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
   20   continue
          ierr=1
          close(kr)
          return
  c
  
negative.code: c     is pretty sloppy, assuming that they were already specified during
  c     the scan phase.
  c
        ierr=izero
        open(kr,file=cfile,status="old",err=20)
        call pskip(kr)
        read(kr,"(a80)") dummy
        i=index(dummy,"=")
        if(i.eq.izero) then
          ierr=ifive
          errstrng="read_coords"
          return
        end if
        call pskip(kr)
  