From 1be7b7f9f66d5dad757537146d564d9ec3ff3f42 Mon Sep 17 00:00:00 2001
From: Julian Heinze <julian.heinze@ufz.de>
Date: Mon, 15 Jan 2024 13:24:01 +0100
Subject: [PATCH] convert to Kelvin instead of Celcicus + change viscosity -
for all temperature related arrays
---
ogstools/feflowlib/_feflowlib.py | 4 +++-
ogstools/feflowlib/tools.py | 2 +-
2 files changed, 4 insertions(+), 2 deletions(-)
diff --git a/ogstools/feflowlib/_feflowlib.py b/ogstools/feflowlib/_feflowlib.py
index b0334d8af..921870f59 100644
--- a/ogstools/feflowlib/_feflowlib.py
+++ b/ogstools/feflowlib/_feflowlib.py
@@ -244,8 +244,10 @@ def _convert_to_SI_units(mesh: pv.UnstructuredGrid):
for to_be_converted in arrays_to_be_converted
):
mesh[data] *= 1 / 86400
- if "4TH" in data or "2ND" in data:
+ elif "4TH" in data or "2ND" in data:
mesh[data] *= -1 / 86400
+ elif "HEAT" in data or "TEMP" in data:
+ mesh[data] = mesh[data] + [273.15] * len(mesh[data])
return mesh
diff --git a/ogstools/feflowlib/tools.py b/ogstools/feflowlib/tools.py
index 67e0e0498..2a879d6bc 100644
--- a/ogstools/feflowlib/tools.py
+++ b/ogstools/feflowlib/tools.py
@@ -478,7 +478,7 @@ def materials_in_HT(
phase_type="AqueousLiquid",
name="viscosity",
type="Constant",
- value=1,
+ value=0.001124,
)
model.media.add_property(
medium_id=material_id,
--
GitLab