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

# hif_editimlist(parameter_file='parameter.15167.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_T13t32.J211420+083000_P15167v1_2017_11_10T19_34_13.099/working/VLASS1.1.sb34607978.eb34635523.58053.05546966435_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T13t32.J211420+083000_P15167v1_2017_11_10T19_34_13.099/rawdata/VLASS1.1.sb34607978.eb34635523.58053.05546966435.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='7252,8254,8257,8256,8251,8250,8253,8252,8259,8258,8264,8789,8788,8260,8265,8783,8782,8787,8266,8785,8609,7861,8267,7599,7598,7597,7596,7595,7594,7593,7592,7591,7590,7199,7198,8262,8188,7193,7192,7191,7190,7197,7196,7195,7194,7801,7254,7803,7604,7805,7804,7807,7255,7859,7256,6848,6849,8261,6842,6843,6840,6841,6846,6847,6844,6845,6796,6797,6794,6795,6792,6793,6790,6791,7606,7605,7253,7258,7259,6798,6799,8793,8187,8796,6855,6854,6857,6856,6851,6850,6853,6852,8589,7603,8239,8238,6859,6858,8604,7602,8263,7601,8805,7600,7858,6781,6780,6783,6782,6785,6784,6787,6786,6789,6788,8255,7249,7248,8804,8790,8806,8807,8800,8801,8802,8803,8784,8808,8809,8208,8209,8206,8207,8204,8205,8202,8203,8200,8201,7866,7865,7864,7250,7862,7788,7789,7786,7787,7784,7785,7782,7783,7780,7781,8798,6778,6779,8799,6774,6775,6776,6777,8248,6772,6773,8249,8812,8811,8810,7863,7852,6833,6832,6831,6837,6836,6835,6834,8794,6839,6838,8786,8608,7799,7798,7850,7851,7856,7857,7854,7855,7791,7790,7793,7792,7795,7794,7797,7796,8192,8193,8190,7579,8196,8197,8194,8195,8246,8198,8199,7261,7260,7263,7262,8602,8791,8611,8244,8598,8599,8596,8597,8593,8245,8592,8583,8590,7860,7242,8601,7842,8605,8268,8606,6802,6800,6801,7802,7849,7848,7845,7844,7847,7846,7841,7840,7843,7589,7838,8595,8189,7853,7839,8185,8184,7251,8186,8181,8180,8183,8182,8191,8588,7257,8600,7585,8581,8247,8582,8585,8584,8587,8586,8594,7587,8792,7800,8179,7175,7174,7177,7176,8610,8243,8795,7179,7178,7836,7837,8607,7577,7576,7203,7202,7201,7200,7204,7777,8603,8591,8797,7779,7778,7247,7246,7806,7245,7236,7237,7234,7235,7233,8242,7243,8240,8241,6860,6861,7238,7239,7588,7241,7584,7244,7586,7240,7580,7581,7582,7583,7188,7189,7578,7180,7181,7182,7183,7184,7185,7186,7187',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 21:14:20.048 +08.30.0.0000',
       scan=['128,132,135,138,141,144,116,119,148,122,125'],
       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_T13t32.J211420+083000_P15167v1_2017_11_10T19_34_13.099/working/VLASS1.1.sb34607978.eb34635523.58053.05546966435_split.ms'],
       imagename='VLASS1.1.ql.T13t32.J211420+083000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='6772,6773,6774,6775,6776,6777,6778,6779,6780,6781,6782,6783,6784,6785,6786,6787,6788,6789,6790,6791,6792,6793,6794,6795,6796,6797,6798,6799,6800,6801,6802,6831,6832,6833,6834,6835,6836,6837,6838,6839,6840,6841,6842,6843,6844,6845,6846,6847,6848,6849,6850,6851,6852,6853,6854,6855,6856,6857,6858,6859,6860,6861,7174,7175,7176,7177,7178,7179,7180,7181,7182,7183,7184,7185,7186,7187,7188,7189,7190,7191,7192,7193,7194,7195,7196,7197,7198,7199,7200,7201,7202,7203,7204,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,7263,7576,7577,7578,7579,7580,7581,7582,7583,7584,7585,7586,7587,7588,7589,7590,7591,7592,7593,7594,7595,7596,7597,7598,7599,7600,7601,7602,7603,7604,7605,7606,7777,7778,7779,7780,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,7836,7837,7838,7839,7840,7841,7842,7843,7844,7845,7846,7847,7848,7849,7850,7851,7852,7853,7854,7855,7856,7857,7858,7859,7860,7861,7862,7863,7864,7865,7866,8179,8180,8181,8182,8183,8184,8185,8186,8187,8188,8189,8190,8191,8192,8193,8194,8195,8196,8197,8198,8199,8200,8201,8202,8203,8204,8205,8206,8207,8208,8209,8238,8239,8240,8241,8242,8243,8244,8245,8246,8247,8248,8249,8250,8251,8252,8253,8254,8255,8256,8257,8258,8259,8260,8261,8262,8263,8264,8265,8266,8267,8268,8581,8582,8583,8584,8585,8586,8587,8588,8589,8590,8591,8592,8593,8594,8595,8596,8597,8598,8599,8600,8601,8602,8603,8604,8605,8606,8607,8608,8609,8610,8611,8782,8783,8784,8785,8786,8787,8788,8789,8790,8791,8792,8793,8794,8795,8796,8797,8798,8799,8800,8801,8802,8803,8804,8805,8806,8807,8808,8809,8810,8811,8812',
       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 21:14:20.048 +08.30.0.0000',
       scan=['128,132,135,138,141,144,116,119,148,122,125'],
       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_T13t32.J211420+083000_P15167v1_2017_11_10T19_34_13.099/working/VLASS1.1.sb34607978.eb34635523.58053.05546966435_split.ms'],
       imagename='VLASS1.1.ql.T13t32.J211420+083000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00052497550700882122, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='6772,6773,6774,6775,6776,6777,6778,6779,6780,6781,6782,6783,6784,6785,6786,6787,6788,6789,6790,6791,6792,6793,6794,6795,6796,6797,6798,6799,6800,6801,6802,6831,6832,6833,6834,6835,6836,6837,6838,6839,6840,6841,6842,6843,6844,6845,6846,6847,6848,6849,6850,6851,6852,6853,6854,6855,6856,6857,6858,6859,6860,6861,7174,7175,7176,7177,7178,7179,7180,7181,7182,7183,7184,7185,7186,7187,7188,7189,7190,7191,7192,7193,7194,7195,7196,7197,7198,7199,7200,7201,7202,7203,7204,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,7263,7576,7577,7578,7579,7580,7581,7582,7583,7584,7585,7586,7587,7588,7589,7590,7591,7592,7593,7594,7595,7596,7597,7598,7599,7600,7601,7602,7603,7604,7605,7606,7777,7778,7779,7780,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,7836,7837,7838,7839,7840,7841,7842,7843,7844,7845,7846,7847,7848,7849,7850,7851,7852,7853,7854,7855,7856,7857,7858,7859,7860,7861,7862,7863,7864,7865,7866,8179,8180,8181,8182,8183,8184,8185,8186,8187,8188,8189,8190,8191,8192,8193,8194,8195,8196,8197,8198,8199,8200,8201,8202,8203,8204,8205,8206,8207,8208,8209,8238,8239,8240,8241,8242,8243,8244,8245,8246,8247,8248,8249,8250,8251,8252,8253,8254,8255,8256,8257,8258,8259,8260,8261,8262,8263,8264,8265,8266,8267,8268,8581,8582,8583,8584,8585,8586,8587,8588,8589,8590,8591,8592,8593,8594,8595,8596,8597,8598,8599,8600,8601,8602,8603,8604,8605,8606,8607,8608,8609,8610,8611,8782,8783,8784,8785,8786,8787,8788,8789,8790,8791,8792,8793,8794,8795,8796,8797,8798,8799,8800,8801,8802,8803,8804,8805,8806,8807,8808,8809,8810,8811,8812',
       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.T13t32.J211420+083000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T13t32.J211420+083000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T13t32.J211420+083000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#