diff --git a/columnphysics/icepack_snow.F90 b/columnphysics/icepack_snow.F90
index 315ce455..1d10d95b 100644
--- a/columnphysics/icepack_snow.F90
+++ b/columnphysics/icepack_snow.F90
@@ -360,6 +360,9 @@ subroutine icepack_step_snow(dt,                   &
                hsn(n)   = vsnon(n)/aicen(n)
                hin(n)   = vicen(n)/aicen(n)
             endif
+            do k = 1, nslyr
+               rsnw (k,n) = max(rsnw_fall, rsnw(k,n))
+            enddo
          enddo
 
          call update_snow_radius (dt,                &