def select_bounds_around_rivermouth(pair,border=10):
    """Select a set of bounds for a plot of a river mouth"""
    imin = max(min(pair[0].get_lat(),pair[1].get_lat()) - border,0)
    imax = max(pair[0].get_lat(),pair[1].get_lat()) + border
    jmin = max(min(pair[0].get_lon(),pair[1].get_lon()) - border,0)
    jmax = max(pair[0].get_lon(),pair[1].get_lon()) + border
    return imin,imax,jmin,jmax