Skip to content

Update ice_area to handle NaNs#237

Merged
pgierz merged 3 commits into
FESOM:mainfrom
Mukulikaa:sia
May 12, 2026
Merged

Update ice_area to handle NaNs#237
pgierz merged 3 commits into
FESOM:mainfrom
Mukulikaa:sia

Conversation

@Mukulikaa
Copy link
Copy Markdown
Contributor

Changes copied from updates to ice_ext. Tested with AWI-ESM data which shows ~17 Mio km2 sea ice area for the Southern Ocean.

Comment thread pyfesom2/diagnostics.py
def ice_area(data, mesh, hemisphere="N", attrs={}):
""" Compute sea ice volume.
def ice_area(data, mesh, hemisphere="N", threshhold=0.15, attrs={}):
""" Compute sea ice area.
Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

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

I am not sure it's a good idea to put threshhold=0.15 - this is something we use for ice_extent, but definition of the area it the total area covered by sea ice, and at leas in my mind this is what I would want to see by default. Adding threshold is totally fine, but I would put 0 as a default.

@pgierz
Copy link
Copy Markdown
Member

pgierz commented May 12, 2026

I've enabled the CI, since this was an outside fork. Apart from @koldunovn comment about the threshold, this looks fine to me.

@Mukulikaa
Copy link
Copy Markdown
Contributor Author

Thanks, updated the branch.

@pgierz pgierz self-requested a review May 12, 2026 11:56
@pgierz
Copy link
Copy Markdown
Member

pgierz commented May 12, 2026

Alright, all tests pass, looks good. I am merging, nice work! :)

@pgierz pgierz merged commit 53f2561 into FESOM:main May 12, 2026
15 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.

3 participants