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

# hif_editimlist(parameter_file='parameter.7471.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_T26t01.J000415+623000_P7471v1_2018_01_01T18_46_59.683/working/VLASS1.1.sb34459038.eb34579917.58037.11761364584_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T26t01.J000415+623000_P7471v1_2018_01_01T18_46_59.683/rawdata/VLASS1.1.sb34459038.eb34579917.58037.11761364584.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='2659,2658,4168,4565,2653,2652,2651,2650,2657,2656,2655,2654,3049,3048,3414,3415,3412,3413,3410,3411,3041,3040,3043,3042,4569,3047,3046,4909,4908,4568,3799,3798,3795,3794,3797,3796,3790,3793,4919,4181,4562,4145,4151,4183,4930,4182,3765,3766,3767,3768,3769,2701,2700,4184,3401,3400,3403,3402,3405,3404,3407,3406,3409,3408,3418,4186,4171,4570,3038,3039,4571,3034,3035,3036,3037,3030,3031,3032,3033,4932,4573,4563,4574,4180,4928,4575,3773,3772,3771,3770,3777,3776,3775,3774,4178,3779,3778,4564,4179,4929,4558,4559,4556,4557,4554,4555,4552,4553,4922,3029,3028,3027,3026,3025,3024,3023,3022,3021,3020,2666,2667,2664,2665,2662,2663,2660,2661,3419,4194,2668,2669,3812,3813,3810,3811,3816,3817,3814,3815,4196,4918,3818,4195,4192,4193,4190,4191,4549,4548,4541,4149,4543,4542,4545,4544,4547,4148,3391,4158,3393,3019,3395,3394,3397,3396,3399,3398,4143,2671,2670,2673,2672,2674,2677,4144,2679,2678,4927,4147,4159,4189,4188,4146,3809,3808,4914,3805,3804,3807,3806,3801,3800,3803,3802,4561,3781,4540,4926,4915,4533,4917,4538,4529,4534,4535,4536,4537,4530,4531,4532,3788,2688,2689,4925,2684,2685,2686,2687,2680,2681,2682,2683,4152,4528,4923,4153,4933,3780,3070,3392,4924,4577,4156,4157,4174,4175,4176,4177,4154,4172,4173,4527,4526,4525,4155,3441,3440,3443,3442,4916,4911,4539,2699,2698,2697,2696,2695,2694,2693,2692,2691,2690,4572,3789,4912,4913,3063,3062,3061,3060,3067,3066,3065,3064,3069,3068,3438,3439,4165,4164,4163,4162,4161,4160,3430,3431,3432,3433,3434,3435,3436,3437,4910,4524,4546,4167,4931,4185,4166,4921,4920,3058,3059,3056,3057,3054,3055,3052,3053,3050,3051,3782,3783,3429,3428,3786,3787,3784,3785,3423,3422,3421,3420,3427,3426,3425,3424,4150,4560,4567,4576,4187,4566',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 00:04:15.607 +62.30.0.0000',
       scan=['65,91,68,72,75,78,81,84,53,87,56,59,62'],
       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_T26t01.J000415+623000_P7471v1_2018_01_01T18_46_59.683/working/VLASS1.1.sb34459038.eb34579917.58037.11761364584_split.ms'],
       imagename='VLASS1.1.ql.T26t01.J000415+623000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='2650,2651,2652,2653,2654,2655,2656,2657,2658,2659,2660,2661,2662,2663,2664,2665,2666,2667,2668,2669,2670,2671,2672,2673,2674,2677,2678,2679,2680,2681,2682,2683,2684,2685,2686,2687,2688,2689,2690,2691,2692,2693,2694,2695,2696,2697,2698,2699,2700,2701,3019,3020,3021,3022,3023,3024,3025,3026,3027,3028,3029,3030,3031,3032,3033,3034,3035,3036,3037,3038,3039,3040,3041,3042,3043,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,3391,3392,3393,3394,3395,3396,3397,3398,3399,3400,3401,3402,3403,3404,3405,3406,3407,3408,3409,3410,3411,3412,3413,3414,3415,3418,3419,3420,3421,3422,3423,3424,3425,3426,3427,3428,3429,3430,3431,3432,3433,3434,3435,3436,3437,3438,3439,3440,3441,3442,3443,3765,3766,3767,3768,3769,3770,3771,3772,3773,3774,3775,3776,3777,3778,3779,3780,3781,3782,3783,3784,3785,3786,3787,3788,3789,3790,3793,3794,3795,3796,3797,3798,3799,3800,3801,3802,3803,3804,3805,3806,3807,3808,3809,3810,3811,3812,3813,3814,3815,3816,3817,3818,4143,4144,4145,4146,4147,4148,4149,4150,4151,4152,4153,4154,4155,4156,4157,4158,4159,4160,4161,4162,4163,4164,4165,4166,4167,4168,4171,4172,4173,4174,4175,4176,4177,4178,4179,4180,4181,4182,4183,4184,4185,4186,4187,4188,4189,4190,4191,4192,4193,4194,4195,4196,4524,4525,4526,4527,4528,4529,4530,4531,4532,4533,4534,4535,4536,4537,4538,4539,4540,4541,4542,4543,4544,4545,4546,4547,4548,4549,4552,4553,4554,4555,4556,4557,4558,4559,4560,4561,4562,4563,4564,4565,4566,4567,4568,4569,4570,4571,4572,4573,4574,4575,4576,4577,4908,4909,4910,4911,4912,4913,4914,4915,4916,4917,4918,4919,4920,4921,4922,4923,4924,4925,4926,4927,4928,4929,4930,4931,4932,4933',
       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 00:04:15.607 +62.30.0.0000',
       scan=['65,91,68,72,75,78,81,84,53,87,56,59,62'],
       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_T26t01.J000415+623000_P7471v1_2018_01_01T18_46_59.683/working/VLASS1.1.sb34459038.eb34579917.58037.11761364584_split.ms'],
       imagename='VLASS1.1.ql.T26t01.J000415+623000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00056110961976155518, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='2650,2651,2652,2653,2654,2655,2656,2657,2658,2659,2660,2661,2662,2663,2664,2665,2666,2667,2668,2669,2670,2671,2672,2673,2674,2677,2678,2679,2680,2681,2682,2683,2684,2685,2686,2687,2688,2689,2690,2691,2692,2693,2694,2695,2696,2697,2698,2699,2700,2701,3019,3020,3021,3022,3023,3024,3025,3026,3027,3028,3029,3030,3031,3032,3033,3034,3035,3036,3037,3038,3039,3040,3041,3042,3043,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,3391,3392,3393,3394,3395,3396,3397,3398,3399,3400,3401,3402,3403,3404,3405,3406,3407,3408,3409,3410,3411,3412,3413,3414,3415,3418,3419,3420,3421,3422,3423,3424,3425,3426,3427,3428,3429,3430,3431,3432,3433,3434,3435,3436,3437,3438,3439,3440,3441,3442,3443,3765,3766,3767,3768,3769,3770,3771,3772,3773,3774,3775,3776,3777,3778,3779,3780,3781,3782,3783,3784,3785,3786,3787,3788,3789,3790,3793,3794,3795,3796,3797,3798,3799,3800,3801,3802,3803,3804,3805,3806,3807,3808,3809,3810,3811,3812,3813,3814,3815,3816,3817,3818,4143,4144,4145,4146,4147,4148,4149,4150,4151,4152,4153,4154,4155,4156,4157,4158,4159,4160,4161,4162,4163,4164,4165,4166,4167,4168,4171,4172,4173,4174,4175,4176,4177,4178,4179,4180,4181,4182,4183,4184,4185,4186,4187,4188,4189,4190,4191,4192,4193,4194,4195,4196,4524,4525,4526,4527,4528,4529,4530,4531,4532,4533,4534,4535,4536,4537,4538,4539,4540,4541,4542,4543,4544,4545,4546,4547,4548,4549,4552,4553,4554,4555,4556,4557,4558,4559,4560,4561,4562,4563,4564,4565,4566,4567,4568,4569,4570,4571,4572,4573,4574,4575,4576,4577,4908,4909,4910,4911,4912,4913,4914,4915,4916,4917,4918,4919,4920,4921,4922,4923,4924,4925,4926,4927,4928,4929,4930,4931,4932,4933',
       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.T26t01.J000415+623000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T26t01.J000415+623000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T26t01.J000415+623000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#