anchor.diff:        end
   c
   c
  -      subroutine td_mat_1(state,dmat,prop,iddmat,nstate,nprop,matchg)
  +      subroutine td_matinit_1(state,dstate,dmat,prop,iddmat,tmax,nstate,
  +     & nprop,matchg,ierr,errstrng)
   c
   c...  subroutine to form the material matrix for an integration point
   c     for the time-dependent solution.  This routine is meant to be
  
positive.code: c
  c...  subroutine arguments
  c
        integer nprop
        integer iddmat(nstr,nstr)
        double precision dmat(nddmat),prop(nprop)
  c
  c...  local variables
  
negative.code:       end
  c
  c
        subroutine elas_strs_1(state,ee,dmat,nstate,ierr,errstrng)
  c
  c...  subroutine to compute stresses for the elastic solution.  For this
  c     material, there are just 2 state variables:  total stress and
  