hi  = hIce(i,j)

     I                  iceMask, hIce, hSnow, tSrf, ageSnw,

      _RL hIce   (siLo:siHi,sjLo:sjHi)

         fsalt(i,j) = (msalt0(i,j) - rhoi*hIce(i,j)*saltIce)/dt

     &         (mwater0(i,j)-(rhos*hSnow(i,j)+rhoi*hIce(i,j)))/dt,

          rqh = rq * hIce(i,j) + rs * hSnow(i,j)

          freshe=(rhos*hSnow(i,j)+rhoi*hIce(i,j))/dt

          salte=(rhoi*hIce(i,j)*saltIce)/dt

           hIce(i,j)  =0. _d 0

     I        iceMask(i,j), icFrac(i,j), hIce(i,j), hSnow(i,j),

      _RL hIce   (1-OLx:sNx+OLx,1-OLy:sNy+OLy)

     U                  icFrac, hIce, hSnow, tSrf, qIc1, qIc2,

         xxMlt = aMlt*(hIce(i,j)-hMlt)

         xxFrz = aFrz*(hIce(i,j)-hFrz)

           IF (hIce(i,j).GT.hThickIce) THEN

           ELSEIF (hIce(i,j).GE.hThinIce) THEN

         mwater0(i,j) = rhos*hSnow(i,j) + rhoi*hIce(i,j)

         msalt0 (i,j) = rhoi*hIce(i,j)*saltIce

         IF ( hIce(i,j).GT.hThickIce .OR. fracEnMelt.EQ.0. _d 0 ) THEN

         ELSEIF (hIce(i,j).LT.hThinIce) THEN

         hnew(i,j,k) = hIce(i,j) * rec_nlyr

         hIce(i,j) = hnew(i,j,1) + hnew(i,j,2)

     &        etop(i,j), ebot(i,j), hIce(i,j), hSnow(i,j)

         IF ( hIce(i,j).LT.hIceMin

     &        .AND. (hIce(i,j)+hSnow(i,j)).GT.0. _d 0 ) THEN

          hIce(i,j)   = 0. _d 0

     &        - (rhos*hSnow(i,j)+rhoi*hIce(i,j)))/dt

         IF ( hIce(i,j) .LE. 0. _d 0 ) THEN

         IF ( hIce(i,j) .GT. 0. _d 0 ) THEN

           IF (hIce(i,j).GT.0. _d 0.AND.evapLoc(i,j).GT.0. _d 0) THEN

         IF (hIce(i,j) .GT. 0. _d 0) THEN

          hIce(i,j) = hnew(i,j,1) + hnew(i,j,2)

          IF ( hIce(i,j).GT.0. _d 0 .AND. hIce(i,j).LT.hIceMin ) THEN

     &          + (rhos*hSnow(i,j) + rhoi*hIce(i,j))/dt

           hIce(i,j)   = 0. _d 0

         IF ( hIce(i,j) .GT. 0. _d 0 ) THEN

          IF ( hSnow(i,j) .GT. hIce(i,j)*floodFac

           dhs = (hSnow(i,j) - hIce(i,j)*floodFac) * rhoi / rhosw

           hIce(i,j)  = hIce(i,j) + dhi

          IF (hIce(i,j).GT.hiMax) THEN

           chi=hIce(i,j)-hiMax

          hIce(i,j) = hnew(i,j,1) + hnew(i,j,2)

          hlyr = hIce(i,j) * rec_nlyr

     &         icFrac(i,j),hIce(i,j), (qicen(i,j,1)+qicen(i,j,2))*0.5,

         IF (hIce(i,j).LE.0. _d 0) icFrac(i,j)=0. _d 0

     &    hIce, hSnow, snowPrc(i,j,bi,bj)*thSIce_deltaT/rhos

     I             iceFrac, compact, hIce, hSnow, qicen,

      _RL compact, hIce, hSnow, qicen(nlyr)

     &               + rhoi*hIce*(qicen(1)+qicen(2))*0.5

     &     - compact*( rhos*hSnow + rhoi*hIce )

     &     - compact* rhoi*hIce*saltice

     U                  icFrac, hIce, hSnow,

          iceThick= hIce(i,j)

          hIce(i,j)   = iceThick

      _RL hIce   (1-OLx:sNx+OLx,1-OLy:sNy+OLy)

     I                  iceMask, hIce, hSnow, tFrz, flxExSW,

          IF ( hIce(i,j).LT.hIceMin ) THEN

          fswocn = fswpen * exp(-ksolar*hIce(i,j))

     &             / (kSnow*hIce(i,j) + 4. _d 0*kIce*hSnow(i,j))

          k32      = 2. _d 0*kIce  / hIce(i,j)

          a10(i,j) = rhoi*cpIce *hIce(i,j)/(2. _d 0*dt) +

     &          k32*( 4. _d 0*dt*k32 + rhoi*cpIce *hIce(i,j) )

     &           / ( 6. _d 0*dt*k32 + rhoi*cpIce *hIce(i,j) )

          b10(i,j) = -hIce(i,j)*

     &               +rhoi*cpIce*hIce(i,j)*tIc2(i,j) )

     &           / ( 6. _d 0*dt*k32 + rhoi*cpIce *hIce(i,j) )

          c10(i,j) = rhoi*Lfresh*hIce(i,j)*Tmlt1 / (2. _d 0*dt)

     &      ') with hIce = ', hIce(ii,jj)

     &       ' BB: not converge: i,j,it,hi=',i,j,bi,bj,myIter,hIce(i,j)

          k32       = 2. _d 0*kIce  / hIce(i,j)

     &                 + rhoi*cpIce*hIce(i,j)*tIc2(i,j))

     &               /(6. _d 0*dt*k32 + rhoi*cpIce*hIce(i,j))

          flxCnB(i,j) = 4. _d 0*kIce *(tIc2(i,j)-tFrz(i,j))/hIce(i,j)

      _RL hIce   (1-OLx:sNx+OLx,1-OLy:sNy+OLy)