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':

                    raise ParseError(

                        f'{block.header} has disorder, pass '

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

                        'remove/ignore disorder'

                    )

                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':

                    raise ParseError(

                        f'{block.name} has disorder, pass '

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

                        'remove/ignore disorder'

                    )

                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]