diff --git a/columnphysics/icepack_fsd.F90 b/columnphysics/icepack_fsd.F90
index 291ee6fb..2e5175af 100644
--- a/columnphysics/icepack_fsd.F90
+++ b/columnphysics/icepack_fsd.F90
@@ -592,7 +592,6 @@ subroutine fsd_lateral_growth (dt,        aice,         &
       endif ! vi0new_lat > 0
 
       ! Use remaining ice volume as in standard model,
-      ! but ice cannot grow into the area that has grown laterally
       tot_latg = SUM(d_an_latg(:))
 
       end subroutine fsd_lateral_growth
diff --git a/columnphysics/icepack_therm_itd.F90 b/columnphysics/icepack_therm_itd.F90
index 3610b7a0..0fcc9673 100644
--- a/columnphysics/icepack_therm_itd.F90
+++ b/columnphysics/icepack_therm_itd.F90
@@ -959,9 +959,6 @@ subroutine lateral_melt (dt,         fpond,      &
 
       if (.not. any(rsiden(:) > c0)) return  ! no lateral melt, get out now
 
-!      write(warnstr,*) 'LM ',rsiden(1)
-!      call icepack_warnings_add(warnstr)
-
       dfhocn   = c0
       dfpond   = c0
       dfresh   = c0
@@ -1027,7 +1024,7 @@ subroutine lateral_melt (dt,         fpond,      &
          ! floe size distribution
          if (tr_fsd) then
             if (rsiden(n) > puny) then
-               if (aicen(n) > puny) then    ! not sure if this should be aicen or aicen_init
+               if (aicen(n) > puny) then
 
                   ! adaptive subtimestep
                   elapsed_t = c0
diff --git a/columnphysics/icepack_therm_vertical.F90 b/columnphysics/icepack_therm_vertical.F90
index b5108986..bf96096e 100644
--- a/columnphysics/icepack_therm_vertical.F90
+++ b/columnphysics/icepack_therm_vertical.F90
@@ -696,8 +696,6 @@ subroutine frzmlt_bottom_lateral (dt,                 &
             rsiden(n) = rsiden(n) * xtmp  ! xtmp is almost always 1 so usually nothing happens here
          enddo ! ncat
 
-!         write(warnstr,*) 'FBM ',rsiden(1), xtmp
-!         call icepack_warnings_add(warnstr)
       endif
 
       if (present(wlat)) wlat=wlat_loc