WRF-CHEM 新手向实践(一):MEIC数据的下载及处理

本文详细描述了如何基于MEIC排放清单处理WRF-CHEM数据,包括数据下载、meic2wrf工具使用、文件名调整以及namelist注意事项。作者分享了实践经验,以帮助其他用户避免初期摸索的困扰。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

背景介绍

基于MEIC排放清单的WRF-CHEM数据处理、运行和一些namelist设置的说明都很少,以至于我在摸索的时候非常痛苦,索性把我自己探索wrf-chem过程中的经验汇总了一下。除了希望能为同志提供一些帮助外,也是想和大家进行更多交流,毕竟本人也是个不知道有没有入门的外行。需说明的是,该系列实践需有一定的WRF基础,仅针对Chem部分详细展开。

MEIC污染源清单(人为排放)

数据下载

数据下载链接:http://meicmodel.org.cn/?page_id=560
在这里插入图片描述
在注册并获取下载资格后,下载最新版(中国污染物排放)即可,为方便后续处理,参考下载选项:在这里插入图片描述

虽然中国污染物排放数据介绍中有CO2,但实际上数据中并不包含,因此还需下载中国碳排放数据(与污染物排放版本号对应):

在这里插入图片描述

解压后的数据如下图所示,是包含了12个月的所有文件(污染物和碳排放),(以一月份为参考案例)选择1月数据放入MEIC文件夹中:

在这里插入图片描述

MEIC数据处理工具及使用

在此使用meic2wrf进行处理,详情可参见此链接:https://cloud.tencent.com/developer/article/1692292,大部分按照其操作手册进行环境配置和运行即可,不过存在一个关键问题需要解决。
由于MEIC数据的更新,部分数据文件名和meic2wrf_GUI.py中的设置不同了,比如PM2.5,PM10等,meic2wrf_GUI.py中的污染物文件名如下所示:

在这里插入图片描述

依据个人经验,如果对代码原理不熟悉的话不建议修改代码,直接修改下载的数据文件名即可,比如将原始文件名中的’PM25’修改为符合代码命名的’PM2.5’。可以不用一个一个对照,可以先运行meic2wrf_GUI.py中的第一步,得到merged文件夹,其中将包含上述污染物种类聚合的结果,每个污染物对应一个文件:

在这里插入图片描述

在报错后只需检查哪个污染物缺失,对照meic2wrf_GUI.py的污染物名称修改原始数据文件的文件名即可。在完成全部污染物聚类后,使用默认剖面,开始meic2wrf_GUI.py的第二步(此处需先按一般WRF运行流程生成所有domain的wrfinput文件),选择对应domain的wrfinput文件,运行即可,每个domain运行一次,都会生成2个对应的wrfchemi文件,各包含12小时的排放结果:

在这里插入图片描述

将wrfchemi排放文件全部复制或者Link到将要运行模拟的文件夹下(即real.exe和wrf.exe所在目录)。至此,MEIC的处理就全部完成了,但后续在设置namelist时仍有几个注意点,可参考此系列第(三)部分的内容,后续将继续更新,有问题也可多交流。

### WRF-Chem 参数化方案与配置方法 WRF-Chem种区域气象-大气化学在线耦合模式,广泛应用于研究气溶胶、臭氧和其他污染物的大气传输及其影响。其参数化方案和配置方法对于实现精确的数值模拟至关重要。 #### 1. 物理过程参数化 WRF-Chem 中涉及的主要物理过程包括微物理过程、辐射过程、陆面模型以及边界层参数化等。这些可以通过 `namelist.input` 文件中的相应选项进行设置: - **微物理过程** 常用的微物理方案包括 Morrison 双矩方案 (Morrison et al., 2009)[^1] 和 Thompson 方案 (Thompson et al., 2008)[^1]。具体选择取决于研究的目标和分辨率。 - **辐射过程** RRTMG(Rapid Radiative Transfer Model for GCMs)是种常用的选择[^2],适用于短波和长波辐射计算。 - **陆面模型** Noah LSM(Land Surface Model)是最常使用的陆面模型,能够较好地描述土壤湿度和温度的变化[^3]。 - **边界层参数化** YSU(Yonsei University Scheme)因其稳定性被推荐用于复杂地形的研究场景[^4]。 #### 2. 大气化学机制 WRF-Chem 支持多种大气化学机制,可通过 `chem_opt` 参数指定。常见的化学机制如下: - **CBMZ/MOSAIC** MOSAIC 提供了种灵活的方式处理多组分气溶胶的动力学行为。它适合于高精度的颗粒物浓度模拟。 - **RADM2/MADE-SORGAM** MADE-SORGAM 主要关注二次有机气溶胶(SOA)形成的过程,而 RADM2 则是个经典的气体相反应机制。 #### 3. 排放数据准备 为了准确反映现实情况下的污染状况,需要准备好详细的排放清单数据。这通常分为人为和自然两部分: - **人为** 包括工业活动、交通运输等活动产生的固定和移动排放。可以利用 EDGAR 或 MEIC 数据集作为输入。 - **自然** 如植被释放的挥发性有机化合物(BVOC),可采用 MEGAN 模型来估算。 #### 4. namelist.input 设置 以下是几个重要的控制参数说明: - `start_year`, `start_month`, `start_day`, `start_hour`: 定义模拟起始时间; - `end_year`, `end_month`, `end_day`, `end_hour`: 定义结束时间; - `interval_seconds`: 输入资料的时间间隔秒数; - `chem_opt`: 选定化学机制编号; - `bio_emis`: 是否开启生物排放开关。 示例代码片段展示如何编辑 `namelist.input` 文件的部分字段: ```fortran &time_control start_year = 2023, start_month = 1, start_day = 1, start_hour = 0, end_year = 2023, end_month = 1, end_day = 7, end_hour = 0, interval_seconds = 21600, / &domains e_we = 100, e_sn = 100, dx = 30000, dy = 30000, grid_id = 1, parent_id = 0, i_parent_start = 1, j_parent_start = 1, parent_grid_ratio= 1, parent_time_step_ratio= 1, max_dom = 1, / &chemistry chem_opt = 200, ! 对应 CBMZ-MOSAIC 组合 bio_emis = .true., / ``` ---
评论 34
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值