Description Usage Arguments Details Value See Also
rmt_create_rct
creates an RMT3DMS
rct object.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 | rmt_create_rct(
btn,
isothm = 0,
ireact = 0,
irctop = 2,
igetsc = 0,
ireaction = 0,
rhob = 157,
prsity2 = 0.1,
srconc = as.list(rep(0, btn$ncomp)),
sp1 = as.list(rep(0, btn$ncomp)),
sp1im = as.list(rep(0, btn$ncomp)),
sp2 = as.list(rep(0, btn$ncomp)),
rc1 = as.list(rep(0, btn$ncomp)),
rc2 = as.list(rep(0, btn$ncomp)),
rc3 = as.list(rep(0, btn$ncomp)),
yld = rep(1, btn$ncomp - 1),
ied = 1,
iea = 2,
f = 1,
rec_filename = "reaction_module.dat",
ned = 1,
nea = 5,
nspecial = 0,
ifesld = 0,
ispec = rep(0, nspecial),
specialispec = rep("SOLID", nspecial),
efcmax = rep(0, nspecial),
hsc = rep(0.5, nea),
ic = rep(0.01, nea),
decayrate = matrix(0.0186, nrow = nea, ncol = ned),
yieldc = matrix(1, nrow = nea + ned, ncol = ned)
)
|
btn |
|
isothm |
integer flag to set type of sorption (or dual-domain mass transfer). See MT3DMS/MT3D-USGS manual. Defaults to 0. |
ireact |
integer flag indicating which type of kinetic rate reaction is simulated. See MT3DMS/MT3D-USGS manual. Defaults to 0. |
irctop |
integer flag indicating how reaction variables are entered. See MT3DMS/MT3D-USGS manual. Defaults to 2 (as 3D-arrays). |
igetsc |
integer flag specifying if initial concentration of nonequilibrium sorbed or immobile phases should be read. See MT3DMS/MT3D-USGS manual. Defaults to 0 (not read). |
ireaction |
MT3D-USGS only. Integer flag specifying a reaction module. See MT3D-USGS manual. Defaults to 0 (no reaction module. |
rhob |
3d array with bulk density of aquifer medium. Defaults to 157 for all cells. |
prsity2 |
3d array with porosity values for the immobile domain. Defaults to 0.1 for all cells. |
srconc |
list with |
sp1 |
list with |
sp1im |
list with |
sp2 |
list with |
rc1 |
list with |
rc2 |
list with |
rc3 |
list with |
yld |
numeric vector of length |
ied |
integer species number representing the electron donor in a EA/ED reaction ( |
iea |
integer species number representing the electron acceptor in a EA/ED reaction ( |
f |
stoiciometric ratio in the simulated equation |
rec_filename |
character, name of the input file that will provide parameter information when |
ned |
number of electron donors when |
nea |
number of electron acceptors when |
nspecial |
number of special cases when |
ifesld |
integer flag specifying if solid phase Fe(III) is simulated. Defaults to 0 (not simulated). |
ispec |
numeric vector of length |
specialispec |
character vector of length |
efcmax |
numeric vector of length |
hsc |
numeric vector of length |
ic |
numeric vector of length |
decayrate |
matrix with |
yieldc |
matrix with |
If irctop < 2
, all 3D variables should be given as 1D vectors with one value per layer.
an object of class rct
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.