diff --git a/tutorial/tutorialproblem_coupled.hh b/tutorial/tutorialproblem_coupled.hh
index bf94e53b5462c52535646273e433fafb989197cc..c7651b62a3c9824df52d71746f0a506a253b2b68 100644
--- a/tutorial/tutorialproblem_coupled.hh
+++ b/tutorial/tutorialproblem_coupled.hh
@@ -57,8 +57,8 @@ SET_PROP(TutorialProblemCoupled, Grid) /*@\label{tutorial-coupled:set-grid}@*/
         Dune::FieldVector<int, 2> cellRes;
         Dune::FieldVector<ctype, 2> lowerLeft(0.0);
         Dune::FieldVector<ctype, 2> upperRight;
-        cellRes[0] = 30;
-        cellRes[1] = 10;
+        cellRes[0] = 100;
+        cellRes[1] = 1;
         upperRight[0] = 300;
         upperRight[1] = 60;
         return new Dune::SGrid<2,2>(cellRes,
@@ -192,10 +192,10 @@ public:
             fvElemGeom.boundaryFace[boundaryFaceIdx].ipGlobal;
         Scalar right = this->bboxMax()[0];
         // extraction of oil on the right boundary for approx. 1.e6 seconds
-        if (pos[0] > right - eps_ && this->timeManager().time() <= 1.e6) {
-            // oil outflux of 0.3 g/(m * s) on the right boundary.
+        if (pos[0] > right - eps_) {
+            // oil outflux of 30 g/(m * s) on the right boundary.
             values[Indices::contiWEqIdx] = 0;
-            values[Indices::contiNEqIdx] = 3e-4;
+            values[Indices::contiNEqIdx] = 3e-2;
         } else {
             // no-flow on the remaining Neumann-boundaries.
             values[Indices::contiWEqIdx] = 0;
diff --git a/tutorial/tutorialproblem_decoupled.hh b/tutorial/tutorialproblem_decoupled.hh
index a9137082646f8697e2af5ddb11e5e4b5c2b919e5..2398b274b19142b3e95f13de30f1bce725b7b3fd 100644
--- a/tutorial/tutorialproblem_decoupled.hh
+++ b/tutorial/tutorialproblem_decoupled.hh
@@ -184,7 +184,7 @@ public:
     bool shouldWriteOutput() const /*@\label{tutorial-decoupled:output}@*/
     {
         return this->timeManager().timeStepIndex() > 0 &&
-        (this->timeManager().timeStepIndex() % 1 == 0);
+        (this->timeManager().timeStepIndex() % 20 == 0);
     }
 
     //! Returns the temperature within the domain.
@@ -262,7 +262,7 @@ public:
         std::vector<Scalar> neumannFlux(2,0.0);
         if (globalPos[0] > this->bboxMax()[0] - eps_)
         {
-            neumannFlux[nPhaseIdx] = 3e-4;
+            neumannFlux[nPhaseIdx] = 3e-2;
         }
         return neumannFlux;
     }
diff --git a/tutorial/tutorialspatialparameters_coupled.hh b/tutorial/tutorialspatialparameters_coupled.hh
index 88c616167081d9d534fca9d05f2b8f49968bdce2..b6f62f29aa7bbeec48d6a0d2c0eb4decb89ac945 100644
--- a/tutorial/tutorialspatialparameters_coupled.hh
+++ b/tutorial/tutorialspatialparameters_coupled.hh
@@ -104,7 +104,7 @@ public:
         materialParams_.setSnr(0.0);
 
         //linear material law
-        materialParams_.setPe(1000.0);
+        materialParams_.setPe(500.0);
         materialParams_.setAlpha(2);
     }
 
diff --git a/tutorial/tutorialspatialparameters_decoupled.hh b/tutorial/tutorialspatialparameters_decoupled.hh
index 997d7dc247eff0833eeea4afe984f006719e2e0d..22479f852bd639e300d846b2ae0679e2f5930872 100644
--- a/tutorial/tutorialspatialparameters_decoupled.hh
+++ b/tutorial/tutorialspatialparameters_decoupled.hh
@@ -92,7 +92,7 @@ public:
 
         // parameters for the Brooks-Corey Law
         // entry pressures
-        materialLawParams_.setPe(1000);
+        materialLawParams_.setPe(500);
 
         // Brooks-Corey shape parameters
         materialLawParams_.setAlpha(2);