_RL Tice (nlyr)

      Tice(1) = 0.5 _d 0 *(-b1 - SQRT(b1*b1-4. _d 0*a1*c1))/a1

      Tice(2) = (Lfresh-qicen(2)) / cpIce

      IF (Tice(1).GT.0. _d 0 ) THEN

     &      ' BBerr: i,j,bi,bj,Tice = ',i,j,bi,bj,Tice

      IF ( Tice(2).GT.0. _d 0) THEN

     &      ' BBerr: i,j,bi,bj,Tice = ',i,j,bi,bj,Tice

     &  'ThSI_SOLVE4T: k, Ts, Tice=',0,Tsf,Tice

     &      (rhoi*cpIce*Tice(1)+rhoi*Lfresh*Tmlt1/Tice(1))

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

         Tice(1) = -(b1 + SQRT(b1*b1-4. _d 0*a1*c1))/(2. _d 0*a1)

         dTsf = (flx0 + k12*(Tice(1)-Tsf)) / (k12-df0dT)

     &        'ThSI_SOLVE4T: k,ts,t1,dTs=', k,Tsf,Tice(1),dTsf

            Tice(1) = (-b1 - SQRT(b1*b1-4. _d 0*a1*c1))/(2. _d 0*a1)

     &     'ThSI_SOLVE4T: k,ts,t1,dTs=', k,Tsf,Tice(1),dTsf

      Tice(2) = (2. _d 0*dt*k32*(Tice(1)+2. _d 0*Tf)

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

     &  'ThSI_SOLVE4T: k, Ts, Tice=',k,Tsf,Tice

      fct    = k12*(Tsf-Tice(1))

      flxCnB(i,j) = 4. _d 0*kIce *(Tice(2)-Tf)/hi

      qicen(1) = -cpWater*Tmlt1 + cpIce *(Tmlt1-Tice(1))

     &            + Lfresh*(1. _d 0-Tmlt1/Tice(1))

      qicen(2) = -cpIce *Tice(2) + Lfresh

      IF (Tice(1) .GE. Tmlt1) THEN

     &   ' BBerr - Bug: IceT(1) > Tmlt',i,j,bi,bj,Tice(1),Tmlt1

      IF (Tice(2) .GE. 0. _d 0) THEN

     &   ' BBerr - Bug: IceT(2) > 0',i,j,bi,bj,Tice(2)

          tIc1(i,j)  = Tice(1)

          tic2(i,j)  = Tice(2)

     &                              Tsf, Tice, dTsf