Multi-Panel EOF Visualization

This section demonstrates how to combine EOF spatial maps and time series panels.

fig = plt.figure(figsize=(12, 8))
ax1 = fig.add_subplot(2, 2, 1, projection=ccrs.LambertConformal())
ax2 = fig.add_subplot(2, 2, 2, projection=ccrs.LambertConformal())
ax3 = fig.add_subplot(2, 1, 2)

# Plot correlation maps
ax1.pcolormesh(lon, lat, cor1, transform=ccrs.PlateCarree(), cmap="coolwarm")
ax2.pcolormesh(lon, lat, cor2, transform=ccrs.PlateCarree(), cmap="coolwarm")

# Plot PCs vs CRWS/PDSI
ax3.plot(years, pc1, label="PC1")
ax3.plot(years, crws, label="CRWS")
ax3.legend()

plt.tight_layout()
plt.savefig(fnFIG)
plt.show()