anchor.diff:        dmat(iddmat(1,2))=f1-f2
         dmat(iddmat(1,3))=dmat(iddmat(1,2))
         dmat(iddmat(2,3))=dmat(iddmat(1,2))
  -      dmat(iddmat(4,4))=three*f2
  +      dmat(iddmat(4,4))=half*three*f2
         dmat(iddmat(5,5))=dmat(iddmat(4,4))
         dmat(iddmat(6,6))=dmat(iddmat(4,4))
  -cdebug      write(6,*) "dmat:",(dmat(idb),idb=1,nddmat)
  -      return
  -      end
  -c
  -c
  -      subroutine td_strs_5(state,dstate,state0,ee,dmat,prop,rtimdat,
  -     & rgiter,ntimdat,iddmat,tmax,nstate,nstate0,nprop,matchg,ierr,
  -     & errstrng)
  +      return
  +      end
  +c
  +c
  +      subroutine td_strs_5(state,dstate,state0,ee,scur,dmat,prop,
  +     & rtimdat,rgiter,ntimdat,iddmat,tmax,nstate,nstate0,nprop,matchg,
  +     & ierr,errstrng)
   c
   c...  subroutine to compute the current values for stress, total strain,
   c     and viscous strain increment.
  
positive.code:       integer iddmat(nstr,nstr)
        character errstrng*(*)
        logical matchg
        double precision state(nstr,nstate),dstate(nstr,nstate)
        double precision state0(nstate0),ee(nstr)
        double precision dmat(nddmat),prop(nprop),tmax
  c
  c...  included dimension and type statements
  c
  
negative.code:           dmat(iddmat(i,j))=od
          end do
        end do
        return
        end
  c
  c
        subroutine elas_strs_5(state,state0,ee,scur,dmat,nstate,nstate0,
       & ierr,errstrng)
  c
  c...  subroutine to compute stresses for the elastic solution.  For this
  c     material, there are 3 state variables:  total stress, total
  c     strain, and viscous strain.  The current total strain is contained
  c     in ee.
  c
  c     state(1:6)   = Cauchy stress
  c     state(7:12)  = linear strain
  c     state(13:18) = viscous strain
  c
  c     The state0 array contains initial stresses.
  c
        include "implicit.inc"
  c
  