Changeset 3271 in Sophya for trunk/Cosmo/SimLSS/cmvobserv3d.cc

Timestamp:

Jun 20, 2007, 6:50:48 PM (18 years ago)

Author:

cmv

Message:

intro evolution bruit avec redshift cmv 20/06/2007

File:

• trunk/Cosmo/SimLSS/cmvobserv3d.cc (modified) (5 diffs)

trunk/Cosmo/SimLSS/cmvobserv3d.cc

@@ -30 +30 @@
      <<" -z nz,dz : size along z axis (redshift axis, npix,Mpc)"<<endl
      <<" -Z zref : redshift for the center of the simulation cube"<<endl
-     <<" -s snoise : gaussian noise sigma in equivalent Msol"<<endl
+     <<" -s snoise,evol : gaussian noise sigma in equivalent Msol"<<endl
+     <<"                  if evol>0 noise evolved with distance (def no)"<<endl
      <<" -2 : compute also 2D spectrum (default: no)"<<endl
      <<" -M schmin,schmax,nsch : min,max mass and nb points for schechter HI"<<endl
@@ -80 +81 @@
  // *** Niveau de bruit
  double snoise= 0.;   // en equivalent MSol
+ int isnoise_evol = 0;
 
  // *** AGN
@@ -125 +127 @@
     break;
   case 's' :
-    sscanf(optarg,"%lf",&snoise);
+    sscanf(optarg,"%lf,%d",&snoise,&isnoise_evol);
     break;
   case 'Z' :
@@ -169 +171 @@
      <<"), schmax="<<schmax<<" ("<<lschmax<<") Msol"
      <<", schnpt="<<schnpt<<endl;
- cout<<"snoise="<<snoise<<" equivalent Msol"<<endl;
+ cout<<"snoise="<<snoise<<" equivalent Msol, evolution="<<isnoise_evol<<endl;
  if(do_agn)
    cout<<"AGN: <log10(Jy)>="<<lfjy_agn<<" , sigma="<<lsigma_agn
@@ -487 +489 @@
 
  if(snoise>0.) {
-   cout<<"\n--- Add noise to HI Flux snoise="<<snoise<<endl;
-   fluct3d.AddNoise2Real(snoise);
+   cout<<"\n--- Add noise to HI Flux snoise="<<snoise<<", evolution="<<isnoise_evol<<endl;
+   fluct3d.AddNoise2Real(snoise,(isnoise_evol>0? true:false));
      nm = fluct3d.MeanSigma2(rm,rs2);
    PrtTim(">>>> End Add noise");