In this paper, we construct a combined upwinding and mixed finite element method for the numerical solution of a two-dimensional mean field model of superconducting vortices. An advantage of our method is that it works for any unstructured regular triangulation. A simple convergence analysis is given without resorting to the discrete maximum principle. Numerical examples are also presented.