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

# hif_editimlist(parameter_file='parameter.14242.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_T20t25.J195832+373000_P14242v1_2018_01_18T22_00_58.959/working/VLASS1.1.sb34555518.eb34684966.58067.79193422454_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T20t25.J195832+373000_P14242v1_2018_01_18T22_00_58.959/rawdata/VLASS1.1.sb34555518.eb34684966.58067.79193422454.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='22504,22505,22506,22500,22501,22502,22503,22489,22488,24047,24046,23619,24044,24043,24042,24041,24040,24794,24049,24048,24397,24411,24412,24413,24414,24415,24416,24417,24418,24792,24791,24790,24797,24796,24399,24398,24033,23275,23274,24052,23276,23271,23270,23273,23272,24429,24788,24430,24431,22849,22848,23638,22847,24045,24039,22870,22871,23266,23267,23264,23265,23262,23263,23260,23261,24422,23268,23269,24025,24024,22487,22486,24021,24020,22858,22859,22856,22857,22854,22855,22852,22853,22850,22851,24798,24400,23662,23663,23660,23661,23664,24803,24802,24801,24800,24807,24799,24805,24804,23253,23252,23251,23250,23257,23256,24410,23259,23258,24423,22496,22497,22494,22495,22492,22493,22490,22491,22863,22862,22861,22860,22867,22866,22865,22864,24421,24441,24440,24427,24426,24806,23659,23658,23657,23656,23655,23654,23653,23652,23651,23650,22874,22875,22876,22877,23248,23249,22872,22873,23244,23245,23246,23247,22878,22879,23242,23243,24793,24009,24008,24050,23648,23649,24428,24037,23640,23643,23644,23645,23646,23647,23239,23238,24035,24032,23233,23235,23234,23237,23236,24406,24018,24019,24014,24015,24016,24017,24010,24011,24012,24013,22889,22888,24808,22881,22880,22883,22882,22885,22884,22887,22886,24424,24789,24402,24425,23639,22498,23240,23635,23634,23637,22499,23631,23630,23633,23632,24028,23636,24442,24036,24438,24439,23241,24432,24433,22890,22891,24436,24437,24434,24435,24031,24027,24026,23628,23629,23626,23627,23624,23625,23622,23623,23620,23621,24407,24038,24405,24404,24403,24023,24401,24051,24034,24022,24409,24408,24787,24795',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 19:58:32.307 +37.30.0.0000',
       scan=['386,389,393,396,418,399,402,405,408,415,412,383'],
       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_T20t25.J195832+373000_P14242v1_2018_01_18T22_00_58.959/working/VLASS1.1.sb34555518.eb34684966.58067.79193422454_split.ms'],
       imagename='VLASS1.1.ql.T20t25.J195832+373000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='22486,22487,22488,22489,22490,22491,22492,22493,22494,22495,22496,22497,22498,22499,22500,22501,22502,22503,22504,22505,22506,22847,22848,22849,22850,22851,22852,22853,22854,22855,22856,22857,22858,22859,22860,22861,22862,22863,22864,22865,22866,22867,22870,22871,22872,22873,22874,22875,22876,22877,22878,22879,22880,22881,22882,22883,22884,22885,22886,22887,22888,22889,22890,22891,23233,23234,23235,23236,23237,23238,23239,23240,23241,23242,23243,23244,23245,23246,23247,23248,23249,23250,23251,23252,23253,23256,23257,23258,23259,23260,23261,23262,23263,23264,23265,23266,23267,23268,23269,23270,23271,23272,23273,23274,23275,23276,23619,23620,23621,23622,23623,23624,23625,23626,23627,23628,23629,23630,23631,23632,23633,23634,23635,23636,23637,23638,23639,23640,23643,23644,23645,23646,23647,23648,23649,23650,23651,23652,23653,23654,23655,23656,23657,23658,23659,23660,23661,23662,23663,23664,24008,24009,24010,24011,24012,24013,24014,24015,24016,24017,24018,24019,24020,24021,24022,24023,24024,24025,24026,24027,24028,24031,24032,24033,24034,24035,24036,24037,24038,24039,24040,24041,24042,24043,24044,24045,24046,24047,24048,24049,24050,24051,24052,24397,24398,24399,24400,24401,24402,24403,24404,24405,24406,24407,24408,24409,24410,24411,24412,24413,24414,24415,24416,24417,24418,24421,24422,24423,24424,24425,24426,24427,24428,24429,24430,24431,24432,24433,24434,24435,24436,24437,24438,24439,24440,24441,24442,24787,24788,24789,24790,24791,24792,24793,24794,24795,24796,24797,24798,24799,24800,24801,24802,24803,24804,24805,24806,24807,24808',
       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 19:58:32.307 +37.30.0.0000',
       scan=['386,389,393,396,418,399,402,405,408,415,412,383'],
       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_T20t25.J195832+373000_P14242v1_2018_01_18T22_00_58.959/working/VLASS1.1.sb34555518.eb34684966.58067.79193422454_split.ms'],
       imagename='VLASS1.1.ql.T20t25.J195832+373000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00045594224003907583, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='22486,22487,22488,22489,22490,22491,22492,22493,22494,22495,22496,22497,22498,22499,22500,22501,22502,22503,22504,22505,22506,22847,22848,22849,22850,22851,22852,22853,22854,22855,22856,22857,22858,22859,22860,22861,22862,22863,22864,22865,22866,22867,22870,22871,22872,22873,22874,22875,22876,22877,22878,22879,22880,22881,22882,22883,22884,22885,22886,22887,22888,22889,22890,22891,23233,23234,23235,23236,23237,23238,23239,23240,23241,23242,23243,23244,23245,23246,23247,23248,23249,23250,23251,23252,23253,23256,23257,23258,23259,23260,23261,23262,23263,23264,23265,23266,23267,23268,23269,23270,23271,23272,23273,23274,23275,23276,23619,23620,23621,23622,23623,23624,23625,23626,23627,23628,23629,23630,23631,23632,23633,23634,23635,23636,23637,23638,23639,23640,23643,23644,23645,23646,23647,23648,23649,23650,23651,23652,23653,23654,23655,23656,23657,23658,23659,23660,23661,23662,23663,23664,24008,24009,24010,24011,24012,24013,24014,24015,24016,24017,24018,24019,24020,24021,24022,24023,24024,24025,24026,24027,24028,24031,24032,24033,24034,24035,24036,24037,24038,24039,24040,24041,24042,24043,24044,24045,24046,24047,24048,24049,24050,24051,24052,24397,24398,24399,24400,24401,24402,24403,24404,24405,24406,24407,24408,24409,24410,24411,24412,24413,24414,24415,24416,24417,24418,24421,24422,24423,24424,24425,24426,24427,24428,24429,24430,24431,24432,24433,24434,24435,24436,24437,24438,24439,24440,24441,24442,24787,24788,24789,24790,24791,24792,24793,24794,24795,24796,24797,24798,24799,24800,24801,24802,24803,24804,24805,24806,24807,24808',
       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.T20t25.J195832+373000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T20t25.J195832+373000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T20t25.J195832+373000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#