dtcenter / MET

Model Evaluation Tools
https://dtcenter.org/community-code/model-evaluation-tools-met
Apache License 2.0
79 stars 24 forks source link

Enhance Ensemble-Stat to support the use of Gen-Ens-Prod output in the -ens_mean command line option. #2139

Open JohnHalleyGotway opened 2 years ago

JohnHalleyGotway commented 2 years ago

Describe the Enhancement

After the introduction of the Gen-Ens-Prod tool in MET version 10.1.0, #1908 removes the ensemble post-processing logic from Ensemble-Stat. With this change, users will likely run Gen-Ens-Prod to generate products, including the ensemble mean, and also run Ensemble-Stat to evaluate the ensemble. The question of handling the ensemble mean arises, since it can be computed by Gen-Ens-Prod but is also needed for computing statistics in Ensemble-Stat. For this issue, consider the two following scenarios and provide solutions for both.

  1. The user runs Gen-Ens-Prod, creates a NetCDF output file including the ensemble mean, and passes that file to Ensemble-Stat using the -ens_mean command line option. As of MET-10.1.0, if -ens_mean is used, it is assumed to be in the same file format at the members. If running Gen-Ens-Prod, that assumption would no longer be valid. Recommend providing enhanced support for the NetCDF output of Gen-Ens-Prod.

    • If -ens_mean specifies a Gen-Ens-Prod output file (need function to determine this), search the NetCDF variables looking for match. For GRIB1/2 inputs, the logic would be pretty easy but if the ensemble members came from NetCDF files originally, it could get messy. Be sure to test both uses.
    • If -ens_mean does NOT specify a Gen-Ens-Prod output file, continue assuming that it follows the same file format/naming conventions as the members provided.
  2. The user runs Ensemble-Stat but does not provide the -ens_mean command line option. Obviously, the ensemble mean can be derived from the input members. Per #1908, if the ensemble_flag.mean = TRUE; write the mean to an output NetCDF file. If not, compute it in memory and use it... but DO NOT write it to an output NetCDF file.

Recommend working the @michelleharrold and the DTC UFS/SRW team to test this functionality.

Time Estimate

2 days.

Sub-Issues

Consider breaking the enhancement down into sub-issues. No sub-issues needed.

Relevant Deadlines

List relevant project deadlines here or state NONE.

Funding Source

Charge 278354 - NOAA R2O RRFS

Define the Metadata

Assignee

Labels

Projects and Milestone

Define Related Issue(s)

Consider the impact to the other METplus components.

Enhancement Checklist

See the METplus Workflow for details.

JohnHalleyGotway commented 1 year ago

@michelleharrold we were reviewing GitHub issues and ran across this one. Can you please update me on the need for this functionality at this time?

michelleharrold commented 1 year ago

This is nice-to-have but not a need-to-have functionality -- so there is no deliverable or imminent project need in relation to this. Thanks!