Skip to content

Conversation

@vanroekel
Copy link
Contributor

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.

@vanroekel
Copy link
Contributor Author

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

@vanroekel
Copy link
Contributor Author

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.
@mark-petersen
Copy link
Contributor

Rebased, tested on grizzly with gnu and intel 19, debug and optimized. Passes nightly regression suite.

@mark-petersen mark-petersen self-assigned this Sep 9, 2020
@mark-petersen mark-petersen force-pushed the ocean/alterRediK33Limiting branch from 30d2c0f to 6b92c01 Compare September 9, 2020 20:44
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.

@vanroekel please check this last commit. Otherwise, this is ready to merge.

@vanroekel
Copy link
Contributor Author

@mark-petersen looks good. Thanks for taking out the extra calculation.

mark-petersen added a commit that referenced this pull request Sep 10, 2020
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.
mark-petersen added a commit that referenced this pull request Sep 10, 2020
Fixes limiting on redi k33 #684
Fixes surface buoyancy forcing calculation #690
Revert GM buoyancy gradient calculation to v1 form #687
@mark-petersen mark-petersen merged commit aa900ab into MPAS-Dev:ocean/develop Sep 11, 2020
mark-petersen added a commit that referenced this pull request Sep 11, 2020
Fixes limiting on redi k33 #684
Fixes surface buoyancy forcing calculation #690
Revert GM buoyancy gradient calculation to v1 form #687
jonbob added a commit to E3SM-Project/E3SM that referenced this pull request Sep 14, 2020
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]
jonbob added a commit to E3SM-Project/E3SM that referenced this pull request Sep 15, 2020
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]
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants