def _get_epoch(obsdate, limits):
    """
    Returns the epoch of the observation with date obsdate 
    given the epochs dates defined in limits.
    """
    epoch = f"e{len(limits) + 1}"
    for i, d in enumerate(obsdate < limits):
        if d:
            epoch = f"e{i + 1}"
            break

    return epoch