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

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

# hif_transformimagedata(datacolumn='corrected', modify_weights=False, clear_pointing=True)
#
mstransform(timebin='0s',
            outputvis='VLASS1.1.sb34900388.eb35074823.58165.26784681713_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T06t16.J101433-163000_P24701v1_2018_03_07T01_11_10.709/rawdata/VLASS1.1.sb34900388.eb35074823.58165.26784681713.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='6074,6312,4570,4571,4572,4573,4574,4575,4576,4577,6084,6082,6083,6080,6081,4707,4706,4705,4704,4703,4702,4701,4700,5684,4709,4708,6327,6308,6325,6067,6323,6322,6321,6320,6060,6329,6061,6310,6332,5923,6306,6305,4563,4562,4561,4560,4567,4566,4565,4564,4569,4568,5687,5686,5342,6328,6319,5681,5343,5680,5692,5693,5690,5079,5919,5683,5694,6315,5915,5667,5917,5682,5911,5913,5666,6316,5665,6302,5664,6326,4558,4559,4556,4557,4554,4555,4552,4553,4550,4551,5098,5099,6313,5090,5091,5092,5093,5094,5095,5096,5097,6064,6065,6066,5691,5689,5688,6062,6063,5685,4966,4965,4964,4963,4962,4961,4960,6314,6307,5671,5548,6317,4549,4548,5672,5549,5674,5675,5676,5677,5678,5679,6318,4547,4546,5326,5089,5088,5940,5327,5083,5082,5081,5080,5087,5086,5085,5084,6077,6076,4958,4959,6073,6072,6071,6070,4952,4953,4950,4951,4956,4957,4954,4955,6069,6304,5673,5351,5546,5547,5331,5330,5333,5332,5335,5334,5337,5336,5339,5338,5669,5668,5670,4945,4944,4947,4946,4941,4940,4943,4942,5542,5543,5540,5541,4949,4948,5544,5545,6324,5328,5329,5918,6079,5108,5109,5551,5550,5102,5103,5100,5101,5106,5107,5104,5105,5344,5345,5346,5347,5340,5341,4936,4937,4938,4939,5348,5349,4689,6059,6058,6055,6054,6057,6056,6068,5521,5522,5523,5524,5525,5526,5527,5528,5529,6303,4691,4690,4693,4692,4695,4694,4697,4696,4699,4698,6075,5924,5925,5926,5927,5920,5921,5922,5914,6330,5912,5928,5929,6331,6309,5941,5356,5355,5530,5533,5532,5531,5354,5537,5536,5535,5534,5353,5539,5538,5352,4714,4715,4716,4717,4710,4711,4712,4713,5934,6311,5350,4718,4719,5937,5936,5935,5916,5933,5932,5931,5930,6078,5939,5938',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 10:14:33.852 -16.30.0.0000',
       scan=['97,100,103,106,109,81,84,87,90,93'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34900388.eb35074823.58165.26784681713_split.ms'],
       imagename='VLASS1.1.ql.T06t16.J101433-163000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='4546,4547,4548,4549,4550,4551,4552,4553,4554,4555,4556,4557,4558,4559,4560,4561,4562,4563,4564,4565,4566,4567,4568,4569,4570,4571,4572,4573,4574,4575,4576,4577,4689,4690,4691,4692,4693,4694,4695,4696,4697,4698,4699,4700,4701,4702,4703,4704,4705,4706,4707,4708,4709,4710,4711,4712,4713,4714,4715,4716,4717,4718,4719,4936,4937,4938,4939,4940,4941,4942,4943,4944,4945,4946,4947,4948,4949,4950,4951,4952,4953,4954,4955,4956,4957,4958,4959,4960,4961,4962,4963,4964,4965,4966,5079,5080,5081,5082,5083,5084,5085,5086,5087,5088,5089,5090,5091,5092,5093,5094,5095,5096,5097,5098,5099,5100,5101,5102,5103,5104,5105,5106,5107,5108,5109,5326,5327,5328,5329,5330,5331,5332,5333,5334,5335,5336,5337,5338,5339,5340,5341,5342,5343,5344,5345,5346,5347,5348,5349,5350,5351,5352,5353,5354,5355,5356,5521,5522,5523,5524,5525,5526,5527,5528,5529,5530,5531,5532,5533,5534,5535,5536,5537,5538,5539,5540,5541,5542,5543,5544,5545,5546,5547,5548,5549,5550,5551,5664,5665,5666,5667,5668,5669,5670,5671,5672,5673,5674,5675,5676,5677,5678,5679,5680,5681,5682,5683,5684,5685,5686,5687,5688,5689,5690,5691,5692,5693,5694,5911,5912,5913,5914,5915,5916,5917,5918,5919,5920,5921,5922,5923,5924,5925,5926,5927,5928,5929,5930,5931,5932,5933,5934,5935,5936,5937,5938,5939,5940,5941,6054,6055,6056,6057,6058,6059,6060,6061,6062,6063,6064,6065,6066,6067,6068,6069,6070,6071,6072,6073,6074,6075,6076,6077,6078,6079,6080,6081,6082,6083,6084,6302,6303,6304,6305,6306,6307,6308,6309,6310,6311,6312,6313,6314,6315,6316,6317,6318,6319,6320,6321,6322,6323,6324,6325,6326,6327,6328,6329,6330,6331,6332',
       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 10:14:33.852 -16.30.0.0000',
       scan=['97,100,103,106,109,81,84,87,90,93'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34900388.eb35074823.58165.26784681713_split.ms'],
       imagename='VLASS1.1.ql.T06t16.J101433-163000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00053667305415146986, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='4546,4547,4548,4549,4550,4551,4552,4553,4554,4555,4556,4557,4558,4559,4560,4561,4562,4563,4564,4565,4566,4567,4568,4569,4570,4571,4572,4573,4574,4575,4576,4577,4689,4690,4691,4692,4693,4694,4695,4696,4697,4698,4699,4700,4701,4702,4703,4704,4705,4706,4707,4708,4709,4710,4711,4712,4713,4714,4715,4716,4717,4718,4719,4936,4937,4938,4939,4940,4941,4942,4943,4944,4945,4946,4947,4948,4949,4950,4951,4952,4953,4954,4955,4956,4957,4958,4959,4960,4961,4962,4963,4964,4965,4966,5079,5080,5081,5082,5083,5084,5085,5086,5087,5088,5089,5090,5091,5092,5093,5094,5095,5096,5097,5098,5099,5100,5101,5102,5103,5104,5105,5106,5107,5108,5109,5326,5327,5328,5329,5330,5331,5332,5333,5334,5335,5336,5337,5338,5339,5340,5341,5342,5343,5344,5345,5346,5347,5348,5349,5350,5351,5352,5353,5354,5355,5356,5521,5522,5523,5524,5525,5526,5527,5528,5529,5530,5531,5532,5533,5534,5535,5536,5537,5538,5539,5540,5541,5542,5543,5544,5545,5546,5547,5548,5549,5550,5551,5664,5665,5666,5667,5668,5669,5670,5671,5672,5673,5674,5675,5676,5677,5678,5679,5680,5681,5682,5683,5684,5685,5686,5687,5688,5689,5690,5691,5692,5693,5694,5911,5912,5913,5914,5915,5916,5917,5918,5919,5920,5921,5922,5923,5924,5925,5926,5927,5928,5929,5930,5931,5932,5933,5934,5935,5936,5937,5938,5939,5940,5941,6054,6055,6056,6057,6058,6059,6060,6061,6062,6063,6064,6065,6066,6067,6068,6069,6070,6071,6072,6073,6074,6075,6076,6077,6078,6079,6080,6081,6082,6083,6084,6302,6303,6304,6305,6306,6307,6308,6309,6310,6311,6312,6313,6314,6315,6316,6317,6318,6319,6320,6321,6322,6323,6324,6325,6326,6327,6328,6329,6330,6331,6332',
       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.T06t16.J101433-163000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T06t16.J101433-163000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T06t16.J101433-163000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#