anchor.diff:  c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   c
   c
  -      subroutine scan_coords(nsd,numnp,kr,ierr,coord_units,cfile)
  +      subroutine scan_coords(numnp,kr,coord_units,cfile,ierr,errstrng)
   c
   c...  subroutine to perform an initial scan of the nodal coordinates
   c     file to determine the number of nodal points and the units of
  
positive.code: c
  c...  open input file
  c
        ierr=0
        numnp=0
        nget=1
        open(kr,file=cfile,status="old",err=20)
  c
  c...  get coordinate 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
        coord_units=units(1)
  c
  c... scan the file, counting the number of entries.
  
negative.code:  20   continue
          close(kr)
  	ierr=1
          errstrng="scan_coords"
          return
  c
  c...  read error
  c
   30   continue
          close(kr)
          ierr=3
          errstrng="scan_coords"
          return
  c
        end
  c
  c version
  c $Id: scan_coords.f,v 1.2 2004/07/12 19:44:14 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  