anchor.diff:  c
         vol=zero
         do l=1,ngauss
  -        call getjac(xl,xs,det(l),shj(1,1,l),nen,nsd,iel,ierr)
  -        if(ierr.ne.0) return
  -        call getder(det(l),sh(1,1,l),shd(1,1,l),xs,nen,nsd)
  +        call getjac(xl,xs,det(l),shj(1,1,l),nen,iel,ierr,errstrng)
  +        if(ierr.ne.izero) return
  +        call getder(det(l),sh(1,1,l),shd(1,1,l),xs,nen)
           det(l)=gauss(4,l)*det(l)
           vol=vol+det(l)
         end do
  
positive.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine getshapn(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.
  
negative.code:       end
  c
  c version
  c $Id: getshapn.f,v 1.4 2004/07/05 20:04:32 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  