Ocean Modeling Discussion

It ran stably overnight. The waves are gone and the ice model supresses submesoscale activity until the eddies are hardly visibly in the vorticity field. I can't attach an animation, so here are a few frames. I'm happy with this, it looks reasonable. I'm not sure whether it was my overhaul of bulk ...

Hi Tale,

I solved the bit for bit problem. This was hard to track down, but I turned the ice model off and still had the issue. Needed to provide the forcing for bulk fluxes in a .nc file instead of all analytical. I couldn't track down exactly what analytical settings were causing the problem, but ...

Thanks for posting this, I will test the new code. I was using the older version of your branch that had the BRY fixes in it. One recent bug I've discovered is that the ROMS ice model is not bit for bit and will seg fault or impose grid scale noise depending on the tiling. I've attached some of the ...

Agreed, I think something like this has to be added to mod_clima.F (likely just the start to actually get CLM support to work):
# ifdef ICE_MODEL
# if defined AICLM_NUDGING || defined MICLM_NUDGING
real(r8), pointer :: aiclm(:,:)
real(r8), pointer :: hiclm(:,:)
real(r8), pointer :: uiclm ...

Hi Tale and others,

I'm working with Brianna on the Beaufort setup. I figured out last night that the code to process BRY files was missing in get_data.F. Varinfo.yaml needs to be updated as well. My impression is that ice nudging and the ability to read in ice CLM files is also broken ...

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