anchor.diff:        integer neq,nextflag,numnp,iter,itmaxp,idebug,idout,kto,kw
         integer id(ndof,numnp),idx(ndof,numnp)
         double precision bextern(nextflag*neq),bintern(neq),bresid(neq)
  -      double precision dispvec(neq)
  +      double precision dispvec(neq),gtol(3),gi(3),gprev(3),gcurr(3)
         logical converge
   c
   c...  included dimension and type statements
   c
  -      include "gtol_dim.inc"
  -      include "gi_dim.inc"
  -      include "gprev_dim.inc"
  -      include "gcurr_dim.inc"
   c
   c...  intrinsic functions
   c
  
positive.code:         write(kto,900) gcurr(1),gcurr(2),gcurr(3),gi(1),gi(2),gi(3),
       &   acc(1),acc(2),acc(3)
        end if
  cdebug      write(6,*) "From end of residu_f!"
  cdebug      write(6,*) "bextern:",(bextern(idb),idb=1,neq)
  cdebug      write(6,*) "bintern:",(bintern(idb),idb=1,neq)
  cdebug      write(6,*) "bresid:",(bresid(idb),idb=1,neq)
  cdebug      write(6,*) "dispvec:",(dispvec(idb),idb=1,neq)
  800   format(/,
       & '   WARNING!  Apparent divergence in iteration #',i5,'!',//)
  810   format(/,
  
negative.code:       logical debug,div(3)
  c
  cdebug      write(6,*) "Hello from residu_f!"
  c
        tmp(3)=zero
        debug=(idebug.eq.1).and.(idout.gt.1)
  