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

# hif_editimlist(parameter_file='parameter.16390.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_T27t06.J080443+653000_P16390v1_2018_01_07T15_11_05.607/working/VLASS1.1.sb34640437.eb34680548.58066.26112458333_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T27t06.J080443+653000_P16390v1_2018_01_07T15_11_05.607/rawdata/VLASS1.1.sb34640437.eb34680548.58066.26112458333.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='11829,11828,11825,11824,11827,11826,11821,11820,11823,11822,11380,11381,11382,11383,11384,11385,11386,11387,11388,11389,10581,10580,10583,10582,10585,10584,10587,10586,10589,10588,12594,12595,12596,12597,12598,12599,11003,11002,11001,11007,11006,11005,11004,11009,11008,4101,4100,4103,4102,4104,5704,5700,5701,5702,5703,12235,12234,12236,12231,12230,12233,12232,11413,12198,12199,12190,12191,12192,12193,12194,12195,12196,12197,11407,11406,11405,11404,11403,11402,11401,6108,6105,6104,6103,6102,11409,11408,6109,6101,6100,11379,11378,11375,11374,11377,11376,5692,10593,10590,10591,5696,10597,10594,5695,10598,5699,4488,4489,10974,10975,10976,10977,10978,10979,11018,11019,11014,11015,11016,11017,11010,11011,11012,11013,10628,10620,10621,10622,10623,10624,10625,10626,10627,12226,12227,12224,12225,12222,12223,12220,12221,12228,12229,10578,10579,12619,12618,12617,12616,12615,12614,12613,12612,12611,12610,11410,11411,11799,11798,11414,11415,11416,11417,11793,11792,11791,11790,11797,11796,11795,11794,4888,4889,4885,4886,4887,11025,11024,11021,11020,11023,11022,10611,10610,10613,10612,10615,10614,10617,10616,10619,10618,11809,11808,11801,11800,11807,11806,11805,11804,4497,4496,4495,4494,4491,4490,4499,4498,10989,10988,10985,10984,10987,10986,10981,10980,10983,10982,11421,11420,11423,11422,11425,11424,5298,5299,5296,5297,5294,5295,5290,5291,10577,12213,12212,12211,12217,12216,12215,12214,12219,12218,12608,12609,12600,12601,12602,12603,12604,12605,12606,12607,11397,11396,11395,11394,11393,11392,11391,11390,11398,10601,11784,11785,11786,11787,11780,11781,11782,11783,11788,11789,11779,11778,11777,11776,4898,4893,4890,4897,4896,4895,4894,10592,5693,4486,10600,5691,10606,10607,10604,10605,10596,10608,10609,5694,10595,10599,10998,10996,10997,10994,10995,10992,10993,10990,10991,11814,11815,11816,11817,11810,11811,11812,11813,11818,11819,11412,4487,12189,12188,12187,12186,12185,12184,12183,6110,6111,6112,6113,11418,11419,5289,5288,4099,5287,5286,12208,12204,12205,12206,12207,12200,12201,12202,12203',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 08:04:43.944 +65.30.0.0000',
       scan=['184,187,190,193,196,200,203,206,79,209,82,212,86,215,89,92,95,98,101,105,108,111'],
       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_T27t06.J080443+653000_P16390v1_2018_01_07T15_11_05.607/working/VLASS1.1.sb34640437.eb34680548.58066.26112458333_split.ms'],
       imagename='VLASS1.1.ql.T27t06.J080443+653000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='4099,4100,4101,4102,4103,4104,4486,4487,4488,4489,4490,4491,4494,4495,4496,4497,4498,4499,4885,4886,4887,4888,4889,4890,4893,4894,4895,4896,4897,4898,5286,5287,5288,5289,5290,5291,5294,5295,5296,5297,5298,5299,5691,5692,5693,5694,5695,5696,5699,5700,5701,5702,5703,5704,6100,6101,6102,6103,6104,6105,6108,6109,6110,6111,6112,6113,10577,10578,10579,10580,10581,10582,10583,10584,10585,10586,10587,10588,10589,10590,10591,10592,10593,10594,10595,10596,10597,10598,10599,10600,10601,10604,10605,10606,10607,10608,10609,10610,10611,10612,10613,10614,10615,10616,10617,10618,10619,10620,10621,10622,10623,10624,10625,10626,10627,10628,10974,10975,10976,10977,10978,10979,10980,10981,10982,10983,10984,10985,10986,10987,10988,10989,10990,10991,10992,10993,10994,10995,10996,10997,10998,11001,11002,11003,11004,11005,11006,11007,11008,11009,11010,11011,11012,11013,11014,11015,11016,11017,11018,11019,11020,11021,11022,11023,11024,11025,11374,11375,11376,11377,11378,11379,11380,11381,11382,11383,11384,11385,11386,11387,11388,11389,11390,11391,11392,11393,11394,11395,11396,11397,11398,11401,11402,11403,11404,11405,11406,11407,11408,11409,11410,11411,11412,11413,11414,11415,11416,11417,11418,11419,11420,11421,11422,11423,11424,11425,11776,11777,11778,11779,11780,11781,11782,11783,11784,11785,11786,11787,11788,11789,11790,11791,11792,11793,11794,11795,11796,11797,11798,11799,11800,11801,11804,11805,11806,11807,11808,11809,11810,11811,11812,11813,11814,11815,11816,11817,11818,11819,11820,11821,11822,11823,11824,11825,11826,11827,11828,11829,12183,12184,12185,12186,12187,12188,12189,12190,12191,12192,12193,12194,12195,12196,12197,12198,12199,12200,12201,12202,12203,12204,12205,12206,12207,12208,12211,12212,12213,12214,12215,12216,12217,12218,12219,12220,12221,12222,12223,12224,12225,12226,12227,12228,12229,12230,12231,12232,12233,12234,12235,12236,12594,12595,12596,12597,12598,12599,12600,12601,12602,12603,12604,12605,12606,12607,12608,12609,12610,12611,12612,12613,12614,12615,12616,12617,12618,12619',
       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 08:04:43.944 +65.30.0.0000',
       scan=['184,187,190,193,196,200,203,206,79,209,82,212,86,215,89,92,95,98,101,105,108,111'],
       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_T27t06.J080443+653000_P16390v1_2018_01_07T15_11_05.607/working/VLASS1.1.sb34640437.eb34680548.58066.26112458333_split.ms'],
       imagename='VLASS1.1.ql.T27t06.J080443+653000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00048425323064730033, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='4099,4100,4101,4102,4103,4104,4486,4487,4488,4489,4490,4491,4494,4495,4496,4497,4498,4499,4885,4886,4887,4888,4889,4890,4893,4894,4895,4896,4897,4898,5286,5287,5288,5289,5290,5291,5294,5295,5296,5297,5298,5299,5691,5692,5693,5694,5695,5696,5699,5700,5701,5702,5703,5704,6100,6101,6102,6103,6104,6105,6108,6109,6110,6111,6112,6113,10577,10578,10579,10580,10581,10582,10583,10584,10585,10586,10587,10588,10589,10590,10591,10592,10593,10594,10595,10596,10597,10598,10599,10600,10601,10604,10605,10606,10607,10608,10609,10610,10611,10612,10613,10614,10615,10616,10617,10618,10619,10620,10621,10622,10623,10624,10625,10626,10627,10628,10974,10975,10976,10977,10978,10979,10980,10981,10982,10983,10984,10985,10986,10987,10988,10989,10990,10991,10992,10993,10994,10995,10996,10997,10998,11001,11002,11003,11004,11005,11006,11007,11008,11009,11010,11011,11012,11013,11014,11015,11016,11017,11018,11019,11020,11021,11022,11023,11024,11025,11374,11375,11376,11377,11378,11379,11380,11381,11382,11383,11384,11385,11386,11387,11388,11389,11390,11391,11392,11393,11394,11395,11396,11397,11398,11401,11402,11403,11404,11405,11406,11407,11408,11409,11410,11411,11412,11413,11414,11415,11416,11417,11418,11419,11420,11421,11422,11423,11424,11425,11776,11777,11778,11779,11780,11781,11782,11783,11784,11785,11786,11787,11788,11789,11790,11791,11792,11793,11794,11795,11796,11797,11798,11799,11800,11801,11804,11805,11806,11807,11808,11809,11810,11811,11812,11813,11814,11815,11816,11817,11818,11819,11820,11821,11822,11823,11824,11825,11826,11827,11828,11829,12183,12184,12185,12186,12187,12188,12189,12190,12191,12192,12193,12194,12195,12196,12197,12198,12199,12200,12201,12202,12203,12204,12205,12206,12207,12208,12211,12212,12213,12214,12215,12216,12217,12218,12219,12220,12221,12222,12223,12224,12225,12226,12227,12228,12229,12230,12231,12232,12233,12234,12235,12236,12594,12595,12596,12597,12598,12599,12600,12601,12602,12603,12604,12605,12606,12607,12608,12609,12610,12611,12612,12613,12614,12615,12616,12617,12618,12619',
       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.T27t06.J080443+653000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T27t06.J080443+653000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T27t06.J080443+653000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#