def load_mean_distances_intact(filein='dists.dat'):
    """Load and return mean distances of the reactant/intact-bond minimum.

    Parameters
    ----------
    filein: str
        Name of the file created with method *save_mean_distances* of
        class *Dissociation*.

    Returns
    -------
    dists: list of floats
        Mean distances.
    forces: list of floats
        External forces.

    """
    dists = []
    forces = []
    fd = open(filein)
    first = True
    for line in fd:
        if first:
            first = False
            continue
        line_arr = line.rstrip().split()
        forces.append(float(line_arr[0]))
        dists.append(float(line_arr[1]))
    fd.close()
    return dists, forces
