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

Add new SOwISC12to30E3r3 ocean and sea-ice mesh #6631

Merged
merged 11 commits into from
Oct 9, 2024

Conversation

xylar
Copy link
Contributor

@xylar xylar commented Sep 20, 2024

Long name: SOwISC12to30kmL80E3SMv3r3

This version of the Southern Ocean Regionally Refined Mesh (SORRM) has resolution that is:

  • 12 km around Antarctica
  • 30 km elsewhere
    It is intended to be similar to the IcoswISC30E3r5 mesh except in the Southern Ocean and around Antarctica.

This is a proposed E3SM v3 (E3) mesh for polar (formerly cryosphere) simulations. This is revision 3 (r3) of the mesh, which includes many changes from revision 2 outlined in MPAS-Dev/compass#807. The mesh is tagged on Compass as: https://github.com/MPAS-Dev/compass/releases/tag/mesh_SOwISC12to30E3r3 to aid reproduction in the future.

The mesh and the G- and B-case results are being reviewed here:
https://acme-climate.atlassian.net/wiki/spaces/OO/pages/4441145345/Review+SOwISC12to30E3r3

A 30-year G- and 50-year B-case have been run, and analysis is posted on the review page. Here are the simulation pages:
https://acme-climate.atlassian.net/wiki/spaces/PSC/pages/4575821825/20240829.GMPAS-JRA1p5-DIB-PISMF.TL319_SOwISC12to30E3r3.mesh-review-test.chrysalis
https://acme-climate.atlassian.net/wiki/spaces/PSC/pages/4575068161/20240904.v3.SORRME3r3.piControl.chrysalis.mesh-review

Relevant discussion can be found at E3SM-Ocean-Discussion#106

Copy link
Member

@darincomeau darincomeau left a comment

Choose a reason for hiding this comment

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

Approved based on test simulations referenced above and discussion in E3SM-Ocean-Discussion#106

@darincomeau
Copy link
Member

I added a default finidat for 1850 conditions with the r05 land grid, which was interpolated from the end of the v3 LR spinup.
10 year test using this finidat is here: https://acme-climate.atlassian.net/wiki/spaces/PSC/pages/4628250625/20240923.v3.SORRME3r3.piControl.lndIC.chrysalis
Notably there is no initial drastic drop in sea level as in the case with a blank finidat.

@xylar
Copy link
Contributor Author

xylar commented Sep 28, 2024

Rebased to fix conflicts with #6643

@xylar
Copy link
Contributor Author

xylar commented Oct 1, 2024

@alicebarthel, @cbegeman and @proteanplanet, could you please review and hopefully approve this week?

@cbegeman
Copy link
Contributor

cbegeman commented Oct 1, 2024

I added a default finidat for 1850 conditions with the r05 land grid, which was interpolated from the end of the v3 LR spinup.
10 year test using this finidat is here: https://acme-climate.atlassian.net/wiki/spaces/PSC/pages/4628250625/20240923.v3.SORRME3r3.piControl.lndIC.chrysalis
Notably there is no initial drastic drop in sea level as in the case with a blank finidat.

I examined the land ice mass balance fields and they look as expected. We hope that the negative river runoff anomaly around Amery shown here will be reduced once we have a configuration with the deep firn model.

Copy link
Contributor

@cbegeman cbegeman left a comment

Choose a reason for hiding this comment

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

Approved on the basis of water column depth around ISC, examination of @proteanplanet's mesh figures, and review of simulation results in MPAS-Analysis pages linked above and separately generated surface mass flux figures in the SO.

Copy link
Contributor

@alicebarthel alicebarthel left a comment

Choose a reason for hiding this comment

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

Approving after examining @proteanplanet's mesh figures, @darincomeau's and @cbegeman's MPAS-Analysis pages.
The Ocean Heat uptake near Zapiola is visible immediately, but I could not find any mesh or topography irregularities in the region, so I don't think this mesh specifically enhances the dynamical problem.

Copy link
Contributor

@proteanplanet proteanplanet left a comment

Choose a reason for hiding this comment

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

I have thoroughly reviewed this mesh and compared it with the IcoswISC30E3r5 mesh, including checking for integrity of the ice shelf component. Differences include a minimum ocean depth of 12 m as compared to 20.1m in the IcoswISC30E3r5 mesh to accommodate greater refinement in ice shelf cavities. Figures are being updated at https://acme-climate.atlassian.net/wiki/spaces/OO/pages/edit-v2/4441145345 for later use in publications.

@proteanplanet
Copy link
Contributor

An addendum to my acceptance of this PR. The final set of plots are now posted on the mesh review page at https://acme-climate.atlassian.net/wiki/spaces/OO/pages/4441145345/Review+SOwISC12to30E3r3. Something I want to make sure is on everyone's minds is the nearly halving of the minimum depth in the Arctic for which we need to keep a close eye on how the sea ice evolves in SORRM runs:

MPAS_OceanIce_SOwISC12to30E3r3_IcoswISC30E3r5_zoom_bathymetry_56

The issue we need to consider in association with this mesh is the spreading functions applied for liquid and frozen runoff, which particularly affect the Arctic and are bathymetrically sensitive. I suggest we consider a similar set of functions to those applied to the E3SM V3 HR mesh, with which @cbegeman is familiar. However, that's not a consideration in whether or not this mesh is ready for implementation. It's now a matter of how we use this excellent mesh. Thank you, @xylar, @darincomeau and @cbegeman for your hard work putting this together.

@darincomeau
Copy link
Member

I pushed an adjustment to the Chrysalis 'M' PE layout to this: https://pace.ornl.gov/exp-details/199117

@jonbob jonbob added the BFB PR leaves answers BFB label Oct 8, 2024
jonbob added a commit that referenced this pull request Oct 8, 2024
Add new SOwISC12to30E3r3 ocean and sea-ice mesh

Long name: SOwISC12to30kmL80E3SMv3r3
This version of the Southern Ocean Regionally Refined Mesh (SORRM) has
resolution that is:
* 12 km around Antarctica
* 30 km elsewhere
It is intended to be similar to the IcoswISC30E3r5 mesh except in the
Southern Ocean and around Antarctica. This is a proposed E3SM v3 (E3)
mesh for polar (formerly cryosphere) simulations.

[BFB]
@jonbob
Copy link
Contributor

jonbob commented Oct 8, 2024

Passes:

  • SMS_P960.T62_SOwISC12to30E3r3.GMPAS-IAF.chrysalis_intel
  • ERP_Ld3.ne30pg2_r05_IcoswISC30E3r5.WCYCL1850.chrysalis_intel.allactive-pioroot1

and copied all new files to the inputdata repo. Merged to next

@xylar
Copy link
Contributor Author

xylar commented Oct 9, 2024

I also synced over (and chmodded) the files that belong in inputdata/shared and diagnostics.

@jonbob jonbob merged commit 9b58110 into E3SM-Project:master Oct 9, 2024
9 checks passed
@jonbob
Copy link
Contributor

jonbob commented Oct 9, 2024

merged to master

@xylar xylar deleted the mesh/add-sowisc12to30e3r3 branch October 10, 2024 00:27
@xylar
Copy link
Contributor Author

xylar commented Oct 10, 2024

Woot, woot! Thanks @jonbob!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
BFB PR leaves answers BFB mpas-ocean mpas-seaice RRM Regionally refined model
Projects
None yet
Development

Successfully merging this pull request may close these issues.

6 participants