rmt_convert_dis_to_btn: Convert RMT3DMS dis to RMT3DMS btn object
In rogiersbart/RMT3DMS: Pre- and post-processing of MT3DMS & MT3D-USGS files

Description Usage Arguments Details Value See Also Examples

Convert RMT3DMS dis to RMT3DMS btn object

rmt_convert_dis_to_btn(
  dis,
  nper = dis$nper,
  perlen = dis$perlen[1:nper],
  tunit = NULL,
  lunit = NULL,
  prj = RMODFLOW::rmf_get_prj(dis),
  ...
)

`dis`	`RMODFLOW` dis object
`nper`	number of stress periods in the btn object. Defaults to `dis$nper`. See details.
`perlen`	vector of stress period lengths in the btn object. Defaults to `dis$perlen[1:nper]`. See details.
`tunit`	character specifying the length unit. The default (NULL), guesses from the dis object. Note that lunit is for identification purposes only does not affect the model outcome.
`prj`	`RMODFLOW` prj object. By default, obtained from the dis object.
`...`	additional arguments passed to `rmt_create_btn`

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.

RMT3DMS btn object

rmt_convert_btn_to_dis, rmt_create_btn

library(RMODFLOW)
dis <- rmf_read_dis(rmf_example_file('example-model.dis'))
rmt_convert_dis_to_btn(dis)
rmt_convert_dis_to_btn(dis, nper = 2, perlen = c(1000, 500), tunit = 'd', ncomp = 2)