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.sb34517306.eb34782681.58097.74992206019.ms'], session=['session_1'])
#
# If required, ASDMs are converted to measurement sets.
#

# hif_editimlist(parameter_file='parameter.12529.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_T19t18.J142142+343000_P12529v1_2018_02_01T00_19_21.949/working/VLASS1.1.sb34517306.eb34782681.58097.74992206019_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T19t18.J142142+343000_P12529v1_2018_02_01T00_19_21.949/rawdata/VLASS1.1.sb34517306.eb34782681.58097.74992206019.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='3710,5316,4479,4478,5307,3659,4471,4470,4473,4472,4475,4474,4477,4476,3258,3259,4906,5277,4466,5271,5272,4509,4901,4900,4903,4902,4905,3255,3256,3257,4868,4869,5305,4866,4867,4888,4889,4467,5266,4880,4881,4882,4883,4884,4885,4886,4887,4515,4898,4875,4874,4877,4876,4871,4870,4873,4872,4503,5301,5302,5303,4879,4878,5269,5268,4484,4075,4486,4487,4480,4481,4482,4074,4105,5275,4488,4077,4104,5289,4899,4076,4109,4502,5288,4893,4071,4890,4897,4896,4895,3686,3276,3277,3274,3275,3272,3273,3270,3271,3306,3673,3670,3305,3676,3677,3278,3279,4485,5313,5312,5311,5310,5317,5306,5315,5314,4497,4496,4495,4494,4493,4492,4063,4108,5283,4499,4498,3665,3664,3667,3666,3661,3660,3663,3662,4483,3669,3668,3261,3260,3263,3262,3265,3264,3267,3266,3269,3268,5308,5309,5298,4908,5273,4114,4489,4506,4112,4068,4113,4069,5299,4110,5278,4111,5270,5282,4904,4101,4100,4103,4102,3298,3299,4107,4106,3294,3295,3296,3297,3290,3291,3292,3293,3698,3699,4465,4066,4067,4064,4065,3690,3691,3692,3693,3694,3695,3696,3697,3704,5297,4508,4894,5304,4907,3289,3288,3287,3286,3285,3284,3283,3282,3689,3688,4079,4078,3683,3682,3681,3680,3687,3678,4073,4072,3679,5279,4080,4081,4082,4083,4084,4085,4086,4087,5296,5267,5294,5295,5293,5290,3672,4070,3304,3671,5276,4909,4512,3302,4510,4511,4516,4514,3303,3674,3675,4093,4092,4091,4090,4097,4096,4095,4094,5281,5280,4099,4098,5285,5284,5287,5286,4505,4504,3708,3709,4501,4500,4468,4469,3702,3703,3700,3701,3706,3707,3300,3705,5274,4916,4917,4914,4915,4912,4913,4910,4911,4507,4513,5300,3301',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 14:21:42.715 +34.30.0.0000',
       scan=['65,68,71,75,78,81,84,87,56,90,59,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_T19t18.J142142+343000_P12529v1_2018_02_01T00_19_21.949/working/VLASS1.1.sb34517306.eb34782681.58097.74992206019_split.ms'],
       imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='3255,3256,3257,3258,3259,3260,3261,3262,3263,3264,3265,3266,3267,3268,3269,3270,3271,3272,3273,3274,3275,3276,3277,3278,3279,3282,3283,3284,3285,3286,3287,3288,3289,3290,3291,3292,3293,3294,3295,3296,3297,3298,3299,3300,3301,3302,3303,3304,3305,3306,3659,3660,3661,3662,3663,3664,3665,3666,3667,3668,3669,3670,3671,3672,3673,3674,3675,3676,3677,3678,3679,3680,3681,3682,3683,3686,3687,3688,3689,3690,3691,3692,3693,3694,3695,3696,3697,3698,3699,3700,3701,3702,3703,3704,3705,3706,3707,3708,3709,3710,4063,4064,4065,4066,4067,4068,4069,4070,4071,4072,4073,4074,4075,4076,4077,4078,4079,4080,4081,4082,4083,4084,4085,4086,4087,4090,4091,4092,4093,4094,4095,4096,4097,4098,4099,4100,4101,4102,4103,4104,4105,4106,4107,4108,4109,4110,4111,4112,4113,4114,4465,4466,4467,4468,4469,4470,4471,4472,4473,4474,4475,4476,4477,4478,4479,4480,4481,4482,4483,4484,4485,4486,4487,4488,4489,4492,4493,4494,4495,4496,4497,4498,4499,4500,4501,4502,4503,4504,4505,4506,4507,4508,4509,4510,4511,4512,4513,4514,4515,4516,4866,4867,4868,4869,4870,4871,4872,4873,4874,4875,4876,4877,4878,4879,4880,4881,4882,4883,4884,4885,4886,4887,4888,4889,4890,4893,4894,4895,4896,4897,4898,4899,4900,4901,4902,4903,4904,4905,4906,4907,4908,4909,4910,4911,4912,4913,4914,4915,4916,4917,5266,5267,5268,5269,5270,5271,5272,5273,5274,5275,5276,5277,5278,5279,5280,5281,5282,5283,5284,5285,5286,5287,5288,5289,5290,5293,5294,5295,5296,5297,5298,5299,5300,5301,5302,5303,5304,5305,5306,5307,5308,5309,5310,5311,5312,5313,5314,5315,5316,5317',
       outframe='LSRK', gridder='mosaic', cycleniter=500, 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 14:21:42.715 +34.30.0.0000',
       scan=['65,68,71,75,78,81,84,87,56,90,59,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_T19t18.J142142+343000_P12529v1_2018_02_01T00_19_21.949/working/VLASS1.1.sb34517306.eb34782681.58097.74992206019_split.ms'],
       imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.0004798826961236045, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='3255,3256,3257,3258,3259,3260,3261,3262,3263,3264,3265,3266,3267,3268,3269,3270,3271,3272,3273,3274,3275,3276,3277,3278,3279,3282,3283,3284,3285,3286,3287,3288,3289,3290,3291,3292,3293,3294,3295,3296,3297,3298,3299,3300,3301,3302,3303,3304,3305,3306,3659,3660,3661,3662,3663,3664,3665,3666,3667,3668,3669,3670,3671,3672,3673,3674,3675,3676,3677,3678,3679,3680,3681,3682,3683,3686,3687,3688,3689,3690,3691,3692,3693,3694,3695,3696,3697,3698,3699,3700,3701,3702,3703,3704,3705,3706,3707,3708,3709,3710,4063,4064,4065,4066,4067,4068,4069,4070,4071,4072,4073,4074,4075,4076,4077,4078,4079,4080,4081,4082,4083,4084,4085,4086,4087,4090,4091,4092,4093,4094,4095,4096,4097,4098,4099,4100,4101,4102,4103,4104,4105,4106,4107,4108,4109,4110,4111,4112,4113,4114,4465,4466,4467,4468,4469,4470,4471,4472,4473,4474,4475,4476,4477,4478,4479,4480,4481,4482,4483,4484,4485,4486,4487,4488,4489,4492,4493,4494,4495,4496,4497,4498,4499,4500,4501,4502,4503,4504,4505,4506,4507,4508,4509,4510,4511,4512,4513,4514,4515,4516,4866,4867,4868,4869,4870,4871,4872,4873,4874,4875,4876,4877,4878,4879,4880,4881,4882,4883,4884,4885,4886,4887,4888,4889,4890,4893,4894,4895,4896,4897,4898,4899,4900,4901,4902,4903,4904,4905,4906,4907,4908,4909,4910,4911,4912,4913,4914,4915,4916,4917,5266,5267,5268,5269,5270,5271,5272,5273,5274,5275,5276,5277,5278,5279,5280,5281,5282,5283,5284,5285,5286,5287,5288,5289,5290,5293,5294,5295,5296,5297,5298,5299,5300,5301,5302,5303,5304,5305,5306,5307,5308,5309,5310,5311,5312,5313,5314,5315,5316,5317',
       outframe='LSRK', gridder='mosaic', cycleniter=500, 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.T19t18.J142142+343000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.pb.tt0')

# hif_makermsimages(pipelinemode="automatic")
#
imdev(zscore=-1, ylength='60arcsec', stretch=False, interp='cubic',
      outfile='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.image.pbcor.tt0.rms',
      imagename='VLASS1.1.ql.T19t18.J142142+343000.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.T19t18.J142142+343000.10.2048.v1.I.iter1.image.pbcor.tt0')
imhead(imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.image.pbcor.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.image.pbcor.tt0.subim',
           imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.image.pbcor.tt0')
imstat(imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imhead(imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.image.pbcor.tt0.rms.subim',
           imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imstat(imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.psf.tt0')
imhead(imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.psf.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.psf.tt0.subim',
           imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.psf.tt0')
imstat(imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.residual.tt0')
imhead(imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.residual.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.residual.tt0.subim',
           imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.residual.tt0')
imstat(imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.pb.tt0')
imhead(imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.pb.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.pb.tt0.subim',
           imagename='VLASS1.1.ql.T19t18.J142142+343000.10.2048.v1.I.iter1.pb.tt0')

# hifv_exportvlassdata(pipelinemode="automatic")
#