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

# hif_editimlist(parameter_file='parameter.8204.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_T21t20.J191816+403000_P8204v1_2017_10_05T20_11_37.206/working/VLASS1.1.sb34446513.eb34497773.58018.931944780095_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T21t20.J191816+403000_P8204v1_2017_10_05T20_11_37.206/rawdata/VLASS1.1.sb34446513.eb34497773.58018.931944780095.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='7077,6156,6926,7070,6157,6161,7071,5476,5477,5474,5475,5472,5473,5470,5471,6474,7073,6619,5478,5479,6006,6007,6004,6005,6002,6003,6000,6001,6473,6472,6471,6470,6160,6476,6008,6009,6956,7392,6957,7395,6618,7397,7396,7393,6010,7078,7399,6169,5461,5460,5463,5462,5465,5464,5467,5466,5469,5468,6953,6011,5309,6013,6012,6015,6014,6017,6016,6019,6018,5692,5693,5690,5691,5696,5697,5694,5695,6604,5698,5699,6605,6606,6607,7098,6608,6143,6142,6141,6147,6146,6145,6144,6149,6148,6927,5319,5318,5313,5312,5311,5310,5317,5316,5315,5314,6478,5689,5688,5685,5684,5687,5686,5681,5680,5683,5682,6475,6479,5326,5327,5324,5325,5322,5323,5320,5321,7092,6938,6939,5328,5480,6932,6933,7401,7422,6930,6942,6931,7419,7421,7420,6950,6936,6937,6954,6934,6477,6935,6941,5331,5330,5333,5332,5335,5334,5337,5336,5339,5338,6163,6162,6165,6164,6167,6166,7096,7097,7094,7095,5708,5709,7090,7091,5704,5705,5706,5707,5700,5701,5702,5703,7412,7089,7410,7411,7416,6170,6929,7415,7418,5329,6168,7417,6482,7093,6469,6947,6171,6613,6612,6488,6489,6617,6616,6615,6614,6158,6483,6480,6481,6486,6487,6484,6485,7081,7080,7083,7082,7085,7084,7087,7086,6946,7088,6609,6611,7099,5710,7409,7408,6159,7072,7405,7404,7407,7406,6952,7400,7403,7402,6491,6490,6493,6492,6610,6928,7413,6951,6465,6943,6949,6948,6628,6629,6626,6627,6624,6625,6622,6623,6620,6621,5458,5459,6940,6944,7079,5450,5451,5452,5453,5454,5455,5456,5457,6028,6029,6945,6020,6021,6022,6023,6024,6025,6026,6027,7414,7398,6955,6631,6630,6633,6632,6635,6634,7394,6468,6150,7069,6151,6462,6463,6464,6152,6466,6467,7074,6153,7075,6154,6030,7076,6155',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 19:18:16.542 +40.30.0.0000',
       scan=['98,102,105,108,111,82,85,89,92,95'],
       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_T21t20.J191816+403000_P8204v1_2017_10_05T20_11_37.206/working/VLASS1.1.sb34446513.eb34497773.58018.931944780095_split.ms'],
       imagename='VLASS1.1.ql.T21t20.J191816+403000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='5309,5310,5311,5312,5313,5314,5315,5316,5317,5318,5319,5320,5321,5322,5323,5324,5325,5326,5327,5328,5329,5330,5331,5332,5333,5334,5335,5336,5337,5338,5339,5450,5451,5452,5453,5454,5455,5456,5457,5458,5459,5460,5461,5462,5463,5464,5465,5466,5467,5468,5469,5470,5471,5472,5473,5474,5475,5476,5477,5478,5479,5480,5680,5681,5682,5683,5684,5685,5686,5687,5688,5689,5690,5691,5692,5693,5694,5695,5696,5697,5698,5699,5700,5701,5702,5703,5704,5705,5706,5707,5708,5709,5710,6000,6001,6002,6003,6004,6005,6006,6007,6008,6009,6010,6011,6012,6013,6014,6015,6016,6017,6018,6019,6020,6021,6022,6023,6024,6025,6026,6027,6028,6029,6030,6141,6142,6143,6144,6145,6146,6147,6148,6149,6150,6151,6152,6153,6154,6155,6156,6157,6158,6159,6160,6161,6162,6163,6164,6165,6166,6167,6168,6169,6170,6171,6462,6463,6464,6465,6466,6467,6468,6469,6470,6471,6472,6473,6474,6475,6476,6477,6478,6479,6480,6481,6482,6483,6484,6485,6486,6487,6488,6489,6490,6491,6492,6493,6604,6605,6606,6607,6608,6609,6610,6611,6612,6613,6614,6615,6616,6617,6618,6619,6620,6621,6622,6623,6624,6625,6626,6627,6628,6629,6630,6631,6632,6633,6634,6635,6926,6927,6928,6929,6930,6931,6932,6933,6934,6935,6936,6937,6938,6939,6940,6941,6942,6943,6944,6945,6946,6947,6948,6949,6950,6951,6952,6953,6954,6955,6956,6957,7069,7070,7071,7072,7073,7074,7075,7076,7077,7078,7079,7080,7081,7082,7083,7084,7085,7086,7087,7088,7089,7090,7091,7092,7093,7094,7095,7096,7097,7098,7099,7392,7393,7394,7395,7396,7397,7398,7399,7400,7401,7402,7403,7404,7405,7406,7407,7408,7409,7410,7411,7412,7413,7414,7415,7416,7417,7418,7419,7420,7421,7422',
       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 19:18:16.542 +40.30.0.0000',
       scan=['98,102,105,108,111,82,85,89,92,95'],
       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_T21t20.J191816+403000_P8204v1_2017_10_05T20_11_37.206/working/VLASS1.1.sb34446513.eb34497773.58018.931944780095_split.ms'],
       imagename='VLASS1.1.ql.T21t20.J191816+403000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00047509318556312791, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='5309,5310,5311,5312,5313,5314,5315,5316,5317,5318,5319,5320,5321,5322,5323,5324,5325,5326,5327,5328,5329,5330,5331,5332,5333,5334,5335,5336,5337,5338,5339,5450,5451,5452,5453,5454,5455,5456,5457,5458,5459,5460,5461,5462,5463,5464,5465,5466,5467,5468,5469,5470,5471,5472,5473,5474,5475,5476,5477,5478,5479,5480,5680,5681,5682,5683,5684,5685,5686,5687,5688,5689,5690,5691,5692,5693,5694,5695,5696,5697,5698,5699,5700,5701,5702,5703,5704,5705,5706,5707,5708,5709,5710,6000,6001,6002,6003,6004,6005,6006,6007,6008,6009,6010,6011,6012,6013,6014,6015,6016,6017,6018,6019,6020,6021,6022,6023,6024,6025,6026,6027,6028,6029,6030,6141,6142,6143,6144,6145,6146,6147,6148,6149,6150,6151,6152,6153,6154,6155,6156,6157,6158,6159,6160,6161,6162,6163,6164,6165,6166,6167,6168,6169,6170,6171,6462,6463,6464,6465,6466,6467,6468,6469,6470,6471,6472,6473,6474,6475,6476,6477,6478,6479,6480,6481,6482,6483,6484,6485,6486,6487,6488,6489,6490,6491,6492,6493,6604,6605,6606,6607,6608,6609,6610,6611,6612,6613,6614,6615,6616,6617,6618,6619,6620,6621,6622,6623,6624,6625,6626,6627,6628,6629,6630,6631,6632,6633,6634,6635,6926,6927,6928,6929,6930,6931,6932,6933,6934,6935,6936,6937,6938,6939,6940,6941,6942,6943,6944,6945,6946,6947,6948,6949,6950,6951,6952,6953,6954,6955,6956,6957,7069,7070,7071,7072,7073,7074,7075,7076,7077,7078,7079,7080,7081,7082,7083,7084,7085,7086,7087,7088,7089,7090,7091,7092,7093,7094,7095,7096,7097,7098,7099,7392,7393,7394,7395,7396,7397,7398,7399,7400,7401,7402,7403,7404,7405,7406,7407,7408,7409,7410,7411,7412,7413,7414,7415,7416,7417,7418,7419,7420,7421,7422',
       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.T21t20.J191816+403000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T21t20.J191816+403000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T21t20.J191816+403000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#