Dtheta: Hydraulic diffusivity

Description Usage Arguments Details Value See Also Examples

Description

Dtheta calculates hydrualic diffusivity

Usage

1
Dtheta(sm, sp, timestep = 86400)

Arguments

sm

Soil water fraction

sp

a list of van Genuchten soil parameters as returned by getsoilparams()

timestep

number of seconds in each time step of model run (default one day)

Details

Applies the van Genuchten equation to estimate the hydrualic diffusivity (van Genuchten 1980 Soil Science Society of America Journal 44:892-898). Paramaters of the van Genuchten equation can be found in the dataset soilparams(). Ksat and the returned value are adjusted so that returned units are in cm per timestep. Dtheta is contrained to be finite by ensuring that sm marginally less than Smax when saturated and is assumed to be zero, when sm = Smin.

Value

Hydraulic conductivity (cm / time-step)

See Also

soilparams()

Examples

1
2
sm <- c(74:422) / 1000
plot(log(Dtheta(sm, getsoilparams("Loam")), 10) ~ sm)

ilyamaclean/ecohydrotools documentation built on June 10, 2019, 5:45 a.m.