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

# hif_editimlist(parameter_file='parameter.14023.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_T24t13.J120319+533000_P14023v1_2018_02_06T18_19_36.291/working/VLASS1.1.sb34541777.eb34701358.58077.63418230324_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T24t13.J120319+533000_P14023v1_2018_02_06T18_19_36.291/rawdata/VLASS1.1.sb34541777.eb34701358.58077.63418230324.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='5453,5418,5451,6157,6161,5419,5078,5079,6159,5072,5073,5070,5071,5076,5077,5074,5075,4707,4706,4705,4704,4703,4702,4679,4678,4677,4676,4675,4674,4673,5809,4709,4708,5808,3999,3998,3997,3996,3995,3994,3993,3992,3991,3990,5065,5064,5067,5066,5061,5060,5821,4336,4337,4334,4335,4332,4333,4330,4331,6150,6151,6152,6153,6154,6155,6156,5802,6158,5429,5819,6162,5804,6145,4329,4328,5792,4321,4320,4323,4322,4325,4324,4327,4326,6143,6142,6141,6140,6147,5423,5768,5769,6139,6149,5422,6138,5414,5415,5416,5417,5410,5411,5412,5413,4725,4724,4726,4721,4720,4723,4722,5818,5793,5789,5797,5409,3959,3958,5791,5820,3953,5790,3957,3956,3955,3954,5779,5778,5408,5796,4358,4359,6160,6146,4354,4355,4356,4357,4350,4351,4352,4353,5407,5406,5089,5088,5403,4701,5401,5083,5082,5081,5080,5087,5086,5085,5084,5449,5799,5798,5452,5038,5039,5036,5037,6137,5771,5816,4349,4348,4347,4346,4345,4344,4343,4342,4341,4340,5432,5433,5430,5431,5436,5437,5434,5435,5775,5438,5439,5774,5784,5785,5786,5787,5780,5777,5782,5783,5770,5454,5788,5776,5781,6148,5405,5069,5404,5773,5402,4686,4687,4684,4685,4682,4683,4680,4681,5425,5424,5426,5421,5420,4688,4689,3966,3967,3964,3965,3962,3963,3960,3961,5810,5811,5812,5813,5814,5815,3968,3969,4365,4364,4361,4360,4363,4362,5447,6144,5446,5058,5059,5054,5055,5056,5057,5050,5051,5052,5053,4691,4690,4693,4692,4695,4694,4697,4696,4000,4001,4002,4003,4004,4005,3971,3970,3973,3972,3975,3974,3977,3976,5803,3978,5801,5800,5807,5806,5805,5448,4312,4313,4314,4315,4316,4317,4318,4319,5772,5068,5817,5049,5048,5047,5046,5045,5044,5043,5042,5041,5040,4714,4715,4716,4717,4710,4711,4712,4713,5443,5442,5441,5440,4718,4719,5445,5444,5450,3988,3989,4698,3984,3985,3986,3987,3981,3982,3983',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 12:03:19.397 +53.30.0.0000',
       scan=['97,100,103,106,109,78,81,84,113,87,116,90,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_T24t13.J120319+533000_P14023v1_2018_02_06T18_19_36.291/working/VLASS1.1.sb34541777.eb34701358.58077.63418230324_split.ms'],
       imagename='VLASS1.1.ql.T24t13.J120319+533000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='3953,3954,3955,3956,3957,3958,3959,3960,3961,3962,3963,3964,3965,3966,3967,3968,3969,3970,3971,3972,3973,3974,3975,3976,3977,3978,3981,3982,3983,3984,3985,3986,3987,3988,3989,3990,3991,3992,3993,3994,3995,3996,3997,3998,3999,4000,4001,4002,4003,4004,4005,4312,4313,4314,4315,4316,4317,4318,4319,4320,4321,4322,4323,4324,4325,4326,4327,4328,4329,4330,4331,4332,4333,4334,4335,4336,4337,4340,4341,4342,4343,4344,4345,4346,4347,4348,4349,4350,4351,4352,4353,4354,4355,4356,4357,4358,4359,4360,4361,4362,4363,4364,4365,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,4698,4701,4702,4703,4704,4705,4706,4707,4708,4709,4710,4711,4712,4713,4714,4715,4716,4717,4718,4719,4720,4721,4722,4723,4724,4725,4726,5036,5037,5038,5039,5040,5041,5042,5043,5044,5045,5046,5047,5048,5049,5050,5051,5052,5053,5054,5055,5056,5057,5058,5059,5060,5061,5064,5065,5066,5067,5068,5069,5070,5071,5072,5073,5074,5075,5076,5077,5078,5079,5080,5081,5082,5083,5084,5085,5086,5087,5088,5089,5401,5402,5403,5404,5405,5406,5407,5408,5409,5410,5411,5412,5413,5414,5415,5416,5417,5418,5419,5420,5421,5422,5423,5424,5425,5426,5429,5430,5431,5432,5433,5434,5435,5436,5437,5438,5439,5440,5441,5442,5443,5444,5445,5446,5447,5448,5449,5450,5451,5452,5453,5454,5768,5769,5770,5771,5772,5773,5774,5775,5776,5777,5778,5779,5780,5781,5782,5783,5784,5785,5786,5787,5788,5789,5790,5791,5792,5793,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,6137,6138,6139,6140,6141,6142,6143,6144,6145,6146,6147,6148,6149,6150,6151,6152,6153,6154,6155,6156,6157,6158,6159,6160,6161,6162',
       outframe='LSRK', gridder='mosaic', cycleniter=100, 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 12:03:19.397 +53.30.0.0000',
       scan=['97,100,103,106,109,78,81,84,113,87,116,90,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_T24t13.J120319+533000_P14023v1_2018_02_06T18_19_36.291/working/VLASS1.1.sb34541777.eb34701358.58077.63418230324_split.ms'],
       imagename='VLASS1.1.ql.T24t13.J120319+533000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00057525230375722451, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='3953,3954,3955,3956,3957,3958,3959,3960,3961,3962,3963,3964,3965,3966,3967,3968,3969,3970,3971,3972,3973,3974,3975,3976,3977,3978,3981,3982,3983,3984,3985,3986,3987,3988,3989,3990,3991,3992,3993,3994,3995,3996,3997,3998,3999,4000,4001,4002,4003,4004,4005,4312,4313,4314,4315,4316,4317,4318,4319,4320,4321,4322,4323,4324,4325,4326,4327,4328,4329,4330,4331,4332,4333,4334,4335,4336,4337,4340,4341,4342,4343,4344,4345,4346,4347,4348,4349,4350,4351,4352,4353,4354,4355,4356,4357,4358,4359,4360,4361,4362,4363,4364,4365,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,4698,4701,4702,4703,4704,4705,4706,4707,4708,4709,4710,4711,4712,4713,4714,4715,4716,4717,4718,4719,4720,4721,4722,4723,4724,4725,4726,5036,5037,5038,5039,5040,5041,5042,5043,5044,5045,5046,5047,5048,5049,5050,5051,5052,5053,5054,5055,5056,5057,5058,5059,5060,5061,5064,5065,5066,5067,5068,5069,5070,5071,5072,5073,5074,5075,5076,5077,5078,5079,5080,5081,5082,5083,5084,5085,5086,5087,5088,5089,5401,5402,5403,5404,5405,5406,5407,5408,5409,5410,5411,5412,5413,5414,5415,5416,5417,5418,5419,5420,5421,5422,5423,5424,5425,5426,5429,5430,5431,5432,5433,5434,5435,5436,5437,5438,5439,5440,5441,5442,5443,5444,5445,5446,5447,5448,5449,5450,5451,5452,5453,5454,5768,5769,5770,5771,5772,5773,5774,5775,5776,5777,5778,5779,5780,5781,5782,5783,5784,5785,5786,5787,5788,5789,5790,5791,5792,5793,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,6137,6138,6139,6140,6141,6142,6143,6144,6145,6146,6147,6148,6149,6150,6151,6152,6153,6154,6155,6156,6157,6158,6159,6160,6161,6162',
       outframe='LSRK', gridder='mosaic', cycleniter=100, 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.T24t13.J120319+533000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T24t13.J120319+533000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T24t13.J120319+533000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#