Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

Chemistry tendency isn't equal to net IRR

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

zxdawn

Member
Dear All,

I'm trying to get the net O3 chemistry production and the contribution of each reaction.

The net O3 chemistry production is just the chemistry tendency variable named chem_o3 and the contributions can be calculated in the equation:
ozone_reaction.jpg
(Equation 1 in https://doi.org/10.1029/JD095iD09p14049)

In my opinion, dO3 should be equal to chem_o3.
So, I output the IRR and chemistry tendency with the same time step.
Here's the summary of used variables in the IRR output file:

Code:
O3 Production:
['NO_HO2_IRR',
'ALKO2_NO_IRR', 'C2H5O2_NO_IRR', 'C3H7O2_NO_IRR', 'CH3CO3_NO_IRR', 'CH3O2_NO_IRR',
'ENEO2_NO_IRR', 'EO2_NO_IRR', 'ISOPO2_NO_IRR', 'MACRO2_NO_IRR', 'MEKO2_NO_IRR',
'PO2_NO_IRR', 'RO2_NO_IRR', 'TERPO2_NO_IRR', 'TOLO2_NO_IRR', 'XO2_NO_IRR']

O3 Loss:
['O1D_CB4_H2O_IRR', 'OH_O3_IRR', 'HO2_O3_IRR']

Substracting Loss from Production yields the Net IRR (named Accumulated IRR in the figure).

However, if I plot the IRR and Chemistry Tendency at the same grid and level together, the values are not equal:
comp_irr_chemdt.jpg

Note that the IRR Change and Chemistry Tendency Change are the differences by the Time dimension.
Although the trend of the accumulated values seems the same, the changes are opposite at some times.

If you have any ideas, please don't hesitate ;)

Regards,
Xin
 
Top