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

# hif_editimlist(parameter_file='parameter.5312.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_T25t10.J131600+603000_P5312v1_2017_12_22T16_15_37.320/working/VLASS1.1.sb34397357.eb34435512.58014.72557283565_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T25t10.J131600+603000_P5312v1_2017_12_22T16_15_37.320/rawdata/VLASS1.1.sb34397357.eb34435512.58014.72557283565.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='13834,13835,13836,13837,13830,13831,13832,13833,15417,15422,13838,13839,14601,14610,13401,15018,13403,13402,13405,13404,13407,13406,13409,13408,15026,14213,15043,15042,14611,14212,15047,15046,15045,15044,15429,15049,15048,15428,15433,13007,15039,13009,13008,13841,13840,13843,13842,13845,13844,13847,13846,13849,13848,15430,13412,13413,13410,13411,13416,13417,13414,13415,14211,14210,13418,13419,14215,14214,14217,14216,15058,15059,15440,15054,15055,15056,15057,15050,15051,15052,15053,13016,13017,13014,13015,13012,13013,13010,13011,14612,13018,13019,15421,15438,13850,13851,15008,15441,15436,15437,14208,14209,14619,14618,13429,15431,14202,13426,13425,13424,13423,13422,13421,13420,15029,15028,14198,14199,14613,15021,15020,15023,15022,15025,15024,15027,14636,13023,13022,13021,13020,13027,13026,13025,13024,13029,13028,15423,14237,14236,14235,14234,14233,14232,14231,14230,15425,15041,14239,14238,13438,13439,15040,15432,14608,14609,13430,13431,13432,13433,13434,13435,13436,13437,15009,14616,15032,15030,15031,15036,15037,15035,13038,13039,13034,13035,13036,13037,13030,13031,14220,14221,14222,14223,14226,14227,14228,14229,15014,15007,14606,14639,14638,15418,15419,14607,14633,14632,14631,14630,14637,14604,14635,14634,14605,13449,13448,13445,13444,13447,13446,13441,13440,13443,13442,13049,13048,15424,13041,13040,13043,13042,13045,13044,13047,13046,14624,14625,14626,14620,14621,14622,14623,14203,14629,15010,15011,15012,15013,13809,13808,15016,15017,13805,13804,13807,13806,13801,13800,13803,13802,14200,14602,13452,13453,13450,13451,13058,15435,13052,13053,13050,13051,13056,13057,13054,13055,14651,14650,14653,14652,14240,13799,13798,14241,15420,15439,14246,14247,14244,14245,14242,14243,13818,13819,13816,13817,13814,13815,13812,13813,13810,13811,14201,15060,14206,14615,14207,14249,14204,14654,14603,14205,14642,14643,14640,14641,14646,14647,14644,14645,14614,15038,14648,14649,13829,13828,15434,14248,13823,13822,13821,13820,13827,13826,15019,14219,14617,14218,15416,15427,15426,14251,15015,14250',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 13:16:0.000 +60.30.0.0000',
       scan=['224,227,246,230,233,237,240,243,214,249,218,252,221'],
       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_T25t10.J131600+603000_P5312v1_2017_12_22T16_15_37.320/working/VLASS1.1.sb34397357.eb34435512.58014.72557283565_split.ms'],
       imagename='VLASS1.1.ql.T25t10.J131600+603000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='13007,13008,13009,13010,13011,13012,13013,13014,13015,13016,13017,13018,13019,13020,13021,13022,13023,13024,13025,13026,13027,13028,13029,13030,13031,13034,13035,13036,13037,13038,13039,13040,13041,13042,13043,13044,13045,13046,13047,13048,13049,13050,13051,13052,13053,13054,13055,13056,13057,13058,13401,13402,13403,13404,13405,13406,13407,13408,13409,13410,13411,13412,13413,13414,13415,13416,13417,13418,13419,13420,13421,13422,13423,13424,13425,13426,13429,13430,13431,13432,13433,13434,13435,13436,13437,13438,13439,13440,13441,13442,13443,13444,13445,13446,13447,13448,13449,13450,13451,13452,13453,13798,13799,13800,13801,13802,13803,13804,13805,13806,13807,13808,13809,13810,13811,13812,13813,13814,13815,13816,13817,13818,13819,13820,13821,13822,13823,13826,13827,13828,13829,13830,13831,13832,13833,13834,13835,13836,13837,13838,13839,13840,13841,13842,13843,13844,13845,13846,13847,13848,13849,13850,13851,14198,14199,14200,14201,14202,14203,14204,14205,14206,14207,14208,14209,14210,14211,14212,14213,14214,14215,14216,14217,14218,14219,14220,14221,14222,14223,14226,14227,14228,14229,14230,14231,14232,14233,14234,14235,14236,14237,14238,14239,14240,14241,14242,14243,14244,14245,14246,14247,14248,14249,14250,14251,14601,14602,14603,14604,14605,14606,14607,14608,14609,14610,14611,14612,14613,14614,14615,14616,14617,14618,14619,14620,14621,14622,14623,14624,14625,14626,14629,14630,14631,14632,14633,14634,14635,14636,14637,14638,14639,14640,14641,14642,14643,14644,14645,14646,14647,14648,14649,14650,14651,14652,14653,14654,15007,15008,15009,15010,15011,15012,15013,15014,15015,15016,15017,15018,15019,15020,15021,15022,15023,15024,15025,15026,15027,15028,15029,15030,15031,15032,15035,15036,15037,15038,15039,15040,15041,15042,15043,15044,15045,15046,15047,15048,15049,15050,15051,15052,15053,15054,15055,15056,15057,15058,15059,15060,15416,15417,15418,15419,15420,15421,15422,15423,15424,15425,15426,15427,15428,15429,15430,15431,15432,15433,15434,15435,15436,15437,15438,15439,15440,15441',
       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 13:16:0.000 +60.30.0.0000',
       scan=['224,227,246,230,233,237,240,243,214,249,218,252,221'],
       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_T25t10.J131600+603000_P5312v1_2017_12_22T16_15_37.320/working/VLASS1.1.sb34397357.eb34435512.58014.72557283565_split.ms'],
       imagename='VLASS1.1.ql.T25t10.J131600+603000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00060931079441543259, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='13007,13008,13009,13010,13011,13012,13013,13014,13015,13016,13017,13018,13019,13020,13021,13022,13023,13024,13025,13026,13027,13028,13029,13030,13031,13034,13035,13036,13037,13038,13039,13040,13041,13042,13043,13044,13045,13046,13047,13048,13049,13050,13051,13052,13053,13054,13055,13056,13057,13058,13401,13402,13403,13404,13405,13406,13407,13408,13409,13410,13411,13412,13413,13414,13415,13416,13417,13418,13419,13420,13421,13422,13423,13424,13425,13426,13429,13430,13431,13432,13433,13434,13435,13436,13437,13438,13439,13440,13441,13442,13443,13444,13445,13446,13447,13448,13449,13450,13451,13452,13453,13798,13799,13800,13801,13802,13803,13804,13805,13806,13807,13808,13809,13810,13811,13812,13813,13814,13815,13816,13817,13818,13819,13820,13821,13822,13823,13826,13827,13828,13829,13830,13831,13832,13833,13834,13835,13836,13837,13838,13839,13840,13841,13842,13843,13844,13845,13846,13847,13848,13849,13850,13851,14198,14199,14200,14201,14202,14203,14204,14205,14206,14207,14208,14209,14210,14211,14212,14213,14214,14215,14216,14217,14218,14219,14220,14221,14222,14223,14226,14227,14228,14229,14230,14231,14232,14233,14234,14235,14236,14237,14238,14239,14240,14241,14242,14243,14244,14245,14246,14247,14248,14249,14250,14251,14601,14602,14603,14604,14605,14606,14607,14608,14609,14610,14611,14612,14613,14614,14615,14616,14617,14618,14619,14620,14621,14622,14623,14624,14625,14626,14629,14630,14631,14632,14633,14634,14635,14636,14637,14638,14639,14640,14641,14642,14643,14644,14645,14646,14647,14648,14649,14650,14651,14652,14653,14654,15007,15008,15009,15010,15011,15012,15013,15014,15015,15016,15017,15018,15019,15020,15021,15022,15023,15024,15025,15026,15027,15028,15029,15030,15031,15032,15035,15036,15037,15038,15039,15040,15041,15042,15043,15044,15045,15046,15047,15048,15049,15050,15051,15052,15053,15054,15055,15056,15057,15058,15059,15060,15416,15417,15418,15419,15420,15421,15422,15423,15424,15425,15426,15427,15428,15429,15430,15431,15432,15433,15434,15435,15436,15437,15438,15439,15440,15441',
       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.T25t10.J131600+603000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T25t10.J131600+603000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T25t10.J131600+603000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#