listobs(vis='A-W49A.ms')
plotms(vis='A-W49A.ms',plotfile='A-uvcov.png',xaxis='Uwave',yaxis='Vwave',spw='0:32',plotrange=[-1000000,1000000,-1000000,1000000])
plotms(vis='B-W49A.ms',plotfile='B-uvcov.png',xaxis='Uwave',yaxis='Vwave',spw='0:32',plotrange=[-1000000,1000000,-1000000,1000000])
plotms(vis='C-W49A.ms',plotfile='C-uvcov.png',xaxis='Uwave',yaxis='Vwave',spw='0:32',plotrange=[-1000000,1000000,-1000000,1000000])
plotms(vis='D-W49A.ms',plotfile='D-uvcov.png',xaxis='Uwave',yaxis='Vwave',spw='0:32',plotrange=[-1000000,1000000,-1000000,1000000])
tclean(vis='A-W49A.ms',imagename='Aonly',cell='0.05arcsec',imsize=[1280,1280],weighting='briggs',
       robust=0,specmode='mfs',phasecenter='J2000 19:10:14  +09.06.13.7') 
tclean(vis='B-W49A.ms',imagename='Bonly',cell='0.05arcsec',imsize=[1280,1280],weighting='briggs',
       robust=0,specmode='mfs',phasecenter='J2000 19:10:14  +09.06.13.7') 
tclean(vis='C-W49A.ms',imagename='Conly',cell='0.05arcsec',imsize=[1280,1280],weighting='briggs',
       robust=0,specmode='mfs',phasecenter='J2000 19:10:14  +09.06.13.7') 
tclean(vis='D-W49A.ms',imagename='Donly',cell='0.05arcsec',imsize=[1280,1280],weighting='briggs',
       robust=0,specmode='mfs',phasecenter='J2000 19:10:14  +09.06.13.7')
plotms(vis='A-W49A.ms',plotfile='A-vis-weight.png',spw='0:32',xaxis='uvwave',yaxis='wt',coloraxis='antenna1')
plotms(vis='B-W49A.ms',plotfile='B-vis-weight.png',spw='0:32',xaxis='uvwave',yaxis='wt',coloraxis='antenna1')
plotms(vis='C-W49A.ms',plotfile='C-vis-weight.png',spw='0:32',xaxis='uvwave',yaxis='wt',coloraxis='antenna1')
plotms(vis='D-W49A.ms',plotfile='D-vis-weight.png',spw='0:32',xaxis='uvwave',yaxis='wt',coloraxis='antenna1')
statwt(vis='A-W49A.ms',datacolumn='data')
statwt(vis='B-W49A.ms',datacolumn='data')
statwt(vis='C-W49A.ms',datacolumn='data')
statwt(vis='D-W49A.ms',datacolumn='data')
plotms(vis='A-W49A.ms',plotfile='adjusted-A-vis-weight.png',spw='0:32',xaxis='uvwave',yaxis='wt',coloraxis='antenna1')
plotms(vis='B-W49A.ms',plotfile='adjusted-B-vis-weight.png',spw='0:32',xaxis='uvwave',yaxis='wt',coloraxis='antenna1')
plotms(vis='C-W49A.ms',plotfile='adjusted-C-vis-weight.png',spw='0:32',xaxis='uvwave',yaxis='wt',coloraxis='antenna1')
plotms(vis='D-W49A.ms',plotfile='adjusted-D-vis-weight.png',spw='0:32',xaxis='uvwave',yaxis='wt',coloraxis='antenna1')
concat(vis=['A-W49A.ms','B-W49A.ms','C-W49A.ms','D-W49A.ms'],concatvis='ABCD-W49A.ms')
plotms(vis='ABCD-W49A.ms',plotfile='ABCD-uvcov.png',xaxis='Uwave',yaxis='Vwave',spw='*:32')
tclean(vis=['ABCD-W49A.ms'], imagename='W49A_combinedABCD',specmode='mfs',cell='0.05arcsec',imsize=[1280,1280],
niter=10000, weighting='briggs',robust=0, phasecenter='J2000 19:10:14  +09.06.13.7',threshold='2mJy')