Index: wflow-py/wflow/wflow_sbm.py
===================================================================
diff -u -r32ef34c05a505cf391a132636c035b99789a626b -rc3e926a4c0a1df8a776e53de1aaeb17fda07b6b0
--- wflow-py/wflow/wflow_sbm.py	(.../wflow_sbm.py)	(revision 32ef34c05a505cf391a132636c035b99789a626b)
+++ wflow-py/wflow/wflow_sbm.py	(.../wflow_sbm.py)	(revision c3e926a4c0a1df8a776e53de1aaeb17fda07b6b0)
@@ -693,9 +693,11 @@
         # Use supplied riverwidth if possible, else calulate
         self.RiverWidth = ifthenelse(self.RiverWidth <= 0.0, W, self.RiverWidth)
 
-        # Only allow reinfiltration in river cells
-        self.MaxReinfilt = ifthenelse(self.River, self.ZeroMap + 999.0, self.ZeroMap)
+        # Only allow reinfiltration in river cells by default
 
+        if not hasattr(self,'MaxReinfilt'):
+            self.MaxReinfilt = ifthenelse(self.River, self.ZeroMap + 999.0, self.ZeroMap)
+
         # soil thickness based on topographical index (see Environmental modelling: finding simplicity in complexity)
         # 1: calculate wetness index
         # 2: Scale the capacity (now actually a max capacity) based on the index, also apply a minmum capacity
@@ -1436,7 +1438,7 @@
         ##########################################################################
 
         self.QCatchmentMM = self.SurfaceRunoff * self.QMMConvUp
-        #self.RunoffCoeff = self.QCatchmentMM/catchmenttotal(self.PrecipitationPlusMelt, self.TopoLdd)/catchmenttotal(cover(1.0), self.TopoLdd)
+        self.RunoffCoeff = self.QCatchmentMM/catchmenttotal(self.PrecipitationPlusMelt, self.TopoLdd)/catchmenttotal(cover(1.0), self.TopoLdd)
         #self.AA = catchmenttotal(self.PrecipitationPlusMelt, self.TopoLdd)
         #self.BB = catchmenttotal(cover(1.0), self.TopoLdd)
         # Single cell based water budget. snow not included yet.