{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Atom-wise distances between matching AtomGroups\n", "\n", "Here we compare the distances between alpha-carbons of the enzyme adenylate kinase in its open and closed conformations. ``distances.dist`` can be used to calculate distances between atom groups with the *same number of atoms* within them.\n", "\n", "**Last updated:** December 2022 with MDAnalysis 2.4.0-dev0\n", "\n", "**Minimum version of MDAnalysis:** 0.19.0\n", "\n", "**Packages required:**\n", " \n", "* MDAnalysis (Michaud-Agrawal *et al.*, 2011, Gowers *et al.*, 2016)\n", "* MDAnalysisTests\n", " \n", "**Optional packages for visualisation:**\n", "\n", "* [matplotlib](https://matplotlib.org)" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "execution": { "iopub.execute_input": "2021-05-19T05:57:15.912303Z", "iopub.status.busy": "2021-05-19T05:57:15.911769Z", "iopub.status.idle": "2021-05-19T05:57:16.659643Z", "shell.execute_reply": "2021-05-19T05:57:16.659995Z" } }, "outputs": [], "source": [ "import MDAnalysis as mda\n", "from MDAnalysis.tests.datafiles import PDB_small, PDB_closed\n", "from MDAnalysis.analysis import distances\n", "\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline\n", "\n", "import warnings\n", "# suppress some MDAnalysis warnings when writing PDB files\n", "warnings.filterwarnings('ignore')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Loading files" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The test files we will be working with here feature adenylate kinase (AdK), a phosophotransferase enzyme. (Beckstein *et al.*, 2009) AdK has three domains: \n", "\n", " * CORE\n", " * LID: an ATP-binding domain (residues 122-159)\n", " * NMP: an AMP-binding domain (residues 30-59)\n", " \n", "The LID and NMP domains move around the stable CORE as the enzyme transitions between the opened and closed conformations." ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "execution": { "iopub.execute_input": "2021-05-19T05:57:16.665302Z", "iopub.status.busy": "2021-05-19T05:57:16.664753Z", "iopub.status.idle": "2021-05-19T05:57:16.889002Z", "shell.execute_reply": "2021-05-19T05:57:16.889497Z" } }, "outputs": [], "source": [ "u1 = mda.Universe(PDB_small) # open AdK\n", "u2 = mda.Universe(PDB_closed) # closed AdK" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Calculating the distance between CA atoms" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We select the atoms named 'CA' of each Universe." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "execution": { "iopub.execute_input": "2021-05-19T05:57:16.894787Z", "iopub.status.busy": "2021-05-19T05:57:16.893567Z", "iopub.status.idle": "2021-05-19T05:57:16.896802Z", "shell.execute_reply": "2021-05-19T05:57:16.897195Z" } }, "outputs": [], "source": [ "ca1 = u1.select_atoms('name CA')\n", "ca2 = u2.select_atoms('name CA')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "`distances.dist`([API docs](https://docs.mdanalysis.org/stable/documentation_pages/analysis/distances.html#MDAnalysis.analysis.distances.dist)) returns the residue numbers of both selections given. The `offset` keyword adds an offset to these residue numbers to help with comparison to each other and other file formats. Here we are happy with our residue numbers, so we use the default offset of 0. (See the documentation of `distances.dist` for more information.)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "execution": { "iopub.execute_input": "2021-05-19T05:57:16.902403Z", "iopub.status.busy": "2021-05-19T05:57:16.901729Z", "iopub.status.idle": "2021-05-19T05:57:16.903560Z", "shell.execute_reply": "2021-05-19T05:57:16.903935Z" } }, "outputs": [], "source": [ "resids1, resids2, dist = distances.dist(ca1, ca2, \n", " offset=0) # for residue numbers" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Plotting\n", "\n", "Below, we plot the distance over the residue numbers and highlight the LID and NMP domains of the protein. The LID domain in particular moves a significant distance between its opened and closed conformations." ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "execution": { "iopub.execute_input": "2021-05-19T05:57:16.922914Z", "iopub.status.busy": "2021-05-19T05:57:16.917665Z", "iopub.status.idle": "2021-05-19T05:57:17.052135Z", "shell.execute_reply": "2021-05-19T05:57:17.052577Z" }, "scrolled": true }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjMAAAGdCAYAAADnrPLBAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAACAlElEQVR4nO3dd3ib5dU/8O+j7SHLe8UjzoSQELIZISQ0CYRRAh2UUgqlpS2rZfb3UlpG+5YU3pbSt9DNbEuhfdkFAgGSEEhC9iDTSezYie14S7Jl7ef3x6PnkWRLtmRr2t/PdfnClmX7DlpH5z73OYIoiiKIiIiI0pQq2QsgIiIiGgkGM0RERJTWGMwQERFRWmMwQ0RERGmNwQwRERGlNQYzRERElNYYzBAREVFaYzBDREREaU2T7AXEm9frRVNTE4xGIwRBSPZyiIiIKAKiKMJqtaK8vBwq1eC5l1EfzDQ1NaGysjLZyyAiIqJhaGxsREVFxaDXGfXBjNFoBCD9z8jJyUnyatKfzWXDpw2fQqfSQa/RJ3s5CeFwO+D0OnFe1XnI1GYmezmUKtw2oO1TQKUD1GPjsRATHgfgdQJF5wEaPp4oPIvFgsrKSuV1fDCjPpiRt5ZycnIYzMSAxqVBljELRp0RBo0h2ctJCLvbDqvTipycHAYz5OfWAPYsQGsE1GPjsRATHjvgsgI5OQxmKCKRlIiwAJiIiIjSGoMZIiIiSmsMZoiIiCitjfqaGSIiongSRRFutxsejyfZS0krarUaGo0mJm1TGMwQERENk9PpRHNzM2w2W7KXkpYyMzNRVlYGnU43ot/DYIaIiGgYvF4v6urqoFarUV5eDp1Ox+asERJFEU6nE21tbairq8PkyZOHbIw3GAYzREREw+B0OuH1elFZWYnMTB4zj1ZGRga0Wi2OHz8Op9MJg2H4LQ5YAExERDQCI8kojHWx+n/HW4CIiIjSGreZiIiIYs3jBER34v6eoAHUIyuiTWcMZoiIiGLJ4wQ6tgDunsT9TU02UDB/zAY03GYiIiKKJdEtBTIqnTS7K94fKp3096LIBN1www1YuXJlyO+NHz8eTzzxRNDXgiBAEARkZGRg/Pjx+OpXv4qPPvpohP+jYofBDBERReSzRife2N+X7GWkD7VeGkIa94/4T23/2c9+hubmZhw6dAgvvPACcnNzsXTpUvziF7+I+9+OBLeZiIhoSC6PiO+82gWrQ8T0Ei0mFvDlYywxGo0oLS0FAFRVVWHRokUoKyvDAw88gC9/+cuYOnVqUtfHzAwREQ1pT4sLVocIADjYlsDCVkpZP/zhDyGKIt54441kL4XBDBERDW1zg1P5/EgHgxkC8vPzUVxcjPr6+mQvhcEMERENbRODGQpBFMWUGOHAYIaIiAblcIvYdjIgmOlkMENAR0cH2traUFNTk+ylMJghIqLB7W52we4GdGrp62Odbni8YnIXRUn329/+FiqVKuwR70RiOToREQ1K3mK6cKIeHx11wOEGTlo8qMrlS8igPI6U/jtmsxm7du0Kuiw/Pz/kda1WK1paWuByuVBXV4e///3v+Otf/4pVq1Zh0qRJw/r7scR7IhERheXxivj0uPRieV61HnWdHhxqd+NoB4OZsASN1JHX3QN4nUNfPxY02dLfjcK6deswa9asoMuuv/76kNd94IEH8MADD0Cn06G0tBRnn302PvzwQyxZsmTYS44l3hOJiGgAt1fELW90Y0O9E30uaUvpnCodNjc4cajdjSMdbiyZGP9mbWlJrZNGC6TwbKbnnnsOzz33XETXTYXTSkNhMENERAMcbnfj/VopI6PXAJdMMWBivhoTC6TCGZ5oGoJaB2BszklKBgYzREQ0gMUuZWPG56rx4XcKoVZJx28n+Tr/8kQTpRKeZiIiogEsDi8AIDdDpQQyQEAw0+GGKPJEE6UGBjNERDSAxTe6IMcQ3BBtQp4GAgCzXUS7zZuElRENxGCGiIgGsNilQCVHH/wyYdAKqMyV6maOdXoSvi6iUBjMEBHRAEpmRj+wVX1+hvTSYbYzM0OpgcEMERENoGRmDANfJrJ0UoAjH9kmSjYGM0RENMBgmZkMrXRZr5PBDKUGHs0mIqIB5NNM/WtmACDLF8zYmJkJy+lxwu1N3PF1jUoDXRRN80YbBjNERDSA3Gem/2kmAMjUcptpME6PE1tObEGPqydhfzNbm435FfMjDmhuuOEGPP/881i1ahX+67/+S7n89ddfx5VXXglRFLFu3TosWbIEubm5aG5uhsFgUK63ZcsWLFiwAACUI/ry9WWFhYWYO3cufvnLX2LmzJmx+GeGxW0mIiIaYLDMjLLNxGAmJLfXjR5XD3QqHYw6Y9w/dCodelw9UWeCDAYDHn30UXR1dQ16PaPRiNdeey3osmeeeQZVVVUhr3/o0CE0Nzfj7bffRldXFy6++GKYzeao1hYtBjNERDRAuD4zAAuAI6XX6GHQGOL+odcMb0bW0qVLUVpailWrVg16veuvvx7PPPOM8nVfXx9eeumlsEMpi4uLUVpaivnz5+PXv/41WlpasHnz5mGtMVIMZoiIaIBwfWYAf2aGNTPpTa1W45FHHsHvfvc7nDhxIuz1rrvuOmzYsAENDQ0AgFdeeQXjx4/H7Nmzh/wbGRkZAACXyxWbRYfBYIaIiIJ4RRHWQU4zyTUzNp5mSntXXnklzjrrLDz44INhr1NcXIwVK1YoU7afeeYZ3HjjjUP+7o6ODjz88MMwGo2YP39+rJYcEoMZIiIK0uMUIYcpoTIzmTpmZkaTRx99FM8//zz2798f9jo33ngjnnvuORw7dgybNm3CtddeG/a6FRUVyM7ORmFhIQ4cOIB///vfKC4ujsfSFQxmiIgoiJyV0akBfYgzr0pmxsUOwKPBokWLcNFFF+HHP/5x2OtccsklsNvt+Pa3v43LL78cBQUFYa+7YcMG7N69G2azGYcPH8ZFF10Uj2UHSWows2rVKsybNw9GoxHFxcVYuXIlDh06FHSdG264AYIgBH2cffbZSVoxEdHoF1gvIwjcZhoLfvnLX+Ktt97Cxo0bQ35frVbjuuuuw7p164bcYqqpqcHEiRORk5MTj6WGlNRgZv369bj11luxefNmrFmzBm63G8uXL0dvb2/Q9S6++GI0NzcrH++8806SVkxENDztvR68tNsGmzP1sxmDdf8FgEyt9NLBbabRY8aMGbj22mvxu9/9Lux1fv7zn6OtrS0hmZZoJbVp3urVq4O+fvbZZ1FcXIzt27dj0aJFyuV6vR6lpaWJXh4RUcw8uakXz+2wwekR8c3ZWclezqDkzIwxxFwmgE3zIuVwO9Lq7/z85z/Hv/71r7Df1+l0KCwsjMnfirWU6gAsN9XJz88PunzdunUoLi5Gbm4uLrjgAvziF78IW0zkcDjgcPhvWIvFEr8FExFF6KTFAwBoso6CzIyOTfMGo1FpkK3NRo+rB06nMyF/M1ubDY0q8pd0+WRSoOrqatjtduXrxYsXK919Q1m5cmXQ94e6fjylTDAjiiLuuusuLFy4ENOnT1cuX7FiBb7yla+guroadXV1+OlPf4oLL7wQ27dvh14/sFHQqlWr8PDDDydy6UREQzL7sh3yf1PZYN1/AWZmhqJT6zC/Yj5nMyVQygQzt912G/bs2YNPPvkk6PKrr75a+Xz69OmYO3cuqqur8fbbb+Oqq64a8Hvuu+8+3HXXXcrXFosFlZWV8Vs4EVEEzL5ZR/LMo1Q22FwmwB/MOD2AyyNCqw59vbFMp9aN6eAi0VIimLn99tvx5ptv4uOPP0ZFRcWg1y0rK0N1dTVqa2tDfl+v14fM2BARJVO3LyMjZz1S2ZCZGZ0/eLG5RJgYzFCSJfU0kyiKuO222/Dqq6/io48+Qk1NzZA/09HRgcbGRpSVlSVghUREsZFe20yDZ2Z0agEa36sHt5ooFSQ1mLn11lvx97//HS+++CKMRiNaWlrQ0tKCvr4+AEBPTw/uuecebNq0CfX19Vi3bh0uv/xyFBYW4sorr0zm0omIImZ3i7D7yifMabHNNHhmBgiYnM1eM5QCkhrM/OEPf4DZbMbixYtRVlamfLz88ssApCY9e/fuxRVXXIEpU6bg+uuvx5QpU7Bp0yYYjcZkLp2IKGKWgGxMWmVmwpxmAoAsFgErknWCZzSI1f+7pNbMDPWPyMjIwHvvvZeg1RARxUd3QABjcYjwiiJUITrrpgolMxOmzwwQONJg7L6Qa7VaAIDNZlOmQ1N0bDYbAP//y+FKiQJgIqLRrLvP/4LvFaWtGeMgWY9kiyQzo2wzjeFgRq1WIzc3F62trQCAzMzMkOMfaCBRFGGz2dDa2orc3Fyo1eoR/T4GM0REcdZ/a8lsF2FM4UOX8mkm0yCZmSwdt5kAKN3p5YCGopObmxuTDv8MZoiI4qx7QDDjRYVpZO9E48UrisrUbKNu6MzMWN5mAgBBEFBWVobi4mK4XK5kLyetaLXaEWdkZAxmiIjirP8JplQuAu51ivD6lhtRzUwaDM5MBLVaHbMXZopeUk8zERGNBf2DF7kmJRXJa9OqAMMgb3flxnljPTNDqYHBDBFRnHWnUWYm8CTTYMWsmVrp5YPBDKUCBjNERHHW3dcvM5PCjfOsEZxkAgK3mVL330JjB4MZIqI4M/tOB2X7tmbMKTyfyRxB91+AfWYotTCYISKKM7Ovz0xVrlQgmsojDeSTV6aMyDIzY/1oNqUGBjNERHEmBwjVSjCTypkZKTgxDZGZYdM8SiUMZoiI4kwOXqpyNUFfpyK5ANgUZmK2jE3zKJUwmCEiiiOvKCrZjkpfo7xUPppt9q1tsO6/AJvmUWphMENEFEdWhwj55b4qhbaZXt3Xh/lPtWJPc3DXWvnkVe4QwYycmenlaSZKAQxmiIjiSA5cMrQCirKkp9xkH812eUT8cp0Vrb1evH/EHvQ9syOybSYWAFMqYTBDRBRH8sRsk15Qtm7Mdi9EMXlBwHu1drT2SkFLk8UT9D15S2ywUQYAt5kotTCYISKKIznTkZuhUrIdLm9yMxrP77ApnzdZ+g3BjHSbiU3zKIUwmCEiiiO5BsVkEJCpFaDxPeuak1QEvO+UC1tP+Otkmq39MjNKAfDg20yBmZlkZpmIAAYzRERxJc9lMvlmHclbTZYkFQH/fZeUlTmrTAtACma8vmBEFEWlxmeo00xyzYwIwO6O02KJIsRghogojuSgRd62kWceJasL8PaTUlbm+wuyIABweoAOm7RGm0uE2xdj5UaYmZF/jiiZGMwQEcVRd78mdDkBRcCJ5vGKqO+S0ijTijUozpbWItfNyFkkrSo4WAlFrRJgkHoAwuZK/lFzGtsYzBARxZEcIMiZGTmoSUbNTJPVA6cH0KmBcTlqlOdIfW/kupnALSZBGDyYATg5m1IHgxkiojjqX4MizzxKRmamrlMKWqpzNVCrBJQbpWDmpKV/MDN0IAMAmTrp38JtJko2BjNERHHU6atHycsIzswkowD4WKe0xVSTLwUx5TnSmpqVYCayUQYyNs6jVMFghogojuQtnFKjrwBYrplJwjZTXZe0lpo8qdilzChvM0mBVaQnmWTK5GxuM1GSMZghIooTj1fEqR4pQJDrU4y+00w9SQhm5MzMBCUz03+bKbIeM7IsdgGmFMFghogoTtp7vXB7AZUAZS5TlpLNSELNjC8zMyFfysyUG0MXAA/V/VeW4wt6+jfeI0o0BjNERHHS3CO9yJdkq6BRSS/8Wb6i2Z4Eb83Y3SJOmuVtpuCamdYeL5weUTlGnhNhZubsSh0A4KOjjlgvlygqDGaIiOKk2de/pdSXAQGALF1y6kwaut0QIW1zFWRKT/0FmSro1FIX31M9nqgLgJdOMgAAtp10oauPvWYoeRjMEBHFSZNv+6U8IJjJTlIwc8x3LHtCnkbpISMI/uPZTRZv1NtMFSY1Ti/WwCsCa48xO0PJw2CGiChOWvqdZAICMjMJLprtX/wrK5Mb51k8UfeZAYBlE/UAgDW19lgsk2hYGMwQEcWJPCagLCgzIz3tJroAWDmW7Sv+lZX5Aq2T1ui3mQD/VtP6Oifsbp5qouRgMENEFCctvgLgshSomamTG+blBWdmKk3S14fb3FFvMwHA9FINSrJVsLlEbG5wxmi1RNFhMENEFCdyZ92ynIHbTE4P4PQkJqARRRFHOuRtpuDMzHnV0jbR+joHLI7o+swAgEoQcEGN9Dt2nGQwQ8nBYIaIKA6CGuaFyMwAicvOnOrxotsuQi0AkwqCg5nZ47QwGQR020V4fcvJiSIzA0hHz4HkdDUmAhjMEBHFRVuvFx4RUAc0zAMAjUqA3hdP9CSobuZgm38mk14TnHXRqAQsmaBXvjZoAIMm8swMAOT4hmcmY94UEcBghogoLuSuuCXZKqhVwcGBvwg4MZmMg20uAMBpRdqQ379woj+Yiab4VyaPaLAwM0NJwmCGiCgO5OGNgQ3zZIkuAj7ky8ycVqQJ+f0LavRQ++Kt4QQz8raU1cHMDCUHgxkiojiQMzNlKRDMHGyXgpmphaGDGZNBhXkVOt/n0W0xAUAOMzOUZAxmiIjiQDnJZBz4NJvILsAuj4gjvmDmtOLQwQwALJssbTWFCr6GYmTNDCVZ+Hs2ERENm7zNJHfYDZTpm5ydiALgui43XF4pgKoIsRbZN2dlQqNCUDFwpOTMjJWZGUoSBjNERHEw2DZTIguA5ZNMUwr9M5lC0aoFXD87a1h/Q6mZcYrweMUBBc9E8cZtJiKiOGjrlbIuxVkDn2YTWTMjBzNTwxT/xoJ8mgkAehLc2ZgIYDBDRBQX8hZSToiCWjmYScQLv3yS6fQ4BjM6tQCD79ezboaSgcEMEVEc9PjqR+QtpUByAbAtAZOz5R4zU8P0mIkVpXEe62YoCRjMEBHFmMMtwuVLUGTrB8vMxDeL4RVFZXJ3/wGTsSZnoCzsNUNJwGCGiCjGArePsrShgpnEFABbHSLkvzCcZnjRMDIzQ0nEYIaIKMZ6fNmJLK0Q8mRPogqA5aPSeg0GzGSKNaVxHmtmKAkYzBARxZjVF6SE2mICElcALG/5yFmTeJJrZthrhpKBwQwRUYz5i39DBzOJ6gAsBxY5YdYRS/6aGQYzlHgMZoiIYkwu7M0OkxHxbzPFd0tGHvxojHO9DMCRBpRcDGaIiGJs6MyM9NQb/20mX2YmzHZXLHHYJCVTUoOZVatWYd68eTAajSguLsbKlStx6NChoOuIooiHHnoI5eXlyMjIwOLFi7Fv374krZiIaGhKzUyYYEbOzNjivs3ka9yX0JoZZmYo8ZIazKxfvx633norNm/ejDVr1sDtdmP58uXo7e1VrvPYY4/h8ccfx5NPPomtW7eitLQUy5Ytg9VqTeLKiYjC61W2mQYPZlxeqSdNvFjs0u82JiAzY2RmhpIoqv7Whw4dwj//+U9s2LAB9fX1sNlsKCoqwqxZs3DRRRfhS1/6EvT6yCeurl69OujrZ599FsXFxdi+fTsWLVoEURTxxBNP4P7778dVV10FAHj++edRUlKCF198Ed/73veiWT4RUULI20fGEN1/Af/UbEAqAo7XsWmlZiYR20wG1sxQ8kSUmdm5cyeWLVuGmTNn4uOPP8a8efNwxx134Oc//zm+8Y1vQBRF3H///SgvL8ejjz4Kh8MxrMWYzWYAQH5+PgCgrq4OLS0tWL58uXIdvV6PCy64ABs3bhzW3yAiijelZiZMEKFR+WcZxbMI2F8zk4htJunfauWgSUqCiDIzK1euxL333ouXX35ZCTRC2bRpE37zm9/g17/+NX784x9HtRBRFHHXXXdh4cKFmD59OgCgpaUFAFBSUhJ03ZKSEhw/fjzk73E4HEHBlMViiWodREQjZZW3mQY5Ep2tU8Hu9sa1CFg+mp2IzIyJmRlKooiCmdraWuh0uiGvd8455+Ccc86B0+mMeiG33XYb9uzZg08++WTA9wQh+IEoiuKAy2SrVq3Cww8/HPXfJyKKlcGGTMqydALabfHtNSM3zctJxNFsnb9mZrDnaKJ4iOgeHkkgM5Lr33777XjzzTexdu1aVFRUKJeXlpYC8GdoZK2trQOyNbL77rsPZrNZ+WhsbIxqLUREI6XUzAySEUlEF2B5m8mYkKZ50suJV4x/M0Ci/qIqAJZt2bIF69atQ2trK7ze4JTi448/HvHvEUURt99+O1577TWsW7cONTU1Qd+vqalBaWkp1qxZg1mzZgEAnE4n1q9fj0cffTTk79Tr9VEVIRMRxdpQfWYCv2dzxXObKXHjDAwaQKuSTmhZnSKy+TRMCRR1MPPII4/gJz/5CaZOnYqSkpKgVGK0acVbb70VL774It544w0YjUYlA2MymZCRkQFBEHDHHXfgkUceweTJkzF58mQ88sgjyMzMxNe//vVol05ElBDWIToAA/7J2T2JKAA2xD8zIwgCjHoBnX0iLHYvyozquP9NIlnUwcxvf/tbPPPMM7jhhhtG/Mf/8Ic/AAAWL14cdPmzzz6r/P4f/ehH6Ovrwy233IKuri4sWLAA77//PoxG44j/PhFRPESSmZGPZ8dzSyaRR7MB6dRUZ5+HvWYo4aIOZlQqFc4777yY/HFRHPoOLwgCHnroITz00EMx+ZtERPHWG0HNTLyHTTrcIhxu6fNEHM0GAoZN8kQTJVjU9/A777wTTz31VDzWQkSU9jxeEb2+OpisIU4zAfErAA4cKzBYhiiW/CMNmJmhxIo6M3PPPffg0ksvxcSJEzFt2jRotdqg77/66qsxWxwRUbrpDSjojaQAOF6ZGWvAVpdalZhgxj/SgJkZSqyog5nbb78da9euxZIlS1BQUMBeAkREAeR6GZ0ag44pkLM28eoAnMiJ2TI5M8OaGUq0qIOZF154Aa+88gouvfTSeKyHiCit9UTQ/RdI3DZTIo5ly+SaGW4zUaJFfS/Pz8/HxIkT47EWIqK0p2zvDBFE5GVI3++0xTczk6iTTNLfkv5NZhYAU4JFHcw89NBDePDBB2Gz2eKxHiKitCZnWobKzJRkS0+/rb3xeeGXMzOJOskE+AO0jjgFaEThRL3N9L//+784evQoSkpKMH78+AEFwDt27IjZ4oiI0k0kPWYAoNgXzJzq8cRlllEih0zKirKkf1M7gxlKsKiDmZUrV8ZhGUREo4NcMzNUrUpxltQh1+GWtoRMMe7SK/d6ScSQSZkczLTFKdtEFE7UwcyDDz4Yj3UQEY0KkW4zGbQCTAYBZruI1h4vTDEOOiwRNO6LNSUz0xufbBNROMMaNAkA27dvx4EDByAIAqZNm6YMgiQiGsuUbaYIgojiLDXMdjdO9XgwuXDYT8chyZmZRAYzhZlSMGN3S0FdIv82jW1RP3paW1vxta99DevWrUNubi5EUYTZbMaSJUvw0ksvoaioKB7rJCJKC/6j2UNnWkqyVajtiE8RsFXpM5O4baZMnQrZOgE9ThFtvd6EHgunsS3qe9rtt98Oi8WCffv2obOzE11dXfj8889hsVjwgx/8IB5rJCJKG5EWAAPBRcCxZknwkEmZnJ1pZ90MJVDUmZnVq1fjgw8+wOmnn65cNm3aNDz11FNYvnx5TBdHRJRurM4otpmypSLg1p7RkZkBgKJsFeq7PSwCpoSK+l7u9XoHHMcGAK1WC6+Xd14iGtsi7QAMAMW+gtl4BjOJzsz4TzTFPttEFE7UwcyFF16IH/7wh2hqalIuO3nyJO6880584QtfiOniiIjSTY8SRERSMyNlZuK5zZTozExhpvRvYq8ZSqSo7+VPPvkkrFYrxo8fj4kTJ2LSpEmoqamB1WrF7373u3iskYgobUR6NBuIXxdgURQDtpmSlJmJQ7aJKJyoa2YqKyuxY8cOrFmzBgcPHoQoipg2bRqWLl0aj/UREaUVeYxAZDUz8ekC3OsU4fXNekzWNpOcmWnt8SDHoIJhkAniRCMVVTDjdrthMBiwa9cuLFu2DMuWLYvXuoiI0lKvkpkZOvEdry7AXX1SIKFTAxnaBJ9mCqiZaex24wtPt+PsKh1e+Ep+QtdBY0tU20wajQbV1dXweFjYRUTUn90lwmyXghl56OJgDFpB2QaKZRFwo1l6jh6Xo054F94iX4DW1uvFJ8edcHqA2nZ3QtdAY0/UNTM/+clPcN9996GzszMe6yEiSltHOtwQAeRnCMjPiCyIkIuAW2N4+kcOZipM6pj9zkj5Rxp4savJBcCfrSKKl2FNzT5y5AjKy8tRXV2NrKysoO9zajYRjVWHfBmIKYXaiDMichfgU3HIzFTlJj6YKfA1zXN5gY/rHQAAm0vkrCaKq6iDmSuuuIJ3SCKiEA77gpmpRZE/tcajC3Bjt/S7KpOQmdFr/AM0m61SgOb2Ak4PoI/t+CkiRdR3rYceeigOyyAiSn8H26RtlSlRDI2MRxdgOTNTaUpO9FCUpYLZHhyc2Vwi9DzRRHESdc3MhAkT0NHRMeDy7u5uTJgwISaLIiJKR3Jm5rRoMjNx6ALckMRtJsDfOC8Q62YonqIOZurr60OeZnI4HDhx4kRMFkVElG7Mdq+yrTI5isxMrAuA+1yiMuQxGdtMgL8IOJDNxSZ6FD8RP+LefPNN5fP33nsPJpNJ+drj8eDDDz9ETU1NbFdHRJQm5OPH5UZVVCME5ILZjhi1/z/hy8oYdULM+tZEK1Qww8wMxVPEwczKlSsBAIIg4Prrrw/6nlarxfjx4/HrX/86posjIkoXykmmKLaYAH+nYHmm00g1mqV1VOYmvseMTG6cJwAoz1HhpMULm4vBDMVPxI86eSJ2TU0Ntm7disLCwrgtiogo3Rxq851kKtRG9XNG3wynnhhlLhqSeJJJJmdmJhdqkK0TcNLiZWaG4irqmpm6uroBgUx3d3es1kNElJYOtUd/kgnwT9e2uUR4vCN/wfefZEpeMLN4gh7nVetwy4IsZPmCNRuDGYqjqIOZRx99FC+//LLy9Ve+8hXk5+dj3Lhx2L17d0wXR0SUDkRRHFaPGQDKiz0Qm+yMEswk6SQTII00+MfV+Vh5RgYyfbOhernNRHEUdTDzpz/9CZWVlQCANWvW4IMPPsDq1auxYsUK3HvvvTFfIBFRquvqE9HVJ71YT8qPLpjRawTofHGHNQZ1M/I2U1WSesz0x8wMJULU9/bm5mYlmPnPf/6Dr371q1i+fDnGjx+PBQsWxHyBRESprtsu1RQa9QIMw5hSnaNXod3mRY/TC2D4GRVRFJXTTMnMzARiZoYSIerMTF5eHhobGwEAq1evxtKlSwFIDyJO0yaiscjikIIZeQJ2tLJ92YuRZma67aKyVVWRkxrBTJbOVxPkZJ8Zip+oMzNXXXUVvv71r2Py5Mno6OjAihUrAAC7du3CpEmTYr5AIqJUJwchxij6ywSK1fFsOStTlKUaVoYoHpiZoUSIOpj5zW9+g/Hjx6OxsRGPPfYYsrOzAUjbT7fcckvMF0hElOosdumFeviZGSkIso4we9HZJ/18YYimdcnCmhlKhKiDGa1Wi3vuuWfA5XfccUcs1kNElHaUbSbD8IIIY4wyM92+YCZ3mOuIB2ZmKBGiDmYCxxoEEgQBBoMBkyZN4lgDIhpTRpqZkYOZkdfMSMFMXkZqbDEBzMxQYkQdzKxcuRKCIEAUg++Y8mWCIGDhwoV4/fXXkZeXF7OFEhGlqpFmZuRtpp4RbjN1+46Hm5iZoTEm6nv8mjVrMG/ePKxZswZmsxlmsxlr1qzB/Pnz8Z///Acff/wxOjo6Qm5FERGNRhZfRiVHN8LMzAizF3JmJpW2mZiZoUSIOjPzwx/+EH/+859x7rnnKpd94QtfgMFgwHe/+13s27cPTzzxBG688caYLpSIKFVZR5qZUbaZRpaZMadgMJOpY2aG4i/qe/zRo0eRk5Mz4PKcnBwcO3YMADB58mS0t7ePfHVERGkgVqeZRloALHchTqmaGS37zFD8RR3MzJkzB/feey/a2tqUy9ra2vCjH/0I8+bNAwDU1taioqIidqskIkphIz7NFKPJ2fI2kymDmRkaW6LeZnr66adxxRVXoKKiApWVlRAEAQ0NDZgwYQLeeOMNAEBPTw9++tOfxnyxRESpSKmZGeFpppFmZswpeDQ7y1cA7HADbq8Y/YsOUQSivl9NnToVBw4cwHvvvYfDhw9DFEWcdtppWLZsGVQq6QG0cuXKWK+TiChlWezyOIPh1sz4muaFqZk50OrCtpMuXHtWBlRC+ICpKwWPZmcGFEXbXCJyGM1QHAzrbiUIAi6++GJcfPHFsV4PEVHaUTIzhhHOZgqzzfTABxZsPeHChHw1zqvWh7yOVxRh9tXupFJmRqcWoFUBLq90oonBDMXDsO5WH374IT788EO0trbC6w1+J/HMM8/EZGFEROnA5RFhc41sNpOyzRQmmGntkZ5nj3d5cF516N9hdYjw+n58uLU78ZKpE2C2i+h1ikBmsldDo1HUwczDDz+Mn/3sZ5g7dy7KysogDJLyJCIa7QIDEOOwa2ak4KPPJcLlEaFVB/8e+W+0WD1hf4d8LDtDK8CgSa3n5SytFMzYWARMcRJ1MPPHP/4Rzz33HK677rp4rIeIKK3IW0xZWgEa1fCCiKyAupJep4jcfjUvci1Nc0/4483KsexhbnXFk3KiiY3zKE6izkU6nc6ghnlERGPZSOtlAKmuRO97a9l/crbDLcLpS8gMlplJxWPZMqXXjIu9Zig+or7Xf+c738GLL74Yj7UQEaUdf8O8kQUR8lZT/+PZgdtYLdbwwUAqTsyWMTND8Rb1NpPdbsef//xnfPDBBzjzzDOh1WqDvv/444/HbHFERKkuFpkZQGqc1947cHJ24HHtSGpmUulYtkzuNSPVzKTe+ij9RR3M7NmzB2eddRYA4PPPP4/1eoiI0oq/Yd7IMiLZShfg4OxLYKbG6hRhdXhDnprqSsGJ2TJmZijeog5m1q5dG7M//vHHH+N//ud/sH37djQ3N+O1114Larh3ww034Pnnnw/6mQULFmDz5s0xWwMR0UjImZThnmSS+RvnBb/gW/o10jvVEzqYScWJ2bLgzAxR7MXkXu/1evHWW29F3fm3t7cXM2fOxJNPPhn2OhdffDGam5uVj3feeWeEqyUiGpk+pwcPv12LjY1izDIzRmVydviaGSD8VlMqTsyWMTND8TaiXoy1tbV45pln8Pzzz6OrqwsXXXRRVD+/YsUKrFixYtDr6PV6lJaWjmSZREQx9dLWBjy76QQ2FQBnV8WmZibcNlP/4KY5TDDTnYITs2XMzFC8RR3M9PX14V//+heefvppbN68GR6PB7/5zW9w4403Ijs7O+YLXLduHYqLi5Gbm4sLLrgAv/jFL1BcXBz2+g6HAw6HQ/naYrHEfE1ENLat/rwFAHC4E6jJj/FpJmf/zExwcBPuRFNXCh/NjjQzI4oiG7HSsER8r9+yZQu++93vorS0FE8++SS+9KUvobGxESqVCkuXLo1LILNixQr84x//wEcffYRf//rX2Lp1Ky688MKgYKW/VatWwWQyKR+VlZUxXxcRjV3tPQ5sre8EAHhF4LMTbgDDn5gtCzc5u//X4TIzqTgxW5alG7rPzPMb6zHz4fexs6ErUcuiUSTizMy5556L22+/HVu2bMHUqVPjuSbF1VdfrXw+ffp0zJ07F9XV1Xj77bdx1VVXhfyZ++67D3fddZfytcViYUBDRDGzZv8pZQYSAHT2ydtMsTnN1L/gN7DDcK9LxKkwXYC7U/hodoZ28MxMm9WBR1cfhM3pwdqDrZhVlZfI5dEoEHEwc+GFF+Lpp59Ga2srrrvuOlx00UUJTweWlZWhuroatbW1Ya+j1+uh14eeKktENFLv+raY8jO16LS5lMtHfJpJN/g206RCDXY3u0JmZryiiO4UnJgtC1Uzs6W+G/tOncLXF1ThqbVHYPO1OT7ZbU/KGim9RRzMvP/++2hsbMSzzz6Lm2++GX19fUrmJFFBTUdHBxobG1FWVpaQv0dEFMjc58LGI+0AgFsuqMZ/v3tE+V6sTjP131aSC4AnF0jBTKjTTKk8MRsYWDPjFUXc8tI+tPc48fLWRhxt61Gu29Tdl5Q1UnqL6l5fWVmJBx54AHV1dfjb3/6G1tZWaDQaXHHFFfjxj3+MHTt2RPXHe3p6sGvXLuzatQsAUFdXh127dqGhoQE9PT245557sGnTJtTX12PdunW4/PLLUVhYiCuvvDKqv0NEFAvrDrXC7RUxpSQbV8wsCfperGpmrI7QTfMmF0jvPTv7RNjdwQFPKk/MBgIyM75gZn8b0N7jBAAcbLHC5RFRmK0DADSbGcxQ9IYdwi9btgz//Oc/0dTUhNtvvx3vvvsu5s2bF9Xv2LZtG2bNmoVZs2YBAO666y7MmjULDzzwANRqNfbu3YsrrrgCU6ZMwfXXX48pU6Zg06ZNMBqNw102EdGwneiSXmhnVuSiKFuH8oBzDyOvmRl8m6nCpIbBl0tv7QnOznSn8MRsICAz49tm2tAg/Xd+TT4WTy1CQZYO/71yBgCgyWyH18sj3BSdEfWZAYC8vDzcfvvtuP3226POzCxevBiiGP5O+9577410eUREMWOxSzUypgxpJt2MEqDJt0My8g7A0s93271BR5QtAR2Gy4xq1HV50Gz1oirX/7OpfCwb8GdmepwivKKIj33BzGVnluGb54yHKIpweUQIAuB0e9HR60SRkbWPFLmI7vkNDQ0R/bLZs2cDAE6ePDn8FRERpShLn3QM22iQgpkzi6UX6QytAJ16ZMFMpUmNDK2Arj4R6+qcyuVypiZbL6AkWw0AaLYEZ2bkbSbTCOt24qXEqIZRJ6DPJeKVfU5sa5IuP39yEQCp7lKnUaHYF8CwboaiFdE9f968ebjpppuwZcuWsNcxm834y1/+gunTp+PVV1+N2QKJiFKF1ZeZMfr2e2b4ghnTCLMygNQ079qZGQCAJzf1KFlruYYmR69CQab0lG3uV1cjF9bKx7tTjUEj4NvzsgAAD3zYB5cXqMg1YHxBZtD1ynOlfz+DGYpWRNtMBw4cwCOPPIKLL74YWq0Wc+fORXl5OQwGA7q6urB//37s27cPc+fOxf/8z/8MOaKAiCgdWey+Bnm+baZzKoCrpmkxpyIjJr//pnlZeGGnDdtPurC50YmzK3VKAXC2TgjbSVf+OitFgxkAuHFuJp7Z1qtsm50/KX/ASdjy3AzsbOjGSQYzFKWIMjP5+fn41a9+haamJvzhD3/AlClT0N7ervR7ufbaa7F9+3Z8+umnDGSIaNTqn5nRqgU8viIL156VOdiPRazEqMZXZ8jZmV7Y3YDHF7dk64SwM47SIZjJ0avw3flZyteLJg9sjFduMgAAms3sNUPRiaoA2GAw4KqrrgrbfZeIaDSzypkZX81MPHxvQRb+vqsPnx53otEs/T0BUqCS0e+Is6zXd+IpVbeZZDfMzsTfd/bC4RZx7oQQwQy3mWiYRnyaiYhorLD0BWdm4qHSpEFptgotPV7saZGCmWy9AEEQlMxL/8yMXCScqUvNAmBZtl6Fd64zwu3sUU6EBWIwQ8OV2vd8IqIUkojMDACMz5OCpc9bpOBJbsiXGWbGUaoXAAfKz1ShOCv0Osf5ghmONKBoMZghIoqAy+NFn0s6Ep2TEd+kdk2edAR77ykpmJEb6snBTF8a1sxEQs7MtPc44HCHng5OFAqDGSKiCMhZGQDI1sc3mJEzM/tb/dtMwMBOujL563QPZvIytTBopZelFhYBUxQYzBARRUA+yZSpU0Ojju9T5/h8KTMjZ2Dk7sL9ZxzJ0qUAeCiCIKDcJG81sW6GIjesR+Tf/vY3nHfeeSgvL8fx48cBAE888QTeeOONmC6OiCgZ/r2tEdc/s0UZXwD4u//Gu14G8G8zyeRtpowwR7OVAmBt+r8/9RcBMzNDkYv6nv+HP/wBd911Fy655BJ0d3fD45H2NXNzc/HEE0/Een1ERAn3h/VHsf5wGzYe6VAu699jJp6qTMF/Q8nM+IIamyu9OgBHozxX6jXDE00UjaiDmd/97nf4y1/+gvvvvx9qtf/dw9y5c7F3796YLo6IKNE8XhGNnTYAgLnPPyNJztLkhDhSHGsGrYByo//p2aiTg5n07QAcqVLfNhMb51E0og5m6urqMGvWrAGX6/V69Pb2xmRRRETJ0tTdB5ev7a65L2CbyS4PmUxMey65CBiQ+rMAobeZRFEcVcFMSY40bLLNymCGIhd1MFNTU4Ndu3YNuPzdd9/FtGnTYrEmIqKkqe/wvynrtgXWzMjbTPHPzADA+IC6mf4FwE4PlICrzyVCDm1GwzZTsVHaZmq1OpK8EkonUb/FuPfee3HrrbfCbrdDFEVs2bIF//znP7Fq1Sr89a9/jccaiYgSpr7DpnwemJnxN8xLTGamJjAzows+mg1I2RmTWlCKfwX4MzfprNgoZWZaLQxmKHJRPyq/9a1vwe1240c/+hFsNhu+/vWvY9y4cfjtb3+Lr33ta/FYIxFRwtS3+zMzoYKZ5GRmpCS6Ti1AowLcXikjYzIE18v0n0Kdjop920ztPQ54vSJUqvT/N1H8Destxk033YSbbroJ7e3t8Hq9KC4ujvW6iIiS4nhH6GDGXwCc+JoZeZsJkLoAWxz+OpnRVC8DAIXZeggC4PaK6LQ5UZitT/aSKA0MqwC4trYWAFBYWKgEMrW1taivr4/p4oiIEq0uIDNjCcrMJLZmpjJXDTk8yQ4YINl/2GTPKAtmtGoV8jN1ALjVRJGLOpi54YYbsHHjxgGXf/bZZ7jhhhtisSYioqSQjmX7+5sEZWb6ElszY9AImFuhRbZOCNpy6n+iSe45MxqKf2VFvrqZUzzRRBGKOpjZuXMnzjvvvAGXn3322SFPORERpYtmcx+cHn9DuqCaGYdvmylBmRkA+MfV+fj0+0UwGQIyM/1GGoym7r+y4hzpRFMbMzMUoajv/YIgwGq1DrjcbDYr3YCJiNJRfbt0kknuJWPuc8HrlYIFa4L7zABSwW9gIAMEDpuUgq7R1P1XViKfaGJmhiIUdTBz/vnnY9WqVUGBi8fjwapVq7Bw4cKYLo6IKJHkHjNnVpgAAF4R6HFKQUyi+8yEI2dg5G2m0VYADPhPNLHXDEUq6rcYjz32GBYtWoSpU6fi/PPPBwBs2LABFosFH330UcwXSESUKPKx7KklOdha3wWn2wuzzQWjXuPvM5Og00zhZIbZZhpVwYzcOI/bTBShqDMz06ZNw549e/DVr34Vra2tsFqt+OY3v4mDBw9i+vTp8VgjEVFCyA3zagozYfLNYDL3udDn8sDt225KdmZGmc+kZGZGXwFwMbeZKErDeotRXl6ORx55JNZrISJKKrnHTHVBFkwZWrRZHbD0uZSsjEoAsnTqwX5F3MmZmb5+20yZutFUAOzfZjL3uXDny7uwbFoJrplfleSVUaoaVjDT3d2NLVu2oLW1FV5v8Cj6b37zmzFZGBFRorX4JjWX52YgNyAzE9hjJtlddjP7Tc4ejQXAgfOZXt95Eh8dbMXmYx249MyyhJ4mo/QRdTDz1ltv4dprr0Vvby+MRmPQA1sQBAYzRJSWXB4vrA4pA5OfpVO2mbr7XDD3Jf4kUzj9j2aPxgJguc+M0+3FW7ubAAA2pwf/3nYC315Yk8ylUYqKOi95991348Ybb4TVakV3dze6urqUj87OzniskYgo7uQJ2YIAmDK0QTUzcmYmFbIC/ZvmjcYCYINWrTQn3Ha8S7n8b5vqlaPyRIGiDmZOnjyJH/zgB8jMzIzHeoiIkqLb5gQgBSxqlYCcgGDGkoQeM+H0H2cgFwBnjYKJ2YFKfI3zAKAwWwejQYP6DhvW17YlcVWUqqIOZi666CJs27YtHmshIkqaLl9mJi9TCmJCZWaSfZIJCDiaPaDPzOgpAAb8RcAAcMGUYnxlTiUA4IWN9UlaEaWyqN9mXHrppbj33nuxf/9+zJgxA1pt8IP7i1/8YswWR0SUKF2+zIzJN+QwMJhR5jIluccMMLDPjHxEezQVAAP+ImAAWDSlEJOKs/HMp3XYHrDtRCSL+pF50003AQB+9rOfDfieIAgcaUBEaUneZuqfmbGkWM1MuHEGo6lmBvD3mhEE4PzJRfCKvrESDjc8XhFq1ej699LIRB3M9D+KTUQ0GnQr20wDMzN1vs7AgXUcyZLlG2fQ5xIhiuKoDWbk/9dnVuQiP0sHl28AqCgCVrsLub7biQgYZp8ZIqLRRq6ZyfVlZuT/dttcaPPNCJpVlZuUtQUK7DNjc4mQz/aMtmBm5axx2FrfiW+eMx4AoFWrkK3XoMfhRreNwQwFG1Yw09vbi/Xr16OhoQFOpzPoez/4wQ9isjAiokTybzMFZ2aauvvg9opQCf4BlMkUWAAsZ2WEgMtHi/wsHf7wjTlBl5kytFIw4xv6SSSLOpjZuXMnLrnkEthsNvT29iI/Px/t7e3IzMxEcXExgxkiSktdYWpm5JlMU0tzkKlLfjJbzsA4PYDZ7t9iSnZn4kTIzdTiZHefEngSyaI+y3fnnXfi8ssvR2dnJzIyMrB582YcP34cc+bMwa9+9at4rJGIKO7820xSZkbuMyNLhS0mwN80DwDaeqUDF6NtiymcwDomokBRBzO7du3C3XffDbVaDbVaDYfDgcrKSjz22GP48Y9/HI81EhHFXf9tJoNWDb3G/xQ5qzI3GcsaQKcW4KsBRluvr2HeGAlm5DomBjPUX9TBjFbrH7RWUlKChoYGAIDJZFI+JyJKN/0LgAF/JgAAZlXlJXxN4chFwK09o7P7bzimDCnQlE+eEcmi3gCeNWsWtm3bhilTpmDJkiV44IEH0N7ejr/97W+YMWNGPNZIRBRXoij6MzNZ/lMypgwtWq0O5Bg0mFCYlazlDZCpFWC2i2hVtplGV/ffcAJPmBEFivoR8Mgjj6CsrAwA8POf/xwFBQW4+eab0draij/96U8xXyARUbz1Oj1weaRi2ryAzIz84jmzMheqFGrSJp9ckreZRlv333BylUnmLACmYFFnZubOnat8XlRUhHfeeSemCyIiSrSuXunFUadRIUOrVi6XtzVSaYsJ8NfIHG73DcDUj41gRikAZmaG+ok6M3PhhReiu7t7wOUWiwUXXnhhLNZERBR3oijirxuOYfXnzQHdf7VBR5yvmV+JBTX5+MqcimQtMyT5RNP+VimYuaBGP9jVRw0WAFM4UWdm1q1bN6BRHgDY7XZs2LAhJosiIoq3t/Y047/fPgCDVoXfXzsbgP8kk+wLp5fgC6eXJGN5g5JHGgBASbYKl5yW/DELiaAUADOYoX4iDmb27NmjfL5//360tLQoX3s8HqxevRrjxo2L7eqIiOLA7vLg0XcP+j734rNjnQCCTy+lssyAGpnrZmVCpx4b20wsAKZwIg5mzjrrLAiC1GUy1HZSRkYGfve738V0cURE8fD8xnqc7O5Tvl5/uA3AwMxMqpILgPUa4JqZmUleTeL4t5mcEEVxTHQ9pshEHMzU1dVBFEVMmDABW7ZsQVFRkfI9nU6H4uJiqNXqQX4DEVHymW0uPLn2CABgQlEWjrX14mCLFQCQl5UemZmSbGmb6cppGSjIHBvHsgF/5szlEWFzepClT/54CUoNEd8TqqurAQBerzduiyEiircdjV2w2t2oLsjEPcun4pZ/7FC+ly6TmL81JwvF2Wp86YyxUSsjy9CqoVOr4PR40d3nYjBDiqhD+ueffx5vv/228vWPfvQj5Obm4txzz8Xx48djujgiolhrszoAAOMLsnBWvxEFgT1mUll+pgrXzcpE5hhplicTBAGmTB7PpoGG1TQvIyMDALBp0yY8+eSTeOyxx1BYWIg777wz5gskIoql9h4pmCnM1qPMZEBhtj8bky6ZmbGMjfMolKiDmcbGRkyaNAkA8Prrr+PLX/4yvvvd72LVqlVRH83++OOPcfnll6O8vByCIOD1118P+r4oinjooYdQXl6OjIwMLF68GPv27Yt2yUREinar9CJYaNRBEAScWZGrfC9dCoDHslxmZiiEqIOZ7OxsdHR0AADef/99LF26FABgMBjQ19c32I8O0Nvbi5kzZ+LJJ58M+f3HHnsMjz/+OJ588kls3boVpaWlWLZsGaxWa7TLJiIC4M/MFGVLjebOrDAp30uXbaaxjL1mKJSoq6eWLVuG73znO5g1axYOHz6MSy+9FACwb98+jB8/PqrftWLFCqxYsSLk90RRxBNPPIH7778fV111FQCpXqekpAQvvvgivve970W7dCKioG0mIDiY4TZT6pNPNLHXDAWKOjPz1FNP4ZxzzkFbWxteeeUVFBQUAAC2b9+Oa665JmYLq6urQ0tLC5YvX65cptfrccEFF2Djxo1hf87hcMBisQR9EBHJBgYzucr3CrIYzKQ6pXFeQM3M1vpOPLX2CLxeMVnLoiSLOjOTm5sbclvo4YcfjsmCZHKH4ZKS4FbiJSUlg56aWrVqVczXQkSjR3uPv2YGkIKa/3fxabC7PMhjMJPyckMMm7z/tb04fKoHs6vycM7EgmQtjZIoomBmz549mD59OlQqVdBYg1DOPPPMmCxM1r/D41BdH++77z7cddddytcWiwWVlZUxXRMRpSe3x4sumy+YyfYPZ7x58cRkLYmi1H/YpNPtxdG2XgBAY5cN54DBzFgUUTBz1llnoaWlBcXFxcpYA1H0p/PkrwVBgMfjicnCSktLAUgZmrKyMuXy1tbWAdmaQHq9Hnr92JggS0TR6ex1QhQBlcCTS+nK5Lvd5JqZhs5eeHzbSy1me9LWRckVUTBTV1enjC+oq6uL64JkNTU1KC0txZo1azBr1iwAgNPpxPr16/Hoo48mZA1ENLq0+epl8rP0UKs41ycdKQXAvszMkdZe5XstFgYzY1VEwYw8yqD/5yPV09ODI0eOKF/X1dVh165dyM/PR1VVFe644w488sgjmDx5MiZPnoxHHnkEmZmZ+PrXvx6zNRDR2KHUy2QzK5Ou/DUz0m15tK1H+R4zM1IpxtOf1CHHoMVX542dEouIgpk333wz4l/4xS9+MeLrbtu2DUuWLFG+lmtdrr/+ejz33HP40Y9+hL6+Ptxyyy3o6urCggUL8P7778NoNEb8N4iIZO2+UQZFRm5Fpyv/aSYpM3O01R/MNDOYwaZjHfjvtw9AJQCXzyxHhm5sDICOKJhZuXJl0NehamZk0dTMLF68OOj39CcIAh566CE89NBDEf9OIqJw+h/LpvQjnzizOT0w21xBmZlT3GbCn9YfAwB4ReB4Zy9OK81J8ooSI6I+M16vV/l4//33cdZZZ+Hdd99Fd3c3zGYz3nnnHcyePRurV6+O93qJiIZNHjLJbab0lWPQYlJxNgBgw5E25SQTIBV4212xOYSSjg40W7D+cJvydX177yDXHl2i7jNzxx134I9//CMWLlyoXHbRRRchMzMT3/3ud3HgwIGYLpCIKFaYmRkdlkwtwpHWHry8tRE9DjfUKgEalQCH24tTFjuqC7KSvcSk+NP6o0Ff17XbkrSSxIu6A/DRo0dhMpkGXG4ymVBfXx+LNRERxYW/AJjBTDpbMrUYALChth0AUJ2fiTKTAcDYLQI+2d2Ht/Y0AwCWni61LxlLmZmog5l58+bhjjvuQHNzs3JZS0sL7r77bsyfPz+miyMiiiUlM8MC4LQ2d3w+sgIKWycUZaNUDmbGaN3M3zcfh8cr4pwJBbjsTKk3W10Hg5mwnnnmGbS2tqK6uhqTJk3CpEmTUFVVhebmZjz99NPxWCMRUUz4t5lYM5POdBoVFk4uVL6eWJyFMlMGgLF5osnu8uClLQ0AgBvOG4/qgkwAwPExFMxEXTMzadIk7NmzB2vWrMHBgwchiiKmTZuGpUuXDjpmgIgomTxeEZ290jYTj2anvyVTi/HevlMAgIlF2RAgvf6MxW2mN3c3ocvmwrjcDCw9vQRWu3Rs/ZTFAZvTjUyd9FJ/pLUHHx44hW8vrIFGHXUuI6VFHcwA0pHp5cuXB020JiJKZZ29TnhFQBCAfI4ySHuLfXUzgBTM9DmlU0xjLZgRRRHPb6wHAFx3TjXUKgG5mTrkZmrRbXOhvt2GaeU5cHu8+M7zW1HfYUNFXiYuPbNs8F+cZkZXaEZEFIa8xZSfqRt170rHolKTAd84uwrnTy7EjHEmpWameYzVzOw5Yca+Jgv0GhWunuvv+Dved6Kr3rfV9NrOk6jvkE43He8cfdtPw8rMEBGlGx7LHn3+e+UM5fPSHCmYOTXGMjOHT1kBAPNr8pWGggBQU5iFXY3dqGvvhdvjxZNr/aODRmP2im9PiGhMkN+VFhq5xTQayUezW612uD3eJK8mceR2A/3rwAKLgF/deRLHO/w9Z0ZjkTSDGSIa9URRxD8/k057LJxUlOTVUDwUZEuT0L2ifzr6WCBnHIv6ZRxrCqVtpu3Hu/DY6oMAgAU1+QBG59iHEQUzfX19sFgsQR9ERKnms7pO7G+2wKBV4Zr5Y2eS8FiiVgko8WUnRuM2Sjjhtk/lmpmjbb1o73HitFIj7r1oKgBmZgAANpsNt912G4qLi5GdnY28vLygDyKiVPPMJ3UAgC/NrkAuTzKNWkoR8Ch8sQ6nLcwk+PEBIx3KTAY89635GO/L1rT3OOAaZVtxUQcz9957Lz766CP8/ve/h16vx1//+lc8/PDDKC8vxwsvvBCPNRIRDVtDhw1rDkj9SL513vjkLobiSs5OdPj6CY0F4TIzpkwtZlbmIj9Lh2e/NQ+lJgPyM3XQqVUQRaDVOrq24qI+zfTWW2/hhRdewOLFi3HjjTfi/PPPx6RJk1BdXY1//OMfuPbaa+OxTiKiYXlrTxNEETh/ciEmFRuTvRyKI1OGFgBg6XMleSWJo8wbC1HY/urN58Lp9iLDN/pBpRJQnKPHia4+tJj7MC43I6FrjaeoMzOdnZ2oqakBAOTk5KCzsxMAsHDhQnz88cexXR0R0Qgda5N6apw9oSDJK6F4k4MZ8xgJZtweL7ps4YenqlWCEsjIykbpVlzUwcyECROU6djTpk3Dv/71LwBSxiY3NzeWayNKCFEUk70EiqPGTulIakXe6HkXSqEpwYxtbAQznb1OiKIUtORFWAtW6pthNdqKpKMOZr71rW9h9+7dAID77rtPqZ258847ce+998Z8gUTxZHd5cdsL7bj3n+3Y2zi69pBJ0uALZqryM5O8Eoo3U+bYyszIdS/5WTqoVZHNRizNGZ0nvqKumbnzzjuVz5csWYKDBw9i27ZtmDhxImbOnBnTxRHF25FTLrSYpZkuD73WhSWnZ+C2ZaYkr4pixe7yoMXXU4PBzOin1MzYx0YwM5yu1nJmZrSNfRjxOIOqqipUVVXFYi1ECdfQ4QYA5GSoYOnzYu2BPnxzoRE5GewnORqc6OoDAGTp1MjP4pHs0S5njNXMKMW/2ZHft+WamdGWmYn4Gfujjz7CtGnTQjbGM5vNOOOMM7Bhw4aYLo4o3hp9wczSMzKUAKazx5PMJVEMNXZJW0yV+ZkQhMjS8JS+cgxjLZgJ3f13MCU5YzyYeeKJJ3DTTTchJydnwPdMJhO+973v4fHHH4/p4ojirbFTCmaqCjTIz/IFM72jq5nUWNbIepkxZaydZgrXMG8wcmbmlMUOr9d/+MFsc8HjTd/DEBEHM7t378bFF18c9vvLly/H9u3bY7IookQQRVHZZqos0CA/WzrCyMzM6NHQwWBmLJGDGavdndYvzJEaTs1MkVEPlQC4vSLae6Wf39HQhbN+/j5WvXMgLutMhIiDmVOnTkGr1Yb9vkajQVtbW0wWRZQInb1e9DpEqARgXC4zM6ORcpKpgMHMWCAHMwBgHQNFwEowE8UkeK1apQQ/p8zSz68/1AZRBF7e2giHOz3fzEUczIwbNw579+4N+/09e/agrKwsJosiSgQ5K1Oeq4ZWIwQEM+n5YKaB5GCmkpmZMUGnUSFDK2VYx8JWU7s1fMO8wfgb50kF8kdaewAAVocbnx5pj+EKEyfiYOaSSy7BAw88ALt9YNFQX18fHnzwQVx22WUxXRxRPDUGbDEBQJ6yzcTMzGggiiJrZsYg/0gDd5JXEn9KAXAUNTPAwIGcta1W5Xvv7m2J0eoSK+Jg5ic/+Qk6OzsxZcoUPPbYY3jjjTfw5ptv4tFHH8XUqVPR2dmJ+++/P55rJYopf/Gv9ORXwMzMqNLZ60Sv0wNBwKiaQUODC1UEbHd58NcNx1B7yhrux9KO2+NF5yCjDAYjT9Q+fMoKl8eLuvZe5Xvv7z+VlhO1I+4zU1JSgo0bN+Lmm2/Gfffdp7SAFwQBF110EX7/+9+jpKQkbgslirWG/pmZLCkz08WamVFB3mIqMRpg0KqHuDaNFqGCmd+sOYw/fXwMnxxpx3Pfmp+spcWUPMpAJSDiUQaysypzAQA7GrpxvMMGl0dEpk6NTJ0a7T1ObD7WgfMnF8Vh1fETVdO86upqvPPOO+jq6sKRI0cgiiImT56MvLy8eK2PKC68oogTAceyAaAgW8rMmG1euD0iNGr2JUlnHGMwNvVvnFff3otnPq0DABz3nW4bDdp65FEG+ohHGchmV0uv2YdaLNjd2A0AmFycjTPGmfDiZw14Z2/LgGDmhU31OHzKip99cTpUUf69RBhWm9O8vDzMmzcP8+fPZyBDaanN4oHdJUKjAkpN0rt2Y4YKGhUgAui2MTuT7hpZ/Dsm5WRIb07kYOYX7xyAyyPtJDR1942awbLD6f4rK8kxoNxkgFcEXtlxAgAwucSIZdOk3ZXNxzqCrm+1u/Czt/bj75sbsK9pYOPcVMCe7TQmNXVLdTHleRrlXY1KEJCbxS7Ao0Vjp3RSozKf9TJjSeA20+ZjHViz/5TyGHe4vehK84nan9S2Y8EjH+D2F3cAiL74VzarSkpEbDwqBS6Ti7MxpcQIADjRZQvq0/PpkQ64fV/LJ6BSDYMZGpPMvsxLXmbwQyDfVzfDXjPpTx6kV25iMDOWBAYzaw+2AgBWnjVOyWCk6otxpN7e24RTFgcsdmmbfFrZwK78kZhVlRv09eSSbJTmGKBVC3B5RGVAKwCsO9SqfH4qRQdUjnjQJFE6stqlYMWY0T+Y4Ymm0eKU79hpie8YKo0NgZOz5Xb/MytNOHzKivYeJ5q77Tij3JTMJY6IPFPp1iUTccGUYsysHN6/ZUAwU2yEWiWgMi8Tx9p7cbyjF+NyMyCKItYd8jfEPWVxDHvt8cTMDI1Jlj4pmOk/HTufvWZGDfmdZWkOg5mxxN9nxoVj7VIzuAmF2QMaxaWrFl8wMXd8PubX5EOvGd5JvTPKTdD6DjkYtCqlfYFcYybXnB1ssQZlaVpSNDPDYIbGJDmYMTEzMyrZXR6lAJTBzNgiBzMdPU5lNldNUZYSzDSl+bToFl8wVjbCjKNBq8Y0X4ZqUnG2ckJJPv0nnwaUszJy3VGqbjMxmKExSQ5mBmwzZbNmZjSQU/EGrUo53UJjgxzM1LZa4faKMGhVKMsxoMyXeWjuTt/MjN3lUQqYYxGkz/L1m5lSbFQuk4MZ+Ri7XC9z4WnFABjMEKWUsNtMPM00KgRuMQlC6vXEoPiRgxn5OHZNoZR1KOvXwj8dtfq2mPQaVdBQzeH63gUTcOWscbh58UTlMnkoa2OnDT0ON7Yf7wIAfG1eJYDUrZnhWxYakyy+AuAcQ+jMDLsApzf53WMJt5jGnP4v8hMKpdb95XJmJo2DGSVIN8UmSC8zZeA3V58VdFngNtO2+k64vSIq8zMwd3w+AOmUmN3lSbmu2szM0Jg0VGbG5hTR52RAk67kbaZSnmQac3L6BzNFUjAjZ2ZazHZ4venZOE8uXo5nkC4XAHfZXFiz/xQA4OyaAuQYNDBopefHVNxqYjCTAkRRxPbjXehzDr210Wzuw6ajHUNej8LzeEX02qUns5yM4Hc3GToVMrTSZczOpC+eZBq7DFo1dBr/S1uNLzNTkmOAIABOjxcdvc5kLW9E5CBipMW/g8nWa1CQJfXkeXN3EwBgwYQCCIKgPJ5aUjC7xWAmBby5uwlf+sNG/Pzt/UNe9+a/78A1f9mMfU3mBKxsdOqxeyECEABkGwY+BEy+RnrmPgYz6YrbTGNb4FbThKJsAIBWrUKRb7p0uh7PbjFL9SrxDtLluhmrrzHfghppi6lYDmaYmaFQ3tnbDAB4f1/LoHNDehxu7D7RDQCoPdWTiKWNSnKQkm0QQg5ok49rmzmfKW3J7xzj+Q6WUldgMCNnZgD4TzSlYGYhEokK0gOHs47LzVC2nuQgqjUFi4AZzCSZy+PFp0ekbaP2HicOnbKGve6exm7IsU5Tmr6zSAXWPul/Yv9j2bIcX2bGwsxM2pJPXLD779gkBzOF2bqgwKZcPtGUpsezAwuA4ykwmDl7QoHyeUmOPmgdqYTBTJJtP96FHodb+fqT2vaw193pG9UOAM3dqXdnShfhGubJlMwMg5m05PWKyjtY1syMTXIAM6EwO+jyMlN6Z2bkjGMiMzNnT8hXPpf/LguAaYD1h6Xuijq1dFNsHKS4d2dDt/J5U5q+s0i0Y60uPLnGjA6rv7haPpZtDFEvA/hPOHGbKT119Drh9ooQhOFPFKb0pgQzRVlBl5fnpm8X4MAgPd7bp+EzMwxmKIz1vlbRN5w3HgCw+VgHXJ6BL6KiKGJXQGYmHR+MyfD27l6sPdCHD/f7g79wx7JlJm4zpTX5ibYwWw+tmk9xY9EZ5dIk6QUBWQXAvz1zqMUCs6+TbrpIZJB+enkOCrN1mFudh4o8/9R5+f9fKjbOY9O8JGq12LG/2QJBAG46fwL+b/sJdPY6sauxG/PGBz8IT3T1ob3HfwdK12r8RJOzK83d/q28IYMZbjOlNaXHDLeYxqxvL6zBZWeWD6gtmVIite0/fKoHi/5nLc4oz8GR1h5MLTXi+W/NV+YTpaJEBuk5Bi0++X8XAkBQc77SgNNMoihCEARsP96FA80WzKnOw+llOXFd12D4tiWJPvbVx8wYZ0KRUY9zJ0rpvE+PDKybkbMy431H5rptroj60ox1Pb5+Mi3mgG2mIYIZ+XJmZtJTC49lj3mCIIQskp1SYsRz35qHKSXZMPe5sPFoB1qtDmyobU/5QxXNCQ7SDVr1gC6/ckbI6fai25fZemdvM37y+uf417bGhKwrHAYzSfSfPVJDosVTpQFe500qBBC6bkYOZhZNKUKWTrqDpfqDLxVYffUxLcPJzLBmJi0pxb8m1svQQIunFuPdHy7CH78xG4996UylQVxnijfSS4Ug3aBVIy9Tqkc6ZZXWU9sqtQmRs17JwmAmSVotdnzsK/69ctY4AMCc6jwAwL6T5qB223KHYAA4qzJXmTHCIuCh9TikgMRiF9Hr+9zaF3ouk0yumbHavfAO0veHUhO3mWgoapWAi6eX4avzKpXgINW7Ap9Kkd5J8v8v+fWn1tdOZHJxdtifSQQGM0ny2s6T8IrA3Oo8panThMIs6DUq9Do9qO/oBSCNfP/R/+3BrsZuCAIwtzo/YJQ9i4AHEzi2AACau6WtJrkWRu4n0598yskr+repKH2c9D3JlpoyhrgmEVCQ7cvM9KR2MJOoHjNDkTMwB5qtsNpdyvbXZGZmxh5RFPF/208AAL48p0K5XKNWKQVUnzdZ4PWKuOHZLfj39hNQCcBPL52GqoJMpfETt5kGZ3OKCAxFWsxuiKI45DaTRi0gWy8VvXGrKb2IooiDLdI7xSklyX2nSOkhP022mfY3WQAA4wuyhrhmfE0fJ71G7T1hVraYSnL0A6aVJ1pKBzMPPfQQBEEI+igtLU32skZsj+9OYNCqcMmZZUHfk48U7msyY+9JMzYf64RBq8Jz35qPGxfWAAgYZc/MzKB67MGBSEu3B3aXCLfv4nDbTAC7AKerUxYHOnudUKuEpO/hU3qQg5lU3mbqdbhxsEUKZmZX5yZ1LdPHmQAAnzeZccQ3VmdycfIfayl/NPuMM87ABx98oHytVqsHuXZ6eH3XSQDARWeUIscQHM3Kd5R9Jy3I0kk3zwVTirBoSpFynTJmZiJi7RfMNJvdyhaTXiNArw1/DNOUoUJTl4fHs9PM/mZpAOvEoqwBJzGIQvEXAKde7xTZ7hPd8IrSOIayJG+fnlEuvUad6OrDlvpOAMDkFMiCpnwwo9FoRkU2JtC+k1KEvcR3iinQ9HJ/1GtzSidwFve7HguAI9O/3qWl26MU/xozBu8n4e8C7AHAF8V0IT+25CdcoqHkZ0mn3hKxzeT1ivjkSDvmjc9Hhi7y5xW5+/ss3yGRZDJlaFFdkInjHTas/rwFQGpkZlJ6mwkAamtrUV5ejpqaGnzta1/DsWPHBr2+w+GAxWIJ+kg1db7i3sBprrIppdnQqAR021zY4bsDB2ZlAH9mptlsH3TK9lgnZ2ayDVLg0mL2DDmXSWZir5m0tL9ZerxPS2LzLkovidxmevqTOnzzmS144sPDUf3cDt9p1tlVyQ9mAP8OgjxXMBXq01I6mFmwYAFeeOEFvPfee/jLX/6ClpYWnHvuuejoCD+/aNWqVTCZTMpHZWVlAlc8NKvdhTarlM4cHyKY0WvUQVXhk4uzMS43OK0oZ2ZsTg8sfW5QaHLNzKRiaSuv2+bF4VNSo6dwxb8yuWaG20zpRQlmyhnMUGSU00wJCGbe3C31Ftt4JPxrWH+iKCpDhmdX5cZhVdGb3i/zyczMEFasWIEvfelLmDFjBpYuXYq3334bAPD888+H/Zn77rsPZrNZ+WhsTG5Xwv6Od9gASPu04aq/zwh4Il48tWjA9w1atfJu4iS3msKSMzPFJrUSvLy2TcqKzZsw+PFGZmbSj8XuUh5fyWyrTulFOc0U56PZTd192HtSquk62GKBwx1ZB/f6Dhs6e53QaVQps306Y5x/HcVGPUyZyT3JBKR4MNNfVlYWZsyYgdra2rDX0ev1yMnJCfpIJXXt4beYZNMDgpkLpgysqwECt5oYzIQj18wYDSqUmaT9aY8XGF+owdIzBi+iYxfg9HOwWTqSXWYyKC9QREORC4CtDnfEAUakunqd+PRIO0RRxAcHTimXuzwiDvjur0ORt5hmjDNBp0mNl2z5eDaQ/M6/stT4PxMhh8OBAwcOoKysbOgrpyg5mAm1xSSbUSFFvRlaNebVhN4jlSvaOT07PKVmRq9Caa6/2O7bF+RAPcRAOc5nSj/7m6R3vWdwi4mikGPQKs8HXb2xnaT90zc+x7V//Qy/+aAWa/ZLwYw8t3HPie6IfseOBrleJjemaxuJ3EydMk17UpI7/8pSOpi55557sH79etTV1eGzzz7Dl7/8ZVgsFlx//fXJXtqw1UeQmZldlYd7lk/Bb66eCb0mdMW7nJk5xWAmLLlmxpihwqQSKQ26aKoB08YN/a7dxJqZtMPiXxoOlUpAXqZcBBy749miKGLzMeno8u8+qlVm7l12ZjkAYHejOaLfIf/cnBQ4yRTonAnSYORUWVdKH80+ceIErrnmGrS3t6OoqAhnn302Nm/ejOrq6mQvbdiORRDMCIKA2y6cPOjvkVtayy2uaaAeh7TNlK0XsHBKJsblaXBGBIEM4N9mstpFeLw8MZYO5LQ962UoWgVZOrT3OGJaBHzK4kB7jxQciSLgEUVMKs7GFTPL8dbupogyM/uaLKhr74VBq8L5kwfWTybTTy6bhstmlmPR5MJkLwVAigczL730UrKXEHP1gxzLjoY8RO8Ug5mwlJ4yBhW0agEzqyKfopxtUEEAIEKqvVGl9COFAOC477E1MUXS3pQ+4jHSQC72nVCUBZ1ahYMtVlx0RgnOrJTKCI609aDH4Ua2PvyTy1t7pNNPXzitBFmDXC8ZTBlaXDAldQKslN5mGm26ep3otkl7siOdr1Ea0GuGQpMnZmcPMrYgHLVKgNHXn8bcx8xMqrPaXbDYpTYF5bkcMEnRyfcdz+6I4Ymmvb7My+yqPLxw43z89LJpuHnxJBQbDSgzGSCKwOcnw281eb0i/rO7GQBw+cz0rRNNFAYzCSQ3yyszGaLq/hiKPIadNTOhebwieh3yaabBi33DycmUbqNUKQL+3w9rcfETHyt9isivyTenzJShHfSdLlEoBXHMzMwYZ0JxjgHfXlij3DfP9B3yGGyraWdjF0529yFbrxnQBZ4GYjCTQHLxbyymnsqZGavDrXRhJD85kAGGl5kB/EFQ4O9Kpn9vb8TBFive9qWeR4OT3X3Y5WsINrLfI/WX6d9gkigSse4CLIoi9vpGa0wfN7A3zJkVuQCAbfVdYX/HW76szPJpJZwzFgEGMwmk9JgpGnkwk63XwOiL8lsCsjNer4gNtW0pP84+3uRj2Zk6Ychj2OFk6qSHhy0FghlRFHHKImVkPomie2gqe29fC5Y9vh5X/v5T7Gsa+mTHYE52Sf2WuMVEwxHrYZNy8a9KCH26Tm6G+uHB1pCNT71eEe/slYKZy7jFFBEGMwmkBDMxyMwA/uxMYBHwf799ANc9vQUPvbkvJn8jXfXYh18vI8vSS0FQKgQz5j4XnG7p37T5WAdcntTY+hquv244hu/9bTtsTg9EEfjnloYR/b6Tvm0mufcFUTRiPWxS3mKaUmIMWVJwRrkJ504sgMcr4plP6uBwe7Dq3QP4x2fHAUinmFqtDmTq1DhvUmqcFkp1DGYSqNH37rGqIDMmv69/EfA/PjuOZz6tAwAcPhVZd8nRSukxM8x6GcAfCPU6kx/MyFkZQBrutjsGWzPJ0uNwY9W7BwH4h6i+sbMJfc7hd1+V391ym4mGoyDGBcBy8W+oLSbZTYsmAABe2tKAW/6+A39afww/ff1znLLY8dHBVgDAwkmFYXuNUTAGMwnU5sugyMW7I1UScDx7R0MXHnjDn40Z7Mh2s7lvxGn9VGf1jTIYSWYmU586NTP9b89PjrQnaSUjt7/JAo9XRJnJgOdumIfK/AxYHW687UurD8fJLqlmhttMNBwFMa6ZCSz+DWfxlCJMLTGi1+nBh77gxSsCr+w4gY8OSV9feBoLfyPFYCZBRFFEm6+BUrEx8n4ngwmcz/R/20/A4xVx3iSpK2OXzQW7K/Q73W89uxVXPPkpjrb1xGQdqSgm20wpVDMzIJipDQ5mGjtt8KZJcz/5OOoZ5SaoVAK+Nq8KAPDy1uFvNcmnmcZxm4mGQS4ANve5YrKFKzdwHGy0hiAI+M75NQAAjUrAlbPGAQD+tum4csppCYOZiDGYSRDpQSK92MgpzZGSMzMtZge21Elts795znhlGFmrZWAxm9nmwsEWK9xeEesPtcVkHalILgA26od/F882pE7NTKvvOLbcQnxnYzesdqln0erPm3H+Y2vx2w/DD2BNJZ/7soLysLovz6mASgC21nehdhjbo063F6esvmCGmRkahtxMnTIzqcs2suyMxe5SOrNPHmII41WzK3DfitPwwrfn479XTkemTo1msx2iKBUOxyqLPxYwmEkQuTdIbqY2ZnugcmbmQLMFR1qlLMv88flKd+BQow4OtFiUzzceTd+tiqHIE7OzM4ZfM5PlC4RSaZtp7vg8VOVnwuMVlQB27UEpKH1xS0NajF7Y3+Q7sloupeBLcgxYenoJAOBPHx+L+ve1+J78dRqVsl1AFA11wHymkRYBy8/FJTl6mDK0Q/7d710wEedOLESWXoNLZ/hPLnGLKToMZhJEDmaKsmOzxQT4MzNy8eOUkmzkZekGDWYONvuDmc+OdcKd5qdiwolFZiZT56uZSYkCYOm2LM4x4OwJ+QD803TlAYttVgc+q4vu2Pa7e5tx0W8+xsGAIDee7C4Pan1P9oHFkbcsmQQAeG3nSTR22qL6nYHFv6phHsMnkp+bR9qU8sgp6f49uXjwrEwoX51XqXzOLaboMJhJEHmboChG9TKAPzMjmzdeepErGWSi9sEWfxrf6nDj86bEvIglmjUmR7NTqWZGuv+UGPVKw609J8xwe7w4FLA1IzfaitQrO07g0Ckr3t4z/OLbaBxsscLjFVGQpUNJjv+xcFZlLs6fXAiPV8Qf1x+N6nfyJBPFQrHv/ngqxPZ8NGpbpcfj5JLoZ4TNrc7DVbPG4dIZZTirMndE6xhrGMwkSFscgpn8LB10av9NOL9GCmZKfQ/K0NtM0gMt09f7YLRuNbVapOLnQmMMamZSIDPTGnASbmZAMHO0rVfpPwNI9TPRFDDKT9y1pxJTDK4U/44zQRCCsyi3+ybF/3vbiaBGkEORG+YxmKGRKInR8F458ziczIwgCHj86rPw1LWzh93sc6xiMJMgsT7JBEh3/OKAd7dyMFMSZpvJ4xVxyLed8JU5FQCATUdHRzfZQB6viDarFMyUmoY/p0fOzDjdUIq3k8HrFZXMXkmOAVNLjdCpVTD3ufDevhYAUmajIEuHLpsLn0ZxbFt+4j6SoJNtckuA6SFOecyvycf88flwerx4ZccJ5fIjrT1wuMP3oGnqZvdfGjn5ubl1pMGMvM00jMwMDR+DmQSJR2YG8G81VeZnoMwkPZmXhtlmOt7RC7vLC4NWhWsWSMdht9Z3DvpCkY7arB54vIBWDeRnD/8unqETIL836kvibMeOXifcXhGCABRm66DTqHB6mfSu7+WtjQCkwXUrZpQCAP4T4ZaR2+NFuy/Irm/vTUhX4c8HmVcDAJefVQ4A+PiwVNT8nz1NWPr4evzP6kNhf6eyzcRj2TQC8pvA1ghrZnodbnzh1+tw3dOfKZnEHodbuT9OKmIwk0gMZhIkXsGM/ACU62UAhC0AlutlppYYMbXEiMJsPewuL3Y3jq4GeqfMUnBWnKOGShh+qlYlCErjvBGe1hwROXtSmK2HxretKNfNyE+cp5fl4BLfSYi1B1sj6jnT3uOEfDW3V8Rx31R3WavVjr0nzNjR0IXeGAwzdXm8OOS7D8onmfpbNFlq3b6joQs9DrcSrK32ZaBCYc0MxUKJUjMTWWZmd2M3jrb1YkNtOy753w3YUNuGo74tpsJsPfJ4si6hGMwkSKuvD0ZRdmz7Blx2ZjmKjXp8fX6VcpnyDsPigCj6X9Tkk0ynleZAEARM86X65Wne6ajN6kFXb3BmqcU88i0mmbzV1JfEYEa+7wQWzJ5ZERwMTCvLwdzqfGTp1OjodWJfBIXd/YNd+UgpANSesuKcVR/h8ic/wVW/34iv//WzkfwTAACHWqxwerwwGjSozA8deFQXZKEqPxMuj4jVn7dgo28b9ERXH050DTzl1G1zKqefwv1OokgUKzUzkWVmGgPuj529Tnz7+W344MApANLJUkosBjMJEq/MzMXTS7Hl/qWYG5CZkYMZp8cb1DNhv68rpbxFUd5vtlOq2XrMjkf/04WGDlfI7/c5vbj7xXbc9Y92dNv8AU2LWcoilOaOvJ9Plk5unDfiXzVs/pNM/kB4ZsBJB5UAqY5Go8K5vqF06w+3RvB7wwcz2493weMVYdBKTxG7G7thc44sO7PTN0/qrMrcAcW/gc73ZWceXX0wqG/OZ8c6B1z3rT3NcHtFnF6Wg4q82Mw8o7FJqZmx2oPeBIbT2CllBL8ypwLnTCiA0+3FU2uPAAAmFzOYSTQGMwngdHvRZZNekGMdzIQS2Dws8N233EvkNN9IernGptk8cAR9soiiCJvDi39utuKX/+nGlmMO/HVd6CzDyS4Peh0iLHYRz23wH08+pWRmRh7MZMrHs5N4oimwx4xsYlG2ciJtQlE2DFrp88VTpcGN6yLo7tw/mKkNCGYafNmOr8ypRKGv/8ZITzzt9PXFmTXEkdPzJ0v/BvkNQI5ByrBtPjawWP3/tkuFwl+aPW5EayOSn5tdHlF5vgakAvyOHseAAEd+jEwuycYvrpwOnVqlbNtOGqLzL8Ueg5kE6OiVnpS1agG5Q3SEjJX+xwz7nB6c8B1hnep7oMnFw00pkpn552Yrrn7qFK77Uyv+b4u09SUA2HfShSOnBmZnmrv9mYINh+zY3SD9f5a3mUpisM0kH89O5jaTkpkJ2GZSqwSliPb0Mv/JoAt8U6h3NHTBbAud0fL/Xul2lzN0R0IEM1X5mUrK/NAIJ7HvkjMzVbmDXu/cSQVBx1LvWDoFALC5X0PAI6092N3YDbVKwBVnMZihkdFr1MqMplMWOzxeEa9sP4GLf/sx5vz3B/j39hNB15e3mSrzMjGhKBvfv2CC8j1mZhKPwUwCyO8wC7P1CetQKp9oajFLf/tkt/TAMxo0SmFaWa5vm6k7+ZmZPqcXb+7ohXygJscg4PZlJiw6TVrj69sH1vU0d0tBi28UFf681gKXR4xpZkZpnJfMmhlf0FHab07L0tOlDqFyAAMAFXmZmFScDa849GRt+b5xnm9r6mhbj1I43OgLfCvzMzHFF/wOZ25SR48Dbo8XZpsLx9qk2/CsyrxBfybHoFWyN6eVGvHVeZVQqwQ0dvYpxb4AlOPbi6cUJSTjSaOff6vJgSc/OoK7/70bh30Zyff3nQq6rr9WS9revGXJJJxWakRhtm7QAZMUHwxmEkAe+JjIJ9z+vWYaQzQWk7eZImlQZnd54np0d1udA063FIC8eHMJnv1uCRafnoEvzs4CAHx21K7UwsjkzMwXZ2chJ0OFFrMH6w/2we4SIQAoNsZgmykVamaUAuDgYObbCyfgw7svGLDFstgX3AxVNyMXFs8bnw+dWgW7y6sEC41BmRkpmDkU5TbTJ7XtmP/Ih7jv1b3Y5ZsCXF2Qqbz7HcwVvgnC18yvQrZeo2ShPvNtNXm9Il7bcRIA8CVfzySikQrMaL/7udTiYIlv61bukQQANqcb7T3SOxw5mDFo1Xj91vPwyf+7EEZDYjLw5MdgZgQaO20RFYrJDfNiOZdpKPK7eLnXjLzFFFgkKW8zWR1uZQJzKK1WO+b94gN87c+b4xbQbDgkrXPhVAP0Wn/2anyhFmdV6eAVgXd3B59mkbeTJhRpcf5U6d/y6lbp3X+BUQWtZuRZsGSdZnJ7vPjF2/tx3i8/Unqz9A+G1SoBE4uyBxTTLvIFM5tDFMwGkoPY8twM1BRKQeOR1h5Y7S6lcLwyPwNTS6WU+eGWyDMzoijif94/BI9XxP/tOIF/bZOOWEfaov0bC6rwyf9bgm+eUw0AyjwquW5mf7MFLRY7svUafOF0zrCh2JAzM3XtvTjsy0T+5LJpEATpoITcl0l+Ps0xaIKGSRq0aqV+jRKLwcww/XtbI5Y+vh5//6xhyOvG6yTTYEpNwSMNTirBjD8zk6XXKMWVg51o+vhwO6x2N7Yf78LTG4b+90bL2ufFLl+9y/lTBh6vXTZdCsB2HA9Oj8iZmdJcNc6fIr+jkgKcshjUywABIw0SmJkRRRH3v/Y5/rKhDie7+yAIwHmTCjC1NLKiQvnIfWOXDXZX+IaIp5QRCXpM8tXF1LZalVMa+Vk6GA1aTPZlZlosdpj7Bq/DkX1ypB27fTUyoghl9tNQxb8yQRBQkZepBGrnTCgAAHx6pAOiKOIz38TwuePzYjaFnkjOzHx0oBVeUXpTOLEoWwn25XEcDR3BW0yUfAxmhsnc54LD7cXP39qPPb4UejjJCWaCTyrJPToq+nVJlVvADxbMbAkovPz9uno0dcX2ZM+mo3Z4vMD4Qg0q8gcGIdMrdFAJQFOXBx090otzj90Lq11aR6lJjUkl2qAamZIY1MsAQKZOzszE/zST1yviSGsPHnpzH17e1giVADz25TOx96GL8I/vnA2tOrKHa0GWDjkGDUQRqO8I3UOoz+mBxS4FgyUmA07zBSx7TpiV4t9K330lx6BVsniR1s387iPpiKp8ukp2VtXg9TLhLKgpgE6jwsnuPhxp7VEyNGf7ghyiWJCL7OVidzmTKDd5lPs3ycW/VQxmUgaDmWH69sIaLJtWAqfHi1v+sWPQkyNyMBPLuUxDmeB7J1HXLg0iDLXNBPgLhQcrAt7iexdcbjLA7RHxtw3eiLbXIvXpYSmQkreK+ss2qDChWErl7m2Utj/kLabcTBUydCoIghD087Eo/gWALF8H4N44bzM1dtpw4a/XYenj6/H8puMAgEeunIGvzq1Etj66LJMgCJjoO01xtDV0MCNnZTK0ahj1GqVP0Za6zgGFjQCUupnDEdTNbK3vxJa6TujUKvzyqjOxYro0ZkGn9o9hiFaGTq0ELh8ebFXukwtq8gf7MaKoFPerS5P7Oc3w1WwpmZkQjxFKLgYzwyQIAn715ZmozM/Aia4+/HL1wbDXVbr/JjCYqcjLgFGvgcsj4lh7T0AwE5yZkYuAwx3PPmWxo77DBkEAnr9xPjK0KtS1AkdOxWaek8Mt4kCTFCmcPTF8d+QZFVLR6N5GKTCUt5jKAhrjnT/V/2+LRfdfQAqkgPjPZvrF2wdQ32GDQavC/Jp8PP7VmfhaQFfnaE30zYU5FmaApBzMlJoMEAQBs6pyoVOr0Gp1YIPvFFRVUDAj/b7Pm8x44I3P8b2/bQu75STPVbr0zDKUmgy4c9kUGPUaLJ1WPKItIbkQ8/mN9TD3uZClU4ed8UQ0HP3fcM6slO5fZ4yTtm4/9xUBy1uxDGZSB4OZETBlavHgZWcACN3QC5Cq3ut9+6uJDGYEQcBpvnfBuxq6lcK1AdtMyhHu0JkZ+R3wtLIcTC4x4qLpUrHlugOxSVXUtbrg8UoZlsG2hmZUSsHMnkYnRFFUjmWX5fqDlnF5Gsyq1kGvETClNDanCeTMTDwLgDcf68DqfS1QqwS8ddtC/Ot75+Cq2SM7oTOhSMrMHQ0TzMi1VPKTt0GrVp64N9RKwUhViMzMi5814IVNx/HevlP47gvbQg4plQuLJ/rWMKXEiM0//gKevGb2iP5Ni6dK9z15S3Tu+PyIt96IIhF4YlAQ/BmZM3zbTI2dfcEjNDjcNGXwmWCE5DRkfUdvyHbv//vhEXT2OjEuN0N5QCTKaaXSu4kPDkhHdLP1wZX3QMA2U0BmxusV8fHhNvQ43EowM9+Xzv/SHGmY4aZaJ/qcIz/ZdLhFenc/pVQ7aIv708p10KqBzl4vTnZ5/JmZfgHQjy7Nw1++XYSCGBzLBoAsX81MrwMx3VqTeb0i/vvt/QCAa+ZXKsW2IyVnZo62hd5mktsFyLc/INWlAFLBLhAczAQWH+s1KmTrNfisrhN3/2v3gP8vLUrWJ7jYfKQ9lmoKs1Bd4F/TggncYqLYCnzDObk4WzlibcrQKo+Hz09a/A3zmJlJGQxmRqjIqEdBlg6iOLDd+6EWK/664RgA4OEvnpHwI3tyZuaTI9I77Yq8jAEBg1wA3BRQM/PCpnp885ktWPnUp1jv2zKQaxNmV5lQYgIcbuDT2pF3Dj7cIqU8hsqk6DUCTivzbzUpwyRzg7eTdBpBOU4dC/LUbK8I9Llifyz9/f0t+PykBUa9Rul0GwuB20yhgrAW5SRTQDDTLzgIfKKeXGxEbqYWmTo1nvvWfPzxG3OgUQn4z55m7GjoDvq5U2Ga/MXC4oAGgSz+pVjTqlUozJaeZ2b6JtPL5CzN+/tbYHN6IAic1J5KGMzEgBw0HOrXh+OBNz6H2yti+bQSLJ1Wkvh1+TIzdt+LcP8tJsDfa6bZLA1X83hFPPNpPQCp54hc6DbPVyAqCALOmyq9wH+4b+Sdg/2ZmaEbqclbTR/s68PJroE1M/Fg0AqQdzKs9pENWgxlzX4pa3b1PP8MpFioys+EWiWg1+kJOQX4VIhgZnZVnjJGQK0SlPsGIBXgrv7hIqy7ZzHOmViAhZMLcc5EKZio6zd1Xd5mktsDxJK81ZSpUysvLkSxVOwb6DqzXxsBuT7rBV+BfonRwJ4yKYTBTAxMLZGChgMt/oGIFrtL6YXxwOXTkrKu0/r1JQn1LkIuALb5juquPdiKhk4bTBla5UTUpOJsFAS80J49SXqBP9ziCjvROhIdVg86erxQCcDEkqELdhdOMSBTJ6C+3Y1eh/9YdjwJgqB0Abb0BQczbVYH/t//7VGyV9ESRVGpT5FfpGNFp1Gh2pdZCVU3I3f6DZz3lKXXKAHCuNwMaPrVo5SaDEGnPcp99x25hxEg1YgpR77jkJlZNKUI379gIlZdNYP1MhQXN5w3HmdPyMclM8qCLv/ynApcdEaJMkk+0gaQlBixOfIxxoXKzNT73q0WZusHHIdOlCy9BtUFmTjeIfeYGbiODJ0auZladNtcaDb34dmNdQCAr82vxHfPn4D//bAWy6aVBv1MTqaAWdVabKtz4cN9ffjWouEV2x72DY+sLtAoTxCDKTFp8MurC/DYf7pwosuD/CzpWHa8ZeoFWO0iLAFdku0uD256YRt2NXZjS30n1t6zOOrfe+iUFa1WBwxaFeaOH17/lcFMKMrCsfZeHG3rUeYvAVJ34QPNvgnqpcEzZBbU5GNXY3dE/TNCbVHKWZksnTouLd3VKgH/teK0mP9eItlX51biq3MrB1xeZNTjT9fNhd3lwf5mCyZxmGRK4VubGJAzIAdbrEp9gpx6rylMboFYYHYm1DYT4M/OvLDpOD490gGVAFx3djUKsvV4+IrpWDi5cMDPLD5d2vL5+GAfXO7hFcbK9TKTyyJ/0RuXJwU0V83NwneXJGaYm3yiSd5mEkUR/++VPcoU6DpfwBCtDYelI9ALagrikq72180EbwPVtvbA7vIiW69Rsm+yry+owszKXHzj7KGPhY/Lk4/1BwQz8vaVKfZZGaJUYNCqMbsqDzmcv5RSGMzEwORiIwQB6Ox1KnOY5GzI+IKswX407gLfeYfLEMnHs1/0jWZYPq10yGzSmVUa5GepYLGL2Fo3vCYstVHUywTK0Klw7blGzJuQmBfM/ttMb+xqwhu7mqBWCcrpmg8PnAr78+F87NtiWjSlaIhrDo//RFNwoLX3hNQrY/q4nAEnjKoLsvDGrefh4unBKfZQyn1T1wMnWcez+JeIKBwGMzGQoVOjxhe0HGyWtprkbabxhckNZgI7ro4Lk5m5cWEN5o3Pw/mTC3HV7HH48SWnD/l71SoBS6ZJv+/DfbYhrh2s3erBn9aacajZfyw7lRXlSA+TfU3SbfvOXmnO0PcWTcCN59UA8B9/j5Td5VGOvS8KkfmKBaXXTGtwMLPbN36j/2mNaI0L2GaSM5ItZt+RbwYzRJRArJmJkamlRhxr78WhFisWTSlCnW8mTrIzMzMqciEIQF6mDnmZoYOG8yYVBtVUROrCaRl4ZWsvdjc40WrxoDhH2irpc3rD1rJY+7y496UOWPqkE1aLphpQHucTSSM1s0qDtfudWH+oA063FxuPSg0SL55eivwsHR58cx+21Xeiq9eJvKzgLJPd5cGeE2bUt/eivkP6aLc6odOo4HB7UZpjiNve++RiKZBtMtvRZnUoPTT2+lqynznCYEbuUWN3edHZ60RBtt5/SorbTESUQAxmYuS00hy8+3mLcqLJn5lJbs3MuNwM/PWbc5GfpRu0Kd1wlJo0mFGhw94TTqzdb8PVZxvx+vYe/O3THvzXZbkht4G21tlh6fOiyKjC7ctzcca46LaYkmF6hRZqFdDQ2YdXdpxAj8ONgiwdppeboFIJOK3UiIMtVvzjs+OwOT3QqlX48pwKHGqx4ievf67UkYSyeGpRzG8XmSlTq6xta30nLplRBofboxT/nlkxsqPNeo0aRUY92qwONHXbUZCt9x/LZmaGiBKIwUyMyB1S9zdZYLa50GWTT+okNzMDAF84PX49br5wRgb2nnDiowN9uPSsLLyyVQriNhyyhwxmtvnqa5acnpEWgQwAZOgETC4FDjYBv3rvEACpzkWuN1l6egkOtljxq/cPKz/z2w9rlc8Ls/U4vcyI8QVZGF+YhSKjHq0WOyx9Lnzj7Oq4rn1BTT4OtlixpU4KZg42W+HyiMjL1IYtCI9GeW4G2qwOnOzuw4wKU8hmfERE8cZgJkZmV+cCkE40bamXaiGKjPqoJx6nmwUTDcjWW9Bu9eLx1d2wOaXaiX0npRlKgVkHl1vE7gbpBNOcmvR6sZtRJeBgk4gO3/jsxVP9RbuXzCjD79cdgVcELphSBK8oYkNtO1QCcNOiCbhz6ZSkNdeaX1OA5zcdV3oe7fFtMUnbjyPPCFXkZmB3Y7dyPDtwgCURUaKM7lfaBCo2GjCzwoTdJ8x49lOpV0tNCmRl4k2nEXD+1Ay8u8emBCoA0G2TZihV5PvvYp+fdMLuEpGXpcKE4vS6602vFPDvzVKgJgjA+ZP9wcy08hy8cetCZOrVygmiE102eLxi0jNz82qk/jUHW6SM4R7fcfKZI9xiksknmpq6++Dximi1sgCYiBKPp5li6MLTpO0cuUA02fUyibL0DP92RX6WCqeXS4XGn58IHjUtbzHNGa+HKk51IvFSYhJQXSD9O8+syEV+v0LfGRUmJZABpGPwyQ5kACnInlCUBVEENh5tx+Y66b4Zq1EAcuO8k9196OhxwOMVoRKgzLchIkoEBjMx9IXTg1vSp8KLWSKML9JiUokUwFw+Kwszq6RTM/tO+oMZURSxvU7agphbE/uZPYmwYrp0+146o3SIa6YWeUjoj1/bi8bOPuRn6bAgRkMaA7sAy/UyRUb9gFEIRETxlF65/hR3RnkOSnL0ymC/miT3mEmkOy8yYe8JJy6clqH0j9l3wl8309DhRpvVC50aOLMyPYOZm5eMx5KpZZg7Pn/oK6eQ+TX5+OeWRqUo/dEvnQlTRmx6+4xTMjP2gAGTnCRMRInFt08xJAiCstUEJL/HTCKV5mqwbLo0qXlyiRY6DWDu8+JEpwcAsPmI9EJ3ZpUeem16bTHJtGoVFkwoUCZLp4v5Nf4szDfOrsKyGE5wl4OZ9h6HMmG9NCc9g1UiSl8MZmLsC6f5t5rkVvdjjVYjYGqZVDPx+QkpS7X5qPTfsyexMDTRxuVm4Jr5lVh6egnuvyS2E9xzM7XI8J3UWntI6oLM4l8iSjRuM8XYwsmFWFCTj3G5Gcga5ceyB3NmhQ57G534cH8fZlbp0dDhhloFzEvTepl0t+qqM+PyewVBQHmuAUfbevHpEam4eGZlblz+FhFROGP31TZODFo1Xv7eOcleRtItnZ6JV7b1oq7Njd9/6BtsWKFDtoHJwNFmfEEWjrb1wqjX4I5lU7DyrHHJXhIRjTEMZigucjJUuPSsTLyytRcHmqTCU24xjU73XXIaZlfn4StzK1Bs5G1MRInHYIbi5vJZWXh3tw02pwgBwIIJ3GIajSYVGzGp2Dj0FYmI4oQ5f4obo0GFy2dJJ7rOqNDBlJna07GJiCg9pUUw8/vf/x41NTUwGAyYM2cONmzYkOwlUYSumpuF7y3JwS1fyEn2UoiIaJRK+WDm5Zdfxh133IH7778fO3fuxPnnn48VK1agoaEh2UujCGjUApbPyESJiTuaREQUHykfzDz++OP49re/je985zs4/fTT8cQTT6CyshJ/+MMfkr00IiIiSgEpHcw4nU5s374dy5cvD7p8+fLl2LhxY8ifcTgcsFgsQR9EREQ0eqV0MNPe3g6Px4OSkuD26yUlJWhpaQn5M6tWrYLJZFI+KisrE7FUIiIiSpKUDmZkghA8C0ceXhjKfffdB7PZrHw0NjYmYolERESUJCldlVlYWAi1Wj0gC9Pa2jogWyPT6/XQ69nPhIiIaKxI6cyMTqfDnDlzsGbNmqDL16xZg3PPPTdJqyIiIqJUktKZGQC46667cN1112Hu3Lk455xz8Oc//xkNDQ34/ve/n+ylERERUQpI+WDm6quvRkdHB372s5+hubkZ06dPxzvvvIPq6upkL42IiIhSQMoHMwBwyy234JZbbkn2MoiIiCgFpXTNDBEREdFQGMwQERFRWmMwQ0RERGmNwQwRERGltbQoAB4JURQBgDOaYsTmsqHX2guXygW9Zmw0J3S4HXB6nbBYLHBr3cleDqUKtw2w9gIqF6AeG4+FmPA4AK8TMFgADR9PFJ78ui2/jg9m1AczVqsVADijiYiIKA1ZrVaYTKZBryOIkYQ8aczr9aKpqQlGozHsPKdoWCwWVFZWorGxETk5OTFYIcUDb6f0wNspPfB2Sg+j7XYSRRFWqxXl5eVQqQavihn1mRmVSoWKioqY/96cnJxRcWcZ7Xg7pQfeTumBt1N6GE2301AZGRkLgImIiCitMZghIiKitMZgJkp6vR4PPvgg9HqeXkhlvJ3SA2+n9MDbKT2M5dtp1BcAExER0ejGzAwRERGlNQYzRERElNYYzBAREVFaYzBDREREaY3BTJR+//vfo6amBgaDAXPmzMGGDRuSvaQx66GHHoIgCEEfpaWlyvdFUcRDDz2E8vJyZGRkYPHixdi3b18SVzw2fPzxx7j88stRXl4OQRDw+uuvB30/ktvF4XDg9ttvR2FhIbKysvDFL34RJ06cSOC/YvQb6na64YYbBjy+zj777KDr8HaKv1WrVmHevHkwGo0oLi7GypUrcejQoaDr8DHFYCYqL7/8Mu644w7cf//92LlzJ84//3ysWLECDQ0NyV7amHXGGWegublZ+di7d6/yvcceewyPP/44nnzySWzduhWlpaVYtmyZMq+L4qO3txczZ87Ek08+GfL7kdwud9xxB1577TW89NJL+OSTT9DT04PLLrsMHo8nUf+MUW+o2wkALr744qDH1zvvvBP0fd5O8bd+/Xrceuut2Lx5M9asWQO3243ly5ejt7dXuQ4fUwBEitj8+fPF73//+0GXnXbaaeJ//dd/JWlFY9uDDz4ozpw5M+T3vF6vWFpaKv7yl79ULrPb7aLJZBL/+Mc/JmiFBEB87bXXlK8juV26u7tFrVYrvvTSS8p1Tp48KapUKnH16tUJW/tY0v92EkVRvP7668Urrrgi7M/wdkqO1tZWEYC4fv16URT5mJIxMxMhp9OJ7du3Y/ny5UGXL1++HBs3bkzSqqi2thbl5eWoqanB1772NRw7dgwAUFdXh5aWlqDbS6/X44ILLuDtlUSR3C7bt2+Hy+UKuk55eTmmT5/O2y7B1q1bh+LiYkyZMgU33XQTWltble/xdkoOs9kMAMjPzwfAx5SMwUyE2tvb4fF4UFJSEnR5SUkJWlpakrSqsW3BggV44YUX8N577+Evf/kLWlpacO6556Kjo0O5TXh7pZZIbpeWlhbodDrk5eWFvQ7F34oVK/CPf/wDH330EX79619j69atuPDCC+FwOADwdkoGURRx1113YeHChZg+fToAPqZko35qdqwJghD0tSiKAy6jxFixYoXy+YwZM3DOOedg4sSJeP7555VCRd5eqWk4twtvu8S6+uqrlc+nT5+OuXPnorq6Gm+//TauuuqqsD/H2yl+brvtNuzZsweffPLJgO+N9ccUMzMRKiwshFqtHhDFtra2DoiIKTmysrIwY8YM1NbWKqeaeHullkhul9LSUjidTnR1dYW9DiVeWVkZqqurUVtbC4C3U6LdfvvtePPNN7F27VpUVFQol/MxJWEwEyGdToc5c+ZgzZo1QZevWbMG5557bpJWRYEcDgcOHDiAsrIy1NTUoLS0NOj2cjqdWL9+PW+vJIrkdpkzZw60Wm3QdZqbm/H555/ztkuijo4ONDY2oqysDABvp0QRRRG33XYbXn31VXz00UeoqakJ+j4fUz5JKz1OQy+99JKo1WrFp59+Wty/f794xx13iFlZWWJ9fX2ylzYm3X333eK6devEY8eOiZs3bxYvu+wy0Wg0KrfHL3/5S9FkMomvvvqquHfvXvGaa64Ry8rKRIvFkuSVj25Wq1XcuXOnuHPnThGA+Pjjj4s7d+4Ujx8/LopiZLfL97//fbGiokL84IMPxB07dogXXnihOHPmTNHtdifrnzXqDHY7Wa1W8e677xY3btwo1tXViWvXrhXPOecccdy4cbydEuzmm28WTSaTuG7dOrG5uVn5sNlsynX4mBJFBjNReuqpp8Tq6mpRp9OJs2fPVo7HUeJdffXVYllZmajVasXy8nLxqquuEvft26d83+v1ig8++KBYWloq6vV6cdGiReLevXuTuOKxYe3atSKAAR/XX3+9KIqR3S59fX3ibbfdJubn54sZGRniZZddJjY0NCThXzN6DXY72Ww2cfny5WJRUZGo1WrFqqoq8frrrx9wG/B2ir9QtxEA8dlnn1Wuw8eUKAqiKIqJzgYRERERxQprZoiIiCitMZghIiKitMZghoiIiNIagxkiIiJKawxmiIiIKK0xmCEiIqK0xmCGiIiI0hqDGSIiIkprDGaIiIgorTGYISIiorTGYIaIiIjSGoMZIiIiSmv/H+rkhbM0HYPVAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(resids1, dist)\n", "plt.ylabel('Ca distance (Angstrom)')\n", "plt.axvspan(122, 159, zorder=0, alpha=0.2, color='orange', label='LID')\n", "plt.axvspan(30, 59, zorder=0, alpha=0.2, color='green', label='NMP')\n", "plt.legend()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Calculating the distance with periodic boundary conditions" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "It is common to want to calculate distances with the minimum image convention. To do this, you *must* pass the unitcell dimensions of the system to the ``box`` keyword, **even if your Universe has dimensions defined**.\n", "\n", "This should have the format: ``[lx, ly, lz, alpha, beta, gamma]``, where the first three numbers are the box lengths along each axis and the last three are the angles of the box." ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "execution": { "iopub.execute_input": "2021-05-19T05:57:17.056695Z", "iopub.status.busy": "2021-05-19T05:57:17.055913Z", "iopub.status.idle": "2021-05-19T05:57:17.058843Z", "shell.execute_reply": "2021-05-19T05:57:17.059530Z" } }, "outputs": [], "source": [ "resids1_box, resids2_box, dist_box = distances.dist(ca1, ca2, \n", " box=[10, 10, 10, 90, 90, 90])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Plotting" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "execution": { "iopub.execute_input": "2021-05-19T05:57:17.070129Z", "iopub.status.busy": "2021-05-19T05:57:17.068098Z", "iopub.status.idle": "2021-05-19T05:57:17.201460Z", "shell.execute_reply": "2021-05-19T05:57:17.200983Z" }, "scrolled": true }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAioAAAGdCAYAAAA8F1jjAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAAC+r0lEQVR4nO29eZwcdZ3//6q+j7knmSSTm3BoCHe4EQUBQVCQVVdFhXUXVxdQwMXfoi54fJeg3xV1BS+8T9SvgLoKGIQAisgVIEAIIQnkzkwyV0/f3VW/P6o+n/pUdVV3VXd1d/XM+/l4zGNmerqrP9NdXfWu1/v1fr8lRVEUEARBEARB+JBAuxdAEARBEARhBwUqBEEQBEH4FgpUCIIgCILwLRSoEARBEAThWyhQIQiCIAjCt1CgQhAEQRCEb6FAhSAIgiAI30KBCkEQBEEQviXU7gU0gizL2L17N7q7uyFJUruXQxAEQRCEAxRFQSqVwvDwMAKB6ppJRwcqu3fvxuLFi9u9DIIgCIIg6mDHjh1YtGhR1ft0dKDS3d0NQP1He3p62ryazidTzOCv2/+KSCCCaCja7uW0hHwpj4JcwKlLTkUinGj3cgi/UMoAo38FAhEgODs+C55QzgNyAZh7KhCizxNhz9TUFBYvXszP49Xo6ECFpXt6enooUPGAUDGEZHcS3ZFuxEKxdi+nJeRKOaQKKfT09FCgQuiUQkAuCYS7geDs+Cx4QjkHFFNATw8FKoQjnNg2yExLEARBEIRvoUCFIAiCIAjfQoEKQRAEQRC+paM9KgRBEATRTBRFQalUQrlcbvdSOopgMIhQKORJ6xAKVAiCIAjCgkKhgD179iCTybR7KR1JIpHAggULEIlEGtoOBSoEQRAEYUKWZWzbtg3BYBDDw8OIRCLUWNQhiqKgUChgdHQU27ZtwyGHHFKzqVs1KFAhCIIgCBOFQgGyLGPx4sVIJKjU2i3xeBzhcBivvfYaCoUCYrH6y/zJTEsQBEEQNjSiBMx2vHrt6B0gCIIgCMK3UOqHIAiCINxQLgBKqXXPJ4WAYGOG1E6GAhWCIAiCcEq5ABx4HChNt+45Q13A4AmzNlih1A9BEARBOEUpqUFKIKLOgmr2VyCiPp8LBeeyyy7DRRddZPm3ZcuW4atf/arhd0mSIEkS4vE4li1bhne/+9144IEHGnyhvIMCFYIgiBmOrCj4/pNpPLen2O6lzByCUXVgZdO/mj+9+/Of/zz27NmDTZs24cc//jH6+vpw1lln4b/+67+a/txOoNQPQRDEDOepXUV8/oEUjhkO4673D7Z7OYTP6O7uxvz58wEAS5Yswemnn44FCxbghhtuwDvf+U4cdthhbV0fKSoEQRAznPGsDAAYy8htXgnRKXz84x+Hoij47W9/2+6lUKBCEAQx08mVFABAuqi0eSVEpzAwMIChoSG8+uqr7V4KBSoEQRAzHRaoZAoUqBDOURTFF2MDKFAhCIKY4eS1gpFsUYGsULBC1ObAgQMYHR3F8uXL270UClQIgiBmOnlNUVEA5Cj9Qzjga1/7GgKBgG2Zcyuhqh+CIIgZDkv9AKpPJTE7+4Z5Sznv6+eZnJzEM888Y7htYGDA8r6pVAp79+5FsVjEtm3b8NOf/hTf/e53sWbNGhx88MF1Pb+XUKBCEAQxw8kLgUqmoADJNi6m05FCaqfY0jQgF1rznKEu9XldsG7dOhxzzDGG2y699FLL+95www244YYbEIlEMH/+fJx00kn485//jDPOOKPuJXsJBSoEQRAzHFFRyVDqpzGCEbWdvY9n/fzwhz/ED3/4Q0f39UNVTy0oUCEIgqiT32/M4gsPpPCNC/uwepF/8yk54ZxKgYoHBCMA/Pt+zzTITEsQBFEn97ycw0haxl9fa1EKoE7E1E+aSpSJDoMCFYIgiDrZPaV2ehVTK36EUj9EJ0OBCkEQRJ3sSZUBUKBCEM2EAhWCIIg6KJYVjEx3nqJCqR+i06BAhSAIog72Tctgp/y8zwOVivJkguggKFAhCIKoA5b2ATpBUdF/ThdpgjLRWVCgQhAEUQe7p4RApdjGhThAVFSy5FEhOgzqo0IQBFEHnaSoUHmytxTKBZTk1jV8CwVCiLho+DbToECFIAiiDvak9BSK3wMVqvrxjkK5gMd3Po7p4nTLnrMr3IUTFp3gOFi57LLL8KMf/Qhr1qzBf/zHf/Db7777brzjHe+AoihYt24dzjjjDPT19WHPnj2IxWL8fo8//jhOPPFEAICiTdtm92fMmTMHq1evxs0334yjjjrKi3/TFkr9EARB1IEh9dNBgQopKo1RkkuYLk4jEoigO9Ld9K9IIILp4rRrBScWi+GLX/wixsfHq96vu7sbd911l+G273//+1iyZInl/Tdt2oQ9e/bgD3/4A8bHx3HuuedicnLS1drcQoEKQRBEHeztoNSPaKYlj4o3RENRxEKxpn9FQ9G61nfWWWdh/vz5WLNmTdX7XXrppfj+97/Pf89ms7jjjjtsBxgODQ1h/vz5OOGEE/DlL38Ze/fuxWOPPVbXGp1CgQpBEEQdiKkfP5cnK4pi9KhQoDIrCAaDuOmmm/D1r38dO3futL3fBz7wATzyyCPYvn07AOA3v/kNli1bhmOPPbbmc8TjcQBAsdhcNzkFKgRBEC7JlxTsz4gelTYupgZFGZCF2IT6qMwe3vGOd+Doo4/GjTfeaHufoaEhnHfeeXza8ve//3186EMfqrntAwcO4HOf+xy6u7txwgkneLVkSyhQIQiCcImY9gH8nU4xp6VIUZldfPGLX8SPfvQjvPjii7b3+dCHPoQf/vCH2Lp1K/72t7/hkksusb3vokWL0NXVhTlz5mDjxo349a9/jaGhoWYsnUOBCkEQhEt2a4FKMiwBUIMBVh3hN8xpKT8HVYT3nH766XjLW96CT33qU7b3eetb34pcLod//ud/xtve9jYMDg7a3veRRx7Bs88+i8nJSbz88st4y1ve0oxlG2hroLJs2TJIklTxdcUVV7RzWQRBEAZ2TpaxcUTPwzN/yvKBIAA1teLXhq8VikrBpwslmsbNN9+M3//+93j00Uct/x4MBvGBD3wA69atq5n2Wb58OVasWIGenp5mLNWStgYqTzzxBPbs2cO/1q5dCwB417ve1c5lEQRBGHjPHWO48CcHMJlTT/J7tNLk5f16Kyq/Vv7kTf6ZXAkoy/5cK9EcjjjiCFxyySX4+te/bnufL3zhCxgdHW2JQuKWtjZ8mzt3ruH3m2++GStWrMAb3/jGNq2IIAjCSKGsYOekGpiMpmX0xgI89bOsPwgJgAIgV1TQU18laVPJaamevpiEiZz6c6aooDsqtXNZHU++lO+o5/nCF76AX/3qV7Z/j0QimDNnjifP5TW+6UxbKBTw05/+FNdeey0kyfoDlM/nkc/rb9rU1FSrlkcQxCyFqSiA3ixtr5b6WdAdRCwsIVtUfFuinC+r6+qNBTCVL0NWVJ9Ktw+Dqk4gFAihK9yF6eI0CoVCS56zK9yFUMD56ZpV8IgsXboUuVyO//6mN72pqq/qoosuMvy91v2biW8ClbvvvhsTExO47LLLbO+zZs0afO5zn2vdogiCmPVM5sSurrJ2m/q9Px5ALARki/4tUWYpqVhIQjIsIVVQPOlOW5IVfPBX4xjqCuCrF/Q1vL1OIRKM4IRFJ9Csnxbim0Dle9/7Hs477zwMDw/b3uf666/Htddey3+fmprC4sWLW7E8giBmKaKiwubkTGsn+mREQiykJn/86lFhqZ9oCEhEtEDFg8qfV8fLeHS7qijcfK6ivQ6zg0gwMqsDh1bji0Dltddew/33348777yz6v2i0SiiUdIrCYJoHRPZytQPC1j0QMXHZlqt5UssJCGhlVN70fSN+XYA4EBGxsKeYMPbJAgrfNFH5Qc/+AGGhoZw/vnnt3spBEEQBibzlQP90oKiEvV5oMLWFRUDFQ8UlV3CUMb9aSp5JppH2wMVWZbxgx/8AJdeeilCIV8IPARBEBxRUdFTP+ptBkXFp43U2LpiIQnJiHeBym5DoFKuck+CaIy2Byr3338/tm/f7mi2AEEQRKsxV/2UZIX3JumKBHyf+jEoKlqg4kXTN4Oikpm5iopfOw53Al69dm2XMM455xzaEQiC8C1i1U+maKyYSYQlxMLqz1mfBiqsPDkWkiArTUr9zMBAJRxW39hMJsOnBBPuyGQyAPTXsl7aHqgQBEH4GaOiIvNAJRxQVQqmqJg7wPoFPfUDMBHd+9TPzAtUgsEg+vr6MDIyAgBIJBK2Pb4II4qiIJPJYGRkBH19fQgGGzNaU6BCEARRhQkbRYX5PXzvUdECqGhIQlBL9jda9VOSFd70DpiZgQoAzJ8/HwB4sEK4o6+vj7+GjUCBCkEQRBUmTB4V0UgLwPceFWPqR72t0YZv+6ZllIVN2KV+fvvMLvzXHzbiWx84Dscu6W/oOduBJElYsGABhoaGUCwWaz+A4ITD4YaVFAYFKgRBEFUwN3zTe6io8kQs7PNARaj6YWQbVH92TRqrfOwUlfs3jmAklcfd63d1ZKDCCAaDnp10Cfe0veqHIAjCz4hm2umCYuhKC/hfURFb6POqn0YDFc2fMphQTyH7M9blybmievvT28cbej5idkOBCkEQhA2KohgVFWFOTleHBSrREJD0qDMtM9IeNV+t5hjPKiiWK7eZL6mv3cY9KaT96jYmfA8FKgRBdCT5UvObjKWLQEkWf7cy06p/82ugwuIDLxWVnVqgcvi8EIJaRmnMwqfCFJWyrODZnRMNPScxe6FAhSCIjmPdphGsuvE+/OLx7U19HrHiBwAyBZk3S2Pt6PWqn6YupW6a0UKfpX4W9wYxoKV/Ri0ClXxRDybXb59o6DmJ2QsFKgRBdBxPb59Asazg8W1jTX0eMe0DGBWVLmam5X1U/KmoWAYqHqV+FvYGdZ+KhaE2L8hRT71GPhWiPihQIQii42AphVSTZQxmpJ2bVA+V+RIwle8sM21eMNOySqVMsf6+J4qiYNeU+vjh7iB/baxKlHOCovL09nHqQk7UBQUqBEF0HOwEOJVrrkGTBSrDPXpp6si0+tzcTKupFH5toW+o+gmzWT/1r3U8q/Dy5gU9QczhikqlZygnBEQTmSK27k/X/bzE7IUCFYIgOg4WqEw3OVBhHpW5iQDC2tFyVEtxMGNq1PedabVAJaxW/gBAsQEfMkv7zEmqAxnnJKulftT7DiQjACj9Q9QHBSoEQXQc7Eo9lW9u6ocFKr1xvWJmJG3uTKutyaeKCqv6iQYlhLUSnaJc/1rHsur/z5SUOQlVbbJO/ai3Hb9Mbfb20p5U3c9LzF4oUCEIouPQPSqtSf30RgO8B8k+nvoxmmn9G6joqR+mChXKqNsvwkYK9MW0QMVGUVEUBTlNUZnfEwMAZArUS4VwDwUqBEF0HDmtmiSVKzXVoMkClb54gCsqTKEwlyf7tZ+ZnvrRFRXA2B/GDROaotIfV7c1x8ZMWywrYG9Nb0JN/WQbyTkRsxYKVAiC6DhyBb2RWDNPfpNahU9fTEIybDxcMjNt3MezfsqyAuZnjQYlRIRAxaqTrBNYOowrKjblyTmhIV9fXO1gmy1QoEK4hwIVgiA6DvEk2ExDLfeoxALck8LohPJksbdLLASEhbl69VYoM0WlN66ePlh58lhWRlnwvrD0nCQBvSxQIUWFqAMKVAiC6DjE/hzNLFGeyOqBSsIUqJhn/ZRkoNSASbUa63cXsCfl/iQvvjSxsISgBLD/olC3osI8KuqWBhIBSABkBRjL6tvMa5FQNBRAIqJGSDkKVIg6oECFIIiOQ+zP0cymb5N57aQcl7iZlpEwmWnVdXkfqOycLOMdPx3D5Xe6L+1lKk8kCAQkCZIkcVWl3piB+Xb6NUUlFJDQHZUMfwP00uRoKIiYFqiQokLUAwUqBEF0HOKVeTMrf8SqH7OikuR9VIR1NZD+2TVVxpceTmGvSTkZ0RqpvTRacq3Y8Pb5gjeF+VTq9aiMs9RPTD99sJ+ZpwfQg8lYOIC4Fh2RR4WoBwpUCILoOLItCFSKZQXTBfXnvnjAoKiEAkBUUyYkSeLBSiNL+fHTGXzjsTTueC5ruF3WxKOSDNfpH+ZRiQpr54pKvR4VU3kyAPRqaaApC0UlFg7yQCXXQOt+YvZCgQpBEB1HXjjhTTep6VuqoP/cHZUMikoyoqZRGF4Yaqe0ACBrSh+JwseOCXeBitg+nxEONDZEUS/Z1rfJFBVx2nRO8KjEwuRRIeqHAhWCIDqKsqygUBY9Ks1RVLKso2tI9WGI5clmv0rMgzb6Ge2xZVN6R/x9e92Bin6brqi4X6uiKLzqx0pRsfKoiIoKeVSIeqBAhSCIjiJfMp7smlX1k9WEmrgWhCQNiorx0OmFosKGGpo3IcYT2yfdpn7U79YeFfdrTBUUrvAYApVoFY9KKIiY9npli2WaoEy4hgIVgiA6CrPPoVlVP0xRYQ3dkqbUj4gngYqmqMgmG4fYQbZuRcXgUanfTDuZZcGHcZs9zExrSP1oVT+CmVZRgHy9LXGJWQsFKgRBdBRmn0OzGr6ZFRXRo9LVxEClbFIcZEVM/bj7X/OWHhX1ez0FOOautAye+snrQQgLSKKhIPeoAORTIdxDgQpBEB2F2efQbI8KV1QEBcFcqhz1YIKynaJSbkBRyVuUJzeiqPCKn7g5ULFXVGLhAMLBAH9e8qkQbqFAhSCIjsJ8RZ5qUtWPOVCppqjo834aeT6mqBhvF3+fyCmYzDlPnbAsWUhonR9poDzZykgLWJcn631Ugobv1EuFcAsFKgRBdBSVHpXmKCqsgsdKUbE10zZQ9WOX+imbAoodLgy1zIgrCCoNKiqVpcmAdcM3vTOt+jfqpULUCwUqBEF0FPkWpX4y2mZjFh6VVpppzYGLm14qrLQ5IImpH/V7PbN+bBWVqFXqx0ZRodQP4RIKVAiC6CjY5ORIUD18Nc2jwsy0XFGp3Uel3iZqgH3qxxy4uPGpsG2FhCM9a/hWT7zAPCq9dmZaG48KICoqFKgQ7qBAhSCIjiJbUE+Wc7ujAJpXnsziH+s+KqZAJdyYolKWFd7zxJzqMQcu2yedB2ZsW0HhSM/7qNTR8I0FKv02qZ9sSVdqxKofAPpgQvKoEC6hQIUgiI6CXZHP0QKVfElGoQm9Ocxm2khQVya8Lk8WH1fpUTH+7k5RqUz9RBqYnjyRtS5PZtOTAWAyr37PVygqetM3gnADBSoEQXQULPUzJxnhtzVDVTGnfiRJQsKiAgjQy5PNc3qcP5cQqNgoKiwYeH5fEdfdM4nvPpGuuV2WNrIy09blUbEpTw4GJL6+yZx6W05ooQ+A2ugTdUOBCkEQHQUzaSajISQ1eaAZPhWzogLoKR9z1U+0gbb0gDFQMWdkWOBy8KAaDY1nFfx6Qxb/58EURqarP2HJIvUTbkhRYR4VqeJveuWP+nteGEoIAPEIeVSI+qBAhSCIjkI0aXbHwgCA6bz3gYrZowIA87vVk+1wt/HQ2YjvA9CNtIB9Z9rh7iC+cHYP3ndUnHeXraXgsG0FAxYN3+pY66RNZ1oA6NUUlSktUDErKtRHhaiXtgcqu3btwvvf/34MDg4ikUjg6KOPxlNPPdXuZREE4VN070MQXdpY4KkmpH7YNGNxps3XLujFD/6hH4fNDRvuy1SKeuOlrLB8u9RPQAI+cEwCN72ll6s8tVqSWPVRYWKQW2FDURTBTGsRqJi60+aKRjMt9VEh6iVU+y7NY3x8HKeeeirOOOMM3HPPPRgaGsKWLVvQ19fXzmURBOFjskKg0q0FKq1K/SzpC2FJX+VhsxGVAgAywsnbLvVjTN9IABSUajyf/WPde1SmCwpPJVkqKnzej/o7b/imyT/UR4Wol7YGKl/84hexePFi/OAHP+C3LVu2rH0LIgjC9/BGYiE99dOUQMU0lLAakQaaqAFArpqiYhFssOqjUo1zPqsY8qIzLVNToqbJyQyzR0V/n8yKCgUqhDvamvr53e9+h9WrV+Nd73oXhoaGcMwxx+D222+3vX8+n8fU1JThiyCI2QX3qERERcX71A/3qISr3w8QPCr1mmlLopnWVJ7MfCZCiTELVGopODxtZPCooK61TtqUJjP01I/6e0XDN+qjQtRJWwOVrVu34pvf/CYOOeQQ3HffffjIRz6Cj33sY/jxj39sef81a9agt7eXfy1evLjFKyYIot3kSvqVeo8WqEw3MfUTc6CoNDI/B9D9MIBeqcOw9Jk4DIzKlh6V+tJU1fwpgFXqx9TwjVI/RJ20NVCRZRnHHnssbrrpJhxzzDH413/9V1x++eX45je/aXn/66+/HpOTk/xrx44dLV4xQRDtJieaabUGJqkmVP2wQMXcM8UKlvrJe1CebI51WOASMKR+JO1v9XhU1O9uhQ1mku2JWr8e5sGEdi30KVAh3NLWQGXBggVYuXKl4bbXv/712L59u+X9o9Eoenp6DF8EQcwurMqTm9rwrQWKijh1WTYFH+z3kJC+CbH0TY0CGt2jYlGe7HKtLCgKB+0CFWPDt7xpKGFcKzcijwrhlrYGKqeeeio2bdpkuO3ll1/G0qVL27QigiD8DjvRxYWqn6kWVf3Y0UhbesCY+jHHD2J5MiPMzbTVgw2eNvKg4ZtV8zgRpqhM5dXgqlA2NXyjPipEnbQ1ULnmmmvw2GOP4aabbsIrr7yCn//85/jOd76DK664op3LIgjCx+SEK/Vezek6ni54+hyKonBFxYlHJdJAW3rAbKY1/s3KZxLiPpPq2y1bBCrMo+J2rUydCUk2ikpUr/rJC0Ybc8O3XK1SJYIw0dZA5fjjj8ddd92FX/ziF1i1ahW+8IUv4Ktf/SouueSSdi6LIAgfwxSVaDiAg+Z2AQBe3peCotQXJFiRL8lgW3OiqIQD3qV+KsuTK7vL8vLkmh4V+9SP20CltqKim2lzQqBCigrRKG3towIAF1xwAS644IJ2L4MgiA5BbM1+2LxuBCRg/3QBo6k8hnpi3jyHIFU4UlS0I2m952Bj6sc8PVn9bkj9aMFG7T4qVo9Vv7ttEMvimlCN1E+mqI80CAUkhILGhm/UmZZwS9tb6BMEQbhBbCQWjwSxQlNVXtjjXV8lVpkSDtibR0XCdaZT9OcTzbTGv8kWAQKv3HHYR8WqM229ZlpR2RHpFqqBRlJqKo4FJ4DQR4XMtIRLKFAhCKKjyGmyBTvxrRxWq/9e3O1loMJ8MM7uH6lTpeDPV3JnpuXlyQ4704oVQ/Uaf5k6Y6eoBAMSuqPqzyMptZkKK00GKPVD1A8FKgRBdBR66kc9fK1c4H2gktFOpgkH/hRAMNOWGldUyrI59WMxAdmpR8WyYqhRRcX+Pv1a+mfXhBqosGZvgLHhm5d+ImLmQ4EKQRAdQ1lW+AmWzZDhioqHqR9eAu3AnwLo6ZSyUhloOKFawzfetK2Oqh/Z4rF62sjdGrmiYlP1AwALe9RTypbRNAB9ICGgK2CAsSqIIGpBgQpBEB2D2CyMXaEzReXVA2lu4mwUt6kfwYpRV/rHSerHqhdKzdSPpx6Vym2ZWdyr/vHlETVQiYmKipAzovQP4QYKVHzE1tFpfPeRrdS50QWjKZKRZxPiZ4OVvQ52RTG/JwZFATbt9UZVybpUVCKCZFEoK1AUBb/ekMGmUWcdc7OFKp1pLYcSOmuhb2WA1T0qdfZRsTHTAsASLVDZPJIBYFRUQsEAIlqU06peKvlS2fMeO0TrcRWobNq0CZ/97Gfx5je/GStWrMCCBQtw5JFH4tJLL8XPf/5z5PP5Zq1zVvClezfh//xhI9a+uK/dS+kIHtmUxUd+MIrfr8+0eylEi2ABRDQUMEwEZumfFzzyqWS1vIiTHiqASVEpA8/uLeK6e6Zw/X3O1lNVUbGc9aM9V63Uj0WzON2j4mhpnJKFOmOGKSqsU7CoqAC6r6hVisoHv/c4Tr75zxhNdd65aevoNDIF7zsudyKOApX169fj7LPPxlFHHYWHH34Yxx9/PK6++mp84QtfwPvf/34oioJPf/rTGB4exhe/+EUKWOpk75Q6JGOMrgAcsWOsZPhOzHxypvkxDK8NtbpHxdn9A5LEDa6FsoKxjLrO8ayzPFC18uSShc8k4jB9U60HS63SZrtt2VX9AHqgwhCrftTfW1eivGMsg79vG0OuKGPneGddzDy/axJnfvkhXPPLZ9q6jkJJRsEHfiJHH8OLLroI1113HX75y19iYGDA9n5/+9vf8JWvfAVf/vKX8alPfcqzRc4WJrWe3WmKoh3BDpxiV09iZmOeyMtgispGjwy1Gd6rxZmiAqgBQFEz+zKrjNO+KtXMtDz1Y9mZtvp2WfM4L8qTa/VRAYAlfcb3JWpSVJihthXp7T9v1JXpejsGt4vtY2pg9cjm/SjLStXXvJn88ont+M/fvoDzj1iA2y45ti1rABwGKps3b0YkEql5v5NPPhknn3wyCgVSBOphIqO+bpl6Z8XPMthAtjwFKrOGvNCVVmRIa+DBgv1GYakJp6kfQEv/FNXutDktT+LE21ssK4YUTkVnWqtZP24VFXHWj/ZYWatQcnoSZNsKV1FUBuMS4iF9oKM5oNR7qTT/Kv3+jSP856J5LoHPYevNFMp4eV8Kr9cUw1aTLujjKtqJo2d3EqQ0cn9CNdCRouIOUlRmH2JXWpGIJjF4deXMrvgTDqt+AONgwjwLVBysJ2fqvWI+p8oW1TZhx4qK9liL8mTAnariRFGRJAlLevXfzQFlq1I/U7kiHtt6gP/eaYGKWL799Pbxtq2D9RNKRto7baeuZ3/88cexbt06jIyMQDYlVG+55RZPFjbbmC6UuPGNFBVnsANnvs4mW0TnwVM/EeMJMKydxb3qz8HKk51W/QDGlAoLUJykfrLF6oFKtc60xRo+E9mqWZxYoSQriMGdolLNowIAi3uATVqMEA3ZKCpNDlQe2jRqqIjqtNSP6AtZv30Cl5y4tC3ryGiSYML0eWs1rgOVm266CZ/5zGdw2GGHYd68eZCEkjmpSiMgojqTGV2yJkXFGaSozD7YCS5mOgEyRaXgUdkrfx4XiorY8ZV7VEqAoihVj43mQEWB8TFWSgYTlGrFZSUrRUV46VwpKg6qfgBgca8EaLOnzYoK96g0uern/o3GyslOU1TE9a5vp6LClMVOU1S+9rWv4fvf/z4uu+yyJixn9jIhBCoZaobkCN2j0lkHIaJ+7Kp+WH+OgkcnpLoUlZCe+mHpHAVqCXG1C1JWmhyUdPWkrADsqauVGJdqKAVWaSNJq1Aqyu6UBt5HpcYF6RLBThE1ByotUFQURcG6TaMAgO5oCKl8qeMCFVFR2TKaxmSmiF43eUiP8Iui4tohEwgEcOqppzZjLbOaiaxuQE571F1zpsOOPWSmnT3YVf1EPfColGUFv3piB7aMTtdnpuXlycZ0ZK35Pxlt/+0Spg+L59WylUeFtcGv2ZlW3bbZVlJPd1onnWkBYHGP/mTm1A8LMJtZ9ZMvydzvd+j8bgCdl/oxB1brd7RHVeEzr6IdFqhcc801uO2225qxllkNKSruYZK42YxIzFz0QMXao1KWlZqzdqZyRTy7Y6Kio/FjWw/gk795DtffuUHoTOt8baKiYghUTB/nv7yax0+EJoUsddkV0U/w4r/QSGdamftKjJGK00BHxKlHpbqZVmv41sRARWwm1xNT38BOVlQA1afSDrI89dNhHpV///d/x/nnn48VK1Zg5cqVCIeNctSdd97p2eJmE2JZJXlUnMHsCIWSejAPkEdqxsPMsnFz6kc4exZKsmEAHmPvZA63PfgKfvP0TmQKZXzjkmPx1iMW8L+Pa+0BXtg1iVXD6pV4rA5FpVg2Bs/myp//uG8KOyfLOHlJBAcPhrhHpSsSAKD+f2qwxYIR9XGBOqp+2DKsFRWlphnXsC0HVT8AsEhI/diWJzcxUGG+ikgwwPeDUocFKnltvYlIEJlCGet3TLRlHUzdj4fb61FxrahcddVVePDBB3HooYdicHAQvb29hi+iPsRAhQZ2OUO8cqb0z+yAfTYqPCpioGJzUrru/z2Lnzz2Glcsd09kDX9nV93pQhlb9quKR8JFoCJ2ixUL98ypn+m8+jx7UuqdLFM/wkOq9lGpEWiUeXBhWqsQVDnFqaKSCEuY06W2qLBt+NbEYxzzVcQjQa60OW285xeK2j5z1KI+AN7NsHILL09uc+rHdZj04x//GL/5zW9w/vnnN2M9sxbW7A0gj4pTxKvJfElBnNr3zHhY6sfcgEpMbdi1/N68bxoAsGQgge1jmQqfRFEIKPZPq59HV6kfiz4q6u/G+7HlHUirPzAzbdIu9WNVYswUlRrne7YdswE2LKzVKSWLDrl2HDQYx/7pAk+9MFrRR4X7KiJBhALMu9RZikqhrP4PQz1qI8N0m1pWiK9lO3GtqAwMDGDFihXNWMusxuxRoYnAtREVFSpRnh2wqbvmhm+SJOklyhYnpUJJxr6UOkvriEWq8mvuuZK3eJyr1I/g+8hXSf2wFMp+bR4Q23eTwnOJ3hNLRcXx9GT1e8B0pA8LPV+cwhUVBy/J9eeuwJVnHIzTD51ruF1P/TQvcBBPrsw3VPTBvBo3sGC7L65aKzKFUsU54f89tRMf+8V63q25GeivZYelfj772c/ixhtvRCbTWUOe/M6EkPopyYpnZZYzGfHzSamf2QFrhmglRUeZzG9xUtozmYWiqFUoC/viACorT6xOZm7Kk8NCOsbgUbHpPHtAC1RY6icelriXROyjWa3qp57pyeJaXSkqDqt+AOCYxb3497ccZt+Ztompn2yRldSGeOrHjRfHD7Aqpd6EKhPLSmVg/eU/bcLvnt2Np1+baNo62PTmdisqrsOk//mf/8GWLVswb948LFu2rMJM+/TTT3u2uNmE2PANUA/I5vwuYYQ8KrOPaS0tmoxWHroioQCQt5b5d46rfpSF/XHeLM584Ld6XNxVC331e6FknPFjDgZKpkAlKwQqQUmbwWPwqFSWGIvN5apRtjHARhw+3mpb5goiN0SrqF5ewVSAeCenfrSdpFfYATOFMg/0MoUS9kyqCiELzLxGlhWh6qfDGr5ddNFFTVjG7OD5XZMY7IpgQW+84m/mYWrpQgn9STJdVEM8z1CJ8uyAVcR1WQQq4SqKyi4tUFnUn+BNyCoUFatAxVULfV1REdM9onhQlhWwv7BAJaedZxJhSU3RyMbBhFYlxk4701qljQDnioyIG0XFDq87CFshpn7CnZr6KevVbZFQAIWSjEyhhAHtnLB1NM3v26x2FrlSGWw37DhF5cYbb2zGOmY8+6ZyuPC2v+KQoS7ce/XpFX8XG74B1EvFCSVSVGYd0yz1Y3GFF7FRSgBgp1bhs7Avzq/qzfezCnDcKCo8nVKyT/2IT7E/bUr9hCStV4piTP1YzvpRvzudnmwOLupJ/Tit+qmG3WvvJVnRo6L947W8PH6D7YvhoIREJIhCSTaky7bu1wOVZqXRxHOQuR1Aq6lbz3nqqaewceNGSJKElStX4phjjvFyXTOOneNZlGUFW/en+RyP0VQeqVwRB83t4mbaYEBCWVao8scBZeHzSWba2UG6VuoH1gHHznHVU7eoP27bHdXqpB2rYyhhQTZ1phW2K6Yr92fU52epn1hY4gGFZXmyoY+KsceKHXrayLrhm9fTk2tR7T3yCqa6xcMhnvrpNM8fe30ioQCSkRAmMkVD4LB1dJr/3Kwuv8wPFg8HEWjgPfcC14HKyMgI3vOe92DdunXo6+uDoiiYnJzEGWecgTvuuANz586tvZFZCMutF0oyckW1IdU/fudv2DWexf3XvpFfYczviWHXRJYUFQeQR2X2wQIVq9QPu3q2SuHoqZ84PwnU8qgEpeozesyIbentypPFVMuBjAxFUSo8KoCxPJn7TCSr1I/TzrTG2+vyqLBSZw9SP61QVJLRzk39sH0xGtKb1omNQMXUT7NKvTNFdlHQfq9kXQ3fpqam8MILL2BsbAzj4+N4/vnnMTU1hY997GPNWOOMIJXTPSgT2QLKsoJt+9PIl2Sse1kdoBUMSJjbzermSVGpBXlUZh+6mbby4BmuqqjogYqdosJODkntxBAPuZsIHxHLk20avpUN+6ya9tk3rd55TiLA1QoxCGePMaZ+WKBRfU0li8cC9XpUKgMmt7ACgWYqKqKZtmNTP2WW+glwf4gx9aMrKtlCc15L1rvFqstzq3GtqNx77724//778frXv57ftnLlStx2220455xzPF3cTGI6pwceE5kioqEgNyr9bct+AGrNPLtSJEWlNuRRmV0oilJVUYnaTFAulWXsnVIrJBb2JXBAa+Zm51FZOdyDJ14dR8zl0dHQmVbYH0VjrfmEeSAj49Vx9bO+vD/EAwpR6JAtlIwIDzTs93tF0Y275nRNuz0qrQhUEuEQD+g6OfXDAhX2fymKgm0tUFS4MtXmih+gjkBFluWKkmQACIfDkOXO2hlaScoUqISFhPPfthwAoJaiJSxkPsIa8qjMLnJFmZ+03XhU9k7lUJYVhIMShrqjgqJiClS0k9nJK+YgVyjhhKFpuCEc0E/+hqof4aNsPj+/Ol7mjd+W9Ad56scwPdmiIyxv+FblHCXGILZVP03qo2KHbqZt5lBCvfcHU9k6LfVjVFTYxav6f42k8kgLF7JeelTypTJeGZnGygU9utfHB4qK613uzDPPxMc//nHs3r2b37Zr1y5cc801ePOb3+zp4mYSKSGVM5ktGFrmj2tG2t5EmB+AM21qmdxJGDwqlPqZ8bC0jyRZl0vqKoHxpMT8KcN9cQQCku3Jkp20++Jh/P7fVuM/T3d3eGQXnvmS0ZdiZ6YFgKd2qceBwUQAPdGAYKatrBQy9lHR1lxFURFfhspZP85SR8btNd5HxWsz7WSmWNGxVUz9hDs19SMoKnGTorJl1BhAe1n188V7NuH8//kL7t844itFxXWgcuuttyKVSmHZsmVYsWIFDj74YCxfvhypVApf//rXm7HGGYGY+hnPFDGWLlTcp09QVCj1UxvDrB9SVGY8vOInErL0jtidBEV/CgDeRyVvUlSKwsmhHpiikjJ5BoypH+Njnt6tXqQs7VPXxKpzxCCDCdVWQwmrne/FoMiuM23rFRXttfcgUHl+1ySO/T9r8bnfv2i4XW9SFhT+z85SVNh6I8EAEmHjOUE00gLepn62j6nbfnlfyhDwtRvXodLixYvx9NNPY+3atXjppZegKApWrlyJs846qxnr6zie3j6OG3/7Aj59/utx0kGD/HaDmTZTtDSk9SUiQqBCqZ9qyIpiqIyg1M/Mp5qRFgAiNkbNXUIPFQCIha0VlYJwcqgHdjxP5Y37ouiLN8cF67VAZVm/+mDLqh/L1I/6vSSDtzswIz6XubxULKV2ipdVPyVZQVlWGip13rBrEmVZwQu7Jw23s4A2EQmhpEV5zfTENANDebKmsmdNgUpXNITpfMnTQIUFkGPpAlcek50WqJRKJcRiMTzzzDM4++yzcfbZZzdrXR3Lfc/vxYZdk/j1kzsNgcq0cLSayBZgZZxXPSrqW0IeleqUTZ9NUlRmPtV6qAB6gGFO/bAeKgv7EgD0q3qzR4VdxYZd9E4xPr/6uGlToCKmfswKxnRB/X1pn/o/8dSPqKhU6aMCqMGKVT8uMdgxBxf1KSpeVP3oCymU5Iau1se19Ln5RC12pmUNvzst9cPel0iwsjyZVfy8fkE3nnh13FOPSkEIVFj7/ninpX5CoRCWLl2KsvksQXDs8ohioDKZKfIPmXhB0RsP86tF8qhUx3zgofLkmQ87UNvlzO0m5TJFhaV+bBUV3g20ztQPC1RMMoXRo2L9WK6oaAcE2cKjYpyejIq/mxFvtzfTWj/WCi+qfiKmQKUR2Hw0s0djJqR++L4Y0lM/ZkXl8OFew+1eICoq+uet/YqK613uM5/5DK6//nqMjY01Yz0dD/uQbBmdNpi8pkxVPxNp9UN25KI+fntfghQVp5iPcaSozHymq0xOBuwVlV3CQEJAV1SKZcXg4yg2mPoJO0j9lBTr/XRpP/OoqL+z2EYRUpwGM60QediVGLP/TUJlP5hII4pKA+maUEDianKjlT+sm7dZGbMy03ZS6kdRFEMaMiG0rFAUBXsm1f350HndALz1qIiKijiKoN24/kT+z//8Dx555BEMDw/jsMMOw7HHHmv4csNnP/tZSJJk+Jo/f77bJfkKttOkciWMTuf57dOCR2U8U8CYpqicdvAcfntfQlBUyExbFXP1BHlUZj7VeqgA9mZadpHABroxRQUwnix5SWidkgFP/RTsUz925+ZlLPVjKk8Wd3MxQBBTPXbnYKuUEX88K6V2kRLxQlGRJMmzeT9sPpo59aGfYEMdWfUjBo/GPiolZApl/vfhvhgAbxUV9hkYSxd4w7eEzeetlbhewYUXXuiqW2MtDj/8cNx///3892Cw/dFbI+SEnWbLSBpD3erOZEj9ZItc2l053IP+RBjjmSJ642G+E1Jn2uqQR2X2UdOjYnMCzGsnMnaCZIqK+jcZCW1IebGk+wLqwe7Cs1Cl6gcAemMS+uLqc+oeFfUxhl4owrICkoSApAYjJTtFpVqgUtesH/vtuSESDCBXlBsOVFhbB1FRUBSFFyJ0aupHVAQjQWPDt0nNdBMOSjzw9tZMq25rLF1Atqi/ju3GdaDy2c9+1tsFhEIdr6KIiDvNltFpnLxCNdSaG76xD05/IoIPnrwM9z6/F8cu6cfzu6YAkKJSC3P/COqjMvOZrhGohG1m/eRKbG6K7gMJByUUywpygsRRFJps1UPEbATRKBiqftT9dDARwAGt0duyfv3/4eXJ2u4sV/GZhAJqvxa7NvhWM4IYbs20stDltpE+KoBWHp4reedRKZZ55VO+pDcFTAipn1Y1fNsyOo3F/Ym6S9wB41pVRUVP/bB0l9gctBmpn2yxjP1aB+dEp5lpAeCggw7CgQMHKm6fmJjAQQcd5HoBmzdvxvDwMJYvX473vOc92Lp1q+198/k8pqamDF9+wxyoAOoBQww8JrIFvsP1J8O45uxDcd81p6vlyVHqTOsE8wUSpX5mPvWkfkplmZ+wxZQP7+chnOW5L6DePio2gUreUPWjfp+bDIDde1mffsXKYiR2si0ZeqGY2uDzCcp2HhX2uMq/RVwqKgZjrgeKCtB4W3uW+lEUXUUTj7Ni6qdaYzyv+PvWA3jzlx/CDb99vqHtsNclGJAQDEiG1A9TVHrjYd5h2dPUj/BGM2+XHxQV17vcq6++aln1k8/nsXPnTlfbOvHEE/HjH/8Y9913H26//Xbs3bsXp5xyimUgBABr1qxBb28v/1q8eLHb5TcdcafZormzp01pnFxR5h6VAaY7a7CKBqr6qY7Zo5IvKhUdKomZBTfT2lX9WJwAxfSCmPJhQYuoqOhVP/UpBlYlwup6KjvTRkMS+uPq8zAjLVDpUTH2QjFuV5+gbP281VM/7jwqYkxRZ/U2J8qqrhpUAtjFHqD7VFjaJxIKcOUMaE3qZ+Me9cL5lRF3oxfMmPfDuCH1o543euNhxFnjwpIM2aNATPy8sPlYfghUHGs6v/vd7/jP9913H3p7e/nv5XIZf/7zn7F8+XJXT37eeefxn4844gicfPLJWLFiBX70ox/h2muvrbj/9ddfb7h9amrKd8GKaOzaou2wrNlbJBhAWVErDdg5tc8UqNCsH2ew80s8LCFbVGXpQhnwge+LaBLpmg3fKhUVY6BSXVFptOrHSeqHZShDAWAwEcRYtmRI/fDpydoBomovlIAEQLFN38gWjeL0x6rfnSsqgrLTYOrHC0UlWygb3ttssYw+oKJSpZWpH9Zt3Hxh6hZz40FxejJTVPoSEUMPmlyp7EmKRvzssKDaD6kfxyu46KKLAKiu7UsvvdTwt3A4jGXLluHLX/5yQ4tJJpM44ogjsHnzZsu/R6NRRKPRhp6j2Yipn10TWWQLZb7j9sRDkBV9h+6KhipkZpZ/zxXlhjs3zmTYgTMZVQMVQFVVxtNl9CeDiDZ62Uf4jpqpHwuPCrtwiAQDhu6s7KpevLBgJ/z6PSrG31kQnbeo+glKwEUrY/h/z2dx2lL9YoV3pmWKihAgBMypnxrpG6sZQfpj3XlUDIpKg6kfuxEGbmBpHwYrUdYnJ5sClRakfvZrx3XRj1gPYldaQFDZhUClNx5GTFAIs4XGA5WyrFimEf2gqDje5WRZhizLWLJkCUZGRvjvsiwjn89j06ZNuOCCCxpaTD6fx8aNG7FgwYKGttNOzPnCrfun+Y7bFQ2hL65Pnu5LVE6hFncKaqNvDy+VDEr8gL1lXxFX/mg/brlnom3rIppHLTOtVdVPnhtpjYc6q5kzXntUeqJaekUwerMTQTgo4YqTu/Dg5XMx1KV/5ll6h8UP1XwmzNRalBXIioL9aeOxx2pGEMO9R0UMmJw9xo6oB4qKmPYB9OOueT5NSEj9NDs1PDbNApVijXtWhwcq2usUF0yzY2k9UBEHbHphqLUzN3dUoMLYtm0b5syZY7htYmKirif/93//dzz00EPYtm0b/v73v+Od73wnpqamKhSbToJF9myuyJbRNB9I2B0Lo1cITlh5mUhUy60CVPlTDVaSGQpKiIbV1+uZ7XkoAPZMUIA3E2HpUDdmWlZuGQ0bD3Uxk6KiKErDnWnNqZ9uLVARFZVqvhFAN8yy1E9Vn4kw7+eGtVM4/rZRbNirnySdeFTyLhWVUKCyeZxb9IBSfe3/97ndWL993NU2KgIV7X1kJbUsmGUne0Wp9LV5zYG02jdrOl9qKCgqmgJmMVAY0Xwjent7Ng6imYFK+1M/rj+RX/ziF/HLX/6S//6ud70LAwMDWLhwIZ599llX29q5cyfe+9734rDDDsPFF1+MSCSCxx57DEuXLnW7LF9QKsv8KmHVwh4Aqk8lJUjW/YInxexPAdSDgHlaJlFJSThwxrQ0z9YR9eDlYbUe4SPSvDOt8/JkduEgGmnV343qi6gYNNqZltEdZcqBfht7npDNyb4i9aOd8KxUDH2CsoIXR0pQAGw+oAfp1TrJsoohp6mfUo0Ayw1RIaDcOZ7BlT9fjyt/vt7VNiYy5tSPSVExpX6A5jd9O6ClfmSlMYXDHDDHQkHezXe31pWWByq88qdxD45dp+CEjSeslbje7b797W9zA+vatWtx//33495778V5552H6667ztW27rjjDuzevRuFQgG7du3Cb37zG6xcudLtknyDuHOu0uYwbBmd5lJgV8yY+hmwSP0A+o5BTd/s4T0iAuCKytZR9fWya4BFdDa1pydbVP0UrRUVfTCh+ncxuGl0KCGDpX7EHj/co2Jz5A0IKgmgByxWvUtYhqpYBjLMpyU8F+9Ma+lR0R7r8PxWrhFguUFM0bFeHSOpnCsVYiJrk/ph3VRNqR+g8XLoWjDvIQCuoteDOQUZCEg8INkzqSoqzDbAAxUPrs7sGvAl7MrZWohrTWfPnj08UPnf//1fvPvd78Y555yDZcuW4cQTT/R8gZ0E21kkCTh0vjqHYftYBkcsVIOW7liIR8KAtaICMPNUnhSVKpSFA3hMe0lZLxUKVGYmmRpmWu59sPSoGA+2+mBC9e/FkklRqWMXMsVC6OGKSmXqx+7YzxUVU+rH0hAr9FHJWgQq3N9iERS5nfXjVVdawKioMB9esawgX5J5b5Ba2KV+9K606j4SFmq6m1n5UyzLhjVN5UoY6qlvW1YpyEQkiEyhzAMVdh6JeRioWAVykVAAIS/e9AZxvYL+/n7s2LEDAHDvvffirLPOAqDmeGf7VOWcJr/Fw0EsHVRHym8fy/Arwe5oCH1xPTix8qgAoKZvDtA9Kqio8KHUz8xDlhWkC9VTP9XKk+3MtExRyZf1i4x6K+2CAcmgXjCPilXqx6pbLFDZmbZasCH2UdEVFf3vvDzZsjOt+t3pZ8WLOT8MUVERiw+mXJhQK6t+tEClaDTTBgISV6OamfoZN6WiGilRtjJ1s8CL7dtmj4oXTd+sqrD8YKQF6lBULr74Yrzvfe/DIYccggMHDvBeKM888wwOPvhgzxfYSbCoNh4OYnG/GqhMZIq8w19XLGSo9Om3S/1Q07ea6Fd4UsVB3M1EWKIzyAhnVDtFRfeo6O8/O4HFbMy0XFERSpMbMYtGQnq5fE9MfQ5ZUU+SoYDEUz92J3x91o/2vVqwwQYLlgVFxWKukLlRnPmxTvBicjJDrLhKCyfYVK6EoW5n22DT5xnsfbaa+BsKSijJSlMnKB+YNgUqDaR+WBoyGjIqKiLsPKK30W/8otZKUbFrrthqXMfHX/nKV3DllVdi5cqVWLt2Lbq6ugCoKaF/+7d/83yBnUSWHxSDSEZDmNOlKiYvah0Lu2NhY6Bio6gkqelbTUSPSixsPHiWZFCX2hkG82sFhZJMM9XLk81mWtbLQ/OomEpC60WMh5hHRV2HlpZU9Go1K9i/JjtQVMKWiopFhVG18mTHHhXj+hpBVL6ywjHOTf8Rs6KSNZlpxUoVuxlQXiL6UwD3Jco7xjL45P97Fq8dSFumfuKmQKXHnPrxwEzLnrdbuBAwP2+7cB0uhcNh/Pu//3vF7VdffbUX6+lozNH8koEE9k8XsFnrUNsVNXpU+m08KmysdobMtLbwqh+hPNn8dx94wAiP4EbaSNBW8dBPgPpVOi9PDtVSVBrrocLXEFS7xQJ61Q+gpn+SEE74NsJEwFSeXM0Qy1Ia6YLC72cw01ar+nHrUanS5dYtYnlyxqCouEj9aH6QnlgIU7kSP1FnLBQVFnw2M/VzwByouDx233zvS/jDc3vQn4hg0YCqxkeC9opKRdWPJ2ZadRvzemNIaeesZKcGKmIrfRFJkhCLxXDwwQe7bqU/U8iZ8qNLBhJ4evsEv/rvjpnLk61TP7qiQqkfO8rMoxIAYhYnl1JZqXtmC+E/anWlBaxbs/Py5LC1osI9Kg3O+dHXoP+cCKuelbKiN31jgYF9HxX1O6v2YSfXQJWqn6m8oCCJfpiqVT/Mt6F6Wcxdb83UCrDcYDTTGlM/TmEdWhf0xjGVS1mYaY2pH/Z8zeLAdN7wu5vUT7ZQxgMbRwCoAdhQj7YvWnhUADU4YfsvC1S87KOSjATRHQshlSt1rqJy0UUXQZKkCmmd3SZJEk477TTcfffd6O/v92yhnYCY+gGAJYNJw9+7TR4VOzNtLGyUpYlKRI+KlaJSLAPxFq+JaB61utIC+pW6qBKwq0RzMGunqNTb7I0hBjqxkMQ9K8w7wpZml/oxd6Z1UmI8lRP/X1FR0R5bJW0EqJ+VWjOyvKz6sTPTulFUmHl1fm8Mm/alKvuoCCfYZqV+pvMlPPrKfrzxsLkVqR83Ztp1m0b4uWO6UKroTAsYAy/xHOKlmVZs3T+YjCCVK3WuR2Xt2rU4/vjjsXbtWkxOTmJychJr167FCSecgP/93//Fww8/jAMHDlimh2Y6WVOzoSWahMfoioYxkIwgIKmyrV3qx9yMiqiE93Ww8KgAVKI806jV7A3QD+xlWeH7R54rKrX6qCiGbdSLGABEQwBr+cLOI6UayoQ+PVkLbBy00DcoKgaPir0RNypsMFeq/VnRP2/emWkbUVRY6mdBbwxAdTNts1I/31z3Cj78k6fwk7+9xvvBMNwEXf+7YQ//eTpXskxDiv+PaB/wsjxZ9HMx/2THKiof//jH8Z3vfAennHIKv+3Nb34zYrEYPvzhD+OFF17AV7/6VXzoQx/ydKGdgFj1A1QGKt2xELpjYXzpnUchHJRsewZYmQIJI1YelYhwUqDKn5mFo9SPcGAvlGTEI0HnfVQ8UlTEpm/RkMQ9KyyAKNeontGnJ0P7bn9/nvoRFRWxZ0vVPirq40sykC0q0M73tjRHUSkbqlWmHAYquaI+OXm+tnBz6ice1vcTPu/H4+Mpq+Z8fNsY7xw71B3FSCrvWFER0z6AqsToior+noupnx4hUPHSoyIqKrGwGqj4RVFxvYotW7agp6eyk01PTw+2bt0KADjkkEOwf//+xlfXYWRNsiPrpcJgB9l3Hreo6nYiwcqBaYSRkpDrZ31UhvtD2J8qo1BWqJfKDKNWV1rAGGSwQMWuPNk8lNA8sbZexNRPNKimfgC9DJhV2dgZvQNcUTF+r9a0zaio6H+v1ixOHdUhYSqvIF1woqio372o+rH3qDhTIZiaEgpImNMVBVA5lFDcT8IW3iUvYPvkhl2TWNSvJpqXDiYwkso7VofEtA+gBuS1FBWxu3k8os2s8qKPCnveYADdsZC2fX8oKq53u+OOOw7XXXcdRkdH+W2jo6P45Cc/ieOPPx4AsHnzZixaVP1kPBMxe1TmdkUN1Qbsza8Fk6mbaf7qdPQDuISDhkIISMDRS6P6/BNSVGYUaQceFTFIYCcl2/Jk01DCAldUGkttiHFULCRxlY8FENVSOeLtvDy5SrBhqagYOtNWN+4mI9rw02Ltz4q3fVR0NSudd5/6Yf6UvkRY6CNSNnxPWHhUvD4msPXumczxys6lmi/R6f/y55dUNeWE5QP8cXmLoDluk/rxtOpHGDdxmNZZffmcZLWHtAzXisr3vvc9XHjhhVi0aBEWL14MSZKwfft2HHTQQfjtb38LAJiensZ//ud/er5Yv2NO/QQCEpYMJAzlyU6wql4gjPDhbkHg0PkR/OjDQ4hHJPzlZVWOtUv9PLkthye25vHPb+zhV7uE/3GS+pEkCZFQAIWSLAQq1uXJZkWlGWbaaEhX+5iiUqrh9eCpH+1+cpXAJmSlqAj7PQt27ObzsEBl2kEPDm4C9lhRkSSxj4o7RaU3HuYXhZVDCcU+KqwUuzmKirimZZqK7jT1w9JHJy4fwOPbxtTUj8W+KM7bEc20zWihHwkG8E+nLsfJKwbx+vl1zgHwGNeBymGHHYaNGzfivvvuw8svvwxFUfC6170OZ599NgKaZf2iiy7yep0dgZ760XcwFqjEw0HHMxN4DpfyF7aIs34AIKH1rBDLLq341d+nsWWkhGOWRnHSwTUS84RvmHZgpgXUeT+FkszVSGamNfvBYiZFxas+KkYzrcTTM3qgov7NtjzZVPVTzaPCsllTeWtFpVpnWgBIhgMAysg4SP3Uav3vBnF4ZFmoHnWqQkxmmaISqVAUWO8pK0XF6ws/q6GxrNLTaXnyfq2smSkXadGjIqZ+hP3eoKg0oYV+JBRAMCDhcG2wrh+oyykjSRLOPfdcnHvuuV6vp6PJmRQVAFisGWq7HKZ9AOGKgxQVW0o2/ShqNbJKawf13ePUTK+TcKKoAFrvibyeNnWqqFiVhNZDxFSeHK1I/Wit+u0UFZvUj+WsH20bkzlrj0q10mYASGiKStpB6sfTzrTMg1eUIcY9TgMVpl70GRQVGYqi8FELrUj9mFWTgATuVXGqqIxqgcoyLVApyQpXluzKk61SP570USlbp0n9QF2Byp///Gf8+c9/xsjICGTZeDL9/ve/78nCOhGzRwXQDbVO/SmA9XA1wohZUWGwA6nd55ZNWN49QWpVJzFd0DvTViNi6pnBPSpmM22FR0Wf9dMIYgASDUkWqR/1b3ZPo5tpa5cns39JPEwYq36q+0qSYb2zbS08nZ4cFhQVWVRUnKV+WLO33njYoCjkSzKYQGPso9Kc1I85sOpPRHgQ4WTAYr5U5kHXUqFCdFybY2Rbniy0tYhHPEz9eGQobwauA5XPfe5z+PznP4/Vq1djwYIFDQ3wmmmYq34A4NB5qilpXrfzNIP5ao+oRPSoiIRrmGmzBRaokKLSSWQcmGmBytJ+3kfFXJ5s9qiwzrQNt9DXf44GYZH60fv/WFFRnsw601ocZ62axjmd9QPoHhVnVT/e9VFhwWS+WDYMRXSqqLBuw/FI0KAoiBVEzZ71UyzLFcfngWSEX5BO50u8AaodbJBhOKj21OqKhjCdL2FMMwsbFRWb1I/2/2e8SP3YqI9+wHWg8q1vfQs//OEP8YEPfKAZ6+lorBznp6wYxH+/6ygcs6TP8XaipKjURKz6EalW9VOW9X4Weyj101GIJ6dqsECVfXZyrDNtDUWl6FHVT2UfFfXnioZvNVI/bP+Wq6R+rEqcjVU/9o8F9NSPs6qf6ttyg6ioZOto+MZOqJFQwOBRYb6QWDhgOC5YTdVuFNGfctDcJLaOpjHYFUF3VA0iFEUNHqoF1syfMqcrikBA0gOVNAtgnJQne99C34+KiusVFQoFQ7M3QsfKoyJJEt553CKsmNvleDtiQyTCGnuPivrd6nObFw7IUzkF0zkKBDuFHL/aq5H6YV1PyzUUlbBRUWEH6UavJtn+Fw6oQTQLXPSGb+rf7Z5Gb6FvbLlvGahYBDvirB+2Dbs5Pm5SP15W/TClIFeUuacEUN8zJydcseQ8phUuZItl7vcYTEYN9w81IfWTEoKiY5f08+cVg6RagddoSg9UAL33ywRTVBx0puWBmoct9P2oqLhe0b/8y7/g5z//eTPW0vFYeVTqgRSV2thdmYaqmGmzpivHPZT+6Rj0gKP6IYvP+6lpptU/Y7KseF6ezLwp7HueN3yr0ZlWCyrMQwmtqm2sXop8SeFz2OQaQZGe+ql9nPGy6ocNiFTTI8a/OVFVxPeUHWsVBdg9oZb6zuk2Bipm35IXTHNzdxhnvX4IAHDMkj5IksQN39P56j4VFqjM1dbbFVMDEKaihV2kfrLFcsX8PbfkfRyouE795HI5fOc738H999+PI488EuGwcQLwLbfc4tniOg3zrJ96ITNtbfiVaYVHRf1uGaiYrhx3T5SxeA55rDoBp/nzqKkUNWdbnqz/XijLnplpmYISC5nGOjhUVMzlyVX7qFgEO7KiBvHhoFCebOtRUZ/MSerH26qfyo0kI0GkC2WkckV+4rZDnN8kHmt3aj1J5nYZZ6g1I/WjBypBnLtqAZ654Wz0aSbX7lgIk9lizaCLpX7maopKtylNJO7rfYkwoqEAIqGAoYV+TNvBZEXdjxup2PFz6sd1oPLcc8/h6KOPBgA8//zzXq+no7Ey09YDmWlrUy5bX5nqZtrKx+SK5kClBCBceUfCd9gFHGbCIaNHpZaiom677OGsH7Z99jsz06q/12z4Zm6hzzrTOvSoAKqqEg5KVXuwAEDCReqnlhLkhsoBkQH0xsNaoOJEUdFTP+FgAKGAhJKsYMd4BoCeSmE0I/XD/DCs7USfUImjKyoOUz/d2lwd03gIY+onhJ/9y4m8xwlDDNRyhQYDFY96CTUD14HKgw8+2Ix1zAjMnWnrhRSV2thVT1Qz02ZNEvee8TIoUOkM7MqMzZi7Otu10A8H1QO+arCWPbuaZIFJ1JQCMntU7Bu+aakflr6pUrljF+zky0AXqhtxAbdVP+w5a961JmZFJREJojsWBiZzrlM/gHq8TeVL2DFmHag0N/VTeQpllT81PSomRaUrajwWmYPm1csGKrYhBmrZYhm9DRzP7PxcfsCT0EmWZfz+97+ftR1pGV55VGh6cm3E6cki1VI/1ooK0QnkHZtpjUF+tQCHnehERSXSYNUP2/8i5tSPeShhrT4q2q5abcaOXdBQERTVKk9ucdVPZaAS4sqEk14qZi8FS3+w1M+clqZ+KgMDrqjUSv2kVNPsXK11hbnXltOg2at5P+JQQr/R0Io2b96M66+/HosWLcK73/1ur9bUkciy4riEshZiZ9pGDVIzFbsrU13mrXwM86gMJNUH7Zlo3IBGtAanRj9e9VOSURQaisUsAhyx8sdqvko9hM0elZAx9VOrCRvvTMtSP1WCDbugKscDFWaAtV4rH0rY4j4qgYBkOBnGI0FdhXDQ0VX3qKjvHztRs7k5ZjNtU1M/FtO8uzVTbK3/hSsqzExrUmeclsrH6mij/8cNe/DbZ3YZbis4VC3bgevUTzabxa9+9St873vfw2OPPYZyuYyvfOUr+NCHPoSuLucluDMNUf3wKvUDNG6QmqmUanlUZIvUj3bluGxuGBOZPHJFBRMZBSHK/vgaWVb4QbSmR4V7QowNuWorKt6aadk5J2oqTy7VSKHo5cnq92rpm5qKCve3VPeoOBlKWKqhzrglEgrw4DDJUj9wV/XDgh12vGXbM6d+mtHwjQUhVqNRnKpDenky86jYm2mr4VZRKZZlXH3HMygrCs543RB6tNfe/Lr6Cccrevzxx/HhD38Y8+fPx6233op/+Id/wI4dOxAIBHDWWWfN6iAFMO4kDad+hB2F0j/W2B3wQzz1U/kYlvrpjkkY6lHvuIda6fseceZVzaofoTxZHOppdfAVFZWiRx4VdoEdM5UnmzvT2npUJGOgXa0XilVnWkAvha7V+6TLVdVP9Y66bhHfR4Oi4iD1w2fSaMGnuZmfnUfFy1k/6Sqpn24HqZ9soczTR3p5sin1E3R2HnE77ydbLPPxBaw7LuDvqh/HKzrllFOQTCbx+OOP44knnsDHP/5xzJs3r5lr6yhYoGJ2ZdeD+CEmQ601/MBZ4VGpbaaNhwNY0Kd+uEcm/fH6TmaKeOjlUcPsE0KFSf2Ag9SPYKbNCQdeK1VBVFQKHuXn37A8ihMWhfHuIxLa9qCtR/3OYi6nQwmr9UIxi0SsgRvLOLDH2h2OEi7MtLpHxRtJxVzR4tSAClT21DFfGM61qfrxcsgrC0KsZrg5qfphpcmxcIDf31yezCrYauF2grIY0LDmcoDY8M1/Cr7jT+WZZ56J733ve/j85z+Pe++9l3L7JrzqoQKo3Wz5AZcCFUtsPSqBKg3ftANyLCLxqzAPGjp6wpp7NuLS7z+OtS/ua/dSfAeTpIMBCaEagURY+NwwRcUuuGG354uyXp7s8ORgx4LuIH71vkGc/zrVIFk560e9n33Vj/qd7d/VeqGY/SJ9cXOaqXqTNhbYFMq1jaZedqYFKhWVHp76cWOmDfLHMyLBAHriZq9Hm6p+qgQqI0JXWjYPyJz6cRo0u039iIE/G4oI6K9rRysqf/rTn/DCCy/gsMMOw0c/+lEsWLAAH//4xwGABhNCj1ITDRppGfwgSoGKJUwxMR+sq7XQZ6mfeFgSTgj+CLj3TOYAAK8dSLd5Jf6D91BxcAAVK+bsSpMZUdFMy4YSepyfN7fQr9lHxVSeXK0XivmaqD/O/ndjabPdy8YUFaB2+qda9VE9GBSVcNCdomJRnswY7IpUnI+akfpJVRmS2eXAb2PuSgtUBj1OAwZ20eU4UBGaTE1krRSVDg5UAGDx4sW44YYbsG3bNvzkJz/ByMgIQqEQLrzwQnzqU5/C008/3ax1+h6veqgwqJdKdfTyZOPtVVM/RV1RYZ9F2ScvL3ufJ7LORt3PJviJycFni7fQF8y0Zg8Dw7o82duDNDv3sBk8Nc20pvJkp31UwgG9isfsUbFqFgeoQRS7rqqV/vGyjwpgDB4TEbeBivF9FY+5Zn8K0KzUj/o5tVJU9PJk+8+yuSstUJlGcho0c6+Vw0AlJygq42lRUfHv9OS6V3T22WfjF7/4BXbv3o2rrroK99xzD44//ngv19ZRsNRPo0ZaBg0mrI5dmaduprXoo8JSP2GJexb8Egeyg6iYMyZU3MwgEQP8XI3Uj2V5sscHaXPqp1aZb0Vn2iqpH/FQEw9LQnM59tja83lY5U9tRUVbn0cvj0FRiYb41GFHqR9TY7JYxKiomGlG6ietRZ5WHhV221SjiorDF1sPuJ39fwaPinZhVCrLPCju6NSPHf39/bjqqquwfv16PPHEE16sqSPJeNQ+n0GDCatj51Gp1kJfN9PqqR+fZH50RSVDiooZN1d6opm2ZupHVFS0dInXikpl6ke9vZZHhad+qlTbiIbcZESqKIWuFuSIjwNqlyjzdXiU5hdfZzepH0VRKvYHsUeOlaISbkLqp5pHZV6PuoZ9kzlbLyfroSKuV0wjhQKSbVm5mZjLqp+cwaOiXhiJalPHBirbt293tLFjjz0WALBr164a95x55Cj101LsPCrVpidzj0pEr8zyi0dlJgYq6XwJr4xMN7wdp3N+AOPnhptpbVI/hvJkjxq+mWHnnoIp9WPfmZbtl+rv1XqhiOcTg6LisDwZEAYT1kj9eK2oiO9JXOijUk2FUNeh8IsL3Uyrb8sqUGFN9zzto6IpP1YeleG+uHqffAlTWev/x0pREdUZN8ECD1Qcqu8Gj4p2vBENth3bR+X444/H5Zdfjscff9z2PpOTk7j99tuxatUq3HnnnZ4tsFPwqn0+g9ro26MoShWPivrd6mXL8kBFqqiuaDfsimZ8BqV+rv7lMzjrloewYedkQ9upR1EpCuXJVl1pxe0ZypOblfoxV+LYmmnV7yzIYN8tPSrCv5UIS7ofhptpaxtgnQ4mZKZer14eg6JiKE+uHqiLx0P2Xhk9KpWpn1CAqWzeXJQoisIVFavUTyISwkBSXcfOiYzlNvZbKCrRUIBfeLnZD1nQl3ec+hE8KiZFxUllXTtw1Jl248aNuOmmm3DuueciHA5j9erVGB4eRiwWw/j4OF588UW88MILWL16Nf7v//2/OO+885q9bt/h1eRkBk1QtkcUQew601YtTw5LFV6AdsMUlckZYqbdP53HnzeqpdbP7ZrAEYt6696Wm2Fphqofh4pKrlgWqn68rWCsUDlqmFLN+yUze1stS0z9JMJi6ke9rVppM8PpYEKv+6iIxuiEUJ7MKrDsTtQFi0BFvDic210t9ePNhz1X1P0cVqkfAFjUH8dYuoBd41kcPly574+n1QBB9NRIkoRkNITJbNGVsscCcbOiMpEp4PO/fxHvXL0Ip6yYI6y/UlHxc8UP4FBRGRgYwH//939j9+7d+OY3v4lDDz0U+/fvx+bNmwEAl1xyCZ566in89a9/rTtIWbNmDSRJwtVXX13X49uNXvXjzRvNrjjITFuJ+JLYTU+umvoJS0LqpzlrdAsLSGeKovKnF/bxg/k+rfS6XnKl6gGHiKGPSo2DL7vqPZAuNK/qRzuH5ksKFEUR0jHVy5OZguF0KKFV6qdaszgGV1RqmGm9rvoxz/oRL/Cqldmy42E4qH+Gxce2IvWTyqsnd0myb0exUEv/sEGJZsa0QKU/YVSAWODjZj/UA27j//eLx3fgzvW7cPvDWw23ixe/rDyZt8/3aaDiatZPLBbDxRdfjIsvvtjTRTzxxBP4zne+gyOPPNLT7bYSrz0q7KAsXkEUyzI++7sXcNJBg3jbUcOePE8nIs7xqeijon3OzPGdrCiCR0WqMC22m4K24FxRrVbxKoXYLv64YQ//ee9UY4FKPYpKoaxX/di9luzqeySV92zWjxmxEkc8j9i30Fe/syCj2qwfUf1JGKp+TOXJVUSQLoeDCb3uoyIGnclICOGghICk/r+5Yhm9cesBXFb7Qs3y5ABLB3rzWecDCSMh2x5iLFDZNVEZqBTLMvfiDCaNgQpLJbnzqOgpTJGHXh4BoBd6MAyKilaezJu9+TDtA3hQ9dMo09PTuOSSS3D77bejv7+/3cupG16e7FHqR6xeYDz40gh+9vft+NqfN3vyHJ2KeGFkHocRsumjkheuGGPhgO8UFfF97nRD7Vi6gL9tPcB/3zuVb2h7XBlxoKhY9VGxU1SGWKAylWtaeXJMUywU6B4pwPlQwmrBRtisqAjqDeAsuEg4rvrR1u3VUEKToiJJkqOZNVbvaS2PitflydNVBhIyFvVrgYqFosJU04AE9JgCsnoUFRa0iYrKdL6Ep14bB1BpHxBTRKl8CcWy7OvJyYAPApUrrrgC559/Ps4666ya983n85iamjJ8+YVmNXwTDVJPbVd3PDfjvGcirFInIFUOa7PzqDB/SkACIqHKVuXtpmAhx3Yqa1/ci7KscLWr0dRPXeXJDjrTDnWrbe5HU3l9IJvHV5Qx4cwu+kBszbSSqTNtlV4oITtFRTs8VGsWx0i2qY+K+F6y9AlTvpykfsTHs8cFA1JFKgXwPvXDFRUbfwoALOxXZz0xM206X0KmoD6ONVnrS0Qq9gNWRVSPoiLaBP625QA/BporR80poslskRSVatxxxx14+umnsWbNGkf3X7NmDXp7e/nX4sWLm7xC53gdqLCDq3il/bRNhDzbYMcxaznceB8GS/vEwhIkyV8t9MVmS4CxW2Qn8scNewEAbz1iAYDGUz+5elI/JZkfuO060w5p/S4OpPXA0OsDdTgAsFORqFrYXbhWmGmrpX6E2yxTPw6CC1ae7Lzqx6PUT5VApVrjMl1d0/cF9riBZMSmjNvj1I8DRYWnfsazKJVlXHjbX3H2LQ8jXyoL/pTK9BbbphtTt96ZVn/dWNoHqPQ5mn+fyBSFycn+TDm3LVDZsWMHPv7xj+OnP/0pYrGYo8dcf/31mJyc5F87duxo8iqd46Z7phPMikqhJONZrcyzMMsNttW6e9qlfsTSZAC+Sv2YW3tPdrii8uIeVel892r1QmIyW3TcjMqKWgGHiEFRqRHgDCQilR4nr3IbGpIk8fTPtBAM2J3vzamfakqGJOnVa/FwZcM3Xp5cpUkbq/ppZ2faOA9U1NuqKcbmyckAsHJBDxb2xfHWVfMtHxMW0oFeUK3ZG2OhlvoZzxTx+LYxvDIyjV0TWewaz/LUz0CyUv3pbkBRYSkdRVHw0Muj/O/mC1tzGfNEpuD7qh9XZloveeqppzAyMoLjjjuO31Yul/Hwww/j1ltvRT6fR9BkQIhGo4hGK81SfsDriJR3pi2rO98Luyf5c8x2RYVfKVq81GIfFUVRuNmNdaVlJw0/pX7M0ux4h3tUmDS+ZCCBeDiIbLGMvZM5LJuTrGt79ZlpKzuYmgkEJMzpihoUH6/NtAAQCwHZoq5ahAL2g1y5mZbP+lEMt5sJBYFySVVG9LlCxi64drN+AL3qZ7rmrB9v+6gYZ/2oC2cBS7XGZVbVKb2JMP7y/51h+5rq6eDWBSq98TC6YyGkciX8+qmd/PZ9U3nbih9AT/3UVZ6sXQy8eiCDHWO6N6Yy9VOpqPh5cjJQp6Lyk5/8BKeeeiqGh4fx2muvAQC++tWv4re//a3jbbz5zW/Ghg0b8Mwzz/Cv1atX45JLLsEzzzxTEaT4Hb2zpbfSKNvJmDEKUK/A7VozzwaqTaAV8/bi8Y55VFgXS31KbbNW6RzzgaSTzbRlWeFp0K5YCPN7VbW0kfSPG7VSL08u6wFOFSWGpX8AtdzUq9SGCPOpTOdrn+zZfsn2cb0Nvs1JmJXoWs36YaXQDvqoOO9M61ELfQszLDvhVhuuZ7cv2AUpABDWIjVZ8SbVm3LgUQGARZpP5Z7n9Qq4kVSO91CxUlTYNt0oG1FTyuwvr+wHoLfyr1BUKi6MCvyC2K+KiutVffOb38S1116Lt771rZiYmEBZ+wf7+vrw1a9+1fF2uru7sWrVKsNXMpnE4OAgVq1a5XZZbafgcURq7kz79HY9UFEUY4nubKNU1aOiH7DEnLTYQwUQvQDtfx3NB45ONtOyq00ASEaD+tyTBgKVWmXGIqy8M10oY0rrclrtcUNCg7BwMFD1hFcvPFDRVL1qwZC5PJl3prU5rLDDjZVHRWam8yrPxxu+tbiPCp/TI1TgOTHT6ikK5xeyYiWXF6qKE48KoPtURM/N3skc90T1W6V+6ihPFjssA7p5/YiFaqM5syfFUlGxSKn5Cder+vrXv47bb78dn/70pw2qx+rVq7FhwwZPF9dJeN0wSlRUFEUxKCrs9tlKdY+K/rP4ecwKZlpA8AL44GU0e1QmOthMm9YO4uGghGgoiPk9mqLSQOWPG0VlqDuKvkQYZVnBc5qnq9rj5nbr/rhmVTywAIL5QKo9TUUL/RreEBaYJyw8KiUHVT96C/3qH4RSleqjemAnYpb2AZyaad1f+RsvXhr/wLN9vLumohKvuG3fVF73qFikfs583RCOWtSLC49e6Hg95qGEaa26qE/bPjuHMNj9klqqbSJbaNoICa9w7VHZtm0bjjnmmIrbo9Eo0ul0Q4tZt25dQ49vJ81SVAolGbsmstg3lUcoIPEDRr4kI+lPu07T4TK0xUVVQKvoKctGQ22Op37Ug1bIR2baCnPbDFBUmIQ9z5PUj/POtJIkYdVwL/7yyn6MaIPfql19GxUV79UUQFRUtKZyVRQOPpSQlSfXMMSylyQuzvoxdaZ1UvVTK/XjvaKiDRQ0VO84MNPW0e8jHBAVlcYVVObBshpIKCIGKqyZ3b5UDlPamAwrReWguV347ZWnuVqPXp6svjYZrT6dpZZkTYFn+zcLBOf1xrB1NI3xTJH7ZdwoVa3E9W63fPlyPPPMMxW333PPPVi5cqUXa+pICh53ttRb6Mv8yvDw4R5+gp3Niord5GQGN88JaR1upuVVP+rtfghUzIpKJ5tpWaDCDuILNEWlkdRPrX4oZg4f7jH8Xq1aSPSoNOtqkp3PmEfFiaLCO9OyYMMmtmFBUHfUyqPivOqnVuqnVCMF5RZdUansMFvVTOvCWM0IBPR2+14oKimXqR8AOO2QuQDU5oJMUTF3pa0XcQq4oijIaIpJn1D+LJ4vWOC/QLuImMz4v4+Ka0XluuuuwxVXXIFcLgdFUfD444/jF7/4BdasWYPvfve7zVhjR1DweFZCVNj5DmiTNof74tg8Mo1SoWwbqGwZncaTr47hncct9sz45jdqXd2FggCKRjOt7lFRHxTgZtr2e1TM7+VkJwcqJqMhN9M2kPrRPSrOPluHLzQOgauuqOipn2ZU/AB6MCFW/djBvVPmzrQ2j7nylC48tr2AoxaEsWlUfe11j4q2zaqKir42sUrOTLV0az2smJtEQAIOm9/Nb+MpDAtFpViWEQ4G6kr9AOrFS1lWPAlUWOO2ZKRGoCIoKu84ZhgPvzyKvVM5/r5YKSr1IHqw8iUZGS2QEquKRAWeKyraRcR4puD7qh/Xgco//dM/oVQq4ZOf/CQymQze9773YeHChfja176G97znPc1YY0fg9awQcSghm6jbEwsjEgogUyjbDiv87O9ewCOb92O4L443aFH8TKNcozW4ejBVDDKv2aPCPo9+EKYqy5M7N/WTNqd+uKJSfxt9t4rKKpOiUu2kJqZ+mnU16cZMy/4kK2p5fbXOtABw8eFxXHy4ekI0DyVkKki12IJ5VEoycM0fJvH3HQX85pJBDPcYX2uv+6gcNLcLj13/ZkPli16ebPw8fP3Pm3Hbulfw6389pe5+VeFAADnInqR+0lpqxW4gIePQed1YMTeJ4b44Vi8dAKB+Dtj7YeVRqYeY8FrkimXuUemKqjOUikKpPqArVsw/JjZ886uZtq4+Kpdffjkuv/xy7N+/H7IsY2hoyOt1dRxev9HiUEI2wKo3Eebbt+ulMqrl5fc02Lbcz7B/3e68pXenFVM/Ro8KV1T8EKholXODyQgOpAuYyBarXt36GXPqhykq+6ZykGWlagWKHW4rEpYNJpGMBJEu1Pa2iKmfZisqzKNSPfWjvz6yUrvqR4R7VHhn2toqCFNUAODuF9Vjxm83ZvHRE7sM9/O6jwoADPUYG33GTNUrjN8+uxu5ooz1O8YtO9M6IRwKAHlvUj/MQ5OooajEwkH8+RNvgqIoPL0rXpT0J60HL7olFFQrp8qyglxR5utLRoOIBAMolo0KPPs8sdTPRKbg++nJrle1bds2bN6sDsWbM2cOD1I2b96MV1991dPFdRJ6HxVvFZVCWeapgN542DAd1gp2oujk9EEtmEfFTlEJ8+60+m08UDEpKn4oT2YHETbNt1CSq5Zo+hlz6ebcrigCkmrmE1vVuyHnUu4PBCSsFFSVakrMnK4oWDzYrIN0ZeqntqICqEEKr/pxEN/pVT+qGuNk1k8ooJtwGQ9sqVS/uKLSxOCZDXQVzbSpXBFbRqcBAFPZIu+x4lb98rLpW6ao7uOJqLNgSZLUCjhRPQoHpZp9WNwgBnksQI+HQwYLAYMFJUztnMgWPS8G8RrXq7rsssvw6KOPVtz+97//HZdddpkXa+pImtZHpSgLqZ+QoUW4FWykdydXjtSitkeFHZQq+6gwMy07IfjBTMsOIv2JCD+gik3fXt6X6phBlDz1o11thoIBzOlqrJcKuwJ00keFcfiw7lOp5m0JBwNcgm9a1U/YhaIiBiqy0JnWkaKiPlgBUJSFZnE1Hjvcrb6uHzlR7Rz81K4ixrPGD4bXVT9W8A6rwrFtw65JMBvZVK5Ud+rHy3k/GYepHzNimrE/EfFUMRUNtRmuagZ1BV6oLGQeFaZ2ZgplfoExY6p+1q9fj1NPPbXi9pNOOsmyGmi2UPC8M60+lJA1ruqJh/ntdqkftsN1cnfTWtQaX8/b6BsCFfX14g3ffNiZNhIK8N4HzKfy1GtjOOcrD+PTd3VGjyKriohGDbX1lKQe7lBRAfQURNNTP3ntGFHlaUS1pSwrNTvTikSFFrT5kqKnjWo89PaL+/Hzf+zHf7yxG6+bG4KsAA9tNaoqJaX6Z84LeMM3IShnFY+Aqq4U6tgXAP3CzxszrRaohN0pIuxzAFh3pW0EsZcKq/pJREKCAi94VLS/z+mK8veTXUTMGEVFkiSkUqmK2ycnJ3mX2tmGmIP0vjOtYKYVUz8WgUqxLPPbJ7IzN1DhV3c25x9dUdFvY6mfmNY3wldmWmHf6YureWuWutsyqvYm2ri38jPnR9L5yh4TTGLeU6+iwlM/9SkqtU5q7Eq3aeXJ2rL1hm9VUj/CEsqCR8XJ0sRMRL6k1DSdMw4eDOGUpeprcOYK9fv9pvRPKxQVNt5CNH4+u2OC/zyVLbk2VjO8Sv3IwogIp6kfxjyhwsxqzk8jsH08VywbFB+zolIqy/xCLx4O8nW8vE9Nr/nVTOt6VW94wxuwZs0aQ1BSLpexZs0anHaau0Y1M4WyrHB5shmdaZmi0lsjUGE7KADeVGgm4rSPijhmwNxCXzfTtl9SERUVduBggSY78TOTtN9hFRFi107WT2LneKaubdbT3vuQeV0YTEbQFQ2hJ1bdtMgDlWZ1pjWlfpyUJwNsNo36sxMhQ5IksGxEvgzHiorIm7VA5aFteZ4mURTF86ofK1jqR1RUDIFKrlh3ebJXqR/RO+Y29TNPMG57rqhor910vsQvfJKRUIUCLyrxsXAQbztqAQDwQYl+VVRcu3m+9KUv4fTTT8dhhx2GN7zhDQCARx55BFNTU3jggQc8X2AnIBpbPeujIgQk7IPbE9OrfgoW6hUrSwNmduqndltx9buoqLBAJerj6cnRYAC9WpMmlvphgcpYOo+yrPi+N07Komvn4gF1ONtOYaKrU2RZVyvdeFTCwQB+e+WpKMtKzcexyp/W9VFxaKaVFb1pm8P3PRqStMnRSs3PiRVHLwhjIC5hLKvgyV0FnLwkakiPNmNoIyNmmp48kspht5AunMoW0aMpjm5TP2zeT7FBCTUjBFExl6rOvCamftjrMSYY1uORYMXMOLGiKhoK4ONvPgR3Pr2Lq/YzRlFZuXIlnnvuObz73e/GyMgIUqkUPvjBD+Kll17qyGGCXlAs6Z9kz6p+tB0mWyzzvH9vPGxpjmKkhYFwM9lMy6cn21wqWplpC9p7FAn5N1ARUz8s0GQOflkxHoT8Slow8jGWaIHK9jH3iop4EeD2ILqoP4Glg8ma91sxVy3FndvdnJkULFCRHaRxJEky9FKp1ZnWjDiY0G2Qw+57/CItHbBffS/Fc3srFZXndqj+FGbPmcqV6upMCwARj1I/emly0HWpvSH10yRFhR0jwkEJkVBAaGdR1r7rHWgDAQl9iQiuPusQvh2/Bip11UcNDw/jpptu8notHUteUze8HBPPPojigK6eeKhqeXJaiPZnsqJSbXoyoL8HLEUkC9I1O2Cxg7cHRQANI3pUWPqBpfvE4HM0lW/aydQrmJm7OyYqKmrqZ0cdqR8xIG/WQfRtRw2jKxrCicsHm7L9WMh4TKh1sg9Keg8V3rTN4b8e5akfxXWQw+iLqU+Wyhv7sQBAqImCHm/4pr3nz+2cAACsGu7Fhl2TSOWK6NH2K/edaau3dXAKU63dpn0A3asFAAMJb3qoMFhl25imxLIeL+a+W0xRERWp95+0FD/7+3a8MjKNBb2VgxT9QF2BysTEBB5//HGMjIxANnXM+uAHP+jJwjoJsSutVyVn5hRSLBxANBSsWp4sntRSuRJKZRkhn85uaIRajaz0qh/1e1F/WcB6NOkzVdofqRSEq5xuzU/BUijTwnu6f9r/PhWuqAjNsBb3q4rKRKaIqVyxpmdEhKUBggGpaftyOBjAOYfPb8q2gcpApdbFTCAAQCsvZvtnyOFxRZz3w9vvuzwkdUfVB7BARRCMm1z1Y7z6f0ar+DntkDnYsGsSU9kSBpN1Kio1GmU6haV+4vUEKr1CeXKTqn7GNUWFTUY2expzFqX+4WAAv/rXk/HSnikctbjP03V5hetA5fe//z0uueQSpNNpdHd3G07MkiTNykBF9Bh4hTlQ6dVSAtU+cGKgAqhSqde5UD9Qy9gXNqV+CsKRVk/9aIqKD1I/4pwNVtbLZuaYFRW/Y1WenIyGeNfdHWMZQ0VOLXgPFZ9K0k4wz66r9a+oTdUUyIpSc9aPGRao5EpCwzeXwUV3VFP1tHJq8TPS1KofU3nyTi1VeNySfgDGVg1uvYC1+k85hXd9rdGV1orBZJR3kG1WeTJL/bBAymymtWueOJCM4JSD53i6Ji9xvdt94hOfwIc+9CGkUilMTExgfHycf42NjTVjjb6H7fxhDz/F5h2JXYVW66OSMTUFm+jgmTHVKNfwqNgFKsGAftDmw998EKiIqR9WLcOUlLRQyeVWUdkxlsH1dz7HO3u2AvOsH8Yizaeyw6Whlld5uGyZ7iei5tRPDXWEBeAl2X3lDutOmxWmIbs9LFUoKoLq2EwvN+8Fok0BZkHJcF+cPy/7DLhN/XjVR4WlfupRVIIBCUu1z4E4WdkLYiYzLTOzR8welTqaJ/oB12fWXbt24WMf+xgSiUQz1tORsJ3fy/LGUECCeDwzKypWVwbTJkVlpvZSYSmdqtOToVf9sCFtEeFoLzZ8U9o8QblQTVEp1K+o/PrJHfjF4zvw879v92il1SnLCg+Wk6ZAZQkPVNz5VOrtROonzKmfWucItpuqZlrnnWkBfd5Pplh/cKEHKkZFJRRAU+dPsZOnOuVYwVRWn3HGUqLs4qPehm9eKSr1eFQA4Nb3HYtvvf9YHDS3q/adXRA1mWkTXFExpX5K7iaR+wXXq33LW96CJ598shlr6VjyXFHx7kOszofQ354eB4FKpmAMVCZnaKBSc3oyn/VjVFQiITFQ0e/fbpuK6FFhSoSlmdalosLM1alca/YDMagyKyqL++sz1HLz3wwKVGqaaXlaUnE9Y8dcCq0+1uFCNfTUD1NU2Lrcbcct4slzIlPgSmNvPIyeuHF/cutRMZ+w6yXjcCChHSuHe3DuqgUNrcEKFrgdSJvNtKY+Kvzz1FmKiutX+/zzz8d1112HF198EUcccQTCYaMx7u1vf7tni+sUmqGosO0x8xNTVMzlZiLTeeNtM3UwYUm4wrPCbKbVAxX9PuJb1e70jzh5m1XLeJH6YfuIOSXYLFhQpQ66M745i+ssUWYH2E6TqkXYrB9GLWMsn0OlQPCZOHsulmYSUz9uPSo9FYqKO0NvvUSCAUgSoCjAvil1Xw8GJCQjQS31racNXad+PKr6yTRQ9dNMWHkyM+EnTGZalvLRzbSdFfi7DlQuv/xyAMDnP//5ir9JkjQr2+jr0r23O6+6PXXHY2V51TvTmlI/M9SjYqWQiJg9Kqzqx6ioiDNVmrFK5xg8KqZAZboBMy3bR1oVqLB0VVcsVJEiqD/10/mKinlIrpPyZEDtoVKvRyVjCFScPZbBFJVUixUVSZIQDweRKZQxklIbvfVo+5K5Uqxej4p3iorPAhVTIJ80lSezBqHs8+S2WV27cR2omMuRCVFR8faKQ/wwco9KlSuDdGF2eFRqBSrmhm81Uz/tDlREj0pUfZ+ncyUoimJI/eyfdhd46oFKqcY9vWHaojSZwUqUd4xnIcuK42ZZuTobfPmJSo+KMzNtWdHVDKevF1Nv0gV9p3af+jGaaVnjuFbEijEtUGGKCkt5d5tKp9yaq70uT6439dMszAoJr/oJ2ykqnfV56tzLFB8hnmi8xMqjUr0zrdEoNVObvuVrKirqd2amLViYaY2tyr1foxt0j0qQm2lLsmomFCsuxtIFV1ULLJjNtiz1oz6P2Z8CAAv6YggGJBRKsiuvjV7107mHqsqqn+r3Z5OSy4rQBt9hsNEdUe84nlX3GwnuDbBMUckUFZREn0wLxjewEmU2zZcpKez4x3Cf+tGn0TeCb1M/ZkVF6/xnvrDtVM9XXWFhOp3GQw89hO3bt6NQMF7lfexjH/NkYZ0E2wm8nhUSsQxU7D9w7Op7uC+OraPpGWumZQqJ+QTACJmGElopKpIkIRhQg5R2d6fNl3WPSjIS5Hn6vcK0YdZ/4cB0wTAuvhotT/3k1f2ty9w4BOpnY0FvDDvHs9gxljF06axGvS3T/URlw7fq97dM/Tg8tPRoXWXHs7Kj57KCKSoAMJ1XWjI5mcECUpb6YUqymPqppwO496kffykq5sCDm2m1AIZ9jngVXYcpKq5f7fXr1+Otb30rMpkM0uk0BgYGsH//fiQSCQwNDc3OQKUVikqsdtUPS/0s1AKVGe9RCdspKsaqHzsFJigBZfgr9SNJErqiIaRyJR6oxMIB9MTCGEnlsX867zhQybc8ULEuTWYs7k9g53gW28cyWL1swNE29QNrZ10BioQDqoLHZ/04TP0UBTXNadUPM8KOaYGK00ZxIpGghGhI7W6byss84G+tosJSPyHDd0A9LrpVifxSntwszIqKuTyZKZNMUek0M63r1V5zzTV429vehrGxMcTjcTz22GN47bXXcNxxx+G///u/m7FG3yO20PeSiJVHpVqgop0oWDOhmepRYbYNu2NFRerHouoH8E932oJ2EGHvLWv6tk+bHNsVDfEZP24MtWwfEUfTN5NprQy62yZQqWc4YadK1SKSJBlUlVpdDNh+KcaXTg8tvUxRyajvvdMAx4xYoswVlebHKfyEa1ZUugVFpR51zatApZGGb83EHHhUmmlnmUflmWeewSc+8QkEg0EEg0Hk83ksXrwYX/rSl/CpT32qGWv0PeYTjVeIH0izmdaqPJl9iIa1QGXGp35sFBUmC+tmWvX2CkWFzftpdx+VslGRY6kTpqgkoyHM6dICFRf+DrZd82iFZpHmzd6sD4JszL2bYGsmlCcDxjb6tZQJ9mdx+rdTMYMpKuM5LVCp85DUE9ENta2q+gEsFBXmUYkZFRW3RD0rT66/hX4zMVfxJKImRYWZaTu0is71asPhMJfd5s2bh+3b1a6Xvb29/OfZRpGbNZvpUWE5R3v3esakqMzUPiq1y5PV7yz1w9+fikDFL4qKsQ8PM6Pu0RSVZKQ+RUXMS5dbEI2xHg52qZ9+bWKsG5P3TChPBow+lVqCQJAHKvptTj0ZzKMymVMM23JLN5+gLAtVP82XVJgywHoGMW+eaKatJw3ImnHO1NSP2XOSsJn106kt9F2HhccccwyefPJJHHrooTjjjDNwww03YP/+/fjJT36CI444ohlr9D2FJjZ8YzgqT84bFZWJbBGKojS17XU7yBedmWkrUz/WiopfAhV2MmYy9whXVIK6ouIm9SP8Y9li2bIax0vY/meX+ulPqIPYxl14p2aCmRYwBSoOZ/0U6lBUemPWwbhbeNO3giIM8qxrU65gJ1w21aLHwkxbV+qHVf3MktQPM9NWTE/mfVQ6K/B3vdqbbroJCxaoLYC/8IUvYHBwEB/96EcxMjKCb3/7254vsBPQhxJ63EdF2/kCki412uVaFUUxmGkBtQeDef7PTKBgo5AwzGZaFqiY+1cwo2HbUz+l2qkfpqi46U4r7iOt6KXCAhU7RaVPU1TGBUVlMlvk82ys6FSp2oxdDx8rWHmy2IHAcdVP1HjHekUQcd4P70zbwtQPg6V8zGZat/A+Kg1elTQyPbmZmIM3s0dFH0rYmUM+Xb/aq1ev5j/PnTsXf/zjHz1dUCeiKyoed6bVjk498TBv+GQ3PTlXlPkJd6ArgmgogHxJxkSmaDCidTqlsm7us1dU1O/mhm/m+4f8kvoxeVS4mdbgUVHViHrMtEBreqmk8npnWiuYosKq0baMTuMtX3kYFx+7EF9651GWj+lUqdqMqKjUHErIFJWS+3k9PVHzPu7scWa6I3p32lb2UTErA9aKSv2BihtFRR2OKBv2PeZR8buiErdL/XDPV2cF/q5Xe+aZZ2JiYqLi9qmpKZx55plerKnjKDZZURE/pHYfOLErbSIc5FevM81QKx68ayoqsvEx5osgrqi0MVCRtUmxQKVHhXWi7RI9KnWYaYHWlCgzRcUuxWRO/Ty3cwIlWcEzOyZstzkTpicDxkClViUO96jI7pu2JSOSIahpXFFpbR8Vs6Ji1UelvtQPO246+xyUyjLO/spDuPDWv1oqk37zqNg2fDOnflh5coelUl3veuvWrato8gYAuVwOjzzyiCeL6jTYCSHquUdF3Zl645VXExWBSl7/AAUCEvri7Op1ZgUqrCeKBPsrU/Osn7xN1Q9XVNqY+hGDiYjJo8JIRkMYTKqBynjaub/DeIBtfqAyXSNQ6Uuq/1euKCNXLGN/Sv1fqgXTM6EzLQCIb2nNhm+sPLnk7P4i6lwc65lWbtDLk4U+Ki3wuplPuCxAEVW6eqorI6Yy3VpMZIvYOprGpn0pPPTyKAB1VAo7pvgt9VPZR8Wc+uns8mTHr/Zzzz3Hf37xxRexd+9e/nu5XMa9996LhQsXeru6DqHY5M60Yn42Yso5MtKmZlu9rMIiO7OavonGWLurTFa2zEy3RYsW+oB+5drO1I9VoGJOnXRFg7xiZjJbRFlWHJ2AWu1Rma7hUemOhniH3YlMkfttqgUqM2HWD2BSVBx2pi3wOT/unqs3GsB4Vj0e1CvyWikqrTi3mU+g7CItGJDQHQ0hlS/VV57sMvWTE0qu7n5mF85eOc8Q7Psu9VPRmdY068fU8K3TFErHgcrRRx8NSVJPDlYpnng8jq9//eueLq5TyJvMkF7BdiYrRUVWVHkypB31WOonGTGqMDNNUdF7otjfJx7Wq36KJcW26scPqR/xwMnkaXPVTCIaQp+WNpEVYCpbRH8yUnW7sqz4LvUjSRL64mEcSBcwninwNFauKCNfKlsGI3zaa6crKqJHpVYfFZNHxa2SISoq9XSmBYxmWpZtbUXVjzlQEYcRdse0QKWOiMmtRyUnOJnvf3EfUrkiD/ZDAcnzY32jhIIBhAISSrKCSDDAL5rZMaVYViDLysxvob9t2zYoioKDDjoIjz/+OObOncv/FolEMDQ0hKDHZtJOoVmdaVmr9MVaR0/AGAwVxEDFdDXLThatGkjXKvI15vwAQDyi/y1btA9U/GCmFXuoMIXIrKgkoyFtsnII0/kSxjOFmoGKWeJuxX4wnaseqABq5Q8LVA4I06Ans0UMdVsFKrNRUTFW/bjthSJW/njTmVar+mlB6kf0qERDAUPg0hMPY/dkrj4zbbB+RSVfknHfC/twzJI+AP7zpzBi4SCm8yXe7A0wBiSFstyxLfQdBypLly4FAMjtHoziQ5rVmfYfjl2EeT1RnLB8kN8m9lYplGRoF9p66oeVMXvUidFv1Gr2BqgycTQkIV9SkC3IenmyjaLSVo+KhRpnHmnfpR14+pNhHqjU3K7pfW+2oiLLitCZ1v6wohpq04bUD6CqREPdlTOM9NRPZx1YzYiBda3GaXrDNzZjx91zsTb6QP2KSo+oqLSwM614Au01TUxmfpWGqn4cHg/NqfXfPrMLr5vfDcB/AwkZ0VAA03m1mIIhni/yRTFQ8WewZYfrd/xHP/oR/vCHP/DfP/nJT6Kvrw+nnHIKXnvtNU8X1yk0szPtma+bZ7hCDQUD3J8gliiz1A9vncxyky2a89IqnCgqgK6qZAuComL6bLLXsVofj2ZjLk0GKhUJFnzyqpl07XSe+crR7FH525YD+NZDW3DTHzfi3uf3olGmhe2bAy2RPqHyRwxU7FKUKW1+ULXgpxMwNHyr1UdF+zv76Lo1xIolyvV7VPTyZL2PSmsVlR5zoMK6czcw64elQGrBAmTWx+Wvr+zHznF1RpWfFRVATRUzwkEJTAjLl8sdW0VXV8O3eFxtKPa3v/0Nt956K770pS9hzpw5uOaaa1xt65vf/CaOPPJI9PT0oKenByeffDLuuecet0tqO82anmyHlYxpTv2Y3d4zBSeKCgAktEAlU1D4cDdzcOOHzrTm9vmAlaJiDFTGnCgqJfvUz97JHN733cdw8z0v4TsPb8XHfrHecnaUG1j7/EgwUPVqjZmCx9OVqR8ziqJgROsbM6SVZ3cqYpzlVFFhjQ3dxgc9BkWl3tRPe2b9iPuOWVFh1XD1VICZU+a1YMrDsjlJLJ+ThKwAT702DgCG1IqfYK9LUgikJEkyzPvp1NlZri9TduzYgYMPPhgAcPfdd+Od73wnPvzhD+PUU0/Fm970JlfbWrRoEW6++Wa+vR/96Ee48MILsX79ehx++OFul9Y2Ck2q+rEjEgogWywbghB9WJaxfn7GBio2AwkZuqIi2wY3XFFpY+rHyojdFTUeoNkV0kDS2DCtGhWKiqCsjaRyUBR1X0kXyiiUZaTz1mZWpzDlw67ZG4N5a149kOFlr4B1oDKZLfL/Y6inswMVNx4V1pmWxZauPSpieXKDikqmqCBdaF1nWjEI6THtS0sHVa/ecG9lirAWhpS5qYmbFbyMNxTEvKEYtu1P49kdkwCARNif6h7rjWKuSIoEA8gVZUzl9M9YpwUqrne9rq4uHDhwAADwpz/9CWeddRYAIBaLIZvNutrW2972Nrz1rW/FoYceikMPPRT/9V//ha6uLjz22GNul9VWWq6oWJQoVyoq1h1sO4myrFQM08vzNE4tRUU70BbszbR+UFTyFv4mq/JkQG9BP+Yk9WP2qAijFNhBeF5PjHsCGp2wzBSVamkfQP8fNo9MG263ClSYmtKXCM8oM22twwT7e6FOj4popq27M62QPvrLa+r7cPBg80/Q1VI/H3njCvz8X07Ee05Y4nq7YqBSdHBM5GW84QAOGeoCADy/Ww1U/FaazIhxRcX4PjFDrfgZ67RZP673vLPPPhv/8i//gmOOOQYvv/wyzj//fADACy+8gGXLltW9kHK5jF//+tdIp9M4+eSTLe+Tz+eRzwsGvKmpup/PS/Q+Kq0Z/mfVE4AHKjYzHjqRz901jpGpEr72/rm8Nwpvh19DUYlZeVTMgYof+qhYpH6SkSAkSR/MxoLPgYRzRSVfNHtU9P0gJ8z7SEZCyBULDZttmaJSK1Bh6astDgIVNkJgnoXJttOIuTHTco8KC1Tq96gE6qzUiQQlRENqs8QndqrvzclLqleaeUG11E8sHMQpB8+pa7uBgIRwUEKxrDhL/ZR00+kh89RAhavWPk39WHlUAP1cMJVVzxHBgMSrRTsF16u97bbbcPLJJ2N0dBS/+c1vMDioVqQ89dRTeO973+t6ARs2bEBXVxei0Sg+8pGP4K677sLKlSst77tmzRr09vbyr8WLF7t+vmbQLkXFEKiYKi46waOye7zED8ZmFEXBxt0FjKZkbNuvn8Ts2uGbYR6VtFC1YJv68UOgIuw7kiQZDLXsPe3T0iZjDrrTFsrGwENM/Ygliizfnm6wIRxXVKLV50oxj4p5WKZ1oKL5Uzo87QMY/VGOhxLWmfoRq34aOR+x9I+sAMmwhFXzmj8zzKCoeDyjzE2JstjB9ZChbsPf4j5N/bBjfsKU1mHHlkmt+ad5TEEn4PoV7+vrw6233lpx++c+97m6FnDYYYfhmWeewcTEBH7zm9/g0ksvxUMPPWQZrFx//fW49tpr+e9TU1O+CFaYouJ11Y8dVqXHeuqHeVS8GWveLP76cha33DuJi1cncckp3RV/L5V178iOAyW8boF6kmbdZmtX/bAPp/7/mz+ffkj9WFX9AGrTN3byZyqZrqjUTv2YA1TRTJvVzoDxcJDn2xvts+I89WN9VW6d+lEVFauy5U7DTcM3Fpjk6zbTioqKu8cathORsD+t/nz84nBLFGNj3xRvA4JIKKB6slykfmKhAFbM7TIpnP480euKinF9LG26fUytWprbgcZ0R3vCc889h1WrViEQCBha6Vtx5JFHulpAJBLhZtrVq1fjiSeewNe+9jV8+9vfrrhvNBpFNOq/F7nVioro4makTePH/a6orH1e9TNt2mOtDuSF4YM7xvSrb3Y+dVr1M5kRW9T7z0xbsCkX7IqFgEk1mGDrZGpEPVU/YnmyOJGYKyqeeVSqXwUzjwojFlaNfpMWwdeIpqjMmwGKihi/1VRUGi5PFj0q9QcX3bEAAHURJy1uftoHMF7tm1M/jeKmwCAv9BuJR4JY1B/HjjH1mOVfj4q6rgqPivZ/v3pADVTm93Re4O8oUDn66KOxd+9eDA0N8Vb6iqIf3dnvkiShXG7sykxRFIMPpRNoddUPi5CrKSp+7qMykSnjhV3qyXY0ZX3QMAQqB8STrLs+KhNaoBIKVB7wfaGoWHhUAP2EL1699TdQ9SMqJjmhLT07qLXao8I4aE4XXtwzVUNRmQGBSlj0qFS/r3l6svuGb41X/QBGr8tJLfCnAOaqn+YEKs48KiyYVx9zyFA3D1T8WvUzTwtA5puqotj/vV0LVBbUUTXVbhy94tu2beMt87dt2+bZk3/qU5/Ceeedh8WLFyOVSuGOO+7AunXrcO+993r2HK3AFx6VDqr6eeyVPFcx9qfKlkP2xEBlp6ioOO6jor5GLFCxur8vAhWb1A/zqIiNznjDt0yRXxjU2i6TrK3MtLFQEMWI+np65VExl5SaMSsqB81N2gYq+7ii0nkHVjPuzLTG6d9ug41uQVGptzOtuh1tpEOkNf4UQL36Z/us14pK2JVHxdjB9ZChLjzw0ggA/6Z+/u2MFThmSR/eeOhcw+1MUXntgJrHMwcynYCjQIW1zzf/3Cj79u3DBz7wAezZswe9vb048sgjce+99+Lss8/27DlaQaHVHhWr8mTWmZa10Hc5hKuVPLpZL2OXFWAsLWOuac5LQQhUxtIy0nkZyWjAcWdadgU7la0WqPgn9VMRqGgnfFHGZSf5sqxgKleqeiBn2+2JhTGZLZoCFX0wmawpo417VJiiUv3kEg0FkYgE+XpWzFUrKqpV/cw0M23NzrS84Zv63e28nlhIr9hp5NqpWwv2j18UaUlXWkBV52OhILLFckV5cqO4M9MaA5WDtRJlwL+pn55YGG85fH7F7eyidUq7mJixisrvfvc7xxt8+9vf7vi+3/ve9xzf18/wFvot9qiIH7hsQf2Z5Xj9Wp48ni7jxV2sLbqEdF7ByFS5MlApGqOHHQdKeN1wxHVn2uqBivq9vX1UbFI/puGSgHrQZE3axtMFR4FKX4IFKmIfFT31I0nq9tmsqHphikqthm+AqgxlCmqwumLIOlAxdqXtvAOrGWPDN5edaes4rPRGAxgpyXWXJwPAEfPD+OWGLM5/XWtf/yMW9mLTvhRv8OYVVsdNAHhlZBqLB+KGXj3mGVOHzNMN/35toW+H2f/WiQqlo0DloosuMvxu5VFhNOpR6TTEpmSt7EwLGNM6edNUTL+ZacfTZfzx2QyefjUPBcAh88KIhSVs2FnA6FQZWGi8v5j6AYDtY2qgkncYqDCPClNLrI4tfi1PBsTUj3HhfYkI0oUsxjMFLEPSdrvsfe+Lh/EazKkf3UzLrpTNs4Dc4rTqB1CDp10TWqAyN8nXmyuW+RXsTOpKC5hb6Fe/L9svp/PacaWOw0pPTMJIujGPyvuPSeC8w2IYTLS258bPLj8R+ZJcdQp3PVh5VH7x+HZcf+cGfPj0g/Cpt76e315NUfHrUEI7zIHKgt54m1ZSP472QFmW+def/vQnHH300bjnnnswMTGByclJ/PGPf8Sxxx7bcd4SLygKO33LZ/0Iz22e4eC31M9PH03hzifTeHW/ekJ7y5EJzO1R1zqaqgxuzYEKM9QWHKZ+ElHje2GenAzoEntbpyfblSdrKRRz8ybWRr/WBGX2vvdqvpZ8SeYBNTfThoKIawfdhj0qeWdVP4DRULtsMMnfhylBVWH+lJnQlRZw15mWXe+MpNX38PVD7lMgrPKn0WunVgcpgHrB53WQAlQeE8fTBdx8z0sAgJf3pQz3zZmOp13REG/d34y1NRPzbKQZ61ERufrqq/Gtb30Lp512Gr/tLW95CxKJBD784Q9j48aNni7Q74iKRas603JFRbsyLst6t8UYT/34y0y7b1I9Ob79mATOOyqBoZ6QqqQA/LtI3pz6GTMGKk5n/TCsFJiQD1I/dorKEYt6AABHLuw13N7Hh/pV76XC9oc+IT2ULZbRFQ0ZUj+MVlX9APr/kIgEkYyG0BMPYyJTxES2iCFNlmYVPzOhKy3g0kxr+vPpy9xX3LCmb25Lm2cyZo/KLWtf5inH6ZwxULf6jFx37mF48KVRHLe0vxXL9QwxrRwOShhMtqaCy0tcBypbtmxBb29vxe29vb149dVXvVhTR2FQVNpUniyqJkzm85tHhVXfHH9QDEM96m43pCkqIxaKCjtvdsUkTOcUrqiwdh+10sQJB4GKn1I/UdO+c+br5uGZG86uaJDmVFFhAWpPPCRU/pRMgUqQT9fNeORRqVX1A+iKymCX+r1PC1QmLRSVmZD2AYBEWMLcZADFslKxb5oRY4twoL7SYFZa3KJrp46AX+CVZby0dwo/+/tr/G/mTsl5oTKO8Y5jFuEdxyxqwUq9JSr0phnqjtU9UbuduD6zHn/88bj66quxZ88eftvevXvxiU98AieccIKni+sE2IkmHJSqlot6iVnCzAm9UnigwvqolGSDn6hdTGgydl9S3+WYgbaaorJirnr1PZGRkcrpk5BrzfqJm+ROqyGGLDbwY9UPYN3FVS9Rdpb6iYaC3GDNKnuYRyWumXOBxlI/iqLwA72T1A9TVOZ0qUEIMwWLTd9mUldaQA2K/3DpIO75pzk1B2qKKsjqRREkI+4vgHo0RaUDz0lNQ+zW/btndkNWgIV9ql/DHKiIPq5OR/SodGLFD1BHoPL9738fIyMjWLp0KQ4++GAcfPDBWLJkCfbs2TNjqnjc0Or2+YBFoKKpJiFh2FQ0qH7AFAUotfNMDPVDn9UCj34h5808KvtTZV4my2ABSW8igDld6mN2jZccV/1EQsaDdDVFxY99VOxggUqtCcpiAMTMf6yyR5wMyyoYGilPzhTK3P/iJPXDAhSW1mFlqKKiMjLDFBUAGOoKYkF37ROfGMecvrw+mX4grm4kVuNzMpsQUz8T2r52zJI+ABaBSqky9dOpiOemTvSnAHWkfg4++GA899xzWLt2LV566SUoioKVK1firLPOapmi4Ce4otLCsdnmtE7eIvoXDVT5ktyyiiQrWNonGpIMHToHuwIISEBJVhWXgS59/WJ1z5zuIPZPyxiflh0HKpIkIRGReOWElf+NlX36wqPi8P3pTzKPSg1FRau+iwSFYKSoHozF1A9TW9INBCos7RMMSI4Gnp1/5AK8sHsS7z1hCQBBUTGkfphHZeYEKk4Ry5FPX1bf//+uIxIYmZbxvqO9LfHtZMQLvJSpp0g6XzI0UTRX/XQy4rmgUxWVuuzLkiThnHPOwTnnnOP1ejqOVjd7Ayr7AbDoX5T4xPXkNRNlu2CBSl8iYAhmgwEJg10BjKZkjKTKxkClqKd4+jQVZjxTdtzwDVANtSxQsTbTMkWlfYoT76PisLLFaeqHBa+qoqJuO2NK/cTCQV7+7LY8OVcs42O/WI83HDoXJy0fAKCqKU4uVuZ0RfGldx7Ff7cKVFgPlU7s+dAorMHbnGQArx+q73O7qDeINedWeglnM/y4WS5z8/d8rVS3WFaQL8k8MNE/IzNLUenUz1PnvwttRveotCH1owVJVopKICDxKiQnsy2aiZU/hWHnU2HNrqIhCX1J9T4HpmWuftRSVAB9gjJgU57sB4+Ky9SP4/Lkst6wKm4OVHh5cmVayCmPbT2AP724D7c+sJl3vHSS9rGimqIyk1I/ThnSPidnHxxtqGEbYURUVFiVjzigTxzMydOjM6A0XjTTdmIPFaBORYXQKfITagsDlaCxPFn/UBnXEA0FUSyXDFOW24GoqJgZ6gnixd3FikBFHD7IjLOsxJndXguxusLKwMhervZW/WgpGof7Dy9Ptpg2bNxupaLCfCiG6cn8b+4UlR3ayPh9U3nsn1bVj65ofS3PrQIVts25XZ15BdgIb3t9DHO7Ajhqfmvm68wWRI8KS/30xsN8pMN0voTBrqg6GLdUefHXqYjnhU71qJCi0iCtHkgI6DlHrqiwCg/Th8ov3WnHM+oJsj9R+aGfY9P0TfSisABnZEo9uEgAnBw/xF4qVgoMu1rtKI9KwtkEZXG78bBxQrKYf2eKSqZYhuxCWtoxrs9renH3FADvFJVcscyl997E7DtZBwMSTl0aRVeUDs9eIirRepVaiKfF2W3i8XJGpH5mY9UPYYRV/bQ09RM0NnOzak4EOOulIssK3vudx/AvP3qiaWXM1VI/Q1rqZ8SsqAhlyKxSiCkqkZCzUnAx9VN11k+bUj9bRqexXVMmzN0j7WBDCotlxdDDxwxP/YRFj4p6IM4K+wvzqCiKnhJyAhsZDwAvaIGKkx4qVpgDFdahVpL0mUcE0Sji6JEpzaPSJQYqOaPZHJgpior6PwQkYG6HmtMbOgpks1kUi0YJuqenp6EFdRr5NigqleXJxgFadvezYnQ6j79tPQAA2De1wvO1AsB42j71M6D5TyYzxjWKqR/mUdGNsc6e15D6sXhMqI0N3x7cNIKP/Xw9UvkSFvTGcPyyAUePEye3Zgpl9Mat9zt92KFomC1DURSDohILBXlDuHS+7HiOyY5xMVCZBOCsh4oVTDVhgQr73hMLd2RzKsKf8JR5yaSoxIxjJJiaFwxIba2W9Ap2XpjTFe3Y/8f1qjOZDK688koMDQ2hq6sL/f39hq/Zhq6otO6Aah5KmLcppXPSRl/0Bby8L+3pOhkTWurHKlDpiqmvWypvClSE1E9/wqwUOXuta6Z+2lSeXJYVXPvLZ5DKl3DC8gH87srTHFdlhYMS7/8iXvmZEVOSYuqnWFa4eZh1pjU3hHMCU4IAYM+kanxtNPUzkTEGKn2zMO1DNA923JzMFMHE4+5omKuUKZOiEmvhxWczOXReN5KRIE47eE67l1I3rt+J6667Dg888AC+8Y1vIBqN4rvf/S4+97nPYXh4GD/+8Y+bsUZfU3BZXuoFenkyq+KwVlT07rT2JyAxUNm8b9rTdTKYmbY/WfkadWkdNKdz1g3foiGgJx6AGGY4qfgBgISQ+rEKbtqV+tm4ZwrjmSK6oyH89J9PdCXHSpKERNhYxWOFtZm2ZEjvsFRhwuVgwslMkR/UReoNVETfjaIofJ/sjVOgQngHC1SYUTsUkBALB3RFJW+qipsBaR9ANdA+9Z9n48vvPqr2nX2K60Dl97//Pb7xjW/gne98J0KhEN7whjfgM5/5DG666Sb87Gc/a8YaW87W0Wl895Gt+O0zu2reV+9M23pFRS9Ptv5gmauDrJhqsqIiKwpP61grKuptuaLCK6gAYaZPSEIoKKE7Xj3osEJUVKwUr3bN+nlMS7Udv3ygrpRh3OQ5sYIFp5GgsTyZXS1Kkr5/uO2lIqZ9ROpN/bBApSQrSOVLFKgQTYHt72Nas0TW90c30zIz98yp+GHEwsGObsjq+ig5NjaG5cuXA1D9KGNjYwCA0047DQ8//LC3q2sTT7w6hv/zh4244/EdNe/rtg+GF5gDkHwNRaVaH5Vmp37SOQUs82Tlp0hEJK6WpIX0j3mmT79gxHWqqMTD1VM/wTalfh7bqn5mTlzuzJdihgUeVVM/5UpFJVMs66XJIf3A5baXCkv7sFb4jHoVlXhE75A7kdaHE/ZQoEJ4CDtGs0CFKSl6oFI5YoLwB67fiYMOOohPSV65ciV+9atfAVCVlr6+Pi/X1jaOXNQHAHh+12TNks12NHyLmcuTa3lUqigqYqCydTTteZdWVprcFZMsm64FAxKS2qRXMf0jmmkBoxoTqTGQkBGvYaZth6IiywqeeFULVA4arGsbcRepn2gowA/EqVzJskLM3Lm2FqyHyonLBwwD9OpVVAC9kd1YpsC9KqSoEF7CLvDY7LNure9P0qbqJzYDmr3NFFyfXf/pn/4Jzz77LADg+uuv516Va665Btddd53nC2wHhwx1IRoKIJUvYduB6ipDO1ro8/Jk1vDNpjmRk/JkMVAplhXsm/R0qUKzN/sPve5T0SMGc6t80VDrNM2W8GF58kt7U5jMFpGMBLFquL4KOSeBhRiosBP+VLZoKWuby5drwRSV5XOSfPos0FgpMW9kly5Q6odoCmbVmymA3dyjYqz6mQk9VGYKro8s11xzDf/5jDPOwEsvvYQnn3wSK1aswFFHda5ZRyQUDODw4R48vX0CG3ZOYsXcLtv7FrUTaiuHEnLpv6Q26bLrTGuuDrJiKms8Oe0eU3DYPO/WynqomCt3RLpjAeydLCOlBSqKolQMH+wTjLhRp4pKtEbqR7uplakf5k85btkAn3TtFkepH8FMawhULIyCSW6mdaioaM3eFg/EsWQgwQOXelM/gKCopAvcN0WBCuEldoFKMsLaH7CGbzPLTDsTaLib0pIlS7BkyRIv1uIrjlzUh6e3T+DZnRO46JiFtvcTp9S2CrFJV1b0HdTRmVZUVABg17iXKxUUFYtmbwxeoqwFKiVZn7/DPSqiolJXHxV7M20rA5W/b1MDlZMOqs+fAqCi06wVokelR2ioxkqQxaA2wcy0eYdmWi0wWTyQwOIBfTpvI6kfcdgiL0+mQIXwkMpARd2/urTv0/nK6eKEP3B8dn3ggQewcuVKTE1NVfxtcnIShx9+OB555BFPF9dOjlykTh7dsLN6LoRVqrTSTBsPq026ALWk1Gp6svq78z4qh83rBgDsGmuOR8Wq4odhLlFmagogKir1VP0Ij7Gs+lG/t8qjsncyJxhp6/OnALVTP7KsV1BFggFD51erg7Abj4osK9jFFJX+BJYYApX6r3v6+QwjSv0QzcF8Mcm8W11Rs6JCqR+/4fid+OpXv4rLL7/csvNsb28v/vVf/xW33HKLp4trJyxQeWH3FErVWpW7nNXiBZIkcbk+k9cVFftZP/YnINZK+gStAmW314pK2oFHRZtpMp1n5mD1JBuQ9MGB4pwg76p+NDNtCzwqz+2cwNtv/Qsms0Us7IvjiIW9dW8rHq6e+hGrvMTUT0lW+NRl8SDM9yUHHpV9qRwKZRmhgIQFvTEsHfQoUOGpnyIFKkRTMF/IdfOqH3U/S5sVFTLT+gbHZ9dnn30W5557ru3fzznnHDz11FOeLMoPLJ/ThWQkiGyxjFdG7RuhFdow6wfQr4Kn8/aKipMW+swPcLwWqOxPAdmCd2fuapOTGSz1w8y0opGWldDWo6gkohIWDYSwoC+IRLSKmbbJispEpoBLbv87RlJ5HDqvC3d8+KSGFLhafVREBY2VJ7NxAfum1GZXcYOi4tyjwmb8DPfFEQoGuKIiSXrAUw/MoyKaaak8mfASu9QPS6XrnWmtL/yI9uH4yLJv3z6Ew/YHjlAohNHRUU8W5QeCAQmrFvbi79vG8NzOSbxuvnWFRjumJwOqbDmSyiNTqOZRcZ76WTqQQH8ijPFMEaNTMvoTtg9xxVRWm4JrM5MGsE/9iKZZQ3myw0AlIEn48nsHoSgwlNEyWhWoPL9rCql8CcO9Mfzmo6c05OUAYGjgZoUYmEaCAUiShN54GAfSBeybUtvdG8y0Uect9Pdqjx/uU6ewHjKvC6+b341F/YmG5vJYeVRIUSG8xHyM7jJX/RTMHhVK/fgFx+/EwoULsWHDBtu/P/fcc1iwYIEni/ILTnwqeh+V1nb9YwbIdBVFhbfQd9CZtjcexlCP2sBrLO3dmTujqTNWigaDp35MiooYkCQiEjfROlVUACAUtO7fArQu9bN1v6rIrRzubThIAcBb6NdK/bAgBdBP+kxREQMVFvikHZhp2XMyFSYaCuKej78Bt3/wONf/hwgLVPZO5Xhg3UuzfggPMafn2bRvsY+KOrRz5nWm7XQcBypvfetbccMNNyCXy1X8LZvN4sYbb8QFF1zg6eLazRFa47fndtkHKqyFvjlIaDZJYT6LnaLCPph2nWmLZZnL/b3xMOb1MPndw0BF850ko9XKk7XUj8mjIgYqkiRxn4tTRaUWoqKiKM2LVraOqr14VsxNerI9p4qKuE+yNMpoiikqVh6V2oqK1bYlSWq4PXd/Ul3fTs2oG5CArgZSSQRhpkJR4WZa9XtJVpAvyXoJP3lUfIPjI8FnPvMZ3HnnnTj00ENx5ZVX4rDDDoMkSdi4cSNuu+02lMtlfPrTn27mWlvOUi3/PjpVGZwx2tGZFtCvAkRFxTztU1dUrE9A4pyf7lgI8zRFZTztzUlbURRdUYlUUVS01E/KnPox7Z39iQBGpsreBSrCybWZqsoWzeN0UIsDlYhFoMIUlWiosuonXSjh03dtwL3P78XvrzoNw0IzN0a+SalO5lFhnZF74+GGUkkEYcbWoyIExNN56+7NRHtxHKjMmzcPjz76KD760Y/i+uuv51egkiThLW95C77xjW9g3jwPO4X5AB4MOOxX0Ur0QKVa1U91jwrzAnRFQwgFAxjqZoGKN4pKrqjwAKBa6qc7Zk79qLebUzxnH5GArGRwxKKIJ+sTY8tqlV2NwhSVg6o0DnRDokbDN6tAhaV+RqcrUz9sXzowXcCvn9yJQlnG2hf34dJTllVsO9+kKjeW+jGvlyC8wrzPMm9KICAhGQkiXShjOmevUBPtw5W2unTpUvzxj3/E+Pg4XnnlFSiKgkMOOQT9/f3NWl9bcTJVtm2KilD5oXdSNPdRqV6ePKW53NlJYZ7HHhWmpgSk6r4SpqhkCgpKZcXSTAsAZ7w+jjNeX3mVXy+iwbbUJEklVyxj96Sazlg+xyNFpUbDN96E0BCoqI9hioW4rzCFhnWYBYDHXx2zDFSaZR6PhdXBhNminookCC+RJAmRYIBfXHYJkm0yGlIDFVJUfEldSeD+/n4cf/zxXq/FdzDDYLGsIF8qG+RyxmhKvUI1XxE2m6Qw8ZOX05nWV6s82VwGOtfj1E8mrxtpq3kYkoLaks4rlmbaZiDGll4PY2Rs25+GoqjGvcGkN/tIrdQPuyIUryDNJ36rFvoiT2wbg6IoFe9boWy9r3nBQDKCXRNqUEelyUQziIT0QKVHMLZ3xdQqSrHdAykq/oFCxiow1QJQG6uZyZfKeFUbWnjIPG9kfae4U1RqBCqaBDrP49RPpqBuRxwOaEUwIHEPy3Rerpic3CwCYuqnSYGKmPZp1HDKqJX6yVukIysCFbGFfqTygDySyhsUFkYzy/GZoRYgRYVoDuJ+2yU0KOwSPX82F35E+6BApQqhYICf7NMW6Z+to2nI2tUy83e0ioT2wZrKFnm7dPMHq5ZHxTz8jVX9TGX1FuyNwBSVZBUjLUOcoGweSNgsApLERxE0T1Hx1kgL6M3a7FKS1TwqfBuRSjMtoAaNbJzC49vGKrbNguKmBCqCKkmBCtEMmMqYiAQNqd8urlBT6seP0DtRA5ZisZLZX96XAgAcMq/bs6tlt+say+iVO+YPltPUDzsp9CXCvGW9F6oKV1SqGGkZ3UJ32ryNR6UZsP/Xi8DMCr002TvFrZ7y5KqpHyFXf/TiPpzxuiEAwBOvVgYqVtv2CgpUiGbDeiqZxz0kLQMVUlT8AgUqNUiYRoCLvDKiXi0f2uK0D6CnfsbSeX5bpaJSw0xr8qhIkoRerSPtWLp2T41apHlpcu3djDV9S+UEM22TFRVANfoCzVNUtuzXUj8eGWkBN1U/+v5g9nyI+0o0FOCvwxsPnYsTlqvm+CderRz81MxAZUDw8PRRszeiCTBFpcvU+6BbaPpGDd/8R1sDlTVr1uD4449Hd3c3hoaGcNFFF2HTpk3tXFIF4vA/M1xRGepu6ZoAY0kpoHbGNbeJZwqLbeonV9mqvE87n45Ne6Co5J0rKjz1k5cNs36ajd70zftARVEUbNV6qCz3MPWTCOsm76JFWbXYmZbREzMrKsaGbX2amnH6oXNx3NIBSJJqBB5JGXsINauPCmAMTkhRIZoBC97NHaLFvlR2nj+ifbT1nXjooYdwxRVX4LHHHsPatWtRKpVwzjnnIJ1Ot3NZBliJspVHZbOmqLTaSAvoAdSBtDYN18L4xT0qQgv9XLGM/+//PYc/b9xnOVOlL6EGB2PTHigqeReKipj6sehM2yxYcNeM1M/+6QJSuRIkCVg26F2gEhNeT6v0j9vUDwD8n4tW4TPnvx5HLepFbzzMfSpPv2ZUVZo5LVxUVChQIZoBC7DNqR9mrE0JZlrqTOsf2tqj+t577zX8/oMf/ABDQ0N46qmncPrpp7dpVUZ0j4oxUMmXynhNmyR76Lx2KCrqh4iflCyif+5REa66/7xxBL98cgce3jzKDZ49cX03YIrKAQ9SP07m/DDE7rStKk8GmquobNPSPgv74p7KyJFgAMGAhLKsIFcsV5zUrQyv5rk55vW89QjjnK4lAwm8tDfFA2EGL09uwtWm6FGh8mSiGUSDNoGKoeqHPCp+w1fDNCYn1Zk6AwMDln/P5/PI53VPxtTUVNPXpHtUjCfuraNplGUF3W2o+AGMBkjAupSOXVGXZQWlsoxQMIBdE2pwtWcyx8eaGxUV9bsnZtq8s/JkwNid1q7hWzNoZqDCeuws6I15ul1JkpAIB5HKl6oqKqLq0RUJISDpowJqydqsy7F5oKXembY5fVQYpKgQzYArKlHj/sUClYlMkbcqoNSPf/DNO6EoCq699lqcdtppWLVqleV91qxZg97eXv61ePHipq9L96gYFZXN3Ejb+oofwCJQsfhQicELO8HsntA9B8wgbO1R8U5RSbrxqIiBSivMtFrqpyR730J/PKOqEX1NaAYYi9iXKFuVJwcCkkGlqCVr2/XgIY8K0cmw/bbLpupHVBBJUfEPvglUrrzySjz33HP4xS9+YXuf66+/HpOTk/xrx44dTV9XgntUjCfuzZqRth0VP4CxGR1gfeIRTybs5MU6f4pYeVQOeGmmddJHJcoavinCrJ+Gl1CToLa0ZjR8m9AClf4mVLBUq/zJ20z0Ft/nWgdh9ljz9pvZ8I0UFaLZRGqkfpj5HWhOZRtRH75I/Vx11VX43e9+h4cffhiLFi2yvV80GkU02to0C/eomBQVVvFzcBsqfgC9vT/DSlEJBiSEAhIfXw4AeyYrAxWxIoQrKumyZQt1N+geFeepn8lMmUuurfSolJpgph3Xetw0Y7yC3vStSuqnaqBSI/Vj0yywoPlfmnEQH+qO4ajFfYiHAxXlowThBf1aMDyvx5iOPXH5AAaTEa6oREOBtijlhDVtPRooioKrrroKd911F9atW4fly5e3czmWsNSPWVFhPVQOGWqPohIJBQwDtuyk/GgogFKhzA2WLPWTiAT5SU5MCbA+KoWSWrXDqnHqgad+HCgqCwdCCAaA0ZSMSIgZNjvbo9LM1E+1pm92gYoYkNZUVHhpu3H7zUz9BAMS7v63UwCAThJEU/j4mw/ByuEeXHj0sOH2/mQEP/in4/He7zyGdKFMaR+f0VZt64orrsBPf/pT/PznP0d3dzf27t2LvXv3IputvOpvF0xiT5sUlcms+vvcNhhpGSwtBdhfIXNTZElGtlDGmHbFwD6o0VDA8KGMhCSehmm06ZveR8WZonLkYvWEzmwXrSxPbk7qhykqrU391FJUJKm2ImLnUWlmebK6tuoDLAmiEeb3xvCBk5ZWKNIAcOSiPnz7A6sRDQVw2Pz2KOWENW0NVL75zW9icnISb3rTm7BgwQL+9ctf/rKdyzJgV57MThDxNkbe4tRbuwFa7IRSKMk87ZOMBPGOY9QUm1kCBYC+JOulUr9PRVEUPfXjQFEBgFMOMa6lNWZa9XtzUj9NVFTC9qMdcryPinGfYMqZE1k7ZlP1wxU8qoggZiCnHTIHj13/Zvzkn09o91IIgbanfvyO3rFQPyEoioIsC1QsJs+2iqQjRUWX8NkV/nBfHCcsH8D/feeRWGrRiGwgGcDOMbkhRSVXVHgprJM+KgBwwkExfDswBXYR35LOtE1sod9MRaVa6mda6zhsNgwyRcWJrG03fqHQxPJkgvAD/UnvLyyIxqDLohokLcpAi2WFn9jamcsUS5TtFBV+winK2K1V/CzoiwMA3rV6MU5YXtmzpj+pPqYRRYWpKQHJecDRFQvgqCV6Kq01HhXr8uTJTBFnfnkdbvjt83VvmykqzTjwJcL2qR9Wdt4dtQlUHHTctDPTNtOjQhAEYQUdbWrAcpniUMKccJXpl9SPraLCTjhlGbu11M/CvuoNyAa6tBLlBhSVTF7vSuvGc3DqofraWlr1Y1JUHt48iq2jafzhuT11bbcsK3xEQTMG7MWr9FFhjfzMvSJ0RaX2x97Ko1KW9QCdAhWCIFoF1QDWgKVXRIk9p/0cDEgIB9tn/EsIaaeoTcAUsVBUhnvjVbfLFJXxhhQV511pRU44KIqhniDiEQmtOBcyj4o59fOUNuNmPFOALCu8MZxTprJFsMxmX7y1VT8sqDaX+LpK/YQr+6gUhKCFekwQBNEqKFCpAVNUxKqfrGCkbWeFgngiitmcOESvAStNZqkfO3jqpwFFhQ0kdFKaLBKPBHDL+wYRDram+iNkY6Zdv10NVGQFmMgWDc3InMDSPl3RUFPUh2qpH6aomD0qqxb2IBIM4IiFvTW3b5X6Ef0qpKgQBNEqKFCpQVdUr65gDdCyfGhVew/WYnmynaIiSvgs9TNcI/XTn2y8Oy1XVBwaaUXiLlWYRghIleXJ2UIZL+zW50iNpfN1BCrNS/sA9oqKoiiComJ87qWDSTz1n2c5aqame5sqFRVJAkIuFSaCIIh6ocuiGrBgoCQrvDQzW/DHdE2jmdb6rYyIgcoE86hUV1QGkqxLrFx32a7uUfH3LmbV8O25nROGwGUsXXS9Xb19fnMqCOwClVxR5v+L2aMCAN2xsCOliu3bBYOiorfmp14nBEG0Cn+fRXxAQghGWIly1gc9VACzmdZOUVFv3zeZQ07riTG/xjTf7rjqD1EATGTqU1WYouI29dNqrMy0T2lpH8ZYOg+3NFtRYf6krClQSeXV55WkynlQbrAy07JAvVnN3giCIKygI04NQsEAP2gzn0rOBz1UAGeKCrt924E0AGBOV9S2lJkRkCT0NehTSXeYoiIGKk+/NmG4jzhR1SlNV1S0wDRr8qhM53QjbSOqh1ULfdb8LeKgvJkgCMIr/H0W8QnJqLELaLbAunO2W1ERG77ZKCraCWfjHtVzUas0mTGQVLdXby8Vt11p2wXvo6KpBYqi4GlNUWGTscem3Qcq402cnAwA8Yh1Z1q7Hipu4WbaYqWiQhU/BEG0EjriOICVKKcLJkWl3YGKE4+K1kF066iqqCwaSDja9kCDioqbOT/tJGDqTPvqgQzG0gVEQgGcfshcAPUpKnrqpzmKip76MfZRmbbpoeIWy9RPiQIVgiBaD1X9OCBpKlH2jUclWltROWqxWoraGw/jtIPn4JqzDnG07YEudXv1Vv50iqLCshgsUHlRq/ZZuaCHe3mYOuKGiWYrKjapn6mcdQ8Vt7BgpFBWzbnBgMTTQFSaTBBEK6FAxQH6BGX1QO0Xj0oiUltRufDohXjjoXPRHQvzNIcTBrtYG/0GFRW/Byraa8KMxsw4O78nxv0lY/UoKlqlULPmhthV/fDS5FhjAZJY7l4oyYhHgrZTmQmCIJoJHXEcYJ6grJcnt/flMzR8q6Lu9CUiroIUQPeojKcbVFR8nvrpiqmvC0vVjAkBxkCXGmQcaMCj0vzUj9lMqw0k9EhRAXRDLaV+CIJoB3TEcQBP/RSM5cntNtMmHJhp64V1p61n3o+iKEjlOqM8uSeuro9NOhZNsIPJ+hWVZk5OBvTUT0lWUCzrwSQ30zboUQkHAzy4ZT4VXp5MgQpBEC2EjjgOYE3f/OZR6XJgpq2Xwa76q36efjWPA9MyIkFgQZ+/s4vdJkWFBSoDyQjvRjuWLkBR3DW+G29RwzfAmP5J2cz5qQdx8rb4nfqoEATRSuiI4wCmqGTMfVTarag4TP3Uw0AXG0qncL+JE8qygp8+Og0AOO+oJLrj/t7F2PqYp0Ss1hlMRgGoSkLaYvifHblimasQzWr4FgkGuL9GnEPlVdUPYJwTBagTuAFSVAiCaC10xHEAV1R4HxWfmGnDQfQnwoiHg3wyrlfEwgFuhB1z4VN5ZFMO2w+UkIxKeMfqpKdragZMUeGpnzRTVMKIR4I8GHXTS4WpKaGA5ImyYYUkSejR3vOpnN7iP+VR1Q9QOZhQ96hQwzeCIFqHv3V5n8AVlYIx9dNuj0ogIOFX/3oy8lpVhtcMdAWRGSthLF3GooHau4qiKPjV31U15aLjkuiO+T8OZoHKVK6EYlmuMMEOJCPYNZHFgXQeSwad9aBh6kxfItLUmTi98TDG0gVMZvRAxSuPClDZnZbKkwmCaAcUqDiAVf3o5cnqlWW7Uz8AcMi87qZteyAZwM4x5z6VTEHBvin1NTrvSGcn9XaTjEqQJEBRVCWEKyqmQMWJofbRLfvxwMYRhDQPR7OMtAxdUbFI/UQbf26W+mH7O5UnEwTRDihQcUAyYmOmbXPqp9kMuOylMpllowUkxCOdcTILBCQko8B0Dtg3mefpPdb/hBlqq3Wn3bBzEl+67yU8snm/4fZmGWkZLN03mRVSP3kvPSos9WMsTyYzLUEQrYQCFQcw06q5hX67+6g0G9ZL5YDmUZEVBdsPlLB4IGTZlyWlBSo9PjfQmunSApWt+9W0VTAgoUc70VcrUd47mcMX/vAi/vDcHgBAOCjhza+bh79u2Y9UroSF/fGmrputUQxUprXpyZ6kfkxVP9yjMsP3e4Ig/AUFKg5ImrqA6g3fZraiopcoq//vuo1Z3Hb/FN5zUhfedUJXxf2ntEClt9MClRiASWDLiBqo9CfC3FvCFJVxi0Dlyp8/jSdfG4ckARcdvRDXnHUolgwmMJUr4sGXRnDKijlNXbeVosJSP402fAP0/ZuZadn3KCkqBEG0EApUHKB7VPzVR6XZsMGErDvtK/vUE+KOsZLl/Vnqx+8lyWaS2kDpLdrgRrGbLO9Omy7g0S37MZUt4dxV8/H8rkk8+do4wkEJd/3bqVi1sJc/picWxoVHL2z6ulmgMqUFKmqjveaVJ5NHhSCIdkCBigP0oYT+mvXTbPTBhOr/O6IZZdM5a3PtFE/9+LsbrRm18kfBllFVURkQAxXt5/Xbx3H3+l0oyQq+d+lq3PfCXgDAeasWGIKUVmIOVPIlGSVtuKIn5cm86sfYmZbKkwmCaCUUqDhA76NinPUzWxSViYw6QXeUBSp56y6tHetR0RSVbfuZoqJXzLDUD1NbAOD6OzfwdMsHT17aolVWYk79pITqn2TEQzMtVf0QBNFGKFBxQLeQ+pFlZdakfnoTAQQkQFbUYGUkpf7fqVqKSgf0TxFhgQpTDgaEiceDXfrP0VAA83pi2D6WAQC8bn43jlva37qFmugxBSrTQvv8gMshlFZUdKalPioEQbQBOuI4gHkWZAXYP52Hpq4jNsNTP8GAhL6Euou8ur8ETVBC2qal/mSHKyoMg0dFa6MPAP/2poPxlX88GiwG+ODJy5ra0K0WZkVl2sOutEBlH5U8lScTBNEGSFFxQCQUQHc0hFS+hF0TWX57bBbk6ge7ghhLy9i0R696SecVyIqCgOkkncqpEVynBSpJzaPCGEjqqZ+FfXGsmJtEMhrCv77xIMTCQdz0jiPw7M5JXHxs8w2z1eg1tdBPeViaDADRsHUfFSpPJgiilVCg4pCBroghUAkGJISDnWUarYd+zafy0m69BFZWgGxBQTJq/P87tTy526SoiI3aIqEA7r/2jSjLCu84+54TluA9J7RyhdbYKipeBSoho5mWFBWCINoBHXEcwk5eu7VAJR4OtlX2bxWs8mfzPmMfEav0z1SHliebUz/mjrKSJPEgxU8wj0quKCNfKns6kBAQ+qiQmZYgiDZCRxyHsA6lu8bVQGWmN3tjDGqKSsHUOmU6Z6z8yZcU5Iqdmfrpihp/7082t/W9V3RHQ2Cx8mS26OlAQsCijwqVJxME0QY664zSRlglyK6JHAB0zCybRmGKihmzosJKk0MBIBHpLKUpGoYhjdfsYYJeEQhIvCJtKlsyVP14gTn1Q4oKQRDtgI44DtEDFT31MxtgHhUGO52bFRUx7dNpKTFJkgzByUCHKCoA0JvQfSopDycnA+JQQuZRKWu302GDIIjWQUcch7CT1+5ZFqgMmhSV4X7192mTotKppckMlu4JSGoL/E5B7E7LBhJ6ZqYNm1I/pKgQBNEG6IjjEHYiYxUWs8WjMmBSVJbNUU+Mdqmfjg1UNGWiLxHxpFlaqxArf1jVT4/HHpWceXoyBSoEQbSQth5xHn74YbztbW/D8PAwJEnC3Xff3c7lVGXQlA6YLYFKIhpALMwmCQfQpwUuKZvUT6cHKp3iT2Ew9WcqV/S86kdP/bDOtKSoEATRetp6xEmn0zjqqKNw6623tnMZjjD7FmZL6gcABrrU3WRuTxBdUfVns6Iy2aHt8xn9SRaodI4/BRAUlUwR4xm1hNzzPipFGbKs8IGH1EeFIIhW0taGb+eddx7OO++8di7BMRWBygxvny8ykAxi93gZQz1BdMVUdcVspmXzfzpVUWGDCDulNJnBApWJbBEv71OnP6+Y2+XJtvXOtDIvTRZvJwiCaAUd1Zk2n88jn8/z36emplr23OZAZbakfgBgkCkq3UEkNUXFbKbt9NTPUYt7IEnAMUv62r0UV7Cmbxt2TWI6X0IkFMDBQx4FKkIfFZb2AUhRIQiitXRUoLJmzRp87nOfa8tzd0VDCAclFMuqkjCbUj9vOSKB6ZyCM14fx54J1QeRzs2sQOUNhwzimRvO4QpFp8DWu377OADgsHndCHsUSMR41Y/MfSoAZsXoCIIg/ENHnVWuv/56TE5O8q8dO3a07LklSTKoKrOl4RsAHLYggk+9vR/D/SF0xZiiYmem7dyTWKcFKYCuqLAAetXCHs+2zc20RdlQmtxpfXIIguhsOkpRiUajiEajte/YJPoTEeybUlNPs0lREUnamGk7XVHpVMzB1crhXs+2bZX6odJkgiBaDR11XDDYpSsqs8mjIsLMtOm8grJWBVKWFW6u7bTJyZ2OOVA5fNh7RUVWgEyeutISBNEe2qqoTE9P45VXXuG/b9u2Dc888wwGBgawZMmSNq7MmoGkrubM1kCFKSoAkCko6I5JmM7JYImgrg4tT+5UxEAlIAGvn+9hoBLW38upnNrokIy0BEG0mrYGKk8++STOOOMM/vu1114LALj00kvxwx/+sE2rsmdAaAY2W1M/4aCEWFhCrqggnZPRHQvwHipdMQnBDurqOhMQA5UVc7s8LZsXg5KUFqhQaTJBEK2mrYHKm970JiiKUvuOPkFUVGZTHxUzyagaqDBD7WRGDVT6KO3TcrqF5m5epn0AdTpzJBRAoSRjSut6S4oKQRCtho46LhhIkqIC6Omdaa1EmRlpexO0O7WacDCApBY0H+6hkZbBPClT2owrap9PEESroaOOC8ijotIVZYZaNUCZyLCKn9n7mrSTwS51v1y1sBmBivqesjlCFKgQBNFqOqo8ud30i4rKrE79GAcTMo9KHykqbeE/L1iJ53ZO4MTlA55vmykqLFChqh+CIFoNBSouGBQ9KrNZUYkZe6kwjwqlftrD2Svn4eyV85qybVb5w6t+KFAhCKLF0FHHBYbOtLM5UNFSP9PmQIXMtDMOPfVD5ckEQbQHOuq4oD8RRkgrv01GZ2+gkmSKiin1Q4rKzEM305JHhSCI9kCpHxeEggHc+PbDMZkpcAPjbKTLNEGZUj8zFzaYcF8qB0BXWAiCIFoFBSou+cBJS9u9hLbDTLNj01qgkqXUz0yFBSZbR9MAgNXL+tu5HIIgZiF0ZiFcM79XPXntnSwhV5SRK6opIKr6mXmIVT6DyQjecczCNq6GIIjZCJ1ZCNfM62UGSwV7J9RhdZEgEAtT+/yZhtgy//0nLZ3V/YMIgmgPFKgQrolHAjzNs2mvWg3SkwhAkihQmWkwRSUSCuADJ1PakyCI1kOBClEXLP2zaU8BANCXoCvtmchgl1qS/w/HLsKcWWwgJwiifZCZlqiLeb1BbNpbxKY9qqJCRtqZyb+cdhAW9MTwrtWL270UgiBmKRSoEHWhG2pVjwqVJs9M5nZHcdmpy9u9DIIgZjF0diHqYn6fMcbtIUWFIAiCaAJ0diHqgikqDCpNJgiCIJoBnV2IujAHKuRRIQiCIJoBnV2IuuiJBwx9U8ijQhAEQTQDOrsQdSFJkkFVoUCFIAiCaAZ0diHqZp4YqFDqhyAIgmgCdHYh6mZ+r1r5I4GqfgiCIIjmQGcXom7m96mKSndMQjBA7fMJgiAI76FAhaibJQOqojK3h9rnEwRBEM2BOtMSdXPYgjA++uYerBgKt3spBEEQxAyFAhWibiRJwlmHJ9q9DIIgCGIGQ6kfgiAIgiB8CwUqBEEQBEH4FgpUCIIgCILwLRSoEARBEAThWyhQIQiCIAjCt1CgQhAEQRCEb6FAhSAIgiAI30KBCkEQBEEQvoUCFYIgCIIgfAsFKgRBEARB+BYKVAiCIAiC8C0UqBAEQRAE4VsoUCEIgiAIwrd09PRkRVEAAFNTU21eycwgU8wgnUqjGCgiGoq2ezktIV/KoyAXMDU1hVK41O7lEH6hlAFSaSBQBIKz47PgCeU8IBeA2BQQos8TYQ87b7PzeDU6OlBJpVIAgMWLF7d5JQRBEARBuCWVSqG3t7fqfSTFSTjjU2RZxu7du9Hd3Q1JkjzZ5tTUFBYvXowdO3agp6fHk20S3kLvUWdA71NnQO9TZzDT3idFUZBKpTA8PIxAoLoLpaMVlUAggEWLFjVl2z09PTNiZ5jJ0HvUGdD71BnQ+9QZzKT3qZaSwiAzLUEQBEEQvoUCFYIgCIIgfAsFKiai0ShuvPFGRKPk9Pcr9B51BvQ+dQb0PnUGs/l96mgzLUEQBEEQMxtSVAiCIAiC8C0UqBAEQRAE4VsoUCEIgiAIwrdQoEIQBEEQhG+hQEXgG9/4BpYvX45YLIbjjjsOjzzySLuXNKv57Gc/C0mSDF/z58/nf1cUBZ/97GcxPDyMeDyON73pTXjhhRfauOLZwcMPP4y3ve1tGB4ehiRJuPvuuw1/d/K+5PN5XHXVVZgzZw6SySTe/va3Y+fOnS38L2Y2td6jyy67rOKzddJJJxnuQ+9R81mzZg2OP/54dHd3Y2hoCBdddBE2bdpkuA99nihQ4fzyl7/E1VdfjU9/+tNYv3493vCGN+C8887D9u3b2720Wc3hhx+OPXv28K8NGzbwv33pS1/CLbfcgltvvRVPPPEE5s+fj7PPPpvPgCKaQzqdxlFHHYVbb73V8u9O3perr74ad911F+644w785S9/wfT0NC644AKUy+VW/RszmlrvEQCce+65hs/WH//4R8Pf6T1qPg899BCuuOIKPPbYY1i7di1KpRLOOeccpNNpfh/6PAFQCEVRFOWEE05QPvKRjxhue93rXqf8x3/8R5tWRNx4443KUUcdZfk3WZaV+fPnKzfffDO/LZfLKb29vcq3vvWtFq2QAKDcdddd/Hcn78vExIQSDoeVO+64g99n165dSiAQUO69996WrX22YH6PFEVRLr30UuXCCy+0fQy9R+1hZGREAaA89NBDiqLQ54lBigqAQqGAp556Cuecc47h9nPOOQePPvpom1ZFAMDmzZsxPDyM5cuX4z3veQ+2bt0KANi2bRv27t1reM+i0Sje+MY30nvWRpy8L0899RSKxaLhPsPDw1i1ahW9dy1k3bp1GBoawqGHHorLL78cIyMj/G/0HrWHyclJAMDAwAAA+jwxKFABsH//fpTLZcybN89w+7x587B37942rYo48cQT8eMf/xj33Xcfbr/9duzduxennHIKDhw4wN8Xes/8hZP3Ze/evYhEIujv77e9D9FczjvvPPzsZz/DAw88gC9/+ct44okncOaZZyKfzwOg96gdKIqCa6+9FqeddhpWrVoFgD5PjI6enuw1kiQZflcUpeI2onWcd955/OcjjjgCJ598MlasWIEf/ehH3PhH75k/qed9ofeudfzjP/4j/3nVqlVYvXo1li5dij/84Q+4+OKLbR9H71HzuPLKK/Hcc8/hL3/5S8XfZvvniRQVAHPmzEEwGKyIPkdGRioiWaJ9JJNJHHHEEdi8eTOv/qH3zF84eV/mz5+PQqGA8fFx2/sQrWXBggVYunQpNm/eDIDeo1Zz1VVX4Xe/+x0efPBBLFq0iN9OnycVClQARCIRHHfccVi7dq3h9rVr1+KUU05p06oIM/l8Hhs3bsSCBQuwfPlyzJ8/3/CeFQoFPPTQQ/SetREn78txxx2HcDhsuM+ePXvw/PPP03vXJg4cOIAdO3ZgwYIFAOg9ahWKouDKK6/EnXfeiQceeADLly83/J0+Txpts/H6jDvuuEMJh8PK9773PeXFF19Urr76aiWZTCqvvvpqu5c2a/nEJz6hrFu3Ttm6davy2GOPKRdccIHS3d3N35Obb75Z6e3tVe68805lw4YNynvf+15lwYIFytTUVJtXPrNJpVLK+vXrlfXr1ysAlFtuuUVZv3698tprrymK4ux9+chHPqIsWrRIuf/++5Wnn35aOfPMM5WjjjpKKZVK7fq3ZhTV3qNUKqV84hOfUB599FFl27ZtyoMPPqicfPLJysKFC+k9ajEf/ehHld7eXmXdunXKnj17+Fcmk+H3oc+TolCgInDbbbcpS5cuVSKRiHLsscfyEjGiPfzjP/6jsmDBAiUcDivDw8PKxRdfrLzwwgv877IsKzfeeKMyf/58JRqNKqeffrqyYcOGNq54dvDggw8qACq+Lr30UkVRnL0v2WxWufLKK5WBgQElHo8rF1xwgbJ9+/Y2/Dczk2rvUSaTUc455xxl7ty5SjgcVpYsWaJceumlFa8/vUfNx+o9AqD84Ac/4Pehz5OiSIqiKK1WcQiCIAiCIJxAHhWCIAiCIHwLBSoEQRAEQfgWClQIgiAIgvAtFKgQBEEQBOFbKFAhCIIgCMK3UKBCEARBEIRvoUCFIAiCIAjfQoEKQRAEQRC+hQIVgiAIgiB8CwUqBEEQBEH4FgpUCIIgCILwLRSoEARBEAThW/5/XCT4xia+6dcAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(resids1_box, dist_box)\n", "plt.ylabel('Ca distance (Angstrom)')\n", "plt.axvspan(122, 159, zorder=0, alpha=0.2, color='orange', label='LID')\n", "plt.axvspan(30, 59, zorder=0, alpha=0.2, color='green', label='NMP')\n", "plt.legend()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## References\n", "\n", "[1] Oliver Beckstein, Elizabeth J. Denning, Juan R. Perilla, and Thomas B. Woolf.\n", "Zipping and Unzipping of Adenylate Kinase: Atomistic Insights into the Ensemble of OpenClosed Transitions.\n", "Journal of Molecular Biology, 394(1):160–176, November 2009.\n", "00107.\n", "URL: https://linkinghub.elsevier.com/retrieve/pii/S0022283609011164, doi:10.1016/j.jmb.2009.09.009.\n", "\n", "[2] Richard J. Gowers, Max Linke, Jonathan Barnoud, Tyler J. E. Reddy, Manuel N. Melo, Sean L. Seyler, Jan Domański, David L. Dotson, Sébastien Buchoux, Ian M. Kenney, and Oliver Beckstein.\n", "MDAnalysis: A Python Package for the Rapid Analysis of Molecular Dynamics Simulations.\n", "Proceedings of the 15th Python in Science Conference, pages 98–105, 2016.\n", "00152.\n", "URL: https://conference.scipy.org/proceedings/scipy2016/oliver_beckstein.html, doi:10.25080/Majora-629e541a-00e.\n", "\n", "[3] Naveen Michaud-Agrawal, Elizabeth J. Denning, Thomas B. Woolf, and Oliver Beckstein.\n", "MDAnalysis: A toolkit for the analysis of molecular dynamics simulations.\n", "Journal of Computational Chemistry, 32(10):2319–2327, July 2011.\n", "00778.\n", "URL: http://doi.wiley.com/10.1002/jcc.21787, doi:10.1002/jcc.21787." ] } ], "metadata": { "kernelspec": { "display_name": "guide", "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.9.15" }, "toc": { "base_numbering": 1, "nav_menu": {}, "number_sections": false, "sideBar": true, "skip_h1_title": false, "title_cell": "Table of Contents", "title_sidebar": "Contents", "toc_cell": false, "toc_position": {}, "toc_section_display": true, "toc_window_display": false }, "vscode": { "interpreter": { "hash": "5edc5d8d8cbc0935a054a8e44024f729bc376180aae27775d15f2ff38c68f892" } } }, "nbformat": 4, "nbformat_minor": 2 }