Abstract
A semi-Lagrangian advection scheme has been implemented for the linear passive advection equation. The advected scalar field is represented in terms of finite-element bases on a self-adaptive grid.
Results concerning the accuracy of the numerical scheme compare well to those with finite-difference schemes on regular meshes found in the literature. Quasi-monotone and quasi-conservative interpolation schemes have been implemented.
A parallelized version of the code shows good performance on a KSR-1 virtual shared memory computer.
