Ocean Modeling Discussion

It makes sense to me now. Thank you John for this very clear explanation.

Thank John and Greg for helping with this. Yes I think "gls" = "epsilon" for my case, so I can directly use the "gls" variable from the output file as dissipation. And I also compared the vertical distribution of gls (tke dissipation) and it is generally consistent with shear production. Below is ...

Hi John, Thanks a lot for this reply. I am well aware of this paper and the GLS method. I am just a bit confused by what ROMS writes in the output file. I think gls = epsilon for GLS k-epsilon closure, but according to the glossary in ocean.in, the output gls is always tke*l even if I am not using ...

Hi all, I wondered what the "gls" is in the output file. It is "tke*l" ("l" is the generic length) according to the glossary I found in the .in file. So then I could calculate TKE dissipation using "epsilon = c_mu0^3*tke^1.5/l = c_mu0^3*tke^1.5/(gls/tke)" - I am using the GLS k-epsilon closure. But ...

Thank you so much! I'll take a look at that section.

Hi all, I wondered how the bottom boundary layer height is calculated in ROMS. I noticed there is an Hbbl variable that can be directly written in the history files, but does anyone know how that was calculated in the code, and is there any way that I can calculate it using other model output ...