qh: Horizontal Darcy Velocity

In CJBarry/DNAPL: DNAPL Source Term Framework

Description

Calculates the horizontal Darcy velocity, or specific discharge, at a specific location in a MODFLOW model.

Usage

qh(x, y = NULL, L, mfdata)

Arguments

x, y	location, read by xy.coords; should be in absolute co-ordinates with the same origin and units as mfdata
L	integer []; a continuous range is insisted upon and this is corrected to min(L):max(L) within function
mfdata	NetCDF object; a MODFLOW dataset; see GW.nc

Details

Uses linear interpolation to calculate the velocity between cell faces. For the layer with the water table, the thickness (and therefore flow normal area) is transient. NA is returned for dry layers.

Value

numeric array [length(L), NTS], dropped, where NTS is the number of time steps in the MODFLOW results; the Darcy velocity in each layer for each time step at the location x,y

Examples

library("RNetCDF")

# get example MODFLOW data set
mfdata <- open.nc(system.file("rflow_mf_demo.nc", package = "Rflow"))
qh(625, 825, 1L, mfdata)

CJBarry/DNAPL documentation built on May 6, 2019, 9:25 a.m.