WRF-CHEM 实战指南（一）：MEIC清单的获取、适配与namelist关键配置

1. MEIC排放清单WRF-CHEM模拟的起点刚开始接触WRF-CHEM时最让我头疼的就是找不到合适的人为排放源数据。国内研究常用的MEIC清单Multi-resolution Emission Inventory for China虽然权威但实际操作中会遇到各种坑。记得第一次下载数据时光是搞明白哪些文件有用就花了半天时间。MEIC官网提供了中国高分辨率人为源排放数据覆盖了SO2、NOx、PM2.5等主要污染物。最新版本的数据按月份分类存储每个压缩包包含12个月的数据文件。这里有个实用建议下载时一定要选择与你的模拟时段匹配的月份比如研究冬季污染就下载1月数据。我遇到过有人下载了全年数据最后只用到一个月份既浪费存储空间又增加处理复杂度。下载后的文件结构通常是这样MEIC_2020/ ├── 01 │ ├── CO_202001.nc │ ├── NOx_202001.nc │ └── ... ├── 02 │ ├── CO_202002.nc │ └── ... └── ...2. 数据下载与预处理避坑指南MEIC官网的注册流程比较简单但有几个细节需要注意。注册时需要提供机构邮箱个人邮箱可能会被拒绝。我帮学生处理过这种情况最后只能用实验室邮箱重新申请。下载权限通常会在1-2个工作日内开通建议提前规划好时间。下载数据时会遇到两个选择中国污染物排放和中国碳排放。这里有个容易忽略的点虽然污染物排放数据介绍中提到包含CO2但实际上需要单独下载碳排放数据。我第一次模拟时就栽在这个坑里导致后续处理时总是报错缺少CO2数据。解压后的文件命名也可能带来麻烦。最新版MEIC数据使用了PM25这样的命名但常用的meic2wrf工具仍要求PM2.5这样的旧格式。我的经验是不要直接修改工具代码而是重命名数据文件更稳妥。可以用这个bash命令批量修改for file in *PM25*; do mv $file ${file/PM25/PM2.5}; done3. meic2wrf工具实战操作meic2wrf是目前最常用的MEIC数据转换工具它能将NetCDF格式的MEIC数据转换为WRF-CHEM需要的wrfchemi文件。安装过程比较简单主要依赖Python3和netCDF4库。建议用conda创建独立环境conda create -n meic2wrf python3.8 conda activate meic2wrf pip install netCDF4 numpy工具运行时最容易出问题的就是污染物种类匹配。这里分享我的调试方法先运行第一步生成merged文件夹检查里面的文件是否齐全。如果缺少某个污染物工具会报错并提示缺失的文件名。这时只需要回到原始数据找到对应的文件并按照要求重命名即可。处理完所有污染物后第二步需要提供对应domain的wrfinput文件。这里要注意domain编号的对应关系特别是当使用嵌套网格时。我建议在处理前先用ncdump检查wrfinput的维度信息ncdump -h wrfinput_d014. namelist.chem关键配置解析完成wrfchemi文件生成后还需要正确配置namelist.chem才能让WRF-CHEM正常读取排放数据。以下几个参数需要特别注意emiss_opt: 这个参数决定了排放源的种类使用MEIC数据时通常设置为2 anthropogenic emissionschem_opt: 根据你的化学机制选择比如202对应RADM2机制emiss_inpt_opt: 对于MEIC数据这个参数应该设置为1读取外部排放文件我遇到过最棘手的问题是时间插值设置。MEIC数据是月均值而WRF-CHEM需要小时数据。正确的做法是在namelist中设置emiss_cycle_yr 1, emiss_cycle_dy 1, emiss_cycle_hr 125. 实战中的常见问题排查即使按照流程操作新手还是会遇到各种问题。最常见的是文件路径错误特别是当使用符号链接时。建议先用绝对路径测试确认无误后再改用相对路径。可以用这个命令检查文件链接ls -l wrfchemi_*另一个常见问题是内存不足。处理高分辨率MEIC数据时meic2wrf可能会消耗大量内存。我的经验是单次处理不要超过3个domain或者分批次处理。如果遇到内存错误可以尝试增加swap空间sudo fallocate -l 4G /swapfile sudo chmod 600 /swapfile sudo mkswap /swapfile sudo swapon /swapfile最后提醒一点不同版本的WRF-CHEM对排放文件格式要求可能不同。如果遇到莫名其妙的报错建议检查WRF版本是否与meic2wrf工具兼容。我在WRF4.3上就遇到过因为时间变量格式不兼容导致的运行错误。