A Bochev–Dohrmann–Gunzburger stabilization method for the primitive equations of the ocean

We introduce a low-order stabilized discretization of the primitive equations of the ocean with highly reduced computational complexity. We prove stability through a specific inf–sup condition, and weak convergence to a weak solution. We also perform some numerical tests for relevant flows.