[WIP] Add capability to run GEOS-Chem at 0.125x0.15625 resolution #1980
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
…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 <[email protected]>
msulprizio
added
category: Feature Request
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 <[email protected]>
…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 <[email protected]>
…tion Signed-off-by: Melissa Sulprizio <[email protected]>
|
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 <[email protected]>
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 <[email protected]>
…iaolin Wang The PS and SPHU fields at 0.125x0.15625 resolution are in the I1dyn collection, while the U and V fields are read from the A1dyn collection. All of these fields should be read hourly, not 3-hourly as is the case for 0.25x0.3125 resolution and coarser. New subroutines FlexGrid_Read_A1dyn, FlexGrid_Read_I1dyn_1, FlexGrid_Read_I1dyn_2, and Copy_I1dyn_Fields have been added to flexgrid_mod.F90 following work by Xiaolin Wang. The code in hco_interface_gc_mod.F90 and main.F90 that calls these routines has also been updated. New I1dyn routines have also been added to time_mod.F90. Signed-off-by: Melissa Sulprizio <[email protected]>
…tion Signed-off-by: Melissa Sulprizio <[email protected]>
yantosca
added
topic: Structural Modifications
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
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, SPHUGEOS.fp.asm.A1dyn: U, VSo far only run directory files have been updated in this PR. The script
createRunDir.shhas 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)