anchor.diff:          read(kr,"(a80)") dummy
           j=index(dummy,"=")
           if(j.eq.0) then
  -          ierr=2
  +          ierr=5
  +          errstrng="read_bc"
             return
           end if
         end do
  
positive.code: c
  20    continue
          ierr=1
          close(kr)
          return
  c
  
negative.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine read_bc(bond,dscale,vscale,fscale,ibond,id,numnp,
       & numbc,neq,kr,kw,idout,bcfile,ofile,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
  