anchor.diff:        end
   c
   c version
  -c $Id: formf_ss.f,v 1.1 2004/06/16 21:14:06 willic3 Exp $
  +c $Id: formf_ss.f,v 1.2 2004/06/21 21:03:39 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
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)
  