这是这个系列的最后一节了,这些经验是自己作为一个半路出家的外行一步步摸索出来的,其中不免有很多理解不到位的地方,但也正是如此,我相信这是目前新手入门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入门系列就完结了,有问题或者建议都可以向我反映。谢谢大家,祝学业有成。