Python 语言顶刊绘图复现 —— 全球海洋变量均值趋势图 + 纬度统计图

def plot_mean_tr_ocean_var(    var_mean,    var_trend,    var_name='Salinity',    units='psu',    trend_units='psu decade$^{-1}$',    cmap={'mean': 'viridis', 'trend': 'RdBu_r'},    mean_norm=mpl.colors.Normalize(vmin=30, vmax=37),    trend_range=(-0.05, 0.05),    save_path=None):    fig = plt.figure(figsize=(30, 15), dpi=600)    fig, axs = plt.subplots(nrows=4, ncols=1,subplot_kw={'projection': ccrs.Robinson(central_longitude=180)})    ccrs_land = feature.NaturalEarthFeature('physical', 'land', '50m',                                        edgecolor='darkgrey',                                        facecolor='darkgrey',                                        linewidth=0.2)    # ================= Mean map =================    im_mean = var_mean.plot(ax=axs[0],transform=ccrs.PlateCarree(),add_colorbar=False,cmap=cmap['mean'],norm=mean_norm,rasterized=True)    axs[0].set_extent([0.1, 359.99, -85, 85], crs=ccrs.PlateCarree())    axs[0].add_feature(ccrs_land)    axs[0].set_title(f"a  {var_name}",loc='left', fontsize=15, fontweight='bold', fontname='Arial')    # colorbar (mean)    cax = fig.add_axes([0.41, 0.7226, 0.005, 0.15])    cbar = fig.colorbar(im_mean, cax=cax, extend='both')    cbar.set_label(units, fontsize=13, fontname='Arial', labelpad=-55)    cbar.ax.tick_params(labelsize=12)    cbar.ax.yaxis.tick_left()    # ================= Trend map =================    ax_tr = plt.axes([0.4183, 0.49, 0.19, 0.168], projection=ccrs.Robinson(central_longitude=180))    im_tr = var_trend.plot(ax=ax_tr, transform = ccrs.PlateCarree(), add_colorbar=False, cmap=cmap['trend'], vmin=trend_range[0], vmax=trend_range[1], rasterized=True)    ax_tr.set_extent([0.1, 359.99, -85, 85], crs=ccrs.PlateCarree())    ax_tr.add_feature(ccrs_land)    ax_tr.set_title(f"d  {var_name} trend", loc='left', fontsize=15, fontweight='bold', fontname='Arial')    cax = fig.add_axes([0.41, 0.50, 0.005, 0.15])    cbar = fig.colorbar(im_tr, cax=cax, extend='both')    cbar.set_label(trend_units, fontsize=13, fontname='Arial', labelpad=-65)    cbar.ax.tick_params(labelsize=12)    cbar.ax.yaxis.tick_left()    # ================= Zonal mean =================    ax_zm1 = plt.axes([0.615, 0.7126, 0.05, 0.167])    zm = var_mean.mean(axis=1)    y = var_mean.lat    ax_zm1.fill_betweenx(y, 0, zm, color=mpl.colormaps[cmap['mean']](0.6))    ax_zm1.set_yticks(range(-80,81,20))    ax_zm1.set_yticklabels(['80$\u00b0$S','60$\u00b0$S','40$\u00b0$S','20$\u00b0$S','0$\u00b0$','20$\u00b0$N','40$\u00b0$N','60$\u00b0$N','80$\u00b0$N'],fontname = 'Arial', fontsize = 12)    ax_zm1.set_title("b  Salinity", loc='left', fontsize=15, fontweight='bold')    ax_zm1.set_xlabel(units)    ax_zm1.axvline(0, color='k')    ax_zm2 = plt.axes([0.615, 0.49, 0.05, 0.167])    zm_tr = var_trend.mean(axis=1)    y = var_trend.lat    cmap_tr = mpl.colormaps[cmap['trend']]    pos = np.ma.masked_where(zm_tr < 0, zm_tr)    ax_zm2.fill_betweenx(y, 0, pos, color=cmap_tr(0.85), linewidth=1,alpha=1)    neg = np.ma.masked_where(zm_tr > 0, zm_tr)    ax_zm2.fill_betweenx(y, 0, neg, color=cmap_tr(0.15), linewidth=1, alpha=1)    ax_zm2.set_yticks(range(-80,81,20))    ax_zm2.set_yticklabels(['80$\u00b0$S','60$\u00b0$S','40$\u00b0$S','20$\u00b0$S','0$\u00b0$','20$\u00b0$N','40$\u00b0$N','60$\u00b0$N','80$\u00b0$N'],fontname = 'Arial', fontsize = 12)    ax_zm2.set_title(f"e  {var_name} trend", loc='left', fontsize=15, fontweight='bold', fontname='Arial')    ax_zm2.set_xlabel(trend_units, fontsize=13, fontname='Arial')    ax_zm2.axvline(0, color='k', linewidth=0.8)    for i in [1, 2, 3]:        fig.delaxes(axs[i])    if save_path:        plt.savefig(save_path, dpi=600, bbox_inches='tight')    plt.show()