anchor.diff:  c
         subroutine formdf_ss(
        & bdisp,neq,                                                       ! force
  -     & x,d,deld,nsd,ndof,numnp,                                         ! global
  +     & x,d,deld,numnp,                                                  ! global
        & s,stemp,                                                         ! stiff
  -     & dmat,ien,lm,lmx,infiel,iddmat,nstr,nddmat,ndmatsz,numelt,nconsz, ! elemnt
  +     & dmat,ien,lm,lmx,infiel,iddmat,ndmatsz,numelt,nconsz,             ! elemnt
        & infmat,numat,                                                    ! materl
  -     & gauss,sh,shj,infetype,netypes,                                   ! eltype
  -     & skew,nskdim,numrot,                                              ! skew
  +     & gauss,sh,shj,infetype,                                           ! eltype
  +     & skew,numrot,                                                     ! skew
        & getshape,bmatrix,                                                ! bbar
  -     & ierr)                                                            ! errcode
  +     & ierr,errstrng)                                                   ! errcode
   c
   c...program to compute forces due to kinematic boundary conditions
   c***  Leave this routine essentially 'as-is' for now.  In the near
  
positive.code:       double precision sh(nsd+1,nenmax,ngaussmax,netypes)
        double precision shj(nsd+1,nenmax,ngaussmax,netypes)
        double precision skew(nskdim,numnp)
  c
  c...  defined constants
  c
        include "nconsts.inc"
        include "rconsts.inc"
  c
  c...  external routines
  c
  
negative.code: c
        include "implicit.inc"
  c
  c...  parameter definitions
  c
        include "ndimens.inc"
        include "nshape.inc"
        include "nconsts.inc"
        include "rconsts.inc"
  c
  c...  subroutine arguments
  c
        integer neq,numnp,ndmatsz,numelt,nconsz,numat,numrot,ierr
        integer ien(nconsz),lm(ndof,nconsz),lmx(ndof,nconsz)
        integer infiel(6,numelt),iddmat(nstr,nstr),infmat(6,numat)
        integer infetype(4,netypes)
        character errstrng*(*)
        double precision bdisp(neq),x(nsd,numnp),d(ndof,numnp)
        double precision deld(ndof,numnp),s(neemax*neemax)
        double precision stemp(neemax*neemax),dmat(nddmat,ndmatsz)
  