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

# hif_editimlist(parameter_file='parameter.12837.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_T31t03.J160254+813000_P12837v1_2017_10_17T15_51_22.121/working/VLASS1.1.sb34535845.eb34554723.58025.888074629635_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T31t03.J160254+813000_P12837v1_2017_10_17T15_51_22.121/rawdata/VLASS1.1.sb34535845.eb34554723.58025.888074629635.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='1218,2859,2858,2851,2850,2853,2852,2855,2854,2857,2856,2268,2269,1019,2262,2263,2260,2261,2266,2267,2264,2265,2103,108,109,2102,102,103,100,101,106,107,104,105,2517,2510,2106,2105,647,649,648,2513,99,98,91,90,93,92,95,94,97,96,2860,2861,2862,2863,2864,2865,2866,2867,2868,2869,1199,1198,1197,1759,1758,1757,1756,1755,1754,1753,1752,1751,1750,2088,2089,2086,2087,2084,2085,1579,1578,2532,2533,1571,1570,2536,2537,2534,2535,494,495,490,491,492,493,88,89,82,83,86,87,84,85,2280,2281,2282,2275,2274,2277,2276,2271,2270,2273,2272,2279,2278,1227,1226,1225,1224,1223,1222,1221,1220,656,657,654,655,652,653,650,651,2110,2111,2112,2113,2114,658,659,2873,2872,2871,2870,2874,1573,1572,2530,1043,1042,1041,2531,1047,1046,1045,1044,1577,1049,1048,1576,1575,1762,1763,1760,1761,669,668,667,666,665,664,663,662,661,660,2095,2094,2097,2096,2091,2090,2093,2092,2099,2098,2525,2524,2527,2526,2521,2520,2523,2522,2529,2528,2516,1588,1589,1582,1583,1580,1581,1586,1587,1584,1585,1038,1039,1032,1033,1030,1031,1036,1037,1034,1035,1212,1213,1210,1211,469,468,1214,1215,465,464,467,466,2109,2108,2518,2519,2514,2515,2101,2100,2107,2511,2512,2104,1574,2848,2849,2846,2847,2844,2845,2843,674,675,676,677,670,671,672,673,2507,1559,2259,2258,2257,2256,2255,2254,2253,2252,111,110,112,1739,1738,1733,1732,1735,1734,1737,1736,1040,1025,1024,1027,1026,1021,1020,1023,1022,1029,1028,1748,1749,1744,1745,1746,1747,1740,1741,1742,1743,1209,1208,1205,1204,1207,1206,1201,1200,1203,1202,1216,1217,1568,1569,2509,2508,1560,1561,1562,1563,1564,1565,1566,1567,1219,489,488,487,486,485,484,483,482,481,480,472,473,470,471,476,477,474,475,478,479',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 16:02:54.793 +81.30.0.0000',
       scan=['34,38,9,12,15,18,22,41,25,28,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_T31t03.J160254+813000_P12837v1_2017_10_17T15_51_22.121/working/VLASS1.1.sb34535845.eb34554723.58025.888074629635_split.ms'],
       imagename='VLASS1.1.ql.T31t03.J160254+813000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,464,465,466,467,468,469,470,471,472,473,474,475,476,477,478,479,480,481,482,483,484,485,486,487,488,489,490,491,492,493,494,495,647,648,649,650,651,652,653,654,655,656,657,658,659,660,661,662,663,664,665,666,667,668,669,670,671,672,673,674,675,676,677,1019,1020,1021,1022,1023,1024,1025,1026,1027,1028,1029,1030,1031,1032,1033,1034,1035,1036,1037,1038,1039,1040,1041,1042,1043,1044,1045,1046,1047,1048,1049,1197,1198,1199,1200,1201,1202,1203,1204,1205,1206,1207,1208,1209,1210,1211,1212,1213,1214,1215,1216,1217,1218,1219,1220,1221,1222,1223,1224,1225,1226,1227,1559,1560,1561,1562,1563,1564,1565,1566,1567,1568,1569,1570,1571,1572,1573,1574,1575,1576,1577,1578,1579,1580,1581,1582,1583,1584,1585,1586,1587,1588,1589,1732,1733,1734,1735,1736,1737,1738,1739,1740,1741,1742,1743,1744,1745,1746,1747,1748,1749,1750,1751,1752,1753,1754,1755,1756,1757,1758,1759,1760,1761,1762,1763,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,2114,2252,2253,2254,2255,2256,2257,2258,2259,2260,2261,2262,2263,2264,2265,2266,2267,2268,2269,2270,2271,2272,2273,2274,2275,2276,2277,2278,2279,2280,2281,2282,2507,2508,2509,2510,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,2843,2844,2845,2846,2847,2848,2849,2850,2851,2852,2853,2854,2855,2856,2857,2858,2859,2860,2861,2862,2863,2864,2865,2866,2867,2868,2869,2870,2871,2872,2873,2874',
       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 16:02:54.793 +81.30.0.0000',
       scan=['34,38,9,12,15,18,22,41,25,28,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_T31t03.J160254+813000_P12837v1_2017_10_17T15_51_22.121/working/VLASS1.1.sb34535845.eb34554723.58025.888074629635_split.ms'],
       imagename='VLASS1.1.ql.T31t03.J160254+813000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00058189766334435259, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,464,465,466,467,468,469,470,471,472,473,474,475,476,477,478,479,480,481,482,483,484,485,486,487,488,489,490,491,492,493,494,495,647,648,649,650,651,652,653,654,655,656,657,658,659,660,661,662,663,664,665,666,667,668,669,670,671,672,673,674,675,676,677,1019,1020,1021,1022,1023,1024,1025,1026,1027,1028,1029,1030,1031,1032,1033,1034,1035,1036,1037,1038,1039,1040,1041,1042,1043,1044,1045,1046,1047,1048,1049,1197,1198,1199,1200,1201,1202,1203,1204,1205,1206,1207,1208,1209,1210,1211,1212,1213,1214,1215,1216,1217,1218,1219,1220,1221,1222,1223,1224,1225,1226,1227,1559,1560,1561,1562,1563,1564,1565,1566,1567,1568,1569,1570,1571,1572,1573,1574,1575,1576,1577,1578,1579,1580,1581,1582,1583,1584,1585,1586,1587,1588,1589,1732,1733,1734,1735,1736,1737,1738,1739,1740,1741,1742,1743,1744,1745,1746,1747,1748,1749,1750,1751,1752,1753,1754,1755,1756,1757,1758,1759,1760,1761,1762,1763,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,2114,2252,2253,2254,2255,2256,2257,2258,2259,2260,2261,2262,2263,2264,2265,2266,2267,2268,2269,2270,2271,2272,2273,2274,2275,2276,2277,2278,2279,2280,2281,2282,2507,2508,2509,2510,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,2843,2844,2845,2846,2847,2848,2849,2850,2851,2852,2853,2854,2855,2856,2857,2858,2859,2860,2861,2862,2863,2864,2865,2866,2867,2868,2869,2870,2871,2872,2873,2874',
       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.T31t03.J160254+813000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T31t03.J160254+813000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T31t03.J160254+813000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#