ax6.set_ylim([0.08, 2000])
ax6.xaxis.set_major_formatter(dates.DateFormatter('%Y'))
ax6.xaxis.set_major_locator(dates.YearLocator(5))
ax6.plot(stimes, d_sflow, linestyle = '-', lw = 0.5, c = 'red',
label = 'simulated')
ax6.plot(otimes, d_oflow, marker = '+', lw = 0, c = 'red',
markersize = 4, label = 'observed')
ax6.set_title('Hydrograph', fontsize = titlesize)
ax6.legend(loc = 2, fontsize = ticksize, frameon = False)
# get runoff area (km2)
comids = postprocessor.get_upstream_comids(comid)
area = sum(postprocessor.get_subbasin_areas(comids))
# convert to mm
conv = 1000
# start and end dates
s, e = process_dates
# get the observed flows
t, oflow = postprocessor.get_obs_flow(tstep = 'daily')
tot_obs = (sum([f * 86400 / conv / area for f in oflow]) /
(e - s).days * 365.25)
ax6.plot(stimes, d_sflow, linestyle='-', lw=0.5, c='red', label='simulated')
ax6.plot(otimes,
d_oflow,
marker='+',
lw=0,
c='red',
markersize=4,
label='observed')
ax6.set_title('Hydrograph', fontsize=titlesize)
ax6.legend(loc=2, fontsize=ticksize, frameon=False)
# get runoff area (km2)
comids = postprocessor.get_upstream_comids(comid)
area = sum(postprocessor.get_subbasin_areas(comids))
# convert to mm
conv = 1000
# start and end dates
s, e = process_dates
# get the observed flows
t, oflow = postprocessor.get_obs_flow(tstep='daily')
tot_obs = (sum([f * 86400 / conv / area
for f in oflow]) / (e - s).days * 365.25)
