anchor.diff:  c     scur
   c
           do l=1,ngauss
  -          indstateg=indstate+(l-1)*nstate*nstr
  -          inddmatg=inddmat+(l-1)*nddmat
  -          call elas_strs(state(1,indstateg),ee(nstr*(l-1)),
  +          call elas_strs(dstate(1,indstateg),ee(nstr*(l-1)),
        &     dmat(1,inddmatg),nstate)
  -          call dcopy(nstr,state(1,indstateg),ione,scur(nstr*(l-1)),ione)
  +          call dcopy(nstr,dstate(1,indstateg),ione,scur(nstr*(l-1)),
  +     &     ione)
  +          indstateg=indstateg+incstate
  +          inddmatg=inddmatg+nddmat
           end do
   c
   c...  compute equivalent nodal loads and add them to global load
  
positive.code:       imatvar=infmat(4)
        ngtest=izero
        if(imatvar.eq.izero) ngtest=ngaussmax
  c
  c...  loop over elements in a material group
  c
  
negative.code:       end
  c
  c version
  c $Id: elas_strs_mat_cmp_ss.f,v 1.3 2004/07/02 17:57:23 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  