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

# hif_editimlist(parameter_file='parameter.4945.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_T24t18.J172349+543000_P4945v1_2017_09_21T17_27_11.805/working/VLASS1.1.sb34382472.eb34389458.58009.054277766205_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T24t18.J172349+543000_P4945v1_2017_09_21T17_27_11.805/rawdata/VLASS1.1.sb34382472.eb34389458.58009.054277766205.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='4448,4449,3557,3556,3555,3554,3553,3552,3551,3550,3559,3558,3856,3857,3858,3859,2859,2858,2851,2850,2853,2852,2855,2854,2857,2856,3333,3332,3331,3330,3337,3336,3335,3334,3203,3339,3338,3202,2829,3208,2989,2988,2983,2982,2981,2980,2987,2986,2985,2984,4215,4214,4217,4216,4219,4218,2634,3869,2619,2618,4101,4100,3868,2639,2638,3863,4455,3862,2637,2636,2631,2630,2633,2632,4589,3207,3206,3205,3204,4079,4078,3201,3200,4075,4074,4077,4076,4071,4073,4072,4588,3348,3349,3342,3343,3340,3341,3346,3347,3344,3345,4460,4461,4228,4229,4224,4225,4226,4227,4220,4221,4222,4223,4592,4593,4590,4591,4597,4594,4595,3562,3563,3560,3561,3566,3564,3565,4583,3744,3742,3743,3740,3741,2839,2838,2837,2836,2835,2834,2833,2832,2831,2830,4435,4434,4437,4436,4431,4433,4432,4439,4438,2622,2623,2620,2621,2626,2627,2624,2625,3870,3871,3872,2629,3874,3875,3876,3877,3885,3884,3886,3881,3880,3883,3882,3878,3879,2628,3873,2978,2979,3351,3350,3352,2972,2973,2970,2971,2976,2977,2974,2975,3209,3715,3714,3717,3716,3713,3719,3718,4574,4575,4576,4577,4578,4579,2488,2489,2486,2487,2484,2485,2482,2483,2480,2481,2969,2968,4239,4238,4237,4236,4235,4234,4233,4232,4231,4230,4454,4440,4441,4442,4443,3728,3729,4446,4447,3724,3725,3726,3727,3720,3721,3722,3723,4080,4081,4082,4083,4084,4085,4086,4087,4088,4089,3861,3860,3867,3866,3865,3864,3535,3537,3536,3539,3538,4459,4458,3184,3185,3186,3187,4596,3181,3182,3183,4598,4599,3188,3189,3544,3545,3546,3547,3540,3541,3542,3543,3548,3549,2848,2849,2846,2847,2844,2845,2842,2843,2840,2841,3322,3323,3324,3325,3326,3327,3328,3329,2998,2498,2499,2999,2491,2490,2493,2492,2495,2494,2996,2496,2990,2991,2992,2648,2993,2994,2479,2995,2497,2997,4445,2642,2643,3210,3211,4604,4605,4602,4603,4600,4601,4244,4242,4243,4240,4241,4453,4452,4451,4450,4457,4456,3739,3738,3737,3736,3735,3734,3733,3732,3731,3730,2635,2509,2508,2507,2506,2505,2504,2503,2502,2501,2500,4093,4092,4091,4090,4097,4096,4095,4094,2640,2641,4099,4098,2644,2645,2646,2647,4584,4586,3197,3196,3195,3194,3193,3192,3191,3190,4585,4444,4587,4582,4581,4580,3199,3198',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 17:23:49.093 +54.30.0.0000',
       scan=['64,67,70,74,77,80,83,52,86,55,89,58,61'],
       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_T24t18.J172349+543000_P4945v1_2017_09_21T17_27_11.805/working/VLASS1.1.sb34382472.eb34389458.58009.054277766205_split.ms'],
       imagename='VLASS1.1.ql.T24t18.J172349+543000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='2479,2480,2481,2482,2483,2484,2485,2486,2487,2488,2489,2490,2491,2492,2493,2494,2495,2496,2497,2498,2499,2500,2501,2502,2503,2504,2505,2506,2507,2508,2509,2618,2619,2620,2621,2622,2623,2624,2625,2626,2627,2628,2629,2630,2631,2632,2633,2634,2635,2636,2637,2638,2639,2640,2641,2642,2643,2644,2645,2646,2647,2648,2829,2830,2831,2832,2833,2834,2835,2836,2837,2838,2839,2840,2841,2842,2843,2844,2845,2846,2847,2848,2849,2850,2851,2852,2853,2854,2855,2856,2857,2858,2859,2968,2969,2970,2971,2972,2973,2974,2975,2976,2977,2978,2979,2980,2981,2982,2983,2984,2985,2986,2987,2988,2989,2990,2991,2992,2993,2994,2995,2996,2997,2998,2999,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,3322,3323,3324,3325,3326,3327,3328,3329,3330,3331,3332,3333,3334,3335,3336,3337,3338,3339,3340,3341,3342,3343,3344,3345,3346,3347,3348,3349,3350,3351,3352,3535,3536,3537,3538,3539,3540,3541,3542,3543,3544,3545,3546,3547,3548,3549,3550,3551,3552,3553,3554,3555,3556,3557,3558,3559,3560,3561,3562,3563,3564,3565,3566,3713,3714,3715,3716,3717,3718,3719,3720,3721,3722,3723,3724,3725,3726,3727,3728,3729,3730,3731,3732,3733,3734,3735,3736,3737,3738,3739,3740,3741,3742,3743,3744,3856,3857,3858,3859,3860,3861,3862,3863,3864,3865,3866,3867,3868,3869,3870,3871,3872,3873,3874,3875,3876,3877,3878,3879,3880,3881,3882,3883,3884,3885,3886,4071,4072,4073,4074,4075,4076,4077,4078,4079,4080,4081,4082,4083,4084,4085,4086,4087,4088,4089,4090,4091,4092,4093,4094,4095,4096,4097,4098,4099,4100,4101,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,4240,4241,4242,4243,4244,4431,4432,4433,4434,4435,4436,4437,4438,4439,4440,4441,4442,4443,4444,4445,4446,4447,4448,4449,4450,4451,4452,4453,4454,4455,4456,4457,4458,4459,4460,4461,4574,4575,4576,4577,4578,4579,4580,4581,4582,4583,4584,4585,4586,4587,4588,4589,4590,4591,4592,4593,4594,4595,4596,4597,4598,4599,4600,4601,4602,4603,4604,4605',
       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 17:23:49.093 +54.30.0.0000',
       scan=['64,67,70,74,77,80,83,52,86,55,89,58,61'],
       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_T24t18.J172349+543000_P4945v1_2017_09_21T17_27_11.805/working/VLASS1.1.sb34382472.eb34389458.58009.054277766205_split.ms'],
       imagename='VLASS1.1.ql.T24t18.J172349+543000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00054305247357478944, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='2479,2480,2481,2482,2483,2484,2485,2486,2487,2488,2489,2490,2491,2492,2493,2494,2495,2496,2497,2498,2499,2500,2501,2502,2503,2504,2505,2506,2507,2508,2509,2618,2619,2620,2621,2622,2623,2624,2625,2626,2627,2628,2629,2630,2631,2632,2633,2634,2635,2636,2637,2638,2639,2640,2641,2642,2643,2644,2645,2646,2647,2648,2829,2830,2831,2832,2833,2834,2835,2836,2837,2838,2839,2840,2841,2842,2843,2844,2845,2846,2847,2848,2849,2850,2851,2852,2853,2854,2855,2856,2857,2858,2859,2968,2969,2970,2971,2972,2973,2974,2975,2976,2977,2978,2979,2980,2981,2982,2983,2984,2985,2986,2987,2988,2989,2990,2991,2992,2993,2994,2995,2996,2997,2998,2999,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,3322,3323,3324,3325,3326,3327,3328,3329,3330,3331,3332,3333,3334,3335,3336,3337,3338,3339,3340,3341,3342,3343,3344,3345,3346,3347,3348,3349,3350,3351,3352,3535,3536,3537,3538,3539,3540,3541,3542,3543,3544,3545,3546,3547,3548,3549,3550,3551,3552,3553,3554,3555,3556,3557,3558,3559,3560,3561,3562,3563,3564,3565,3566,3713,3714,3715,3716,3717,3718,3719,3720,3721,3722,3723,3724,3725,3726,3727,3728,3729,3730,3731,3732,3733,3734,3735,3736,3737,3738,3739,3740,3741,3742,3743,3744,3856,3857,3858,3859,3860,3861,3862,3863,3864,3865,3866,3867,3868,3869,3870,3871,3872,3873,3874,3875,3876,3877,3878,3879,3880,3881,3882,3883,3884,3885,3886,4071,4072,4073,4074,4075,4076,4077,4078,4079,4080,4081,4082,4083,4084,4085,4086,4087,4088,4089,4090,4091,4092,4093,4094,4095,4096,4097,4098,4099,4100,4101,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,4240,4241,4242,4243,4244,4431,4432,4433,4434,4435,4436,4437,4438,4439,4440,4441,4442,4443,4444,4445,4446,4447,4448,4449,4450,4451,4452,4453,4454,4455,4456,4457,4458,4459,4460,4461,4574,4575,4576,4577,4578,4579,4580,4581,4582,4583,4584,4585,4586,4587,4588,4589,4590,4591,4592,4593,4594,4595,4596,4597,4598,4599,4600,4601,4602,4603,4604,4605',
       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.T24t18.J172349+543000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T24t18.J172349+543000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T24t18.J172349+543000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#