| [626] | 1 | set user "JEC" |
|---|
| [595] | 2 | set toppath "/sps/baoradio/AmasNancay/${user}" |
|---|
| [626] | 3 | set source "NGC4383" |
|---|
| 4 | set srclower "ngc4383" |
|---|
| 5 | set date "20111119" |
|---|
| [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 |
|---|
| [626] | 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 |
|---|
