From a17dabfa946594320f70ad4289feac667c95a399 Mon Sep 17 00:00:00 2001
From: robin-zweigel <57904608+robin-zweigel@users.noreply.github.com>
Date: Thu, 28 Nov 2019 14:32:53 +0100
Subject: [PATCH] Update GROUND_base_class.m

---
 modules/@GROUND_base_class/GROUND_base_class.m | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/modules/@GROUND_base_class/GROUND_base_class.m b/modules/@GROUND_base_class/GROUND_base_class.m
index 371ed518..b046ff8e 100644
--- a/modules/@GROUND_base_class/GROUND_base_class.m
+++ b/modules/@GROUND_base_class/GROUND_base_class.m
@@ -120,7 +120,7 @@
             E_frozen = -Lf.*ground.STATVAR.waterIce;
             
             ground.STATVAR.T = double(ground.STATVAR.energy < E_frozen) .* (ground.STATVAR.energy - E_frozen) ./ (c_i.*ground.STATVAR.waterIce + c_m.*ground.STATVAR.mineral + c_o.*ground.STATVAR.organic) + ...
-                double(ground.STATVAR.energy >0) .* ground.STATVAR.energy ./ (c_i.*ground.STATVAR.waterIce + c_m.*ground.STATVAR.mineral + c_o.*ground.STATVAR.organic);
+                double(ground.STATVAR.energy >0) .* ground.STATVAR.energy ./ (c_w.*ground.STATVAR.waterIce + c_m.*ground.STATVAR.mineral + c_o.*ground.STATVAR.organic);
             ground.STATVAR.ice = double(ground.STATVAR.energy <= E_frozen) .*ground.STATVAR.waterIce + double(ground.STATVAR.energy > E_frozen & ground.STATVAR.energy < 0) .* ground.STATVAR.energy ./ (-Lf);
             ground.STATVAR.water = double(ground.STATVAR.energy >= 0) .*ground.STATVAR.waterIce + double(ground.STATVAR.energy > - Lf.*ground.STATVAR.waterIce & ground.STATVAR.energy < 0) .* (ground.STATVAR.energy + Lf.*ground.STATVAR.waterIce) ./ Lf;