For F2, the high-resolution grid has 221 points from R = (1.1 to 2.2) Å, a spacing of 0.005 Å. The coarse, irregular grid consists of the nine points with R = 1.28, 1.30, 1.32, 1.36, 1.42, 1.48, 1.54, 1.58 and 1.63 Å, which are close to the energy minimum and the classical turning points (on the high-level potential) for v = 0, 1, 2, and 3. Low-level methods are HF/6-31G*, B3LYP/6-31G*, fcMP2/cc-pVDZ and fcCCSD/cc-pVDZ, all spin-restricted, where the prefix "fc" indicates that core electrons are left uncorrelated ("frozen"). The high-level method is RHF-CCSDT/aug-cc-pwCVTZ with all electrons correlated. We are testing procedures, so a high-resolution grid is also computed at the high level to provide access to "correct" values. All calculations are done using Gaussian09 except for CCSDT, which is done using the CFOUR programs. For CO, the high-resolution grid has 201 points from R = (0.8 to 1.8) Å, a spacing of 0.005 Å. Guiding potentials are computed at: CASPT3(10,8)/aug-cc-pwCVTZ (third-order multireference perturbation theory ), icMRCI(10,8)+Q /aug-cc-pwCV5Z (internally contracted multireference singles and doubles configuration interaction with Davidson correction for additional correlation), and UHF-fcCCSD/aug-cc-pVQZ. The CCSD calculations are done using Gaussian09 and the multireference calculations are done using MOLPRO. For C2, the high-level data from Boschen et al. span a wide range, from 0.9 to 20 Å. Here, only distances up to 6 Å are considered. The high-resolution grid is from 0.9 to 6.0 Å in steps of 0.005 Å (1021 points). Inexpensive guiding potentials are CASSCF(8,8)/cc-pVTZ and icMRCI +Q/aug-cc-pwCVTZ using three-state-averaged orbitals. The X state is identified as the lower of the  = 0 states in the CASSCF, and as the icMRCI state dominated by that CASSCF reference. Abrams and Sherrill obtained a correctly shaped potential from RHF-fcFCI (full configuration interaction) using a small basis set. However, this is costly (about 100 times more than the present MRCI and 10 4 times more than the present CASSCF), so for this method we use a superposition of three grids: from 0.9 to 1.9 (R = 0.005 Å), from 1.9 to 3.2 (R = 0.02 Å), and from 3.2 to 6.0 (R = 0.1 Å) for a total of 294 points. We solve for the three lowest states (of irrep Ag in D2h), to get beyond the curve crossing near 1.7 Å while obtaining values for X Σ . The X state in the FCI calculations is identified by requiring smoothness in the energy and in the components of the quadrupole moment, although it is not always clear. These calculations, including the FCI, are done using MOLPRO. Data analysis, including splining and polynomial fitting, is done in Python using the standard Scipy and Numpy libraries.