Numerical approximation of the three‐dimensional ocean primitive equations