日志 - 如何增加WRFout文件中的化学诊断变量

Legendary

在模拟大气臭氧污染时，一些诊断量比如臭氧的净生成率对污染过程的研究具有十分重要的意义。然而，wrf-chem并没有输出这些变量到wrfout文件中。实际上只要简单一步操作，就能轻松实现将化学反应诊断变量输出到文件。

1.打开namelist.input
2.找到文件中化学过程参数的设置部分 &chem
3.寻找chemdiag这一参数，如果没有，则添加在&chem 这一部分参数的最后
4.chemdiag=0 turns off chemical tendency diagnostics
  chemdiag= 1 turns on chemical tendency diagnostics for equation budget analysis
  所以我们此时应该令chemdiag=1

注意：
1. 打开chemdiag=1后，计算过程会变得异常漫长，文件体积变大。
2. 打开后运行过程可能会莫名出错。如果你使用的mpirun wrf.exe np=X，那么尝试减小X的数值，会增加运行成功的概率。
3. 更改namelist.input 中debug_level的设置，有可能反馈更多运行错误的信息。
常见debug_level = 0,50,100,200,300,数值越大越详细。

模式报错：
WARNING: When using an adaptive time-step the boundary layer time-step should be 0 (i.e., equivalent to model time-step). In order to proceed, for surface calculations, the boundary layer time-step will be rounded to the nearest minute, possibly resulting in inaccurate results.
运行结果：
关闭chemdiag
（待补充）
打开chemdiag

相关链接：
该链接描述了新wrfout里变量的物理意义
https://svn-dares-dart.cgd.ucar. ... FIXES/registry.chem

同时提到了deposition的相关变量（化学物质的沉降和去除还不知道怎么开启模拟，求回帖）

相关应用文献：
Elucidating the ozone pollution in Yangtze River Delta region during the 2016 G20 summit for MICS-Asia III

IN`VIZEL

您好，
想请问一下，我模拟时打开了chemdiag=1，但是生成的wrfout文件中没有找到相应的变量，如conv_o3, chem_o3...... 我是用ARWpost来处理wrfout文件的，请问是不是在namelist.ARWpost中除了把对应的fields写上还需要其他处理？或者在编译WRF-Chem时需要额外设定来打开chemdiag？感谢

诶诶诶诶诶诶呀

IN`VIZEL 发表于 2021-3-22 12:48
您好，
想请问一下，我模拟时打开了chemdiag=1，但是生成的wrfout文件中没有找到相应的变量，如conv_o3, c ...

还需要在WRF/Registry文件夹下面的registry.chem里进行修改，找到chemdiag的有关选项，如果要让conc_o3输出出来，就把这个量从r改成rh，然后重新编译WRF和WPS，重跑一边，就可以看到了

棉花糖grace

诶诶诶诶诶诶呀 发表于 2022-9-5 14:22
还需要在WRF/Registry文件夹下面的registry.chem里进行修改，找到chemdiag的有关选项，如果要让conc_o3输 ...

你好，请问registry里chemdiag的变量是在哪里啊，我用的3.9版本，没有找到呢

诶诶诶诶诶诶呀

棉花糖grace 发表于 2023-2-6 16:42
你好，请问registry里chemdiag的变量是在哪里啊，我用的3.9版本，没有找到呢

我用的4.2版本的，打开这个文件搜conc_o3就能搜到了，可能是版本不同的原因

棉花糖grace

诶诶诶诶诶诶呀 发表于 2023-2-7 15:06
我用的4.2版本的，打开这个文件搜conc_o3就能搜到了，可能是版本不同的原因

你好，请问当前时刻化学过程对O3生成的贡献=当前时刻chem_o3 减去 上一时刻chem_o3，是这样计算吗

tyf

棉花糖grace 发表于 2023-2-16 21:22
你好，请问当前时刻化学过程对O3生成的贡献=当前时刻chem_o3 减去 上一时刻chem_o3，是这样计算吗

也想知道这个！

诶诶诶诶诶诶呀

棉花糖grace 发表于 2023-2-16 21:22
你好，请问当前时刻化学过程对O3生成的贡献=当前时刻chem_o3 减去 上一时刻chem_o3，是这样计算吗

对的，前后两个时刻做减法

诶诶诶诶诶诶呀

tyf 发表于 2023-2-21 17:20
也想知道这个！

对的，是当前时刻的chem_o3减去上一时刻的chem_o3

tyf

诶诶诶诶诶诶呀 发表于 2023-3-21 13:00
对的，是当前时刻的chem_o3减去上一时刻的chem_o3

谢谢！！！