Skip to content

WIP: Add support for Omega's PseudoThickness variable#440

Draft
xylar wants to merge 12 commits intoE3SM-Project:mainfrom
xylar:support-omega-pseudo-thickness
Draft

WIP: Add support for Omega's PseudoThickness variable#440
xylar wants to merge 12 commits intoE3SM-Project:mainfrom
xylar:support-omega-pseudo-thickness

Conversation

@xylar
Copy link
Copy Markdown
Collaborator

@xylar xylar commented Dec 19, 2025

Checklist

  • User's Guide has been updated
  • Developer's Guide has been updated
  • API documentation in the Developer's Guide (api.md) has any new or modified class, method and/or functions listed
  • Documentation has been built locally and changes look as expected
  • Testing comment in the PR documents testing used to verify the changes
  • New tests have been added to a test suite

@xylar xylar self-assigned this Dec 19, 2025
@xylar xylar added ocean Related to the ocean component utility Utility script(s) in the utils directory Omega PR required The polaris changes won't work with the current Omega submodule and require an update labels Dec 19, 2025
@xylar
Copy link
Copy Markdown
Collaborator Author

xylar commented Feb 18, 2026

Much of this work is now included in #464

@xylar xylar force-pushed the support-omega-pseudo-thickness branch from b465575 to 8b0bc09 Compare April 27, 2026 14:06
Comment thread utils/omega/ctest/omega_ctest.py Outdated
Comment on lines +131 to +133
'ocean.QU.240km.omega_vars.251219.nc',
'PlanarPeriodic48x48.omega_vars.251219.nc',
'cosine_bell_icos480.omega_vars.250219.nc',
Copy link
Copy Markdown
Collaborator Author

Choose a reason for hiding this comment

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

Need to be updated again.

@xylar xylar force-pushed the support-omega-pseudo-thickness branch 2 times, most recently from d225583 to 7267693 Compare May 4, 2026 11:21
xylar added 9 commits May 5, 2026 13:50
This is needed so we can compute the reference pseudo-thickness
for p* from `restingThickness`.

This merge also adds support for multiple iterations when computing
pseudo-thickness using TEOS-10.
Switch layerThickness --> GeometricLayerThickness

Add a few vertical coordinate variables
The new meshes have as many variables as possible mapped to
the new Omega names, and include PseudoThickness (if they
previously had layerThickness) as well as GeometricLayerThickness.
The PseudoThickness is computed using TEOS-10 for the specific
volume.
We no longer assume that `LayerThickness` exists (and that it
is pseudo-thickness).  Instead, we explicitly compute only
`PseudoThickness` and `SpecVol` when writing a model dataset
and compute `GeomLayerThickness` if needed when reading it back.
This is in anticipation of other code besides just the vertical
coordinate diagnostics needing this capability.
@xylar xylar force-pushed the support-omega-pseudo-thickness branch 3 times, most recently from 8b2cadb to 8e99646 Compare May 6, 2026 15:55
xylar added 3 commits May 6, 2026 18:14
Following the renaming of fields, trying to leverage the same
infrastructure is unnecessarily complicated and poteintially
limiting.
@xylar xylar force-pushed the support-omega-pseudo-thickness branch from f4f15be to 782681e Compare May 6, 2026 16:14
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

ocean Related to the ocean component Omega PR required The polaris changes won't work with the current Omega submodule and require an update utility Utility script(s) in the utils directory

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant