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.14255.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.J192452+393000_P14255v1_2018_01_26T23_09_15.488/working/VLASS1.1.sb34555518.eb34684966.58067.79193422454_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T20t25.J192452+393000_P14255v1_2018_01_26T23_09_15.488/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='20447,20446,20445,20444,20443,20442,20441,20440,20449,20448,21567,21566,21569,21568,21083,21084,21085,21101,21100,21103,21102,21105,21086,21107,21106,21109,21108,21087,19973,19972,19971,19970,19977,19976,19975,19974,19979,19978,21208,21209,21202,21203,21200,21201,21206,21207,21204,21205,21198,21199,21192,21193,21190,21191,21196,21197,21194,21195,20821,20820,20823,20822,20825,20824,20827,20826,20829,19968,19680,19689,19688,19980,19981,20708,20709,20704,20705,20706,20707,21097,21096,21095,21094,21093,21092,21091,21090,21099,21098,21484,21485,21486,21487,21480,21481,21482,21483,21488,21489,20450,20451,20452,20453,20454,20455,20456,20457,20458,20459,21570,21571,21572,21573,21574,21575,21576,21577,21578,21579,21479,21478,21471,21470,21473,21472,21475,21474,21477,21476,20357,20356,20355,20354,20353,20352,20351,20350,20429,20428,20058,20059,20054,20055,20056,20057,20052,20053,19704,19705,19706,19707,19700,19701,19702,19703,21112,21113,21110,21111,19685,19684,19687,19686,19681,19969,19683,19682,19964,19965,19966,19967,19960,19961,19962,19963,20731,20730,20733,20732,20735,20734,21104,20832,20833,20830,20831,20836,20837,20834,20835,20083,20082,20081,20080,20065,20064,20067,20066,20061,20060,20063,20062,20069,20068,19955,19954,19957,19956,19951,19950,19953,19952,21088,21089,19959,19958,21589,21588,21585,21584,21587,21586,21581,21580,21583,21582,21468,21469,21463,21466,21467,21464,21465,20340,20341,20342,20343,20432,20345,20430,20431,20348,20349,20438,20439,20436,20437,20434,19678,19679,20435,20344,20433,20346,20347,19696,19697,19694,19695,19692,19693,19690,19691,19698,19699,20335,20334,20337,20336,20331,20330,20333,20332,20339,20338,20722,20723,20720,20721,20726,20727,20724,20725,20728,20729,20809,20808,20807,20806,20828,19677,20076,20077,20074,20075,20072,20073,20070,20071,20078,20079,21211,21210,21213,21212,21215,21214,21216,21189,21188,21187,21186,21596,21594,21595,21592,21593,21590,21591,20326,20327,20328,20329,20719,20718,20717,20716,20715,20714,20713,20712,20711,20710,20818,20819,20810,20811,20812,20813,20814,20815,20816,20817,21493,21492,21491,21490',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 19:24:52.056 +39.30.0.0000',
       scan=['355,358,361,364,367,336,339,342,345,348,351'],
       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.J192452+393000_P14255v1_2018_01_26T23_09_15.488/working/VLASS1.1.sb34555518.eb34684966.58067.79193422454_split.ms'],
       imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='19677,19678,19679,19680,19681,19682,19683,19684,19685,19686,19687,19688,19689,19690,19691,19692,19693,19694,19695,19696,19697,19698,19699,19700,19701,19702,19703,19704,19705,19706,19707,19950,19951,19952,19953,19954,19955,19956,19957,19958,19959,19960,19961,19962,19963,19964,19965,19966,19967,19968,19969,19970,19971,19972,19973,19974,19975,19976,19977,19978,19979,19980,19981,20052,20053,20054,20055,20056,20057,20058,20059,20060,20061,20062,20063,20064,20065,20066,20067,20068,20069,20070,20071,20072,20073,20074,20075,20076,20077,20078,20079,20080,20081,20082,20083,20326,20327,20328,20329,20330,20331,20332,20333,20334,20335,20336,20337,20338,20339,20340,20341,20342,20343,20344,20345,20346,20347,20348,20349,20350,20351,20352,20353,20354,20355,20356,20357,20428,20429,20430,20431,20432,20433,20434,20435,20436,20437,20438,20439,20440,20441,20442,20443,20444,20445,20446,20447,20448,20449,20450,20451,20452,20453,20454,20455,20456,20457,20458,20459,20704,20705,20706,20707,20708,20709,20710,20711,20712,20713,20714,20715,20716,20717,20718,20719,20720,20721,20722,20723,20724,20725,20726,20727,20728,20729,20730,20731,20732,20733,20734,20735,20806,20807,20808,20809,20810,20811,20812,20813,20814,20815,20816,20817,20818,20819,20820,20821,20822,20823,20824,20825,20826,20827,20828,20829,20830,20831,20832,20833,20834,20835,20836,20837,21083,21084,21085,21086,21087,21088,21089,21090,21091,21092,21093,21094,21095,21096,21097,21098,21099,21100,21101,21102,21103,21104,21105,21106,21107,21108,21109,21110,21111,21112,21113,21186,21187,21188,21189,21190,21191,21192,21193,21194,21195,21196,21197,21198,21199,21200,21201,21202,21203,21204,21205,21206,21207,21208,21209,21210,21211,21212,21213,21214,21215,21216,21463,21464,21465,21466,21467,21468,21469,21470,21471,21472,21473,21474,21475,21476,21477,21478,21479,21480,21481,21482,21483,21484,21485,21486,21487,21488,21489,21490,21491,21492,21493,21566,21567,21568,21569,21570,21571,21572,21573,21574,21575,21576,21577,21578,21579,21580,21581,21582,21583,21584,21585,21586,21587,21588,21589,21590,21591,21592,21593,21594,21595,21596',
       outframe='LSRK', gridder='mosaic', cycleniter=200, 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:24:52.056 +39.30.0.0000',
       scan=['355,358,361,364,367,336,339,342,345,348,351'],
       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.J192452+393000_P14255v1_2018_01_26T23_09_15.488/working/VLASS1.1.sb34555518.eb34684966.58067.79193422454_split.ms'],
       imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00051577667606252358, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='19677,19678,19679,19680,19681,19682,19683,19684,19685,19686,19687,19688,19689,19690,19691,19692,19693,19694,19695,19696,19697,19698,19699,19700,19701,19702,19703,19704,19705,19706,19707,19950,19951,19952,19953,19954,19955,19956,19957,19958,19959,19960,19961,19962,19963,19964,19965,19966,19967,19968,19969,19970,19971,19972,19973,19974,19975,19976,19977,19978,19979,19980,19981,20052,20053,20054,20055,20056,20057,20058,20059,20060,20061,20062,20063,20064,20065,20066,20067,20068,20069,20070,20071,20072,20073,20074,20075,20076,20077,20078,20079,20080,20081,20082,20083,20326,20327,20328,20329,20330,20331,20332,20333,20334,20335,20336,20337,20338,20339,20340,20341,20342,20343,20344,20345,20346,20347,20348,20349,20350,20351,20352,20353,20354,20355,20356,20357,20428,20429,20430,20431,20432,20433,20434,20435,20436,20437,20438,20439,20440,20441,20442,20443,20444,20445,20446,20447,20448,20449,20450,20451,20452,20453,20454,20455,20456,20457,20458,20459,20704,20705,20706,20707,20708,20709,20710,20711,20712,20713,20714,20715,20716,20717,20718,20719,20720,20721,20722,20723,20724,20725,20726,20727,20728,20729,20730,20731,20732,20733,20734,20735,20806,20807,20808,20809,20810,20811,20812,20813,20814,20815,20816,20817,20818,20819,20820,20821,20822,20823,20824,20825,20826,20827,20828,20829,20830,20831,20832,20833,20834,20835,20836,20837,21083,21084,21085,21086,21087,21088,21089,21090,21091,21092,21093,21094,21095,21096,21097,21098,21099,21100,21101,21102,21103,21104,21105,21106,21107,21108,21109,21110,21111,21112,21113,21186,21187,21188,21189,21190,21191,21192,21193,21194,21195,21196,21197,21198,21199,21200,21201,21202,21203,21204,21205,21206,21207,21208,21209,21210,21211,21212,21213,21214,21215,21216,21463,21464,21465,21466,21467,21468,21469,21470,21471,21472,21473,21474,21475,21476,21477,21478,21479,21480,21481,21482,21483,21484,21485,21486,21487,21488,21489,21490,21491,21492,21493,21566,21567,21568,21569,21570,21571,21572,21573,21574,21575,21576,21577,21578,21579,21580,21581,21582,21583,21584,21585,21586,21587,21588,21589,21590,21591,21592,21593,21594,21595,21596',
       outframe='LSRK', gridder='mosaic', cycleniter=200, 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.J192452+393000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T20t25.J192452+393000.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.J192452+393000.10.2048.v1.I.iter1.image.pbcor.tt0.rms',
      imagename='VLASS1.1.ql.T20t25.J192452+393000.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.J192452+393000.10.2048.v1.I.iter1.image.pbcor.tt0')
imhead(imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.image.pbcor.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.image.pbcor.tt0.subim',
           imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.image.pbcor.tt0')
imstat(imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imhead(imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.image.pbcor.tt0.rms.subim',
           imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imstat(imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.psf.tt0')
imhead(imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.psf.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.psf.tt0.subim',
           imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.psf.tt0')
imstat(imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.residual.tt0')
imhead(imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.residual.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.residual.tt0.subim',
           imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.residual.tt0')
imstat(imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.pb.tt0')
imhead(imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.pb.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.pb.tt0.subim',
           imagename='VLASS1.1.ql.T20t25.J192452+393000.10.2048.v1.I.iter1.pb.tt0')

# hifv_exportvlassdata(pipelinemode="automatic")
#