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:       end
  c
  c version
  c $Id: getjac.f,v 1.6 2004/08/12 01:28:23 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.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)
  cdebug      do idb=1,nen
  cdebug        write(6,*) idb,(x(jdb,idb),jdb=1,nsd)
  cdebug      end do
        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
  