diff --git a/column/ice_therm_vertical.F90 b/column/ice_therm_vertical.F90
index f2f3cc4427d1..d2ad23be29f3 100644
--- a/column/ice_therm_vertical.F90
+++ b/column/ice_therm_vertical.F90
@@ -964,13 +964,20 @@ subroutine init_vertical_profile(nilyr,    nslyr,    &
          endif                  ! tice_high
 
          if (tice_low .and. heat_capacity) then
+
                write(warning,*) ' '
                call add_warning(warning)
-               write(warning,*) 'Starting thermo T < Tmin, layer', k
+               write(warning,*) 'Starting thermo, zTin < Tmin, layer', k
+               call add_warning(warning)
+               write(warning,*) 'k:', k
+               call add_warning(warning)
+               write(warning,*) 'zTin =',zTin(k),', Tmin=',Tmin
                call add_warning(warning)
-               write(warning,*) 'zTin =', zTin(k)
+               write(warning,*) 'zSin =',zSin(k)
+               call add_warning(warning)
+               write(warning,*) 'hin =',hin
                call add_warning(warning)
-               write(warning,*) 'Tmin =', Tmin
+               write(warning,*) 'zqin =',zqin(k)
                call add_warning(warning)
                l_stop = .true.
                stop_label = "init_vertical_profile: Starting thermo, zTin < Tmin, layer"