diff --git a/dune/fem-dg/models/modelwrapper.hh b/dune/fem-dg/models/modelwrapper.hh
index 71d5e337abe7cd91ba37f6a5af63416e171aa2d9..de5f5ac459d831b4c17760363a58e0657e935fd3 100644
--- a/dune/fem-dg/models/modelwrapper.hh
+++ b/dune/fem-dg/models/modelwrapper.hh
@@ -374,10 +374,6 @@ namespace Fem
                      RangeType& jump) const
     {
       jump = AdditionalType :: jump( it, x, uLeft, uRight );
-      // take pressure as shock detection values
-      //const RangeFieldType pl = pressure( uLeft );
-      //const RangeFieldType pr = pressure( uRight );
-      //jump  = (pl-pr)/(0.5*(pl+pr));
     }
 
     // calculate jump between left and right value
diff --git a/dune/fem-dg/solver/dg.hh b/dune/fem-dg/solver/dg.hh
index e9d4c9f59b3f83a72219234cbf8a466dfa4157ec..5f151ab65239492ae962b83719cbec926c2f753c 100644
--- a/dune/fem-dg/solver/dg.hh
+++ b/dune/fem-dg/solver/dg.hh
@@ -30,7 +30,7 @@ namespace Fem
              class Additional,
              Solver::Enum solverId             = Solver::Enum::fem,
              Formulation::Enum formId          = Formulation::Enum::primal,
-             AdvectionLimiter::Enum limiterId  = AdvectionLimiter::Enum::unlimited,
+             AdvectionLimiter::Enum limiterId  = AdvectionLimiter::Enum::limited,
              AdvectionFlux::Enum advFluxId     = AdvectionFlux::Enum::llf,
              DiffusionFlux::Enum diffFluxId    = DiffusionFlux::Enum::primal
              >
diff --git a/pydemo/euler/shock_tube.py b/pydemo/euler/shock_tube.py
index 3d82c379987a1f86aef1d953e2d9f90d763bf564..b3838c97f681e6e07fb4c1026df91c2046e71cdc 100644
--- a/pydemo/euler/shock_tube.py
+++ b/pydemo/euler/shock_tube.py
@@ -11,7 +11,7 @@ from euler import Model,Sod
 
 parameter.append({"fem.verboserank": -1})
 
-grid = structuredGrid([-1, 0], [1, 0.1], [200, 10])
+grid = structuredGrid([-1, 0], [1, 0.1], [20, 5])
 dimR      = 4
 t = 0
 dt = 1e-5
@@ -58,7 +58,7 @@ def useGalerkinOp():
 def useODESolver():
     global count, t, dt, saveTime
     spaceName = "dgonb"
-    polOrder = 0
+    polOrder = 2
     space = create.space(spaceName, grid, order=polOrder, dimrange=dimR)
     u_h   = space.interpolate(Sod(), name='u_h')
     u_h_n = u_h.copy(name="previous")
@@ -71,6 +71,7 @@ def useODESolver():
         u_h_n.assign(u_h)
         operator(u_h_n, u_h)
         dt = operator.deltaT()
+        print('dt = ', dt )
         t += dt
         if t > saveTime:
             count += 1