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.sb34641006.eb34808190.58115.235410231486.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_T14t14.J090647+143000_P16458v1_2018_04_30T17_28_34.285/working/VLASS1.1.sb34641006.eb34808190.58115.235410231486_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T14t14.J090647+143000_P16458v1_2018_04_30T17_28_34.285/rawdata/VLASS1.1.sb34641006.eb34808190.58115.235410231486.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='1700,1701,1702,1703,1704,1705,1706,1707,1708,1709,3683,3682,3681,3680,3687,3686,3678,3679,3685,3672,3673,3670,3671,3676,3677,3674,3675,3029,3028,3027,3026,3025,3024,3023,3022,3021,3020,1457,1849,1848,1843,1842,1841,1840,1847,1846,1845,1844,2617,2616,2615,2614,2613,2612,2611,2610,2619,2618,3294,3295,3290,3291,3292,3293,1434,3019,3397,3396,3399,3398,1438,3015,2228,2229,2226,2227,2224,2225,2222,2223,2220,2221,1696,1697,1694,1695,1692,1693,1690,1691,1698,1699,1456,2088,2089,2083,2086,2087,2084,2085,1304,1305,1306,1307,1300,1301,1302,1303,1308,1309,3416,3417,3414,3415,3412,3413,3410,3411,3418,3419,1832,1833,1830,1831,1836,1837,1834,1835,1838,1839,2503,2901,2900,1713,1712,1711,1710,1717,1716,1715,1714,1719,1718,3665,3664,3667,3666,3661,3663,3662,3669,3668,1436,1437,3018,1435,1432,1433,1430,1431,3012,3013,3010,3011,3016,3017,3014,1439,3031,2882,2883,2880,2881,2886,2887,2884,2885,3690,3691,2888,2889,2110,2111,2112,2113,2622,2623,2620,2621,2626,2627,2624,2625,2628,2629,3289,3288,3287,3286,3285,3284,3283,3282,3281,3280,2873,2872,2871,2877,2876,2875,2874,2879,2878,2239,2238,2231,2230,2233,2232,2235,2234,2237,2236,1689,2095,2094,2097,2096,2091,2090,2093,2092,2099,2098,2488,2489,2486,2487,2484,2485,2482,2483,2480,2481,2244,2240,2241,2242,2243,1317,1316,1315,1314,1313,1312,1311,1310,1319,1318,3401,3400,3403,3402,3405,3404,3407,3406,3409,3408,1829,1828,1825,1824,1827,1826,1821,1820,1823,1822,3276,3277,3274,3275,3272,3273,3270,3271,3278,3279,3005,3004,3007,3006,3003,3002,3009,3008,1427,1426,1429,1428,3689,3688,2899,2898,2895,2894,2897,2896,2891,2890,2893,2892,1454,1455,2109,2108,1450,1451,1452,1453,2103,2102,2101,2100,2107,2106,2105,2104,2638,2635,2634,2637,2636,2631,2630,2633,2632,1298,1299,2499,2498,2491,2490,2493,2492,3030,2494,2497,2496,1449,1448,1322,1323,1320,1321,1326,1327,1324,1325,2495,2479,2478,3032,1850,1297,2608,2609,3265,3267,3266,3269,3268,2477,3684,2217,2216,2215,2214,2219,2218,1447,1446,1445,1444,1443,1442,1441,1440,2507,2506,2505,2504,1296,2502,2501,2500,3423,3422,3421,3420,3426,3425,3424',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 09:06:47.757 +14.30.0.0000',
       scan=['32,67,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_T14t14.J090647+143000_P16458v1_2018_04_30T17_28_34.285/working/VLASS1.1.sb34641006.eb34808190.58115.235410231486_split.ms'],
       imagename='VLASS1.1.ql.T14t14.J090647+143000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='1296,1297,1298,1299,1300,1301,1302,1303,1304,1305,1306,1307,1308,1309,1310,1311,1312,1313,1314,1315,1316,1317,1318,1319,1320,1321,1322,1323,1324,1325,1326,1327,1426,1427,1428,1429,1430,1431,1432,1433,1434,1435,1436,1437,1438,1439,1440,1441,1442,1443,1444,1445,1446,1447,1448,1449,1450,1451,1452,1453,1454,1455,1456,1457,1689,1690,1691,1692,1693,1694,1695,1696,1697,1698,1699,1700,1701,1702,1703,1704,1705,1706,1707,1708,1709,1710,1711,1712,1713,1714,1715,1716,1717,1718,1719,1820,1821,1822,1823,1824,1825,1826,1827,1828,1829,1830,1831,1832,1833,1834,1835,1836,1837,1838,1839,1840,1841,1842,1843,1844,1845,1846,1847,1848,1849,1850,2083,2084,2085,2086,2087,2088,2089,2090,2091,2092,2093,2094,2095,2096,2097,2098,2099,2100,2101,2102,2103,2104,2105,2106,2107,2108,2109,2110,2111,2112,2113,2214,2215,2216,2217,2218,2219,2220,2221,2222,2223,2224,2225,2226,2227,2228,2229,2230,2231,2232,2233,2234,2235,2236,2237,2238,2239,2240,2241,2242,2243,2244,2477,2478,2479,2480,2481,2482,2483,2484,2485,2486,2487,2488,2489,2490,2491,2492,2493,2494,2495,2496,2497,2498,2499,2500,2501,2502,2503,2504,2505,2506,2507,2608,2609,2610,2611,2612,2613,2614,2615,2616,2617,2618,2619,2620,2621,2622,2623,2624,2625,2626,2627,2628,2629,2630,2631,2632,2633,2634,2635,2636,2637,2638,2871,2872,2873,2874,2875,2876,2877,2878,2879,2880,2881,2882,2883,2884,2885,2886,2887,2888,2889,2890,2891,2892,2893,2894,2895,2896,2897,2898,2899,2900,2901,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,3027,3028,3029,3030,3031,3032,3265,3266,3267,3268,3269,3270,3271,3272,3273,3274,3275,3276,3277,3278,3279,3280,3281,3282,3283,3284,3285,3286,3287,3288,3289,3290,3291,3292,3293,3294,3295,3396,3397,3398,3399,3400,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,3661,3662,3663,3664,3665,3666,3667,3668,3669,3670,3671,3672,3673,3674,3675,3676,3677,3678,3679,3680,3681,3682,3683,3684,3685,3686,3687,3688,3689,3690,3691',
       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 09:06:47.757 +14.30.0.0000',
       scan=['32,67,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_T14t14.J090647+143000_P16458v1_2018_04_30T17_28_34.285/working/VLASS1.1.sb34641006.eb34808190.58115.235410231486_split.ms'],
       imagename='VLASS1.1.ql.T14t14.J090647+143000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.0005521264640443405, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='1296,1297,1298,1299,1300,1301,1302,1303,1304,1305,1306,1307,1308,1309,1310,1311,1312,1313,1314,1315,1316,1317,1318,1319,1320,1321,1322,1323,1324,1325,1326,1327,1426,1427,1428,1429,1430,1431,1432,1433,1434,1435,1436,1437,1438,1439,1440,1441,1442,1443,1444,1445,1446,1447,1448,1449,1450,1451,1452,1453,1454,1455,1456,1457,1689,1690,1691,1692,1693,1694,1695,1696,1697,1698,1699,1700,1701,1702,1703,1704,1705,1706,1707,1708,1709,1710,1711,1712,1713,1714,1715,1716,1717,1718,1719,1820,1821,1822,1823,1824,1825,1826,1827,1828,1829,1830,1831,1832,1833,1834,1835,1836,1837,1838,1839,1840,1841,1842,1843,1844,1845,1846,1847,1848,1849,1850,2083,2084,2085,2086,2087,2088,2089,2090,2091,2092,2093,2094,2095,2096,2097,2098,2099,2100,2101,2102,2103,2104,2105,2106,2107,2108,2109,2110,2111,2112,2113,2214,2215,2216,2217,2218,2219,2220,2221,2222,2223,2224,2225,2226,2227,2228,2229,2230,2231,2232,2233,2234,2235,2236,2237,2238,2239,2240,2241,2242,2243,2244,2477,2478,2479,2480,2481,2482,2483,2484,2485,2486,2487,2488,2489,2490,2491,2492,2493,2494,2495,2496,2497,2498,2499,2500,2501,2502,2503,2504,2505,2506,2507,2608,2609,2610,2611,2612,2613,2614,2615,2616,2617,2618,2619,2620,2621,2622,2623,2624,2625,2626,2627,2628,2629,2630,2631,2632,2633,2634,2635,2636,2637,2638,2871,2872,2873,2874,2875,2876,2877,2878,2879,2880,2881,2882,2883,2884,2885,2886,2887,2888,2889,2890,2891,2892,2893,2894,2895,2896,2897,2898,2899,2900,2901,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,3027,3028,3029,3030,3031,3032,3265,3266,3267,3268,3269,3270,3271,3272,3273,3274,3275,3276,3277,3278,3279,3280,3281,3282,3283,3284,3285,3286,3287,3288,3289,3290,3291,3292,3293,3294,3295,3396,3397,3398,3399,3400,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,3661,3662,3663,3664,3665,3666,3667,3668,3669,3670,3671,3672,3673,3674,3675,3676,3677,3678,3679,3680,3681,3682,3683,3684,3685,3686,3687,3688,3689,3690,3691',
       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.T14t14.J090647+143000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T14t14.J090647+143000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T14t14.J090647+143000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#