{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Initializing session" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "import matplotlib\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "from classy_sz import Class\n", "from scipy.optimize import fsolve\n", "from scipy.interpolate import interp1d\n", "import math\n", "\n", "font = {'family':'STIXGeneral'}\n", "axislabelfontsize='large'\n", "matplotlib.rc('font', **font)\n", "\n", "plt.rcParams.update({\n", " \"text.usetex\": True,\n", " \"font.family\": \"sans-serif\",\n", " \"font.sans-serif\": [\"Helvetica\"]})\n", "\n", "\n", "cosmology = {\n", " # LambdaCDM parameters\n", " 'H0':67.556,\n", " 'omega_b':0.022032,\n", " 'omega_cdm':0.12038,\n", " # 'sigma8': 0.831,\n", " 'ln10^{10}A_s': 3.047,\n", " 'n_s': 0.9665,\n", " 'tau_reio':0.0925,\n", "}\n", "\n", "\n", "import os \n", "path_to_class_sz = os.environ['PATH_TO_CLASS_SZ_DATA']" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Galaxy kernels" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0.9999999999999999\n", "1.0000000000000002\n" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhQAAAGcCAYAAACBcTnNAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAABJR0lEQVR4nO3deXxcZ2Hv/8/MaJcsjeR9S+JRNrJbskNIgACWgaRsASmGFnqBYKlACy1tLdzLLVsvxoZS2ltK7QChLZSfLbEUwiolQICQxJaAbCQkGifxEsebRpK1jJY5vz+eOTOSrGX2c0b6vl8vvWY7Z+bR8bHOd57VY1mWhYiIiEgavE4XQERERPKfAoWIiIikTYFCRERE0qZAISIiImlToBAREZG0KVCIiIhI2hQoREREJG0KFCIiIpK2glx9UCQS4fjx4yxZsgSPx5OrjxUREZE0WJbFwMAAa9asweudvR4iZ4Hi+PHjrF+/PlcfJyIiIhl05MgR1q1bN+vrOQsUS5YsiRWosrIyVx8rIiIiaejv72f9+vWx6/hschYo7GaOyspKBQoREZE8M193BXXKFBERkbQpUIiIiEjaFChEREQkbTnrQyEislBNTEwwNjbmdDFEUlJYWIjP50v7fRQoRERSZFkWJ06cIBQKOV0UkbT4/X5WrVqV1jxRChQiIimyw8SKFSsoKyvTpH2SdyzLYmhoiJMnTwKwevXqlN9LgUJEJAUTExOxMLF06VKniyOSstLSUgBOnjzJihUrUm7+UKdMEZEU2H0mysrKHC6JSPrs8zidvkAKFCIiaVAzhywEmTiPFShEREQkbQoUIiKLjMfj0cgUyTgFChEREUmbAkWG9Y+M8cunThOJWE4XRURyzLIshkbHHfmxrNT+5oRCIZqamqiurqa6uprW1tYpr3s8Hrq7u6mvr8fj8VBfXz+ldmPfvn3U1tbGXuvu7k74vefaV/KPho1m2D/c/TgHDh3lnTdexMfecKXTxRGRHBoem+CKv/+xI5/9+CdeQ1lR8n/St2zZQiAQoKurKxYAAHbv3h3bprW1lba2NgC2bt1Ka2sre/fuJRgM0tLSQltbG3V1dezevZumpiZ6enrmfe/59pX8o0CRQZZl8dMnTwHw1fufYeMFft543VqHSyUiMrPOzk66u7u555578Pv9AOzdu5etW7dOCRQtLS0EAoHY/Y6ODgCCwSAADQ0N+P1+du/ezaFDhxJ677n2lfykQJFBwdODnBoIxx5/5DuP8opLV1BVVuhgqUQkV0oLfTz+idc49tnJsi/qGzZsiD03U2fNurq62H07HIAJA4FAgOrqahobG9m2bRuNjY0Jvfdc+0p+Uh+KDHogeAaA6y+q4bKVSxgYGefLvww6XCoRyRWPx0NZUYEjP6nMI1BTU0NdXR29vb2xH8uyzuuPUVNTM+t79PT00NHRQU1NDdu3b6e2tjbh955tX8lPChQZ9GDwLAAvqV3KXzZcAsBXfvUMoaFRAJ7vG+aP73yAD3zjNwyPTjhWThERgEAgcF5HyPb29lhfh/l0dnayb98+Ghoa2Lt3L729vQSDQTo7O+d977n2lfykQJEhlmXFaihuCCzlNVeu4vJVSzgXHueO/zjETx47QeMXf839PWf47u+O857/PMjImEKFiDinrq6Ouro6tm7dGruYb9++nc2bNyf8Hi0tLXR2dhIKhdi3bx9ggkoi7z3bvpKf1IciQ46cHebkQJgin5eNF/jxej3sabyGt3/pQbqe7aX5v7oAuHBpGacHwvzq6TN84u7H+dRtVztcchFZzO655x62b99OfX09AM3NzezYsSOhfRsaGtixYwdNTU2EQiECgQBtbW2xUDDXe8+3r+Qfj5Xq4OUk9ff3U1VVRV9fH5WVlbn4yJx6MHiGbfseYMOycn76N6+IPf/EiX7u+OohJiIWr7lyJR/Ycgm/f36At3/5QQD+493Xc/Olyx0qtYikamRkhMOHD7NhwwZKSkqcLo5IWuY6nxO9fquGIkPODpp+EjXlRVOev3xVJb9sfeWUDlMvvaSYd954EV+9/xn+97cf4Rc7XqkFhkREJK+pD0WGnJklUMDMq7jteO1llBR6Odo7TM+pc1kvn4iISDYpUGSIXUOxdIZAMZOyogI2rq8G4MHDZ7NWLhERkVxQoMiQM+fMhFYz1VDM5voNZmz3QwoUIiKS5xQoMsRu8lhaUZzwPi+OBooHg2dTXthHRETEDRQoMiTZJg+AjRdUU+jzcKJ/hKO9w9kqmoiISNYpUGTIbKM85lJa5OOadX4gPm23iIhIPlKgyJC5RnnM5Zp1VQA8fVIjPUREUhUKhdi6davTxUiYx+MhFAoRCoWyMm1AZ2dnbG2Uzs5OWltbM/4Z0ylQZEAkYtEb60ORXKBYV10GwNGQmjxERFK1a9cuWlpanC5G0uyl27OpoaEhNsV5NilQZED/yBjjEdOpMtkairV+MyPZcQUKEcmRffv2UVtbi8fjob6+fsoiXsFgkK1bt+LxeKitrWXPnj2x12b6Nr1nz54pi4l5PB6CwSC1tbWx9Tmmv2d7e/uU92xqaqK6uprq6uop36TnKudkoVCI9vb2Kcufezweuru7qa+vj+0/+YI61+850+/h8Xhi3/qrq6tpaWmZ8v7Ta0cmv2ZvPxu73J2dnXg8nvN+7JqM2Y7T9M/bunXreceqpaUl+7UUVo709fVZgNXX15erj8yZp08OWBe23m1d9fc/Snrfh4+ErAtb77Y2/0NHFkomItkyPDxsPf7449bw8LDTRUlKT0+PBVhtbW1WT0+P1dzcbAUCgdjrgUDAamxstHp6eqyOjg7L7/dbe/futSzLsnp7e63pl43du3dbjY2NsceA1dDQYHV0xP+m+f1+q7m52erp6bH27t1rAVZXV5dlWZZVV1cX+7yuri4rEAhYO3bsmLec08vQ3Nw85Tm7HD09PVZPT48VCASmbDPX7znT72E/7u3ttdra2izA8vv9sXIDU/af/DvP9v69vb0zHlNbc3OzVVdXN+dxmvzvMv0YTz9ec13y5zqfE71+K1BkwEOHz1gXtt5tvXzPvUnve3pgxLqw9W7rwta7rZGx8SyUTkSyYc5AET43+8/ocBLbDiW2bRI6OjpiFzPLMhcj+6JpX/gm2717d+zClGigaGtriz1ua2s77z2bm5uttra288oyuXxzlXO6urq6KRfrmcqxe/duq6GhIaHfc6b97TLZ/H5/7IJuWZbV0NAw5fH08kx/fb5AYR83+/ee7TjZnzU9PMwUwPx+/6zHMBOBQmt5ZMCZc6l1yLT3KSn0MjIW4UTfCBcuLc908UQk1z61ZvbXLnk1/Elb/PFnLoaxoZm3vfCl8K7vxx9//moYmmFE2Mf6Ei5aQ0MDgUCA6upqGhsb2bZtW6zKvbu7m02bNp23fbJV5XV1dbH7Bw8epKGhYcrre/fuBYg1iWzYsCH2mt0sMVc5pwsGgzOuUjq5HH6/P3Y/0d9z8v7AlH1qamqmLMU++f3BrKy6b98+enp66O7uprOz87z3m00oFGL79u3ceeed+P1+gsEgMPNxAujq6jrvGNfX19PZ2Xle+e33ygb1ociAVOagsHk8Htb4SwE4prkoRCQHenp66OjooKamhu3bt8dGA6TizJnzA05NTU1C+9bU1FBXV0dvb2/sxzI150mV017+PNVyzFW+VIRCIWpra2lra6O2tpbdu3fPGoZmsmXLFm6//fbYPvMdp0TL7vf76enpSf4XSpBqKDLg7KCZdntpeeKzZE621l9K8NSgRnqILBR/d3z21zy+qY//9uk5tp32ne8vH0m9TFGdnZ0Eg0Gam5tpaGhg7969sQ6HdXV17Nq1a8r2hw4dOu9iHQqFYt/I5/vGu3nz5imdMAGamprYunUrmzZtOq/zYHt7O/v376elpWXWck7/Nm5/i58pVMwk0d8zVZ2dnZw9ezali3drayuhUChWiwMQCARmPU5tbW3U19efN1Kko6PjvPe2g062qIYiA07bTR5JDhm1ras2NRQa6SGyQBSVz/5TWJLEtqWJbZuklpaW2DBCu9khEAjQ0NBATU0NTU1NBIPB2PwFdlOAHSJ27doVG1kxvVp9usbGRs6ePUtrayvBYJB9+/bR3t7Opk2bqKuro66ujq1bt8Y+b/v27bGmhNnKOV0gEEiqKn++3zNdfr+fUChEd3d3rOzt7e3zDtvs7u5mz549tLW1TXl+vuN0++23EwwGYyGsvb2dAwcOnPf+mQxNM1ENRQak0+QBsKZKTR4ikhsNDQ3s2LGDpqamWFNBW1tb7ELT0dFBS0sLtbW1BAIBdu7cSXNzc2z/trY2tm/fzr59+2hoaKCtrW3W4Zy2rq6uKe/Z1tYW609wzz33sH37durr6wHT92DHjh0Ac5Zzsm3bttHV1ZXUcZjv90xHQ0MDzc3NbNmyBTC/k33ctm7dOmvzh10rYR+LyWWd6zj5/X66urpizUJ1dXW0tbWdF0xCodB5tTuZ5LHmaoTJoP7+fqqqqujr66OysjIXH5kz77rrIX765Cn2vOUabt+8Pun9v9V9lA8d+B03XbyUr7/nhiyUUEQybWRkhMOHD7NhwwZKSkrm30GyJhQKUV9fn9X+Afmuvb2djo6OKU0pk811Pid6/VaTRwYMjU4AUF6cWoWPOmWKiKTO7/fT2Nh4Xl8Nidu1a1fWZ+RUoMgAO1CUFfnm2XJma/12H4oRLWMuIpKCnTt3zvrte7GzO7JOH9qaaepDkQGDo+NA6oFi+RIzOmR0IkL/yDhVpYUZK5uIyGLg9/tnHNkgpk9HNvtO2FRDkQHDsRqK1PJZSaGPimhzyZlz4YyVS0REJFcUKDJgMBytoShOrYYC4quU2sugi4iI5BMFigyIdcpMsYYC4kNOVUMhIiL5SIEiTaPjkdjS5aUp9qEAWFph+lHYk2SJSH5QR2pZCDJxHitQpGko2iETUu+UCbDMbvJQoBDJC4WFpvP00NAsC3uJ5BH7PLbP61RolEea7OaOogIvhb7U85m9DsiZQTV5iOQDn8+H3+/n5MmTAJSVleHxeBwulUhyLMtiaGiIkydP4vf78flS/2KsQJGmoTSHjNrUKVMk/6xatQogFipE8pXf74+dz6lSoEjTYDj9DpkQ70OhTpki+cPj8bB69WpWrFjB2NiY08URSUlhYWFaNRM2BYo0pTtLpm1ZufpQiOQrn8+XkT/IIvlMnTLTlLkmD7sPhQKFiIjkHwWKNA2mOUumrSZaQ9E7NMr4RCTtcomIiOSSAkWahqM1FOVpzJIJUF1WiMcDlgW9Q2qLFRGR/JJyoGhtbSUUCmWwKPnJ7pRZmmYNRYHPS3WZPdJDHTNFRCS/pBQouru72bNnT6bLkpeGx+xRHul3yFqqjpkiIpKnUgoUwWCQQCCQ6bLkpdjCYGnWUEB8LorTGjoqIiJ5JulA0d7eTmNjYzbKkpcyNWwU4iM9zmqkh4iI5JmkvlaHQiH8fn9C24bDYcLh+Dft/v7+pAqWL2LDRtPslAlq8hARkfyVVA3FgQMHaGhoSGjbXbt2UVVVFftZv359SgV0u9iw0cL0A0VVqVmUpW9YozxERCS/JBwoOjs7uf322xN+4507d9LX1xf7OXLkSEoFdLshuw9Fcfp9KBQoREQkXyV1FTxw4EDsfjAYZNeuXWzbto26urrzti0uLqa4uDj9Erqc3Yci3bU8QIFCRETyV8JXwelNHS0tLbS0tCz60R6Z7JSpQCEiIvkq6VEeoVAoNgfF7t276e7uznih8kmm1vKAeKDoV6AQEZE8k3Q9vd/vZ8eOHezYsSMb5ck7sSaPTPShKDOBIqRAISIieUZreaTJntiqNMNNHpZlpf1+IiIiuaJAkab41NuZ65Q5EbFiw1FFRETygQJFGkbHI4xNmJqETExsVVroo8hn/knUMVNERPKJAkUa7A6ZkJmJrTweD5V2s4eWMBcRkTyiQJEGu0NmUYGXAl9mDmVVqWk6UQ2FiIjkEwWKNNg1FJlYutymuShERCQfKVCkYXg0AkBJBpo7bPFAoQXCREQkfyhQpGFk3DR5ZCdQqIZCRETyhwJFGsJjpoaiuCBzh9FfZpYwV6AQEZF8okCRhnC0hqI4gzUUlaqhEBGRPKRAkYaRaA1FSQZrKOJNHuPzbCkiIuIeChRpyEYNhfpQiIhIPlKgSEN2aygUKEREJH8oUKQhqzUUQxo2KiIi+UOBIg0jWRnloRoKERHJPwoUaQjH5qHIfJNH/8i4ljAXEZG8oUCRhngNRQaHjZbElzAf0hLmIiKSJxQo0pCNGoqSQi8+rweAgRENHRURkfygQJGG8Hjmayg8Hg9LSsyKo+fC6kchIiL5QYEiDSNj0VEeGeyUCVBRbAJFv2ooREQkTyhQpMGuocjk4mAAS6L9KNTkISIi+UKBIg3hLNVQxJo8FChERCRPKFCkIWs1FNEmj4ER9aEQEZH8oECRhmz1oYh3ylQNhYiI5AcFijRkuw+FOmWKiEi+UKBIQzgLU28DVJSoyUNERPKLAkUaRmKLg6lTpoiILG4KFGkIZ2HqbZjcKVOBQkRE8oMCRRpGsjD1NsT7UKhTpoiI5AsFijRkrYZCfShERCTPKFCkyLKsrPWhqFCTh4iI5BkFihSNTVhYlrmf+RqK6NTbavIQEZE8oUCRInvpcshGHwo1eYiISH5RoEjRSLT/hMcDRb7sBIqRsQhjE5GMvreIiEg2KFCkyK6hKC7w4vF4Mvredh8K0FwUIiKSHxQoUjSSpREeAAU+L6XR6bzVMVNERPKBAkWKJtdQZEOsH0VY/ShERMT9FChSZNdQZHphMFu8Y6ZqKERExP0UKFKU7RqKCnvoqAKFiIjkAQWKFGVr6XJbpb1AmJo8REQkDyhQpCg8lqM+FKqhEBGRPKBAkaJs11Bo+m0REcknChQpGslyDUVFsfpQiIhI/lCgSJFdQ5HphcFsFdEmj0Gt5yEiInlAgSJFdg1FSRYmtgKoKDbvq0AhIiL5QIEiReGx7NZQlNt9KBQoREQkDyhQpCjW5JG1Ggo1eYiISP5QoEhRrFNmtvpQKFCIiEgeUaBIUbZrKNTkISIi+USBIkWxTpmqoRAREVGgSFXu+lBMZOX9RUREMkmBIkXZrqGwmzwGR8eJRKysfIaIiEimKFCkyK6hKPJldy0Py4KhMdVSiIiIuylQpGjUDhRZmnq7uMCLz+sB1I9CRETcT4EiRaMT2e1D4fF4KC8y731OgUJERFxOgSJFo7FOmdk7hHbHzHNaIExERFxOgSJF2W7yAC0QJiIi+UOBIkV2k0c2A4U90kNNHiIi4nYKFCkKR0deZGuUB0xq8lCgEBERl1OgSFEuaig0W6aIiOQLBYoUhXPQKTPe5KF5KERExN0UKFKUk06ZsUAxlrXPEBERyYSCZDbu7OwEIBQKcfDgQbZt20ZdXV1WCuZmlmXluMlDNRQiIuJuSV0Nm5qaqKmpobGxkdraWpqamrJVLlcbj1hY0eU1in3ZmdgKNMpDRETyR1KBoq2tbUqNhN/vz3R58oLdfwKyXUNhwoo6ZYqIiNsldTVsaGiI3W9ra6OlpSXjBcoHo7kKFCWqoRARkfyQVB8KgO7ubvbv38/WrVtpbm6edbtwOEw4HI497u/vT62ELmQHigKvJ7aAVzaUFylQiIhIfkj663VdXR07d+6kp6eH9vb2WbfbtWsXVVVVsZ/169enVVA3ycUID9A8FCIikj9SuiL6/X6amppoamoiFArNuM3OnTvp6+uL/Rw5ciSdcrrK6ER0lsxsB4oSLQ4mIiL5IeErYmdnJ9XV1bHHgUAAgGAwOOP2xcXFVFZWTvlZKOxOmdmcdhs0ykNERPJHwlfEmpqaKZ0yu7u78fv9i3Ieipw3eYxOYNnjVEVERFwo4U6ZdXV1bNu2jX379gHQ0dFBV1dX1grmZuEcBQq7hmIiYjEyFqG0KHtzXoiIiKQjqVEejY2NsftzjfBY6EZj63hk9wJfVujD4wHLMs0eChQiIuJWWssjBblq8vB6PbGhoxrpISIibqZAkQJ7HY/iLHfKBCiPzpapjpkiIuJmChQpyFUNBUxecVSBQkRE3EuBIgVOBAo1eYiIiJspUKQgPB6d2ConTR6qoRAREfdToEhBroaNgpo8REQkPyhQpCDWKVNNHiIiIoACRUpy2Yci3uQxkfXPEhERSZUCRQpy2ilTC4SJiEgeUKBIgUZ5iIiITKVAkQK7U2ZOJraKTrd9blSBQkRE3EuBIgWO9KFQk4eIiLiYAkUK4qM8sr9Y15ISNXmIiIj7KVCkwJlRHgoUIiLiXgoUKcjlxFYKFCIikg8UKFJgN3nkYurtJRrlISIieUCBIgXhsehaHjmsoRjUxFYiIuJiChQpiNVQ5DBQjE5EYouSiYiIuI0CRQrsTpm5XMsDVEshIiLupUCRglyO8vB5PZQWmuGp6kchIiJupUCRglyuNgoa6SEiIu6nQJGCWA2FL/sTW0F8cisFChERcSsFihTkch4KgPLi6HoeChQiIuJSChQpyGUfCoDyIs1FISIi7qZAkYJcjvKA+EgPLRAmIiJupUCRJMuycjoPBUCF+lCIiIjLKVAkyQ4TkMs+FJotU0RE3E2BIkl2cwfkZi0PmNTkER7LyeeJiIgkS4EiSWFHA4VqKERExJ0UKJJk11AU+jx4vZ6cfGa5VhwVERGXU6BIUnyER24mtQKo0DwUIiLicgoUScr1CA+AiuJCQIFCRETcS4EiSfFpt3N36OyZMtXkISIibqVAkaTwuOkYWVyYyxoK9aEQERF3U6BIUngst7Nkgia2EhER91OgSJI9bLSkMHedMu21PBQoRETErRQokjQyFm3yyGmnzILoZ0cYnzRTp4iIiFsoUCQp7MCwUXseCtD02yIi4k4KFEmyO2WW5LBTZlGBNzZM9dyomj1ERMR9FCiSNDKW+xoK0EgPERFxNwWKJMWGjeawDwXE56IYGFGgEBER91GgSFJs2GgOmzwgPlumaihERMSNFCiSNBKroch1k4dmyxQREfdSoEiSUzUU9kiPAQUKERFxIQWKJDkxbBTUKVNERNxNgSJJ9sRWuRw2CgoUIiLibgoUSXKqhsJu8jinia1ERMSFFCiS5NSw0YpYoBjL6eeKiIgkQoEiSU4sDgaTmzxUQyEiIu6jQJEkJxYHg8lNHupDISIi7qNAkaR4HwpnZso8p5kyRUTEhRQokmTPQ5HrJo8lJdEmDy0OJiIiLqRAkaQRp9byKFKTh4iIuJcCRZLiM2U6NGxUTR4iIuJCChRJcqoPRazJQzUUIiLiQgoUSQrHZsp0poZicHSCSMTK6WeLiIjMR4EiSU7VUNjzUIA6ZoqIiPsoUCQhErEYnXAmUBQXeCn0eQAYUD8KERFxGQWKJNi1E5D7Jg+Px8OSkkJAgSKrzvTAZy6G/3g9nHrS6dKIiOQNBYok2Ot4QO5rKCDeMVPreWTQkz+E7v8CK9ovpcQPg6fg8H1w6C5HiyYikk8K5t9EbHYNhc/rocDnXKDoVw1Fak49CWeDUFgKBaUwEoL/eT8MnQFvAVz3Niithlf8HfzsU2ZbERFJSFKBoru7m87OTgAOHjzInXfeid/vz0a5XMlex6PEgdoJiHfMVJNHCgZPw96Xw/jI+a+tuhqueou57/XCBS829xUoREQSltSVsbOzkx07drBjxw42b97Mli1bslUuV4qN8Mhx/wlbvA+FmjySNjEK9e8y91dcAdUXwYorTZDY9jUoKIpvWxMwt73PQESru4qIJCLhQNHd3c2uXbtijxsbG+nu7iYYXDzf4mKzZDpUQ2E3eaiGIgWVa+CWT8PH+uB9v4YP/g7edz80fsWEiynbrgVfEUTGoO+oI8UVEck3CV8Z6+rquPPOO2OPQ6EQADU1NRkvlFvZnTJzPcLDVhmtodD021nm9cVDhpo9REQSklQfisbGxtj9/fv309DQMGsfinA4TDgcjj3u7+9PrYQuMuJwDUW8D4WaPJJiWfDs/bD6Gihektg+r/o/4PHAqmuyWzYRkQUipVEeoVCI9vZ2urq6Zt1m165dfPzjH0+5YG4UdmilUZuaPFIUeha+eiv4imHn0an9JWZzxRuyXy4RkQUkpStja2srHR0dc47w2LlzJ319fbGfI0eOpFpG14hPu+1sp0wNG03SsW5zu+JFiYUJERFJWtI1FHv27KG1tZVAIBDrRzFTsCguLqa4uDjd8rmKPWy0uNDpGgo1eSTleDRQrK1LfJ+Rfgj+DIZ7of5/ZaVYIiILSVJXxvb2durq6mJh4sCBA4tqHgrnayjsmTJVQ5GU/uPmduklie8zdBoOvAN+uAMikfm3FxFZ5BKuoQgGgzQ1NU15zu/309zcnPFCuVXYNTUUChRJCQ+Y25LKxPepWGVux0dgbAiKKzJfLhGRBSThQBEIBLDs9Q4WqZFoDUWJw30o1OSRJDtQJDrCA6CgJH5/PKxAISIyDy0OloTYxFYuqKFY7OEuKakECq/XTG4FM0/XLSIiU2hxsCQ4P2zU1FCMRyzC4xHHJtjKO5vebWa8tKfUTlRBiZmyW4FCRGReChRJsCe2cupCXlbow+Mx8zT1j4wpUCRq8x2p7VdQDGEUKEREEqAmjyQ4XUPh9Xq04mgu2f0oFChEROalGookOD1sFMx6HgMj4woUiYpMwIlHTP+JmoCZTjtRt34GJsagekP2yiciskAoUCTBntiqxKFOmaDJrZI2dBb23Wzu/31vcoHisluyUyYRkQVITR5JcEMNRWxyK9VQJGY0OsKjaIkZuSEiIlmhGookxKbedqgPBaA+FMmKDRlNYR6JZ+83s2yu2xRfzlxERGakr2xJGB41gaK82MkaCnuBMDV5JCSVOShs930WvnkHPPvrzJZJRGQBUqBIwlA0UJQWOVexo+m3k5ROoNAoDxGRhClQJGFo1FzEy4ocHOVRqhqKpITPmduUAkV0tdzxcObKIyKyQClQJCFWQ+HghFLVZSZQ9A6OOlaGvBLuN7dFKfShUA2FiEjC1CkzCfE+FM4dtppy86357JBqKBKy6hp46V/BssuS31c1FCIiCVOgSJBlWQy6oMmjplw1FElZv9n8pEI1FCIiCVOTR4LC4xEi0QU+Sx0MFNVlZgXMswoU2acaChGRhKmGIkF2cweYRbqcstRu8lCgSEzoCETGoGIlFJUnt+8Vb4Tll8OKy7NTNhGRBUQ1FAkaik5qVVTgpcDn3GGrjjZ5DI9NTAk5Mosf74R/2Qi/+0by+66tg+veBms2Zr5cIiILjAJFgobCzvefADNTZqHPrEfRO6RainnF5qGodLYcIiILnAJFguwho042dwB4PB71o0hGOhNb9R2DJ38Izz2Y2TKJiCxAChQJis+S6WygAKgpV6BIWDqB4vDP4Rtvhfs+k9kyiYgsQAoUCRoeM00eTs5BYbMDhZo8EpDW1Nv2KA8NGxURmY8CRYIGw87PkmmrVg1F4uxAkdZMmRo2KiIyHwWKBNkjKpzulAmwVIEiMZEJGLXX8kihU6ZqKEREEuZ8/X2eiC8M5vwhU6fMBEUmzLTb4QEoSSVQqIZCRCRRzl8d88Sgi2oo1IciQQVF0PCxNPbX1NsiIolSk0eC3NTkoT4UOaImDxGRhKmGIkHxYaPOHzL1oUjQC4+BtxCWXgzeFLJz5Vq49bOaFEtEJAGqoUiQPWzUFTUUsT4UWsJ8Tvf+A3xhMzy0L7X9y2rg+u1w7bbMlktEZAFSoEiQPWzUDYFich8Ky7IcLo1LWRYcPWTuay0OEZGsU6BIUGzqbRc0edgLhE1ELPqGVUsxo/5jMHgSPD5YfU1q7xGZgMO/gKc6zX0REZmVAkWC3NTkUVzgw19mQsXJAQ1pPM/QWXjuAXN/5ZVQWJra+0TG4T9eB19/S3w+CxERmZHzX7fzxJCLRnkArFhSTGhojJP9YS5dmcK00gvZF2+CgePm/tr61N/HVwR4AEtzUYiIzEM1FAkaCrunyQNgxRIzR8IL/RrSOIVlwUgo/viyW1N/L49Hc1GIiCTIHVfHPDAUbfJww2qjACsqzRwJavIADt0FE6NwdZMZmbHzGAz3mtfKl6b33gXFMD4MYwoUIiJzUaBIkJsmtgLVUMREJszy4v3HoGwpXN1o5pxIN0jYVEMhIpIQNXkkyO5DUe6SJo+V0RqKU4u9hqLnpyZMlPjh8tdl/v1js2Uu8uMsIjIPd1wdXS4SsSbNlKkaCsd1/D2EnjMh4rFvmeeu2QaFJZn/LNVQiIgkRIEiASPj8TkIXNPksZj7UDx9L7zwSPzxyqvMqqLZcNMHzWqlNRuy8/4iIguEAkUC7NoJgNJCdwSKldEaipMDI1iWhcfjcbhEOfSq/w1nD8O5E7DyarjqzeDN0r/Lxj/JzvuKiCwwChQJsIeMlhb68HrdceG2ayhGxiL0j4xTVVrocImyLBKBcy9A5Wq47BanSyMiItOoU2YChlw0S6atpNBHZYnJg6cGFkH7/vFu+Nzl8NXXmbkmcuVMDzz7a+h/PnefKSKShxQoEmAvDOaWDpm2FZV2x8xF0I/iDz8yt+XLzIRTudL5UbjrtfDk93P3mSIieUiBIgEDI2YBrsoSdzUrrFhid8xcBDUUx7rN7Yabc/u5sVEeiyC0iYikQYEiAf0jpsmjstRdXU5WLqYaCnvmy4qVuf3c2DwUiyC0iYikQYEiAf3D7q6hWBRzUdjrc5T6c/u5qqEQEUmIAkUC+u0mD5eNpFhVtYgmtxoOmdsSf24/1w4Uo4O5/VwRkTyjQJGA/mHT5OG2oZmrok0ez/ct8EBhWTDSZ+7nuoaicq257TuS288VEckz7uoU4FJ9Lm3yiNVQLPRAMTEG177N9KMorc7tZ9cEzO3ZYG4/V0QkzyhQJCDe5OGuw7W6qhSAFwbCTEQsfC6ZdCvjCorgTV9w5rNXXQ03t8Lyy535fBGRPOGuK6RLubVT5vIlxfi8HiYiFqfPhWOjPiSDqtbCK//O6VKIiLie+lAkID5s1F2Bwuf1xEZ6LOh+FOOjplNkLmfIzIYnvg9dX4XnHnS6JCIiGadAkYCBaA2F2zplQrwfxYm+YYdLkkVP/Rg+tQbuutWZz+8/Dk/fA6f+kN773P+v8L0PwldeA0cPZaZsIiIuoUCRgFinTJf1oQBYXbUIRnrYQ0aLK5z5/Ps+A197Mzy8P733ueAGWHElYMEPW82CZyIiC4QCxTwsy4p3ynRZHwqIz5Z5YiHPRWFPapXrOShsqY70CD0HX341PN1pmmsaPgrv+BYUVcCxQ/BIW+bLKiLiEAWKeYyMRRibMG33butDAfEaihOLoYaipMqZz081UDz6TTjyIPzy8/EFzZasgpf8ubn/+P9krIgiIk5zXx2+y9i1Ez6vh3KXrTYKsCo6dHRBN3k4Ne22zQ4Up540fSCKl0DglXDxlrn3633G3F5409Tn12w0t5osS0QWENVQzCM+ZLQATy6XzU7Q4qqh8Dvz+dUbTDPF+LAZpfFwW2K1JaHnzK1//dTn7cf9xzJaTBERJ6mGYh7xDpnua+6A+PTbJ/pGsCzLlaEnbU7XUBSWwDvvhucegPAArK2HdZvm3y8UrYGomhYoll4Cf9Edn9ZbRGQBUKCYh5s7ZILplOnxwOhEhDODoyyrKHa6SJm3/gbwFUNNrXNlWLMx3lSRCMuKN2n4L5j6WkERLHXwdxERyQIFinnYC4O5ccgoQFGBl+UVxZwcCHM8NLwwA8XNf+t0Cc53pgd+8Y8mONz2xfNfHzwF4yOARzURIrIoqA/FPNxeQwGwttp0zDzau4Ant3KbyAT89uvw+Hdmnk+i/zh4vLBktamRmO7hA/Cd98EffpL1ooqI5IICxTz6XTxLpm2t3wSKYwsxUFiWO6fdrgmYZpixIeg9fP7ra66Dj5yE5p/OvP9zvzaB5OjBrBZTRCRXkgoU3d3d1NfXZ6ssruTWdTwmW1ddBsCx0AIMFKODZtrtTy43993CVwDLLzP3Tz4+yzaFZt6JmVStM7d9RzNfNhERByQcKNrb2wETKhaTvqH4sFG3WtBNHiN90TsWFJY5WpTzrLzK3L4wS6CYiz3yQ3NRiMgCkfBVsrGxMZvlcK1YHwo311DYTR4LsYZi8rTbbhsSu/IKc3vysfNf6/goDDwPL24xw0ynUw2FiCww6kMxj9BQPnXKHHK4JFlwLFojNn3opRusiAaKmWoonvqJWUxsuHfmfe1A0X9Mi4SJyIKQtXr8cDhMOByOPe7v78/WR2XV2cFRAJZWzNBT3yXsTpkDI+P0j4y5Ovwk7ckfmNtLX+tsOWay8kpzu26TCQXeaD4fOgtnox01/RfOvO+S1WYUyMSoqcmo0tBSEclvWauh2LVrF1VVVbGf9evXz7+TC50ZNKFoabl753coLy7AX2ZCxIIa6TE6BD3RURKX3+psWWZSsRL+1/fgtn+PhwmAez5upuleccXsk3H5Cs3rJX71oxCRBSFrNRQ7d+7kQx/6UOxxf39/3oWKiYgVq6FY5uIaCoB11aWEhsY41jvMi1ZXOl2c9Njf9g//3FyYqy6Id4B0E48HNrw8/nj/2+H003DqCfP4j/5xatCY7ra9Zl0Pp1ZRFRHJoKwFiuLiYoqL3futPhGhoVEi0ekPqsvdHSjW+kt59Fh/fnfMtCz40U4zEdTWT8Dyy+HmD0Nxhfs6ZM7kTBBO/d7cv+5P4MIb595+1aSQ9HAb/KgVCsth2cVw2a2w6Y65A4mIiIukFChCoRB+vz/DRXGfM9HaCX9ZIYU+d/9hX+tfAHNR/PKf4MEvwsZ3mHBRswFeudPpUiXuTf8WHeZqwQUvSW7fcB8MnQHOQN9z0HOvGSa78U+yUVIRkYxLOFB0dnbS0dEBmP4RmzdvXvBDSU+fs/tPuLt2AhbASI9TT5q+BwCrrs6PGonp1lyX+r5XNcKFN8FIPzxyAA5+CTr+Hi7/I+dWWRURSULCgaKhoYGGhgZ2796dzfK4yplz9ggP9zfd5P3028/eb24vepmZu2GxKfXHg8OajXD4Pjj9B/jyq+F9vwavz7x24tGpTSUiIi7h7np8h52J1lC4vUMmmE6ZkMdNHiceNrczTQK12BQUwes+b/pTlC+PhwnLgm+8DY7/1snSiYjMyL3zSbtAbA4KFw8Ztdk1FKfPjTIyNkFJoc/hEiXp+WigWHW1s+Vwi4tugg89bpZBtz20zwwx/cOP0mteERHJAtVQzOF0HkxqZfOXFVJeZEJE3tVSRCbghej01auvdbYsblLqh2WXxB8XlAAWPH2PUyUSEZmVAsUc7CaPfOhD4fF4Yh0z864fxdBZ0y+gfPnsE0EJXLzF3B47NPuU3iIiDlGgmIPdKXNZHozygHizR96tOlqxHN7TCX/zlOZdmEvVOlh2GVgRCP7c6dKIiEyhv95zODOYP6M8ID509FgoT4eO5uNQ0Vyzayl61OwhIu6iQDGH2DwUedCHAiZNbpVPNRSDp+HQXU6XIn+s22xuTz3pbDlERKbRKI9ZhMcnGBgZB2BZHozygMk1FC4OFCN9cO//hd7D4PGZuRbO9pg+AS/70Pz7L3bly8ztcMjRYoiITKdAMQt7yGiB10NlaX4cpnVu75QZicA3t8NTP576fGkNvOj1zpQp36zbDB/6PZRWO10SEZEp8uNK6YD4LJlFePKkbX9dtFPmif4RxiYi7lt/5JefM2GioMQs/lUQrfm5eCtUrXW2bPmisNT8iIi4jALFLE4OjAD5MamVbVlFMcUFXsLjEY72DrNhWbnTRZqqoMSEhyveCHXvcLo0IiKSQS77CuseL/SbDpmrqkocLknivF5PLEQET51zuDQzuPHP4e3tChPp6vgofPvP4Nyp+bcVEckRBYpZvNBvaihWVuZPoAC4eEUFAD1uDBSSGb/7hvkZeN7pkoiIxChQzMKuoVhZmT9NHgC1y02gePpkjgJFZCKx7fqP6xt1ptgdMjVbpoi4iALFLPK1hqI2VkMxmJsP/OJN8A8rof0OExpmc99n4bMXw89256ZcC5kChYi4kALFLOKBIr9qKC6eVENhWVbmP+DIQ2b4p210EMZH4NF2+LcboPeZmfd74VFzu1RrdaRNgUJEXEiBYhbxJo/8qqEILC/H44G+4bHY1OEZMz4KX2+C738IJsykX/zZfXBHJ6y4wkxa9et/O3+/SCS+mujKKzNbpsVIgUJEXEiBYgZjExHODOZnoCgp9MUmuOrJdD+Kwz+HkRA8cXf8udJqWL8ZXvMp8/g3Xzt/FsfQMzB6DnzFsPQSJE0KFCLiQgoUMzg1EMayzCyZNWX5sY7HZHbHzIz3o3j8O+b2ijeCb9oUJoFXwPIXwdgg3P8v8ectC366y9xfddX5+0nySv3mVoFCRFxEgWIGdv+JFUuK8XrzY5bMyexA8dTJgcy96cQYPPF9c/+KN53/uscDL3mfuV+1Lv78z3bBIwfMuh1bPpq58ixmm+6ADz0Bt37W6ZKIiMTo6+IMYh0y82hSq8kuX7UEgMeO92fmDcdG4Oe7zTfi8hVw4Y0zb3fd22F0CDZOmrjKHlZ66x4I3JyZ8ix2ZTVOl0BE5DwKFDOIdchckp+B4pp1fgAeO9ZHJGKlX8vy7ZZ4c8cNfwZe38zbeb3m9SmFuR2uvM00d4iIyIKlJo8Z5OuQUVvt8nJKC30Mjk4QPJ1kPwrLgl98DkLPxZ8bH4HiKrhtL7zsr5N7v+WXKUxk2sAL8KO/gx9+2OmSiIjEKFDMwK6hWJFnIzxsBT4vV66pBOCRY6Hkdn7yh3DPx+FLDTASbTJpvAt29MC1b81sQSU14yPwwBeg66tOl0REJEaBYgb5OkvmZFevqwLg4aN9ye348H5ze+WbocSEEorKwFeYwdJJWuxho+PDMDbsbFlERKIUKGbwfJ/5I706TztlAlwTDRSPJBMoxkbg6c7oGzRloVSSEcVLzKgZOH/ODxERhyhQTGNZFs/3mRqKfA4UV6/1A2akx/hEZO6NbcGfmQmoKtfCmrqslU3S5PFocisRcR0Fimn6h8cZGjVDHVdXlTpcmtQFlpVTUVzA8NgET5xIcD6K33/P3F7+R+aiJe5VttTc9h1xthwiIlEKFNMcjzZ3VJcVUlo0y/DIPOD1eqi70HyLPfjM2fl3eOZXZoEvgMtfl8WSSUZc8GJz23Ovs+UQEYlSoJgm3n8if2snbC/eYCZAmjFQHP4F/L9NZgZMgOIKKKqAS18LF70sh6WUlFzyanN75mlnyyEiEqWJraZZCP0nbJsvMoHiocO9WJaFx27GOP0U7P8TMyw0Mm5GcKy+Ft7TAUtWmwmqxN1qt8AHfgM1AadLIiICqIbiPM+HooHCn/+B4pp1VRT5vJw+F+awPcGVZUHbO81S46uuAs+kU6AmAIX5XzOzKBSVxcPE878zP0MJNG2JiGSJAsU0xxdQk0dJoY9r15vho7FmjzM98MKj4CuCt38LCvJzNlCZ5H/eD3tfDns2mPuRBEf1iIhkkALFNHYNxZoFUEMBcH20H8UDwWigsDvxXXADVKxwqFSSURWroHy5uf+br8E9H3O0OCKyOClQTHOi3+5Dkf81FAA31S4D4BdPnSYSsaDnHvNC7ascLJVk1Nvb4W+fNmutAPzqn+HoIWfLJCKLjgLFJJZlcTxkmjzWLJBAUX9RNaWFPk6fC/PEsTNmdAcoUCxE174VXvR6c7/np86WRUQWHQWKSXqHxgiPm/bnlVULo29BcYGPl9SaSZAeeuIZuOo2WHUNrLza2YJJdlz0cnP77K+cLYeILDoKFJPYtRPLKoooLsjfSa2me9klptnjJ8+Mwxu/AH/2Cw0NXagCN8PGd0DdO5wuiYgsMpqHYpKjvSZQrPUvgOYOe6Gv3md4/fAYId8jfPGZJgbD45QX6599wVp+GbzxX50uhYgsQrqyTPLMGTNXw0XLyh0uSZp6fhqdayIEwDLgrwrh30dez31/OMUtV692snQiIrIAKVBM8mw0UFy4NI8DRegItL/LhInKtWZ4qGXx0OlCip4d50ePnVCgWOgmxuHE7+DUk3Dt27TQm4jkhALFJM+cHgLgoqVlDpckQZYFx38Dyy41a3EAfKvZLGm9ZiO8+8exiat8z/Yy8MX7uff3JwmPTyyoPiIyzZmn4EtbwZqA/uNQ/y4oqTRTrIuIZIl65k2SV00eE2OmJuLOV8LZnvjzta80S1s3fXXKLJgb1/tZsaSYgfA49/ecyX15JXdWvAhe+2lz/95PwmcC8Mll8JmL4Qc7IPScs+UTkQVJgSJqZGwitjDYRW5v8rAs+OZ74LFvR5+YVKV9/XZ4/0GovmjKLl6vh9dcuQqA7z/8fG7KKc65fju8/G+huDL+3OApeGgvHLrLuXKJyIKlJo+o586a5o4lJQVUl7m8avj538Lj3wFvIbztG7D6mvhrpdWz7vaG69bwXw88yw8eeZ6PveFKKjTaY+HyeOBVHzE/E+MQ7odj3eb2qjc7XToRWYBUQxFlr8a5YVl5fJlvt/rN18ztFW+ES7YmvNumC6sJLC9naHSCu393PEuFE9fxFUBZDVzSMDVMnPoDPPcAnHwCep81zWgiIinSV9SovBnhMTYMj7SZ+xvfntSuHo+HbZvWs+uHT/D/HTzCW6+/IAsFlLwwNgIH3gGnnog/V1AC618Mt+w2/TBERJKgGoqow9ERHhvcPsKj+z9hpA+q1sOGm5Pe/c116yjwevjtkRCPHO3LQgElLww8bzrvVq2H0hrwFcP4CBz+OXypAZ7qiG9rWc6VU0TyhgJF1OHT5wAX1lAc64JH2uOP119v/vhvendK02cvX1LM664x81D828+ezlQpJd/UbIB3/QD+6lFoPQwfeQHe/xBseDmMnjM/tke/Cf90NTz+P86VV0RcT4ECmIhYsW/rV6ypnGfrHLEsuO+zcOcWuOcT8efXbIR3/whu+mDKb/3eV1wMwI8eO8HTJwfSLaksBB6Pmbb77d82Q04n136FnoW+5+DAn8K9/9e5MoqIqylQAH94YYDB0Qkqigu4dOUSp4tj/PRTZg4BLDNx1eRq57V14E19YqrLVi1h6xUrsSz4fOdT6ZdVFg5fAdzwXtOJ07bp3XDjB8z9X3xW81iIyIwUKIDu53oBuHZ9FT6vC0Z4PPeA+cMNcMseeHt7xqdP/suGS/B44O6Hn+eBoCa6kjmUVsOrPwmBV4AVgQf3Ol0iEXEhBQqg+9kQAHUXzD6HQ05Ylul02fZO84f72rfBi1uy8lFXrqnibdFRHh/77mOMT0Sy8jmygNzwPnPb/V8QVlOZiEylQAH8JlpDkfNAceIROHoo/nhsCL77F6YHfk3ADN/Lor959WVUlRbyxIkB/vshVWPLPC7eCksvgXAffOd9Gv0hIlMs+kDROzhKMDqp1XXr/bn74GfvN8PzvrXdzC0BMDpk/mC/6v9Ay31QUpXVItSUF/HXr74UgH/8yR84Ozia1c+TPOf1mpDrKzbNH3Yz3NP3mJEgChgii9qiDxS/6jkNQGB5OdXlRbn50Gfvh2+81Yz7918YDxQVy+EvDsHL/waKc9M59I+vv4DLVy2hb3iMT//w9zn5TMljF2+BD3TD5jvizz1xN7S/2/zY57KILDqLPlB8s+soALdetTrLH/QeuOuPYO/NcNctZnKq9TfAW/97ao/6HCvwefnkm64C4MCho/zyqdOOlUXyRNW6qY+XXQreAnjsW/Crf3GmTCLiuEUdKE72j/DzP5wC4C3162bf0LLgDz82Sz8HfxZ/fmwE+o7NvP305488BM/+0izs5fFC3Z/CH++HIudn5tx8UQ3vuOFCAD78rYfpG9aaDpKEG94Lr48Gie7/gMiEs+URWWx+tx/+X/3UOYscsKjX8vj2b44RscyiWRuWzTJD5sQ4/HAHHPqyeXzRTfHXgj81TRc1tbD8ciiuMM0YLzwG507C+x+EyjVm21t2m+pgbwGsusp0unSR1lsu594nTnK0d5g//+9u7nrnZgp8izpvSjKuegv85H9D/zF4uhMufY3TJRJZPM4G4czTMOTsFACL9ooxMDLGV351GJijdsKy4Nst0TDhgfp3wupr46+ffNw8f7YHnvw+PLzfTE985mmIjEPPT+PbXnaLWenxije4LkwAVBQXsPcd9ZQW+vjFU6f58LceYSKiTnaSoMISuPaPzf2DX1IHTZFc6jtibqvWO1qMRVtDsedHT/JCf5iLlpZx28a1M2/0y3+CR9tNrULTV+FFr5/6+sv+GurfZZoxzvSYGoiCEqhYYTqv5ahjZaZctbaKz7/1Ot77tS7au44SHo/wmcZrKClMfVZOWUTq3wkP/Bs89RN44vvwoteZGr7hs2Zp9JEQnHvB1N6NDpq5VlZdDRfcYPYfPA33/4tpSvQVQokfqi+E2ldB+TIHfzERl7Nnr1WgyKHnHoTBU9x73MePHzxDERV86rarz79gRibgvs/Azz5tHt+y5/wwYSurMX/wal+V3bLnyGuuXMW//nEdH/jGb/je747zxPP9fO7267h6XXaHsMoCsPxSeM2noPcwXHarea73MPzrptn3ecmfxwPF6CD86p/P38ZbANe+FV73eRM0RGQqu4bCn0eBIhgM0t7eTiAQIBgM0tzcjN/vz1LRMuDcKXPBj657MfrwNyk6tJdXAQ8VR7f570JTk1BcAe/4DiytNTUNj7QBlvmDN3mI3CJw69Wr8ZcW8sH9v+Wpk+d4wxd+yZs3ruOdN16kYCFze0l0wqvpU8V7C8y8KhUrTQ1eUYX5f7nyqvg2pdVww/tN88nEmKnZeP5hOPEw/OZrUHUBvKI1t7+PiNtFIvFBAA7XUHgsK/HGzvr6erq6ugATLlpbW2lra0to3/7+fqqqqujr66OyMosrevYdg9/+Nxz+OTz3ayJv3c/j5Zv53sPHGXnoa9w28UNWec6ywtOHl2nTTf/lo/GE98JjZibLa9+avbK63JlzYT5x9+P8z2+Px567ck0ljfXruLF2GZesqMDrhrVPZGF7/Ltm9MhbvmRCx3gYhkNmtFRk3EwD3veceS4yDms3wTKzoi6hI/D775pax8j4pNtxE2hqXzW1huSFx6GkEoorzW1hWcbX0RHJqP7j8LkXgccHHzlpFvjL9EckeP1O+JODweCUx4FAgM7OztRLmEHj4xMMBR/Ae/DfKXv6B3it8dhrd33j63xyOBx9dCP3LdtK62sv57VXLDd/iEbPmdvwOfPtybbySvOziC2tKOaf37qRd954EXf96hl+9OgJHjvez2PHHwegqrSQjRf4uXh5BbUrKriwpozq8iJqyovwlxVSXKC+F5IBV7zBNDnaF/bDv4Cvv2X27f/oc/FAcTYIP/672bctKo8HilNPwJcbpr5e4jcdqgOvNP2i7L4co4MwcMIElIlRiIyZWpWJMfN42aVQFe2bdfIJ02Hb6zM1NR6fuV9SCWXLYPU18Y7aYyOmn0lBCZT6oaCYOUUipi9KLPR4zH2FoMVjdAhWX2fOgyyEiWQk/OmdnZ3U1EydgKmmpobu7m7q6uoyXrBEvekz3+VjAx/lOm888DwQeRHfnbiR+yNX8Iy1mrIiHzddvIzG+nW86vIVFNrDIUv95kfmtPGCajZeUE3v4Cjf+s0x7n3iBbqfDdE3PMbPnjzFz548NeN+ZUU+yosLKCn0UlLgo6TQZ+4X+uI/BV4KfB7A/AH0eMw9c+uJPTavTd5m+mvmdQ/mrWZ8ffJz9rYzvGa/F3O9Pqm8dtnir016btrnTd93+nOOcejjr1hdyVVrE2hGm3KBtDAFtkwtRdESc/EuW2r6WNhDtQGWrIKrm8yF3L6g2z+RcVizMb5tJAL+C2CkH8L95g/0SAh+9w3z86YvwnXRkSyHfwHf2DZ7eV/3edj0LnO//xj88nOzb7v1k3BTdHn4538LX5k05LawzNSWeKJ/s172Ibh+u7n/3IPwlVfP/r43fxheudPcP/kE/PtN5rgVlJgm3uIl0cDiMb/XDe812w6cgP1vj36mHVAm3b/8dXDDn5lth0PwzTtm3s7jhQ0vjy9wOB6G77wXvIVQUAS+omggG4OJMKzbHN92Ygz+843mdnJYs+9vuBne9IX47/rFl4I1YX63ghLTbFZQam5XXmk60Nvu+YQpi11GK2L+7SPj5t//xj+Pb/uDv4Whs+a9I9Ef+77/AnjdpH/XA39qasQi4+Y9rUj8nKtaB9u+Nul9d5jzoqjc/Bvb/w4ej2kafMWH49s+dKfpdOnxxLexz/+CErNty89nPw9yKOFAEQqFZnz+7NmzMz4fDocJh8Oxx/39/cmVLEFDvkqKmCBsFfADXsrdpW9gbMXVXLS0jHcuK+fqdX6uWlupb8sZUF1exB0v3cAdL93A2ESEx4738+ixPnpOnSN4apCjvUOEhsboHRolYsHQ6ARDo5rkSGb3wS2XJBYoJrtkK3wsFB+aOte38eWXmaaSRKzfDH/5iLlvWaYW4sTDZsTKkQfNH3+bx2su9F6fuUD6isy3Q/t+yaRq4aW18OL3Ri9c4/GfkT4zb8DS2vi2kQlzIRwfASyzYODYUPz10cH4/cKSxH4v8wuZzwQYHYPRAbMIoW3glfH7Y8Nw9ODsb7X88vj98bCZd2Q2Jf74/Ykxs+bLbCLj8UDh8cGzv5p928GTUx+ffNxc6GcydGZqoDj4JXPsZ7Lu+qmB4vffm3qcJlsxrQb7xKNmCoGZjEy7/j17P7zwyMzbLlk9NVA8fACOPjTztsXTwofD0q4fmS1o7Nq1i49//OPpvv28vvruF1N29it4l63ntsqV3Jb1TxSAQp+X69b7Z1xQLRKxGBgZp3dolKHRCYbHJgiPTTAyPsHIWITh0fj9kbGJ2HwXlgUWVvSW2AXDmuE1+1piYZ4wz03ad9o+RF+f7/2wn5u+76R9Jn/e1NeiZYi/VfwzZnk/Jj2XSdmYBSKJ7lYJm3VCuURks1rf4zHf4i+80fxMd+mrYeeRxN6r+iK45dOJbXvRTfCRE+Ybc7gPhnunXoyWTFoiYMUV8DdPmyATP8kxJ54FhaXxbZdeDB/6vXl+ImyaeUf6zX0LMzzXVrHCLAlgWSYE2e9n3588j05Jpam5iX1uZOq2Sy+Ob+srgtd+Ol4jMTFmvr37oiFs6SXxbb1eM1TfWxh9vXDq/dJJK0NbFrzrhzA+bJqMJt+Oh6ceMzDhbnw4XlaPJ94UNX1a+ZtboxMS+kyInNxsNbkMYGorxkai20S3t2s1pjddvXKnadoaHTQ/E6PxYzh9uoGrm2D99ef/+3q8SYbK7Eu4U+a+ffvYu3dvrFMmQHV1NW1tbTQ0NJy3/Uw1FOvXr89+p0wRERHJmEQ7ZSY8U+ZMoQFg06aZx5gXFxdTWVk55UdEREQWpoQDRSAwdbroYDDIpk2b3D0PhYiIiOREUn0o2traaG1tZfPmzRw8eDDhOShERERkYUtqYqt05GxiKxEREcmYjPehEBEREZmNAoWIiIikTYFCRERE0qZAISIiImlToBAREZG0KVCIiIhI2hQoREREJG0KFCIiIpI2BQoRERFJW9rLlyfKnpCzv79/ni1FRETELezr9nwTa+csUAwMDACwfv36XH2kiIiIZMjAwABVVVWzvp6ztTwikQjHjx9nyZIleDyejLxnf38/69ev58iRI1ofJAE6XonTsUqOjlfidKwSp2OVnGwdL8uyGBgYYM2aNXi9s/eUyFkNhdfrZd26dVl578rKSp1sSdDxSpyOVXJ0vBKnY5U4HavkZON4zVUzYVOnTBEREUmbAoWIiIikLa8DRXFxMR/96EcpLi52uih5QccrcTpWydHxSpyOVeJ0rJLj9PHKWadMERERWbjyuoZCRERE3EGBQkRERNKmQCEiIiJpy9k8FOkIBoO0t7cTCAQIBoM0Nzfj9/vT3nYhSub37+7uBqCuro5gMEgoFKKuri6HpXVWd3c327dvp6ura87tFvs5ZUv0eC328wrMMejs7ATg4MGD3HnnnfqbNYdkjtdiP7/s4xQKhTh48CDbtm2b9ffP+bll5YG6urrY/Z6eHquxsTEj2y5Eyfz+zc3NFmABVkNDg9Xb25uDErpDW1ub1dXVZSXyX2Cxn1OWldzxWsznlW337t1T7k8+h6bT+ZXc8Vrs55ff77e6urosy7KsvXv3WoFAYNZtc31uuT5Q9PT0nHdy+f3+tLddiJL9/ffu3Wv19vYuuv+Qk813gVzs59R0iQSKxX5edXV1TTlHenp6LMDq6ek5b1udX8kdL8vS+dXR0RG7v3fv3lnDlxPnluv7UHR2dlJTUzPluZqamli1V6rbLkSp/P5+v3/RVa8mY7GfU6lazOdVXV0dd955Z+xxKBQCOO88Ap1fkNzxsi3m86uhoSF2v62tjZaWlhm3c+Lccn0fCvvkmu7s2bNpbbsQJfv7h0Ih2tvbAdNu2dLSQiAQyFbx8tJiP6dSofMKGhsbY/f3799PQ0PDjBdAnV9GoscLdH6B6Ueyf/9+tm7dSnNz84zbOHFuuT5QzGa2g5XutgvRbL//5A46gUCArVu30tPTk7uC5bHFfk7NRedVnH3xm68j60z7LUaJHC+dX6ZWJxAI0NraSnt7+5RANp9snluub/Lw+/3nJaqzZ8/OmF6T2XYhSvb3DwaDsft2L+DJz4nOqVTovIprbW2lo6Nj1vNF59dU8x0v0Pll8/v9NDU10dTUNGNIcOLccn2gmNxeNNmmTZvS2nYhSub37+7uZsuWLec9P1e75WK02M+pZOm8ituzZw+tra0EAgFCodCMf/R1fsUlcrwW+/nV2dlJdXV17LHd1DNToHLi3HJ9oJjeNhYMBtm0aVMsZXV3d8cO5nzbLnTJHqvdu3fHtu3s7KSxsXHRHKvJpv/h0jk1t/mOl84raG9vj1VLh0IhDhw4oL9Zc0jmeC3m86umpmZKUOju7sbv98fmoXD63MqLxcGCwSB79+5l8+bNHDx4kJ07d8YOSlNTE5s3b2bHjh3zbrsYJHOs7Mlk/H4/PT09U/6jLnSdnZ10dHSwZ88eduzYwebNm2PtkDqnzpfM8VrM5xWY86W2tnbKc36/n97eXkDn13TJHq/Ffn61t7fHmjI6OjrYvXt3LDw4fW7lRaAQERERd3N9k4eIiIi4nwKFiIiIpE2BQkRERNKmQCEiIiJpU6AQERGRtClQiIiISNoUKERERCRtChQiIiKSNgUKERERSZsChYiIiKRNgUJERETSpkAhIiIiafv/AdFo1zK1jPRbAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "z , nz_lenses = np.loadtxt(path_to_class_sz+'/class_sz/class-sz/class_sz_auxiliary_files/includes/nz_lens_bin1.txt',unpack=True)\n", "plt.plot(z,nz_lenses,label='lenses')\n", "print(np.trapz(nz_lenses,z))\n", "z , nz_sources = np.loadtxt(path_to_class_sz+'/class_sz/class-sz/class_sz_auxiliary_files/includes/nz_source_normalized_bin4.txt',unpack=True)\n", "plt.plot(z,nz_sources,label='sources (normalized)',ls='--')\n", "\n", "print(np.trapz(nz_sources,z))\n", "plt.legend()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Compute" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/Users/boris/Work/CLASS-SZ/SO-SZ/mcfit/mcfit/mcfit.py:130: UserWarning: use backend='jax' if desired\n", " warnings.warn(\"use backend='jax' if desired\")\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "CPU times: user 25.5 s, sys: 2.17 s, total: 27.7 s\n", "Wall time: 3.05 s\n" ] } ], "source": [ "%%time\n", "\n", "M = Class()\n", "M.set(cosmology)\n", "M.set({\n", "'output':'gallens_lens_1h,gallens_lens_2h',\n", "'ell_max' : 5e3,\n", "'ell_min' : 2,\n", "'dlogell':0.1,\n", " \n", "\n", "'ndim_redshifts':80,\n", " \n", "'mass_function':'T08M200c',\n", "\n", "\n", "'M_min':1e11, # Msun/h\n", "'M_max':1e15, # Msun/h\n", "# 'mass_epsrel': 1e-3,\n", "# 'mass_epsabs': 1e-40,\n", " \n", "'z_min':1e-5, \n", "'z_max': 2.,\n", "# 'redshift_epsrel': 0.5e-3,\n", "# 'redshift_epsabs': 1e-40,\n", "\n", "'delta_for_galaxies':'200c',\n", "'delta_for_matter_density':'200c',\n", " \n", "'concentration_parameter':'fixed',\n", " \n", "\n", "'Delta_z_source':0.00,\n", " \n", "'galaxy_sample' : 'custom',\n", "'full_path_to_source_dndz_gal' : path_to_class_sz+'/class_sz/class-sz/class_sz_auxiliary_files/includes/nz_source_normalized_bin4.txt', # source galaxies\n", "\n", "'N_samp_fftw':1024, #precision parameter for the bessel transform to theta space\n", "'l_min_samp_fftw' : 1e-8,\n", "'l_max_samp_fftw' : 1e8,\n", " \n", "'hm_consistency' : 1,\n", "\n", "})\n", "M.compute_class_szfast()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "ell = np.asarray(M.cl_kg_k()['ell'])\n", "cl1h = M.cl_kg_k()['1h']\n", "cl2h = M.cl_kg_k()['2h']" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Plot" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAk4AAAGGCAYAAACNCg6xAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAACt0ElEQVR4nOx9d3wcxdn/d0+nU9ed5IIxzT6Zl1Bjn2xaaMaSIUDKiyWbhFACWAJCgPwCEoZATIBgmV5jyZDQwZapyUux5ARCAiFYokOCpcMG0wzWFZ10ujq/P9a73jttmd3bvSLP9/O5j6Rnn3nmmZmvZudmnpnhCCEEDAwMDAwMDAwMmrDl2gEGBgYGBgYGhkIBGzgxMDAwMDAwMFCCDZwYGBgYGBgYGCjBBk4MDAwMDAwMDJRgAycGBgYGBgYGBkqwgRMDAwMDAwMDAyXYwImBgYGBgYGBgRJs4MTAwMDAwMDAQAl7rh1g2IlkMokvvvgCVVVV4Dgu1+4wMDAwMDDsEiCEYHh4GNOnT4fNpj6nxAZOeYQvvvgCe+21V67dYGBgYGBg2CXx2WefYc8991TVYQMnC9Db24u5c+fC5XLpSldVVQWAb7jq6moLPOMxb948vPnmmzmzoycdja6WjpHnwWAQe+21l+VtYRbMalOr7TPO5A+s5oxZeRQKZ7R0GGeyl4cRG0JbCO9hNbCBkwXo6OhAb2/vOLnH40FfX59iOmF5rqqqyrJ/IkIIOI7LeDnQqB096Wh0tXSMPheucLSyLcyCWW1qtX3GmfyB1ZwxK49C4YyWDuNM9vLIhDMAqNKwgZMFcLlc6OnpSZF1d3ejtbU1Rx4xMDAwMDAwmAE2cLIAS5YsQUNDg/i31+tFfX09PB5PDr1iYGBgYGBgyBQTcuDk9XrR3NyM1atXqw5Wurq64Pf74XK5MDg4iMbGxpQBj1E0NTWl/N3Z2YmOjo6M7TIwMDAwMDDkFhNq4NTa2oqhoSG43W709/dr6tbX16OtrU2UNTc3A4ApgycBXV1dWLJkiWn2zIDWVkur7ehJR6OrpZPp80KA1WVgnDHuT74iG2UwI49C4YyWDuNM9vKw2s8JNXDq7OwEAPj9fqxcuVJRr7+/H729vaK+gI6ODtTX18Pn85nqU0tLi640Vp7hxHEcLrzwwozzMGpHTzoaXS0do89LS0txzTXXoLS0VNPPXMOsNrXaPuNM/sBqzpiVR6FwRkuHcSZ7eWTCGWpdIoSSTyD4/X7U1NSgr69PdqlOCNJOHzgBQE1NDbq7u8VZp3Xr1mHNmjWaeXZ0dMDtdqfIurq60NfXJ5uPHILBIJxOJ/x+P5xOJ1UavSCEIBAIwOl0ZrxDyogdPelodLV0jD43q56yAat9ZZyhe844Y34ehcIZLR3GmezlYdRGIBCAy+VCIBDQ3OE4oWacaNHb2zsuDkmA2+1GT0+POHBqampS1NVCd3c3CwhnYGBgYGCYQCj8RVcD8Hq9qKurk31WW1urGR9Fi97eXkyaNMkUWwwMDAwMDAy5xy45cNLC0NBQxjb8fj8A6D49nIGBgYGBgSF/scst1QkDGiW4XC54vd6M83G5XHC73ePinmgQDAZT/uY4DoQQlJSUoKSkZJw8HWpyAOL6rfSk1HR9GttSO3p8qa6uptaXrjUr6aeXR8lPJRtCPFn6MyW5kXq3Sq6nTY3KzbJPayefOKPEbSV/lOT5xBlBrlUHWnas7CMEeboNPT7S8IC2j6CtLzU7ejkjZ9+IPJ84oyUHMn8/admIRCKIRCLj5OnvXTXscgMnGpg1SzQ4OGgo3d577y0rb29vxxVXXAEAKCkpQVlZGUZGRhCPx0WdsrIylJSUYHh4GMlkUpRXVFSguLgYgUAAiUQCNpsNHMcfS2+z2RAIBFLycjqdSCaTGB4eFmUcx8HpdCIejyMUCiGZTMJms6GoqAjV1dWIRqMIh8Oivt1uR2VlJcbGxkSiEkJEeTgcRjQaFfXTy0QIQTKZREVFBUpLS2XLZLfb4ff7wXGc+E8jLZNgw2azweVyjSuToJ9IJDA6OppS1srKSsRiMYyNjamWCQAcDgfKy8s1y0TbTsFgMKVjUGunRCKBQCAgtqm0nUZGRkRdm81G3U7SMo2OjmJsbEy0n23uEUJACEFNTY1mmYT2djgcqKqqki2T1HeBM9IyxWIxkTPl5eXjykQIEcsqbSdCiMjJ9E5Y6//JjHbSwz2HwyH6LtSBEe5Z0UdIy2SUewIPqqqqUsoqVyatPkIYMAi+C/Wlh3tK/R4hRCyrlX2EGdyrqKgQfRfqwOx+z4z3UywWw/DwsGgjvUwrVqzI/FxFMgHh8/kIANLX16f4rLOzUzZtQ0MD8Xg8Vrsoi0AgQACQLVu2EL/fL34CgQDx+/0kHA6TZDIpfgghKX/TyBOJBBkaGiKJREJVX8t2uh1aX4R0crrp+tI8lOwnk8lx5VHyU85GIpEgPp9vXHoludF6t0qup02t5oxV3MslZ/T4o+ZnPnGGpg5ouaema7SPMIN7NDyg6SNWrFhBPB4PAUCcTidZsGAB2bhxY4oNj8dDli5dqmpHeAaAABBtpHNm7dq1xOPxEJfLRVwuF/F4PGTt2rXjfO/u7pbVo20ntXwIIWRgYED0Ve7jdrup2mPp0qUp9ZVMJonL5SLr16/PqL0z4Uw4HB73bvX7/WTLli0EAAkEAkQLu1yMkzCbpLZkV1tbmx1nFFBdXS37kS7TMTAwMDBYh4ULF2LFihVYvHgx1q9fj9/+9rfw+XyYO3euoQ1Ezz77rPi73BE3XV1dWLx4MRYsWIDe3l6sXbsWDQ0NWLx4ccq5hF1dXWhubk7RW7BgwTg9JQj5NDQ0KKZ3u93YuHHjuM+tt94KAFS7xfv7+7F69WpTYobNRElJieI7lha75FKdx+PB9u3bZZ95vd6cX8Y7d+5c2O180whLL4QQnHPOOVi6dCmAzJaAQqGQaDuTaXjBjt6lOuF3mqW6UCiEsrIy1aU6YapZaalO8FNpqQ4A4vH4uKU6AIhGo3m/VJdMJlPa1IolIKn9bHOPEIKRkRG4XC6qpbpQKAS73a66VBeJRBAIBBSX6gQ/lZbqEokEAIxbqkvuWCYqhKU6ge9WLtUZ6SPSy2SEewIPKioqqJbq0vuIu+66C729vXjrrbcwc+ZMVFdXY86cOTj77LMxf/58XHbZZXj22WdRXV0NQgii0SgCgYAi9wgh+NOf/gSPx4PZs2ejq6sLV155pZgWAFatWoWzzjoLV155JQDgoIMOQmNjIyKRCG688caU8wdbWlpw1VVXgRCCuro6zJs3D4QQ3HjjjeMOXE5vpz/84Q84++yz0dHRgVgshpGREcybN0/Mp729HdFoNGXnubBUt2jRIsyYMQOrV69GIBCQ5d7LL7+Ma6+9Fm+//TYAYGRkROROWVkZAGB0dDSFT1a8nwRZ+v/T6tWrcf/99wPg/9fIjlAAKZ80oTknVYBQW6ojhJC2tjbS1NQk+0wtndUQlup8Pp9l0/ZyS1By+lq20+3Q+iKkk9OVm24X8lCyn0wmZZfU5PyUs5FIJIjf75ddqpOTm90emcr1tKnVnLGKe7nkjB5/1PzMJ87Q1AEt99R0jfYRZnCPhgdqfYTL5SKXX365bH1t3LiRrFixQrQhXapTynP79u0EAFm1ahVZv349AUDWr1+fwhmXy0UaGhrG+Tg0NJSytCWnl0zyS689PT2a9SukT5cL+SjVV1NTEwFANm3apGp/YGCArF27lnR2dorlTM9/7dq1pKmpSVwqXLFiRUZ9hBmcEcYNNEt1u+TAqa+vj7hcLmp5tiAMnGgajoGBgYHBfAjvj56eHip9j8dDWlpaVHWEQYTP5yOEEOJyucal6ejoIACIx+MhHR0diu8vWj0lGEkvpOnu7qbOR6kehcFSZ2cn6evrI21tbTmdsBCg5/07IQdOg4ODmg3R0tIyLkC8oaGB+p/FCmRrxikSiZgy4yS1Q+uLkE5OV+5bo5CHkv1kMjmuPEp+ytlIJBIkGo3KzjjJyc1uj0zletrUas5Yxb1cckaPP2p+5hNnaOqAlntqutJ6LC9PkpISQkpKkjIfa+Xl5fR+EkLIxo0bCQAyMDBAVV/SGScl7nk8HnL88ceLfwuzN+mcWb9+PVm6dClxu91iIHZTU9M4X3p6emT1BgcHqdpp/fr1pKWlRTafdH1hhuzyyy/XxRlh4CQ347R06dJxsu7ublPeT0b7Kz0zThMqxqm9vR1+vx+9vb0AgObmZvHqlPT74jo7O7Fy5UqsXLkSLpcLg4ODaG9vF/VzCatjnIaHh1FZWWlKjFNlZaWhGKfddtuNOsZpypQpqjFO27ZtQ0VFhWqMU2VlpWqMk7D1WVpWQghKS0sLIsbp66+/FtvUitiZoaEh0X62uUd2xDjtscce1DFOLpdLNcZpaGgIJSUlqjFOlZWVqjFOkyZNko1xqqmpKYgYp2+++Qbl5eWWxzhVVlbiyy/5PiISkS9TOCxfppGRUVnuhULy3AsGU2OcdtttNzgcDgQC2jFOQh8xY8YMAMB7772HSZMmiccRBINBFBUVgeM4bN68GW63WzbGKZ17mzdvRn9/P37xi1/A6/VibGwMBx10ENatW4eHHnoI3//+9zF16lQx1mjevHkA+POGnn32Wdx4442YNWsWNm/eDKfTiYqKCjQ0NIhxTYFAAE8//TTuuusu1NXViXpq7XTooYeisbERsVgMX3zxBZ5++mnccccdmDVrFnw+H8rKyhAOhxEIBNDc3Iz58+ejo6MD27Ztg8PhoDqOQIBcjNNRRx01jk8ATHs/SftDFuM0wZGtGSd2HEGqP+w4AnM4YxX3cskZPf6o+ZlPnKGpA1ruqeka7SPM4B4ND9T6CLkYJ0Fn7dq1BAAZGhoihGgfR7BixQrFbf2LFi0iQ0NDZGBggDQ1NY3jDyFEPMZg1apVop5cLJ1UT6lepPmk14uQvrOzU5QJRxYIfunhjNqMk5xMmHGy+jgCJTt6Zpx2ueMICgHCrIH0UMd0mdXyXOTJysTKlG++sDLtmmVatmwZbrrpJnzyySfjdFesWAG3242amhoIULPd1dWFhoYG+Hw+cTaMEIKWlhY8+eSTCAaDqK2txbp169Dd3T3Ojs/nAwDMmzcPkyZNwrp167B27dpxeUr1lHwR0nd3d4/zW0g/d+5ccByH888/H/39/diwYYNYVr31S8sbJd1sc48WE2qpbqJA+McSwHHmHsef/lxOX8tGuh1aX4R0wu9q+tI8lOxr+a9lgxAiTtfSyPWUNRtyPW2aidws+0a4l0vO6PFHzc984ozaMyPcU9I12keky41wj4YHan3E5ZdfjjVr1qCurg4rVqyAx+PBhx9+iIceegj9/f3YuHFjSnt7vV709fVhZGQEFRUVKX2H1+vFjTfeOI4bLS0t6OrqwjPPPIOLL74YK1asQGtrKzZu3IjGxka4XC588skn6OjoQENDA+bMmQMAsnperxcrV64U9ZTqy+l0iun7+vrQ0NAwLr3H40F3dze6urqwdOlSEELQ19cHAOIRDxzHiWc5qXFMKK9c/yqXxqz3k1HO0IIjerQZLEUwGITT6cSsWbNQVFQEYOdomBD5GKdQKKSw1h9UPNJe2uRGj7SXi8mIRCIK8QthxesU5GIyWJlYmViZWJnyoUx33HEHnn76abzzzjtwuVyYPXs2brvtNsyYMUMsk8fjwVtvvQU5nHbaaXjhhRewefNm2TLNmTMHbrcbGzZsQCgUQm9vL+644w689dZbCAQCcLvdOOuss3DxxRenlOmVV17BDTfcIOrNmDEDra2taG9vTymTdFZMgM/nwyuvvIK7774bGzduhN/vx4wZM3D22WfjV7/6Faqrq3HeeeeJcUBymDFjBt566y3VdopEIqipqcHTTz+N4447TmynadOm4YEHHsCxxx4r6s+cOROrV69GY2Oj5dxTinFKJBIYGBhAIBDQPAyTDZzyCMLAyefzpQT3mfntM5lMIhqNpgT4Gf02KbWjZ8YpGo0qnoKe/g1AyENp9gfggyil5VHyU2lmKRaLobi4OCW9kpymbrIp19OmRuVm2TfKvVxyRo7bSv6o+ZlPnBGgVgc0dqzqI6Ryo9yj4QFtHyFAq77U7BjhDG1ZteT5xBktuRnvJ6OcCQQCqKmpoRo4saW6PITQyaTLlHT1ysfGxlJ2FCnpa9lOt0Pri1z+SvqCrpJ9QoiivXQ/lWyEw+FxnYGaXMlOruR62tSo3Cz7RrmXK84ocVspvZqf+cQZmjqgsWNVHyGVG+WeFg/k/FTSpa0vNTtGOEPjO408nzijJTfj/WSEM0rP5cCCwxkYGBgYGBgYKMFmnPIQ6YFzZk6FCrYzDb5Lt0PrizQNzXS78FGbPk4vj5KfcjaUfKH1MddyPW2aidws+0a4l0vO6PFHzc984oycv0bsWNVHKNnR4yMND2j7CAFa9aVmxwhnaMuqJc8nztDIzXw/6fFFqcxyYAOnPITVB2AKwZkcl9kBmIIdvQdgCv7SHIAZDoc1L/kVDqgUplpTDrdLJDD27bfgvvoKzpISJEdHMbp9OzA2Bi4SASIRlBCCBMchGg4DiQSQSIAjBKV2O+J2O2LxOGCzATYbbMXFKHO5EOU4RG02oLQUpKQExdXVKJsyBWG7HVHJP2C2DsCUtqm0ncw6WFFqP9vcI4SIbUxzAGY4HNa85Dcej6te8iv4qXQAppAu/QBMIbalEA7ATL/o2IoDMI30EellMsI9gQe0l/wKecgdkstx/AGYHMel1Jce7in1e4QQkVdW9xGZcq+iogLJZDKlDqw4+NeM95PUhhUHYLLg8DyCEBzu9/tTgtNyOZORzTwV5ckkOL8f5OuvgW++AbZtA7ZvB3w+Xr59O+D3g3w7BAQD4ELDiPuGYRsZRlF4ZFweViNZUopkRRVQXQ2bqxq2SbUgkycDkycDU6aIP7nddwfZYw9g+nRAEhxasO1kkjyffDFLnk++mCXPJ1/MkueTL2bJ88kXs+RW2A4Gg3C5XCw4vJAhjOiV/jYqB/hgxNLS0hQdOX0128IsgNQOjS8p6eJx4IsvgK1b+Z9ffgl8+SW4L79E8vMvEd/6JbBtG+z+b8ElE1AqEZf2EwCK03QIxyFmL0OsqBTx4jLYykpRNaUU/rESjMWKALsdRcVFmDKtCMGRIoyEbUgmE7ARYNrUJEZHkgj6k+AScXCxKCZVjIGMjmHUNwZ7fAzF8TAcCX5WxBYZgy0yBgx9M85HJd8Tk6agaO89gT33BPbeG5zbDdTVAW43/6moUK1j2jY1IjfLvh47ipzRyFOqq2RbzR4Nt5XSa/lpVnuYIdeqA1o7VvQR6TDCPRoeyPmppEtTX2p2jHKGpqw08nzijJocyPz9pMeGlh0lsIHTLohIJCL+c1tmhxDA7wc2b+Y/W7bwPz/9DBjcgtgXX6B4+1fgFCY8bQAcabJEtQtBxxTEaqaATJoM54xaFE91wWergHPGNBRPnQQ4nUBVlfghlZUIJJNw7rYbHDbbOJtOQoBAAE6nU/wHqgZQRfg7oAR5BYAKGT+d0j/icWB4GAgGd34CAZDtQwgOfoOxrd8i9uU3qAx/C1fsG/g+/BIV/q1wJCMo2v4NsP0bQOE8mMSU3WD7n1ngDjgAOPBA4IAD+M/06dptYQKywhmT0tHoaukYfW51O5iJbPhqRh6FwhktHcaZ7OVhtZ9s4MRgHNEosGkT4PUCg4OA14vkwCCSg5/A/tlmfuCQBg5AmVTgcGBs8h4IVu0BTNsdxXvvDuf+u8M2fXeQadMwXF6OKrcb3JQpKHI4UJNmjxACeyAAu9MJyH1z2DEwkn1mNux2oKaG/0jAgR9gOdPUawT/hob4WbetWxH66DME39uC2MdelH/lxZTAIODzoeibr4Fvvgb++c8UG8lqJ7iDDkTZAQcAhx8OzJ3LD6iK0+fcGBgYGBjMABs45SGM7AjQI9e1ayEe52eKNm0CNm0Ct2kT4h9tQuK/m1D9xRZwJJmSzobUMy5ClbthbOo+sLlnoGb23iB77Y3RSbWo2G8/YK+9gMmTUWKzYUq6jzvqIREIgDid4FQOpksvj1xZhYDedD1pYC+NXKm+DMsBoLaW/xxyCCpOAiol+gQAfD7A60XkvU0IvP4Bpn77IYbf+BAVXw6gKBgAXnsNJa+9Btx3HwAgWewA991DgDlzgCOOAI48Eth3X4AzdgghLWfMspMuk6bRylOrveVsKvkoPKf1R83PfIoPkfPXiB3atlaqR7129PhIwwPaPkKAVn2p2THCGdqyasnziTM0ciN9hBEb6XKlMsuBBYfnEYTg8Lq6OupddXqvHvD7/SmHgwm7FoKffgrbxx+jaNMm2DZtQsnmzUh+9F9gYABFiZiizzFHOWyzZiA4eV/AvQ+K998TmDkDyX32gW3GDFROnZpynQIh/C6S2tpaxZ0YQpkIIYhEInC5XCgtLZUtk91ux7Zt21IOqkzfMROJRFBSUiK7Ywbgd4bY7XaMjo6KMo7jUFxcDJvNJu7oAvLsiohIBM6vv0botbdhe+d1jP7zQ5T/5x1UxFJ3oABAonYSkocfhuJjjkF03jyMHnigGJCuVaaRkREMDw+LnDGbezS76qLRKKZOnaq4C0i4TkFo74qKClRVVcmWqaysDENDQ7DZbCJnpGWKxWIiZ4RdddIyEUJQVFSEqqqqcbvqhN1HenbV5eJ6EofDgW+++SblZHyzr1wJhUJiPQq76vSWySj3BB7U1tbC4XColkmrj+A4flfd8PAwEomEWF96uKfU7wkv+JqaGsUdaPlyNU5FRQV8Pp/4TpK2k1n9ntE+QlqmaDSKoaEh0QbtlSvxeByDg4NUweEgDHmDQCBAABC/30+SyaT4IYSk/G1YnkyS5FdfkeRf/0qSd99Nkr/4BYkefTxJ7jaNEH7RSPYTtZeS7XscTIYXnkpIeztJdK0myVdeIcnPPyfJRMJcHy2U55MvlpUpkSDJgQESemAtSV52GfEddBSJ2ErGtWmitIwkFy4kyRUrSPLNNwmJx/O3TBNAnk++sDLRl2nFihXE4/EQAMTlcpGGhgaycePGFH2Px0OWLl1KlSf4SeQUG1L9tWvXEo/HQ1wuF3G5XMTj8ZC1a9eOs9Pd3S2rR9tOavkQQsjAwIDoq9zH7Xar2r/88suJ2+0mAIjH4yHr168Xn7lcrpS/84V7fr+fACCBQIBogQ2c8gjSgVPG2L6dkFdeIeTeewm58EISPfJYEnFOVh0gBZ17kG9mLyD+M35BErffSchLLxGyZQshicQ488lkkoyMjIjko4WedDS6WjpGnxstXy6g6evYGIn/43Wy9de3kOT/nkrCzqnj27+mhpBTTyWks5OQrVv12TfLTxPSMc7QIRu+mpFHLjnT0NBAXC4X6ejoIOvXryd33XWXOIjq6+sT9TweD2lpaVHNM5lMkkceeUQceLS1tY3Lu7OzU3zW19dHenp6SFtbGwFAOjo6RFu0ekqgTd/X15fy2bhxI7nrrrsIANLU1KRov6mpibhcLtLZ2Un6+vpIS0tLSp25XC7S09Ojqy30wKgNNnAqUBgaOI2OErJxIyF/+hMhv/41iR2/kEQm7a44OEpyHInsNYN8c+QPSOiiNpL80wOE/PvfhASDunxNJpPE5/MZ6tBo09HoaukYfW60fLmAbl+TSZJ89z3yedvtZODAH5BkVdV4nsyeTchVVxHy2mskGYuZUheMM/mDbPhqRh654owwuBgcHByn4/F4SENDg6grHTgp5ZlMJslxxx0n6rpcLlm7LS0t49K2tbWJ+tL8tPSUYDR9MpkkM2bMIG63W1HH5/MRAKS7u1sxT62BU644o2fgxILDCwWEAJ9+Crz7LvDOO0i8/S5ife/CsWUTbJIAbWmDJvbaB8nvHIjYfgehfN6BwEEHAfvth9FYDJMk2+8ZdjFwHLiDD8L0joOAjkuAeByx1zfi8wd6wL34PPb64g3Y3n4bePtt4IYbgMmTUbZwIfCTnwCNjSmHdTIwTES0t7ejra0Nbrd73LPVq1ejt7dXlz2/34+XX34Zq1atQl1dHbq6utDb24sFCxaIOl6vF7W1tePSLlu2DI2Njbr1lGA0/eLFi7F582YMDAwo6gwNDcHj8aChoSFFXltbC6/XK/7t9/vR3Nws1uOyZcvQ1tam6XvewOiojsF8iDNOX37JzwKtXk3IRReRyOFHk2iFU3EWaaRiMkkcO5/EL/wlSXZ2EfL664QojJrN+qbJZg/yB6b7um0bIQ8+SL48djEZtjtT+VZdTcjPfkbIM88QEg5nxU/GGfPBZpyUdYVZE+msiJo9mhmnVatWEQBkaGiIEMLPuqSn6ejoEGOCOjo6UpYDpaDVU4KR9EKaBx54QHdbDA4OpiwDCnFVwlKesEzY19dXMDNObOCURxAGTj7+WMhxn5itmGzf+7sk+bMzSHLlSkJeeokkv/hCDNAWPoSQlL+l8kQiQUZGRkhCkkZOX81GMpkcZ0dLX/gI6eR00/WleSjZTyaT48qj5KecjUQiQUZHR8elV5LrKWs25Hra1JA8GiXDz/WSTSdeQJK7py4BJysrSfKnPyXJ//s/koxGLeNeLjmjxx81P/OJMzR1QMs9NV1pPf7ud0lSUkJISUmSlJQkyeefJ8njjyd3/M3LX345Sfr6kqJOSQkhnZ0Jsm3biESWJJdfztt3u3fKfvQjPt8TT9wp+853EiSRSJBt2+j9JISQjRs3EgBkYGCAqr6E4HA17nk8HjJ//nzx76amJgJgHGfWr19Pli5dKgZWY0c8UbovPT09snqDg4NU7bR+/XrS0tIim0+6/vr16wkAcvnll+vmzPr168Xgc0HucrnGBdO7XC7S3d1tyvvJaH8oDJhpBk7sOII8gnAcQQBAhKvER44qbKqI4OPSUvzH4cAxLS0494ILAFi8zR35t32alSkPypRMIvLyv7G983lMffVpTA5vFXWTU6YiedoS2M8+G8G6OiQlvud1mSZiO7EyZVSmoaEhTJ48GU8//TSOO+442TJt3rwZbrcb1dXV8Hg8OOSQQ3D77bfLlmnz5s2YM2cOfvWrX+Giiy5COBzGunXrsHz5cjzwwANYvHixbJkikQieffZZ3Hjjjdi8eTM2b94Mp9M5rkyBQABPP/007rrrLni9XlGPtp2++OILPP3007jjjjuwefNm+Hw+lJWVIRwOIxAI4Lvf/S48Hg/++te/UrdTOBzGmWeeid7eXlx88cW49tprxXaaNm0a7rjjDvzwhz8U9WfOnInVq1ejsbHRcu4pHUeQSCQwMDDAjiMoNIgzTv/5j2XfPhOJBAkGg6bMOEnt0PoipJPTTdeX5qFkP5lMjiuPkp9yNhKJBBkeHpadcZKTm90emcr1tKmpnEkkSOzV10jyootIsDR1t2Zy//1J8ve/54+rMIF7ueSMHn/U/MwnztDUAS331HSN9hFmcJuGB2p9hMvlIpdffrlsfa1duzZl2U064ySX54oVKxS39S9atIgEg0EyMDBAmpqayNDQ0Dgfh4aGCACyatUqUU9YipJ+pHpK9SLNJ71ehPSdnZ2iTDiyQPCLhjMbN24Uj25InykT6jb9OALpjFOm7yejnNEz42STH04x5BLctGniIWPSyzHTP0blwiFuavo0tqV29PiSSCSo9aUHzinpp5dHyU8lG/F4XJfc7PbIRK6nTU3ljM0G+1FHgLvrLlQFv8A3f/wzPjp4MWL2UnAffQTuyiv5y4pPPRXc888DO9Ib5V6uOKPXHyV5PnGGtg6MtJPeejTEPRN5oKS7bNky3HTTTfjkk0/G1deKFSvgdrtRI7laSS3Prq4uNDQ0wOfziYdzEkLQ0tKCJ598EkNDQ6itrcW6devQ3d09zkefzwcAmDdvHiZNmoR169Zh7dq148oj1VMqp5C+u7t7XL0I6efOnQuO43D++eejv78fGzZsEMtKw5nFixdj8eLF6OnpQV1dHTWX9LS3VZyhBdtVx8DAkBmKizHl56dgys9P4S817l6HT675I9xfvgY88wz/2XNP4Oc/B9fUBBx8cK49ZmBQRVtbG9asWYO6ujp0dHRgzpw5+OCDD/Dwww+jv78ffX19Kfperxf9/f0IhUKorKxMeQl7vV6sWLFiXB6tra3o6urC008/jUsuuQQdHR1obW1FX18fGhsb4XK54PV60dHRgYaGBng8HgCg1pODy+WiSr9u3Tp0dXWhpaUFANDf3w9CSEr55PLp7e2F1+tFfX39uJ2HtbW1qr4VFDTnpBiyBlMPwFSAMCUpTFNm246edDS6WjpGn5tVT9mA1b4atT/y7/fJew2XkrHK2pRzxJI//CEhPT2EUNpjnDEf2fDVjDxy3c90dHSIh146nU7S0NAgnu0kQHgu9znzzDOJy+VSzMPtdpPjjjtOlPf09IgHb2LHCd1yh1rS6sn5RJNeOLRS6aN0lpNw/pXcRzg0U+4cJ2GpLpec0bOrjgWH5xGE4HC/358S3GcmyI57v6R3u2XTjp50NLpaOkafm1VP2YDVvmZsf2wMI48+g63XdGG/L/62U77//sBFFwFnnglUVpqSP+MMHbLhqxl5FEo/o6XDOJO9PIzaCAQCcLlcVMHhbOCURxB31dFE9TMwFCBCb36EwV/fg4P7HoRtNAQAINXV4M4+G7jkEkDmwEEGBgYGq6Hn/cuCw/MQVo5lCSEpN7pn246edDS6WjpGn5tVT9mA1b6ayZnkfnvgkFfugu3Lz/F52x3YWr4vuGAQuPNOkH33BRYvBt54w3D+jDN0yIavZuRRKP2Mlg7jTPbyyIQztGADpzwE2bHrQvjIyTKRJxIJTX0a21I7enwRdh7R6At5qOmnl0fJTyUb0p0uNHKz2yMTuZ42tZozuuxUVWH6il9iz+H/YHPnC/jPjBPBJZNAdzdw+OEgxxwD8uyzIDv0c8UZvRxWkucTZ2jrgJZ7mdSjldyj4YGZnNGyo5czZrV3PnHGrPa2ijO0YLvq8hBz586F3c43jbBNkhCCc845B0uXLgWw8yC4kZER2YPghoeHZQ+CCwaDCIVCom0jB4zF43GEQiHRTlFREaqrqxGNRmUPtxsbGxMPTSOEiL+Hw2HZw+2EMhHC7+IoKytDaWmpbJnsdrt42Jmwni0tk2AD4HeUpJdJQDwex+joaEpZASAajWJsbEy1TMDOg+C0yqSnnaT/yFrtJG1TaTvJHQRH007pZZLaN5t7rsWHA4sPR/CDD7Dl0lXYb+MaOF59FXj1VSRmzcLYL3+JkVNOgcvl0iyT0N52ux1VVVWyZSorK0MkEkEgEBDbWVqmWCwm+lleXj6uTNKXnbSdCNk52A4Gg7LtpPT/ZFY70XLP4XCIfBfqwCj3zO4j0stkhHsCDyoqKuBwOFTLpNVHcByH6upqJBKJFM7o4Z5Sv0cIEX+3uo/IlHsVFRWIRqMpdWBFv2fW+0mQpZdJ6QBMqZ9aYDFOeQRhjdXn86UEhwsDp3QYkSeTSQQCATgll/zK6WvZJoSk2KH1RUjncrlk60CqL81DILgchGD69PNA0v2UsyG85Kqrq8cFbcrJaeomm3I9bWpUbpZ9Wjsjm77EhxfciYP/uQplEb7DTO6xB7i2NuC884CyMsU8zeaMHLeVOKzG7XzijAC1OqCxY1UfIZUb5R4ND2j7CAFa9aVmxwhnaMuqJc8nzmjJzXg/GeVMIBBATU0NCw4vNGRrV108Hofdbs94V50RO3rS0ehq6Rh9blY9ZQNW+5ozzgwPI35vF3y/uRlT4l/xNqZOBferXwEXXgjIdG6MM3TIhq9m5FEo/YyWDuNM9vIwaoPtqitQsF11DAzjkRgZw1uXPIBpD3Vgz9hmAABxucBdcglw6aWAwrd1BgYGBlqwXXUFDivHssLUcKZ5GLWjJx2NrpaO0edm1VM2YLWvueaMrbwE+97yE+we/C8Sf3wQw3t+B5zfD1x7LciMGcC11wI7YiAYZ+iQDV/NyKNQ+hktHcaZ7OWRCWdowQZOuyDMIr5RO3rS0XZ8VjwvhM5MgNW+5gNnbCXFKPr5majc/AH62tZgoPRAcIEAsHw5MGMGcN11/JUvjDNUyIavZuRRKP2Mlg7jTPbysNpPNnBiYGAoKHBFNtR3LEZd6F1sun4NcMABgN8PXHMN4Haj5JZbgB27ahgYMsHKlStRX18Pm82GGTNmYOHChejv70/Rqa+vR2trK5U9m80GjuPG2RCwbt061NfXo6amBjU1Naivr8e6desM6ylBK73X6xUD+oWPzWZDTU0NbDYb6urqqPJpbW1VLKsW2tvbxUuC0+++q6mpGXcXXjbBBk4MDAwFCa7Ihn2vWgy8+y7+dekTGHDsD87nQ9n11wOzZgF33AFIjpJgYNCDxsZG3HjjjViyZAnWr1+P5cuXY2hoCPX19YYGA88++6z4+5o1a8Y97+rqQnNzMxoaGrBhwwZ0d3ejoaEBzc3NWLlypW49JdCkd7vd6OvrS/ls3LgRt912GwBQXdbb39+Prq4uDA0NaeqmY/Hixejq6kJ7ezv6+vowd+5cNDY2Gh6EmQ7CkDfI1iW/8XjclEt+jdjRk45GV0vH6HOz6ikbsNrXQuFMbCxOXm55hAR2m0UIwH/22ouQ++4jJBajtsc4kz955KqfES6rFS70lep4PB7S0NAg6no8HtLS0qKaZzKZJAsWLBB1XS6XrN2WlpZxadva2kR9aX5aekowmj6ZTBK32614wa+Anp6elIuP0y/0VUMymSTffvstAUC6u7sV/Za7KFhqwwhn9Fzyy2ac8hBWbksVplwzzcOoHT3paHS1dIw+N6uesgGrfS0UzthLinBs5+mo+vQDfH1dF76w7QF89hl/9tMBBwBr1gDJJOMMsuOrGXnkqp9pb29HW1sb3DvuTpTqrF69Go2NjbryDAQC2LBhA1pbW9Hc3Ay/34/e3t6UNF6vF16vd1zaZcuWobu7W/ybVk8JRtMvXrwYXq8XPT09qvbdbjeWLVuGzs5OTV/SZ5A4joPf74fH40FDQ0PKs9ra2hS//X4/mpubxeVGYbYsE85QQ9eQjMFSZGvGyefzmTJ7YMSOnnQ0ulo6Rp+bVU/ZgNW+FipnAl+HyQsLbyXfcJN3zkDV15Nkby/jTBZ8TSaTxDc0RJLDw4SEQoY+yeFh4tu6VbeNcel09CE+n2/cTIlafUlnnJR0Vq1aRQCQoaEhQgg/Y5KepqOjgwAgHo+HdHR0kL6+Pll/afWUYCS9kOaBBx6g5oxcPco9l0KpDgcHBwkA0tHRQQjh68/lcpHOzk7S19dH2traCADS19dnmNt6ZpzYwCmPwAZO+nXZwIkNnLR0A1uDhFx7LRlzVIoDqOiCBST59tuG8mOc0ZHH1q07B625/IRC1HXR19eXskynVV80AyePx0OOO+448VlTUxMBMC5NT08PaWlpIW63W1zqampqSvFFj54S9KTv6ekhAMjll1+uizNmDZx6enqIy+UiHo9HlAkDTylcLhfp7u7OysCJLdXlIQix9hJFGnku8mRlYmWywnbV9EqQq6/G1pcH8Lz7F4jBjuING4A5c0DOOgtky5aCK1M+yGl18wF6yjRz5kwAwODgoKKu8EzJvlQ2ODiI/v5+HHLIIfB6vRgcHMTcuXMBAN3d3Sn6CxYswKpVqzAwMIChoSF0dnaiv78fdXV18Pl8os2GhoYUvVWrVsnqKZVzwYIF6OzslE3v9/tFPZ/PJwaSd3R06OaMWt1s37495Xm6rs/nQ2NjIxobG7F06VJs3LgxRbehoYEqPyu4yi75zUOwS37ZJb+70iW/WtwjhGBkZCTjS34n7V+C495bibefvQSTbrkc7r5ngYceAlm7FvGLLkLx1VdjxGZjl/yaeclvMgl89pnYR0QiEdkyhcNh2TKNjIzA5/OhoqIihXuhUEiWe8FgUGyDkZERTJ06lb/kNxYDkfiv1kdUV1fD6XTiL3/5C+bNmweOS73k97nnnsPZZ5+NLVu2YO+99wYhRLz8Vo57jzzyCADgzjvvxJ133plSh48++iiOPfZYeL1eXHbZZbj99tvF67bKy8vR0tKCE088Efvssw8efPBBHHfccbjuuutw//33i+8FjuOwZMkSNDU1YfLkyXjwwQdx9tlny7bT5s2bsXz5ctxxxx3YZ599EI/HxfQnnngiZsyYgbVr1+Kss85COBzG/PnzQQjBn/70JwDQdcmvgJGREZE7//nPf3DEEUek1IHNtnP+ZvPmzbDZbHjnnXfw4x//GHPnzsXAwAAmTZo0jn/V1dUpMoGz2bjkF4QhbyBdqksmk+KHEJLyd6byRCKhqU9jW2pHjy+JRIJaX8hDTT+9PEp+KtnQ+zG7PTKR62lTqzljJfes4MyXz/6LvF56LBGXc6ZOJclVq0gyGtXktpI/SvJ84oy0DszgnprcaB9hBvdoeKCku2LFCgKADAwMjKsvj8dD3G63KPd4PGTp0qWKebrdbtLQ0DCOG0uXLiUAyPbt28Vlq1WrVo3zcWBggAAgGzduJENDQwQA6ezsHFceqZ5SOYX0q1atGlcvQnohTkjwT7CnlzNCmdavXy+rL/gi195ut1usUzn7LpdrnF3pUp0RzuhZqmMzTnkKuR07NHpacgDit5T0G6712CY7RulSOzS+EMm0qJZ+up7Shaty5ZHzU84G2TFLkL4LQ0kuV6Z8kNO2qRG5Wfb12MkGZ3b7waGo8f0NT7Y8h9mPXY66bZvAnX8+yN13g6xcCe7EE2XtKPmj5Wc+cUbr/4bWjhV9hJyverlHwwM5PwVZe3s71q5di1mzZqGjowNz5szB4OAgVq9ejf7+fvT19aXY/OSTT9Df3y/bZ3i9XqxYsWJcPueffz5Wr16NtWvXorW1FR0dHTj//PPR39+PxsZGuFwueL1edHR0oKGhAfX19QCAjo4OtLa2oq+vT1VPrl5qamo08/F4PFi3bh1Wr16NlpYWcByHt956a1yfKD3LSa39lDim1DY9PT3wer1ob2/Hhg0bUp7V1taK+crZFWBWf6UIzaEVQ9bAgsP162rpGH1uVj1lA1b7uitwxr8tQsgdd5DRshpxBir5/e8T8tFH1HkwzpifR64509HRIZ5J5HQ6SUNDw7gAaumZRemfM888k7hcLsU83G53StB4T08PaWhoIC6XiwAgbrdb3EkmBa2enE806VtaWhTLJOhrwWhw+G233aaYb1NTEyFE/hynbAaHc4ToiIhisBTC7cx+v19c5zYbhPAXIDqdTmMj7Qzt6ElHo6ulY/S5WfWUDVjt667Emf7eIXx85nVY9NU9KCYxwG4HLr6Yv85lx/8k40z28igEztDoMM5kLw+jNgKBAFwuFwKBAKqrq1V12a46BgYGhh3wNNRi8dZb8fR1b+LT2T8A4nHg1luB//kf4P77AUnQOwMDw64JNnDaBWHWtwWjdvSko9HV0jH6PN+/AUphta+7Emc4DjjhohnYq/9ZDN79AgaK9gO2beNPID/0UOC11xhnsphHIXCGRodxJnt5WN4fsqW6/IGwVEczVcjAwJAdfPRuDC833Y0zvMtRmdhxzMBZZwErVwJTp+bWOQYGBlOg5/3LZpzyEFaOZQkhiMViGedh1I6edDS6WjpGn5tVT9mA1b7u6pzZ/5BinP/fX+E/z20CzjsPhOOABx8E2W8/4N57gR1nOjHOmJ9HoXBGS4f1M9nLIxPO0IINnHZBSA/by4UdPelodLV0jD43q56yAat93dU5w3HA3JOmAqtX4w9nvIZ37R5wfj/wi1/wy3dvvKHbz1wjG76akUehcEZLh/Uz2cvDaj/ZwImBgYFBBy544DAEenpx04y7ES13Av39wBFHAC0t4Hy+XLvHwMBgMdjAySJ0dXVh5cqV6OrqQnt7O/x+f65dYmBgMAkHz+bwq00XwrbpY2w57iyAEHD33Yeqww4DHn+cPwmKgYFhQoKdHG4Burq6sHjxYrhcLgCA3+/H0qVL0d3dnVvHdkB6N1Au7OhJR6OrpWP0uVn1lA1Y7SvjzPjnRUUAN30qHI89gN+efi5+9toF2PebD4DTTwcefpiPf9pxYWw+Ihv8NiOPQuGMlg7rZ7KXh9V+sl11FqCxsRE9PT2asnSwXXUMDIWL/3s6Cu6mlTip7zogGgXKy4FrrwUuvZQ/SJOBgSFvscvvqvN6vaivr0d/f7+qXvpyWm9vryn5Dw0NobW1dZyMFlbvWohEIqbskDJiR086Gl0tHaPPzaqnbMBqXxln6J6f9ONiLPjb5fi65x38u/xYYHQUuPxyPnhcoy/KNrLBbzPyKBTOaOmwfiZ7eWTCGVpMqIFTa2srmpub0dnZqTloEgY2bW1taGlpQUdHBzo7O00ZPHV0dKCrq0scvLW3t+fNMh0AhMPhnNrRk45GV0vH6HOz6ikbsNpXxhm65+FwGFOP3g+Bp/+G9sn3Y6SkBnjrLX7wdPXVQCSi6Vu2kA1+m5FHoXBGS4f1M9nLw2o/J9TAqbOzE93d3Vi2bJmqXn9/P3p7e9HS0pIi7+joQHNzc8Z+NDQ0oLu7G/39/eJN1W63O2O7DAwMhYHGhRyu/uQcdF36ERJNi/mznq6/Hpg7F9i4MdfuMTAwZIBdcuG9s7MTDQ0N4+TC4Ka3t1d8vm7dOqxZs0bTZkdHR0p6r9cLn8+HG2+8EStXrkR/f79mjBMDA8PEQWUl8KsVu4GQNbj6gGa0bbkQVe+/Dxx+ONDezl8cXFKSazcZGBh0YpccOPX29qKpqUn2mdvtRk9PjzhwampqUtSVg9/vR3t7O/r6+gDwA6rW1lbU19dj3bp1umxZBbtJgapG7ehJR6OrpWP0uVn1lA1Y7SvjDN1zOTnHAZe/0YTfnH8sFv7lYpw8/ATw+98DzzwDPPggPwuVA2SD32bkUSic0dJh/Uz28rDazwm1VEcLr9eLuro62We1tbWa8VFatuemdYRutxsdHR3UAeJWXlDIcRwqKyszzsOoHT3paHS1dIw+N6uesgGrfWWcoXuulq66GrjjsSlIPPI4tnc9ieSUqcCHH/IHZ153HRCPa5bNTGSD32bkUSic0dJh/Uz28siEM7TYJQdOWtCzAy4dHo8HXq93nLyvrw+LFy+msmH1roVwOGzKDikjdvSko9HV0jH63Kx6ygas9pVxhu45Td4//CEwaempWHnWh+ipWcwPmK65BjjqKGDTJs3ymYVs8NuMPAqFM1o6rJ/JXh6ZcIYWhTNHaBK0TvB2uVyyAx896OzsRHt7OyZNmiTK2tvbxQMxtRAIBFL+5jgOhBCUlJSgRBITIcjToSUfGxvTtKNlI90OrS9CutLSUk19aR5K9gGMK4+Sn3I2hK2r6emV5HrKmg25njbNRG6WfSPcyyVn9Pij5me6jbaOWvzxfx7HhZf9ELfHfwHHG2+AzJ4N3HorsHQpOJvNUi5p1QGNHav6iHS5Ee7R8IC2jxCgVV9qdoxwhrasWvJ84gyN3Iz3k5qNSCSCyI7drVJ5+ntXDbvcwIkGtAMcJQhLc0axzz77yMrb29txxRVXAABKSkpQVlaGkZERxCXT/GVlZSgpKcHw8DCSyaQor6ioQHFxMYLBIEKhEACeNFVVVbDZbONI43Q6kUwmMTw8LMo4joPT6UQ8HkcoFBLtFBUVobq6GtFoNGUbqN1uR2VlJcbGxkSiCgMSgN8yGo1GRf30MhFCEAqFUFZWhtLSUtky2e128UJHYapVWibBBsC3a3qZBMTjcYyOjqaUFQCi0SjGxsZUywQADocD5eXlmmXS007SjkGrnaRtKm0n6WWXNpuNup3SyyS1n23uEUIwMjICl8ulWSahve12O6qqqmTLVFZWhkgkgkAgILaztEyxWEz0s7y8fFyZCCFIJBIAkNJOhBAkk0kQQhAMBlXLBADNizkcdfTpeP8/R+DAm89GyT9fBc4/H7FnnkHxAw8g6nLpbida7jkcDpHvQh0Y5Z7ZfUR6mYxwT+BBRUUFHA6Hapm0+giO41BdXY1EIpHCGT3cU+r3CCHi71b3EUrtRNtHVFRUIBqNptSBFf2eWe8nQZZephUrVmT0fgYm6Mnhfr8fNTU16Ovrg8fjkX3W2dk57jgCgD/he2hoSAzuziaEk0u3bNkCp9Mpys2ccUomkwgEAnA6nSL5jX6blNrRM+MUCAQUB6fp3+CEPGwK38ABvk2l5VHyU86G8JKrrq5OSa8kp6mbbMr1tKlRuVn2jXIvl5yR47aSP2p+atVX95okBi66HVcEr4QtFgWmTgV54AHgxBN12aGVa9UBjR2r+gip3Cj3aHhA20cI0KovNTtGOENbVi15PnFGS27G+0nLhtqM0z777LPrnhyuBun9cUqora3NjjMKqK6ulv3ILRsZgcPhyKkdPelodLV0jD43q56yAat9ZZyhe260fM1LbDjoj/8PC5wb8c30Q4Bt28CddBJw2WX89S0WIBv8NiOPQuGMlg7rZ7KXh5qNkpISxXcsLXbJpTqPx4Pt27fLPvN6veOuS8k25s2bh6KiIgBI+aZ2zjnnYOnSpQCMT4UODw+DEIJYLAYg8+ndWCxmaAmI4ziMjo5STe/a7XbVMsXj8ZSlEbkyxWIx1TLFYjHZKetIJGLZcolZ0/BCewptavZS3djYWIr9XHGP4zjFdkovEyFEtUwcx6VwRq5MsVhMtUwcx4lLPdJ2AsbHS9D8Px1zDDDrlb3xct/f0PTGb8Hdcze4W29F/K9/xeh998G2336mci+RSKTUgVVLQEb7CDO453A4YLPZqMqk1UcIS34CjHBPqd/jOA6hUCivl+oqKyvFukxvJ7P6PTPeT4lEIoUz6WVavXo17r//fvGZMOsoLL/TYJdbqgP4WCGv1yt7DQrHcYrprIawVOfz+WSX6tJhdCo0HA6jrKws46U6qR09S3XhcBjl5eWydZA+9S3koTaFPjo6mlIeJT+VluqE4Mz0pTo5OU3dZFOup02Nys2yb5R7ueSM0lKdnD9qfuqtrzff5LD6B8/i3rFzUBwcAqmoAO65B9yZZ0KuREaWXdTqgMaOVX2EVG6UezQ8oO0jBGjVl5odI5yhLauW3Mylukw5oyU34/1klDOBQAA1NTVsqU4JS5Yskb2Trr+/Hy6XKyeDJimETkb4yMkykcdiMU19GttSO3p8kc6MaOkLeajpp5dHyU8lG9FoVJfc7PbIRK6nTa3mjJXcyxVn9PqjJNdbX4ceCpz6wI9wWOm7+Hzf48CNjIA7+2zgjDP4301oV606MNJOeuvRSu7R8MBMzmjZ0cuZTOrMbBtmccas9raKM7SYkAMnrXOYPB4PFi9ejK6urhR5vlzGK3wLEj5yMqvluciTlYmVKd98yWWZTjyR4Kk3pmPN0h6Q664DKSoCHn0UZO5ckHffLcgy5bs8n3xhZcp+mWgxoWKc2tvb4ff7xdmk5uZm8eqUzs7OFN3Ozk6sXLkSK1euhMvlwuDgINrb22XvsMs25s6dKx4ZL4yECTEnxokdR8COI2DHEeTHcQQ07eRyhXHOUuDtTy7GrfcfjT+GT0fxf/8LcthhCN90E6Knnw5HSQk7joAdR8COI+DojiNQinGS+qmFCRnjVKjIVoxTeuyOnL6WbUJSY4BofRHSlZWVydaBVF+ah1rsQTgcHheLJOenUoyTcNClNL2SnKZusinX06ZG5WbZN8q9XHJGjttK/qj5aUZ7PP44cMOl32LDnmdi2lsv8nmeeSYf+1RZqcu+Vh3Q2LGqj5DKjXKPhge0fYQArfpSs2OEM7Rl1ZKb1RcAmXNGS27G+8koZ/TEOLGBUx5BGDjRNBwDA8Ouh3ffBW69OYk/fqcDtqt/AySTwP77A93dwIEH5to9BoaChZ7374SMcSp0CFPNVqznCtOb0jzk9LVsp9uh9UVIJ6ebri/NQ8k+IWRceZT8lLMhLHOlp1eSm90emcr1tKnVnLGKe7nkjB5/1Pw0qz0OOQT404McvjzrCly/4K9ITtsd+Ogj4NBDQR5+mNqOVh3Qck9N12gfYQb3aHhA20fQ1peaHSOcMYs3ZnHPDM6Y1d5WcYYWEyrGaaLA6hin4eFhxONxU2Kc4vG4oRinyspK6hin4uJi1RinQCAglgeQj1+Ix+O6Y5wIIQUT4+T3+8U6sCLGSWo/29wjhI9xqqyspI5xEupMKcZJ8F0txikej6vGOFVWVloe46TUTqWlEWw/8DAc/nEfeg47A843NoA780xE/vpXhG+8ESgt1YxxCgaDKXVgVYyT3j5CaKdMuCfwoKSkhDrGSamPEGKcIpFISn3p4Z5Svyf0MZWVlQUR4zQ6OppSB1b0e2a8n6ScYTFOExzZinFiV67shPCSY1euZM4Zs+zkE2eUYpzMvnLFqLy7G7hlZQKvn3Q9uOt/B44QkPp6/sGMGarxKuzKFbo+QgC7cmViXLmiZEdPjBObccpDCJ1MukxJV69c+lHT17KdbofWFz36WvYJIbLlkfNTK0+lOjBSN9mU62lTo3Kz7BvlXq44o2RHr5zWd73yxYuBU0+1IxJbjr/5D8f3H/sZuL4+oL4eePhh4OSTDdcBjT9W9RFSuVHu6cnXDM5o2THCGbPkZtgwizNacjPeT0Y4o/RcDizGaReE0g6ObNnRk45GV0vH6HOz6ikbsNpXxhm657ngjN0OjI0BN713Ii44tB+JeYcBPh9wyinAVVcBCksQ2fDVjDwKhTNaOvnEGaOY6JyhBVuqyyMIS3V+vz9lqpDjcrMElC/LUWbJ88kXs+T55ItZ8nzyxSx5NvKMRgkuvRR4d2MUfz/sMtjuvhsAQI4/HnjsMWDq1IIr00RsJ1am/CxTMBiEy+VixxEUGoSBU11dHXVwuNLFkMFgUDb4zu/3Y2RkRLzo1EjwnRA8Ozo6ivLycjHwU+lC3HA4nBIkGYvFMGXKFMWAQqFMhBCMjo5i0qRJKC0tlS2T3W7HF198IZYHGB/4KfipFBxus9lQWlo6Ljic4zg4HA7Z4HBpmYCdQZJKFxfrbSe5y2OV2imRSODrr78W60DaTnouLlYq08jICPx+v2g/29wT2nD69OmKwaxCmQTd6upqVFVVyZaprKwM33zzDYqLi0XOSMsUi8VEzgjB4dIyCXExNTU144LDBf/1BIeb1U4C9wYHbXC7kwh0PocZ118IjIwgOX06Rh54AIl588Tg8K+++irlrBsj3LOij5CWySj3BB5MmTIFDodDtUxafQTH8cHhPp8vZflHD/eU+j1C+KDkyZMnKwZSm9FHmMG9iooKfPvtt7Db7WIdmN3vmfF+ikaj+Oabb0Qb6WVSCw4fHBxkA6dCAwsOT9WX5sGCw5XbiQWHW8cZOW4r+aPmZ7a58dlnwOGHA3e0foSmx08F99//ghQXA7feClx4IcBxLDicso8QwILDWXC4ABYcnoeQfquRypR09cqlHzV9Ldvpdmh90aOvZZ8QFhyup02Nys2yb5R7ueKMkh29clrfzZLvvTfwz38CP/jBAeg7/t+48eBzwa1bB/zyl8C//gWyapVmHdDka1UfIZUb5Z6efM3gjJYdI5wxS26GDdo6yFRuxvvJCGeUnsuBBYczMDAwTEDMnAm89howy1MNbu1aJG++BdhxUTCOOAK2wcFcu8jAUJBgA6ddEBUVFTm1oycdja6WjtHnZtVTNmC1r4wzdM/zjTPV1cB55wGbBjjUP/L/8OWjfwWmTQP3/vuoOv544NlnLc3fjPooFM5o6RQKZ9SQDV9zyRlasBinPALbVcfKpFeeT76YJc8nX8yS54MvXV3A734HPH3vlzj05iXAP/4BACDLlvEPiooKrkxmy/PJF7Pk+eSLWXIrbLNddQWKbO2qC4VCqKioAMdltqtuZGQEFRUVunfVRaNRTJ06lWpX3cjICCZPnqy6q+7zzz8XywOM3zEj+Km0q47jOJSVlY3bVSf4VAi76r766iuxDqTtZNauOp/PJ9rPNveENtxjjz2odtWNjIzA6XSq7qrbtm0bHA5HSjtLd9UJnFHaVZdMJlFbWztuVx0hBC6XK6e76gTItdPf/27HI4+U4YmHExj9xYWo3LG7KDZ/PvDooyjefXdTd9UZ6SOkZTLKPYEHU6dOpdpVp9ZHcBy/q25oaAg2m03kjB7uKfV7hOzcbVcIu+rSd6Nasasu0/dTNBrFtm3bRBtsV90EB9tVl6ovzYPtqlNuJ7arzjrOyHFbyR81P/OJMwCweTPwwANh/PZ/ngbX0gJudBRk773BPfkkf2ULhR2r+gipnO2qy5/ZGYDtqhPAdtXlIYROJl2mpKtXLv2o6WvZTrdD64sefS37hLBddXra1KjcLPtGuZcrzijZ0Sun9T1b8qoqgt5eBzZt+in+9MohKPnJqeAGBoDvfQ/cvfcC555LZceqPkIqN8o9PfmawRktO0Y4Y5bcDBu0dZCp3Iz3kxHOKD2XAwsOZ2BgYNjFMGkS8NRTIdjtwJk3HQxs3Aj88IdANMpHk59/PiBZOmNgYNgJtlSXR5AGh0uX6syEEJ8hXavPph096Wh0tXSMPjernrIBq31lnKF7Xoic4Tgbvv6ag9MJfLE1iVndNwJXXw0QAhx2GLBuHbDnnhnlkUl9FApntHQmEmes9DWXnBGWTGmW6tiMUx7Cyn8gjuNMIb5RO3rS0ehq6Rh9blY9ZQNW+8o4Q/e8EDljs3HYfXd+wumY42x4+XtXAc8/D9TUAG+8AdTXA3//e0Z5ZFIfhcIZLZ2JxJl8fz9lwhlasIFTHkLYFSJ8AKT8nYk8mUzC7/en5CGnr2U73Q6tL0I6Od10fWkeSvYJIePKo+SnnA0hkDA9vZLc7PbIVK6nTa3mjFXcyyVn9Pij5mc+cUauDo4+mmDtWoLTTwce/uYEkDffBPnud4Ft20AWLADuvBNEZ79ktI8wg3s0PKDtI2g4o2XHCGfMam+rOGMFV814P2XCGVqw4PA8xNy5c6mPI1Dawjo8PCy73TMYDCIUCom2jWz3jMfjCIVCoh1hq3E0GpXdajw2NpayLVf4Xek4AqFMhBCEQiGUlZWhtLRUtkx2u13cTit8Y0jfaiz4qXQcAQDE43HZ4wii0ajscQTSMgE7t+VqlUlPOwn/2OllkmsnaZtK20luqzFNO6WXSWo/29wjhN/m7XK5NMsktLfdbkdVVZVsmcrKyhCJRBAIBMR2lpZJ2EoPQDyOQFomQggSiQQAjDuOQOis9RxHYGY70XLP4XCIfBfq4LDDKvDyy8X4979H4a+dBPzf86i49BI41q0DLrkEsX/8A6O33w6Ul1OVyWgfkV4mI9wTeFBRUQGHw6H6/6TVR3AcfxxBIpFI4Ywe7in1e4QQ8Xer+4hMuVdRUYFoNJpSB1b0e2a9nwRZepnUjiOgBYtxyiOw4whS9aV5sOMIlNuJHUfAjiPQKwfU6+DBB/nzMe+9h6D4D3eCu+wyIJEAmTMHeOopYJ99LOsjpHJ2HAE7jkBvmdhxBLsohE4mXaakq1cu/ajpa9lOt0Prix59LfuEsOMI9LSpUblZ9o1yL1ecUbKjV07re7bkWnWweDHw3HPAj37MYe3aS1E1Zw7Q3AzurbeAuXOBtWuB44+3rI+Qyo1yT0++ZnBGy44RzpglzwZnMrUvyM14PxnhjNJz2TRsxil/IMw40Yx4M4HwD5ArO3rS0ehq6Rh9blY9ZQNW+8o4Q/d8InEmkQD+3//jN9VdfjmATz8FTj0V6OvjLwu+6Sbg0kuBDOrRDD/NSGcGZ7R0dgXO5EseRmzoef+y4PA8hJVjWWkcRi7s6ElHo6ulY/S5WfWUDVjtK+MM3fOJxpmiIuCOO4Bf/xp4+WXgo5G9gVdfBc44Y+eo6owzAElsoN48zPAz03RmcEZLZ1fhTD7kkQlnaMEGTrsg5IKjs2lHTzoaXS0do8/NqqdswGpfGWfonk9EzthswPbtwIIFwGtvlQEPPsiPqIqKgEcfBY46ip+NyiAPM/zMJJ0ZnNHS2ZU4k+s8rPaTDZwYGBgYGFSxaBE/Rlq0CHj5FQ64+GKgtxeYPBkQ4p5efTXXbjIwZAVs4MTAwMDAoIn584H164FDDtmxOnfccfzJmbNnA998Axx/PLBqVY69ZGCwHmxXXR5C2CIrgOPM35qsZV/LRrodWl+UflfzUQj004ov0PJTzoZeuZ6yZkOup00zlZtlXy/3cskZPf6o+ZlPnJFCr58HHcT/feqpQH09hyuv3BvcP/4BnHsuuDVrgAsuAN5+G+SOO0CKi1PyyCb3aHhA20dIQfP/J2fHCGdoy6olzzVnjMgzfT/R2FDT1wLbVZdHEKL6Z82ahaKiIgA7t1YSIn8AZigUkj1gLBgMyh4wJhz4JsDIAWOxWEz20LRIJCJ7uF04HJY93G50dFT20DRWJlYmVqb8LtO333I4/fRq1NcTXHddEEU2gpLbb0fpddeBIwTJ730Pw3/8I8jUqQVTponYTqxM48ukdABmIpHAwMAA1a46NnDKI2TrAMx4PA673Z5yfoWRGSepHT0zTvF4HMU7vo2q6UvzUDuYLhaLpZRHyU+lAzATiQSKiopS0ivJaeomm3I9bWpUbpZ9o9zLJWfkuK3kj5qf+cQZAWp1QGNndJTDOecQXHcdsO++O4TPPw/u9NOBQABkzz0R7+5G0aGH6uojpHKj3KPhAW0fQVtfanaMcIa2rFryfOKMltyM95NRzug5AJPFOOUhhE5G+MjJMpGPjo5q6tPYltrR44v0qgctfSEPNf308ij5qWRjZGREl9zs9shErqdNreaMldzLFWf0+qMkzyfO0NaBlryiAlizhkNdHYebbuIwOsqBO/lk/nLg/fYDt3Ur7PPng3v88Zxwj4YHZnJGy45ezpjV3vnEGbPa2yrO0IINnBgYGBgYMsInnwALFwI+H4D99gPeeAPkpJPAjY2B+9nPgPZ2/uwnBoYJADZwYmBgYGAwjKIi4N57gfnz+U8kAsDpBJ59FmO/+hWvtHIlcMopgN+fS1cZGEwBGzjtgrDZzGl2o3b0pKPR1dIx+tysesoGrPaVcYbu+a7KGY4Drr8euO8+oKQEGBkBUFSE6PLlII89BpSVAS++CBx6KPCf/2TFz2xzRkuHcSZ7eVjtJwsOzyNk6646BgYGBqvwySf8EU/PPw8ceOAOYX8/8OMfA599BlRX86dpnnJKDr1kYEgFu6uuwGHlWJYQgkgkknEeRu3oSUejq6Vj9LlZ9ZQNWO0r4wzdc8YZHjNnArfeCjQ2EvzrX1E+D4+HPyzz6KOBYBD44Q+B3/8e0Mi/UDijpcM4k708MuEMLdjAaReE9HyOXNjRk45GV0vH6HOz6ikbsNpXxhm654wzPBYt4pftVq2SvIymTuWvabngAn7AdNVVwJIlO9b1zPcz25zR0mGcyV4eVvvJBk4MDAwMDKbj+98Hbr89jPfeA/761x1Ch4OPJO/qAoqLge5u4Mgj+fU9BoYCAfXAKRgMWukHAwMDA8MExOgocPrpwF/+IhEuXQr87W/AbrsB774LzJsnGV0xMOQ3qAdOb775Ju677z68/fbbFrrDkA3Y7eZcUWjUjp50NLpaOkafm1VP2YDVvjLO0D1nnBmfx2GHAS+8ALS2As89J3n4ve/xcU9z5wLbt/MHQd1++7i4p0LhjJYO40z28rDaT9276j755BNs2LABtbW1aGhoYLu/TIQQ1e/3+1PqleNye3nsRJHnky9myfPJF7Pk+eSLWfJ88sUsuV4b//kPh9JSgr324s9+EvXHxkBaW8E9/DAAgJxxBrBqFbjy8rwvUyHI88kXs+RW2A4Gg3C5XFS76nQPy2bOnImZM2fC6/XivPPOw6RJk1BXVwe32w2Px4MZM2boNcmQhvr6enHELBwFT4j8Jb8jIyOylygODw8rXqIYiUTgcDjAcZyhSxTj8ThCoRCi0SgcDgeKiopQXV2NaDQqezHk2NiYeDEkIfzdUDU1NQiHw7IXQwplIoQgGo3C6XSitLRUtkx2ux3btm0TywOkXgwp2HA4HHC5XOPKJORbVFQkXn0glFWwOTY2plomYOdll1plom2nYDCY8o+u1k6JRALbt28X/ZW2k9xllzTtJC3T6OgoQqGQaD/b3BPacOrUqZplEnTLy8tRVVUlW6aysjL4fL6UaxakZYrFYiJnysvLx5WJECLWmbSdCCEoLi5GeXn5uNAGrf8nM9pJD/ccDge+/fbblPu8jHBPTx+x//7ViESiOP10DkceGcdZZ0V3lmnVKpADDkDZb34D7uGHkXzvPXDPPovRSZMMcU/gQU1NDRwOh2qZtPoIjuNQXV0tXlgr1Jce7in1e4JPLpfL0j7CDO5VVFTAv+MAU6EOzO73zHg/xWIx+Hw+0UZ6mZQu+ZX6qQXqGadgMIilS5di3rx5aGlpGTcie+utt9Db24s333wThx56KC677DJqJxh4ZOuS30AgAKfTmXJXj94RPSEkxQ6tL0I6l8slWwdSfWkeapdv+v3+lPIo+al0yW8wGER1dXVKeiU5Td1kU66nTY3KzbJvlHu55Iwct5X8UfMznzgjQK0OaOwY7SO2bAEaGoBLLwV+8Ys0O3/7G7B4Mbjt24EpU5BcswaB2bN1c4+GB7R9BG19qdkxwhnasmrJ84kzWnIz3k9G+ys9l/xSzzgtWLAADQ0NigOiOXPmYM6cObTmGFQg/SYslSnp6pVLP2r6WrbT7dD6okdfyz4hRLY8cn5q5alUB0bqJptyPW1qVG6WfaPcyxVnlOzoldP6ni05TR3Q2DHSR8yYAbzyCr/r7uST+XOfRDvHHw/09QH/+7/AW2+Ba2xEyQ03gLv8ct3c09M3mcEZLTtGOGOWPJ84oyU34/1kpL9Sei4H6uDwwcFBLFmyhNowAwMDAwODHPbYgz9MfMYM4OWX0x7usw/wj38Ap58OLpFA+RVXAGefDRTQeUcMExvUA6fFixdj48aNVvrCkCU4HI6c2tGTjkZXS8foc7PqKRuw2lfGGbrnjDP0edjt/FEFl1wCLF+etpmuvBx4+GGQm28GKSriA8e/9z1g82ZT8jaim6kO40z28rDaT1276hYvXowrr7wSs2fPttClXRfsrjoGBoZdDdu3AyecADQ28rewjFsx2RH3hG+/BSZNAh5/nFdmYDARlt1Vt3btWgwODmKzjlE/g37oGMsasj06OppxHkbt6ElHo6ulY/S5WfWUDVjtK+MM3XPGGWN5TJrE38QydaqCjcMOA5Ge93TiiUBHh+o9d9nmjJYO40z28sikn6GF7itXFi1aJB45wA7DLExIt43mwo6edDS6WjpGn5tVT9mA1b4yztA9Z5wxlofLBfzqV8BrrwFXXpk6JopGo8DeewOvvgqccw6QTAJXXAE0NwNpR4sYyZtWN1Mdxpns5WG1n9QDpyVLloybaXI6nVi9ejX+yo7KZ2BgYGDIEAcfzAeL//rXMhNKpaX8zcF/+AN/z92TTwKHHQZ89FEuXGXYhUE9cPL5fONkM2fOxNKlS8VDsRgYGBgYGIyiuhp48UXgjTf421fGgeOA88/nzzOYPp0fNB16KLBuXbZdZdiFQT1wam9vR21tLQD+oKhZs2ahqKgI++67L958803LHGQwHyUlJTm1oycdja6WjtHnZtVTNmC1r4wzdM8ZZzLPQxg8nX024PMBDoeMjSOO4M8zOO44IBTil+0uuwyQnP6cbc5o6TDOZC8Pq/3Utavuvvvuw3nnnYfVq1fD7XZj7ty5AJByyjWDcbBddQwMDAw7ce65/Ard3XfL7LYD+IHSlVcCN93E/33MMcCaNcC0aVn1k6Hwoef9q/uS3yeffBIcx+HUU0+lTvPUU0+BEIKhoSHU1dXh+OOP15NlQaK9vR2TJk0CAGzfvh0dHR2aaaSX/Fo1GCWEYGRkBBUVFbpOSjXLjp50NLpaOkafm1VP2YDVvjLO0D1nnDE/j1CI4MQTk5gzx4Y77+TkB08AH+/085/zweK77w6yZg1G5szJGme0dBhnspeHURvC1TemHkdw/vnn4+abbwbHcfB6vdTOAPyp44sWLcLSpUtRX1+P1atXyx5p8OSTT06Iow5aW1sBAG1tbWhra0NdXZ0oywfouczQCjt60tHoaukYfW5WPWUDVvvKOEP3nHHG3DwqKoAnnhjG5s3Ali0qiosWAW++CRxwAPDll8D8+Si69VbVIwv0+pmpDuNM9vKw2k/qgVNzczMGBgbEwcAJJ5yAm2++mWpHncfjwZIlS/DOO+/A6XRi6dKl8Pl849IuWrQIvb29+kuRR/D7/ejq6kq5nmbx4sXo6upiQfQMDAwMOlFZCTz3HLDnnsBjj6mMhfbbj48q/+lPwSUSKPvtb/k772Q2NjEwZALqgdOCBQuwatUqDAwMwOfzoaWlBQMDA2hpacG+++6rmbarqws33HAD9t13Xyxbtgwcx8HtduO+++5L0RUC0DOB1+tFfX09+vv7VfW6urqwcuVKdHV1ob293ZRBm3AtjfS2a+H3Qh8UMjAwMOQKY2PAXXcBbW0qg6fKSuCRR0DuvRfE4QD33HNAfT1/cTADg0mwG0nkdDqxaNEiLFq0SFeatWvXwuv1YuXKlTj++OMRCATgdrvR09MDt9sNn8+HWbNmGXEJAL9ENjQ0BLfbrTloam1tRX19Pdra2kRZc3MzAKChocGwD263GwBSZpeE3/UucVqFsrKynNrRk45GV0vH6HOz6ikbsNpXxhm654wz5uch2Kis5HfbLVwILFsGrFihkGDHkQWx2bNRfPrp4D75BDjySP58g/PPl40yN4MzWjqMM9nLw2o/dQeHA8BNN92EuXPnYv78+Rll/sknn6C/vx9DQ0Po6enBfffdZ8puMr/fj5qaGvT19cHj8Yx73t/fj+bmZgwODqbIhZkquTOr9KC5uRnz5s0TB2W9vb1obGxER0dHykAtHWxXHQMDA4M6/H7ghReAn/yEQtnn4881eO45/u+f/ARYtYo/84CBQQI9719DM04NDQ1YsGABAoEAGhoa0NjYiIaGBt2X/86cORMzZ84UbQ4NDWVlwNDZ2Sk7qyTMFvX29orP161bhzVr1mja7OjoENN3d3ejvb0dXV1dqK2tFeXCTy1YfRfQ8PAwqqqqMt4hZcSOnnQ0ulo6Rp+bVU/ZgNW+Ms7QPWecMT8PORsuFz/+eeop4IMPgKuvVknncoF75hngllv4a1oef5wPIl+zBtjxpdoMzmjpMM5kL49M+hlaGBo4zZkzBxs3bsSsWbPw5ptvghCCf//731i4cCHOO+88IybFAVQ20Nvbi6amJtlnwtKhMHBqampS1FWD9PgBYalObvYrF0gmkzm1oycdja6WjtHnZtVTNmC1r4wzdM8ZZ8zPQ8nGsccCy5cDDgfQ3q6SjuP4wzGPPJIfcQ0M8Ado3nIL8ItfUPuZqQ7jTPbysNpP3Zf8CnjqqafQ2dmJoaEhrF+/HmvXrh03aPrkk0+wbNmyvDtiwOv1oq6uTvZZbW2tZnyUFnp7e1PimYQZLNoZJwYGBgYGdUyaBPT0AA88ADzzDEWCI48E3noL+NGPgGgU+OUv+WMM2K47Bp0wNON03333YebMmZrB4TNnzkRLSwtWrFiBuXPnGp6NyjaGhoYySt/a2or29na0tLQAAG688UasXr2aOn0gEEj5m+M4EEJQUlKScpS8IE+Hllz4qOnT2hB0aH2RptHSl+ahZD/dppqfcjaUfKH1MddyPW2aidws+0a4l0vO6PFHzc984oycv0bsWNVHKNmR0586lb8U2OUCAgEO1dUaPKip4df47roLXFsb8PTTQH8/bJ2dII2NGdWXGveMcCaTOjPbhpy/Ztu34v0kpx+JRBCJRMbJ09+7ajA0cOrs7KS+n27mzJlYtWoVNmzYgJtvvhmXXXaZkSxNg9ZZSi6XK+Pdb52dnfB6vejq6sLg4CC6u7t1zTbts88+svL29nZcccUVAPi7eMrKyjAyMpJy2FdZWRlKSkowPDycMl1ZUVGB4uJiBINBJBIJBAIBcByHqqoq2Gy2caRxOp1IJpMYHh4WZRzHwel0Ih6PIxQKiXaKiopQXV2NaDSKcDgs6tvtdlRWVmJsbEwkKiEERUVFAIBwOIxoNCrqp5eJEIJEIoFoNIrS0lLZMtntdiSTSbE8AFLKJNgQToVNL5NgJx6PY3R0NKWsFRUViEajGBsbUy0TADgcDpSXl2uWSU87Sf/xtdpJ2qbSdhoZGRF1bTYbdTull0lqP9vcI4SI9rTKJLT3yMgIqqqqZMtUVlaGoqKiFM5IyxSLxUQ/y8vLx5WJECJ+gZG2EyEEZWVlIIQgGAyqlgmAJe1Eyz2HwwFCSEodGOWe2X1Eepm0uFdSAhBShhNPLMHPfjaG00+PiDyIx+NwOBzjy3TRRbB973tILl6Mok8+QdXJJ2OsvR2l116LJMeNK1N1dTVKSkpS6ksP95T6PUII7Hb+NWx1H5Ep94T+VloHVvR7ZryfpDbSy7RixQqqmzzUYGhX3cKFC7F+/XrdmT355JOoqamx/MoVtV11wrPOzk5xRkiK5uZmeL1e9OXg3A8hqv/TTz9NCZI3e8aJRp5v35BZmViZClmeT76YJc8nXwBg82YO8+cT3HAD8NOfUtoJBIALLgD3+OO88JhjQB5+GNhrr7wo00Rsp1yXSWnGKRgMYu+99zb3yhUpjG7XX7RokaEBVzaRD6d7V1dXw+l0ih/h79LSUnFGQRjxS/+mkQMQvw2r6WvZTrdD64uQTpjOVtOX5qFkX/rtXstPpXIK32Zo5Ebr3Sq5nja1mjNWcS+XnNHjj5qf+cQZmjqg5Z6arlY9ms29mTOBnh4OnZ0cIhFtHgAA53QCjzyC0VWrQCorgb//Hdx3vwtux52seupLjXtGOGNWe+cTZ8xq70xslJaWjnu3Cr/TwtDAycAklYglS5bg5ptvNpw+UwineKsNkMw4vTwTSNdopcsAZsmTyaSmPo1tqR09vkiXPLT0hTzU9NPLo+Snkg29H7PbIxO5nja1mjNWci9XnNHrj5I8nzhDWwdG2klvPZrNvVmzCF5+mSAaJXj9dfn4o3F2AEQWLwbp7wcOPZQ/KKq5GeS880BCIV31pcY9vZwxq73ziTNmtbeZnEn/mwaGz3G67777DAV7z5kzB1dccUVOY508Hg+2b98u+8zr9eb8Qt65c+eK697S0f4555yDpUuXAshsDTkUCom2M4lfEOzojXESfqeJcQqFQigrK1ONcRLW6IVvFukxToKfSjFOAGRjnAAUTIyTtE2l7WRW7IzUfra5Rwh/27nL5aKKcQqFQrDb7aoxTpFIRDXGSfBTKcYpkUgAGB/jJHTYhRDjJPBdqAMrYpyM9BHpZTLCvQ8+4HDmmZV47LEEGhpsqmUS+4jJk+F69VWQ3/4WXEcHuPvvR+KVVxDu7ETF/PkpcTNy7aTGPaV+jxAi/l4IMU7RaNTyGCez3k+CLL1Mq1evxv333y+WVxg06bkY2FCM01tvvYXFixdj06ZNepMCAE444QS89NJLhtLSQOvk8Pb2dni9XnR3d497xnGcYjqrIcQ4+Xw+OJ3OFJ/kmsmIXAikdjqdKVOY6fpatgkhKXZofRHSSe/yU9KX5iEQXA5+vz+lPEp+ytkQXnLV1dUp6ZXkNHWTTbmeNjUqN8u+Ue7lkjNy3FbyR83PfOKMALU6oLFjVR8hlRvlHiEEL70Uws9/XoVnn+Vw6KF0foqcefll4IwzwH3+OUhREfCb38B/4YVwTp6sWF9q3DPCGdqyasnziTNacjPeT0Y5EwgEUFNTY12M05w5czBz5kxceeWVRpJnvN0/UyxZskT2wt3+/n64XK68OajSKlRVVeXUjp50NLpaOkafm1VP2YDVvjLO0D1nnDE/D6M2GhvL8fDDBF99ZSCP444D3n0X5Cc/AZdIgLv2WrhOOQX4+GPDvjLOZC8Pq/00tFQHAKtWrcKsWbNQV1eHc889V1daA5NcuqA1MPN4PFi8eDG6urpSdta1t7fLzkJlG/PmzRO3rpq9VCfI5Za1pKCZ3iWEyE6FClCahi8uLkZFRQVGR0c1p3cJIbLLJdIyCdOyAtLLJPipVKbq6mrZKeuqqipLl0vMmoYnhL9iQPrtyswlIEEut6yVLe4B/FJrLBbTLBMhBMXFxaplikQiiMViiu0kcEapTOXl5eLWbGk7VVZWAhh/Jky+LdWVlJSI/8Pp7WTmEpDRPsIM7hFCcNxxlSgutuG668Zwwgkx/M//JGXLJNtH2GzAvffCsWAByi+7DNybbwIeD8K/+x2i554Lm8zyoxr31Po9u92OUCiU10t1lZWViMVilvZ7ZryfEokEQqGQ4pKq0lKdsPxOA0NLdQK6urpw/vnni4MQmqj0m2++GQMDA1i1apXRbBXR3t4Ov98vntztdrvFq1M6OzvH6a9cuRIA3yEPDg6Kd+7lCmypji3V6ZWzpTq2VGdEDkz8pTopD554gqC9Hfjb3wDhdi/aPgIAsHUr4mecgeJXXuHtn3ACsHo1uL32Ykt1efZ+ysZSneEZJwBoaWmB3+/HFVdcgd7eXixbtgxLly5VzPTJJ5/EjTfeaNkZSXoPtWpra7PEj0whdDLpMiVdvXLpR01fy3a6HVpf9Ohr2Re+KarVGa2PSnVgpG6yKdfTpkblZtk3yr1ccUbJjl45re/ZktPUAY0dq/oIqdwo96RpTjuNQzgMnHAC8MYbgLBxmraPIHvuiZGnnoLzkUfAtbeDe+kl4OCDgdtuA3f22QCFHSOcMUueT5zRkpvxfjLCGaXncsho4ATwgw+Xy4Xzzz8fbW1taGtrQ0NDAzweD+bNmweA36nW09OD3t5edHR0YMaMGZlmy8DAwMDAQI2f/xyoq+NvXSFEHOvQw2bj77dbuBA4+2x+BHbOOcC6dUBXFzB9uhVuM+QhMh44AfzMU0NDA1pbW7FhwwZxkCRAmArr7OwUY3QYlCFMHwvgOHOn7Wns09oQdGh9kaahmW4XPkr2022q+SlnQ8kXWh9zLdfTppnIzbJvhHu55Iwef9T8zCfOyPlrxI5VfYSSHT0+KvHg6KOBTZuA1lbguec4VFZq9xECRD/22w949VXg1lvBXXMN8PzzIAceCHLbbSA/+pGsHSOcyaTOzLYh56/Z9q14P+nxRanMcsgoxkkOgUAAnZ2d2LhxI/x+P1wuF+bNm4empibMFBaXGWQhxDjNmjWLOjhcKaAwGAzKBt8FAoGMg++EAF2hgxCC7yKRiGxAYTgc1hUcLi0TITuDw5XKlH6YaTaCw9PLJARJ0pSJtp2k/5pa7RQMBlOmnKXtJC2TnnaSlkkuQDeb3APUg8OlZSJkZ4CuUpmE85qU2kngjFKZysvL4XA4ZIPDhXvwtMpkRTvp4V6670a5Z0UfYQb3CCGorKwcVyZCgOXLq/H++xwefTSA0lLtPsLpdCIajaac9Waz2VC9dSuSZ50F28aNAIDowoWI3XknKvbbL6N+T6lMVrVTSpmywD2r309yweFCmdSCwwcGBqhinEwfODEYR7aCw5PJJGw2W8qL1siIXmpHz4xTMpkUB4Zq+tI81AJ9E4lESnmU/FQKDhc6TWl6JTlN3WRTrqdNjcrNsm+Ue7nkjBy3lfxR8zOfOCNArQ5o7FjVR0jlRrmnxYNkEjjvPA6TJxN0dGTImVgMuPlmYPlycNEoSFUVuJUrQZYu5Zf3YIwztGXVkucTZ7TkZryfjHLG8nOcGKyF0MlIX9rpskzkwmhcTZ/GttSOHl/ST7lW05d+c1DSTy+Pkp9KNoSt/LRys9sjE7meNrWaM1ZyL1ec0euPkjyfOENbB0baSW89Wsk9NR4UFXFYvRq4+moOQ0McgsEMOFNcDG7ZMqCvD/F588AND/MXBx9/PLiPPzbMGbPaO584Y1Z7W8UZWrCBEwMDAwPDLge7HaiuBq64AmhvL0PGay8HHojQCy+A3H47UFEB/P3vwHe/C/z+94BkWZih8MEGTnkIaXCbdOo5m/Jc5MnKxMqUb76wMk38Mt16K0F/fxF+9zsTbNtsIL/8JfD++/xZT5EIcNVVwKGHinFQrJ3yt0y0MGVXHYO5sPqSX+mluJlcoijY0XvJrxBESHPJ78jISMaX/ArPlS755ThO9pJfjuMK5pJfaR1I28msE6ml9rPNPaENaS/5HRkZ0bzkN/2y0vRLfoU8lC75FX5Pv+RX+BTCJb+CbaEOrLjk10gfkV4mI9wTeFBRUQGHw6F5ye8DD4TR0jIF55/PYcoU+Q0kwsGKQn2pcm/GDIw99RSSDz2EsiuvhO2dd1C1cCFwwQUI/+Y3iJaVAeA5I2xSKIRLfmOxmOWX/JrxfpLayJtLfhmsgRAc7vf7U4LThIFTOqyU5yJPq+X55ItZ8nzyxSx5PvliljyffDFLnk++mCEnBEgmObzxBsGRR5poe9s24PLLwT38MJ/P7rsDt90GNDcDkjgb1k65LVMwGITL5WK76goN0oGTdFedmRBG1na7XVcwnFl29KSj0dXSMfrcrHrKBqz2lXGG7jnjjPl5ZJszn37K4cgjgfvuA77/fX321HQIIUisX4+iX/4S3KZNvPDEE0HuvhvxvfdmnDExD6M2hKtv2K46BllIp2ZzYUdPOhpdLR2jz82qp2zAal8ZZ+ieM86Yn0c2ObPPPsBzz/EHgr/2mn57ajqhww4D3nkHWL4ccDiAF18EDjoI8WuvBSThAPmMic4ZWrCBUx5CGishFz+RDXku8mRlYmXKN19YmXa9Mnk8BI89RvDeexbkWVoKcs01IO++C3L88eDGxlB2ww3AIYeAPP88a6ccl4kWLDg8D2F1cLj03JBMAj8FO3qDw4XfaYLDQ6FQxsHhgp9KweEAZIPDARRMcHj6WTBmBx1L7Webe4ToCw4PhUKaweGRSEQ1OFzwUyk4PJFIABgfHC4EJhdCcLjAdyuDw430EellMsI9gQe0weHpfcScOcOYMwfo6irGkUcmcNBBVUgkEprB4UrcG9fv7bYbsG4d7E8+ibLf/AZFAwPAyScjevLJGPv975Hce++8DA5P31RhRb9n1vtJkLHg8AmObJ0cPjw8jKqqKpH8RoLvCCEpdugDMPl0SmvIUn1pHmon+gaDwZTyKPmpdHJ4KBRCZWXluHgVOTlN3WRTrqdN9cgJAUIhwO/nsH17El99NYpEohyhEIdQCAiFOASDBCMj/CpDOCz85BAOE0Qi/NE10ajwk0M0ShCNJkGIDYkEf3pzIsEhkSAQ+lDeBblyEthsHDiOP4xZ+FlUBBQVcbDZCIqKeJnNlkRxsQ12Owe7nZfb7akfII7SUjuKi/m/i4uB4mJe3+EACImiosIBhwMoKdkpdziA4mKCZHIMTmcZSkogyh0Ognh8FLW1FSgpQcqntJRDSQnZkXd2uKEmB9T/b2jsWNVHSOVGuU3Td9D0EatXA7feyl9N53Co15eaHaV+jxCC4c8/R/XttwN33AEukQApLQWuvBK47DJwZWWm9BH5whktuRnvJ6Oc0XNyOBs45RGEgRNNwzEw6AEhQCAAfP458MUX/EafbduAb75J/X37dsDnA4aHgfJy/ib5mhrA6QSqqoDKSv6n8CkvB8rKUj+lpfyHH2Sk/rTbhcEOxIGO8FPoi9N/EpL6SSZ3/tw5ANv5M5EA4vGdP4VPIsEP4uLxnT+F34WPMNBLl+0cAEIcFEYiO59FIjv/Fn5P/4yN8fkVF++sI+Ej1Jv0p1C30p/ST0VF6u8VFXz7CL8XF+eGaxMN118P/PnPwIYNfP1ahg8+AC66CHj5Zf5vtxu45RbgRz/a+c/AYBn0vH/ZwCmPkK1dddFoFA6HI+MdUkbs6ElHo6ulY/S5WfWUDRBCEAhE8dVXDmzZwmHzZmDzZuDTT4GtW/mB0pdf8i/fPfYAdt8dmDYNmDoVmDKF/yn8XlvLD5Sqq8VrthhnKJ/T5c0PnsbGUj87Z+x2/j46uvPvcJj/W/iMjKT+Lv2EQvzPRGLnYEoY8O78SVBRkYTLZYPTyaG6mm/zqip+kOx0Ai7Xzt/TZ8poYAZv8oEzhAAPPkiwaFEUlZXGeEXNGUKAJ54ALruM/8cFgMZG4PbbgQMOoC6/FchGn5hLzujZVcdinHZBhMNhOByOnNnRk45GV0vH6HOz6skshELAwACwaVPqZ2AAGBlxYJ99+F1BM2bwn9mzgT333DlYKikxnjfjDN1zrXQcJywL8oMUK5FM8gOrUIifQZT+DASAbdsiiMXKEAoBg4NAMMh/AgH+4/fv/N1m4wdSLtfOWUjhU1vLfyZNGv87IZnzJtec4TjgrLOAjRsjeOghB+64Y+eXCj12qDjDccBPfgL84AfAjTfylwf39ACHHMLPRi1fzjdCjpCNPtGMPKz2kw2cGBjyDH4/8OGHOz8ffMD/HB4GZs0C9t2X/xx/PNDays/oFxcH4HI52Yw+gwibbeds07Rpqc/4pdsonM4yKs5EIvwAyufj+enzpX68XmDjRn6pd2iI/7l9OxAKOVFTw89mTp7Mf4QZzqlTgd122/n7tGkWL4VliL33TuLDD/nJoFtusXj1rLISuOEG4NxzgV//GnjmGeCOO4BHH+XXDs8919g0IIMpYDWfh0jfGml2oCiNfVobgo6e4HBpMKWavjQPtYDF9PIo+SlnQ8kXWh8zkSeTHDZtInjnHf54l3ffBd5+mw+wPvBAgv3352fnv/99/uf06RwAueBwgkAg95wxy04+cUaPP2p+5tOGAjl/tfRLSjhMmUIwZQp9noQQDA0FkEw68e23wPbtHLZtI2Jc3bvvCjF2HL7+muDrr/kYsWnT+BnS3XfnsNtuQE1NCerqiDh7uueeHCoqMucBbR8hwOEgePJJguOPB+69F7jwQvo8jXAGADi3G+Spp4DeXuDSS8F9+CFw/vnAPfeA3HIL0NCgWAe07UQrT/fXbPtWvJ/0+KJUZjmwgVMewurjCKR3VGWy3VOwo/c4AsFfmuMIwuGw5nEEwnEBwnp2+lZjwSel4wjsdrvscQR2u9204wiCwRF8+GESb79dhLfesuPdd4vx4Yccpk9P4sADEzjooATOOiuBe+4pwd57y281BpTbSdqm0nYya5u71H62uUcIEduA5jiCcDiseRxBPB5XPY5A8FPpOAIhXfpxBMJuqkI4jiD9SAYrjiOIRsMoKwP23LMIBxxQjUhEvkzhMF+msTF+IDU0VAKfrxRbtkSxeXMc69fb8NVXNnz1VRG++IJDMglMn57EHnskseeeScyYUQS3247a2lHssUcC06cnkEiEqY8jUOsjOI6/q47vhwN44gkbYjEgGOTgdNJxT6nfI2TnvYeq/0/z5gEvvwzHn/6Eso4OcO+9B27hQsROOAHh665Dct99s3IcQfp9fVYcR2DG+0lqgx1HMMHB7qqbOGXauhV4/XXgX/8CNm7k8M47BHvvDcydC9TX8z+/+10OpaWFU6ZsyfPJF7Pk+eSLWfJc+hIKEXz+OfDZZ/wmiE8/BT77jMOnnxJs2cLvHp00CZg5E5g5k8OMGQQzZwJ1dfxy99SpgM2WWZmamoDzz+fQ0JDlOhgaAnfddSD33AMuHgex2/npr2uuATdpUl61UyFxj91VV6DI1q66sbExlJaWit8asmlHTzoaXS0do8/1+BmNAv39/BUNr7/OfwDgiCP4z2GH8YHaFRWqZgzDrDa12j7jTP4gG76akYdRG4kEwebNY/jyy1Js3szvNPV6+SD4gQE+XlAYRNXVEcycGcNBBxVj//051NbS+fHvfwM//CHw/POAx6PuqyWc+e9/+YCrv/yF/7umBrjmGn4QZUFg9ETnDNtVx6CKSCSC0tLSnNnRk45GV0vH6HMleSjEzyS9+irw97/zg6b99weOPBJYvJi/+HzPPVVdNh1mtanV9ndVzuQjsuGrGXkYsWGzAZMmReB2l+Koo8Y/HxnhB1IDA8DHHwP//CfBI48A//kPwHHAfvsB3/kO/3990EF8jGFFRaofhx4K/PGP/P/8Bx/wYxU1X03nzH778QdM9fYCv/oV8P77/M977gE6OoD//V++MCZiInNGD9jAiYFBAyMjwD/+Afztb/znv//ll9qOPhq4+mp+Rsmq2SQGBgbzUVEBHHww/yEECATCcDod4Dh+N+B//8sPoj78kD8N4P33geHhahx0EHDggXy62bOBo47iZ5gdDr6fyAkaGvhdJX/8I98hDQwAixbxHdQttwDz5uXIsYkLNnBiYEhDJAL885/ACy+U4vXX+U7z0EOB+fP5HcH19exUZgaGiYpJk/jZ4yOP3CkjBPj00yA+/9yJDz4A3nsPWLuW3xk4aRK/E3DLFuD66+044QR+V2BWUVQELF0KnHYasHIlP2B69VW+4/rpT4Hf/54/5I3BFLAYpzxCtmKchJ1qmcarGLGjJx2NrpYOzfPR0TAGBsrQ28uhp4dfhps9m+Doo2NYuLAYhx3GWREyYBrMalOr7U8kzsg9t7odzEQ2fDUjj3znDCH8Kf1vvQX8/vcEn35KUFTEoaiIw9y5ED/19QQVFVnkzNatwFVXAQ89xP9dUgJceimwbBl/FLwBTHTO6IlxYgOnPAK7qy572LYNWL8eePFF/g6qSZP4mw0aG4Fjjsnvg/gYGBjyD4kEcOqpwM9+Bnzve/yBoMLnzTf52O0jjgAOP5z/edBBWTjDsr+fDyD/29/4vydNAn77W/4sKDZtngJ2V12Bgs046delnR1wOMrw739zePFF4IUX+G+JDQ3ACScQHHXUGOrqxu92YbMH5tovNM6wGafc51EonBF0gsEwqqvLsGULhxkzdj5LJgnef38Mb71Vin/9i8Prr/N9UH09wRFHxHD88cU44ggOZWXURaQHIfzWv8svBz76iJftuy8fQP7jH1MHkE90zuiZcVK4cYchlyBk/OnF6bJM5JFIRFOfxrbUjh5fpIfCaekLeajpp5dHkG/fTvDoowRnnmnHtGn8TDVAcOedBF9+SfDYYwRnnkkwefL49ITwF0WaWe9WyPW0qdWcsZJ72eJMug29/ijJ84kztHVgpJ301qOV3KPhgZmcIYQgmYwgFiM4+WT+YuCdegR77TWGM88E7r2X4K23+GW99naCaDSB3/0OmDaN4KijCJYtI3jxRYJQyKQ+guNATjoJ5J13QP7wB5CpU/lLLk89FeSYY0DeeCNrnDGrva3iDC1YcHgewuqTw0OhkGg7k1OBBTt6Tw4Xfqc5OTwUCmmeHL7z5FsO//mPDS+/XInnnwc++IDg2GPjOPbYYdx4YxFmzdpZph2ui5A7ORyAaSeHG2kn6T+yVjtJ21TaTmadSC21n23uEUIwMjICl8tFdXJ4KBTSPDk8/dTs9JPDBT+VTg5PJBIAxp8cnkwmQUhhnBwu8F2oAytODjfSR6SXyQj3BB7Qnhwu5KF2cngikUjhjBr3nn22CsccAzidIzjuuLhivzdvHsF3vxvF1KmlGBoawb//TfDPf9rR0WHHaafZMXs2cPTRERxzTAxz5iRgt2fQTiMjfPD4ySej9M47UXrPPeD+8Q/g8MMRPfVUjF1zDTBzpurJ4dFo1PKTw816PwkydnL4BIewVOfz+VKW6sw8JVU4Mt/pdIrkN3IyKyEkxQ6tL0I6l8IN31J9aR4CwdMRjwMvvhjChg2VeO45PmbglFOAH/yAw/e+R2C3q9sQXnLClQpacpq6yaZcT5salZtl3yj3zOYMAHE5XNq2tNxW8kfNz3zijAC1OqCxY1UfIZUb5R4ND+T8NMoZuTz/9S+C224DnnjCGGfGxoDXX+fQ00OwYQN/ysAxxwAnnACcfDKHvfbKsJ0+/xzk6quBBx8ERwiIwwH88pfgrroKROF/LVPOaMnNeD8Z5UwgEEBNTQ2LcSo0sJPD6XSHh4GXXgKefRZ48UWCWbOS+PGPbfjxjznst5++/JSem1VP2YDVvk4EzujRYZzJnzwKhTNqOps3AxxHMHVqZpwZGuJjvF94gQ9ZmjQJOPlk4KST+KMTDAeav/02H0C+YQP/d20tfwL5BReknEA+0TnDdtUVKNiuOmV8+y3w3HPA00/zh1EefTTwox/xs0u77ZZr7xgYGBjk8cc/ArfeyvdbChM5ukEIP955/nng//6PP/184UL+sPDvf9/ArmBC+C3Gl1/OH4MO8AHkK1fyHW2efxEwA3revyw4PA9h5VhWWIfPNA+jdvSk+/RTgptuimD+fILvfAf461+Bs8/mjyh57jng3HOBqVPV7Wnlp/TcrHrKBqz2tZA4Q6NrlBNazxlnzM+jUDijpnPOOcBJJxH8+MdxRCLmcIbjgDlz+KOaXnuNP+V84ULgkUeAvfbiv1Defz/wzTc6DH7/+/xorLOTvwV50yZ+JHbcccDGjbsEZ2jBBk67IPQEwVlhRy2d1wvcdBN/jclhh/FXHlxxBfDFF3ynsGjR+OtNtPww+tysesoGrPY1nzljRJdxJju+mpFHoXBGTefGG4H99ktg61ZjdrUweTL/pfLZZ/kvlj//Ob+s953vAMcfD3R18bP2mrDbgZYWPqDqqquA0lL+Qs5584Azz0Ri8+aMfdVCLjlDCzZwYsg5Pv6YvxHA4+GDH7duBW6+GfjsM+C228JYuNCSy74ZGBgYsgKbDejoCGP33YF166zNq6KC/4L5yCPAl1/y4Uv//Cd/J/BJJwEPPgikbVIbj6oq4Prr+c75Zz8DAHCPPILqQw8Fli/P4cV8+QE2cGLICT7+mP+//O53+SnmoSHg3nuBTz/l74M7+mj++iUGBgaGiYLRUeDKK3fehGI1HI6dg6XPP+evs3v+eWDmTH5w9dxzQCymYmCvvYCHHwbefBPk6KPBhcPgfvc7fhT28MOA5EiBXQls4LQLosyk42n12vn4Y+CGG4BjjqnGCSfw33pWrwY++YSfYTr8cP6bmd48tHSMPjernrIBq33NFWeMpGOcoUM2fDUjj0LhjJZOWVkZJk3iBy7LlvExm3r9zASlpXzI0po1/BfUH/8YuPNOfmz0q1/x4U2KmDsXePllxB57DGTGDH4UduaZfKf92mum+plLztCC7arLI0zEXXVeL3+L+Jo1/KzS4sVAczO/ZL4LbNRgYGBgGIfXXwe2bOHPo8w1Pv2Unzx68EGgrIyPjzrzTP5UAlmMjfHLAjfcwJ8NAwA/+Qm/A2/PPbPmt9lgu+oKHFbvWkg/cdZsO599BtxyC3DoofyS21dfAX/4A3+eycqVBN/5ThAA3W4XLV+1dIw+N6uesgGrfc0GZ8xKxzhDh2z4akYehcIZLZ30Z0ccwQ+a/vAHgk2bhnPKmb335uPA//tfYNUq4P33gf/5H+CMM/jYKME1sQwlJUB7O7/r7rzz+G/Ajz/OL99dfz0/sDKIXHOGFmzglIcgxNq7gBKJhKY+jW2pnW3bgLvv5u9amjePYHCQ4KabyI4Ab4LDDiMAdqaj9V3IQ00/vTxKfirZEK5ooJWb3R6ZyPW0qdWcsZJ7ueKMXn+U5PnEGdo6MNJOeuvRSu7R8MBMzmjZkePG8DDBkiVlGB21vo/QkgMEhx9OcN99wMAA//sFFxAcfDB/v6ffn1YHU6eCdHXxRxV873t8ANfVV4Psvz/IU08BMnma1d5WcYYW7K66PESh3FX3+ecjePRRB556qgTvvmvHKackcdllYRx1VBx2O3+/kc1WiXA4O3fVCafE6rmHSgC7q47dVcfuqmN31Zl1V10695T6vfPOI3j//RKccUYZHnhgBMmkdX2EHu6VlUVx+ulh/PSnwJtvFuHBB0vxu9/ZcfLJJbjggmF85ztkZzt5PAi/9BLI44+j7Le/hW3zZj7yfMECjK5YgdisWbrKxO6qY9CFQrirLhzm8Oc/Ezz2GPDqqwQnnMDhtNOAE0/k4HBo+0JI7u4dY3fVsbvqBLC76thddVbfVUfLmW+/DWDdOhdaWsbvJDarjzCDM99+C9x11xj+9KdSHHggcPHF/L15RUUS/ZERYMUK4OabwUUiIHY7cMkl/BUuVVWa+bK76hh0I1t31cXjcdjtdup7fGIxoKeHX8Z+/nnge98DfvITgu9/Pw6nk96O3vxpdLV0jD43Uk+5gtW+mmXfqB3GGfORDV/NyKNQOKOlQ8OZf/yDw6ZN/Enj+QjBV0LsePppDrffDvh8/ADq7LOB8nKJ8ief8Fv1nn2W/3v6dD7wdckS1V1BueQMu6uuQJFPu+qSSX6X6WOPAU8+CRx4IL9xYtEild0WDAwMDAyGsGUL/6X0/vv5mZxCwL//zd/D98orwC9+wX9qaiQKzz/Pj6wGB/m/588H7rqLf6HkGdiuugKHlWNZYcpYKY933+WvOJk5E7j0UmDWLKC/nz9zZOnSnYMmLTtG89erq6Vj9LnR8uUCVvtqln3GmfxBNnw1I49C4YyWDg1n9tmHv8T87LP5nW35BrkyHHoo8MQT/AXGn3/O78b79a+x82qZk07iC/O73/EHSf3tb8Ds2UBbm+zp47nmDC3YwGkXRDpBtmzh71I6+GB+RqmkBFi/Hti4Efh//w/YYw86O0bzz1RXS8fo80J4AQqw2lez7DPO5A+y4asZeRQKZ7R0aDgzbx4/y7/77ppZ5QRKZair44+cee89/v3h8fCXsA8MgB8wXX01f/Hoj34ExOP8haQHHgj85S/UeZjhp1lgA6ddFNu380Q/6ij+8Nevv+aniD/+GLj2Wv5IDgYGBgaG7GL+fP6u3bPO4nf4FxKmTePvHd20Cdh3X/798vOf71ipmzkTeOYZ4M9/5g+P2rIF+MEP+G/rcrcf5zHYwGkXwugoP6162mkV2G8/4I03gN/+lj+w8vbb+WnXPI9pZWBgYJjwqK7m75k766zCvA7O6eRDPjZt4pfvvvc9fgbK6wVwyin87NNll/HbCJ96Cth/f/7+lx3HfOQ7WHB4HsGKXXXxOLBhA/Doo8D//R9w9NEEp51G8MMfcigvz2yHVDKZhM1m073bhTYdja6WjtHnRsuXC1jtq1n2GWfyB9nw1Yw8CoUzWjpGOBOLASeeyN++sHw5RaEtRibtOTzMx4TfeSc/brrmGn7SCe+8A5x/PvCvf/F5zJ2LZGcnbHPmZJ0zenbVsRmnPESmHRkh/G6HSy7hL3D8/e/5KdNNm4BnnuGwZElmgybBRyP/QHrS0ehq6Rh9brR8uYDVvppln3Emf5ANX83Io1A4o6VjhDPFxcC6dfzFvPmATNqzqgq48ko+FGSfffg7gy+7DNi+53f5e11WrQKcTnAbN8J22GHgfvMbw1e3ZMIZWrCBUx7C6CTgpk38N5P99uOnRadP55fjXnkFaGnhd8SxHVJ0z9kOKfPtM87kD9iuOn262dhVJ4eaGuC73wV++Uvg73/XLIqlMKM9q6t3xoknEsABBwA3dtgwekYr8NFHIKeeCi4e57/tz57Nb9fLkp9sV90uhK+/5i+qPvRQYMECIBzmv6W89x5/D+Pee+faQwYGBgaGTNDczF8KPDCQa0/MweTJwG238V/sP/qIj4PqfG53xJ9Yh5GHHgKZNo2/dfjoo/nDodKuMMo12MDJIPx+P7q6utDY2Cj7fOXKlVi3bh26urrQ3t5uat7BIPDQQ/whaQcdBHzwAb+7c/NmoKMDOOQQU7NjYGBgYMghjjmGn4T5wQ8AyZV+BY8ZM/h32fPP84eMz54NvFTxv/xL7dxzeaV77+VfdD09uXQ1BeySXwPo7e2F1+vF4OAgvF7vuOcrV64EADQ1NYn6zc3N6O7uNpxnJAK8+CJ/xkdPD9DQAFx4IR88WFJi2CwDAwMDQwHg7LP5gUVJCR/HWgChdNQ45BB+8PTii8Cll5bhT/tW4JZb7sP//PSn/MnLXi+wcCHQ2srPEuy49y5XYLvqMkBXVxc6OjowKBwnvwM1NTXo6+uD2+0WZRzHwefzKV5UCow/8j2Z5Ne1H3uMP1F2zhzg9NP5YMFMbmQhhJgSFGrUjp50NLpaOkafm1VP2YDVvjLO0D1nnDE/j0LhjJaOWZx55hn+Joc776ROYhqywZlolKCzk8MNNwA//Slwza9H4OpYxm/LA/jo8vvv52NTTPRzwl+54vV6UV9fj/7+flW9rq4urFy5Ulwu6+3ttdy3/v5++P1+1KZd6OZyubB27VoqG2+/TXDZZXx80hVX8LOU773Hn+Z91lmZD5qSyaQpgb5G7OhJR6OrpWP0uVn1lA1Y7SvjDN1zxhnz8ygUzmjpmMmZhgbg1VeBe+6hTmIKssWZoqIkLrqI4IMP+CMZ9p9bgdUH34nkhr/xh2hu2cJXwi9+AYRCpvk5YYPDW1tb0dzcjM7OTs1BU2trKwCgra0NLS0t6OjoQGdnp+WDp6GhIQAYN7NUW1sLv99PZePcc/kDxF5+mT/e4uKL+RNZzcLw8HBO7ehJR6OrpWP0uVn1lA1Y7SvjDN1zxhnz8ygUzmjpmMWZykrgueeAFSuAN9/UlTRjZJMzkybxk0y9vfyqy5FXHod3Hn6Xj1EB+NinQw6R3XlntZ8FFePU2dkJgA/MFuKI5NDf34/e3l5RX0BHRwfq6+vh8/ks81FtcLR9+3YqG2++Cais6DEwMDAw7MLYay/+mJl99pl48U7pOPBAfmnysceAkxZXorn5Htzw7KmouPhc4JNPgGOP5beQL1/OH7eeBRTUwIkWnZ2daGhoGCcXYo56e3vF5+vWrcOaNWs0bXZ0dKTELClBLYaJFhP5n4CBgYGBIXO43cA33wAnnwy88AI/QzNRwXF8fO/JJ/PnQH3nFwtw++/examvXAzuwQf5W+pfegl45BHgO9+x3J8JOXDq7e0Vd7Slw+12o6enRxw4NTU1KeoagRDb5Pf7UwZRQ0NDqKurMy2fTGBWcJ9RO3rS0Z78a8XzQgnyBaz3lXGG7jnjjPl5FApntHSs4MyUKfzO6qYmftxg9YRLrjnjcvHLd2efDVxwQTW6ah7Aw/f+AFN/0wL09wMeD9DRAe6MMyz1saBinGjh9XoVBym1tbWa8VGZwOPxwOVyibFOAvx+v+wsmBysvgbB6XSacn2GETt60tHoaukYfW5WPWUDVvvKOEP3nHHG/DwKhTNaOlZyZvly/kDJ664zbIIK+cSZ+nrg9df5e+8OXr4IXRe/D3LCicDYGLhLLoHztNPAffWV7rxpMSEHTlpIH9SYjZaWFqxbt078W1gapFnqA4xfuUJrOxaLmbJDyogdPelodLV0jD43q56yAat9ZZyhe844Y34ehcIZLR0rOWOzAQ88AFx6KX+NiVXIN84UFfFX0bzxBvDka7vjCN/z+PKqu0FKS4H160GEw6F05E2LCTdw0tq55nK5qHe3KaG/vx8rV65EZ2cnvF4v2tvb0dXVJT7v6OgAwB+H0NXVhZ6eHl2HXwYCgZRPMBhEIBDA2NgYCCHiB0DK37TyUCikqU9jW2pHjy+hHVtIafSFPNT008uj5KeSjZGREV1yo/VuhVxPm1rNGSu5lyvO6PVHSZ5PnKGtAyPtpLcereQeDQ/M5IyWHb2c0VPW8nKC2lqCU04BenvzlzNmtbdUts8+BC+8QHDBhYDn/l/gnnP7EDvwIHDffgucfDLIpZcCkYioPzY2Nu7dKnxoMSFjnLSQaQC3x+OBx+NBW1uboo7aMy3ss88+svL29nZcccUVAICSkhKUlZVhZGQE8Xhc1CkrK0NJSQmGh4eRTCZFeUVFBYqLixEMBsV/VI7jUFVVBZvNNo40TqcTyWQyZVunMI0aj8cRCoVEO0VFRaiurkY0GkU4HBb17XY7KisrMTY2hsiOewIIIeLv4XAY0WhU1E8vk/BPVFZWhtLSUtky2e12jIyMiP4BSCmTtGNyuVzjyiQgHo9jdHQ0pawAEI1GMSa5pVuuTADgcDhQXl6uWSY97SR0EOllkmsnaZtK20moGwCw2WzU7ZReJqn9bHOPEH4Q63K5NMsktLfdbkdVVZVsmcrKyhCJRBAIBMR2lpYpFouJfpaXl48rEyEEiR1f7aXtRMjO82OCaXdraf0/mdVOtNxzOBwi34U6MMo9s/uI9DIZ4Z7Ag4qKCjgcDtUyafURHMehuroaiUQihTN6uKfU7xFCxN/N6COuvtqGU08l+L//C6GuLknVTrTcq6ioQDQaTakDK/o9o++nH/+Yw4knOnHRRfuhK/pPPL/oCuz55B/A3XEH8MoriD30EEb32gsrVqwQJzeMoiBPDvf7/eLp3B6PR/ZZZ2cnWlpaxqVtbGzE0NAQ+vr6suUuNYSTS7ds2QKn0ynKOY4DIQQlJSUokdyvIsjToSZPJpMIBAIp68hy+lq2CSEpdmh9EdIpDV6l+tI8bDabrH2Ab/P0dXE5P+VsCC+56urqlPRKcpq6yaZcT5salZtl3yj3cskZOW4r+aPmZz5xRoBaHdDYsaqPkMqNco+GB7R9hACt+lKzY4QztGVNlz/wAME99/BLWByXX5zRkpvxfkomk3jooVFceWUFbpn/F5z20s/Bbd8OUl4O3HUXIj/5CSI7BnlSO4FAAPvss8/EPTlcDQIB1Zbj0k/1zjdUV1fLfkpMupTOZjOn2Y3a0ZOORldLx+hzs+opG7DaV8YZuueMM+bnUSic0dLJFmfOOou/MJcQ82OeCoUz//u/CfT1ETwROgUn7fEORg6bD250FNy556KkpQXVO2bW0j+0mJBLdR6PR/GwSa/XK54qnq+YN28eioqKACDlm9o555yDpUuXAjA+FTo8PJyylGB0Gl6Y3g0Gg4aWFjiOw+joKNX0biwWUy2T4IcAuTIJs3lKZYrFYrJT1pFIxLLlErOW6oRvzUIdmL0EJMTWCfZzxT2O4xTbKb1Mo6OjqmUS6litnYLBoGqZOI4Tl3qk7QTA8P9TJu2kl3scx6XUgdlLdZn2EWZwL5FIwGazUZVJq48oLy9PqS8j3FPq9ziOQygUMq2PKC8Hrr++BN9+y+HOO0tM456w7JneTmb1e2a8nxKJBAghKCsL4oEHgDVrarHvb3vwxAk34uje5eAeeQSDa9fikunT4S0rE/vPhI5R5oRbqgP4WCCv1ysbkM1xnGK6XEP4x/X5fLJLdekwOhUajUbFf2IlfZppeKkdPUt10WhUcfYsfepbyENtCj0SiaSUR8lPpaW6WCyG4uLicUt1cnKausmmXE+bGpWbZd8o93LJGaWlOjl/1PzMJ84IUKsDGjtW9RFSuVHu0fCAto8QoFVfanaMcIa2rEryQAA48kigrY3DWWflB2e05Ga8n+RsfPYZh3POIdj/21dx21c/gf3rL0DKysD94Q8gZ54JgP+yU1NTs2su1QHAkiVLZO+k6+/vh8vlystBkxRCJyN85GSZyMfGxjT1aWxL7ejxRQi2ptEX8lDTTy+Pkp9KNsLhsC652e2RiVxPm1rNGSu5lyvO6PVHSZ5PnKGtAyPtpLcereQeDQ/M5IyWHb2cybS9XS4Ozz3H4corgc8/zw/OmNXeem3svTewfj2HWT8/BrOTb+GrQxaCC4eBs88Gd9554Hb09bQoyIGT1jlMHo8HixcvTjkiAOBnovQcC5ArCN+ChI+czGp5LvJkZWJlyjdfWJlYmQpZXldH8M47wB57EESj+eljttqJ4wguvpjg0Z4paIg9j6fm/A7EZgP++EeQww4D+fhj0KKgYpza29vh9/vF2aTm5mbxNO70C307OzuxcuVKrFy5Ei6XC4ODg2hvb6c+vTuXmDt3Lux2vmmEETMh5sQ4seMI2HEE7DgCdhwBO45gYh5HINdOkyc78cYbSfzylwTPPRdCWZkx7uX7cQTp3BNk6WVav341ovFHcO4nbVhlewSP4ueY8v77ePPYY0GLgoxxmqjIVozTyMiIGNyqpK9lW3iZCXZofRHSVVZWytaBVF+ah1rsgdBBSqda5fxUinEaHR1FeXl5SnolOU3dZFOup02Nys2yb5R7ueSMHLeV/FHzM584I0CtDmjsWNVHSOVGuUfDA9o+QoBWfanZMcIZ2rJqyfl6JDjnHH6X3YMPApyBYwpo6iBT3814P+nhzJ//zOHqpV9i7ZSLMOXhazDJM5sqxokNnPIIwsCJpuEYGBgYGBhoMTYGHHcc8OtfA83NufYmf/DFF8C11wLXXhvE7rvTvX/ZwCmPkK0Zp7GxMZSWlmY84yS1o2fGaWxsDGVlZbJ1kP4NTshD7ZtgOBxOKY+Sn0ozTpFIRNwSrCWnqZtsyvW0qVG5WfaNci+XnFGacZLzR83PfOKMALU6oLFjVR8hlRvlHg0PaPsIAVr1pWbHCGdoy6oll8p8PqC6GkgmAYcj+5zRkpvxfjLKGT276goqxmlXgdUxTsPDw6isrATHZR7jVFlZaSjGqaysjDrGacqUKaoxTkNDQynTsnLxC5WVlaoxTkVFReNinAgh4DiuIGKctm/fLraptJ3Mip0ZGhoS7Webe8Jyxh577EEd4+RyuVRjnILBICKRiMgZuRinyspK1RgnwY70xZlMJlFaWloQMU4+ny9lKdqqGCe9fUR6mYxwT+DBbrvtRh3jpNRHcBwf4zQ6OprCGT3cU+r3COFjnLLRRwwPD8NmA7ZtA44/vgp//nMRZszQF+M0PDycUgdW9HtmvJ+k/WF6mVavXo37779fLK8weJb6qQU245RHyNaME7tyZScIYVeuaMnZlSt0/qj5mU+cEcCuXKHrIwQUypUrNO308MPAihUcXn+dIH1yxUrOaMnNunLFCGfYjFOBQ+hk0mVKunrl0o+avpbtdDu0vujR17IvzAqp1Rmtj0p1YKRusinX06ZG5WbZN8q9XHFGyY5eOa3v2ZLT1AGNHav6CKncKPf05GsGZ7TsGOGMWfJ02ZlnAm+/DVx5JYd77qGzYRZntORmvJ+McEbpuRzYwGkXhMPhyKkdPelodLV0jD43q56yAat9ZZyhe844Y34ehcIZLZ1848zKlcDwMBCPA3bKkcBE5wwt2FJdHkFYqvP7/SlThRyXmyWgfFpaYGViZSpkeT75YpY8n3wxS55Pvpgl19L90Y+Ac84BfvSj/PPdaJmMyIPBIFwuF1uqK1TU19dbFhweCAQwNjYm7hbLJPBT2HWmNzg8mUyitraWKjg8EonA5XKpBodv27Yt5V6i9MBPwU+l4HCHwwG73T4uOLy4uBg2my3vg8MTiQS+/fZbsU2l7WRG0PHo6CiGh4dF+9nmHiF8AO3UqVOpgsMjkQgqKipUg8OHhoZgs9lEzqQHhwucUQoOLyoqQlVV1bjgcCGIthCCw7/55puUuxitCA430kdIy2SUewIPamtrqYLD1foIjuODw4eHh5FIJMT60sM9pX5PiLGqqanJSnC4XDtdfnkYTU2VmDYthAMOgGpwuM/nS1kCM7vfM+P9FIvFMDQ0JNpgweETHCw4PFVfmgcLDlduJxYcbh1n5Lit5I+an/nEGQEsOJyujxAwkYLD0+Vr1vBLdxs3AjabdZzRkrPgcAbDkI7opTIlXb1y6UdNX8t2uh1aX/Toa9knhAWH62lTo3Kz7BvlXq44o2RHr5zW92zJaeqAxo5VfYRUbpR7evI1gzNadoxwxiy5lu5ppwHHH8/LEgmgqIgFh6uhIC/5ZWBgYGBgYDAPU6fys05XX51rT/IfbOC0C6KkpCSndvSko9HV0jH63Kx6ygas9pVxhu4544z5eRQKZ7R0CoEzra1Adzfw5JPyzyc6Z2jBYpzyCGxXHSuTXnk++WKWPJ98MUueT76YJc8nX8yS55MvZsn12nj/fQ4LFxK88w4wZcrEKBONnO2qK3BYvatudHQUZWVl4LjMdtWFw2GUlZXp3lUXj8cxefJkql114XAYtbW1qrvqvvrqq5R7idJ3zAh+Ku2qs9vtcDgc43bVFRUVwW63F8Suum3btoltKm0ns3bVBQIB0X62uUcIf5/XtGnTqHbVhcNhVFVVqe6q+/bbb2G320XOpO+qEzijtKtOqOP0XXU2m03cbadWJgCmt5PeXXVff/11yl2MVuyqM9JHSMtklHsCDyZPnky1q06tj+A4fledoCvUlx7uKfV7hOzcbZerXXXp3Dv44Gq89loUDkcYPh9gs+3cVZe+G9WKXXWZvp9isRi+/fZb0QbbVTfBwXbVpepL82C76pTbie2qs44zctxW8kfNz3zijAC2q46ujxAwkXfVycl7e4G77waeeoofPNHUQab5sl11DIYhdDLpMiVdvXLpR01fy3a6HVpf9Ohr2Re+AarVGa2PSnVgpG6yKdfTpkblZtk3yr1ccUbJjl45re/ZktPUAY0dq/oIqdwo9/TkawZntOwY4YxZciM25s/ng8Wvuw5YvpztqpOCBYczMDAwMDAwpMBuB554AnjoIeAf/8i1N/kFNuO0C6KsrCyndvSko9HV0jH63Kx6ygas9pVxhu4544z5eRQKZ7R0CpEzkyYBr7wCTJ8OEDLxOUMLFuOURxBinGjWWBkYGBgYGLKBL78EFi0CXnoJqKrKtTfWQM/7ly3V5SGEXSHCB0DK35nIhcA5aR5y+lq20+3Q+iKkk9NN15fmoWRfCK5UqjMtG8lkEsFgcFx6JbnZ7ZGpXE+bWs0Zq7iXS87o8UfNz3ziDE0d0HJPTddoH2EG92h4QNtH0NaXmh0jnDGrvc2wMW0aQX09wRlnxJBI5Pf7KRPO0IIt1eUh5s6da9lxBMFgUHzGcZkdRxAKhZBMJnUfRxCJRFBdXU11HEEoFEJJSYnqcQTS8gDjtxoLfiodRwAA8Xh83HEEhPCXy+b7cQTSQZ7AF7O3uUvtZ5t7hBCMjIygurqa6jiCUCgkHgugdBxBOBxO4Uz6cQQCZ5SOI0gkEgAw7jgCobMuhOMIhDIKdWDFcQRG+oj0MhnhnsADoaxaxxGo9REcxx9HEIvFEAgExPrSwz2lfk/oY6qrq/PmOAKldrr55goccwzB7beHce65sXHtZFa/Z8b7ScoZdhzBBAc7joB+e68U7DgCdhyBVZyR47aSP2p+5hNnBLDjCNhxBHo58/HHAUyeXI3aWnYcAUOeoRC2e6bb0uOLHn0t+4Sw4wj0tKlRuVn2jXIvV5xRsqNXTut7tuQ0dUBjx6o+Qio3yj09+ZrBGS07Rjhjltwszuy2G1BezmHhQg5//COw117m2RfkZryfjHBG6bkcWIzTLoiKioqc2tGTjkZXS8foc7PqKRuw2lfGGbrnjDPm51EonNHSmSicKS4GfvADoKkJkKyumppHPthQA1uqyyOwXXUMDAwMDPkOQoCf/Yw/ruDOO3PtjTlgu+oKHFaOZYU19UzzMGpHTzoaXS0do8/NqqdswGpfGWfonjPOmJ9HoXBGS2eicYbjgNWrgQsusC6PbNvQo88GTrsgzPonNWpHTzrajs+K54XQmQmw2lfGGbrnjDPm51EonNHSmWicKS8H9t+fHzz191uTRy5tqIENnBgYGBgYGBgM4aSTgOZmYGgo155kD2xXXR5C2CIrgOPM3ZpMY5/WhqBD64s0jZa+NA8l++k21fyUs6HkC62PuZbradNM5GbZN8K9XHJGjz9qfuYTZ+T8NWLHqj5CyY4eH2l4QNtHCNCqLzU7RjhDW1YtudWc+cEPOLz+OsGZZwLPPQcIm9Py4f2kp6xKZZYDGzjlIaw+AFM454LjMjsAU7Cj9wBMoWw0B2Amk0lEo1HVAzAJISkH06Ufbif4qXQAZlVVlewBmFVVVQVzAKa0TaXtZNbBilL72eaetBOkOQAzmUxiZGRE9QBMu92ewpn0AzAFP5UOwBTuwko/ALOiogKEFMYBmBzHpdSBFQdgGukj0stkhHsCD+LxONUBmGp9BMfxB2CWlZVpHoCpxD2lfo8QAofDAQB5fwBmRUUFHA5HSh0I7bRsWRgvv5xEMMj7n+v3k9QGOwBzgkOI6vf7/SlR/VbMHqSfb2HkW4rUjh5fCCGah9JJdaWDRzl96enHan6qfZPSg3yaPdDTppnIzbJvlHu54oxef5Tk+cQZQa5VB1pyK/sIqdwo92h4QNtH0NaXmh29nNFTVjV5Njlz993AvvsCJ5yQ2/eTEc4Eg0G4XC62q66QIfzzSTud9I8ROQDx27CavpbtdDu0vgjppB2Mkr40DyX70m/3Wn4qlVP4NkMjN7s9MpXraVOrOWMV93LJGT3+qPmZT5yhqQNa7qnpatWjldyj4YGcn5nUl5odI5wxq72zyZl58ziccw6HTz/N3fspE87Qgg2cGBgYGBgYGDLG4YcDV13FB4tbcThmvoDFODEwMDAwMDCYggsvBKZOBewTeHTBZpwYGBgYGBgYTAHH8TNOf/sb8MQTufbGGrDg8DxCtq5cSQ+cy7YdPelodLV0jD43q56yAat9ZZyhe844Y34ehcIZLZ1djTP//S9w7LFAby9w0EHW5GGmDXblSoHDyrEsIUTcrpsLO3rS0ehq6Rh9blY9ZQNW+8o4Q/ecccb8PAqFM1o6uyJn9tsPuOMO/jLgtBNgTMvDTBt69NnAaRdE+jlG2bajJx2NrpaO0edm1VM2YLWvjDN0zxlnzM+jUDijpbMrcmbJEuDyy4FEwro8rLKhBjZwYmBgYGBgYLAE554LbNsGPPZYrj0xD2zgxMDAwMDAwGAZSkuByy4D3ngj156YAzZw2gVhViCiUTt60tHoaukYfV4oAZuA9b4yztA9Z5wxP49C4YyWzq7Mmb33Bv74R37pbvt2a/Iw24aqfbarLn+QrStXaOS5yNNqeT75YpY8n3wxS55PvpglzydfzJLnky9myfPJF7Pk+eTLvfcCCxdymDUr/8qk58qVCXxEVeGivr7eskt+A4EA4vE4ioqKwHGZXfKbSCRQVFSk+5LfoqIiVFVVUV3ym0gkUFlZqXrJ79DQEGw2m/gtI/0CT8FPpUt+hYtZ0y/5LS8vRyKRyPtLfhOJBPx+v9im0nYy4/LY0dFRhMNh0X62uSfskqmtraW65DeRSKCkpET1kt/h4WEkEgmRM+mX/AqcUbrkV9BPv+RXuEC3EC759fl8Il+k7WTmJb9G+ghpmYxyT+BBdXU11SW/an0Ex/GX/IbDYUQiEbG+9HBPqd8jhIh1UAiX/IZCIcTjcbEO9HLvvPPKMDZWgkcfDeOUU3bqm/l+isViCAaDog12ye8EhzDj5PP54HQ6RbmZo2vh1min0ymS38iInhCSYofWFyGdy+WSrQOpvjQPpUswAcDv96eUR8lPpQs2g8EgqqurU9IryWnqJptyPW1qVG6WfaPcyyVn5Lit5I+an/nEGQFqdUBjx6o+Qio3yj0aHtD2EQK06kvNjhHO0JZVS55PnAGAb77hUF9P8OCDwPz5qfpmvJ+MciYQCKCmpobNOBUqpN8CpTIlXb1y6UdNX8t2uh1aX/Toa9knZPylmUp+auWpVAdG6iabcj1talRuln2j3MsVZ5Ts6JXT+p4tOU0d0Nixqo+Qyo1yT0++ZnBGy44RzpglzyfOTJ0KPPIIh5/9DNi4Edhtt1R9M95PRjij9FwOLDicgYGBgYGBIWs49ljg178GPvww154YAxs47YKw2cxpdqN29KSj0dXSMfrcrHrKBqz2lXGG7jnjjPl5FApntHQYZ1Jx6aXAMccAr7xifh5W1ymLccojZOuuOgYGBgYGhlzD5wNmzwZWrwYWLsytL+yuugKHlWNZQggikUjGeRi1oycdja6WjtHnZtVTNmC1r4wzdM8ZZ8zPo1A4o6XDOCOPmhrg8ceBn/8c+OKL3HOGFmzgtAtCus00F3b0pKPR1dIx+tysesoGrPaVcYbuOeOM+XkUCme0dBhn5HHkkfyp4i+9ZF4eVtcp21XHwMDAwMDAkDP86lf8z//+N3WXXb6CzTgxMDAwMDAw5BRjY8D3vw9s2JD/8zls4LQLQjiVPFd29KSj0dXSMfrcrHrKBqz2lXGG7jnjjPl5FApntHQYZ9RRWgo88gjwy1+W44svMrNldZ2yXXUG4ff7sXbtWnR3d6Onp0f3czmwXXUMDAwMDLsybroJCAaB667Lbr5sV53F6O3txdq1azE4OAiv16v7uRas3uki3KuUCzt60tHoaukYfW5WPWUDVvvKOEP3nHHG/DwKhTNaOowz9HlceGEYy5cT+HzGbRjlDC3YwMkAGhoa0NLSgrq6OkPPcw3pZZq5sKMnHY2ulo7R52bVUzZgta+MM3TPGWfMz6NQOKOlwzhDh1iMz+Poo4HeXmM2rPazcBZXJfB6vWhubsbq1avh8XgU9bq6uuD3++FyuTA4OIjGxkY0NDRk0VMGBgYGBgYGPbDZgFWrgNNO4++zmzYt1x6loqBmnFpbW9Hc3IzOzk709/dr6gJAW1sbWlpa0NHRgc7OTvQaHcJOIKxevTqndvSko9HV0jHyPBKJYMWKFQXzbdCsNrXaPuNM/sBqzpiVR6FwRkuHcUZfHkcdBVx0EfCb3xi3YRlIAcLn8xEApK+vT/Z5X18fcbvd4+SDg4PE5XKZ5kdnZ6dsPrTP0xEIBAgA4vf7zXBPFslkknznO98hyWQyJ3b0pKPR1dIx+tzv91veFmbBrDa12j7jTP7Aas6YlUehcEZLh3HGWB6JBCHDw4TEYsZt0EJoi0AgoKlbkEt1Wujs7JRdknO73QD44G3h+bp167BmzRpNmx0dHWJ6q8FxnKW2hU8u7OhJR6OrpWP0ufC3lW1hFsxqU6vtM87kD6zmjFl5FApntHQYZ4zlYbMBFRXAggXAtdfycU9W+alHf0IOnHp7e9HU1CT7zO12o6enRxw4NTU1KermCsTiXQvCJ5N/AKN29KSj0dXSMfpcaAMr28IsmNWmVttnnMkfWM0Zs/IoFM5o6TDOGM+D44CrrgJ+9jOgvx+YNMkaP/W0QUHFONHC6/Uq7mirra3VjI+a6DDrn9SoHT3paHS1dDJ9XgiwugyMM8b9yVdkowxm5FEonNHSYZwxnseCBcAZZ/AxT0ZtmIkJOeOkhaGhoVy7wMDAwMDAwECJ5cuBrVtz7QWPCTdw8vv9qs9dLpehQyml6O/vR29vL9asWQOv14v29nbU1dWhpaWF6rkShFHy1q1bEQwGRTnHcSCEoKSkBCUlJePk6VCTJ5NJJJNJBAKBlPX1dH0t24SQFDu0vgjppOVT0pfmYbPZFL9FpJdHyU85G4QQJBKJcekDgUDKTz11k025njY1KjfLvlHu5ZIzctwWOJPuj8AVOT/ziTMC1OqAxo5VfYRUbpR7NDyg7SMEyPUTtHkqcViNM7Rl1ZKbyRmtOsg0X632rqkBTj8daGnhMHu2MRuRSETcxSiVC21AM1tVkFeu+P1+1NTUoK+vb9w5TsKzzs5O2YFKc3MzvF4v+vr6suUuNbZu3Yq99tor124wMDAwMDDskvjss8+w5557qupMuBknLWjNSOUS06dPx+DgIIqLi2WD2tJnnBgYGBgYGBjoIZ1xkoIQglgshunTp2vamHADJ5fLBUB9gFRbW5sdZ3TCZrNl7cgDBgYGBgYGBv2YkLvqPB4Ptm/fLvvM6/WisbExyx4xMDAwMDAwTARMyIFTQ0ODYgC41+tl99UxMDAwMDAwGMKEW6oDgCVLlmDBggXj5P39/XC5XKoXAzMwMDAwMDAwKKEgZ5y0zmHyeDxYvHgxurq6UuTt7e3o7u620jUGBgYGBgaGCYyCOo6gvb0dfr8fvb298Hq9cLvd4rJbZ2fnOP2VK1cC4APGBwcH0djYyJbpGBgYGBgYGAyjoAZODAwMDAwMDAy5REEu1TEwMDAwMDAw5AJs4MTAwMDAwMDAQIkJuauOgcFq+P1+rF27Ft3d3ejp6cm1OwwFgPb2dgD8kSi1tbXo6OgQD+xlYEiH3+9HV1eXGKPr9XrR0dHBDknOA7CBEwODTgibE4TOjIFBC62trSkDpebmZixYsCAv78xkyA8sXboUjY2N4p2rra2taGxsxODgYI49Y2BLdQwMOtHQ0ICWlhbU1dXl2hWGAoB0J7CAZcuWob+/H/39/Tn0jCHfIR1Y19XVsS9qeQI248SwS8Dr9aK5uRmrV69WPQC1q6sLfr+fHWHBYCpnhoaG4PV6RTvCzJNUxlD4MJMz6WcO9vT0oKmpyRK/GfSBDZwYJjRaW1sxNDQEt9ut+e2+tbUV9fX1aGtrE2XNzc0AwAZPuxDM5ozL5YLP50tJJ9hlvJoYsLqf6erqgtvtlj2vkCEHIAwMuwB8Ph8BQPr6+mSf9/X1EbfbPU4+ODhIXC6XbJrOzk7ZNAwTA1ZwRoDb7SYdHR2m+MmQP7Cqn+no6CAtLS1kcHDQVH8ZjIHFODEwgD95Xu7bnrCDpbe3N9suMeQ5jHKmvb0dTU1NKTMODLsGjHCmpaUFbW1tqKurY3FOeQI2cGJgAN9hKW0Nd7vd7MgBhnEwwpmuri5MmjQJHR0dFnvHkI+g5Yzf70ddXV3KQEoYcK1bt85yPxnUwQZODAzggzqVdsnV1tay3U8M46CXM729vaitrRVnmoTddgy7Dmg54/V6x11mL8w0sc0EuQcbODEwUCC9E2Ng0IKUM/39/eju7obL5UJvby96e3tx4403/v/27u8ocTaKA/DZma+ADHYAHUStAOhApAMoITVgB9gBDh1IB47pgHQApgP2wpHR1U9fFhU3PM9lcnMuzmR+vH8OhhnywlPP5Hkel5eXL7b1ZrNZ9Ho9Fwp+ALfqOHp1Xb/7PsuyF+cKyrKMxWIRs9ksqqqKoiii0+lsB9XRfLv0TF3X0e12t5Ogn7Nldzx2/c5MJpMoiiJOTk5itVpFq9V6NaKAwxCcIMHzcwl5nkee5w738q6nnnlrHAG85fl3JssywfqHslUHH/jolyL8Sc+wKz3z7xCcOHpPv/Le+3C1Wq3vKYZ/gp5hV3qmOQQniMftt9Vq9ea7qqqi3+9/c0X8dHqGXemZZhCcIB5npPzfYLmqqtxk4RU9w670TDM4HA4RMRwOo9vtvnpelmVkWWZ2Cq/oGXalZ5rBihNH4aM5TE9zU/68Ll4UhSvAR0rPsCs9cxx+bTabzaGLgK9SFMV2QnNVVdFut7fL4W/90/jV1VVEPB7kXC6X0e/3LZ8fGT3DrvTMcRGcAAAS2aoDAEgkOAEAJBKcAAASCU4AAIkEJwCARIITAEAiwQkAIJHgBACQSHACAEgkOAEAJBKcAAAS/XfoAgCapK7rKIoisiyL8/PzuLi4OHRJwCey4gTwibrdbhRFEZPJJO7u7g5dDvDJBCeATzIYDKLX60W73T50KcAXsVUH8AnKsoz5fB4PDw+HLgX4Qr82m83m0EUA/Ov6/X602+2YTqcREXF9fR1nZ2eR5/mBKwM+k606gD2VZRmLxSIGg8H22XK5FJqggQQngD1Np9PIsix6vV5EPK42jcfjA1cFfAXBCWBPNzc329BU13XUde2AODSUw+EAe6iqKuq6juFwGBERRVFszzkBzWPFCWAPi8UiIiLyPN/ObwKaS3AC2MP9/X1kWRbT6TSGw2FkWXbokoAvJDgB7GG9Xm+36tyig+YTnAD2UFVVjEYjoQmOhOAE8JeKojh0CcA3E5wA/sJ8Po/xeBzD4TDW6/Wr90+HxoFm8ZcrAHuo6zpOT09juVxun43H4xiPx7bvoIEEJ4A9zefzuL29jU6nExERo9HI7TpoKMEJACCRM04AAIkEJwCARIITAEAiwQkAIJHgBACQSHACAEgkOAEAJBKcAAASCU4AAIkEJwCARIITAEAiwQkAIJHgBACQSHACAEj0G44KWhDsoy3bAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "label_size = 17\n", "title_size = 18\n", "legend_size = 13\n", "handle_length = 1.5\n", "fig, (ax1) = plt.subplots(1,1,figsize=(6,4))\n", "\n", "ax = ax1\n", "ax.tick_params(axis = 'x',which='both',length=5,direction='in', pad=10)\n", "ax.tick_params(axis = 'y',which='both',length=5,direction='in', pad=5)\n", "ax.xaxis.set_ticks_position('both')\n", "ax.yaxis.set_ticks_position('both')\n", "plt.setp(ax.get_yticklabels(), rotation='horizontal', fontsize=label_size)\n", "plt.setp(ax.get_xticklabels(), fontsize=label_size)\n", "ax.grid( visible=True, which=\"both\", alpha=0.2, linestyle='--')\n", "ax.set_xlabel(r'$\\ell$',size=title_size)\n", "ax.set_ylabel(r'$C_\\ell^{\\kappa\\gamma}$',size=title_size)\n", "ax.set_xscale('log')\n", "ax.set_yscale('log')\n", "ax.set_ylim(1e-11,1.5e-7)\n", "ax.set_xlim(2.,5e3)\n", "fac = ell*(ell+1.)/2./np.pi\n", "ax.plot(ell,cl1h/fac,label='CLASS_SZ 1h',c='b',ls='-',lw=0.7)\n", "ax.plot(ell,cl2h/fac,label='CLASS_SZ 2h',c='b',ls='--',lw=0.7)\n", "ax.plot(ell,cl1h/fac+cl2h/fac,label='CLASS_SZ 1+2h',c='r')\n", "\n", "\n", "ax.legend(loc=1,ncol = 1,frameon=False,fontsize=11)\n", "\n", "fig.tight_layout()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.12.8" }, "widgets": { "application/vnd.jupyter.widget-state+json": { "state": {}, "version_major": 2, "version_minor": 0 } } }, "nbformat": 4, "nbformat_minor": 4 }