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

# hif_editimlist(parameter_file='parameter.13833.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_T23t16.J150259+483000_P13833v1_2018_01_18T22_44_28.357/working/VLASS1.1.sb34539903.eb34684984.58068.81223880787_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T23t16.J150259+483000_P13833v1_2018_01_18T22_44_28.357/rawdata/VLASS1.1.sb34539903.eb34684984.58068.81223880787.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='8880,8890,8507,9320,9295,8919,7302,7303,7300,7301,7306,7307,7304,7305,9285,9284,7308,7309,9281,9280,8099,8098,7694,7695,7696,7697,8497,8496,7692,7693,8905,8904,8499,8498,7698,7699,9276,9277,9278,9279,8920,9310,9286,9326,9296,8877,8876,8879,8878,8495,9328,7687,7686,7685,7684,7683,7682,7681,7680,9327,8928,7689,7688,9312,8104,9322,9308,9309,9323,8523,8522,8521,8520,8527,8526,8525,8524,9293,8888,8889,8884,8885,8886,8887,8107,8881,8882,8883,9324,8924,9325,8516,8517,8514,8515,8512,8513,8510,8511,9282,9294,8072,8095,8518,8519,9288,8094,9319,9318,9305,8093,9315,9314,8123,8122,8121,8092,7706,7707,7704,7705,7702,7703,7700,7701,8101,8896,8895,8090,8893,8892,7708,7709,9306,9287,8503,8502,8505,8504,8100,8506,8509,8508,9317,8118,8119,9283,9316,8112,8113,8110,8111,8116,8117,8114,8115,7711,7710,7713,7712,7715,7714,7717,7716,9313,8491,9290,8490,9311,8909,8120,8073,8908,7290,7291,7292,7293,7294,7295,7296,7297,7298,7299,8494,8109,8108,9291,8479,8478,8475,8474,8477,8476,8907,9297,8906,8899,8901,8898,8900,8897,9292,9289,8079,8078,8894,7287,7286,7285,7284,8071,8070,7289,7288,8075,8074,8077,8076,9304,8106,8891,7678,7679,7676,7677,7674,7675,7672,7673,7670,7671,8912,9307,8103,8913,9300,8102,8910,9301,8927,8926,8925,8911,8923,8922,8921,8493,8916,8917,9321,8914,8915,7669,7668,8492,7665,7664,7667,7666,8084,8085,8086,8087,8080,8081,8082,8083,8105,9298,9299,8088,8089,8480,8481,8482,8483,8484,8485,8486,8487,8488,8489,8918,8091',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 15:02:59.337 +48.30.0.0000',
       scan=['160,163,166,138,141,144,147,170,151,154,157'],
       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_T23t16.J150259+483000_P13833v1_2018_01_18T22_44_28.357/working/VLASS1.1.sb34539903.eb34684984.58068.81223880787_split.ms'],
       imagename='VLASS1.1.ql.T23t16.J150259+483000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='7284,7285,7286,7287,7288,7289,7290,7291,7292,7293,7294,7295,7296,7297,7298,7299,7300,7301,7302,7303,7304,7305,7306,7307,7308,7309,7664,7665,7666,7667,7668,7669,7670,7671,7672,7673,7674,7675,7676,7677,7678,7679,7680,7681,7682,7683,7684,7685,7686,7687,7688,7689,7692,7693,7694,7695,7696,7697,7698,7699,7700,7701,7702,7703,7704,7705,7706,7707,7708,7709,7710,7711,7712,7713,7714,7715,7716,7717,8070,8071,8072,8073,8074,8075,8076,8077,8078,8079,8080,8081,8082,8083,8084,8085,8086,8087,8088,8089,8090,8091,8092,8093,8094,8095,8098,8099,8100,8101,8102,8103,8104,8105,8106,8107,8108,8109,8110,8111,8112,8113,8114,8115,8116,8117,8118,8119,8120,8121,8122,8123,8474,8475,8476,8477,8478,8479,8480,8481,8482,8483,8484,8485,8486,8487,8488,8489,8490,8491,8492,8493,8494,8495,8496,8497,8498,8499,8502,8503,8504,8505,8506,8507,8508,8509,8510,8511,8512,8513,8514,8515,8516,8517,8518,8519,8520,8521,8522,8523,8524,8525,8526,8527,8876,8877,8878,8879,8880,8881,8882,8883,8884,8885,8886,8887,8888,8889,8890,8891,8892,8893,8894,8895,8896,8897,8898,8899,8900,8901,8904,8905,8906,8907,8908,8909,8910,8911,8912,8913,8914,8915,8916,8917,8918,8919,8920,8921,8922,8923,8924,8925,8926,8927,8928,9276,9277,9278,9279,9280,9281,9282,9283,9284,9285,9286,9287,9288,9289,9290,9291,9292,9293,9294,9295,9296,9297,9298,9299,9300,9301,9304,9305,9306,9307,9308,9309,9310,9311,9312,9313,9314,9315,9316,9317,9318,9319,9320,9321,9322,9323,9324,9325,9326,9327,9328',
       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:02:59.337 +48.30.0.0000',
       scan=['160,163,166,138,141,144,147,170,151,154,157'],
       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_T23t16.J150259+483000_P13833v1_2018_01_18T22_44_28.357/working/VLASS1.1.sb34539903.eb34684984.58068.81223880787_split.ms'],
       imagename='VLASS1.1.ql.T23t16.J150259+483000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00052360813007471733, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='7284,7285,7286,7287,7288,7289,7290,7291,7292,7293,7294,7295,7296,7297,7298,7299,7300,7301,7302,7303,7304,7305,7306,7307,7308,7309,7664,7665,7666,7667,7668,7669,7670,7671,7672,7673,7674,7675,7676,7677,7678,7679,7680,7681,7682,7683,7684,7685,7686,7687,7688,7689,7692,7693,7694,7695,7696,7697,7698,7699,7700,7701,7702,7703,7704,7705,7706,7707,7708,7709,7710,7711,7712,7713,7714,7715,7716,7717,8070,8071,8072,8073,8074,8075,8076,8077,8078,8079,8080,8081,8082,8083,8084,8085,8086,8087,8088,8089,8090,8091,8092,8093,8094,8095,8098,8099,8100,8101,8102,8103,8104,8105,8106,8107,8108,8109,8110,8111,8112,8113,8114,8115,8116,8117,8118,8119,8120,8121,8122,8123,8474,8475,8476,8477,8478,8479,8480,8481,8482,8483,8484,8485,8486,8487,8488,8489,8490,8491,8492,8493,8494,8495,8496,8497,8498,8499,8502,8503,8504,8505,8506,8507,8508,8509,8510,8511,8512,8513,8514,8515,8516,8517,8518,8519,8520,8521,8522,8523,8524,8525,8526,8527,8876,8877,8878,8879,8880,8881,8882,8883,8884,8885,8886,8887,8888,8889,8890,8891,8892,8893,8894,8895,8896,8897,8898,8899,8900,8901,8904,8905,8906,8907,8908,8909,8910,8911,8912,8913,8914,8915,8916,8917,8918,8919,8920,8921,8922,8923,8924,8925,8926,8927,8928,9276,9277,9278,9279,9280,9281,9282,9283,9284,9285,9286,9287,9288,9289,9290,9291,9292,9293,9294,9295,9296,9297,9298,9299,9300,9301,9304,9305,9306,9307,9308,9309,9310,9311,9312,9313,9314,9315,9316,9317,9318,9319,9320,9321,9322,9323,9324,9325,9326,9327,9328',
       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.T23t16.J150259+483000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T23t16.J150259+483000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T23t16.J150259+483000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#