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

+*Non-Boussinesq revision of MOM_sponge code #441

Merged
merged 2 commits into from
Aug 23, 2023

Conversation

Hallberg-NOAA
Copy link
Member

Use tv%SpV convert thicknesses to vertical distances in apply_sponge when it is allocated to replace multiplication by GV%H_to_Z, thereby eliminating the dependence on GV%RHo_0 in this modue in non-Boussinesq mode. The new internal array dz_to_h is used to reduce the code duplication as a result of these changes. All answers in Boussinesq test cases are bitwise identical, but answers change in fully non-Boussinesq mode. In semi-Boussinesq mode answers are mathematically equivalent but change at roundoff unless RHO_0 is an integer power of 2.

  Use tv%SpV convert thicknesses to vertical distances in apply_sponge when it
is allocated to replace multiplication by GV%H_to_Z, thereby eliminating the
dependence on GV%RHo_0 in this modue in non-Boussinesq mode.  The new internal
array dz_to_h is used to reduce the code duplication as a result of these
changes.  All answers in Boussinesq test cases are bitwise identical, but
answers change in fully non-Boussinesq mode.  In semi-Boussinesq mode answers
are mathematically equivalent but change at roundoff unless RHO_0 is an integer
power of 2.
@codecov
Copy link

codecov bot commented Aug 6, 2023

Codecov Report

Merging #441 (e96ac0f) into dev/gfdl (d60c2e0) will decrease coverage by 0.02%.
The diff coverage is 0.00%.

❗ Current head e96ac0f differs from pull request most recent head 47b109a. Consider uploading reports for the commit 47b109a to get more accurate results

@@             Coverage Diff              @@
##           dev/gfdl     #441      +/-   ##
============================================
- Coverage     38.03%   38.01%   -0.02%     
============================================
  Files           269      269              
  Lines         77554    77576      +22     
  Branches      14319    14323       +4     
============================================
  Hits          29494    29494              
- Misses        42709    42731      +22     
  Partials       5351     5351              
Files Changed Coverage Δ
src/parameterizations/vertical/MOM_sponge.F90 0.00% <0.00%> (ø)

📣 We’re building smart automated test selection to slash your CI/CD build times. Learn more

@Hallberg-NOAA Hallberg-NOAA added enhancement New feature or request answer-changing A change in results (actual or potential) labels Aug 6, 2023
@Hallberg-NOAA Hallberg-NOAA changed the title +(*)Use tv%SpV in MOM_sponge code +*Non-Boussinesq revision of MOM_sponge code Aug 6, 2023
Copy link
Member

@marshallward marshallward left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Code appears to preserve the existing Boussinesq solution. To the best of my limited knowledge, the nonBoussinesq ones using specific volume seem correct.

@marshallward
Copy link
Member

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

@marshallward marshallward merged commit bd5fe0c into NOAA-GFDL:dev/gfdl Aug 23, 2023
10 checks passed
@Hallberg-NOAA Hallberg-NOAA deleted the nonBous_sponge branch September 27, 2023 16:17
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
answer-changing A change in results (actual or potential) enhancement New feature or request
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants