c
c
      subroutine lift
     i              (liftf,nlay,
     i               dflat,droot,dblwgnd,
     i               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    liftf
      real    dflat(mnbpls),dblwgnd(mnbpls,mnsz),
     &         droot(mnbpls,mnsz)
c
c
c     + + + ARGUMENT DEFINITIONS + + +
c
c     dblwgnd   - (decomp) below ground residue / layer and decomp
c                 pool (kg / m^2)
c     dflat     - (decomp) flat residue pools (kg / m^2)
c     droot     - (decomp) root mass / layer and decomp pool
c     liftf     - fraction of buried material lifted to the surface
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    liftlay(mnbpls,mnsz),lifttot(mnbpls)
c
c     + + + LOCAL VARIABLE DEFINITIONS + + +
c
c     i         - biomass pools (1-3)
c     lay       - number of layers in a specified subregion
c     liftlay   - buried material lifted to the surface in each layer
c     lifttot   - total buried material lifted to the surface
c
c     + + + END SPECIFICATIONS + + +
c
C *** debugging fix, may or may not be good
      lifttot(1) = 0.0
C *** eodf

c     perform the lifting of biomass

      do 101 i=2,mnbpls-1
       if (BTEST(bflg,i))then
          do 100 lay=1,nlay
            liftlay(i-1,lay)=dblwgnd(i-1,lay)*liftf
            lifttot(i-1)=lifttot(i-1)+droot(i-1,lay)*
     &       liftf+liftlay(i-1,lay)
100       continue
       endif
101   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) -
     &      liftlay(i-1,lay)
           droot(i-1,lay) = droot(i-1,lay)*(1.0-liftf)
200       continue
           dflat(i-1)=dflat(i-1)+lifttot(i-1)
         endif
201     continue
      end