Index: wflow/wflow_sbm.py =================================================================== diff -u -ra6b939e08b980c412cc05d2dc0bf707e7f49d2c4 -r11471548f856592485a77ca37df4ad838eed2eaf --- wflow/wflow_sbm.py (.../wflow_sbm.py) (revision a6b939e08b980c412cc05d2dc0bf707e7f49d2c4) +++ wflow/wflow_sbm.py (.../wflow_sbm.py) (revision 11471548f856592485a77ca37df4ad838eed2eaf) @@ -2430,7 +2430,7 @@ # Calculate simple recharge to groundwater self.Recharge = self.Transfer - self.CapFlux - self.ActEvapSat - + # Now add capflux to the layers one by one (from bottom to top) for n in np.arange(len(self.UStoreLayerThickness) - 1, -1, -1): @@ -2913,7 +2913,9 @@ self.InflowKinWaveCell = pcr.upstream( self.TopoLdd, self.OldSurfaceRunoff ) - deltasup = float(pcr.mapmaximum(pcr.abs(oldsup - self.SurfaceWaterSupply))) + deltasup = float( + pcr.mapmaximum(pcr.abs(oldsup - self.SurfaceWaterSupply)) + ) if deltasup < self.breakoff or self.nrit >= self.maxitsupply: break