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

# hif_editimlist(parameter_file='parameter.7387.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_T26t04.J051707+613000_P7387v1_2018_01_01T18_46_23.386/working/VLASS1.1.sb34459038.eb34579917.58037.11761364584_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T26t04.J051707+613000_P7387v1_2018_01_01T18_46_23.386/rawdata/VLASS1.1.sb34459038.eb34579917.58037.11761364584.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='17760,17761,17762,17763,17764,17765,17766,17767,17768,17769,18554,18555,18552,18553,18550,18551,18158,18159,18152,18153,18150,18151,18156,18157,18154,18155,18952,18953,18950,18951,18956,18558,18954,18955,18958,18559,19716,19317,18556,19315,19314,19313,19312,19311,18557,19711,19319,19318,18167,18166,18165,18164,17799,17798,18160,17795,17794,17797,17796,17791,17790,17793,17792,17399,17398,17409,17408,17405,17404,17407,17406,17401,17400,17403,17402,19718,18962,18961,18960,19308,19309,18168,18170,18171,18172,18173,18174,18175,17788,17789,17786,17787,17784,17785,17782,17783,17780,17781,19336,17418,17419,17416,17417,17414,17415,17412,17413,17410,17411,19704,19705,19706,19707,19708,19709,19353,19352,19351,19350,19355,19354,18185,18184,18181,18180,18183,18182,18543,18566,19721,18957,18562,18927,18570,18571,18572,19723,19710,18959,19348,19349,19344,19345,19346,19347,19340,19341,19342,19343,17801,17800,19712,19720,18560,18176,19316,18177,19724,18178,18918,18179,18916,18917,18915,19713,18567,19338,19330,19333,18919,19335,19334,19337,19327,17759,17758,19714,17755,17757,17756,19726,19717,18929,18928,19310,18923,18922,18921,18920,18565,18926,18925,18924,19339,19328,19329,19322,19323,19320,19321,19326,18569,19324,19325,18568,18139,18163,18564,18934,18935,18936,18937,18930,18931,18932,18933,19715,18563,18529,18528,18527,18526,18525,19719,18561,19725,17776,17775,17774,17773,17772,17771,17770,18545,18544,18547,18546,18541,18540,17779,18542,18149,18148,18169,18141,18140,18143,18142,18145,18144,18147,18146,18941,18940,18943,18942,18945,18944,18947,18946,18949,18948,19722,18538,18539,18530,18531,18532,18533,18534,18535,18536,18537',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 05:17:7.679 +61.30.0.0000',
       scan=['323,326,294,297,300,304,307,310,329,313,316,319'],
       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_T26t04.J051707+613000_P7387v1_2018_01_01T18_46_23.386/working/VLASS1.1.sb34459038.eb34579917.58037.11761364584_split.ms'],
       imagename='VLASS1.1.ql.T26t04.J051707+613000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='17398,17399,17400,17401,17402,17403,17404,17405,17406,17407,17408,17409,17410,17411,17412,17413,17414,17415,17416,17417,17418,17419,17755,17756,17757,17758,17759,17760,17761,17762,17763,17764,17765,17766,17767,17768,17769,17770,17771,17772,17773,17774,17775,17776,17779,17780,17781,17782,17783,17784,17785,17786,17787,17788,17789,17790,17791,17792,17793,17794,17795,17796,17797,17798,17799,17800,17801,18139,18140,18141,18142,18143,18144,18145,18146,18147,18148,18149,18150,18151,18152,18153,18154,18155,18156,18157,18158,18159,18160,18163,18164,18165,18166,18167,18168,18169,18170,18171,18172,18173,18174,18175,18176,18177,18178,18179,18180,18181,18182,18183,18184,18185,18525,18526,18527,18528,18529,18530,18531,18532,18533,18534,18535,18536,18537,18538,18539,18540,18541,18542,18543,18544,18545,18546,18547,18550,18551,18552,18553,18554,18555,18556,18557,18558,18559,18560,18561,18562,18563,18564,18565,18566,18567,18568,18569,18570,18571,18572,18915,18916,18917,18918,18919,18920,18921,18922,18923,18924,18925,18926,18927,18928,18929,18930,18931,18932,18933,18934,18935,18936,18937,18940,18941,18942,18943,18944,18945,18946,18947,18948,18949,18950,18951,18952,18953,18954,18955,18956,18957,18958,18959,18960,18961,18962,19308,19309,19310,19311,19312,19313,19314,19315,19316,19317,19318,19319,19320,19321,19322,19323,19324,19325,19326,19327,19328,19329,19330,19333,19334,19335,19336,19337,19338,19339,19340,19341,19342,19343,19344,19345,19346,19347,19348,19349,19350,19351,19352,19353,19354,19355,19704,19705,19706,19707,19708,19709,19710,19711,19712,19713,19714,19715,19716,19717,19718,19719,19720,19721,19722,19723,19724,19725,19726',
       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 05:17:7.679 +61.30.0.0000',
       scan=['323,326,294,297,300,304,307,310,329,313,316,319'],
       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_T26t04.J051707+613000_P7387v1_2018_01_01T18_46_23.386/working/VLASS1.1.sb34459038.eb34579917.58037.11761364584_split.ms'],
       imagename='VLASS1.1.ql.T26t04.J051707+613000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00045899868975920978, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='17398,17399,17400,17401,17402,17403,17404,17405,17406,17407,17408,17409,17410,17411,17412,17413,17414,17415,17416,17417,17418,17419,17755,17756,17757,17758,17759,17760,17761,17762,17763,17764,17765,17766,17767,17768,17769,17770,17771,17772,17773,17774,17775,17776,17779,17780,17781,17782,17783,17784,17785,17786,17787,17788,17789,17790,17791,17792,17793,17794,17795,17796,17797,17798,17799,17800,17801,18139,18140,18141,18142,18143,18144,18145,18146,18147,18148,18149,18150,18151,18152,18153,18154,18155,18156,18157,18158,18159,18160,18163,18164,18165,18166,18167,18168,18169,18170,18171,18172,18173,18174,18175,18176,18177,18178,18179,18180,18181,18182,18183,18184,18185,18525,18526,18527,18528,18529,18530,18531,18532,18533,18534,18535,18536,18537,18538,18539,18540,18541,18542,18543,18544,18545,18546,18547,18550,18551,18552,18553,18554,18555,18556,18557,18558,18559,18560,18561,18562,18563,18564,18565,18566,18567,18568,18569,18570,18571,18572,18915,18916,18917,18918,18919,18920,18921,18922,18923,18924,18925,18926,18927,18928,18929,18930,18931,18932,18933,18934,18935,18936,18937,18940,18941,18942,18943,18944,18945,18946,18947,18948,18949,18950,18951,18952,18953,18954,18955,18956,18957,18958,18959,18960,18961,18962,19308,19309,19310,19311,19312,19313,19314,19315,19316,19317,19318,19319,19320,19321,19322,19323,19324,19325,19326,19327,19328,19329,19330,19333,19334,19335,19336,19337,19338,19339,19340,19341,19342,19343,19344,19345,19346,19347,19348,19349,19350,19351,19352,19353,19354,19355,19704,19705,19706,19707,19708,19709,19710,19711,19712,19713,19714,19715,19716,19717,19718,19719,19720,19721,19722,19723,19724,19725,19726',
       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.T26t04.J051707+613000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T26t04.J051707+613000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T26t04.J051707+613000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#