Changeset 3973 in Sophya for trunk/Cosmo/RadioBeam/plpkn.pic


Ignore:
Timestamp:
Apr 18, 2011, 5:30:44 PM (14 years ago)
Author:
ansari
Message:

Corrections diverses: choix lobe gaussien/triangle et specif DishDiameter au lieu de DoL ds applobe/calcpk2, possibilite application lobe freq.independante ds applobe, Reza 18/04/2011

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/Cosmo/RadioBeam/plpkn.pic

    r3947 r3973  
    22### Script de trace de PNoise(k) et reponse dans le plan (u,v)
    33###  de l'interferometre a partir du fichier PPF produit par pknoise.cc
    4 ###     Fev - Avril 2010 ,  BAORadio/Reza
     4###    Janvier 2011 ,  BAORadio/Reza
    55########################################################################
    66
    7 if ( $# < 1 ) then
    8   echo ' Usage: exec plpkn PPFName_pknoise'
    9   return
    10 endif
    11 
    12 echo "---> openppf $1 "
    13 openppf $1
     7echo ' -----> plpknew.pic : opening PPF files ... "
     8delobjs *
     9# openppf ../../PkNoise/cmvhpkz.ppf
     10# openppf ../pkz0p25.ppf
     11openppf ../pkz0p7.ppf
     12rename hpkz hpkz0p7
     13openppf ../pkz1p0.ppf
     14rename hpkz hpkz1p0
     15
     16defscript 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
     26endscript
     27
     28pknopen pknf11x11.ppf A121d
     29pknopen pknnan128.ppf B128d
     30pknopen pknconfC.ppf  C129d
     31pknopen pknD50m.ppf   D50m
     32pknopen pknD75m.ppf   D75m
     33pknopen pknD100m.ppf   D100m
     34pknopen pknD200m.ppf   D200m
     35pknopen pknf20x20.ppf E400d
     36pknopen pknf4cyl.ppf F4cyl
     37pknopen pknf4cylp.ppf F4cylp
     38pknopen pknf8cyl.ppf G8cyl
     39pknopen pknf8cylp.ppf G8cylp
     40
     41
     42listobjs
     43
     44#  A z = 0.7
     45z = 0.7
     46c = 3.e5
     47H = 102
     48Da = 2488
     49nu21 = 1.42e9
    1450
    1551echo '----> Executing anapkn.pic'
    1652exec anapkn.pic
    1753setup5
    18 scalewz 0.7 2500 100 1
     54# scalewz 0.7 2488 102 1
    1955# scalewz 0.25 989.13 80.26 1
    20 y1 = 500*$cct21
    21 y2 = 9e4*$cct21
    22 xyl = "xylimits=0.01,0.5,$y1,$y2 logx logy minorticks"
     56scalewz 1 3300 120.5 1
     57Da = 3300
     58
     59
     60y1 = 400*$cct21
     61y2 = 7e4*$cct21
     62xyl = "xylimits=0.005,0.5,$y1,$y2 logx logy minorticks"
     63# H = 71.9 km/s/Mpc = 1.0271 x 70 km/s
     64h70 = 1.0271
     65h70cube = $h70*$h70*$h70
     66set kk pow(10.,x)/$h70
     67
    2368defscript plpklss
    24   n/plot hpkz.val*${cct21}%$kk ! ! "same notit nsta connectpoints black $xyl line=solid,2"
    25   addtext 0.03 2000 '*** P(k)-LSS ***' 'font=helvetica,bolditalic,16 black'
     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
    2679endscript 
    2780
    … …  
    3386Delnu = 1.e6
    3487
    35 #  A z = 0.7
    36 z = 0.7
    37 c = 3.e5
    38 H = 100
    39 Da = 2500
    40 nu21 = 1.42e9
    4188nu = $nu21/(1+$z)
    4289pi23 = 8.*Pi*Pi*Pi
    4390
    44 Lambda = 0.357
     91Lambda = 0.21*(1+$z)
    4592Lam2 = $Lambda*$Lambda
    4693
    … …  
    5299  Dyol = $Dy/$Lambda
    53100
    54   FOV = (1.2*1.2*$Lam2/$Dx/$Dy)
     101#   FOV = (1.2*1.2*$Lam2/$Dx/$Dy)
     102  FOV = ($Lam2/$Dx/$Dy)
     103
    55104  FOVDEG = $FOV*$R2D2
    56105  NPointing = 10000/$FOVDEG
    57106  tinteg = 365*86400/$NPointing
    58   PNOISE = $Tsys*$Tsys/$tinteg/$Dxol/$Dyol
     107  PNOISE = 2.*$Tsys*$Tsys/$tinteg/$Dxol/$Dyol
    59108  PNOISE = $PNOISE*$Da*$Da*$c/$H*(1+$z)/$nu
    60109#   PNOISE = $PNOISE*$Da*$Da*$c/$H*(1+$z)/$nu/$pi23
    61110  PNOISE = $PNOISE*1.e6
     111  PNOISE = $PNOISE*1.05
    62112  echo " FOV = $FOV deg^2  NPointing= $NPointing"
    63113  echo " tinteg= $tinteg sec  PNOISE= $PNOISE mK^2"
    64114endscript
    65115
    66 defscript plnoisedish
     116
     117defscript 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'
    67122  Dx = 100
    68123  Dy = 100
    69124  calcul
    70   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'
     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'
     135endscript
     136
     137defscript 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 violetred nsta notit'
     145  plot2d pknG8cylp k/$Da pnoise*$PNOISE fracmodok>0.2 'same line=solid,2 cpts violet 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 turquoise 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 ( turquoise forestgreen green violetred violet  )
     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'
     156endscript
     157
     158defscript 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'
     175endscript
     176
     177defscript 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'
     183endscript
     184
     185defscript 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
     193endscript
     194
     195defscript 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
    71212  Dx = 200
    72213  Dy = 200
    73214  calcul
    74   plot2d noiseD2 x/$Da val*$PNOISE nb>10 'same line=solid,2  cpts grey nsta notit'
    75   set lines ( '100mDish' '200mDish'  )
    76   set cols ( navyblue grey  )
    77   textdrawer lines cols 'frame font=helvetica,bold,16 inset=0.1,0.3,0.7,0.8'
    78   setaxelabels 'k (Mpc^-1) ' 'PNoise(k) mk^2 Mpc^3' 'font=helvetica,bolditalic,16'
    79 endscript
    80 
    81 defscript plnoiseA
    82   Dx = 5*0.95
    83   Dy = 5*0.95
    84   calcul
    85   plot2d noisemdf64 x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts red nsta notit'
    86 #  plot2d noisemds x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts red nsta notit'
    87   plot2d noisemdsB x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts siennared nsta notit'
    88   plot2d noisemdsC x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts violetred nsta notit'
    89   Dx = 25*0.3
    90   Dy = 0.5*0.9
    91   calcul
    92   plot2d noise2cyl x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts forestgreen nsta notit'
    93   plot2d noise2cylP x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts green nsta notit'
    94   set lines (  'FilledA:64x5mD' 'SparseB:72x5mD' 'SparseC:129x5mD' 'Pitts2Cyl(=64C)' 'PerfPitts2Cyl(=64C)' )
    95   set cols ( red siennared magenta  forestgreen green )
     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 )
    96219  textdrawer lines cols 'frame font=helvetica,bold,16  inset=0.1,0.3,0.15,0.35'
    97   settitle ' PNoise(k) : Dishes/Cylinders, 64/72/129 channels' ' ' 'font=helvetica,bold,16'
    98 endscript
    99 
    100 defscript plnoiseB
    101   Dx = 5*0.95
    102   Dy = 5*0.95
    103   calcul
    104   plot2d noisemdsC x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts red nsta notit'
    105   plot2d noisemdf x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts gold nsta notit'
    106   Dx = 10*0.95
    107   Dy = 0.5*0.9
    108   calcul
    109   plot2d noise3cyl x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts violet nsta notit'
    110   plot2d noise3cylP x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts magenta nsta notit'
    111 
    112   Dx = 12*0.95
    113   Dy = 0.5*0.9
    114   calcul
    115   plot2d noisefcyl x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts blue nsta notit'
    116   plot2d noisefcylP x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts royalblue nsta notit'
    117 
    118   set lines ( 'SparseC:129x5mD' 'Filled400x5m' '3xCyl10x64m(=384C)' 'Perf3xCyl10x64m(=384C)' )
    119   set lines ( $lines '8Cyl12x96m(=1536C)' 'Perf8Cyl12x96m(=1536C)' )
    120   set cols ( red gold  violet magenta  blue royalblue )
    121   textdrawer lines cols 'frame font=helvetica,bold,16 inset=0.25,0.55,0.6,0.95'
    122   settitle ' PNoise(k) : Dishes/Cylinders, 129/384/400/1536 channels' ' ' 'font=helvetica,bold,16'     
    123 endscript
    124 
    125 defscript xxx
    126 newwin
    127 disp mfill 'h2disp=img colbr128 notit'
    128 setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16'
    129 settitle ' u-v coverage , 400 x 5m Dishes - No Pointing' ' ' 'font=helvetica,bold,16'
    130 
    131 disp dish 'h2disp=img colbr128 notit'
    132 setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16'
    133 settitle ' u-v coverage , 100 m Dish' ' ' 'font=helvetica,bold,16'
    134 
    135 disp msparsfp 'h2disp=img colbr128 notit'
    136 setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16'
    137 settitle ' u-v coverage , 63 x 5m Dishes T-config - No Pointing' ' ' 'font=helvetica,bold,16'
    138 
    139 disp mspars 'h2disp=img colbr128 notit'
    140 setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16'
    141 settitle ' u-v coverage , 63 x 5m Dishes T-config - Pointing ~Pi/4' ' ' 'font=helvetica,bold,16'
    142 
    143 disp mcylf 'h2disp=img colbr128 notit'
    144 setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16'
    145 settitle ' u-v coverage , Filled Cylinder Array 8 Cyl 12mx96m ' ' ' 'font=helvetica,bold,16'
    146 
    147 disp mcylfP 'h2disp=img colbr128 notit'
    148 setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16'
    149 settitle ' u-v coverage , Perfect Filled Cylinder Array 8 Cyl 12mx96m ' ' ' 'font=helvetica,bold,16'
    150 
    151 disp m2cyl 'h2disp=img colbr128 notit'
    152 setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16'
    153 settitle ' u-v coverage , Pittsburgh 2 Cyl 16mx8m , dist=25m ' ' ' 'font=helvetica,bold,16'
    154 
    155 disp m2cylP 'h2disp=img colbr128 notit'
    156 setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16'
    157 settitle ' u-v coverage , Perfect-Pitts. 2 Cyl 16mx8m , dist=25m ' ' ' 'font=helvetica,bold,16'
    158 
    159 
    160 disp noiseD 'logy nsta'
    161 disp noiseD2 'same grey nsta'
    162 disp noisemdf 'same cyan nsta'
    163 disp noisemdf64 'same brown nsta'
    164 disp noisemds 'same red nsta'
    165 disp noisemdsC 'same orange nsta'
    166 
    167 disp noisemdsB 'same yellow nsta'
    168 disp noisemdsfp 'same yellow nsta'
    169 # disp noisemdsd7 'same gold nsta'
    170 disp noisefcyl 'same blue nsta'
    171 disp noisefcylP 'same skyblue nsta'
    172 disp noise3cylP 'same magenta nsta'
    173 disp noise3cyl 'same violet nsta'
    174 disp noise2cyl 'same forestgreen nsta'
    175 disp noise2cylP 'same green nsta'
    176 
    177 endscript
    178 
    179 defscript AA
    180   plnoisedish
    181   plpklss
    182   plnoiseA
     220  settitle ' PNoise(k) : Dish D=50m,75m,100m,200m' ' ' 'font=helvetica,bold,16'
     221endscript
     222
     223
     224
     225defscript 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
     248endscript
     249
     250defscript ABCD
     251#  plnoisedish
     252#  zone 1 2
     253#  newwin 1 1 1000 600
     254  plpklss
     255  plnoiseABCD
     256#  plfracABCD
     257endscript
     258
     259defscript EFGH
     260#  plnoisedish
     261#  zone 1 2
     262#  newwin 1 1 900 600
     263  plpklss
     264  plnoiseEFGH
     265#  plfracEFGH
     266endscript
     267
     268defscript DDDD
     269#  plnoisedish
     270#  zone 1 2
     271  plpklss
     272  plnoiseDishes
     273#  plfracABCD
    183274endscript
    184275
    … …  
    189280endscript
    190281
     282defscript CC
     283  plnoisedish
     284  plpklss
     285  plnoiseC2
     286endscript
     287
     288defscript nancay
     289  plnoisedish
     290  plpklss
     291  plnoisenancay
     292endscript
     293
    191294defscript POSCOV
    192   disp posspB red
     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 '
    193310  setaxelabels ' X (meters) ' ' Y (meters) ' 'font=helvetica,bolditalic,16'
    194   settitle ' Config B dish positions - 72 dishes  ' ' ' 'font=helvetica,bold,16'
    195 #  w2ps
    196    
    197   disp posspC red
     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'
    198313  setaxelabels ' X (meters) ' ' Y (meters) ' 'font=helvetica,bolditalic,16'
    199   settitle ' Config C dish positions - 129 dishes  ' ' ' 'font=helvetica,bold,16'
    200 #  w2ps
    201   disp mfill64 'h2disp=img colbr128 notit'
    202   setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16'
    203   settitle ' u-v coverage , Filled 8x8 - 64 x 5m Dishes' ' ' 'font=helvetica,bold,16'
    204 #  w2ps
    205   disp m3cyl 'h2disp=img colbr128 notit'
    206   setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16'
    207   settitle ' u-v coverage , 3 Cylinders 10mx64m ' ' ' 'font=helvetica,bold,16'
    208 #  w2ps
    209   disp mfill 'h2disp=img colbr128 notit'
    210   setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16'
    211   settitle ' u-v coverage , Filled 20x20 - 400 x 5m Dishes' ' ' 'font=helvetica,bold,16'
    212 #  w2ps
    213   disp msparsC 'h2disp=img colbr128 notit'
    214   setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16'
    215   settitle 'u-v coverage, Sparse-C: 129x5mD Over 80mx80m (Rot~Pi/4)' ' ' 'font=helvetica,bold,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'
    216340 
    217341endscript 
Note: See TracChangeset for help on using the changeset viewer.