| 1 | ########################################################################
|
|---|
| 2 | ### Script de trace de PNoise(k) et reponse dans le plan (u,v)
|
|---|
| 3 | ### de l'interferometre a partir du fichier PPF produit par pknoise.cc
|
|---|
| 4 | ### Janvier 2011 , BAORadio/Reza
|
|---|
| 5 | ########################################################################
|
|---|
| 6 |
|
|---|
| 7 | echo ' -----> plpknew.pic : opening PPF files ... "
|
|---|
| 8 | delobjs *
|
|---|
| 9 | # openppf ../../PkNoise/cmvhpkz.ppf
|
|---|
| 10 | # openppf pkz0p25.ppf
|
|---|
| 11 | openppf pkz0p7.ppf
|
|---|
| 12 | rename hpkz hpkz0p7
|
|---|
| 13 | openppf pkz1p0.ppf
|
|---|
| 14 | rename hpkz hpkz1p0
|
|---|
| 15 |
|
|---|
| 16 | defscript pknopen
|
|---|
| 17 | set f $1
|
|---|
| 18 | set nm $2
|
|---|
| 19 | echo '------ pknopen File= ' $f ' name=' $nm
|
|---|
| 20 | openppf $f
|
|---|
| 21 | rename h1dnoise h1dn$nm
|
|---|
| 22 | rename h1drep h1dr$nm
|
|---|
| 23 | rename fracmodok fracok$nm
|
|---|
| 24 | mv dtnoise pkn$nm
|
|---|
| 25 | rename h2drep h2dr$nm
|
|---|
| 26 | endscript
|
|---|
| 27 |
|
|---|
| 28 | pknopen pknf11x11.ppf A121d
|
|---|
| 29 | pknopen pknnan128.ppf B128d
|
|---|
| 30 | pknopen pknconfC.ppf C129d
|
|---|
| 31 | pknopen pknD50m.ppf D50m
|
|---|
| 32 | pknopen pknD75m.ppf D75m
|
|---|
| 33 | pknopen pknD100m.ppf D100m
|
|---|
| 34 | pknopen pknD200m.ppf D200m
|
|---|
| 35 | pknopen pknf20x20.ppf E400d
|
|---|
| 36 | pknopen pknf4cyl.ppf F4cyl
|
|---|
| 37 | pknopen pknf4cylp.ppf F4cylp
|
|---|
| 38 | pknopen pknf8cyl.ppf G8cyl
|
|---|
| 39 | pknopen pknf8cylp.ppf G8cylp
|
|---|
| 40 |
|
|---|
| 41 |
|
|---|
| 42 | listobjs
|
|---|
| 43 |
|
|---|
| 44 | # A z = 0.7
|
|---|
| 45 | z = 0.7
|
|---|
| 46 | c = 3.e5
|
|---|
| 47 | H = 102
|
|---|
| 48 | Da = 2488
|
|---|
| 49 | nu21 = 1.42e9
|
|---|
| 50 |
|
|---|
| 51 | echo '----> Executing anapkn.pic'
|
|---|
| 52 | exec anapkn.pic
|
|---|
| 53 | setup5
|
|---|
| 54 | # scalewz 0.7 2488 102 1
|
|---|
| 55 | # scalewz 0.25 989.13 80.26 1
|
|---|
| 56 | scalewz 1 3300 120.5 1
|
|---|
| 57 | Da = 3300
|
|---|
| 58 |
|
|---|
| 59 |
|
|---|
| 60 | y1 = 400*$cct21
|
|---|
| 61 | y2 = 7e4*$cct21
|
|---|
| 62 | xyl = "xylimits=0.005,0.5,$y1,$y2 logx logy minorticks"
|
|---|
| 63 | # H = 71.9 km/s/Mpc = 1.0271 x 70 km/s
|
|---|
| 64 | h70 = 1.0271
|
|---|
| 65 | h70cube = $h70*$h70*$h70
|
|---|
| 66 | set kk pow(10.,x)/$h70
|
|---|
| 67 |
|
|---|
| 68 | defscript plpklss
|
|---|
| 69 | scalewz 1 3300 120.5 1
|
|---|
| 70 | # scalewz 0.7 2488 102 1
|
|---|
| 71 | convpk2t21
|
|---|
| 72 | cct21h70 = $cct21*$h70cube
|
|---|
| 73 | n/plot hpkz1p0.val*${cct21h70}%$kk ! ! "notit nsta connectpoints black $xyl line=solid,2 "
|
|---|
| 74 | # n/plot hpkz0p7.val*${cct21h70}%$kk ! ! "notit nsta connectpoints black $xyl line=solid,2 "
|
|---|
| 75 |
|
|---|
| 76 | addtext 0.01 4000 '[ P(k)-LSS z=1.0 ]' 'font=helvetica,bolditalic,16 black'
|
|---|
| 77 | setaxelabels 'k_comov (h70 Mpc^-1)' 'P21(k) mK^2 x (Mpc/h70)^3' 'font=helvetica,bolditalic,16'
|
|---|
| 78 |
|
|---|
| 79 | endscript
|
|---|
| 80 |
|
|---|
| 81 | setaxesatt 'minorticks font=helvetica,bold,16 autofontsize'
|
|---|
| 82 |
|
|---|
| 83 | Rad2Deg = 180/3.141596
|
|---|
| 84 | R2D2 = $Rad2Deg*$Rad2Deg
|
|---|
| 85 | Tsys = 50
|
|---|
| 86 | Delnu = 1.e6
|
|---|
| 87 |
|
|---|
| 88 | nu = $nu21/(1+$z)
|
|---|
| 89 | pi23 = 8.*Pi*Pi*Pi
|
|---|
| 90 |
|
|---|
| 91 | Lambda = 0.21*(1+$z)
|
|---|
| 92 | Lam2 = $Lambda*$Lambda
|
|---|
| 93 |
|
|---|
| 94 | Dx = 100
|
|---|
| 95 | Dy = 100
|
|---|
| 96 |
|
|---|
| 97 | defscript calcul
|
|---|
| 98 | Dxol = $Dx/$Lambda
|
|---|
| 99 | Dyol = $Dy/$Lambda
|
|---|
| 100 |
|
|---|
| 101 | # FOV = (1.2*1.2*$Lam2/$Dx/$Dy)
|
|---|
| 102 | FOV = ($Lam2/$Dx/$Dy)
|
|---|
| 103 |
|
|---|
| 104 | FOVDEG = $FOV*$R2D2
|
|---|
| 105 | NPointing = 10000/$FOVDEG
|
|---|
| 106 | tinteg = 365*86400/$NPointing
|
|---|
| 107 | PNOISE = 2.*$Tsys*$Tsys/$tinteg/$Dxol/$Dyol
|
|---|
| 108 | PNOISE = $PNOISE*$Da*$Da*$c/$H*(1+$z)/$nu
|
|---|
| 109 | # PNOISE = $PNOISE*$Da*$Da*$c/$H*(1+$z)/$nu/$pi23
|
|---|
| 110 | PNOISE = $PNOISE*1.e6
|
|---|
| 111 | PNOISE = $PNOISE*1.05
|
|---|
| 112 | echo " FOV = $FOV deg^2 NPointing= $NPointing"
|
|---|
| 113 | echo " tinteg= $tinteg sec PNOISE= $PNOISE mK^2"
|
|---|
| 114 | endscript
|
|---|
| 115 |
|
|---|
| 116 |
|
|---|
| 117 | defscript plnoiseDishes
|
|---|
| 118 | Dx = 75
|
|---|
| 119 | Dy = 75
|
|---|
| 120 | calcul
|
|---|
| 121 | plot2d pknD75m k/$Da pnoise*$PNOISE/5 fracmodok>0.2 'same line=solid,2 cpts red nsta notit'
|
|---|
| 122 | Dx = 100
|
|---|
| 123 | Dy = 100
|
|---|
| 124 | calcul
|
|---|
| 125 | plot2d pknD100m k/$Da pnoise*$PNOISE/20 fracmodok>0.2 'same line=solid,2 cpts blueviolet nsta notit'
|
|---|
| 126 |
|
|---|
| 127 | plot2d pknD100m k/$Da pnoise*$PNOISE/100 fracmodok>0.2 'same line=solid,2 cpts blue nsta notit'
|
|---|
| 128 |
|
|---|
| 129 |
|
|---|
| 130 |
|
|---|
| 131 | set lines ( '(a) 75m Dish, 5 beams' '(b) 100m Dish, 10 beams' '(b) 100m Dish, 100 beams' )
|
|---|
| 132 | set cols ( red blueviolet blue )
|
|---|
| 133 | textdrawer lines cols 'frame font=helvetica,bold,16 inset=0.1,0.3,0.15,0.35'
|
|---|
| 134 | settitle ' PNoise(k) : 5/10/100 beams/polar @z=1' ' ' 'font=helvetica,bold,16'
|
|---|
| 135 | endscript
|
|---|
| 136 |
|
|---|
| 137 | defscript plnoiseEFGH
|
|---|
| 138 | Dx = 12*0.9
|
|---|
| 139 | Dy = 0.844248*0.8
|
|---|
| 140 | calcul
|
|---|
| 141 | plot2d pknF4cyl k/$Da pnoise*$PNOISE fracmodok>0.2 'same line=solid,2 cpts forestgreen nsta notit'
|
|---|
| 142 | plot2d pknF4cylp k/$Da pnoise*$PNOISE fracmodok>0.2 'same line=solid,2 cpts green nsta notit'
|
|---|
| 143 |
|
|---|
| 144 | plot2d pknG8cyl k/$Da pnoise*$PNOISE fracmodok>0.2 'same line=solid,2 cpts orangered nsta notit'
|
|---|
| 145 | plot2d pknG8cylp k/$Da pnoise*$PNOISE fracmodok>0.2 'same line=solid,2 cpts siennared nsta notit'
|
|---|
| 146 |
|
|---|
| 147 | Dx = 5*0.9
|
|---|
| 148 | Dy = 5*0.9
|
|---|
| 149 | calcul
|
|---|
| 150 | plot2d pknE400d k/$Da pnoise*$PNOISE fracmodok>0.2 'same line=solid,2 cpts darkviolet nsta notit'
|
|---|
| 151 |
|
|---|
| 152 | set lines ( '(e) 20x20:400xD=5m' '(f) 4Cyl-12mx85m, 400 rec/pol' '(fp) 4Cylp-12mx85m, 400 rec/pol' '(g) 8Cyl-12mx105m, 960 rec/pol' '(gp) 8Cylp-12mx105m, 960 rec/pol' )
|
|---|
| 153 | set cols ( darkviolet forestgreen green orangered siennared )
|
|---|
| 154 | textdrawer lines cols 'frame font=helvetica,bold,16 inset=0.1,0.3,0.15,0.35'
|
|---|
| 155 | settitle ' PNoise(k) : Dishes/Cylinders, 400/400/960 recv/pol @z=1' ' ' 'font=helvetica,bold,16'
|
|---|
| 156 | endscript
|
|---|
| 157 |
|
|---|
| 158 | defscript plnoiseABCD
|
|---|
| 159 | Dx = 5*0.9
|
|---|
| 160 | Dy = 5*0.9
|
|---|
| 161 | calcul
|
|---|
| 162 | plot2d pknA121d k/$Da pnoise*$PNOISE fracmodok>0.2 'same line=solid,2 cpts magenta nsta notit'
|
|---|
| 163 | plot2d pknB128d k/$Da pnoise*$PNOISE fracmodok>0.2 'same line=solid,2 cpts red nsta notit'
|
|---|
| 164 | plot2d pknC129d k/$Da pnoise*$PNOISE fracmodok>0.2 'same line=solid,2 cpts orange nsta notit'
|
|---|
| 165 |
|
|---|
| 166 | Dx = 75
|
|---|
| 167 | Dy = 75
|
|---|
| 168 | calcul
|
|---|
| 169 | plot2d pknD75m k/$Da pnoise*$PNOISE/100 fracmodok>0.2 'same line=solid,2 cpts blueviolet nsta notit'
|
|---|
| 170 |
|
|---|
| 171 | set lines ( '(a) 11x11:121xD=5m' '(b) 128xD=5m' '(c) 129xD=5m' '(d) 75m Dish, 100 beams' )
|
|---|
| 172 | set cols ( magenta red orange blueviolet )
|
|---|
| 173 | textdrawer lines cols 'frame font=helvetica,bold,16 inset=0.1,0.3,0.15,0.35'
|
|---|
| 174 | settitle ' PNoise(k) : 121,128,129,100 receiver/polar @z=1' ' ' 'font=helvetica,bold,16'
|
|---|
| 175 | endscript
|
|---|
| 176 |
|
|---|
| 177 | defscript plfracABCD
|
|---|
| 178 | xylf = "xylimits=0.005,0.5,0.,1. minorticks"
|
|---|
| 179 | plot2d pknA121d k/$Da fracmodok fracmodok>0. '$xylf line=solid,2 cpts magenta nsta notit'
|
|---|
| 180 | plot2d pknB128d k/$Da fracmodok fracmodok>0. 'same line=solid,2 cpts red nsta notit'
|
|---|
| 181 | plot2d pknC129d k/$Da fracmodok fracmodok>0. 'same line=solid,2 cpts orange nsta notit'
|
|---|
| 182 | plot2d pknD75m k/$Da fracmodok fracmodok>0. 'same line=solid,2 cpts violetblue nsta notit'
|
|---|
| 183 | endscript
|
|---|
| 184 |
|
|---|
| 185 | defscript plfracEFGH
|
|---|
| 186 | xylf = "xylimits=0.005,0.5,0.,1. minorticks"
|
|---|
| 187 | plot2d pknF4cyl k/$Da fracmodok fracmodok>0. '$xylf line=solid,2 cpts forestgreen nsta notit'
|
|---|
| 188 | plot2d pknF4cylp k/$Da fracmodok fracmodok>0. 'same line=solid,2 cpts green nsta notit'
|
|---|
| 189 | plot2d pknG8cyl k/$Da fracmodok fracmodok>0. 'same line=solid,2 cpts gold nsta notit'
|
|---|
| 190 | plot2d pknG8cylp k/$Da fracmodok fracmodok>0. 'same line=solid,2 cpts yellow nsta notit'
|
|---|
| 191 | plot2d pknE400d k/$Da fracmodok fracmodok>0. 'same line=solid,2 cpts magenta nsta notit'
|
|---|
| 192 |
|
|---|
| 193 | endscript
|
|---|
| 194 |
|
|---|
| 195 | defscript pldishes
|
|---|
| 196 |
|
|---|
| 197 | Dx = 50
|
|---|
| 198 | Dy = 50
|
|---|
| 199 | calcul
|
|---|
| 200 | plot2d pknD50m k/$Da pnoise*$PNOISE/100 fracmodok>0.2 'same line=solid,2 cpts skyblue nsta notit'
|
|---|
| 201 |
|
|---|
| 202 | Dx = 75
|
|---|
| 203 | Dy = 75
|
|---|
| 204 | calcul
|
|---|
| 205 | plot2d pknD75m k/$Da pnoise*$PNOISE/100 fracmodok>0.2 'same line=solid,2 cpts blue nsta notit'
|
|---|
| 206 |
|
|---|
| 207 | Dx = 100
|
|---|
| 208 | Dy = 100
|
|---|
| 209 | calcul
|
|---|
| 210 | plot2d pknD100m k/$Da pnoise*$PNOISE/100 fracmodok>0.2 'same line=solid,2 cpts navyblue nsta notit'
|
|---|
| 211 |
|
|---|
| 212 | Dx = 200
|
|---|
| 213 | Dy = 200
|
|---|
| 214 | calcul
|
|---|
| 215 | plot2d pknD200m k/$Da pnoise*$PNOISE/100 fracmodok>0.2 'same line=solid,2 cpts blueviolet nsta notit'
|
|---|
| 216 |
|
|---|
| 217 | set lines ( 'D:Dish50m' 'D:Dish75m' 'D:Dish100m' 'D:Dish200m' )
|
|---|
| 218 | set cols ( skyblue blue navyblue blueviolet )
|
|---|
| 219 | textdrawer lines cols 'frame font=helvetica,bold,16 inset=0.1,0.3,0.15,0.35'
|
|---|
| 220 | settitle ' PNoise(k) : Dish D=50m,75m,100m,200m' ' ' 'font=helvetica,bold,16'
|
|---|
| 221 | endscript
|
|---|
| 222 |
|
|---|
| 223 |
|
|---|
| 224 |
|
|---|
| 225 | defscript plnoisenancay
|
|---|
| 226 | # Dx = 5*0.95
|
|---|
| 227 | # Dy = 5*0.95
|
|---|
| 228 | # calcul
|
|---|
| 229 | # plot2d pknnan24 x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts red nsta notit'
|
|---|
| 230 | # plot2d pknnan25 x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts orange nsta notit'
|
|---|
| 231 | # Dx = 3.5*0.95
|
|---|
| 232 | # Dy = 3.5*0.95
|
|---|
| 233 | # calcul
|
|---|
| 234 | # plot2d pknnan36 x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts red nsta notit'
|
|---|
| 235 | # plot2d pknnan40 x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts violet nsta notit'
|
|---|
| 236 | # Dx = 7*0.9
|
|---|
| 237 | # Dy = 2*$Lambda*0.9
|
|---|
| 238 | # calcul
|
|---|
| 239 | # plot2d pknpit2cyl x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts green nsta notit'
|
|---|
| 240 | Dx = 9*0.9
|
|---|
| 241 | Dy = 2*$Lambda*0.9
|
|---|
| 242 | calcul
|
|---|
| 243 | plot2d pknpit2cylw x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts red nsta notit'
|
|---|
| 244 | set lines ( '2 Cylinders 20mx9m' )
|
|---|
| 245 | set cols ( red )
|
|---|
| 246 | textdrawer lines cols 'noframe font=helvetica,bold,16 inset=0.4,0.4,0.4,0.4'
|
|---|
| 247 |
|
|---|
| 248 | endscript
|
|---|
| 249 |
|
|---|
| 250 | defscript ABCD
|
|---|
| 251 | # plnoisedish
|
|---|
| 252 | # zone 1 2
|
|---|
| 253 | # newwin 1 1 1000 600
|
|---|
| 254 | plpklss
|
|---|
| 255 | plnoiseABCD
|
|---|
| 256 | # plfracABCD
|
|---|
| 257 | endscript
|
|---|
| 258 |
|
|---|
| 259 | defscript EFGH
|
|---|
| 260 | # plnoisedish
|
|---|
| 261 | # zone 1 2
|
|---|
| 262 | # newwin 1 1 900 600
|
|---|
| 263 | plpklss
|
|---|
| 264 | plnoiseEFGH
|
|---|
| 265 | # plfracEFGH
|
|---|
| 266 | endscript
|
|---|
| 267 |
|
|---|
| 268 | defscript DDDD
|
|---|
| 269 | # plnoisedish
|
|---|
| 270 | # zone 1 2
|
|---|
| 271 | plpklss
|
|---|
| 272 | plnoiseDishes
|
|---|
| 273 | # plfracABCD
|
|---|
| 274 | endscript
|
|---|
| 275 |
|
|---|
| 276 | defscript BB
|
|---|
| 277 | plnoisedish
|
|---|
| 278 | plpklss
|
|---|
| 279 | plnoiseB
|
|---|
| 280 | endscript
|
|---|
| 281 |
|
|---|
| 282 | defscript CC
|
|---|
| 283 | plnoisedish
|
|---|
| 284 | plpklss
|
|---|
| 285 | plnoiseC2
|
|---|
| 286 | endscript
|
|---|
| 287 |
|
|---|
| 288 | defscript nancay
|
|---|
| 289 | plnoisedish
|
|---|
| 290 | plpklss
|
|---|
| 291 | plnoisenancay
|
|---|
| 292 | endscript
|
|---|
| 293 |
|
|---|
| 294 | defscript POSCOV
|
|---|
| 295 | openppf hdt_repnan128.ppf
|
|---|
| 296 | rename mdish mdB128d
|
|---|
| 297 | rename h2rep uvB128d
|
|---|
| 298 |
|
|---|
| 299 | openppf hdt_repf11x11.ppf
|
|---|
| 300 | rename mdish mdA121d
|
|---|
| 301 | rename h2rep uvA121d
|
|---|
| 302 |
|
|---|
| 303 | openppf hdt_repconfC.ppf
|
|---|
| 304 | rename mdish mdC129d
|
|---|
| 305 | rename h2rep uvC129d
|
|---|
| 306 |
|
|---|
| 307 | newwin 2 1 800 400
|
|---|
| 308 | zone 2 1
|
|---|
| 309 | nt2d mdB128d posx posy - - - - 'xylimits=-10,90,-10,90 marker=circle,15 notit nsta red '
|
|---|
| 310 | setaxelabels ' X (meters) ' ' Y (meters) ' 'font=helvetica,bolditalic,16'
|
|---|
| 311 | settitle '(b) 128 D=5m dishes in 8 rows ' ' ' 'font=helvetica,bold,16'
|
|---|
| 312 | nt2d mdC129d posx posy - - - - 'xylimits=-10,90,-10,90 marker=circle,15 notit nsta red'
|
|---|
| 313 | setaxelabels ' X (meters) ' ' Y (meters) ' 'font=helvetica,bolditalic,16'
|
|---|
| 314 | settitle '(c) 129 D=5m dishes ' ' ' 'font=helvetica,bold,16'
|
|---|
| 315 | pssetfilename configab.ps
|
|---|
| 316 | w2ps
|
|---|
| 317 | w2eps configab.eps
|
|---|
| 318 | psclosefile
|
|---|
| 319 |
|
|---|
| 320 | newwin 2 2 800 800
|
|---|
| 321 | disp uvA121d 'h2disp=img colbr128 h2dyn=1,80 notit nsta'
|
|---|
| 322 | setaxelabels 'u (Radian^-1) u=1000->~21 arcmin' ' v (Radian^-1) ' 'font=helvetica,bolditalic,12'
|
|---|
| 323 | settitle 'u-v coverage, (a) 11x11 D=5m dishes Over 55mx55m' ' ' 'font=helvetica,bold,12'
|
|---|
| 324 | disp uvB128d 'h2disp=img colbr128 h2dyn=1,80 notit nsta'
|
|---|
| 325 | setaxelabels 'u (Radian^-1) u=1000->~21 arcmin' ' v (Radian^-1) ' 'font=helvetica,bolditalic,12'
|
|---|
| 326 | settitle 'u-v coverage, (b) 8 row of 16xD=5m dishes Over 80mx80m' ' ' 'font=helvetica,bold,12'
|
|---|
| 327 | disp uvC129d 'h2disp=img colbr128 h2dyn=1,80 notit nsta'
|
|---|
| 328 | setaxelabels 'u (Radian^-1) u=1000->~21 arcmin' ' v (Radian^-1) ' 'font=helvetica,bolditalic,12'
|
|---|
| 329 | settitle 'u-v coverage, (c) 129 D=5m dishes Over 80mx80m' ' ' 'font=helvetica,bold,12'
|
|---|
| 330 | disp h2drD75m 'h2disp=img colbr128 h2dyn=0.02,0.8 notit nsta'
|
|---|
| 331 | setaxelabels 'u (Radian^-1) u=1000->~21 arcmin' ' v (Radian^-1) ' 'font=helvetica,bolditalic,12'
|
|---|
| 332 | settitle 'u-v coverage, (d) D=75 m dish with 100 beams' ' ' 'font=helvetica,bold,12'
|
|---|
| 333 |
|
|---|
| 334 | pssetfilename uvcovabcd.ps
|
|---|
| 335 | w2ps
|
|---|
| 336 | w2eps uvcovabcd.eps
|
|---|
| 337 | psclosefile
|
|---|
| 338 |
|
|---|
| 339 | # imag uvB128d 'colbr128 stdaxes showcmap=right'
|
|---|
| 340 |
|
|---|
| 341 | endscript
|
|---|
![(please configure the [header_logo] section in trac.ini)](/PI/chrome/site/your_project_logo.png){kind=logo}
