c
c
      subroutine cbury
     i              (buryf,nlay,dflat,
     i               dblwgnd,bflg)


c     + + + PURPOSE + + +
c
c     This subroutine performs the biomass manipulation process of transfering
c     the above ground biomass into the soil or the inverse process of bringing
c     buried biomass to the surface.  It deals only with the biomass
c     pools (ie no live crop is involved)
c
c
c     + + + KEYWORDS + + +
c     bury, lift, biomass manipulation

*$noreference
      include 'p1werm.inc'
*$reference
c
c     + + + ARGUMENT DECLARATIONS + + +
      integer nlay, bflg
      real    buryf
      real    dflat(mnbpls),dblwgnd(mnbpls,mnsz)
c
c
c     + + + ARGUMENT DEFINITIONS + + +
c
c     buryf     - fraction of flat material buried
c     dblwgnd   - (decomp) below ground residue / layer and decomp
c                 pool (kg / m^2)
c     dflat     - (decomp) flat residue pools (kg / m^2)
c     nlay      - number of soil layers used in the operation(s)
c
c     + + + ACCESSED COMMON BLOCK VARIABLE DEFINITIONS + + +
c
c     mnbpls        - max number of biomass pools
c     mnsz          - max number of soil layers
c
c     + + + PARAMETERS + + +
c
c     + + + LOCAL VARIABLES + + +
c
      integer  lay,i
      real    bury(mnbpls)
c
c     + + + LOCAL VARIABLE DEFINITIONS + + +
c
c     bury      - mass of biomass that is buried
c     i         - biomass pools (1-3)
c     lay       - number of layers in a specified subregion
c
c     + + + END SPECIFICATIONS + + +
c
c     perform the burying of biomass

      do 100 i=2,mnbpls-1
       if (BTEST(bflg,i)) then
        bury(i-1)=dflat(i-1)*buryf
       endif
100   continue
c     Now let's update the 4 pool types using the temporary variables
c     we calculated above.
        do 201 i=2,mnbpls-1
         if (BTEST(bflg,i)) then
          do 200 lay=1,nlay
           dblwgnd(i-1,lay) = dblwgnd(i-1,lay)+bury(i-1)/
     &      nlay
200       continue
           dflat(i-1)=dflat(i-1)-bury(i-1)
         endif
201     continue

      end