Abstract
A technique for following the circulation of a tracer in a turbulent fluid is developed from the integral form of the mass continuity equation. Numerical methods based on this technique are shown to be stable, to ensure that the total tracer mass is conserved, and that the mass in any region is always nonnegative. As an illustration of the utility of the technique, a numerical method is developed for a two-dimensional model of the stratosphere.