anchor.diff:        imatvar=infmat(4)
         ngtest=izero
         if(imatvar.eq.izero) ngtest=ngaussmax
  +      incstate=nstr*nstate
   c
   c...  loop over elements in a material group
   c
  
positive.code:       end
  c
  c version
  c $Id: td_strs_mat_cmp_ss.f,v 1.1 2004/06/25 15:51:51 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code: c     scur
  c
          do l=1,ngauss
            call td_strs_mat(state(1,indstateg),dstate(1,indstateg),
       &     ee(nstr*(l-1)),dmat(1,inddmatg),prop,rtimdat,rgiter,ntimdat,
       &     nstate,matchg)
            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
  