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(vis=['VLASS2.1.sb38619162.eb38631105.59096.06071193287.ms'], session=['session_1'], nocopy=True)
#
# If required, ASDMs are converted to MeasurementSets.
#

# hif_editimlist(parameter_file='parameter.list')
#
# No comment registered for hif_editimlist
#

# hif_transformimagedata(modify_weights=False, datacolumn='corrected', clear_pointing=True)
#
# No comment registered for hif_transformimagedata
#
mstransform(vis='VLASS2.1.sb38619162.eb38631105.59096.06071193287.ms',
            outputvis='VLASS2.1.sb38619162.eb38631105.59096.06071193287_split.ms',
            field='875,872,873,870,871,24,25,26,27,1836,1837,1078,1079,1076,1077,1074,1075,28,29,1070,1071,818,401,400,810,811,812,813,814,815,816,817,1838,1490,1491,1492,1493,1494,1495,1496,1497,1498,1499,1069,1068,1067,379,1857,1087,1086,1085,1084,1083,1082,1081,1080,1849,1088,393,392,391,390,397,396,395,394,399,398,874,797,1850,1851,1852,1853,1854,1855,799,798,449,448,443,442,447,446,445,444,380,381,382,383,384,385,386,387,388,389,1072,1436,1437,39,38,1432,1433,1430,1431,33,32,31,30,37,36,35,34,1846,1073,458,459,1839,450,451,452,453,454,455,456,457,1840,1421,1423,1422,1425,1424,1427,1426,1429,1428,859,1502,1503,1500,1501,1504,464,461,460,463,462,1435,1434,869,868,861,860,863,862,865,864,867,866,1843,1842,1487,1841,1486,880,1485,1847,1484,1844,1845,1848,1438,1439,1489,1488,878,879,46,877,44,45,42,43,40,41,1442,1441,1440,1856,1483,876,809,808,803,802,801,800,807,806,805,804',
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            datacolumn='corrected', chanbin=1, timebin='0s', reindex=False)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(vis=['VLASS2.1.sb38619162.eb38631105.59096.06071193287_split.ms'],
       field='24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,379,380,381,382,383,384,385,386,387,388,389,390,391,392,393,394,395,396,397,398,399,400,401,442,443,444,445,446,447,448,449,450,451,452,453,454,455,456,457,458,459,460,461,462,463,464,797,798,799,800,801,802,803,804,805,806,807,808,809,810,811,812,813,814,815,816,817,818,859,860,861,862,863,864,865,866,867,868,869,870,871,872,873,874,875,876,877,878,879,880,1067,1068,1069,1070,1071,1072,1073,1074,1075,1076,1077,1078,1079,1080,1081,1082,1083,1084,1085,1086,1087,1088,1421,1422,1423,1424,1425,1426,1427,1428,1429,1430,1431,1432,1433,1434,1435,1436,1437,1438,1439,1440,1441,1442,1483,1484,1485,1486,1487,1488,1489,1490,1491,1492,1493,1494,1495,1496,1497,1498,1499,1500,1501,1502,1503,1504,1836,1837,1838,1839,1840,1841,1842,1843,1844,1845,1846,1847,1848,1849,1850,1851,1852,1853,1854,1855,1856,1857',
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       antenna=['0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25'],
       scan=['9,12,15,18,21,25,28,31,34'], intent='OBSERVE_TARGET#UNSPECIFIED',
       datacolumn='data',
       imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter0',
       imsize=[5760, 5760], cell='1.0arcsec', phasecenter='J2000 15:51:34.069
       +24.30.0.0000', stokes='I', specmode='mfs', reffreq='3.0GHz', nchan=-1,
       outframe='LSRK', perchanweightdensity=False, gridder='mosaic',
       chanchunks=-1, wprojplanes=1, mosweight=False, conjbeams=False,
       usepointing=False, rotatepastep=360.0, pblimit=0.2, deconvolver='mtmfs',
       scales=[0], nterms=2, restoration=False, restoringbeam='common',
       pbcor=False, weighting='briggs', robust=1.0, npixels=0, niter=0,
       threshold='0.0mJy', nsigma=0, cycleniter=75, cyclefactor=1.0,
       interactive=0, fastnoise=True, savemodel='none', parallel=False)
copytree(src='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter0.residual.tt0',
         dst='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.residual.tt0')
copytree(src='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter0.psf.tt0',
         dst='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.psf.tt0')
copytree(src='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter0.weight.tt2',
         dst='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.weight.tt2')
copytree(src='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter0.model.tt0',
         dst='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.model.tt0')
copytree(src='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter0.pb.tt0',
         dst='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.pb.tt0')
copytree(src='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter0.residual.tt1',
         dst='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.residual.tt1')
copytree(src='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter0.psf.tt2',
         dst='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.psf.tt2')
copytree(src='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter0.sumwt.tt2',
         dst='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.sumwt.tt2')
copytree(src='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter0.model.tt1',
         dst='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.model.tt1')
copytree(src='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter0.sumwt.tt0',
         dst='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.sumwt.tt0')
copytree(src='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter0.weight.tt0',
         dst='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.weight.tt0')
copytree(src='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter0.sumwt.tt1',
         dst='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.sumwt.tt1')
copytree(src='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter0.weight.tt1',
         dst='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.weight.tt1')
copytree(src='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter0.psf.tt1',
         dst='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.psf.tt1')
tclean(vis=['VLASS2.1.sb38619162.eb38631105.59096.06071193287_split.ms'],
       field='24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,379,380,381,382,383,384,385,386,387,388,389,390,391,392,393,394,395,396,397,398,399,400,401,442,443,444,445,446,447,448,449,450,451,452,453,454,455,456,457,458,459,460,461,462,463,464,797,798,799,800,801,802,803,804,805,806,807,808,809,810,811,812,813,814,815,816,817,818,859,860,861,862,863,864,865,866,867,868,869,870,871,872,873,874,875,876,877,878,879,880,1067,1068,1069,1070,1071,1072,1073,1074,1075,1076,1077,1078,1079,1080,1081,1082,1083,1084,1085,1086,1087,1088,1421,1422,1423,1424,1425,1426,1427,1428,1429,1430,1431,1432,1433,1434,1435,1436,1437,1438,1439,1440,1441,1442,1483,1484,1485,1486,1487,1488,1489,1490,1491,1492,1493,1494,1495,1496,1497,1498,1499,1500,1501,1502,1503,1504,1836,1837,1838,1839,1840,1841,1842,1843,1844,1845,1846,1847,1848,1849,1850,1851,1852,1853,1854,1855,1856,1857',
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       antenna=['0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25'],
       scan=['9,12,15,18,21,25,28,31,34'], intent='OBSERVE_TARGET#UNSPECIFIED',
       datacolumn='data',
       imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1',
       imsize=[5760, 5760], cell='1.0arcsec', phasecenter='J2000 15:51:34.069
       +24.30.0.0000', stokes='I', specmode='mfs', reffreq='3.0GHz', nchan=-1,
       outframe='LSRK', perchanweightdensity=False, gridder='mosaic',
       chanchunks=-1, wprojplanes=1, mosweight=False, conjbeams=False,
       usepointing=False, rotatepastep=360.0, pblimit=0.2, deconvolver='mtmfs',
       scales=[0], nterms=2, restoration=True, restoringbeam='common',
       pbcor=False, weighting='briggs', robust=1.0, npixels=0, niter=20000,
       nsigma=4.5, cycleniter=75, cyclefactor=2.0, interactive=0,
       fastnoise=True, restart=True, savemodel='none', calcres=False,
       calcpsf=False, parallel=False)
move(src='summaryplot_1.png',
     dst='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.summaryplot.png')

# hifv_pbcor(pipelinemode="automatic")
#
# No comment registered for hifv_pbcor
#
impbcor(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.tt0',
        pbimage='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.pb.tt0',
        outfile='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.pbcor.tt0',
        mode='divide', cutoff=-1.0, stretch=False)
impbcor(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.residual.tt0',
        pbimage='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.pb.tt0',
        outfile='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.residual.pbcor.tt0',
        mode='divide', cutoff=-1.0, stretch=False)

# hif_makermsimages(pipelinemode="automatic")
#
# No comment registered for hif_makermsimages
#
imdev(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.pbcor.tt0',
      outfile='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.pbcor.tt0.rms',
      overwrite=True, stretch=False, grid=[10, 10], anchor='ref',
      xlength='60arcsec', ylength='60arcsec', interp='cubic', stattype='xmadm',
      statalg='chauvenet', zscore=-1, maxiter=-1)

# hif_makecutoutimages(pipelinemode="automatic")
#
# No comment registered for hif_makecutoutimages
#
imhead(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.tt0')
imsubimage(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.tt0',
           outfile='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.tt0.subim',
           box='1020.0,1020.0,4741.0,4741.0')
imhead(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.residual.tt0')
imsubimage(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.residual.tt0',
           outfile='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.residual.tt0.subim',
           box='1020.0,1020.0,4741.0,4741.0')
imhead(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.pbcor.tt0')
imsubimage(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.pbcor.tt0',
           outfile='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.pbcor.tt0.subim',
           box='1020.0,1020.0,4741.0,4741.0')
imhead(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imsubimage(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.pbcor.tt0.rms',
           outfile='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.pbcor.tt0.rms.subim',
           box='1020.0,1020.0,4741.0,4741.0')
imhead(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.psf.tt0')
imsubimage(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.psf.tt0',
           outfile='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.psf.tt0.subim',
           box='1020.0,1020.0,4741.0,4741.0')
imhead(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.residual.pbcor.tt0')
imsubimage(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.residual.pbcor.tt0',
           outfile='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.image.residual.pbcor.tt0.subim',
           box='1020.0,1020.0,4741.0,4741.0')
imhead(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.pb.tt0')
imsubimage(imagename='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.pb.tt0',
           outfile='VLASS2.1.ql.T17t22.J155134+243000.10.2048.v1.I.iter1.pb.tt0.subim',
           box='1020.0,1020.0,4741.0,4741.0')

# hifv_flagdata(quack=False, edgespw=False, clip=False, scan=False, autocorr=False, hm_tbuff='manual', template=False, online=False, baseband=False)
#
# No comment registered for hifv_flagdata
#
flagdata(vis='VLASS2.1.sb38619162.eb38631105.59096.06071193287_split.ms',
         mode='summary', name='before')
flagdata(vis='VLASS2.1.sb38619162.eb38631105.59096.06071193287_split.ms',
         mode='list',
         inpfile='VLASS2.1.sb38619162.eb38631105.59096.06071193287_split.flagcmds.txt',
         tbuff=[0.0, 0.0], action='apply', flagbackup=False, savepars=False)

# hifv_exportvlassdata(pipelinemode="automatic")
#
# No comment registered for hifv_exportvlassdata
#