anchor.diff:  c...  local variables
   c
         integer nmatel,nstate,ind,iel,indien,ietype,indstate,inddmat
  -      integer ngauss,nen,nee,l,indstateg,inddmatg
  +      integer ngauss,nen,nee,l,indstateg,inddmatg,incstate
         double precision dl(60),xl(60),scur(162),ee(162),p(60),det(27)
   c
   c...  included variable definitions
  
positive.code:         ngauss=infetype(1,ietype)
          nen=infetype(2,ietype)
          nee=infetype(4,ietype)
  c
  c...  localize coordinates and displacements
  c
  
negative.code: c     scur
  c
          do l=1,ngauss
            call td_strs(state(1,indstateg),dstate(1,indstateg),
       &     ee(nstr*(l-1)),dmat(1,inddmatg),prop,rtimdat,rgiter,ntimdat,
       &     nstate)
            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
  