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

# hif_editimlist(parameter_file='parameter.7513.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_T25t04.J041852+583000_P7513v1_2018_01_01T18_47_18.103/working/VLASS1.1.sb34459038.eb34579917.58037.11761364584_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T25t04.J041852+583000_P7513v1_2018_01_01T18_47_18.103/rawdata/VLASS1.1.sb34459038.eb34579917.58037.11761364584.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='23994,23995,23996,23997,23998,23999,24111,23693,23692,23691,23690,23697,23696,23695,23694,23699,23698,24554,23669,24883,24882,24881,24880,24885,24884,23149,23148,23147,23146,23145,24120,24121,24122,24123,23248,23249,24126,24127,23246,23247,24450,24451,24452,24453,22434,22435,22436,22437,22430,22431,22432,22433,24018,24019,24014,24015,24016,24017,24010,24011,24012,24013,22751,22750,22753,22752,22755,22754,24526,24527,24528,24529,24872,24873,24870,24871,24876,24877,24874,24875,24878,24879,22838,22839,23592,23593,23594,23595,23596,23597,22830,22831,22832,22833,22834,22835,22836,22837,23275,23274,23276,23271,23270,23273,23272,24119,24118,23688,23689,22845,22844,22847,22846,22841,22840,22843,24116,22409,22408,22407,22406,24025,24024,24021,24020,24023,24022,22424,23687,23158,23159,23150,23151,23152,23153,23154,23155,23156,23157,23680,24096,24097,24098,24099,23267,24429,24428,24423,24425,24424,24427,24426,23683,22748,22749,22742,22743,22740,22741,22746,22747,22744,22745,23575,24100,24101,22852,22853,24551,24550,24553,24552,24555,24104,24556,24105,24859,23569,23568,23590,23591,23161,23160,23163,23162,23165,23164,23167,23166,23169,23168,23266,23598,23264,23265,23262,23263,23260,23261,23268,23269,24108,24109,24102,24103,22854,22855,24106,24107,22850,22851,22418,22419,22416,22417,22414,22415,22412,22413,22410,22411,22733,22732,22731,22730,22737,22736,22735,22734,22739,22738,24438,24439,24432,24433,24430,24431,24436,24437,24434,24435,24857,24542,24543,24540,24541,24546,24547,24544,24545,22849,24548,24549,22848,24858,23684,23685,24124,23686,24856,24125,23671,23670,23673,23672,23675,23674,23677,23676,23679,23678,23681,23682,23578,23579,23576,23577,23574,22842,23572,23573,23570,23571,24110,23172,23173,23170,23171,23174,23175,24113,23253,23252,23251,23250,23257,23256,23255,23254,23259,23258,24112,24855,24443,24442,24441,24440,24447,24446,24445,24444,24449,24448,24115,22829,22828,22429,22428,22423,22422,22421,22420,22427,22426,22425,24114,24003,24002,24001,24000,24007,24006,24005,24004,24009,24008,24117,22728,22729,22725,22726,22727,24539,24538,24533,24532,24531,24530,24537,24536,24535,24534,24861,24860,24863,24862,24865,24864,24867,24866,24869,24868,23587,23586,23585,23584,23583,23582,23581,23580,22827,22826,22825,22824,23589,23588',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 04:18:52.248 +58.30.0.0000',
       scan=['386,389,392,405,395,399,408,402,373,376,380,383'],
       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_T25t04.J041852+583000_P7513v1_2018_01_01T18_47_18.103/working/VLASS1.1.sb34459038.eb34579917.58037.11761364584_split.ms'],
       imagename='VLASS1.1.ql.T25t04.J041852+583000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='22406,22407,22408,22409,22410,22411,22412,22413,22414,22415,22416,22417,22418,22419,22420,22421,22422,22423,22424,22425,22426,22427,22428,22429,22430,22431,22432,22433,22434,22435,22436,22437,22725,22726,22727,22728,22729,22730,22731,22732,22733,22734,22735,22736,22737,22738,22739,22740,22741,22742,22743,22744,22745,22746,22747,22748,22749,22750,22751,22752,22753,22754,22755,22824,22825,22826,22827,22828,22829,22830,22831,22832,22833,22834,22835,22836,22837,22838,22839,22840,22841,22842,22843,22844,22845,22846,22847,22848,22849,22850,22851,22852,22853,22854,22855,23145,23146,23147,23148,23149,23150,23151,23152,23153,23154,23155,23156,23157,23158,23159,23160,23161,23162,23163,23164,23165,23166,23167,23168,23169,23170,23171,23172,23173,23174,23175,23246,23247,23248,23249,23250,23251,23252,23253,23254,23255,23256,23257,23258,23259,23260,23261,23262,23263,23264,23265,23266,23267,23268,23269,23270,23271,23272,23273,23274,23275,23276,23568,23569,23570,23571,23572,23573,23574,23575,23576,23577,23578,23579,23580,23581,23582,23583,23584,23585,23586,23587,23588,23589,23590,23591,23592,23593,23594,23595,23596,23597,23598,23669,23670,23671,23672,23673,23674,23675,23676,23677,23678,23679,23680,23681,23682,23683,23684,23685,23686,23687,23688,23689,23690,23691,23692,23693,23694,23695,23696,23697,23698,23699,23994,23995,23996,23997,23998,23999,24000,24001,24002,24003,24004,24005,24006,24007,24008,24009,24010,24011,24012,24013,24014,24015,24016,24017,24018,24019,24020,24021,24022,24023,24024,24025,24096,24097,24098,24099,24100,24101,24102,24103,24104,24105,24106,24107,24108,24109,24110,24111,24112,24113,24114,24115,24116,24117,24118,24119,24120,24121,24122,24123,24124,24125,24126,24127,24423,24424,24425,24426,24427,24428,24429,24430,24431,24432,24433,24434,24435,24436,24437,24438,24439,24440,24441,24442,24443,24444,24445,24446,24447,24448,24449,24450,24451,24452,24453,24526,24527,24528,24529,24530,24531,24532,24533,24534,24535,24536,24537,24538,24539,24540,24541,24542,24543,24544,24545,24546,24547,24548,24549,24550,24551,24552,24553,24554,24555,24556,24855,24856,24857,24858,24859,24860,24861,24862,24863,24864,24865,24866,24867,24868,24869,24870,24871,24872,24873,24874,24875,24876,24877,24878,24879,24880,24881,24882,24883,24884,24885',
       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:18:52.248 +58.30.0.0000',
       scan=['386,389,392,405,395,399,408,402,373,376,380,383'],
       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_T25t04.J041852+583000_P7513v1_2018_01_01T18_47_18.103/working/VLASS1.1.sb34459038.eb34579917.58037.11761364584_split.ms'],
       imagename='VLASS1.1.ql.T25t04.J041852+583000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.0004941820257041399, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='22406,22407,22408,22409,22410,22411,22412,22413,22414,22415,22416,22417,22418,22419,22420,22421,22422,22423,22424,22425,22426,22427,22428,22429,22430,22431,22432,22433,22434,22435,22436,22437,22725,22726,22727,22728,22729,22730,22731,22732,22733,22734,22735,22736,22737,22738,22739,22740,22741,22742,22743,22744,22745,22746,22747,22748,22749,22750,22751,22752,22753,22754,22755,22824,22825,22826,22827,22828,22829,22830,22831,22832,22833,22834,22835,22836,22837,22838,22839,22840,22841,22842,22843,22844,22845,22846,22847,22848,22849,22850,22851,22852,22853,22854,22855,23145,23146,23147,23148,23149,23150,23151,23152,23153,23154,23155,23156,23157,23158,23159,23160,23161,23162,23163,23164,23165,23166,23167,23168,23169,23170,23171,23172,23173,23174,23175,23246,23247,23248,23249,23250,23251,23252,23253,23254,23255,23256,23257,23258,23259,23260,23261,23262,23263,23264,23265,23266,23267,23268,23269,23270,23271,23272,23273,23274,23275,23276,23568,23569,23570,23571,23572,23573,23574,23575,23576,23577,23578,23579,23580,23581,23582,23583,23584,23585,23586,23587,23588,23589,23590,23591,23592,23593,23594,23595,23596,23597,23598,23669,23670,23671,23672,23673,23674,23675,23676,23677,23678,23679,23680,23681,23682,23683,23684,23685,23686,23687,23688,23689,23690,23691,23692,23693,23694,23695,23696,23697,23698,23699,23994,23995,23996,23997,23998,23999,24000,24001,24002,24003,24004,24005,24006,24007,24008,24009,24010,24011,24012,24013,24014,24015,24016,24017,24018,24019,24020,24021,24022,24023,24024,24025,24096,24097,24098,24099,24100,24101,24102,24103,24104,24105,24106,24107,24108,24109,24110,24111,24112,24113,24114,24115,24116,24117,24118,24119,24120,24121,24122,24123,24124,24125,24126,24127,24423,24424,24425,24426,24427,24428,24429,24430,24431,24432,24433,24434,24435,24436,24437,24438,24439,24440,24441,24442,24443,24444,24445,24446,24447,24448,24449,24450,24451,24452,24453,24526,24527,24528,24529,24530,24531,24532,24533,24534,24535,24536,24537,24538,24539,24540,24541,24542,24543,24544,24545,24546,24547,24548,24549,24550,24551,24552,24553,24554,24555,24556,24855,24856,24857,24858,24859,24860,24861,24862,24863,24864,24865,24866,24867,24868,24869,24870,24871,24872,24873,24874,24875,24876,24877,24878,24879,24880,24881,24882,24883,24884,24885',
       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.T25t04.J041852+583000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T25t04.J041852+583000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T25t04.J041852+583000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#