模式利器 | MEIC污染源清单向WRF-Chem模式网格插值分配工具——meic2wrf

2020-09-04 11:08:20 浏览数 (1)

开源工具地址:

https://gitee.com/jinfan0931/meic2wrf (推荐国内仓库√)

https://github.com/jinfan0931/meic2wrf

欢迎大家star!

工具简介

本工具针对清华大学MEIC源排放清单[1](0.25°×0.25°,2016),实现了清单中各污染物排放量向WRF-Chem模式网格的插值和分配,可窗口界面运行或脚本运行。

获取地址(已从Github搬迁至Gitee,Gitee为国内代码托管平台,稳定,下载速度快):https://gitee.com/jinfan0931/meic2wrf.git[2](求star,欢迎pull request共同改进?)

使用方法请参考Gitee该项目仓库中的《操作手册》。

注意: 本工具要求WRF-Chem用户设定的最大模拟范围处于MEIC清单覆盖范围内。

基本原理

以对2016年MEIC清单1月份NH3排放量的处理为例:

  1. 提取MEIC清单中某一部门(如农业部门)NH3排放量,将排放量单位转换为WRF-Chem模式要求的单位后,插值到模式模拟区域
  1. 按照部门排放特征,确定逐小时和逐层排放因子后,将模拟区域的NH3排放量进行时空分配:
  1. 对其余四部门(工业、电力、民用、交通)的NH3排放量重复1、2步骤后,再将所有部门NH3排放量同时刻、同层进行加和,即完成将MEIC清单NH3排放量插值并分配到WRF-Chem模拟区域网格上。对MEIC清单中其余物种的处理过程同NH3。
运行环境

Python 3.7.7[3]

PyNIO 1.5.5[4]

可通过Conda——环境管理及包安装工具,搭建上述运行环境。

窗口界面运行——meic2wrf_GUI.py

搭建好运行环境后,通过命令python meic2wrf_GUI.py启动窗口界面。工具的主界面由1. 整合MEIC清单中各污染物不同部门的排放数据2. 向WRF-Chem模式网格插值分配”两部分组成。 如图所示:

窗口界面运行方法请参考Gitee该项目仓库中的《操作手册》[5]

终端脚本运行——meic2wrf_noGUI.py

若以脚本模式运行,终端下用文本编辑器打开meic2wrf_noGUI.py(如vim meic2wrf_noGUI.py),在脚本末尾if __name__ == '__main__':后,指定相应文件和目录路径,保存退出后,终端输入python meic2wrf_noGUI.py,等待脚本运行结束,则会在save_dir路径下生成WRF-Chem源排放文件。

WRF-Chem多层嵌套网格

若用户设定的WRF-Chem网格为多层嵌套,即有多个wrfinput文件(如对于三层嵌套,有:wrfinput_d01,wrfinput_d02,wrfinput_d03三个初始场文件),只需对每个wrfinput文件重复上述操作,即可得到对应WRF-Chem嵌套网格的MEIC清单插值分配结果。

其它说明

引用: 樊晋,周永龙,徐旋烨,蒋沛雅,李卓.MEIC污染源清单向WRF-Chem模式网格插值分配程序软件V1.0[Z].国家版权局,2020SR0414107.

感谢中国气象局气象科学研究院的张磊博士在工具开发期间提供的帮助。

感谢Github开源平台吕浩工程师对meic2wrf_noGUI.py脚本的贡献。

感谢《气象学家》公众号的支持和推广。

感谢成都信息工程大学大创项目的支持(项目编号:202010621016)。

感谢成都兰蓝环境科技有限公司的协助。

参考资料

[1]

清华大学MEIC源排放清单: http://www.meicmodel.org/

[2]

工具地址: https://gitee.com/jinfan0931/meic2wrf.git

[3]

Python 3.7.7: https://www.python.org/downloads/release/python-377/

[4]

PyNIO 1.5.5: https://www.python.org/downloads/release/python-377/

[5]

《操作手册》: https://gitee.com/jinfan0931/meic2wrf.git

0 人点赞