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.14033.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.J120324+543000_P14033v1_2018_02_06T18_21_03.332/working/VLASS1.1.sb34541777.eb34701358.58077.63418230324_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T24t13.J120324+543000_P14033v1_2018_02_06T18_21_03.332/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='3967,3433,2716,2717,2714,2715,3434,2718,2719,3435,3422,3436,3974,3437,3418,3419,3971,4679,4678,4677,4676,4363,4674,4673,3964,3991,2406,2407,2404,2405,2402,2403,2400,2401,3997,3975,3995,3994,3993,3992,2408,2409,3977,4336,4337,4334,3965,4332,4333,4330,3420,3976,3649,3648,3647,3646,3645,3644,3643,3642,3641,3640,4680,3962,4313,2396,2410,2394,2395,2415,2414,2416,2398,2399,4329,3629,4335,4691,3638,3628,4321,3639,4323,4322,3966,4312,4327,4326,3978,4000,4688,4314,3981,3989,4364,2411,4344,3470,3471,3961,2397,3425,4315,4317,2413,3953,2412,2392,4001,2393,3968,3959,3958,3960,3424,2763,2762,2761,2760,3957,3956,2765,2764,3996,4358,3627,4354,4355,4356,3626,4350,4351,4352,4353,3467,3466,3465,3464,3463,3462,3461,3460,4003,3469,3468,3092,3093,3096,3097,3094,3095,3098,3099,4675,3954,4359,4365,3969,4328,4004,4349,4348,4347,4346,4345,4002,4343,4342,4341,4340,3452,3453,3450,3451,3456,3457,3454,3455,4005,3458,3459,4690,3089,3088,4361,3637,3085,3084,3087,3086,3081,3080,3083,3082,3999,4006,3632,3998,4362,4324,3070,3071,3072,3073,3074,3075,3076,3077,3078,3079,3631,4686,4687,4684,4685,3449,3448,3955,4681,4682,3447,3446,3441,3440,3443,3442,2756,2757,2754,2755,2752,2753,2750,2751,4357,3963,3990,4325,3421,2758,2759,2734,2735,2736,2737,2730,2731,2732,2733,3636,3427,3634,3635,2738,3633,3630,3426,4698,3067,3066,3065,3064,3069,3068,3438,3439,4693,4692,3108,3109,4697,4696,3104,3105,3106,3107,3100,3101,3102,3103,2741,3970,2743,2742,2745,2744,2747,2746,2749,2748,4689,3973,2727,2726,2725,2724,2723,2722,2721,2720,4318,4319,3625,3624,2729,2728,4316,4320,4695,4331,3972,4694,4683,3429,3428,4360,3117,3116,3115,3114,3113,3112,3111,3110,3430,3988,3423,3431,3984,3985,3986,3987,3432,3982,3983',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 12:03:24.156 +54.30.0.0000',
       scan=['65,68,71,75,78,81,84,53,87,56,90,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_T24t13.J120324+543000_P14033v1_2018_02_06T18_21_03.332/working/VLASS1.1.sb34541777.eb34701358.58077.63418230324_split.ms'],
       imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='2392,2393,2394,2395,2396,2397,2398,2399,2400,2401,2402,2403,2404,2405,2406,2407,2408,2409,2410,2411,2412,2413,2414,2415,2416,2714,2715,2716,2717,2718,2719,2720,2721,2722,2723,2724,2725,2726,2727,2728,2729,2730,2731,2732,2733,2734,2735,2736,2737,2738,2741,2742,2743,2744,2745,2746,2747,2748,2749,2750,2751,2752,2753,2754,2755,2756,2757,2758,2759,2760,2761,2762,2763,2764,2765,3064,3065,3066,3067,3068,3069,3070,3071,3072,3073,3074,3075,3076,3077,3078,3079,3080,3081,3082,3083,3084,3085,3086,3087,3088,3089,3092,3093,3094,3095,3096,3097,3098,3099,3100,3101,3102,3103,3104,3105,3106,3107,3108,3109,3110,3111,3112,3113,3114,3115,3116,3117,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,3446,3447,3448,3449,3450,3451,3452,3453,3454,3455,3456,3457,3458,3459,3460,3461,3462,3463,3464,3465,3466,3467,3468,3469,3470,3471,3624,3625,3626,3627,3628,3629,3630,3631,3632,3633,3634,3635,3636,3637,3638,3639,3640,3641,3642,3643,3644,3645,3646,3647,3648,3649,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,4006,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',
       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:24.156 +54.30.0.0000',
       scan=['65,68,71,75,78,81,84,53,87,56,90,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_T24t13.J120324+543000_P14033v1_2018_02_06T18_21_03.332/working/VLASS1.1.sb34541777.eb34701358.58077.63418230324_split.ms'],
       imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00057536866673726765, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='2392,2393,2394,2395,2396,2397,2398,2399,2400,2401,2402,2403,2404,2405,2406,2407,2408,2409,2410,2411,2412,2413,2414,2415,2416,2714,2715,2716,2717,2718,2719,2720,2721,2722,2723,2724,2725,2726,2727,2728,2729,2730,2731,2732,2733,2734,2735,2736,2737,2738,2741,2742,2743,2744,2745,2746,2747,2748,2749,2750,2751,2752,2753,2754,2755,2756,2757,2758,2759,2760,2761,2762,2763,2764,2765,3064,3065,3066,3067,3068,3069,3070,3071,3072,3073,3074,3075,3076,3077,3078,3079,3080,3081,3082,3083,3084,3085,3086,3087,3088,3089,3092,3093,3094,3095,3096,3097,3098,3099,3100,3101,3102,3103,3104,3105,3106,3107,3108,3109,3110,3111,3112,3113,3114,3115,3116,3117,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,3446,3447,3448,3449,3450,3451,3452,3453,3454,3455,3456,3457,3458,3459,3460,3461,3462,3463,3464,3465,3466,3467,3468,3469,3470,3471,3624,3625,3626,3627,3628,3629,3630,3631,3632,3633,3634,3635,3636,3637,3638,3639,3640,3641,3642,3643,3644,3645,3646,3647,3648,3649,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,4006,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',
       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.J120324+543000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T24t13.J120324+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.T24t13.J120324+543000.10.2048.v1.I.iter1.image.pbcor.tt0.rms',
      imagename='VLASS1.1.ql.T24t13.J120324+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.T24t13.J120324+543000.10.2048.v1.I.iter1.image.pbcor.tt0')
imhead(imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.image.pbcor.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.image.pbcor.tt0.subim',
           imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.image.pbcor.tt0')
imstat(imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imhead(imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.image.pbcor.tt0.rms.subim',
           imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imstat(imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.psf.tt0')
imhead(imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.psf.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.psf.tt0.subim',
           imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.psf.tt0')
imstat(imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.residual.tt0')
imhead(imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.residual.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.residual.tt0.subim',
           imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.residual.tt0')
imstat(imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.pb.tt0')
imhead(imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.pb.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.pb.tt0.subim',
           imagename='VLASS1.1.ql.T24t13.J120324+543000.10.2048.v1.I.iter1.pb.tt0')

# hifv_exportvlassdata(pipelinemode="automatic")
#