anchor.diff:  c
   c...  local variables
   c
  -      integer l,i,nee
  +      integer l
         double precision shd(4,nenmax,ngaussmax),b(6,3*nenmax)
         double precision shbar(4,nenmax),vol
   c
  
positive.code: c     being used.
  c
        include "implicit.inc"
        include "nshape.inc"
  c
  c...  subroutine arguments
  c
        integer n,nen,nee,nsd,ndof,nstr,ngauss,ierr
        double precision xl(nsd,nen),dl(ndof*nen)
        double precision sh(nsd+1,nenmax,ngaussmax)
        double precision shj(nsd+1,nenmax,ngaussmax),ee(nstr,ngauss)
  
negative.code: c
  cdebug      write(6,*) "Hello from bdeld_ss_f!"
  c
        call getshape(xl,sh,shj,shd,shbar,det,gauss,vol,iel,nen,nsd,
       & ngauss,ierr)
        if(ierr.ne.0) return
  c
  c...construct b matrix and compute strains at gauss point(s)
  c
        do l=1,ngauss
          call bmatrix(b,shd(1,1,l),shbar,nsd,ndof,nen,nstr)
  