diff --git a/test_problems/diamondSurf/runDiamond.cpp b/test_problems/diamondSurf/runDiamond.cpp
index a380a5c40..738490ed9 100644
--- a/test_problems/diamondSurf/runDiamond.cpp
+++ b/test_problems/diamondSurf/runDiamond.cpp
@@ -49,18 +49,18 @@ int main(int argc, char** argv) {
       string path = findInputFile(infile);
       ctml::get_CTML_Tree(xc, path);
 
-      XML_Node *xg = findXMLPhase(xc, "gas");
+      XML_Node * const xg = (XML_Node *) findXMLPhase(xc, "gas");
       ThermoPhase *gasTP = newPhase(*xg);
       int nsp = gasTP->nSpecies();
       cout << "Number of species = " << nsp << endl;
 
-      XML_Node *xd = findXMLPhase(xc, "diamond");
+      XML_Node * const xd = (XML_Node *) findXMLPhase(xc, "diamond");
       ThermoPhase *diamondTP = newPhase(*xd);
       int nsp_diamond = diamondTP->nSpecies();
       cout << "Number of species in diamond = " << nsp_diamond << endl;
 
 
-      XML_Node *xs = findXMLPhase(xc, "diamond_100");
+      XML_Node * const xs = (XML_Node *) findXMLPhase(xc, "diamond_100");
       ThermoPhase *diamond100TP = newPhase(*xs);
       int nsp_d100 = diamond100TP->nSpecies();
       cout << "Number of species in diamond_100 = " << nsp_d100 << endl;
diff --git a/test_problems/diamondSurf/runDiamond_blessed.out b/test_problems/diamondSurf/runDiamond_blessed.out
index 5f7df8c72..b874439cc 100644
--- a/test_problems/diamondSurf/runDiamond_blessed.out
+++ b/test_problems/diamondSurf/runDiamond_blessed.out
@@ -2,11 +2,11 @@ Number of species = 4
 Number of species in diamond = 1
 Number of species in diamond_100 = 8
 Number of reactions = 20
-0  1  -8.95682e-05
-1  2  4.48368e-05
-2  3  -3.51474e-08
+0  1  -8.95751e-05
+1  2  4.48403e-05
+2  3  -3.51539e-08
 3  4   nil
-4  0  3.51474e-08
+4  0  3.51539e-08
 5  2   nil
 6  1   nil
 7  1   nil
@@ -15,14 +15,14 @@ Number of reactions = 20
 10  3   nil
 11  3   nil
 12  2   nil
-sum = -1.07461e-20
-growth rate = 0.43175 microns per hour
+sum = -2.53573e-21
+growth rate = 0.43183 microns per hour
 Coverages:
-0   c6HH   0.462324
-1   c6H*   0.0370432
-2   c6*H   0.474239
-3   c6**   0.0219335
-4   c6HM   0.00174798
-5   c6HM*   2.56191e-05
-6   c6*M   0.00264844
-7   c6B   3.81558e-05
+0   c6HH   0.462262
+1   c6H*   0.037052
+2   c6*H   0.474283
+3   c6**   0.0219445
+4   c6HM   0.00174648
+5   c6HM*   2.56272e-05
+6   c6*M   0.00264858
+7   c6B   3.8171e-05