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

# hif_editimlist(parameter_file='parameter.15312.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_T13t30.J192606+113000_P15312v1_2017_11_08T16_49_19.653/working/VLASS1.1.sb34609053.eb34623434.58051.07626767361_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T13t30.J192606+113000_P15312v1_2017_11_08T16_49_19.653/rawdata/VLASS1.1.sb34609053.eb34623434.58051.07626767361.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='5453,6233,6251,6249,5445,6572,6573,6570,6193,5379,5378,6574,6200,5375,5374,5377,5376,5373,4981,4980,4983,4982,4985,4984,4987,4986,4989,4988,6242,6243,6244,6245,6246,6247,5065,5064,5382,5066,5061,5060,5063,5062,5388,5389,5384,5460,5463,5462,5465,5464,5466,5838,6583,4648,4649,6569,6568,4642,4643,4640,4641,4646,4647,4644,4645,5854,5855,5856,5857,4978,4979,5852,5839,4974,4975,4976,4977,5858,5859,6255,4973,5393,5392,5391,5390,5397,5396,5395,5394,6598,6599,5399,5398,6581,4655,4654,4657,4656,4651,4650,4653,4652,6239,4659,4658,6586,6194,6262,6188,6260,5796,6196,5847,5846,5845,5844,5843,5842,5841,5840,5850,6195,5849,5848,5003,5002,5001,5000,6587,6192,6585,6584,6201,5790,6580,5779,5778,6571,6258,6178,6179,6176,6177,5773,5772,5775,5774,5777,5776,6253,6250,5403,5402,5401,5400,6591,6257,6170,6575,6199,6252,6197,6238,6578,6590,6579,5038,5039,5036,5037,6592,6256,6232,6593,4637,4636,6594,4639,4638,6174,6595,6248,6175,5436,5437,5444,6172,5438,5439,6173,5784,5785,5786,5787,5780,5781,5782,5783,5865,5864,6171,5788,5789,5863,5862,6596,5837,5385,6254,5834,6576,6240,5835,6241,6236,5861,6577,5860,5797,5443,5795,5794,5793,5792,5791,5442,6597,6198,5441,5799,5798,5380,5381,6589,5447,5446,5383,5058,5059,6261,5054,5055,5056,5057,5050,5051,5052,5053,5458,5459,5386,6588,5450,5451,5452,5387,5454,5455,5456,5457,6187,6186,6185,6184,6183,6182,6181,6180,5803,5802,5801,5800,6189,5851,6582,6259,6190,6263,5049,5048,5047,5046,5045,5044,5043,5042,5041,5040,5836,5461,5449,5448,6191,4664,4665,4666,5440,4660,4661,4662,4663,4996,4997,4994,4995,4992,4993,4990,4991,6235,6234,6237,5853,4998,4999',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 19:26:6.738 +11.30.0.0000',
       scan=['96,99,102,105,108,111,83,86,89,92'],
       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_T13t30.J192606+113000_P15312v1_2017_11_08T16_49_19.653/working/VLASS1.1.sb34609053.eb34623434.58051.07626767361_split.ms'],
       imagename='VLASS1.1.ql.T13t30.J192606+113000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='4636,4637,4638,4639,4640,4641,4642,4643,4644,4645,4646,4647,4648,4649,4650,4651,4652,4653,4654,4655,4656,4657,4658,4659,4660,4661,4662,4663,4664,4665,4666,4973,4974,4975,4976,4977,4978,4979,4980,4981,4982,4983,4984,4985,4986,4987,4988,4989,4990,4991,4992,4993,4994,4995,4996,4997,4998,4999,5000,5001,5002,5003,5036,5037,5038,5039,5040,5041,5042,5043,5044,5045,5046,5047,5048,5049,5050,5051,5052,5053,5054,5055,5056,5057,5058,5059,5060,5061,5062,5063,5064,5065,5066,5373,5374,5375,5376,5377,5378,5379,5380,5381,5382,5383,5384,5385,5386,5387,5388,5389,5390,5391,5392,5393,5394,5395,5396,5397,5398,5399,5400,5401,5402,5403,5436,5437,5438,5439,5440,5441,5442,5443,5444,5445,5446,5447,5448,5449,5450,5451,5452,5453,5454,5455,5456,5457,5458,5459,5460,5461,5462,5463,5464,5465,5466,5772,5773,5774,5775,5776,5777,5778,5779,5780,5781,5782,5783,5784,5785,5786,5787,5788,5789,5790,5791,5792,5793,5794,5795,5796,5797,5798,5799,5800,5801,5802,5803,5834,5835,5836,5837,5838,5839,5840,5841,5842,5843,5844,5845,5846,5847,5848,5849,5850,5851,5852,5853,5854,5855,5856,5857,5858,5859,5860,5861,5862,5863,5864,5865,6170,6171,6172,6173,6174,6175,6176,6177,6178,6179,6180,6181,6182,6183,6184,6185,6186,6187,6188,6189,6190,6191,6192,6193,6194,6195,6196,6197,6198,6199,6200,6201,6232,6233,6234,6235,6236,6237,6238,6239,6240,6241,6242,6243,6244,6245,6246,6247,6248,6249,6250,6251,6252,6253,6254,6255,6256,6257,6258,6259,6260,6261,6262,6263,6568,6569,6570,6571,6572,6573,6574,6575,6576,6577,6578,6579,6580,6581,6582,6583,6584,6585,6586,6587,6588,6589,6590,6591,6592,6593,6594,6595,6596,6597,6598,6599',
       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 19:26:6.738 +11.30.0.0000',
       scan=['96,99,102,105,108,111,83,86,89,92'],
       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_T13t30.J192606+113000_P15312v1_2017_11_08T16_49_19.653/working/VLASS1.1.sb34609053.eb34623434.58051.07626767361_split.ms'],
       imagename='VLASS1.1.ql.T13t30.J192606+113000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00052856442226948276, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='4636,4637,4638,4639,4640,4641,4642,4643,4644,4645,4646,4647,4648,4649,4650,4651,4652,4653,4654,4655,4656,4657,4658,4659,4660,4661,4662,4663,4664,4665,4666,4973,4974,4975,4976,4977,4978,4979,4980,4981,4982,4983,4984,4985,4986,4987,4988,4989,4990,4991,4992,4993,4994,4995,4996,4997,4998,4999,5000,5001,5002,5003,5036,5037,5038,5039,5040,5041,5042,5043,5044,5045,5046,5047,5048,5049,5050,5051,5052,5053,5054,5055,5056,5057,5058,5059,5060,5061,5062,5063,5064,5065,5066,5373,5374,5375,5376,5377,5378,5379,5380,5381,5382,5383,5384,5385,5386,5387,5388,5389,5390,5391,5392,5393,5394,5395,5396,5397,5398,5399,5400,5401,5402,5403,5436,5437,5438,5439,5440,5441,5442,5443,5444,5445,5446,5447,5448,5449,5450,5451,5452,5453,5454,5455,5456,5457,5458,5459,5460,5461,5462,5463,5464,5465,5466,5772,5773,5774,5775,5776,5777,5778,5779,5780,5781,5782,5783,5784,5785,5786,5787,5788,5789,5790,5791,5792,5793,5794,5795,5796,5797,5798,5799,5800,5801,5802,5803,5834,5835,5836,5837,5838,5839,5840,5841,5842,5843,5844,5845,5846,5847,5848,5849,5850,5851,5852,5853,5854,5855,5856,5857,5858,5859,5860,5861,5862,5863,5864,5865,6170,6171,6172,6173,6174,6175,6176,6177,6178,6179,6180,6181,6182,6183,6184,6185,6186,6187,6188,6189,6190,6191,6192,6193,6194,6195,6196,6197,6198,6199,6200,6201,6232,6233,6234,6235,6236,6237,6238,6239,6240,6241,6242,6243,6244,6245,6246,6247,6248,6249,6250,6251,6252,6253,6254,6255,6256,6257,6258,6259,6260,6261,6262,6263,6568,6569,6570,6571,6572,6573,6574,6575,6576,6577,6578,6579,6580,6581,6582,6583,6584,6585,6586,6587,6588,6589,6590,6591,6592,6593,6594,6595,6596,6597,6598,6599',
       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.T13t30.J192606+113000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T13t30.J192606+113000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T13t30.J192606+113000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#