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

# hif_editimlist(parameter_file='parameter.5662.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_T28t05.J094602+693000_P5662v1_2017_09_26T22_41_52.086/working/VLASS1.1.sb34417859.eb34434020.58014.453630601856_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T28t05.J094602+693000_P5662v1_2017_09_26T22_41_52.086/rawdata/VLASS1.1.sb34417859.eb34434020.58014.453630601856.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='4446,4447,4448,4449,6099,6098,6152,6153,6154,6155,6156,6091,6090,6093,6092,6095,6094,6097,6096,5300,5301,5302,5670,5671,5672,5673,5674,5675,4813,4812,4811,6508,4817,4816,4815,4814,4819,4818,6507,6506,6077,6076,6075,6074,6079,6078,5708,5709,5704,5705,5706,5707,5700,5701,5702,5703,5229,5228,5223,5222,5221,5227,5226,5225,5224,6104,6103,6102,6101,6100,4468,4466,4467,4464,4465,4462,4463,4460,4461,6150,6151,6572,6573,6570,6571,6576,6577,6574,6575,6578,6579,4868,4869,4862,4863,4860,4861,4866,4867,4864,4865,6569,6568,5663,5662,4437,5660,5696,5697,5667,5698,5699,6533,5665,6561,5664,6560,6143,6142,6141,6140,6147,6146,6145,6144,6149,6148,5250,5251,6583,6582,6581,6580,6587,6586,6585,6584,6589,6588,6536,6537,5661,6535,6532,5666,6530,6531,5669,5668,4439,4438,4828,4829,4826,4827,4824,4825,4822,4823,4820,4821,5719,5718,5717,5716,5715,5714,5713,5712,5711,5710,6527,6526,6529,6528,6138,6139,6132,6133,6130,6131,6136,6137,6134,6135,4880,4881,4882,4883,4875,4874,4877,4876,4871,4870,4873,4872,6565,6564,6567,6566,4879,4878,6563,6562,6129,6128,5249,5248,5245,5244,5247,5246,5241,5240,5243,5242,6509,6521,6520,6523,4810,6525,6524,5658,5659,5656,5657,5654,5655,5652,5653,5650,5651,4453,4452,4831,4830,4833,4832,6534,5298,5299,5296,5297,5294,5295,5292,5293,5290,5291,5722,5723,5720,5721,5726,6559,5724,5725,4858,6125,6127,6126,6088,6089,6086,6087,6084,6085,6082,6083,6080,6081,5274,5275,5276,5277,5272,5273,5278,5279,6522,6518,6519,4802,4803,4804,4805,4806,4807,4808,4809,6512,6513,6514,6515,6516,6517,5238,5239,5230,5231,5232,5233,5234,5235,5236,5237,5649,5648,4451,4450,4457,4456,4455,4454,4459,4458,5645,5647,5646,5289,5288,5281,5280,5283,5282,5285,5284,5287,5286,4856,4440,4441,4442,6510,4443,4859,6511,4857,4444,4855,4854,4853,4852,4445',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 09:46:2.176 +69.30.0.0000',
       scan=['97,100,103,107,110,81,84,113,88,91,94'],
       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_T28t05.J094602+693000_P5662v1_2017_09_26T22_41_52.086/working/VLASS1.1.sb34417859.eb34434020.58014.453630601856_split.ms'],
       imagename='VLASS1.1.ql.T28t05.J094602+693000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='4437,4438,4439,4440,4441,4442,4443,4444,4445,4446,4447,4448,4449,4450,4451,4452,4453,4454,4455,4456,4457,4458,4459,4460,4461,4462,4463,4464,4465,4466,4467,4468,4802,4803,4804,4805,4806,4807,4808,4809,4810,4811,4812,4813,4814,4815,4816,4817,4818,4819,4820,4821,4822,4823,4824,4825,4826,4827,4828,4829,4830,4831,4832,4833,4852,4853,4854,4855,4856,4857,4858,4859,4860,4861,4862,4863,4864,4865,4866,4867,4868,4869,4870,4871,4872,4873,4874,4875,4876,4877,4878,4879,4880,4881,4882,4883,5221,5222,5223,5224,5225,5226,5227,5228,5229,5230,5231,5232,5233,5234,5235,5236,5237,5238,5239,5240,5241,5242,5243,5244,5245,5246,5247,5248,5249,5250,5251,5272,5273,5274,5275,5276,5277,5278,5279,5280,5281,5282,5283,5284,5285,5286,5287,5288,5289,5290,5291,5292,5293,5294,5295,5296,5297,5298,5299,5300,5301,5302,5645,5646,5647,5648,5649,5650,5651,5652,5653,5654,5655,5656,5657,5658,5659,5660,5661,5662,5663,5664,5665,5666,5667,5668,5669,5670,5671,5672,5673,5674,5675,5696,5697,5698,5699,5700,5701,5702,5703,5704,5705,5706,5707,5708,5709,5710,5711,5712,5713,5714,5715,5716,5717,5718,5719,5720,5721,5722,5723,5724,5725,5726,6074,6075,6076,6077,6078,6079,6080,6081,6082,6083,6084,6085,6086,6087,6088,6089,6090,6091,6092,6093,6094,6095,6096,6097,6098,6099,6100,6101,6102,6103,6104,6125,6126,6127,6128,6129,6130,6131,6132,6133,6134,6135,6136,6137,6138,6139,6140,6141,6142,6143,6144,6145,6146,6147,6148,6149,6150,6151,6152,6153,6154,6155,6156,6506,6507,6508,6509,6510,6511,6512,6513,6514,6515,6516,6517,6518,6519,6520,6521,6522,6523,6524,6525,6526,6527,6528,6529,6530,6531,6532,6533,6534,6535,6536,6537,6559,6560,6561,6562,6563,6564,6565,6566,6567,6568,6569,6570,6571,6572,6573,6574,6575,6576,6577,6578,6579,6580,6581,6582,6583,6584,6585,6586,6587,6588,6589',
       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 09:46:2.176 +69.30.0.0000',
       scan=['97,100,103,107,110,81,84,113,88,91,94'],
       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_T28t05.J094602+693000_P5662v1_2017_09_26T22_41_52.086/working/VLASS1.1.sb34417859.eb34434020.58014.453630601856_split.ms'],
       imagename='VLASS1.1.ql.T28t05.J094602+693000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00066763238494713724, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='4437,4438,4439,4440,4441,4442,4443,4444,4445,4446,4447,4448,4449,4450,4451,4452,4453,4454,4455,4456,4457,4458,4459,4460,4461,4462,4463,4464,4465,4466,4467,4468,4802,4803,4804,4805,4806,4807,4808,4809,4810,4811,4812,4813,4814,4815,4816,4817,4818,4819,4820,4821,4822,4823,4824,4825,4826,4827,4828,4829,4830,4831,4832,4833,4852,4853,4854,4855,4856,4857,4858,4859,4860,4861,4862,4863,4864,4865,4866,4867,4868,4869,4870,4871,4872,4873,4874,4875,4876,4877,4878,4879,4880,4881,4882,4883,5221,5222,5223,5224,5225,5226,5227,5228,5229,5230,5231,5232,5233,5234,5235,5236,5237,5238,5239,5240,5241,5242,5243,5244,5245,5246,5247,5248,5249,5250,5251,5272,5273,5274,5275,5276,5277,5278,5279,5280,5281,5282,5283,5284,5285,5286,5287,5288,5289,5290,5291,5292,5293,5294,5295,5296,5297,5298,5299,5300,5301,5302,5645,5646,5647,5648,5649,5650,5651,5652,5653,5654,5655,5656,5657,5658,5659,5660,5661,5662,5663,5664,5665,5666,5667,5668,5669,5670,5671,5672,5673,5674,5675,5696,5697,5698,5699,5700,5701,5702,5703,5704,5705,5706,5707,5708,5709,5710,5711,5712,5713,5714,5715,5716,5717,5718,5719,5720,5721,5722,5723,5724,5725,5726,6074,6075,6076,6077,6078,6079,6080,6081,6082,6083,6084,6085,6086,6087,6088,6089,6090,6091,6092,6093,6094,6095,6096,6097,6098,6099,6100,6101,6102,6103,6104,6125,6126,6127,6128,6129,6130,6131,6132,6133,6134,6135,6136,6137,6138,6139,6140,6141,6142,6143,6144,6145,6146,6147,6148,6149,6150,6151,6152,6153,6154,6155,6156,6506,6507,6508,6509,6510,6511,6512,6513,6514,6515,6516,6517,6518,6519,6520,6521,6522,6523,6524,6525,6526,6527,6528,6529,6530,6531,6532,6533,6534,6535,6536,6537,6559,6560,6561,6562,6563,6564,6565,6566,6567,6568,6569,6570,6571,6572,6573,6574,6575,6576,6577,6578,6579,6580,6581,6582,6583,6584,6585,6586,6587,6588,6589',
       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.T28t05.J094602+693000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T28t05.J094602+693000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T28t05.J094602+693000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#