Skip to content

Added frazil to ice shelf#985

Merged
Hallberg-NOAA merged 2 commits into
NOAA-GFDL:dev/gfdlfrom
alex-huth:frazil_pr
Dec 2, 2025
Merged

Added frazil to ice shelf#985
Hallberg-NOAA merged 2 commits into
NOAA-GFDL:dev/gfdlfrom
alex-huth:frazil_pr

Conversation

@alex-huth
Copy link
Copy Markdown

The frazil mass flux to the ice-shelf base is calculated by multiplying frazil energy [J m-2] by the inverse of the timestep times the latent heat of fusion [kg J-1 s-1].

This frazil mass flux is incorporated as a negative water flux from the ice shelf, which then acts to add the frazil mass to the ice shelf base
(MOM_ice_shelf.F90/change_thickness_using_melt) and remove it from the ocean surface as evaporation (MOM_ice_shelf.F90/add_shelf_flux).

Note frazil is reset to zero at the start of each therm timestep in MOM.F90/step_MOM.

Some additional changes were also made to how the ice-shelf flux factor is implemented, so that is only scales ice-shelf melt without affecting the frazil mass flux.

Comment thread src/ice_shelf/MOM_ice_shelf.F90
Comment thread src/ice_shelf/MOM_ice_shelf.F90 Outdated
Copy link
Copy Markdown
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 am convinced that this code is correct as written, and should correctly address the issue of how to handle the frazil heat fluxes under ice shelves.

At a later date, we may need to revisit how frazil is handled when there is partial coverage of the ocean by ice shelves, especially when that coverage is evolving with time.

The frazil mass flux to the ice-shelf base is calculated by
multiplying frazil energy [J m-2] by the inverse of the timestep times
the latent heat of fusion [kg J-1 s-1].

This frazil mass flux is incorporated as a negative water flux from
the ice shelf. This negative water flux then acts to add the frazil
mass to the ice shelf base
(MOM_ice_shelf.F90/change_thickness_using_melt) and remove it from
the ocean surface as evaporation (MOM_ice_shelf.F90/add_shelf_flux).

Note frazil is reset to zero at the start of each therm timestep in
MOM.F90/step_MOM.

Some additional changes were also made to how the ice-shelf flux
factor is implemented, so that is only scales ice-shelf melt without
affecting the frazil mass flux.
@Hallberg-NOAA
Copy link
Copy Markdown
Member

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

@Hallberg-NOAA Hallberg-NOAA merged commit 96e186f into NOAA-GFDL:dev/gfdl Dec 2, 2025
52 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