Skip to content
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

Fixes MPAS-Ocean Richardson number calculation #5946

Conversation

vanroekel
Copy link
Contributor

@vanroekel vanroekel commented Sep 21, 2023

This changes the calculation of the gradient and bulk richardson numbers to utilize the normal and tangential velocities to avoid the factor of two issue and also fixes a number of places where the loop was over maxLevelCell and should have been to maxLevelEdgeTop.

Testing is documented in -- E3SM-Ocean-Discussion#66

[NML]
[CC]

vanroekel and others added 6 commits September 21, 2023 16:22
This changes the calculation of the gradient and bulk richardson numbers
to utilize the normal and tangential velocities to avoid the factor of
two issue and also fixes a number of places where the loop was over
maxLevelCell and should have been to maxLevelEdgeTop
@mark-petersen
Copy link
Contributor

Three ocean PRs were merged locally and pass the nightly test suite, as described in #5945 (comment).

Copy link
Contributor

@mark-petersen mark-petersen left a comment

Choose a reason for hiding this comment

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

I carefully reviewed the code and the simulation results. This is a straightforward substitution, from delU2=normalVelocity**2 to delU2=normalVelocity**2 + tangentialVelocity**2. I agree with these changes.

@rljacob rljacob added this to the v3.0alpha05 milestone Oct 12, 2023
jonbob added a commit that referenced this pull request Oct 24, 2023
… next (PR #5946)

Fixes MPAS-Ocean Richardson number calculation

This changes the calculation of the gradient and bulk richardson numbers
to utilize the normal and tangential velocities to avoid the factor of
two issue and also fixes a number of places where the loop was over
maxLevelCell and should have been to maxLevelEdgeTop.

[CC]
@jonbob
Copy link
Contributor

jonbob commented Oct 24, 2023

merged to next

jonbob added a commit that referenced this pull request Oct 25, 2023
@jonbob
Copy link
Contributor

jonbob commented Oct 25, 2023

Re-merged to next after fix for threading issues

@jonbob
Copy link
Contributor

jonbob commented Oct 26, 2023

@philipwjones -- if you get a chance, could you please check the acc update I made in the last commit?

@jonbob jonbob merged commit 510fca6 into E3SM-Project:master Oct 26, 2023
2 checks passed
@jonbob
Copy link
Contributor

jonbob commented Oct 26, 2023

merged to master

xylar added a commit to xylar/compass that referenced this pull request Dec 3, 2023
This merge updates the E3SM-Project submodule from [894b5b2](https://github.com/E3SM-Project/E3SM/tree/894b5b2) to [5d5f15c](https://github.com/E3SM-Project/E3SM/tree/5d5f15c).

This update includes the following MPAS-Ocean and MPAS-Frameworks PRs (check mark indicates bit-for-bit with previous PR in the list):
- [ ]  (ocn) E3SM-Project/E3SM#5945
- [ ]  (ocn) E3SM-Project/E3SM#5946
- [ ]  (ocn) E3SM-Project/E3SM#5947
- [ ]  (ocn) E3SM-Project/E3SM#5999
- [ ]  (ocn) E3SM-Project/E3SM#6037
- [ ]  (ocn) E3SM-Project/E3SM#5989
- [ ]  (ocn) E3SM-Project/E3SM#6035
- [ ]  (ocn) E3SM-Project/E3SM#6077
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
CC PR is climate changing mpas-ocean NML
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants