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

(*)Offset G%bathyT by -G%Z_ref #1490

Merged
merged 2 commits into from
Sep 11, 2021

Conversation

Hallberg-NOAA
Copy link
Collaborator

Added code to adjust the value of G%bathyT in the ocean_grid_type by -G%Z_ref,
with changes scattered throughout the code to compensate. This does not impact
the value of bathyT in the dyn_hor_grid. Using a non-zero value of
REFERENCE_HEIGHT leads to mathematically equivalent solutions that differ at
round-off but are demonstrably similar for values of REFERENCE_HEIGHT ranging
from -1e4 m to 1e4 m, and in cases that do not use the less exact
..._2018_ANSWERS algorithms the solutions are qualitatively equivalent for
values ranging from -1e8 m to 1e8 m. (For these more extreme values, the 64-bit
roundoff in the free surface height calculation starts to become physically
significant at of order 0.1 mm.) This new capability is useful for developing
wetting and drying and for identifying algorithmic shortcomings, but because
answers are not bitwise identical for various reference heights, it is not such
a good candidate for fully automated testing. By default, all answers are
bitwise identical and all output files are the same.

  Added code to adjust the value of G%bathyT in the ocean_grid_type by -G%Z_ref,
with changes scattered throughout the code to compensate.  This does not impact
the value of bathyT in the dyn_hor_grid.  Using a non-zero value of
REFERENCE_HEIGHT leads to mathematically equivalent solutions that differ at
round-off but are demonstrably similar for values of REFERENCE_HEIGHT ranging
from -1e4 m to 1e4 m, and in cases that do not use the less exact
..._2018_ANSWERS algorithms the solutions are qualitatively equivalent for
values ranging from -1e8 m to 1e8 m.  (For these more extreme values, the 64-bit
roundoff in the free surface height calculation starts to become physically
significant at of order 0.1 mm.)  This new capability is useful for developing
wetting and drying and for identifying algorithmic shortcomings, but because
answers are not bitwise identical for various reference heights, it is not such
a good candidate for fully automated testing.  By default, all answers are
bitwise identical and all output files are the same.
@codecov
Copy link

codecov bot commented Aug 31, 2021

Codecov Report

Merging #1490 (f44b9c4) into dev/gfdl (a316416) will increase coverage by 0.00%.
The diff coverage is 9.39%.

Impacted file tree graph

@@            Coverage Diff            @@
##           dev/gfdl    #1490   +/-   ##
=========================================
  Coverage     29.01%   29.02%           
=========================================
  Files           236      236           
  Lines         71506    71515    +9     
=========================================
+ Hits          20750    20759    +9     
  Misses        50756    50756           
Impacted Files Coverage Δ
src/ALE/MOM_regridding.F90 21.45% <0.00%> (ø)
src/ALE/coord_adapt.F90 0.00% <0.00%> (ø)
src/core/MOM_PressureForce_Montgomery.F90 6.36% <0.00%> (-0.02%) ⬇️
src/core/MOM_grid.F90 61.68% <0.00%> (ø)
src/core/MOM_open_boundary.F90 20.00% <ø> (ø)
src/core/MOM_transcribe_grid.F90 24.63% <0.00%> (ø)
src/diagnostics/MOM_PointAccel.F90 3.08% <0.00%> (ø)
src/diagnostics/MOM_sum_output.F90 60.64% <0.00%> (ø)
src/diagnostics/MOM_wave_speed.F90 23.68% <0.00%> (ø)
src/framework/MOM_diag_remap.F90 79.57% <0.00%> (ø)
... and 23 more

Continue to review full report at Codecov.

Legend - Click here to learn more
Δ = absolute <relative> (impact), ø = not affected, ? = missing data
Powered by Codecov. Last update a316416...f44b9c4. Read the comment docs.

@marshallward
Copy link
Collaborator

Gaea regression: https://gitlab.gfdl.noaa.gov/ogrp/MOM6/-/pipelines/13555 ✔️

This did eventually pass, despite many snags and reported failures during testing. I believe they were all related to Lustre problems and hopefully do not reflect any problems with this PR or any other recently merged PRs, such as the rotational read_data support. But let's keep an eye on it.

@marshallward marshallward merged commit 1623085 into mom-ocean:dev/gfdl Sep 11, 2021
@Hallberg-NOAA Hallberg-NOAA deleted the offset_bathyT_by_Z_ref branch November 22, 2021 09:07
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants