Open xylar opened 2 weeks ago
This has been moved from https://github.com/E3SM-Ocean-Discussion/E3SM/pull/90, which includes relevant discussion.
@mark-petersen, you are the best person I could think of to sanity check changes to init mode. Let me know if you can do that.
@mark-petersen, this is a high priority for the FAnSSIE project. Can you review when you can?
Sorry, I forgot I already pinged you...
@sbrus89, @mark-petersen is away and I need someone to review this that knows about MPAS-Ocean's global ocean init mode. You're the best candidate. Would you be up for reviewing this? If not, remove yourself as a reviewer and I'll find someone else.
@xylar, sure I will take a look.
This merge updates MPAS-Ocean init mode for
global_ocean
in several ways that are critical for initialization that is compatible with planned coupling to MALI.The fields
landIcePressure
,landIceDraft
andssh
are now read from a topography file. This means thatlandIcePressure
is no longer computed based on the density at the ocean surface in init mode. Similarly,ssh
is not required to be the same aslandIceDraft
.Support is removed for reading land-ice topography fields from a dataset on a lat-lon grid. We no longer use this capability and maintaining it is not worth the effort.
The land-ice topography data set is now required to have
landIceFloatingFracObserved
as a field, and this is used to setlandIceFloatingFraction
. Previously,landIceFracObserved
was being used to compute bothlandIceFraction
andlandIceFloatingFraction
(assuming no grounded ice in the domain).