Abstract
The planetary geostrophic equations with inviscid balance equation are reformulated in an alternate form, and a fourth-order finite difference numerical method of solution is proposed and analyzed in this article. In the reformulation, there is only one prognostic equation for the temperature field and the velocity field is statically determined by the planetary geostrophic balance combined with the incompressibility condition. The key observation is that all the velocity profiles can be explicitly determined by the temperature gradient, by utilizing the special form of the Coriolis parameter. This brings convenience and efficiency in the numerical study. In the fourth-order scheme, the temperature is dynamically updated at the regular numerical grid by long-stencil approximation, along with a one-sided extrapolation near the boundary. The velocity variables are recovered by special solvers on the 3-D staggered grid. Furthermore, it is shown that the numerical velocity field is divergence-free at the discrete level in a suitable sense. Fourth order convergence is proven under mild regularity requirements.