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.

Improving Met Office Weather and Climate Forecasts with Bespoke Multigrid Solvers

Multigrid preconditioners for the mixed finite element dynamical core of the LFRic atmospheric model

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

Precision of the ENDGame: Mixed-precision arithmetic in the iterative solver of the Unified Model

Hybridised multigrid preconditioners for a compatible finite element dynamical core

Mixed-precision arithmetic in the ENDGame dynamical core of the Unified Model, a numerical weather prediction and climate model code

Mixed-Precision Computing in the GRIST Dynamical Core for Weather and Climate Modelling

A Scalable Helmholtz Solver in GRAPES over Large‐scale Multicore Cluster

A directive based hybrid met office NERC cloud model

Computationally-Efficient Climate Predictions using Multi-Fidelity Surrogate Modelling

Helmholtz Solving and Performance Optimization in Global/Regional Assimilation and Prediction System

Productive Performance Engineering for Weather and Climate Modeling with Python

Optimizing the Weather Research and Forecasting Model with OpenMP Offload and Codee

10M-core scalable fully-implicit solver for nonhydrostatic atmospheric dynamics.

A Barotropic Solver for High-Resolution Ocean General Circulation Models

Statistical Modelling of 2m Temperature and 10m Wind Speed Forecast Errors

Twelve Times Faster yet Accurate: A New State‐Of‐The‐Art in Radiation Schemes via Performance and Spectral Optimization

On Producing Reliable and Affordable Numerical Weather Forecasts on Public Cloud-Computing Infrastructure

Improving the Scalability of the Ocean Barotropic Solver in the Community Earth System Model

Advances in Land Surface Model-based Forecasting: A comparative study of LSTM, Gradient Boosting, and Feedforward Neural Network Models as prognostic state emulators