Ocean Modeling Discussion

Not 100% sure, but I think there is a coding error in the application of gls_Kmin and gls_Pmin to limit values produced by horizontal and vertical advection of tke and gls for models using any of the GLS-wrapped versions of second moment closure (and hence entering ROMS/Nonlinear/gls_prestep.F &...

I had some exchange with Hans Burchard on this with CC to John Warner -- ought to have cc'd you too, Hernan. I summarize that exchange here: Burchard: Around 2005, we had a major update in the turbulence routines coordinated by Lars, and I cannot find any pre-2005 GOTM code where we use the limiting...

Here is a Matlab program (with apologies to open-source folk) to demonstrate the problem with the published & ROMS-coded equation for the smoothed Gh limiter, based on 3 different closure options that invoke it in ROMS within gls_corstep.F . High resolution of Gh is needed to see the problem in ...

Hi. I've been working on putting different turbulence closures into the COAWST version of ROMS and into GOTM, and I think I have found a longstanding problem in the published account of an aspect of the GLS turbulence closure model, and in the resulting implementation in ROMS. The problem involves a...

I've been working on putting in a MY2.5-variant parameterization of Langmuir turbulence (Harcourt,2013; 2015; + newer mods) that entails a component of the vertical momentum flux directed down the gradient of the surface waves' Stokes drift, i.e.: <u'w'> = - K_M du/dz - K_W du^s/dz where u is horizo...