Convergence of a numerical scheme for a mixed hyperbolic-parabolic system in two space dimensions