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

# hif_editimlist(parameter_file='parameter.4472.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_T27t13.J204415+663000_P4472v2_2017_09_14T17_37_20.602/working/VLASS1.1.sb34346984.eb34356943.58004.0160307176_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T27t13.J204415+663000_P4472v2_2017_09_14T17_37_20.602/rawdata/VLASS1.1.sb34346984.eb34356943.58004.0160307176.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='11389,2911,2912,2913,2914,2915,2916,2917,2918,2919,3333,3332,3331,3330,3337,3336,3335,3334,3339,3338,10209,10199,11003,11002,11001,11000,11006,11005,11004,4101,4100,4103,4102,4105,4104,4107,4106,4109,4108,3296,3297,4079,4078,4075,4074,4077,4076,4073,4072,2949,2948,2947,2946,2945,2944,2943,2942,2941,2940,3340,3341,9442,9441,9440,4505,4504,3708,4506,4501,4500,4503,4502,3702,4467,3700,3701,4509,4508,10206,4912,4913,4910,4911,4507,10205,10202,10203,10200,10201,4868,4869,4866,4867,9061,9060,4484,4485,4486,4487,4480,4481,4482,4483,10598,10599,4488,4489,3308,3309,3306,3307,3304,3305,3302,3303,3300,3301,10204,2578,2579,2576,2577,2574,2575,2572,2573,2570,2571,3709,4468,4469,3703,3698,3699,3690,3691,3692,3693,3694,3695,3696,3697,9058,9059,9050,9051,9054,9055,9056,9057,3715,3714,3717,3716,3711,3710,3713,3712,3719,3718,4479,4478,4471,4470,4473,4472,4475,4474,4477,4476,4909,4908,4901,4900,4903,4902,4905,4904,4907,4906,4888,4880,4881,4882,4883,4884,4885,4886,4887,9825,9824,9827,9826,9821,9820,9823,9822,4875,4874,4877,4876,4871,4870,4873,4872,4879,4878,10216,10211,10210,10213,10212,4497,4496,4495,4494,4493,4492,4499,4498,2938,2939,2936,2937,2932,2933,2930,2931,3319,3311,3310,3313,3312,3315,3314,3316,10989,10215,2569,2568,2561,2560,2563,2562,2565,2564,2567,2566,10214,3728,3724,3725,3726,3727,3720,3721,3722,3723,3689,3688,3683,3682,3681,3687,3686,3685,3684,9438,9439,9436,9437,9435,9432,9430,9431,4080,4081,4082,4083,4084,4085,4086,4087,4088,4089,9814,9815,9816,9817,9810,9811,9812,9813,9049,9048,9047,9046,9045,9044,11396,11395,11394,11393,11392,11391,11390,3320,3321,3322,3323,3324,3325,3326,3327,3328,3329,2557,2558,2559,3298,3299,4899,4898,4893,4892,4891,4897,4896,4895,4894,3294,3295,10592,10593,10602,10603,10606,10607,10604,10605,10596,10608,10609,10597,10594,10595,4116,4117,4114,4115,4112,4113,4110,4111,4118,4119,2929,2928,2921,2920,2923,2922,2925,2924,2927,2926,10999,10996,10994,10995,10992,10993,10990,10991,2958,2954,2955,2956,2957,2950,2951,2952,2953,3706,4512,4513,4510,4511,4514,9429,9428,3707,9425,9427,9426,4093,4092,4091,4090,4097,4094,4099,4098',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 20:44:15.281 +66.30.0.0000',
       scan=['156,160,163,166,169,172,175,179,52,182,55,185,58,188,61,191,65,194,68,71,74,77,80,84,87,90'],
       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_T27t13.J204415+663000_P4472v2_2017_09_14T17_37_20.602/working/VLASS1.1.sb34346984.eb34356943.58004.0160307176_split.ms'],
       imagename='VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='2557,2558,2559,2560,2561,2562,2563,2564,2565,2566,2567,2568,2569,2570,2571,2572,2573,2574,2575,2576,2577,2578,2579,2911,2912,2913,2914,2915,2916,2917,2918,2919,2920,2921,2922,2923,2924,2925,2926,2927,2928,2929,2930,2931,2932,2933,2936,2937,2938,2939,2940,2941,2942,2943,2944,2945,2946,2947,2948,2949,2950,2951,2952,2953,2954,2955,2956,2957,2958,3294,3295,3296,3297,3298,3299,3300,3301,3302,3303,3304,3305,3306,3307,3308,3309,3310,3311,3312,3313,3314,3315,3316,3319,3320,3321,3322,3323,3324,3325,3326,3327,3328,3329,3330,3331,3332,3333,3334,3335,3336,3337,3338,3339,3340,3341,3681,3682,3683,3684,3685,3686,3687,3688,3689,3690,3691,3692,3693,3694,3695,3696,3697,3698,3699,3700,3701,3702,3703,3706,3707,3708,3709,3710,3711,3712,3713,3714,3715,3716,3717,3718,3719,3720,3721,3722,3723,3724,3725,3726,3727,3728,4072,4073,4074,4075,4076,4077,4078,4079,4080,4081,4082,4083,4084,4085,4086,4087,4088,4089,4090,4091,4092,4093,4094,4097,4098,4099,4100,4101,4102,4103,4104,4105,4106,4107,4108,4109,4110,4111,4112,4113,4114,4115,4116,4117,4118,4119,4467,4468,4469,4470,4471,4472,4473,4474,4475,4476,4477,4478,4479,4480,4481,4482,4483,4484,4485,4486,4487,4488,4489,4492,4493,4494,4495,4496,4497,4498,4499,4500,4501,4502,4503,4504,4505,4506,4507,4508,4509,4510,4511,4512,4513,4514,4866,4867,4868,4869,4870,4871,4872,4873,4874,4875,4876,4877,4878,4879,4880,4881,4882,4883,4884,4885,4886,4887,4888,4891,4892,4893,4894,4895,4896,4897,4898,4899,4900,4901,4902,4903,4904,4905,4906,4907,4908,4909,4910,4911,4912,4913,9044,9045,9046,9047,9048,9049,9050,9051,9054,9055,9056,9057,9058,9059,9060,9061,9425,9426,9427,9428,9429,9430,9431,9432,9435,9436,9437,9438,9439,9440,9441,9442,9810,9811,9812,9813,9814,9815,9816,9817,9820,9821,9822,9823,9824,9825,9826,9827,10199,10200,10201,10202,10203,10204,10205,10206,10209,10210,10211,10212,10213,10214,10215,10216,10592,10593,10594,10595,10596,10597,10598,10599,10602,10603,10604,10605,10606,10607,10608,10609,10989,10990,10991,10992,10993,10994,10995,10996,10999,11000,11001,11002,11003,11004,11005,11006,11389,11390,11391,11392,11393,11394,11395,11396',
       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 20:44:15.281 +66.30.0.0000',
       scan=['156,160,163,166,169,172,175,179,52,182,55,185,58,188,61,191,65,194,68,71,74,77,80,84,87,90'],
       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_T27t13.J204415+663000_P4472v2_2017_09_14T17_37_20.602/working/VLASS1.1.sb34346984.eb34356943.58004.0160307176_split.ms'],
       imagename='VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.iter1',
       nterms=2, threshold=0.00048821114961647254, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='2557,2558,2559,2560,2561,2562,2563,2564,2565,2566,2567,2568,2569,2570,2571,2572,2573,2574,2575,2576,2577,2578,2579,2911,2912,2913,2914,2915,2916,2917,2918,2919,2920,2921,2922,2923,2924,2925,2926,2927,2928,2929,2930,2931,2932,2933,2936,2937,2938,2939,2940,2941,2942,2943,2944,2945,2946,2947,2948,2949,2950,2951,2952,2953,2954,2955,2956,2957,2958,3294,3295,3296,3297,3298,3299,3300,3301,3302,3303,3304,3305,3306,3307,3308,3309,3310,3311,3312,3313,3314,3315,3316,3319,3320,3321,3322,3323,3324,3325,3326,3327,3328,3329,3330,3331,3332,3333,3334,3335,3336,3337,3338,3339,3340,3341,3681,3682,3683,3684,3685,3686,3687,3688,3689,3690,3691,3692,3693,3694,3695,3696,3697,3698,3699,3700,3701,3702,3703,3706,3707,3708,3709,3710,3711,3712,3713,3714,3715,3716,3717,3718,3719,3720,3721,3722,3723,3724,3725,3726,3727,3728,4072,4073,4074,4075,4076,4077,4078,4079,4080,4081,4082,4083,4084,4085,4086,4087,4088,4089,4090,4091,4092,4093,4094,4097,4098,4099,4100,4101,4102,4103,4104,4105,4106,4107,4108,4109,4110,4111,4112,4113,4114,4115,4116,4117,4118,4119,4467,4468,4469,4470,4471,4472,4473,4474,4475,4476,4477,4478,4479,4480,4481,4482,4483,4484,4485,4486,4487,4488,4489,4492,4493,4494,4495,4496,4497,4498,4499,4500,4501,4502,4503,4504,4505,4506,4507,4508,4509,4510,4511,4512,4513,4514,4866,4867,4868,4869,4870,4871,4872,4873,4874,4875,4876,4877,4878,4879,4880,4881,4882,4883,4884,4885,4886,4887,4888,4891,4892,4893,4894,4895,4896,4897,4898,4899,4900,4901,4902,4903,4904,4905,4906,4907,4908,4909,4910,4911,4912,4913,9044,9045,9046,9047,9048,9049,9050,9051,9054,9055,9056,9057,9058,9059,9060,9061,9425,9426,9427,9428,9429,9430,9431,9432,9435,9436,9437,9438,9439,9440,9441,9442,9810,9811,9812,9813,9814,9815,9816,9817,9820,9821,9822,9823,9824,9825,9826,9827,10199,10200,10201,10202,10203,10204,10205,10206,10209,10210,10211,10212,10213,10214,10215,10216,10592,10593,10594,10595,10596,10597,10598,10599,10602,10603,10604,10605,10606,10607,10608,10609,10989,10990,10991,10992,10993,10994,10995,10996,10999,11000,11001,11002,11003,11004,11005,11006,11389,11390,11391,11392,11393,11394,11395,11396',
       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.T27t13.J204415+663000.10.2048.v2.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.iter1.pb.tt0')

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

# hifv_exportdata(pipelinemode="automatic")
#
exportfits(dropstokes=False, minpix=0,
           fitsimage='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T27t13.J204415+663000_P4472v2_2017_09_14T17_37_20.602/products/VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.tt0.fits',
           imagename='VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.iter1.image.tt0',
           overwrite=True, bitpix=-32, maxpix=-1, velocity=False, optical=False,
           stokeslast=True)
exportfits(dropstokes=False, minpix=0,
           fitsimage='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T27t13.J204415+663000_P4472v2_2017_09_14T17_37_20.602/products/VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.tt1.fits',
           imagename='VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.iter1.image.tt1',
           overwrite=True, bitpix=-32, maxpix=-1, velocity=False, optical=False,
           stokeslast=True)
exportfits(dropstokes=False, minpix=0,
           fitsimage='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T27t13.J204415+663000_P4472v2_2017_09_14T17_37_20.602/products/VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.alpha.fits',
           imagename='VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.iter1.alpha',
           overwrite=True, bitpix=-32, maxpix=-1, velocity=False, optical=False,
           stokeslast=True)
exportfits(dropstokes=False, minpix=0,
           fitsimage='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T27t13.J204415+663000_P4472v2_2017_09_14T17_37_20.602/products/VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.alpha.error.fits',
           imagename='VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.iter1.alpha.error',
           overwrite=True, bitpix=-32, maxpix=-1, velocity=False, optical=False,
           stokeslast=True)
exportfits(dropstokes=False, minpix=0,
           fitsimage='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T27t13.J204415+663000_P4472v2_2017_09_14T17_37_20.602/products/VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.pb.tt0.fits',
           imagename='VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.iter1.pb.tt0',
           overwrite=True, bitpix=-32, maxpix=-1, velocity=False, optical=False,
           stokeslast=True)
exportfits(dropstokes=False, minpix=0,
           fitsimage='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T27t13.J204415+663000_P4472v2_2017_09_14T17_37_20.602/products/VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.mask.fits',
           imagename='VLASS1.1.ql.T27t13.J204415+663000.10.2048.v2.I.iter1.mask',
           overwrite=True, bitpix=-32, maxpix=-1, velocity=False, optical=False,
           stokeslast=True)

# hifv_exportvlassdata(pipelinemode="automatic")
#