IF( EVAP(CHNO) .GT. EMAXRT(CHNO) ) THEN

        DEA = EM(CHNO) - EA(CHNO)

     &   (Q0 + ALHX(CHNO)*(ETURB(ChNo)-EMAXRT(CHNO)) - DHSDEA(CHNO)*DEA)

     &          /  ( A11 - ALHX(CHNO)*DEDTC(ChNo) )

        EVAP(CHNO) = EMAXRT(CHNO)

        SMELT(CHNO)=0.

      ESATNW = ESATTC(CHNO)+DESDTC(CHNO)*DTC

      EANEW = EA(CHNO) + DEA

      IF( EVAP(CHNO) .LT. 0. .AND. EM(CHNO).LT.ESATNW )  THEN

        DEA = EM(CHNO) - EA(CHNO)

        DTC = ( Q0 + ALHX(CHNO)*ETURB(ChNo) - DHSDEA(CHNO)*DEA ) /

     &            ( A11 - ALHX(CHNO)*DEDTC(ChNo) )

        EVAP(CHNO) = 0.

        IF(TC(CHNO)+DTC .GT. TF .AND. SNOW(CHNO).GT.0.) THEN

          Q0X=Q0-ALHM*SNOW(CHNO)/DTSTEP

          SMELT(CHNO)=SNOW(CHNO)/DTSTEP

          DTC = ( Q0X + ALHX(CHNO)*ETURB(ChNo) - DHSDEA(CHNO)*DEA ) /

     &              ( A11 - ALHX(CHNO)*DEDTC(ChNo) )

          IF(TC(CHNO)+DTC .LT. TF) THEN

            DTC=TF-TC(CHNO)

            SMELT(CHNO)=(Q0-Q0SNOW)/ALHM

      IF( CHOKE .AND. ABS(DEA) .GT. 0.5*EA(CHNO) ) THEN

        DEA = SIGN(.5*EA(CHNO),DEA)

        IF(TC(CHNO)+DTC .GT. TF .AND. SNOW(CHNO).GT.0.) THEN

          SNLEFT=SNOW(CHNO)-EVAP(CHNO)*DTSTEP*ESNFRC(CHNO)

          SMELT(CHNO)=SNLEFT/DTSTEP

          IF(TC(CHNO)+DTC .LT. TF) THEN

            DTC=TF-TC(CHNO)

            SMELT(CHNO)=(Q0-Q0SNOW)/ALHM

      HLWUP(CHNO) = HLWTC + DHLWTC*DTC

      SNOW(CHNO)=SNOW(CHNO)-SMELT(CHNO)*DTSTEP

      EA(CHNO) = MAX(EA(CHNO), 0.0 _d 0)

      VPDSTR(CHNO) = 1.

      FTEMP(CHNO) = MIN( 1. _d 0, MAX( FTEMP(ChNo), 1. _d -10 ) )

      HESAT = ESATTC(CHNO) * MIN( 1. _d 0, WET(CHNO)*2. _d 0)

      IF( EA(CHNO) .LT. HESAT ) THEN

          RSURF=RSURF*( 1. + (ESATTC(CHNO)-HESAT)/(HESAT-EA(CHNO)) )

      RX1(CHNO)=RC(CHNO)

      RX2(CHNO)=RSURF

      RC(ChNo) = RC(CHNO) * RSURF / ( RC(ChNo) + RSURF )

      IF(SATCAP(CHNO).GT..001) THEN

        RC(ChNo)=RC(ChNo)*(1.-POTFRC(CHNO))/

     &                ( 1.+POTFRC(CHNO)*RC(ChNo)/RA(ChNo) )

      IF(SATCAP(CHNO) .LE. .001) THEN

        RC(ChNo)=RC(ChNo)*(1.-SNWFRC(CHNO))/

     &                  ( 1.+SNWFRC(CHNO)*RC(ChNo)/RA(ChNo) )

      RAMPFC=ETURB(CHNO)/ETCRIT

      IF ( RAMPFC .GE. 0. ) RC(CHNO) = RC(CHNO)*(1.-RAMPFC)

      ESNO(CHNO)=ESNFRC(CHNO)*EVAP(CHNO)*DTSTEP

      EINT(CHNO)=EIRFRC(CHNO)*EVAP(CHNO)*DTSTEP

      EGRO = EVAP(CHNO)*DTSTEP - ESNO(CHNO) - EINT(CHNO)

      IF(ESNO(CHNO) .GT. SNOW(CHNO)) THEN

        EINT(CHNO)=EINT(CHNO)+(ESNO(CHNO)-SNOW(CHNO))

        ESNO(CHNO)=SNOW(CHNO)

      EGROI=EGRO+EINT(CHNO)

      IF(EGROI .GT. CAPAC(CHNO)+WTOTAL) THEN

        ESNO(CHNO)=ESNO(CHNO)+EGROI-(CAPAC(CHNO)+WTOTAL)

        EGROI=CAPAC(CHNO)+WTOTAL

      EINT(CHNO)=EGROI-EGRO

      IF(EINT(CHNO) .GT. CAPAC(CHNO)) THEN

        EGRO=EGRO+EINT(CHNO)-CAPAC(CHNO)

        EINT(CHNO)=CAPAC(CHNO)

        EINT(CHNO)=EINT(CHNO)+EGRO-WTOTAL

      IF( RX1(CHNO)+RX2(CHNO) .NE. 0. ) THEN

          ESOI(CHNO)=EGRO*RX1(CHNO)/(RX1(CHNO)+RX2(CHNO))

          EVEG(CHNO)=EGRO - ESOI(CHNO)

          ESOI(CHNO)=EGRO/2.

          EVEG(CHNO)=EGRO/2.

          FWS = EVEG(CHNO) / WTOTAL

      ELAY1 = WLAY1*FWS + ESOI(CHNO)

      IF(EVAP(CHNO) .LT. 0.) THEN

          EINT(CHNO)=(1.-ESNFRC(CHNO))*EVAP(CHNO)*DTSTEP

          ESOI(CHNO)=0.

          EVEG(CHNO)=0.

          ESNO(CHNO)=ESNFRC(CHNO)*EVAP(CHNO)*DTSTEP

      SNOW(ChNo) = SNOW(ChNo) - ESNO(CHNO)

      CAPAC(ChNo) = CAPAC(ChNo) - EINT(CHNO)

      SWET(ChNo,ROOTLY) = (WLAY2 - ELAY2) / VGWMAX(ROOTLY,ITYP(CHNO))

        FWSOIL(CHNO) = FWSOIL(CHNO) + (THRU-DEWRUN)/DTSTEP

      IF( SWET(CHNO,SFCLY) .GT. WETEQ1(CHNO) ) THEN

          FAREA=(SWET(CHNO,SFCLY) - WETEQ1(CHNO)) / (1. - WETEQ1(CHNO))

          DHDZ = 1. + 2.*(PHSAT(CHNO)-PHLAY(CHNO,ROOTLY))/ZDEP12

          AKAVE = AKSAT(CHNO)

      HALFMX=0.5*ABS( SWET(CHNO,SFCLY)-WETEQ1(CHNO) )

     &                               * WSMAX(SFCLY,ITYP(CHNO))

      FRAMP=(TC(CHNO)-TFM2)/2.

      HALFMX=0.5*ABS( SWET(CHNO,ROOTLY)-WETEQ2(CHNO) )

     &                               * WSMAX(ROOTLY,ITYP(CHNO))

      FRAMP=(TC(CHNO)-TFM2)/2.

      GDRAIN(CHNO)=GWFLUX/DTSTEP

      FRAMP=(TC(CHNO)-TFM2)/2.

      BOMB(CHNO) = 0.

      DEDQA(CHNO)  = MAX(  DEDQA(CHNO), 500. _d 0/ALHE )

      DEDTC(CHNO)  = MAX(  DEDTC(CHNO),   0. _d 0)

      DHSDQA(CHNO) = MAX( DHSDQA(CHNO),   0. _d 0)

      DHSDTC(CHNO) = MAX( DHSDTC(CHNO), -10. _d 0)

      EM(CHNO)     = QM(CHNO) * PSUR(CHNO) / EPSILON

      EA(CHNO)     = QA(CHNO) * PSUR(CHNO) / EPSILON

      ESATTC(CHNO) = QSATTC(CHNO) * PSUR(CHNO) / EPSILON

      DESDTC(CHNO) = DQSDTC(CHNO) * PSUR(CHNO) / EPSILON

      DEDEA(CHNO)  = DEDQA(CHNO) * EPSILON / PSUR(CHNO)

      DHSDEA(CHNO) = DHSDQA(CHNO) * EPSILON / PSUR(CHNO)

      PAR(CHNO)    = (PARDIR(CHNO) + PARDIF(CHNO) + 1.E-20)

      PDIR(CHNO)   = PARDIR(CHNO) / PAR(CHNO)

      RA(CHNO)     = ONE / ( CD(CHNO) * UM(CHNO) )

      CSOIL(CHNO)  = CSOIL0(ITYP(ChNo))

      CAPAC(CHNO)       = max( min(CAPAC(ChNo),SATCAP(CHNO)), 0. _d 0)

      SNWFRC(CHNO) = SNOW(CHNO) / ( SNOW(CHNO) + SNWMID(ITYP(CHNO)) )

      FCAN(CHNO) = MIN( 1. _d 0, MAX(0. _d 0,CAPAC(ChNo)/SATCAP(ChNo)) )

      POTFRC(CHNO)=1.-(1.-SNWFRC(CHNO))*(1.-FCAN(CHNO))

      EMAXRT(CHNO) = ( SNOW(CHNO) + CAPAC(CHNO) + WSOI12 ) / DTSTEP

      RUNOFF(CHNO) = 0.

      RUNSRF(CHNO) = 0.

      FWSOIL(CHNO) = 0.

      ESATTX(ChNo) = ESATTC(ChNo) + DESDTC(CHNO) * DELTC

      DO 200 CHNO=1,NCH

      RTBS=RX1(CHNO)*RX2(CHNO)/(RX1(CHNO)+RX2(CHNO))

      EPFRC=POTFRC(CHNO) * ( RA(CHNO) + RTBS ) /

     &                    ( RA(CHNO) + POTFRC(CHNO)*RTBS )

      ESNFRC(CHNO)=EPFRC*SNWFRC(CHNO)/

     &                      (SNWFRC(CHNO)+FCAN(CHNO)+1.E-20)

      EIRFRC(CHNO)=EPFRC*FCAN(CHNO)/(SNWFRC(CHNO)+FCAN(CHNO)+1.E-20)

      ALHX(CHNO) = (1.-ESNFRC(CHNO))*ALHE + ESNFRC(CHNO)*ALHS

      DO 300 CHNO=1,NCH

      IF(EVAP(CHNO) .GT. 0.) THEN

        EADJST=(ESNO(CHNO)/DTSTEP) - ESNFRC(CHNO)*EVAP(CHNO)

        SHFLUX(CHNO)=SHFLUX(CHNO)-EADJST*(ALHS-ALHE)

      QA(CHNO) = EA(CHNO) * EPSILON / PSUR(CHNO)

      HLATN (CHNO) = EVAP (CHNO) * ALHX(CHNO)

      IF(SWET1(ChNo).LT.1.E-10) SWET1(CHNO) = 0.0

      IF(SWET2(ChNo).LT.1.E-10) SWET2(CHNO) = 0.0

      IF(SWET3(ChNo).LT.1.E-10) SWET3(CHNO) = 0.0

      IF(CAPAC(ChNo).LT.1.E-10) CAPAC(CHNO) = 0.0

      IF(SNOW (ChNo).LT.1.E-10) SNOW (CHNO) = 0.0

      IF(RUNOFF(ChNo).LT.1.E-10) RUNOFF(CHNO) = 0.0

      EINT(CHNO) = EINT(CHNO) * ALHE / DTSTEP

      ESOI(CHNO) = ESOI(CHNO) * ALHE / DTSTEP

      EVEG(CHNO) = EVEG(CHNO) * ALHE / DTSTEP

      ESNO(CHNO) = ESNO(CHNO) * ALHS / DTSTEP

      DEFCIT = EVAP(CHNO) * ( RC(CHNO) + RA(CHNO) )

      IF( SNOW(CHNO).GT.0. .AND. TC(CHNO).GT.TF ) THEN

        IF( SNOWM .EQ. SNOW(CHNO) ) THEN

            SNOW(CHNO)=0.

            TC(CHNO)=TF

        SMELT(CHNO)=SMELT(CHNO)+SNOWM/DTSTEP

        XTCORR= (1.-TIMFRL) * MIN( 1. _d 0,(CAPAC(CHNO)/SATCAP(CHNO))/

      WATADD = TRAINL(ChNo)*DTSTEP + SMELT(CHNO)*DTSTEP

      CAVAIL = ( SATCAP(CHNO) - CAPAC(CHNO) ) * FWETL

      WETINT = CAPAC(CHNO)/SATCAP(CHNO)

      CAPAC(CHNO)=CAPAC(CHNO)+WATADD-THRU1-THRU2

      XTCORR= (1.-TIMFRC) * MIN( 1. _d 0,(CAPAC(CHNO)/SATCAP(CHNO))/

      CAVAIL = ( SATCAP(CHNO) - CAPAC(CHNO) ) * FWETC

      WETINT = CAPAC(CHNO)/SATCAP(CHNO)

      CAPAC(CHNO)=CAPAC(CHNO)+WATADD-THRU1-THRU2

      IF(SWET1(CHNO).GT.WETEQ1(CHNO) .AND. WETEQ1(CHNO).NE.1.) THEN

          WETFRC=(SWET1(CHNO)-WETEQ1(CHNO))/(1.-WETEQ1(CHNO))

      CAVAIL = ( 1.-WETFRC)*FWETL*WMAX(1,ITYP(ChNo))*(1-WETEQ1(CHNO))

          SWET1(CHNO) = SWET1(CHNO) + CAVAIL / WMAX(1,ITYP(ChNo))

      IF(SWET1(CHNO).GT.WETEQ1(CHNO) .AND. WETEQ1(CHNO).NE.1.) THEN

          WETFRC=(SWET1(CHNO)-WETEQ1(CHNO))/(1.-WETEQ1(CHNO))

      CAVAIL = (1.-WETFRC)*FWETC*WMAX(1,ITYP(ChNo))*(1-WETEQ1(CHNO))

          SWET1(CHNO) = SWET1(CHNO) + CAVAIL / WMAX(1,ITYP(ChNo))

      FWSOIL(CHNO) = FWSOIL(CHNO) + (THRUC+THRUL-WRUNC-WRUNL)/DTSTEP

      WET0  = MAX(WET(CHNO),0.01 _d 0)

      PHEQ = PHR(CHNO) - DELZ(ITYP(CHNO))

      WETEQ(CHNO) = ( PHEQ/VGPSAX(CHNO) ) ** ( -1/VGBEEX(CHNO) )

      HLWTC = ALWRAD(CHNO) + BLWRAD(CHNO) * TC(CHNO)

      DHLWTC = BLWRAD(CHNO)

     &        ALHX(CHNO)*DEDTC(ChNo) +

      A12 = DHSDEA(ChNo) + ALHX(CHNO) * DEDEA(ChNo)

     &        ALHX(CHNO) * ETURB(ChNo) -

          STRDG6(ChNo)=ALHX(CHNO)*ETURB(ChNo)/A11

      IF( RC(CHNO) .GT. 0.) THEN

              EHARMN = 2.*EPLANT*ETURB(CHNO) / (EPLANT + ETURB(ChNo))

      SMELT(CHNO)=0.

      SNLEFT=SNOW(CHNO)-EVAP(CHNO)*DTSTEP*ESNFRC(CHNO)

      IF(TC(CHNO)+DTC .GT. TF .AND. SNLEFT.GT.0.) THEN

        SMELT(CHNO)=SNLEFT/DTSTEP

        IF(TC(CHNO)+DTC .LT. TF) THEN

          DTC=TF-TC(CHNO)

          SMELT(CHNO)=(Q0-Q0SNOW)/ALHM