[exercise-coupling-ff-pm]{turbulence] Correct boundary values

exercises/exercise-coupling-ff-pm/README.md

@@ -357,7 +357,8 @@ Instead of computing the whole cross-section of a channel, you can use symmetric
 values.setAllSymmetry();
 ```
 
-In addition, you have to remove the condition `onUpperBoundary_(globalPos)` from the `isOnWall(globalPos)` method.
+In addition, you have to remove the condition `onUpperBoundary_(globalPos)` from the `isOnWall(globalPos)`
+and `initialAtPos(globalPos)` method.
 
 __Task C__:
 

exercises/exercise-coupling-ff-pm/turbulence/ex_turbulence_ffproblem.hh

@@ -276,7 +276,7 @@ public:
         values[Indices::velocityXIdx] = refVelocity();
         values[Indices::temperatureIdx] = refTemperature();
 
-        if(onUpperBoundary_(globalPos))
+        if(onUpperBoundary_(globalPos) || onLowerBoundary_(globalPos))
             values[Indices::velocityXIdx] = 0.0;
 
         return values;

exercises/solution/exercise-coupling-ff-pm/turbulence/ex_turbulence_ffproblem.hh

@@ -307,8 +307,13 @@ public:
         values[Indices::velocityXIdx] = refVelocity();
         values[Indices::temperatureIdx] = refTemperature();
 
-        if(onUpperBoundary_(globalPos))
+#if EXNUMBER >= 2
+        if(onLowerBoundary_(globalPos))
+            values[Indices::velocityXIdx] = 0.0;
+#else
+        if(onUpperBoundary_(globalPos) || onLowerBoundary_(globalPos))
             values[Indices::velocityXIdx] = 0.0;
+#endif
 
         return values;
     }