implementation of DG scheme for MHD equation on 2D Cartesian mesh, locally divergence-free space included