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.8565.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.J153000+393000_P8565v1_2017_11_05T04_31_20.175/working/VLASS1.1.sb34450715.eb34635774.58054.85418905092_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T20t20.J153000+393000_P8565v1_2017_11_05T04_31_20.175/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='6619,7002,6848,6849,6842,6843,6840,6841,6846,6847,6844,6845,6989,6988,6987,6986,6985,6984,6983,6982,6981,6980,8516,8517,8514,8515,8512,8513,8510,8511,8518,8519,8129,8128,8127,8126,8125,8124,8123,8122,8121,8120,6613,6612,6611,6610,6617,6616,6615,6614,7001,7000,7003,6618,7005,7004,7007,7006,8381,8380,7221,7220,7223,7222,7225,7224,7227,7226,7229,7228,7380,7381,7386,7387,7384,7385,7599,7598,7597,7596,7595,7594,7593,7592,7759,7758,8495,8494,7755,7754,7757,7756,7751,7750,7753,7752,6855,6854,6857,6856,6851,6850,6853,6852,6859,6858,7621,7620,7622,8524,8501,8500,8503,8502,8505,8504,8507,8506,8509,8508,8118,8119,8116,8117,8114,8115,7980,7981,7982,7983,7984,7985,7986,7987,7988,7989,8376,8377,8374,8375,8372,8373,8370,8371,8378,8379,7357,7356,7359,7358,7218,7219,7215,7216,7217,7764,7765,7760,7761,7762,7763,6628,6629,6626,6627,6624,6625,6622,6623,6620,6621,6868,6869,6864,6865,6866,6867,6860,6861,6862,6863,7618,7619,7614,7615,7616,7617,7610,7611,7612,7613,7382,7383,7368,7369,7245,7244,7243,7242,7241,7240,7360,7361,7362,7363,7364,7365,7366,7367,8497,8496,8499,8498,7993,7992,7991,7990,7997,7996,7995,7994,7999,7998,8361,8360,8363,8362,8365,8364,8367,8366,8369,8368,8141,8140,8143,8142,8144,6631,6630,6633,6632,6635,6634,6998,6999,6994,6995,6996,6997,6990,6991,6992,6993,7609,7608,7607,7606,7605,7604,7603,7602,7601,7600,8134,7979,7978,8523,7974,8521,8520,7971,8131,7973,7972,8135,7739,7738,7737,7736,7735,8136,8358,8359,8351,8352,8353,8354,8355,8356,8357,7379,7378,8137,7373,7372,7371,7370,7377,7376,7375,7374,6839,6838,6605,6606,6607,6608,6609,7010,8138,8139,7975,8522,7977,7976,8130,8132,8133,7236,7237,7234,7235,7232,7233,7230,7231,8000,8001,7238,7239,7009,7008,7742,7743,7740,7741,7746,7747,7744,7745,7748,7749',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 15:30:0.878 +39.30.0.0000',
       scan=['130,133,137,140,143,146,149,118,121,124,127'],
       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.J153000+393000_P8565v1_2017_11_05T04_31_20.175/working/VLASS1.1.sb34450715.eb34635774.58054.85418905092_split.ms'],
       imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='6605,6606,6607,6608,6609,6610,6611,6612,6613,6614,6615,6616,6617,6618,6619,6620,6621,6622,6623,6624,6625,6626,6627,6628,6629,6630,6631,6632,6633,6634,6635,6838,6839,6840,6841,6842,6843,6844,6845,6846,6847,6848,6849,6850,6851,6852,6853,6854,6855,6856,6857,6858,6859,6860,6861,6862,6863,6864,6865,6866,6867,6868,6869,6980,6981,6982,6983,6984,6985,6986,6987,6988,6989,6990,6991,6992,6993,6994,6995,6996,6997,6998,6999,7000,7001,7002,7003,7004,7005,7006,7007,7008,7009,7010,7215,7216,7217,7218,7219,7220,7221,7222,7223,7224,7225,7226,7227,7228,7229,7230,7231,7232,7233,7234,7235,7236,7237,7238,7239,7240,7241,7242,7243,7244,7245,7356,7357,7358,7359,7360,7361,7362,7363,7364,7365,7366,7367,7368,7369,7370,7371,7372,7373,7374,7375,7376,7377,7378,7379,7380,7381,7382,7383,7384,7385,7386,7387,7592,7593,7594,7595,7596,7597,7598,7599,7600,7601,7602,7603,7604,7605,7606,7607,7608,7609,7610,7611,7612,7613,7614,7615,7616,7617,7618,7619,7620,7621,7622,7735,7736,7737,7738,7739,7740,7741,7742,7743,7744,7745,7746,7747,7748,7749,7750,7751,7752,7753,7754,7755,7756,7757,7758,7759,7760,7761,7762,7763,7764,7765,7971,7972,7973,7974,7975,7976,7977,7978,7979,7980,7981,7982,7983,7984,7985,7986,7987,7988,7989,7990,7991,7992,7993,7994,7995,7996,7997,7998,7999,8000,8001,8114,8115,8116,8117,8118,8119,8120,8121,8122,8123,8124,8125,8126,8127,8128,8129,8130,8131,8132,8133,8134,8135,8136,8137,8138,8139,8140,8141,8142,8143,8144,8351,8352,8353,8354,8355,8356,8357,8358,8359,8360,8361,8362,8363,8364,8365,8366,8367,8368,8369,8370,8371,8372,8373,8374,8375,8376,8377,8378,8379,8380,8381,8494,8495,8496,8497,8498,8499,8500,8501,8502,8503,8504,8505,8506,8507,8508,8509,8510,8511,8512,8513,8514,8515,8516,8517,8518,8519,8520,8521,8522,8523,8524',
       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:30:0.878 +39.30.0.0000',
       scan=['130,133,137,140,143,146,149,118,121,124,127'],
       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.J153000+393000_P8565v1_2017_11_05T04_31_20.175/working/VLASS1.1.sb34450715.eb34635774.58054.85418905092_split.ms'],
       imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00048478307756231516, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='6605,6606,6607,6608,6609,6610,6611,6612,6613,6614,6615,6616,6617,6618,6619,6620,6621,6622,6623,6624,6625,6626,6627,6628,6629,6630,6631,6632,6633,6634,6635,6838,6839,6840,6841,6842,6843,6844,6845,6846,6847,6848,6849,6850,6851,6852,6853,6854,6855,6856,6857,6858,6859,6860,6861,6862,6863,6864,6865,6866,6867,6868,6869,6980,6981,6982,6983,6984,6985,6986,6987,6988,6989,6990,6991,6992,6993,6994,6995,6996,6997,6998,6999,7000,7001,7002,7003,7004,7005,7006,7007,7008,7009,7010,7215,7216,7217,7218,7219,7220,7221,7222,7223,7224,7225,7226,7227,7228,7229,7230,7231,7232,7233,7234,7235,7236,7237,7238,7239,7240,7241,7242,7243,7244,7245,7356,7357,7358,7359,7360,7361,7362,7363,7364,7365,7366,7367,7368,7369,7370,7371,7372,7373,7374,7375,7376,7377,7378,7379,7380,7381,7382,7383,7384,7385,7386,7387,7592,7593,7594,7595,7596,7597,7598,7599,7600,7601,7602,7603,7604,7605,7606,7607,7608,7609,7610,7611,7612,7613,7614,7615,7616,7617,7618,7619,7620,7621,7622,7735,7736,7737,7738,7739,7740,7741,7742,7743,7744,7745,7746,7747,7748,7749,7750,7751,7752,7753,7754,7755,7756,7757,7758,7759,7760,7761,7762,7763,7764,7765,7971,7972,7973,7974,7975,7976,7977,7978,7979,7980,7981,7982,7983,7984,7985,7986,7987,7988,7989,7990,7991,7992,7993,7994,7995,7996,7997,7998,7999,8000,8001,8114,8115,8116,8117,8118,8119,8120,8121,8122,8123,8124,8125,8126,8127,8128,8129,8130,8131,8132,8133,8134,8135,8136,8137,8138,8139,8140,8141,8142,8143,8144,8351,8352,8353,8354,8355,8356,8357,8358,8359,8360,8361,8362,8363,8364,8365,8366,8367,8368,8369,8370,8371,8372,8373,8374,8375,8376,8377,8378,8379,8380,8381,8494,8495,8496,8497,8498,8499,8500,8501,8502,8503,8504,8505,8506,8507,8508,8509,8510,8511,8512,8513,8514,8515,8516,8517,8518,8519,8520,8521,8522,8523,8524',
       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.J153000+393000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T20t20.J153000+393000.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.J153000+393000.10.2048.v1.I.iter1.image.pbcor.tt0.rms',
      imagename='VLASS1.1.ql.T20t20.J153000+393000.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.J153000+393000.10.2048.v1.I.iter1.image.pbcor.tt0')
imhead(imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.image.pbcor.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.image.pbcor.tt0.subim',
           imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.image.pbcor.tt0')
imstat(imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imhead(imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.image.pbcor.tt0.rms.subim',
           imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imstat(imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.psf.tt0')
imhead(imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.psf.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.psf.tt0.subim',
           imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.psf.tt0')
imstat(imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.residual.tt0')
imhead(imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.residual.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.residual.tt0.subim',
           imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.residual.tt0')
imstat(imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.pb.tt0')
imhead(imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.pb.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.pb.tt0.subim',
           imagename='VLASS1.1.ql.T20t20.J153000+393000.10.2048.v1.I.iter1.pb.tt0')

# hifv_exportvlassdata(pipelinemode="automatic")
#