From c91c62dd6ef5da593d81bac49049e15b17d9ed43 Mon Sep 17 00:00:00 2001
From: HenrikAsmuth <henrik.asmuth@geo.uu.se>
Date: Mon, 3 Oct 2022 14:53:22 +0200
Subject: [PATCH] Add domain dimensions to config

---
 apps/gpu/LBM/BoundaryLayer/BoundaryLayer.cpp | 10 ++++++----
 1 file changed, 6 insertions(+), 4 deletions(-)

diff --git a/apps/gpu/LBM/BoundaryLayer/BoundaryLayer.cpp b/apps/gpu/LBM/BoundaryLayer/BoundaryLayer.cpp
index 969d3381a..779027761 100644
--- a/apps/gpu/LBM/BoundaryLayer/BoundaryLayer.cpp
+++ b/apps/gpu/LBM/BoundaryLayer/BoundaryLayer.cpp
@@ -114,11 +114,12 @@ void multipleLevel(const std::string& configPath)
 
     LbmOrGks lbmOrGks = LBM;
 
-    const real H = 1000.0; // boundary layer height in m
 
-    const real L_x = 6*H;
-    const real L_y = 4*H;
-    const real L_z = 1*H;
+    real H = 1000.0; // boundary layer height in m
+    const real L_x = config.contains("L_x")? config.getValue<real>("L_x"): 6*H;
+    const real L_y = config.contains("L_y")? config.getValue<real>("L_y"): 4*H;
+    const real L_z = config.contains("L_z")? config.getValue<real>("L_z"): 1*H;
+    H = L_z;
 
     const real z0  = 0.1; // roughness length in m
     const real u_star = 0.4; //friction velocity in m/s
@@ -179,6 +180,7 @@ void multipleLevel(const std::string& configPath)
     bool useStreams = (nProcs > 1 ? true: false);
     para->setUseStreams(useStreams);
     para->setMainKernel("TurbulentViscosityCumulantK17CompChim");
+    // para->setMainKernel("CumulantK17CompChimRedesigned");
     para->setIsBodyForce( config.getValue<bool>("bodyForce") );
 
     para->setTimestepStartOut(uint(tStartOut/dt) );
-- 
GitLab