NZ = ea.panel_beater(storage.sum(axis=1), inflows.sum(axis=1), 365)
Taupo = ea.panel_beater(storage["Taupo"], inflows["Taupo"], 365)
Tekapo = ea.panel_beater(storage["Tekapo"], inflows["Tekapo"], 365)
Pukaki = ea.panel_beater(storage["Pukaki"], inflows["Pukaki"], 365)
Hawea = ea.panel_beater(storage["Hawea"], inflows["Hawea"], 365)
TeAnauManapouri = ea.panel_beater(
storage["TeAnau"] + storage["Manapouri"], inflows["TeAnau"] + inflows["Manapouri"], 365
)
NZ_HRC_SO, SI_HRC_SO = ea.get_SO_HRC(so_hrc)
NZ_HRC, SI_HRC = ea.get_EA_HRC(ea_hrc)
SI_HRC2 = SI_HRC_SO["6/2013":].append(SI_HRC["12/2012":"5/2013"])
NZ_HRC2 = NZ_HRC_SO["6/2013":].append(NZ_HRC["12/2012":"5/2013"])
history_repeats = ea.calc_mean_percentile(storage.sum(axis=1), 80, future=True)
history_repeats = history_repeats["12/" + str(dt.datetime.now().year - 1) : "2/1/" + str(dt.datetime.now().year + 1)]
path = """/home/dave/python/reporting/"""
os.chdir(path)
p = ea.hrc_plot(
1,
SI_HRC2["12/" + str(dt.datetime.now().year - 1) : "2/1/" + str(dt.datetime.now().year + 1)],
history_repeats["Actual"],
history_repeats["mean"],
"figures/hrc.pdf",
)
p = ea.hrc_plot(
1,
