sedimentation_matrix: Create Sedimentation and Migration Matrix for SALMO-1D

Description Usage Arguments Details Value Examples

Description

Note: state variable indexes still hard coded

Usage

1
2
sedimentation_matrix(parms, nstates = NULL, nphy = NULL, depth = NULL,
  focussing = c(1, 1))

Arguments

parms

SALMO parameter list

nstates

total number of state variables (deprecated !!!)

nphy

number of phytoplankton groups (deprecated !!!)

depth

vector of individual depths (m) (deprecated !!!)

focussing

optional two-valued vector for a linear sediment focussing heuristics (e.g. c(1, 2) - experimental!)

Details

ToDo: pass parameters regularly to this function

Value

matrix with sedimentation / migration velocities

Examples

1
2
3
4
5
nstates <- 8 # nOfVar$numberOfStates
nphy    <- 3 # nOfVar$numberOfAlgae
parms   <- get_salmo_parms(nlayers=140, macrophytes=TRUE)
depth  <- seq(0, 70, 0.5) # unique(inputs[,3])
vmat   <- sedimentation_matrix(parms, nstates, nphy, depth)

rSALMO documentation built on May 2, 2019, 6:12 p.m.