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

Retire ExtState%WLI (land-water-ice index) #170

Merged
merged 3 commits into from
Jan 5, 2023

Conversation

lizziel
Copy link
Contributor

@lizziel lizziel commented Oct 11, 2022

This PR is 1 of 2 for super-projects. See instructions here:
geoschem/GCHP#264
geoschem/GCClassic#26

This pull request retires use of external meteorology field land-water-ice (LWI) index and the ExtState array it is stored in (ExtState%WLI). We are removing it as a dependency because it will not be in the GEOS-IT dataset and determining whether a grid cell is over ocean, ice, or land is possible through other input meteorology fields, specifically GMAO meteorology constants FRLAND, FRLANDIC, FRLAKE, and FROCEAN, and time-varying fields FRSEAICE and FRSNO.

The ExtState%WLI field is currently used in HEMCO to retrieve land type via subroutine HCO_LANDTYPE in hco_geotools_mod.F90. That subroutine is called in the following extension modules, which will be impacted by this update:

  1. hcox_dustdead_mod
  2. hcox_gc_POPs_mod
  3. hcox_iodine_mod
  4. hcox_lightnox_mod
  5. hcox_seaflux_mod
  6. hcox_soilnox_mod
  7. hcox_tomas_dustdead_mod
  8. hcox_tomas_jeagle_mod

Note that the LWI input field used for WLI in HEMCO is currently computed during GEOS-Chem input pre-processing from the other GMAO surface type met-fields. Changing to using these fields at run-time will result in differences in the model. In addition, the new way for determining if mostly land, ice, or ocean is to evaluate which surface type fraction is the greatest. This is a slightly different algorithm than what was used to compute LWI and this will also cause small differences.

There is an accompanying pull request in GEOS-Chem that should be merged at the same time:
geoschem/geos-chem#1435.

LWI is the land-water-ice index included in input meteorology. It is
derived from fraction land types which are also available as input.
It is therefore retired in favor of using the fractions which are
already read into in HEMCO.

Signed-off-by: Lizzie Lundgren <[email protected]>
@lizziel lizziel added the topic: Input Data Related to input/emissions data, or disk read/write operations label Oct 11, 2022
@lizziel lizziel added this to the 3.5.1 milestone Oct 11, 2022
@lizziel lizziel self-assigned this Oct 11, 2022
@lizziel lizziel changed the base branch from main to dev October 11, 2022 20:37
@msulprizio msulprizio changed the base branch from dev to main October 18, 2022 17:29
@msulprizio msulprizio modified the milestones: 3.5.1, 3.6.0 Oct 18, 2022
@msulprizio msulprizio changed the base branch from main to dev/3.6.0 October 18, 2022 17:32
@lizziel
Copy link
Contributor Author

lizziel commented Nov 4, 2022

@christophkeller, this PR changes the subroutine HCO_LANDTYPE in hco_geotools_mod.F90. Since you are the original author of that subroutine could you review the changes I made to that file? We are swapping out use of input LWI with using the other GMAO surface type fields.

@lizziel lizziel marked this pull request as ready for review November 4, 2022 18:36
@lizziel lizziel requested a review from yantosca November 4, 2022 18:36
@lizziel
Copy link
Contributor Author

lizziel commented Nov 4, 2022

Results of a 1-month GEOS-Chem benchmark comparison with and without this update show differences in the land/ocean/ice boundaries, as shown here:
Screen Shot 2022-10-17 at 3 57 03 PM
Screen Shot 2022-10-17 at 3 56 51 PM

@stale
Copy link

stale bot commented Dec 16, 2022

This issue has been automatically marked as stale because it has not had recent activity. If there are no updates within 7 days it will be closed. You can add the "never stale" tag to prevent the Stale bot from closing this issue.

@stale stale bot added the stale No recent activity on this issue label Dec 16, 2022
@msulprizio msulprizio removed the stale No recent activity on this issue label Dec 16, 2022
@yantosca yantosca added never stale Never label this issue as stale topic: Structural Modifications Related to HEMCO structural modifications (as opposed to scientific updates) labels Dec 16, 2022
Copy link
Contributor

@yantosca yantosca left a comment

Choose a reason for hiding this comment

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

This PR now sends uses the GMAO land fraction fields as arguments to the HCO_Landtype_Sp and HCO_LandType_Dp routines. Good to merge!

@msulprizio msulprizio merged commit 7de865f into dev/3.6.0 Jan 5, 2023
@msulprizio msulprizio deleted the feature/retire_wli_metfield branch January 5, 2023 23:41
zmoon added a commit to noaa-oar-arl/NEXUS that referenced this pull request Apr 26, 2023
This was removed for HEMCO v3.6

In geoschem/HEMCO#170
they just removed the relevant block from the standalone mod,
as we do here
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
never stale Never label this issue as stale topic: Input Data Related to input/emissions data, or disk read/write operations topic: Structural Modifications Related to HEMCO structural modifications (as opposed to scientific updates)
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

3 participants