Hi, all -
I've got a simple setup - periodic channel, flat bottom, sloping front, geostrophically-balanced along-channel jet. No forcing, no mixing. In fully prognostic setup the model maintains this state, as it should.
When I turn TS_FIXED on, it breaks down to something that looks like baroclinic instability in a matter of days. The final state that it seems to try to achieve, consists of two counter-flowing almost barotropic jets, which are obviously not in geostrophic balance.
The problem doesn't seem to be sensitive to temporal or spatial resolution.
Any idea of what may be going on?
Thanks!
- Andrey
Model details below:
Domain size: 20km (across channel), 2km(along channel), 150m (depth)
horizontal spacing: .4km, vertical spacing: 1-2m
AKT_BAK = AKV_BAK = 1e-5
VISC2 = RDRG = 0
CPP options:
# define TS_FIXED /*No evoluteion of tracers!*/
# define UV_ADV
# define UV_COR
# define UV_LDRAG
# define UV_VIS2
# define TS_DIF2
#define UV_U3HADVECTION
#define TS_U3HADVECTION
#define TS_C4VADVECTION
#define MIX_S_UV
#define MIX_S_TS
# define SOLVE3D
# define SALINITY
# define NONLIN_EOS
# define SPLINES
# define DJ_GRADPS
# define ANA_VMIX
# define ANA_BSFLUX
# define ANA_BTFLUX
# define ANA_M2CLIMA
# define ANA_M3CLIMA
# define ANA_TCLIMA
# define EASTERN_WALL
# define WESTERN_WALL
# define NS_PERIODIC