Ocean Modeling Discussion

Hi Wilkin,

Thanks. I'll try it and will update the result in this thread.

Hi Wilkin,

Do you mean GLS might not be a good choice?

Also should I try to switch back the advection scheme to U3 and C4 and see if it helps?

Hi Wilkin, I just check the roms.in and I have the advection scheme as this: ! The user has the option of specifying the full Keyword or the first ! two letters, regardless if using uppercase or lowercase. If nested ! grids, specify values for each grid (see glossary below). Hadvection == A4 ...

Hi Wilkin, I believe by default the horizontal advection is 3rd-order upstream bias and vertical advection is 4th-order centered. I didn't specify advection schemes myself and so I think it uses the default setting. I set the horizontal mixing coefficients as 1 m^2/s and I attach the code below ...

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

Ah! ok , yes i gave you the k-e. And we dont see those wiggles with k-e. I dont typically use the k-gen settings. Wonder if you can adjust any parameters to get rid of that noise. try this: Generic length-scale turbulence closure parameters. GLS_P == 2.0d0 ! GLS_M == 1.0d0 GLS_N == -0.67d0 GLS_Kmin ...

Hi John, These are the one used to represent the k-epsilon turbulent closure model, not the gls one right? I've tried the k-epsilon one before and there is no noise for the vertical mixing. But k-epsilon model tends to overmix and not very appropriate for stratified estuary situation so that's why ...

Hi John, I show the gls setting below here. Please check it. ! Generic length-scale turbulence closure parameters. GLS_P == 2.0d0 ! GLS_M == 1.0d0 GLS_N == -0.67d0 GLS_Kmin == 1.0d-8 GLS_Pmin == 1.0d-8 GLS_CMU0 == 0.5544d0 GLS_C1 == 1.0d0 GLS_C2 == 1.22d0 GLS_C3M == 0.1d0 GLS_C3P == 1.0d0 GLS_SIGK ...

Hi John, Thanks for your reply! It took me a while to play around with all the options and found out the best for our study. For the white box, it's just due to the vertical grid is evenly layered so the grid size is bigger in the middle compared to the side of the channel because of depth change ...

Hi Everyone, Thanks for reading this post. I'm doing my research on circulation and stratification in microtidal estuaries. I use an idealized straight-channel estuary grid (show below) and I'm looking at the wind effects on the circulation and stratification. here is my model domain. The grid ...

If you have a linear salt profile then you should have constant d(rho)/dz and then the vertical diffusion term (d/dz*AKT*constant) would be identically zero. So, I doubt that explicit vertical diffusion is the the issue. Besides, AKT_BAK is just a background minimum to the chosen vertical ...

Hello everyone, Thanks for reading this post. I'm doing my research on circulation and stratification in microtidal estuaries. I use an idealized straight-channel estuary grid (show below) and I want to check the numerical circulation isn't significant in my grid situation. here is my model domain ...