diff --git a/columnphysics/icepack_therm_vertical.F90 b/columnphysics/icepack_therm_vertical.F90
index 1b0fa027..53f082f6 100644
--- a/columnphysics/icepack_therm_vertical.F90
+++ b/columnphysics/icepack_therm_vertical.F90
@@ -375,12 +375,11 @@ subroutine thermo_vertical (dt,          aicen,     &
          einter = einter + hilyr * zqin(k)
       enddo ! k
 
-      Tsnice = c0
-      if ((hslyr+hilyr) > puny) then
+      if (hilyr > puny) then
          if (hslyr > puny) then
-            Tsnice = (hslyr*zTsn(nslyr) + hilyr*zTin(1)) / (hslyr+hilyr)
+            Tsnice = Tsnice + aicen*((hilyr*zTsn(nslyr) + hslyr*zTin(1)) / (hslyr+hilyr))
          else
-            Tsnice = Tsf
+            Tsnice = Tsnice + aicen*Tsf
          endif
       endif
 
diff --git a/configuration/driver/icedrv_flux.F90 b/configuration/driver/icedrv_flux.F90
index bde38972..b0509286 100644
--- a/configuration/driver/icedrv_flux.F90
+++ b/configuration/driver/icedrv_flux.F90
@@ -688,6 +688,7 @@ subroutine init_history_therm
       congel (:) = c0
       frazil (:) = c0
       snoice (:) = c0
+      Tsnice (:) = c0
       dsnow  (:) = c0
       meltt  (:) = c0
       melts  (:) = c0