########################################################################
### Script de trace de PNoise(k) et reponse dans le plan (u,v) 
###  de l'interferometre a partir du fichier PPF produit par pknoise.cc
###     Fev - Avril 2010 ,  BAORadio/Reza
########################################################################

if ( $# < 1 ) then 
  echo ' Usage: exec plpkn PPFName_pknoise'
  return 
endif

echo "---> openppf $1 "
openppf $1

echo '----> Executing anapkn.pic'
exec anapkn.pic
setup5 
scalewz 0.7 2500 100 1
y1 = 500*$cct21
y2 = 9e4*$cct21
xyl = "xylimits=0.01,0.5,$y1,$y2 logx logy minorticks"
defscript plpklss
  n/plot hpkz.val*${cct21}%$kk ! ! "same notit nsta connectpoints black $xyl line=solid,2"
  addtext 0.03 2000 '*** P(k)-LSS ***' 'font=helvetica,bolditalic,16 black'
endscript  

setaxesatt 'minorticks font=helvetica,bold,16 autofontsize'

Rad2Deg = 180/3.141596
R2D2 = $Rad2Deg*$Rad2Deg
Tsys = 50 
Delnu = 1.e6

#  A z = 0.7
z = 0.7
c = 3.e5
H = 100 
Da = 2500 
nu21 = 1.42e9
nu = $nu21/(1+$z)
pi23 = 8.*Pi*Pi*Pi

Lambda = 0.357 
Lam2 = $Lambda*$Lambda 

Dx = 100
Dy = 100

defscript calcul 
  Dxol = $Dx/$Lambda 
  Dyol = $Dy/$Lambda 

  FOV = (1.2*1.2*$Lam2/$Dx/$Dy) 
  FOVDEG = $FOV*$R2D2 
  NPointing = 10000/$FOVDEG
  tinteg = 365*86400/$NPointing
  PNOISE = $Tsys*$Tsys/$tinteg/$Dxol/$Dyol
  PNOISE = $PNOISE*$Da*$Da*$c/$H*(1+$z)/$nu
#   PNOISE = $PNOISE*$Da*$Da*$c/$H*(1+$z)/$nu/$pi23
  PNOISE = $PNOISE*1.e6
  echo " FOV = $FOV deg^2  NPointing= $NPointing"
  echo " tinteg= $tinteg sec  PNOISE= $PNOISE mK^2"
endscript 

defscript plnoisedish 
  Dx = 100
  Dy = 100
  calcul 
  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' 
  Dx = 200
  Dy = 200
  calcul 
  plot2d noiseD2 x/$Da val*$PNOISE nb>10 'same line=solid,2  cpts grey nsta notit' 
  set lines ( '100mDish' '200mDish'  )
  set cols ( navyblue grey  )
  textdrawer lines cols 'frame font=helvetica,bold,16 inset=0.1,0.3,0.7,0.8' 
  setaxelabels 'k (Mpc^-1) ' 'PNoise(k) mk^2 Mpc^3' 'font=helvetica,bolditalic,16' 
endscript 

defscript plnoiseA
  Dx = 5*0.95
  Dy = 5*0.95
  calcul 
  plot2d noisemdf64 x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts red nsta notit'
#  plot2d noisemds x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts red nsta notit'
  plot2d noisemdsB x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts siennared nsta notit'
  plot2d noisemdsC x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts violetred nsta notit'
  Dx = 25*0.3
  Dy = 0.5*0.9
  calcul 
  plot2d noise2cyl x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts forestgreen nsta notit' 
  plot2d noise2cylP x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts green nsta notit' 
  set lines (  'FilledA:64x5mD' 'SparseB:72x5mD' 'SparseC:129x5mD' 'Pitts2Cyl(=64C)' 'PerfPitts2Cyl(=64C)' )
  set cols ( red siennared magenta  forestgreen green )
  textdrawer lines cols 'frame font=helvetica,bold,16  inset=0.1,0.3,0.15,0.35' 
  settitle ' PNoise(k) : Dishes/Cylinders, 64/72/129 channels' ' ' 'font=helvetica,bold,16' 
endscript 

defscript plnoiseB
  Dx = 5*0.95
  Dy = 5*0.95
  calcul 
  plot2d noisemdsC x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts red nsta notit'
  plot2d noisemdf x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts gold nsta notit' 
  Dx = 10*0.95
  Dy = 0.5*0.9
  calcul 
  plot2d noise3cyl x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts violet nsta notit' 
  plot2d noise3cylP x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts magenta nsta notit' 

  Dx = 12*0.95
  Dy = 0.5*0.9
  calcul 
  plot2d noisefcyl x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts blue nsta notit' 
  plot2d noisefcylP x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts royalblue nsta notit' 

  set lines ( 'SparseC:129x5mD' 'Filled400x5m' '3xCyl10x64m(=384C)' 'Perf3xCyl10x64m(=384C)' )
  set lines ( $lines '8Cyl12x96m(=1536C)' 'Perf8Cyl12x96m(=1536C)' )
  set cols ( red gold  violet magenta  blue royalblue )
  textdrawer lines cols 'frame font=helvetica,bold,16 inset=0.25,0.55,0.6,0.95' 
  settitle ' PNoise(k) : Dishes/Cylinders, 129/384/400/1536 channels' ' ' 'font=helvetica,bold,16'	
endscript

defscript xxx
newwin
disp mfill 'h2disp=img colbr128 notit'
setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 
settitle ' u-v coverage , 400 x 5m Dishes - No Pointing' ' ' 'font=helvetica,bold,16' 

disp dish 'h2disp=img colbr128 notit'
setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 
settitle ' u-v coverage , 100 m Dish' ' ' 'font=helvetica,bold,16' 

disp msparsfp 'h2disp=img colbr128 notit'
setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 
settitle ' u-v coverage , 63 x 5m Dishes T-config - No Pointing' ' ' 'font=helvetica,bold,16' 

disp mspars 'h2disp=img colbr128 notit'
setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 
settitle ' u-v coverage , 63 x 5m Dishes T-config - Pointing ~Pi/4' ' ' 'font=helvetica,bold,16' 

disp mcylf 'h2disp=img colbr128 notit'
setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 
settitle ' u-v coverage , Filled Cylinder Array 8 Cyl 12mx96m ' ' ' 'font=helvetica,bold,16' 

disp mcylfP 'h2disp=img colbr128 notit'
setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 
settitle ' u-v coverage , Perfect Filled Cylinder Array 8 Cyl 12mx96m ' ' ' 'font=helvetica,bold,16' 

disp m2cyl 'h2disp=img colbr128 notit'
setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 
settitle ' u-v coverage , Pittsburgh 2 Cyl 16mx8m , dist=25m ' ' ' 'font=helvetica,bold,16' 

disp m2cylP 'h2disp=img colbr128 notit'
setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 
settitle ' u-v coverage , Perfect-Pitts. 2 Cyl 16mx8m , dist=25m ' ' ' 'font=helvetica,bold,16' 


disp noiseD 'logy nsta' 
disp noiseD2 'same grey nsta'
disp noisemdf 'same cyan nsta'
disp noisemdf64 'same brown nsta'
disp noisemds 'same red nsta'
disp noisemdsC 'same orange nsta'

disp noisemdsB 'same yellow nsta'
disp noisemdsfp 'same yellow nsta'
# disp noisemdsd7 'same gold nsta'
disp noisefcyl 'same blue nsta'
disp noisefcylP 'same skyblue nsta'
disp noise3cylP 'same magenta nsta'
disp noise3cyl 'same violet nsta'
disp noise2cyl 'same forestgreen nsta'
disp noise2cylP 'same green nsta'

endscript 

defscript AA
  plnoisedish
  plpklss
  plnoiseA
endscript 

defscript BB
  plnoisedish
  plpklss
  plnoiseB
endscript 

defscript POSCOV
  disp posspB red
  setaxelabels ' X (meters) ' ' Y (meters) ' 'font=helvetica,bolditalic,16' 
  settitle ' Config B dish positions - 72 dishes  ' ' ' 'font=helvetica,bold,16' 
#  w2ps
   
  disp posspC red
  setaxelabels ' X (meters) ' ' Y (meters) ' 'font=helvetica,bolditalic,16' 
  settitle ' Config C dish positions - 129 dishes  ' ' ' 'font=helvetica,bold,16' 
#  w2ps
  disp mfill64 'h2disp=img colbr128 notit'
  setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 
  settitle ' u-v coverage , Filled 8x8 - 64 x 5m Dishes' ' ' 'font=helvetica,bold,16' 
#  w2ps
  disp m3cyl 'h2disp=img colbr128 notit'
  setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 
  settitle ' u-v coverage , 3 Cylinders 10mx64m ' ' ' 'font=helvetica,bold,16' 
#  w2ps
  disp mfill 'h2disp=img colbr128 notit'
  setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 
  settitle ' u-v coverage , Filled 20x20 - 400 x 5m Dishes' ' ' 'font=helvetica,bold,16' 
#  w2ps
  disp msparsC 'h2disp=img colbr128 notit'
  setaxelabels 'kx (Radian^-1)  k=1000 -> ~21 arcmin ' ' ky (Radian^-1) ' 'font=helvetica,bolditalic,16' 
  settitle 'u-v coverage, Sparse-C: 129x5mD Over 80mx80m (Rot~Pi/4)' ' ' 'font=helvetica,bold,16' 
  
endscript