

import sys

from casa_tools import myclean

'''
Cleans an MS with a single channel given a mask and a model
'''

vis_1 = sys.argv[-5]
vis_2 = sys.argv[-4]
model = sys.argv[-3]
mask = sys.argv[-2]
out_root = sys.argv[-1]

if model == "None":
    model = None
else:
    model = [model] * 2
if mask == "None":
    mask = None
else:
    mask = [mask] * 2

myclean(vis=[vis_1, vis_2], imagename=out_root+'.clean', field='M33*',
        restfreq='1420.40575177MHz',
        mode='channel', width=1, nchan=1, start=1,
        cell='1.5arcsec', multiscale=[0, 4, 8, 20, 40, 80, 160],
        threshold='2mJy/beam', imagermode='mosaic', gain=0.1,
        imsize=[4096, 4096], weighting='natural', robust=0.0, niter=200000,
        pbcor=True, minpb=0.2, interpolation='linear', usescratch=False,
        phasecenter='J2000 01h33m50.904 +30d39m35.79', veltype='radio',
        modelimage=model, mask=mask)



import sys

from casa_tools import myclean

'''
Cleans an MS with a single channel given a mask and a model
'''

vis_1 = sys.argv[-5]
vis_2 = sys.argv[-4]
model = sys.argv[-3]
mask = sys.argv[-2]
out_root = sys.argv[-1]

if model == "None":
    model = None
else:
    model = [model] * 2
if mask == "None":
    mask = None
else:
    mask = [mask] * 2

myclean(vis=[vis_1, vis_2], imagename=out_root+'.clean', field='M33*',
        restfreq='1420.40575177MHz',
        mode='channel', width=1, nchan=1, start=1,
        cell='1.5arcsec', multiscale=[0, 4, 8, 20, 40, 80],
        threshold='2mJy/beam', imagermode='mosaic', gain=0.1,
        imsize=[4096, 4096], weighting='natural', robust=0.0, niter=200000,
        pbcor=True, minpb=0.2, interpolation='linear', usescratch=False,
        phasecenter='J2000 01h33m50.904 +30d39m35.79', veltype='radio',
        modelimage=model, mask=mask)

--- 
+++ 
@@ -26,7 +26,7 @@
 myclean(vis=[vis_1, vis_2], imagename=out_root+'.clean', field='M33*',
         restfreq='1420.40575177MHz',
         mode='channel', width=1, nchan=1, start=1,
-        cell='1.5arcsec', multiscale=[0, 4, 8, 20, 40, 80],
+        cell='1.5arcsec', multiscale=[0, 4, 8, 20, 40, 80, 160],
         threshold='2mJy/beam', imagermode='mosaic', gain=0.1,
         imsize=[4096, 4096], weighting='natural', robust=0.0, niter=200000,
         pbcor=True, minpb=0.2, interpolation='linear', usescratch=False,
Add larger scale to match C config only map
