The model utilizes
the modified Mellor and Yamada (1982) level 2.5
(MY-2.5) and Smagorinsky turbulent closure schemes
as default setups for vertical and horizontal mixing,
respectively, with an option for various turbulent parameterizations
using the Generalized Ocean Turbulence
Model (GOTM) (Burchard, 2002).