convert to Kelvin instead of Celcicus + change viscosity

- for all temperature related arrays

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
 
 

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,