@@ -121,6 +121,28 @@ mKGE_score = Streamfall.mKGE(dam_obs[366:end], dam_sim[366:end])
121121@info " Scores: " rmse_score nse_score mKGE_score
122122
123123f2 = quickplot (dam_obs, dam_sim, climate, " Modelled - 406000" , false ; burn_in= 366 )
124+
125+ temporal_cross_section (sim_dates, calib_data[:, " 406000" ], sn2[3 ]. level)
124126```
125127
128+ The last two lines produces the plots below
129+
126130![ ] ( ../assets/calibrated_example.png )
131+
132+ The ` quickplot() ` function creates the figure displayed above which shows dam levels on the
133+ left (observed and modelled) with a [ Q-Q plot] ( https://en.wikipedia.org/wiki/Q%E2%80%93Q_plot )
134+ on the right.
135+
136+ ![ ] ( ../assets/temporal_xsection_historic_calibrated.png )
137+
138+ The above shows a "cross-section" of model predictions for each month-day across simulation
139+ time. It is useful to gain an understanding on when models may underperform and give a
140+ sense of a models predictive uncertainty. The units of the y-axis are the same as for the
141+ node (in this case, meters).
142+
143+ Ideally, the median error would be a straight line and the confidence intervals would
144+ be as thin and consistent as possible for all month-days.
145+
146+ Here, we see that while performance is generally good (mean of Median Error is near zero),
147+ the model can under-estimate dam levels in late-April to May and displays a tendency to
148+ over-estimate dam levels between January and June, relative to other times.
0 commit comments