[Feature]: Support for p_levs instead of hyam and hybm - re EAMxx remapped files
Opened this issue · 1 comments
Is your feature request related to a problem?
EAMxx introduced a new capability which allows the user to remap output onto a static pressure level vertical coordinate. This capability does not use the hyam
and hybm
hybrid coordinate system because we do not need to build the pressure from the surface pressure, since pressure levels are static in the data. E3SM Diags assumes that hyam
and hybm
will be in the source data and as a result is unable to produce any pressure based diagnostics for this remapped data.
Describe the solution you'd like
An option to use p_levs
(which is the fixed set of pressure levels in the data) rather than build on the fly from hyam
and hybm
would be great. The solution could also be extended to use the p_mid
and p_int
variables to determine the pressure levels of data rather than rely on hyam
and hybm
.
Describe alternatives you've considered
So far I haven't tried alternative approaches, although I may try "hacking" an hyam and hybm with a hacked PS
to ensure that e3sm-diags always builds the appropriate pressure levels for the data.
Additional context
No response
@AaronDonahue thanks for filing this issue! I just got back from vacation and am catching up. Actually e3sm_diags doesn't request the vertical coordinate to be hybrid, it also supports pressure levels, such as CMIP datasets and those reanalysis datasets. I'm wondering if the program can't recognize specs of the p_lev
coordinate that EAMxx data has. Could you share the input data and I can probably figure out the code change needed for supporting the file.