Changeset 3283 in Sophya for trunk/Cosmo

Timestamp:

Jul 25, 2007, 2:41:39 PM (18 years ago)

Author:

cmv

Message:

scale du spectre final cmv 24/07/2007

Location:

trunk/Cosmo/SimLSS

Files:

• cmvobserv3d.cc (modified) (3 diffs)
• genefluct3d.cc (modified) (2 diffs)
• genefluct3d.h (modified) (1 diff)

trunk/Cosmo/SimLSS/cmvobserv3d.cc

@@ -208 +208 @@
  GrowthFactor growth(om0,ol0);
  // GrowthFactor growth(1.,0.); // D(z) = 1/(1+z)
+ double growth_at_z = growth(zref);
+ cout<<"...Growth factor at z="<<zref<<" = "<<growth_at_z<<endl;
 
  PkSpectrum0 pk0(pkini,tf);
@@ -439 +441 @@
    double varr;
    fluct3d.VarianceFrReal(R,varr);
+   cout<<"...Computed variance = "<<varr
+       <<" , Theorical variance at (z=0) = "<<pow(sigmaR,1.)
+       <<" , at (z="<<zref<<") = "<<pow(sigmaR*growth_at_z,1.)<<endl;
    PrtTim(">>>> End Check variance sigmaR in real space");
  }
@@ -511 +516 @@
  }
 
+ double scalecube = -1., offsetcube=0.;
  if(snoise>0.) {
    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);
+     if(rs2>0. && sr2int>0.) scalecube = sqrt(sr2int)/sqrt(rs2);
+     offsetcube = -rm;
    PrtTim(">>>> End Add noise");
  }
 
+ if(scalecube>0.) {
+   cout<<"\n--- Scale and offset scale="<<scalecube
+       <<" off="<<offsetcube<<endl;
+   fluct3d.ScaleOffset(scalecube,offsetcube);
+ }
+ 
  if(wfits) {
    fluct3d.WriteFits("!cmvobserv3d_rf.fits");

trunk/Cosmo/SimLSS/genefluct3d.cc

@@ -450 +450 @@
 // On ecrit 3 tranches du cube selon chaque axe
 {
- cout<<"--- GeneFluct3D::WriteSlicePPF: Writing Cube Slices"<<cfname<<endl;
+ cout<<"--- GeneFluct3D::WriteSlicePPF: Writing Cube Slices "<<cfname<<endl;
  try {
 
@@ -996 +996 @@
                <<" because out of range=]"<<vmin<<","<<vmax<<"["<<endl;
  return nbad;
+}
+
+void GeneFluct3D::ScaleOffset(double scalecube,double offsetcube)
+// Replace  "V"  by  "scalecube * ( V + offsetcube )"
+{
+ if(lp_>0) cout<<"--- ScaleCube "<<scalecube<<endl;
+
+ for(long i=0;i<Nx_;i++) for(long j=0;j<Ny_;j++) for(long l=0;l<Nz_;l++) {
+   int_8 ip = IndexR(i,j,l);
+   data_[ip] = scalecube * ( data_[ip] + offsetcube );
+ }
+
+ return;
 }
 

trunk/Cosmo/SimLSS/genefluct3d.h

@@ -86 +86 @@
   int_8 NumberOfBad(double vmin=-1.e+150,double vmax=1.e+150);
   int_8 SetToVal(double vmin, double vmax,double val0=0.);
+  void  ScaleOffset(double scalecube=1.,double offsetcube=0.);
 
   void TurnFluct2Mass(void);