anchor.diff:  c
   c
         subroutine unlock(b,btot,id,iwink,idhist,ibond,bond,histry,
  -     & nstep,ndof,numnp,nwink,nhist,neq,numdif,lastep,iflag)
  +     & nstep,numnp,nwink,nhist,neq,numdif,lastep,iflag)
   c
   c       program to remove the winkler forces and leave the applied
   c       forces at a specified step.  code uses load history factors
  
positive.code: c
        do i=1,ndof
          do j=1,numnp
            if(id(i,j).ne.0) then
              do k=1,nwink
                idk=iwink(2,k)
                if(id(i,j).eq.idk) then
                  mode=iwink(1,k)
                  if(mode.lt.0) then
                    ihist=-mode
                    diff=histry(ihist,nstep+1)-histry(ihist,nstep)
                    if(diff.lt.-cutoff) then
                      if(iflag.eq.1) then
                        idforc=ibond(i,j)
                        forc=bond(i,j)
                      else if(iflag.eq.2) then
                        idforc=idhist(j)
                        forc=zero
                        if(numdif.ne.0) forc=bond(i,j)
                      end if
                      if(idforc.eq.0) then
                        b(idk)=-btot(idk)+forc
                      else
                        b(idk)=-btot(idk)+forc*histry(idforc,nstep+1)
  
negative.code:       end
  c
  c version
  c $Id: unlock.f,v 1.2 2004/07/12 21:04:56 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  