ax6.set_ylabel('Daily Flow (m\u00B3/s)') 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]) /
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