Closed 1kyne-256 closed 2 years ago
Hi @1kyne-256 the model works in molecule cm-3
(concentrations) and seconds
(time). The units of the reaction rates are calculated accordingly.
productionRates
and lossRates
are simply the reaction rates of each reaction that forms and consumes a given molecule. The values should be the same as those in the reactionRates/
directory, only in a more user-friendly format for selected species of interest. Hope this helps.
Thank you very much .
There's another problem.
If I wanted to calculate the net photochemical production rate of ozone, what data would I use to calculate it(PO3)?
The calculation formula is as follows:
Photochemical reaction is the dominant source of near-surface ozone, and its process can be represented by (R1)-(R3) :
During the whole process, peroxic radicals (HO2 and RO2) oxidize NO to produce NO2(φ describes the yield of NO oxidized by RO2 to produce NO2, and its value is less than 1), which then photolyzes to produce O3.
The total photochemical generation rate of ozone can be expressed as:
At the same time, ozone can also be removed by some photochemical reactions, such as (R4)-(R8):
The photochemical removal rate of ozone can be expressed as:
The net photochemical production rate of ozone (P(O3)) is the difference between G(O3) and D(O3).
@1kyne-256 you can set up the model to output the production and loss rates of OH, HO2, O1D, O3, NO and the reaction rates will be the terms in your G(O3) and D(O3). It will take quite a bit of processing of the output files I am afraid. Also be careful that the model must NOT be constrained to O3 and/or NO otherwise the results won't make sense.
@rs028 Thank you very much . So that means that the output of O3 production rates is not G (O3), and the output of O3 loss rates is not D (O3). I have to do a lot of processing to get to my final goal, right? At one time, I thought, all I had to do was set the productivity and loss rate for O3, and then I could calculate P (O3).Now it seems that I misunderstood.
That's correct: productionRates
and lossRates
will give you the reaction rates of each individual reaction that forms and destroys the selected species, not their sum.
What are the units of the rate variable in the productionRates.output and lossRates.output files? How is productionRates and lossRates calculated? What's the formula for the rate of formation? For example, what about the productionRate of O3 and the lossRates of O3 consumption?![图片1](https://user-images.githubusercontent.com/64003532/141946589-aeaa85ec-cecd-459e-84d7-a756bae38cd1.png)