1 | set user "AST" |
---|
2 | set toppath "/sps/baoradio/AmasNancay/${user}" |
---|
3 | set source "Abell85" |
---|
4 | set srclower "abell85" |
---|
5 | set date "20110812" |
---|
6 | set saveplot "1" |
---|
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 |
---|
16 | objaoper gain_${date}_${srclower} row 0 g0 |
---|
17 | objaoper gain_${date}_${srclower} row 1 g1 |
---|
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)" |
---|
22 | if ( $saveplot == "1" ) then |
---|
23 | echo "saving plot" |
---|
24 | w2eps ${toppath}/${source}/eps/gain/gain_${date}_${srclower}.eps |
---|
25 | endif |
---|
26 | |
---|
27 | #################################################### |
---|
28 | # Control of the gains |
---|
29 | #################################################### |
---|
30 | newwin 1 1 |
---|
31 | graphicatt "xylimits=0,25,0,0.006" |
---|
32 | del gainmoni gainmoni0 gainmoni1 |
---|
33 | |
---|
34 | openppf ${toppath}/${source}/${date}${srclower}/gain_monitor_${date}_${srclower}.ppf |
---|
35 | #attention les "monitor" ont ete utilise a une epoque ou il n'y avait pas de |
---|
36 | #reordering... |
---|
37 | ls |
---|
38 | #plot2d monitor n val n>0 "black cpts nsta notit" |
---|
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" |
---|
42 | settitle "Gain monitor $source ${date} Ch 0 (blue) Ch 1 (red)" |
---|
43 | if ( $saveplot == "1" ) then |
---|
44 | # w2eps ${toppath}/${source}/eps/gain/gain_monitor_${date}_${srclower}.eps |
---|
45 | endif |
---|
46 | |
---|
47 | del g0 g1 gainmoni gainmoni0 gainmoni1 |
---|
48 | del gain_${date}_${srclower} |
---|
49 | ls |
---|