anchor.diff:  c
   c...  subroutine arguments
   c
  -      integer nprop,ierr
  +      integer ipauto,nprop,ierr
         integer iddmat(nstr,nstr)
         character errstrng*(*)
  -      double precision autofac,dmat(nddmat),prop(nprop)
  +      double precision tpois,tyoungs,dmat(nddmat),prop(nprop)
   c
   c...  local variables
   c
         double precision ptmp(10)
   c
         call dcopy(nprop,prop,ione,ptmp,ione)
  -      ptmp(2)=autofac*ptmp(2)
  +      if(ipauto.eq.ione) then
  +        ptmp(2)=tyoungs
  +        ptmp(3)=tpois
  +      end if
         call elas_mat_5(dmat,ptmp,iddmat,nprop,ierr,errstrng)
         return
         end
   c       
   
   c version
  -c $Id: mat_5.f,v 1.6 2005/01/18 19:40:36 willic3 Exp $
  +c $Id: mat_5.f,v 1.7 2005/02/24 00:14:45 willic3 Exp $
   
   c Generated automatically by Fortran77Mill on Tue May 18 14:18:50 2004
   
  
positive.code: c
  c...  subroutine arguments
  c
        integer nstate,nprop,ierr
        integer iddmat(nstr,nstr)
        character errstrng*(*)
        logical matchg
        double precision state(nstr,nstate),dstate(nstr,nstate),ee(nstr)
        double precision dmat(nddmat),prop(nprop),tmax
  c
  c...  included dimension and type statements
  
negative.code:       pr=prop(3)
        vis=prop(4)
        rmu=half*e/(one+pr)
        rvis=one/vis
        tmax=two*vis/rmu
        fac1=half*deltp*rvis
        fac2=(one+pr)/e
        f1=fac1*(one-alfap)
        f2=one/(fac2+fac1*alfap)
        emean=third*(ee(1)+ee(2)+ee(3))
        smean=e*emean/(one-two*pr)
        smean0=third*(state0(1)+state0(2)+state0(3))
        smeantp=third*(state(1,1)+state(2,1)+state(3,1))
  cdebug      write(6,*) "e,pr,vis,rmu,tmax,deltp,alfap,f1,f2,emean,smean"
  cdebug      write(6,*) "smean0,smeantp,emean:"
  cdebug      write(6,*) e,pr,vis,rmu,tmax,deltp,alfap,f1,f2,emean,smean,smean0,
  cdebug     & smeantp,emean
        do i=1,nstr
          eet=ee(i)-diag(i)*emean-state(i,3)
          sdevtp=state(i,1)-diag(i)*smeantp
          sdev0=state0(i)-diag(i)*smean0
          sdev=f2*(eet-f1*sdevtp+fac2*sdev0)
          dstate(i,1)=sdev+diag(i)*(smean+smean0)
          stau=(one-alfap)*sdevtp+alfap*sdev
          dstate(i,3)=half*deltp*stau*rvis
          dstate(i,2)=ee(i)
  cdebug        write(6,*) "i,ee,eet,sdevtp,sdev0,sdev,ds1,ds2,ds3,s1,s2,s3,",
  cdebug     &   "stau:",
  cdebug     &   i,ee(i),eet,sdevtp,sdev0,sdev,dstate(i,1),dstate(i,2),
  cdebug     &   dstate(i,3),state(i,1),state(i,2),state(i,3),stau
        end do
  c
        return
        end
  c
  c
        subroutine td_strs_mat_5(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.
  