anchor.diff:  c
         tmp(3)=zero
         debug=(idebug.eq.1).and.(idout.gt.1)
  -      if(debug) call printv(b,btot,id,idx,neq,ndof,numnp,iout,idout,kw)
  +      if(debug) call printv(b,btot,id,idx,neq,numnp,iout,idout,kw)
         ftmp=dnrm2(neq,bres,ione)
         do i=1,neq
           en=deld(i)*bres(i)
  
positive.code:       end
  c
  c version
  c $Id: residu.f,v 1.1 2004/04/14 21:18:30 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code: 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
        include "implicit.inc"
  c
  c...  parameter definitions
  c
        include "ndimens.inc"
        include "nconsts.inc"
        include "rconsts.inc"
        integer iout
        parameter(iout=1)
  c
  c...  subroutine arguments
  c
        integer neq,numnp,iter,itmaxp,idebug,idout,kto,kw
        integer id(ndof,numnp),idx(ndof,numnp)
        double precision b(neq),bres(neq),btot(neq),deld(neq)
        logical converge
  