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

# hif_editimlist(parameter_file='parameter.10096.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_T24t09.J080944+523000_P10096v1_2017_10_10T20_38_32.637/working/VLASS1.1.sb34499142.eb34512738.58022.410128460644_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T24t09.J080944+523000_P10096v1_2017_10_10T20_38_32.637/rawdata/VLASS1.1.sb34499142.eb34512738.58022.410128460644.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='7077,5467,6694,5466,7071,6125,6124,6127,6126,5472,5473,5470,5471,5749,7119,7073,6129,6128,7059,5821,5820,5823,5822,5825,5824,5827,5826,5829,5828,6674,6740,7056,7057,7054,7055,6327,6326,6325,6324,6323,6322,6321,6320,7121,6329,6328,7125,7128,6382,7072,5461,5460,5751,5462,5465,5464,5755,5754,5469,5468,5759,5758,7049,7048,7120,6687,7113,7065,7112,6392,6393,6385,6391,7111,6394,6395,6312,6313,6310,6311,6316,6317,6314,6315,6121,6318,6319,6390,6120,6123,6122,6757,6143,6142,6141,6140,6147,6146,5768,5769,5766,5767,5764,5765,5762,5763,5760,5761,7050,7118,7063,6137,6698,6309,6308,5840,6305,6304,6307,6306,6303,7130,7114,7070,7067,6748,6134,7123,7058,5779,5778,6677,6676,6702,6692,5771,5770,5773,5772,5775,5774,5777,5776,6675,6381,6690,7133,6135,6697,7062,7140,7141,6696,6378,6379,6705,6695,6766,6370,6371,6372,6373,6374,6375,6376,6377,6380,7051,7132,7061,7138,7117,6762,6704,7116,7127,7126,6760,6383,6753,7135,5752,7060,6767,6679,6765,6764,6763,6369,6761,6678,6688,6689,6139,6368,7075,6680,6681,6682,6683,6367,6366,6132,6364,7134,7069,5837,7115,5834,6693,5835,7124,6742,6744,6138,6758,6759,6118,6119,6384,6752,6684,6750,6751,6756,6387,6754,6685,5818,5819,6699,6386,6365,7066,5810,5811,5812,5813,5814,5815,5816,5817,7064,6145,6136,6144,6691,5446,6745,6701,6747,6746,6741,5838,6743,6700,5839,6749,6703,5458,5459,6148,7078,7079,5450,5451,5452,5453,5454,5455,5456,5457,7053,7137,7139,7136,7068,7131,6755,6686,6389,7076,6736,6737,6738,6739,7122,5836,5753,5449,5448,5832,5833,5830,5831,5443,6133,6130,6131,5447,5463,5445,5444,6388,7074,5750,6330,6331,6332,6333,5757,7129,7052,5756',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 08:09:44.737 +52.30.0.0000',
       scan=['128,131,134,106,109,112,115,119,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_T24t09.J080944+523000_P10096v1_2017_10_10T20_38_32.637/working/VLASS1.1.sb34499142.eb34512738.58022.410128460644_split.ms'],
       imagename='VLASS1.1.ql.T24t09.J080944+523000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='5443,5444,5445,5446,5447,5448,5449,5450,5451,5452,5453,5454,5455,5456,5457,5458,5459,5460,5461,5462,5463,5464,5465,5466,5467,5468,5469,5470,5471,5472,5473,5749,5750,5751,5752,5753,5754,5755,5756,5757,5758,5759,5760,5761,5762,5763,5764,5765,5766,5767,5768,5769,5770,5771,5772,5773,5774,5775,5776,5777,5778,5779,5810,5811,5812,5813,5814,5815,5816,5817,5818,5819,5820,5821,5822,5823,5824,5825,5826,5827,5828,5829,5830,5831,5832,5833,5834,5835,5836,5837,5838,5839,5840,6118,6119,6120,6121,6122,6123,6124,6125,6126,6127,6128,6129,6130,6131,6132,6133,6134,6135,6136,6137,6138,6139,6140,6141,6142,6143,6144,6145,6146,6147,6148,6303,6304,6305,6306,6307,6308,6309,6310,6311,6312,6313,6314,6315,6316,6317,6318,6319,6320,6321,6322,6323,6324,6325,6326,6327,6328,6329,6330,6331,6332,6333,6364,6365,6366,6367,6368,6369,6370,6371,6372,6373,6374,6375,6376,6377,6378,6379,6380,6381,6382,6383,6384,6385,6386,6387,6388,6389,6390,6391,6392,6393,6394,6395,6674,6675,6676,6677,6678,6679,6680,6681,6682,6683,6684,6685,6686,6687,6688,6689,6690,6691,6692,6693,6694,6695,6696,6697,6698,6699,6700,6701,6702,6703,6704,6705,6736,6737,6738,6739,6740,6741,6742,6743,6744,6745,6746,6747,6748,6749,6750,6751,6752,6753,6754,6755,6756,6757,6758,6759,6760,6761,6762,6763,6764,6765,6766,6767,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,7079,7111,7112,7113,7114,7115,7116,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',
       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 08:09:44.737 +52.30.0.0000',
       scan=['128,131,134,106,109,112,115,119,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_T24t09.J080944+523000_P10096v1_2017_10_10T20_38_32.637/working/VLASS1.1.sb34499142.eb34512738.58022.410128460644_split.ms'],
       imagename='VLASS1.1.ql.T24t09.J080944+523000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00051338888602949731, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='5443,5444,5445,5446,5447,5448,5449,5450,5451,5452,5453,5454,5455,5456,5457,5458,5459,5460,5461,5462,5463,5464,5465,5466,5467,5468,5469,5470,5471,5472,5473,5749,5750,5751,5752,5753,5754,5755,5756,5757,5758,5759,5760,5761,5762,5763,5764,5765,5766,5767,5768,5769,5770,5771,5772,5773,5774,5775,5776,5777,5778,5779,5810,5811,5812,5813,5814,5815,5816,5817,5818,5819,5820,5821,5822,5823,5824,5825,5826,5827,5828,5829,5830,5831,5832,5833,5834,5835,5836,5837,5838,5839,5840,6118,6119,6120,6121,6122,6123,6124,6125,6126,6127,6128,6129,6130,6131,6132,6133,6134,6135,6136,6137,6138,6139,6140,6141,6142,6143,6144,6145,6146,6147,6148,6303,6304,6305,6306,6307,6308,6309,6310,6311,6312,6313,6314,6315,6316,6317,6318,6319,6320,6321,6322,6323,6324,6325,6326,6327,6328,6329,6330,6331,6332,6333,6364,6365,6366,6367,6368,6369,6370,6371,6372,6373,6374,6375,6376,6377,6378,6379,6380,6381,6382,6383,6384,6385,6386,6387,6388,6389,6390,6391,6392,6393,6394,6395,6674,6675,6676,6677,6678,6679,6680,6681,6682,6683,6684,6685,6686,6687,6688,6689,6690,6691,6692,6693,6694,6695,6696,6697,6698,6699,6700,6701,6702,6703,6704,6705,6736,6737,6738,6739,6740,6741,6742,6743,6744,6745,6746,6747,6748,6749,6750,6751,6752,6753,6754,6755,6756,6757,6758,6759,6760,6761,6762,6763,6764,6765,6766,6767,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,7079,7111,7112,7113,7114,7115,7116,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',
       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.T24t09.J080944+523000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T24t09.J080944+523000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T24t09.J080944+523000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#