From a6726928f50f10268d9659130fd4f349d9d1a34c Mon Sep 17 00:00:00 2001 From: ivanderkelen Date: Thu, 27 Aug 2020 13:50:41 -0600 Subject: [PATCH 1/2] Adjust ice to liquid density in ice mass calculation, as lake layer is not adjusted --- src/biogeophys/TotalWaterAndHeatMod.F90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/biogeophys/TotalWaterAndHeatMod.F90 b/src/biogeophys/TotalWaterAndHeatMod.F90 index b774fe9761..87d3256bbf 100644 --- a/src/biogeophys/TotalWaterAndHeatMod.F90 +++ b/src/biogeophys/TotalWaterAndHeatMod.F90 @@ -518,7 +518,7 @@ subroutine AccumulateLiqIceMassLake(bounds, num_c, filter_c, & c = filter_c(fc) ! calculate lake liq and ice content per lake layer first h2olak_liq = dz_lake(c,j) * denh2o * (1 - lake_icefrac(c,j)) - h2olak_ice = dz_lake(c,j) * denice * lake_icefrac(c,j) + h2olak_ice = dz_lake(c,j) * denh2o * lake_icefrac(c,j) ! use water density of liquid water as layer depth is not adjusted liquid_mass(c) = liquid_mass(c) + h2olak_liq ice_mass(c) = ice_mass(c) + h2olak_ice From ffc879f88203767e8fa4d70d69e0b192d73deb8d Mon Sep 17 00:00:00 2001 From: ivanderkelen Date: Fri, 28 Aug 2020 00:50:23 -0600 Subject: [PATCH 2/2] 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