Numerically Converged Solutions of the Global Primitive Equations for Testing the Dynamical Core of Atmospheric GCMs

Solutions of the dry, adiabatic, primitive equations are computed, for the first time, to numerical convergence. These solutions consist of the short-time evolution of a slightly perturbed, baroclinically unstable, midlatitude jet, initially similar to the archetypal LC1 case of Thorncroft et al. The solutions are computed with two distinct numerical schemes to demonstrate that they are not dependent on the method used to obtain them.

These solutions are used to propose a new test case for dynamical cores of atmospheric general circulation models. Instantaneous horizontal and vertical cross sections of vorticity and vertical velocity after 12 days, together with tables of key diagnostic quantities derived from the new solutions, are offered as reproducible benchmarks. Unlike the Held and Suarez benchmark, the partial differential equations and the initial conditions are here completely specified, and the new test case requires only 12 days of integration, involves no spatial or temporal averaging, and does not call for physical parameterizations to be added to the dynamical core itself.