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.sb34697045.eb34742816.58084.0767340625.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_T13t02.J005813+093000_P18524v1_2018_03_17T19_48_15.410/working/VLASS1.1.sb34697045.eb34742816.58084.0767340625_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T13t02.J005813+093000_P18524v1_2018_03_17T19_48_15.410/rawdata/VLASS1.1.sb34697045.eb34742816.58084.0767340625.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='1718,3714,3712,2913,2914,2915,2916,2917,2918,2919,3718,1700,1701,1702,1703,1704,1705,1706,2319,1708,1709,2314,2315,2312,2313,2310,2311,2708,3140,3141,3142,3143,2541,2540,1519,1518,3128,1515,3129,1514,1517,1516,1510,1536,1535,1534,1533,1532,1531,1530,1512,3121,3126,2528,3298,3299,1908,1909,3296,3297,3519,2128,3513,3517,3516,3515,3514,2709,2943,2942,2941,2940,1696,1697,1694,1695,1692,1693,1690,1691,1698,1699,3708,3709,3702,3703,3700,3701,3706,3707,3704,3705,2121,2120,2123,2122,2125,2124,3119,3118,3117,3116,3115,3114,3113,1288,2538,2539,2532,2533,2530,2531,2536,2537,2534,2535,1304,1305,1306,1307,1300,1301,1302,1303,2718,2719,1308,1309,3518,2318,1707,2316,2317,3308,3309,3306,3307,3304,3305,3302,3303,3300,3301,1713,1712,1711,1710,1717,1716,1715,1714,2301,2300,1719,2302,2305,2304,2307,2306,1513,2308,2714,3698,3699,3696,3697,2118,2119,1506,1507,2110,2111,2112,1509,2114,2115,2116,2117,1287,1919,1918,3524,3525,3522,3523,3520,3521,1911,1910,1913,1912,1915,1508,1917,1916,2113,2127,2126,2129,2699,2698,2697,2696,2695,2694,1289,3526,3527,1689,3715,2712,3717,3716,3711,3710,3713,2713,3719,2710,2711,2716,2717,2525,2524,2527,2526,2521,2520,2523,2522,3122,3123,3120,2715,2529,3127,3124,3125,3528,3529,1317,1316,1315,1314,1313,1312,1311,1310,2705,2704,2707,2706,2701,2700,2703,2702,1920,1921,1922,1923,1924,1925,1926,1927,1928,1929,2938,2939,2936,2937,2934,2935,2932,2933,2930,2931,3319,3318,3311,3310,3313,3312,3315,3314,3317,3316,2299,2298,2293,2292,2297,2296,2295,2294,3724,3725,3726,3720,3721,3722,3723,2109,2518,2519,2514,2515,2516,2517,1511,2511,2512,2513,3531,3530,3533,3532,3535,3534,3537,3536,3539,3538,3540,3541,3542,3543,3320,3321,3322,3323,3324,3325,3326,2322,2321,2320,3135,3134,3137,3136,3131,3130,3133,3132,3139,3138,1914,1524,1525,1526,1527,1520,1521,1522,1523,1528,1529,1933,1932,1931,1930,1937,1936,1935,1934,1938,2929,2928,2921,2920,2923,2922,2925,2924,2927,2926,2136,2137,2134,2135,2132,2133,2130,2131,2138,2139,1298,1299,1292,1293,1290,1291,1296,1297,1294,1295,2309,2724,2723,2722,2721,2720,2303',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 00:58:13.462 +09.30.0.0000',
       scan=['32,64,67,35,38,42,45,61,48,51,54,58,29'],
       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_T13t02.J005813+093000_P18524v1_2018_03_17T19_48_15.410/working/VLASS1.1.sb34697045.eb34742816.58084.0767340625_split.ms'],
       imagename='VLASS1.1.ql.T13t02.J005813+093000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='1287,1288,1289,1290,1291,1292,1293,1294,1295,1296,1297,1298,1299,1300,1301,1302,1303,1304,1305,1306,1307,1308,1309,1310,1311,1312,1313,1314,1315,1316,1317,1506,1507,1508,1509,1510,1511,1512,1513,1514,1515,1516,1517,1518,1519,1520,1521,1522,1523,1524,1525,1526,1527,1528,1529,1530,1531,1532,1533,1534,1535,1536,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,1908,1909,1910,1911,1912,1913,1914,1915,1916,1917,1918,1919,1920,1921,1922,1923,1924,1925,1926,1927,1928,1929,1930,1931,1932,1933,1934,1935,1936,1937,1938,2109,2110,2111,2112,2113,2114,2115,2116,2117,2118,2119,2120,2121,2122,2123,2124,2125,2126,2127,2128,2129,2130,2131,2132,2133,2134,2135,2136,2137,2138,2139,2292,2293,2294,2295,2296,2297,2298,2299,2300,2301,2302,2303,2304,2305,2306,2307,2308,2309,2310,2311,2312,2313,2314,2315,2316,2317,2318,2319,2320,2321,2322,2511,2512,2513,2514,2515,2516,2517,2518,2519,2520,2521,2522,2523,2524,2525,2526,2527,2528,2529,2530,2531,2532,2533,2534,2535,2536,2537,2538,2539,2540,2541,2694,2695,2696,2697,2698,2699,2700,2701,2702,2703,2704,2705,2706,2707,2708,2709,2710,2711,2712,2713,2714,2715,2716,2717,2718,2719,2720,2721,2722,2723,2724,2913,2914,2915,2916,2917,2918,2919,2920,2921,2922,2923,2924,2925,2926,2927,2928,2929,2930,2931,2932,2933,2934,2935,2936,2937,2938,2939,2940,2941,2942,2943,3113,3114,3115,3116,3117,3118,3119,3120,3121,3122,3123,3124,3125,3126,3127,3128,3129,3130,3131,3132,3133,3134,3135,3136,3137,3138,3139,3140,3141,3142,3143,3296,3297,3298,3299,3300,3301,3302,3303,3304,3305,3306,3307,3308,3309,3310,3311,3312,3313,3314,3315,3316,3317,3318,3319,3320,3321,3322,3323,3324,3325,3326,3513,3514,3515,3516,3517,3518,3519,3520,3521,3522,3523,3524,3525,3526,3527,3528,3529,3530,3531,3532,3533,3534,3535,3536,3537,3538,3539,3540,3541,3542,3543,3696,3697,3698,3699,3700,3701,3702,3703,3704,3705,3706,3707,3708,3709,3710,3711,3712,3713,3714,3715,3716,3717,3718,3719,3720,3721,3722,3723,3724,3725,3726',
       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 00:58:13.462 +09.30.0.0000',
       scan=['32,64,67,35,38,42,45,61,48,51,54,58,29'],
       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_T13t02.J005813+093000_P18524v1_2018_03_17T19_48_15.410/working/VLASS1.1.sb34697045.eb34742816.58084.0767340625_split.ms'],
       imagename='VLASS1.1.ql.T13t02.J005813+093000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.0005925693632681941, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='1287,1288,1289,1290,1291,1292,1293,1294,1295,1296,1297,1298,1299,1300,1301,1302,1303,1304,1305,1306,1307,1308,1309,1310,1311,1312,1313,1314,1315,1316,1317,1506,1507,1508,1509,1510,1511,1512,1513,1514,1515,1516,1517,1518,1519,1520,1521,1522,1523,1524,1525,1526,1527,1528,1529,1530,1531,1532,1533,1534,1535,1536,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,1908,1909,1910,1911,1912,1913,1914,1915,1916,1917,1918,1919,1920,1921,1922,1923,1924,1925,1926,1927,1928,1929,1930,1931,1932,1933,1934,1935,1936,1937,1938,2109,2110,2111,2112,2113,2114,2115,2116,2117,2118,2119,2120,2121,2122,2123,2124,2125,2126,2127,2128,2129,2130,2131,2132,2133,2134,2135,2136,2137,2138,2139,2292,2293,2294,2295,2296,2297,2298,2299,2300,2301,2302,2303,2304,2305,2306,2307,2308,2309,2310,2311,2312,2313,2314,2315,2316,2317,2318,2319,2320,2321,2322,2511,2512,2513,2514,2515,2516,2517,2518,2519,2520,2521,2522,2523,2524,2525,2526,2527,2528,2529,2530,2531,2532,2533,2534,2535,2536,2537,2538,2539,2540,2541,2694,2695,2696,2697,2698,2699,2700,2701,2702,2703,2704,2705,2706,2707,2708,2709,2710,2711,2712,2713,2714,2715,2716,2717,2718,2719,2720,2721,2722,2723,2724,2913,2914,2915,2916,2917,2918,2919,2920,2921,2922,2923,2924,2925,2926,2927,2928,2929,2930,2931,2932,2933,2934,2935,2936,2937,2938,2939,2940,2941,2942,2943,3113,3114,3115,3116,3117,3118,3119,3120,3121,3122,3123,3124,3125,3126,3127,3128,3129,3130,3131,3132,3133,3134,3135,3136,3137,3138,3139,3140,3141,3142,3143,3296,3297,3298,3299,3300,3301,3302,3303,3304,3305,3306,3307,3308,3309,3310,3311,3312,3313,3314,3315,3316,3317,3318,3319,3320,3321,3322,3323,3324,3325,3326,3513,3514,3515,3516,3517,3518,3519,3520,3521,3522,3523,3524,3525,3526,3527,3528,3529,3530,3531,3532,3533,3534,3535,3536,3537,3538,3539,3540,3541,3542,3543,3696,3697,3698,3699,3700,3701,3702,3703,3704,3705,3706,3707,3708,3709,3710,3711,3712,3713,3714,3715,3716,3717,3718,3719,3720,3721,3722,3723,3724,3725,3726',
       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.T13t02.J005813+093000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T13t02.J005813+093000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T13t02.J005813+093000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#