The atoms in disordered structures such as liquids and glasses are arranged in a complex manner. Their atomic arrangements are often studied using the Voronoi tessellation method. This method divides the space containing atoms into regions called Voronoi polyhedra. Each polyhedron contains one atom. From the structure of the polyhedron, we can know the bond network, or the topology, of the cluster composed of the atoms in the polyhedron and its neighbors. Thus, the topological order can be characterized by identifying the most dominant Voronoi polyhedron. For this purpose, the Voronoi index has long been used. However, it does not specify how polygons are arranged in the polyhedron. Therefore, polyhedra with different graph structures can have the same index. This problem can lead to failure in characterizing the topological order. In addition, the Voronoi index does not tell anything about the chemical arrangement, namely how different types of atoms are arranged in the cluster. To characterize both of the topological and chemical order of disordered structures, a polyhedron-naming system called the polyhedron code has been developed and implemented in the Vorotis software. This paper presents the essential features of Vorotis. Program summaryProgram Title: VorotisCPC Library link to program files:https://doi.org/10.17632/4548pnvzpy.1Developer's repository link:https://doi.org/10.5281/zenodo.6069726Code Ocean capsule:https://doi.org/10.24433/CO.1700993.v1Licensing provisions: MITProgramming language: Fortran, Csh, PythonNature of problem: Characterizing the arrangement of atoms in substances.Solution method: The polyhedron code is combined with the Voronoi tessellation method.Additional comments including restrictions and unusual features: The current version of Vorotis assumes that every vertex of a Voronoi polyhedron meets three faces.