ICBMModel(SoilR)
ICBMModel()所属R语言包：SoilR
Implementation of the Introductory Carbon Balance Model (ICBM)
执行的介绍碳平衡模型（ICBM）
译者：生物统计家园网 机器人LoveR
描述----------Description----------
This function is an implementation of the Introductory Carbon Balance Model (ICBM).
是一个实现此功能的介绍性碳平衡模型（ICBM）。
用法----------Usage----------
参数----------Arguments----------
参数：t
A vector containing the points in time where the solution is sougth.
一个矢量点时间的解决方案是sougth的。
参数：ks
A vector of length 2 with the decomposition rates for the young and the old pool.
一个向量长度为2的分解速率为年轻人和老池。
参数：h
Humufication coefficient (transfer rate from young to old pool).
Humufication系数（传输速率从年轻到老池）。
参数：r
External (environmental or edaphic) factor.
外部因素（环境或土壤）。
参数：c0
A vector of length 2 with the initial value of carbon stocks in the young and old pool.
一个向量长度为2的初始值碳储量的年轻人和老年人的游泳池。
参数：In
Mean annual carbon input to the soil.
平均每年的碳输入到土壤中。
（作者）----------Author(s)----------
Carlos A. Sierra <csierra@bgc-jena.mpg.de>, Markus Mueller <mamueller@bgc-jena.mpg.de>
参考文献----------References----------
参见----------See Also----------
TwopSeriesModel
TwopSeriesModel
实例----------Examples----------
#This example reproduces the simulations presented in Table 1 of Andren and Katterer (1997).[此示例再现的Andren，Katterer的（1997）的表1中提出的模拟。]
#First, the model is run for different values of the parameters representing different field experiments. [首先，模型运行的参数为不同的值，代表不同的田间试验。]
times=seq(0,20,by=0.1)
Bare=ICBMModel(t=times) #Bare fallow[休耕]
pNpS=ICBMModel(t=times, h=0.125, r=1, c0=c(0.3,4.11), In=0.19+0.095) #+N +Straw[+ N +秸秆]
mNpS=ICBMModel(t=times, h=0.125, r=1.22, c0=c(0.3, 4.05), In=0.19+0.058) #-N +Straw[-N +秸秆]
mNmS=ICBMModel(t=times, h=0.125, r=1.17, c0=c(0.3, 3.99), In=0.057) #-N -Straw[-N-稻草]
pNmS=ICBMModel(t=times, h=0.125, r=1.07, c0=c(0.3, 4.02), In=0.091) #+N -Straw[+ N-稻草]
FM=ICBMModel(t=times, h=0.250, r=1.10, c0=c(0.3, 3.99), In=0.19+0.082) #Manure[肥料]
SwS=ICBMModel(t=times, h=0.340, r=0.97, c0=c(0.3, 4.14), In=0.19+0.106) #Sewage Sludge[污水处理厂污泥]
SS=ICBMModel(t=times, h=0.125, r=1.00, c0=c(0.25, 4.16), In=0.2) #Steady State[稳态]
#The amount of carbon for each simulation is recovered with the function getC[与函数GETC回收的碳的量为每个模拟]
CtBare=getC(Bare)
CtpNpS=getC(pNpS)
CtmNpS=getC(mNpS)
CtmNmS=getC(mNmS)
CtpNmS=getC(pNmS)
CtFM=getC(FM)
CtSwS=getC(SwS)
CtSS=getC(SS)
#This plot reproduces Figure 1 in Andren and Katterer (1997)[此图再现图1 Andren Katterer（1997）]
plot(times, rowSums(CtBare),type="l",ylim=c(0,8),xlim=c(0,20),ylab="Topsoil carbon mass (kg m-2)",xlab="Time (years)")
lines(times,rowSums(CtpNpS),lty=2)
lines(times,rowSums(CtmNpS),lty=3)
lines(times,rowSums(CtmNmS),lty=4)
lines(times,rowSums(CtpNmS),lwd=2)
lines(times,rowSums(CtFM),lty=2,lwd=2)
lines(times,rowSums(CtSwS),lty=3,lwd=2)
#lines(times,rowSums(CtSS),lty=4,lwd=2)[线（倍，rowSums（CTSS），LTY = 4，LWD = 2）]
legend("topleft",c("Bare fallow", "+N +Straw", "-N +Straw", "-N -Straw", "+N -Straw", "Manure","Sludge"),lty=c(1,2,3,4,1,2,3),lwd=c(1,1,1,1,2,2,2),bty="n")
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-9-30 11:34:20
