#include "sopnamsp.h" #include "machdefs.h" #include #include #include #include #include #include #include "fioarr.h" #include "skymap.h" // set rname = gsm820948 // set fname = ( `cat ${rname}.list | awk '{print $3}' | sed 's/data/ppf/'` ) // cmvgsm2cube -o ${rname}_cube.ppf -p 360,0.1666667,160. -t 360,0.1666667,30. $fname void usage(void) { cout <<"cmvgsm2cube [options] liste_des_spheres_healpix"< phi0 + nphi*dphi (degres)"< theta0 + ntheta*dtheta (degres)"< sphname; for(int i=optind;i Theta(ntheta); int nthgood = 0; float thmin = 1.e30, thmax = -1.e30; for(int it=0;itM_PI) continue; // [0,Pi] Theta(it) = t; nthgood += 1; if(tthmax) thmax = t; } cout<<"Theta: ngood="< sph; pis>>PPFNameTag("sph")>>sph; Freq(is) = (r_4)sph.Info()["FMHz"]; for(int it=0;it