[MITgcm-support] Vertical diffusion for non-hydrostatic ocean model with a stretched grid

[Burns, Paul]

Hi,

I was hoping that someone could give me some advice about the best way to define the vertical viscosity (and diffusivity) in MITgcm for a non-hydrostatic model using a stretched vertical grid.  I've used a stretched vertical grid to resolve very small features of interest (~0.2 m) in the upper ocean and to keep the simulation practical.

The viscosity should vary with the grid resolution (e.g., Smagorinsky, 1993), however, the MITgcm documentation says that the Smagorinsky scheme isn't implemented in the vertical. So would a reasonable approach be to compute a vertical profile of viscosities using a vertical Reynolds number and use that to set viscArNr?


Many thanks,
Paul




Dr Paul Burns MInstP FRMetS

Research Fellow

Geophysical and Astrophysical Fluid Dynamics Group

Department of Mathematics

Harrison Building

University of Exeter

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20230201/66386c39/attachment.html>