![(please configure the [header_logo] section in trac.ini)](/PI/chrome/site/your_project_logo.png){kind=logo}
Changeset 3769 in Sophya for trunk/Cosmo/RadioBeam/plpkn.pic
- Timestamp:
- May 7, 2010, 6:44:43 PM (15 years ago)
- File:
-
- 1 edited
-
trunk/Cosmo/RadioBeam/plpkn.pic (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/Cosmo/RadioBeam/plpkn.pic
r3756 r3769 5 5 ######################################################################## 6 6 7 openppf pknoise.ppf 7 if ( $# < 1 ) then 8 echo ' Usage: exec plpkn PPFName_pknoise' 9 return 10 endif 11 12 echo "---> openppf $1 " 13 openppf $1 14 15 echo '----> Executing anapkn.pic' 16 exec anapkn.pic 17 setup5 18 scalewz 0.7 2500 100 1 19 y1 = 500*$cct21 20 y2 = 9e4*$cct21 21 xyl = "xylimits=0.01,0.5,$y1,$y2 logx logy minorticks" 22 defscript plpklss 23 n/plot hpkz.val*${cct21}%$kk ! ! "same notit nsta connectpoints black $xyl line=solid,2" 24 addtext 0.03 2000 '*** P(k)-LSS ***' 'font=helvetica,bolditalic,16 black' 25 endscript 8 26 9 27 setaxesatt 'minorticks font=helvetica,bold,16 autofontsize' 10 disp noiseD 'logy nsta'11 disp noiseD2 'same grey nsta'12 disp noisemdf 'same orange nsta'13 disp noisemds 'same red nsta'14 disp noisemdsfp 'same yellow nsta'15 disp noisemdsd7 'same gold nsta'16 disp noisefcyl 'same blue nsta'17 disp noisefcylP 'same skyblue nsta'18 disp noise2cyl 'same forestgreen nsta'19 disp noise2cylP 'same green nsta'20 28 21 29 Rad2Deg = 180/3.141596 … … 55 63 endscript 56 64 57 Dx = 100 58 Dy = 100 59 calcul 60 plot2d noiseD x/$Da val*$PNOISE nb>10 'line=solid,2 logy logx xylimits=0.002,0.8,1,1e5 grid cpts nsta notit' 61 Dx = 200 62 Dy = 200 63 calcul 64 plot2d noiseD2 x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts grey nsta notit' 65 Dx = 5*0.95 66 Dy = 5*0.95 67 calcul 68 plot2d noisemdf x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts orange nsta notit' 69 plot2d noisemds x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts red nsta notit' 70 Dx = 7.55*0.95 71 Dy = 7.5*0.95 72 calcul 73 plot2d noisemdsd7 x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts siennared nsta notit' 74 Dx = 12*0.9 75 Dy = 0.5*0.9 76 calcul 77 plot2d noisefcyl x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts blue nsta notit' 78 plot2d noisefcylP x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts skyblue nsta notit' 79 Dx = 8*0.9 80 Dy = 0.5*0.9 81 calcul 82 plot2d noise2cyl x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts forestgreen nsta notit' 83 plot2d noise2cylP x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts green nsta notit' 84 85 set lines ( '100mDish' '200mDish' '400x5mDish' '63x5mDish' '63x7.5mDish' '8Cyl12x96m' 'Perf8Cyl12x96m' ) 86 set lines ( $lines 'Pitts2Cyl' 'PerfPitts2Cyl' ) 87 set cols ( black grey orange red siennared blue skyblue forestgreen green ) 88 textdrawer lines cols 'frame font=helvetica,bold,16' 89 setaxelabels 'k (Mpc^-1) ' 'PNoise(k) mk^2 Mpc^3' 'font=helvetica,bolditalic,16' 90 91 65 defscript plnoisedish 66 Dx = 100 67 Dy = 100 68 calcul 69 plot2d noiseD x/$Da val*$PNOISE nb>10 'line=solid,2 logy logx xylimits=0.002,0.8,1,1e5 navyblue grid cpts nsta notit' 70 Dx = 200 71 Dy = 200 72 calcul 73 plot2d noiseD2 x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts grey nsta notit' 74 set lines ( '100mDish' '200mDish' ) 75 set cols ( navyblue grey ) 76 textdrawer lines cols 'frame font=helvetica,bold,16 inset=0.1,0.3,0.7,0.8' 77 setaxelabels 'k (Mpc^-1) ' 'PNoise(k) mk^2 Mpc^3' 'font=helvetica,bolditalic,16' 78 endscript 79 80 defscript plnoiseA 81 Dx = 5*0.95 82 Dy = 5*0.95 83 calcul 84 plot2d noisemdf64 x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts red nsta notit' 85 # plot2d noisemds x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts red nsta notit' 86 plot2d noisemdsB x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts siennared nsta notit' 87 plot2d noisemdsC x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts violetred nsta notit' 88 Dx = 25*0.3 89 Dy = 0.5*0.9 90 calcul 91 plot2d noise2cyl x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts forestgreen nsta notit' 92 plot2d noise2cylP x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts green nsta notit' 93 set lines ( 'FilledA:64x5mD' 'SparseB:72x5mD' 'SparseC:129x5mD' 'Pitts2Cyl(=64C)' 'PerfPitts2Cyl(=64C)' ) 94 set cols ( red siennared magenta forestgreen green ) 95 textdrawer lines cols 'frame font=helvetica,bold,16 inset=0.1,0.3,0.15,0.35' 96 settitle ' PNoise(k) : Dishes/Cylinders, 64/72/129 channels' ' ' 'font=helvetica,bold,16' 97 endscript 98 99 defscript plnoiseB 100 Dx = 5*0.95 101 Dy = 5*0.95 102 calcul 103 plot2d noisemdsC x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts red nsta notit' 104 plot2d noisemdf x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts gold nsta notit' 105 Dx = 10*0.95 106 Dy = 0.5*0.9 107 calcul 108 plot2d noise3cyl x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts violet nsta notit' 109 plot2d noise3cylP x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts magenta nsta notit' 110 111 Dx = 12*0.95 112 Dy = 0.5*0.9 113 calcul 114 plot2d noisefcyl x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts blue nsta notit' 115 plot2d noisefcylP x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts royalblue nsta notit' 116 117 set lines ( 'SparseC:129x5mD' 'Filled400x5m' '3xCyl10x64m(=384C)' 'Perf3xCyl10x64m(=384C)' ) 118 set lines ( $lines '8Cyl12x96m(=1536C)' 'Perf8Cyl12x96m(=1536C)' ) 119 set cols ( red gold violet magenta blue royalblue ) 120 textdrawer lines cols 'frame font=helvetica,bold,16 inset=0.25,0.55,0.6,0.95' 121 settitle ' PNoise(k) : Dishes/Cylinders, 129/384/400/1536 channels' ' ' 'font=helvetica,bold,16' 122 endscript 123 124 defscript xxx 92 125 newwin 93 126 disp mfill 'h2disp=img colbr128 notit' … … 122 155 setaxelabels 'kx (Radian^-1) k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 123 156 settitle ' u-v coverage , Perfect-Pitts. 2 Cyl 16mx8m , dist=25m ' ' ' 'font=helvetica,bold,16' 157 158 159 disp noiseD 'logy nsta' 160 disp noiseD2 'same grey nsta' 161 disp noisemdf 'same cyan nsta' 162 disp noisemdf64 'same brown nsta' 163 disp noisemds 'same red nsta' 164 disp noisemdsC 'same orange nsta' 165 166 disp noisemdsB 'same yellow nsta' 167 disp noisemdsfp 'same yellow nsta' 168 # disp noisemdsd7 'same gold nsta' 169 disp noisefcyl 'same blue nsta' 170 disp noisefcylP 'same skyblue nsta' 171 disp noise3cylP 'same magenta nsta' 172 disp noise3cyl 'same violet nsta' 173 disp noise2cyl 'same forestgreen nsta' 174 disp noise2cylP 'same green nsta' 175 176 endscript 177 178 defscript AA 179 plnoisedish 180 plpklss 181 plnoiseA 182 endscript 183 184 defscript BB 185 plnoisedish 186 plpklss 187 plnoiseB 188 endscript 189 190 defscript POSCOV 191 disp posspB red 192 setaxelabels ' X (meters) ' ' Y (meters) ' 'font=helvetica,bolditalic,16' 193 settitle ' Config B dish positions - 72 dishes ' ' ' 'font=helvetica,bold,16' 194 # w2ps 195 196 disp posspC red 197 setaxelabels ' X (meters) ' ' Y (meters) ' 'font=helvetica,bolditalic,16' 198 settitle ' Config C dish positions - 129 dishes ' ' ' 'font=helvetica,bold,16' 199 # w2ps 200 disp mfill64 'h2disp=img colbr128 notit' 201 setaxelabels 'kx (Radian^-1) k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 202 settitle ' u-v coverage , Filled 8x8 - 64 x 5m Dishes' ' ' 'font=helvetica,bold,16' 203 # w2ps 204 disp m3cyl 'h2disp=img colbr128 notit' 205 setaxelabels 'kx (Radian^-1) k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 206 settitle ' u-v coverage , 3 Cylinders 10mx64m ' ' ' 'font=helvetica,bold,16' 207 # w2ps 208 disp mfill 'h2disp=img colbr128 notit' 209 setaxelabels 'kx (Radian^-1) k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 210 settitle ' u-v coverage , Filled 20x20 - 400 x 5m Dishes' ' ' 'font=helvetica,bold,16' 211 # w2ps 212 disp msparsC 'h2disp=img colbr128 notit' 213 setaxelabels 'kx (Radian^-1) k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 214 settitle 'u-v coverage, Sparse-C: 129x5mD Over 80mx80m (Rot~Pi/4)' ' ' 'font=helvetica,bold,16' 215 216 endscript
Note:
See TracChangeset
for help on using the changeset viewer.
