'bowl1': bowl1_p, 'bowl2': bowl2_p, 'labBowl': labBowl} varAtl = {'name': 'Atlantic', 'var_change': var_change_a, 'var_mean': None, 'bowl1': bowl1_a, 'bowl2': bowl2_a, 'labBowl': labBowl} varInd = {'name': 'Indian', 'var_change': var_change_i, 'var_mean': None, 'bowl1': bowl1_i, 'bowl2': bowl2_i, 'labBowl': labBowl} # ------------------------------------ # Plot # ------------------------------------ fig, axes = plt.subplots(nrows=2, ncols=3, figsize=(17,5)) if name == 'Durack & Wijffels': levels = np.linspace(minmax[0], minmax[1], minmax[2]) cmap = custom_div_cmap() cnplot = zonal_2D(plt, 'total', axes[0,0], axes[1,0], 'left', lat, density, varAtl, domrho, cmap, levels, clevsm, clevsm_bold) cnplot = zonal_2D(plt, 'total', axes[0,1], axes[1,1], 'mid', lat, density, varPac, domrho, cmap, levels, clevsm, clevsm_bold) cnplot = zonal_2D(plt, 'total', axes[0,2], axes[1,2], 'right', lat, density, varInd, domrho, cmap, levels, clevsm, clevsm_bold) else: levels = np.linspace(minmax[0], minmax[1], minmax[2]) cmap = custom_div_cmap() # plt.get_cmap('bwr') cnplot = zonal_2D(plt, 'total_mme', axes[0, 0], axes[1, 0], 'left', lat, density, varAtl, domrho, cmap, levels, clevsm, clevsm_bold) cnplot = zonal_2D(plt, 'total_mme', axes[0, 1], axes[1, 1], 'mid', lat, density, varPac, domrho, cmap, levels, clevsm, clevsm_bold)
'bowl2': bowl2z[3, :], 'labBowl': labBowl, 'density': density_z[3, :, :] } # ------------------------------------------------------------------------------- # Plot # ------------------------------------------------------------------------------- domzed = [0, 500, 2000] # -- Create figure and axes instances fig, axes = plt.subplots(nrows=2, ncols=3, figsize=(17, 5)) # -- color map cmap = custom_div_cmap() # -- levels levels = np.linspace(minmax[0], minmax[1], minmax[2]) # Change #levels = np.arange(33.5,35.5,0.2) # Mean salinity #levels = np.arange(-2,30,2) # Mean temperature #levels = np.arange(0,701,50) # Mean volume ext_cmap = 'both' contourDict = { 'cmap': cmap, 'levels': levels, 'ext_cmap': ext_cmap, 'isopyc': False }
plt.subplots_adjust(hspace=.0001, wspace=0.05, left=0.04, right=0.86) cb = plt.colorbar(cnplot1, ax=axes.ravel().tolist(), ticks=levels[::3]) plt.suptitle('dS/dy (%s)' %(name,), fontweight='bold', fontsize=14, verticalalignment='top') #plt.close() # ==== dS/dy change ==== fig2, axes = plt.subplots(nrows=2, ncols=3, figsize=(17,5)) levels = np.linspace(-0.1,0.1,16) cmap = custom_div_cmap() cnplot2 = zonal_2D(plt, 'spiciness_change', axes[0,0], axes[1,0], 'left', lat, density, varAtl, domrho, cmap, levels) cnplot2 = zonal_2D(plt, 'spiciness_change', axes[0,1], axes[1,1], 'mid', lat, density, varPac, domrho, cmap, levels) cnplot2 = zonal_2D(plt, 'spiciness_change', axes[0,2], axes[1,2], 'right', lat, density, varInd, domrho, cmap, levels) plt.subplots_adjust(hspace=.0001, wspace=0.05, left=0.04, right=0.86) cb = plt.colorbar(cnplot1, ax=axes.ravel().tolist(), ticks=levels[::3])