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.sb34477887.eb34601198.58042.00525377315.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.sb34477887.eb34601198.58042.00525377315_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T10t32.J205400-013000_P9613v1_2018_04_25T22_59_26.844/rawdata/VLASS1.1.sb34477887.eb34601198.58042.00525377315.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='8465,7810,7811,6969,8605,8621,8622,8623,7388,7389,8871,8882,7382,7383,7380,7381,7386,7387,7384,7385,8866,8867,8864,8476,8862,8863,8860,8861,8868,8869,8211,8473,8213,7395,7394,7397,7396,6998,6999,7393,7392,6994,6995,6996,6997,6990,6991,6992,6993,8451,8872,8877,8875,8874,7254,7255,7256,7257,7250,7251,7252,7253,8879,8878,8876,7258,7259,8450,8618,8461,8469,8466,8873,8467,8608,8609,8464,6989,6988,6987,6986,6985,6984,6983,6982,6981,6980,7791,8883,7790,8468,7247,7246,7245,7244,7243,7242,7241,7240,8884,8885,8886,8209,8880,8881,7249,7248,8206,8612,8207,8044,8045,8046,8047,8048,8049,8204,8205,8202,8203,8200,8201,8457,8456,8455,8454,8453,8452,7788,7789,7786,7787,7784,7785,7782,7783,8459,7781,7379,7378,7639,7398,7668,8865,7377,7376,7375,8463,8053,8052,8051,8050,8057,8056,8055,8054,8074,8210,8059,8058,8215,8214,8217,8216,7650,7651,7652,7653,7654,7655,7656,7657,7658,7659,7793,7792,7795,7794,7797,7796,8192,8193,8190,8191,8196,8197,8194,8195,8198,8199,7261,7260,8462,7262,8615,8602,8617,8616,8611,7667,8598,8599,8212,8597,8594,7666,8619,8593,8610,8601,8613,8068,8069,8066,8067,8064,8065,8062,8063,7799,8061,7643,7642,7641,7640,7647,7646,7645,7644,8475,8474,7649,7648,8471,8470,8458,8472,8595,8189,8188,8870,8187,7798,8596,7805,8600,7405,7404,8607,7401,7400,7403,7402,7801,8460,7800,8071,8070,8073,8072,7803,7804,8620,7802,8060,7391,7390,7807,7806,7809,8606,8603,7808,8208,6976,6977,6974,6975,6972,6973,6970,6971,7399,8614,7638,6978,6979,7236,7237,7234,7235,7232,7233,7661,7660,7663,7662,7665,7664,7238,7239,8477,8604,8480,8859,8858,8479,8857,8856,8478',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 20:54:0.128 -01.30.0.0000',
       scan=['160,163,166,169,141,144,147,150,153,157'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34477887.eb34601198.58042.00525377315_split.ms'],
       imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='6969,6970,6971,6972,6973,6974,6975,6976,6977,6978,6979,6980,6981,6982,6983,6984,6985,6986,6987,6988,6989,6990,6991,6992,6993,6994,6995,6996,6997,6998,6999,7232,7233,7234,7235,7236,7237,7238,7239,7240,7241,7242,7243,7244,7245,7246,7247,7248,7249,7250,7251,7252,7253,7254,7255,7256,7257,7258,7259,7260,7261,7262,7375,7376,7377,7378,7379,7380,7381,7382,7383,7384,7385,7386,7387,7388,7389,7390,7391,7392,7393,7394,7395,7396,7397,7398,7399,7400,7401,7402,7403,7404,7405,7638,7639,7640,7641,7642,7643,7644,7645,7646,7647,7648,7649,7650,7651,7652,7653,7654,7655,7656,7657,7658,7659,7660,7661,7662,7663,7664,7665,7666,7667,7668,7781,7782,7783,7784,7785,7786,7787,7788,7789,7790,7791,7792,7793,7794,7795,7796,7797,7798,7799,7800,7801,7802,7803,7804,7805,7806,7807,7808,7809,7810,7811,8044,8045,8046,8047,8048,8049,8050,8051,8052,8053,8054,8055,8056,8057,8058,8059,8060,8061,8062,8063,8064,8065,8066,8067,8068,8069,8070,8071,8072,8073,8074,8187,8188,8189,8190,8191,8192,8193,8194,8195,8196,8197,8198,8199,8200,8201,8202,8203,8204,8205,8206,8207,8208,8209,8210,8211,8212,8213,8214,8215,8216,8217,8450,8451,8452,8453,8454,8455,8456,8457,8458,8459,8460,8461,8462,8463,8464,8465,8466,8467,8468,8469,8470,8471,8472,8473,8474,8475,8476,8477,8478,8479,8480,8593,8594,8595,8596,8597,8598,8599,8600,8601,8602,8603,8604,8605,8606,8607,8608,8609,8610,8611,8612,8613,8614,8615,8616,8617,8618,8619,8620,8621,8622,8623,8856,8857,8858,8859,8860,8861,8862,8863,8864,8865,8866,8867,8868,8869,8870,8871,8872,8873,8874,8875,8876,8877,8878,8879,8880,8881,8882,8883,8884,8885,8886',
       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 20:54:0.128 -01.30.0.0000',
       scan=['160,163,166,169,141,144,147,150,153,157'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34477887.eb34601198.58042.00525377315_split.ms'],
       imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00067597185500562764, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='6969,6970,6971,6972,6973,6974,6975,6976,6977,6978,6979,6980,6981,6982,6983,6984,6985,6986,6987,6988,6989,6990,6991,6992,6993,6994,6995,6996,6997,6998,6999,7232,7233,7234,7235,7236,7237,7238,7239,7240,7241,7242,7243,7244,7245,7246,7247,7248,7249,7250,7251,7252,7253,7254,7255,7256,7257,7258,7259,7260,7261,7262,7375,7376,7377,7378,7379,7380,7381,7382,7383,7384,7385,7386,7387,7388,7389,7390,7391,7392,7393,7394,7395,7396,7397,7398,7399,7400,7401,7402,7403,7404,7405,7638,7639,7640,7641,7642,7643,7644,7645,7646,7647,7648,7649,7650,7651,7652,7653,7654,7655,7656,7657,7658,7659,7660,7661,7662,7663,7664,7665,7666,7667,7668,7781,7782,7783,7784,7785,7786,7787,7788,7789,7790,7791,7792,7793,7794,7795,7796,7797,7798,7799,7800,7801,7802,7803,7804,7805,7806,7807,7808,7809,7810,7811,8044,8045,8046,8047,8048,8049,8050,8051,8052,8053,8054,8055,8056,8057,8058,8059,8060,8061,8062,8063,8064,8065,8066,8067,8068,8069,8070,8071,8072,8073,8074,8187,8188,8189,8190,8191,8192,8193,8194,8195,8196,8197,8198,8199,8200,8201,8202,8203,8204,8205,8206,8207,8208,8209,8210,8211,8212,8213,8214,8215,8216,8217,8450,8451,8452,8453,8454,8455,8456,8457,8458,8459,8460,8461,8462,8463,8464,8465,8466,8467,8468,8469,8470,8471,8472,8473,8474,8475,8476,8477,8478,8479,8480,8593,8594,8595,8596,8597,8598,8599,8600,8601,8602,8603,8604,8605,8606,8607,8608,8609,8610,8611,8612,8613,8614,8615,8616,8617,8618,8619,8620,8621,8622,8623,8856,8857,8858,8859,8860,8861,8862,8863,8864,8865,8866,8867,8868,8869,8870,8871,8872,8873,8874,8875,8876,8877,8878,8879,8880,8881,8882,8883,8884,8885,8886',
       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.T10t32.J205400-013000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T10t32.J205400-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.T10t32.J205400-013000.10.2048.v1.I.iter1.image.pbcor.tt0.rms',
      imagename='VLASS1.1.ql.T10t32.J205400-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.T10t32.J205400-013000.10.2048.v1.I.iter1.image.pbcor.tt0')
imhead(imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.image.pbcor.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.image.pbcor.tt0.subim',
           imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.image.pbcor.tt0')
imstat(imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imhead(imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.image.pbcor.tt0.rms.subim',
           imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imstat(imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.psf.tt0')
imhead(imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.psf.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.psf.tt0.subim',
           imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.psf.tt0')
imstat(imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.residual.tt0')
imhead(imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.residual.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.residual.tt0.subim',
           imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.residual.tt0')
imstat(imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.pb.tt0')
imhead(imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.pb.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.pb.tt0.subim',
           imagename='VLASS1.1.ql.T10t32.J205400-013000.10.2048.v1.I.iter1.pb.tt0')

# hifv_exportvlassdata(pipelinemode="automatic")
#