From ce7324779f331c33d2d02808bee659d06c6e78df Mon Sep 17 00:00:00 2001
From: Julian Heinze <julian.heinze@ufz.de>
Date: Wed, 7 Feb 2024 11:22:49 +0100
Subject: [PATCH] bug fix - remove bulk_node/element_ids after assignment
---
ogstools/feflowlib/tools.py | 6 ++++++
1 file changed, 6 insertions(+)
diff --git a/ogstools/feflowlib/tools.py b/ogstools/feflowlib/tools.py
index 6a72d53cd..6c8c95aea 100644
--- a/ogstools/feflowlib/tools.py
+++ b/ogstools/feflowlib/tools.py
@@ -110,6 +110,9 @@ def extract_point_boundary_conditions(
dict_of_point_boundary_conditions[
str(out_mesh_path / point_data) + ".vtu"
] = filtered_points
+ # Remove bulk node/element ids from bulk mesh, as they are not needed anymore.
+ mesh.point_data.remove("bulk_node_ids")
+ mesh.cell_data.remove("bulk_element_ids")
return dict_of_point_boundary_conditions
@@ -168,6 +171,9 @@ def extract_cell_boundary_conditions(
topsurf.point_data.remove(pt_data)
# correct unit for P_IOFLOW, in FEFLOW m/d in ogs m/s
topsurf.cell_data["P_IOFLOW"] = topsurf.cell_data["P_IOFLOW"]
+ # Remove bulk node/element ids from bulk mesh, as they are not needed anymore.
+ mesh.point_data.remove("bulk_node_ids")
+ mesh.cell_data.remove("bulk_element_ids")
return (
bulk_mesh_path.with_stem("topsurface_" + bulk_mesh_path.stem),
topsurf,
--
GitLab