Ocean Modeling Discussion

Sorry, but those two prestep routines have errors in the new code changes still. In each of the prestep routines there is a pair of lines, one for each horizontal direction of third-order upstream advection, that says: IF (cff.gt.0.0_r8) THEN The first of these should read IF (XF(i,j).gt.0.0_r8 ...

There is a problem with the advection schemes for turbulence if you are not specifying one of the 2 centered schemes, which may be contributing to the noise you are dealing with. Please see my modified code posted for GLS and MY25 here: https://www.myroms.org/forum/viewtopic.php?p=25663#p25663 I ...

Attached are revisions of 4 code modules for the predictor and corrector steps of GLS and MY25 second moment closure options. They correct an incongruity in the default cpp selection of predictor and corrector advection schemes, one that seems to go back at least 2 decades to version 2.0 per ...

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 ...