Ocean Modeling Discussion

I'm working with a student who has a regional model of the Beaufort Shelf. We plan on doing some advection scheme experiments to look at numerical mixing and I'm curious whether HSIMT can support negative temperatures. Maybe I'm misunderstanding, but there appears to be contradictory info about this ...

Thanks for sending this. I was unaware of the bug - but I am using COAWST v3.8, which has ROMS 4.1. I will need to see how difficult it is to git cherry pick the commit and deal with merge conflicts.

That said, I don't think GLS is overmixing in this case (recall I paired it with harmonic lateral ...

Just for an update, I found significant improvement with Richardson number smoothing paired with biharmonic mixing. See attached for a vorticity comparison with the same GLS simulation (that uses harmonic).

kpp_w_biharmonic_RI_smooth.png

RI_HORAVG Smooth Richardson number horizontally
RI ...

Hi All,

I'm doing a model-model comparison of the Texas-Louisiana shelf with ROMS and MPAS-O. Our ROMS model is validated (Hetland 2017 JPO, Kobashi & Hetland 2020 OD) and normally uses k-omega or k-epsilon for vertical mixing with great success. The mean resolution is ~1.5 km. 30 vertical layers ...

Hi All,

I'm doing some experiments with Smagorinsky viscosity/diffusion applied biharmonically along geopotential surfaces to explore its impacts on numerical mixing. In hmixing.F, I see that Bsmag is controlled by PecletCoef, which is equal to 1/12. Where does the 1/12 come from? It's unclear to ...

A little late, but I believe you're correct. Parker MacCready posted a thread about this about a decade ago, it gets at this issue as well. Took me too long to figure out myself until I discovered that thread, posting this to help aggregate.


I figured out the problem with my salt_rate test. salt ...

Is it possible to calculate the vertical tracer diffusion coefficients "AKt" or "AKs" offline with model output? I'm using the GLS k-epsilon vertical mixing scheme for context and would appreciate if there's a way to compute them using the model equations.

I'd like to construct an offline ...