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

openppf pknoise.ppf

setaxesatt 'minorticks font=helvetica,bold,16 autofontsize'
disp noiseD 'logy nsta' 
disp noiseD2 'same grey nsta'
disp noisemdf 'same orange nsta'
disp noisemds 'same red nsta'
disp noisemdsfp 'same yellow nsta'
disp noisemdsd7 'same gold nsta'
disp noisefcyl 'same blue nsta'
disp noisefcylP 'same skyblue nsta'
disp noise2cyl 'same forestgreen nsta'
disp noise2cylP 'same green nsta'

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 

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 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' 
Dx = 5*0.95
Dy = 5*0.95
calcul 
plot2d noisemdf x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts orange nsta notit' 
plot2d noisemds x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts red nsta notit'
Dx = 7.55*0.95
Dy = 7.5*0.95
calcul 
plot2d noisemdsd7 x/$Da val*$PNOISE nb>10 'same line=solid,2 cpts siennared nsta notit' 
Dx = 12*0.9
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 skyblue nsta notit' 
Dx = 8*0.9
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 ( '100mDish' '200mDish' '400x5mDish' '63x5mDish' '63x7.5mDish'  '8Cyl12x96m' 'Perf8Cyl12x96m' )
set lines ( $lines 'Pitts2Cyl' 'PerfPitts2Cyl' ) 
set cols ( black grey  orange red siennared  blue skyblue forestgreen green )
textdrawer lines cols 'frame font=helvetica,bold,16' 
setaxelabels 'k (Mpc^-1) ' 'PNoise(k) mk^2 Mpc^3' 'font=helvetica,bolditalic,16' 


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'