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

# hif_editimlist(parameter_file='parameter.12513.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_T19t18.J135727+333000_P12513v1_2018_02_01T00_19_11.324/working/VLASS1.1.sb34517306.eb34782681.58097.74992206019_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T19t18.J135727+333000_P12513v1_2018_02_01T00_19_11.324/rawdata/VLASS1.1.sb34517306.eb34782681.58097.74992206019.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='3382,1948,1949,1942,1943,1940,1941,1946,1947,1944,1945,3386,3387,3384,3385,2547,2546,3380,2544,2981,2548,3388,3389,2381,2380,2383,2382,2385,2384,2387,2386,2389,2388,3580,3581,2769,2768,1339,3585,3586,3587,3588,1537,1536,1535,1534,1533,1532,1531,1530,2165,2164,2166,2161,2160,1539,1538,3809,3808,3805,3804,3807,3806,3801,3800,3803,3802,3559,3558,1759,1758,1757,1756,1755,1754,1753,1752,1751,1750,2377,3799,3798,3795,3794,3797,3796,3791,3790,3793,3792,2154,2155,2156,2157,2150,2151,2152,2153,2158,2159,1955,1954,1957,1956,1951,1950,1953,1952,1959,1958,3562,3563,3560,3561,3566,3567,3564,3565,3568,3569,2798,2799,3153,3157,3156,3155,3154,3159,3158,3391,3390,3393,3392,3395,3394,3397,3396,3399,3398,2574,2572,2573,2570,2571,2378,2379,2374,2375,2376,2795,2370,2371,2372,2373,1348,1349,1340,1341,1342,1343,1344,1345,1346,1347,1960,1961,1962,2978,2979,2972,2973,2970,2971,2976,2977,2974,3383,2545,3381,1768,1769,1762,1763,1760,1761,1766,2549,1764,1765,2980,2954,2955,2956,1546,1547,1544,1545,1542,1543,1540,1541,2950,1548,1549,2952,2953,3401,3400,3403,3402,3405,3404,3407,3406,2951,2147,2146,2145,2144,2143,2142,2141,2140,2792,2793,2790,2791,2796,2797,2149,2148,3582,3583,3584,3579,3578,3575,3574,3577,3576,3571,3570,3573,3572,3168,3169,3166,3167,3164,3165,3162,3163,3160,3161,2569,2568,2561,2560,2563,2562,2565,2564,2567,2566,2369,2368,2367,2366,2365,2364,2363,2362,2361,1359,1358,1353,1352,1351,1350,1357,1356,1355,1354,2794,3184,3180,3181,3182,3183,2163,2162,2965,2964,2967,2966,2961,2960,2963,2962,2969,2968,1775,1774,1777,1776,1771,1770,1773,1772,1779,1778,1551,1550,1552,2550,2551,2552,2553,2554,2555,2556,2557,2558,2559,2390,2391,2778,2779,2770,2771,2772,2773,2774,2775,2776,2777,2785,2784,2787,2786,2781,2780,2783,2782,3172,2789,2788,1524,1525,1526,1527,1522,1523,1528,1529,3810,3811,1933,1932,1937,1936,1935,1934,1939,1938,2975,3179,3178,2958,2959,3171,3170,3173,2957,3175,3174,3177,3176,3379,3378,3377,1749,2136,2137,2138,2139,3782,3783,3781,3786,3787,3784,3785,3788,3789,1767,1366,1367,1364,1365,1362,1363,1360,1361,1368,1369',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 13:57:27.752 +33.30.0.0000',
       scan=['33,59,36,65,40,43,46,49,52,62,56,27,30'],
       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_T19t18.J135727+333000_P12513v1_2018_02_01T00_19_11.324/working/VLASS1.1.sb34517306.eb34782681.58097.74992206019_split.ms'],
       imagename='VLASS1.1.ql.T19t18.J135727+333000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='1339,1340,1341,1342,1343,1344,1345,1346,1347,1348,1349,1350,1351,1352,1353,1354,1355,1356,1357,1358,1359,1360,1361,1362,1363,1364,1365,1366,1367,1368,1369,1522,1523,1524,1525,1526,1527,1528,1529,1530,1531,1532,1533,1534,1535,1536,1537,1538,1539,1540,1541,1542,1543,1544,1545,1546,1547,1548,1549,1550,1551,1552,1749,1750,1751,1752,1753,1754,1755,1756,1757,1758,1759,1760,1761,1762,1763,1764,1765,1766,1767,1768,1769,1770,1771,1772,1773,1774,1775,1776,1777,1778,1779,1932,1933,1934,1935,1936,1937,1938,1939,1940,1941,1942,1943,1944,1945,1946,1947,1948,1949,1950,1951,1952,1953,1954,1955,1956,1957,1958,1959,1960,1961,1962,2136,2137,2138,2139,2140,2141,2142,2143,2144,2145,2146,2147,2148,2149,2150,2151,2152,2153,2154,2155,2156,2157,2158,2159,2160,2161,2162,2163,2164,2165,2166,2361,2362,2363,2364,2365,2366,2367,2368,2369,2370,2371,2372,2373,2374,2375,2376,2377,2378,2379,2380,2381,2382,2383,2384,2385,2386,2387,2388,2389,2390,2391,2544,2545,2546,2547,2548,2549,2550,2551,2552,2553,2554,2555,2556,2557,2558,2559,2560,2561,2562,2563,2564,2565,2566,2567,2568,2569,2570,2571,2572,2573,2574,2768,2769,2770,2771,2772,2773,2774,2775,2776,2777,2778,2779,2780,2781,2782,2783,2784,2785,2786,2787,2788,2789,2790,2791,2792,2793,2794,2795,2796,2797,2798,2799,2950,2951,2952,2953,2954,2955,2956,2957,2958,2959,2960,2961,2962,2963,2964,2965,2966,2967,2968,2969,2970,2971,2972,2973,2974,2975,2976,2977,2978,2979,2980,2981,3153,3154,3155,3156,3157,3158,3159,3160,3161,3162,3163,3164,3165,3166,3167,3168,3169,3170,3171,3172,3173,3174,3175,3176,3177,3178,3179,3180,3181,3182,3183,3184,3377,3378,3379,3380,3381,3382,3383,3384,3385,3386,3387,3388,3389,3390,3391,3392,3393,3394,3395,3396,3397,3398,3399,3400,3401,3402,3403,3404,3405,3406,3407,3558,3559,3560,3561,3562,3563,3564,3565,3566,3567,3568,3569,3570,3571,3572,3573,3574,3575,3576,3577,3578,3579,3580,3581,3582,3583,3584,3585,3586,3587,3588,3781,3782,3783,3784,3785,3786,3787,3788,3789,3790,3791,3792,3793,3794,3795,3796,3797,3798,3799,3800,3801,3802,3803,3804,3805,3806,3807,3808,3809,3810,3811',
       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 13:57:27.752 +33.30.0.0000',
       scan=['33,59,36,65,40,43,46,49,52,62,56,27,30'],
       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_T19t18.J135727+333000_P12513v1_2018_02_01T00_19_11.324/working/VLASS1.1.sb34517306.eb34782681.58097.74992206019_split.ms'],
       imagename='VLASS1.1.ql.T19t18.J135727+333000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00046007942630831163, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='1339,1340,1341,1342,1343,1344,1345,1346,1347,1348,1349,1350,1351,1352,1353,1354,1355,1356,1357,1358,1359,1360,1361,1362,1363,1364,1365,1366,1367,1368,1369,1522,1523,1524,1525,1526,1527,1528,1529,1530,1531,1532,1533,1534,1535,1536,1537,1538,1539,1540,1541,1542,1543,1544,1545,1546,1547,1548,1549,1550,1551,1552,1749,1750,1751,1752,1753,1754,1755,1756,1757,1758,1759,1760,1761,1762,1763,1764,1765,1766,1767,1768,1769,1770,1771,1772,1773,1774,1775,1776,1777,1778,1779,1932,1933,1934,1935,1936,1937,1938,1939,1940,1941,1942,1943,1944,1945,1946,1947,1948,1949,1950,1951,1952,1953,1954,1955,1956,1957,1958,1959,1960,1961,1962,2136,2137,2138,2139,2140,2141,2142,2143,2144,2145,2146,2147,2148,2149,2150,2151,2152,2153,2154,2155,2156,2157,2158,2159,2160,2161,2162,2163,2164,2165,2166,2361,2362,2363,2364,2365,2366,2367,2368,2369,2370,2371,2372,2373,2374,2375,2376,2377,2378,2379,2380,2381,2382,2383,2384,2385,2386,2387,2388,2389,2390,2391,2544,2545,2546,2547,2548,2549,2550,2551,2552,2553,2554,2555,2556,2557,2558,2559,2560,2561,2562,2563,2564,2565,2566,2567,2568,2569,2570,2571,2572,2573,2574,2768,2769,2770,2771,2772,2773,2774,2775,2776,2777,2778,2779,2780,2781,2782,2783,2784,2785,2786,2787,2788,2789,2790,2791,2792,2793,2794,2795,2796,2797,2798,2799,2950,2951,2952,2953,2954,2955,2956,2957,2958,2959,2960,2961,2962,2963,2964,2965,2966,2967,2968,2969,2970,2971,2972,2973,2974,2975,2976,2977,2978,2979,2980,2981,3153,3154,3155,3156,3157,3158,3159,3160,3161,3162,3163,3164,3165,3166,3167,3168,3169,3170,3171,3172,3173,3174,3175,3176,3177,3178,3179,3180,3181,3182,3183,3184,3377,3378,3379,3380,3381,3382,3383,3384,3385,3386,3387,3388,3389,3390,3391,3392,3393,3394,3395,3396,3397,3398,3399,3400,3401,3402,3403,3404,3405,3406,3407,3558,3559,3560,3561,3562,3563,3564,3565,3566,3567,3568,3569,3570,3571,3572,3573,3574,3575,3576,3577,3578,3579,3580,3581,3582,3583,3584,3585,3586,3587,3588,3781,3782,3783,3784,3785,3786,3787,3788,3789,3790,3791,3792,3793,3794,3795,3796,3797,3798,3799,3800,3801,3802,3803,3804,3805,3806,3807,3808,3809,3810,3811',
       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.T19t18.J135727+333000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T19t18.J135727+333000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T19t18.J135727+333000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#