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

# hif_editimlist(parameter_file='parameter.list')
#

# hif_transformimagedata(datacolumn='corrected', modify_weights=False, clear_pointing=True)
#
mstransform(timebin='0s',
            outputvis='VLASS1.1.sb34739674.eb34813174.58124.37509868055_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T10t16.J100159-003000_P19398v1_2018_04_10T23_01_45.486/rawdata/VLASS1.1.sb34739674.eb34813174.58124.37509868055.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='7300,6513,6914,6881,6921,4877,6120,6088,6089,6086,6087,6084,6085,6082,6083,6080,6081,5274,5275,5276,5277,5270,5271,5272,5273,4901,4900,7308,4902,5278,5279,6475,6474,4868,4869,6915,6888,4862,4863,4860,4861,4866,4867,4864,4865,4888,4889,7310,4874,4880,4881,4882,4883,4884,4885,4886,4887,6099,6098,6916,5303,6910,6911,6912,6913,6091,4871,6092,6095,6919,6097,4870,5267,5266,5265,5264,5263,5262,5261,4872,5300,5301,5302,5260,4879,4878,5269,5268,5692,5304,5690,5691,5696,5697,5694,5695,6929,5698,5699,6116,6928,4899,4898,6879,4893,4892,4891,4890,4897,4896,4895,4894,6907,6904,6903,5306,6901,6900,5307,6927,6909,6908,6518,6519,6889,6926,5256,5257,5255,6510,6511,5258,5259,6514,6515,6516,6517,6478,6067,5689,6931,6493,5685,5684,5686,5681,5680,5683,5682,7302,7303,6479,5670,5671,5672,5673,5674,5675,5676,5677,5678,5679,4873,7307,6932,6509,6090,6930,7304,6503,6502,6501,7305,6507,6506,6505,6504,6077,6076,6075,6074,6073,6072,6071,6070,6079,6078,7297,7309,7286,5308,6492,6115,5663,5662,5661,6496,5667,5666,5665,5664,6925,6924,5669,5668,6884,6920,6923,6922,5708,5709,6512,5704,5705,5706,5707,5700,5701,5702,5703,6887,6885,7289,7285,6882,7290,7291,7292,7288,7294,7295,7296,6473,7298,7299,6117,6498,6521,6520,6488,6489,6525,6524,6880,6526,6482,6483,6480,6481,6486,6487,6484,6485,6917,7301,6118,6119,7306,6114,6110,5305,5714,5713,5712,5711,5710,6918,6883,7293,6902,6495,6494,6497,6068,6491,6490,5298,5299,5296,5297,5294,5295,5292,5293,5290,5291,6891,6890,6893,6892,6895,6894,6897,6896,6899,6898,6069,6105,6112,6109,6108,6107,6106,6096,6104,6103,6102,6101,6100,4849,7287,6111,6477,5289,5288,6113,5280,5283,5285,5284,5287,5286,5693,6476,6886,6508,6523,6522,4859,4858,4857,4856,4855,4854,4853,4852,4851,4850',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 10:01:60.000 -00.30.0.0000',
       scan=['97,101,104,107,110,113,82,116,85,88,120,91,94'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34739674.eb34813174.58124.37509868055_split.ms'],
       imagename='VLASS1.1.ql.T10t16.J100159-003000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='4849,4850,4851,4852,4853,4854,4855,4856,4857,4858,4859,4860,4861,4862,4863,4864,4865,4866,4867,4868,4869,4870,4871,4872,4873,4874,4877,4878,4879,4880,4881,4882,4883,4884,4885,4886,4887,4888,4889,4890,4891,4892,4893,4894,4895,4896,4897,4898,4899,4900,4901,4902,5255,5256,5257,5258,5259,5260,5261,5262,5263,5264,5265,5266,5267,5268,5269,5270,5271,5272,5273,5274,5275,5276,5277,5278,5279,5280,5283,5284,5285,5286,5287,5288,5289,5290,5291,5292,5293,5294,5295,5296,5297,5298,5299,5300,5301,5302,5303,5304,5305,5306,5307,5308,5661,5662,5663,5664,5665,5666,5667,5668,5669,5670,5671,5672,5673,5674,5675,5676,5677,5678,5679,5680,5681,5682,5683,5684,5685,5686,5689,5690,5691,5692,5693,5694,5695,5696,5697,5698,5699,5700,5701,5702,5703,5704,5705,5706,5707,5708,5709,5710,5711,5712,5713,5714,6067,6068,6069,6070,6071,6072,6073,6074,6075,6076,6077,6078,6079,6080,6081,6082,6083,6084,6085,6086,6087,6088,6089,6090,6091,6092,6095,6096,6097,6098,6099,6100,6101,6102,6103,6104,6105,6106,6107,6108,6109,6110,6111,6112,6113,6114,6115,6116,6117,6118,6119,6120,6473,6474,6475,6476,6477,6478,6479,6480,6481,6482,6483,6484,6485,6486,6487,6488,6489,6490,6491,6492,6493,6494,6495,6496,6497,6498,6501,6502,6503,6504,6505,6506,6507,6508,6509,6510,6511,6512,6513,6514,6515,6516,6517,6518,6519,6520,6521,6522,6523,6524,6525,6526,6879,6880,6881,6882,6883,6884,6885,6886,6887,6888,6889,6890,6891,6892,6893,6894,6895,6896,6897,6898,6899,6900,6901,6902,6903,6904,6907,6908,6909,6910,6911,6912,6913,6914,6915,6916,6917,6918,6919,6920,6921,6922,6923,6924,6925,6926,6927,6928,6929,6930,6931,6932,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,7310',
       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 10:01:60.000 -00.30.0.0000',
       scan=['97,101,104,107,110,113,82,116,85,88,120,91,94'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34739674.eb34813174.58124.37509868055_split.ms'],
       imagename='VLASS1.1.ql.T10t16.J100159-003000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00069449914560318576, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='4849,4850,4851,4852,4853,4854,4855,4856,4857,4858,4859,4860,4861,4862,4863,4864,4865,4866,4867,4868,4869,4870,4871,4872,4873,4874,4877,4878,4879,4880,4881,4882,4883,4884,4885,4886,4887,4888,4889,4890,4891,4892,4893,4894,4895,4896,4897,4898,4899,4900,4901,4902,5255,5256,5257,5258,5259,5260,5261,5262,5263,5264,5265,5266,5267,5268,5269,5270,5271,5272,5273,5274,5275,5276,5277,5278,5279,5280,5283,5284,5285,5286,5287,5288,5289,5290,5291,5292,5293,5294,5295,5296,5297,5298,5299,5300,5301,5302,5303,5304,5305,5306,5307,5308,5661,5662,5663,5664,5665,5666,5667,5668,5669,5670,5671,5672,5673,5674,5675,5676,5677,5678,5679,5680,5681,5682,5683,5684,5685,5686,5689,5690,5691,5692,5693,5694,5695,5696,5697,5698,5699,5700,5701,5702,5703,5704,5705,5706,5707,5708,5709,5710,5711,5712,5713,5714,6067,6068,6069,6070,6071,6072,6073,6074,6075,6076,6077,6078,6079,6080,6081,6082,6083,6084,6085,6086,6087,6088,6089,6090,6091,6092,6095,6096,6097,6098,6099,6100,6101,6102,6103,6104,6105,6106,6107,6108,6109,6110,6111,6112,6113,6114,6115,6116,6117,6118,6119,6120,6473,6474,6475,6476,6477,6478,6479,6480,6481,6482,6483,6484,6485,6486,6487,6488,6489,6490,6491,6492,6493,6494,6495,6496,6497,6498,6501,6502,6503,6504,6505,6506,6507,6508,6509,6510,6511,6512,6513,6514,6515,6516,6517,6518,6519,6520,6521,6522,6523,6524,6525,6526,6879,6880,6881,6882,6883,6884,6885,6886,6887,6888,6889,6890,6891,6892,6893,6894,6895,6896,6897,6898,6899,6900,6901,6902,6903,6904,6907,6908,6909,6910,6911,6912,6913,6914,6915,6916,6917,6918,6919,6920,6921,6922,6923,6924,6925,6926,6927,6928,6929,6930,6931,6932,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,7310',
       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.T10t16.J100159-003000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T10t16.J100159-003000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T10t16.J100159-003000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#