{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Problem 14.4" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Fundamentals of Solar Cells and Photovoltaic Systems Engineering**\n", "\n", "**Solutions Manual**\n", "\n", "**Problem 14.4**\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Determine the dominant frequencies for solar PV and wind power fluctuations by calculating the Fast Fourier Transform (FFT) of the solar PV and wind power time series corresponding to the country in Problem 14.3 and plotting the power spectra.**" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We downloaded solar PV and wind time series for Spain from https://model.energy/ and saved it in the data folder." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We will use the packages [pandas](https://pandas.org/), [numpy](https://numpy.org/), and [matplotlib.pyplot](https://matplotlib.org/stable/index.html) to plot the results" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import pandas as pd\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import matplotlib.gridspec as gridspec" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "First, we retrieve the time series. " ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
onwindsolar
2011-01-01 00:00:000.0210.0
2011-01-01 01:00:000.0210.0
2011-01-01 02:00:000.0200.0
2011-01-01 03:00:000.0180.0
2011-01-01 04:00:000.0160.0
\n", "
" ], "text/plain": [ " onwind solar\n", "2011-01-01 00:00:00 0.021 0.0\n", "2011-01-01 01:00:00 0.021 0.0\n", "2011-01-01 02:00:00 0.020 0.0\n", "2011-01-01 03:00:00 0.018 0.0\n", "2011-01-01 04:00:00 0.016 0.0" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "CF=pd.read_csv('data/time-series-SPAIN.csv',\n", " sep=',', index_col=0)\n", "CF.head()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We calculate the Fast Fourier Transform (FFT) of the solar PV capacity factor time series. " ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "C:\\Users\\34620\\.conda\\envs\\env_pv_textbook\\lib\\site-packages\\ipykernel_launcher.py:7: RuntimeWarning: divide by zero encountered in double_scalars\n", " import sys\n" ] } ], "source": [ "t_sampling=1 # sampling rate, 1 data per hour\n", "x = np.arange(1,8761, t_sampling) \n", "y = CF['solar']\n", "n = len(x)\n", "y_fft = np.fft.fft(y)/n #n for normalization \n", "frq = np.arange(0,1/t_sampling,1/(t_sampling*n)) \n", "period = np.array([1/f for f in frq]) " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We plot the power spectra for solar CF time series." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(740, 0.95, 'month')" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAzoAAAFzCAYAAAD2eXw5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAA9hAAAPYQGoP6dpAABgKUlEQVR4nO3dd3gU5fo38O9syqYHQiAhkCbNQEKRgIQiIBhFwS5YDlgA5YceFBCVwzmivig2ONjAgsDRI4oNLAeRoPQiNSgkEAiEBEgI6RDSZ94/HnaTTSG7IZvZmf1+ritXdmdnZ+/dZ57duecpIymKooCIiIiIiEhHDGoHQERERERE1NyY6BARERERke4w0SEiIiIiIt1hokNERERERLrDRIeIiIiIiHSHiQ4REREREekOEx0iIiIiItIdJjpERERERKQ7rmoHYA1ZlnH27Fn4+vpCkiS1wyEiIiIiIpUoioILFy4gJCQEBkPD7TaaSHTOnj2L0NBQtcMgIiIiIiIHkZGRgY4dOzb4uCYSHV9fXwDizfj5+akcDRERERERqaWoqAihoaHmHKEhmkh0TN3V/Pz8mOgQOSlZlpGamgoA6NSp0xWbqomaE/c9/WLZElnPEetLY0NaNJHoEBEBQFlZmdohkJPivqdfLFsi62mtvqifihERERERETUzJjpERERERKQ7uum6pigKKioqUFlZqXYoVIurqyvc3Nw4NTgRERERtRhdJDplZWVIS0vDxYsX1Q6FGuDj44OIiAgYjUa1QyEiIiIiJ6D5REeWZSQlJcHV1RWRkZEwGo1sOXAgiqKgrKwMZ86cweHDh9GjRw8mO0RERERkd5pPdEpLSyHLMiIjI+Hj46N2OFQPb29vuLu74+jRo/jxxx8xYsQIBAQEqB0WaZCbm5vaIZCT4r6nXyxbIutprb7YPBnBli1bMGbMGISEhECSJKxZs6bR52zevBl9+/aFh4cHrrnmGnz44YdNifWKHGEub2qYqXyysrLw888/o6ioSOWISGsMBgO6deuGbt26sb5Ti+K+p18sWyLrabG+2BxlcXExevXqhffff9+q9U+ePIlbb70VQ4YMwYEDB/CPf/wD06ZNw3fffWdzsKR9wcHBOHv2LLKystQOhYiIiIh0zOaua6NGjcKoUaOsXv/DDz9EWFgYFi1aBACIiorC3r178fbbb+Oee+6x9eVJ40xnAEpKSlSOhDTp5BdA4SGg+3OAe2u1oyEiIiIHZvd2p507dyI+Pt5i2c0334y9e/eioqKi3ueUlZWhqKjI4o+qpaWlQZIkJCYmqh0KUYuR8w8j9bwLUst7Qd43q0nbyM3NxdSpU/HZZ581c3SkZ7IsIzU1FampqZBlWe1wrtrOnTsxdepU7Ny5U+1QVKe3sq3t559/xtSpU5GSkqJ2KKQDWqwvdk90srKyEBQUZLEsKCgIlZWVyMnJqfc58+fPh7+/v/kvNDTU3mHSZStWrIAkSea/9u3bY+zYsTh58iTKy8sRGBiIefPm1fvc+fPnIzAwEOXl5S0cNTmF89tQ4hWNEq9oIG2l2tGQkykpKWFLtE5puWx58oZamtbqS4uMJKo93bOiKPUuN5k9ezYKCwvNfxkZGXaP0dlcKRnx8/NDZmYmzp49i5UrVyIxMRG33347XFxc8Le//Q0rVqwwl2FNy5cvx/jx4+Hu7m7P0ImIiIiIGmX3RCc4OLjOwPPs7Gy4urqiTZs29T7HaDTCz8/P4k9vvv32W8TExMDT0xNt2rTByJEjUVxcDEA0Db7yyivo2LEjjEYjevfujXXr1jW4raqqKkycOBGRkZHw9PREt27d8M4771is88gjj+DOO+/E/PnzERISgq5duza4PUmSEBwcjPbt22P48OGYO3cuDh06hOPHj2PixIlITU3Fli1bLJ6zdetWHDt2DBMnTryKT4WIiIiIqHnY/To6cXFx+OmnnyyWrV+/HrGxsfabi3tdLFCiwqxensHALXsbXS0zMxMPPPAA3nzzTdx11124cOECtm7dam4leeedd7BgwQJ89NFH6NOnD5YtW4bbb78dhw8fRpcuXepsT5ZldOzYEV9//TUCAwOxY8cOPP744+ZuZya//fYb/Pz8kJCQUG+LTINvy9MTAFBRUYGYmBj069cPy5cvx9ChQ83rLFu2DP3790d0dLTV2yWyF1mWsWHDBmzfvh35+flo3bo14uLiEBsbW2fdo0ePYvfu3UhNTUVhYSEA0b128ODBGDx4sHm90tJSzJ49GwEBAfjXv/5V72vOmTMHVVVVeO211+DqqvnLlFETFRcX47nnnkNMTAymTJliXp6WloY333wTAPDqq6+idevqCTXefPNNZGZm4u2334aLiwsA4ODBg9i0aRPS09NRUVGBtm3bYsCAARgxYkS9U7vaun5teXl5eO+995CXl4dJkyYhJibmaj8KzUpJScGiRYtw6623IioqCj/++CPS09Ph5uaGPn364J577oG7uzsOHz6MX375BadPn4bRaMSAAQNwxx13WHzeVVVV2LRpE3bt2oXs7Gy4uLggNDQUI0eOrPMZ79y5E59//jnGjx+P1q1b43//+x8yMjLg5uaG6Oho3HPPPeZrBprWBYBdu3Zh165d5u0888wzdU5o7tu3D+vXr0dWVhY8PT1x3XXX4c4772QvDNI1m3+JL168iOPHj5vvnzx5EomJiQgICEBYWBhmz56NM2fOmPuLTpkyBe+//z5mzJiByZMnY+fOnfj000/x5ZdfNt+7qK0kCyg5Y7/tX6XMzExUVlbi7rvvRnh4OABYfNm9/fbbeP7553H//fcDAN544w1s3LgRixYtwgcffFBne25ubnj55ZfN9yMjI7Fjxw58/fXXFomOt7c3li5datOX2unTp/HWW2+hY8eO5i/Nxx57DM8++yzef/99+Pj44OLFi/jmm2+wcOFC2z4IIjv54osvsHPnTrRp0wZDhw5FRUUFfv/9d5w4caLOugkJCTh//jwiIiLQunVrXLp0CUlJSVi5ciXOnTtnnh3Sw8MDsbGx2L59O1JTU9GpUyeL7Rw6dAiFhYUYMWIEkxwn5+3tjQ4dOuDYsWOQZdl80Hvs2DHzOikpKbj++usBiCQ6PT0dUVFR5iTnhx9+wK+//opWrVqhT58+8PDwwPHjx7F69WqkpaVh8uTJFq9p6/q1ZWZm4r333kN5eTmmTZtWZ/92VmlpaUhISEBUVBQGDx6MlJQUbN26FaWlpejVqxf+85//oGfPnoiIiMChQ4eQkJAADw8P8+y0iqLg008/RWJiItq1a4cbbrgB5eXl2L9/P5YsWYL77rsPw4cPr/O6f/31Fw4dOoSYmBgMGTIEx48fxx9//IGcnBzMnDkTANCxY0cMHz4cGzduRMeOHdGzZ0/z82v3mNm8eTOSkpLQs2dPdOnSBUlJSdi0aROKi4vx6KOP2vETJFKXzb/Ge/futaiUM2bMAAA8/PDDWLFiBTIzM5Genm5+PDIyEmvXrsX06dPxwQcfICQkBO+++659p5b2DLbftpvhdXv16oURI0YgJiYGN998M+Lj43HvvfeidevWKCoqwtmzZzFo0CCL5wwaNAgHDx5scJsffvghli5dilOnTqGkpATl5eXo3bu3xToxMTFWJTmFhYXw8fGBoii4dOkSrrvuOnz//ffm5z7wwAOYMWMGVq1ahYkTJ2LVqlVQFMWcmBGpKSUlBTt37kTHjh0xc+ZMGI1GAMAtt9yC1157rc76999/PwIDAy2WVVVV4YMPPsDGjRsxfPhwBAQEAAAGDx6M7du3Y8eOHXUOBHfs2AEAdeouOaeuXbvi999/x+nTpxEWFgZA7JshISG4cOGCRaJz/PhxyLJsPpmUnJyMX3/9FT169MDkyZPN372KouCrr77C1q1bceDAAfTp06dJ69d24sQJLFmyBK6urpgxYwZCQkLs+tloSVJSEp544gn06tULgPhueP3117F3714kJydj+vTpiIiIAACMHj0ac+fOxe+//474+Hi4uLhg9+7dSExMRJcuXfD3v//dfBLklltuweuvv47vv/8eMTExdb6D/vzzT0yfPt38PSPLMt555x0cO3YMJ0+eRGRkJEJDQ+Hl5WVOdEaPHt3g+zhy5AheeOEF8+RQ5eXleO2117B3717cddddaNWqVTN/ckSOweZEZ9iwYVfs9rRixYo6y4YOHYr9+/fb+lJNZ0X3MTW5uLggISEBO3bswPr16/Hee+9hzpw5+OOPP8xnYeqbwKGhyRu+/vprTJ8+HQsWLEBcXBx8fX3x1ltv4Y8//rBYz9vb26r4fH19sX//fhgMBgQFBdV5nr+/P+69914sX74cEydOxPLly3HvvffqciwVORaXyrxG1zHt96NGjTInOQDQqlUrDB8+vE5X2toHGICoo0OGDMGRI0eQkpKCAQMGAADCw8MRGhqK/fv347777oOHhwcAcXLg0KFD6NSpE4KDVTrRQnZlammxlinROXr0KMLCwlBVVYXU1FTExcWhsLAQR48eNa9rmvrX1DV58+bNAMRJpZonpyRJwp133olt27Zh79695sTF1vVrOnToEJYuXYpWrVrh73//e4NjZ/XsSmXbpUsXc5JjWrdPnz44c+YMYmJizEkOIFp9o6OjsXPnThQUFKBNmzbm7mR33XWXRUtv69atceONN+KHH37Anj176lyfsF+/fhYnUwwGAwYMGIBjx44hLS0NkZGRNr3H4cOHW8yA6+7ujtjYWKxduxbp6elMdMhqtn4Xqo39K1QiSRIGDRqEQYMG4cUXX0R4eDhWr15tPpu2bds23HDDDeb1d+zYgf79+9e7ra1bt2LgwIGYOnWqeVlqamqTYzMYDOjcufMV15k4cSKGDRuGn3/+Gdu3b6/3TDlRczJIMqKShja63unTpwGg3n24vu44paWl2LBhAw4ePIicnByUlZVZPF5QUGBxf8iQIVi5ciX27t1rHsOza9cuyLLM1hydMhgMiIqKsuk5nTt3hiRJSElJwU033YT09HSUlpaia9euKCwsxP79+5GTk4PAwECkpKTAw8PD3PJz8uRJGI1GcythbW5ubhaT/Ni6vsn+/fuRnJyM0NBQPPnkk+axH86ksbKt7/IW/v7+AETXsYYeMyU6pvE1NRMiE1MLnuk7q7HXNY3pasrUvs29PXJOTfkuVBsTHRX88ccf+O233xAfH4927drhjz/+wPnz5807z6xZszB37lx06tQJvXv3xvLly5GYmIgvvvii3u117twZn332GX799VdERkbi888/x549e2w+42OLoUOHonPnzpgwYQI6d+5skZQRqam0tBSSJNV70Fa71bGyshL//ve/kZGRgdDQUPTv3x/e3t4wGAzIy8vDrl27UFlZafGc2NhYfPfdd9i+fbs50dm5c6d5cC8RAHh5eSE0NBSpqamoqqrC0aNHIUkSunTpggsXLgAQLTne3t7IyMhAjx49zGN5iouLIcsy1q5d2+D2a14iwNb1TU6ePAlZltG5c2enTHKsYWq1rclUTld6rKqqCoD4Pqo56URNpu+j+hIN0yRA9W27KRdqbO7tEWkFEx0V+Pn5YcuWLVi0aBGKiooQHh6OBQsWmJuup02bhqKiIsycORPZ2dno3r07fvzxx3pnXAPEhA+JiYkYN24cJEnCAw88gKlTp+KXX36x6/t47LHH8I9//AOzZjXtKvVE9uDh4QFFUXDx4kX4+vpaPFZUVGRx/88//0RGRgYGDRqEhx56yOKxvXv3WsxiVHP7/fr1w7Zt23D69GlcunQJ2dnZuOGGGzh7EVno2rUr0tPTkZ6ejmPHjqFjx47w8vKCl5cX/P39kZKSYh4P2a1bN/PzTAelb731llWvY+v6Jrfffjv+/PNP/PbbbzAYDLjrrrtsej41zsPDo873jolpeX0JExE1jxa5YChZioqKwrp165CdnY3S0lIcPXoUTz31lPlxg8GAF198EadPn0Z5eTkSExNxyy23mB+PiIiAoijmyQaMRiOWL1+OgoIC5OfnY/HixZg/fz4SExPNz1mxYgXWrFnTaGyPPPJIna46DZk9ezYURcHs2bOtWp/oasiKASeuWYYT1yyDLBkbXM/UnaTm7JAmtbt0nj9/HgDqnUa3vuebmFpyduzYwUkInIAsyzhx4gROnDhh09lv08mp5ORknDhxwiKZ6dq1K1JSUuqMzwHEd3xxcTGys7Oteh1b1zdxc3PDE088gR49eiAhIQHff/+9Tc/Xg6aWrbVCQ0NRUVGBtLS0Oo+ZZuGrrwuctUxjd9kqQy3B3vXFHpjoEJFmXPLph0s+/XClry7TTFa//PKLxXibgoICbNy40WJd02xqtROglJQUbN++vcHXCAsLQ1hYGHbv3o0DBw4gLCys3j7wpB+XLl3CpUuXbHpO586dYTAYsGnTJpSVlVlc16Rr164oKCjA7t274enpaXGwO2zYMADA559/josXL9bZbmFhITIzM5u8fk1ubm54/PHHER0djQ0bNjhlstOUsrWW6fvohx9+MHdnA8T30e+//w6DwdDg+FtreHl5QZIk5OfnX3WsRNawZ32xB3ZdIyJd6dq1K+Li4rBz507MmzcPvXv3RmVlJfbt22e+1oVJTEwM2rRpg4SEBJw9exYhISE4d+4cDh06hF69euHAgQMNvs7gwYOxcuVKAMDAgQPt/r5Iezw9PREWFoa0tLQ6k7yYWncuXryInj17WlxgskePHhg1ahR++eUXzJ07Fz169EBAQIC51SY1NRVjxoxB+/btm7R+baZk5+OPP8aGDRugKIp9LwHhRK6//nokJibizz//xLx58xATE4OysjLs378fxcXFuOeee+qd+dFaHh4eCA8Px/Hjx/HZZ5+hXbt2kCQJ/fr1M5/IIXJmTHSISHceeughtGvXDtu3b8fmzZvRqlUr3Hjjjejbt69FouPh4YGnn34aq1evxrFjx3Ds2DG0b98ejzzyCPz8/K6Y6MTGxmLVqlVwcXFBv379WuJtkQZ17doVaWlpCAsLsxiLERgYiICAAOTl5dW5gj0AjBkzBl26dMHGjRtx5MgRlJSUwNvbG23atMFtt91WZ5+zdf3aXF1d8fjjj+OTTz7Bb7/9BkVRcO+99zbPh+DEJEnC5MmTsXHjRuzatQubNm2Cq6srQkNDMWLECIuLfDbVww8/jG+//RYHDx5EaWkpFEVBZGQkEx0iAJJypYviOIiioiL4+/ujsLCwzqxJly5dQnJyMqKiouDl5aVShNQYUzmlpaXh2LFjGDVqlMW1CYgaI6d8hKRyMQ6m+1/9YXhA3abztLQ0vPnmm4iLi8P48eNVjYXsS5ZlJCUlAQC6d+9u0fpC2sayJbKeI9WXK+UGNemmRmtlUJSzYvmQ3mzYsAGAuK4OEREROR7Nd10zTed68eJFXgfAgZkGyFZUVKgcCVHT5eXlYc+ePcjMzMT+/fvRvXv3ei8ESEREROrTfKLj6uqKwMBAnDlzBgDg4+PDpmcHIssyLl68iDNnzqCgoIAtO3RVJFnd7mo5OTn44YcfYDQa0bNnTzzwwAOqxkMtxzSNL+kPy5bIelqrL5pPdAAx1assy+ZkhxxPQUEBzp07Z55es76rNBNdiUGS0ePQ9arG0LVrVyxevFjVGKjlGQwG9OjRQ+0wyA5YtkTW02J90UWiI0kSIiMjcerUKRw6dAje3t48kHYQiqKgvLwcsiyjqqoK2dnZ6NixI4KDg9UOjYiIiIh0TBeJjsmQIUMgyzIOHz6M/Px8aGBCOafi4uKCsLAw3HrrrVecIYOIiIiI6GrpKtFxcXHBsGHD0K9fP01dtdVZuLm5wdfXF25ubmqHQhokKxLSIz4AAISdmq6fKSPJ4cmyjPT0dACiqzTHgeoHy5bIelqsL7pKdADRf9DX1xe+vr5qh0JEzUrCRb8bLt92UTUScj6mmSNJf1i2RNbTWn1x/FSMiIiIiIjIRkx0iIiIiIhId5joEBERERGR7jDRISIiIiIi3WGiQ0REREREusNEh4iIiIiIdEd300sTkT4ZJBnRf8aoHQY5IYPBgOjoaLXDIDtg2RJZT4v1hS06RERERESkO0x0iIiIiIhId9h1jYg0QVYknA5bAADomDGbZ2moxciyjNOnTwMAOnbsCIOBe59esGyJrKfF+sJEh4g0QkJRq3hxM+Of6oZCTqeoqEjtEMhOWLZE1tNafXH8VIyIiIiIiMhGTHSIiIiIiEh3mOgQEREREZHuMNEhIiIiIiLdYaJDRERERES6w0SHiIiIiIh0h9NLE5EmSJDR/a/+4rZSonI05EwkSUL37t3Nt0k/WLZE1tNifWGiQ0SaIElMcEgdkiRp5kedbMOyJbKeFusLu64REREREZHusEWHiDRBViSc7TgPABBy5mWepaEWI8syzp49CwAICQmBwcC9Ty9YtkTW02J9cfwIiYgAABIKAu5AQcAd4DkaamkFBQUoKChQOwyyA5YtkfW0Vl+Y6BARERERke4w0SEiIiIiIt1hokNERERERLrDRIeIiIiIiHSHiQ4REREREekOEx0iIiIiItIdztFKRJogKVW49vANl2+XqBwNORNJknDttdeab5N+sGyJrKfF+tKkFp3FixcjMjISHh4e6Nu3L7Zu3XrF9b/44gv06tULXl5eaN++PR599FHk5uY2KWAick6SBLhW5cO1Kh/a+HolvZAkCa6urnB1ddXMjztZh2VLZD0t1hebE51Vq1bhmWeewZw5c3DgwAEMGTIEo0aNQnp6er3rb9u2DRMmTMDEiRNx+PBhfPPNN9izZw8mTZp01cETkRPRyJcqEREROQabE52FCxdi4sSJmDRpEqKiorBo0SKEhoZiyZIl9a6/a9cuREREYNq0aYiMjMTgwYPxxBNPYO/evVcdPBE5D1mRcDZkDs6GzIEsuakdDjkRWZZx9uxZnD17FrIsqx0ONSOWLZH1tFhfbEp0ysvLsW/fPsTHx1ssj4+Px44dO+p9zsCBA3H69GmsXbsWiqLg3Llz+Pbbb3Hbbbc1+DplZWUoKiqy+CMiZychL/B+5AXeDw4vpJaWl5eHvLw8tcMgO2DZEllPa/XFpkQnJycHVVVVCAoKslgeFBSErKysep8zcOBAfPHFFxg3bhzc3d0RHByMVq1a4b333mvwdebPnw9/f3/zX2hoqC1hEhERERGRk2vSZAS1ByApitLgoKSkpCRMmzYNL774Ivbt24d169bh5MmTmDJlSoPbnz17NgoLC81/GRkZTQmTiIiIiIiclE39PwIDA+Hi4lKn9SY7O7tOK4/J/PnzMWjQIMyaNQsA0LNnT3h7e2PIkCGYN28e2rdvX+c5RqMRRqPRltCIiIiIiIjMbGrRcXd3R9++fZGQkGCxPCEhAQMHDqz3OZcuXYLBYPkyLi4uAERLEBERERERUXOzuevajBkzsHTpUixbtgzJycmYPn060tPTzV3RZs+ejQkTJpjXHzNmDL7//nssWbIEJ06cwPbt2zFt2jT0798fISEhzfdOiIiIiIiILrN56qJx48YhNzcXr7zyCjIzMxEdHY21a9ciPDwcAJCZmWlxTZ1HHnkEFy5cwPvvv4+ZM2eiVatWuPHGG/HGG28037sgIiIiIiKqQVI00H+sqKgI/v7+KCwshJ+fn9rhEJEKlJQlqDg4DwDgVpEJ6UFtzOFP2qcoCioqKgAAbm5umrkiODWOZUtkPUeqL9bmBrwYBRFpgiQB7hVn1Q6DnJAkSXB3d1c7DLIDli2R9bRYX5o0vTQRUYtz/MZnIiIiciBs0SEiTZAVCdntZwAA2mW9y7M01GJkWUZ2djYAoF27dnVmEiXtYtkSWU+L9cXxIyQiAgDJgJy2jyKn7aMA3NSOhpxMTk4OcnJy1A6D7IBlS2Q9rdUXJjpERERERKQ7THSIiIiIiEh3mOgQEREREZHuMNEhIiIiIiLdYaJDRERERES6w0SHiIiIiIh0h9fRISJNkCCj89E7xW2lVN1gyKlIkoTOnTubb5N+sGyJrKfF+sJEh4g0QZIAj7JUtcMgJyRJEjw8PNQOg+yAZUtkPS3WF3ZdIyIiIiIi3WGLDhFpgqxIOB/0fwCAttmf8CwNtRhZlnH+/HkAQNu2bWEwcO/TC5YtkfW0WF+Y6BCRNijA+aCpAIC22SvUjYWcTs0fd9IXli2R9bRWXxw/FSMiIiIiIrIREx0i0gaNzPBCREREjoGJDhERERER6Q4THSIiIiIi0h0mOkREREREpDtMdIiIiIiISHc4vTQRaYIEGdccu1/cVspUjoaciSRJuOaaa8y3ST9YtkTW02J9YaJDRJogSYBXyWG1wyAnJEkSvLy81A6D7IBlS2Q9LdYXdl0jIiIiIiLdYYsOEWmCrEjIbfsIAKBNzn95loZajCzLyM3NBQC0adMGBgP3Pr1g2RJZT4v1hYkOEWmEhHPtZwIA2uSsUjkWcjbnzp0DIH7cSV9YtkTW01p9cfxUjIiIiIiIyEZMdIhIGxRF7QiIiIhIQ5joEBERERGR7jDRISJt0Mic/UREROQYmOgQEREREZHuMNEhIiIiIiLd4fTSRKQJEmREpD4qbitlKkdDzkSSJERERJhvk36wbImsp8X6wkSHiDRBkgCf4r1qh0FOSJIk+Pj4qB0G2QHLlsh6Wqwv7LpGRERERES6wxYdItIERZGQ1+Z+AEBA7rfQRqM56YGiKMjLywMABAQEaKbLBjWOZUtkPS3WFyY6RKQJigJkdpgDAGid9wMTHWoxiqIgMzMTANC6dWtN/LiTdVi2RNbTYn1h1zUiIiIiItIdJjpEpA0aOHNEREREjoOJDhERERER6Q4THSLSBkVROwIiIiLSECY6RERERESkO0x0iEgbOEaHiIiIbMDppYlIEyTICD85VdxWylWOhpyJJEkIDw833yb9YNkSWU+L9YWJDhFpgiQBvhe2qh0GOSFJkuDr66t2GGQHLFsi62mxvjSp69rixYsRGRkJDw8P9O3bF1u3Xvngo6ysDHPmzEF4eDiMRiM6deqEZcuWNSlgIiIiIiKixtjcorNq1So888wzWLx4MQYNGoSPPvoIo0aNQlJSEsLCwup9ztixY3Hu3Dl8+umn6Ny5M7Kzs1FZWXnVwROR81AUCQWt7wAAtMr/H7TRaE56oCgKCgoKAACtWrXSTJcNahzLlsh6WqwvNic6CxcuxMSJEzFp0iQAwKJFi/Drr79iyZIlmD9/fp31161bh82bN+PEiRMICAgAAERERFxd1ETkdBRIOBM6DwDgX7CeiQ61GEVRcObMGQCAv7+/Jn7cyTosWyLrabG+2NR1rby8HPv27UN8fLzF8vj4eOzYsaPe5/z444+IjY3Fm2++iQ4dOqBr16549tlnUVJS0uDrlJWVoaioyOKPiIiIiIjIWja16OTk5KCqqgpBQUEWy4OCgpCVlVXvc06cOIFt27bBw8MDq1evRk5ODqZOnYq8vLwGx+nMnz8fL7/8si2hEZHe8YKhREREZIMmTUZQu6lKUZQGm69kWYYkSfjiiy/Qv39/3HrrrVi4cCFWrFjRYKvO7NmzUVhYaP7LyMhoSphEREREROSkbGrRCQwMhIuLS53Wm+zs7DqtPCbt27dHhw4d4O/vb14WFRUFRVFw+vRpdOnSpc5zjEYjjEajLaERkd5poC8wEREROQ6bWnTc3d3Rt29fJCQkWCxPSEjAwIED633OoEGDcPbsWVy8eNG8LCUlBQaDAR07dmxCyERERERERFdmc9e1GTNmYOnSpVi2bBmSk5Mxffp0pKenY8qUKQBEt7MJEyaY13/wwQfRpk0bPProo0hKSsKWLVswa9YsPPbYY/D09Gy+d0JE+sYxOkRERGQDm6eXHjduHHJzc/HKK68gMzMT0dHRWLt2LcLDwwEAmZmZSE9PN6/v4+ODhIQE/P3vf0dsbCzatGmDsWPHYt68ec33LohI9yTICD01U9xWylWOhpyJJEkIDQ013yb9YNkSWU+L9UVSFMc/TVpUVAR/f38UFhbCz89P7XCISA3HlgB7plbff9Dhv7qIiIjIDqzNDZo06xoREREREZEjs7nrGhGRGhQFKPIXFyv2K/wN2mg0Jz1QFMV84Wo/Pz/NdNmgxrFsiaynxfrCFh0i0gQFBmSEL0BG+AIokrva4ZATURQFGRkZyMjIgAZ6e5MNWLZE1tNifWGiQ0REREREusNEh4iIiIiIdIeJDhERERER6Q4THSLSBo30ByYiIiLHwESHiIiIiIh0h4kOEWmDBqaxJCIiIsfB6+gQkSZIUNAh45/itlKhcjTkTCRJQocOHcy3ST9YtkTW02J9YaJDRJogSQpa5/+gdhjkhCRJQuvWrdUOg+yAZUtkPS3WF3ZdIyIiIiIi3WGLDhFpgqIAF32HAAB8LuyANhrNSQ8URcHFixcBAD4+PprpskGNY9kSWU+L9YUtOkSkCQoMOBW5GKciF0OR3NUOh5yIoig4deoUTp06BYXTnOsKy5bIelqsL0x0iIiIiIhId5joEBERERGR7jDRISIiIiIi3WGiQ0REREREusNEh4iIiIiIdIeJDhFpg0ZmeCEiIiLHwOvoEJEmSFDQ/syr4rZSoXI05EwkSUL79u3Nt0k/WLZE1tNifWGiQ0SaIBmANrlfqR0GOSFJktCmTRu1wyA7YNkSWU+L9YVd14iIiIiISHfYokNEmqAoQLF3LADAu3g/tNFoTnqgKAqKi4sBAN7e3prpskGNY9kSWU+L9YUtOkSkCQoMSOu0HGmdlkORjGqHQ05EURSkpaUhLS0NCifF0BWWLZH1tFhfmOgQEREREZHuMNEh0ovSbCB5IVBwSO1IiIiIiFTHMTpEerFtHJC9CTgA4AEZ0EDfWSIiIiJ7YYsOkV5kb6q+rVSqFgYRERGRI2CiQ0REREREusNEh4iIiIiIdIdjdIhIIxQEZS64fLtC1UjI+QQFBakdAtkJy5bIelqrL0x0iEgTDJDR9vwKtcMgJ2QwGNC2bVu1wyA7YNkSWU+L9YVd14iIiIiISHfYokNEmqBAQolnDwCAZ0kyOHk2tRRFUVBSUgIA8PT0hMSp23WDZUtkPS3WF7boEJEmKDDgRJevcKLLV1Ako9rhkBNRFAUnTpzAiRMnoCiK2uFQM2LZEllPi/WFiQ4REREREekOEx0iIiIiItIdJjpERERERKQ7THSIiIiIiEh3mOgQ6ZFGBgkSERER2QsTHSLSBiZvREREZANeR4eINEJB23OLL9+uUDUScj5auxo4WY9lS2Q9rdUXJjpEpAkGAxB0bonaYZATMhgMCAoKUjsMsgOWLZH1tFhf2HWNSJfYzYuIiIicW5MSncWLFyMyMhIeHh7o27cvtm7datXztm/fDldXV/Tu3bspL0tETkyRFZQaO6HU2AkKJLXDISeiKApKS0tRWlqqmauBk3VYtkTW02J9sTnRWbVqFZ555hnMmTMHBw4cwJAhQzBq1Cikp6df8XmFhYWYMGECRowY0eRgich5KTDgeLc1ON5tDRTJQ+1wyIkoioLjx4/j+PHjmvlxJ+uwbImsp8X6YnOis3DhQkycOBGTJk1CVFQUFi1ahNDQUCxZcuW+80888QQefPBBxMXFNTlYInJiEltxiIiIyHo2JTrl5eXYt28f4uPjLZbHx8djx44dDT5v+fLlSE1Nxdy5c616nbKyMhQVFVn8ERERERERWcumRCcnJwdVVVV1ZlwICgpCVlZWvc85duwYXnjhBXzxxRdwdbVukrf58+fD39/f/BcaGmpLmETEyQiIiIjIyTVpMgKpVhcSRVHqLAOAqqoqPPjgg3j55ZfRtWtXq7c/e/ZsFBYWmv8yMjKaEiYRERERETkpm66jExgYCBcXlzqtN9nZ2fXOq33hwgXs3bsXBw4cwFNPPQUAkGUZiqLA1dUV69evx4033ljneUajEUaj0ZbQiIiIiIiIzGxq0XF3d0ffvn2RkJBgsTwhIQEDBw6ss76fnx/++usvJCYmmv+mTJmCbt26ITExEddff/3VRU9ERERERFQPm1p0AGDGjBkYP348YmNjERcXh48//hjp6emYMmUKANHt7MyZM/jss89gMBgQHR1t8fx27drBw8OjznIiakYamfbRNgoCzy+/fLtC1UjI+QQGBqodAtkJy5bIelqrLzYnOuPGjUNubi5eeeUVZGZmIjo6GmvXrkV4eDgAIDMzs9Fr6hAR2coAGcGZC9UOg5yQwWBAcHCw2mGQHbBsiaynxfoiKRq44k9RURH8/f1RWFgIPz8/tcMhckwra0wIMvYS4OqpXiz2kLIY2Ptk9f0HHf6ri4iIiOzA2tzA5hYdIiI1KJBQ4RYCAHCryAQvH0otRVEUVFSI7pJubm71zjJK2sSyJbKeFutLk6aXJiJHp7/WDgUGpET9ipSoX6FIHmqHQ05EURSkpKQgJSUFGugEQTZg2RJZT4v1hYkOEWmDRr5UiYiIyDEw0SEiIiIiIt1hokNE2qCBvsBERETkOJjoEOkSu3kRERGRc2OiQ0REREREusNEh4iIiIiIdIfX0SEijVAQkPPV5duVqkZCzicgIEDtEMhOWLZE1tNafWGiQ0SaYJAUhJx9Ve0wyAkZDAaEhISoHQbZAcuWyHparC/sukakR7zmDBERETk5tugQkSYosoIql9YAAJeqfHCyaWopiqKgqqoKAODi4gKJU53rBsuWyHparC9s0SEiTVBgwJEeW3CkxxYokqfa4ZATURQFR44cwZEjR6CwtVRXWLZE1tNifWGiQ0TaoIEzR0REROQ4mOgQ6ZI2zrQQERER2QsTHSLSBo00kxMREZFjYKJDRERERES6w0SHiLSBY3SIiIjIBkx0iHSJ3byIiIjIufE6OkSkEQpa5f1w+XalqpGQ82nVqpXaIZCdsGyJrKe1+sJEh4g0wSAp6Hj6n2qHQU7IYDCgY8eOaodBdsCyJbKeFusLu64REREREZHusEWHiDRBUQBF8gQASEoJODUBtRRFUcxXAZckCRInxtANli2R9bRYX5joEOmRDq85o8CApJjdAIDuf/VnokMtRlEUJCUlAQC6d++uiR93sg7Llsh6Wqwv7LpGRNqgw+SNiIiI7IeJDhERERER6Q4THSLSBg00kRMREZHjYKJDpEvs5kVERETOjYkOEWkDx+gQERGRDZjoEBERERGR7nB6aSLSBgnwK1h/+U6VqqGQ8/Hz81M7BLITli2R9bRWX5joEOmS/rp5GSQFYekz1Q6DnJDBYEBYWJjaYZAdsGyJrKfF+sKua0REREREpDtMdIhIGzgZAREREdmAXdeISBNkxYCknn8BALr/1Z9naajFyLKMpKQkAED37t1hMHDv0wuWLZH1tFhfHD9CIrKdHls/eMFQIiIisgETHSIiIiIi0h0mOkSkDXpspSIiIiK7YaJDRERERES6w0SHiLSBY3SIiIjIBkx0iHSJ3byIiIjIuXF6aSLSBkWGT9GWy3eqVA2FnI+Pj4/aIZCdsGyJrKe1+sJEh4g0wSApiEh7Uu0wyAkZDAZERESoHQbZAcuWyHparC/sukZE2sAxOkRERGQDJjpEesSpmImIiMjJNSnRWbx4MSIjI+Hh4YG+ffti69atDa77/fff46abbkLbtm3h5+eHuLg4/Prrr00OmIick6wYcDj6DxyO/gOy5Kl2OOREZFnG4cOHcfjwYciyrHY41IxYtkTW02J9sTnRWbVqFZ555hnMmTMHBw4cwJAhQzBq1Cikp6fXu/6WLVtw0003Ye3atdi3bx+GDx+OMWPG4MCBA1cdPBE5EUWBYvCCYvBSOxJyQoqiQGFLqS6xbImsp7X6YnOis3DhQkycOBGTJk1CVFQUFi1ahNDQUCxZsqTe9RctWoTnnnsO/fr1Q5cuXfDaa6+hS5cu+Omnn646eCIiIiIiovrYlOiUl5dj3759iI+Pt1geHx+PHTt2WLUNWZZx4cIFBAQENLhOWVkZioqKLP6IyBbaOdtiNU5GQERERDawKdHJyclBVVUVgoKCLJYHBQUhKyvLqm0sWLAAxcXFGDt2bIPrzJ8/H/7+/ua/0NBQW8IkIiIiIiIn16TJCKRaZ1YVRamzrD5ffvklXnrpJaxatQrt2rVrcL3Zs2ejsLDQ/JeRkdGUMIlITzTUJ5iIiIjUZ9MFQwMDA+Hi4lKn9SY7O7tOK09tq1atwsSJE/HNN99g5MiRV1zXaDTCaDTaEhoREREREZGZTS067u7u6Nu3LxISEiyWJyQkYODAgQ0+78svv8QjjzyClStX4rbbbmtapERkAx22fkgSvC7ugdfFPQC0Ma0l6YeXlxe8vDjjnx6xbImsp7X6YlOLDgDMmDED48ePR2xsLOLi4vDxxx8jPT0dU6ZMASC6nZ05cwafffYZAJHkTJgwAe+88w4GDBhgbg3y9PSEv79/M74VItIzgyTjmhOPqR0GOSGDwYBrrrlG7TDIDli2RNbTYn2xOdEZN24ccnNz8corryAzMxPR0dFYu3YtwsPDAQCZmZkW19T56KOPUFlZiSeffBJPPvmkefnDDz+MFStWXP07ICLnwDE6REREZANJ0cBVf4qKiuDv74/CwkL4+fmpHQ6RY1pZY0KQu88BHg1P+KFJKYuBvdUnS/Cgw391ERERkR1YmxvY3KJDRKQGWTHgaPfNAIBuybc0bcpIoiaQZRlHjx4FAHTr1g0GA/c+vWDZEllPi/WFiQ6RHjl+Q22TVLk2fKFhInuqqqpSOwSyE5YtkZUKkqC16uL4qRgREREREanjYhqw82Hg1+url5VkqxaOLdiiQ0Qaoc9WKiIiIodUcg44/Cpw/ENArgAkz+rHilMB72D1YrMSEx0iIiIiIhLKC4EjC4AjC4HK4url7q2qbwfGtXhYTcFEh0iX2PpBRERENqgqFTOcJr0GlOVWL3fxAq6dDnSbARw/q158TcBEh4iIiIjIWcmVwMn/AH+9BFw6Xb3c4AZ0fgLoMQfwDAZkGQATHSKi5qco8Lx06PIdWdVQyPl4eno2vhJpEsuWnJaiABnfA3/OAYqO1nhAAiIeAnq+DPhcY/EUrdUXJjpEpAkGSUan4w+oHQY5IYPBgE6dOqkdBtkBy5acVtZvQOILQN5ey+Uho4FerwKte9Z5ihbrCxMdIl3iGB0iIiKqJXcPkDgbOPeb5fK2g4HerwNtB6kTl50w0SEiIiIi0rPCI8Cf/wQyvrNc3qon0Gs+EDIKkCR1YrMjJjpEpAmyIuHYtesAAF2O3smrHVOLkWUZx44dAwB06dIFBgP3Pr1g2ZLuFWcAh14GTiwHlBrjW32uAXr+PyD8fkCybr/XYn1hokNEmlHh3uHyLf2ddSLHVlFRoXYIZCcsW9Kl0hwgaT6Q8gEgl1Uv9wgGYl4ErpkIuLjbvFmt1RcmOkR6pHCMDhERkdOpuAgc+TeQ/BZQeaF6uZs/0P15oNs0wNVbvfhaGBMdIiIiIiItqyoDjn8EHJoHlJ2vXu7iAXR7Goh6DjAGqBefSpjoEBERERFpkVwFpH0B/PUiUHyqernkAnSaBET/C/Dq0PDzdY6JDhFpg8KLhBIREQEQXdTP/AgcnAMUHrZ8LGycmGjAr4s6sTkQJjpEpA1KldoREBERqe/cZnGxz9xdlsvb3wL0eg0I6KNOXA6IiQ6RLulwMgKlAsbS46Y7qoZCzsdoNKodAtkJy5Y0I+8AcPAfQOY6y+VtBgC95wNBw+wegtbqCxMdItIEg1KOLil3qR0GOSGDwYAuXdgFRI9YtqQJF44Df/4LOPWV5XL/HkCvV4EOt7fIxT61WF+Y6BCRNiiVakdARETUci6dBQ69AqR+avkb6B0OxLwCRDwEGFzUi08DmOgQkTbITHSIiMgJlOcDSW8AR98FqkqqlxvbAtH/BDo/AbhoqwuZWpjoEOmS/sawyIqM1K6rAQCdjj0Ag8rxkPOQZRmpqakAgE6dOsFg4N6nFyxbciiVxSK5SXoDqCisXu7qC0TNAq59BnDzVS08LdYXJjpEpA1yFco8Ol++Y/++yEQ1lZWVqR0C2QnLllQnVwCpS4G/XgFKs6qXG4xA1yeB7rMBj0D14qtBa/WFiQ4RaQPH6BARkZ4oMnBqlZho4GJq9XLJAEQ+AsTMBbzDVAtPD5joEJE2cIwOERHpgaIAZ38RU0UXHLR8LPQecbFP/yh1YtMZJjpEeqTob4wOW3SIiEjzzm8HEmcD57daLg8aIS72Gdhfnbh0iokOOY+0r4DKIuCaxwADd33NYaJDRERaVfAXcHAOcOYny+UBseJin8Ej1YlL53i0R84h63dgxwPV9zs/rl4s1DTsukZERFpz8QTw51wg7QtYzIjq1w3o+SoQeneLXOzTWTHRIedw/MPq23+9zERHi+QKuJWfuXxHh13zyKG5ubmpHQLZCcuW7KLkHHB4HnD8IzGrmolXRyDmJSDyYU32LtFafdHeJ0x01exw5qSqHMjZDrQZALh6Nv/2CQalFN2O3KJ2GOSEDAYDunXrpnYYZAcsW2p25YVA8lvA0UXiujgm7gFAj38AXaZq9jhBi/WFiQ5Rc/hjEpD2ORByGzDsZ7WjgS5bPDhGh4iIHJFcBeTtAzJ/ERf8LM+rfszVG7h2BnDtTMDdX70YnRQTHXJCdkgC0j4X/8/+r/m3TQLH6BARkSNQFKDwMHDudyDrNyB7M1BRaLmOwQ3oPAXoMQfwDFInTmKiQ86CA/20TpYVnOz8JQAgMvURGFSOh5yHLMs4efIkACAyMhIGA/c+vWDZktUunhATG537TSQ4pdkNrCgBkePFOByfyJaM0O60WF+Y6BCRNigySryiL99x/C9X0peSkhK1Q3AeBYeA9G8A/+5A0HDAo51dX45lS/UqyRIJjanVpjit4XWNgUDQjUDwCCA4HvCJaKkoW5zW6gsTHXISOhyzckU6fL8co0Okf5fOAhtuAMrzq5e1ihEXU2wfD7S/hVPxkn2UFwDnNl1Obn4DCpMaXtfVF2g3VCQ2QTcCraIBiSfgHBETHXJCOvyRVHSY2JhUXgJcvSwTHdMPSnm+mNXGq6M6sRE1B0UWFxP07SwGLjsrRRETu9RMcgDx2RT8JWaxChsLDPqSB5V09SovAee3i6Qm63cgf5+oi/UxGIG2g6pbbQJiNTk1tDNiKRHpgV5bO/Y/CxxZAPT4Z933eOkM8HMUUFUC3LQdCOyvToxEV+vQq8BfLwL+PYBRic57AJW6VMxaBQAewWKcw7mNYjYrUyt1+teA37VAz5dVC9Os8hKQuU50awocCEQ+pHZEDSs+BeQfFPuYbye1o1GHXAHk7hblde53IGcnIJfXv65kAAL6XW6xGQEExml2Smhn56TfpqQZlSVA6ieAdyTQcYza0VhHUVq+a0XNi5HpyZEF4v/heUDr2OrligIcmAVUXhD3/3oJGL5W3C44BJz5WRwkeXVo0XCJbKbIwPEl4nbhYSA/EWgTe8Wn6NLFk8D+GdX3r18KdLhN3C7PB9K/BXY/AUABDr0iugqF3dfycVZcAM78D8j4Dji7Fqi6JJYfWwwUHQF6vuI4XevkCuD0j+I3NHM9zMmidyQQPBJof5NooTC2UTVMu1FkkdyZxtic32J5XZvaTF0kg24E2t3AqaB1gokOObak18WPGgDclgT4R6kbjzWUKkBq4apVO9HRYle2knPAzgmARxAwYDlgcLF8vOAwEGq6owCnvqx+zJTwAMC6vuIs3Zkfgfgd9o6a9EquqrsP2kPuHqAks/p+9mbnS3QUGdj1KFB5UdzvNLE6yQEA99ZA58li+t4Ds8SynQ8DPp2BgD72j6+8UHyfZHwHnF0HyGX1r3d4noix7yJ1u9YVHROtYydX1D8zWPFJkfykfgJAAgKuE4lP8Eig7WDAxaN54qgoEi1JF9PEQP7iNLGvSy7iNVw8RJewmv9dTPdr3rZyHckFuJBSndic22h5PZvafDpVd0VrgUkvSB1MdMixmZIcADi9pumJjj0P/OWqWq9VhRavWnUSnar613Nke58CstaL28EjgIi63UBcKhv40TK4i/9yVXVXhJyddgiSNEuRgQupottOYwehqcuBvU8CEQ+KlgUALi52SnpOr7G8n70FiJppn9dyVEffEQkeAHiHA9ctrH+9a2eKsTonPxNdVrfcAdy856qvUXLFss36Ddh6T91rpABipq2Od4kD5MOvimUp74l1r/+0ZbsgVpUCGauB4x8D2ZvqPu4dAXS4HSj4E8jZXuM3QxFdA/P2AUlviKSh7WAg+CaR+LTu3XB9qSiyTGJq3i4+deUkwx4kQ8NjbADRHdLUFS34RrGvkc3s9l1oJ0x0iK6W4gBJRsmZWjFoZMxO/kFxQBk0UpwtNcnbB4TeZbGqQSlBVNLQ+rdj+iGu70yrXAnk7Rdnfg1uzRQ42cXxj0U3pWtnNH9Z7fk/sf3Qe4Ah31553T8eE/9TPwV6/AMGn2sQFWWn1uTaic75reJg7UrJWN4BcSbbv7t1r1FVDpz6Smw7b6/oYuUfA3R/Thyot0TLVW0l54C8PaJFK+mN6uUDlgNufvU/R5KA/h8BRUeB3D+ASxnA1ruBEb+Lz6MJDAZDw2V7bjOweYxIqkw8goHQu4Gwe4G2Q6qTGd+uYr9RqkQiVlEEDPqqyXFZrTAJOP6JeM3aiYXBDehwh2gNCx5ZvU9VFgPZ24CsBCBrA1BwsPo5VaViWdYGcd/YRiQGAdeJ6ZZNSUxxWt1JI9RWO8lxayVaakytNn7XOk63Qo26Yn1xUEx0qGUUJovuXH5dmr6Nq+kKYOpHDTT/F51cK6lo6UQnawPw+02Wy36+Fui3RIxTqSx2jCZ5RRGxuPmI29vGAhmXDzjPb7dcN3szcGSR9dvO2gB84y8OLmrb9SiQ9l9xRtO3s+jX3/nx6sdTl4sxPt2fA7o+aeObogZdOiuS19C7rJsV79ymy2MwAJTlAX3euOLqNsdy/GNxO+M7cYDm3rrWOqeBw6+Jfvo1nf4BuHZ688VSU+ERkXTUVJ4vWi1a96r/Occ/rv6c2gwAuvyfGL/o3lp0sTq1UoxTC4gFQu8UY372PAkUJVtuJ28PsO0+wK+bGHQtV4gTJG2HiG26uDfPezy3UYxpKc8FynLFwXjxKfF519btaXFgeiUuHsANq4F1/cQJnpwdYlISrxBxYOvmD7i3Ajzbi7ru161pcZ/fAWy+rTrJaX+zmBQlMK7+xPCaCYCbL7D9ftGqfHoNsHk0cMOaxmfSKz4FnPoayPgeKMsR67v5AC7e4nZ9f5KLaMHJqad7rm9XoNMk4JqH6//ud/UGQm4Wf4BIOs/9fjnBSRAJpElZrpgAIv1raz41QXIBvELFd65PBOAVLv57R4jvAkURJ6WqSqv/W9wuA+TL/2s/VmedWreNbauTm9Z91EniyaFIiuL4nfmLiorg7++PwsJC+Pk1cKaHWk7KB8CJ5UCv+WIwY2Py9gHrYgFIwOhk2354VtZISnq/IQ5Gm2JtL9FkD4juBvecb9p26lOeD3wbUH3/3ry6B1H2UFUmzmyu69twC457gDj4v2krEDjA/jHVVp4vEpaA64Cj74oEZuAX4od2cwtMLvGgYrkPmdx3QRxIAJaPP6iIH/2U94GgYeIsINWlKACUuicfKovFWV/vcGBtjDiIDxwIxG+vdzMWEv8BJM2vvj/mePPNDnXoVeDPf1bfH/oT0GF09f2qcuDX/pZntk3aDQVGbmqeOGoqzxfv+fiH4r5PJ+Biqrjd+w1xsFjwl2jd8OkEdHpMXOfjf92rx7KYSAagVS/xeddsfZBc6343SAYxDq7muKDa/KKA/h+KAdmKApSdF2fLPYOtf3+KAhyaJ2aTs0b7W4Ah31s/s1XuXmDDEHFweyVBw4HOT4iWK2uTt5zdwMabqk+ctB8lkitrWmeyNgCb76g+uRYYBwz7X93fhJJMcVHUU181Tzdbg1G0MnWaLMqtqSf0FAW4cKw66Tm3sW63PckF8AqrTl68wy//jxDLPDs478yB1GKszQ24J5JQXii+1IKGX3kGFkURYykAYGO8ODBszO7/Mz0Z2D8TGPZzE4O8ipYYY9sam2nmMzy1x8dUNTBdZXPJ/xM49LI4+9cYU1eGrfcAt58UPz4XTwDHPhTdL9oOvLpYGhqwLVcBZ/8n9pWaZwcBccbTM8T2l5KMSIsUs1NFnPw/GJQGBgTX1NDB3IWjQOvr6k4teums+GyPfywGFkc+DPSeL8Z2+FwjzhoXHhFnD1v1FPWh5vsvyRQteqYWDEURB4gGF6DiojjLGxh39QfwchWQ/KY4uLl2upiM4dIZsU8Ej7BMaqtKRbcWvyig3ZB6tlUpDoZrDkDO2SUOvjo/bnk2Wq4A0r4QiUPVJaDfh6LrZuk5IPQ+4NdYcXa6/ajqloqcHWL9wOsB326iO1jAdUCHMWKykeytYkxGfqJlXCc/a3wKYblKnIn2DgMSnxcX8Ws7WEz5695ajPMKGi5es6ZTq8QZf1cfYO+06rFh9cneDHnPNKT5TQTcfBHhkQqDpIhpehub1a8sT5RT9hYxo5p7a3G2veSMOElRs/W39xvAtnvF7YOz63bDSV0qzs6bkhz31tVdhxQZyD9Q9/VrJjltrgd6vSr2DRcv8X1/+LXqsTE1FSUDG4YC/tHApfQaB/y3AH3eEjOeAeIkVvLbYn8JHwd0f17EJVcB+/4OHFtS/+fi1kq0WLXpJ1qT2vQTB8m2HJy3iRVdw/Y9I/a/mgleTec2ij+PdsA1j4p92uca88OyLCMtLQ0AEBERAUNBIrDx5ur3HDwSuOF767ugBY8EbtwAbLoVqCgQ9WjDMGD4+sutMN8B6atEC2Z9F3Z2ayXqVkPTHtfmHy26pkX8DTAGNL5+YyQJ8Osq/rpOvdz1d6+YFc+rgygnzxAmMk6qTn0xOP71rJrUorN48WK89dZbyMzMRI8ePbBo0SIMGVLPD+hlmzdvxowZM3D48GGEhITgueeew5QpU6x+Pd236OTsAk5+Lr6Aa3dXOLYE2DMV6LMAiJpR97mpy8TZ56E/1f+jW5YHrA4RAwrjdzb8Q7LtfvHl23awOPvfkMoS4Guv6vvWJDpre1efKQ0aLvpTW8uiRedNoPss659b0283ih87QByg3lrPmdumunQaWBNafX/MMdFFqqlqT099fofYP/L2iR/K3F1N3zYgEtmyXHG7dR9xgOJ3rTgw9e0qznRXlYiDxvSvxcF0h9HiYCrje6BVD+DMT8DBOdXb9AoV3TYqLtRNbJqJLHkiKWY3AKD7X/1hUBo4sLFF3/fEAZk1DEYxSPzwa5bL2w0D+r4D/DVXJDKSC3DzH+Iz/eU6MQuQX5Q4ADR1Mxl9VBzoVl4U+0/nyUDJWfH5G1zF9TmOLhIHPWFjxcFR8MjLswjli4PPlPfFtjza1Z1ZKew+IO6/4uB6yx3VZXLdQlGOYWMB71CR8O6fWX32uc/b4mzu8Y8ub2ecOLBVqsTBzaZbxcX1moObX/1dDU0MbsDILeJ1q0rEwWlVCZCZILqTyaVijNeFlCu/Ts2EoIka3Pc63inK3thGHHCnfy3qVPBNohXm4D/EDFeNCb1XHLR/16b+Qe+1ebYXs1AWHgbSvxOJWuFh0W0r4m9AyC1iEH3GdyI57fmKmMmsvu6/l06L73UXd5Ho739GXGukIZJBnAAoTq+7L7j5i2Qnb6/liZie84COd4jPyb1N83WLq6mqXHx2FYWi5ev8VrF/17d/tB0MhN8PhN0H2T0QSUlJAIDu7WUYNg6vPkHUbphojXH1qruNxuQfFCcDTXXTGCj2w/q6NreKuRzPuOqTIHKlaCE1/VUVi5MlptuVly53O4zluBNqMbIsV9eX7t1VTXSszQ1sTnRWrVqF8ePHY/HixRg0aBA++ugjLF26FElJSQgLC6uz/smTJxEdHY3JkyfjiSeewPbt2zF16lR8+eWXuOeee5r1zWiW6WDe4A7cX1b/Y0DdpKI0G/g+qOHHAWBNuDgoAsQPcrdpV46hoe2YlOUC3wVaty4gDtrXXVd9tjYwzvopfyuLga99qu+3HQLctKWR15PFD0nREfFDXHQEOPBs3fU63ikSnm7TxAG9wU08T66w7DphmlFNMohkTZHFdQmChomuWKdXN/4+2g0VB2i5u8XryBXiILjL/4kDtnM2JH5OzC6JjiNqFSO6LF2t1tcB+fvrf8ytFdD7dWCP9SecnJMEQLnyvicZRBfRspwrb8o7XJx4qrwgvgf8rhWDvMPHidYWSQI2jRYtoYBICuL+K1oS9vyfSIxNBn8LhNX6/awoAlw8607iYOt1veQqIPVjMW6tLEeMr/DtIr5LTb8l1pJcgQHLxFhBNSiKmIHs2EfA6e/rtr5LBshBo5DU7nUAQPeUW2EovXxSoO0gYNi66i6uTVGUAvw+sv6TP75dqpObVj2a/hpELcgpEp3rr78e1113HZYsqW6SjoqKwp133on58+fXWf/555/Hjz/+iOTk6oGQU6ZMwcGDB7Fzp3X9Up0m0dE6r471DzAlagZOk+g4OoM7cP0yIPNXIO1ztaOpn2eIaFnY93T1Mr8oMXbB1FplEjRCTNxhcBUnHf6YJJZHThAtYzl/QE5aYLnvhY4RXb5Kz1luq77pbdv0BwZ/LRIdRREnqIwB9c8qd2YtsOV2Mf5h2M/Vs6qVFwDbHwQyfxFxDVhh/7P4iiK6vpnirCwBUt4VLZqmljifTkDULNFSn/SGuG6L6f27eouELOQW+8ZprdJs0QPi5GcWEzPU+73S5nrgxvUNz/5mi+J0YOMt4jW9wkRyE37/5WmbdfLbT05Di4mOTZ0sy8vLsW/fPrzwwgsWy+Pj47FjR/1n6Xfu3In4+HiLZTfffDM+/fRTVFRUwM2t7pd9WVkZysqqWzaKiq7QvYEcB5MculrdngEuHAc6TwK23Nnweq2igfw9lssGrhRd7LaNFd1XeB2dppNcRV/8i8frPtb7DSDyISB8rOjSk5UgWqFO/1C9TpsBYkxOj3+Is9rFJ8XFYA3uYqzCnidEC4VfN8D7GtHyVFkM+EQCPf+fmDK4oethtBsqumelfVH3eiEel1u4b/hBjOEoy62+FlePf4iz9Lm7xetcv0y8pmf76gNOrwnigpAXjor1/bqJfSo4HjDNXxJ6FzDov2IfS3oDOPMzUHhIHLgO+kp0XczZKbouGduIMVSm8U+SdOVrvnS4FbgnR4wdqjkGwr0VMHytmCjDo13LHCBLEiDV+H129RTJ4zWPia60PpHiuiymMWoDPgWinhWJUMkZsZ+06Wf/OK3l0Q7o8YJ4DwV/iUkATn0JXKqVrLa+Dhi+rnmSHECMHxuVKFp1fK5hckPUwmxKdHJyclBVVYWgIMsv6qCgIGRlZdX7nKysrHrXr6ysRE5ODtq3b1/nOfPnz8fLLzcyCJWIxDVBio6KAy1IYla645+IPuYht4mDjbXRzfd6nu3rDvD36Wx5QBw8UszY4xUqxiqUZomB0q6+1f3lvcPFD7+LlzigdvcXF9irOSPfg4qYOOHiSeDIQqDTE4Bp+MKN64GDL4hJJnpcPvFiGjQ//Bfxvzgd+PNFcXDd4wVxkHHqa3FQLbmIA1FTN6F2NwA9XxUH7Oe3i/suHqLLpc814sAz5T0g+S1xIJS1XnR1DL9fjJsxthVn7ovTxCDi8gIxDqDnK+K1KgqBXY+Jg5yw+8R06x1Gi88z9VPxOfT4h+gqlP61eO1Ok8RrnPlJdHOqKBKffVEycM0jor//7ikAJHFdmPyDQLvBlweDV4pE4MzPQMitYprZrN/ELE+h9wBl2aJFIv1b8Z5v/E3MrHXhuDgwc28tyit7kzjgy90tXqdVtPhsAHGmv/+H1eVVcFh8th3GiANzE492AAaLMncPENsYcXmAeIcxYmxXbTfvFZ9ZxmoxHiH0XjGpQofR1QlAp4li//AOF+MxvMMtBpkDAGLmigNyGMTEBJIEjKrRnc/d33J9gxsw5Ju68bQdBJwXZzFx/VKxHfdWYqKK3vMvj3ExVo+B8e3c9O5aNT+72q7ywpjNwqNt/eNFAXFB54EO2spnIklA657ir9erwPndgGmIW/BIYOCKK5dBU7i4N98MgkRkE5u6rp09exYdOnTAjh07EBcXZ17+6quv4vPPP8eRI0fqPKdr16549NFHMXv2bPOy7du3Y/DgwcjMzERwcN3pKutr0QkNDdVv1zVA/FCW54mDJIObOPCRXCCmcHURBzwGd9GVoLxQ/NhXlV0+IPECSs+LH+3KkstdDWQxeFqpEgeAlcWXZ3ERfc7NywExwNHdX5y9k0svv5aH+EGoLBavKxmq57539RZnfOXy6rOuBtfq/s/mM7GXdy2lUhw0SQbx+nKZeK2KQhG75CIOYlw8qgdqSi6X378k4jBtU5Evb9cg+rq7eIn/rj7Vz5HLL39mlSJOU3xVZSIGyVVsQ66onoGt5oXUJFex3OB6+bNwuzz2R7Fc1xRP7VnHTH3iTTNuKZXiMzQxndFTZJhnkmvsLJ8iX34t18YvJtiYmjOBXY2rjcNGjtRkrhtyJWdPsgL3Pf1i2RJZz5Hqi126rgUGBsLFxaVO6012dnadVhuT4ODgetd3dXVFmzb1T2NsNBphNNr5asKOxtUTcG1kqlKT+prUvS5P19vQhcnqO2tqfu0azzHUen6dM1s1BurXnm7T2iuZm7px1JwK0/zcBrZhOqCueWBtis2l1n5kjsu1geUAINU/XWjtz6nm69aXi9Q3VbVUI3mRXAA0kFDYkiRIhvo/g6Ywx3WVWjDJMb8ku300LyY5VuO+p18sWyLraa2+2PQr5+7ujr59+yIhIQF33XWXeXlCQgLuuOOOep8TFxeHn376yWLZ+vXrERsbW+/4HCKi+hgMBvTowdmJqOVx39Mvli2R9bRYX2w+JTtjxgwsXboUy5YtQ3JyMqZPn4709HTzdXFmz56NCRMmmNefMmUKTp06hRkzZiA5ORnLli3Dp59+imefrWfKXyIiIiIiomZgc7+FcePGITc3F6+88goyMzMRHR2NtWvXIjw8HACQmZmJ9PTqufYjIyOxdu1aTJ8+HR988AFCQkLw7rvvWn0NHSIiIiIiIlvZfB0dNej+OjpE1ChZls0nUcLCwjhomFoM9z39YtkSWc+R6otdJiMgIlLTxYsX1Q6BnBT3Pf1i2RJZT2v1hacuiIiIiIhId5joEBERERGR7jDRISIiIiIi3WGiQ0REREREusNEh4iIiIiIdEcTs66ZZsAuKipSORIiUossy+bZXoqKijgNLLUY7nv6xbIlsp4j1RdTTtDYVXI0kejk5uYCAEJDQ1WOhIiIiIiIHMGFCxfg7+/f4OOaSHQCAgIAAOnp6Vd8M1rWr18/7NmzR9cxNNf2r2Y7TXmuLc+xZt3G1ikqKkJoaCgyMjJ0e4Fc7u8tsx177u/NtR73d33EoMX9vbm/261Zj/u7PmLg/t74evbe1xVFwYULFxASEnLF9TSR6Jiaxvz9/XX7xeDi4qL6e7N3DM21/avZTlOea8tzrFnX2u35+fmpvk/YC/f3ltmOPff35l6P+7u2Y9Di/t7c3+22rMf9XdsxcH+3fj177uvWNH6wM6qDePLJJ9UOwe4xNNf2r2Y7TXmuLc+xZl1HKGu1OcJnwP396p7T3OvpmSN8Btzfr2597u/Wc4TPgPv71a2vp/1dUhobxeMAioqK4O/vj8LCQtXPEhDZG/d3cibc38mZcH8nZ+Eo+7omWnSMRiPmzp0Lo9GodihEdsf9nZwJ93dyJtzfyVk4yr6uiRYdIiIiIiIiW2iiRYeIiIiIiMgWTHSIiIiIiEh3mOgQEREREZHuMNEhIiIiIiLdcehEZ8uWLRgzZgxCQkIgSRLWrFmjdkhEzaKxfVtRFLz00ksICQmBp6cnhg0bhsOHD6sTLJGNmmP/Lisrw9///ncEBgbC29sbt99+O06fPt2C74Korpbat/Pz8zF+/Hj4+/vD398f48ePR0FBgZ3fHTk7R9q/09PTMWbMGHh7eyMwMBDTpk1DeXm5ze/JoROd4uJi9OrVC++//77aoRA1q8b27TfffBMLFy7E+++/jz179iA4OBg33XQTLly40MKREtmuOfbvZ555BqtXr8ZXX32Fbdu24eLFixg9ejSqqqpa6m0Q1dFS+/aDDz6IxMRErFu3DuvWrUNiYiLGjx9v9/dHzs1R9u+qqircdtttKC4uxrZt2/DVV1/hu+++w8yZM21/U4pGAFBWr16tdhhEza72vi3LshIcHKy8/vrr5mWlpaWKv7+/8uGHH6oQIVHTNWX/LigoUNzc3JSvvvrKvM6ZM2cUg8GgrFu3rsViJ7oSe+3bSUlJCgBl165d5nV27typAFCOHDli53dFJKi5f69du1YxGAzKmTNnzOt8+eWXitFoVAoLC216Hw7dokPkjE6ePImsrCzEx8eblxmNRgwdOhQ7duxQMTKiq2fN/r1v3z5UVFRYrBMSEoLo6GjWAXJYzbVv79y5E/7+/rj++uvN6wwYMAD+/v7c/0k1Lbl/79y5E9HR0QgJCTGvc/PNN6OsrAz79u2zKW4mOkQOJisrCwAQFBRksTwoKMj8GJFWWbN/Z2Vlwd3dHa1bt25wHSJH01z7dlZWFtq1a1dn++3ateP+T6ppyf07Kyurzuu0bt0a7u7uNtcBJjpEDkqSJIv7iqLUWUakVU3Zv1kHSAuaY9+ub33u/+QIWmr/bq46wESHyMEEBwcDQJ2zFtnZ2XXOcBBpjTX7d3BwMMrLy5Gfn9/gOkSOprn27eDgYJw7d67O9s+fP8/9n1TTkvt3cHBwndfJz89HRUWFzXWAiQ6Rg4mMjERwcDASEhLMy8rLy7F582YMHDhQxciIrp41+3ffvn3h5uZmsU5mZiYOHTrEOkAOq7n27bi4OBQWFmL37t3mdf744w8UFhZy/yfVtOT+HRcXh0OHDiEzM9O8zvr162E0GtG3b1+b4na1/a22nIsXL+L48ePm+ydPnkRiYiICAgIQFhamYmREV6exffuZZ57Ba6+9hi5duqBLly547bXX4OXlhQcffFDFqImsc7X7t7+/PyZOnIiZM2eiTZs2CAgIwLPPPouYmBiMHDlSrbdF1CL7dlRUFG655RZMnjwZH330EQDg8ccfx+jRo9GtW7eWf9PkNBxl/46Pj0f37t0xfvx4vPXWW8jLy8Ozzz6LyZMnw8/Pz7Y3ZdMcbS1s48aNCoA6fw8//LDaoRFdlcb2bVmWlblz5yrBwcGK0WhUbrjhBuWvv/5SN2giKzXH/l1SUqI89dRTSkBAgOLp6amMHj1aSU9PV+HdEFVrqX07NzdXeeihhxRfX1/F19dXeeihh5T8/PwWepfkrBxp/z516pRy2223KZ6enkpAQIDy1FNPKaWlpTa/J0lRFMW21IiIiIiIiMixcYwOERERERHpDhMdIiIiIiLSHSY6RERERESkO0x0iIiIiIhId5joEBERERGR7jDRISIiIiIi3WGiQ0REREREusNEh4hIxzZt2gRJklBQUAAAWLFiBVq1aqVqTFerud6DJElYs2bNFdfJzc1Fu3btkJaWBqDu5+lIfv75Z/Tp0weyLKsdChGRQ2CiQ0TkRMaNG4eUlBS1w7gqLfke5s+fjzFjxiAiIqJFXu9qjB49GpIkYeXKlWqHQkTkEJjoEBE5EU9PT7Rr107tMJqsoqKixd5DSUkJPv30U0yaNMnur9UYRVFQWVnZ6HqPPvoo3nvvvRaIiIjI8THRISJSmSzLeOONN9C5c2cYjUaEhYXh1VdfBQDceOONeOqppyzWz83NhdFoxO+//w4AKCsrw3PPPYfQ0FAYjUZ06dIFn376ab2vVbvb10svvYTevXvj888/R0REBPz9/XH//ffjwoUL5nUuXLiAhx56CN7e3mjfvj3+/e9/Y9iwYXjmmWcafE+m7X700UcIDQ2Fl5cX7rvvvjpdvpYvX46oqCh4eHjg2muvxeLFi82PpaWlQZIkfP311xg2bBg8PDzw3//+t96ua0uWLEGnTp3g7u6Obt264fPPP7d4/NixY7jhhhvg4eGB7t27IyEhocHYTX755Re4uroiLi6uzmP79u1DbGwsvLy8MHDgQBw9etTqeEzvKzEx0bysoKAAkiRh06ZNAKq7yP3666+IjY2F0WjE1q1bcfDgQQwfPhy+vr7w8/ND3759sXfvXvN2br/9duzevRsnTpxo9P0REekdEx0iIpXNnj0bb7zxBv71r38hKSkJK1euRFBQEABg0qRJWLlyJcrKyszrf/HFFwgJCcHw4cMBABMmTMBXX32Fd999F8nJyfjwww/h4+Nj9eunpqZizZo1+Pnnn/Hzzz9j8+bNeP31182Pz5gxA9u3b8ePP/6IhIQEbN26Ffv37290u8ePH8fXX3+Nn376CevWrUNiYiKefPJJ8+OffPIJ5syZg1dffRXJycl47bXX8K9//Qv/+c9/LLbz/PPPY9q0aUhOTsbNN99c53VWr16Np59+GjNnzsShQ4fwxBNP4NFHH8XGjRsBiETy7rvvhouLC3bt2oUPP/wQzz//fKPxb9myBbGxsfU+NmfOHCxYsAB79+6Fq6srHnvsMavjscVzzz2H+fPnIzk5GT179sRDDz2Ejh07Ys+ePdi3bx9eeOEFuLm5mdcPDw9Hu3btsHXrVptfi4hIdxQiIlJNUVGRYjQalU8++aTex0tLS5WAgABl1apV5mW9e/dWXnrpJUVRFOXo0aMKACUhIaHe52/cuFEBoOTn5yuKoijLly9X/P39zY/PnTtX8fLyUoqKiszLZs2apVx//fXm+Nzc3JRvvvnG/HhBQYHi5eWlPP300w2+r7lz5youLi5KRkaGedkvv/yiGAwGJTMzU1EURQkNDVVWrlxp8bz/9//+nxIXF6coiqKcPHlSAaAsWrTIYp3a72HgwIHK5MmTLda57777lFtvvVVRFEX59ddf640FgLJ69eoG38Mdd9yhPPbYYxbLTJ/nhg0bzMv+97//KQCUkpISq+Ixva8DBw6YH8/Pz1cAKBs3brR4nTVr1lhsx9fXV1mxYkWDMSuKovTp08e8fxAROTO26BARqSg5ORllZWUYMWJEvY8bjUb87W9/w7JlywAAiYmJOHjwIB555BHzfRcXFwwdOrTJMURERMDX19d8v3379sjOzgYAnDhxAhUVFejfv7/5cX9/f3Tr1q3R7YaFhaFjx47m+3FxcZBlGUePHsX58+eRkZGBiRMnwsfHx/w3b948pKamWmynoVYVk+TkZAwaNMhi2aBBg5CcnGx+vL5YGlNSUgIPD496H+vZs6f5dvv27QHA/Jk1Fo8tar/3GTNmYNKkSRg5ciRef/31Op8VIMZhXbp0yebXIiLSGyY6REQq8vT0bHSdSZMmISEhAadPn8ayZcswYsQIhIeHW/38xtTs+gSIaZdNUxQrimJeVpNpuS1M26i5/U8++QSJiYnmv0OHDmHXrl0Wz/P29rZ62zXjMy2rL9ba69cnMDAQ+fn59T5W8zMzbavmtM5XisdgMNSJq6Kiot7Xqf3eX3rpJRw+fBi33XYbfv/9d3Tv3h2rV6+2WCcvLw9t27a94nsjInIGTHSIiFTUpUsXeHp64rfffmtwnZiYGMTGxuKTTz7BypUrLcaDxMTEQJZlbN682S7xderUCW5ubti9e7d5WVFREY4dO9boc9PT03H27Fnz/Z07d8JgMKBr164ICgpChw4dcOLECXTu3NniLzIy0qYYo6KisG3bNotlO3bsQFRUFACge/fu9cbSmD59+iApKcmmWKyJx5SEZGZmmh+vOTFBY7p27Yrp06dj/fr1uPvuu7F8+XLzY6WlpUhNTUWfPn1sjpuISG9c1Q6AiMiZeXh44Pnnn8dzzz0Hd3d3DBo0COfPn8fhw4cxceJE83qTJk3CU089BS8vL9x1113m5REREXj44Yfx2GOP4d1330WvXr1w6tQpZGdnY+zYsVcdn6+vLx5++GHMmjULAQEBaNeuHebOnQuDwdBoq4iHhwcefvhhvP322ygqKsK0adMwduxYBAcHAxCtE9OmTYOfnx9GjRqFsrIy7N27F/n5+ZgxY4bVMc6aNQtjx47FddddhxEjRuCnn37C999/jw0bNgAARo4ciW7dumHChAlYsGABioqKMGfOnEa3e/PNN2P27NnIz89H69atmy0eT09PDBgwAK+//joiIiKQk5ODf/7zn41ut6SkBLNmzcK9996LyMhInD59Gnv27ME999xjXmfXrl0wGo1Wdc0jItI7tugQEansX//6F2bOnIkXX3wRUVFRGDdunHm8h8kDDzwAV1dXPPjgg3XGjSxZsgT33nsvpk6dimuvvRaTJ09GcXFxs8W3cOFCxMXFYfTo0Rg5ciQGDRpknhL6Sjp37oy7774bt956K+Lj4xEdHW0xffSkSZOwdOlSrFixAjExMRg6dChWrFhhc4vOnXfeiXfeeQdvvfUWevTogY8++gjLly/HsGHDAIiuYqtXr0ZZWRn69++PSZMmmafvvhJTS9rXX3/drPEAwLJly1BRUYHY2Fg8/fTTmDdvXqPbdXFxQW5uLiZMmICuXbti7NixGDVqFF5++WXzOl9++SUeeugheHl52RQzEZEeSUpTOloTEVGLysjIQEREBPbs2YPrrrtO1ViKi4vRoUMHLFiwwKLVqaaXXnoJa9assalLliNau3Ytnn32WRw6dMg8tsZRnT9/Htdeey327t1rc7JIRKRH7LpGROTAKioqkJmZiRdeeAEDBgxQJck5cOAAjhw5gv79+6OwsBCvvPIKAOCOO+5o8Vha2q233opjx47hzJkzCA0NVTucKzp58iQWL17MJIeI6DImOkREDmz79u0YPnw4unbtim+//Va1ON5++20cPXoU7u7u6Nu3L7Zu3YrAwEDV4mlJTz/9tNohWKV///4W04ATETk7dl0jIiIiIiLdcewOx0RERERERE3ARIeIiIiIiHSHiQ4REREREekOEx0iIiIiItIdJjpERERERKQ7THSIiIiIiEh3mOgQEREREZHuMNEhIiIiIiLdYaJDRERERES68/8BZhfwcvNhZ9YAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(10, 4))\n", "gs = gridspec.GridSpec(1, 1)\n", "ax0 = plt.subplot(gs[0,0])\n", "ax0.semilogx(period[1:n//2],\n", " abs(y_fft[1:n//2])/np.max(abs(y_fft[1:n//2])), \n", " color='orange',\n", " linewidth=2, label='solar PV') \n", "ax0.legend(loc='best', shadow=True, fancybox=True)\n", "ax0.set_xticks([1, 10, 100, 1000, 10000])\n", "ax0.set_xticklabels(['1', '10', '100', '1000', '10000'])\n", "ax0.set_xlabel('cycling period (hours)')\n", "\n", "#We add lines indicating day, week, month \n", "plt.axvline(x=24, color='lightgrey', linestyle='--')\n", "plt.axvline(x=24*7, color='lightgrey', linestyle='--')\n", "plt.axvline(x=24*30, color='lightgrey', linestyle='--')\n", "plt.axvline(x=8760, color='lightgrey', linestyle='--') \n", "plt.text(26, 0.95, 'day', horizontalalignment='left', color='dimgrey', fontsize=14)\n", "plt.text(24*7+20, 0.95, 'week', horizontalalignment='left', color='dimgrey', fontsize=14)\n", "plt.text(24*30+20, 0.95, 'month', horizontalalignment='left', color='dimgrey', fontsize=14)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We calculate the Fast Fourier Transform (FFT) of the onshore wind capacity factor time series. " ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "C:\\Users\\34620\\.conda\\envs\\env_pv_textbook\\lib\\site-packages\\ipykernel_launcher.py:7: RuntimeWarning: divide by zero encountered in double_scalars\n", " import sys\n" ] } ], "source": [ "t_sampling=1 # sampling rate, 1 data per hour\n", "x = np.arange(1,8761, t_sampling) \n", "y = CF['onwind']\n", "n = len(x)\n", "y_fft=np.fft.fft(y)/n #n for normalization \n", "frq=np.arange(0,1/t_sampling,1/(t_sampling*n)) \n", "period=np.array([1/f for f in frq]) " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We plot the power spectra for onshore wind CF time series." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(740, 0.95, 'month')" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAzoAAAFzCAYAAAD2eXw5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAA9hAAAPYQGoP6dpAACUY0lEQVR4nO3dd3wUZf4H8M9uekISCCUQCARpho4BFFCKBcXeTk89sAAnZ0FAUTl/HsrZFcVyIGcB9eyn2A5LVECqEiAIJLRQQgmEmgLpO78/Hmdndna219n9vF+vvJhsZmefzc6G+ez3KSZJkiQQERERERFFEHOoG0BERERERORvDDpERERERBRxGHSIiIiIiCjiMOgQEREREVHEYdAhIiIiIqKIw6BDREREREQRh0GHiIiIiIgiDoMOERERERFFnNhQN8AdFosFBw8eRGpqKkwmU6ibQ0REREREISJJEqqqqpCVlQWz2XHdxhBB5+DBg8jOzg51M4iIiIiIKEzs27cPHTp0cPhzQwSd1NRUAOLJpKWlhbg1REREREQUKpWVlcjOzrZmBEcMEXTk7mppaWkMOkRRymKxoKSkBADQpUsXp6VqIn/iuRe5+NoSuS8c3y+uhrQYIugQEQFAXV1dqJtAUYrnXuTia0vkPqO9X0IfxYiIiIiIiPyMQYeIiIiIiCJOxHRdkyQJDQ0NaGxsDHVTyABiY2MRFxfH6cqJiIiIIlREBJ26ujrs2bMH1dXVoW4KGUizZs2Qk5ODhISEUDeFiIiIiPzM8EHHYrGgqKgIsbGx6Ny5MxISEvgpPTklSRLq6upw4MABFBUVoU+fPoiNNfxbgYiIiIhUDH91V1tbC4vFgs6dO6NZs2ahbg4ZREpKCuLj47Ft2zYsX74cw4cPR0xMTKibRS7ExcWFugkUpXjuRS6+tkTuM9r7xePJCH755RdcccUVyMrKgslkwhdffOHyPsuWLUNeXh4SExNxxhln4PXXX/emrU6Fw1zeZCzyOVNUVIRffvklxK0hV8xmM3r06IEePXrw/U5BxXMvcvG1JXKfEd8vHrfy1KlT6NevH1577TW39t+9ezcuvfRSnHfeediwYQP+/ve/Y/Lkyfjss888bixRICQlJWHr1q1oaGgIdVOIiIiIyE887ro2ZswYjBkzxu39X3/9dXTs2BFz5swBAOTm5qKgoAAvvPACrrvuOk8fnsjv4uLiUF1djdraWsOVZImIiIiC4emngSVLgEGDgHvuAdq1C3WLXAt43Wn16tUYPXq0zW0XX3wxCgoKHH6CXldXh8rKSpsv8szChQvRvHnzUDfDK/5q+8iRIzFlyhSfj0PhwWKxoKSkBCUlJbBYLF4d49ixY7jrrrvw7rvv+rl1FMn8ce6Fk9WrV+Ouu+7C6tWrQ92UkIu011brm2++wV133YXt27eHuikUAX75xYI77ihBnz7Geb8EPOgcOnQImZmZNrdlZmaisbERR48e1b3P008/jfT0dOtXdnZ2oJtJYeTGG2/kH2XSVVNTg5qamlA3g6IQz73IZeTXlh/eULBYLMD69UDv3jXo3bvGENUcIAhBB4DddM+SJOneLpsxYwYqKiqsX/v27Qt4G8k9wRjHkpSUhDZt2gT8cYiIiIjIte3bgaqqULfCcwEPOm3btsWhQ4dsbisvL0dsbCxatmype5+EhASkpaXZfEWauro6TJ48GW3atEFiYiLOPfdcrF271vrzpUuXwmQy4aeffsLAgQORnJyMoUOHYtu2bdZ9Nm7ciFGjRiE1NRVpaWnIy8tDQUGBzeN8//33yM3NRbNmzXDJJZegrKzM+jOLxYJZs2ahQ4cOSEhIQP/+/fHdd99Zf75nzx6YTCZ88sknGDlyJBITE/Gf//wHALBgwQLk5uYiMTERZ555JubOnevwuX799ddo3ry5tcxZWFgIk8mE6dOnW/e58847cdNNNwGw77r22GOPoX///njvvfeQk5OD9PR0/PnPf0aV6h136tQpjBs3Ds2aNUO7du0we/Zst14HIiIiInLut99C3QLvBHwdnSFDhuDrr7+2ue2HH37AwIEDAzbwe+BAQJOtgqJtW0CTMxx68MEH8dlnn+Gdd95Bp06d8Nxzz+Hiiy/Gzp07kZGRYd3vkUcewezZs9G6dWtMmjQJd9xxB1auXAkAuOWWWzBgwADMmzcPMTExKCwstPmdnj59Gi+88ALee+89mM1m/OUvf8EDDzyA999/HwDw8ssvY/bs2Zg/fz4GDBiAt99+G1deeSW2bNmCbt26WY/z0EMPYfbs2ViwYAESEhLwxhtvYObMmXjttdcwYMAAbNiwARMnTkRKSgpuvfVWu+c6fPhwVFVVYcOGDcjLy8OyZcvQqlUrLFu2zLrP0qVLMXXqVIe/r5KSEnzxxRf45ptvcOLECdxwww145pln8OSTTwIApk+fjiVLlmDRokVo27Yt/v73v2PdunXo37+/ey8IRRSLxYIff/wRK1euxIkTJ9CiRQsMGTIEAwcOtNt327Zt+O2331BSUoKKigoAonvtueeei3PPPde6X21tLWbMmIGMjAw8+uijuo/5yCOPoKmpCU899RQXoY1ip06dwoMPPog+ffpg0qRJ1tv37NmD5557DgDw5JNPokWLFtafPffccygrK8MLL7xgXdNr48aNWLp0KUpLS9HQ0IDWrVvjnHPOwQUXXKA7taun+2sdP34cr776Ko4fP44JEyagT58+vv4qDGv79u2YM2cOLr30UuTm5uKrr75CaWkp4uLiMGDAAFx33XWIj4/Hli1b8O2332L//v1ISEjAOeecg6uuusrm993U1ISlS5dizZo1KC8vR0xMDLKzs3HhhRfa/Y5Xr16N9957D2PHjkWLFi3wv//9D/v27UNcXBx69+6N6667zrpmoLwvAKxZswZr1qyxHmfKlCno3r27zbHXrVuHH374AYcOHUJSUhLOOussXH311YiPjw/Ur5EiSNQEnerqauzcudP6/e7du1FYWIiMjAx07NgRM2bMwIEDB6z9RSdNmoTXXnsN06ZNw8SJE7F69Wq89dZb+PDDD/33LDQOHQIOHAjY4X126tQpzJs3DwsXLrTOYPfGG28gPz8fb731lk2l48knn8SIESMAAA8//DAuu+wy1NbWIjExEaWlpZg+fTrOPPNMALAJJ4DoZvb666+jS5cuAIB77rkHs2bNsv78hRdewEMPPYQ///nPAIBnn30WS5YswZw5c/Cvf/3Lut+UKVNw7bXXWr//5z//idmzZ1tv69y5M4qKijB//nzdoJOeno7+/ftj6dKlyMvLs4aaxx9/HFVVVTh16hS2b9+OkSNHOvydWSwWLFy4EKmpqQCAsWPH4qeffsKTTz6J6upqvPXWW3j33Xdx0UUXAQDeeecddOjQwdnLQBHs/fffx+rVq9GyZUuMGDECDQ0N+Pnnn7Fr1y67ffPz83HkyBHk5OSgRYsWOH36NIqKivDBBx/g8OHD1tkhExMTMXDgQKxcuRIlJSXW95Vs8+bNqKiowAUXXMCQE+VSUlLQvn177NixAxaLxXrRu2PHDus+27dvx9lnnw1AhOjS0lLk5uZaQ86XX36J77//Hs2bN8eAAQOQmJiInTt3YtGiRdizZw8mTpxo85ie7q9VVlaGV199FfX19Zg8ebLd+R2t9uzZg/z8fOTm5uLcc8/F9u3bsXz5ctTW1qJfv35455130LdvX+Tk5GDz5s3Iz89HYmKi9f92SZLw1ltvobCwEG3atMHw4cNRX1+P9evXY968efjTn/6EUaNG2T3upk2bsHnzZvTp0wfnnXcedu7ciV9//RVHjx7F/fffDwDo0KEDRo0ahSVLlqBDhw7o27ev9f7aHjPLli1DUVER+vbti27duqGoqAhLly7FqVOncPvttwfwN0iR4tdfQ90C73j8v3FBQYHNm3LatGkAgFtvvRULFy5EWVkZSktLrT/v3LkzFi9ejKlTp+Jf//oXsrKy8MorrwR0aum2bQN2aL88bklJCRoaGjBs2DDrbXFxcRg8eDCKi4tt9lX/4Wr3x8iv8vJydOzYEdOmTcOECRPw3nvv4cILL8Sf/vQnm/+ckpOTbb5v164dysvLAQCVlZU4ePCgTRsAYNiwYdi4caPNbepPwY8cOYJ9+/Zh/PjxNv9xNjY2Ij093eFzHjlyJJYuXYpp06Zh+fLleOKJJ/DZZ59hxYoVOHnyJDIzM62BTU9OTo415GifS0lJCerr6zFkyBDrzzMyMtCjRw+Hx6PItX37dqxevRodOnTA/fffj4SEBADAJZdcgqeeespu/z//+c9o1aqVzW1NTU3417/+hSVLlmDUqFHWKuu5556LlStXYtWqVXYXgqtWrQIAu/cURafu3bvj559/xv79+9GxY0cA4tzMyspCVVWVTdDZuXMnLBaL9RP44uJifP/99+jVqxcmTpxo/cRdkiR89NFHWL58OTZs2IABAwZ4tb/Wrl27MG/ePMTGxmLatGnIysoK6O/GSIqKinDnnXeiX79+AMTfhmeeeQYFBQUoLi7G1KlTkZOTAwC4/PLLMXPmTPz8888YPXo0YmJi8Ntvv6GwsBDdunXDvffea/0Q5JJLLsEzzzyDzz//HH369LH7G/T7779j6tSp1r8zFosFL7/8Mnbs2IHdu3ejc+fOyM7ORnJysjXoXH755Q6fx9atW/Hwww9bJ4eqr6/HU089hYKCAlxzzTWGnaWVgqO2Fti4ETDiZ3gej9EZOXIkJEmy+1q4cCEAMb5i6dKlNvcZMWIE1q9fj7q6OuzevdumlB8IBQXA/v3B/3K325qjyRgkSbK7Td0VTf6ZPNblsccew5YtW3DZZZfh559/Rs+ePbFo0SLd+8r3lx9be0xnbUhJSbFuy4/9xhtvoLCw0Pq1efNmm7K51siRI7F8+XJs3LgRZrMZPXv2xIgRI7Bs2TIsXbrUWrVyRO+5yG3RPieKXDExMdZPvB359Y+PncaMGWMNOQDQvHlz3U9OtRcY8uOcd955sFgsNjMAdurUCdnZ2Vi/fj1qa2utt1dUVGDz5s3o0qUL2obqkxYKKHfOPTU5tMjjKpuamlBSUoIePXqgW7duNuMt5XNMrsrL3Xpvuukmm25FJpMJV199NUwmk814TE/3V9u8eTNeeeUVpKSk4IEHHojKkOPste3WrZs15Mj7DhgwAJIkoU+fPtaQA4iqb+/evXHq1CmcPHkSAKz/L15zzTU2ld4WLVrg/PPPR1NTk834XNmgQYNsPkwxm80455xzAIgqk6dGjRplMwNufHw8Bg4cCEmSbD6cJtKzcSMgz0V1+rRnfwtDzYDZzPi6du2K+Ph4rFixAjfffDMA0c2soKDA43Vfunfvju7du2Pq1Km46aabsGDBAlxzzTUu75eWloasrCysWLECw4cPt96+atUqDB482OH9MjMz0b59e+zatQu33HKL2+2Ux+nMmTMHI0aMgMlkwogRI/D000/jxIkTuO+++9w+llbXrl0RFxeHNWvWWD85PXHiBLZv3+4yQJFxmM1m5Obmutxv//79AMR5oaXXHae2thY//vgjNm7ciKNHj6Kurs7m5/IFi+y8887DBx98gIKCAusYnjVr1sBisbCaE6HcPffUunbtCpPJhO3bt+Oiiy5CaWkpamtr0b17d1RUVGD9+vU4evQoWrVqhe3btyMxMdH692v37t1ISEiwVgm14uLibCb58XR/2fr161FcXIzs7Gzcfffd1rEf0cTVa6u3vIXce0Gve7T8s5MnT6Jly5bW8TXqQCSTw7D8N8vV48pjuryZCtvfx6PoIo/PqakxY926XPztb6FtjycYdEIgJSUFf/vb3zB9+nTr2KbnnnsOp0+fxvjx4906Rk1NDaZPn47rr78enTt3xv79+7F27VqPugROnz4dM2fORJcuXdC/f38sWLAAhYWF1skKHHnssccwefJkpKWlYcyYMairq0NBQQFOnDhh7cqoJY/T+c9//oOXX34ZgAg/f/rTn9DQ0OB0fI4rzZo1w/jx4zF9+nS0bNkSmZmZeOSRR9wafEuRp7a2FiaTSfeiTTuDY2NjI1566SXs27cP2dnZGDx4MFJSUmA2m3H8+HGsWbMGjY2NNvcZOHAgPvvsM6xcudIadFavXm0d3EsEiK7D2dnZKCkpQVNTE7Zt2waTyYRu3bpZZ4zcvn07UlJSsG/fPvTq1cv6N+vUqVOwWCxYvHixw+PX19dbtz3dX7Z7925YLBZ07do1KkOOOxITE+1uk18nZz9ramoCIP4eqSedUJP/HukFjaSkJIfH9mahRn8fj6KLeiICJ5+FhyUGnRB55plnYLFYMHbsWFRVVWHgwIH4/vvvHf5B1IqJicGxY8cwbtw4HD58GK1atcK1116Lxx9/3O02TJ48GZWVlbj//vtRXl6Onj174quvvrKb1EBrwoQJSE5OxvPPP48HH3wQKSkp6NOnj8tq1KhRo7B+/XprqGnRogV69uyJgwcPevxpqdbzzz+P6upqXHnllUhNTcX9999vnUGLoktiYiIkSUJ1dbXNuC5AjE1T+/3337Fv3z4MGzbMrkJZUFCg2x0zMTERgwYNwooVK7B//36cPn0a5eXlGD58OGcvIhvdu3dHaWkpSktLsWPHDnTo0AHJyclITk5Geno6tm/fjmbNmkGSJJsxhfJF6fPPP+/W43i6v+zKK6/E77//jp9++glms9mt3gDkmcTERLu/OzL5dr3ARBRO5KCTkAAYbTJGfuQdIomJiXjllVdw5MgR1NbWYsWKFRg0aJD15/JYKPUAwf79+0OSJOTk5CA+Ph4ffvghSktLUVdXhwMHDuDVV1+1/sG87bbb7LrcXH311TbjWcxmM/7xj39g//79qK+vR2FhIS655BLrz3NyciBJku4UzTfffDM2bNiAuro6HD9+HMuWLXP5n+QLL7wASZLQq1cv622FhYUoLy+3GRekbftjjz2GwsJCm2NNmTLFpp9ys2bN8N577+HUqVM4dOgQpk+fjqVLl2LOnDlO20TGYbFYsGvXLuzatcvpJ5BydxL17JCykpISm++PHDkCALrT6OrdXyZXclatWsVJCKKAu+eelvyhUXFxMXbt2mUTZrp3747t27fbjc8BxN/eU6dOWSdcccXT/WVxcXG488470atXL+Tn5+Pzzz/36P6RwNvX1l3Z2dloaGjQHVcjz8Lnywyh2rG7RP524oRYLBQABg+2YP/+wL1fAoFBh4gM4/Tp0zh9+rTTfeSZrL799lub8TYnT57EkiVLbPaVZ1PTBqDt27db16vS07FjR3Ts2BG//fYbNmzYgI4dO+r2gafI4c65p9W1a1eYzWYsXboUdXV1NuuadO/eHSdPnsRvv/2GpKQkm4tduer93nvvobq62u64FRUVNos/e7q/WlxcHP7617+id+/e+PHHH6My7Hjz2rpL/nv05ZdfWruzAeLv0c8//wyz2ex0XKwrycnJMJlMOHHihM9tJdKjnsdk0KDAvl8CgV3XiCiidO/eHUOGDMHq1avxxBNPoH///mhsbMS6deusa13I+vTpg5YtWyI/Px8HDx5EVlYWDh8+jM2bN6Nfv37YsGGDw8c599xz8cEHHwAAhg4dGvDnRcaTlJSEjh07Ys+ePTCbzTYTZMjVnerqavTt29dmTGGvXr0wZswYfPvtt5g5cyZ69eqFjIwMa9WmpKQEV1xxhXXJAU/315LDzr///W/8+OOPkCQpoEtARJOzzz4bhYWF+P333/HEE0+gT58+qKurw/r163Hq1Clcd911ujM/uisxMRGdOnXCzp078e6776JNmzYwmUwYNGiQzeLjRN5Sj8/RWXM77DHoEFHEueWWW9CmTRusXLkSy5YtQ/PmzXH++ecjLy/PJugkJibivvvuw6JFi7Bjxw7s2LED7dq1w2233Ya0tDSnQWfgwIH4+OOPERMTY9PtlEite/fu2LNnDzp27GgzFqNVq1bIyMjA8ePH7VawB4ArrrgC3bp1w5IlS7B161bU1NQgJSUFLVu2xGWXXWZ3znm6v1ZsbCz++te/4o033sBPP/0ESZJw/fXX++eXEMVMJhMmTpyIJUuWYM2aNVi6dCliY2ORnZ2NCy64wGatPG/deuut+O9//4uNGzeitrYWkiShc+fODDrkF+qgk5cH6MxrEtZMkgEWIamsrER6ejoqKirsZk06ffo0iouLkZubi+Tk5BC1kIxIPneKiopw4sQJ3H777XaD1yl8WCwWFBUVAQB69uwZ8ln19uzZg+eeew5DhgzB2LFjQ9oWCqxwO/fIf/jaEjkmSUC7dsDhw0Dz5sDRoxYUF4fH+8VZNlCLmHe0UQZFUfjgOUO++PHHHwGIdXWIiIgizf79IuQAYlppzXryhmD4rmvydK7V1dVcB4A8Ig/a1VtfgkjP8ePHsXbtWpSVlWH9+vXo2bOn7kKARERERmfk9XNkhg86sbGxaNWqFQ4cOABATDPM0jM5Y7FYUF1djQMHDuDkyZM2M+FQeDOF+OOko0eP4ssvv0RCQgL69u2Lm266KaTtoeAJ9blHgcPXlkifXtAx2vvF8GN0AECSJOzZswfHjx8PQevIqE6ePInDhw/j2LFjiI2NxYQJExAXFxfqZhERERGF3KhRwNKlYrusDGjbNqTNseHuGB3DV3QAkS47d+6MvXv3YvPmzUhJSbGuFE2kJUkS6uvrYbFYUFNTg1OnTuHss89myCEiIiIC0NSkrKHTsWN4hRxPRETQkZ133nmwWCzYsmULTpw4AQMUqyiETCYT4uLiMGjQIAwfPjzUzSEiIiIKC1u3AvL6w0YdnwNEWNCJiYnByJEjMWjQIEOt2kqhk5ycjJSUFI7rMgCLxYLS0lIAQMeOHfmaUdDw3ItcfG2J9OmNzzHi+yWigg4AmM1mpKamcj0Uoggkz5RHFGw89yIXX1sie45mXDPa+yX8oxgREREREQWNHHTMZiAvL7Rt8QWDDhERERERAQBqaoDffxfbPXsCRl6mkkGHiIiIiIgAAIWFQGOj2DbyRAQAgw4REREREf3B0fgcI2LQISIiIiIiAAw6REREREQUgeSgk5gI9O4d2rb4yiQZYFXNyspKpKeno6KiAmlpaaFuDhERERFRxDl+HGjZUmwPHQqsXBna9jjibjZgRYeIiIiIiLB2rbJt9G5rAIMOEREREREhssbnAEBsqBtAROQOi8WC/fv3AwA6dOgAs5mf01Bw8NyLXHxtiWw5CzpGfL8w6BCRYVRWVoa6CRSleO5FLr62RIIkKUGnZUvgjDPs9zHa+yX8oxgREREREQXU3r1AebnYHjwYMJlC2x5/YNAhIiIiIopykTY+B2DQISIiIiKKegw6REREREQUcdRBZ9Cg0LXDnxh0iIiIiIiiWGMjsG6d2O7cGWjdOrTt8RcGHSIiIiKiKFZUBJw+LbYjpdsawOmlicggTCYTevbsad0mChaee5GLry2R4M74HCO+Xxh0iMgQTCaTYf6wUmThuRe5+NoSCe4GHaO9X9h1jYiIiIgoislBJyYGGDAgtG3xJ1Z0iMgQLBYLDh48CADIysqC2czPaSg4eO5FLr62RMCpU8DmzWK7d28gJUV/PyO+X8K/hUREfzh58iROnjwZ6mZQFOK5F7n42lK027ABaGoS264mIjDa+4VBh4iIiIgoSkXiQqEyBh0iIiIioijFoENERERERBFHDjrJycAfs0dHDAYdIiIiIqIodOQIsHu32M7LA2IjbJoyBh0iIiIioii0dq2yHWnd1gAGHSIiIiKiqBTJ43MArqNDRAZhMplw5plnWreJgoXnXuTia0vRzpOgY8T3i1cVnblz56Jz585ITExEXl4eli9f7nT/999/H/369UNycjLatWuH22+/HceOHfOqwUQUnUwmE2JjYxEbG2uYP7AUGXjuRS6+thTNJEkJOq1bA506Od/fiO8Xj4POxx9/jClTpuCRRx7Bhg0bcN5552HMmDEoLS3V3X/FihUYN24cxo8fjy1btuDTTz/F2rVrMWHCBJ8bT0REREREntu9G5DrDoMHAwbJLh7xOOi8+OKLGD9+PCZMmIDc3FzMmTMH2dnZmDdvnu7+a9asQU5ODiZPnozOnTvj3HPPxZ133omCggKfG09E0cNiseDgwYM4ePAgLBZLqJtDUYTnXuTia0vRzNPxOUZ8v3gUdOrr67Fu3TqMHj3a5vbRo0dj1apVuvcZOnQo9u/fj8WLF0OSJBw+fBj//e9/cdlllzl8nLq6OlRWVtp8EREdP34cx48fD3UzKArx3ItcfG0pWnkzEYHR3i8eBZ2jR4+iqakJmZmZNrdnZmbi0KFDuvcZOnQo3n//fdx4442Ij49H27Zt0bx5c7z66qsOH+fpp59Genq69Ss7O9uTZhIRERERkRPqoDNoUOjaEUheTUagHYAkSZLDQUlFRUWYPHky/vGPf2DdunX47rvvsHv3bkyaNMnh8WfMmIGKigrr1759+7xpJhERERERaTQ0AOvXi+0uXYCWLUPbnkDxaHrpVq1aISYmxq56U15eblflkT399NMYNmwYpk+fDgDo27cvUlJScN555+GJJ55Au3bt7O6TkJCAhIQET5pGRERERERu2LIFqKkR25G4fo7Mo4pOfHw88vLykJ+fb3N7fn4+hg4dqnuf06dPw2y2fZiYmBgAohJERERERETBE+kLhco87ro2bdo0vPnmm3j77bdRXFyMqVOnorS01NoVbcaMGRg3bpx1/yuuuAKff/455s2bh127dmHlypWYPHkyBg8ejKysLP89EyIiIiIicilago5HXdcA4MYbb8SxY8cwa9YslJWVoXfv3li8eDE6/bHKUFlZmc2aOrfddhuqqqrw2muv4f7770fz5s1x/vnn49lnn/XfsyAiIiIiIrfIQScmBhgwILRtCSSTZID+Y5WVlUhPT0dFRQXS0tJC3RwiCgFJktDQ0AAAiIuLM8yqzGR8PPciF19bikbV1UB6OmCxiJAjT0rgSji9X9zNBh5XdIiIQsFkMiE+Pj7UzaAoxHMvcvG1pWi0bp0IOYBn3daM+H7xanppIiIiIiIynmgZnwOwokNEBmGxWFBeXg4AaNOmjd1sjkSBwnMvcvG1pWjkbdAx4vsl/FtIRPSHo0eP4ujRo6FuBkUhnnuRi68tRRs56KSkALm5nt3XaO8XBh0iIiIioihw6BAgT448cKCYdS2SMegQEREREUWBtWuV7UgfnwMw6BARERERRQX1+Jyzzw5dO4KFQYeIiIiIKApE04xrAIMOEREREVHEkyQl6LRtC3ToENr2BAODDhERERFRhNu5Ezh5UmwPHgyYTCFtTlBwHR0iMgSTyYSuXbtat4mChede5OJrS9HE125rRny/MOgQkSGYTCYkJiaGuhkUhXjuRS6+thRN/BF0jPZ+Ydc1IiKiKLFtG3DttcDcuaFuCREFmzroDBwYunYEE4MOERmCxWLB4cOHcfjwYVgsllA3h6JIJJ17zz4LLFoE3Hef0lc/mkXSa0vkTH09sGGD2O7eHWjRwvNjGPH9wqBDRIZx5MgRHDlyJNTNoCgUKefe/v3i38ZGoKIitG0JF5Hy2hI5s2kTUFcntn2ZVtpo7xcGHSIioiihDjcNDaFrBxEFV7StnyNj0CEiIooSDDpE0YlBh4iIiCJaZaWy3dgYunYQUXDJQScuDujXL7RtCSYGHSIioijBig5R9KmsBIqLxXa/foDBZoj2CYMOERFRFGhoAE6ftv2eiCLfunWAJIntaOq2BjDoEBERRYWqKtvvGXSIokO0js8BgNhQN4CIyB0mkwlnnHGGdZsoWCLl3NNOJ80xOpHz2hI546+gY8T3C4MOERmCyWRCcnJyqJtBUShSzj1t0GFFJ3JeWyJn5KCTmgr06OH9cYz4fmHXNSIioiignnENYNAhigYHDyoLBQ8aBJij7MqfFR0iMgSLxYJjx44BAFq2bAlztP21ppCJlHOPXdfsRcprS+TI2rXKtq/jc4z4fmHQISLDOHz4MADxB5YomCLh3GPXNX2R8NoSOeLviQiM9n4J/yhGREREPmPXNaLoE80zrgEMOkRERFGBFR2i6GKxKF3XsrKA9u1D255QYNAhIiKKAhyjQxRdduxQ3vfRWM0BGHSIiIiiAruuEUWXaO+2BjDoEBERRQV2XSOKLr/+qmwz6BAREVHEYtc1ouiirugMHBi6doQSp5cmIkMwmUzIycmxbhMFS6Sce+y6Zi9SXlsirbo6oLBQbJ95JpCe7vsxjfh+YdAhIkMwmUxo1qxZqJtBUShSzj12XbMXKa8tkdbGjcp73F/d1oz4fmHXNSIioijAoEMUPTgRgcCKDhEZgiRJOH78OAAgIyPDMGVzMr5IOfe0Xdc4RidyXlsirUAEHSO+Xxh0iMgQJElCWVkZAKBFixaG+ANLkSESzj1J4hgdPZHw2hLpkYNOfDzQt69/jmnE9wu7rhEREUW4U6eApibb2xh0iCLTyZPAtm1iu39/ICEhlK0JLQYdIiKiCKet5gDsukYUqQoKlO1oHp8DMOgQERFFPO1EBAArOkSRihMRKBh0iIiIIhyDDlH0UAeds88OXTvCAYMOERFRhNPrusagQxR5JAn49Vex3bw50LVrSJsTcgw6REREEU6vosMxOhRJNmwQi2RGuwMHgEOHxPagQYA5yq/0Ob00ERmCyWRCp06drNtEwRIJ5x4rOvoi4bUlYM0aYMgQwGQSA/HPOivULQqdQI7PMeL7hUGHiAzBZDIhNTU11M2gKBQJ5x7H6OiLhNeWgMWLxb+SBHz5JYOOLBBBx2jvF68KWnPnzkXnzp2RmJiIvLw8LF++3On+dXV1eOSRR9CpUyckJCSgS5cuePvtt71qMBEREXmGXdcokm3erGyvXRu6doQDddAZNCh07QgXHld0Pv74Y0yZMgVz587FsGHDMH/+fIwZMwZFRUXo2LGj7n1uuOEGHD58GG+99Ra6du2K8vJyNPIvLBF5QJIknDx5EgDQvHlzw5TNyfgi4dxj1zV9kfDakm3QKSgQlZ1ofCmbmpQ1dLKzgXbt/Ht8I75fPA46L774IsaPH48JEyYAAObMmYPvv/8e8+bNw9NPP223/3fffYdly5Zh165dyMjIAADk5OT41moiijqSJOHAgQMAgPT0dEP8gaXIEAnnHruu6YuE1zba1dQAJSXK90eOAKWlwB9DSaLKtm1AVZXYDsT6OUZ8v3jUda2+vh7r1q3D6NGjbW4fPXo0Vq1apXufr776CgMHDsRzzz2H9u3bo3v37njggQdQU1Pj8HHq6upQWVlp80VERETeCZeg8/PPwLRpwO7dwX9sikxbtwIWi+1t0dp9jQuF2vOoonP06FE0NTUhMzPT5vbMzEwckuey09i1axdWrFiBxMRELFq0CEePHsVdd92F48ePOxyn8/TTT+Pxxx/3pGlERETkgN7nhcHuQV5TA1x7rQhdBw8CH30U3MenyKTutiYrKACuvz74bQk1Bh17Xk1GoC1VSZLksHxlsVhgMpnw/vvvY/Dgwbj00kvx4osvYuHChQ6rOjNmzEBFRYX1a9++fd40k4iIiKBUdNRragS7orNundKOnTuD+9gUufSCTrRXdEwmIC8vtG0JFx4FnVatWiEmJsauelNeXm5X5ZG1a9cO7du3R3p6uvW23NxcSJKE/fv3694nISEBaWlpNl9ERETkHTlgpKUBcXFiO9hBR93DXa8rHZE31EEnOVn8W1Bg350t0tXWKgum9uwJGGwW6IDxKOjEx8cjLy8P+fn5Nrfn5+dj6NChuvcZNmwYDh48iOrqautt27dvh9lsRocOHbxoMhEREXlC7rqWlgbE/tFpPdhd1xh0KBDkoJOWBlx4odiurIy+qmFhofKeZrc1hcdd16ZNm4Y333wTb7/9NoqLizF16lSUlpZi0qRJAES3s3Hjxln3v/nmm9GyZUvcfvvtKCoqwi+//ILp06fjjjvuQFJSkv+eCREREemSg0V6emgqOpIErF5t3x49mzcDTzwB7N0b+HaRsVVWihnWAKB3b9sL/GjrvsbxOfo8nl76xhtvxLFjxzBr1iyUlZWhd+/eWLx4MTr9MY9fWVkZSuWzDkCzZs2Qn5+Pe++9FwMHDkTLli1xww034IknnvDfsyCiiGcymZCdnW3dJgoWo597DQ1iIgAgdF3Xdu0CysuV7+vrRVebxET7fW+8ESgqAtasAb75JrDtMvprG+22bFG2e/UCBg5Uvl+7FrjlluC3KVSCEXSM+H7xOOgAwF133YW77rpL92cLFy60u+3MM8+06+5GROQJk8lkM9aPKFiMfu6pZ1wLVUVHbwWKigr7oNPUJEIOABQXB75dRn9to5066PTubR90ookcdBISgD59AvMYRny/eDXrGhERERmDuptYenpoxuiou63J9LqvqW+TFz4kckQ9EUHv3kDLlsAZZ4jvN2wI/ji0UDl+HNixQ2yfdZbyYQYx6BCRQUiSZJ1yXpKkUDeHoojRzz11eAhV1zVHFR2tEyeU7WAEHaO/ttFOG3QAYNAg8W9NjVIdjHQFBcp2IMfnGPH9wqBDRIYgSRL27duHffv2GeYPLEUGo597oe66VlUFbNpkf7uroFNbG/hP5I3+2kY7Oei0bg20aSO2o7H7WrAmIjDi+4VBh4iIKIKFuuvab78pa5rExOi3S6YOOgCgWpmCyMaRI8Dhw2JbruYASkUHsK10RDLOuOYYgw4REVEEC3XXNXW3tSFD9Nsl0wYdjtMhR7QTEcjOOguQJwSLhoqOJAG//iq2W7QAunQJbXvCDYMOERFRBAt11zX1RARjxijbDDrkC73xOQCQmgqceabY/v13oK4uuO0KttJSZer2wYOVkEcCgw4REVEE03Zdk4OOxaJ0KQsUi0UJOm3aAAMG6LdLxqBD7lIHnV69bH8md19raAA2bgxem0KB3dacY9AhIiKKYOqKTmqqMkYHCPw4na1bgZMnxfbQoSJo6bVLxqBD7tIuFqoWTeN0GHScY9AhIiKKYDU1ynZysu0aG4HuvqbutjZkiG3QYUWHvCVJSkWnQwegeXPbn6uDTqSP01EHHfXzJiHW9S5ERKFnMpnQvn176zZRsBj93KutVbYTE4MbdNQTEQwdKiZDkIVD0DH6axutDh5UKoXq8Tmyfv1E5bKxMbKDTmOjUrHq1AnIzAzs4xnx/cKgQ0SGYDKZ0KJFi1A3g6KQ0c89ddBJSgpu1zU56MTFAXl5tsEqXIKOkV/baOVoIgJZYiLQpw+wYQNQXCymKW/WLHjtC5biYuD0abEdjG5rRny/sOsaERFRBAtVRef4cTFGBxCTECQliYtN+YPgcAg6ZEyugg6gdOOyWETgiUQcn+Magw4RGYIkSaiqqkJVVZVhVmSmyGD0cy9UQWfHDmU7L0/8azYr3dfCIegY/bWNVu4EnYEDle1I7b4W7KBjxPcLgw4RGYIkSdi7dy/27t1rmD+wFBmMfu6FKuioJ0FQT0Igb/sadFavBoYPB1591fs2Gv21jVZy0DGZgNxc/X2iYeY1OeiYzWKh1EAz4vuFQYeIiCiCaYNOsMboaB9X5ijoWCz2t1VXOz7+k08Cy5cDDzygjFOgyGexAEVFYrtLFzGToJ5evZTzLhIrOqdPA5s2ie1evSJzDJI/MOgQERFFMDlwmEyimhOsio6roFNbC9TXK7dXVIhpg9WcVXQOHhT/1tcDR4741lYyjj17lGCrXT9HLS5OWaB25077aqHRbdgANDWJbY7PcYxBh4iIKILJgSMxUQk7slAGHcC2gqN3Ieos6MjTCwMMOtHEnfE5MvU4nUjrvsaJCNzDoENERBTB1EEHCJ+ua4BvQUe9P4NO9PAk6ETyOB0GHfcw6BAREUUwbdCJhIqOdjzP0aPetZGMx9ugE2njdOSgk5Tk+vcQzRh0iIgoKjQ1AV99FblrajgSiUGnstJ2PE+0VXQKC4GSklC3IjTkoBMbC3Tv7nzf7t2B1FSxHUlB5+hRYNcusZ2XZ1ulJVsMOkRkCCaTCe3atUO7du1gklccJPLAwoXAVVeJaVjlgezuMPq55yzoGLXrmnp8DuB90DHia/vJJ2KQfdeuwIUXAl9+qQxKj3QNDcoitD16APHxzvc3m5U1nPbvBw4dCmz7gkUd2oLZbc2I7xcGHSIyBJPJhJYtW6Jly5aG+QNL4WXCBGX7nXfcv5/Rzz1nY3QCWdGpq1O2wznoGO21zc9Xtn/6Cbj6ahF6nn8eOH48ZM0Kih07lHPW3e5akThOJ1Tjc4z4fmHQISIiilAWizKFcyi7riUkKNueBJ36etspqB3tG01jdEpL7W/bswd48EGgfXsR6DduDHqzgsKT8Tky9cxrkdJ9jRMRuI9Bh4gMQZIkVFdXo7q62jArMlP48uTDSCOfe3pVFSOM0WnXTtnWq+r4q6JjxNdWDjrJycD//geMGaP8rLYWeOstoH9/YPhw0c0tkK9xsG3Zomx7U9GJhKAjSUrQadUKyMkJ5mMb7/3CoENEhiBJEvbs2YM9e/YY5g8sRQYjn3t6YcMI00t37KhsV1fbH9ufQcdIr60kKUGnY0fg0kuBxYuB7duB++4D0tKUfZcvB268EejcGXjiCaC8PDRt9id1RcfZYqFqOTlAy5Ziu6DAflFao9mzR6lgDh7s2Yc2vjLa+wVg0CEiIvLaihXiQnLixFC3RJ9e2DBCRUcddPQqOuHUda2pKXgXz8ePA6dPi23176hbN2DOHODAAWDuXKBnT+VnBw4Ajz4KZGeLCTmMTA46iYnAGWe4dx+TSanqHDmi3/XPSNhtzTMMOkRERF664QbxCeubbwLr1oW6NfaMGHTMZiArS7ndna5rJ06EpovW118DzZqJykowwo76Il0ddGTNmgF/+5sIBD/+KGYZNP9xpVdfL362e3fg2xkINTXAzp1iu2dPICbG/ftG0jgdBh3PMOgQkeGoL46IvOGv7h5lZcq2evxAuHAVdMKx61rz5rZdsNyp6ADAsWNeNdNrTU3AlCnieX73HVBcHPjHdBV0ZCYTcMEFwBdfiPV2/vQncXttLTB1akCbGDBbt4rJNQDPF8iMpJnX1EFH/bxIH4MOERnC9u3K9v/9X+jaQeSIJ2vzBIurMTqhqOikpipBUy/otGihLPIIuFfRAYLffe3zz5VFGwHb7UBxN+io5eSIimPbtuL7L78Evv3W700LOG9mXJNFSkWnsVGpHJ9xhpiMgJxj0CEiQ1i2TNlesCB07SBy5MCBULfAXk2Nsh0uXdfMZiXIyEHHYlHCS/Pm3gUdbyck8IYkAc89Z3tbMLqE7dunbLsbdABRIXvhBeX7yZNtZ+QzAl+CTlaW0h2yoECpDBnNli3Ke5rd1tzDoENEhmCQtckoyrRpo2zv2ROyZjgULmN01OvoAEr3NTnoVFYqY1xatBBjTWTudl0LZtBZtsy+C1SwKzrZ2Z7d9+abxZTTgBjrog4+RuBL0AGUbl6VlWLhUSPi+BzPMegQkSFIEjB7diZmz86MqHUhKDQ8Dc6ZmZnIzMy0u715c2VbHigdTsJheumEBPvftzboqIOLt13XvA06jl5bZ7TVHCD4QadDB8/uazIBr72mDOJ/8klg717/tS3Q5DFwaWmeP3cgMsbphEPQ8eb9EkoMOkRkCJJkxsKFrbFwYWs0NvJPFwWP2WxG69at0bp1a5jNtueeOkjs2iUGqIeTcKjoqLutyeSgU1Mj2uCPoOPNGB1nr60jmzYpY1w6dQLi48V2MINOZqb+79WVPn2Ae+4R2zU1wLRp/mtbIFVWKqGsVy/vKvyRsHCoHHRiYoABA4L/+N68X0LNGK0kIiKCWPTw++99r0T4K5Cox8DU14ffOJ1QBh15DIi22xpgP/Oap0EnlF3X1F2+pk1TVqbftSuwU0w3NCgTXngyPkfr8cdFUALEhAo//OB72wKtqEjZ9qbbGgDk5SnbRgw61dVK970+fYDk5NC2xygYdIjIICT06nUavXqdhtlsjBWZyb8sFmDIEOCSS4B//tO3Y3kSlCRJwunTp3H69Gm71cDVQQcQU/mGk3CYXtpZRQcQQUddoXEVdOrrlUUz1evteBN0nL22evbtAz74QGxnZADjxysLV54+LYJ4oBw4oAQpX4JOerpt17t77w3/iQl8HZ8DAC1bKq/Vhg2BPfcDYf16ZRKFUHVb8/T9Eg4YdIjIEMxmCR99tAsffbQLCQnG+ANL/rVjh9I9aNYs347lSSVDkiTs2rULu3btsvnPXZKUC26ZEYJOsKeXdifoOKvoVFfb3lcdirp1U7a96brm6LV15OWXlQvku+8GUlKUi2cgsN3XvJla2pGxY4Fhw8T29u3ASy/5drxA80fQAZTuazU1tlUiIwiH8Tmevl/CAYMOERkCZ10jf54D/rjAb2iwn6bW2YQE9fXAkiXBXfA23MfoAJ53XVMHnQ4dlK5xge66dvIkMH++2E5MVMa6GDHomEzAv/4lpvoGRIVUPXV1uPFX0DHyejrhEHSMiEGHiIiijj8u8LXd1gDnFZ0ZM4DzzwdGjAjsWA41ddBJShL/BiPoSJJvQcfZ9NLabm7yoomBDjrz5yvVpdtuU6YWVwedQK6l48+gAwD9+gF33SW2T58G7r/f92MGihx0Wre2ndLdU0aeeU0OOikpQM+eoW2LkTDoEBGRIfizouOP/vnabmuA86Dzv/+JfzduDF5VJ1TTSzc0KGHOm6BjNosLOsA+6Gj3bd1abB89GrgAWVcHzJkjtk0m29nKOndWto1S0ZH985/K7+/TT4Eff/TPcf3p6FHg8GGx7Us1BwDOOkv5O2Kkis7hw8qsc3l5yhTh5BqDDhEZAruukT8FqqJTWam/b2Oj7UWwfOEWaKHquqb3uGqugg6gdF9zVtFp3ly5UG9sDFyAfP994NAhsX3ttbZjg4wcdJo3B559Vvn+3ntFF8twIq+fA/gedFJTgdxcsf377+E/CYNMHcrYbc0zDDpERGQI4TZGRy/o6N0GiE9j1Y8ZyNm51IwcdOTua666rslBBwhM9zWLBXj+eeX76dNtf56eLmb0AoITdBISbJ+zr269FTjnHLG9dauYcCGc+Gt8jkwep9PQIKqrRsDxOd5j0CEiIkPQDvz3RaCCjvoCX237dtvvwyXoBKrrmrdBx2RSfqau6Ki7pKlDUfPmyhgdIDBB55tvRAAAgOHDgbPPtt9HHqezf3/gqgTyZAEdO/o39JvNYmIC+ZiPPx5e60Gpg06vXr4fz4jjdBh0vMegQ0SGMXdua8yd2zqgM0VR+PLXIp+A+xf4R48Cr74KSJJYDVxNb4yOo6CzY4ft96EMOsGYXtrboJOerswEJgedxkbb8OCo6xrg3RTT8krvjqirOQ8+qL+PHHQkSRlL4U8VFUq3yOxs/x//rLOASZPE9qlTwAMP+P8xvBXIoGOEcTqSpASdNm38123RW67eL+GGQYeIDEGSzJg3LxPz5mWisZF/uqKRP4NOQ4P4hHzQIOCKKxwHn0mTgClTzBgxIhOtWmXCbFbOPUcVHb0B8eFa0Qm3oCN3WwMcTzHtaDICwPOKjtlsRmZmJjIzbV9b2a+/AitWiO2ePYExY/SPE+gppgMxPkfriSeULngffSSmQg81SVKCTocOItj6ql8/JewbIeiUlCjn/ODBoR2v6ur9Eo6M0UoiIop6vnRd0963oUGsal9QILomvfmm/v0++0z8e+KE/UB3vaBjseiHB0cVHe1imP4WDkFHXudGTR10Tp5UqjTuBB1tRSeQXdfk1x8QM605urYL9BTTwQg6GRnAM88o399zT2DXWXLHwYPK6+2P8TmAeB/06SO2i4sD/x70Fbut+caroDN37lx07twZiYmJyMvLw/Lly92638qVKxEbG4v+/ft787BEFMVMJgldutSiS5damEzGWJGZ/MuXio72gq2hAfjpJ+X74mLH95XPvaqqWpvVwPW6rgH63de0FZ3Dh4HXXhMX/Lfd5rztvgjV9NLqrmZ6FR11iDlwQHltHQUd9cWoP7uuSZKE2tpa1NbW6q70/vPPyvZVVzk+TiRUdADgjjuUi+miItFtM5T8PRGBTO6+ZrEAGzb477iBEE5Bx9X7JRx5HHQ+/vhjTJkyBY888gg2bNiA8847D2PGjEGp+l2oo6KiAuPGjcMFF1zgdWOJKHqZzRK++GInvvhiJxITjfEHlvzLn0GnsdH1mityFSgxUZx7VVU7bf5zdzTDmjbo1Nbaj9soLxdT+VoswDvv+LdbnqO2hFPXtZgYJcjs2aPc7knXNZMJSEvzreuaJEnYuXMndu7caXfhduIEsH692O7Xz7ZypBXoKaaDFXS0ExM89hhQVha4x3Ml0EEHCP/ua+qgo253KDh7v4Qrj4POiy++iPHjx2PChAnIzc3FnDlzkJ2djXnz5jm935133ombb74ZQ4YM8bqxREQUvfxd0VHT6/d+7JjzY7obdHbtsg9V8posro7lK7ktZrNSyQmHoAMo3dfUlTFPuq7JExcEanrppUuV1+38853vm52tLOJo5KADiOmXJ04U21VV9tNpB5M/19BRk6eYBsI76DQ0KGG7a1fRvZA841HQqa+vx7p16zB69Gib20ePHo1Vq1Y5vN+CBQtQUlKCmTNnuvU4dXV1qKystPkioujGBUPJ30HH1QeS2jCipQ4nKSnKtjboaLutAWKAsZqjbnC+ktuSmKi8h4LRdc2ToKPmSdCRB6a3aKE8N38GHXW3NVdBJzYW6NRJbOsFW1+pg04gZl3Teuop5aL6/feBX34J/GPqkSs6JpOy0Kc/9OqlnJfhPMX0pk1KN9BQd1szKo+CztGjR9HU1ITMzEyb2zMzM3HIwf8IO3bswMMPP4z3338fseq/rk48/fTTSE9Pt35lB+NdTUREYc2XyQhcVXT0HD7s/OeOKhHa6oxe0NE+l2AEHZnJpFQfwqGio+Yq6EiS/QxtMTHKbGHeTC/tiBx0YmLE+jmuyON0KiuB48f91w5ACTqtWgHJyf49tp6WLUXYkd1zT+BCsSMWi1LR6dLFv887Lg4YMEBs79xpO5NfOAmn8TlG5dVkBCbNR6uSJNndBgBNTU24+eab8fjjj6N79+5uH3/GjBmoqKiwfu2TV8kiIqKo5e8xOmp6FUNXQUcdaNRdSrQVHfWMax066B8rmEEHULqvhXPQadZM2ZaDzunTymunnmpY7r7mr4pOWZkYjA+IcRFpaa7vE6gJCRoblQU8g7mGyoQJQF6e2N60CZg/P3iPDYixW/L7wh/r52ipu6+Fa1WHQcd3HgWdVq1aISYmxq56U15eblflAYCqqioUFBTgnnvuQWxsLGJjYzFr1ixs3LgRsbGx+FldF1ZJSEhAWlqazRcRRTd2XaNAjtFRO3oUePll21nZ9LgbdA4eVLYdDSYOVNCR22jEoKNX0dGuoSOTJwo4dco/453Ua8i46rYmC1TQKStTzv1gBp2YGDExgew//wneYwOBm4hApn4vhnvQiY0FOGGxdzwKOvHx8cjLy0N+fr7N7fn5+Rg6dKjd/mlpadi0aRMKCwutX5MmTUKPHj1QWFiIs88+27fWExFR1AjWZAQPPABMmQIsWOD8mOoLavVFtzboqKdD7tFD/1juBp3du0WXm6uucu/34aiiI/ckN9oYHe3U0jJfppjWo/4c1t3JYgO1lk4wJyLQOvts4MwzxXZBgQiSwRLMoBOOExJUVSlVxb59gaSk0LbHqNwbNKMybdo0jB07FgMHDsSQIUPw73//G6WlpZg0aRIA0e3swIEDePfdd2E2m9Fbc3a2adMGiYmJdrcTETkjScCCBeJj21AvYkf+MXMmkJ8PzJ3r3qeVwarovPOO/b4LFrTC8OFAz57K7e6O0ZEvzlNSgPbt9R/T3aDz7rtAYaH4+uUXYNQo5/sbueuat0HnyBHPBuy30pk3Wq7mJSQA7k4WG6iKTiiDDgCMHAls3SpC8apVwEUXBedxAx10uncX51hVVXgGnXXrlEktwqkuoPd+CWceB50bb7wRx44dw6xZs1BWVobevXtj8eLF6PTHdCNlZWUu19QhIvKUJJnx4ottQ90M8pOSEmDWLLE9apR7g4H9ORmBJ5WMxkZx7iUn2y4a6W7XNfUsYW3a6D+Gu0FHXa1wNStcY6PyPCMl6DjquubtFNNmsxlt29r+Xdm9W1nbZ+hQ9z9JD1TQUQ9TDsXcTCNGAK+/LraXLg1+0ImNFaHE38xmMQZp6VJg/37xfmobRv/FhOP4HL33S7jzajKCu+66C3v27EFdXR3WrVuH4arpSBYuXIilS5c6vO9jjz2GwsJCbx6WiIgihHrcivpTemeCvY6OVnW17ff+Djo1NcC//23bbUpLPc2yq3V+5GlpAcdBJxhd1xIS9PfxZ0VH/SGzr13XPJlWWq1FC+U5RVJFZ8QIZXvZsuA8ZkODqCIBortnfHxgHiecx+mEY9AxIq+CDhFR8EnIyqpHVlY9TCZjrMhM/hWsrmvanhkmkzj3zOZ6m9XAHXVdU1/k19cr+7kKOq+9Btx5pxgTcscd+uMhPAk6zqoq8hidcKvoqG9TBx05ZLrbdc1dkiShvr4e9fXKa6uehMLd8TkyuapTWuq/322og067dkpF5bffAjdxhtrOncrvL5AjHcJ5nI4cdFJTHY/tCza990u4Y9AhIkOIiZHw/ffb8f3325GYaIw/sOSYN7PoaYOOJ//PelLRUYcWAEhMFOfe7bdvt/nPXV3RUV90qy/yKyps99GZoBSAuHhUj0lYsED0y9eux6IOOq7WanEWNgLddc1ZNUmmDTrp6cr6PoAY0yS/LoHquiZJErZv347t28VrK0lKRadZM9spiN0hB52mJtsuZ76Qg05cXOi6VslVnYYGYM2awD9eoMfnyNSvbzgFnbIy5fwZOND2fRFK2veLETDoEBFR0Pkj6HhS4fFkjI67F/9y0ElKsh3HoQ5A2qDTooX+Rcvp07YhBhCLJX7zje1t/qrohOMYHW3ANJmUtXQ86brmy1o6xcXK+knDhyu/J3cFYpyOHHQ6dBDjSkJB3X3NyegEvwlW0MnJURabLSjw7MOTQFKHLnZb8w2DDhERGYJ2MgJPxpe4quioj+1N0FFfzKsv8rUX5mazbfVBdvq0/RggwLbbEuD/oBPK6aW1S+Rpgw6gdF/TCzqOKjq+jNHxZlppNX8HnaoqpYoVim5rsmCP01EHnUAsFiozmZTua0eO2L/fQoXjc/yHQYeIDIELhkaWUFd0nH3v7sW/PFZBW9FxFnQA/e5rehUdADhwwPZ7f4/RaWwMzKfY/qjoAPYVHXXXNX+M0dFSj8/xZCICmb/X0lF3fwtl0OnQAejSRWz/+qv9hBv+JgedxETb32kghGP3NQYd/2HQISIiQ9AGG39WdNTHUv9swADHx/SmogMAkyaJiop6fRZHQWf/ftvv/T1GBwhMVcdfQUdd0ZEkx13XEhKUfb0NOk1NSresjAyxSKOn/F3RCfVEBGpyVaeuLrDjdGpqxGQEgFi3KtDjU8Jt5jWLRQlc7do5XnuL3MOgQ0REQefNmjj+rOhoqxjqn8vbubnA+vX6Xc0AJegkJ9tezKvH6KgvzOUL+0mTxIW7emFSdde11q2Vios66EiS/7uuAYEPOo6ml/ak65rFIn6vckUnPt5+fRt5nI63Xdd+/115vUaN8m48TMeOyv0iNegAge2+tnWr8vchGGvLh9vMazt3Kufh4MHszeArBh0iMgT+sY8s3kwV7c+go6VX0ZHDQHKy/f4Wi3Ix72lFBxAX/+oLdXVFJz1d+RRXHXTq6mzbefKk85DiTtc1IDATEsiPHRtr+1jaNqSkKN87CzqA+P2o1yTS/k2QA+nx496dX+pB9t6MzwFEAOvQQWwz6HgnWBMRyNRVk4IC3xYm9odff1W22W3Ndww6RGQIkgR89FEGPvooI2ADqCl4vHkN/TkZgbOfy8eVL9CTksS59+mnyqqg6hDh6RgdmTpAqYNOaqpysXz0qHI8va5t6jErWto2qqkrOoEMOo66rcnU3dc8DTpactCRJNfVLrWMjAxkZGR4vVColtx97fhx9xfDdSScgk6nTmKWMgBYvdp2CnF/2rJF2Q5G0AGUcTqVlcCOHcF5TEfCfXyO/H4xCgYdIjIESTLjySez8OSTWWho4J8uowu3io78c0lSgo4cBhISxLk3a1YWJEmce+ruadqua94EnRMnlDY0a6YEHUCZkEAv6Di7oHe365pRgk5FhTJdt96+6imm3e2+ZjabkZWVhVatsrB0qXhts7KUBTK94c8JCdRBJzvbt2P5g1zVqa21vSD3p2BXdIDwGqej/r16uo5ToMnvl6ysLJhDNde5h4zRSiKKeuy6Flm8qej4czICR+1RH1MOA/LMX4Ay05o66HjTdQ0Q3dfk87q8XLk9NdV2ALLcfU0v6DibkCAcxug4Gp8j8yToqGegc1bRATyfkODXX5XX9oILfPt7488JCeSg07y57e8iVILRfU0OOmlptoE/kMJlnE5dHVBYKLZ79NA/z8kzDDpEZBASWrRoRIsWjQDCZFU38po/gk4gKjrq/eQwkJKinHvV1eLcky+KAfug42gyAu1Fi8mkVHXkRSoB265rgPOg421FJ1hjdPxZ0VFPteyvoCNJEhobG7FihfJ3xZdua4D/KjoWi/Lah7rbmmzkSGU7EEGnshLYu1ds9+oVvA+48vKU7VAGnd9/B+rrxXY4dluT3y+NjY2QwmV1VRcYdIjIEMxmCb/8shW//LIVSUnG+ANLjjkLKStWAD/+aD8zWrArOnIYaNFCOffkoKOt6MTGKrNtOaroaKdTBpSgc+qUclswgo4Ru66pg47evt4Gna1bt+KKK5S/K6NGuXdfR/xV0Tl8WHltwiXo5OQoXehWrlQuyv2lqEjZDla3NQBo2VJ53TZsCNxCuq6E+/gc+f2ydetWBh0iIn9i17XI4uhCYv164LzzgIsusl2lHrCfjEAbfF56CZg2TXwqrOVLRUc7aQBgP0bHZFIG/OsFneRkMSOXlt6MbtoxOoEOOoG4qJMHqrsKOv37i39TUoCuXe1/rg466vEqehUdb8boaHXpIgbd+8JfQSecJiKQmUxK97WaGv+PZwnF+ByZ3H2tpsY2cAVTuAcdI2LQISIiv1u1Cvjb34B16/R/7qii88gjyvbddzu/j/r7n38WIeell4CnnrI/ri9BRz1GR668aCs6gHJRrw4Y8uB5R33t9YKOtqLjbDICb8foBLLrWmOjEp5cBZ377gPeeAP46Sf9Ko36dx+Iio6Wt9NKq7VqpbQ70oIOENhxOuEQdIDQdV+Tg05cHNCvX2jaEGkYdIiIyClJAsaOBYYMAfbscb3/ihXAsGHA66/bhxWZO1UET7quffGFsv3ss/bHcrfrmquKjhx0tGN0AOWiXm+MjqdBp21bpSuc0bquqacddhV0kpOBCROAs8/W/7knFR1/BB1fu60BouohV3X27PFuhkEgfINOIMfphDLoqGc4C0XQqagQi6UCotLpaiIPcg+DDhEROfX558B//gOsWQPcdJO44HZ2ITlnjrKtXvxOzd+TETi7qAdcjyWQL/b1xuioF7X0pKLT2AhUV4ttT4NOXJwIO4BnQWfnThEwjx61bWMwg46r18ITnsy65o+ua0OGeHc/rc6dxb8NDbbt9kS4Bp0uXcQU3ID4UMOf54+8hk7r1kCbNv47rjvOOkvpIh2KKabVj8lua/7DoENEhsAxOqGzbZuyvWaN6FbVoQOwcaP+/uqLTEcXuo4+5Va/zq4qOnv3Kvu4urj2paKjF3TUFR05rGjH6Mjd1gDPgw6gdF87dEi0y1XQsViAMWNEl8G77grdGJ1ABR3166/XdS0tTXle3lR00tP9Fyj8MU4nXIOOepzOqVNiXJ0/HD0qznUg+NUcQJxrubli+/ffA7cgqiMcnxMYDDpEZAgGmeAlIqnHcwBisH99PTBxov7+6gtnR2vKeXNxrZ2M4PrrgUmTxLavQcf59NLKbXLQUc+SJv9cW9FxNrW0zNFkBIASdCQJKCtzHXSKikRFBxAr14dqjE6ggo6a3u/TZFK6r7kbdNShvG9f/32g4s+gExMDtGvne5v8KRDd1+RqDhCaoAMo3dcaGhx/kBMoDDqBEet6FyKi0JMk4MsvmwMI3dSf0SomRv92+dNXLfWFszacyBy9hs4uNPWqQJ9+Csyf77yblrZNztrjKOjI597jj4vAoQ7e2qDT1CSO42pqacC9ig4guq+5mozgl1+U7YMHbWefi4Sua2qOgmPr1uK5Hz0qXiNXweX334FffxUH69XL62ba8cdaOnLQad/e/sOGUFNPSLB0KfDgg74fM5Tjc2SDBgHvviu2CwqCGzjkoJOWBnTvHrzH9VRzg61iGmZvHSIiR8z4v/8L0jLZZMPRRZaj29UXzo4uorWhxWIR1R9Puq4B4mJekmwv+PUG8bpb0VEHMGUyAttz74EHgDvuUPaTw4r6or621vuKjnxx3769cps26LRsKao56oqO+pN1i0Wp7mjbBhiz65qaXtc1QBmnU18vfl9pac6PX1iovLYLFnjZSB2+VnROn1aqTeHUbU3WvTuQmSnW+lmxQpxDvoYxddDxZ+j0RKhmXjtwQAR0uQ2OKuGhZjab0aGDsf4fDtNfJRERhQtHFzCOKj3qUNHUpB9QtBfX7sxMpbdPU5O4KFR3VVJfaOu1SY9e1zX5eaunOJaVlCjb2oqO3AZfg462oiNPbACIhRsBUcmqqRFhT9uFyFnQCVbXNV9njkpK0r/oc1Qh83TmtcJCZVte08cf5NcH8C7oqKfSlhfoDCfqcTpVVba/R2+FQ9Dp1095bwQz6LDbWuAw6BCRQUhISrIgKckCgAN2gsmXio7e94DjqaKdVXQcdYOrrATKy22/d9UmLWdd11q0sD/3Dh9W9pODjjwZAeBb0NGO0QFsKzomk+3F77FjwI4dtm3SPhejdl0zmeyrOikptu1X8zTobNwoXttmzSw480z//V1JTFQqct4EnXCdiEDNn+N0JEkJOh06OH6/BFpiItCnj9guLrb9cCGQjBJ0JEmCxWKBxWKBZJCBsww6RGQIMTESfvutCL/9VoSkJGP8gY0Ujio37lR0AP2pnbUVHfk+no7RAcQFrXqGM7k7m7M2aTmbjCA31/7cUwerQHVd0y4aKgedZs1sp1I+ftz5hWZMjH0oNUrQAeyDjqNua4BnU0zX1gK7d4vXdvXqIsTH+/fvitx9rbzc8wtmdUUnXIOOdpyOLw4eVN4voRqfI5O7r1kswIYNwXlMIwWdoqIiFBUVMegQEVFkCERFRxt09MaJuDNGB7D/xLypyb77mrsVHb0xOnrhSz0mSK/rWk2N7xUdea0SwLaik5oKZGQoPzt2zHnQ0QsbRhmjA9h3HXT2ab8nFZ2iIsdVQn/wZUICI1R0cnOV3/fy5d4vjAqEx4xrsmCP07FYlMdp3972fU++Y9AhIiKnPA062gtnvYqOo65rau4GHfVYFJm2+5ovY3RccTRGR11lcnfWteRkpVKWmKhcSGqDTsuWyn3UQUdvTIxe2AjkGB31+iPBruh4EnT8Ma7EGXnRUMDz7mtGCDomEzB8uNiuqBAz2HkrHGZck8lTTAPBCTrbtinv7XCu5hgVgw4RETkVjK5r7lQVgh10HI0D0dIuGAqIoKNug6PZv9T3AeyrF3L3tQMHlO5P2qCzbp0IQgBw7rn2FQ9XFR2jdV1zVtHxpOtaoINOpFd0ANtxOr50XwunoNOrl3LeFhQE/vGM0m3NqBh0iIjIKUfjZhoagClTgL/8xbablquua5JkH058mYxAPQOaTLvmjC+TEbgihxVtRUfdBkdBR1vR0V7Uy0GnqUl5/tqg89//KtsjRtjO+KVtl8xIXdc8CTqeVHQCvSCkL1NMy0EnNdVxNTAcqMfp+DIhgRx0TCbRJS6U4uKAAQPE9s6dwIkTgX08Bp3AYtAhIiKnHAWMDRuAl18G3n8fmDFDud1VRefmm4GPP7a9TS/oaDmq6OgFHU8rOk1NIljpjdFxJjlZmf5YO0bHnYqOu0FHu496jI46NI4YAXTqZLt/sLuuhXIygg4dlHPIWcVGkoJb0fEk6EiSEnQ6dnS96Gko9eqlnIvLl3s35sliUcbodOmiP24t2NTjdAJd1ZGDjskE5OUF9rGiEYMOERE5dPq0ewspfvSRsu0q6Kj3lckBQ13F2bMHmD5duXB2FHTKyuxv8zToyG3wdIyO+qJMW9FRt0Eex+Ps/oD7QUdd0ZFdcw1w3nnuBZ1I7bqWnq58Gl9YaDs7ntqePfrTkPtT27bK8/ck6Bw5ooxzCudua4AI+XJV5/hx2y5o7tqzR/ydAUK3fo6WepxOIINOba1SWTzzzPCu3hkVgw4RGYLFAvzwQxp++CHNp9l9yDOPPAL8+KPr/eTxOpLk3hgdLb1ZzwDghReAF18U245ed70FQj3tuibv46jrWlpaGgD7c08dYBwFndRUxyuda4OOdoyOvBaLWmoq0KaN7W133w188on4VNjTrmuRFHQA4KKLlO2ff9bfR764bGoCSkvT/nh9/ctkUqo6u3e7X+0wyvgcma/TTIfT+BxZsGZe27hRef8ZpdtaWlpg3i+BwqBDRIYgSWbcf39H3H9/R9TX809XICxbBjzxhO0g7jlz3LuvfCGvdzHnbjUF0A9F8+eLfz0JuN5UdBwFHbPZjI4dOyI93f7cUwcd7WQEcthydk3gbUWndWtg0iRRNXjlFeDVV5UKlKcVHaON0XHWdQ0ALrxQ2c7P199H7rZWX2/G6dMd0bFjR5gdpVEfyEGnthY4dMi9+xg56HgzTiccg0737sp5F8igY7TxOfLfwkC9XwLBGK0koqhnkLXJDKuiQsyg9OijwF//6vn95YqOXqDwpKKjd3/5tfek/7+3XdecjdHR637mTte1QAQdAJg3Tyy0eO+9tuM43KnoBGuMjt50157ytKJz7rnKc87P1//boZ6IoF8/n5rnlDfjdIwWdPr0UV6TX37xfJxOOK2hIzOblfEy+/e7H1I9pQ46Z58dmMeIdgw6RGQIDDqBtW2bsr1okef3lz/ccxV0HF0EuRN0gl3R0Y7R0XYrAxx3XTt1ynbdG0dcBR1HXddkegPVI32MjquKTmKiCDsAsG+f/vTjckUnJUUMgA8Ub9bSUQed7Gz/ticQYmKU9XSOHhULsXpCrujExopKSrgIxoQEctBJSBCBkfyPQYeIDMFstmDTps3YtGkzkpICuJx5lPI1SLobdByFFWdd17wJOt6M0dFORiCHAYvFgs2bN2P79s1o1sz23HMUdNTd/3yp6DRrZl/BcBacADELljqUadfqAYzVdU0bMF1VdADbcTra7msnT4oB8AAwaJAFRUWbsXnzZli8mTLMBW/W0jFaRQfwvvtaQwOwdavY7tEDiI/3b7t8EehxOidOANu3i+0BA8LruTsi/y0M1PslEBh0iMgQWNEJLF9/v866rp08CcycCfz7366DTjhVdPSml9ZedKuDijpQqNdw8STo6FWNtN3XXI0DNplsqzqRNr20O0FHPU5HO5nG778r24H+FF0ddPSmQdcjBx2TSb+iF468DTo7dyofboRLtzVZoIOOukpkhPE5RuXG5JlERKHHoBNYgQw6zz8PbNokttu1079/uAQdV+voaIOJo4qOelpjZ8EkLk6EDvlx9ao1HTrYDth2VdEBxDgdeexDtHVdA4D+/cUU3MeOiZnXGhuVcKdeP6dvX9/b58wZZyiv7/ffi9+Nq9/Jvn3i33btjPEpPyB+32lp4n23bJl4z7qz/k84TkQg69RJOYcKCtx/Tu4y2kQERsWKDhER+cxZ1zU55ADAwoX693fWdU3uIeFJTwntlNNyu5xdOLqzjo624uJo1jX12j6ugok6PDkKOmruBB1XFZ1ABh15DRhHj+0p9fM1m/WrXlpmM3DBBWK7ogJYt075WTCDTnIycO21Yru8XH8NKbW6OmXgu1G6rQHig47zzhPb5eVKdzRX1EEnXNbQkZlMSlXnyBHbLoX+8OuvyjaDTuAw6BCRIbCiE1j+qui4Gu/haAFHZxUdmScVHfXFtvq4rVs7vo87Xde0M6+pQ4p6EU/1eAxXXc0CEXTUM69F0vTS6emO1yTSctR9TZ5xzWwGevb0vX2uTJmibM+Z4/y9tn+/sm2koAN4130tnCs6QOC6r0mSUtFp3hzo2tV/xyZbDDpEZAgMOoEVyK5raq6Cji+TEahDiHamN7ka1KqV4/s7mozA0WNov1cfu7pa2fYk6LgzRsedoKOe7UtvWuxgjNExmfR/h55SP193uq3J9CYkaGhQLq67d7fvihgI55yjfGK/caPzEGDEiQhkvgSdxETb8UzhYuBAZdufM6/t2wccPiy2Bw/2b5c4ssWgQ0QUgSQJ+Oc/gcmTbS+6ne3vC2dd19S8qejIIcVV0FEHCnXQUR8zPd3xxbc7Y3ScBZ2EBP1Q42vXNe2AdHeCzpgxoltW+/bAddfZ/zwYY3QSE/1zAad+vu5MRCDLyVGmjl61Skz5vW2bcm4Ecv0cNZPJtqrz8suO9zVy0DnrLCWoL13q+m9Kba0y9XfPnsqHJeEkUBUdjs8JHgYdIjIEiwX45Zdm+OWXZh51YYpWn38O/OMfwKuvisDjSGWlmAo5WBWdkyf1b/dH17XkZGUMjqOgExfn+GLZ2RidZs2aoVmzZnZTNWsrAnpd40LRdS0lRYxF2bsX6NbN/ufB6Lrmj25r8nHkdXHUVRp3yPs3NIjFLNXjc/r3F//Kr20gXX89kJUltr/80vGaOkYOOrGxyut06BCwY4fz/bduVT7ECMdua4CYEEL+oKGgwPPFUB0xctAJxvvFnxh0iMgQLBYz7r47B3ffnYP6ev7pcuXrr5Xt+fNFkJk2DbjySuDgQXF7eblYkLB9e9spd73hbtBxxFnXNflnri4ykpLcCzrp6fr3dzRGx2w2IycnBzk5OUhKsj33tBUeva5x/gw6cXGicuQOk8nxp+TOuq5JkpgK/L333HscLX8HHUB0PSsoAJ56yrP7acfpaIOO+rU1uzv4xwtxccDdd4ttSRIfQOgxctABPOu+Fu7jc2Ry97XKStfhzV3qoKOuGoW7YL1f/MkYrSQiIo+oQ4HZDHz3HfDSSyIATZwobn/8cfGfd309cM89vj2e/H9eTY1393dW0Tl1Svwb6IqON5MRaL/3pqKTna0cKyPD/ufp6crjuFPNcYezrmvz5gF33gmMGwd8843nx5aDjruBzB2JiUBenvsTEchGjVK6z+XnKxMRAMHruib761+V8PfWW/ZToAPGDzojRyrbkRJ01EHEH+N0mpqU43TsCLRt6/sxyTEGHSIyBE5G4Bn178tstp1ed/Fi8W9Vlf8eLyZGVIjkKX091dgo2qzXjUoOIK6CjrsVHWdd11yN0XE2vTSgH3RchZOZM4Hx44F33tEfHG8yKV2CBgxwfix3OQo6TU3ACy8o37/7rufHDkRFx1sZGcon8ps2AWvWiO02bYJ/gdmqFTB2rNiuqtKfal0OOklJ+qE33OXlKe8JV+N0jBh0/DFOp7hY+fDGaN3WjMiroDN37lx07twZiYmJyMvLw/Llyx3u+/nnn+Oiiy5C69atkZaWhiFDhuD777/3usFEFJ1MJgt+/XULfv11C5KS/NRROoKpKzrHjgGPPmq/jz9mxJLFxAAPP+z9/bXjY7ROn/Zf0HG365rcvctisWDLli3YsmUL0tNtzz1/jNHp0gV48039SQNk774rgtAHHzg/lrscjdFZtMh2auxvvnFvMgu1cAo6gG33Nfm59O8vAqT6tbX4awCGE/fdp2y/8ortOS1JStDp2NGYM3HFxQFDh4rtAwccj0UClKCTlmY/Di2cqGde80fQMfL4nGC/X/zB46Dz8ccfY8qUKXjkkUewYcMGnHfeeRgzZgxKHayk9Msvv+Ciiy7C4sWLsW7dOowaNQpXXHEFNmzY4HPjiSi6JCdLSE5maccd7vwf5M+gYzbbLgzqKVdB59Qp97quyd2l/DEZgfr3I0kSJEmyCzb+6LrmjjZtRFeyNm18PxbgeIzO7Nm2+9XUAP/7n/vHtViU3324BB29CQzU3dbk1zYYevVS2lNSYvu7PX5cBHrAmN3WZO6M06msFBNlAOJ3Es6hLiNDmfp6wwbfJ+8wctABgvt+8QePg86LL76I8ePHY8KECcjNzcWcOXOQnZ2NefPm6e4/Z84cPPjggxg0aBC6deuGp556Ct26dcPX6pGyREQuGOjvalgIdtCxWHzrCtfQoD8RgezUKf9NRuDLGB1vuq75I+j4m9lsPyX4qlVK1y511euTT9w/rnqh1nAJOkOHwm62PHnGtVDQLiAqM/r4HJl6nM7Spfr7FBUp2+HcbU0md1+rqbFtuzfkoGM2i65+FFgeBZ36+nqsW7cOo0ePtrl99OjRWLVqlVvHsFgsqKqqQoaTzqd1dXWorKy0+SKi6Mag45lgB53GRs+7OKlNmWI7UFwrWF3XXI3R0VZ0tN9rZ12LjfXvoHx/kp+f/JzV1ZzZs5Xq0eLF7odYudsaED5BJyEBGD7c9rZgT0SgdsklYrFSAFiyRJnxcN8+ZR8jB51Bg5Rg6aiis2WLsm2koAP41n2tpkZ5vXv21F8gmPzLo6Bz9OhRNDU1ITMz0+b2zMxMHDp0yK1jzJ49G6dOncINN9zgcJ+nn34a6enp1q9seUoaIiJyS7CDTkODb0EHAK64wvHP3O26pg46cjhWhx5vu67JPK3opKWFb7ccuftaQ4PoRrVokfi+bVvgL38Ra78AIry4O/taOFZ0ANtxOgkJQI8eoWuL2SwW8pXJC4hGSkUnPh4YMkRsl5YCe/bY72OUiQhk/hqns2GD8nfMiN3WjMiryQhMmr/akiTZ3abnww8/xGOPPYaPP/4YbZx0NJ4xYwYqKiqsX/vUH3MQUVRiRccz7vy+5FDgD75WdADn93cn6KgrOoC4gF+0CPj0U+U2T7quqcexyDwdoxOO3dZkcpBraABef105Z+69V4QB9eeR7nZfC8eKDmAbdPr00X9tg+nWW5XK4vvvixkLIyXoAK7H6Rgt6Jx1lvKBhS9TTBt9fI4ReRR0WrVqhZiYGLvqTXl5uV2VR+vjjz/G+PHj8cknn+BC9V8cHQkJCUhLS7P5IqLoxqDjmWBXdNzpWuaLceNE1cEZdUUHAL76Crj2WuD555XbPKno6F0Mays62ot5owad774T22azWO8FEFNay1Mwf/ut/rovkgScPKl8H65Bp29f5eL7jjtC2xZAnEfyelZ1dWJR30gKOq7G6chBp3Vr/02wEUipqUBurtj+/XfbyqUnGHSCz6OgEx8fj7y8POTn59vcnp+fj6HyfII6PvzwQ9x222344IMPcNlll3nXUiKKemvXJmPt2mS3LuKjnTu/I392qdq/33/H0uNO72htRWf6dPt93B2jExtr+/tJTk7+40t9m/0Clikpthf4/lrgMxDkoFNerlx45uUp44xiYoA//Uls19WJxWa1/vIXoEUL4MUXxffhGnTMZuDnn8Vz/dvfbH8mv7bBds89yvkzdy6wc6fys3CebtkdgwcrY9O0FZ2jR5X3sxGqOTJ5nE5Dg/PxhM7IQScx0VjPXS1U7xdvedx1bdq0aXjzzTfx9ttvo7i4GFOnTkVpaSkmTZoEQHQ7GzdunHX/Dz/8EOPGjcPs2bNxzjnn4NChQzh06BAqKir89yyIKOI1NZlxxx1n4I47zkBdHdc6dsWdoOPrNKlq3n7C6St1VUobdPRmcXO365r6uGazGWeccQbOOOMMpKYq55622xogwpG6qhPOFR25YiUvXggA559vu4+6+5p2nE5NjbKuz9tvi3/VQSfcJmEwm+0rburX1qxNrQHWqRNwzTVi+9AhZVHfzMzwConeSEwEzjlHbO/ebVutMtpEBDL1OB1vuq8dO6ZUpc86y78V9WAJ5fvFWx638sYbb8ScOXMwa9Ys9O/fH7/88gsWL16MTp06AQDKysps1tSZP38+Ghsbcffdd6Ndu3bWr/vUq2YREbnArmueCXbQ0SMPsg4kdXVG23XN06Cj7rrmaAyHuuuaow81jRJ09C60LrjA9nt19xrtoPKDB5XtvXvFezRcKzrhSj3VtMzo3dZk6u5r6qqO0cbnyHydeU19H3ZbCx6vhuPddddduOuuu3R/tnDhQpvvlzqaRJ2IiLxSXS0GM8fHAwsW6F9QOluTBhDjaZwt0An4XqUJxoSZ6emiKwxgX9HRa7+rrmt6FR01dRVHr6ID2E4xbYSua7L4eGDYMPvbWrcGjhyxDTaAbXfF6mqgooJBx1PDhonugnI1B4icoKOdkGDsWLGtDjq9egW3Tb7o1098ANLY6F3Q4fic0DBG3YmIop7JZMGyZcVYtqwYSUnRPUhn5kzg88+Bjz6yX8lepr7g1NOjhxIQ9EiSWGjRFw0NgZ/dSl0xcTfoOAofjY1KlUsdAiwWC4qLi1FcXIy4OAvkZeDatdM/jlEqOtrXZsgQ/SpV+/bi34MHbSuFBw7Y7rd3r/GCjvq1tYRg8J/JZF/ViZQVNc45R3k/qis66q5rRgo6iYlixj4AKC72fJbJSAg6oX6/eINBh4gMQZKAjIwmZGQEcGovD4WqO508QxagXEBUVQHvvqsMaHYVdEpKRDXIkYMHgfXrfWtnTY3vQefGG4E//9nxz9VBIjnZdlyIXlUrPt5+AgGZs4pOU1MTmpqaYDKJGbKuuAJ48kn94xgl6Gifo3Z8jiwrS/zb2GgbjrVBp7TUeEEHUF7bULnhBmV2OyByKjpJScoF/c6d4nyRJKWi06GD426k4UruvmaxiDVx3CVJStDJyADOOMP/bQuWUL9fPMWgQ0SGEG4fHs2aJS5o33kn1C0RHn1UdGcbNUpckLoKOq744/8xd4NOs2b6a20MHw7Mmwd8+KHjNX+cVXT0OBsA7M4YHUAspPnVV8DZZ+v/3KhBRzs+RyZXdADbcBMJFZ1wEB8P3H238n3PnqFri79px+mUlQEnTojvjTQ+R+btOJ29e0X3T0CEv3BdRDgSMegQkSGE2wdIM2eKWXRuuy34j62uJMn/YcoD//fvFzM4+Rp0/BUs3Qk6R4+KUKOutIwdKy6MWrQQ3zsaD6Ot6LgbdN58U1zA33yz8jN3xui4Q931KJzXCFE/x5QU24s4NbmiA9iO09FOKW7Uik44ePBB4KGHgMceAy66KNSt8R/tOB2jTkQg8zboREK3NaNi0CEiQwinoBPqGeD0go5adbXriQZc8XUighYtxBor7gQdubuZekzQ6dO2+zgKOnKf+cREICfH/aAzfry4UJ86VfmZeh0dX4LOddcBV18NXHWVMn1wOFK/NsOHO/7dqYOOs4oOg4734uOBZ54RH6AYZNZetwwZopxnkRB0evZUzmtPpphm0AmdCHo7EVEkC/RUyJ5wNaNZoGmDjrZ68+WXtutWeMOXoPPaa6JPfrNmno3ReestccFnMolPuNX0gs7YscC994rFFvPzgZYtPe+6pm6fuuuaL0EnORlYtAj44gvHU1CHA/VzdDQ+B7Dtuqau6Oh1XVOfNww6lJKiXNhv2ybepzIjBp24OGDAALG9c6fSDc8VddBxVDmlwGDQISJDCKeKjq/dwnylDTraLkQPP+z7Y/gSdM46C9aZyTwJOt27A9u3A1u32n/qqQ46nToBn34qAk5ysljp/txzxc88DTrq79Vd1wI9W1w4UE/c4Czo6FV0LBYx3kKNFR3So+6+9v334l+TCcjNDU17fKUOKu5UdRoblenDc3LCuztrJGLQISJDaGwENm9OwubNSSGfmMDXbl2+0gYd7fom/uBLmFNXMTytjHTqJAKPs2O2bCkmBFAv3inzpaLT0KAEau1+SUlJSEpKcn5wg7nuOtFNasQIoH9/x/vpVXTKy+2rrGVlQGWl8r1Rgk4kvrbhRB105L9dZ5wR3tVOZwYOVLbdCTpFRUpX3Ejotma090sUfGZFRJGgvt6Mm27qEupmAAi/oBOICpN2jIwn1NUXV5WRRx9175jq/1ed/f59qejU1Ojfbjab0aVLeJx7/nT77cCVV4rqm7NZoFq1Er+PhgaloqOtIgLivJSnNwdsK0bhKlJf23AydCgQE2NblTditzWZpxMSRNL4HCO+X1jRISJDCKcxOuHWdS0Qv5tgBJ2JE8UsU+5QVwec/f59qeg4CjqRrGVL11Pdms3K4qhyRUc9Pkc9eH77dmXbKBUdCqzUVNsqCGDsoNO9uzLbY7QFHSNi0CEiQwinMTrhVtHxdYY1PeqLfk+pu6Q4Czrnnef+DFPqi2ZnbfM06KiPW1GhbEfDGB1PyON0jhwR57866PTtq2yXlCjbDDokU3dfA4wddMxmIC9PbMvT+Tvz66/i35gYMX6RgotBh4gMwoLvvtuG777bhsTE0A7SCbeKjidBx90LeF8qOu4GnZgY94/pbkXHVXcpbdCR1+kBxLgTvf0sFgu2bduGbdu2wRLqAWIhoh6nc+iQbdBRTwuuPheNEHT42gZHJAUdwP1xOqdOKVNq9+7teJp8ozDi+4VBh4gMoakJaN++Ae3bN4R8VelAV3TmzQMuvhjYsEH/574EHfUCm874EnTUISHYQcebMTqpqWL78GHH+zU0NKAhEKUzg9DOvKYeo6MOOmpGCDoAX9tgOPdcpXobG6s/4YiRuDtOZ/16ZfHlSOm2ZrT3C4MOERlCOHVdC2RFp6ICuOsu4Icf9C8gDx8G9uxRvvc06KSnu7efL0FHLRBBx5+TEQDKVNjqrmvRMkbHXdqZ19QVnSFD9O9jlKBDgZeWBlx6qdi+8ELX79Nw527Q4fic0GPQISJDCKfJCAJZ0Tl5UtnWC1RTpth+7yroaMfAuFvRcTVG529/Azp3dn0cZ2HG26DjLPR6E3RatrS/jWN0bGkrOnLQSUkR54FelxwGHVL74ANg8WLgk09C3RLfdeqk/N0oKLCtsqsx6IQegw4RGUK0VHRcdctbv95+f2dBRztmxV9d1265Bdi1S38tG237HPE26DjjS0XH1X7RzFFFp3178Rp36mR/HyNML03Bk5oKjBmjdBU1MpNJqeocOSIWy9UjB53kZKBnz+C0jWwx6BCRIYRT0An1rGtqngYdf3Vdk0OK9rE9WQTQk6oJg05oqSs6W7cCVVViWw5AHTva7h8f7/6MekRG5Kr7Wnm50s04L49V4lDhnyEiMoRwCjraio63k88sWSLG4fz738ptrrroabtIuAo62oDgr4qOo6CzcqV7x1cfwx2BDDp6XdcYdGypKzrq7jiOgg67rVGkcxV01Lex21roMF8SkSE0NgI7d4ryhKP+0IGyYwdw4oT4j81ksq/oNDa6N7i2shJ48kkgJ0eMcTn/fHH76tVi8Ux3JhZQD5gHgL17xcQFjmgrOu5Ob+pqjI4cUtQhb9o0oH9/946vPoY7gl3R0X76mhDl/bBSU0U3xepq23VDOnQQ/2q7rhkp6ET7a0vecTXFdKSOzzHa+4VBh4gMob7ejGuu6Rb0x92xQ/StbmwEvvtOTPusreioq00lJcDVV4tPuL/6yvZi/tFHgVdeEdt9+tgfo6ICGDDA9vZ33xXjIaZMEfdXr/UCuF6ZW3vB6e4FqLsVHTVPZ1IK56Cj3s9sNqNbt+Cfe+GmfXtg2zb72wD7io5RroX42pK32rUT5/+BAyLoWCy23TUjMegY8f3CrmtEZAjedF07eBA4dsz7x1yxQqz3IHcnGztW/Lt9u+1+6u5mN90kFohbvBj48EPb/eSQA4jQpFZfD9x/v30l5dZbgb//HXjjDWD2bM+fg/aC019BR6+/uV6IMMJkBHpd19wdyxRN1ON0ZHLQMXJFh8hbcve1ykrxoZhMkpSg07q1/mQdFBwMOkRkCJ6Og9m4UXzK3KGD7ZofnrjgAtvv5S5z8+bZ3q4OOuoKy969jo+t7X7X0ACsWuV4/0cecfwzZ7QXnO5+0u5uRWfcOOU27ernrngyOLdVK/f281dFp2tX9x4vmqjH6Whv4xgdikaOuq/t2gUcPy62Bw92PZsmBQ67rhGRIcTEWLBoUQkA4KabusDV5zTjxokqUFOTqIgsXCiCRPPmQK9e+vf59lvg11/FOJaLLhJVFrXERP3KkqNqk3phTy1t0Kmvdz5tdZs2YnyEpwJV0ZGDzuzZIjyccYZ9MAT8V9G58EJxUfH778Dnnzvez1WQczfoqHtnWCwWlJSIc69Lly4wR+l0YnoVHXmMTvv2otuO/IGEUYIOX1vyhXZCgltuEduR2G0NMOb7hUGHiMLaSy+JikxlJdC1q5gFwGQCRo4EZswQY2b0qBferK4G8vOVfYcNE+FnzhyxBsLXX4ugIa/c7Uhiov7U0nJF5+uvbW9/800xo5rexb42RC1YAOzb5/ixMzLEp4Se8naMjruTEbRqJZ6nNzwJOmazCKFVVc67lbmq6Oj9v6zXdU3bDb0unOYUDxFtRScmBsjMFNuxseLn8jlslKAD8LUl76krOupqfqQGHcB47xcGHSIKWz/+KGbyAoCkJNufLVsmvrSVkcZG/S5RN92kbK9cCVx2mfL9+PHAlVe6bk9ion7VpalJzJymd4zKSv0Lc/XMVQDw0EPOH1v7/N0V6IqOLzw9htnseuyMpxMiAPYVnZYtgRYtPD9OpNNWdNq2tX0NO3UyZtAh8lZGhqhm79oFbNig/P+jDjrqqg8FX/jXnIjILfv3AzNn2v6BNbpvv3Vvv8ZG4OhRYOpU0TXJZLJfqdpZF6qlS90bA7R5s/6n/42NwFNP6d9H7qet9f77rh9PzZtua4D+GJ1//EP53tECn+EYdNzhTdDRhhqDTSoUNNqgI3dbk6nH6TDoULSQg0xNDVBUJMZbrl8vbuvSRf//DAoeBh2iCHHVVcCsWcDZZ3u/gGW4cWcAZ2Mj0K+fmNlmzhz9fRYtsu3KpncMX9bmuf9+YOdO/Z/Js75pKzie2rDBvf20Y1D0KjqPPy66A1osYkrriy6yP86pU84fx92Q4q8xOu7yZvVx7e+MQUeftuua9nsGHYpG2nE6mzcrlf9I67ZmRAw6RH6wbJkIGosXh64N8idIgO0sYOHk99+BTz91vSimzJ2g88034lM0ZyTJ+fTUFot9BcgTn30GbN2q/7NBg4DXXxfVtmDQzginDTryxWhWlvj9xsbqX5S66obtTaAIxDG0TCbvqjpqDDr62rWz/V4bdNRT6DLoULTQBp1IHp9jRAw6RH4wcqRYHFI97iOUvFlzJtCOHROVlxtuAF5+2X/HddQ1zFPerFHjrr/9TVROgsHV5AM9etjfx5upT72pxmgnAghERQfwPeh06eKfdkSa+HhROZWxokMkFnmW/4YWFDDohBsGHaIIoJ3GOByDzpdfKtvTp7t3H/UFuCQBBw7E4cCBOJtuZu5Wh0Lt44+D8zjaSQu0FZ02bezv46wC6GhQvjdd17QVnHANOvJMYmpxcXGI05ubOsqox+lox+j066e8xkaqivG1JV+kpgK5uWL799+B5cvFdkyMCEGRxmjvF866RhQBRo60/T4cg4720/x//lOsZ3Pttfb7SpK4QFZfJNfWmnHJJfblCO00zdFOG3S0n6zrVW+chcUWLYATJ+xv9yakxMbavl7hFHQefBB47jkxcPicc2x/Zjab0UOvFBaF2rcXi/HK29qfff216Er6178Gv23e4GtL/jBokDIRwY4d4ra+fb2fLTNcGfH9wooOUQTYu9f2+3AJOvv2KRfR2ovaf/wDuO465T8F2aRJYm2Wr75y7zGMUtHxp2uucfwzbbCJiwOef16ESkez2LkKOnq8qehoPwQMp6Dz6KNiPaCffxYLxpI+ebHd2Fj9bpCXXCKmhG/WLLjtIgol9Xo6MnZbCw8MOkQRKByCzrvvij77eXliauTycv39Hn1UmfGsvByYP1+Mu7nqKjEzmCuVlf5rs1E4G1Oj/QQxNhZ44AExE9All+jfJ5Bd19SCFXS8GWPTrJlYT6lvX/+3J5JMny7Op/fft5+cgCha6a2Vw6ATHhh0iCLQxImBOa6ztVXmzQOmTFEmB7j1VvHvpk2iD/MDD+jf7+OPRbe2O++0XytGvdZMQoIFH35Ygg8/LEFCgjJ/9uOPe/48jM7Z66AXdFxxVtHRLqbpyXG1tEEnELOuAcDw4f49nsViQUlJCUpKSmCJlLnbvdS6tagQ3nBDqFviH3xtyR/U49NkkRh0jPh+YdAhikBffqn0o/eXa68Vn+5/+KH9z1avBu66S8ymNmuWqOZ46t//BkpKHP/cbAZ6965B7941duN9ok3Pno5/pu265mvQcVTR8eY1CNZkBCNG+P+YNTU1qKmp8f+BKeT42pKvEhNtq8EpKcoEBZHGaO+XKL9cIIpchw/7foz9+0VXlbfeEotu1tcDN99sv596PM3LLyvVHE85W9SThPbtgf/7P8c/145PcSdMaGfPUmve3K1muSVYQefsswNzXCIiR9TjdAYODNzfN/IMgw5FBUkCFi4EXn01+gavL1smuncdOeL5fa+/HnjhBWDCBOf7abucecsglfCQ2rHDvsrSqRNw5ZWiG6A33cPmzAHS0vR/plfRufRSt5oKIDTTSycnB+a4RESOqMfpRGK3NaNi0KGokJ8P3H47MHmy6CIVSZx196qsFFNPP/aYWD+ltlbcvnw58PDDQGmp82P/+qt7bfBXV7KHH/bPcSKZ3kKMf/mL6K54ww32YcKdoNO5s5j44Ykn7H+Wnm77/ZYttmsiuRKKoAMAb78duGMTEWldey2QkyM+HLrjjlC3hmQMOhQV3npL2X7qqdC1wx+OHwd+/BH47DNg1y7gzDP193v+eWDnTtvb5s8XFa3hw4FnnxXTO3ujb1/g00+V73ft8u44WtqFT41qyhTgoouA88/3/7H1ZlxThx9teHB3wH+zZsDFF9vfrl1gtGdPzyYRUI8n0i6eF8igc/vtwO7dwOefB+4xiIhkGRni/9xDhxz/v0zBx6BDZCDHjwPdu4uL6OuvF9PoOpoa+McfxdTOalOmAGPHKt8XFAA1NcDBg0BZmaj0uNN9bNMmUT144QVgyRLgm2+8fkphT299BFeeew744QcxE48znToBmZnetUutrk7Z9qVqMnAgMHeuWMcoM1NM8X3ZZb61bdYsoH9/sf7KSy953zZv5ORwTRwiCp6YGO/W8aLACdDknkQUCK+/Dhw75tsxPv7Y9nvteIZ//lN0h3Nn8ePp031ri6eOHw/clXFysu20ze+8A5xxBnDOOfbjXlyR909IcL5f164iEPXpI1bV9kRKCnDqlNjevVu53Zuua2p/+5v48pf0dGD9erGtrUYFY/Y8bdc7b8VwZHHE4mtL5D6jvV9Y0SEAoirw22/OFw701okTwKhRoluMs/U/jEiSlHEvsunTbWchc9fJk+ICuW9fUV2R7d8PzJ4tBqFXVfnUXLc8+qiYuGHGjMA/lidqaswYMSIXI0bkoqbG/3+6nn/e9vvUVODcc31b60VvPI1a27biYn/1as8Hr15zjbKtXiDT16Cj59NPgdGjRfXOGyaT80VOA2nwYOCCC8T6Qt5WHs1mM3Jzc5GbmwtztM9tHmH42hK5z4jvF2O0kgJu3DgxJettt/n/2H//O7B0qfjkWm+wszO//iou4vzZLkny33EuucR+gcYXXhBdfuSLu7173TteixYiaG7aJLoLvfQS8MorQHa2WGyze3fgmWf803ayV19vGza6dnV9n1deAVq2tL3t6aeV7WbNnN9fXlk+Lc12xh53wskrrwBnnSW6hd13n3K79sM2f1Q0rr8e+P57MbGF0ZhMYjKSo0d974ZHRETGwq5rBjd1KrB4sRhsf+653h9HXgTy/feB//zH+b4VFeIT6JEjXX9iDQArVyrb7s7iJRsxQow/2LULmDRJdCPyVE0NsGKF5/dzRJLEmjLuDuTPyRFr2rRuLULMl18Cw4aJT8k7dhRdxPr3t73Phg3ii9yXmGhfXVN379IzcCCwbp2oqlx2maiSTJ0qQkefPq4fc/Bg8drGxIhq29attuOiXK1Bk52tbDc1KduxsWLCiJ9/tt1f/QFaixai7VrRsDq3p0wmTjlNRBSNGHQMbNcusf4FAJx3nv8qFa5cfrkIDuPHA2++6Xr/o0eVbU/XsFEPsn7qKTEI/5JLgG7dHN9HkpRuMlVV9uuDlJWJgdH5+eJicd48EULGjhUXyomJ4hgnT4qB+W+8oXTjSkiwbZO7/DHgPNolJFgwb94eAMDf/paDujrlqv+ZZ8TaLn/9K7BmjXKfSZNEtz9Hvv5aBJWEBOWc0pula84cMZFDVhYwcaJYl6hdOzGLmFxBSU21rcoAjtemkXXvrmyr3xuxsWJ65KefBj74QOmyqK0e6tH2JtDOmkaes1gs2PPHlIA5OTmG6bJBrvG1JXKfEd8vXrVw7ty56Ny5MxITE5GXl4fly5c73X/ZsmXIy8tDYmIizjjjDLz++uteNTYa1deLmbD0Lq6PH3d9/9OnxaxHvqwpoQ5QkqRUR9RTNsvq6oDCQtv7qMebLF8OzJwpLhw3bhThZfx4sf/+/cCLL4oKx/PP2/fp//prsQ5O9+4ilKSni8AidxGTv8xmZdvRhWZRkVhhPjkZuP9+0XXPZBIXkvIxMjLE7FPqsSrehBzyD7MZGDToNAYNOm1zMT9/PvDQQ6ICo36tzjlHdJXs1cv+WH/5i6j0tW0rZkZzNRXovfeKc7ewEHjkEXEurlnjenYdV1UEdWAfPVrZnjBBzMj2+uvATTcpt7sbWv7xDxGu5Uot+e706dM4HWmDDAkAX1siTxju/SJ56KOPPpLi4uKkN954QyoqKpLuu+8+KSUlRdq7d6/u/rt27ZKSk5Ol++67TyoqKpLeeOMNKS4uTvrvf//r9mNWVFRIAKSKigpPm2sI338vSeJSn1/84hcgSX362N+WlNQkbdq0Sdq0aZOUlNRkvX35cuW9ZLFI0h13SNKwYZK0a5dy+08/2R7rxIngvLd//ln/+XXoIEn33GO7r8UiSX//uyTdeackVVYqt+/ZI0mxseJ+q1YFp92BNGGCeC7XXhvqlrivqUk595qamkLdHPIjvrZE7gun94u72cAkSZLkSTA6++yzcdZZZ2HevHnW23Jzc3H11VfjafUo3D889NBD+Oqrr1BcXGy9bdKkSdi4cSNWr17t1mNWVlYiPT0dFRUVSHPVF8SAQjUbEUWXmBjg99/FmC5/TAvdr5+oyjnSqZNYn0fbXXHuXDGz3AcfiG09o0bZzvA1dCiwYYMFv/0m5mD++uue+PvfzejfX6wF5M5sl/36iecPiEqpp1NGe2PNGmDIEPvbPfurC5SWigqUO1N+h7umJjH+rF+/4LwG/mCxWFD0x/zfPXv2NER3DXIPX1si94XT+8XdbOBRC+vr67Fu3TqMVvexADB69GisWrVK9z6rV6+22//iiy9GQUEBGhwM2Kirq0NlZaXNFxE517y5fRfFNm3ErG+SJGZ069lTzOD23Xeie9Q334jJETx1661iQdLCQuW2wkIx8P/rr0X3xj17xAX6rl1iXJfsiivEZAz/+pdo0/bt4iJenoHs0UdFl624ODF2ZtcuMaGFesruhx4CiovFpBjuTun/1VfiuS9ZErwLbO2EBsnJ3k1xLE9aEQliYsQkEEYJOUREZFweTUZw9OhRNDU1IVMzsjozMxOHDh3Svc+hQ4d0929sbMTRo0fRTr66UXn66afx+OOPe9I0IpfatwcOHPD/cXNygAULRKXko4/EBAYHDogppu+8UwyS1+rWTQyCv/BCYOdOcTFfXS2mC77nHnHR/8wzwLJlYhzJhx+KQJCTI6od27cDublinFbLlmJ8U2amGDNy++3iMeSqgV7F8OKLxZe839q1YkKGTp2UfQ4eFGOXTp8W46oWLhRjnp58UoxdAsS/J06ISRzathW3qUONfNv8+eL5nHce0KGD8vOYGGWcysqV4pP+yy4TAWf/fjEuR36s88+3XVTT1bgarU6d7NfKCbSUFDEz2qpVYlxZQoJ7MxUSERGR77yadc2kuXKSJMnuNlf7690umzFjBqZNm2b9vrKyEtnqeVgjjNxzv6lJGQjf0CC+j40VP4uLEwPhY2OVqWzlT0Tl+8TGKvc1m8WXJIkLb0kS3zc1KccymcSFscWiPFZjo9hP/Sl5dbX4JLqxUexvMonjqV++mhrlZyaTOJ58DPX+DQ3icbTHkNvozqfz2sfWft/Q4N9Pi9VtdmbkSOC55+xvl19fT7soXnKJ/fFl8gB7OQTovT3cfTyTSXzCrpWVJao2MjkYabmaQlk+1iuvON+nc2fxJYuU2cLOOkt8ERERUXB5FHRatWqFmJgYu+pNeXm5XdVG1rZtW939Y2Nj0VK70t4fEhISkJCQ4EnTDM9ksl3/Qm82J/Unwc6GKqkv8uUQJJODhPpYMTHK7XqPm5pqe1/5uGraaW+1a3nI+2vbJvOkm6f2sbXf+7tLjD+Ox3FY/uHsAxWiQOK5F7n42hK5z2jvF4+CTnx8PPLy8pCfn49rrrnGent+fj6uuuoq3fsMGTIEX3/9tc1tP/zwAwYOHIg4dtImIjeZzWb00psrmijAeO5FLr62RO4z4vvF4+kSpk2bhjfffBNvv/02iouLMXXqVJSWlmLSpEkARLezcePGWfefNGkS9u7di2nTpqG4uBhvv/023nrrLTzwwAP+exZEREREREQqHo/RufHGG3Hs2DHMmjULZWVl6N27NxYvXoxOf4xiLisrQ2lpqXX/zp07Y/HixZg6dSr+9a9/ISsrC6+88gquu+46/z0LIiIiIiIiFY/X0QmFSF9Hh4hcs1gs1g9ROnbsyPUuKGh47kUuvrZE7gun94u72cCrWdeIiEKhuro61E2gKMVzL3LxtSVyn9HeL/zogoiIiIiIIg6DDhERERERRRwGHSIiIiIiijgMOkREREREFHEYdIiIiIiIKOIYYtY1eQbsysrKELeEiELFYrFYZ3uprKzkNLAUNDz3IhdfWyL3hdP7Rc4ErlbJMUTQOXbsGAAgOzs7xC0hIiIiIqJwUFVVhfT0dIc/N0TQycjIAACUlpY6fTJGNmjQIKxduzai2+Cv4/tyHG/u68l93NnX1T6VlZXIzs7Gvn37InaBXJ7vwTlOIM93f+3H8z0y2mDE893ff9vd2Y/ne2S0gee76/0Cfa5LkoSqqipkZWU53c8QQUcujaWnp0fsH4aYmJiQP7dAt8Ffx/flON7c15P7uLOvu8dLS0sL+TkRKDzfg3OcQJ7v/t6P57ux22DE893ff9s92Y/nu7HbwPPd/f0Cea67U/xgZ9Qwcffdd4e6CQFvg7+O78txvLmvJ/dxZ99weK1DLRx+BzzffbuPv/eLZOHwO+D57tv+PN/dFw6/A57vvu0fSee7SXI1iicMVFZWIj09HRUVFSH/lIAo0Hi+UzTh+U7RhOc7RYtwOdcNUdFJSEjAzJkzkZCQEOqmEAUcz3eKJjzfKZrwfKdoES7nuiEqOkRERERERJ4wREWHiIiIiIjIEww6REREREQUcRh0iIiIiIgo4jDoEBERERFRxAnroPPLL7/giiuuQFZWFkwmE7744otQN4nIL1yd25Ik4bHHHkNWVhaSkpIwcuRIbNmyJTSNJfKQP87vuro63HvvvWjVqhVSUlJw5ZVXYv/+/UF8FkT2gnVunzhxAmPHjkV6ejrS09MxduxYnDx5MsDPjqJdOJ3fpaWluOKKK5CSkoJWrVph8uTJqK+v9/g5hXXQOXXqFPr164fXXnst1E0h8itX5/Zzzz2HF198Ea+99hrWrl2Ltm3b4qKLLkJVVVWQW0rkOX+c31OmTMGiRYvw0UcfYcWKFaiursbll1+OpqamYD0NIjvBOrdvvvlmFBYW4rvvvsN3332HwsJCjB07NuDPj6JbuJzfTU1NuOyyy3Dq1CmsWLECH330ET777DPcf//9nj8pySAASIsWLQp1M4j8TntuWywWqW3bttIzzzxjva22tlZKT0+XXn/99RC0kMh73pzfJ0+elOLi4qSPPvrIus+BAwcks9ksfffdd0FrO5EzgTq3i4qKJADSmjVrrPusXr1aAiBt3bo1wM+KSAjl+b148WLJbDZLBw4csO7z4YcfSgkJCVJFRYVHzyOsKzpE0Wj37t04dOgQRo8ebb0tISEBI0aMwKpVq0LYMiLfuXN+r1u3Dg0NDTb7ZGVloXfv3nwPUNjy17m9evVqpKen4+yzz7buc8455yA9PZ3nP4VMMM/v1atXo3fv3sjKyrLuc/HFF6Ourg7r1q3zqN0MOkRh5tChQwCAzMxMm9szMzOtPyMyKnfO70OHDiE+Ph4tWrRwuA9RuPHXuX3o0CG0adPG7vht2rTh+U8hE8zz+9ChQ3aP06JFC8THx3v8HmDQIQpTJpPJ5ntJkuxuIzIqb85vvgfICPxxbuvtz/OfwkGwzm9/vQcYdIjCTNu2bQHA7lOL8vJyu084iIzGnfO7bdu2qK+vx4kTJxzuQxRu/HVut23bFocPH7Y7/pEjR3j+U8gE8/xu27at3eOcOHECDQ0NHr8HGHSIwkznzp3Rtm1b5OfnW2+rr6/HsmXLMHTo0BC2jMh37pzfeXl5iIuLs9mnrKwMmzdv5nuAwpa/zu0hQ4agoqICv/32m3WfX3/9FRUVFTz/KWSCeX4PGTIEmzdvRllZmXWfH374AQkJCcjLy/Oo3bGeP9Xgqa6uxs6dO63f7969G4WFhcjIyEDHjh1D2DIi37g6t6dMmYKnnnoK3bp1Q7du3fDUU08hOTkZN998cwhbTeQeX8/v9PR0jB8/Hvfffz9atmyJjIwMPPDAA+jTpw8uvPDCUD0toqCc27m5ubjkkkswceJEzJ8/HwDw17/+FZdffjl69OgR/CdNUSNczu/Ro0ejZ8+eGDt2LJ5//nkcP34cDzzwACZOnIi0tDTPnpRHc7QF2ZIlSyQAdl+33nprqJtG5BNX57bFYpFmzpwptW3bVkpISJCGDx8ubdq0KbSNJnKTP87vmpoa6Z577pEyMjKkpKQk6fLLL5dKS0tD8GyIFME6t48dOybdcsstUmpqqpSamirdcsst0okTJ4L0LClahdP5vXfvXumyyy6TkpKSpIyMDOmee+6RamtrPX5OJkmSJM+iERERERERUXjjGB0iIiIiIoo4DDpERERERBRxGHSIiIiIiCjiMOgQEREREVHEYdAhIiIiIqKIw6BDREREREQRh0GHiIiIiIgiDoMOEVEEW7p0KUwmE06ePAkAWLhwIZo3bx7SNvnKX8/BZDLhiy++cLrPsWPH0KZNG+zZsweA/e8znHzzzTcYMGAALBZLqJtCRBQWGHSIiKLIjTfeiO3bt4e6GT4J5nN4+umnccUVVyAnJycoj+eLyy+/HCaTCR988EGom0JEFBYYdIiIokhSUhLatGkT6mZ4raGhIWjPoaamBm+99RYmTJgQ8MdyRZIkNDY2utzv9ttvx6uvvhqEFhERhT8GHSKiELNYLHj22WfRtWtXJCQkoGPHjnjyyScBAOeffz7uuecem/2PHTuGhIQE/PzzzwCAuro6PPjgg8jOzkZCQgK6deuGt956S/extN2+HnvsMfTv3x/vvfcecnJykJ6ejj//+c+oqqqy7lNVVYVbbrkFKSkpaNeuHV566SWMHDkSU6ZMcfic5OPOnz8f2dnZSE5Oxp/+9Ce7Ll8LFixAbm4uEhMTceaZZ2Lu3LnWn+3ZswcmkwmffPIJRo4cicTERPznP//R7bo2b948dOnSBfHx8ejRowfee+89m5/v2LEDw4cPR2JiInr27In8/HyHbZd9++23iI2NxZAhQ+x+tm7dOgwcOBDJyckYOnQotm3b5nZ75OdVWFhove3kyZMwmUxYunQpAKWL3Pfff4+BAwciISEBy5cvx8aNGzFq1CikpqYiLS0NeXl5KCgosB7nyiuvxG+//YZdu3a5fH5ERJGOQYeIKMRmzJiBZ599Fo8++iiKiorwwQcfIDMzEwAwYcIEfPDBB6irq7Pu//777yMrKwujRo0CAIwbNw4fffQRXnnlFRQXF+P1119Hs2bN3H78kpISfPHFF/jmm2/wzTffYNmyZXjmmWesP582bRpWrlyJr776Cvn5+Vi+fDnWr1/v8rg7d+7EJ598gq+//hrfffcdCgsLcffdd1t//sYbb+CRRx7Bk08+ieLiYjz11FN49NFH8c4779gc56GHHsLkyZNRXFyMiy++2O5xFi1ahPvuuw/3338/Nm/ejDvvvBO33347lixZAkAEyWuvvRYxMTFYs2YNXn/9dTz00EMu2//LL79g4MCBuj975JFHMHv2bBQUFCA2NhZ33HGH2+3xxIMPPoinn34axcXF6Nu3L2655RZ06NABa9euxbp16/Dwww8jLi7Oun+nTp3Qpk0bLF++3OPHIiKKOBIREYVMZWWllJCQIL3xxhu6P6+trZUyMjKkjz/+2Hpb//79pccee0ySJEnatm2bBEDKz8/Xvf+SJUskANKJEyckSZKkBQsWSOnp6dafz5w5U0pOTpYqKyutt02fPl06++yzre2Li4uTPv30U+vPT548KSUnJ0v33Xefw+c1c+ZMKSYmRtq3b5/1tm+//VYym81SWVmZJEmSlJ2dLX3wwQc29/vnP/8pDRkyRJIkSdq9e7cEQJozZ47NPtrnMHToUGnixIk2+/zpT3+SLr30UkmSJOn777/XbQsAadGiRQ6fw1VXXSXdcccdNrfJv88ff/zRetv//vc/CYBUU1PjVnvk57Vhwwbrz0+cOCEBkJYsWWLzOF988YXNcVJTU6WFCxc6bLMkSdKAAQOs5wcRUTRjRYeIKISKi4tRV1eHCy64QPfnCQkJ+Mtf/oK3334bAFBYWIiNGzfitttus34fExODESNGeN2GnJwcpKamWr9v164dysvLAQC7du1CQ0MDBg8ebP15eno6evTo4fK4HTt2RIcOHazfDxkyBBaLBdu2bcORI0ewb98+jB8/Hs2aNbN+PfHEEygpKbE5jqOqiqy4uBjDhg2zuW3YsGEoLi62/lyvLa7U1NQgMTFR92d9+/a1brdr1w4ArL8zV+3xhPa5T5s2DRMmTMCFF16IZ555xu53BYhxWKdPn/b4sYiIIg2DDhFRCCUlJbncZ8KECcjPz8f+/fvx9ttv44ILLkCnTp3cvr8r6q5PgJh2WZ6iWJIk621q8u2ekI+hPv4bb7yBwsJC69fmzZuxZs0am/ulpKS4fWx1++Tb9Nqq3V9Pq1atcOLECd2fqX9n8rHU0zo7a4/ZbLZrV0NDg+7jaJ/7Y489hi1btuCyyy7Dzz//jJ49e2LRokU2+xw/fhytW7d2+tyIiKIBgw4RUQh169YNSUlJ+Omnnxzu06dPHwwcOBBvvPEGPvjgA5vxIH369IHFYsGyZcsC0r4uXbogLi4Ov/32m/W2yspK7Nixw+V9S0tLcfDgQev3q1evhtlsRvfu3ZGZmYn27dtj165d6Nq1q81X586dPWpjbm4uVqxYYXPbqlWrkJubCwDo2bOnbltcGTBgAIqKijxqizvtkUNIWVmZ9efqiQlc6d69O6ZOnYoffvgB1157LRYsWGD9WW1tLUpKSjBgwACP201EFGliQ90AIqJolpiYiIceeggPPvgg4uPjMWzYMBw5cgRbtmzB+PHjrftNmDAB99xzD5KTk3HNNddYb8/JycGtt96KO+64A6+88gr69euHvXv3ory8HDfccIPP7UtNTcWtt96K6dOnIyMjA23atMHMmTNhNptdVkUSExNx66234oUXXkBlZSUmT56MG264AW3btgUgqhOTJ09GWloaxowZg7q6OhQUFODEiROYNm2a222cPn06brjhBpx11lm44IIL8PXXX+Pzzz/Hjz/+CAC48MIL0aNHD4wbNw6zZ89GZWUlHnnkEZfHvfjiizFjxgycOHECLVq08Ft7kpKScM455+CZZ55BTk4Ojh49iv/7v/9zedyamhpMnz4d119/PTp37oz9+/dj7dq1uO6666z7rFmzBgkJCW51zSMiinSs6BARhdijjz6K+++/H//4xz+Qm5uLG2+80TreQ3bTTTchNjYWN998s924kXnz5uH666/HXXfdhTPPPBMTJ07EqVOn/Na+F198EUOGDMHll1+OCy+8EMOGDbNOCe1M165dce211+LSSy/F6NGj0bt3b5vpoydMmIA333wTCxcuRJ8+fTBixAgsXLjQ44rO1VdfjZdffhnPP/88evXqhfnz52PBggUYOXIkANFVbNGiRairq8PgwYMxYcIE6/TdzsiVtE8++cSv7QGAt99+Gw0NDRg4cCDuu+8+PPHEEy6PGxMTg2PHjmHcuHHo3r07brjhBowZMwaPP/64dZ8PP/wQt9xyC5KTkz1qMxFRJDJJ3nS0JiKioNq3bx9ycnKwdu1anHXWWSFty6lTp9C+fXvMnj3bpuqk9thjj+GLL77wqEtWOFq8eDEeeOABbN682Tq2JlwdOXIEZ555JgoKCjwOi0REkYhd14iIwlhDQwPKysrw8MMP45xzzglJyNmwYQO2bt2KwYMHo6KiArNmzQIAXHXVVUFvS7Bdeuml2LFjBw4cOIDs7OxQN8ep3bt3Y+7cuQw5RER/YNAhIgpjK1euxKhRo9C9e3f897//DVk7XnjhBWzbtg3x8fHIy8vD8uXL0apVq5C1J5juu+++UDfBLYMHD7aZBpyIKNqx6xoREREREUWc8O5wTERERERE5AUGHSIiIiIiijgMOkREREREFHEYdIiIiIiIKOIw6BARERERUcRh0CEiIiIioojDoENERERERBGHQYeIiIiIiCIOgw4REREREUWc/wdtFxUc/hc/iQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(10, 4))\n", "gs = gridspec.GridSpec(1, 1)\n", "ax0 = plt.subplot(gs[0,0])\n", "ax0.semilogx(period[1:n//2],\n", " abs(y_fft[1:n//2])/np.max(abs(y_fft[1:n//2])), \n", " color='blue',\n", " linewidth=2, label='onshore wind') \n", "ax0.legend(loc='best', shadow=True, fancybox=True)\n", "ax0.set_xticks([1, 10, 100, 1000, 10000])\n", "ax0.set_xticklabels(['1', '10', '100', '1000', '10000'])\n", "ax0.set_xlabel('cycling period (hours)')\n", "\n", "#We add lines indicating day, week, month \n", "plt.axvline(x=24, color='lightgrey', linestyle='--')\n", "plt.axvline(x=24*7, color='lightgrey', linestyle='--')\n", "plt.axvline(x=24*30, color='lightgrey', linestyle='--')\n", "plt.axvline(x=8760, color='lightgrey', linestyle='--') \n", "plt.text(26, 0.95, 'day', horizontalalignment='left', color='dimgrey', fontsize=14)\n", "plt.text(24*7+20, 0.95, 'week', horizontalalignment='left', color='dimgrey', fontsize=14)\n", "plt.text(24*30+20, 0.95, 'month', horizontalalignment='left', color='dimgrey', fontsize=14)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.12" } }, "nbformat": 4, "nbformat_minor": 4 }