diff --git a/columnphysics/icepack_therm_vertical.F90 b/columnphysics/icepack_therm_vertical.F90
index 98acbdcb9..78adc7f26 100644
--- a/columnphysics/icepack_therm_vertical.F90
+++ b/columnphysics/icepack_therm_vertical.F90
@@ -1711,21 +1711,12 @@ subroutine thickness_changes (nilyr,     nslyr,    &
                fhocnn = fhocnn &
                       + zqsn(k)*hsn/(real(nslyr,kind=dbl_kind)*dt)
                zqsn(k) = -rhos*Lfresh
-!tcx, tcraig, in columnphysics, this is 
-! is it correct that now everything is "if snwgrain"?
-!               if (tr_snow) then
-!                 meltsliq = meltsliq + smicetot(k)  ! add to meltponds
-!                 smice(k) = rhos
-!                 smliq(k) = c0
-!               endif
-!               if (tr_rsnw) rsnw(k) = rsnw_fall
                if (snwgrain) then
                   meltsliq = meltsliq + massice(k)  ! add to meltponds
                   smice(k) = rhos
                   smliq(k) = c0
                   rsnw(k) = rsnw_fall
                endif
-!tcx
             enddo
             melts = melts + hsn
             hsn = c0