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

# hif_editimlist(parameter_file='parameter.14595.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_T28t03.J045013+693000_P14595v1_2017_11_12T21_11_05.704/working/VLASS1.1.sb34589996.eb34600351.58041.24583186343_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T28t03.J045013+693000_P14595v1_2017_11_12T21_11_05.704/rawdata/VLASS1.1.sb34589996.eb34600351.58041.24583186343.ms',
            datacolumn='corrected', reindex=False,
            spw='0,1,10,11,12,13,14,15,2,3,4,5,6,7,8,9',
            field='5821,5820,5988,5989,5825,5824,5982,5983,5980,5981,5986,5987,5984,5985,4936,6659,4937,4938,4526,5797,5796,5823,5822,5799,5798,6437,6436,6435,6434,6433,6432,6431,6430,4770,5182,5183,5180,5181,5186,5187,5184,5185,5188,5189,5818,5819,5795,5810,5811,5812,5813,5814,5815,5816,5817,5357,5356,5355,5354,5353,5352,5351,5350,5359,5358,4505,4504,4507,4506,4501,4500,4503,4502,4509,4508,4918,4919,4916,4917,4915,5173,5172,5171,5170,5177,5176,5175,5174,5179,5178,5586,5587,5584,5585,5582,5583,5581,5588,5589,6248,6249,6240,6241,6242,6243,6244,6245,6246,6247,5995,5994,5997,5996,5991,5990,5993,5992,5999,5998,6668,6669,6662,6663,6660,6661,6666,6667,6664,6665,4945,4944,4941,4940,4943,4942,6406,6407,6408,6409,4761,4760,4763,4762,4765,4764,4767,4766,4769,4768,5809,5808,5803,5802,5801,5800,5807,5806,5805,5804,5362,5363,5360,5361,5364,5365,6004,6005,6002,6003,6000,6001,5591,5590,5593,5592,5595,5594,5597,5596,5599,5598,5610,5611,6253,6252,6251,6250,6255,6254,4499,6675,6674,6677,6676,6671,6670,6673,6672,6679,6678,4758,4759,4750,4751,4752,4753,4754,4755,4756,4757,5979,5978,5977,5976,5975,5974,6688,6689,6680,6681,6682,6683,6684,6685,6686,6687,4930,4931,4932,4933,4934,4935,4529,4528,4527,4939,4525,4524,4523,4522,4521,4520,6419,6418,6415,6414,6417,6416,6411,6410,6413,6412,6226,6227,6225,6228,6229,4518,4920,4929,4511,5159,5158,5605,5604,5607,5606,5601,5600,5603,5602,5609,5608,6424,6425,6426,6427,6420,6421,6422,6423,6428,6429,5335,5337,5336,5339,5338,4749,4748,4743,4742,4741,4740,4747,4746,4745,4744,5344,5345,5346,5347,5340,5341,5342,5343,5348,5349,6690,4923,4922,4921,4519,4927,4926,4925,4924,4512,4513,4510,4928,4516,4517,4514,4515,5160,5161,5162,5163,5164,5165,5166,5167,5168,5169,6239,6238,6235,6234,6237,6236,6231,6230,6233,6232',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 04:50:13.662 +69.30.0.0000',
       scan=['97,101,113,104,107,110,81,85,88,91,94'],
       spw=['0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15'],
       vis=['/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T28t03.J045013+693000_P14595v1_2017_11_12T21_11_05.704/working/VLASS1.1.sb34589996.eb34600351.58041.24583186343_split.ms'],
       imagename='VLASS1.1.ql.T28t03.J045013+693000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='4499,4500,4501,4502,4503,4504,4505,4506,4507,4508,4509,4510,4511,4512,4513,4514,4515,4516,4517,4518,4519,4520,4521,4522,4523,4524,4525,4526,4527,4528,4529,4740,4741,4742,4743,4744,4745,4746,4747,4748,4749,4750,4751,4752,4753,4754,4755,4756,4757,4758,4759,4760,4761,4762,4763,4764,4765,4766,4767,4768,4769,4770,4915,4916,4917,4918,4919,4920,4921,4922,4923,4924,4925,4926,4927,4928,4929,4930,4931,4932,4933,4934,4935,4936,4937,4938,4939,4940,4941,4942,4943,4944,4945,5158,5159,5160,5161,5162,5163,5164,5165,5166,5167,5168,5169,5170,5171,5172,5173,5174,5175,5176,5177,5178,5179,5180,5181,5182,5183,5184,5185,5186,5187,5188,5189,5335,5336,5337,5338,5339,5340,5341,5342,5343,5344,5345,5346,5347,5348,5349,5350,5351,5352,5353,5354,5355,5356,5357,5358,5359,5360,5361,5362,5363,5364,5365,5581,5582,5583,5584,5585,5586,5587,5588,5589,5590,5591,5592,5593,5594,5595,5596,5597,5598,5599,5600,5601,5602,5603,5604,5605,5606,5607,5608,5609,5610,5611,5795,5796,5797,5798,5799,5800,5801,5802,5803,5804,5805,5806,5807,5808,5809,5810,5811,5812,5813,5814,5815,5816,5817,5818,5819,5820,5821,5822,5823,5824,5825,5974,5975,5976,5977,5978,5979,5980,5981,5982,5983,5984,5985,5986,5987,5988,5989,5990,5991,5992,5993,5994,5995,5996,5997,5998,5999,6000,6001,6002,6003,6004,6005,6225,6226,6227,6228,6229,6230,6231,6232,6233,6234,6235,6236,6237,6238,6239,6240,6241,6242,6243,6244,6245,6246,6247,6248,6249,6250,6251,6252,6253,6254,6255,6406,6407,6408,6409,6410,6411,6412,6413,6414,6415,6416,6417,6418,6419,6420,6421,6422,6423,6424,6425,6426,6427,6428,6429,6430,6431,6432,6433,6434,6435,6436,6437,6659,6660,6661,6662,6663,6664,6665,6666,6667,6668,6669,6670,6671,6672,6673,6674,6675,6676,6677,6678,6679,6680,6681,6682,6683,6684,6685,6686,6687,6688,6689,6690',
       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 04:50:13.662 +69.30.0.0000',
       scan=['97,101,113,104,107,110,81,85,88,91,94'],
       spw=['0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15'],
       vis=['/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T28t03.J045013+693000_P14595v1_2017_11_12T21_11_05.704/working/VLASS1.1.sb34589996.eb34600351.58041.24583186343_split.ms'],
       imagename='VLASS1.1.ql.T28t03.J045013+693000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00048485624266031663, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='4499,4500,4501,4502,4503,4504,4505,4506,4507,4508,4509,4510,4511,4512,4513,4514,4515,4516,4517,4518,4519,4520,4521,4522,4523,4524,4525,4526,4527,4528,4529,4740,4741,4742,4743,4744,4745,4746,4747,4748,4749,4750,4751,4752,4753,4754,4755,4756,4757,4758,4759,4760,4761,4762,4763,4764,4765,4766,4767,4768,4769,4770,4915,4916,4917,4918,4919,4920,4921,4922,4923,4924,4925,4926,4927,4928,4929,4930,4931,4932,4933,4934,4935,4936,4937,4938,4939,4940,4941,4942,4943,4944,4945,5158,5159,5160,5161,5162,5163,5164,5165,5166,5167,5168,5169,5170,5171,5172,5173,5174,5175,5176,5177,5178,5179,5180,5181,5182,5183,5184,5185,5186,5187,5188,5189,5335,5336,5337,5338,5339,5340,5341,5342,5343,5344,5345,5346,5347,5348,5349,5350,5351,5352,5353,5354,5355,5356,5357,5358,5359,5360,5361,5362,5363,5364,5365,5581,5582,5583,5584,5585,5586,5587,5588,5589,5590,5591,5592,5593,5594,5595,5596,5597,5598,5599,5600,5601,5602,5603,5604,5605,5606,5607,5608,5609,5610,5611,5795,5796,5797,5798,5799,5800,5801,5802,5803,5804,5805,5806,5807,5808,5809,5810,5811,5812,5813,5814,5815,5816,5817,5818,5819,5820,5821,5822,5823,5824,5825,5974,5975,5976,5977,5978,5979,5980,5981,5982,5983,5984,5985,5986,5987,5988,5989,5990,5991,5992,5993,5994,5995,5996,5997,5998,5999,6000,6001,6002,6003,6004,6005,6225,6226,6227,6228,6229,6230,6231,6232,6233,6234,6235,6236,6237,6238,6239,6240,6241,6242,6243,6244,6245,6246,6247,6248,6249,6250,6251,6252,6253,6254,6255,6406,6407,6408,6409,6410,6411,6412,6413,6414,6415,6416,6417,6418,6419,6420,6421,6422,6423,6424,6425,6426,6427,6428,6429,6430,6431,6432,6433,6434,6435,6436,6437,6659,6660,6661,6662,6663,6664,6665,6666,6667,6668,6669,6670,6671,6672,6673,6674,6675,6676,6677,6678,6679,6680,6681,6682,6683,6684,6685,6686,6687,6688,6689,6690',
       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.T28t03.J045013+693000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T28t03.J045013+693000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T28t03.J045013+693000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#