Skip to content

Add geopotential height at 700 hPa diagnostic (processed__geopot)#605

Open
Zubair Maalick (zmaalick) wants to merge 1 commit into
MetOffice:mainfrom
zmaalick:604_diagnostics_geopot
Open

Add geopotential height at 700 hPa diagnostic (processed__geopot)#605
Zubair Maalick (zmaalick) wants to merge 1 commit into
MetOffice:mainfrom
zmaalick:604_diagnostics_geopot

Conversation

@zmaalick

@zmaalick Zubair Maalick (zmaalick) commented Jul 8, 2026

Copy link
Copy Markdown

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.

  • 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.

PR Summary

Sci/Tech Reviewer:
Code Reviewer: Benjamin Went (@MetBenjaminWent)

Code Quality Checklist

  • I have performed a self-review of my own code
  • My code follows the project's style guidelines
  • Comments have been included that aid understanding and enhance the readability of the code
  • My changes generate no new warnings
  • All automated checks in the CI pipeline have completed successfully

Testing

  • I have tested this change locally, using the LFRic Apps rose-stem suite
  • If any tests fail (rose-stem or CI) the reason is understood and acceptable (e.g. kgo changes)
  • I have added tests to cover new functionality as appropriate (e.g. system tests, unit tests, etc.)
  • Any new tests have been assigned an appropriate amount of compute resource and have been allocated to an appropriate testing group (i.e. the developer tests are for jobs which use a small amount of compute resource and complete in a matter of minutes)

trac.log

Security Considerations

  • I have reviewed my changes for potential security issues
  • Sensitive data is properly handled (if applicable)
  • Authentication and authorisation are properly implemented (if applicable)

Performance Impact

  • Performance of the code has been considered and, if applicable, suitable performance measurements have been conducted

AI Assistance and Attribution

  • Some of the content of this change has been produced with the assistance of Generative AI tool name (e.g., Met Office Github Copilot Enterprise, Github Copilot Personal, ChatGPT GPT-4, etc) and I have followed the Simulation Systems AI policy (including attribution labels)

Documentation

  • Where appropriate I have updated documentation related to this change and confirmed that it builds correctly

PSyclone Approval

  • If you have edited any PSyclone-related code (e.g. PSyKAl-lite, Kernel interface, optimisation scripts, LFRic data structure code) then please contact the TCD Team

Sci/Tech Review

  • I understand this area of code and the changes being added
  • The proposed changes correspond to the pull request description
  • Documentation is sufficient (do documentation papers need updating)
  • Sufficient testing has been completed

(Please alert the code reviewer via a tag when you have approved the SR)

Code Review

  • All dependencies have been resolved
  • Related Issues have been properly linked and addressed
  • CLA compliance has been confirmed
  • Code quality standards have been met
  • Tests are adequate and have passed
  • Documentation is complete and accurate
  • Security considerations have been addressed
  • Performance impact is acceptable

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>
@github-actions github-actions Bot added the cla-required The CLA has not yet been signed by the author of this PR - added by GA label Jul 8, 2026
@github-actions

github-actions Bot commented Jul 8, 2026

Copy link
Copy Markdown

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).

📄 You can read the CLA here.

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 iboutle left a comment

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

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

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

cla-required The CLA has not yet been signed by the author of this PR - added by GA

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants