dirflow {raem} | R Documentation |
Compute flow in the direction of a given angle
Description
dirflow()
computes a flow variable at the given points in the direction of the supplied angle.
Usage
dirflow(
aem,
x,
y,
angle,
flow = c("discharge", "darcy", "velocity"),
as.grid = FALSE,
...
)
Arguments
aem |
|
x |
numeric x coordinates to evaluate |
y |
numeric y coordinates to evaluate |
angle |
numeric, angle of the direction to evaluate |
flow |
character specifying which flow variable to use. Possible values are |
as.grid |
logical, should a matrix be returned? Defaults to FALSE. See details. |
... |
additional arguments passed to |
Details
The x and y components of flow
are used to calculate the directed value using angle
.
The z
coordinate in discharge()
, darcy()
or velocity()
is set at the aquifer base. Under Dupuit-Forchheimer,
the x and y components of the flow vector do not change along the vertical axis.
Value
A vector of length(x)
(equal to length(y)
) with the flow values at x
and y
in the direction of angle
.
If as.grid = TRUE
, a matrix of dimensions c(length(y), length(x))
described by
marginal vectors x
and y
containing the directed flow values at the grid points.
See Also
Examples
rf <- constant(-1000, 0, hc = 10)
uf <- uniformflow(TR = 100, gradient = 0.001, angle = -45)
w <- well(10, -50, Q = 200)
m <- aem(k = 10, top = 10, base = 0, n = 0.2, rf, uf)
dirflow(m, x = c(0, 100), y = 50, angle = -45)
m <- aem(k = 10, top = 10, base = 0, n = 0.2, rf, uf, w, type = 'confined')
dirflow(m, x = c(0, 50, 100), y = c(0, 50), angle = -90,
flow = 'velocity', as.grid = TRUE)