In this work, we consider a boundary value problem involving Caputo derivatives defined in the plane. We develop a fast locally refined finite volume method for
variable-coefficient conservative space-fractional diffusion equations in the plane to resolve boundary layers of the solutions. Numerical results are presented to show the
utility of the method.