I'm having trouble with what looks like a numerical instability that starts in the vertical velocity field.
My setup consists of an idealized flat bottom shelf + slope at the southern end of a beta-plane channel with a west to east (along-isobath) barotropic inflow (2cm/s). There is an anticyclonic eddy in the channel that moves southwards towards the topography since it is on the beta plane. The southern boundary (i.e. the coast) is closed and the remaining three are open. I have sponge layers along all open boundaries. The boundary conditions are
- * Shc or RadNud for 2D momentum
* Che for free surface
* RadNud for 3D momentum
* RadNud on tracers at the western (inflow boundary) and Rad everywhere else
* VolCons is not enabled
* RADIATION_2D is enabled
The instability happens at the edge of the continental slope (where the vertical levels are stretched more) but not over the flat bottom region. For this case, bottom slope = 0.02, N = 3e-3, f = 5e-5 i.e., a slope Burger number of 1.2. The trouble is
- * results in vertical velocity growing with time to O(10^-4) m/s. I haven't run it long enough to see if it results in blow up eventually.
* Similar patterns are seen in u,v,rho also. magnitude of noise is O(1cm/s) in u.
* Independent of where the eddy is in the domain. It is seen at locations pretty far away from the eddy.
* Has a wavelength of 10-15km in the horizontal = 9-10 grid points. So, my horizontal viscosity is not effective at killing it.
* Has a time period of 20 inertial periods
* Looks like grid-scale noise in the vertical (see figure).
* It is always seen where the vertical grid spacing is maximum. Changing the stretching parameters only makes it appear at different depth levels in the vertical i.e., wherever the dz is large.
* At a horizontal resolution dx=dy=1.5, Doubling the number of vertical levels to 60 delays the onset of instability by approx. 1.5 times but it still appears at a time where I'm interested in the dynamics.
* Increasing horizontal resolution to dx=dy=1km doesn't change it by much.
Have any of you seen anything like this before? Any suggestions on what might be causing it or what might help mitigate it?