anchor.diff:        iel=infiel(4,matgpt)
         inddmat0=infiel(6,iel)
         ietype=infiel(3,iel)
  -      call elas_mat(dmat(1,inddmat0),prop,iddmat,nprop)
  +      call elas_mat(dmat(1,inddmat0),prop,iddmat,nprop,ierr,errstrng)
  +      if(ierr.ne.izero) return
         ngauss=infetype(1,ietype)
         ng=ngauss
         if(imatvar.eq.izero) ng=ngaussmax
  
positive.code:       subroutine elas_mat_cmp(
       & dmat,infiel,iddmat,ndmatsz,numelt,                               ! elemnt
       & prop,infmat,nprop,matgpt,elas_mat,                               ! materl
       & infetype)                                                        ! eltype
  c
  c...  compute subroutine to form the d-matrix for the elastic solution
  c
  
negative.code:       end
  c
  c version
  c $Id: elas_mat_cmp.f,v 1.3 2004/06/21 20:14:47 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  