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

# hif_editimlist(parameter_file='parameter.7795.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_T27t15.J233459+653000_P7795v1_2017_10_05T16_52_09.012/working/VLASS1.1.sb34461205.eb34469117.58017.220875416664_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T27t15.J233459+653000_P7795v1_2017_10_05T16_52_09.012/rawdata/VLASS1.1.sb34461205.eb34469117.58017.220875416664.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='1263,1267,1266,1265,1264,1669,1668,1667,1666,1665,1664,1663,270,271,272,273,274,275,276,277,278,279,2169,2168,2165,2164,2167,2166,2161,2160,2163,2162,153,559,558,555,554,557,556,551,553,552,152,1692,1693,1690,1691,155,1278,1279,1175,1174,1173,1172,1171,1170,1270,1271,1272,1273,1274,1275,1179,1178,1285,1284,1287,1286,1281,1280,1283,1282,1289,1288,1579,1578,688,1573,1572,1571,1570,1577,1576,1575,1574,156,2154,2155,2156,2157,2153,2158,2159,1955,1957,1956,1959,1958,1177,1176,144,145,142,143,141,580,581,1276,149,249,172,171,170,977,976,975,974,973,972,971,970,979,978,658,659,1292,1968,1969,1964,1965,1966,1967,1960,1961,1962,1963,1293,685,684,687,686,681,680,1681,1680,1683,1682,1685,1684,1687,1686,1689,1688,682,1269,1268,669,668,667,666,665,664,663,662,661,660,148,151,579,578,150,573,572,571,570,577,576,575,574,258,259,252,253,250,251,256,257,254,255,1588,1162,1582,1583,1580,1581,1586,1587,1584,1585,168,169,1675,164,165,166,167,160,161,162,163,964,965,966,967,960,961,962,963,968,969,1159,1158,1277,683,1979,1978,1977,1976,1975,1974,1973,1972,1971,1970,1678,1679,1674,280,1676,1677,1670,1671,1672,1673,678,679,674,675,676,677,670,671,672,673,263,262,261,260,267,266,265,264,269,268,1290,1291,1559,1558,982,983,980,981,986,987,984,985,2172,2173,2170,2171,2176,2177,2174,2175,2178,2179,146,147,1986,1984,1985,1982,1983,1980,1981,568,569,560,561,562,563,564,565,566,567,2183,2182,2181,2180,1188,1186,1187,1184,1185,1182,1183,1180,1181,1164,1165,1166,1167,1160,154,157,1163,159,158,1168,1169,1161,959,958,957,956,1568,1569,1560,1561,1562,1563,1564,1565,1566,1567',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 23:34:59.163 +65.30.0.0000',
       scan=['32,42,35,38,10,13,16,19,22,26,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_T27t15.J233459+653000_P7795v1_2017_10_05T16_52_09.012/working/VLASS1.1.sb34461205.eb34469117.58017.220875416664_split.ms'],
       imagename='VLASS1.1.ql.T27t15.J233459+653000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='141,142,143,144,145,146,147,148,149,150,151,152,153,154,155,156,157,158,159,160,161,162,163,164,165,166,167,168,169,170,171,172,249,250,251,252,253,254,255,256,257,258,259,260,261,262,263,264,265,266,267,268,269,270,271,272,273,274,275,276,277,278,279,280,551,552,553,554,555,556,557,558,559,560,561,562,563,564,565,566,567,568,569,570,571,572,573,574,575,576,577,578,579,580,581,658,659,660,661,662,663,664,665,666,667,668,669,670,671,672,673,674,675,676,677,678,679,680,681,682,683,684,685,686,687,688,956,957,958,959,960,961,962,963,964,965,966,967,968,969,970,971,972,973,974,975,976,977,978,979,980,981,982,983,984,985,986,987,1158,1159,1160,1161,1162,1163,1164,1165,1166,1167,1168,1169,1170,1171,1172,1173,1174,1175,1176,1177,1178,1179,1180,1181,1182,1183,1184,1185,1186,1187,1188,1263,1264,1265,1266,1267,1268,1269,1270,1271,1272,1273,1274,1275,1276,1277,1278,1279,1280,1281,1282,1283,1284,1285,1286,1287,1288,1289,1290,1291,1292,1293,1558,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,1663,1664,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,1955,1956,1957,1958,1959,1960,1961,1962,1963,1964,1965,1966,1967,1968,1969,1970,1971,1972,1973,1974,1975,1976,1977,1978,1979,1980,1981,1982,1983,1984,1985,1986,2153,2154,2155,2156,2157,2158,2159,2160,2161,2162,2163,2164,2165,2166,2167,2168,2169,2170,2171,2172,2173,2174,2175,2176,2177,2178,2179,2180,2181,2182,2183',
       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 23:34:59.163 +65.30.0.0000',
       scan=['32,42,35,38,10,13,16,19,22,26,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_T27t15.J233459+653000_P7795v1_2017_10_05T16_52_09.012/working/VLASS1.1.sb34461205.eb34469117.58017.220875416664_split.ms'],
       imagename='VLASS1.1.ql.T27t15.J233459+653000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00057206724878051311, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='141,142,143,144,145,146,147,148,149,150,151,152,153,154,155,156,157,158,159,160,161,162,163,164,165,166,167,168,169,170,171,172,249,250,251,252,253,254,255,256,257,258,259,260,261,262,263,264,265,266,267,268,269,270,271,272,273,274,275,276,277,278,279,280,551,552,553,554,555,556,557,558,559,560,561,562,563,564,565,566,567,568,569,570,571,572,573,574,575,576,577,578,579,580,581,658,659,660,661,662,663,664,665,666,667,668,669,670,671,672,673,674,675,676,677,678,679,680,681,682,683,684,685,686,687,688,956,957,958,959,960,961,962,963,964,965,966,967,968,969,970,971,972,973,974,975,976,977,978,979,980,981,982,983,984,985,986,987,1158,1159,1160,1161,1162,1163,1164,1165,1166,1167,1168,1169,1170,1171,1172,1173,1174,1175,1176,1177,1178,1179,1180,1181,1182,1183,1184,1185,1186,1187,1188,1263,1264,1265,1266,1267,1268,1269,1270,1271,1272,1273,1274,1275,1276,1277,1278,1279,1280,1281,1282,1283,1284,1285,1286,1287,1288,1289,1290,1291,1292,1293,1558,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,1663,1664,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,1955,1956,1957,1958,1959,1960,1961,1962,1963,1964,1965,1966,1967,1968,1969,1970,1971,1972,1973,1974,1975,1976,1977,1978,1979,1980,1981,1982,1983,1984,1985,1986,2153,2154,2155,2156,2157,2158,2159,2160,2161,2162,2163,2164,2165,2166,2167,2168,2169,2170,2171,2172,2173,2174,2175,2176,2177,2178,2179,2180,2181,2182,2183',
       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.T27t15.J233459+653000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T27t15.J233459+653000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T27t15.J233459+653000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#