diff --git a/ogstools/meshplotlib/plot_setup_defaults.py b/ogstools/meshplotlib/plot_setup_defaults.py
index 816693760306ff4846eeb52fb5b31acded0cdccb..29810fa350d783dfc11ac55f1a018040c8d492c1 100644
--- a/ogstools/meshplotlib/plot_setup_defaults.py
+++ b/ogstools/meshplotlib/plot_setup_defaults.py
@@ -34,12 +34,15 @@ setup_dict = {
     "custom_cmap": None,
     "cmap_dict": {
         "displacement": "Greens",
-        "temperature": ["Blues", "plasma"],
+        "temperature": "plasma",
         "pressure": "Blues",
         "velocity": "coolwarm",
         "MaterialIDs": "tab20",
     },
-    "cmap_dict_if_bilinear": {"displacement": "PRGn"},
+    "cmap_dict_if_bilinear": {
+        "displacement": "PRGn",
+        "temperature": ["Blues", "plasma"],
+    },
     "cmap_if_mask": ["lightgrey", "green"],
     "rcParams": {
         "font.weight": "normal",
diff --git a/ogstools/propertylib/presets.py b/ogstools/propertylib/presets.py
index 8712f01803d986453516f34c3626e025c0eeaa07..93bc44b6bdf9a0a6dc56e7fe0b968415237a83a3 100644
--- a/ogstools/propertylib/presets.py
+++ b/ogstools/propertylib/presets.py
@@ -22,7 +22,7 @@ hydraulic_height = Scalar("pressure", "m", "m", "hydraulic_height", H_mask)
 qp_ratio = Scalar("sigma", "Pa", "percent", "QP_ratio", M_mask, v2s.qp_ratio)
 strain = Matrix("epsilon", "", "percent", "strain", M_mask)
 stress = Matrix("sigma", "Pa", "MPa", "stress", M_mask)
-temperature = Scalar("temperature", "K", "°C", mask=T_mask)
+temperature = Scalar("temperature", "K", "°C", mask=T_mask, bilinear_cmap=True)
 velocity = Vector("velocity", "m/s", "m/s", "darcy_velocity", H_mask)
 von_mises_stress = Scalar("sigma", "Pa", "MPa", "von_Mises_stress", M_mask, v2s.von_mises)
 # fmt: on