[676] | 1 | ############################################################### |
---|
| 2 | ### Estimations pour PAON-1 (Interfero de test a Nancay ) #### |
---|
| 3 | ####################### Avril 2012 ########################## |
---|
| 4 | ############################################################### |
---|
| 5 | |
---|
| 6 | setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks' |
---|
| 7 | |
---|
| 8 | ##### 1ere partie ############## |
---|
| 9 | # Trace des radio sources a partir du catalog NVSS |
---|
| 10 | readfits nvss.fits |
---|
| 11 | plot2dw nvss C_RAJ2000/15. C_DEJ2000 S1_4 S1_4>5000 'nsta orange marker=circle,7' |
---|
| 12 | plot2dw nvss C_RAJ2000/15. C_DEJ2000 S1_4 S1_4>10000 'nsta orange marker=fcircle,9 same' |
---|
| 13 | plot2dw nvss C_RAJ2000/15. C_DEJ2000 S1_4 S1_4>20000 'nsta red marker=fcircle,11 same' |
---|
| 14 | plot2dw nvss C_RAJ2000/15. C_DEJ2000 S1_4 S1_4>50000 'nsta blueviolet marker=fcircle,13 same' |
---|
| 15 | |
---|
| 16 | settitle 'NVSS radio sources (flux>5 Jy)' ' ' 'font=helvetica,bold,16' |
---|
| 17 | setaxelabels 'RA (alpha) - hours' 'DEC (delta)' 'font=helvetica,bold,16' |
---|
| 18 | |
---|
| 19 | set lines ( 'flux > 10 Jy' 'flux > 20 Jy' 'flux > 50 Jy' ) |
---|
| 20 | set cols ( orange red blueviolet ) |
---|
| 21 | textdrawer lines cols 'font=helvetica,bold,16 frame' |
---|
| 22 | addtext 20. 62. 'CasA ~500 Jy' 'font=helvetica,bold,14 blueviolet' |
---|
| 23 | addtext 17. 36. 'CygA ~800 Jy' 'font=helvetica,bold,14 blueviolet' |
---|
| 24 | |
---|
| 25 | # Latitude 47 degres |
---|
| 26 | addline 0 44 24 44 'line=solid,1' |
---|
| 27 | addline 0 50 24 50 'line=solid,1' |
---|
| 28 | |
---|
| 29 | addline 0 67 24 67 'line=dotted,1' |
---|
| 30 | addline 0 27 24 27 'line=dotted,1' |
---|
| 31 | |
---|
| 32 | w2ps |
---|
| 33 | w2eps nvss_bright.eps |
---|
| 34 | |
---|
| 35 | ############################################# |
---|
| 36 | ########## 2eme partie ############## |
---|
| 37 | ### Calcul bruit et franges d'interference |
---|
| 38 | ############################################# |
---|
| 39 | # Calcul bruit |
---|
| 40 | WDISH = 3. |
---|
| 41 | lambda = 0.21 |
---|
| 42 | # Vitesse angulaire Omega en degre/second = 15./3600.*cos(delta) |
---|
| 43 | OMEGA = cos(Pi*47./180)*15./3600. |
---|
| 44 | set ANG $OMEGA*x |
---|
| 45 | set T (2*M_PI*$WDISH/$lambda*sin(M_PI*$ANG/180.)) |
---|
| 46 | set XX ((fabs($T)<1e-8)?1.:(sin($T)/$T)) |
---|
| 47 | func $XX*$XX -1800 1800 240 |
---|
| 48 | |
---|
| 49 | # Niveau de bruit en Jy |
---|
| 50 | TSYS = 100 |
---|
| 51 | k_B = 1.38e-23 |
---|
| 52 | A = 0.75*Pi*$WDISH*$WDISH/4. |
---|
| 53 | DNU = 5.e6 |
---|
| 54 | TINT = 15 |
---|
| 55 | NOISE = 2*$k_B*$TSYS/$A/sqrt($DNU*$TINT)*1.e26 |
---|
| 56 | echo 'NoiseLevel: ' $NOISE 'Jansky ' |
---|
| 57 | |
---|
| 58 | |
---|
| 59 | ##### Calcul N-tuple avec Visib(temps) |
---|
| 60 | |
---|
| 61 | # Executer pour FLXSRC=500, NOISE=10 |
---|
| 62 | del visi |
---|
| 63 | c++execfrf paon1.cc |
---|
| 64 | mv visi visi500 |
---|
| 65 | # Executer pour FLXSRC=20, NOISE=5 (4-5 d'integrations/transits) |
---|
| 66 | del visi |
---|
| 67 | c++execfrf paon1.cc |
---|
| 68 | |
---|
| 69 | mv visi visi20 |
---|
| 70 | |
---|
| 71 | #### Traces |
---|
| 72 | set titres ( '500 Jy - single transit (Noise=10 Jy)' '20 Jy - 4 transits (Noise=3 Jy)' ) |
---|
| 73 | set O ( visi500 visi20 ) |
---|
| 74 | for i 0:2 |
---|
| 75 | newwin 1 2 |
---|
| 76 | plot2de $O[i] time vI 0. err*1.4 1 'marker=fcircle,5 cpts notit nsta' |
---|
| 77 | settitle "$titres[i]" ' ' 'font=helvetica,bold,16' |
---|
| 78 | setaxelabels 'Time (s)' 'Vis=I - Jy ' ' font=helvetica,bolditalic,14' |
---|
| 79 | plot2de $O[i] time vre 0. err 1 'marker=fcircle,5 red cpts notit nsta' |
---|
| 80 | setaxelabels 'Time (s)' 'Real[Vis] - Jy ' ' font=helvetica,bolditalic,14' |
---|
| 81 | w2ps |
---|
| 82 | end |
---|
| 83 | |
---|
| 84 | |
---|
| 85 | |
---|