pastas.recharge.FlexModel.get_root_zone_balance#

static FlexModel.get_root_zone_balance(pe, ep, srmax=250.0, lp=0.25, ks=100.0, gamma=4.0, dt=1.0)[source]#

Method to compute the water balance of the root zone reservoir.

Parameters
  • pe (numpy.Array) – Effective precipitation flux in mm/d.

  • ep (numpy.Array) – Potential evaporation flux in mm/d.

  • srmax (float, optional) – Maximum storage capacity of the root zone.

  • lp (float, optional) – Parameter determining when actual evaporation equals potential.

  • ks (float, optional) – Saturated hydraulic conductivity in mm/d.

  • gamma (float, optional) – Parameter determining the nonlinearity of outflow / recharge.

  • dt (float, optional) – time step for the calculation of the recharge. Only dt=1 is possible now.

Returns

  • sr (numpy.Array) – Storage in the root zone reservoir.

  • r (numpy.Array) – Recharge flux in mm/d

  • ea (numpy.Array) – Evaporation flux in mm/d. Consists of transpiration and soil evaporation. Does not include interception evaporation.

  • q (numpy.Array) – surface runoff flux in mm/d.

  • pe (numpy.Array) – Incoming infiltration flux in mm/d.

Notes

If Numba is available, this method is significantly faster.