anchor.diff:  cdebug      call flush(6)
         if(det.le.zero) then
           ierr=113
  +cdebug        write(6,700) iel,det
           write(errstrng,700) iel,det
         end if
   c
  
positive.code: c...calculate jacobian matrix for (x,y,z) to (r,s,t) transformation
  c
  cdebug      call dgemm("n","t",nsd,nsd,nen,one,shj,nsd+1,x,nsd,zero,xs,nsd)
        call dgemm("n","t",nsd,nsd,nen,one,x,nsd,shj,nsd+ione,zero,xs,nsd)
  c
  c...form determinant of jacobian matrix and check for error condition
  
negative.code:       end
  c
  c version
  c $Id: getjac.f,v 1.7 2005/02/24 00:00:13 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  