Skip to content
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

[WIP] Change units of Loss_Ox diagnostic from molec/cm3/s to unitless (to fix an issue in tagO3 simulations) #2096

Merged
merged 20 commits into from
Feb 9, 2024

Conversation

yantosca
Copy link
Contributor

@yantosca yantosca commented Jan 3, 2024

Name and Institution (Required)

Name: Bob Yantosca
Institution: Harvard + GCST

Confirm you have reviewed the following documentation

Describe the update

This is a a placeholder PR and should not be merged yet.

As described in #2094 (by @ncolombi and @xpye98), the units of the Loss_Ox diagnostic need to be changed from molec/cm3/s to unitless in order to prevent a problem in the tagO3 simulation. Please see the discussion in #2094 for the time being.

What needs to be changed:

  • Units of Loss_Ox as archived from GeosCore/fullchem_mod.F90. Specifically, we need to divide the C(ind_LO3) by the mass of all species in the LOx family.
  • Routine Chem_Tagged_O3 in GeosCore/tagged_o3_mod.F90, to ingest Loss_Ox in the new units and apply the loss rates.

Caveat: We will also need to make sure that this does not affect other specialty simulations that may use archived Loss_Ox data.

Tagging @msulprizio @lizziel @ncolombi @xpe98 for your reference.

Expected changes

This should be a zero-diff update, as we are only changing the diagnostic output. However, because we are touching the full chemistry module as well, we should benchmark this just to be sure. I recommend that this update be added to the next Y version.

Reference(s)

N/A

Related Github Issue(s)

lizziel and others added 20 commits July 21, 2023 13:31
With this commit the GEOS model builds and runs using GEOS-Chem transport
tracer simulation with new config files stored in
geos-chem/run/GEOS/TransportTracers. The following is currently required
to run successfully:
1. CEDSv2 must be turned off in HEMCO_Config.rc
2. ExtData2G must be turned off in CAP.rc

Signed-off-by: Lizzie Lundgren <[email protected]>
This commit addresses issue #2055, by using a newer Hg restart file

run/shared/download_data.yml
- Change the remote restart directory from v2012-12 to v2023-12

CHANGELOG.md
- Updated accordingly

Signed-off-by: Bob Yantosca <[email protected]>
This commit addresses issue #2029 by @jessicahaskins:

KPP/fullchem/commonIncludeVars.H
- Update units (in comment) for SR_MW in comment line
- Add units (in comment) for PRESS, RELHUM, SUNCOS, PHOTOL, H2O,
  K_STRAT, K_TROP, TROP

Signed-off-by: Bob Yantosca <[email protected]>
This merge brings PR #2061 (Add missing units in comments of
KPP/fullchem/commonIncludeVars.sh, by @yantosca) into the
"no-diff-to-benchmark" development stream.

This is a zero-diff update which updates comments in the fullchem
mechanism include file commonIncludeVars.sh.  The KPP mechanisms were
not rebuilt, as this file is #include-d into the mechanism.

Signed-off-by: Bob Yantosca <[email protected]>
The timers used in benchmarking to track model performance and generate this
chart https://wiki.seas.harvard.edu/geos-chem/index.php/GEOS-Chem_1-month_benchmark_timing_results
include:

  - GEOS-Chem,
  - HEMCO
  - All chemistry
  - Gas-phase chem
  - Photolysis
  - Aerosol chem
  - Linearized chem
  - Transport
  - Convection
  - Boundary layer mixing
  - Dry deposition
  - Wet deposition
  - Diagnostics
  - Unit conversion

The remaining timers are superfluous. They are removed here to clean up
output and reduce unnecessary time spent computing those timers.

Signed-off-by: Melissa Sulprizio <[email protected]>
This merge brings PR #2059 (Now read the Hg restart file from
GEOSCHEM_RESTARTS/v2023-12, by @yantosca) into the GEOS-Chem
"no-diff-to-benchmark" development stream.

This PR now copies updated Hg simulation restart files (that
with well-spun-up initial conditions) to GCClassic and GCHP run
directories, as was requested by @angothelene in #2055.

Signed-off-by: Bob Yantosca <[email protected]>
Specific humidity will be available in transport tracer sim checkpoints
for purposes of post-processing unit conversions.

Signed-off-by: Lizzie Lundgren <[email protected]>
This avoids an out-of-bound issue when vertically flipping
diagnostics with 73 levels in GEOS.

Signed-off-by: Lizzie Lundgren <[email protected]>
This merge brings PR #2069 (Reduce timers to essential list only, by
@msulprizio) into the "no-diff-to-benchmark" development stream.

This removes GEOS-Chem timers that are no longer needed.

Signed-off-by: Bob Yantosca <[email protected]>
This merge brings PR #2058 (GEOS-Chem Transport Tracers simulation GEOS,
by @lizziel) into the GEOS-Chem "no-diff-to-benchmark" development stream.

This adds several updates to get the GEOS-Chem Transport Tracers simulation
working in the NASA GEOS ESM.  It will not affect the GEOS-Chem Classic or
GCHP implementations of GEOS-Chem.

Signed-off-by: Bob Yantosca <[email protected]>
This commit fixes the issue raised in PR #2083 by @lizziel.  Integration
and parallel test execute scripts now use the run directory name and not
the absolute path to determine the executable file name.

test/integration/GCClassic/integrationTestExecute.sh
test/integration/GCHP/integrationTestExecute.sh
test/parallel/GCClassic/parallelTestExecute.sh
- Pass ${runDir} instead of ${runAbsPath} to the exe_name function.
  This will prevent integration test root paths with e.g. "rrtmg"
  from returning an incorrect executable name.

Signed-off-by: Bob Yantosca <[email protected]>
CHANGELOG.md
- Added note that we now use the run directory and not the absolute
  path to determine the executable file to use in integration and
  parallel tests.

Signed-off-by: Bob Yantosca <[email protected]>
This merge brings PR #2084 (Use rundir (not absolute path) to select
the executable file for integration and parallel tests, by @yantosca)
into the GEOS-Chem "no-diff-to-benchmark" development stream.

This PR fixes the issue reported in #2083 by @lizziel.

Signed-off-by: Bob Yantosca <[email protected]>
@yantosca yantosca added topic: Diagnostics Related to output diagnostic data category: Bug Fix Fixes a previously-reported bug topic: Math and Units Issues pertaining to math expressions and/or unit conversions labels Jan 3, 2024
@yantosca yantosca self-assigned this Jan 3, 2024
@msulprizio msulprizio merged commit def47a0 into main Feb 9, 2024
@msulprizio
Copy link
Contributor

@yantosca Were there any commits added to this branch (bugfix/tago3-pl-units)? It looks like it was just a branch made off of dev/no-diff-to-benchmark but nothing was actually added to it. If so, can we now delete this branch? If the issue is not resolved, I recommend creating a new branch/pull request (and reopening Nadia's original issue).

Copy link
Contributor Author

yantosca commented Feb 9, 2024

@msulprizio, yes... this was a branch I used to open a draft PR. I can close this out.

@msulprizio msulprizio deleted the bugfix/tago3-pl-units branch February 9, 2024 21:34
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
category: Bug Fix Fixes a previously-reported bug topic: Diagnostics Related to output diagnostic data topic: Math and Units Issues pertaining to math expressions and/or unit conversions
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Issue with units of O3 loss in tagged simulation
3 participants