I                         hi, hs, Tsf, age, 

      _RL  hi                  ! ice height

      albice = albIceMax + (albIceMin-albIceMax)*exp(-hi/hAlbIce)

     I                     useBlkFlx, flxExcSw, Tf, hi, hs, 

      if ( hi.lt.himin ) then

      fswocn = fswpen * exp(-ksolar*hi)

      k12 = 4. _d 0*kice*ksnow / (ksnow*hi + 4. _d 0*kice*hs)

      k32 = 2. _d 0*kice  / hi

      a10 = rhoi*cpice *hi/(2. _d 0*dt) +

     &      k32 * (4. _d 0*dt*k32 + rhoi*cpice *hi)

     &       / (6. _d 0*dt*k32 + rhoi*cpice *hi)

      b10 = -hi*

     &      - k32 * (4. _d 0*dt*k32*Tf+rhoi*cpice *hi*Tice(2))

     &       / (6. _d 0*dt*k32 + rhoi*cpice *hi) - fswint

      c1 = rhoi*Lfresh*hi*Tmlt1 / (2. _d 0*dt)

            write (6,*) 'BB: thermw conv err, iceheight ', hi

     &        + rhoi*cpice *hi*Tice(2))

     &         /(6. _d 0*dt*k32 + rhoi*cpice *hi)

      flxCnB = 4. _d 0*kice *(Tice(2)-Tf)/hi

      _RL hi

     U                     compact, hi, hs, Tsf, qicen, qleft, 

        if (hi.gt.hihig.and. compact.le.iceMaskmax) then

      mwater0 = rhos*hs + rhoi*hi

      msalt0  = rhoi*hi*saltice 

      if (hi.lt.himin0) then

      if (hi.gt.hihig) then

      hlyr = hi / rnlyr

      hi = 0. _d 0

         hi = hi + hnew(k)

     &                                         etop, ebot, hi, hs    

      if ( hi.lt.himin .AND. (hi+hs).gt.0. _d 0 ) then

         hi = 0. _d 0

      fresh = (mwater0 - (rhos*hs + rhoi*hi))/dt

      IF ( hi .LE. 0. _d 0 ) THEN 

      if (hi.gt.0. _d 0.and.evap.gt.0. _d 0) then

      hi = 0. _d 0

         hi = hi + hnew(k)

      if ( hi.gt.0. _d 0 .AND. hi.lt.himin ) then

         fresh = fresh + (rhos*hs + rhoi*hi)/dt

         hi = 0. _d 0

      IF ( hi .le. 0. _d 0 ) GOTO 200

      if ( hs .gt. hi*rhoiw/rhos ) then

         dhi = (hs * rhos - hi * rhoiw) / rhosw

         hi = hi + dhi

      if (hi.gt.hiMax) then

         chi=hi-hiMax

      hi = 0. _d 0

         hi = hi + hnew(k)

       hlyr = hi/rnlyr

     &                  compact,hi,(qicen(1)+qicen(2))*0.5, hs

      if (hi.le.0. _d 0) compact=0. _d 0

      fsalt = (msalt0 - rhoi*hi*saltice)/dt

     &   (mwater0-(rhos*hs+rhoi*hi))/dt,evap,fresh,fsalt

         rqh = rq * hi + rs * hs

         freshe=(rhos*hs+rhoi*hi)/dt

         salte=(rhoi*hi*saltice)/dt

           hi=0. _d 0 

      _RL  hi