real function timepk
     i         (timpkd,k10)
c
c     + + + PURPOSE + + +
c     Calculates the time to peak from the 12 interval time to peak
c     accummulated distribution parameters TIMPKD input for each station
c     location.
c
c     + + + ARGUMENT DECLARATIONS + + +
      real timpkd(0:12)
      integer k10(4)
c
c     + + + ARGUMENT DEFINITIONS + + +
c     timpkd      - The 12 interval time to peak accummulated distribution
c                   parameters for the station.  Cumulative distribution of 
c                   computed time to peak rainfall intensity values based on
c                   NWS 15-minute rainfall data (section 2.1.4 WEPP tech 1995)
c     k10         - Seed for random number generation.
c
      include 'crandom.inc'
c        read: zx
c     + + + LOCAL VARIABLES + + +
      integer i
      real    diff1, diff2, ratio
c
c     + + + END SPECIFICATIONS + + +
c
      z=zx(dax)
      i=0
 10   continue
        i=i+1
      if(timpkd(i).lt.z .and. i.lt.12) goto 10
c
      diff1= timpkd(i)-z
      diff2=timpkd(i) - timpkd(i-1)
      ratio= diff1/diff2
      timepk=0.08333*i - ratio*0.08333
c
      return
      end