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

# hif_editimlist(parameter_file='parameter.17039.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_T32t01.J032041+863000_P17039v1_2018_01_18T23_18_35.820/working/VLASS1.1.sb34666958.eb34700016.58074.035491539355_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T32t01.J032041+863000_P17039v1_2018_01_18T23_18_35.820/rawdata/VLASS1.1.sb34666958.eb34700016.58074.035491539355.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='4022,4020,4021,4149,4148,4145,4144,4147,4146,4141,4140,4143,4142,3999,3998,3997,3996,3995,3994,3993,3992,3658,3659,3665,2289,3084,4655,4654,4657,4656,4651,4650,4653,4652,3672,3673,3670,3671,4659,4658,3674,3675,3650,3651,2381,2380,2268,2269,2385,2384,2387,2386,2389,2263,2266,2267,2264,2265,2769,2768,2763,2284,2767,2766,2765,2285,3092,3093,3090,3091,3096,3097,3094,3095,4672,3519,3518,3517,3516,3515,3514,2680,3207,3206,3205,3204,3203,3202,3201,3200,3209,3208,3067,3066,3069,3068,4505,4504,4507,4506,4501,4500,4503,4502,4509,4508,2088,2089,4668,4669,2082,2083,4666,2081,2086,2087,2084,2085,2659,2658,2653,2652,2651,2650,2657,2656,2655,2654,4138,4139,4131,4132,4133,4134,4135,4136,4137,2288,3664,3667,3666,3661,3660,3663,3662,2280,2281,2282,2283,3669,3668,2286,2287,2275,2274,2277,2276,2271,2270,2273,2272,2279,2278,2378,2379,2374,2375,2376,2377,2370,2371,2372,2373,4530,3089,3088,3085,2111,3087,3086,3081,3080,3083,3082,3544,4150,3526,2110,3524,4151,3522,3523,3520,3521,3540,3528,3529,3541,3543,4664,4665,4011,4159,4660,4661,4662,4663,4008,4009,4161,4160,4000,4001,4002,4003,4004,4005,4006,4007,3527,3525,2095,2094,2097,2096,2091,2090,2093,2092,3654,3655,3656,3657,2099,2098,3652,3653,4678,4677,4676,4675,4674,4673,2383,4671,4670,2382,2388,4499,2792,2793,2790,2791,2794,2666,2667,2664,2665,2662,2663,2660,2661,2668,2669,2764,2293,2292,2291,2290,2369,2368,4667,4529,4528,4527,4526,4525,4524,4523,4522,4521,4520,2109,2108,2103,2102,2101,2100,2107,2106,2105,2104,3531,3530,3533,3532,3535,3534,3537,3536,3539,3538,3189,4152,4153,4019,4018,4156,4157,4154,4155,4013,4012,4158,4010,4017,4016,4015,4014,3649,3648,3647,3646,3645,4648,4649,4647,2396,2397,2394,2395,2392,2393,2390,2391,2398,2399,2778,2779,2770,2771,2772,2773,2774,2775,2776,2777,2785,2784,2787,2786,2781,2780,2783,2782,3542,2789,2788,2671,2670,2673,2672,2675,2674,2677,2676,2679,2678,3214,3215,3216,3217,3210,3211,3212,3213,3218,3219,3070,3071,3072,3073,3074,3075,3076,3077,3078,3079,4518,4519,4512,4513,4510,4511,4516,4517,4514,4515,2649,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 03:20:41.893 +86.30.0.0000',
       scan=['99,102,105,76,111,80,108,83,86,89,92,115,95'],
       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_T32t01.J032041+863000_P17039v1_2018_01_18T23_18_35.820/working/VLASS1.1.sb34666958.eb34700016.58074.035491539355_split.ms'],
       imagename='VLASS1.1.ql.T32t01.J032041+863000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='2081,2082,2083,2084,2085,2086,2087,2088,2089,2090,2091,2092,2093,2094,2095,2096,2097,2098,2099,2100,2101,2102,2103,2104,2105,2106,2107,2108,2109,2110,2111,2263,2264,2265,2266,2267,2268,2269,2270,2271,2272,2273,2274,2275,2276,2277,2278,2279,2280,2281,2282,2283,2284,2285,2286,2287,2288,2289,2290,2291,2292,2293,2368,2369,2370,2371,2372,2373,2374,2375,2376,2377,2378,2379,2380,2381,2382,2383,2384,2385,2386,2387,2388,2389,2390,2391,2392,2393,2394,2395,2396,2397,2398,2399,2649,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,2675,2676,2677,2678,2679,2680,2763,2764,2765,2766,2767,2768,2769,2770,2771,2772,2773,2774,2775,2776,2777,2778,2779,2780,2781,2782,2783,2784,2785,2786,2787,2788,2789,2790,2791,2792,2793,2794,3066,3067,3068,3069,3070,3071,3072,3073,3074,3075,3076,3077,3078,3079,3080,3081,3082,3083,3084,3085,3086,3087,3088,3089,3090,3091,3092,3093,3094,3095,3096,3097,3189,3190,3191,3192,3193,3194,3195,3196,3197,3198,3199,3200,3201,3202,3203,3204,3205,3206,3207,3208,3209,3210,3211,3212,3213,3214,3215,3216,3217,3218,3219,3514,3515,3516,3517,3518,3519,3520,3521,3522,3523,3524,3525,3526,3527,3528,3529,3530,3531,3532,3533,3534,3535,3536,3537,3538,3539,3540,3541,3542,3543,3544,3645,3646,3647,3648,3649,3650,3651,3652,3653,3654,3655,3656,3657,3658,3659,3660,3661,3662,3663,3664,3665,3666,3667,3668,3669,3670,3671,3672,3673,3674,3675,3992,3993,3994,3995,3996,3997,3998,3999,4000,4001,4002,4003,4004,4005,4006,4007,4008,4009,4010,4011,4012,4013,4014,4015,4016,4017,4018,4019,4020,4021,4022,4131,4132,4133,4134,4135,4136,4137,4138,4139,4140,4141,4142,4143,4144,4145,4146,4147,4148,4149,4150,4151,4152,4153,4154,4155,4156,4157,4158,4159,4160,4161,4499,4500,4501,4502,4503,4504,4505,4506,4507,4508,4509,4510,4511,4512,4513,4514,4515,4516,4517,4518,4519,4520,4521,4522,4523,4524,4525,4526,4527,4528,4529,4530,4647,4648,4649,4650,4651,4652,4653,4654,4655,4656,4657,4658,4659,4660,4661,4662,4663,4664,4665,4666,4667,4668,4669,4670,4671,4672,4673,4674,4675,4676,4677,4678',
       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 03:20:41.893 +86.30.0.0000',
       scan=['99,102,105,76,111,80,108,83,86,89,92,115,95'],
       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_T32t01.J032041+863000_P17039v1_2018_01_18T23_18_35.820/working/VLASS1.1.sb34666958.eb34700016.58074.035491539355_split.ms'],
       imagename='VLASS1.1.ql.T32t01.J032041+863000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00054194105717142514, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='2081,2082,2083,2084,2085,2086,2087,2088,2089,2090,2091,2092,2093,2094,2095,2096,2097,2098,2099,2100,2101,2102,2103,2104,2105,2106,2107,2108,2109,2110,2111,2263,2264,2265,2266,2267,2268,2269,2270,2271,2272,2273,2274,2275,2276,2277,2278,2279,2280,2281,2282,2283,2284,2285,2286,2287,2288,2289,2290,2291,2292,2293,2368,2369,2370,2371,2372,2373,2374,2375,2376,2377,2378,2379,2380,2381,2382,2383,2384,2385,2386,2387,2388,2389,2390,2391,2392,2393,2394,2395,2396,2397,2398,2399,2649,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,2675,2676,2677,2678,2679,2680,2763,2764,2765,2766,2767,2768,2769,2770,2771,2772,2773,2774,2775,2776,2777,2778,2779,2780,2781,2782,2783,2784,2785,2786,2787,2788,2789,2790,2791,2792,2793,2794,3066,3067,3068,3069,3070,3071,3072,3073,3074,3075,3076,3077,3078,3079,3080,3081,3082,3083,3084,3085,3086,3087,3088,3089,3090,3091,3092,3093,3094,3095,3096,3097,3189,3190,3191,3192,3193,3194,3195,3196,3197,3198,3199,3200,3201,3202,3203,3204,3205,3206,3207,3208,3209,3210,3211,3212,3213,3214,3215,3216,3217,3218,3219,3514,3515,3516,3517,3518,3519,3520,3521,3522,3523,3524,3525,3526,3527,3528,3529,3530,3531,3532,3533,3534,3535,3536,3537,3538,3539,3540,3541,3542,3543,3544,3645,3646,3647,3648,3649,3650,3651,3652,3653,3654,3655,3656,3657,3658,3659,3660,3661,3662,3663,3664,3665,3666,3667,3668,3669,3670,3671,3672,3673,3674,3675,3992,3993,3994,3995,3996,3997,3998,3999,4000,4001,4002,4003,4004,4005,4006,4007,4008,4009,4010,4011,4012,4013,4014,4015,4016,4017,4018,4019,4020,4021,4022,4131,4132,4133,4134,4135,4136,4137,4138,4139,4140,4141,4142,4143,4144,4145,4146,4147,4148,4149,4150,4151,4152,4153,4154,4155,4156,4157,4158,4159,4160,4161,4499,4500,4501,4502,4503,4504,4505,4506,4507,4508,4509,4510,4511,4512,4513,4514,4515,4516,4517,4518,4519,4520,4521,4522,4523,4524,4525,4526,4527,4528,4529,4530,4647,4648,4649,4650,4651,4652,4653,4654,4655,4656,4657,4658,4659,4660,4661,4662,4663,4664,4665,4666,4667,4668,4669,4670,4671,4672,4673,4674,4675,4676,4677,4678',
       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.T32t01.J032041+863000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T32t01.J032041+863000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T32t01.J032041+863000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#