Skip to content
Snippets Groups Projects
Commit 0966d8e3 authored by Felix Kleinert's avatar Felix Kleinert
Browse files

update rank hist plot

parent ed4bf2ba
Branches
No related tags found
2 merge requests!474Draft: Resolve "DataHandler with multiple stats per variable",!466Draft: Resolve "Include CRPS analysis and other ens verif methods or plots"
Pipeline #109553 passed
......@@ -1333,11 +1333,14 @@ class PlotRankHistogram(AbstractPlotClass):
pdf_pages = matplotlib.backends.backend_pdf.PdfPages(plot_path)
for ah in value[self.ahead_dim]:
fig, ax = plt.subplots()
ax.bar(value.sel({"rank_hist_type": "bins", self.ahead_dim: ah.values}),
value.sel({"rank_hist_type": "freq", self.ahead_dim: ah.values}),
align='center')
ax.set_xlabel(f"{self.target_var} (in {self.target_unit})")
ax.set_ylabel("freq")
bins = value.sel({"rank_hist_type": "bins", self.ahead_dim: ah.values})
freq = value.sel({"rank_hist_type": "freq", self.ahead_dim: ah.values})
relfreq = freq/freq.sum()
ax.bar(bins, relfreq, align='center', color="gray")
ax.hlines(1. / len(bins), 0, len(bins), color="black", linestyle="dashed")
ax.set_xlabel(r"Verification Rank ($n_{ens}+1$)")
ax.set_ylabel("Relative Frequency")
plt.title(f"{ah.values}")
pdf_pages.savefig()
plt.close('all')
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment