| [595] | 1 | set user "AST"
|
|---|
| 2 | set toppath "/sps/baoradio/AmasNancay/${user}"
|
|---|
| [597] | 3 | set source "Abell1205"
|
|---|
| 4 | set srclower "abell1205"
|
|---|
| 5 | set date "20110924"
|
|---|
| [595] | 6 | set saveplot "1"
|
|---|
| [525] | 7 |
|
|---|
| 8 | ####################################################
|
|---|
| 9 | # Gains
|
|---|
| 10 | ####################################################
|
|---|
| 11 | newwin 1 1
|
|---|
| 12 | graphicatt "xylimits=1250,1500,0,0.008"
|
|---|
| 13 |
|
|---|
| 14 | openfits ${toppath}/${source}/${date}${srclower}/gain_${date}_${srclower}.fits
|
|---|
| 15 | del g0 g1
|
|---|
| [540] | 16 | objaoper gain_${date}_${srclower} row 0 g0
|
|---|
| 17 | objaoper gain_${date}_${srclower} row 1 g1
|
|---|
| [525] | 18 | plot2d g1 (n/8192)*250+1250 val n>0 "red cpts nsta notit"
|
|---|
| 19 | plot2d g0 (n/8192)*250+1250 val n>0 "same blue cpts nsta notit"
|
|---|
| 20 | settitle "Gain $source ${date} Ch 0 (blue) Ch 1 (red)"
|
|---|
| 21 | setaxelabels "Freq. (MHz)" "I (a.u)"
|
|---|
| [595] | 22 | if ( $saveplot == "1" ) then
|
|---|
| 23 | echo "saving plot"
|
|---|
| [597] | 24 | # w2eps ${toppath}/${source}/eps/gain/gain_${date}_${srclower}.eps
|
|---|
| [595] | 25 | endif
|
|---|
| [525] | 26 |
|
|---|
| 27 | ####################################################
|
|---|
| 28 | # Control of the gains
|
|---|
| 29 | ####################################################
|
|---|
| 30 | newwin 1 1
|
|---|
| 31 | graphicatt "xylimits=0,25,0,0.006"
|
|---|
| [595] | 32 | del gainmoni gainmoni0 gainmoni1
|
|---|
| [525] | 33 |
|
|---|
| 34 | openppf ${toppath}/${source}/${date}${srclower}/gain_monitor_${date}_${srclower}.ppf
|
|---|
| [540] | 35 | #attention les "monitor" ont ete utilise a une epoque ou il n'y avait pas de
|
|---|
| 36 | #reordering...
|
|---|
| [595] | 37 | ls
|
|---|
| 38 | #plot2d monitor n val n>0 "black cpts nsta notit"
|
|---|
| [525] | 39 | plot2d gainmoni n val n>0 "black cpts nsta notit"
|
|---|
| 40 | plot2d gainmoni0 n val n>0 "same blue cpts nsta notit"
|
|---|
| 41 | plot2d gainmoni1 n val n>0 "same red cpts nsta notit"
|
|---|
| [595] | 42 | settitle "Gain monitor $source ${date} Ch 0 (blue) Ch 1 (red)"
|
|---|
| 43 | if ( $saveplot == "1" ) then
|
|---|
| [597] | 44 | w2eps ${toppath}/${source}/eps/gain/gain_monitor_${date}_${srclower}.eps
|
|---|
| [595] | 45 | endif
|
|---|
| 46 |
|
|---|
| 47 | del g0 g1 gainmoni gainmoni0 gainmoni1
|
|---|
| 48 | del gain_${date}_${srclower}
|
|---|
| 49 | ls
|
|---|
