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

# hif_editimlist(parameter_file='parameter.4729.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_T23t19.J184640+503000_P4729v1_2017_09_15T20_53_07.350/working/VLASS1.1.sb34361741.eb34367683.58005.085977442126_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T23t19.J184640+503000_P4729v1_2017_09_15T20_53_07.350/rawdata/VLASS1.1.sb34361741.eb34367683.58005.085977442126.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='10518,10519,10109,10108,10105,10104,10107,10106,10101,10100,10103,10102,10899,10898,10893,10892,10891,10890,10897,10896,10895,10894,10581,10580,10583,10582,10585,10584,10587,10586,10589,10588,10206,10207,10204,10205,10202,10203,10200,10201,10208,10209,10196,10197,10194,10195,10192,10193,10190,10191,10198,10199,10512,10513,10510,10511,10516,10517,10514,10515,8666,9436,8580,9435,9410,9411,9412,9413,9414,9415,9416,9417,9418,9419,9731,9730,9733,9732,9735,9734,9737,9736,9739,9738,9328,9329,9322,9323,9321,9326,9327,9324,9325,8941,8940,8943,8942,8945,8944,8947,8946,8949,8948,9029,9028,9025,9024,9027,9026,9023,9022,10099,10098,8935,8936,8937,8938,8939,9044,8648,8649,8646,8647,8644,8645,8642,8643,8640,8641,10521,10520,9341,9709,10886,10887,10118,10119,10116,10117,10114,10115,10112,10113,10110,10111,10592,10593,10590,10591,10596,10597,10594,10595,10598,10599,10905,10904,10907,10906,10901,10900,10903,10902,10909,10908,9043,9042,10503,9040,9047,9046,9045,9409,9408,9335,9334,9337,9336,9331,9330,9333,9332,9339,9338,8574,8575,8576,8577,8570,8571,8572,8573,8578,8579,10129,9718,10499,10498,10497,10496,10495,10494,10493,10492,10491,10490,9713,10579,10122,9050,9051,9052,10121,10120,10127,10126,10125,10124,8659,8658,8651,8650,8653,8652,8655,8654,8657,8656,9825,9824,9827,9826,9821,9820,9823,9822,9719,10128,10123,9712,9711,9710,9717,9716,9715,9714,8963,8962,8961,8960,8965,8964,10914,10912,10913,10910,10911,8660,8661,8662,8663,8664,8665,9438,8667,8581,9437,9434,8582,9432,9433,9430,9431,9340,8668,9342,9343,9344,9345,9346,9347,9348,9349,8567,8566,8565,8564,8563,8562,8561,8560,8569,8568,9818,9819,9814,9815,9816,9817,9810,9811,9812,9813,10509,10508,9049,9048,10501,10500,9041,10502,10505,10504,10507,10506,10884,10885,10215,10214,10217,10216,10211,10210,10213,10212,10187,10186,10888,10189,10188,10889,8639,8638,10602,10603,10600,10601,10606,10607,10604,10605,10608,10609,9726,9727,9724,9725,9722,9723,9720,9721,9728,9729,8956,8957,8954,8955,8952,8953,8950,8951,8958,8959,9797,9796,9799,9798,9032,9033,9030,9031,9036,9037,9034,9035,9038,9039,9429,9428,9421,9420,9423,9422,9425,9424,9427,9426,9351,9350,8552,8553,8551,8556,8557,8554,8555,8558,8559,9809,9808,9807,9806,9805,9804,9803,9802,9801,9800',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 18:46:40.303 +50.30.0.0000',
       scan=['162,165,169,172,194,175,178,188,181,184,191,156,159'],
       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_T23t19.J184640+503000_P4729v1_2017_09_15T20_53_07.350/working/VLASS1.1.sb34361741.eb34367683.58005.085977442126_split.ms'],
       imagename='VLASS1.1.ql.T23t19.J184640+503000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='8551,8552,8553,8554,8555,8556,8557,8558,8559,8560,8561,8562,8563,8564,8565,8566,8567,8568,8569,8570,8571,8572,8573,8574,8575,8576,8577,8578,8579,8580,8581,8582,8638,8639,8640,8641,8642,8643,8644,8645,8646,8647,8648,8649,8650,8651,8652,8653,8654,8655,8656,8657,8658,8659,8660,8661,8662,8663,8664,8665,8666,8667,8668,8935,8936,8937,8938,8939,8940,8941,8942,8943,8944,8945,8946,8947,8948,8949,8950,8951,8952,8953,8954,8955,8956,8957,8958,8959,8960,8961,8962,8963,8964,8965,9022,9023,9024,9025,9026,9027,9028,9029,9030,9031,9032,9033,9034,9035,9036,9037,9038,9039,9040,9041,9042,9043,9044,9045,9046,9047,9048,9049,9050,9051,9052,9321,9322,9323,9324,9325,9326,9327,9328,9329,9330,9331,9332,9333,9334,9335,9336,9337,9338,9339,9340,9341,9342,9343,9344,9345,9346,9347,9348,9349,9350,9351,9408,9409,9410,9411,9412,9413,9414,9415,9416,9417,9418,9419,9420,9421,9422,9423,9424,9425,9426,9427,9428,9429,9430,9431,9432,9433,9434,9435,9436,9437,9438,9709,9710,9711,9712,9713,9714,9715,9716,9717,9718,9719,9720,9721,9722,9723,9724,9725,9726,9727,9728,9729,9730,9731,9732,9733,9734,9735,9736,9737,9738,9739,9796,9797,9798,9799,9800,9801,9802,9803,9804,9805,9806,9807,9808,9809,9810,9811,9812,9813,9814,9815,9816,9817,9818,9819,9820,9821,9822,9823,9824,9825,9826,9827,10098,10099,10100,10101,10102,10103,10104,10105,10106,10107,10108,10109,10110,10111,10112,10113,10114,10115,10116,10117,10118,10119,10120,10121,10122,10123,10124,10125,10126,10127,10128,10129,10186,10187,10188,10189,10190,10191,10192,10193,10194,10195,10196,10197,10198,10199,10200,10201,10202,10203,10204,10205,10206,10207,10208,10209,10210,10211,10212,10213,10214,10215,10216,10217,10490,10491,10492,10493,10494,10495,10496,10497,10498,10499,10500,10501,10502,10503,10504,10505,10506,10507,10508,10509,10510,10511,10512,10513,10514,10515,10516,10517,10518,10519,10520,10521,10579,10580,10581,10582,10583,10584,10585,10586,10587,10588,10589,10590,10591,10592,10593,10594,10595,10596,10597,10598,10599,10600,10601,10602,10603,10604,10605,10606,10607,10608,10609,10884,10885,10886,10887,10888,10889,10890,10891,10892,10893,10894,10895,10896,10897,10898,10899,10900,10901,10902,10903,10904,10905,10906,10907,10908,10909,10910,10911,10912,10913,10914',
       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 18:46:40.303 +50.30.0.0000',
       scan=['162,165,169,172,194,175,178,188,181,184,191,156,159'],
       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_T23t19.J184640+503000_P4729v1_2017_09_15T20_53_07.350/working/VLASS1.1.sb34361741.eb34367683.58005.085977442126_split.ms'],
       imagename='VLASS1.1.ql.T23t19.J184640+503000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00053653568283698399, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='8551,8552,8553,8554,8555,8556,8557,8558,8559,8560,8561,8562,8563,8564,8565,8566,8567,8568,8569,8570,8571,8572,8573,8574,8575,8576,8577,8578,8579,8580,8581,8582,8638,8639,8640,8641,8642,8643,8644,8645,8646,8647,8648,8649,8650,8651,8652,8653,8654,8655,8656,8657,8658,8659,8660,8661,8662,8663,8664,8665,8666,8667,8668,8935,8936,8937,8938,8939,8940,8941,8942,8943,8944,8945,8946,8947,8948,8949,8950,8951,8952,8953,8954,8955,8956,8957,8958,8959,8960,8961,8962,8963,8964,8965,9022,9023,9024,9025,9026,9027,9028,9029,9030,9031,9032,9033,9034,9035,9036,9037,9038,9039,9040,9041,9042,9043,9044,9045,9046,9047,9048,9049,9050,9051,9052,9321,9322,9323,9324,9325,9326,9327,9328,9329,9330,9331,9332,9333,9334,9335,9336,9337,9338,9339,9340,9341,9342,9343,9344,9345,9346,9347,9348,9349,9350,9351,9408,9409,9410,9411,9412,9413,9414,9415,9416,9417,9418,9419,9420,9421,9422,9423,9424,9425,9426,9427,9428,9429,9430,9431,9432,9433,9434,9435,9436,9437,9438,9709,9710,9711,9712,9713,9714,9715,9716,9717,9718,9719,9720,9721,9722,9723,9724,9725,9726,9727,9728,9729,9730,9731,9732,9733,9734,9735,9736,9737,9738,9739,9796,9797,9798,9799,9800,9801,9802,9803,9804,9805,9806,9807,9808,9809,9810,9811,9812,9813,9814,9815,9816,9817,9818,9819,9820,9821,9822,9823,9824,9825,9826,9827,10098,10099,10100,10101,10102,10103,10104,10105,10106,10107,10108,10109,10110,10111,10112,10113,10114,10115,10116,10117,10118,10119,10120,10121,10122,10123,10124,10125,10126,10127,10128,10129,10186,10187,10188,10189,10190,10191,10192,10193,10194,10195,10196,10197,10198,10199,10200,10201,10202,10203,10204,10205,10206,10207,10208,10209,10210,10211,10212,10213,10214,10215,10216,10217,10490,10491,10492,10493,10494,10495,10496,10497,10498,10499,10500,10501,10502,10503,10504,10505,10506,10507,10508,10509,10510,10511,10512,10513,10514,10515,10516,10517,10518,10519,10520,10521,10579,10580,10581,10582,10583,10584,10585,10586,10587,10588,10589,10590,10591,10592,10593,10594,10595,10596,10597,10598,10599,10600,10601,10602,10603,10604,10605,10606,10607,10608,10609,10884,10885,10886,10887,10888,10889,10890,10891,10892,10893,10894,10895,10896,10897,10898,10899,10900,10901,10902,10903,10904,10905,10906,10907,10908,10909,10910,10911,10912,10913,10914',
       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.T23t19.J184640+503000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T23t19.J184640+503000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T23t19.J184640+503000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#