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.sb34903570.eb35034891.58160.909213252315.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.sb34903570.eb35034891.58160.909213252315_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T05t07.J043916-213000_P24497v1_2018_04_02T20_29_17.853/rawdata/VLASS1.1.sb34903570.eb35034891.58160.909213252315.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='12482,12483,12480,12481,10299,10298,12484,12485,10295,10294,10297,10296,10291,10290,10293,10292,10699,10698,12079,12078,10691,10690,10693,10692,10695,10694,10697,10696,11119,11118,11111,11110,11113,11112,11115,11114,11117,11116,12488,12489,10288,10289,10725,10724,10721,10720,10723,10722,10688,10689,12477,10686,10687,10684,10685,11108,11109,12516,11100,11101,11106,11107,11104,11105,10321,10320,10323,10322,10325,10324,10327,10326,10329,10328,12514,11883,11882,11881,11880,11887,11886,11885,11884,11889,11888,12503,12502,12501,12500,12507,12506,12505,12504,12509,12508,12095,12094,12096,12091,12090,12093,12092,12099,12486,10709,10708,10703,10702,10701,10700,10707,10706,11120,11121,11089,11088,11083,11082,11896,11897,11087,11086,11085,11084,11519,11518,11515,11514,11517,11516,11511,11122,11513,11512,12086,12087,12084,11123,12082,12083,12080,12081,12088,12089,11094,11095,11096,11097,11090,11091,11092,11093,10307,10306,10305,10304,10303,10302,10301,10300,11894,11895,11508,11509,11506,11507,11504,11505,11502,11503,11891,11892,11893,10318,10319,12077,12517,10310,10311,10312,10313,10314,10315,10316,10317,11098,12518,12107,12106,12105,12104,12103,12102,12101,12100,12479,12478,12109,12108,10714,10715,10716,10717,10710,10711,10712,10713,12515,10718,10719,11099,12511,12510,11510,11520,11521,11489,11488,11487,11486,11485,11484,11483,11482,11481,11480,12110,12111,12112,12113,12114,12115,12116,12117,12118,12085,11890,11878,11879,12487,12513,12491,12490,12493,12492,12495,12494,12496,12499,11498,11499,11490,11491,11492,11493,11494,11495,11496,11497,12512',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 04:39:16.837 -21.30.0.0000',
       scan=['192,163,198,166,169,173,176,179,182,201,185,188,195'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34903570.eb35034891.58160.909213252315_split.ms'],
       imagename='VLASS1.1.ql.T05t07.J043916-213000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='10288,10289,10290,10291,10292,10293,10294,10295,10296,10297,10298,10299,10300,10301,10302,10303,10304,10305,10306,10307,10310,10311,10312,10313,10314,10315,10316,10317,10318,10319,10320,10321,10322,10323,10324,10325,10326,10327,10328,10329,10684,10685,10686,10687,10688,10689,10690,10691,10692,10693,10694,10695,10696,10697,10698,10699,10700,10701,10702,10703,10706,10707,10708,10709,10710,10711,10712,10713,10714,10715,10716,10717,10718,10719,10720,10721,10722,10723,10724,10725,11082,11083,11084,11085,11086,11087,11088,11089,11090,11091,11092,11093,11094,11095,11096,11097,11098,11099,11100,11101,11104,11105,11106,11107,11108,11109,11110,11111,11112,11113,11114,11115,11116,11117,11118,11119,11120,11121,11122,11123,11480,11481,11482,11483,11484,11485,11486,11487,11488,11489,11490,11491,11492,11493,11494,11495,11496,11497,11498,11499,11502,11503,11504,11505,11506,11507,11508,11509,11510,11511,11512,11513,11514,11515,11516,11517,11518,11519,11520,11521,11878,11879,11880,11881,11882,11883,11884,11885,11886,11887,11888,11889,11890,11891,11892,11893,11894,11895,11896,11897,12077,12078,12079,12080,12081,12082,12083,12084,12085,12086,12087,12088,12089,12090,12091,12092,12093,12094,12095,12096,12099,12100,12101,12102,12103,12104,12105,12106,12107,12108,12109,12110,12111,12112,12113,12114,12115,12116,12117,12118,12477,12478,12479,12480,12481,12482,12483,12484,12485,12486,12487,12488,12489,12490,12491,12492,12493,12494,12495,12496,12499,12500,12501,12502,12503,12504,12505,12506,12507,12508,12509,12510,12511,12512,12513,12514,12515,12516,12517,12518',
       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:39:16.837 -21.30.0.0000',
       scan=['192,163,198,166,169,173,176,179,182,201,185,188,195'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34903570.eb35034891.58160.909213252315_split.ms'],
       imagename='VLASS1.1.ql.T05t07.J043916-213000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00051139583626309099, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='10288,10289,10290,10291,10292,10293,10294,10295,10296,10297,10298,10299,10300,10301,10302,10303,10304,10305,10306,10307,10310,10311,10312,10313,10314,10315,10316,10317,10318,10319,10320,10321,10322,10323,10324,10325,10326,10327,10328,10329,10684,10685,10686,10687,10688,10689,10690,10691,10692,10693,10694,10695,10696,10697,10698,10699,10700,10701,10702,10703,10706,10707,10708,10709,10710,10711,10712,10713,10714,10715,10716,10717,10718,10719,10720,10721,10722,10723,10724,10725,11082,11083,11084,11085,11086,11087,11088,11089,11090,11091,11092,11093,11094,11095,11096,11097,11098,11099,11100,11101,11104,11105,11106,11107,11108,11109,11110,11111,11112,11113,11114,11115,11116,11117,11118,11119,11120,11121,11122,11123,11480,11481,11482,11483,11484,11485,11486,11487,11488,11489,11490,11491,11492,11493,11494,11495,11496,11497,11498,11499,11502,11503,11504,11505,11506,11507,11508,11509,11510,11511,11512,11513,11514,11515,11516,11517,11518,11519,11520,11521,11878,11879,11880,11881,11882,11883,11884,11885,11886,11887,11888,11889,11890,11891,11892,11893,11894,11895,11896,11897,12077,12078,12079,12080,12081,12082,12083,12084,12085,12086,12087,12088,12089,12090,12091,12092,12093,12094,12095,12096,12099,12100,12101,12102,12103,12104,12105,12106,12107,12108,12109,12110,12111,12112,12113,12114,12115,12116,12117,12118,12477,12478,12479,12480,12481,12482,12483,12484,12485,12486,12487,12488,12489,12490,12491,12492,12493,12494,12495,12496,12499,12500,12501,12502,12503,12504,12505,12506,12507,12508,12509,12510,12511,12512,12513,12514,12515,12516,12517,12518',
       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.T05t07.J043916-213000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T05t07.J043916-213000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T05t07.J043916-213000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#