Just pre mDGsweep.f90 modification
Migrated modifications made to 1D code to 2D. Should be more conservative now.
Commiting just before correcting mistakes in DG updates-- update should use (rho*u)*phi rather than u*(rho*phi) for fluxes!
Going to revert changes back to pre-mass consistency
Still working on mass consistent split. Having trouble with deformation flows when velocities reverse sign