The function grd.Magkap
computes the maximum capacity for each saturation level
grd.Magkap(NoL, Res_prob, CapacityUpperLevel)
The output is a scalar
## Not run: ------------------------------------ # grd.Magkap() ## ---------------------------------------------