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)
# simulated runoff
times, vols = postprocessor.get_reach_timeseries('ROVOL', comid,
tstep = 'hourly')
vols = postprocessor.aggregate_hourly_daily(vols)
tot_sim = (sum([v * conv / area for v in vols]) /
(e - s).days * 365.25)
pct = (tot_sim - tot_obs) / tot_obs
t = ('total observed runoff: {:.0f} mm/yr\n'.format(tot_obs) +
'total simulated runoff: {:.0f} mm/yr\n'.format(tot_sim) +
'percent error: {:.1%}'.format(pct))
ax6.text(0.995, 0.02, t, transform = ax6.transAxes,
ha = 'right', va = 'bottom', color = 'black', size = ticksize)
# 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)
# simulated runoff
times, vols = postprocessor.get_reach_timeseries('ROVOL',
comid,
tstep='hourly')
vols = postprocessor.aggregate_hourly_daily(vols)
tot_sim = (sum([v * conv / area for v in vols]) / (e - s).days * 365.25)
pct = (tot_sim - tot_obs) / tot_obs
t = ('total observed runoff: {:.0f} mm/yr\n'.format(tot_obs) +
'total simulated runoff: {:.0f} mm/yr\n'.format(tot_sim) +
'percent error: {:.1%}'.format(pct))
ax6.text(0.995,
0.02,
t,
