SRHD: Fix FV slope limiter method to use up-to-date primitives in evolution
Due to the cons2prim methods being tied to D, S, tau registered as dynamic_variables, we need some juggling to get the primitives of the current step.
This inconvenience is due to
- not (yet) separating the time stepping from the flux computation for the FV methods,
- and the 'hard coding' of the primitives to the
cons2primmethods.
Side note: We face the same issue with the GRHD version.
This PR also adds some more reftests to cover the blastwave1 exercise.
To this end we also updated the FV solver to work with non-periodic bdry conditions with bdry values taken from initial data.