anchor.diff:  c
             call fill(p,zero,nee)
             call dsymv("u",nee,one,s,nee,dld,ione,zero,p,ione)
  -          call addfor(bintern,p,lm(1,indien),lmx(1,indien),neq,nee)
  +          call addfor(bintern,p,lm(1,ielg),lmx(1,ielg),neq,nee)
    150      continue
           end do
  -        matgpt=matgpt+nmatel
         end do
  -cdebug      write(6,*) "bintern:"
  -cdebug      write(6,*) (bintern(idb),idb=1,200)
         return
         end
   c
   c version
  -c $Id: formdf_ss.f,v 1.7 2005/02/23 23:57:20 willic3 Exp $
  +c $Id: formdf_ss.f,v 1.8 2005/03/21 19:12:43 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code: c
  c...  local variables
  c
        integer matgpt,imat,matmodel,nmatel,imatvar,ngtest,ielg,iel,indien
        integer ietype,nen,inddmat,ngauss,nee,ngaussdim,i
        double precision p(60),dld(60)
  cdebug      integer idb
  c
  cdebug      write(6,*) "Hello from formdf_ss_f!"
  c
        matgpt=1
  c
  c...  loop over material groups
  c
        do imat=1,numat
          matmodel=infmat(1,imat)
          nmatel=infmat(2,imat)
          imatvar=infmatmod(4,matmodel)
          ngtest=0
          if(imatvar.eq.izero) ngtest=ngaussmax
          do ielg=matgpt,matgpt+nmatel-1
            iel=infiel(4,ielg)
            indien=infiel(1,iel)
            ietype=infiel(3,iel)
            nen=infetype(2,ietype)
  c
  c...  localize displacement BC
  c
            call ldisbc(dld,deld,ien(indien),lm(1,indien),nen,numnp)
            do i=1,ndof*nen
              if(dld(i).ne.zero) go to 100
            end do
  
negative.code: c
  c...  form element stiffness matrix
  c
            call formes_ss(
       &     x,numnp,                                                     ! global
       &     s,stemp,                                                     ! stiff
       &     dmat(1,ielg),ien(1,ielg),lm(1,ielg),ielg,                    ! elemnt
       &     gauss,sh,shj,nen,ngauss,nee,                                 ! eltype
       &     skew,numrot,                                                 ! skew
       &     getshape,bmatrix,                                            ! bbar
       &     ierr,errstrng)                                               ! errcode
  