def PlotProps(pars):
import numpy as np
import pylab as pl
import vanGenuchten as vg
psi = np.linspace(-10, 2, 200)
pl.figure
pl.subplot(3, 1, 1)
pl.plot(psi, vg.thetaFun(psi, pars))
pl.ylabel(r'$\theta(\psi) [-]$')
pl.subplot(3, 1, 2)
pl.plot(psi, vg.CFun(psi, pars))
pl.ylabel(r'$C(\psi) [1/m]$')
pl.subplot(3, 1, 3)
pl.plot(psi, vg.KFun(psi, pars))
pl.xlabel(r'$\psi [m]$')
pl.ylabel(r'$K(\psi) [m/d]$')
def PlotProps(pars): import numpy as np import pylab as pl import vanGenuchten as vg psi=np.linspace(-10,2,200) pl.figure pl.subplot(3,1,1) pl.plot(psi,vg.thetaFun(psi,pars)) pl.ylabel(r'$\theta(\psi) [-]$') pl.subplot(3,1,2) pl.plot(psi,vg.CFun(psi,pars)) pl.ylabel(r'$C(\psi) [1/m]$') pl.subplot(3,1,3) pl.plot(psi,vg.KFun(psi,pars)) pl.xlabel(r'$\psi [m]$') pl.ylabel(r'$K(\psi) [m/d]$')
q[j] = -Kmid * ((psi[i + 1] - psi[i]) / dz + 1.0)
if q[n] != 0:
if q[n] < q[n - 1]: # because here q is a negative value
print('Water accumulated at the top, time = ', t / 60, 'min(s)')
KTop = vg.KFun(psiTop, p)
dhdz = (psiTop - psi[n - 1]) * 2 / dz + 1
q[n] = -KTop * dhdz
# Continuity
i = np.arange(0, n)
dpsidt = (-(q[i + 1] - q[i]) / dz) / C
return dpsidt
psi = np.linspace(-10, 5)
theta = vg.thetaFun(psi, p)
C = vg.CFun(psi, p)
K = vg.KFun(psi, p)
pl.figure()
pl.rcParams['figure.figsize'] = (5.0, 10.0)
pl.subplot(311)
pl.plot(psi, theta)
pl.ylabel(r'$\theta$', fontsize=20)
pl.subplot(312)
pl.plot(psi, C)
pl.ylabel(r'$C$', fontsize=20)
pl.subplot(313)
pl.plot(psi, K)
pl.ylabel(r'$K$', fontsize=20)
pl.xlabel(r'$\psi$', fontsize=20)