Performance Improvement: shock capturing on curvilinear meshes

There is a performance improvement if the metric terms of the subcell FV method are only computed once at the beginning of the simulation, instead of at every Runge-Kutta stage after #721 is merged.