From bc4f2b71f19c9481b91dc1c6a4952d045082bfdb Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Max=20J=C3=A4schke?= <max.jaeschke@stud.htwk-leipzig.de>
Date: Thu, 14 Mar 2024 07:40:08 +0000
Subject: [PATCH] fix missing point in bhe mesh example
---
ogstools/meshlib/gmsh_meshing.py | 16 ++++++++--------
1 file changed, 8 insertions(+), 8 deletions(-)
diff --git a/ogstools/meshlib/gmsh_meshing.py b/ogstools/meshlib/gmsh_meshing.py
index 51862b685..29cd28ed7 100644
--- a/ogstools/meshlib/gmsh_meshing.py
+++ b/ogstools/meshlib/gmsh_meshing.py
@@ -142,15 +142,15 @@ def bhe_mesh(
delta_2 = 0.866 * delta
delta_3 = 0.5 * delta
- geo.addPoint(x_BHE, y_BHE, 0, delta)
- d1 = geo.addPoint(x_BHE, y_BHE - delta, 0, delta)
- d2 = geo.addPoint(x_BHE, y_BHE + delta, 0, delta)
- d3 = geo.addPoint(x_BHE + delta_2, y_BHE + delta_3, 0, delta)
- d4 = geo.addPoint(x_BHE - delta_2, y_BHE + delta_3, 0, delta)
- d5 = geo.addPoint(x_BHE + delta_2, y_BHE - delta_3, 0, delta)
- d6 = geo.addPoint(x_BHE - delta_2, y_BHE - delta_3, 0, delta)
+ d1 = geo.addPoint(x_BHE, y_BHE, 0, delta)
+ d2 = geo.addPoint(x_BHE, y_BHE - delta, 0, delta)
+ d3 = geo.addPoint(x_BHE, y_BHE + delta, 0, delta)
+ d4 = geo.addPoint(x_BHE + delta_2, y_BHE + delta_3, 0, delta)
+ d5 = geo.addPoint(x_BHE - delta_2, y_BHE + delta_3, 0, delta)
+ d6 = geo.addPoint(x_BHE + delta_2, y_BHE - delta_3, 0, delta)
+ d7 = geo.addPoint(x_BHE - delta_2, y_BHE - delta_3, 0, delta)
geo.synchronize()
- model.mesh.embed(0, [d1, d2, d3, d4, d5, d6], 2, 1)
+ model.mesh.embed(0, [d1, d2, d3, d4, d5, d6, d7], 2, 1)
soil_1 = geo.extrude([(2, 1)], 0, 0, -depth / 2, [6], [1], True)
soil_2 = geo.extrude([soil_1[0]], 0, 0, -depth / 2, [6], [1], True)
--
GitLab