Index: wflow-py/Sandbox/wflow_topoflex/reservoir_Su.py
===================================================================
diff -u -r5f72a83356683168a37f127593b6f33c4eccd3b5 -r878b78ee5bcc03356fd85011b4da54a9d9636813
--- wflow-py/Sandbox/wflow_topoflex/reservoir_Su.py	(.../reservoir_Su.py)	(revision 5f72a83356683168a37f127593b6f33c4eccd3b5)
+++ wflow-py/Sandbox/wflow_topoflex/reservoir_Su.py	(.../reservoir_Su.py)	(revision 878b78ee5bcc03356fd85011b4da54a9d9636813)
@@ -65,6 +65,22 @@
         name = 'unsatZone_LP_beta_Ep_cropG'
     return name
 
+def unsatZone_no_reservoir(self, k):
+    """
+    This function is used when no unsaturated zone reservoir is used and only
+    passes fluxes from the upper reservoirs to the lower
+    Qu = Pe
+    Eu = 0.
+    Perc = 0.
+    Cap = 0.
+    Storage in unsaturated zone = 0.
+    """
+    self.Qu_[k] = max(self.Pe_[k], 0)
+    self.Eu_[k] = 0.
+    self.Perc_[k] = 0.
+    self.Su_[k] = 0.
+    self.Cap_[k] = 0.
+    self.wbSu_[k] = self.Pe - self.Eu - self.Qu - self.Perc + self.Cap - self.Su[k] + self.Su_t[k]
 
 def unsatZone_LP_beta(self,k):
     """