subroutine soilinit(isr)
C ***************************************************************** wjr
C Contains init code from main
C
C       Edit History
C       04-Mar-99       wjr     created
C
      include 'p1werm.inc'
      include 's1layd.inc'
      include 's1layr.inc'
      include 's1dbc.inc'
C
      integer isr
      integer idx
C
c initialize soil depth to bottom of layers (mm) from layer thickness (mm)
      aszlyd(1, isr) = aszlyt(1, isr)
      do idx = 2, nslay(isr)
        aszlyd(idx,isr) = aszlyt(idx,isr) + aszlyd(idx-1, isr)
      end do
C
C This should go away, possibly becoming a data statement   
      asmno3(isr) = 0.

      end