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

# hif_editimlist(parameter_file='parameter.11892.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_T18t21.J153019+313000_P11892v1_2017_11_19T19_29_20.588/working/VLASS1.1.sb34516171.eb34560849.58028.92610690973_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T18t21.J153019+313000_P11892v1_2017_11_19T19_29_20.588/rawdata/VLASS1.1.sb34516171.eb34560849.58028.92610690973.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='12828,13460,12707,13467,12706,13477,13466,12833,12824,12831,12830,12837,12699,12835,12825,12698,12829,13200,12679,12678,11661,11660,13085,13084,13087,13086,13081,13080,12822,13082,11911,11910,11913,11912,11915,11914,11917,11916,11919,11918,13483,12823,13485,13484,13487,13469,13218,13211,13196,13212,13465,12437,12436,12435,12434,12433,12432,12431,12430,12709,12708,12826,12827,12820,12821,12439,12438,12315,13098,12832,12701,13097,13094,13095,13092,12694,13090,13091,11902,11903,11900,11901,11906,11907,11904,11905,11908,11909,12309,12308,12836,12700,13069,13096,13204,12834,13093,12420,12421,12422,12423,12424,12425,12426,12427,12428,12429,13468,12305,13068,12304,12307,12316,12317,12314,12306,12312,12313,12310,12311,13221,12301,13223,13222,13225,13224,13227,12300,13207,12702,12303,12302,12444,13073,12445,11898,11899,11895,11896,11897,13099,13088,12705,12806,13201,12807,13070,13476,13089,13488,12287,11920,11921,11922,11923,11924,11925,12704,13471,13472,13226,13215,13470,13481,13461,13083,13473,13220,13208,13478,12446,13206,12039,12038,12442,12443,12440,12441,12033,12032,12031,12030,12037,12036,12035,12034,13475,13479,13213,12680,13459,12681,13072,12024,12025,12026,12027,13480,13074,13458,12028,12029,13075,13210,12815,13199,13198,13214,13076,13216,12814,12810,13219,13077,13197,12817,13078,12816,11638,11639,13079,11634,11635,11636,11637,11631,11632,11633,13071,13205,12051,12050,12053,12052,12811,12054,12813,12812,12819,12818,12697,12696,12299,12298,12693,12692,12691,12690,12293,12292,12291,12290,12297,12296,12295,12294,11649,11648,13217,13474,13202,11641,11640,11643,11642,11645,11644,11647,11646,13482,13464,13203,12695,12042,12043,12040,12041,12046,12047,12044,12045,12048,12049,12288,12289,12808,12809,12419,12418,12682,12683,12684,12685,12686,12687,12688,12689,12703,12417,12416,13209,13463,11658,11659,13486,13462,11652,11653,11650,11651,11656,11657,11654,11655',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 15:30:19.948 +31.30.0.0000',
       scan=['193,196,199,203,206,209,212,215,218,190'],
       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_T18t21.J153019+313000_P11892v1_2017_11_19T19_29_20.588/working/VLASS1.1.sb34516171.eb34560849.58028.92610690973_split.ms'],
       imagename='VLASS1.1.ql.T18t21.J153019+313000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='11631,11632,11633,11634,11635,11636,11637,11638,11639,11640,11641,11642,11643,11644,11645,11646,11647,11648,11649,11650,11651,11652,11653,11654,11655,11656,11657,11658,11659,11660,11661,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,11920,11921,11922,11923,11924,11925,12024,12025,12026,12027,12028,12029,12030,12031,12032,12033,12034,12035,12036,12037,12038,12039,12040,12041,12042,12043,12044,12045,12046,12047,12048,12049,12050,12051,12052,12053,12054,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,12316,12317,12416,12417,12418,12419,12420,12421,12422,12423,12424,12425,12426,12427,12428,12429,12430,12431,12432,12433,12434,12435,12436,12437,12438,12439,12440,12441,12442,12443,12444,12445,12446,12678,12679,12680,12681,12682,12683,12684,12685,12686,12687,12688,12689,12690,12691,12692,12693,12694,12695,12696,12697,12698,12699,12700,12701,12702,12703,12704,12705,12706,12707,12708,12709,12806,12807,12808,12809,12810,12811,12812,12813,12814,12815,12816,12817,12818,12819,12820,12821,12822,12823,12824,12825,12826,12827,12828,12829,12830,12831,12832,12833,12834,12835,12836,12837,13068,13069,13070,13071,13072,13073,13074,13075,13076,13077,13078,13079,13080,13081,13082,13083,13084,13085,13086,13087,13088,13089,13090,13091,13092,13093,13094,13095,13096,13097,13098,13099,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,13227,13458,13459,13460,13461,13462,13463,13464,13465,13466,13467,13468,13469,13470,13471,13472,13473,13474,13475,13476,13477,13478,13479,13480,13481,13482,13483,13484,13485,13486,13487,13488',
       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 15:30:19.948 +31.30.0.0000',
       scan=['193,196,199,203,206,209,212,215,218,190'],
       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_T18t21.J153019+313000_P11892v1_2017_11_19T19_29_20.588/working/VLASS1.1.sb34516171.eb34560849.58028.92610690973_split.ms'],
       imagename='VLASS1.1.ql.T18t21.J153019+313000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00063378752936332598, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='11631,11632,11633,11634,11635,11636,11637,11638,11639,11640,11641,11642,11643,11644,11645,11646,11647,11648,11649,11650,11651,11652,11653,11654,11655,11656,11657,11658,11659,11660,11661,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,11920,11921,11922,11923,11924,11925,12024,12025,12026,12027,12028,12029,12030,12031,12032,12033,12034,12035,12036,12037,12038,12039,12040,12041,12042,12043,12044,12045,12046,12047,12048,12049,12050,12051,12052,12053,12054,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,12316,12317,12416,12417,12418,12419,12420,12421,12422,12423,12424,12425,12426,12427,12428,12429,12430,12431,12432,12433,12434,12435,12436,12437,12438,12439,12440,12441,12442,12443,12444,12445,12446,12678,12679,12680,12681,12682,12683,12684,12685,12686,12687,12688,12689,12690,12691,12692,12693,12694,12695,12696,12697,12698,12699,12700,12701,12702,12703,12704,12705,12706,12707,12708,12709,12806,12807,12808,12809,12810,12811,12812,12813,12814,12815,12816,12817,12818,12819,12820,12821,12822,12823,12824,12825,12826,12827,12828,12829,12830,12831,12832,12833,12834,12835,12836,12837,13068,13069,13070,13071,13072,13073,13074,13075,13076,13077,13078,13079,13080,13081,13082,13083,13084,13085,13086,13087,13088,13089,13090,13091,13092,13093,13094,13095,13096,13097,13098,13099,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,13227,13458,13459,13460,13461,13462,13463,13464,13465,13466,13467,13468,13469,13470,13471,13472,13473,13474,13475,13476,13477,13478,13479,13480,13481,13482,13483,13484,13485,13486,13487,13488',
       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.T18t21.J153019+313000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T18t21.J153019+313000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T18t21.J153019+313000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#