{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# SISO for Event History" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "from mdof.utilities.printing import *\n", "import mdof\n", "from mdof import modal, transform\n", "from mdof.utilities.config import Config" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "## Upload a Dataset and Identify Spectral Quantities" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "inputs = np.loadtxt(\"./uploads/opensees_sdof/Accgrd_set1.csv\", delimiter=\",\")\n", "outputs = np.loadtxt(\"./uploads/opensees_sdof/Accrsp_set1.csv\", delimiter=\",\")\n", "# inputs = np.loadtxt(\"./uploads/opensees_sdof/Accgrd_set2.csv\", delimiter=\",\")\n", "# outputs = np.loadtxt(\"./uploads/opensees_sdof/Accrsp_set2.csv\", delimiter=\",\")\n", "t = inputs[0,:]\n", "dt = t[1] - t[0]" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "### Method Inputs\n", "\n", "See [`01_SISO_Intro`](./01_SISO_Intro.ipynb) for parameter definitions" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set parameters\n", "conf = Config()\n", "conf.m = 300\n", "conf.horizon = 140\n", "conf.nc = 140\n", "conf.order = 2\n", "conf.a = 0\n", "conf.b = 0\n", "conf.l = 10\n", "conf.g = 3\n", "conf.damping = 0.0001\n", "conf.period_band = (0.1, 3.0)" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "### Perform System Identification" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "100%|█████████▉| 885/886 [00:00<00:00, 281528.94it/s]\n", "100%|█████████▉| 885/886 [00:00<00:00, 174254.02it/s]\n", "100%|█████████▉| 885/886 [00:00<00:00, 115929.89it/s]\n", "100%|█████████▉| 885/886 [00:00<00:00, 165130.08it/s]\n", "100%|█████████▉| 885/886 [00:00<00:00, 237474.19it/s]\n", "100%|█████████▉| 885/886 [00:00<00:00, 217175.23it/s]\n", "100%|█████████▉| 885/886 [00:00<00:00, 161995.24it/s]\n", "100%|█████████▉| 885/886 [00:00<00:00, 417919.28it/s]\n", "100%|█████████▉| 885/886 [00:00<00:00, 213958.10it/s]\n", "100%|█████████▉| 885/886 [00:00<00:00, 191288.79it/s]\n", "100%|█████████▉| 885/886 [00:00<00:00, 159874.19it/s]\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ " period(s) damping\n", " [[0.41068031 0.02342356]\n", " [0.41083068 0.02249289]\n", " [0.4107063 0.00068529]\n", " [0.41201371 0.02015356]\n", " [0.41051848 0.02368659]\n", " [0.4112227 0.02544649]\n", " [0.40933146 0.01991049]\n", " [0.41281326 0.0191038 ]\n", " [0.40928958 0.02311352]\n", " [0.41883737 0.0023086 ]\n", " [0.40943662 0.02034821]]\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAqkAAAFyCAYAAADWN0O2AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8fJSN1AAAACXBIWXMAAA9hAAAPYQGoP6dpAACjmklEQVR4nOzdeVzUdf7A8dcc3NcA4oVHDp6pqaB2n0LnRlmQS7W/aktpt3N3S7L2aK8Ut7a22i2wttraJZVqs+0SNLu9QK28ZfAakWsYbpjr+/tjnC8MDAiKcvh++uDBdz7f73zn850ZnPd8jvdHoyiKghBCCCGEEH2ItrcrIIQQQgghRFsSpAohhBBCiD5HglQhhBBCCNHnSJAqhBBCCCH6HH1vV+B0mDx5MnFxcR3uN5vNxMbG9tjjyfkG7vn6ct3kfHI+OV//PV9frpucT853us5XVFTE9u3bWwqUM8D1119/Uvt7+vHkfP33fH25bnI+OZ+cr/+ery/XTc4n5ztd52t7P+nuPwXS0tLkfH1IT9avrz93Z9JrcSaer6f19evt6+frSX39Wvv6+XpaX7/evn6+nqJRlIGfJzUhIYHY2FjS0tJ8vhDJycmsWrWqF2omfJHXo++Q16Jvkdejb5HXo++Q16Jv6e7rkZOTQ05ODmazmYKCArX8jBiTGhsb2+mT1Ve/QZyp5PXoO+S16Fvk9ehb5PXoO+S16Fu6+3p4GhGTk5O9ys+IllT5hiWEEEII0be1jddkTKoQQgghhOhzJEgVQgghhBB9jgSpQgghhBCiz5EgVQghhBBC9DkSpAohhBBCiD5HglQhhBBCCNHnSJAqhBBCCCH6HAlSRY/QaNw/QgghhBA9QYJUcdJ++KG3ayCEEEKIgUaCVHHS1q3r7RoIIYQQYqA5I4JUs9lMcnIyOTk5vV2VAamgoLdrIIQQQoj+Kicnh+TkZMxms1e5vpfqc1rFxsZ6rQUrelZxcW/XQAghhBD9VVpaGmlpaSQnJ3uVnxEtqeLUKi3t7RoIIYQQYqCRIFWctOrqlm2LpffqIYQQQoiBQ4JUcdIaGlq2t27ttWoIIYQQYgCRIFWctObmlu3du3uvHkII0VdYrVZSU1NJSEhAo9Gg0WhITU31+klISCAyMhKNRkN6enqv1zcuLo6lS5f2aj2EaK1XJ05ZrVYWL16M1WrFZDJhsVhYtGgRKSkpPo9funQplZWVREdHU1RURFJSUofHitPH7rJDcDU0DGLfvt6ujRBC9D6DwcDKlSuxWq1ERkZiNBpZuXKlz2PT09PZvHnzaa6hN4vFgslkYtOmTb1aDyFa67Ug1Wq1kpGRQWZmJgaDAYDCwkISEhJISUlp98ecnp5OXFwcmZmZallSUhIWi4UFCxaczqqLNpzX3QPT/gVPKjLTXwghWvF8vnl++5KVlUVCQsLpqVAHjEYjiqL0ah2EaKvXuvsXL17sFaACxMfHk5mZSW5uLvn5+Wp5YWEh2dnZLFy40OscmZmZvd5FIoC4TzmrCt7kdhSTRKlCCNFdiYmJvV0FIfqcXgtSc3NzfX5z9Pyhtm5JzcrKIj4+vt2xnrLc3NxTVEvRJS4dkY1wO/9GsVT1dm2EEKJfMJlM6nZcXBxWq7X3KiNEH9RrQarRaMTiI1+Rp2W19b78/HyMRqPP8xgMBvLy8k5JHUUXKTp1s6mpF+shhBD9SFZWlrq9YMEC9fNv6dKl6k9GRgbZ2dle9/MMjYuMjPRq7OmoPD8/Xy3PyMigsLCQpUuXkpSURG5uLiaTiaSkpHb387BaraSnp6v1SU1NpbCwsMvnF+JE9dqY1I4CS88bf9asWWqZyWTqsCskKiqq1wecn/EUTct2c2Pv1UMI0ac0NMCuXb1di+6bOBGCg0/tY7Qe0tZaQkJCuwnEnsDQ08MYHx9PQUFBu4DSU+6Zr+GRmJhIQUGB2lqbn5/PwoULycrKYvny5aSkpJCXl9fufuD+/E1ISGDNmjVq72Xbsq6cX4gT0eeWRc3KysJgMKiToY7X/WEwGKSLpLe5WlpSsdl6rx5CiD5l1y7o5flAJ6SgAHyMMDspJpOJ9PR0LBaLVyDXWnp6OlFRUe2CuszMTCIjI8nNzfXaZzQavYYMeBgMBp89lUajkRUrVlB8bIZrQUGB17wQXz2c6enpJCYmeg25MxqNJCYmkpGR4dXgdLzzC9FdfSpIzc/PJz8/n5UrV8obuz9RdJSEwZO6DPYro3q7NkKIPmLiRHfA199MnNjz5zQajWr3vie7TVvZ2dleQwBa8wSFrYPUqKgon0FqR+WeenQl44Cnnvn5+V5ZdTySkpJ8XkN3zi/E8fSpIDU1NZWsrCyvP8Ku/BEdj9lsJjk5Wb2dlpZGWlraiVZTtOXSczQMfh/wCAHOQb1dGyFEHxEc3PMtkgOBwWAgMzOTxYsXq2WeoW5RUVE+79NZ4NkdHZ3fF89QuqKionbjYgGfwWt3zi9ETk4OOTk56m2z2ey1v88EqampqSxatKjbOU8tFkuHk6o8YmNjWbVq1clUT3TGpSesCc535bHFcSUQ3ds1EkKIPs1gMJCUlKTe9tU9f6oet6s8AWd3Fs6R1lPRHW0bDVs3KEIfWRY1IyODWbNmtRuf49HR+Bpwt6TOnDnzVFZPHI9Lx1gLfNp0K+Nc/XCWhBBC9ILWE4I9n2MdtZZ2pUGm9bE9ofVEKSF6Q68HqdnZ2URHR7cLUFt3Ldxyyy2d/pG0/jYqekGrFFSyXokQQnSfwWAgJSWF5cuX+9yfm5vrcwyoryFvhYWFPTah2DNL3xdZTEecar0apObn52O1Wn22oLb+A/PkZGv7R+dJ4SErdfQeRVG8ZvcrrdNRCSGEALo2f2LZsmVYrdZ24z89M+zbDodLTU1t14CTm5vbYR7yrtSj7X7PypBLly5t9zipqandPr8Q3dFrY1I96Tg8Mxah5c3t2eeRmJhISkqKupSqR2ZmpmQC6GWNjkavllQXEqQKIYRnBr9n8pEnYb7RaCQ1NdVn44rBYKCoqIiMjAwyMjKIjo6mqKiIhIQEn62ZiYmJZGZmkpqaquYWj4+Px2g0snnzZjXnqtFoZPHixWrDjud4TwORyWQiIyNDbTjyzBHxdPcXFBR41cnz2J79hYWFnZ5fiBOlURSlV3po4+LiOu3CLygoaLcU6tKlS6msrFT/cLs6mDs5OVkmTp0i5hozIx6ex9lBX/N+1ghudn7ANmV6b1dLCCGEEP1M23it11pSi4qKun0f+VbW91ibrKDo2DEYxg1fDoem93aVhBBCCDEA9PrEKdG/FVcVo06X0tp7tS5CCCGEGDgkSBUnZfn25aCBqUeh7NANnMfXvV0lIYQQQgwAEqSKk6LRaMCpR++CGFc1IdT1dpWEEEIIMQBIkCpOiqIobfKkyltKCCGEECdPIgpxUuwuu1eeVJe8pYQQQgjRA86IiMJsNpOcnExOTk5vV2XAsTlt4GpJEiF5UoUQQgjRHTk5OSQnJ2M2m73Key0F1ekUGxsreVJPEXdLqp490XD+uEXs2Dujt6skhBBCiH4kLS2NtLQ0kpOTvcrPiCBVnDp2hx3QUx8A68NHA5G9XSUhhBBCDABnRHe/OHUcigNQiK2Gvx5+mzHs7e0qCSGEEGIAkCBVnJQmRxNoXAyuh1+UrmMYpb1dJSGEEEIMABKkipNSb6sHjaLeVmTilBBCCCF6gASp4qR4WlKFEEIIIXqSBKnipNicNglShRBCCNHjZHa/OCnNzmbQKFQEw0sjpnH08JDerpIQQvQJqampmEwmCgsLAUhJSfHab7VaSU9Pb1cuek5GRgZWq1W9HRcXx8KFC1m6dCkpKSkYjcbeq5w4LglSxUlxt6Q6OWSAn0+9GA6P7e0qCSFEn7By5UoANBoN8fHx6m0Pq9XKnDlzWLx4MWvWrMFgMPRCLQeuhIQEFi1a5PUlwGQykZqaSm5ubp/9cpCRkUFmZmZvV6NPkO5+cVIcLgdoFIJsEF97mCDqertKQgjRLxgMBtasWUNhYSFz5szp7eoMKEuXLsVoNLYLRI1GY58PAE0mU29Xoc+QIFWcFHeQ6mJiBRR89V+msL23qySEEP2GwWDAYDBQWFjo1S0tTk5eXl6HXfm+gte+Ijc3V94HrUiQKk6KU3GCxtmqRFJQCSFEd3iCEunu71n5+fkd7ktKSjqNNekak8nE/Pnze7safcoZMSbVbDaTnJysrg0reo7T6fKa3a90cqwQQghvubm5AB12QVutVjIyMoiLi6OyshKTycSiRYuIj49XjzGZTGRlZREXF6cGvEVFRaSnpxMfH09+fj4ZGRmYTCYSExNJT08nLy9PvW9nXeBLly5VtysrK4mLi2PBggVqWWFhIfPnz1fPvWzZMrKzszEYDGprZttzH6++3bn2jqSmppKenk56ejqZmZntvgC0voaTeX66U8fs7GwKCgrUuiQlJZGYmAi43wfLly8HYPPmzaSmpgItwxNa13HBggXMmzeP/Px88vLySE9Px2g0qq+D0WikoKCg3evTUbmv1y0qKoqsrCysVivZ2dkAbNq06ZQNl8jJySEnJwez2ey9QzkDXH/99b1dhQHL78lAhZ9eoMxYgKKAksCG3q6SEEL0KYASHx/vVVZUVKRkZmYq8fHxysqVK33er6ioSDEYDEpBQUGnZW3PrSiKkpiY6HWM5ziDwaBkZWV5laekpPg8h6+6LVy4UElJSfH5eImJiUpmZqZXucFgaFfWlfp29do7k5iYqOBuO1GMRqOyYMGCDp9rT7268/x0p44pKSnKggULvMoWLlyo5OXltatzYmJih3X0XIfnOTUajV6vR3x8fIfPb0flHb1uCxcu7NLr2ZPaxmtnREuqOHXsdgW0TlwaqPHTYrBXQ2kpDJFUVEKc6RrsDeyq2NXb1ei2iYMmEuwX3KPntFqtXq2SmzZtwmQysWzZsg5bBtPT00lMTPTabzQaSUxMJCMjg7y8vA7Hsqanp7cri4qKArxbEQGWLVtGZGSk16zy9PR0oqKi2o3dzMzMJDIyst3seIPBQG5uLllZWV7Hz5w5k+XLl7Nw4UKALte3K9d+PHl5eWRnZ5OVlUVhYSHZ2dlqa+GaNWvaPe/deX66U8fs7Gxyc3OpqqryOm9ubq7aktlVRqORFStWUFxcDODVMuvZ72vilcFgwGKx+Czv6HVbunQpRUVF7cpbv56nmgSp4uTYQkHjYtswiFgwHv6eCENkXKoQAnZV7CIhO6G3q9FtBQsKiB92/C7l7jAYDO0+2HNzc0lISGDlypU+c6jm5+f77FpNSkoiIyMDaAlKPGWegKejiUGeQKxt3VJSUsjOzlYfzxPc+eIJwlo/RlRUFAaDod1kJaPR6BUcdaW+Xb32rliwYIEadHq6x5cuXUpCQgJVVVXthgF09fnpTh09z1Xbx4qPj2fevHldvhYPo9GonstX/X0FqZ2Vd/S6dVR+OrMPSJAqTk5zeMuYVJ2td+sihOhTJg6aSMGCgt6uRrdNHDTxtDxOSkoK8fHxzJ8/v11QuXnzZsA9VtMzJrA1T3BkMBhYuXIl8+fPVycDxcfHk5mZ2e0WOqvVitVqVYMQXwGbp7yjgOd4ulLfrl57dyUmJpKYmEhSUhJJSUksXry4y+dq/fwYDIYu19FzH1+ZBtrmze2qrjzPPXG+nn6cEyFBqjg5tmDQOJlUBivfPsiDrIEdw+Hss3u7ZkKIXhbsF9zjLZIDjdFopLCwUJ3Y4uEJEJKSko6bLiklJYWUlBS1pTA3N5ekpCTy8vK6Fah6+OoW7oquZic4Xn27c+0d6SwhvidY9awEdiK6WkfPcxkdHX1Cj+MJilvr6SwQHZ2vL2SbkBRU4uTYg0HjItABky0OAmmGpqberpUQQvQrbQMmzzjH43Wt5ufnq6mWEhMTyczMpKioiJSUlG611BUWFqo5W2fOnNnpY1sslhNeTrQr9e3qtXfmeAFofHx8t4Kw1s+P5/5w/Dp6nqe2Yzu7avPmzT3SvX6iXzx6mwSp4uQ4grxSUNmlcV4IIbrM0yLXOhDxBFgLFy7scFxo64lGvoLR9PR0n4GJrzLP+MpFixYBLWMwPSmR2srNze3WuNC2ulLfrl57RywWi89ueI/8/Hyf40G78vx0t46eVmNf2taxbeBstVq73e3ua2Jaf10sQoJUcXJswaBtSeZvxw+7vRfrI4QQfYQnuAF3EGoymdoFCp5gpvVsdU9w6Mnv2TorALiDRE8eTXAHOm1b2/Ly8nwGYSaTqV3ANH/+fBITE70mdi1btswrR2br+iYmJrabAW+xWHwGQb7Ku1Lfrl57Z/Ly8sjIyGj3+EuXLvWZuQC6/vx0p47Lli0DaBfYFxYWtgtAk5KS1PGunvr4Clw7kpqa2u65zc3NbTeBzaM7r5un/HTSKIoy4POvJycns2rVqt6uxoCkSZ0LQ3cyw7abwmy4mM9579NQBl0p49CEEGc2X62N0dHR7YIdzyxxTxeyJzl72/N4xjW2TnuUn5+vdkW3Dip8Lf3pmajkCdxMJhNFRUUkJCS0Czp9PbavY00mExkZGeTn52O1WklJSVHr37rckyTfswRsV+p7vGvvTHp6OllZWerCAZ7HslgszJo1y2cKpRN5frpTR8954+Li1JnzvsYMZ2RkUFhYSFJSEvHx8er42cWLF6uLP6SkpHR4HUuXLmXTpk3MmjULcA9NWLlyJStWrMBoNKqLDXTldTve69nTy8u2jdckSBUnRXNzGsQWEBG8l0uKtXy1oow1a3XMuNzQ21UTQgjRiicI60qO0TORPD+9r228Jt394uQ4AkDjpNpfxweTXFQRibne0Nu1EkIIIUQ/J0GqOGGKooAjEDQuhtRqeOwrGKI9TOm2o71dNSGEEEL0c2dEkGo2m0lOTiYnJ6e3qzKgOFwOcAaAxsXwWg2L18BwXTGWHSW9XTUhhBBtWCyWfpuK6HSQ56f35OTkkJycjNls9io/I4LU2NhYVq1aRVpaWm9XZUCxu+zg9HfP7leOLYWqs1FZKcuiCiFEX1FYWEhqaiqFhYXqdkcpkc5E8vz0vrS0NFatWkVsbKxXuSS1FCfM7rSD088rTyr6ZixWCVKFEKKv8MzuFr7J89N3nREtqeLUsDltane/St9MTa28rYQQQghxciSaECfM7rKDyw80TqyBsPJssOr9KW2K6O2qCSGEEKKfkyBVnLAmR5N7TKrGRXGklltugeKgaIrso3q7akIIIYTo5yRIFSfMHaS6x6T6OZ3EVoOfrhZnk6yLKoQQQoiTI0GqOGFNjiZwuWf3TylXOPwsTHHuIrZxX29XTQghhBD9nASp4oQ12Zsg/hUIqG2ZPKWzYXfK7H4hhBBCnBwJUsUJq7fXw/BC9w1PkKpvxiFBqhBCCCFOkgSp4oTV2epabnjiUo0Tp0uCVCGEEEKcHAlSxQlrsDe0L9RKkCqEEEKIkycrTokT1roldetQCHgC7B8PQ0tcL9ZKCCGEEAPBGRGkms1mkpOTSUtLIy0trberM2DU2mrVbUULNg2gd6JIA70QQgCQkZGB1WpVb8fFxbFw4UKWLl1KSkoKRqOx9yrXAwb69YnTIycnh5ycHMxms1f5GRGkxsbGsmrVqt6uxoBT01Sjbo+rgOwPYIHhCLWuEmBY71VMCCH6gISEBBYtWkRKSopaZjKZSE1NJTc316vcIyMjg8zMzB55/J48ly8ncn1C+OJpRExOTvYqlyYvccJqmluC1FAbXHYAQpUGAmjuxVoJIUTvW7p0KUajsV2gZjQaOw0cTSZTj9WhJ8/V1olenxDdIUGqOGHVzdXqtkY5tqF1oiATp4QQZ7a8vLwOu7p9BXcAubm5Xl3nJ6Mnz+XLiVyfEN0lQao4Ya0nTnneSFqNTYJUIYQA8vPzO9yXlJTkddtkMjF//vweedyePFdnunN9QpyIM2JMqjg1Wqeg8uTy12oduCRIFUKc4VJTU0lPTyc9PZ3MzEwMBoPX/gULFqjbubm5LF++HIDNmzeTmpoKtO86N5lMZGVlER0dTWVlJUC7rvWunstqtZKRkUFcXByVlZWYTCYWLVpEfHx8j19ffn4+GRkZmEwmEhMTSU9PJy8vT72mjoYIdOV6W8vOzqagoECtS1JSEomJiV7HnOx1i9NLoyiKcvzD+rfk5GSZOHUKXP765aw7sA6AQXVw8054t/Zn1H/xF+qVkN6tnBCibygpcf+0FhkJY8ZAUxPs2NH+Pp6AYfduqK/33nfWWRAVBeXlcOiQ976wMBg3DpxO2Lat/XmnTgU/Pygqgupq732xsTBkSLcu7XiSkpLU1kaj0UhiYiJJSUkddoV7Wh89AVxb6enpZGVlqbczMjLIz8+noKCgW+cymUwkJCSwZs0aNTjzVdbT15eQkIDJZCIzM9MriE1NTcVkMrW7ju5cb2pqKlFRUe2Obx2o9tR1i1OnXbym9LKqqiolJSVFWblyZYfHJCYmKgsXLlQKCgoURVGUoqIiZeHChcqCBQu69BjXX399j9RVeJudPVvhSbx/rnpQ6f13lRCiz/jd7xQFvH9uu829b+/e9vta/wdy3nnt9735pnvfiy+233flle591dW+z1tW5t5//fXt9z3zzCm5/KysLCU+Pl4B1B+DwaB+nrWWmJioJCYm+jxPXl6eAih5eXlqWVVVlQL4/Pzs7FyJiYlKSkpKu/KUlJQO79OR7l5ffHx8u3LPdSxcuFAt6871ZmVlKYBSVVXlVW40Gr2usyevW5wabeO1Xuvu93zrAXf3xLx58zo81mKxsHTpUpYuXaqWJSYmdvhtU5wejY5GdXtwLdy0E95zWKmhHpCWVCEEkJ4ObdLKEBnp/j1iBPhoFVO9/rrvllSAW26B88/33hcW5v4dEuL7vJ4u6WefhSef9N4XG9txPU7CggUL1FbD/Px88vLyWLp0KQkJCVRVVbXrJu+Ip6Wy9WQlz327M4vfarWSn5/vs9s8KSmJjIyMLp8Lun99ns/91gwGAykpKWRnZ6v16s71ZmRkkJKS0u6x4uPj1diip69bnB69FqSuXLkScL/ZsrOzOz3WaDQyb948Nm3ahNFo9DnORJx+TfYmddtYBS99BD+cX8VeapEgVQgBwLBh7h9fAgNbuvZ9mTCh430xMe4fX3S6zs8b1zur4iUmJqpd4klJSSxevLjL6ZqMRqPaMGO1WjGZTGqw5hmv2RWbN28GoKioyOdn78mkjzrZ67NarVitVgwGQ5ev13MfX5kGPHEGnNrrFqdOv5g4FRUVxcKFC3u7GqKNZmdLPlTPVCmt1i6z+4UQZ7zOEul7grnCwsJOz+EJ2Dxyc3PJyspSWwh9tR4e71yelszOxo52RU9c3/F05XotFgsA0dHRnZ6rp65bnF79IkgVfZPNaVO3PXlSNRKkCiHEcQO0+Pj443bTb968GaPRiNFoJDs7m4yMDAoKCk5oqVHPuVpPGDoZPXF9bc9nMBjUILSr1+vZV1RUdNz6wKld4ED0vH6VJ7WwsJDs7OyT/nYmeobdZVe31SBV45AgVQhxxrNYLJ0OZcvPz283F6NtK6HValVbADMyMrjlllvaBWytE/a3nrfR2bkWLlzoNQu+tfT09A7r3NqJXJ/nfm15xosuWrRILevO9aakpHSYs7V1HXviusXp1S+CVIvFQkZGBhaLhQULFmCxWNRUFqL3OJwOdbvWH74dAdUBCjb8e7FWQgjRN+Tl5ZGRkdFu5aelS5cSFRXVrts5KSlJHTsJ7lY/T7AZFRXV7jMvNzeXxMREnytLdXYuT17T1kGt53yevKqn4vo89WgbUM6fP5/ExESvYX3dud5ly5YBtJv8VFhY6DVRq6euW5w+vZ4n1WQyERcXx8qVKzscJ5Kdne2VUw3cfwRZWVnHbeIHd2622FYzN9PS0khLSzu5igvCngqjzl7nXbjvSnjrUwZ+9l0hhOiYJ8enJyG9J7CyWCzMmjWrw3kWGRkZFBYWkpSURHx8vDpJ2Gq1Mn/+fKxWq5oD1fOZmZ6eTnx8PIsWLfJqQe3oXK33Q8t4zsTExC7nCj2R6/PU2xPYmkwmioqKSEhIaPcZf6LXa7VaiYuLUydf+ZpkfTLXLXpWTk4OOTk56m2z2eyVB7dfBKm+FBYWkpCQ0KX7STL/UyP4z8FqGiqNC/wdYDt0Kcqbn6Eo0uUvhBCixfEWKxCibbzWL7r7ffGMU5E3e+9xupzq9kUHoOkpuOJoFQasvVcpIYQQQgwIfT5ITU1NJSEhocP9vgZhi9PDpbjalWm0DlwycUoIIYQQJ6nPp6AqLCz0mX7CE5zOmjXrdFdJHKPQMlJEDUu1TpndL4QQoh1pVBLd1edbUlNSUnx26efm5gK0G2wtTp/Ww5k9KajQSJAqhBCiRWFhIampqRQWFqrbHaWMEqK1Xm9JbT0j0JdFixapswg9CgsLWbx4MStXruzyahui5ykooGjgywyChi5xF2odyMR+IYQQHvHx8V5LlArRVb0WpGZkZGAymdTE/BkZGeTl5REVFeUVkBoMBjIzM0lPT8dgMKh509asWSMpI/oCpz+UTmPdeTDh53BA0WHTBPV2rYQQQgjRz/VakNrRmr++GAyGDleJEL1HQQGnH+jsNPnDnsFAlQuku18IIYQQJ6lbY1K3bt3KVVddRXR0NK+88opafu+997J27doer5zoB1x+oGtm1iH47DU4p6IJP5p7u1ZCCCGE6Oe6HKRu2bKFK664goiICJYsWeK17+WXX6aqqoqtW7f2dP1EX+f0B62dYXVw2QGIabSjlVGpQgghhDhJXQ5SlyxZQkFBAStWrGD+/Pnt9t98880yW+9M5PQDnQ1PXKpoXTK7XwghhBAnrctB6pgxYxgzZsyprIvoR9T0Uy53kKq+kSRIFUIIIUQP6HKQOmjQIK/brXNkelRWVp58jU4Bs9lMcnIyOTk5vV2VAcPusrs3nP6gs6lhqUvrQOn76XeFEEII0Ufk5OSQnJyM2Wz2Ku9yNLFv3z62bdum3tZovFvLnn766ZOs4qkTGxvLqlWrSEtL6+2qDBiN9kb3htMP9Db2RMPr50BxJDjR9W7lhBBCCNFvpKWlsWrVKmJjY73Ku5yCasmSJRiNRpKSkpg1axZFRUVERUVhMpnIysrCYDCwadOmHq+46JuaHE3gQp3dv30I3HUTYFMALXY7+Pn1ciWFEEII0W91uSXVYDCwefNmLBYLCxcuJCsri5SUFBYuXMjNN98sAeoZptHRCC6d2t0/ygKPfQ7RDU5QtJSW9nYNhRBCCNGfdSuZv9FoJC8vj+rqajZv3kxUVBQzZsw4VXUTfVijvRFcenV2/3mHYfFn8OVIJ1+jYcsWGDGit2sphBBCiP7qhGa4REREMGfOnHYBquRJPXM0OZpA0auz+z0jlBWte0Ld9u29VzchhBBC9H89Og3bV/5UMTB5d/c3t5rd7wKgqKj36iaEEEKI/s9nd/+8efO6dRKr1YrFYsFkMvVIpUTf12hvBEV7bHZ/M8diU5wad0tqmywSQgghhBDd4jNIzcvLw2g0EhUV5VW+efNmn+WeAHXWrFmnrqaiT3F392vV7v46HVQGQb2/O0jtoylzhRBCCNFP+AxSjUYjmzdv9ipbs2YNAHPmzPF5omXLlkmQegZpdDQCGnd3v18tqybAoLOP7dS4sFolob8QQgghTpzPSCIzM7NdWXFxcYcBKrjHoy5fvrznaib6tLbJ/L1onNTXn/46CSGEEGLg8Bmk+gpGrVbrcU/WdhiAGLga7A3ujWPd/bdtheY/wLQSQOugqak3ayeEEEKI/q7LfbJFXZiu3XapVDFw1dpqQaOos/v9nODvAq0CaB3YbMc9hRBCCCFEh7ocpBqNRhYtWtTh/nfffZeKiooeqVRPM5vNJCcnk5OT09tVGTBqm2vdG04/0Nq930h+jdjtvVErIYQQQvQ3OTk5JCcnY26TGqjLK049+uijJCQkkJ2dzbx58zAajYC7hTU/Px+DwdBnl0aNjY1l1apVvV2NAaXOVuduSXW5U1AdyzzlFliFs2Zwr9VNCCGEEP1HWloaaWlpJCcne5V3a1nUgoICMjIyWLZsmdcY1YULF7JkyZIeqajoH+psdYCrJZl/6yA1yIJS3Vs1E0IIIcRA0K0gFdwz/zMzMykuLgZgzJgxPV4p0ffVNNeAxuXu7tc5WDUB4iywaxAQXInT2ds1FEIIIUR/dsLJLMeMGdMuQH333XdPukKif6i317d092vtlIXDoqug0R8IsvR29YQQQgjRz/VoxvXFixf35OlEH9bSkuoPOjuJe6DuTxBXCQRVoijHPYUQQgghRIe63N3f2WpSVqsVi0Vaz84kdc11oHGqs/snWCDEAeMroSi4vLerJ4QQQoh+rstBalFRETNnzsRgMHgl7bdYLBQWFpKQkKDO+BcDX729HrSuY8n87erEqfBmILhvpiITQgghRP/R5SDVaDSyevXqDve/8847REZG9kilRN/nHpPqdHf3a+3uJP5ARDNogiqQ3n4hhBBCnIwuj0ldtmxZp/tvvvlmCgsLT7pCon9otDe6W1KdvlpSZeiHEEIIIU5Ol4PUGTNmnMp6iH6m2dHs3nDpQOtiZwzU6WF3NCiB1l6tmxBCCCH6v27nSe1MX11xSvQ8u+vYuqeKBoDV4yHyt3ocLgdUNPVizYQQQggxEPTY7H6TyURmZmaPVEr0fQ6Xw73hjlEZUwl/znfwuytgr19j71VMCCGEEANCt2b3G41GjEaj1+x+j5dffpk5c+b0aOV6itlsJjk5WV0bVpw8l+Jybxwbi5qyHdJ2gjUE7ru0SSZOCSGEEKJLcnJyyMnJwWw2e5V3a3b/5s2be7xip0NsbCyrVq3q7WoMKC5c3gXHWlQDHKDom3FHr5rTXS0hhBBC9DOeRsTk5GSv8h6b3S/ObE3Hvu4EOAF9c6/WRQghhBD930nP7t+/fz81NTU9ViHR96ld/YCnv3/EsbdAgAPQ2aBtS6sQQgghRDd0OUhdtGhRu7Lq6mqKiorIy8vj6aefZu3atT1aOdE3NTlazd4/1qPv1LrDVYcGd/5UnbSmCiGEEOLEdWviVFsRERFek6Wefvpprrjiip6pmeizGu2tZ++7W1IfT3L/qPzrgeDTWS0hhBBCDCBdbknVaI4/CSYvL++kKiP6B68gVdPSrX9e7Hkt5X71p7FGQgghhBhoOmxJfeyxxzCZTFRXVwOwefNmrrrqqg5PtHnzZhYsWNDzNRR9TnljecuNY+uhPrkGfvvlej4YBzfcBvhJQn8hhBBCnLgOg9QlS5YAkJuby4IFC9BoNCiK7+yXBoOBJUuWMH/+/FNTS9GnlNaVttw41pJqaHIPTx1Re6xcLwn9hRBCCHHijjsmNSUlhfj4eB577DFWrFhxOuok+rjy+tYtqU73r2PfXwKPLUSFXlpShRBCCHHiujQm1Wg0Mm/evFNdF9FPlDe0DlK9U035O49tyOx+IYQQQpyELk+cuvnmm497zNatW0+mLqKfqGioaLlxLEj1vJEC1JZU22mtkxBCCCEGli4HqV2xePHinjyd6KMsjZaWG8eC1NenQcEw+MWVx8p1zTid7e8rhBBCCNEVPsekPv3002zevJm3335bLZs1a1anJ7JarZhMpp6tXQ8xm80kJyera8OKk1PVWNVyQ+uORAtHavj5H2ax0bzRXa6389130MFCZUIIIYQQAOTk5JCTk4PZbPYq9xmkvvTSS+zfv98rSC0qKmLmzJkYjUafD6AoChaLxee+3hYbG8uqVat6uxoDRnWTOy0ZTj1o3f37ad8pZP1pI/+eCvfeAOgbWbtWglQhhBBCdM7TiJicnOxV7jNILSwsbBdwzpw5k9WrV3f6IFVVVZ3uFwODtdnq3nAGgM499vSyYgh1wPwtniC1icLCXquiEEIIIfo5n0FqREQEERERXmVZWVnHPVlmZmbP1Er0adXNx1pSHYGgd8/i97yRtIDGpUXxa2T/QTvg1xtVFEIIIUQ/1+WJU2PGjDnuMcXFxSdVGdE/1DXXuTccAWqqKZ3S8lbSuLSgb+JoheRKFUIIIcSJ6dHZ/V1pbRX9X4O9wb3hCAKtHQCdS6Pu1yga0DdSXSdpqIQQQghxYnx29x9vJr8vfXl2v+hZzc5jifodAaBzB6nmwQE4NA1sHg4oGvBroKFRclAJIYQQ4sT4DFKPN5PflxOd3W+1Wpk/fz7z5s0jJSWlw+OWLl1KZWUl0dHRFBUVkZSU1Onx4tSxu+zgwj0m9djs/o9um0XGRZ8TExyDs6Ya/BqwOzo/jxBCCCFER3wGqUaj8bgz+X3pzuz+1NRUoqKiAMjNze102dX09HTi4uK8JmYlJSVhsVhYsGBBt+spTo7D6QCXn1eQmlgdzXmfw76xDt4drKD41+NyKL1cUyGEEEL0Vz6D1JUrV57Qybozu9/zGCaTiezs7A6PKywsJDs7G0XxDngyMzNJSEiQIPU0cykuXLjAFXAsSHV391+SvZpLdgCfVRGaoaPerwGXy+fbSwghhBDiuHxOnOrKTP6amhpqamq6fb/uysrKIj4+vl25pyw3N7fHH1N0rNHe6N5w6d15Uo8FqfpW3yF0OMGvHlw9Oi9PCCGEEGeQbkcRixYtIjo6msjISCIjI4mOjuaZZ545FXUDID8/v8OxsQaDgby8vFP22KK9Wlute0PRuFtSde7ufn2rt5IWwL9OglQhhBBCnLBu9cfOnDmTwsJCUlJS1MCxoKCARx99lLy8PD755JMer6DJZCIxMdHnvqioKDZv3tzjjyk6Vmc7liNV0R1L5u9OR6VVWlJQaRXAvxZ0MrtfCCGEECemy0HqY489htFoZM2aNe1Wo7JarSxYsIBXXnmFe+65p8cqZ7VaO91vMBiOe4zoWS1BqtYdpAZYAdC1C1LrQSPT+4UQQghxYrrcH2symVixYkW7ABXcweKKFSsoKCjo0cr1FLPZTHJysvqTk5PT21Xqt9QgFY1XntTtK15A8yQ88/XTVITgDlK1rlbHCyGEEEK0yMnJ8YrPzGaz1/4ut6R2JWdqd/KqdoXBYOh0f1dbUWNjY1m1atXJV0hQ23xsTCrKsTGp7lWlJsdMBiBQF+jerW8EjYPSulJCo0J7oaZCCCGE6MvS0tJIS0tTbycnJ3vt79GZLb6C1LVr1/bkQ3ixWCzHDWRFz1JbRjXuIHWkLZGGxxsYc+1tuJ6EO6/4BYPqAb8G0Do5Wne0N6srhBBCiH6qy0FqUlISr7zySof7165dS2RkZLvyrKysE6vZMQaDocOVrKxWKzNnzjyp84vuqW6qdm9oXOAIZETkMIL8gvCvqEIDhDTYCbQD+ibQOiipK+nN6gohhBCin+pyd39WVhZr1qwhKytLXSnKw2KxYDKZmDlzpldCf4vFQmFh4UlV8JZbbul0Bn9SUtJJnV90T0VjhXtD6wBnAJON7i8mOtpMnNI5QNcsLalCCCGEOCFdbknNz89nzJgxREZGoiiK109kZCQJCQntytuuEnUiUlNTKSwsbDf+ND8/H6DD9FTi1KhpPraAg84OjkCuvjQaAL2rTZAKcGeiBKlCCCGEOCHdmjh1IjlJb7nllk73e4LPjrr0ExMTSUlJYfHixV6ttJmZmaxcuVLGpJ5m6phUrR2ceiaNGgyArtUxapAaViJBqhBCCCFOSJeD1NYBYncsWrTIZ3lGRgYmk0kdDpCRkUFeXh5RUVHtxrGuXLmSpUuXkpGRQXR0NEVFRaSnp5OSknJCdRInrt5W797QukDREBMc4779k5+wNzuTb6ZFUhZSpR7fYG/ohVoKIYQQor/rcpA6Z86c4x7z7rvvctNNN3mVzZgxw+ex3Q16Fy5c2K3jxanhFXQ6A4kMOjZZbskSrh/xXw7XHKbe3nJIs6P59FZQCCGEEANCj6agOtmZ/KLvq7fXt9xoDkGvPfY952c/4/OFu8l8t57oesDpB0Cjo/H0V1IIIYQQ/V6XW1IBnn76aZYvX95hEn2TydQTdRJ9mNcKUrZWSfpXrWJII9xXAK8mQOW+W2Dav2lyNJ3+SgohhBCi3+tykPrYY4+RnZ3NzJkzSUhIaLe/srKyw8lPYuBoWXEKsIV1fKDGAcDRQ8GnuEZCCCGEGIi6HKSaTKbjBqHHm8kv+r9aW+sgtZPlTrVOAMpMg09xjYQQQggxEHV5TOqsWbOOe8yJZgA41cxmM8nJyeTk5PR2Vfo9dUyqPQgUXccHatxBam1FJ62tQgghhDjj5eTkkJycjNls9irv0YlTxcXFPXm6HhMbG8uqVatIS0vr7ar0e+oY06ZwoNViDUOH4gJ2R0JtAO48qoC9MfB0V1EIIYQQ/UhaWhqrVq0iNjbWq7zLQeqCBQt4+umn2b9/f4fHyOz+gU8NUpsjoNVSqBQU8MuPH2LiQ3Agxg/07uMUZ7fm5gkhhBBCAN0YkxoREUFFRQVxcXEYDAaioqK8VnuyWq0yu/8M4HC5J0TRZACl1XecgwcZnLeeGAUcg4Ko8js2LODkV8YVQgghxBmoy0Hqvffey4oVK5gzZw5Go7Hdfpndf2ZwOFsFqa5WY1LPP5/HjxzhceCqh/1ZHVzuLtfbTncVhRBCCDEAdDlItVgsMrtf4MLl3rCFgMvP5zGGQAMEHXLf0MuKU0IIIYTovi4HqUlJScc9pq/O7hc9x6W0ClLtQT6PCQsIA82x4FRr93mMEEIIIURnujxxqqNVplrrq7P7Rc+wOVt13dvCwO47UX+IX0jLDa3jFNdKCCGEEAORzO4XXVbVWNVywxbmbk31IUjfuoVVQVFk9pQQQgghuqfL3f0LFizAarWSkZEhs/vPUNYmq3tDAZrDQWn19vngA3j5ZZg3D0v1W/A94PAHjYLD5cBP53v8qhBCCCGEL10OUvPy8pg5cyY333wzUVFR7fbL7P6BT11tCg04grzSpBIfD9nZAAR//IG7zOkPQLOjWYJUIYQQQnRLl4NUo9HI6tWrOz1GZvcPbI32xmNbCtiD0GlbdeOfdx5s2AAGA0OfTXGXOf0ALaYqE+cMPed0V1cIIYQQ/ViXx6QuW7bsuMfI7P6BrdFxLEjVAI4g/P1aNaUeOpZyymplJOHubZcfKBq2HN1yWusphBBCiP6vy0HqjBkzjntMX53dbzabSU5OJicnp7er0q+1tKQC9iBCQny/fWJCYtwbLj2g4ZN9n5z6ygkhhBCiX8rJySE5ORmz2exV3uUgtSv66uz+2NhYVq1aRVpaWm9XpV9TW1IBHIEYDL7fPv5a91hUd5CqSJAqhBBCiA6lpaWxatUqYmNjvcp7LEh95513yM/P76nTiT7IqyXVEcTIkb6PUydJKVrQuKhuqm1ZBEAIIYQQogtOKkh99913mTdvHjqdjltuuaVLCf9F/+XVkmoPYtasVjuvvNL9e+ZMXMOHubddegisQakZQnVT9WmrpxBCCCH6v24HqWvXrlUD09TUVFauXMmYMWN49NFHvfKmioGn3lbfcsMRxNy5rXa+9hooCmzahHbosSC10QBDt0L1aO/VqoQQQgghjqNLKajWrl3LypUrWbFiBVarVV1BaOHChcybN0+dVCXJ/Ae2Oltdyw17IOee22rn734Hf/4zzJ1L0JP3ustsITBiMxydIkGqEEIIIbqlwyB169atZGVleQWmBoOB+fPnk5qaSlZWFkuWLPG6z4oVK055hUXv8QpSHUFoWifzf+UVcDohN5fIe3/sLtM63L8VjQSpQgghhOiWdkHqX/7yF5YsWeIzMJ0zZ456XG5u7mmtqOh9DfaGlhuOwA6Piw6Kdm+ElkD1SDAcpNnZfIprJ4QQQoiBpN2YVIPBQHx8PIqisHTpUiwWCy+//LJXgCrOTC3LogL2oA6PCwsIQ4sWDPuh4B4Y+p13K6wQQgghxHG0C1Lnz59PXl4eVVVVREREMG/ePBYtWsTWrVt7oXqiL2nb3d8RjUZDiH8I6F1giQO9jcqGytNQQyGEEEIMFB2OSY2IiGD+/PnMnz+f6upqsrOzeeqpp4iLi2PevHlYLJbTWU/RB6hBqqIBR4D3ToMBjhyBgAAICiIiIIJaW606LrWmueb0VlYIIYQQ/VqXUlBFRETw6KOPsmLFChYsWEBeXh4FBQUsWrSI/fv3q8fN8kqcKQYaNUh1BAAa753bt7tTUDU1waRJRAVFuct17glT1ibraaunEEIIIfq/LqWgas2TE/XRRx+luLiYl19+meLiYqqqqigsLDwVdRR9RG1zrXvDFnLcY9Ug1d89jrW0vvRUVUsIIYQQA9BJrTg1ZswYlixZwvLly0lJSempOvU4s9lMcnIyOTk5vV2Vfk2d3d8Y3X7nqFGg0bh/tmwhJiTGXe7vDmzL6stOUy2FEEII0Z/k5OSQnJyM2Wz2Kj+pILW1BQsW9NkMALGxsaxatYq0tLTerkq/pgapDdGA4r3T6WzZVhTGRI5xbwe4x6K+v+v9U19BIYQQQvQ7aWlprFq1itjYWK/yHgtSATIzM3vydKKPabQ3ujcaYo577PQh090bgVYALI0y0U4IIYQQXdejQapneVQxMNmcNnBpoGEQ7SZOtXF2zNnujSB3cNrYoDvFtRNCCCHEQNKjQaoY2NxLm2qhMRqttvMg1RBocG8EV4JTj2ILklypQgghhOgyCVJFlzkUB6BAQzT+/m12vvgijBsHzz0HkyYR6h/qLg+sAnsIaJwM+ssgNh/ZfJprLYQQQoj+qNspqMSZS0EBrQKN0YSFtdk5d67755gQx7GW1tCjYA+GIHcr6u6K3cwcPvM01VgIIYQQ/ZW0pIouUZRWs/kbBjFiRJsDEhPd6adiY+HAAQJ0AWjRQkgZ2IJB707qH6gPPH2VFkIIIUS/JUGq6JImR1PLjYZoZs9uc8DOne7fR45AZSUajcbd5a8FXC1jA7zOI4QQQgjRAQlSRZd4pZBqjOKnPz3+fUIDjo1LdbW8zQ7XHu7hmgkhhBBiIJIgVXTJ4ZpWwaWvllQf1MlTzpYu/iJLUQ/XTAghhBADkUycEp1SFAUFBXNNq6XKGqO6dF81SLWFqGUHqw/2ZPWEEEIIMUBJS6ro1IX/vJDzXjmPkroSd4EtGFx+7Q/0NK2OGwdDhgAQHhDe7rAjtUdOVVWFEEIIMYCcES2pZrOZ5ORk0tLSSEtL6+3q9CvfHv4WgAB9gLugyeD7wPfea1c0NHSoe6N+kFpW3lDek9UTQgghRD+Xk5NDTk4OZrPZq/yMaEmNjY1l1apVEqB2k3uFKbevDn7l3mgY5Pvgv/8dQkPhoYegthaA4aHDW+7z8hZoiMbSYMGluE5ltYUQQgjRj6SlpbFq1SpiY2O9ys+IIFWcGFOVqX1hfYzvg596Curr4fnnYe9eAIaFDXPvG/EtHJ0OjQZsLhtbj27F4XKcmkoLIYQQYkCQIFV0aG/l3vaF9YMBTZfuPzzsWEvq4O9BawdbKFqNlvNeOY+ozCjvjAFCCCGEEK1IkCo6tM+yjwBdQEuBQsfd/T4MCz3WkqoDAqqhKZLIwEjsLju1tlqWfLWEt394m6N1R3u03kIIIYTo/yRIFR06UH2AwSGDW5VooDEaTdcaUlu6+wGCy6FuCMF+wXz/s++5Z8Y9vLT5JdLeSWPksyOpt9X3aN2FEEII0b9JkCo6dLD6IJGBkd6FTZH4+chAReCxhP06HZ4D1JZUgIhDUD2K6uZqpgyegjHSqE6gcrgc/H3T30/BFQghhBCiv+o3QWpSUhIZGRkUFhYCYDKZyMjIID09vZdrNnAdqD7QsrQpgEaBJgNRvnL5FxWBooDDAVOnAhARGEHZI2Vo0YKhGKyjqWuuQ1EURkaMBCA6KJr0hHSWfr2UktqS03BVQgghhOgP+k2QarFYWLp0KQkJCWg0GuLi4igsLCQrK6u3qzZgHaw+SJA+yLuwMZJp07p+jpiQGAL1gRCzC5qicOGiwd7AyHB3kDoqYhRPXPwEDpeDuOfjeGfHOz14BUIIIYTor/pNkGo0GsnMzCQlJYWFCxeSl5dHXl5eb1drwGqwN1DRUIGm7Uz+JgMPPeTjDnFxoNGAXg/ff++1Kzo4GoYUgC0MgNe2vqa2pI6KGMXIiJEUPVjEhaMu5A9f/OFUXI4QQggh+pl+s+JUVFQUCxcu7O1qnDEOWA8AUFpf6r2jycA11/i4Q1OT+7fTCXa7165BQYM4NPIb+MLdKvv+7veZHz8fgNERowF3IHvDhBt44OMHSM5J5vlrnucsw1k9dj1CCCGE6F/6TUuqOL32VO4B4Psy71bRDpdF7cTg0MHg54RAKwCldaV8ceALxkeN57wR56nHTRk8BYAP9nzA7e/eLgn/hRBCiDNYvwtSCwsLyc7OVidQiVNjd+VuwvzD2u9o9DVrqnOTB092b4SUAVBWX8aVb13JHsse0qa2LFU7OWayuv3t4W95ceOL3X4sIYQQQgwM/SZItVgsZGRkYLFYWLBgARaLhYSEBEwmH0t3ipO2p3IP46LGeRc6/cAe3O1zzRo+y70R6p69X1pbru7bWb5T3Y4JiWFc1Dhevu5lLhx5IZuPbO5+xYUQQggxIPSbMalJSUksWLBAvZ2YmMi8efNISkqiqKio0/uazWaSk5PV22lpaaSlpXVyD7G7cjejIkZReLRVi3WTgQ6XRH38ccjIgLvvhnHewe05Q85xbxwLUqkcS2jsEepsdcxdPpdd9+9Sj93zgHuYwSdFn1DeUM6h6kPqJCshhBBCDBw5OTnk5OSot81ms9d+jaIoyumuVE8pLCwkISGBlStXkpKS0uFxycnJrFq16jTWrH9TFIWopVHcMOEG3tj2RsuOivHw4m66+445UnuE2L/GgnkGDP0B9iURfs5X1DTX4Kf1w/YbW7v7/PT9n/La1tcA2HP/HsZFj2t3jBBCCCEGjrbxWr/p7vfFaDQCSCqqHmaqMmFtstLsaPbeUTOi4yVR5851p6AaPx7afBNSx7YO3QLNYVA3lNrmWsC92lSTo4mxz4/16t6PCmoZ+7r4q8UnfU1CCCGE6F/6RZCamppKQkJCh/stFstprM3At+nIJgC+K/3Oe8ehC9XVT9vZuNH9e+9eKPVOWxXiHwJAgH8AOP1B0aPgbo5VUNhdsZuiqiKe+vIp9T6tl2NdtXuVuoSqEEIIIc4M/SJILSwsJMrHWpye4HTWrFmnu0oD2kbzRkZHjOZo3VHvHebZTJjQ/fNpNVpC/UOJDoo+VmCHVkMGnljzBADVzdVqWWSQO0jVaXRUNlaypWRL9x9YCCGEEP1WvwhSU1JSfHbp5+bmAnhNqBInL8+Ux+VjLqfOXtdSWHY27L+Uxx8/sXNGB0VzpO4IuHQQXAnVo9R9H+77EIDqppYg1dPdP8c4h0B9IF8e/JImR9OJPbgQQggh+p1+EaQuWrSI9PR0r7LCwkIWL17MypUrMRgMvVOxAchcY+aHsh+4Ou5q7M5WK0e9kQ+2UFJTT+y8I8JHuDdcOhj1NRye3e6YioYKddvT3R/qH0pMcAy/XvtrZmbPpB/P8xNCCCFEN/SLFFQGg4HMzEzS09MxGAxqbtQ1a9YQHx/fy7UbWApKCgBIGJagjhsFoDGSDtNPAUyaBEeOwPDhEB3dbvfIiJFwCIg86C4ougqm5HodY641c6T2CMPDhqstqf46fwyBBg7VHGJ7+Xb2WfbJTH8hhBDiDNAvglRwB6pZWVm9XY0Br6qxCgBLY5vJaC7/zu+Yn9/p7pHhrXKdOvWgb253jMPlID4rnr0P7FUnW/lp/dTxqQD5pnwJUoUQQogzQL/o7henXpOjiW8PfUt5QzlB+iC+OfRNmyM6aUUFeO89d/qpv/0NGhvb7Va7+wEsY2HkN1B6tlqk1WiZOXwm5Q3lPPDxA4yKGIVOo+PemfdiCDSox+UXdx4MCyGEEGJgkCBVAPDalte44J8X8NZ3bxEZFMkXB79o2enowtvk/vvd6acefhh27my3O0AXAMDzVz8PukYY/SUU3gO4A1SX4uKcwecwInwEb2x7A1OVCcdvHVww8gI1SNWgYW3xWpwu58lerhBCCCH6OAlSBQD7rfsBKG8oxxBo4IeyH1p2VkwENPj5nfj5E4a789yeP/J8woZYIOIQ1LhbV/207hM32BsI0AWg1+r56uBX6n09k6guGX0J1iarOm5WCCGEEAOXBKkCgJK6EgDqbHUYAg0crjncstPsnok/duyJn3/m8JnUP17PzOEzmRQzCapHwPgPwalFURT0Gj0ldSU02BuYMXSG13ADT0vqFWOuAGB3xe4Tr4gQQggh+oUzIkg1m80kJyeTk5PT21XpszyJ+xvtjQTrg2l0tBpXeuhCAP7+95N7jGC/YABuO+c2OHgBTPgAtv0Em8tGZFAkeyr3UG+v57wR57HRvFG9X6h/KOCefBUREIG51uzz/EIIIYTof3JyckhOTsbcZln1MyJIjY2NZdWqVaSlpfV2VfqskroSgv2CsbvsOJU2Yz7NMwG4/PJOTqDTtWxrOp9klTYlDYyfQ7AFz4Ss8oZySupKqLfVMzZqLAerD6o5UT2/Q/1DiQ2PxVwjQaoQQggxUKSlpbFq1SpiY2O9ys+IIFUcX0ltCecMOQeAelu9985q4/FPcPAgKIr7Z8aMTg+NCYnh3kvmQlMExL8OipZzh58LgN1l559b/kmzs5m5y+dy34f3senIJsA9wSo2LNa9cpUQQgghBjQJUgU2p43KxkqmDZkGtHT9txwQ1uOPuSRxCVjGuG9s+QmBfoHqvj0VewAorS9l3YF15O5wJ/0vqy+TllQhhBDiDCFBqlCDUk9L6sGagy07XVqOmyMVYPJkdzd/YKDPFFRtRQRGoBviDkbZeD/bSrap+26ZfAsAT176JNt/vp0DDx8gYVgCD33yEDVNNTImVQghhDgDSJB6httZvpO3f3gbaEn15OVYmqiw4zWmWq3u383NPpP5+zJv2o3wbBEcTcC6uWXA613T7wLg+pzrufd/9zIyYiT3zboPu8vOJ0WfUFJbIrlShRBCiAFOgtQz2JHaI5z9j7PJyM9Aq9GSXZDd/qAS9/jSX/6y5x9/jnEOXPWI+8aqN9XyiKAIwv3DsbvsZBVkYW2y8uz6ZwF3LlWn4qSsvkw93qW4eOrLp/jHxn/0fCWFEEII0SskSD2D/fWbv6rboyJG+Z6QVDERgCef7PnHHxQ8CM5+D1DAEQx297CCqsYqrzGqq4tW833Z97xy/Stojg09+Ms3fwHg032fMunvk3hi7RPc9/F93Pj2jXx98GsqGyp7vsJCCCGEOG0kSD2DWZos6nZcZByV9T4CO+voU/b4w0KHuTfOew7QwKfPA7D4q8WU1ZepK1H95/v/EBMcw/9N+z/+fdO/AXh2/bP8cd0fuebf17Cnco96zvd3v89Fr13EU18+dcrqLYQQQgx0Tpez14fW6Xv10UWvKq4qVreHhw3n8wOftz+oYXDXTnbPPfDnP8PcuTBmTJfuMnP4TH5zyW/4I7+CDQ/CnmTgAfJMeeox/jp/Pj/wObdPvR0/nR+Xj7kcLVq0Gi2//fy3AOg0Om6deitvfvcm46LGsdeylyO1R/iu9DuKq4opLCnkvtn3MTiki9cihBBC9AGKomB32WlyNNFob3T/djT6vN3lYzrY57ntKXO4HHxy2ydcNfaqXrt+CVLPYDsrWmbhRwVG4XA51NtatLiqh8PBi9B35V3y+9+7f7pBo9Hwh8v/wIPnPkjMvz6Dgxd77Xe4HIyJHIOpyoQx0sgfPv8Dv1v3OwYFDaLZ2UytrRYAp+Lk06JPAZg9fDZ7LXv5ZN8nvL39bfVcr255lRWpK7hg5AXdqqMQQgjhi9PlpNZWS21zLbW2Wmqaa9Tt2mb37UZHY/tgsIOgsKOg0qW4ulUvf50/QfogAvWBBPkd+93mdrBfMNFB0T73eW4H6gOZPHjyKXr2ukaC1DNUbXMtpfWl6m1Lo8Vrv4ICq7Khfgg/vr0LJ7zrLnj9dZg5E/73PxgypMt1GRQ8iP99bONHE3WgaMAWCAGNKCgYAgwA/HJ1y8ytWlstzc5m7p5xN69ueRVw51DVarT46/0BsDZbeXPum5w34jyC9EHMy53Hpa9fyn/n/Zfrxl8HwDs73iEuKo7pQ6d3ua5CCCH6r2ZHszuYbBVcegLKDgPOTgLQzgToAgjxD1EDPl/BoCHQ4L6tC+w0qGwdOHZ2TIA+AK1m4IzklCD1DLWjfIfX7datjnAsSC26GoA33+T4Vq92/968GczmbgWpANeMuwZmPQfNYRBYo5a3/ga594G93PnfO/n60NcMCh7E3Ilz1SBVgwaX4mJ7+XY0aPDT+XH7OS3R9Wd3fEbqylRuePsGbpx4IzqNjvzifFyKi5eve5mSuhI2H9mMv86fX57/S6YMngK4g/ftZdsJCwhj2pBpaI6z5KsQQoie5XA5qG6qprq5GmuTleqm6s5bMG3tA0rPtt1l7/SxQv1DCfMPIywgjDD/MMIDwgkLCGNk+Ei1PDwg3Ocxbff76fxO0zM0cJ0RQarZbCY5OZm0tDTS0tJ6uzp9wg9lP3jdPt4f7qmm1Wi55VcbWPHy/TBoN9RHw+xstpZuBdxB6NiosQwNHQqAn9aPvZa96v1D/UOps9WxybzJHWC34afz48VrX2T94fW8s/MdtVyDhh+/82M1sLU5bRSWFKLRaCirL+NIbUvGg9mxs1kyZwmXj7m83fmFEEL4ZnfasTZZqWqqoqqxSt32BJzWJmtLAHrsd+t99fb6Ds+t0+h8BohhAWEMDR3a5YAyLCCMUP/QAdUK2Z/k5OSQk5OD2ey9WM8ZEaTGxsayatWq3q5Gn7K9fDt6jR6H4mi3T4MGzdGZdG8UzMl7be5rrPhyEqz4ABoGwbS3CApRaHS4u/6dLicjw0cCUGer4/vS79X7esan6jQ6nIoTu7N90P3mtjeps9URHRRNZWOlevyUwVPYWrqV1LNT+WDPB2wr3caw0GHcNf0upg6ZypTBUzhUfYgnP3+SK/51BbdOvZUXrnmBqKCo0/CsCCFE71IUhQZ7Q7sgs23A2dH+joJMrUZLREAEhkADhkADEYHu7XFR49y3feyLCIggIjBCDSyD9EHSwzUAeBoRk5OTvcrPiCD1THak9ggf7vmQ+Qnzvco3mjf6DFDB/R+Hc3kOoGHkyNNQyWOC/YK5K3EOr1U+DivegU+fpfH6dHVV1u9Kv+MPl/+BUP9Q/vTln/jm8DcMDR2qLuvqr/PH5rQB7uEKjfZGgvyC1PMXlBQQPyyeLw9+qZY5FIfaWvvoBY9S3lDOtqPb+PzOzxkXPU49bsrgKVw99mre+u4tHvzkQaxNVj689cNT/IyInlRWX4ZOoyM6OLq3qyLEaeNwOWi0N1Jnq6POVke9vV7drm2u9Q4yG6t8BpzWJmuHvW1B+iAigyIxBBqIDIwkMiiS0YbRTA+c7lXm2TYEGtTbYf5hEmCKTkmQOsDd/u7tfLb/M+6acRcvb36ZT/Z9wvKU5aw/vL7D+7gUF1S500ht3drFBxo5Eo4cAYMBIiJOuL5LEpfw2tYhMOMV2HwfDNsCs14GID47nktGXcKMYe5VsHZX7GZ89Hg1SJ0SM4XCo4XuzAS4mLt8Lt+VfkdmYiY/mfYTCksKuWz0ZV5B6pyz5rBm/xrAvaDBiLARVDdVewWoHhqNhp9M+wlbj27lw70SoPYXiqLw4sYXeTTvUZyKkzfnvsmPp/y4t6slzkCVDZWU1pfS5GjC5rR5ze5uPaPb5rRx/fjriYuKO6HH+Xz/58z/YD5FVUXHnRmuQUNEYES7AHJE+Ih2Zb4CzgB9wAnVUYiukCB1gPN0g89eNpstR7cA8FnxZziVjhP0usd0ur/dRnW1R3t9x0FvdwwOGcy8yfNYzv2wIxU+egEG7YQx7hyuXxz8Qg0yFdxdUENDhnK0/ig3TryRwqOFXD32aj7a9xGfFn3KuKhx3PHfO4iLiqPYWqwGJ9eOvZb84nxSJqeoQaoh0EB0cMtQAIA9lXsw15i9xqEaAg1Ym6w9cr19gcPlQK8dmP8VlNWXcdf7d/HR3o9IT0inqqmK2969jV0Vu/jZzJ8R7BfMPss+pg+dLi06okMltSWsLlrNrNhZnB1z9nGPb3I08e2hb1ldtJpvDn/D7ordWBotXRr7r9fo0Wg0/OmLP7Hhng3dDlTrbHXc8d87GBo6lF+e/0t1RniYv3vMZah/KCH+Iep2eEB4nx2HWdlQyfrD6yksKVRTMWk1WiKDIokKiiI6KJrRhtFMHDSRQH3g8U8o+p2B+cl0hvDkNe0swPD85+MJUAFe2fIKeq3eKy9qa/7ls7F1tzKFhfDyyzBvHlxyCfid+KzGf9/0b7J+lEW8XxKm366DNz6Daf+EufcAeE2MOlxzmEcueIQpg6cQExwDwBzjHD7a9xGXjr6UYWHDqGioIOnNJLQaLWEBYQA8fN7DfHDrB2w2b1bPpdFoGBQ8iIqGCrVs8VeL+de2f/HOLe9w48QbgYEVpL79w9ukvZNG6SOlA2qxA0VReH3r6zya9yhajZb/pf2P68Zfh9PlZHzUeJ768il+/3lLXl9Pwmqb08aGwxuos9WRMDyh28/Jx3s/5tZ3b6XZ0UygPpAAfYD7ty5ADRY8Za3Lj7tfH4Cf1g+dVodWo0Wn0aHT6rx+azQatBqte0y5RoOGY7d9bI+KGMWwsGE9/bQPSOv2r2Pu8rnq3/x5I87jqrirmDRoElqNFmuTlX1V+9hdsZudFTsx15jbjcE0BBiIi4xjRPgIhoUNw1/rj06rA407HWBlYyXmGjP7rfvd91WgsrGSCS9OYPLgyYyOGM3I8JGMCB9BqH8oAfoA9Fo9Nc01OFwObj/ndoaGDsXpcnLX+3dR2VjJmv9bc8Itsb3B6XKyq2IX6w+v55tD3/DN4W/YVbELgOigaMICwtxD0VxOqpqqqGluyQKj1WgZGzWWKYOnMCVmCpMHT2bq4KmMjx7vfp5FvyVBaj/2wEcPUN1czX9u/k+Hx9Tb2g9Yf3/3+52e17blJkDDlVd2ozLXX+/u7l+2DAoKID6+G3f2ptPqiAiMYN+iDdykWcp/f/1L2PZTaI6CH9+kHqdBg4LCm9veZP096wn1D8Vf509pnTv/q2cFrVnDZ7HpyCYuGX0J7+18D3B/8CTFJXGo9pDXY0cHRWNtsqqti4eqD6FBw7zceWy4ZwPTh7rHWTU7m2lyNLX79m532jHXmhkVMarPtk54KIrCojWLAMgryuO2c27r5Rr1jO9Lv+f+j+/niwNfcNvU23jmymcYEupOiabT6vjjFX/k/tnu/TanjZ++/1N+u+63/GPzP1hbvJY6W516rrNjzubS0Zdy2VmXcenoS9Xz+LL16FZSV6Zy4agLuW7cdTQ5mmh2NKtduZ73TOvtqqYqn8e0LetpWo2WF695kZ/N+lmPn7u7FEU54VZsRVFodjbTaG/E5rTR7Gx2/3a4n7MgvyCC/YIJ0rt/d5YSqKy+jEPVhxgfPZ41xWvYaN5Ig72BZQXLuGDUBbxw9Qv8dt1vyTfls/HwRlw+ppZqNVr8tH4E6YOwu+xqQ4C12Yq12cquyl0dPr5eqycyMJKhoUMJ8gvC5rCx17KX/VX7qbfVs/nIZmqaa7A77dhc7maEQH0giqKw+MvFvDn3TZ7f+Dz5pnxWpq5sF6C6FBd2p7tOnro5XA61zNdP6+MURcFP54ef1k/97a/zb1fmpztWrvVDr9V3+No6XU7WH17Pqt2r+OrQV2w9upUGewNajZapg6dy+VmX8+uLf80FIy/gLMNZ7c5jd9qxNFoothbzQ9kPfF/6PdvLt5NVkKXmAA/2C2b60OnED41n5vCZnDviXMZFjcOpONW/s9Z/l8dbran1bbvTjlNxLxvq+e1SXO5txb2t1+rRa/TotXr8dO7nw/Pjr/MnxC/Eq2Xb1+0zPeuABKn92A/lP6gThTpS3lDe5fNFBERQ3VwNtbEAfPrpSVXvpGk0Gt5blMGo2ss59JcPYdeN8OL3sGAq+Le0qGo0Gi7854V8cvsnzI6dzevbXic6KJro4GhmD5/NW9+/xTmDz+HpxKe58LULAbA5bXxX+h0pK1LcaaiODQMYFDwIgDv+ewdGgxFzrZn0hHS+OPgFP/vwZ3z9068xBBoAsDZZ1ZRYHk9/8zSPr32cYL9gpg6eyvSh05k2ZBrThk5j6uCpaktuT9tv3c9Lm14iJiSGu2fcTWRQ5HHvs610G/ut+wG4e9XdKCikTUnrtZYHl+Jio3kj9bZ69cPR7rSr2zqNjqlDpnbY3VpkKeJ3637Hf77/D2OjxpL/k3zmGOcA7sD1w70fotPocCku9lv3s6tyFzvLd2Jz2dho3ojRYOSy0ZcxcdBEgvRBHKw5SJGliHd3vstLm18CIDYslrNjzmZyzGTOGXIOI8JHqIHPAx8/wOCQwdybcC9+Oj+vD6+uXr8nIHAqTvdvlxO7y64GXw6nw32c4sDlcnl9OLoUlzsLhuJEURS1rPWPorj3bynZwgMfP8DXh75Ws1QoSvvUbYDa+tp6u6MywGt/Z2UuxcXKHSsJ8Qvhhgk3qOWN9kbq7fU02Buot9dTb6tvue1juzur8eg0OjVw1Wl0hAeEE6APwNpk5XDNYa9zeSZiatBQcKSAs/9xts/0djqNjpiQGGLDYhkaOpQhIUMYEjqEISFDiAyK9Ep/FKgPREFRXx+Hy6FOUrI0WtSJS2X1ZZTWl+JwOTBZTdTZ69BpdO2GCzQ5mgBodjZzXc51an1uf/d29Dq9VxDqq+6ng16rVwNavVaPS3G5x+M6GnEpLjWFVFRQFGdFnEVEYAQB+gD2WvZSbC3m7e1vo9O4/0/yvMc97+W2t52Kk7ioOIaGDqW6uZq65jq2Ht3KhsMbOh3idjxajdYrcX6ALgB/nX+73ozWvz2tvh19IWh2NlNvc09iO94wEA0awgLCiAiIIDwgnIjAiJbtY9kOfO0LDwgn2C9Y/QnyCyJIH9SvWpclSO3HjtYd7fTblaIo3eqWnjpkKl8f/AbF1J0m1FPvwR9dx6PaCPjbfqiYAn+pgwXxELMHcC8EsK54Hck5yRRbi9Fr9eg0OhwuBzk/5ADwXdl3PLb2MewuO35aP8oaytQFDRQU1pjc41I9LR//+d7dOu2n9cPSaGHOWXPI2Z7DP7f8kwnREwB3ov/vSr8j1D9UXW71o30fccnoS7h+/PVsK93GN4e+4dUtr6rnjYuMY9rQae7A9VjwOjpitPoB7QkUWrca2Jw27E47If4hPp+fioYKrnjjCmqaa6i317PkqyX85+b/cGVc56/jv7b9i8Ehg/nV+b/imW+f4f/e+z/yTfn884Z/9sq39l988gue3/h8p8eE+oVy/YTryUzMZGSEO/XEoepD/PGLP/La1teICY7hH9f9g7tn3K0Gj9vLtnPeq+fhcDnUwKA1z0S7gzUHOVJ3hLX71wKoHzCtAxdzrRlzrZk8U16HdbxpxU0d7usuT6uLTqNDr9Wj1WjVbntPzOEJelr/1qBBq3UPC/AaHnDsQ9RTtnLHSsZGjVWDgLatVZ73Y9vHaFvmObY7+y2NFrVb/EjtEfU9F+QXRIhfCMF+wYT4hxDiF0JMcAwhfiHq7db7PB++AboAAvTu4MFf5155rtHeSKOjkQZ7A+X15Wwu2czWkq3srNhJs7OZkroS9VonRE9gWNgwrI1WbC4buyt2qxOEGuwNnBt7LteNv47x0eO9AlFDoOGEWoIrGirYWb6TwzWH2WvZyz7LPvZZ9lFsLfbqygbULxld4VScNDga0Dl1amtdgD7A67Vv/aPT6Nq9V7RaLXqNXi33vA+D/YK9hqWE6EOICGwJiDwtfnaXXf2CaXfaKasvY8X2FeyocP+fGx0UzZTBUzBGGokJjvFuwVV8t+h6vuBoNVp3C22rISzq3wUaAvQBXsNoPLcVFI7WHeVwzWEOVB+gyFKkztkYGjKUiTETmTRoEpNjJjNt6DTGRo1VW+E7axHuCTanTf3S5QlcPVkYPIsRVDdXU91U3bLdXE1ZfRn7LPvUfdXN1cdtuAL3SlievxvPMKLWQbZeq1e3n7v6uV5dTlyC1H6spLakw8AFYJ9lX5fesB5fHfwK9l8G9YOZ3LvL9XoZHz0e/BxMzryK7QvzoX4I/H0XJLwE19/HP7f8kwBdgJob1eFy4MD9TdUQYCD7+mxuyb2FtcVr0Wq0GCONlNWXsbtiN4NDBjNp0CQ+P/A5q4tW85P3fgLAogsX8fb2tym2FoMCy7cvJzYsloc+eYiHz30YgL9++1d1xatLR1/KIxc8wvrD63n+6ue9ulGbHc3srNjJtqPb2Fbq/nl+w/PqBK2IgAjOGXIO04ZM4/MDn7OjfAcTB01k671b0Wq0zFo2i+9Kv1NTvXgmPATqA/HT+vFd6Xc02Bu4Y9odDA8fzns73+Oqt67i7hl389iFjzE4dDBh/mE0OZqobKxkv3U/+yz7eGnzS1ww4gIa7Y3cPvV2thzdwhvb3sDhcpD1oywaHY08lv8YDfYG9Zv4eSPOI21KWo//h11cVczzG5/nT5f/ibSpaWq3oV6r542tb/BI3iME6gOps9fxvz3/w+60syJ1BZWNlcxaNgun4mTxnMXcN+s+Ne1Ys6OZxV8t5ulvnqbB3kCgPpArxlzBVXFXcfXYq5kQPUG9jtvfvZ18Uz7f3v0tof6hRAdHq0FT6xZOT0tIeX05G8wb2FyymT0Ve9hg3uDuheiGIH2QGmiF+ocS7BfsntwScOz11QWqXah2l/uDvryhnLL6Msrqyzr8AhrmH0aAPoA6W53a0tYZa6OVnfftJDwwnCO1RyiyFAEtLaQKClWNVeSb8skz5RHsF0zq2anMjp3N8LDhDAsb1qVUQk6Xk+3l29lweAMvbnyRQzWHGB42nJLaEqqbqtXlI6OCohgZPpKzDGcxY+gMEoYnMCF6glfrj91pZ1vpNr4++DXbSrcRGxbLmMgxapBV2VyJtclKRUMFuyt3s610GxvNG3EpLmYOn8miixYxxziH/db9uBQXl46+lNGG0TQ7mvnD539g8VeLuX7C9eTcnEOwX3DXX9Q2FEXhUM0hdpbvZGfFzpbfFTvVse+eMcLjosdx/ojzuW3qbQwNHcrgkMEMDhlMeEA4flo//vzln8kuzCYuMo45Y+bQaG/kq0NfUWwt5tqx13JP/D0sWrOIA9UH+OV5vyQ2PNYr8PF8SWv94+mWbn3b00LZYG9gd8VurM1WnC4nwX7BjIwY6ZW+qqqpqt1nzKDgQUyInsDEQROZED0Bu8vO8xufJ1AXyLLrl5FoTOQsw1kn/Jz2FEVR2G/dzwbzBjYc3sC20m2s2L5C/X851D+U8dHjMUYaiYuMIy4yDmOkkdjwWAaHDCYyMLLH/h/01/njH+TfpR6w42l2NHsFtI2ORhrt7i9qDfYG9Uub50ftvWnV89P6d3hAeA9c4YnTKB318QwgycnJAy6Zf21zLeFLwvHX+dP0RJPPP5bXt77OXe/f1a7ck/C+tQBdABeNuog1D7wGNSPp9rsiNtY9JhVOekxqW9VN1fzh8z/wxCVPEBkYSdyMYoq3jQE04FcHt14DY75Sj/d8sHpsnL+R2ctmq7evNF5JRWMF46PHY64xk3FhBj/K+REA58aey/by7Tx24WOs3LGSbaXbOqxXgC6AB2Y/wEWjLuKPX/yRgpICAO6bdR8/m/kzJg/uONJXFIUjtUfcQeux4HXr0a3srtytHjMkZAgxITH8UPYDN068kaEhQwnUB6LT6Ki319PsbKbYWsy6/es4Z/A5NDgaMNeYfa4n7es1Bwjxa2kJ0Wq07CjfQYAuwN1tGTqEPZV7mDF0Bg32Bmqaa9hZsZOUs1N48ZoXOx2f2V0vbXqJBz95kMqFlV7/KZbVlzH2+bHcOvVWnrnyGaKWumf0ltSVMMYwhiC/II7WHeWHn/3QbiLQb9b+hsyvM7G77Pxs5s945spnCPILwuFysNG8kc/3f8620m2UN5RzuOYweyr3qPf11/kzO3Y2l42+jMvOuozzR57fYbBSb6tnxLMjWBC/gMcuekwNIFt32Xv+07e77DTYG9QP+s5+bE6b2hLlp/MjJjhGDVy8tkNiiAmOUbuW2wZzdbY6am21amDRejnJfFM+b373pju1UGAkJqup09dJr9XjdDnbdR0H+wW7A9bQYQwLG8bw0OGEB4RTZ6vjaN1RiqqK2F62nTp7y3jfqYOncvGoiykoKWBb6TYyLsxQh2J4smocrj0MuMdeDg8dTkRgBM3OZkwWE03OJgJ0AYyLHseR2iNYGi3t6huoD2R89Hgmx0zm8rMu57rx1zE8bHi7447WHuW1ra/x901/p6S2hMcvfpz7Zt2HVqtVu5JbdyvbXXZsTps6/rXeXk95fTkVDRWU1Zex17KXXRW72FWxS20tDtQHMiF6ApNiJjFp0LGfmEmMjRrb5Znp6w+v5++b/s6GwxvQa/WcN+I87p5xNxeOcg9jqm2uJfntZDaaN/L81c9zx/Q70Gv16v83uyt3s7dyL/ut+9V8qA32BvXLUJi/e5WmisYKcr7PweFycOnoSwn0C+T9Xe9z94y7eflHL3v1/FgaLZTUlVBSW0JJXQkHrAfYXbmbXRW72F6+Xf2i5Kf1Iy4qTg1eW//uieCsJyiKQkldCd+Xfs93pd+x17IXU5WJoqoiDlYf9GrN9tP6eX2R8OpW1weh1Wi9huG4FBd+Wj9mDp/JPfH3nNZsIqYqE0+sfYJP9n3CxaMuJvXsVK4ae1WfnCzbNl6TILWfqG2uxeFyqH/Meyr3MOFFd7dzw+MNXknrPRZ8sIA3v3vTqzVlWOgwry6u1q7ctY/Vb7sH23f7XXHwIKxYAddcAxMngu7UjnnJ/iKX9MQrwX5sjGdkEdx4GxPiq9lTuUft8lRQeOqKp1j81WK1aydtShqfFn1KZGAk5Q3lPH7R4/xu3e/QarQULCjgoU8eIkAfQHFVMdvLt5P9o2wSjYnc//H97Cjbwf7q/Wo9Xr/hde6YfgeKonDHf+9gxfYVBOgCqLHVED8snqeTnsbmtHFl3JXH/U9p3f51XP7G5bx47Yvc/9H93DzpZjYf2Uxpfak6jgvcgWXC8ARmD5/Np0Wf4lJcfHXXVxiCDCiKQnVzNeX15Xx18Ct++9lvOVx7mCkxU/jVBb9iSMgQNpo38uTnT/L81c/zwLkPqI/f7Ghm+F+Hc8OEGyiuKmbdgXUE64O5cNSFDAkdQpOjiU3mTZhrzei1epYmLuX+2ff3yH+2yTnJHKw+yMPnPax2d5XVl/H29repaKjA9KCJ2PBYZmbPpKCkgDD/MEL8Qzg75myeufIZpg+d7nW+RnsjI58dybCwYZTWlbL+nvV8vPdj8kx5fLb/M2qaawjzD2P60OnEhscSFRhFra2WN797k79d/Tf3c3rwKz4/8DkVDRX4af2YPnQ6M4fPZNKgSYwIH0FseCx6rZ5/bf0XL2x8gVdveBVDoIEmRxOTYyYzdcjUk35eepLnefUEV4H6QFYXreZnH/6MZmczOo2O26bexn2z7iMiMAIFhW8PfcsvPv0F1c3V/Hjyj8m+PpvvSr/jotcuYmniUs4dcS7mGjN7Kvewq2IX+yz7OFRzqNN0S3qtnoRhCcQPi2dk+EiC/YP51ae/4vGLH2fm8JncvepurwwbHp4MBp7hGrFhsVx61qVcNPIizhlyDqMjRqPX6dWWf0OgQZ1YVN5Q7g42LEWYqkyYrCaKq4oprS91t+R2sxW8I3qtnkHBgzBGGr0C0YmDJjI6YvRpGQvYYG8g/X/pvPXdW+4UTRGjOVJ3RM0nrdPoGBkxkuigaAyBBvWLm91p53DNYfZb93tN1gvxC2HSoEn46fz49vC3pE1O44lLniDIL0j9Atb2yxi4P6N++v5PmTV8FmlT09znrt7PAesB9lv3e30ORQVGcVbkWYyOGM1ZhrO8fjwrSnU0LrqjfW2/WHiCRV9jWTvb59m2OW2U1JZQ3lBOZWMlVY1VWBosWJosNNgbvCc8OpvUCYGef2jcvXymKhMXjryQ+fHzCfYL7nSMbUd1rGmuwWQ1YXfaGR89Xh3O0Po4p8vJ0bqjfHP4GwqOFBDiH8K0IdPYb93PoRr3hOFQv1B1optnSIdXYO1y8der/kra1NO3nLwEqf2MJ5XOg588SJ2tDuV37pfr8/2fc9kblwFw5JdHfKaTOfvvZ7O3cq/XylKDQwZTVl/m+8GedJ/70Udh6dIevYxTosHewOQbP2D/RynAsTGUhiK4eR6MLGRU+CgO1hwkJjiGOlud2sJ4buy5bDBv8DpXeEA4aVPSePlHL/O7z37Hkq+WYHO5J000/boJf50/v1/3e17a/BKXjr6UFTtWAO5v0zWLagjUB3LVW1fhp/XjnVveYeX2ldz74b1qK4onmO3MHz7/A3/99q9YMiwMeXoID8x+gJU7VnLhyAv5+7V/50D1AXZV7OKHsh/YdGQT3x76FnNtyzrHEQHuAfMKCpNjJnPn9DsZFzWOT/Z9whOfPUHJr0oYGjqUYc8MY86YObw59812AeYTa57guQ3PccOEG3hv13s8fO7D7KrcRWVDJX46P8ZGjqWysZL/7vovTsXJ7NjZfHbHZyfVJVpSW8LIZ0eqLb1+Wj9C/EJA456cptfqeeT8R1icuJiblt/Ee7ve48+X/5knPnuCL+78gotHX9zunP/c8k/uWeVOWXZP/D28+d2bOFwOLhh5AUnGJJKMScwcPtMraGiwNxCVGUVmYiYPnfcQ4O7q31m+k3X717HxyEY2mTdRVFV03GE0gfpA9ty/Rx032xlFUaiz1Z2ySXUAv177a/785Z997tOgwU/np16TFi0jw0cSEhDCropdjIoYRaIxEUOAQc1N/MHuDyhvKGdc1Dj2V+9XMyL46/yZPnQ658aey3kjzmPakGnsqdxD6spUAvQBnD/ifMYYxlBaX8rB6oMcrD5IVVOVV338df6MCBuhJo0PDwgn3N89IUSn1VFvq+dQzSEOVx/mSJ27BdXTshuoDyRQF6gGKTaXTc2v6eGn9XP3SGjdk+hqm91fXtu2DmvQqGMvPWN2AbX11DPD3nPO4aHDGRkxktjwWIaFDiPYL7jD7vXWM8E7Le/CbbvLro69rWmuoa7ZPZ6x7fWok5c0ehSNot7fM1GvOxPQRO9o20vo0XqMMaBOInMpLjRo1HHenmM8WTGcihOUY+PHjwW5bYP/Z696lrtmtO+RPVXOyCA1ISGB2NhYdW3Y08nhclBnq1NnhHfXweqDjH5utHrbE6Qu/2E5P37HPSN9x893MClmUrv7an7vHYDoNDo0aNothxodFM21sXfw5u3PuB/jRN4RCQnuXKlhYe7u/nHtV2w6VRRFYeIFu9lTOAxsx7qKA6rh8t/CuS+oy6p6tP5Dz7o2i/SP0tV9d8+4G3ONmU+KPsFf509UUBRpU9KwOW0YI438avWv+FnCz3ip4CX1PuvuWEfC8ASGPj2UX1/ya9IT0jnvlfOobKwkUB/IkdojRARGUPpIqTqhw5eZ2TOJi4pjecpyLn/jcoL0QXxa9CkvXfcSCxIWtLvmX6/9NYu/WsylZ13K5iObqbPVqTNpW3f5e6535rCZxA+LJ7swm99d8jvunXUvQ0KGeAWqVY1VTP7HZErqSvj1xb/mj1f80Wddt5RsIT47Xv1PzBPUnYiFeQt5+punuWDkBay9Yy37LPt4ZPUjfLzvY55OehpzrZm3f3ibw788zOWvX84G8wbmTpyLyWrC5rSxef5mr2BTURSmZ02nyFJEgD4Am9PGxaMu5t83/fu43YqJ/0okQB/QbslbT/oeu8uOzWGjpK6EQzWH+GjPR7yw6QVeuu4lzoo4C61WS6O9kTvfv5PZsbN5cPaDXt32NqdNDQqcLic2p43Xtr7G92XfMypiFKMjRqvd+VFBUei0Ou/UNseCitbpltqmXvKUeR63uqkac62Z8IBw/LX+NDma1G53TwuM5zFOxwxwT6vXqXys4z2GZ9KN55oDdAHuD+9jLVcul8trQo6vCTqevyvPa+N5f3iCPb1Wr07gCdIHEewfrE6C02q0XjlvO7qtKAp2pzu7Q7OzWW2p84wz9KRG8vDX+TMsdBgjwkcwKtydC7fZ2czuit3sqdxDSV2JzxbuqKAoJg2axOzY2UyOmdwyYa/NhBqny8k/t/6TD/d+SGxYLFFBURysPuizJVqn0TEkdAgRARHqWOtAXaB7/KXeH3+tv9dEN89kQJfiotHRqC7HWt1cTW1zrdewFV85vv11/upkKU9A3noSlWdmfri/e5KXIdCAIchAdFA0kYGRhAW4e2eC9cFotS2vs69cw63zEbc9rrN9nu3DNYdZ8vUSVvywgkRjIm/d9BbBfsGdPtZXB7/i+pzrOctwFr88/5dcOvpStBot7+58l68Pfc3Oip1UN1UToHcP15oxdAZXj72aOWPm9ItVwXJycsjJycFsNlNQUKCWnxFBam+2pC74YAHLCpepwWV3bTJvYvYrLeMpnb91otVoeWT1IzzzrTuo/OV5v+SZq55pd9+2QWqHjzF/E7NHzERRYPBgKC09gYqewjGp3XFd6lE+yo3B3bKqgWsegHNfBFA/lFq70nglq02rfX5DNUYaMVW1jNEL8w9ThwwABOmCcCgOHr/4cUaEj2DBBwv48qdfkrIihaN1R3n/x++TZEzi7vfvJmd7DvHD4vn6p1/7HH92qPoQo54bxX9u+g9pU9NY8tUSNYfp5vmbSRiewO6K3Tya9yj3xN/D8xueZ02xOyNBwrAErhl7DdeMu4ZZw2fhUlyMfWEsc8bMIT0hnR0VO/j5/35OXFQcdbY6tasH3K2vEwZNYFzUOEZHjGa0YTSDggZR1VTFj6f8uNOJeck5yazbvw6Hy8FD5z7EL8//JTEhMV17oXB/gdtdsZtpL0/DqTi5Y9odFJYU8n3Z9wwPG87PZ/6cmcNnsunIJn7z2W/4xXm/4Ln1zzFz2Ew2lWwicUwi+cX5JI5JZOKgiTQ7mzlYfZAvDnyhBul+Wj8igyIZEjJEnSXcNuhrncapwd5As9OdiN+z71S0MHk++D2PEaAL8Goh89Br9WqOz2C/YDWJuycFjueD2LPtr/NX9/lp/XAqTv7z/X8YFDyIeZPn8fWhr/ny4JecHXM2N068EUOAwWtGr91l54D1AHaXnXd3votWo+WPl/+RUP/Qdt2JOd/n8PG+j/n95b9nZPhIr9n7nmN2lO3g2Q3P8otzf+H1RVpB4YD1AP/Y9A/q7fVcN+46Khsr2V25mw/TPiQ8MFxN+ROoD6TB3kCtrVZN0eSZJAnuv2tDoIGooCj1R6fVYWm08K9t/+LvG/9OsbWYZ658hu/LvqeioYLCkkIO1Rxiduxsll2/jHOGnNMjr6uiKByoPsCGwxvcE3LMG9TVkrQaLWfHnE3CsASmDZnG0NChNDubqWqsUtNQVTVVUdlQibnWzOGaw+3G2UYFRTEifISa2L/tz7iocZ0OKXC6nBRVFVFWX4bdaSc8IJzx0eO73Yq//IflPPjJg1Q0VHBu7LlEBERQa6tlW+k26mx1XDzqYuaMmUODvUG9NmuTtV2+0dY5SD3v+9YhSev/k9V8usqxbn1PhguNerBXi+DJ/O16/sY8WSTCA8KJDIxU31+etGKh/qHqSl6eNGOe7db7OsrPu7poNTe+fSNTBk8h95ZcRkWManeMw+nghY0vsGjNImYMncGzVz+LBo06+99XVgC17NgxnuEITXb38906H7M6tENpn90D4LUbXuPms28+oefxRJyRLam9GaTG/CWGioYK7L+xd2npyeqmaj7c+yG3Tr0VgA92f0Dy28nqfk/X/gWvXsC3h79Vyx+Y/QDnjzifW9+9levHX8/dM+7mxuU3qvtD/UO9kpR76DQ6jtxXz5BB7m9aJ/xu6CNBqsfjua+y+M5rQdHDlY9AwyDIfwpiN8A1v4AhW0DX/mJjw2K5a/pd/OnLP3mVL7t+Ge/sfIdP9n3iVR6gCyA2PBYUMFlNxIbFqmM2o4Oi+eyOz5gwaAIxf4nB0mjh0tGXsu7OdQDsLN/J6qLV3Hz2zXx76Ftuyb2FskfKiAlxHxu9NJqLR11MkjGJ9eb1rNu/Tl0GttHRiMPl4L5Z95GZlNnuOn772W95bv1zHH3kKMF+wVz0z4sYFjaMy0ZfxsOfPExheiFFVUXsqtjF7srd7LPs44D1AOZas9d/7IOCB6nB6/DQ4fjp/NRv9vsq9/HfPf8FUMtGRYwiIiDCnSZGp8fmtHmN1Wr9n+OJ5i30BJAOl4NAfSBNjiaMkUYiAiL4oewHtRt2dPhoLhh1AUH6ILVFxU/n12lewz2Ve3h1y6v89tLfMjh4sFfanrYJzL899C1Lvl7CC9e8QMKwBK/E5ho0pKxModZWS85NOUwZMkW9r07rTo/2zDfP8Niax1gyZwkZF2Wo11dnq2Nn+U6+L/uebw99y5riNe4sE8BNk27i1eRXu9QzU9lQyY9yfsTeyr0Upheyu2I3V751JX++4s8sumjRcccS76rYxfmvnk9cZBy5t+S2m5Hd7Ghm1rJZAHz90699BjpXvHEFtbZaNt6z0evx1u1fx5VvXsns2NmsSF3B8LDh7CzfyeR/TOal614ifWZ6u3P5suHwBh7+9GEUReHO6Xdic9rYZ9nH14e+ZuvRreg0Om4++2YWXrCQGcNmqPdzKS6aHE0nNUylq+xOO9vLt7P5yGYKjhSwuWQz28u2q1+kgv2C3RPXgiLVQCg2LLZdABobHnta6ttV1iYrK7evJM+UR4O9AYAZQ2eQcnYK04ZO6+XatfBMdLM77TQ5mrA2WdU8tOX15ZQ3uCe8ldeXU1pfSmVDJZYmC9VN1dTaatu12mrQqP9fKIrSpRy0bdM6ee6vQYNTcVLTVIMLFwG6AHVcaHeHYbRtsVWvv80Xx+569fpX+Wn8T7t9vxMlQepp8MWBL7jjv3fww89+YOwLYzlad1QdD3g89314H//Y/A/MvzQzPGw42QXZpP+v5T/s9Xev59wR5xL3fBwHrAe8PuhbL3V67dhr+WjfR4D7zTs6YrTXhB9wtzKdO+JcvrnnS1wuCAyExvYTw7umjwWpHgetB5mb8V8Kl6UDOlB0gAZCS+Csz2D0lxBqhrH/Az/3n8KkQZPYWbETcKcJanQ08sLVL5A+Mx3/P7m76z1pgTry0OyHeH3b69icNgaHDGbRRYv4w+d/4EjdEe6ecTeN9kZyfshBQWFwyGCuNF7JmuI1HPmV+zlsdjTz5Lon+cemf2B32ZljnEP80HimD53OTStuUr90FCwoIH6Y93OtKAr/2vYv7nz/Tv5147/4ybSf8Ocv/szirxZz+ZjLOVR9iK33bvVZb8+KWQesBzhQfaDld/UBjtYdxely0mBvoKS2RA0G9Rr3f6rDw4dTVl/mNV4zMjCSOWPmMC56XLtlPutsdSzKX6QOP/nrlX9lduxsn8uJbji8gUfyHuHxix7n/6b/H3W2Oi59/VJ+KP1BnfE9KnwUP5T/AHTcu3A8ZfVlDHl6CMtTlnPL5Fs6PK7B3sDEFycyY9gM3v+x7xXciquKueqtq9hr2cv0odN59qpnueysy9hweAO3v3c7+yz7ePyix/nTFX86bsBorjGzumg1v1r9K/RaPXfPuJtbJt9Cna0Oc62ZuRPnenXprdu/jgUfLKCqqYqPb/uYKYOnMPWlqYwIH8Ha/1vb5cluW49uJTknmfKGcm6edDMLL1zo1er4fen3XPTaRUyOmczylOVeY3DXFq9lzr/m8M4t73DTpJbcsXsq93DeK+cRPyyej277yGsITNo7aXxz6BuKHixq96W+zlbH/A/ms+HwBs4dcS5l9WWsLV7L5JjJNNgbOFh9EL1Wz5jIMcyOnc0loy7hR+N/1KNZKHqKoijU2+u9crqKvsmToaLtT0ltCUfrj/2uO4q1yeo16cizQpcaoB7LP9t26WJA7cXxtThFiF8IUUFR7uVwdQH46/3VISRB+iB3gv5WCwx4vgx7vjS37mVpvSJY67LWx7a9PSRkiM+J2aeKBKmnQUJ2AoUlhTw4+0He3v42ZfVlhAeEc/+s+/nzHN+TFzzmr5rPK1te4d83/Zu5E+ey9OulPPn5k9w69Vb+8/1/WJGygtTJqUQsiSAyMJID1QcAGBc1jr2Wvep5jhdEnR97Pt+Xf8/cwb/kzXvca5i7XHDCE7X7aJDa1rp1Cok/3oWzOgZc/mAPORa4KhC5B6a/CcO2QPQOiNqvdiUNChqE+VdmIhZH0ORs8gpkPS4YcQHfHP5Gvd36NfA1nOCy0Zfx5k1vkroilfXm9WjQcNnoy7jp7JtYuWMlXx/8mrQpaSy6aBETYyai1WhpcjQRtjiMISFDuGv6Xe3GjJbVl3H7u7eTZ8pjWOgwxkePZ92d6zBVmYh73p25oe2s/u7YdnQbl79xOaMiRqmTXsIDwokKjCLYP5gPb/2QEL8Qdab3a1tf48uDXzI8bDgXjbqIt29+W53U8qOcH/HlgS+xu+zEBMew/+H9x11EwKW4MFWZ2FKyhU3mTazYsYJDNYe8/mP/0+V/4olLnjih6wMw/s3IteOu5cVrX+zwmEX5i/jr+r+y4+c7Ol0fvcHewEd7P+JvG/7GVwe/YnbsbApLCkkYlsDLP3q5XVaC4zlgPcAz3z7Dv7b9y2sc4Mjwkbya/CpJcUl8uOdDblpxEzOHz+S1G15jfPR4Hvz4QV7e/DLb7t3mc/x6Z2qaa3hx44u8tvU1iquK+cV5v+A3l/5GTRW20byRucvnUt1Uzbo71zFz+EwqG9z5a0eEj2DdnevU19XaZGX2stnotDq+vfvbdi3CnrHO797yLnMnzVXLLY0WrnjjCoqqirgq7ioqGiqICIzghgk38JNzftLpcqdC9BeKorC7cjcVDRUE+wUzcdDEPtV6fjpIkHqKOVwOIpZEqN0fbXlaVNcfXk9EQES7D4zH8h8j82t39+1d0+8iQBfAevN6ChcUErEkgoUXLuTOaXcy+m+jOTf2XGbHzuZvG/5GxgUZZH7Tvtv38Yse56mvnmpX/tkdn3H5G5fDHxvAGcTQoVDiOzNV1/zsZ/Dyy3DRRfDuuxDT9bGJvenhP+zmH0uHYG8KcAetGge4PC1Sivv2pHdhzGegb4Dhm2DwrnaTscB7zKqvILY1DRr0Wj03TLiBd3e96xVkeVZAak2LFn+9u8WlydHkswX1k32f8NP3f4pTcTJz+Ew+3vsxCgp7H9jL2KixvLH1DVabVpP9o2xC/EPajWHaUb6DwSGDmTZ0GqH+oe1aeA5YD3DuK+cSExxDsH8wuyt2c9Okm3ht62s8fN7D/HfXf2m0N5L3kzw1/ZJLcfHezvfI+SGHd3a+w79v+je3Tr2V3B25pK5MZXjocGqaa7hz+p28cO0LPp+rktoSXt/6Omv3r2X94fXqsJUhIUPUZNWTBk0iflg8D537ELNiZ3X4vHfFAx89wPu73+fAwwd8tjiu2r2KG96+gcVzFvPYRY+12+90OXk071FW7V6ltq60nhQxJWYKWT/KYkTEiBOuY4O9gW1Ht+Gv8yfIL4iHP3mYfFM+EwdNZGfFTm6YcAMrUlfgr/NnjWkNiW8m8vdr/87PZ/38hB/T5rTxzDfP8Icv/kCALoCfnPMTrhhzBdeOu5ZmZzNJbyZxpPYIf73yrzz11VMcrjnMxns2MiZyjHqO+z+6n39t+xdb0rd0GNyf/+r52Jw21t+9Hj+dH02OJq5880p2lO8g7yd5OFwOtpVuI9Q/lItHXewebiNEH9Z2RbjWKaY6SoHVlWM7Sld1Iufq6NizY84+rXlsJUg9xVrnL/XFMwnGM6mp9YSqioYKfr3212QVZKlliWMS8dP58dFtH3HuK+dysPqgmu/u9qm3syx5GUF/DuLJS5/kyc+fVNeb9njl+le454N7vOqg1+jJuCiDJc+X43z/ZUBz4mNRB5jyqkYuvrqUPd8ZUJrCcUejnifnWCuffy1E7wH/OrAYYfz/4Ef3gUaB5iAIaOSiERfx16v+yp3v36kuBdhaeEB4u6UPj8czezVYH8wbc99Qlz1ttDfy+89/z9Kvl3LRqItYnrIcl+JixLMj3Omc/EMwRhp548Y3mDJ4Coqi8NAnD/HCRt9BoYdOo1NXt9JqtJTUlqitwYZAA2+nvM0VZ11B5NJI6m31fHv3tyS/nUx5fTmvJL/CndPv9DrfZa9fRpBfEC9f9zLnv3o+UwZPUZcX/eS2T7hq7FU+6zHnX3NYf3g9l511GReNvIijdUfJM+VR3lBOs6OZq8dezfrD6yl+qLhH8lDmm/JJejNJHVrTWll9GeNeGEeiMZHc1Nx2Qazdaeenq35Kzvc5LEhYwKiIUeokh9K6UvZY9rD5yGamD53Ol3d92WEd1u1fR5OjSV0KNCooyr00aAcT2RwuB29ue5NNRzZx3ojzuG3qbers8Etev4RmRzMb7tnQIzltzTVm/rbhb6zYvoID1QcwRhp5a+5bjDaMJvFfieys2Mn46PHkpuZ65Yr9357/cX3O9Tx31XOdZoPYZN7EBf+8gHFR45gxdAafH/ic0vpSJsdMZk/lHq/MFVFBUbw37z0uGX3JSV/XmajJ0cThmsMcqj7k/l1ziKN1R93jLDsKXHowAOqtwOt0P25/9tGtH3HNuGtO2+O1jddkWdQetrPc3Xp22VmXsW7/unb7zbVmEkjwed8f/edH7fJ3mqpMXHbWZYB7hZaN5o3qvktGX+IOWvyC1W7/tvkb2wao4J5E9dS6v6C83whoyMhod0j3PfaYO7nqddfBG29AVFQPnPT0i4kMYteGs9qUanjhtcM8/riGuho9OAKg/Gxw+oOihYKfQel0GF4A390OTRF8hcJsrQP034BfI+iaIdIEoUdh2hvUjF3tjn81Ha8GBe4PYc8MXz+tH5NjJvPt4W+56q2rKH+knBU7VvCLT3+hvu6Haw7T5GgiNixWHbdqbbJSWFLIRf+8iJd/9DKfFX9GdmF2u8e6edLNTB08lU+LPuXbw9/y81k/Z2T4SDaaN5K7M1c9LkAXQGVjJTe+fSPXj7+euMg4thzd4pWF4qfv/xRDoIEbJ97I1qNb1aDr8wOfE/e3OIaEDuHc2HNZU7yGYaHDuHzM5T6vv7qpms/3f84L17xAdFA0v/7s1+y17CVlUgr7q/erAd/KHSvJN+V3GOh2x+VnXc4YwxieXf8sb930ljo2srSulP977/9odjQzadAknvryKVyKiymDpzB30lyKLEXc/t7tbDJv4t83/Zt5U+b5PH/O9znc+u6tFFcVe7Uyeny671Ou/vfV7cr9tH784fI/EB0UzYHqA/zu0t+p3dx6rZ67ZtzVLp/huv3r+OrgV/wv7X89tsJNbHgsS5OWsjRpKT+U/cCCDxZw8WsX88DsB/j09k+pt9d7zTKvaqzi39//m0dWP8INE25oN9TEk2x/b+Ve9lTu4bvS7zhnyDl8X/q92hsxOHgw46LH8X/T/o8LRl7AjKEzsDZZSXsnjcR/JfLyj17mrul3ndZVfHpSZUMlG80bOctwFnFRcSc9TrXZ0awun1taX0ppXSml9aWYqkzsrNjJ4erDlDeUq3mcPQyBBkaEu7/ctp6E01Fqpbbpkrqy3zNG80TOpaXrj9OTde7KsUVVRbyx7Y0OlyvuTHhAOI9e8Cjjo8ef0jqeyHPTlVzPp5IEqT3s60NfE+oXyjVjr2kXpOo0Oo7UHlGTR3s4XA7yivLaBagAVU1V6oSrqYO9V7DxvHmigqIoriruch0VFJS/HAG06PWwZEmX79qxN990pwb43/9g//5+G6R25IG7RvCAj3zGLsXFG2+X8Zsn4yjZdg4ufS0EusAZAA5/sIUdy92qgZpj+W63p0FADUx7A87OxTn6S3Dp3fepHQbR++DARTBoFxZbPfgD1cOprxvCt7b1aNChaJzEPO09pEKDhmJrMRNenMDQ0KHqkJPIgEiqmquobq4m7Z2WPMFxkXHst+4nOjgaa5OV/+76L3mmPHQaHUH6IDYc3sC0hGm8u+tdABKGJnDd+Ot4bsNzNDvdS0G+vf1tJkRPQKvREqQP4s25b/LS5pfYaN7Izctv5vIxl6upsjycODlSd4S/b/o7iqLwy/N/ib/OH4fLwT7LPs4ynKWm6fq06FOcipNiazH3fXSf2irxwZ4PaHa6g8UXNr7A5JjJvLrlVSYMmsDwsOHd+pDfZN7EW9+9pbaO1DTX4KfzY/n25ZhrzDx52ZO8u/NdXt3yKs3OZoL0QWQXZKspzcobyvn9pb/nT1/+ieFhw/nyri85f+T5HT7e9ROuJ9gvmNe3vs7vL/99u/dTRn4GF468kLdT3qbeVk+trRZLo4XXtr7GojWL1MdNGJbgNW4T3AHf7z//Pf/5/j/85JyfkLszl5nDZ3LtuGu7/Hx0x5TBU/j8zs9Z+vVS/vLNX8gqyOLns37OteOuZUf5Dv6x6R9qoJk2JY3bz7mdVwtfpdhaTLG1mL2Ve9lr2evVqzAuahzThk7jpok3MX3odHVVsLaGhA7hk9s/4b4P7+PuVXezumg1ryS/Qqh/6Cm51lOhsqGS59Y/x3MbnlOHseg0OsZEjmF89HjGR41nwqAJjI8eT3RQtDtN1bE0XJ7UVZZGC2UNZZTWlapBqa9AKdQ/lAZ7gzq0KFAfyKDgQYT4hRAWEEZlQyXNzmZW/XiVzy9PA53T5SS7IJsfyn5Ao9HgdDnV9Fh2l51pQ6bxwOwH2nV919vqeW79c/x909+5aNRFPHbhYx0uc+vrS5TT5eSJtU/w1JdP8Z+b/8ONE288FZfXb0l3fw8b/JfBNDoaeeX6V9Rk+x4jwkdw1/S7uGXyLUx9yR1wKr9T2uUznTtxLvW2elabVqPT6Hju6ue4f/b9rN63mqv+3dJStDV9K9HB0Uz+x2RC/UM5UnvkuPUL8QtB85dS6qzubsMee/X7ycSp3vLJNwdZ/MdAthbqqbMG4HJqwaUDjQtCSkHfDBf+1T3m9avHYEcKoAGNE+bdDCO+hSALbHzAnUprwiqY9F/3ogX2IAgpg6HbwL8ediXD2e+5H/jDv4F5tntIwpQVEGUCrROOTnFv65sIbTibuuCd7rpoWpaTbU2n0fHMlc9wzpBzWLljJaMiRvHlgS/Vda3PGXIOBSUFBOmDGBQ8yCsXK0CwPpgGhzto1uAel/l9+feA++/CX+dPeX05tbZaxhjG8P6P32fqkKlc/sblbDi8AZvThlNx8tkdnwHw313/JWFYAlMGTyE+O5758fP555Z/otFomDl8JitTV9Jgb+C9ne+xrHAZlY2V2B12dNpjrTjH8izanXaqm6vVRRD0Wj3RwdGMDB/JD2U/qKshRQVF0WhvZHz0eDbcs0GdSa8oCle/dTWfH/icmcNn8untn3aaW9bjwY8f5N/f/5vih4rVCUgA//7u39z+3u18/dOvuWDkBV73URSF3RW72XRkE09/+zQjw0fyv1v/57X/d+t+xx+/+CPTh05nR/kObE4bE6IncNOkm7h09KVMGzpNXcDB5rSxcvtKPi36lP3W/RyuOYzdZefS0ZcyPGy4Ols41D+UiMAI98pPAeFEBLi3Y0JiiA6KVj94qxqreHb9s/xtw9+oaa5Bq9Eyd+Jc5k6c615LfuPzgDtd2YjwEZxlOIuxkWMZFz2OcVHjGBc9jrFRY09oosjyH5Zzzwf3MCVmCuePPB+X4k7nc8PEG9o9j6eC0+Xko70fUVZfpqZZ8/yoS2UeWy6z3lbP4ZrDHKg+QGldKRqNhuigaDWI1Gg0TBsyDb1Wz17LXoosRT57WsL8w9S0VYNDBjMkdAiDg4/9DhnMkJAh6naofyjTXp6GMdLIs1c9yxjDGCICI7zOZ2m0MGvZLLQaLV/c+YXPVQwHql0Vu7jr/bvYcHiDOkRFp9G5Z87rg9BqtHxz6BucipNrx13LOYPPQavRsrtyNx/v+5ja5loeOvchnprz1AlN4mu0N3LHf+9g5Y6V3D/rfhYnLu5XX7Z6knT3n2Jjo8by7eFvvbpHLxhxAT+f9XOe3/g85hp3eh+Ptq2qABGBEeoqIk7FyZcHv+S2qbe1WzUiUB/I7z77HTXNNe4PBR8TbtqK2baE/ccC1M8+O+HLFN109QWjuPrjjvZ6VhR7WS1xOBRWrHTxzzdr+WH9K1RVBmBzNruHx2qbYdeNsOsGd8CpcblbYl3Hxswqehj1NZzzFmy5BxzBcOQ8WHcs9+uMVyFpIRy6AL58nLr9rbvaj4WnBhMMK4SGGLCMxVk3jId/73I/Fue7x9/ygPu3xkmBxgWxG2m84BkOhZRD2aVwJN59jrGraQgtg6Zw+PgZlIkf8n1xFQweBP51HN41CIKsENEAWiguP8o5f0gFNBBVhE7jT6hiJDTIjxKLFVPtLtYWf8b2su1cMPICzh50Nqt2r0LBvXLQnso9jHy28y4qT7ejQ3EQGRjJbVNvo8HewOHaw3x54Ev2W/cTpG9Ju2JptJBkTGJ5ynL179CluCg4UsDgkME0O5u5dPSlXgGqS3FxqPoQn+3/DJvTRmldKeUN5YT4uZOD19vqueHtG3j9htcZbRhNg72B33z2G26Y0BJYOV1OdpTvYOqQqWg0Gp7f+DwvbX6JAF0A35V+x3PfPse9s+6lvL6cO967g88OfMbciXOJHxbP3sq9TIiewFmGs3il8BUWf7UYOBbcBEZS3lBOo6NR/VIS4heCzWnj7R/eVlsumx3N6uQ6XwyBBsZHj2dc1DiGhg4lKiiKP13+JxwuB6MjRjMyYiQFRwp4fuPzPHHxE9w1/S5GRozsckv3ropdrNu/ju1l2ymtL2VE+AhmDZ/FuSPOZYxhjBogz5syD41Gw49zf8ymI5s4y3AWjY5Gln6zlOwfZWOMNLL16Fa2HN3ClqNbqGyo5JELHuGRCx7pUj06Y3faufY/15Jvygdall31pPPxfCFyuBw0OZqobf7/9u4+uI3yzgP4V7Id20nsrCUHEgw0kSkEAg1ISfrCAb1G4trpEW6mUjy56fXaGyJx7QwtHSrVf9BeS7mM3DKU+6MzEp2WY3o1YLflrXeARKcUrgFiqaXHe9GGl9gJBMtrO/FL9LL3x2rXerVlaWXL9vcz47G1u3r07D56rJ+efV4mIUNG27o2XLXlKlzSeQk6Wjq0hQteH30dv37t17j1E7fi8X98HIlUAsekY5BmJG0+VaFFKGve7bScxqunXsUdz9yhDWqbTc7itQ9fK9qn8gd/+wPc8j+3YO9P9+I7134HV5x7BWxbbRUFXuqiBOrqVdk/6i3//B/1NnM1/jr6V9x15C60NrbC0mHBJZ2XYOfmnUVb498bfw+BSAA/+uOPcOGmC/HcvzxX8kvNydMncf9L9+PRNx7FT979iXYr/GbbzbjJelNVrc+tTa14wPkArnnxGnjDXjz06kO47ZO34ctXfnlRC6ToTZqRsKl507J2o2FLqs5kWYa5z5yzHvVtn7oNL518Cal0Cusa1+GGi2/A1/77awCAZ7/yLK75ee7a49YtVhyfPI4PznygbbvvxvvQaGzEF3/zRW3btz71Lfzwjz8sO28dk1dj7D+eBBIb8IUvAIODCz+nbGxJrVvvvQcM/mYWT/3hA7wRm0H8zBRmUtNInjYjPdsMOW0EZjYpfWybM9MapZuULgkzwtxSs0hnglMAyKz2ok1zUOyfWGZb0xRw0RPAyV3AWOnpmtA4o7QYb/8dYHka6HwdiDmA3/4EmDZnDir335Ws/BiyH6sr1MiZVuM0IBsAGJVZHJC1D7LS31g7VgaQmnueAZnfKSWNdCO0QXaG9NyXh40nlT7JhpSSXroBMCaVc208m8kDgNNblb7O6QYg1aJsWzelpGdMAKlWJY8tcaApAZzzZ6BpVtlmkJUvAK2jymsZk4AxrTyvIQkYzyqvn2pWXiPVlMm7AUZDo7Kizro2NDYZkcYsTp45jsRsCz513jXYvGEzDE3TmJUn0dA8BWNzAqkGCee2t2NdSyPiiRF8MH0c708NYzz1Ic4kJUwkJMiGs1lloPxsbN6ItuYNOJ2YwKWbL8Eh2yHs2nIFNrasx0xyGqNTo2hqaEIylcQjbz6CX736K4ycHkGDoQEXmy9G5/pOvCO9g3cn3gWgBMg7Ondgh3kHtndsR//L/RifHsfG5o14K/4Wtndsx/DEsLaijtFgRNu6Nmxq3oSUnMLw5DC+ffW3ceOOG3FB+wXYsnELDAYD4tNxnDpzCps3bEbn+s6i7650Oo33z7yPoyNH8b1nvofoiSg+ft7H0dbShpHJEbw3/l7O6nQGGHCR6SLs2rILV59/Nd4efxsvHH8Bn972aXSu78y5Bd/S2IKnxafx4KsP4vZrbseVW6/Ex875GDZv2Iz25nbtVvSJ0ydwfvvcDBGpdArx6TgmZieQTCdx8vRJ9D7dm7PoSyUuP+dyPNLzCJJyUlnvPZ3CNmFbwfRhsizjlVOv4PE3H8fjbz6OI8ePVLTiU6kAtsHYoN3tyP5pamjSAtyRyRG8Lb1dNF2jwYjmhmZthbbZ5CzGZ8fRZGyC7TwbPrPtM2htai0ImidmJzA+M47x2czPzDgS6QQ2NG1Q0jKuU1aYalJWmmpvbtdWnFIHPm5o2qA8zvy9Yd0GtDa2Fg3+3pHewfef+T5+8X+/QDKVhHWrFXu79mJr21a0rWtT+gw3NGEqMYULN12IHZ070NHSoWsgOTY9hruO3IV7XrgHD/c8jH2WfbqlvRCO7q8xdVAEAOzcvBOxeAwzKaVVdO95ezGdnMaVW67U+sDd89l78PUnFl73/Pru67Fv+z74wsoop+aG5rnlzIpoNjZjNj23v+3tHkze1w/AAKtViSN1tWsX8Je/AOvXA9EocEnpGQ5o9ZmYnci5bQ0AE1NTiL07hdhbDXj55SReeyeOI8ePYGxiBtsm/gnxmTimp1M4e7oNyaQRcqIVKUwj1fIB0DQFw8Q2yM1jSmB1dr0StMkNSpA1uykTYMpKYLbutPI3GpQW5cR65VjZoARJskGZYixHJpgElACyIPDOVskHwMocxFO9lfKRkvWFq2iei23LOtaQd6gh/zn5B2R+q8ep7838Ywz5z8nPU9Y6oDBmnUaR1wAyAzwrfC8aZKBRGWA7l0ZmxGf2l5Ds7GYHS0XfCvl5UdOQc89Bff38hNS6mp94Q0L5gpZoAZItc8/XvjSmMr+Tym/ZqBynfhnNzkfONcz/O/MjGzNfTvPOpdi5FlwH9Qtt/k4D5hoDsl6r2IU0pJRzBjJlDMCY/R6QlbQA5UurIalcH/V1ZGPmHIyZ/5UNuX9n8n7hlW/inf/9m8LXrxHe7q8hWZbxsz//THv8pV1fwlf3fBW3PXUbApEAWhpb8MboGxiZHIHb5kYgEsgJUIUWoeTIwKdiT+H4+HHtcXtzO05NnSqZFzVAvaD9Arx3rBmT//lfAAz46EdrEKACwEsv1SBRWinyA1QAaF+/HlftWI+rdgDOvweAzQCyv7wU67vZAiB7YEKpv7O1YTa5DpETEYzPjMO61YpzNy7cLxRogCwbcXziOM5vP7+ilghZljGbmkWTce52aPY0WKenZzE+kcbkRAPGxmSMxWWMjQPj40mMT6YxPiEjPbMeI5Pv49j4Wzg1NYyZZAIbkhegffZSTMxOYursaUwnZnE2ASTONkA+uxGNqXYY0y0wyDLONkwgYTiNtOEskGqGMdEOJFogpxsBOdN9Q25UWmnlzMed3ACkGyAbknOtvFpQn4LcOJX50E/MtTan1eDEiJxWddmQea4a6Gc++LJblpF5mGpU9huzW6WVLiPacalmYLatSACgMuT9LSsfvg1nM63GMpBsnQtUSpee0ofbkJ5r5VY/uNUAUst7VtCUVq+BmgW15T1z/WCY26bFUpntaqCnXe8FsqjmU/2dE9TkBWran9mB0fyLY5T10omNQHY3soLYu5ZfxkqlLS+wnypTeD3fe71wOfWltCZaUm02G7q6unDw4EEcPHhw4SdUIZFK4Pnjz+O6+67Dj//ux7jlE7cAAK75+TVIpVParZfnvvIcnow9iTv+oKwYdMPFN+DE6RMYGhkqSNN5mRNPvPWENvqzXE3GJqTuOIN0QrkVed11wO9/X9XpERERUYXSchrSjIRTZ06hpbEF57efX9H8zrIsI5lOluwvnEwntfEq7c3tEFoEpNIpjM+OIz6lzBCRTKUyXRXa0NbchvWNGzCbmsXYzBgmZk5DlmVYTB+BsL6t2tNeUH9/P/r7+zE8PIxIVkvammhJ7erqWrLb/U0NTdh5zk50ru+EN+zFJy/4JPZ07UFXWxde//B17bhdW3Zhb9deLUjdt30fHnvzsYL0LjFfogwQqOAbY+LRe4CE8gbu7QX+vXDhKf2wTyoREdG8jAYjTK0mmFqrm6bRYDDMO6Ct0dhY8DrGBiM613eW7GsNABvQDNOGwjtjtaY2Iu7fvz9n+5oIUpfa3UfuxumzpzGbmsWND9yId77xDrrauvD88ee1Yy68+8KcFXnUVX3yuS5zob25PacT/nwcFgeeuf9anBU/AYh2AMCpU0Bn6fckERERUd1hkKqzU2dOIRAJwG6xw2gw4pE3HsFnf/FZdG7o1FaFAoDZ1Czufv5u7fFNj92Ej2z6iPZ473l78eLIi/jtX3+LP538U1mvffHvjyD0bx+H2q/kgguAd9/V57yIiIiIlhKDVB3Jsox99+/D5NlJPPbmY9o8i797+3doNORe6ust1yNyIpIz6Xl2EPviyItob27HS++XHpBkNBjRe3UvnviuF5FIEm/OzDXrHzsGbNum04kRERERLTEGqTpb17BOm4i/o7UD8pSMmdQMknIy57iH33g453EDGpBC7qoi2UsFFpAB4+CvcOd3/yFn8113Ad/8ZsXZJyIiIqoLKy5I7evrw+joKMxmM2KxGBwOB5xO53JnC4DSkfki00WInFBGpnW1dRVdqrTYspP5AWopW48GcSJiBd7/GJLyXKfpYBA4dKiKzFfrpz8Fbr8duPlm4LLLljEjREREtBqsqCDV4/Ggu7sbfr9f2+ZwOBCPx+F2u5cxZ3P2bd+HkckRPPvuszg6crToMfkB6nyM6Wak//gN4GUX8MEVOJE1IXlzs7Ka0OblWzVtzuc+p/wQERER6WDFBKnRaBTBYBD507r6/X7YbLa6CVJvffLWkutclyXZBLziAmY6gJf+GemR3cieYNdgAO6+G/j6wotULa1rrwWefVaJmF94Adhe+TrGRERERCsmSA0EArAWmXtT3TY4OLjst/3VyXXL0WRsQmKqGTh1OXDyciB6CBi7SFkrHUZkL3/X1gYcOQLs3FmrnOsgFlN+nzoFjI0xSCUiIqKqrJggNRwOFw1SAUAQBIRCoWUPUgHg0s5L8eqHrwIycDaZAM6cC0jbgPCdwPu7lKX62k4gMdMBTJszz8pd4q29HfjlLw34/OeX4wyIiIiIlt+KCVJFUYTdbi+6z2QyYWiocDnRpWYwGPDnf42qj0ofOHaR9qfZDNx5pwEeT23zRkRERLSSrIggVZKkefcLgrDgMUsnNzhtbATOOQc4eBC44w6gtXWZskVERES0gqyIILVaw8PDOevBqmvE1oJc/sD91eXaa4EHHlA6zm7duty5ISIiojrX39+P/v5+7fHw8HDO/hURpAqCMO/+hVpRu7q68Oijj+qXISrU36/8EBEREZUhv9Ewu0ERUIaRr3jxeHzBQJaIiIiIVo4VE6QKgoB4PF50nyRJ2L179xLniIiIiIhqZcUEqQcOHIAoiiX3OxyOJcwNEREREdXSiglSXS4XotFoQf/TcDgMACWnpyIiIiKilWfFBKl2ux1OpxOHDx/O2e73+zEwMMA+qURERESryIoY3a8aGBhAX18ffD4fzGYzYrEYPB5PXaw0RURERET6WVFBKgB4vd7lzgIRERER1diKud1PRERERGsHg1QgZ7UDWn4sj/rBsqgvLI/6wvKoHyyL+qJXeTBIBd/c9YblUT9YFvWF5VFfWB71g2VRXxik1jG9K8taS09veuav3q/dWiqLtZie3ur9fOs9PT3V+7nWe3p6q/fzrff09LImgtTh4WHs379/yQqh3t889Z6e3hik1o96P996T09v9X6+9Z6enur9XOs9Pb3V+/nWe3qVvP7+/fsxPDycs90gy7K8THlaMjt37kR3d3fJ/cPDw+jq6tLt9Zje6k2vnvPG9Jge01u56dVz3pge01uq9GKxGF555RXt8ZoIUomIiIhoZVkTt/uJiIiIaGVhkEpEREREdYdBKhERERHVnRW3LGq5+vr6MDo6CrPZjFgsBofDAafTueRpkKLaa+lwOGC1WtHT0wOr1QpRFBEIBCBJEgKBQA1zvrpIkoRDhw6hp6enovcy64S+qikP1gl9SZKEw4cPQ5IkiKKIeDyO3t5efm4sAz3KgvVDX2qZmM1mjI6OIhqNwuFwLHqp+kXXD3kVcrvdst/vz9lmt9vlQCCwpGmQQo9rabVaZQA5P3a7Xe+srlpOp1N2u92y2+2WAcgDAwOLToN1Qj96lAfrhH7GxsZkt9stj42NadsikYgMQHY6nWWnwzpSPb3KgvVDP2qZ5G8TBKHm9WPVBanqm7nc7bVKgxR6XUun0yn7/X7Z6XTKXq9XDoVCemZzzYjFYhUFRawTtVFpecgy64SevF5vTlCk8vv9MoCyri3riD70KAtZZv3QUyAQKHrtnU6nDECOxWILplFp/Vh1NcftdstWq7XovnI/DPRIgxR6Xcv8b3FUmUqDItaJ2qgmSGWd0I/FYpEtFkvBdvUDtJxrzTqiDz3KQpZZP/QUiURkQRAKglT1TlAkElkwjUrrx6obOBUOh2GxWIruEwQBoVBoSdIgBa/l6sBypNXMYrEgHo8XbBcEAQCK7svHOqIPPcqC9GW1WjE2Nga73Z6zXX3PW63WBdOotH6suiBVFEWYTKai+0wmE4aGhpYkDVLofS2j0SiCwSCi0age2aMysU7UL9aJ6oVCIYyNjRVsV6/pnj17FkyDdUQfepRF/vNYP/Q3ODiIeDyOgYGBso6vtH6sqiBVkqR59wuCsOAxeqRBCj2vZTweh8/nQzweh9vtRjweh81mgyiK1WeU5sU6UZ9YJ2ovEAhAEAS43e55j2Mdqb1yy0LF+qE/URTR19cHl8uFQCCAY8eOldWKWk39WLVTUNHq4nA4cv452e129PT0wOFwIBaLLWPOiJYH60RthcNhhMNhDAwMaLeaaXlUUhasH/qzWCzwer2QJAnBYBA2mw0DAwNlBaqVWlUtqQu9ecv5JqtHGqTQ81oW+/Zst9shiiIGBwcXmTNaDNaJ+sQ6UVtqa1E5c3OyjtTWYspCxfpRO4IgwOv1wm63l9U6XU39WFVB6kLi8XjV34j1SIMU1V5LtRM2ByQsL9aJ+sE6oQ+Xy4Xe3t6yby0vhHWkcnqWBeuHvlwuFwDA5/NVlc589WPVBamCIJQc/SdJEnbv3r0kaZBCj2vpcrlgs9lK7udoz9pjnagvrBO14/P5sGfPnkWvpMM6or9Ky4L1Q182mw0Oh6NguzoQqpxBaZXWj1UXpB44cGDepudiF7oWaZBCj2sZjUaLjgpU3/CLHe1Ji8c6UV9YJ2ojGAzCbDYXBEXBYHDB57KO6KuasmD90Fc0Gi363lavZzl9UiutH6suSHW5XIhGowV9HMLhMAAUzPNVqzRIoce1dDqdRW/PqP2K9LolR6WxTtQX1gn9hcNhSJJUtNWunD6lrCP6qbYsWD/0ZbfbEYlECrar17inp2fBNCquH4tfe6D+qcugZbPb7UVXNCi1usVi0qD5VVsexdYNVlfAYHksjrpqy3xrJbNOLJ1Ky4N1Ql+xWEy2WCyy2+2WvV6v7PV6ZbfbLbvd7qLvcdaR2tGjLFg/9BWLxWSn05mzXK26Wl6xlb30rB8GWZblRQTUK0ZfXx9GR0dhNpsRi8XgcDiKjgy02WwwmUxFv3WVmwYtrNrykCQJPp8PgiBotwx6e3trOvXFauLz+SCKonbbRhAE2O12mEwmBAKBnGNZJ2pPj/JgndBPd3f3vLciI5FIznVlHakdvcqC9UNf2ddTkiSIogiPx1PzuGrVBqlEREREtHKtuj6pRERERLTyMUglIiIiorrDIJWIiIiI6g6DVCIiIiKqOwxSiYiIiKjuMEglIiIiorrDIJWIiIiI6g6DVCIiIiKqOwxSiYhoQfOtAlRMNBqtUU6IaK1gkEpERPOKRqMFy7WWw+fz1SA3RLRWMEglojXN5XLBZrPBYDDAYDDA5XLl/DgcDng8HkiSVJPXlyQJ3d3d6Ovr0zXdcDgMm82Gjo4OuFyuitNR1+z2+/2Lep7VaoXZbEYwGKz4tYlobWtc7gwQES2ngYEBAIDBYIDVatUeZ/N4POjo6EAoFILdbtf19ePxOERRxNGjR3VN1263IxKJwGazVZWOy+WqqBUVALxeL2w2G+x2OywWS1X5IKK1hy2pREQLCAQCsFgscLlcureoWiwWyLJcNDjWg8lkqvi54XAYAKoKMHt7e3nbn4gqwiCViKgMdrsdkiRpgdtaEAgEqg4wnU4nwuFwzbpLENHqxSCViIiKCofDunRvsNvt7JtKRIvGPqlERGVQW1CzgzZ1UFF3dzdGR0chiiJ6e3thtVq15/h8PoiiCLfbjZ6eHoTDYYRCIXg8HlitVng8HgwNDcFisSASiRS8bvaAqtHRUXR3d8PtdhfNY19fH2KxGLq7uyEIQlW36QcHB0s+XxRFBAIBdHd3ay2ksVhMO6d8DocDAwMD8Hq9FeeHiNYeBqlERAvo6+uDKIoIhUIQBAGAEqjZbDY8/fTTWmCWv00dvKQGc+FwGF6vF4FAAA8++CCcTidCoRAcDgfi8XjB69psNvT29sLpdGrbfD4fXC5XQR9WdYBS9iCncDiMoaGhilpDjx49WvJ5LperIKB2OBwl09q9ezc8Hs+i80BEaxuDVCKiDFEUC/pgiqIIk8mEWCyW07Lo8Xhgt9tzWg4tFgvsdjt8Ph9CoVDO9oceegjHjh0DAEQiES3YVffnB6kejwcmkyknQAUAv9+Pjo4ODA4OavvU1tr8aaKqGVUviiL27NlTsD0ajRbtXzpfEKrmQZKknPMmIpoPg1QiogyLxVLWfKBqq2ixYx0OR9HBRhaLRQvQygnUgsFgyamf1EBYDVKDwSAOHDhQ9NhKR/eXCigtFgtEUdTOU21tzQ+ms6npxONxBqlEVDYGqUREizQ0NARA6YdZbEBQseB1McGiuqRoqeeYTCZtmVJRFLUFAfSktiDnEwQBAwMDOHTokHaL32q1wu/3L9itgCP8iWgxGKQSES2SGrw5HI55WxCzLaYFsVj/1KUmCELJfDidTm1qqVAohMHBQTgcjgUXO+CE/kS0GJyCiohokbIHStXC7t27500/Ho9rAZ/6OxaL6ZoHi8VStOUzHA7nzHTg9/sRi8XgdDpLLkigpsNb/US0GAxSiYgqoI7SL6bakeyCIMDpdOLBBx8sun9wcDCn36vX6y25yEA8Hq/oNrvJZMLo6GjRfaWWji3V8iqKIgNUIlo0BqlERBmLCeb8fj8EQciZxxRQAkiXy7XotPP333vvvZAkqaDPqzqrQPZcqWpeBgcHc44NBoMQRbGi7gM2m61k4Kummy0UCqGnp6fo8UNDQ1rrMBFRuQyyLMvLnQkiouWiTqavDlZSp5UqZ5Q/AK1F02w25zwfUAZAHT58WAsenU4n9uzZo01qr055pS4b6nQ6cxYDyE8/FovBZrOVnMzf5/Np+RgdHUVPTw98Pp+2WMC9995bdLL9YkRRRHd3N/I/IsLhMKLRKARByAmsLRZLyf65Ho9n3nwTERXDIJWIiIqy2WxljdpfSHd3t+59Zolo9ePtfiIiKqq3t7dkv9tyDQ4Olt16S0SUjUEqEREV5XQ6IUlSVbMYBAKBsrtOEBFlY5BKREQlBQKBimcrCAaDcDgcnB+ViCrCIJWIiEqyWCzw+XxFl3qdTzQaRSwW0waJEREtFgdOERHRgqLR6KL6li72eCKifAxSiYiIiKju8HY/EREREdUdBqlEREREVHcYpBIRERFR3WGQSkRERER15/8BV4LBZkVnAZ4AAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "mode_predictions = np.empty((inputs.shape[0]-1, 2))\n", "plt.rc\n", "fig, ax = plt.subplots(figsize=(8,4))\n", "for i,motion in enumerate(inputs[1:,:]):\n", " conf.decimation = 8 # transfer function decimation\n", " A,B,C,D = mdof.system(method=\"srim\", inputs=motion, outputs=outputs[i+1,:], **conf)\n", " ss_modes = modal.system_modes((A,B,C,D),dt,decimation=conf.decimation)\n", " mode_predictions[i,:] = [[1/v[\"freq\"], v[\"damp\"]] for v in ss_modes.values()][0] # save predicted period and damping\n", " conf.decimation = 1 # transfer function decimation\n", " periods, amplitudes = transform.fourier_transfer(inputs=motion, outputs=outputs[i+1,:], step=dt, **conf)\n", " ax.plot(periods, amplitudes, color=\"blue\", label=[\"Fourier\" if i==0 else None][0])\n", " periods, amplitudes = transform.response_transfer(inputs=motion, outputs=outputs[i+1,:], step=dt, periods=periods, threads=10, **conf)\n", " ax.plot(periods, amplitudes, color=\"green\", label=[\"Response Spectrum\" if i==0 else None][0])\n", " ax.vlines([1/v[\"freq\"] for v in ss_modes.values()], 0, max(amplitudes), colors=\"r\", linestyles=\"--\", label=[\"State Space\" if i==0 else None][0])\n", "plt.legend()\n", "plt.xlabel(\"Period (s)\")\n", "plt.ylabel(\"Amplitude\")\n", "print(\" period(s) damping\\n\", mode_predictions)" ] } ], "metadata": { "kernelspec": { "display_name": "base", "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.12.0" }, "orig_nbformat": 4 }, "nbformat": 4, "nbformat_minor": 2 }