anchor.diff:          close(kr)
           return
   c
  -1000  format(1x,///,'  n o d a l   f o r c e s   a n d   d i s p l a',
  -     1 ' c e m e n t s',//,
  -     2 '      key to boundary condition codes:',//,
  -     3 '          free = unconstrained degree of freedom; the',/,
  -     4 '                 numerical value is meaningless.',/,
  -     5 '          disp = fixed displacement',/,
  -     6 '          velo = constant velocity',/,
  -     7 '          forc = applied (constant) force',//,
  -     8 '      hfac = load history factor applied (0 if none)',///,
  -     9 ' node # ',3(4x,a4,6x,'type  hfac  '))
  -2000  format(' ')
  -3000  format(1x,i7,2x,3(1pe12.5,2x,a4,1x,i3,4x))
         end
   c
   c version
  -c $Id: read_bc.f,v 1.4 2005/03/12 01:59:39 willic3 Exp $
  +c $Id: read_bc.f,v 1.5 2005/04/12 22:33:27 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code: c
  c                             Charles A. Williams
  c                       Rensselaer Polytechnic Institute
  c                        (C) 2004  All Rights Reserved
  c
  c  Copyright 2004 Rensselaer Polytechnic Institute.
  c  All worldwide rights reserved.  A license to use, copy, modify and
  c  distribute this software for non-commercial research purposes only
  c  is hereby granted, provided that this copyright notice and
  
negative.code: c
  c
        subroutine read_bc(bond,dscale,vscale,fscale,ibond,numnp,
       & numbc,nconcforce,kr,bcfile,ierr,errstrng)
  c
  c...  subroutine to read in boundary conditions.  The bc types are
  c     stored in the ibond array and the bc value is stored in the bond
  