New Mesh: SOwISC12to30E3r3

[xylar]

Long name: SOwISC12to30kmL80E3SMv3r3

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

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

This mesh differs from r2 in #753 in that:

• holes in land-ice mask are filled and isolated "islands" of land ice are removed (#800, #844, #845, #846, #848)
• the landIceFraction is capped at 1 (#805)
• it has 80 (instead of 64) vertical levels with a highest vertical resolution of 2 m (instead of 10)
• it specifies a minimum of 5 layers with a minimum thickness of 2 m (usring #812)
• the sea surface height is iteratively adjusted for hydrostatic balance (following #813)
• Omega and pi have been corrected in https://github.com/E3SM-Project/E3SM/pull/6481, brought in through #839
• Bug fix to the haney-number-constrained vertical coordinate code in https://github.com/E3SM-Project/E3SM/pull/6535, updated in the E3SM-Project submodule in https://github.com/MPAS-Dev/compass/pull/849
• Bug fix related to masking of ice thickness: #854
• Bug fix related to constraining ice thickness to flotation thickness: #855

Mesh, initial condition, dynamic adjustment and files for E3SM are on Chrysalis at:

/lcrc/group/e3sm/ac.xylar/compass_1.4/chrysalis/e3smv3-meshes/sowisc12to30e3r3

Checklist

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

[xylar]

We will switch to 80 levels

[xylar]

Probably sticking with 12 km highest res, but up for further discussion.

[xylar]

@cbegeman and I discussed switching to having landIcePressure be supplied from observations and ssh be adjusted to be dynamically consistent. (At the moment, Compass does the opposite.) This would be more consistent with plans for MALI coupling in the near future.

[xylar]

My preliminary testing shows that the new ssh adjustment technique from #813 works just fine and converges sufficiently in 4 iterations (change of about 10 cm). So I think this is ready to include in the next version of the mesh.

[xylar]

Here are the figures from the initial_state step, showing the water-column thickness > 10 m in both the open ocean and in cavities:

[xylar]

Needs https://github.com/E3SM-Project/E3SM/pull/6535