This paper deals with the development of efficient numerical solvers for the primitive equations of the ocean in turbulent regime. We derive the numerical approximation of a reduced model by the Smagorinsky turbulence model that includes stabilization of the pressure discretization by a penalty...