-
Notifications
You must be signed in to change notification settings - Fork 388
Fixes limiting on redi k33 #684
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 limiting on redi k33 #684
Conversation
|
@mark-petersen we had originally put k33 before limiting slopes to have stronger diffusion for stability, but testing in E3SM has shown it creates too much mixing and overly deep boundary layers. This PR moves the k33 calculation post limiting to reduce the associated vertical mixing. I have not run into any stability issues with this change (30 years of G-case and 40 years of B-case). |
|
I will run the nightly soon and post results. |
on the vertical redi term, the surface scaling is currently being applied twice. This PR removes that. It also reorganizes the code to increase the strength of slope limiting on k33.
|
Rebased, tested on grizzly with gnu and intel 19, debug and optimized. Passes nightly regression suite. |
30d2c0f to
6b92c01
Compare
mark-petersen
left a comment
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.
@vanroekel please check this last commit. Otherwise, this is ready to merge.
|
@mark-petersen looks good. Thanks for taking out the extra calculation. |
Fixes limiting on redi k33 #684 on the vertical redi term, the surface scaling is currently being applied twice. This PR removes that. It also reorganizes the code to increase the strength of slope limiting on k33.
This PR brings in a new mpas-source submodule with changes only to the ocean core, plus scripts updates corresponding to Registry changes in the ocean. It includes three changes to improve ocean heat uptake: * Fixes limiting on redi k33 (MPAS-Dev/MPAS-Model/pull/684); * Fixes surface buoyancy forcing calculation (MPAS-Dev/MPAS-Model/pull/690); and * Revert GM buoyancy gradient calculation to v1 form (MPAS-Dev/MPAS-Model/pull/687). [NML] [non-BFB]
Update mpas-source: GM/Redi and surface buoyancy fixes This PR brings in a new mpas-source submodule with changes only to the ocean core, plus scripts updates corresponding to Registry changes in the ocean. It includes three changes to improve ocean heat uptake: * Fixes limiting on redi k33 (MPAS-Dev/MPAS-Model/pull/684); * Fixes surface buoyancy forcing calculation (MPAS-Dev/MPAS-Model/pull/690); and * Revert GM buoyancy gradient calculation to v1 form (MPAS-Dev/MPAS-Model/pull/687). [NML] [non-BFB]
on the vertical redi term, the surface scaling is currently being
applied twice. This PR removes that. It also reorganizes the code to
increase the strength of slope limiting on k33.