anchor.diff:  c
   c
         subroutine gload_drv(
  -     & b,bres,gvec1,gvec2,grav,neq,                                     ! force
  +     & bgravity,ngravflag,grav,neq,                                     ! force
        & x,d,numnp,                                                       ! global
        & dx,numslp,                                                       ! slip
        & tfault,numfn,                                                    ! fault
  
positive.code: c   positions (gvec2) and current body force (gvec1) and update
  c   body force vector and global load vector by this amount.
  c
        do i=1,neq
          dif=gvec2(i)-gvec1(i)
          gvec1(i)=gvec2(i)
          b(i)=b(i)+dif
          bres(i)=bres(i)+dif
        end do
  c
        return
        end
  c
  c version
  c $Id: gload_drv.f,v 1.2 2004/07/08 21:25:16 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code: c
        integer i,imat,matgpt,matmodel,nmatel,nprop,indprop
        logical matchg
        double precision gsum,dens
        double precision ptmp(100)
  c
  c...  included variable definitions
  