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

# hif_editimlist(parameter_file='parameter.15973.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_T12t33.J213004+073000_P15973v1_2017_11_19T22_17_07.397/working/VLASS1.1.sb34612269.eb34618429.58046.08591032407_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T12t33.J213004+073000_P15973v1_2017_11_19T22_17_07.397/rawdata/VLASS1.1.sb34612269.eb34618429.58046.08591032407.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='12808,12392,12393,11988,11989,12390,12708,12391,11986,11987,11984,11985,12408,12409,12809,12402,12403,12400,12401,12406,12314,12404,12405,11911,11910,12315,13213,12822,12312,12819,13225,12313,12310,13211,12311,13098,12701,12700,12703,12702,12705,12704,12707,12706,12709,12389,12820,12821,12799,13223,12307,12388,12805,13099,13096,12796,13094,13095,12694,12797,11902,11903,11900,11901,11906,11907,11904,11905,11908,11909,11995,12308,11478,11479,12807,13114,13204,13093,11883,11882,11881,12711,11887,11886,11885,11884,12818,11889,11888,12288,12804,12305,12806,12398,11477,12396,12014,12394,12395,12011,12010,12013,12012,11605,11604,11607,11606,11601,11600,11603,11602,11997,12301,12289,13222,11609,11608,12794,12300,12303,13226,12302,12306,11898,11899,11894,11895,11896,11897,11890,11891,11892,11893,11999,13201,12008,11998,11994,12007,12004,12005,12002,12003,12000,12001,11595,11594,11597,11596,11591,11590,11610,11592,13221,11599,11598,12803,13212,12407,11991,13224,12812,13209,11990,13220,13208,11993,13207,13206,13205,11992,13203,13202,12718,13200,11506,11507,11504,11505,11502,11503,11500,11501,11586,11587,11584,11585,11582,11583,11580,11581,13122,11996,11588,11589,12719,13120,12418,13121,12792,12285,13210,12712,13199,13198,13214,13215,12795,12713,13218,13219,12009,13197,12710,13113,12006,13111,13110,13117,13116,13115,12689,12716,13119,13118,13097,12717,12412,12714,12813,12715,12815,12814,12817,12816,12811,12810,11489,11488,11487,11486,11485,11484,11483,11482,11481,11480,12697,12696,12299,12298,12693,12692,12691,12690,12293,12292,12291,12290,12297,12296,12295,12294,13217,13108,13109,12699,13104,13105,13106,13107,13100,13101,13102,13103,11593,12798,12695,12304,13112,11498,11499,12286,12287,12802,12698,12800,12801,11490,11491,11492,11493,11494,11495,11496,11497,12309,12793,13196,12411,12410,12413,12399,12415,12414,12417,12416,13216,12397,13123',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 21:30:4.506 +07.30.0.0000',
       scan=['192,195,198,202,205,208,211,214,217,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_T12t33.J213004+073000_P15973v1_2017_11_19T22_17_07.397/working/VLASS1.1.sb34612269.eb34618429.58046.08591032407_split.ms'],
       imagename='VLASS1.1.ql.T12t33.J213004+073000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='11477,11478,11479,11480,11481,11482,11483,11484,11485,11486,11487,11488,11489,11490,11491,11492,11493,11494,11495,11496,11497,11498,11499,11500,11501,11502,11503,11504,11505,11506,11507,11580,11581,11582,11583,11584,11585,11586,11587,11588,11589,11590,11591,11592,11593,11594,11595,11596,11597,11598,11599,11600,11601,11602,11603,11604,11605,11606,11607,11608,11609,11610,11881,11882,11883,11884,11885,11886,11887,11888,11889,11890,11891,11892,11893,11894,11895,11896,11897,11898,11899,11900,11901,11902,11903,11904,11905,11906,11907,11908,11909,11910,11911,11984,11985,11986,11987,11988,11989,11990,11991,11992,11993,11994,11995,11996,11997,11998,11999,12000,12001,12002,12003,12004,12005,12006,12007,12008,12009,12010,12011,12012,12013,12014,12285,12286,12287,12288,12289,12290,12291,12292,12293,12294,12295,12296,12297,12298,12299,12300,12301,12302,12303,12304,12305,12306,12307,12308,12309,12310,12311,12312,12313,12314,12315,12388,12389,12390,12391,12392,12393,12394,12395,12396,12397,12398,12399,12400,12401,12402,12403,12404,12405,12406,12407,12408,12409,12410,12411,12412,12413,12414,12415,12416,12417,12418,12689,12690,12691,12692,12693,12694,12695,12696,12697,12698,12699,12700,12701,12702,12703,12704,12705,12706,12707,12708,12709,12710,12711,12712,12713,12714,12715,12716,12717,12718,12719,12792,12793,12794,12795,12796,12797,12798,12799,12800,12801,12802,12803,12804,12805,12806,12807,12808,12809,12810,12811,12812,12813,12814,12815,12816,12817,12818,12819,12820,12821,12822,13093,13094,13095,13096,13097,13098,13099,13100,13101,13102,13103,13104,13105,13106,13107,13108,13109,13110,13111,13112,13113,13114,13115,13116,13117,13118,13119,13120,13121,13122,13123,13196,13197,13198,13199,13200,13201,13202,13203,13204,13205,13206,13207,13208,13209,13210,13211,13212,13213,13214,13215,13216,13217,13218,13219,13220,13221,13222,13223,13224,13225,13226',
       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 21:30:4.506 +07.30.0.0000',
       scan=['192,195,198,202,205,208,211,214,217,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_T12t33.J213004+073000_P15973v1_2017_11_19T22_17_07.397/working/VLASS1.1.sb34612269.eb34618429.58046.08591032407_split.ms'],
       imagename='VLASS1.1.ql.T12t33.J213004+073000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00072092760336381616, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='11477,11478,11479,11480,11481,11482,11483,11484,11485,11486,11487,11488,11489,11490,11491,11492,11493,11494,11495,11496,11497,11498,11499,11500,11501,11502,11503,11504,11505,11506,11507,11580,11581,11582,11583,11584,11585,11586,11587,11588,11589,11590,11591,11592,11593,11594,11595,11596,11597,11598,11599,11600,11601,11602,11603,11604,11605,11606,11607,11608,11609,11610,11881,11882,11883,11884,11885,11886,11887,11888,11889,11890,11891,11892,11893,11894,11895,11896,11897,11898,11899,11900,11901,11902,11903,11904,11905,11906,11907,11908,11909,11910,11911,11984,11985,11986,11987,11988,11989,11990,11991,11992,11993,11994,11995,11996,11997,11998,11999,12000,12001,12002,12003,12004,12005,12006,12007,12008,12009,12010,12011,12012,12013,12014,12285,12286,12287,12288,12289,12290,12291,12292,12293,12294,12295,12296,12297,12298,12299,12300,12301,12302,12303,12304,12305,12306,12307,12308,12309,12310,12311,12312,12313,12314,12315,12388,12389,12390,12391,12392,12393,12394,12395,12396,12397,12398,12399,12400,12401,12402,12403,12404,12405,12406,12407,12408,12409,12410,12411,12412,12413,12414,12415,12416,12417,12418,12689,12690,12691,12692,12693,12694,12695,12696,12697,12698,12699,12700,12701,12702,12703,12704,12705,12706,12707,12708,12709,12710,12711,12712,12713,12714,12715,12716,12717,12718,12719,12792,12793,12794,12795,12796,12797,12798,12799,12800,12801,12802,12803,12804,12805,12806,12807,12808,12809,12810,12811,12812,12813,12814,12815,12816,12817,12818,12819,12820,12821,12822,13093,13094,13095,13096,13097,13098,13099,13100,13101,13102,13103,13104,13105,13106,13107,13108,13109,13110,13111,13112,13113,13114,13115,13116,13117,13118,13119,13120,13121,13122,13123,13196,13197,13198,13199,13200,13201,13202,13203,13204,13205,13206,13207,13208,13209,13210,13211,13212,13213,13214,13215,13216,13217,13218,13219,13220,13221,13222,13223,13224,13225,13226',
       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.T12t33.J213004+073000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T12t33.J213004+073000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T12t33.J213004+073000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#