anchor.diff:  c...  subroutine arguments
   c
         integer iel,nen,nee,ngauss,ierr
  +      character errstrng*(*)
         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)
  
positive.code: c
  c
        subroutine bdeld_ss(xl,dl,sh,shj,ee,det,gauss,iel,nen,nee,ngauss,
       & ierr,getshape,bmatrix)
  c
  c...  Subroutine to compute strain increments in each element.
  c     Strain is evaluated at the number of gauss points specified by
  
negative.code: c   points, and then compute average dilatational component.
  c
        call getshape(xl,sh,shj,shd,shbar,det,gauss,vol,iel,nen,ngauss,
       & ierr,errstrng)
        if(ierr.ne.0) return
  c
  c...construct b matrix and compute strains at gauss point(s)
  