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

# hif_editimlist(parameter_file='parameter.12488.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_T19t21.J161212+353000_P12488v1_2017_12_13T03_32_11.174/working/VLASS1.1.sb34516741.eb34573529.58031.90910875_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T19t21.J161212+353000_P12488v1_2017_12_13T03_32_11.174/rawdata/VLASS1.1.sb34516741.eb34573529.58031.90910875.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='12921,12920,12625,12922,11829,11828,12927,12926,11825,11824,11827,11826,11821,11820,11823,11822,13307,12522,12924,12525,12213,12527,12526,12521,12520,12212,12136,13303,12529,12211,13005,12929,13007,12210,11708,11709,13003,12928,12217,11706,11707,13009,13008,12216,12932,12930,12931,12523,12215,11838,11839,11836,11837,11834,11835,11832,11833,11830,11831,13034,12535,12231,12209,13319,13318,13324,12534,12134,12532,12533,12530,12531,13315,13314,13317,13316,13311,13309,13313,13312,13016,13017,13014,13015,13012,13013,13010,13011,13320,13018,13019,12135,11809,13305,12909,12908,12907,12906,12905,12636,12903,12902,13310,13323,12108,13325,12226,12507,12506,12505,12227,13306,13030,12509,12508,13302,12224,13300,13301,13023,13022,13021,11439,13027,13026,13025,13024,12222,13029,13028,12223,12918,12919,12913,12220,12524,12910,12911,12912,12221,12914,12915,12916,12917,13020,13326,13327,12518,12519,12614,12514,12515,12230,12517,12510,12511,12512,12513,12129,12128,12228,12121,12120,12123,12122,12125,12124,12127,12126,12925,12229,13004,12923,12528,12631,12630,12633,12632,12635,12634,12637,12239,13328,13299,13006,11429,13321,12236,11428,12516,12237,11814,11815,11816,11817,11810,11811,11812,11813,11421,11420,11423,11422,11818,11819,11427,11426,12904,12233,12232,12623,12622,12132,12620,12621,12626,12627,12624,12133,12628,12629,12130,13297,12107,12106,12131,11425,12214,12109,11424,11735,11734,11438,11736,11731,11730,11733,11732,11432,11433,11430,11431,11436,11437,11434,11435,13298,12619,12618,12617,12616,12238,12225,12613,12612,12611,12610,12110,12111,12112,12113,12114,12115,12116,12117,12118,12119,11409,12235,11726,11727,11724,11725,11722,11723,11720,11721,12208,12219,12218,11728,11729,13308,12234,13033,12608,12609,13032,12606,12607,11410,11411,11412,11413,11414,11415,11416,11417,11418,11419,12504,12615,13031,11713,11712,11711,11710,11717,11716,11715,11714,11719,11718,13304,13322',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 16:12:12.465 +35.30.0.0000',
       scan=['193,196,200,203,206,209,212,215,187,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_T19t21.J161212+353000_P12488v1_2017_12_13T03_32_11.174/working/VLASS1.1.sb34516741.eb34573529.58031.90910875_split.ms'],
       imagename='VLASS1.1.ql.T19t21.J161212+353000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='11409,11410,11411,11412,11413,11414,11415,11416,11417,11418,11419,11420,11421,11422,11423,11424,11425,11426,11427,11428,11429,11430,11431,11432,11433,11434,11435,11436,11437,11438,11439,11706,11707,11708,11709,11710,11711,11712,11713,11714,11715,11716,11717,11718,11719,11720,11721,11722,11723,11724,11725,11726,11727,11728,11729,11730,11731,11732,11733,11734,11735,11736,11809,11810,11811,11812,11813,11814,11815,11816,11817,11818,11819,11820,11821,11822,11823,11824,11825,11826,11827,11828,11829,11830,11831,11832,11833,11834,11835,11836,11837,11838,11839,12106,12107,12108,12109,12110,12111,12112,12113,12114,12115,12116,12117,12118,12119,12120,12121,12122,12123,12124,12125,12126,12127,12128,12129,12130,12131,12132,12133,12134,12135,12136,12208,12209,12210,12211,12212,12213,12214,12215,12216,12217,12218,12219,12220,12221,12222,12223,12224,12225,12226,12227,12228,12229,12230,12231,12232,12233,12234,12235,12236,12237,12238,12239,12504,12505,12506,12507,12508,12509,12510,12511,12512,12513,12514,12515,12516,12517,12518,12519,12520,12521,12522,12523,12524,12525,12526,12527,12528,12529,12530,12531,12532,12533,12534,12535,12606,12607,12608,12609,12610,12611,12612,12613,12614,12615,12616,12617,12618,12619,12620,12621,12622,12623,12624,12625,12626,12627,12628,12629,12630,12631,12632,12633,12634,12635,12636,12637,12902,12903,12904,12905,12906,12907,12908,12909,12910,12911,12912,12913,12914,12915,12916,12917,12918,12919,12920,12921,12922,12923,12924,12925,12926,12927,12928,12929,12930,12931,12932,13003,13004,13005,13006,13007,13008,13009,13010,13011,13012,13013,13014,13015,13016,13017,13018,13019,13020,13021,13022,13023,13024,13025,13026,13027,13028,13029,13030,13031,13032,13033,13034,13297,13298,13299,13300,13301,13302,13303,13304,13305,13306,13307,13308,13309,13310,13311,13312,13313,13314,13315,13316,13317,13318,13319,13320,13321,13322,13323,13324,13325,13326,13327,13328',
       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 16:12:12.465 +35.30.0.0000',
       scan=['193,196,200,203,206,209,212,215,187,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_T19t21.J161212+353000_P12488v1_2017_12_13T03_32_11.174/working/VLASS1.1.sb34516741.eb34573529.58031.90910875_split.ms'],
       imagename='VLASS1.1.ql.T19t21.J161212+353000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00048139693587836513, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='11409,11410,11411,11412,11413,11414,11415,11416,11417,11418,11419,11420,11421,11422,11423,11424,11425,11426,11427,11428,11429,11430,11431,11432,11433,11434,11435,11436,11437,11438,11439,11706,11707,11708,11709,11710,11711,11712,11713,11714,11715,11716,11717,11718,11719,11720,11721,11722,11723,11724,11725,11726,11727,11728,11729,11730,11731,11732,11733,11734,11735,11736,11809,11810,11811,11812,11813,11814,11815,11816,11817,11818,11819,11820,11821,11822,11823,11824,11825,11826,11827,11828,11829,11830,11831,11832,11833,11834,11835,11836,11837,11838,11839,12106,12107,12108,12109,12110,12111,12112,12113,12114,12115,12116,12117,12118,12119,12120,12121,12122,12123,12124,12125,12126,12127,12128,12129,12130,12131,12132,12133,12134,12135,12136,12208,12209,12210,12211,12212,12213,12214,12215,12216,12217,12218,12219,12220,12221,12222,12223,12224,12225,12226,12227,12228,12229,12230,12231,12232,12233,12234,12235,12236,12237,12238,12239,12504,12505,12506,12507,12508,12509,12510,12511,12512,12513,12514,12515,12516,12517,12518,12519,12520,12521,12522,12523,12524,12525,12526,12527,12528,12529,12530,12531,12532,12533,12534,12535,12606,12607,12608,12609,12610,12611,12612,12613,12614,12615,12616,12617,12618,12619,12620,12621,12622,12623,12624,12625,12626,12627,12628,12629,12630,12631,12632,12633,12634,12635,12636,12637,12902,12903,12904,12905,12906,12907,12908,12909,12910,12911,12912,12913,12914,12915,12916,12917,12918,12919,12920,12921,12922,12923,12924,12925,12926,12927,12928,12929,12930,12931,12932,13003,13004,13005,13006,13007,13008,13009,13010,13011,13012,13013,13014,13015,13016,13017,13018,13019,13020,13021,13022,13023,13024,13025,13026,13027,13028,13029,13030,13031,13032,13033,13034,13297,13298,13299,13300,13301,13302,13303,13304,13305,13306,13307,13308,13309,13310,13311,13312,13313,13314,13315,13316,13317,13318,13319,13320,13321,13322,13323,13324,13325,13326,13327,13328',
       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.T19t21.J161212+353000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T19t21.J161212+353000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T19t21.J161212+353000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#