Index: wflow-py/wflow/wflow_sbm.py =================================================================== diff -u -r3736503d1154e7e76dcf4e1a9c11d7d8cc277d83 -re2a41c70e2142cd8157fdcc4264039e194e0522d --- wflow-py/wflow/wflow_sbm.py (.../wflow_sbm.py) (revision 3736503d1154e7e76dcf4e1a9c11d7d8cc277d83) +++ wflow-py/wflow/wflow_sbm.py (.../wflow_sbm.py) (revision e2a41c70e2142cd8157fdcc4264039e194e0522d) @@ -719,7 +719,10 @@ # Used to seperate output per LandUse/management classes OutZones = self.LandUse - self.QMMConv = self.timestepsecs / (self.reallength * self.reallength * 0.001) #m3/s --> mm + self.QMMConv = self.timestepsecs / (self.reallength * self.reallength * 0.001) #m3/s --> actial mm of water over the cell + #self.QMMConvUp = 1000.0 * self.timestepsecs / ( catchmenttotal(cover(1.0), self.TopoLdd) * self.reallength * self.reallength) #m3/s --> mm over upstreams + temp = catchmenttotal(cover(1.0), self.TopoLdd) * self.reallength * 0.001 * 0.001 * self.reallength + self.QMMConvUp = cover(self.timestepsecs * 0.001)/temp self.ToCubic = (self.reallength * self.reallength * 0.001) / self.timestepsecs # m3/s self.KinWaveVolume = self.ZeroMap self.OldKinWaveVolume = self.ZeroMap @@ -878,7 +881,7 @@ self.KinWaveVolume = self.WaterLevel * self.Bw * self.DCL self.OldKinWaveVolume = self.KinWaveVolume - self.SurfaceRunoffMM = self.SurfaceRunoff * self.QMMConv + self.QCatchmentMM = self.SurfaceRunoff * self.QMMConvUp self.InitialStorage = self.FirstZoneDepth + self.UStoreDepth + self.CanopyStorage + self.LowerZoneStorage self.CellStorage = self.FirstZoneDepth + self.UStoreDepth @@ -912,6 +915,7 @@ :var self.ActEvap: Actual evaporation (transpiration + Soil evap + open water evap) [mm] :var self.SurfaceRunoff: Surface runoff in the kinematic wave [m^3/s] :var self.SurfaceRunoffDyn: Surface runoff in the dyn-wave resrvoir [m^3/s] + :var self.SurfaceRunoffCatchmentMM: Surface runoff in the dyn-wave reservoir expressed in mm over the upstream (catchment) area :var self.WaterLevelDyn: Water level in the dyn-wave resrvoir [m^] :var self.ActEvap: Actual EvapoTranspiration [mm] (minus interception losses) :var self.ExcessWater: Water that cannot infiltrate due to saturated soil [mm] @@ -1249,6 +1253,7 @@ self.MassBalKinWave = (self.KinWaveVolume - self.OldKinWaveVolume) / self.timestepsecs + self.InflowKinWaveCell + self.Inwater - self.SurfaceRunoff Runoff = self.SurfaceRunoff + self.QCatchmentMM = self.SurfaceRunoff * self.QMMConvUp # Updating # --------