anchor.diff:          mode=iwink(1,i)
           k=iwink(2,i)
           if(mode.eq.ione) then
  -          btot(k)=btot(k)-wink(i)*deld(k)
  +          bintern(k)=bintern(k)+wink(i)*deld(k)
           else
             ihist=-mode
  -          btot(k)=btot(k)-wink(i)*deld(k)*histry(ihist,nstep+1)
  +          bintern(k)=bintern(k)+wink(i)*deld(k)*histry(ihist,nstep+1)
           end if
         end do
         return
         end
   c
   c version
  -c $Id: winklf.f,v 1.2 2004/07/12 21:09:09 willic3 Exp $
  +c $Id: winklf.f,v 1.3 2005/01/05 22:35:38 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine winklf(btot,deld,iwink,wink,histry,nwink,nhist,nstep,
       & neq,lastep)
  c
  c       program to compute winkler restoring forces from the
  c       displacements and add them to the global force matrix btot
  c       note: minus sign in expression insures that these are
  c       *restoring* forces that oppose displacment increments,
  c       deld.
  c
        include "implicit.inc"
  c
  
negative.code: c
        integer nwink,nhist,nstep,neq,lastep
        integer iwink(2,nwink)
        double precision bintern(neq),deld(neq),wink(nwink)
        double precision histry(nhist,lastep+1)
  c
  c...  local variables
  