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

In SoilR: Models of Soil Organic Matter Decomposition

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.

See Also

RothCModel. There are other predefinedModels and also more general functions like Model.

Examples

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

SoilR documentation built on Oct. 13, 2023, 5:06 p.m.