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

# hif_editimlist(parameter_file='parameter.5220.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_T18t22.J161519+313000_P5220v1_2018_01_12T21_06_00.709/working/VLASS1.1.sb34396276.eb34742211.58081.44143899306_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T18t22.J161519+313000_P5220v1_2018_01_12T21_06_00.709/rawdata/VLASS1.1.sb34396276.eb34742211.58081.44143899306.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='15168,15169,15417,14757,15160,15025,15162,15163,14382,15165,15166,15024,15154,15027,15172,14264,14265,14266,14267,14260,14261,14262,14263,14776,14759,14774,14775,14772,14773,14770,14771,15049,15048,14778,14779,13979,13978,15173,14381,14391,14380,13847,13846,13849,13848,15430,14386,15175,15424,15427,15426,15421,14385,15423,15422,15152,14384,15429,15428,13603,13602,13601,13600,13607,13606,13605,13604,13609,13608,15431,15174,14368,14783,15043,13852,13853,13850,13851,13856,13857,13854,13855,13858,13859,15158,15436,15414,15434,15435,15432,15433,15044,15415,15416,15420,15439,13993,13992,13991,13990,13997,13996,13995,13994,13614,13615,13999,13998,13610,13611,13612,13613,15177,15412,15176,15413,14648,14237,14236,14649,13869,13868,13867,13866,13865,13864,13863,13862,13861,13860,15040,14239,15031,15047,14238,15149,15046,15038,15039,13988,13989,15170,13984,13985,13986,13987,13980,13981,13982,13983,15042,14389,14388,13975,13974,13977,13976,13870,13871,13872,13873,13874,13875,13876,13877,14639,14638,15418,15419,14758,14633,14632,14631,14630,14637,14636,14635,14634,14657,15023,15425,15164,15034,15438,15151,15150,15153,14627,15155,15022,15157,15156,15159,14786,14785,14784,14628,14629,14781,14780,14255,14254,14257,14256,14251,14250,14253,14252,15019,14259,14258,13597,13596,13595,13594,13593,13592,13591,13590,15161,13599,13598,15032,14369,15033,15167,15030,15440,14651,14650,14653,14652,14655,14654,14387,14656,15036,15148,15045,15028,15037,14246,14247,14244,14245,14242,14243,14240,14241,15035,14365,14364,14248,14249,15442,13585,13586,13587,13588,13589,14394,14395,14390,14367,14392,14393,15437,14366,15026,14642,14643,14640,14641,14646,14647,14644,14645,14000,14001,14002,14003,14004,14005,15171,14782,14378,14379,14372,14373,14370,14371,14376,14377,14374,14375,14765,14764,14767,14766,14761,14760,14763,14762,15441,15041,14769,14768,14383,15021,15178,15020,15029,14756,14777',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 16:15:19.948 +31.30.0.0000',
       scan=['226,229,232,235,238,242,245,248,251,223'],
       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_T18t22.J161519+313000_P5220v1_2018_01_12T21_06_00.709/working/VLASS1.1.sb34396276.eb34742211.58081.44143899306_split.ms'],
       imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='13585,13586,13587,13588,13589,13590,13591,13592,13593,13594,13595,13596,13597,13598,13599,13600,13601,13602,13603,13604,13605,13606,13607,13608,13609,13610,13611,13612,13613,13614,13615,13846,13847,13848,13849,13850,13851,13852,13853,13854,13855,13856,13857,13858,13859,13860,13861,13862,13863,13864,13865,13866,13867,13868,13869,13870,13871,13872,13873,13874,13875,13876,13877,13974,13975,13976,13977,13978,13979,13980,13981,13982,13983,13984,13985,13986,13987,13988,13989,13990,13991,13992,13993,13994,13995,13996,13997,13998,13999,14000,14001,14002,14003,14004,14005,14236,14237,14238,14239,14240,14241,14242,14243,14244,14245,14246,14247,14248,14249,14250,14251,14252,14253,14254,14255,14256,14257,14258,14259,14260,14261,14262,14263,14264,14265,14266,14267,14364,14365,14366,14367,14368,14369,14370,14371,14372,14373,14374,14375,14376,14377,14378,14379,14380,14381,14382,14383,14384,14385,14386,14387,14388,14389,14390,14391,14392,14393,14394,14395,14627,14628,14629,14630,14631,14632,14633,14634,14635,14636,14637,14638,14639,14640,14641,14642,14643,14644,14645,14646,14647,14648,14649,14650,14651,14652,14653,14654,14655,14656,14657,14756,14757,14758,14759,14760,14761,14762,14763,14764,14765,14766,14767,14768,14769,14770,14771,14772,14773,14774,14775,14776,14777,14778,14779,14780,14781,14782,14783,14784,14785,14786,15019,15020,15021,15022,15023,15024,15025,15026,15027,15028,15029,15030,15031,15032,15033,15034,15035,15036,15037,15038,15039,15040,15041,15042,15043,15044,15045,15046,15047,15048,15049,15148,15149,15150,15151,15152,15153,15154,15155,15156,15157,15158,15159,15160,15161,15162,15163,15164,15165,15166,15167,15168,15169,15170,15171,15172,15173,15174,15175,15176,15177,15178,15412,15413,15414,15415,15416,15417,15418,15419,15420,15421,15422,15423,15424,15425,15426,15427,15428,15429,15430,15431,15432,15433,15434,15435,15436,15437,15438,15439,15440,15441,15442',
       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 16:15:19.948 +31.30.0.0000',
       scan=['226,229,232,235,238,242,245,248,251,223'],
       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_T18t22.J161519+313000_P5220v1_2018_01_12T21_06_00.709/working/VLASS1.1.sb34396276.eb34742211.58081.44143899306_split.ms'],
       imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00049353784652485049, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='13585,13586,13587,13588,13589,13590,13591,13592,13593,13594,13595,13596,13597,13598,13599,13600,13601,13602,13603,13604,13605,13606,13607,13608,13609,13610,13611,13612,13613,13614,13615,13846,13847,13848,13849,13850,13851,13852,13853,13854,13855,13856,13857,13858,13859,13860,13861,13862,13863,13864,13865,13866,13867,13868,13869,13870,13871,13872,13873,13874,13875,13876,13877,13974,13975,13976,13977,13978,13979,13980,13981,13982,13983,13984,13985,13986,13987,13988,13989,13990,13991,13992,13993,13994,13995,13996,13997,13998,13999,14000,14001,14002,14003,14004,14005,14236,14237,14238,14239,14240,14241,14242,14243,14244,14245,14246,14247,14248,14249,14250,14251,14252,14253,14254,14255,14256,14257,14258,14259,14260,14261,14262,14263,14264,14265,14266,14267,14364,14365,14366,14367,14368,14369,14370,14371,14372,14373,14374,14375,14376,14377,14378,14379,14380,14381,14382,14383,14384,14385,14386,14387,14388,14389,14390,14391,14392,14393,14394,14395,14627,14628,14629,14630,14631,14632,14633,14634,14635,14636,14637,14638,14639,14640,14641,14642,14643,14644,14645,14646,14647,14648,14649,14650,14651,14652,14653,14654,14655,14656,14657,14756,14757,14758,14759,14760,14761,14762,14763,14764,14765,14766,14767,14768,14769,14770,14771,14772,14773,14774,14775,14776,14777,14778,14779,14780,14781,14782,14783,14784,14785,14786,15019,15020,15021,15022,15023,15024,15025,15026,15027,15028,15029,15030,15031,15032,15033,15034,15035,15036,15037,15038,15039,15040,15041,15042,15043,15044,15045,15046,15047,15048,15049,15148,15149,15150,15151,15152,15153,15154,15155,15156,15157,15158,15159,15160,15161,15162,15163,15164,15165,15166,15167,15168,15169,15170,15171,15172,15173,15174,15175,15176,15177,15178,15412,15413,15414,15415,15416,15417,15418,15419,15420,15421,15422,15423,15424,15425,15426,15427,15428,15429,15430,15431,15432,15433,15434,15435,15436,15437,15438,15439,15440,15441,15442',
       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.T18t22.J161519+313000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T18t22.J161519+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.T18t22.J161519+313000.10.2048.v1.I.iter1.image.pbcor.tt0.rms',
      imagename='VLASS1.1.ql.T18t22.J161519+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.T18t22.J161519+313000.10.2048.v1.I.iter1.image.pbcor.tt0')
imhead(imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.image.pbcor.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.image.pbcor.tt0.subim',
           imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.image.pbcor.tt0')
imstat(imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imhead(imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.image.pbcor.tt0.rms.subim',
           imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imstat(imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.psf.tt0')
imhead(imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.psf.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.psf.tt0.subim',
           imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.psf.tt0')
imstat(imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.residual.tt0')
imhead(imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.residual.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.residual.tt0.subim',
           imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.residual.tt0')
imstat(imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.pb.tt0')
imhead(imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.pb.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.pb.tt0.subim',
           imagename='VLASS1.1.ql.T18t22.J161519+313000.10.2048.v1.I.iter1.pb.tt0')

# hifv_exportvlassdata(pipelinemode="automatic")
#