-
Notifications
You must be signed in to change notification settings - Fork 24
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
[FEATURE REQUEST] Add OH metrics output for benchmarks #34
Comments
The calculations for MCF and CH4 lifetimes require the time-in-troposphere diagnostic output from GEOS-Chem. This is currently only saved out to bpch diagnostics and needs to be implemented in History. @yantosca is working on that now. |
The time in troposphere netCDF diagnostic is now implemented in commit geoschem/geos-chem@c01871c. A 1-day zonal mean of ND54 (bpch) vs. FracOfTimeInTrop (netCDF) looks like this: |
This is now pushed in 509cb4f. We still need to apply fraction of time in troposphere to calculation to include only purely tropospheric boxes. This will be added once the 12.6.0 1-month benchmark with the new FracOfTimeInTrop diagnostic is available. |
This calculation now also works for GCHP but the values are off. For example,
Further investigation is required. |
We should probably implement this as a HISTORY diagnostic in GEOS-Chem, so that we can obtain the output for GCHP as well. This is documented in GEOS-Chem pull request https://github.com/geoschem/geos-chem/242. As a stopgap measure, I also have implemented a GCPy function (mean_oh_from_logs.py) to strip the 1-year mean OH from the log files and to print them. But this will not work for GCHP. I will close out this issue, as it seems that further work will need to be done as a GEOS-Chem diagnostic as described in the issue above. |
Reproduce the GAMAP metrics output in GCPy so that we do not need to generate with IDL. Ideally, we would compute and compare global mean OH, MCF lifetime, and CH4 lifetime.
The text was updated successfully, but these errors were encountered: