THE COST-EFFECTIVENESS OF SEMI-LAGRANGIAN ADVECTION

The purpose of this paper is to ascertain the cost-effectiveness of se mi-Lagrangian advection schemes for a wide variety of geophysical flow s at all scales. The approach used is first to determine the minimum c omputational overhead associated with these schemes and then to examin e temporal Variability in the Lagrangian and Eulerian frames by employ ing simple turbulent cascade phenomenologies. The goal is to evaluate whether the Lagrangian variability is sufficiently slower than that of the Eulerian frame to overcome the computational overhead. It is foun d that the most efficient semi-Lagrangian schemes require a factor of 5-10 times more floating point operations per grid point per time step than the classic second-order leapfrog scheme. In the enstrophy casca de of 2D or quasigeostrophic turbulence, evolution of flow quantities is considerably slower in the Lagrangian frame and semi-Lagrangian adv ection schemes can be very cost-effective. In an energy cascade such a s the Kolmogorov range of 3D turbulence or the inverse cascade of QG o r 2D turbulence, the Lagrangian evolution remains slower than the Eule rian evolution. However, the difference is very much less than in the enstrophy cascade. Since the computational overhead of semi-Lagrangian schemes is considerable, they are at best marginally cost-effective a t current resolutions for these flows, which prevail in the atmosphere at scales below 300-400 km. In the presence of stationary forcing fie lds in the Eulerian frame, the time step must respect the advective ti mescale even in the Lagrangian frame, at length scales where the forci ng is significant. Hen semi-Lagrangian schemes are not recommended.