From 2c40807a2cb8b65ee1ba6b807bf4c0d29e5fc637 Mon Sep 17 00:00:00 2001
From: Wenqing Wang <wenqing.wang@ufz.de>
Date: Wed, 10 Aug 2016 18:25:25 +0200
Subject: [PATCH] [GW] Corrected the sign of the groundwater flow velocity

---
 ProcessLib/GroundwaterFlow/GroundwaterFlowFEM.h | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/ProcessLib/GroundwaterFlow/GroundwaterFlowFEM.h b/ProcessLib/GroundwaterFlow/GroundwaterFlowFEM.h
index edf90943fdf..45ddbeb8dfd 100644
--- a/ProcessLib/GroundwaterFlow/GroundwaterFlowFEM.h
+++ b/ProcessLib/GroundwaterFlow/GroundwaterFlowFEM.h
@@ -100,7 +100,7 @@ public:
                                  sm.detJ * wp.getWeight();
 
             // Darcy velocity only computed for output.
-            auto const darcy_velocity = (k * sm.dNdx *
+            auto const darcy_velocity = -(k * sm.dNdx *
                 Eigen::Map<const NodalVectorType>(local_x.data(), ShapeFunction::NPOINTS)
                 ).eval();
 
-- 
GitLab