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

# hif_editimlist(parameter_file='parameter.5540.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_T27t08.J120605+663000_P5540v1_2017_11_17T01_19_22.291/working/VLASS1.1.sb34401340.eb34620284.58047.7253216088_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T27t08.J120605+663000_P5540v1_2017_11_17T01_19_22.291/rawdata/VLASS1.1.sb34401340.eb34620284.58047.7253216088.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='10109,10108,11275,10105,10104,10107,10106,10101,10100,10103,10102,11289,10899,10898,11288,10893,10892,10891,10890,10897,10896,10895,10894,9939,9938,9933,9932,9931,9930,9937,9936,9935,9934,11281,11280,11283,11282,11285,11274,11284,11270,11287,11286,11271,8934,9943,8936,8937,9155,9154,10709,9156,9528,9150,9153,9152,10703,10702,9526,9527,9159,10706,10705,10704,9551,9550,9553,9552,9555,9554,9557,9556,9328,9329,9322,9323,9320,9321,9326,9327,9324,9325,8941,8940,8943,8942,8945,8944,8947,8946,8949,8948,10093,10092,10091,10090,10097,10096,10095,10094,10099,10098,9942,8935,9940,9941,8931,8932,8933,8938,8939,11279,9718,10699,10698,9168,9169,9164,9165,9166,9167,9160,9161,9162,9163,9539,9538,9537,9536,10305,9534,9533,9708,9709,10886,10887,10880,10881,10882,10883,9700,9701,9702,9703,9704,9705,9706,9707,10879,10878,10871,10870,10873,10872,10875,10874,10877,10876,11120,11121,11122,11123,11124,11125,10901,10900,10309,10308,10307,10306,9535,10304,10303,9532,9531,9530,9529,9335,9334,9337,9336,9331,9330,9333,9332,10330,9339,9338,10714,10715,10716,10717,10710,10711,10499,10713,10497,10496,10495,10494,10493,10719,10491,10490,10475,10474,10477,10476,9712,10479,10478,10088,10089,10084,10085,10086,10087,10083,11119,11118,11111,11110,11113,11112,11115,11114,11117,11116,10725,10724,10727,10726,10721,10720,10723,10722,9175,9174,9173,9172,9171,9170,10332,10333,9157,11300,9719,10331,9713,9151,9711,9710,9717,9716,9715,9714,9913,9915,9914,9917,9916,9919,9918,10697,8961,8960,10701,10700,10712,10707,10498,9158,10718,10492,11298,11299,11292,11293,11290,11291,11296,11297,11294,11295,11095,11096,11097,11098,11099,10318,10319,10310,10311,10312,10313,10314,10315,10316,10317,9340,9341,10112,10488,10489,10480,10481,10482,10483,10484,10485,10486,10487,10501,10500,10503,10502,10505,10504,10884,10885,11276,9928,9929,11277,9920,9921,9922,9923,9924,9925,9926,9927,11272,11273,10888,10889,10110,10111,11108,11109,11102,11103,11100,11101,11106,11107,11104,11105,10321,10320,10323,10322,10325,10324,10327,10326,10329,10328,9148,9149,9146,9147,9144,9145,10113,9548,9549,9546,9547,9544,9545,9542,9543,9540,9541,9724,9725,9722,9723,9720,9721,9319,9318,9317,9316,9315,9314,9313,9312,9311,8956,8957,8954,8955,8952,8953,8950,8951,8958,8959,11278,10708,9698,9699,9696,9697,9695',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 12:06:5.343 +66.30.0.0000',
       scan=['160,188,194,163,166,169,172,175,179,182,185,156,191'],
       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_T27t08.J120605+663000_P5540v1_2017_11_17T01_19_22.291/working/VLASS1.1.sb34401340.eb34620284.58047.7253216088_split.ms'],
       imagename='VLASS1.1.ql.T27t08.J120605+663000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='8931,8932,8933,8934,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,9144,9145,9146,9147,9148,9149,9150,9151,9152,9153,9154,9155,9156,9157,9158,9159,9160,9161,9162,9163,9164,9165,9166,9167,9168,9169,9170,9171,9172,9173,9174,9175,9311,9312,9313,9314,9315,9316,9317,9318,9319,9320,9321,9322,9323,9324,9325,9326,9327,9328,9329,9330,9331,9332,9333,9334,9335,9336,9337,9338,9339,9340,9341,9526,9527,9528,9529,9530,9531,9532,9533,9534,9535,9536,9537,9538,9539,9540,9541,9542,9543,9544,9545,9546,9547,9548,9549,9550,9551,9552,9553,9554,9555,9556,9557,9695,9696,9697,9698,9699,9700,9701,9702,9703,9704,9705,9706,9707,9708,9709,9710,9711,9712,9713,9714,9715,9716,9717,9718,9719,9720,9721,9722,9723,9724,9725,9913,9914,9915,9916,9917,9918,9919,9920,9921,9922,9923,9924,9925,9926,9927,9928,9929,9930,9931,9932,9933,9934,9935,9936,9937,9938,9939,9940,9941,9942,9943,10083,10084,10085,10086,10087,10088,10089,10090,10091,10092,10093,10094,10095,10096,10097,10098,10099,10100,10101,10102,10103,10104,10105,10106,10107,10108,10109,10110,10111,10112,10113,10303,10304,10305,10306,10307,10308,10309,10310,10311,10312,10313,10314,10315,10316,10317,10318,10319,10320,10321,10322,10323,10324,10325,10326,10327,10328,10329,10330,10331,10332,10333,10474,10475,10476,10477,10478,10479,10480,10481,10482,10483,10484,10485,10486,10487,10488,10489,10490,10491,10492,10493,10494,10495,10496,10497,10498,10499,10500,10501,10502,10503,10504,10505,10697,10698,10699,10700,10701,10702,10703,10704,10705,10706,10707,10708,10709,10710,10711,10712,10713,10714,10715,10716,10717,10718,10719,10720,10721,10722,10723,10724,10725,10726,10727,10870,10871,10872,10873,10874,10875,10876,10877,10878,10879,10880,10881,10882,10883,10884,10885,10886,10887,10888,10889,10890,10891,10892,10893,10894,10895,10896,10897,10898,10899,10900,10901,11095,11096,11097,11098,11099,11100,11101,11102,11103,11104,11105,11106,11107,11108,11109,11110,11111,11112,11113,11114,11115,11116,11117,11118,11119,11120,11121,11122,11123,11124,11125,11270,11271,11272,11273,11274,11275,11276,11277,11278,11279,11280,11281,11282,11283,11284,11285,11286,11287,11288,11289,11290,11291,11292,11293,11294,11295,11296,11297,11298,11299,11300',
       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 12:06:5.343 +66.30.0.0000',
       scan=['160,188,194,163,166,169,172,175,179,182,185,156,191'],
       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_T27t08.J120605+663000_P5540v1_2017_11_17T01_19_22.291/working/VLASS1.1.sb34401340.eb34620284.58047.7253216088_split.ms'],
       imagename='VLASS1.1.ql.T27t08.J120605+663000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00049046056983233937, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='8931,8932,8933,8934,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,9144,9145,9146,9147,9148,9149,9150,9151,9152,9153,9154,9155,9156,9157,9158,9159,9160,9161,9162,9163,9164,9165,9166,9167,9168,9169,9170,9171,9172,9173,9174,9175,9311,9312,9313,9314,9315,9316,9317,9318,9319,9320,9321,9322,9323,9324,9325,9326,9327,9328,9329,9330,9331,9332,9333,9334,9335,9336,9337,9338,9339,9340,9341,9526,9527,9528,9529,9530,9531,9532,9533,9534,9535,9536,9537,9538,9539,9540,9541,9542,9543,9544,9545,9546,9547,9548,9549,9550,9551,9552,9553,9554,9555,9556,9557,9695,9696,9697,9698,9699,9700,9701,9702,9703,9704,9705,9706,9707,9708,9709,9710,9711,9712,9713,9714,9715,9716,9717,9718,9719,9720,9721,9722,9723,9724,9725,9913,9914,9915,9916,9917,9918,9919,9920,9921,9922,9923,9924,9925,9926,9927,9928,9929,9930,9931,9932,9933,9934,9935,9936,9937,9938,9939,9940,9941,9942,9943,10083,10084,10085,10086,10087,10088,10089,10090,10091,10092,10093,10094,10095,10096,10097,10098,10099,10100,10101,10102,10103,10104,10105,10106,10107,10108,10109,10110,10111,10112,10113,10303,10304,10305,10306,10307,10308,10309,10310,10311,10312,10313,10314,10315,10316,10317,10318,10319,10320,10321,10322,10323,10324,10325,10326,10327,10328,10329,10330,10331,10332,10333,10474,10475,10476,10477,10478,10479,10480,10481,10482,10483,10484,10485,10486,10487,10488,10489,10490,10491,10492,10493,10494,10495,10496,10497,10498,10499,10500,10501,10502,10503,10504,10505,10697,10698,10699,10700,10701,10702,10703,10704,10705,10706,10707,10708,10709,10710,10711,10712,10713,10714,10715,10716,10717,10718,10719,10720,10721,10722,10723,10724,10725,10726,10727,10870,10871,10872,10873,10874,10875,10876,10877,10878,10879,10880,10881,10882,10883,10884,10885,10886,10887,10888,10889,10890,10891,10892,10893,10894,10895,10896,10897,10898,10899,10900,10901,11095,11096,11097,11098,11099,11100,11101,11102,11103,11104,11105,11106,11107,11108,11109,11110,11111,11112,11113,11114,11115,11116,11117,11118,11119,11120,11121,11122,11123,11124,11125,11270,11271,11272,11273,11274,11275,11276,11277,11278,11279,11280,11281,11282,11283,11284,11285,11286,11287,11288,11289,11290,11291,11292,11293,11294,11295,11296,11297,11298,11299,11300',
       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.T27t08.J120605+663000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T27t08.J120605+663000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T27t08.J120605+663000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#