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

# hif_editimlist(parameter_file='parameter.11861.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_T18t21.J152454+283000_P11861v1_2017_11_19T19_28_57.739/working/VLASS1.1.sb34516171.eb34560849.58028.92610690973_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T18t21.J152454+283000_P11861v1_2017_11_19T19_28_57.739/rawdata/VLASS1.1.sb34516171.eb34560849.58028.92610690973.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='8307,8306,8305,8304,8303,8302,8301,8300,8309,8308,7302,7303,7300,7301,7306,7307,7304,7305,7308,7309,8099,8098,8909,8908,8905,8904,8907,8906,8901,8900,8903,8902,7889,7888,7881,7880,7883,7882,7885,7884,8730,7886,8097,8096,8095,8094,7900,7901,7902,7903,7904,7905,7906,8092,7706,7707,7704,7705,7702,7703,7700,7701,8090,7708,7709,6891,6890,6893,6892,6895,6894,6897,6896,6898,7078,7079,8728,7074,7075,7076,7077,7071,7072,7073,8493,8492,8491,7697,8497,8496,8495,8494,8499,8498,7698,7699,8330,6716,6717,6714,6715,6712,6713,6710,6711,7319,7318,7315,7314,7317,7316,7311,7310,7313,7312,8879,7713,7712,8729,7492,7493,8721,8720,8723,8722,8725,8724,8727,8726,8501,8500,8503,8502,8505,8504,8507,8506,8509,8508,8299,7711,7710,7498,7499,7715,7714,7717,7716,7719,7718,7490,7491,7496,7497,7494,7495,8479,8325,6700,8327,7324,8326,7320,7321,7322,7323,6705,8320,8323,8322,7100,7101,6868,6869,8710,8711,8712,8713,8714,8715,8716,8717,8718,8719,7887,7720,7721,7722,7723,7724,7725,7726,7727,6709,6708,8329,8328,6701,8324,6703,6702,8321,6704,6707,6706,8888,8889,8884,8885,8886,8887,8880,8881,8882,8883,7096,7097,7094,7095,7092,7093,7090,7091,7098,7099,6688,6689,8490,6687,8105,8104,8107,8106,8101,7295,7296,8102,7298,7299,7485,7484,7487,7486,7481,7480,7483,7482,7489,7488,7478,7479,7474,7475,7476,7477,8314,8315,8316,8317,8310,8311,8312,8313,8318,8319,8084,8085,8086,8087,8080,8081,8082,8083,8088,8089,7504,7500,7501,7502,7503,6879,6878,6877,6876,6875,6874,6873,6872,6871,6870,8703,8702,8701,8700,8707,8706,7898,7899,7896,7897,8709,7895,7892,7893,7890,7891,7876,7877,7878,7879,8705,8704,7894,8708,8899,8898,8897,8896,8895,8894,8893,8892,8891,8890,6886,6887,6884,6885,6882,6883,6880,6881,6888,6889,7081,7080,7083,7082,7085,7084,7087,7086,7089,7088,6699,6698,6693,6692,6691,6690,6697,6696,6695,6694,8079,8078,8077,7294,8100,8103,7297,8093,8480,8481,8482,8483,8484,8485,8486,8487,8488,8489,8091',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 15:24:54.992 +28.30.0.0000',
       scan=['129,133,136,139,109,142,112,145,117,120,123,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_T18t21.J152454+283000_P11861v1_2017_11_19T19_28_57.739/working/VLASS1.1.sb34516171.eb34560849.58028.92610690973_split.ms'],
       imagename='VLASS1.1.ql.T18t21.J152454+283000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='6687,6688,6689,6690,6691,6692,6693,6694,6695,6696,6697,6698,6699,6700,6701,6702,6703,6704,6705,6706,6707,6708,6709,6710,6711,6712,6713,6714,6715,6716,6717,6868,6869,6870,6871,6872,6873,6874,6875,6876,6877,6878,6879,6880,6881,6882,6883,6884,6885,6886,6887,6888,6889,6890,6891,6892,6893,6894,6895,6896,6897,6898,7071,7072,7073,7074,7075,7076,7077,7078,7079,7080,7081,7082,7083,7084,7085,7086,7087,7088,7089,7090,7091,7092,7093,7094,7095,7096,7097,7098,7099,7100,7101,7294,7295,7296,7297,7298,7299,7300,7301,7302,7303,7304,7305,7306,7307,7308,7309,7310,7311,7312,7313,7314,7315,7316,7317,7318,7319,7320,7321,7322,7323,7324,7474,7475,7476,7477,7478,7479,7480,7481,7482,7483,7484,7485,7486,7487,7488,7489,7490,7491,7492,7493,7494,7495,7496,7497,7498,7499,7500,7501,7502,7503,7504,7697,7698,7699,7700,7701,7702,7703,7704,7705,7706,7707,7708,7709,7710,7711,7712,7713,7714,7715,7716,7717,7718,7719,7720,7721,7722,7723,7724,7725,7726,7727,7876,7877,7878,7879,7880,7881,7882,7883,7884,7885,7886,7887,7888,7889,7890,7891,7892,7893,7894,7895,7896,7897,7898,7899,7900,7901,7902,7903,7904,7905,7906,8077,8078,8079,8080,8081,8082,8083,8084,8085,8086,8087,8088,8089,8090,8091,8092,8093,8094,8095,8096,8097,8098,8099,8100,8101,8102,8103,8104,8105,8106,8107,8299,8300,8301,8302,8303,8304,8305,8306,8307,8308,8309,8310,8311,8312,8313,8314,8315,8316,8317,8318,8319,8320,8321,8322,8323,8324,8325,8326,8327,8328,8329,8330,8479,8480,8481,8482,8483,8484,8485,8486,8487,8488,8489,8490,8491,8492,8493,8494,8495,8496,8497,8498,8499,8500,8501,8502,8503,8504,8505,8506,8507,8508,8509,8700,8701,8702,8703,8704,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,8879,8880,8881,8882,8883,8884,8885,8886,8887,8888,8889,8890,8891,8892,8893,8894,8895,8896,8897,8898,8899,8900,8901,8902,8903,8904,8905,8906,8907,8908,8909',
       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 15:24:54.992 +28.30.0.0000',
       scan=['129,133,136,139,109,142,112,145,117,120,123,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_T18t21.J152454+283000_P11861v1_2017_11_19T19_28_57.739/working/VLASS1.1.sb34516171.eb34560849.58028.92610690973_split.ms'],
       imagename='VLASS1.1.ql.T18t21.J152454+283000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00055380240594713533, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='6687,6688,6689,6690,6691,6692,6693,6694,6695,6696,6697,6698,6699,6700,6701,6702,6703,6704,6705,6706,6707,6708,6709,6710,6711,6712,6713,6714,6715,6716,6717,6868,6869,6870,6871,6872,6873,6874,6875,6876,6877,6878,6879,6880,6881,6882,6883,6884,6885,6886,6887,6888,6889,6890,6891,6892,6893,6894,6895,6896,6897,6898,7071,7072,7073,7074,7075,7076,7077,7078,7079,7080,7081,7082,7083,7084,7085,7086,7087,7088,7089,7090,7091,7092,7093,7094,7095,7096,7097,7098,7099,7100,7101,7294,7295,7296,7297,7298,7299,7300,7301,7302,7303,7304,7305,7306,7307,7308,7309,7310,7311,7312,7313,7314,7315,7316,7317,7318,7319,7320,7321,7322,7323,7324,7474,7475,7476,7477,7478,7479,7480,7481,7482,7483,7484,7485,7486,7487,7488,7489,7490,7491,7492,7493,7494,7495,7496,7497,7498,7499,7500,7501,7502,7503,7504,7697,7698,7699,7700,7701,7702,7703,7704,7705,7706,7707,7708,7709,7710,7711,7712,7713,7714,7715,7716,7717,7718,7719,7720,7721,7722,7723,7724,7725,7726,7727,7876,7877,7878,7879,7880,7881,7882,7883,7884,7885,7886,7887,7888,7889,7890,7891,7892,7893,7894,7895,7896,7897,7898,7899,7900,7901,7902,7903,7904,7905,7906,8077,8078,8079,8080,8081,8082,8083,8084,8085,8086,8087,8088,8089,8090,8091,8092,8093,8094,8095,8096,8097,8098,8099,8100,8101,8102,8103,8104,8105,8106,8107,8299,8300,8301,8302,8303,8304,8305,8306,8307,8308,8309,8310,8311,8312,8313,8314,8315,8316,8317,8318,8319,8320,8321,8322,8323,8324,8325,8326,8327,8328,8329,8330,8479,8480,8481,8482,8483,8484,8485,8486,8487,8488,8489,8490,8491,8492,8493,8494,8495,8496,8497,8498,8499,8500,8501,8502,8503,8504,8505,8506,8507,8508,8509,8700,8701,8702,8703,8704,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,8879,8880,8881,8882,8883,8884,8885,8886,8887,8888,8889,8890,8891,8892,8893,8894,8895,8896,8897,8898,8899,8900,8901,8902,8903,8904,8905,8906,8907,8908,8909',
       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.T18t21.J152454+283000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T18t21.J152454+283000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T18t21.J152454+283000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#