anchor.diff:  c
   c...  subroutine arguments
   c
  -      integer nprop
  +      integer nprop,ierr
         integer iddmat(nstr,nstr)
  +      character errstrng*(*)
         double precision dmat(nddmat),prop(nprop)
   c
   c...  local variables
  
positive.code:       end
  c
  c
        subroutine td_mat_1(state,dmat,prop,iddmat,nstate,nprop,matchg)
  c
  c...  subroutine to form the material matrix for an integration point
  c     for the time-dependent solution.  This routine is meant to be
  
negative.code:       end
  c
  c
        subroutine elas_mat_1(dmat,prop,iddmat,nprop,ierr,errstrng)
  c
  c...  subroutine to form the material matrix for an integration point
  c     for the elastic solution.  The material matrix is assumed to be
  