210: Plots the output of 'getC' or 'getReleaseFlux' for each pool...

Description Usage Arguments Author(s) Examples

Description

This function produces a plot with the C content or released C for each pool over time. Needs as input a matrix obtained after a call to getC or getReleaseFlux.

Usage

1
plotCPool(t, mat, col, ...)

Arguments

t

A vector containing the time points for plotting.

mat

A matrix object obtained after a call to getC or getReleaseFlux

col

A color palette specifying color lines for each pool (columns of mat).

...

Other arguments passed to link{plot}.

Author(s)

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

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
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)
Rt=getReleaseFlux(Ex)

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

Example output

Loading required package: deSolve
Loading required package: parallel
Loading required package: RUnit
[1] "format= Delta14C"
[1] "format= Delta14C"

SoilR documentation built on May 1, 2019, 8:06 p.m.