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

# hif_editimlist(parameter_file='parameter.17025.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_T32t01.J041225+853000_P17025v1_2018_01_18T23_06_05.509/working/VLASS1.1.sb34666958.eb34700016.58074.035491539355_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T32t01.J041225+853000_P17025v1_2018_01_18T23_06_05.509/rawdata/VLASS1.1.sb34666958.eb34700016.58074.035491539355.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='5821,5820,5823,5822,4149,3999,3998,3997,3996,3995,3994,3993,3992,3991,3990,4176,4177,5612,5613,4178,4179,5797,5796,5799,5817,6432,6431,6430,6204,6205,6206,6207,6200,6201,6202,6203,6208,6209,4174,4175,5229,5228,4170,4171,4172,4173,5223,5222,5221,5220,5227,5226,5225,5224,5818,5819,5795,5794,5793,5810,5811,5812,5813,5814,5815,5816,5798,4505,4504,4507,4506,4501,4500,4503,4502,4509,4508,4668,4669,4667,5586,5587,5588,5589,14101,14100,14103,14102,14105,14104,14107,14106,14109,14108,4484,4485,4486,4487,4481,4482,4483,4488,4489,5038,5039,5036,5037,5034,5035,5032,5033,5030,5031,14022,14023,14020,14021,14024,6402,6403,6401,6406,6407,6404,6405,6408,6409,6214,6213,6212,6211,6210,5216,5217,14095,5215,14093,14092,5218,5219,4167,4166,4165,4164,4163,4162,4161,4160,4000,4001,4002,4003,4004,4169,4168,5809,5808,3975,3974,3977,3976,5803,5802,5801,5800,5807,5806,5805,5804,4679,4678,4677,4676,4675,4674,4673,4672,4671,4670,5591,5590,5593,5592,5595,5594,5597,5596,5599,5598,5610,5611,5616,5617,5614,5615,14112,14113,14110,14111,14116,14117,14114,14115,14118,14119,13997,13996,13995,13994,13999,13998,4497,4496,4495,4494,4493,4492,4491,4490,4499,4498,5245,5244,5241,5240,5243,5242,5029,5028,5021,5020,5023,5022,5025,5024,5027,5026,4686,4687,4684,4685,4682,4683,4680,4681,4688,4689,6194,6195,6196,6197,6190,6191,6192,6193,6198,6199,14019,14018,14017,14016,14015,14014,14013,14012,14011,14010,6419,6418,6415,6414,6417,6416,6411,6410,6413,6412,14099,14098,4152,4153,4150,4151,4156,4157,4154,4155,4158,4159,3988,3989,3984,3985,3986,3987,3980,3981,3982,3983,14097,14096,3979,3978,14094,5019,5605,5604,5607,5606,5601,5600,5603,5602,5609,5608,14122,14121,14120,6424,6425,6426,6427,6420,6421,6422,6423,6428,6429,5238,5239,5230,5231,5232,5233,5234,5235,5236,5237,4691,4690,4693,4692,4695,4694,4697,4696,4510,4511,6187,6186,6185,6184,6189,6188,14008,14009,14000,14001,14002,14003,14004,14005,14006,14007,5049,5048,5047,5046,5045,5044,5043,5042,5041,5040',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 04:12:25.301 +85.30.0.0000',
       scan=['128,131,134,243,105,108,111,240,115,118,121,124'],
       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_T32t01.J041225+853000_P17025v1_2018_01_18T23_06_05.509/working/VLASS1.1.sb34666958.eb34700016.58074.035491539355_split.ms'],
       imagename='VLASS1.1.ql.T32t01.J041225+853000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='3974,3975,3976,3977,3978,3979,3980,3981,3982,3983,3984,3985,3986,3987,3988,3989,3990,3991,3992,3993,3994,3995,3996,3997,3998,3999,4000,4001,4002,4003,4004,4149,4150,4151,4152,4153,4154,4155,4156,4157,4158,4159,4160,4161,4162,4163,4164,4165,4166,4167,4168,4169,4170,4171,4172,4173,4174,4175,4176,4177,4178,4179,4481,4482,4483,4484,4485,4486,4487,4488,4489,4490,4491,4492,4493,4494,4495,4496,4497,4498,4499,4500,4501,4502,4503,4504,4505,4506,4507,4508,4509,4510,4511,4667,4668,4669,4670,4671,4672,4673,4674,4675,4676,4677,4678,4679,4680,4681,4682,4683,4684,4685,4686,4687,4688,4689,4690,4691,4692,4693,4694,4695,4696,4697,5019,5020,5021,5022,5023,5024,5025,5026,5027,5028,5029,5030,5031,5032,5033,5034,5035,5036,5037,5038,5039,5040,5041,5042,5043,5044,5045,5046,5047,5048,5049,5215,5216,5217,5218,5219,5220,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,5586,5587,5588,5589,5590,5591,5592,5593,5594,5595,5596,5597,5598,5599,5600,5601,5602,5603,5604,5605,5606,5607,5608,5609,5610,5611,5612,5613,5614,5615,5616,5617,5793,5794,5795,5796,5797,5798,5799,5800,5801,5802,5803,5804,5805,5806,5807,5808,5809,5810,5811,5812,5813,5814,5815,5816,5817,5818,5819,5820,5821,5822,5823,6184,6185,6186,6187,6188,6189,6190,6191,6192,6193,6194,6195,6196,6197,6198,6199,6200,6201,6202,6203,6204,6205,6206,6207,6208,6209,6210,6211,6212,6213,6214,6401,6402,6403,6404,6405,6406,6407,6408,6409,6410,6411,6412,6413,6414,6415,6416,6417,6418,6419,6420,6421,6422,6423,6424,6425,6426,6427,6428,6429,6430,6431,6432,13994,13995,13996,13997,13998,13999,14000,14001,14002,14003,14004,14005,14006,14007,14008,14009,14010,14011,14012,14013,14014,14015,14016,14017,14018,14019,14020,14021,14022,14023,14024,14092,14093,14094,14095,14096,14097,14098,14099,14100,14101,14102,14103,14104,14105,14106,14107,14108,14109,14110,14111,14112,14113,14114,14115,14116,14117,14118,14119,14120,14121,14122',
       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 04:12:25.301 +85.30.0.0000',
       scan=['128,131,134,243,105,108,111,240,115,118,121,124'],
       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_T32t01.J041225+853000_P17025v1_2018_01_18T23_06_05.509/working/VLASS1.1.sb34666958.eb34700016.58074.035491539355_split.ms'],
       imagename='VLASS1.1.ql.T32t01.J041225+853000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00048479786160460204, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='3974,3975,3976,3977,3978,3979,3980,3981,3982,3983,3984,3985,3986,3987,3988,3989,3990,3991,3992,3993,3994,3995,3996,3997,3998,3999,4000,4001,4002,4003,4004,4149,4150,4151,4152,4153,4154,4155,4156,4157,4158,4159,4160,4161,4162,4163,4164,4165,4166,4167,4168,4169,4170,4171,4172,4173,4174,4175,4176,4177,4178,4179,4481,4482,4483,4484,4485,4486,4487,4488,4489,4490,4491,4492,4493,4494,4495,4496,4497,4498,4499,4500,4501,4502,4503,4504,4505,4506,4507,4508,4509,4510,4511,4667,4668,4669,4670,4671,4672,4673,4674,4675,4676,4677,4678,4679,4680,4681,4682,4683,4684,4685,4686,4687,4688,4689,4690,4691,4692,4693,4694,4695,4696,4697,5019,5020,5021,5022,5023,5024,5025,5026,5027,5028,5029,5030,5031,5032,5033,5034,5035,5036,5037,5038,5039,5040,5041,5042,5043,5044,5045,5046,5047,5048,5049,5215,5216,5217,5218,5219,5220,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,5586,5587,5588,5589,5590,5591,5592,5593,5594,5595,5596,5597,5598,5599,5600,5601,5602,5603,5604,5605,5606,5607,5608,5609,5610,5611,5612,5613,5614,5615,5616,5617,5793,5794,5795,5796,5797,5798,5799,5800,5801,5802,5803,5804,5805,5806,5807,5808,5809,5810,5811,5812,5813,5814,5815,5816,5817,5818,5819,5820,5821,5822,5823,6184,6185,6186,6187,6188,6189,6190,6191,6192,6193,6194,6195,6196,6197,6198,6199,6200,6201,6202,6203,6204,6205,6206,6207,6208,6209,6210,6211,6212,6213,6214,6401,6402,6403,6404,6405,6406,6407,6408,6409,6410,6411,6412,6413,6414,6415,6416,6417,6418,6419,6420,6421,6422,6423,6424,6425,6426,6427,6428,6429,6430,6431,6432,13994,13995,13996,13997,13998,13999,14000,14001,14002,14003,14004,14005,14006,14007,14008,14009,14010,14011,14012,14013,14014,14015,14016,14017,14018,14019,14020,14021,14022,14023,14024,14092,14093,14094,14095,14096,14097,14098,14099,14100,14101,14102,14103,14104,14105,14106,14107,14108,14109,14110,14111,14112,14113,14114,14115,14116,14117,14118,14119,14120,14121,14122',
       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.T32t01.J041225+853000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T32t01.J041225+853000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T32t01.J041225+853000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#