# ============================================================================= # Plot # ============================================================================= # Setup vmax = 10 levels_tau = levels_olr = np.linspace(-vmax, vmax, 21, endpoint=True) levels_r2 = np.linspace(0, 1, 21, endpoint=True) # Generate the base plot fig, axs = plt.subplots(2, 1, figsize=(8, 4)) lon2, lat2 = np.meshgrid(y_tau.lon, y_tau.lat) m = [] for i in range(2): m.append( Basemap(projection='merc',llcrnrlat=-30,urcrnrlat=30,\ llcrnrlon=100,urcrnrlon=300,lat_ts=5,resolution='c', ax=axs[i])) x, y = m[i](lon2, lat2) m[i].drawparallels(np.arange(-90., 120., 15.), labels=[1, 0, 0, 0], color='grey') m[i].drawmeridians(np.arange(0., 360., 30.), labels=[0, 0, 0, 1], color='grey') m[i].drawmapboundary(fill_color='white') m[i].drawcoastlines() axs[0].text(9.2e5, 5.3e6, 'a', weight='bold',