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
