anchor.diff:          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(/,
  
positive.code:       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)
        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
  
negative.code:       end
  c
  c version
  c $Id: residu.f,v 1.4 2005/03/29 21:43:32 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  