{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "f07fa1f2-187e-4ce0-af95-31d6120977fe", "metadata": { "ExecuteTime": { "end_time": "2024-11-20T18:02:26.764056Z", "start_time": "2024-11-20T18:02:16.122104Z" }, "pycharm": { "name": "#%%\n" } }, "outputs": [], "source": [ "from mimic.utilities.utilities import set_all_seeds\n", "from mimic.utilities.utilities import plot_CRM\n", "from mimic.model_infer import *\n", "from mimic.model_simulate import *\n", "\n", "import pandas as pd\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from scipy.integrate import odeint\n", "\n", "import pymc as pm\n", "import arviz as az\n", "from pymc.ode import DifferentialEquation\n", "from mimic.model_simulate.sim_CRM import sim_CRM\n", "import pytensor.tensor as at\n" ] }, { "cell_type": "markdown", "id": "82eb9f01", "metadata": { "pycharm": { "name": "#%% md\n" } }, "source": [ "## Simulate some time course data from the CRM\n", "\n", "The McArthur Consumer Resource Model takes the form\n", "\n", "\n", "$$ dN_i/dt = 1/\\tau_i N_i (\\sum_a(c_{ia} w_a R_a - m_i)) $$\n", "$$ dR_a/dt = 1/(r_a K_a) (K_a - R_a) R_a - \\sum_i(N_i c_{ia} R_a) $$\n", "\n", "where:\n", "- $N_i$ is the concentration of a species\n", "- $R_a$ is the concentration of a resource\n", "- $c_{ia}$ is the preference of species $i$ for resource $a$\n", "- $w_a$ is the quality of resource $a$\n", "- $m_i$ is the mortality rate of species $i$\n", "- $K_a$ is the carrying capacity of resource $a$\n", "- $\\tau_i$ is the timescale of species $i$\n", "- $r_a$ is the timescale of resource $a$" ] }, { "cell_type": "markdown", "id": "fbc11bbc", "metadata": { "pycharm": { "name": "#%% md\n" } }, "source": [ "### Generate parameters for model with two species" ] }, { "cell_type": "code", "execution_count": 3, "id": "3e0845a5", "metadata": { "ExecuteTime": { "end_time": "2024-11-20T18:02:29.987899Z", "start_time": "2024-11-20T18:02:29.970540Z" }, "collapsed": false, "pycharm": { "name": "#%%\n" } }, "outputs": [], "source": [ "set_all_seeds(123)\n", "\n", "num_species = 2\n", "num_resources = 2\n", "times = np.arange(0, 10, 0.1)\n", "\n", "\n", "tau = np.random.uniform(0.1, 0.9, num_species) # species timescales\n", "w = np.random.uniform(0.1, 0.9, num_resources) # resource quality\n", "c = np.random.uniform(0.1, 0.9, (num_species, num_resources)) # relative resource preferences\n", "m = np.random.uniform(0.3, 0.7, num_species) # mortality rates\n", "r = np.random.uniform(0.1, 0.9, num_resources) # resource timescales\n", "K = np.random.uniform(1.0, 10.0, num_resources) # resource carrying capacities\n", "\n", "\n", "# write the parameters to a dictionary and pickle\n", "params = {'num_species': num_species, 'num_resources': num_resources, 'tau': tau, 'w': w, 'c': c, 'm': m, 'r': r, 'K': K}\n", "pd.to_pickle(params, 'params-s2-r2.pkl')\n" ] }, { "cell_type": "code", "execution_count": 4, "id": "44dae2c8", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "{'num_species': 2, 'num_resources': 2, 'tau': array([0.65717535, 0.32891147]), 'w': array([0.28148116, 0.54105182]), 'c': array([[0.67557518, 0.43848517],\n", " [0.88461136, 0.64786379]]), 'm': array([0.49237276, 0.45684701]), 'r': array([0.37454241, 0.68323977]), 'K': array([4.9471502 , 1.53710107])}\n" ] } ], "source": [ "print(params)" ] }, { "cell_type": "markdown", "id": "09cde050", "metadata": {}, "source": [ "### Simulate single time course" ] }, { "cell_type": "code", "execution_count": 5, "id": "0ead4c25", "metadata": { "ExecuteTime": { "end_time": "2024-11-20T18:02:33.149100Z", "start_time": "2024-11-20T18:02:32.677261Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA9wAAAJOCAYAAABFiQ/hAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAD4MElEQVR4nOzdd3wUdf7H8dfuppAeAoEAoYMQ0NBRinQEO4ooVhBPT0WsqAcqYDl7Oe9QTk+KDVQsyE9sgCCIjV4UkN47JAHSd+f3x5AlIQmTZJNMsnk/H499ZHd2dvYz2Xc2+9mZ+Y7DMAwDERERERERESlVTrsLEBEREREREfFHarhFREREREREyoAabhEREREREZEyoIZbREREREREpAyo4RYREREREREpA2q4RURERERERMqAGm4RERERERGRMqCGW0RERERERKQMqOEWERERERERKQNquEVECjB8+HAcDgfbt2+3uxQppoULF+JwOJgwYYLdpXhVtjxVtnqLo1GjRjRq1MjuMio0f379RUTKmxpukXK2fPlybrvtNpo3b05YWBghISE0bdqUm2++mblz5+aZd8KECTgcjjyX0NBQzj33XB577DFSUlIKfI6ceYODgzly5EiB8xw7doyQkBDvvMW1ZcsWHnzwQdq2bUv16tUJDAwkNjaWHj168NRTT7Fjx45iL7M82dGU/fXXX4waNYrWrVsTGRlJcHAw9evX55prruGzzz7D4/GUWy2VzZIlSxgyZAj16tUjKCiI6tWr07JlS2644Qbeffddu8srF9OmTcPhcDBt2jS7S6kwFixYwHXXXUf9+vUJDg4mJiaG7t2789prr5Genm53eSIiIgTYXYBIVeHxeBg9ejSvvfYaAQEB9OnThyuuuILAwEC2bt3KnDlz+OCDD3jqqad44okn8jx28ODBnHvuuQAcOHCAr7/+mmeffZavvvqK33//neDg4HzPFxAQQGZmJh9++CH33ntvvvs//PBD0tPTCQgIIDs7u1jr8uqrr/Loo4+SnZ3NBRdcwE033URkZCRHjx5l2bJlTJgwgWeeeYYlS5bQqVOnYi3bX73yyis8+uijeDweunfvTv/+/QkNDWXXrl3MmzePzz77jBEjRjB58mS7S61wpk2bxogRIwgICOCSSy6hefPmOBwONm7cyNdff82iRYsYNmyYd/7OnTuzfv16atasaWPVUpays7MZOXIkb7/9NmFhYVx88cU0a9aM5ORkvv/+ex588EH++9//MmfOHJo1a2Z3uSIiUoWp4RYpJ48//jivvfYabdu25dNPP6Vp06Z57k9LS2PixIkFbpG+5pprGDp0qPd2eno6F1xwAatXr2b69Onceuut+R7TtGlTDMNg6tSpBTbcU6ZMoUWLFgBs3LixyOvx1ltv8dBDD9G4cWM++eQTOnbsmG+ev/76i3HjxhW6Bb6qefvttxk9ejSNGjXis88+o3379nnuz87O5t1332Xx4sU2VVhxpaamcu+99xIREcHPP/9M69at89yflZXFwoUL80wLDQ2lZcuW5VillLcxY8bw9ttv06lTJ7744gvq1avnvc/tdvPUU0/x1FNPMXDgQFasWEFkZKSN1YqISJVmiEiZ27Rpk+FyuYwaNWoY+/fvP+u86enp3uvjx483AGPGjBn55nvxxRcNwBg5cmS++wCjRYsWxnPPPWcAxvLly/Pcv2rVKgMwnn/+eaNFixZGUd8Kjh49akRGRhrBwcHGhg0bLOfPysrKV1fPnj2N3bt3GzfffLNRu3Ztw+FwGAsWLPDOM2XKFKNz585GWFiYERYWZnTu3NmYOnVqvjqcTqdx6aWX5pm+cuVKAzAAY9OmTXnu69mzp1GtWjUjPT3d+3st6LJt2zbDMAxj2LBhBmBs3brVeP31140WLVoYQUFBRoMGDYwJEyYYbre7SL+zY8eOGZGRkUZQUJDxxx9/nHXe3K99zvPn1JNbTv25f28LFiwwAGP8+PHGkiVLjP79+xtRUVEGYGzfvt1wOBxG7969C3zezMxMo0aNGkZ8fHye9crIyDBeeeUVo127dkZoaKgRHh5udO/e3fjyyy/zLSMpKcl44oknjISEBCMsLMyIiIgwmjZtatxyyy3G9u3bLX5Lhfvtt98MwLjyyiuL/Jjcv4vcGjZsaDRs2NBISkoy7rzzTiMuLs4IDQ01LrzwQu/fyJ49e4wbb7zRiI2NNapVq2b079/f+Ouvv/IsZ9u2bQZgDBs2rMDnz8l5bgW9nhkZGca///1v46KLLjLi4+ONoKAgIzY21rjqqquMFStWFPj4gi65paSkGOPGjTNatWplVKtWzYiKijIuuugiY/HixQXWum7dOuPSSy81wsPDjcjISOPiiy821q5de9b85VaSbPmalY0bNxpOp9OIiYk56/vpDTfcYADGE088kWd6Tg6OHTtm3HHHHUbt2rWN4OBgo23btsb06dPzLSctLc14+eWXjcTERCMyMtIIDQ01GjZsaAwZMsRYtWpVvvlnzZpl9OnTx4iOjjaCg4ON1q1bGy+99JKRnZ2dZ76pU6cagDF16lRj9uzZRteuXY3w8HCjYcOGxqJFiwzAuPXWWwtctwMHDhgBAQFG165d80wv79dfRESsaQu3SDmYNm0abrebv//979SuXfus8xa0e/jZBAQU/mc8bNgwHn/8caZOnZpnq+rkyZNxuVzccsstTJ06tcjP9emnn5KSksJNN93k3Tpe3NqOHDlCly5diImJYejQoaSnp3u3Pt1777385z//oV69etx2220AfPbZZ9x6662sXLmS119/HYDq1avTpk0bFi9ejNvtxuVyAebxnDkWLFjg3ZU0PT2dX3/9la5duxIcHEyvXr3Yvn077777Lj179qRXr17ex0VHR+ep9+GHH+bHH3/ksssuY8CAAcyaNYsJEyaQmZnJP//5zyL/zm644QZatWp11nmL+9oX5Oeff+bZZ5+ld+/e3HHHHezcuZOGDRvSo0cPfvzxR3bv3k18fHyex3z99dccOXKERx99FKfTHNojIyODgQMHsnDhQtq2bcttt91GVlYWc+bM4corr+Q///kP99xzDwCGYTBgwAB+++03unXrxsCBA3E6nezYsYPZs2dz880307BhQ+/z5YwZYBiG5frUqFEDgK1bt+Z5rUsqMzOT/v37k56eznXXXceBAwf45JNP6NevHz///DMDBgygTp063HTTTWzevJn/+7//49JLL2X9+vU+P/eZjh49yv3338+FF17IJZdcQvXq1dm6dSuzZ8/mm2++YdGiRd5DMgYNGkRSUhJffvklV155JW3bti1weT169OCPP/6gW7du3HnnnaSkpPDll1/Su3dvZs6cyaBBg7zzr1u3jm7dunHixAmuvvpqmjdvzu+//063bt1o06ZNkdahuNkqblYK8u677+LxeLjjjjvO+n76xBNPMH36dKZMmcJTTz2V577MzEz69evHiRMnuPnmmzl58iSffPIJN9xwA4cPH2bUqFHeeYcNG8Ynn3xCYmIit956K8HBwezatYsFCxawdOnSPL+rMWPG8Pzzz1OvXj2uvvpqoqKiWLx4MQ8//DC//fYbM2fOzFfnzJkz+f7777nsssu4++67SUlJoXv37t49Yt58802qVauW5zEzZswgOzubm2++2TvNjtdfRESKwOaGX6RK6NWrlwEY8+bNK9bjCtvCnZaWZrRp08YAjJkzZ+Z7HKe2cBuGYVx22WVGTEyMd+tpenq6ERMTY1x++eWGYRjF2sJ96623GoAxefLkYq1H7ro4tdXmzK09P/74owEYCQkJRlJSknf60aNHjXPOOccAjEWLFnmnP/jggwZg/Pbbb95pl19+uXHOOecY9evXN66//nrv9Pnz5xuA8dRTT3mnFbYVNEfOFp7GjRsbe/fu9U4/dOiQER0dbURERBgZGRmW6zx8+HADMN555x3LeQt6/uJu4QaMKVOm5HvMO++8YwDGCy+8kO++wYMHG4Cxbt0677SxY8d6tw56PB7v9JSUFKNjx45GUFCQsWfPHsMwDGPNmjUGYAwaNCjfstPT043jx4/nmZZTZ1F4PB6jQ4cOBmB0797d+N///mesXbs2X35yO9sWbsAYMmRInr0vXnjhBQMwoqOjjQceeCDP+t51110GYHz22WfeaaW1hTs9Pd3YvXt3vsevW7fOCA8PN/r165dneu4togXJ2aL7v//9L8/0AwcOGPXr1zdiY2ONtLQ07/SePXsagPHBBx/kmX/MmDH59vg4m+Jkq7hZKUjO++ncuXMt561bt64BGDt37vROy8lBjx498vwN79q1y6hZs6YRHBzsfV2SkpIMh8NhdOjQIV/msrOzjWPHjnlvf//99wZgDBgwwDhx4oR3usfjMe68804DMD799FPv9JzX0+l0Frgujz/+uAEYH3/8cb77OnToYAQFBRlHjhzxTrPr9RcRkbPTKOUi5WD//v0A+bb+FNWnn37KhAkTmDBhAnfffTctWrRg9erVXHXVVVx99dVnfeyIESM4evQos2bNAmDWrFkcPXqUESNGFLuOnPWoW7duvvtWrVrlrTHnkvOcuQUFBfHiiy/m21qYM9L0hAkTiIqK8k6vXr0648ePB8gzOnPv3r0B+OGHHwDzuM1FixbRu3dvevfunW9rN5BnS3ZRPfHEE9SpU8d7u2bNmlx55ZUcP368SMe++/raF1f79u0LPKb/mmuuoVq1anzwwQd5piclJfHVV1/Rtm1b7/HRHo+HSZMm0bRpU5588sk8o9hHREQwbtw4MjMz+fzzz/MsKyQkJN/zBgcHEx4enmfa+vXrWb9+fZHWx+Fw8Omnn9KtWzd++uknbr/9ds477zwiIyPp16+fd++R4nj55Zfz7H1x/fXXA+ax9M8880ye9c25b/Xq1cV6jqIIDg7Oc+xxjtatW9O7d28WLVpEVlZWkZZ1+PBhPv74Y/r06cPf/va3PPfVqlWLhx9+mEOHDjFv3jwAdu7cyY8//khiYiI33nhjnvnHjh2bb0+PsylOtnIUNSsFyfmbql+/vuW8OfPs27cv333PPvssQUFB3tvx8fHcd999ZGRk8NFHHwFm/gzDoFq1at69P3K4XK48v6eJEycCeAdyy+FwOHj++edxOBzMmDEjXx1XXnkl/fr1yzc9Z+v1mb/X9evXs3z5ci655BJiYmIAe19/ERE5O+1SLlIJfPbZZ3z22Wd5pg0ZMoSPP/7Y8pRel112GbVq1WLKlClcd911TJkyhVq1anHZZZeVao2rVq3iySefzDNt2LBheXZhBGjcuHGBo0evXLkSKLgpzmmuV61a5Z3Wo0cPXC4XCxYs4B//+AcrV64kOTmZPn36kJqaynvvvcf69etJSEhgwYIFhISEcP755xd7vTp06JBvWk7znJSUVOzllbXCRoWPioriiiuu4JNPPmH16tXeXUZnzpxJRkZGnl1TN27cyLFjx6hbt26+1xTg0KFDAGzYsAGAhIQEEhMTmTFjBrt372bQoEH06tWLtm3b5mtSgGIPaNaoUSN++uknVq1axbx581i2bBlLlixh/vz5zJ8/n/fee49vvvmmSLvkV69enQYNGuSZlvOFSvPmzQkNDS3wvr179xar5qJatWoVL774Ij/99BP79+/P12AfPnw4zxc+hVm6dClut5uMjIwCT3W3adMmwHzNLrvsMu8XCN27d883b3h4OG3bts03GF1hipOt4malrAQEBNClS5d80y+88ELg9PtRZGQkl1xyCV9//TXt27dnyJAh9OrVi06dOhEYGJjnsb/++ithYWFMmTKlwOcMCQnx/s3k1rlz5wLnP+ecc+jcuTPffvsthw8f9r5v5jTguX+vdr7+IiJydmq4RcpBXFwcGzZsYM+ePUU69vlMM2bMYOjQoWRnZ7Nx40ZGjx7NzJkzadGiBU8//fRZHxsYGMhNN93Ev/71L37++WfmzZvHAw88cNZjvwuTc7xkQc3H8OHDGT58OGB+8Czow2zuZZwpJSUFp9NJbGxsgY9xOBx5Rj2PjIykffv2LFmyhKysLBYsWIDD4aB3796kpqYC5pbthg0b8vvvv9OzZ888W7OKqqDRjXN+d0XZshoXFwfAnj17iv3cJXG2Y1pvvvlmPvnkEz744ANvU/T+++/jcrm44YYbvPMdPXoUgD/++IM//vij0OWdPHkSMH8fP/zwAxMmTOCzzz7joYceAiA2NpZ77rmHxx57rFSOf27btm2eY5cXLlzITTfdxIIFC3jzzTd54IEHLJdxttfzbPcVdUtzcfz888/06dMHgIsuuojmzZsTHh6Ow+Fg1qxZrF69moyMjCItK+c1W7JkCUuWLCl0vpzXLDk5GTC3fhbEaqyJMxU1W6WRlZz30127dlm+n+7atQsg35cWNWvWLLDBz1nvnN8PmF8cPPvss0yfPp3HHnsMMLNy66238uyzz3q/pDl69CjZ2dkFfkmVI+f3X9BzFuTmm2/m999/5+OPP2bkyJEYhsGHH35I9erVufTSS73z2f36i4hI4bRLuUg56NatGwDz58/3aTkBAQG0bt2aL774gmbNmvHPf/6TFStWWD7utttuw+PxcO211+LxeLwDkhVX165dgbyDkxVXYVvkIyMj8Xg83q2nuR08eBDDMPI1RL179+bkyZP8/vvvLFy4kNatWxMbG0vDhg1p3LgxCxYs8DbkOVvJy1tJX/ucZqCgc6TnbgbOdLY9HgYOHEhsbCwzZszA4/Gwfft2fvrpJ/r16+f9YgBON56DBw/GMIxCL7kH3KtRowb/+c9/2LNnD3/++ScTJ04kJiaG8ePH8+KLLxZr3YuqV69e3i+ccg4tKA8lfW3O9M9//pOMjAzmzZvH7NmzeeWVV3jyySeZMGFCntejKHJes4ceeuisr1nO4Rk5h20cPHiwwOUdOHCgWM9f1GyB71nJeR+y+pvasGEDe/fupV69evl2Pz98+DAejyffY3LWO/dhLaGhoTzzzDNs3bqVrVu3MnnyZFq0aMHrr7+e50ueyMhIatSocdbf/7Zt2/I959n+ZocOHUpgYKB3q/aiRYvYsWMH1157bZ49Oux+/UVEpHBquEXKwfDhw3G5XLz99tsFNpS5FWWLVrVq1Xj55ZcxDIN//OMflvO3atWK888/nz179nDBBReQkJBQ5Npzu+aaa4iIiGDmzJneXRRLS7t27QAK3I0xZ9qZIzPnNNHff/89ixcv9m4tBOjTpw8LFy70NmJn7qqesxWtuMf/Ftc111xDZGQkn332WYG7k+aW+7WvXr06UPCW8ZzdXYsrICCAoUOHsmfPHhYsWMCHH36IYRjcdNNNeeZLSEggMjKSZcuWFXvLrsPhICEhgZEjRzJ37lwAZs+eXaJ6i6Iox/yWtpzjW319bbZs2UJMTEy+3XpTU1ML/CLtbJnt1KkTDoeDX375pUjPnbMV+qeffsp334kTJ/IcvlEURc1WbiXNyi233ILT6eR///vfWd9Pc84iUNB4FdnZ2QX+rhYvXgycfj86U+PGjRkxYgQ//vgj4eHheeo9//zzOXLkSKm+N9asWZOBAwfy66+/snnzZm/jfebv1e7XX0RECqeGW6QcNGvWjEceeYTDhw9z8cUXF7iVIz09nVdffbXA4+8KcuWVV9K+fXvmzp3r/ZB4NlOmTOGLL75g8uTJxS3fq3r16rz00ktkZGRw8cUXs3z58gLnK8mxzcOGDQPgySefzLPreHJysncXzZx5cnTv3p2AgAAmTZrE8ePH8zTcvXv35vDhw0yePJmwsLB8xzbnDDaUs8tpWYmOjvb+zi699NICP8i63W7effdd7rzzTu+0nHpzDxQH5gB6P/74Y4nryTnu8/333+f9998nLCyMq666Ks88AQEB3HXXXezYsYPRo0cX2HSvW7fOu3Vs+/btbN++Pd88OVvJzjyl0YYNGyy/fMixbds2Jk6cyPHjx/Pdl5qa6j1VXEHHopaVyMhIWrRowU8//cTmzZu9048fP86YMWOKvJyGDRty7NixPLvtu91uRo8eXWAjebbMxsXFce211/Lzzz/z0ksvFXjKtd9++817uEWDBg3o0aMHa9as4cMPP8wz37PPPluiv+GiZKu4WSlIixYtuO+++zhy5AiXX355vgHRPB4PTz/9NB988AFNmzZl9OjRBS5n7NixZGZmem/v3r2b119/neDgYIYOHQqY4xWsW7cu32OPHTtGRkZGnnrvvfdewGzwjxw5ku8x+/fvL/Jggbnl/F7feecdZs6cSePGjb17zuSoCK+/iIgUTMdwi5STZ555hvT0dF577TVatGhBnz59OPfccwkMDGTbtm3MmzePI0eO8MwzzxR5mRMmTOCKK65g3Lhxlrt5t2rVyvI80EXx97//nRMnTvDoo4/SsWNHunTpQocOHYiMjOTIkSNs2LCBRYsWERgYWKxBynr06MGoUaP4z3/+w7nnnuvdnfmzzz5j9+7d3HvvvfTo0SPPY8LDw+nUqRO//PILTqeTnj17eu/L2fp96NAhBgwYkG+Ao5YtW1K3bl0++ugjgoODiY+Px+FwMGrUqDy7k5aGO+64g5SUFP7xj3/Qvn17evToQbt27QgJCWHPnj3Mnz+fPXv25Bld+Morr6Rp06ZMmzaNXbt20a5dO9avX88PP/zgHcSpJDp16kSLFi2YPn06WVlZ3HzzzXlGVM7x5JNPsmLFCv79738zZ84cevToQa1atdizZw9r165l9erV/PLLL9SqVYtVq1Zx9dVX07lzZ1q1akVcXBx79uxh1qxZOJ3OfMdW5+xhUVBTcKbk5GRGjRrFww8/TPfu3Tn33HO9v7c5c+Zw5MgROnTokOe8yeXhoYce4o477qBLly4MGTIEj8fDN998U+igdQUZNWoU33//Pd27d+faa6+lWrVqLFy4kD179tCrV698e3t06dKFkJAQ/vWvf3Hs2DHveAePP/44AG+++SYbN27kkUce4f3336dLly5ER0eza9culi1bxqZNm9i3b5/3mOM33niDbt26ccsttzBr1izveZiXLl3KhRdeWKQv8nIrSraKm5XCvPjiiyQnJzNlyhSaN2/OpZdeStOmTUlJSeH7779n06ZNNG/enK+//rrAY/Pr1KnDyZMnSUxM5PLLL/eeh/vIkSP8+9//9o4ev2fPHtq1a0ebNm1ITEykXr16HDlyhC+//JKsrKw8zfzAgQN54oknePrpp2nWrBkDBw6kYcOGHDlyhM2bN7N48WKeeeaZYu9hdPnllxMVFcWrr75KVlYW9957b4G7odv9+ouISCHK8pxjIpLf0qVLjREjRhjNmjUzQkJCjODgYKNRo0bGDTfckO9crIWdhzu3jh07GoAxf/587zRynYfbSnHOw53bpk2bjPvvv99ITEw0IiMjjYCAAKNGjRpG9+7djfHjxxvbt2/P9xgKOD/xmaZMmWJ06tTJCA0NNUJDQ41OnToVeF7pHDnni+7QoUO++3LO3/3cc88V+Nhff/3V6NmzpxEREZHvvLPFPQ92UWzYsMG45557jFatWhnh4eFGYGCgUa9ePWPQoEHGp59+muf8z4Zhnu950KBBRkREhBEWFmb07dvXWLp06VnPw13YecVze+aZZ7zr+9133xU6X3Z2tvHWW28Z3bp1MyIjI43g4GCjQYMGxsCBA41JkyZ5zzW8a9cu4x//+IdxwQUXGLVq1TKCgoKMBg0aGFdffbXxyy+/5FtuznMXRXp6uvHZZ58Zd9xxh9GmTRujZs2ahsvlMqpXr250797dePXVV/OcW/hsv4uGDRsaDRs2LPB5Csvm2c65/cYbbxjNmzc3AgMDjQYNGhjjxo0zMjMzi3websMwjE8//dRo3769ERoaatSsWdO49tprjS1bthQ6/5w5c4xOnToZISEhBf4eU1NTjRdffNHo0KGDERYWZoSEhBiNGzc2Bg0aZLz33nt5zj9uGIaxdu1a45JLLjHCw8ONiIgI4+KLLzbWrl171vyfjVW2ipsVK3PnzjWGDBli1K1b1wgMDDSio6ONLl26GK+88oqRmppa4GNycnD06FHjjjvuMGrXrm0EBwcbbdq0MaZPn55n3mPHjhkTJkwwevToYdSpU8cICgoy6tatawwcOND45ptvCq3p8ssvN2JjY43AwEAjLi7O6NKli/H000/nOR+41XnVc/vb3/7m/b1u3Lix0Pnsfv1FRCQ/h2EUYRODiIiIiIiIiBSLjuEWERERERERKQNquEVERERERETKgBpuERERERERkTKghltERERERESkDKjhFhERERERESkDarhFREREREREykCA3QWUN4/Hw969e4mIiMDhcNhdjoiIiIiInGIYBsePH6du3bo4ndo2KJVflWu49+7dS/369e0uQ0RERERECrFr1y7i4+PtLkPEZ1Wu4Y6IiADMP+LIyEjb6jAMg+TkZKKiorSlXQqkjIgVZUSsKCNiRRkRK+WdkZSUFOrXr+/9zC5S2VW5hjvnjSIyMtLWhjs7O5v169dTt25dAgKq3MsgRaCMiBVlRKwoI2JFGRErdmVEXwCJv9CBESIiIiIiIiJlQA23iIiIiIiISBlQw20Th8NBSEiIdpeRQikjYkUZESvKiFhRRsSKMiLiG4dhGIbdRZSnlJQUoqKiSE5OtvUYbhERERERyasyf1Z3u91kZWXZXYaUg8DAQFwuV5Hm1egYNvF4PBw+fJiaNWvqHINSIGVErCgjYkUZESvKiFhRRqwZhsH+/ftJSkqyuxQpR9HR0cTFxVnu/aGG2yYej4etW7cSExOjNy8pkDIiVpQRsaKMiBVlRKwoI9Zymu1atWoRGhqq3e/9nGEYpKamcvDgQQDq1Klz1vnVcIuIiIiIiJSA2+32Nts1atSwuxwpJyEhIQAcPHiQWrVqnXX3cn1NJSIiIiIiUgI5x2yHhobaXImUt5zX3Oq4fTXcNnE4HERFRWmXEymUMiJWlBGxooyIFWVErCgjRaPfT9VT1Ndcu5TbxOVykZCQYHcZUoEpI2JFGREryohYUUbEijIi4htt4baJx+Nh9+7deDweu0uRCkoZESvKiFhRRsSKMiJWlBEpSw6Hg1mzZtldRplSw20TvXmJFWVErCgjYkUZESvKiFhRRvzToUOHuOuuu2jQoAHBwcHExcUxYMAAlixZUq517Nu3j4svvrjMlv/Pf/6Trl27EhoaSnR0dJk9z9lol3IREREREZEqZPDgwWRmZvLuu+/SpEkTDhw4wPz58zly5Ei51hEXF1emy8/MzGTIkCF06dKFyZMnl+lzFUZbuEVERERERKqIpKQkFi9ezAsvvEDv3r1p2LAhnTt3ZsyYMVxxxRXe+RwOB5MmTeLiiy8mJCSEJk2a8Omnn+ZZ1q5du7j22muJjo4mJiaGK6+8ku3bt+eZZ8qUKbRu3Zrg4GDq1KnDPffck+c5cu9SbrW8hQsX0rlzZ8LCwoiOjqZbt27s2LGj0HV98skneeCBBzjvvPNK9ssqBWq4beJ0OomNjcXp1EsgBVNGxIoyIlaUEbGijIgVZaR4DMMgNTPblothGEWqMTw8nPDwcGbNmkVGRsZZ533iiScYPHgwq1ev5sYbb2To0KGsX78eME+HNWDAACIiIli8eDFLliwhPDycgQMHkpmZCcCkSZMYOXIkd9xxB2vXrmX27Nk0a9aswOeyWl52djaDBg2iZ8+erFmzhl9++YU77rijwo8Q7zCK+sr4iZSUFKKiokhOTiYyMtLuckRERERE5JTK9lk9PT2dbdu20bhxY6pVq0ZqZjatxn1nSy1/PjWA0KCiHTH82Wefcfvtt5OWlkb79u3p2bMnQ4cOJTEx0TuPw+HgzjvvZNKkSd5pF1xwAe3bt+fNN9/kgw8+4JlnnmH9+vXepjczM5Po6GhmzZrFRRddRL169bj11lt55plnCqzD4XDwxRdfMGjQIMvldezYkRo1arBw4UJ69uxZrN/NtGnTuP/++0lKSirW487mzNe+MPqqyiYej4ctW7ZoAAoplDIiVpQRsaKMiBVlRKwoI/5p8ODB7N27l9mzZzNw4EAWLlxI+/btmTZtWp75unTpku92zhbu1atXs3nzZiIiIrxbzWNiYkhPT2fLli0cPHiQvXv30rdv3yLVZLW8mJgYhg8fzoABA7j88st5/fXX2bdvX6n8PsqSBk2zicfj4dChQzRs2FC76EiBlBGxooyIFWVErCgjYkUZKZ6QQBd/PjXAtucujmrVqtG/f3/69+/PE088wd/+9jfGjx/P8OHDi/T4EydO0KFDBz788MN895XkMASr5QFMnTqVe++9l2+//ZaPP/6Yxx9/nLlz53LBBRcU67nKkxpuEREREZGqzp0N2ekQHG53JZWaw+Eo8m7dFU2rVq3ynRP7119/5ZZbbslzu127dgC0b9+ejz/+mFq1ahW6+3+jRo2YP38+vXv3tnz+oiwPoF27drRr144xY8bQpUsXpk+froZbRERERERslHoUdi+DpB1w4gAc3wfHD8CJ/ebPk4eg3U1w5US7K5UyduTIEYYMGcKIESNITEwkIiKCZcuW8eKLL3LllVfmmXfmzJl07NiR7t278+GHH/L77797T69144038tJLL3HllVfy1FNPER8fz44dO/j888955JFHiI+PZ8KECdx5553UqlWLiy++mOPHj7NkyRJGjRqVry6r5WVlZfH2229zxRVXULduXTZu3MimTZvyfCFwpp07d3L06FF27tyJ2+1m1apVADRr1ozw8PL5ckkNt02cTifx8fHaNUcKpYyIFWVErCgjYkUZqYA8bti/Frb/BCl7oFo0hMZASPVTP2PMn9WiISgMnAXsRmwYcHQr7PwVdv4Cu36Dw39ZP/eJg/kmKSP+Jzw8nPPPP5/XXnuNLVu2kJWVRf369bn99tsZO3ZsnnmffPJJPvroI+6++27q1KnDjBkzaNWqFQChoaEsWrSIRx99lKuvvprjx49Tr149+vbt691CPWzYMNLT03nttdcYPXo0NWvW5JprrimwLqvlpaWlsWHDBt59912OHDlCnTp1GDlyJH//+98LXddx48bx7rvvem/nbJ1fsGABvXr18uXXWGQapVxERERExC4eN+xfYzbY23+CHb9ARnLRHx8YajbeQWEQFG7ePrbN3GJ9pprnmJeIOAiPg4jap39G1IHQGgU38OWosn1WL+pI1ZVR7hHEJb+ivvbawm0Tt9vNX3/9xTnnnIPLZe8bm1RMyohYUUbEijIiVpSRUpK8G1Z+aDa5WamQeQIyT56+ZKWC4TG3PGPk+gmkHYPM43mXFxwJDbpA7DmQngJpRyH12KmfR82fbvM8x2SlmpczG2xXENRtDw0uMC/1zze3jBeTMiLiGzXcNjEMg+Tk5CKfoF6qHmVErCgjYkUZESt+nZHkPXBk86lG14O3wc1peN0ZkHECMo6bDW/GidONckxjaHMDRNU7+3OcOAQ/vQpL3zndAJdEcCQ07AqNukPDbhCXCK6zfEw3DMjOONXQHzd/eus/YW61rtsWAoJLXpP3qfw4IyLlQA23iIiIiJS9lH3m8cTVIqFOWwirWbrLzzgO25fAlh9g64KiHbN8NguehWb9ocMwaD4gbwOclgS/TIRf3oSsk+a0ht3Mpjn3rt3e69XA4QKHAxxOwGFexwGBIRDboni7cjsc5jIDq0FYDd/WU6QQ+pKldKjhFhEREZHSl3nSbIC3LoAtC+DQ+rz3R9aDOm1OXdqaPyPrFH35hgH7VsOmueZz7PoNPNmn73c4oUYzcAaa1x2Qp9F1BUFwhHkarKCI09cDqpn17vgJNn1nXsLjoN2NcN618Nc38NO/ID3JfJ667aDvOGjS+9SyRUROU8NtE6fTSZMmTTTioxRKGREryohYUUbESoEZSdoFm+fCpnnm6aPqtTeP/61/PkTXL3hBhmEex7x/7ekBwHb+Cp6sXDM5IO48yEozd/VO2WNeNn59epboBtCoh7lrdeMLISo+7/NkpcP2xeZj/vrOfHxu1RtB0z5m89v4QnNk75LoMRoOb4YV78Kq6eapsxa/Yl5yxLaEPo9Dy8v8utHW+4iIbzRKuYiIiEhVlp1p7uq9ea65tfjQhsLnjaxnNt4NLjC3CO9fZzbY+9ee3uKbW3QDs/lt2hsa9zw9aFfGcfOx+1afvhzaAIY77+OrNzab71qtYMcSc8tzzi7cAIFh5rKb9TWfJ6axz7+OfLIzzQZ/xbvm80fXh15jIfFa20f09keV7bO6P49SLmenUcorOLfbzbp16zj33HM14qMUSBkRK8qIWFFGBMMwR8FO2QPH90PKXvPncfOnkbIXz+FNuLLTTj/G4YT4TtC8P0Q3gj3LzK3V+9eay/njc/NyJmcAxCZA3LlQr4O5pTmmScFbf4MjoGEX85Ij4wTs+hW2LTa3Yu9daZ7e6ti2vI+NqAstLjYvjS40j2MuSwFB0HqQeclMNQciq0KNtt5HRHyjhtsmhmGQlpamwQikUMqIWFFGxIoyUkVkZ8LBP80tzUk7zdG5U3abzXXyHsjdTJ/BAbgAI6wWjmb9oHk/c0tx7tNHJQ4xf2aehD3LYedv5vHS2elQ+1xzN/G488yBv3wZFTs4HJr1My9gng5r5y+wbREc2gjxHc0mOy7Rvl24g0LteV4b6X1ExDdquEVERETKy8kjsPt32L301OmnmkKNpubgXlH1weo4WY/bHH17zwpzC/DeFeau2e6Msz8utKY5IFlErktkHdyhtVi3O4XWva4mIDDo7MsICoPGPcxLeagWCecMMC8iIpWUGm4RERGRspCZag4Otvt32LXU/Hl0a+Hzu4JPNd9NISDEPM45I8Xc0puRcvr6mcc5A1SLNkfLjmlinjs6Mh4i65rXI+oWutu1kZ1NWvKyU6eqEhEpXw6Hgy+++IJBgwbZXUqZUcNtE5fLRcuWLXUsjBRKGREryohYUUbKiGHAiYNwZBMc3gRJO+DkITh5ONfPw3kH98qtZguo38kcQfvIVrMpP7rV3Ep98E/zcjaBYVC3rdlg121njiJevXGJdrNWRsSKMuKfDh06xLhx45gzZw4HDhygevXqtGnThnHjxtGtW7dyq2Pfvn1Ur17CswlY2L59O08//TQ//PAD+/fvp27dutx000089thjBAVZ7NFTitRw28ThcBAdHW13GVKBKSNiRRkRK8pIMRz6C9Z/aZ6yyuHEPF+z8/TFkwVHt5lN9pEt5tbmogiKgPgOEN8Z6nc2BxPLfXx0Dnc2JO8yl31ks/l8wZHm4GLVIiE46vT18NqlNmiXMiJWlBH/NHjwYDIzM3n33Xdp0qQJBw4cYP78+Rw5cqRc64iLiyuzZW/YsAGPx8Nbb71Fs2bNWLduHbfffjsnT57k5ZdfLrPnPZNOC2aT7OxsVq5cSbt27QgI0Pcekp8yIlaUEbFS6TOSnQknD0JYbNEG43Jnm6NvJ+82m9PYBHCdZb2z0mH9/8HyabDjp+LV5nCap7yq0dw8FVV4bQiradYaFguhNcyfwREV+hzNlT4jUubKOyMV5bN6UVXG04IlJSVRvXp1Fi5cSM+ePQudz+Fw8OabbzJ79mwWLlxInTp1ePHFF7nmmmu88+zatYuHHnqI77//HqfTyYUXXsjrr79Oo0aNvPNMmTKFV155hc2bNxMTE8PgwYOZOHGi9zly71JutbyFCxfyyCOP8McffxAYGEjr1q2ZPn06DRs2LNK6v/TSS0yaNImtW89yeE8R6bRglYDbXcAxWCK5KCNiRRkRK5UuI4ZhDii2egas+/z0uZ1Da5wa7Cvu1KWuOT15lzkyd9Iu85RVuY9vDgwzd7eO72ie5qpeR4iobW7NXvEurJoOaUfNeR1OaH6RuWu24cl/cTgguiHUbG4OcBbTxLcRuSuQSpcRKXfKSDEYBmSl2vPcgaFF+oIvPDyc8PBwZs2axQUXXEBwcOHvZU888QTPP/88r7/+Ou+//z5Dhw5l7dq1JCQkkJWVxYABA+jSpQuLFy8mICCAZ555hoEDB7JmzRqCgoKYNGkSDz74IM8//zwXX3wxycnJLFmypMDnslqe0+lk0KBB3H777cyYMYPMzEx+//13HMX4UjM5OZmYmAL2MipDarhFRETEfse2w5pPzEY7z8BiDsCA1CPm5cC6sy/HGWgOFJZ61Nzte/upczrnCK8NJw6cvh0ZD+1vgXY3mY8TEfFFVio8W9ee5x671zybgIWAgACmTZvG7bffzn//+1/at29Pz549GTp0KImJiXnmHTJkCH/7298AePrpp5k7dy7/+c9/ePPNN/n444/xeDy888473qZ36tSpREdHs3DhQi666CKeeeYZHnroIe677z7vMjt16lRgXVbL69ixI8nJyVx22WU0bdoUgISEhCL/ejZv3sx//vOfct2dHNRwi4iISFnJTDUb6aNbzC3Q2WnmbuLZ6eDOhOwMc6CwI1vM8y3nCAyDVldA4nXQ6EKzcT6+D1L2mT+P7zd3HTcMiK5vbnmOqm9eD48zT63l8cDhjebW8t1LYfcyOLjebLYdTjhnIHQYbp7zuZSOhxYRqSwGDx7MpZdeyuLFi/n111/55ptvePHFF3nnnXcYPny4d74uXbrkeVyXLl1YtWoVAKtXr2bz5s1ERETkmSc9PZ0tW7Zw8OBB9u7dS9++fYtUk9XyLrroIoYPH86AAQPo378//fr149prr6VOnTqWy96zZw8DBw5kyJAh3H777UWqp7So4baJy+UiMTFRIz5KoZQRsaKMiJViZ8Qw4MAfsOEr83zPddtC3fbm+ZvP9pgTB2D/WvOxRzabg4sd3Wo2xUXmgCa9oM1QaHkZBIefvis0xrzUbl30xTmdUCvBvLS/xZyWnmLWWL2hecos0fuIWFJGiikw1NzSbNdzF0O1atXo378//fv354knnuBvf/sb48ePz9Nwn82JEyfo0KEDH374Yb77YmNjcTqLd7pBq+WBucX73nvv5dtvv+Xjjz/m8ccfZ+7cuVxwwQWFLnfv3r307t2brl278vbbbxerptKghttG5TkcvVROyohYUUbESpEycnA9/PGFeTn8V/77w+NOn4KqThtIT4YDa80me/86SD1c+LKDo6BGE3MrdHC4ea7pgGoQEHTqejBUi4IWF5d9E1wtEhp2sZ6vitH7iFhRRorB4SjSbt0VUatWrZg1a1aeab/++iu33HJLntvt2rUDoH379nz88cfUqlWr0AHuGjVqxPz58+ndu7fl8xdleQDt2rWjXbt2jBkzhi5dujB9+vRCG+49e/bQu3dvOnTowNSpU4v9JUBpUMNtE7fbzbJly+jYsaNGBZUCKSNiRRkRKwVmxOM2d8lO2gnbFsEfn8OhDacf5AqG5v2hWjTsXQmH1sOJ/fDXN+alIA6nOVp33LlQ8xxzQLGcS0j1Cj1Kd1Wn9xGxooz4nyNHjjBkyBBGjBhBYmIiERERLFu2jBdffJErr7wyz7wzZ86kY8eOdO/enQ8//JDff/+dyZMnA3DjjTfy0ksvceWVV/LUU08RHx/Pjh07+Pzzz3nkkUeIj49nwoQJ3HnnndSqVYuLL76Y48ePs2TJEkaNGpWvLqvlZWVl8fbbb3PFFVdQt25dNm7cyKZNm/J8IZDbnj176NWrFw0bNuTll1/m0KFD3vvK8nRkZ9JfjYiIiD/KzsSxZQHxm77EuWeKOYJ38k5I2Que7LzzuoKgaV8492rz2OZqubYsZKaaW7L3rjQv+9eYW6Rrnwtx55lNdq1WEBhSvusnIiIlEh4ezvnnn89rr73Gli1byMrKon79+tx+++2MHTs2z7xPPvkkH330EXfffTd16tRhxowZtGrVCoDQ0FAWLVrEo48+ytVXX83x48epV68effv29W6hHjZsGOnp6bz22muMHj2amjVr5jmtWG5Wy0tLS2PDhg28++67HDlyhDp16jBy5Ej+/ve/F7i8uXPnsnnzZjZv3kx8fHye+8rzzNg6D7dNsrOz9W2hnJUyIlaUEcknOxO2/Qh/zDKPw845pdaZnAHm7tu1WkGrK6HFJRASXY6FSkWh9xGxUt4ZqSif1YuqMp6Hu6jOPEe25KXzcIuIiFQ2WenmOZ8Dgs8+crbHY556JvMkZJ4wR/n+88t8TbYRVotD0W2pcU5nXNVzjeQdUUcjc4uIiJQDNdw2cblcdOzYUSM+SqGUEbGijPiBk4dhx8+wY4l52b8OOLXjmTPAHFzMFXTqZ4DZkGeehKyThS8zrJZ5Sq1Wg6BBF2IMcLpcOo5aCqT3EbGijIj4Rg23jTIzMwkJ0TFvUjhlRKwoIxXQiYOw8Wvzp9NlNs7OwFM/TzW++9eZDXbuwcrO5Mk2t16fjcMJQeHmwGTN+5tNdsOup7deGwaZaWnKiJyV3kfEijJSNVWxI4/LjBpum7jdbtasWaNjpqRQyohYUUYqkJR9sP7/zN26dyzBu5W6KGIToFE3s1Fu0AWCI8xjsbPTwZ1x+ron2xyYLCjMbLKDwswt32fZcq2MiBVlRKwoIyK+0V+NiIhISSTvMRvsP7+EXb/mva9eB3MEb0+2eRouTza4s07frt7IbLIbdIWwGvmXHVwuayAiIiJlTA23iIhIUaUlmQ322pmw/SfybMmuf7454nfC5RDdwK4KRUREpAJRw20jDT4hVpQRsaKMFFF6ChzfZ/7MSIGM43kv7gwIj4PIOubpsiLqQngt81jo7Az46ztY+4n50515erkNupjHTSdcDlH1bFu9s1FGxIoyIlaUEZGS03m4RUSkcslKM0f3zko9dWqs1NPXM06YjXXybkjZY/5M3gMZycV/HofLPH1WxvG8j49NgMRr4bxrtCVbRKSUVbbP6v58Hm45O52Hu4IzDIPk5GSioqJw6FQtUgBlRKxUuYxkpsKS1+Hnf5vNdXEFR0FIFARHmgOTeS+R5gjiJ/abg5+l7DWvG25I2W0+NqKu2WAnXgu1z600p9iqchmRYlNGxIoyIuIbNdw2cbvdbNiwQSM+SqGUEbFSZTJiGLD2U5g33txqDea5qQNDzZG6A0NPj94dGAoRcRBZD6Lizd28I0/9DI4o+nO6s+HkQbP5Bqjb7vSptiqRKpMRKTFlRKwoIyK+0V+NiIj4zp0FrsDSX+6e5fDNP2D37+btqAZw0dPm4GRluaXFFWAeyx1Zt+yeQ0RExCbDhw/n3XffBSAgIID4+HiGDBnCU089VSV3jf/888/573//y/Llyzl69CgrV66kbdu2pbJsNdwiIlJymakwexT88bm5JTm+A9TrCPEdoU4bc8tzSSTvhh/+Caunm7cDw+DCB6DLPSVfpoiIiHgNHDiQqVOnkpWVxfLlyxk2bBgOh4MXXnjB7tIKlZmZSVBQUKkv9+TJk3Tv3p1rr72W22+/vVSX7SzVpUmRORwOQkJCdCyMFEoZESu2ZyRlH0y9GNZ9CoYHknfCH1/A94/BlAHwXDy81cNsyBe9DKs/hu1LIGmnucs2mLuLp+yDjd/Cwhdgxg3wamt4rfXpZrvN9TBqOfR4WM12MdmeEanwlBGxooz4r+DgYOLi4qhfvz6DBg2iX79+zJ0713u/x+Phueeeo3HjxoSEhNCmTRs+/fRT7/3Hjh3jxhtvJDY2lpCQEJo3b87UqVO9969du5Y+ffoQEhJCjRo1uOOOOzhx4oT3/l69enH//ffnqWnQoEEMHz7ce7tRo0Y8/fTT3HLLLURGRnLHHXcAsGTJEnr16kVoaCjVq1dnwIABHDt2rEh1F+Tmm29m3Lhx9OvXr9i/Ryvawm0Tl8tFmzZt7C5DKjBlRKzYmpG9q2DG9XB8L4TEwOB3zIHH9iyD3cvNnycOwL7V5uVMDqc5EJk70zxWuiANusJFz5hbzaVE9D4iVpQRsaKMlEzqWQb3dDldBLuCizSv0+GkWkA1y3lDA0NLUOVp69at4+eff6Zhw4beac899xwffPAB//3vf2nevDmLFi3ipptuIjY2lp49e/LEE0/w559/8s0331CzZk02b95MWloaYG4xHjBgAF26dGHp0qUcPHiQv/3tb9xzzz1MmzatWLW9/PLLjBs3jvHjxwOwatUq+vbty4gRI3j99dcJCAhgwYIFuN3uItVd3tRw28Tj8XD48GFq1qyJ06kdDSQ/ZUSs2JaR9V/B57ebI4XXbAE3fAQxTcz7mpz6R2YY5m7he5bBwfXm9aSdp0/X5c48PQK4w2kup06b05e486BaxT8dTEWn9xGxooyIFWWkZM6ffn6h911Y70Le7Pem93avT3qRlp1W4Lwda3dk6sDTW40HfjaQYxnH8s23dtjaYtf41VdfER4eTnZ2NhkZGTidTiZOnAhARkYGzz77LPPmzaNLly4ANGnShJ9++om33nqLnj17snPnTtq1a0fHjh0Bc2t0junTp5Oens57771HWFgYABMnTuTyyy/nhRdeoHbt2kWus0+fPjz00EPe2zfccAMdO3bkzTdP/w5bt25d5LrLmxpum3g8HrZu3UpMTIzevKRAyohYKfeMGAYs+RfMexIwoGkfGDINqkXln9fhgOj65qX1VWcWbm7VTtplzlerFQT59s28FEzvI2JFGREryoj/6t27N5MmTeLkyZO89tprBAQEMHjwYAA2b95Mamoq/fv3z/OYzMxM2rVrB8Bdd93F4MGDWbFiBRdddBGDBg2ia9euAKxfv542bdp4m22Abt264fF42LhxY7Ea7pyGPseqVasYMmRIgfMWpe7ypoZbRETMrc/rv4ItP5hNcLVos5EOOfWzWjTsWAKrZ5jzd7odBj5vjuZdXE6neequiLhSXAEREZGK47cbfiv0PtcZp5lceO3CQud1OvJ+yfHt4G99qiu3sLAwmjVrBsCUKVNo06YNkydP5rbbbvMeaz1nzhzq1auX53HBwebu8BdffDE7duzg66+/Zu7cufTt25eRI0fy8ssvF+n5nU4nhmHkmZaVlVVgnbmFhBQ+nktR6i5varhFRPxBWhIc/gsObYTDG+HwJvOY6lqtoFYC1G4NMU3zNsiHNsL62bD+/wo+zrogDicMfAHOv6NMVkNERMQfFOeY6rKatzicTidjx47lwQcf5IYbbqBVq1YEBwezc+fOs+6GHRsby7Bhwxg2bBgXXnghDz/8MC+//DIJCQlMmzaNkydPehvmJUuW4HQ6adGihfex+/bt8y7L7Xazbt06evfufdZaExMTmT9/Pk8++WS++4pad3lSw20Th8NBVFSURnyUQikjclaZJ3H++hbnrZuD66dd5gBlBdnw1enrrmCIPQdqngP715oNeg6H0xykrOUlEBwB6clmE5+edPq64YEuI6FZ3zJcMSlNeh8RK8qIWFFGqo4hQ4bw8MMP88YbbzB69GhGjx7NAw88gMfjoXv37iQnJ7NkyRIiIyMZNmwY48aNo0OHDrRu3ZqMjAy++uorEhISALjxxhsZP348w4YNY8KECRw6dIhRo0Zx8803e3cn79OnDw8++CBz5syhadOmvPrqqyQlJVnWOWbMGM477zzuvvtu7rzzToKCgliwYAFDhgyhZs2alnUX5OjRo+zcuZO9e/cCsHHjRgDi4uKIi/Ntjzw13DZxuVzeQIoURBmRArmzYdWHsOBZnCf2k2cnq4i6ENvCvNQ8xxyY7MAfcPBPOLgBsk6ajfb+UwOrOAOhSS9IuBxaXALhsTaskJQlvY+IFWVErCgjVUdAQAD33HMPL774InfddRdPP/00sbGxPPfcc2zdupXo6Gjat2/P2LFjAQgKCmLMmDFs376dkJAQLrzwQj766CMAQkND+e6777jvvvvo1KkToaGhDB48mFdffdX7fCNGjGD16tXccsstBAQE8MADD1hu3QY455xz+P777xk7diydO3cmJCSE888/n+uvvx7Asu6CzJ49m1tvvdV7e+jQoQCMHz+eCRMmFPt3mZvDOHPHeT+XkpJCVFQUycnJREbaNwKux+Nh79691K1bVwNQSIGUEcnDMGDT9zB3PBxab06KbkBSq2FEJfTEGdvi7KN6ezyQtMNsvg9thKj6cM5FBQ94Jn5D7yNiRRkRK+WdkYryWb2o0tPT2bZtG40bN6ZatWrWDxC/UdTXXlu4beLxeNi9ezdxcXH6BycFUkaqmKw0wAGuQDhjMBX2LDcb7e2LzdvVoqHnI7jbDWfjqrV0rNMOZ4DF27nTCTGNzUvLS8tiDaQC0vuIWFFGxIoyIuKbCttwP//884wZM4b77ruPf/3rX4D5LcJDDz3ERx99REZGBgMGDODNN98s1rDyIiIVRuZJWPc5LJ9mnq86h8Np7u7tOnVJO3W+TVcwXHAndH8AQqpDdrYtZYuIiIhI0VTIhnvp0qW89dZbJCYm5pn+wAMPMGfOHGbOnElUVBT33HMPV199NUuWLLGpUhGREti3xmyy13wCmcfz3294wJ1hXgBwQJuh0Psx87zWIiIiIlIpVLiG+8SJE9x4443873//45lnnvFOT05OZvLkyUyfPp0+ffoAMHXqVBISEvj111+54IIL7Cq5RJxOJ7Gxsdo1RwqljPgRjxuObDHPY73iPdi74vR9MU2g/TBIvA6CwsCTbQ525s4yf3qyzeOsCzhntTIiVpQRsaKMiBVlRMQ3Fa7hHjlyJJdeein9+vXL03AvX76crKws+vXr553WsmVLGjRowC+//FJow52RkUFGRob3dkpKCgDZ2dlkn9od0+l04nQ68Xg8eDwe77w5091ud56Tshc23eVy4XA4vMvNPR3Mc8vl1qRJE28tuQUEBGAYRp75HQ4HLpcrX42FTbdrnQqbrnUq2To1atQoz33+sE6V/nU6vBnnmo8gKxVCY3CG1cQdHA0h1TFCa0BQBM7knTgP/Yln/1ocB/6AQxtwZKd7l2c4AzFaXo7R/haMht1wBQRar1Ou30POOnk8Hho2bIjH48EwDL1OWqcCpzds2NB7Oh9/WaecGv3pdbJznZo2bYrH48nzvJV9nfzxdbJrnRwOh/d/Tc5zlOU65b4u4g8qVMP90UcfsWLFCpYuXZrvvv379xMUFER0dHSe6bVr12b//v2FLvO5554r8KToK1eu9J6EPTY2lqZNm7Jt2zYOHTrknSc+Pp74+Hj++usvkpOTvdObNGlCrVq1WLduHWlpad7pLVu2JDo6mpUrV+Z5o0lMTCQoKIhly3IdownUqFGDOnXqsG7dOu80l8tFp06dSE5OZsOGDd7pISEhtGnThsOHD7N161bv9KioKBISEti7dy+7d+/2TrdrnTp27EhmZiZr1qzROvm4TnXr1uW3337L841yZV+nSvs61axJU8dO0n54hdA9S3CQ9+QOZwxx5pV7W4DbFYwR24qAxMGs8rQkIyACDgOHV/i0TqmpqYSGhup10joVuk6pqam0bduWmJgYv1kn8L/Xya51ioyMJDg4mMDAQO/5Zyv7Ovnj62TnOgUEBPDTTz8RGhpaLuukkb7F31SY04Lt2rWLjh07MnfuXO+x27169aJt27b861//Yvr06dx66615tlYDdO7cmd69e/PCCy8UuNyCtnDXr1+fI0eOeE81YMc3nm63mxUrVtCxY0fvloccleUbzzPX6WzTtU7FXyePx8PSpUtp3769t4bKvk6V7nXKPIFj7Sc4l76D4/BG73yepv2gVgKOtGM40o5hpB6B1CPm4GbpSRBZD0ftc/HUaoVRqxVG7dYQ3QhXYFCprlNmZiYrVqygffv2BAQEVN3XSetU6PTc/2sCAwP9Yp1y1+gvr5Od6+R2u1m5ciXt27fP8wVvZV4nf3yd7Fyn7Oxsli1blufzSFmu04kTJ6hevbpOCyYVXqU7Ldjy5cs5ePAg7du3905zu90sWrSIiRMn8t1335GZmUlSUlKerdwHDhwgLi7/sY05goODCQ4Ozjc9ICCAgDNOo5PzJnSmnDeVok4/c7mFTXc4HDgcjgLnL2x6YTUWd3pZrdPZpmudijfd4/F4/xmd+RyVdZ2gErxOhoFz30qcaz6BVTMg49SWgaAIaHcjdLodZ81meZdfYCV5t3CXtPbCpudkI+fnmV/K5KvF316nEkyvquuU87+muLUXNr0irJNVjcWdrnUypxdn/sqyTv74OpX3OuW8hxT0eaQs1qmgeUQqswrTcPft25e1a9fmmXbrrbfSsmVLHn30UerXr09gYCDz589n8ODBAGzcuJGdO3fSpUsXO0oWEX9y6C9YOxPWfQpHT+/yRkwT6Px3aHsDVKv437SLiIiISMVRYRruiIgIzj333DzTwsLCqFGjhnf6bbfdxoMPPkhMTAyRkZGMGjWKLl26VLoRysH89i4+Pl7f4kmhlJESMgzIToeM46cuKZCZmnPn6Xlybu9dCWs/hf2nj0MjIARaXAxtrodm/aCCvgbKiFhRRsSKMiJWlBER31SYhrsoXnvtNZxOJ4MHDyYjI4MBAwbw5ptv2l1WieS8eYkURhkpovRkWPI6rPvcPH4647h5Kq3icgZA075w3hCz2Q4OL/VSS5syIlaUEbGijIgVZcQ/DR8+nHfffRcwDw2Ij49nyJAhPPXUU1XuWPSsrCwef/xxvv76a7Zu3UpUVBT9+vXj+eefp27duj4vv0I33AsXLsxzu1q1arzxxhu88cYb9hRUitxuN3/99RfnnHNOocfRSNWmjFjIzoTlU+HHF8wBy/JxQHCEeQmoBg4neAcoPPXT4YDwWtD6Kki4EsJqlFf1pUIZESvKiFhRRsSKMuK/Bg4cyNSpU8nKymL58uUMGzYMh8NR6GDUFUFmZiZBQUGluszU1FRWrFjBE088QZs2bTh27Bj33XcfV1xxRb6zDpSE9g2xiWEYJCcn5xkhUiQ3ZaQQhgF/zII3z4dvHjGb7ZrnwJBpMPJ3eHA9jNkN447CmF3w4J9w7woYtQzuWXrq8rt5GfkbDPs/6Dii0jXboIyINWVErCgjYkUZ8V/BwcHExcVRv359Bg0aRL9+/Zg7d673fo/Hw3PPPUfjxo29p3f79NNPvfcfO3aMG2+8kdjYWEJCQmjevDlTp0713r927Vr69OlDSEgINWrU4I477uDEiRPe+3v16sX999+fp6ZBgwYxfPhw7+1GjRrx9NNPc8sttxAZGckdd9wBwJIlS+jVqxehoaFUr16dAQMGcOzYsSLVfaaoqCjmzp3LtddeS4sWLbjggguYOHEiy5cvZ+fOnSX63eZWobdwi4h4GQbs/BXmPgG7l5rTwmpB7zHQ7hZw6e1MREREKgZPamrhd7pcOHOdRems8zqdOHPt4l3YvM5c50kviXXr1vHzzz/TsGFD77TnnnuODz74gP/+9780b96cRYsWcdNNNxEbG0vPnj154okn+PPPP/nmm2+oWbMmmzdv9p7f/eTJkwwYMIAuXbqwdOlSDh48yN/+9jfuuecepk2bVqzaXn75ZcaNG8f48eMBWLVqFX379mXEiBG8/vrrBAQEsGDBAu/p6KzqLork5GQcDkees2OVlD6hikjFlLIP9q0yBzXLuZw8ZN4XGApd74Wu95i7jIuIiIhUIBvbdyj0vrCePWjw1lve2391645xqlE9U2inTjR8/z3v7c19++E+tSU3t4QN64td41dffUV4eDjZ2dlkZGTgdDqZOHEiABkZGTz77LPMmzfPe0aoJk2a8NNPP/HWW2/Rs2dPdu7cSbt27ejYsSNgbo3OMX36dNLT03nvvfcICwsDYOLEiVx++eW88MIL1K5du8h19unTh4ceesh7+4YbbqBjx455xvJq3bp1keu2kp6ezqOPPsr1119fKueCV8NtE6fTSZMmTTTioxTKrzPizoYTB+DEfjh+xs+UvbB/LRzfl/9xzkDz9Fy9x0JEXPnXXcH4dUakVCgjYkUZESvKiP/q3bs3kyZN4uTJk7z22msEBAR4T7+8efNmUlNT6d+/f57HZGZm0q5dOwDuuusuBg8ezIoVK7jooosYNGgQXbt2BWD9+vW0adPG22wDdOvWDY/Hw8aNG4vVcOc09DlWrVrFkCFDCpy3KHWfTVZWFtdeey2GYTBp0qQi13g2arht4nQ6qVWrlt1lSAXmdxnJOA6bvoc/Z8OmuZB18uzzO5wQ2xLqtjMvddpC3LkQGFIu5VYGfpcRKXXKiFhRRsSKMlIyLVYsL/zOMwafO2fJT4XPe8YXHc3mz/OlrDzCwsJo1qwZAFOmTKFNmzZMnjyZ2267zXus9Zw5c6hXr16exwWf2h3+4osvZseOHXz99dfMnTuXvn37MnLkSF5++eUiPb/T6cw3NkBWVlaBdeYWElL4Z8Gi1F2YnGZ7x44d/PDDD6WydRvUcNvG7Xazbt06zj33XI34KAXyi4ykHYON38L62bB5PrgzTt/ncEF4bYioDRF1Tl2PM3/WSoC48yAorPBli39kRMqUMiJWlBGxooyUTHGOqS6reYvD6XQyduxYHnzwQW644QZatWpFcHAwO3fuPOtu2LGxsQwbNoxhw4Zx4YUX8vDDD/Pyyy+TkJDAtGnTOHnypLdhXrJkCU6nkxYtWngfu2/f6T0ac7LWu3fvs9aamJjI/PnzefLJJ/PdV9S6z5TTbG/atIkFCxZQo0bpDaarhtsmhmGQlpamER+lUJUmI+nJkLwbknZB8i7zevJuSNoJe1fkPSd2TBNIuAJaXQF12uX71laKp9JkRGyjjIgVZUSsKCNVx5AhQ3j44Yd54403GD16NKNHj+aBBx7A4/HQvXt3kpOTWbJkCZGRkQwbNoxx48bRoUMHWrduTUZGBl999RUJCQkA3HjjjYwfP55hw4YxYcIEDh06xKhRo7j55pu9u5P36dOHBx98kDlz5tC0aVNeffVVkpKSLOscM2YM5513HnfffTd33nknQUFBLFiwgCFDhlCzZk3Lus+UlZXFNddcw4oVK/jqq69wu93s378fgJiYGJ9PQ6aGW0RKZtfv5mm59q48+3yxCWaDnXAF1G6d61zYIiIiIlJRBAQEcM899/Diiy9y11138fTTTxMbG8tzzz3H1q1biY6Opn379owdOxaAoKAgxowZw/bt2wkJCeHCCy/ko48+AiA0NJTvvvuO++67j06dOhEaGsrgwYN59dVXvc83YsQIVq9ezS233EJAQAAPPPCA5dZtgHPOOYfvv/+esWPH0rlzZ0JCQjj//PO5/vrrASzrPtOePXuYPXs2AG3bts1z34IFC+jVq1dxf5V5OIwq9nVVSkoKUVFRJCcnl9p++SWRnZ3NsmXL6NixIwEB+t5D8quwGTlxCOaNh1Ufnp4WEgNR8RBVH6Lrn7oeD7XPg5rN7KvVz1XYjEiFoYyIFWVErJR3RirKZ/WiSk9PZ9u2bTRu3JhquU7fJf6vqK+93llt4nK5aNmypY6FkUJVuIy4s2HZFPjhGchINqe1uwn6jDOPw5ZyV+EyIhWOMiJWlBGxooyI+EYNt01K60Tq4r8qVEZ2/gpzRsOBtebtuES49BWo39neuqq4CpURqZCUEbGijIgVZUTEN2q4bZKdnc3KlStp166dduGSApV5RtxZsHupOXr4lh9g/xpzujPAHEHcGWAOauZwQeph875qUdB3HHS4FZz6pttueh8RK8qIWFFGxIoyIuIb/dXYyO12212CVHClmhHDgGPbzOZ68w+wbRFkHs8/X+5RxXNrdzP0mwBhNUuvJvGZ3kfEijIiVpQRsaKMiJScGm4Rf+Vxw8E/zd3Bd/xs/jy+N+88oTWgSW9o1hcadgVXsNlwe7LB8Jy+Xi0aourZshoiIiIiFV0VG4daKPprroZbxJ+cOGSOHr79J/O0XTmDm+VwBkD986FpH7PJjmujc2GLiIiIlFBgYCAAqamphISE2FyNlKfU1FTgdAYKo9OC2cQwDNLS0ggJCcGh8xJLAYqVkSNb4JeJsGo6ZKefnh4Ubg5s1qCLeanXAYJCy7ZwKTd6HxEryohYUUbESnlnpKJ8Vi+Offv2kZSURK1atQgNDdXfkp8zDIPU1FQOHjxIdHQ0derUOev82sJto6CgILtLkArOMiN7VsCS12H9bHMXcDCb6vOGmA127XPBpT9zf6b3EbGijIgVZUSsKCNnFxcXB8DBgwdtrkTKU3R0tPe1Pxt9EreJ2+1m2bJldOzYUSM+SoHOmpEtC2DxK7B98elpzS+Cbvebx2Lrm9UqQe8jYkUZESvKiFhRRqw5HA7q1KlDrVq1yMrKsrscKQeBgYFFPje9/mpEKpOTR+Cbh2HdZ+ZtZ4C5NbvrKKjd2t7aRERERKowl8tV5CZMqg413CKVxZ+zYc6DcPKQeW7sTn8zG+3o+nZXJiIiIiIiBVDDLVLRnTwM34+BPz43b8cmwKA3oV57e+sSEREREZGz0ijlNjEMA7fbjcvl0kiGUiDDMPCs+wLnNw/jSD1sbtXu/gD0fAQCgu0uTyoAvY+IFWVErCgjYqW8M1JRPquLlBZt4bZRZmamztcnBdu7Cha9iGvDHPN2rVbmVu267WwtSyoevY+IFWVErCgjYkUZESk5p90FVFVut5s1a9bgdrvtLkUqCsOAzfPg3cvh7Z44NszBcDjxdH8I7lioZlvy0fuIWFFGxIoyIlaUERHfaAu3iN3cWeao4z//Bw6sM6c5XHhaX8XaqP607nUNTp2GQ0RERESk0tGneJGyZBiw9B3480twOMEVCM5AcAWAK8i8vv0nSNltzh8YBh2GwwV34QmvQ9qyZbaWLyIiIiIiJaeG20Y6T5+fyzgOs+6G9bOt5w2rBRfcCR1HQEh1c1p2tjIilpQRsaKMiBVlRKwoIyIlp1HKRcrCoY3w8U1w+C9zK3avf0B0Q/BkmbuQuzPBk23+jKgDCVdAYDW7qxYRERGxlT6ri7/RFm6bGIZBcnIyUVFROg2Hv1n3OXx5D2SdhIi6cO17UL9TsRejjIgVZUSsKCNiRRkRK8qIiG80SrlN3G43GzZs0IiP/sSdBd89Bp/eajbbjS6Evy8qUbMNyohYU0bEijIiVpQRsaKMiPhGW7hFSkPKPvjsNtixxLzd7X7o84Q5OJqIiIiIiFRJ6gZEfOHxwPKpMG8CZKRAUARcNQkSLre7MhERERERsZkabps4HA5CQkJ0LExlduBP+Op+2PWbebtue7j6bajZvFQWr4yIFWVErCgjYkUZESvKiIhvNEq5SHFlpcGil2DJ6+ZI40Hh0HccdPobOHXaDBEREZGS0md18TcaNM0mHo+HgwcP4vF47C5FimPrQpjUFRa/YjbbLS+Dkb/D+X8v9WZbGREryohYUUbEijIiVpQREd+o4baJx+Nh69atevOqLI5sgY9vhveuhKNbzXNnX/cBDP0QouqVyVMqI2JFGREryohYUUbEijIi4hsdwy1yNicOwY8vmAOjebLB4YSOt5m7kFfTbk4iIiIiIlI4NdwiBck8Cb+8CUv+BZknzGnNB0C/CVC7lZ2ViYiIiIhIJaGG2yYOh4OoqCiN+FjRGAas/AAW/BOO7zOn1WkLFz0NjXuUaynKiFhRRsSKMiJWlBGxooyI+EajlIvkMAz49h/w23/N29ENzV3HW18NTg13ICIiIlLW9Fld/I26CJt4PB52796tASgqkh+ePt1s9x0P9yyF866xrdlWRsSKMiJWlBGxooyIFWVExDdquG2iN68KZvEr5gXgkpfhwgchINjWkpQRsaKMiBVlRKwoI2JFGRHxjRpukd/egvlPmdf7PwWdb7e3HhERERER8QtquKVqW/kBfPOIeb3HI9DtPnvrERERERERv6GG2yZOp5PY2FicGozLPus+h9mjzOsX3A29x9pbzxmUEbGijIgVZUSsKCNiRRkR8Y1GKZeqaeO38PGN4MmG9rfA5f8Gne5CRERExFb6rC7+Rl9V2cTj8bBlyxYNQFHeknfDVw+cbrbPGwKX/atCNtvKiFhRRsSKMiJWlBGxooyI+EYNt008Hg+HDh3Sm1d5SdkLc0bDv9vBsilms33uYBg0CZwuu6srkDIiVpQRsaKMiBVlRKwoIyK+CbC7AJEydXw//PQaLJsK7gxzWqMLodcYaNTN3tpERERERMSvqeEW/+TOggX/hF8nQXa6Oa1BV+g9Bhr3sLc2ERERERGpEtRw28TpdBIfH68RH8tCWhLMHA5bF5i3659vbtFu0qtCHqtdGGVErCgjYkUZESvKiFhRRkR8o1HKxb8c3QbTr4PDGyEwFAa9Ca0GVapGW0RERKSq0md18Tf6qsombreb9evX43a77S7Ff+z8Fd7pazbbEXVhxLfQ+qpK22wrI2JFGREryohYUUbEijIi4hvtUm4TwzBITk6miu1gUHbWfAJfjgR3JtRpC9d/BJF17K7KJ8qIWFFGxIoyIlaUEbGijIj4Rg23VG6GAQufhx+fN2+3vAyufhuCwuytS0REREREqjw13FJ5eTzmVu3V083b3e6HvuNBg3qIiIiIiEgFoIbbJk6nkyZNmmjER1/Mf9Jstp0BcNm/oP3NdldUqpQRsaKMiBVlRKwoI2JFGRHxjUYpl8pp2VT46n7z+lVvQZuhtpYjIiIiIr7TZ3XxN/qqyiZut5vVq1drxMeS2DQP5jxkXu81xm+bbWVErCgjYkUZESvKiFhRRkR8o4bbJoZhkJaWphEfi2v/Wpg5HAw3tLkeej5qd0VlRhkRK8qIWFFGxIoyIlaUERHfqOGWyiNlL3x4LWQeh0YXwuX/rrTn2BYREREREf+nhlsqh4zjMP1aOL4Xap4D170PAUF2VyUiIiIiIlIoNdw2cblctGzZEpfLZXcpFZ87G2beau5OHhYLN86EkOp2V1XmlBGxooyIFWVErCgjYkUZEfGNTgtmE4fDQXR0tN1lVA7f/gM2z4WAELj+Y6jeyO6KyoUyIlaUEbGijIgVZUSsKCMivtEWbptkZ2ezdOlSsrOz7S6lYlv+Liz9H+CAq9+G+A52V1RulBGxooyIFWVErCgjYkUZEfGNGm4b6fQKFnb9Dl+PNq/3fgxaXWFvPTZQRsSKMiJWlBGxooyIFWVEpOTUcEvFlLIPPr4Z3JmQcDlc+JDdFYmIiIiIiBSLGm6peLIz4JOb4cR+iE2AQZPAqaiKiIiIiEjl4jCq2FnsU1JSiIqKIjk5mcjISNvqMAyDtLQ0QkJCcOhc0qcZBsweBSvfh2pRcPsCqNHU7qpsoYyIFWVErCgjYkUZESvlnZGK8lldpLRos6GNgoJ0Hul8lk02m22HE66ZUmWb7RzKiFhRRsSKMiJWlBGxooyIlJwabpu43W6WLVumQShy2/EzfPOoeb3veGjWz956bKaMiBVlRKwoI2JFGREryoiIb9RwS8WQvBs+uQU82dD6auh2n90ViYiIiIiI+EQNt9jPMGDWXXDyENQ+D66cCDqOTEREREREKjk13GK/tTNh2yIIqAbXvgtBYXZXJCIiIiIi4jONUm4TwzBwu924XK6qPSpo2jGY2Mncut3ncejxsN0VVRjKiFhRRsSKMiJWlBGxUt4ZqSif1UVKi7Zw2ygzM9PuEuw3/2mz2a7RHLrea3c1FY4yIlaUEbGijIgVZUSsKCMiJaeG2yZut5s1a9ZU7REfdy+HZVPM65e9CgHB9tZTwSgjYkUZESvKiFhRRsSKMiLiGzXcYg93Nnx1P2BA4nXQuIfdFYmIiIiIiJQqNdxij6XvwP41UC0KLnrG7mpERERERERKnRpuG7lcLrtLsEfKPvjhVJPddzyE17K3ngqsymZEikwZESvKiFhRRsSKMiJSchqlXMrfzOHwxxdQryPcNhec+t5HRERERPRZXfyPOh2bGIZBUlISVez7Dtg8z2y2HU5zoDQ124WqshmRIlNGxIoyIlaUEbGijIj4Rt2OTdxuNxs2bKhaIz5mpcGc0eb1zn+HOm3sraeCq5IZkWJRRsSKMiJWlBGxooyI+EYNt5Sf7x+HY9sgog70Hmt3NSIiIiIiImVKDbeUj9//Z45MjgMu+xdU0zE5IiIiIiLi39Rw28ThcBASEoLD4bC7lLK35Qf45lHzer/x0GKgvfVUElUqI1IiyohYUUbEijIiVpQREd9olHIpW4f+gnf6QUYytLkeBk0CvWGLiIiISAH0WV38jbZw28Tj8XDw4EE8Ho/dpZSd1KMw4zqz2a5/AVz+uprtYqgSGRGfKCNiRRkRK8qIWFFGRHyjhtsmHo+HrVu3+u+bV3YmfHILHN0K0Q3gug8gINjuqioVv8+I+EwZESvKiFhRRsSKMiLiGzXcUvoMA74eDdsXQ1A4XP8xhMfaXZWIiIiIiEi5qlAN96RJk0hMTCQyMpLIyEi6dOnCN998470/PT2dkSNHUqNGDcLDwxk8eDAHDhywsWIp0G//hRXvAg64ZgrUbmV3RSIiIiIiIuWuQjXc8fHxPP/88yxfvpxly5bRp08frrzySv744w8AHnjgAf7v//6PmTNn8uOPP7J3716uvvpqm6suGYfDQVRUlP+N+LjjF/ju1Dm2L3oGzhlgbz2VmN9mREqNMiJWlBGxooyIFWVExDcVfpTymJgYXnrpJa655hpiY2OZPn0611xzDQAbNmwgISGBX375hQsuuKBIy9PIh2XI44a3e8L+tZA4FK76rwZJExEREZEi02d18TcBdhdQGLfbzcyZMzl58iRdunRh+fLlZGVl0a9fP+88LVu2pEGDBmdtuDMyMsjIyPDeTklJASA7O5vs7GwAnE4nTqcTj8eTZ0CInOlut5vc30sUNt3lcuFwOLzLzT09Z51yeDweDhw4QN26dTnzO4+AgAAMw8gzv8PhwOVy5auxsOl2rJNj5fu49q/FCI7E3fdJyHVfZV2ns00v63UC2L17N3Fxcd7blX2d/PF1snOdsrKy2L9/P3FxcbhcLr9YJ398nexcJ4/Hw/79+6lXrx4ul8sv1il3jf7yOtm5ToZhcODAAeLi4vIsozKvkz++Tnauk8fjYc+ePXk+j5TlOmlwNvE3Fa7hXrt2LV26dCE9PZ3w8HC++OILWrVqxapVqwgKCiI6OjrP/LVr12b//v2FLu+5557jySefzDd95cqVhIWFARAbG0vTpk3Ztm0bhw4d8s4THx9PfHw8f/31F8nJyd7pTZo0oVatWqxbt460tDTv9JYtWxIdHc3KlSvzvNEkJiYSFBTEsmXLvNNy3uSio6O9u8yD+QbWqVMnkpOT2bBhg3d6SEgIbdq04fDhw2zdutU7PSoqioSEBPbu3cvu3bu908t7nVzZJ2mzaAIuIKvrg6xYvx3YXqnXKUfHjh3JzMxkzZo15fo6xcXF8ddff7Fnzx7vblyVfZ388XWyc53Wr19PUlISe/bsITQ01C/WyR9fJzvXyTAMkpKSCAsLo0aNGn6xTv74Otm5ThERERw/fhy3282+ffv8Yp388XWyc51cLhd//PFHns8jZblO1apVQ8SfVLhdyjMzM9m5cyfJycl8+umnvPPOO/z444+sWrWKW2+9Nc/WaoDOnTvTu3dvXnjhhQKXV9AW7vr163PkyBHvbip2fOPpdrtZsWIFHTt2zHdMTGX5xjP3OjnnT8D5y38wajSDu37GjavSr5PV9LJeJ4/Hw9KlS2nfvr23hsq+Tv74Otm5TpmZmaxYsYL27dsTEBDgF+vkj6+TneuU+39NYGCgX6xT7hr95XWyc53cbjcrV66kffv23q2XlX2d/PF1snOdsrOzWbZsWZ7PI2W5TidOnKB69erapVz8RoXbwh0UFESzZs0A6NChA0uXLuX111/nuuuuIzMzk6SkpDxbuQvaDSq34OBggoPzn/85ICCAgIC8q5/zJnSmnDeVok4/c7mFTXc4HDgcjgLnL2x6YTUWd3qprtORLebI5IBjwLMQEFxgsCrVOhVxelmuk8fj8f4zOvM5Kus6gf+9TmDfOuV8QHK5XPm+lClq7RVtnfzxdbJ7nXL+1xS39sKmV4R1sqqxuNO1Tub04sxfWdbJH1+n8l6nnPeQgj6PlMU6FTSPSGVW4RPt8XjIyMigQ4cOBAYGMn/+fO99GzduZOfOnXTp0sXGCkvG6XQSGxvrH28q3z8Onixo2heaX2R3NX7DrzIiZUIZESvKiFhRRsSKMiLimwq1hXvMmDFcfPHFNGjQgOPHjzN9+nQWLlzId999R1RUFLfddhsPPvggMTExREZGMmrUKLp06VLkEcorEqfTSdOmTe0uw3dbFsDGr8HhgoHPaVTyUuQ3GZEyo4yIFWVErCgjYkUZEfFNhfqq6uDBg9xyyy20aNGCvn37snTpUr777jv69+8PwGuvvcZll13G4MGD6dGjB3FxcXz++ec2V10yHo+HLVu2VO6RGN3Z8O0Y83rn2yG2hb31+Bm/yIiUKWVErCgjYkUZESvKiIhvKtQW7smTJ5/1/mrVqvHGG2/wxhtvlFNFZcfj8XDo0CEaNmxYeXfRWT4VDq2HkOrQ81G7q/E7fpERKVPKiFhRRsSKMiJWlBER3+ivRkom9SgseNa83vsxCI2xtx4REREREZEKRg23lMyPL0DaUYhNgA632l2NiIiIiIhIhaOG2yZOp5P4+PjKuWvOoY3w+//M6wOfBVeFOjLBb1TqjEi5UEbEijIiVpQRsaKMiPjGYeQ+630VkJKSQlRUFMnJyURGRtpdTuVjGPDB1bDlB2hxCVw/w+6KRERERMRP6LO6+Bt9VWUTt9vN+vXrcbvddpdSPH99azbbriC46Bm7q/FrlTYjUm6UEbGijIgVZUSsKCMivlHDbRPDMEhOTqZS7WCQnXH6NGAX3A01dE7GslQpMyLlShkRK8qIWFFGxIoyIuIbNdxSdL++Cce2QXgc9BhtdzUiIiIiIiIVmhpuKZrj+2HRy+b1fhMgOMLWckRERERERCo6Ndw2cTqdNGnSpPKM+DhvAmSegHodIfE6u6upEipdRqTcKSNiRRkRK8qIWFFGRHyj8znZxOl0UqtWLbvLKJrdy2D1qdHIL34R9IZbLipVRsQWyohYUUbEijIiVpQREd+oc7KJ2+1m9erVFX/ER48HvnnEvN7mBojvYG89VUilyYjYRhkRK8qIWFFGxIoyIuIbNdw2MQyDtLS0ij/i45qPYM9yCAqHfuPtrqZKqTQZEdsoI2JFGREryohYUUZEfKOGWwqXcdw8dhugx8MQEWdrOSIiIiIiIpWJGm4p3KKX4cQBiGkCF9xldzUiIiIiIiKVihpum7hcLlq2bInL5bK7lIId3WqedxtgwHMQEGxvPVVQhc+I2E4ZESvKiFhRRsSKMiLiG59GKf/uu++YPHkyW7du5dixY/mO7XA4HGzZssWnAv2Vw+EgOjra7jIK9/0T4M6Epn3gnAF2V1MlVfiMiO2UEbGijIgVZUSsKCMivilxw/3SSy/xj3/8g9q1a9O5c2fOO++80qzL72VnZ7Ny5UratWtHQEAFOzvbtkWw4StwuGDAs+Bw2F1RlVShMyIVgjIiVpQRsaKMiBVlRMQ3Jf6ref311+nTpw9ff/01gYGBpVlTlVEhT6/gccO3Y83rHUdArQR766niKmRGpEJRRsSKMiJWlBGxooyIlFyJj+E+duwY11xzjZptf7PyfTiwFqpFQa8xdlcjIiIiIiJSaZW44e7cuTMbN24szVrEbunJMP9p83rPf0BYDXvrERERERERqcRK3HC/+eabfP7550yfPr0066kyXC4XiYmJFWvEx8WvQOphqNEcOt9udzVVXoXMiFQoyohYUUbEijIiVpQREd84jDOHFi+ixMREjh49yr59+wgPDyc+Pj7fH6LD4WD16tWlUmhpSUlJISoqiuTkZCIjI22rwzAM3G43LpcLR0UYlOzoVnjjfHNk8hs+0cjkFUCFy4hUOMqIWFFGxIoyIlbKOyMV5bO6SGkp8RbumJgYmjdvTo8ePWjfvj21atWiRo0aeS4xMTGlWatfcbvdLFu2rOIMQjF33OnTgDW/yO5qhAqYEalwlBGxooyIFWVErCgjIr4p8SjlCxcuLMUyxFbbFsP6/wOHU6cBExERERERKSUl3sItfsLjhm9PjUau04CJiIiIiIiUGp/OXu92u/nggw+YM2cOO3bsAKBhw4Zcdtll3HjjjRpcoTJY+YF5GrDgKOg11u5qRERERERE/EaJB01LTk5mwIABLF26lIiICJo0aQLAtm3bSElJoXPnznz33XcVbrCDijIQQ4UYpCQrHV5PhBMHzF3Ju4y0pw4pUIXIiFRoyohYUUbEijIiVjRomohvSrxL+WOPPcby5cv5z3/+w6FDh1ixYgUrVqzg4MGDTJw4kWXLlvHYY4+VZq1+JzMz094C1nxkNtuR9aCTTgNWEdmeEanwlBGxooyIFWVErCgjIiVX4ob7iy++4O677+buu+8mMDDQOz0wMJC77rqLu+66i88++6xUivRHbrebNWvW2Dfio8cNP//HvN5lJAQE2VOHFMr2jEiFp4yIFWVErCgjYkUZEfFNiRvuI0eO0KJFi0Lvb9myJUePHi3p4qWsbZgDRzZDtShof4vd1YiIiIiIiPidEjfczZo1Y/bs2YXeP3v2bJo2bVrSxUtZMgxY8i/zeqfbITjC1nJERERERET8UYkb7rvvvpvvv/+eSy65hO+//57t27ezfft2vvvuOy699FLmzp3LPffcU5q1+h3bRnHfsQT2LAdXMJz/d3tqkCLRSP9iRRkRK8qIWFFGxIoyIlJyJR6lHGDChAk8//zzZGVl5ZkeGBjImDFjGD9+vM8FljaNfAh8OAQ2fW+ed/uy1+yuRkREREQE0Gd18T8+NdwAhw8fZt68eXnOw92vXz9q1qxZKgWWtoryR2wYBsnJyURFRZXvaTgO/AGTuoLDCfcsgxra7b+isi0jUmkoI2JFGREryohYKe+MVJTP6iKlJcDXBdSsWZOhQ4eWRi1VitvtZsOGDXTs2JGAAJ9fhqJb8m/zZ8IVarYrONsyIpWGMiJWlBGxooyIFWVExDdF/qvZuXMnAA0aNMhz20rO/FIBJO2CdZ+a17vdZ28tIiIiIiIifq7IDXejRo1wOBykpaURFBTkvW1F5+yrQH59EzzZ0LgH1GtvdzUiIiIiIiJ+rcgN95QpU3A4HAQGBua5LSXjcDgICQkpv99h6lFY/q55XVu3K4Vyz4hUOsqIWFFGxIoyIlaUERHf+DxoWmVTZQdi+PElWPAM1D4P7lwMetMUERERkQqmyn5WF79V4vNwjxgxgt9++63Q+3///XdGjBhR0sX7PY/Hw8GDB/F4PGX/ZFlp8Nt/zevd7lOzXUmUa0akUlJGxIoyIlaUEbGijIj4psQN97Rp09iyZUuh92/bto133323pIv3ex6Ph61bt5bPm9eqDyH1MEQ1gNZXlf3zSako14xIpaSMiBVlRKwoI2JFGRHxTYkbbit79+4lJCSkrBYvRWUY8Ntb5vWu94BLp3MQEREREREpD8Xqvr788ku+/PJL7+23336befPm5ZsvKSmJefPm0alTJ98rFN/s+h0O/wWBodDmerurERERERERqTKK1XD/+eefzJw5EzBHLPztt99Yvnx5nnkcDgdhYWH06NGDV199tfQq9TMOh4OoqKiyH/Fx5fvmz9ZXQTUNPFGZlFtGpNJSRsSKMiJWlBGxooyI+KbEo5Q7nU4++OADbrjhhtKuqUxVqZEPM07AKy0g8wTc+g007Gp3RSIiIiIihapSn9WlSijxMdwej6fSNdsVicfjYffu3WU7AMWfX5rNdkxTaNCl7J5HykS5ZEQqNWVErCgjYkUZESvKiIhvymzQNDm7cnnzytmdvN1NOhVYJaR/cGJFGREryohYUUbEijIi4hufGu5vvvmG/v37U6NGDQICAnC5XPkuYpPDm2HnL+BwarA0ERERERERG5S44f7ss8+47LLLOHDgAEOHDsXj8XD99dczdOhQQkJCSExMZNy4caVZqxTHqg/Mn836Q2Qde2sRERERERGpgkrccD/33HN07tyZlStX8uSTTwIwYsQIPvzwQ9atW8e+ffto3LhxqRXqb5xOJ7GxsTidZbBXvzsbVs0wr7e7qfSXL+WiTDMifkEZESvKiFhRRsSKMiLimxL/5fz5558MHToUl8tFQIB5drGsrCwAGjVqxN13380LL7xQOlX6IafTSdOmTcvmzWvzPDixH0JrwjkDS3/5Ui7KNCPiF5QRsaKMiBVlRKwoIyK+KfFfTmhoKEFBQQBER0cTHBzMvn37vPfXrl2bbdu2+V6hn/J4PGzZsqVsBqDIGSytzVAICCr95Uu5KNOMiF9QRsSKMiJWlBGxooyI+KbEDXeLFi34888/vbfbtm3L+++/T3Z2Nunp6UyfPp0GDRqUSpH+yOPxcOjQodJ/8zpxEP761rze9sbSXbaUqzLLiPgNZUSsKCNiRRkRK8qIiG9K3HBfddVVfPnll2RkZADw2GOPsXDhQqKjo4mNjWXx4sX84x//KLVCpYjWfAyebKjXAWq3srsaERERERGRKiugpA8cPXo0o0eP9t6+7LLLWLhwIZ9//jkul4tLL72U3r17l0qRUkSGAStPjU7e7mZ7axEREREREaniStRwZ2Rk8N1339GoUSMSExO90y+88EIuvPDCUivOnzmdTuLj40t3AIo9y+HQBggIgXOvLr3lii3KJCPiV5QRsaKMiBVlRKwoIyK+KdFfTlBQEEOGDOHnn38u7XqqjDJ581rxnvmz1ZVQLar0liu20D84saKMiBVlRKwoI2JFGRHxTYn+chwOB82bN+fw4cOlXU+V4Xa7Wb9+PW63u3QWmHkS1n1uXm+v3cn9QalnRPyOMiJWlBGxooyIFWVExDcl/qpq7NixTJw4kY0bN5ZmPVWGYRgkJydjGEbpLPDP2ZB5HKo3hobdSmeZYqtSz4j4HWVErCgjYkUZESvKiIhvSjxo2q+//kqNGjU499xz6dWrF40aNSIkJCTPPA6Hg9dff93nIqUIVn1o/mx7Izgc9tYiIiIiIiIiJW+4J06c6L0+f/78AudRw11Oju2A7YsBB7S93u5qREREREREBB8abo/HU5p1VDlOp5MmTZqUzgAUaz42fzbuAVHxvi9PKoRSzYj4JWVErCgjYkUZESvKiIhvSvyXs3PnTtLS0gq9Py0tjZ07d5Z08X7P6XRSq1Yt39+8DANWTTevt73B98Kkwii1jIjfUkbEijIiVpQRsaKMiPimxH85jRs35osvvij0/tmzZ9O4ceOSLt7vud1uVq9e7fuIjzt/hWPbICgcEi4vneKkQii1jIjfUkbEijIiVpQRsaKMiPimxA231UiFWVlZ+ibsLAzDIC0tzfcRH1ef2rrd6koICvO9MKkwSi0j4reUEbGijIgVZUSsKCMivinWMdwpKSkkJSV5bx85cqTA3caTkpL46KOPqFOnjs8FyllkpcEfs8zrbTRYmoiIiIiISEVSrIb7tdde46mnngLMEcjvv/9+7r///gLnNQyDZ555xucC5Sw2zIGMFIhqoHNvi4iIiIiIVDDFargvuugiwsPDMQyDRx55hOuvv5727dvnmcfhcBAWFkaHDh3o2LFjqRbrT1wuFy1btsTlcpV8ITmDpbUZCtp93++USkbErykjYkUZESvKiFhRRkR8U6yGu0uXLnTp0gWAkydPMnjwYM4999wyKczfORwOoqOjS76AlH2wdYF5vc3QUqlJKhafMyJ+TxkRK8qIWFFGxIoyIuKbEm8WHT9+vJptH2RnZ7N06VKys7NLtoA1H4PhgfoXQI2mpVucVAg+Z0T8njIiVpQRsaKMiBVlRMQ3xdrCfaZjx44xY8YMtm7dyrFjx/KNXuhwOJg8ebJPBfqzEp9ewTBg9QzzelsNlubPdAoOsaKMiBVlRKwoI2JFGREpuRI33N999x3XXHMNJ0+eJDIykurVq+ebx+Fw+FScFGLvSji0AQKqQeur7K5GREREREREClDihvuhhx4iLi6Ozz//nPPOO680axIrOVu3W14G1aLsrUVEREREREQKVOJjuDdv3sy9996rZruEXC4XiYmJxR/xMTsT1n5qXtfu5H6txBmRKkMZESvKiFhRRsSKMiLimxI33M2bN+f48eOlWUuVExQUVPwHbfoO0o5CRB1o0rv0i5IKpUQZkSpFGREryohYUUbEijIiUnIlbrifeeYZ3nzzTbZv316K5VQdbrebZcuWFX8QilWndidPvBac+qbRn5U4I1JlKCNiRRkRK8qIWFFGRHxT4mO458+fT2xsLAkJCfTv35/69evn29XE4XDw+uuv+1yknHLysLmFG6DNDfbWIiIiIiIiImdV4oZ74sSJ3utfffVVgfOo4S5lf3wBnmyo2w5qtbS7GhERERERETmLEjfcHo+nNOuQoti91PzZ4hJ76xARERERERFLDsMwDLuLKE8pKSlERUWRnJxMZGSkbXUYhoHb7cblchX9fOVvnG+ef/uGmXDORWVboNiuRBmRKkUZESvKiFhRRsRKeWekonxWFyktJR40Lcevv/7Kc889xwMPPMCmTZsASE1NZcWKFZw4ccLnAv1ZZmZmMWZOhcN/mdfrJJZNQVLhFCsjUiUpI2JFGREryohYUUZESq7EDXdmZiZXX3013bp147HHHuPf//43u3btMhfqdHLRRRfp+O2zcLvdrFmzpugjPh74AwwPhNeGiLiyLU4qhGJnRKocZUSsKCNiRRkRK8qIiG9K3HA/8cQTfPXVV0yaNImNGzeSe8/0atWqMWTIEL788stSKVKAfavMn3Xa2FqGiIiIiIiIFE2JG+4ZM2Zw1113cccddxATE5Pv/oSEBLZu3epTcZLLvtXmTzXcIiIiIiIilUKJG+6DBw9y3nnnFXq/y+UiNTW1pIuvEs48b/lZqeGukoqVEamSlBGxooyIFWVErCgjIiVX4tOC1a9fnw0bNhR6/5IlS2jWrFlJF+/3AgIC6NSpU9Fmzs6Ag+vN62q4q4xiZUSqJGVErCgjYkUZESvKiIhvSryF+4YbbuCtt97il19+8U7LOVXA//73Pz755BNuueUW3yv0U4ZhkJSURJHOynZwPXiyIKQ6RNUv++KkQihWRqRKUkbEijIiVpQRsaKMiPimxA33Y489RteuXenRowe9e/fG4XDwwAMP0KBBA/7+978zcOBAHnjggdKs1a+43W42bNhQtBEfc3Ynj0sEnSOzyihWRqRKUkbEijIiVpQRsaKMiPimxA13UFAQ3377LVOnTqVJkya0bNmSjIwMEhMTmTZtGv/3f/+n4z1Ky/415k/tTi4iIiIiIlJplPgYbjB3Ib/pppu46aabSqseKYgGTBMREREREal0SryF++jRo6xZs6bQ+9euXcuxY8eKtcznnnuOTp06ERERQa1atRg0aBAbN27MM096ejojR46kRo0ahIeHM3jwYA4cOFCidbCTw+EgJCTEe9x7odzZsH+deb1O2zKvSyqOImdEqixlRKwoI2JFGREryoiIbxxGCUdAGDZsGBs3buTXX38t8P6uXbuSkJDA5MmTi7zMgQMHMnToUDp16kR2djZjx45l3bp1/Pnnn4SFhQFw1113MWfOHKZNm0ZUVBT33HMPTqeTJUuWFOk5UlJSiIqKIjk5mcjIyCLXZpuD6+HNCyAoHP6xC5wl/o5ERERERKRCq3Sf1UUslLh7++GHH7jiiisKvf/yyy9n3rx5xVrmt99+y/Dhw2ndujVt2rRh2rRp7Ny5k+XLlwOQnJzM5MmTefXVV+nTpw8dOnRg6tSp/Pzzz4U2/hWVx+Ph4MGDeDyes8+Ye8A0NdtVSpEzIlWWMiJWlBGxooyIFWVExDcl7uAOHTpEzZo1C72/Ro0aHDx4sKSLB8wGGyAmJgaA5cuXk5WVRb9+/bzztGzZkgYNGuQ5PVll4PF42Lp1a9Ebbh2/XeUUOSNSZSkjYkUZESvKiFhRRkR8U+JB0+rUqcPKlSsLvX/58uXExsaWdPF4PB7uv/9+unXrxrnnngvA/v37CQoKIjo6Os+8tWvXZv/+/QUuJyMjg4yMDO/tlJQUALKzs8nOzgbA6XTidDrxeDx53kxyprvd7jznHixsusvlwuFweJebezqQ53QKOY81DCPf/AEBARiGgdvtxrV3NQ7AU/tcnKd+L7lrdDgcuFyuQmsvz3U62/Tc62RVu9bJnA7kW35lXyd/fJ3sXqecn/60Tla1a52Kvk65/9cAfrFOuWv0l9fJznXKmcfj8eR53sq8Tv74Otm5TjnvIbnrLMt1UmMv/qbEDfegQYN44403uPjii/PtWv7ll18ydepU7rrrrhIXNnLkSNatW8dPP/1U4mWAORDbk08+mW/6ypUrvceFx8bG0rRpU7Zt28ahQ4e888THxxMfH89ff/3l3doO0KRJE2rVqsW6detIS0vzTm/ZsiXR0dGsXLkyzxtNYmIiQUFBLFu2zDst500uPT2dP/74wzvd5XLRqVMnkpOT2bD+TzruXUkAsOlEKC2Aw4cPs3XrVu/8UVFRJCQksHfvXnbv3u2dbsc6AXTs2JHMzMw8A+rlWacNG7zTQ0JCaNOmjdapkHWKi4vj5MmTrFixwjtQSWVfJ398nexcp/Xr15OUlMSKFSsIDQ31i3Xyx9fJznUyDIOkpCRSUlKoUaOGX6yTP75Odq5TREQEAPv27WPfvn1+sU7++DrZuU4ul8v7vybn80hZrlO1atUQ8SclHjQtOTmZ7t278+eff9KmTRvvVuh169axevVqEhIS+Omnn/JtjS6Ke+65hy+//JJFixbRuHFj7/QffviBvn37cuzYsTzLbdiwIffffz8PPPBAvmUVtIW7fv36HDlyxDsQgx3feHo8HjZv3sw555yTr2bvt4OHNhHwZieMgGp4Ht2JKzC4Qn3jeeY6nW16ZfkWtyKtk2EYbNiwgebNm3u3eFf2dfLH18nOdcrKymLTpk00b94cl8vlF+vkj6+Tnevk8XjYtGkTLVq0ICAgwC/WKXeN/vI62blOHo+HLVu20KxZszyjUFfmdfLH18nOdXK73WzcuDHP55GyXKcTJ05QvXp1DZomfqPEDTfAyZMnefHFF/n888/ZsmULAE2bNmXw4ME8/PDD3i3IRWUYBqNGjeKLL75g4cKFNG/ePM/9ycnJxMbGMmPGDAYPHgzAxo0badmyJb/88gsXXHCB5XNUqpEP130On94K9TrA7T/YXY2IiIiISJmqVJ/VRYqgxLuUA4SFhfHkk08WuMt2SYwcOZLp06fz5ZdfEhER4T0uOyoqipCQEKKiorjtttt48MEHiYmJITIyklGjRtGlS5ciNdsVicfjYe/evdStW9f7bWE+GjCtSitSRqRKU0bEijIiVpQRsaKMiPimQv3VTJo0ieTkZHr16kWdOnW8l48//tg7z2uvvcZll13G4MGD6dGjB3FxcXz++ec2Vl0yHo+H3bt3n31gCDXcVVqRMiJVmjIiVpQRsaKMiBVlRMQ3Pm3hXr9+PVOnTmXr1q0cO3aMM/dOdzgczJ8/v8jLK8re7dWqVeONN97gjTfeKHa9lYphqOEWERERERGpxErccL///vvceuutBAYG0qJFC6pXr55vHh8OD5eUPZB2FJwBUKuV3dWIiIiIiIhIMZW44Z4wYQLt2rXjm2++oWbNmqVZU5XgdDqJjY21Pn47NgECgsuvMKkwLDMiVZ4yIlaUEbGijIgVZUTENyX+y9m7dy8jRoxQs11CTqeTpk2basA0KZRlRqTKU0bEijIiVpQRsaKMiPimxH85iYmJ7N27tzRrqVJyzntZ6AAUarirPMuMSJWnjIgVZUSsKCNiRRkR8U2JG+5XX32VyZMn8/PPP5dmPVWGx+Ph0KFDarilUJYZkSpPGREryohYUUbEijIi4psSH8P9wgsvEBUVxYUXXkirVq1o0KABLpcrzzwOh4Mvv/zS5yKrnOMH4Pg+wAFx59pdjYiIiIiIiJRAiRvuNWvW4HA4aNCgASdOnODPP//MN4/D4fCpuCpr/xrzZ83mEBRmby0iIiIiIiJSIiVuuLdv316KZVQ9TqeT+Pj4ggeg0O7kgkVGRFBGxJoyIlaUEbGijIj4psQNt/gm582rQGq4BYuMiKCMiDVlRKwoI2JFGRHxjc9fVf3444888sgjXHfddVx33XU88sgj/Pjjj6VRm19zu92sX78et9ud/0413IJFRkRQRsSaMiJWlBGxooyI+KbEW7gzMzO5/vrrmTVrFoZhEB0dDUBSUhKvvPIKV111FTNmzCAwMLC0avUrhmGQnJyMYRh570g7Bkk7zOtxieVfmFQYhWZE5BRlRKwoI2JFGREryoiIb0q8hfvJJ5/kiy++4KGHHmLfvn0cPXqUo0ePsn//fkaPHs3nn3/OU089VZq1Vg37Tg2YVr0RhETbWYmIiIiIiIj4oMQN9/Tp0xk2bBgvvvgitWvX9k6vVasWL7zwArfccgvvv/9+qRRZpeRs3a7R3N46RERERERExCclbrj37dvH+eefX+j9559/Pvv37y/p4v2e0+mkSZMm+Ud8PHnY/BkWW/5FSYVSaEZETlFGxIoyIlaUEbGijIj4psR/OfHx8SxcuLDQ+3/88UeNaHgWTqeTWrVq5X/zSj1i/gyrWf5FSYVSaEZETlFGxIoyIlaUEbGijIj4psR/OcOGDeOTTz7hzjvvZOPGjbjdbjweDxs3buSuu+5i5syZDB8+vBRL9S9ut5vVq1fnH/HRu4VbDXdVV2hGRE5RRsSKMiJWlBGxooyI+KbEo5SPHTuWLVu28Pbbb/O///3P+62Xx+PBMAyGDRvG2LFjS61Qf2MYBmlpaflHfEw91XCHquGu6grNiMgpyohYUUbEijIiVpQREd+UuOF2uVxMmzaNBx98kK+//podO8zBvho2bMgll1xCYqJOaVUi2sItIiIiIiLiF4rVcKenp3P//ffTunVrRo0aBUBiYmK+5vrf//43//3vf3n99dd1Hu7iyjmGW1u4RUREREREKrViHcP99ttvM23aNC699NKzznfppZcyZcoU3nnnHZ+K82cul4uWLVvicrlOTzQMOHnIvB5Ww57CpMIoMCMiuSgjYkUZESvKiFhRRkR8U6yG+5NPPmHw4ME0adLkrPM1bdqUIUOGMGPGDJ+K82cOh4Po6GgcDsfpiZknITvdvK4t3FVegRkRyUUZESvKiFhRRsSKMiLim2I13GvXrqV79+5Fmrdr166sWbOmREVVBdnZ2SxdupTs7OzTE3MGTAuoBkFh9hQmFUaBGRHJRRkRK8qIWFFGxIoyIuKbYjXcmZmZBAUFFWneoKAgMjIySlRUVZH/lGC5jt/Wt4hCARkROYMyIlaUEbGijIgVZUSk5IrVcNetW5d169YVad5169ZRt27dEhVVZeVs4dbx2yIiIiIiIpVesRrufv368d5773Hw4MGzznfw4EHee+89+vfv71NxVc5JnYNbRERERETEXxSr4X700UdJT0+nT58+/PbbbwXO89tvv9G3b1/S09N5+OGHS6VIf+RyuUhMTMw74mOqzsEtpxWYEZFclBGxooyIFWVErCgjIr4p1nm4mzRpwieffML1119P165dadKkCeeddx4REREcP36cdevWsWXLFkJDQ/noo49o2rRpWdXtF/IdD68t3HKGoo6ZIFWXMiJWlBGxooyIFWVEpOSKtYUbzHNsr1mzhjvuuIP09HRmzZrF+++/z6xZs0hNTeX2229n9erVXH755WVRr99wu90sW7Ys7yAUqacGTdMx3EIhGRHJRRkRK8qIWFFGxIoyIuKbYm3hztGoUSMmTZrEpEmTOH78OCkpKURGRhIREVHa9VUt2sItIiIiIiLiN0rUcOcWERGhRru06BhuERERERERv1HsXcqlDGkLt4iIiIiIiN9wGIZh2F1EeUpJSSEqKork5GQiIyNtq8MwDNxuNy6XC4fDYU58th5knoBRK6CGBpyr6grMiEguyohYUUbEijIiVso7IxXls7pIadEWbhtlZmaevpGVbjbbAKEaNE1MeTIiUgBlRKwoI2JFGREryohIyanhtonb7WbNmjWnR3zMOX7bGQjVouwrTCqMfBkROYMyIlaUEbGijIgVZUTEN2q4K4qTuQZM0y5dIiIiIiIilZ4a7ooiVQOmiYiIiIiI+BM13DZyuVynb5w8Yv4M0/HbclqejIgUQBkRK8qIWFFGxIoyIlJyGqW8ovjlDfhuLJx7DVwz2e5qRERERETKXYX9rC5SQtrCbRPDMEhKSsL7fUfuY7hFKCAjImdQRsSKMiJWlBGxooyI+EYNt03cbjcbNmzIP0q5juGWU/JlROQMyohYUUbEijIiVpQREd+o4a4ovFu4dQy3iIiIiIiIP1DDXVGc1BZuERERERERf6KG2yYOh4OQkBAcOefcTtUx3JJXvoyInEEZESvKiFhRRsSKMiLiG41SXlE81wAykmHkUog9x+5qRERERETKXYX9rC5SQtrCbROPx8PBgwfxeDyQnWk226At3OKVJyMiBVBGxIoyIlaUEbGijIj4Rg23TTweD1u3bjXfvFKPmBMdLqgWbWtdUnHkyYhIAZQRsaKMiBVlRKwoIyK+UcNdEXhPCRYDTr0kIiIiIiIi/kDdXUWgEcpFRERERET8jhpumzgcDqKioswRH3N2Kdfx25JLnoyIFEAZESvKiFhRRsSKMiLimwC7C6iqXC4XCQkJ5g3vFu4a9hUkFU6ejIgUQBkRK8qIWFFGxIoyIuIbbeG2icfjYffu3acGTdM5uCW/PBkRKYAyIlaUEbGijIgVZUTEN2q4bZLnzUvHcEsB9A9OrCgjYkUZESvKiFhRRkR8o4a7ItAWbhEREREREb+jhrsiOKlB00RERERERPyNGm6bOJ1OYmNjcTqduc7DrYZbTsuTEZECKCNiRRkRK8qIWFFGRHyjUcpt4nQ6adq0qXnjpHYpl/zyZESkAMqIWFFGxIoyIlaUERHf6Ksqm3g8HrZs2YInOwvSjpkTtYVbcvFmRIOUSCGUEbGijIgVZUSsKCMivlHDbROPx8OhQ4fwnDwMGIADQmPsLksqEG9G9A9OCqGMiBVlRKwoI2JFGRHxjRpuu6WeGjAtpDo4XfbWIiIiIiIiIqVGDbfNHDp+W0RERERExC+p4baJ0+kkPj4eR9qpLdw6flvOkJMRjQoqhVFGxIoyIlaUEbGijIj4Rn85NvG+eaUdNSeE1bC3IKlw9A9OrCgjYkUZESvKiFhRRkR8o78cm7jdbtavX4/nxEFzgrZwyxlyMuJ2u+0uRSooZUSsKCNiRRkRK8qIiG/UcNvEMAySk5N1Dm4pVE5GDMOwuxSpoJQRsaKMiBVlRKwoIyK+UcNtt1Qdwy0iIiIiIuKPAuwuoMrLabiLuIV788ETHDqeQVpWNqmZbtIy3aRlub3X07PMS1qWm/Qsz6mfbgKcDto3qE7nxjG0qR9NtUCdgkxERERERKQsqeG2idPppEmTJjiWn9qlPPTsg6YZhsE/56znnZ+2lfg5F2w8BECQy0mb+lF0bhxD58Y1aBMfRVRIIA6Ho9DHpmW6WbsnmRU7j7FixzHW7kmmRngQvVvUonfLWrSJj8blLPzxUnw5GdEgJVIYZUSsKCNiRRkRK8qIiG8cRhU7ICMlJYWoqCiSk5OJjIy0uxx4qRmcPAR3/gRx5xU627/m/cW/5m0CoGlsGKFBAYQE/X979x0fR3ktfPw3s12r3VVvlovcezcGm2YwGEIAE3oCF9K4NyGFkISQRhohhSQvIRAI3CSkwCUJCZCE6phqum25N9kqVq+r1fYyM+8fK60tS/IaXHZtn68/+5E8enb2zO7ZmTkzzzxjwmExkWM1pX53WJK/2/t/t1tMOKwqfeEE79b38G5dD53+6JD528wqxS4bJS5b/087JS4bXYEo6/f2sr21j4Q+cqoUOK2cPbmYZVNLOHNyMR6H5fDfGyGEEEIIcVLJun11IQ6TnOHOEE3T2LJ5E7NDPShw0Gu4f7umLlVsf/fi6dy4tOoDveYNS8ZhGAb13SHeq+vhnboe3qnrpskbJprQafKGafKGR3x+icvG/DH5zBuTx5zReTR5w7y8o4PXdnXSE4zxj+pm/lHdjElVmFLqYnqFmxkVbmZUeJhW7sJllyL8/dA0jS1btjBz5kxMJrkEQAwlOSLSkRwR6UiOiHQkR4Q4PFJwZ4hhGMT7OlGM/lssjNCl/K9rG/nBv7cB8OXzJn/gYnuAoihUFTmpKnJy1aLRAETiGp3+KB3+KJ3+CB3+KB19UTr8EZw2c6rIHpXnGNLt/IoFlcQ1nbX1Xl7e2cFLOzrY3RFgW2sf21r7eGLdvrZjC3OYUeHm4tkVnDe9FLNJuiYdjGEYhMNhGRVUjEhyRKQjOSLSkRwR6UiOCHF4pODOIHOsN/mLzQNm65C/P7e5ldv/vgmAT59RxefOmXhU4rBbTIwuyGF0Qc4Her7FpHLahEJOm1DINz40jZbeMJubfWxt6WNbSx/bWny0+CI0dIdo6A7x7OY2Kjx2PnbqWK49ZQwFzqHLLoQQQgghhBDHOym4M8gS8yV/cQ49u/3ark6+8Hg1ugFXLxzNNz407aCDmmWTijwHFXkOVswoS03rCcbY1tLHm3u6+Mt7jbT4Itz9wk5+ubqGS+dUcMOSccwc5cEwDLoCMXa1+9nZ5k/+bPcTjml8aFY5ly+oZFSeI4NLJ4QQQgghhBCHRgZNyxDDMAitexznv/8HKk+BT61K/W1dQw/X/e+7hOMaF80q595r551QI4BH4hr/3tTKH96sZ3OzLzV9Ukku3cEYPcHYiM9VFDh9YhFXLRzNedNLT+jbmxmGgc/nw+PxHDcHW8SxJTki0pEcEelIjoh0jnWOZMu+uhBHipzhzhBFUXDSP0CZszg1vbEnxI2/f49wXOOsycX8v6vnnlDFNiS7sF+xoJLL549i/d5e/vBmPc9ubqWmIwAki+pxhU4ml+YypdTF5DIX0bjOE+uaeKu2m9druni9pguPw8LKuRVcuXA0M0d5MrxUR56iKOTl5WU6DJHFJEdEOpIjIh3JEZGO5IgQh0cK7gxJJBK07dpIJQzqUv7oO3vxRxLMqfTw4HULsJpP3IHFFEVhwdh8FozN51sXTWP9Xi+V+TlMKM7FYR165vryBZXs7Q7xxLpGnljXRIsvwh/eauAPbzUwa5SHqxeN5tK5FSfMaOiJRILq6mrmzZuH2SxfVTGU5IhIR3JEpCM5ItKRHBHi8Mi3JoNM0d7kL/23BNN1g6c3NAPwmbMnDFt0nqhK3HYumFmett2YwhxuPX8KX1w+mTd2d/GXtY2s2trO5mYfm5t9/PCZ7Vw0u5xrFo1mwdj84757nKZpmQ5BZDnJEZGO5IhIR3JEpCM5IsQHJwV3BlkGRil3Jgvud+p6aPVFcNvNnD2lJHOBHQdMqsKZk4s5c3IxPcEYT1Y38/i7e6npCPDEuiaeWNfExJJcLps3iqUTi5g1ynPCdc0XQgghhBBCZDcpuDPIPDBKef8Z7qeqk2e3PzSr/IQeDOxIK3Ba+eTpVXxi6TjW7+3l8Xf38u9NrezuCHD3Czu5+4WduO1mTh1fyNKJRSydWMiE4tzj/uy3EEIIIYQQIrvJKOUZYhgGxgNLUTu2wnV/JzJ2GYvu/A/+aILHbzqVU8cPvVWYOHT+SJx/b2rl5R0dvFXbjT+SGPT3UreNU6oKmT8mj3lj8ple7s666+UNwyAcDuNwOOTggBiW5IhIR3JEpCM5ItI51jmSLfvqQhwpcoY7g5Rwd/KXnCJe2tGBP5qgwmPnlHEFmQ3sBOCyW7j2lDFce8oYNN1gS7OPNbu7eHNPF+/Ve2nvi/KvjS38a2MLAFazyqxRnlQBPmuUh1F5DtQMd0O3Wq0ZfX2R/SRHRDqSIyIdyRGRjuSIEB+cFNwZoiUSqIEuFABnEU/9J9md/NJ5ozJe5J1oTKrCnNF5zBmdx83LJhKJa6xv8LKuwUt1Yy/Ve714Q3HW9U+DOgBybWamlLmYWuZiarmbaWUuJpW6sJpUNMNANwx03UDTDTTDwGY24XEcuRHSNU1j7dq1LFy4UEYFFcOSHBHpSI6IdCRHRDqSI0IcHvnWZEq0D9VIdnPuVdy8vHMTAJfNG5XJqE4KdouJJROLWDIxee28YRjUd4dY3+ClutFL9d5eatoDBKKJ/YrwQ3PBjDK+uHwS08qlC5QQQgghhBAnOym4MyXUBYBhdfLMdi9xzWBauZvJpa4MB3byURSFqiInVUVOLl9QCUBc06nrCrK9tY8dbX529P9s9UVGmAcYBjy/tY3nt7Zx4cxk4T21TApvIYQQQgghTlZScGdKcN/12wOjk182ryKDAYn9WUwqk0tdTC51cel+00OxBLoBJkVBVZM/TaqCoijsavfzy9U1PLu5lee2tPHcljYumlXOF86dxJQyOZAihBBCCCHEyUZGKc8QY/u/Uf7yMaKl85jS8FUUBd66/VzKPPaMxSSOjJ1tfu5dXcMzm1uB5Nnvi2aV840PTaMiz3HI8zEMA03TMJlMMnKsGJbkiEhHckSkIzki0jnWOZIt++pCHCnZdR+kk0koeYa7Ne4EYMmEQim2TxBTylzc/7H5PH/LGXxoVhmGAf/e1Mp5v3iV362pQ9MP/RhXLBY7ipGKE4HkiEhHckSkIzki0pEcEeKDk4I7Q/RABwA7+pK3WVg5VwZLO9FMLXPz648t4NkvnMHCsfkEYxrf//c2Lvv1G2xp9qV9vqZpbNq0CU3TjkG04ngkOSLSkRwR6UiOiHQkR4Q4PFJwZ4jSf4a7PpKDzaxywcyyDEckjpbpFW7++t+n8cPLZuKym9nU5OPS+9/grme3E4olhrQPxRJsafbxzOZWanriGYhYCCGEEEIIcSTIoGmZEkyOUt5juFg+oxSX/cjdv1lkH1VV+NjisZw3rZTv/Xsbz2xq5aHXanlmUyvXnzaWNl+EPZ0BajuDNPeGBz13Ves6vv6haXKrMSGEEEIIIY4zUnBnSv8Z7h7cXCbdyU8aJW479390PlfM7+BbT22huTfMj5/bMaRdgdPK6HwHW5p9vFbTxev3vs5H5lXy5fMnv6+B18SJz2QyZToEkeUkR0Q6kiMiHckRIT64rBql/LXXXuPuu+9m3bp1tLa28uSTT7Jy5crU3w3D4Dvf+Q4PP/wwvb29LF26lAceeIBJkyYd8mtky8iH/nuX4OrZyueV2/n5N7+G1Sy9+082oViCB1/Zw442P1VFTiYU5zK+2Mn44lwKnMlr++u7gtz94k6e2ZQc8dxqVvn40nF89uyJeBzSK0IIIYQQJ5Zs2VcX4kjJqjPcwWCQOXPm8IlPfIKPfOQjQ/7+05/+lHvvvZc//OEPVFVV8e1vf5sVK1awbds27Pbja4TvRF8nANMnjpdi+ySVYzVz6/lTRvy7YRjkmePcd+08Pn3GeH707HbeqevhN6/W8vi7jVw6t4JJpS4mFucyqTSXQqdVbulykjEMA5/Ph8fjkc9eDEtyRKQjOSLSkRwR4vBkVcF94YUXcuGFFw77N8MwuOeee/jWt77FpZdeCsAf//hHSktLeeqpp7jmmmuOZaiHJRSN44h7QYHTZk3OdDgiS2maxo4dO1i4cCFzR+fx+E2n8vLODn783A52tQf441sNg9rn5ViYVJLLxJJcppW7mTnKw7QyNw6rdAM7Ue2fI2ZzVq3ORZaQHBHpSI6IdCRHhDg8x823pq6ujra2NpYvX56a5vF4WLx4MW+99dZxVXC/urWeC5Xk6NMzJlRlOBpxvFAUhXOmlnLW5BJe2NrGxsZedncEqOkI0OgN0RuK8169l/fqvannmFSFicW5zBjlZtYoD7MrPcwbnY+qyhFqIYQQQgghjrbjpuBua2sDoLS0dND00tLS1N+GE41GiUajqf/39fUBkEgkSCSSt2RSVRVVVdF1HV3XU20Hpmuaxv6Xuo803WQyoShKar77TwdS9y88pzJZ7EQVK6rVOaS92WzGMIxB9ztUFAWTyTQkxpGmH+tlSjddlun9LxMwZP6qqmJSVVZML+H8acWp6XEdaruC7GrtY1dHgG2tfWxp7qM7GGNnu5+d7X7+sb4ZgDmVHn542UxmVHjkczoBlmng54m0TOlil2U69GUa+DnQ5kRYpv1jPFE+p0wu00AbXdcHve7xvEwn4ueUyWUaWIfsH+fRXKb9fxfiRHDcFNwf1I9+9CO+973vDZleXV2N0+kEoLi4mAkTJlBXV0dnZ2eqTWVlJZWVlezatQufz5eaPn78eEpKStiyZQvh8L5bOE2dOpW8vDyqq6sHrWhmz56N1Wpl7dq1AAS63uau0RXEVRPPRaNs2bIl1dZkMrFo0SJ8Ph87duwbvdrhcDBnzhy6urqora1NTfd4PEybNo2WlhaamppS04/1Mg1YuHAhsViMTZs2yTId5jKVl5cTjUZZv3596pqpgy3TjIo8Qs27KM/XOCsfjGm5lFVNZVdXmBff20G9L8HWrjgbm3xcct8bXH9KJWcXhXBYFPmcjtNl2r59O36/n/Xr15OTk3NCLNOJ+DllcpkMw8Dv9+P3+ykoKDghlulE/JwyuUxutxuHw0FbWxstLS0nxDKdiJ9TJpfJbDantjUD+yNHc5mOt3GZhEgnq0Yp35+iKINGKa+trWXChAlUV1czd+7cVLuzzjqLuXPn8stf/nLY+Qx3hnv06NF0d3enRj481kc8vVvW8Optn0Ex4EPPrseiDB5t+ng54rn/MqWbLsuUHcvU3hfhh8/u5Nkt/T1G3DbuuGga508vQVGU43KZ0k2XZZJlkmWSZZJlkmWSZTp+likQCJCfny+jlIsTxnFTcBuGQUVFBV/5ylf48pe/DCSL55KSEh555JFDvoY7G241EGtvZ89ZZ6Mp4HrzOcbmj8tIHCK76bpOV1cXRUVFqOqRHcn+lZ0d3PH0Vvb2hAA4d2oJ371kBqMLco7o64ij62jmiDgxSI6IdCRHRDrHOkeyYV9diCMpq7qUBwIBdu/enfp/XV0dGzZsoKCggDFjxnDLLbdw5513MmnSpNRtwSoqKgbdq/t4YMnPB8BkQGd7vRTcYli6rlNbW0tBQcER38CdPaWEF79UyH0v7eY3r+1h9Y4OXqvpZHRBDqPyHKlHRZ6DUfnJ30vddrmFXZY5mjkiTgySIyIdyRGRjuSIEIcnqwrutWvXsmzZstT/b731VgBuuOEGHnnkEW677TaCwSA33XQTvb29nH766Tz//PPH3bUeitVKxK5ij+h4OxpgaqYjEicju8XEV1ZMYeW8Cr755BbeqeuhtjNIbWdw2PaKAkW5Nio8dso8dso9Diry7Ewv97B0YqHcm1MIIYQQQogDZFXBffbZZ3OwHu6KovD973+f73//+8cwqqMj6rJjj4QIdDRnOhRxkptY4uLxm06lsSdMkzdEU2+Ylt4wzd4wLb6BnxFiCZ1Of5ROf5SNTb5B85hS6uKzyyZw0axyzCY5+i2EEEIIIQRkWcF9MtHcOdAZItQ58i3NxMlNURQ8Hs8xOXOsKApjCnMYUzj8NdyGYdAdjNHmi9DSG6bVF0kV46/s7GRnu58vPr6Bn7+4i5vOHM8VCyqxW0xHPe6T3bHMEXF8khwR6UiOiHQkR4Q4PFk7aNrRki0DMbx7/Upc7+2k6bMXc94XfpqxOIQ4XL5QnD++Vc/v36ynJxgDoNhl45OnV/GxxWNw2S1p5iCEEEIIkZQt++pCHCnS9zNDKsfNRsnPZ0HJ/EyHIrKUrus0NTUNulVGNvLkWPj8uZNY87VlfOfi6ZR77HT6o/z4uR0s+fFL/Oi57bT6wulnJN634yVHROZIjoh0JEdEOpIjQhweKbgzpOQ7d9D7/36B+4orMh2KyFLH2wYux2rm40urePWry/jpFbMZX+zEH0nwm1drOeMnL/Olv2xga4sv/YzEITveckQce5IjIh3JEZGO5IgQh0eu4c4wTdcwy8cgTiBWs8pVC0dzxfxKVu/o4OHXa3m3rocnq5t5srqZpRML+fQZ4zltQiEWVUVV5ZowIYQQQghxYpJKL0MC8QBf2f0V/Lv8vPuxd7GZbJkOSYgjSlUVzpteynnTS9nY2MvDr9fy3JY23tjdzRu7u1PtFAXMqoJJVTCrKhaTQrnHwYSSXCYUO5lQnMuE4lzGFztlIDYhhBBCCHFckYI7Qyy7G7n1D156nQZdH+liVO6oTIcksoyqqhQXF6Oqx/+VH3NG53HfR+fT5A3x+zfq+ct7jQSiCQAMA+KaQVwzgGR3NW8ozrbWvkHzUBQYU5DDZ8+ewFULR8toqZxYOSKODskRkY7kiEhHckSIwyOjlGdIeNMm6q+6mi43uP79f8wtmZuxWIQ41hKaTiiuoWkGCd1A0w3imo6mG0QTOo09IfZ0BvofQXZ3BPCF46nnXzyngrsumykjoAshhBAnmGzZVxfiSJEz3Bmi5OUB4A5BR6gjs8GIrKTrOnV1dVRVVZ1wR5XNJhW3aeRlmlLmYjmlqf8bhkFPMMZf1jby8xd38a+NLWxq6uW+a+czq9JzLELOSidyjogjQ3JEpCM5ItKRHBHi8Mi3JkNUT7JIsCag29uS4WhENtJ1nc7OThkVFFAUhcJcG589eyJ//e/TGJXnoKE7xEceeIPfrqnjJOuokyI5ItKRHBHpSI6IdCRHhDg8UnBniJKTQ8KcfPt9HY0ZjkaI48eCsfk8+4UzWDGjlLhm8IN/b+PTf1yHNxjLdGhCCCGEEEIMIl3KM0RRFKJOG2ZfmGBHa6bDEeK44smx8OB1C/jT2w3c+e/t/Gd7O8t/8SoVeQ4MDAwjORibQbI7uttu4czJRZw3vYzJpbky4JoQQgghhDgmpODOEFVVUfM94AszxsjPdDgiC6mqSmVlpVwvNQJFUfiv08Yxf0w+n/+/auq6gnQf5Cz3u/U9/OzFXYwpyEndrmzh2HzMB1xLrukGfeE4vnAcp81MsSt7b9knOSLSkRwR6UiOiHQkR4Q4PDJKeQY1/s9niGzdStn3vovrnHMyGosQx7NwTOO9+h4Suo6iKCiw309o8oZZta2dNbu7iCX2XYOWn2NhapkbfzRZYPeG4vgjiUHzLnRamVLmYkqZi6llLqaUuZlcmkuOVY5XCiGEEEdaNu2rC3EkSMGdIZqmsWvXLiZPnozJZMpYHCJ7SY4cecFogtdrOnlxWzsv7eigNxQfsa3TaiIU1xhuDakocPrEIm49bzLzxmSuh4rkiEhHckSkIzki0jnWOZIt++pCHClyiiZDDMPA5/NhGAYJPYFJMcl1pWKQ/XNEHBlOm5kLZpZzwcxyEprO2gYvrb4weQ4rboeFvBwLHkfyYTGphGMaNR1+drT52dn/2NHmpysQ5fWaLl6v6eK86aV8+fzJTC079jsFkiMiHckRkY7kiEhHckSIwyMFdwYZhsGFT15Ie6idFy5/gfLc8kyHJMRJw2xSOXV84UHbOKwmZlfmMbsyb9D0vd0hfvVSDX9f38Sqbe38Z3s7l86p4JblkxlX5DyKUQshhBBCiOOJjH6QQZat2/jC/3Zy44sJOsIdmQ5HCHGIxhTmcPeVc3jxS2dx0axyDAOe2tDC8l+8yjee3MyOtj45EyCEEEIIIeQMd6aoqkppngdnQxQtAV2hrkyHJLKMqqqMHz9eRgXNYhNLcrn/Y/P5TLOPn724k1d2dvLYO3t57J29lHvsnD2lhGVTilk6sQinbfDqNq7p7O4IsKXZx9aWProCUc6cXMwFM8tw2y2H9PqSIyIdyRGRjuSISEdyRIjDI4OmZVBo7Voarrue1nxo//0dXDv12ozGI4Q4PO/V9/DgK3tYs7uL6H6joVtNKqdUFXBKVQGtvjBbW/rY0eYfNGJ6qq1Z5bxppVw6t4Kzp5RgNcsOjhBCiJNHNu2rC3EkyBnuDNE0jZqODqyAOwRbQp2ZDklkGU3T2LJlCzNnzpSRY48Ti8YVsOjGAiJxjbdqu3llRwcv7+xkb0+INbu7WLN7cE8Wl83MtAo3Mys85NrNPLu5ld0dAZ7Z3Mozm1vxOCxcNLuclXNHsXBsPqo6eGBFyRGRjuSISEdyRKQjOSLE4ZGCO0MMwyBitWIFnFHo8rdnOiSRZQzDIBwOy7XAxyG7xcSyKSUsm1LCdw2D2q4gr+zsZENjL6PzHcyo8DCjws2YgpxBRfSXlk9ia0sfT29o5ukNLXT4o6ku6qPyHFw6t4LL5o1iUqkLkBwR6UmOiHQkR0Q6kiNCHB4puDPIyMnBUBUU3SDQ1ZLpcIQQR4GiKEwozmVCce4htZ05ysPMUR5uv3Aab9d282R1M89vaaO5N8yvX9nDr1/Zw4wKNyvnjuJDM0uOwRIIIYQQQogPSgruTFJVDHcuSq+fGabRmY5GCJFFTKrC0olFLJ1YxJ0rZ/Kf7e08Vd3CKzs72NrSx9aWPu56bjvTCs1cEqnj3GllTCzJRVGU9DMXQgghhBDHhAyaliGGYeDz+ei58Ua07h4qf3UvjrlzMxaPyD4DOeLxeKSIEik9wRjPbG7l6epm1jZ4B/1tVJ6Ds6YUs2xKCUsmFA4ZGV2cfGQ9ItKRHBHpHOscyZZ9dSGOFCm4M8wwDNnACSE+kMaeEP/Z3s7LOzt5u7Z70KjnVpPK4vEFrJhRxvnTSylx2zMYqRBCCHFosm1fXYjDJQV3hiQSCaqrq5k3bx5ms5m4HkdFxaTK6I8i6cAcEeJA++dIXFd4q7aLV3Z28vLODhp7wql2igLzx+SzYkYpK2aUMbbQmcGoxbEk6xGRjuSISOdY50i27KsLcaTImjWDNE0D4Kp/XcWOnh08+qFHmVU8K8NRiWwykCNCjGQgRxxWE+dMLeWcqaUYhsGeziCrtrXzwtY2NjT2sq7By7oGL3c9u4OpZS7GFuag6aDpOpqR/JnQDAygwmNnQnEuE0uSj7GFTrkf+HFM1iMiHckRkY7kiBAfnBTcGeZ/9lk++UAtb1VqdCzryHQ4QogTgKIoqWL5M2dPoM0X4cVtbbywtY23a3vY0eZnR5v/kOdnVhXGFOYwqSSXKWVuppcnH5X5jiH3BhdCCCGEEPtIwZ1hmtfLmLogzTaFrlBXpsMRQpyAyjx2/uu0cfzXaePwBmO8VtNJIJrApCioqoJZVTD1PwwDGr0hdncE2NMRYHdHgGBMo7YzSG1nkBe2tqfmm2szM7XMxfQKNzMq3JwxqZiKPEcGl1QIIYQQIrtIwZ0hJpOJ2bNnE+9KFtmuELSF5Qy32GcgR0wmua5fDO+D5Ei+08qlc0cdcnvDMGjri7C7I8Cu9gA7WvvY1tpHTXuAQDTB2gbvoNHSp5e7OXdaCedOK2X2KI+cAc8wWY+IdCRHRDqSI0IcHim4M8hqtaLn5wPgDhlsCcsZbjGY1WrNdAgiyx3tHFEUhXKPg3KPgzMmFaemxzWd2s4g21v72N7ax9oGL9V7vWzrL8h/9dJuil02zplSwllTiplcmrxu3GKSa8GPNVmPiHQkR0Q6kiNCfHBScGeIpmmsXbuWWZ48ANwh6Ax1ZjYokVUGcmThwoUycqwYViZzxGJSmVLmYkqZi5XzkmfMe4IxXt7Rweod7by2q4tOf5S/rG3kL2sbgeS14GMLc5hQnMuEktzkz2InVUVO8nJkZ+5okPWISEdyRKQjOSLE4ZFvTYaZCpJnuF1h6ApKl3IhxPGrwGnl8gWVXL6gklhC5526blZv72BtQw+1nUFCMY09nUH2dAZhW/ug5+bnWKgqclJVlEtVUQ5VRbnMrvQwuiAnQ0sjhBBCCHH4pODOMFN/l3KTAYtyZ2Q4GiGEODKsZpUzJhWnuqHvfy34no4AezqD7O4IUNcVpK0vgjcUx7u3l/V7ewfNZ3JpLudNL2X5tFLmVObJNeFCCCGEOK5IwZ1hisWCqaAAxWzmi1NuynQ4QghxVIx0LThAMJqgvjtIfVeIuq4AtV3Js+Bbmn3sak8O1nb/y3sodtlYPq2E5dNKOW1CITlW2YQJIYQQIrsphmEYmQ7iWOrr68Pj8eDz+XC73RmLwzAMNE1LjfioKHLWRgy2f45IfojhnOg50huK8crOTlZtb+fVnclbmQ0wqQpTy1zMH5PP/LF5zBudz9jCnCHvg2EYBGMaPYEY4bjGxJJcTCfRWfITPUfE4ZMcEekc6xzJln11IY4UKbgzxDAMwuEwDocjtfKK63EALKolY3GJ7DFcjgixv5MpR6IJjXdqe1i1rZ2XdnTQ3Bse0qbQaWV2pQeA7mCM7kCMrkCUaEJPtZlS6uLW8ydz/vTSE/49g5MrR8QHIzki0jnWOZIt++pCHClScGdIIpEYNOLjf6/6b95qeYt7lt3DOWPOyVhcInscmCNCHOhkzpGW3jDr93pZ39DL+r1etrb4iGsjb87sluTtyCLxZPE9p9LDV1ZM4fSJRcPuQPrCcVZvb+e5LW34QnGuOWU0l8ypwHyc3dbsZM4RcWgkR0Q6xzpHsmVfXYgjRdasWaDn0Ue56rFNOCZpdJ0q9+IWQoh0KvIcVOQ5+PDsCgAicY2tLX1sa/FhM5sozLVSmGuj0GmlMNdKjtWMLxzn4ddq+d0bdWxs8nH9b99lcVUBX10xhYXjCugJxli1rY1nN7fx5p6uQQX8u/U93Lu6hpuXTWTlvFFyP3EhhBBCHBIpuLNAorWVij29VBYodIblXtxCCPF+2S0mFozNZ8HY/BHbeBwWvrJiCjcuHcevX97Dn99u4J26Hq548C0ml+aypzOIpu8rsieV5HLhzDIsJpXfvVFHfXeIrz6xiV+9tJubl03gsnmVWM1SeAshhBBiZFJwZ9DAgGmm/AIA3CFoCUnBLfYZyBEhRiI58v4V5dq44+LpfPKMKn61uoa/rWtiV3sAgBkVbi6cWcYFM8uZWJKbes4nTq/iz2838NBrteztCfG1v2/m3tW7+ejiMVTmOyhx2Sl12yh123Hahm5aYwmdvkicvnAcfyTB6IIcCpzWY7K8kiMiHckRkY7kiBAfnFzDnQV6n3yK1q9/nQ1VCmu+soz7z70/0yEJIcRJo64ryIZGLwvGFDCmMOegbUOxBI+9s5ffvFZLpz86bJtcm5kStw2TovQX2QnCcW1Iu5mj3Jw+sZgzJhWxYGw+dovs0AohRDbuqwtxOOQMd4YYhoHP58Pj8WDKzwPAHTbolDPcot/+OSIjx4rhSI4cGVVFTqqKnIfUNsdq5lNnjOe6U8fyt7WNrGvw0t4Xpd0foaMvSiCaSD46E8M+32Uz47Ca6PBH2dLcx5bmPh58dQ92i8opVYWcMbGI+WPzmV7uxmE9/AJcckSkIzki0pEcEeLwSMGdIZqmsWPHjuSIjwX7upTLNdxiwKAckZFjxTAkRzLHbjFx/WnjuP60cYOmB6IJOvoitPdFMQwDt8OCx2HBbbeQazen7gHe4Y/wxu4uXq9JPjr9UV7b1clru5LbAJOqMKkkl5mjPMyu9DBzlIfp5e73fRZcckSkIzki0pEcEeLwyLcmC5j6C25PWGVJ+WkYhiFHEIUQ4jiUazOTW5zL+OLcg7Yrcdm5bF4ll82rxDAMdrUHeL2mkzf3dLOpqZeuQIwdbX52tPl5Yl0TAGZV4azJxVw2fxTLp5VKF3QhhBDiOCAFdxYw5eWDxYIjP58fLLpDim0hhDiJKIrClDIXU8pcfOqM8RiGQVtfhM1NPrY0+9jUnPzZFYixekcHq3d04LKbuWhWOZfNG8WicQWoqmw3hBBCiGwkBXeGKIqCw+FAURRUZw5TN22UQlsMsn+OCDEcyZETk6IolHsclHscnD+jDEheQ7mnM8CT1c08Vd1Cc2+Yx99r5PH3GqnMT96PfEKxkzKPnXKPnTKPg1ybeUiOxBI6PcEY3cEoPcEYLruF2aM8UrCfxGQ9ItKRHBHi8Mgo5VkmrsXR0bGZbJkORQghRBbSdYN36np4srqJZze3EYiOPEBbmSd5mzJvKEZPIIZ/mLZlbjsXzCzjotnlLBiTf9DiO67pNHvDlHns0qVdCHFUZPu+uhDvlxTcGaLrOl1dXRQVFaGqKgC3v347z9Q+wx2n3cGVk6/MWGwiOwyXI0LsT3JEhGMaq7a380ZNF619Edp8YVp9EfyR4YtwSA7Ilp9jpcBpoaU3MqhgL3XbuHBmOR+aVU5RrpVd7QFq2v3sbPdT0x6gtitAXDPIz7Hw8aVV3HDaODw5lmOxqOIokfWISOdY50i27KsLcaRIl/IM0XWd2tpaCgoKUFWVzvvu56IX3qBjtk7nHBmpXAzNESEOJDkiHFYTl8yp4JI5FYOmB6IJ2nwRmnuCbNy2g0Wzp1PsdlDotOJxWFJnsSNxjddrunh2cyv/2dZOe1+UR96s55E360d8TbOq4A3F+cWqXTz0Wi3XnTqWT55eRbFr5J5Zmm7QHYxS5LRJ9/UsI+sRkY7kiBCHRwruLBGrr6eopovS0arcGkwIIcRhybWZmViSy7gCO7ZeGwvH5Q97Ox+7xcR500s5b3op0YTGmpountncyqpt7cQ1nUklLiaV5jK51MWU0uTvZW47z25p49cv72ZHm58HX93D79+o45pFo7nprAnYzWpqhPUdrX3sbPezq91PJK6Tl2NhwZh8Fo4rYNG4fGaO8kjXdCGEECc0KbizhKkgHwB3yKA5JAW3EEKIY8tmNnHutFLOnVaKrievNhvpbPQlcyq4eHY5q7d3cN/Lu9nQ2Msf3mrgj283cLAL1XpD8dRI6wBWk8rsSg/zx+YzodjJmAIn44pyKHXZ5Uy4EEKIE4IU3BmiKAoejyc14qO5/17c7hBskDPcgqE5IsSBJEdEOh80Rw6l2FUUheXTSzl3Wglv7enmvpd38+aebhQFxhbkMKXMxdQyN1PLXEwtd1PusbOjzc/a+h7W1ntZ29BDVyDG2gYvaxu8g+ZtM6uMLcxhbKGTCcW5zBuTx/wx+Qftti4+GFmPiHQkR4Q4PDJoWpbw/uWvtH3nO7w3SeGR68t46aqXMh2SEEII8b6090XItZlx2tIfzzcMg4buEO/V97C52Ud9d4i93UEavWE0ffhdkzEFOSwYm8/8sfnMH5PH5FIXFpNcUyrEiSRb99WF+KDkDHeG6LpOS0sLFRUVqKo6qEt5d6QbTdcwqXJd28nswBwR4kCSIyKdY50jpW77IbdVFIVxRU7GFTm5cuHo1PS4ptPSG6a+O0RDd5DtrX2sb+hlV4efvT0h9vaEeLK6OdXeZlZx2szkWE04rWacNhNOm5lJJS5Wzqtg1ig5M3cwsh4R6UiOCHF4pODOEF3XaWpqoqysDFVVMecnC+6SmJ2Lqi4kqkXJUXMyHKXIpANzRIgDSY6IdI7HHLGYVMYWOhlb6ASKU9P7InE27O1lXYOX9Xu9bNjbiz+aIJrQiSZi9AQHz+f1mi5+90YdE0tyuWzeKFbOG8WoPMexXZjjwPGYI+LYkhwR4vBIwZ0lTAUFKFYrZXmV3HXGXZkORwghhMgqbruFMycXc+bkZBGu6wa94TihWIJgVCMYSxDq/9kXjvPqrk5WbWtnd0eAu1/Yyd0v7OTU8QV8ZF4lc0bnkWs347KbybWah1yzHogmaOkN0+wN09QbpqU3THtfhIRmoBkGum6g6Qa6kfxpt5hYOrGI86aXvq+z/EIIIU58UnBnCWtVFVM2bpBub0IIIcQhUFWFAqeVAqd12L9fuXA0fZE4z29u4x/VTbxd25N6HCjXliy+HVYT3YEYvnD8fcfz3JY2vvXUFuZUevpvtVbG5NJc2a4LIcRJTgruDFFVleLi4lTXnP03yHEtTsJI4DBL17eT2YE5IsSBJEdEOid7jrjtFq5aNJqrFo2muTfMU9XNPLu5lTZfBH8kQUzTgeQZ7UA0Mei5HoeFijwHo/IcjMqzU+qxYzObMCnJYl9VFEyqgklR6AxE+c/2djY09rKxycfGJh8/e3EXYwpyOGNSEWVuO0UuG0W5NopyrRTl2ih22bLiHuQne46I9CRHhDg8Mkp5lvnJuz/hz9v/zGfnfpbPzPlMpsMRQgghTliRuIY/kiy2/ZE4wahGYa6Vco8dl93yvufX4Y+wensHq7a1s2Z3F7GEftD2VrOK1aRiNatYTEr/z+Q0gFhCJ6bpxDWduGYQSyR/n1Say/nTy1gxQ86iixNPtu+rC/F+ScGdIbquU1dXR1VVVeqIYdsP76LunVXcfWonc5ddxR2n3ZGx+ETmDZcjQuxPckSkIzmSOcFogtdrOtnU5KM7EKMrEO1/xOgMRNMW44dqbGEO508v5fwZZcwfk4/pEO6hvj/JEZHOsc6RbNlXF+JIkS7lGaLrOp2dnYwdOza18oru2EHerjaKp6k8X/88/pifOcVzmFM8h6kFU7GY3v/RdnH8Gi5HhNif5IhIR3Ikc5w2MxfMLOeCmeVD/mYYRv9Z9UTqrHV04Gx2/08FZdBZb1v/T4B363p4YWsbr+/uoqE7xMOv1/Hw63UUOq2cObmYJRMKOX1SEeWe9JemSY6IdCRHhDg8UnBnEVNBAQCFUSv+mJ/n65/n+frnAfj+ku9z2aTLMhmeEEIIIY4ARVFw2S0fqNs6wLgiJ1ctGk0wmuC1XZ28uK2d1dvb6Q7GeLK6OXWf8vHFTk6fWMSSCUWcNr4QT86hvZ6uG/gjCbyhGAndoMRtw2UzS9d1IYT4AKTgziKm/DwAbh53PZecv4SNHRvZ2LmRTV2bmFMyJ9Xu+frnebb2WT4969PMKp6VoWiFEEIIkUlOm5kLZ5Vz4axy4prOe/U9vLG7izd2d7OpqZfaziC1nUH++FYDkLxmPNdmxmkz4bSacdrM5FhVQv4+lOp38Ibi9Ibi9IbjaPrgKw5zrCbK3HbKPHbK3MlB5MYU5DCjws3kUldWDAAnhBDZSAruDFFVlcrKykFdc8z9Z7hVn59Ty0/l1PJTgWTXswGGYfC7zb9je892Xm58mcXli/n0rE9zStkpcuT5BDNcjgixP8kRkY7kyMnDYlJZMiF5NvurK8AXjvN2bTdv7u5ize4u9nQGiSV0ehIxeoLDzaF3yJQcqwmTouCPJgjFNGq7gtR2DX2yWVWYWJLLjAoPMyrczKhwU+5xYDYpmNXkaO5mk5r63RuK0eaL0N4XodUXoa0vQrsvQk8ozrzReVw+v5IxhTlH/D0SH4ysR4Q4PDJoWhbp+dOfaf/hD3GtWEHlL+8ZsV2tr5bfbf4dz9Q+Q8JI3sZkdvFsPjnzkywdtRSbyXaMIhZCCCHE8aAvEqcvHCcU0whEEwT7H4GoRlzT8Tgs5OVYyM+xkp9jJS/HkjprHYolaBsojPsitPmitPdF2NMZYEuzD2/o/d+3PJ1TxhVw+YJRfGhW+Qfuei+OT9m8ry7EByEFd4ZomsauXbuYPHkyJlNyg+Z75hlavvwVchYtYuyf/ph2Hi2BFh7Z+gj/qPkHUS0KwOWTLue7S74LgG7oxLQYdrP9qC2HOHqGyxEh9ic5ItKRHBHpHG6OGIZBqy/ClmYfW1v62NriY1tLH73hOAndIKHpHNA7HbOqUOKyUeqxU+6xU+pOdlPPsZp4sf+WagN7p3aLygUzylg5bxQTinMpzLWSY5UOmsfSsV6PZMu+uhBHiqyxMsQwDHw+36Du4ub8fBSrFQ5xZVaRW8E3Fn+Dm2bfxJ+2/Yl/7vkn80vnp/6+vWc71z97PbOKZjGjaAZmNflxK/3/Lqy6kCkFUwDY1LmJv+78K+FEmHAiTESLEI6HKc8tZ+XElSypWJJ6vjg2hssRIfYnOSLSkRwR6RxujiiKQkWeg4o8B+fPKBu2ja4baIaBphvENR2n1Yw6wu3Lrj9tHK2+ME9WN/P3dU3s6Qzy1IYWntrQkmpjt6gUOm0U5lopcFopyLGSa09ek57b/0j+bsJuMY14yZ3NrPaf0beQl2PFapYu08OR9YgQh0cqqCySc9ppTNm44X1fi13kKOJLC77ELfNvQTO01PRNnZuI63HWd6xnfcf6Ic+bWjA1VXC3Bdt4es/TQ9ps6d7CqoZV3H7K7Xxs2sfe5xIJIYQQ4mSnqgoqChYThzS4WrnHwWfPnshnzprAxiYff1/XxMs7O+j0R4kmdCJxnebeMM294SMap9NqIi/HSr5zX9f6Ameye33yZ7K4r8izM7ogJ3WbtqPBH4lT13+9/KxRHhmnR4jjmBTcWeRwV6aKomBW9n2k10y5hiUVS1jbtpY6Xx0GySOTAz/Husem2k4pmMIt82/BYXakHlaTlbXta3m+7nkuGHdBqu1bLW/hi/o4s/JMciwyqIkQQgghjjxFUZg7Oo+5o/OA5JnWYEyjJxCjOxilJxijOxDDG4oRjCbwp65N11K/R+LaiPMPx7XkqOyhGLoBwZhGMHZohbxJVRiV52BckZNxhTmMK3QyriiHEpc9debdZh754EI0oeENxukKROnwR6jtDLKnM0htZ4DariCd/miq7bjCHK5eNIYrFlRS7JJxeoQ43sg13Bmi6zpdXV0UFRVl/aiPmq5hUvdtNG58/kbWta8DoMJZQVVeFRM8ExjvGc/4vPHMLZ571I/EtgXbCMVDhBIhVEWlzFlGvi3/hDoCfDzliMgMyRGRjuSISEdyZN99x3tCyeK9NxTDG4zj7f+/NxTHG0z+3hOM0dgTJnyQQn6Ay2ZOFd9uhwVfOE5PMEZPIIY/mkj7/KJcK+GYRjCWfC2zqrB8WinXnDKaMyYVYxqhW/6RdqxzJFv21YU4UqTgzjLNX72NeGMjFXf/FOvo0ZkOZwjd0Hlg4wM8vftpWoOtQ/5e7CjmpateSv3/d1t+h27ojHWPZax7LBXOCnqjvbQGWzEpptQ155qucc0z1xDVouiGTkJPoBkauq6TMBLMLJrJ/efen5rv2X85m+5I96DXtplslDnLmF44nZ+e+dPU9Hdb38VislBoL6TAXoDT4syKwjyuxekIdzAqd1SmQxFCCCHEITIMgw5/lLquIA3dQeq6QtR3BWnoCdEVSJ55P/A+5sMxqwr5TitFuTbGFzmpKnIyvtjJ+OJcqoqceBwWgtEE/97Uwv+928iGxt7Uc0flOVg+rYRil438/uvY851WCp3Jnx6H5ah2eT+asn1fXYj3S7qUZ4imaWzZsoWZM2cOGvExvHEj8b17SXR0ZGXBrSoqN8+9mZvn3ow34qXWV5t89CZ/5tnyBrV/bPtjtIfah53XgtIFPHLBIwCYVBM94R46wh3Dtu0MdQ76v8fmIWEkyDHnENNidEe6iWpRGvoahsTwrTe+NejggFW1UugoJNeay+T8yfz4jB+n/nbv+nsJxAOYlMHdwBRFodhRzMdnfjw17eFND9Md6cYwDAwMbCYbp1Wcxillpxx0gLlQPMQTu57gD9v+QLGjmP+76P9SBwA2d25masFULCbLiDkC0NDXQE+kB6vJik21YTPZsJn7f/Y/jvZBhbgWp6GvgUA8gNPiJNeSS641F6fFiaoM3siH4iHaQm14I168ES+90V6qPFXMK5k3pO1wDMNgbftaphdOx2lxvq84G/oaqO2tZUHZAtzWE2vDfbAcGY5hGLzc+DJ/2PoHzKqZRWWLuHTCpZTnlh+DaDMrGA+yrXsb49zjKM4pznQ4x8z7zZF0EnoCwzCwmOQ2TZnQFmzjhfoXKHOWsWLciiMyzyOdIycDRVEodSdHVz91fOGQv+u6QV8kTnd/l/eeYJS+cAK3w5I6413otOK2W+iJdtMWbGNK/pRhv1dOm5mrF43h6kVj2NHWx+PvNvKP9U0094b5w1sNB43TaTXhcVjw5FjxOMzJ3x0WdAMicY1IXCea0FK/R+IaZ08p5psXTR80H8kRIQ6PFNwZYhgG4XB4yIiP5vz8ZMHd05OhyA5dvj2fBfYFLChdMOzfDcPgyslXUt9Xz96+vdT11eGP+bGZbJQ7y6lwVgxq/4tlvyCmxTCrZlRFxayYMakmVEXFYXYMavv0ysEDvMW0GO2hdtqCbSgMLjRHu0ZjUkz0RHoIJULE9FiyAA8yZL7/3PPPEQ8QTMybOKjg/ueef1LfVz+ozSNbHyHPlseKcSv45uJvDip6fVEfj+14jMe2P0ZvtLf/TYL2UDtlzjJ8UR83Pn8judZcLp1wKZeMvwRf0MfGzo3s6t3FBeMuIM+el1z+3U/z8OaHh40T4N5l97JszDIA1rWv45naZyh0JM/wW1QLJiX5vppUE4vLFqcKkLZgG7W+WiyqJfUwq2YiWoQ6Xx1nV56diuHBTQ/y0KaHhn19p8XJPcvu4dTyUwH4vx3/xz3r7xnSrsJZwUXjL+LqKVdT6iwd8vf2YDv/3PNPntz9JI3+Rr572ne5fPLlQLKAimvxVDwj2dK1hdtfvx2zYuaU8lM4d8y5LBu97LgvunxRH//13H8RCoe4iqu4csqVB30vtnVv42drf8Z7be+lpr3b9i5njDojVXDXeGsIxoMoikJPuIfuSDc9kR66w90E4gF+ePoPU88NxUM4zI6jcmCnN9JLfV89jf5GWoOttARaaA22UuQoGhTDX3f+lZKcEk6rOA2baeTrGu/fcD+/3/L71O0Tx7jGML90PgtKk+uvytzKo36AStM1OsOdNPobaQu2Mdo1mtnFsw/pgNPhGGlbcyjP6wx3ssu7ixpvTfLRW0Ntby13nn4nF1ZdCCS/h3t69zCjcMagS48G5rH/+/rztT+nyFHEyokr8dg8h79wJwlN11jTvIa/7fobrze/jm7oXDLhklTBndATfO+t7zElfwrTCqcxu3g2FvXQDog0+Zt4eNPDvNv4Lrfl3cayscsOO17DMOiJ9OCNeDEwmJQ/KfW3x7Y/Rme4k3HucZxVeVba9feRsH8edoQ6eK3pNbrCXXSFu1AVlbHusYxxjWGcexzlueVH5I4sqqqQl5McZG3CAZuarnAXa9vX8Nju93iv/T3qfHVAcru5pGIJd51+14i3c51a5ua7l8zg9gun8sLWNra3+vEGY3T3d3f3BmP0hGL09t8XPdjfHb3FF0kTsQFKAjCYVJo79K8fcD0ihEiSgjvLmAoKANC83gxHcvgUReG/5/x36v+GYRBOhEfcSZ9TPOcDv5bVZGW0azSjXUN7Bfx2xW9Tv4cTYXoiPfSEewjEA0MK7htn3Ig36k2dtR7YuBgYFNoHH8W+aspVeCPe1LJ2h7t5ufFleiI9NAeaU8uoGzo/W/sz/r7r74QSISB5EOATMz/BJRMuwWqyAlDfV4/b5qYr3MXvt/6e32/9PSoqeo0OJG8Dd2blmQDML53PmPoxRLUoMS1GRIukuuMDFDr2xbq1ayt/2/W3Ed+73yz/Tar4XNO8hu+99b2Dtl0yagkA4z3jcVqc5NvyCSVC+GN+4nr/Rj4exKpaU88rsBfgsrjIs+eRb8/HZXGxoXMDLcEWHt78MBdUXZAquGNajFcaX+HJ3U/yZsubqWVyWpz4Y/7UPP++6+/8cv0vWTFuBVdOuZJZRbN4q+Utnt7zNHOK53D99OsBOHfMuVTmVtIUaOLNljd5s+VN7nz7TmYXz+bcMedy1ZSrUmfNDywQDkVcj9Pga2C3bze7vbvZ27eXkpwSbpx5I0WOovc1r5H4oj7eaH6DnkgP102/DgC31U1cj9MWa+PeDffym82/4aLxF/HRqR9N3X0AwBvx8rO1P+Nfe/6FgYFVtXL99Ospd5azoXMDUwumptr+adufeHL3kyPGccdpd6QK2x+8/QPWNK9hZtFMZhXNYkLeBELxEL6oD5NqSr3/ALe9dhvBeJACewH59nwK7YXk2/PJt+XjtrkHffcvfvLiIQeyBux/+UVUi/LztT8nlAjhtDg5c9SZLB+7nCUVS9javZVphdNSPRo8Vg9RLUqRo4jucDd7/XvZ69/LU7ufAmDVFasocyZvZ/Sz937G5q7NRLVo6qGgUOYso9xZzh2n3ZH6zvqiPuJ6nL5YH/6Yn75o/89YH6qictWUq1LxXvvMtWzv2T5oeUpzSlkxbgUrxq1gVtGsYXMvkoiwtXsrGzo2sLNnJ5dMvITTR50OwC7vLu5dfy+6oaMoCk6zk1xrbqq3ySllpzCzYCYArcFWXtj7AoFYAH/Mjz/uxx/zE4qHiGpRrpl6DSsnrgRga/dWrn/2+tT3+UA13ppUwf1my5vc+sqtuK1uTi0/lYVlC+kMdbKpaxM7enaw6opVqfVsVIvys7U/477q+7ho/EVcM/WaQfm3v+5wN2+3vs2bLW/SHmpntGs049zjGOsey5ziOeTb8we1H7hEpyO076GgsHzs8tRnG9fiycFFj9EtLkPxELt7dzO7eHZq2qbOTdT31acODNvN9kHfC7fVncqD9mA7/9j9D/5R8w/agm2pecwtnss5Y85J/b/OV5fKZUiuby+deCmXT7p80OCo+2vyN/Hw5of55+5/kjCS1xKXO/f1dNnTu4eEnmBy/uQheXngevLR7Y9S56ujJdBCc6CZlkALES1Z4E3Jn8ITlzyRavvYjsdo6EuelTUpJuaVzOOcMeewbPQyKl2VQ+LUDZ32YDt7/XtpD7XjMDtwWV2McY2hInffQfuEnqAl0EJ9Xz31vnrq++pp6Gug3lfPTbNv4uqpVwPQ6G886Dbuptk38fl5nweS373qjmrmlswdsq9woLgWZ5d3F9t7thOKh9AMjYSeIGEkSOgJbpxxIy6rC4CvvfY13m17N/VcBYVcay7+mJ8ab82gYvs/Df/BZXURjAfxRX30xfrojfbii/qYWjCV2y+8akgsAAlNxx9J0BuO4zvg8UbHk9QEXydhhEkYEWJ6mJgeQkfnzNKL+dL85QddViHE+ycFd5YxFSR3IrSe47/gPpCiKBkf1dxhdjAqd9SI100PFDOHYv+CYsC39G+xvn19aocckjtNf9r2JwAm50/mU7M+xXljzxuy0zeneA4vXvEia5rW8I+af/Ba82vohk6BvYDphdOxm/ZthE8fdTrPfOSZIa+f0BNEEhFs5n1n++aUzOEzcz5Dd7gbb9RLXI+j6Rq6oaMZ2qAzDC6riyn5U4jr8X0PLY5JNTHOPW5QzCvGreBDVR8atOMV1aIEYgEC8QAlOSWp6SsnruSySZcNijWSiPBK4ytUd1QzOX9yavqXX/kyrzS9kvr/wtKFXDbpMpaPWT4ofzZ2biSmx/hX7b/4V+2/cJgdhBPJkWV39uzkumnXoSgKdrOd5y5/jnpfPav3rualvS+xqWsTGzs3srFzI8vHLE8V3PdW38uj2x/d1z3e7MRhcWA32XGYHXx10VdTO++PbHmEp/c8Tb2vPrXDur/9e0P8bdff2Nmzk0l5yTM9oURywL9QPERCT/D1xV9Ptf3Ful/wbuu7RLUokUSEiBbBG/GiGRpOi5Orp1yNxWRBURR+tPRHPFf9HOti69jh3cE/apI75/NL5vO5eZ9jUdkirCYra5rXYGDwoaoP8cX5X0ztqA7shA6wmWwUOYpSl10U2Av2/bQXDjq7saNnB73RXtY0r2FN85pB8ylxlAz6frQGWtnQuWHIewRQmVvJc5c/l/r/wCUhpTmljHGPodxZzqjcUZQ7ywftkIfjYVZOXMl/9v6HjlAHz9U/x3P1++bzvSXf4yOTPgLAh8d/mHml85heMB1/3M+Gjg2sa1/H+vb1eKPe1GcKUNNbM+xtFPf69+IwO/jB0h+kpn37jW/zcuPLwy5XjjmHKydfmfp+lDnLqPHWUJ5bTmlOKdt7ttMeaueP2/7Ifxr+w/OXPw+AP+bnndZ32NCxgerOarZ1byOh78uvuSVzoX/15Yv6eLXp1WFfH+BLC76UKrjbQ+38cv0vR2y7/2U3FtVCXI+nzv5NypvEpPzkY3L+5EHrT2/Ei8vioi/Wx4sNL/Jiw4uD5ru9e3tqrI5lo5exrn0du7y7+HvN3/l7zd+ZVzKPa6Zcw3ljz0t1p/37rr/z3be+O2g+77S+k/r9l8t+mSo4H93+KA9teoieyPC9wqYWTE19vn+v+Tt3vXMXHpsHl9WF0+Ikx5xDjiWHHHMO/z3nv1ProteaXuOp3U+lvoORRISEnqA0p5TR7tF8ZOJHmJg/MfU6hmHQ0NeQXLd0bGRT1yZqvDVohsYb176ROvjzrz3/4vGdj4/4OTx96dOMzxsPwM2rb2andyeQ/F5cOuFSrph8BeM84wY9x2V18dm5n2V793Y2dm6kJ9LD77f8nt9v+T2LyhbxhXlfSOYNyYLzfzf/76BC+7Ty05igT2Bi3r7l+c2m3/Bc3XNM8ExgQekCvFEv7cF22kJtFDuKefzD+5bh8R2PDzlIpqDgtrnJtQ4+W3rR+IvoDndT3VHNLu8u1ravZW37Wn763k9ZPmY5/2/Z/wOS6/DbXruNJn8TMT025H36xMxP8KUFXwKSBw8+/OSHB90WdX+d4X2XpFU4Kzir8iyKHEUUOgqJ63H29u2loa+BvX17Bx2g2Ni5kZtW3YRZNTOzcCaLyhaxsHQhc0vm0hfrozSnNPX9/vqar/NC/QvDvj7AFZOuSBXcp5SdQl+sj4WlC1lUtogFpQtwWV1s7dqKL+ZLPSeqRfnGmm+ktmsHGjhgDMleENc9ex3TC6ezuHwxi8oWETJCrO95m7db3+am2TdxVn9PA8+e8XxzzW+HnWdujsbEEteIyyGE+GCk4M4Qk8nE1KlTh1wLY87vL7i92d+lXAxlVpPdlve317+X88aex8qJKzlj1BkHPYNqUS0sG7OMZWOWJbv09nYzoXTCIY8KalbNQ3Zw5hTPOeTeAwNn2w71tQ5kM9mwOWyDzrDD8Le8s5vtXFB1ARdU7bvlXDgR5p22dyhxlHDpxEtZOXElY9xjhn39n531M7Z0beFvu/7Gc3XPEU6E8dg8XFR1ESsnrhzymuM84/jkrE/yyVmfpD3YzsuNL7N672rctn3XdftjfsKJMOFEeNBO2oAvzv9i6ve+WB+7e3cDybPvE/ImMClvEmPcY/DH/IPeg9UNq3mj5Y1hl0NB4WunfC3VtbjJ38TW7q1D2k3Mm8hZlWcR0SKpwmR2yWzGnjoWt9vNxq6NPLb9Mf7T8B/Wd6ynL9qXiu17S75Hob2QWcWzho1hwDdP/SbfPPWbB20z4C8f/gu7vLvY1LmJLV1b2OvfS641l3xb/pDu+rcuvJXa3lq8UW/qwE9PuAdv1ItZMQ+6E8LdZ92N2+pOe3Auz57H1xd/na+d8jU2d21mdcNqVjWsoinQhNvqJhQPDWo7cGDJbXVzZuWZqd4imj54J/1Tsz7F5ZMux262p8ZDiOtx2oJtqe72A7rD3amzU26rO/VwWV24rC6iWjR1turO0+8kx5yT+t5EtShrmtfwQv0LVLmrUvNt6GvgS698aVBMRY4i5pXMY2bRTBaWLUxNr/JU8b0l30NBwcAgGA8SiAcIxAIE40GmFkxNbWtCphCXTrg0FVuuJReX1UWOJQeH2UGVu2rffN1VvHj5ixQ4Cg7aVR+SPX0+MukjbO3eypstb7KxYyMlOSXMKp7FrKJZg4q40ypO44mLn2B9x3oe3/E4/2n4D9Ud1VR3VBNMBLly8pUATC1MnvWekj+FJaOWMN4zniZ/Ew19DTT0NVDl2RerqqipYtuiWijJKUk9gEFnQXsiPRgY9EZ7913Ws5+PTvto6vdGfyOrGlYNaTPQS+HMyjOZSHLZ/rXnX/zg7R8MWxiV5pTSGmjFXZBcz4zPG8/SiqWp23OG4qHU98Ef91NgL0g997JJl/Fi/YtcOeVKzht73oifRZmzjM/M+QyQ7HHzWuNrPFHzBG80v8F7be+lDtj4Y34u/+flqTiXVCzhM3M+w5ziOfh8vkH7I2bFjFW1sse3hz2+PYNeb/8DQACXTryUUDxERW5F6oB2mbNs0IHnAQNxQnJd90rjK7zU+BLr29cPKnZdVhe1vtpkLKqZytxKypxlxLQYfbG+QWfj+2J9aIaGzWRLDdA6zj2OKk8VY91jB+VLeW45951737Dv48BB6P3nW5pTSnuonQ2dG9jQuYGHNz+c+r49e9mzjHYne9XNLJzJWy1vMbNoJnm2PMyqOfnovzRu/zPkN82+aVDvvwEHrp+9ES8LShfQ5G/CbXPjsXrw2PofVs+grvo7vDvY0r2FLd1b+Ouuvw6Z98yiman2Z1WexS+X/RKnxZk86GTJwWlO/j7SmfyR9lmFEIdGRinPMt2//R0dd9+N++KLGXX3T9M/QYgTjC/qI9eSO+R60IPpi/VR56tjWsG0YXfyDlUgFsAb9abO0ofioVQBHk6EuXTipamz4Xt699ASaGFi3kTKnGUHPZDy0t6XqO6ops5Xh1k1Dzqr5rA4+PSsT6cKsc2dm/FGvdhMNuxmO3aTnXx7/qAeAwfTHmznX7X/4sYZNx6zrrPZwjAM2kPtFNoLj8mgXgOXnhzJ67Djepzrn72e6YXTmVcyj7klc4/JNeaZ0Bnq5ImaJ3hi5xMsGbUk1XtAN3R6Ij2HdEnGwLW4JTklaW8NmdAT9EZ76Yn0EIwHCcVDyZ/9vU3OH3d+6jV39uxkfcf6VO8Wu9mOqqi0Blpp9DcOumTkV9W/4qFND2Ez2ZheOJ05xXOYXTybWUWzBvWeSCemxbColtQyfJBLXPbXGmhlVcMqrp9+fWo+P3jrBzQFmvjMnM+kznqPxB/z89Lel2joa6DIUURpTimlztJBBzSOlN5IL5qhpQ5UarrGO23vMNY9lrKcsoNuD2JaDG/ES3FO8REfE8EwDJoCTaxtS56Jf6/tvdRdVn657JecNfqsVAz7f3bHWige4p3Wd3i37V3eaXuHGm8NZsXMrOJZnFp+KueOOXfQpUbZLtv31YV4v6TgzpBEIkF1dTXz5s3DbN63U9z71FO03fEdSr58KwU33JCx+ETmjZQjQgyQHBHpHA85Etfj9EZ6j9uBDL0RL13hLsa5x2X96O1xPT5kQLXjIUeySWeok1xrbtrrujPJF/VhUS1H7DK+Y50j2bKvLsSRImvWDNK0odcb5S5dyoTnnsUySu7NLIbPESH2Jzki0sn2HLGoluO22IbkHTsOHMQtW400enm250g2OR5y9WjcBUByRIgPTgruLGMuzv4VuRBCCCGEEEKI9I7uDUDFYQlVV9P34ovpGwohhBBCCCGEyDpyDXeGGIZBOBzG4Rj+ntSBNW/Q+KlPYcrLY8ILz2PyHPnuQSK7pcsRISRHRDqSIyIdyRGRzrHOkWzZVxfiSJEz3BlktY48mrLz1MVYJ05A6+2l68HfHMOoRDY5WI4IAZIjIj3JEZGO5IhIR3JEiA9OCu4M0TSNtWvXjjgIhWI2U3rbbQB4//xnYo2NxzI8kQXS5YgQkiMiHckRkY7kiEhHckSIwyMFdxZznnEGziVLMOJxOn7+i0yHI4QQQgghhBDifZCCO4spikLJ124DRcH//POE1q/PdEhCCCGEEEIIIQ6RFNxZzj5lCnlXXA5A+09+wkk2xp0QQgghhBBCHLdklPIMMQwDTdMwmUxpR3xMdHbScP1/UfjpT+G57DIUVY6TnAzeT46Ik5PkiEhHckSkIzki0jnWOZIt++pCHCnmTAdwMovFYjgcjrTtzMXFjH/2GSm0T0KHmiPi5CU5ItKRHBHpSI6IdCRHhPjgjssK7v7772fcuHHY7XYWL17Mu+++m+mQ3jdN09i0adMhj/i4f7Hd8f/uwfu3v6H5/UcrPJEF3m+OpNN25w/ZufhUGm/+HMF33pXLE04ARzpHxIlHckSkIzki0pEcEeLwHHdnuP/yl79w66238uCDD7J48WLuueceVqxYwc6dOykpKcl0eEddoquL7v/9X9A02u/8Ia5zz8Wz8lKcS5agmI+7j1McRXowiGKzpfLCXFSI7vMRWL2awOrV2KZMoeC//gv3hy9CtdkyHO3JxUgkTujvq2EYaL29JDo6SHR0oAdDmPLzsVSUYx09OtPhHfcMwyDe3EJk8ybCm7cQ2bSJeEcH7hUrKPnyrZkO76iLNzfjX/0S4Q0bMBIJchYuIGfxYmyTJh3RnmBGIkG4uhr/Sy8Tb26m8t5f7ouhowNzQcEJ/T0+kRmGQaK1FUPTsIwaJT0IhRBH1XG3pfjFL37Bpz/9aT7+8Y8D8OCDD/LMM8/wu9/9jttvvz3D0R19isVC8S1fxPfU08T27KHv2Wfpe/ZZTMVFuFdcQN6VV2CfMgUALRAg0dGJKT8Pk9MJB153YzKlNjKGrsPBjlyqKorJtK9tIoGh6xiJRPL31EPD5HFjcrkAiDU04Hv6nyS6ukh0daH19mLyeDCXlGAuKSZ36VIcc+cm5xuLoQUCyZhUFVQTiin5u6IoyXj3j+FQ443FiHd0YoRD6JEIeiiMEQmDakLNcWApL8dSUZGarxGPj7zx3X++hgGJROpPhq6jB4PogQB6MIiam5sqLrRAAO+fH0Xz+ZKPPh96b/J3FAXX8nMp/sIXUvMKb9mKUlyUfG/j8cFnoxUltZNnGAZ6X1/qvUdLkPB68T35FL4nn6T8rrtwrzgfgLyrr8Y6bhzBd97B99TTRHfupPWb36Tj5z+n4PrrKPrMZ4ZdrnTvgRGJoIfDqfd1/9/NZeXYp0wGQA+F8P3zX/2xxlHtdlSXC5M7mS/m0lIs5eWHFsMB78GhtgUw4vGj2tZIJIjW1hLdsZPIzh1Et+/ANm0qZd/4RqptzVlno9psWMeNxTJ2LNaxY7GUlaOHw1jKy3CeeioAejRK69e/jpHQMBcVYS4pxlxcnPz+FBejFBWl5qlHIvQ9+xyJjnYSHR3EOzognkBxOFAdDhyzZ5F/7bXJGGMxun/7W7Q+P3rAj+YPoPf1ofXnbs7ChZR/77upeYfWrsVUUIi5pBjd7yfR1U2iuwutuxtzWRm5S5cCoPX2UnfFlSQ6OjBisSFvWe655zL6/vuSMRgGdR+5HJPbjaWsDHNZGZayUsylZVjKy7CUl2PKy0suWyhE33PPpZYr0dmZ/J4qKigKOYsWUfiJj6dep+2Hd2FyuzEXFyXft6IizMXFmAoLUczmYXPHiMWS66mB1+joxDq6Etfy5an5dv/2d5g87n2fQUlJcl2nKEcszwzDQA8E0Hw+9L4+FJsN+9SpyffX76f5K18hsnkLWk/PkNnqweCgeXb+6j7spywCTRu6Htnve6wFAsR27ybh9aL1eNF6vWheLygqpoICnItPwT59enK+iQR6JIpqsw6/bAduK/ZfTxtGapqhaShWK6rdnnr/4x2doCUwND35s38bk+jowDpuHLaqKgAi27fTftddqdn6X3wRAFN+fioXUtsVTQNdH/mz2H87qGnowSDBN98k8PLLBF55NbmO7hdvbU2to1q+8lWiu3aRe9aZ5C47h5xFC9EDARIdHahOZ+r90gJBWr785RFf3jZlCiW3fin1/8Abb2DKzcVUWITqzEltDxVVBbN50MHRhNeLEYthRKMYsRh6/0/FbMHkdmEdOzbVVg8GwWIZvL2OJyARJ3FAbyfv439JHizz9qD1eDEikWQsOTlYK0dT9D//ve+9f+UV0DRUpxPVmZv8mevE5HSiWK0oFkuq7UG/EzCorR4KJXNW05Jxaonkd7SnByMaxXnaaam2rd/5LrHa2uT3cv/1SFkp5rIyzCUlqeuO277/AyLbthGtqUl9XxSHA9uECTjmzqXsW99MH29/Hg/kLkBkx45kzNEoHPB+qk4njjlzUv+P1taCYaDY7Kh2G4rdjmq1pvbRDnzPtECAWH09sYYGYg0NxBsaiNU3YK4oZ/R996Xaev/2NxTVhHVccptiaDrx5mbizc2oTieuc5al2u79xCcw4glM+fn9jzzM+fnYpkzFeerig35OQoj377gquGOxGOvWrePrX/96apqqqixfvpy33npr2OdEo1Gi0Wjq/319fQAkEgkS/Ttaqqqiqiq6rqPvt2EemK5p2qAdlZGmDwwmkThg5980sFOz346Hpmmoqpo8ynpAe7PZnBqgYoCiKMn5u1zkffzjeG68kejWbfj/9U/8zzyL1tmF989/xjZ7NuYJE1BVleBbb9H8+S8wktK77sJ9ycUAhF55hebPfX7EtsXf/jZ5V1+VbPv22zR/6tMjti25/XY8130MgEhLC12//vWIbVW7A9vs2ei6Tnj9eppu/PiIbYtu/RLFN92EruuENm2m8ZprRmxb8D//Q/EXPo+qqoR21bD3iitGbJv3X9dTfNttmEwm4m1t7Dnn3BHbeq69hpJv9m+Qe3upOf2MEdu6r7yC0u98BwA9GqPznntGbGufPSuVB3ooRH1/vB5g94ExrFxJ6Q/vTO6YB4PsWXzqiPP1r16N6/zzkrnkcpGzfDk5y5dT/IUv0PvE3/E++iiJtjbCNTUkEgnMZjPR3bupu/iSEeeZd+ONFH/lyyiKgtbaxp79CpIDeT56LSX9habu99P23e+O2NZ1ycWU9e9Ea3191C5ZOmJb90UXUfqTHwPJHZLd8+aP2Db33HOpvO9Xqe9Tzbz5IxboOaedytjf/z61LtizZCn6CJdu2OfMYez/PZZaF+w5dzlaR8ewbbVAgEQigclkSha23d1oQLylBd4cvO7KXbEC28KFAKgmE33PPjfisjmWLkW56dPJdVE8Tut+Rf2QGIJBXFdeCSQLoc5f3jtiW9ukSan1m6Hr7P34J0bc+XRdsILcpUtJJBIYDgfx1tZUkWXKz8dUXIzqdKL19mKuHJVaZ8Z8PqLbt48Yg/Pssxl1/32YTCa0aJTWb35rxLZKjoNEIoGqqhCL4f3Tn0Zsm7tiBeU//1nyecCuWbNHbnveeTjOPhtI5lnH3XeP2NZ1/vmU/eLnqf/XHGS+zjPOoPI3D6a2N7tPWYwRDg/bNmfJEsb+7rfJnLTbCVdvQO/rA7MZ+5Qp2GbOxD5rJqb+YkPXdVRVJbB+Pd0PPQQPPTTseqT4m98g79prMZlMhNato+m//2fEeItu+yrmyckDZ/GtW6m/euR1b9Gtt5Lff/AjumULe6+5dsS2hZ+7mYL/Sb5udPdu9n7k8pHn+7nPUfCZ/8EwDCwzZ5JzxhnkzJ+HYjITfPttwtXVaF4v/hdfJO+qq1Lb0J6HHqb73pFzvfKPf8R1yiIMw6Drj3+i6yc/GfR31ePBecYZOJedjZGbi6ZpKIkE0d270Xp78T39T3xP/3PQcwblmN1GYM2aEQ8QG/vtMxiGQdPNn8OIRIZt61y6lIrfPJj6/+6zlyULvGE45s9n7KN/Tq33ai+4EK2zc9i2tpkzMX3tNnRdJ5FI0Hn//SO3nTqVgps+nVrvtf/wLuKNjcO2tY4fz9h/Pp3aN6q97CPEdh+YiUnm8nKqVvUfODGZaLj+v4hs3Tp829ISqlavBpL7RtHduwmvWzdsWzU3lwlv71vHBt97j1hNTf+MzCiKghEOE9myBcVsHrRft2fpUvS+4df/joULGffnP6U+u72f+OSwB8Eg+Z5VPfkPILnv1/jZm4nX1w/b1jJuHFXP/Du5rdB16j5y+b54D3wfenpS6z1VVem6/9ck2tpGjNdxZnJ/RVVVQtUbhl3neC69FNvCBYOmmUwmDMNIbvP3y+Ph9m3h4PuwB+5njzRdP9hBMiGOQ8dVwd3V1YWmaZSWlg6aXlpayo4dO4Z9zo9+9CO+973vDZleXV2N0+kEoLi4mAkTJlBXV0fnfhuZyspKKisr2bVrF779jnKPHz+ekpIStmzZQni/FdbUqVPJy8ujurp60Ipm9uzZWK1W1q5dOyiGhQsXEovF2LRpU2qayWRi0aJF+Hy+QcvkcDiYM2cOXV1d1NbWpqZ7rriCqV/7GvVPPoXvtVepiUbQ166luLiYolgMnE7Y76zH/urq6oj3x1Q5wkZlQGNjI7v725p27SJ3mDaGyYRqsRBPxFPLqni9OM45h/Lp04k5c2j1B1CDAZTeXiz+APaZM1PLZNq+Y9j5DvD2eCnuj7tn+7aDtm1pbUHp6qKkpITdzU1YrFYMmxWsNqwuF2ank2DAD5EorZEoDWvXMnv27BHfqwEdHR3sXbsWk8nEvP6zLQdSc3IwHA46+/w09r8PntxcPJd/hJBqwm/oGE4nujMXz6gKyktKaQ0F971nXV14ioowuruHHCkfkMq9WAzPwESTCUNVwWQiMW0qsXOXM+aaq9E0bdjcy/nYR9k9fRrmdevxl5fRVl3NokWL8AcCB30P2tvaaFi7Fo/HwwTXvtFDDYsFw2ZDdTiwulzEVZW2hMbe/tcelZ9P7rnnEgiFiBk6SiyOEg5hS2io4TC9ikrTQJzB0L7lGoZu6PuWKZE4aFtg0PfJbRiMNMZqKBQCoKWlhaamJtyJxIhtA8EALS0tqXWEHoulBsVQXC4c06fjLywgWl6Bf+wYOtauZerUqXg8HoL3/QqjpQW1vR21vZ2CSAS9uwd/PE5Xbi7N/cu2cOFCCm7/Gs1NTSg+H2qvD9XnwxmLEW9vx2dKHrSrrq7G4XBQsPxc4lYrPpMZIy8Pw2ImR1Updnvwu3IH5UH+RReRV1KCLxEnoBsYOQ6MnByKKispmjhx33ovGCS3pARTTw+EwxgmE4bbjeF2o3vcMH4CQGq9Z/rmN9BdbmaecTq23MGv2QUUa1pyvbdlC6bbvora68Pk66XMZCbU1ERw717Unh56VBXvli3MmTOHnliM+Ny5GC4Xel4e9vIyyseOpaenB29PD4HiYpr613tVZWWoH72WcEsrqs+H4vNhCQahpwcjHqfH25N6f6vGjRv0mRpWK3p+Po6KCuzl5XQWFqTaEg4z6rKVyc+poR7F24vi96P0f0cP/J6ly8n91+VuXR+SZ4bdjuF0Euk/2zWwfTLfcANGnoey005j9PjxbN++nb0D2yevl/H96726tjaMpUswb9mKut/2a8Dehr3s7l/vKfn56IWF6G4XRq4Lw+WiZOIEtHic7rp66g2DPf3rvWn935GRNDY1sWdgvTfCAagB3u4eavvbqi2tuGw2VLMZHdBVBZTk+szi8aDmOgdvcz/5idQ2t27eXLRoFFNdHebtO2Da1NTnYWtuxj5yCOyq2cWCU5Lb3MbGvTgArbQUFi1k/NVX4x81irq9e5ONt27F4/Ewbdo0HP/3GG0vv4ylegPmDdWY2jtQHA6M/Dy6E/FU3lRWVlLx4x/RUt9AKLzvvSsuLsbtdrM3EqFlIG9icYrGV2H09hLv7EI54CCXoTAox9yqmuyxYbWimVSwWDBMZhQ92SNm//WeKxQaMmjPwPYimkiwaNEimpqaaGpqwj5/HsTiOEpLKZ4wga6An0BXF0okSjjXOWi9FysvR7VaUSIRLPE4RCJo/d+LcCTC2v71Xl5eHuFwGNMIn0MsFkst28KFC9GVoXGaLBaU/Hwiubmsfe89UBQcDgcTb/4sPfX1dNbUoPZ4UXp6sPb1Ye7rQ3c4Br1nBR+5jFGlpbTbbHTbbKAoqB0dFIcjFJaVDdrfcye0Edf/iVByX2FgvefMy0OxmLG7XCiqadC+od/tQhtY723ahFNRUJ1OlFhsyGcciUTY0r/e6+rqSr1nekE+asUo8qZPJ5yXR4/dhl5WRnf/em/8+PHoCxcQr6vD1NaO0tODoihYysuJ5+fTW1xM28B6r6qK0fffR+2GjcR7elACfhR/AI+q4pg3b8R92IFtzYAjsg/b/30a2OYOsNsP9q0V4vhzXN0WrKWlhVGjRvHmm29y2n7diW677TZeffVV3nnnnSHPGe4M9+jRo+nu7k7daiATZ7gNwyAQCODxeIYcyTsSRwf3XyYtFkPv30naP3bFbk91XVI0DaLREZdJM5uTXZ5InulRYjEwmdAVJdl9sL/d4R7xNHQdxTBQFQU9kUDv77qOrmNyODA7HAddpoHYFasVc07O+/6cDF0n7vcP6oK4/zIpFgtq/yidJlVF8/v3LZOqYsrJwWy1fqDP6cDpxOP0tLWRm5ub6g6nqiommw2jf+NnGAZoGiarFVVVDyn3DjbdbDajJxLE+3uCDBf7QDdQRVFQFQUtGMSw2VLdMo/E98kwDJRQaMTPSTGbMfq7VhqGge73j7hMFrsdxeFITdf6+kb+PMwWLG5Xarrm94NhDLtMismE2eVKxZ7o60seIFEUzG538szsYawj0n1OhmEQj8fx+/243e5kbnzAdcShfE6KoqBEo+jWwV2Jj/QyDUw3+t/3I7VMiqIk81pRUHNyUvMhGEy+psmE4nDs6010CMtkxOPo4fCQnIRkno203jNbrdC/Lhtou/8yYbenupwf7vdJ13V8LS3k5eUNWkcMdGP9IJ+TFgyS2G+7un/shsWyr9uzpiVz5oDvjWqxYDKb0VWV/XdAjsY2V49GU2eBh1smNScHi92e/D6FwxixGKb+de6h5p5hGCjxOCaHI/WeH+4yxePx5PpE1zF0HVP/ZQvG/pciaFoyl0b4nPbPPT0SgXgcs9WKYTKhD3RT7xcIBHC5XMPG+H5zL5FIJLtX63qye3z/MkV7e0fctqIoqUvRTCYTeiyGbhipfYvhlmmkz2P/6SPFeCjLpPn9qIoy7LKaHQ5MDsdhr/dMJhN6JEJioEBXVcwuV2qZ4n19yf0Om+19rff0WAyT2YzJaj0i3yfDMPB6vbjd7tTncTTPcAcCAfLz8+W2YOKEcVwV3LFYjJycHJ544glWrlyZmn7DDTfQ29vL008/nXYe2XJvv0Qiwdq1a1m4cCFmGXRFDENyRKQjOSLSkRwR6UiOiHSOdY5ky766EEfKcTUso9VqZcGCBazuv3YHktd5rF69etAZbyGEEEIIIYQQItOOu0OZt956KzfccAMLFy7klFNO4Z577iEYDKZGLRdCCCGEEEIIIbLBcVdwX3311XR2dnLHHXfQ1tbG3Llzef7554cMpJbtlP7BPgauhRHiQJIjIh3JEZGO5IhIR3JEpCM5IsThOa6u4T4S5LoQIYQQQgghspPsq4sTzXF1DfeJRNd1Ojo65F6DYkSSIyIdyRGRjuSISEdyRKQjOSLE4ZGCO0N0Xae2tlZWXmJEkiMiHckRkY7kiEhHckSkIzkixOGRglsIIYQQQgghhDgKpOAWQgghhBBCCCGOAim4M0RRFDwej4z4KEYkOSLSkRwR6UiOiHQkR0Q6kiNCHB4ZpVwIIYQQQgiRFWRfXZxo5Ax3hui6TlNTkwxAIUYkOSLSkRwR6UiOiHQkR0Q6kiNCHB4puDNEVl4iHckRkY7kiEhHckSkIzki0pEcEeLwSMEthBBCCCGEEEIcBVJwCyGEEEIIIYQQR4EU3BmiqirFxcWoqnwEYniSIyIdyRGRjuSISEdyRKQjOSLE4ZFRyoUQQgghhBBZQfbVxYlGDlVliK7r7NmzRwagECOSHBHpSI6IdCRHRDqSIyIdyREhDo8U3Bmi6zqdnZ2y8hIjkhwR6UiOiHQkR0Q6kiMiHckRIQ6PFNxCCCGEEEIIIcRRYM50AMfawCXrfX19GY0jkUgQDAbp6+vDbD7pPgZxCCRHRDqSIyIdyRGRjuSISOdY58jAPvpJNsyUOIGddGtWv98PwOjRozMciRBCCCGEEGI4fr8fj8eT6TCEOGwn3Sjluq7T0tKCy+VCUZSMxdHX18fo0aNpbGyUERjFsCRHRDqSIyIdyRGRjuSISOdY54hhGPj9fioqKuRWZOKEcNKd4VZVlcrKykyHkeJ2u2UDJw5KckSkIzki0pEcEelIjoh0jmWOyJltcSKRw0ZCCCGEEEIIIcRRIAW3EEIIIYQQQghxFEjBnSE2m43vfOc72Gy2TIcispTkiEhHckSkIzki0pEcEelIjghxeE66QdOEEEIIIYQQQohjQc5wCyGEEEIIIYQQR4EU3EIIIYQQQgghxFEgBbcQQgghhBBCCHEUSMGdAffffz/jxo3DbrezePFi3n333UyHJLLIj370IxYtWoTL5aKkpISVK1eyc+fOTIclstSPf/xjFEXhlltuyXQoIss0Nzdz3XXXUVhYiMPhYNasWaxduzbTYYksoWka3/72t6mqqsLhcDBhwgR+8IMfIEP7nLxee+01Lr74YioqKlAUhaeeemrQ3w3D4I477qC8vByHw8Hy5cupqanJTLBCHEek4D7G/vKXv3Drrbfyne98h/Xr1zNnzhxWrFhBR0dHpkMTWeLVV1/l5ptv5u2332bVqlXE43HOP/98gsFgpkMTWea9997jN7/5DbNnz850KCLLeL1eli5disVi4bnnnmPbtm38/Oc/Jz8/P9OhiSzxk5/8hAceeID77ruP7du385Of/ISf/vSn/OpXv8p0aCJDgsEgc+bM4f777x/27z/96U+59957efDBB3nnnXdwOp2sWLGCSCRyjCMV4vgio5QfY4sXL2bRokXcd999AOi6zujRo/n85z/P7bffnuHoRDbq7OykpKSEV199lTPPPDPT4YgsEQgEmD9/Pr/+9a+58847mTt3Lvfcc0+mwxJZ4vbbb+eNN97g9ddfz3QoIkt9+MMfprS0lN/+9repaZdffjkOh4M///nPGYxMZANFUXjyySdZuXIlkDy7XVFRwZe//GW+8pWvAODz+SgtLeWRRx7hmmuuyWC0QmQ3OcN9DMViMdatW8fy5ctT01RVZfny5bz11lsZjExkM5/PB0BBQUGGIxHZ5Oabb+aiiy4atD4RYsA///lPFi5cyJVXXklJSQnz5s3j4YcfznRYIossWbKE1atXs2vXLgA2btzImjVruPDCCzMcmchGdXV1tLW1DdrmeDweFi9eLPuwQqRhznQAJ5Ouri40TaO0tHTQ9NLSUnbs2JGhqEQ203WdW265haVLlzJz5sxMhyOyxOOPP8769et57733Mh2KyFK1tbU88MAD3HrrrXzjG9/gvffe4wtf+AJWq5Ubbrgh0+GJLHD77bfT19fH1KlTMZlMaJrGD3/4Qz72sY9lOjSRhdra2gCG3Ycd+JsQYnhScAuRxW6++Wa2bNnCmjVrMh2KyBKNjY188YtfZNWqVdjt9kyHI7KUrussXLiQu+66C4B58+axZcsWHnzwQSm4BQB//etfefTRR3nssceYMWMGGzZs4JZbbqGiokJyRAghjiDpUn4MFRUVYTKZaG9vHzS9vb2dsrKyDEUlstXnPvc5/v3vf/Pyyy9TWVmZ6XBElli3bh0dHR3Mnz8fs9mM2Wzm1Vdf5d5778VsNqNpWqZDFFmgvLyc6dOnD5o2bdo09u7dm6GIRLb56le/yu23384111zDrFmzuP766/nSl77Ej370o0yHJrLQwH6q7MMK8f5JwX0MWa1WFixYwOrVq1PTdF1n9erVnHbaaRmMTGQTwzD43Oc+x5NPPslLL71EVVVVpkMSWeTcc89l8+bNbNiwIfVYuHAhH/vYx9iwYQMmkynTIYossHTp0iG3E9y1axdjx47NUEQi24RCIVR18G6gyWRC1/UMRSSyWVVVFWVlZYP2Yfv6+njnnXdkH1aINKRL+TF26623csMNN7Bw4UJOOeUU7rnnHoLBIB//+MczHZrIEjfffDOPPfYYTz/9NC6XK3VtlMfjweFwZDg6kWkul2vI9fxOp5PCwkK5zl+kfOlLX2LJkiXcddddXHXVVbz77rs89NBDPPTQQ5kOTWSJiy++mB/+8IeMGTOGGTNmUF1dzS9+8Qs+8YlPZDo0kSGBQIDdu3en/l9XV8eGDRsoKChgzJgx3HLLLdx5551MmjSJqqoqvv3tb1NRUZEayVwIMTy5LVgG3Hfffdx99920tbUxd+5c7r33XhYvXpzpsESWUBRl2Om///3vufHGG49tMOK4cPbZZ8ttwcQQ//73v/n6179OTU0NVVVV3HrrrXz605/OdFgiS/j9fr797W/z5JNP0tHRQUVFBddeey133HEHVqs10+GJDHjllVdYtmzZkOk33HADjzzyCIZh8J3vfIeHHnqI3t5eTj/9dH79618zefLkDEQrxPFDCm4hhBBCCCGEEOIokGu4hRBCCCGEEEKIo0AKbiGEEEIIIYQQ4iiQglsIIYQQQgghhDgKpOAWQgghhBBCCCGOAim4hRBCCCGEEEKIo0AKbiGEEEIIIYQQ4iiQglsIIYQQQgghhDgKpOAWQgghhBBCCCGOAim4hRBCZK0bb7yRcePGZToMIYQQQogPxJzpAIQQQpxcFEU5pHYvv/zyUY5ECCGEEOLoUgzDMDIdhBBCiJPHn//850H//+Mf/8iqVav405/+NGj6eeedR0FBAbquY7PZjmWIQgghhBBHhBTcQgghMupzn/sc999/P7I5EkIIIcSJRq7hFkIIkbUOvIa7vr4eRVH42c9+xv3338/48ePJycnh/PPPp7GxEcMw+MEPfkBlZSUOh4NLL72Unp6eIfN97rnnOOOMM3A6nbhcLi666CK2bt16DJdMCCGEECcDuYZbCCHEcefRRx8lFovx+c9/np6eHn76059y1VVXcc455/DKK6/wta99jd27d/OrX/2Kr3zlK/zud79LPfdPf/oTN9xwAytWrOAnP/kJoVCIBx54gNNPP53q6moZpE0IIYQQR4wU3EIIIY47zc3N1NTU4PF4ANA0jR/96EeEw2HWrl2L2ZzcvHV2dvLoo4/ywAMPYLPZCAQCfOELX+BTn/oUDz30UGp+N9xwA1OmTOGuu+4aNF0IIYQQ4nBIl3IhhBDHnSuvvDJVbAMsXrwYgOuuuy5VbA9Mj8ViNDc3A7Bq1Sp6e3u59tpr6erqSj1MJhOLFy+WkdGFEEIIcUTJGW4hhBDHnTFjxgz6/0DxPXr06GGne71eAGpqagA455xzhp2v2+0+onEKIYQQ4uQmBbcQQojjjslkel/TB0ZA13UdSF7HXVZWNqTd/mfHhRBCCCEOl+xZCCGEOGlMmDABgJKSEpYvX57haIQQQghxopNruIUQQpw0VqxYgdvt5q677iIejw/5e2dnZwaiEkIIIcSJSs5wCyGEOGm43W4eeOABrr/+eubPn88111xDcXExe/fu5ZlnnmHp0qXcd999mQ5TCCGEECcIKbiFEEKcVD760Y9SUVHBj3/8Y+6++26i0SijRo3ijDPO4OMf/3imwxNCCCHECUQxBkaSEUIIIYQQQgghxBEj13ALIYQQQgghhBBHgRTcQgghhBBCCCHEUSAFtxBCCCGEEEIIcRRIwS2EEEIIIYQQQhwFUnALIYQQQgghhBBHgRTcQgghhBBCCCHEUSAFtxBCCCGEEEIIcRRIwS2EEEIIIYQQQhwFUnALIYQQQgghhBBHgRTcQgghhBBCCCHEUSAFtxBCCCGEEEIIcRRIwS2EEEIIIYQQQhwFUnALIYQQQgghhBBHwf8Huh5cGZQB9t0AAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# initial conditions\n", "init_species = 10 * np.ones(num_species+num_resources)\n", "\n", "# instantiate simulator\n", "\n", "simulator = sim_CRM()\n", "\n", "simulator.set_parameters(num_species = params['num_species'],\n", " num_resources = params['num_resources'],\n", " tau = params['tau'],\n", " w = params['w'],\n", " c = params['c'],\n", " m = params['m'],\n", " r = params['r'],\n", " K = params['K'])\n", "\n", "simulator.print_parameters()\n", "\n", "observed_species, observed_resources = simulator.simulate(times, init_species)\n", "observed_data = np.hstack((observed_species, observed_resources))\n", "\n", "# add Guassian noise to the data\n", "observed_species = observed_species + np.random.normal(loc=0, scale=0.1, size=observed_species.shape)\n", "observed_resources = observed_resources + np.random.normal(loc=0, scale=0.1, size=observed_resources.shape)\n", "\n", "# Replace negative values with 0\n", "observed_species = np.maximum(observed_species, 0.0001)\n", "observed_resources = np.maximum(observed_resources, 0.0001)\n", " \n", "# plot simulated species and resouce dynamics\n", "plot_CRM(observed_species, observed_resources, times)\n", "\n", "\n", "# Write the data to a csv file\n", "df_species = pd.DataFrame(observed_species, columns=[f'species_{i + 1}' for i in range(observed_species.shape[1])])\n", "df_resources = pd.DataFrame(observed_resources, columns=[f'resource_{i + 1}' for i in range(observed_resources.shape[1])])\n", "\n", "# Add the time column to ensure combination matches up\n", "df_species['time'] = times\n", "df_resources['time'] = times\n", "\n", "# Combine the species and resources into one DataFrame\n", "df_combined = pd.concat([df_species.drop(columns=['time']), df_resources], axis=1)\n", "\n", "# Ensure time is the first column\n", "cols = df_combined.columns.tolist()\n", "cols = cols[-1:] + cols[:-1]\n", "df_combined = df_combined[cols]\n", "\n", "# Write the combined data to a CSV file\n", "df_combined.to_csv('data-s2-r2.csv', index=False)\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "id": "ce932174", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "MIMIC_env", "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.10.12" } }, "nbformat": 4, "nbformat_minor": 5 }