diff --git a/columnphysics/icepack_therm_vertical.F90 b/columnphysics/icepack_therm_vertical.F90
index 8f264155..ea7ddca0 100644
--- a/columnphysics/icepack_therm_vertical.F90
+++ b/columnphysics/icepack_therm_vertical.F90
@@ -2932,6 +2932,8 @@ subroutine icepack_step_therm1(dt,                        &
 
          endif
 
+         if (present(dsnown      )) dsnown(n)    = l_dsnown
+
       enddo                  ! ncat
 
       !-----------------------------------------------------------------
@@ -2957,6 +2959,7 @@ subroutine icepack_step_therm1(dt,                        &
 
       if (present(meltsliqn   )) meltsliqn    = l_meltsliqn
       if (present(meltsliq    )) meltsliq     = l_meltsliq
+      if (present(dsnow       )) dsnow        = l_dsnow
 
       !-----------------------------------------------------------------
       ! Calculate ponds from the topographic scheme