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

Updates from GMAO #125

Merged
merged 24 commits into from
Jan 21, 2022
Merged

Conversation

lizziel
Copy link
Contributor

@lizziel lizziel commented Jan 11, 2022

This PR brings in changes made to HEMCO for use in GEOSgcm with GEOS-Chem 13.3.4.

Changes include:

  1. Skip check for PBLHEIGHT field if not running GEOS. There are HEMCO applications in GEOS that don't require PBLHEIGHT to be defined.
  2. Change criteria for assuming restart field is all fill values. The new criteria is max value is negative. Previously it was global sum equal to zero.
  3. Force the soil NOx code to grab the restart fields from the internal state on each call if using GEOS, rather than rely on whatever is in local memory.
  4. Rework how to determine when to read volcano table. The new method is to keep track of the emissions data day that is in memory (new variable YmdOnFile) and then compare the current simulation day against that.
  5. Edits to prevent undefined variable when calculating vertical scale factor. This updates adds handling for when HcoState%Grid%BXHEIGHT_M%Val is associated and HcoState%Options%VertWeight is false.

Christoph Keller and others added 19 commits April 19, 2021 10:05
We now have two copies of files in `src/Shared`. This PR removes the
older versions (without `hco_` prefix)
…files

Remove redundant files in src/Shared
# Conflicts:
#	CMakeLists.txt
#	run/createRunDir.sh
#	src/Core/CMakeLists.txt
#	src/Core/hco_config_mod.F90
#	src/Core/hcoio_read_std_mod.F90
#	src/Core/hcoio_util_mod.F90
#	src/Extensions/hcox_gc_POPs_mod.F90
#	src/Extensions/hcox_seasalt_mod.F90
#	src/Interfaces/GCHPctm/HEMCO_GridCompMod.F90
#	src/Shared/GeosUtil/hco_julday_mod.F90
#	src/Shared/NcdfUtil/hco_ncdf_mod.F90

Signed-off-by: Lizzie Lundgren <[email protected]>
…SMF. For many HEMCO applications within GEOS, the PBL height is not needed and the code works without that information. There is an explicit check for PBL height fields if it is indeed needed.
…mco_3.2.2

# Conflicts:
#	src/Core/hco_calc_mod.F90
#	src/Extensions/hcox_megan_mod.F90

Signed-off-by: Lizzie Lundgren <[email protected]>
@lizziel lizziel added the category: Feature Request New feature or request label Jan 11, 2022
@lizziel lizziel self-assigned this Jan 11, 2022
@lizziel
Copy link
Contributor Author

lizziel commented Jan 12, 2022

This new line in CMakeLists.txt, added for 3km GEOS-Chem runs, is problematic:

# Enable OpenMP
target_link_libraries(${this} PRIVATE OpenMP::OpenMP_Fortran)

Enabling OpenMP in HEMCO by default was added for using GEOS with
3 km runs. Rather than explicitly set OpenMP in HEMCO CMakeLists.txt, the
OMP setting should be used higher up in the model. If set to true then
HEMCO will build with OpenMP enabled.

Signed-off-by: Lizzie Lundgren <[email protected]>
@christophkeller
Copy link
Contributor

No need to include the new line in CMakeLists.txt that was added for the 3km test runs.

Copy link
Contributor

@christophkeller christophkeller left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

There was a merge error in hco_calc_mod.F90, see my comment there. Everything looks good otherwise.

src/Core/hco_calc_mod.F90 Outdated Show resolved Hide resolved
@lizziel lizziel changed the title [WIP] Updates for HEMCO in GEOSgcm Updates for HEMCO in GEOSgcm Jan 19, 2022
@lizziel lizziel changed the title Updates for HEMCO in GEOSgcm Updates from GMAO Jan 19, 2022
@lizziel lizziel added this to the 3.4.0 milestone Jan 19, 2022
@lizziel lizziel marked this pull request as ready for review January 19, 2022 17:44
src/Core/hco_calc_mod.F90 Outdated Show resolved Hide resolved
src/Extensions/hcox_soilnox_mod.F90 Show resolved Hide resolved
src/Extensions/hcox_volcano_mod.F90 Show resolved Hide resolved
Copy link
Contributor

@christophkeller christophkeller left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks good now.

This update was originally added by C. Keller for GEOS only. With this
updatet is now standard handling when calculating dilution factor. It is
only executed when HcoState%Grid%BXHEIGHT_M%Val is associated and
HcoState%Options%VertWeight is false.

Signed-off-by: Lizzie Lundgren <[email protected]>
@lizziel lizziel force-pushed the feature/ewlundgr/hemco_3.2.2_in_geos branch from 89f4f9e to 570eb98 Compare January 19, 2022 20:00
Signed-off-by: Lizzie Lundgren <[email protected]>
@yantosca yantosca merged commit f31ffd6 into geoschem:dev Jan 21, 2022
@yantosca
Copy link
Contributor

FYI, this fix has been tagged in GC-Classic and GCHP as 13.4.0-rc.2. I have a GCHP benchmark pending on Cannon and an GC-Classic benchmark running on the AWS cloud.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
category: Feature Request New feature or request
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

6 participants