raise Error('The casa commands log is not executable!')

# # This file contains CASA commands run by the pipeline. Although all
# commands required to calibrate the data are included here, this file cannot
# be executed, nor does it contain heuristic and flagging calculations
# performed by pipeline code. This file is useful to understand which CASA
# commands are being run by each pipeline task. If one wishes to re-run the
# pipeline, one should use the pipeline script linked on the front page or By
# Task page of the weblog. Some stages may not have any commands listed here,
# e.g. hifa_importdata if conversion from ASDM to MS is not required.

# hifv_importdata(nocopy=True, vis=['VLASS1.1.sb34396276.eb34742211.58081.44143899306.ms'], session=['session_1'])
#
# If required, ASDMs are converted to measurement sets.
#

# hif_editimlist(parameter_file='parameter.5154.list')
#

# hif_transformimagedata(datacolumn='corrected', modify_weights=False, clear_pointing=True)
#
mstransform(timebin='0s',
            outputvis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T18t17.J120217+293000_P5154v1_2018_03_22T15_43_39.639/working/VLASS1.1.sb34396276.eb34742211.58081.44143899306_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T18t17.J120217+293000_P5154v1_2018_03_22T15_43_39.639/rawdata/VLASS1.1.sb34396276.eb34742211.58081.44143899306.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='3715,3714,3717,3716,3711,3713,3712,4721,3719,3718,4720,4348,4756,4723,5111,5119,4722,5148,5123,3128,3129,4149,4148,3122,3123,3120,3121,3126,3127,3124,3125,4741,3740,3761,3762,3763,3764,4338,4339,4336,4337,5144,5145,4346,5143,5140,5141,3135,3134,3137,3136,3131,3130,3133,3132,4748,3139,3138,3741,4138,4139,4730,4731,4736,4737,4130,4131,4132,4133,4134,4137,5129,5126,5136,4738,5112,5159,4358,4739,4142,4349,5161,2770,2771,2772,2773,5155,5154,5157,5156,3140,3141,3144,3145,3146,3147,3148,3149,4746,4345,4109,4129,4128,3743,4728,4343,4123,4122,4121,4120,4127,4126,4125,4124,5138,3746,4113,4341,3747,4340,3744,4740,2769,2768,4342,4713,3748,3749,2763,2762,2761,2760,2767,2766,2765,2764,3153,3152,3151,3150,3157,3156,3155,3154,4354,4355,3159,3158,4350,4351,4118,4119,4758,4116,4356,4117,4750,4711,4752,4753,4754,4114,4111,4743,5128,4145,4115,5132,4357,4144,4112,4147,3745,3759,3758,4759,4146,4110,4140,3751,3750,3753,3752,3755,3754,3757,3756,5124,5131,4143,3168,3169,3166,3167,3164,3165,3162,3163,3160,3161,4749,4344,5127,3364,3365,3366,3367,3360,3361,3362,3363,5130,5151,4352,5150,4353,5153,4714,5152,5120,4715,5110,4732,5133,5114,4716,5108,4717,3728,3729,3724,3725,3726,3727,3720,3721,3722,3723,4347,5121,4154,4733,4155,2756,2757,2754,2755,2752,2753,2750,2751,5158,4708,4726,4751,4757,2758,2759,5115,3760,5117,5116,4361,4360,5113,4747,4745,4712,5118,4755,4761,4760,5109,3742,3739,4141,3736,3735,3734,3733,3732,3731,3730,3348,3349,4162,4161,4160,3342,3343,5149,3346,3347,3344,3345,5146,5147,2749,2748,4709,5160,4744,5139,5137,5142,3355,3354,3357,3356,3351,3350,3353,3352,4729,3359,3358,5125,4152,4153,4150,4151,4156,4157,3119,3118,3117,3116,4158,4159,4718,4719,5122,4725,4742,4359,4724,4710,4727',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 12:02:17.202 +29.30.0.0000',
       scan=['65,68,71,74,87,78,81,52,55,84,58,90,62'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T18t17.J120217+293000_P5154v1_2018_03_22T15_43_39.639/working/VLASS1.1.sb34396276.eb34742211.58081.44143899306_split.ms'],
       imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='2748,2749,2750,2751,2752,2753,2754,2755,2756,2757,2758,2759,2760,2761,2762,2763,2764,2765,2766,2767,2768,2769,2770,2771,2772,2773,3116,3117,3118,3119,3120,3121,3122,3123,3124,3125,3126,3127,3128,3129,3130,3131,3132,3133,3134,3135,3136,3137,3138,3139,3140,3141,3144,3145,3146,3147,3148,3149,3150,3151,3152,3153,3154,3155,3156,3157,3158,3159,3160,3161,3162,3163,3164,3165,3166,3167,3168,3169,3342,3343,3344,3345,3346,3347,3348,3349,3350,3351,3352,3353,3354,3355,3356,3357,3358,3359,3360,3361,3362,3363,3364,3365,3366,3367,3711,3712,3713,3714,3715,3716,3717,3718,3719,3720,3721,3722,3723,3724,3725,3726,3727,3728,3729,3730,3731,3732,3733,3734,3735,3736,3739,3740,3741,3742,3743,3744,3745,3746,3747,3748,3749,3750,3751,3752,3753,3754,3755,3756,3757,3758,3759,3760,3761,3762,3763,3764,4109,4110,4111,4112,4113,4114,4115,4116,4117,4118,4119,4120,4121,4122,4123,4124,4125,4126,4127,4128,4129,4130,4131,4132,4133,4134,4137,4138,4139,4140,4141,4142,4143,4144,4145,4146,4147,4148,4149,4150,4151,4152,4153,4154,4155,4156,4157,4158,4159,4160,4161,4162,4336,4337,4338,4339,4340,4341,4342,4343,4344,4345,4346,4347,4348,4349,4350,4351,4352,4353,4354,4355,4356,4357,4358,4359,4360,4361,4708,4709,4710,4711,4712,4713,4714,4715,4716,4717,4718,4719,4720,4721,4722,4723,4724,4725,4726,4727,4728,4729,4730,4731,4732,4733,4736,4737,4738,4739,4740,4741,4742,4743,4744,4745,4746,4747,4748,4749,4750,4751,4752,4753,4754,4755,4756,4757,4758,4759,4760,4761,5108,5109,5110,5111,5112,5113,5114,5115,5116,5117,5118,5119,5120,5121,5122,5123,5124,5125,5126,5127,5128,5129,5130,5131,5132,5133,5136,5137,5138,5139,5140,5141,5142,5143,5144,5145,5146,5147,5148,5149,5150,5151,5152,5153,5154,5155,5156,5157,5158,5159,5160,5161',
       outframe='LSRK', gridder='mosaic', cycleniter=100, niter=0,
       datacolumn='data', savemodel='none', reffreq='3.0GHz', restoration=False,
       intent='OBSERVE_TARGET#UNSPECIFIED', cyclefactor=2.0, robust=1.0,
       parallel=False, stokes='I', nchan=-1, conjbeams=False,
       deconvolver='mtmfs', weighting='briggs', pbcor=False, pblimit=0.2,
       restoringbeam='common', specmode='mfs', chanchunks=-1, interactive=0)
tclean(phasecenter='J2000 12:02:17.202 +29.30.0.0000',
       scan=['65,68,71,74,87,78,81,52,55,84,58,90,62'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T18t17.J120217+293000_P5154v1_2018_03_22T15_43_39.639/working/VLASS1.1.sb34396276.eb34742211.58081.44143899306_split.ms'],
       imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00052679368935932864, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='2748,2749,2750,2751,2752,2753,2754,2755,2756,2757,2758,2759,2760,2761,2762,2763,2764,2765,2766,2767,2768,2769,2770,2771,2772,2773,3116,3117,3118,3119,3120,3121,3122,3123,3124,3125,3126,3127,3128,3129,3130,3131,3132,3133,3134,3135,3136,3137,3138,3139,3140,3141,3144,3145,3146,3147,3148,3149,3150,3151,3152,3153,3154,3155,3156,3157,3158,3159,3160,3161,3162,3163,3164,3165,3166,3167,3168,3169,3342,3343,3344,3345,3346,3347,3348,3349,3350,3351,3352,3353,3354,3355,3356,3357,3358,3359,3360,3361,3362,3363,3364,3365,3366,3367,3711,3712,3713,3714,3715,3716,3717,3718,3719,3720,3721,3722,3723,3724,3725,3726,3727,3728,3729,3730,3731,3732,3733,3734,3735,3736,3739,3740,3741,3742,3743,3744,3745,3746,3747,3748,3749,3750,3751,3752,3753,3754,3755,3756,3757,3758,3759,3760,3761,3762,3763,3764,4109,4110,4111,4112,4113,4114,4115,4116,4117,4118,4119,4120,4121,4122,4123,4124,4125,4126,4127,4128,4129,4130,4131,4132,4133,4134,4137,4138,4139,4140,4141,4142,4143,4144,4145,4146,4147,4148,4149,4150,4151,4152,4153,4154,4155,4156,4157,4158,4159,4160,4161,4162,4336,4337,4338,4339,4340,4341,4342,4343,4344,4345,4346,4347,4348,4349,4350,4351,4352,4353,4354,4355,4356,4357,4358,4359,4360,4361,4708,4709,4710,4711,4712,4713,4714,4715,4716,4717,4718,4719,4720,4721,4722,4723,4724,4725,4726,4727,4728,4729,4730,4731,4732,4733,4736,4737,4738,4739,4740,4741,4742,4743,4744,4745,4746,4747,4748,4749,4750,4751,4752,4753,4754,4755,4756,4757,4758,4759,4760,4761,5108,5109,5110,5111,5112,5113,5114,5115,5116,5117,5118,5119,5120,5121,5122,5123,5124,5125,5126,5127,5128,5129,5130,5131,5132,5133,5136,5137,5138,5139,5140,5141,5142,5143,5144,5145,5146,5147,5148,5149,5150,5151,5152,5153,5154,5155,5156,5157,5158,5159,5160,5161',
       outframe='LSRK', gridder='mosaic', cycleniter=100, niter=20000,
       datacolumn='data', savemodel='none', reffreq='3.0GHz', restoration=True,
       intent='OBSERVE_TARGET#UNSPECIFIED', cyclefactor=2.0, robust=1.0,
       calcpsf=False, restart=True, stokes='I', nchan=-1, conjbeams=False,
       deconvolver='mtmfs', weighting='briggs', pbcor=False, calcres=False,
       pblimit=0.2, restoringbeam='common', specmode='mfs', chanchunks=-1,
       interactive=0)

# hifv_pbcor(pipelinemode="automatic")
#
impbcor(cutoff=-1.0, stretch=False,
        outfile='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.pb.tt0')

# hif_makermsimages(pipelinemode="automatic")
#
imdev(zscore=-1, ylength='60arcsec', stretch=False, interp='cubic',
      outfile='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.image.pbcor.tt0.rms',
      imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.image.pbcor.tt0',
      grid=[10, 10], overwrite=True, maxiter=-1, xlength='60arcsec',
      anchor='ref', stattype='xmadm', statalg='chauvenet')

# hif_makecutoutimages(pipelinemode="automatic")
#
imstat(imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.image.pbcor.tt0')
imhead(imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.image.pbcor.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.image.pbcor.tt0.subim',
           imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.image.pbcor.tt0')
imstat(imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imhead(imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.image.pbcor.tt0.rms.subim',
           imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imstat(imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.psf.tt0')
imhead(imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.psf.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.psf.tt0.subim',
           imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.psf.tt0')
imstat(imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.residual.tt0')
imhead(imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.residual.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.residual.tt0.subim',
           imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.residual.tt0')
imstat(imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.pb.tt0')
imhead(imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.pb.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.pb.tt0.subim',
           imagename='VLASS1.1.ql.T18t17.J120217+293000.10.2048.v1.I.iter1.pb.tt0')

# hifv_exportvlassdata(pipelinemode="automatic")
#