dwiddo / average-minimum-distance

        remove_sites.extend((i for i, n in enumerate(asym_types) if n == 1))

    # Remove atoms with fractional occupancy or raise ParseError

    if disorder != 'all_sites':

        for i, dis in enumerate(has_disorder):

            if i in remove_sites:

                continue

            if dis:

                if disorder == 'skip':

                    msg = f"{block.header} has disorder, pass " \

                            "disorder='ordered_sites' or 'all_sites' to " \

                            "remove/ignore disorder"

                    raise ParseError(msg)

                elif disorder == 'ordered_sites':

                    remove_sites.append(i)

    # Asymmetric unit

    asym_unit = [c for i, c in enumerate(asym_unit) if i not in remove_sites]

        remove_sites.extend(i for i, num in enumerate(asym_types) if num == 1)

    # Remove atoms with fractional occupancy or raise ParseError

    if disorder != 'all_sites':

        for i, dis in enumerate(has_disorder):

            if i in remove_sites:

                continue

            if dis:

                if disorder == 'skip':

                    msg = f"{block.name} has disorder, pass " \

                           "disorder='ordered_sites'or 'all_sites' to " \

                           "remove/ignore disorder"

                    raise ParseError(msg)

                elif disorder == 'ordered_sites':

                    remove_sites.append(i)

    # Asymmetric unit

    asym_unit = [c for i, c in enumerate(asym_unit) if i not in remove_sites]