We present an explicit sixth order compact finite difference scheme for fast high accuracy numerical solutions of the two dimensional convection diffusion equation with variable coefficients. The sixth order scheme is based on the well-known fourth order compact scheme, the Richardson extrapolation technique, and an operator interpolation scheme. For a particular implementation, we use multisca...