1 #! /bin/csh -f 2 # 3 # 3/7/07 4 # 5 # Script to reduce CARMA OPHB11 data 6 # -- using 2 calibrators: 7 # offsetcal to calibrate offset in phase b/w wide and narrow band 8 # cal to calibrate phases 9 # 10 # pre-work: 11 # set vis=c0014.OphB11_C.1.miriad 12 # carmadata $vis.tar.gz 13 # puthd in=$vis/restfreq value=115.271204 type=double 14 # uvcat vis=$vis out=1751+096W.raw 'select=source(1751+096),-auto,time(11:29:00,11:38:20)' 15 # uvcat vis=$vis out=1751+096N.raw 'select=source(1751+096),-auto,time(11:08:00,11:30:00)' 16 # uvcat vis=$vis out=1625-254.raw 'select=source(1625-254),-auto' 17 # uvcat vis=$vis out=ophb11.raw 'select=source(ophb11),-auto' 18 19 umask 002 20 21 set task=flag # task to perform 22 set vis==c0014.OphB11_C.1.miriad # base project vis file 23 set offsetcal=1751+096 # calibrator for phase offset b/w Wide and Narr. bands 24 set cal=1625-254 # primary phase cal name 25 set source=ophb11 # source name 26 set refant=8 # refant for selfcal 27 set interval=25 # interval for selfcal 28 set scselect= # select for selfcal on phase cal 29 set nxy=3,3 # number of windows for plotting 30 set flux=1.6 31 32 set ointerval=9999 # interval for offset cal selfcal 33 set olinecal=chan,1,46,45,45 34 set linenarrow=chan,1,109,63,63 35 set axis=time,phase 36 37 set line= 38 39 set dev=1/xs # plotting device 40 set select= # select for uvplt 41 set options= # options for uvplt 42 43 set cutoff= 44 45 foreach _arg ($*) 46 set $_arg 47 end 48 49 goto $task 50 51 logs: 52 listobs vis= 53 54 55 renewoffset: 56 57 # Secondary Calibrator split into files with with differen correlator configs: 58 rm -rf "$offsetcal"N.red 59 cp -r "$offsetcal"N.raw "$offsetcal"N.red 60 rm -rf "$offsetcal"W.red 61 cp -r "$offsetcal"W.raw "$offsetcal"W.red 62 63 flagoff: 64 uvplt vis="$offsetcal"W.red device=$dev \ 65 nxy=$nxy select=$select options=$options axis=$axis 66 67 uvplt vis="$offsetcal"N.red device=$dev \ 68 nxy=$nxy select=$select options=$options axis=$axis 69 70 exit 0 71 72 selfoff: 73 74 # if applying BP cal to 1751, also apply to sec phase, and to source, depending on 75 # window ... 76 # 77 # or just include USB (calibrating Wide band) : 78 79 selfcal vis="$offsetcal"W.red refant=$refant \ 80 line=$olinecal options=noscale interval=$ointerval 81 82 gpplt vis="$offsetcal"W.red yaxis=phase yrange=-180,180 \ 83 device=$dev nxy=$nxy options=w 84 85 # Apply Wide Band calibration to Narrow Band: 86 gpcopy vis="$offsetcal"W.red out="$offsetcal"N.red mode=copy options=nopass 87 88 rm -rf "$offsetcal"Nref.red 89 uvcat vis="$offsetcal"N.red out="$offsetcal"Nref.red 90 91 # Another selfcal to get offset between Wide and Narrow bands: 92 selfcal vis="$offsetcal"Nref.red refant=$refant select=$select \ 93 line=$linenarrow options=phase interval=$ointerval 94 95 gpplt vis="$offsetcal"Nref.red yaxis=phase yrange=-180,180 \ 96 device=$dev nxy=$nxy options=w 97 98 exit 0 99 renew: 100 # CAREFUL: This will erase $cal.red to start from scratch 101 # uvcat vis=c0014.OphB11_C.1.miriad select=source(1625-254),-auto out=1625-254.raw 102 103 rm -rf $cal.red 104 cp -r $cal.raw $cal.red 105 106 flag: 107 uvplt vis=$cal.red device=$dev line= \ 108 nxy=$nxy select=$select options=$options 109 110 #---Flag data here--- 111 exit 0 112 113 selfcal: 114 115 #selfcal phase calibrator 116 117 selfcal vis="$cal".red flux=$flux refant=$refant select=$scselect options=amp,apriori,noscale interval=$interval line=$olinecal 118 119 gpplt: 120 gpplt vis="$cal".red yaxis=phase yrange=-180,180 \ 121 device=$dev nxy=$nxy options=w 122 123 sleep 2 124 125 rm -rf $cal.dm $cal.bm 126 invert vis="$cal".red map=$cal.dm beam=$cal.bm \ 127 imsize=237,237 cell=1,1 line=$olinecal options=systemp select=$scselect 128 129 rm -rf $cal.cln 130 131 clean map=$cal.dm beam=$cal.bm niters=10000 out=$cal.cln cutoff=$cutoff 132 133 rm -rf $cal.restor 134 restor model=$cal.cln map=$cal.dm beam=$cal.bm out=$cal.restor 135 136 cgdisp in=$cal.restor device=$dev 137 138 chmod -R 775 $cal.* 139 140 exit 0 141 142 renewsource: 143 # CAREFUL: This will erase $source.red to start from scratch 144 # uvcat vis=c0014.OphB11_C.1.miriad select=source(OPHB11),-auto out=ophb11.raw 145 rm -rf $source.red 146 cp -r $source.raw $source.red 147 148 flagsource: 149 uvplt vis=$source.red device=$dev line=$line \ 150 nxy=$nxy select=$select options=$options 151 152 #---Flag source data here--- 153 exit 0 154 155 156 invert: 157 158 # Copy and apply gains table from phase calibrator 159 gpcopy vis="$cal".red/ out=$source.red mode=copy 160 161 rm -rf "$source"GT.red 162 uvcat vis=$source.red out="$source"GT.red 163 164 # Copy and apply gains table from offset calibrator (narrow band only) 165 gpcopy vis="$offsetcal"Nref.red out="$source"GT.red mode=copy 166 167 rm -rf "$source"cal.red 168 uvcat vis="$source"GT.red out="$source"cal.red select='window(5)' 169 170 rm -rf $source.dm $source.bm 171 invert vis="$source"cal.red map=$source.dm beam=$source.bm imsize=237,237 \ 172 cell=1,1 options=systemp,double,mosaic line=$line sup=0 select=$select 173 174 #rm -rf $source.cln $source.restor 175 #clean map=$source.dm beam=$source.bm out=$source.cln \ 176 # niters=10000 cutoff=$cutoff 177 #minmax in=$source.cln 178 #restor map=$source.dm beam=$source.bm \ 179 # model=$source.cln out=$source.restor 180 181 chmod -R 775 $source.*