U                     CAPAC, SNOW, SWET, RUNOFF,

     &       CAPAC(NCH),    SNOW(NCH),    SWET(nch,NLAY),

      WLAY1 = SWET(ChNo,SFCLY) * VGWMAX(SFCLY,ITYP(ChNo))

      WLAY2 = SWET(ChNo,ROOTLY) * VGWMAX(ROOTLY,ITYP(ChNo))

      SWET(ChNo,SFCLY) = (WLAY1 - ELAY1) / VGWMAX(SFCLY,ITYP(ChNo))

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

      SWET(ChNo,SFCLY)  = MIN( 1. _d 0, MAX( 0. _d 0, SWET(ChNo,SFCLY) )

      SWET(ChNo,ROOTLY) = MIN( 1. _d 0, MAX( 0. _d 0, SWET(ChNo,ROOTLY)

        SWET(ChNo,SFCLY) = SWET(ChNo,SFCLY) +

        IF ( SWET(ChNo,SFCLY) .GT. 1. ) THEN

          DEWRUN = ( SWET(ChNo,SFCLY) - 1. ) * VGWMAX(SFCLY,ITYP(ChNo))

          SWET(ChNo,SFCLY) = 1.

     U                     SWET,   RUNOFF,  GDRAIN

     &                  DTSTEP,    SWET(nch,NLAY),

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

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

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

     &               SWET(ChNo,SFCLY) * WSMAX(SFCLY,ITYP(ChNo)) )

     &              SWET(ChNo,ROOTLY) * WSMAX(ROOTLY,ITYP(ChNo)) )

     &                SWET(ChNo,ROOTLY) * WSMAX(ROOTLY,ITYP(ChNo)) )

     &                SWET(ChNo,SFCLY) * WSMAX(SFCLY,ITYP(ChNo)) )

      SWET(ChNo,SFCLY) = SWET(ChNo,SFCLY) -

      SWET(ChNo,ROOTLY) = SWET(ChNo,ROOTLY) +

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

     &                    SWET(ChNo,ROOTLY) * WSMAX(ROOTLY,ITYP(ChNo)) )

     &                    SWET(ChNo,RECHLY) * WSMAX(RECHLY,ITYP(ChNo)) )

     &                    SWET(ChNo,RECHLY) * WSMAX(RECHLY,ITYP(ChNo)) )

     &                    SWET(ChNo,ROOTLY) * WSMAX(ROOTLY,ITYP(ChNo)) )

      SWET(ChNo,ROOTLY) = SWET(ChNo,ROOTLY) -

      SWET(ChNo,RECHLY) = SWET(ChNo,RECHLY) +

     &                 SWET(ChNo,RECHLY) * WSMAX(RECHLY,ITYP(ChNo)) )

      SWET(ChNo,RECHLY) = SWET(ChNo,RECHLY) -

      SWET(ChNo,SFCLY ) = max( min(SWET1(ChNo),1. _d 0), 0. _d 0)

      SWET(ChNo,ROOTLY) = max( min(SWET2(ChNo),1. _d 0), 0. _d 0)

      SWET(ChNo,RECHLY) = max( min(SWET3(ChNo),1. _d 0), 0. _d 0)

      WSOI12=SWET(ChNo,SFCLY ) * VGWMAX(SFCLY ,ITYP(ChNo)) +

     &       SWET(ChNo,ROOTLY) * VGWMAX(ROOTLY,ITYP(ChNo))

     I             NCH, UM, U2FAC, Z2, RDC, SWET(1,SFCLY),

      CALL RSURFP (NCH, UM, U2FAC, Z2, RDC, SWET(1,SFCLY),

      CALL RSURFP (NCH, UM, U2FAC, Z2, RDC, SWET(1,SFCLY),

     U                 TC,     TD,     EA,   SWET,   SNOW,

     U               CAPAC, SNOW, SWET, RUNOFF,

     I             NCH,          ITYP,        SWET (1, SFCLY), VGBEEX,

     I             NCH,          ITYP,        SWET(1,ROOTLY), VGBEEX,

     I             NCH,          ITYP,        SWET(1,RECHLY), VGBEEX,

     U             SWET, RUNOFF, GDRAIN

     I             NCH,          ITYP,        SWET(1,ROOTLY), VGBEEX,

     U             TC, CAPAC, SNOW, SWET, RUNOFF, RUNSRF,

      SWET1(ChNo) = SWET(ChNo,SFCLY)

      SWET2(ChNo) = SWET(ChNo,ROOTLY)

      SWET3(ChNo) = SWET(ChNo,RECHLY)

      _RL     SWET(nch,NLAY),   VGPSAT(NTYPS),   VGCSAT (NTYPS),