anchor.diff:  c
   c...  subroutine arguments
   c
  -      integer n,nen,nsd,ngauss,ierr
  +      integer iel,nen,nsd,ngauss,ierr
         double precision xl(nsd,nen),sh(nsd+1,nenmax,ngaussmax)
         double precision shj(nsd+1,nenmax,ngaussmax)
         double precision shd(nsd+1,nenmax,ngaussmax),shbar(nsd+1,nenmax)
  
positive.code:       end
  c
  c version
  c $Id: getshapb.f,v 1.2 2004/06/15 19:46:30 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine getshapb(xl,sh,shj,shd,shbar,det,gauss,vol,iel,nen,nsd,
       & ngauss,ierr)
  c
  c...  Subroutine to compute shape functions and derivatives at
  c     Gauss points.
  