庄伟
Century 模型 R代码分享
2019-8-30 13:45
阅读:8089

Century 是传统的土壤碳过程模型,这里分享本人最近找到的三个相关版本:


1. Fortran 版本(linux)


该版本的下载地址:

https://github.com/bsulman/INTERFACE-model-experiment-synthesis


DAYCENT.zip

Sulman et al (2018)发表了多个土壤碳模型比较的文章,并分享了部分模型的代码,century的源代码也在其中。其中source文件的代码,通过makefile命令,可以生成DailyDayCent_srad_SOM文件,这个文件是在linux中的可执行文件。将DailyDayCent_srad_SOM文件复制到【RUN】文件夹中。

Linux系统(ubuntu或者WIN10中的ubuntu子系统)中,cd 到DailyDayCent_srad_SOM所在的文件夹【RUN】执行如下命令:

./DailyDayCent_srad_SOM  -s exp_HQ500  -n exp_HQ500.high_protect -c experiment_25C.cdi.txt

即可得到土壤碳变化的输出结果文件:

TIM截图20190830132705.png


2. R版本  by Shoji Hashimoto


可以在https://github.com/LukeEcomod/SoilCarbon 下载到


CENTURY_supplement.BG.R


这个文件直接把Century从Fortran翻译成R,可以直接在R里面执行,比较Fortran容易看懂,初次接触,可以先看看这个版本,容易上手。


   3. SOILR版本  


可以在https://github.com/MPIBGC-TEE/SoilR-exp 下载到。

其中的CENTURY模型是简化版本,可以做参数优化,感兴趣的同学可以尝试。



转载本文请联系原作者获取授权,同时请注明本文来自庄伟科学网博客。

链接地址:https://wap.sciencenet.cn/blog-526092-1195878.html?mobile=1

收藏

分享到:

当前推荐数:0
推荐到博客首页
网友评论0 条评论
确定删除指定的回复吗?
确定删除本博文吗?