!$Author: fredfox $
!$Date: 2006-02-25 00:21:16 $
!$Revision: 1.2 $
!$Source: /weru/cvs/wepp/wepp.watbal/parest.for,v $

      SUBROUTINE PAREST(DEPSTO, KS, SM, POR, SAT, CC, SC)

C     + + + PURPOSE + + +
C
C     THIS SUBROUTINE GETS KS AND SUCTION EITHER FROM USER OR 
C     ESTIMATES THE VALUES FROM SOIL TEXTURE AND CANOPY AND 
C     GROUND COVER.
C
C     CALLED FROM MAIN
C     AUTHOR(S): D. FLANAGAN, J. ASCOUGH
C     VERSION: THIS MODULE TAKEN FROM ASCOUGH STANDALONE IRS CODE
C     DATE CODED:  3-22-2005
C     CODED BY: D. FLANAGAN
C
C     + + + ARGUMENT DECLARATIONS + + +

      REAL DEPSTO, KS, SM, POR, SAT, CC, SC

C     + + + ARGUMENT DEFINITIONS + + +
C
C     DEPSTO  -
C     KS      -
C     SM      -
C     KSM     -
C     POR     -
C     SAT     -
C     CC      -
C     SC      -
C
C     + + + LOCAL VARIABLES + + +
C
C     INTEGER I, J, L, M, N, POND, IU, NTP, IC
      INTEGER IC , KSFLG, KCOV, KDEP
      REAL SF
C     REAL X, Y, Z
C     INTEGER LB, UB
C     CHARACTER*1 ANS
C
C     + + + END SPECIFICATIONS + + +
C
      REAL K1(12), K2(12), K3(12)
      DATA K1 /118., 30., 11., 6.5, 3.4, 2.5, 1.5, 1., .9, .6, .5, .4/
      DATA K2 /49., 63., 90., 110., 173., 190., 214., 210., 253., 260.,
     1    288., 310./
      DATA K3 /.4, .4, .41, .43, .49, .42, .35, .31, .43, .32, .42, .39/
C
C-----------------------------------------------------------------------
C
C     WRITE (*,1100)
C     
C     WRITE (*,1200)
C  
C     READ FLAG FOR ESTIMATING KS,SF, AND POROSITY BY SOIL TEXTURE
C     (KSFLG = 1 FOR YES; KSFLG = 0 FOR NO)       
C
      READ (3,*) KSFLG
C     READ (*,1300) ANS
C     IF (ANS.EQ.'Y'.OR.ANS.EQ.'Y') THEN
      IF (KSFLG.EQ.1) THEN
C        
C        
C       WRITE (*,1400)
C        
C       WRITE (*,1500)
C       LB = 1
C       UB = 12
C       CALL GETI(IC,LB,UB,IEND)
C
C       READ THE SOIL TEXTURE CLASS CODE (1-12)
C
C                       KS      SF    EFF.
C     TEXTURE        (MM/HR)   (MM)   POR.  CODE
C     __________________________________________
C     SAND             118      49    0.40   1
C     LOAMY SAND       30       63    0.40   2
C     SANDY LOAM       11       90    0.41   3
C     LOAM             6.5      110   0.43   4
C     SILT LOAM        3.4      173   0.49   5
C     SILT             2.5      190   0.42   6
C     S. CLAY LOAM     1.5      214   0.35   7
C     CLAY LOAM        1.0      210   0.31   8
C     SL. CLAY LOAM    0.9      253   0.43   9
C     SANDY CLAY       0.6      260   0.32   10
C     SILTY CLAY       0.5      288   0.42   11
C     CLAY             0.4      310   0.39   12
C     __________________________________________

        READ (3,*) IC
C        
C       KS (MM/HR) AND SUCTION (MM) ESTIMATED 
C       FROM SOIL TEXTURE
C        
        KS = K1(IC)
        SF = K2(IC)
        POR = K3(IC)
C       WRITE (*,1600) KS
C        
C       GET SUCTION
C        
C       WRITE (*,1700) SM
C     
      ELSE
C        
C       USER INPUT OF KS AND SUCTION
C        
C       WRITE (*,2100)
C       CALL GETR(XKKS,.0,1000.,IEND)
C       KS = XKKS
C       WRITE (*,2200)
C       CALL GETR(SSM,0.,10000.,IEND)
C       SM = SSM
C       WRITE (*,2300)
C       CALL GETR(PPOR,0.,100.,IEND)
C
C       USER INPUT OF SATURATED CONDUCTIVITY (MM/HR), AVE.
C       SUCTION ACROSS THE WETTING FRONT (MM), AND EFFECTIVE
C       POROSITY (%)
C
        READ (3,*) KS, SF, POR
        POR = POR / 100.0
      END IF
C     
C     
C  10 WRITE (*,2400)
C     CALL GETR(SSAT,0.,100.,IEND)
C     READ THE INITIAL SATURATION (%)
      READ (3,*) SAT       
C     SAT = SSAT
C     IF (SAT.LT.1.) THEN
C       WRITE (*,2500)
C       READ (*,1300) ANS
C       IF (ANS.EQ.'N'.OR.ANS.EQ.'N') GO TO 10
C     END IF
      SM = (1.-SAT/100.) * POR * SF
C     
C     GET CANOPY AND COVER
C     
C     WRITE (*,1000)
C     READ (*,1300) ANS
C     READ FLAG FOR ADJUSTING KS BASED ON CANOPY AND
C     GROUND COVER
C 
      READ (3,*) KCOV
C     IF (ANS.EQ.'Y'.OR.ANS.EQ.'Y') THEN
      IF (KCOV.EQ.1) THEN
C       WRITE (*,1800)
C       WRITE (*,1900)
C       CALL GETR(CCC,0.,100.,IEND)
C       CC = CCC
C       READ THE TOTAL CANOPY COVER AND GROUND COVER
C
        READ (3,*) CC, SC
C        
C       WRITE (*,2000)
C       CALL GETR(SSC,0.,100.,IEND)
C       SC = SSC
C       
        KS = KS * EXP(8.9999999E-03*SC) * EXP(.0105*CC)
      END IF
C     
C     READ FLAG FOR ADJUSTING RAINFALL EXCESS BASED ON
C     DEPRESSION STORAGE
C     
      READ (3,*) KDEP
C      
      IF (KDEP.EQ.1) THEN
C     
C       READ THE DEPRESSION STORAGE (MM)       
C
        READ (3,*) DEPSTO
      ELSE
C
C       DEFAULT VALUE OF DEPRESSION STORAGE IS 0.0
C      
        DEPSTO = 0.0
      END IF
C     
C     WRITE(*,200)
C200  FORMAT(//' DO YOU WANT RAINFALL EXCESS ADJUSTED FOR DEPRESSION',
C     1          ' STORAGE? (Y/N) ' $)
C     READ(*,40)ANS
C     
C     IF (ANS.EQ.'Y'.OR.ANS.EQ.'Y') THEN
C       WRITE(*,210)
C210    FORMAT(//' ENTER DEPRESSION STORAGE (MM) ---> ' $)
C       CALL GETR(DEPSTO,0.,1000.,IEND)
C     END IF
C     
      RETURN
      END