[MITgcm-support] NLFS=4 with z coordinates.. and sea ice

[Jean-Michel Campin]

Hi Dan,

It should work, not used much without z*, precisely for the reason you mentioned,
see also the Ocean modelling paper, doi:10.1016/j.ocemod.2008.05.005 :
"Sea ice-ocean coupling using a rescaled vertical coordinate z*"

Note that you could always bring back an inconsistent "levitating" sea ice with
a recently added param (in "data", first namelist): sIceLoadFac=0., 

Cheers,
Jean-Michel

On Tue, Sep 22, 2026 at 11:18:26AM +0000, Daniel Goldberg wrote:
> Hello!
> 
> I am curious if anyone has used nonlinear free surface (NLFS=4) with z-coordinates and with the sea ice model?
> 
> I suspect that it will be limited by the thickness of the top cell relative to the sea ice thickness, and am prepared to deal with surf_adjustment errors if they occur.
> 
> But aside from this, Im wondering if there are any numerics/packages that simply won't work, of if i will get a RT error simply for trying this combination of parameters?
> 
> Thanks!
> dan
> 
> --
> Daniel Goldberg, PhD (he/him)
> Reader in Glaciology
> School of Geosciences, University of Edinburgh
> Geography Building, Drummond Street, Edinburgh EH8 9XP
> https://dngoldberg.github.io/
> 
> Note: My institution's virus protection makes any attachments of code source files unreadable. Please zip, rename, or share via cloud.
> The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336. Is e buidheann carthannais a th' ann an Oilthigh Dhùn Èideann, clàraichte an Alba, àireamh clàraidh SC005336.

> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support