Most numerical weather prediction models rely on a terrain-following coordinate framework. The computational mesh is thus characterized by inhomogeneities with scales determined by the underlying topography. Such inhomogeneities may affect the truncation error of numerical schemes. In this study, a new class of terrainfollowing coordinate systems for use in atmospheric prediction models is proposed. Unlike conventional systems, the new smooth level vertical (SLEVE) coordinate yields smooth coordinates at mid-and upper levels. The basic concept of the new coordinate is to employ a scale-dependent vertical decay of underlying terrain features. The decay rate is selected such that small-scale topographic variations decay much faster with height than their large-scale counterparts. This generalization implies a nonlocal coordinate transformation. The new coordinate is tested and compared against standard sigma and hybrid coordinate systems using an idealized advection test. It is demonstrated that the presence of coordinate transformations induces substantial truncation errors. These are critical for grid inhomogeneities with wavelengths smaller than approximately eight grid increments, and may overpower the regular-grid truncation error of the underlying finite-difference approximation. These results are confirmed by a theoretical analysis of the truncation error. In addition, the new coordinate is tested in idealized and real-case numerical experiments using a nonhydrostatic model. The simulations using the new coordinate yield a substantial reduction of small-scale noise in dynamical and thermodynamical model fields.

No abstract

The Global Environmental Multiscale (GEM) model is the Canadian atmospheric model used for meteorological forecasting at all scales. A limited-area version now also exists. It is a gridpoint model with an implicit semi-Lagrangian iterative space-time integration scheme. In the ''horizontal,'' the equations are written in spherical coordinates with the traditional shallow atmosphere approximations and are discretized on an Arakawa C grid. In the ''vertical,'' the equations were originally defined using a hydrostatic-pressure coordinate and discretized on a regular (unstaggered) grid, a configuration found to be particularly susceptible to noise. Among the possible alternatives, the Charney-Phillips grid, with its unique characteristics, and, as the vertical coordinate, log-hydrostatic pressure are adopted. In this paper, an attempt is made to justify these two choices on theoretical grounds. The resulting equations and their vertical discretization are described and the solution method of what is forming the new dynamical core of GEM is presented, focusing on these two aspects.

No abstract

EudraCT Number: 2012-000232-25; clinicaltrials.gov Identifier: NCT02184819.