# ALMA Data Reduction Script

# Imaging
tgts=['UDS.0351','UDS.0697','UDS.0021','UDS.0570','UDS.0571','UDS.0572','UDS.0573','UDS.0332','UDS.0577','UDS.0186','UDS.0232','UDS.0087','UDS.0546','UDS.0702','UDS.0662','UDS.0075','UDS.0602','UDS.0601','UDS.0629','UDS.0079','UDS.0052','UDS.0036','UDS.0568','UDS.0545','UDS.0016','UDS.0111','UDS.0684','UDS.0296','UDS.0682','UDS.0030','UDS.0197','UDS.0249','UDS.0619','UDS.0679','UDS.0497','UDS.0599','UDS.0673','UDS.0712','UDS.0670','UDS.0390','UDS.0395','UDS.0227','UDS.0559','UDS.0517','UDS.0222']

# These are the parameters you are most likely to want to change. ? means fill in for yourself.

thevis = ['?']          # Calibrated MS to be imaged
active=True            # whether to clean interactively or not)
c = '?arcsec'       # pixel size  
isz = ?              # image size in pixels
r = 0.5                # robust parameter value  
sig = '250.e-6Jy'      # expected stopping threshold

#for t in range(len(tgts)):
for t in range(1):               # just the first field as a test starter
  os.system('rm -rf '+tgts[t]+'_sci.spw0_1_2_3.mfs.I.manual*')
  tclean(vis = thevis,
         imagename = tgts[t]+'_sci.spw0_1_2_3.mfs.I.manual',
         field = tgts[t],
         stokes = 'I',
         spw = '0,1,2,3',
         specmode = 'mfs',
         imsize = [isz, isz],
         cell = c,
         deconvolver = 'hogbom',
         niter = 100,                   # increase for bright sources
         cycleniter=10,                 # increase for bright sources
         weighting = 'briggs',
         robust = r,
         mask = '',
         gridder = 'standard',
         pbcor = True,
         threshold = sig,
         interactive = True
         )