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

# hif_editimlist(parameter_file='parameter.5598.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_T28t11.J215517+713000_P5598v1_2017_09_26T20_03_35.777/working/VLASS1.1.sb34402674.eb34433625.58014.15711049769_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T28t11.J215517+713000_P5598v1_2017_09_26T20_03_35.777/rawdata/VLASS1.1.sb34402674.eb34433625.58014.15711049769.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='10518,10519,10512,10513,10516,10517,10514,10515,3557,3556,3555,3554,3553,3552,3551,3550,12068,12069,5465,5467,5466,5469,5468,12062,12063,12066,12067,4732,4733,4736,4737,4734,4735,4738,4739,3572,4329,4328,4321,4323,4322,4325,4324,4327,4326,11689,11688,11684,11687,11681,11680,11683,11682,3580,3581,3582,4358,4359,4354,4355,4356,4357,4350,4351,4352,4353,5135,5134,5133,5132,5131,5130,3966,3967,3964,3965,3962,3963,3960,3961,3968,3969,12056,12057,3538,12059,12058,11308,11309,12061,11300,11301,11304,11305,11306,11307,10527,10526,10525,10524,10523,3938,3935,3934,3937,3936,3931,3930,3933,3932,10118,10119,10116,10117,10114,10115,3562,3563,3560,3561,3566,3567,3564,3565,3568,3569,5498,5499,5494,5495,5496,5497,5490,5491,5492,5493,4729,4728,4725,4724,4727,4726,4721,4720,4723,4722,3159,11694,11692,11693,11690,11691,10905,10907,10906,10909,10908,4349,4348,4347,4346,4343,4342,4341,4340,5108,5109,5102,5103,5100,5101,5106,5107,5104,5105,3971,3970,3973,3972,3975,3974,3977,3976,4714,4715,4716,4717,4710,4711,4712,4713,4718,4719,5090,5091,5092,5093,3571,3570,3573,5508,5509,5506,5507,5504,5505,5502,5503,5500,5501,12437,12436,12435,12434,12433,12432,12431,12430,12060,3948,3949,3940,3941,3942,3943,3944,3945,3946,3947,10129,10128,10121,10120,10127,10126,10125,10124,5098,5099,3579,3578,3575,3574,3577,3576,5094,5095,5096,5097,12071,12070,5487,5486,5485,5484,5483,5482,5481,5480,4750,4751,4752,4753,4754,11298,11299,11296,11297,11294,11295,3168,3169,3166,3167,3164,3165,3162,3163,3160,3161,10916,10917,10915,10912,10910,10911,10918,10919,5115,5114,5117,5116,5110,5113,5119,5118,3535,3537,3536,3539,3939,10522,3180,3181,3182,10529,10528,3544,3545,3546,3547,3540,3541,3542,3543,3548,3549,5476,5477,5474,5475,5472,5473,5470,5471,12073,12072,5478,5479,4707,4709,4708,4338,4339,4336,4337,4334,4335,4332,4333,4330,4331,5511,5510,5512,11677,11678,11679,3959,3958,10130,10131,3952,3951,3950,3957,3956,3955,5128,5129,5124,5125,5126,5127,5120,5121,5122,5123,5089,5088,4749,4748,4743,4742,4741,4740,4747,4746,4745,4744,3179,3178,3171,3170,3173,3172,3175,3174,3177,3176,10922,10921,10920,4365,4364,4367,4366,4361,4360,4363,4362,4368,11311,11310',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 21:55:17.278 +71.30.0.0000',
       scan=['194,197,165,169,172,175,200,178,181,184,188,191,66,203,69,72,75,79,56,82,85,88,91,94,60,63'],
       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_T28t11.J215517+713000_P5598v1_2017_09_26T20_03_35.777/working/VLASS1.1.sb34402674.eb34433625.58014.15711049769_split.ms'],
       imagename='VLASS1.1.ql.T28t11.J215517+713000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='3159,3160,3161,3162,3163,3164,3165,3166,3167,3168,3169,3170,3171,3172,3173,3174,3175,3176,3177,3178,3179,3180,3181,3182,3535,3536,3537,3538,3539,3540,3541,3542,3543,3544,3545,3546,3547,3548,3549,3550,3551,3552,3553,3554,3555,3556,3557,3560,3561,3562,3563,3564,3565,3566,3567,3568,3569,3570,3571,3572,3573,3574,3575,3576,3577,3578,3579,3580,3581,3582,3930,3931,3932,3933,3934,3935,3936,3937,3938,3939,3940,3941,3942,3943,3944,3945,3946,3947,3948,3949,3950,3951,3952,3955,3956,3957,3958,3959,3960,3961,3962,3963,3964,3965,3966,3967,3968,3969,3970,3971,3972,3973,3974,3975,3976,3977,4321,4322,4323,4324,4325,4326,4327,4328,4329,4330,4331,4332,4333,4334,4335,4336,4337,4338,4339,4340,4341,4342,4343,4346,4347,4348,4349,4350,4351,4352,4353,4354,4355,4356,4357,4358,4359,4360,4361,4362,4363,4364,4365,4366,4367,4368,4707,4708,4709,4710,4711,4712,4713,4714,4715,4716,4717,4718,4719,4720,4721,4722,4723,4724,4725,4726,4727,4728,4729,4732,4733,4734,4735,4736,4737,4738,4739,4740,4741,4742,4743,4744,4745,4746,4747,4748,4749,4750,4751,4752,4753,4754,5088,5089,5090,5091,5092,5093,5094,5095,5096,5097,5098,5099,5100,5101,5102,5103,5104,5105,5106,5107,5108,5109,5110,5113,5114,5115,5116,5117,5118,5119,5120,5121,5122,5123,5124,5125,5126,5127,5128,5129,5130,5131,5132,5133,5134,5135,5465,5466,5467,5468,5469,5470,5471,5472,5473,5474,5475,5476,5477,5478,5479,5480,5481,5482,5483,5484,5485,5486,5487,5490,5491,5492,5493,5494,5495,5496,5497,5498,5499,5500,5501,5502,5503,5504,5505,5506,5507,5508,5509,5510,5511,5512,10114,10115,10116,10117,10118,10119,10120,10121,10124,10125,10126,10127,10128,10129,10130,10131,10512,10513,10514,10515,10516,10517,10518,10519,10522,10523,10524,10525,10526,10527,10528,10529,10905,10906,10907,10908,10909,10910,10911,10912,10915,10916,10917,10918,10919,10920,10921,10922,11294,11295,11296,11297,11298,11299,11300,11301,11304,11305,11306,11307,11308,11309,11310,11311,11677,11678,11679,11680,11681,11682,11683,11684,11687,11688,11689,11690,11691,11692,11693,11694,12056,12057,12058,12059,12060,12061,12062,12063,12066,12067,12068,12069,12070,12071,12072,12073,12430,12431,12432,12433,12434,12435,12436,12437',
       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 21:55:17.278 +71.30.0.0000',
       scan=['194,197,165,169,172,175,200,178,181,184,188,191,66,203,69,72,75,79,56,82,85,88,91,94,60,63'],
       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_T28t11.J215517+713000_P5598v1_2017_09_26T20_03_35.777/working/VLASS1.1.sb34402674.eb34433625.58014.15711049769_split.ms'],
       imagename='VLASS1.1.ql.T28t11.J215517+713000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00056888473961480481, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='3159,3160,3161,3162,3163,3164,3165,3166,3167,3168,3169,3170,3171,3172,3173,3174,3175,3176,3177,3178,3179,3180,3181,3182,3535,3536,3537,3538,3539,3540,3541,3542,3543,3544,3545,3546,3547,3548,3549,3550,3551,3552,3553,3554,3555,3556,3557,3560,3561,3562,3563,3564,3565,3566,3567,3568,3569,3570,3571,3572,3573,3574,3575,3576,3577,3578,3579,3580,3581,3582,3930,3931,3932,3933,3934,3935,3936,3937,3938,3939,3940,3941,3942,3943,3944,3945,3946,3947,3948,3949,3950,3951,3952,3955,3956,3957,3958,3959,3960,3961,3962,3963,3964,3965,3966,3967,3968,3969,3970,3971,3972,3973,3974,3975,3976,3977,4321,4322,4323,4324,4325,4326,4327,4328,4329,4330,4331,4332,4333,4334,4335,4336,4337,4338,4339,4340,4341,4342,4343,4346,4347,4348,4349,4350,4351,4352,4353,4354,4355,4356,4357,4358,4359,4360,4361,4362,4363,4364,4365,4366,4367,4368,4707,4708,4709,4710,4711,4712,4713,4714,4715,4716,4717,4718,4719,4720,4721,4722,4723,4724,4725,4726,4727,4728,4729,4732,4733,4734,4735,4736,4737,4738,4739,4740,4741,4742,4743,4744,4745,4746,4747,4748,4749,4750,4751,4752,4753,4754,5088,5089,5090,5091,5092,5093,5094,5095,5096,5097,5098,5099,5100,5101,5102,5103,5104,5105,5106,5107,5108,5109,5110,5113,5114,5115,5116,5117,5118,5119,5120,5121,5122,5123,5124,5125,5126,5127,5128,5129,5130,5131,5132,5133,5134,5135,5465,5466,5467,5468,5469,5470,5471,5472,5473,5474,5475,5476,5477,5478,5479,5480,5481,5482,5483,5484,5485,5486,5487,5490,5491,5492,5493,5494,5495,5496,5497,5498,5499,5500,5501,5502,5503,5504,5505,5506,5507,5508,5509,5510,5511,5512,10114,10115,10116,10117,10118,10119,10120,10121,10124,10125,10126,10127,10128,10129,10130,10131,10512,10513,10514,10515,10516,10517,10518,10519,10522,10523,10524,10525,10526,10527,10528,10529,10905,10906,10907,10908,10909,10910,10911,10912,10915,10916,10917,10918,10919,10920,10921,10922,11294,11295,11296,11297,11298,11299,11300,11301,11304,11305,11306,11307,11308,11309,11310,11311,11677,11678,11679,11680,11681,11682,11683,11684,11687,11688,11689,11690,11691,11692,11693,11694,12056,12057,12058,12059,12060,12061,12062,12063,12066,12067,12068,12069,12070,12071,12072,12073,12430,12431,12432,12433,12434,12435,12436,12437',
       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.T28t11.J215517+713000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T28t11.J215517+713000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T28t11.J215517+713000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#