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~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine bdeldss(xl,dl,sh,shj,ee,det,gauss,n,
       & nen,nee,nsd,ndof,nstr,ngauss,ierr,getshape,bmatrix)
  c
  c...  Subroutine to compute strain increments in each element.
  
negative.code: c     being used.
  c
        include "implicit.inc"
  c
  c...  dimension parameters
  c
        include "nshape.inc"
  c
  c...  subroutine arguments
  c
        integer iel,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)
  