RothCModel(SoilR)
RothCModel()所属R语言包:SoilR
Implementation of the RothCModel
实施的RothCModel
译者:生物统计家园网 机器人LoveR
描述----------Description----------
This function implements the RothC model of Jenkinson et al. It is a wrapper for the more general function GeneralModel.
此功能实现了RothC模型的詹金森等。这是一个更普遍的功能GeneralModel的包装。
用法----------Usage----------
k.IOM = 0), C0 = c(0, 0, 0, 0, 2.7), In = 1.7, DR = 1.44,
参数----------Arguments----------
参数:t
A vector containing the points in time where the solution is sought.
一种向量,包含的点的时间寻求溶液。
参数:ks
A vector of lenght 5 containing the values of the decomposition rates for the different pools
一个向量的长度包含的值的分解速率不同的池
参数:C0
A vector of length 5 containing the initial amount of carbon for the 5 pools.
长度为5的一种向量,含有碳的初始量为5个泳池。
参数:In
A scalar or data.frame object specifying the amount of litter inputs by time.
一个标量或数据框对象,指定的凋落物输入量。
参数:DR
A scalar representing the ratio of decomposable plant material to resistant plant material (DPM/RPM).
标量的比例抗性的植物材料分解的植物材料(DPM / RPM)。
参数:clay
Percent clay in mineral soil.
矿质土壤的粘土的百分比。
参数:xi
A scalar or data.frame object specifying the external (environmental and/or edaphic) effects on decomposition rates.
一个标量或数据框对象指定的外部(环境和/或土壤)的分解率的影响。
参数:solver
A function that solves the system of ODEs. This can be euler or ode or any other user provided function with the same interface.
的功能,解决了系统的微分方程。这可能是euler或ode或任何其他用户提供的函数具有相同的接口。
值----------Value----------
A Model Object that can be further queried
一个模型对象,可以进一步查询
(作者)----------Author(s)----------
Carlos A. Sierra <csierra@bgc-jena.mpg.de>, Markus Mueller <mamueller@bgc-jena.mpg.de>
参考文献----------References----------
参见----------See Also----------
ICBMModel
ICBMModel
实例----------Examples----------
t=0:500
Ex=RothCModel(t)
Ct=getC(Ex)
Rt=getReleaseFlux(Ex)
plot(t,Ct[,1],type="l",col=1, ylim=c(0,25),ylab=expression(paste("Carbon stores (Mg C", ha^-1,")")),xlab="Time (years)",lwd=2)
lines(t,Ct[,2],col=2,lwd=2,lty=2)
lines(t,Ct[,3],col=3,lwd=2,lty=3)
lines(t,Ct[,4],col=4,lwd=2,lty=4)
lines(t,Ct[,5],col=5,lwd=2,lty=5)
lines(t,rowSums(Ct),lwd=2)
legend("topright",c("Pool 1, DPM", "Pool 2, RPM", "Pool 3, BIO","Pool 4, HUM","Pool 5, IOM","Total Carbon"),lty=c(1:5,1),lwd=rep(2,5),col=c(1,2,3,4,5,"black"),bty="n")
plot(t,Rt[,1],type="l",ylim=c(0,2),ylab="Respiration (Mg C ha-1 yr-1)",xlab="Time")
lines(t,Rt[,2],col=2)
lines(t,Rt[,3],col=3)
lines(t,Rt[,4],col=4)
lines(t,Rt[,5],col=5)
lines(t,rowSums(Rt),lwd=2)
legend("topright",c("Pool 1, DPM", "Pool 2, RPM", "Pool 3, BIO","Pool 4, HUM","Pool 5, IOM","Total Respiration"),lty=c(1,1,1,1,1,1),lwd=c(1,1,1,1,1,2),col=c(1,2,3,4,5,1),bty="n")
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|