anchor.diff:  c
    20   continue
           ierr=1
  +        errstrng="read_coords"
           close(kr)
           return
   c
  
positive.code:       end
  c
  c version
  c $Id: read_coords.f,v 1.1 2004/04/14 21:18:30 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code: c...  output coordinates to plot file
  c
        if(idsk.eq.0) then
          open(kp,file=pfile,err=50,status="old",access="append")
  	do i=1,numnp
  	  write(kp,"(i7,3(1pe16.8))",err=60) i,(x(j,i),j=1,nsd)
          end do
        else
          open(kp,file=pfile,err=50,status="old",access="append",
       &   form="unformatted")
  	write(kp,err=60) x
        end if
        close(kp)
  c
  c...  output coordinates to human-readable file
  c
        if(idout.gt.izero) then
          open(kw,file=ofile,err=50,status="old",access="append")
  	do i=1,numnp
  	  if(i.eq.1.or.mod(i,npage).eq.0) write(kw,1000)
       &     (labelc(j),j=1,nsd)
  	  write(kw,"(6x,i7,10x,3(1pe20.8))",err=60) i,(x(j,i),j=1,nsd)
  	end do
          close(kw)
        end if
  