anchor.diff:  c...  subroutine arguments
   c
         integer nen,iel,ierr
  -      character errstrng
  +      character errstrng*(*)
         double precision x(nsd,nen),xs(nsd,nsd),det,shj(nsd+1,nenmax)
   c
   c...  local variables
  
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: 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
  