anchor.diff:  c
   c
         subroutine winklr(alnz,iwink,wink,histry,nstep,nwink,nhist,nnz,
  -     & lastep,idout,kto,kw)
  +     & lastep,ierr,errstrng)
   c
   c...program to implement winkler restoring forces on nodes
   c   designated by iwink.  this program adds a constant given
  
positive.code: c
        include "implicit.inc"
  c
  c...  subroutine arguments
  c
        integer nstep,nwink,nhist,nnz,lastep,idout,kto,kw
        integer iwink(2,nwink)
        double precision alnz(nnz),wink(nwink),histry(nhist,lastep+1)
  c
  c...  local variables
  c
  
negative.code:       do i=1,nwink
          mode=iwink(1,i)
          k=iwink(2,i)
          if(mode.eq.ione) then
            alnz(k)=alnz(k)+wink(i)
          else
            ihist=-mode
            if(ihist.gt.nhist) then
              ierr=100
              errstrng="winklr"
              return
            end if
            alnz(k)=alnz(k)+wink(i)*histry(ihist,nstep+1)
          end if
        end do
        return
        end
  c
  c version
  c $Id: winklr.f,v 1.2 2004/07/12 21:12:54 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  