##################  Script plac8.pic #################
## Script pour tracer les 8 auto-correlations de PAON4 
if ( $# < 1 ) then 
  echo 'plac8.pic/missing argument '
  echo 'Usage: plac8 input_base_name   (without .ppf)' 
  return
endif
set flnm $1 
set titre 'PAON4 - transit'
#  Pour les axes 
setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
#  Ouvrir le fichier 
openppf $flnm.ppf 
set NT $flnm 
print $NT 

## Facteur de normalisation entre les voies x 1.e-4
## 1  1/3.48  1.06  0.92  1   0.98   0.377  1.32 
set FNORM ( 1.e-4 0.287e-4 1.06e-4 0.92e-4 1.e-4 0.96e-4 0.377e-4 1.32e-4  )   
##  couleurs pour les traces 
set COLS ( black red orange brown blue green skylbue forestgreen ) 

###### Trace de la puissance large bande    
##  On trace p1 en fct du temps 
newwin 1 1 800 600 
plot2d $NT time/3600 p1*$FNORM[0]  1 " nsta cpts notit $COLS[0] "
### trace de p2 ... p8
for i 1:8
   j = i+1 
   plot2d $NT time/3600 p$j*$FNORM[i]  1 "nsta cpts same notit $COLS[i] "
end

### Axes et titres 
setaxelabels 'Time (hours)' 'P1...P8 (arb.units)' 'font=helvetica,bolditalic,16 black'
settitle "$titre" ' ' 'font=helvetica,bold,18 black'

newwin 2 2 
for i 0:4
 j = i+1 
 plot2d $NT time/3600 p$j*$FNORM[i]  1 "nsta cpts notit $COLS[i] "
 settitle "Transit P$j" ' ' 'font=helvetica,bold,14 black'

end
newwin 2 2 
for i 4:8
 j = i+1 
 plot2d $NT time/3600 p$j*$FNORM[i]  1 "nsta cpts notit $COLS[i] "
 settitle "Transit P$j" ' ' 'font=helvetica,bold,14 black'
end

###### Trace de la puissance a +/-1 MHz @ 1420 MHz    
newwin 1 1 800 600 
plot2d $NT time/3600 (p1HI-p1HI5)*$FNORM[0]  1 " nsta cpts notit $COLS[0] "
### trace de p2 ... p8
for i 1:8
   j = i+1 
   plot2d $NT time/3600 (p${j}HI-p${j}HI5)*$FNORM[i]  1 "nsta cpts same notit $COLS[i] "
end

### Axes et titres 
setaxelabels 'Time (hours)' 'P1...P8 excess  @1420' 'font=helvetica,bolditalic,16 black'
settitle "$titre 1420+/-1 MHz - baseline (pjHI-pjHI5)" ' ' 'font=helvetica,bold,18 black'

###########################################################