A reduced grid for a parallel global ocean general circulation model

A limitation of many global climate models with explicit finite-difference numerics is the timestep restriction caused by the decrease in cell size associated with the convergence of meridians near the poles. To keep the longitudinal width of model cells as uniform as possible, we apply a “reduced” grid to a three-dimensional primitive equation ocean-climate model. With this grid the number of cells in the longitudinal direction is reduced at high latitudes. The grid consists of subgrids which interact at interfaces along their northern and southern boundaries, where the resolution changes by a factor of three. We extend the finite-difference techniques to these interfaces, focusing on the conservation required to perform long time integrations, while preserving the staggered spatial arrangement of variables and the numerics used on subgrids. The common alternative used to reduce the timestep restriction caused by the spherical grid is the filtering of high-frequency modes from the high-latitude solution. The reduced grid allows an increased timestep while eliminating the need for filtering and reduces execution time per model step by roughly 20%. We implement the reduced grid model for parallel computer architectures with two-dimensional domain decomposition and message passing, with speedup results similar to those of the original model. We present results of model runs showing small effects on the solution and sizable improvements to the execution time.