D50_to_porosity: Calculate sediment porosity

In Jack-H-Laverick/MiMeMo.tools: Tools for extracting data for StrathE2E, during MiMeMo

Calculate sediment porosity

Description

This function takes median grain size at a location and returns a measure of porosity. The default parameterisation comes from Matt Pace's thesis for fine sediments. The numerical order of parameters has been changed from Matt's to be consistent with StrathE2E.

Usage

D50_to_porosity(D50, p1 = -1.035, p2 = -0.314, p3 = -0.435, p4 = 0.302)

Arguments

D50	The median grain size at a location.
p1	A parameter in the relationship.
p2	A parameter in the relationship.
p3	A parameter in the relationship.
p4	A parameter in the relationship.

Value

A numeric vector containing estimates of sediment porosity in

See Also

Other Sediment properties: mud_to_permeability()

Examples

data <- data.frame(D50 = seq(0, 1, length.out = 1000)) %>%
mutate(Porosity = D50_to_porosity(D50))

ggplot(data) +
geom_line(aes(x = D50, y = Porosity)) +
scale_x_continuous(trans = "log", breaks = c(0.001,0.01,0.1,1))

Jack-H-Laverick/MiMeMo.tools documentation built on March 6, 2023, 3:44 p.m.