Changeset 3756 in Sophya for trunk/Cosmo/RadioBeam/anapkn.pic

Timestamp:

Apr 28, 2010, 3:41:15 PM (15 years ago)

Author:

ansari

Message:

Ajout des programmes de calcul de la sensibilite de l'interfero (plan(u,v), PNoise(k)) , Reza 28/04/2010

File:

• trunk/Cosmo/RadioBeam/anapkn.pic (modified) (5 diffs)

trunk/Cosmo/RadioBeam/anapkn.pic

@@ -54 +54 @@
 
 ##  ---- VOL1,VOL2,VOL3 
-##  VOL1: 10deg x 10deg x Delta_z=0.2 : 335 x 335 x 680 Mpc^3 = 76 10^6 Mpc^3
-##  VOL2: 20deg x 20deg x Delta_z=0.2 : VOL2 = 300 10^6 Mpc^3
-##  VOL3: 60deg x 60deg x Delta_z=0.2 : VOL4 = 2700 10^6 Mpc^3
+##  VOL0: 10deg x 10deg x Delta_z=0.2 : 335 x 335 x 680 Mpc^3 = 76 10^6 Mpc^3
+##  VOL1: 20deg x 20deg x Delta_z=0.2 : VOL1 = 300 10^6 Mpc^3
+##  VOL2: 360deg x 10deg x Delta_z=0.2 : VOL2 = 2700 10^6 Mpc^3
+##  VOL3: 360deg x 30deg x Delta_z=0.2 : VOL3 = 10800 10^6 Mpc^3
 ##  VOL4: 40deg x 40deg x Delta_z=0.2 : VOL3 = 1200 10^6 Mpc^3
-set angcov ( '10x10 deg^2' '20x20 deg^2' '60x60 deg^2' '40x40 deg^2' )
-VOL1 = 76e6
-VOL2 = 300e6
-VOL3 = 2700e6
+set angcov ( '20x20 deg^2' '360x10 deg^2' '360x30 deg^2' '40x40 deg^2' )
+VOL0 = 76e6
+VOL1 = 300e6
+VOL2 = 2700e6
+VOL3 = 10800e6
 VOL4 = 1200e6
 
@@ -104 +106 @@
 endscript
 
+#  Setup 64 x D=5m sur 100mx100m , 10000 m^2
+defscript setup5
+  vpixz05 = 100
+  vpix = 100
+  kappa = 100*100/1200
+  kappa = $kappa*$kappa
+  echo "setup-5: VPix=$vpix  Kappa=$kappa"
+endscript
 
 
@@ -112 +122 @@
 
 #  --- Temps d'integration  a z=0.5 avec FOV=4x4 degres^2  (D=5m) 
-#  120 jours pour VOl1 en ~ 7 passes (7 couvertures elementaires)
-#  120 jours pour VOl2 en ~ 4x7 passes (28 couvertures elementaires)
-#  120 jours pour VOl3 en ~ 16x7 passes (225 couvertures elementaires)
-set ndays ( '120 days' '120 days'  '120 days' '120 days' ) 
-TINT1 = 1.5e6
-TINT2 = 3.7e5
-TINT3 = 4.5e4
+#  120 jours pour VOl0 en ~ 7 passes (7 couvertures elementaires)
+#  360 jours pour VOl2 en ~ 16x7 passes (225 couvertures elementaires)
+#  360 jours pour VOl1 en ~ 4x7 passes (28 couvertures elementaires)
+#  360 jours pour VOl2 en ~ 90x3 passes (270 couvertures elementaires)
+#  360 jours pour VOl3 en ~ 90x9 passes (810 couvertures elementaires)
+
+set ndays ( '360 days' '360 days'  '360 days' '360 days' ) 
+TINT1 = 1.1e6
+TINT2 = 1.1e5
+TINT3 = 3.8e4
 TINT4 = 1.e5
 
 #  Bruit Tsys = 100 K = 10^5 mK 
 TSYS = 1e5
+#  Bruit Tsys = 50 K = 5.10^4 mK 
+TSYS = 5e4
 
 DELNU = 1.e6
@@ -151 +166 @@
 set vols ( $VOL1 $VOL2 $VOL3 )
 set tints ( $TINT1 $TINT2 $TINT3 )
-set cols ( green blue red )
+set cols ( red orange blue )
 
 #############################################################
@@ -239 +254 @@
   xyl = "xylimits=0.01,0.5,1.e-3,0.4 logx logy minorticks"
   set vols ( $VOL1 $VOL2 $VOL3 )
-  set cols ( green blue red )
+  set cols ( red orange blue )
   csigpk $vols[0]
   n/plot hpkz.$CSPK/$kk%$kk ! ! "notit nsta connectpoints $cols[0] $xyl line=solid,2"