rmf_create_dis: Create an 'RMODFLOW' dis object

In rogiersbart/RMODFLOW: Pre- and post-processing of MODFLOW files

Create an RMODFLOW dis object

Description

rmf_create_dis creates an RMODFLOW dis object.

Usage

rmf_create_dis(
  nlay = 3,
  nrow = 10,
  ncol = 10,
  nper = 1,
  itmuni = 1,
  lenuni = 2,
  laycbd = rep(0, nlay),
  delr = 100,
  delc = 100,
  top = matrix(0, nrow = nrow, ncol = ncol),
  botm = array(rep(seq(0, -10 * (nlay + sum(laycbd != 0)), length = nlay + sum(laycbd !=
    0) + 1)[2:(nlay + sum(laycbd != 0) + 1)], each = nrow * ncol), dim = c(nrow, ncol,
    nlay + sum(laycbd != 0))),
  perlen = rep(1, nper),
  nstp = rep(1, nper),
  tsmult = rep(1, nper),
  sstr = c("SS", rep("TR", nper - 1)),
  prj = NULL
)

Arguments

nlay	number of layers; defaults to 3
nrow	number of rows; defaults to 10
ncol	number of columns; defaults to 10
nper	number of stress periods; defaults to 1
itmuni	time unit; defaults to 1 (seconds)
lenuni	length unit; defaults to 2 (metres)
laycbd	vector of Quasi-3D confining bed flags; defaults to 0 for all layers
delr	vector of cell widths along rows; defaults to 100 for all columns
delc	vector of cell widths along columns; defaults to 100 for all rows
top	matrix with the top elevation of layer 1; defaults to 0 for all nrow x ncol cells
botm	3D array with the bottom elevations of all layers and Quasi-3D confining beds; defaults to nlay layers, equally spaced between 0 (top first layer) and -100 (bottom last layer)
perlen	vector of stress period lengths
nstp	vector of stress period time steps
tsmult	vector of successive time step length multipliers
sstr	character vector with steady state ('SS') or transient ('TR') stress period indicator
prj	optional RMODFLOW prj object as created by rmf_create_prj to set the projection information. Defaults to NULL

Value

Object of class dis

See Also

rmf_read_dis, rmf_write_dis and http://water.usgs.gov/nrp/gwsoftware/modflow2000/MFDOC/index.html?dis.htm

rogiersbart/RMODFLOW documentation built on Jan. 14, 2023, 4:21 a.m.