anchor.diff:        end
   c
   c version
  -c $Id: getshapb.f,v 1.2 2004/06/15 19:46:30 willic3 Exp $
  +c $Id: getshapb.f,v 1.3 2004/06/17 18:09:43 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code: c
        vol=zero
        do l=1,ngauss
          call getjac(xl,xs,det(l),shj(1,1,l),nen,nsd,n,ierr)
          if(ierr.ne.0) return
          call getder(det(l),sh(1,1,l),shd(1,1,l),xs,nen,nsd)
          det(l)=gauss(4,l)*det(l)
  
negative.code: c
  c...  subroutine arguments
  c
        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)
  