!
!$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_mproc/resinit.for $
!
! This routine takes a residue amount and depth and distributes it
! into the array of soil layers

      subroutine resinit(resmass, resdepth, nlay, resarray, laythick)

!     + + + INPUT VARIABLE DECLARATIONS + + +
      real resmass
      real resdepth
      integer nlay
      real resarray(nlay)
      real laythick(nlay)

!     + + + INPUT VARIABLE DEFINITIONS + + +
!     resmass - residue mass (Kg/m^2)
!     resdepth - Depth residue is distributed in soil (mm)
!     nlay - number of soil layers
!     resarray(nlay) - soil residue array by layer (Kg/m^2)
!     laythick(nlay) - soil layer thickness (mm)

!     + + + LOCAL VARIABLE DECLARATIONS + + +
      integer ilay
      real    depth

!     + + + LOCAL VARIABLE DEFINITIONS + + +
!     ilay - array index
!     depth - accumulator for depth
!     thick - thickness of slice to which residue is to be added

      depth = resdepth
      do ilay = 1, nlay
          if (depth.gt.0.0) then
              resarray(ilay) = resmass                                  &
     &                       * min( depth,laythick(ilay)) / (resdepth)
              depth = depth - laythick(ilay)
          else
              resarray(ilay) = 0.0
          end if
      end do

      return
      end