尝试在CMAQ中添加关于干沉降物的新反应

Fortuna

我用的是CMAQv5.3.2，沉降方面是m3dry模块，现在想要添加dNO3 = HONO + NO2这样一个新反应，dNO3代表干沉降NO3，反应速率是调用现有的HNO3光解速率，所以这本身也是一个非均相光解反应。
查了一些资料，试图使用CHEMMECH工具把这个反应加入mech.def、RXNS_DATA和RXNS_FUNC文件，在这之前需要把dNO3这个新物种写入AE.nml或者GC.nml文件，问题在于模型计算时要怎么识别到dNO3代表的是干沉降NO3，计算时用到的是dNO3的生成速率还是浓度？
如果有其他能实现这个反应的修改思路也请告诉我，我想多做一些尝试。

zsx1997

您好，我也在优化反应机理，是否可以反应速率乘以倍数来表明部分HNO3沉降发生了反应。
我现在在添加HNO3置换反应，但也遇到了一些问题，能交流交流吗？

Fortuna

您好，您的意思是直接使用HNO3光解速率乘上HNO3光解反应与dNO3反应之间的倍数，从而直接得到dNO3反应生成的额外HONO和NO2吗？考虑到我并不关注NO3最终的干沉降量，这个思路似乎是可行的。我的想法是利用CHEMMECH工具，直接在mech.def中添加反应方程式的右侧部分，RXNS_FUNC文件中计算的photolysis rates变量“RKI”的单位是ppm/min，最终效果可能类似于凭空生成HONO和NO2，会对沉降N一定的高估。

另外，您有使用WRF-CMAQ耦合模型的经验吗？我最近在耦合模型编译方面遇到问题。