def add_location(ax, name, lon, lat, crs, text_color=None):
    """
    Add a geodesic location to the plot axis
    """
    import matplotlib.patheffects as pe
    xy_point = xy_text = lon, lat
    ax.plot(*xy_point, marker='o', markersize=8, color='white')
    color = 'white'
    path_effects = None
    if text_color:
        color = text_color
        ax.plot(*xy_point, marker='o', markersize=4, color=text_color)
        patheffects = [pe.withStroke(linewidth=3.0, foreground='white')]
    ax.annotate(name, xy=xy_text, xycoords=crs._as_mpl_transform(ax), color
        =color, ha='left', va='center', path_effects=patheffects)
    return ax
