anchor.diff:  c
   c...  open input file
   c
  -      ierr=0
  -      numbc=0
  -      nget=3
  +      ierr=izero
  +      numbc=izero
  +      nget=ithree
         open(kr,file=bcfile,status="old",err=20)
   c
   c...  get units, returning error 2 if they aren't found.
   c
  -      call get_units(kr,ierr,nget,units_defined,units,def)
  -      if(ierr.eq.2) return
  +      call get_units(kr,nget,units_defined,units,def,ierr,errstrng)
  +      if(ierr.ne.izero) return
         displacement_units=units(1)
         velocity_units=units(2)
         force_units=units(3)
  
positive.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine scan_bc(ndof,numbc,kr,ierr,displacement_units,
       & velocity_units,force_units,bcfile)
  c
  c...  subroutine to perform an initial scan of the boundary condition
  c     file to determine the number of boundary condition entries and the
  
negative.code: c     Error codes:
  c         0:  No error
  c         1:  Error opening input file
  c         3:  Read error
  c         5:  Units not specified
  c
        include "implicit.inc"
  c
  c...  parameter definitions
  c
        include "ndimens.inc"
        include "nconsts.inc"
  c
  c...  subroutine arguments
  c
        integer numbc,kr,ierr
        character displacement_units*(*),velocity_units*(*)
        character force_units*(*),bcfile*(*),errstrng*(*)
  c
  c...  local constants
  c
        character def(3)*18
        data def/"displacement_units","velocity_units","force_units"/
  