anchor.diff:  c     being used.
   c
         include "implicit.inc"
  +c
  +c...  dimension parameters
  +c
         include "nshape.inc"
   c
   c...  subroutine arguments
   c
  -      integer n,nen,nee,nsd,ndof,nstr,ngauss,ierr
  +      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)
  
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:       end
  c
  c version
  c $Id: bdeld_ss.f,v 1.2 2004/06/16 16:57:13 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  