anchor.diff:  c
         tmax=big
         call dcopy(nstr,ee,ione,dstate(1,2),ione)
  -      call dspmv("u",nstr,one,dmat,dstate(1,2),ione,zero,dstate(1,1),
  +      call dcopy(nstr,state0,ione,dstate(1,1),ione)
  +      call dspmv("u",nstr,one,dmat,dstate(1,2),ione,one,dstate(1,1),
        & ione)
   c
         return
         end
   c
   c
  -      subroutine td_strs_mat_1(state,dstate,ee,dmat,prop,rtimdat,rgiter,
  -     & ntimdat,iddmat,tmax,nstate,nprop,matchg,ierr,errstrng)
  +      subroutine td_strs_mat_1(state,dstate,state0,ee,dmat,prop,rtimdat,
  +     & rgiter,ntimdat,iddmat,tmax,nstate,nstate0,nprop,matchg,ierr,
  +     & errstrng)
   c
   c...  subroutine to compute the current stress and updated material
   c     matrix for the time-dependent solution.  Since this is a purely
  
positive.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
  
negative.code: c       
  
  c version
  c $Id: mat_1.f,v 1.15 2005/02/24 00:07:03 willic3 Exp $
  
  c Generated automatically by Fortran77Mill on Tue May 18 14:18:50 2004
  
  