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

# hif_editimlist(parameter_file='parameter.5184.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_T18t22.J154715+283000_P5184v1_2018_01_12T20_55_50.449/working/VLASS1.1.sb34396276.eb34742211.58081.44143899306_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T18t22.J154715+283000_P5184v1_2018_01_12T20_55_50.449/rawdata/VLASS1.1.sb34396276.eb34742211.58081.44143899306.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='19490,19491,18495,18494,18493,18492,18491,18490,18093,18092,18091,18090,19492,18095,18094,19493,19055,19494,19863,19495,20076,20074,20075,20072,19496,19856,20071,19849,19497,20079,19498,19499,18480,18481,18482,18483,18484,18485,18486,18487,18488,18489,19467,19466,19461,19460,19463,19462,19869,19868,19861,19860,19447,19862,19865,19864,19867,19866,19061,19060,19063,19062,19065,19064,19067,19066,19069,19068,20093,20455,18479,18478,19478,19479,18475,18474,18477,18476,18471,18470,18473,18472,18079,18078,20058,18071,18070,18073,18072,18075,18074,18077,18076,19072,19073,19076,19077,19074,19075,20054,20055,19078,19079,19465,20052,20053,18084,18085,18086,18087,18080,18081,18082,18083,18088,18089,19446,19448,18466,18467,18464,18465,18462,18463,18460,18461,19489,20480,20459,20104,19093,20102,20103,20100,20073,18862,18863,18860,18861,18866,18867,18864,18865,19047,19046,18868,19044,19459,19050,19475,19051,19488,19052,19850,19851,18459,18458,19854,19855,19853,19857,18453,18452,18451,18450,18457,18456,18455,18454,20460,19056,20473,19057,19058,19059,18859,18858,18857,18856,18855,18854,18853,18852,18851,18850,19476,20098,19477,20101,19474,19464,20057,20089,20088,20087,20086,20085,20084,20083,20082,20081,20080,20099,19470,20475,19471,18848,18849,19049,20468,18843,18844,18845,18846,18847,19454,19455,19456,19457,19450,19451,19452,19453,19449,19458,19045,19083,19082,19081,19080,19087,19086,19085,19084,20090,20091,19089,19088,20094,20095,20096,20097,20469,19872,19873,20461,19469,20463,20462,20465,19870,20467,19468,20051,19871,20456,20457,20070,20464,18444,18445,18446,18447,19874,18442,18443,20466,18448,18449,19487,19486,19485,19484,19483,19482,19481,19480,19094,19095,19096,19097,19090,19091,19092,19048,20478,20479,20472,19053,20470,20471,20476,20477,20474,20059,20065,20064,20067,20066,20061,20060,20063,20062,20092,19852,20069,20068,19858,19054,19859,20458,20056',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 15:47:15.908 +28.30.0.0000',
       scan=['328,321,331,325,296,299,302,305,309,312,315,318'],
       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_T18t22.J154715+283000_P5184v1_2018_01_12T20_55_50.449/working/VLASS1.1.sb34396276.eb34742211.58081.44143899306_split.ms'],
       imagename='VLASS1.1.ql.T18t22.J154715+283000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='18070,18071,18072,18073,18074,18075,18076,18077,18078,18079,18080,18081,18082,18083,18084,18085,18086,18087,18088,18089,18090,18091,18092,18093,18094,18095,18442,18443,18444,18445,18446,18447,18448,18449,18450,18451,18452,18453,18454,18455,18456,18457,18458,18459,18460,18461,18462,18463,18464,18465,18466,18467,18470,18471,18472,18473,18474,18475,18476,18477,18478,18479,18480,18481,18482,18483,18484,18485,18486,18487,18488,18489,18490,18491,18492,18493,18494,18495,18843,18844,18845,18846,18847,18848,18849,18850,18851,18852,18853,18854,18855,18856,18857,18858,18859,18860,18861,18862,18863,18864,18865,18866,18867,18868,19044,19045,19046,19047,19048,19049,19050,19051,19052,19053,19054,19055,19056,19057,19058,19059,19060,19061,19062,19063,19064,19065,19066,19067,19068,19069,19072,19073,19074,19075,19076,19077,19078,19079,19080,19081,19082,19083,19084,19085,19086,19087,19088,19089,19090,19091,19092,19093,19094,19095,19096,19097,19446,19447,19448,19449,19450,19451,19452,19453,19454,19455,19456,19457,19458,19459,19460,19461,19462,19463,19464,19465,19466,19467,19468,19469,19470,19471,19474,19475,19476,19477,19478,19479,19480,19481,19482,19483,19484,19485,19486,19487,19488,19489,19490,19491,19492,19493,19494,19495,19496,19497,19498,19499,19849,19850,19851,19852,19853,19854,19855,19856,19857,19858,19859,19860,19861,19862,19863,19864,19865,19866,19867,19868,19869,19870,19871,19872,19873,19874,20051,20052,20053,20054,20055,20056,20057,20058,20059,20060,20061,20062,20063,20064,20065,20066,20067,20068,20069,20070,20071,20072,20073,20074,20075,20076,20079,20080,20081,20082,20083,20084,20085,20086,20087,20088,20089,20090,20091,20092,20093,20094,20095,20096,20097,20098,20099,20100,20101,20102,20103,20104,20455,20456,20457,20458,20459,20460,20461,20462,20463,20464,20465,20466,20467,20468,20469,20470,20471,20472,20473,20474,20475,20476,20477,20478,20479,20480',
       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 15:47:15.908 +28.30.0.0000',
       scan=['328,321,331,325,296,299,302,305,309,312,315,318'],
       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_T18t22.J154715+283000_P5184v1_2018_01_12T20_55_50.449/working/VLASS1.1.sb34396276.eb34742211.58081.44143899306_split.ms'],
       imagename='VLASS1.1.ql.T18t22.J154715+283000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00047316229444222696, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='18070,18071,18072,18073,18074,18075,18076,18077,18078,18079,18080,18081,18082,18083,18084,18085,18086,18087,18088,18089,18090,18091,18092,18093,18094,18095,18442,18443,18444,18445,18446,18447,18448,18449,18450,18451,18452,18453,18454,18455,18456,18457,18458,18459,18460,18461,18462,18463,18464,18465,18466,18467,18470,18471,18472,18473,18474,18475,18476,18477,18478,18479,18480,18481,18482,18483,18484,18485,18486,18487,18488,18489,18490,18491,18492,18493,18494,18495,18843,18844,18845,18846,18847,18848,18849,18850,18851,18852,18853,18854,18855,18856,18857,18858,18859,18860,18861,18862,18863,18864,18865,18866,18867,18868,19044,19045,19046,19047,19048,19049,19050,19051,19052,19053,19054,19055,19056,19057,19058,19059,19060,19061,19062,19063,19064,19065,19066,19067,19068,19069,19072,19073,19074,19075,19076,19077,19078,19079,19080,19081,19082,19083,19084,19085,19086,19087,19088,19089,19090,19091,19092,19093,19094,19095,19096,19097,19446,19447,19448,19449,19450,19451,19452,19453,19454,19455,19456,19457,19458,19459,19460,19461,19462,19463,19464,19465,19466,19467,19468,19469,19470,19471,19474,19475,19476,19477,19478,19479,19480,19481,19482,19483,19484,19485,19486,19487,19488,19489,19490,19491,19492,19493,19494,19495,19496,19497,19498,19499,19849,19850,19851,19852,19853,19854,19855,19856,19857,19858,19859,19860,19861,19862,19863,19864,19865,19866,19867,19868,19869,19870,19871,19872,19873,19874,20051,20052,20053,20054,20055,20056,20057,20058,20059,20060,20061,20062,20063,20064,20065,20066,20067,20068,20069,20070,20071,20072,20073,20074,20075,20076,20079,20080,20081,20082,20083,20084,20085,20086,20087,20088,20089,20090,20091,20092,20093,20094,20095,20096,20097,20098,20099,20100,20101,20102,20103,20104,20455,20456,20457,20458,20459,20460,20461,20462,20463,20464,20465,20466,20467,20468,20469,20470,20471,20472,20473,20474,20475,20476,20477,20478,20479,20480',
       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.T18t22.J154715+283000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T18t22.J154715+283000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T18t22.J154715+283000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#