Skip to content

Static ice shelves with NUOPC - #112

Draft
angus-g wants to merge 5 commits into
2026.05from
angus-g/ice-shelf-update
Draft

angus-g wants to merge 5 commits into
2026.05from
angus-g/ice-shelf-update

Conversation

@angus-g

@angus-g angus-g commented Oct 25, 2026 •

Copy link
Copy Markdown
Collaborator

This is a minimal set of cherry-picked commits from #38, excluding a bespoke flux override and commits that will already come in through GFDL.

#38 also had a modification to support (dynamic) shelf restarts with NUOPC. When #110 lands, this is no longer required for static shelves since they have no restartable fields. In general, this will need to be addressed for dynamic shelves with NUOPC.

claireyung and others added 5 commits October 23, 2026 15:07
- turn FATAL error when ESMF mask and MOM6 mask differ intol WARNING, so that we can set up a MOM6 domain including ice shelves but avoid ESMF and the mediator transferrring atmospheric data into the ice shelves
- use ocean_state%fluxes%frac_shelf_h when computing the So_omask for the mediator from MOM6 data, so that the coupler does not get sent fields in regions it does not expect. A ceiling is used for frac_shelf_h (fraction of ice shelf values) - partial ice shelf coverage may need to be better thought out.
- for above to work, the ocean_state_type is made public.
- Also initialise ice shelf fluxes and forces separately. Otherwise, shelf_sfc_mass_flux field is initialised which should only be used with dynamic ice shelf (and results in crash when it can't find the file if you're running from a restart.)
This parameter scales the friction velocity computed due to prescribed tides
either through UTIDE = constant value or through TIDEAMP = True and associated
TIDEAMP_FILE. Jourdain et al. 2019 suggest a value of 0.66 to estimate tidal ice shelf
boundary layer velocity from total barotropic tide amplitude. Default is set to 1 so
should not change answers
Parameters READ_TIDEAMP and UTIDE not only affect INT_TIDE_DISSIPATION but also ice shelf melt rates
via the ice shelf basal melt rate friction velocity calculation. In this commit the effect on
ice shelves is added to the description of the parameters to make this clearer.
@angus-g
angus-g requested a review from dougiesquire October 25, 2026 22:35
@angus-g
angus-g marked this pull request as draft October 25, 2026 22:42
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