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

# hif_editimlist(parameter_file='parameter.8543.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_T20t20.J151930+373000_P8543v1_2017_11_05T04_31_09.833/working/VLASS1.1.sb34450715.eb34635774.58054.85418905092_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T20t20.J151930+373000_P8543v1_2017_11_05T04_31_09.833/rawdata/VLASS1.1.sb34450715.eb34635774.58054.85418905092.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='11542,11543,11540,11541,11546,11547,11544,11545,11548,11549,9939,9938,9933,9932,9931,9930,9937,9936,9935,9934,10332,10333,10330,10708,10336,10337,10334,10335,10703,10702,10338,10700,10707,10706,10705,10704,9559,9558,9551,9550,9553,9552,9555,9554,9557,9556,9616,9617,9614,10387,10386,10385,10384,10383,10382,10381,10380,9630,9612,10790,10389,9613,9610,9611,11539,9942,9943,9940,9941,9946,9947,9944,9945,9948,9949,11168,11169,11164,11165,11166,11167,11160,11161,11162,11163,10699,10698,10340,10697,9539,9560,9561,9562,9563,9564,9565,9566,9567,9568,9569,11089,11088,11087,11086,11551,11550,10309,11559,11558,11155,11154,11157,11156,11151,11150,11153,11152,11159,11158,10400,10401,10402,10390,10391,9629,9628,10394,10395,10396,10397,10013,10012,10011,10010,9627,10016,10015,10014,10714,10715,10716,10717,10710,10711,10712,10713,10718,10719,9622,11489,11488,11487,11486,11485,11484,11483,11482,11481,11480,9954,9951,9950,9953,9952,10392,10393,9626,10398,10399,9621,9620,11111,11110,11113,11112,11115,11114,11116,9625,9624,10728,10725,10724,10727,10726,10721,10720,10723,10722,9988,9989,9986,9987,9985,10759,10709,11568,11569,10331,11560,11561,11562,11563,11564,11565,11566,11567,10701,10339,11094,11095,11096,11097,11090,11091,11092,11093,11098,11099,10318,10319,10310,10311,10312,10313,10314,10315,10316,10317,11149,10004,10005,10006,10007,10000,10001,10002,10003,10008,10009,9618,9619,10761,10760,10763,10762,10765,10764,10767,10766,10769,10768,11498,11499,11553,11552,11555,11554,11557,11556,11490,11491,11492,11493,11494,11495,11496,11497,9928,9929,9923,9924,9925,9926,9927,11108,11109,11102,11103,11100,11101,11106,11107,11104,11105,11476,11477,11478,11479,10321,10320,10323,10322,10325,10324,10327,10326,10329,10328,9615,9999,9998,9991,9990,9993,9992,9995,9994,9997,9996,9548,9549,9546,9547,9544,9545,9542,9543,9540,9541,10783,10782,10781,10780,10787,10786,10785,10784,10789,10788,11506,9623,11504,11505,11502,11503,11500,11501,11173,11172,11171,11170,11177,11176,11175,11174,11179,11178,10388,9601,9600,9603,9602,9605,9604,9607,9606,9609,9608,10378,10379,10376,10377,10374,10375,10372,10373,10371,10772,10773,10770,10771,10776,10777,10774,10775,10778,10779',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 15:19:30.769 +37.30.0.0000',
       scan=['194,165,168,171,175,200,178,181,184,187,190,197'],
       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_T20t20.J151930+373000_P8543v1_2017_11_05T04_31_09.833/working/VLASS1.1.sb34450715.eb34635774.58054.85418905092_split.ms'],
       imagename='VLASS1.1.ql.T20t20.J151930+373000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='9539,9540,9541,9542,9543,9544,9545,9546,9547,9548,9549,9550,9551,9552,9553,9554,9555,9556,9557,9558,9559,9560,9561,9562,9563,9564,9565,9566,9567,9568,9569,9600,9601,9602,9603,9604,9605,9606,9607,9608,9609,9610,9611,9612,9613,9614,9615,9616,9617,9618,9619,9620,9621,9622,9623,9624,9625,9626,9627,9628,9629,9630,9923,9924,9925,9926,9927,9928,9929,9930,9931,9932,9933,9934,9935,9936,9937,9938,9939,9940,9941,9942,9943,9944,9945,9946,9947,9948,9949,9950,9951,9952,9953,9954,9985,9986,9987,9988,9989,9990,9991,9992,9993,9994,9995,9996,9997,9998,9999,10000,10001,10002,10003,10004,10005,10006,10007,10008,10009,10010,10011,10012,10013,10014,10015,10016,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,10334,10335,10336,10337,10338,10339,10340,10371,10372,10373,10374,10375,10376,10377,10378,10379,10380,10381,10382,10383,10384,10385,10386,10387,10388,10389,10390,10391,10392,10393,10394,10395,10396,10397,10398,10399,10400,10401,10402,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,10728,10759,10760,10761,10762,10763,10764,10765,10766,10767,10768,10769,10770,10771,10772,10773,10774,10775,10776,10777,10778,10779,10780,10781,10782,10783,10784,10785,10786,10787,10788,10789,10790,11086,11087,11088,11089,11090,11091,11092,11093,11094,11095,11096,11097,11098,11099,11100,11101,11102,11103,11104,11105,11106,11107,11108,11109,11110,11111,11112,11113,11114,11115,11116,11149,11150,11151,11152,11153,11154,11155,11156,11157,11158,11159,11160,11161,11162,11163,11164,11165,11166,11167,11168,11169,11170,11171,11172,11173,11174,11175,11176,11177,11178,11179,11476,11477,11478,11479,11480,11481,11482,11483,11484,11485,11486,11487,11488,11489,11490,11491,11492,11493,11494,11495,11496,11497,11498,11499,11500,11501,11502,11503,11504,11505,11506,11539,11540,11541,11542,11543,11544,11545,11546,11547,11548,11549,11550,11551,11552,11553,11554,11555,11556,11557,11558,11559,11560,11561,11562,11563,11564,11565,11566,11567,11568,11569',
       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 15:19:30.769 +37.30.0.0000',
       scan=['194,165,168,171,175,200,178,181,184,187,190,197'],
       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_T20t20.J151930+373000_P8543v1_2017_11_05T04_31_09.833/working/VLASS1.1.sb34450715.eb34635774.58054.85418905092_split.ms'],
       imagename='VLASS1.1.ql.T20t20.J151930+373000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00051680446365317038, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='9539,9540,9541,9542,9543,9544,9545,9546,9547,9548,9549,9550,9551,9552,9553,9554,9555,9556,9557,9558,9559,9560,9561,9562,9563,9564,9565,9566,9567,9568,9569,9600,9601,9602,9603,9604,9605,9606,9607,9608,9609,9610,9611,9612,9613,9614,9615,9616,9617,9618,9619,9620,9621,9622,9623,9624,9625,9626,9627,9628,9629,9630,9923,9924,9925,9926,9927,9928,9929,9930,9931,9932,9933,9934,9935,9936,9937,9938,9939,9940,9941,9942,9943,9944,9945,9946,9947,9948,9949,9950,9951,9952,9953,9954,9985,9986,9987,9988,9989,9990,9991,9992,9993,9994,9995,9996,9997,9998,9999,10000,10001,10002,10003,10004,10005,10006,10007,10008,10009,10010,10011,10012,10013,10014,10015,10016,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,10334,10335,10336,10337,10338,10339,10340,10371,10372,10373,10374,10375,10376,10377,10378,10379,10380,10381,10382,10383,10384,10385,10386,10387,10388,10389,10390,10391,10392,10393,10394,10395,10396,10397,10398,10399,10400,10401,10402,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,10728,10759,10760,10761,10762,10763,10764,10765,10766,10767,10768,10769,10770,10771,10772,10773,10774,10775,10776,10777,10778,10779,10780,10781,10782,10783,10784,10785,10786,10787,10788,10789,10790,11086,11087,11088,11089,11090,11091,11092,11093,11094,11095,11096,11097,11098,11099,11100,11101,11102,11103,11104,11105,11106,11107,11108,11109,11110,11111,11112,11113,11114,11115,11116,11149,11150,11151,11152,11153,11154,11155,11156,11157,11158,11159,11160,11161,11162,11163,11164,11165,11166,11167,11168,11169,11170,11171,11172,11173,11174,11175,11176,11177,11178,11179,11476,11477,11478,11479,11480,11481,11482,11483,11484,11485,11486,11487,11488,11489,11490,11491,11492,11493,11494,11495,11496,11497,11498,11499,11500,11501,11502,11503,11504,11505,11506,11539,11540,11541,11542,11543,11544,11545,11546,11547,11548,11549,11550,11551,11552,11553,11554,11555,11556,11557,11558,11559,11560,11561,11562,11563,11564,11565,11566,11567,11568,11569',
       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.T20t20.J151930+373000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T20t20.J151930+373000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T20t20.J151930+373000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#