      subroutine xyzko (i,xyz)                                          XYZK001
c ----------------------------------------------------------------------XYZK002
c     returns atom coordinate string                                    XYZK003
c ----------------------------------------------------------------------XYZK004
      character xyz*8                                                   XYZK005
      write(*,'(/a,a,a$)') '  ',char(87+i),' coordinate : '             XYZK006
      read(*,'(a)') xyz                                                 XYZK007
      if (xyz.eq.' ') xyz='0'                                           XYZK008
      return                                                            XYZK009
      end                                                               XYZK010
                                                                        XYZK011
                                                                        XYZK012
