Changeset 4028 in Sophya for trunk/Cosmo/RadioBeam/pknoise.cc

Timestamp:

Oct 18, 2011, 10:56:28 AM (14 years ago)

Author:

ansari

Message:

Debug de calcul Pnoise3D(k) avec D_A(z), Reza 18/10/2011

File:

• trunk/Cosmo/RadioBeam/pknoise.cc (modified) (4 diffs)

trunk/Cosmo/RadioBeam/pknoise.cc

@@ -90 +90 @@
   double H_z=120.5;  // Hubble_param(z) 
   int box3dsz[3]={512,512,256};
-  double cellsz[3]={3.,3.,1.5};
+  double cellsz[3]={3.,3.,3.};
   int NBINPK=256;
   int prtlev=0;
@@ -159 +159 @@
     redshift2da.ReadXYFromFile(interpfilename);
     SLinInterp1D redshift2hz;
-    interpfilename="redshift_hz.txt";
+    interpfilename="redshift_Hz.txt";
     redshift2hz.ReadXYFromFile(interpfilename);
     comov_dA_z=redshift2da(z_Redshift);
@@ -211 +211 @@
     Vector angscales(box3dsz[2]);
     angscales=angscale;
-    if (NMAX<0)  { box3dsz[0]=256; box3dsz[1]=256;  box3dsz[3]=128; }
+    if (NMAX<0)  { box3dsz[0]=256; box3dsz[1]=256;  box3dsz[2]=128; }
     else {
       angscale=comov_dA_z;
@@ -225 +225 @@
          << box3dsz[2] << ")" << endl;
     Four3DPk m3d(rg,box3dsz[0]/2,box3dsz[1],box3dsz[2]);
-    cout << " pknoise[2.d]: m3d.SetCellSize(" << dkxmpc << "," << dkympc << "," << dkzmpc << endl;
+    cout << " pknoise[2.d]: m3d.SetCellSize(" << dkxmpc << "," << dkympc << "," << dkzmpc 
+         << ") cell size (Mpc) : " << cellsz[0] << "x" << cellsz[1] << "x" << cellsz[2] << endl;
     m3d.SetCellSize(dkxmpc, dkympc, dkzmpc);
     m3d.SetPrtLevel(prtlev,prtmod);