-
Notifications
You must be signed in to change notification settings - Fork 375
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
Fix split-explicit upwinded advection #6047
Fix split-explicit upwinded advection #6047
Conversation
TestingThis PR introduces machine-precision-level changes in the tests included in the compass PR suite, which was run on Chrysalis with intel, openmpi. This PR was tested with upwinding on for the |
Tested with gnu on chrysalis with both debug and optimized. Passes nightly suite with default flags and matches master bfb as expected. I also added these flags to compass:
and this branch passes nightly suite with both gnu debug and optimized. When I compare to master using upwind, I get small differences, as expected: diffs:
|
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 also reviewed visually. Based on testing, I am approving. Note that this does not affect E3SM mainline simulations.
@jonbob I think this should be labeled 'BFB' as it pertains to tested E3SM simulations. We could even add |
Thanks for the testing @mark-petersen. I'll change the label but test in E3SM as well before merging |
@mark-petersen Thank you for your review and testing! |
notes: ready but low priority |
…6047) Fix split-explicit upwinded advection This PR introduces a new variable to track the baroclinic water column thickness for the barotropic subcycles. The method of storing the baroclinic water column thickness in a bottomDepthEdge variable was problematic for upwinding because the upwinded cell can change between the baroclinic cycle and the barotropic subcycle. We now store the baroclinic (upwinded or centered) water column thickness at edges and then apply it along with the barotropic (upwinded or centered) ssh on edges to determine the barotropic water column thickness at edges. Fixes #6048 [BFB]
passes:
merged to next |
merged to master |
This PR introduces a new variable to track the baroclinic water column thickness for the barotropic subcycles. The method of storing the baroclinic water column thickness in a
bottomDepthEdge
variable was problematic for upwinding because the upwinded cell can change between the baroclinic cycle and the barotropic subcycle. We now store the baroclinic (upwinded or centered) water column thickness at edges and then apply it along with the barotropic (upwinded or centered) ssh on edges to determine the barotropic water column thickness at edges.Fixes #6048
[BFB]