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

Update and bugfix to movement of the ice-shelf front #475

Merged
merged 1 commit into from
Oct 27, 2023

Conversation

alex-huth
Copy link

Previously, when an ice-shelf cell transitioned from partially-full to full, it was assigned a "reference thickness" equal to the average thickness of surrounding filled cells. In this PR, this average is instead weighted by the flux from each filled cell into the partially-filled cell. This flux-weighted average is more accurate.

Also fixed the incorrect computational bounds used during advection of the ice-shelf front.

The PR has been tested by running the flow-band ice shelf test for 300 years, where the modeled ice-shelf evolution sufficiently matched the analytical solutions for ice-shelf thickness, velocity, and front advection.

@codecov
Copy link

codecov bot commented Aug 29, 2023

Codecov Report

Merging #475 (e2b716d) into dev/gfdl (f514529) will not change coverage.
The diff coverage is 0.00%.

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

@@            Coverage Diff            @@
##           dev/gfdl     #475   +/-   ##
=========================================
  Coverage     37.54%   37.54%           
=========================================
  Files           270      270           
  Lines         79030    79030           
  Branches      14635    14635           
=========================================
  Hits          29670    29670           
  Misses        43894    43894           
  Partials       5466     5466           
Files Coverage Δ
src/ice_shelf/MOM_ice_shelf_dynamics.F90 0.00% <0.00%> (ø)

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

… a partially-filled cell, add the reference thickness contribution from a neighboring filled cell proportionate to its flux into the partially-filled cell. This is more accurate than simply taking the average thickness of all neighboring filled cells. Also fixed incorrect bounds.
Copy link
Member

@Hallberg-NOAA Hallberg-NOAA left a comment

Choose a reason for hiding this comment

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

I agree that these changes seem sensible.

@Hallberg-NOAA
Copy link
Member

This PR has passed pipeline testing at https://gitlab.gfdl.noaa.gov/ogrp/MOM6/-/pipelines/21119.

@Hallberg-NOAA Hallberg-NOAA merged commit 503a9f4 into NOAA-GFDL:dev/gfdl Oct 27, 2023
10 checks passed
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.

2 participants