IC1 = IC + 1

         DO 3703 L=KM1,IC1,-1

      TEM1      = TX2(I) + PRJ(I,IC1) * POI(I,IC) * CP

      TX3(I   ) = (PRJ(I,IC1) - PRH(I,IC)) * POI(I,IC)

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

      if (IC1 .lt. K) then

         IF (TX2(I) .LE. 0.0 .AND. (HOL(I,K) .GT. HST(I,IC1))

      TX2(I)   = 0.5 * (PRS(II,IC) + PRS(II,IC1))

      IF (HST(II,IC1) .LT. HST(II,IC)) THEN

         TX6(I) = (HST(II,IC1)-HOL(II,K))/(HST(II,IC1)-HST(II,IC))

      TX2(I) = 0.5 * (PRS(II,IC1)+PRS(II,IC1+1)) * (1.0-TX6(I))

         WFN(I) = - GAM(II,IC) * (PRJ(II,IC1) - PRH(II,IC))

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

      IF (IC1 .LE. KM1) THEN

         DO 380 L=KM1,IC1,-1

         WFN(I) = WFN(I) + TX1(I) * GAM(II,IC)*(PRJ(II,IC1)-PRH(II,IC))

      TEM    = TX6(I) * (HOL(II,IC)-HOL(II,IC1)) + HOL(II,IC1)

      TEM1   = TX6(I) * (QOL(II,IC)-QOL(II,IC1)) + QOL(II,IC1)

      TX2(I)     = 0.5 * (PRS(II,IC) + PRS(II,IC1))

      IF (IC1 .LE. KM1) THEN

         DO 750 L=KM1,IC1,-1

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

      IF (IC1 .LE. KM1) THEN

         DO 840 L=KM1,IC1,-1

      IF (IC1 .LE. KM1) THEN

         DO 3725 L=KM1,IC1,-1

         IF (L .EQ. IC1) THEN

      TX2(I) = 0.5 * (PRS(II,IC) + PRS(II,IC1))

      TX1(I) = PRS(II,IC1)

      TX5(I) = 0.5 * (PRS(II,IC1) + PRS(II,IC+2))

         TEM        = PRI(II,IC1) / PRI(II,IC)

         HOL(I,IC1) = HOL(I,IC1) + HOL(I,IC) * TEM

         GMH(I,IC1) = GMH(I,IC1) + GMH(I,IC) * TEM

      DO 1020 L=KM1,IC1,-1