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

# hif_editimlist(parameter_file='parameter.15694.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_T25t02.J020023+573000_P15694v1_2017_11_18T00_58_34.952/working/VLASS1.1.sb34614479.eb34620102.58046.99344417824_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T25t02.J020023+573000_P15694v1_2017_11_18T00_58_34.952/rawdata/VLASS1.1.sb34614479.eb34620102.58046.99344417824.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='26124,24168,24169,24164,24165,24166,24167,24160,24161,24162,24163,25899,25699,25698,25691,25690,25693,25692,25695,25694,25697,25696,24830,25888,25889,26126,24379,24378,26127,24375,24374,24377,24376,24371,24370,24373,24372,26125,26122,26123,24812,26121,25233,25237,25236,25235,25234,25239,25238,25905,25904,25907,25906,24591,24590,25903,24592,24599,24598,24155,24154,24157,24156,24151,24153,24152,24159,24158,25019,25018,25017,25016,25015,25014,25475,25474,25477,25476,25471,25470,25473,25472,25479,25478,24799,26360,25688,25689,25682,25683,25680,25681,25686,25687,25684,25685,24809,24808,24803,24802,24801,24800,24807,24806,24805,24804,26140,26141,24368,24369,26135,26134,26137,26136,26131,26130,26133,26132,26139,26138,25916,25917,25914,25915,25912,25913,25910,25911,24588,25919,26351,26350,26353,26352,26355,26354,26357,26356,26359,26358,25480,25260,25261,25262,25263,25264,25466,25467,25464,25465,25462,25463,25460,25461,25468,25469,24397,24396,24395,24394,24393,24392,24391,24390,24398,24586,24587,24584,24585,24582,24583,25043,25040,25041,25042,24581,25044,24605,24604,24607,24606,24601,24600,24603,24602,24589,24609,24608,25039,25038,25035,25034,25037,25036,25031,25030,25033,25032,25259,25258,25251,25250,25253,25252,25255,25254,25257,25256,24814,24815,24816,24817,24810,24811,26120,24813,24818,24819,26128,26129,25679,25678,25673,25672,25671,25677,25676,25675,25674,26342,26343,26340,26341,26346,26347,26344,26345,26348,26349,24173,24172,24171,24170,24177,24176,24175,24174,24179,24178,24597,24829,24828,24825,24824,24827,24826,24821,24820,24823,24822,25453,25452,25451,25450,25457,25456,25455,25454,25459,25458,24380,24381,24382,24383,24384,24385,24386,24387,24388,24389,25891,25890,25897,25896,25895,25894,25893,25892,24182,24180,24181,25701,25700,24595,26339,26338,24594,26333,26332,26331,26330,26337,26336,26335,26334,24596,24612,24610,24611,25901,25900,24593,25902,25028,25029,25026,25027,25024,25025,25022,25023,25020,25021,25909,25908,25248,25249,25242,25243,25240,25241,25246,25247,25244,25245,25918,26113,26112,26111,26117,26116,26115,26114,26119,26118,25898',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 02:00:23.624 +57.30.0.0000',
       scan=['418,421,424,427,430,399,402,405,408,411,414'],
       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_T25t02.J020023+573000_P15694v1_2017_11_18T00_58_34.952/working/VLASS1.1.sb34614479.eb34620102.58046.99344417824_split.ms'],
       imagename='VLASS1.1.ql.T25t02.J020023+573000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='24151,24152,24153,24154,24155,24156,24157,24158,24159,24160,24161,24162,24163,24164,24165,24166,24167,24168,24169,24170,24171,24172,24173,24174,24175,24176,24177,24178,24179,24180,24181,24182,24368,24369,24370,24371,24372,24373,24374,24375,24376,24377,24378,24379,24380,24381,24382,24383,24384,24385,24386,24387,24388,24389,24390,24391,24392,24393,24394,24395,24396,24397,24398,24581,24582,24583,24584,24585,24586,24587,24588,24589,24590,24591,24592,24593,24594,24595,24596,24597,24598,24599,24600,24601,24602,24603,24604,24605,24606,24607,24608,24609,24610,24611,24612,24799,24800,24801,24802,24803,24804,24805,24806,24807,24808,24809,24810,24811,24812,24813,24814,24815,24816,24817,24818,24819,24820,24821,24822,24823,24824,24825,24826,24827,24828,24829,24830,25014,25015,25016,25017,25018,25019,25020,25021,25022,25023,25024,25025,25026,25027,25028,25029,25030,25031,25032,25033,25034,25035,25036,25037,25038,25039,25040,25041,25042,25043,25044,25233,25234,25235,25236,25237,25238,25239,25240,25241,25242,25243,25244,25245,25246,25247,25248,25249,25250,25251,25252,25253,25254,25255,25256,25257,25258,25259,25260,25261,25262,25263,25264,25450,25451,25452,25453,25454,25455,25456,25457,25458,25459,25460,25461,25462,25463,25464,25465,25466,25467,25468,25469,25470,25471,25472,25473,25474,25475,25476,25477,25478,25479,25480,25671,25672,25673,25674,25675,25676,25677,25678,25679,25680,25681,25682,25683,25684,25685,25686,25687,25688,25689,25690,25691,25692,25693,25694,25695,25696,25697,25698,25699,25700,25701,25888,25889,25890,25891,25892,25893,25894,25895,25896,25897,25898,25899,25900,25901,25902,25903,25904,25905,25906,25907,25908,25909,25910,25911,25912,25913,25914,25915,25916,25917,25918,25919,26111,26112,26113,26114,26115,26116,26117,26118,26119,26120,26121,26122,26123,26124,26125,26126,26127,26128,26129,26130,26131,26132,26133,26134,26135,26136,26137,26138,26139,26140,26141,26330,26331,26332,26333,26334,26335,26336,26337,26338,26339,26340,26341,26342,26343,26344,26345,26346,26347,26348,26349,26350,26351,26352,26353,26354,26355,26356,26357,26358,26359,26360',
       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 02:00:23.624 +57.30.0.0000',
       scan=['418,421,424,427,430,399,402,405,408,411,414'],
       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_T25t02.J020023+573000_P15694v1_2017_11_18T00_58_34.952/working/VLASS1.1.sb34614479.eb34620102.58046.99344417824_split.ms'],
       imagename='VLASS1.1.ql.T25t02.J020023+573000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00052179050192353318, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='24151,24152,24153,24154,24155,24156,24157,24158,24159,24160,24161,24162,24163,24164,24165,24166,24167,24168,24169,24170,24171,24172,24173,24174,24175,24176,24177,24178,24179,24180,24181,24182,24368,24369,24370,24371,24372,24373,24374,24375,24376,24377,24378,24379,24380,24381,24382,24383,24384,24385,24386,24387,24388,24389,24390,24391,24392,24393,24394,24395,24396,24397,24398,24581,24582,24583,24584,24585,24586,24587,24588,24589,24590,24591,24592,24593,24594,24595,24596,24597,24598,24599,24600,24601,24602,24603,24604,24605,24606,24607,24608,24609,24610,24611,24612,24799,24800,24801,24802,24803,24804,24805,24806,24807,24808,24809,24810,24811,24812,24813,24814,24815,24816,24817,24818,24819,24820,24821,24822,24823,24824,24825,24826,24827,24828,24829,24830,25014,25015,25016,25017,25018,25019,25020,25021,25022,25023,25024,25025,25026,25027,25028,25029,25030,25031,25032,25033,25034,25035,25036,25037,25038,25039,25040,25041,25042,25043,25044,25233,25234,25235,25236,25237,25238,25239,25240,25241,25242,25243,25244,25245,25246,25247,25248,25249,25250,25251,25252,25253,25254,25255,25256,25257,25258,25259,25260,25261,25262,25263,25264,25450,25451,25452,25453,25454,25455,25456,25457,25458,25459,25460,25461,25462,25463,25464,25465,25466,25467,25468,25469,25470,25471,25472,25473,25474,25475,25476,25477,25478,25479,25480,25671,25672,25673,25674,25675,25676,25677,25678,25679,25680,25681,25682,25683,25684,25685,25686,25687,25688,25689,25690,25691,25692,25693,25694,25695,25696,25697,25698,25699,25700,25701,25888,25889,25890,25891,25892,25893,25894,25895,25896,25897,25898,25899,25900,25901,25902,25903,25904,25905,25906,25907,25908,25909,25910,25911,25912,25913,25914,25915,25916,25917,25918,25919,26111,26112,26113,26114,26115,26116,26117,26118,26119,26120,26121,26122,26123,26124,26125,26126,26127,26128,26129,26130,26131,26132,26133,26134,26135,26136,26137,26138,26139,26140,26141,26330,26331,26332,26333,26334,26335,26336,26337,26338,26339,26340,26341,26342,26343,26344,26345,26346,26347,26348,26349,26350,26351,26352,26353,26354,26355,26356,26357,26358,26359,26360',
       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.T25t02.J020023+573000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T25t02.J020023+573000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T25t02.J020023+573000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#