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 | |
---|