Abstract
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.