Convert cbc object fluxes to darcy velocities

Source: R/utils.R

Convert cbc object fluxes to darcy velocities

rmf_convert_cbc_to_darcy(
  cbc,
  dis,
  hed = NULL,
  porosity = NULL,
  include_bc = TRUE
)

Arguments

cbc

RMODFLOW cbc object
dis

RMODFLOW dis object
hed

RMODFLOW hed object; optional; if specified, the saturated cell thickness is used
porosity

optional 3d array with porosity values. If used, the average linear groundwater flow velocities are returned.
include_bc

logical, should fluxes from boundary conditions be included in the calculation of Darcy fluxes ? Defaults to TRUE.

Value

list of 4d arrays: right, front, lower, left, back, upper, qx, qy, qz and q; all represent Darcy velocities (or average linear groundwater flow velocities if porosity is specified): the first six at the different cell faces, the last four represent the components and magnitude at the cell center