anchor.diff:  c
   c...  localize displacement BC
   c
  -          call ldisbc(dld,deld,ien(indien),lm(1,indien),ndof,nen,numnp)
  +          call ldisbc(dld,deld,ien(indien),lm(1,indien),nen,numnp)
             do i=1,ndof*nen
               if(dld(i).ne.zero) go to 100
             end do
  
positive.code: c
        subroutine formdf_ss(
       & bdisp,neq,                                                       ! force
       & x,d,deld,nsd,ndof,numnp,                                         ! global
       & s,stemp,                                                         ! stiff
       & dmat,ien,lm,lmx,infiel,iddmat,nstr,nddmat,ndmatsz,numelt,nconsz, ! elemnt
       & infmat,numat,                                                    ! materl
       & gauss,sh,shj,infetype,netypes,                                   ! eltype
       & skew,nskdim,numrot,                                              ! skew
       & getshape,bmatrix,                                                ! bbar
       & ierr)                                                            ! errcode
  c
  c...program to compute forces due to kinematic boundary conditions
  c***  Leave this routine essentially 'as-is' for now.  In the near
  
negative.code:           ngaussdim=max(ngauss,ngtest)
  c
            call formes_ss(
       &     x,numnp,                                                     ! global
       &     s,stemp,                                                     ! stiff
       &     dmat(1,inddmat),ien(indien),lm(1,indien),iddmat,iel,         ! elemnt
       &     gauss(1,1,ietype),sh(1,1,1,ietype),shj(1,1,1,ietype),        ! eltype
       &     ngauss,ngaussdim,nen,nee,                                    ! eltype
       &     skew,numrot,                                                 ! skew
       &     getshape,bmatrix,                                            ! bbar
       &     ierr,errstrng)                                               ! errcode
  c
            if(ierr.ne.izero) return
  c
  