Abstract
A second order numerical method for the primitive equations
(PEs) of large-scale oceanic flow formulated in mean vorticity is proposed and
analyzed, and the full convergence in $L^2$ is established. In the reformulation of
the PEs, the prognostic equation for the horizontal velocity is replaced by evolutionary equations for the mean vorticity field and the vertical derivative of
the horizontal velocity. The total velocity field (both horizontal and vertical)
is statically determined by differential equations at each fixed horizontal point.
The standard centered difference approximation is applied to the prognostic
equations and the determination of numerical values for the total velocity field
is implemented by FFT-based solvers. Stability of such solvers are established
and the convergence analysis for the whole scheme is provided in detail.