methods in the numerical treatment of advection-diffusion problems, while also conserving
mass and handling general boundary conditions. This paper discusses some ideas that
could extend earlier formulations to variable coefficients and multiple space dimensions.