An Efficient Two‐Time‐Level Semi‐Lagrangian Semi‐Implicit Integration Scheme

The semi-implicit semi-Lagrangian integration technique enables numerical weather prediction models to be run with much longer timesteps than permitted by a semi-implicit Eulerian scheme. the choice of timestep can then be made on the basis of accuracy rather than stability requirements. to realize the full potential of the technique, it is important to maintain second-order accuracy in time; this has previously been achieved by applying it in the context of a three-time-level integration scheme. In this paper we present a two-time-level version of the technique which yields the same level of accuracy for half the computational effort. Unlike other efficient two-time-level schemes, ours does not rely on operator splitting. We apply this scheme to a variable-resolution barotropic finite-element regional model with a minimum gridlength of 100 km, using timesteps of up to three hours. the results are verified against a control run with uniformly high resolution, and are shown to be of similar accuracy to those of a semi-implicit Eulerian integration with a timestep of 10 minutes.

[1]  J. R. Bates,et al.  Improving the Estimate of the Departure Point Position in a Two-Time Level Semi-Lagrangian and Semi-Implicit Scheme , 1987 .

[2]  Herschel L. Mitchell,et al.  A Semi-Implicit Finite-Element Barotropic Model , 1977 .

[3]  A. Robert A Semi-Lagrangian and Semi-Implicit Numerical Integration Scheme for the Primitive Meteorological Equations , 1982 .

[4]  André Robert,et al.  A stable numerical integration scheme for the primitive meteorological equations , 1981 .

[5]  Janusz A. Pudykiewicz,et al.  Preliminary results From a partial LRTAP model based on an existing meteorological forecast model , 1985 .

[6]  Janusz A. Pudykiewicz,et al.  Some properties and comparative performance of the semi‐Lagrangian method of Robert in the solution of the advection‐diffusion equation , 1984 .

[7]  A. Staniforth,et al.  A Variable-Resolution Finite-Element Technique for Regional Forecasting with the Primitive Equations , 1978 .

[8]  A. Mcdonald Accuracy of Multiply-Upstream, Semi-Lagrangian Advective Schemes , 1984 .

[9]  J. Bates,et al.  Multiply-Upstream, Semi-Lagrangian Advective Schemes: Analysis and Application to a Multi-Level Primitive Equation Model , 1982 .

[10]  J. R. Bates An Efficient Semi-Lagrangian and Alternating Direction Implicit Method for Integrating the Shallow Water Equations , 1984 .

[11]  Harold Ritchie,et al.  Eliminating the Interpolation Associated with the Semi-Lagrangian Scheme , 1986 .

[12]  Monique Tanguay,et al.  Elimination of the Helmholtz Equation Associated with the Semi-Implicit Scheme in a Grid Point Model of the Shallow Water Equations , 1986 .

[13]  Dick Dee,et al.  A Fully Implicit Scheme for the Barotropic Primitive Equations , 1985 .

[14]  A. McDonald A Semi-Lagrangian and Semi-Implicit Two Time-Level Integration Scheme , 1986 .

[15]  Clive Temperton,et al.  Semi-Implicit Semi-Lagrangian Integration Schemes for a Barotropic Finite-Element Regional Model , 1986 .

[16]  S. Cohn,et al.  A Factored Implicit Scheme for Numerical Weather Prediction with Small Factorization Error , 1985 .

[17]  André Robert,et al.  Accuracy and Stability Analysis of a Fully Implicit Scheme for the Shallow Water Equations , 1986 .

[18]  D. Purnell Solution of the Advective Equation by Upstream Interpolation with a Cubic Spline , 1976 .

[19]  A. Robert,et al.  A semi-Lagrangian and semi-implicit numerical integration scheme for multilevel atmospheric models , 1985 .

[20]  A. Staniforth,et al.  Reply to comments on and addenda to “some properties and comparative performance of the semi‐lagrangian method of Robert in the solution of the advection‐diffusion equation” , 1985 .