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.sb34905760.eb35113377.58166.981752928245.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='VLASS1.1.sb34905760.eb35113377.58166.981752928245_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T03t06.J035459-313000_P23819v1_2018_03_26T13_38_16.511/rawdata/VLASS1.1.sb34905760.eb35113377.58166.981752928245.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='818,819,810,811,812,813,814,815,816,817,1700,1701,1702,1703,1669,1668,1667,1666,1665,449,448,443,442,441,440,447,446,445,444,902,900,38,33,32,31,30,37,36,35,34,1241,1240,1243,1242,1245,1244,438,439,1249,1248,1696,1697,1694,1695,1692,1693,1690,1691,1757,1756,1755,1754,1753,1752,1698,1699,1270,1271,1272,1273,1274,1275,2121,2120,2123,2122,2125,2124,2127,2126,2129,2128,1304,1305,1306,1307,1300,1301,1302,1303,1308,1309,24,25,26,27,20,21,22,23,28,29,407,406,405,404,403,402,401,400,409,408,1758,1751,829,828,825,824,827,826,821,820,823,822,393,392,391,390,397,396,395,394,399,398,1247,1750,458,459,450,451,452,453,454,455,456,457,2118,2119,2110,2111,2112,2113,2114,2115,2116,2117,869,868,13,12,15,14,17,16,19,18,865,864,867,866,883,882,881,880,887,886,885,884,889,888,1681,1680,1683,1682,1685,1684,1687,1686,1689,1688,1269,1268,1263,1262,1261,1260,1267,1266,1265,1264,2095,2097,2096,2099,2098,1317,1316,1315,1314,1313,1312,1311,1310,1319,1318,414,415,416,417,410,411,412,413,418,419,1237,832,833,830,831,836,837,834,835,838,839,1238,1330,383,384,385,386,387,388,389,1728,1729,1726,1727,1724,1725,1722,1723,1720,1721,469,468,465,464,467,466,461,39,463,462,901,2109,2108,2103,2102,2101,2100,2107,2106,2105,2104,878,879,876,877,874,875,872,873,870,871,890,891,892,893,894,895,896,897,898,899,1246,2130,2131,809,1678,1679,1674,1675,1676,1677,1670,1671,1672,1673,50,1322,1323,1320,1321,1326,1327,1324,1325,1328,1329,1256,1257,1254,1255,1252,1253,1250,1251,421,420,1258,1259,847,846,845,844,843,842,841,840,1739,1738,1731,1730,1733,1732,1735,1734,1737,1736,1239,1748,1749,1744,1745,1746,1747,1740,1741,1742,1743,2132,2133,48,49,46,47,44,45,42,43,40,41,1298,1299,1292,1293,1296,1297,1294,1295,460,472,473,470,471,476,474,475',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 03:54:59.940 -31.30.0.0000',
       scan=['34,37,9,12,15,18,21,24,28,31'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34905760.eb35113377.58166.981752928245_split.ms'],
       imagename='VLASS1.1.ql.T03t06.J035459-313000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,383,384,385,386,387,388,389,390,391,392,393,394,395,396,397,398,399,400,401,402,403,404,405,406,407,408,409,410,411,412,413,414,415,416,417,418,419,420,421,438,439,440,441,442,443,444,445,446,447,448,449,450,451,452,453,454,455,456,457,458,459,460,461,462,463,464,465,466,467,468,469,470,471,472,473,474,475,476,809,810,811,812,813,814,815,816,817,818,819,820,821,822,823,824,825,826,827,828,829,830,831,832,833,834,835,836,837,838,839,840,841,842,843,844,845,846,847,864,865,866,867,868,869,870,871,872,873,874,875,876,877,878,879,880,881,882,883,884,885,886,887,888,889,890,891,892,893,894,895,896,897,898,899,900,901,902,1237,1238,1239,1240,1241,1242,1243,1244,1245,1246,1247,1248,1249,1250,1251,1252,1253,1254,1255,1256,1257,1258,1259,1260,1261,1262,1263,1264,1265,1266,1267,1268,1269,1270,1271,1272,1273,1274,1275,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,1318,1319,1320,1321,1322,1323,1324,1325,1326,1327,1328,1329,1330,1665,1666,1667,1668,1669,1670,1671,1672,1673,1674,1675,1676,1677,1678,1679,1680,1681,1682,1683,1684,1685,1686,1687,1688,1689,1690,1691,1692,1693,1694,1695,1696,1697,1698,1699,1700,1701,1702,1703,1720,1721,1722,1723,1724,1725,1726,1727,1728,1729,1730,1731,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,2095,2096,2097,2098,2099,2100,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,2127,2128,2129,2130,2131,2132,2133',
       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 03:54:59.940 -31.30.0.0000',
       scan=['34,37,9,12,15,18,21,24,28,31'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['VLASS1.1.sb34905760.eb35113377.58166.981752928245_split.ms'],
       imagename='VLASS1.1.ql.T03t06.J035459-313000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00046400396049903062, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,383,384,385,386,387,388,389,390,391,392,393,394,395,396,397,398,399,400,401,402,403,404,405,406,407,408,409,410,411,412,413,414,415,416,417,418,419,420,421,438,439,440,441,442,443,444,445,446,447,448,449,450,451,452,453,454,455,456,457,458,459,460,461,462,463,464,465,466,467,468,469,470,471,472,473,474,475,476,809,810,811,812,813,814,815,816,817,818,819,820,821,822,823,824,825,826,827,828,829,830,831,832,833,834,835,836,837,838,839,840,841,842,843,844,845,846,847,864,865,866,867,868,869,870,871,872,873,874,875,876,877,878,879,880,881,882,883,884,885,886,887,888,889,890,891,892,893,894,895,896,897,898,899,900,901,902,1237,1238,1239,1240,1241,1242,1243,1244,1245,1246,1247,1248,1249,1250,1251,1252,1253,1254,1255,1256,1257,1258,1259,1260,1261,1262,1263,1264,1265,1266,1267,1268,1269,1270,1271,1272,1273,1274,1275,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,1318,1319,1320,1321,1322,1323,1324,1325,1326,1327,1328,1329,1330,1665,1666,1667,1668,1669,1670,1671,1672,1673,1674,1675,1676,1677,1678,1679,1680,1681,1682,1683,1684,1685,1686,1687,1688,1689,1690,1691,1692,1693,1694,1695,1696,1697,1698,1699,1700,1701,1702,1703,1720,1721,1722,1723,1724,1725,1726,1727,1728,1729,1730,1731,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,2095,2096,2097,2098,2099,2100,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,2127,2128,2129,2130,2131,2132,2133',
       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.T03t06.J035459-313000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T03t06.J035459-313000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T03t06.J035459-313000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#