c$Author: fredfox $
c$Date: 2001-11-21 23:05:24 $
c$Revision: 1.3 $
c$Source: /weru/cvs/weps/weps.src/mproc/loosn.for,v $

      subroutine loosn 
     i                 (u,tillf,nlay,density,sbd,laythk)


c     + + + PURPOSE + + +
     
c     This subroutine reads in the array(s) containing the components 
c     that need to be loosen/compact(ed). 

c     + + + KEYWORDS + + +
c     loosen/compact, tillage 

      include 'p1werm.inc'

c     + + + ARGUMENT DECLARATIONS + + +
      integer nlay
      real    u,tillf,density(mnsz),laythk(mnsz),sbd(mnsz)

c     + + + ARGUMENT DEFINITIONS + + +

c     nlay     - number of soil layers used
c     u        - loosening coefficient
c     tillf    - fraction of soil area tilled by the machine
c     density  - present soil bulk density
c     sbd      - settled soil bulk density
c     laythk   - layer thickness

c     + + + ACCESSED COMMON BLOCK VARIABLE DEFINITIONS + + +
c     mnsz        - max number of soil layers

c     + + + PARAMETERS + + +

c     + + + LOCAL VARIABLES + + +
      integer  i 
      real dum(mnsz)       

c     + + + LOCAL VARIABLE DEFINITIONS + + +

c     dum = dummy variable used in calculating the mass in a subregion
c     i = loop variable on layers 

c     + + + END SPECIFICATIONS + + + 

c     perform the loosen/compact process on the layers in a subregion 

      do 300 i=1,nlay
         dum(i)= density(i)-((density(i)-(2.0/3.0)*sbd(i))*u*tillf)
         laythk(i)=laythk(i)*(density(i)/dum(i))
         density(i)=dum(i)
300   continue
      end