From 0eeaafec43c6fccf197f2b7e9eb3827a79b3e3b7 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Tobias=20M=C3=BCller?= Date: Fri, 8 Mar 2024 00:52:01 -0800 Subject: [PATCH] Pdf fixes and other recent issues (#73) * min width in tour template * slide changes Jon requested * updated sample data * removed header for sample data, so it can be imported quickly * fixed tooltip for rating error rounding incorrectly * added bold line in error plot to clearly show optimization target * made compare curve input smaller * added notes and more info by Jons request * hope to fix problem of 0 Q from autofit. needs testing * increased sig figs on tooltip * auto formatting and fixing PDF plots * unabbridged equation in pdf --- HydRA Test Curves_20240115/RC_LCR2122.csv | 20 -- HydRA Test Curves_20240115/RC_Lizard.csv | 6 - HydRA Test Curves_20240115/RC_cook1.csv | 20 -- HydRA Test Curves_20240115/RC_cook2.csv | 10 - HydRA Test Curves_20240115/RC_grand.csv | 24 -- HydRA Test Curves_20240115/RC_random.csv | 22 -- HydRA Test Curves_20240115/RC_shelly.csv | 11 - .../rctool/develop/develop_table_rating.html | 42 ++- .../rctool/rctool/develop/rctool_develop.html | 2 + .../rctool/develop/rctool_develop_js.html | 39 +- .../rctool/rctool/export/rctool_export.html | 3 +- .../rctool/export/rctool_export_pdf.html | 4 +- .../rctool/rctool/tour/rctool_tour_intro.html | 33 +- frontend/rctool/views.py | 340 ++++++++++++------ frontend/static/images/RMSE.JPG | Bin 21493 -> 25698 bytes frontend/static/sample_data/sample_data.csv | 35 +- 16 files changed, 336 insertions(+), 275 deletions(-) delete mode 100644 HydRA Test Curves_20240115/RC_LCR2122.csv delete mode 100644 HydRA Test Curves_20240115/RC_Lizard.csv delete mode 100644 HydRA Test Curves_20240115/RC_cook1.csv delete mode 100644 HydRA Test Curves_20240115/RC_cook2.csv delete mode 100644 HydRA Test Curves_20240115/RC_grand.csv delete mode 100644 HydRA Test Curves_20240115/RC_random.csv delete mode 100644 HydRA Test Curves_20240115/RC_shelly.csv diff --git a/HydRA Test Curves_20240115/RC_LCR2122.csv b/HydRA Test Curves_20240115/RC_LCR2122.csv deleted file mode 100644 index 2917e17..0000000 --- a/HydRA Test Curves_20240115/RC_LCR2122.csv +++ /dev/null @@ -1,20 +0,0 @@ -datetime,discharge,stage,uncertainty,comments -2021-11-24 10:27,0.9085,0.492,8,"23 panels, max 9.2%" -2021-12-07 11:15,1.751439049,0.501,3.11, -2021-12-16 10:59,1.1704,0.457,8,"22 panels, max 7.4%" -2022-01-14 11:05,4.891958877,0.73,2.67, -2022-02-07 11:38,1.127406583,0.459,3.27, -2022-03-23 10:16,2.319489087,0.55,2.65, -2022-03-30 14:45,0.843988357,0.515,2.94, -2022-04-22 13:56,0.569,0.397,,"20 panels, max 13.8%" -2022-05-03 17:29,0.465,0.377,,"20 panels, max 13.6%" -2022-05-18 14:37,0.981925337,0.437,3.23, -2022-07-19 9:53,0.131,0.277,,"20 panels, max 14.3%" -2022-07-19 10:09,0.1472,0.278,3, -2022-08-04 9:54,0.152,0.283,,"20 panels, max 11.9%" -2022-08-04 12:30,0.183566024,0.286,3.17, -2022-08-31 11:51,0.111,0.263,,"20 panels, max 9.2%" -2022-09-09 10:02,0.101913626,0.261,2.59, -2022-10-20 8:40,0.0813,0.243,2.9, -2022-12-12 10:55,1.672120912,0.4875,3.13, -2022-12-12 11:02,1.631,0.485,2.7, diff --git a/HydRA Test Curves_20240115/RC_Lizard.csv b/HydRA Test Curves_20240115/RC_Lizard.csv deleted file mode 100644 index 4a2e620..0000000 --- a/HydRA Test Curves_20240115/RC_Lizard.csv +++ /dev/null @@ -1,6 +0,0 @@ -datetime,discharge,stage,uncertainty,comments -2021-08-19 11:45,0.364,0.302,0.05,"old control, do not use MMT" -2022-07-20 8:54,0.612,0.225,0.05,suggested PZF -0.150 -2022-07-27 15:27,0.482,0.193,0.05, -2022-09-07 11:00,0.383,0.125,0.05, -2022-10-18 11:33,0.201,0.097,0.05, diff --git a/HydRA Test Curves_20240115/RC_cook1.csv b/HydRA Test Curves_20240115/RC_cook1.csv deleted file mode 100644 index 49fc7fe..0000000 --- a/HydRA Test Curves_20240115/RC_cook1.csv +++ /dev/null @@ -1,20 +0,0 @@ -datetime,discharge,stage,uncertainty,comments -2023-07-25 10:50:08,0.024789663,9.0635,0.099,approx PZF 9.0 -2023-05-08 13:23:17,0.468793188,9.193,0.031, -2023-04-18 12:45:02,1.0337748,9.263,0.042, -2023-02-06 12:00:00,1.080884219,9.27,0.032, -2022-08-18 10:15:00,0.028,9.037,0.049, -2022-07-19 10:11:16,0.093839378,9.065,0.03, -2022-07-19 10:45:00,0.088,9.065,0.033, -2022-04-13 13:00:00,0.557,9.177,0.03, -2022-03-09 11:00:00,0.406728145,9.16,0.03, -,,,, -,,,, -,,,, -,,,, -,,,, -,,,, -,,,, -,,,, -,,,, -,,,, diff --git a/HydRA Test Curves_20240115/RC_cook2.csv b/HydRA Test Curves_20240115/RC_cook2.csv deleted file mode 100644 index 9525aee..0000000 --- a/HydRA Test Curves_20240115/RC_cook2.csv +++ /dev/null @@ -1,10 +0,0 @@ -datetime,discharge,stage,stage_raw,uncertainty,comments -2023-07-25 10:50,0.024789663,0.3085,9.0635,0.099,approx PZF .245 -2023-05-08 13:23,0.468793188,0.438,9.193,0.031, -2023-04-18 12:45,1.0337748,0.508,9.263,0.042, -2023-02-06 12:00,1.080884219,0.515,9.27,0.032, -2022-08-18 10:15,0.028,0.282,9.037,0.049, -2022-07-19 10:11,0.093839378,0.31,9.065,0.03, -2022-07-19 10:45,0.088,0.31,9.065,0.033, -2022-04-13 13:00,0.557,0.422,9.177,0.03, -2022-03-09 11:00,0.406728145,0.405,9.16,0.03, diff --git a/HydRA Test Curves_20240115/RC_grand.csv b/HydRA Test Curves_20240115/RC_grand.csv deleted file mode 100644 index 5caaf45..0000000 --- a/HydRA Test Curves_20240115/RC_grand.csv +++ /dev/null @@ -1,24 +0,0 @@ -datetime,stage,discharge,Uncertainty,Comments -2023-07-26 12:57:30,0.186,0.002772154,,"PZF 0.018 , BP @ 0.395M" -2023-05-26 13:02:30,0.198,0.005309356,,minor debris -2023-04-19 13:51:15,0.2965,0.039120759,6.8,section -2023-03-03 15:42:40,0.442,0.161655436,3.7,submerged -2023-02-12 15:15:00,0.368,0.066844777,3.683071598,section -2022-11-04 13:00:00,0.24,0.016708148,8, -2022-09-19 16:00:00,0.188,0.00292457,8, -2022-08-17 10:00:00,0.19,0.003312321,5, -2022-07-18 15:00:00,0.21,0.004,, -2021-08-17 09:00:00,0.285,0.002114971,,debris -2019-09-19 13:45:00,0.23,0.01125648,10, -2019-07-29 11:15:00,0.196,0.003352529,10, -2019-06-10 10:40:00,0.209,0.00646972,8, -2018-06-05 12:00:00,0.218,0.007823333,8, -2017-11-27 13:30:00,0.54,0.4602,3.6, -2017-07-12 13:30:00,0.205,0.004951,8, -2017-06-13 16:00:00,0.215,0.007517333,8, -2014-12-19 11:30:00,0.498,0.2883,,submerged -2014-02-17 13:00:00,0.658,1.063,8,submerged -2013-10-03 12:00:00,0.318,0.0388,5, -2013-06-25 09:30:00,0.247,0.0162,8, -2012-10-11 15:00:00,0.203,0.00382,5, -2012-09-12 13:00:00,0.198,0.002749,8, diff --git a/HydRA Test Curves_20240115/RC_random.csv b/HydRA Test Curves_20240115/RC_random.csv deleted file mode 100644 index 8d8e49e..0000000 --- a/HydRA Test Curves_20240115/RC_random.csv +++ /dev/null @@ -1,22 +0,0 @@ -date,time,DateTime,discharge,stage,uncertainty,comments -2022-03-25,12:00:00,2022-03-25 12:00,2.796,0.759,0.004, -2021-06-25,0:48:23,2021-06-25 0:48,2.827,0.774,0.008, -2020-11-04,14:45:54,2020-11-04 14:45,3.816,2.792,0.060, -2020-09-23,6:43:37,2020-09-23 6:43,3.530,1.838,0.095, -2020-08-22,3:12:13,2020-08-22 3:12,4.587,5.577,0.041, -2020-03-24,11:10:33,2020-03-24 11:10,4.328,3.894,0.001, -2019-05-30,3:34:14,2019-05-30 3:34,3.805,2.830,0.098, -2018-08-31,12:10:31,2018-08-31 12:10,3.588,1.945,0.022, -2017-09-15,2:11:41,2017-09-15 2:11,3.808,2.590,0.023, -2017-06-13,5:36:18,2017-06-13 5:36,4.358,4.321,0.004, -2016-06-18,14:06:14,2016-06-18 14:06,5.582,8.625,0.029, -2016-02-14,8:54:33,2016-02-14 8:54,5.537,9.343,0.021, -2015-06-18,1:21:39,2015-06-18 1:21,5.468,8.690,0.063, -2015-02-16,13:05:51,2015-02-16 13:05,2.377,0.349,0.047, -2014-09-30,0:54:20,2014-09-30 0:54,5.054,9.185,0.080, -2014-03-14,22:27:11,2014-03-14 22:27,5.047,9.827,0.001, -2013-03-26,21:35:10,2013-03-26 21:35,4.690,6.033,0.064, -2013-03-04,0:46:13,2013-03-04 0:46,4.009,3.816,0.046, -2012-05-21,21:36:24,2012-05-21 21:36,5.416,9.181,0.018, -,,,,,, -,,,,,, diff --git a/HydRA Test Curves_20240115/RC_shelly.csv b/HydRA Test Curves_20240115/RC_shelly.csv deleted file mode 100644 index 64f1f98..0000000 --- a/HydRA Test Curves_20240115/RC_shelly.csv +++ /dev/null @@ -1,11 +0,0 @@ -datetime,discharge,stage,uncertainty,comments -2021-11-22 10:00,0.13,0.33,0.05,"21 panels, some snow on ground" -2021-04-06 14:00,0.01,0.24,0.05,19 panels -2021-03-08 9:00,0.05,0.28,0.03,"35 panels, recent high flow event in December" -2021-02-24 13:00,0.03,0.27,0.05,"21 panels, rain overnight" -2021-01-28 14:00,0.12,0.32,0.04,"21 panels, good mmt" -2021-01-11 10:00,0.18,0.39,0.04,"21 panels, recent high flows after consistent rain. High turbidity" -2020-12-30 11:30,0.36,0.45,0.03,"24 panels, recent high flows and heavy rains" -2020-11-24 13:30,0.03,0.26,0.04,21 panels. Rain overnight -2021-12-13 8:00,0.29,0.4,0.04,20 panels -2020-01-12 17:00,0.67,0.55,0.03,"23 panels, recent rain on snow event, creek rising during mmt" diff --git a/frontend/rctool/templates/rctool/rctool/develop/develop_table_rating.html b/frontend/rctool/templates/rctool/rctool/develop/develop_table_rating.html index eebdab3..6640bd6 100644 --- a/frontend/rctool/templates/rctool/rctool/develop/develop_table_rating.html +++ b/frontend/rctool/templates/rctool/rctool/develop/develop_table_rating.html @@ -9,8 +9,8 @@ coefficient exponent offset - start point - end point + start point
(stage in m) + end point
(stage in m) RMSE

compare

@@ -40,32 +40,52 @@ {% endif %} - + - + - + - + - + - + -

+

- +
diff --git a/frontend/rctool/templates/rctool/rctool/develop/rctool_develop.html b/frontend/rctool/templates/rctool/rctool/develop/rctool_develop.html index ac85361..74f99b2 100644 --- a/frontend/rctool/templates/rctool/rctool/develop/rctool_develop.html +++ b/frontend/rctool/templates/rctool/rctool/develop/rctool_develop.html @@ -207,7 +207,9 @@

|
", onEnd: function (tour) { sessionStorage.setItem("tourStatus", false); $('#tour_request_status_id').val("1") diff --git a/frontend/rctool/templates/rctool/rctool/develop/rctool_develop_js.html b/frontend/rctool/templates/rctool/rctool/develop/rctool_develop_js.html index f095bea..ea5f282 100644 --- a/frontend/rctool/templates/rctool/rctool/develop/rctool_develop_js.html +++ b/frontend/rctool/templates/rctool/rctool/develop/rctool_develop_js.html @@ -215,6 +215,13 @@ rcDict = {'data':rcData, 'parameters':rcParam}; // 3.5 Update HTML equation parameters and user input fields + // make sure all bounds for Q are at minimum 0.001: + if (rcParam[0].seg_bounds[0][1] < 0.001) { + rcParam[0].seg_bounds[0][1] = 0.001; + } + if (rcParam[0].seg_bounds[1][1] < 0.001) { + rcParam[0].seg_bounds[1][1] = 0.001; + } var segNumStr = String(idx + 1); document.getElementById('segment' + segNumStr + '-const').innerHTML = C.toFixed(2); document.getElementById('segment' + segNumStr + '-exp').innerHTML = slope.toFixed(2); @@ -501,7 +508,6 @@ ----------------------------------------------------------------------------------------------------------------------------------------------*/ function setEndpoints(segEP, pointEP, coordEP, valueEP) { - // 1. COLLECT AND INITIATE PARAMETERS if (rcParam.length == 2) { var newSegBounds = [rcParam[0].seg_bounds, rcParam[1].seg_bounds]; @@ -1289,12 +1295,12 @@ } } } - label.push(['H: ' + h_field.toFixed(2) + ' m']); - label.push(['Q: ' + q_field.toFixed(2) + ' m³/s']); + label.push(['H: ' + h_field.toFixed(3) + ' m']); + label.push(['Q: ' + q_field.toFixed(3) + ' m³/s']); label.push(['date: ' + fielddatacsv.datetime[idx]]); if (h_model !== null) { - label.push(['H shift: ' + (h_model - h_field).toFixed(2) + ' m']) - label.push(['H shift: ' + (h_model - h_field).toFixed(2) + ' m']) + label.push(['H shift: ' + (h_model - h_field).toFixed(3) + ' m']) + label.push(['H shift: ' + (h_model - h_field).toFixed(3) + ' m']) } //if (q_model !== null){ // label.push(['Q model: ' + q_model.toFixed(2) + ' m³/s']); @@ -1457,14 +1463,20 @@ display: true, text: 'Q rated vs Q field (% change)' }, - beginAtZero: false + beginAtZero: true, + // make zero tick thicker: + grid: { + color: function(context) { + return context.tick.value === 0 ? 'black' : 'rgba(0,0,0,0.1)'; + }, + }, }, y: { title: { display: true, text: 'Stage H (m)' }, - beginAtZero: false + beginAtZero: false, } }, @@ -1502,21 +1514,17 @@ let q_model = calc_q_model(h_field, rcParam[seg].const, rcParam[seg].exp, rcParam[seg].offset) // r_percent = 100 * (q_field - q_model) / q_field let q_field = - 100 * q_model / (r_percent - 100) - let q_field_rounded = Math.round(q_field * 100) / 100 + let q_field_rounded = Math.round(q_field * 1000) / 1000 let h_model = calc_h_model(q_field, rcParam[seg].const, rcParam[seg].exp, rcParam[seg].offset) let r_volumetric = r_percent/100 * q_field - for (var idx in fielddatacsv.stage) { if (fielddatacsv.stage[idx] == h_field && fielddatacsv.discharge[idx] == q_field_rounded) { - label.push(['H: ' + h_field.toFixed(2) + ' m']); - label.push(['Q: ' + q_field.toFixed(2) + ' m³/s']); + label.push(['H: ' + h_field.toFixed(3) + ' m']); + label.push(['Q: ' + q_field.toFixed(3) + ' m³/s']); label.push(['date: ' + fielddatacsv.datetime[idx]]); if (h_model !== null) { - label.push(['H shift: ' + (h_model - h_field).toFixed(2) + ' m']) + label.push(['H shift: ' + (h_model - h_field).toFixed(3) + ' m']) } - //if (q_model !== null){ - // label.push(['Q model: ' + q_model.toFixed(2) + ' m³/s']); - //} if (r_percent !== null){ if (r_percent > 0){ // show + sign for positive values @@ -1527,6 +1535,7 @@ label.push(['R volume:' + r_volumetric.toFixed(2) + ' m3/s']); } } + break; } } } diff --git a/frontend/rctool/templates/rctool/rctool/export/rctool_export.html b/frontend/rctool/templates/rctool/rctool/export/rctool_export.html index 5017638..72da451 100644 --- a/frontend/rctool/templates/rctool/rctool/export/rctool_export.html +++ b/frontend/rctool/templates/rctool/rctool/export/rctool_export.html @@ -18,7 +18,8 @@
Export
-
+
+

Note: The rating curve generated by HydRA produces an equation ONLY. All input (stage) and output (discharge) datasets including any discrete measurements of stage and discharge must be reviewed and graded according to the provincial hydrometric RISC standards.

{% csrf_token %} diff --git a/frontend/rctool/templates/rctool/rctool/export/rctool_export_pdf.html b/frontend/rctool/templates/rctool/rctool/export/rctool_export_pdf.html index 9b46e60..898b4db 100644 --- a/frontend/rctool/templates/rctool/rctool/export/rctool_export_pdf.html +++ b/frontend/rctool/templates/rctool/rctool/export/rctool_export_pdf.html @@ -109,7 +109,7 @@
Model Parameters
{{ rc.parameters.0.label }} - Q = {{ rc.parameters.0.const|floatformat:3 }} ( H {{ offsets_val.0 }} ){{ rc.parameters.0.exp|floatformat:3 }} + Q = {{ rc.parameters.0.const }} ( H {{ offsets_val.0 }} ){{ rc.parameters.0.exp }} {{ rc.parameters.0.seg_bounds.0.0|floatformat:3 }} , {{ rc.parameters.0.seg_bounds.1.0|floatformat:3}} {{ rc.parameters.0.rmse|floatformat:3 }} @@ -117,7 +117,7 @@
Model Parameters
{% if rc.parameters.1 %} {{ rc.parameters.1.label }} - Q = {{ rc.parameters.1.const|floatformat:3 }}( H {{ offsets_val.1 }}){{ rc.parameters.1.exp|floatformat:3 }} + Q = {{ rc.parameters.1.const }}( H {{ offsets_val.1 }}){{ rc.parameters.1.exp }} {{ rc.parameters.1.seg_bounds.0.0|floatformat:3 }} , {{ rc.parameters.1.seg_bounds.1.0|floatformat:3 }} {{ rc.parameters.1.rmse|floatformat:3 }} diff --git a/frontend/rctool/templates/rctool/rctool/tour/rctool_tour_intro.html b/frontend/rctool/templates/rctool/rctool/tour/rctool_tour_intro.html index 8be589b..34ab337 100644 --- a/frontend/rctool/templates/rctool/rctool/tour/rctool_tour_intro.html +++ b/frontend/rctool/templates/rctool/rctool/tour/rctool_tour_intro.html @@ -32,7 +32,8 @@
Let's Get Started
- Hydrometric field work involves regular visits to a hydrometric gauging station across a range of streamflow conditions. A hydrometric gauging station typically consists of a stage sensor and benchmarks or reference gauge from which water level is referenced. Each field visit should collect three basic observations: reference of stage, measurement of discharge, and control observation. These observations are required to convert stage into discharge through the development of a hydrometric rating curve. (Note: Safety should be considered when collecting these observations, do not attempt any unsafe field work!). The Hydrometric Rating Application (HydRA) is a tool for developing hydrometric rating curves based on observations of stage and discharge. HydRA allows the user to upload observations of stage and discharge and automatically fits a rating curve. Users can adjust the fit according to onsite observations. It then generates a rating equation that can be used for converting stage datasets into discharge. + Hydrometric field work involves regular visits to a hydrometric gauging station across a range of streamflow conditions. A hydrometric gauging station typically consists of a stage sensor and benchmarks or reference gauge from which water level is referenced. Each field visit should collect three basic observations: reference of stage, measurement of discharge, and control observation. These observations are required to convert stage into discharge through the development of a hydrometric rating curve. (Note: Safety should be considered when collecting these observations, do not attempt any unsafe field work!). +

The Hydrometric Rating Application (HydRA) is a tool for developing hydrometric rating curves based on observations of stage and discharge. HydRA allows the user to upload observations of stage and discharge and automatically fits a rating curve. Users can adjust the fit according to onsite observations. It then generates a rating equation that can be used for converting stage datasets into discharge.
@@ -50,22 +51,6 @@
Uploading Dat
-
The Rating Equation
@@ -112,7 +97,7 @@
Time for a Br
Performance
- Some common metrics used to assess the performance of hydrometric rating curves are included below: + The HydRa app offers the option to auto-fit rating curves by finding the best equation describing the uploaded field data. This is done by calculating common metrics used to assess the performance of hydrometric rating curves:
Root Mean Square Error
@@ -132,6 +117,18 @@
Performance
+
+
+
+
Data Review
+ HydRA is designed to ONLY produce a rating curve equation and summary report. The complete hydrometric data review process should include the review, correction, and grading of the input stage dataset and the output discharge dataset according to the 2018 BC Hydrometric RISC standards +

Applying RISC grades to the output discharge dataset may depend on the performance of the rating equation in predicting discharge or the range of discharge observed during the rating period. All hydrometric data should be reviewed and graded by a qualified professional according to the RISC standards. For more information, consult RISC Table 1-1. +
+ +
+
diff --git a/frontend/rctool/views.py b/frontend/rctool/views.py index ab9a91c..865ebd1 100644 --- a/frontend/rctool/views.py +++ b/frontend/rctool/views.py @@ -42,7 +42,9 @@ def rctool_import(request, tour_request_id=0): if request.method == "POST": if request.POST.get("tour_request_status_id"): - context["tour_request_status_id"] = request.POST.get("tour_request_status_id") + context["tour_request_status_id"] = request.POST.get( + "tour_request_status_id" + ) # If user is on the tour, upload test data if tour_request_id == 1: @@ -544,9 +546,9 @@ def rctool_develop_initialize(request): "headings": field_df_raw.columns.values, "data": field_df_raw.values.tolist(), } - context[ - "n_seg" - ] = 1 # this will be initialized as None and make the user enter a value, for now use 2 + context["n_seg"] = ( + 1 # this will be initialized as None and make the user enter a value, for now use 2 + ) # init_offsets = [None, None] # initialize at None, then optimize in autofit init_offsets = [0.0, 0.0] context["breakpoint1"] = None @@ -739,45 +741,100 @@ def rctool_export_initialize(request): def create_export_rc_img(field_data, rc_data): # Initialize and prepaire pallet = ["#80B7AB", "#CC6677", "#003466"] - df_field = field_data.copy() df_field_active = field_data[field_data["toggle_point"] == "checked"] df_field_inactive = field_data[field_data["toggle_point"] == "unchecked"] + + exponent = rc_data["parameters"][0]["exp"] + + log_base_n = lambda x: np.log(x) / np.log(exponent) + # create plot obj fig1, ax1 = plt.subplots(figsize=(11, 5), num=1) # scales etc - plt.yscale("log") - plt.xscale("log") + plt.yscale("function", functions=(log_base_n, np.exp)) + plt.xscale("function", functions=(log_base_n, np.exp)) ax1.grid(True, which="both") ax1.set_ylabel("Stage H (m)") ax1.set_xlabel("Discharge Q (m$^3$/s)") ax1.xaxis.set_major_formatter(ticker.FormatStrFormatter("%.3f")) - ax1.xaxis.set_major_locator(plt.MaxNLocator(10)) - ax1.yaxis.set_major_formatter(ticker.FormatStrFormatter("%.3f")) + + # set bounds to field data plus/minus 10% + x_min = df_field_active["discharge"].min() * 0.9 + x_max = df_field_active["discharge"].max() * 1.1 + y_min = df_field_active["stage"].min() * 0.9 + y_max = df_field_active["stage"].max() * 1.1 + + if x_min > rc_data["parameters"][0]["seg_bounds"][0][1]: + x_min = rc_data["parameters"][0]["seg_bounds"][0][1] * 0.9 + if x_max < rc_data["parameters"][0]["seg_bounds"][1][1]: + x_max = rc_data["parameters"][0]["seg_bounds"][1][1] * 1.1 + if y_min > rc_data["parameters"][0]["seg_bounds"][0][0]: + y_min = rc_data["parameters"][0]["seg_bounds"][0][0] * 0.9 + if y_max < rc_data["parameters"][0]["seg_bounds"][1][0]: + y_max = rc_data["parameters"][0]["seg_bounds"][1][0] * 1.1 + + # if there are two segments, check if the bounds are within the field data + if len(rc_data["parameters"]) > 1: + if x_min > rc_data["parameters"][1]["seg_bounds"][0][1]: + x_min = rc_data["parameters"][1]["seg_bounds"][0][1] * 0.9 + if x_max < rc_data["parameters"][1]["seg_bounds"][1][1]: + x_max = rc_data["parameters"][1]["seg_bounds"][1][1] * 1.1 + if y_min > rc_data["parameters"][1]["seg_bounds"][0][0]: + y_min = rc_data["parameters"][1]["seg_bounds"][0][0] * 0.9 + if y_max < rc_data["parameters"][1]["seg_bounds"][1][0]: + y_max = rc_data["parameters"][1]["seg_bounds"][1][0] * 1.1 + ax1.set_xlim(x_min, x_max) + ax1.set_ylim(y_min, y_max) + + # spread ticks evenly, calculate with log_base_n + ticks_x = np.array( + [ + exponent**i + for i in np.linspace( + math.floor(log_base_n(x_min)), + math.ceil(log_base_n(x_max)), + 10, + ) + ] + ) + # if there are too few y ticks, add more + ticks_y = np.array( + [ + exponent**i + for i in np.linspace( + math.floor(log_base_n(y_min)), + math.ceil(log_base_n(y_max)), + 10, + ) + ] + ) + ax1.xaxis.set_major_locator(ticker.FixedLocator(ticks_x)) + ax1.yaxis.set_major_locator(ticker.FixedLocator(ticks_y)) + ax1.plot( - df_field_active["stage"], df_field_active["discharge"], + df_field_active["stage"], "o", color=pallet[0], label="field data", ) - if len(df_field_inactive["stage"]) > 0: - ax1.plot( - df_field_inactive["stage"], - df_field_inactive["discharge"], - "o", - color="#661100", - edgecolors='#661100', - label="field data (inactive)", - facecolors='none', - ) + # add data points for inactive + ax1.plot( + df_field_inactive["discharge"], + df_field_inactive["stage"], + "o", + color="#661100", + label="field data (inactive)", + ) + i = 1 for segment_parameters in rc_data["parameters"]: seg_bounds = segment_parameters["seg_bounds"] y_seg = [seg_bounds[0][0], seg_bounds[1][0]] x_seg = [seg_bounds[0][1], seg_bounds[1][1]] ax1.plot( - y_seg, x_seg, + y_seg, ls="-", color=pallet[i], label=segment_parameters["label"], @@ -794,14 +851,21 @@ def create_export_rc_img(field_data, rc_data): def create_export_res_img(field_data, rc_data): - res_df = field_data.copy() + df_field_active = field_data[field_data["toggle_point"] == "checked"] # create plot obj fig2, ax2 = plt.subplots(figsize=(11, 5), num=2) ax2.grid(True, which="both") # scales etc ax2.set_ylabel("Stage H (m)") ax2.set_xlabel("Discharge Error (%)") - ax2.plot(res_df["Discharge Error (%)"], res_df["stage"], "o", color="#80B7AB") + ax2.plot( + df_field_active["Discharge Error (%)"], + df_field_active["stage"], + "o", + color="#80B7AB", + ) + # make 0 line thick + ax2.axvline(x=0, color="black", lw=0.5) tmpfile2 = io.BytesIO() fig2.savefig(tmpfile2) b64_2 = base64.b64encode(tmpfile2.getvalue()).decode() @@ -846,7 +910,7 @@ def export_calculate_discharge_error(field_data_output_df, rc_output_dict): loc=4, column="Discharge Error (%)", value=error_perc_dat ) - field_data_output_df["Discharge Error (%)"] = field_data_output_df[ + field_data_output_df["Discharge Error (%)"] = -field_data_output_df[ "Discharge Error (%)" ].round(decimals=2) return field_data_output_df @@ -854,25 +918,29 @@ def export_calculate_discharge_error(field_data_output_df, rc_output_dict): def rctool_export_output(request): context = {} - floatlist = lambda strlist : list(map(float, strlist.split(','))) + floatlist = lambda strlist: list(map(float, strlist.split(","))) if request.method == "POST": export_form = export_rc_data(request.POST) - field_data_output_json = request.POST.get('fielddatacsv-to-output') + field_data_output_json = request.POST.get("fielddatacsv-to-output") field_data_output_df = pd.read_json(field_data_output_json) field_data_output_dict = field_data_output_df.to_dict() - field_data_output_df['datetime'] = field_data_output_df['datetime'].apply(str) - field_data_output_df['stage'] = field_data_output_df['stage'].round(decimals=3) - field_data_output_df['discharge'] = field_data_output_df['discharge'].round(decimals=3) - field_data_output_df['uncertainty'] = field_data_output_df['uncertainty'].round(decimals=3) + field_data_output_df["datetime"] = field_data_output_df["datetime"].apply(str) + field_data_output_df["stage"] = field_data_output_df["stage"].round(decimals=3) + field_data_output_df["discharge"] = field_data_output_df["discharge"].round( + decimals=3 + ) + field_data_output_df["uncertainty"] = field_data_output_df["uncertainty"].round( + decimals=3 + ) field_data_values = field_data_output_df.values.tolist() - rc_output = request.POST.get('rc_output') + rc_output = request.POST.get("rc_output") rc_output_dict = ast.literal_eval(rc_output) - - rc_output_dict['data'].append(field_data_values) - rc_output_df = pd.DataFrame.from_dict(rc_output_dict, orient='index') + + rc_output_dict["data"].append(field_data_values) + rc_output_df = pd.DataFrame.from_dict(rc_output_dict, orient="index") rc_output_df = rc_output_df.transpose() if export_form.is_valid(): @@ -880,145 +948,209 @@ def rctool_export_output(request): fname = export_form.cleaned_data["export_filename"] if export_form.cleaned_data["export_station_name"]: stname = export_form.cleaned_data["export_station_name"] - context['stname'] = stname + context["stname"] = stname if export_form.cleaned_data["export_comments"]: comments = export_form.cleaned_data["export_comments"] - context['comments'] = comments + context["comments"] = comments if export_form.cleaned_data["export_date_applic_init"]: app_period_start = export_form.cleaned_data["export_date_applic_init"] - context['app_period_start'] = app_period_start + context["app_period_start"] = app_period_start if export_form.cleaned_data["export_date_applic_final"]: app_period_end = export_form.cleaned_data["export_date_applic_final"] - context['app_period_end'] = app_period_end - - context['current_time'] = datetime.utcnow().strftime("%d/%m/%Y %H:%M") - context['rc'] = rc_output_dict + context["app_period_end"] = app_period_end + + context["current_time"] = datetime.utcnow().strftime("%d/%m/%Y %H:%M") + context["rc"] = rc_output_dict # Round output parameters - for idx in range(len(rc_output_dict['parameters'])): - rc_output_dict['parameters'][idx]['const'] = round(rc_output_dict['parameters'][idx]['const'], 4) - rc_output_dict['parameters'][idx]['exp'] = round(rc_output_dict['parameters'][idx]['exp'], 4) + for idx in range(len(rc_output_dict["parameters"])): + rc_output_dict["parameters"][idx]["const"] = round( + rc_output_dict["parameters"][idx]["const"], 4 + ) + rc_output_dict["parameters"][idx]["exp"] = round( + rc_output_dict["parameters"][idx]["exp"], 4 + ) # Format offset values for output equation - offsets = [0,0] - offsets[0] = rc_output_dict['parameters'][0]['offset'] - if len(rc_output_dict['parameters']) > 1: - offsets[1] = rc_output_dict['parameters'][1]['offset'] - context['offsets'] = offsets + offsets = [0, 0] + offsets[0] = rc_output_dict["parameters"][0]["offset"] + if len(rc_output_dict["parameters"]) > 1: + offsets[1] = rc_output_dict["parameters"][1]["offset"] + context["offsets"] = offsets # Format equation according to + or - offset (ie. if user enters a - offset, the equation will show (H + offset) instead of (H -- offset)) offsets_val = offsets.copy() if offsets_val[0] < 0: - offsets_val[0] = ' + ' + str(abs(offsets_val[0])) - context['offsets_val'] = offsets_val + offsets_val[0] = " + " + str(abs(offsets_val[0])) + context["offsets_val"] = offsets_val else: - offsets_val[0] = ' - ' + str(abs(offsets_val[0])) - context['offsets_val'] = offsets_val + offsets_val[0] = " - " + str(abs(offsets_val[0])) + context["offsets_val"] = offsets_val if offsets_val[1] < 0: - offsets_val[1] = ' + ' + str(abs(offsets_val[1])) - context['offsets_val'] = offsets_val + offsets_val[1] = " + " + str(abs(offsets_val[1])) + context["offsets_val"] = offsets_val else: - offsets_val[1] = ' - ' + str(abs(offsets_val[1])) - context['offsets_val'] = offsets_val + offsets_val[1] = " - " + str(abs(offsets_val[1])) + context["offsets_val"] = offsets_val - - if ftype == 'session settings': + if ftype == "session settings": try: # output HTTP response response = HttpResponse( - content_type='text/csv', - headers={'Content-Disposition': 'attachment; filename="{0}.dat"'.format(fname)}, + content_type="text/csv", + headers={ + "Content-Disposition": 'attachment; filename="{0}.dat"'.format( + fname + ) + }, ) # write to csv rc_output_df.to_csv(path_or_buf=response) return response except Exception as e: - messages.error(request,"Unable to process request. "+repr(e)) - - elif ftype == 'session results (csv)': + messages.error(request, "Unable to process request. " + repr(e)) + + elif ftype == "session results (csv)": export_calculate_discharge_error(field_data_output_df, rc_output_dict) - field_table = {'headings': field_data_output_df.columns.values, 'data': field_data_output_df.values.tolist()} - - context['rc_img'] = create_export_rc_img(field_data_output_df, rc_output_dict) - context['res_img'] = create_export_res_img(field_data_output_df, rc_output_dict) + field_table = { + "headings": field_data_output_df.columns.values, + "data": field_data_output_df.values.tolist(), + } + + context["rc_img"] = create_export_rc_img( + field_data_output_df, rc_output_dict + ) + context["res_img"] = create_export_res_img( + field_data_output_df, rc_output_dict + ) # create paramters df for output param_to_df = [] - for segment_parameter in rc_output_dict['parameters']: - segment_parameter['stage initial'] = segment_parameter['seg_bounds'][0][0] - segment_parameter['stage final'] = segment_parameter['seg_bounds'][1][0] - segment_parameter.pop('seg_bounds') + for segment_parameter in rc_output_dict["parameters"]: + segment_parameter["stage initial"] = segment_parameter[ + "seg_bounds" + ][0][0] + segment_parameter["stage final"] = segment_parameter["seg_bounds"][ + 1 + ][0] + segment_parameter.pop("seg_bounds") param_to_df.append(segment_parameter) df_parameters = pd.DataFrame(param_to_df) - parameter_table = {'headings': df_parameters.columns.values, 'data': df_parameters.values.tolist()} + parameter_table = { + "headings": df_parameters.columns.values, + "data": df_parameters.values.tolist(), + } # prepaire as output try: # output HTTP response response = HttpResponse( - content_type='text/csv', - headers={'Content-Disposition': 'attachment; filename="{0}.csv"'.format(fname)}, + content_type="text/csv", + headers={ + "Content-Disposition": 'attachment; filename="{0}.csv"'.format( + fname + ) + }, ) # # write to csv # rc_output_df.to_csv(path_or_buf=response) # return response writer = csv.writer(response) - writer.writerow(['RATING CURVE OUTPUT SUMMARY']) - writer.writerow(['date created (utc): {}'.format(context['current_time'])]) + writer.writerow(["RATING CURVE OUTPUT SUMMARY"]) + writer.writerow( + ["date created (utc): {}".format(context["current_time"])] + ) if export_form.cleaned_data["export_station_name"]: - writer.writerow(['location: {}'.format(export_form.cleaned_data["export_station_name"])]) - if export_form.cleaned_data["export_date_applic_init"] and export_form.cleaned_data["export_date_applic_final"]: - writer.writerow(['period of applicability: {0} to {1}'.format(context['app_period_start'], context['app_period_end'])]) + writer.writerow( + [ + "location: {}".format( + export_form.cleaned_data["export_station_name"] + ) + ] + ) + if ( + export_form.cleaned_data["export_date_applic_init"] + and export_form.cleaned_data["export_date_applic_final"] + ): + writer.writerow( + [ + "period of applicability: {0} to {1}".format( + context["app_period_start"], + context["app_period_end"], + ) + ] + ) if export_form.cleaned_data["export_comments"]: - writer.writerow(['comments: {}'.format(export_form.cleaned_data["export_comments"])]) - writer.writerow('') - writer.writerow(['MODEL PARAMETERS']) - writer.writerow(parameter_table['headings']) - for row in parameter_table['data']: + writer.writerow( + [ + "comments: {}".format( + export_form.cleaned_data["export_comments"] + ) + ] + ) + writer.writerow("") + writer.writerow(["MODEL PARAMETERS"]) + writer.writerow(parameter_table["headings"]) + for row in parameter_table["data"]: writer.writerow(row) - writer.writerow('') - writer.writerow(['MODEL INPUTS']) - writer.writerow(field_table['headings']) - for row in field_table['data']: + writer.writerow("") + writer.writerow(["MODEL INPUTS"]) + writer.writerow(field_table["headings"]) + for row in field_table["data"]: writer.writerow(row) return response except Exception as e: - messages.error(request,"Unable to process request. "+repr(e)) + messages.error(request, "Unable to process request. " + repr(e)) else: try: # write to pdf option selected # calculate residuals for output figures and table - field_data_output_df = export_calculate_discharge_error(field_data_output_df, rc_output_dict) - context['field_table'] = {'headings': field_data_output_df.columns.values, 'data': field_data_output_df.values.tolist()} + field_data_output_df = export_calculate_discharge_error( + field_data_output_df, rc_output_dict + ) + context["field_table"] = { + "headings": field_data_output_df.columns.values, + "data": field_data_output_df.values.tolist(), + } # create figure for pdf - context['rc_img'] = create_export_rc_img(field_data_output_df, rc_output_dict) - context['res_img'] = create_export_res_img(field_data_output_df, rc_output_dict) + context["rc_img"] = create_export_rc_img( + field_data_output_df, rc_output_dict + ) + context["res_img"] = create_export_res_img( + field_data_output_df, rc_output_dict + ) # prepaire and return output pdf - template = get_template('rctool/rctool/export/rctool_export_pdf.html') + template = get_template( + "rctool/rctool/export/rctool_export_pdf.html" + ) html = template.render(context) - pdf = render_to_pdf('rctool/rctool/export/rctool_export_pdf.html', context) - response = HttpResponse(pdf, content_type='application/pdf') - content = "inline; filename='%s'" %(fname) + pdf = render_to_pdf( + "rctool/rctool/export/rctool_export_pdf.html", context + ) + response = HttpResponse(pdf, content_type="application/pdf") + content = "inline; filename='%s'" % (fname) download = request.GET.get("download") if download: - content = "attachment; filename='%s.pdf'" %(fname) - response['Content-Disposition'] = content + content = "attachment; filename='%s.pdf'" % (fname) + response["Content-Disposition"] = content return response except Exception as e: - messages.error(request,"Unable to process request. "+repr(e)) + messages.error(request, "Unable to process request. " + repr(e)) - context['form'] = export_form - context['rc_output'] = rc_output + context["form"] = export_form + context["rc_output"] = rc_output return render(request, "rctool/rctool/export/rctool_export.html", context) else: export_form = export_rc_data() - return render(request, "rctool/rctool/export/rctool_export.html", {"form":export_form}) \ No newline at end of file + return render( + request, "rctool/rctool/export/rctool_export.html", {"form": export_form} + ) diff --git a/frontend/static/images/RMSE.JPG b/frontend/static/images/RMSE.JPG index 18dabdd0c2b27cacfe2af99b511a45653685446a..80a5c9ff84abe19647bfe9afcbb2bcd63d4bdca1 100644 GIT binary patch literal 25698 zcmeFZ2UHZ>(kR?B3^_*;P=}oJ5GCiFvw#dRL(YSO0uDI~l5@^c0RaIO$siy}B#C4Y z5lM=m=o>tqb3FIE=f3-W>#cwNZ@s(eo?TtLs&?(#wL@1o{e9;95Ozj-7J763qg!r~bCUJ@9zojie)B~k zNkKu3$4LspBt=C2Iqp9m`Tw-v{rwztKI0OrI-*dXlKlMcK74R{4_gGEorfF0AKa5) zkWYXgkdpQDgxk3wP*7WhgOj^7(}(t_Oi(9#X(m$cG|_mT9IX8Ivq5(A&3`I(?UNKh`)OlFt~R?)*70fq6w z_yl+{3&0H_CX5Vkd!(d+lFCmBm@8?fpQ_%xd6Vy^FrSCF1HYhzgap5U5WkQRFNT8G z$KM?V_v3ZFe$KqXG7I{0LWs8-~mW zV{yShsybhH{-8c@u!EDE=MN5yW&fbUz<&_^Q{~UA#F$L{v-=CL%0q#|sw{x8;Q)5w^VIVj{x4a6x;RowxuJDI{$7JHM{C6UL3;uD|MY zF2^1t2O%gZAP5r|=Cu<*U}_gc3h>&BBLsNG#f3!dg<)_B7!sq0y`7|rhqoIX(+^H= za0djxvnRrV33_fXNqIeWX(k~)fuEOpu5c6*Lo3as<>c<`_Y-U6l-{)RF zVKIqc2tB-wJUm>bnf};6{8pxoaejL^3jU{<-?kG*zi%ga1;ix#Za?~1&B+Jl;qCu(=olirenPHJ&>xN~3AZ~R2GUGE@EZtwrk_YVN4UEK0yCg7 zPW2m^(?98B5ka_+m>8yS?8Pxd4JL$v;7ClDi`dzV3W_5Uf^ae9@BBUGXr@QHhIJfZ+QwsA6i2ZH@(VsTp$LP<0K4yPw3jaU065|PW z=Y8Vg{ZE=K^OtSeZ*p9n{z2kjB>brAm+9hnTFe}D4*oHv$w2>;{$qjvSl~Yv_>Tqt zV}buz;Qu`f{CQJ^aK~(%ZenhEzOR7qt1BvA*E7&nQP)z&#IyhaT;OWw<^z!f05^A( zw}Gkx)Xdxhio1om)5im_0UiJhxAXCo*VEJbaXiS02r zKbUY96hk6s=jn~Yz&03I^d`#l96pbWQDOFf_89m)F2>`HQ4j;qoa1eOgLlqhhad2H zcnzTPHZoAewCy~u26gxaw*3XRbM$t@(7-S>EcR~h82Mnc-(dT5IN%(1bM?jO_9LCg z_J|Pf#)g=j1#>b0s(?D64d?+-z!vZYoB&q<1@K~WcMQe{Fu;^6{0%?d4}Luimo0|N z39!R(C;=XT8vy^o51jV_h6a=Vw5<BXA_u_aAS`mw z_a5L9#__Sv574e3~Kyw_GbI4`I5zKV-GvS%XkdE0MeqLZn z`H2Zo@#DA6=2PtWo#Z^frQ!MRBIFDrGe%GmPRB9Ix`rYhIXt>Q$kIBFDV7Zh<7cAN z)nz{*Wly|C0DGeK4dg);JRpL9azE(a6&pqttyc889BRe|_6<{B?=s$l=jYeICkf>n z*VrtXRR(~Mi3g5)*cECCFN%F=QH!k5JN5%B`4U&5%0os7Nk_ZEx8fycNl6BL7QI>~ zWsiZ&1i*5>Keb4uQr_xk?uiZN{^=WP9m(QhKUtQMh-R0TAFlg~4MPp;1BnTNlWcF6`&w z>wjR1oW7%+XrSC@+*@H7p6E%U&&s<5;8v7s%6*A(%}ZQ8;VG)(VIsOQ8a%_QOkxo! zP#fD*M3U)9mNbhIFUVOi-0GGT#nGs9vF7FGCj^$yB}&HRqGvpzU$|{4K)3Ss0O~;~ zODDBp8tIRL@30V%*ABjlF^I_WVhO>>+Spf-0+0i99>_E~WsdgbTg zr0_J@;Cayljk*N@LB<1;3S99)oB`nC*I$@_1^Fh+K)3?qO+@$IiUlbRIt5NnW$Y=u z6^X@zZe`QPXE87s=73nGI{}0Bp4(`l=n|5>U)lc(dy|+mGefK|Nj zs!49@ol_ZQ%Kqkp zs?axe){@eeCjai$+tBzsw@P68{NN9hnPlewW))VHkKQ#{u ze=584^7dyc)ra1^EqX$liQULKrsNg>>$9}7!yh8vw{E?pb{cp-vf6^3l+iFp#Qyb5 z=~8rYNlBEFXxd>ocZ$bSv2NuZ;;jAGwRkCtd=(1R^KLG79`|meL+{VX-2zG^y_kJ* zqOtpSS=@RK>^EL|JE3B}%}=On&CKxH?GLB06v%nyZ*%Czl&n2WPT&)d2&sr^QO$nM zpEQx#cp-Fwy5%!R^?af#L!3mN>VN47SIMMy3L~@wRNHJog(=IJt<^_#H6goIEW1(0 zY-Kw1lf#GafYXA{WbKBTYp#2aXEW`f`D*rTOQZD2n{96GT4%>X_a4HybR9ew^QVu- z)bYdYQ-W}%-q+(LTyGr-YTWaUb$Uc4tX*PS>Fse5zp>ymp6<7hvXPeia#^(nLbZ!l z98@iQE0y>2pVVETM=gEp{qU`-Bs@A-BgOWfu=z^reBGOs`h)jROnP_{D;leQcdHc0 zVme;s106HGFunQLvyGQ5rrl;%Jy~^Qlom-2G(tgl8P=@So*ad(LZ@j5c<#Tu5&9iy z&RiURU4O~yqE^Js;-qhoygPWV<_T^6mGJZrtJG4TU&@F*Z0NoJQe#rxP)MsqG^;9L z)k&qAV~pUH>*Kgqy^t%lg80-EJU0?ob9E1)935eXg73?li+SVh^Fx;PIj&{;NRf%4 z>XT@2YO9YNmC}^-99;NsR@wJY#tr$BGUP2c zH<>CJj53v;={w&`<<%tH8i-9_`()I^ZM#HPbta>tf7?+*k^|{BqWxj#Vg7U7OvqQ$ z>@4{8^{n{hp5R7W_b=^zltDa;8$_F8$+ZigbzT-HsB*(4gsZ7tM9djfws?)xlp|s~IK6&5TynQS$QFy{zm&<`$Y|P8%ZMnOS=~w1-hc0*3tZ(#u zja~DCiLnci6?oPe=;$f9fTr%`(Z%ziM+fhbPC4@+o+wt@! zM}F$Sn8woZx$*%rrf-zJl(y#ZL9ijzMP}SQ}tTZRS^@V zD?7uvTd@9;C<2FuIJ7R~eWP}2e7UDqZ1M3eLoL^4v3o}@ky~tegv_uLs`Sazcr^LQ zHOSfQ`g$;4x+lq-H^W2DKaxiCg4&Cgnkx)j37={0Rur=)L(PT4hDn7ft)+OBc%1v) z<_x_KnNm1PwPfP1>%$2V)z(#Q973^~EC{dv!-Jc|ZjjzqnyB7ZK2K~&d`I-ew8oBH z1ov|yk>c?}CbLKj5ocy%i_xMBJ#h*{Bl5D1LS%Ju^>Tzh-j1mvKa7@0%=Xa`+{U>kGyk>hl>*3{G0i+72lJQ}mtDW5_t zru(Me20|S_{JzeZ#82xNuO~J>X|{8@-sSE$z_Vb*G)wZh=1|rZ<{XC4KPsMr415h$ zb(-rg?BYx%s{IZePO{YHFG({$`?AoNtKlvFK_Q(AuSz~1W}C<5>A|pEez2>*QpZkx zpVEkTzEg_0U0F=@Q><;f@&l`)nc8b!udCNn8jj7n1wJbVnDS0YFbD9KsNNau9x(pD z_ArvgrW^BXx3r*#XJ#9+qXsy3w9{$DVQ4Q3t`r5@nzHb^S8IL`raT`&cy2t`bq>_I zZ9Qww5P9F_<4can0Uo1a6lvPiiiF>&=Efv!;<5H6(*!-3j>u(1<{Nd`s z^K9upciL;#m-#8SukZm~lp*w|lFMq*uS&X%KQ5;|RSzq5G)$aI?ck+3wdB3XjbXhL z%S{s>=;$7!nafL;n07fC-#Qszt-qAMepln^0DD}@-hiYDmMqQv`e#E`1CXm^kESR( z?nQj1C$|;ix)*dy<6974`Q?b#jP;g|st;WZsn#He&UrusM zxqO**vaL^y=j*K{eLq(j)>1g$^nI25?UQMs(n!&u`1fA8}~7Fd{oi6FxxbC_o>Z?53)IuAhE*7gMnB1BM#pWpzmbg2nzTP zyiBUT?Ws{a_xMUHrp>0mANs^@tTr$s5~vM`awJvpPBdaJ~UH~dsa@nq;o z|0Zp4QU_>q<8mRI`#nuB$(>-Hrfmuc>ODQ$U+fcipt*bgnqOZSkL*&IcZM?~g81+;?A&-)WZPOUWEIUES9B$Otl8b$_6oOC@0=*mL7liKZ3>HkoL)## z@iL8~p?5>q%HBS!ohQmCwdCi{6-O4cbA6`OC_2?0<0}>+&k1|WB*ya{K*swmG@fON zta`urZ*}O*5{QJ_WaEC(b{*Tm2*$)Q?Dt*_V$h<74TiaQT zZ#n#8>XgH=5Mnr?7xeP8#=6vWmi8TlEa}6b3lg`Qyh=_}_x5We%5stoD%ULz(UZDy zQoSKMhMf0w^m}M$fp(Ba~Xj8@QyUWuO9eilQ~;&MD4gsu17r&?1eUVO--9sMXIp(exd;MH?R1f#Py z`153!__t1>%96H~#X4Xm{8864)QlxFw{v>AUOwb@{i!%ABS7ff``CfcB-`Co_ zt_N}!z!UMRR9aDwFydldA!mQ4xAM!9cxTjomX&TjDI8c@L~P#1$Fe*%3sZS18P5JW<5y`y+IiqwY3gWY!=d35-#p93D!)7g*ViqX^rm!VGzMXiT^m~J03sF> z(2SpYDo~sgFU=O32R7iqp(1y-8u4e&TXCVeb~liz)Pk#?c>I*R;RD+josvRzN%>go zPJh$5>Tt&m#vHMHwUvL6NyNxp#rVo+LEcm@nucM1eS_kH6*MQn zBrq@e$cdlFdr3NiZIg<#!K%ow>rSR@mUe_|^Fi;LaDbLIy|dlM9&F})kBvy=kc z1&-lvM;o-}TkylBLfH)CPO2&!{1QUi*(ij3Z=nn6oAw(IPdlu6|I5s}9DRbSC3%RQH*_vkv*)@siWQV}-OvUwmqoEwb$~shXY9u{-`|curv_wIxs= zQikE9LVNLn4(}~4h^LK;pqQ}KMLfme<^%OISmk7&tjzF2ULU7na8Y`qdM$ouXfKhj z&RTapwVo&k2eI@b=Nh|k`nyqdMViAp7Ta@{ielEFlQEyTnWMh5@4(5OSCm9MkJHnv zPSwc$M)X~((;J9I-VXIHMsA;G5 z+B27qd$5dJ0(NTfn(LmW<%%wmWD~pl!FH*GrrI|J*0uHiB=`j-jp2Z7zMkrBJw-zT z))mT7{W`@5iy%^2c8=q$ z4J)qehN=>wqN~9a%NG@vM{Zm`q3I5qNOySQuM^S6ou9T6>;%L;?{ZV-1_5oDFdg~sBNm$il^KrB!PCTEIR$Zs-{qZoGRk6_*L|3dC9%s-Km zK?3j2!p>mLaT4Y__2q;a7%HN1IxfDY@a8w0r_&R)oJAtIxGO8>h*bRC4OWwrcMaC6 z!Bs2u5Dj&UPu$Hn21|R(;C!2~zQnbK)kkdfqZPM;44<7u)-*et3biDK%I8H-(&Uys zFh~hHczNCHWLO(U(9E4`b~$=-C;Bb-k`}_dkYz=gb2rj8Q!Ht1a^+nw&Lfx;Lm9c7 zj+j@8yM2C5AfH-gBxcWm%X%?DY;$CUW9`IqobbL|G}1f=C)(RGJ8YuOPYDM$+>tUF z<@UIX{6IOg;lT7wR_4vE$-D5V4h6YdVa2;`6XuzSn+zJlmj!X2JwTY>aAFCUh825~ zae0!tY8semalovHhEfeaQ2Q-6D>~=o-`-ju3q-Oa-(RG^YS&Ty2q&ec$Y009Rp@;t z6T^_^PSV$*x~7>@rwCR{y?3@RC`BJ{i}pNfOC?o%HG;LsG;w46fhHohFpb`|P$x+z zm5xl@R%5sc!P4L(CfcyAZvak}kM8Oi7Zl@W=F1eNAv1*)8LgvAImd_Ku}Gf-6QcoN zxm6uTO3G&Cd`^C*myTCTDOC6(&qNiodX z8jL=iPu7S`*I4_Yw|T{ah1G|~T*LUeR)9&*gqNa(c8j>`_BCVw-Jv$8FE~68z5}j# zsgx_6dR1zLI6`rxr#>6CSL9B~?{B1k2lUculh2$3`K+(xia%{tmLr#pfw#GHDucSA*Mk{_2plY)G)7ip_AsS1dC{gcxLc$86-{`p z$odu6>u04oXVeUPnT*Hj|G_#B}<9(>`@g{c_GZ>3QD>b4Y$cgVrnwq z^)MY5Ys0!&-xCH+X&nuTj%YN+Mt7ce)|fZj>P+vNm4_vxZj)+V<%nS)a5b5HOTiJi zaV3-aV&sTyo4}`6n^YCXEPA1oqVd^0W*H7cym)MHb#h>yhYFLfqPOy35s;S+GV8b3iGWA*z+sFT_!oIXT9s`6L8-fqnjyf1LHdKF3MtI2&iMtRJC zFIH5ZHbvA_M}+Q9+T|4GF;SUbxudB|T*aZsVqs2%m_WNvI-+_7) zb+i3ScD0qrW($>}bVe>Orr-da96_OojZX~Sx;&!Wb0k)SF$63Jw^y`3;MJDb;+k}1 z6H6(Om3mr{J62x3HnL`byEylSNc#oViR$atS^iH0$Z(UUXyqp%a8p+4(wZ zj=Frh0H~*I5gfF7VD5BTLhg( zYw0=hyjrfCQ$1F|t1;>xH9u?=n!fG8G7Im8awm{nv6pbtv6vjy)T~=EH{9m*GMY5Z zFjDH-c~-1dwqKUVVc8Tn)5w~{pOZhO%ja#kj`NizC)Wvr9nP!p-Y+KJRs$Kh^0m9O zY9D|bYAxuqN^m;qXtXX5(pdPFLWo@C%h>Y1CFu?prb0qfEIwS4HSxPnp%(*EH2HCg+ZXOV;E8;_f}Q25im~AyT#a__r7DnJip9 zD%oP>I~#(_~7;t#&dpeBCm)pnqJuag%)oS~Qh`bULbbey7&% z51SD8`J7bQJe8iraoNepSR5+Ttr28{ykUb`@W)rOHLmbL+8QpaurfYSx;F|D5q6&& zojc(dtVwjb1mi7krOfn!EN!BP8jF9=agbWc`t)Ma+wSe!9;{bL#vgO8rzOWF^RNPgM(Yw+ z)moa1DYNv@XgSThHligt-Jl{?B;txw3!REq;VD|#B}Z|0&M*Y=3XPqX?3>nkggWf5 zH#2Dzn%IS}D;rS`DiBHT^sF|A9UeW`68G5yV0b6UpG(a3vaB-HLb7c}KmFPPFqrz_3)PO$MordLhW zV1fv&gkB^c-%2^BK;FTQtx2MHds-A@)R_29hU+=SI1|aaH->Zdt5Rkv>hXH?4mpVP#KWa zJgZH8H}GNLSf&$=qT0pL)XNm8&a|wKS%rGh72%RD#HQQeV7s(LDWwr}@~8l7wi8D> zek;9NCkYcRK}iF?(nmD~RA&|LL&;UO zgo=Y<<>ERX*6taJ55^TLkscj|r7u?k<(acsoH@#?4jzqOu2R~jt+e24HGrX=gQUhv zpKx=Xx{DhK>VnV#aYxy$E8wv@M@4PXv~0B^yd4<10r|8f9)5k=7Pz zhw%s%L$sUE2AEn{3@P>?>4dmNP}RfmovpyAlAeu}*25nG7gj~4lzb=h|8j#Sx*hPc zgu`^xqU4T)P7>3Q%I$_)f^VHC>2rYYeXjANCp%P+{WQeAeTeVH1RZmHU2!O7FYjJj z^><_pfWG$$dYof;+fGpEnnzNL;ckQjOTPc~(1AN{g;NABAkKs|(`Xu=gjoCH_fF-! z+}an-m#Bo0zSh3eQ;}`@sY2m9fFec7E{hc8P_C8UE6_72PP{c3{jsH)Z#Zp!wXiMU zCfeWLKdw1#asXS(kft;xCng))Lz5I1&nm@`I2FEvU>nrP9;4QxlH^=zO3S@ZVsEJ1 z{_t2u!`f2BX$?(D=!<%)%~W!$NJrk(s5!ZXbJ$2*x213XM%Z+}VMygg5w7A4m6PuP zXvFM|?2<$^{pMF(5jY9-B{wy0?4||Q1aHC$8>_yCl@n7AZ9qQ~c_A55vS^S1(Vlyu z8Nt)^$X~5|;u8LGe)iMmJ9;&Y-5uxSOVlea=!|j8E4lTb`!D z)@i>kJCVhtZY@W4cBLDacvjT`dW)=Qt$jymz=>BPh3z=Ui`J92$SBR=xw*uRqgtab zo>vYGcyQ8E7mK^#jMW#^v{6kAxq%;N!(Hq}{8C!(65wxq%*eErf2pG-GwbR`m6#P{ zG^M~B-Iwy580#J(=20cVWXPrLOfTL-U}FC;0a%Sz&ph^a#YexZE|hdBqG?9H6k4K% zrKRCZl^F24=(xjQ3K@~jUsOz1_JT2(`At*~a{8CD=(BDlwW0bB&A$WTs$F7PFTQfg zkatYqS=Qe}r^K)j05U#ZO&!~~d?LpR2Er_c8)8fAY4wBI*V9o*76+-&Y5&4ZPNUc7 z(U7E|9(JVq0~qDT`aEH|mWN3GXQKdUvjG8Kp*pR>Pti$P>CT9njxV_gQo_1YwjJGBZ|^cK$V(Hx`PQSH70%M<;wU>3aMiFjiDQ^;G+QVqhLL z`VkVDn5>tcO$Z$RRdq6plr(4K1`jGBbg(sp>{aMuI{Tb18^KA)^LzTn#W2pM(K+3Y ztpqM^nKP<`R+$fC+bZeYE~8D+Dc5-lKeBR^(Sj%^C>!~PJ2(!V{kw{_dhA9fNk^W) z&{Dt3qwoy=Jbk~ER4brf+G%9Cbbi1o(%K+dE?@^Wb>-DJ(#E{AlEqh%XRrF!R@)IL zwPJ6pLm2k*rmGX$s)w=fBrM3Te&i~-(41~-+3eRG_Nt~L&2Ukd?Ja!2lVs}1VIyKa z!au!pgZ7IP4XT!-Y(~xWxH6&s<98s$Z5sU|{D|svPT+dPtSdFA#xUfY8Ykg;er9x9mkq)dFq{bEpDj^(njjC{>rK| zjQvSjBH}{bUuIL^_2)^Erm&q|)487Z@xr#U|F=#*rWm8hOQw~wq!cykwu-!DNb0xT z3jVk)NsyYjnRFWp_ccAQbtbXrpNY5oPVRsno4w&*604?P|BCa`=J}{@<+Kh8QVJZr z2J?zZvl~d_9-hg#aiiFr;SWW-Z>fu$X5ByJZ6lgCtDQ0>M`LB?3aLhM~?umsusxGT9ERdUbfilKkJ|3@XKil=sRzs%;ch z_PuBR=`i#o4&SJ+;4Qky?p1tN`RUcPNg2tFhOmPWy@3>s6!lJ|4#QrFhDIY-h_k^E z7$o(kykulhbA)PxAC@{kQ(Bx5Lc4YQEB4+TwVE= zOs9GDx6u~{Z5ol$Z;PJSm!`ZixLCUFG^(UCYTj($O_{>&4PQ^b)*rGTA1^DB81(+j zBi3pX_-{)^rMme}O?O0;in8rvCOV|b?X{y+e+MvE`_r4~o|-HMueBDLzSpb9Tvh5l zhiAo--XjPhm)P)EORiZ<3>fFMc<%XZ`3^Tbc- zWxoS;l8N0WF5b5Os)K85F(mHvt$x|xfj5|Ya`rjwwu6aDHEM%tM8#5VEY$|Jg?hr89v^Hhh_8NRiJM_dU6v^3 z7I7FFZr8Z(pHeZI67`$S6|8NCli^b*b)sDcs?v!b4Y3IB#?a@4vzu zkd(F|H0aMk+APK<=!=`|{Ou_$Fg6zU`DZWAA8e5`0)q1Ta5Bu}tU{>G1I**AXM5wz zoisXZA0qe%KKJ)3`)N5cUggwLZsRy6#3AS6oENus zY?q0%p?&Vk&6f2Y7`nM0w7&F_q~XrZSN#v*c+^26idJQj@>A7v@(L@A?GBOFgVNxMrU5Y)sH_PO1kCSK0R*@4uE3_?#zMM~IirJZaRG z$I=KwAhx$gU>ogthSD$_K!BFCtZn$}y{qnIS`|gu=Js`H@GAt)*z{oE%l0eW-vO(5 zh)CnLY|{Z+^aW|v5E>7 zeVh5+JA%j7GDcBipTgUU4Ddpfx*z3gCr@KxZ=$k&7x5~(*S3ySqPrx)yRAIsfJE=L zn<(SD3h)67vG7^kN&&l3y?OthU%CZ@UOwHMlwl{PyO(7io_AEWbM+shro$U+eSB7y4DvbKUC9br}t#96+_dPq4Dbezs?GB73F% zh00sV&Hr{mOjv_)271krMl^huV?d zbl2Gxv|S{!axte32WvjD_{hjCzGS*Zj%BVtKA=28W$~%Nm9#P)9yVZp((D7`WAT|w zps{zTN9i8D7IXDinT>iQ-tp2;LHR`yslDJ9ESJi<(dR(@s?p+IR_Pdr;c3NZs*1%P z)O(yc*4&w!`z8GBVNMMj(n0ywxpr?T2%4x~`xHyA?0RLXv89fE&F@!D~E~i})8n}qvxjg=W&10#pwB6C|m1Gd!w+^8pr%wV{%_IYRZI=_Ci1h|`JU^mC zF-g6{^c5a7(QDZ^t!_8qavRl!8}5)3l9ZYCtI1!+8{^|ikCGk<7u;h_Z6(ol^_c?aCS)3a*=`N#bCeR zcPYA*Olu2sU+~D3Oo1HVuce%nXVN9xxo{zr-%W2vkcC%HO=`b--AA8d;CbVaH6$Pw z-u7&2#d%7+#9F7`w7Evjts4qP8Z5rnQWrKI0ERjox79iqpv&2p%=JEG4*B1HqvgDB z7E2-%cEsIR??o01W9a8gnNnm7&D3S`||#_DP}*3eNon1)f~mOj@1-w z9!S&=QYGe_drA}8(1^!Wg3nxn-GM^@P^V?Dnz23}Q*3=1$t#kY{SI3W9k4f_QYQ79 ztcXNA$u#b*FmawmWg6d#PCGJUoIaeoEbyZy)9&QSb^O5u^8r%ta=Xu#n%h&_!UnC6 zW$*9d5n(Icq|DX3RoZu&H=nCh@h-?g3@fzN=16Mwt*T>m$713YZS8fp^nf1IRzs}i zjBp;E6#isw!cpm*qG*aToi_IrTzigrKFbF&-70ntY_Ts{Z;374!ZE>{Ww`x!$@X4IYxZHsvUw#V z!y(Tg&w%p)+=QnitS2sBjchGscHmww=uh%bsTbcizoaw-JxxS+8mW3oOiIerW%4Cn zWu!^Zcm&v&`|S!gSjL5qgbZBNy1mlJ{H{w7?9njAP2Iru&QtMfugz(BJO)Uwe^`c8t&r(_2h!{4ZsD!z zHzSq9+rV$Bh3y_xDJDpa&jxHU7M+x0CCC%h_`8&9Q<$M2j~(MNXtCK8Su5kZA*}jt zqN0y)?^nVUL16-~gZVT|Bhv*zAa`q#_@D^lP=o1m@jhb1Th)36kEFfcmZxy?=4NnL zjTBRe$P)=pz2#f!)XU~@u**f#tyC;{b#S;mT|!#Z?B3qb?v9% zs?PQ03|1=jSrCQ>w~n_=nGqEfT}@Vt9!r2iuac&Pa}2zw(yDt`$M8u{m2J*jWFhk+ zHk69q>6D}4xicDMhAZt-_);=Vq57JpsOeReII#txm+MR)Ck1Y(xB{!Lj$fCvxo8zF zY2RqqTzs91?6e^rYNM)nh|j`AUyf>?!fc`uX6r+b%RF?F$iaf%M^m?X7Gnm6Kg!U{ zOh<>!F)Uwu$9rS(eHByL+C}H-%hARbyZ*Xv)_!6?wzAlqb`OkIuMb!N~a+yj%O@hbUB} zd2ph~wYwUq{E`&;teQ9lLe;D?Z)x&VcxC0Ph=_==uLwz2R#sRtU8J_+=d{}+x5$$FL>VB8*5N zwBZU)(z~b%XzWWYzMeU8=c?d`iJnVscW+41bX`nMCKZus)6arIOn__13pg1{$5q%= zR8@u^^qK+of{U5&j;xWE2yX0D*S-dJSa>Y)RtbSC3lz#9VYXXU{(Dm>wS;RjL0aY#=dNEG)2Q_ zh3a5}>Ly2*RMBtb5u~7?f)K^f=I|*wPa}x9K-pYfyknKF;VU*Kxi(1fZ02FMN2yL>Awyi}qnie+l}gqhpqGs@6qY^CU?6+fmIdqPCBAy#bfe8`xA6vJl$ zZOMEF7nh;{Fs?Q({`@!pae=SFi{+sdpPVVsJ`(#%(T@W!qHVCz*8SV92O_tUz6M@` zD zX|m4DoSOb3efk6__Y874IUPa2GYbm1XkrC1M$K*VUF$Ean9&T$+oT8NNPts31BXuI zPcU;_8|on4H`O*=?Kr8261h z$95z1`^LO4x!$*Tx(3Py9uLIUKaIMS>w%ixPd)oAH^qR!DpU5ld#Jcf$*Xe1{heI)p^eYswj|c^f;k3hgl_Y zGxJ^zsKq6GQaSc@rh0;g3~JTvByX|d+Z$9YzDzXv&fVtC#2(}!M_a=uwi})7C4nz* zi?lkN@^@F8e+%QP^QiBCtpg}B-+TF0`@|~>-dlg?n-}IT`p5qS1LyCMk|n9L{I`$u zd-T6a`2T%r>MRuh10DWPu!N{V|IM}oe^;ISU)u{}Nx}pEl!^Cy=|5D2F`hgW{HM%6 zSpGpD2*xsGKdMNSV?pCX$?&1S*$Q*{XDt|-b23N>07(KtlGHH=kc1CSQe&a`LxCT4 zCE}k~lqe4+Q#;o)M4bhYX9>ZBQplf+MXN(Y)G)j(Noo+ZI!lN=3m607K|??+iI^gJ zNRs>y_0%8~=NJH^e+Wj-Ipg^k3dr!!;lCWfc@Y|8lYi00YAhl6EMR#6EQiqzlJsAN z{2@~g2mwJ7e4`1-@Z}*8wC#;=F13}m*{aZfk%u3V-g4;XB3nt(}zQC&O?uTn7E_B zgYo74KO>K1-+{{Y?L@TM1KZwM8lL}YMsCMKs|5^vp$@x$1Yc|x1KsJN|y^5uO~lNJOF@Lt1LL?V<9{PXY4~W3Sc>Kq5qyJ7=8K5 zbtoDTf}){U!Jus-GLvTlSb6ejn*<>6E)_xQu7?r4d5-F;4fYR%dVhv{0>BvU^9xY9 z#3+}6uZmN9(t%}qcx3cLcb7;KA@9Zs!4NDBH`Rf&d=ss$=uQByH4&!>m=mhE*bgoY7l+j>F@U#N z0LppERthJ4W78E9Oz3;pzLf?thNg697za*)mgqS>0B7~13ZTmm0L5Xb%z8I79pVcf&Gy`7_LA0%L$ z*9c4s8b<0f+sQa>o9ROyP<3X2TvU}*6egjgWcO2N3MvFe*;uzY-At86R6oQ};OxCe zyYS$3Z)Mgm@^F@O-g$JRDiy z1*s5w2;He3Bc2q%$QW}i4f9L`SU(}(4aOVClV0@!?LdOIX!kI#jmO&wLCc}#3hrS) z!+gsCB&ZZIj0f7nu}-=4l>@quf1Qi-+Y%XF;>k^+ID8uhk!n8{Ksoq^%~RVq6L?VB zj4$ZVF5Tc3IE-{w?Sf3E1%a*9We`kzi@=p=IA~^`c>xFFi>){Z0VopykdV&fs0v&( z1R`+2E{GQDtbBE+81wZC+(?7Ug5o7M`GVNQ3epLH8ad6|fJdgEl|?$Bf5Ef8%W)#) z$=2TLYGl%uBbF__7;l1VA?3AP*GCa4*k46l;0^dup;U?_vdIcn$O>+T?YKwSdM{Eu99(@y0Eee+55W2k z$X^Ktb&8cer~>@#{Q}TvNZn-pc(9y@MZ{|e8)(z@5$4lOq;bUq*xSs@J_QgPS1(v- z4;V`)xPade>qu_@{q0A`NlON+6cjDTyC!0f?h##czZV5~Vlf03 z6aZr|oMzUsF}VdPd07zChQ z&nzMM4nXOuuJt;G05o4?%YZXEpkeqM@R3#N2`$bC8Sx{uPlit6I~*zMvmpuq%LrV2 z+td+M)ZY4<3ef{Fp*x`v``6 zq7_Zqx%S?on^f65Bc?9tlhPJzfj|vxjCGPjNJ9%81A!|G>X=L4A8H)H$#5#)V8xGc zE<%nyz})qV`mZ|5#cOUXfx!a9S8?`w4L{w*&Z3%&0STY7E%uq*o3|K5|F25U{~788 z4&(37XFngizOCIIme$5%T^7qK$0FTov2rvh*(;7Eac1r>xSyWqFL<7L=9ziLYLUYg-IJjm<_1N@O*lA2 zD(mjX37RQ3Q@iSaq139ntTYuuPSVT7nX@5Pt8;?nWgILSSjYwE*eD)`2aIgkV^d{V(wTxjcB$bcqi+NW$S0r{c{I^rM{o{OzBOFZ zj)J9v@}7OEAEySnyFFnp1>!ij9o)9k0a{tinUKW^Rv+hUL>3Y$Eb3PcVgaIcpOjNW z&h~1K&o|pWBzams!|r;YrxXlYo8w4`XI$IoWQ2Neq{LU3h^jF-w>#a&Dc{pSQMUrC zk?B@V#UrQ+#p#rr4>u6TbMjW=S~JgRptF|FZx`J-aDRAb!3m|Xk%lyws%BwFzimTS zuW5Y)e=0<1tViaaUv#k!-8q=uha9D~_;;G;?nk4AR5w~>HF zw5cfDsVpT`K4F<&@6s8$~#F}D(z$OLzZE$ zcb#4XVVRCuu4wMk(d;{~6KfDwDs8$mf_bCrE@N{p+wzpZgg)U}q@2)3Gwt+-aTS4Q>#)vTfE7OhIJU|A*_i<}sv=^{SC6hyFnCS@%s8P=g4J#+*_{(TShh=g6Z{zV${l9A>Rpbqy=lXePs^W#-Q vecqb2iy8~fGSuc$vs?dImEnlnti0JBozXPRh;^K!ktF5GB5u#vE|2{SeRDCI literal 21493 zcmeIa2Ut_VwlBWvASI$8ozN5o>AeP2nurJrp(6rPLocC*A|OZ+5Ky`zDpDf7H>rw9 zCsL#&C>;`Ngb@Dm)cejo-*@jj=e+m*-+O11A(NfiGqcvLnfa|XYvuIY=>l-}rna6o zKt@IeJSF`Br;Fr!T5y+#0AOGMhyVaU3s91&0OX_`69ACijB^1fNqI8T%PcqV_dCE1 z0OZ4K+yD*f9s}vKIw=9%Yovd_zy0?G+zWDd^x?bV$amAz!wCSWo>KgQpDZN4y?ILc zmwZ!N>L2nHGKv6@miE{DUwKQG{hOWt=IvjFN?KA{R#{S7SxSmeT1Ht)R#{FC_T@ow<#d-lO@ z_DT-I5Ou~gs&HjE%nRn|XU7MJxqJ93!_@@-P`ffI|NFLtfV!%W!$W1`8`^&=LHeX7 z@Ta~61_p`;%7}aVI7vt;DJe-vN=ryfi;+r*`38CT*}=s;d4#De z1pd^#gZ*Dx_wx5~|3ljj_7aZnjxa|LKVK3zq$C6+{@Ce%RhuLh{iV^r%YG7#q%z9a zeH`ul9B+`4`tKnYlawYAbMMc?E-s^>q$=?@( z)8~)P-*xo<=SJJzMg4zRRQ^aN>Fyshr20w1M@L7eD)Aqi{?EMl!yS_Bk)HdDq>`TG z``>W=|GAOB&UQx+(i9IQ&E?Y>;5tA@qnP*G8mws5jA()R!*Gu1hMDNXA0Ms_p;-Yn9QsV`{-uT^%k8joOwWFGiL z(Ve-##?HZc@sjZ6D_|VDO#nso(KOhhu z6de61CN}PId_r3Kvy9BF=P$Bf6}>Kg^Y&fI`>N`i+PeCN#-@*5-95d1{hvOMj*U-D zqP|T{FQPF^KYlK+tghj9cK7xV4)MQ^e$zz;Q2ddszfkr!bTO0YBB!LJprrXt7a4gV zDN!&}Qt?YspVKs=vGYDJARS4|axL{`Whb4Wj4_t=fzQa93qrDs7jeHy`vYbFJ;I{? zA{;<`=&Mt^!0u>h1&L*|yCNLYg6 ztVncnpz>iM(@_FW0`J;Ri5M2`2IE0ArvPsekrUYo-aCV%AyRan0)6i&W}K1dPJ!4` ztj-B}i{>d1fkoifAdBw_%%^|`@+@LO6-)?@P^oxymoBwQ||5rTEx|A*ctKCh4 zSvMf29!oH?Z=A)}AqrM)>TR?ZKGcJ2=DS5Rv`@a5zqU;A{H~LVe>SwJ!P~6Ga)<2y zwq&D)?KpXhiR+<{=CkMNd2MBXTm>3?D2j6HCo(?tQnZ8+KIcJEV!)o+P- zMv{n!@HqXZNk8*4cYn*gOXfGD{aso5OH&)JsL^q7H6aNuGoGh)T888tmXe3?UdhjpiYW8F<6 zL0xv8SGd`#M#|97Z0E*U8)WZ7T7(o)0d9fmun`*GO7EV2Q&RE$+98{Q??9{$10aQY zc6ax_1BP|~XSN4ZhO?*dLHe1XDR5tn$BD14D8Z;P`rgU=vQ5Hu)qr5fyq;Wb%QR+6 zzN<-TaG+Xo=ynkD4M7F14`OlIv2{`RoxVEA;`@nNE0-Y@mwxzkId;X%nZn{9hF1R! zEc^4bqKWt?xa&9YE>14uwX74}sGgph3ze}F)gBvZ8i1$mP=qr@^{W=QZZNm1&Y~q3 z%&;56DB)g{^M;;hU`^M`!GSuFe;TeT#k$f)hkKKbnP=!JUoqs5%?HmQJMB z_R?I)r!8cP`d3ie@??SEe5MZ?`Qi)lV#BYm;dY(d#QosSAmNP_9s;Pssvqh<^jN%@ zU*}etJ-8SM{?7VyP8t7Wu|2L$h}bi^cVc6ia`7R^jln&}QnAy8)#1p!Zazq7>J_X1YP65m~qRRv?NQ|`MI|Vx65d1z6JXy4uM{m}WunirNT|a)V zEs(u7auMF^+a)RFyJ~UseF7_BViW;%khd8AEfW4aU7Z%~ZvPu3+)A@}*ljJ0zMbM* z;v8^Thn+E8Ql1kKy)exs=qB0A9$3g!l^`08A3>Ld@?+ncwADQA3sQcf??gjob3Ll# zC}vitfrV_qOMr1##j7l~E+2Tjc$?SM)>3q}CACW8dMA*qM;Qft+DoCuYj<@v z>3>@tv$5hRmKpL6@*2~VeFRvg4?it>L*)YjROQ&>!kJytHprLnQ6D+=iM4NqlaUjh zCR@cLBnh?{MsNp_H_HsQrJCcXV#MP3(|6=-pfEvi0UB))>V{9p$#kc=*fQ>wS}U`1 zC5{$myu5`8O`hjFq7wcv=mi0XdMi8$7y3S>bw$}DdDF!sf>7MqowB7*F0=BgEgdhm z01sux&h-jE7B32a5*WM`I`hpzc;sGWidFs|*^lt6xW3+$N2m*d>Z0Sa?R>iRK3-0Abqo8!qcsL)$D1c7fJBxm3)pYyzEXVWD1 z`mJSG?^$vmokF=C@W+BPpKxGwd3nYAZEv!`)^VYL0o??jDafXjiSWsL-~bGAEdZ+` zUMG3p9G7i9^VwwHBi-`Vuk~+*&P(pW=QNbIO;j%)UxUZpz~$iG&@pA;=7wWgls8|} z=Dgx;M`&4|5<4c}Z3U0fEiyx!+|Xb((I$@^CbK;o3eQvuSHLH|1TJS@RD^W%Cn~Hl z0%0GSv~>g$IfVs??iZ3^VKe=8ru@~RdlehFh_cc%7%RRm& z=5NvFL_w@hT&QKpoA3B1x7n&)d^lk2g95bQGK3b*sDH4<_G4yIPs6`*mVZq_!rw0L z6Uqf$0zqGf=0>l(dX*+3(m7@!Q}3atdfhKo(4jkJaGn#~*DOIwy*pJ!(~=;db89oUA5y zbRJw=o>AO1D_|n9hLNUl5jwD2k%Gf(TK@U$Yi1^+ezh>2l8em=`+49VPWkqCuBH}h zJJv(RpZQzvnRegS#2?*!RW z!)zbUEAIpj9k#(5Njkna9{PHHH>{v5IJfg-P(WH6X2k2!ayDrNb$a{=LE$G)nJze6 zd2kBoK1JLFaxfBIwBdsAJ140|a?NK)&2Pfz*Uy8s`xd;ye0PSrJbBh2^r8L-FCOG< z91qqNKbRMlj*jkJS@Y2jS@OHowt>t&1zbQe>Y+aHWsPWQfs{?%vtd8D=^_pm6zV_^ zoGDhs6pU*oWt(9-26^IQ;zj0m86bw&4xF{NSBPH`tkf|j$msp4gVp+~hrLh_v7eC6 zjS9r2qjX!-KY%)#+)y_xRa{~lcunSs?pLH8h1Il`2pRunY5O5#c5R(eFGLbf>sIOexQXB>V^mhboVTb2%R>d&GH%cr$(>hSXw;Bw;VT9_*xqMhxo$Fa5)E18%ry?0rvhi2bEC$y**R6KEf+WtF1rm8L8_kc%>A@2Cy%0R*K zx--F`YI9J$& z>@sAs;U9YOrn6ZXBB{Yu@J_Bd#I?ST?LyVLbMS0=!4GyR>l$D(xXNwM^PyTJWge$M z*ujbR#oT(8j9fNh>|l5#>f;Ns+PML z7)>V;oNhMhxu{F-5zpEg5+eo7YS=2V;w5FDb8f#4I45(m^q9Hecy%_Mb@2mf%3@M_^h-8h3b?J-4_>0zaiU%D zP){f%TMMr}D2S^Oh0xU)gw;b9>iyyq;-&kSldPWUtYa-XGU zE_8G<5vGdO)?6Ud;9b=62Gh0igW}cctZghEYMWyN>Y8GuZ^i-K6+xX0pn5@JjtXCE z{6J6`<_*@V%Yc{1&E)Zi#^%On#kFl?v0s($gD97LJ3f)H@;zvx3-|)T{0+gj4hv59 z`_8FXkr2}~GN)0Qd2Jw-UudrrEH0BDa!_!Nc!{)Q%CfYM8&!C0(dp|eexrU7v6M~_ zQT(+0Zef_luvBhsl})U85JolKYjH)q$VI?d=~sri9{4(u6+KVcCXR1L>m+(!-GA=4 zQgdvnoi|v+P~KAYzCs^t3W>RnhCGI_WG%uD6KWahH7iHN%j?7%(>&|6(;vvDSBy%p zd3P=szB>gRPPqE9ZOvs>_rF^gSEbKt=OjZo5QeV6-h2Z1>`T1cBI{&v)p2N=RdI{y zTwq4)^K;>_@)iRE+3eSI_Y9%4MPqe_bryqHIj_pB?dWii%#?oD zJ2WyAw>oC{86h@eBV~_E@3w*94CYd>Z*Np)&kY6|ymw<-XOf8s5dI<)ovf#t%fLB| z*TX4v&nZr}Trb8)>yH)Y+)%elGp~LxTRe4*(9`gXtH9$K*~)bLV?S#$fv^zhhp`I} zZooo(p2i7QC@|+bm_{jo;%sr}VJ9zp%a>wgD;R3ljuY;NBu(<3Z>S$DYN*P971k## z{-`(gs@NE1F%@KIjA5VYD>-ud7}t~Pf?~wf4W~JJ&YlAH@i*Mqxih7AbC0|efE8yX zPYUKjBL;E7+8IGt4JYVRqyKrr}# zx4eRjO)s4K<`x5^4Dh_tSZlw#xN|f2y(H%Css<0DqadDF7@xvh?TcpLD5|TH(g97q zylg(>?fxT0gKZaSgA_W+uJPO$ipE-3!NP8Dxuxef&44-t}&IoAWIV zSnHU!m*sDUsx7TPn9Ym7tkRdbaWx#f6u4?*c?!_!Y%-n#*M#Fvf$by2j!6iNg$!+r zeHK1Ezb61k5Kn=(r@%#{zh7`56w;2%EGjgV=vKQlbaQ-iHS@CjthB7ky;#N6U$a*J z8;bjD&aWa=K?fhkg`RXCUC4j^sG%*$lIP$O4KcB-eM&~xRV5r zPPoYkaa^6NX`+S9vD939lr&Xlv-0jax!;zw#PTD^4PEFtGhiw-hbfmpTqNkKmr7KIFd00ZQaU6*nXoXuz);&5|5<;B!>$`v3JJkxcv52i$ zvo@CCTCdw=tN~7ehdq3JpO@?h`wlgPd0#Jo=%*QN3VRmp@xslX=BU)!m<7Wv5&_9D zd;RcC@$(4alW;A3n&31^aKh-PFcB_~Y#@!t&&!vV^rL41?`FEMqj2t}qx%md)G`|( z@lC|?MO$Evb1UXz!q9_PP|Fnyt<+!`m$lcFAmvg(pAxribeq)LZs_fT4t2G$t7ZF+ z4^pjGx90TgzAxCt^GWY8g;fsf92>eprcQw~7dj1xO6Ikgt@dwxc|>$>&jXwZ^!P~9 zT9tjNv;>}XbyPIuQGNaKRW0exKy{x!Rq62BHUC|w^ZyV<<{k0+u*#mBdtJT_#`K^2GiBnv?S zuNn|%@qOLkZvk00XNT<@6qG#=p3JRD{|cB~tqqj%DF(G!<_=AxQEfc z?)f|^cW}B{u=1qn=a>wsbipj)0TV^a^iPfy-zdIcMqb7$q=XB))`jym;cTO9sK>h9 z9Gzh9SUtauFSH=hO-Czovj8QIRkNn7UI+)?5?z4EpVSo17;LgnyPhk4{>u}HRB>d2!EZ9oNqg>4L=||cZ zd9m{GuQR3(0;IG^1_pP&S%m83cOZ0r(f0r}bKN|O8`f~%v!ib=#sTsUc?(~Sty@RO zROC6g3BoRo4$kR3vIzR<#7AYnTOps+;eXuQGPGEA3VhT&Hp7>i#_Ma^AEf@|4Py9^ z^2H&=&Q|vh&0ssWe|}#D@k96&sE3viF5zuS6R%yjy^@F`jvqU$%bx<(K`N3By$W^V zweuqyR4hh+O{}|dn_`FxyrtEszZj^eY2xYl=#3Dr5LCitLScdQBf$!yZ@wFN5&DP_ zf|KbA3~t2dp&3OK#SL1VJ%n?<@SDUBWjBB&NpoW@SagXouy_7^8xwrBCP}Be+I97; z`8^w@e3j=piV$`4p#Fs@+BL1PyNj@JFbm%B7D}whxv0Hf?-A_NH679I=+tar@>PrF zL3)B5`2ubVHE?=R*EyHI|Z>Qj^HfjdQp^Wl3T@!3>$<>)XU zv)iSV3i5~0UFi(z7)J-Op}U<2k|3lj0XV*t@T>hGpc-fU1o;*x-ikg64ucpJd`ceG|7D}{ zLf#>cp_}#2XYbk6)uA?cPPe=|3wH3J@qMMG<>kbb`cmFVokk(DEoTbjgkSdBxWn%G zq;Tn+TiBmvZQNt%>sJe$7BqP^IJ^|9!M7rURz#M&ZK*@`NU*CNuUoEuahC(!MixI^||>u(a#QkCMpUaGAnt zpmK?03faiGa=#?|>p1!9H?vl{u>8@{U+9<-P1y0hlpG=QkSMjw6!Vvn&(zQ1L;BiT z+~q`j@mZu1<+@@l1aZkn$UVN;{^LD{!m}bUpicsiTtqyEF!+Ut8HHbRYbq}%sk%<~ z2aSu)vn6(Y*;ACFfjpNf);^2oW2ss<>$?W2)hvF7lw$19&+7YLNz+b!CgYeQOo_o( z#Qt5Qfuw4P!c-dX!jc5n#g#pdW z$Nk}Xnp?DL>qUP&(yt7)TWm(K;URa2J-&Y1E|iMv4S?M9?2LfKz>G9_FCIhSv0;}O zuP?U8d%Kl6S2xzje^vb!KRe@#qRc#mQhG6lz`Vm%@b<;ofNfMeCq z<7R}~1V_`{o}*l?Z~Ro!*qeRr5W<~+NA|SwoC1%8(LBT?1cjBRT~N1=l&v~?WSuGfPb&|1hW==Y)O>y4J`P3JHIdF?={>vD)P&t_Dw5gIDH1t(3f!*A zJ|IH&1ED)~ybynYie{eU5P7W-l!n6)8UO!+IXOwIkuzN#zB<0%$66Isc{Dhu9zYF3Yv>kc2vXF~aJLd?haK1T$Tx zIDb%pgVp6Mn5jxpzW0<}a3fDb}E&(+AdBh%c#r%kvPhZU zOcxnZ<&3bzJ_`Pv_qk51)@Adc8p{H}V38hKGrPK5L?_gc`$fV6Wo#Z7&IAFkGW38+K__#f=_L zXjstf;MrZfttYhOs{xiGoH4}TYIioU?;cAUINR7}@I4#)tq2SC{M?#_fj*{s+sPAt zsN0~|GHz)oPs32CTFR9ey4g?=r&8SZ2Rv$u>SB;DLWdtdDz+=VRw zWP-UqHy6^N$#-v`8N|3bT|xrdGXcGzHjb}qwEGofJ+{?*%JbbEs?eA4H-{RH`G)dF zunSg7><^PQ*6s=_t&bm|U*dE&P=)KkaxoQI@2y-8Gs^{s5USUH>i33TaACTa`56}- zncf2gdq(Yf;&oc+U%4iAcu#}r_&H6taPsr>(7hF-B(D!QFPgK(l@l`$nZeNPHP+>WQwhaF?3C_4v&hO@zkLMKHaTjY}oD9AWC~*N3=dYS~VRX_(FGX zVB(CWZ=M7%v(AT}!Dqb-)xx$;p#(JFtyWz%9^#BK?;LQ=;)`Gy@#wW~R;9P8u*FQoC{x}}-+ zuvgl<1MmWf4%xu3&-yUKKHJ)WV#YTWvQK6#r!pshx|uZWw0$jxx_fwq|I|6zpKQHu ziia*Pk}&Ho8kw~1zy4zRtl9LO?uYAF)EZRK>Yr=b?UsU@r50oumnB6 zogkrO``idaToKd|6v}~Gdcyny7T$NRu>B*Mk!!%m9(lB-LlNu>;PDbg_NR($P2>k18cDzBrE7 z{OXx(Mle$e_L0kZYz{E^?lPZD_Bw6q`=6vPv@^g%46idLrYq6rUN$o_VOu@^N^`6} zu^E0I8y>%q$1K~F7o zfhw^D(Zmb0^Kh!7hB4zgH)KuRF6ZLNZGB#)&lqm9y?jCc@UuvJIQcoGsZ#V z{nvYCY&02+;<+dN!l;y(Vx@$FU4fC9dT4>gTZ{mQweWgZSM;(lV( z4w5J4UD~VV>axnL)d~`8-CzRiKe9i53o^VK^tJD)JbRdN4HDrrpr{(;hPI_yABB6D zr?SEptMsB^i8>iqz37ZS-JcT51`arv+AiXMj&<nN&ZUH0unDF5qF;q2pXEWQA(UP;0S z7fL;RUEC{!0s=ypmkiN@$%be;V3-=Q@oPaJc?K>HvozN>fLi%CWH`_S-;k-agmeeW zfDa>76ZYT#+41DB@1RtOE05>#b#zSI3ob|50TyjrJ?<Y1+ws%p@iPpkV;gA!yT1DT=f9iot=F5x;l{dD`CYZ6s`)AkPLgd7!wD_RIBz;1VC z2BYG(ksjsK-k!cZQ`#ctdZtt+@L#&UAfrLxe(0@9CwwC=*=SM4+|ww^e2Vh=3+CX) zJ*!M5H8zf(u&{+q3L>*7)v7Z1#n0NL_t$oq5<5b>L!9K6AXQt&6qkr)8VOjFV1=m161`Zhaouk*T3NV4!+bsi5IJ(g$h@wHVhEcK8p^X0?sqw zN`l+?4+`=vT4C_v<{jcYXNuW&XArlQGm)nrn;+R`;M*`hZc+La%UV%R*>vC7B~`fy|Qx0W7rhE~LEr zwp!Og>%NJpPb3>Lx|T9*QY9@J#6pH+>oknNa0)PmI$rE8Ij{Z#!ItVXQeh+RQ$Z8T zTQr=%0G*O*i@``Rl}6f^yRap&(`aRLXz(SL^eggt1v9USepHtuP*xhqANSbMKeqVf zk)-xBw|#92xR@_R5MsKyC-??L6MK@U3VY{l%M9}aDHZ5XaulVJ`;>@}nriW+d}cX{ z9@5kR+fG|pelzcJv-f1Jtr0nQ*-j&#Gdbe&0?5u*EPSGAxC!?rsi?@&-97tZ=gLl- z)PA=sg7fHo+XehUX~9{71CA1#F4E_Jpc}PYaG4f=lp%~O1o`0E&`27*=)$;<8?V6V zrpM20w}B^j17f|u_VAji#ET?oP^hAO7YdkN+Sv)Y6nkG>f(bqRdifi`#1yE!YlR@l*+xNDSkWK?4aw0d1^NnKk2ZWnVD$I&MCU;_ zq-S5kvo=OQHEMP$X(=vGHg4mcA+eA>8N^b|&?&%G$mIE7dN!>>77K~&r@$~t^wZcN z1do)zsRMMqx=DYXT^!9Im@k)lok_zkwIJNpHF*xLQfgKqHWPZ6&(EDXLy&Lb_}(+X zNQ>onBKsKX{^BKG3~j1GC;w7u&*=^#&jWr2Tpg-_9p|bKvM_k>(JS-#-or+TX-<=r zg@+AgOb%|2?EA+zoQ|MDh1U-djH^$d(fyddk)xyaku3Q##kJp#Z7pLM+m9=4{CKm$ z`w2}y^NTmZZ?l!3$M31Ng2qwB6Ka;GRl_wQ<#$c;U>lo)(gUJ?CScS|D?mEEe1evn z?*ehw1V2voYc8SW%&5Ndor8w=K2IS35Kphp8};wVXN2>Ma`;7D++CbMF}p6DJ9bG= zNYncu7Jwxlr2$TYs9DaVsz2Ghe#c5YvOS1B1zv8GLN5*%A-^vCDL6ulw2XKj@e&}X zK=}DnAl3PwEEy&Jnj*z}x<#MJjF{#k#)=Ik?rJV`$`B&vp>Nb6155$ zht3cvhgEsd@Y+o-YaIEUpd(@6vAL76Pe0}}Wl5iE^or)sR_*^zRUp`_e;50Yy!|&r z{zlmT&5-|>8S-u=)m74&Gt^=UH@oIN-CLi2Lo)3a3ytG~m!*)GdKB`Hj`vn6-AHa1 zA|=*(>t&E%s`V$zczZg}Y?nGhAydmaI<*-8iM&#Lpcan2yXhdsuB1%pdCcZ}gCENr zy$jyCO!?2qu&9sC<>Gy4WZc&NDyFWozQx?e)s-B|IF;-aIYX}bshS`S&Yf(sHK&6C#5g$fYyQWH%ShZsq-OkH%iJW>@^?48Iw z#*QI@p_Gjij3tOVR7b4bTCAerokeX;x~^Ah+jT*O!p|8)u!((QD_TsnOPQNQYY0zy zqG-A(xHUgg2Y3D(Q*9pUd=5j3r@MqloB~u{2%KDH$ttHuz0-k#lVR<+SUvTQ$AbA? z!I#q=iNiF5^VRi_{Z1^=+X-?8H_?i6239K5J%Yd}&9hGe$P%)93jw?nD8d6ge;Gj= zH<4UfY`c1pS>nP|+uz`HF0@MNYT*@{A9kIB#liUn3&=C8Y`A$B8_tmp5A;x$H_VR3 z=aME3!+Y9(p)f|+Fby->+EWv+h3y+i(A^ZMu90+^Tq2BpX4wS%h3bWl)7&uav4_h3 z?L8pS_jME<<1+=^pzLTJ>&?85%Dq`?Ut51yB6y@0k2VO@{V+ay?v`7Kdpf0X30Gv# zL&V1hQXmjr8~7CHv;q?h06*8yYaFzWc270?1==h|+Z3#v!C?OxB(v?B3_MtGxmw>- zX-?%nH)ZOBuQEHnk+~|+@eZmxEcclJo>GW0mY7tQ=z3(3L^=o@El;yOq2&Z{Lyy~Q z6rx{aZx`7uJ}L3HJYILI(-2#V7nm#7m-6T6Q0n-4TiL%t=fDqd@}6|Ess~D}ol$b% zC^t$EJO2&R4U;g&sm-T4G(}n4+%hVNI37Af zrzFVDv`?}TzI5T-Z_=z-4X}77RkxSC^Zk;EyFH(zqN19cD3P^~1m_G$do}S4$>kF5 zOp!7Hah{JewQxVYz**Cj^zn++bdQs?Y|=LmzccEuiT*>xp?p3ZIEkQ%PC36moiUqP zF3<&Y@Zg^LXqI{Z@onk}TA9S1FoWPYf@lgUxHE@j#3Hbgq7%C(>oQ)gS}>MctOPse z=Xh*6IG1fsvOMG?@4if-2snMaL$yY7I+fn{WN_vrRDP;h3HEp^bkkI!r6&17#6Ivr zk+tAI3G}Fu2KW5mHZ=10iZLe3#$HSox|hTv@+y&vqLAXl+yCH+9i*^%MIi(GU{fJx z*JI%p82zi>{uJODTC^!=(5GbXew2}PEh>0{0s1U_&!lAJ2-cX90J2mHO`8o}VM>KLHS4qVT|C-W64` z#d;^yzT>xBs;w=zyi=0{>$SGWk4)z~J!6>(%Gg4;ans3tDsu^kgL5KlmQndaUglE`ZHSE`K|s zYFlBcRWkEy(20?d+DB$3-2Yw!z(A~aHNba9Y}JWh+Rw^~1z0<2;FwAR-+kw^9B}az z=Y&puRF3M-rTBn3DY72wgCo&+>yT0xi>>0K+112qQL|b0(2WlU6;R5Q9Krmy8nRVc zXngqBBWS!F4+FlvyY2F%9PJjNU+iiHQYe!zp@W0=K{)JlM=-xjlR%Fp9gsfWpGk=)$}_6{ zSSr65byGP%)}1#3ltFV}8uAD=?AOniFgd3vR;2qfo?M4hoRx(Ept(<_{zs)jq^$_y zS?rRf9J3cC>SmVhVnv!g?><4EBJ)oH$K$NYP$Nu|tJf>b7yPNMnyEj1hsN*FjQk@} z41bP#VBIx*;@97HDSSeDMS)Mx{_X0JC%h3&39ZX00w&AoAm`rFXQU1~;TK8WqYf;{;0CfA}bX>m$Gv5n}|q0ADk zUs9yVZ1t-*$8O_!-0;!e5Ns{U8jaKm+6v}2CRHbM z5k32Q9ZC@eu-~DHPoOj!-(eLh)m0+7KO@b(D?=>jvoS>zIx^ zECr2vnL-YsQDK(g@5^68<#2AjZJ?2?o7kTpGUab>5G|%Ox+~Z8rOz4O)>{EqL^t*4 zDd76Xt(m(<*tREu{xTkEL3{YT*`rL+G`FYqr$CHKAjS!sR2Q@OORj$^qa1pT-MF-| zifH;3*kA9nxlFMPCMo;z^M5cHkjitBei8tL4m@F)|HThG1$x)lJ@HUpT9SgRKvHn= zq9hwce9?d2#Fh&F6}4c%=Rz%TeNSrY#{5sh-E!plTCOm7m10tx2TW$wl0xD(UXa-` z;uA2rdV%NsFh8yZlOnxJn93z0S+jo=Q zYEq2H#T@D@Bfxo6v7;Tq;u~pj%sgv9aUjqFN)XLyRtR|2^8au9q24w>|{)#Y zSMx6Bj)3m_VO|@EGeGyv;X8O!T zZD@t-HDvt(rjsj(Ym>xBmaf~J!}h+kfB#+;rEr*D?J3U?78b)y*52ZL3aV~Jcs+X53+ih>(xEW@Ak%a8#Z)329ARU}e zk0&);0?O$W?_W@9BYVw`ZR%nAl2C-tV;%ds%J#ppyyao$&l7}42si5H&x-Ap$lWDx z%4y}Vu>CBw`Kue&d)+(xCt^|cBp*Nua#=@r1;+N~5pF!R>r6gagp~CBni0Kf6k*R_}!g<95(zbp$e_1nAko&$Humu-e!fMb8&pH$wVxz;Bd6<6}<5lW&zSSkOVT$ z=#76T{c*IN?4|EZs?U6*Vb33TqpR^6-q+BAd%r@O#n==_dM VN7+euuIXP9g#W*&4$$e;{|DpGL`nbv diff --git a/frontend/static/sample_data/sample_data.csv b/frontend/static/sample_data/sample_data.csv index 64f1f98..0a2e489 100644 --- a/frontend/static/sample_data/sample_data.csv +++ b/frontend/static/sample_data/sample_data.csv @@ -1,11 +1,24 @@ -datetime,discharge,stage,uncertainty,comments -2021-11-22 10:00,0.13,0.33,0.05,"21 panels, some snow on ground" -2021-04-06 14:00,0.01,0.24,0.05,19 panels -2021-03-08 9:00,0.05,0.28,0.03,"35 panels, recent high flow event in December" -2021-02-24 13:00,0.03,0.27,0.05,"21 panels, rain overnight" -2021-01-28 14:00,0.12,0.32,0.04,"21 panels, good mmt" -2021-01-11 10:00,0.18,0.39,0.04,"21 panels, recent high flows after consistent rain. High turbidity" -2020-12-30 11:30,0.36,0.45,0.03,"24 panels, recent high flows and heavy rains" -2020-11-24 13:30,0.03,0.26,0.04,21 panels. Rain overnight -2021-12-13 8:00,0.29,0.4,0.04,20 panels -2020-01-12 17:00,0.67,0.55,0.03,"23 panels, recent rain on snow event, creek rising during mmt" +datetime,stage,discharge,Uncertainty,Comments +2024-01-26 11:16:54,0.502,0.341,0.064,"clear, control submerged, wading MMT 24 panels, max 6.4%" +2023-10-13 11:30:00,0.198,0.005,,"clear, section control, Volumetric bucket MMT x3" +2023-08-14 14:07:05,0.184,0.002,0.17,"clear, section control, wading MMT 13 panels, min spacing, max 18.9%" +2023-08-14 13:18:22,0.184,0.002,,"clear, section control, Volumetric bucket MMT x3" +2023-07-26 12:57:30,0.186,0.003,,"clear, section control, Volumetric bucket MMT x3, approx offset 0.170m, breakpoint at 0.386m" +2023-05-26 13:02:30,0.198,0.005,,"clear, section control, Volumetric bucket MMT x3" +2023-04-19 13:51:15,0.297,0.039,0.136,"clear, section control, wading MMT 20 panels, min spacing, max 11.9%" +2023-03-03 15:42:40,0.442,0.162,0.074,"clear, control submerged, wading MMT 24 panels, max 9.2%" +2023-02-12 15:15:00,0.368,0.067,0.074,"clear, control almost submerged, wading MMT 17 panels, max 11.9%" +2022-11-04 13:00:00,0.24,0.017,0.16,"clear, section control, Volumetric bucket MMT x3" +2022-10-01 10:00,0.285,0.002,,"debris at outlet (cleared), section control, volumetric bucket MMT x5" +2022-09-19 16:00:00,0.188,0.003,0.16,"clear, section control, Volumetric bucket MMT x3" +2022-08-17 10:00:00,0.19,0.003,0.1,"clear, section control, Volumetric bucket MMT x4" +2017-11-27 13:30:00,0.54,0.46,0.72,"clear, control submerged, wading MMT 25 panels, max 6.3%" +2014-12-19 11:30:00,0.498,0.288,,"clear, control submerged, wading MMT 21 panels, max 10.5%" +2014-02-17 13:00:00,0.658,1.06,0.16,"clear, control submerged, panels unknown" +,,,, +,,,, +,,,, +,,,, +,,,, +,,,, +,,,,