Create Sedimentation and Migration Matrix for SALMO-1D

Share:

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)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.