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.sb34475744.eb34497997.58020.06194271991.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_T11t36.J233000+033000_P9868v1_2018_03_14T21_07_54.739/working/VLASS1.1.sb34475744.eb34497997.58020.06194271991_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T11t36.J233000+033000_P9868v1_2018_03_14T21_07_54.739/rawdata/VLASS1.1.sb34475744.eb34497997.58020.06194271991.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='14787,15168,14883,14891,15169,14884,15165,15166,14885,13368,14079,14078,13360,13361,13362,13363,13364,13365,13366,13367,14479,14478,14774,14775,14772,14773,14770,14771,14471,14470,14473,14472,14475,14474,14477,14476,14383,14382,14892,14381,14353,14354,14355,13847,14357,13849,13848,14077,14068,14069,14066,14067,14064,14065,14062,14063,14060,14061,15195,15185,14468,14469,15191,15190,15193,14788,14462,14463,14460,14461,14466,14467,14464,14465,14076,13764,14887,14368,13765,15167,15186,13766,13852,13853,13850,13767,13856,13857,13854,13855,13760,13858,13859,14864,13761,14457,14075,14059,14058,13762,14074,14053,14052,14051,13763,14057,14056,14055,14054,13755,13754,13757,13756,13751,13750,13753,13752,15179,14071,15172,15178,13759,13758,15188,14070,14866,14073,15184,14072,14362,14776,13869,13868,13867,13866,13865,13864,13863,13862,13861,13860,14867,15171,14869,13746,13747,13339,13338,14759,15170,13748,13749,14380,14379,13038,13039,14364,13035,13036,13037,14356,14862,15189,15194,13870,13871,13872,13873,13874,13875,13876,13877,15182,14863,13851,14778,15180,14779,15192,13773,13772,13771,13770,13449,13448,13774,13445,13444,13447,13446,13441,13443,13442,13049,13048,13041,13040,13043,13042,13045,13044,13047,13046,14789,13745,14456,14786,14785,14784,14783,13744,14781,14780,13351,13350,13353,13352,13355,13354,13357,13356,13359,13358,14369,13768,13769,13458,13459,13456,13457,13454,13455,13452,13453,13450,13451,13058,13059,13052,13053,13050,13051,13056,13057,13054,13055,14875,14874,14877,14876,14871,14870,14873,14872,14879,14878,13342,13343,13340,13341,13346,13347,13344,13345,14361,14360,13348,13349,14365,14363,14367,14366,14888,14889,13065,13064,13063,13062,13061,13060,13463,13462,13461,13460,13467,13466,13465,13464,14358,14459,14359,14886,15183,14484,14458,14050,14485,14486,14865,14480,14481,14482,14483,15177,15176,15175,15174,15173,13469,14868,14782,14080,14378,13468,13470,13471,14372,14373,14370,14371,14376,14377,14374,14375,14765,14764,14767,14766,14761,14760,14763,14762,15181,14769,14768,14880,15187,14881,14882,14890,14777',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 23:30:0.823 +03.30.0.0000',
       scan=['225,228,231,234,238,241,244,247,250,253,222'],
       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_T11t36.J233000+033000_P9868v1_2018_03_14T21_07_54.739/working/VLASS1.1.sb34475744.eb34497997.58020.06194271991_split.ms'],
       imagename='VLASS1.1.ql.T11t36.J233000+033000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='13035,13036,13037,13038,13039,13040,13041,13042,13043,13044,13045,13046,13047,13048,13049,13050,13051,13052,13053,13054,13055,13056,13057,13058,13059,13060,13061,13062,13063,13064,13065,13338,13339,13340,13341,13342,13343,13344,13345,13346,13347,13348,13349,13350,13351,13352,13353,13354,13355,13356,13357,13358,13359,13360,13361,13362,13363,13364,13365,13366,13367,13368,13441,13442,13443,13444,13445,13446,13447,13448,13449,13450,13451,13452,13453,13454,13455,13456,13457,13458,13459,13460,13461,13462,13463,13464,13465,13466,13467,13468,13469,13470,13471,13744,13745,13746,13747,13748,13749,13750,13751,13752,13753,13754,13755,13756,13757,13758,13759,13760,13761,13762,13763,13764,13765,13766,13767,13768,13769,13770,13771,13772,13773,13774,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,14050,14051,14052,14053,14054,14055,14056,14057,14058,14059,14060,14061,14062,14063,14064,14065,14066,14067,14068,14069,14070,14071,14072,14073,14074,14075,14076,14077,14078,14079,14080,14353,14354,14355,14356,14357,14358,14359,14360,14361,14362,14363,14364,14365,14366,14367,14368,14369,14370,14371,14372,14373,14374,14375,14376,14377,14378,14379,14380,14381,14382,14383,14456,14457,14458,14459,14460,14461,14462,14463,14464,14465,14466,14467,14468,14469,14470,14471,14472,14473,14474,14475,14476,14477,14478,14479,14480,14481,14482,14483,14484,14485,14486,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,14787,14788,14789,14862,14863,14864,14865,14866,14867,14868,14869,14870,14871,14872,14873,14874,14875,14876,14877,14878,14879,14880,14881,14882,14883,14884,14885,14886,14887,14888,14889,14890,14891,14892,15165,15166,15167,15168,15169,15170,15171,15172,15173,15174,15175,15176,15177,15178,15179,15180,15181,15182,15183,15184,15185,15186,15187,15188,15189,15190,15191,15192,15193,15194,15195',
       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 23:30:0.823 +03.30.0.0000',
       scan=['225,228,231,234,238,241,244,247,250,253,222'],
       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_T11t36.J233000+033000_P9868v1_2018_03_14T21_07_54.739/working/VLASS1.1.sb34475744.eb34497997.58020.06194271991_split.ms'],
       imagename='VLASS1.1.ql.T11t36.J233000+033000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00055356193333227566, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='13035,13036,13037,13038,13039,13040,13041,13042,13043,13044,13045,13046,13047,13048,13049,13050,13051,13052,13053,13054,13055,13056,13057,13058,13059,13060,13061,13062,13063,13064,13065,13338,13339,13340,13341,13342,13343,13344,13345,13346,13347,13348,13349,13350,13351,13352,13353,13354,13355,13356,13357,13358,13359,13360,13361,13362,13363,13364,13365,13366,13367,13368,13441,13442,13443,13444,13445,13446,13447,13448,13449,13450,13451,13452,13453,13454,13455,13456,13457,13458,13459,13460,13461,13462,13463,13464,13465,13466,13467,13468,13469,13470,13471,13744,13745,13746,13747,13748,13749,13750,13751,13752,13753,13754,13755,13756,13757,13758,13759,13760,13761,13762,13763,13764,13765,13766,13767,13768,13769,13770,13771,13772,13773,13774,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,14050,14051,14052,14053,14054,14055,14056,14057,14058,14059,14060,14061,14062,14063,14064,14065,14066,14067,14068,14069,14070,14071,14072,14073,14074,14075,14076,14077,14078,14079,14080,14353,14354,14355,14356,14357,14358,14359,14360,14361,14362,14363,14364,14365,14366,14367,14368,14369,14370,14371,14372,14373,14374,14375,14376,14377,14378,14379,14380,14381,14382,14383,14456,14457,14458,14459,14460,14461,14462,14463,14464,14465,14466,14467,14468,14469,14470,14471,14472,14473,14474,14475,14476,14477,14478,14479,14480,14481,14482,14483,14484,14485,14486,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,14787,14788,14789,14862,14863,14864,14865,14866,14867,14868,14869,14870,14871,14872,14873,14874,14875,14876,14877,14878,14879,14880,14881,14882,14883,14884,14885,14886,14887,14888,14889,14890,14891,14892,15165,15166,15167,15168,15169,15170,15171,15172,15173,15174,15175,15176,15177,15178,15179,15180,15181,15182,15183,15184,15185,15186,15187,15188,15189,15190,15191,15192,15193,15194,15195',
       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.T11t36.J233000+033000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T11t36.J233000+033000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T11t36.J233000+033000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#