anchor.diff:  c
         integer i,imat,matgpt,matmodel,nmatel,nprop,indprop
         logical matchg
  -      double precision gsum,dens,dif
  +      double precision gsum,dens
         double precision ptmp(100)
   c
   c...  included variable definitions
  
positive.code: c
  c
        subroutine gload_drv(
       & b,bres,gvec1,gvec2,grav,neq,                                     ! force
       & x,d,numnp,                                                       ! global
       & dx,numslp,                                                       ! slip
       & tfault,numfn,                                                    ! fault
  
negative.code:         gsum=gsum+grav(i)*grav(i)
        end do
        if(gsum.eq.zero) return
        call fill(bgravity,zero,neq)
        matgpt=1
  c
  c...  loop over material groups
  