Implement exact Riemann solver as hydro option

[maxpkatz]

It would be interesting to have an option to run the exact Riemann solver in the main Castro code, so we could more easily compare both accuracy and performance for any particular problem setup between the approximate/HLL solvers and the exact solver.