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

# hif_editimlist(parameter_file='parameter.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_T11t16.J100200+013000_P19419v1_2018_04_10T23_02_07.590/working/VLASS1.1.sb34739674.eb34813174.58124.37509868055_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T11t16.J100200+013000_P19419v1_2018_04_10T23_02_07.590/rawdata/VLASS1.1.sb34739674.eb34813174.58124.37509868055.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='10141,9735,9742,9368,10140,9366,9367,9364,9365,9362,9363,8150,9361,9319,10131,10142,8549,8548,8545,8544,8547,8546,8541,8540,8543,8542,10128,9737,9736,10148,9739,10147,8909,9738,10537,8922,8149,8148,9734,8141,8140,8143,8142,8145,8144,8147,8146,8538,8539,9360,9722,8531,8532,8533,8534,8535,8536,8537,10557,10554,10543,10555,10164,9726,10533,10534,10535,10536,9727,10538,10539,10551,9724,10129,8528,8138,8139,9774,8523,8522,8521,8137,8130,8131,8132,8133,10132,9771,9721,9327,9763,10134,9773,9762,8962,8961,8960,10172,8516,8517,8514,8515,8512,8513,8510,8511,9733,10138,10139,8518,8519,10545,10544,10547,10546,10541,9728,8129,8128,8127,8126,8125,9729,10542,10548,10180,8956,8957,8954,8955,8952,8953,8950,8951,9760,8958,8959,10549,9761,10169,9766,10151,10540,9731,9730,8503,9732,8505,8504,8507,8506,8509,8508,9318,10173,10145,10556,9317,9328,9329,10552,10553,10550,9316,9322,9323,9320,9321,9326,9315,9324,9325,8941,8940,8943,8942,8945,8944,8947,8946,8949,8948,10168,10167,10127,10157,10166,9335,9334,9337,9336,9331,9330,9333,9332,9764,10143,10159,9339,9338,10162,9744,9745,9746,9740,9741,9725,9743,10130,9749,10133,9765,9756,9755,8934,10155,8937,8930,8931,8932,8933,8938,8939,9340,10152,10558,9343,9344,9345,9346,9347,9348,9349,10165,9750,10163,10150,10161,10160,9757,8134,9723,9754,9753,9752,9751,8135,10156,10135,8136,9759,9758,8520,10144,8527,10137,9770,10178,8526,8927,8926,8925,8924,8923,8525,8921,8920,10136,8524,10158,8929,8928,9353,9352,9351,9350,9357,9356,9355,9354,10170,10171,9359,9358,10174,10175,10176,10177,8552,8553,8550,8551,8556,9767,8554,8555,9768,9769,10146,10179,9772,8912,8913,8910,8911,8916,8917,8914,8915,8918,8919,10149',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 10:02:0.018 +01.30.0.0000',
       scan=['161,164,135,170,139,142,145,148,151,154,167,158'],
       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_T11t16.J100200+013000_P19419v1_2018_04_10T23_02_07.590/working/VLASS1.1.sb34739674.eb34813174.58124.37509868055_split.ms'],
       imagename='VLASS1.1.ql.T11t16.J100200+013000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='8125,8126,8127,8128,8129,8130,8131,8132,8133,8134,8135,8136,8137,8138,8139,8140,8141,8142,8143,8144,8145,8146,8147,8148,8149,8150,8503,8504,8505,8506,8507,8508,8509,8510,8511,8512,8513,8514,8515,8516,8517,8518,8519,8520,8521,8522,8523,8524,8525,8526,8527,8528,8531,8532,8533,8534,8535,8536,8537,8538,8539,8540,8541,8542,8543,8544,8545,8546,8547,8548,8549,8550,8551,8552,8553,8554,8555,8556,8909,8910,8911,8912,8913,8914,8915,8916,8917,8918,8919,8920,8921,8922,8923,8924,8925,8926,8927,8928,8929,8930,8931,8932,8933,8934,8937,8938,8939,8940,8941,8942,8943,8944,8945,8946,8947,8948,8949,8950,8951,8952,8953,8954,8955,8956,8957,8958,8959,8960,8961,8962,9315,9316,9317,9318,9319,9320,9321,9322,9323,9324,9325,9326,9327,9328,9329,9330,9331,9332,9333,9334,9335,9336,9337,9338,9339,9340,9343,9344,9345,9346,9347,9348,9349,9350,9351,9352,9353,9354,9355,9356,9357,9358,9359,9360,9361,9362,9363,9364,9365,9366,9367,9368,9721,9722,9723,9724,9725,9726,9727,9728,9729,9730,9731,9732,9733,9734,9735,9736,9737,9738,9739,9740,9741,9742,9743,9744,9745,9746,9749,9750,9751,9752,9753,9754,9755,9756,9757,9758,9759,9760,9761,9762,9763,9764,9765,9766,9767,9768,9769,9770,9771,9772,9773,9774,10127,10128,10129,10130,10131,10132,10133,10134,10135,10136,10137,10138,10139,10140,10141,10142,10143,10144,10145,10146,10147,10148,10149,10150,10151,10152,10155,10156,10157,10158,10159,10160,10161,10162,10163,10164,10165,10166,10167,10168,10169,10170,10171,10172,10173,10174,10175,10176,10177,10178,10179,10180,10533,10534,10535,10536,10537,10538,10539,10540,10541,10542,10543,10544,10545,10546,10547,10548,10549,10550,10551,10552,10553,10554,10555,10556,10557,10558',
       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 10:02:0.018 +01.30.0.0000',
       scan=['161,164,135,170,139,142,145,148,151,154,167,158'],
       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_T11t16.J100200+013000_P19419v1_2018_04_10T23_02_07.590/working/VLASS1.1.sb34739674.eb34813174.58124.37509868055_split.ms'],
       imagename='VLASS1.1.ql.T11t16.J100200+013000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00053919082078100473, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='8125,8126,8127,8128,8129,8130,8131,8132,8133,8134,8135,8136,8137,8138,8139,8140,8141,8142,8143,8144,8145,8146,8147,8148,8149,8150,8503,8504,8505,8506,8507,8508,8509,8510,8511,8512,8513,8514,8515,8516,8517,8518,8519,8520,8521,8522,8523,8524,8525,8526,8527,8528,8531,8532,8533,8534,8535,8536,8537,8538,8539,8540,8541,8542,8543,8544,8545,8546,8547,8548,8549,8550,8551,8552,8553,8554,8555,8556,8909,8910,8911,8912,8913,8914,8915,8916,8917,8918,8919,8920,8921,8922,8923,8924,8925,8926,8927,8928,8929,8930,8931,8932,8933,8934,8937,8938,8939,8940,8941,8942,8943,8944,8945,8946,8947,8948,8949,8950,8951,8952,8953,8954,8955,8956,8957,8958,8959,8960,8961,8962,9315,9316,9317,9318,9319,9320,9321,9322,9323,9324,9325,9326,9327,9328,9329,9330,9331,9332,9333,9334,9335,9336,9337,9338,9339,9340,9343,9344,9345,9346,9347,9348,9349,9350,9351,9352,9353,9354,9355,9356,9357,9358,9359,9360,9361,9362,9363,9364,9365,9366,9367,9368,9721,9722,9723,9724,9725,9726,9727,9728,9729,9730,9731,9732,9733,9734,9735,9736,9737,9738,9739,9740,9741,9742,9743,9744,9745,9746,9749,9750,9751,9752,9753,9754,9755,9756,9757,9758,9759,9760,9761,9762,9763,9764,9765,9766,9767,9768,9769,9770,9771,9772,9773,9774,10127,10128,10129,10130,10131,10132,10133,10134,10135,10136,10137,10138,10139,10140,10141,10142,10143,10144,10145,10146,10147,10148,10149,10150,10151,10152,10155,10156,10157,10158,10159,10160,10161,10162,10163,10164,10165,10166,10167,10168,10169,10170,10171,10172,10173,10174,10175,10176,10177,10178,10179,10180,10533,10534,10535,10536,10537,10538,10539,10540,10541,10542,10543,10544,10545,10546,10547,10548,10549,10550,10551,10552,10553,10554,10555,10556,10557,10558',
       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.T11t16.J100200+013000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T11t16.J100200+013000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T11t16.J100200+013000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#