diff --git a/source/SourcePopulation.cc b/source/SourcePopulation.cc
index 3d85ca8448b561b234e5863c72bd21a3ae55c123..82c572dc9ca29c9e6f878a8293a527d7d90b9916 100644
--- a/source/SourcePopulation.cc
+++ b/source/SourcePopulation.cc
@@ -45,6 +45,8 @@ int  SourcePopulation::init()
   IDL_plot_control =  1; //AWS20200623
   IDL_psym         = -6; //AWS20200623
 
+  print_ctl        =  0; //AWS20200720
+
   strcpy(galdef_ID, "sourcepop_needs_galdef_ID" ); //AWS20200701
   strcpy(psfile_tag,"sourcepop_needs_psfile_tag"); //AWS20200701
 
@@ -790,12 +792,12 @@ int SourcePopulation::analyse_sample()
           pointing pointing_=pointing(theta,phi);
 
 
-          if(verbose!=4)cout<<"source statistics sample i="<<i<<" Fermi sensitivity  l="<<sample_l[i] <<" b="<<sample_b[i]
-				<<" theta="<<theta<<" phi="<<phi;  //AWS20200716
+          if(print_ctl==1)cout<<"source statistics sample i="<<i<<" Fermi sensitivity  l="<<sample_l[i] <<" b="<<sample_b[i]
+				<<" theta="<<theta<<" phi="<<phi;  //AWS20200720
 	  //          cout<<" interpolated value="<<healpix_skymap_Fermi_sensitivity.interpolated_value(pointing_)<<endl;
           flux_detection_limit_sensitivity_file = healpix_skymap_Fermi_sensitivity.interpolated_value(pointing_); //AWS20200713 was flux_detection_limit
 
-          if(verbose!=4)cout<<" flux detection limit_sensitivity_file="<<flux_detection_limit_sensitivity_file<<" sample flux="<<sample_flux[i]<<endl;//AWS20200716
+          if(print_ctl==1)cout<<" flux detection limit_sensitivity_file="<<flux_detection_limit_sensitivity_file<<" sample flux="<<sample_flux[i]<<endl;//AWS20200720
          }
 
         dlnN_dlnS  [ii]++;
@@ -1065,12 +1067,12 @@ int SourcePopulation::analyse_sample()
           double phi  =    sample_l[i] *dtr;
           pointing pointing_=pointing(theta,phi);
 
-          if(verbose!=4)cout<<"skymap generation source sample i="<<i<<" Fermi sensitivity  l="<<sample_l[i] <<" b="<<sample_b[i]
-				<<" theta="<<theta<<" phi="<<phi; //AWS20200716
+          if(print_ctl==1)cout<<"skymap generation source sample i="<<i<<" Fermi sensitivity  l="<<sample_l[i] <<" b="<<sample_b[i]
+				<<" theta="<<theta<<" phi="<<phi; //AWS20200720
 	  //          cout<<" interpolated value="<<healpix_skymap_Fermi_sensitivity.interpolated_value(pointing_)<<endl;
           flux_detection_limit_sensitivity_file = healpix_skymap_Fermi_sensitivity.interpolated_value(pointing_); //AWS20200713
 
-          if(verbose!=4)cout<<" flux detection limit_sensitivity_file="<<flux_detection_limit_sensitivity_file<<" sample flux="<<sample_flux[i]<<endl;//AWS20200716
+          if(print_ctl==1)cout<<" flux detection limit_sensitivity_file="<<flux_detection_limit_sensitivity_file<<" sample flux="<<sample_flux[i]<<endl;//AWS20200720
          }
 
        ii=int( (log10(sample_flux[i])-lnS_min)/dlnS  );
diff --git a/source/SourcePopulation.h b/source/SourcePopulation.h
index 305d4b71b68d0b08a37cef9e539cbe009a4113e2..e851c6f7d2ffcaeb417144ecd73fdcf1276dfe08 100644
--- a/source/SourcePopulation.h
+++ b/source/SourcePopulation.h
@@ -162,6 +162,7 @@ class SourcePopulation
   char* Fermi_sensitivity_file;                          //AWS20170112
 
   int verbose;
+  int print_ctl;                                         //AWS20200720
 
   ofstream *txt_stream;