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.sb34481252.eb34557673.58027.0348240625.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='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T12t01.J000600+043000_P9063v1_2018_03_09T03_16_25.401/working/VLASS1.1.sb34481252.eb34557673.58027.0348240625_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T12t01.J000600+043000_P9063v1_2018_03_09T03_16_25.401/rawdata/VLASS1.1.sb34481252.eb34557673.58027.0348240625.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='20674,20675,20676,20677,20670,20671,20672,20673,21901,21509,21903,21902,20678,20679,21094,21508,20270,20271,20272,20273,20274,20275,20276,20277,20278,20279,21161,21160,21077,21165,21164,21076,21908,21482,21082,21483,19948,21900,21562,21089,21479,21912,21488,21910,21911,21166,21914,21489,19946,19947,19944,19945,20289,20288,19940,19941,20285,20284,20287,20286,20281,20280,20283,20282,20689,20688,21081,20681,20680,20683,20682,20685,20684,20687,20686,21484,21083,21560,21084,21073,21913,21085,21101,21100,21103,21102,21086,21091,21566,21087,21497,21088,20759,20354,20353,20352,20351,20350,20753,20752,20751,20750,20757,20756,20290,20291,21505,21504,21507,21506,21501,21500,21503,21502,20692,20693,20690,20691,20696,20697,20694,20695,21569,21499,21568,20744,20745,20746,20747,20740,20741,20742,20743,21495,20748,20749,21493,20340,20341,20342,20343,20344,20345,20346,20347,20348,20349,21097,21898,21494,21096,21892,20760,20669,21561,20668,21485,21151,21093,21092,20731,20730,20733,20732,20735,20734,20737,20736,20739,20738,21159,19919,19918,21486,21492,21909,21890,21491,20335,20334,20337,20336,20331,20330,20333,20332,21152,21885,20339,20338,21496,21487,21099,21915,21098,21564,21138,21139,21080,21889,21886,21136,21137,21079,21563,21075,21490,21894,20326,20327,20324,20325,21887,21481,20328,20329,21552,20758,21572,21565,21893,21906,21498,19937,19936,19935,19934,19933,19932,19931,19930,21907,21905,21480,19939,19938,21078,21149,21148,21145,21144,21147,21146,21141,21140,21143,21142,21549,21548,21567,21162,21897,21896,21543,21542,21545,21544,21547,21546,20755,21163,21891,20754,21904,21150,19920,19921,19922,19923,19924,19925,19926,19927,19928,19929,20267,20266,20265,20264,20263,20262,20261,21090,21156,21157,21154,21155,20667,21153,20269,20268,21888,21570,21558,21559,21899,21074,21571,21158,21553,21550,21551,21556,21557,21554,21555,21095,19942,21895,19943',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 00:06:0.879 +04.30.0.0000',
       scan=['322,325,328,332,335,338,341,344,347,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_T12t01.J000600+043000_P9063v1_2018_03_09T03_16_25.401/working/VLASS1.1.sb34481252.eb34557673.58027.0348240625_split.ms'],
       imagename='VLASS1.1.ql.T12t01.J000600+043000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='19918,19919,19920,19921,19922,19923,19924,19925,19926,19927,19928,19929,19930,19931,19932,19933,19934,19935,19936,19937,19938,19939,19940,19941,19942,19943,19944,19945,19946,19947,19948,20261,20262,20263,20264,20265,20266,20267,20268,20269,20270,20271,20272,20273,20274,20275,20276,20277,20278,20279,20280,20281,20282,20283,20284,20285,20286,20287,20288,20289,20290,20291,20324,20325,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,20667,20668,20669,20670,20671,20672,20673,20674,20675,20676,20677,20678,20679,20680,20681,20682,20683,20684,20685,20686,20687,20688,20689,20690,20691,20692,20693,20694,20695,20696,20697,20730,20731,20732,20733,20734,20735,20736,20737,20738,20739,20740,20741,20742,20743,20744,20745,20746,20747,20748,20749,20750,20751,20752,20753,20754,20755,20756,20757,20758,20759,20760,21073,21074,21075,21076,21077,21078,21079,21080,21081,21082,21083,21084,21085,21086,21087,21088,21089,21090,21091,21092,21093,21094,21095,21096,21097,21098,21099,21100,21101,21102,21103,21136,21137,21138,21139,21140,21141,21142,21143,21144,21145,21146,21147,21148,21149,21150,21151,21152,21153,21154,21155,21156,21157,21158,21159,21160,21161,21162,21163,21164,21165,21166,21479,21480,21481,21482,21483,21484,21485,21486,21487,21488,21489,21490,21491,21492,21493,21494,21495,21496,21497,21498,21499,21500,21501,21502,21503,21504,21505,21506,21507,21508,21509,21542,21543,21544,21545,21546,21547,21548,21549,21550,21551,21552,21553,21554,21555,21556,21557,21558,21559,21560,21561,21562,21563,21564,21565,21566,21567,21568,21569,21570,21571,21572,21885,21886,21887,21888,21889,21890,21891,21892,21893,21894,21895,21896,21897,21898,21899,21900,21901,21902,21903,21904,21905,21906,21907,21908,21909,21910,21911,21912,21913,21914,21915',
       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 00:06:0.879 +04.30.0.0000',
       scan=['322,325,328,332,335,338,341,344,347,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_T12t01.J000600+043000_P9063v1_2018_03_09T03_16_25.401/working/VLASS1.1.sb34481252.eb34557673.58027.0348240625_split.ms'],
       imagename='VLASS1.1.ql.T12t01.J000600+043000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00057029727907429329, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='19918,19919,19920,19921,19922,19923,19924,19925,19926,19927,19928,19929,19930,19931,19932,19933,19934,19935,19936,19937,19938,19939,19940,19941,19942,19943,19944,19945,19946,19947,19948,20261,20262,20263,20264,20265,20266,20267,20268,20269,20270,20271,20272,20273,20274,20275,20276,20277,20278,20279,20280,20281,20282,20283,20284,20285,20286,20287,20288,20289,20290,20291,20324,20325,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,20667,20668,20669,20670,20671,20672,20673,20674,20675,20676,20677,20678,20679,20680,20681,20682,20683,20684,20685,20686,20687,20688,20689,20690,20691,20692,20693,20694,20695,20696,20697,20730,20731,20732,20733,20734,20735,20736,20737,20738,20739,20740,20741,20742,20743,20744,20745,20746,20747,20748,20749,20750,20751,20752,20753,20754,20755,20756,20757,20758,20759,20760,21073,21074,21075,21076,21077,21078,21079,21080,21081,21082,21083,21084,21085,21086,21087,21088,21089,21090,21091,21092,21093,21094,21095,21096,21097,21098,21099,21100,21101,21102,21103,21136,21137,21138,21139,21140,21141,21142,21143,21144,21145,21146,21147,21148,21149,21150,21151,21152,21153,21154,21155,21156,21157,21158,21159,21160,21161,21162,21163,21164,21165,21166,21479,21480,21481,21482,21483,21484,21485,21486,21487,21488,21489,21490,21491,21492,21493,21494,21495,21496,21497,21498,21499,21500,21501,21502,21503,21504,21505,21506,21507,21508,21509,21542,21543,21544,21545,21546,21547,21548,21549,21550,21551,21552,21553,21554,21555,21556,21557,21558,21559,21560,21561,21562,21563,21564,21565,21566,21567,21568,21569,21570,21571,21572,21885,21886,21887,21888,21889,21890,21891,21892,21893,21894,21895,21896,21897,21898,21899,21900,21901,21902,21903,21904,21905,21906,21907,21908,21909,21910,21911,21912,21913,21914,21915',
       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.T12t01.J000600+043000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T12t01.J000600+043000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T12t01.J000600+043000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#