def Plot_Coherence_Average(): fig = plt.figure(figsize=(10, 4)) cohere = [] freq = Coherence_Frequency() for ibasin in xrange(0, 10): # Plot coherence cohere_basin = load('%s/coherence_obs_5var_good_station_%s' % (workspace, basinlongs[ibasin])) # cohere_basin = load('%s/coherence_obs_5var_test2048_%s' %(workspace, basinlongs[ibasin])) cohere.append(cohere_basin) # for all stations average ax = fig.add_subplot(1, 1, 1) Plotting.CoherenceVarPlot(ax, mean(vstack(cohere), axis=0), sampling_frequency, freq) plt.show() # savefig('%s/fig2_coh_obs_5var_all_good_station.tif' %figdir, dpi=300) return