-
Notifications
You must be signed in to change notification settings - Fork 66
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Updates MOM_kappa_shear averaging and diagnostics. #846
Updates MOM_kappa_shear averaging and diagnostics. #846
Conversation
breichl
commented
Mar 3, 2025
- Update to horizontal averaging procedure in MOM_kappa_shear for VERTEX_SHEAR = True.
- Adds geometric and thickness weighted averages for moving diffusivity from vertices to tracer points.
- Adds new diagnostics of N2, S2, Kd, and TKE both with and without VERTEX_SHEAR, useful for debugging.
- Bitwise reproduces OM5 config w/ VERTEX_SHEAR = True.
- Update to horizontal averaging procedure in MOM_kappa_shear for VERTEX_SHEAR = True. - Adds geometric and thickness weighted averages for moving diffusivity from vertices to tracer points. - Adds new diagnostics of N2, S2, Kd, and TKE both with and without VERTEX_SHEAR, useful for debugging.
ed026dc
to
c79c5f6
Compare
- Fix OMP directives for new diagnostics - Remove quasi-2d fields that broke OMP threading - Reworking horizontal averaging algorithm.
The code now passes all OMP tests. I've moved the PR to draft until I finish evaluating that the new features still work as intended. |
Revised the Calc_kappa_shear_vertex code where it averages the diffusivities at the vertices back to the tracer points when there is either a geometric mean being used or thickness weighting. With these changes, the code should now pass dimensional consistency testing, work properly with openMP threading enabled, and give sensible values of the diffusivity for massless layers. Also corrected the indenting of several recently added lines to follow the MOM6 2-space indenting convention. This commit can change answers slightly in some cases where VERTEX_SHEAR = True, especially if VERTEX_SHEAR_GEOMETRIC_MEAN = True or VERTEX_SHEAR_THICKNESS_MEAN = True. The recently added runtime parameter VERTEX_SHEAR_KD_MEAN_ANSWER_DATE (which has not yet been added to the dev/gfdl branch of the code) was no longer needed and has been eliminated.
*+Revise horizontal averaging in Calc_kappa_shear_vertex
…help with line length issues.
Codecov ReportAttention: Patch coverage is
Additional details and impacted files@@ Coverage Diff @@
## dev/gfdl #846 +/- ##
============================================
- Coverage 38.07% 38.05% -0.02%
============================================
Files 298 298
Lines 87830 87955 +125
Branches 16467 16516 +49
============================================
+ Hits 33440 33470 +30
- Misses 48352 48425 +73
- Partials 6038 6060 +22 ☔ View full report in Codecov by Sentry. 🚀 New features to boost your workflow:
|
I did a layout test in OM5 (with VERTEX_SHEAR = True, VERTEX_SHEAR_GEOMETRIC_MEAN = True, and VERTEX_SHEAR_THICKNESS_MEAN = True), which reproduced. I also checked that the effects of turning the new averaging on in a 1 month run are consistent with my developmental code and tests. I think this is now ready. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I have carefully examined these changes, and agree that they are correct.
This PR has passed pipeline testing at https://gitlab.gfdl.noaa.gov/ogrp/mom6ci/MOM6/-/pipelines/26630 with the expected warnings about new runtime parameters and additional diagnostics. |