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

# hif_editimlist(parameter_file='parameter.4865.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_T25t14.J182600+593000_P4865v1_2017_10_10T20_52_52.309/working/VLASS1.1.sb34395681.eb34403276.58010.036995381946_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T25t14.J182600+593000_P4865v1_2017_10_10T20_52_52.309/rawdata/VLASS1.1.sb34395681.eb34403276.58010.036995381946.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='9779,9778,9775,9777,9776,8909,8908,8905,8904,8907,8906,8901,8900,8903,8902,10206,10207,10204,10205,10202,10203,10200,10201,10208,10209,10196,10197,10194,10195,10192,10193,10190,10191,10198,10199,9726,9727,10136,10137,10130,9723,9720,10133,10977,10639,10638,10637,10636,10635,10634,10633,10632,10631,10630,9731,9730,10143,9732,10145,10144,10147,10146,9317,9320,10563,10562,10561,10560,10564,8891,8959,9368,9369,9366,9367,9364,9365,9363,8984,8549,8548,8983,8982,8547,8493,8492,8491,8490,8497,8496,8495,8494,8499,8498,10967,10966,10965,9708,9709,10964,10118,10119,10116,10117,9702,9703,9704,9705,9706,9707,10961,10960,8978,8979,9711,8970,8971,8972,8973,8974,8975,8976,8977,9710,10969,10968,9384,9385,9386,9387,9380,9381,9382,9383,10978,10979,9388,9389,9717,10123,9716,10618,9715,8501,8500,8503,8502,8505,8504,8506,9714,10628,10629,10620,10621,10622,10623,10624,10625,10626,10627,8574,8575,8576,8577,8570,8571,8572,8573,8475,8477,8476,8578,10142,10129,10128,9304,9712,9306,10120,10127,10126,10125,10124,8479,8478,8969,9783,10973,9379,9378,9371,9370,9373,9372,9375,9374,9377,9376,8963,8962,10611,10610,10613,10612,10615,10614,10617,10616,10619,8981,8980,10534,10535,10536,10537,10538,10539,9719,9718,9308,9309,9713,10122,10121,9307,9300,9301,9302,9303,8888,8889,8884,8885,8886,8887,8881,8882,8883,9780,9781,9782,8968,9784,9785,9786,9787,9788,9789,8961,8960,8967,8966,8965,8964,10545,10544,10547,10546,10541,10540,10543,10542,10220,10548,9393,9392,9391,9390,10985,10984,10981,10980,10983,10982,10549,10221,8567,8566,8565,8564,8563,8562,8561,8560,8569,8568,9292,9293,9290,9291,9296,9297,9294,9295,9298,9299,8910,8911,10141,10140,10215,10214,10217,10216,10211,10210,10213,10212,10219,10218,10958,10959,10956,10957,10955,10962,10609,10134,10135,9724,9725,9722,10131,10132,9721,10138,10139,9728,9729,9319,9318,9305,9316,9315,9314,9313,9312,9311,9310,8956,8957,8954,8955,8899,8898,8897,8896,8895,8894,8893,8892,8958,8890,9793,9792,9791,9790,9797,9796,9795,9794,9799,9798,10556,10557,10554,10555,10552,10553,10550,10551,10558,10559,10963,10970,10971,10972,8552,8553,8550,8551,8556,8557,8554,8555,10974,8558,8559,10975,10976,9805,9804,9803,9802,9801,9800,8480,8481,8482,8483,8484,8485,8486,8487,8488,8489',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 18:26:0.872 +59.30.0.0000',
       scan=['162,165,168,137,171,140,143,146,149,152,175,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_T25t14.J182600+593000_P4865v1_2017_10_10T20_52_52.309/working/VLASS1.1.sb34395681.eb34403276.58010.036995381946_split.ms'],
       imagename='VLASS1.1.ql.T25t14.J182600+593000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='8475,8476,8477,8478,8479,8480,8481,8482,8483,8484,8485,8486,8487,8488,8489,8490,8491,8492,8493,8494,8495,8496,8497,8498,8499,8500,8501,8502,8503,8504,8505,8506,8547,8548,8549,8550,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,8881,8882,8883,8884,8885,8886,8887,8888,8889,8890,8891,8892,8893,8894,8895,8896,8897,8898,8899,8900,8901,8902,8903,8904,8905,8906,8907,8908,8909,8910,8911,8954,8955,8956,8957,8958,8959,8960,8961,8962,8963,8964,8965,8966,8967,8968,8969,8970,8971,8972,8973,8974,8975,8976,8977,8978,8979,8980,8981,8982,8983,8984,9290,9291,9292,9293,9294,9295,9296,9297,9298,9299,9300,9301,9302,9303,9304,9305,9306,9307,9308,9309,9310,9311,9312,9313,9314,9315,9316,9317,9318,9319,9320,9363,9364,9365,9366,9367,9368,9369,9370,9371,9372,9373,9374,9375,9376,9377,9378,9379,9380,9381,9382,9383,9384,9385,9386,9387,9388,9389,9390,9391,9392,9393,9702,9703,9704,9705,9706,9707,9708,9709,9710,9711,9712,9713,9714,9715,9716,9717,9718,9719,9720,9721,9722,9723,9724,9725,9726,9727,9728,9729,9730,9731,9732,9775,9776,9777,9778,9779,9780,9781,9782,9783,9784,9785,9786,9787,9788,9789,9790,9791,9792,9793,9794,9795,9796,9797,9798,9799,9800,9801,9802,9803,9804,9805,10116,10117,10118,10119,10120,10121,10122,10123,10124,10125,10126,10127,10128,10129,10130,10131,10132,10133,10134,10135,10136,10137,10138,10139,10140,10141,10142,10143,10144,10145,10146,10147,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,10218,10219,10220,10221,10534,10535,10536,10537,10538,10539,10540,10541,10542,10543,10544,10545,10546,10547,10548,10549,10550,10551,10552,10553,10554,10555,10556,10557,10558,10559,10560,10561,10562,10563,10564,10609,10610,10611,10612,10613,10614,10615,10616,10617,10618,10619,10620,10621,10622,10623,10624,10625,10626,10627,10628,10629,10630,10631,10632,10633,10634,10635,10636,10637,10638,10639,10955,10956,10957,10958,10959,10960,10961,10962,10963,10964,10965,10966,10967,10968,10969,10970,10971,10972,10973,10974,10975,10976,10977,10978,10979,10980,10981,10982,10983,10984,10985',
       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:26:0.872 +59.30.0.0000',
       scan=['162,165,168,137,171,140,143,146,149,152,175,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_T25t14.J182600+593000_P4865v1_2017_10_10T20_52_52.309/working/VLASS1.1.sb34395681.eb34403276.58010.036995381946_split.ms'],
       imagename='VLASS1.1.ql.T25t14.J182600+593000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00050167542150020381, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='8475,8476,8477,8478,8479,8480,8481,8482,8483,8484,8485,8486,8487,8488,8489,8490,8491,8492,8493,8494,8495,8496,8497,8498,8499,8500,8501,8502,8503,8504,8505,8506,8547,8548,8549,8550,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,8881,8882,8883,8884,8885,8886,8887,8888,8889,8890,8891,8892,8893,8894,8895,8896,8897,8898,8899,8900,8901,8902,8903,8904,8905,8906,8907,8908,8909,8910,8911,8954,8955,8956,8957,8958,8959,8960,8961,8962,8963,8964,8965,8966,8967,8968,8969,8970,8971,8972,8973,8974,8975,8976,8977,8978,8979,8980,8981,8982,8983,8984,9290,9291,9292,9293,9294,9295,9296,9297,9298,9299,9300,9301,9302,9303,9304,9305,9306,9307,9308,9309,9310,9311,9312,9313,9314,9315,9316,9317,9318,9319,9320,9363,9364,9365,9366,9367,9368,9369,9370,9371,9372,9373,9374,9375,9376,9377,9378,9379,9380,9381,9382,9383,9384,9385,9386,9387,9388,9389,9390,9391,9392,9393,9702,9703,9704,9705,9706,9707,9708,9709,9710,9711,9712,9713,9714,9715,9716,9717,9718,9719,9720,9721,9722,9723,9724,9725,9726,9727,9728,9729,9730,9731,9732,9775,9776,9777,9778,9779,9780,9781,9782,9783,9784,9785,9786,9787,9788,9789,9790,9791,9792,9793,9794,9795,9796,9797,9798,9799,9800,9801,9802,9803,9804,9805,10116,10117,10118,10119,10120,10121,10122,10123,10124,10125,10126,10127,10128,10129,10130,10131,10132,10133,10134,10135,10136,10137,10138,10139,10140,10141,10142,10143,10144,10145,10146,10147,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,10218,10219,10220,10221,10534,10535,10536,10537,10538,10539,10540,10541,10542,10543,10544,10545,10546,10547,10548,10549,10550,10551,10552,10553,10554,10555,10556,10557,10558,10559,10560,10561,10562,10563,10564,10609,10610,10611,10612,10613,10614,10615,10616,10617,10618,10619,10620,10621,10622,10623,10624,10625,10626,10627,10628,10629,10630,10631,10632,10633,10634,10635,10636,10637,10638,10639,10955,10956,10957,10958,10959,10960,10961,10962,10963,10964,10965,10966,10967,10968,10969,10970,10971,10972,10973,10974,10975,10976,10977,10978,10979,10980,10981,10982,10983,10984,10985',
       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.T25t14.J182600+593000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T25t14.J182600+593000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T25t14.J182600+593000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#