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

# hif_editimlist(parameter_file='parameter.16364.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_T27t05.J065825+663000_P16364v1_2018_02_08T17_04_28.290/working/VLASS1.1.sb34640437.eb34680548.58066.26112458333_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T27t05.J065825+663000_P16364v1_2018_02_08T17_04_28.290/rawdata/VLASS1.1.sb34640437.eb34680548.58066.26112458333.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='2828,2829,2824,2825,2826,2827,2820,2821,2822,2823,3580,4214,3582,3583,3584,3585,4213,4212,3588,3589,5009,5008,4219,4218,4358,4359,4354,4355,4356,4357,4351,4352,4353,3452,3453,3450,3451,3456,3457,3454,3455,3458,2688,2689,2684,2685,2686,2687,2680,2681,2682,2683,4618,4777,4774,4775,4772,4773,4770,4771,4778,4779,3207,3206,3205,3204,3203,3202,3201,3200,3209,3208,3966,3967,3964,3965,3962,3963,3960,3961,3968,3969,5007,4619,3062,3061,3060,3067,3066,3065,3064,4611,4610,3069,4612,4615,4614,4617,4616,5006,3068,5005,5004,4228,4229,4224,4225,4226,4227,4220,4221,4222,4223,4381,4380,3834,3835,3836,3837,3830,3831,3832,3833,3838,3839,4613,3567,3568,3569,3841,3840,3843,3842,3845,3844,3846,4620,4621,4622,4623,4624,4625,4626,4627,4628,4629,2830,3593,3592,3591,3590,3597,3596,3595,3594,5034,3598,5032,5033,5030,5031,4215,3449,3448,3445,3444,3447,3446,3441,3440,3443,3442,3063,2694,2693,2692,2691,2690,4761,4760,4763,4762,4765,4764,4767,4766,4769,4768,3971,3970,3973,3972,3975,3974,3977,3976,3979,3978,3058,3059,3056,3057,3054,3055,3052,3053,3050,3051,4239,4238,4237,4236,4235,4234,4233,4232,4231,4230,3827,3826,3825,3824,3823,3822,3821,3820,4217,3829,3828,4216,4211,4210,3579,3578,3575,3574,3577,3576,3571,3570,3573,3572,2666,2667,2664,2665,2663,2668,2669,4758,4759,4750,4751,4752,4753,4754,4755,4756,4757,4633,4632,4631,4630,4635,4634,2802,2803,2800,2801,2806,2807,2804,2805,3581,5029,5028,5021,5020,5023,5022,5025,5024,5027,5026,4378,4379,4372,4373,4370,4371,4376,4377,4374,4375,3438,3439,3430,3431,3432,3433,3434,3435,3436,3437,3184,3185,3186,3187,3181,3182,3183,3188,3189,3586,3587,3984,3985,3986,3987,3980,3981,3982,3983,3049,3048,3045,3044,3047,3046,2808,2809,5010,5011,5012,5013,5014,5015,5016,5017,5018,5019,4208,4209,3959,3958,3816,3817,3818,3819,3957,2671,2670,2673,2672,2675,2674,2677,2676,2679,2678,4749,4748,3210,3211,3212,2815,2814,2817,2816,2811,2810,2813,2812,2819,2818,3070,3071,3072,3073,3074,4608,4609,4606,4607,4605,4365,4364,4367,4366,4361,4360,4363,4362,4369,4368,4776,3429,3428,5035,3197,3196,3195,3194,3193,3192,3191,3190,3199,3198',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 06:58:25.218 +66.30.0.0000',
       scan=['67,86,70,73,76,79,92,82,54,89,57,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_T27t05.J065825+663000_P16364v1_2018_02_08T17_04_28.290/working/VLASS1.1.sb34640437.eb34680548.58066.26112458333_split.ms'],
       imagename='VLASS1.1.ql.T27t05.J065825+663000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='2663,2664,2665,2666,2667,2668,2669,2670,2671,2672,2673,2674,2675,2676,2677,2678,2679,2680,2681,2682,2683,2684,2685,2686,2687,2688,2689,2690,2691,2692,2693,2694,2800,2801,2802,2803,2804,2805,2806,2807,2808,2809,2810,2811,2812,2813,2814,2815,2816,2817,2818,2819,2820,2821,2822,2823,2824,2825,2826,2827,2828,2829,2830,3044,3045,3046,3047,3048,3049,3050,3051,3052,3053,3054,3055,3056,3057,3058,3059,3060,3061,3062,3063,3064,3065,3066,3067,3068,3069,3070,3071,3072,3073,3074,3181,3182,3183,3184,3185,3186,3187,3188,3189,3190,3191,3192,3193,3194,3195,3196,3197,3198,3199,3200,3201,3202,3203,3204,3205,3206,3207,3208,3209,3210,3211,3212,3428,3429,3430,3431,3432,3433,3434,3435,3436,3437,3438,3439,3440,3441,3442,3443,3444,3445,3446,3447,3448,3449,3450,3451,3452,3453,3454,3455,3456,3457,3458,3567,3568,3569,3570,3571,3572,3573,3574,3575,3576,3577,3578,3579,3580,3581,3582,3583,3584,3585,3586,3587,3588,3589,3590,3591,3592,3593,3594,3595,3596,3597,3598,3816,3817,3818,3819,3820,3821,3822,3823,3824,3825,3826,3827,3828,3829,3830,3831,3832,3833,3834,3835,3836,3837,3838,3839,3840,3841,3842,3843,3844,3845,3846,3957,3958,3959,3960,3961,3962,3963,3964,3965,3966,3967,3968,3969,3970,3971,3972,3973,3974,3975,3976,3977,3978,3979,3980,3981,3982,3983,3984,3985,3986,3987,4208,4209,4210,4211,4212,4213,4214,4215,4216,4217,4218,4219,4220,4221,4222,4223,4224,4225,4226,4227,4228,4229,4230,4231,4232,4233,4234,4235,4236,4237,4238,4239,4351,4352,4353,4354,4355,4356,4357,4358,4359,4360,4361,4362,4363,4364,4365,4366,4367,4368,4369,4370,4371,4372,4373,4374,4375,4376,4377,4378,4379,4380,4381,4605,4606,4607,4608,4609,4610,4611,4612,4613,4614,4615,4616,4617,4618,4619,4620,4621,4622,4623,4624,4625,4626,4627,4628,4629,4630,4631,4632,4633,4634,4635,4748,4749,4750,4751,4752,4753,4754,4755,4756,4757,4758,4759,4760,4761,4762,4763,4764,4765,4766,4767,4768,4769,4770,4771,4772,4773,4774,4775,4776,4777,4778,4779,5004,5005,5006,5007,5008,5009,5010,5011,5012,5013,5014,5015,5016,5017,5018,5019,5020,5021,5022,5023,5024,5025,5026,5027,5028,5029,5030,5031,5032,5033,5034,5035',
       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 06:58:25.218 +66.30.0.0000',
       scan=['67,86,70,73,76,79,92,82,54,89,57,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_T27t05.J065825+663000_P16364v1_2018_02_08T17_04_28.290/working/VLASS1.1.sb34640437.eb34680548.58066.26112458333_split.ms'],
       imagename='VLASS1.1.ql.T27t05.J065825+663000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00051268261225450261, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='2663,2664,2665,2666,2667,2668,2669,2670,2671,2672,2673,2674,2675,2676,2677,2678,2679,2680,2681,2682,2683,2684,2685,2686,2687,2688,2689,2690,2691,2692,2693,2694,2800,2801,2802,2803,2804,2805,2806,2807,2808,2809,2810,2811,2812,2813,2814,2815,2816,2817,2818,2819,2820,2821,2822,2823,2824,2825,2826,2827,2828,2829,2830,3044,3045,3046,3047,3048,3049,3050,3051,3052,3053,3054,3055,3056,3057,3058,3059,3060,3061,3062,3063,3064,3065,3066,3067,3068,3069,3070,3071,3072,3073,3074,3181,3182,3183,3184,3185,3186,3187,3188,3189,3190,3191,3192,3193,3194,3195,3196,3197,3198,3199,3200,3201,3202,3203,3204,3205,3206,3207,3208,3209,3210,3211,3212,3428,3429,3430,3431,3432,3433,3434,3435,3436,3437,3438,3439,3440,3441,3442,3443,3444,3445,3446,3447,3448,3449,3450,3451,3452,3453,3454,3455,3456,3457,3458,3567,3568,3569,3570,3571,3572,3573,3574,3575,3576,3577,3578,3579,3580,3581,3582,3583,3584,3585,3586,3587,3588,3589,3590,3591,3592,3593,3594,3595,3596,3597,3598,3816,3817,3818,3819,3820,3821,3822,3823,3824,3825,3826,3827,3828,3829,3830,3831,3832,3833,3834,3835,3836,3837,3838,3839,3840,3841,3842,3843,3844,3845,3846,3957,3958,3959,3960,3961,3962,3963,3964,3965,3966,3967,3968,3969,3970,3971,3972,3973,3974,3975,3976,3977,3978,3979,3980,3981,3982,3983,3984,3985,3986,3987,4208,4209,4210,4211,4212,4213,4214,4215,4216,4217,4218,4219,4220,4221,4222,4223,4224,4225,4226,4227,4228,4229,4230,4231,4232,4233,4234,4235,4236,4237,4238,4239,4351,4352,4353,4354,4355,4356,4357,4358,4359,4360,4361,4362,4363,4364,4365,4366,4367,4368,4369,4370,4371,4372,4373,4374,4375,4376,4377,4378,4379,4380,4381,4605,4606,4607,4608,4609,4610,4611,4612,4613,4614,4615,4616,4617,4618,4619,4620,4621,4622,4623,4624,4625,4626,4627,4628,4629,4630,4631,4632,4633,4634,4635,4748,4749,4750,4751,4752,4753,4754,4755,4756,4757,4758,4759,4760,4761,4762,4763,4764,4765,4766,4767,4768,4769,4770,4771,4772,4773,4774,4775,4776,4777,4778,4779,5004,5005,5006,5007,5008,5009,5010,5011,5012,5013,5014,5015,5016,5017,5018,5019,5020,5021,5022,5023,5024,5025,5026,5027,5028,5029,5030,5031,5032,5033,5034,5035',
       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.T27t05.J065825+663000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T27t05.J065825+663000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T27t05.J065825+663000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#