Intro.md

@@ -320,11 +320,7 @@ void SWE::SWE_ADERDG::boundaryValues(const double* const x,const double t,const
   out.p(0)   = in.p(0);
   out.p(1)   = in.p(1);
   out.p(direction)= -in.p(direction);
-
-  fluxout.h() = fluxin.h();
-  fluxout.p(0) = fluxin.p(0);
-  fluxout.p(1) = fluxin.p(1);
-}
+  \\....
 ```
 
 To get the new flux we apply the inverted velocity vector to the flux term (in this example for normalNonZero=0):
@@ -345,9 +341,9 @@ which results in the implementation
   ReadOnlyVariables f_in (fluxIn) ;
   Variables f_out(fluxOut);
 
-  f_out.h() = -f_in.h();
-  f_out.p(normalNonZero) =  f_in.p(normalNonZero);
-  f_out.p(normalZero)    = -f_in.p(normalZero);
+  f_out.h()  =  f_in.h();
+  f_out.p(0) =  f_in.p(0);
+  f_out.p(1) =  f_in.p(1);
 }
 ```