diff --git a/source/SourcePopulation.cc b/source/SourcePopulation.cc
index e3eba4a7b41fd22b2b020968d2a29241334d1e50..650342ffa2b184be8f097def04529fe4acae43bd 100644
--- a/source/SourcePopulation.cc
+++ b/source/SourcePopulation.cc
@@ -1579,6 +1579,25 @@ int SourcePopulation::analyse_sample()
    *txt_stream<<endl;
 
 
+   *txt_stream<<"NS_soplimit = ["; 
+   for (ii=0;ii< n_dlnN_dlnS;ii++)
+       { *txt_stream<<dlnN_dlnS_soplimit[ii];
+           if(ii<n_dlnN_dlnS-1)    *txt_stream<<",";
+       }
+   *txt_stream<<"]";
+   *txt_stream<<";  idl format";
+   *txt_stream<<endl;
+
+
+   *txt_stream<<"NS_sublimit = ["; 
+   for (ii=0;ii< n_dlnN_dlnS;ii++)
+       { *txt_stream<<dlnN_dlnS_sublimit[ii];
+           if(ii<n_dlnN_dlnS-1)    *txt_stream<<",";
+       }
+   *txt_stream<<"]";
+   *txt_stream<<";  idl format";
+   *txt_stream<<endl;
+
 
    //   *txt_stream<<"device,file='"<<title<<".ps'          ;idl"<<endl;
 
@@ -1601,6 +1620,16 @@ int SourcePopulation::analyse_sample()
      *txt_stream<<"oplot,S,NS,psym="<<IDL_psym<<"   ; idl"<<endl;
 
 
+   if(IDL_plot_control==3)                     // overplot on start plot
+     *txt_stream<<"oplot,S,NS,psym="<<IDL_psym<<"   ; idl"<<endl;
+
+   if(IDL_plot_control==3)                      // overplot on start plot
+     *txt_stream<<"oplot,S,NS_soplimit,psym="<<IDL_psym<<"   ; idl"<<endl;
+
+   if(IDL_plot_control==3                     ) // overplot on start plot
+     *txt_stream<<"oplot,S,NS_sublimit,psym="<<IDL_psym<<"   ; idl"<<endl;
+
+
    // python format
 
   *txt_stream<<"title='" <<title<< "' # python format"<<endl;
diff --git a/source/gen_source_population.cc b/source/gen_source_population.cc
index 35f6b153ded06be24446bab86e4081dde6c7ead4..b14bf22392a9a32c707215d7219cdab0b17828c3 100644
--- a/source/gen_source_population.cc
+++ b/source/gen_source_population.cc
@@ -375,8 +375,8 @@ int Galplot::gen_source_population()
 
    sourcepop5.print();
 
-   sourcepop5.IDL_plot_control=2;                                                                      //AWS20200623 oplot on first plot
-   sourcepop5.IDL_psym        =0;                                                                      //AWS2020623  just line
+   sourcepop5.IDL_plot_control=3;                                          //AWS20200714 oplot on first plot model NS, soplimit, sublimit
+   sourcepop5.IDL_psym        =0;                                          //AWS20200714 just line
 
    sourcepop5.analyse_sample();