def get_atom_indices(labels, atom):
    """
    labels - a list of coordinate labels ("Elements")
    atom - the atom whose indices in labels are sought
    Returns a list of all location of [atom] in [labels]
    """
    indices = []
    for i in range(len(labels)):
        if labels[i] == atom:
            indices.append(i)
    return indices