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

Timestamp:

Jul 10, 2007, 7:23:38 PM (18 years ago)

Author:

cmv

Message:

ecriture fichiers de tranches du cube cmv 10/07/2007

File:

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

trunk/Cosmo/SimLSS/cmvobserv3d.cc

@@ -39 +39 @@
      <<" -W : write cube in FITS format (complex cube is coded as real cube)"<<endl
      <<" -P : write cube in PPF format"<<endl
+     <<" -S : write cube slices in PPF format"<<endl
      <<" -V : compute variance from real space (for check, default: no)"<<endl
+     <<" -T nth : nombre de threads (si compil multi-thread, default: 0)"<<endl
      <<endl;
 }
@@ -91 +93 @@
  bool computefourier0=false;
  bool use_growth_factor = false;
- unsigned short nthread=4;
+ unsigned short nthread=0;
 
  // *** What to do
@@ -97 +99 @@
  bool wfits = false;
  bool wppf = false;
+ bool wslice = false;
  bool compvarreal = false;
 
@@ -106 +109 @@
 
  char c;
- while((c = getopt(narg,arg,"ha0PWV2Gx:y:z:s:Z:M:A:")) != -1) {
+ while((c = getopt(narg,arg,"ha0PWSV2Gx:y:z:s:Z:M:A:T:")) != -1) {
+  int nth = 0;
   switch (c) {
   case 'a' :
@@ -150 +154 @@
   case 'P' :
     wppf = true;
+    break;
+  case 'S' :
+    wslice = true;
+    break;
+  case 'T' :
+    sscanf(optarg,"%d",&nth);
+    nthread = (nth<1)? 0: nth;
     break;
   case 'h' :
@@ -410 +421 @@
    PrtTim(">>>> End WritePPF");
  }
+ if(wslice) {
+   fluct3d.WriteSlicePPF("cmvobserv3d_s_r0.ppf");
+   PrtTim(">>>> End WriteSlicePPF");
+ }
 
  int_8 nm;
@@ -448 +463 @@
    nm = fluct3d.MeanSigma2(rm,rs2,0.,1e200,true,0.);
  PrtTim(">>>> End Apply poisson on galaxy number");
+ if(wslice) {
+   fluct3d.WriteSlicePPF("cmvobserv3d_s_rn.ppf");
+   PrtTim(">>>> End WriteSlicePPF");
+ }
 
  cout<<"\n--- Convert Galaxy number to HI mass"<<endl;
@@ -479 +498 @@
    PrtTim(">>>> End WritePPF");
  }
+ if(wslice) {
+   fluct3d.WriteSlicePPF("cmvobserv3d_s_r.ppf");
+   PrtTim(">>>> End WriteSlicePPF");
+ }
 
  if(do_agn) {
@@ -502 +525 @@
    fluct3d.WritePPF("cmvobserv3d_rf.ppf",true);
    PrtTim(">>>> End WritePPF");
+ }
+ if(wslice) {
+   fluct3d.WriteSlicePPF("cmvobserv3d_s_rf.ppf");
+   PrtTim(">>>> End WriteSlicePPF");
  }
 
@@ -591 +618 @@
 
 ######################################################
+openppf cmvobserv3d_s_r0.ppf
+openppf cmvobserv3d_s_rn.ppf
+openppf cmvobserv3d_s_r.ppf
+openppf cmvobserv3d_s_rf.ppf
+
+######################################################
 openppf cmvobserv3d.ppf