anchor.diff:            if(numslp.ne.0) call rdisp(dxcur,skew,numnp)
           end if
   c
  -c...compute contribution to global load vector from winkler boundary
  +c...  zero internal force vector prior to reconputing it.
  +c
  +        call fill(bintern,zero,neq)
  +c
  +c...compute contribution to internal force vector from winkler boundary
   c   conditions
   c
  -        if(nwink.ne.0) call winklf(btot,gvec2,iwink,wink,histry,nwink,
  -     &   nhist,nstep,neq,lastep)
  -        if(nwinkx.ne.0) call winklf(btot,gvec2,iwinkx,winkx,histry,
  +        if(nwink.ne.0) call winklf(bintern,dispvec,iwink,wink,histry,
  +     &   nwink,nhist,nstep,neq,lastep)
  +        if(nwinkx.ne.0) call winklf(bintern,dispvec,iwinkx,winkx,histry,
        &   nwinkx,nhist,nstep,neq,lastep)
   c
   c...update the total displacement and the displacement increment
  
positive.code: c
  c...  included dimension and type statements
  c
        include "nsysdat_dim.inc"
        include "npar_dim.inc"
        include "rtimdat_dim.inc"
  
negative.code: c
  c...  included variable definitions
  c
        include "nforce_def.inc"
        include "nsysdat_def.inc"
        include "npar_def.inc"
        include "rtimdat_def.inc"
  