# ncol=2, mode="expand", borderaxespad=0.)
#hfmt = dates.DateFormatter('%H:%M')
#ax.xaxis.set_major_formatter(hfmt)
y_formatter = mpl.ticker.ScalarFormatter(useOffset=False)
ax.yaxis.set_major_formatter(y_formatter)
ax.grid(True)
f.suptitle(r"Schnellabschaltung KKW Gundremmingen 25.03.2015, $\Delta P_a = -1290MW$")
f.autofmt_xdate()
plt.savefig("images/gundremmingen-frequenzverlauf.png", bbox_inches='tight')
plt.clf()
print "Drawing momentum graph"
#fig, ax = plt.subplots(figsize=(16, 9), dpi=75)
fig, ax = plt.subplots(2, figsize=(16, 9), dpi=75)
grundremmingen_momentum_df = datatool.resample_add_freq_momentum(grundremmingen, 19000)
lower_momentum_limit = np.min(grundremmingen_momentum_df.momentum)
upper_momentum_limit = np.max(grundremmingen_momentum_df.momentum)
ax[0].set_xlabel("Zeit [Sekunden seit Mitternacht UTC]")
ax[0].set_ylabel("Frequenz [Hz]")
ax[0].set_xlim((np.min(grundremmingen_momentum_df.s_since_midnight.astype(int)),
np.max(grundremmingen_momentum_df.s_since_midnight.astype(int))))
ax[0].plot(grundremmingen_momentum_df.s_since_midnight.astype(int),
grundremmingen_momentum_df.freq, 'r', label=r"Frequenzmessung")
ax[0].plot(grundremmingen_momentum_df.s_since_midnight.astype(int),
grundremmingen_momentum_df.freq_sg, 'b', label=r"Savitzky-Golay (7, 2)")
ax[0].legend(loc="lower right", fontsize="small")
#hfmt = dates.DateFormatter('%H:%M')
#ax[1].xaxis.set_major_formatter(hfmt)
y_formatter = mpl.ticker.ScalarFormatter(useOffset=False)
#ax.legend(bbox_to_anchor=(0., 1.02, 1., .102), loc=3,
# ncol=2, mode="expand", borderaxespad=0.)
#hfmt = dates.DateFormatter('%H:%M')
#ax.xaxis.set_major_formatter(hfmt)
y_formatter = mpl.ticker.ScalarFormatter(useOffset=False)
ax.yaxis.set_major_formatter(y_formatter)
ax.grid(True)
f.suptitle(r"Schnellabschaltung KKW M\"uhleberg 06.07.2015, $\Delta P_a = -1290MW$")
f.autofmt_xdate()
plt.savefig("images/muehleberg-frequenzverlauf.png")#, bbox_inches='tight')
plt.clf()
print "Drawing momentum graph"
f, ax = plt.subplots(2)
muehleberg_momentum_df = datatool.resample_add_freq_momentum(muehleberg, 19000)
lower_momentum_limit = np.min(muehleberg_momentum_df.momentum)
upper_momentum_limit = np.max(muehleberg_momentum_df.momentum)
ax[0].set_xlabel("Zeit [Sekunden seit Mitternacht UTC]")
ax[0].set_ylabel("Frequenz [Hz]")
ax[0].set_xlim((np.min(muehleberg_momentum_df.s_since_midnight.astype(int)),
np.max(muehleberg_momentum_df.s_since_midnight.astype(int))))
ax[0].plot(muehleberg_momentum_df.s_since_midnight.astype(int),
muehleberg_momentum_df.freq, 'r', label=r"Frequenzmessung")
ax[0].plot(muehleberg_momentum_df.s_since_midnight.astype(int),
muehleberg_momentum_df.freq_sg, 'b', label=r"Savitzky-Golay (7, 2)")
ax[0].legend(loc="lower right", fontsize="small")
#hfmt = dates.DateFormatter('%H:%M')
#ax[1].xaxis.set_major_formatter(hfmt)
y_formatter = mpl.ticker.ScalarFormatter(useOffset=False)
