High-latitude truncation errors of box-type primitive equation models

The 'box-type' finite-difference method includes a weighted average of the pressure gradient with weights proportional to the surface of the grid walls. It is shown that this averaging introduces first-order truncation errors near the poles. An example is shown in which the relative error is of zero order and the scheme produces large distortions in the solution at high latitudes.