diff --git a/test_problems/surfSolverTest/results2_blessed.txt b/test_problems/surfSolverTest/results2_blessed.txt
index 12b175a33..415432bfc 100644
--- a/test_problems/surfSolverTest/results2_blessed.txt
+++ b/test_problems/surfSolverTest/results2_blessed.txt
@@ -32,13 +32,13 @@ Surface Phase:  soot_interface (11)
 Surface Temperature = 1.4e+03
 Surface Pressure    = 1.01e+05
                        Name         Coverage         SrcRate 
-   0                  Csoot-*        0.0184677    -8.195044e-20
-   1                  Csoot-H         0.981532    -8.195044e-20
+   0                  Csoot-*        0.0184677     0.000000e+00
+   1                  Csoot-H         0.981532     0.000000e+00
 Sum of coverages = 1
 Surface Phase:  soot_interface (11)
 Surface Temperature = 1.4e+03
 Surface Pressure    = 1.01e+05
                        Name         Coverage         SrcRate 
-   0                  Csoot-*        0.0184677    -8.195044e-20
-   1                  Csoot-H         0.981532    -8.195044e-20
+   0                  Csoot-*        0.0184677     0.000000e+00
+   1                  Csoot-H         0.981532     0.000000e+00
 Sum of coverages = 1
diff --git a/test_problems/surfSolverTest/surfaceSolver2.cpp b/test_problems/surfSolverTest/surfaceSolver2.cpp
index 4dbaff689..7794fd926 100644
--- a/test_problems/surfSolverTest/surfaceSolver2.cpp
+++ b/test_problems/surfSolverTest/surfaceSolver2.cpp
@@ -145,9 +145,13 @@ void printSurf(ostream &oooo, ThermoPhase *surfPhaseTP,
   int nspSurf = surfPhaseTP->nSpecies();
   for (int k = 0; k < nspSurf; k++) {
     kstart = iKin_ptr->kineticsSpeciesIndex(0, 2);
+    double srcK = src[kstart];
+    if (fabs(srcK) < 1.0E-8) {
+      srcK = 0.0;
+    }
     sprintf(buf, "%4d %24s   %14g   %14e\n", 
 	   k, surfPhaseTP->speciesName(k).c_str(),
-	   x[k], src[kstart]);
+	   x[k], srcK);
     oooo << buf;
     sum += x[k];
   }