【发布时间】：2022-12-04 14:35:34
【问题描述】：
I am trying to run the RothC model using the SoilR package. One of the steps involves the estimation of the effects of moisture on decomposition rates.
According to the model description (Page 14), the parameter b is calculated as:
if acc. TSMD < 0.444 max.TSMD,
b = 1.0
otherwise,b = 0.2 + (0.8) * (max.TSMD - acc.TSMD)/(max.TSMD - 0.444 max.TSMD)
but the soilR implementation, fW.RothC.R, is:
b=ifelse(Acc.TSMD > 0.444*Max.TSMD, 1, (0.2+0.8*((Max.TSMD-Acc.TSMD)/(Max.TSMD-0.444*Max.TSMD))))
As you can see, the condition has > instead of <
How to explain this?
I don\'t have enough reputation to create a tag, but it would be useful to add the tags soilr and rothc.