IC1 = IC + 1

      IF (IC .LT. KM1) THEN

      HOL(I,IC) = TX2(I)

      QS1(I)    = ALF(I,IC) + BET(I,IC)*POI(I,IC)

      QOL(I,IC) = MIN(QS1(I)*RHMAX,QOI(I,IC)) 

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

      HOL(I,IC) = TEM1 + QOL(I,IC) * ALHL

      HST(I,IC) = TEM1 + QS1(I)    * ALHL

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

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

      TX2(I) = HOL(I,K)  - HST(I,IC)

      DO 160 L=IC,KM1

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

      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))

      DO 260 L=KM1,IC,-1

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

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

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

      DO 520 L=IC,K

      TX8(I) = HST(II,IC) - HOL(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)

      TX3(I) = HOL(II,IC)

      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)

      DO 635 L=KM1,IC,-1

      DO L=KM1,IC,-1

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

      TX1(I) = HOL(I,IC)

      TX1(I) = TX1(I) + (TX5(I)-TX3(I)+QOL(II,IC)*ALHL)*(PRI(II,IC)/CP)

      HOL(I,IC) = TX1(I) - TX2(I) * ALBCP * PRI(I1(I),IC)

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

     *                                * TX1(I) / PRH(II,IC)

     *       * (TX7(I) - TX8(I) - TX5(I) - QOL(I1(I),IC)*ALHL)

      GMH(I,IC) = TX1(I) + PRI(I1(I),IC) * ONEBCP

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

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

     *       + 0.5*(PRS(II,IC+2) - PRS(II,IC)) * (ONE-TX6(I))

      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

         HOL(I,IC)  = ZERO

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

         GMH(I,IC)  = ZERO

      DO L = IC,K

      DO 990 L=IC,K

         TEM    = half*TX5(I) * PRI(II,IC)

      ucu(II,IC,nt) = TX1(I)