{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# MHKiT Loads Module\n",
"The following example will help familiarize you with some of the functions in the [MHKiT loads module](https://mhkit-software.github.io/MHKiT/mhkit-python/api.loads.html) that you can use to assist you in your loads analysis. \n",
"\n",
"Start by importing the necessary python packages and MHKiT module."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import pandas as pd \n",
"import numpy as np \n",
"from mhkit import utils\n",
"from mhkit import loads \n",
"import matplotlib.pyplot as plt "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Import Loads Data \n",
"\n",
"This example data comes from a land based wind turbine, since there is limited availablality of loads data for MHK devices. The example uses a subset of data from a database of 331 files, each containing 10 minutes of data sampled at 50Hz. \n",
"\n",
"As a start, lets look at the data for one of these files to figure out what formatting we need to apply. We utilize pandas to read in the csv and peek into what the data looks like. "
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
Timestamp
\n",
"
Time
\n",
"
uWind_80m
\n",
"
WD_ModActive
\n",
"
WD_Nacelle
\n",
"
WD_NacelleMod
\n",
"
LSSDW_Tq
\n",
"
LSSDW_My
\n",
"
LSSDW_Mz
\n",
"
TTTq
\n",
"
TT_ForeAft
\n",
"
TT_SideSide
\n",
"
TB_ForeAft
\n",
"
TB_SideSide
\n",
"
BL3_FlapMom
\n",
"
BL3_EdgeMom
\n",
"
BL1_FlapMom
\n",
"
BL1_EdgeMom
\n",
"
ActivePower
\n",
"
yawoffset
\n",
"
\n",
" \n",
" \n",
"
\n",
"
0
\n",
"
42795.061586
\n",
"
0.00
\n",
"
3.226754
\n",
"
1.0
\n",
"
157.302829
\n",
"
157.279582
\n",
"
-41.380694
\n",
"
-234.487436
\n",
"
-6.207381
\n",
"
-70.130726
\n",
"
-936.247028
\n",
"
-12.605151
\n",
"
-330.410779
\n",
"
1024.816867
\n",
"
470.774738
\n",
"
-165.541786
\n",
"
33.427748
\n",
"
-59.452360
\n",
"
-5.279680
\n",
"
0.023247
\n",
"
\n",
"
\n",
"
1
\n",
"
42795.061586
\n",
"
0.02
\n",
"
3.221099
\n",
"
1.0
\n",
"
157.302829
\n",
"
157.279582
\n",
"
-38.614459
\n",
"
-233.715870
\n",
"
-8.886200
\n",
"
-66.916338
\n",
"
-942.675906
\n",
"
-24.350452
\n",
"
-315.445562
\n",
"
873.214212
\n",
"
469.244736
\n",
"
-163.588005
\n",
"
32.697822
\n",
"
-62.300637
\n",
"
-5.671178
\n",
"
0.023247
\n",
"
\n",
"
\n",
"
2
\n",
"
42795.061586
\n",
"
0.04
\n",
"
3.223492
\n",
"
1.0
\n",
"
157.302829
\n",
"
157.279582
\n",
"
-39.717967
\n",
"
-234.341966
\n",
"
-7.970862
\n",
"
-67.860011
\n",
"
-922.971018
\n",
"
-22.485796
\n",
"
-292.252115
\n",
"
876.299461
\n",
"
468.736474
\n",
"
-166.018111
\n",
"
35.495810
\n",
"
-61.733604
\n",
"
-5.551847
\n",
"
0.023247
\n",
"
\n",
"
\n",
"
3
\n",
"
42795.061586
\n",
"
0.06
\n",
"
3.223274
\n",
"
1.0
\n",
"
157.302829
\n",
"
157.279582
\n",
"
-41.415143
\n",
"
-235.645598
\n",
"
-10.451794
\n",
"
-72.371951
\n",
"
-939.515265
\n",
"
-33.030892
\n",
"
-274.812769
\n",
"
763.833828
\n",
"
467.373452
\n",
"
-164.645639
\n",
"
37.952455
\n",
"
-64.390050
\n",
"
-4.626557
\n",
"
0.023247
\n",
"
\n",
"
\n",
"
4
\n",
"
42795.061587
\n",
"
0.08
\n",
"
3.223927
\n",
"
1.0
\n",
"
157.302829
\n",
"
157.279582
\n",
"
-38.614459
\n",
"
-234.755991
\n",
"
-8.648988
\n",
"
-76.530014
\n",
"
-924.771486
\n",
"
-29.228398
\n",
"
-310.213400
\n",
"
704.537757
\n",
"
466.318754
\n",
"
-161.233863
\n",
"
37.430668
\n",
"
-65.974766
\n",
"
-4.708621
\n",
"
0.023247
\n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" Timestamp Time uWind_80m WD_ModActive WD_Nacelle WD_NacelleMod \\\n",
"0 42795.061586 0.00 3.226754 1.0 157.302829 157.279582 \n",
"1 42795.061586 0.02 3.221099 1.0 157.302829 157.279582 \n",
"2 42795.061586 0.04 3.223492 1.0 157.302829 157.279582 \n",
"3 42795.061586 0.06 3.223274 1.0 157.302829 157.279582 \n",
"4 42795.061587 0.08 3.223927 1.0 157.302829 157.279582 \n",
"\n",
" LSSDW_Tq LSSDW_My LSSDW_Mz TTTq TT_ForeAft TT_SideSide \\\n",
"0 -41.380694 -234.487436 -6.207381 -70.130726 -936.247028 -12.605151 \n",
"1 -38.614459 -233.715870 -8.886200 -66.916338 -942.675906 -24.350452 \n",
"2 -39.717967 -234.341966 -7.970862 -67.860011 -922.971018 -22.485796 \n",
"3 -41.415143 -235.645598 -10.451794 -72.371951 -939.515265 -33.030892 \n",
"4 -38.614459 -234.755991 -8.648988 -76.530014 -924.771486 -29.228398 \n",
"\n",
" TB_ForeAft TB_SideSide BL3_FlapMom BL3_EdgeMom BL1_FlapMom \\\n",
"0 -330.410779 1024.816867 470.774738 -165.541786 33.427748 \n",
"1 -315.445562 873.214212 469.244736 -163.588005 32.697822 \n",
"2 -292.252115 876.299461 468.736474 -166.018111 35.495810 \n",
"3 -274.812769 763.833828 467.373452 -164.645639 37.952455 \n",
"4 -310.213400 704.537757 466.318754 -161.233863 37.430668 \n",
"\n",
" BL1_EdgeMom ActivePower yawoffset \n",
"0 -59.452360 -5.279680 0.023247 \n",
"1 -62.300637 -5.671178 0.023247 \n",
"2 -61.733604 -5.551847 0.023247 \n",
"3 -64.390050 -4.626557 0.023247 \n",
"4 -65.974766 -4.708621 0.023247 "
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"loads_data_file = './data/loads/data_loads_example.csv'\n",
"\n",
"# Import csv data file\n",
"raw_loads_data = pd.read_csv(loads_data_file)\n",
"raw_loads_data.head()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Format Loads Data with datetime index\n",
"\n",
"To use MHKiT it is important to index your DataFrame by datetime. The example loads data has two references to time, but neither are in the right format. The `'Timestamp'` column is what will give us the datetime index that we need, but we first need to convert it from microsoft excel format to pd.Datetime. "
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
uWind_80m
\n",
"
WD_ModActive
\n",
"
WD_Nacelle
\n",
"
WD_NacelleMod
\n",
"
LSSDW_Tq
\n",
"
LSSDW_My
\n",
"
LSSDW_Mz
\n",
"
TTTq
\n",
"
TT_ForeAft
\n",
"
TT_SideSide
\n",
"
TB_ForeAft
\n",
"
TB_SideSide
\n",
"
BL3_FlapMom
\n",
"
BL3_EdgeMom
\n",
"
BL1_FlapMom
\n",
"
BL1_EdgeMom
\n",
"
ActivePower
\n",
"
yawoffset
\n",
"
\n",
"
\n",
"
Timestamp
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
" \n",
" \n",
"
\n",
"
2017-03-01 01:28:40.999987200
\n",
"
3.226754
\n",
"
1.0
\n",
"
157.302829
\n",
"
157.279582
\n",
"
-41.380694
\n",
"
-234.487436
\n",
"
-6.207381
\n",
"
-70.130726
\n",
"
-936.247028
\n",
"
-12.605151
\n",
"
-330.410779
\n",
"
1024.816867
\n",
"
470.774738
\n",
"
-165.541786
\n",
"
33.427748
\n",
"
-59.452360
\n",
"
-5.279680
\n",
"
0.023247
\n",
"
\n",
"
\n",
"
2017-03-01 01:28:41.020032000
\n",
"
3.221099
\n",
"
1.0
\n",
"
157.302829
\n",
"
157.279582
\n",
"
-38.614459
\n",
"
-233.715870
\n",
"
-8.886200
\n",
"
-66.916338
\n",
"
-942.675906
\n",
"
-24.350452
\n",
"
-315.445562
\n",
"
873.214212
\n",
"
469.244736
\n",
"
-163.588005
\n",
"
32.697822
\n",
"
-62.300637
\n",
"
-5.671178
\n",
"
0.023247
\n",
"
\n",
"
\n",
"
2017-03-01 01:28:41.039990400
\n",
"
3.223492
\n",
"
1.0
\n",
"
157.302829
\n",
"
157.279582
\n",
"
-39.717967
\n",
"
-234.341966
\n",
"
-7.970862
\n",
"
-67.860011
\n",
"
-922.971018
\n",
"
-22.485796
\n",
"
-292.252115
\n",
"
876.299461
\n",
"
468.736474
\n",
"
-166.018111
\n",
"
35.495810
\n",
"
-61.733604
\n",
"
-5.551847
\n",
"
0.023247
\n",
"
\n",
"
\n",
"
2017-03-01 01:28:41.060035200
\n",
"
3.223274
\n",
"
1.0
\n",
"
157.302829
\n",
"
157.279582
\n",
"
-41.415143
\n",
"
-235.645598
\n",
"
-10.451794
\n",
"
-72.371951
\n",
"
-939.515265
\n",
"
-33.030892
\n",
"
-274.812769
\n",
"
763.833828
\n",
"
467.373452
\n",
"
-164.645639
\n",
"
37.952455
\n",
"
-64.390050
\n",
"
-4.626557
\n",
"
0.023247
\n",
"
\n",
"
\n",
"
2017-03-01 01:28:41.079993600
\n",
"
3.223927
\n",
"
1.0
\n",
"
157.302829
\n",
"
157.279582
\n",
"
-38.614459
\n",
"
-234.755991
\n",
"
-8.648988
\n",
"
-76.530014
\n",
"
-924.771486
\n",
"
-29.228398
\n",
"
-310.213400
\n",
"
704.537757
\n",
"
466.318754
\n",
"
-161.233863
\n",
"
37.430668
\n",
"
-65.974766
\n",
"
-4.708621
\n",
"
0.023247
\n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" uWind_80m WD_ModActive WD_Nacelle \\\n",
"Timestamp \n",
"2017-03-01 01:28:40.999987200 3.226754 1.0 157.302829 \n",
"2017-03-01 01:28:41.020032000 3.221099 1.0 157.302829 \n",
"2017-03-01 01:28:41.039990400 3.223492 1.0 157.302829 \n",
"2017-03-01 01:28:41.060035200 3.223274 1.0 157.302829 \n",
"2017-03-01 01:28:41.079993600 3.223927 1.0 157.302829 \n",
"\n",
" WD_NacelleMod LSSDW_Tq LSSDW_My \\\n",
"Timestamp \n",
"2017-03-01 01:28:40.999987200 157.279582 -41.380694 -234.487436 \n",
"2017-03-01 01:28:41.020032000 157.279582 -38.614459 -233.715870 \n",
"2017-03-01 01:28:41.039990400 157.279582 -39.717967 -234.341966 \n",
"2017-03-01 01:28:41.060035200 157.279582 -41.415143 -235.645598 \n",
"2017-03-01 01:28:41.079993600 157.279582 -38.614459 -234.755991 \n",
"\n",
" LSSDW_Mz TTTq TT_ForeAft TT_SideSide \\\n",
"Timestamp \n",
"2017-03-01 01:28:40.999987200 -6.207381 -70.130726 -936.247028 -12.605151 \n",
"2017-03-01 01:28:41.020032000 -8.886200 -66.916338 -942.675906 -24.350452 \n",
"2017-03-01 01:28:41.039990400 -7.970862 -67.860011 -922.971018 -22.485796 \n",
"2017-03-01 01:28:41.060035200 -10.451794 -72.371951 -939.515265 -33.030892 \n",
"2017-03-01 01:28:41.079993600 -8.648988 -76.530014 -924.771486 -29.228398 \n",
"\n",
" TB_ForeAft TB_SideSide BL3_FlapMom \\\n",
"Timestamp \n",
"2017-03-01 01:28:40.999987200 -330.410779 1024.816867 470.774738 \n",
"2017-03-01 01:28:41.020032000 -315.445562 873.214212 469.244736 \n",
"2017-03-01 01:28:41.039990400 -292.252115 876.299461 468.736474 \n",
"2017-03-01 01:28:41.060035200 -274.812769 763.833828 467.373452 \n",
"2017-03-01 01:28:41.079993600 -310.213400 704.537757 466.318754 \n",
"\n",
" BL3_EdgeMom BL1_FlapMom BL1_EdgeMom \\\n",
"Timestamp \n",
"2017-03-01 01:28:40.999987200 -165.541786 33.427748 -59.452360 \n",
"2017-03-01 01:28:41.020032000 -163.588005 32.697822 -62.300637 \n",
"2017-03-01 01:28:41.039990400 -166.018111 35.495810 -61.733604 \n",
"2017-03-01 01:28:41.060035200 -164.645639 37.952455 -64.390050 \n",
"2017-03-01 01:28:41.079993600 -161.233863 37.430668 -65.974766 \n",
"\n",
" ActivePower yawoffset \n",
"Timestamp \n",
"2017-03-01 01:28:40.999987200 -5.279680 0.023247 \n",
"2017-03-01 01:28:41.020032000 -5.671178 0.023247 \n",
"2017-03-01 01:28:41.039990400 -5.551847 0.023247 \n",
"2017-03-01 01:28:41.060035200 -4.626557 0.023247 \n",
"2017-03-01 01:28:41.079993600 -4.708621 0.023247 "
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Use the datetime conversion from the utils module\n",
"datetime = utils.excel_to_datetime(raw_loads_data['Timestamp'])\n",
"\n",
"# Replace the 'Timestamp' column with our newly formatted datetime\n",
"raw_loads_data['Timestamp'] = datetime \n",
"\n",
"# Set this as our index for our DataFrame\n",
"loads_data = raw_loads_data.set_index('Timestamp')\n",
"\n",
"# Remove the 'time' column since it will not be used\n",
"loads_data.drop(columns='Time',inplace=True)\n",
"loads_data.head()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Loads Analysis\n",
"Now that we have our loads data in the correct format for MHKiT, we do some analysis.\n",
"\n",
"### Damage Equivalent Loads\n",
"Lets say that we wanted to investigate fatigue. We can do this by calculating short-term damage equivalent loads (DELs). In this instance, we calculate the DELs on the tower base moment `'TB_ForeAft'`, and on blade 1 root flap moment `'BL1_FlapMom'`. Our tower is steel while our blade is composite so they will have different material slopes. We will run the function `damage_equivalent_load` on each data signal. \n",
"\n",
"We call our function and apply the default inputs of using at least 100 bins for the load ranges and we let data_length=600 seconds so that we get an equivalent 1Hz DEL for our 10 minute file. "
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"tags": []
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"DEL TB_ForeAft: 3912.6390862773\n",
"DEL BL1_FlapMom: 1435.8222478714554\n"
]
}
],
"source": [
"# Calculate the damage equivalent load for blade 1 root momement and tower base moment\n",
"DEL_tower = loads.general.damage_equivalent_load(loads_data['TB_ForeAft'],4,\n",
" bin_num=100,data_length=600)\n",
"DEL_blade = loads.general.damage_equivalent_load(loads_data['BL1_FlapMom'],10,\n",
" bin_num=100,data_length=600)\n",
"print('DEL TB_ForeAft: '+ str(DEL_tower))\n",
"print('DEL BL1_FlapMom: '+ str(DEL_blade))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Calculate Statistics\n",
"\n",
"Another important part of loads analysis is looking at statistics. Here, we use another function to help us calculate the mean, max, min, and std for this 10 minute file. Per standards, a valid statistical window has to be consecutive in time with the correct number of datapoints. If this 10 minute file did meet this criteria, then no stats would be generated and a warning message would appear. \n",
"\n",
"NOTE: Sometimes individual files may contain enough data for multiple statistical windows. This function can still handle this scenario as long as the correct inputs are specified. "
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
uWind_80m
\n",
"
WD_ModActive
\n",
"
WD_Nacelle
\n",
"
WD_NacelleMod
\n",
"
LSSDW_Tq
\n",
"
LSSDW_My
\n",
"
LSSDW_Mz
\n",
"
TTTq
\n",
"
TT_ForeAft
\n",
"
TT_SideSide
\n",
"
TB_ForeAft
\n",
"
TB_SideSide
\n",
"
BL3_FlapMom
\n",
"
BL3_EdgeMom
\n",
"
BL1_FlapMom
\n",
"
BL1_EdgeMom
\n",
"
ActivePower
\n",
"
yawoffset
\n",
"
\n",
" \n",
" \n",
"
\n",
"
2017-03-01 01:28:41
\n",
"
7.773325
\n",
"
1.0
\n",
"
178.612256
\n",
"
178.602595
\n",
"
127.244191
\n",
"
-252.23813
\n",
"
3.50322
\n",
"
7.032573
\n",
"
-846.663367
\n",
"
271.446574
\n",
"
3785.034515
\n",
"
7.199176
\n",
"
-494.858287
\n",
"
266.790368
\n",
"
-452.652744
\n",
"
21.259999
\n",
"
234.578289
\n",
"
0.009661
\n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" uWind_80m WD_ModActive WD_Nacelle WD_NacelleMod \\\n",
"2017-03-01 01:28:41 7.773325 1.0 178.612256 178.602595 \n",
"\n",
" LSSDW_Tq LSSDW_My LSSDW_Mz TTTq TT_ForeAft \\\n",
"2017-03-01 01:28:41 127.244191 -252.23813 3.50322 7.032573 -846.663367 \n",
"\n",
" TT_SideSide TB_ForeAft TB_SideSide BL3_FlapMom \\\n",
"2017-03-01 01:28:41 271.446574 3785.034515 7.199176 -494.858287 \n",
"\n",
" BL3_EdgeMom BL1_FlapMom BL1_EdgeMom ActivePower \\\n",
"2017-03-01 01:28:41 266.790368 -452.652744 21.259999 234.578289 \n",
"\n",
" yawoffset \n",
"2017-03-01 01:28:41 0.009661 "
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Calculate the means, maxs, mins, and stdevs for all data signals in the loads data file\n",
"means,maxs,mins,stdevs = utils.get_statistics(loads_data,50,period=600)\n",
"\n",
"# Display the results, indexed by the first timestamp of the corresponding statistical window\n",
"means"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"At this point, it would be nice to start visualizing some of this data. In order to do this, we need to calculate the stats and DELs for all the files in our database. In this case, it would be done through a loop that imports each file and applies all the functions we just saw. At the end of the loop, we store the result by appending to a list (this is computationally more efficient than concatenating dataframes). Finally, we can convert our lists to dataframes so that its easier to play with the data. To speed things up, this was already done so we just need to import the resulting dataframes. \n",
"\n",
"As a reference, an example code block is shown of how to create a loop that calculates all the means for each file which are then stored into a dataframe. \n",
"\n",
"```python\n",
"# pre-allocate lists for storage\n",
"means = []\n",
"time = []\n",
"\n",
"# start loop\n",
"for f in os.listdir(pathOut):\n",
" if f.endswith('.csv'):\n",
" # import csv file\n",
" raw_loads_data = pd.read_csv(pathOut+'/'+f)\n",
" # replace the timestamp column with formatted datetime\n",
" datetime = utils.excel_to_datetime(raw_loads_data['Timestamp']) \n",
" raw_loads_data['Timestamp'] = datetime \n",
" # set this as our index for our dataframe\n",
" loads_data = raw_loads_data.set_index('Timestamp')\n",
" # remove the \"time\" column as its unnecessary \n",
" loads_data.drop(columns='Time',inplace=True)\n",
" # get stats\n",
" fmean, fmax, fmin, fstd = utils.get_statistics(loads_data,freq=50,period=600)\n",
" means.append(fmean.values.tolist())\n",
" time.append(fmean.index.values)\n",
"\n",
"# convert lists into a dataframe\n",
"loads_means = pd.DataFrame(np.squeeze(means),columns=loads_data.columns.values,index=time)\n",
"```"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
"text/plain": [
""
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"loads.graphics.plot_statistics(means['uWind_80m'],\n",
" means['BL1_FlapMom'],\n",
" maxs['BL1_FlapMom'],\n",
" mins['BL1_FlapMom'],\n",
" y_stdev=std['BL1_FlapMom'],\n",
" xlabel='Wind Speed [m/s]',\n",
" ylabel='Blade Flap Moment [kNm]',\n",
" title = 'Blade Flap Moment Load Statistics')\n",
"\n",
"loads.graphics.plot_statistics(means['uWind_80m'],\n",
" means['TB_ForeAft'],\n",
" maxs['TB_ForeAft'],\n",
" mins['TB_ForeAft'],\n",
" y_stdev=std['TB_ForeAft'],\n",
" xlabel='Wind Speed [m/s]',\n",
" ylabel='Tower Base Moment [kNm]',\n",
" title = 'Tower Base Moment Load Statistics')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Another common step is to bin the statistical data. This can easily be done with the bin_stats function from the loads module shown below. A warning message will show if there are any bins that were not filled. "
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"tags": []
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Warning: some bins may be empty!\n",
"Warning: some bins may be empty!\n",
"Warning: some bins may be empty!\n"
]
},
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
uWind_80m
\n",
"
WD_ModActive
\n",
"
WD_Nacelle
\n",
"
WD_NacelleMod
\n",
"
LSSDW_Tq
\n",
"
LSSDW_My
\n",
"
LSSDW_Mz
\n",
"
TTTq
\n",
"
TT_ForeAft
\n",
"
TT_SideSide
\n",
"
TB_ForeAft
\n",
"
TB_SideSide
\n",
"
BL3_FlapMom
\n",
"
BL3_EdgeMom
\n",
"
BL1_FlapMom
\n",
"
BL1_EdgeMom
\n",
"
ActivePower
\n",
"
yawoffset
\n",
"
\n",
" \n",
" \n",
"
\n",
"
0
\n",
"
3.582788
\n",
"
1.000000
\n",
"
179.116583
\n",
"
179.215605
\n",
"
29.890756
\n",
"
-252.235192
\n",
"
-2.931479
\n",
"
-26.271654
\n",
"
-932.953738
\n",
"
126.807011
\n",
"
994.638671
\n",
"
-102.050137
\n",
"
-206.591291
\n",
"
129.390923
\n",
"
-238.085263
\n",
"
-14.874877
\n",
"
21.370858
\n",
"
-0.099022
\n",
"
\n",
"
\n",
"
1
\n",
"
4.498995
\n",
"
1.000000
\n",
"
178.968041
\n",
"
178.963826
\n",
"
77.559257
\n",
"
-248.081262
\n",
"
-3.389182
\n",
"
-51.890684
\n",
"
-856.012087
\n",
"
154.470046
\n",
"
2235.672810
\n",
"
55.864708
\n",
"
-327.707803
\n",
"
162.318746
\n",
"
-354.474414
\n",
"
-7.655706
\n",
"
76.033527
\n",
"
0.004216
\n",
"
\n",
"
\n",
"
2
\n",
"
5.525786
\n",
"
0.980392
\n",
"
180.016225
\n",
"
179.996997
\n",
"
171.459516
\n",
"
-230.697510
\n",
"
-2.588070
\n",
"
-50.333653
\n",
"
-761.047719
\n",
"
213.700239
\n",
"
3822.643018
\n",
"
137.424215
\n",
"
-495.457116
\n",
"
191.503291
\n",
"
-599.658658
\n",
"
1436.266991
\n",
"
167.928433
\n",
"
0.019228
\n",
"
\n",
"
\n",
"
3
\n",
"
6.520641
\n",
"
0.970588
\n",
"
180.835469
\n",
"
180.788600
\n",
"
219.138550
\n",
"
-237.979417
\n",
"
-5.688640
\n",
"
-42.827109
\n",
"
-647.310249
\n",
"
302.591017
\n",
"
6122.425105
\n",
"
192.772384
\n",
"
-708.321569
\n",
"
237.072176
\n",
"
-777.545681
\n",
"
599.754365
\n",
"
317.629187
\n",
"
0.046869
\n",
"
\n",
"
\n",
"
4
\n",
"
7.534458
\n",
"
1.000000
\n",
"
181.373937
\n",
"
181.474371
\n",
"
285.072650
\n",
"
-240.838000
\n",
"
-8.602511
\n",
"
-56.222891
\n",
"
-571.093365
\n",
"
376.755167
\n",
"
7551.899048
\n",
"
296.723677
\n",
"
-861.049383
\n",
"
252.077699
\n",
"
-874.336546
\n",
"
88.957621
\n",
"
454.388684
\n",
"
-0.100434
\n",
"
\n",
"
\n",
"
5
\n",
"
8.483451
\n",
"
0.966667
\n",
"
180.006753
\n",
"
180.002341
\n",
"
405.504595
\n",
"
-233.796744
\n",
"
-7.596248
\n",
"
-32.233513
\n",
"
-430.681550
\n",
"
502.070942
\n",
"
9938.523640
\n",
"
335.667789
\n",
"
-1073.556266
\n",
"
313.335072
\n",
"
-1167.966930
\n",
"
1302.506123
\n",
"
679.886822
\n",
"
0.004413
\n",
"
\n",
"
\n",
"
6
\n",
"
9.635533
\n",
"
0.880000
\n",
"
178.894596
\n",
"
178.885289
\n",
"
581.709245
\n",
"
-218.693601
\n",
"
-8.801577
\n",
"
-72.657764
\n",
"
-330.525678
\n",
"
668.243228
\n",
"
11477.966048
\n",
"
588.384954
\n",
"
-1209.842191
\n",
"
342.343342
\n",
"
-1526.574097
\n",
"
4686.436750
\n",
"
942.635185
\n",
"
0.009307
\n",
"
\n",
"
\n",
"
7
\n",
"
10.552790
\n",
"
0.892857
\n",
"
179.473587
\n",
"
179.487138
\n",
"
653.097969
\n",
"
-224.998487
\n",
"
-5.756869
\n",
"
-62.714160
\n",
"
-283.620699
\n",
"
750.186980
\n",
"
12002.102796
\n",
"
740.834502
\n",
"
-1234.723513
\n",
"
354.710593
\n",
"
-1380.119025
\n",
"
2065.459821
\n",
"
1101.930764
\n",
"
-0.013551
\n",
"
\n",
"
\n",
"
8
\n",
"
11.431282
\n",
"
0.814807
\n",
"
179.157452
\n",
"
179.173400
\n",
"
738.418861
\n",
"
-209.685631
\n",
"
-1.320960
\n",
"
-22.419344
\n",
"
-295.116952
\n",
"
874.826602
\n",
"
12237.843754
\n",
"
668.581523
\n",
"
-1248.802900
\n",
"
362.942457
\n",
"
-1291.630102
\n",
"
173.010056
\n",
"
1251.222715
\n",
"
-0.015947
\n",
"
\n",
"
\n",
"
9
\n",
"
12.447532
\n",
"
0.888889
\n",
"
178.668483
\n",
"
178.658239
\n",
"
791.696059
\n",
"
-188.741923
\n",
"
-2.588168
\n",
"
-29.769987
\n",
"
-291.021580
\n",
"
902.150404
\n",
"
11648.909089
\n",
"
743.454123
\n",
"
-1173.606189
\n",
"
328.218911
\n",
"
-1262.103186
\n",
"
738.232546
\n",
"
1316.246001
\n",
"
0.010243
\n",
"
\n",
"
\n",
"
10
\n",
"
13.406650
\n",
"
0.833333
\n",
"
180.438968
\n",
"
180.452124
\n",
"
817.388823
\n",
"
-198.989294
\n",
"
-7.197755
\n",
"
-24.899471
\n",
"
-331.564307
\n",
"
992.711011
\n",
"
10733.428547
\n",
"
791.711119
\n",
"
-1074.877395
\n",
"
300.227189
\n",
"
-1114.404530
\n",
"
121.228715
\n",
"
1384.312354
\n",
"
-0.013155
\n",
"
\n",
"
\n",
"
11
\n",
"
14.348479
\n",
"
0.542676
\n",
"
179.354862
\n",
"
179.346215
\n",
"
875.198293
\n",
"
-59.267478
\n",
"
-1.437998
\n",
"
4.153358
\n",
"
-374.601981
\n",
"
907.552054
\n",
"
9896.827494
\n",
"
677.680004
\n",
"
-971.840066
\n",
"
285.582399
\n",
"
-1036.798595
\n",
"
102.258425
\n",
"
1350.838415
\n",
"
0.008648
\n",
"
\n",
"
\n",
"
12
\n",
"
15.579021
\n",
"
0.833333
\n",
"
178.013520
\n",
"
178.001359
\n",
"
880.603457
\n",
"
-164.536103
\n",
"
-1.851375
\n",
"
-26.848873
\n",
"
-330.960964
\n",
"
1059.024803
\n",
"
9234.857465
\n",
"
926.368719
\n",
"
-899.553762
\n",
"
227.524012
\n",
"
-923.851217
\n",
"
74.708350
\n",
"
1491.638595
\n",
"
0.012160
\n",
"
\n",
"
\n",
"
13
\n",
"
16.400735
\n",
"
1.000000
\n",
"
178.314956
\n",
"
178.311169
\n",
"
866.227371
\n",
"
-185.673240
\n",
"
-12.999453
\n",
"
-78.397677
\n",
"
-287.228505
\n",
"
1071.363839
\n",
"
8983.182063
\n",
"
984.403585
\n",
"
-860.376373
\n",
"
200.755492
\n",
"
-864.637118
\n",
"
53.896186
\n",
"
1499.435045
\n",
"
0.003787
\n",
"
\n",
"
\n",
"
14
\n",
"
17.602406
\n",
"
1.000000
\n",
"
179.171230
\n",
"
179.157532
\n",
"
886.425220
\n",
"
-137.146930
\n",
"
-8.231045
\n",
"
-3.156350
\n",
"
-383.416299
\n",
"
1063.362197
\n",
"
8563.213120
\n",
"
873.449181
\n",
"
-762.282362
\n",
"
225.211156
\n",
"
-810.106716
\n",
"
48.225543
\n",
"
1476.244726
\n",
"
0.013698
\n",
"
\n",
"
\n",
"
15
\n",
"
18.982808
\n",
"
1.000000
\n",
"
179.199942
\n",
"
179.421965
\n",
"
863.882390
\n",
"
-158.281553
\n",
"
-23.757701
\n",
"
-74.334173
\n",
"
-270.882688
\n",
"
1056.399578
\n",
"
8206.587839
\n",
"
1349.955119
\n",
"
-699.125685
\n",
"
205.630765
\n",
"
-703.231687
\n",
"
14.362069
\n",
"
1514.168487
\n",
"
-0.222023
\n",
"
\n",
"
\n",
"
16
\n",
"
19.519921
\n",
"
1.000000
\n",
"
178.361870
\n",
"
178.356581
\n",
"
881.739431
\n",
"
-169.696939
\n",
"
1.018761
\n",
"
-61.347061
\n",
"
-284.548611
\n",
"
1113.928290
\n",
"
7807.792946
\n",
"
981.873207
\n",
"
-695.129807
\n",
"
151.699529
\n",
"
-709.755338
\n",
"
20.152316
\n",
"
1514.549099
\n",
"
0.005289
\n",
"
\n",
"
\n",
"
17
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
\n",
"
\n",
"
18
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
\n",
"
\n",
"
19
\n",
"
22.062534
\n",
"
1.000000
\n",
"
178.524018
\n",
"
178.507149
\n",
"
854.648338
\n",
"
-142.928299
\n",
"
-43.759922
\n",
"
-111.305482
\n",
"
-262.456906
\n",
"
1067.274340
\n",
"
7566.500039
\n",
"
1369.827084
\n",
"
-572.629544
\n",
"
186.356598
\n",
"
-575.251752
\n",
"
4.025634
\n",
"
1500.616164
\n",
"
0.016870
\n",
"
\n",
"
\n",
"
20
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
NaN
\n",
"
\n",
"
\n",
"
21
\n",
"
24.759699
\n",
"
1.000000
\n",
"
179.912515
\n",
"
179.904301
\n",
"
334.771556
\n",
"
-206.441986
\n",
"
-11.694631
\n",
"
-32.524912
\n",
"
-679.532690
\n",
"
492.819279
\n",
"
3419.356234
\n",
"
752.087159
\n",
"
-226.330646
\n",
"
137.392161
\n",
"
-187.876792
\n",
"
-33.225278
\n",
"
615.054202
\n",
"
0.008214
\n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" uWind_80m WD_ModActive WD_Nacelle WD_NacelleMod LSSDW_Tq \\\n",
"0 3.582788 1.000000 179.116583 179.215605 29.890756 \n",
"1 4.498995 1.000000 178.968041 178.963826 77.559257 \n",
"2 5.525786 0.980392 180.016225 179.996997 171.459516 \n",
"3 6.520641 0.970588 180.835469 180.788600 219.138550 \n",
"4 7.534458 1.000000 181.373937 181.474371 285.072650 \n",
"5 8.483451 0.966667 180.006753 180.002341 405.504595 \n",
"6 9.635533 0.880000 178.894596 178.885289 581.709245 \n",
"7 10.552790 0.892857 179.473587 179.487138 653.097969 \n",
"8 11.431282 0.814807 179.157452 179.173400 738.418861 \n",
"9 12.447532 0.888889 178.668483 178.658239 791.696059 \n",
"10 13.406650 0.833333 180.438968 180.452124 817.388823 \n",
"11 14.348479 0.542676 179.354862 179.346215 875.198293 \n",
"12 15.579021 0.833333 178.013520 178.001359 880.603457 \n",
"13 16.400735 1.000000 178.314956 178.311169 866.227371 \n",
"14 17.602406 1.000000 179.171230 179.157532 886.425220 \n",
"15 18.982808 1.000000 179.199942 179.421965 863.882390 \n",
"16 19.519921 1.000000 178.361870 178.356581 881.739431 \n",
"17 NaN NaN NaN NaN NaN \n",
"18 NaN NaN NaN NaN NaN \n",
"19 22.062534 1.000000 178.524018 178.507149 854.648338 \n",
"20 NaN NaN NaN NaN NaN \n",
"21 24.759699 1.000000 179.912515 179.904301 334.771556 \n",
"\n",
" LSSDW_My LSSDW_Mz TTTq TT_ForeAft TT_SideSide TB_ForeAft \\\n",
"0 -252.235192 -2.931479 -26.271654 -932.953738 126.807011 994.638671 \n",
"1 -248.081262 -3.389182 -51.890684 -856.012087 154.470046 2235.672810 \n",
"2 -230.697510 -2.588070 -50.333653 -761.047719 213.700239 3822.643018 \n",
"3 -237.979417 -5.688640 -42.827109 -647.310249 302.591017 6122.425105 \n",
"4 -240.838000 -8.602511 -56.222891 -571.093365 376.755167 7551.899048 \n",
"5 -233.796744 -7.596248 -32.233513 -430.681550 502.070942 9938.523640 \n",
"6 -218.693601 -8.801577 -72.657764 -330.525678 668.243228 11477.966048 \n",
"7 -224.998487 -5.756869 -62.714160 -283.620699 750.186980 12002.102796 \n",
"8 -209.685631 -1.320960 -22.419344 -295.116952 874.826602 12237.843754 \n",
"9 -188.741923 -2.588168 -29.769987 -291.021580 902.150404 11648.909089 \n",
"10 -198.989294 -7.197755 -24.899471 -331.564307 992.711011 10733.428547 \n",
"11 -59.267478 -1.437998 4.153358 -374.601981 907.552054 9896.827494 \n",
"12 -164.536103 -1.851375 -26.848873 -330.960964 1059.024803 9234.857465 \n",
"13 -185.673240 -12.999453 -78.397677 -287.228505 1071.363839 8983.182063 \n",
"14 -137.146930 -8.231045 -3.156350 -383.416299 1063.362197 8563.213120 \n",
"15 -158.281553 -23.757701 -74.334173 -270.882688 1056.399578 8206.587839 \n",
"16 -169.696939 1.018761 -61.347061 -284.548611 1113.928290 7807.792946 \n",
"17 NaN NaN NaN NaN NaN NaN \n",
"18 NaN NaN NaN NaN NaN NaN \n",
"19 -142.928299 -43.759922 -111.305482 -262.456906 1067.274340 7566.500039 \n",
"20 NaN NaN NaN NaN NaN NaN \n",
"21 -206.441986 -11.694631 -32.524912 -679.532690 492.819279 3419.356234 \n",
"\n",
" TB_SideSide BL3_FlapMom BL3_EdgeMom BL1_FlapMom BL1_EdgeMom \\\n",
"0 -102.050137 -206.591291 129.390923 -238.085263 -14.874877 \n",
"1 55.864708 -327.707803 162.318746 -354.474414 -7.655706 \n",
"2 137.424215 -495.457116 191.503291 -599.658658 1436.266991 \n",
"3 192.772384 -708.321569 237.072176 -777.545681 599.754365 \n",
"4 296.723677 -861.049383 252.077699 -874.336546 88.957621 \n",
"5 335.667789 -1073.556266 313.335072 -1167.966930 1302.506123 \n",
"6 588.384954 -1209.842191 342.343342 -1526.574097 4686.436750 \n",
"7 740.834502 -1234.723513 354.710593 -1380.119025 2065.459821 \n",
"8 668.581523 -1248.802900 362.942457 -1291.630102 173.010056 \n",
"9 743.454123 -1173.606189 328.218911 -1262.103186 738.232546 \n",
"10 791.711119 -1074.877395 300.227189 -1114.404530 121.228715 \n",
"11 677.680004 -971.840066 285.582399 -1036.798595 102.258425 \n",
"12 926.368719 -899.553762 227.524012 -923.851217 74.708350 \n",
"13 984.403585 -860.376373 200.755492 -864.637118 53.896186 \n",
"14 873.449181 -762.282362 225.211156 -810.106716 48.225543 \n",
"15 1349.955119 -699.125685 205.630765 -703.231687 14.362069 \n",
"16 981.873207 -695.129807 151.699529 -709.755338 20.152316 \n",
"17 NaN NaN NaN NaN NaN \n",
"18 NaN NaN NaN NaN NaN \n",
"19 1369.827084 -572.629544 186.356598 -575.251752 4.025634 \n",
"20 NaN NaN NaN NaN NaN \n",
"21 752.087159 -226.330646 137.392161 -187.876792 -33.225278 \n",
"\n",
" ActivePower yawoffset \n",
"0 21.370858 -0.099022 \n",
"1 76.033527 0.004216 \n",
"2 167.928433 0.019228 \n",
"3 317.629187 0.046869 \n",
"4 454.388684 -0.100434 \n",
"5 679.886822 0.004413 \n",
"6 942.635185 0.009307 \n",
"7 1101.930764 -0.013551 \n",
"8 1251.222715 -0.015947 \n",
"9 1316.246001 0.010243 \n",
"10 1384.312354 -0.013155 \n",
"11 1350.838415 0.008648 \n",
"12 1491.638595 0.012160 \n",
"13 1499.435045 0.003787 \n",
"14 1476.244726 0.013698 \n",
"15 1514.168487 -0.222023 \n",
"16 1514.549099 0.005289 \n",
"17 NaN NaN \n",
"18 NaN NaN \n",
"19 1500.616164 0.016870 \n",
"20 NaN NaN \n",
"21 615.054202 0.008214 "
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Create array containing wind speeds to use as bin edges\n",
"bin_edges = np.arange(3,26,1)\n",
"bin_against = means['uWind_80m']\n",
"\n",
"# Apply function for means, maxs, and mins \n",
"[bin_means, bin_means_std] = loads.general.bin_statistics(means,bin_against,bin_edges)\n",
"[bin_maxs, bin_maxs_std] = loads.general.bin_statistics(maxs,bin_against,bin_edges)\n",
"[bin_mins, bin_mins_std] = loads.general.bin_statistics(mins,bin_against,bin_edges)\n",
"\n",
"bin_means"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Now lets make some more plots with the binned data. Here we use the binned data and corresponding standard deviations as inputs to the plot_bin_statistics function. "
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAfAAAAFgCAYAAABEyiulAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABkCUlEQVR4nO3dd3gU1dfA8e9NLySQSgsQAqFXpasUEcQGqKiIFVREQP3xqmDvDUXFAoIFBAugoBRBAVGKSgu99wChpYeEkH7fP2YTAiakbTJbzud59kl2dnbm7GQ3Z29XWmuEEEIIYV9czA5ACCGEEGUnCVwIIYSwQ5LAhRBCCDskCVwIIYSwQ5LAhRBCCDskCVwIIYSwQ5LAhbARSqkpSqmXTDjvg0qpv6v6vIXOX67XrZSqr5RKU0q5VkZcQtg6SeBCVBGlVLRS6rwl6SQppRYrperlP661HqG1fsPMGIuilHpIKbVXKZWqlDqjlFqilPKzPPaNUurNMhzrP18WSvu6LdfvukLPO6a1rqa1zi3L6xHCUUgCF6Jq3aK1rgbUBs4An5ocz2UppXoAbwN3a639gObAHHOjEkKAJHAhTKG1zgDmAi3ytxUuzSqleiqlYpRSTymlYpVSp5RSQy/Zd5KlFJ+qlFqvlGpU6PFmSqnlSqlEpdQ+pdSdhR4LUkotVEqdVUptAAqeV4SOwFqt9RZL3Ila6xla61Sl1HDgHmCspVZhkeX4zyqlDlni2q2UutWyvTkwBehq2T+5iNcdrJT6VSmVbIl9jVLKRSn1LVAfWGR57lilVLhSSiul3CzPDVRKTVdKnbTUcMy/3DHL9YcTwobIm1gIEyilfIC7gHWX2a0WUB2oCzwETFJKBRR6fDDwGhAAHATeshzbF1gO/ACEWvabrJTK/7IwCcjAqAUYZrkVZz1wvVLqNaXUVUopz/wHtNZfAN8D71mqsm+xPHQIuMYS+2vAd0qp2lrrPcAIjC8E1bTWNYo431NADBAC1ASeN06l7wOOYanB0Fq/V8RzvwV8gJaW1/3R5Y55mdcshF2QBC5E1ZpvKXmmAH2A9y+zbzbwutY6W2u9BEgDmhZ6/Bet9QatdQ5GIm1n2X4zEK21nq61zrGUnucBd1g6fN0OvKy1Pqe13gnMKC4ArfUa4DbgCmAxkKCU+vByHce01j9prU9qrfO01nOAA0Cny7zOS19zbaCB5XWv0aVYsEEpVRu4ARihtU6yPHdVRY4phK2TBC5E1RpoKXl6AaOBVUqpWsXsm2BJzvnSgWqF7p8u5rEGQGdLlXGy5QvDPRgl+hDADThe6LlHLxew1vo3S+k6EBgAPAg8XNz+Sqn7lVJbC527FRB8uXMU8j5GbcIypdRhpdSzpXxePSBRa51kxWMKYdMkgQthAq11rtb6ZyAXuNrKhz8OrNJa1yh0q6a1fgyIA3IwEl6++qWMOU9rvQL4EyMpwyVV0UqpBsCXGF9OgixfVnYCqqj9izhHqtb6Ka11BNAf+D+lVO9SPPc4EKiUqlHGYwphtySBC2ECZRiA0X69x8qH/xVoopS6Tynlbrl1VEo1twy5+hl4VSnlY2kXf+AycQ5QSg1WSgVYYu4E9OBC2/0ZIKLQU3wxEm2c5flDuZDs8/cPU0p5FHO+m5VSjZVSCqOZIRfIK+ZcBbTWp4DfMNr6AyyvuXspjimE3ZIELkTVWqSUSgPOYnQ6e0BrvcuaJ9BapwJ9MTqvncSoah8P5HdAG41R3X4a+AaYfpnDJQGPYLRjnwW+A97XWn9vefxroIWluny+1no38AGwFiPhtgb+KXS8P4FdwGmlVHwR54sE/sBo718LTNZa/2V57B3gRcu5ni7iufdhtHfvBWKB/5XimELYLSV9OYQQQgj7IyVwIYQQwg5JAhdCCCHskCRwIYQQwg5JAhdCCCHskJvZAVRUcHCwDg8PNzuMcsvKysLDo8gRNaIKyPU3l1x/88i1N1dZrv+mTZvitdYhl263+wQeHh5OVFSU2WGUW3R0NPb8BcTeyfU3l1x/88i1N1dZrr9SqsjZEqUKXQghhLBDksCFEEIIOyQJXAghhLBDdt8GXpTs7GxiYmLIyMgwO5QS5eTksGeP9abC9vLyIiwsDHd3d6sdUwghhO1xyAQeExODn58f4eHhGOsX2K7MzEw8PT1L3rEUtNYkJCQQExNDw4YNrXJMIYQQtskhq9AzMjIICgqy+eRtbUopgoKC7KLmQQghRMU4ZAIHnC5553PW1y2EEM7GIavQS+uj5fv5eMWBYh9/snckY/o0qcKIhBBCiNJx6gQ+pk+TggR919S1AMx5tKuZIQkhhBCl4rBV6GWVlZPH3tOpxKZK+7EQQgjbJwnc4lTKedIyc5i66rBVjhcdHU2zZs148MEHadKkCffccw9//PEHV111FZGRkWzYsIENGzbQo0cP2rdvT7du3di3bx8AH330EcOGDQNgx44dtGrVivT0dKvEJYQQwjE4fBX6a4t2sfvk2cvuk5WTx+mzmfh4uDLj32g2H03Cw6347zYt6vjzyi0tSzz3wYMH+emnn5g2bRodO3bkhx9+4O+//2bhwoW8/fbbzJw5kxUrVuDr68sff/zB888/z7x583jyySfp2bMnv/zyC2+99RZTp07Fx8enzK9dCCGE43L4BF4ap1LO4+6qcHVRBFfz4FTKeRoE+Vb4uA0bNqR169YAtGzZkt69e6OUonXr1kRHR5OSksLo0aM5dOgQSimys7MBcHFx4ZtvvqFNmzY8+uijXHXVVRWORQhnJB1VhSNz+AReUkk59mwGfT5aTas61fFwc+HTu9vTd+JqPh3SnlA/rwqdu/AELS4uLgX3XVxcyMnJ4aWXXqJHjx4sWLCA6OhoevbsWbD/gQMHqFatGidPnqxQDEI4M+moKhyZ07eBT1l1mNuvCCuoMg/19+K29mFWawu/nJSUFOrUqQPAN998c9H2J554gtWrV5OQkMDcuXMrPRYhhBD2xeFL4CXZHpNM1NGkgvvhzy4GoEODgEo/99ixY7n//vsZP348N910U8H2MWPGMGrUKJo0acLXX39Nr1696N69O6GhoZUekxBCCPvg9Al87mPdAOtXr4WHh7Nz586C+4VL2IUf27lzZ0HV+ptvvgnAtGnTCvatV68eBw8etEpMQgghHIdTJ/CiOrjkl8BBOrgIYQbpeCZE6Th1Ai/cwUWIymLNhGSryc2acUnHMyFKx6kTuBBVwZoJyVaTm63GJYQjkwQuRBFstaQrhBD5JIELUQQpUZZfVk4eh+PPEZuaUeG5FIQQxXPuBP7XO7Dq3eIf7/Es9Hqu6uIRwgEUXlfgpZtbmB2OEA7LuRN4r+cuJOjplnHYQxcXv78QNsTWSrq5eZo/95zhTGomfp5u/LD+GNc0DuaqyGDcXZ1+zighrM65E3i+HXPh5GbITodJXaD709B6kNlRiTJytnbripZ0rXG9zmfl8vfBeJbvPs2KPbEknMsCICM7j6zcPB78ZiM+Hq60r1+DDg0C6RgeSPv6NfD1lH89QlSUfIp2zIUVr0NgBHj6GyXyBaONxyqQxKOjo+nXrx9dunTh33//pWPHjgwdOpRXXnmF2NhYvv/+e1q2bMnIkSPZs2cP2dnZvPrqqwwYMIDo6Gjuu+8+zp07B8Bnn31Gt27dWLlyJa+++irBwcHs3LmTK6+8ku+++w6llDWuhN2zh3bripSac3LzSMvMITUjh2OJ6cSmZhIR7Mu8TTE82iOizMcr7/VKSMtkxd5Ylu8+w5oDcWRk5+Hn5UbXiCD+PhhP4xBfvD3ceOnmFtw5dS03t6nNzhNn+eTPA2gNri6KlnX86dAgkHDfHHyDMgnx8yzxvEKIizl+Av/tWTi9o/jHT242knfyMeP+yvHg6QcLH4eo6UU/p1ZruOEybecWJS0n2qJFC3r27MmMGTNITk6mU6dOXHfddYSGhrJ8+XK8vLw4cOAAd999N1FRUQBs2bKFXbt2UadOHa666ir++ecfrr766rJeFWGCnNw8jiemk5qRzfM/7+SmNrVIy8ghNTOHtIwc0iw/L7pvSdhpmdlkZOf955gH487hoqDfxDX0bBJCs9p+NKvlT7PafoRU87Tal7sj8edYvvs0y3efYdPRJPI01KnuxV0d6tGnRS06NQzk3d/2MrhjfXadTAGgVd3qDO5YH6VgyZPXcDYjmy3Hktl4JJGN0Yl8v/4omTl5vLzsOA2DfekYHkCHcKOUHh7kg1LK5poJhLAljp/AS5KdbpS8C/P0N7ZXUEnLicbExLBgwQI+/vhjADIyMjh27Bh16tRh9OjRbN26FVdXV/bv319wzE6dOhEWFgZAu3btiI6OlgRuY7TWxKZmsvd0KvtOn2XvqVT2nk7lwJlUsvM0AH/sOcMfe84UPMfVReHn5UY1T+Pm5+VGcDUPwoN9C+77ebpRzcuNvDzNB8v3U6e6Fxro1iiYWRuOseZAPD9vOVFwzCBfjwsJvZYfzWv70zi0Gl7urv+J+dJEmZen2RqTzPLdZ1i++wwHY9MAaFHbn8evjaRPi5q0rON/0ReEktYV8Pdyp0eTEHo0CSk45/JNe4nJ8GRjdBLLdp/hx6gYAIKredIxPICDsWmkZ+Xw+cpDJa4sWBWcrZlG2DbHT+AllZQndTGqzVeON+4PXQxHVsOSsRXu0FbScqKurq7Mnj27IMnne/XVV6lZsybbtm0jLy8PLy+vIo/p6upKTk5OhWIUJbtcKTA9K4d9p1PZd9pI0ntPn2Xf6VSS0rML9qnp70mzWv64usDRhHTq1PCmSU0/PNxcGNevGX5ebni6uZS6tPz6ot0XlXTfGNgKd1cXlILRvRoXxGF8cTjL9+uPFpTeXRREhFQrSOjNavnRrLY/Jy3t6S/8spMgXw/+2BNLfFombi6KzhGB3Nu5Pte1qElYgE+xcZV1XQEPNxda1vThpvBwHu0BeXmaQ3FpbIxOIio6kbWHE0jNNN7f0/+JZkdMCt0aBdHB0o7u5+VequtlTfbQTCOch+Mn8JJ0f9po8/b0M0reR1Yb93u/XOmnvv7665k8eTKTJ09GKcWWLVto3749KSkphIWF4eLiwowZM8jNza30WETx8juLvff7Xq5tVtNIkKfOsu9MKscS09FGoRofD1ea1PSjX6taNK1pJMZmtfyo4eNRsO58o5BqeLi58OJNzek7cTVj+zUtskR8OZcr6Qb4etC1URBdGwUVPJ6bpzmacK4g7j2nU9kWk8yv209ddFwFLN99Bh8PV3o1C6Vvi5r0bBpKde+qSZQuLorImn5E1vRjSOf6vL5oN4u2naS6txsebq4ciT/H5mNG9b2Lgqa1/OkYHsCVDYyq97o1vKskTiFshSTw/I5qCx83qs2XjDWSdxX0Qn/ppZd4/PHHadOmDXl5eTRs2JBff/2VkSNHcvvttzNz5kz69euHr69vpcdiBlutjtRaE52QztbjSfxzMJ7TZzNRwNxNJ5i76QQuCsKDfWlZx5/brwijaS0/mtfyJyzAGxeXokvR+evO55eaC687X9Ye5GUt6bq6KCJCqhERUo0bW9cu2J6akc3+M6l8sGw/W44l4+/tRtt6NahT3YtX+7cqU0zWXhgo9mwG8zbHEBHsi4ebC5/e3Z6+E1fz51M9OZ6UTlR0ElFHE5m7KYaZa48CRpv8leGBdGgQQIfwAEutx4W/h7SnC0cjCRyMZJ3fYc1K48BLu5zopEmTLqoWB4iMjGT79u0F98ePN6r3e/bsSc+ePQu2f/bZZ1aJ1Sy2Uh2ZnJ7F1uPJbD2ezJZjyWyLSSbZUgXu5qJwd1UE+XoQGepHoK8H4we1sWqp2Sx+Xu7UC/Bh18mzNKtlVOm/OaAVfSeuZmSvxmVKctZeGKi4LzzfrjvKSze34JpIox09JzePvadT2RidSNTRJDYcSWDRtpMAVPN0Kxi+1iE8gBPJMsGMcCzOncCLmont1eoXfpeZ2OxWcaWt7Nw89p5KZcvxJLYeS2bj4TiOp+wCQCloEurH9S1q0a5+DeoH+TDyu80FpcAP7mxL34mrOZuRXeYEXlnrzleUNWsGrKm0X3jcXF1oVbc6repWZ+hVDdFaE5N0nk1Hk9gYncimo0lMXLG/oJnDy92FORuPl2vYnbAdsWczeHrudibc0cap/47OncALz8QmHEp+u/UHS/dxTZMQth5LZsvxZHaeSCEzx+jQFVzNk6bBngzu0pD29WrQOqz6RR2jXl+02yaTmzXZYs0AlP8Lj1KKeoE+1Av0YWD7ugCknM/mmZ+28c/BeLLzNOcycxj53WZ+fLRrsU0eVcFWm5DswZRVh9l2PNmhPovl4dwJXDikNQfiOHM2E1cXmBMVw5yoGDzdjJLavV0a0L5+DdrVq0HdGt4cPXqU8PDwIo9jq8nNmmy1ZsCaMrNzWX8kkSY1/XBRUN3HnT/3xjHkq/V8NqQ9wdXMmUTGVpqQ7E1+/4jvH+7MvV+vd+raFEngwmGkpGfzwfJ9zFx7FAXU8PEgMrQaoX5evH9HWzzcyjYfty0mN2t3FnMGlzYTfP1AR4Z8tZ71hxO48eM1fDy4/UW99oVty/97tqpb3eFqxMpKEriwe3l5mp82HWf87/tIOpeFh6sLzWpVw9vDjY8HG72Xk89nOcS3dGt2FnOWLwOX1qQ0fG4JYExKk5Gdyz1freOJ3pE8fm3kRb3Whe3JL30vH9MdgBE9Iug7cbXTlsKdOoFP3jqZz7d9Xuzjj7V9jJHtRlZhRLbP1trtth1P5uWFu9h2PJmO4QHUjQwm0NfTptqtbTVRWrvnuK26XE3KucwcXpy/k4l/HGDd4QQ+Htyemv7OlwjsRX7pO9TyN7KFz7eZnDqBj2w3siBBD/19KADT+xUz/3klWLhwIdu3b+fFF1+ssnNWlK202yWey+L9pXuZvfE4wdU8+eiutgxsV5c7pqxl/taTBfvZQru1syRKe+Tr6caHd7alW6MgXl6wixs/XsOHd7UrmO5V2Jb82pRp/xy5aLsj9UspC6dO4AD3/3Y/W2K3FNxvPcOY1rR9aHtm3jCzUs/dv39/rr/++ko9h6PJzdP8sOEYE5buIy0zh4euasiT10UW9B63xXZrYduUUtzRoR7t6tVg9A9beGDaBkb0aMRTfZtU2TrmMslM6eR/voXB6RN4y6CWtAxqyd7EvYBRAh+/YXyFV3EqzXKiu3fvZv369Xz++ec8+OCD+Pv7ExUVxenTp3nvvfcYNEjWJC9s09FEXpq/i92nztI1IojXBrSkSU0/s8MSVmZWk0NkTT/mj7qK13/dxZRVh9hwJIFPh1xRJVO0VnRtd+GcHD6Bj98wviA5FyUrN4tdCbvwdPXERblwz+J72JWwi1ZBrQqq1S/VLLAZ4zqNK/HcJS0nOnDgwIv2P3XqFH///Td79+6lf//+ksAtYlMzePe3vfy8+QS1q3vx2ZD23NS6tqyD7qDMbHLw9nDlndva0LVRMM//vIMbP17D+4Pa0LdlrUo5n9aaLceSiU3NJDK0GvM2l29td+GcHD6Bl8TD1YMg7yCSMpLwdPXkdPppgr2DcXet+AIOJS0neqmBAwfi4uJCixYtOHPmzH8edzbZuXnMXHuUicv3k5GTy8iejRjVqzG+npX/trXVjmeiavRvW4c2daszetZmhn+7iaFXhfPsDc3wdCvbDHxFOZ6YzrrDCaw9nMD6w4mcSD4PwKG4c/RqFiKlcFFqDp/AS1NSjkuPo8/cPtT3q8+JcyeYP2A+wd7BFT53ScuJXm5/nT/3o5O4tA1w7aEEXl24i31nUuneJIRXb2lBREi1KotHOp6J8GBf5j3WjXeW7GX6P9FERSfx2ZD2NAgq2+JCJ5LPs+6QkbDXHU4gJslI2IG+HrSrV53Ec1nUreHFyZQMlu+OxdPNhYeubkgdWV1NlMDhE3hphPiEEOQdxL6kfQxpPsQqyVuUTX4b4IfL9nMuK5dF204SFuDNF/ddSZ8WNaW6XJjC082VV/u3pGujIJ75aRs3ffI379zWmlva1im249nplAzWHo5n3SFjTfNjiekA1PBxp3PDQB6+uiFdGgXRJNSPNxfvITyoGrtOphDg40FkLT9+WH+MgZP+4ZdRV8kSqeKyJIFb1PKpRUZOBsNaDTM7FKcTezaDuLQsgn09mL3xOO6uiid7R/JYz0ZlXjREiMpwfctatKzjz+OztvD4rC38eyiBk5bVzT5avp8uEUFGtfihBKITjITt7+VG54ggHugWTteIIJrV8vvP3OuXTjKz0fJ7fFomN32yhg/uaEvv5jWr7oUKuyIJ3MLD1YMmAU2sVvou7XKid999938eB0hLS7NKHLZOa82zP+8ArTmTmkm9QB+6NgwsV/W1tFuLyhQW4MOPj3ZlwrJ9TF11GAAPN8WsDceZteE4fl5udG4YyL1dGtAlIojmtf1LnNmtuGGPR+LPMer7zTw0I4rh3SN45vqmVTakTdgPp07gRc3Elj8OHGQmtsq29XgyL8/fyfYTKXi5udA4yIdvH+pM34mreboc42Gl3VpUNndXF567oTlH4s7x595YfNzdaBtWjTrVvfnwrnZWm4q1YbAvP4/sxluL9/DF6sNsjE7k07vbExbgY5XjC8fg1Am88ExsourEJKXz/tJ9LNh6Ei93F65qHER2Th5KKaefGlHYvtizGaw/kkjrutXxcHPh07uN+fYTzmVadfiXl7srbwxsReeIQJ6dt4ObPvmbD+5oy3UtpEpdGBw2gWutnbLjky33Xk/NyGbyykN8/fcRFDC6V2P+ORTPPwcTCvaxhalPhbicS1c3q+wvnTe3qUOrOtUZ9cNmHp4ZxcNXN2Rsv2ZlXl1POB6HTOBeXl4kJCQQFBTkVElca01CQgJeXrY1CURObh6zNx7no+X7STiXxa3t6/LM9U2pU8Obp2kKyNSnwn6YsU58/pC2t5fs4au/jxB11BjSJlXqzs0hE3hYWBgxMTHExcWZHUqJcnJycHOz3p/By8uLsLAwqx2vIrTWrNwfx9uL93AgNo1O4YFMH9qcNmE1zA5NiHKz1nz7Ze106eXuyusDWtElIohxc7dz48drmHBH20qbJU7YvlJnDqXUNOBmIFZr3cqy7VXgESA/Uz6vtV5ieew54CEgF3hCa73Usr0f8DHgCnyltX7Xsr0hMBsIAjYB92mts8rzotzd3WnYsGF5nlrloqOjCQ8PNzsMq9t7+ixvLd7DmgPxhAf5MOXeK7m+pYznFiJfeTtd3ti6Ni3r+DP6hy0M/3YTw65qyLM3SJW6MyrLX/wboF8R2z/SWrez3PKTdwtgMNDS8pzJSilXpZQrMAm4AWgB3G3ZF2C85ViNgSSM5C/sTGxqBs/OM0oH22NSeOnmFiwb04N+rWpJ8hbCShoE+TL3sa482C2caf8c4Y4p/3LcMmGMcB6lLoFrrVcrpcJLufsAYLbWOhM4opQ6CHSyPHZQa30YQCk1GxiglNoDXAsMsewzA3gVuHiMlyiXoqrqCivv+OjCM1H5ebrz1ZrDfL7qENm5eTzYrSFP9G5MDR+PioQuRIU48twA+bPEdYkI5Jm527nxkzW8P6gtV9SvwdNztzPhjjayKIqDs0bj62il1P1AFPCU1joJqAusK7RPjGUbwPFLtnfGqDZP1lrnFLH/fyilhgPDAerWrVvkwiD2IiEhoeSdKujWSA9ujWwJwJMLjgDw8YCLmxjKcw2PJ6aRlpHD6JnrOJKYSdy5HLo39GN4l5qEVfckOfYkyaU8VkZGRrnjqIiquP6ieJV9/Qu/94tT3vecWe/ZSzWrBlNva8hry48z4rtNNA7y5ExaNu8v2sqobsW3j8t731zWuP4VTeCfA28A2vLzA6DS5yLVWn8BfAHQoUMHbe9tyFUZv5fXKaucM/ZsBknpu3F3dWHD8XM0r+3HZ/e0pHNEkKlxlYe9v3/snb1efzPfs5cKBxa2aMTLC3YyZ2MM10QGs/xgCs/c0u6ypXBbiN3RxJ7NKHUNSEWvf4V6PWitz2itc7XWecCXXKgmPwHUK7RrmGVbcdsTgBpKKbdLtgsb9cmfB0BpcvI03ZsE0zUiqNzJWwhRcZ5urvh6uNO7WSjj+jUrGJsuqtaUVYfZdjy5Sq59hRK4Uqp2obu3AvmTfy8EBiulPC29yyOBDcBGIFIp1VAp5YHR0W2hNmYf+QsYZHn+A8CCisQmKs+ZlPPM2Xic3DyIrFmNCYPa8vOWE8SmZpgdmhBOK/ZsBvM2x/DOba1pVbc6I3pEMG9zjHwuq1D+3+D7hztXybUvdQJXSs0C1gJNlVIxSqmHgPeUUjuUUtuBXsAYAK31LuBHYDfwOzDKUlLPAUYDS4E9wI+WfQHGAf9n6fAWBHxtlVcorO7xWVvJztXUD/TG38v9opmohBDmyJ8hLtTfqLaVz2XVy/8btKpbvUqufVl6od9dxOZik6zW+i3grSK2LwGWFLH9MBeq4IWNWn84gQ3RiQAcSzzPscTzMv2pEDYgf4a4af8cuWi7fC6rRn7pe/mY7gCM6BFB34mrebRHRKWNBnDImdhE5Yg9m8GoH7YQEeLLglFX8fCMKECmPxXCFuTPECfMcbkakMpamEkSuCiV7Nw8Rv2wmXOZOfzwSGf8vNzNDkkIIWyGGTUgksBFqYz/bS8bo5P4eHA7mtT0s8oxHXmSDSGEczGjBkQSuCjR4u2n+OrvIzzYLZwB7YqdX6fMyjsXtBBCiAoOIxOO72BsGmPnbuOK+jV4/sbmZocjhBDCQhK4KNa5zBxGfLcJL3dXJt1zhax2JIQQNkSq0EWRtNaMm7edw3FpfPdQZ2pX9zY7JCFMJ/02hC2RBC6K9M2/0fy6/RRj+zWlW+Ngs8MRwiZIvw1hS6ROVPxHVHQiby3eQ58WNXmsRyOzwxFCCFEEKYHbqMpaw7skcamZjPphM3UDvJlwR1uUUlY/hxBCiIqTBG6jClfV3TV1LVD5M57l5Obx+KzNpJzP5puhnajuLZO1CCGErZIELgpMWLafdYcT+eCOtjSv7W92OEIIIS5D2sAFAEt3nWbKqkPc07k+t18ZZnY4QgghSiAJXHAk/hxP/7iNtmHVefmWypl0XwghhHVJFbqTS8/K4bHvNuHmqph875V4urmaHVL5/fUOrHq3+Md7PAu9nqu6eIQQohJJAndiWmte+GUn+86kMmNoJ+rWsPPJWno9dyFBT7/J+Dl0cfH7CyGEHZME7sS+W3+MX7ac4P/6NKF7k5BSPcdpZqKS0rwQwsZJAndSW48n8/qiXfRqGsLoXo1L/TynmYnKmqV5+TIghKgEksCdTFZOHofi0nh0ZhQ1/b346K52uLjIZC2VSqr2hRCVQBK4kzmZcp6zGTmkZeawcPTV1PDxMDskIYQQ5SDDyJxI7NkMYlMzAfBydyXU39PkiIQQQpSXJHAn8vqvu9Eagqt5MLhjfaauOmx2SEIIIcpJEriTOByXxuIdp/BwVTQI8mVEjwjmbY4hNjXD7NCEEEKUgyRwJzH82yi0hsah1XBzUYT6e3Fb+zAphVelHXPh5GY4+jdM6mLcF0KIcpJObE5g0baTHIw9B8DuU6nAhbHbHRoEmBaXU9kxF1a8DoER4Olv9EpfMNp4rPUgc2MTQtglSeAO7mTyeV74ZQft69fgp0e7cs9X64HKX5pUXGLVeLj6f7BuKmSfh4bdYcBnsGRs2RK4jCkXQlhIAndguXmaMXO2kpunmXhXO9xcpcWkzPKrvbPTjWrv7k8Xn3CzzkHiYUg4BImHIOGw5echOBcLv465sO/kbsZx4veVLR4ZUy6EsJAE7sC+WH2Y9UcSeX9QGxoE+Zodjv0pstp7FJw9AYGNLiTn/KSdevLi51eraezXpC/sXwodH4Z9v0FmKnj4wIrXjP1m9Ie2g6H5LeDpV/WvUwhhlySBO6idJ1L4cPk+bmxdi0Gyvnf5rH4fOj0CaydDSgzMHwUpx2D5yxf28QkyknRED+NnkOUWGHFxMs7/MuDpZyT2Xs/BLyOg7hVwegfMfwx+/T9ofjO0GQwRPcFVPp5CiOLJfwgHdD4rlydmbyHI15O3b22NUjJVaqnlZMLhVbBnAcTthWUvgnIBdx+o3xlq3AVrJsDDKyAoArxL2Qkwv9p94eNGdfySsdDndWO71nB8A2yfDTt/hh0/gW8otL4D2twJtduC/A2FEJeQBO6A3ly8myPx5/j+oc4yVWppZKXDwT9gz0KjqjvzLHj4GdXmXUfCoVXg4gq3fwVHVsPeXyHsyrKfp/UgiJpu/F643Vop48tB/c7Q7104sAy2zYYNX8C6SRDSDNrcZSTz6pbalLK0zQshHJIkcAfzx+4zfL/+GMO7R9CtcbDZ4diujBTYv4yQTbPgxL+Qcx68A6FFf2g+wKgS37PoQrW3p7+RvBeMht4vl3z88nLzNNrCm98C6Ymwez5sm2O0l694HcKvhqDGcHCFDEkTwslJAncgsakZjJu3nRa1/XmqrxMs+Vmc4kqn5xJg3xKjpH14JeRm4ekdAu3vgeb9ocFVF7c7F1Xt3fvlqkuSPoHQYZhxSzwC2380qtmj14CLuxETqvxD0oQQdk0SuIPQWvPMT9tJy8xhzt3t8HRzNTskc1zac7zTcFj8f0aHtPgDoHOhRn1je/P+xOSGEN4wovjjFVftXdUCG0LPcdBjLLweCO3vhy0zID0eZt8D175U9iFpQgi7JgncQcxce5RV++N4fUBLGofayVCkypiUZPUEuGUi/PoUJB+FuUMBDVlpxkQqzftf3CksOrp8sZtFKQhuCq0GQuxuY+ja4ZXG8DSvGnA+qfQd64QQdk0SuAPYfyaVt5bsoVfTEO7r0sDscEqvMiYlid8H67+ApMPg7gu9XoAmN8AX11Ru23VV6v600ebt6Qf+daHPG7BgJJxPhE/aG198Oj4Eru5mRyqEqESSwO1cZk4uT8zagp+nG+8NauvcQ8ayzoGbF+z/zahC96sDPZ4xOp8FNzUvrqJqGl6tfuH3stY0XNo2fz4JbvkYQprC0hfg93Gw8Uvo+yY06SdD0IRwUJLA7dyEpfvYezqVaQ92IMTP0+xwzHM+GX6405hn3CfIGLet86qm53hJCtc0WEtxbfP3LzCGwi17EWYNhoY94Pq3oFbr/x5D5lUXwq5JArdjfx+I58s1R7ivSwOubVbT7HDMcy4evr0VYvfAHd9AXo55PcfNphQ07QeNexsJfuXbMOUaaH+v0dHNr9D7ROZVF8KuSQK3U0nnsnjqp600CvHl+Rubmx2Oec6egpkDjA5rd8+CyD7GdlvoOW4mV3foPBza3GF07Fs/1Zjl7Zox0HU0uHubHaEQooIkgdshrTXP/7KDxHNZfP1AR7w9nHTIWFK0kbzPxcO984xJTsTFvAOMKvQOw+CPV+DPNyHqG7juVWh1O7hYaYU6qY4XospJArdDP0XF8NvO0zx3QzNa1a1e8hMcUdx+I3lnp8P9C8s3takzCWoEd30H0X/D0ufh54dh/edw/dvGQi0VnZZVquOFqHKSwO1MdPw5Xl20i64RQTxyzWUmIHFkp7Ybbd7KBYYugZotzY6o6lS0R3v41fDISthumZ512vVGhz//MPANlmlZhbAjksDtSHZuHk/O2Yq7qwsf3NkWFxcnHB50fAN8P8hYbOT+BRDcuHLOY+2hX9ZijR7tLi7Q7m5j3veJbYz1yRMOQPY5CG0p07IKYSckgduRT1ccYNvxZCYNuYI6NZywE9LhVTDrbqgWCg8sNKZErSyVMfTL1nj4GpO/PLEVvroOzp6Aj9tAp0dlWlYh7IAkcDuQlZPH/thUNhxJZNCVYdzUprbZIVW9fb/Dj/cbE7TcPx/8apkdkWMIbmr04A9uYlSjB9SHvz8wmidWT4DOI8CzmtlRCiGKYKUuqKIynUg+z7nMXKp5uvFqfydq7823cx7MuQdqtjDavCV5W0/+tKwZyeDuZUzB6lcbaraGP9+Aj9vCv58ZE+QIIWyKJHAbF3s2g7jUTAA0kJ6VY25AVW3zTJj7EIR1Mnqb+wSaHZFjaT3ImOgm8TAc+9do++77JoxYDQ+vMGZwW/YCfNwONnwJOZlmRyyEsJAEbuM++fMAGgjx8+TODvWYuuqw2SFVnXWfGzOqNepljPP28jc7IsfUehDUuQIaXA2j1l3ovBbWwWiueHCx0XSx5Gn4tANs/hZyneyLpBA2SNrAbVjs2QzmbToBQKifJyN6RNB34moe7RFBqJ9XqY/z0fL9fLziwEXbwp+9MEb3yd6RjOnTxDpBW4PWRvvrX29Cs5th0DRwc+J53s0WfrXRdHHoT/jrLVg4Gv7+EHo+D61uAxcnnUhICJOVOoErpaYBNwOxWutWlm2BwBwgHIgG7tRaJyljSayPgRuBdOBBrfVmy3MeAF60HPZNrfUMy/YrgW8Ab2AJ8KTWWlfw9dm1KasOU83T+BP5ergS6u/Fbe3DmLrqMC/d3KLUxxnTp4ltJehL7Zh78UQiNerDgaXQZjAMmASu8j3TdEoZ86s3uhb2/w5/vmVMBrPmA+j1PDS/RVY9E6KKlaUK/Rug3yXbngVWaK0jgRWW+wA3AJGW23DgcyhI+K8AnYFOwCtKqQDLcz4HHin0vEvP5XTWH0kgLi2T89m5bIhOIvzZxUz75wjbjiebHZr17JgLK143qmjrdYWAcCN5N+wBAz+X5G1rlIKmN8Cjq42FY3Qu/HgfTO1urIK2/Sfjy9jRv40vYzvmmh2xEA6r1P8dtdarlVLhl2weAPS0/D4DWAmMs2yfaSlBr1NK1VBK1bbsu1xrnQiglFoO9FNKrQT8tdbrLNtnAgOB38rzohxF14gg9p9JpXXd6ri7ujDn0a5mh2R9qycYE4f89S4k7IdzcdDyVojda715ukXRKjJZjYuL8Xdq3h92/AQr3zGWc3X1NHqx+9WGa5+XWd2EqEQVLd7U1Fqfsvx+Gshfq7AucLzQfjGWbZfbHlPE9iIppYZjlOypW7cu0dHR5X8FJktISChye3ZuHj9FHaNrAz9SzmeRm41dv87iNIjfx1Fdm7CEw7idiyOp/ShSWj1Ig+86c7SMr/eH6B+YdXSWcSc/98+4sA723Q3uZkj4kIueU9z1dwoN7zZul1Oav0H1LnDLXMJ+ugGVl4NrcjR55+KJj95DXsfnCFzxNif9OhT5VKe+/iaTa28ua1x/q9VPaq21UqpK2qy11l8AXwB06NBBh4eHV8VpK01R8f++8xQpGbkM7dGUKSsPFbuf3QtuSnjiGjh3CnxDCBjwNgFHVhvby/h6nw9/nud5HoCh3xgJY/qDUSU+zyGvqxkyEmDcMZjaHZfUE4SuGgs+IZAeT3h1BQENinyaXH/zyLU3V0Wvf0XrKM9Yqsax/Iy1bD8B1Cu0X5hl2+W2hxWx3WnN2XicWv5edI8MMTuUynXVE8bYYxc3qNEQjqw2ql27P212ZKKsgpvCqS3gXwfqdIB75kFQBKCNCWG+v9OYUS8v1+xIhXAIFU3gC4EHLL8/ACwotP1+ZegCpFiq2pcCfZVSAZbOa32BpZbHziqlulh6sN9f6FhO53RKBqv2xzHoyjBcHX3Bktg9kJcNyhVObDCSee+Xpc3UHhWe1Q0Nbh6QehpueA+6PwOntsGsu4xkvvp9SD1jdsRC2LWyDCObhdEJLVgpFYPRm/xd4Eel1EPAUeBOy+5LMIaQHcQYRjYUQGudqJR6A9ho2e/1/A5twEguDCP7DSfuwDZvcwx5Gu7oEFbyzlWtqI5PhZVlla6ja+HfT+HKoRBvGadekTWkC8U2PX+bLawg5izyv3QtfNwYEnjpl7EeY2HfbxD1Nfz5Jqx8l5B6vaDH49CwuwxDE6KMytILvbjeLr2L2FcDo4o5zjRgWhHbo4BWpY3HUWmt+THqOF0iAmkQ5Gt2OP9VeJWu6TcZP8uTdDPTYP4IY8x33zfgh8FWjW3o70ONEPtNv9wzhLW1HgRRlmt+6fvC1d1YwrRFf0g4BFHT8Nr8LczsD0GNjS9y7YbIdLlClJIMsrUx648kcjQhnSd7R5odSuX64xVIOgoP/gqefmZHU6kmb53M59s+L/bxx9o+xsh2I6swIhsQ1Aiuf4uYRvfT4NxWo1S+7AVjToDgSDizs/jnSk2KEIAkcJvzY9Rx/DzduKGVAy8ZeuhP2PgVdBllTNPp4Ea2G1mQoKVm4GLazRPa3mXcTu+EqGmwfY7xoKu70Yu9ehj0fsloX5f+EUIUkJkybMjZjGyW7DjFLe3q4O3hoPNLZ6QY/4iDmxj/lIXIV6sV3PwhPLXXWDLWPwxST0LMRljyjNG7/bexsHUWxGwy3ktCODEpgduQX7edIiM7j7s61Ct5Z3v127NGz+SHloO7t9nRCFvk6QdpsfDCGfjqOkhPgMBGELfP+H3+iAv7VqtpfBkMamz8DG4CwY2her2LF1nZMdeY9S9+nzHcrfvTUpIXdk8SuA2ZE3WcpjX9aBNWveSd7dHexbDtB2NIUdiVZkdzWdJubbLgpnB8nZHMPf3g7h+MOQKWPAN3fQfx+42RC/EHjN93/WIZvmbh5mUk9aDGkJcDx9bBtS9Ai1vhzA6Z4lU4BEngVlTUsp2FXW7Zzn2nU9l2PJmXbm6BcsThNOcSYNGTULM1dB9rdjQlknZrE/31DsTtgRm3XNhWeDjgjrn/7cSmtVE6L0js+yHhIBxZBeeTjH1+HWPc8s17CM7sguteqbzXIkQlkgRuRYWX7bxr6lqAUi9A8mPUcdxdFbe2L3YKePulNSweA+eT4b75xgQfQhQnfzhgWaq9lQLfYOPWoNuF7SXNW/DPRDi21hiH3rA7hHWUteeF3ZAEbgOycvL4ZcsJ+rSoSaCvAya3nfNg9wLo/YrRUUmI0mg9qOJV3L2eM957N75nJOh8+5cbXypbD4LDq4yZ4VaNBzdvqN/FktB7QO22Fy9pa82JjISoIEngNmDFnjMknsviDkfsvJZ6GhY/ZZRsuj1hdjTCGeVP8TrgM6jf1ShxL3kKrnv1wheE88lw9F+jnf3IaljxmrHdszqEX3WhhN5j3IXagfwZ50KaS6c4YQpJ4DbgxygHXbhEa+OfXE4mDJxycUlGOJaKrC1e2fIT65KxF6rjLx1P7l0Dmt1o3MDoBR+95kJC37fE2O4TDIENIfGwMczNJ0jWPRemkf+oJotLy2bV/jhG9mzseAuXbPkWDiyDfuONoT3CcRWeYtcWlbU6vlootLrduAEkH7+QzHfONXq2pydAkjus+9yobl/xGjTpB57VKuc1CHEJSeAmW7o/2XYXLqmIpKPw+3MQfg10Gm52NEJUTI160P4e47bjR3hsHcy5BzJTjV7vCZYS+rv1ILQlhHWAep2MpqPARuAic2YJ65MEbqK8PM2Svcm2u3BJeeXlwYJRgIIBk+Sfl3AswU0hPQ78ahu3oYthz2L4fZyxGEvMRtj5M2yyDDv0qmEk9LCOxq3ulUaVfT6ZZEaUkyRwE22ITuTk2SyeucHBOq9t/NJoP+z/KQQ0MDsaIawrv1Ocpx94+hvV6kufu7hTXF6eMRY9ZqPlFgUr3wW08XhwUyOZKwUHV8DAz43OcsfWSnu6KDVJ4Cb6ceNxfD1c6NfSgRYuiT8Iy1+ByL7Q/r7L72vLHZ+EKE5J656DUesU2sy4XWH5HGSchZObLyT0XfMg+7zx2LcDLj7HvIeMqnl5/4vLkARukrMZ2SzZeYq+kdUdZ+GS3Bxjnmo3T7jlE6N0cTm23vFJiEsV9aUzbo+RcOc9dPkvnV7+ENHTuIExSuP1QOj/Gax4AzJTjC8E+Y7+Y5TaG1xlVMHL2gHiEpLATbJo20kysvO4sVmA2aFYz7+fGKWL278GfweqVRAinzW/dK58F3QeLChmTv34/RD9N6DB1cNoO29wlVHVXq8zeDhQvxlRLpLATfJjVAzNavnRNMTL7FCs48wu+OttaDHgwtAbIUTxej0HwZGw4vWLJ5kpvO75+WRjIZajf0P0P/D3R7BmAri4QZ32loR+tZHQvfylQ5yTkQRuAodbuCQnC3551OhZe9NHJVedCyEMJU0y410DmvYzbmAMWzu23qheP/oPrJ1kzOeuXIwlVM8nGTMeXvkAxO2VDnEOThK4CQovXHI27qTZ4VTc6vfg9A4YPAt8g8yORgj7UpZJZjz9IPI64waQlQ4xG4zS+drPIDcb/nrTaJcfNM0o2S8Z61gJXGoZCkgCr2KXLlxyNs7siMppx1yjR212ulG9V7/bhWkohRBVw8PnQse4NRNg3DE4vd3YDka1fPw+MyO0rh1zi25yAKdM4jLDRhXLX7jkTnteuCT/Q1QjHFy9jPmhU44b24UQ5ghuCqe2GJ3c6rQ3th1ba2x3FKsnGMm7YXdwdTd+DvjM2O6EpARexeZEHad2dS+uqcqFS6y9BOLqCXDj+0a7d24G3PaD8WFytKo6IexJUauu5XeIcxTx+4zXVpij1TKUgSTwKnQq5Tyr98cxqlcVL1xSeOjL9JuMn0MXl/948fuMHufnk4x5nhv3NtrenPRDJIRNKM2qa/YuuKnxxaTw2u6OVstQBpLAq9C8TTHGwiVX2nH1efwBUK4Qu8dYB9nH0mnNiT9EQtiMsq66Zm+coZahDCSBV5G8PM2PUTF0jQiifpCP2eGUz7F1MGswuHkZSya6uBoTURxZ7dQfIiFEFXGGWoYykAReRdYfSeRYYjpj+kSaHUr57F4A8x6B6mFw71w4sfnyc0ELIURlcPRahjKQBF5Ffow6jp+XGze0Kt0Uox8t38/HKw5ctC382Qvt1k/2jmRMnyZWjbFYayfD0ueN1ZPunm2M9Q6MgCjLcokVaU8XQghRLpLAq8DZjGyW7DjFoCvD8HIv3cIlY/o0qboEXZy8PFj2IqybBM1uhtu/kgUVhBDCRkgCrwILt54kMyePuzraUee17AxjmNju+dB5BFz/ttHmLYQQwiZIAq8CP0Udp1ktP1rXrV7yzrYgPRFmDzF6ePZ9C7qOkvnNhRDCxkgCr2R7T59lW0wKL9vLwiVJ0fDdIEg+CoOmQ6vbzI5IFDJ562Q+3/Z5sY8/1vYxRrYrZnlKIYRDkQReyX7cGIO7q2Jg+7pmh1KyE5vhh7sgNwvuXwANupkdkbjEyHYjCxL00N+HAjC933QzQxJCmETmQq9EeVrzy5aYgoVLbNr+ZfDNzcYY74eWSfIWQggbJyXwSpSUnk1SerbtL1yy6Rv49f+gVisY8iP41TI7IlEFpDpeCPsmCbwSxaVmVv3CJWWhNfz1Fqx+HxpfB3fMMGZYE05BquOFsG+SwCtJWmY2Keezuf2KulW7cMnlFF7De1Jn8A4wepq3vw9u/shYUUwIIYRdkAReSY4mpAOQlpFjciQW+Wt4B0aAuw/kZhrJu8VA6P+pDBMTQgg7Iwm8Emw+mkhaZi6+Hq4s23OGp1MzCPXzMjeo1ROMFXz+eB3O7DR6mnd7Eg4sk+QtbIq0zQtROpLArSwuNZMHpm9EAY1CqtEhPJCpqw7z0s0tzA0sfh/U7QBxeyAnE+75yVhTd+2n5sYlxCWkbV6I0pEEbkVpmTnc+/U6UjNyiAythreHKyN6RNB34moe7RFhbik8uCn8+j/ISoPgZtC4t7EMqIlreEtJSwghyk8SuJVk5eTx2Heb2H86jT4tQjl73mj7DvX34rb2YeaXwlv0h1XjwasG+ATaxBreUtISQojykwRuBXl5mmfmbmPNgXjCg3xYvju24LH8JUA7NAgwKzzIOAvbZoFPMGSdg2P/yhrepXD/b/ezJXZLwf3WM1oD0D60PTNvmGlWWEIIAUgCt4q3l+xhwdaTPHN9U0b1agzAXVPXAjDn0a5mhmb4bSykxMDQ342e6CBreJdCy6CWtAxqyY74HbgqV2bcMIPxG8bbx5z2QgiHJwm8gr5cfZiv/j7Cg93CGdmzkdnh/NfOn43Sd49xUL+z2dHYlbub3c2tC24lKy8LDxcPZu6ayaJDi5g/cL7ZoQkhhMyFXhG/bInhrSV7uKl1bV6yxdXGUk4YHdfqXgndnzE7Gruy7tQ6hi8fTlZeFl6uXri5uPF+1Pvkkce6U+vIzcs1O0QhhJOTBF5Oq/bH8cxP2+kaEcSHd7W1ndnW8uXlwfwRkJsDt30ps6yVUmpWKq/++yqPLHsEdxd3Pun1Cdl52TTwa4CPmw8h3iE8t+Y5bl94O8uil5Gn88wOWQjhpKQKvRy2xyTz2HebiKzpx9T7r8TTzdXskP5r7WdGT/P+n0KQFav2/3oHVr178bZXq1/4vcez0Os5652vCq2OWc3ra18n7nwcQ1sNZWTbkXi5eRHkHcT+5P0MaT6Epzs8zfKjy5m8dTJPrXqKpgFNGdVuFD3r9bS9GhghhENz+gT+0fL9fLziQLGPP9k7kjF9mhTcPxJ/jqHTNxLo68GMoR3x97LBku2p7UZntWY3G/OcW1Ov5+w2QRcnJTOF9za+x8JDC2lcozETe02kVXCrgsdr+dQiIyeDYa2G4aJcuD78eq6rfx1LjixhyrYpPPHXE7QKasXo9qPpVqebJHIhRJWwSgJXSkUDqUAukKO17qCUCgTmAOFANHCn1jpJGf/dPgZuBNKBB7XWmy3HeQB40XLYN7XWM6wR3+WM6dOkIEGX1HM8NjWD+6etRwMzh3Ui1N/k6VGLkn0efn4EfILglk9kmtQSrI1fy9T1U0nJTGFE2xE80voRPFwvXrvdw9WDJgFNCPYOLtjm6uLKLY1u4YaGN7Do0CKmbJvCiD9G0D60PY+3f5yOtTpW9UsRQjgZa5bAe2mt4wvdfxZYobV+Vyn1rOX+OOAGINJy6wx8DnS2JPxXgA6ABjYppRZqrZOsGGO5pWZkM3T6RuJTs5g1vAsRIZW87GZRVdWFFVdVvfwViNsL9/4MvkGVF5+dS8xI5O31b7M0einNA5szpc8UmgU2K/Nx3FzcuDXyVm6OuJmfD/zMF9u/YNjSYXSu1ZnR7UfTLrSd9YMXQggqtwp9ANDT8vsMYCVGAh8AzNRaa2CdUqqGUqq2Zd/lWutEAKXUcqAfMKsSYyyVzJxcRny3iX2nU/nqgQ60q1ej8k9auKp6+k3Gz5LGbh/4AzZMhc6PGVOliv/QWvN79O+8s/4d0rLTuC/8PsZcMwZ3l4o1hbi7unNXs7sY0HgAP+3/ia92fMV9v93H1XWvZnS70YzfOF4mhRFCWJW1ErgGlimlNDBVa/0FUFNrfcry+GmgpuX3usDxQs+NsWwrbrup8vI0T/24jX8OJvDBHW3p2TTU7JCKdi4e5j8GIc3hulfMjsYmxaXH8ca6N/jr+F+0Dm7N691exy3FrcLJuzAvNy/ua3Eft0fezux9s5m+czqDFw8mrFoYNzW8iTPpZwBjyliZFEYIURHWSuBXa61PKKVCgeVKqb2FH9Raa0tytwql1HBgOEDdunWJjo62ynEzMjIACo6nteazf0/z645EHu1ckyuDckp9rkuPVZyEhIQSj1XLcqzTxR1La0L/HIP3+SRO9p5E9okz5T+WSUp7vcpzPK01f575ky8PfUl2XjZDI4YyIGwArimupbr+5Y3t2mrX0qVDFxadWMS8Y/NYfGQx1dyqEewZzLo961hwYAGTOk4q92u25jWz9vUvrZKuv1lxOYPSvPdF5bHG9bdKAtdan7D8jFVK/QJ0As4opWprrU9ZqsjzJwg/AdQr9PQwy7YTXKhyz9++spjzfQF8AdChQwcdHh5ujZeBl5dRYZB/vCmrDjFvRyLDrmrIszc3L1Np6dJjXU6J+3h5XX6/Td/A8b+g71vUvaJvxY5lEq+91onr0vnLb1l1S8HvV4RewWvdXiO8+sXnKOmcFY2tReMWjMgcwYg/RrAzfidpOWk8uuFRgn2C+Tf9Xzr6daRNSBu83MrWKdJa18wax6rIynKXO6c1X6P4L7mu5qro9a9wAldK+QIuWutUy+99gdeBhcADwLuWnwssT1kIjFZKzcboxJZiSfJLgbeVUvmrfvQFTBuvNG9TDO/+tpdb2tbhxZvKlryrVPxB+P05aNgDusjSm/nzl+9J2EP8+XjiM+LJyMngitAr+Or6r3BR5sxdVN2zOp/0+oQ+c/tQx7cOsedjCfAM4IvtXzBFT8HdxZ02IW3oWKsjHWuWL6GbSVaWE6LqWaMEXhP4xZLg3IAftNa/K6U2Aj8qpR4CjgJ3WvZfgjGE7CDGMLKhAFrrRKXUG8BGy36v53doq2p/7Ytl7LztXNU4iAl3tMHF1mZZy5ebbQwZc/WAgZ+Di0ysN6zVMPrP749Gcy77HO1D2nMw5SDv9XjPtOSdL8QnhCDvIE6eO8mQ5kMY23EsqVmpbIndQtTpKDae3nhRQm8d3NpI6LU60jakbUFCl1XShBBghQSutT4MtC1iewLwn67Qlt7no4o51jRgWkVjKq+snDwOxKby2LebaFbLjyn32ugsa/lWjYeTm+GOGVD9Mv39HHj2tMJSMlP4asdXnMs+B0ADvwa0DG5Jq5BWF43hNlPhSWEA/Dz86B7Wne5h3QH+k9C/3PElU7dPvSihB3kFcXezuzmQZExAJB3ihHBOTj8TW2HHk9JJy8zFz9ONb4Z2ws8WZ1nLd3QtrPkA2t0DLQdefl8HnD2tsOy8bH7a9xOTt00mNSuVWxrdwuLDi/F192XR4UXMHzDf7BALFDUpTGGlSej586/7uPkQ6hNK7LlYm3udQojKJwncYv/pVOLTsnBVxuRlGqt1mre+jBT4ZTjUqA83jDc7GlP9c+If3tv4HodTDtO5dmfGdhxLk4AmrDu1jn1J+xjSfIjNlL7Lo7iEPmnLJPYl7SP6bDTDlw/npoY32fXrFEKUnSRwi1kbj+Hr4Up4sC+dGwYxddVhXrq5hdlhFW3JWEiJgaG/g6dflZ66Ir2NrelIyhEmRE1gdcxq6vvV55Nen1y0oMilVdWOIj+hNw9sTp+5fQj2CuZQyiFclAsn0k5Qt5rpUycIIaqIJHAg9mwGP28+QZOafni4uTCiRwR9J67m0R4RhPrZWE/gnfNg+2zoMQ7qd67y05vd2zglM4Up26Ywe+9svNy8eOrKpxjSfEip5i93JPkd4hLOJ9CrXi+iTkdx1693Mf6a8VxV9yqzwxNCVAHptgxMWXWY268Iw8PNuByh/l7c1j6MqasOmxzZJVJi4NcxULcDdH/G7GiqVE5eDnP2zuHmX27m+z3fMzByIL/e+isPtnrwP8nbWdTyqYWvuy8vd32Z2TfPJtQnlMf+eIwp26bIOuVCOAEpgWOs7x119MKaKeHPGnOOd2gQUNxTqp7W8MsIyM2B274AVxvuYGdla0+u5b2N73Ew+SAda3VkXMdxNA1sanZYpruolsEbvr/xe15b+xqTtk5iZ/xO3r7mbfw9/M0OUwhRSSSBA3Mf6waUvJyoKXbMNYaKZacb9694AIIamRtTFTl69igToiaw8vhKwqqFMbHnRK6tf60MlyqGt5s371z9Dm1D2vLehvcY/OtgPur5kXzZEcJBSQK3ZTvmworXwa82JB6Cel3g8Epje+tBpT6MrXQ8K0pxk5KE+oSSmJGIh4sH/7vif9zb4l48XT1NidGeKKW4u9ndNA9szlMrn+LeJffycteXuaXRLSU/WQhhVySB27LVE6DXC7DwcWO2tcE/QOwuoxd6GRK42R3PLid/6tO9iXvRWnNTo5sYv2E8semx3BZ5G4+3f9xhO6JVpnah7ZhzyxyeWfUMz//9PNvjtjO241jcnajpRQhHJwnclsXvg+Uvg86Fmq3ANwjqdzW2O4hhrYYxcMFA6lSrw9GzR9kUuwlX5crU66bSrW43s8Oza8HewXzZ90smbprIjN0z2JO4hw96fEBN35olP1kIYfOkF7qtit1rzCiTcx5qtgaPasb2Y2sh2HHaNEN8Qrix4Y0cSDpAVm4W19S9hiHNhkjythI3Fzee7vg0E3pMYH/Sfu789U42nt5Y8hOFEGW3Yy5M6gKvBRg/d8yt1NNJArdFZ3bBNzeBuy94+EJeNug8OLIaFoyG7k+bHaHVaK2JTY8lV+dSx7cO2+O3M6y1Y02+YguuD7+e2TfNxt/Dn0eWPcKMXTMwliUQQlhFfp+lG9+DF2ONnyter9QkLgnc1pzaBt/cbLR5P/IX9HkDEg/DsX+Ntu/eL5ep/dvW/bD3B/48/ie+7r6cPHeS/o36S5t3JYmoEcGsm2Zxbf1rmRA1gadWPVWw8IsQooJWT4ABn0HD7sYw34bdjfurJ1TaKaUN3Jac2ATf3gqe/vDAQgiMgODGEGXpcDZ0sbnxWVnU6SgmbJzAtfWuJT49nujUaIeb+tTWVPOoxgc9PmDGrhl8tPkjVsesJjM3s+BxWZpUiHKK32f0USqskvssOX0C/2j5fj5eceCibfkTuQA82TuSMX2aVP6xjm+A724H7wB4YBEENCjlK7BPZ86d4elVTxPmF8ZbV7/F438+7tBTn9oSpRQPtnqQFkEtGLliJG7KjQb+DQjwCpClSYUor+CmRh+lht0vbKvkPktOn8DH9GlS6gRd4rHc5jHG693id3B7FihiWc+ja+H7QVAt1Eje1cOsEo+tysrN4v9W/R/nc87z9fVfUy2/g56oUp1qd2JGvxncs+QeDqUcIjQrlCMpRxxqaVJbngNBOJjuTxt9lAZ8ZpS8j6017vd+udJO6fQJ3KoKr7s9/SbjZwnV3l6nNsKfT4B/XaPa3L9OJQdpvvc2vsf2uO182PNDGtVwjlnlbFXL4Jbc1fQufjnwC7HnY7lj0R20D21PgKe50wgXN8FPWav2bXkOBOFg8vsmLRlrVJsHN630PkuSwM106C9C/xgNgeFw/0Lwc/zxub8c+IU5++YwtNVQ+jToY3Y4Ani49cPM2TeHhv4NOZZ6jHWn1nHXr3cxrtM4OtbqaEpMhSf4AWymal9K9OKyWg+q0k7GksDNcuAPmD2EHP8GeDy4GHwdv+13V8Iu3lz3Jp1rd+aJ9k+YHY6wyF+a9FjqMYY0G0Kb0DZ8FPURw5YOo3f93jzV4Snq+dWr0pgebPkg/ef3x0W5kJ6Tzr1L7mVv4l4m9JhAdl427i7mzCgnJXphSySBm2Hfb/Dj/RDSjNM9P6W+EyTvpIwkxvw1hiDvIN7r/h5uLvLWsyW1fGqRkZPBsNbDCPYOpmdYT2bunslXO75i9fzV3NviXoa3Hl7p/RXSs9NZcGgBP+z5gfScdFxwIdA7kONnj5OZm8njfz6Ot5s3V9a8ks61OtOpdieaBjTF1cW1UuMSwhbJf9GqtnshzB0KtdrAfT+TdybF7IgqXU5eDs+sfoaE8wnMvGEmgV6BZockLnHR0qSAl5sXw9sMZ2DjgXy8+WOm75zOwoMLeeKKJxjQaIDVE2ZMagyz9s7ilwO/kJqdSuvg1rzQ+QXe3fAuod6hnDh3gvn953Pk7BHWn1rPhtMb+GDTBwD4e/jTsVZHOtXqROfanYmoHmF6VbsQVUESeFXaOQ/mPQJ1r4R754JXdcDxE/inWz5l/an1vN7tdVoGtzQ7HFEGoT6hvHX1WwxpNoTxG8fzyr+vMGvvLMZ2HFvh9nGtNTuSd/DRnx+xMmYlLrjQp0Ef7mlxD21D2gLw5Y4v2Ze0jyHNh9AooBGNAhpxXYPrAIhLj2PD6Q1sOL2B9afWs+LYCgCCvILoVLtTQQn9hTUvsCWu4h3ihLA1ksCryrY5MH+EsSToPT+Cp5/ZEVWJ5UeXM23nNO5scie3Rt5qdjiinFoGt2RGvxksjV7Kh5s+ZNjSYfRp0If/u/L/CPMr27DHzNxMlhxewvd7vmdf0j5qeNbgoVYPcVfTu/6z0EpB1X4RE/yE+IRwU8RN3BRhjPiISY1h4+mNrD+9nvWn1vPbkd8A8HXzpVH1RuTpPPw8/Pj+pu9tokOcEBUlCbwqbPkeFoyC8KthyBxjfnMncCj5EC/+/SJtQtowrtM4s8MRFaSUol/DfvSs15MZu2bw9c6vWXl8Jfe3uJ9H2jyCr/vl39ex6bHM2TeHufvnkpiRSOMajRndZDQPdHwALzevIp9zadX+5YT5hRHmF8atkbeiteZIyhHWn17Pmpg1rDmxpmC/b3d/61Bj3YXzkgReGXbMhZObITsdPmgKqachopexnreHj9nRVYnUrFT+99f/8HLz4sMeH+Lh6mF2SMJKvNy8eLTtowxsPJBPtnzC1zu/Zv7B+fi6+3Is9VjBfvlV1U0CmtC4RmOWRS8jV+fSo14P7m1+L51qdeLo0aPFJu+KUEoRUSOCiBoR3N3sbsZvGM+Sw0vI0Tm8t/E9utXpJrP+CbsnCdza8lekCYyArHOQdATcvKHNnWVK3pN/uZvPz+407uQvOWP5hwjwmH8rRt46y4qBW0+ezuOFv1/geOpxvur7law/7aBq+tbkravfYnDTwYzfOJ5tcdsI9AwkxCcEX3dfBjcbzHsb32N/0n5OpJ1gcLPBDGk2hHr+VTskDYx152ftnUVkjUjSs9P59+S/LDq0iFsa3VLlsQhhLZLArW31BOj/KSx8ApKjoemN0PFhWPoCtBtS6sOMvHUW+dNB2Nt40693fM1fx/9iXMdxdKjVwexwRCVrHdKab2/4ljn75vD2+rdJzEzETbmxOXYzLrgwut1o7ml+j6lT5uaPdT+QfIA7m97J4ZTDvPD3CwCSxIXdkgRubfH7IGqakbx9guGOGaBUpa5IY0v+OfEPn275lBsb3sg9ze8xO5xyK2rGrdaFa0Bkxq2LKKUY3Gwwh1IOsejgIjxcPWgX2o661eryaNtHzQ4PuNAhbnib4fi6+/L4isd58Z8XUUpxc8TNZocnRJlJAremM7vBxQ12L4AaDcA/DNw84MjqSl2RxlbEpMYwdvVYIgMieaXrK3bdy7fwjFui9Ia3Hs5P+36igX8DNsdu5uWulbeQQ1ld2iHu096fMnrFaF74+wUUqqA3uxD2QhK4tWybDYv+B25expKgnn6ANpJ3Ja9IYwvO55xnzMoxaDQTe07Ex905OutVNVuvGcivqs4fu23LHcW83bz59NpPGf3naJ7/+3lclAs3NLzB7LCEKDVJ4BWVnQG/jYXNM6DB1TDoa4j+GxY+bvRCXzK20lekMZvWmjfWvsG+xH1M6j3JlE5KzsIeagYuN3bb1vi4+/DZtZ8xcsVInl3zLApjqJwQ9kASeEUkHoYfH4DT2+Hq/4NeL4Crm5GsoywdzkpYTrSqWGt5xqKO1WZmGwBq+9bmmrBrrBCtsGdlGbttC3zcfZjcezKP/fEYz655FhT0C5ckLmyfJPDy2vMrzB9pdFC7ew40te0PvDWXZyx8rNSsVA4lH6KWby161utp5aiFqJjSNjn4uPvw+XWfG0l8tVESvz78+qoOV4gykQReVrnZsOI1+PdTqNPe6GUe0MDsqEo0rNUwBi4YiJ+HHymZKdy/5H52xO+gQ60OPPnnk7i6uOLm4oabciv43VW54u7ijquy3Lds93T15Ls93xHoGUh8Rjw1fWpyNussD7V+yOyXKcRFytLk4OPuw+TrjJL4uNXjcFEusma9sGmSwMvi7EmYOwyOrTXGdl//Nrh5Fjx80bf9oiZfMbGDUbB3MPX96rMzYSderl6cSDtBDc8anM06S2JGIjl5OeTqXHLycgpu+fdzdS7Zednk5OVcdMxT6adwVa60r9meAK8Au6kyFaI4vu6+fH7d54xYPoKxq8byfo/3CxZPEcLWSAIvrcMrYe5DkH0ebv+6yE5phb/tD/3GmMBk+oNRVRllkXLycnh97evsTNiJQtHQvyEn008yf8D8Mifd3LxccnUup8+dpv/8/kTWiGTNiTUyr7RwGAVJ/I8RPLPqGSb0mEDvBr3NDkuI/5AE/tc7sOrd4h/vPg5c3eGvtyCkKdw50/hpJ87nnOeZVc+wKmYVI9qO4Of9P7M/eX+5h/i4urjiiiv1/esXzGxl68OFhCirah7VmHLdFB7941GeXvU0E3pOoHd9SeLCtriUvIuD6/UcvJpi3Bpcbdzy7z9zGE5ugr/ehNZ3wCN/2lXyTs5I5uFlD7M6ZjUvdXmJUe1GUcu3Fr7uvlYZ4lPLx3rHEsLW5CfxFkEteHrl0/x57E+zQxLiIpLAi3N8I0ztbkzEcvNHcNsXdrUM6Mm0k9z3233sTdjLhz0/5M6mdwLWHeJjb8OFhCgrPw8/pvSZQvOg5jy16ilWHl9pdkhCFJAq9HyFlwCd0ATOxUP1MHhomdHbvDQKVccXLDvyavULj/d41ijxV7L9Sft5bPljnM89zxd9v+DKmldW+jnLw9ZnFRMCLiTxR5c9ypiVY5jYcyI96vUwOywhJIEDF5YArdEA0mIh7Qy4e8M1/1f65A1GcrYkaLNWENt4eiNP/vkk3u7ezOg3g8iAyCo9f1nYw6xiQgD4e/gzte9Uhi8bzuN/Po5GFzxWkUmRhKgIqUIHYwnQ698x1u4+nwB9XjcmZ1k3xezIymT50eWMWD6CEJ8QvrvhO5tO3kLYG38Pf6b2mUqAVwAuuNC4RmM61OzAjgd2cG/ze2kV3MrsEIWTkQQOxlKfkX3AOxBqtoarnoQG3exqCdDZe2fz1MqnaB7UnJk3zKR2tdpmhySEw6nuWZ1p108DBQeTD5KUkURcehyLDi+SzpyiykkCB2Opz+PrIDACvCxt1sfW2sUSoFprPtn8CW+tf4se9XrwZd8vqe5ZveQnCiHKpVGNRgyKHISbcuNQyiGeW/Mc/Rv1l86cospJAgfo/rSx5GdGMui8C0uAdn/a7MguKycvh1f+fYUvd3zJ7ZG381HPj/B28zY7LCEc3oi2I9BoqrlVY/3p9eTl5aG1LvmJQliRdGKDC7Oq2dESoJdO0DKy7chyLUwihCi7/HXPE84n0LhGY77f+z2p2am82u1V3F3czQ5POAlJ4PlscAnQ4iRnJDPqz1HsiNvBS11eKhjjLYSoOvnrnn/R5wvmHpjL5K2TSTifwAc9P8DXveQ5I4oaRlmYDKMUJZEEXtRUqiaM3S6tk2kneXT5o5xMO8mHPT+UhRaEMEn+REYhPiE81vYxavrU5PW1rzP096FMvm5yiW3iF62dYIVhp/f/dj9bYrf8Z7sMb3NcksALjd22NZd+IPPHm7oqV76+/mubnaBFCGd0W+RtBHsH8/Sqp7l3yb1MuW4K4dXDq+z8LYNa0jKoJeM6jSvYNn7DeGlac2DSic2GtQxqyb3N76VDzQ40DWiKn7sf3m7e3NjwRkneQtig7mHd+brv16Rnp3Pfb/exLW5blZ17WKthLDi0gPEbxvPbkd9keJsTsLkErpTqp5Tap5Q6qJR61ux4zDS05VDmH5xPTGoM+5P2E+AVgJuLG//X4f/MDk0IUYzWIa357sbv8PPw4+GlD1fJ/OlpWWnMOzCPjJwMvtvzHTvidzBt5zQZ3ubgbKoKXSnlCkwC+gAxwEal1EKt9W5zI6s6p8+dZt2pdaw7tY4NpzaQlp1GWnYafu5+dKndBU83T/lA2hmZ89351Pevz7c3fMuoFaN48q8neaHzC5XS2TQ9O51Ze2cxfdd0UjJT6FqnK9vjtnNLxC08svwR5g+Yb/VzCtthUwkc6AQc1FofBlBKzQYGAA6bwFOzU1l+dDnrT61n3al1HD17FIBAr0A61+pMi6AWTNw8kXp+9Vh6dKl8IO2QzPnunIK8g5h2/TSeWvUUb6x7g9j0WEa1G2WVNumMnAx+3PcjX+/8msSMRK6qexWj242mVXArxm8Yz8PLHmZA4wHyZd/B2VoCrwscL3Q/BuhsUizlUlzHs/yeoOdzzrMldgvrTq1j/an17EnYg0bj4+ZDh1oduLPJnXSu3ZnIgEhclNHC8e2eb9mXtI8hzYfIB1IIO+Lj7sMn137CG2vfYOr2qZxJP8PLXV8u91jxrNws5h2Yx1fbvyL2fCyda3dmdLvRtAttV7DPsFbDOJJyRNq+nYCtJfBSUUoNB4YD1K1bl+joaHMDKqSeez3q1a3H4bTDALzZ9k0+2PMBcRlxDJk/hD1n95Cjc3BTbjT1b8rA0IF0qdOFJn5NcHOx/DnOwrGzxwqOGeAaQLprOn2q96nwa83IyACwyjWz5rGsrbSxJSQkVEE0ts+s90VJ198R3mMAQ+sMxSPLgzkH53A88TjjWozD2/XCrIklHSsnL4cVZ1Yw5+gc4jLjaOHfgv81+R+ta7SG9P8+79nIZ0k7k0YaacXGJO99c1nj+ttaAj8B1Ct0P8yy7SJa6y+ALwA6dOigw8PDqyS40vhf6P8YuGAggV6BJGcmM+TfIaTnpAPQLLAZ9zS/hy51unBF6BX4uPsQHR1NSfH77fWjmU8z2jcpw9KmxfDa6wVQ4jmr+ljWVpbYbDH+qmbm++Jy+znKewzgxYYv0nR/U95c9yav7XmNSb0nEeQddNlj5eTlsPjwYqZsm0JMWgxtgtvwZvc36Vq7q1Wq4m3xujqTil5/W0vgG4FIpVRDjMQ9GBhibkhlE+ITQtc6XVkavRQPFw/q+9entm9t3rz6TQK9As0OT4gC0rmu6t3R5A6CvYIZu3os9/12H1Oum0J9//r/2S9P5/H7kd/5fNvnRJ+Npnlgcyb1nsQ1da+Rcd2igE0lcK11jlJqNLAUcAWmaa13mRxWmZzNOsvWM1sBCPcP58z5M3zZ90tJ3sLm2GrnOkf/YtGrfi++uv4rRq8YTf/5/cnVuQWP5b9OL1cvMnIzaFyjMRN7TuTa+tdK4hb/YVMJHEBrvQRYYnYc5aG15tV/XyUhI4EAzwAOpRySjmfCKVgz6drqFwtrahvSlm9v+JbBiweTkZNBw+oNAXBzcWNv4l48XD1446o36Bvet6AzqxCXsrkEbs/mHpjL8qPLGXPlGFYcXUH02WjpCSqcgjMkXWsLrx7OjH4zuPPXOzmYfBCAOr518Hbz5uf+P1PTt6bJEQpbJ1/trORA0gHGbxhPtzrdeLDlgwULHUjpWwhRnKaBTbmzyZ14u3rTwL8BPev1ZFCTQZK8RalICdwKzuecZ+zqsVRzr8ZbV79lc1Vejt6mKIQ9G95mOD/t/wlfN18WH1kskzWJUpMEbgXvb3yfg8kHmXrdVJsscTtL9aZ8URH2KMQnhCDvIJmsSZSZJPAKWhq9lJ/2/8SwVsPoVreb2eE4NWf5oiIcTy2fWmTkZEifGVEmtlXXa2dOpJ3gtX9fo01wG0a3H212OEIIOyV9ZkR5SAIvp+y8bMatHodGM777+HLPbSyEEEKUh1Shl9PkrZPZFreN97u/T5hfmNnhCCGEcDKSwMth3al1fL3ja26PvJ1+DfuZHU6Vko5iQghhG5w+gReVkAq7NCElnE/guTXP0bB6Q8Z1GlcVIdoU6SgmhBC2wekTeOGENPT3oQBM7ze9yH3zdB4v/vMiZzPPMrXPVLzdvIvcTwjhmKQGStgSp0/gZfHt7m/5+8TfvNj5RZoENDE7HCFEFZMaKGFLpBd6Ke2K38XEzRPpXb83dza90+xwhBBCODlJ4KWQlpXGM6ufIdg7mNe6vSbL+gkhhDCdVKGXQGvNG+ve4GTaSab3m051z+pmhySEEEJIAi/JgkMLWHJkCaPbjaZ9aHuzwxFCOAjpECcqShL4ZRxJOcLb69+mY62OPNz6YbPDEUI4EOkQJypK2sCLkZmbyTOrnsHT1ZN3rn4HVxdXs0MSQgghCkgJvBgfRn3IvqR9TOo9iZq+Nc0ORwghhLiIJPAi/HXsL37Y+wP3Nr+X7mHdTYlB2seEEEJcjiRw4P7f7mdL7JaC+1FnogDYEb+jTMexZtKV9jEhhBCXIwkcaBnUkpZBLdmTsId9SfvI03n0adAHf0//Mh1Hkq4QQoiqIp3YgGGthrHw0EJiUmNIy07jiSueYGXMSoa1GmZ2aEIIIUSRJIEDIT4hXB9+PbHnYwnyCiImNYb+jfoT7B1sdmhCCCFEkSSBWzzW9jFclAvBXsEsOrxISt9CCCFsmiRwixCfEIK8gziYclBK30IIIWyedGIrpJZPLTJyMqT0LYQQwuZJCbwQD1cPmgQ0kdK3EEIImyclcCGEEHbh0jk78rUPbc/MG2aaEJG5JIELIYSwC/lzdozrNK5g2/gN41FKmRiVeaQKXQghhF3In7MjLj0OgLj0OKceNSQJXAghhF0I8Qmhf6P+TNs5DYBpO6c59aghSeBCCCHsRn4pfE/CHqcufYMkcCGEEHYkvxT+8LKHnbr0DdKJTQghhJ0Z1moYR1KOOHXpGySBCyGEsDMhPiFM6TPF7DBMJ1XoQgghhB2SBC6EEELYIaevQp+8dTKfb/v8om2tZ7Qu+P2xto8xst3Iqg5LCCGEuCynT+Aj242UBC2EEMLuSBW6EEIIYYckgQshhBB2SBK4EEIIYYckgQshhBB2SBK4EEIIYYckgQshhBB2SBK4EEIIYYckgQshhBB2SBK4EEIIYYckgQshhBB2qEJTqSqlXgUeAeIsm57XWi+xPPYc8BCQCzyhtV5q2d4P+BhwBb7SWr9r2d4QmA0EAZuA+7TWWRWJTwghhKgK9/92P1tit/xne/vQ9sy8YWalnNMac6F/pLWeUHiDUqoFMBhoCdQB/lBKNbE8PAnoA8QAG5VSC7XWu4HxlmPNVkpNwUj+F68yIoQQQtiglkEtaRnUknGdxhVsG79hPEqpSjtnZVWhDwBma60ztdZHgINAJ8vtoNb6sKV0PRsYoIxXeC0w1/L8GcDASopNCCGEsKphrYax8NBC4tKNCum49DgWHV7EsFbDKu2c1iiBj1ZK3Q9EAU9prZOAusC6QvvEWLYBHL9ke2eMavNkrXVOEfv/h1JqODAcoG7dukRHR1vhZZgjISHB7BCcmlx/c8n1N49ce+vrGdKTif9O5JHGj/DlwS/pGdyTtDNppJH2n32tcf1LTOBKqT+AWkU89AJGFfcbgLb8/ACovK8bFlrrL4AvADp06KDDw8Mr+5SVyt7jt3dy/c0l1988cu2t63+h/2PggoHc2/5eVsWvYv6A+QR7Bxe7f0Wvf4kJXGt9XWkOpJT6EvjVcvcEUK/Qw2GWbRSzPQGooZRys5TCC+8vhBBC2LwQnxD6N+rPw8seZkDjAZdN3tZQoTZwpVTtQndvBXZafl8IDFZKeVp6l0cCG4CNQKRSqqFSygOjo9tCrbUG/gIGWZ7/ALCgIrEJIYQQVW1Yq2G0Dm5dqW3f+SraBv6eUqodRhV6NPAogNZ6l1LqR2A3kAOM0lrnAiilRgNLMYaRTdNa77IcaxwwWyn1JrAF+LqCsQkhhBBVKsQnhCl9plTJuSqUwLXW913msbeAt4rYvgRYUsT2wxi91IUQQghRApmJTQghhLBDksCFEEIIOyQJXAghhLBDksCFEEIIOyQJXAghhLBDksCFEEIIOyQJXAghhLBDksCFEEIIO6SMWUztl1IqDjhqdhwVEAzEmx2EE5Prby65/uaRa2+uslz/BlrrkEs32n0Ct3dKqSitdQez43BWcv3NJdffPHLtzWWN6y9V6EIIIYQdkgQuhBBC2CFJ4Ob7wuwAnJxcf3PJ9TePXHtzVfj6Sxu4EEIIYYekBC6EEELYIUngQgghhB2SBG4ipVS0UmqHUmqrUirK7HgcnVJqmlIqVim1s9C2QKXUcqXUAcvPADNjdGTFXP9XlVInLJ+BrUqpG82M0VEppeoppf5SSu1WSu1SSj1p2S7v/0p2mWtf4fe+tIGbSCkVDXTQWstkClVAKdUdSANmaq1bWba9ByRqrd9VSj0LBGitx5kZp6Mq5vq/CqRprSeYGZujU0rVBmprrTcrpfyATcBA4EHk/V+pLnPt76SC730pgQunobVeDSResnkAMMPy+wyMD5aoBMVcf1EFtNantNabLb+nAnuAusj7v9Jd5tpXmCRwc2lgmVJqk1JquNnBOKmaWutTlt9PAzXNDMZJjVZKbbdUsUsVbiVTSoUD7YH1yPu/Sl1y7aGC731J4Oa6Wmt9BXADMMpSxShMoo32JGlTqlqfA42AdsAp4ANTo3FwSqlqwDzgf1rrs4Ufk/d/5Sri2lf4vS8J3ERa6xOWn7HAL0AncyNySmcsbVT5bVWxJsfjVLTWZ7TWuVrrPOBL5DNQaZRS7hgJ5Hut9c+WzfL+rwJFXXtrvPclgZtEKeVr6dCAUsoX6AvsvPyzRCVYCDxg+f0BYIGJsTid/ORhcSvyGagUSikFfA3s0Vp/WOghef9XsuKuvTXe+9IL3SRKqQiMUjeAG/CD1votE0NyeEqpWUBPjGX8zgCvAPOBH4H6GMvS3qm1lo5WlaCY698TowpRA9HAo4XaZIWVKKWuBtYAO4A8y+bnMdpi5f1fiS5z7e+mgu99SeBCCCGEHZIqdCGEEMIOSQIXQggh7JAkcCGEEMIOSQIXQggh7JAkcCGEEMIOSQIXQggh7JAkcCGEEMIO/T975myWnABDWAAAAABJRU5ErkJggg==\n",
"text/plain": [
"