Skip to content

Commit

Permalink
Adjust ice density to liquid density in AccumulateHeatLake routine
Browse files Browse the repository at this point in the history
  • Loading branch information
Ivanderkelen committed Aug 28, 2020
1 parent a672692 commit ffc879f
Showing 1 changed file with 2 additions and 1 deletion.
3 changes: 2 additions & 1 deletion src/biogeophys/TotalWaterAndHeatMod.F90
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down

0 comments on commit ffc879f

Please sign in to comment.