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.sb34617180.eb34646687.58058.06139675926.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.sb34617180.eb34646687.58058.06139675926_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T08t32.J205810-083000_P16710v1_2018_04_25T14_25_14.034/rawdata/VLASS1.1.sb34617180.eb34646687.58058.06139675926.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='6501,5310,5317,5316,6723,6091,6721,6720,6727,5315,5508,5509,5506,5507,5504,5314,5502,5503,5500,5501,4909,4908,4705,4704,4703,4702,4701,4700,4901,4900,4903,4902,4905,4904,4907,4906,6327,5903,6325,6324,5906,5907,5904,5905,5908,5909,6329,6328,6716,5711,6109,6728,6515,6112,5505,6098,6718,6719,5511,5510,5513,5512,6095,6094,6097,6096,5308,5309,5710,5719,6319,5300,5301,5302,5303,5304,5305,5306,5307,6312,6114,6310,6311,5919,5918,6107,6115,5915,5914,5917,5916,5911,5910,5913,5912,4899,4898,6722,6088,6110,6089,6488,5729,6111,6100,6086,5498,5499,6087,5494,5495,5496,5497,5490,5491,5492,5493,5098,5099,5319,5318,6313,5090,5091,5092,5093,5094,5095,5096,5097,6309,6308,6305,6316,6307,6306,6712,6105,4707,6314,4706,6315,6317,6512,5326,5327,5324,5325,5322,5323,5320,5321,5487,5486,5485,5484,5483,6104,5328,5329,6714,5089,5088,4679,5083,5082,5081,6500,5087,5086,5085,5084,6733,6090,6708,6734,6103,6726,6735,6737,6335,5489,6717,6511,5488,6106,6707,5330,6489,6093,4709,6102,6725,4708,6490,5722,5708,5709,6326,5704,5705,5706,5707,6508,6116,5702,5703,6092,6318,6101,6724,6323,6322,6510,5723,6715,6513,6321,6491,6320,6509,5108,5109,4924,5102,5103,5100,5101,5106,5107,5104,5105,4686,4687,4684,4685,4682,4683,4680,4681,5717,5716,5715,5714,5713,5712,4688,4689,6503,5718,6502,5721,6713,6495,6494,6497,6496,5111,5110,6493,6492,6507,6709,6499,6498,6506,6729,6505,5728,6504,6710,6516,4691,4922,4693,4692,4695,4926,4925,4696,4699,4698,5720,4928,5726,5727,5724,5725,5924,5925,5926,5927,5920,5921,5922,5923,6330,6711,5928,5929,6331,4923,6113,6332,4690,6517,5930,4921,6730,6731,6732,4920,5731,5730,6736,5732,6514,4927,4694,6108,6099,4697,4918,4919,4916,4917,4914,4915,4912,4913,4910,4911,6334,5313,6333,5933,5932,5931,5312,6518,5311',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 20:58:10.614 -08.30.0.0000',
       scan=['98,102,105,108,114,111,82,86,89,92,95'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34617180.eb34646687.58058.06139675926_split.ms'],
       imagename='VLASS1.1.ql.T08t32.J205810-083000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='4679,4680,4681,4682,4683,4684,4685,4686,4687,4688,4689,4690,4691,4692,4693,4694,4695,4696,4697,4698,4699,4700,4701,4702,4703,4704,4705,4706,4707,4708,4709,4898,4899,4900,4901,4902,4903,4904,4905,4906,4907,4908,4909,4910,4911,4912,4913,4914,4915,4916,4917,4918,4919,4920,4921,4922,4923,4924,4925,4926,4927,4928,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,5110,5111,5300,5301,5302,5303,5304,5305,5306,5307,5308,5309,5310,5311,5312,5313,5314,5315,5316,5317,5318,5319,5320,5321,5322,5323,5324,5325,5326,5327,5328,5329,5330,5483,5484,5485,5486,5487,5488,5489,5490,5491,5492,5493,5494,5495,5496,5497,5498,5499,5500,5501,5502,5503,5504,5505,5506,5507,5508,5509,5510,5511,5512,5513,5702,5703,5704,5705,5706,5707,5708,5709,5710,5711,5712,5713,5714,5715,5716,5717,5718,5719,5720,5721,5722,5723,5724,5725,5726,5727,5728,5729,5730,5731,5732,5903,5904,5905,5906,5907,5908,5909,5910,5911,5912,5913,5914,5915,5916,5917,5918,5919,5920,5921,5922,5923,5924,5925,5926,5927,5928,5929,5930,5931,5932,5933,6086,6087,6088,6089,6090,6091,6092,6093,6094,6095,6096,6097,6098,6099,6100,6101,6102,6103,6104,6105,6106,6107,6108,6109,6110,6111,6112,6113,6114,6115,6116,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,6333,6334,6335,6488,6489,6490,6491,6492,6493,6494,6495,6496,6497,6498,6499,6500,6501,6502,6503,6504,6505,6506,6507,6508,6509,6510,6511,6512,6513,6514,6515,6516,6517,6518,6707,6708,6709,6710,6711,6712,6713,6714,6715,6716,6717,6718,6719,6720,6721,6722,6723,6724,6725,6726,6727,6728,6729,6730,6731,6732,6733,6734,6735,6736,6737',
       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 20:58:10.614 -08.30.0.0000',
       scan=['98,102,105,108,114,111,82,86,89,92,95'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34617180.eb34646687.58058.06139675926_split.ms'],
       imagename='VLASS1.1.ql.T08t32.J205810-083000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00059003723060559415, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='4679,4680,4681,4682,4683,4684,4685,4686,4687,4688,4689,4690,4691,4692,4693,4694,4695,4696,4697,4698,4699,4700,4701,4702,4703,4704,4705,4706,4707,4708,4709,4898,4899,4900,4901,4902,4903,4904,4905,4906,4907,4908,4909,4910,4911,4912,4913,4914,4915,4916,4917,4918,4919,4920,4921,4922,4923,4924,4925,4926,4927,4928,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,5110,5111,5300,5301,5302,5303,5304,5305,5306,5307,5308,5309,5310,5311,5312,5313,5314,5315,5316,5317,5318,5319,5320,5321,5322,5323,5324,5325,5326,5327,5328,5329,5330,5483,5484,5485,5486,5487,5488,5489,5490,5491,5492,5493,5494,5495,5496,5497,5498,5499,5500,5501,5502,5503,5504,5505,5506,5507,5508,5509,5510,5511,5512,5513,5702,5703,5704,5705,5706,5707,5708,5709,5710,5711,5712,5713,5714,5715,5716,5717,5718,5719,5720,5721,5722,5723,5724,5725,5726,5727,5728,5729,5730,5731,5732,5903,5904,5905,5906,5907,5908,5909,5910,5911,5912,5913,5914,5915,5916,5917,5918,5919,5920,5921,5922,5923,5924,5925,5926,5927,5928,5929,5930,5931,5932,5933,6086,6087,6088,6089,6090,6091,6092,6093,6094,6095,6096,6097,6098,6099,6100,6101,6102,6103,6104,6105,6106,6107,6108,6109,6110,6111,6112,6113,6114,6115,6116,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,6333,6334,6335,6488,6489,6490,6491,6492,6493,6494,6495,6496,6497,6498,6499,6500,6501,6502,6503,6504,6505,6506,6507,6508,6509,6510,6511,6512,6513,6514,6515,6516,6517,6518,6707,6708,6709,6710,6711,6712,6713,6714,6715,6716,6717,6718,6719,6720,6721,6722,6723,6724,6725,6726,6727,6728,6729,6730,6731,6732,6733,6734,6735,6736,6737',
       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.T08t32.J205810-083000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T08t32.J205810-083000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T08t32.J205810-083000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#