WRF-CHEM 新手向实践（三）：生物质燃烧排放及化学初始/边界场

头发凌乱的在读博士

已于 2024-08-22 19:07:34 修改

阅读量943

点赞数 17

分类专栏： WRF-CHEM 文章标签：经验分享学习笔记

于 2024-08-22 19:00:31 首次发布

本文链接：https://blog.csdn.net/qq_42335815/article/details/135902989

版权

WRF-CHEM 专栏收录该内容

3 篇文章 6 订阅

订阅专栏

这几个月实在太忙，但因为有不少同学询问，抽空更新一章，最后的第四个部分也会尽快完成。

在WRF-CHEM 新手向实践（二）：生物排放中已完成了Biogenic Emissions的处理，接下来介绍：（1）Biomass Burning Emissions即生物质燃烧排放的处理，（2）使用mozbc处理CAM-Chem以获得化学初始和边界场。

1.生物质燃烧排放的处理

资源下载：FINN数据及其处理工具，下载链接：FINN: Data

具体的，首先填写自己的信息：

接下来可以选择下载不同版本的FINN数据，较新的v2.5数据需要转到另一个链接进行下载：

本例选择下载1.5版本的FINN数据，如果自己的案例是2020年之前的都可以选择这个老版本，选择相应的时间和化学物种（需与后续选择的化学方案相对应）：

在勾选FINNv1.5数据后点击网页最下方的Submit，将获得数据下载链接。随后再勾选FINN处理工具，然后Submit，将获得数据处理工具的下载链接（注意该网页一次只能勾选一个）：

FINN工具编译及运行：

FINN工具的压缩包中包含run和src两个文件夹，其中run文件夹中是几个inp文件的示例，src是FINN工具源码位置，在此文件夹下进行程序编译。首先，执行命令设置fortran编译器，在这我使用的是gfortran，随后执行make命令：

export FC=gfortran
./make_fire_emiss

运行成功的话，会得到fire_emis可执行程序。

FINN排放数据处理：

对于WRF-CHEM模拟，从run文件中复制example_finn15_wrf.inp文件至src下并修改，具体如下：

fire_directory和fire_filename(1)分别是之前下载的FINN数据存放的路径和文件名，wrf_directory是wrfinput文件的所在路径，start_date和end_date是所需要的数据时间，domains是模拟域的数量， FinnVers是所使用FINN数据的版本，其余保持默认。

这部分内容包括编译部分都可以参考src文件夹中的README_GRID_FINN_v20200609.pdf，还是比较详细的。

在修改完成后，执行以下命令：

./fire_emis < example_finn15_wrf.inp > fire_emis.out

如果成功，会在src文件夹下生成每个模拟域规定时间范围内的每小时的wrffirechemi文件：

至此，生物质排放输入文件就准备好了。

2.初始/边界场的生成

CAM-CHEM数据下载：CAM-Chem Download。基本流程和前面都差不多，填写信息然后提交，在邮箱中接收数据链接。

mozbc工具下载：WRF-Chem，一样的流程，不赘述：

mozbc工具编译

解压mozbc压缩包，修改其中的Makefile文件，在设置NETCDF环境变量的情况下，只需添加以下命令：

FC = gfortran
AR_FILES = -lnetcdff -lnetcdf

此处设置类似于WRF-CHEM 新手向实践（二）：生物排放，可参考。随后，执行make_mozbc命令，得到mozbc可执行程序。

mozbc运行

需修改对应的inp文件，在mozbc的文件夹中包含了几个示例的inp文件，比如我们在WRF-CHEM中使用CBMZ-MOSAIC_8bins的话，就可以选择CBMZ-MOSAIC_8bins.inp文件，但需注意的是，对于CAM数据，无法直接使用CBMZ-MOSAIC_8bins.inp示例文件中的物种，需要映射CAM至WRF-CHEM气相物种和气溶胶模型。

暂时不看这部分，先修改上半部分的设置，在CBMZ-MOSAIC_8bins.inp文件中：

&control
do_bc = .true.
do_ic = .true.
domain = 1
dir_wrf = '/glade/scratch/pfister/WRFreal_WACCM/'
dir_moz = './'
fn_moz = 'output_WACMM_0001.nc'
moz_var_suffix = ''
def_missing_var = .true.

其中的domain并非是模拟域数量，而是模拟域代号，1就代表最外层，在运行嵌套模拟时，从最外层到最内层，每次修改domain后都需重复运行一次mozbc，有几个wrfinput文件就运行几次；dir_wrf是wrfinput文件所在位置， dir_moz是之前下载的CAM数据所在目录，fn_moz是CAM数据的名称，其与参数保持默认。随后修改spc_map参数，可参考官方的CESM-WRFchem_aerosols_20190822.pdf文档：

spc_map = 'o3 -> O3', 'n2o -> N2O', 'no -> NO',
'no2 -> NO2', 'nh3 -> NH3', 'hno3 -> HNO3', 'hno4 -> HO2NO2',
'n2o5 -> N2O5', 'h2o2 -> H2O2',
'ch4 -> CH4', 'co -> CO', 'ch3ooh -> CH3OOH',
'hcho -> CH2O', 'ch3oh -> CH3OH', 'c2h4 -> C2H4',
'ald -> CH3CHO', 'acet -> CH3COCH3', 'mgly -> CH3COCHO',
'pan -> PAN', 'mpan -> MPAN', 'macr -> MACR',
'mvk -> MVK', 'c2h6 -> C2H6', 'c3h6 -> C3H6', 'c3h8 -> C3H8',
'c2h5oh -> C2H5OH', 'c10h16 -> MTERP',
'isopr -> ISOP','acetol -> HYAC', 'mek -> MEK',
'bigene -> BIGENE', 'bigalk -> BIGALK',
'tol -> TOLUENE', 'benzene -> BENZENE', 'xylenes -> XYLENES',
'cres -> CRESOL', 'dms -> DMS', 'so2 -> SO2',

再加上气溶胶部分，选择与我们的CBMZ-MOSAIC_8bins化学方案对应的气溶胶模型：