open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

Attempting to setup a baroclinic model which is forced UV,H, Salinity and Temp on the boundaries. The shore parallel boundary has prescribed tracers, velocities and water levels while the two shore normal boundaries have only prescribed velocities and tracers. Binary boundary files are provided with W=0 and some custom modifications have been made to sum the contribution of TPXO tides.

The model has been running relatively stable with a limited number of layers and we've been increasing the number of horizontal layers to better reproduce the stratification. I am however, noticing a strong vertical velocity across the full depth developing at a single node on one of the lateral boundaries. I cannot attribute to this anything obvious - bathymetry is relatively smooth within this location etc.

My understanding is that the vertical velocity remains free regardless of if the U,V components are prescribed but would like to constrain vertical velocities on the boundaries to be 0 in the hope to suppress this issue from developing. Setting LIWBOL%I(IBORD) = KENTU (6) with a prescribed value of WBORL%R(IBORD) = 0.D0 within BORD3d does not seem to have any effect though.

Wondering if anybody has any thoughts on the issue?

Thanks,
Toby

Managed to resolve the issue by setting LIWBOL%I(IBORD) = KENTU (i.e. prescribed) within the LIMI3D subroutine which has resolved the vertical instabilities on some of my boundary nodes.

Given this is on ocean boundaries where I expect the vertical velocities to be rather small it seems likes a reasonable assumption to make. Wondering if anybody has any thoughts?