Changeset 3252 in Sophya for trunk/Cosmo/SimLSS

Timestamp:

May 22, 2007, 11:00:18 AM (18 years ago)

Author:

cmv

Message:

details cmv 22/05/2007

Location:

trunk/Cosmo/SimLSS

Files:

• cmvdefsurv.cc (modified) (3 diffs)
• cmvobserv3d.cc (modified) (1 diff)

trunk/Cosmo/SimLSS/cmvdefsurv.cc

@@ -33 +33 @@
       <<" -P : on donne -x -y -z en Mpc au lieu d\'angles et de redshift"<<endl
       <<" -L lobewidth : taille du lobe d\'observation (FWHM) en arcmin (def= 1\')"<<endl
+      <<"                Si <=0 alors angle solide du lobe = celui du pixel"<<endl
       <<" -O surf,tobs : surface effective (m^2) et temps d\'observation (s)"<<endl
       <<" -N Tsys : temperature du system (K)"<<endl
@@ -275 +276 @@
  //     = NuHi*dz/(1.+z0)^2 * 1/[1-(dz/(2*(1+z0)))^2]
  double dnuhiz = Fr_HyperFin_Par *dred/(unplusz*unplusz)
-               / (1.- (dred/.2/unplusz)*(dred/.2/unplusz));
+               / (1.- pow(dred/.2/unplusz,2.));
  cout<<"\n--- Observation:"<<endl
      <<"      surf_eff="<<surfeff<<" m^2, tobs="<<tobs<<" s"<<endl
@@ -322 +323 @@
  // Elargissement doppler de la raie a 21cm: dNu = vrot/C * Nu(21cm) / (1+z)
  double doplarge = vrot / SpeedOfLight_Cst * nuhiz;
- cout<<"     Doppler width="<<doplarge*1.e3<<" MHz  for rotation width of "<<vrot<<" km/s"<<endl;
- if(doplarge>dnuhiz) {
+ double dzvrot = vrot / SpeedOfLight_Cst * unplusz;
+ cout<<"     Doppler width="<<doplarge*1.e3<<" MHz  for rotation width of "<<vrot<<" km/s"<<endl
+     <<"                 dx= "<<dzvrot<<" a z="<<redshift<<endl;
+ if(doplarge>dnuhiz)
    cout<<"Warning: doppler width "<<doplarge<<" GHz > "<<dnuhiz<<" GHz redshift bin width"<<endl;
- }
 
  // Synchrotron

trunk/Cosmo/SimLSS/cmvobserv3d.cc

@@ -396 +396 @@
    cout<<"...D(z=0)="<<growth(0.)<<"  D(z="<<zref<<")="<<growth(zref)<<endl;
    fluct3d.ApplyGrowthFactor();
-   rmin,rmax; rgen.MinMax(rmin,rmax);
+   rgen.MinMax(rmin,rmax);
    cout<<"rgen.Min = "<<rmin<<" , Max="<<rmax<<endl;
    PrtTim(">>>> End Applying growth factor");