WRF-CHEM 新手向实践（四）：化学初始/边界场及namelist设置

这是这个系列的最后一节了，这些经验是自己作为一个半路出家的外行一步步摸索出来的，其中不免有很多理解不到位的地方，但也正是如此，我相信这是目前新手入门wrf-chem最友好的教程之一，希望能对大家有实质性的帮助。

本节继续WRF-CHEM 新手向实践（三）：生物质燃烧排放及化学初始/边界场未完成的内容。在mozbc生成化学初始和边界条件之前，还需要先在wrfinput文件中开启化学部分。

namelist.input设置

主要涉及到 &time_control和 &chem 这两部分的修改。这里只介绍各排放源输入的必要设置，更具体的namelist选项设置可自行查询用户手册。

针对MEIC人为排放的输入，添加以下设置：

&time_control
auxinput5_inname                    = 'wrfchemi_<hour>z_d<domain>',
auxinput5_interval_m                = 60, 60, 60,
frames_per_auxinput5                = 12, 12, 12,
io_form_auxinput5                   = 2,
 
&chem
chem_opt                            = 10,
io_style_emissions                  = 1,
emiss_opt                           = 4,         4,         4,
emiss_inpt_opt                      = 101,      101,       101,

auxinput5_inname是我们在第（一）节中生成的12小时人为排放文件名，对于d01即为wrfchemi_00z_d01和wrfchemi_12z_d01，每个domain两个；

io_form_auxinput5 = 2 以WRF Netcdf格式开启人为排放输入；

auxinput5_interval_m是数据更新的时间间隔，单位是分钟，对于我们使用的wrfchemi_<hour>z_d<domain>文件，间隔为每小时也就是60分钟，可按自己的需求整数倍调整；

frames_per_auxinput5 = 12表示每个wrfchemi_<hour>z_d<domain>文件中包含12个子时刻。

对于&chem部分，chem_opt  = 10表示开启化学部分，并选择选项为10的CBMZ chemical mechanism。此前的wrfinput都在chem_opt  = 0的情况下生成，只包含气象场，在这一节中，我们需要将化学部分开启。

io_style_emissions = 1 这个设置非常重要，这代表我们使用的是12小时的循环文件，emiss_opt 代表人为排放的化学机制，需要与chem_opt = 10相对应，使用的是CBMZ/MOSAIC机制（选项4）。然而，我们在第（一）节中处理的MEIC人为排放实际上是 RADM2机制的。这就需要使用emiss_inpt_opt  = 101进行转换，其作用是将RADM2排放数据文件进行调整，以遵循CBMZ/MOSAIC框架。

如果使用其他形式（非12小时循环）的人为排放文件，比如特定日期/时间排放数据文件，需设置

auxinput5_inname = wrfchemi_d<domain>_<date>,
io_style_emissions = 2, 
frames_per_auxinput5 = 1

针对Megan生物排放的输入，添加以下设置：

&time_control
auxinput6_inname                    = 'wrfbiochemi_d<domain>',
auxinput6_interval_m                = 5760, 5760, 5760,
io_form_auxinput6                   = 2,
 
&chem
chem_opt                            = 10,
bio_emiss_opt                       = 3,         3,         3,
ne_area                             = 41，

auxinput6_inname 是第（二）节中处理的Megan生物排放文件名；

auxinput6_interval_m 是数据输入更新的时间间隔，由于Megan生物排放输入并没有时间变化，所以可以直接设置一个较长的间隔；

bio_emiss_opt = 3 生物排放选项3即使用MEGAN生物排放；

ne_area = 41，MEGAN生物排放使用的物种总数，这个值也可以越大越好。

针对FINN生物质燃烧排放的输入，添加以下设置：

&time_control
auxinput7_inname                    = 'wrffirechemi_d<domain>_<date>',
auxinput7_interval_m                = 60, 60, 60,
io_form_auxinput7                   = 2,
 
&chem
chem_opt                            = 10,
biomass_burn_opt                    = 3,         3,         3,
scale_fire_emiss                    = .false.,    .false.,   .false.,

auxinput7_inname 是第（三）节中处理的FINN生物质燃烧排放文件名；

auxinput7_interval_m 是数据输入更新的时间间隔，FINN可以做到每小时更新，可根据需要按整数倍调整；

biomass_burn_opt = 3 该选项代表使用MOZART的生物质燃烧排，这在FINN数据下载时确定；

scale_fire_emiss 表示是否要将FINN数据分解为焖烧和燃烧两个部分，只对MOZART or MOZCART机制有效，如果使用这两个chem_opt的话，则必须设置为true。

完成namelist.input设置后，确保所有排放文件都在运行目录下，随后运行real.exe，生成包含化学变量的wrfinpu文件。需观察自己的wrfinput文件内容是否发生变化，若无变化，需在bashrc中添加export WRF_CHEM=1，随后重新编译WRF-CHEM。

化学初始/边界场生成

在mozbc可执行文件目录下，执行以下命令：

./mozbc < CBMZ-MOSAIC_8bins.inp > OUT.log

有几个domain就运行几次，从外层向内层运行，每次都需要修改CBMZ-MOSAIC_8bins.inp中的domain编号。mozbc的编译和配置参考第（三）节内容。完成后，可得到包含化学初始场和边界场的wrfinput文件，随后就可以开始执行WRF程序了。

&chem部分还有很多自定义设置就需要大家自己去了解和使用了。至此，WRF-CHEM入门系列就完结了，有问题或者建议都可以向我反映。谢谢大家，祝学业有成。