subroutine updres(isr)
C *****************************************************************

C     Update geometric properties of the decomp residue pools

      include 'p1werm.inc'
      include 'p1const.inc'
      include 's1layr.inc'
      include 'c1db1.inc'
      include 'd1glob.inc'
      include 'c1glob.inc'
      include 'decomp/decomp.inc'
C
C arguments
C
      integer isr
C
C local variables
C
      integer idx,jdx
C      
c     + + + FUNCTIONS CALLED + + +

c *****************************************************************

c     For each decomp pool compute
c       residue SAI 
c       residue LAI (currently zero)
c       flat residue surface cover
c       standing residue surface cover
c       total surface residue cover

      do 10 idx=1,mnbpls
         adrsai(idx,isr) = adzht(idx,isr) * addstm(idx,isr) *
     &                     adxstm(idx,isr)
         adrlai(idx,isr) = 0.0
         adffcv(idx,isr) = 1.0 - exp(-covfact(idx,isr)*admf(idx,isr))
         adfscv(idx,isr) = addstm(idx,isr) * pi *
     &                     (adxstm(idx,isr)/2.0)*(adxstm(idx,isr)/2.0)
         adftcv(idx,isr) = (1.0 - adffcv(idx,isr)) * adfscv(idx,isr) +
     &                     adffcv(idx,isr) 
10    continue

c     Update SAI and LAI by height in the future
c     Currently SAI by height "adrsaz" is not being used
c     Currently LAI by height "adrlaz" is not being used

    ! print *, 'before summation: ', 
    !&         acmrt(isr), admrt(1,isr), admrt(2,isr), admrt(3,isr)
    !&         admbg(1,isr), admbg(2,isr), admbg(3,isr)

      !update crop root mass here for now - LEW
      acmrt(isr) = 0.0
      do 15  jdx=1,nslay(isr)
         acmrt(isr) = acmrt(isr) + acmrtz(jdx,isr)
15    continue

c     Root biomass and below ground biomass distributed by layer
      do 30 idx=1,mnbpls
         admbg(idx,isr) = 0.0
         admrt(idx,isr) = 0.0
         do 20 jdx=1,nslay(isr)
            admbg(idx,isr) = admbg(idx,isr) + admbgz(jdx,idx,isr)
            admrt(idx,isr) = admrt(idx,isr) + admrtz(jdx,idx,isr)
20       continue
30    continue
    ! print *, 'after summation: ', 
    !&         acmrt(isr), admrt(1,isr), admrt(2,isr), admrt(3,isr)
    !&         admbg(1,isr), admbg(2,isr), admbg(3,isr)

      do 40 idx=1,mnbpls
         adm(idx,isr) = admbg(idx,isr) + admrt(idx,isr) +
     &              admst(idx,isr) + admf(idx,isr)
40    continue
      end