anchor.diff:        end
   c
   c version
  -c $Id: eforce.f,v 1.2 2004/06/18 15:29:21 willic3 Exp $
  +c $Id: eforce.f,v 1.3 2004/06/21 20:11:19 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine eforce(xl,sh,shj,det,gauss,evp,p,iel,nen,ngauss,ierr,
       & getshape,bmatrix)
  c
  c...this subroutine computes the effective forces at each node
  c   within an element: p=(b)t*evp, where p is the local force vector,
  
negative.code: c...  subroutine arguments
  c
        integer iel,nen,ierr
        character errstrng*(*)
        double precision xl(nsd,nen),sh(nsd+1,nenmax,ngaussmax)
        double precision shj(nsd+1,nenmax,ngaussmax),det(ngaussmax)
        double precision gauss(nsd+1,ngaussmax),evp(nstr,ngauss)
  