set user "AST"
set toppath "/sps/baoradio/AmasNancay/${user}"
set source "Abell1205"
set srclower "abell1205"
set date   "20110924"
set saveplot "1"

####################################################
# Gains
####################################################
newwin 1 1
graphicatt "xylimits=1250,1500,0,0.008"

openfits ${toppath}/${source}/${date}${srclower}/gain_${date}_${srclower}.fits
del g0 g1
objaoper gain_${date}_${srclower} row 0 g0
objaoper gain_${date}_${srclower} row 1 g1
plot2d g1 (n/8192)*250+1250 val n>0 "red cpts nsta notit"
plot2d g0 (n/8192)*250+1250 val n>0 "same blue cpts nsta notit"
settitle "Gain $source ${date} Ch 0 (blue) Ch 1 (red)"
setaxelabels "Freq. (MHz)" "I (a.u)"
if ( $saveplot == "1" ) then
  echo "saving plot"
#  w2eps ${toppath}/${source}/eps/gain/gain_${date}_${srclower}.eps
endif

####################################################
# Control of the gains
####################################################
newwin 1 1
graphicatt "xylimits=0,25,0,0.006"
del gainmoni gainmoni0 gainmoni1

openppf ${toppath}/${source}/${date}${srclower}/gain_monitor_${date}_${srclower}.ppf
#attention les "monitor" ont ete utilise a une epoque ou il n'y avait pas de
#reordering...
ls
#plot2d monitor n val n>0 "black cpts nsta notit"
plot2d gainmoni n val n>0 "black cpts nsta notit"
plot2d gainmoni0 n val n>0 "same blue cpts nsta notit"
plot2d gainmoni1 n val n>0 "same red cpts nsta notit"
settitle "Gain monitor $source ${date} Ch 0 (blue) Ch 1 (red)"
if ( $saveplot == "1" ) then
  w2eps ${toppath}/${source}/eps/gain/gain_monitor_${date}_${srclower}.eps
endif

del g0 g1 gainmoni gainmoni0 gainmoni1
del gain_${date}_${srclower}
ls