CNPDIAvars: Functions to retrieve variables for the CNPDIA and MPBDIA...

In CNPDIA: Diagenetic Model with Simple Nitrogen, Carbon, Phosphorus Dynamics

Description

CNPDIA0D, CNPDIA1D, CNPDIAsvar retrieve the (0-dimensional or 1-dimensional) output variables or the state variables of CNPDIA model solutions. When called with out not specified, will return the names of these variables and their units.

MPBDIA0D, MPBDIA1D, MPBDIAsvar retrieve the (0-dimensional or 1-dimensional) output variables or the state variables of MPBDIA model solutions. When called with out not specified, will return the names of these variables and their units.

Usage

  CNPDIA0D(out, as.vector = FALSE, which = NULL) 
  CNPDIA1D(out, which = NULL) 
  CNPDIAsvar(out, which = NULL) 
  MPBDIA0D(out, as.vector = FALSE, which = NULL) 
  MPBDIA1D(out, which = NULL) 
  MPBDIAsvar(out, which = NULL) 

Arguments

out	an output object returned by CNPDIAsolve, MPBDIAsolve, or CNPDIAdyna, MPBDIAdyna. If NULL, CNPDIAparms or MPBDIAparms will return the default (parameter) values.
as.vector	if TRUE will return the parameter vector, else a data.frame that also contains the units.
which	if not NULL, a vector with names of the variables/parameters to return.

Value

CNPDIA0D and CNPDIA1D return the output variables of the solution as a vector or data.frame. For dynamic runs, the output is averaged over the mean of the run.

CNPDIA1D always returns the sediment depth and the porosity as the first two columns.

Author(s)

Karline Soetaert

References

Soetaert K, PMJ Herman and JJ Middelburg, 1996a. A model of early diagenetic processes from the shelf to abyssal depths. Geochimica Cosmochimica Acta, 60(6):1019-1040.

Soetaert K, PMJ Herman and JJ Middelburg, 1996b. Dynamic response of deep-sea sediments to seasonal variation: a model. Limnol. Oceanogr. 41(8): 1651-1668.

Examples

# defaults
  CNPDIAsvar()
  head(CNPDIA0D())
  CNPDIA1D()
  
# some runs to work with  
  defsteady <- CNPDIAsolve()
  defdyn    <- CNPDIAdyna()

# altered steady-state run
  out <- CNPDIAsolve(parms = list(Cflux = 1000))

# 0-D outputs
  cbind(steady = CNPDIA0D(defsteady, as.vector = TRUE), 
        dyna = CNPDIA0D(defdyn, as.vector = TRUE), 
        ref = CNPDIA0D(out))

# 1-D outputs
  head(CNPDIA1D(out))
  head(CNPDIA1D(defdyn, which = c("O2", "TOC")))

CNPDIA documentation built on April 15, 2021, 3 p.m.