set toppath "/sps/baoradio/AmasNancay/JEC"
set source "Abell85"
set srclower "abell85"
set date   "20110507"

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

####################################################
# Control of the gains
####################################################
newwin 1 1
graphicatt "xylimits=0,25,0,0.006"

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...
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"