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

# hif_editimlist(parameter_file='parameter.14900.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_T23t17.J161553+513000_P14900v1_2018_01_05T15_14_10.803/working/VLASS1.1.sb34596438.eb34603649.58042.87761814815_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T23t17.J161553+513000_P14900v1_2018_01_05T15_14_10.803/rawdata/VLASS1.1.sb34596438.eb34603649.58042.87761814815.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='11988,11989,11986,11987,11985,344,13242,13241,13240,340,341,342,343,13162,13163,13160,13161,13164,13165,13544,13540,13541,13542,13543,12015,12014,12011,12010,12013,12012,13904,13905,13906,13907,13900,13901,13902,13903,13908,13909,39,12088,12089,339,338,335,334,337,336,331,330,333,332,13243,13247,13246,13245,13244,13645,13644,13643,13642,13641,13640,13249,13248,12758,12759,12756,12757,12754,12755,12753,12491,12490,12493,12492,12495,12494,12497,12496,12499,12498,12884,12885,12882,12883,12880,12881,12370,12371,12372,12373,12374,12375,12376,12377,12378,12379,12400,13157,13156,13155,13154,13153,13152,13151,13150,13159,13158,13258,13259,13254,13255,13256,13257,13250,13251,13252,13253,13618,13619,13615,13616,13617,12886,13531,13530,13533,13532,13535,13534,13537,13536,13539,13538,12008,12009,12006,12007,12004,12005,12002,12003,12000,12001,12877,12876,12875,12874,12873,12872,12871,12870,12879,12878,12473,12477,12476,12475,12474,12479,12478,62,13896,13897,13894,13895,13892,13893,13891,13898,13899,13265,13264,320,13266,13261,13260,13263,13262,328,329,12503,64,12502,13142,13143,13625,13624,12769,12768,13621,13620,13623,13622,12763,12762,12761,12760,12767,12766,12765,12764,12482,12483,12480,12481,12486,12487,12484,12485,12488,12489,13920,13921,12389,12388,12385,12384,12387,12386,12381,12380,12383,12382,319,318,313,317,316,315,314,13140,13141,12501,12500,13144,13145,13146,13147,13148,13149,13637,13631,13135,13137,13136,13139,13138,13917,13522,13523,13520,63,13526,13527,13524,13525,68,69,13528,13529,12860,12861,12862,12863,12864,12865,12866,12867,12868,12869,12107,12106,12105,12104,12103,12102,12101,12100,12109,12108,11999,11998,11991,11990,11993,11992,11995,11994,11997,11996,12116,12117,12118,12774,12775,12776,12777,12770,12771,12772,12773,12778,12779,13636,322,13634,13635,13632,13633,13630,323,13638,13639,321,326,327,324,325,60,59,58,61,55,54,57,56,51,50,53,52,13521,12781,12780,12783,12782,12784,65,66,67,12859,12858,13919,13918,13913,13912,13911,13910,12855,13916,12857,12856,12398,12399,12396,12397,12394,12395,12392,12393,12390,12391,12095,12094,12097,12096,12091,12090,12093,12092,12099,12098,13238,13239,13236,13237,13627,13626,13519,13518,13517,13516,13515,13514,13629,13628,12110,12111,12112,12113,12114,12115,48,49,46,47,44,45,42,43,40,41,13915,13914',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 16:15:53.409 +51.30.0.0000',
       scan=['227,230,9,12,239,242,211,236,214,233,217,220,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_T23t17.J161553+513000_P14900v1_2018_01_05T15_14_10.803/working/VLASS1.1.sb34596438.eb34603649.58042.87761814815_split.ms'],
       imagename='VLASS1.1.ql.T23t17.J161553+513000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,313,314,315,316,317,318,319,320,321,322,323,324,325,326,327,328,329,330,331,332,333,334,335,336,337,338,339,340,341,342,343,344,11985,11986,11987,11988,11989,11990,11991,11992,11993,11994,11995,11996,11997,11998,11999,12000,12001,12002,12003,12004,12005,12006,12007,12008,12009,12010,12011,12012,12013,12014,12015,12088,12089,12090,12091,12092,12093,12094,12095,12096,12097,12098,12099,12100,12101,12102,12103,12104,12105,12106,12107,12108,12109,12110,12111,12112,12113,12114,12115,12116,12117,12118,12370,12371,12372,12373,12374,12375,12376,12377,12378,12379,12380,12381,12382,12383,12384,12385,12386,12387,12388,12389,12390,12391,12392,12393,12394,12395,12396,12397,12398,12399,12400,12473,12474,12475,12476,12477,12478,12479,12480,12481,12482,12483,12484,12485,12486,12487,12488,12489,12490,12491,12492,12493,12494,12495,12496,12497,12498,12499,12500,12501,12502,12503,12753,12754,12755,12756,12757,12758,12759,12760,12761,12762,12763,12764,12765,12766,12767,12768,12769,12770,12771,12772,12773,12774,12775,12776,12777,12778,12779,12780,12781,12782,12783,12784,12855,12856,12857,12858,12859,12860,12861,12862,12863,12864,12865,12866,12867,12868,12869,12870,12871,12872,12873,12874,12875,12876,12877,12878,12879,12880,12881,12882,12883,12884,12885,12886,13135,13136,13137,13138,13139,13140,13141,13142,13143,13144,13145,13146,13147,13148,13149,13150,13151,13152,13153,13154,13155,13156,13157,13158,13159,13160,13161,13162,13163,13164,13165,13236,13237,13238,13239,13240,13241,13242,13243,13244,13245,13246,13247,13248,13249,13250,13251,13252,13253,13254,13255,13256,13257,13258,13259,13260,13261,13262,13263,13264,13265,13266,13514,13515,13516,13517,13518,13519,13520,13521,13522,13523,13524,13525,13526,13527,13528,13529,13530,13531,13532,13533,13534,13535,13536,13537,13538,13539,13540,13541,13542,13543,13544,13615,13616,13617,13618,13619,13620,13621,13622,13623,13624,13625,13626,13627,13628,13629,13630,13631,13632,13633,13634,13635,13636,13637,13638,13639,13640,13641,13642,13643,13644,13645,13891,13892,13893,13894,13895,13896,13897,13898,13899,13900,13901,13902,13903,13904,13905,13906,13907,13908,13909,13910,13911,13912,13913,13914,13915,13916,13917,13918,13919,13920,13921',
       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:53.409 +51.30.0.0000',
       scan=['227,230,9,12,239,242,211,236,214,233,217,220,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_T23t17.J161553+513000_P14900v1_2018_01_05T15_14_10.803/working/VLASS1.1.sb34596438.eb34603649.58042.87761814815_split.ms'],
       imagename='VLASS1.1.ql.T23t17.J161553+513000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00049800908351535327, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,313,314,315,316,317,318,319,320,321,322,323,324,325,326,327,328,329,330,331,332,333,334,335,336,337,338,339,340,341,342,343,344,11985,11986,11987,11988,11989,11990,11991,11992,11993,11994,11995,11996,11997,11998,11999,12000,12001,12002,12003,12004,12005,12006,12007,12008,12009,12010,12011,12012,12013,12014,12015,12088,12089,12090,12091,12092,12093,12094,12095,12096,12097,12098,12099,12100,12101,12102,12103,12104,12105,12106,12107,12108,12109,12110,12111,12112,12113,12114,12115,12116,12117,12118,12370,12371,12372,12373,12374,12375,12376,12377,12378,12379,12380,12381,12382,12383,12384,12385,12386,12387,12388,12389,12390,12391,12392,12393,12394,12395,12396,12397,12398,12399,12400,12473,12474,12475,12476,12477,12478,12479,12480,12481,12482,12483,12484,12485,12486,12487,12488,12489,12490,12491,12492,12493,12494,12495,12496,12497,12498,12499,12500,12501,12502,12503,12753,12754,12755,12756,12757,12758,12759,12760,12761,12762,12763,12764,12765,12766,12767,12768,12769,12770,12771,12772,12773,12774,12775,12776,12777,12778,12779,12780,12781,12782,12783,12784,12855,12856,12857,12858,12859,12860,12861,12862,12863,12864,12865,12866,12867,12868,12869,12870,12871,12872,12873,12874,12875,12876,12877,12878,12879,12880,12881,12882,12883,12884,12885,12886,13135,13136,13137,13138,13139,13140,13141,13142,13143,13144,13145,13146,13147,13148,13149,13150,13151,13152,13153,13154,13155,13156,13157,13158,13159,13160,13161,13162,13163,13164,13165,13236,13237,13238,13239,13240,13241,13242,13243,13244,13245,13246,13247,13248,13249,13250,13251,13252,13253,13254,13255,13256,13257,13258,13259,13260,13261,13262,13263,13264,13265,13266,13514,13515,13516,13517,13518,13519,13520,13521,13522,13523,13524,13525,13526,13527,13528,13529,13530,13531,13532,13533,13534,13535,13536,13537,13538,13539,13540,13541,13542,13543,13544,13615,13616,13617,13618,13619,13620,13621,13622,13623,13624,13625,13626,13627,13628,13629,13630,13631,13632,13633,13634,13635,13636,13637,13638,13639,13640,13641,13642,13643,13644,13645,13891,13892,13893,13894,13895,13896,13897,13898,13899,13900,13901,13902,13903,13904,13905,13906,13907,13908,13909,13910,13911,13912,13913,13914,13915,13916,13917,13918,13919,13920,13921',
       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.T23t17.J161553+513000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T23t17.J161553+513000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T23t17.J161553+513000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#