anchor.diff:  c
         integer nmatel,nstate,ind,iel,indien,ietype,indstate,inddmat
         integer ngauss,nen,nee,l,indstateg,inddmatg,imatvar,ngtest
  -      integer ngaussdim
  +      integer ngaussdim,incstate
         double precision dl(60),xl(60),scur(162),ee(162),p(60),det(27)
   c
   c...  included variable definitions
  
positive.code: 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)),
       &     dmat(1,inddmatg),nstate)
            call dcopy(nstr,state(1,indstateg),ione,scur(nstr*(l-1)),ione)
          end do
  c
  c...  compute equivalent nodal loads and add them to global load
  
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
  