Changeset 3789 in Sophya for trunk/Cosmo/RadioBeam/calcpk.cc

Timestamp:

Jun 27, 2010, 4:57:36 PM (15 years ago)

Author:

ansari

Message:

Corrections diverses, Reza 27/06/2010

File:

• trunk/Cosmo/RadioBeam/calcpk.cc (modified) (3 diffs)

trunk/Cosmo/RadioBeam/calcpk.cc

@@ -20 +20 @@
 #include "qhist.h" 
 #include "lobe.h" 
+#include "cubedef.h" 
 
 #include "histinit.h"
@@ -86 +87 @@
     RandomGenerator rg;
     Four3DPk pkc(four3d, rg);
-    
+    double dkxmpc = DeuxPI/(double)inmap.SizeX()/XCellSizeMpc;
+    double dkympc = DeuxPI/(double)inmap.SizeY()/YCellSizeMpc;
+    double dkzmpc = DeuxPI/(double)inmap.SizeZ()/ZCellSizeMpc;
+    pkc.SetCellSize(dkxmpc, dkympc, dkzmpc);
+
     HProf hp = pkc.ComputePk(0.,256);
     {
@@ -92 +97 @@
       POutPersist po(outname);
       po << hp;
+      outname = "f3d_" + outname;
+      cout << "calcpk[4.b] : writing four3d to PPF file  " << outname << endl;
+      POutPersist pof(outname);
+      pof << four3d;
     }
     tm.Split(" Done CalcP(k) ");