Convert RMT3DMS btn to RMODFLOW dis object

Source: R/utils.R

Convert RMT3DMS btn to RMODFLOW dis object

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

Value

RMODFLOW dis object

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.

See also

rmt_convert_dis_to_btn

Examples

btn <- rmt_create_btn()
rmt_convert_btn_to_dis(btn)

#> RMODFLOW Discretization File object with: 
#> 10 rows, 10 columns and 3 layers totalling 300 cells 
#> 
#> Time units: seconds 
#> Length units: meters 
#> 
#> Spacing along rows (DELR):  100 (constant) 
#> Spacing along columns (DELC):  100 (constant) 
#> 
#> Summary of top elevations: 
#>       Top   
#>  Min.   :0  
#>  1st Qu.:0  
#>  Median :0  
#>  Mean   :0  
#>  3rd Qu.:0  
#>  Max.   :0  
#> 
#> Summary of bottom elevations: 
#>         Layer 1 Layer 2 Layer 3
#> Min.        -10     -20     -30
#> 1st Qu.     -10     -20     -30
#> Median      -10     -20     -30
#> Mean        -10     -20     -30
#> 3rd Qu.     -10     -20     -30
#> Max.        -10     -20     -30
#> 
#> Information for 1 stress-period: 
#>  Period Length Timesteps Multiplier         Type
#>       1      1         1          1 Steady-state

# 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)

#> RMODFLOW Discretization File object with: 
#> 10 rows, 10 columns and 3 layers totalling 300 cells 
#> 
#> Time units: seconds 
#> Length units: meters 
#> 
#> Spacing along rows (DELR):  100 (constant) 
#> Spacing along columns (DELC):  100 (constant) 
#> 
#> Summary of top elevations: 
#>       Top   
#>  Min.   :0  
#>  1st Qu.:0  
#>  Median :0  
#>  Mean   :0  
#>  3rd Qu.:0  
#>  Max.   :0  
#> 
#> Summary of bottom elevations: 
#>         Layer 1 Layer 2 Layer 3
#> Min.        -10     -20     -30
#> 1st Qu.     -10     -20     -30
#> Median      -10     -20     -30
#> Mean        -10     -20     -30
#> 3rd Qu.     -10     -20     -30
#> Max.        -10     -20     -30
#> 
#> Information for 1 stress-period: 
#>  Period Length Timesteps Multiplier         Type
#>       1      1         1          1 Steady-state

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

#> RMODFLOW Discretization File object with: 
#> 10 rows, 10 columns and 3 layers totalling 300 cells 
#> 
#> Time units: seconds 
#> Length units: meters 
#> 
#> Spacing along rows (DELR):  100 (constant) 
#> Spacing along columns (DELC):  100 (constant) 
#> 
#> Summary of top elevations: 
#>       Top   
#>  Min.   :0  
#>  1st Qu.:0  
#>  Median :0  
#>  Mean   :0  
#>  3rd Qu.:0  
#>  Max.   :0  
#> 
#> Summary of bottom elevations: 
#>         Layer 1 Layer 2 Layer 3
#> Min.        -10     -20     -30
#> 1st Qu.     -10     -20     -30
#> Median      -10     -20     -30
#> Mean        -10     -20     -30
#> 3rd Qu.     -10     -20     -30
#> Max.        -10     -20     -30
#> 
#> Information for 2 stress-periods: 
#>  Period Length Timesteps Multiplier         Type
#>       1   1000         1          1 Steady-state
#>       2    500         1          1    Transient