Hello Giulia,
If I were you, I would read the recommendation appendix (E) and would try a few tips (maybe not suitable to your model, but let's try):
- IMPLICITATION FOR DEPTH = 0.55 (default), 0.6 or a value up to 0.9
- IMPLICITATION FOR VELOCITIES = 0.55 (default), 0.6 or a value up to 0.9
I would also change the vertical turbulence model rather than constant viscosity along the vertical:
e.g. VERTICAL TURBULENCE MODEL = 3 (k-epsilon). You can also try = 2 (mixing length) and try different mixing length models (perhaps with damping function). I do not know if they could suit your model, but you can also try = 5 (Spalart-Allmaras) or = 7 (k-omega). You can also try what is available in the GOTM library.
Hope this helps,
Chi-Tuan