diff --git a/columnphysics/icepack_therm_vertical.F90 b/columnphysics/icepack_therm_vertical.F90
index b162fda67..92a035037 100644
--- a/columnphysics/icepack_therm_vertical.F90
+++ b/columnphysics/icepack_therm_vertical.F90
@@ -1431,6 +1431,9 @@ subroutine thickness_changes (nilyr,     nslyr,    &
          ebot_mlt = ebot_mlt - dhs*zqsn(k)
          ebot_mlt = max(ebot_mlt, c0)
 
+         ! Add this to the snow melt (J. Zhu)
+         melts = melts - dhs
+
       enddo                     ! nslyr
 
       !-----------------------------------------------------------------