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.sb34688349.eb34756507.58091.09663136574.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='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T13t06.J032201+103000_P17956v1_2018_05_04T01_35_25.868/working/VLASS1.1.sb34688349.eb34756507.58091.09663136574_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T13t06.J032201+103000_P17956v1_2018_05_04T01_35_25.868/rawdata/VLASS1.1.sb34688349.eb34756507.58091.09663136574.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,2580,2980,2616,1780,1781,1782,1783,1784,1785,1786,1787,3416,3417,3414,2614,3412,3021,3410,3411,2613,3020,3418,2612,1775,1774,1777,1776,2611,1773,2610,1779,1778,1373,1375,1374,1377,1376,1379,1378,2997,1798,1793,1792,1791,1790,1797,1796,1795,1794,3401,1403,3403,3402,3405,3024,3407,1396,3409,3408,1401,2998,3005,3396,2179,3379,3017,3004,1825,1407,2992,1826,1821,1820,1823,1406,1418,1419,3006,1388,1391,1824,1410,1411,1412,1390,1414,1415,1416,1417,2993,3003,2994,3007,3002,2991,2995,3022,2996,1822,3419,3386,3387,2989,2988,3382,3383,3380,3381,2983,2982,2981,2603,2987,2986,2985,2984,2173,1409,1408,2176,2177,2174,1398,1397,1402,1395,1394,1393,1392,1405,1404,2604,2605,2606,2607,2601,2602,2201,3385,2608,2202,2609,1384,2203,1385,2204,1386,1389,3391,3390,3018,1413,3395,3394,2578,2579,2576,1380,2574,2207,3016,2175,3014,1381,1382,3019,3398,1383,2594,2595,2596,2597,2590,2591,2592,2593,3023,2598,2619,2618,2190,2191,2192,2193,2194,2195,2196,2197,2198,3406,3009,3008,1421,1420,1423,1422,1425,1424,2206,1426,2208,2209,3388,3392,2587,2586,2585,2584,2583,2582,2581,2178,3413,3389,2589,2588,2183,2182,2181,2180,2187,2186,2185,2184,3393,2189,2188,2213,2212,2211,2210,2217,2216,2215,2214,3377,2990,2219,2218,3373,3375,3012,2577,2621,2222,3010,2622,2575,2205,2623,3374,2620,2617,2573,1387,3425,3415,2626,2223,3015,1814,1815,1816,1817,1810,1811,1812,1813,2226,2224,2225,1818,1819,2220,2221,3376,3378,3424,3013,2615,3384,3397,2978,2979,2624,2973,2976,2977,2974,2975,1807,1806,1805,1804,1803,1802,1801,3011,3423,3422,3421,3420,3426,1809,1808,2625,3026,3001,3025,3404',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 03:22:1.851 +10.30.0.0000',
       scan=['32,35,38,41,44,48,51,54,57,60,29,63'],
       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_T13t06.J032201+103000_P17956v1_2018_05_04T01_35_25.868/working/VLASS1.1.sb34688349.eb34756507.58091.09663136574_split.ms'],
       imagename='VLASS1.1.ql.T13t06.J032201+103000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='1373,1374,1375,1376,1377,1378,1379,1380,1381,1382,1383,1384,1385,1386,1387,1388,1389,1390,1391,1392,1393,1394,1395,1396,1397,1398,1401,1402,1403,1404,1405,1406,1407,1408,1409,1410,1411,1412,1413,1414,1415,1416,1417,1418,1419,1420,1421,1422,1423,1424,1425,1426,1773,1774,1775,1776,1777,1778,1779,1780,1781,1782,1783,1784,1785,1786,1787,1788,1789,1790,1791,1792,1793,1794,1795,1796,1797,1798,1801,1802,1803,1804,1805,1806,1807,1808,1809,1810,1811,1812,1813,1814,1815,1816,1817,1818,1819,1820,1821,1822,1823,1824,1825,1826,2173,2174,2175,2176,2177,2178,2179,2180,2181,2182,2183,2184,2185,2186,2187,2188,2189,2190,2191,2192,2193,2194,2195,2196,2197,2198,2201,2202,2203,2204,2205,2206,2207,2208,2209,2210,2211,2212,2213,2214,2215,2216,2217,2218,2219,2220,2221,2222,2223,2224,2225,2226,2573,2574,2575,2576,2577,2578,2579,2580,2581,2582,2583,2584,2585,2586,2587,2588,2589,2590,2591,2592,2593,2594,2595,2596,2597,2598,2601,2602,2603,2604,2605,2606,2607,2608,2609,2610,2611,2612,2613,2614,2615,2616,2617,2618,2619,2620,2621,2622,2623,2624,2625,2626,2973,2974,2975,2976,2977,2978,2979,2980,2981,2982,2983,2984,2985,2986,2987,2988,2989,2990,2991,2992,2993,2994,2995,2996,2997,2998,3001,3002,3003,3004,3005,3006,3007,3008,3009,3010,3011,3012,3013,3014,3015,3016,3017,3018,3019,3020,3021,3022,3023,3024,3025,3026,3373,3374,3375,3376,3377,3378,3379,3380,3381,3382,3383,3384,3385,3386,3387,3388,3389,3390,3391,3392,3393,3394,3395,3396,3397,3398,3401,3402,3403,3404,3405,3406,3407,3408,3409,3410,3411,3412,3413,3414,3415,3416,3417,3418,3419,3420,3421,3422,3423,3424,3425,3426',
       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 03:22:1.851 +10.30.0.0000',
       scan=['32,35,38,41,44,48,51,54,57,60,29,63'],
       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_T13t06.J032201+103000_P17956v1_2018_05_04T01_35_25.868/working/VLASS1.1.sb34688349.eb34756507.58091.09663136574_split.ms'],
       imagename='VLASS1.1.ql.T13t06.J032201+103000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00054977381914795838, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='1373,1374,1375,1376,1377,1378,1379,1380,1381,1382,1383,1384,1385,1386,1387,1388,1389,1390,1391,1392,1393,1394,1395,1396,1397,1398,1401,1402,1403,1404,1405,1406,1407,1408,1409,1410,1411,1412,1413,1414,1415,1416,1417,1418,1419,1420,1421,1422,1423,1424,1425,1426,1773,1774,1775,1776,1777,1778,1779,1780,1781,1782,1783,1784,1785,1786,1787,1788,1789,1790,1791,1792,1793,1794,1795,1796,1797,1798,1801,1802,1803,1804,1805,1806,1807,1808,1809,1810,1811,1812,1813,1814,1815,1816,1817,1818,1819,1820,1821,1822,1823,1824,1825,1826,2173,2174,2175,2176,2177,2178,2179,2180,2181,2182,2183,2184,2185,2186,2187,2188,2189,2190,2191,2192,2193,2194,2195,2196,2197,2198,2201,2202,2203,2204,2205,2206,2207,2208,2209,2210,2211,2212,2213,2214,2215,2216,2217,2218,2219,2220,2221,2222,2223,2224,2225,2226,2573,2574,2575,2576,2577,2578,2579,2580,2581,2582,2583,2584,2585,2586,2587,2588,2589,2590,2591,2592,2593,2594,2595,2596,2597,2598,2601,2602,2603,2604,2605,2606,2607,2608,2609,2610,2611,2612,2613,2614,2615,2616,2617,2618,2619,2620,2621,2622,2623,2624,2625,2626,2973,2974,2975,2976,2977,2978,2979,2980,2981,2982,2983,2984,2985,2986,2987,2988,2989,2990,2991,2992,2993,2994,2995,2996,2997,2998,3001,3002,3003,3004,3005,3006,3007,3008,3009,3010,3011,3012,3013,3014,3015,3016,3017,3018,3019,3020,3021,3022,3023,3024,3025,3026,3373,3374,3375,3376,3377,3378,3379,3380,3381,3382,3383,3384,3385,3386,3387,3388,3389,3390,3391,3392,3393,3394,3395,3396,3397,3398,3401,3402,3403,3404,3405,3406,3407,3408,3409,3410,3411,3412,3413,3414,3415,3416,3417,3418,3419,3420,3421,3422,3423,3424,3425,3426',
       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.T13t06.J032201+103000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T13t06.J032201+103000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T13t06.J032201+103000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#