rmt_convert_btn_to_dis: Convert RMT3DMS btn to RMODFLOW dis object

In rogiersbart/RMT3DMS: Pre- and post-processing of MT3DMS & MT3D-USGS files

Description

Convert RMT3DMS btn to RMODFLOW dis object

Usage

rmt_convert_btn_to_dis(
  btn,
  nper = btn$nper,
  perlen = btn$perlen[1:nper],
  nstp = btn$nstp[1:nper],
  tsmult = btn$tsmult[1:nper],
  sstr = c("SS", rep("TR", nper - 1)),
  prj = RMODFLOW::rmf_get_prj(btn),
  ...
)

Arguments

btn	RMT3DMS btn object
nper	number of stress periods in the dis object. Defaults to btn$nper. See details.
perlen	vector of stress period lengths in the dis object. Defaults to btn$perlen[1:nper]. See details.
nstp	vector of stress period time steps in the dis object. Defaults to btn$nstp[1:nper]. See details.
tsmult	vector of stress period time step multipliers in the dis object. Defaults to btn$tsmult[1:nper]. See details.
sstr	character vector with steady state ('SS') or transient ('TR') stress period indicator; defaults to first period 'SS' and all others 'TR'
prj	RMODFLOW prj object. By default obtained from the btn object.
...	additional arguments passed to RMODFLOW::rmf_create_dis

Details

Conversion works best when number of stress-periods is equal in the dis and btn objects. nper, perlen, nstp & tsmult should only differ from the defaults when the flow model is steady-state with only 1 stress period. In that case, MT3DMS allows for differing stress periods between the flow and transport model.

Value

RMODFLOW dis object

See Also

rmt_convert_dis_to_btn

Examples

btn <- rmt_create_btn()
rmt_convert_btn_to_dis(btn)

# BTN with 2 stress periods but flow model is single stress period steady-state
btn <- rmt_create_btn(nper = 2, perlen = c(1000, 500))
rmt_convert_btn_to_dis(btn, nper = 1, perlen = 1, nstp = 1, tsmult = 1)

# otherwise converts to dis object with only first stress period steady-state
rmt_convert_btn_to_dis(btn)

rogiersbart/RMT3DMS documentation built on Oct. 16, 2021, 9:45 a.m.