diff --git a/landlab/components/river_flow_dynamics/river_flow_dynamics.py b/landlab/components/river_flow_dynamics/river_flow_dynamics.py
index 7955229c72..5d09966cb2 100644
--- a/landlab/components/river_flow_dynamics/river_flow_dynamics.py
+++ b/landlab/components/river_flow_dynamics/river_flow_dynamics.py
@@ -118,7 +118,7 @@
 
 >>> flow_depth = np.reshape(grid['node']["surface_water__depth"],(nRows,nCols))[10,:]
 
->>> np.round(flow_depth,2)
+>>> np.round(flow_depth,3)
 array([0.5  , 0.491, 0.48 , 0.473, 0.467, 0.464, 0.46 , 0.458, 0.455,
        0.454, 0.452, 0.45 , 0.449, 0.448, 0.446, 0.445, 0.443, 0.442,
        0.441, 0.439, 0.438, 0.437, 0.435, 0.434, 0.433, 0.431, 0.43 ,
@@ -1617,7 +1617,7 @@ def run_one_step(self):
             left_hand_side,
             right_hand_side,
             M=Mi,
-            tol=self._pcg_tolerance,
+            #rtol=self._pcg_tolerance,
             maxiter=self._pcg_max_iterations,
             atol=0,
         )