1D tested and working code, 2D updated Riemann solver 2D