integer, dimension(num_symmetry_operations,
current_cell%max_ions_in_species,current_cell%num_species)
:: symmetry_equiv_atoms

This contains the indexes, within each species, of ions related to one another by the symmetry operations in symmetry_operations.