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

# hif_editimlist(parameter_file='parameter.7170.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_T16t23.J165333+213000_P7170v1_2017_09_29T21_21_34.655/working/VLASS1.1.sb34441820.eb34483341.58017.956089363426_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T16t23.J165333+213000_P7170v1_2017_09_29T21_21_34.655/rawdata/VLASS1.1.sb34441820.eb34483341.58017.956089363426.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='9674,9675,10518,9677,9670,9671,9672,9673,10105,10513,10107,10106,9678,9679,10103,10102,9472,9473,8626,8627,9476,9474,9475,10964,10963,10962,10961,10960,9676,10519,10109,10108,8608,8609,8602,8603,8606,8607,8604,8605,10512,10104,10510,10511,10101,10100,10514,10515,10332,10333,10330,10331,10336,10337,10334,10335,9256,9257,9254,9255,9252,9253,9250,9251,10097,10096,10095,10099,10098,8835,8834,8837,8836,8831,8830,8833,8832,10738,10938,10939,10934,10935,10936,10937,10933,9449,9448,9447,9446,10739,10736,8389,8388,8387,8386,8385,8384,8383,8382,8381,8380,9689,9688,10749,10748,9681,9680,9683,9682,9685,9684,9687,9686,10527,10526,10525,10524,10523,10522,10521,10520,10529,10528,10118,10119,10116,10117,10114,10115,10112,10113,10110,10111,9902,9900,9901,8806,8807,8808,8809,9249,9248,10309,10308,10307,9243,9245,9244,9247,9246,9878,9879,9876,9877,9874,9875,9872,9873,8615,8614,8617,8616,8611,8610,8613,8612,8619,8618,8378,8379,10737,10516,9058,9059,9050,9051,9052,9053,9054,9055,9056,9057,10517,9894,9895,9896,9897,9890,9891,9892,9893,9898,9899,10949,10948,10941,10940,10943,10942,10945,10944,10947,10946,8628,8629,8624,8625,9470,9471,8620,8621,8622,8623,9270,9271,9272,9273,10758,10759,10750,10751,10752,10753,10754,10755,10756,10757,10530,10531,10532,10533,10534,10535,10536,10537,10538,10539,10123,10122,10121,10120,10125,10124,9258,9259,10338,10540,8817,8816,8815,8814,8813,8812,8811,8810,8819,8818,10318,10319,10310,10311,10312,10313,10314,10315,10316,10317,10747,10746,10745,10744,10743,10742,10741,10740,9060,8408,8409,8400,8401,8402,8403,8404,8405,8406,8407,10761,10760,10762,10509,9049,9048,10322,9043,9042,9041,9040,9047,9046,9045,9044,9887,9886,9885,9884,9883,9882,9881,9880,9889,9888,10958,10959,9669,10952,10953,10950,10951,10956,10957,10954,10955,8633,8632,8631,8630,10321,9262,9261,9260,9267,9266,9265,9264,10329,10328,9269,9268,10732,10733,10731,9465,9464,9467,9466,9461,9460,9463,9462,9469,9468,9263,10320,10323,9032,9033,9030,9031,9036,9037,9034,9035,10325,9038,9039,10324,10327,10326,8822,8823,8820,8821,8826,8827,8824,8825,8828,8829,10734,10735,9458,9459,9454,9455,9456,9457,9450,9451,9452,9453,8398,8399,8394,8395,8396,8397,8390,8391,8392,8393,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 16:53:33.911 +21.30.0.0000',
       scan=['172,163,166,134,137,140,143,147,150,169,153,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_T16t23.J165333+213000_P7170v1_2017_09_29T21_21_34.655/working/VLASS1.1.sb34441820.eb34483341.58017.956089363426_split.ms'],
       imagename='VLASS1.1.ql.T16t23.J165333+213000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='8378,8379,8380,8381,8382,8383,8384,8385,8386,8387,8388,8389,8390,8391,8392,8393,8394,8395,8396,8397,8398,8399,8400,8401,8402,8403,8404,8405,8406,8407,8408,8409,8602,8603,8604,8605,8606,8607,8608,8609,8610,8611,8612,8613,8614,8615,8616,8617,8618,8619,8620,8621,8622,8623,8624,8625,8626,8627,8628,8629,8630,8631,8632,8633,8806,8807,8808,8809,8810,8811,8812,8813,8814,8815,8816,8817,8818,8819,8820,8821,8822,8823,8824,8825,8826,8827,8828,8829,8830,8831,8832,8833,8834,8835,8836,8837,9030,9031,9032,9033,9034,9035,9036,9037,9038,9039,9040,9041,9042,9043,9044,9045,9046,9047,9048,9049,9050,9051,9052,9053,9054,9055,9056,9057,9058,9059,9060,9243,9244,9245,9246,9247,9248,9249,9250,9251,9252,9253,9254,9255,9256,9257,9258,9259,9260,9261,9262,9263,9264,9265,9266,9267,9268,9269,9270,9271,9272,9273,9446,9447,9448,9449,9450,9451,9452,9453,9454,9455,9456,9457,9458,9459,9460,9461,9462,9463,9464,9465,9466,9467,9468,9469,9470,9471,9472,9473,9474,9475,9476,9669,9670,9671,9672,9673,9674,9675,9676,9677,9678,9679,9680,9681,9682,9683,9684,9685,9686,9687,9688,9689,9690,9691,9692,9693,9694,9695,9696,9697,9698,9699,9872,9873,9874,9875,9876,9877,9878,9879,9880,9881,9882,9883,9884,9885,9886,9887,9888,9889,9890,9891,9892,9893,9894,9895,9896,9897,9898,9899,9900,9901,9902,10095,10096,10097,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,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,10334,10335,10336,10337,10338,10509,10510,10511,10512,10513,10514,10515,10516,10517,10518,10519,10520,10521,10522,10523,10524,10525,10526,10527,10528,10529,10530,10531,10532,10533,10534,10535,10536,10537,10538,10539,10540,10731,10732,10733,10734,10735,10736,10737,10738,10739,10740,10741,10742,10743,10744,10745,10746,10747,10748,10749,10750,10751,10752,10753,10754,10755,10756,10757,10758,10759,10760,10761,10762,10933,10934,10935,10936,10937,10938,10939,10940,10941,10942,10943,10944,10945,10946,10947,10948,10949,10950,10951,10952,10953,10954,10955,10956,10957,10958,10959,10960,10961,10962,10963,10964',
       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 16:53:33.911 +21.30.0.0000',
       scan=['172,163,166,134,137,140,143,147,150,169,153,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_T16t23.J165333+213000_P7170v1_2017_09_29T21_21_34.655/working/VLASS1.1.sb34441820.eb34483341.58017.956089363426_split.ms'],
       imagename='VLASS1.1.ql.T16t23.J165333+213000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00051544366429643555, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='8378,8379,8380,8381,8382,8383,8384,8385,8386,8387,8388,8389,8390,8391,8392,8393,8394,8395,8396,8397,8398,8399,8400,8401,8402,8403,8404,8405,8406,8407,8408,8409,8602,8603,8604,8605,8606,8607,8608,8609,8610,8611,8612,8613,8614,8615,8616,8617,8618,8619,8620,8621,8622,8623,8624,8625,8626,8627,8628,8629,8630,8631,8632,8633,8806,8807,8808,8809,8810,8811,8812,8813,8814,8815,8816,8817,8818,8819,8820,8821,8822,8823,8824,8825,8826,8827,8828,8829,8830,8831,8832,8833,8834,8835,8836,8837,9030,9031,9032,9033,9034,9035,9036,9037,9038,9039,9040,9041,9042,9043,9044,9045,9046,9047,9048,9049,9050,9051,9052,9053,9054,9055,9056,9057,9058,9059,9060,9243,9244,9245,9246,9247,9248,9249,9250,9251,9252,9253,9254,9255,9256,9257,9258,9259,9260,9261,9262,9263,9264,9265,9266,9267,9268,9269,9270,9271,9272,9273,9446,9447,9448,9449,9450,9451,9452,9453,9454,9455,9456,9457,9458,9459,9460,9461,9462,9463,9464,9465,9466,9467,9468,9469,9470,9471,9472,9473,9474,9475,9476,9669,9670,9671,9672,9673,9674,9675,9676,9677,9678,9679,9680,9681,9682,9683,9684,9685,9686,9687,9688,9689,9690,9691,9692,9693,9694,9695,9696,9697,9698,9699,9872,9873,9874,9875,9876,9877,9878,9879,9880,9881,9882,9883,9884,9885,9886,9887,9888,9889,9890,9891,9892,9893,9894,9895,9896,9897,9898,9899,9900,9901,9902,10095,10096,10097,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,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,10334,10335,10336,10337,10338,10509,10510,10511,10512,10513,10514,10515,10516,10517,10518,10519,10520,10521,10522,10523,10524,10525,10526,10527,10528,10529,10530,10531,10532,10533,10534,10535,10536,10537,10538,10539,10540,10731,10732,10733,10734,10735,10736,10737,10738,10739,10740,10741,10742,10743,10744,10745,10746,10747,10748,10749,10750,10751,10752,10753,10754,10755,10756,10757,10758,10759,10760,10761,10762,10933,10934,10935,10936,10937,10938,10939,10940,10941,10942,10943,10944,10945,10946,10947,10948,10949,10950,10951,10952,10953,10954,10955,10956,10957,10958,10959,10960,10961,10962,10963,10964',
       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.T16t23.J165333+213000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T16t23.J165333+213000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T16t23.J165333+213000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#