Index: doc/release_notes.rst
===================================================================
diff -u -rfa4daa7879ffb68b5ef714c2604a0f0bbc691dea -ra9bbe9f0c2f701510e8c85b3b5319f63551ff29f
--- doc/release_notes.rst	(.../release_notes.rst)	(revision fa4daa7879ffb68b5ef714c2604a0f0bbc691dea)
+++ doc/release_notes.rst	(.../release_notes.rst)	(revision a9bbe9f0c2f701510e8c85b3b5319f63551ff29f)
@@ -13,7 +13,7 @@
 + Added bmi support
 + wflow_sbm s-curve sharpness of height distribution now based on upper and lower
   half of the distribution (average)
-+ Added seperate soil evap to wflow_sbm
++ Added separate soil evap and RootZoneSoilMoisture to wflow_sbm
 
 
 
Index: wflow-py/wflow/wflow_sbm.py
===================================================================
diff -u -rfa4daa7879ffb68b5ef714c2604a0f0bbc691dea -ra9bbe9f0c2f701510e8c85b3b5319f63551ff29f
--- wflow-py/wflow/wflow_sbm.py	(.../wflow_sbm.py)	(revision fa4daa7879ffb68b5ef714c2604a0f0bbc691dea)
+++ wflow-py/wflow/wflow_sbm.py	(.../wflow_sbm.py)	(revision a9bbe9f0c2f701510e8c85b3b5319f63551ff29f)
@@ -1221,6 +1221,8 @@
         else:
             self.reinfiltwater = self.ZeroMap
 
+
+        self.RootZonSoilMoisture = self.UStoreDepth * max(1.0, self.RootingDepth/self.zi)
         # The MAx here may lead to watbal error. Howvere, if inwaterMMM becomes < 0, the kinematic wave becomes very slow......
         self.InwaterMM = max(0.0,self.ExfiltWater + self.ExcessWater + self.SubCellRunoff + self.SubCellGWRunoff + self.RunoffOpenWater + self.BaseFlow - self.reinfiltwater - self.ActEvapOpenWater)
         self.Inwater = self.InwaterMM * self.ToCubic  # m3/s