#!/usr/bin/env python
# This script describes the imaging for the ALMA LBC SV B6 data on Juno

# Requires CASA 4.2.2 or higher. Meant to be run interactively.

###################################################################
###################################################################
# This script assumes you have downloaded the ALMA Science Portal SV
# data product "CalibratedData". You can either apply any flagdata 
# commands and run the final applycals for each dataset (found in the 
# calibration subdirectory) or follow the script end to end.
###################################################################
###################################################################

# For each of 10 uv datasets, we create an initial model image, then
# run two iterations of phase-only selfcal, followed by an amplitude
# selfcal.  Finally, we create a final image with multiscale clean.

# First dataset: Create initial model
os.system('rm -rf Juno0.ms.cal.R0.*')
clean('Juno0.ms.cal', imagename='Juno0.ms.cal.R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy')
# 0.0331 x 0.0268"

# First dataset: phase self-cal first iteration
os.system('rm -rf Juno0.pcal1')
gaincal('Juno0.ms.cal', caltable='Juno0.pcal1',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='inf',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno0.pcal1',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno0.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno0.pcal1'],gainfield='',calwt=F,flagbackup=F)

os.system('rm -rf Juno0.pcal1_R0.*')
clean('Juno0.ms.cal', imagename='Juno0.pcal1_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno0.ms.cal.R0.mask')

# First dataset: phase self-cal second iteration
gaincal('Juno0.ms.cal', caltable='Juno0.pcal2',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='15.36s',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno0.pcal2',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno0.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno0.pcal2'],gainfield='',calwt=F,flagbackup=F)

os.system('rm -rf Juno0.pcal2_R0.*')
clean('Juno0.ms.cal', imagename='Juno0.pcal2_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno0.pcal1_R0.mask')

plotms('Juno0.ms.cal', xaxis='uvdist', yaxis='amp', coloraxis='spw',
       avgtime='1000', plotfile='Juno0.ms.cal.amp_vs_uvdist.png',
       ydatacolumn='corrected', overwrite=True)

# First dataset: amplitude self-cal 
gaincal('Juno0.ms.cal', caltable='Juno0.apcal',gaintype='T',
        refant='DV09',calmode='a',combine='spw,scan',
        gaintable=['Juno0.pcal2'],spwmap=[0,0,0,0],
        solint='300s',minsnr=3.0,minblperant=6)

#plotcal(caltable='Juno0.apcal',xaxis='time',yaxis='amp',timerange='',
#        spw='',iteration='antenna',subplot=421,plotrange=[0,0,0.6,1.4])

os.system('rm -f Juno0.apcal.png')
plotcal(caltable='Juno0.apcal',xaxis='time',yaxis='amp',timerange='',
        spw='',subplot=111,plotrange=[0,0,0.6,1.4],figfile='Juno0.apcal.png')

applycal('Juno0.ms.cal', spwmap=[[0,0,0,0],[0,0,0,0]],spw='',field='',
         gaintable=['Juno0.pcal2','Juno0.apcal'],gainfield='',calwt=F,
         flagbackup=F)

#os.system('rm -rf Juno0.apcal_R0.*')
#clean('Juno0.ms.cal', imagename='Juno0.apcal_R0',
#      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
#      interactive=True, weighting='briggs',robust=0.0,
#      niter=10000,threshold='0.2mJy',mask='Juno0.pcal1_R0.mask')

# First dataset: final imaging
os.system('rm -rf Juno0.apcal_R0ms.*')
clean('Juno0.ms.cal', imagename='Juno0.apcal_R0ms',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      multiscale=[0,5,15], negcomponent=10,
      niter=10000,threshold='0.15mJy',mask='Juno0.pcal2_R0.mask')

#####  Repeat everything for second dataset (and subsequent datasets)

os.system('rm -rf Juno1.ms.cal.R0.*')
clean('Juno1.ms.cal', imagename='Juno1.ms.cal.R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy')

os.system('rm -rf Juno1.pcal1')
gaincal('Juno1.ms.cal', caltable='Juno1.pcal1',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='inf',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno1.pcal1',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno1.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno1.pcal1'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno1.pcal1_R0.*')
clean('Juno1.ms.cal', imagename='Juno1.pcal1_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno1.ms.cal.R0.mask')

gaincal('Juno1.ms.cal', caltable='Juno1.pcal2',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='15.36s',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno1.pcal2',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno1.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno1.pcal2'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno1.pcal2_R0.*')
clean('Juno1.ms.cal', imagename='Juno1.pcal2_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno1.pcal1_R0.mask')

gaincal('Juno1.ms.cal', caltable='Juno1.apcal',gaintype='T',
        refant='DV09',calmode='a',combine='spw,scan',
        gaintable=['Juno1.pcal2'],spwmap=[0,0,0,0],
        solint='300s',minsnr=3.0,minblperant=6)

#plotcal(caltable='Juno1.apcal',xaxis='time',yaxis='amp',timerange='',
#        spw='',iteration='antenna',subplot=421,plotrange=[0,0,0.6,1.4])

os.system('rm -f Juno1.apcal.png')
plotcal(caltable='Juno1.apcal',xaxis='time',yaxis='amp',timerange='',
        spw='',subplot=111,plotrange=[0,0,0.6,1.4],figfile='Juno1.apcal.png')

applycal('Juno1.ms.cal', spwmap=[[0,0,0,0],[0,0,0,0]],spw='',field='',
         gaintable=['Juno1.pcal2','Juno1.apcal'],gainfield='',calwt=False,
         flagbackup=False)

os.system('rm -rf Juno1.apcal_R0ms.*')
clean('Juno1.ms.cal', imagename='Juno1.apcal_R0ms',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      multiscale=[0,5,15], negcomponent=10,
      niter=10000,threshold='0.15mJy',mask='Juno1.pcal2_R0.mask')

########  Third dataset
os.system('rm -rf Juno2.ms.cal.R0.*')
clean('Juno2.ms.cal', imagename='Juno2.ms.cal.R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy')

os.system('rm -rf Juno2.pcal1')
gaincal('Juno2.ms.cal', caltable='Juno2.pcal1',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='inf',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno2.pcal1',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno2.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno2.pcal1'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno2.pcal1_R0.*')
clean('Juno2.ms.cal', imagename='Juno2.pcal1_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno2.ms.cal.R0.mask')

gaincal('Juno2.ms.cal', caltable='Juno2.pcal2',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='15.36s',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno2.pcal2',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno2.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno2.pcal2'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno2.pcal2_R0.*')
clean('Juno2.ms.cal', imagename='Juno2.pcal2_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno2.pcal1_R0.mask')

gaincal('Juno2.ms.cal', caltable='Juno2.apcal',gaintype='T',
        refant='DV09',calmode='a',combine='spw,scan',
        gaintable=['Juno2.pcal2'],spwmap=[0,0,0,0],
        solint='300s',minsnr=3.0,minblperant=6)

#plotcal(caltable='Juno2.apcal',xaxis='time',yaxis='amp',timerange='',
#        spw='',iteration='antenna',subplot=421,plotrange=[0,0,0.6,1.4])

os.system('rm -f Juno2.apcal.png')
plotcal(caltable='Juno2.apcal',xaxis='time',yaxis='amp',timerange='',
        spw='',subplot=111,plotrange=[0,0,0.6,1.4],figfile='Juno2.apcal.png')

applycal('Juno2.ms.cal', spwmap=[[0,0,0,0],[0,0,0,0]],spw='',field='',
         gaintable=['Juno2.pcal2','Juno2.apcal'],gainfield='',calwt=False,
         flagbackup=False)

os.system('rm -rf Juno2.apcal_R0ms.*')
clean('Juno2.ms.cal', imagename='Juno2.apcal_R0ms',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      multiscale=[0,5,15], negcomponent=10,
      niter=10000,threshold='0.15mJy',mask='Juno2.pcal2_R0.mask')


#####  Fourth dataset

os.system('rm -rf Juno3.ms.cal.R0.*')
clean('Juno3.ms.cal', imagename='Juno3.ms.cal.R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy')

os.system('rm -rf Juno3.pcal1')
gaincal('Juno3.ms.cal', caltable='Juno3.pcal1',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='inf',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno3.pcal1',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno3.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno3.pcal1'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno3.pcal1_R0.*')
clean('Juno3.ms.cal', imagename='Juno3.pcal1_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno3.ms.cal.R0.mask')

gaincal('Juno3.ms.cal', caltable='Juno3.pcal2',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='15.36s',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno3.pcal2',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno3.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno3.pcal2'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno3.pcal2_R0.*')
clean('Juno3.ms.cal', imagename='Juno3.pcal2_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno3.pcal1_R0.mask')

gaincal('Juno3.ms.cal', caltable='Juno3.apcal',gaintype='T',
        refant='DV09',calmode='a',combine='spw,scan',
        gaintable=['Juno3.pcal2'],spwmap=[0,0,0,0],
        solint='300s',minsnr=3.0,minblperant=6)

#plotcal(caltable='Juno3.apcal',xaxis='time',yaxis='amp',timerange='',
#        spw='',iteration='antenna',subplot=421,plotrange=[0,0,0.6,1.4])

os.system('rm -f Juno3.apcal.png')
plotcal(caltable='Juno3.apcal',xaxis='time',yaxis='amp',timerange='',
        spw='',subplot=111,plotrange=[0,0,0.6,1.4],figfile='Juno3.apcal.png')

applycal('Juno3.ms.cal', spwmap=[[0,0,0,0],[0,0,0,0]],spw='',field='',
         gaintable=['Juno3.pcal2','Juno3.apcal'],gainfield='',calwt=False,
         flagbackup=False)

os.system('rm -rf Juno3.apcal_R0ms.*')
clean('Juno3.ms.cal', imagename='Juno3.apcal_R0ms',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      multiscale=[0,5,15], negcomponent=10,
      niter=10000,threshold='0.15mJy',mask='Juno3.pcal2_R0.mask')

##### Fifth dataset

os.system('rm -rf Juno4.ms.cal.R0.*')
clean('Juno4.ms.cal', imagename='Juno4.ms.cal.R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy')

os.system('rm -rf Juno4.pcal1')
gaincal('Juno4.ms.cal', caltable='Juno4.pcal1',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='inf',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno4.pcal1',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno4.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno4.pcal1'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno4.pcal1_R0.*')
clean('Juno4.ms.cal', imagename='Juno4.pcal1_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno4.ms.cal.R0.mask')

gaincal('Juno4.ms.cal', caltable='Juno4.pcal2',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='15.36s',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno4.pcal2',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno4.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno4.pcal2'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno4.pcal2_R0.*')
clean('Juno4.ms.cal', imagename='Juno4.pcal2_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno4.pcal1_R0.mask')

gaincal('Juno4.ms.cal', caltable='Juno4.apcal',gaintype='T',
        refant='DV09',calmode='a',combine='spw,scan',
        gaintable=['Juno4.pcal2'],spwmap=[0,0,0,0],
        solint='300s',minsnr=3.0,minblperant=6)

#plotcal(caltable='Juno4.apcal',xaxis='time',yaxis='amp',timerange='',
#        spw='',iteration='antenna',subplot=421,plotrange=[0,0,0.6,1.4])

os.system('rm -f Juno4.apcal.png')
plotcal(caltable='Juno4.apcal',xaxis='time',yaxis='amp',timerange='',
        spw='',subplot=111,plotrange=[0,0,0.6,1.4],figfile='Juno4.apcal.png')

applycal('Juno4.ms.cal', spwmap=[[0,0,0,0],[0,0,0,0]],spw='',field='',
         gaintable=['Juno4.pcal2','Juno4.apcal'],gainfield='',calwt=False,
         flagbackup=False)

os.system('rm -rf Juno4.apcal_R0ms.*')
clean('Juno4.ms.cal', imagename='Juno4.apcal_R0ms',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      multiscale=[0,5,15], negcomponent=10,
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.15mJy',mask='Juno4.pcal2_R0.mask')

##### Sixth dataset

os.system('rm -rf Juno5.ms.cal.R0.*')
clean('Juno5.ms.cal', imagename='Juno5.ms.cal.R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy')

os.system('rm -rf Juno5.pcal1')
gaincal('Juno5.ms.cal', caltable='Juno5.pcal1',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='inf',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno5.pcal1',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno5.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno5.pcal1'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno5.pcal1_R0.*')
clean('Juno5.ms.cal', imagename='Juno5.pcal1_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno5.ms.cal.R0.mask')

gaincal('Juno5.ms.cal', caltable='Juno5.pcal2',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='15.36s',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno5.pcal2',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno5.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno5.pcal2'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno5.pcal2_R0.*')
clean('Juno5.ms.cal', imagename='Juno5.pcal2_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno5.pcal1_R0.mask')

gaincal('Juno5.ms.cal', caltable='Juno5.apcal',gaintype='T',
        refant='DV09',calmode='a',combine='spw,scan',
        gaintable=['Juno5.pcal2'],spwmap=[0,0,0,0],
        solint='300s',minsnr=3.0,minblperant=6)

#plotcal(caltable='Juno5.apcal',xaxis='time',yaxis='amp',timerange='',
#        spw='',iteration='antenna',subplot=421,plotrange=[0,0,0.6,1.4])

os.system('rm -f Juno5.apcal.png')
plotcal(caltable='Juno5.apcal',xaxis='time',yaxis='amp',timerange='',
        spw='',subplot=111,plotrange=[0,0,0.6,1.4],figfile='Juno5.apcal.png')

applycal('Juno5.ms.cal', spwmap=[[0,0,0,0],[0,0,0,0]],spw='',field='',
         gaintable=['Juno5.pcal2','Juno5.apcal'],gainfield='',calwt=False,
         flagbackup=False)

os.system('rm -rf Juno5.apcal_R0ms.*')
clean('Juno5.ms.cal', imagename='Juno5.apcal_R0ms',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      multiscale=[0,5,15], negcomponent=10,
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.15mJy',mask='Juno5.pcal2_R0.mask')

##### Seventh dataset

os.system('rm -rf Juno6.ms.cal.R0.*')
clean('Juno6.ms.cal', imagename='Juno6.ms.cal.R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy')

os.system('rm -rf Juno6.pcal1')
gaincal('Juno6.ms.cal', caltable='Juno6.pcal1',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='inf',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno6.pcal1',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno6.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno6.pcal1'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno6.pcal1_R0.*')
clean('Juno6.ms.cal', imagename='Juno6.pcal1_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno6.ms.cal.R0.mask')

gaincal('Juno6.ms.cal', caltable='Juno6.pcal2',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='15.36s',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno6.pcal2',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno6.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno6.pcal2'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno6.pcal2_R0.*')
clean('Juno6.ms.cal', imagename='Juno6.pcal2_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno6.pcal1_R0.mask')

gaincal('Juno6.ms.cal', caltable='Juno6.apcal',gaintype='T',
        refant='DV09',calmode='a',combine='spw,scan',
        gaintable=['Juno6.pcal2'],spwmap=[0,0,0,0],
        solint='300s',minsnr=3.0,minblperant=6)

#plotcal(caltable='Juno6.apcal',xaxis='time',yaxis='amp',timerange='',
#        spw='',iteration='antenna',subplot=421,plotrange=[0,0,0.6,1.4])

os.system('rm -f Juno6.apcal.png')
plotcal(caltable='Juno6.apcal',xaxis='time',yaxis='amp',timerange='',
        spw='',subplot=111,plotrange=[0,0,0.6,1.4],figfile='Juno6.apcal.png')

applycal('Juno6.ms.cal', spwmap=[[0,0,0,0],[0,0,0,0]],spw='',field='',
         gaintable=['Juno6.pcal2','Juno6.apcal'],gainfield='',calwt=False,
         flagbackup=False)

os.system('rm -rf Juno6.apcal_R0ms.*')
clean('Juno6.ms.cal', imagename='Juno6.apcal_R0ms',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      multiscale=[0,5,15], negcomponent=10,
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.15mJy',mask='Juno6.pcal2_R0.mask')

##### eighth dataset

os.system('rm -rf Juno7.ms.cal.R0.*')
clean('Juno7.ms.cal', imagename='Juno7.ms.cal.R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy')

os.system('rm -rf Juno7.pcal1')
gaincal('Juno7.ms.cal', caltable='Juno7.pcal1',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='inf',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno7.pcal1',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno7.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno7.pcal1'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno7.pcal1_R0.*')
clean('Juno7.ms.cal', imagename='Juno7.pcal1_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno7.ms.cal.R0.mask')

gaincal('Juno7.ms.cal', caltable='Juno7.pcal2',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='15.36s',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno7.pcal2',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno7.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno7.pcal2'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno7.pcal2_R0.*')
clean('Juno7.ms.cal', imagename='Juno7.pcal2_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno7.pcal1_R0.mask')

gaincal('Juno7.ms.cal', caltable='Juno7.apcal',gaintype='T',
        refant='DV09',calmode='a',combine='spw,scan',
        gaintable=['Juno7.pcal2'],spwmap=[0,0,0,0],
        solint='300s',minsnr=3.0,minblperant=6)

#plotcal(caltable='Juno7.apcal',xaxis='time',yaxis='amp',timerange='',
#        spw='',iteration='antenna',subplot=421,plotrange=[0,0,0.6,1.4])

os.system('rm -f Juno7.apcal.png')
plotcal(caltable='Juno7.apcal',xaxis='time',yaxis='amp',timerange='',
        spw='',subplot=111,plotrange=[0,0,0.6,1.4],figfile='Juno7.apcal.png')

applycal('Juno7.ms.cal', spwmap=[[0,0,0,0],[0,0,0,0]],spw='',field='',
         gaintable=['Juno7.pcal2','Juno7.apcal'],gainfield='',calwt=False,
         flagbackup=False)

os.system('rm -rf Juno7.apcal_R0ms.*')
clean('Juno7.ms.cal', imagename='Juno7.apcal_R0ms',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      multiscale=[0,5,15], negcomponent=10,
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.15mJy',mask='Juno7.pcal2_R0.mask')

##### Ninth dataset

os.system('rm -rf Juno8.ms.cal.R0.*')
clean('Juno8.ms.cal', imagename='Juno8.ms.cal.R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy')

os.system('rm -rf Juno8.pcal1')
gaincal('Juno8.ms.cal', caltable='Juno8.pcal1',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='inf',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno8.pcal1',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno8.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno8.pcal1'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno8.pcal1_R0.*')
clean('Juno8.ms.cal', imagename='Juno8.pcal1_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno8.ms.cal.R0.mask')

gaincal('Juno8.ms.cal', caltable='Juno8.pcal2',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='15.36s',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno8.pcal2',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno8.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno8.pcal2'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno8.pcal2_R0.*')
clean('Juno8.ms.cal', imagename='Juno8.pcal2_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno8.pcal1_R0.mask')

gaincal('Juno8.ms.cal', caltable='Juno8.apcal',gaintype='T',
        refant='DV09',calmode='a',combine='spw,scan',
        gaintable=['Juno8.pcal2'],spwmap=[0,0,0,0],
        solint='300s',minsnr=3.0,minblperant=6)

#plotcal(caltable='Juno8.apcal',xaxis='time',yaxis='amp',timerange='',
#        spw='',iteration='antenna',subplot=421,plotrange=[0,0,0.6,1.4])

os.system('rm -f Juno8.apcal.png')
plotcal(caltable='Juno8.apcal',xaxis='time',yaxis='amp',timerange='',
        spw='',subplot=111,plotrange=[0,0,0.6,1.4],figfile='Juno8.apcal.png')

applycal('Juno8.ms.cal', spwmap=[[0,0,0,0],[0,0,0,0]],spw='',field='',
         gaintable=['Juno8.pcal2','Juno8.apcal'],gainfield='',calwt=False,
         flagbackup=False)

os.system('rm -rf Juno8.apcal_R0ms.*')
clean('Juno8.ms.cal', imagename='Juno8.apcal_R0ms',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      multiscale=[0,5,15], negcomponent=10,
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.15mJy',mask='Juno8.pcal2_R0.mask')

##### Tenth dataset

os.system('rm -rf Juno9.ms.cal.R0.*')
clean('Juno9.ms.cal', imagename='Juno9.ms.cal.R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy')

os.system('rm -rf Juno9.pcal1')
gaincal('Juno9.ms.cal', caltable='Juno9.pcal1',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='inf',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno9.pcal1',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno9.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno9.pcal1'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno9.pcal1_R0.*')
clean('Juno9.ms.cal', imagename='Juno9.pcal1_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno9.ms.cal.R0.mask')

gaincal('Juno9.ms.cal', caltable='Juno9.pcal2',gaintype='T',
        refant='DV09',calmode='p',combine='spw',
        solint='15.36s',minsnr=3.0,minblperant=6)

plotcal(caltable='Juno9.pcal2',xaxis='time',yaxis='phase',timerange='',
        spw='',iteration='antenna',subplot=421,plotrange=[0,0,-180,180])

applycal('Juno9.ms.cal', spwmap=[0,0,0,0],spw='',field='',
         gaintable=['Juno9.pcal2'],gainfield='',calwt=False,flagbackup=False)

os.system('rm -rf Juno9.pcal2_R0.*')
clean('Juno9.ms.cal', imagename='Juno9.pcal2_R0',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      niter=10000,threshold='0.2mJy',mask='Juno9.pcal1_R0.mask')

gaincal('Juno9.ms.cal', caltable='Juno9.apcal',gaintype='T',
        refant='DV09',calmode='a',combine='spw,scan',
        gaintable=['Juno9.pcal2'],spwmap=[0,0,0,0],
        solint='300s',minsnr=3.0,minblperant=6)

#plotcal(caltable='Juno9.apcal',xaxis='time',yaxis='amp',timerange='',
#        spw='',iteration='antenna',subplot=421,plotrange=[0,0,0.6,1.4])
os.system('rm -f Juno9.apcal.png')
plotcal(caltable='Juno9.apcal',xaxis='time',yaxis='amp',timerange='',
        spw='',subplot=111,plotrange=[0,0,0.6,1.4], figfile='Juno9.apcal.png')

applycal('Juno9.ms.cal', spwmap=[[0,0,0,0],[0,0,0,0]],spw='',field='',
         gaintable=['Juno9.pcal2','Juno9.apcal'],gainfield='',calwt=False,
         flagbackup=False)

os.system('rm -rf Juno9.apcal_R0ms.*')
clean('Juno9.ms.cal', imagename='Juno9.apcal_R0ms',
      mode='mfs',imagermode='csclean', imsize=1000,cell='0.005arcsec',
      interactive=True, weighting='briggs',robust=0.0,
      multiscale=[0,5,15], negcomponent=10,
      niter=10000,threshold='0.15mJy',mask='Juno9.pcal2_R0.mask')

# Output all images as FITS files
for i in range(10):
    exportfits('Juno%d.apcal_R0ms.image'%i, 'Juno%d.apcal_R0ms.image.fits'%i)