Note
Go to the end to download the full example code.
FLOW: Vadose: van Genuchten#
Shows the water retention curve and the hydraulic conductivity function for a number of soil types.
For more information about the parameters used see:
The RETC code for quantifying the hydraulic functions of unsaturated soils, Van Genuchten, M Th, Leij, F J, Yates, S R
import matplotlib.pyplot as plt
import discretize
from simpeg.flow import richards
def run(plotIt=True):
mesh = discretize.TensorMesh([10])
VGparams = richards.empirical.VanGenuchtenParams()
leg = []
for p in dir(VGparams):
if p[0] == "_":
continue
leg += [p]
params = getattr(VGparams, p)
k_fun, theta_fun = richards.empirical.van_genuchten(mesh, **params)
theta_fun.plot(ax=plt.subplot(121))
k_fun.plot(ax=plt.subplot(122))
plt.legend(leg, loc=3)
if __name__ == "__main__":
run()
plt.show()
Total running time of the script: (0 minutes 0.618 seconds)
Estimated memory usage: 163 MB