Implement FV scheme for SRHD in 1d
I added a simple wave ref test. If we evolve that till the shock forms than this very vanilla FV method also struggles to nicely resolve it. So I guess we won't get around implement a WENO scheme to have a nice reference for comparison.