R: Implementation of a general model for linear non-autonomous...

linearScalarModel {SoilR}

R Documentation

Implementation of a general model for linear non-autonomous systems with scalar modifiers

Description

This function implements a linear model with scalar modifier for inputs and compartmental matrix.

Usage

linearScalarModel(
  t,
  A,
  C0,
  u,
  gamma,
  xi,
  xi_lag = 0,
  solver = deSolve.lsoda.wrapper
)

Arguments

`t`	A vector containing the points in time where the solution is sought.
`A`	A square (n x n) matrix with compartmental structure
`C0`	A vector of length n containing the initial amount of carbon for the n pools.
`u`	A vector of length n with constant mass inputs for the n pools.
`gamma`	A scalar or data.frame object specifying the modifier for the mass inputs.
`xi`	A scalar, data.frame, function or anything that can be converted to a scalar function of time `ScalarTimeMap` object specifying the external (environmental and/or edaphic) effects on decomposition rates.
`xi_lag`	A time shift/delay for the automatically created time dependent function xi(t)
`solver`	A function that solves the system of ODEs. This can be `euler` or `deSolve.lsoda.wrapper` or any other user provided function with the same interface.

Value

A Model Object that can be further queried

References

C.A., M. Mueller, S.E. Trumbore. 2012. Models of soil organic matter decomposition: the SoilR package version 1.0. Geoscientific Model Development 5, 1045-1060.

Examples

t=seq(0,52*200,1) # Fix me! Add an example.

[Package SoilR version 1.2.105 Index]