diff --git a/columnphysics/icepack_therm_vertical.F90 b/columnphysics/icepack_therm_vertical.F90
index 7ec9bf47..a0aa6b84 100644
--- a/columnphysics/icepack_therm_vertical.F90
+++ b/columnphysics/icepack_therm_vertical.F90
@@ -2602,7 +2602,6 @@ subroutine icepack_step_therm1(dt,                        &
          congeln(n) = c0
          snoicen(n) = c0
          l_dsnown   = c0
-         if (present(dsnown)) dsnown(n) = c0
 
          Trefn  = c0
          Qrefn  = c0
diff --git a/columnphysics/version.txt b/columnphysics/version.txt
index 7efd5d79..f46731dc 100644
--- a/columnphysics/version.txt
+++ b/columnphysics/version.txt
@@ -1 +1 @@
-ICEPACK 1.4.1
+ICEPACK 1.5.0
diff --git a/doc/source/conf.py b/doc/source/conf.py
index 49c6e2c4..ba2e0471 100644
--- a/doc/source/conf.py
+++ b/doc/source/conf.py
@@ -63,9 +63,9 @@
 # built documents.
 #
 # The short X.Y version.
-version = u'1.4.1' 
+version = u'1.5.0' 
 # The full version, including alpha/beta/rc tags.
-version = u'1.4.1' 
+version = u'1.5.0' 
 
 # The language for content autogenerated by Sphinx. Refer to documentation
 # for a list of supported languages.
diff --git a/doc/source/user_guide/interfaces.include b/doc/source/user_guide/interfaces.include
index b41f6055..b28eff1a 100644
--- a/doc/source/user_guide/interfaces.include
+++ b/doc/source/user_guide/interfaces.include
@@ -2758,12 +2758,14 @@ icepack_step_therm1
            melttn      , & ! top ice melt                           (m)
            meltbn      , & ! bottom ice melt                        (m)
            congeln     , & ! congelation ice growth                 (m)
-           snoicen     , & ! snow-ice growth                        (m)
-           dsnown          ! change in snow thickness (m/step-->cm/day)
+           snoicen         ! snow-ice growth                        (m)
   
         real (kind=dbl_kind), dimension(:), intent(in) :: &
            fswthrun        ! SW through ice to ocean            (W/m^2)
   
+        real (kind=dbl_kind), dimension(:), intent(inout), optional :: &
+           dsnown          ! change in snow thickness (m/step-->cm/day)
+  
         real (kind=dbl_kind), dimension(:), intent(in), optional :: &
            fswthrun_vdr , & ! vis dir SW through ice to ocean   (W/m^2)
            fswthrun_vdf , & ! vis dif SW through ice to ocean   (W/m^2)