source: Sophya/trunk/SigPredictor/lightcmbprim.cc@ 3930

                        // Dominique YVON, CEA/DAPNIA/SPP 02/2000

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "lightcmbprim.h"

#include "fitsspherehealpix.h"

#ifdef __MWERKS__
   #include "macenvvariables.h"
#endif


LightCMBPrim::LightCMBPrim(int_4 nside)
        :nlat(nside) {

  resolution=2*3.14159/(nlat*4.);                       // Radian
  sprintf(Name, "CMB Primordial");
  cout<<"Objet LightCMBPrim :"<<endl;
  cout<<"Resolution de travail nlat= "<<nlat<<endl;
  
  // On reserve la mŽmoire pour la carte gorski
  try {
        pMap =new SphereHEALPix<r_4> (nlat);
  }
  catch (bad_alloc) {
        cerr<<" Memory booking error in LightCMBPrim"<<endl;
        exit(-1);
  }
  cout<<"Vous avez rŽservŽ: "<<12*nlat*nlat*sizeof(float)<<" Octets"<<endl;

  // On lit le fichier gŽnŽrŽ par synfast
                        
  char filename[150];  

#ifndef __MWERKS__
  char* PATHDataLScr=getenv("PATHDataLScr"); 
  sprintf(filename,"%s/cmb_res%04i.fits",PATHDataLScr,nlat);
#else
  sprintf(filename,"%scmb_res%04i.fits",PATHDataLScr,nlat); 
#endif


  FITS_SphereHEALPix<r_4> FSphere(*pMap);
  FSphere.Read(filename);                               // Units:  Microkelvin
        
/*      
  FitsIoServer FitsServer;
  FitsServer.load(*pMap,filename);
*/
}

LightCMBPrim::~LightCMBPrim() {
        delete pMap;
}