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.sb34897149.eb35035193.58161.27141606482.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.sb34897149.eb35035193.58161.27141606482_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T03t25.J162509-293000_P25193v1_2018_04_23T14_48_32.579/rawdata/VLASS1.1.sb34897149.eb35035193.58161.27141606482.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='23838,22469,22770,24239,24058,24059,24066,23855,23854,23857,23856,23851,22771,23853,23852,23859,23858,23846,24220,22760,22761,22762,22763,22764,22765,22766,22767,22768,22769,22472,22473,22474,22475,22476,22477,24050,24051,24052,24053,23569,23568,24056,24057,23565,23564,23567,23566,23561,23560,23563,23562,23836,23847,23844,23845,23842,23843,23840,23841,23848,23849,24230,24216,24217,22965,24214,24215,23578,23579,23576,23577,23574,23575,23572,23573,23570,23571,22485,22484,22487,22486,22481,22480,22483,22482,22489,22488,23424,24234,24212,23415,23414,23417,23416,23411,23410,23413,23412,23419,23418,23850,22980,22981,22494,22983,22984,22985,22490,22987,22988,22989,23860,23585,23584,23149,23148,23147,23146,23145,23144,23143,23142,23141,23140,24231,23129,23128,23125,23124,23127,23126,23406,23407,23404,23405,23402,23403,23400,23401,24218,24219,23408,23409,22979,22978,22975,22974,22977,22976,22971,22970,22973,22972,23150,23151,23398,23399,22496,23138,23139,23136,23137,23134,23135,23132,22982,23130,23131,24227,22495,24229,24228,24213,22492,24225,24224,24223,24222,24221,22493,23837,22986,22968,22969,22966,22967,22964,22491,22962,22963,24226,24235,22478,22479,22759,22758,23835,23834,23833,24236,23425,22751,22750,22753,22752,22755,22754,22757,22756,24069,24068,24233,24061,23420,24063,24062,24065,24064,24067,23421,22470,24237,23422,22471,23423,24232,23133,23587,23586,22748,22749,23583,23582,23581,23580,22743,22746,22747,22744,22745,23839,24072,24073,24070,24071,24076,24077,24074,24075,24238,24060,24054,24055',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 16:25:9.227 -29.30.0.0000',
       scan=['357,360,364,367,370,373,376,380,383'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34897149.eb35035193.58161.27141606482_split.ms'],
       imagename='VLASS1.1.ql.T03t25.J162509-293000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[5760, 5760], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='22469,22470,22471,22472,22473,22474,22475,22476,22477,22478,22479,22480,22481,22482,22483,22484,22485,22486,22487,22488,22489,22490,22491,22492,22493,22494,22495,22496,22743,22744,22745,22746,22747,22748,22749,22750,22751,22752,22753,22754,22755,22756,22757,22758,22759,22760,22761,22762,22763,22764,22765,22766,22767,22768,22769,22770,22771,22962,22963,22964,22965,22966,22967,22968,22969,22970,22971,22972,22973,22974,22975,22976,22977,22978,22979,22980,22981,22982,22983,22984,22985,22986,22987,22988,22989,23124,23125,23126,23127,23128,23129,23130,23131,23132,23133,23134,23135,23136,23137,23138,23139,23140,23141,23142,23143,23144,23145,23146,23147,23148,23149,23150,23151,23398,23399,23400,23401,23402,23403,23404,23405,23406,23407,23408,23409,23410,23411,23412,23413,23414,23415,23416,23417,23418,23419,23420,23421,23422,23423,23424,23425,23560,23561,23562,23563,23564,23565,23566,23567,23568,23569,23570,23571,23572,23573,23574,23575,23576,23577,23578,23579,23580,23581,23582,23583,23584,23585,23586,23587,23833,23834,23835,23836,23837,23838,23839,23840,23841,23842,23843,23844,23845,23846,23847,23848,23849,23850,23851,23852,23853,23854,23855,23856,23857,23858,23859,23860,24050,24051,24052,24053,24054,24055,24056,24057,24058,24059,24060,24061,24062,24063,24064,24065,24066,24067,24068,24069,24070,24071,24072,24073,24074,24075,24076,24077,24212,24213,24214,24215,24216,24217,24218,24219,24220,24221,24222,24223,24224,24225,24226,24227,24228,24229,24230,24231,24232,24233,24234,24235,24236,24237,24238,24239',
       outframe='LSRK', gridder='mosaic', cycleniter=75, 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 16:25:9.227 -29.30.0.0000',
       scan=['357,360,364,367,370,373,376,380,383'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34897149.eb35035193.58161.27141606482_split.ms'],
       imagename='VLASS1.1.ql.T03t25.J162509-293000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00063322571957542091, imsize=[5760, 5760],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='22469,22470,22471,22472,22473,22474,22475,22476,22477,22478,22479,22480,22481,22482,22483,22484,22485,22486,22487,22488,22489,22490,22491,22492,22493,22494,22495,22496,22743,22744,22745,22746,22747,22748,22749,22750,22751,22752,22753,22754,22755,22756,22757,22758,22759,22760,22761,22762,22763,22764,22765,22766,22767,22768,22769,22770,22771,22962,22963,22964,22965,22966,22967,22968,22969,22970,22971,22972,22973,22974,22975,22976,22977,22978,22979,22980,22981,22982,22983,22984,22985,22986,22987,22988,22989,23124,23125,23126,23127,23128,23129,23130,23131,23132,23133,23134,23135,23136,23137,23138,23139,23140,23141,23142,23143,23144,23145,23146,23147,23148,23149,23150,23151,23398,23399,23400,23401,23402,23403,23404,23405,23406,23407,23408,23409,23410,23411,23412,23413,23414,23415,23416,23417,23418,23419,23420,23421,23422,23423,23424,23425,23560,23561,23562,23563,23564,23565,23566,23567,23568,23569,23570,23571,23572,23573,23574,23575,23576,23577,23578,23579,23580,23581,23582,23583,23584,23585,23586,23587,23833,23834,23835,23836,23837,23838,23839,23840,23841,23842,23843,23844,23845,23846,23847,23848,23849,23850,23851,23852,23853,23854,23855,23856,23857,23858,23859,23860,24050,24051,24052,24053,24054,24055,24056,24057,24058,24059,24060,24061,24062,24063,24064,24065,24066,24067,24068,24069,24070,24071,24072,24073,24074,24075,24076,24077,24212,24213,24214,24215,24216,24217,24218,24219,24220,24221,24222,24223,24224,24225,24226,24227,24228,24229,24230,24231,24232,24233,24234,24235,24236,24237,24238,24239',
       outframe='LSRK', gridder='mosaic', cycleniter=75, 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.T03t25.J162509-293000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T03t25.J162509-293000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T03t25.J162509-293000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#