# =============================================================================
# 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',