def ob_readXYZs(filename):
    """
    Ready stacked xyz's from file.

    Args:
        filename (str): Filepath.

    Returns:
        tuple: [elements,coords]
        
        - elements (list): Coordinate list of shape (Molecules,Atoms,3).
        - coords (list): Molecule list of shape (Molecules,Atoms).
    """
    infile = open(filename, "r")
    coords = [[]]
    elements = [[]]
    for line in infile.readlines():
        if len(line.split()) == 1 and len(coords[-1]) != 0:
            coords.append([])
            elements.append([])
        elif len(line.split()) == 4:
            elements[-1].append(line.split()[0].capitalize())
            coords[-1].append([float(line.split()[1]), float(line.split()[2]), float(line.split()[3])])
    infile.close()
    return elements, coords