Added Roes Riemann solver

We only seem to use `q*q`; why are we taking the square root?