Sweep interpolation: a cost-effective semi-Lagrangian scheme in the Global Environmental Multiscale model

Abstract

The interpolation process is the most computationally expensive step of the semi-Lagrangian (SL) approach for solving advection and is commonly used in numerical weather prediction (NWP) models. It has a significant impact on the accuracy of the solution and can potentially be the most expensive part of model integration. The sweep algorithm, which was first described by Mortezazadeh and Wang (2017), performs SL interpolation with the same computational cost as a third-order polynomial scheme but at the accuracy of the fourth order. This improvement is achieved by using two third-order backward and forward polynomial interpolation schemes in two consecutive time steps. In this paper, we present a new application of the sweep algorithm within the context of global forecasts produced with Environment Climate Change Canada's Global Environmental Multiscale (GEM) model. Results show that the SL scheme with sweep interpolation is computationally more efficient compared to a conventional SL scheme with fourth-order polynomial interpolation, especially when a very large number of passive tracers are advected. An additional advantage of this new approach is that its implementation in a chemical and weather forecast model requires minimum modifications of the interpolation weighting coefficients. An analysis of the computational performance for a set of theoretical benchmarks as well as a global ozone forecast experiment show that up to 15 % reduction in total wall clock time is achieved. Forecasting experiments using the global version of the GEM model and the new interpolation show that the sweep interpolation can perform very well in predicting ozone distribution, especially in the tropopause region, where transport processes play a significant role.