anchor.diff:        pr=prop(3)
         vis=prop(4)
         f1=third*e/(one-two*pr)
  -      f2=one/((one+pr)/e+half*alfap*deltp/vis)
  -      dmat(iddmat(1,1))=f1+two*third*f2
  +      f2=third/((one+pr)/e+half*alfap*deltp/vis)
  +      dmat(iddmat(1,1))=f1+two*f2
         dmat(iddmat(2,2))=dmat(iddmat(1,1))
         dmat(iddmat(3,3))=dmat(iddmat(1,1))
  -      dmat(iddmat(1,2))=f1-third*f2
  +      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))=half*f2
  +      dmat(iddmat(4,4))=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)
  
positive.code:       end
  c
  c
        subroutine elas_strs_5(state,ee,dmat,nstate,ierr,errstrng)
  c
  c...  subroutine to compute stresses for the elastic solution.  For this
  c     material, there are 3 state variables:  total stress, total
  
negative.code: c
  c...  subroutine arguments
  c
        integer ipauto,nprop,ierr
        integer iddmat(nstr,nstr)
        character errstrng*(*)
        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)
        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.7 2005/02/24 00:14:45 willic3 Exp $
  
  c Generated automatically by Fortran77Mill on Tue May 18 14:18:50 2004
  
  