satthick {raem} | R Documentation |
Compute the saturated thickness
Description
satthick()
computes the saturated thickness of the aquifer from an aem
object
at the given x and y coordinates.
Usage
satthick(aem, x, y, as.grid = FALSE, ...)
Arguments
aem |
|
x |
numeric x coordinates to evaluate at. |
y |
numeric y coordinates to evaluate at. |
as.grid |
logical, should a matrix be returned? Defaults to |
... |
additional arguments passed to |
Details
If the aquifer is confined at x
and y
, the saturated thickness equals the aquifer thickness.
For flow with variable saturated thickness (aem$type = 'variable'
), if the aquifer is unconfined at x
and y
,
the saturated thickness is calculated as the hydraulic head at x
and y
minus the aquifer base.
Value
A vector of length(x)
(equal to length(y)
) with the saturated thicknesses at x
and y
.
If as.grid = TRUE
, a matrix of dimensions c(length(y), length(x))
described by
marginal vectors x
and y
containing the saturated thicknesses at the grid points.
See Also
Examples
uf <- uniformflow(100, 0.001, 0)
rf <- constant(-1000, 0, 11)
m <- aem(k = 10, top = 10, base = 0, n = 0.2, uf, rf, type = 'confined')
satthick(m, x = c(-200, 0, 200), y = 0) # confined
s <- satthick(m, x = seq(-500, 500, length = 100),
y = seq(-250, 250, length = 100), as.grid = TRUE)
str(s)
mv <- aem(k = 10, top = 10, base = 0, n = 0.2, uf, rf, type = 'variable')
satthick(mv, x = c(-200, 0, 200), y = 0) # variable