diff --git a/ogstools/meshlib/gmsh_meshing.py b/ogstools/meshlib/gmsh_meshing.py
index 51862b6852226dccfb7f7d29687a80c1e8160ab2..29cd28ed76535768cdc2b0d939106678fb38be34 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)