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

# hif_editimlist(parameter_file='parameter.15182.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_T13t32.J205412+103000_P15182v1_2017_11_10T19_34_20.974/working/VLASS1.1.sb34607978.eb34635523.58053.05546966435_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T13t32.J205412+103000_P15182v1_2017_11_10T19_34_20.974/rawdata/VLASS1.1.sb34607978.eb34635523.58053.05546966435.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='10109,10108,11544,11545,10105,10104,10107,10106,10101,10100,10103,10102,10899,10898,10893,10892,10891,10890,10897,10896,10895,10894,10370,11546,11547,10512,10513,11548,10511,11689,11688,10516,11687,10517,10514,10515,11574,11573,11572,11948,11949,11946,11947,11944,11945,11571,11570,10093,10092,10091,10090,10097,10096,10095,10094,10099,10098,9946,9947,9944,9945,9948,9949,11168,11169,11164,11165,11166,11167,11160,11161,11162,11163,11713,11712,11711,11710,11717,11716,11715,11714,11308,11309,10349,10348,11300,11301,11302,11303,10347,10346,10345,10344,10488,10489,9689,9688,10762,10749,10748,10747,10746,10745,10744,9687,9708,9709,10887,10116,10117,9702,10115,10112,10113,9706,10111,11289,11288,11287,10888,10509,11698,11699,11696,11697,11694,11695,11692,11693,11690,11691,10905,10904,10907,10906,10901,10900,10903,10902,10909,10908,10889,10501,10500,11155,11154,11157,11156,11151,11150,11153,11152,11159,11158,11549,11973,11972,11971,11970,11974,10510,10499,10498,10497,10496,10495,10494,10493,10492,10491,10490,11560,11561,9711,9710,10088,10089,11564,10087,11565,9959,9715,9955,9954,9957,11567,9951,9950,9953,9952,11708,11709,11704,11705,11706,11707,11700,11701,11702,11703,10354,10355,10356,10357,10350,10351,10352,10353,10358,10359,10758,10759,10750,10751,10752,10753,10754,10755,10756,10757,11568,11569,9713,9712,11562,11563,9717,9716,11566,9714,11304,11305,10508,11306,11307,9964,9965,11298,11299,9960,9961,9962,9963,11292,11293,11290,11291,9968,11297,11294,11295,10916,10917,10914,10915,10912,10913,10910,10911,11146,11147,11144,11145,11148,11149,10369,10368,10365,10364,10367,10366,10361,10360,10363,10362,9966,10761,10760,10763,9967,10765,10764,10767,10766,10769,10768,10487,11551,11550,11553,11552,11555,11554,11557,11556,11559,11558,10503,10502,10505,10504,10507,10506,11296,9969,11964,11965,11966,11967,11960,11961,11962,11963,11968,11969,9700,9701,10114,9703,9704,9705,10110,9707,9958,9956,11959,11958,11955,11954,11957,11956,11951,11950,11953,11952,10374,10372,10373,9974,9973,9972,9971,9970,10371,11173,11172,11171,11170,11174,9698,11317,11316,11315,11314,11313,11312,11311,11310,10772,10773,10770,10771,10774,9699,9696,9697,9694,9695,9692,9693,9690,9691',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 20:54:12.958 +10.30.0.0000',
       scan=['192,195,164,198,167,170,173,176,179,183,186,189'],
       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_T13t32.J205412+103000_P15182v1_2017_11_10T19_34_20.974/working/VLASS1.1.sb34607978.eb34635523.58053.05546966435_split.ms'],
       imagename='VLASS1.1.ql.T13t32.J205412+103000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='9687,9688,9689,9690,9691,9692,9693,9694,9695,9696,9697,9698,9699,9700,9701,9702,9703,9704,9705,9706,9707,9708,9709,9710,9711,9712,9713,9714,9715,9716,9717,9944,9945,9946,9947,9948,9949,9950,9951,9952,9953,9954,9955,9956,9957,9958,9959,9960,9961,9962,9963,9964,9965,9966,9967,9968,9969,9970,9971,9972,9973,9974,10087,10088,10089,10090,10091,10092,10093,10094,10095,10096,10097,10098,10099,10100,10101,10102,10103,10104,10105,10106,10107,10108,10109,10110,10111,10112,10113,10114,10115,10116,10117,10344,10345,10346,10347,10348,10349,10350,10351,10352,10353,10354,10355,10356,10357,10358,10359,10360,10361,10362,10363,10364,10365,10366,10367,10368,10369,10370,10371,10372,10373,10374,10487,10488,10489,10490,10491,10492,10493,10494,10495,10496,10497,10498,10499,10500,10501,10502,10503,10504,10505,10506,10507,10508,10509,10510,10511,10512,10513,10514,10515,10516,10517,10744,10745,10746,10747,10748,10749,10750,10751,10752,10753,10754,10755,10756,10757,10758,10759,10760,10761,10762,10763,10764,10765,10766,10767,10768,10769,10770,10771,10772,10773,10774,10887,10888,10889,10890,10891,10892,10893,10894,10895,10896,10897,10898,10899,10900,10901,10902,10903,10904,10905,10906,10907,10908,10909,10910,10911,10912,10913,10914,10915,10916,10917,11144,11145,11146,11147,11148,11149,11150,11151,11152,11153,11154,11155,11156,11157,11158,11159,11160,11161,11162,11163,11164,11165,11166,11167,11168,11169,11170,11171,11172,11173,11174,11287,11288,11289,11290,11291,11292,11293,11294,11295,11296,11297,11298,11299,11300,11301,11302,11303,11304,11305,11306,11307,11308,11309,11310,11311,11312,11313,11314,11315,11316,11317,11544,11545,11546,11547,11548,11549,11550,11551,11552,11553,11554,11555,11556,11557,11558,11559,11560,11561,11562,11563,11564,11565,11566,11567,11568,11569,11570,11571,11572,11573,11574,11687,11688,11689,11690,11691,11692,11693,11694,11695,11696,11697,11698,11699,11700,11701,11702,11703,11704,11705,11706,11707,11708,11709,11710,11711,11712,11713,11714,11715,11716,11717,11944,11945,11946,11947,11948,11949,11950,11951,11952,11953,11954,11955,11956,11957,11958,11959,11960,11961,11962,11963,11964,11965,11966,11967,11968,11969,11970,11971,11972,11973,11974',
       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 20:54:12.958 +10.30.0.0000',
       scan=['192,195,164,198,167,170,173,176,179,183,186,189'],
       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_T13t32.J205412+103000_P15182v1_2017_11_10T19_34_20.974/working/VLASS1.1.sb34607978.eb34635523.58053.05546966435_split.ms'],
       imagename='VLASS1.1.ql.T13t32.J205412+103000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00047264802713543755, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='9687,9688,9689,9690,9691,9692,9693,9694,9695,9696,9697,9698,9699,9700,9701,9702,9703,9704,9705,9706,9707,9708,9709,9710,9711,9712,9713,9714,9715,9716,9717,9944,9945,9946,9947,9948,9949,9950,9951,9952,9953,9954,9955,9956,9957,9958,9959,9960,9961,9962,9963,9964,9965,9966,9967,9968,9969,9970,9971,9972,9973,9974,10087,10088,10089,10090,10091,10092,10093,10094,10095,10096,10097,10098,10099,10100,10101,10102,10103,10104,10105,10106,10107,10108,10109,10110,10111,10112,10113,10114,10115,10116,10117,10344,10345,10346,10347,10348,10349,10350,10351,10352,10353,10354,10355,10356,10357,10358,10359,10360,10361,10362,10363,10364,10365,10366,10367,10368,10369,10370,10371,10372,10373,10374,10487,10488,10489,10490,10491,10492,10493,10494,10495,10496,10497,10498,10499,10500,10501,10502,10503,10504,10505,10506,10507,10508,10509,10510,10511,10512,10513,10514,10515,10516,10517,10744,10745,10746,10747,10748,10749,10750,10751,10752,10753,10754,10755,10756,10757,10758,10759,10760,10761,10762,10763,10764,10765,10766,10767,10768,10769,10770,10771,10772,10773,10774,10887,10888,10889,10890,10891,10892,10893,10894,10895,10896,10897,10898,10899,10900,10901,10902,10903,10904,10905,10906,10907,10908,10909,10910,10911,10912,10913,10914,10915,10916,10917,11144,11145,11146,11147,11148,11149,11150,11151,11152,11153,11154,11155,11156,11157,11158,11159,11160,11161,11162,11163,11164,11165,11166,11167,11168,11169,11170,11171,11172,11173,11174,11287,11288,11289,11290,11291,11292,11293,11294,11295,11296,11297,11298,11299,11300,11301,11302,11303,11304,11305,11306,11307,11308,11309,11310,11311,11312,11313,11314,11315,11316,11317,11544,11545,11546,11547,11548,11549,11550,11551,11552,11553,11554,11555,11556,11557,11558,11559,11560,11561,11562,11563,11564,11565,11566,11567,11568,11569,11570,11571,11572,11573,11574,11687,11688,11689,11690,11691,11692,11693,11694,11695,11696,11697,11698,11699,11700,11701,11702,11703,11704,11705,11706,11707,11708,11709,11710,11711,11712,11713,11714,11715,11716,11717,11944,11945,11946,11947,11948,11949,11950,11951,11952,11953,11954,11955,11956,11957,11958,11959,11960,11961,11962,11963,11964,11965,11966,11967,11968,11969,11970,11971,11972,11973,11974',
       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.T13t32.J205412+103000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T13t32.J205412+103000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T13t32.J205412+103000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#