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

# hif_editimlist(parameter_file='parameter.16305.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_T29t09.J215412+733000_P16305v1_2018_01_05T19_32_53.301/working/VLASS1.1.sb34639297.eb34672679.58065.06301421297_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T29t09.J215412+733000_P16305v1_2018_01_05T19_32_53.301/rawdata/VLASS1.1.sb34639297.eb34672679.58065.06301421297.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='7170,6723,6722,6721,6720,6727,5456,6725,6724,6842,6729,5358,7393,6836,5821,5820,5823,5822,5825,5824,5827,5826,5829,5828,6834,5926,5902,5903,5900,5901,5906,5907,5904,5905,7193,7192,5908,5909,7197,7196,7195,7194,7395,7394,7397,7396,7183,7390,7174,7392,7169,7399,5927,5461,5460,5463,5462,7401,6718,6719,6830,6843,6840,6841,6846,6844,6845,7406,6392,5928,6390,6391,7178,5929,5919,5918,5915,5914,5917,5916,5911,5910,5913,5912,6293,6292,6385,6384,6387,6386,6381,6380,6383,6382,6389,6388,6266,6267,6264,6265,6263,6268,6269,6375,7419,7402,6820,6821,6822,6823,6824,6825,6826,6827,6828,6829,7405,7404,6288,7179,7407,6289,6271,6270,6273,6272,6275,6274,6277,6276,6279,6278,7400,6378,6379,7403,7421,7420,6370,6371,6372,6373,6374,6286,6376,6377,6833,6832,6831,6287,6837,6728,6835,6818,6839,6838,5331,7190,5333,5332,5335,5334,5337,5336,5339,5338,5433,5436,5437,5434,5435,6284,6285,5438,5439,6280,6281,6282,6283,7412,7413,6369,6368,7416,7417,7414,7415,6363,6362,7418,7171,6367,6366,6365,6364,5449,5361,7198,7191,5344,5345,5346,5347,5340,5341,5342,5343,5348,5349,6291,6290,5818,5819,7175,7176,7409,7408,7410,5810,5811,5812,5813,5814,5815,5816,5817,6819,6726,7173,7185,7172,6817,6816,6745,6744,6747,6746,6741,6740,6743,6742,7391,6749,6748,5458,5356,5355,5354,5353,5352,5351,5350,5450,5451,5452,5453,5454,5455,5359,5457,5924,5925,5809,5808,5920,5921,5922,5923,5803,5802,5801,5800,5807,5806,5805,5804,5357,5459,7411,7398,6730,6731,6732,6733,6734,6735,6736,6737,6738,6739,5362,5360,5448,5830,5443,5442,5441,5440,5447,5446,5445,5444,7188,7189,5931,5930,7180,7181,7182,7177,7184,7168,7186,7187',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 21:54:12.182 +73.30.0.0000',
       scan=['129,132,135,138,142,113,116,119,122,126'],
       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_T29t09.J215412+733000_P16305v1_2018_01_05T19_32_53.301/working/VLASS1.1.sb34639297.eb34672679.58065.06301421297_split.ms'],
       imagename='VLASS1.1.ql.T29t09.J215412+733000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='5331,5332,5333,5334,5335,5336,5337,5338,5339,5340,5341,5342,5343,5344,5345,5346,5347,5348,5349,5350,5351,5352,5353,5354,5355,5356,5357,5358,5359,5360,5361,5362,5433,5434,5435,5436,5437,5438,5439,5440,5441,5442,5443,5444,5445,5446,5447,5448,5449,5450,5451,5452,5453,5454,5455,5456,5457,5458,5459,5460,5461,5462,5463,5800,5801,5802,5803,5804,5805,5806,5807,5808,5809,5810,5811,5812,5813,5814,5815,5816,5817,5818,5819,5820,5821,5822,5823,5824,5825,5826,5827,5828,5829,5830,5900,5901,5902,5903,5904,5905,5906,5907,5908,5909,5910,5911,5912,5913,5914,5915,5916,5917,5918,5919,5920,5921,5922,5923,5924,5925,5926,5927,5928,5929,5930,5931,6263,6264,6265,6266,6267,6268,6269,6270,6271,6272,6273,6274,6275,6276,6277,6278,6279,6280,6281,6282,6283,6284,6285,6286,6287,6288,6289,6290,6291,6292,6293,6362,6363,6364,6365,6366,6367,6368,6369,6370,6371,6372,6373,6374,6375,6376,6377,6378,6379,6380,6381,6382,6383,6384,6385,6386,6387,6388,6389,6390,6391,6392,6718,6719,6720,6721,6722,6723,6724,6725,6726,6727,6728,6729,6730,6731,6732,6733,6734,6735,6736,6737,6738,6739,6740,6741,6742,6743,6744,6745,6746,6747,6748,6749,6816,6817,6818,6819,6820,6821,6822,6823,6824,6825,6826,6827,6828,6829,6830,6831,6832,6833,6834,6835,6836,6837,6838,6839,6840,6841,6842,6843,6844,6845,6846,7168,7169,7170,7171,7172,7173,7174,7175,7176,7177,7178,7179,7180,7181,7182,7183,7184,7185,7186,7187,7188,7189,7190,7191,7192,7193,7194,7195,7196,7197,7198,7390,7391,7392,7393,7394,7395,7396,7397,7398,7399,7400,7401,7402,7403,7404,7405,7406,7407,7408,7409,7410,7411,7412,7413,7414,7415,7416,7417,7418,7419,7420,7421',
       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 21:54:12.182 +73.30.0.0000',
       scan=['129,132,135,138,142,113,116,119,122,126'],
       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_T29t09.J215412+733000_P16305v1_2018_01_05T19_32_53.301/working/VLASS1.1.sb34639297.eb34672679.58065.06301421297_split.ms'],
       imagename='VLASS1.1.ql.T29t09.J215412+733000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00052836648411557689, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='5331,5332,5333,5334,5335,5336,5337,5338,5339,5340,5341,5342,5343,5344,5345,5346,5347,5348,5349,5350,5351,5352,5353,5354,5355,5356,5357,5358,5359,5360,5361,5362,5433,5434,5435,5436,5437,5438,5439,5440,5441,5442,5443,5444,5445,5446,5447,5448,5449,5450,5451,5452,5453,5454,5455,5456,5457,5458,5459,5460,5461,5462,5463,5800,5801,5802,5803,5804,5805,5806,5807,5808,5809,5810,5811,5812,5813,5814,5815,5816,5817,5818,5819,5820,5821,5822,5823,5824,5825,5826,5827,5828,5829,5830,5900,5901,5902,5903,5904,5905,5906,5907,5908,5909,5910,5911,5912,5913,5914,5915,5916,5917,5918,5919,5920,5921,5922,5923,5924,5925,5926,5927,5928,5929,5930,5931,6263,6264,6265,6266,6267,6268,6269,6270,6271,6272,6273,6274,6275,6276,6277,6278,6279,6280,6281,6282,6283,6284,6285,6286,6287,6288,6289,6290,6291,6292,6293,6362,6363,6364,6365,6366,6367,6368,6369,6370,6371,6372,6373,6374,6375,6376,6377,6378,6379,6380,6381,6382,6383,6384,6385,6386,6387,6388,6389,6390,6391,6392,6718,6719,6720,6721,6722,6723,6724,6725,6726,6727,6728,6729,6730,6731,6732,6733,6734,6735,6736,6737,6738,6739,6740,6741,6742,6743,6744,6745,6746,6747,6748,6749,6816,6817,6818,6819,6820,6821,6822,6823,6824,6825,6826,6827,6828,6829,6830,6831,6832,6833,6834,6835,6836,6837,6838,6839,6840,6841,6842,6843,6844,6845,6846,7168,7169,7170,7171,7172,7173,7174,7175,7176,7177,7178,7179,7180,7181,7182,7183,7184,7185,7186,7187,7188,7189,7190,7191,7192,7193,7194,7195,7196,7197,7198,7390,7391,7392,7393,7394,7395,7396,7397,7398,7399,7400,7401,7402,7403,7404,7405,7406,7407,7408,7409,7410,7411,7412,7413,7414,7415,7416,7417,7418,7419,7420,7421',
       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.T29t09.J215412+733000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T29t09.J215412+733000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T29t09.J215412+733000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#