anchor.diff:        include "rconsts.inc"
         integer iout
         parameter(iout=1)
  +      double precision alpha
  +      parameter(alpha=-1.0d0)
   c
   c...  subroutine arguments
   c
  -      integer neq,numnp,iter,itmaxp,idebug,idout,kto,kw
  +      integer neq,nextflag,numnp,iter,itmaxp,idebug,idout,kto,kw
         integer id(ndof,numnp),idx(ndof,numnp)
  -      double precision b(neq),bres(neq),btot(neq),deld(neq)
  +      double precision bextern(nextflag*neq),bintern(neq),bresid(neq)
  +      double precision dispvec(neq)
         logical converge
   c
   c...  included dimension and type statements
  
positive.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine residu(b,bres,btot,deld,gtol,gi,gprev,gcurr,id,idx,neq,
       & numnp,iter,itmaxp,idebug,idout,kto,kw,converge)
  c
  c...computes the deviation from overall force and energy balance
  c
  
negative.code:       end
  c
  c version
  c $Id: residu.f,v 1.3 2005/01/05 22:33:11 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  