diff --git a/src/fluid/RiemannSolver/riemannSolver.hpp b/src/fluid/RiemannSolver/riemannSolver.hpp
index 4b306f38..abb70eed 100644
--- a/src/fluid/RiemannSolver/riemannSolver.hpp
+++ b/src/fluid/RiemannSolver/riemannSolver.hpp
@@ -62,6 +62,8 @@ class RiemannSolver {
   template<int dir>
   ExtrapolateToFaces<Phys, dir>* GetExtrapolator();
 
+  std::unique_ptr<ShockFlattening<Phys>> shockFlattening;
+
  private:
   template <typename P, int dir, PLMLimiter L, int O>
   friend class ExtrapolateToFaces;
@@ -75,8 +77,6 @@ class RiemannSolver {
 
   Solver mySolver;
 
-  std::unique_ptr<ShockFlattening<Phys>> shockFlattening;
-
   // Because each direction is a different template, we can't use
   std::unique_ptr<ExtrapolateToFaces<Phys,IDIR>> slopeLimIDIR;
   std::unique_ptr<ExtrapolateToFaces<Phys,JDIR>> slopeLimJDIR;