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

# hif_editimlist(parameter_file='parameter.15709.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_T25t02.J023142+583000_P15709v1_2017_11_18T00_58_43.972/working/VLASS1.1.sb34614479.eb34620102.58046.99344417824_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T25t02.J023142+583000_P15709v1_2017_11_18T00_58_43.972/rawdata/VLASS1.1.sb34614479.eb34620102.58046.99344417824.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='24047,24046,24045,24044,24043,24042,24041,24040,24049,24048,22818,22760,22761,22762,22763,22764,22765,22766,22767,22768,22769,22819,24499,24902,24903,24900,24901,24906,24907,24904,24905,24908,24909,23662,23663,23660,23661,23666,23667,23664,23665,23668,23669,24883,24887,24886,24885,24884,24889,24888,24089,24088,23248,23244,24082,24081,23247,24087,23241,23242,23243,24458,24459,24454,24455,24456,24457,24451,24452,24453,23183,23182,23228,22811,22816,22817,22814,22815,22759,22758,22753,22752,22755,23189,22757,22756,23188,22754,24524,24525,24526,24527,24520,24521,24522,24523,24528,24482,24481,24480,23200,23201,22379,23203,23619,23618,23596,23597,23598,23599,23611,23610,23617,23616,23615,23614,24050,24051,24052,24053,24465,24464,24467,24466,24461,24460,24463,24462,24469,24468,22799,22798,22797,22409,22408,22405,22404,22407,22406,22401,22400,22403,22402,24025,24024,24027,24026,24023,24022,24029,24028,23659,23658,23657,23656,23655,23654,23653,23652,23651,23650,24898,24899,24894,24895,24896,24897,24890,24891,24892,24893,24519,24518,24515,24514,24517,24516,24511,24510,24513,24512,24094,24095,24096,24097,24090,24091,24092,24093,24098,24099,22801,22800,22803,22802,22805,22804,22807,22806,23231,23230,23233,23232,23235,23234,23237,23236,24074,24069,24068,23239,23238,23202,24498,22809,22808,24497,23604,23605,23606,23607,23600,23601,23602,23603,23608,23609,22399,22398,23612,22391,22390,22393,22392,22395,22394,22397,22396,24476,24477,24474,24475,24472,24473,24470,24471,24478,24479,22782,22780,22781,24036,24037,24034,24035,24032,24033,24030,24031,24038,24039,23648,23649,23641,23642,23643,23644,23645,23646,23647,24508,24509,24506,24507,24504,24505,24502,24503,24500,24501,24085,22812,22813,23181,23180,23187,23186,23185,23184,23222,23223,23220,23221,23226,23227,23224,23225,24084,22825,24078,24079,24072,24073,24070,24071,24076,24077,22824,24075,22827,22777,22776,22775,22774,22773,22772,22771,22770,22823,22779,22778,22822,22821,22820,24911,24910,24913,24912,23671,23670,23672,24083,23245,23246,24080,23240,24086,23173,23176,23177,23174,23175,23178,23179,22388,22389,22382,22383,22380,22381,22386,22387,22384,22385,23194,22826,23624,23197,23622,23191,23192,23193,22810,23229,23219,23218,23198,23199,23626,23195,23196,23625,23190,23623,23620,23621,23613',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 02:31:42.542 +58.30.0.0000',
       scan=['386,389,392,405,395,399,408,402,373,376,380,383'],
       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_T25t02.J023142+583000_P15709v1_2017_11_18T00_58_43.972/working/VLASS1.1.sb34614479.eb34620102.58046.99344417824_split.ms'],
       imagename='VLASS1.1.ql.T25t02.J023142+583000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='22379,22380,22381,22382,22383,22384,22385,22386,22387,22388,22389,22390,22391,22392,22393,22394,22395,22396,22397,22398,22399,22400,22401,22402,22403,22404,22405,22406,22407,22408,22409,22752,22753,22754,22755,22756,22757,22758,22759,22760,22761,22762,22763,22764,22765,22766,22767,22768,22769,22770,22771,22772,22773,22774,22775,22776,22777,22778,22779,22780,22781,22782,22797,22798,22799,22800,22801,22802,22803,22804,22805,22806,22807,22808,22809,22810,22811,22812,22813,22814,22815,22816,22817,22818,22819,22820,22821,22822,22823,22824,22825,22826,22827,23173,23174,23175,23176,23177,23178,23179,23180,23181,23182,23183,23184,23185,23186,23187,23188,23189,23190,23191,23192,23193,23194,23195,23196,23197,23198,23199,23200,23201,23202,23203,23218,23219,23220,23221,23222,23223,23224,23225,23226,23227,23228,23229,23230,23231,23232,23233,23234,23235,23236,23237,23238,23239,23240,23241,23242,23243,23244,23245,23246,23247,23248,23596,23597,23598,23599,23600,23601,23602,23603,23604,23605,23606,23607,23608,23609,23610,23611,23612,23613,23614,23615,23616,23617,23618,23619,23620,23621,23622,23623,23624,23625,23626,23641,23642,23643,23644,23645,23646,23647,23648,23649,23650,23651,23652,23653,23654,23655,23656,23657,23658,23659,23660,23661,23662,23663,23664,23665,23666,23667,23668,23669,23670,23671,23672,24022,24023,24024,24025,24026,24027,24028,24029,24030,24031,24032,24033,24034,24035,24036,24037,24038,24039,24040,24041,24042,24043,24044,24045,24046,24047,24048,24049,24050,24051,24052,24053,24068,24069,24070,24071,24072,24073,24074,24075,24076,24077,24078,24079,24080,24081,24082,24083,24084,24085,24086,24087,24088,24089,24090,24091,24092,24093,24094,24095,24096,24097,24098,24099,24451,24452,24453,24454,24455,24456,24457,24458,24459,24460,24461,24462,24463,24464,24465,24466,24467,24468,24469,24470,24471,24472,24473,24474,24475,24476,24477,24478,24479,24480,24481,24482,24497,24498,24499,24500,24501,24502,24503,24504,24505,24506,24507,24508,24509,24510,24511,24512,24513,24514,24515,24516,24517,24518,24519,24520,24521,24522,24523,24524,24525,24526,24527,24528,24883,24884,24885,24886,24887,24888,24889,24890,24891,24892,24893,24894,24895,24896,24897,24898,24899,24900,24901,24902,24903,24904,24905,24906,24907,24908,24909,24910,24911,24912,24913',
       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 02:31:42.542 +58.30.0.0000',
       scan=['386,389,392,405,395,399,408,402,373,376,380,383'],
       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_T25t02.J023142+583000_P15709v1_2017_11_18T00_58_43.972/working/VLASS1.1.sb34614479.eb34620102.58046.99344417824_split.ms'],
       imagename='VLASS1.1.ql.T25t02.J023142+583000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00047077937354242746, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='22379,22380,22381,22382,22383,22384,22385,22386,22387,22388,22389,22390,22391,22392,22393,22394,22395,22396,22397,22398,22399,22400,22401,22402,22403,22404,22405,22406,22407,22408,22409,22752,22753,22754,22755,22756,22757,22758,22759,22760,22761,22762,22763,22764,22765,22766,22767,22768,22769,22770,22771,22772,22773,22774,22775,22776,22777,22778,22779,22780,22781,22782,22797,22798,22799,22800,22801,22802,22803,22804,22805,22806,22807,22808,22809,22810,22811,22812,22813,22814,22815,22816,22817,22818,22819,22820,22821,22822,22823,22824,22825,22826,22827,23173,23174,23175,23176,23177,23178,23179,23180,23181,23182,23183,23184,23185,23186,23187,23188,23189,23190,23191,23192,23193,23194,23195,23196,23197,23198,23199,23200,23201,23202,23203,23218,23219,23220,23221,23222,23223,23224,23225,23226,23227,23228,23229,23230,23231,23232,23233,23234,23235,23236,23237,23238,23239,23240,23241,23242,23243,23244,23245,23246,23247,23248,23596,23597,23598,23599,23600,23601,23602,23603,23604,23605,23606,23607,23608,23609,23610,23611,23612,23613,23614,23615,23616,23617,23618,23619,23620,23621,23622,23623,23624,23625,23626,23641,23642,23643,23644,23645,23646,23647,23648,23649,23650,23651,23652,23653,23654,23655,23656,23657,23658,23659,23660,23661,23662,23663,23664,23665,23666,23667,23668,23669,23670,23671,23672,24022,24023,24024,24025,24026,24027,24028,24029,24030,24031,24032,24033,24034,24035,24036,24037,24038,24039,24040,24041,24042,24043,24044,24045,24046,24047,24048,24049,24050,24051,24052,24053,24068,24069,24070,24071,24072,24073,24074,24075,24076,24077,24078,24079,24080,24081,24082,24083,24084,24085,24086,24087,24088,24089,24090,24091,24092,24093,24094,24095,24096,24097,24098,24099,24451,24452,24453,24454,24455,24456,24457,24458,24459,24460,24461,24462,24463,24464,24465,24466,24467,24468,24469,24470,24471,24472,24473,24474,24475,24476,24477,24478,24479,24480,24481,24482,24497,24498,24499,24500,24501,24502,24503,24504,24505,24506,24507,24508,24509,24510,24511,24512,24513,24514,24515,24516,24517,24518,24519,24520,24521,24522,24523,24524,24525,24526,24527,24528,24883,24884,24885,24886,24887,24888,24889,24890,24891,24892,24893,24894,24895,24896,24897,24898,24899,24900,24901,24902,24903,24904,24905,24906,24907,24908,24909,24910,24911,24912,24913',
       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.T25t02.J023142+583000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T25t02.J023142+583000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T25t02.J023142+583000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#