listobs
2011-04-11 15:56:54 INFO listobs Fields: 3
2011-04-11 15:56:54 INFO listobs ID Code Name RA Decl Epoch SrcId nVis
2011-04-11 15:56:54 INFO listobs 0 none J0006-063 00:06:13.8920 -06.23.35.3350 J2000 0 582
2011-04-11 15:56:54 INFO listobs 1 none J0108+015 01:08:38.7710 +01.35.00.3170 J2000 1 4854
2011-04-11 15:56:54 INFO listobs 2 none NGC253 00:47:33.1300 -25.17.17.7997 J2000 2 23112
2011-04-11 15:56:54 INFO listobs (nVis = Total number of time/baseline visibilities per field)
2011-04-11 15:56:54 INFO listobs Spectral Windows: (3 unique spectral windows and 2 unique polarization setups)
2011-04-11 15:56:54 INFO listobs SpwID #Chans Frame Ch1(MHz) ChanWid(kHz)TotBW(kHz) Ref(MHz) Corrs
2011-04-11 15:56:54 INFO listobs 1 128 TOPO 229207.813 15625 2000000 229200 XX YY
2011-04-11 15:56:54 INFO listobs 2 1 TOPO 230184.375 1796875 1796875 229200 XX YY
2011-04-11 15:56:54 INFO listobs 0 4 TOPO 184550 1500000 7500000 183300 I
2011-04-11 15:56:54 INFO listobs Antennas: 6:
2011-04-11 15:56:54 INFO listobs ID Name Station Diam. Long. Lat.
2011-04-11 15:56:54 INFO listobs 0 DV02 J504 12.0 m -067.45.17.0 -22.53.23.0
2011-04-11 15:56:54 INFO listobs 1 DV03 N605 12.0 m -067.45.17.4 -22.53.23.9
2011-04-11 15:56:54 INFO listobs 2 DV04 T702 12.0 m -067.45.18.6 -22.53.24.1
2011-04-11 15:56:54 INFO listobs 3 DV05 J505 12.0 m -067.45.18.0 -22.53.22.8
2011-04-11 15:56:54 INFO listobs 4 DV07 J510 12.0 m -067.45.17.8 -22.53.23.5
2011-04-11 15:56:54 INFO listobs 5 PM02 T701 12.0 m -067.45.18.8 -22.53.22.2
Both calibrators are about the same brightness, use the one observed for more time (field '1' = 9.3 minutes, field '0' = 1.0 minutes, field '2'=NGC253=25 minutes):
---+ script
split('uid___A002_X1288b6_X273.ms',outputvis='X273spw1.ms',datacolumn='data',spw='1')
gaincal('X273spw1.ms',caltable='bpphase.gcal',field='1',spw='0:55~60',refant='DV02',calmode='p',solint='10s',minsnr=2.0)
bandpass('X273spw1.ms',caltable='bandpass.bcal',field='1',refant='DV02',solint='inf',solnorm=T,gaintable=['bpphase.gcal'])
plotcal(caltable='bpphase.gcal',xaxis='time',yaxis='phase',spw='0',iteration='antenna',subplot=221,plotrange=[0,0,-180,180])
plotcal(caltable='bp.gcal',xaxis='time',yaxis='phase',spw='0',iteration='antenna',subplot=221,plotrange=[0,0,-180,180])
gaincal('X273spw1.ms',caltable='intphase.gcal',field='',spw='0:10~120',combine='',refant='DV02',calmode='p',solint='int',minsnr=2.0,gaintable=['bandpass.bcal'])
gaincal('X273spw1.ms',caltable='scanphase.gcal',field='',spw='0:10~120',combine='',refant='DV02',calmode='p',solint='inf',minsnr=2.0,gaintable=['bandpass.bcal'])
gaincal('X273spw1.ms',caltable='amp.gcal',field='',spw='0:10~120',combine='',refant='DV02',calmode='ap',solint='inf',minsnr=2.0,gaintable=['bandpass.bcal','intphase.gcal'])
applycal(vis='X273spw1.ms',calwt=F,gaintable=['bandpass.bcal','scanphase.gcal','amp.gcal'],field='1',gainfield=['1','1','1'])
applycal(vis='X273spw1.ms',calwt=F,gaintable=['bandpass.bcal','scanphase.gcal','amp.gcal'],field='0',gainfield=['1','1','1'])
applycal(vis='X273spw1.ms',calwt=F,gaintable=['bandpass.bcal','scanphase.gcal','amp.gcal'],field='2',gainfield=['1','1','1'])
plotms(vis='X273spw1.ms',field='1',ydatacolumn='corrected',xaxis='frequency',yaxis='amp',correlation='XX',avgtime='1e8',avgscan=T,spw='0:10~120',antenna='DV02')
split(vis='X273spw1.ms',combine='scan',timebin='100000s',outputvis='X273spw1_timeavg.ms',antenna='*&&&')
--
ToddHunter - 2011-04-11