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.sb34504025.eb34557675.58027.36083796296.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='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T11t13.J083800+003000_P10398v1_2018_03_21T19_29_52.899/working/VLASS1.1.sb34504025.eb34557675.58027.36083796296_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T11t13.J083800+003000_P10398v1_2018_03_21T19_29_52.899/rawdata/VLASS1.1.sb34504025.eb34557675.58027.36083796296.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='17478,17479,18259,18280,17472,17473,17474,17475,17476,17477,16684,16685,16680,16681,16682,16683,18712,18288,18299,18289,18709,16671,16670,16673,16672,16675,16674,16677,16676,16679,16678,18679,18678,18681,18298,18708,17485,17484,17487,17486,17481,17480,17483,17482,17856,17857,17854,17855,17489,17488,17850,17851,17081,17080,17083,17082,17085,17084,17087,17086,17089,17088,18710,18290,16662,16663,16660,16661,16666,16667,16664,16665,18664,17852,16668,16669,18301,18662,17853,18672,17900,17901,17902,17903,18302,17496,17497,17494,17495,17492,17493,17490,17491,17863,17862,17861,17860,17867,17866,17865,17864,18714,17090,17091,18715,18705,18257,17869,17038,17039,18256,18295,18711,17874,17875,18692,17870,17871,17872,17873,18308,17878,17879,18707,18696,18294,18683,18699,18279,18278,18277,18276,18275,18274,18273,18272,18271,18270,18703,18666,18702,18297,18701,17045,17044,17047,17046,17041,17040,17043,17042,18668,17049,17048,18669,18691,18268,18693,18309,18695,18694,18697,18269,18296,18698,18260,17889,17888,18281,18293,18665,17881,17880,17883,17882,17885,17884,17887,17886,18306,18667,17445,17444,17447,17446,17449,17448,17056,17057,17054,17055,17052,17053,17050,17051,18291,18663,18262,17058,17059,18682,18267,18680,18263,18686,18687,18684,18685,18264,18670,18700,18265,17898,17899,18266,17892,17893,17890,17891,17896,17897,17894,17895,17452,17453,17450,17451,17456,17457,17454,17455,17458,17459,18671,17063,17062,17061,17060,17067,17066,18690,18677,17069,17068,18676,18292,18675,18706,18674,17858,18713,18307,18261,18305,18304,18303,17859,17469,17468,17467,17466,17465,17464,17463,17462,17461,17460,17868,17078,17079,17074,17075,17076,17077,17070,17071,17072,17073,18286,18704,18287,18284,18285,18300,18258,18673',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 08:38:0.041 +00.30.0.0000',
       scan=['353,343,321,337,296,324,330,299,302,305,308,334,311,327,340,346,315,349,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_T11t13.J083800+003000_P10398v1_2018_03_21T19_29_52.899/working/VLASS1.1.sb34504025.eb34557675.58027.36083796296_split.ms'],
       imagename='VLASS1.1.ql.T11t13.J083800+003000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='16660,16661,16662,16663,16664,16665,16666,16667,16668,16669,16670,16671,16672,16673,16674,16675,16676,16677,16678,16679,16680,16681,16682,16683,16684,16685,17038,17039,17040,17041,17042,17043,17044,17045,17046,17047,17048,17049,17050,17051,17052,17053,17054,17055,17056,17057,17058,17059,17060,17061,17062,17063,17066,17067,17068,17069,17070,17071,17072,17073,17074,17075,17076,17077,17078,17079,17080,17081,17082,17083,17084,17085,17086,17087,17088,17089,17090,17091,17444,17445,17446,17447,17448,17449,17450,17451,17452,17453,17454,17455,17456,17457,17458,17459,17460,17461,17462,17463,17464,17465,17466,17467,17468,17469,17472,17473,17474,17475,17476,17477,17478,17479,17480,17481,17482,17483,17484,17485,17486,17487,17488,17489,17490,17491,17492,17493,17494,17495,17496,17497,17850,17851,17852,17853,17854,17855,17856,17857,17858,17859,17860,17861,17862,17863,17864,17865,17866,17867,17868,17869,17870,17871,17872,17873,17874,17875,17878,17879,17880,17881,17882,17883,17884,17885,17886,17887,17888,17889,17890,17891,17892,17893,17894,17895,17896,17897,17898,17899,17900,17901,17902,17903,18256,18257,18258,18259,18260,18261,18262,18263,18264,18265,18266,18267,18268,18269,18270,18271,18272,18273,18274,18275,18276,18277,18278,18279,18280,18281,18284,18285,18286,18287,18288,18289,18290,18291,18292,18293,18294,18295,18296,18297,18298,18299,18300,18301,18302,18303,18304,18305,18306,18307,18308,18309,18662,18663,18664,18665,18666,18667,18668,18669,18670,18671,18672,18673,18674,18675,18676,18677,18678,18679,18680,18681,18682,18683,18684,18685,18686,18687,18690,18691,18692,18693,18694,18695,18696,18697,18698,18699,18700,18701,18702,18703,18704,18705,18706,18707,18708,18709,18710,18711,18712,18713,18714,18715',
       outframe='LSRK', gridder='mosaic', cycleniter=200, 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 08:38:0.041 +00.30.0.0000',
       scan=['353,343,321,337,296,324,330,299,302,305,308,334,311,327,340,346,315,349,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_T11t13.J083800+003000_P10398v1_2018_03_21T19_29_52.899/working/VLASS1.1.sb34504025.eb34557675.58027.36083796296_split.ms'],
       imagename='VLASS1.1.ql.T11t13.J083800+003000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00052949032887700015, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='16660,16661,16662,16663,16664,16665,16666,16667,16668,16669,16670,16671,16672,16673,16674,16675,16676,16677,16678,16679,16680,16681,16682,16683,16684,16685,17038,17039,17040,17041,17042,17043,17044,17045,17046,17047,17048,17049,17050,17051,17052,17053,17054,17055,17056,17057,17058,17059,17060,17061,17062,17063,17066,17067,17068,17069,17070,17071,17072,17073,17074,17075,17076,17077,17078,17079,17080,17081,17082,17083,17084,17085,17086,17087,17088,17089,17090,17091,17444,17445,17446,17447,17448,17449,17450,17451,17452,17453,17454,17455,17456,17457,17458,17459,17460,17461,17462,17463,17464,17465,17466,17467,17468,17469,17472,17473,17474,17475,17476,17477,17478,17479,17480,17481,17482,17483,17484,17485,17486,17487,17488,17489,17490,17491,17492,17493,17494,17495,17496,17497,17850,17851,17852,17853,17854,17855,17856,17857,17858,17859,17860,17861,17862,17863,17864,17865,17866,17867,17868,17869,17870,17871,17872,17873,17874,17875,17878,17879,17880,17881,17882,17883,17884,17885,17886,17887,17888,17889,17890,17891,17892,17893,17894,17895,17896,17897,17898,17899,17900,17901,17902,17903,18256,18257,18258,18259,18260,18261,18262,18263,18264,18265,18266,18267,18268,18269,18270,18271,18272,18273,18274,18275,18276,18277,18278,18279,18280,18281,18284,18285,18286,18287,18288,18289,18290,18291,18292,18293,18294,18295,18296,18297,18298,18299,18300,18301,18302,18303,18304,18305,18306,18307,18308,18309,18662,18663,18664,18665,18666,18667,18668,18669,18670,18671,18672,18673,18674,18675,18676,18677,18678,18679,18680,18681,18682,18683,18684,18685,18686,18687,18690,18691,18692,18693,18694,18695,18696,18697,18698,18699,18700,18701,18702,18703,18704,18705,18706,18707,18708,18709,18710,18711,18712,18713,18714,18715',
       outframe='LSRK', gridder='mosaic', cycleniter=200, 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.T11t13.J083800+003000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T11t13.J083800+003000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T11t13.J083800+003000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#