From f1e415266d726e6c6333f56ccc0991d25afe89e8 Mon Sep 17 00:00:00 2001
From: Andrew Heather <a.heather@opencfd.co.uk>
Date: Thu, 19 Jan 2017 13:30:59 +0000
Subject: [PATCH] BUG: wave models - corrected setting of U for
 shallowWaterAbsorption.  Fixes #388

---
 .../derived/shallowWaterAbsorption/shallowWaterAbsorption.C | 6 +++++-
 1 file changed, 5 insertions(+), 1 deletion(-)

diff --git a/src/waveModels/waveAbsorptionModels/derived/shallowWaterAbsorption/shallowWaterAbsorption.C b/src/waveModels/waveAbsorptionModels/derived/shallowWaterAbsorption/shallowWaterAbsorption.C
index 471ac164bcf..d6c7de60407 100644
--- a/src/waveModels/waveAbsorptionModels/derived/shallowWaterAbsorption/shallowWaterAbsorption.C
+++ b/src/waveModels/waveAbsorptionModels/derived/shallowWaterAbsorption/shallowWaterAbsorption.C
@@ -61,7 +61,11 @@ void Foam::waveModels::shallowWaterAbsorption::setVelocity
     const scalarField& level
 )
 {
-    U_ = vector::zero;
+    // Apply zero-gradient condition to z-component of velocity only
+    const volVectorField& U = mesh_.lookupObject<volVectorField>(UName_);
+    U_ = U.boundaryField()[patch_.index()].patchInternalField();
+    U_.replace(0, 0);
+    U_.replace(1, 0);
 }
 
 
-- 
GitLab