Skip to content

Add sea ice recentering#1479

Merged
guillaumevernieres merged 3 commits into
feature/lowres-bfrom
feature/ice_recenter
Feb 6, 2025
Merged

Add sea ice recentering#1479
guillaumevernieres merged 3 commits into
feature/lowres-bfrom
feature/ice_recenter

Conversation

@shlyaeva
Copy link
Copy Markdown
Collaborator

@shlyaeva shlyaeva commented Feb 6, 2025

Description

Adds recentering for sea ice by adding the recentering increment to the sea ice background, and calling Soca2Cice variable change.

Requires copying background files and ice restarts, I wasn't familiar with how to code this, so please review.

Companion PRs

Needs a change to global-workflow as well: guillaumevernieres/global-workflow#10

Issues

Resolves #1454

Automated CI tests to run in Global Workflow

  • atm_jjob
  • C96C48_ufs_hybatmDA
  • C96C48_hybatmaerosnowDA
  • C48mx500_3DVarAOWCDA
  • C48mx500_hybAOWCDA
  • C96C48_hybatmDA

Comment thread ush/soca/marine_recenter.py Outdated
Comment on lines +187 to +192
fname_in = os.path.abspath(os.path.join(mem_dir_real, f00))
fname_out = os.path.realpath(os.path.join(self.task_config.ens_dir,
"cice_model.res."+str(mem)+".nc"))
ens_cice_list.append([fname_in, fname_out])
fname_out = os.path.realpath(os.path.join(self.task_config.ens_dir,
"cice_model.res.output."+str(mem)+".nc"))
Copy link
Copy Markdown
Collaborator Author

Choose a reason for hiding this comment

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

Note: this currently copies restarts twice (to keep input and output separate), we can change this before the merge to only copy the file once, and change it in place. (I coded input and output separately so I can easily see what the recentering does).

Comment on lines +53 to +67
variable change name: Soca2Cice
arctic:
seaice edge: 0.4
shuffle: true
rescale prior:
rescale: true
min hice: 0.5
min hsno: 0.1
antarctic:
seaice edge: 0.4
shuffle: true
rescale prior:
rescale: true
min hice: 0.5
min hsno: 0.1
Copy link
Copy Markdown
Collaborator Author

Choose a reason for hiding this comment

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

note: these parameters should be the same as in soca2cice used in the marineanlckpt; ideally we would reference the same yaml section in both places.

@shlyaeva
Copy link
Copy Markdown
Collaborator Author

shlyaeva commented Feb 6, 2025

After one cycle (once the deterministic background includes the large analysis increment from the previous cycle):
CICE restart for member 5 before recentering:
Screenshot 2025-02-06 at 2 42 43 PM
and after recentering:
Screenshot 2025-02-06 at 2 42 50 PM

@shlyaeva shlyaeva changed the base branch from develop to feature/lowres-b February 6, 2025 21:51
@guillaumevernieres guillaumevernieres marked this pull request as ready for review February 6, 2025 21:57
@guillaumevernieres guillaumevernieres merged commit be41879 into feature/lowres-b Feb 6, 2025
@guillaumevernieres guillaumevernieres deleted the feature/ice_recenter branch February 6, 2025 21:58
guillaumevernieres pushed a commit that referenced this pull request Apr 14, 2025
#1614)

# Description

In #1479 I introduced a
duplicate copy of ensemble members to marine recentering task by
mistake. This PR removes that.

I didn't test this as thoroughly as I normally would. Below is what I
ran to test (for both of the low-res tests I made sure that both rundir
and comdir were completely wiped out before I started the tests):
- all tests in C48mx500_hybAOWCDA up to and including ocnanalecen --
tests passed
- all tests in C48mx500_hybAOWCDA up to and including ocnanalecen except
for the letkf one (to mimic retro runs) -- tests passed
- introduced this change in my high-res 30-member ensemble experiment
that runs LETKF -- task completed successfully, 33 min instead of 38 min
for the previous cycle.

# Issues

Resolves #1613 

# Automated CI tests to run in Global Workflow
<!-- Which Global Workflow CI tests are required to adequately test this
PR? -->
- [ ] atm_jjob <!-- JEDI atm single cycle DA !-->
- [ ] C96C48_ufs_hybatmDA <!-- JEDI atm cycled DA !-->
- [ ] C96C48_hybatmaerosnowDA  <!-- JEDI aero/snow cycled DA !-->
- [ ] C48mx500_3DVarAOWCDA <!-- JEDI low-res marine 3DVar cycled DA !-->
- [x] C48mx500_hybAOWCDA <!-- JEDI marine hybrid envar cycled DA !-->
- [ ] C96C48_hybatmDA <!-- GSI atm cycled DA !-->

Co-authored-by: Anna Shlyaeva <anna.v.shlyaeva@noaa.gov>
jspsl pushed a commit to NOAA-PSL/GDASApp that referenced this pull request Apr 30, 2025
NOAA-EMC#1614)

# Description

In NOAA-EMC#1479 I introduced a
duplicate copy of ensemble members to marine recentering task by
mistake. This PR removes that.

I didn't test this as thoroughly as I normally would. Below is what I
ran to test (for both of the low-res tests I made sure that both rundir
and comdir were completely wiped out before I started the tests):
- all tests in C48mx500_hybAOWCDA up to and including ocnanalecen --
tests passed
- all tests in C48mx500_hybAOWCDA up to and including ocnanalecen except
for the letkf one (to mimic retro runs) -- tests passed
- introduced this change in my high-res 30-member ensemble experiment
that runs LETKF -- task completed successfully, 33 min instead of 38 min
for the previous cycle.

# Issues

Resolves NOAA-EMC#1613 

# Automated CI tests to run in Global Workflow
<!-- Which Global Workflow CI tests are required to adequately test this
PR? -->
- [ ] atm_jjob <!-- JEDI atm single cycle DA !-->
- [ ] C96C48_ufs_hybatmDA <!-- JEDI atm cycled DA !-->
- [ ] C96C48_hybatmaerosnowDA  <!-- JEDI aero/snow cycled DA !-->
- [ ] C48mx500_3DVarAOWCDA <!-- JEDI low-res marine 3DVar cycled DA !-->
- [x] C48mx500_hybAOWCDA <!-- JEDI marine hybrid envar cycled DA !-->
- [ ] C96C48_hybatmDA <!-- GSI atm cycled DA !-->

Co-authored-by: Anna Shlyaeva <anna.v.shlyaeva@noaa.gov>
DavidNew-NOAA pushed a commit that referenced this pull request Jan 16, 2026
#1614)

# Description

In #1479 I introduced a
duplicate copy of ensemble members to marine recentering task by
mistake. This PR removes that.

I didn't test this as thoroughly as I normally would. Below is what I
ran to test (for both of the low-res tests I made sure that both rundir
and comdir were completely wiped out before I started the tests):
- all tests in C48mx500_hybAOWCDA up to and including ocnanalecen --
tests passed
- all tests in C48mx500_hybAOWCDA up to and including ocnanalecen except
for the letkf one (to mimic retro runs) -- tests passed
- introduced this change in my high-res 30-member ensemble experiment
that runs LETKF -- task completed successfully, 33 min instead of 38 min
for the previous cycle.

# Issues

Resolves #1613 

# Automated CI tests to run in Global Workflow
<!-- Which Global Workflow CI tests are required to adequately test this
PR? -->
- [ ] atm_jjob <!-- JEDI atm single cycle DA !-->
- [ ] C96C48_ufs_hybatmDA <!-- JEDI atm cycled DA !-->
- [ ] C96C48_hybatmaerosnowDA  <!-- JEDI aero/snow cycled DA !-->
- [ ] C48mx500_3DVarAOWCDA <!-- JEDI low-res marine 3DVar cycled DA !-->
- [x] C48mx500_hybAOWCDA <!-- JEDI marine hybrid envar cycled DA !-->
- [ ] C96C48_hybatmDA <!-- GSI atm cycled DA !-->

Co-authored-by: Anna Shlyaeva <anna.v.shlyaeva@noaa.gov>
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.

Recenter sea ice ensembles

2 participants