Changeset 3973 in Sophya for trunk/Cosmo/RadioBeam/cubedef.h


Ignore:
Timestamp:
Apr 18, 2011, 5:30:44 PM (14 years ago)
Author:
ansari
Message:

Corrections diverses: choix lobe gaussien/triangle et specif DishDiameter au lieu de DoL ds applobe/calcpk2, possibilite application lobe freq.independante ds applobe, Reza 18/04/2011

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/Cosmo/RadioBeam/cubedef.h

    r3829 r3973  
    77// Definition tailles et bornes du cube angles (X,Y), frequence (Z)
    88
    9 static sa_size_t NTheta=360;
    10 static sa_size_t NPhi=360;
     9static sa_size_t NTheta=600;
     10static sa_size_t NPhi=800;
    1111static sa_size_t NFreq = 256;
    1212
     13// Carte 30x50 deg, centre en delta = +10 deg , alpha= 150 deg = 10h00
    1314static double Theta0Degre = 65.;  // -> Delta = +30 deg
    14 static double Phi0Degre = 135.;   // -> alpha = 9h00
     15static double Phi0Degre = 130.;   // -> alpha = 8.66 heures
    1516static double ThetaSizeDegre = 30.;  // Taille de la carte en degre selon alpha (axe X)
    16 static double PhiSizeDegre = 30.;    // Taille de la carte en degre selon delta (axe Y)
     17static double PhiSizeDegre = 40.;    // Taille de la carte en degre selon delta (axe Y)
    1718
    1819/* ---  Parametres avec generation Reza from HASLAM + NVSS
    … …  
    3435
    3536// z = 0.6, freq_center=884 MHz , dA_comov ~= 2185 Mpc , 
    36 //       5 arcmin -> 3 Mpc , 0.550 MHz -> 3 Mpc  (H(ze) ~ 96 km/s/Mpc)
     37//       5 arcmin -> 3 Mpc (3.17 Mpc) , 0.550 MHz -> 3 Mpc  (H(ze) ~ 96 km/s/Mpc)
     38//       3 arcmin -> 1.9 Mpc, 0.5 MHz -> 2.8 Mpc
    3739// Taille de cellule pour le calcul du spectre de puissance 3D
    3840static double ComovDA = 2185.;
    39 static double XCellSizeMpc = 3.;
    40 static double YCellSizeMpc = 3.;
    41 static double ZCellSizeMpc = 3.;
    42 
     41static double XCellSizeMpc = 1.9;
     42static double YCellSizeMpc = 1.9;
     43static double ZCellSizeMpc = 2.8;
     44// Nb de bin de HProf pour calcul P(k)
     45static int HPk_NBin = 384;
    4346
    4447
    … …  
    5255static double sigPLidx2 = 0.15;
    5356
    54 static double PLidxSrc = -2.2;  // index de la loi de puissance des sources 
     57// Generation de la loi de puissance des sources :
     58static double PLidxSrc = -2.0;  // index de la loi de puissance des sources 
    5559static double sigPLidxSrc = 0.15;  // Sigma de la variation (gaussienne) de l'index 
    5660
Note: See TracChangeset for help on using the changeset viewer.