UnsHydCond Function

public function UnsHydCond(ksat, theta, thetas, thetar, psdi) result(k)

Compute hydraulic conductivity of partially saturated soil (m/s)

Arguments

Type IntentOptional Attributes Name
real(kind=float), intent(in) :: ksat

saturated hydraulic conductivity (m/s)
real(kind=float), intent(in) :: theta

volumetric water content (m3/m3)
real(kind=float), intent(in) :: thetas

saturated volumetric water content (m3/m3)
real(kind=float), intent(in) :: thetar

residual volumetric water content (m3/m3)
real(kind=float), intent(in) :: psdi

Brooks & Corey pore size distribution index (-)

Return Value real(kind=float)


Source Code

FUNCTION UnsHydCond &
!
 (ksat, theta, thetas, thetar, psdi) &
 !
 RESULT (k)

IMPLICIT NONE

!Arguments with intent in
REAL (KIND = float), INTENT(IN) :: ksat !!saturated hydraulic conductivity (m/s)
REAL (KIND = float), INTENT(IN) :: theta !!volumetric water content (m3/m3)
REAL (KIND = float), INTENT(IN) :: thetas !!saturated volumetric water content (m3/m3)
REAL (KIND = float), INTENT(IN) :: thetar !!residual volumetric water content (m3/m3)
REAL (KIND = float), INTENT(IN) :: psdi !!Brooks & Corey pore size distribution index (-)


!local declarations:
REAL (KIND = float) :: k


!-------------------------end of declarations----------------------------------

k = ksat * ( (theta - thetar) / (thetas - thetar) ) ** (2./psdi + 3.)


RETURN

END FUNCTION UnsHydCond