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

# hif_editimlist(parameter_file='parameter.13143.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_T27t12.J190642+683000_P13143v1_2017_10_17T14_53_25.018/working/VLASS1.1.sb34538714.eb34596436.58039.97741366898_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T27t12.J190642+683000_P13143v1_2017_10_17T14_53_25.018/rawdata/VLASS1.1.sb34538714.eb34596436.58039.97741366898.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='11191,11190,11193,11192,11195,11194,11197,11196,11199,11198,12609,11542,11543,11540,11541,11546,11547,11544,11545,11561,11548,12613,12279,12278,11913,11912,10817,11914,11917,11916,11919,11918,10819,10818,12275,12274,12277,12276,11182,11183,11180,11181,11186,11187,11188,11189,11539,12608,11576,11575,11574,11573,11572,11571,11570,11902,11903,11900,11901,11906,11907,11904,11905,12262,12263,11908,11909,12266,12264,12265,11201,11200,11203,11202,11205,11204,11207,11206,11209,11208,11562,11563,11564,11565,11566,11567,12257,12256,12255,12254,12253,12252,12251,12250,12259,12258,12269,12260,11936,11935,11934,11933,11932,11931,11930,12595,12596,12261,12598,12290,12599,11210,11211,11898,11899,11894,11895,11896,11897,11890,11891,11892,11893,12284,12242,12243,12244,12245,12246,12247,12248,12249,12286,11928,11929,12287,11920,11921,11922,11923,11924,11925,11926,11927,12282,11568,12283,12600,12601,12597,12604,11159,11158,12605,11915,12272,11888,12271,12270,12273,10840,10841,10842,11538,11569,11533,11532,11531,11530,11537,11536,11535,11534,11887,11886,11885,12618,11524,11525,11526,11527,12619,11523,12617,12616,12615,12614,11528,11529,12611,12610,11173,11172,11171,11170,11177,11176,11175,11174,12293,12292,11179,11178,12280,12612,10835,10834,10837,10836,10831,10830,10833,10832,10839,10838,12281,11551,12285,11553,11552,11555,11554,11557,11556,11559,11558,12602,12603,12288,12289,12606,12607,11168,11169,11889,11164,11165,11166,11167,11160,11161,11162,11163,10826,10827,10824,10825,10822,10823,10820,10821,11560,12291,10828,10829',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 19:06:42.688 +68.30.0.0000',
       scan=['194,197,201,204,207,210,213,216,188,191'],
       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_T27t12.J190642+683000_P13143v1_2017_10_17T14_53_25.018/working/VLASS1.1.sb34538714.eb34596436.58039.97741366898_split.ms'],
       imagename='VLASS1.1.ql.T27t12.J190642+683000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='10817,10818,10819,10820,10821,10822,10823,10824,10825,10826,10827,10828,10829,10830,10831,10832,10833,10834,10835,10836,10837,10838,10839,10840,10841,10842,11158,11159,11160,11161,11162,11163,11164,11165,11166,11167,11168,11169,11170,11171,11172,11173,11174,11175,11176,11177,11178,11179,11180,11181,11182,11183,11186,11187,11188,11189,11190,11191,11192,11193,11194,11195,11196,11197,11198,11199,11200,11201,11202,11203,11204,11205,11206,11207,11208,11209,11210,11211,11523,11524,11525,11526,11527,11528,11529,11530,11531,11532,11533,11534,11535,11536,11537,11538,11539,11540,11541,11542,11543,11544,11545,11546,11547,11548,11551,11552,11553,11554,11555,11556,11557,11558,11559,11560,11561,11562,11563,11564,11565,11566,11567,11568,11569,11570,11571,11572,11573,11574,11575,11576,11885,11886,11887,11888,11889,11890,11891,11892,11893,11894,11895,11896,11897,11898,11899,11900,11901,11902,11903,11904,11905,11906,11907,11908,11909,11912,11913,11914,11915,11916,11917,11918,11919,11920,11921,11922,11923,11924,11925,11926,11927,11928,11929,11930,11931,11932,11933,11934,11935,11936,12242,12243,12244,12245,12246,12247,12248,12249,12250,12251,12252,12253,12254,12255,12256,12257,12258,12259,12260,12261,12262,12263,12264,12265,12266,12269,12270,12271,12272,12273,12274,12275,12276,12277,12278,12279,12280,12281,12282,12283,12284,12285,12286,12287,12288,12289,12290,12291,12292,12293,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 19:06:42.688 +68.30.0.0000',
       scan=['194,197,201,204,207,210,213,216,188,191'],
       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_T27t12.J190642+683000_P13143v1_2017_10_17T14_53_25.018/working/VLASS1.1.sb34538714.eb34596436.58039.97741366898_split.ms'],
       imagename='VLASS1.1.ql.T27t12.J190642+683000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00054917689928224715, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='10817,10818,10819,10820,10821,10822,10823,10824,10825,10826,10827,10828,10829,10830,10831,10832,10833,10834,10835,10836,10837,10838,10839,10840,10841,10842,11158,11159,11160,11161,11162,11163,11164,11165,11166,11167,11168,11169,11170,11171,11172,11173,11174,11175,11176,11177,11178,11179,11180,11181,11182,11183,11186,11187,11188,11189,11190,11191,11192,11193,11194,11195,11196,11197,11198,11199,11200,11201,11202,11203,11204,11205,11206,11207,11208,11209,11210,11211,11523,11524,11525,11526,11527,11528,11529,11530,11531,11532,11533,11534,11535,11536,11537,11538,11539,11540,11541,11542,11543,11544,11545,11546,11547,11548,11551,11552,11553,11554,11555,11556,11557,11558,11559,11560,11561,11562,11563,11564,11565,11566,11567,11568,11569,11570,11571,11572,11573,11574,11575,11576,11885,11886,11887,11888,11889,11890,11891,11892,11893,11894,11895,11896,11897,11898,11899,11900,11901,11902,11903,11904,11905,11906,11907,11908,11909,11912,11913,11914,11915,11916,11917,11918,11919,11920,11921,11922,11923,11924,11925,11926,11927,11928,11929,11930,11931,11932,11933,11934,11935,11936,12242,12243,12244,12245,12246,12247,12248,12249,12250,12251,12252,12253,12254,12255,12256,12257,12258,12259,12260,12261,12262,12263,12264,12265,12266,12269,12270,12271,12272,12273,12274,12275,12276,12277,12278,12279,12280,12281,12282,12283,12284,12285,12286,12287,12288,12289,12290,12291,12292,12293,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.T27t12.J190642+683000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T27t12.J190642+683000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T27t12.J190642+683000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#