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

Floodfill landIceFraction, not landIceMask #848

Merged
merged 1 commit into from
Aug 3, 2024

Conversation

xylar
Copy link
Collaborator

@xylar xylar commented Jul 30, 2024

During the cull_mesh step of creating a global ocean mesh, if ice-shelf cavities are present, we need to flood fill to make sure all land ice is connected (no islands of isolated land ice). This merge changes the flood fill to be on all locations where landIceFraction > 0, not all locations where landIceFraction > 0.5. It then masks all landIce* fields to only be nonzero within the flood-filled region. This prevents islands of nonzero landIceFraction, landIcePressure, landIceDraft, etc.

Checklist

  • Document (in a comment titled Testing in this PR) any testing that was used to verify the changes

@xylar xylar added bug Something isn't working ocean labels Jul 30, 2024
@xylar xylar self-assigned this Jul 30, 2024
@xylar xylar mentioned this pull request Jul 30, 2024
1 task
@xylar xylar force-pushed the fix-land-ice-frac-floodfill branch from 4b66796 to c207ac9 Compare July 31, 2024 13:27
During the `cull_mesh` step of creating a global ocean mesh,
if ice-shelf cavities are present, we need to flood fill to make
sure all land ice is connected (no islands of isolated land ice).
This merge changes the flood fill to be on all locations where
`landIceFraction > 0`, not all locations where `landIceFraction > 0.5`.
It then masks all `landIce*` fields to only be nonzero within the
flood-filled region.  This prevents islands of nonzero
`landIceFraction`, `landIcePressure`, `landIceDraft`, etc.
@xylar xylar force-pushed the fix-land-ice-frac-floodfill branch from c207ac9 to 1d6ef70 Compare July 31, 2024 18:29
@xylar
Copy link
Collaborator Author

xylar commented Aug 1, 2024

Testing

I tested this in conjunction with #807. Without this fix, there were still regions where landIceFraction was nonzero around islands, which I think we probably don't want.

Having nonzero landIceFraciton around islands was also caused confusion since @proteanplanet's sanity checking was using landIceFraction (rather than, e.g., landIceMask or landIceFloatingFraction) to determine the ice-shelf extent.

@xylar xylar merged commit df18be6 into MPAS-Dev:main Aug 3, 2024
4 checks passed
@xylar xylar deleted the fix-land-ice-frac-floodfill branch August 3, 2024 16:17
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working ocean
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant