A fast converging and concise algorithm for computing the departure points in semi‐Lagrangian weather and climate models

Abstract

AbstractAccurate calculation of departure points is crucial for the overall accuracy of a semi‐Lagrangian advection scheme. The iterative method used in the semi‐Lagrangian advection scheme of the ECMWF model IFS (Integrated Forecast System) is based on SETTLS (Stable Extrapolation Two‐Time‐Level Scheme). The scheme converges slowly due to the long timesteps that the spectral semi‐implicit, semi‐Lagrangian formulation permits: five iterations are needed for accurate determination of the departure points. Spherical coordinates increase the expense and complexity of the relevant code while computations become sensitive to arithmetic precision and require further approximations in the polar high‐curvature region. In this article, an accurate alternative, derived by reformulating the current SETTLS scheme on a geocentric Cartesian framework, is proposed. The new scheme simplifies considerably the computations by avoiding the update of the Lagrangian trajectory equations in spherical coordinates, which involve additional expensive metric terms, eliminates approximations, and is advantageous for single precision. Methods to accelerate the convergence of the departure points were also investigated. The most accurate and efficient approach was to start the SETTLS iterations from the departure point of the previous timestep. Testing the proposed method in high‐resolution global forecast and 4D‐Var analysis experiments yields equally accurate results as the existing method, while reducing the cost of the iterative scheme.