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

Icepack area underflow patch #184

Merged

Conversation

MJHarrison-GFDL
Copy link
Contributor

@MJHarrison-GFDL MJHarrison-GFDL commented Sep 29, 2022

New parameter RIDGE_AREA_UNDERFLOW is introduced in order to avoid the occurrence of undetectably small fractional area and negative ice volume following ridging adjustments from Icepack. The default is 0.0, which does not introduce changes in existing runs. A reasonable value would be 10^-26 which would correspond to an Angstrom-scale ice patch within a km-scale grid cell. This addresses #178 and #183

…ero, which

preserves previous behavior which resulted in extremely small negative ice volumes in
subsequent calls to ice_continuity, resulting in model failure. This parameter can be
set to a reasonable value (e.g. 10^-30) to avoid this scenario.
@MJHarrison-GFDL MJHarrison-GFDL changed the title Icepack compile message Icepack area underflow patch Sep 29, 2022
Copy link

@MitchBushuk MitchBushuk left a comment

Choose a reason for hiding this comment

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

Thanks @MJHarrison-GFDL for putting this ice area underflow fix in.

src/ice_ridge.F90 Outdated Show resolved Hide resolved
@marshallward
Copy link
Member

Gaea regression: https://gitlab.gfdl.noaa.gov/ogrp/SIS2/-/pipelines/16960 ✔️

I will squash-merge this PR.

@kshedstrom
Copy link
Contributor

I'm afraid you have to check for negative ice volume as well. I get:

WARNING from PE    63: Negative ice thickness at:    233   235     7 -6.2085E-61

FATAL from PE    63: Negative thickness input to ice_continuity().

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.

4 participants