Saturation layers The function updates the saturation levels. It computes the runoff event that occurs in each of layer, shift the layer by one timestep and then update the layer value with the runoff value
do.groundwater.Layers(NoL, Layers, ddist, X, layerUH, nbStepsDelay)
The output is a matrix describing the groundwater. the number of column represent the delay-steps, the row represent the levels (1st is the fastest, the last is the slowest)
## Not run: ------------------------------------ # do.groundwater.Layers() ## ---------------------------------------------