anchor.diff:  c...  form element stiffness matrix
   c
           call formes_ss(
  -     &   x,nsd,ndof,numnp,                                              ! global
  +     &   x,numnp,                                                       ! global
        &   s,stemp,                                                       ! stiff
  -     &   dmat(1,inddmat),ien(indien),lm(1,indien),iddmat,nstr,nddmat,   ! elemnt
  +     &   dmat(1,inddmat),ien(indien),lm(1,indien),iddmat,               ! elemnt
        &   iel,                                                           ! elemnt
        &   gauss(1,1,ietype),sh(1,1,1,ietype),shj(1,1,1,ietype),          ! eltype
        &   ngauss,ngaussdim,nen,nee,                                      ! eltype
  -     &   skew,nskdim,numrot,                                            ! skew
  +     &   skew,numrot,                                                   ! skew
        &   getshape,bmatrix,                                              ! bbar
  -     &   ierr)                                                          ! errcode
  +     &   ierr,errstrng)                                                 ! errcode
   c
           if(ierr.ne.izero) return
   c
  -        call lflteq(dl,dfault(1,i),nfault(1,i),ien(indien),nen,ndof)
  +        call lflteq(dl,dfault(1,i),nfault(1,i),ien(indien),nen)
   	call dsymv("u",nee,one,s,nee,dl,ione,zero,p,ione)
           call addfor(bdisp,p,lm(1,indien),lmx(1,indien),neq,nee)
         end do
  
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
        subroutine formf_ss(
       & bdisp,neq,                                                       ! force
       & x,numnp,                                                         ! global
       & s,stemp,                                                         ! stiff
       & dmat,ien,lm,lmx,infiel,iddmat,ndmatsz,numelt,nconsz,             ! elemnt
       & infmat,numat,                                                    ! materl
       & gauss,sh,shj,infetype,                                           ! eltype
       & skew,numrot,                                                     ! skew
       & nfault,dfault,tfault,numfn,                                      ! split
       & getshape,bmatrix,                                                ! bbar
       & ierr,errstrng)                                                   ! errcode
  c
  c      generates forces due to faulted nodes
  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,numfn,ierr
        integer ien(nconsz),lm(ndof,nconsz),lmx(ndof,nconsz)
        integer infiel(6,numelt),iddmat(nstr,nstr),infmat(6,numat)
        integer infetype(4,netypes),nfault(3,numfn)
        character errstrng*(*)
        double precision bdisp(neq),x(nsd,numnp),s(neemax*neemax)
        double precision stemp(neemax*neemax),dfault(ndof,numfn)
        double precision tfault(ndof,numfn),dmat(nddmat,ndmatsz)
  