From ffc879f88203767e8fa4d70d69e0b192d73deb8d Mon Sep 17 00:00:00 2001
From: ivanderkelen <inne.vanderkelen@vub.be>
Date: Fri, 28 Aug 2020 00:50:23 -0600
Subject: [PATCH] Adjust ice density to liquid density in AccumulateHeatLake
 routine

---
 src/biogeophys/TotalWaterAndHeatMod.F90 | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/src/biogeophys/TotalWaterAndHeatMod.F90 b/src/biogeophys/TotalWaterAndHeatMod.F90
index 87d3256bbf..79ee0a77d6 100644
--- a/src/biogeophys/TotalWaterAndHeatMod.F90
+++ b/src/biogeophys/TotalWaterAndHeatMod.F90
@@ -1086,7 +1086,8 @@ subroutine AccumulateHeatLake(bounds, num_c, filter_c, &
                 heat_liquid = lake_heat_liquid(c), &
                 latent_heat_liquid = lake_latent_heat_liquid(c))
             ! ice heat
-            h2olak_ice = dz_lake(c,j) * denice * lake_icefrac(c,j)      
+            ! use water density as lake layer does not adjust
+            h2olak_ice = dz_lake(c,j) * denh2o * lake_icefrac(c,j)      
             lake_heat_ice(c) = lake_heat_ice(c) + & 
                 TempToHeat(temp=t_lake(c,j), cv = (h2olak_ice * cpice))
         end do