The function grd.S() computes the deficit S for all sub surface layers (not the overland flow layer)

grd.S(ddistx, X, ddist)

Arguments

ddistx
distance distribution of each saturation level
X
excess water, scalar
ddist
states (in mm) of each saturation level

Value

The output is a scalar

Details

It follows the expression dS/dt = -dD/dt from Skaugen,Peerebom and Nilsson(2015)

Examples

## Not run: ------------------------------------ # grd.S() ## ---------------------------------------------