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

# hif_editimlist(parameter_file='parameter.13836.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_T23t16.J152055+483000_P13836v1_2018_01_18T22_44_55.995/working/VLASS1.1.sb34539903.eb34684984.58068.81223880787_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T23t16.J152055+483000_P13836v1_2018_01_18T22_44_55.995/rawdata/VLASS1.1.sb34539903.eb34684984.58068.81223880787.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='9230,9231,7632,7633,7630,7631,7634,8040,8041,7346,7347,7344,7345,7342,7343,7340,7341,7348,7349,8839,8838,8835,8834,8837,8836,8831,8830,8833,8832,9379,8561,9378,8560,9371,8569,9370,8568,8163,8162,8161,8160,8167,8166,8165,8164,8169,8168,8426,8427,8424,8425,8422,8423,8420,9369,9366,9367,9364,9365,9362,9363,8428,8429,8985,8984,8987,8986,8981,8980,8019,8018,8017,8016,8015,8014,8013,8012,8011,8010,7759,7758,7755,7754,7757,7756,7751,7750,7753,7752,8978,8979,8970,8971,8972,8973,8974,8975,8976,8977,9384,9385,9386,9387,9380,9381,9382,9383,7625,7624,7627,7626,7621,7620,7623,7622,7629,7628,9248,9241,9240,9243,9242,9245,9244,9247,9246,7351,7350,7353,7352,9238,9239,7357,8573,7359,7358,9236,9237,8578,8579,9232,9233,8181,8180,8182,7764,7765,7766,7767,7760,7761,7762,7763,7768,7769,8840,8841,8842,8843,8844,8845,8846,8847,7356,8574,8575,8576,8577,7355,7354,8572,8158,8159,8156,8157,8154,8155,8152,8153,9234,9235,7367,7618,7619,7614,7615,7616,7617,7610,7611,7612,7613,8431,8430,8433,8432,8435,8434,8437,8436,8439,8438,9373,9372,9375,9374,9377,9376,8028,8029,9360,8022,8023,8020,8021,8026,8027,8024,8025,8421,8983,8982,8444,8445,8440,8441,8442,8443,7368,7369,7360,7361,7362,7363,7364,7365,7366,8988,8969,8968,8963,8962,8961,8960,8967,8966,8965,8964,8817,8819,8818,8581,8580,8583,8582,8585,8584,8586,8567,8566,8565,8564,8563,8562,9229,9228,9227,9226,9225,9224,9223,9222,9221,9220,7777,7776,7775,7774,7773,7772,7771,7770,8570,8571,7339,9368,7609,7608,7607,7606,7605,7604,8035,8034,8037,8036,8031,8030,8033,8032,8039,8038,8957,8419,8958,8959,8418,9361,8416,8822,8823,8820,8821,8826,8827,8824,8825,8828,8829,8170,8171,8172,8173,8174,8175,8176,8177,8178,8179,9357,9359,9358,9218,9219,8556,8557,8558,8559,8417,9217,8415,7746,7747,7748,7749',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 15:20:55.360 +48.30.0.0000',
       scan=['160,163,166,138,141,144,147,170,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_T23t16.J152055+483000_P13836v1_2018_01_18T22_44_55.995/working/VLASS1.1.sb34539903.eb34684984.58068.81223880787_split.ms'],
       imagename='VLASS1.1.ql.T23t16.J152055+483000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='7339,7340,7341,7342,7343,7344,7345,7346,7347,7348,7349,7350,7351,7352,7353,7354,7355,7356,7357,7358,7359,7360,7361,7362,7363,7364,7365,7366,7367,7368,7369,7604,7605,7606,7607,7608,7609,7610,7611,7612,7613,7614,7615,7616,7617,7618,7619,7620,7621,7622,7623,7624,7625,7626,7627,7628,7629,7630,7631,7632,7633,7634,7746,7747,7748,7749,7750,7751,7752,7753,7754,7755,7756,7757,7758,7759,7760,7761,7762,7763,7764,7765,7766,7767,7768,7769,7770,7771,7772,7773,7774,7775,7776,7777,8010,8011,8012,8013,8014,8015,8016,8017,8018,8019,8020,8021,8022,8023,8024,8025,8026,8027,8028,8029,8030,8031,8032,8033,8034,8035,8036,8037,8038,8039,8040,8041,8152,8153,8154,8155,8156,8157,8158,8159,8160,8161,8162,8163,8164,8165,8166,8167,8168,8169,8170,8171,8172,8173,8174,8175,8176,8177,8178,8179,8180,8181,8182,8415,8416,8417,8418,8419,8420,8421,8422,8423,8424,8425,8426,8427,8428,8429,8430,8431,8432,8433,8434,8435,8436,8437,8438,8439,8440,8441,8442,8443,8444,8445,8556,8557,8558,8559,8560,8561,8562,8563,8564,8565,8566,8567,8568,8569,8570,8571,8572,8573,8574,8575,8576,8577,8578,8579,8580,8581,8582,8583,8584,8585,8586,8817,8818,8819,8820,8821,8822,8823,8824,8825,8826,8827,8828,8829,8830,8831,8832,8833,8834,8835,8836,8837,8838,8839,8840,8841,8842,8843,8844,8845,8846,8847,8957,8958,8959,8960,8961,8962,8963,8964,8965,8966,8967,8968,8969,8970,8971,8972,8973,8974,8975,8976,8977,8978,8979,8980,8981,8982,8983,8984,8985,8986,8987,8988,9217,9218,9219,9220,9221,9222,9223,9224,9225,9226,9227,9228,9229,9230,9231,9232,9233,9234,9235,9236,9237,9238,9239,9240,9241,9242,9243,9244,9245,9246,9247,9248,9357,9358,9359,9360,9361,9362,9363,9364,9365,9366,9367,9368,9369,9370,9371,9372,9373,9374,9375,9376,9377,9378,9379,9380,9381,9382,9383,9384,9385,9386,9387',
       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:20:55.360 +48.30.0.0000',
       scan=['160,163,166,138,141,144,147,170,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_T23t16.J152055+483000_P13836v1_2018_01_18T22_44_55.995/working/VLASS1.1.sb34539903.eb34684984.58068.81223880787_split.ms'],
       imagename='VLASS1.1.ql.T23t16.J152055+483000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00049806657461595022, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='7339,7340,7341,7342,7343,7344,7345,7346,7347,7348,7349,7350,7351,7352,7353,7354,7355,7356,7357,7358,7359,7360,7361,7362,7363,7364,7365,7366,7367,7368,7369,7604,7605,7606,7607,7608,7609,7610,7611,7612,7613,7614,7615,7616,7617,7618,7619,7620,7621,7622,7623,7624,7625,7626,7627,7628,7629,7630,7631,7632,7633,7634,7746,7747,7748,7749,7750,7751,7752,7753,7754,7755,7756,7757,7758,7759,7760,7761,7762,7763,7764,7765,7766,7767,7768,7769,7770,7771,7772,7773,7774,7775,7776,7777,8010,8011,8012,8013,8014,8015,8016,8017,8018,8019,8020,8021,8022,8023,8024,8025,8026,8027,8028,8029,8030,8031,8032,8033,8034,8035,8036,8037,8038,8039,8040,8041,8152,8153,8154,8155,8156,8157,8158,8159,8160,8161,8162,8163,8164,8165,8166,8167,8168,8169,8170,8171,8172,8173,8174,8175,8176,8177,8178,8179,8180,8181,8182,8415,8416,8417,8418,8419,8420,8421,8422,8423,8424,8425,8426,8427,8428,8429,8430,8431,8432,8433,8434,8435,8436,8437,8438,8439,8440,8441,8442,8443,8444,8445,8556,8557,8558,8559,8560,8561,8562,8563,8564,8565,8566,8567,8568,8569,8570,8571,8572,8573,8574,8575,8576,8577,8578,8579,8580,8581,8582,8583,8584,8585,8586,8817,8818,8819,8820,8821,8822,8823,8824,8825,8826,8827,8828,8829,8830,8831,8832,8833,8834,8835,8836,8837,8838,8839,8840,8841,8842,8843,8844,8845,8846,8847,8957,8958,8959,8960,8961,8962,8963,8964,8965,8966,8967,8968,8969,8970,8971,8972,8973,8974,8975,8976,8977,8978,8979,8980,8981,8982,8983,8984,8985,8986,8987,8988,9217,9218,9219,9220,9221,9222,9223,9224,9225,9226,9227,9228,9229,9230,9231,9232,9233,9234,9235,9236,9237,9238,9239,9240,9241,9242,9243,9244,9245,9246,9247,9248,9357,9358,9359,9360,9361,9362,9363,9364,9365,9366,9367,9368,9369,9370,9371,9372,9373,9374,9375,9376,9377,9378,9379,9380,9381,9382,9383,9384,9385,9386,9387',
       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.T23t16.J152055+483000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T23t16.J152055+483000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T23t16.J152055+483000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#