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


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


*$noereference
      include 'p1werm.inc'
*$reference

c
c     + + + ARGUMENT DECLARATIONS + + +

      integer nlay
      real    u,tillf,density(mnsz),laythk(mnsz),sbd(mnsz)
c
c
c     + + + ARGUMENT DEFINITIONS + + +
c
c     nlay     - number of soil layers used
c     u        - loosening coefficient
c     tillf    - fraction of soil area tilled by the machine
c     density  - soil density
c     sbd      - soil bulk density
c     laythk   - layer thickness
c
c     + + + ACCESSED COMMON BLOCK VARIABLE DEFINITIONS + + +
c     
c     mnsz        - max number of soil layers
c
c     + + + PARAMETERS + + +
c
c     + + + LOCAL VARIABLES + + +
c
      integer  i 
      real dum(mnsz)       
c
c     + + + LOCAL VARIABLE DEFINITIONS + + +
c
c     dum = dummy variable used in calculating the mass in a subregion
c     i = loop variable on layers 
c
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