在WRF-CHEM 新手向实践(一)中已完成了MEIC人为排放源的处理,根据官方的教程,还应输入Biogenic Emissions(生物排放),Biomass Burning Emissions(生物质燃烧排放),Dust Emissions(粉尘排放)、GOCART Background Data(背景场数据)volcanic emissions(火山灰排放)等等,但根据我的理解,这些都是可选项,本节仅介绍生物排放和生物质燃烧排放这两个最经常使用的部分。(注:MEIC中有粉尘了,一般不需要额外Dust Emissions,GOCART Background Data只在特定化学机制下使用chem_opt = 300,301,火山排放一般情况下都不使用)
1. 生物排放采用Megan及其处理程序
数据及工具下载链接:WRF-Chem
进入网页后,填写自己的信息,输入要下载的数据范围,建议直接填写最大范围,因为数据量不大:
填写完成后,分别点击bio_emiss和bio_emiss_input files按钮下载megan处理工具及数据 (其中数据需要等待一段时间后在邮箱中确认下载链接):
下载完成并解压后,开始megan_bio_emiss的编译
在megan_bio_emiss文件夹中的Makefile中修改如下命令:
FC = gfortran
AR_LIBS = -lnetcdff -lnetcdf
其中FC为你的Fortran编译器类型,此外,如果没有提前将netcdf路径设置到环境变量中,此处还要额外加入一行:
NETCDF_DIR = 你的netcdf路径
完成后在该路径下执行make命令(即在终端中输入make),我在该过程遇到了gfortran报错,是因为高版本不兼容问题,因此需要使用较低的gfortran版本(如不存在这个问题可跳过)。可输入以下命令查看自己的编译器版本:
gfortran --version
我之前用的是gfortran-11,之后额外安装了gfortran-9用来编译megan和FINN,不用卸载之前的11版本,因为WRF4.3的编译仍然需要gfortran-11,我们可以使用如下命令在多个gfortran版本中来回切换:
sudo update-alternatives --config gfortran
切换至gfortran-9后再次make,成功后生成megan_bio_emiss可执行文件。
bio_emiss数据处理
在megan_bio_emiss文件夹中创建megan_bio_emiss.inp文件,其中内容的格式可参照README.bio_emiss文件中的介绍,相关参数的含义也在其中:
&control
domains = 3,
start_lai_mnth = 5,
end_lai_mnth = 6,
wrf_dir = '/home/stacy/megan/wrf_files',
megan_dir = '/home/stacy/megan/30sec'
/
因为是新手向,这里再额外多嘴两句:domain是你WRF的domain数量,start_lai_mnth代表处理的起始月份,另一个是终止月份,需要注意的是,除了要模拟时段外,megan选项还额外要求模拟时段前一个月的输入,举例来说,如果要模拟6月份,则start_lai_mnth=5,但是要模拟1月份,则必须要输入12个月的结果,即:
start_lai_month = 1,
end_lai_month = 12,
此外,wrf_dir是wrf_input文件所在的目录(在不开启chem选项的情况下事先正常运行real.exe生成),megan_dir是之前下载的megan.data的目录。这些内容在README.bio_emiss中都有更详细的介绍。
完成megan_bio_emiss.inp的创建后,就可以运行megan_bio_emiss生成生物质排放文件了,执行以下命令:
./megan_bio_emiss < megan_bio_emiss.inp > megan_bio_emiss.out
如果顺利运行,会生成与wrfinput文件数量对应的wrfbiochemi文件:
到这megan生物质排放处理就做完了,其实大部分内容都可以参考README.bio_emiss,希望这篇博客能更好地帮助理解。如何设置namelist.input加入生物质排放将在后续内容中统一说明 。下一节将继续介绍FINN生物质燃烧排放的处理。