Add getting and converting ERA5 ml analyses to MPAS ICs#418
Merged
Add getting and converting ERA5 ml analyses to MPAS ICs#418
Conversation
jim-p-w
reviewed
Feb 20, 2026
jim-p-w
reviewed
Feb 25, 2026
jim-p-w
approved these changes
Mar 5, 2026
Collaborator
Author
|
Thanks for your reviews and tests, @jim-p-w and @junmeiban! |
jim-p-w
reviewed
Mar 5, 2026
| base = 'UngribExternalAnalysis' | ||
| queue = 'UngribExternalAnalyses' | ||
| if base in self['PrepareExternalAnalysisOuter']: | ||
| if base in self['PrepareExternalAnalysisOuter'] and self['externalanalyses__UngribPrefixOuter'] != 'ERA5': |
Collaborator
There was a problem hiding this comment.
When I try to run test/testinput/3dvar_OIE120km_ColdStart.yaml this line is throwing an exception. It looks like the key externalanalyses__UngribPrefixOuter doesn't exist for that scenario.
Collaborator
Author
There was a problem hiding this comment.
Thanks for pointing this out @jim-p-w! I should have checked this more carefully. I found a solution by doing if self.__ungribtask is not None: instead. I will add this fix into my open PR.
11 tasks
ibanos90
added a commit
that referenced
this pull request
Mar 11, 2026
### Description The graphics package can already handle multiple reference datasets for model-space verification (see [JCSDA-internal/mpas-jedi#1096]). However, this capability has not been integrated into the MPAS-Workflow. This PR introduces support for specifying the reference dataset (-R {GFS, ERA5, EC}) at the statistics generation stage, i.e., in the verifymodel task. This change is needed because the statistics files themselves contain the reference dataset definition. For example, `"log_mogfsan", "mmgfsan"`, indicates GFS as the reference dataset used in the statistics. Here, it was added a new configurable variable `anaRef` in the VerifyModel component that will be passed to the verifymodel.csh script as an argument (-R $anaRef). The option was added to `ForecastFromGFSAnalyses` and it works as expected. Additionally, the `script directory` variable is returned to be user-configurable in both VerifyModel and VerifyObs. The path was updated to use the tip of the develop branch similar to the build. ### Issue closed Closes #419 Fixes the issue mentioned by @jim-p-w in #418 (comment) ### Tests completed #### Tier 1: - [x] 3dvar_OIE120km_WarmStart - [x] 3denvar_OIE120km_IAU_WarmStart - [x] 3dvar_OIE120km_ColdStart - [x] 3dvar_O30kmIE60km_ColdStart - [x] 3denvar_O30kmIE60km_WarmStart - [x] eda_OIE120km_WarmStart - [x] getkf_OIE120km_WarmStart - [x] ForecastFromGFSAnalysesMPT - [x] 4denvar_OIE120km_WarmStart - [x] 4dhybrid_OIE120km_WarmStart #### Tier 2 (: - [x] ForecastFromGFSAnalyses Under `diagnostic_stats/model`, `myCommand` successfully shows: `python DiagnoseModelStatistics.py 2018041500 -n 128 -r /glade/derecho/scratch/ivette/pandac/ivette_eda_OIE120km_WarmStart_TEST/ExternalAnalyses/120km/2018041500/x1.40962.init -R GFS -rd /glade/derecho/scratch/ivette/pandac/ivette_eda_OIE120km_WarmStart_TEST/ExternalAnalyses/120km/2018041500/diag -m 5`
5 tasks
This file contains hidden or 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
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
3 participants
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.
Description
This PR adds the capability to get and convert ERA5 model level analyses to MPAS format. It relies on era5_to_int code to get the ERA5 ml data on GDEX and the NetCDF files to WPS intermediate. Then,
mpas_initis run to create the MPAS init files. A new scenario is added to test this functionality.Issue closed
None
Tests completed