Add geopotential height at 700 hPa diagnostic (processed__geopot)#605
Add geopotential height at 700 hPa diagnostic (processed__geopot)#605Zubair Maalick (zmaalick) wants to merge 1 commit into
Conversation
Implements issue MetOffice#604: geopotential height at 700 hPa (zg700), output 6-hourly. Reuses the already-computed multi-level plev__geopot by slicing its 700 hPa level (index 5 of the pressure_levels axis) via an XIOS axis zoom, mirroring the existing processed__pstar two-step extraction pattern. - grid_def_main.xml: new pressure_levels_700_face_zoom grid. - field_def_diags.xml: zoomed__geopot (level slice) and processed__geopot (face field, name=zg700, standard_name=geopotential_height, unit=m). - file_def_diags_oper_nwp_gl.xml: new 6-hourly file group requesting it. Registry-only change; no Fortran/kernel/algorithm edits. Co-Authored-By: Claude Opus 4.8 <noreply@anthropic.com>
|
Hello Zubair Maalick (@zmaalick)! 👋 Thank you for your contribution. Since this is your first time contributing to this repository, we ask that you sign our Contributor Licence Agreement (CLA). To agree to the CLA, please add your details (GitHub username, Real Name, Affiliation, and Date) to the CONTRIBUTORS.md file (create one, if required) in the development branch for this PR. After signing the CLA, you won't need to do this again for future PRs. |
iboutle
left a comment
There was a problem hiding this comment.
Hi Zubair,
I'm a bit puzzled for what the motivation for this diagnostic is - are you able to give any further details about why we need a specific diagnostic of the 700hPa geopotential that is separate from the general geopotential height on pressure levels diagnostic already available?
In particular, you're adding this diagnostic to the operational suite diagnostics list, but I can't see anything equivalent in the operational UM diagnostics list?
As an aside, I don't think this is the correct way to produce this diagnostic if indeed we do need it (we can't assume that 700hPa will always be level 5, because the pressure level axis is user defined), but it will be easier to establish how to approach it when I understand the requirement better.
Thanks
Ian
Implements issue #604: geopotential height at 700 hPa (zg700), output 6-hourly. Reuses the already-computed multi-level plev__geopot by slicing its 700 hPa level (index 5 of the pressure_levels axis) via an XIOS axis zoom, mirroring the existing processed__pstar two-step extraction pattern.
Registry-only change; no Fortran/kernel/algorithm edits.
PR Summary
Sci/Tech Reviewer:
Code Reviewer: Benjamin Went (@MetBenjaminWent)
Code Quality Checklist
Testing
trac.log
Security Considerations
Performance Impact
AI Assistance and Attribution
Documentation
PSyclone Approval
Sci/Tech Review
(Please alert the code reviewer via a tag when you have approved the SR)
Code Review