for tick in ax.yaxis.get_major_ticks():
tick.label.set_fontsize(16)
fig = plt.savefig('%s.flow.season.png' %cfg['OUTPUT']['output_plot_basename'], format='png', dpi=dpi)
#============== plot flow duration curve (based on weekly data) ===============#
s = list_s_to_plot_weekly[0]
fig = my_functions.plot_duration_curve(\
list_s_data=list_s_to_plot_weekly_noNAN, \
list_style=list_plot_style, \
list_label=list_plot_label, \
figsize=(10,10), xlog=False, ylog=True, \
xlim=None, ylim=None, \
xlabel='Exceedence', ylabel='Flow (thousand cfs)', \
title='{}, WY {}-{}\n'.format(cfg['PLOT_OPTIONS']['plot_title'], \
plot_start_date.year+1, \
plot_end_date.year), \
fontsize=18, legend_loc='upper right', \
add_info_text=True, model_info=model_info, \
stats='Flow duration curve based on weekly data', show=False)
ax = plt.gca()
for tick in ax.xaxis.get_major_ticks():
tick.label.set_fontsize(16)
for tick in ax.yaxis.get_major_ticks():
tick.label.set_fontsize(16)
fig = plt.savefig('%s.flow_duration_weekly.png' %cfg['OUTPUT']['output_plot_basename'], format='png', dpi=dpi)
#-------------------------------------------
df_to_plot = pd.DataFrame(s_TVA_to_plot.values, index=s_TVA_to_plot.index, columns=['tva'])
df_to_plot['usgs'] = s_usgs_to_plot
df_to_plot = df_to_plot.dropna()
s_usgs_to_plot = df_to_plot['usgs']
s_TVA_to_plot = df_to_plot['tva']
#-------------------------------------------
#=== Plot flow duration curves (daily data) ===#
fig = my_functions.plot_duration_curve(\
list_s_data=[s_TVA_to_plot, s_usgs_to_plot], \
list_style=['k-', 'b-'], \
list_label=['TVA pass-through daily', 'USGS daily'], \
figsize=(10,10), xlog=False, ylog=True, \
xlim=None, ylim=None, \
xlabel='Exceedence', ylabel='Flow (thousand cfs)', \
title='Dam {}, daily flow duration, WY {}-{}'.format(dam_number, \
plot_start_date.year+1, \
plot_end_date.year), \
fontsize=18, legend_loc='upper right', \
add_info_text=True, model_info='Compare USGS flow data with TVA pass-through', \
stats='Flow duration curve based on daily data', show=False)
ax = plt.gca()
for tick in ax.xaxis.get_major_ticks():
tick.label.set_fontsize(16)
for tick in ax.yaxis.get_major_ticks():
tick.label.set_fontsize(16)
plt.savefig(os.path.join(out_plot_dir, \
'dam{}.flow_duration_daily.png'.format(dam_number)), \
fig = plt.savefig('%s.flow.season.png' % cfg['OUTPUT']['output_plot_basename'],
format='png',
dpi=dpi)
#============== plot flow duration curve (based on weekly data) ===============#
s = list_s_to_plot_weekly[0]
fig = my_functions.plot_duration_curve(\
list_s_data=list_s_to_plot_weekly_noNAN, \
list_style=list_plot_style, \
list_label=list_plot_label, \
figsize=(10,10), xlog=False, ylog=True, \
xlim=None, ylim=None, \
xlabel='Exceedence', ylabel='Flow (thousand cfs)', \
title='{}, WY {}-{}\n'.format(cfg['PLOT_OPTIONS']['plot_title'], \
plot_start_date.year+1, \
plot_end_date.year), \
fontsize=18, legend_loc='upper right', \
add_info_text=True, model_info=model_info, \
stats='Flow duration curve based on weekly data', show=False)
ax = plt.gca()
for tick in ax.xaxis.get_major_ticks():
tick.label.set_fontsize(16)
for tick in ax.yaxis.get_major_ticks():
tick.label.set_fontsize(16)
fig = plt.savefig('%s.flow_duration_weekly.png' %
cfg['OUTPUT']['output_plot_basename'],