anchor.diff:  c
             else if(itype.eq.ithree) then
               k=id(i,j)
  -            if(nstep.eq.izero) then
  -              b(k)=bond(i,j)
  -              if(ihist.ne.izero) b(k)=bond(i,j)*histry(ihist,1)
  -            else if(nstep.ge.ione) then
  -              if(ihist.ne.izero) b(k)=bond(i,j)*(histry(ihist,nstep+1)
  -     &                                 -histry(ihist,nstep))
  -            end if
  +            bconcforce(k)=bond(i,j)
  +            if(ihist.ne.izero) bconcforce(k)=bond(i,j)*
  +     &       histry(ihist,nstep+1)
  +c*            if(nstep.eq.izero) then
  +c*              b(k)=bond(i,j)
  +c*              if(ihist.ne.izero) b(k)=bond(i,j)*histry(ihist,1)
  +c*            else if(nstep.ge.ione) then
  +c*              if(ihist.ne.izero) b(k)=bond(i,j)*(histry(ihist,nstep+1)
  +c*     &                                 -histry(ihist,nstep))
  +c*            end if
             end if
           end do
         end do
  
positive.code:       end
  c
  c version
  c $Id: load.f,v 1.3 2004/08/12 01:51:31 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code: c
  c...  subroutine arguments
  c
        integer numnp,neq,nconcflag,nhist,nstep,lastep,ierr
        integer id(ndof,numnp),ibond(ndof,numnp)
        character errstrng*(*)
        double precision bond(ndof,numnp),bconcforce(neq*nconcflag)
        double precision d(ndof,numnp),deld(ndof,numnp)
        double precision histry(nhist,lastep+1),deltp
  c
  c...  local variables
  c
  