R Graphical Manual

Browse All

Last data update: 2014.03.03

R: Add lines with the output of 'getC14', 'getC', or...

linesCPool

R Documentation

Add lines with the output of `getC14`, `getC`, or `getReleaseFlux` to an existing plot

Description

This function adds lines to a plot with the C content, the C release, or Delta 14C value of each pool over time. Needs as input a matrix obtained after a call to getC14, getC, or getReleaseFlux.

Usage

linesCPool(t, mat, col, ...)

Arguments

`t`	A vector containing the time points for plotting.
`mat`	A matrix object obtained after a call to `getC14`, `getC`, or `getReleaseFlux`.
`col`	A color palette specifying color lines for each pool (columns of `mat`).
`...`	Other arguments passed to `plot`.

Author(s)

Carlos A. Sierra <csierra@bgc-jena.mpg.de>, Markus Mueller <mamueller@bgc-jena.mpg.de>

Examples

years=seq(1901,2009,by=0.5)
LitterInput=700 

Ex=ThreepFeedbackModel14(t=years,ks=c(k1=1/2.8, k2=1/35, k3=1/100),
                         C0=c(200,5000,500),F0=c(0,0,0), In=LitterInput, 
                         a21=0.1,a12=0.01,a32=0.005,a23=0.001,inputFc=C14Atm_NH)
Ct=getC(Ex)

pal=rainbow(3)
plotCPool(t=years,mat=Ct,col=pal,xlab="Time (yrs)",
          ylab="Carbon stocks",ylim=c(min(Ct),max(Ct)))

LitterInput2=350 

Ex2=ThreepFeedbackModel14(t=years,ks=c(k1=1/2.8, k2=1/35, k3=1/100),
                          C0=c(200,5000,500),F0=c(0,0,0), In=LitterInput2, 
                          a21=0.1,a12=0.01,a32=0.005,a23=0.001,inputFc=C14Atm_NH)
Ct2=getC(Ex2)

linesCPool(t=years,mat=Ct2,col=pal,lwd=2)