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

# hif_editimlist(parameter_file='parameter.5419.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_T26t07.J080424+633000_P5419v1_2017_09_30T19_18_22.718/working/VLASS1.1.sb34398871.eb34498197.58021.421088125004_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T26t07.J080424+633000_P5419v1_2017_09_30T19_18_22.718/rawdata/VLASS1.1.sb34398871.eb34498197.58021.421088125004.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,2659,2658,1780,1781,1782,1783,1784,1785,1786,1787,3049,3048,3057,3415,3412,3413,3410,3411,3041,3040,3043,3042,3047,3046,1775,1774,1777,1776,1771,1770,1078,1772,1076,1077,1074,1075,1072,1073,1070,1071,2154,2653,2150,2151,2152,2153,3414,2670,2651,2466,2467,2660,3019,2490,1792,1791,1790,2468,2469,3055,2656,3401,1403,3403,3402,3405,2655,3407,1402,3409,3408,2488,1401,1069,1068,2649,1400,3038,3039,3400,1061,1060,1063,1062,1065,1064,1067,1066,1418,1419,2145,2144,1388,1391,2141,2140,1410,1411,1412,1413,1414,1415,1416,1417,2477,2476,2475,2474,2473,2472,2471,2470,3030,2694,2479,2478,2104,3032,2486,3033,2684,2147,2149,3036,3029,2146,3027,2693,3025,1773,3023,3022,3021,3020,2666,2667,1409,2665,2662,2663,1399,1398,1397,1396,1395,1394,1393,1392,1405,1390,3037,3051,2657,2674,1384,2692,1385,3034,2686,2652,1386,3404,1389,3391,3390,3018,1387,1432,1433,1430,1431,3399,2107,2489,2106,2671,1406,2673,2672,1382,2677,2487,2679,2678,1383,2484,3031,2143,2485,1779,3071,2687,1778,2695,2480,2696,2133,3406,2481,2142,1421,1420,1423,1422,1425,1424,1427,1426,1429,1428,1404,3053,2118,2119,2688,2689,2130,2110,2111,2112,2113,2114,2115,2116,2117,3035,3024,3393,2650,3070,2483,2482,2685,3395,2661,3394,2664,3397,2680,2691,2126,2681,2683,2109,2108,1748,1749,2699,2698,1744,1745,1746,1747,1740,1741,1742,1743,2668,2669,3398,3063,3062,3061,3060,3067,3066,3065,3064,2654,3069,3068,2136,2137,2134,2135,1054,1055,1056,1057,1058,1059,2138,2139,3392,1759,1758,1757,1756,1755,1754,1753,1752,1751,1750,2682,2103,2491,2690,2701,2700,2148,3058,3059,3056,2132,3054,2702,3052,3028,3050,2102,2121,2120,2123,2122,2125,2124,2697,2105,2129,2131,3396,1768,1769,3026,1762,1763,1760,1761,2101,1764,1765',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 08:04:24.323 +63.30.0.0000',
       scan=['65,34,59,37,40,43,46,49,53,56,27,62,31'],
       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_T26t07.J080424+633000_P5419v1_2017_09_30T19_18_22.718/working/VLASS1.1.sb34398871.eb34498197.58021.421088125004_split.ms'],
       imagename='VLASS1.1.ql.T26t07.J080424+633000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='1054,1055,1056,1057,1058,1059,1060,1061,1062,1063,1064,1065,1066,1067,1068,1069,1070,1071,1072,1073,1074,1075,1076,1077,1078,1382,1383,1384,1385,1386,1387,1388,1389,1390,1391,1392,1393,1394,1395,1396,1397,1398,1399,1400,1401,1402,1403,1404,1405,1406,1409,1410,1411,1412,1413,1414,1415,1416,1417,1418,1419,1420,1421,1422,1423,1424,1425,1426,1427,1428,1429,1430,1431,1432,1433,1740,1741,1742,1743,1744,1745,1746,1747,1748,1749,1750,1751,1752,1753,1754,1755,1756,1757,1758,1759,1760,1761,1762,1763,1764,1765,1768,1769,1770,1771,1772,1773,1774,1775,1776,1777,1778,1779,1780,1781,1782,1783,1784,1785,1786,1787,1788,1789,1790,1791,1792,2101,2102,2103,2104,2105,2106,2107,2108,2109,2110,2111,2112,2113,2114,2115,2116,2117,2118,2119,2120,2121,2122,2123,2124,2125,2126,2129,2130,2131,2132,2133,2134,2135,2136,2137,2138,2139,2140,2141,2142,2143,2144,2145,2146,2147,2148,2149,2150,2151,2152,2153,2154,2466,2467,2468,2469,2470,2471,2472,2473,2474,2475,2476,2477,2478,2479,2480,2481,2482,2483,2484,2485,2486,2487,2488,2489,2490,2491,2649,2650,2651,2652,2653,2654,2655,2656,2657,2658,2659,2660,2661,2662,2663,2664,2665,2666,2667,2668,2669,2670,2671,2672,2673,2674,2677,2678,2679,2680,2681,2682,2683,2684,2685,2686,2687,2688,2689,2690,2691,2692,2693,2694,2695,2696,2697,2698,2699,2700,2701,2702,3018,3019,3020,3021,3022,3023,3024,3025,3026,3027,3028,3029,3030,3031,3032,3033,3034,3035,3036,3037,3038,3039,3040,3041,3042,3043,3046,3047,3048,3049,3050,3051,3052,3053,3054,3055,3056,3057,3058,3059,3060,3061,3062,3063,3064,3065,3066,3067,3068,3069,3070,3071,3390,3391,3392,3393,3394,3395,3396,3397,3398,3399,3400,3401,3402,3403,3404,3405,3406,3407,3408,3409,3410,3411,3412,3413,3414,3415',
       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 08:04:24.323 +63.30.0.0000',
       scan=['65,34,59,37,40,43,46,49,53,56,27,62,31'],
       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_T26t07.J080424+633000_P5419v1_2017_09_30T19_18_22.718/working/VLASS1.1.sb34398871.eb34498197.58021.421088125004_split.ms'],
       imagename='VLASS1.1.ql.T26t07.J080424+633000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.0005800063064569258, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='1054,1055,1056,1057,1058,1059,1060,1061,1062,1063,1064,1065,1066,1067,1068,1069,1070,1071,1072,1073,1074,1075,1076,1077,1078,1382,1383,1384,1385,1386,1387,1388,1389,1390,1391,1392,1393,1394,1395,1396,1397,1398,1399,1400,1401,1402,1403,1404,1405,1406,1409,1410,1411,1412,1413,1414,1415,1416,1417,1418,1419,1420,1421,1422,1423,1424,1425,1426,1427,1428,1429,1430,1431,1432,1433,1740,1741,1742,1743,1744,1745,1746,1747,1748,1749,1750,1751,1752,1753,1754,1755,1756,1757,1758,1759,1760,1761,1762,1763,1764,1765,1768,1769,1770,1771,1772,1773,1774,1775,1776,1777,1778,1779,1780,1781,1782,1783,1784,1785,1786,1787,1788,1789,1790,1791,1792,2101,2102,2103,2104,2105,2106,2107,2108,2109,2110,2111,2112,2113,2114,2115,2116,2117,2118,2119,2120,2121,2122,2123,2124,2125,2126,2129,2130,2131,2132,2133,2134,2135,2136,2137,2138,2139,2140,2141,2142,2143,2144,2145,2146,2147,2148,2149,2150,2151,2152,2153,2154,2466,2467,2468,2469,2470,2471,2472,2473,2474,2475,2476,2477,2478,2479,2480,2481,2482,2483,2484,2485,2486,2487,2488,2489,2490,2491,2649,2650,2651,2652,2653,2654,2655,2656,2657,2658,2659,2660,2661,2662,2663,2664,2665,2666,2667,2668,2669,2670,2671,2672,2673,2674,2677,2678,2679,2680,2681,2682,2683,2684,2685,2686,2687,2688,2689,2690,2691,2692,2693,2694,2695,2696,2697,2698,2699,2700,2701,2702,3018,3019,3020,3021,3022,3023,3024,3025,3026,3027,3028,3029,3030,3031,3032,3033,3034,3035,3036,3037,3038,3039,3040,3041,3042,3043,3046,3047,3048,3049,3050,3051,3052,3053,3054,3055,3056,3057,3058,3059,3060,3061,3062,3063,3064,3065,3066,3067,3068,3069,3070,3071,3390,3391,3392,3393,3394,3395,3396,3397,3398,3399,3400,3401,3402,3403,3404,3405,3406,3407,3408,3409,3410,3411,3412,3413,3414,3415',
       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.T26t07.J080424+633000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T26t07.J080424+633000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T26t07.J080424+633000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#