Added Roes Riemann solver

Rewrite as `sqrt(ul[0])*sqrt(ur[0]))` as we use both of these further down.