!$Author: joelevin $
!$Date: 2011-03-24 11:33:26 -0500 (Thu, 24 Mar 2011) $
!$Revision: 11724 $
!$HeadURL: https://svn.weru.ksu.edu/weru/weps1/trunk/weps.src/src/lib_soil/depthini.for $

      subroutine depthini(nlay, bszlyt, bszlyd)

      integer nlay
      real    bszlyt(*), bszlyd(*)

      integer idx

!     nlay - number of soil layers
!     bszlyt - soil layer thickness (mm)
!     bszlyd - depth to bottom of soil layer (mm)

      bszlyd(1) = bszlyt(1)
      do idx = 2, nlay
        bszlyd(idx) = bszlyt(idx) + bszlyd(idx-1)
      end do

      return
      end