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.sb34905760.eb35113377.58166.981752928245.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='VLASS1.1.sb34905760.eb35113377.58166.981752928245_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T05t09.J055610-203000_P23962v1_2018_03_23T20_04_59.822/rawdata/VLASS1.1.sb34905760.eb35113377.58166.981752928245.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='7986,7993,7992,7998,9328,9329,9326,9327,9325,8941,8940,8943,8942,8945,8197,8947,8946,8949,8948,8198,7553,7552,7551,7550,7557,7556,7555,7554,8955,8953,8934,8935,8936,8937,8930,8931,8932,8933,8938,8939,7175,7174,7177,7176,7171,7170,7173,7172,7179,7178,8167,8166,8169,8168,7188,7189,7180,7181,7182,7183,7184,7185,7186,7187,7957,7956,7955,7954,7953,7952,7951,7950,8545,8544,8547,8546,8541,7968,7959,7958,7969,7597,7596,7595,7594,7593,7592,7591,7590,8978,8979,7526,7527,8970,8971,7528,7529,8974,8975,8976,8977,7944,7945,7940,7980,7981,7982,7983,7984,7985,8598,7987,8596,7989,8594,8595,8592,8593,8590,8591,9335,9334,9337,9336,9331,9330,9333,9332,9339,9338,8574,7988,8576,8577,8570,8571,8572,8573,8578,8579,7566,7567,7565,8189,8188,8185,8184,8187,8186,8181,8180,8183,8182,7941,8927,8929,8928,8575,7193,7192,7191,7190,7195,7194,8538,8539,7966,7967,8530,8531,8532,8533,8534,8535,8536,8537,8985,8984,7531,8987,7530,8986,8549,8548,8983,8982,7568,7569,8989,8988,8969,8968,7539,7538,7535,7534,7537,7536,8967,8966,7533,7532,7931,7930,7933,7932,7935,7934,7937,7936,7939,7938,7927,7928,7929,8589,8588,7991,7990,7997,7996,7995,7994,8581,8580,8583,8582,8585,8584,8587,8586,9340,9341,9342,9343,9344,9345,9346,9347,9348,9349,8567,8566,8569,8568,7571,7570,7573,7572,7575,7574,7577,7576,7579,7578,8529,8528,7979,7978,7975,7974,7977,7976,8527,7970,7973,7972,8192,8193,8190,8191,8196,8972,8944,8956,8957,8954,8973,8952,8194,8950,8951,8195,7946,8958,8959,7548,7549,7947,7540,7541,7542,7543,7544,7545,7546,7547,7942,7943,7587,7580,7168,7169,7166,7167,7164,7165,7584,7585,8540,7586,8170,8171,8172,8173,8174,8175,8176,8177,8178,8179,8597,7581,8994,8543,7583,7971,8981,8542,9353,9352,9351,9350,9357,9356,9355,9354,8980,8552,8553,8550,8551,8556,8557,8554,8555,8558,8559,7948,7949,8998,7588,7589,8992,8993,8990,8991,8996,8997,7582,8995',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 05:56:10.923 -20.30.0.0000',
       scan=['131,147,134,137,140,143,115,118,121,124,127'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34905760.eb35113377.58166.981752928245_split.ms'],
       imagename='VLASS1.1.ql.T05t09.J055610-203000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='7164,7165,7166,7167,7168,7169,7170,7171,7172,7173,7174,7175,7176,7177,7178,7179,7180,7181,7182,7183,7184,7185,7186,7187,7188,7189,7190,7191,7192,7193,7194,7195,7526,7527,7528,7529,7530,7531,7532,7533,7534,7535,7536,7537,7538,7539,7540,7541,7542,7543,7544,7545,7546,7547,7548,7549,7550,7551,7552,7553,7554,7555,7556,7557,7565,7566,7567,7568,7569,7570,7571,7572,7573,7574,7575,7576,7577,7578,7579,7580,7581,7582,7583,7584,7585,7586,7587,7588,7589,7590,7591,7592,7593,7594,7595,7596,7597,7927,7928,7929,7930,7931,7932,7933,7934,7935,7936,7937,7938,7939,7940,7941,7942,7943,7944,7945,7946,7947,7948,7949,7950,7951,7952,7953,7954,7955,7956,7957,7958,7959,7966,7967,7968,7969,7970,7971,7972,7973,7974,7975,7976,7977,7978,7979,7980,7981,7982,7983,7984,7985,7986,7987,7988,7989,7990,7991,7992,7993,7994,7995,7996,7997,7998,8166,8167,8168,8169,8170,8171,8172,8173,8174,8175,8176,8177,8178,8179,8180,8181,8182,8183,8184,8185,8186,8187,8188,8189,8190,8191,8192,8193,8194,8195,8196,8197,8198,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,8553,8554,8555,8556,8557,8558,8559,8566,8567,8568,8569,8570,8571,8572,8573,8574,8575,8576,8577,8578,8579,8580,8581,8582,8583,8584,8585,8586,8587,8588,8589,8590,8591,8592,8593,8594,8595,8596,8597,8598,8927,8928,8929,8930,8931,8932,8933,8934,8935,8936,8937,8938,8939,8940,8941,8942,8943,8944,8945,8946,8947,8948,8949,8950,8951,8952,8953,8954,8955,8956,8957,8958,8959,8966,8967,8968,8969,8970,8971,8972,8973,8974,8975,8976,8977,8978,8979,8980,8981,8982,8983,8984,8985,8986,8987,8988,8989,8990,8991,8992,8993,8994,8995,8996,8997,8998,9325,9326,9327,9328,9329,9330,9331,9332,9333,9334,9335,9336,9337,9338,9339,9340,9341,9342,9343,9344,9345,9346,9347,9348,9349,9350,9351,9352,9353,9354,9355,9356,9357',
       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 05:56:10.923 -20.30.0.0000',
       scan=['131,147,134,137,140,143,115,118,121,124,127'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34905760.eb35113377.58166.981752928245_split.ms'],
       imagename='VLASS1.1.ql.T05t09.J055610-203000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00049595500642517438, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='7164,7165,7166,7167,7168,7169,7170,7171,7172,7173,7174,7175,7176,7177,7178,7179,7180,7181,7182,7183,7184,7185,7186,7187,7188,7189,7190,7191,7192,7193,7194,7195,7526,7527,7528,7529,7530,7531,7532,7533,7534,7535,7536,7537,7538,7539,7540,7541,7542,7543,7544,7545,7546,7547,7548,7549,7550,7551,7552,7553,7554,7555,7556,7557,7565,7566,7567,7568,7569,7570,7571,7572,7573,7574,7575,7576,7577,7578,7579,7580,7581,7582,7583,7584,7585,7586,7587,7588,7589,7590,7591,7592,7593,7594,7595,7596,7597,7927,7928,7929,7930,7931,7932,7933,7934,7935,7936,7937,7938,7939,7940,7941,7942,7943,7944,7945,7946,7947,7948,7949,7950,7951,7952,7953,7954,7955,7956,7957,7958,7959,7966,7967,7968,7969,7970,7971,7972,7973,7974,7975,7976,7977,7978,7979,7980,7981,7982,7983,7984,7985,7986,7987,7988,7989,7990,7991,7992,7993,7994,7995,7996,7997,7998,8166,8167,8168,8169,8170,8171,8172,8173,8174,8175,8176,8177,8178,8179,8180,8181,8182,8183,8184,8185,8186,8187,8188,8189,8190,8191,8192,8193,8194,8195,8196,8197,8198,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,8553,8554,8555,8556,8557,8558,8559,8566,8567,8568,8569,8570,8571,8572,8573,8574,8575,8576,8577,8578,8579,8580,8581,8582,8583,8584,8585,8586,8587,8588,8589,8590,8591,8592,8593,8594,8595,8596,8597,8598,8927,8928,8929,8930,8931,8932,8933,8934,8935,8936,8937,8938,8939,8940,8941,8942,8943,8944,8945,8946,8947,8948,8949,8950,8951,8952,8953,8954,8955,8956,8957,8958,8959,8966,8967,8968,8969,8970,8971,8972,8973,8974,8975,8976,8977,8978,8979,8980,8981,8982,8983,8984,8985,8986,8987,8988,8989,8990,8991,8992,8993,8994,8995,8996,8997,8998,9325,9326,9327,9328,9329,9330,9331,9332,9333,9334,9335,9336,9337,9338,9339,9340,9341,9342,9343,9344,9345,9346,9347,9348,9349,9350,9351,9352,9353,9354,9355,9356,9357',
       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.T05t09.J055610-203000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T05t09.J055610-203000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T05t09.J055610-203000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#