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