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

# hif_editimlist(parameter_file='parameter.list')
#

# hif_transformimagedata(datacolumn='corrected', modify_weights=False, clear_pointing=True)
#
mstransform(timebin='0s',
            outputvis='VLASS1.1.sb34915423.eb35019637.58158.66329216435_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T04t33.J213951-253000_P23282v1_2018_04_17T22_10_41.886/rawdata/VLASS1.1.sb34915423.eb35019637.58158.66329216435.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='10062,10063,10060,10061,10066,10067,10064,10065,10068,10069,11300,11301,11302,11289,11288,10893,10888,10891,10890,10889,11911,10665,10666,10667,11708,11709,11917,11916,11704,11705,11706,11707,10668,10669,11702,11703,10059,10058,10057,10056,10055,10054,10053,10052,10280,10281,10884,10885,10688,10689,10880,10881,10882,10883,10682,10683,10680,10681,10686,10687,10684,10685,11282,11285,11284,11902,11903,11900,11901,11906,11287,11904,11905,11908,11286,10071,10261,11298,10073,10072,10075,10074,11689,11688,10076,10879,10878,10871,10870,10873,10872,10875,10874,10877,10876,11910,11280,11089,11088,11898,11899,11083,11082,11081,11080,11087,11086,11085,11084,11498,11909,11698,11499,11696,11697,11694,11695,11692,11693,11690,11691,10868,10869,11489,10892,11292,11293,11290,11291,11296,11297,11294,11295,11094,11095,11096,11097,11090,11091,11092,11093,11098,10679,11907,10259,10258,10678,11488,11895,11896,10257,11897,11508,11509,11506,11507,11504,11505,11502,11503,11500,11501,11912,11496,11915,11497,11914,11919,10675,11918,10674,11278,11700,11699,11913,11701,11099,11279,10277,10276,10275,10274,10273,10272,10271,10270,10673,10672,10671,10670,10677,10676,10279,10278,10475,10474,10477,10476,10471,10470,10473,10472,11487,11486,11485,10479,10478,11078,11079,11710,11281,11076,11077,11074,11075,11299,10480,10481,10482,10483,10484,10260,10070,10262,10263,10264,10265,10266,10267,10268,10269,10466,10467,10464,10465,10462,10463,10460,10461,11490,11491,11492,11493,11494,11495,10468,10469,11283,10886,11712,11711,10887',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 21:39:51.648 -25.30.0.0000',
       scan=['194,165,168,172,175,178,181,184,188,191'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34915423.eb35019637.58158.66329216435_split.ms'],
       imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[5760, 5760], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='10052,10053,10054,10055,10056,10057,10058,10059,10060,10061,10062,10063,10064,10065,10066,10067,10068,10069,10070,10071,10072,10073,10074,10075,10076,10257,10258,10259,10260,10261,10262,10263,10264,10265,10266,10267,10268,10269,10270,10271,10272,10273,10274,10275,10276,10277,10278,10279,10280,10281,10460,10461,10462,10463,10464,10465,10466,10467,10468,10469,10470,10471,10472,10473,10474,10475,10476,10477,10478,10479,10480,10481,10482,10483,10484,10665,10666,10667,10668,10669,10670,10671,10672,10673,10674,10675,10676,10677,10678,10679,10680,10681,10682,10683,10684,10685,10686,10687,10688,10689,10868,10869,10870,10871,10872,10873,10874,10875,10876,10877,10878,10879,10880,10881,10882,10883,10884,10885,10886,10887,10888,10889,10890,10891,10892,10893,11074,11075,11076,11077,11078,11079,11080,11081,11082,11083,11084,11085,11086,11087,11088,11089,11090,11091,11092,11093,11094,11095,11096,11097,11098,11099,11278,11279,11280,11281,11282,11283,11284,11285,11286,11287,11288,11289,11290,11291,11292,11293,11294,11295,11296,11297,11298,11299,11300,11301,11302,11485,11486,11487,11488,11489,11490,11491,11492,11493,11494,11495,11496,11497,11498,11499,11500,11501,11502,11503,11504,11505,11506,11507,11508,11509,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,11895,11896,11897,11898,11899,11900,11901,11902,11903,11904,11905,11906,11907,11908,11909,11910,11911,11912,11913,11914,11915,11916,11917,11918,11919',
       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 21:39:51.648 -25.30.0.0000',
       scan=['194,165,168,172,175,178,181,184,188,191'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34915423.eb35019637.58158.66329216435_split.ms'],
       imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00064002174047345756, imsize=[5760, 5760],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='10052,10053,10054,10055,10056,10057,10058,10059,10060,10061,10062,10063,10064,10065,10066,10067,10068,10069,10070,10071,10072,10073,10074,10075,10076,10257,10258,10259,10260,10261,10262,10263,10264,10265,10266,10267,10268,10269,10270,10271,10272,10273,10274,10275,10276,10277,10278,10279,10280,10281,10460,10461,10462,10463,10464,10465,10466,10467,10468,10469,10470,10471,10472,10473,10474,10475,10476,10477,10478,10479,10480,10481,10482,10483,10484,10665,10666,10667,10668,10669,10670,10671,10672,10673,10674,10675,10676,10677,10678,10679,10680,10681,10682,10683,10684,10685,10686,10687,10688,10689,10868,10869,10870,10871,10872,10873,10874,10875,10876,10877,10878,10879,10880,10881,10882,10883,10884,10885,10886,10887,10888,10889,10890,10891,10892,10893,11074,11075,11076,11077,11078,11079,11080,11081,11082,11083,11084,11085,11086,11087,11088,11089,11090,11091,11092,11093,11094,11095,11096,11097,11098,11099,11278,11279,11280,11281,11282,11283,11284,11285,11286,11287,11288,11289,11290,11291,11292,11293,11294,11295,11296,11297,11298,11299,11300,11301,11302,11485,11486,11487,11488,11489,11490,11491,11492,11493,11494,11495,11496,11497,11498,11499,11500,11501,11502,11503,11504,11505,11506,11507,11508,11509,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,11895,11896,11897,11898,11899,11900,11901,11902,11903,11904,11905,11906,11907,11908,11909,11910,11911,11912,11913,11914,11915,11916,11917,11918,11919',
       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.T04t33.J213951-253000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.pb.tt0')

# hif_makermsimages(pipelinemode="automatic")
#
imdev(zscore=-1, ylength='60arcsec', stretch=False, interp='cubic',
      outfile='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.image.pbcor.tt0.rms',
      imagename='VLASS1.1.ql.T04t33.J213951-253000.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.T04t33.J213951-253000.10.2048.v1.I.iter1.image.pbcor.tt0')
imhead(imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.image.pbcor.tt0')
imsubimage(box='1020,1020,4741,4741',
           outfile='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.image.pbcor.tt0.subim',
           imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.image.pbcor.tt0')
imstat(imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imhead(imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imsubimage(box='1020,1020,4741,4741',
           outfile='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.image.pbcor.tt0.rms.subim',
           imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imstat(imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.psf.tt0')
imhead(imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.psf.tt0')
imsubimage(box='1020,1020,4741,4741',
           outfile='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.psf.tt0.subim',
           imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.psf.tt0')
imstat(imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.residual.tt0')
imhead(imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.residual.tt0')
imsubimage(box='1020,1020,4741,4741',
           outfile='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.residual.tt0.subim',
           imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.residual.tt0')
imstat(imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.pb.tt0')
imhead(imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.pb.tt0')
imsubimage(box='1020,1020,4741,4741',
           outfile='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.pb.tt0.subim',
           imagename='VLASS1.1.ql.T04t33.J213951-253000.10.2048.v1.I.iter1.pb.tt0')

# hifv_exportvlassdata(pipelinemode="automatic")
#