pyDICE:DICE模型的python版本

Posted by CryoECNU on June 18, 2022

初步根据DICE模型的GAMS代码,翻译了一版Python的代码,主要用了Pyomo包+ipopt数学规划求解器。

源GAMS代码可以从网上下载:http://www.econ.yale.edu/~nordhaus/homepage/homepage/DICE2016R-091916ap.gms

参考文献:”Revisiting the social cost of carbon: Estimates from the DICE-2016R model”, https://www.pnas.org/doi/10.1073/pnas.1609244114

验证了GAMS代码输出的33个变量,挑几个放这里展示一下。

基线运行

优化解

Pyomo这个优化包很人性化,与GAMS和AMPL之类的软件风格相似,会GAMS的很容易上手将GAMS代码改写成python代码。pyomo也提供直接转写成其他语言的源程序【https://pyomo.readthedocs.io/en/stable/library_reference/solvers/gams.html?highlight=write#gams-writer】,试了一下GAMS的转写结果,可以用GAMS运行。不过暂时没体会(发现)这功能有啥特别用处,待挖掘。

DICE模型的区域版——RICE模型——正在路上……