Abstract:The demand for substantial increases in the spatial resolution of global weather and climate prediction models makes it necessary to use numerically efficient and highly scalable algorithms to solve the equations of large‐scale atmospheric fluid dynamics. For stability and efficiency reasons, several of the operational forecasting centres, in particular the Met Office and the European Centre for Medium‐Range Weather Forecasts (ECMWF) in the UK, use semi‐implicit semi‐Lagrangian time‐stepping in the dynamical core of the model. The additional burden with this approach is that a three‐dimensional elliptic partial differential equation (PDE) for the pressure correction has to be solved at every model time step and this often constitutes a significant proportion of the time spent in the dynamical core. In global models, this PDE must be solved in a thin spherical shell. To run within tight operational time‐scales, the solver has to be parallelized and there seems to be a (perceived) misconception that elliptic solvers do not scale to large processor counts and hence implicit time‐stepping cannot be used in very high‐resolution global models. After reviewing several methods for solving the elliptic PDE for the pressure correction and their application in atmospheric models, we demonstrate the performance and very good scalability of Krylov subspace solvers and multigrid algorithms for a representative model equation with more than 10 10 unknowns on 65 536 cores on the High‐End Computing Terascale Resource (HECToR), the UK's national supercomputer. For this, we tested and optimized solvers from two existing numerical libraries (the Distributed and Unified Numerics Environment (DUNE) and Parallel High Performance Preconditioners (hypre)) and implemented both a conjugate gradient solver and a geometric multigrid algorithm based on a tensor‐product approach, which exploits the strong vertical anisotropy of the discretized equation. We study both weak and strong scalability and compare the absolute solution times for all methods; in contrast to one‐level methods, the multigrid solver is robust with respect to parameter variations.

Semi-Lagrangian shallow water equations solver on the cubed-sphere grid as a prototype of new-generation global atmospheric model

Improved Local Time Step for 2D Shallow-Water Modeling Based on Unstructured Grids

Parallel Multilevel Methods for Implicit Solution of Shallow Water Equations with Nonsmooth Topography on the Cubed-Sphere

A parallel well-balanced finite volume method for shallow water equations with topography on the cubed-sphere

A global nonhydrostatic dynamical core on cubed sphere using multi-moment finite volume method: formulation and preliminary test

A Fully Implicit Domain Decomposition Algorithm for Shallow Water Equations on the Cubed-Sphere.

2d Well-Balanced Positivity Preserving Shallow Flow Model With Wetting And Drying On Triangular Grids

ExoCubed: A Riemann-Solver based Cubed-Sphere Dynamic Core for Planetary Atmospheres

A Well-Balanced Numerical Model for the Simulation of Long Waves over Complex Domains

A global non-hydrostatic dynamical core using the spectral element method on a cubed-sphere grid

ExoCubed: A Riemann-solver-based Cubed-sphere Dynamic Core for Planetary Atmospheres

Review on development of a scalable high-order nonhydrostatic multi-moment constrained finite volume dynamical core

A Global [math]-Staggered Composite Model for Shallow Water Equations with Latitude–Longitude Grid and Reductions in the Polar Regions

Scalability Studies of an Implicit Shallow Water Solver for the Rossby-Haurwitz Problem

Development of A Cell-Centered Godunov-Type Finite Volume Model for Shallow Water Flow Based on Unstructured Mesh

Massively parallel solvers for elliptic partial differential equations in numerical weather and climate prediction

The Spectral Element Method for the Shallow Water Equations on the Sphere

Newton-Krylov-Schwarz Method for a Spherical Shallow Water Model*

A mixed finite-element discretisation of the shallow-water equations

Development of a Hybrid Parallel MCV-based High-Order Global Shallow-Water Model

Analysis and improvement of a semi-Lagrangian exponential scheme for the shallow-water equations on the rotating sphere