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

# hif_editimlist(parameter_file='parameter.5720.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_T28t06.J104300+713000_P5720v1_2017_12_20T15_55_08.830/working/VLASS1.1.sb34417859.eb34434020.58014.453630601856_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T28t06.J104300+713000_P5720v1_2017_12_20T15_55_08.830/rawdata/VLASS1.1.sb34417859.eb34434020.58014.453630601856.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='9076,9077,9074,9075,9072,9073,9070,9071,9078,9079,10109,10108,10105,10104,10107,10106,10101,10100,10103,10102,8097,8099,8098,9854,8667,9843,9842,9841,9840,9847,9846,9845,9844,8695,8126,8125,9848,8691,8690,8121,8692,8350,8349,8348,8351,8343,8342,9320,9321,9326,9327,8345,9325,8941,8940,8353,8354,9312,9311,8357,9856,10248,10249,10242,10243,10240,10241,10246,10247,10244,10245,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,9322,9323,9838,9839,9833,9836,9837,9834,9835,8493,8492,8491,8490,8497,8496,8495,8494,8499,8498,9061,9060,9063,9062,9065,9064,9067,9066,9069,9068,9708,9709,10118,9700,9701,8338,8339,9704,9705,9706,9707,8685,10233,10232,9324,10237,10236,10235,10234,10239,10238,8501,8500,8503,8502,8505,8504,8118,8119,8112,8113,8110,8111,8116,8117,8114,8115,7980,7981,7982,7983,7984,7985,7986,7987,7988,7989,8479,8478,8475,8477,8476,9058,9304,9305,9050,9051,9052,9053,9054,9306,9056,9057,9307,10088,10089,9300,10087,9301,8927,8926,8925,8924,8923,8922,8921,8920,9303,8929,8928,7962,7963,7966,7967,7964,7965,7968,7969,9719,9718,9308,9309,9713,9712,9711,9710,9717,9716,9715,9714,8687,8127,10115,8694,9849,8124,8123,8122,8693,8120,8105,8104,8107,8106,8101,8100,8103,8102,8109,8108,9861,9860,9863,9862,10114,9059,7993,7992,7991,7990,8665,8666,9439,8668,8669,8361,8360,8363,8362,8365,8364,8367,8366,8369,8368,9055,10260,10261,10262,9049,9298,9299,8912,8913,8910,8911,8916,8917,8914,8915,8918,8919,10116,10117,9702,9703,10112,10113,10110,10111,9328,9329,8682,9851,8680,8681,8686,7978,8684,9857,7975,7974,7977,7976,7971,7970,7973,7972,8340,8347,8346,9465,9464,9467,9466,9461,9460,9463,9462,8344,9469,9468,9720,9721,8352,9080,8358,8359,9319,9318,9317,9316,9315,9314,9313,8355,8356,9310,10259,10258,10251,10250,10253,10252,10255,10254,10257,10256,9850,8683,9852,9302,9853,7979,9855,9458,9459,9454,9455,9456,9457,9450,9451,9452,9453,9858,9859,8688,8689,8673,8672,8671,8670,8677,8676,8675,8674,8679,8678,8486,8487,8488,8489,8341,8480,8481,8482,8483,8484,8485,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 10:43:0.105 +71.30.0.0000',
       scan=['160,163,170,135,138,167,141,173,144,148,151,154,157'],
       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_T28t06.J104300+713000_P5720v1_2017_12_20T15_55_08.830/working/VLASS1.1.sb34417859.eb34434020.58014.453630601856_split.ms'],
       imagename='VLASS1.1.ql.T28t06.J104300+713000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='7962,7963,7964,7965,7966,7967,7968,7969,7970,7971,7972,7973,7974,7975,7976,7977,7978,7979,7980,7981,7982,7983,7984,7985,7986,7987,7988,7989,7990,7991,7992,7993,8097,8098,8099,8100,8101,8102,8103,8104,8105,8106,8107,8108,8109,8110,8111,8112,8113,8114,8115,8116,8117,8118,8119,8120,8121,8122,8123,8124,8125,8126,8127,8338,8339,8340,8341,8342,8343,8344,8345,8346,8347,8348,8349,8350,8351,8352,8353,8354,8355,8356,8357,8358,8359,8360,8361,8362,8363,8364,8365,8366,8367,8368,8369,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,8665,8666,8667,8668,8669,8670,8671,8672,8673,8674,8675,8676,8677,8678,8679,8680,8681,8682,8683,8684,8685,8686,8687,8688,8689,8690,8691,8692,8693,8694,8695,8910,8911,8912,8913,8914,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,9049,9050,9051,9052,9053,9054,9055,9056,9057,9058,9059,9060,9061,9062,9063,9064,9065,9066,9067,9068,9069,9070,9071,9072,9073,9074,9075,9076,9077,9078,9079,9080,9298,9299,9300,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,9439,9440,9441,9442,9443,9444,9445,9446,9447,9448,9449,9450,9451,9452,9453,9454,9455,9456,9457,9458,9459,9460,9461,9462,9463,9464,9465,9466,9467,9468,9469,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,9720,9721,9833,9834,9835,9836,9837,9838,9839,9840,9841,9842,9843,9844,9845,9846,9847,9848,9849,9850,9851,9852,9853,9854,9855,9856,9857,9858,9859,9860,9861,9862,9863,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,10114,10115,10116,10117,10118,10232,10233,10234,10235,10236,10237,10238,10239,10240,10241,10242,10243,10244,10245,10246,10247,10248,10249,10250,10251,10252,10253,10254,10255,10256,10257,10258,10259,10260,10261,10262',
       outframe='LSRK', gridder='mosaic', cycleniter=100, 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 10:43:0.105 +71.30.0.0000',
       scan=['160,163,170,135,138,167,141,173,144,148,151,154,157'],
       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_T28t06.J104300+713000_P5720v1_2017_12_20T15_55_08.830/working/VLASS1.1.sb34417859.eb34434020.58014.453630601856_split.ms'],
       imagename='VLASS1.1.ql.T28t06.J104300+713000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00072544545782147899, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='7962,7963,7964,7965,7966,7967,7968,7969,7970,7971,7972,7973,7974,7975,7976,7977,7978,7979,7980,7981,7982,7983,7984,7985,7986,7987,7988,7989,7990,7991,7992,7993,8097,8098,8099,8100,8101,8102,8103,8104,8105,8106,8107,8108,8109,8110,8111,8112,8113,8114,8115,8116,8117,8118,8119,8120,8121,8122,8123,8124,8125,8126,8127,8338,8339,8340,8341,8342,8343,8344,8345,8346,8347,8348,8349,8350,8351,8352,8353,8354,8355,8356,8357,8358,8359,8360,8361,8362,8363,8364,8365,8366,8367,8368,8369,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,8665,8666,8667,8668,8669,8670,8671,8672,8673,8674,8675,8676,8677,8678,8679,8680,8681,8682,8683,8684,8685,8686,8687,8688,8689,8690,8691,8692,8693,8694,8695,8910,8911,8912,8913,8914,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,9049,9050,9051,9052,9053,9054,9055,9056,9057,9058,9059,9060,9061,9062,9063,9064,9065,9066,9067,9068,9069,9070,9071,9072,9073,9074,9075,9076,9077,9078,9079,9080,9298,9299,9300,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,9439,9440,9441,9442,9443,9444,9445,9446,9447,9448,9449,9450,9451,9452,9453,9454,9455,9456,9457,9458,9459,9460,9461,9462,9463,9464,9465,9466,9467,9468,9469,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,9720,9721,9833,9834,9835,9836,9837,9838,9839,9840,9841,9842,9843,9844,9845,9846,9847,9848,9849,9850,9851,9852,9853,9854,9855,9856,9857,9858,9859,9860,9861,9862,9863,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,10114,10115,10116,10117,10118,10232,10233,10234,10235,10236,10237,10238,10239,10240,10241,10242,10243,10244,10245,10246,10247,10248,10249,10250,10251,10252,10253,10254,10255,10256,10257,10258,10259,10260,10261,10262',
       outframe='LSRK', gridder='mosaic', cycleniter=100, 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.T28t06.J104300+713000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T28t06.J104300+713000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T28t06.J104300+713000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#