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

# hif_editimlist(parameter_file='parameter.15444.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_T28t08.J141644+693000_P15444v1_2018_01_12T18_19_26.940/working/VLASS1.1.sb34613347.eb34742769.58081.77578660879_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T28t08.J141644+693000_P15444v1_2018_01_12T18_19_26.940/rawdata/VLASS1.1.sb34613347.eb34742769.58081.77578660879.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='9076,8255,9074,9075,9072,9073,9070,9071,9078,9079,8307,8306,8305,8304,8303,8302,8301,8300,8251,8309,8250,8253,8252,6796,6797,6794,6795,6792,6793,6790,6791,6798,6799,8233,8232,8231,8230,8237,8236,8235,8234,8239,8238,7889,7888,7881,7880,8734,8308,7885,7884,8730,8731,7867,7866,7865,7864,7863,7869,8289,8710,8658,6808,6809,6806,6807,6804,6805,6802,6803,6800,6801,8718,8719,7175,7174,7177,7176,7171,7170,7173,7172,7179,7178,8248,8249,8242,8243,8240,8241,8246,7826,8244,7824,8648,8649,8646,8647,7182,8645,7184,7185,7186,7187,7221,7220,7223,7222,7225,7224,7227,7226,7229,7228,7456,7457,7454,7455,7452,7453,7450,7451,7458,7459,8735,9061,9063,9062,9065,9064,9067,9066,9069,9068,8732,8290,8296,6783,6785,6784,6787,6786,6789,6788,7829,8729,8728,7828,8721,8720,8723,8722,8725,8724,8727,8726,8291,7823,8293,8292,8295,8294,8297,7822,8299,8298,7821,7820,7827,8247,7825,8245,7218,7219,6811,6810,6813,6812,7180,7181,7183,9077,7882,7816,7817,7814,7815,7812,7813,7810,7811,7818,7819,8659,8711,8712,8713,8714,8715,8716,8717,8651,8650,8653,8652,8655,8654,8657,8656,7883,7441,7440,7443,7442,7445,7444,7447,7446,7449,7448,8733,7887,7886,9090,9091,7247,7246,7245,7244,7243,7242,7241,7240,7249,7248,7434,7435,7436,7437,7438,7439,8286,8287,7159,7158,7157,7156,8288,7868,8660,8661,8662,8663,8664,8665,8666,8667,8668,8669,8226,7806,8228,8229,8254,8314,8315,8316,8317,8310,8311,8312,8313,8256,7803,7802,7805,7804,7807,8227,7809,7808,8707,8706,8705,8709,8708,7892,7893,7890,7891,7874,7875,7876,7877,7870,7871,7872,7873,7878,7879,9089,9088,9087,9086,9085,9084,9083,9082,9081,9080,7168,7169,7162,7163,7160,7161,7166,7167,7164,7165,7830,7831,7832,8673,8672,8671,8670,8676,8675,8674,7236,7237,7234,7235,7232,7233,7230,7231,7238,7239,7463,7462,7461,7460,7464',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 14:16:44.554 +69.30.0.0000',
       scan=['129,132,135,138,141,145,116,119,148,122,126'],
       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_T28t08.J141644+693000_P15444v1_2018_01_12T18_19_26.940/working/VLASS1.1.sb34613347.eb34742769.58081.77578660879_split.ms'],
       imagename='VLASS1.1.ql.T28t08.J141644+693000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='6783,6784,6785,6786,6787,6788,6789,6790,6791,6792,6793,6794,6795,6796,6797,6798,6799,6800,6801,6802,6803,6804,6805,6806,6807,6808,6809,6810,6811,6812,6813,7156,7157,7158,7159,7160,7161,7162,7163,7164,7165,7166,7167,7168,7169,7170,7171,7172,7173,7174,7175,7176,7177,7178,7179,7180,7181,7182,7183,7184,7185,7186,7187,7218,7219,7220,7221,7222,7223,7224,7225,7226,7227,7228,7229,7230,7231,7232,7233,7234,7235,7236,7237,7238,7239,7240,7241,7242,7243,7244,7245,7246,7247,7248,7249,7434,7435,7436,7437,7438,7439,7440,7441,7442,7443,7444,7445,7446,7447,7448,7449,7450,7451,7452,7453,7454,7455,7456,7457,7458,7459,7460,7461,7462,7463,7464,7802,7803,7804,7805,7806,7807,7808,7809,7810,7811,7812,7813,7814,7815,7816,7817,7818,7819,7820,7821,7822,7823,7824,7825,7826,7827,7828,7829,7830,7831,7832,7863,7864,7865,7866,7867,7868,7869,7870,7871,7872,7873,7874,7875,7876,7877,7878,7879,7880,7881,7882,7883,7884,7885,7886,7887,7888,7889,7890,7891,7892,7893,8226,8227,8228,8229,8230,8231,8232,8233,8234,8235,8236,8237,8238,8239,8240,8241,8242,8243,8244,8245,8246,8247,8248,8249,8250,8251,8252,8253,8254,8255,8256,8286,8287,8288,8289,8290,8291,8292,8293,8294,8295,8296,8297,8298,8299,8300,8301,8302,8303,8304,8305,8306,8307,8308,8309,8310,8311,8312,8313,8314,8315,8316,8317,8645,8646,8647,8648,8649,8650,8651,8652,8653,8654,8655,8656,8657,8658,8659,8660,8661,8662,8663,8664,8665,8666,8667,8668,8669,8670,8671,8672,8673,8674,8675,8676,8705,8706,8707,8708,8709,8710,8711,8712,8713,8714,8715,8716,8717,8718,8719,8720,8721,8722,8723,8724,8725,8726,8727,8728,8729,8730,8731,8732,8733,8734,8735,9061,9062,9063,9064,9065,9066,9067,9068,9069,9070,9071,9072,9073,9074,9075,9076,9077,9078,9079,9080,9081,9082,9083,9084,9085,9086,9087,9088,9089,9090,9091',
       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 14:16:44.554 +69.30.0.0000',
       scan=['129,132,135,138,141,145,116,119,148,122,126'],
       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_T28t08.J141644+693000_P15444v1_2018_01_12T18_19_26.940/working/VLASS1.1.sb34613347.eb34742769.58081.77578660879_split.ms'],
       imagename='VLASS1.1.ql.T28t08.J141644+693000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00052288770331651905, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='6783,6784,6785,6786,6787,6788,6789,6790,6791,6792,6793,6794,6795,6796,6797,6798,6799,6800,6801,6802,6803,6804,6805,6806,6807,6808,6809,6810,6811,6812,6813,7156,7157,7158,7159,7160,7161,7162,7163,7164,7165,7166,7167,7168,7169,7170,7171,7172,7173,7174,7175,7176,7177,7178,7179,7180,7181,7182,7183,7184,7185,7186,7187,7218,7219,7220,7221,7222,7223,7224,7225,7226,7227,7228,7229,7230,7231,7232,7233,7234,7235,7236,7237,7238,7239,7240,7241,7242,7243,7244,7245,7246,7247,7248,7249,7434,7435,7436,7437,7438,7439,7440,7441,7442,7443,7444,7445,7446,7447,7448,7449,7450,7451,7452,7453,7454,7455,7456,7457,7458,7459,7460,7461,7462,7463,7464,7802,7803,7804,7805,7806,7807,7808,7809,7810,7811,7812,7813,7814,7815,7816,7817,7818,7819,7820,7821,7822,7823,7824,7825,7826,7827,7828,7829,7830,7831,7832,7863,7864,7865,7866,7867,7868,7869,7870,7871,7872,7873,7874,7875,7876,7877,7878,7879,7880,7881,7882,7883,7884,7885,7886,7887,7888,7889,7890,7891,7892,7893,8226,8227,8228,8229,8230,8231,8232,8233,8234,8235,8236,8237,8238,8239,8240,8241,8242,8243,8244,8245,8246,8247,8248,8249,8250,8251,8252,8253,8254,8255,8256,8286,8287,8288,8289,8290,8291,8292,8293,8294,8295,8296,8297,8298,8299,8300,8301,8302,8303,8304,8305,8306,8307,8308,8309,8310,8311,8312,8313,8314,8315,8316,8317,8645,8646,8647,8648,8649,8650,8651,8652,8653,8654,8655,8656,8657,8658,8659,8660,8661,8662,8663,8664,8665,8666,8667,8668,8669,8670,8671,8672,8673,8674,8675,8676,8705,8706,8707,8708,8709,8710,8711,8712,8713,8714,8715,8716,8717,8718,8719,8720,8721,8722,8723,8724,8725,8726,8727,8728,8729,8730,8731,8732,8733,8734,8735,9061,9062,9063,9064,9065,9066,9067,9068,9069,9070,9071,9072,9073,9074,9075,9076,9077,9078,9079,9080,9081,9082,9083,9084,9085,9086,9087,9088,9089,9090,9091',
       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.T28t08.J141644+693000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T28t08.J141644+693000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T28t08.J141644+693000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#