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

# hif_editimlist(parameter_file='parameter.list')
#

# hif_transformimagedata(datacolumn='corrected', modify_weights=False, clear_pointing=True)
#
mstransform(timebin='0s',
            outputvis='VLASS1.1.sb34898226.eb35042893.58162.29820586806_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T05t17.J111916-213000_P25121v1_2018_04_01T02_33_51.391/rawdata/VLASS1.1.sb34898226.eb35042893.58162.29820586806.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='1788,1789,3592,1785,1786,1787,3613,3614,3615,2383,2406,2407,2404,2382,2019,2018,2400,2401,2015,2014,2017,2016,2011,2010,2013,2012,1799,1798,2387,1793,1792,1791,1790,1797,1796,1795,1794,1397,1396,1395,3602,2798,1400,1825,1393,1826,1821,1820,1823,1822,1418,1419,2394,2395,1388,1389,2390,1824,1410,1411,1412,1413,1414,1415,1416,1417,3581,2411,2805,2408,2820,2821,3585,2785,2396,2787,2786,2781,2780,2783,2397,3575,3574,3577,3576,2789,2413,3610,1409,1408,2385,2384,1399,1398,1403,1402,1401,1394,1407,1392,1391,1390,3580,3214,3582,3583,3584,2417,3586,3215,3588,3589,2391,3611,3216,2419,3217,1385,3210,1386,3211,1387,2796,3612,3213,2797,3605,2398,2792,2399,2386,3593,2415,3591,3590,3597,3596,3599,3598,3219,2392,2808,2809,2793,2802,2803,2389,2806,2807,2804,2414,1421,1420,1423,1422,1425,1424,3212,1426,2412,3218,1404,2393,2405,2790,2388,2815,2788,2817,2816,2811,2810,2813,2812,2420,2421,2422,2423,2819,2818,3179,3178,2791,2416,3207,3206,3205,3204,3203,3202,3201,2799,3609,3209,3208,2024,2025,2020,2021,2022,2023,3200,2409,3603,2814,3601,3600,2418,3607,3606,1814,1815,1816,1817,1810,1811,1812,1813,3604,1818,1819,3184,3185,3186,3187,3180,3181,3182,3183,3193,3188,3189,3192,3608,2784,2795,3579,2794,3587,3578,1807,1804,1803,1802,1801,1800,2782,1809,1808,3197,3196,3195,3194,2006,2007,3191,3190,2008,2009,2410',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 11:19:16.837 -21.30.0.0000',
       scan=['33,59,36,65,40,43,46,49,52,62,55,27,30'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34898226.eb35042893.58162.29820586806_split.ms'],
       imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='1385,1386,1387,1388,1389,1390,1391,1392,1393,1394,1395,1396,1397,1398,1399,1400,1401,1402,1403,1404,1407,1408,1409,1410,1411,1412,1413,1414,1415,1416,1417,1418,1419,1420,1421,1422,1423,1424,1425,1426,1785,1786,1787,1788,1789,1790,1791,1792,1793,1794,1795,1796,1797,1798,1799,1800,1801,1802,1803,1804,1807,1808,1809,1810,1811,1812,1813,1814,1815,1816,1817,1818,1819,1820,1821,1822,1823,1824,1825,1826,2006,2007,2008,2009,2010,2011,2012,2013,2014,2015,2016,2017,2018,2019,2020,2021,2022,2023,2024,2025,2382,2383,2384,2385,2386,2387,2388,2389,2390,2391,2392,2393,2394,2395,2396,2397,2398,2399,2400,2401,2404,2405,2406,2407,2408,2409,2410,2411,2412,2413,2414,2415,2416,2417,2418,2419,2420,2421,2422,2423,2780,2781,2782,2783,2784,2785,2786,2787,2788,2789,2790,2791,2792,2793,2794,2795,2796,2797,2798,2799,2802,2803,2804,2805,2806,2807,2808,2809,2810,2811,2812,2813,2814,2815,2816,2817,2818,2819,2820,2821,3178,3179,3180,3181,3182,3183,3184,3185,3186,3187,3188,3189,3190,3191,3192,3193,3194,3195,3196,3197,3200,3201,3202,3203,3204,3205,3206,3207,3208,3209,3210,3211,3212,3213,3214,3215,3216,3217,3218,3219,3574,3575,3576,3577,3578,3579,3580,3581,3582,3583,3584,3585,3586,3587,3588,3589,3590,3591,3592,3593,3596,3597,3598,3599,3600,3601,3602,3603,3604,3605,3606,3607,3608,3609,3610,3611,3612,3613,3614,3615',
       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 11:19:16.837 -21.30.0.0000',
       scan=['33,59,36,65,40,43,46,49,52,62,55,27,30'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34898226.eb35042893.58162.29820586806_split.ms'],
       imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.0005596410670663317, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='1385,1386,1387,1388,1389,1390,1391,1392,1393,1394,1395,1396,1397,1398,1399,1400,1401,1402,1403,1404,1407,1408,1409,1410,1411,1412,1413,1414,1415,1416,1417,1418,1419,1420,1421,1422,1423,1424,1425,1426,1785,1786,1787,1788,1789,1790,1791,1792,1793,1794,1795,1796,1797,1798,1799,1800,1801,1802,1803,1804,1807,1808,1809,1810,1811,1812,1813,1814,1815,1816,1817,1818,1819,1820,1821,1822,1823,1824,1825,1826,2006,2007,2008,2009,2010,2011,2012,2013,2014,2015,2016,2017,2018,2019,2020,2021,2022,2023,2024,2025,2382,2383,2384,2385,2386,2387,2388,2389,2390,2391,2392,2393,2394,2395,2396,2397,2398,2399,2400,2401,2404,2405,2406,2407,2408,2409,2410,2411,2412,2413,2414,2415,2416,2417,2418,2419,2420,2421,2422,2423,2780,2781,2782,2783,2784,2785,2786,2787,2788,2789,2790,2791,2792,2793,2794,2795,2796,2797,2798,2799,2802,2803,2804,2805,2806,2807,2808,2809,2810,2811,2812,2813,2814,2815,2816,2817,2818,2819,2820,2821,3178,3179,3180,3181,3182,3183,3184,3185,3186,3187,3188,3189,3190,3191,3192,3193,3194,3195,3196,3197,3200,3201,3202,3203,3204,3205,3206,3207,3208,3209,3210,3211,3212,3213,3214,3215,3216,3217,3218,3219,3574,3575,3576,3577,3578,3579,3580,3581,3582,3583,3584,3585,3586,3587,3588,3589,3590,3591,3592,3593,3596,3597,3598,3599,3600,3601,3602,3603,3604,3605,3606,3607,3608,3609,3610,3611,3612,3613,3614,3615',
       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.T05t17.J111916-213000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T05t17.J111916-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.T05t17.J111916-213000.10.2048.v1.I.iter1.image.pbcor.tt0.rms',
      imagename='VLASS1.1.ql.T05t17.J111916-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.T05t17.J111916-213000.10.2048.v1.I.iter1.image.pbcor.tt0')
imhead(imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.image.pbcor.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.image.pbcor.tt0.subim',
           imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.image.pbcor.tt0')
imstat(imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imhead(imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.image.pbcor.tt0.rms.subim',
           imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imstat(imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.psf.tt0')
imhead(imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.psf.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.psf.tt0.subim',
           imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.psf.tt0')
imstat(imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.residual.tt0')
imhead(imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.residual.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.residual.tt0.subim',
           imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.residual.tt0')
imstat(imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.pb.tt0')
imhead(imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.pb.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.pb.tt0.subim',
           imagename='VLASS1.1.ql.T05t17.J111916-213000.10.2048.v1.I.iter1.pb.tt0')

# hifv_exportvlassdata(pipelinemode="automatic")
#