geoschem / geos-chem

GEOS-Chem "Science Codebase" repository. Contains GEOS-Chem science routines, run directory generation scripts, and interface code. This repository is used as a submodule within the GCClassic and GCHP wrappers, as well as in other modeling contexts (external ESMs).
http://geos-chem.org
Other
170 stars 166 forks source link

Tracking OH changes in simulations #2543

Open fang159753 opened 3 weeks ago

fang159753 commented 3 weeks ago

Your name

fangjiong

Your affiliation

shanghai university

Please provide a clear and concise description of your question or discussion topic.

I'm new to using GEOS-Chem and have encountered some issues during actual operation. I ran simulations by combining fullchemistry with CH4 chemistry. First, I obtained SpeciesConcVV_OH and OHconcAfterChem from the full-chemistry run. In the metrics.py script, I used the LossOHbyCH4columnTrop variable to calculate CH4 lifetime. I'm curious to know if the OH referenced in LossOHbyCH4columnTrop corresponds to SpeciesConcVV_OH, OHconcAfterChem, or another variable?

In fullchem_mod.F90, I found that LossOHbyCH4 = LossOHbyCH4 + ( Ktrop OHconc_mcm3 airMass_m ), where OHconc_mcm3 = State_Chm%Species(id_OH)%Conc(I,J,L). However, I'm not sure which OH variable in which file State_Chm%Species(id_OH)%Conc(I,J,L) refers to.

Additionally, when I use SpeciesConcVV_OH from the fullchemistry simulation as input for CH4 chemistry, should I look at OHconcAfterChem for the OH output? Interestingly, I’ve noticed that OHconcAfterChem values are higher than either SpeciesConcVV_OH or OHconcAfterChem obtained in full-chemistry. Could you help me understand where I might have gone wrong?

Your help would be greatly appreciated! Thank you so much! Snipaste_2024-10-31_15-48-50