  This article concerns a scalar multidimensional conservation law where the
flux is of Panov type and may contain spatial discontinuities. We define a
notion of entropy solution and prove that entropy solutions are unique. We
propose a Godunov-type finite volume scheme and prove that the Godunov
approximations converge to an entropy solution, thus establishing existence of
entropy solutions. We also show that our numerical scheme converges at an
optimal rate of $\mathcal{O}(\sqrt{\D t}).$ To the best of our knowledge,
convergence of the Godunov type methods in multi-dimension and error estimates
of the numerical scheme in one as well as in several dimensions are the first
of it's kind for conservation laws with discontinuous flux. We present
numerical examples that illustrate the theory.
