Merge pull request #2 from Ivanderkelen/dynlakes_master_notools

Adjust ice to liquid density in ice mass calculation Adjust ice to liquid density in ice mass calculation, as lake layer is not adjusted for change in density.
glemieux · Aug 28, 2020 · 6ae44c3 · 6ae44c3
2 parents 2a8c23c + ffc879f
commit 6ae44c3
Showing 1 changed file with 3 additions and 2 deletions.
diff --git 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
@@ -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