anchor.diff:        end
   c
   c
  -      subroutine elas_strs_1(state,ee,dmat,nstate)
  +      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
  
positive.code: c
  c...  subroutine arguments
  c
        integer nstate
        double precision state(nstr,nstate),ee(nstr),dmat(nddmat)
  c
        call dcopy(nstr,ee,ione,state(1,2),ione)
  
negative.code:       end
  c
  c
        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
  