Created
October 27, 2022 20:45
-
-
Save kandersolar/672763bb0dc8c23432a2072ba011ce9f to your computer and use it in GitHub Desktop.
Showing how calculating interval energy from cumulative energy depends on how the cumulative energy was recorded
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
{ | |
"cells": [ | |
{ | |
"cell_type": "code", | |
"execution_count": 1, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"import pvlib\n", | |
"import pandas as pd\n", | |
"import matplotlib.pyplot as plt\n", | |
"\n", | |
"times = pd.date_range('2019-06-01', '2019-06-04', freq='1s', tz='Etc/GMT+5')\n", | |
"location = pvlib.location.Location(40, -80)\n", | |
"cs = location.get_clearsky(times)\n", | |
"\n", | |
"# \"real\" and \"intantaneous\" values according to the meter itself\n", | |
"# GHI looks similar enough to a power signal that we can just use that:\n", | |
"power_1s = cs.ghi # kW\n", | |
"# and the accompanying cumulative energy, again according to the meter itself\n", | |
"cumulative_energy_1s = power_1s.cumsum() / 3600 # kWh" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 2, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"# interval energy is usually a straight average:\n", | |
"interval_energy = power_1s.resample('15T').mean() / 4 # left-labeled average, kWh\n", | |
"# but the historian/DAS has a few options for reporting aggregated\n", | |
"# cumulative energy values. simulate each in turn:\n", | |
"resampler = cumulative_energy_1s.resample('15T')\n", | |
"data_cases = {\n", | |
" 'cumulative_energy_eoi': resampler.last(), # end-of-interval, kWh\n", | |
" 'cumulative_energy_aoi': resampler.mean(), # left-labeled average, kWh\n", | |
"}" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 3, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"# two simple ways of estimating the interval values corresponding to a cumulative energy signal\n", | |
"calculation_cases = {\n", | |
" 'simple diff': lambda kWh: kWh.diff(),\n", | |
" 'average diff': lambda kWh: 0.5 * (kWh.diff().shift(-1) + kWh.diff())\n", | |
"}" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 4, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAk0AAAJOCAYAAACqbjP2AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAA9hAAAPYQGoP6dpAAC4h0lEQVR4nOzdd3gUVRfA4d9JhRBaQhUMKGJFBQuCNBVpIr2DFAuCgqIIUgQBQdqnKCoqICq9l1CkKyJKUQQBC9JDLwmEkpB6vz9mE5clkA0kmWxy3ufZh+zMnZmTTXI5c9uIMQallFJKKXV9XnYHoJRSSinlCTRpUkoppZRygyZNSimllFJu0KRJKaWUUsoNmjQppZRSSrlBkyallFJKKTdo0qSUUkop5QZNmpRSSiml3KBJk1JKKaWUGzRpUqkSESMig9P5nOtEZF16nlNlHSJyUES+tTsOlbVoXaLSKqvVJZo0qQwjIveKyGARKW13LEopz6V1icoqfOwOQGVr9wKDgHXAQZd9tTM7GJWp7gIS7Q5CZRtal+RcWaou0aRJ2cIYE2t3DFmNiAQYY6LsjiM9GGNi7I5B5Qxal1xN65KMo91zmUxESojIJBE5JiIxInJARL4QET9H87NJ4ZhOjrEApZ22HRSRpSLyhIj8JiLRIrJTRJ5w7G/qeH9ZRLaKSAWXc6Y4DkBEvhWRg6l8D6VE5HMR2e24briIzHWJrxMw1/H2B0f8xim+5OuLSFERiReRQSlc6y7Hcd2dthUQkY9F5LDjM9wrIn1EJM2/zyLymIisEJFIEYkSkR9FpIpLmcGOGO5wfD7nHOW/EZGAFM75nOMzjxaRCBGZJSK3upRZJyK7RORhEVkvIlHAcMe+YBGZKiLnHdeaLCIPOmLo5CjzvON9hRSu319EEkSkRBo+B7c+UxHJIyIfOpXbLSK9RERcymWpcQjZkdYlyfFpXaJ1SabRlqZMJCK3AFuAAsAE4B+gBNAcuOoPxg13ADOA8cA0oBewRES6Yv3RfO4o1w+YIyJ3GWPSo5nzUeBxYBZwBCgNvAKsE5F7HXc464FPgNcdsfztOPZv15MZY06KyI9AS2CIy+5WQAKOStNRsfyI9bmNB8IcsYwAigNvuPtNiMhTwHJgq+O6icDzwPciUs0Ys8XlkDnAAazP8yHgJeAU0MfpnO8AQx1lvwIKA68B60WkgjHmnNP5gh3Xn4X18zvpqFiWABWBL7B+RxoBk11imQeMA9oB21z2tQPWGWOOuvk5uPWZOiqzxcCTwCRgO1AH+J/j2DfduZ66eVqXAFqXnHM6n9YlmcUYo69MemH9siYAj6SwT4DB1o/kqn2dAAOUdtp20LGtstO22o5tUUCI0/aXHdufcNq2DuuPwfVa3wIHXbYZYLDT+9wpHFfJUa6907bmrte91vWdYiznUu5PYK3T+wHARaCsS7kRQDxwq5s/CwH+BVYA4vy9AfuBVU7bBjtim+RyjgXAGaf3pRwx9HcpVw6Ic97u+P4N0MWlbFPH9h5O27yAtY7tnZy2zwCOAl5O2yq4lnPjs3DrM8WqcA3wjku5uVj/SZRx+f381q6/tez+QusSrnV9tC5J2q51SQa8tHsukziy/sbAEmPMb677jeO3I43+MsZsdHq/2fHv98aYsBS2334D17iKMSY66WsR8RWRYGAvcA7rrulGLMD6o2rldO5yWANAZzuVawH8BJwVkUJJL2AN4A1Ud/N65YGyWJVFsNN58mBVKtVTaKL/0uX9T45j8zneN8WqlOa4xHYC2IN1V+UsBvjGZVtdrEpxYtIGY93Rj0vhe5gC3OJy3nZANDA/pW/6Gtz9TJ/B+o/6E5fjP8T6j6NeGq6pbpDWJanSusSidUkG0O65zFMYyAfsSsdzOldmGGMiHd3Bh13KRTr+LZgeFxWR3FjNys9jNaU690Hnv5FzGmPOiMharGb1gY7NrbAqvwVORcsCDwCnr3GqIm5esqzjX9emamf5gbNO78Nc9iftKwicd5xTsCq1lMS5vD9qrh7EWgo4bq4exLk3hfOtBo5jVW5rHRVzGyDUGHPhGjGkxN3PtBRwLIVz/+20X2U8rUuuQ+uSZFqXZABNmrKWa90hel9je0IatztXSMblfWrXcvYpViX3MbARqyI1WP3pN9N6OQv4RkTKG2O2Y1V6a40xZ5zKeGH9gY++xjn+dfNaSXH2xupPT8lFl/epfa5eWJ9DvWuUdT1fdApl3GaMSRCRGUBnEXkVqIJ1tzgtjadKr89UZR1al2hd4jatS9ynSVPmOY11B1HuOmXOgjX7wFw5yC8jsu6zpNzE7s61mgOTjTFvJW0QkVxYg1KdpbWbYBHW4MFWjrvcO7H6wp3tAwKNMWvSeG5X+xz/nk+HczmfU4ADxpgbrRwOAU/K1VOG77hG+SnAW0ADrAr2NLAyjdd09zM9BDwtInld7hDvdtqvMp7WJalbhNYlWpdkAB3TlEkcfcmLgAYi8ojrfsdsgqQ/vupO2/MAHTMgpH3A3SJS2OlaD2LdYaQmgavvLF/j6jvLS45/C7gTkKNyX4l1V9gaiMX6zJzNASqLSB3X4x1TXd29EdiK9Rn0EpHAFM5V+OpDUrUA67MZlMK0WXGM10jNSsAX6Ox0rBfQLaXCxpgdwA6s2TfNgFnGmPg0xu3uZ/od1s+4u0uxN7H+U1uexuuqG6B1Seq0LgG0LskQ2tKUufpjzUr5UUQmYPXfFscaPFcVWIXV1z1JRP6H9UfzAlbGH5LOsXwN9ARWisgkrL7mrlgzTPJd70BgKdBeRCKBv4DKwNNAuEu57VjfQx8RyY81WPF7Y8yp65x7NlaT8KvASpe7ZLCmpDYEloq1dsdWrAGX92PdtZYGzpAKY0yiiLyE9cf5p4h8gzV7pATWYMjzWHdcbjPG7BORAVh3tKVFZBFwAbgNaII1NfyDVE6zCGsq+YcicgfWNOGGQFDSZVI4ZorTedPanA7uf6ZLgB+A98VaR+cPrN/nRsDHxph9ridWGUbrEq1LtC6xg93T93LaC6vCmoy1JsdlrDuUzwA/x/6HgE1YlcIhrMy7EylPE16awvkN8JnLttKO7b1ctrdzXD8Ga32O2rg3TbgAVkV5GusPeQXWUvcHcZkainXXsg9rEGbylGGuPU05L9Y0ZwO0u8ZnGIi1XsseR+yngZ+xmpZ90/jzKI81O+SM4+dxEKuyfcqpzGBHPIVcjr3q5+LY3hRrBslFx+tvx8/4Tqcy64Bd14ipEDAdq7I9hzUr5nHHtVqlUL6Y4/PdfRO/l259po5yY7D+U4jFGqPQC6ep1k6/n9/eaDz6cutnpnWJ0boErUsy9SWOoJRSWZiINAYWAlWNMT+77CuENfPlPWPMUBvCU0p5CK1Lbo6OaVIqi3FMw3Z+7401zuM88HsKh3TCGh8wNcODU0p5DK1L0p+OaVLZjogEAX7XKZJgjLnWOiJZwaeOym4j4I/VRP841irAzosBPoW1YN87wCJjzEHnkzjOkdpaNxFGH3iqVIq0Lkner3WJg3bPqWxHrId31rhOkUPGmNKZE03aiUhbrP7/O4BcWIvRfWGM+cyl3DqsCvBn4Dnj8nwosR7I6bpKsKsnjTHr0iVwpbIZrUuS93dC6xJAkyaVDYnIw1x/xeJo17787EhEigP3pVJsqzHmbCpllMqRtC6xaF3yH02alFJKKaXcoAPBlVJKKaXcoAPBSV5B9xasdUKUUp4nL9YDQG1rOtd6RCmPl2o9okmT5RbgiN1BKKVuSkmshfLsovWIUp7vuvWIJk2WCwCHDx8mX77UVv1XStltzZo1tGvXjsuXL/Pwww+zdetWsL+FR+sRpTzIxYsXadu2LT/++CO+vr7ExcVBKvWIJk1O8uXLp5WdUlnc3LlzadeuHXFxcdStW5dvvvmG4sWL2x1WMq1HlMr6IiIiaNq0KZs3byYwMJCZM2fSoEHqjwjUgeBKKY/x1Vdf0bp1a+Li4mjVqhWhoaEEBATYHZZSyoMcO3aM6tWrs3nzZoKCgli7di3Vq1d361hNmpRSHuGDDz6gc+fOJCYm8vLLLzN9+nT8/K63WLNSSl1p//79VKtWjT///JNbbrmF9evXU7FiRbeP16RJKZWlGWPo378/vXv3BqBPnz58+eWXeHt72xyZUsqT7Nq1i6pVq7J//37KlCnDhg0buO++1NbsvJKOaVJKZVmJiYl0796dL774AoCRI0fSp08fm6NSSnmazZs3U69ePc6ePcv999/PypUrb2gspCZNSqksKS4ujo4dOzJz5kxEhC+//JKXX37Z7rCUUh5m7dq1NGrUiEuXLlGpUiWWLVtGUFDQDZ1LkyalVJYTFRVFy5YtWbZsGT4+PkybNo1WrVrZHZZSysMsXLiQ1q1bExsbS61atVi4cCF58uS54fPpmCalVJbSe/ovFLq7IsuWLcMvVy4WL16sCZNSKk1Ctx3l9qa9adqsObGxsTRr1owlS5bcVMIE2tKklMpCag1fzPpPehJ7ch/iF0DBZoOILfaA3WEppTxE6LajjF27h+0rZnB27UQA8txfi7vavIO/v/9Nn19bmpRSWULjUaGs+7AbsSf34RWQn2JtR5Cr5H1M3xJmd2hKKQ8Quu0or8/axu8LJyQnTHkfbUxwvdcJ3XEyXa6hLU1KKdu9P2MtS9/vTMKF03jnLUzR1sPwDSoBQLuKITZHp5TK6kK3HaX/wj84u3YCF7YuAaBAtfbkq9wSEaFJhRLpch1NmpRStvpo5koGdW5FQlQkPkElKdpqKD75CpPbVxjZ9EEapVNlp5TKnkK3HeX1mVsJ/+5jLv35AwBBtbqS96FnAahSJpihje9Pl2tp0qSUsk3H4ZOZOuRVTGwUfkXLUKTFELzzFKBKmWCmd65kd3hKKQ8w5ec9nF44nOi9m0G8CK7/JoH3PQlA+ZL507Uu0aRJKZXpQrcdZdC4aeyYPAgTH4N/yfso0vxdcucJ5O06d/NitdvtDlEplYWFbjvK9C1hNL0viD2T3yF672bEx49CjfoQcMdjyeX61rsnXa+rSZNSKlOFbjvKS4M/4czSDyExgdy3P0Khxn3x8s3F6GbaHaeUur7QbUfpMXs7CdHnCR3Skdjje/D2D6BQ04HkCrG64QoF+vFZm4eoVCY4Xa+tSZNSKlP1eX8MZxZ/BBgC7qlBofpvIt4+tK8UogmTUipV07eEEX/hDKdmv0tceBjeufPR8b3xrD2TN7nMwPr3pnvCBJo0KaUyUcce/dk9fwwAgeXrEVSrK2WK5KNHzbKaMCml3PJU8UQW9OlDQuRJvAODKdpqGGvP5KV9pRB2n7xIu4oZdwOmSZNSKsMt+v0Iffv3Y/fKaQDkq9SCAtU7UKZwIN/3esLe4JRSWV7SGKYqBS4y8o0OJESexC/oFgq3HIpP/qIA7D55kTldKmdoHJo0KaUy1ILfwujYuQsXt68AoMATz5P/sWYA9KhZ1s7QlFIeIGkMU8zRv5k/dzCJMZd44IEHePODbxm89lhyucxY002TJqVUhomNjeWVlzpx8Y8fACGobnfK12xKcF7/DG1CV0plH9O3hBF9YBunFw7DxMWQu+S9vP3pTNpVv5eChYowfUtYptUnmjQppTJEVFQU1Wo/y6k/fgAvHwo16EWeu6vyuo5fUkqlwW0XdjJ33hBIjCfXbQ9RqHF/3vnuAIF589OoQolMrU80aVJKpbvp6//ijRfbcmbvH4iPP4Wb9Cf37Q9ze6E8mjAppa4rafxSu4ohhG9byQd9XoHERPLdW40Cz/REvH0BqwUqs+sTTZqUUulq8to/6Ny2KXGn9iP+eSjSfBC5St4L6BgmpdT1JY1fAlgzexJnf5gEwEsvvUS9LgPpOW9nclk7nkupSZNSKt189d0WunVoSlz4UbwCClC01Xvcc+/9OoZJKeWW6VvCMMYQ+dM0IjfOBqB3796MGjUKEcHb2ztTxzC50qRJKZUuxi34kR6dmpNw4Qze+YpQtNVQfINK6BgmpZTb2jxSkhUT3ufC78sAuL3eSzze+nVEBCDTxzC50qRJKXXTfv/9d3p2bEzCxXP4BJWkaKthFChclPcb368Jk1LKLXFxccz/qJ+VMIkQVOsVEh54hjfm/IGIZIm6xMvuAJRSnmvTvnCaDPyKx6vVIPbiOfyK3UGxdqPwyVdIEyalVKpCtx2l5fiNzNm4l2bNmjF9+nR8fHx47IVB5K3wTHK56VvCbIzyP9rSpJS6IZv2hdOw76ecWTQCEx9L7pD7GT1+KusPRen4JaVUqpIGfSfGRLH4/ZeJObyLXLlyMW/ePOJvKZ88IBzsGfSdEk2alFJpFrrtKG+OGMfpBaMhMYHcZR6lUKO+FCsczJy6D9odnlLKA0zfEkZCVCSn5g4i9sRefHLlYeXK76hevfoVZbLSTZgmTUqpNAnddpQX+g4nYtUXgCHPvU8Q/Mwb5PL3o2LpILvDU0p5iDqlfFjYvw/xEUfwyp2PCi//j7N5yyTvt3vQd0o0aVJKuc0Yw4D3hhGx6ksA8j5Un9sbdKfNY6VpUzGEkOAAmyNUSnmCPXv2MLRrC+IjjuBfoAjBzYdwKleJ5C65rJYsJdGB4Eoptxhj6NOnD7sWWQlT/sqtKPh0V4Y0up8+9e7WhEkplaKkwd6h244C8Mcff1CtWjUOHTrEnXfeydNvf4lv8K3J5bPKoO+U2Jo0iUg/EflVRC6IyCkRWSQid7mUySUi40QkXEQuish8ESnqUiZERJaJSJTjPP8TEW1FUyqd/PzvKe57qin/+9//AHi+57vU7tiDT1pXyLJ3hEop+yUN9t5yIIIes7cz8ptF1KhRg5MnT1K+fHl++uknOteteMUxWWXQd0rsTixqAOOAXx2xDAdWici9xphLjjIfAfWBFkAk8BmwAKgCICLewDLgBPA4UByYAsQB/TPtO1Eqm/rp7+PUbdKSqN0bQLwo+szrDO7fR1uWlFKpcm41it6/lQFjhpMQF0PVqlVZsmQJBQoUoFGR/8pmpUHfKbE1aTLG1HV+LyKdgFPAw8B6EckPvAi0NcZ87yjzPPC3iFQyxmwCagP3Ak8bY04C20VkIDBKRAYbY2Iz7ztSKnuZ/cseXu7Ylqi9v4GXD4Ua9ibXXVXYcjBCkyalVKraVQxhy4EILv2zgTNLPoDEeOrWrcv8+fMJCPivDsmKg75TktXGNOV3/Bvh+PdhwBdYk1TAGPMPEAZUdmyqDOx0JExJVgL5gPtSuoiI+ItIvqQXkDf9vgWlsofp6/+iY8tGnN/7G+LrT5Hmg8hzVxX8fbx0lhxajyjljkYVSvCs31+cWTwaEuOpWrsheZ7tx+rdZ+0O7YbY3T2XTES8gI+Bn40xuxybiwGxxphzLsVPOvYllTmZwn6cyrjqBwy6mXiVys5OnjxJt3aNiTm6By//PBRpMZhiZR+g9aMhOkvuP1qPKJWKDz74gHFD3wagdtN2/FOmJYcPX+C3LD5L7lqyUkvTOKAc0DoTrjUCq1Ur6VUyE66plEc4dOgQVatWJfLIHrzyFKBo25H4l7iHwQ3u01lyV9J6RKlrMMbQv39/evfuDUCfPn3IV+tVxMs7uUxWniV3LVkiaRKRz4BngSeNMUecdp0A/ESkgMshRR37ksoUTWE/TmWuYIyJMcacT3oBF24mfqWyi8/mr+Oe8o+yd+9eSpcuzbjpS6j62MOMbVXe4+4IM5rWI0r9x3lZgcTERLp168aIESMAGDlyJCNHjuS5x0pdcUxWniV3LbZ2z4mIAJ8CTYAnjDEHXIpsxZoFVxOY7zjmLiAE2OgosxF4R0SKGGNOObbVAs4Df2Xsd6BU9vHh9O94u3MbEqPP4xscwsAv5/JCnUfoandgSqksLWlZAYDNe0/xYf9v+GlFKCLCl19+ycsvvwz81xXnCbPkrsXuMU3jgLZAI+CCiCSNQYo0xkQbYyJFZBIwRkQisBKhT4GNjplzAKuwkqOpIvI21jimYcA4Y0xMZn4zSnmqdevW0felFiRejsKveFmKtBjCioNxvGB3YEqpLC+pmy0x7jJnQkcRtu9XfHx8mDZtGq1atbqirKfMkrsWu7vnXsEaC7AOOO70cv6U3wSWYrU0rcfqcmuatNMYk4DVtZeA1eo0DWudpnczPHqlsoF3Pv6WmrXqEH85Cv+QByja6n28c+fzyKZzpVTma1cxhMSYS5yaM4jofb/i55+Ltz/8mspPN7A7tHQnxhi7Y7CdY7pwZGRkJPny5bM7HKUyTc/3P+WjgW+ASSR32Ur0GPYp+87GeVTT+fnz58mfPz9AfsfYIltoPaJyqtOnT1Opek32/7OT3HkCKdRsEF7F78Hfx4vVb9bwiMkj7tYjdrc0KaVsMm7cOD4a8DqYRPLc9ySFG/dj39k45nSp7DEJk1LKXocPH6ZatWrs/2cnhQsX5qWRk/Eqfg8AMfGJzPTAGXLXo0mTUjmMMYZhw4bRvXt3API+3IDg+m8iXt7aJaeUuibXB+/++++/VKlShd27d3PrrbeyYcMGEgqWvuKYsIgoGyLNOHYPBFdKZaJFvx+hT5/e/LtmFgDvvvsuFRp1Zsavhz2qS04plbmcZ8htORDB/n92MaJHe06fPs1dd93F6tWrufXWW2nvHc6ynceTj2tfqdQ1zuiZNGlSKodY8OshOrzYmUs7VwPwwluDGTLEWtC68UO6LqNS6tqcF6K8fORP3h77HvGXL/HQQw+xYsUKChcuDEClMsHM6lyJhduP0qR8CSqVCbYr5AyhSZNSOUBMTAyvd+nEpZ3rQLwIrvc6F8rWtjsspZSHSHrwbvS+3zi9aAQmPobq1auzePHipAHUySqVCc52yVISHdOkVDY36+d/KVWhGke3rQNvHwo37kvg/U/r+CWllNsaVShBiwJhnF4wFBMfQ/369VmxYsVVCVN2py1NSmVj03/8kxfaNiP22G7ENxdt3vmEuGLldPySUipNJkyYwIf9u2GMoU2bNjR7Yzgdp2zPcXWJJk1KZVPHjx+nW7tGxB7bh1euQIq0GEJcsXLM6VLZ7tCUUh5k1KhR9O3bF4BXXnmFWi/24825OwBrUDiQYxIn7Z5TKhs6cOAA1apVI/LoPrwDgyjadiT+t9ylXXJKKbcZY+jbt29ywtSvXz/GjRvHzN+OXFFuejZbi+l6NGlSKhsJC4/io7nfU/nxKuzbt4/bbruNcTOWULXiQ4xtVT7H3A0qpdLGdQ2mhIQEunbtyqhRowAYPXo0w4cPR0SuuvnKSTdj2j2nVDYRFh5Ftd4TOTJzIImXL3Dn3ffww9o13HLLLXSxOzilVJblugZTXFws88b0Y/bs2YgIEyZM4KWXXkoun3TzNX1LmI5pUkp5pmET53J4ej9MbDR+xe+k1eBJ3HLLLXaHpZTK4py71xLjLvPaC+048edGfH19mT59Oi1atLjqmEYVSuSoZCmJds8plQ2Ehoby9bsvY2KjyVXqQYq2fp/weH+7w1JKeYCk7rXEyxc5NeddTvy5kdy5c7NkyZIUE6acTJMmpTzcG0PH0qRpMxLiYsldthJFmg/Cyy93tnt8gVIqYzSqUIIhtUpyYcFAYo78Rf78+Zk6bwkXCt1HWHj2enbczdLuOaU8WOc+Q/lq9LsA5ClXk4GjPiYiOjFbPr5AKZUxwsLCGP5qS84d3kORIkWYPCeUHqvPEhP/B/4+Xqx+swYhwQF2h5klaEuTUh7IGMN7772XnDDlfbghwc/0YOvhC4xq9oAmTEopt/zzzz9UqVKFf//9l1KlSrFhwwb+iCpATHwiADHxiczMQUsKpEaTJqU8TGJiIj179mTQIOthu/mrtqNgzc6IeOWoqb9KqZvz+++/U61aNY4cOcLdd9/Nhg0bKFu2LGERV3bJub7PybR7TikPsv/keV7s/BLrlswFYOzYsZSq1ixHTv1VSrln075wpm46xK1BuWlbsRQhwQGsX7+eBg0acP78eR5++GFWrFhBoUKFAGhfqRTLdh5PPl7HR/5HjDF2x2A7EckHREZGRpIvXz67w1EqRf8ejeDhmg24uPsXEC/GfPYlb77a2e6wbHf+/Pmkh4bmN8actysOrUdUVrRpXzitJ25Kfu/nLfQrd5lXX2jH5cuXqVGjBosXL77qd3bTvnAWbj+aY8ZHuluPaEuTUh7g4sWL1H+2Phd3bwJvHwo36kvsbdXsDksplcUt3H70ivdnd66j86iPSEiIp0GDBsyePZvcuXNfdVylMsE5IllKKx3TpFQWFxERwdNPP83e7ZsQv9wUaTGEgLKVdJyBUipVTcr/12V/Ydt3nFnyAQkJ8bRr14758+enmDCpa9OkSaks7JtVW7nt/kfYvHkz+QoUpGirYeQu9SCg4wyUUqmrVCaYmS89RrEDy4lY9Tlg6NatG1OmTMHX19fu8DyOds8plUWNX/IL3Ts0I/7cCbwDgxg2cS6PVngwR40zUErdHGMMCyeMZvOccQAMGDCAh5t0ofXEzTp55AZo0qRUFrRr1y7eaN+I+Mgz+BQoRpFWw/gpPDev6TgDpZSbEhIS6Nq1K1999RUAH374IWWebHXFw3kBTZzSQLvnlMpiNm/eTPXq1bkceQbfwqUp2m40vgWK6RpMSim3xcTE0KZNG7766iu8vLyYNGkSPXv2vOLhvMBV79X1adKkVBaydu1aatasydmzZ6lUqRJfz1nC4/ffwdhW5fVuMAOIyOsiksvxdYiIiN0xKZUWoduO0nL8RkK3/TdL7tKlSzRs2JC5c+fi6+vLnDlzeOGFFwCuuvnSm7G00e45pbKAsPAo3vl4EjNGvkVifBy1atViwYIFBAYG8twTdkeXrY0BZgGXgQNAceCUrREp5abQbUev6mqrXjqAZ599ll9++YWAgAAWLVpErVq1ko9JuvnSBXFvjCZNStksLDyKip0GcHLZWDCJBNz5OAM+/pbAwEC7Q8sJjgHNROQ7QICSSS1Prowx2o+hshTXrrVJa7bz7rT+7NixgwIFCvDdd99RuXLlq45rVKGEJks3SLvnlLLZoBH/4+TSj8Akkuf+pynUqA/L/j5jd1g5xTDgY2A/YIBfsVqcnF8HHf8qlaU4d63FR57kt0+7s2PHDooWLcqPP/6YYsKkbo62NCllE2MMQ4YM4dsPhwCQ99HGFHzyRUTkigXpVMYxxkwQkZlAKWAH8DQQbm9USrknqbXoi9D1/PbNO4SfOkHp0qVZvXo1d9xxh83RZU+aNCllg8TERBq2e4lls74BoFf/QQRVacmRs5dpX6mULiuQSUTkdWCCMWaXiDwPbDTGRNsdl1LuKpFwnN8+e43w8HDuvfdevpkdyvbIXPiFRxESHGB3eNmOW0mTiOSz80GYSmUn8fHx1G7Shh+WzgMgqFZXqjbvrGMM7OE8EPxrYDmgSZPyCOvWraNhw4ZcuHCBRx99lInT59Ny8i5i4hPx9/Fi9Zs1NHFKZ+6OaTorIkUAROR7ESmQcSEplX1dvnyZ5s2bWwmTeBH87FvkfehZXSvFPkkDwUvx30DwkJReNsep1BWWLFlC3bp1uXDhAk8++SRr165lxd5LxMQnAhATn8hMrVfSnbtJ00Ugqb/gCSBdHlgjItVFZImIHBMRIyKNXfZ/69ju/FrhUiZIRKaLyHkROScik0REpx2pLOfChQs888wzhIaG4uvnT+Gm7xB435OArpViIx0IrjzO9OnTadKkCTExMTRq1IjvvvuOvHnzXvUQb32od/pzd0zTGuAHEfnb8X6hiMSmVNAY81Qarp8H+AOrWXzBNcqsAJ53eh/jsn861toqtbCSuW+ACUDbNMShVIYKDw+nXr16/Prrr+TNm5fFixcTmb+srpViMx0IrrKq0G1HU6wfxo0bR/fu3QFo3749X3/9NT4+1n/l7SuVYtnO48ll9aHe6c/dpOk5oCNQBqgB/AncdAprjFmONYaA6yzEG2OMOZHSDhG5B6gLPGqM+c2x7TXgOxHpZYw5drMxKnWzvl75G28+35Lzxw8QHBzMihUreOSRRwB95lNWYIy5ACQNBP/ZGON6Y6ZUpkpp0cqG5W/h/fffZ+DAgQC89tprfPzxx3h5/ddhVKlMMLM6V9KHemcgd5MmX2PMlwAi8gjQxxhzLsOiutITInIKOAt8DwwwxiTdCVYGziUlTA5rgETgMWBhSicUEX/A32lT3nSPWingy8U/071DMxIiT+IdGMzg8fOSEyaVtRhjJqelvNYjKqO4jnGctvkQ66eNYcyYMQAMGjSIQYMGpdjYUEkf6p2h0jwQHKvfP7OsADoANYE+WK1cy0XE27G/GC6PPDDGxAMRjn3X0g+IdHodSd+wlYIdO3bwZvtGJESexKdgcYo9N5r1Z/xTP1BlGhGJEJFCjq/POt6n+ErhcK1HVIZwHuNoEhM4u+KT5ITp448/ZvDgwdfrnVEZyN2WpqSB4KewEpd0GQieGmPMLKe3O0VkB7APazD62ps49QisqcZJ8qIVnkpHGzdu5JlnnuHy+XP4Fi5N0ZZD8Q4sqAO+s543gQtOX6flplDrEZUhkrrtp/y8l2MLR7H2++V4eXkxadIkOnXqZG9wOdyNDAQX0m8geJoYY/aLyBngDqyk6QRQxLmMiPgAQY591zpPDE4DyjVjV+lpyBczGfrGCyTEXubxxx/n1eETCP3nvA74zoKcu+SMMd+m8VitR1SGqVk2P+P6vM+m71fj5+fHrFmzaNKkid1h5Xi2DgRPKxEpidXilTQ9YCNQQEQeNsZsdWx7CqvbcXNmx6dUnw8mMLpPN0iMJ1fpCrz+v29p9XhZ2tWwOzKVEhHJ525ZXeBXZZazZ8/yzDPPsGnTJvLkyUNoaCh5Spen2/TfCQkKoE3FEF200iZuJU2Oxwqk+0Bwx3pKzg/IuU1EymONSYoABgHzsVqNygCjgb3ASkdcfzvWbZooIl2xug0/A2bpzDmV2b7++mtGv/0KmEQC7qpKoWffYv7OM7R6vKzdoalrO4f7XXLeqRdR6uYcP36cOnXqsHPnTgoWLMjy5csxhe6g9cRNyWW+/vmArvZtkzQ/e84Y82Q6Xv8R4Aen90njAyYDrwAPYLVwFcBauXcVMNBlSnA7rERpLdasufnA6+kYo1KpGjNmDG+99RYAgQ/UJqhON8TLW8cwZX3O9VlpYCTwLVYrNlgzdDtiDfpWKkMdOHCAWrVqsW/fPooXL86qVasoV64cfebvuKJcTHwiWw5GaNJkA3efPTcGK1m55Pj6mowxPd29uDFmHdYYqWup48Y5ItCFLJVNjDG8++67DBs2DIBevXpRpU0PZvx6WMcweQBjzI9JX4vIu0BPY8xMpyKLRWQn8DLWzZxSGeLPP/+kdu3aHDt2jNtuu401a9Zw++23A9CkfAlm/3o4uay/jxcVSwfZFWqO5m5LUwX+mzFX4TrlMnM5AqVs9cue03R//TW2rZgNwPDhw+nbty8iQuOHStocnboBlYGuKWz/Dfgqk2NR2VRKK33/+uuv1K1bl4iICO677z5WrVrFLbfcknxM0qKVUzcd0jFNNnN3TNOTKX2tVE614Z8T1G3elkt//gAIvYeMol+/3naHpW7OYaAz8LbL9pcc+5S6KSmt9J3v3L80bNiQixcvUrFiRb777juCg69enFIXrcwa0jymSamcLjo6mvZtW3Hpz/Xg5U2h+j25VKam3WGpm/cmMF9E6vHf7NuKQFmgmW1RqWzDdaXv0ROmsfWbQcTExFCzZk0WLVpEYKA+bz4rc3dM07UepnsVY0zTGw9Hqazt/PnzNGzYkIPb1iM+fhRq3I+AMo9ya1Buu0NTN8kY852IlMWahHKPY/MS4EtjjLY0qZvWrmJIcgvTxV3fs2nFWBITEmjcuDEzZ84kV65cNkeoUuNuS1Ok09cCNHFsS3rm28NYM9zcTq6U8jRTv99JjxdacvbQPwQG5qVA4wF4l7gPP2+hbUV9mnh2YIw5Arxjdxwqe0oawzR41IccWvYRAB07duSrr77Cx0c7fjyBu2Oank/6WkRGAXOArsaYBMc2b+BzQBd/U9nSpBW/8spzzYgLP4xX7nwMGT+b5nWqs+VgBBVLB+mgTKVUqowx/LFkEttnWwlTjx49eOOdYUz6+RAA9coV17oki7uR1PYFoGpSwgRgjElwLEXwC6CjYVW2smfPHl5v15C4iBN45y1E0VZD2XQhPz2DA7SCU0q5JTExkbfeeouPP/4YgCFDhtDx1beo9dGPxCZYE88/Wv0vq3TRyizN6waO8QHuTmH73Td4PqWyrD/++INq1aoRFXECn4K3UKzdaHyDb9VFK5VSbouPj+fFF19MTpjGjh3Lu+++y4o/TyQnTACXHYtWqqzrRlqavgEmiUgZYItj22NAX8c+pbKFX375hfr163Pu3DnKly/PGx98y7K9UbpopVLKbZcvX6Zt27YsXLgQb29vvv76azp06JBiWV9v0UUrs7gbSZp6YT0L7i2guGPbceB/wIfpFJdStlq5ciVNmjQhOjqaqlWrsmTJEgoUKEBHXVkgWxMRH+AJrGddzjDGXBCRW4DzxpiLtgansqxN+8JTXHjy4sWLNG7cmLVr1+Lv78/s2bNp1KhR8nH1yhVnzOp/iYlPxMdLmPrCY9o1l8WJMTe+iHfSE8I9/enfju8jMjIyknz53H7oucqGNu0LZ9DYr1j7xUAS4uOoW7cu8+fPJyBAK7Ks6vz58+TPnx8g/83URSJSClgBhAD+wJ3GmP0iMhbwN8aktFq48/Faj+RAm/aFX/EwXX8fL1a/WYNAucwzzzzD5s2bCQwMZPHixTz55NVrQ4eFR+mEkizA3XrkpuY4enqypJSzTfvCqd99CBErx4FJJPCeanz69UxNmHKOsVjLqDwIhDttXwhMtCUileUt3H70ivcx8Yms+PUvPunViT///JOgoCCWL19OxYoVUzw+RCeUeBQduK2UQ78hw4lY8amVMD1Yl6Bne7H9mPbI5CDVgGHGmFiX7QcBHcSmUtSk/JW/GnLhJO+/0pI///yTW265hfXr118zYVKeR5MmleMZY+jfvz/rpo4BIN9jzQmq041cfr46KDNn8QK8U9heEriQybEoD5H0MN369xenSakEoua/Q9jBA5QpU4YNGzZw33332R2iSke6BKnK0RZuPUzvnj3Yt34hAK/2fhcebKRPEs+ZVgFvAC873hsRCQSGAN/ZFZTK+iqVCUbO7KVevU6cPXuW+++/n5UrV1K8ePHUD1Ye5aaSJhHJZYy5nF7BKJWZ5m85SIdOnYj6+0dAeKX/CMa938fusJR93gJWishfQC5gBtbDes8AbewMTGVta9eupVGjRly6dIlKlSqxbNkyfjoUzfTFG3WJkmwmzUmTiHhhPZupK1BURJJmmAwFDhpjJqV3kEqlt6ioKF57qR1Rf/8CXt4UevYtzoRUtzssZSNjzBEReRBoDTwABAKTgOnGmGhbg1NZ1sKFC2ndujWxsbHUqlWLBQsWsHZPJD1mbwdIfkCvJk7Zw42MaRoAdALeBpwHTO4CXkqHmJTKUJGRkdStW5fjO39BfPwp0nQgee6prqt8K4wx8caYacaYt40xrxpjvtKESV3Lt99+S/PmzYmNjaVp06YsWbKEwMBApm8Ju6Kc63vluW6ke64D8LIxZq2IfOm0/Q9SfryKUlnG6dOnqVu3Lr///jv58uWjz5hv2B5fXJvQFSLS8Bq7DHAZ2GuMOZCJIaksbOzYsbzxxhsAPP/880yYMAEfH+u/1HYVQ5JbmJLeq+zhRpKmEsDeFLZ7Ab43F45SGefw4cPUqlWL3bt3U7hwYVauXEmFChXsDktlHYuwEiRx2Z60zYjIBqCxMeZsJsemMllYeBQzthzicEQ07SuVolKZYMCabTt48GDee+89AHr27MkHH3yAyH+/Nkk3YNO3hOkNWTZzI0nTX1jrmRxy2d4c2HbTESmVAeas2UznNk04f+Y4t956K2vWrOHOO++0OyyVtdQC3scas5n0XM2KwFBgGBAJjAc+AF60I0CVOcLCo3h6zLrkh+ku23mcWZ0rUfG2grzxxht8+umnAAwbNoz+/ftfkTAlaVShhCZL2dCNJE3vAZNFpARW61JTEbkLq9vu2fQMTqn08NGsVfR6sRWJUefwCSrJp9OXasKkUjIWa+jBL07b1orIZWCCMeY+EXkD+NqW6FSm2XIwIjlhSjJ/6yE+H/ImU6dOBeCzzz6jW7dudoSnbJTmgeDGmFCgAfA0cAkriboHaGCMWZ2+4Sl1c+YvW8NbnZqQGHUOv6JlKNZ2JL+csjsqlUWVAVJ6NNR54HbH13uAQpkWkbJFxdJB+Hn/13pk4mPZMnEAU6dOxdvbm2nTpmnClEPd0DpNxpifsJqylcqyli9fTttmTTExl/EveR9Fmr+Ll3+eqx57oJTDVuB/ItLBGHMaQEQKA6OBXx1lygKHbYpPZZKQ4ADW9HyCGVsOse/oGbZN6s/6TRvw9/dn7ty5NGjQwO4QlU10RXCV7WzaF87wcZP47tN3SIiPJ88djxLUsA++frkY07J88oBOpVy8CIQCR0QkKTG6FdgPNHK8D8Qa36SyuZDgADpXLEK9eh3Z+uuv5M2bl8WLF/PEE0/YHZqykVtJk4icxZpBkipjjD6sS9lm075w6ncbTMTKcYChVoOmfDHha7YdvUDF0kH6WBR1TcaY3SJyL1AbSBr0thtYbYxJdJRZZFN4KpMdPXqU2rVr89dffxEcHMyKFSuIL3gbfebvoEn5EnrzlUO529L0RkYGoVR6efOdIUSs/AyAwPL1KNOiL2WK5adMsfw2R6Y8gSM5WuF4qRxq3759PP300xw8eJASJUqwevVqIv2K0HriJgBm/3qYWZ0raeKUA7mVNBljJmd0IErdDGMM/fr1Y9NsaypwvkotKFC9A6UKBdocmfIkIpIHqAGEAH7O+4wxn9gSlMpUO3bsoE6dOpw4cYI77riD1atXU7p0abpN//2KclM3HdKkKQe6oTFNIuINNMGaNQfW2k2hxpj49ApMKXclJCTQrVs3xo8fD0Chp14gz6NN8fMW2lYsZXN0ylOISAXgOyAAyANEYM2UiwJOAZo0ZXMbN27kmWee4dy5czzwwAOsXLmSYsWKAXBrUO4ryrq+VznDjTyw9z5gMVAMq78foA9wWkQaGGN2pWN8Sl1XbGwsHTp0YPbs2YgI48ePp07Tdmw5GKFjmFRafQQswXoYeSRQCYgDpmGt4aSysdWrV9O4cWOioqJ4/PHHWbZsGQUKFEje37ZiKb7ecIDYBKM3ZDnYjbQ0fQX8CTyS9CgBESkIfAtMAB5Pt+iUuo7dh8/QpnVLtv3yA76+vkyfPp0WLVoAaLKkbkR5oIsxJlFEEgB/Y8x+EXkbmAwssDU6ddNCtx1N8dEm8+fPp02bNsTFxVGnTh3mz59Pnjx5rjg2aRkCvSHL2W4kaSqPU8IEYIw5KyLv8N9aJkplqNXb9tOwYUMuH/kT8fHnq2mzaNGisd1hKc8WByQ6vj6FNa7pb6xWp1vtCkqlj9BtR+kxeztA8sN0G1Uowddff03nzp1JTEykRYsWTJs2DT8/vxTPERIcoMlSDncjSdO/QFGs1iZnRUj5Qb5Kpavf/znIs/VqEXtyP+KfhyLNBxFd9H67w1KebxvwKNaq3z8C74lIIaA9oMMOPNz0LWFXvd/3w2zeeustAF588UXGjx+Pt7e3HeEpD5Hmx6gA/YBPRKS5iJR0vJoDHwN9RCRf0iu1E4lIdRFZIiLHRMSISGOX/SIi74nIcRGJFpE1IlLWpUyQiEwXkfMick5EJomITpnKpsLCwni2Tk1iT+7HK6AAxdqOIFfJe+0OS2UP/YHjjq/fAc4CXwCFgZftCkqlj3YVQ5K/NsbAb/8lTL169WLixImaMKlU3UhL01LHv3P4b8HLpIf0LHF6b4DUfgPzAH9gPQAzpfECbwOvAx2BA1hPG18pIvcaYy47ykwHimM91sUX+AZrbFVb978l5Qn++ecfatWqxfEjR/DJV5girYbhG1QCP2+hXrnidoenPJhYj6k/haNFyRhzCqhra1AqXSWNYZq26SBR6yYyd461ks7w4cPp27cv1q+AUtd3I0nTk+l1cWPMcmA5cNUvrKMSewMY5nhIMCLSATgJNAZmicg9WBXbo8aY3xxlXgO+E5Fexphj6RWrsteYGSvo16UNsRfPcffdd/PtnMVsj7AaSuuVK67jDNTNEqzhBfdhdc+pbOiZckWY+2Ef5s2Zjogwbtw4XnnlFbvDUh4kzUmTMebHjAgkBbdhLWuwxunakSKyGagMzHL8ey4pYXJYgzWY8zFgYUonFhF/wN9pU970DV2lp+FfzWdAtw6Y2Cj8it1Bv3Gzeez+sjxmd2Aq23DMmNsDBONm0qT1iGeJjo6mVatWLFmyBB8fH6ZMmUKbNm2uOaNOqZTc6OKWuYAHsAZ/XzEuyhizOB3iAithAqtlydlJp33FsJrUna8fLyIRTmVS0g8YlB5Bqoz13XffMfDVtpi4WPxvLUeRZu+ydM8lOjxld2QqG+oL/E9EXnFzvTmtRzzE+fPnadiwIT/++CO5cuVi3rx51K9f/5oz6pS6lhtZ3LIuMAVrpVxX7oxjygpGAGOc3ucFjtgUi7qGmTNn0qFDBxLj48ld5lEKNeqLl6//FQM6lUpHU7BWA/9DRGKBaOedKTyMXOsRD3DmzBnq1q3L1q1byZs3L0uXLqV69epAyjPqNGlS13MjLU2fAnOB94wxrq1A6emE49+i/DejJen9dqcyRZwPEhEfIMjp+KsYY2KAGKdjbj5ala6+/PJLXn31VYwxtGvXjiavD2P2tuPahK4y0htpKaz1SNZ35MgRatWqxT///EOhQoVYuXIlDz30UPL+dhVDkluYkt4rdT03kjQVBcZkcMIE1my5E0BNHEmSYxmDx7CmAQNsBAqIyMPGmK2ObU9hdRluzuD4VAYwxjBy5Ej69+8PQLdu3fjkk0/w8vKiWcXS9gansjV9MLlnutaYpD179lCrVi0OHTpEyZIlWb16NXffffcVxyaV1zFNyl03kjTNA54A9t3sxR3rKd3htOk2ESkPRBhjwkTkY2CAY4Bm0pIDx4BFAMaYv0VkBTBRRLpiLTnwGTBLZ855no17z/DqG2+xfdkUAAYMGMB7772nd/Aq04hIGeB5oAzQwxhzSkTqAWHGGNcFfZXNrjUm6Y8//qBOnTqcPHmSsmXLsnr1akqVSvlZcY0qlNBkSbntRpKm7sBcEakG7MR69EAyY0xangT+CPCD0/uk8QGTgU7AaKy1nCYABYANQF2nNZoA2mElSmuxZs3Nx1rbSXmQn/89Rd0W7bm4YxUAhZ9+ic5v9NOESWUaEamBtQTKz0B1rAUuTwEPAi8Cze2LTqUkpTFJhaMPUb9+fc6dO0f58uVZuXIlRYoUucYZlEqbG0ma2gC1gctYLU7GaZ8B3E6ajDHr+G9hzJT2G+Bdx+taZSLQhSw9WmxsLC92ciRM4kVw3e4EPFCbLQcjdP0llZlGAgOMMWNE5ILT9u+xbhZVFuM6June+P08/XRnoqOjqVq1KkuWLKFAgQL2BaiynRtJmt7HmmY70hiTmFphpa7n0qVLNGvWjN0bV4GXD4Ua9ibPXVXw9/GiYmnXyUpKZaj7SfkG7BQpzxZWNnMek1Q6cgfD33yNuLg46taty/z58wkI0Jsulb5uJGnyA2ZrwqRu1upt+3m+TTOO7t5OQEAAIz6fwm6f2wgJCqBNxRBtZVKZ7RzWI5kOuGyvABzN9GiUWxpVKMHprcvp0u9VEhMTadWqFVOmTMHPz8/u0FQ2dCNJ02SgFTA8nWNROcjkNdvp3K4pcacO4OWfh4+/nUfnFvqoL2WrWcAoEWmBNdTAS0SqAB9greGksqAPPviA3r17A/Dyyy/z+eef8+vBcyzcfpQm5UtQqUywzRGq7ORGkiZv4G0RqQPs4OqB4D3TIzCVff2y/R9eavkM8WeP45WnAEVbDmWv1y12h6VUf2AccBirnvvL8e8MYJiNcakUGGN45513GDFiBAB9+vRhxIgRLN5+LHlG3exfDzOrcyVNnFS6uZGk6X5gm+Prci77DEpdx99//02DOk8Rf/YE3vmLUrTVUHwL3kKT8jrlV9nLGBMLdBaRoVh1WyCwzRijD/DNYhITE+nevTtffGEt2Tdy5Ej69OlDWHgUPedsv6Lswu1HNWlS6eZGHtj7ZEYEorK/3377jbp16xIRHo5foVsp3HIo/vkKMaZlea3UlO1EpKoxZoMxJgwIS/UAZYu4uDg6duzIzJkzERG++OILunTpAsCWgxEkuNy66w2ZSk839MBepdJq3bp1NGzYkAsXLvDoo48yYdo89l6wZsjpgG+VRXwvIkeBmcA0Y8xfdgekrhQVFUXLli1ZtmwZPj4+TJs2jVatWiXvr1g6iNy+3kTHJeAt6A2ZSnduJU0isgDoZIw57/j6mowxTdMlMpUthIVH8fmU2Xzc7xViYmJ48sknCQ0NJW/evJS3OzilrnQL0BprLbq+IrIDmA7MNMbog3gzSVh4FDO2HOJwRDTtK5VKTnoiIyNp0KABP/30E7lz52b+/PnUq1fvimNDggNY+UZ1thyM0BsylSHcbWmK5L/xSpEZFIvKZsLCo3jsxcGcWPwhmERq1a3P4oXzyJUrl92hKXUVY8wZrKcLfCYit2Gt2dQRGCEi640xT9kaYA4QFh7F02PWEevoY1u28zizOleiTL5E6taty++//06+fPlYunQp1apVS/EcIcEBmiypDONW0mSMeT6lr5W6ntcHjuBE6P8AyHPfkzzRbaQmTMojGGMOiMhI4A+sZ17WsDmkHGHLwYjkhCnJ5DVb+eGj19m9ezeFCxdm5cqVVKhQwaYIVU6X5jFNIpIbEGNMlON9KaAJ8JcxZlU6x6c8kDGG4cOHE/qFNUs778MNKFizM0cjY22OTKnUOdZmaof1rLlcQCjQz9agcoiKpYPw85bkxCku4igLp7zHyeNHufXWW1mzZg133nmnzVGqnMzrBo4JBToAiEgBYAvwFhAqIq+kX2jKExlj6N27NwMGDAAg/+NtKFjzZUS8aF8p5aeMK5UViMgIETmA9ay5EKAHUMwY094Ys8Le6HKGkOAA1vR8gq41bqdi3kii57/DyeNHueuuu/j55581YVK2u5HZcw8Bbzq+bg6cwHrMQDPgPeCL9AlNeZoFv4XxVo9uHPxlKQAfffQRlRq015V5laeoDvwPmOMY36RsEBIcQNW84YwY/Qrnz5+nQoUKrFixgiJFitgdmlI3lDQFAElPAK8NLDDGJIrIJkCbEnKoeZv307FDe6L+/QXEi9cGfcAbb7wBoMmS8gjGmCp2x6Bg+fLlNGvWjOjoaKpVq8aSJUv4+0w8Q1b/rs+lVLa7kaRpL9BYRBYCdYCPHNuLAOfTKzDlOS5dukT351sT9e+v4O1D4YZvc6JYJbvDUuqGiMi9WN1zVzzx1Riz2J6Ico7Zs2fz3HPPER8fzzPPPMPcuXPZcTya1hM3JZf5+ucDrH6zhiZOyhY3kjS9h/Uspo+AtcaYjY7ttfnv8Soqhzh79iz169fn5N+/Ir65KNx0ALlLl6ddxRC7Q1MqTUTkdmAh1qOiDCCOXUnTubztiCunmDBhAl27dsUYQ5s2bZg8eTK+vr4s3L73inIx8YlsORihSZOyRZoHghtj5mHdhT0COD+Wfi3/jXVSOcCJEyeoUaMGGzdupGDBgoz6ag41nnyKsa3K06iCPrpAeZyxwAGsVvMo4D6scU6/AU/YF1b2N2rUKLp06YIxhq5duzJ16lR8fX2Bqx+D4u9jPUlAKTvc0GNUjDEnsAaAO2/bki4RKY9w4MABatWqxb59+yhevDirVq2iXDnX5zcr5VEqA08ZY86ISCKQaIzZICL9gE+wJryoGxAWHpXiKt3GGPr168eoUaMA6NevH++//z4iklymUplgZnWuxNRNh3RMk7KdPntOpUlYeBSTlq5nbK9ORJ45yW233caaNWu4/fbb7Q5NqZvlzX+TXM5gPVZlN3AIuMuuoDxdWHgUtT76kZj4RPx9vJLHIyUkJNCtWzfGjx8PwOjRo+ndu3eK56hUJlgnlKgsQZMm5baw8Ciq9prA0Vnvknj5An6FSzF78WpNmFR2sQt4EKuLbjPwtojEAi8D++0MzJPN3BJGTHwiYI1HmrkljDdr3k6HDh2YPXs2IsKECRN46aWXbI5UqdRp0qTc9vHUBRyZ0R8TG41f8Tsp0mIwh2L8edTuwJRKH8OAPI6v3wWWAj8B4UAru4LydGERUVe833c8nMaNe7J8+XJ8fX2ZPn06LVq0sCk6pdJGkyblltDQUD59+0VMXCy5Sj1A4SYD8A/IowMyVbZhjFnp9PVe4G4RCQLOGmPMtY9U19O+UimW7TwOQOLli2z5Yhh//LaJ3Llzs3DhQurUqWNzhEq570Yeo6JymKlTp9KsWTPi42LJc2dlijQfjF/uPEx94TEdkKmyNWNMhCZMNydpIPezdwaQa9X7/PHbJvLnz8/q1as1YVIeR1ua1HV17jOUr0a/C0DHjh15d9Qn/H7k/FWzYJRS6lpu8b3E6lFd2fPvvxQpUoRVq1ZRsGRZ5m09onWJ8iiaNKkUGWNo92pvZn75IQB5H25Io9eGcnvRfNxeNJ/N0SmlPMU///xDrVq1OHLkCCEhIaxZswb/oBIpzqhTKqvT7jl1lcTERHr27JmcMOWv2o6CNTsz87cjNkemlPIkv//+O9WqVePIkSPcfffd/Pzzz5QtWzbFGXVKeQJtaVJXiI+Pp3Pnznz77bcAFKz5MvkeaQigj0ZRSrlt/fr1NGjQgPPnz/Pwww+zfPlyChcuDFw9o871vVJZlSZNKllMTAxt2rRh4cKFeHt78/XXX5P//ppM3xJGu4oh+mgUpZRbli1bRvPmzbl8+TI1atRg8eLF5Mv3X7e+84y6pPdKeQJNmhQAFy9epHHjxqxduxY/Pz/mzJlDo0aNADRZUkpdJXTb0RRvqGbOnEmHDh2Ij4/n2WefZc6cOeTOnfuKY5Nm1C3cfpQm5Uvoat/KY2jSpJi2bhc9XmhNxIE/CQwMJDQ0lKeeesrusJRSWVTotqP0mL0dgC0HIgDr5uqLL76gW7du1kSSdu345ptvkh+860ofjaI8kSZNOdw3q36nS7umxJ05hFeuvLz7+QxNmJRS1zXdZeD2tM2H+GvFFPr37w9At27d+OSTT/Dy0rlGKnvR3+gcbP/+/bzWrgFxZw7hHRhE0bYj+TVK7/yUUtfnPCnEGEP8xqnJCdOAAQP49NNPNWFS2VKW/q0WkcEiYlxe/zjtzyUi40QkXEQuish8ESlqZ8yeYteuXVStWpVLZ47hU6AYRduNxq9wKZ0hp5RKVaMKJRjbqjyPlsrPXXtmsGjKlwB8+OGHDB06FBGxOUKlMkaWTpoc/gSKO72qOu37CGgAtABqALcACzI7QE8Suu0oNft+ReUqVTl+/Dj3338/E2Yvo0qFexnbqrwO+lZKuaXefYWRHz5h9cKZeHl5MWnSJHr27Gl3WEplKE8Y0xRvjDnhulFE8gMvAm2NMd87tj0P/C0ilYwxmzI5ziwvdNtRXh71LacXDMPEXeau+x9i3brVBAUF8bzdwSmlPMalS5do1qwZK1euxNfXl5kzZ/LoE/UY/+M+AOqVK64rfKtsyROSprIicgy4DGwE+hljwoCHAV9gTVJBY8w/IhIGVAY0aXIxevxUTs0bDAnx5CpdgXtfHE1QUJDdYSmlPMi5c+eoX78+v/zyCwEBASxatIi7HqrC02PWEZtgPdv4o9X/skofjaKyoayeNG0GOgG7sbrmBgE/iUg5oBgQa4w553LMSce+axIRf8DfaVPedIo3y5o8eTKbJr4DiYkE3Pk4hRr0pmP1u+0OSymPlRPrkZMnT1KnTh3++OMPChQowHfffUflypUZ/+O+5IQJ4HJ8IlsORmjSpLKdLJ00GWOWO73dISKbgUNASyD6Jk7dDysByxHGjh3LG2+8AUDNhq0oUKcb7SvfruOXlLo5OaoeOXToEE8//TR79+6laNGirFq1igceeCDFsr7eQsXS2oqtsh9PGAiezNGq9C9wB3AC8BORAi7Fijr2Xc8IIL/Tq2S6BppFGGMYPHhwcsL05ptvsmrhDOa9Wk0TJqVuXrarR8LCo5i39Qhh4Vc+C+7vv/+mSpUq7N27l9KlS7Nhw4YrEqZ65Yrj72P9d+LjJUx94TFtZVLZUpZuaXIlIoFAGWAqsBWIA2oC8x377wJCsMY+XZMxJgaIcTpvBkVsn1/2nOa1Hj34fflMAIYOHco777yTLb9XpeyQ3eqRsPAoan30IzHxifj7eLHaMSbpt99+o27duoSHh3PPPfewevVqSpS48qYrJDiA1W/WYMvBCCqWDtKESWVbWTppEpEPgCVYXXK3AEOABGCmMSZSRCYBY0QkAjgPfApszOkz5zbsPknd5m25tOt7AHq+O5IBA/rYHJVSKiubuSWMmPhEAGLiE5m5JYzHcp+gYcOGXLhwgUcffZTvvvuOQoUKpXh8SHCAJksq28vq3XMlgZlYA8HnAOFAJWPMacf+N4GlWC1N67G65ZraEGeWcfnyZTq0a2UlTOJF8LNvcfnOWnaHpZTK4sIiruySW79mBXXr1uXChQs8+eSTrF279poJk1I5RZZuaTLGtE5l/2Wgm+OV4124cIFGjRpxYOuP4O1L4cZ9CbjjMW4Nyp36wUqpHK19pVIs23kcgIt//sDK5R+TkJBAw4YNmT17Nrly5bI5QqXsl9VbmpSbwsPDqVmzJj/88AN58gRSovV7BNzxGH7eQtuKpewOTymVxVUqE8yszpUoe+Znwpd+SEJCAu3bt2f+/PmaMCnlkKVbmpR7vl75G28+35Lzxw8QHBzMihUrKHLbvTooUynlNmMMa2Z+wZpJIwB47bXX6DlwOJN+PgToKt9KgSZNHu/LxT/TvUMzEiJP4h0YzODxc3nkkUcAtIJTSrnFGEOvXr0YM2YMAO+++y4vvPY2tT76UVf5VsqJds95sJ07d/Jm+0YkRJ7Ep2Bxij03mvVntBldKeW++Ph4XnrppeSE6aOPPmLIkCGs+PNEiqt8K5WTaUuTh9q0aRPPPPMMl8+fxbdwaYq2HIp3YEHaVQyxOzSllIeIiYmhbdu2LFiwAC8vLyZNmkSnTp1SLKurfCulLU0eafXq1dSsWZOzZ8/y+OOP883cpVS+vwxjW5XXlb6VUle41irfFy9epEGDBixYsAA/Pz/mzZt3RcKkq3wrdTVtafIgm/aFM+rLySwb24+4uFhq167NggULyJMnD+1q2B2dUiqrudYq32fPnuWZZ55h06ZN5MmTh0WLFvH0009fcayu8q3U1TRp8hCb9oXz7GtDCV/xKZhEnqrXiMULZ+Pv75/6wUqpHCmlVb47lM9PnTp12LlzJwULFmT58uU89thjKR6vq3wrdSXtnvMQvQYNJ3z5WDCJBD5QmzvbvKMJk1LqulxX+d71zx6qVavGzp07KV68OOvXr79mwqSUupomTVmcMYaBAwfy83RrZku+ik0JqvsapQoF2hyZUiqra1/pv4VtY8+EsWLUy+zbt4/bbruNDRs2UK5cORujU8rzaNKUhf2y5zSP1G/LsGHDAAh+oiMFnngefx8vXeVbKZWqpFW+qxeM5NL8dzhz8gT33XcfGzZs4Pbbb7c7PKU8jo5pyqI2/HOCus3bcunPHwCh95BRdO/WTQdlKqXSJDpsBwuHdeHixYtUrFiRCdPmMeWPSA7/cIL2lUpRqUyw3SEq5TE0acqCoqOjad+2FZf+XA9e3hSq35NLZWrqoEylVJqEhobSqlUrYmJiqFmzJp9OmkGjCb8lL1q5bOdxZnWupImTUm7S7rks5vz589SrV4+D29YjPn4UbjqAPPfW4Nag3HaHppTyIFOnTqVZs2bExMTQuHFjli5dyp9nYq9Y5Rtg4fajNkWolOfRpCkLOXPmDDVr1uTHH38kMDAvt7R6j4Ayj+LnLTqGSSnltk8//ZQOHTqQkJBAx44dmTt3Lrly5aJi6SD8vOWKsk3K64K4SrlLu+eyiCNHjlC7dm3+/vtvChUqxIoVKyhc+h4dw6SUSlFYeBQzthzicER08tgkYwxDhw5l0KBBAPTo0YMxY8bg5WXdH4cEB7Cm5xNXHaeUco8YY1Ivlc2JSD4gMjIyknz58mX69b8I/YleL7QkKuIEJUuWZPXq1dx9992ZHodSnuj8+fPkz58fIL8x5rxdcWRmPRIWHsXTY9Zd0dU248WKzP18OB9//DEAQ4YMYeDAgYjINc6ilEribj2iLU02+3j2at56sSWJl87hU/AW3v1yniZMSqnr2nIw4oqEySQm0KXLS+z6YTEAY8eO5fXXX7crPKWyLR3TZKNffvmFt19oRuKlc/gWuZ1i7Uax+kii3WEppbI457FJJj6W06Ej2fXDYry9vZk8ebImTEplEE2abLBpXzgtB37BUzWfJi7qAv4l7qVYm+F45ylIu4ohdoenlMriksYmPV+xKIkrRhL970b8/PyYP38+HTp0sDs8pbItTZoy2aZ94TR46wPmDn+dmMvRVKpekylzQql0TwhjW5WnUQWdyaKUSl2gXGb56O4c+XMLgYGBfDtrIX/73Um36b+zaV+43eEplS3pmKZM1m/kWM4sHg0mkYC7q/HgC+/T8vE7aPn4HXaHppTyEMeOHaN27dr8+eefBAUF8c2shfRad5HYhP2ALlqpVEbRlqZM9MEHH7Duq6FgEgl8sC6FGvSidJHMn62nlPJc+/fvp1q1avz555/ccsstrF+/ntig23XRSqUygSZNmcAYwzvvvEPv3r0BKFi5OUF1uuHv66OLViqlUhS67Sgtx28kdNt/yc+uXbuoWrUq+/fv5/bbb2fDhg3cd999umilUplEu+cy2MKth+ndswf71i8EYOTIkbR56TVdtFIpdU2h247SY/Z2ALYciACgWOwR6tWrx9mzZ7n//vtZuXIlxYsXB3TRSqUyiyZNGWj+loN0eL4TUX/9CAiv9B9Bnz59ADRZUkpd0/QtYVe8HzNlAVsn9uPSpUtUqlSJZcuWERQUdEWZkOAA+ta7JzPDVCrH0e65DBIdHc1rL7WzEiYvbwo17M2ZkOp2h6WU8gDOS49E/fsLv4zrxaVLl6hVqxarV6++KmFSSmUOTZoywIz1fxPy4OMc3/kL4uNPkaYDyXNPdV2DSSnllkYVSjC2VXmKHN9IeOhI4uNieaLOs9z/wnB2nYyxOzylciztnktnU77fQee2TYk9uQ/xC6D9oHFEB99Fu4ohugaTUsptB9fP49cp7wNQs2Er9tzVlgPbTzJ/+0ldTkApm2jSlI4OHz5M97YNiD0ZhldAfoq2fI/o4LuY06Wy3aEppTyEMYYhQ4YwZMgQAF565TXW5quN8N/suIXbj2rSpJQNtHsunfz7779UrVqVCyfD8M5bmGJtR+FXtIx2ySmlrst5aYHExER69OiRnDANHTqU2i/0JhFdTkCprEBbmtLBR7NW0e/l1sRcOMudd95J74+msupwgnbJKaWuy3lpgc37TjP23an8sHQeAJ9++indu3cnLDyK3L7eRMcl4C0wpmV5bWVSyiaaNN2kEV8v5J1XnsPERuFXtAz9P59Dx5oP8pLdgSmlsrykpQVMfCynF48mbM8mvL29+fbbb3nuuecAaymBlW9U17XdlMoCsk33nIh0E5GDInJZRDaLSMWMvua7n03lnS5tMLFR+Je8j6JthrNsb1RGX1YplU20qxhCYkwUJ+cOJnrPJrx9/KjT4wPuqFzvinIhwQE0f7ikJkxK2SxbJE0i0goYAwwBHgL+AFaKSJGMumbvkV8wtMfzmPgYct/+CEVaDsHLP4+OYVJKpWrTvnC6Tf+drf+G4b1iGDFhOxC/3BRqMZg/fcvSeuImNu0LtztMpZSL7NI91xOYaIz5BkBEugL1gReAkel5odBtR+nz/hh2z/8IMATcU51C9d8kb0Au3m98v45hUkpd16Z94bSeuIn4C2c4Nftd4sLDyJOvIHkbv4t/8bLJ5XSGnFJZj8cnTSLiBzwMjEjaZoxJFJE1QLrO9Q/ddpRObw7g3I/fAhBYvh5BtboiXt6aMCml3LJw+1Hizh7n5OwBJESexDswmMbvjGdLZG5iE0xyOZ0hp1TW4/FJE1AI8AZOumw/Cdyd0gEi4g/4O23Km9pFjDF0eOUNzm+2Zrbkq9SCAtU7kDeXjyZMSuVAN1KPAExdtp5Tc94l4dJZfAoWp2irYVwIKMqaFx/SB+4qlcVlh6TpRvQDBrlbOCEhgW7duiUnTAWeeJ78jzUD0IRJqZwrTfUIwMaNGzk5oy+JMZfwLVyaoi2H4h1YkBeq3KYP3FXKA2SHgeBngASgqMv2osCJaxwzAsjv9Cp5rZPHxsbSrl07xo8fDwhBdbonJ0yAJkxK5Vxu1yMAq1ev5umnnyYx5hL+Je6haNuReAcWBLQeUcpTeHzSZIyJBbYCNZO2iYiX4/3GaxwTY4w5n/QCLqRULioqisaNGzN79mx8fX2ZM2c2ecvXTd5/cGT99PxWlFIexN16BGD+/PnUr1+fqKgoateuTZGWQ/HOFQhoPaKUJ8ku3XNjgMki8huwBXgDyAN8c6MnjIyM5Nlnn2XDhg3kzp2bBQsWULduXVq0SJ+AlVI5w9dff03nzp1JTEykRYsWTJ06FX9//9QPVEplOdkiaTLGzBaRwsB7QDFgO1DXGOM6ONwtp06dok6dOmzfvp38+fOzbNkyqlSpko4RK6VygjFjxvDWW28B8OKLLzJ+/Hi8vb1tjkopdaM8vnsuiTHmM2NMKWOMvzHmMWPM5hs5T1hYGNWqVWP79u0UKVKEH3/8URMmpVSaGGMYOHBgcsLUq1cvJk6cqAmTUh4uW7Q0pZc9e/bQpEkTDh8+TEhICGvWrKFs2bKpH6iUUg6JiYm89tprjBs3DoDhw4fTt29fRMTmyJRSN0uTJie1a9cmIiKCu+++m9WrV1Oy5HUnwyil1FW6dOnCnDlzEBHGjRvHK6+8YndISql0IsaY1EtlcyKSD4gEePjhh1m+fDmFCxe2OSqllDvOnz9P/vz5AfI7ZrHZwrke8fHxYcqUKbRp08aucJRSaeBuPaItTU4qVqzI/Pnz8ff35/x52+pepVQaZLW/VT8/P6ZNm0adOnWyXGxKqZS5+7eqLU2AiJQAjtgdh1LqppQ0xhy16+JajyiVLVy3HtGkCRBrhOYtXGdxOid5sSrGkm6WV+lDP3d7eMrnnhc4Zmys0LQe8Qj6udvDUz73VOsR7Z4DHB+QW3eoTjNgLtg5fiKn0c/dHh70udsem9YjWZ9+7vbwoM891diyzTpNSimllFIZSZMmpZRSSik3aNKUdjHAEMe/KvPo524P/dwzhn6u9tDP3R7Z5nPXgeAqWxCRb4ECxpjGNoeilPJQWo+o1GhLk8o0IjJYRLbbHYdSynNpPaLspEmTUhlIRPzsjkEp5dm0Hsk6NGlSbhORuiKyQUTOiUi4iCwVkTIuZUqKyEwRiRCRSyLym4g8JiKdgEHAgyJiHK9OIlLa8XV5p3MUcGx7wvHeW0QmicgBEYkWkd0i0uMG4q8qIj85znFYRD4RkTxO+w+KSH8R+VpELohImIi87HKOW0VkjuMziBCRUBEp7bT/WxFZJCLviMgxYLdj++Misl1ELjs+k8ZJ37dY9opIL5drlXeUuSOt36tSWZXWI1qPeDJNmlRa5AHGAI8ANYFEYKGIeAGISCDwI1ACaAg8CIzG+j2bDXwI/AkUd7xmu3ldL6yF0VoA9wLvAcNFpKW7gTsq5RXAfOABoBVQFfjMpehbwG9ABeBz4AsRuctxDl9gJdbibNWAKsBFYIVceSdYE7gLqAU8K9YzyZYAO4GHgIHAqKTCjvV9vgaed4nleWC9MWavu9+nUh5A6xGtRzyXMUZf+rqhF1AIMEA5x/uXsRYHC7pG+cHAdpdtpR3nKO+0rYBj2xPXufZnwDyn998Ci65T/itgvMu2qkACkMvx/iAw1Wm/ACeBro73zwH/4JhA4djmB0QBtZ3iOAH4OZXpCpxJuo5j20vO3zfWStLxQEXHe1/gNNDR7p+zvvSVkS+tR7Qe8aSXtjQpt4lIWUeT+X4ROY9VOQCEOP4tD2wzxkRkwLW7ichWETktIhexKtaQ1I5z8iDQSUQuJr2w7va8gNucyu1I+sJYtc4JoIjTOe4ALjidIwLIBTh3L+w0xsQ6vb8L2GGMuey0bYtzcMaYY8Ay4AXHpgaAPzA3Dd+jUlme1iNaj3gyfYyKSoslwCGgM3AMq6LYhXWXBBB9A+dMdPwrTtt8nQuISGvgA6wm741Yzdq9gcfScJ1AYDzwSQr7wpy+jnPZZ/ivGzsQ2Aq0S+Ecp52+vpSGuJx9BUwVkTexmtRnG2OibvBcSmVVWo9oPeKxNGlSbhGRYKw7nc7GmJ8c26q6FNsBvCQiQde4S4wFvF22JVUSxYFtjq/Lu5SpAvxijPncKZ4ypM3vwL3m5vr1f8caw3DKpO35SbuB50TE3xiTtLjboymU+w6ronwFqAtUv4lYlcpytB5JPofWIx5Ku+eUu84C4cDLInKHiDyFNZjT2UysZuhFIlJFRG4XkWYiUtmx/yBwm2M2RyHHH380sAnoKyL3iEgNYJjLefcAj4hIHRG5U0SGknJlcT2jgMdF5DPH9cuKSCMRcR3AeT3TscYUhIpINRG5TUSecMyeKXmd42Zg/a1NcHyPdYCkGS7Jq8saYxKwxjKMAPYYYzamITalPIHWI1qPeDRNmpRbjDGJQGvgYaym9I+wmrady8QCtYFTWHc7O4G+WIMkwZpxsgL4AevOsI1j+wtYrZ5bgY+BAS6XHw8swJolsxkIxpqRkpb4dwA1gDuBn7DuRt/D6h5w9xxRWHdtYY54/gYmYY1FuOYdo+NusgHWne924H3HtQEuuxSfhNVN8Y27cSnlKbQe0XrE0+ljVJSygYi0w6rQ8jvukpO2VwPWArcaY07aFZ9SKuvTeiTz6ZgmpTKBiHQA9gNHsWbPjALmJFV0IuIPFMaaTj1XKzqllCutR+yn3XNKZY5iwDSspviPsKYAO68S3AZrRlEB4O3MDk4p5RG0HrGZds8ppZRSSrlBW5qUUkoppdygSZNSSimllBs0aVJKKaWUcoMmTUoppZRSbtCkSSmllFLKDZo0KaWUUkq5QZMmpZRSSik3aNKklFJKKeUGTZqUUkoppdygSZNSSimllBs0aVJKKaWUcoMmTUoppZRSbtCkSSmllFLKDZo0KaWUUkq5QZMmlSoRMSIyOJ3PuU5E1qXnOVXWkRG/M8rzaV2i0iqr1SWaNKkMIyL3ishgESltdyxKKc+ldYnKKnzsDkBla/cCg4B1wEGXfbUzOxiVqXID8XYHobINrUtyrixVl2jSpGxhjIm1O4asRkQCjDFRdseRHowxl+2OQeUMWpdcTeuSjKPdc5lMREqIyCQROSYiMSJyQES+EBE/R/OzSeGYTo5+3dJO2w6KyFIReUJEfhORaBHZKSJPOPY3dby/LCJbRaSCyzlTHAcgIt+KyMFUvodSIvK5iOx2XDdcROa6xNcJmOt4+4MjfuMUX/L1RaSoiMSLyKAUrnWX47juTtsKiMjHInLY8RnuFZE+IpLm32cReUxEVohIpIhEiciPIlLFpcxgRwx3OD6fc47y34hIQArnfM7xmUeLSISIzBKRW13KrBORXSLysIisF5EoYLhjX7CITBWR845rTRaRBx0xdHKUed7xvkIK1+8vIgkiUsLNzyBIRD5w/L5cdFx3uYg8mELZIo7f35OO360/RKRjCuWy1DiE7EjrkuT4tC7RuiTTaEtTJhKRW4AtQAFgAvAPUAJoDlz1B+OGO4AZwHhgGtALWCIiXbH+aD53lOsHzBGRu4wxiTfzPTg8CjwOzAKOAKWBV4B1InKv4w5nPfAJ8Lojlr8dx/7tejJjzEkR+RFoCQxx2d0KSMBRaToqlh+xPrfxQJgjlhFAceANd78JEXkKWA5sdVw3EXge+F5EqhljtrgcMgc4gPV5PgS8BJwC+jid8x1gqKPsV0Bh4DVgvYhUMMacczpfsOP6s7B+ficdlfUSoCLwBdbvSCNgskss84BxQDtgm8u+dsA6Y8xRNz+K24HGWJ/xAaAo0AX40fHzPOb43nJjdY/cAXzmKNsC+FZEChhjxrp5PXWTtC4BtC4553Q+rUsyizFGX5n0wvplTQAeSWGfAIOtH8lV+zoBBijttO2gY1tlp221HduigBCn7S87tj/htG0d1h+D67W+BQ66bDPAYKf3uVM4rpKjXHunbc1dr3ut6zvFWM6l3J/AWqf3A4CLQFmXciOw+r1vdfNnIcC/wApAnL83YD+wymnbYEdsk1zOsQA44/S+lCOG/i7lygFxztsd378BuriUberY3sNpmxew1rG9k9P2GcBRwMtpWwXXcm58Fv7O53BsKw1cBgY6bevhOHc7p22+wC/ABSDvtX5n9JW+L7Qu4VrXR+uSpO1al2TAS7vnMokj628MLDHG/Oa63zh+O9LoL2PMRqf3mx3/fm+MCUth++03cI2rGGOik74WEV8RCQb2Auew7ppuxAKsSqKV07nLYQ0Ane1UrgXwE3BWRAolvYA1gDdQ3c3rlQfKYlUWwU7nyYNVqVRPoYn+S5f3PzmOzed43xSrUprjEtsJYA/wpMvxMcA3LtvqYlWKE5M2GOuOflwK38MU4BaX87YDooH5KX3TKTHGxDiugYh4O36eF4HdXPnzfMbxvcx0OjYOqxUgEKjh7jXVjdO6JFVal1i0LskA2j2XeQoD+YBd6XhO58oMY0ykiAAcdikX6fi3YHpc1NG02g+r+bkE1p1Wkvw3ck5jzBkRWYvVrD7QsbkVVuW3wKloWeAB4PQ1TlXEzUuWdfzr2lTtLD9w1ul9mMv+pH0FgfOOcwpWpZaSOJf3R83Vg1hLAcfN1YM496ZwvtXAcazKba2jYm4DhBpjLlwjhqs4jusBvArchvUfRpJwl9j2mKu7Zf522q8yntYl16F1STKtSzKAJk1Zy7XuEL2vsT0hjdudKyTj8j61azn7FKuS+xjYiFWRGqz+9JtpvZwFfCMi5Y0x27EqvbXGmDNOZbyw/sBHX+Mc/7p5raQ4ewPbr1Hmosv71D5XL6zPod41yrqeLzqFMm4zxiSIyAygs4i8ClTBuluclsZT9ccaO/E11n8yEVhjMj5GJ4t4Kq1LtC5xm9Yl7tOkKfOcxrqDKHedMmfBmtFhrhzklxFZ91lSbmJ351rNgcnGmLeSNohILqxBqc7S2k2wCGtAZivHXe6dWOMLnO0DAo0xa9J4blf7HP+eT4dzOZ9TgAPGGHcrXFeHgCfl6inDd1yj/BTgLaABVgV7GliZxms2B34wxrzovFFECgDO/8kcAh4QES+XO8S7nfarjKd1SeoWoXWJ1iUZIFtkfp7A8YuxCGggIo+47hfrLzvpj6+60/Y8QMcMCGkfcLeIFHa61oNYdxipSeDqO8vXuPrO8pLj3wLuBOSo3Fdi3RW2BmKxPjNnc4DKIlLH9Xixpg+7eyOwFesz6CUigSmcq/DVh6RqAdZnM8jx83Q+nzj691OzEmtAZGenY72AbikVNsbsAHZgzb5pBswyxqR1Ibirfp4i0gKru8TZd0Axrhwr4oP1s7+INRNJZTCtS1KndQmgdUmG0JamzNUfa1bKjyIyAav/tjjWgMSqwCqsvu5JIvI/rF/AF7Ay/pB0juVroCewUkQmYfXfd8WaYZLvegcCS4H2IhIJ/AVUBp7myj5rsJqqE4A+IpIfa7Di98aYU9c592ysJuFXgZUud8kA/wMaAktF5FusCisPcD/WXU5prryjSZExJlFEXsKapvuniHyDNXukBNZgyPNYd1xuM8bsE5EBWHe0pUVkEdZMkNuAJlhTwz9I5TSLsKaSfygid2BNE24IBCVdJoVjpjidN63N6WD9PN91fAa/YH2W7bBm/jibgDV9+FsReRhr1lVzrP8c30jL2Ad107Qu0bpE6xI72D19L6e9sCqsyVhrclzGukP5DPBz7H8I2IRVKRwC3uTa04SXpnB+A3zmsq20Y3svl+3tHNePwVqfozbuTRMugFVRnsb6Q14B3OWI6VuXY19yXCMepynDXHuacl6sac5XTEd1KROItV7LHkfsp4GfsZqWfdP48yiPNTvkjOPncRCrsn3KqcxgRzyFXI696ufi2N4UazbMRcfrb8fP+E6nMuuAXdeIqRAwHauyPYc1K+Zxx7VapVC+mOPz3X2Dv5P+WBXlMcdnvwFr2vdVPyOs/xCTfvYxWHemna7xezj4RuLRl9s/N61LjNYlaF2SqS9xBKWUysJEpDGwEKhqjPnZZV8hrJkv7xljhtoQnlLKQ2hdcnN0TJNSWYxjGrbze2+svv7zwO8pHNIJawzI1AwPTinlMbQuSX86pkllOyISBPhdp0iCMeZaa7NkBZ86KruNWM3dTbGa1PubKxcDfAprwb53gEXGmIPOJ3GcI7W1biKMPvBUqRRpXZK8X+sSB+2eU9mOWA/vvN6KsoeMMaUzJ5q0E5G2WGMq7gByYS1G94Ux5jOXcuuwKsCfgeeMy/OhxHogp+sqwa6eNMasS5fAlcpmtC5J3t8JrUsATZpUNuSYjXG9FYujXfvysyMRKQ7cl0qxrcaYs6mUUSpH0rrEonXJfzRpUkoppZRygw4EV0oppZRygw4EJ3kF3Vuw1glRSnmevMAxY2PTudYjSnm8VOsRTZostwBH7A5CKXVTSmKtxGwXrUeU8nzXrUc0abJcADh8+DD58qW26r9Syk4JCQn06NGDqVOtpWTeffdd3nvvPbC/hUfrEaU8yMaNG2nRogUXLlzg3nvv5a+//oJU6hEdCA6ISD4gMjIyUis7pbKwmJgYnnvuOebNm4eXlxcTJ06kefPm5M+fHyC/Mea8XbFpPaKU51i+fDnNmjUjOjqaqlWrMmPGDEJCQiCVekQHgiulPMKlS5do2LAh8+bNw8/Pj7lz5/LCCy/YHZZSysPMnj2bhg0bEh0dTb169Vi5cmXSjVeqNGlSSmV5Z8+epVatWqxatYo8efKwdOlSmjZtandYSikPM2HCBNq0aUN8fDytW7dm0aJFBAQEuH28Jk1KqSztxIkTPPHEE2zcuJECBQqwZs0aatWqZXdYSikPM2rUKLp06YIxhq5duzJt2jT8/K73lJyradKklMqyDh48SLVq1dixYwfFihVj/fr1VKpUye6wlFIexBhD37596du3LwD9+vXj888/x9vbO83n0tlzSqks6a+//qJWrVocO3aM2267jdWrV1OmTBm7w1JKeZCEhAS6devG+PHjAau16e23377h82lLk1IqSwnddpRHe3xJ+Ucrc+zYMe69915++uknTZiUUmnSf+7vBD34FOPHj0dEmDBhwk0lTKAtTUqpLCIsPIohi3exbPVaTs0fiomNxq94Wfp+NpsSJUrYHZ5SyoPUHLmCXyYO4PL+reDlQ5OeI+jcufNNn1eTJqWU7cLCo3h6zDrO/bOJ06EjISGOXKUeoHCTASz59wLtn7Q7QqWUp2gxdjU/j3uLmCN/IT7+FG7Sn90B5dLl3Jo0KaVsFbrtKB+v3UPEju8JX/YRmERyl61E4YZvIz5+tKsYYneISikP0fPbdYQO70rcqf2Ifx6KNH+XXCXvo0mF9Gmt1qRJKWWb0G1H6TF7O+e3LuHsGmugZp5yTxFcrwfi5U37SiE0SqfKTimVvb0+cRVf9H2B+IijeAUUoGjL9/Arejt3FsnD0Mb3p8s1NGlSStlm2uZDnPt5JpEbpgOQ9+EGFKzZGREv2lcKSbeKTimVvb3y2WIm9n+RhAtn8M5XmKKthuEbVIIqZYKZ3jn9linRpEkplelCtx1l2qaDHF4+Pjlhyl+lLfmrtEFEuL1Q+t0ZKqWytzEzVjDh7fYkRp/HJ6gkRVsNxSdf4Qy58dKkSSmVqUK3HeX1mVsJX/4pl3atAeDWZ17B6/76yWV61CxrV3hKKQ8y/Kv5DOzegcSYKPyKlqFIy/fwDsifYS3VmjQppTLVlJ/3ciZ0FFH//gLiRXC9HnjdX5P2lULYffIi7SrqOCalVOoGfjKF99/qjImPxf/WchRp9i5e/gEZ2rWvSZNSKsOFbjvK9C1hNC0XzL5pA62EyduHwg37EHBnZQB2n7zInC6VbY5UKeUJ3hoxjjEDekBiArnLPEqhRn3JFxjA+43vz9CbLk2alFIZKmmGXEL0BUKHPE/s8d14++WiTNshxBS9L7mcLi2glHLHK++M5Mvh/QFDwL01KPTMm4i3T4YnTKBJk1Iqg03fEkb8xQhOzR5I3JlDeOUKpHCLIcQUvUu75JRSbjPG0OH1/kz7bCQAgRXqE1SrC3lz+WZKwgSaNCmlMtjTJWBh37eJP3cC78AgirR8D7/CpQHtklNKuWfR70d4vWcvDv84G4B8lVtRoNpziEimJUygSZNSKgMkjWGqFhzNiNefI/7cCfIUuoVOQ79i6cHE5HLaJaeUSs2C38Lo+GJnLu5YBUDBJ18gX8WmBPp7Z2rCBJo0KaXSWdIYpphju5k/dzCJly9Qrlw5Vq1aRfHixanlSKi0S04plZrY2Fh6dHmeizu+B/EiqE538j5YGyDTEybQpEkplc6mbwkj+tAfnF4wDBMbTb6Qeyjz/AdsOZFIo+LQqEIJTZaUUqma/cseenRuz8m/NoOXD4Ua9CLP3VUpUygPr9csa0s9okmTUipdlbn0F3PnDoKEeHKVepD8TQew/VQ8PWZvB9CESSmVqunr/+LFts2JOfo34utP6/5jiS/+gO0t1Jo0KaVu2KZ94UzddIiQoADaVAxh3bJ5/O/tLpCQQInyNbijVX8OnotLLj99S5gmTUqp65q8Zjtdn2tGzMn9ePnnoXDzwcQXfyBLTBrRpEkpdUM27Qun9cRNye8//PhjTq8aD0CnTp2YOHEiy3aeTG5hAh34rZS6vonLNvNq+ybEnz2OV54CFG05FL8it2WZukOTJqXUDVm4/ShgrZ0S+fNMIn+eAUDZp1rRsPt7+Pj4JLcq6cBvpVRq/v77b3o814D4c6fxzleEoq2HUbDYrbYM+L4WL7sDUEp5piblS2BMImfXTkhOmPJXbUfMI8/x5twdhG6zkqpGFUowp0vlLFPpKaWyng+nf0f5RysTfe40vsG3Uuy50fgWvCVLJUygSZNS6gY9Uio/5fbO4MLWJQCUa96DAlXaICKA1bqklFKpGTZxLr1faE7spUj8ipXltTHTePyBOxnbqnyWSphAu+eUUjfg8uXLtGnThu8WLcLb25tvvvmGfOWe0vFLSqk0eefjbxnR62VMQhz+IfdTpOlADkf7ZolB3ynRpEkplSrnWXIN7inIay+25fvvv8ff3585c+bQsGHD5LI6fkkp5Y6e73/KRwPfAJNI7jsqUqhhH7x8/bP0DZcmTUqp63KeJZcQfZ7BLw3m8rF/CQwMZPHixTz55JPJZXXhSqWUO8aNG8dHA14HIM99TxJcrwd5A/yz3BgmV7aOaRKRfiLyq4hcEJFTIrJIRO5yKZNLRMaJSLiIXBSR+SJS1KVMiIgsE5Eox3n+JyKaECqVDpJmycVfCOfkjL5cPvYvufMW4LHuH3O+wJ02R6eU8iTGGIYNG0b37t0ByPvQswTXfxPx9snyCRPYPxC8BjAOqATUAnyBVSKSx6nMR0ADoIWj/C3AgqSdIuINLAP8gMeBjkAn4L2MD1+p7K9J+RLEnT3OyelvE3cmDO/AYAq0HM5eU4wes7cnz5JTSqnrMcbQuENXBg4cCMDAgQOZMvELHru9UJYc9J0SW1tjjDF1nd+LSCfgFPAwsF5E8gMvAm2NMd87yjwP/C0ilYwxm4DawL3A08aYk8B2ERkIjBKRwcaY2Mz7jpTKfvJEHSNm4TvER54iuHgI5bt8wN7ogOT9usq3Uio1CQkJ1GnWjrWhswEo+NRLPNykC40qlKDxQyVtjs59drc0ucrv+DfC8e/DWK1Pa5IKGGP+AcKApKH1lYGdjoQpyUogH3BfhkarVDa3adMmatSoQfjpU9x///3s+n0zrzV8/IoyWXnQplLKfvM276fUo09bCZN4EVyvB/kebeyRy5JkmXE/IuIFfAz8bIzZ5dhcDIg1xpxzKX7SsS+pzMkU9uNUxvVa/oC/06a8Nxa1UtlL6LajybPf8oT/TePGjbl06RKVK1dm2bJlFCxYkEaOv6qcPktO6xGlUjf7lz10ateKywe3gbcPhRu8TcBd1o2XJ95wZZmkCWtsUzmgaiZcqx8wKBOuo5THCN12NHmdpXUrlnJ22QfEx8VSq1YtFi5cSJ48/w011FlygNYjSl3X2bNn6da+KZcP7kJ8c1G46QDue6QKhfL6e+wNV5bonhORz4BngSeNMUecdp0A/ESkgMshRR37ksoUTWE/TmVcjcDqCkx6eU6HqlIZJKmp/OKO1ZwOHUl8XCyVaz5D3obvsObfc/YGlzVpPaLUNZw4ccLq2t+/C69cgRRtNYzcpcvTo2ZZj36skt1LDogjYWoCPGWMOeBSZCsQB9R0OuYuIATY6Ni0EbhfRIo4HVcLOA/8ldJ1jTExxpjzSS/gQrp8Q0p5sHYVQzj/6yLCl48Fk8iDNZtw9KEubD1yUWfJpUDrEaVSNmHpRu544BF27txJsWLFGDN5IdWqPu4xM+Sux+7uuXFAW6ARcEFEksYgRRpjoo0xkSIyCRgjIhFYidCnwEbHzDmAVVjJ0VQReRtrHNMwYJwxJiYzvxmlPJUxhq0Lx3P2+68AaNS+C76Pd+DXg2eTy+gsOaVUaj6b/wNvdGpOwsUIfPIXZdD4eXRtWIUedgeWTuzunnsFq1l7HXDc6dXKqcybwFJgPrAeq8utadJOY0wCVtdeAlar0zRgCvBuhkevVDaQmJhIjx49GDp0KADDhg1j4eQveO6xUleU88RBm0qpzPPrr7/yVsfGJFyMwDc4hKLtRvP9cbvTjPRl9zpN4kaZy0A3x+taZQ4Bz6RjaEplS86z4xpVKEFcXBwvvPAC06ZNA6xHG7z66qsAya1KOX2WnFIqdcMmzGHI688THxOFX/E7KdJiMN6582W7my27u+eUUpnEeXbclgMRxMZcZtqInixevBhvb28mT55Mu3btrjhGZ8kppVLT/6OvGdG7KyTEkavUA7w2YjwHz5tsebOlSZNSOYTzQnKJMVF069iS0//+Tq5cuZg7dy7PPvusjdEppTzR1KlTGdnrZUhMIHfZShRu+DYHzxvmdKmc+sEeKHt1NiqlrimpmTwhKpKTs97h9L+/kzdvXgZ8OpUpR4N1dpxSKk0+/fRTOnTogElMIE+5pyjcuB/i45ftuuScaUuTUjlEowolCD91nDc6dSf2xEGCg4PpN3Yqn+5MBCLYciAiuZxSSl2LMYZ2r/Zm5pcfAvD666/zRIdezPztSLbsknOmLU1K5RB79+7lvS7NuXDiICVLluSnn35i88UCV5TxxGdBKaUyT2JiIo2eezk5YcpfpS1PduxNk4dv9ehFK92lSZNSOcCOHTuoWrUqhw4d4o477mDDhg3cc889VzWjZ+dmdaXUzYmPj+fFF19kyQxrPbeCNTtToGpbZvx62ObIMo92zymVjWzaF87C7UdpUr4ElcoEA/DLL79Qv359zp07x4MPPsjKlSspWtR60pAuK6CUcsfly5d58pkmbPphBeLlTVC91wksZz2sIyfdbGnSpFQ2sWlfOK0nWgvlz/71MLM6V+L8vq00adKEqKgoqlSpwtKlSylQoMAVx+myAkqp67l48SJVatZjx5YN4O1DoYZ96NKhFbtPXsxxN1uaNCmVTUzddOiK94M/mcT3XwwgLi6OunXrMm/ePPLkyWNTdEopTxQREcEzzzzDji2bEd9cFG46gNyly7P75MVsu6zA9eiYJqWyiVuDcid/feGPVaz6rB9xcXE8XutZ8jzbjzX/nrMvOKWUxzl+/Dg1atRg8+bN5M1fgKKt3yd36fJAzuqSc6YtTUplE20rluLrDQc4s3E+Z3/4GoBaTdqw+47WHD18gd8cq4HnpKZ0pdSNGb/kF956oQWXzhyjePHirFq1in1xBXP8+EdNmpTKJm4Nyk3t6HWMcyRMvXv35sDtjZGDZ5PLTN8SlmMrO6WUez6Zu5aeL7Qk4WIEPgWKMWj8PMqVK0c59KZLu+eUygYSExPp3r074z76HwAjRoxg9OjRPPdYqSvK5dQmdaWUezZv3kyvTk1IuBiBb6FSFG03mrXHxO6wsgxtaVLKw8XFxdGpUydmzJiBiPD555/TtWtXQJcUUEq5J3TbUcZMWcCW8X2Ji47Cr/hdFGkxGO/cefVmy4kmTUp5kNBtR69IgKKjo2nRogXLli3Dx8eHKVOm0KZNmyuO0SUFlFLXE7rtKJ2Hfs7pxaMhIZ4HH6tG3w+/YsGucL3ZcqFJk1IeInTbUXo4BnNvORBB1MULfDmgC+vXrydXrlzMmzeP+vXr2xukUsrjvD/2S04vGgkmkYA7H6fMc0NpXeVOWlexO7KsR5MmpTyE83PhEqIiefW5xpwL202+fPlYunQp1apVszE6pZQnGjt2LL9OHgZAnnJPE1zvNTpUucPmqLIuTZqU8hDtKoaw5UAE8edPc3L2QOIjjlC4cGH6fjyFT//yISLwqDajK6XcYoxhyJAhDBkyBICG7TrjV6Ujz1UqrfXIdWjSpJSHaFShBMfC9tPrhReJjzjJrbfeSt+x0xi9+QIQwZYDEcnllFLqWhZuPczbvXuy94d5AAwdOpR33nkHEZ0llxpdckApD7F9+3YGv9ycqIiT3HnnnWzYsIF1p3yvKOPchaeUUq4W/HqI9p06JSdMnd8eyoABAzRhcpMmTUp5gJ9//pknnniCU6dOUb58eX766SdCQkKumgqsU4OVUtdy+fJlXuvcnku7vgfxIvjZtzh3e027w/Io2j2nVBa3YsUKmjZtSnR0NFWrVmXp0qXkz58f0HWYlFLuuXDhAo0aNeLYHz+Bty+FG/UloOxjeqOVRpo0KZVFbNoXzsLtR2lSvgSVygQDMGfOHJ577jni4uKoV68e8+bNIyAg4IrjdB0mpdT1rPxtDx1aNeHU/j/JkyeQvmMmscPcqjdaN0CTJqWyAOc1mGb/ephZnSux8/sFdOnSBWMMrVq1YsqUKfj5+dkbqFLKoyz+eRfNG9UnLjwMr9z5CGr5Hh2aPUtIcEDqB6ur6JgmpWwWFh5Fzznbr9jWd/AwXn75ZYwxdOnShenTp2vCpJRKk3379tGhaR3iwsPwDgymaNuReBW5gy0HI+wOzWNpS5NSNttyMIIEY31tjOHc+skc2mTNbOnbty/Dhw/XmS1KqTTZsWMHderUIfLUCXwKFKdo62H45C+Kn7dQsXSQ3eF5LG1pUspmFUsHkdvXG5OYwNnVn3PekTCNGjWKSi2702rCJkK3HbU5yuxJRF4XkVyOr0NEs1OVDWzcuJEaNWpw4sQJHnjgARYtX0OjquXpWuN21vR8QrvmboK2NClls5DgAJa8WolOnToStm05IsL48eMp8sgzVzxrDnThygwwBpgFXAYOAMWBU7ZGpNRNGPLFTIa+8QIJsZepXLkyy5Yto2DBgtSvZHdk2YO2NClls6ioKF5/sS0/r1qMr68vs2bNonPnzlctVKkLV2aIY0AzESkFCFDS0eJ01cvmOJVKVd8PJjL4tQ4kxF4mV+kK9PhgMgULFrQ7rGxFkyalbBQZGUndunVZvnw5uXPnZvHixbRs2RK4eqFKXU8lQwwDPgb2Awb4FavFyfl10PGvUlnWN998w6i3u0JCPAF3VaFIs3eZv/OM3WFlO9o9p5RNTp06Rd26ddm2bRv58+dn6dKlVK1aNXm/LlyZ8YwxE0RkJlAK2AE8DYTbG5VS7gvddpRBw0fzx7xPAAh8oDZBdbohXt56o5UB3EqaRCSfMeZ8RgejVHYUuu3oVYlPWFgYtWrV4t9//6Vw4cKsWrWK8uXLX3WsLlyZsUTkdWCCMWaXiDwPbDTGRNsdl1LuWPT7ETp170XkxtkANGrfhU5vDGDGr4f1RiuDuNvSdFZEihtjTonI90BTY8y5DIxLqWzBedHKpMHcdwdcpFatWhw+fJiQkBBWr17NnXfeaWOUOZrzQPCvgeWAJk0qy0tMTKT3W28SudGabVugegd8H+9A44dK0vihkjZHl325O6bpIhDs+PoJwPfaRd0nItVFZImIHBMRIyKNXfZ/69ju/FrhUiZIRKaLyHkROScik0QkMD3iU+pmuQ7e/mzeaqpVq8bhw4e566672LBhgyZM9tKB4MrjxMXF0bFjR/aumwcIQbVeIX/lljz3WCm7Q8v23G1pWgP8ICJ/O94vFJHYlAoaY55Kw/XzAH9g3eEtuEaZFcDzTu9jXPZPx5omXAsrmfsGmAC0TUMcSmWIdhVDkluYLh/exS+L3yfq4gUeeughVqxYQeHChW2OMMcbBnwKfMZ/A8FdiWOfdybGpVSKoqOjadWqFUuWLMHb25vX3/uYI8EPa3dcJnE3aXoO6AiUAWoAfwJRN3txY8xyrObw6614HGOMOZHSDhG5B6gLPGqM+c2x7TXgOxHpZYw5drMxKnUzkiqxDybNYsv8QcTGxFC9enUWL17Muv0Xmb5go1Z2NtKB4MqTnD9/noYNG/Ljjz+SK1cu5s6dy7PPPmt3WDmKu0mTrzHmSwAReQTok4ljmp4QkVPAWeB7YIAxJqlSqwycS0qYHNYAicBjwMJMilGpa4re/RObxvclPj6e+vXrM3fuXFb9E6ELV2YRxpgLQNJA8J+NMa6t2UrZbsWv/9KhVRNOH/iLwMC8LF26hBo1atgdVo7j7pimsyJSxPG1yahgUrAC6ADUBPpgtXItF5GkZvJiuKzea4yJByIc+1IkIv4iki/pBeTNiOCVGj9+PG3btiU+Pp62bduycOFCcufOrQtXZkHGmMlpSZi0HlGZZfGGnTSs9zSnD/yFV+58BLUcxm3lHrU7rBzpRgaC1yCdBoKnxhgzyxiz2Biz0xizCHgWeBRrMPrN6AdEOr2O3OT5lLrKyJEj6dq1K8YYXnnlFaZOnYqvr/WnowtXZg0iEiEihRxfn3W8T/GVwuFaj6gMt2fPHto3rUtc+GG88xaiWLtRSOEybDmY0q+kymg3MhBcSL+B4GlijNkvImeAO4C1wAmgiHMZEfEBghz7rmUE1lTjJHnRCk+lE2MMffv2ZfTo0QD079+fYcOGXTFuTxeuzDLeBC44fZ2WlnStR1SG+uOPP6hTpw7nT5/Ep+AtFG01DJ/8RfD38aJi6SC7w8uRbB0InlYiUhKrxeu4Y9NGoICIPGyM2erY9hRWC9rma53H0QSf3AyvDzZXN2LTvnAWbj9Kk/IlqFTGaohNSEjglVdeYeLEiQD873//o1evXikerwtX2s8YM9np62/TeKzWIyrDjPo2lIHd2hMXdYEHH3yQ4eNnsmxvNCFBAbSpGEJIcIDdIeZIbiVNjhVy030guGM9pTucNt0mIuWxxiRFAIOA+VitRmWA0cBeYKUjrr8d6zZNFJGuWN2GnwGzdOacykib9oXTeuImAGb/ephZnSvx0K15ad++PXPmzMHLy4sJEybw4osv2hypuh7HWCS36FMRVGYZPG46773xIiY+Bv8S99L7kxk889g9PPOY3ZGpND97zhjzZDpe/xHgB6f3SU3dk4FXgAewWrgKYC1CtwoY6DJYsx1WorQWa9bcfOD1dIxRqatM3XToivdf//gPQ+YOY8WKFfj6+jJjxgyaN29uU3QqDc7hfpecrtOkMtzcuXN5r0cnTEI8uW57mMJN+hH6dyTtqtsdmQL3nz03BitZueT4+pqMMT3dvbgxZh3WGKlrqePGOSLQhSxVJrs1KHfy14mXL7Lqw3c59NfvBAQEsGDBAurUSfVXV2UNzjeBpYGRwLdYXf9gLWvSEWvQt/p/e3ceZ1P9P3D89Z4xixlbM3YaJPUlRSn0zVI0luxF1iJLKcqWohVJRELaKEL2ZBvbDL9vpGxhlKSyjn2ZwWDMYubz++Pc0TUNrjEz596Z9/PxOI+595xz733Pwcf7fs7n8/6oLPXVV1/xwgsvYFJSCPhPbQo37Y94++hEETfiak/T/fwzY+7+65yXneUIlLJNh+plmLJ+P5diYzg57x0ST+6nUKFCLFu2jFN5y/D0l1q00hMYY9amPhaRd4D+xpjZTqcsEZHfgOexesCVyhJjxoxh4MCBAPTo0YNGz7/FnK1HtB1xM66OaXosvcdK5VYhwQF80+YOWrfoTeLJ/RQrVozw8HD2Jwdr0UrP9TDQM539vwBfZXMsKpfYsOc0vQe8zrYlUwB47bXXGDlyJCLCkw9qD5O7cbVOk1LKye7du2nbLJRjUfspU6YMP/74I/fdd58WrfRsh4Ae6ezv7jimVKb6+e9TNHi685WEqfBjXej12rs6E9ONuTqm6VqL6f6LMebJjIejlPvbunUrjRo14vTp01SsWJGIiAhKlbJ6k5wX6E19rjxGP2CBiDTmn5Il1YEKwFO2RaVypKSkJLp37cKF7csBIajhSwRWbczmAzFaTsCNuTqm6ZzTYwFaOfalrvlWDWuGm8vJlVKeaO3atTRr1ozz589TrVo1Vq5cSeHCha8c16KVnssYs1xEKmDN3K3o2L0U+MIYoz1NKtNcunSJNm3a8Mf65eDlTeEm/QmsVFeLVnoAMebmxm6LyCisits9jTHJjn3ewGdArDFmYKZHmcUctVrOnTt3jgIFXC7bonKZsLAw2rRpQ3x8PHXr1mXJkiX698UNxMbGUrBgQYCCdtZS0nZEuWLWuj/o070Dp/+OxN/fn/cnTmWv/11atNJmrrYjN12nCegK1EpNmACMMcmOUgQ/Ax6XNCmVauPeaGZsPPivBmzWrFl07tyZy5cv06xZM+bOnUvevHlv8G5KKfWP6f/3Kz06PEniib2IbwBvT5xB/246osWTZCRpygP8B/gzzf7/oAPLlQdzrvINMOWn/UT0q0vY3G/o3bs3xhg6derElClTriy8q5RSrjh06BC9OzQj8UQUXgEFKdpmKJGXS9gdlrpJGUmapgJfi0h5YLNjXw1gkOOYUh5pYeSRq57HJyXz+ttDmPP5aAB69erFhAkT8PLS7wZKKdf99ddfhIaGcv5EFN75i1Cs7Xv4BJfWiSIeKCNJ06tYa8ENAFLT5GPAaOCjTIpLqWzXqmop5m6xxvsaYzi/dipzNllzG95++22GDh2qU4GVUjclMjKShg0bcvLkSe666y4GfjyD8EPJOlHEQ930QPCrXuxY7NLTF7LUAZwq1ca90Uz7eR+/fDuKX8IXADB27Fj69evH4u1HdFacG8rMgeAikgd4FGuB8FnGmPMiUhJrksuFG7xW2xF1lfXr19OkSRNiY2O5//77WblyJUWLFrU7LJWOrBwIfoWnJ0tKpXV/6XyMWWwlTF5eXkyePJmuXbuyePsRrfSdw4lIGWAlEAL4ARHAeeB1x/P0qoUrla53Js5gRP8eJCclUKtWLcLCwlL/U1YeTAdnKOVw8eJFmjdvzoIFC/D19WX+/Pl07doV+Hdlb630nSONx6o9dxtwyWn/QqC+LREpjzRw5Oe81+c5kpMS8L+jGi+PmqoJUw6hSZNSwJkzZwgNDSU8PJzAwEDCwsJ48sl/pgKnHbCpAzhzpNrAcGNMYpr9BwDtVlQumTRpEmPe6AUpyQRUrEPRJ9/iu99O2R2WyiS3dHtOqZzg+PHjNGzYkF9//ZVChQqxYsUKatasedU5Wuk7V/ACvNPZXxrrNp1S1zVq1CgGDRoEQL6qjQkK7Yl4eeuXrBzklpImEfE3xsRnVjBKZbcDBw4QGhrKnj17KF68OOHh4dx7773pntvi/lKaLOVs4UBf4HnHcyMi+YChwHK7glLub8Oe0/Tq9yrbw6YBMHjwYKq3folZWw7pl6wc5qaTJhHxAt7EGhRZTETuMsbsE5H3gAPGmK8zO0ilbkVUdByzNh/kUMwlnqlZhprlgwHYtWsXoaGhHD16lHLlyhEREUH58uVtjlbZaACwSkR2Af7ALKzFek8D7e0MTLmvn/46SaO2nbkQuRKAwvWeo+eAtwgJDqDlA6Vtjk5ltoz0NL0FdAZeAyY77d+J9S1NkyblNqKi43h87A8kJlulNZb9dow5PWriHbOPxo0bEx0dTaVKlQgPD6dUKf02mJsZYw6LSBWgHXAfkA+rPZtpjLl03RerXCkxMZHuzz3LhchVgBDUsBeBVRux+UCMriGXQ2UkaXoWeN4Ys0ZEvnDavwNrKRWl3MbmAzFXEqZU479dRNiYvly4cIGHHnqIFStWEBwcbFOEyp0YYy4D39odh3J/cXFxtG7dmt0/rwKvPBRuOoDAirXxy+NF9bJBdoenskhGkqZSwJ509nsBuiCXcivVywbh6y1XEqe4vzfxfdiHJCYmUK9ePRYtWkT+/PltjlK5AxFpfo1DBogH9hhj9mdjSMpNnTt3jqZNm7J+/Xry5s3L+xO/YY/fnf9a6FvlPBlJmnZhTc09mGZ/a2D7LUekVCYKCQ5gdf9HmbX5IKuXLOCHxSNITk6mRYsWzJkzh1V/RDNz804drKkAFmElSGnXykndZ0RkPdDSGHMmm2NTbuLkyZM0bNiQyMhIChYsyLJly3jkkUfsDktlk4zUaRoGTBSR1x2vf1JEJmMNDh+WmcEplRlCggPIt3cNa754m+TkZJ599lm+++47Vv0RTZ+5kWzeH0OfuZEs3n7kxm+mcrJQYIvjZ0HHFgpsApoCdYBgYIxdASp7fbV8M+Xve4jIyEiKFi3KDz/8oAlTLnPTSZMxZjHQDHgcuIiVKFUEmhljIjI3PKVujTGG9957j5dffhmAV155halTp5InTx6t8q3SGg/0N8asMcacd2xrgIHAaGPMT1iTXULtDFLZ47OF6+jZrikXTkThXaAI7375HVWrVrU7LJXNMlQR3BjzozEm1BhT1BgTYIypZYwJz+zglLoVKSkpDBgwgHfeeQeAd999l3HjxuHlZf211yrfKo3yQHrracYCdzge/w0UzraIlFvYtm0b/Z5tQfL5U+QJKk3xjqP54YTWhs6N9E9d5UiXL1+mR48efPPNNwCMGzeOPn36XHWOVvlWaWwFRovIs8aYUwAiUgT4EOu2HVh1mw7ZFJ+ywbp162jWrBmJF2LxLX4nRdsMxTugoH7JyqVcSppE5AzWYMgbMsboXEtlq4SEBDp06MD333+Pl5cXU6ZMoXPnzumeq1W+lZNuwGLgsIikJka3A/uAFo7n+YDhNsSmsllUdBzvfDqdb4f3JTkpgTp16tBz+Jcs3HVGv2TlYq72NPXNyiCUulnXqvJ94cIFWrVqxerVq/H19WXu3Lm0bNnS3mCVRzDG/CkilYAGwF2O3X8CEcaYFMc5i2wKT2WjqOg4anQbwvGlYyElmbzlH2LoZzN59J7StK9td3TKTi4lTcaYaVkdiFKuulaV77tuE5544gk2bdpEYGAgixcvpn79+jZHqzyJIzla6dhULjV8zHiOLx4DGAIq1aXwE/1YsTuGR+/RZVFyuwyNaRIRb6AV1qw5sGo3LXZU01UqS6VX5XvG/yJZN74vO3fu5LbbbmPFihXUqFHDpgiVpxKRQKAuEAL4Oh8zxkywJSiVbYwxjBw5kskj3wAg3/1NCAp9AREvWlXV23EqYwv23gMsAYpjdV0DvA6cEpFmxpidmRifUv+Stsp30tnjLB42nCOHDlCiRAnCw8OpXLmyzVEqTyMi9wPLgQAgEIjBmikXB5wENGnKwYwxvP7664wePRqA3v1fo2T9zhw+E3/VEACVu2Wkp+kr4HfgwdSquCJyG/ANMAn4b6ZFp1Q6nKt87/h1J6unvsnpkye44447iIiI4Ldzfgz7coMO1lQ362NgKdATOAfUBJKw1qIbb2NcKoslJyfTs2dPvvrqKwDGjBnDgAEDbI5KuaOMJE1VcUqYAIwxZ0TkTf6ZlqtUlgoJDuCxoFhGf9iTmJgYKleuTHh4OJuPp9BnbiQAm/fHAGjipFxVFXjBGJMiIsmAnzFmn4i8BkwDvrc1OpUlfvzjGD26deHPDeF4eXkxadIkunXrZndYyk1lpLjlX0CxdPYXJf2FfJXKdGvWrKF+/frExMRQo0YN1q5dS4kSJbTKt7oVSUCK4/FJrHFNYPU63W5LRCpL/e+3KBo0acqfG8LBKw/vjf9KEyZ1XRlJmgYDE0SktYiUdmytgXHA6yJSIHW70RuJSB0RWSoiR0XEiEjLNMdFRIaJyDERuSQiq0WkQppzgkRkpojEishZEflaRPJl4PdSHmLRokU88cQTXLx4kfr167N69WqCgqzyYFrlW92C7cBDjsdrgWEi0hGrbdOxmjnM2bNn6fBUc+L3b0N8/Cja+h2O3FbF7rCUm8vI7bkwx895/FPwMnVV8KVOzw3gfYP3CgR2AFNIv+v7NeAVoDOwH3gPWCUilYwx8Y5zZgIlsNaD8gGmYo2t6uD6r6Q8xfTp0+natSvJycm0atWK2bNn4+fnd+W4VvlWt+ANIL/j8ZvAdOBzrKVTutoVlMp8J06coGHDhhz/ewdefoEUaT0E/9IVuT0or92hKTeXkaTpscz6cGPMCmAFgIhcdUysHX2B4Y5FghGRZ4ETQEtgjohUBBoBDxljfnGc8zKwXEReNcYczaxYlf0mTJhwZSmULl26MHnyZPLk+fdfYa3yrW6Wo705iaNHyRhzEqttUTnMwYMHefzxx9mzZw9FihYlb/N3kOCy+HoLHaqXsTs85eZuOmkyxqzNikDSUQ6rrMFqp88+JyKbgIeBOY6fZ1MTJofVWOMSagALsylWlQmiouNYsfMYAI0rlyAkOACwpgIPGzaMIUOGANC3b18++uijKwvvKpUJBGtM5j1YPUsqB5q44Ade6/40l86eokyZMtbKAbeVZPOBGKqXDbrS5ih1LRktbukP3Ic1+Puq/7mMMUsyIS6wEiawepacnXA6Vhzr26Hz518WkRinc/5FRPwAP6dd+a91rsoeaat8fxzxF+H96lL6Nn/69evHhAlWiZxhw4bx1ltv/atnUqlb4Zgx9zcQjItJk7YjnuWjmSt4rUc7Ui7F4hN8O+98+R133nkngCZLymUZKW7ZCOtef+F0DrsyjskdDAbetTsI9Y8VO49dVeU7/nIKG/acZMXnQ5g2zVrFZ8KECbz88st2hahyvkHAaBF50cUivdqOeIi1a9cyqHtrUuLj8C1egaJthrDyQJIOVFM3LSP3Nz4B5gMljDFeabbMTJiOO36mLW9QzOnYcazeritEJA8Q5HROej4ACjptuqCQm/E2SUwd9grTpk3D29ub6dOna8Kkstp0oDqwwzFbN8Z5S+d8bUc8QFhYGI0aNeJyfBx+IfdSrN37eAcU1Jm1KkMycnuuGDDWGJP2tllm24+V+NQHIgEcZQxqYM1oAdgAFBKRasaYrY599bCSwU3XemNjTAKQkPpcb/XYr3HlEoyN+IuEyyl4JV3itg0TWLXhR/z8/Jg3bx7Nmze3O0SV8/W9mZO1HXF/M2fOpHPnziQnJ9O8eXM6DR7L/B0ndWatyrCMJE3fAY8Ce2/1wx31lO502lVORKoCMcaYKBEZB7zlGGuQWnLgKLAIwBjzh4isBCaLSE+skgMTgTk6c86zhAQHENGvLmsi9zLhta5s3fYL+fLlY8mSJTz22GMs3n5EywioLGWMmWZ3DCrzfPbZZ/Tu3RtjDJ06dWLKlCn4+PjQpmZ5u0NTHiwjSVNvYL6I1AZ+w6qie8VNrgT+IPA/p+djHT+nAV2AD7FqOU0CCgHrgUZONZoAOmIlSmuwZs0twKrtpDxMnoSzfNSnPb///jtBQUGsXLmShx56iMXbj+jSKCpbiEh54DmgPNDHGHNSRBoDUcaY3+2NTrni4OmLPD/gDcKnW/8V9e7dm/Hjx+tsW5UpMpI0tQcaAPFYPU7G6ZjhJlYCN8b8wD+FMdM7boB3HNu1zolBC1l6vL179xIaGsr+/fspWbIkERERVKpUCfj3UigzN0dp0qQynYjUxaob9xNQB6vA5UmgCtANaG1fdMoVB09fpGrTZzi7yao2U/C/7enQ511NmFSmycjfpPexZowUNMaUNcaUc9ruyOT4VC7w22+/UatWLfbv30/58uVZv379lYQJdGkUlW1GAm8ZY0KBRKf9/wfUtCck5ark5GSe69btSsJ0W70eFKrdkUU7dKSGyjwZ6WnyBeYaY1JueKZSN7Bx40aeeOIJzpw5w7333suqVasoUaLEVefo0igqm9xL+r3WJ0m/xIpyEwkJCXTs2JH/LVkA4kVwo5fJd18oAK2qanuhMk9GkqZpQFtgRCbHonKZ1atX07JlSy5evMjDDz/MsmXLuO2229I9V5dGUdngLNY6lvvT7L8fOJLt0SiXzP35b/p078iJP7bg6+vLhEnfcKZoVQ7FXOKZmmWoWT7Y7hBVDpKRpMkbeE1EGgK/8u+B4P0zIzCVM1xr1tv3339P+/btSUxMJDQ0lIULFxIYGGhjpEoxBxglIm2wxmd6icgjwBisGk7Kzcxc+zvdOrQm4ehuxMefN8dP5YXO7ewOS+VgGUma7gW2Ox5XTnPMoJTDtWa9TZ06le7du5OSksJTTz3FzJkz8fPzu847KZUt3gA+BQ5hfTnc5fg5CxhuY1wqHcePH6dXxxYkHN2Ll38+irYewk7RBXdV1srIgr2PZUUgKudJb9bbvh/m0b+/1RnZtWtXvvzyS/LkydASiEplKmNMItBDRN7D+kKYD9hujNEFfN3M/v37CQ0N5dyRvXgH3kbRtu/hW6SsThJRWU7/t1JZpmP1kCs9TMYYZOtc+k8eD8CAAQMYPXq0VlFWbkNEahlj1htjooCoG75A2WLXrl2EhoZy9OhRypUrx2vjZvB/x7x0kojKFi4lTSLyPdDFGBPreHxNxpgnMyUy5fFSG7BvNx4gbu1XzJv7DQDDhw/njTfeYNO+GBZGHqFV1VI6WFO5g/8TkSPAbOBbY8wuuwNSV9uyZQuNGjUiJiaGSpUqERERQcmSJelpd2Aq13C1TtM5/hmvdO4Gm1JXPFG5KP4/f8FyR8L06aef8uabb7Ik8ijtJm9k7pZDtJu8kY17o+0NVCkoCXwE1AV2ikikiAwUEV2I12ZR0XH0H/cttes+RkxMDNWrV2fdunWULFnS7tBULuNST5Mx5rn0Hit1PfHx8Tz99NMsXboUb29vpk2bRseOHa0GcF7kVecujDyivU3KVsaY01hLMk0UkXJYNZs6Ax+IyDpjTD1bA8yloqLjqNlzJMcWjoTkJALKVmHq3CUEB2t7obLfTVcEF5G8IhLg9LyMiPQVkQaZG5ryZOfPn6dx48YsXboUf39/Fi1aRMeOHQHYfCCG5DTzLLUAnXInxpj9WBXCB2GtsVnX3ohyr3fHfsaxBe9DchJ5K9Sk8FPvsis66cYvVCoLZGQZlcXAswAiUgjYDAwAFovIi5kXmvJUp0+fpl69evzwww/kz5+flStX0rRp0yvHq5cNIq+PNwDeAuPbVtVeJuU2ROQREfkMOIZVbmAn0MTeqHKnTz75hG9GDASTQmDlehRpORhfPz+qlw2yOzSVS2Vk9twDQD/H49bAcayKuU8Bw4DPMyc05YmOHDlCaGgof/zxB8HBwaxatYpq1apddU5IcACr+tZh84EYqpcNIiQ44BrvplT2EZEPgHZYY5sigD7AYmNMnK2B5ULGGIYPH84771hrtRd6qDkFHuuOj7c3M7rW0DZD2SYjSVMAcN7xuAHwvTEmRUQ2AlpZLBeIio5LN+HZs2cPjz/+OAcPHqR06dKEh4dTsWLFdN8jJDhAGz7lbuoAo4F5jvFNygYHTl2g0/O9+GmRVYR9yJAhdOn1KlsOntEvWcp2GUma9gAtRWQh0BD42LG/KBCbWYEp97R4+xH6z99Bcoohr483q/rWISQ4gF9//ZUGDRpw4sQJ7rzzTlavXk2ZMppDK89hjHnE7hhyu30nYqna4CnO/7oagNvq96Bhp96UKRxImcK6zJKyX0aSpmFY9/k/BtYYYzY49jfgn+VVVA60cW/0lWVRAC4lJbP5QAyH/4ykSZMmnD17lipVqrBq1SqKFStmX6BK3QIRqQSEAL7O+40xS+yJKHdISEigffu2VsIkXgQ/0Yd8levrzFrlVjKyjMp3IrIeazXwHU6H1gALMysw5X4WRl690Lu3lxC/fxuhndsTFxfHI488QlhYGIUKFbInQKVugYjcgdWG3YtVly61XH3qXE9vO+LKDS5cuEDLli3Z/L81iHceCjd/nYC7HgZ0Zq1yLxlaRsUYcxxrALjzvs2ZEpFyW62qlmLulkNXnre+7RBdO7xMUlISDRs2ZMGCBQQGahe68ljjgf1AfcfP6kAwVsHLV22MK0eLiYnhiSeeYNOmTQQGBjJ5xjwO+pfjUMwlnqlZRnuZlFvRteeUy2qWD2ZOj5osjDyC91//Y9Sg/qSkpNCmTRu+/fZbVvx+ipmbf9U1oJSnehioZ4w5LSIpQIoxZr2IDAYmYM0SVpno2LFjNGjQgJ07dxIUFMSKFSuoXr263WEpdU0ZqdOkcrGa5YMpeiCCD97oS0pKCt27d2f27Nms+P0UfeZGsnl/DH3mRrJ4+5Ebv5lS7sWbf2YGn8YqPQBwELjblohysC+X/kyF+x5k586dlChRgnXr1mnCpNyeJk3KZcYY3nzzTV591bpTMXDgQCZNmoS3tzczN1+9KHza50p5gJ1AFcfjTcBrIvII8A6wz7aocqAJ89fQq0MzLp4+Sp5CxRkyaQH33HOP3WEpdUOaNCmXpKSk0Lt3b0aMGAHABx98wIcffoiINVa2Y/WQq85P+1wpDzCcf9rEd4BywI/AE8ArdgWV02zatIlXu7Qi+UIMPoXLUKzjh6zWjmnlIXRMk7qhpKQkunTpwqxZsxARPvvsM3r27HnVOaljmGZujtIxTcojGWNWOT3eA/xHRIKAM8YYc+1XKletWbOGFi1akBR3Ed+Sd1O09RC88+bXL1nKY2jSpK7r0qVLtGnThmXLlpEnTx6mT59O+/bt0z23xf2lNFlSOYoxJsbuGHKCjXujGfXldJaNG0RSUiKPP/443YZ8yvc7o/VLlvIomjQpNu6NZmHkEVpVLXXV9N7Y2FiaNWvGunXr8Pf357vvvqNJE123VCnluo17o2n6yvtErxgPJoVHGzYlbPF3+Pn50U5rsCsPo0lTLrdxbzTtJm8EYO6WQ8zpUZOa5YM5deoUjRo1Ytu2bRQoUICwsDBq165tc7RKKU/z6pAPiF5urbYVeO/j/KfD2/j5+dkclVIZo0lTLjdj48F/PS/lG0eDBg3YvXs3RYoUYeXKlTzwwAM2RaiU8kTGGIYOHcpP334EQP4HW3BbvW6UKZLf5siUyjhNmnK524PyXvXcP+4EtWq1Iioqittvv52IiAjuvltL1CilXJeSkkK/fv2YMGECAMF1niGw5tP45fGiQ3VdyFt5Lk2acrkO1cswZf1+EpMN5vQ+pk95j9OnTnHXXXcRERFBSIjOalFKue7y5ct069aN6dOnA/DJJ5/QvH1XNh+IoXrZIEKCA2yOUKmM0zpNuVxIcACr+z9K9wqJnJ3/NqdPnaJq1ar8+OOPhISEsHj7EZ7+coNW+FZK3dBfR2L47+NNmD59Ot7e3syYMYPevXsTEhxA62qlNWFSHk97mhS7tqxjxMsduXTpErVq1SIsLIyCBQuyePsR+syNBGDzfmvmtU4NVkqlZ876P3muU1viD+5AvH348puZdOrUxu6wlMpU2tOUy82bN4/mzZtz6dIlGjduzKpVqyhYsCDw76VQdGkUpVR6Vv3yN52fbmYlTL55KdpmKAUrPmx3WEplOk2acrHJkyfTrl07kpKSaNu2LYsWLSIg4J/uc10aRSl1I0eOHKFjq0YkHvsbr7wFKNZuBIHlqlC9bJDdoSmV6dw6aRKRISJi0my7nY77i8inIhItIhdEZIGIFLMzZk/x4Ycf8vzzz2OM4YUXXmDmzJn4+vpedU6L+0sxvm1VqpcLYnzbqnprTil1lb1791KrVi2iD+/DO18wxTqMxK9EBca2qaLjl1SO5Aljmn4HHnd6ftnp8cdAE6ANcA6YCHwPaJ3ZazDG8MYbbzBy5EgABg0axIgRI64svJuWLo2ilErP+HmrGdSjHfGx0ZQvX57RX89n42nvf60soFRO4glJ02VjzPG0O0WkINAN6GCM+T/HvueAP0SkpjFmYzbH6VYWbz/yr8Vzk5OT6d27N1988QUAo0aN4rXXXrMzTKWUB/pw2hIG9+xISvwFfIqU5a3P59Oq7v20sjswpbKYJyRNFUTkKBAPbAAGG2OigGqAD7A69URjzG4RiQIeBnJt0pTerLfG9xShc+fOzJkzBxHhiy++4Pnnn7cxSqWUJ4qIiOCN59uSkhiPX8n/UKTNEJbvi6eL3YEplQ3cPWnaBHQB/gRKAO8CP4pIZaA4kGiMOZvmNSccx65JRPwA58WPclRd/7Sz3Kat/4sv3+zBihUr8PHxYcaMGbRt29am6JTKGXJ6O5Ke77//nvbt25OcmIh/2fsp0upNvHz9dZKIyjXceiC4MWaFMWa+MeZXY8wq4AmgEPD0Lb71YKwxUKnb4Vt8P7fi3IClJFzkz6mDWbFiBXnz5mXJkiWaMCmVOXJ0O5LW1KlTadOmDYmJibRu3ZoZc76j5t0ldZKIylXcOmlKy9Gr9BdwJ3Ac8BWRQmlOK+Y4dj0fAAWdttKZGqjNUme93RcMPivfY9f2TRQsWJDw8HASit2rFb6Vyhw5uh1x9vHHH9O1a1dSUlLo1q0bc+bMoXWNO5j3wsOaMKlcxd1vz11FRPIB5YEZwFYgCagPLHAcvxsIwRr7dE3GmAQgwel9syhi+9wfnMyfX/Vn319/UaRIEcLDwzloimiFb6UySW5oRzbsOc3Lrw5i6+KvAXj11Vf58MMPc+TvqpQr3LqnSUTGiEhdESkrIv8FFgLJwGxjzDnga2CsiDwmItWAqcCG3D5z7s8//6RWrVr89ddfhISEsH79eqpWraoVvpVSLvv571M0bPvclYQp+NHO9H59iCZMKldz66QJq7t7NtZA8HlANFDTGHPKcbwfEIbV07QO67bckzbE6Ta2bdtG7dq1OXToEHfffTfr16/nrrvuArTCt1LKNUlJSfTo/hznt4UBQlCDl8hXow1bDp6xOzSlbOXWt+eMMe1ucDwe6OXYcr0ff/yRpk2bEhsbywMPPMDKlSspUqTIleOpt+LS1m9SSqlUly5dom3btuxatwy8vCncpB+BlR7FL4+XLo2icj23TpqU65YvX85TTz1FfHw8derUYcmSJVcW3nWmFb6VUtcy+8fd9OnWgVN/b8ff35/hn0xhX967CQkKoH31EF0aReV6mjR5kPSqfAPMmTOHZ555hsuXL9OkSRPmz59P3rx5bYxUKeVpZvzfb3Tv+BSJx/9GfPPy1ifTGdC9td1hKeVWNGnyEOlV+W5xfym++OILXnrpJYwxdOjQgW+++QYfHx8bI1VKeZrDhw/Tq2NzEo8fwCtvAYo+PYwdydojrVRamjR5iPRmvv2xagaDBw8G4MUXX2TixIl4ebn72H6llDvZs2cPjz/+OOePH8Q7f2GKtX0Pn+DbdaKIUunQ/2E9hHMDZozh8oZ/EqY33niDTz/9VBMmpdRN2bFjB7Vq1eLgwYNUqFCBz2eH8ciDVbTKt1LXoD1NHiK1Aft2437ORXzOwoWzABg9ejQV6ren7aSNOiNOKeWyn3/+mSZNmnD27FmqVKnCqlWrKFasGD3sDkwpN6ZJkwdpfE8RZo3sT8TCeXh5eTFp0iQKP9BIq3wrpVwWFR3HkM9nMmP4K1xOiOeRRx4hLCyMQoUK2R2aUm5P7+d4iLi4OFq0aMG8efPw8fFh7ty5dOvWTat8K6VcFhUdR40ew5n67ktcTojHv1w1hn8xWxMmpVykSZMHOHv2LA0aNGDlypUEBASwdOlSWre2pgJrlW+llKs+GPcZxxeNgpTLBPynNkWfeosVf2qVb6Vcpbfn3NyJEydo2LAhO3bsoFChQixbtoz//ve/V45rlW+llCvGjBnDF8MHApCvSkOCGryEeHnTqqq2GUq5SpMmN3bw4EFCQ0P5+++/KVasGOHh4dx3333/Ok+rfCulrmXRtsMMfvNNdq+cDkDPl/sR0qg7h8/E80zNMtQsH2xzhEp5Dk2a3EB6lb53795NaGgohw8fpkyZMkRERFChQgWbI1VKeZKFWw/xbI+eXNi+HIBnXh7E5xM+sDkqpTyXJk02S6/Sd+mU4zRq1IjTp09TsWJFwsPDKV26tI1RKqU8TVJSEn1f7M6F7eGAENTwJeLvaW53WEp5NB0IbrO0s93GfbuYxx57jNOnT1OtWjXWrVunCZNS6qZcunSJVq1aEbUlHLy8KdzsVfJXbawTRZS6RdrTZLOO1UOu9DDF7dnMz2GjSExIoG7duixZsoQCBQrYHKFSypOcO3eO5s2bs27dOvz9/Xn1wy/507eCThRRKhNo0mSz1Ebsg4lf8cviESRfvkyzZs2YO3cu4btjmLn5d23slFIuWbH5T55r34oT+/6gQIEChIWFUbt2bbvDUirH0KTJDRzZsJjNU4dijKFTp05MmTKF5TtPaqVvpZTLvlqxhRc7PcnlmMN4BRRkwoyFmjAplcl0TJONjDG8//779OrVC2MMvXr1Ytq0afj4+Gilb6WUy37YtIMX2zXhcsxhvPMXoXiHUexK0lICSmU2TZpsYoxh4MCBvPXWWwC8/fbbfPLJJ3h5WX8kWulbKeWKyMhIWj7xOJdjT5EnqBTFO43CJ7i0Fq1UKgvo7bksFhUdx+YDMVQvG0RIcAAAycnJvPDCC3z99dcAjB07ln79+l31Oq30rZS6nqjoOEZM/Z6p775IYtwF/IqXp0jrofjmK8TYp6tq0UqlsoAmTVkoKjqO0I/XknA5Bb88XkT0q0uxfN507NiRBQsW4OXlxeTJk+natWu6r9dK30qp9ERFx/FwrzEcWzACczkBv9L3MGXmfHwD81/1BU0plbk0acpCszdHkXA5BYCEyylMW7eb9V8MJjw8HF9fX2bPns2TTz5pc5RKKU/z8aRvODr/PUi5TN47HqRwy0HsOJXEqDpa002prKRJUxaKiom78jg5/gKfvv4cJ/7+lcDAQBYuXEhoaKiN0SmlPNHkyZMZ/2ZvMIaAinUo3KQf4u2jY5iUygaaNGWhZ2qWYdlvx0i+cIYT894m6dQBChUqxOBx05i8Lx9x24/o7TellEuiouN4c+j7fPvJCAA6delOxaf6cORcoi68q1Q20aQpC9UsH8zHjUvStd2LJJ06RPHixRk0fgYfb0sAYrT+klLKJRv2nOaJZ1/i7Ib5APTq+yqfjP0QEbE5MqVyFy05kIV27dpF7w5NOXv8EOXKlWP9+vX8dCbwqnO0/pJS6nr2nzxPo6efvZIwFXr0Oe5t2VMTJqVsoD1NWWTLli00btyY6OhoKlWqRHh4OKVKlaJjrP+VHibQ+ktKqfRt3BvNd1sOEPH528RuXwEIQQ17kb9qI7tDUyrX0qQpE2zcG83CyCO0qlqKmuWD+eGHH2jWrBkXLlzgoYceYsWKFQQHW+MNtP6SUupGNu6N5unPfuDUog+I37cVvPJQuOkAAivWxtdbaFy5hN0hKpUradJ0izbujabd5I0AzN1yiF7lzvB2n+4kJCRQr149Fi1aRP78+a96jdZfUkpdz1f/t5OT894h4fAuJI8f9V8eResWTQFoXLmE1mFSyiaaNN2iGRsPXnl84ff/8frocZiUZFq0aMGcOXPw9/e3MTqllKc5efIky0f2JOHwbsQvkKKt36V+aANeqFve7tCUyvV0IPgtuj0oLwCxW5cSHfYRJiWZx5q2Jk+DAaz6I9rm6JRSniIqOo5B09Zwd5XqHNu3G++AQhRv/wEFytxDh+pl7A5PKYX2NN2y9g+FMHrkCM6s+xaAek89y57yrdkXFcsvUZGAlhRQSl1fVHQcdd78lsOz3iL5/Cm8CxRhwrQFFL29nC6LopQb0Z6mDIiKjuO7rYc5cOoC40e8TYwjYeo78A2CHn8BkX8uq5YUUEpdz+LtR3jyvW85NH0gyedPkSeoNMU7juZgckFaVyutCZNSbiTH9DSJSC9gIFAc2AG8bIzZnNmfk7oIb3xiEmdXfULsr6sBGDduHH369GHx9iNsOXDmyvlaUkApdS2Ltx/hhTHfcvK7YZjEOHyL30nRNkPxDiioy6Io5YZyRNIkIm2BsUBPYBPQF1glIncbY05m5mfN3hxFfHwCp5eOJu6vn0G8eGXIR/Tp0wfQkgJKKdcs3n6EV0Z/zcn572MuJ+J3e2UqdxlOjbtv12VRlHJTOSJpAvoDk40xUwFEpCfQBOgKjMzMD1qyZS8nvxtK/MFI8M5DkeavsziuAvWc1pHTkgJKqetpOmEdm1aHcXrZWEhJJm/5hyjcYhBvt3pQ2w6l3JjHj2kSEV+gGrA6dZ8xJsXx/OHM/KxXZ6znly8GEH8wEvHxp2jrIQTcZX2Ejl1SSrmix7QtbFg2l9NLx0BKMgGV6lKm7Tt80qmGJkxKubmc0NNUGPAGTqTZfwL4T3ovEBE/wM9pV/70znN27NgxJvTvRNLpg3j556Nom6H4lbz7ynEdu6RU7pKRdgTgu6mfcnbtNADy3d+EoNAXeLJaGU2YlPIAHt/TlEGDgXNO2+Hrnbxv3z5q1apF0umDeOcLoliHkVcSpurlghjftqo2eErlPjfVjhhjeP31168kTAUebktQaE9EvLRwpVIeIif0NJ0GkoFiafYXA45f4zUfYA0cT5WfazR4v//+O6GhoRw7dow8hYpTtO1wfAoVv3J83guZegdQKeU5XG5HkpOTefHFF5k8eTIAtz3WlQLVn7xyXMsKKOUZPL6nyRiTCGwF6qfuE6tQUn1gwzVek2CMiU3dgPPpnbdp0ybq1KnDsWPHqFy5MlG7tl2VMB0Y2SQzfxWllAdxtR1JTEykffv2TJ48GS8vL77++uurEiZtR5TyHDmhpwmsb3vTROQXYDNWyYFAYGpG33DNmjW0aNGCixcvUqNGDZYvX05QUJA2cEopl128eJGnnnqKVatW4ePjw6xZs2jdujVd7Q5MKZUhOSJpMsbMFZEiwDCs4paRQCNjTNrB4S5ZtGgRbdu2JTExkfr167No0SLy5cuXiRErpXK6s2fP0rRpU3766ScCAgJYuHAhDRo0sDsspdQtyBFJE4AxZiIw8VbfZ9q0aXTt2pWUlBRatWrF7Nmz8fPzu/ELlVLK4cSJEzRs2JAdO3ZQqFAhli9fzsMP6/hHpTydx49pykyff/45Xbp0ISUlhS5dujBv3jxNmJRSNyUqKopatWqxY8cOihUrxtq1azVhUiqHEGOM3THYTkQKYE0ZBqBv37589NFHeHlpTqmUu4uNjaVgwYIABR0Dsm2R2o6UKFGCY8eOUbZsWSIiIrjzzjvtCkkp5SJX2xHNCtIYNmwYY8eO1YRJKZUhx44do2LFiqxfv14TJqVymBwzpikzDB06lD59+nD+fLozh5VSbig21rbOpXTdd999LFq0iPz587tdbEqp9Ln6b1VvzwEiUoobVPNVSrm90saYI3Z9uLYjSuUI121HNGkCRESAklyjOF0aqVV/S7t4vsocet3t4SnXPT9w1NjYoGk74hH0utvDU677DdsRvT0HOC6QS99QrXYRgPN2DjrNbfS628ODrrvtsWk74v70utvDg677DWPT0c5KKaWUUi7QpEkppZRSygWaNN28BGCo46fKPnrd7aHXPWvodbWHXnd75JjrrgPBlVJKKaVcoD1NSimllFIu0KRJKaWUUsoFmjQppZRSSrlAk6abICK9ROSAiMSLyCYRqW53TJ5MROqIyFIROSoiRkRapjkuIjJMRI6JyCURWS0iFdKcEyQiM0UkVkTOisjXIpIvW38RDyMig0Vki4icF5GTIrJIRO5Oc46/iHwqItEickFEFohIsTTnhIjIMhGJc7zPaBHR2m8u0LYkc2lbkv1yazuiSZOLRKQtMBZrBsADwA5glYgUtTUwzxaIdR17XeP4a8ArQE+gBnAR65r7O50zE7gHCAWaAnWASVkVcA5RF/gUqIl13XyAcBEJdDrnY6AZ0MZxfkng+9SDIuINLAN8gf8CnYEuwLCsD9+zaVuSJbQtyX65sx0xxujmwgZsAiY6PffCqv47yO7YcsIGGKCl03MBjgGvOu0rCMQD7RzPKzpe96DTOY2AFKCk3b+Tp2xAEcd1rON0nROB1k7n/MdxTk3H88ZAMlDM6ZyewDnA1+7fyZ03bUuy/PpqW2LPdc8V7Yj2NLlARHyBasDq1H3GmBTH84ftiiuHKwcU5+prfg7rP5zUa/4wcNYY84vT61ZjNXQ1sinOnKCg42eM42c1rG+Nztd+NxDF1df+N2PMCaf3WQUUwPq2rtKhbYkttC3JHrmiHdGkyTWFAW/gRJr9J7D+MarMl3pdr3fNiwMnnQ8aYy5j/aPVPxcXiIgXMA74yRiz07G7OJBojDmb5vS01z69PxvQa3892pZkP21LslhuakfcdrCVUipbfApUBmrZHYhSymPlmnZEe5pccxrHfdc0+4sBx7M/nFwh9bpe75ofB64aPOuYdRGE/rnckIhMxBrw+pgx5rDToeOAr4gUSvOStNc+vT8b0Gt/PdqWZD9tS7JQbmtHNGlygTEmEdgK1E/d5+iOrA9ssCuuHG4/1j8a52teAGt8Qeo13wAUEpFqTq+rh/X3elM2xelxHNOvJwKtgHrGmP1pTtkKJHH1tb8bCOHqa39vmhlfoUAssCurYvd02pbYQtuSLJBr2xG7R6J7yga0xZpt0RlrpsWXwBmcRv3rdtPXNB9Q1bEZoJ/jcYjj+OuOa9wcuBdYBOwD/J3eYwWwDagOPAL8Bcyy+3dz5w34DDiLNQW4uNOW1+mcz4GDwGNYAzp/Bn52Ou4N/IY1aLMK0BBrTMgIu38/d9+0LcmSa6ptSfZf81zZjtgegCdtQG/HX4AErG8fNeyOyZM34FFHA5d2+8ZxXLDqdRx3/CezGrgrzXsEAbOA81jTVKcA+ez+3dx5u8Y1N0AXp3P8scYpxGDVtPkeKJ7mfcoAy4E44BQwBshj9+/nCZu2JZl+PbUtyf5rnivbEXEErZRSSimlrkPHNCmllFJKuUCTJqWUUkopF2jSpJRSSinlAk2alFJKKaVcoEmTUkoppZQLNGlSSimllHKBJk1KKaWUUi7QpEkppZRSygWaNKkcQUS+EZFFdsehlPJc2o6oG9GkSWUbERkiIpF2x6GU8lzajig7adKkVBYSEV+7Y1BKeTZtR9yHJk3KZSLSSETWi8hZEYkWkTARKZ/mnNIiMltEYkTkooj8IiI1RKQL8C5QRUSMY+siImUdj6s6vUchx75HHc+9ReRrEdkvIpdE5E8R6ZOB+GuJyI+O9zgkIhNEJNDp+AEReUNEpojIeRGJEpHn07zH7SIyz3ENYkRksYiUdTr+jYgsEpE3ReQo8Kdj/39FJFJE4h3XpGXq7y2WPSLyaprPquo4586b/V2Vclfajmg74sk0aVI3IxAYCzwI1AdSgIUi4gUgIvmAtUApoDlQBfgQ6+/ZXOAj4HeghGOb6+LnegGHgTZAJazVykeIyNOuBu5olFcCC4D7gLZALWBimlMHAL8A9wOfAZ+LyN2O9/ABVmGtgl4beAS4AKyUq78J1gfuBkKBpiJSAFgK/AY8ALwNjEo92VirZk8BnksTy3PAOmPMHld/T6U8gLYj2o54LmOMbrplaAMKAwao7Hj+PBALBF3j/CFAZJp9ZR3vUdVpXyHHvkev89kTge+cnn8DLLrO+V8BX6bZVwtIBvwdzw8AM5yOC3AC6Ol43gnYDYjTOb5AHNDAKY7jgK/TOT2B06mf49jX3fn3BkoCl4Hqjuc+wCmgs91/zrrplpWbtiPajnjSpj1NymUiUsHRZb5PRGKxGgeAEMfPqsB2Y0xMFnx2LxHZKiKnROQCVsMacqPXOakCdBGRC6kb1rc9L6Cc03m/pj4wVqtzHCjq9B53Aued3iMG8Aecby/8ZoxJdHp+N/CrMSbead9m5+CMMUeBZUBXx65mgB8w/yZ+R6XcnrYj2o54sjx2B6A8ylLgINADOIrVUOzE+pYEcCkD75ni+ClO+3ycTxCRdsAYrC7vDVjd2gOBGjfxOfmAL4EJ6RyLcnqclOaY4Z/b2PmArUDHdN7jlNPjizcRl7OvgBki0g+rS32uMSYug++llLvSdkTbEY+lSZNyiYgEY33T6WGM+dGxr1aa034FuotI0DW+JSYC3mn2pTYSJYDtjsdV05zzCPCzMeYzp3jKc3O2AZXMrd3X34Y1huGkMSb2Jl73J9BJRPyMMQmOfQ+lc95yrIbyRaARUOcWYlXK7Wg7cuU9tB3xUHp7TrnqDBANPC8id4pIPazBnM5mY3VDLxKRR0TkDhF5SkQedhw/AJRzzOYo7PjHfwnYCAwSkYoiUhcYnuZ9/wYeFJGGInKXiLxH+o3F9YwC/isiEx2fX0FEWohI2gGc1zMTa0zBYhGpLSLlRORRx+yZ0td53Sysf2uTHL9jQyB1hotJPckYk4w1luED4G9jzIabiE0pT6DtiLYjHk2TJuUSY0wK0A6ohtWV/jFW17bzOYlAA+Ak1red34BBWIMkwZpxshL4H9Y3w/aO/V2xej23AuOAt9J8/JfA91izZDYBwVgzUm4m/l+BusBdwI9Y30aHYd0ecPU94rC+tUU54vkD+BprLMI1vzE6vk02w/rmGwm87/hsgPg0p3+NdZtiqqtxKeUptB3RdsTTiWOEvVIqG4lIR6wGraDjW3Lq/trAGuB2Y8wJu+JTSrk/bUeyn45pUiobiMizwD7gCNbsmVHAvNSGTkT8gCJY06nna0OnlEpL2xH76e05pbJHceBbrK74j7GmADtXCW6PNaOoEPBadgenlPII2o7YTG/PKaWUUkq5QHualFJKKaVcoEmTUkoppZQLNGlSSimllHKBJk1KKaWUUi7QpEkppZRSygWaNCmllFJKuUCTJqWUUkopF2jSpJRSSinlAk2alFJKKaVc8P+T/H56VHP/YwAAAABJRU5ErkJggg==\n", | |
"text/plain": [ | |
"<Figure size 600x600 with 4 Axes>" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"fig, axes = plt.subplots(len(data_cases), len(calculation_cases),\n", | |
" sharex=True, sharey=True, figsize=(6, 6), dpi=100)\n", | |
"\n", | |
"for i, (case_name, cumulative_energy) in enumerate(data_cases.items()):\n", | |
" for j, (calc_name, function) in enumerate(calculation_cases.items()):\n", | |
" calculated_energy = function(cumulative_energy)\n", | |
" axes[i, j].scatter(interval_energy, calculated_energy, s=3)\n", | |
" axes[i, j].axline((0, 0), slope=1, c='k')\n", | |
" axes[i, j].set_title(case_name)\n", | |
" axes[i, j].set_ylabel(calc_name)\n", | |
" axes[i, j].set_xlabel('actual energy')\n", | |
"\n", | |
"fig.tight_layout()" | |
] | |
} | |
], | |
"metadata": { | |
"kernelspec": { | |
"display_name": "Python 3", | |
"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.7" | |
} | |
}, | |
"nbformat": 4, | |
"nbformat_minor": 4 | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment