anchor.diff:        end
   c
   c version
  -c $Id: read_bc.f,v 1.2 2004/07/07 20:52:31 willic3 Exp $
  +c $Id: read_bc.f,v 1.3 2005/01/06 01:15:31 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code: c
  c...  open input file and define scaling factors
  c
        ierr=0
        open(kr,file=bcfile,status="old",err=20)
        scale(1)=zero
        scale(2)=dscale
  
negative.code:           ihist(j)=imode/10
            itype(j)=imode-10*ihist(j)
            bond(j,n)=scale(itype(j)+1)*bond(j,n)
            if(itype(j).eq.ithree) nconcforce=nconcforce+ione
          end do
        end do
        close(kr)
  