-
Notifications
You must be signed in to change notification settings - Fork 62
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
Ice-shelf FMAs #729
Ice-shelf FMAs #729
Conversation
Tested in ice-ocean mode and ice-shelf-only mode using ISOMIP and MISOMIP+, respectively
…s. Note that the changes here are to the MOM_ice_shelf_dynamics.F90 function 'quad_area', which is only called by the subroutine 'bilinear_shape_functions'; this subroutine is currently unused so that rotationally-consistent use of FMAs in both 'quad_area' and 'bilinear_shape_functions' has not been tested
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.
These changes all look correct to me, and I agree with them.
It should be noted that the mathematically equivalent changes on lines 3361 to 3364 of MOM_ice_shelf_dynamics.F90 are likely to change answers for cases that use a dynamic ice sheet, regardless of whether FMAs are being used.
Gaea regression: https://gitlab.gfdl.noaa.gov/ogrp/mom6ci/MOM6/-/pipelines/25072 ✔️ |
* Fixes to preserve rotational symmetry when using Fused-multiply-adds Tested in ice-ocean mode and ice-shelf-only mode using ISOMIP and MISOMIP+, respectively * Added more parentheses to preserve rotational symmetry when using FMAs. Note that the changes here are to the MOM_ice_shelf_dynamics.F90 function 'quad_area', which is only called by the subroutine 'bilinear_shape_functions'; this subroutine is currently unused so that rotationally-consistent use of FMAs in both 'quad_area' and 'bilinear_shape_functions' has not been tested
Fixes to preserve rotational symmetry when using Fused-multiply-adds. Tested in ice-ocean mode and ice-shelf-only mode using ISOMIP and MISOMIP+, respectively