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: c
  c...  subroutine arguments
  c
        integer ngravflag,neq,numnp,numslp,numfn,numelt,nconsz,numat
        integer npropsz,nhist,lastep,numrot,ierr
        integer ien(nconsz),lm(ndof,nconsz),lmx(ndof,nconsz),lmf(nconsz)
        integer infiel(6,numelt),mhist(npropsz),infmat(3,numat)
        integer infmatmod(5,nmatmodmax),infetype(4,netypes)
        character errstrng*(*)
        double precision bgravity(ngravflag*neq),grav(ndof)
        double precision x(nsd,numnp),d(ndof,numnp),dx(ndof,numnp)
        double precision tfault(ndof,numfn),prop(npropsz)
        double precision gauss(nsd+1,ngaussmax,netypes)
  