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

# hif_editimlist(parameter_file='parameter.17139.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_T29t03.J075716+753000_P17139v1_2018_01_18T17_04_39.994/working/VLASS1.1.sb34667861.eb34684334.58067.0700902199_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T29t03.J075716+753000_P17139v1_2018_01_18T17_04_39.994/rawdata/VLASS1.1.sb34667861.eb34684334.58067.0700902199.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='20447,20446,20445,20444,20443,20442,20441,20440,21901,21900,21903,21902,21488,21486,20449,20448,21080,21081,21082,21083,21084,21085,21086,22677,21480,22289,21481,22288,21908,22672,22673,23079,23078,22693,21912,21905,21910,21911,23077,23076,23075,21489,21079,21078,22678,22679,21075,21074,21077,21076,21071,21070,21073,21072,21479,21478,22689,21471,21470,21473,21472,21475,21474,21477,21476,22276,22277,22278,22279,22688,22676,23065,22705,22711,22710,21068,21069,21064,21065,21062,21063,21060,21061,21468,21469,21501,21500,21502,21463,21466,21467,21464,21465,22674,22675,23067,20436,20437,20434,20435,20432,20433,21484,20431,21053,21052,21051,21050,21057,21056,21055,21054,22706,23071,22707,21059,22311,22310,22313,21058,21485,22314,23062,23063,22702,23068,23069,22703,23070,21047,22700,22701,21048,21049,21493,21492,23080,21491,22708,22302,21490,21487,22303,22300,22301,22306,20438,22304,22305,22308,20439,22709,22307,22698,22283,22282,22281,22280,20649,20648,22285,22284,20645,20647,20646,22699,22309,22686,21904,22683,23074,21906,22682,22681,21907,22680,22687,22294,22685,22684,22298,22299,20658,20659,20656,20657,20654,20655,20652,20653,20650,20651,23064,22297,21899,21898,21897,21896,21895,21894,23066,21891,21890,23073,22312,22292,23072,22704,22315,22293,22290,21875,21874,21877,21876,21873,22291,21879,21878,22690,21499,21498,22694,22695,22696,22697,20663,20662,20661,20660,21497,21496,21495,21494,21888,21889,21909,21880,21881,21882,21883,21884,21885,21886,21887,22287,22286',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 07:57:16.444 +75.30.0.0000',
       scan=['385,388,401,392,395,398,369,404,373,407,376,379,382'],
       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_T29t03.J075716+753000_P17139v1_2018_01_18T17_04_39.994/working/VLASS1.1.sb34667861.eb34684334.58067.0700902199_split.ms'],
       imagename='VLASS1.1.ql.T29t03.J075716+753000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='20431,20432,20433,20434,20435,20436,20437,20438,20439,20440,20441,20442,20443,20444,20445,20446,20447,20448,20449,20645,20646,20647,20648,20649,20650,20651,20652,20653,20654,20655,20656,20657,20658,20659,20660,20661,20662,20663,21047,21048,21049,21050,21051,21052,21053,21054,21055,21056,21057,21058,21059,21060,21061,21062,21063,21064,21065,21068,21069,21070,21071,21072,21073,21074,21075,21076,21077,21078,21079,21080,21081,21082,21083,21084,21085,21086,21463,21464,21465,21466,21467,21468,21469,21470,21471,21472,21473,21474,21475,21476,21477,21478,21479,21480,21481,21484,21485,21486,21487,21488,21489,21490,21491,21492,21493,21494,21495,21496,21497,21498,21499,21500,21501,21502,21873,21874,21875,21876,21877,21878,21879,21880,21881,21882,21883,21884,21885,21886,21887,21888,21889,21890,21891,21894,21895,21896,21897,21898,21899,21900,21901,21902,21903,21904,21905,21906,21907,21908,21909,21910,21911,21912,22276,22277,22278,22279,22280,22281,22282,22283,22284,22285,22286,22287,22288,22289,22290,22291,22292,22293,22294,22297,22298,22299,22300,22301,22302,22303,22304,22305,22306,22307,22308,22309,22310,22311,22312,22313,22314,22315,22672,22673,22674,22675,22676,22677,22678,22679,22680,22681,22682,22683,22684,22685,22686,22687,22688,22689,22690,22693,22694,22695,22696,22697,22698,22699,22700,22701,22702,22703,22704,22705,22706,22707,22708,22709,22710,22711,23062,23063,23064,23065,23066,23067,23068,23069,23070,23071,23072,23073,23074,23075,23076,23077,23078,23079,23080',
       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 07:57:16.444 +75.30.0.0000',
       scan=['385,388,401,392,395,398,369,404,373,407,376,379,382'],
       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_T29t03.J075716+753000_P17139v1_2018_01_18T17_04_39.994/working/VLASS1.1.sb34667861.eb34684334.58067.0700902199_split.ms'],
       imagename='VLASS1.1.ql.T29t03.J075716+753000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00050919232424712383, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='20431,20432,20433,20434,20435,20436,20437,20438,20439,20440,20441,20442,20443,20444,20445,20446,20447,20448,20449,20645,20646,20647,20648,20649,20650,20651,20652,20653,20654,20655,20656,20657,20658,20659,20660,20661,20662,20663,21047,21048,21049,21050,21051,21052,21053,21054,21055,21056,21057,21058,21059,21060,21061,21062,21063,21064,21065,21068,21069,21070,21071,21072,21073,21074,21075,21076,21077,21078,21079,21080,21081,21082,21083,21084,21085,21086,21463,21464,21465,21466,21467,21468,21469,21470,21471,21472,21473,21474,21475,21476,21477,21478,21479,21480,21481,21484,21485,21486,21487,21488,21489,21490,21491,21492,21493,21494,21495,21496,21497,21498,21499,21500,21501,21502,21873,21874,21875,21876,21877,21878,21879,21880,21881,21882,21883,21884,21885,21886,21887,21888,21889,21890,21891,21894,21895,21896,21897,21898,21899,21900,21901,21902,21903,21904,21905,21906,21907,21908,21909,21910,21911,21912,22276,22277,22278,22279,22280,22281,22282,22283,22284,22285,22286,22287,22288,22289,22290,22291,22292,22293,22294,22297,22298,22299,22300,22301,22302,22303,22304,22305,22306,22307,22308,22309,22310,22311,22312,22313,22314,22315,22672,22673,22674,22675,22676,22677,22678,22679,22680,22681,22682,22683,22684,22685,22686,22687,22688,22689,22690,22693,22694,22695,22696,22697,22698,22699,22700,22701,22702,22703,22704,22705,22706,22707,22708,22709,22710,22711,23062,23063,23064,23065,23066,23067,23068,23069,23070,23071,23072,23073,23074,23075,23076,23077,23078,23079,23080',
       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.T29t03.J075716+753000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T29t03.J075716+753000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T29t03.J075716+753000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#