Changeset 3565 in Sophya for trunk/Cosmo/SimLSS

Timestamp:

Jan 30, 2009, 10:18:23 AM (17 years ago)

Author:

ansari

Message:

Ajout options -d/-D pour ecriture d_rho/rho en PPF/fits et sortie du programme apres - Reza 30/01/2009

File:

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

trunk/Cosmo/SimLSS/cmvobserv3d.cc

@@ -56 +56 @@
      <<" -V : compute variance from real space (for check, default: no)"<<endl
      <<" -T nth : nombre de threads (si compil multi-thread, default: 0)"<<endl
+     <<" -D fileName : write only real cube d_rho/rho in FITS format and exit "<<endl
+     <<" -d fileName : write only real cube d_rho/rho in PPF format and exit "<<endl
+
      <<endl;
      ////<<" -A <log10(S_agn in Jy at 1.4 GHz)>,sigma,powlaw :"<<endl
@@ -126 +129 @@
  bool wslice = false;
  bool compvarreal = false;
+ //Ajoute par Reza : Ecriture d_rho/rho en fits ou PPF et 
+ bool fgrzproc = false;
+ string rzppfname = "";
+ string rzfitsname = "";
+ 
  unsigned long ntnent = 10000;  // 0 = do not fill NTuple
 
@@ -136 +144 @@
 
  char c;
- while((c = getopt(narg,arg,"ha0PWSV2U:G:F:x:y:z:s:Z:M:A:T:N:Q:R:8:")) != -1) {
+ while((c = getopt(narg,arg,"ha0PWSB2U:G:F:x:y:z:s:Z:M:A:T:N:Q:R:8:D:d")) != -1) {
   int nth = 0;
   switch (c) {
@@ -204 +212 @@
   case 'S' :
     wslice = true;
+    break;
+  case 'd' :
+    fgrzproc = true;
+    rzppfname = optarg;
+    break;
+  case 'D' :
+    fgrzproc = true;
+    rzfitsname = optarg;
     break;
   case 'T' :
@@ -551 +567 @@
    PrtTim(">>>> End Check variance sigmaR in real space");
  }
-
+ if(fgrzproc) {
+   cout << "--- Flag -d/-D specified, saving d_rho/rho and exiting ..." << endl;
+   if (rzfitsname.length()>0) {
+         cout<<" Writing d_rho/rho to fits file " << rzfitsname << endl;
+     fluct3d.WriteFits(rzfitsname);
+   }
+   if (rzppfname.length()>0) {
+         cout<<" Writing d_rho/rho to PPF file " << rzppfname << endl;
+     fluct3d.WritePPF(rzppfname);
+   }
+   PrtTim(">>>> End WriteFits/PPF_d_rho/rho");
+   cout<<"-------------------- END ------------------"<<endl;
+   return 0;
+ }
  //-----------------------------------------------------------------
  if(no_poisson_type!=0) {