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.4728.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.J184026+503000_P4728v1_2017_09_15T20_53_53.901/working/VLASS1.1.sb34361741.eb34367683.58005.085977442126_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T23t19.J184026+503000_P4728v1_2017_09_15T20_53_53.901/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='9072,9070,9071,10109,10108,10105,10104,10107,10106,10101,10100,10103,10102,10893,10892,10891,10890,10894,10207,10208,10209,9438,9439,9843,9842,9841,9840,9847,9846,9845,9844,9328,9329,9322,9323,9320,9321,9326,9327,9324,9325,8941,8940,8943,8942,8945,8944,8946,10093,10092,10091,10090,10097,10096,10095,10094,10099,10098,8934,8935,8936,8937,8930,8931,8932,8933,8938,8939,9449,9448,9447,9446,9445,9444,9443,9442,9441,9440,9825,9824,8549,8548,8545,8544,8547,8546,8541,8540,8543,8542,9838,9839,9832,9829,9830,9831,9836,9837,9834,9828,9689,9061,9060,9063,9062,9065,9064,9067,9066,9069,9068,9708,9709,10886,10887,10880,10881,10882,10883,9700,9701,9702,9703,9704,9705,9706,9707,10599,10233,10232,10231,10230,10237,10236,10235,10234,10879,10878,10871,10870,10873,10872,10875,10874,10877,10876,9716,9715,10500,10628,10629,10620,10621,10622,10623,10624,10625,10626,10627,9331,9330,10499,10498,10497,10496,10495,10494,10493,10492,10491,10490,10475,9304,10477,10476,10471,10470,10473,9712,9050,9051,9052,9053,10479,9306,9056,9057,9710,10088,10089,10084,10085,10086,10087,10080,10081,10082,10083,8927,8926,8925,8924,8923,8922,8921,8920,9714,8929,8928,9055,8659,8658,8657,8538,8539,9827,9826,9821,9820,9823,9822,8532,8533,8534,8535,8536,8537,10611,10610,10613,10612,10615,10614,10617,10616,10619,10618,9719,9718,9308,9309,9713,9305,9711,9307,9717,9301,9302,9303,10228,10229,10224,10225,10226,10227,10220,10221,10222,10223,10868,10869,9833,10866,10867,10864,10865,9835,9058,9059,8660,8661,8662,8663,8664,8665,8666,8667,9436,9437,9434,9435,9432,9433,9430,9431,8668,10472,8562,8561,8560,8669,9818,9819,9054,9816,9817,10478,10488,10489,10480,10481,10482,10483,10484,10485,10486,10487,9049,9048,9043,9042,9047,9046,9045,9044,10884,10885,8916,8917,8915,8918,8919,10215,10214,10217,10216,10211,10210,10213,10212,10219,10218,9429,8672,10888,10889,8682,8683,8680,8681,8686,8687,8684,8685,8688,10602,10603,10600,10601,10606,10607,10604,10605,10608,10609,9319,9318,9317,9316,9315,9314,9313,9312,9311,9310,9458,9454,9455,9456,9457,9450,9451,9452,9453,10474,8673,9428,8671,8670,8677,8676,8675,8674,8679,8678,10079,10078,8552,8553,8550,8551,8556,8557,8554,8555,8558,8559,9698,9699,9696,9697,9694,9695,9692,9693,9690,9691',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 18:40:26.929 +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.J184026+503000_P4728v1_2017_09_15T20_53_53.901/working/VLASS1.1.sb34361741.eb34367683.58005.085977442126_split.ms'],
       imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='8532,8533,8534,8535,8536,8537,8538,8539,8540,8541,8542,8543,8544,8545,8546,8547,8548,8549,8550,8551,8552,8553,8554,8555,8556,8557,8558,8559,8560,8561,8562,8657,8658,8659,8660,8661,8662,8663,8664,8665,8666,8667,8668,8669,8670,8671,8672,8673,8674,8675,8676,8677,8678,8679,8680,8681,8682,8683,8684,8685,8686,8687,8688,8915,8916,8917,8918,8919,8920,8921,8922,8923,8924,8925,8926,8927,8928,8929,8930,8931,8932,8933,8934,8935,8936,8937,8938,8939,8940,8941,8942,8943,8944,8945,8946,9042,9043,9044,9045,9046,9047,9048,9049,9050,9051,9052,9053,9054,9055,9056,9057,9058,9059,9060,9061,9062,9063,9064,9065,9066,9067,9068,9069,9070,9071,9072,9301,9302,9303,9304,9305,9306,9307,9308,9309,9310,9311,9312,9313,9314,9315,9316,9317,9318,9319,9320,9321,9322,9323,9324,9325,9326,9327,9328,9329,9330,9331,9428,9429,9430,9431,9432,9433,9434,9435,9436,9437,9438,9439,9440,9441,9442,9443,9444,9445,9446,9447,9448,9449,9450,9451,9452,9453,9454,9455,9456,9457,9458,9689,9690,9691,9692,9693,9694,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,9816,9817,9818,9819,9820,9821,9822,9823,9824,9825,9826,9827,9828,9829,9830,9831,9832,9833,9834,9835,9836,9837,9838,9839,9840,9841,9842,9843,9844,9845,9846,9847,10078,10079,10080,10081,10082,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,10207,10208,10209,10210,10211,10212,10213,10214,10215,10216,10217,10218,10219,10220,10221,10222,10223,10224,10225,10226,10227,10228,10229,10230,10231,10232,10233,10234,10235,10236,10237,10470,10471,10472,10473,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,10599,10600,10601,10602,10603,10604,10605,10606,10607,10608,10609,10610,10611,10612,10613,10614,10615,10616,10617,10618,10619,10620,10621,10622,10623,10624,10625,10626,10627,10628,10629,10864,10865,10866,10867,10868,10869,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',
       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:40:26.929 +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.J184026+503000_P4728v1_2017_09_15T20_53_53.901/working/VLASS1.1.sb34361741.eb34367683.58005.085977442126_split.ms'],
       imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00049151342786317399, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='8532,8533,8534,8535,8536,8537,8538,8539,8540,8541,8542,8543,8544,8545,8546,8547,8548,8549,8550,8551,8552,8553,8554,8555,8556,8557,8558,8559,8560,8561,8562,8657,8658,8659,8660,8661,8662,8663,8664,8665,8666,8667,8668,8669,8670,8671,8672,8673,8674,8675,8676,8677,8678,8679,8680,8681,8682,8683,8684,8685,8686,8687,8688,8915,8916,8917,8918,8919,8920,8921,8922,8923,8924,8925,8926,8927,8928,8929,8930,8931,8932,8933,8934,8935,8936,8937,8938,8939,8940,8941,8942,8943,8944,8945,8946,9042,9043,9044,9045,9046,9047,9048,9049,9050,9051,9052,9053,9054,9055,9056,9057,9058,9059,9060,9061,9062,9063,9064,9065,9066,9067,9068,9069,9070,9071,9072,9301,9302,9303,9304,9305,9306,9307,9308,9309,9310,9311,9312,9313,9314,9315,9316,9317,9318,9319,9320,9321,9322,9323,9324,9325,9326,9327,9328,9329,9330,9331,9428,9429,9430,9431,9432,9433,9434,9435,9436,9437,9438,9439,9440,9441,9442,9443,9444,9445,9446,9447,9448,9449,9450,9451,9452,9453,9454,9455,9456,9457,9458,9689,9690,9691,9692,9693,9694,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,9816,9817,9818,9819,9820,9821,9822,9823,9824,9825,9826,9827,9828,9829,9830,9831,9832,9833,9834,9835,9836,9837,9838,9839,9840,9841,9842,9843,9844,9845,9846,9847,10078,10079,10080,10081,10082,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,10207,10208,10209,10210,10211,10212,10213,10214,10215,10216,10217,10218,10219,10220,10221,10222,10223,10224,10225,10226,10227,10228,10229,10230,10231,10232,10233,10234,10235,10236,10237,10470,10471,10472,10473,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,10599,10600,10601,10602,10603,10604,10605,10606,10607,10608,10609,10610,10611,10612,10613,10614,10615,10616,10617,10618,10619,10620,10621,10622,10623,10624,10625,10626,10627,10628,10629,10864,10865,10866,10867,10868,10869,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',
       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.J184026+503000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T23t19.J184026+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.J184026+503000.10.2048.v1.I.iter1.image.pbcor.tt0.rms',
      imagename='VLASS1.1.ql.T23t19.J184026+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.J184026+503000.10.2048.v1.I.iter1.image.pbcor.tt0')
imhead(imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.image.pbcor.tt0')
imsubimage(box='1785.0,1785.0,5506.0,5506.0',
           outfile='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.image.pbcor.tt0.subim',
           imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.image.pbcor.tt0')
imstat(imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imhead(imagename='VLASS1.1.ql.T23t19.J184026+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.J184026+503000.10.2048.v1.I.iter1.image.pbcor.tt0.rms.subim',
           imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imstat(imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.psf.tt0')
imhead(imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.psf.tt0')
imsubimage(box='1785.0,1785.0,5506.0,5506.0',
           outfile='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.psf.tt0.subim',
           imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.psf.tt0')
imstat(imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.residual.tt0')
imhead(imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.residual.tt0')
imsubimage(box='1785.0,1785.0,5506.0,5506.0',
           outfile='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.residual.tt0.subim',
           imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.residual.tt0')
imstat(imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.pb.tt0')
imhead(imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.pb.tt0')
imsubimage(box='1785.0,1785.0,5506.0,5506.0',
           outfile='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.pb.tt0.subim',
           imagename='VLASS1.1.ql.T23t19.J184026+503000.10.2048.v1.I.iter1.pb.tt0')

# hifv_exportvlassdata(pipelinemode="automatic")
#