Created
March 12, 2026 18:18
-
-
Save adrn/b815f9cd01394aa48478d2adcc284ac5 to your computer and use it in GitHub Desktop.
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": 33, | |
| "id": "07cbe3b0", | |
| "metadata": {}, | |
| "outputs": [], | |
| "source": [ | |
| "import arviz as az\n", | |
| "import astropy.table as at\n", | |
| "import astropy.units as u\n", | |
| "import corner\n", | |
| "import matplotlib as mpl\n", | |
| "import matplotlib.pyplot as plt\n", | |
| "import numpy as np\n", | |
| "import pymc as pm\n", | |
| "import thejoker as tj\n", | |
| "import thejoker.units as xu\n", | |
| "from astropy.time import Time\n", | |
| "\n", | |
| "%matplotlib inline" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": null, | |
| "id": "f76adbfe", | |
| "metadata": {}, | |
| "outputs": [], | |
| "source": [ | |
| "rng = np.random.default_rng(42)\n", | |
| "fake_t = 55555.0 + np.sort(rng.uniform(0, 1000, size=3))\n", | |
| "\n", | |
| "fake_rv = np.zeros(len(fake_t)) * u.km / u.s\n", | |
| "fake_rv_err = np.full(len(fake_t), 0.5) * u.km / u.s\n", | |
| "fake_rv = rng.normal(fake_rv.value, fake_rv_err.value) * fake_rv.unit\n", | |
| "\n", | |
| "data = tj.RVData(fake_t, fake_rv, fake_rv_err)" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 78, | |
| "id": "6ae09c13", | |
| "metadata": {}, | |
| "outputs": [ | |
| { | |
| "data": { | |
| "text/plain": [ | |
| "<Axes: xlabel='time [BMJD]', ylabel='RV [$\\\\mathrm{km\\\\,s^{-1}}$]'>" | |
| ] | |
| }, | |
| "execution_count": 78, | |
| "metadata": {}, | |
| "output_type": "execute_result" | |
| }, | |
| { | |
| "data": { | |
| "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3AAAANwCAYAAABu6ojXAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAWJQAAFiUBSVIk8AAAYAVJREFUeJzt3QeYXGX9P+wnBQIBEjqE0DF0EEhQQKoSitIFaRFFgyhiCSjiDyTSi4AIKkUjSpMOKogakRbA0IIQQgkiSegQkgCBEBL2vb7nemf/s5vdzfaZZ+e+r2uyM3PKPDM5O3s+52m96urq6hIAAABVr3elCwAAAEDrCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZKJvpQtA3tZcc830zjvvpLXWWqvSRQEAgCz873//SwMGDEgvvvhim7cV4OiQCG9z5sxJuZs8eXLxc8iQIZUuCjXA8UZ3c8zRnRxvdKfJmR5vHTl/FuDokFLN26OPPppyNnz48OLn2LFjK10UaoDjje7mmKM7Od7oTsMzPd6GDh3a7m31gQMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmehVV1dXV+lCkK+hQ4cWPx999NFKFwUAAHr8ObQaOAAAgEz0rXQBaJ2LLrooDRw4MB122GGdsr8XX3wxXXPNNemll15Kiy66aJo/f37aaqut0v7775/69evXKa8BAAB0LjVwGXjjjTfSv/71r07b33333Ze+853vpGHDhqVf//rX6YILLkhnnnlmeu6554rnZ82a1WmvBQAAdB4BropF98SJEyem448/Pr333nudss9Jkyals846K+28885pl112qX9+ySWXLF7nzTffTKNHj05z587tlNcDAAA6jyaUVejGG29MDz74YFETNnv27PTWW291yn6jmeQ555xThLPdd999geWLL7542mGHHdLtt9+ebrjhhnTooYd2yusCAACdQw1cFYp+aOedd1767W9/m37+85932n7vvPPO9PLLL6f+/funddZZp8l1Nt100+LnTTfdVIRHAACgeghwNSRq1sJqq62W+vTp0+Q6a621VvHz3XffTffcc0+3lg8AAGiZAFcjog/dM888U9xfeeWVm12vfNkjjzzSLWUDAABaRx+4GvHUU0+ljz/+uLi/9NJLN7te9INbbLHF0pw5c9KTTz7Zqn1Pnjw5DR8+vN1lGzt2bLu3BQCASujI+W+cPw8ZMqRd26qBqxHTp09vENJaUloeg6jEwCcAAEB1UANXI8rnduvbt+X/9tLymMYgtlt22WVbXD+uHqhFAwCgloztwPnv0KFD272tGrga0a9fv/r7paaUzSlfXr4dAABQWQJcjVhmmWXq73/00UctrluaxHvRRRdNSyyxRJeXDQAAaB0BrkbE1AEl77zzTou1b6X531ZdddVuKRsAANA6AlyNiIm7S6NPlg9o0tjbb79d34Ry2LBh3VY+AABg4QS4GtGrV6+09dZbF/dffvnlZtd76aWX6u9vs8023VI2AACgdQS4GnLggQemPn36pNdee63BqJTlSpN9b7bZZmmjjTbq5hICAAAtEeB6kFdffTUdccQRaa+99kpXX331AssHDx6c9txzz2J6gPHjxze5jwcffLCYRuDrX/96N5QYAABoCwGuyn344YcLjA7ZnHHjxqUXX3wxffDBB+n6669vcp0jjzyyqF279tpri/XKPfzww2nSpElp1KhRaf3110+5iPcUNwAA6OlM5F2FbrzxxqImLLz11lv1z//pT38qAlb0Z4tasrPPPrvBdttuu236xz/+kV5//fV0wAEHNLnv2O7kk09Ov/jFL9L3v//9tO+++xZTDDz99NPFtscdd1waPnx4F79DAACgPXrVRXs6atLUqVPT448/XkwbMGjQoLTlllu2ed630izyjz76aKqUUu3bpZdeWrEyAABAd5xDq4GrYauvvnpxAwAA8qAPHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcDXmo48+SjNmzGhxnTfeeKPbygMAALRe3zasSzd76qmn0rXXXlsErj59+hTP7bTTTmnPPfesf9xWN9xwQ7r++uvTLrvskoYNG5YGDx6cevfuXbzG5MmT0z333JM23HDDNHLkyE5+NwAAQEcJcFXq1ltvTWPGjEknnXRS2nLLLYvnpk+fnk4++eR0//33p9NPPz0tuuiibd7vvHnz0uzZs9Mtt9xS3Motssgi6ctf/nI6+OCDO+19AAAAnUcTyio0bty4dPHFF6dDDz20PryF5ZZbLp1wwgnpiSeeSGeddVaqq6tr1/433njjtN5666V+/foVoS1q4fbaa690ySWXCG8AAFDF1MBVmagdO//884twtttuuy2wfKWVVkqbb755uu+++9Kdd96Zdt555za/Rmx/2GGHdVKJAQCA7qIGrsrcfPPN6d13302rr756WnrppZtcZ9NNNy1+XnXVVWn+/PndXEIAAKBSBLgqc/vttxc/11hjjWbXWWuttYqfL7/8cnr88ce7rWwAAEBlaUJZRf73v/8VA5WElVdeudn1Bg0aVH//kUceSUOHDm3zaz399NPpjjvuSDNnzkyvv/56euedd9IGG2xQjHAZTSwBAIDqI8BVkSeffLL+fnPNJxsvmzhxYptf5+67707vvfdeOvroo+tHsnzhhRfS6NGj03HHHZf23XffdNRRR7V6fzH9wPDhw1N7jR07tt3bAgBAJXTk/DfOn4cMGdKubTWhrCKl2rew+OKLN7te//79m9ymtZZYYon0rW99q8E0BGuvvXb60Y9+VNyP6QViGgMAAKC6qIGrIrNmzaq/37dv8/815cuiCWRbxLQBW2+9derVq1eT0wusttpqadq0aenKK69Mu+++ezHVwMLE1QO1aAAA1JKxHTj/bU8XqBI1cFWkPCx9/PHHza5XvmyxxRZr02tE88gDDjig2eXrr79+8TP6xMV8cwAAQPUQ4KrIMsssU3//o48+ana9uXPntqqvXFMi8EUtXGvKEDVxOdRaPvPMM+mxxx5LY8aMaVCLCQAAPY0AV0VWXXXV+vtRA9ac8mXR5LEt5syZk8aPH5+mTp3a5PLevXs3GRSr0emnn54GDx6c7r333mI0zpEjRxaP43kAAOiJ9IGrIjF8f58+fYrJuVsanKR8WVvaz3744YfFyJNTpkwpXufSSy9dYL658j51K664YqpWEdJOPPHEBZ6fPXt2/fMnnHBCBUoGAABdRw1cFYnRIT/5yU/WT9LdnFLTxqgt22qrrVq9/6ipivAWIiQ2VQv39ttv199v79CmXS2aSZ555pktrhPLW6rFBACAHAlwVeaQQw6pnxuiuYFMos9Xae6JttSSlYezLbfcMn3qU59qsDz63U2aNKm4v8kmm7S5eWZ3ufHGG4uatpbE8lgPAAB6EgGuykQNXNSqRQApn9i7JGrOHnrooaK2bsSIEQ2Wvfrqq+mII45Ie+21V7r66qsX2LbUXHLYsGFFE8TGUwTEQCAxwXcMcvLNb34zVat4n525HgAA5EKAq0LHH398WnPNNYu52CKwlfvrX/+aZsyYkU466aS08sorN1g2bty49OKLL6YPPvggXX/99QvsN2rdollkDE7SeB64eO6KK64o5pgbNWpUWnfddVO1GjRoUKeuBwAAuehVV1dXV+lC0HRzx3POOacIVl/4whdS//7904QJE9KDDz6Yvvvd7xZhrKkapwh2r7/+evrSl760QA1diPD3s5/9LC266KJpu+22SwMHDkxvvfVWuvXWW4s+dd/+9rfTRhtt1OpylgZRefTRR1N39oGL0SZbakYZNZSvvPJKGjBgQLeVCwAAuvocWoCrcs8991yaOHFiEeRWX331ovljhK+OigFMYr8RhmLut7XXXrtdtW6VCHAtjUJZctpppxmFEgCAqtSRc2jTCFS5CFVd0ZwxwmDcclUKZzHaZHlNXNS8/fjHPxbeAADokdTA0SGVqoEriakC9txzz/T++++nb33rW2n//ffXbBIAgKqmBo6aFWFt/fXXL+5/7Wtfq3RxAACgSxmFEgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkIm+XbHTK664IlWbww47rNJFAAAAqL4A99WvfjX16tUrVRMBDgAAyF2XBLiSurq6VA2qLUwCAABUVYBbYYUV0kEHHZSWWWaZVClTp05NV111VZo3b17FygAAAFD1Ae7HP/5x+t73vpcqbdasWemWW26pdDEAAACqdxTKTTbZJFWDT37yk5UuAgAAQPUGuOhzttFGG6VqsPHGG1e6CAAAANUb4JZeeum00korpWoQQbJaBlMBAACougD3+OOPp2oxZMiQ9Pzzz1e6GAAAANUZ4FZbbbVULaI559prr13pYgAAAFTvICYAAAB0LgEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMlE1AW733XevdBEAAACqWtUEuOeee67SRQAAAKhqVRHg5syZk1566aVKFwMAAKCq9e3Ixtddd1368MMPO1yI+++/P82bN6/D+wEAAOjJOhTgbrjhhnTzzTenXr16dV6JAAAA6PwAd+ONN6aTTjopnXbaaWn55ZdPSyyxRLv2M3369DR79uyOFAUAAKDH61CAC6ecckp66KGH0mc+85n0k5/8pF37ePrpp9PGG2/c0aIAAAD0aJ0yiMkPf/jDdN9997V7+w022CAtu+yynVEUAACAHqtTAty2226bXn311Q7tY4UVVuiMogAAAPRYHW5CGfr165dGjhzZoX2cffbZnVGUHuWpp55K1157bZoxY0bq06dP8dxOO+2U9txzz/rH7fXiiy+ma665ppi+YdFFF03z589PW221Vdp///2L/08AAKCHBrjwve99r0PbRyjh/7n11lvTmDFjikFittxyy/rBXk4++eRi2oXTTz+9CF7tEc1dzznnnPSd73wn/d///V/x3HvvvZd+9rOfFc/Fz4EDB3bq+wEAAHrIRN40NG7cuHTxxRenQw89tD68heWWWy6dcMIJ6YknnkhnnXVWqqura/O+J02aVGy78847p1122aX++SWXXDIdf/zx6c0330yjR49Oc+fO7bT3AwAAdA4BrsrEdArnn39+Ec522223BZavtNJKafPNNy9q0e6888427TuaSUbNW4Sz3XfffYHliy++eNphhx2Kppsxxx8AAFBdBLgqExOjv/vuu2n11VdPSy+9dJPrbLrppsXPq666qghlrRWB7+WXX079+/dP66yzTov7vummm8zNBwAAtRzgHn744e58uSzdfvvtxc811lij2XXWWmut4meEsccff7zN+15ttdWaHQSltO8Ikffcc0+byg4AAPSgAHf00Ud358tl53//+18xUElYeeWVm11v0KBB9fcfeeSRVu07Bil55plnFrrv8mWt3TcAAJDZKJSt8fHHH3fny2XnySefrL/fXPPJxssmTpzYqn1Hv7bS59/SvqMf3GKLLZbmzJnToDwtmTx5cho+fHhqr7Fjx7Z7WwAAqISOnP/G+fOQIUOqvwauV69e3fly2SnVvpWCVHOiD1tT23TGvsuXz5o1q0197AAAgB5UA0fLIjCV9O3b/H9N+bKZM2d26r7Ll8dImLHdsssu2+L6cfVALRoAALVkbAfOf4cOHdrubY1CWUX69evXquam5cuiuWNn7rvx8vLtAACAyhLgqsgyyyxTf/+jjz5qdr3ySbZb6s/Wnn2X73/RRRdNSyyxRKv2DwAAdD0Broqsuuqq9fffeeedZtcrXxZTArRG+Xot7Ttq30rzv5WXBwAAqDwBropsvvnm9fOztTQ4Sfmy1rafjYm7S7V1Le377bffrm9COWzYsFaXHQAA6HoCXBWJ5oqf/OQn6yfpbs60adOKn717905bbbVVq0cA3XrrrRe675deeqn+/jbbbNPqsgMAAF1PgKsyhxxySP3cEM0NNlKakDvmnlhxxRVbve8DDzywqOF77bXXGoxK2dS+N9tss7TRRhu14x0AAABdRYCrMlEDF7Vq0Q+tqYm0Y162hx56qKitGzFiRINlr776ajriiCPSXnvtla6++uoFth08eHDac889i+kBxo8f3+TrP/jgg8U0Al//+tdTLi699NLiBgAAPZ0AV4WOP/74tOaaa6Yrr7xygYm0//rXv6YZM2akk046Ka288soNlo0bNy69+OKL6YMPPkjXX399k/s+8sgji9q1a6+9tliv3MMPP5wmTZqURo0aldZff/0ueGcAAEBH9KqL6phu8qlPfaqoPWLhYjCRc845pxjS/wtf+ELq379/mjBhQlFD9t3vfjdtueWWC2wTNXAR7F5//fX0pS99aYEaupL3338//eIXvyjC3r777ltMMfD000+nf/zjH+nwww8vmma2VmkQlUcffbQD7xYAAGrH0A6cQwtwVe65555LEydOLILc6quvXowMGfOzdYapU6emxx9/vGiuOWjQoCIUtnXeNwEOAABSt51D923zFnSrddddt7h1hQiEcQMAAPLQrX3gurGyDwAAoMfp1gA3cuTI7nw5AACAHqVbA1yMgAgAAED7mEYAAAAgEwIcAABAJgQ4AACATGQzjcDHH3+c7rjjjmLC6Zir7DOf+UzadNNNK10sAACA2glwzz77bBo/fnyD53r16pW+8IUvpGWXXbZ4HBNZf/GLX0zPP/98g/X22GOPdOWVV6YBAwZ0a5kBAABqMsC99NJL6Ygjjkjz5s0rHi+11FJpww03TJtvvnkR4GL5jjvumGbMmFHMI7faaqulT3/60+mVV15Jt912W9p3333TnXfeWem3AQAA0PP7wH3uc59Lp5xySlpsscXSxRdfnKZPn54eeOCBtPHGGxfLR40ald5+++3i/sEHH5wmT56crr/++jRu3Lh0yy23pHvvvTf96U9/qvC7AAAA6Hq96qJaq4Li5TfaaKN08sknpwMOOKDBsilTpqS11167uL/mmmump556qgh65Q488MDi53XXXdeNpaZk6NChxc9HH3200kUBAIAefw5d8Rq4u+++O82fP3+B8BZuv/32IuCFkSNHLhDewogRI9IjjzzSLWUFAACopIoHuMcff7w+gTb28MMPN2hq2ZSooXv11Ve7rHwAAADVouIBbu7cualPnz5NLnvyySeLn7179252yoColWtuewAAgJ6k4gFuueWWK0aUbOyDDz4o+rzFlALrrbdek80nw6RJk9KQIUO6oaQAAAA1HuA++clPFvPARWArFyNLfvjhh8X9nXbaqdntf/Ob3xRzxAEAAPR0FQ9wMd9bzP122WWX1T83a9asNHr06PrHhxxySJPbnn/++cUk39/73ve6pawAAAA1PZF337590zHHHJOOPfbY9O9//zsNHjw4/fnPf07//e9/i+aTu+++e9p6660bbPPcc88VAS/mg4tJwJdccsmKlR8AAKBmAlz44Q9/WNSkXXnllUVoK00dsMUWW6Q//OEP9ev94he/KNaZMGFC/TrRhDL6x11wwQUVKz8AAEDNBLgQQe3www9Pd955Z9GEctiwYUXTyaihK5kzZ07aZJNNilu5Ul85AACAnqxXXakqC7p5FnkAAKhFQztwDl3xQUwAAABoHQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZCL7APf222+nr33ta5UuBgAAQJfLPsC99tpr6Q9/+EOliwEAANDlsg9wEydOrHQRAACAjBx55JHFLUd9UxV588030z//+c/0yiuvpNmzZy90/blz56arrrqqW8oGAABQaVUR4CKIHXvssenSSy9N8+fPb9O2dXV1qVevXl1WNgAAgGpRFQHuuOOOS7/+9a+LMAYAAECVBrjp06enX/3qV6l3797pM5/5TNphhx3SCiuskJZaaqlW1dzFtk8++WS3lBUAAKCmA1z0eevTp0+6++6701ZbbdXm7QcNGpT23nvvLikbAABANan4KJQvv/xy2nnnndsV3sK6666r6SUAAFATKh7g5s2bl5Zddtl2b7/yyiun0aNHd2qZAAAAqlHFA9yGG26Ynn/++XZvP2DAAAEOAACoCRUPcMOHD08vvPBCevrpp9u1/TvvvJNOOeWUTi8XAABAtal4gOvXr186//zz04gRI4ow1lavvfZaOvnkk7ukbAAAANWk4gEuHHLIIemAAw5IW2+9dXrooYfatO3kyZO7rFwAAADVpOLTCJQcf/zx6VOf+lTabbfd0sYbb1yMTLnqqqumvn37tjgP3C9/+ctuLScAAECq9QA3ZsyYdOqpp6aZM2em+++/v7i1Rkwh0KtXry4vHwAAQKVVRYC78MIL06hRo+ofm9cNAACgCgPc+++/XwxCEqEtmkzusMMOaaWVVkpLLbXUQreNJpRXX311mjZtWreUFQAAoKYD3F133ZVmzZpVBLGDDz64zdtvttlm6aCDDuqSsgEAAFSTio9COWXKlKLWrT3hLcSAJwAAALWg4gHu3XffTauttlq7t4/mlocddlinlgkAAKAaVTzArbPOOmnq1Knt3n7ZZZdNl19+eaeWCQAAoBpVPMDFvG9PPPFEevXVV9s9CMoVV1zR6eUCAACoNhUPcEsuuWQaPXp0Ovzww9O8efPavP0rr7xSbAsAANDTVTzAhe985ztpvfXWS7vuumt6+eWX2zwICgAAQC2o+DQCJb/4xS+Kvmxbb7112mmnndLOO++cBg8enPr27dviPHBnnHFGt5YTAACgZgPcueeem370ox/VP44Jva+66qri1hqxfq9evbqwhAAAANWh4gFu+eWXL0JYucaPAQAAqJIAF6IWbZNNNikm5l5kkUVatW00ofzb3/6WZs6c2cWlBAAAqLyqCXCf//zn01/+8pc2bx/95kaOHNkFJQMAAKguvaslwG233Xbt2n7LLbfs5BIBAABUp4oHuJVWWintsMMOad11123X9iussELafvvtO71cAAAA1abiTSiXWmqpdNddd3UoAHZkewAAgFxUvAauo95///10xRVXVLoYAAAAtRHgTjnllDRr1qx2bfvCCy+kww8/vNPLBAAAUG2qIsCdfPLJ6ZlnnmnXtg8//HCnlwcAAKAaVUWAi4m777nnnjZv9+KLL6af/vSnXVImAACAalMVAa40n9vHH3/c6vVvu+229KlPfSpNmzatS8sFAABQLaomwD333HPpl7/85ULXmz17djriiCPS3nvvnaZPn94tZQMAAKgGVRPgwvHHH58mTJjQ7PJ77703bbrppul3v/td0ewybgAAALWiagJchLE5c+akvfbaK7388ssNls2dOzcde+yx6bOf/WzR7y3WHTx4cLrsssuK0AcAAFALKj6Rd3mfts022ywdeuihadddd033339/GjhwYHrsscfSYYcdlp5++un6GrcRI0YUzS0HDBiQ3nrrrfTqq69WuvgAAAC1EeDKBy/561//mnbaaaeiJi5q3M4444w0b968Irwtv/zy6ZJLLkn77bdf/frxXAyAAgAA0NNVTRPKksUXX7wIca+//noxwfdHH31UhLc999wzTZw4sUF4AwAAqCVVF+DCsssum/7+97+nlVdeOfXt2zeNGTMm/elPf0orrrjiAutGf7m11167IuUEAABItR7gwhprrJHuuOOOtMQSS7TYx+3NN99MU6ZM6dayAQAA1GwfuObElAG33npr+vznP59WW2219OUvf3mBde65556KlA0AAKDHBLiYaPvJJ5/stP2NHDky/eY3v1lgeoGW5o0DAADoSboswMWcbjFnW69evYrHHZl0O/YR28fUAo3F86XXAAAA6Mn6duVAJKWAteSSSxaPOztoRQ2cOeAAAIBa0WUBbrnllisC26mnnpqOP/741Lt314yXcvLJJxfTDQAAAPR0vbu6Bu6oo47qsvAWvvCFL3TZvgEAAKpJlyWrZZZZpmg+OXDgwNSVll9++Q71rwMAAEi1HuBGjRqVPv744y4fYGTNNdcsXgcAAKCnq9qJvAEAAGhIgAMAAKjlADd+/PhUTaqtPAAAAFUT4HbZZZdUTaqtPAAAAFUT4N599900bdq0VA2mTJmS3nvvvUoXAwAAoHr7wE2aNClVg4kTJ1a6CAAAANUd4B5//PFUDSZMmFDpIgAAAHSKvqmLnHLKKenf//53MaF3pcyYMSP94x//qNjrAwAAZBHg5syZk/785z+nSqurq+vyycQBAACyDnARnAAAAKjyAHfXXXd1xW4BAABqWpcEuB122KErdgsAAFDTumwUSgAAADqXAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAajnAbbDBBl2xWwAAgJrWJQHuhRde6IrdAgAA1LQuCXAfffRRmjp1alfsGgAAoGZ1WR+4n/zkJ121awAAgJrUZQHuqquuSp/73OfS/fff31UvAQAAUFO6dBTKu+66K22//fbpE5/4RDruuOPSAw880JUvBwAA0KN1WYCrq6tLQ4cOTcOHD0/9+vVL5513Xtpuu+3SKqusko466qg0duzYNG/evK56eQAAgB6nywLc5Zdfnh5++OH0t7/9LT311FNp+vTpRbPKbbfdtvi52267pRVXXDEddthh6ZZbbkkffPBBVxUFAACgR+iSANenT5906KGHNnhu6aWXTgcffHC6/vrr05tvvlmEtr333jvdcccd6Ytf/GJafvnl07777puuvPLKNHPmzK4oFgAAQNb6dsVOF1lkkdS3b/O7jiaVe+21V3GbP39+uueee9LNN9+c/vSnPxW32D76zu23335pn332SYMGDeqKYgIAAGSlS2rgVl999TbV1n32s59Nv/zlL9O0adPSgw8+mEaNGlXc//a3v51WW221tPXWW6ef/exn6fnnn++K4gIAANRugHvmmWfave2nP/3pdNZZZxX7ePLJJ9Po0aPThx9+mH70ox+l9dZbL22yySbFc48//ninlhkAAKCmpxHoqI022qiYEPyxxx5Lt956a1pppZWKAVFOO+20YoTLtddeOx177LFp3LhxxaiXtM4bb7zR4vIZM2YYIRQAAGqlD1xnefHFF9M111yT/vjHP6ZJkybVP18Ka1OmTEkXXHBBcTv77LPTD37wg9RTzJo1K1177bXpP//5T1p00UXT3Llz07rrrlsMDrPCCiu0e79vvfVWGjFiRBo2bFjacccd05AhQ1L//v3T+++/XzRbjZFD//vf/xa1oAMGDOjU9wQAAPSwABe1Q9ddd10R2saPH98gsPXq1au4lT+3zTbbpEMOOSQdeOCBqad46aWX0g9/+MOi71+E0whwMdhLfCZHHnlkUQO54YYbtmvfUbMWn10Etbg1tsUWWxT7F94AgJ4izp/CpZdeWumiQHUGuDPOOCP93//9X6vXf/fdd9NNN91U1LbdfffdRVhpHNzKH2+22WbFlAQHHXRQMchJT/LOO++kE044IS255JLp6KOPTr17964f7CVqzkr9Ai+88MJ2j8651FJLpc033zxNmDAhzZ49u5jiYZ111km77757MU9f6fMGAABqIMBFDc7CAlwMTHLbbbcVoS3mgovHTdW2xeO4RVO/CG1xi8FMeqrLLrssvfLKK+kb3/hGfXgrt+uuu6Yzzzwz/fznP0/nnHNOu14jmkxG30IAACAvXRLgmhtQ5OOPP07//Oc/i9AWg5JEzVtLoW3w4MFF08hoIhlN+3q6CG7x+YRNN920yXVKz0ft2RNPPNHsegAAQM/TJQHuo48+KprmLbHEEsXjBx54oOi/dcMNN6Q333yzxdC23HLLpf3337+oadtuu+1qqjnf3/72t6L5aLzn5ubSW3755YsmkBF+b7/9dgEOAABqSJcEuKhp22uvvYpmj//4xz+K0SJbCm3R32ufffYpQtvw4cNT375VN7ZKt3jkkUeKnwMHDkyLL754s+tF37cIcDG9Qnx+7Qm50dfu5ptvLgZMiZEpX3/99bTiiium7bffPu2xxx6pX79+HXovAABA5+uypBSDkcStudAWASEGzYjQtueee6bFFlss1bIPPvggPf/888X9GFSkJaXlM2fOTFOnTk1rrLFGm8NbTBNw1FFHpVVXXbX+9S+66KJ0ySWXFDWBMRBNa6crmDx5chG822vs2LHt3hYAACqhI+e/cf4clV1VOZF3+QiSMShHvNHf/e53RY1P1AAdcMABNR/ewttvv10fdluqfSsNQlIyffr0Nr9WhLUYTrcU3kqvGZOix6ieMf/eySefXNSkAgAA1aNL2yqWmvfFXG1R0xZhLZrp0fTE3SULa0Javrx8u9ZYZJFFisFhmqq1i6kKdtlllzRmzJj07LPPpvvuuy/tsMMOC91nXD1QiwYAQC0Z24Hz36FDh1ZfgIsBTI455ph0+OGHt7mJXy0q73O2sJqv8uVt7asWg8REU8nmrL/++vX3x40b16oABwAAdI8uC3DXXXdd0cetJ4j51rqyhin2vcwyyzQYxbMlc+fOrb9fvl1rxSiWzSnf37Rp09q8bwAAILMAF038dtttt9RTRH+xww47rEtfIwYmKU0PEIOMtKS0PPoURnPItormkbGPYcOGLTCCZfnk4eVBEQAA6MEBrifN3xbD+setK0Vw2myzzYp+ZzNmzGhxeoDSwCXR92zAgAFtep0rr7wyXXHFFcX9mCA9mriWi5EtS1o7CiUAANA9umQUytJoirTNZz7zmeLnhx9+WD/heWPRvPK1114r7sfgMG0Rk4THhOolL7zwQpOjYZasu+66bdo/AACQYYCLybtpux133DGtssoq9c0cm5szIoJYNLfce++927T/aDZZ6l8XTTZHjBixwDr/+c9/ip9R+9dT+jACAEBP0SUBbrvttuuK3fZ4MYz/yJEji/v//ve/m1yn9HyErxjps7Hzzz+/mBh91KhR6b333muwLEJb3KK55q9+9au03nrrLTC65f3331/c32+//erDJAAAUB26fCLvzhKTS//lL38pAsbUqVNTTxXh99BDD0133313+t///tdgWTSrvP3229Ouu+5aBKymatjuuOOONGfOnDRx4sQ0YcKEBsujVi3m44ugFrV4jd10001FE8qtttqqPkgCAAA1MpF3Z3r11VfTtddeW4SYJ598shgoJZocfv/73+/QRHjV6Ktf/Woxv9uJJ55YNJNcc801i/d/8803FxNtNxeuYkCTaPZ41113pU984hNp8803X2CdCH7RjPKkk04q1o1RLCPMPfjgg8Xt61//ejrwwAN71CA0AADQU/Sqy3TEkRtuuKGoTYpmh7fddlsaPnx46mliNMqHH344vfXWW0XTxxj2f8UVV+yUfc+ePTs98sgj6Y033ig+wwhyEfgWXXTRNu2nFJ4fffTRTikXAEBXTAkVLr300koXhSpxZIWPiY6cQ3d5Ddy8efOKiapj8I24H0PfRxPAtgaFxg444ICin9crr7ySRo8e3SMDXEyqHTVuXSH6z+2www5dsm8AAKBrdGmAi/5a0SSv8ZD4yy23XDGIRoSw9njuuefSZZddVoS3EE0qAQAAerouC3AxIfU+++zT5GAZ0SQwmj9+8MEH6bDDDmt2H5MmTSqqFZ9++ulizrK4RXh79913G6y3yCKLdMl7AAAAqIkA9+1vf7sIb6XBMEpd7UqPYyTEo446Km277bZp7bXXrt8uRkGMofB///vfFwN3lGvcXa+0ry233LKr3gYAAEDPDnD33HNPMYx9BKwIXdFkMkJazD/28ssvp5deeqlYL2rgov/alVdeWTweP358UWsXA2s0FdYaj4wY60RfulNOOaUr3gYAAEDPnwcu5msL6667bvrnP/9Z9IGLcBbD1MccbjG/2fe+970ikF1//fVp5syZxfMxrH2Et6ZEWCu/RXCLIfFjlMZPf/rTXfE2AAAAen4NXPRbi7nLHnjggWIkxcbWWGON9POf/zxtv/32xUAmMUrlLbfcUgS5ENtsuummaaWVViqGze/fv38R2OL5QYMGFSNZbrLJJh0eyRIAACDVeoCLKQPOPPPMJsNbuX333beYNDqaXN54443pc5/7XDrmmGOKKQFiom4AAAD+ny5JSVGTFs0hW+PLX/5y2nvvvdOXvvSldNVVV3VFcQAAAHqELukDFxN2L7/88q2ehfyjjz5KJ598clcUBQAAoMfokgAXo022VoxQ2adPn7TOOuu0+/V22223dm8LAABQ0wGuTQXo3bvDE3Hfe++9nVYeAACAalXxAAcAAEDrCHAAAAC1PArlnDlz0u23317M1bYwpYm5p02bVvxsqwkTJqQPP/ywnSUFAADIR5dNtrbXXnu1af2Y+BsAAIDmdels2a2tUevVq1e7at/KtwcAAOjpqqIPXEfCGwAAQK3oshq4ffbZJ2266aapqz3++OPpz3/+c5e/DgAAQI8McDGv24033tgtTRvnz5+f+vfv3+WvAwAA0CObUPbp06fb+qXFa8UNAACgp+vdE/q06UMHAADUgi4JcCeeeGJX7LZqXg8AAKDHBLgTTjihK3ZbNa8HAABQs9MIAAAAsHACHAAAQCYEOAAAgEwIcAAAQM2YNWtWeuaZZ9Jjjz2WxowZUzzOiQAHAADUhNNPPz0NHjw43XvvvemRRx5JI0eOLB7H87noW+kCAAAAdLXTTz+9yenHZs+eXf98DqPbq4EDAAB6tFmzZqUzzzyzxXVi+TvvvJOqnQAHAAD0aDfeeGNR09aSWB7rVTsBDgAA6NFeffXVTl2vkgQ4AACgRxs0aFCnrldJAhwAANCj7b///mmJJZZocZ1YHutVOwEOAADo0QYOHJh+/OMft7hOLB8wYECqdqYRAAAAerwT/v8pAmK0yfIBTaLmLcJbDlMIBDVwAABATTjhhBPSK6+8krbffvs0bNiwNGbMmOJxLuEtqIEDAABqxoABA9L6669f3P/a176WcqMGDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAD3WrFmz0jPPPJMee+yxNGbMmOIx5EyAAwCgRzr99NPT4MGD07333pseeeSRNHLkyOJxPA+56lvpAgAAQGeLkHbiiScu8Pzs2bPrnz/hhBMqUDLoGDVwAAD0KNFM8swzz2xxnVj+zjvvdFuZoLOogcvAs88+m0499dR01VVXddo+P/jgg3TjjTem8ePHp759+6a5c+em1VdfPR188MFpjTXW6LTXAQDobnGOEzVtLYnlsd7Xvva1bisXdAYBLgPXXnttp+5v5syZ6Qc/+EFaddVV09lnn52WWGKJ4vm//e1v6eijj04//vGP0zbbbNOprwkA0F1effXVTl0PqokmlFVsxowZ6dxzz03jxo3rtH1GTdvo0aOLfR9//PH14S3stttuaccdd0xnnHFGMVoTAECOBg0a1KnrQTUR4KrMlClT0rHHHpu+853vpO9+97vpX//6V6fu/7rrrkuTJk1KO+20U1psscUWWL777runDz/8MJ1zzjlp/vz5nfraAADdYf/9929wkbopsTzWg9wIcFUm+p+dd9556aKLLkpXXnll2mCDDTpt39HW++abby7ub7LJJk2us9566xXBbtq0aZ0eHgEAusPAgQOLLiEtieUDBgzotjJBZxHgasg999yT3nvvveL+mmuu2eQ6ffr0Sauttlpx/7bbbuvW8gEAdJaYIuC0005boCYuHsfzphAgVwJcDXn44Yfr76+88srNrldaFv3gSoEPACA3EdJeeeWVtP3226dhw4alMWPGFI+FN3JmFMoaMnHixOLn4osvnvr169fseksvvXTx8+OPPy76y33qU59qcb+TJ09Ow4cPb3e5xo4d2+5tAQBaEs0k119//eK+KQPoTB05/43z5yFDhrRrWzVwNWLevHnFpJalANeS/v3719+fPn16l5cNAABoHTVwNSLCW11dXXE/Ju5uSfnyUuhrSVw9UIsGAEAtGduB89+hQ4e2e1s1cDWifMqAaBrZkvLlLTW1BAAAupcauFaIOdG6soapO2qvYsSlRRZZJH300UfFbWGTfZcss8wyXV42AACgdQS4VjjyyCPTYYcdlnIX0wO88MIL6d133y2aU/bq1avJ9d555536+6uuumo3lhAAAGiJANfKySDjlrtoaxsBLppIzpgxIy277LJNrlcauCRq39ZZZ51uLiUAANAcfeBqyDbbbFN//+WXX252vZdeeqn4ufXWWzdbSwcAAHQ/Aa6GbLzxxmnTTTetn6S7KW+//XZ64403Up8+fdKBBx7YzSUEAABaIsD1MOeff37ac88906hRo9J77723wPKRI0cW0wSMHz++ye3//e9/Fz/32muvtMoqq3R5eQEAgNYT4KpcaUTI8pEhmxODj9xxxx1pzpw5aeLEiWnChAkLrLPBBhuk7373u+k///lPeuyxxxose//999P111+ftthii2LgFgAAoLoYxKTKTJkyJV144YXF/Xnz5qXnnnuuuB+DjnzrW99K/fv3Lx6PGDEibb755g22HTBgQNp9993TXXfdlT7xiU8ssLwk1om+beeee27adddd0/rrr1/s/+abb04bbrhh+s53vlM0oQQAAKqLAFdl1lhjjXTeeee1e/tjjjmmuC3Mbrvtlrbbbrv00EMPFSNTxjxxJ510kmkDAACgiglwNSxC20477VTpYgAAAK2kDxwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIRN9KFwAAAKA7XXrppSlXauAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAABAJgQ4AACATAhwAAAAmRDgAAAAMiHAAQAAZEKAAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAFeD3njjjRaXz5gxI82bN6/bygMAALRO31auRwU9++yz6dRTT01XXXVVh/f11ltvpREjRqRhw4alHXfcMQ0ZMiT1798/vf/++2natGnp4YcfTv/973/TWWedlQYMGNAp5QcAqKRLL7200kWATiPAZeDaa6/ttH1FzVpdXV0R1OLW2BZbbJFOO+004Q0AAKqQAFfFoinjmDFj0rhx49JKK63Uaftdaqml0uabb54mTJiQZs+enZZeeum0zjrrpN133z1tu+22qVevXp32WgAAQOcR4KrMlClT0oUXXpjmzp2bZs6cmaZPn97prxFNJn/yk590+n4BAICuJcBVmTXWWCOdd9559Y+PPfbY9MQTT1S0TAAAQHUwCiUAAEAm1MDVqHfeeSfdfPPN6aWXXipGpnz99dfTiiuumLbffvu0xx57pH79+rV6X5MnT07Dhw9vd1nGjh3b7m0BAKASOnL+G+fPMRp8ewhwNRreYpqAo446Kq266qrFcx988EG66KKL0iWXXJL+9re/pTPOOCOtsMIKlS4qAABQRoCrQRHWjjzyyPrwFhZffPGiv90zzzyTXnzxxXTyyScXg6n07r3wVrZx9UAtGgAAtWRsB85/hw4d2u5t9YGrMYssskgaPHhwMVhKY3369Em77LJL/eTh9913XwVKCAAANEeAqzHLLbdc0VSyOeuvv379/Zh/DgAAqB6aULbCOeec06VNBLu7+WFM5N2cZZZZpv7+tGnTuqlEAABAawhwrRD9xQ477LDUU0TzyBjIZNiwYalXr14NlpX3eYvJxAEAgOohwLXCwIEDi1tPcOWVV6YrrriiuH/IIYekww8/vMHymTNn1t83CiUAAFQXfeBqyPz589Mf//jH+scvvPDCAuu8/fbb9ffXXXfdbisbAACwcAJcDYlmkx999FFxf+mll04jRoxYYJ3//Oc/xc9oWrn77rt3exkBAIDmCXA9zPnnn5/23HPPNGrUqPTee+81WBahLW7Rz+1Xv/pVWm+99Ros//jjj9P9999f3N9vv/3SKqus0q1lBwAAWibAVbnSQCKtGVAkatjuuOOONGfOnDRx4sQ0YcKEBsujVu3ggw8uglo0p2zspptuKppQbrXVVmnkyJGd+C4AAIDOYBCTKjNlypR04YUXFvfnzZuXnnvuueL+jBkz0re+9a3Uv3//4nE0f9x8880bbDtgwICi2eNdd92VPvGJTyywvFSzFs0oTzrppGLdmNQ7wtyDDz5Y3L7+9a+nAw88cIHRKQEAgMrrVVdXV1fpQtD9Zs+enR555JH0xhtvpD59+hRBLgLfoosu2qb9DB06tPj56KOPdlFJAQCgZxnagXNoNXA1aokllkg77LBDpYtRNYYPH16RSdWpTY43uptjju7keKM7Da/B400fOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQiV51dXV1lS4E+Vp22WXTnDlz0gYbbJByNnny5OLnkCFDKl0UaoDjje7mmKM7Od7oTpMzPd6efvrptNhii6W33367zdsKcHTImmuumd5555201lprVbooAACQhf/9739pwIAB6cUXX2zztgIcAABAJvSBAwAAyIQABwAAkAkBDgAAIBMCHAAAQCYEOAAAgEwIcAAAAJkQ4KiI2bNnF7eWvPnmm00+HzNfvPXWWy1u+8Ybb3SofPQsjje6k+ON7uaYozs53iqvb6ULQG267LLL0r///e+06667ps033zytvPLK6eOPPy5mo3/mmWfSv/71r/T5z38+7bnnngtse9ddd6Xzzjsvffazn01bbbVVWmONNVKfPn2KCcVjUsRx48al/v37p//7v/9baDnmzZuX/vKXv6R777039e7du/hiWXzxxdOhhx6aNtxww4VuP3bs2HTHHXek+fPnF+WP191///3Tlltu2arPISZvvOaaa9JLL72UFl100WI/8Z5iH/369WvVPqid463cEUcckY4++uj0yU9+stXbON66R0843saPH1+81nvvvZdeffXV4vWjLNtvv33aZZddijItjOOt++R+zH300Ufp7rvvTg8++GBxnES5IwCsv/76aa+99kpbbLFFqz6Hp556Kl177bVpxowZ9cfoTjvtVLzv1hyz1Mbx1pz//Oc/6Yc//GG64oorivdUzd9xJvKmIs4555wi/DQlfnGPOuqo4ouhKX//+9/Tueee2+SyXr16pX322ScdeeSRC/2yji/4E044Ia2++urF6w0YMKB4furUqcWJ8U9/+tNm/2jEF9WZZ56Znn/++XTaaaelwYMHF88/++yzxT5333339PWvf73F17/vvvuKz+E73/lOcUIU4mTpZz/7WXHCFD8HDhzY4j6ojeOtsQceeCCNHj26KFdrA5zjrfvkfLzFKcGFF16YhgwZUpSx9DpxjES5nnjiibTmmmsW33srrbRSs6/veOteOR9z77//fvrBD36Qtt1223TggQfWv06cHB9//PFp+vTpxQlxlKElt956axozZkw66aST6i+ixrYnn3xycTJ9+umnFyfZ1Pbx1pw5c+YUr/vKK6+kK6+8cqEBrtLfcZpQUjHDhg0rTgQWWWSR4ss1rsIcdNBB6be//W2zv/gla6+9dtp0003TYostlvr27Vv8og0fPjxddNFFxS/ywn7x4xf1Rz/6UbF9/Cz94ocPPviguBp44403Nrv9pZdeWlwtPO644+rDW1hvvfXSN7/5zeIKYPwxac6kSZPSWWedlXbeeef6X/yw5JJLFn+w4spjnKDPnTu3xfdBbRxv5Vcb449m/NFoC8db98v1eLvuuuuK77G4el7+OoMGDUqnnnpqWmGFFYoT6x//+MfFfprieKuMXI+566+/Pk2ePHmBmt14LyNHjizux7bNBYYQtTYXX3xxUfNS3gJmueWWK07y48JDHJPqLDpPrsdbcyL8R3hrjWr4jtOEkorZcccdF/pL3px11lmnCE/tdfXVVxdV9b/4xS+KKz4l8UsXVwLjRDmu3DUl/hDcfPPNxRfXBhtssMDyHXbYofgSipD36U9/ujjxKRdV7HECHr/YUVPXWFT/xz5uv/32dMMNNxR/kKjd4y2cffbZ6fXXXy/WiSt8C+t7UM7xVhm5Hm/x3faHP/yh2SvrcQxF86Jp06alO++8M+22224N1nG8VU6ux9xzzz1X/Dz22GPTiSeeWNT+lkTzx/h7GrV0f/7zn4uT/Mbi+/D8888vwlnj4zFETXE084sakzhm46Sb2j3emvLkk0+mP/3pT6k1quU7Tg0cNWfWrFnppptuKq74NG4jHU0jQ7SlbqrtdoiTl7DJJps0uTyuRkW7/fgCiS+ZxuIPyMsvv1ycDMWXWFPiylSIcrblZJ2ed7yFuMIYJyhxYv2Nb3yjTa/veKstHTneYttolhQnus3VVJQfQ3ExqzHHW+3p6Hdc6TiJ2o844S0XNTGlprpR89vcRYd33323aEq39NJLt3jMXXXVVcUJOLX9N7VxbV406Wx8sb051fIdpwaOmhMdaKN6vanas/hDEQEtvuCXX375BZZHR9Xo5BqiBq45a621VnrssceKjrzf/va3iysyJXFVJqy22mrNNhOI7UP8UbrnnnuK5kzU3vHWGRxvtaUjx1t8T8WJT/TdiOZE3/rWt4qr7OWiuVNJhL3GHG+1p6PfcV/+8peLE+/ouxT9nxorHXNxcbQppWMumvA1p3TMxYn3448/noYOHdrKd0dP/5v6u9/9rgiDG220UdH3bWGq5TtOgKOiHn744aIv2cyZM9Nrr72WPvzww6Jma999903rrrtui9tOmTKlqPKOUY9iyNmoLv/EJz5RNKGIztDl1erlYpSr0i9fU5ZZZplmX/ORRx6pv99SB9fSlZz4konAF6MShWj+FiM0LWz78mXxmk5wavN46yjHW2XldrzFAA9bb711uv/++4vX/c1vfrNAgCsf/ru8/29wvFVebsdc6biLUXWbEifiMSBEaCp0RTO6UlO51vxNLh1zAlztHm/lJk6cmP7xj38Uo2r+7W9/SwtTTd9xAhwVEx1MYxjZY445pv4qRoSdGDEqqqjjC/2AAw5otr1yXLGL0X+WWGKJ4rn48jjllFOK23bbbVcMQVt+tbhxe/v4JY8/DnE1JYagLV2BjhOWuDV1ZSVet6S5phqNl8U2pQAXQxyXqvhb2j7KEZ1zo2q//DWpreOtoxxvlZPr8RbljRPq6KTfVF+hOGEL8dqNmyg53ior12OuJTFUfZw0L7XUUulrX/tak+Vu69/kOGmn43I/3j788MOi6WQMlLPiiiu26j1X03ecPnBUTDR3OPjggxv8ksWQ6NHkMPpfxBWR0pWWxuIXKL40Sr/4pSseP/nJT4pf+OjDESMhNRZXj+OPQYirOzFaUDSpiCGx4xc55puJvkbRCbapdsvlnWLLm0U2Fm2jm9qmtduXL4/23trs1+bx1lGOt8rJ9XiL0dzixCn6Cn31q19tsCz6IMUJdTSzHDVq1AJN1hxvlZXrMVeqSYvXj36+MWBJiP3EYGBRKxP7aKqPUkf/JlObx1u4/PLLi+C2xx57pNaqpu84AY6KiCu822yzTZPLYvSe0pdtU7/A0Q4+hq9taj6X+IKP0aZCVM3HiETl4hepJP5QxJdH9POIKyUhhh+Oq8pxhS6GJG6sfPumrgw1tax8m9ZuX748vgjLt6N2jreOcrxVRk883qIvR8x9WTppKh86u6nXd7x1r9yPuTjZjgsEMY1F1IjEUPLxM/o0xZQVzfU5b8/f5GjuR20fb0899VTRZDJGPm2LavqOE+CoiKhab9y3ovygL43sE52aY7jqcttvv30x11pzYg6jEKNAjh8/vsGymB+kvAP+iBEjFtg+2l6HaBcd7bLLxVwnJaVq9KaULyvfprXbt7QPaud46yjHW2X0tOMt+qj88Ic/LPYfw3bH0O5NcbxVTu7H3MYbb1yMJvnXv/41XXPNNenXv/51UQsctRjx3i644IKiyVtj7fmbXDrZpzaPt7lz56bzzjsvff3rX68f4bS1quk7ToCjIuKqSTTDac6yyy5bf7/xL398ObRUdd3StuXNKKL9cjTNaCymAChdNWn85VHeOba5SWxD+eSN5du0dvvyfcRVqvJmBtTO8dZRjrfK6EnH2/PPP5+OPvroYsCSuKJdPkdXY463yulJx1z5JNzRpG7VVVct+jlF36jGJ83t+ZvcUt8lev7xdvnllxfHVluaTlbjd5wAR0VEc5wHHnigmJy4KeWjDzWeyT7aEkdH+lJH1rZsG/07Spq78hJfLqVfttLoVyXxh6TknXfeaXL7xsvKR0oqv9/S9vFHqtR+u/w1qa3jraMcb5XRU463eA/RXPKQQw4pTqQXdhLieKucnnLMNdXcrjRZ9EMPPVSMeNiZf5OpreNt0qRJRS1v9L9rbpTLXL7jjEJJt4vmONFmOX7G1ZSYd6P8l7JxG/UVVlih/n5cUYkTipjHJZx11lkLDAfc3LalqyfxmtFJuqn21+XV3fHL1/gKU7TbvvXWWxfaEbp8WXn5ollBXDWKMra0fXw2pSuN8ZrU5vHWUY637tdTjrfbbrut6IMSc8KV5jRqPOhEHDPlE9k63ioj92Nu7NixxYAX0bQxLhQ0Hn6+fLCcGBQjRj4sif5SMYhGhIL2/E2mto63+++/v5gfLgY8aaw8jJ5++un1+4/+vqWLCNX0HacGjm4XTSHi4A7xS9jUFZzSBLHxC1TeeTmGYy394pdOIhor7Ts0bu4TX/QxR0lLV0/il670BbLKKqs0WBYjLJWu7MTIR82JCb9DNDsq/+MTV3xinqXWbh+a6yhMzz/eOsrx1v16wvEWE8/G+4jR3JoKb6Ur2aXR4Eocb5WR8zEXNTnRr7I0j9g///nPFt9r435w8fc4/i4v7JgrNcWLE/rStD7U3vF2xBFHpDFjxhR94BrfygdmOuGEE+qfL4W3avuOE+DodqVf7BC/GI3bMMcvb+nLNkYzKm+2U75tXKUbPnz4AvuPeUhKQ9I2daWt1ME1RjdqamjXV155pfgCaOqLPq4Q7rfffsX9Z599ttn3WJro8dBDD11g2YEHHlh8CcUfq+ZGJiptv9lmm6WNNtqo2dehZx9vncHx1r1yP97ipCiCWwwyEHNvNac0H1xjjrful/MxFye6pVAW/aKaOibLm8E1nroiRBPfMHny5GYHligdc/H+WjvnFz3veOsM1fIdJ8DR7VZfffXi5xe+8IViro7G7ZDHjRtX/PLFycNXvvKVJreNKzoxKtXAgQMbLJ8yZUox6lGIKv6mmmvEF0YMVRtXjuIqcmPRrjvE1ZjyzrQlX/ziF4tq/SeeeKLBiEjlf5DiyyuuHJU39SiJWrkY5ralDt3RTCTacccoSdT28dZRjrfulfvxFhPwxrZxEhUhrfEttr/kkkuKpkhNcbx1v5yPudLjgw46KN1yyy1piy22WGD70rEW/eEaTx4fogYuTtSjyVxTkybHSX70n4sg0dSohdTO8dYZquU7ToCj25WugDXunFq6+hvzwEQb55hMtnEn1WjOE9XRsW3jX+z4woiq8fgyOfzww5utto4rJ8cdd1zRRvovf/lLg2Vz5swp+n5Ep9OYg6Yp8Ufg1FNPLX55r7322gWWx9wkEfBixKzyCS7LHXnkkcWVmdi+cQiMk6T4UoqJckujKVG7x1tj5e+jqffUFMdb98n9eHv00UeLn9GvJcrY+DZ69Oh00003tfgZON66V87HXJQn5u6KWram/l4+8sgjRZO7WPajH/2o2dqz6FcVoSCO28a1MjFoRdT8nHTSSUWtDrV7vLW2ee7C/rZWw3dcr7o4C4VuFrVU55xzTtEcIjp4RiiKau84MYgvhhi2uqmmEiF+WS688MKi+jy+SOIKS7Sjjy/pqLqPNs6lNsoteeyxx9LZZ59ddIKO+UxivzfeeGNacskliy+HGGa2JVFFHttH9Xi8XsxZEu33o+11TDy6sP5McfUo2v7H5KX77rtv0Tn36aefLuYuiS+vppoWUJvHW2xXms8mrk6W2vdHWUojXMUfiihLcxxv3Sfn4y32H8dIa8RAAKX+R4053rpXzsdc1Jz96le/Kl7zM5/5TLFe/D2NWuA4QY95weKEuTQ/WHOirPEZxMl31A5FiJgwYUJRG/Ld7363CIp0jpyPt3Ixv2UEx4hC0QQ3AmCIoF+6WBDHZKnrTDV9xwlwVFQMIxt9yeILPH7Zok10c7/0jcUJbXw5x5W1+IWNq28RptoyNGx80cfVkmjyGF9A8foL+yNRLq70xZfICy+8UFxN2mCDDYoJSdsiTsjjCmN8BtEsIP7ImBepa+R+vHUGx1v3cbw53rpbzsfcW2+9VdRclAbFiPLHgBWNRyJszWcwceLEoizRZC8CRksjFlKbx1vu33ECHAAAQCb0gQMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQCQEOAAAgEwIcAFTQmmuumXr16rXArW/fvun73//+AuvPmjUrLbbYYvW3prZdZJFFGqyz+OKLp6WXXjqtt956ae+9905jxoxJb731VqvKd9ppp9XvJ/bb1Ou9+OKLbXrPP/7xjxfYR+/evetfJ8rZWrvsskv9dn369GnTfn/5y18u9L2Vf47lt379+qVBgwalnXfeOZ1++ulp6tSpqa2GDh3a5GvG7atf/Wqb9wfUBgEOgC7x/PPPp2HDhqWllloqfeMb30h1dXWVLlLV23///dM111yTHn/88TR79ux0wQUXLLDOwIED05w5c+pvN998c4Pl++67b4Plcfvggw+KwHbjjTem1VdfPR155JFprbXWKsJZLG/JiSee2GBfsf/G7rjjjja9z7FjxzZ4HAFz2rRp9a/x7LPPtnpf//jHP+q3mzlzZvH+Svr3759efvnlZvd79NFHt/je4rNt/FmWv9a//vWvIsBdcsklxed5xBFHpDfeeKPVZX/00UeL/TzwwAPpN7/5TRHoABZGgAOgVaKWJWoGWmv06NHFCep7771XnJw2PmlnQZdddlk6+OCD0yc/+cmihqc1Nt988waPN91006ImqrGo0dtkk03SRRddlG699dYihPzkJz9Jn/vc54pavdaI/cb+G/vrX/+aWmv69OlpwoQJabnllqt/bsUVV0yDBw9OHRUXCyJIlaywwgpFLVl731vjz7Zx6Nxggw3S8ccfn5566qnic/ztb39b/N/de++9rS5zBPKtt946jRw5Mp1wwgmt3g6oXQIcAK3y8ccft2n9xmGvLeGvVi2zzDLd8jp77LFHEbBD1P7svvvuad68eW3ax4ABA+rv33XXXenDDz9s1Xb//Oc/i8C20UYbpZ4iPovbb789bbHFFum1115Lu+66a/r73/9etf//QN4EOABa5d13323T+qecckrRhHLJJZcsmuxFUzOqx6hRo4omhuHBBx9Mv//979u0ffQ9K4nmnvfcc0+rtoua2OHDh/e4QB996C6//PKipjNqN7/4xS+mp59+utLFAnogAQ6AVvdpa4u11147Pfzww0Xwiz5CPe2EPXdLLLFE+vSnP13/+LrrrmvT9nvuuWe7+sFFgIsaqp4ommCW+tFFqD3ooIPS/PnzK10soIcR4ABolbvvvrvSRaCTlfcPe+GFF9q0bQwWEn3q2tIP7rnnnisGK4kauJ7q29/+dv39J554Il166aUVLQ/Q8whwACxU9G9qaw0N1a+831rUyLVVeS1chLOFhcAYMXKzzTYr+sD1VNtvv31afvnl6x+fffbZauGATiXAAbBQP/vZz9Kbb75Z6WLQyWKI/ZIYPbGtGjejXFgtXDSfLO871xNFU+Hy/p4xP1xbRukEWBgBDoBmvf322+mss86qH7GQnuP9998vmviVfOUrX2nzPj71qU81qE1rqR9cjHIZzXB7av+3cttss02Dx3/+858rVhag5xHgAGjSeuutl1ZZZZX005/+tMEUAosttliDW6xXLgYsiedjXq2ojSjdmvLDH/6wWLd3797166255poNmtzFEPcrr7xy/Wv96Ec/WmBEzAiaxxxzTLFtrBd9u770pS8VE2K31Z/+9KdimP0IJjEX20orrVTUNP3lL39JPcn1119fhLgwYsSIdo0SGv9vn//85xtMJ9DcxODjx48vmhJ+5jOfST3dkCFDGjyOCb8BOosAB0CTnn322eJkPG5rrLFG/fOl50q3WK/cN7/5zeL5mTNnFpMUL6xpZqw7ZcqUBZZFUIvwePLJJ6eXXnqpGNEyAtU555xTTJo8d+7c+iZqEQoibD7yyCNFmIsh8m+44Yaihqi1J88ffPBBEdT22WefYoLmqJ2KkQSvvvrqYr977bVX+trXvtbm+fCqUXxmpUmjDzzwwGIC6vYqb0YZn2Fzg91EGN9hhx3Soosumnq6ddZZp8HjF198sfhsADqDAAdAl1hqqaXS0ksv3ap1V1tttQaPY+S+mJsswkCEsJhbK0Y8jFqjqKWLMBfBLmqQooYuguAPfvCDYvCImNvsuOOOK4Zw/+ijj9JXv/rVhU4yHfuJmqTbbrutqMmL/UWtX7xu1Ez98Y9/LNaLeb6++93vphxFgHjmmWfSeeedl7bccsvivcXcb9dee20RjNsr+rSVh7Lm+nv15OkDGmt83Efob+oiBUB7CHAAVJUYLOUnP/lJuuaaaxaorVl//fWLWpxSyDv//PPTjjvuWDR5bCxqAkMMW3/LLbe0+JpRGxVhMUJn1Po1Fq+x4YYbFvd//etfp8ceeyxVq9NOO22BZq5xi2AbNYvx/o4++ug0ceLEdvV7aywmao/Pp6V+cLNmzUoPPfRQjx/ApPwzaaxxs1+A9hLgAKgqURv2xS9+Ma266qpNLt9oo42Kn9OnTy/CSKkpYHPrLWxwjQgyv/zlL4v7MQlzhLimlGqP6urq0umnn56q1YknnrhAM9e4xecaISpC2ymnnFJ8vtEHccaMGR1+zfJmlDHh++TJkxssj75xgwcPLgJ4LVhkkUUWeK7U5BegowQ4AKpOBKnmlM+xFc0ro+/bwtabNGlSs/uLGrUYITGU1yQ1Vj74RjQHjOaZOVl88cWLppMRVkvlP/fcc9Pmm29e9PHriMbTCTQOzNH/rVZq30o1jq2plQNoDwEOgKqz7bbbtmq97bbbrlXrNVfLFH2Tbr755vrH6667brP7WGGFFRo0hysfgj83EVTPOOOM4n70zYraxYVNwt2SGORm4403brYfXC31f2suwC233HIVKQvQ8whwAFSVAQMGFP21WiMGGmmNUg1bYzGC5uuvv17/eJlllml2H41PwF955ZWUsyOOOKK+VihG7oxpFzoywmZ5Ldw999xTP+pijMD4v//9rxg5tFbECKzl4nOOJqQAnUGAA6CqtHbkyrau25SXX365weNoTtjUACBxGzp0aDG3XdyWWGKJ4rmcxXsYNmxY/eNHH3202REk2xrgos9dafqGqH2LppstheOeJkb7bHxcNTcXIkBbCXAAVJW2nOh29KT4jTfeaPD4P//5T5MDgJRuUZMXt/feey8NHz485a7xoCK33357u/f16U9/ukEz01I/uM7s/xY1hFHzWe39DxtPIF9L/f+ArifAAVCzGjfVrLXJlhvXYEZTx/bq3bt3MZdeeYCLwBU1cZ3R/y32t9ZaaxVNEaM5a8zZtzCVqvVqHOD233//ipQD6JkEOABqVvlIlaU56GpZNA/tiPJmlDEoytVXX53mz59f1M51xJNPPlmMTDp16tT6QWSamq+vqZE3u1vU1MYk9CU77bRTzUyfAHQPAQ6ADvcjiwm1c7TpppsWNUdtrYGK8HDhhRemnjZa4uqrr96h/UVTwfLJ12OOvhi8pKPBMCZt//DDD9s8LP+gQYNSd4uawmhiWxJz7gF0JgEOgA73I8u15ipGvNxiiy0aDOTRGnfffXdaZ511Uk8bbCNqizoiJkHfYYcd6h9PmzatU/p/xeTg5WL+v1/96ldNrhsTlpcPHtLdLrroovr7Bx10UKunxABoLQEOgIXq169f/f25c+c2WDZ58uQGg1fk5pvf/Gb9/X/+85+prq6uxfVfeuml9Nhjj3U47FRaBJ3yCbxXXXXVtPfeey90u9Ln09zn1HhS787o/1aqbYv+bw8//HAaP3582mijjZpc97///W/9/f322y91pxgE5q677qqfG+/iiy/u1tcHaoMAB8BCxcloeWAr9+c//7nBcPS5+fKXv1w0pSw1oSwNf9+cE088sZhDrbVz1VWr3/3ud0VfstC3b9/0m9/8pkFQb05p/rvy+fPK7bHHHvX3hwwZktZcc80OlzVq3ELUZrV0rN13333FnHYhagK32Wab1F1ee+219I1vfKO4v+KKK6a///3vHZ7mAqApAhwAC7XXXnvV3z/11FOLgSnC9ddfn/7yl7+krbfeOuUq+mxdd9119bWI3/ve9xr0YSp39tlnp3HjxqXRo0ennD3wwAPp+OOPL+5HH8DLL7887bbbbgvdLqZQuPPOO+trK5sSI0WWasc6o/YtfPWrX00DBw5M99577wKTZJfEMRl97kIEpz/84Q+pu0SojfcaPzfccMPiGFlvvfW67fWB2iLAAbBQUeNUqvmIsBMnyDFAxIEHHljUOiyyyCILbDN79uwGJ9ulmpvmmiWWxDaxbWsm3248EXdz68VgHS3tM0YJjFCz2Wabpaeeeiptv/32RUCJaQVizrE4IY+apTFjxhTzmkXfuWoRzTkbD2FfCtjl4n3EPHfHHHNM0fwzPo+YYDuaI44YMaJVfR3jOCgN9HLFFVek3//+9wsMLlLejLK5/m/x/1He1DH2XX4MNBY1WldddVWaPn16UQt34403NpgLLvpgHnLIIUUNXATxW2+9tUGtcWvEZxafT0ufbbl4/WeffTadfvrpRWCN+xGKH3rooaLmEaDL1AFAK8ycObPuqKOOqltxxRXr+vXrV7fuuuvWnXbaaXXz5s1rsN7FF19cLO/du3d0kqq/xeN4/sQTT6xf9wc/+EHxXK9evZpc96KLLqp/7Xi8yCKLNFgvbosuumjdSiut1OZ9NmX+/Pl11157bd0Xv/jFutVXX71u8cUXrxswYEDdsGHD6s4555y6999/v9M/1zXWWKO+jK1R+ixKt8afR9z69u3bYJ3SZ7fKKqsU72XUqFF199xzT/F+F+bUU09t9nXi1qdPn2J5ufvvv794vXfffbf+ueHDh9eXpfGxEbf4/yotj2OrKRMmTKjbfvvti/UXW2yxunXWWadugw02KN5vbL/33nvXTZ06ta614lgovWbso6n31/hzLL8NGjSobtddd60799xz69544426jrrrrrvqX/crX/lKh/cH9Ey94p+ui4cAQEuij9iUKVOK+/4kt07MMXf//fcXtboxWXhMf7DjjjsWg5zkLEY3LQ2O85WvfKWo4QRorO8CzwAAVLG11167uAHUIn3gAAAAMiHAAQAAZEKAAwAAyIQABwBVYsaMGZUuAhU0derUShcByIBBTACgSsScenvvvXfaeOONi7npFltssUoXiS4Uc99NnDixmHswft52222VLhKQAQEOAKpETFAdt9CnT5909NFHpwsuuKDSxaKLbL311i1OYA7QFPPAAQAAZEIfOAAAgEwIcAAAAJkQ4AAAADIhwAEAAGRCgAMAAMiEAAcAAJAJAQ4AACATAhwAAEAmBDgAAIBMCHAAAACZEOAAAAAyIcABAABkQoADAADIhAAHAACQ8vD/AVAgLLufN+o9AAAAAElFTkSuQmCC", | |
| "text/plain": [ | |
| "<Figure size 432x432 with 1 Axes>" | |
| ] | |
| }, | |
| "metadata": { | |
| "image/png": { | |
| "height": 440, | |
| "width": 440 | |
| } | |
| }, | |
| "output_type": "display_data" | |
| } | |
| ], | |
| "source": [ | |
| "data.plot()" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 79, | |
| "id": "83e852be", | |
| "metadata": {}, | |
| "outputs": [], | |
| "source": [ | |
| "prior = tj.JokerPrior.default(\n", | |
| " P_min=1.0 * u.day,\n", | |
| " P_max=1000.0 * u.day,\n", | |
| " sigma_v=100 * u.km / u.s,\n", | |
| " sigma_K0=30 * u.km / u.s,\n", | |
| ")" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 80, | |
| "id": "bed7f4f2", | |
| "metadata": {}, | |
| "outputs": [], | |
| "source": [ | |
| "joker = tj.TheJoker(prior)\n", | |
| "prior_samples = prior.sample(1_000_000, rng=rng)" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 81, | |
| "id": "3977e65b", | |
| "metadata": {}, | |
| "outputs": [], | |
| "source": [ | |
| "samples = joker.iterative_rejection_sample(\n", | |
| " data, prior_samples, n_requested_samples=1024\n", | |
| ")" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 82, | |
| "id": "d5f5718f", | |
| "metadata": {}, | |
| "outputs": [ | |
| { | |
| "data": { | |
| "text/plain": [ | |
| "Text(0, 0.5, '$K$')" | |
| ] | |
| }, | |
| "execution_count": 82, | |
| "metadata": {}, | |
| "output_type": "execute_result" | |
| }, | |
| { | |
| "data": { | |
| "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3AAAAJQCAYAAADPMYZVAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAWJQAAFiUBSVIk8AAAn9ZJREFUeJzt3Qm4FNWd//+DICIuLCKKbGrELQQCSIJJFBSIBtzDaEBhfiYzv5ioGedvNMw46kTDuEw0iUtMZqLOgAmauTGoUX8R9w1BARFUhCgoi4psIqAIev/P54ynU7duVXV3dXV3VfX79Tz36Xtvd1efqq7qrm99z/meNs3Nzc0GAAAAAJB6O9W7AQAAAACA0hDAAQAAAEBGEMABAAAAQEYQwAEAAABARhDAAQAAAEBGEMABAAAAQEYQwAEAAABARhDAAQAAAEBGEMABAAAAQEYQwAEAAABARhDAAQAAAEBGEMABAAAAQEa0q3cDUH/777+/2bRpkznggAPq3RQAAAAg95YtW2b23HNPs3z58rKfSwYONnj76KOPAu/79NNPzbp16+xtMeU8dunSpfYH5W23RmhftV4vqeVWspw4z63GMcjx1xLHYG1eL4nlVroMjsF04hhsjGOw2sdf1o5BnXvrHDyWZjS8wYMHN3/xi19sfu+991r9vPbaa82jRo1qfuutt4ouR4/XY3VbjB6nH5S33RqhfdV6vaSWW8ly4jy3nOeU+liOv5Y4Bmvzekkst9JlcAymE8dgYxyD1T7+0ngMbtmyJfD8Wj8DBgyw599x0IUS1oYNG8z48eND77///vvNOeecU9M2AQAAAFnV1NRkpk2bFtqFsmvXrrGWSwAHq0uXLmb69Omt/r9+/Xpz7rnnmrFjx9alXQAAAEAWjRs3zowZMybwvpEjR5qddoo3mo0ADpZ2oG7duoXev+uuu9a0PQAAAECWdezY0f4EadcufhhGERMAAAAAyAgCOAAAAADICAI4AAAAAMgIAjhEUr/diRMnhvbfjftYZGe71bp91Xq9pJZbyXLiPJdjsPrSvt04BpNbBsdgOqV9u3EMJrOMah9/WdiXktJGcwnUuxGoryFDhtjbuXPn1uw1R48ebW9nzpxZs9cE8L84/oD64hgE6mt0Co7BSs6/ycABAAAAQEYQwAEAAABARhDAAQAAAEBGEMABAAAAQEYQwAEAAABARhDAAQAAAEBGMI0A6jKNAAAAANCohjCNAAAAAADkHwEcAAAAAGQEARwAAAAAZES7ejcAyXj//ffN9OnT7e3q1avNxo0bzbBhw8yECRNMp06d6t08AAnZsGGDmTFjhlmxYoXp3bu3OeWUU0yXLl3q3SwAAFAjBHA5sG3bNvPTn/7UnHPOOaZnz572f4sWLTIXX3yxefbZZ82NN95oOnfuXO9mAkggeLvwwgvNsmXLzJ577mm2bNlig7mzzz673k0DAAA1QhfKHHjqqadsJRsXvEn//v3NcccdZ9555x1z991317V9AJKhYE3Bm4oHK9u+ePFim4kDAACNgwxclSn7pS6MkyZNKvk5L7/8srnzzjvt1fa2bdva/x1zzDHmxBNPLPzttX79evPwww+bPfbYw4wYMaLwmAMOOMDevv7664mtD4D6UbCmzJuCN9m0aZPtRgkAABoHGbgqWrNmjXn00UfLvsI+efJkc9JJJ5mbbrrJ/OIXvzCXXXaZXY7+//HHH7d6zuDBg82bb75prr76avP0008X/r9161Z726NHjwTWBkC9KVg79NBD7UWhNm3a2Is0GgMHAAAaBwFcFah7k8agKeDavHlzyc9T8HXLLbeYM8880wwdOrTw/7322stccskl5qWXXrJBmn/u9YMOOsj8/Oc/N+eee25hUkCZNWuWvVp/xhlnJLRmAOpJwZoCOB3nyupfd911FDABAKDB0IUyQU1NTTZoUvcmFRdYu3Ztyc/V46+//nobnB1//PGt7t9nn33MoEGD7Hi3Rx55xIwaNarF/Ycddpj9cf74xz/a11dxk7333rvCNQOQBgrWKFgCAEBjI4BL0Lhx4+yPqHjIxIkTS36uCo188MEHpm/fvqEVIwcMGGDmzp1r7rjjDjsmLmg8nGj8nMbRKZunDBwAAACAfKALZUrcf//99lYBXBhXlGTVqlXmxRdfbHX/9u3bbVfK9u3bmyuvvJLgDQAAAMgZArgUUFnwdevW2d/33Xff0Md5i5G88MILLe5TcZMbbrjBjB071px22mkt7ps5c2bibQYAAABQewRwKbBw4cLC71ETbnvvU5EUr9/85jdm+PDhpl+/fi3+v3LlStudEwAAAED2MQYuBVz2TXbdddfQx3Xs2DHwOR9++KH505/+ZN5++21bRMX55JNPzLx581oVPAmydOlSM3r06JhrQJYPAAAAjWV0BefOOvf2J15KRQCXAm5SXmnXLvwt8d63cePGwu9vvfWWHf/23HPPVbGVAAAAAOqNAC4Fdtlll8Lvn376aejjvPd16NCh8PshhxxScQZMVwDIogEAAAClqeTc2Tt3c7kYA5cC3ol4lUkLo0IlpYyVAwAAAJBPZOBSoFevXoXfN23aFPo47329e/dOtA3K7kVNPK7xd94xeAAAAADCbd261f4E2bFjh9lpp3i5NAK4FBg0aJCdlFtFR7zFSfy891WSdg2yYcMGM378+ND7NSn5pEmTEn1NAAAAIK+amprMtGnTQqcR69q1a6zlEsClwG677WYGDhxoK0Zqku4wK1assLeK1ocNG5Z4N87p06eH3k/2DQAAACjduHHjzJgxYwLvGzlyJBm4rJswYYIN4FRSVN0Zg97QxYsXF0qWdu/ePdHX1+t169Yt0WUCAAAAjapjxBCkqMrzxVDEJCWUgVNWbcuWLS0m9nbUvXLOnDk2W3fWWWfVpY0AAAAA6osMXJVs27YtsHpklMmTJ5sLLrjA9pXt37+/HRfnPPDAA3ac2pQpU8y+++6beHspYgIAAACkv4hJm+bm5uYK2wbPQMVZs2bZ3xUMrV69ujBnm+ZZa9OmjU2XXnPNNaHLWL9+vbn22mtt0Dd27FgbNM2fP98u9wc/+IEZOnRo4u1WQRQVSImaDZ4iJkg7XeCYMWOGHSuqKq2nnHJKiyk6AAAAamnq1KmhRUxmz55ti5gsX7687OUSwKXUkiVLzKJFi2wg16dPH3PEEUeY9u3bV+W1FMApAxc1GSEZOKTd7bffbt58883C33379jVnn312XdsEAAAa92Lr1ogMnCtiokRNuehCmVIHH3yw/akVipgg61yV1rC/AQBAfSl4cxdbdau/83yxtSNFTAAgnH9y+6QnuwcAAJXhYmsyCOAA5IK6YajbpLLJutXfAAAgPbjYmgy6UMKiCiWyTn3o89wNAwCArNPFVf8YuDzbShVKVAtVKAEAAFAveS1uMpUqlKgWqlACAACgXoFZXitJb6UKJaqJKpQAAACoR9XJvBY36UgVSgAAAABZUWpgRnGT8hDAAQAAAEhcqYEZlaTLQxdKWFShBAAAQD2qTua1kvRWqlCiWqhCCQAAACSLKpSoGqpQpkdey+gCAIB84xymNapQoqqoQpmtak0AAABpwjlM7apQEsAhVRr96k1ey+gCAIB84xymdgjgkCqNfvVGQat3IkvK6AIAGlmjX9jNkjSfw2zI2X7ENAJIlUa/ekMZXQAAWl/Y1Vh9d2EX+TyHUZB1++23myuuuMLe6u+kzMjZfkQGDqmS5qs3tZDXMroAAMTR6Bd2G+kcppq9sFbkbD8igEOq5oErdb4QAACQf41+YbeRVDPI6l2n/aha88ARwMFSmnr8+PF1nweODFTt5a1fOAAgP7iw2ziqGWSdUqf9qKmpKXQeuGXLltl54OJgHjgwD1yDUz9z7wem+q0TRKcTwTYAIK/fVf7njRgxwjz++OOZ/s7byjxwqCbmgWtceesXnmeNXqUVAJDd76pigZ2/F5b3ArP/O29DCUFiGi56VmseOKpQAg3O30WB8QXpRbANAMjqd5WCqcWLF5u5c+eaqVOnmgsvvDCy0mTUd96MEqpK5q3ypBcBHNDgmLogOwi2AQD1VEqp/7DvKgVgCuDef/99oxFcGgMWFVRFfeetKOGCZp4vetKFEmhwFI7JDgbzAwBqzdsVcfny5Wbvvfc2u+66a2hX/rDvKv3+2GOPFR635557RgZVUd95vUsoeJLnCqYEcACQEQTbAIByVToWzDumTVmzdevWmUGDBtm/gwKwsO8qve5TTz1ll6Hg7dBDDy0pqPrwww/t81566SWbvdtnn31Mz549be2G9evXB17Q1Dpv2bLFLFiwoFAwJE8XPQngAKRioC8AAEhfASxvkKbAa9OmTYW//QFY1PmEbq+77rqSe5K4dr/yyit2rmI9R8M9unTpYr7+9a/bAPC8884Lfa6eM3DgQPv3brvtlqvzGgI4pGoib9QH1Q0BAMinSseCebsiKmh67733bCAVFIAVO58o1pNEAeBvf/tb88gjj9jxcgcccIANGHWOquzbzjvvbDZu3GgeeuihyEmw/eu4ZMkSO2av1heqmcgbDTGRN+ojzwN9AQBoZJWOBfOORVMAN3ny5NDgp9LzCb2O5iV2hU5ef/11m0DYtm1bIehp166dPW999913S15nPXaXXXap+YXqak3kTQAHSwfi9OnTQ+8n+5ZveR7oCwBAI6u0AFY5468rPZ9QG10Xzb322susWbPGjmVTMKdslTJwHTp0sG3SWLhS19kFgN7XqYVx48aZMWPGRE7kHQcBHCwm8m5sVDcEACCfalkAq9LzCT1H4+xcd8ndd9/d9OnTx2b+HnzwQdt9s2PHjqZHjx62kEmYqEnB3etkeSJvAjgAVDfMEArOAEA68flc+fmEtpkybhoD54qPKHjTtAUK2D7++GMbvCm4a+QL1W2alZNEQxsyZIi9nTt3br2bAqAI/1VETb5O8A0A9cfnc/W2qaYSuOeee8wnn3xis3IK7DQf3W233ZbZILmS828ycEglrmIBwSg4AwD5+HzO27lONdbHZc7uvfdem4XbvHmzfR0FdBoj16hVswngkEqUtQeCUXAGANIZBJX7+ZzUuU7YNoi7beI+rxrnbq5LptqisXC///3vbfCmYG7//fcv6yJmmvaVSsUrfQJUGVkGIJi+cNQtR4WHdJv1fvwAUAkXNGg+Wxc0ZOXzOalzHa2z5kxTV7ypU6eaCy+8sBCsxNk2cZ9XzXM3ldt/8sknzfbt203btm1tF8rly5cXDZK1HdQN84orrrDbRdspDftKpcjAIZXIMgDBKDgDAOm84Fvu53O55zphGST9rcBElRvd/GLucXG2TdznRa1PqdmvqMcp86aKjtu3bzdt2rSxgXKxINmbFdR2WbdunRk0aFBZ65VGBHBIpbxVCwIAAI19wdcfnIwYMcI8/vjjoec6/serOuPatWtbdVHUfY899ljheSrD754TZ9vEfV7UuVtU90rvei5ZssSsX7/eBmtaD63zeeedZ/93yCGHFILUNm3amJNOOqloF0hvkKbluTnmylmvNCKAQyqRZQAAAHm64OsPYhS8RZ3r+B+/YMECM3DgwFbBidb5qaeeshkmBSkqu++2RZxtE/d5UeduUVk973q+8MILdtJuTRWgYE3TCSiAcwGsMo2bNm0yBxxwQNGAV/d7g1FtF80jp8xd2veVYgjgYKk/sLuqU+5EhAAAAPWQpQu+5XZNLHa/yyBpG1x33XWBXQ/jbJtqbNOorJ5/Pbdt2xYaVGrsW++QLphBWT5vMKoAbvLkyTUtXLJ161b7E2THjh02mIyDAA6Fqxbjx48PvX/ixIlm0qRJNW0TAABAXpTbNdH/+JEjR9oAJigz5oIul4W68cYbU1VpMSqr511PZdbeeecd20VS2UStc7GgcsNn6zxt2jSzxx57FCb+1mvVO8Bvamqy7QqijKmKs8TBRN6wEwkqAzdz5szQx5CBAwAAiK/cMvZxyt6ncTLxYuvhvd8FNBrzVmqxk+XLl9tJvV0hl06dOtlCJWlY96gMnIJTZeDmz59f9nLJwMHSDtStW7d6NwMAACCXys0GxckepakqZ6ldG/2BmjJTV155pR0jqLF9l156qc3MFasueeihh9og7oMPPgidxqHWc8FFJUDatYsfhhHAAciFPE3QCQCobzXIUr8/avndU8prpbEqp9qrqpKuAEnnzp1DK2qKgrc33njD/q5b/X3bbbe1Wqa/uqS6TbpCJbpfy/S/t1GvmyV0oYTtQimaABLIqjR2GwEApP/74q233jJ9+vSJ9f1Ry++eUl6rlCAvTpfFUgW9vv7Wj5sCQF0c1fPLW1FTf1922WX292OPPdYWMlFW7aOPPrKB2bPPPtuiXbd7toWCQwVt+++/f6E7pZ4T9N76K3l6XzdL59/xSp8AQMqksdsIACB9/N8POumPuj8t3z2lvJbrdqmgRLdBwZjrfujqH+hHv7uMVCW8y/Z2l1SgpMIkCt6UJXMUfGkMmH4UlCkAVCDmgjfp0KFDq3adcsopNoDVcrU8VeHUOuu5LnjzvrfudZQF1K3+TkuGMg66UALIhTR2GwEApP/7Qif9/vtLzWol+d1TLHtW7muFLc8b+Hkntk4iAA0KMvWamnQ7qKLmvffea/8++OCDCwGfxrx985vftFk4dbf88pe/bB9XyhQJvX3baL/99rMB2yuvvGLHxuk1lZVTCf9x48Zldi44MnAAcsF7NS5s8HIe6AtZVymvuOKKwtVKAED87wsFDFHfH0FZpWp890S9TtRrhX0veJenzNOFF15oH6OslMtAafyYfqKCwnK+d/zPd3/7237mmWfaAEzVIvWjrJnapEBNZfe/+MUvmtNPP92cfPLJ5u2337brUEqW8BTf6wwYMMD+X+Pe1G4Fb0uXLrXbIMtj5RkDB8bAARnCWD8AeZClwlMKXBQ8VHvclF5HRTa8xT5UvCPu1AHedisLpWUOHz68xZgx3f/nP//Zzr3Wq1cvO3+cgidXCVKBzvbt283gwYML7fB+71RaEMbbdldOXwGd2rhq1arC9lBbFYwp0NupjO3vtsHPfvYzuyxX/XGXXXYxP//5z+v6/ckYOABoEIz1A5AHxbJNaRKWVUq6p4SW6+YyU35F26aU7RL2veBtp4I3l2lTEKSASEGQMlIKZhSUtW3b1txwww0tKkGqDVrek08+Gfh6/vdRwVuxMXhhGTM3ns21UQGXCo6oe6WrYulfr2Lb3hVpUdCmMXjuR39n+fuTMXAAkCGM9QOQB1m6GBU0Z1m5856VkunRcpuammw2TDR/2JIlS2J/L3jbrXnUVJ3R/xh1J1yzZo0db6ZA7pNPPmlV2EX/37hxY6vnBr1vaq8Cp1KrYHof488kOm5+N1XKVMC5bds2+9igCcHVTVRtUPCn8W4KDhUEuoBQwZt+19i4LH9/EsABQIaUeyIBAGmUpYtRpU6o7QITjeHaY489bOChYKHU4FSvs88++xS6+m3dutW8++67RZ+nLofqarly5UrbDXLKlCktgqjzzz/fPi7su0MBkfdWlKFzc7F169atkCFzz9W6/ud//qe59dZbbXsVEB1zzDG2UIgCvqjgNSzA9X6/KXOmAPaJJ56wmUM3v5ubEuBN37Jd8PbYY4/Z7aZun2qvxr4paFNRFN23evVqG9hNnDgx09+fBHAAkCGlnkigMcbmAFmVxYtRxT4bXGCi4E3dIJU1UnBVTnCqAE4l9F2XR/1djLo9qvujsk0KXP7u7/7O9OzZ0z5fwZULdIK+O4444gizefNmm2HTmDv9LSrs4sbAHXLIIfZvZfEcBYjTp08vjLFTYDRv3jzzta99rcXyg4LXsOyr9/tNy1cwqrZpWyh4c9viw8+6U6pbp/7//PPPm6efftoGj5pyQM9RJlHbRF0lZ82aZfr372/22msvuw6HH364ufjii02WEcABQIYQXFRf3O5PAPJ9MarYZ4MLRFyXPwUU5Vam1Hgvl8ESPb8Yb3dHZZwU1Cgj5oJIZZzCqDCI5mbzv54CHWX1wmhdFfQpSNp9991twKTX1uspqHPZx6DgtZTsq5av5ysAFpf90/MWL15sX0uZNhUn0a2oDQruNEWAaByh/lZAt++++5qdd97ZbpNSspppRxETAMiQLA38z6osjc0BkJ7PBheIuMBD3fRKKeThVWxagqACKd557NQNUhlAR9mrqAxg3GkQtExl7JTxUzCloEhj6VyWTOPQwpZXymsGFY5xz/vggw9s0KYulgrQtM4K3hSgKWhr3769/V3ZON0qyNR7pe9MBXDegDWryMABQIYQXFRflsbmAKheLwb/YzUuS8FK2GdDVLfQUl+3WGYyKAvo7e6oZauLoF5HY9gU5KgbpV4/zuuFbQ8FRxqPpmyYC570v9mzZ9u52xSceZer5/72t781jzzySGEib43NC9v2QdvS29apU6eajz/+uFCtU21Q4KZsowJavU8KLl2XShcsioK4rCOAg6Wd3/uh5KcPAP0AqC+Ci+rL4tgcAMl2kXZFMTQfmhtLps8DBQFhnw1hwVDQsuJ2zVZm65VXXimMkVP2Sctx3R1dkKVsmMauqSujzu/Kfb2ggNO77Z555hl7q66JCqI++ugjG8S5apX+7yY9d+bMmYXgSb8r2AprU9S23LJli12OJvhWFlCv6zKP5557rvmbv/kbG9C67a3soO5XMZNSxxUmRZlC18XTT109FVjGQQCHwgExfvz40PvVDWDSpEk1bRPAeK/WCC6qL4tjcwAk24tBn7MKABSceMeSnXfeeWW/Ztiy4lAw4oKgoPFc7vNL6+WdfLzc3hpBga53GQogRQGUvqtdBk4X+4O6Req57jnu+d7lBWXozjzzzFbf+WrH2rVrzejRo80999xj/6eJxxWoapu64iTXXXddoc3KTGoKBXVtLXVcYVI0LYSqkgbRPuHmqSsXARwsHSCqJhSG7BvqgWISrRFcAED1ezHoxF/ZGhcsFRtLFiXJZQVVqQy62Flpb42gQNdNNL5w4cLCpNrK8imTpJ8jjzyyVbVK7+t7t4F+1/9c2++9914baOl8U8FgWIbOtWvXXXe1r63t4AqdeNfRfVcGBYa1vPA5btw4M2bMmMD71BYycKiIdiDN84HGk+YsF+O9svV+AUBeejHoPnXVU6CiIEFBSdwT/ySXFVSlMuhiZ6W9NYICQC1D49bmz59vuySqRP+qVavMUUcdZTNeUd9Feq62gT+Qcm1X10t9v6m7YY8ePVpl6ET3K8hb9lnXSI11U0ET7/x0fi5jp8m8RUFhLb8zo4YgafvFRQAHNLg0Z7kY7xX9fulkQOMq9CVGMAcAyfVicMGFTvijPl9LuagWtCzxTrZd6ud3UGB24403tniM7vOuZ5wLf2FFRBSEKYB0QaQqPOo7qNjydL+6n/q7oLogTQGZumBqLJ2KjyiQU6CobeReW+1RV8h1n2UgFbz5A0f/umrMoDfgzcuFYAI4oMGlOcvFeK/WvO+Pu5qrSmBpC74BoBECvbgXQeM+T+3yfjfqtlh1zDivFbX+CoiUgQt7vTBRXT0VBC5dutQ+Rj+aW06vo8craFTg558bbtu2bZGFVnSrMYL6jiy3rWlHAAc0uDRnuRjvFf1+uTEQaQy+AaARlHIRtFhBkLDnSdAYLnEBm5anITBR1TGTvFCr11fgpCkKZMiQISVfXI3q6rlgwQKbzevevbvtlqmxcup6qKyc1l0BnP985d133y1k1/R/t53UHVPfjSpsojGCUdsmqwjggAZHliu775fGUag7SRqDbwBoBKVcBA0rCBL2PG+mSmO+3nnnnUIpehX30JgvN6ZL1q9fH1kdM+6FWtcOdUNUsKRgSBdW9bPffvvZrJmKlkR1n/Sui7pEagyfqwbp7eqp31999VXbfVJUFEVdJbWuem11pRwxYoR5/PHHC9vPmwUUBW+qvOmt9qnvzDxeCCaAAxocWa7sGjx4cOHLm+AbQNLyWjQpyfUq5SJoUAClYMRNvq1A6G//9m8DM1Uq2KH5zjTfmut5odL9/uVX2sYgrh2ad04BkQIqUZClOdfUNq1DVAET77oouFJQ5bpAquunGweo7aDJuEVZN03SrR/RxUotQ8Gb93zl9ttvb7FdRVk3N7xAwV/QZOregFQBZRb36zbNClPR0JT+lrlz59a7KQCK8H9hqWtIJQF4Xk/QAKTv86ZR16vYpNgffvihee+992xAocBCY8EU3Lhg5LXXXisEcJ06dbJzoKkgSrU/t6+44gobdD3xxBM2q6WJsMUbUOp/mic4bPu5Zbj1VPCkAM5V5nRdQXXfk08+abeVAjll6VTQ5LDDDrPbQX8rILvssstCt6t3eaIumNp+LkjWclxxFAWk2pZqS73260rOv8nAAUADF51JcxVSAPWV5iJXWVqvoJ4uQQWplOFSYLF582Yb9Dz00EM2UFPGTbft27cPneC6kgt2Yf93mUM3f5sbc60AztH/orafmzvOO4WCpiLQ8v/pn/7JZvfceO5hw4aZAQMGhAZkxTKNY8eObdHFUl0qV65cae/TmL0XX3zRnHTSSYUJxd1tFvdrAjgAaOCiM3k9QQOQ7yJXtVyvavRUCCpI5QIKBWsKXlSMQxUUlTnSj5uYOm5bwi7Yhf3fZQo11kyvp9cWdZ9UoREVHVGWS9nCMFqGprtx66hx2275yja6ib11qx9voOtdV72GtomCvnc/6/6oWzcO3E0GrgDNBYjTpk0rLEvdPhV4Kpuo7ap2u/mPs7hfE8ABQAMXncnrCRqAyuW1yFW565VUTwV/QKIAQmOYXUEqZakUxOy1115m5513tsGbGy/mLq5V0pawC3Zh//dmDr3dTufMmWOLjCgzqHZG0TLUfdGtn56rbpTa5grC3JxuCu7UpdE/N54LWhUEasydC77WffY8Ny5P201t8gZymhBcgaYok6c26HUUzKnr5+GHH267T2ZxvyaAA5BqjNGqbtGZvJ6gAahcXotclbteSfVU8AZfCigUPKh6pPuecxUX/dklb8EPZZWUAXPjwoLaUqxLpHfcncaoaYyYXstVh9Rr3XTTTS2mLlAg5Er267nKErqsoR4X1a1Tr6v2uGybuofqbxUQ8U6y/dZbbwUGp5oe4Nlnn7Vj4tS9tHPnzjYAU0ZObdcyFPRqfFu7du3sfXq+KnXqPj1G7T3qqKNsgRQFkpoEPKsFTIQADkCqMUaruvJ6ggYgP+p9IS+pngqlZLrC1lnBlqYQWL16tS0ooizSF77whUIQ5t0u3u9NBSvKXikL5s366XkK2hRM6VZ/6zFu7JleywVc+l0BkpsQW4Gbgi0FS96ALOy7RO1SVkyBlZufTeulro7edfRPC+C6bd52222FsXAqcLJx40YbxCoT5wJKBZIK8PR4bQP9TwGxnisKfrUtNK7QPUZ/Z/WcYqd6NwAAojBGCwAamwtIFCi4C3m1pABE2TIFBJV0ufMHflGBoAvqFOR4AxB1M9R2UCVFBVaaI04Bl3e7+AukqOuhnqMgSJkoVXJUsOYybrpV1k/tUfdGBT1avsaMKVBUcOQmxNY2UBVMzQOngEztcQFZ1LqoS+Pw4cPtYxcuXGjb+u1vf9u2XeuodVVGzL999Dhl1tRGVwlzt912M7vvvrsNRr/+9a+bk08+2Xbl1Jg8ZdwU2Ok1lHlTcOvmkFOQquDPPUbbJqvnFGTgAKQaY7QAoLHV6kJeUKZPvP9TIBA3Gxiny7oLXhW8KEOlAERBlAIn/Si4cXOruYyVAhcFbcp2aYyYsmf+bef/blWXTQU2qgqpoEpzsOl1FPQp6FJw5c1UKYjyPt87p1vQdnHrrkycHqMMmYJEtVXUlTRo+9x4443mwAMPNK+//roN3hRQ9u/f3waU3m6fCi7VRlftUoGtMpQu6Fd1SgWtqnLpMot6XFbPKcjAAUi1pK58AgCyqZzMVdKZPv0oKNBcXVOnTrVznunvONlABTT6DlP7FaToua6yYxgXcKnQiZvgWoGM/nZl/V31SpexcsU69H89VnO23XPPPea///u/bdZOwZ37blVQqO6QmntOAZUCPmW21EVRr+cybv7vXv93s0RlSV1GUYGmgjAVHBG9njJ+ypTpOVquyzoqeFOA169fPxtEKmumdenfv3+h26d7fY3TUzCn5SvTp/+54M5tR20fLUeBr9uGWT2nIAMHINXch767MqoPdIqZAEDjqFWxpbBMn6sOKZpXzGW8vI8pdZxesXHd/uUos6UsmDJHmgdOAYuW6wqcqC3eC5z6jnSBjChAe/rpp+1ylWFTgHXppZfagEev4YqluDnolH1TgKOgSd0StWxl2/z84/YUgDnKrrlMm3++OQVkbjyb2qlgTllEf+DntpECyXnz5tm2KIBTIZJdPwvMlFFzE3tr2Sp24gqvqM1qh3usdwJ13Zf18wgCuJx6/vnnzauvvmqvFAF5QDETAGhMtSq2FNZl/7HHHiv8TxUQXcbL+5hSv6OKdQf1L0fZMAVnepwCDn+3Tvc/F4j410HdCmfPnm2DHUfZxIMOOsj+rmycgjdlphQMqkiIuiwqsHPZPmUen3rqKXPddddFVpr0Fk4Rb1Dmqkmqa6PGqilQ1K0yhXquy/4pmFOApcBVbVLAp2IpRx99dKHLZRfPujruf3pNV9DErb83WMvLeQMBXA6pQtEvf/lLc8wxx9S7KUBiKGYCpLNCH5CXfTos06fgxY0pU2l6laBXZsr7mFK/o4qN6/Y/T+X7/cU9ogKRoHVQ+9944w17v8t4aUJrrY8yVAp4XNZOwaKWrYyaAj2XedT6u0As6P3xvq62jbfNbp2UHdNrK3hUsKbg8Mgjj7T3K+jT/xSsde/e3QZ3bjyb6yqqgE5FVnbybXtHy3ftVfZN7fR2V83TZyMBXI7cddddZsGCBXZH1RUOIE8oZgIEIzuNRtynXRChE3o3d1ql83qFBUbKPBULKEv9jirWHTSsuEjUtghbB7eN1CVSmTVR8KPukRqDpmBHgZvGgvmDIv3uzTwqiCo2mXjQpN9B20LdITWWTW3VtlV3TjeFgO5T8KU2K1B2k5yL2qpqlmeH7AsK+FTsxU2Eriqa+t9LL71kfve739nupd/61rfsc2bNmmVv9dpRc9ilFQFcDWiH0Y5YTnfGl19+2dx55512p9RVCFFG7cQTTyz87XfGGWfYHwVxP/zhDxNrP2qLq+nBmHAaCEZ2Go24T7sgQlUTFYgoe6NAJyrAifv9WkrXu1K/o4oty78cF9goqFEgoov1Gl9WSuDhtpHOQfV41xVT3RdddkuFSoK6RrrMncs8KvvlArFi70/YtlAbVETFZdX0t15X7dLfCri8c8Gp3YMHD7a/a+66sO0647P1VLCnapXaFw455BDbFdNNo6Aulfq54447CpU8tb/o9fS4rF30IoCrsjVr1phHH33UnHrqqSU/RzvirbfeagdmDh061P5PO+OPf/xj88wzz5gpU6bYgazIJ66mR3/hqjoVAS3wV2Sn0Yj7tAsa3Hg0dxt1AaOa369Jja/yL8dlsxSIqFLk5s2bW5XfD/vO1OOUlXSFPNz3qJbnipwoeAr6TnXZsaBArNj7E7YtFHAqWAoL7BRgaRiQ3ktVxtS4NxU8URv96+ldXwW0Csq0np/73OfscrRsBap/+MMfCuPhlJVTMKdupJpHTllJdS3N4kUvphGoEqVwFy1aZCZPnmwPtlKpUtAtt9xid3IXvInSwZdccom9+nL11Vfb5SOfuJqerglcgbRjqg34T2h10u8mMC5Wpj6r+7QLGtz4KHcbdQEji9+vbluoO6Er/KHzSgUiKr/vf3+935n6cQVFxAVNpX5euEBMCQXdukAv6c+csWPH2oBK7XbBlSpOKkh1FS3D9ucZM2bY9dR5scbXiQq2uPe2Z8+eNsuobJseox+Ns1MXUi3bbZesIQOXsKamJtuv1pVidaVSS6HHX3/99XbnOv7441vdr/7dumKilLYGao4aNSrh1iMNuJqe/S9coJbyVFkNlctDL45yuiyqC5wbA1csmMji96t3W1xzzTX2vFLBhwI5BSz+99f7HRlU9COJz4uwZRTrourfN1WVUhk5BWh/+ctf7HupYEzrpQIn6m0WVtHS0Wu5Cpr+Cby1rYYNG2Z/dN6sx2gfUeDmMrZDhgzJ5EUvAriEjRs3zv6IdqKJEyeW/Ny7777bXmHRzqUytUE0D4iqAqkPr8bEhY2HQ3Yx1iv7X7gAUC+NctErTiCS9PdrNQqphNEyH374YfPnP//ZZpB0nqgxX/731/udqS6F6p7oui2WWokx7ljBYhcPXFsVQCmY0nbT9tIwIY1/07QFWjeX1HDnwnqMe86TTz5Z2B5qk79LqOpAeCfwdlU8db8er3nllNlTAKftFzU1QpoRwKXI/fffb2/djPZBtLO5HfLFF1+0Vw6QL1xNb4mAFgDye9GrWLCQZGGvpL9fXcCioEHFM5Q1UvEMBR8asxXW9jjrpPtvuukmc+GFFwYWFgn7zlRbys3IlpPFVVuuvPJKGxRpHJom2nbr4g8uNbebxrq99tprtreZsmSux5ooU+gCuB49etguoq4tGs+nx6qwyZueNvnXV/PmeXu/eat4qprle++9Z8+ds14kjgAuJdxEiqLSqmG0QzsvvPBCbgI4Ki8iDAEtAOT3olexYCHNXUJdgKJCGK56ooIMdddTABfWdvd/ZZX0u7oQqjx+sXOfqMIi3sd4t4/GjgW12a9Y8ZMwCt7cHHMKxJQhO/nkkyMvHqiLpLJtCuIUtKnNOr9V4KUATROJjxkzxv7f7c9aru5T0Optk399dT79L//yL2bevHm2/QrovFk3ja1TZcuszw1HAJcSCxcuLPwe1n3Sf5+KpORFmj+gUX0E8NHYPqgX9r3sqcVFryT3i2JdPtPYJdStvwIdZZHcHGsqllFK292tugQq4FPmznvuE7V9y31/9Xy9jhsjpp5cWr7//fKeh2mddF6qqo1RzxFXDVMULKncf9hE27pPXRlVYVLBmgI+JS9UJfLcc8+1wZV/ugDv+paSWX788cdtkRd1ZxVl3PSj4E8FUpQp1Y/LlGb1fJMALiVc9k28fXf9tHMHPadSKk07evTo2M9XSrwSafyARu0QwEdj+6Be2PdQ7f2iWJfPNHYJdeuvLJUCI52bKWjRWDM351pU293/XSENV0HTO1G2C7o0mbaK18Udq6VASN0u3dxr6kYY9H75i5/ovE7rEvUcV/HRZeAU8B155JG2amUQt97aThrzpm2moE8/eo2wCbrVNnW/1HZ2E3B369atUJHSG+xqXJ3brqLuk3pNde9UoK3gVAGjAmdtX7WlEpWcO+vcW5Oqx0EAlxLakRzX/zeI9z7tjF733Xef3XHdQagUvq466IqQrmykWRo/oFE7BPDR2D6oF/Y9VHu/KNblM41dQt366oK7MkqHH364Deb8bQxru25VgVHnfupOqMBE3Sm93QNdds51C4wbJCvoU5DVp08f+7cry+/P7vnPw0QZK2XLlDDQ+WVQpvXSSy8tjIHT6+jvMG69tc6qpKllaZos/URN1i7K2KktAwcOtOvwwAMP2DFyomybivzpsRrzpqDTbTt1bdXQJAWjCjRVLFAFALV9VeVSrx2WXUwzAriUcAMsRVcHwnjv86fqTzzxxNivrysAlWbRKpHGD2ikJ4Bv9G5cXOBAvbDvodr7RbEuge5+9z1w44031v17wL/+Ct6C1iFs3fR/ZX6UvXFdG3XBXXMHu+Ur8+YoIKkkSPa2N6wsv/c8TG1R0KOuhgqAFDAp0AkKItW90gVSXnq/FKwpmSDKSrqJvLXeDz30kH2MgrmgYiziX2cFiQpEFy9eXGibaPutXLnSBpxaLwV5LrOm11CgrYBO66HAUT/KwimrF5VdLEUl586V1LFgIu+U8H4IaQcMo0GfpYyVy5qwySLRGIpNCtrok3nnaaLmPEwy3EjytO8h2/tFmr4Hklh/BScugzd8+HCbvfJOlK3ASGPjXOEOF+DE+Qz1tlc/LtPn2uE/D1NbNJeaXt/Rc8oJIvX+KLhR+9asWWNuvfVW8+1vf9tm/+TrX/+6DbT0ulp21GTtjtrlArZt27bZ5Id+9BrKUmr6Lt2nsW1uWi89RwkPbUetuwqaqKunsp7a/vrJYs8CMnAp0atXr8Lv3r67ft77krwS6iY8DKMd3Tv+DkhSsSuwjd6NK0+VOBlTlS152veQ7f1CXfheeeWVwlgul30pp5dGUr05klj/qCxmVLXJOJ+h3vYq6CuWPXVtU6EPV7pfgY4yVnq+a9OIESNs0ZCg7an/uXNWnV/q/VLiQe+dMmgKXPXjAuCg98X//7/927+1r9e5c2d7vzsvVeCm81h1jVRwprFlChTVdgV4rouqAjhVuHTtcXMpV7NngTKM+gmiTKDaFAcBXEpoJ9aOpAGWUcVJvPclOYWAPtTGjx8fer8mJJ80aVJirweUg25c+dHowTiAeDS2yY1r0q3+dkoNampxASkoSHSv7f1fsaEjYUFipZ+hpQxZcY9RkKPHaZ001526Jfbs2dMGc9p+Gvvmxtb5t6eW7caiuWyZ/lZGTEVZNBZNvysoCyvaIsqmzZ492z6mqanJfOMb3zDXX3+9+Z//+R9z55132sBMy1W2UBUs9Xj9rnbqtd15s9qp11Mwp2BPbdLzVAylmhlktXnatGmB9ylrqKA4DgK4lNBBoVSy5q3QJN1h3IGqiH3YsGGJvb4+KKZPnx56P9k31BNjJPODYBxAHCpUoZNxl4FzZeJLCWpcUKUTaXWhU5e9anWdCwoS3e/+QMcFO+VkBkv5DI0zDUFY4KkKlsogqYiexqAp8FLSwTsmzfFuTz1fwZTGwOk905g6bfeXXnrJzv+mIUGahkHz5amYiAIqvY6GET377LO2u6XaqmqVbp49LWf16tV2mbpPY+n0Ps6ZM8c+R4GZ6yb6l7/8xWb39HyNd1Nbjj76aHuerX1n9913t+fe+ltj86pF3Tg1p10QjQskA5cDEyZMsDuSUr9uJ/Rzg0+103bv3j2x13alXIE0ohtXfhCMA4hDhUK8Bd90cl5qUOOCKgVvrny8gpAkLyCVEyT6/1dOZtB9hqpLqctKKthRUKJt5O4vJ9OotitQU0ZI48Wee+452wVR3716HQVveh0FQeqi6B+TFtYNVMGZfrzBodqspIAKpSgo05xtCqZUYVLLVmCngE73q9q6AjDXXVZ/q9th586d7e8K6PU+arze888/b/+vfUS92dRNUsvT+uj/rtul/qf3x1V1r6TCZymihiBFVZ0vhgAuRZSBU1ZNB44mUNTfXtohdZVBVwzOOuusurUT1dfoVReRXwTjAJK++FPswpALmBRUKXhTFinp4ivFgsSoALOcbpHuM1Rj0RSsKIvlugrqb7cdvBSE6fEu6PMGe1qenqNARpkxd7vffvvZIElBm8bAic4/VVzFTdTtxqS5edr0WBVXicr6qb06x3FBmQI4JSQUmOk8V7d6HS1L2T0FcFov93gFZHt+Nm+eG2OnQFnvrRvTpiye2qiulHqsMnuqUnnPPffY13DbR49zQWrWEMBVkXeArbd6ZBSVkL3gggvsFZz+/fsXdkbRnBfa6adMmWJT0UmiiEm6UOgBAIDSLv4UuzDkMnSu6qOCt3K/U4tdWC0WJHpL6rtJqN3z/RlEBSnKqnnnVlPg5OVezwUx7ta1z7s8BW0KglQExh/saTvoOW68mrobiguSFORouymw03ng4MGDbZdDtd1bAVM9yFSS342PCztvUbdBVafUOipQVEZMP5qmQIGW1kOvqfNf3Sobp0DMjaHT2LZDP6uiqSydCyY1vcQNN9xgt5mKr9x88802WNQYNHW7VKCo90TLU1VMBXI6l9Y6ecdTJo0iJhmhHcXNEu8NiBT168DRTqiU6TXXXBP4fF110H3XXnutueiii8zYsWPtzqWdUMvVlQ0dPEmjiEm6UOgBAID0dN0udmG1WJDoah2480Pv8/3tU6CnjJFo/JeKhfjnWnOv5wIvF3C59fMuzyUUFLgp6HHr4f6vxyjjpcBTgZXLaIker3NEBUUKJhV8qX0nnXSSfY4713Vl/BWMKQhTl8WgKpVuLjiXDdSPzov1PwVs6k6pwEbLU2ZORUvUnqBqlZMnTy4Ekvq/2nfUUUcVgusvfvGL9rFKUjzxxBM2kFMwpYyiG9enYNo7njIrRUzaNCv8RSpp5160aJHN3imNfMQRR9grB0lTNUvt3FGTEZKBqy1/md84VwsBAMiKtA8d0AV0nSs5ypxozrRS26/nuyDJBTgKyoLW8dhjj231Wo8++miLx7jXC+sWGXROoWSCnqdsVo8ePWxJfbXB23YXUKiio1sPZbfUHtddU8kIzV2nypQuKNV9r732WqGHmIIkZdW8RU6CzmXca6sdaqPCEgVzCrDcVAPFtvfy5csL2T//62jdtc01cbgb46cgURk+dRPVelTzHCsqA+eKmGjblYsMXIrpQNRPLVDEJF0o9AAAaCRBGa6w+cHqoVihlFK6cWpd3FQICojCuhkqk6TMmzJSynApoFEgUko1ST8FOgocFWzpVpk6/SjAcpmnoDnXXGVGdf3U/xWEqNuhsmJB54tuom837kx/q7BIWJXKoO02depUG8CJtlPU3MczPPuLglgFQW7snAI0UQZQfytB4Yql6HxXAaxulTFU8RSN5asWipgADYRCD8jiFWoASHLoQJrGg1d6YVWPV2VHb4ATNjxCY97UbVJDZ3Tyr26Bcddfz1EQqEyZy+q5rJa3kmfY9AcKfpTZUgDngiDXvVIZJAVNWg/9T9vFO3woqkpl0PbR/G8KxhS8abyagitNrTV06NBWQdYKz7ZTkKbvR7XRdQv1zlOnYFXbXNtSXTz1GgcddFChUqi6efrHGKYdARwAZESaTmYAoNoZrrjjwatxsavcC6tBbdC4sVLmwVQwoS6F/m6bccbDe5+jgEXBi4IcdbvUrcvshW1rZcEU9KhrpDJv6qqooC1ouwatsyss4ibtDqPlaPJuTWegOd1UTEUBlwJeBY7+IKtr1662uJ8ylSpK4ibo1q3aJ5qWS+P+VMlSmT21fdSoUTb75q30nsU6AwRwaMgqlGQykEXeL5kPP/zQXs1lH64ePieA+ma41H1PGSA3ebfmwC2FlqNxT/p57LHHbGZHwYG7rxbHdFSXUDduzRtABbWjlEm7i/EuQ9kmZc1UMVLjwRTcuK6cQa+lzJamEnDj5j73uc/ZILTUaqBaN2XA9KPvLGXEFMh5M5j+90P3q23Kpun1lYlTEBdU6v/tzybpVn0IjW3TFA7KtrlpD0TZPI3F0zmulqmsowJR7zlvkvMB+lGFElXVaFUoyWQgi7xfsDoxEX35sg9XB58TQDaHDigg0GekG2/mJmuWWh3TQRktt4433XSTrUyudulWgYomvC6n22bUBSbvfQpslIlSsKPgSBfj3Vxvbr46ZazOP//8VsGlujAq2FHwJW6MWdBcb8W2gSveomDO20XT/35ouQrWVcBPY+7cODF/qf9Vq1YV2qVunQrydKvCLBrbp0BNxf+UodPrKvt4+OGH2+0ftO2yVoWSAA6WDsDp06eH3p+n7JtQph9Z5L16q33WTTESNZYC8fE58VdkI1EPqoToxmu5v0uhfVSZN0cBQdDxW84xXe4xEJU9Uxl+F1zqVn8HBXBRQW3UBSbvfc8884y9Pfroo+2txtW5qQdEwY3a5p8gXFR2X4/91re+Zf9W18Zbb73VZu5E5426+B+WTfNuA72OMoD6znJVOPXdpddSIKYA78knnzSnnXaa+epXv2peffXVQlCmrpP+Uv/vvvtuoeCJMllqv+bO82cAlY1zy9ePy3jW6mLcuHHjzJgxYyKrUMYR71nIHVeFMuwnbwGcP11ezfQ5qk9frPpQ1lVB3XonF80T9wWr6rTaZ3Vl0l1BZR9OHp8Tf+VOCL0ZXyCtx6BO0HXSrwBAJ/CuwEYlx3S5x4DaoCyQzq+8E3rX4gKT93cFS95qjuKKd6gKo75DXBEQ/3MVvHmfq2kCXFVL/X/hwoV2W6irq8auKbj74x//aG/1t6pAum2g98Nl/XSr7eiyai5bqi6Q6tqoc86vfe1rhayZvvP0WO93/D777GO7dGp8nrKFapN3Pbzvget+qeXU+vNL6xJ2bq3vcLpQAmWgTH++NFpXN1fxy3VJ0RcA+3Dy+Jz4K7KRyNIx6ApiBD037jFd7jEQlT1T5uX++++3XeikV69eNvAoNautx6owiJ7vKlq6ypDizXy5Cb69r+0m0lZ3Q/nDH/5gx8Vpm3mfq2VqIm99x+j/ypa5ro2iaQ5EGUSNldNjFUxpzJcCKxUecdtAbVaGTGPWXJtdcKvMmwJKLU8TbqsoiSho03KUCXRz47nv+IMPPtj+ra6hep6rhun9/nfvgZvMu9T3LgsI4NCQKNOfL412cum+YL2loOnOlry8fU5U0g0yiWIKQC2PwbDnljp3WlR3QNHfcY8pzbH23HPPFQI4dQ1VFsvfjTJs+fqfKiqqK6Mu4ilwmjx5cmDg6wq/eCfm1jJ039y5c1uME1QbRHPGuWBPy3Xr9PDDD9vnKPBSIKeATl0i1YVRQZlbHwVx/m6rWoa/CqfLTKrIjLpnKjOnzNTKlStt4Obmc9PfbpJuWbFihR2zpyyfAkeN81Omy43n8/O+d+pOqe1V6ji+tCKAQ0NWoUS+NNrJJZkh1DpTzT6HRlLqpOJxjykFDOrup9L8ooyVfxycgjcFKC7Lpse75asNCmjcRTwFVN4gpJTAN2icoNqg8vquxL6CIe9yf/KTn9hKkk8//bTNpGlSb7VPwaHa46oqquqj67bqFbYN1XVSgaSCK62nluPmc/MGbt62d+nSxRZlGTBgQKvxfH7e11XwpuC3VgXAqEKJqmq0KpTIl0Y7ucxbZgjpz1SzzyFPimXOoipI+v8f9Xcl1L6gapFJXbR0mS9vN0x/qX7/+vjnp1P2Te1SEKIfBVsKSjV+TV0co6pmanyc/laFRt2ndVNGT3O+tW3b1o5fVICj7pFf/vKXbUVM/3xyvT+b6sANJ1D79JpB769775KYW68cVKFEVTVaFUrkS95PLvVlpK4tujrqurWoC04Wu32gfhotUw2EKZY5CztW/IGBTr7jziemz3HvHHf620uvEVQtUqLmk3PrV6xbp/536aWX2oyaAiNlpo488sgW2aKw9XHbR23Saytw69Gjhw26jjnmGHPZZZcVHuu2meYtVeCkQFHP1etqSgFl31TqX10m1RVSAZwyZPpd66Zlap43/c8/n1zXrl1tWxT8uUDoxhtvtOujx6tdCu6UyXTzz1XynqWpCiUBHFpUoQSQPvry0xe9+yLX7/rCynPQiuQ1WqYaCFMscxZ2rPgDP503aRxXnGNKF+Fc4Q2XkVIQ5g0OvcWqXHZJgkr+R82tFvZdoSIjbqJtd7FegU9Ym1ww6CZZ1xg8BWg6h9TvCrYUjM2ePdt84xvfsOvotpnGqimbqPVR108FWXpdraOeo2Ik2p777bef6dmzZyFoUxCmapf6OeSQQ+zYu48++sg+V4GYtr+6nmpuPX03qh2vv/66HT+nCpbex1f6nsURNQTJzXEXBwEcAKScvmi8pZz1e60KtTD/V37kPVON7Kr150yxbLQ7Vly7lNXRY5TxcgGTaOxX0PxtcY5HBUr+4FDBjYK8sG1SShdO7//829m7PspsPfjggzZwEgWQqpTpslVuugAXdCoAUyl/PU/BlzJ3WpYyZ1quKke6XiPqTqn/q8ugW0dl7EQBmpsWR3PD6W+9tp6vKSDUhVKBnzJ9aouCvU6dOrWaA9XNradMooI2rav+1oTfCuZU3VIZTQV1V111lck6AjgAqUYA8b8nF96uNPq9Vt3fGm2KBiBtGuEzsNafM6Vmo/3tUpc+l62SJD+HvYGWgiIFJMpURb3nus9NUO2ydIMHD7aBjpu8Wtkx170yan30WBc0ibJZeq4raCJqk4IxPfall16ymS1NRaAgSxk20fg17bP6nwvIHnroIds18oMPPrBdJkXLddk+BX5f//rXCwVL9LrqXqg2KvBSMKvnKnjT7aefjb/TuLjvfOc7LbaJAj1lvDR+ThlBBXP6240ldHPPZR0BHIBUI4D435MNDdT2joGrVfe3RpuiAUibRvgMrPXnTKnZaH87lHFyE0P7x56VG1RHjadTgCSu8qTGj7kxct7pAPSj8V1uHJ0CIlEb9Rz/5NVh66P/K2hSdk2Bl7oh6jlaJz3fWwnSTbqtoEvBlNrsMm8KmEQBloI7ZeY0tk2/K6BScOgCxvvuu8/069fPtsG/7bwVKhUEankKwhRgKhDbsmWLbadeV9MxaJvofr0faocomFXxE/1f2T23jfR6ebgoQgAHINUIIP73ZEPddOJ21alEIxe+yMOXPBrvMzCL+21aP2f87VIwEzb2rNygWs9x2TOV89e4r2HDhtkuf3rvFHgoc6VbdWtURkyUlfO+popzeLOCrltn0OTVYesjWict87XXXrMBl9ZPFSXVJVFBjwJWZbM0vkzTBOhHwZkCM2Xh1E4Fdgro9D8FUnqM1kf7n9qo+9UOLVcZOWXUFCyqK6Qep8BMXE0GtU3j8FS0RPPEKWhThm+nz6peqmuknnfQQQfZydDVZmX2lJlTRlDz2em1tW1dEKqANQ8XRQjgYDEPHNIqrV/sjaKRC1/k4UsejfcZmIX9NqikvApqhH3OVDMojVq2//PPFfVQWXh1BVTgocAgzoVFPcdls0RjyhSoqT3KMulHv2tsmQIkBVZubjTv2K+w/SPo/2Hro7+VAfzqV79qs4sKvvSaX/jCFwpj8BQcKcOnLJiCLhUU0XQBrgqkd/tde+215pVXXrFZL21L3e8KsigQVCZOQZ7oMdqGCt5cd02dj7r91hVaUSZOr6vHtW3b1gaPCgLVVdSNh1MwqwqY3vGE2oYK/vR4NwXBf/zHfxTap6yc267VwDxwqKq8zwOXxSuS+F+NHECkQSMXviD7iyx+BmZhv/UHmTpJj/qcKRaUVvIdH7XssEIjbuJpV1FRr1luG/wTaSswUZEP0QVznfS77oCua6ICPPe6bj8I2z+C/h9VOMVVdLz44otbBX5uH1KgpbFqCvIUMCkYW7Roke2i6KYH0Gsqs+ctjqIsmZ6jYKV///52XV3SQF0dg7jXdLfaFnq+bjt27GiDOGXgXObPbRO9D+o+qgycgjNl7dS+o48+2i5H+5ra7QLnao+LYx44VFXe54HLwhVJBGvkAAL1lcfsLxez8v8ZWOp+W899odwgs9jjK/mO9y9LgUZQ6XzvY102SWPAFPQEFQgp1gb/RNrKxujWBRXKOClYUhdGZecU6CiQU0CkQMY/pUCc/SZou55//vmt9gtvd09XLEXdEhWkucnGVeZft3fddZd9nraNAieNcxswYEAhCzVnzpzClAEK5LQu2ob6O2h+Nrc/K2unLpzaBocccojtKvqVr3yl0Ca3TVQpU9tL+7fapx5mqjzpXUd1CVU2zzsurlqYBw5Vlfd54LJwRRJAuuQx+6uTG+/kwTqBqsfYSpjUVVis5YXNci+OKEvh3W9Hjx6d2He8vy3KxoSNb3OPVSZKmTcFHu6+ctugAOu6664rvFdu8mlx1SMVKE6ePLnFeyV63XKC7bBgPeh9CAr8goqlqKiWujy6CslLly61XRqV8dLrKTv2+c9/3nZ/nDVrln2sgk+Ne1OXRf1fXTZdAKMsnba9gikFh/psuuKKK+x7r/NTvaYynwrgNnmm1dH74N0mapeSDnotjeVTUKnt6K3Kqce57q/uudXCPHBAg1xJ5wo5kA55zP66uZJEt/qbAC5f4lZYrOWFzaQvjlTyHe9vi388lHe7RLU7Thu875X3u99VmBT9z1W9VHCj7onlbi9vAKggRsGYxoO54Mhb2TKsnf5iKY7LRmqMnoI3BUXKvimj6A203Gv7u4Kq0IkCZv1o+S6Y8nftvO2222y7XcZy//33t9lJ3frbrvFuKpQibvydtyqnhD03Kwjg0BCydCWd7p4AgDxf2Cz34ogCDJdpcXOk+cvpx/2OjxoX5t8uUe2O04ZiF2xdW7zBTZzzAW8Q6rocankuOCrlIo5/zrn99tuvEPyp3QrYVKFSmTP9ru6O6qY4e/Zsm/VSZUg91t8V1N99MWwy8qAgcq+99rJFS/zdEpWtdQVPFMCJgjjvJOyu4ElWEcChIWTpSjrdPQFUizu5cV2h3NxSpaB3QL7U+8JmOfuTN9h0WRyNbfJe5EzqOz7udolznlHsgm1S5wPe7eeO/bBlhr0v/m6UCt7mzZtnAyEZPny4fU9efPFFmwGTlStX2mIip556ql2GspvqZun9/FGAFRQwh1XW9I/F27BhQ4v95swzz7TLVMZNAZ7ap/VR5k6FT3Sf/u/vhps1BHBAymSpuyeyhRNwuJObOCft9A7Il3pf2Cxnf/IGVcraqBucG9P05JNPFh7j/zwr9pnn7g/qoljtz0Z/tUTv1ABJnw94t5+CHjfWLmiZYe+Ly4DpudruDzzwgA3GPve5z9mASOPZTj/9dBukad007kvj2kTj+/R8rat/3FdYwBxWWdM/Fm+Gb79x+7WbB09TDKi7pvYXjYd75513bOZR4+n0eZjV70ACOCBl6n1V1IsT/nzhBByVnLTTOwBJKmd/8u63rluhOzFXdUL3eeb//vSWqw/6zHPVFTVhtsvk6KS/Fp+Neg29loIK/9QATtC8bTfddJPtQirKYCkIccsL+64OG2vnPcdw/4+a406P12PUXj1etH0VlKlQicapucqaqjqpTJyyXm4snMa+ubneRF0qwz6Twiprum6UH34WwD///POFbeVdZxf8utdWO7St3QTlCi6z/B1IAAeLibzTo95XRb044c8XTsBRCXoHIA37kwtqlHlT8KZAw32e+b+zdHKv4EIn8Zo/TH+7SasVYCgYVFZIWSNlphWUKChw46aSEhQ06XdXAMRbBj/qwqnaqy7QrhCRfndtLfW7Ouwcw2077xx3ap+KfagapAsglTXUGDYFZ9q2qvToXYayYpoCQeuhLJkqSmrbR00VEGe/WfzZJOjeAN67Xm5KAT1G77WbJN2dy6pttfgOZCJvVFXeJ/JGPJzw5wsn4MhL7wBkX6VFP1x3PlcK3k047TIzCopUsr5Xr1422FBWSI/VBWsFPiqqod/1eHcirQBDz0v6szHoYqj7PPaXwfcWUfEHJlo/b2VH/a7un1rfqK6YpQia407Bm7axG2+oSbBPOumkwlQKKlqi+xQkuXFn6oKq+evcMtSt0gWhQcFpuT19ogJ4Ly1DbdJYN7VFRU0UbOr/aq/eczdfXDV7FjGRN6oq7xN5Ix5O+POFE3DkpXcAsq/Soh8KLPyl4HW/6+Inyroo+6GMkLhslQIfBTyqpOjGRynIUwCnMWJJfzaWOmF22GMdPc7Nuyb6XWP3FJhocmwFBApMv/Od75TdxqA57tw4sqB2K2uproh6fWXiXOCnHz3fP09e2HseFbAG8S7jTc/5iQKha6+91tx55532/ezfv7/dTpraQNvJzTmn9rqpCILG0CWNibxRVXmfyBvxcMKfL5yAA8jymGlvMKNAw18KXq/puvjpBN1N6KwKiRoz5+g+BXUK2lSEQ0GelqXsUjXaXc7F0KjHqm0a1+cdA6exXCrMofFdovV13UPDCoO4rI+bikFdI7XcBQsW2PL/2jYK0BTs+LOcjv531FFH2QygK8+vIE5/65yy1HMG73uqtuv9K2UfOiVgvKMSEW5snsbiKYA/8MADW4wxVLu9UxFUu2dRtSbybtPs3nE0rCFDhtjbuXPnmkZAYQ7grzgegOxLw3Hsnz8t7pxllb6G9zEKCFyWzhu0uDFwKrzhLQZSrW3mf38UMF155ZWFTJACH/1oXUp5L72PUTdABU0uGFC3QgVQ3mIh/smxXTDrsmbutdzfogyad/t5s5ze90AZP29AVO777n2//O3SepQaVF9xxRV2WeoaqfdWAZ0Cp2HDhtn10O/f+MY3bPCrrKIr0lKN/bQW599k4NBwKMwB/BXHA5B9aTiOazFmupReId7H6CR98uTJJQdm1ZpWIGiycAVv3gqUrohKKYG4vyvp0qVLC1lHlwWLei/cODpXCMQVHHFjB9XtUMGUxo0pGPJOtu1flu5z3S3j9NQJmiJi4cKFdnydsoEKthSM+Scb3+DbVgrKlTlUQPnJJ5/YgFYZOGVZv/SlL9n/q4CKlq/11DZywWEWEcCh4STxJZOGq51AEihUA+T/OK7Fd1YtxkyX0g28kq7iLjB65ZVXbGCjAEbdA6MC4jjbVo/1jmVzhVNKDcT9XUm/9rWv2cDEdY9Ut8onnniiRXbPLVP0fwVpr732mh0/pyBJf7v51RTsKABSwKMgU9vDbQctX0Vg3GNVJKSSiwVBU0So+Iy6cKp6pbaRsmb+AG6Gb1tpGJCK8f3sZz+z66LxjaqAqUIqCjC1PLe9lF1UsJjli5UEcGg4SXzJpOFqZ7WkMThNY5vygkI1yING/4wodhzX4jsrD2OmXWDkMlTuNurCVpxt68ZsuYyXK5xy4403BrbHzx9EqUCHM2/ePDuxtgIX3a8ukMpCurZqmQq6XPdRBWbKXCnT1r17d9ud9MEHHyyUvneVOV17FFQpwFMRGPe8pI5Ptw8pm6jX0XqEWeHbNuoaq/GQykgGdbP1d7/N+ncdARwaThJfMnnOWqQxOE1jm/IiDyddQLU+I7ISGBY7jmvxnZW2Iklx3js9zpXkd+XlNX7KZbCCxNm27v1St0k3Jk5/q9uid3xWsSBDmSYFU+qOqSBQz3Pj6tz0BMo0ufX2vj8KeDTFggsilbHSRNyu/L53snQtz20fFRnRODv9iAsEkzg+3T6k4NYboCqoLPWixSmfzf/mCr0oANV7mbfvOgI4NJwkvmTynLVIY3CaxjblRdpOuoA0fUZk5eJRseO4mt9Z1Q5yg5YvxV6z2HsXtFz9XHjhhYXpBnTrzWBFBX3ebFqxucWCxsSpjaWOz1LwpQBtzpw5tn1r1qyx48X0o4DHO1ecsnX+ipR6fbdPeAuXaHu47eCmCtBYQAVrChT1t5tqQGPLSn2/yj0+VVBG297/nnuX3bVrVxucuUqa3m2loPLVV1+1v7uuk+qCmcZjN654kw8ADU4fFErL68qWbrN+JcfL/8WehuA0jW0CkB7V+ozI6sUjnejqpN1V5lOGp1rfWS5QcpM96+8kBS2/lNcs9t75l6Gsjf6nqo7quqfJqEeNGmUrMBabWFpBlMsWubnFyuHa5sZn6UfBRtjruv3bjRVTN0iNZVP2TQGZgkj3XkvQtnLnMa7giYJHd78LMK+66iqbldNUCxoP+Nxzz9nXcePU9KNsmfYxV74/6L0p9/h0r68ukd7t4F32ihUrbHdRPz1GVR3VTv2oGIrLxuUJGTgghjxnLdLYzSCNbQKQHtX6jMhqbwt/9unxxx+v2ndWtYPcUpYf9L9i753/OTrJV+n9PfbYo1AZUoFUKcGGgrxK5hbzttWV7lfwHZZd9I8VUyVIdfvU3/7Km1pO0Hq785igybr9tG1cwRV121SVx913390WOFGxE2+Gs9ik5cqcKeiLWr8w3mUv/izjqe3ufX0FpG5ydrXNu255QgAHIPXBaRrbBCD/nxFZvXhUy8xhtYPcoKqHboxW1GsWe+/87XZcl0JXvVDPK9ZNNGgblNO11NtWVZBUt0A3li2ohH7YWLHNmzebadOm2UmsL730UpuJK/b+lPv+KYBTl8oePXrYAE6BkncfC1peUKVJPU/rrDF1pc715l32ps/W2XGvry6f2j/UnVRzwknQGLqsYyJv2IkEdYVCHwJxZpIHAADpUYtJtWs1Bu6mm25qFcBpjFSlr+lvt4Ihlc13NL7KjcNSt0p1jVQXx6DtGbQN/BNel/oejB071vzlL38pdI086KCDzP333190HbSNdC6nsWkKstSFUtvKO4F50LYq5f3zvgfK9Cl40+so06VgTplKt35Bk5YrA+z+VoZM6+UKpChrOHz48BbbJ6xN/gnM9w54T/7pn/7JvPjiizYAFnX/vOOOO+pWfEhVOsMqdSqw1PvkJjAvBwEcbACnvs39+vULfczEiRPNpEmTatouANmsiAegvvL0WaGudt5ucDrh1diopEUFdG5ONW9lx2JtiNvuwYMH22BL1FVRz/vOd75T9H089thjbdCn8znXbj1PAU4pwWPUPuMvHiLKFEZNdu6eowybt7KmiqGo26O2qUIQFwB6t08pFyA2hLS3lhcvSjF16lSbFQ0ye/Zsuz0VjJaLLpSwtNNPnz499H6yb0D9ZaUiHoD6ylO381qNQ/RvM+/YMX9lx1LaELfdym4pY6MCHOoCqLFm3oIgYe+rxuE9++yzhedpbJwbx1dKF9qo75c4+5NbnjJ2CtRcWxTwKajSHHIusPNvn6guwP7A7fzzz28ROKat2/O4cePMmDFjIjNwcRDAwdIOpO4CANIrqxXxACCuuCfklWYhvQGYggwVFtG5Uilt0Gsrg7dgwYLCiXqp7f7GN75R6K749ttvmwMPPLCkz3yNefvmN79ps3AKOBUIuqCzlOAx6e8X93y1RV0lXVuUrVMwGPW+RgW/es7iz6ZteOyxx+x4v+uuu67w3qbt4kXUECQVWYmLAA4AMiKrFfEAIK64J+SV9ljwBhj+yo7FgsampiazevVqO4ZOUxK49SiFdw40N86rlM98FSxRNsoVCFGAU860EXG/X8ICZf88c/62RL2vUcGd/rd48eJCVUyNdWvE3igEcACQEWnrGoLK5GmcUtqwbfMlzmTelWaUKuk2qLnH1JVRXSGVZVEZfn8lyVJeN2y9w/Zv73eEm5C71P1ej9VceG7ONPXKCpqQPGqsoDdQDmuLm6Mw6tj0r4ubYNwFho899ljhscrwNWJvFAI4AMiItHUNQXrHNNY6gElbwFTP8aJp2xZ5WN+g99P9HvYe16PHgj+QUHfGSrdj0H6rQMtbmVPdO5WpCxsXVgo9Xpk/zYUnCsqCjhv/e6Fuou453m0Q1vZSj82wx2m7PPXUU4VpFpTha8TeKARwAADkbExjrQOYtBXYqed40bRti2oHYbVY3ziTedejx4ILGtWd0ZWxV5XFUuYhK2c7eifW1u2dd95p51Ir5blRvNtQ3TBVQdK/TxQ7looFU6Uem2GPUxuuu+66ou9tUOVMVfcsNq1CVhDAAQBQB9XMENQ6gElbgZ16jhdN27aodhBWi/UNez+j3uN69FhwQaMybyoi4i2xX4x/u2m+tGJdDR03mXbYsuJsZ40zE38FTP97oeDUjdmLCpRdQKU5z7zTCoQdm5UewzM8+7Sb51hVML2/Z+UCSxACOAAA6qCaGYJaBzBpK7BTz/GiadsW1Q7CarG+Ye9nLd7jcrKTUUFjseV4t6MCMk0J8Nxzz9lughpn5g00FDR5u1D27NmzxWvFfQ+821lFRxR8uoIoTz75pH2Mf2LuUjNYLqDSMrU8BajKGoa9b1HH8IyQSpTuPj1HgaJeS0GidwoI7+9ZucAShAAOAIA6qGaGoNYBTNoK7NRzvGjatkW1g7Ak1rdYcBP2ftbiPU6qi2ix5Xi3o8a0ae40/aiLpAIVZbmCKlVqewUFVXF4t7ObEFuBkNqgrqD6W69TafdMcV1AS2lL0LIW+ypRalzgvHnzCmPjtm/fXph7Tn873t+zcoElCAEcAKChizDkUa0DGArsZHtbFAvCoj4TkljfNI8bTKqLaLHleLejJhFft25dIUBR1sgbaARtc427S/Lz3e0TyrwpeHMTblfaPdMFXi4gjPNe9w6oRKlxgW7ScC1fmTdlEfUzevTowrg37+9ZucAShAAOAJCZkykAySsWhFX7MyHN4waT6iJaznJceX4FOwreFJwlFWhEvZdh0xZUsv5umeoy+e677xYKiVQSEJ4SUIly0aJFNoDT63Xo0MFOg3D66afn9ruLAA5AppAdqv72KPVkivcCyK5yjt9qB1hpGjfo3y7qnnj//fcXnR+tGJfRcoGMCp2om2LQeD69pqibZNKfrVHvZVBwF6eLrHcbusnIlRHr06ePvd/dxn2vuwRUoly1apUN4DT3nubga9OmTWaza6UggAOQKWSHqr89Sj2Zyvp7QQCKtKrX3Gphx6+/wIbGaKmbX1JtS9O4Qf920ZivUuZHK/beuSyngrZddtklck67UsaZlbOPhAVU/s/3oOCukgnNRVkydQfVWDRRZc6+ffuW9F6X03V3yZIl5oMPPjDt27e3WbnDDz8815/nO9W7AQCQl642edke+pLUF6zGDug27As26++FO8nwlskGGmXfLOf49X4muEmjk2ybOxm/7LLL7G3cE2+d8CtAUnCpW/2dxHYpZ1sVe+8qXX6prxP22D322MNWuNTYtrfeequQ7Qu6WBc3E+ptv4Ipb+VHf8XMpNbx4IMPtkHi8OHD7a3+zjMCOACZktQXTF5UY3uUejKV9fci6wEo8qtWc6tF/R32mbD//vsXsjfValutAt+ggC9ou5SzrYq9d5Uu37tcZUNVJfKJJ56wk26HBayuDXq8AreVK1faQE6BuLJ95V68K8bbfo1P0xg+t0wp9T2Ke5GhbwVtzwoCOFg6kNQtIOxHfYqBNGi0D+k0b4+svxdZD0CRX7XYN+Mev2k+bsoNfIMCvqDtoh+NfVuwYIH9UYGRsGCp2PZRxkuZL28GLM57oeW6io6qvKh1CAuGXBv0eLVb3Qzd9ATebZRUJtS7PgrgNF7NLVNFTEp9j+JeZDi7grYnTefPYefWO3bssO9bHG2a9a6joQ0ZMsT2T+7Xr1/oYyZOnGgmTZpU03YBQDUxBg5pEFb5r9i+Wa/9N83HjZu7zFEQETV+S5k37wm0Ag4FAJUsu9j2iVpOuePavv3tb9vCHa4So8bqBbXfLXfatGn2944dO5qdd97ZFvrQuV0txy6X8x6leV8r1dSpU+12DzJ79mxbkVPjEstFAAcbwOkDbObMmaGP0cGuHwAAUL+go9LnZUHcE/dyn1fONiwn2IsStZxy31M9Xlk073QDynaFrXO5j6+GPARl5WbgwnqxjRw50r7/6gZbLqpQwtIOpO4BAAAg/ePd8jyGM06F26jAIOy+cqpf1mI+uHLfU3W/vO222+yYts6dO9txbVHbyq1vNaYnyPNE95WISoC0axc/DCOAAwAAqKKo4CJuYBD2vDxkOCqpyhgU9IXdV04wkdRUB1HLKXdfUAESBQGuOIi64u21116R+0EjBU95RgAHAABQRVHBRdzAIOx5WZ+fMW5Q66oyuu6Byki5oCWJbGVSwU/UcsrdF/Q4jX9TQRLRerttVav9IA8XDLKIAA4AAKCKJ7RRAUTcwCDseXnoWhknqNXj9BwXzLiqjNpGSXV/LFWc7pxx9gU9XxUxFy5caN544w3bVc9VyKzVfpCHCwZZRAAHADnFlVGg9rwntMoGXXjhhfZYVEChSoGaQ62aAUStg5VqiBPU6vNNc6GpsqKryuiCFhcQLlmyxLz77rtm27ZttqBHqRU/y/2sVRClMvGi9+K3v/2tHXem+9XNUfOvaT+oNODxrtchhxxi11mvq//Vaj+IGyjGDXLxv5gHDkBZk5wiO8qd1BZA5bwnsArgli1bZg4++GD7t060qz1nYtScYrX4TK/X94ZO8E866SQzfPhwM2jQoBaBsgsI9T706dPH7LLLLoXPxCQ+J/3LeOSRR1rcr7/d/doftF8kkRlz66X1devsllmreTrjzgsYtd357iqOAA5AKD5Esy0PXanwV1xQyQbvCazGJCkbpBNrd5Jd7UmGoyY0rsVnej2/N4oFr8rQPfHEE7Zsu8bL6TMxic/Jcp6j/UH7RZKZMf8yNLdYrTJYcQPFqO3Od1dxdKEEEIoP0WzLQ1cq/BVjTbI3fkvzbKm7XFqOwVp8pvuXqayjLjgk2UUxbBlRXS/1fAWVmv5Y4+SUBXPBRqWfk/7PWs3vJS4Tp+6TChgVyKub43vvvWcDnkqqWUaNGVQXzmp9VgS9F3GWHfX9xHdXcWTgACTeNQLpUKsuNKgNLqhkgzcDpkmSdcKelmOwnM/0uBlf/zI15izpLopxlqHjRe9Fp06d7Dg5vSd6P5L4nPQv48wzz7RB28CBA+1Pz549C0Gb2qD9IihDmlTWdf369a3WPan3OakMa9R257urODJwAKo+7w3qgzl/kpOGQfV5uyqdhm3aaMdgOZ/p5WR8ve+luu9169bNBhF6DRUMSbqLYpxluONH3VhFgYHb3yp9j4LeZ28blXnbf//9bYBVi2PEraubVkGBkCva4n++CqzMnDmz0N1X2bvzzjuv6heSoo6NtB03aUQAByAUH6JAerov5u2CShq2aaMp5zNdXR9feeWVwom9PxALey9VBVHBkQsCFDgk3UUxzjKSPn7Cgiv3f42zq7TyaNxjxK2rxvyJircEPV9tve2222yw3aFDB/Pxxx/bbp9RAVy570UjXKipBwI4AAAy0H0xbxdU0rBNa63YSX+aTnLV9dHNqaZb/R3nvUwicEpiGUkfP2HBlfu/giZlvxQIqzpmnDa77egyaU8++aT9u9j+4dZVz1cQ6V+edx22bt1qf//oo4/MunXrTPfu3RN9L7hQUx0EcAAANFj3xTRoxG1a7KTf//962meffewJvcvA6e8472USgVMaL16EBa3u1lUeVffFuG1321XBm4Jojd8rZ/8odoyprQceeKB5/fXXCxlWV4AlqfciyQs1abzQUS8EcECF+EAB8i9v3RfTII/btNj3QbGT/rDHVas9UZRB0nxpjrpFhr2Gxk0tWLCgEABk6b2sdJyZ9++o/1dyjCjztvPOO9s55F566SXbvhEjRtgqp/71cJOVK+BW8RTveET/++IqVmrZCtS1PBVgSXKbeicud68ZVxovdNQLARxQIT5QGhsBfGNIYwYg6/K4TYt9H9TipL+c9iQRYOsxGvemaoui6otp/gz0f2YrgFH7444z82+fJC9MeI+RG264oTB/nLo9XnnllXb8mv+91rhFZeuUPVUA7h2PGLYOyhK6oE9/V/o95t3vFLypAqeKuFS6PRqx23UYAricUHr997//ve3moH7MuupyxhlnmPbt29e7abnHB0pjI4AHWmrkixrFvg9qcdKf1JxspQbYWfsO9H9mK3Pogs9y2h+2fapxYULv07XXXmt/V7GRvfbay2a2vFy7XZDnbqPWx7VV+4TLtlbyPeaO/WnTppk99tijUMAlTgXOIBRQ+SsCuByYN2+eueqqq8zNN99cGHx6yy23mB//+Md2Xo+2bdvWu4m51ojjOJDdkxcgCVEnRo18UaPY90EtT/qD2qMsSxIn6lGvUWxuuUpOqJM4IS/2GV2r7/By1kX/P/LII80bb7xhtm/fbjOGynR6pwZw74Mu5CsDp9tS1yep7zF37Ct4c5Olaxxg3G3q30bqNvr4449TQIWJvLNPlYmuvvpqe2B7Kwd985vfNHPmzIk9wSJKx4STyYs7gWw9MNk5GlHUZL6NfFEjbd8H/vb4C5Ek8d6Us86VTgKdxCTS/s9ojdmrx3tW7rpceumlplevXjaIU3ZNWTgFSO557n04/PDDbWES3Za6Pkl9j7n9yU2Y/sEHH1S0Tf3bSMGbd8LyLkWC9zx/FpGBq6Ibb7zR7sCTJk0q+Tkvv/yyufPOO+0Jq8ucHXPMMebEE08MzKQ9+OCD9rE6WLwUzGkyTS3r1FNPtR9MqI48juOot1peNav0im4eCzEAxUSdGGWxV0JSXa3S9n3gb08Sc7IVe41qnlCX8/yw99T7ma3zJLecWnexK3dbqMCIgs1Vq1aZ5uZmm4lTAKdMnHjXzc07V+r6JPU95o59V4FTwVslx0Ol+0vvDH4WlYqz+ipZs2aNefTRR8t6jg6eyZMn2/lCbrrpJvOLX/zCXmXQcvR/TbDoN2vWLHvbp0+fVvfpfxs3brQDWpHdDE8jquVVs0qv6LqTl1KvCKI0HKPpFnXFPm1ZqFpldrKwr9f7vak001PO88PeU+9ntoIfdUesx/seZ1vou9B1jRRl4rzPi7sfJ/U9ltT+5Y4JTYauH/U2i7O/nJLBz6JSkYFLmK6KKIt2/fXXm82bN5f8vKefftqOW9OBM3To0ML/lSK/5JJLbBZPXSWVQm/Tpk3hfjeI1V2B8XL/02P69+9f4ZrlS577RedBLa+a5bmLRZZxjKZb1BX7WmShki5OUM/PgVru6/XOEJaS6Yl6b8vJFJXyntbzfY+T9XJVM5V5c2X/vc+r9/dZUvtXUpOhd0lZRjxJBHAJaWpqstkwDdr0lqQthR6vgE/B3/HHH9/qfvVZVyr6qaeeMo888ogZNWpU4T4XJHrnanFUqUjUBxkt1ftDDunplpjnLhZZxjGabvU+MUo66Knn50Aj7eul7DdR7205+10p72k93/c4x5D7btQF+qALF3n5PktyMvS8IoBLyLhx4+yPvPPOO2bixIklP/fuu+8uDPTs3Llz4GMGDBhg5s6da+644w47Js6Nh9t9991tN0lNHeDn/qdqQGgpLx9yeVXLk0PGsKUTx2jjKiW7lnTQU43PgVKzhEnt61ksmR7U5qTe21Le06x9/hf7bsza+oSp5JjYkMHjIA4CuBS4//777a0CuDBKk4sGr7744otmyJAhhf+rf7CyeH7uf+65yN+HHLKfSUAwjtHGVUp2LekAvxqfA6VmCZPa17PY7TiozUm9t6W8p2n+/I8TiKR5fcpRyTExI4PHQRwEcHW2bNkys27dOvv7vvvuG/q4Hj16FH5/4YUXCgGcpg9QAKed/Atf+EKL57z11lu2wtJhhx1WtfZnVV4+5IC84hhtXKVkYLIQ4JeaSUpqX89iV8ygNp9//vmpf29rEaA1SiCS9DGxIoPHQRwEcHW2cOHCwu9h3Sf99y1atKjwu8bMaaqAV1991YwZM6ZFFcz169eb733ve5mYQqBRUt5IDvtM7bHNUYv9pZQMTBYC/KD1iNomlR5ftZxQOylBbc7Ce5uEYgFapYFIJe9xWvaPOHo3SPd7Arg6c9k3N1gzTMeOHUOfoykGpkyZYsfeuSyeiqp8+ctfNieffHJJ7Vi6dKkZPXp0zLUwZubMmaYSjXylCfGwz9Qe2xy12F+ykF0rRdB6RG2TSo+vcrZbWo7luO91vQKMJF+3WIBWaSBSyXtczf2j2u/dKTX+/Kjk3Fnn3v369attAKfS9goaUBlVrXTatQt/O7z3qWiJlyr0/OQnPzG//vWvbVETFS/RTqv5PIIm/06jRkl5IznsM7XHNkct9pe8ZGCC1iNqm1R6fNVyQu2klPteu5P/e++91851duihh9Y0AE0ysPEGaJrn7L333rNzAbqgo9JApJL3uJr7R7UvHnTx7VNuTrksZhOrEsBpTjJ1z+vVq1eyLWow3vL/+jAK473PTQ/gpQ+xyy+/PHY7dAWg0ixaJRol5Y3ksM/UHtsc5WB/KW+b1HJ7ZfW9cSf/upCtqZc0R5guYtcqAE0ysPEGaAre9t577xYTcCsIqdfUGNXcP2p58WDDhg3mwgsvtPUmNAG6ivslHTBWcu7s6lnEEXtwlDtwUBnvVYDt27eHPu7jjz8uaaxcVumDTFU4NV5Pt1ntMpMX7oqVrgbqVn+nDftM7bHNkdX9JS2faVHbpBbby20HTYysQmfbtm2r+3tTDneyr5Nx0WTWtQxA/a9Tyeu6TJF6S+2///4thtEkEdRUsj9Vc19MchsWo2BNwZtiFvV4U9ySl54jFY2B0/xl3kmlUT5vBtN9EAXx3leNnV1XfaImH9cYPO84vKTlpctMXqRlfEQU9pnaY5sjq/tLWj7TorZJLbaX2w7q/dOnTx97cp6W96iczJB6HelkvNYXB6o1vqoaGa+43VOr3dWwlmPUVqxYYYN9N1xJ59K1zjZv3brV/gTZsWNH7EKDFQVwmlT6X//1X0337t0rWUxDU+pf49Q++eSTFsVJ/Lz3VZJyjTpwx48fH3q/JiafNGlS4q+LdErL+AgASAKfafnYDt6TfzdOrJbjmaoVZKehcE+tLnLU8sJO7969bbdJBfsK3jQvcq23rYoKTps2LfA+ZQc13VfNAzhtlEsvvdQWz0A8u+22mxk4cKCZN2+enaQ7jPuQVaQ+bNiwqhxQ06dPD72/mtk31F6xK21ZHR+RZVku29yI64Vs4TMtH9shTVndvK1X1oP7qMBY59r1+v4ZN25ci2m+vEaOHBk7A1fRBGHqU/qb3/zG/OlPfzLVovTiD37wA5NnEyZMKJQTDStk4sYbqlxpNTKe2oG6desW+kMAly/uSpt3wHRax640imLvSVbldb2QLXym5Xc7VDq+MS3jI+utlmPTaqWLZ5yhbutx8VDnz2Hn1qowX5culGeddZadLFobRZNL77PPPiZpqjR08803mxtuuMHklTJwyqo999xzdmJv/e2l7pVz5syxVxC0zYFqX2lLw9XARpPHq595Xq9aaZQMZrXXk8+05LZDpe9V0u91pV3/qtV1MGvHbrndOLO2fnkTOwOnUvYKqpR9O/PMM+2E0d5KiUlZsGCByRpVdXJK3SaajFtViNRPVgGb1wMPPGAPFF1BcBN1V6uISdhP2ABMZFMer7RlXV7fk7yuV600SgazUdYzDyp9r5J+r91FIc2lNn/+fHseVU4mrVoXmbK2T5ebrcra+tWLzp/Dzq3VyzBqCrGqZOD0Zrly9j/72c/M2LFjzd///d+b//7v/zZJmTt3rjnnnHNMFmiQ4qxZs+zv3mqO99xzj3nllVdMmzZtbKr0mmuuCXy+smu679prrzUXXXSR3Z5Ku+rDSMtVan/w4MFVaz9FTBpLGgZMozHek7yuV600SgazGutJhiCd71XS77Ub16ehJqo22KlTp7IyadUaF1jKeia9j9Zyn2+Uz6a0FjFp06yBbAnQQaNugH/3d39nJ82rJHt111132cIo6lKo5in48Wel8kzzs6hLqrJ3KvN7xBFHmPbt21ft9VTVUlcAoiYjrPY0AgCA1pRJ8J5cZq3seznrqRNwb7W46667rqKTz0bZdrXknxhZ5fz1U852Tfp9cUGLTpL32GMP2x7NqaaxRcom1SvoKWU9k94Wtdzn/a+lMV1KRnDBpPRpBFwREyVrajoGzktXPJRt+spXvmL69+9vjjvuuLKe/+qrr9qgTQegxr1JQrFl5hx88MH2p5ZcERMAQHo0SgZT66XAQMGbAoO999674rFIZAiSp/dE742mNtJ79d5779khIPXcp73j+uJk0uo5NUDS+2gt93n/+qkyfRrmWkybqASIeubFlVgAJwo6NDecxsSp21+/fv0iH68Mk1KLv/rVr8wzzzzTImhT1g1AetE9Cai+WhXfqPfxrNfSOHD1Okm6a533b1RG74myW5rD1l38LXc/yfNcauWuZ9L7aC33ef/6aaiPFxdMqquiaQSCHH/88eZHP/qROfHEE+3VmbAugj/84Q9Nz5497dgqBW/+bJuuwGnuBD0GQPowgBnIjzQcz0kXvMljyfxGLUpUSqn/NJSMr/c+Ws99noJVtZXYGDg/lbtXil0VFJVNU6WVP/zhD7ab5BNPPNEq2+Z+13gsFS5RQQ2lHPXYY489tqHGwNUaY+AQh75IvdWTSh1vACB90nA81zsLiPS+R4xnzO++kffjfmuVxsBVLYBTMZKjjjrKVk7U+Lj/+q//KlRndIVJ3O8KDL71rW/ZwE0FO7w0z5zGZsUts4nSAjgF21FdXqlCCT++UNMr71+ISB7HM9J+gUFjrFyRG1VBv+2220I/19L8GZjmtuXps2dDSrbz1KlTQ6tQzp4921ahXL58ee0CuFGjRpk///nPpm3btqGPWbVqlRk6dKh59913A7Nthx9+uPnud79rAwMFeWGmTJliLrnkkjjNRAnIwCHLH45ojZNxlIvjGWn/TNP+qYrnonNG7aNhn2tp/gxMc9vylP2/PSXbOXVVKB977DF7NURj1cJo/Nrdd99tjjnmGFuwRIHbzjvvbE477TTzve99z2boSkHwVn1UoURaiysgX9X3CBTSuZ04npFm2v/vvfdemwRw0xdEfa6l+TOw3m1L22dwqYVXym33ipTsA9WqQhm7iImCsc2bNxd9nOaGu+WWW+zjL7jgArsBf/e735UcvAEA8jOYPA3FMrKA7QT8lU7UTzrpJDN8+HBbAVOVMKM+19L8GVirtoUVfknbZ0uphVfKbXfvFO8DSahoGgEFY/vtt1/Rx/2f//N/zIIFC8yzzz4b2VUSAJC/ctppvCqadrXcTmm7It9I2PbV+VxL82dgrdrmAh7/vGxp+wwuNftfbrtPSfE+UPcA7vLLLzfXXnutGTBgQNHHXnfddXaKAY1505UApIuuaLgiM0EYAwdkR5q7wzE3V/q2U9iJHqovTds+7cFkOZ9raf4MrFXbwgKerH4Gl9vuLinZB6LGwKlCvzKPNQ/gVPRCqWxVUPna175mf9Q1UkUx/MVN1MA777zTfPnLXzY/+9nPzD/+4z9W8tKowge3pm4IQxVKAEnI+1XRLG6ntF2RbyRp2vZpCiZrIe0Ba7UCnqx+Bme13U1NTaFVKJctW2ZjqDhiV6FUQOatKOmmBRD1TdbYNwVz+jnyyCPt/+Tll1+2//v9739vK1mW4vHHHzcjRoyI00yUgCqUANC4ilVrq9WJbt5PqNNcKS8tcwE26ravhkY8ntIodfPABaX8wgI6VVn54he/WAjo1FVv8uTJdv6DAw88MPJ1VOlS4+aUZkT1AjiZO3duvZsCAEjZiV6tTnRr9TppOrFNU1vyHtA0esCadmk6FrJw/l1RF8qrr77aThXw9NNPm6eeesq8+uqrhfvcZN263b59u3nhhRfsj7pPOmPHjrWB3Fe/+lVz0EEHhU4IziTeAABUR7GxIrXq5ler10lTV8FyxulU+wQ3q13U4srqWLAsKmXfTdNxmQUVBXATJkwwvXr1MmeeeWbhDXrmmWdsMKegThGl5n/zB3Si31977TXz7W9/2/7dvXt385WvfKUwlm7w4MF2HN26detaZPMAAED+TnSTfJ2oE8Y0jTsrR7VPcNNS9KFWshawZjlDVcq+m9XjMpMBnH8Sb+1IJ5xwgv2Rjz76yMyZM6cQ0M2aNcts2rSpRUDnfn/33XftG+rmddCYuS996Utk3wAAaIAT3SRfJ+qEMauZF05wk5W1gDXLGapS9t2sHpeZC+D69etnOnToEPkY3X/00UfbH1Ew9tJLLxUCOv28/fbbhcd7s3Qa8PfEE0+0CPQAAEA+T3STfJ2oE8asZV4cTnAbW5YD+FL23awel/USu4hJUt54441CQKfbJUuWtHqMArhPPvmkLu1rBFShBADkSR4LcmS5Cx0ae59u5H13a9qqUFaLKlS6YE7TB2ilCOCqH8BprKGyqmGYBw4AkBWNfMKIfGKfzqapU6eGzgOnavyaB2758uXZD+D8rr/+enPRRRcRwFURGTgAALKh3ifySb1+vdcDqIWGycD5LV682Hz+858ngKsi5oEDACAb6t2VLqnXr/d6JIEgFPU6/249G3fK7L333vVuAgAAQCrUu5hFUq9f7/VIsjKkejG5ypBA6qcRqIVddtmlMHccAACoHjIK6VfvapRJvX691yMJ9QpCOU6R+gzc7rvvbn7729/WuxkAAOQeGYX008m6uhtq7Ixuq11uXcGCujteccUV9nbEiBGJvH456+Fvg/5OA3/QWasgtJGP07TuC7WW+gycjB8/vt5NAJATXLkE8t2tLe9qPQG1fwJpVQhP4vXLWQ/Xhg8//ND+fu+995qTTjqp7p/f9Zq7rJGP0yxPaN5wARwAJIUPfyDf3dqQv2DBvaYK273//vt2eqk0fH7XOph2Gvk4TcP+mAap70IJAEniwx9IT/c8pF+9ugkGveamTZvs7Z577tnq87uRutbV+zit57ZOw/6YBmTgYKkftSZRD8M8cMiLRr5yCaQ1o4D0qlc3waA2dO7c2Z6vHHrooa0+v2vRuyItXfDrfZzWsydLGvbHpOaB27Fjhw3C40j9PHCoPs1DsW7dOtOvX7/Qx0ycONFMmjSppu0CqiEtX8AAgOQ+v5UNUnDn6MT4sssuS/T18zB3XRJqsa3zYurUqWbatGmB982ePdt07drVLF++vOzlkoGDpQ/A6dOnh95P9g15Ue8rlwAQV6NfgIr6/K5F7wq64P8verKUbty4cWbMmDGB940cOTJ2Bo4xcLC0A3Xr1i30hwAOAID6auTy8WkYF8b4q3SMwcsSnT+HnVu3a9cudgBHBg4AACADyADVt3dF1sZfVQs9WeqPAA4AAI9G76aG9KLrWn0RuCAt6EIJAIAH3dSQVnRdAyBk4AAA8KCbGtKKDBAAIYADAMCDbmoA/OhajTShCyUAAB50UwPgR9dqpAkZOAAAPOimlj5kP1BvdK1Ovw0N9DlBBg4AAKQa2Q/UG3PApd+MBvqcIAMHSzv72rVrIyciZDJvAEA9kP1AvTEHXPqtSOHnxNatW+1PkB07djCRNypPO48fPz70/okTJ5pJkybVtE0AAAiFZVBvdK1Ov94p/Jxoamoy06ZNC7xv2bJlpmvXrrGW26a5ubm5wrYh44YMGWIzcDNnzgx9DBk4AEC9NNLYFgD5+ZzYGpGBGzlypM3AzZ8/v+zlEsDBBnAyd+7cejcFAAAAyL0hFZx/U8QEAAAAADKCAA4AAAAAMoIiJgAANPg4DQBAdhDAAQBQh7mKxM1VRHU7ANXCRaP8IYADcooPbCCd0jhXEYD84qJR/jAGDsj5B7amiHAf2ADqzz83URrmKtIFn9tvv91cccUV9lZ/A8gHLhrlDwEckFN8YAPppGx437597fw/utXf9cYFHyC/FzzSeNEIlaELJZBT+oB2XSbc3wDqT12Z09Z9iQs+QH67JOoikX9IBbKNAA7IKT6wAZSKCz5Afi94pPGiESpDAAfkFB/YAErFBR+gdFzwQL0RwAEAUolKqrXDBR+gdFzwQL0RwAF1xkkqkI9xJgAaAxc80mVDA55HUYUSqDOqvwH5GGcCAKi9GQ14HkUGDpZ2+rVr14be37FjR/uD5HGSCgRjnAkAIMvnUVu3brU/QXbs2GGnk4mDAA6F9PP48eND7584caKZNGlSTdvUKDhJBYIxzgQAkOXzqKamJjNt2rTA+5YtW2a6du0aa7ltmpubmytsGzJuyJAhNgM3c+bM0MeQgaueRuy7DQAAkPfzqK0RGbiRI0faDNz8+fPLXi4BHGwAJ3Pnzq13UwAAAIDcG1LB+TddKAEgxdJ8ZREAovD5BVQHVSgBIMUasboWgHzg8wuoDgI4AEixNFfXAoAofH4B1UEABwAp5q+mlabqWgAQhc8voDoI4AAgxTRmpG/fvrZSlW4ppQ8gK/j8AqqDKpSgCiUAAACQkfNvMnAAAAAAkBEEcAAAAACQEQRwAAAAAJARTOQNAACAxDGRN1AdZOAAAACQOCbyBqqDAA4AAACJYyJvoDoI4AAAAJA4JvIGqoMxcAAAACmV5XFkaqu/7QAqRwAHAACQ8nFk4saRnX322SYLFGhmpa1ItyxfyKgGulDm1PPPP2+mTp1a72YAANDiJOz22283V1xxhb3V34jGODKAgjh+BHA5tHr1avPLX/6y3s0AAKCuJ2F5CBgZRwZwIcOPAC5H7rrrLvPP//zP5kc/+pFZuXJlvZsDAEBdT8LycNVeXcX69u1rdtppJ3vLODI0Ii5ktMQYuBA33nij6dSpk5k0aVLJz3n55ZfNnXfeaa/wtW3b1v7vmGOOMSeeeGLh72o644wz7M+CBQvMD3/4w6q/HgAA5dBJlxvP5f6upjxctWccGUBBHD8CuABr1qwxjz76qDn11FNLfo52qltvvdVcdtllZujQofZ/69atMz/+8Y/NM888Y6ZMmWLat29fxVYDAJButT4Jq3XACKA6uJDREl0oPZqbm82iRYvM5MmTzebNm0t+3tNPP21uueUWc+aZZxaCN9lrr73MJZdcYl566SVz9dVX2+UDANDoJ2G62KnbaleRo/shgDwiA2eMaWpqMrNmzTLvv/++2bJli1m7dm3Jz9Xjr7/+ehucHX/88a3u32effcygQYPMU089ZR555BEzatSohFsPAACCcNUeQB4RwBljxo0bZ3/knXfeMRMnTiz5uXfffbf54IMP7JW9zp07Bz5mwIABZu7cueaOO+6wY+I0Hm779u3mvPPOq6jdJ5xwgh1fB+QJc70AAACEI4Cr0P33329vFcCFOeCAA+ztqlWrzIsvvmiGDBlidt55Z/PrX/+6Zu0EsiLLk9YCAABUG2PgKrBs2TJbqET23Xff0Mf16NGj8PsLL7xQk7YBWZWHqnEAAADVQgauAgsXLiz8HtZ90n+fiqSk0dKlS83o0aNjP3/mzJmJtgeNi6pxAAAgC0ZXcO6sc+9+/frFei4ZuAq47JvsuuuuoY/r2LFj4HMAtEbVOAAAgHBk4CqgqpVOu3bhm9J738aNG6vWnvvuu88sWbKk0OVMVS/fe+8906FDB3PuuedGPldXAMiiIQ2oGgcAALJgZgXnzqqJERcBXAV22WWXwu+ffvpp6OO89ymYqhYqUgIAAAD5RgBXAW9pc00LEObjjz8uaaxcPSnIjJr/Tt1AvV1BAQBIE6YgAZA2W7dutT9BduzYYYeLxEEAV4FevXoVft+0aVPo47z3pbUgg774xo8fH3q/5sabNGlSTdsEAECpmIIEQNo0NTWZadOmhVaz79q1a6zlEsBVYNCgQXZS7k8++SSyOIn3vkr6u1aTrlJOnz499H6ybwCQP3nKWjEFCYC0GTdunBkzZkzgfSNHjiQDVw+77babGThwoJk3b56dpDuM+xLRmzRs2DCTRmpbt27d6t0MAEAN5SlrxRQkANImaghSVAHEYphGoEITJkwozOUQVshk8eLFhbkiunfvXtP2AQDQCFkrpiAB0CjIwFVIGThl1Z577jk7sbf+9lL3yjlz5ths3VlnnWXSiiImANB48pS1YgoSAGlDEZMa2bZtW2D1yCiTJ082F1xwgR2k2L9/fzsuznnggQfsGIMpU6aYfffd16QVRUwAoPEoS+UfAwcASHcRkzbNzc3NpsFp486aNcv+rizU6tWrC3O2aYLrNm3a2H6q11xzTegy1q9fb6699lob9I0dO9Zmq+bPn2+X+4Mf/MAMHTrUpJUKqygDFzUZIRk4AAAAIJkMnCtionihXARwCVuyZIlZtGiRDeT69OljjjjiCNO+fXuTZq4y5ty5c+vdFAAAACD3hlRw/k0XyoQdfPDB9gcAAAAAkkYVSgAAAADICDJwsKhCCQAAACSHKpSoKqpQAgAAAMmhCiWqhiqUAAAAQDaqUJKBg6UdqFu3bqbRso7++Y80ESwAAABQqagEiKYoi4siJmhYCt7efPNNm33Urf4GAAAA0owADg1LmbeovwEAAIC0IYBDw1K3yai/AQAAgLRhDBwadhoBjXnzj4EDAADIM2oAZH8aAapQwlahXLdunenXr1/oY5hGAAAAIPtuv/12O/bf6du3rzn77LPr2qa8mjp1aug0ArNnz7bTCCxfvrzs5ZKBg6UrL9OnTw+9P2/ZNwAAgEZEDYDaGTdunBkzZkzkNAJxEMChYacRAAAAaDTqNunNwFEDoHqYRgAAAABARTTmTd0mdfFet9QAyB4ycAAAAEADDZthzFu2kYEDAAAAgIwgA4eGnUYAAAAAyNo0AgRwKMwJMn78+ND7mUYAAAAAKF1TU1PoNALLli2z0wjEwTxwsPPAKQM3c+bM0MeQgQMAAACSycC5aQTmz59vykUGDhbTCAAAAADJYRoBAAAAAGhwBHAAAAAAkBEEcAAAAACQEQRwAAAAAJARBHAAAAAAkBEEcAAAAACQEUwjAEvzwK1duzb0fuaBAwAAAJKZB27Hjh12Gq84COBgbdiwwYwfPz70/okTJ5pJkybVtE0AAABAVjU1NZlp06YF3rds2TLTtWvXWMtt09zc3Fxh25BxQ4YMsRm4mTNnhj6GDBwAAACQTAZu5MiRNgM3f/58Uy4ycLC0A3Xr1q3ezQAAACnqnTNjxgyzYsUK07t3b3PKKaeYLl261LtZQGZEJUDatYsfhlHEBAAAAK0oeHvzzTdtLx3d6m8A9UcABwAAgFaUeYv6G0B90IUSAAAArajbpDJv3r8BoXttfZGBAwAAQCs6Ke/bt68dJ69b/Q0I3WvriwwcAAAAWlFG5eyzz653M5DCDBrda+uLDBwAAADQQCrNoPm709K9trYI4AAAAIAGUmkGje619UUXSgAAAKCBVFqghu619UUAB0sp9LVr18aaiBAAAADZoYyZfwwckrd161b7E2THjh02gxlHm+bm5uYK24aMGzJkiFm3bp3p169f6GMmTpxoJk2aVNN2AQAAAFk1depUM23atMD7Zs+ebbp27WqWL19e9nIJ4GADOGXgZs6cGfoYMnAAAABAMhm4kSNH2gzc/PnzTbnoQglLO1C3bt3q3QwAAAAgFzpGJEDatYsfhlGFEgAAAAAyggAOAAAAADKCAA4AAAAAMoIADgAAAAAyggAOAAAAADKCAA4AAAAAMoIADgAAAAAyggAOAAAAADKCAA4AAAAAMoIADgAAAAAyggAOAAAAADKCAA4AAAAAMoIADgAAAAAyol29G4B0+PTTT83atWtD7+/YsaP9AQAAAFDc1q1b7U+QHTt2mJ12ipdLI4CDtWHDBjN+/PjQ+ydOnGgmTZpU0zYBAAAAWdXU1GSmTZsWeN+yZctM165dYy23TXNzc3OFbUPGDRkyxGbgZs6cGfoYMnAAAABAMhm4kSNH2gzc/PnzTbnIwMHSDtStW7d6NwMAAADIhY4RCZB27eKHYRQxAQAAAICMIIADAAAAgIwggAMAAACAjCCAAwAAAICMIIADAAAAgIwggAMAAACAjCCAAwAAAICMIIADAAAAgIwggAMAAACAjCCAAwAAAICMIIADAAAAgIwggAMAAACAjCCAAwAAAICMIIADAAAAgIwggAMAAACAjGhX7wYgGe+//76ZPn26vV29erXZuHGjGTZsmJkwYYLp1KlTvZsHAAAAIAEEcDmwbds289Of/tScc845pmfPnvZ/ixYtMhdffLF59tlnzY033mg6d+5c72YCAAAAqBBdKHPgqaeeMkOGDCkEb9K/f39z3HHHmXfeecfcfffddW0fAAAAgGSQgQugjJW6HU6aNKnk57z88svmzjvvNBs2bDBt27a1/zvmmGPMiSeeWPi7WtavX28efvhhs8cee5gRI0YUXu+AAw6wt6+//npVXx8AAABAbZCB81mzZo159NFHy3rOjBkzzOTJk81JJ51kbrrpJvOLX/zCXHbZZXY5+v/HH39sqmnw4MHmzTffNFdffbV5+umnC//funWrve3Ro0dVXx8AAABAbRDAfaa5udmOG1PAtXnz5pKfp4DplltuMWeeeaYZOnRo4f977bWXueSSS8xLL71kAystv1oOOugg8/Of/9yce+65tiulM2vWLLPnnnuaM844o2qvDQAAAKB2Gr4LZVNTkw10VL1xy5YtZu3atSU/V4+//vrrbXB2/PHHt7p/n332MYMGDbJj1B555BEzatQoUy2HHXaY/XH++Mc/2nVRcZO99967aq8LAAAAoHYaPoAbN26c/REV/Jg4cWLJz1VxkA8++MD07ds3tMrjgAEDzNy5c80dd9xhx8RpfNr27dvNeeedV1G7TzjhBDu+LojG4mlMnjKDysABAAAAyIeGD+Aqcf/999tbBXBhXCGRVatWmRdffNF2cdx5553Nr3/968Tbo8Dw5ptvNn369DFXXnll4ssHAAAAUF+MgYtp2bJlZt26dfb3fffdN/Rx3gIiL7zwQtXao0IpN9xwgxk7dqw57bTTWtw3c+bMqr0uAAAAgNohAxfTwoULC79HTZLtvU9FUqrlN7/5jRk+fLjp169fi/+vXLnSdg0tZunSpWb06NGxX58gEQAAAI1kdAXnzjr39p+3l4oALiaXfZNdd9019HEdO3YMfE6SPvzwQ/OnP/3JvP3227Ygi/PJJ5+YefPmVbV4CgAAAIDaIYCLSVUrnXbtwjej976NGzdWpS1vvfWWHf/23HPPxV6GrgCQRQMAAABKU8m5s3fqr3IRwMW0yy67FH7/9NNPQx/nva9Dhw5VacshhxxC8AUAAAA0AIqYxNSlS5fC78p+RRUXKWWsHAAAAAAUQwYupl69ehV+37RpU+jjvPf17t3bpJUyhVGTmGssn3c8HwAAAIBwW7dutT9BduzYYXbaKV4ujQAupkGDBtlJuVUoJKo4ife+Svq6VtuGDRvM+PHjQ+/XBOeTJk2qaZsAAACArGpqajLTpk0LnZKsa9eusZZLABfTbrvtZgYOHGirPGqS7jArVqywt4qwhw0bZtLcJXT69Omh95N9AwAAAEo3btw4M2bMmMD7Ro4cSQauHiZMmGADOM3joC6IQW/C4sWLC/NEdO/e3aSV2t6tW7d6NwMAAADIhY4RQ5CiqtgXQxGTCigDp6zali1bWkzs7ah75Zw5c2y27qyzzqpLGwEAAADkBxk4j23btgVWj4wyefJkc8EFF9j+rf3797fj4pwHHnjAji2bMmWK2XfffU2aUcQEAAAASH8RkzbNzc3NpsEHF86aNcv+rgBm9erVhTnbNLl1mzZtbIrzmmuuCV3G+vXrzbXXXmuDvrFjx9pAZ/78+Xa5P/jBD8zQoUNNmqm4ioqtaH3DUMQEAAAAKN3UqVNDi5jMnj3bFjFZvny5KVfDB3BJWrJkiVm0aJEN5Pr06WOOOOII0759e5N2CuCUgYuaDJwMHAAAAJBMBs4VMVHSp1x0oUzQwQcfbH+yiCImAAAAQHIoYgIAAAAADY4MHCyKmAAAAADpL2JCAAdL1TLHjx8fej9FTAAAAIDyiiWGFTFZtmyZLWISB0VMQBETAAAAIGEUMUFVUcQEAAAASA5FTAAAAACgwRHAAQAAAEBGEMABAAAAQEYwBg4W0wgAAAAAyWEaAVQV0wgAAAAAyWEaAVQN0wgAAAAAyWIaAVQV0wgAAAAAyWEaAQAAAABocARwAAAAAJARBHAAAAAAkBGMgYPFNAIAAABAcphGAFXFNAIAAABAcphGAFXDNAIAAABAsphGAFXFNAIAAABAcphGAAAAAAAaHAEcAAAAAGQEARwAAAAAZAQBHAAAAABkBAEcAAAAAGQEARwAAAAAZATTCMDSPHBr164NvZ954AAAAIBk5oHbsWOHncYrDgI4WBs2bDDjx48PvX/ixIlm0qRJNW0TAAAAkFVNTU1m2rRpgfctW7bMdO3aNdZy2zQ3NzdX2DZk3JAhQ2wGbubMmaGPIQMHAADSfCF6xowZZsWKFaZ3797mlFNOMV26dKl3s9DgtkZk4EaOHGkzcPPnzy97uWTgYGkH6tatW72bAQAAUDYFb2+++ab9Xbf6++yzz653s9DgOkYkQNq1ix+GUcQEAAAAmabMW9TfQJ4QwAEAACDT1G0y6m8gTwjgAAAAkGka89a3b187JES3+hvIK8bAAQAAINNUsIQxb2gUZOAAAAAAICMI4AAAAAAgIwjgUBejR4+2PwBqj+MPqC+OQaC+Rmf8GCSAAwAAAICMoIgJrE8//dSsXbs21kSEAAAAAFraunWr/QmyY8cOWzU1DgI4WBs2bDDjx48PvX/ixIlm0qRJNW0TAAAAkFVNTU1m2rRpgfctW7bMdO3aNdZyCeBQKL87ffr00PvJvgEAAAClGzdunBkzZkzgfSNHjiQDh8poB+rWrVu9mwEAAADkQseIIUjt2sUPwyhiAgAAAAAZQQAHAAAAABlBAIdIqpwzderU0Ao6cR+L7Gy3WrevWq+X1HIrWU6c53IMVl/atxvHYHLL4BhMp7RvN47BZJZR7eMvC/tSUgjgEEkHgKrnlPrFVepjkZ3tVuv2Vev1klpuJcuJ81yOwepL+3bjGExuGRyD6ZT27cYxmMwyqn38ZWFfSgoBHAAAAABkBAEcAAAAAGQEARwAAAAAZAQBHAAAAABkBAEcAAAAAGREm+bm5uZ6NwL11bVrV/PRRx+Zww47rNV9O3bsMMuWLTMHHHBA0Rnjy3ns0qVL7W2/fv1MoytnuzVC+6r1ekktt5LlxHluNY5Bjr+WOAZr83pJLLfSZXAMphPHYGMcg9U+/rJ2DL766qumQ4cOZv369WU/N31HCWpuzz33DL1vp512sgGebosp57F8acXbbo3Qvmq9XlLLrWQ5cZ5bjWOQ468ljsHavF4Sy610GRyD6cQx2BjHYLWPv6wdgwreos7Bo5CBAwAAAICMSOelDgAAAABAKwRwAAAAAJARBHAAAAAAkBEEcAAAAACQEVShRKYsXrzY/P73v7dVezT1Qc+ePc0ZZ5xh2rdvX++mAQ3n+eeft2WQJ02aVO+mAA3h/fffN9OnT7e3q1evNhs3bjTDhg0zEyZMMJ06dap384DcW7p0qXnooYdMmzZt7HG4YsUKM3jwYHP66afHrigZB1UokRnz5s0zV111lbn55ptN9+7d7f9uueUWs3LlSnPFFVeYtm3b1ruJQMPQyeMll1xijjnmGAI4oAa2bdtmfvKTn5hzzjnHXryURYsWmYsvvtjstdde5sYbbzSdO3eudzOB3Hr//fft992vf/1rs++++9r/vfnmm+b73/++6d27tz0/rdW5KF0okQkffvihufrqq82RRx5ZCN7km9/8ppkzZ46ZMWNGXdsHNIq77rrL/PM//7P50Y9+ZC+eAKiNp556ygwZMqQQvEn//v3NcccdZ9555x1z991317V9QN6tXbvWbN261Vx00UVmw4YN9n99+/a1GbjXX3/d/tQKXShRMV31U9eNcq7Cv/zyy+bOO++0B4C7WqEr+SeeeGLg1YsHH3zQPvbQQw9t8X8Fc5qwUcs69dRTUzsJKJCH41DUZVk/CxYsMD/84Q8Taz+QZbU4/tavX28efvhhs8cee5gRI0YUHnPAAQfY21qePAKNeAx+7nOfM//wD/9gVq1aZXbdddfC/z/99FN7/tmlSxdTKwRwqMiaNWvMo48+aoOnUilbduutt5rLLrvMDB061P5v3bp15sc//rF55plnzJQpU1qNaZs1a5a97dOnT6vl6X8vvviieeWVV+zVSKDR1Oo4BFC/409X+fUc9UZp166dGT58uP2/MgLSo0ePRNcLyIo1NfwOPOGEE1r8rQsn6sqsC5p77723qRXSFYhFQye1w06ePNls3ry55Oc9/fTTdtzamWeeWThgRP33NZ7mpZdesl9O/qGZy5cvt7e77bZbq2W6/7nHAI2i1schgPodfwcddJD5+c9/bs4991zbldJ7gVPFE5QZBxpJcx2+AzXmTceaXvOnP/2pfbyGFIwePdrUEhk4lKWpqcl+WWgg55YtW2x/4FLp8ddff709II4//vhW9++zzz5m0KBBtp//I488YkaNGlW4zx2Yu+yyS6vndejQwd5+8MEHMdcKyJZ6HYcA6nv8HXbYYfbH+eMf/2hfXyeStbz6DzTqMdi3b19bxMQVDHrttdfMhRdeaJd13nnnmVohgENZxo0bZ39Eg6YnTpxY8nM1wFpBlnb+sEpZAwYMMHPnzjV33HGH7Yfs+iDvvvvutlyypg7wc//TuACgEdTrOASQnuNPY3c0hkeZhFqWLwca/Rjs7HneIYccYkaOHGnuuecee3FFv9cCXShRM/fff7+91UETxg3G1gBRjWvz/19XTvzc/9xjAFTnOARQ/+Nv+/bttiulxudceeWVBG9ADY9BJQ10DHrtv//+9lZZu1ohgENNLFu2zA4OFTd3RhDvIOwXXnih8LumDxBNmOj31ltv2UqU3m4lAJI/DgHU9/j7+OOPzQ033GDGjh1rTjvttBb3zZw5M/E2A3myrMJjUPMRa/qqs846q0WPMFWhlFpWQieAQ00sXLiw8HvURKPe+zQw1VHfYgVpr776aqvKQyqtrAGlTCEAVPc4BFDf4+83v/mNrT7Zr1+/Fv/XnIzqSgagesegMmy6iKIqsG3atCn8f8mSJfb2iCOOMLXCGDjUhLviId65M/w6duwY+hxV/FFZV31JuSsnGsj65S9/2Zx88slVazuQF5UehwDqd/x9+OGH5k9/+pN5++23C1PryCeffGIzAxQcAqp7DB5++OG2qJ4qTzqawkqBnc5FjzvuOFMrBHCoCVUKcnTlIoz3PhUt8VJVoJ/85Ce2+o+Kmih93bt3bzuHB0UWgNoch/fdd5+92ui6M6tK13vvvWerwaq8OYDqHH8aLqCxN88991wVWwnk1/sVHoOaKkDdJdWNWRdOdJ8Konzve98zY8aMqem5KAEcasJb/t/1FQ7ivc9ND+B16KGHmssvv7wKLQTyL4nj8MQTT6xS64B8q/T4U7U7xrkB8SXxHagsWy0zbWEYNISa6NKlS+F3f/UeL/UtLqV/MoDycRwC9cPxB9RXlxwdgwRwqIlevXoVft+0aVPo47z3qXskgORwHAL1w/EH1FevHB2DBHCoCY1fc32Do4oieO8bMmRITdoGNAqOQ6B+OP6A+hqUo2OQAA41sdtuu5mBAwcWJkYM4wojaEqAYcOG1ax9QCPgOATqh+MPqK/dcnQMEsChZiZMmGBvly5dGjp4dPHixYVKP927d69p+4BGwHEI1A/HH1BfE3JyDBLAoWZ01UNXMrZs2dJiMkVHJVnnzJljr5BolnsAyeM4BOqH4w+or4E5OQYJ4BDbtm3bAiv2RNFk3Pvvv7+ZNm2aPUi8HnjgAbNhwwY7r5ubqBtANI5DoH44/oD62tagx2Cb5ubm5no3AtnR1NRkZs2aZX9fu3atWb16dWGejH79+pk2bdrYCRCvueaa0GWsX7/eXHvttfZAGzt2rJ3xfv78+Xa5P/jBD8zQoUNrtj5AFnEcAvXD8QfUVxPHIAEc6mfJkiVm0aJF9uDp06ePOeKII0z79u3r3SygoXAcAvXD8QfU15KMHoMEcAAAAACQEYyBAwAAAICMIIADAAAAgIwggAMAAACAjCCAAwAAAICMIIADAAAAgIwggAMAAACAjCCAAwAAAICMIIADAAAAgIwggAMAAACAjCCAAwAAAICMIIADAAAAgIwggAMAAACAjCCAAwAAVderVy/ToUMH+9OmTZvCz4gRI+rdNADIFAI4AABy5sorrzSdO3c2hxxyiJk1a5ZJg5UrV5qPPvrI/txxxx31bg4AZBYBHAAAObJkyRJz2WWXmffff9/+fuGFF5q06dmzZ72bAACZRQAHAECO7LRTy692dVMEAOQHARwAADly0EEHmSlTpphOnTqZgw8+2PzsZz+rd5MAAAkigAMAIGf++Z//2WzcuNG89tpr5ktf+lK9mwMASBABHAAAAABkBAEcAAAAAGQEARwAAAAAZAQBHAAAAABkBAEcAKChDBkyxHTo0MG0bdvWlth3P5s3b7b333XXXeYrX/mKnQhbP2PGjDEPPfRQWa9xzz33mBNOOMF0797d7LLLLmafffYxJ554ornvvvtCn6P71K527dq1aNfy5cvt/WrfRRddZPr06WN23313M3ToUHPvvfe2WEavXr3s63mf/6//+q8ltXnbtm3mlltuMcOHD7frveuuu5oDDzzQ/N//+3/NwoULTbmamprMsccea/bcc0/TrVs3c9xxx6VmUnEAyDICOABAQ5k7d6756KOPzKZNm2yg4vX973/f/P73vze33XabWbt2rfl//+//mbfeessGH6VMiP3hhx/aQO2UU04xhx12mHnppZfMli1bzG9/+1vzwgsvmJNOOsl8+9vfNp9++mmr5+p5apeWceihh7a4TxUljz76aBsIvfzyy2bVqlX2f3qdRx55pPC4lStX2kDsiiuuKGubLF261AwcOND86Ec/MmeddZYNGj/44ANz991329caNGiQufLKK0ta1o4dO8yECRPM3/zN35hPPvnEPPvss+add94xV111la2O+Z//+Z9ltQ0A4NMMAECD6tu3b7O+CvVzzTXXNI8fP775008/bfGYVatWNe+xxx72MRdddFHosrZs2dI8YsQI+7j/7//7/1rd/9hjjxVe69xzz41s1/DhwwuPXbZsWfOxxx7b/Ktf/apw/0MPPdTcoUMHe//EiRNbPf/2228vPP/yyy+PfK0FCxY0d+/evblNmzbNjz/+eKv7d+zYUViv73//+83FfOtb37KP1bbdvHlzi/u2bdtm1+3f/u3fCu3T3wCA0pGBAwDAGJshuvnmm223Q6/99tvPZubk3//9321GKcgll1xiHn/8cbPHHnsEdlscMWKEOfzww+3vv/zlL828efNKapcygrvttpv57ne/a/9WNkuZPGXrRBm+uJStO/30082aNWtsBlDdJ/3U1fTaa68ttPu//uu/Qpc3ffp0c+edd9rfJ0+ebNvt1b59e/Mf//Ef5uqrr47dZgBodARwAAAYYyZOnGi6dOkSeJ+CHOeCCy5odf+iRYvMTTfdZH8/9dRTbRAXRF0xpbm52UyZMqWkdikYvPzyywt/q2vjxx9/bH/faaedzDnnnGPi+slPfmIn+xZ1ewyj8XYHHXSQ/V1dSd97773AYPDiiy+2vysIVhfKIAcffLDtDgoAiIcADgAAY8zXv/710Ps0Pqxjx4729+eff9688sorLe5XZkpjv1ymLcxXv/rVwu8zZ84027dvL9quL3zhC7bwitOvXz87Nk8ZP41/Gz16tIlDRVF+8YtfBLYtyFFHHWVv169fb4ud+D388MN2DJ4L0vbaa6/QZblAFgBQPgI4AACMsUVHwqgboSoyOn/84x8Lv6sgiYp9OApewuy9994tMmkqclLMqFGjWv1PQZuyZ1HBYjEPPvigbYOoqqUqWJa6fdSt08+7TaK2gfiLtAAASkcABwCAMZEZI+nZs2eLSpaOuiC+++67hb/DumEGvcbq1auLtuvII4801eAdy+ddtzDeAE+VMN9///0W93u3SbHlqZomACCedjGfBwBArmj+tSjeghxvv/124XeV2fdSyX1/IRR/Nk8055t+itl3331NNbzxxhuF35WBK8b/GD1f6xq0TfzFS8rd1gCAcHyCAgBQAk2Q7WgcmKMKjl4LFixItItg165dTTW4ictddchy1j+o+qV3m5SyPABAPHShBACgBKqyGBRUueImjibiTpIqTVaDN6PmqlqWuv5BWTbvNilleQCAeAjgAAAwplBFMszWrVsLv/fo0SN0PFdQif008hZl8Wbjwvgf432+v6tnsbnpim1rAEA4AjgAAIwx69ati7zfO9bNW9Z/wIABLbJky5YtK3l+txtuuMHUy1e+8pXQcXxB3BQB8vnPf9506tSpxf3ebVKsOMvatWvLbC0AwCGAAwDAGPPqq6+G3vfJJ5+0KPqhybqdPffc0wwePDiwGmOUxx9/3Hzuc58z9XL88ccXJhxXds0boBXbPkGTdJ988smF393k4GEWL14co8UAACGAAwDAGPPQQw+F3qf52ly3wKFDh5rDDz+8xf3nnHNOiwmtm5ubI19LwdK8efPMMcccY+pFwdv5559f+PuZZ56JfPzTTz9dmCbhe9/7XuBE6K5r6ZIlS1oUNfH785//XEHLAaCxEcABAGCMaWpqMhs2bAi8z01crekBfvGLX7S6f+LEibYrpetC+eijj0a+1r/8y7+Yv//7v29VAKXWLr30UtOvXz/7+/Tp00Mf9/zzz5ulS5fa36+77jrTvXv3Vo/RlAhXX321/V0BrLZnEAV3BHAAEB8BHAAAxpjzzjvP/vhpfrNf/vKX9vcf/ehHgRNrq2z+XXfdZfbee2/79z/8wz+EFga55pprbDbr8ssvN/WmoEvBqdp97733mieeeCKw+6jWW7773e+as88+O3R5CmS/+c1v2t8VzHkLv8j27dvtMv7t3/4t8XUBgEZBAAcAgDFmwoQJplevXub000+3Y7hUKXH27Nm2a+CmTZvMD3/4Q3PVVVeFPl9zvz377LPmi1/8onn55ZfN0UcfbbtTaloBBS4K2k444QRz66232u6aGjsXRI999913Sy4IEsQ7nq3Y89Xep556yo7HO/HEE81vfvMbs3HjRhu4aU47jW1TYHfZZZeZX/3qV5HLUoZSmbzTTjvNZiLHjBljt4W25fz58+24O/1Pt47m0dM6AwBK06a5WEd9AAByav/99zdvvvmm/V0Bh/6eMWOG+fd//3cbeCgg+epXv2r+8R//0YwcObKkZX766afmf/7nf+yPuh5qWoGdd97ZHHzwwTY4VJZv1113bfW8++67zxYHUbCj4MlRG5ThU1D0u9/9LvK1FYDq9fzzsOn5Ws7777/fakJu56OPPrLBmzKJGvOned80pk3rrbFyAwcONOVQIKeAT0Vd9JqDBg0ykydPNqNGjTKLFi0yX/jCFwqPbdu2rWnXrp0N7u6+++6yXgcAGg0BHACgYQUFcAAApBldKAEAAAAgIwjgAAAAACAjCOAAAAAAICMI4AAAAAAgIwjgAAAAACAjCOAAAAAAICMI4AAADWnLli12XjTn7bffrmt7AAAoBfPAAQAaypAhQ+wk3du3b7eTbvsnzHaTWgMAkEYEcAAAAACQEXShBAAAAICMIIADAAAAgIwggAMAAACAjCCAAwAAAICMIIADAAAAgIwggAMAAACAjCCAAwAAAICMIIADAAAAgIwggAMAAACAjCCAAwAAAICMIIADAAAAgIwggAMAAACAjCCAAwAAAACTDf8/nfzsFr9u0zcAAAAASUVORK5CYII=", | |
| "text/plain": [ | |
| "<Figure size 432x288 with 1 Axes>" | |
| ] | |
| }, | |
| "metadata": { | |
| "image/png": { | |
| "height": 296, | |
| "width": 440 | |
| } | |
| }, | |
| "output_type": "display_data" | |
| } | |
| ], | |
| "source": [ | |
| "plt.figure(figsize=(6, 4))\n", | |
| "plt.loglog(\n", | |
| " samples[\"P\"].value,\n", | |
| " np.abs(samples[\"K\"].value),\n", | |
| " alpha=0.5,\n", | |
| " marker=\".\",\n", | |
| " linestyle=\"none\",\n", | |
| ")\n", | |
| "plt.xlabel(\"period\")\n", | |
| "plt.ylabel(\"$K$\")" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 86, | |
| "id": "b2121adf", | |
| "metadata": {}, | |
| "outputs": [ | |
| { | |
| "data": { | |
| "text/plain": [ | |
| "np.float64(0.8994140625)" | |
| ] | |
| }, | |
| "execution_count": 86, | |
| "metadata": {}, | |
| "output_type": "execute_result" | |
| } | |
| ], | |
| "source": [ | |
| "P_min = 1.0 * u.day\n", | |
| "K_min = 1.0 * u.km / u.s\n", | |
| "\n", | |
| "mask = (samples[\"P\"] > P_min) & (np.abs(samples[\"K\"]) > K_min)\n", | |
| "mask.sum() / len(mask)" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": null, | |
| "id": "8ec2f8f6", | |
| "metadata": {}, | |
| "outputs": [], | |
| "source": [] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": null, | |
| "id": "a848ad35", | |
| "metadata": {}, | |
| "outputs": [], | |
| "source": [] | |
| } | |
| ], | |
| "metadata": { | |
| "kernelspec": { | |
| "display_name": "demo-notebooks (3.12.10)", | |
| "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.10" | |
| } | |
| }, | |
| "nbformat": 4, | |
| "nbformat_minor": 5 | |
| } |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment