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

Seaice over land after the analysis insertion #1131

Closed
guillaumevernieres opened this issue Feb 4, 2025 · 7 comments
Closed

Seaice over land after the analysis insertion #1131

guillaumevernieres opened this issue Feb 4, 2025 · 7 comments
Assignees
Labels
bug Something is not working SOCA Sea-ice, Ocean, and Coupled Assimilation

Comments

@guillaumevernieres
Copy link
Contributor

Description

This was mentioned here:
ufs-community/ufs-weather-model#2562 (comment)
But the discussion above is based on a ~2 yo version of the soca to cice variable change. It's possible that it has been fixed, but better double check.

To Reproduce:

Run soca2cice and check that we have no sea ice over land.

@guillaumevernieres guillaumevernieres added bug Something is not working SOCA Sea-ice, Ocean, and Coupled Assimilation labels Feb 4, 2025
@DeniseWorthen
Copy link

The other item to check is whether any ice volume exists where there is no associated ice coverage (aicen=0 but vicen /= 0)

@shlyaeva
Copy link
Contributor

shlyaeva commented Feb 4, 2025

I'm hoping this is resolved now by calling icepack's cleanup_itd in initialization, but we'll need to check.

@guillaumevernieres
Copy link
Contributor Author

I'm hoping this is resolved now by calling icepack's cleanup_itd in initialization, but we'll need to check.

🤞

@shlyaeva
Copy link
Contributor

I started working on it, and couldn't make any ice over land in my test as hard as I tried. Looks like someone smart already added code to make sure we're outputting zero over land:

if ( pe == mpp_root_pe()) glb%aicen(1:self%ni,1:self%nj,n) = var2d*glb%iceumask
.

I'll double check that the second concern

The other item to check is whether any ice volume exists where there is no associated ice coverage (aicen=0 but vicen /= 0)

is covered by cleanup_itd before closing this.

@DeniseWorthen
Copy link

@shlyaeva Thanks for the followup. I found a third condition also, which is to check for whether there is snow volume (vsnon) where there is no aicen. The opposite is fine I think---ice can have no snow coverage.

@shlyaeva
Copy link
Contributor

I just confirmed that if I change the code that puts analysis into the CICE restart files and artificially set vsnon and vicen to be non-zero when aicen is zero, cleanup_itd that is called at the end resets those values back to zero. So I think the current code covers all the issues reported here. OK with everyone to close?

@DeniseWorthen
Copy link

From my perspective, it seems like whatever was previously the problem has been fixed in the intervening years. So, all is good.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something is not working SOCA Sea-ice, Ocean, and Coupled Assimilation
Projects
None yet
Development

No branches or pull requests

3 participants