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