ODE-States consisting of several elements/isotopes.
This package supports the use of array based state variables with package deSolve.
This often occurs when several isotopes of the same element need to be calculated. The matrix of elements in rows and isotopes in columns, is converted to state vector with entries named: row_col
An example modMetaICBMDemo with the simple ICBM carbon turnover model is provided using radiocarbon data.
setting up a structure about the models Metadata twCreateModMeta
adding auxiliary outputs to the Metadata, twSetModMetaAuxGroups
initializing the state, initStateModMeta
converting state vector to state matrix, twStateMatODERes
constructing entry names twModElementNames
radiocarbon constants: c14Constants
radiocarbon atmospheric record of the bomb peak: delta14Catm
and function fmAtmosphere
.
Thomas Wutzler <twutz@bgc-jena.mpg.de>
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.