Century 模型 R代码分享
2019-8-30 13:45
阅读:8893
Century 是传统的土壤碳过程模型,这里分享本人最近找到的三个相关版本:
1. Fortran 版本(linux)
该版本的下载地址:
https://github.com/bsulman/INTERFACE-model-experiment-synthesis
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 |
即可得到土壤碳变化的输出结果文件:
2. R版本 by Shoji Hashimoto
可以在https://github.com/LukeEcomod/SoilCarbon 下载到
这个文件直接把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 条评论