The combined Lagrangian advection method

We present and test a new hybrid numerical method for simulating layerwise-two-dimensional geophysical flows. The method radically extends the original Contour-Advective Semi-Lagrangian (CASL) algorithm [5] by combining three computational elements for the advection of general tracers (e.g. potential vorticity, water vapor, etc.): (1) a pseudo-spectral method for large scales, (2) Lagrangian contours for intermediate to small scales, and (3) Lagrangian particles for the representation of general forcing and dissipation. The pseudo-spectral method is both efficient and highly accurate at large scales, while contour advection is efficient and accurate at small scales, allowing one to simulate extremely fine-scale structure well below the basic grid scale used to represent the velocity field. The particles allow one to efficiently incorporate general forcing and dissipation. (C) 2010 Elsevier Inc. All rights reserved.