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

# hif_editimlist(parameter_file='parameter.12634.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_T19t20.J155205+323000_P12634v1_2018_01_29T22_43_48.799/working/VLASS1.1.sb34518430.eb34756479.58089.71818545139_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T19t20.J155205+323000_P12634v1_2018_01_29T22_43_48.799/rawdata/VLASS1.1.sb34518430.eb34756479.58089.71818545139.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='6718,8158,8159,8156,8157,8154,6719,8152,8153,8150,8087,6723,6722,6721,6720,6727,6726,6725,6724,8783,8782,7119,6728,8787,6714,8785,8784,8097,8096,7058,7059,8093,8092,8091,8090,7052,7122,7050,7051,7056,7057,7054,7055,8493,8492,8491,8490,7759,6710,8495,8494,7755,7754,7757,6711,7751,7750,7753,7752,8149,8148,8155,8496,8147,8146,8760,7128,7129,7126,7127,7124,6715,7049,7048,7120,7121,8498,7473,7737,8789,7736,8767,8788,8779,7735,7117,8586,7529,8500,8766,7139,7138,8515,8763,7131,7130,6703,7132,7135,6704,6707,6706,7118,8573,7765,8578,8761,8579,8786,7739,7738,7539,7538,7535,7534,7537,7536,7531,7530,7533,7532,8516,8517,8514,8168,8512,8513,8510,8511,7763,7758,8780,8769,8094,8762,7740,8499,6716,7145,7146,7147,7140,7141,7142,7143,7548,7549,7053,7540,7541,7542,7543,7544,7545,7546,7547,8501,8170,8503,8502,8505,8504,6717,8171,8509,8508,8507,8173,8099,8098,8488,7490,8506,7559,7558,7553,7552,7551,7550,7557,7556,7555,7554,8105,8104,8107,8106,8101,8100,8103,8102,8772,8773,8770,8771,8776,8777,8774,8775,7485,7484,7487,7486,7481,7480,7483,7482,7489,7488,8575,6712,6709,7756,8082,6708,8570,8083,8565,7764,7144,7125,8487,7760,7761,7762,7072,8581,8580,8583,8582,8576,8584,6713,8562,8079,8078,8172,8095,8174,8175,8176,8585,8781,8088,8778,8077,8089,8567,8566,6729,8564,8563,7073,8561,8560,8768,8577,7133,8569,8568,8765,7078,8764,7074,6705,7076,7077,7070,7071,7123,7134,7478,7479,8151,7741,7137,7470,7471,7472,7136,7474,7475,7476,7477,8163,8162,8161,8160,8167,8166,8165,8164,8169,7745,8572,8574,6730,6731,6732,6733,8556,8557,8571,8489,8790,8497,8558,8559,7748,8084,8085,8086,7749,8080,8081,7069,7068,7067,7066,7065,7064,7063,7062,7061,7060,7742,7743,7469,7468,7746,7747,7744,7075,7463,7462,7461,7460,7467,7466,7465,7464',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 15:52:5.524 +32.30.0.0000',
       scan=['130,133,136,139,146,142,114,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_T19t20.J155205+323000_P12634v1_2018_01_29T22_43_48.799/working/VLASS1.1.sb34518430.eb34756479.58089.71818545139_split.ms'],
       imagename='VLASS1.1.ql.T19t20.J155205+323000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='6703,6704,6705,6706,6707,6708,6709,6710,6711,6712,6713,6714,6715,6716,6717,6718,6719,6720,6721,6722,6723,6724,6725,6726,6727,6728,6729,6730,6731,6732,6733,7048,7049,7050,7051,7052,7053,7054,7055,7056,7057,7058,7059,7060,7061,7062,7063,7064,7065,7066,7067,7068,7069,7070,7071,7072,7073,7074,7075,7076,7077,7078,7117,7118,7119,7120,7121,7122,7123,7124,7125,7126,7127,7128,7129,7130,7131,7132,7133,7134,7135,7136,7137,7138,7139,7140,7141,7142,7143,7144,7145,7146,7147,7460,7461,7462,7463,7464,7465,7466,7467,7468,7469,7470,7471,7472,7473,7474,7475,7476,7477,7478,7479,7480,7481,7482,7483,7484,7485,7486,7487,7488,7489,7490,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,7558,7559,7735,7736,7737,7738,7739,7740,7741,7742,7743,7744,7745,7746,7747,7748,7749,7750,7751,7752,7753,7754,7755,7756,7757,7758,7759,7760,7761,7762,7763,7764,7765,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,8146,8147,8148,8149,8150,8151,8152,8153,8154,8155,8156,8157,8158,8159,8160,8161,8162,8163,8164,8165,8166,8167,8168,8169,8170,8171,8172,8173,8174,8175,8176,8487,8488,8489,8490,8491,8492,8493,8494,8495,8496,8497,8498,8499,8500,8501,8502,8503,8504,8505,8506,8507,8508,8509,8510,8511,8512,8513,8514,8515,8516,8517,8556,8557,8558,8559,8560,8561,8562,8563,8564,8565,8566,8567,8568,8569,8570,8571,8572,8573,8574,8575,8576,8577,8578,8579,8580,8581,8582,8583,8584,8585,8586,8760,8761,8762,8763,8764,8765,8766,8767,8768,8769,8770,8771,8772,8773,8774,8775,8776,8777,8778,8779,8780,8781,8782,8783,8784,8785,8786,8787,8788,8789,8790',
       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:52:5.524 +32.30.0.0000',
       scan=['130,133,136,139,146,142,114,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_T19t20.J155205+323000_P12634v1_2018_01_29T22_43_48.799/working/VLASS1.1.sb34518430.eb34756479.58089.71818545139_split.ms'],
       imagename='VLASS1.1.ql.T19t20.J155205+323000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00048444338304234403, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='6703,6704,6705,6706,6707,6708,6709,6710,6711,6712,6713,6714,6715,6716,6717,6718,6719,6720,6721,6722,6723,6724,6725,6726,6727,6728,6729,6730,6731,6732,6733,7048,7049,7050,7051,7052,7053,7054,7055,7056,7057,7058,7059,7060,7061,7062,7063,7064,7065,7066,7067,7068,7069,7070,7071,7072,7073,7074,7075,7076,7077,7078,7117,7118,7119,7120,7121,7122,7123,7124,7125,7126,7127,7128,7129,7130,7131,7132,7133,7134,7135,7136,7137,7138,7139,7140,7141,7142,7143,7144,7145,7146,7147,7460,7461,7462,7463,7464,7465,7466,7467,7468,7469,7470,7471,7472,7473,7474,7475,7476,7477,7478,7479,7480,7481,7482,7483,7484,7485,7486,7487,7488,7489,7490,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,7558,7559,7735,7736,7737,7738,7739,7740,7741,7742,7743,7744,7745,7746,7747,7748,7749,7750,7751,7752,7753,7754,7755,7756,7757,7758,7759,7760,7761,7762,7763,7764,7765,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,8146,8147,8148,8149,8150,8151,8152,8153,8154,8155,8156,8157,8158,8159,8160,8161,8162,8163,8164,8165,8166,8167,8168,8169,8170,8171,8172,8173,8174,8175,8176,8487,8488,8489,8490,8491,8492,8493,8494,8495,8496,8497,8498,8499,8500,8501,8502,8503,8504,8505,8506,8507,8508,8509,8510,8511,8512,8513,8514,8515,8516,8517,8556,8557,8558,8559,8560,8561,8562,8563,8564,8565,8566,8567,8568,8569,8570,8571,8572,8573,8574,8575,8576,8577,8578,8579,8580,8581,8582,8583,8584,8585,8586,8760,8761,8762,8763,8764,8765,8766,8767,8768,8769,8770,8771,8772,8773,8774,8775,8776,8777,8778,8779,8780,8781,8782,8783,8784,8785,8786,8787,8788,8789,8790',
       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.T19t20.J155205+323000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T19t20.J155205+323000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T19t20.J155205+323000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#