get_cmap("terrain") ocean_color = cmap(0.18) cmap, norm = colors.from_levels_and_colors(topo_clevs, cmap(np.linspace(0.3, 1, len(topo_clevs) - 1))) add_rectangle(ax, xx, yy, margin=20, edge_style="solid") add_r