*,         ETA(lng,K), HOI(lng,K)

     *,           ALF, BET, GAM, PRH, PRI, HOI, ETA

     *,                 ALF, BET, GAM, PRH, PRI, HOL, ETA

      _RL HOL(LENC,K), ETA(LENC,K), HST(LENC,K), QOL(LENC,K)

      ETA(I,K) = ZERO

         ETA(I,L) = ETA(I,L+1) + TX1(I)*CPBG

      ETA(I,IC) = ETA(I,IC1) + CPBG * TX3(I)

      TX1(I) = TX1(I) + (HST(I,IC) - HOL(I,L)) * (ETA(I,L) - ETA(I,L+1))

      ETA(I,K) = 1.0

      TX1(I) = ETA(IA(I),L)

      ETA(I,L) = 1.0 + ALM(I) * TX1(I)

     *                         *  HST(II,IC) * ETA(I,IC1)

         TEM = TX1(I) + (ETA(I,L) - ETA(I,L+1)) * HOL(II,L)

         TEM1   = ETA(I,L+1) * PCU(I)

         TEM1   = (TEM1 + ETA(I,L) * PCU(I)) * HST(II,L)

      TX1(I) = ETA(I2(I),L)

      ETA(I,L) = TX1(I)

      WLQ(I) = QOL(II,K) - QS1(I)     * ETA(I,IC)

      UHT(I,NT) = UOI(II,K+nltop-1,NT)-UOI(II,IC+nltop-1,NT) * ETA(I,IC)

      TEM    = ETA(I,L) - ETA(I,L+1)

      TEM    = ETA(I,L) - ETA(I,L+1)

         TEM      = (POI(II,L-1) - POI(II,L)) * ETA(I,L)

     *            * (ETA(I,L)   * (PRH(II,L)-PRJ(II,L)) +

     *               ETA(I,L+1) * (PRJ(II,L+1)-PRH(II,L))) / PRH(II,L)

         AKM(I) = AKM(I) - ETA(I,IC1) * (PRJ(II,IC1) - PRH(II,IC)) 

         TX3(I) = (QOL(II,L-1) - QOL(II,L)) * ETA(I,L)

     *          * (TX3(I)*(ALHL*HALF) + ETA(I,IC) * TX8(I))

         TX1(I) = TX1(I) + (ETA(I,L) - ETA(I,L+1)) * GMH(I,L)

       CMASS(II,L) =  ETA(I,L+1) * WFN(I) * ONEBG

      TX1(I) = (UOI(II,L-1+nltop-1,nt)-UOI(II,L+nltop-1,nt)) * ETA(I,L)

     &       ETA,ZETA,DRAGS,DRAGA,AMASS,FORCEX,FORCEY,UICEC,VICEC

      _RL ETA        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,  nSx,nSy)

          ETA(I,J,bi,bj)=ECM2*ZETA(I,J,bi,bj)

      _EXCH_XY_R8(ETA, myThid)

          ETA(I,J,bi,bj)=ECM2*ZETA(I,J,bi,bj)

      _EXCH_XY_R8(ETA, myThid)

          ETAMEAN(I,J,bi,bj) =QUART*(ETA(I,J-1,bi,bj)+ETA(I-1,J-1,bi,bj)

     &                 +ETA(I,J,bi,bj)+ETA(I-1,J,bi,bj))

          AA1=((ETA(I,J-1,bi,bj)+ZETA(I,J-1,bi,bj))

     &           +(ETA(I,J,bi,bj)+ZETA(I,J,bi,bj))

          AA2=((ETA(I-1,J-1,bi,bj)+ZETA(I-1,J-1,bi,bj))

     &           +(ETA(I-1,J,bi,bj)+ZETA(I-1,J,bi,bj))

          AA3=ETA(I-1,J,bi,bj)+ETA(I,J,bi,bj)

          AA4=ETA(I-1,J-1,bi,bj)+ETA(I,J-1,bi,bj)

          AA5=-(ETA(I-1,J,bi,bj)+ETA(I,J,bi,bj)-ETA(I-1,J-1,bi,bj)

     &           -ETA(I,J-1,bi,bj))*TNGICE(I,J,bi,bj)

     4            -HALF*(ETA(I,J,bi,bj)*(VICEC(I+1,J+1,bi,bj)

     4            +ETA(I,J-1,bi,bj)*(VICEC(I+1,J,bi,bj)

     4            -VICEC(I,J-1,bi,bj))+ETA(I-1,J,bi,bj)

     4            +ETA(I-1,J-1,bi,bj)*(VICEC(I,J-1,bi,bj)

     5            +HALF*(ETA(I,J,bi,bj)*RECIP_CSTICE(I-1,J,bi,bj)

     5            +ETA(I-1,J,bi,bj)*RECIP_CSTICE(I-1,J,bi,bj)

     5            +ETA(I,J-1,bi,bj)*RECIP_CSTICE(I-1,J-1,bi,bj)

     5            +ETA(I-1,J-1,bi,bj)*RECIP_CSTICE(I-1,J-1,bi,bj)

     6            +(ETA(I,J,bi,bj)+ETA(I,J-1,bi,bj)-ETA(I-1,J-1,bi,bj)

     6            -ETA(I-1,J,bi,bj)))

          UVRT1(I,J,bi,bj)=(ETA(I-1,J-1,bi,bj)+ETA(I,J-1,bi,bj))

          UVRT2(I,J,bi,bj)=(ETA(I-1,J,bi,bj)+ETA(I,J,bi,bj))

           AA2=((ETA(I-1,J-1,bi,bj)+ZETA(I-1,J-1,bi,bj))

     &            +(ETA(I-1,J,bi,bj)+ZETA(I-1,J,bi,bj))

           AA1=((ETA(I,J-1,bi,bj)+ZETA(I,J-1,bi,bj))

     &            +(ETA(I,J,bi,bj)+ZETA(I,J,bi,bj))

          AA1=ETA(I-1,J,bi,bj)+ZETA(I-1,J,bi,bj)+ETA(I,J,bi,bj)

          AA2=ETA(I-1,J-1,bi,bj)+ZETA(I-1,J-1,bi,bj)+ETA(I,J-1,bi,bj)

          AA3=(ETA(I,J-1,bi,bj)*RECIP_CSUICE(I,J,bi,bj)+ETA(I,J,bi,bj)

          AA4=(ETA(I-1,J-1,bi,bj)+ETA(I-1,J,bi,bj))

          AA5=((ZETA(I-1,J,bi,bj)-ETA(I-1,J,bi,bj))

     &          +(ZETA(I,J,bi,bj)-ETA(I,J,bi,bj))

     &          -(ZETA(I-1,J-1,bi,bj)-ETA(I-1,J-1,bi,bj))

     &          -ETA(I,J-1,bi,bj)))*TNGICE(I,J,bi,bj)

     4        *(ETA(I-1,J,bi,bj)+ETA(I,J,bi,bj)

     4        -ETA(I-1,J-1,bi,bj)-ETA(I,J-1,bi,bj))*DELXY(I,J,bi,bj)

     5        +HALF*(ETA(I,J,bi,bj)*(UICEC(I+1,J+1,bi,bj)

     5        -UICEC(I+1,J,bi,bj)-UICEC(I,J,bi,bj))+ETA(I,J-1,bi,bj)

     5        -UICEC(I,J-1,bi,bj))+ETA(I-1,J,bi,bj)

     5        +ETA(I-1,J-1,bi,bj)*(UICEC(I,J-1,bi,bj)

     6        +(ETA(I,J,bi,bj)+ETA(I,J-1,bi,bj)-ETA(I-1,J-1,bi,bj)

     6        -ETA(I-1,J,bi,bj))

          UVRT1(I,J,bi,bj)=(ETA(I-1,J-1,bi,bj)*RECIP_CSUICE(I,J,bi,bj)

     &            +ETA(I-1,J,bi,bj)*RECIP_CSUICE(I,J,bi,bj))

          UVRT2(I,J,bi,bj)=(ETA(I,J-1,bi,bj)*RECIP_CSUICE(I,J,bi,bj)

     &            +ETA(I,J,bi,bj)*RECIP_CSUICE(I,J,bi,bj))

            AA2=ETA(I-1,J-1,bi,bj)+ZETA(I-1,J-1,bi,bj)

     &         +ETA(I,J-1,bi,bj)+ZETA(I,J-1,bi,bj)

            AA1=ETA(I-1,J,bi,bj)+ZETA(I-1,J,bi,bj)+ETA(I,J,bi,bj)

          ETA(I,J,bi,bj)=ZETA(I,J,bi,bj)/4.0 _d 0