anchor.diff:  c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   c
   c
  -      subroutine read_coords(x,cscale,nsd,numnp,kr,kw,kp,idout,idsk,
  -     & ierr,cfile,ofile,pfile)
  +      subroutine read_coords(x,cscale,numnp,kr,kw,kp,idout,idsk,
  +     & cfile,ofile,pfile,ierr,errstrng)
   c
   c...  subroutine to read in nodal coordinates and store them in the x
   c     array.  Note that the current implementation would allow a 3D
  
positive.code: c     Error codes:
  c         0:  No error
  c         1:  Error opening input file
  c         2:  Units not specified
  c         3:  Read error
  c         4:  Wrong number of nodes
  c
        include "implicit.inc"
  c
  c...  subroutine arguments
  c
        integer nsd,numnp,kr,kw,kp,idout,idsk,ierr
        double precision cscale
        double precision x(nsd,numnp)
        character cfile*(*),ofile*(*),pfile*(*)
  c
  c...  included dimension and type statements
  c
        include "labelc_dim.inc"
  c
  c...  defined constants
  c
        integer npage
        parameter(npage=50)
  c
  c...  intrinsic functions
  c
  
negative.code: c
   20   continue
          ierr=1
          errstrng="read_coords"
          close(kr)
          return
  c
  