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

# hif_editimlist(parameter_file='parameter.7166.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_T16t23.J163625+213000_P7166v1_2017_09_29T21_21_33.757/working/VLASS1.1.sb34441820.eb34483341.58017.956089363426_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T16t23.J163625+213000_P7166v1_2017_09_29T21_21_33.757/rawdata/VLASS1.1.sb34441820.eb34483341.58017.956089363426.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='8909,8908,8905,8904,8907,8906,8901,8900,8903,8902,10604,9549,8468,10608,10609,8552,8550,9528,9529,9526,9527,9551,9550,9553,9552,9555,9554,9557,9556,8462,8463,8460,8461,8466,8467,10028,10029,10026,10027,10024,10025,10022,10023,10020,10021,8949,10030,8459,10589,10248,10249,10242,10243,10240,10241,10246,10247,10244,10245,8950,10680,8951,9182,9183,8958,11037,8959,9593,9592,11044,11043,11042,11041,11040,8530,8549,8548,8531,8545,8544,8547,8546,8541,8540,8543,8542,8533,8534,9172,10664,10655,10666,10667,10660,9980,10662,10654,9170,10668,10657,10656,10651,11026,10653,10652,8978,8979,10590,10591,10596,10597,10594,10595,8970,8971,8972,8973,8974,8975,8976,8977,10233,10232,10231,10230,10237,10236,10235,10234,10239,10238,11035,10618,10659,10658,9539,9538,9537,9536,9535,9534,9533,9532,9531,9530,11018,11019,11014,11015,11016,11017,11013,10620,10599,10019,10018,8479,8478,8475,8474,9588,9589,10017,10016,8473,10014,9959,9958,9955,9954,9957,9956,9953,9952,10682,9982,8477,8476,8471,8470,10015,8472,9981,8538,8539,9179,9178,9177,9176,8532,9174,9173,9983,8536,8537,11025,11024,11027,8960,11021,11020,11023,11022,11029,11028,10611,10610,10613,10612,10615,10614,10617,10616,10619,9168,9169,9164,9165,9166,9167,8888,8889,8887,8969,8968,9162,8963,8962,8961,9163,8967,8966,8965,8964,10228,10229,10227,9964,9965,9966,9967,9960,9961,9962,9963,9968,9969,8535,10665,10661,10663,10669,9616,9617,9614,9615,9612,9613,9610,9611,9618,9619,9186,9187,9184,9185,9599,9598,9180,9181,9595,9594,9597,9596,9591,9590,9188,9189,8912,8913,8910,8911,8916,8917,8914,8915,8464,10040,10041,10042,8465,10044,10045,8529,8528,8523,8522,8527,8526,8525,8524,11036,10592,11034,8469,11032,11033,11030,10593,11038,11039,10602,10603,10600,10601,10606,10607,9548,10605,9546,9547,9544,9545,9542,9543,9540,9541,10039,10038,10035,10034,10037,10036,10031,10598,10033,10032,8956,8957,8954,8955,8952,8953,8899,8898,8897,8896,8895,8894,8893,8892,8891,8890,8551,9175,10258,10681,10251,10250,10253,10252,10255,10254,10257,10256,9977,9976,9975,9974,9973,9972,9971,9970,9979,9978,10043,9171,10679,10678,10673,10672,10671,10670,10677,10676,10675,10674,9601,9600,9603,9602,9605,9604,9607,9606,9609,9608,8980,9191,9190,9193,9192,11031,8480,8481,8482,8483,8484,8485,8486,8487,8488,8489',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 16:36:25.612 +21.30.0.0000',
       scan=['172,163,166,134,137,140,143,147,150,169,153,156,159'],
       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_T16t23.J163625+213000_P7166v1_2017_09_29T21_21_33.757/working/VLASS1.1.sb34441820.eb34483341.58017.956089363426_split.ms'],
       imagename='VLASS1.1.ql.T16t23.J163625+213000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='8459,8460,8461,8462,8463,8464,8465,8466,8467,8468,8469,8470,8471,8472,8473,8474,8475,8476,8477,8478,8479,8480,8481,8482,8483,8484,8485,8486,8487,8488,8489,8522,8523,8524,8525,8526,8527,8528,8529,8530,8531,8532,8533,8534,8535,8536,8537,8538,8539,8540,8541,8542,8543,8544,8545,8546,8547,8548,8549,8550,8551,8552,8887,8888,8889,8890,8891,8892,8893,8894,8895,8896,8897,8898,8899,8900,8901,8902,8903,8904,8905,8906,8907,8908,8909,8910,8911,8912,8913,8914,8915,8916,8917,8949,8950,8951,8952,8953,8954,8955,8956,8957,8958,8959,8960,8961,8962,8963,8964,8965,8966,8967,8968,8969,8970,8971,8972,8973,8974,8975,8976,8977,8978,8979,8980,9162,9163,9164,9165,9166,9167,9168,9169,9170,9171,9172,9173,9174,9175,9176,9177,9178,9179,9180,9181,9182,9183,9184,9185,9186,9187,9188,9189,9190,9191,9192,9193,9526,9527,9528,9529,9530,9531,9532,9533,9534,9535,9536,9537,9538,9539,9540,9541,9542,9543,9544,9545,9546,9547,9548,9549,9550,9551,9552,9553,9554,9555,9556,9557,9588,9589,9590,9591,9592,9593,9594,9595,9596,9597,9598,9599,9600,9601,9602,9603,9604,9605,9606,9607,9608,9609,9610,9611,9612,9613,9614,9615,9616,9617,9618,9619,9952,9953,9954,9955,9956,9957,9958,9959,9960,9961,9962,9963,9964,9965,9966,9967,9968,9969,9970,9971,9972,9973,9974,9975,9976,9977,9978,9979,9980,9981,9982,9983,10014,10015,10016,10017,10018,10019,10020,10021,10022,10023,10024,10025,10026,10027,10028,10029,10030,10031,10032,10033,10034,10035,10036,10037,10038,10039,10040,10041,10042,10043,10044,10045,10227,10228,10229,10230,10231,10232,10233,10234,10235,10236,10237,10238,10239,10240,10241,10242,10243,10244,10245,10246,10247,10248,10249,10250,10251,10252,10253,10254,10255,10256,10257,10258,10589,10590,10591,10592,10593,10594,10595,10596,10597,10598,10599,10600,10601,10602,10603,10604,10605,10606,10607,10608,10609,10610,10611,10612,10613,10614,10615,10616,10617,10618,10619,10620,10651,10652,10653,10654,10655,10656,10657,10658,10659,10660,10661,10662,10663,10664,10665,10666,10667,10668,10669,10670,10671,10672,10673,10674,10675,10676,10677,10678,10679,10680,10681,10682,11013,11014,11015,11016,11017,11018,11019,11020,11021,11022,11023,11024,11025,11026,11027,11028,11029,11030,11031,11032,11033,11034,11035,11036,11037,11038,11039,11040,11041,11042,11043,11044',
       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 16:36:25.612 +21.30.0.0000',
       scan=['172,163,166,134,137,140,143,147,150,169,153,156,159'],
       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_T16t23.J163625+213000_P7166v1_2017_09_29T21_21_33.757/working/VLASS1.1.sb34441820.eb34483341.58017.956089363426_split.ms'],
       imagename='VLASS1.1.ql.T16t23.J163625+213000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00052342423837908425, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='8459,8460,8461,8462,8463,8464,8465,8466,8467,8468,8469,8470,8471,8472,8473,8474,8475,8476,8477,8478,8479,8480,8481,8482,8483,8484,8485,8486,8487,8488,8489,8522,8523,8524,8525,8526,8527,8528,8529,8530,8531,8532,8533,8534,8535,8536,8537,8538,8539,8540,8541,8542,8543,8544,8545,8546,8547,8548,8549,8550,8551,8552,8887,8888,8889,8890,8891,8892,8893,8894,8895,8896,8897,8898,8899,8900,8901,8902,8903,8904,8905,8906,8907,8908,8909,8910,8911,8912,8913,8914,8915,8916,8917,8949,8950,8951,8952,8953,8954,8955,8956,8957,8958,8959,8960,8961,8962,8963,8964,8965,8966,8967,8968,8969,8970,8971,8972,8973,8974,8975,8976,8977,8978,8979,8980,9162,9163,9164,9165,9166,9167,9168,9169,9170,9171,9172,9173,9174,9175,9176,9177,9178,9179,9180,9181,9182,9183,9184,9185,9186,9187,9188,9189,9190,9191,9192,9193,9526,9527,9528,9529,9530,9531,9532,9533,9534,9535,9536,9537,9538,9539,9540,9541,9542,9543,9544,9545,9546,9547,9548,9549,9550,9551,9552,9553,9554,9555,9556,9557,9588,9589,9590,9591,9592,9593,9594,9595,9596,9597,9598,9599,9600,9601,9602,9603,9604,9605,9606,9607,9608,9609,9610,9611,9612,9613,9614,9615,9616,9617,9618,9619,9952,9953,9954,9955,9956,9957,9958,9959,9960,9961,9962,9963,9964,9965,9966,9967,9968,9969,9970,9971,9972,9973,9974,9975,9976,9977,9978,9979,9980,9981,9982,9983,10014,10015,10016,10017,10018,10019,10020,10021,10022,10023,10024,10025,10026,10027,10028,10029,10030,10031,10032,10033,10034,10035,10036,10037,10038,10039,10040,10041,10042,10043,10044,10045,10227,10228,10229,10230,10231,10232,10233,10234,10235,10236,10237,10238,10239,10240,10241,10242,10243,10244,10245,10246,10247,10248,10249,10250,10251,10252,10253,10254,10255,10256,10257,10258,10589,10590,10591,10592,10593,10594,10595,10596,10597,10598,10599,10600,10601,10602,10603,10604,10605,10606,10607,10608,10609,10610,10611,10612,10613,10614,10615,10616,10617,10618,10619,10620,10651,10652,10653,10654,10655,10656,10657,10658,10659,10660,10661,10662,10663,10664,10665,10666,10667,10668,10669,10670,10671,10672,10673,10674,10675,10676,10677,10678,10679,10680,10681,10682,11013,11014,11015,11016,11017,11018,11019,11020,11021,11022,11023,11024,11025,11026,11027,11028,11029,11030,11031,11032,11033,11034,11035,11036,11037,11038,11039,11040,11041,11042,11043,11044',
       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.T16t23.J163625+213000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T16t23.J163625+213000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T16t23.J163625+213000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#