WRF /化学:编辑NetCDF (wrfchemi文件)的排放

Question

这里有人有经验的编辑wrfchemi输出prep_chem效用,这样我覆盖文件用自己的排放?

我在python中知识渊博,但我不是那么好面向NetCDF文件的操作。

我想做的,是用我自己的排放在WRF /化学模型

1. 我有一个shapefile排放(网格)的每个网格都有排放值。
2. 我要将其转换为可读的格式WRF /化学;然而,我认为编辑wrfchemi文件由prep_chem容易吗?(告诉我如果不是)

我读过WRF /化学排放指南,但它并没有阐述这尤其是关于下面的详细信息

• 如何处理物种形成(我的shapefile只有PM10和PM2.5)。哪里我把这个新的排放文件(编辑wrfchemi)我想要的。有许多点物种(如pm2.5成核模式,点,等等)。我把这些不确定权重。
• 我的排放不与时间有关的,我做什么当我编辑wrfchemi文件。(wrfchemi文件我现在也没有时间相关)

希望有人可以帮忙。

谢谢!

Answer 1

你可以试试R包eixport,wrf_put

#示例#读取数组排放,有限公司<——wrf_get(文件=“Path_to_WRFCHEMI”,名称=“E_CO”) #改变的值,这里你应该使用你的数据有限公司[]= rnorm(长度(CO)) # Inyect排放到wrfchemi wrf_put(文件=“Path_to_WRFCHEMI name = " E_CO "波尔=有限公司)

如何处理物种形成(我的shapefile只有PM10和PM2.5)。哪里我把这个新的排放文件(编辑wrfchemi)我想要的。有许多点物种(如pm2.5成核模式,点,等等)。我把这些不确定权重。

你需要知道你当地的物种形成。大学在我的部门。圣保罗,我们使用下面的物种形成的PM 2.5 (g / h /公里^ 2)

e_so4i = 0.0077, e_so4j = 0.0623, e_no3i = 0.00247, e_no3j = 0.01053, e_pm2.5i = 0.1, e_pm2.5j = 0.3, e_orgi = 0.0304, e_orgj = 0.1296, e_eci = 0.056, e_ecj = 0.024,水= 0.277

检查这

我的排放不与时间有关的,我做什么当我编辑wrfchemi文件。(wrfchemi文件我现在也没有时间相关)

首先需要编辑名称列表。wps和名称列表。输入所需的时间长度。然后,在运行。/真实的。exe wrfinput_d0x你会。在这个阶段,您可以使用wrf_create创建一个wrfchemi文件和0。阅读手册,看看你想要两个0-12z 12-0z文件或一个文件的长度。然后您可以使用wrf_put。只是效仿和阅读手册。

你也可以尝试R包EmissV

引用:

1. Ibarra-Espinosa et al ., (2018)。eixport: R包出口排放大气模型。开源软件学报,3(24),607年,https://doi.org/10.21105/joss.00607Schuch et al ., (2018)。
2. EmissV: R包来创建车辆和其他排放空气质量模型。开源软件学报,3(30),662年,https://doi.org/10.21105/joss.00662

Answer 2

作为处理netCDF文件而言,您可以使用这个例子。我所做的在过去通过维度,循环变量,和文件属性,if语句更改的部分我想要编辑的文件。

在你的情况中,你可以做你想做的是最好的。或许你可以把平均排放在次网格平滑wrfchemi网格和shapefile之间变化。只要确保单位是一致的。关于时间变化,你可以做你认为是合适的,指定你做什么当你沟通方法。如果它对你来说更容易输入常数每小时平均每天排放的排放,那么这样做。如果你认为更容易保持满满当当的模式,那么你可以这样做。

物种形成而言(见页14日至15日),你可以把PM10 e_pm10变量。PM2.5排放取决于你,但请记住这两个的区别是什么(请参阅这和这)。