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] Add capability to run GEOS-Chem at 0.125x0.15625 resolution #1980

Draft
wants to merge 5 commits into
base: main
Choose a base branch
from

Conversation

msulprizio
Copy link
Contributor

Name and Institution (Required)

Name: Melissa Sulprizio
Institution: Harvard / GCST

Describe the update

NASA GMAO has provided 1-hourly mass flux fields for GEOS-FP meteorology at c720 resolution to drive GCHP. The mass flux fields have been converted to derived winds fields and further regridded from c720 to 0.125x0.15625 lat-lon resolution. The new files include the following collections and met fields:

  • GEOS.fp.asm.I1dyn : PS, SPHU
  • GEOS.fp.asm.A1dyn : U, V

So far only run directory files have been updated in this PR. The script createRunDir.sh has been updated to include this resolution as an option when running with GEOS-FP meteorology. The use of this resolution and these fields is still being evaluated for it's impact on model output and performance.

Expected changes

This update will not impact the full-chemistry benchmark simulations which are run at 4x5 resolution.

Related Github Issue(s)

…lution

NASA GMAO has provided 1-hourly mass flux fields for GEOS-FP meteorology
at c720 resolution to drive GCHP. The mass flux fields have been converted
to derived winds fields and further regridded from c720 to 0.125x0.15625
lat-lon resolution. The new files include the following collections and met
fields:

- GEOS.fp.asm.I1dyn : PS, SPHU
- GEOS.fp.asm.A1dyn : U, V

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
@msulprizio msulprizio added category: Feature Request New feature or request topic: Input Data Related to input data topic: Nested Grid Simulation Related to GEOS-Chem nested model simulations labels Oct 6, 2023
@msulprizio msulprizio self-assigned this Oct 6, 2023
Several IF statements have been expanded to include State_Grid%GridRes =
0.125x0.15625 resolution so that grid parameters and scale factors are
properly determined.

Also updated createRunDir.sh for GCClassic to correctly set transport and
chemistry timesteps for 0.125x0.15625 resolution. For now the same default
timesteps as 0.25x0.3125 and 0.5x0.625 will be used (transport=5 min,
chemistry=10 min). This recommendation may later changed based on proper
evaluation of the 0.125x0.15625 simulations.

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
…s for non-dynamic fields

The 0.125x0.15625 files are only available for dynamic collections A1dyn
and I1dyn. The remaining met fields when running at this resolution should
be read at 0.25x0.3125. This file is currently hardcoded to use the North
America (NA) files at 0.25x0.3125 resolution for testing, but this should
be made more flexible in the future.

The 0.125x0.15625 files have also been updated here so that daily files are
used instead of hourly files. The use of daily 24-hour met field files has
already been implemented in GEOS-Chem, but it does not currently know how
to handle hourly files.

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
…tion

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
@msulprizio
Copy link
Contributor Author

There seems to be an issue in input_mod.F90 cutting off the longitude resolution string to 4 decimal places and changing 0.15625 to 0.1563. This causes differences in the longitude edge and center values saved to netCDF output. I will push a fix shortly and post difference plots to show the impact on model output.

Tagging @xlwangg

Xiaolin Wang identified the following bugs in the 12km simulations:

1. In GeosCore/pjc_pfix_window_mod.F90 for 12km code, we need to change the
   parameter dp in correspond with 12km resolution.
2. In the input_mod.F90, we need to improve floating point precision,
   otherwise the longitude resolution is 0.1563 rather than 0.15625.

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
@msulprizio
Copy link
Contributor Author

The above issue in the longitude string length was already fixed by @xlwangg and is now also fixed in commit cb71dfa

msulprizio added a commit to geoschem/integrated_methane_inversion that referenced this pull request Apr 24, 2024
…n Wang

The grid resolution option is expanded to include 0.125x0.15625. These
modifications were obtained from Xiaolin Wang. These simulations require
using the feature/0.125x0.15625_resolution in the GCClassic/src/GEOS-Chem
repository (see also geoschem/geos-chem#1980).

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
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 topic: Input Data Related to input data topic: Nested Grid Simulation Related to GEOS-Chem nested model simulations
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

1 participant