anchor.diff:  c*     &   (dfault(idb,l),idb=1,ndof)
         end do
         return
  -1000  format(//' fatal split node error!'//
  -     & ' attempt to use undefined load history # ',i5,
  -     & ' at split node # ',i7)
         end
   c
   c version
  -c $Id: loadf.f,v 1.1 2004/04/14 21:18:30 willic3 Exp $
  +c $Id: loadf.f,v 1.2 2004/07/07 17:52:46 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code: c
        include "implicit.inc"
  c
  c...  subroutine arguments
  c
        integer nstep,ndof,numfn,nhist,lastep,idout,kto,kw
        integer nfault(3,numfn)
        double precision fault(ndof,numfn),dfault(ndof,numfn)
        double precision histry(nhist,lastep+1),deltp
  c
  c...  defined constants
  c
        include "rconsts.inc"
  c
  c...  local variables
  c
  cdebug      integer idb
        integer l,j,ihist
  c
  cdebug      write(6,*) "Hello from loadf_f!"
  c*      write(6,"(8i7)") nstep,ndof,numfn,nhist,lastep,idout,kto,kw
  c
        call fill(dfault,zero,ndof*numfn)
        do l=1,numfn
          ihist=nfault(3,l)
          if(ihist.gt.nhist) then
            if(idout.gt.1) write(kw,1000) ihist,nfault(2,l)
            write(kto,1000) ihist,nfault(2,l)
            stop
          end if
          do j=1,ndof
            if(ihist.eq.-1) then
              dfault(j,l)=deltp*fault(j,l)
            else if((ihist.eq.0).and.(nstep.eq.0)) then
              dfault(j,l)=fault(j,l)
            else if(ihist.gt.0) then
              dfault(j,l)=fault(j,l)*histry(ihist,nstep+1)
            end if
          end do
  
negative.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine loadf(fault,dfault,histry,deltp,nfault,nstep,
       & numfn,nhist,lastep,ierr,errstrng)
  c
  c...program to compute current split node displacements
  c
  