{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Quick Start\n", "\n", "This is a quick start reference for an overview of what the code can do. For more details, please visit the other parts of the documentation." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "# preamble\n", "from __future__ import annotations\n", "\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import pyccl as ccl\n", "from init_halo_model import ( # halo model\n", " a_arr,\n", " bM,\n", " cM_relation,\n", " cosmo,\n", " hmc,\n", " hmd,\n", " k_arr,\n", " r_arr,\n", ")\n", "\n", "import glasz\n", "\n", "# Initialize the halo model\n", "\n", "# CMASS PARAMETERS\n", "z_lens = 0.55 # Mean z for CMASS\n", "a_sf = 1 / (1 + z_lens)\n", "\n", "# constituent fractions\n", "fb = cosmo[\"Omega_b\"] / cosmo[\"Omega_m\"] # Baryon fraction\n", "fc = cosmo[\"Omega_c\"] / cosmo[\"Omega_m\"] # CDM fraction\n", "\n", "# define 2-halo term\n", "xi_mm_2h = glasz.profiles.calc_xi_mm_2h(\n", " cosmo, hmd, cM_relation, hmc, k_arr, a_arr, r_arr, a_sf\n", ")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Profiles" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Halo Mass\n", "M_halo = 3e13\n", "\n", "# 2-halo Amplitude\n", "A_2h = 0.7\n", "\n", "# GNFW Parameters\n", "alpha = 1.0\n", "beta = 3.0\n", "gamma = 0.2\n", "x_c = 0.5\n", "\n", "Rb = 10 * (hmd.get_radius(cosmo, M_halo, a_sf) / a_sf)\n", "\n", "# COMPUTE GNFW AMPLITUDE\n", "prof_nfw = ccl.halos.HaloProfileNFW(\n", " mass_def=hmd, concentration=cM_relation, truncated=False, fourier_analytic=True\n", ")\n", "\n", "prof_baryons = glasz.profiles.HaloProfileGNFW(\n", " hmd,\n", " rho0=1.0,\n", " alpha=alpha,\n", " beta=beta,\n", " gamma=gamma,\n", " x_c=x_c,\n", ")\n", "\n", "prof_baryons.normalize(cosmo, Rb, M_halo, a_sf, prof_nfw)\n", "\n", "\n", "# COMPUTE 3D DENSITY PROFILES\n", "def rho_2h(r):\n", " return (\n", " xi_mm_2h(r)\n", " * bM(cosmo, M_halo, a_sf)\n", " * ccl.rho_x(cosmo, a_sf, \"matter\", is_comoving=True)\n", " * A_2h\n", " )\n", "\n", "\n", "prof_baryons.rho_2h = rho_2h # add 2-halo term to baryon profile\n", "\n", "prof_matter = glasz.profiles.MatterProfile(\n", " mass_def=hmd, concentration=cM_relation, rho_2h=rho_2h\n", ")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Galaxy-Galaxy Lensing (GGL)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "R = np.geomspace(\n", " 1e-2, 6e1, 100\n", ") # Mpc/h comoving units (like observations of Amon & Robertson et al. 2022)\n", "\n", "# COMPUTE ΔΣ PROFILE\n", "ds_b = (\n", " fb\n", " * glasz.GGL.calc_ds(\n", " cosmo,\n", " R / cosmo[\"h\"], # convert from Mpc/h to Mpc\n", " M_halo,\n", " a_sf,\n", " prof_baryons,\n", " )\n", " / cosmo[\"h\"]\n", ") # convert from Msun/pc^2 to h Msun/pc^2\n", "\n", "\n", "ds_dm = (\n", " fc\n", " * glasz.GGL.calc_ds(\n", " cosmo,\n", " R / cosmo[\"h\"], # convert from Mpc/h to Mpc\n", " M_halo,\n", " a_sf,\n", " prof_matter,\n", " )\n", " / cosmo[\"h\"]\n", ") # convert from Msun/pc^2 to h Msun/pc^2\n", "\n", "ds = ds_b + ds_dm" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### kinetic Sunyaev-Zeldovich (kSZ) Effect" ] }, { "cell_type": "code", "execution_count": 19, "metadata": {}, "outputs": [], "source": [ "# COMPUTE kSZ PROFILE\n", "theta = np.geomspace(0.5, 6.5, 20)\n", "\n", "\n", "def rho_gas_3D(r):\n", " \"\"\"\n", " We need to be careful about comoving units here. The kSZ code assumes\n", " that the density profile is in physical units. The array which it will\n", " be feeding into this function is in physical units r = aχ. The CCL\n", " profile assumes comoving units so we need to convert the input to\n", " comoving units before passing it to the profile. Then we need to convert\n", " the density profile back into physical units before returning it by\n", " dividing by the scale factor a^3.\n", " \"\"\"\n", " return (fb * prof_baryons.real(cosmo, r / a_sf, M_halo, a_sf)) / a_sf**3\n", "\n", "\n", "T_kSZ_150 = glasz.kSZ.create_T_kSZ_profile(theta, z_lens, rho_gas_3D, \"f150\", cosmo)" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAABKYAAAEhCAYAAABFk6ghAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAthBJREFUeJzs3XlcjdkfB/DPbS8qSlJUMtaiVLasZYnsjGXsIVsIWbLvZBuMlG0kO2MJM+M39mxZWywVQmQpWdv3e35/nHGJMi339txu3/fr9by4z733Od8eOs9zv/ec7xExxhgIIYQQQgghhBBCCClhSkIHQAghhBBCCCGEEELKJkpMEUIIIYQQQgghhBBBUGKKEEIIIYQQQgghhAiCElOEEEIIIYQQQgghRBCUmCKEEEIIIYQQQgghgqDEFCGEEEIIIYQQQggRBCWmCCGEEEIIIYQQQoggKDFFCCGEEEIIIYQQQgRBiSlCCCGEEEIIIYQQIghKTBFCCCGEEEIIIYQQQVBiihBCCCGEEELycenSJXTr1g3GxsYQiUQ4duxYvq8dM2YMRCIR1q9fX2LxEUJIaUeJKUIIIYQQQgjJR0pKCqytrbFx48Yfvu7YsWO4ceMGjI2NSygyQghRDCpCB0AIIYQQQggh8srZ2RnOzs4/fM2rV68wYcIEnDp1Cl26dCmhyAghRDEofGKqV69eCAwMRLt27XD48GHJ/ujoaIwYMQJv3ryBsrIyrl+/jnLlyhXomGKxGK9fv4a2tjZEIpGsQieEkBLFGENSUhKMjY2hpEQDaqWBrheEEEVD14rvicViDBkyBNOnT4elpWWB3pORkYGMjIxcx/jw4QP09fXpekEIUQiFuV4ofGLK3d0dI0aMwM6dO3Ptd3FxwdKlS9GqVSt8+PAB6urqBT7m69evYWJiIu1QCSFELrx48QLVqlUTOgyFQNcLQoiiomvFFytXroSKigrc3d0L/B4vLy8sWrRIhlERQoh8KMj1QuETU46OjggMDMy1Lzw8HKqqqmjVqhUAQE9Pr1DH1NbWBsBPsI6OjlTiJIQQoSUmJsLExETSx5Hio+sFIUTR0LUit+DgYPz2228ICQkp1EinWbNmwcPDQ/I4ISEBpqamdL0ghCiMwlwv5DoxdenSJaxevRrBwcGIjY1FQEAAevbsmes1vr6+WL16NWJjY2FpaYn169dLEk75iYqKQvny5dG9e3e8fPkSffr0wezZswsc1+eLjo6ODl04CCEKh6YQSA9dLwghioquFdzly5cRHx8PU1NTyb6cnBxMnToV69evx7Nnz/J8n7q6ep4zNuh6QQhRNAW5Xsj1xPD/WgHj4MGDmDx5MubMmYPQ0FC0atUKzs7OiImJ+eFxs7KycPnyZfj4+ODatWs4c+YMzpw5k+/rMzIykJiYmGsjhBBCCCGElG1DhgzB3bt3ERYWJtmMjY0xffp0nDp1SujwCCGkVJDrEVP/tQLG2rVrMXLkSLi6ugIA1q9fj1OnTmHTpk3w8vLK933VqlVD48aNJXU/OnfujLCwMHTo0CHP19MccEIIIYQQQsqm5ORkPH78WPI4OjoaYWFh0NPTg6mpKfT19XO9XlVVFVWqVEGdOnVKOlRCCCmV5Dox9SOZmZkIDg7GzJkzc+13cnJCUFDQD9/buHFjvHnzBh8/foSuri4uXbqEMWPG5Pv6b+eAf54rSQghhBBCCFFst2/fhqOjo+Tx588Fw4YNg7+/v0BREVJ25OTkICsrS+gwyDdUVVWhrKwslWOV2sTUu3fvkJOTA0NDw1z7DQ0NERcXJ3ncsWNHhISEICUlBdWqVUNAQAAaN26M5cuXo3Xr1mCMwcnJCV27ds23rfzmgBNCCCGEEEIUm4ODAxhjBX59fnWlCCGFwxhDXFwcPn36JHQoJB8VKlRAlSpVil13sNQmpj779gQwxnLty29u939NE5SF9PR0qKioQEWl1J92QgghhBCiCBgDqJA5IUQOfU5KVa5cGVpaWrToghxhjCE1NRXx8fEAACMjo2Idr9RmSCpVqgRlZeVco6MAID4+/rtRVPLi999/x5IlS9C/f38MHjwYjRs3pl8uQgghhBAiW8nJQFQU3x49yv1nQADwHytaE0JIScvJyZEkpb6t40bkg6amJgCeg6lcuXKxpvWV2sSUmpoa7OzscObMGfTq1Uuy/8yZM+jRo4eAkeXv77//Rnx8PLy9veHt7Y1atWph8ODBGDRoEH766SehwyOEEEIIIaVVejrw5EneyafY2PzfFxVFiSlCiNz5XFNKS0tL4EjIj3z+98nKylLcxNR/rYDh4eGBIUOGoFGjRrC3t8fWrVsRExODsWPHChh1/k6cOIGzZ89iz549CAgIQFRUFBYsWIAFCxagefPmGDJkCPr16wc9PT2hQyWEEEIIIfImKwt49izv5FNMDJ+Wlx8DA6BWLaB27dx/1qpVYuETQkhh0Qwj+Satfx+5Tkz91woY/fv3x/v377F48WLExsaifv36OHnyJMzMzIQK+YdUVVUlta2Sk5MREBCAPXv24OzZswgKCkJQUBAmTZqErl27YujQoXB2doaamprQYRNCCCGEkJIiFgMvXuSdfIqOBrKz83+vjg5POOWVfKpQocR+BEIIIaQwRKwwS0wQAEBiYiJ0dXWRkJAAHR2dYh8vNjYW+/fvx65du3Dnzh3Jfn19ffzyyy9wcXGBnZ0dZYsJITIl7b6N0DklhOSDMSAuLu/k05MnfFpefjQ1v086ff67gYHMC5lTvyYbdF4JyS09PR3R0dEwNzeHhoaG0OGQfPzo36kw/Zpcj5gqK4yMjODh4QEPDw/cu3cPu3btwp49exAXFwcfHx/4+PjAwsICLi4uGDx4cLEr3hNCCCGEkBKQmAhERHyffIqK4gXJ86OqCvz0U95T74yNASWlkvsZCCGEFIqDgwMaNmyI9evXCx1KqUGJKTnToEEDrF69Gl5eXjh79ix27tyJY8eOISIiAjNmzMCsWbPQqVMnDB8+HN26daOpfoQQQgghQktJ4Qmo8HC+3b/P/3zxIv/3KCkB1avnnXwyNQVU6DadEEJI2UBXPDmloqKCTp06oVOnTkhISMAff/wBf39/BAUF4e+//8bff/8NfX19DBo0CCNHjoSVlZXQIRNCCCGEKLa0NCAy8ksC6vMWHZ3/e4yMgLp1v08+mZsD6uolFzshhJBSiTGGnJwcqCjwFxY0DrgU0NXVxahRo3D16lU8fPgQM2fOhLGxMd6/f48NGzbA2toajRo1wqZNm5CQkCB0uIQQQgghpVtGBnD3LrB/PzBnDtCzJ08mlS8P2NkBQ4cCK1cCf/31JSlVuTLg6AhMnAhs3gxcvgx8+AC8fg2cP8/3TZ0KdOvGE1WUlCKEkMJhjI9QFWIrZGnu7OxsTJgwARUqVIC+vj7mzp2Lz+W99+zZg0aNGkFbWxtVqlTBwIEDER8fL3lvYGAgRCIRTp06hUaNGkFdXR27d++GkpISbt++nasdb29vmJmZSY598eJFNGnSBOrq6jAyMsLMmTOR/dWiGQ4ODnB3d8eMGTOgp6eHKlWqYOHChbmOuXDhQpiamkJdXR3GxsZwd3cv1M9eFIqbclNQtWvXhpeXF5YsWYIzZ87Az88Px48fR3BwMIKDgzF16lT069cPrq6uaNGiBRVMJ4QQQgjJT1YWr/v09ein+/eBx4+BnJy836OvD1haftnq1+d/VqpUsrETQkhZk5rKvyAQQnIyUK5cgV++c+dOjBw5Ejdu3MDt27cxevRomJmZYdSoUcjMzMSSJUtQp04dxMfHY8qUKXBxccHJkydzHWPGjBlYs2YNatSogQoVKqB9+/bYsWMHGjVqJHnNjh074OLiApFIhFevXqFz585wcXHBrl278ODBA4waNQoaGhq5kk87d+6Eh4cHbty4gWvXrsHFxQUtWrRAhw4dcPjwYaxbtw4HDhyApaUl4uLici3QJiu0Kl8RyNuqGW/fvsWePXuwfft2hIeHS/bXq1cPo0aNwtChQ6Gvry9ghISQ0kDe+jZFQOeUEDmRnc1Xu/u6/lN4OE9KZWXl/R5d3dyJp8+boaHMV76TZ9SvyQadV0Jyy3O1t5SUUpGYcnBwQHx8PMLDwyUDRWbOnIkTJ04gIiLiu9ffunULTZo0QVJSEsqXL4/AwEA4Ojri2LFj6NGjh+R1f/zxB8aOHYvY2Fioq6vjzp07sLGxwdOnT1G9enXMmTMHR44cQWRkpKRdX19feHp6IiEhAUpKSnBwcEBOTg4uX74sOW6TJk3Qtm1brFixAmvXrsWWLVtw//59qKqq/ufPKq1V+WgqnwIwMDDAlClTcO/ePVy7dg0jRoxAuXLlEBkZCQ8PD1StWhWDBw/G5cuXQXlIQggpnF69eqFixYro06eP0KEQQgri3Tvg9GlgxQpg0CDA2pp/kKlbF/j5Z2DBAuCPP3hiKiuLP9e0KTBiBPDrr8CpU8DLl8DHj8DVq8CWLYC7O9CuHVClSplOShFCiKC0tHiCSIhNS6tQoTZr1izX7CV7e3tERUUhJycHoaGh6NGjB8zMzKCtrQ0HBwcAQExMTK5jfD0yCgB69uwJFRUVBAQEAAD8/Pzg6OiI6tWrAwAiIyNhb2+fq90WLVogOTkZL1++lOz7tj61kZGRZCph3759kZaWhho1amDUqFEICAjINRVQVmgqnwIRiURo1qwZmjVrhnXr1mH//v3YunUrQkJCsHfvXuzduxcWFhYYN24chgwZAl1dXaFDJoQQuefu7o4RI0Zg586dQodCCPnWmzdASAgQHPzlz29u7CW0tIB69b4fAWVqSskmQggpDUSiQk2nk0fp6elwcnKCk5MT9uzZAwMDA8TExKBjx47IzMzM9dpy3/ysampqGDJkCHbs2IHevXtj3759WL9+veR5xth3pXw+D0z5ev+3I6FEIhHEYjEAwMTEBA8fPsSZM2dw9uxZuLm5YfXq1bh48WKBRlAVFSWmFJSOjg7GjBmDMWPG4Pbt29iyZQv27duHiIgITJw4EZ6enhg0aBDGjx8Pa2trocMlhBC55ejoiMDAQKHDIKRsYwyIjc2dgAoJAV69yvv1NWvyIuXW1l8SUdWrA0o0WYAQQojsXb9+/bvHtWrVwoMHD/Du3TusWLECJiYmAPBdQfMfcXV1Rf369eHr64usrCz07t1b8pyFhQWOHDmSK0EVFBQEbW1tVK1atcBtaGpqonv37ujevTvGjx+PunXr4t69e7C1tS3wMQqLrs5lQKNGjbBt2za8fv0aGzduhKWlJVJTU7Ft2zY0bNgQLVu2xIEDB5CVX40FQggppS5duoRu3brB2NgYIpEIx44d++41vr6+knnxdnZ2uebcE0IEwBgf9XTsGDBvHtClC2BkBFStCnTvDixcCPz5J09KiUR8it6gQXwa3oULwKdPQFQUcOAAMGsWXwWvRg1KShFCCCkxL168gIeHBx4+fIj9+/fD29sbkyZNgqmpKdTU1ODt7Y2nT5/ixIkTWLJkSYGPW69ePTRr1gyenp4YMGAANDU1Jc+5ubnhxYsXmDhxIh48eIDjx49jwYIF8PDwgFIBr4H+/v7Yvn077t+/j6dPn2L37t3Q1NSEmZlZoc9BYdCIqTJEV1cX48ePh5ubGy5fvgwfHx8cPXoUV69exdWrV2FkZIQxY8Zg7NixMDQ0FDpcQggptpSUFFhbW2P48OH4+eefv3v+4MGDmDx5Mnx9fdGiRQts2bIFzs7OiIiIgKmpaaHby8jIQEZGhuRxYmJiseInROExBjx79v1IqHfvvn+tkhJgYQHY2vLRULa2QMOGwhXCJYQQQvIxdOhQpKWloUmTJlBWVsbEiRMxevRoiEQi+Pv7Y/bs2diwYQNsbW2xZs0adO/evcDHHjlyJIKCgjBixIhc+6tWrYqTJ09i+vTpsLa2hp6eHkaOHIm5c+cW+NgVKlTAihUr4OHhgZycHDRo0AB//vmnzBdTo1X5ikCRVs2IjY3F1q1bsXnzZsTFxQHgc1f79++PSZMmwc7OTuAICSElRZH6tryIRCIEBASgZ8+ekn1NmzaFra0tNm3aJNlXr1499OzZE15eXpJ9gYGB2LhxIw4fPvzDNhYuXIhFixZ9t19RzykhhSIW85Xxvk5AhYTwIuPfUlHh0+8+J6Ds7AArq0IXnyXSp+jXCqHQeSUktx+t9lbWLVu2DAcOHMC9e/eEDkVqq/LRiKkyzsjICAsWLMCsWbNw5MgRbNiwAdevX8fu3buxe/dutGrVClOmTEH37t2hrKwsdLiEECI1mZmZCA4OxsyZM3Ptd3JyQlBQUJGOOWvWLHh4eEgeJyYmSuoHEFLmfPwIBAXxle2uXeNJqLxGEaqpAQ0a5B4J1aABQB9ECCGEEInk5GRERkbC29u7UNP/SgOabE8A8FFSAwYMwLVr13Djxg0MHDgQKioquHz5Mnr37o3atWvD29sbycnJQodKCCFS8e7dO+Tk5Hw3ddnQ0FAyghQAOnbsiL59++LkyZOoVq0abt26le8x1dXVoaOjk2sjpExgDIiOBnbvBsaO5QXH9fSArl0BLy8gMJAnpdTVgSZNgHHjgG3beLIqKQm4fRvYuhUYMwZo3JiSUkSu/KheYVZWFjw9PdGgQQOUK1cOxsbGGDp0KF6/fi1cwIQQhTRhwgS0bNkSbdq0+W4aX2lHI6bId5o0aYK9e/di1apV8PHxwebNm/H06VO4u7tjwYIFGDduHNzd3akOFSFEIeS1rO7X+06dOlXSIREi/7KzgTt3gCtX+IioK1f4qnnfql0baNkSaN6cJ5zq1QNkuNw0IbLwo3qFqampCAkJwbx582BtbY2PHz9i8uTJ6N69e6FW2iKEkP/i7+8Pf39/ocOQCUpMkXxVrVoVy5cvx5w5c7Bz506sXbsWT548wfLly/Hrr79i2LBhmDZtGmrVqiV0qIQQUmiVKlWCsrJyrtFRABAfH1/sxLuPjw98fHyQk5NTrOMQIjeSkoDr178koq5fB1JScr9GVZVPxWvZEmjRgiejKlcWJl5CpMjZ2RnOzs55Pqerq4szZ87k2uft7Y0mTZogJiamSAtpEEJIWUNT+ch/KleuHNzc3PDw4UMcOXIETZs2RUZGBrZu3Yq6deuif//+CAkJETpMQggpFDU1NdjZ2X33geLMmTNo3rx5sY49fvx4RERE/HDaHyFy7eVL4OBBYOJEXvOpQgXAyQlYvBg4d44npSpUADp3BpYvBy5eBBISeC2p1auBnj0pKUXKrISEBIhEIlSoUCHf12RkZCAxMTHXRgghZRWNmCIFpqysjN69e6NXr164cuUKVq5cib///ht//PEH/vjjD3Ts2BFz585Fy5YthQ6VEEIA8CKRjx8/ljyOjo5GWFgY9PT0YGpqCg8PDwwZMgSNGjWCvb09tm7dipiYGIwdO1bAqAkpYTk5QHj4lyl5V68Cz59//7rq1b+MhmrZErCwAJToO05Cvpaeno6ZM2di4MCBP6wz6OXllecqroQQUhZRYooUmkgkQqtWrdCqVSvcu3cPK1euxIEDB3Dq1CmcOnUKrVu3xrx589CuXbvvarcQQkhJun37NhwdHSWPP6+YN2zYMPj7+6N///54//49Fi9ejNjYWNSvXx8nT56EmZmZUCETInuMAZGRwOnTwJkzPBGVkJD7NUpKQMOGXxJRLVoAVasKEi4hpUVWVhZ++eUXiMVi+Pr6/vC1tIorIYR8QYkpUiwNGjTAnj17sGjRIqxcuRL+/v64dOkSOnTogObNm2P+/PlwcnKiBBUhRBAODg5gjP3wNW5ubnBzc5Nqu1Rjisidjx+Bs2eBU6d4QurFi9zPly8PNGv2JRHVtCmgrS1MrISUQllZWejXrx+io6Nx/vz5/1yVVV1dHerq6iUUHSGEyDcaf02k4qeffsLWrVvx9OlTTJgwAerq6ggKCkKnTp3QvHlznDlz5j8/HBJCiKKgGlNEcNnZvN7TwoWAvT1QqRLQrx+wfTtPSqmrAx06AGvWAMHBPHF15gywYAHQvj0lpQgphM9JqaioKJw9exb6+vpCh0QIIaUKjZgiUlWtWjV4e3tj9uzZWL16NTZt2oTr16/DyckJLVu2xJIlS+Dg4CB0mIQQQojiefGCj4g6dYqPjvr0Kffz9eoBHTvyrXVrQEtLkDAJKW1+VK/Q2NgYffr0QUhICP766y/k5ORIVnvV09ODmpqaUGETQkipQSOmiEwYGRlh7dq1iI6OxuTJk6Guro4rV67A0dERHTp0wI0bN4QOkRBCCCndUlOB//0PmDyZJ51MTYFRo4DDh3lSqkIFoE8fYNs2ICYGiIgA1q0DOnWipBQhhXD79m3Y2NjAxsYGAK9XaGNjg/nz5+Ply5c4ceIEXr58iYYNG8LIyEiyBQUFCRw5IUQILi4uEIlEEIlEUFVVhaGhITp06AA/Pz+IxWLJ66pXrw6RSIQDBw58dwxLS0uIRCL4+/uXYOTCocQUkakqVapg3bp1ePr0Kdzc3KCqqoqzZ8+iWbNm6NWrF8LDw4UOkRBCCCkdGAPu3ePT7zp0APT0gM6dgd9+Ax484AXL7e359L1r14C3b4FDhwBXV4CKKhNSZJ/rFX67+fv7o3r16nk+xxijWQKElGGdOnVCbGwsnj17hv/9739wdHTEpEmT0LVrV2RnZ0teZ2Jigh07duR67/Xr1xEXF4dy5cqVdNiCoal8pEQYGxvDx8cH06dPx+LFi7Fz504cO3YMJ06cwJAhQ7BkyRJaiYQQojCo+DmRmrQ0Pirqzz950fLXr3M/b2LyZXpeu3ZAxYrCxEkIIYTIGGN8sLAQtLSAwqznpa6ujipVqgAAqlatCltbWzRr1gzt2rWDv78/XF1dAQCDBg3CunXr8OLFC8nnYT8/PwwaNAi7du2S+s8hr2jEFClR1atXh5+fH+7fv4/evXtDLBZj586dqFWrFjw9PfHp23oYhBBSClHxc1IsaWnA0aPAgAGAgQHw88+Avz9PSmlq8ql469bxqXnPn/Open36UFKKEEKIQktN5YvICrFJIyHWtm1bWFtb4+jRo5J9hoaG6NixI3bu3Pnvz5iKgwcPYsSIEcVvsBShxBQRRL169XDkyBFcv34dbdq0QUZGBlatWoWffvoJGzZsQGZmptAhEkIIISUnr2TUgQNASgqvHTVlCh8x9eFD7rpShfn6lhBCCCGCqlu3Lp49e5Zr34gRI+Dv7w/GGA4fPoyffvoJDRs2FCQ+odBUPiKopk2b4sKFCzh58iQ8PT0RHh6OSZMmYePGjVi9ejW6d+8OEd10E0IIUUSfp+kdOsSn6qWkfHnO1BTo25dvTZpQAoqQPJw4caLQ7+nQoQM0NTVlEA0hRNa0tIDkZOHalgbG2Hefb7t06YIxY8bg0qVL8PPzK3OjpQBKTBE5IBKJ0KVLF3Ts2BF+fn6YN28eoqKi0LNnT7Rt2xbr1q2DlZWV0GESQgghxUfJKEKkpmfPnoV6vUgkQlRUFGrUqCGbgAghMiUSAaW9HnhkZCTMzc1z7VNRUcGQIUOwYMEC3LhxAwEBAQJFJxyaykfkhoqKCkaPHo3Hjx9j9uzZUFdXx/nz52FjY4Nx48bh3bt3QodICCGEFN5/TdObOhW4fh149oyvuNe0KSWlCCmguLg4iMXiAm1a0hryQAghRXD+/Hncu3cPP//883fPjRgxAhcvXkSPHj1QsQzWjKTEFJE72traWLZsGR48eIC+fftCLBZj8+bNqF27NjZu3JhreU1CCJFHPj4+sLCwQOPGjYUOhQiFklGEyNywYcMKNS1v8ODB0NHRkWFEhBDCZWRkIC4uDq9evUJISAiWL1+OHj16oGvXrhg6dOh3r69Xrx7evXuHHTt2CBCt8GgqH5Fb1atXxx9//IGLFy/C3d0dd+/excSJE7Ft2zb4+PigZcuWQodICMkD1fzgq/KNHz8eiYmJ0NXVFTocUpJCQoDNm4H9+3MXwqBpeoRIXWE/wG3atElGkRBCSG7//PMPjIyMoKKigooVK8La2hobNmzAsGHDoKSU9/ggfX39Eo5SfogYY0zoIGSpV69eCAwMRLt27XD48OFcz6WmpqJevXro27cv1qxZU+Bjfv6gkZCQQN+6lJDs7Gxs27YNc+bMwcePHwHwb8lWrVqFypUrCxwdIYpBWn1bfhfb/ChyzQ+6XpQRKSnAwYM8IXXr1pf9lIwiCkhe+rW0tDR8+PABVatWzbU/PDwclpaWAkVVdPJyXgmRF+np6YiOjoa5uTk0NDSEDofk40f/ToXp1xR+Kp+7uzt27dqV53PLli1D06ZNSzgiUhQqKioYN24cHj16BFdXVwDAzp07UbduXWzZsgVisVjgCAkhX6OaH6RMCA8H3N2BqlWBkSN5UkpVlU/fCwykaXqEyMjhw4dRu3ZtdO7cGVZWVrhx44bkuSFDhggYGSGEkKJQ+MSUo6MjtLW1v9sfFRWFBw8eoHPnzgJERYqqUqVK2LZtG65fvw4bGxt8/PgRY8eORYsWLXD37l2hwyOEgGp+EAWXkQHs2we0bg3Urw94ewMJCUCNGsDKlcDLl/z5Nm0oGUWIjCxduhQhISG4c+eOZGn1ffv2AeBLsRNCCCld5DoxdenSJXTr1g3GxsYQiUQ4duzYd6/x9fWVDBuzs7PD5cuXC3TsadOmwcvLS8oRk5LStGlT3Lx5E7/99hu0tbVx/fp12NnZYebMmUhNTRU6PELKtB07duT5hUB+Nm3ahEqVKskwIkKkICoKmD4dqFYNGDQIuHwZUFYGevUCTp3iz8+YAdD0ckJkLisrCwYGBgCARo0a4dKlS9iyZQsWL14MESWECSGk1JHrxFRKSgqsra2xcePGPJ8/ePAgJk+ejDlz5iA0NBStWrWCs7MzYmJifnjc48ePo3bt2qhdu3aB4sjIyEBiYmKujQhPRUUF7u7uiIyMRO/evZGdnY2VK1eiQYMGOH/+vNDhEUIIKe2ysoAjR4AOHYDatfm0vHfveHJq8WIgJoavvOfkBBSythohpOgqV66ca6S8vr4+zpw5g8jISBpBTwghpZBc30U5Oztj6dKl6N27d57Pr127FiNHjoSrqyvq1auH9evXw8TE5D9X3Lh+/ToOHDiA6tWrY9q0adi2bRsWL16c7+u9vLygq6sr2UxMTIr1cxHpqlq1Ko4cOYLjx4+jWrVqePr0Kdq1a4eRI0dKCqUTQkpGWloaXr169d3+8PBwAaIRjo+PDywsLNC4cWOhQyFFERMDzJvHi5f36QOcPcun5XXuDJw4AURH8+eNjYWOlJAyaffu3d8tfqOmpob9+/fj4sWLAkVFCCGkqOQ6MfUjmZmZCA4OhpOTU679Tk5OCAoK+uF7vby88OLFCzx79gxr1qzBqFGjMH/+/HxfP2vWLCQkJEi2Fy9eSOVnINLVvXt3hIeHw83NDQDg5+cHS0vLIi1dTwgpPCpG+8X48eMRERGBW1+v0EbkX3g40L8/YG4OLF0KxMUBhobA7NnA06fA338D3boBKipCR0pImVatWjVUqVIlz+datGhRwtEQQggprkLfWRXlQ36HDh0KVQi3IN69e4ecnBwYGhrm2m9oaIi4uDjJ444dOyIkJAQpKSmoVq0aAgICCv0Ntrq6OtTV1aUSN5EtHR0d+Pj4YMCAARg5ciQePXqEHj16YMCAAdiwYQPVsSFEhj4XozUwMMDt27cxbNgwzJkzBwMHDqRitES+3b/Pp+YdPgx8/r/arh0wdizQvTugpiZsfISQH4qLi8s3UUUIIUT+FTox1bNnz0K9XiQSISoqCjVq1ChsUwU+/tcYY7n2nTp16ofvd3FxkUVYRGAtW7ZEWFgYFi1ahNWrV2P//v04d+4ctmzZUuj/w4SQgsmrGG3v3r3x+PFjKkZL5NO9e18SUp/16cOn6VlZCRcXIaRQnJycqLYUIYSUYkWayhcXFwexWFygTUtLS9oxAwAqVaoEZWXlXKOjACA+Pv67UVSkbNLU1MSKFStw/fp1WFpaIj4+Hr169cLgwYPx4cMHocMjROFQMVpSaty9yxNQVlZfklJ9+/L9hw5RUoqQUoZG5RJCSOlW6MTUsGHDCjUtb/DgwdDR0SlsM/9JTU0NdnZ2OHPmTK79Z86cQfPmzaXeHim9GjdujODgYMycORNKSkrYu3cvGjRogNOnTwsdGiEKhYrRErl35w7w88+AtTVfbU8kAvr14yOn/vgDaNBA6AgJIUVAo3IJIWWVi4uLQswIKnRiaseOHdDW1i7w6zdt2lTkuj7JyckICwtDWFgYACA6OhphYWGIiYkBAHh4eOD333+Hn58fIiMjMWXKFMTExGDs2LFFao8oLnV1dXh5eeHq1auoXbs2Xr9+jY4dO8LNzQ0pKSlCh0eIQqBitERuhYUBvXoBDRsCR4/yhFT//jwhdfAgUL++0BESQgghREE4ODhg8uTJMn+PIpHrZWVu374NR0dHyWMPDw8AfNSWv78/+vfvj/fv32Px4sWIjY1F/fr1cfLkSZiZmQkVMpFzzZo1Q2hoKGbOnAlvb29s2rQJ586dw549e2hZd0Jk5Ny5czh37hzi4+MhFotzPefn5ydQVKRMCA0FFi0Cjh/njz8npObNAywshI2NEFJqXLp0CatXr0ZwcDBiY2MREBCQa4QCYwyLFi3C1q1b8fHjRzRt2hQ+Pj6wtLQULmhCCClFCjVi6uPHj5LaPG/fvsWRI0dw//59mQQG8KwhY+y7zd/fX/IaNzc3PHv2DBkZGQgODkbr1q1lFg9RDFpaWtiwYQPOnDmDqlWr4tGjR2jevDmWLFmC7OxsocMjRKEsWrQITk5OOHfuHN69e4ePHz/m2hSVj48PLCwsKOEtlPBwoGdPwNaWJ6VEImDAAL763v79lJQiRMGoyXjlzJSUFFhbW2Pjxo15Pr9q1SqsXbsWGzduxK1bt1ClShV06NABSUlJMo2LkLKGMYaUlBRBtoLWsnNxccHFixfx22+/QSQSQSQS4dmzZ7h48SKaNGkCdXV1GBkZYebMmZLPnvm9JycnByNHjoS5uTk0NTVRp04d/Pbbb7I8xcJhBbRt2zZWo0YNVr16debr68tatGjBxo0bxywtLdmWLVsKehiFkJCQwACwhIQEoUMhxfT+/XvWr18/BoABYC1atGDR0dFCh0WIIGTRt1WpUoXt2rVLascrbeh6UcKSkxmbMYMxFRXGAMZEIsYGDmQsIkLoyAhRGGW9XwPAAgICJI/FYjGrUqUKW7FihWRfeno609XVZZs3b873OOnp6SwhIUGyvXjxokyfV0K+lZaWxiIiIlhaWppkX3JysuRzW0lvycnJBYr706dPzN7eno0aNYrFxsay2NhY9vLlS6alpcXc3NxYZGQkCwgIYJUqVWILFizI9z3Z2dksMzOTzZ8/n928eZM9ffqU7dmzh2lpabGDBw9K2hs2bBjr0aOHNE99oeT17/RZYa4XBZ7K5+3tjfDwcKSmpsLU1BTR0dEwMDBAYmIiWrdujdGjRxcpMUaIkPT09HDgwAF069YNbm5uuHr1KqytrbF582YMGDBA6PAIKfUyMzNpQQpSMo4fB9zdgX/rUKJHD2DFCqBuXWHjIoSUCC8vLxgaGmLEiBG59vv5+eHt27fw9PSUSbvR0dGIi4uDk5OTZJ+6ujratGmDoKAgjBkzJt94Fy1aJJOYCCHC0dXVhZqaGrS0tCS1V+fMmQMTExNs3LgRIpEIdevWxevXr+Hp6Yn58+fn+R4AUFZWztVPmJubIygoCH/88Qf69etX4j+bLBV4Kp+ysjI0NDSgp6eHmjVrwsDAAACgo6NDK2GQUk0kEmHw4MG4c+cO7O3tkZiYiIEDB2L48OFITk4WOjxCSjVXV1fs27dP6DCIInv+HOjenU/di4kBzMyAEyeAY8coKUVIGbJlyxbUzeN33tLSEps3b5ZZu3FxcQAAQ0PDXPsNDQ0lz+Vl1qxZSEhIkGwvXryQWYyEKAotLS0kJycLsmlpaRU57sjISNjb2+fKm7Ro0QLJycl4+fLlD9+7efNmNGrUCAYGBihfvjy2bdsmWQxOkRR4xJSKigrS09OhoaGRa9lvmjtNFIW5uTkuXbqEJUuWYOnSpfD398fVq1dx4MAB2NraCh0eIaVSeno6tm7dirNnz8LKygqqqqq5nl+7dq1AkZFSLzMTWLsWWLwYSEsDVFSAadOAuXOBcuWEjo4QUsLi4uJgZGT03X4DAwPExsbKvP1vv6hnjP3wy3t1dXWoq6vLOixCFIpIJEK5UniNz6s/YP/WrPpRP/HHH39gypQp+PXXX2Fvbw9tbW2sXr0aN27ckGm8QihwYur8+fOSzlNXV1eyPy0tDdu3b5d+ZIQIQEVFBYsWLULbtm0xePBgREVFwd7eHqtXr8bEiRNpdCAhhXT37l00bNgQAL5bLIN+n0iRXbwIuLkBERH8cZs2gK8vFTUnpAwzMTHB1atXYW5unmv/1atXYWxsLLN2P0+7+TYxFh8f/90oKkJI2aCmpoacnBzJYwsLCxw5ciRXgiooKAja2tqoWrVqnu8BgMuXL6N58+Zwc3OT7Hvy5EkJ/AQlr8CJqfLly+d6HBcXhypVqqBy5cqoXLmy1AMjREht2rTBnTt3MGLECBw/fhyTJk3C+fPn4efnBz09PaHDI6TUuHDhgtAhEEUSHw9Mnw7s2sUfGxgAa9YAQ4bwlfcIIWWWq6srJk+ejKysLLRt2xYAcO7cOcyYMQNTp06VWbvm5uaoUqUKzpw5AxsbGwC8vuLFixexcuVKmbVLCJFf1atXx40bN/Ds2TOUL18ebm5uWL9+PSZOnIgJEybg4cOHWLBgATw8PKCkpJTnez6XUNq1axdOnToFc3Nz7N69G7du3fouAa8IClxj6ltfF/gjBRQYyL/RPXwYuHQJePAA+PABKODSk6Rk6enpISAgABs2bICamhqOHz8OGxsbXL9+XejQCCGkbBGLga1bec2oXbt4EmrMGH4dHTqUklKEEMyYMQMjR46Em5sbatSogRo1amDixIlwd3fHrFmzinXs5ORkhIWFISwsDAAveB4WFoaYmBiIRCJMnjwZy5cvR0BAAO7fvw8XFxdoaWlh4MCBUvjJCCGlzbRp06CsrAwLCwsYGBggKysLJ0+exM2bN2FtbY2xY8di5MiRmDt3br7viYmJwdixY9G7d2/0798fTZs2xfv373ONnlIkIsaKlhVp0KAB7t27J+14SoXExETo6uoiISEBOjo6BX+jmxuwadP3+1VVgcqVAUNDvlWp8v1mZMS3b0aukZIREhKC/v374/Hjx1BRUcGKFSvg4eFBU5GIQily3/YNDw8PLFmyBOXKlYOHh8cPX6voNaakdU7LtPv3gVGjgM9fCjRsyK+lzZoJGhYhZZW892vJycmIjIyEpqYmatWqJZU6ToGBgXB0dPxu/7Bhw+Dv7w/GGBYtWoQtW7bg48ePaNq0KXx8fFC/fv0CtyHv55WQkpaeno7o6GiYm5tDQ0ND6HBIPn7071SYfq3AU/m+RR/Ii8DGBujdm09FiI8H3rwBEhKArCzg1Su+/RdtbZ6gMjb+slWtyrdq1fhmZMSLwBKpsbW1RXBwMFxdXXHo0CFMmzYNly5dgr+/PypWrCh0eITIldDQUGRlZUn+nh+6jpAfYgzw8wMmTADS0/kXM0uW8Md0jSOE5KN8+fJo1KgRAOldZxwcHPCj7/JFIhEWLlyIhQsXSqU9Qggpa+jOriSNGsW3r2VkfElSvXkDxMV9+fPzFhvLt5QUICmJb48e5d+OkhIfZWViwjdT0y+bmRlQvTpQsSJNfSgkHR0dHDx4EA4ODpgyZQpOnDgBOzs7HDp0CHZ2dkKHR4jc+LquFNWYIkWSksJHGX+uJeXsDGzbxr+EIYSQfGzfvh3r1q1DVFQUAKBWrVqYPHkyXF1dBY6MEELIj1BiSmjq6l8SSP8lKYknqF6/5n9+HmX17ZadzV/z+jWQ31KS5cvzBJW5OVCjxvcbDZfMk0gkgpubG5o2bYq+ffsiOjoazZs3x4YNGzB69GgaAUJIHtLT03H37l3Ex8dDLBZL9otEInTr1k3AyGTHx8cHPj4+362uQgogMhLo04evuKekBCxbBsyYwf9OCCH5mDdvHtatW4eJEyfC3t4eAHDt2jVMmTIFz549w9KlSwWOkBBCSH6KXGOqUaNGuH37trTjKRXkeg64WMxHXL18ybcXL4CYGOD5c77FxPDnf0Qk4omymjWBWrX4Vrs2Lzprbk5TKP718eNHuLi44MSJEwCAoUOHYtOmTdDS0hI4MkKKRhZ92z///IMhQ4bg/fv33z0nEokUPnEj19cLebRnDy9qnprKp6Xv3w+0aSN0VISQr8hrv1apUiV4e3tjwIABufbv378fEydOxLt37wSKrGDk9bwSIhSqMVU6CF5jqqwmpeSektKXQumNG+f9mrQ0nqCKjubb06dftidP+MismBi+nT+f+72qqjxRZWEB1KsHWFoC9evzxJWqqux/PjlSsWJFBAQEYPXq1Zg9ezZ27dqF0NBQHD16FDVr1hQ6PELkwoQJE9CvXz/Mnz8fhoaGQodD5FVaGjBpEp+uBwDt2gF79/IFQQghpABycnIktaW+Zmdnh+zsbAEiIoRIw9ej7Yn8kda/Dw19KYs0NYE6dfj2LcaAt295gurRIyAqiv/56BHw8CEvQBsRwbevqaryEVUNGgBWVnzVJBsbvtqgAlNSUoKnpyeaNm2KX375Bffu3UOjRo2wa9cudO/eXejwCBFcfHw8PDw8KClF8hcVBfTtC9y5w0fszp8PzJsHKCsLHRkhpBQZPHgwNm3a9N1qr1u3bsWgQYMEiooQUlRqampQUlLC69evYWBgADU1NSqbIkcYY8jMzMTbt2+hpKQENTW1Yh2vyFP5AMDLywuGhoYYMWJErv1+fn54+/YtPD09ixWcvCqzQ23FYj41MDKSb+HhPEF1/z4fZZUXY2PA1hawswMaNeKjuBT0A+rr16/Rt29fBAUFAQBmz56NxYsXQ5k+XJFSQhZ924gRI9CiRQuMHDlSKscrbcrs9aKgDh0CRo7k1xADAz5KqkMHoaMihPyAvPZrEydOxK5du2BiYoJmzZoBAK5fv44XL15g6NChUP1qZP+3ySt5IK/nlRAhZWZmIjY2FqmpqUKHQvKhpaUFIyOjPBNThenXipWYql69Ovbt24fmzZvn2n/jxg388ssviI6OLuqh5RpdOL7BGJ/2d+8ecPcu/9Y7LIx/C57Xfy8zM6BJE8Denm82NrwIvALIzMzE9OnTsWHDBgCAk5MT9u3bB319fYEjI+S/yaJvS01NRd++fWFgYIAGDRrk+mAAAO7u7lJpR17R9SIfGRnAtGnAxo38ccuWwIEDtOoeIaWAvPZrjo6OBXqdSCTC+W9LVcgBeT2vhAiNMYbs7GyFr0taGikrK0NFRSXfkWwllpjS0NBAZGQkzM3Nc+1/+vQpLCwskJ6eXtRDyzW6cBRQUhJPUgUH8+32beDBg++TVRoaPFHVsiXQqhXQogWgrS1MzFKyf/9+jBw5EmlpaahevTqOHj0KGxsbocMi5Idk0bf9/vvvGDt2LDQ1NaGvr5/rwiUSifD06VOptCOv6HqRh7g4oHt34NYt/tjTE1i6lBbWIKSUoH5NNui8EkIUTYkUPwcAExMTXL169bvE1NWrV2FsbFycQxNFoK3Nk00tW37Zl5DAk1TXr/Pt2jXg3Tvg0iW+AbyuSKNGgKMjL4DbogWvi1WKDBgwAJaWlujduzeePHmC5s2bY+vWrRgyZIjQoRFSoubOnYvFixdj5syZUFJSEjocIrSYGKB9ez6itmJFYNcuoGtXoaMihBBCCCECKlZiytXVFZMnT0ZWVhbatm0LADh37hxmzJiBqVOnSiVAomB0dYG2bfkG8NFTjx4BV64Aly8DFy8Cz54BN27wbcUKPs2vZUugY0egUye+CmApKHxnZWWFW7duYfDgwTh58iSGDh2K27dvY82aNd9NZyJEUWVmZqJ///6UlCI8GdW+PU9OmZoCZ8/yVV4JIYQQQkiZVqypfIwxzJw5Exs2bEBmZiYAPr3P09MT8+fPl1qQ8oaG2srY8+dAYCBw7hzfXr/O/Xy1akDnzvxb9nbtAC0tQcIsKLFYjIULF2LJkiUAgNatW+PQoUOorOArFpLSRxZ925QpU2BgYIDZs2dL5XilDV0v/nXvHi9q/uYNULs2T0qZmAgdFSGkCKhfkw06r4QQRVNiNaY+S05ORmRkJDQ1NVGrVi2oK0gh6/zQhaMEMcbrUp0+DZw6xRNWaWlfntfUBJycgB49eM0SOS4yfuzYMQwdOhRJSUmoVq0aAgIC0KhRI6HDIkRCFn2bu7s7du3aBWtra1hZWX03WlAeV0aSJrpeALh5k492/fgRsLLi/bmCrs5KSFlA/Zps0HklhCiaEk9MAXz0FIB8K7IrErpwCCgtjSen/v4b+PNPPiXkM2VlXpeqb1+gd2+gUiXBwsxPZGQkevbsiUePHkFdXR1btmzBsGHDhA6LEACy6dt+tEqSvK6MJE1l/npx8SIf3ZqcDDRtCpw8CejpCR0VIaQYyny/JiN0XgkhiqZEE1Pbt2/HunXrEBUVBQCoVasWJk+eDFdX1+IcVq7RhUNOMMZX/Tt+HDh6FLh798tzysq8lsmAATxJJUer/CUkJGDw4MH466+/APARJVR3isgD6tukr0yf0//9j/e/6en8S4Pjx+WqLyaEFE2Z7tdkiM4rIUTRFKZfK1Y12nnz5mHSpEno1q0bDh06hEOHDqFbt26YMmUK5s6dW5xDE/LfRCKgYUNgwQKeoIqK4sXSbW2BnBw+9c/FhU8ZGTCAf0jKyRE6aujq6uL48eNYsGABAGDDhg3o0KED3r59K3BkhBAiJYcP8ynW6elAly58lCslpQghhBBCSB6KNWKqUqVK8Pb2xoABA3Lt379/PyZOnIh3794VO0B5RN9olAJRUcCBA8CePXzVv8+MjYGhQ4ERI+RiNahjx45hyJAhSE5OhqmpKY4dOwYbGxuhwyJllCz6Ni8vLxgaGmLEiBG59vv5+eHt27fw9PSUSjvyqkxeL3bu5H2sWAz06wfs3g2oqQkdFSFESkpjv/bkyRMAwE8//SRwJPkrjeeVEEJ+pMRGTOXk5ORZvNnOzg7Z2dnFOTQhxVOrFjBvHi+cfusWMHEir2vy+jUfVVW7NtC2LU9e/buipBB69uyJGzduoFatWoiJiUGLFi2wb98+weIhRNq2bNmCunXrfrff0tISmzdvFiCikuHj4wMLCws0btxY6FBKlo8PH6kqFvPk1L59lJQihAjmw4cP6Nq1K7y9vbFhwwZ06dIFHz9+FDosQggh3yhWYmrw4MHYtGnTd/u3bt2KQYMGFefQhEiHSAQ0agRs2MCTUocPA87OfP+FC3yKn4kJMGdO7kLqJcjCwgI3b96Es7Mz0tLSMGjQIEyfPp2Su0QhxMXFwcjI6Lv9BgYGiI2NFSCikjF+/HhERETg1q1bQodScnbtAiZM4H+fNAnYto3X+yOEEIHMmjULHh4eWL9+PX777Td4eHhg2rRpQodFCCHkG8VKTAG8+Hn9+vXh6uoKV1dX1K9fH9u2bYOSkhI8PDwkGyGCU1cHfv6Zrwr17BmvTWVsDMTHA8uXA+bmQJ8+wOXLvLB6CapQoQL+/PNPzJw5EwCwZs0adO7cGR8+fCjROAiRNhMTE1y9evW7/VevXoWxsbEAERGZCAoCRo3if582DVi3DlAq9i0GIYQUS3h4ONq2bYvFixfjw4cPaNeuHSIjI4UOixBCyDeKddd4//592NrawsDAAE+ePMGTJ09gYGAAW1tb3L9/H6GhoQgNDUVYWJiUwiVESkxNgYULeYLqyBE+rU8s5n9v3ZqPstq/H8jKKrGQlJWV4eXlhYMHD0JLSwtnzpxB48aNce/evRKLgRBpc3V1xeTJk7Fjxw48f/4cz58/h5+fH6ZMmYJRnxMZpHR7/hzo2ZNPi+7VC1i5ko9KJYQQOfF1SV2xWCz142dnZ2Pu3LkwNzeHpqYmatSogcWLF8ukLUIIUUQqxXnzhQsXpBUHIcJQVeXLmffuDdy/z6f87d4NhIQAAwcCM2cCU6bwkQDlypVISP369UOdOnXQs2dPPH36FPb29ti5cyd+/vnnEmmfEGmaMWMGPnz4ADc3N2T+W89NQ0MDnp6emDVrlsDRkWJLTga6dwfevuWrpO7eTSOlCCFyo3nz5ti3b59kJeT9+/ejRYsWUm9n5cqV2Lx5M3bu3AlLS0vcvn0bw4cPh66uLiZNmiT19gghpCRkZmbi1q1buHv3Lu7du4f79+/jyZMnEIvFUFJSgpKSErS0tGBhYYEGDRrAysoKVlZWqFmzJpQKeT9YrFX5yipaNUPBvXsHbNoEeHvzD1sAoK/Pa6ZMmABUrFgiYbx//x79+/fHuXPnAABz5szBokWLoEw1W4iMyLJvS05ORmRkJDQ1NVGrVi2oq6tL9fjySqGvF2IxHyF14gRgaMgXmjAxEToqQoiMlaZ+LS0tDRMmTMCnT58AALq6uvD19YWGhoZU2+natSsMDQ2xfft2yb6ff/4ZWlpa2L17d4GOUZrOKyFEMTHGEBkZidOnT+PMmTMIDAxEampqoY+jra2Nhg0bShY7Kki/VqTE1LfLfufHz8+vsIeWul69eiEwMBDt2rXD4cOHAQAvXrzAkCFDEB8fDxUVFcybNw99+/Yt8DHpwlFGpKXxYr6rVwP/LjMMHR2+wt+UKTxZJWPZ2dmYOXMmfv31VwBA586dsXfvXlSoUEHmbZOyR9p9W1ZWFpycnLBlyxbUrl1bChGWPgp9vZg5k0/bU1cHAgOBZs2EjogQUgJKY7+Wnp4Oxhg0NTVlcvwVK1Zg8+bNOH36NGrXro07d+7AyckJ69evx4ABA/J8T0ZGBjIyMiSPExMTYWJiUqrOKyGk9GOM4c6dOzh8+DAOHz6Mhw8f5nrewMAAjRs3Rv369dGgQQPUqVMHampqEIvFEIvFSEhIwP3793H37l3JyKr09PRcx5BZYkpJSQlmZmawsbHBj94eEBBQ2ENL3YULF5CcnIydO3dKElOxsbF48+YNGjZsiPj4eNja2uLhw4coV8CpWqXxgkyKITubr+a3fDnwud5T+fI8QTVtGqCnJ/MQ9u7dC1dXV6Snp6NmzZo4duwYLC0tZd4uKVtk0bcZGBggKCgItWrVksrxShuFvV7s2gUMG8b/vncvn/pMCCkTFLZfKwbGGGbPno2VK1dCWVkZOTk5WLZs2Q+nrC9cuBCLFi36bj+dV0JISXj37h38/f2xbds2PHr0SLJfTU0Nbdq0gZOTEzp06IAGDRoUalpednY2Hjx4gJCQEAQFBWHLli2yS0y5ubnhwIEDMDU1xYgRIzB48GDolcCH86IKDAzExo0bJYmpb1lZWeHvv/+GSQGnINAFuYwSi/mUlcWLgdBQvk9Hh4+e8vDgf5eh0NBQ9OrVC8+fP0e5cuXg7++PPn36yLRNUrbIom+bOnUqVFVVsWLFCqkcr7RRyOtFUBDg6MiLnc+eDSxbJnREhJASpJD9WjEdOHAA06dPx+rVq2FpaYmwsDBMnjwZa9euxbDPSfxv0IgpQogQrl27JsmNfF3/1dnZGX369EHXrl2l1gcV5npR5BpTGRkZOHr0KPz8/BAUFIQuXbpg5MiRcHJygkhKq/FcunQJq1evRnBwMGJjYxEQEICePXvmeo2vry9Wr16N2NhYWFpaYv369WjVqlWu1/woMXX79m24uLjg/v37BY6LLshlHGPAn38C8+cDd+7wffr6/AOamxsg5boFX3v37h369++P8+fPAwA8PT2xbNkyqjtFpEIWfdvEiROxa9cu1KxZE40aNfpuZOratWul0o68UrjrxfPnQJMmQHw8ry91+DAVOyekjJG3fm3hwoWwtbWFnZ0dqlatKkgMJiYmmDlzJsaPHy/Zt3TpUuzZswcPHjwo0DHk7bwqkrdv3yI6OhovX77Ey5cv8erVK6SkpCAjIwPp6enIysqClpYWtLW1oa2tjQoVKsDMzAzVq1eHubk59PX1pfb5lhChXLp0CYsXL5bULwYAW1tbjBkzBgMGDIC2trbU2yxMv1bkVfnU1dUxYMAADBgwAM+fP4e/vz/c3NyQlZWFiIgIlC9fvqiHlkhJSYG1tTWGDx+e54pkBw8exOTJk+Hr64sWLVpgy5YtcHZ2RkREBExNTf/z+O/fv8fQoUPx+++///B1eX2jQcowkYivQtW1K3D0KDBvHvDgATB1KrB+PbBkCTB4MCCDZFGlSpVw6tQpzJo1C2vWrMHKlSsRHByM/fv3o1KlSlJvj5Diun//PmxtbQEg1zBhAHSTV9p8XoEvPp5W4COEyI3FixdLrieVKlWCnZ0dbG1tJckqMzMzmceQmpr63VQXZWVliMVimbdNcktOTsbFixdx48YNhIaGIiQkBK9fvy7WMXV0dGBjYwM7OzvY2dmhSZMm+Omnn+g+hpQKV65cwdy5c3Hx4kUAgIqKCgYPHozx48ejUaNGAkf3RZETU18TiUQQiURgjEm1A3Z2doazs3O+z69duxYjR46Eq6srAGD9+vU4deoUNm3aBC8vrx8eOyMjA7169cKsWbPQvHnzH77Wy8srzzngpIxTUgL69AF69uT1VhYuBF68AFxcgLVredF0JyepN6uiooLVq1ejUaNGGDFiBM6ePQs7OzscPXoUdnZ2Um+PkOK4cOGC0CEQaWAMGDoUuHuXr8B34gRQwLqMhBAiS40bN0ZsbCyGDx+OKlWqICQkBCdPnsTq1auRnZ2NihUrwtbWFqdPn5ZZDN26dcOyZctgamoKS0tLhIaGYu3atQVeMIoUT0REBE6cOIHTp0/jypUryMrKyvW8SCRC1apVUa1aNVSrVg1Vq1aFjo4O1NXVoaGhARUVFaSmpiIpKQlJSUn48OEDnj9/jmfPnuH169dITEzExYsXJR/sAcDMzAzt27eXbPQFMZE3sbGxmD59Ovbu3QsAUFVVxciRIzFz5swSSdgXllSm8l25cgVdu3bF8OHD0alTp0IVxyookUiUaypfZmYmtLS0cOjQIfTq1UvyukmTJiEsLCxXx/HtVD7GGAYOHIg6depg4cKF/9k2zQEnBZKeDnh783orCQl8n7Mz8OuvQL16Mmny/v376N27N6KioqCurg4fHx+MHDlSJm0RxSfLaQQRERGIiYmRzGUHeL/erVs3qbYjbxRmasbnYudqasDFi7QCHyFlmDz2a/7+/pg9ezZsbGywbt061K5dG1lZWbh79y5CQkIQGhoKX19fmbWflJSEefPmISAgAPHx8TA2NsaAAQMwf/58qKmpFegY8nhe5VlsbCz279+PPXv2IPRz7dd/mZubw8HBAba2trCxsYG1tXWRZ/NkZGTg4cOHCAkJQXBwMIKDg3H79u1cyS9lZWU4ODjg559/Rq9evVClSpVi/WyEFEdWVha8vb2xcOFCJCUlQSQSwdXVFfPmzStwTW1pKVS/xopg3LhxrGLFisza2pqtX7+evXv3riiHKRQALCAgQPL41atXDAC7evVqrtctW7aM1a5dW/LYycmJVapUiWlqarKqVauymzdvssuXLzORSMSsra0l2927dwscS0JCAgPAEhISiv1zEQX07h1jU6YwpqrKGMCYsjJj48cz9v69TJr79OkT6969OwPAALCRI0ey1NRUmbRFFJss+rYnT54wKysrJhKJmJKSEhOJRJK/KykpSa0deaUQ14vYWMYqVuT9mZeX0NEQQgQmr/1aUlISmz59OtPW1mbTpk1jSUlJQodUKPJ6XuWJWCxm58+fZ927d2dKSkqSe18VFRXWuXNn5u3tzR49esTEYrFM40hOTmb/+9//2LRp05iVlZUkDgBMJBKx9u3bs3379tH9OClxISEhzNraWvL/sWnTpuzWrVuCxVOYfq1II6aUlJRgamoKGxubH86tPXr0aGEPna9vR0y9fv0aVatWRVBQEOzt7SWvW7ZsGXbv3l3gQoNFQd9okAKJigI8PYGAAP5YT4/Xnxo9GlCRyixaCbFYjBUrVmDevHkQi8WwsbHB4cOHUaNGDam2QxSbLPq2bt26QVlZGdu2bUONGjVw8+ZNvH//HlOnTsWaNWu+W6xC0SjE9aJvX17k3NYWuHFD6v0XIaR0kfd+LSoqCh4eHrh16xaWL19eaqbTyft5FVJmZib279+P9evXIywsTLK/efPmGDx4MPr27SvoVLqnT5/iyJEjOHz4MG7evCnZX6FCBQwcOBBjxoyBlZWVYPERxZeRkYElS5ZgxYoVyMnJgb6+PlauXInhw4fLZDZbQRWmXytSlEOHDoWjoyMqVKgAXV3dfDdZqlSpEpSVlREXF5drf3x8PAwNDWXaNiEFUqsWL45+/jzQoAHw4QMwfjxgZwdcvizVppSUlDB79mycOnUKlSpVQmhoKGxtbXHs2DGptkNIYV27dg2LFy+GgYEBlJSUoKSkhJYtW8LLywvu7u5Ch0f+y9GjPCmlogJs305JKUKIXMvKykJaWhp++eUXmJqaYtSoUfjw4YPQYZEiys7Ohr+/P+rUqQMXFxeEhYVBS0sL48aNw4MHD3D16lWMGzdO8PpONWrUwPTp03Hjxg08ffoUCxYsgKmpKT59+gRfX19YW1vD0dERx44dQ05OjqCxEsUTEhICW1tbLFu2DDk5Oejbty/Cw8MxcuRIQZNShVWkO0x/f38ph1F4ampqsLOzw5kzZ3LVmDpz5gx69OghYGSEfMPREQgJAbZs4Sv43b0LtG7NV+5bvRqQ4jz09u3bIyQkBP3798e1a9fQq1cvTJ06FV5eXlBVVZVaO4QUVE5OjqSuQ6VKlfD69WvUqVMHZmZmePjwocDRkR/6+BFwc+N/9/TkK/ERQoicWbZsGe7du4d79+7h0aNHKFeuHKysrNC0aVOMGTNG5l+WE+ljjOHw4cOYP3++ZBZMlSpVMGnSJIwePRp6enrfvSclBXj4kC+UHRkJvH4NvH/Pvxf+8AHIzuYLZn/e9PT4LXiVKoCREVCnDmBpCVSvXrwFZ83NzbFw4ULMnz8f58+fx9atW3H06FEEBgYiMDAQ5ubm8PDwwMiRI6GpqVn0hkiZxxjD+vXr4enpiaysLFSuXBm+vr74+eefhQ6taGQ8rbBYkpKSWGhoKAsNDWUA2Nq1a1loaCh7/vw5Y4yxAwcOMFVVVbZ9+3YWERHBJk+ezMqVK8eePXsm07hoDjgpsrdvGRs1ijGRiNdr0dVlbONGxrKzpdpMZmYmmzJlimR+sb29veT3hpD8yKJva9mypaQ+4IABA1inTp3YlStX2NChQ5mlpaXU2pFXpfp64eLC+6m6dRlLSxM6GkKInJC3fk0kEjFzc3O2aNEi9ujRI6HDKTJ5O69CCQ0NZS1atJDcw+rp6bFVq1axlJSUXK+LiWFszx7GRo9mrE4dfrmSxqalxViTJrxc7LFjvHRsccXExDBPT0+mp6cn+bkqV67MVqxYUeb/vUnRxMfHs86dO0v+P/Xs2bNE6n4XlsxrTJWUwMBAODo6frd/2LBhklFbvr6+WLVqFWJjY1G/fn2sW7cOrVu3lmlcNAecFNvNm3wkQnAwf9yoEbB5M5/mJ0UBAQEYPnw4EhISULFiRfj7+6N79+5SbYMoDln0badOnUJKSgp69+6Np0+fomvXrnjw4AH09fVx8OBBtG3bVirtyKtSe704fRro2BEQiYCrV4GvajkSQso2eevXWrdujTt37iApKQmampqwsrKCnZ0dbG1tYWdnh/r160NZWVnoMP+TvJ3Xkvbx40fMmzcPmzZtglgshpaWFqZPn44pU6ZAV1cXjAEREXx2+eHDwP373x/DwACoW5dvZmaAvj7f9PQAVVUgJ4dv2dl8NFVcHN9evuTHfvAA+GrxYAkrK6B7d6BnT15u8Qclln8oNTUVO3bswOrVq/H8+XMAvA6Vh4cHJk2aVCb/3UnhBQYGYuDAgYiNjYW6ujrWrl2LcePG/bD2t1AK06/JdWJKXpX1CweRkpwcnoyaMwdISODjht3deYH0Ii5pm5fo6Gj0798ft27dAgC4u7tj5cqV0NDQkFobRDGUVN/24cMHVKxYUS4voNJWKq8XyclA/frA8+fApEnA+vVCR0QIkSPy2q9FRUUhODgYISEhCA4ORmhoKD59+gR1dXU0aNAgV1FqeSSv57Uk/PHHH5gwYQLevn0LAOjfvz/WrFmDatWq4e1bwM8P8PfniaPPlJV5kqh1a6BNG/79SXFLTWVnA48f8wocly7xLTIy92uqVQN69QKGDSt6kiorKwv79++Hl5eXZKqinp4eZsyYgQkTJqBcuXLF+0GIQsrJyYGXlxcWLFgAsViMevXq4cCBA3JdWJ8SUzJWli8cRAbi4gAPD2D/fv7YxATw9QW6dpVaE5mZmZg1axbWrl0LALC2tsaBAwdQt25dqbVBSj/q26SvVJ5Td3fA25sX2rh3T6qJckJI6Vea+rXo6Gjcvn0boaGhWL58udDh/FBpOq/SEh8fj/Hjx+Pw4cMAAAsLC3h7e8PRsS2uXeO3w4cOfRnFpKbGB/P26QN06wZUrFgSMfJBxMeOAf/8w2tZfWZpCbi4AEOGAEVZe0ssFuOPP/7AwoULJXU3DQ0NMX/+fIwaNYrqwxKJ+Ph4DB48GGfOnAEAuLi4YOPGjXKfxCxUvybTSYUKiuaAE5n43/8YMzf/Msm9Xz/GYmOl2sTff//NDAwMGACmqanJNm/ezMRisVTbIKUX9W15+/PPP1nt2rVZzZo12bZt2wr13lJ3Tq9c+VID78wZoaMhhMihUtevlRJl7bwePHiQ6evrMwBMRUWFzZ8/n6WnZ7DTpxlr2TJ33afGjRnbvp0xoU9NWhpjf/7J2C+/MKau/iU+VVXGBg1i7MaNoh03KyuL7dy5k9WoUUNSM6hWrVrs0KFDdJ9O2MWLF5mRkZHk89uOHTuEDqnACtOvlZ71AwlRdJ068Qnz06fz8cl//AHUq8fHLktpYGPnzp1x584dtGvXDmlpaRg7dix69uwpGTpNCMktOzsbHh4eOH/+PEJCQrBy5UrFXXo8PR0YOZL3NyNGAO3bCx0RIYQQBZOcnIzhw4ejf//+eP/+PaysrHDjxk3Y2y+Co6ManJyAK1cAdXVg+HDg1i1emnXECEDogWQaGnxCw/79fMLDli1A06ZAVhawdy//e7Nm/BY+J6fgx1VRUcHQoUMRGRmJjRs3wsDAAFFRUejbty/s7e1x7do12f1QRG4xxrBq1Sq0bdsWsbGxsLCwwK1bt+Di4iJ0aDIhs8TUkydP8OTJE1kdnhDFpKUFrFrFr8K2tsCnT/yq3LEj8OyZVJowMjLC6dOnsWbNGqipqeHEiRNo0KAB/vrrL6kcnxBFcvPmTVhaWqJq1arQ1tZG586dcerUKaHDkg0vL77WtpER8OuvQkdDCCFEwYSFhaFRo0bw9/eHkpIS5s6dix07bsHT0wbOzsC1azz5M2kS8PQpry3VqJHQUeetQgVg9Gjg+nV+2z50KJ9qeOMG0L8/YGEB7NiRdzH1/KipqWH8+PF48uQJ5s+fDy0tLdy4cQPNmzfHL7/8IimYThTfx48f0bNnT3h6eiInJweDBw+W3JMqKqknpj58+ICuXbvC29sbGzZsQJcuXfDx40dpN0OIYrOx4Ve2Vav4FfrMGV6MeONGQCwu9uGVlJQwdepU3Lx5ExYWFnjz5g26deuGUaNGISkpSQo/ACHy4dKlS+jWrRuMjY0hEolw7Nix717j6+sLc3NzaGhowM7ODpcvX5Y89/r1a1StWlXyuFq1anj16lVJhF6y3r79koz67Td+x00IIYRIAWMMPj4+aNq0KR4+fIiqVavi6NHzePduCRo3VsPZszyp4+EBREfzNTeMjYWOuuAaNQJ27gRiYoAFC3jtq0eP+CivmjX5WkdZWQU/nra2NhYtWoTHjx9jxIgREIlEOHjwIOrUqYM5c+YgOTlZdj8MEVxISAjs7Oxw4sQJqKurY8uWLdi1a5fc15MqLqknpmbNmgUPDw+sX78ev/32Gzw8PDBt2jRpN0OI4lNR4dP67t7lS46kpAATJ/KlRx49kkoT1tbWCA4OhoeHB0QiEX7//XdYWVnh/PnzUjk+IUJLSUmBtbU1Nm7cmOfzBw8exOTJkzFnzhyEhoaiVatWcHZ2RkxMDAB+M/2tH60mmJGRgcTExFxbqbBqFe9j7Ox4VVlCCCFEClJTUzF06FBMmDABmZmZ6NatG6ZODYOLSxts3sy/b+3bl6+49+uvQJUqQkdcdIaGwMKFfFHbVav44xcvgHHjgLp1gT17CjfFz8jICNu3b0dwcDAcHByQkZGB5cuXo06dOtizZ0+e9yik9GKMYevWrWjevDmio6Nhbm6OoKAgjB49ukysZC314uctWrRgjDG2aNEi9v79e8YYY/b29tJuRlBlrTghkQM5OYz5+DBWrhyvsqihwdivvzKWnS21Ji5cuMDMzMwkRRfd3NxYUlKS1I5P5J8s+zZ7e3tWuXJlqR+3MACwgICAXPuaNGnCxo4dm2tf3bp12cyZMxljjF29epX17NlT8py7uzvbu3dvvm0sWLBA8jv09SbX14vYWMY0NXnf8vffQkdDCJFz8nofvHz5crZ9+/bv9m/fvp2tWLFCgIgKR17Pa3E8ffqUNWzYkAFgysrKbObMX1mbNmJJ0XAbG8YuXRI6StlJTWXst98Yq1z5S6F0S0teQL2wNc3FYjELCAjIVSDd3t6e3b59WzbBkxKVkpLChg4dKvm37d69O/v48aOgMWVlMRYczNimTYxNnsxYly6M1arFP4aqqfEFANTVGatQgTFbW75u15w5jO3cydjNm3yxgsL0azJLTC1cuFCSmGratKm0mxGUIl44SCkRHc1Y+/Zfrm729ow9eCC1wycmJrKxY8dKOkUzMzP2zz//SO34RL7Jsm9bv349mzt3rtSPWxjfJqYyMjKYsrIyO3r0aK7Xubu7s9atWzPG+Eo5NWvWZC9fvmSJiYmsZs2a7N27d/m2kZ6ezhISEiTbixcv5P964e7O+5NmzQp/p0wIKXPk9T7YzMyMXb169bv9169fZ9WrVxcgosKR1/NaVKdPn2Z6enoMADMwMGBjxlxgGhr8cqOlxdjatVL9flWuJScztnw5/wD/+Ra+bVvGQkIKf6y0tDS2fPlyVq5cOQaAiUQiNmrUKBYfHy/9wEmJiIyMZPXr12cAmJKSEluxYgXLyckp8TgyMhgLDOTJJQcH/nv69eqYRdmMjARMTE2fPj3Xt8n79u1jHh4e0m5GUIp24SCljFjM2NatjGlrfxk9tWaNVK/uZ8+eZdWrV5ckqIYMGcLevn0rteMT+aTofdu3ialXr14xAN99kFm2bBmrXbu25PHx48dZrVq12E8//cS2bNlSqDbl/py+eMG/9gIYO3NG6GgIIaWAvPZr6urq7OnTp9/tf/LkCVNXVxcgosKR1/NaFD4+PkxZWZkBYNbWjVnTpjGSD6odOjCWxz9TmfDxI2OennyUCcCYSMTYsGGMvXpV+GO9evWKDR48WHKvXqFCBebt7c2ysrKkHTaRob1790qSjIaGhuzChQsl2v7r14xt3MhY165fJuZ8venqMubkxJiHB2ObNzN27hxjT54wFhPD2PPnfAsPZ+z4cf5xdMwYxhwdGatS5fMxBExMpaamshEjRrDevXuz3r17s+HDh7O0tDRpNyMoRbpwkFLs+XPeU3zuOZo3Z+zhQ6kdPjk5mU2ePJmJRCIGgOnr6zN/f38mphEVCkvR+7b8ElNBQUG5Xrd06VJWp04dqbQp9+d07Fjef7RuTaOlCCEFIq/9Ws2aNdnu3bu/279r1y5mbm4uQESFI6/ntTCysrLYxIkTJcmSVq2GMh2dNAYwVr48Y9u306WGMT4BYsCAL7fw5csztnIlY+nphT/WlStXJNMlATArKyt2SZHnRyqItLQ0Nnr0aMm/m4ODA4uNjS2RtuPiGPP25rd+IlHuRFTlyowNGsTHQNy/z6vJFNWHD4ydOSNgYuqztLQ0lpqaKqvDC0oRLhxEQYjFjG3blnv01Lp1xetFvnH9+nVmZWUl6TjbtGnDwsPDpXZ8Ij8UvW8rylS+otq4cSOrV68eq127tvye0+hoxlRVed9x8aLQ0RBCSgl5vVasWLGC6evrMz8/P/bs2TP27Nkztn37dqavr8+WL18udHj/SV7Pa0ElJCSwTp06Se4XbWy8GCCWzBR//FjoCOXPjRv83HxOCtSqVbRSj9nZ2WzTpk2sYsWKkvM/cOBA9qooQ7GIzEVGRjJra2vJVMx58+axbBnPaxWL+Winvn0ZU1HJnYyyt+eJ0bAwqX6EZIwJXGOqLCjtFw6igJ4/52OjP/cwrVoxFhUltcNnZmayVatWMU1NTQaAqaiosOnTp1NxdAWj6H1bfsXPx40bl2tfvXr1JMXPi0uuz+nIkby/aN9e6EgIIaWIvPZrYrGYzZgxg2loaDAlJSWmpKTEtLS02KJFi4QOrUDk9bwWxMuXLyVfYmpoaDITkyMMYExJibEFC3gRZZK3nBxeLPrL1CfGevRg7Nmzwh/r3bt3bMyYMZLZDuXLl2crV65kGRkZUo+bFJ5YLGY7duxgWlpaktprp06dkmmbaWl8dFTt2rmTUY0b83W0nj+XafOUmJK10nzhIApMLGZsyxY+Hhjgq2xt2CDV1Hd0dDTr0aOH5NsYY2NjtmfPHprepyBKum8riXaSkpJYaGgoCw0NZQDY2rVrWWhoKHv+75X4wIEDTFVVlW3fvp1FRESwyZMns3LlyrFnRbkjzIPcXi+iohhTVuZ9xbVrQkdDCClF5LZf+1dSUhK7efMmu3fvHksvytwogcj7ec1PeHg4MzEx+bfOURWmqXmLATzRQoNxCy4hgbFp076MZtHU5AXTi5JTun37NrO3t5fcr9euXZudPHlS+kGTAktMTMxVE6xt27bs9evXMmsvJYVPojEy+pKM0tZmbNw4PjKqpMg8MbVgwQJ2/Phx9vLly6K8vdQrrRcOUkZER/OlPj73Qm3a8Cp1UvTXX3/lWq62efPmtFytApBm37Z69er/bKskVmy9cOGC5P/p19uwYcMkr/Hx8WFmZmZMTU2N2drasotSvJOW2+vFkCG8f+jcWehICCGljNz2a18Ri8Ul/qXZy5cv2aBBg5ienh7T1NRk1tbWhbo3Kg3n9VuXL19mFSpUYABYxYq1GfCUAXxFrxIql6Nw7t/ntX8+38bXrctXSiusnJwc5u/vzwwNDSX3Pt26dWNRUpxRQQomKChI8rlJWVmZLV26VGZT99LTGVu/nteK+vx/yMSEFzhPTJRJkz8k88SUSCSSDJGtXLkyc3Z2ZnPmzGFHjhyR2rfM8qw0XjhIGZOTw5iPz5d1PsuVY8zXV6qjp9LS0tiyZcskw1FFIhEbNmxYmU1YKwJp9m0aGhrMz88vz+eSkpJYs2bNWL169YrdjryS6xpTkZF8fgXAGCWUCSGFJM/3wb///juztLRkampqTE1NjVlaWrJt27bJvN0PHz4wMzMz5uLiwm7cuMGio6PZ2bNn2eNCFFaS5/Oal4CAAKaurs4AMF1dewa8ZQBjM2fS1L3iEosZ27Urd3Jh+HDGirJAdkJCAps6dSpTUVFhAJiqqiqbMWNGqfl/VpplZWWxBQsWSFaoNDU1ZZcvX5ZJWzk5jO3dy1j16l/+z5ib81LEQs7klHliqkmTJszExITNnz+f+fr6MldXV2ZjY8PU1NSYkpIS09fXZx06dCjKoUuF0nbhIGXYkyd8xNTnHqptWz6iSopevHjBBg0aJPk2RktLi82fP5/qT5VC0uzbDh06xDQ0NL6r6ZSUlMTs7e1Z7dq1WVxcXLHbkXdyeb3o35/3Bz17Ch0JIaQUkst+jTE2d+5cVq5cOTZz5kx2/Phxdvz4cTZz5kxWvnx5NmfOHJm27enpyVq2bFmsY8jrec2Ln58fU1JSYgBYuXLdGZDCtLQYO3RI6MgUy4cPjI0Z8+U2Xl+fMX//oq1sGBERwTp27Ci5Xzc0NGTbt2+XedHtsioqKirXdMqBAweyjx8/yqStS5cYs7X98v/EyIhXd8nMlElzhVIiNaZ27NjBjIyMWOfOndnDf5eoz8zMZLdv32Zbt279rpisIilNFw5CWE4OrzWlqSmz0VOMMXbjxg3WokWLXBe8TZs2sUx56BVJgUi7b9u2bRvT1NRk58+fZ4zxpFSLFi1YrVq1ZDqvXp7I3fXi/v0vawPfuSN0NISQUkju+rV/6evrs3379n23f9++fUxfX1+mbderV49NnjyZ9enThxkYGLCGDRuyrVu3/vA96enpLCEhQbK9ePFCLs/rt3799VfJvZ6q6ggGZDFTU8ZCQ4WOTHFdvcpY/fq5v2d+9KjwxxGLxeyvv/5itWrVkvwbNmzYkJ07d076QZdROTk57LfffpMsGKWjo8P27t0rk7ZiY79UZvhcQ2rpUsaSk2XSXJGUWPHzpKQkNn36dKatrc2mTZtWZkZIyOsFmZAfioriq/V9fVV7+lSqTYjFYnbo0CFWs2bNXAUX//jjD5Yj7fVHidTJom9buXIl09HRYRcuXGAtW7ZkP/30U5ma7il314tx4/jvf+/eQkdCCCml5K5f+1eFChXYozw+rT98+JDp6urKtG11dXWmrq7OZs2axUJCQtjmzZuZhoYG27lzZ77vWbBgQZ41EOXtvH4mFovZnDlzJHGKRNMYIGatWjH25o3Q0Sm+zEzGVqz48j2zujpPQhRlmlZGRgZbs2YN09XVlfx7du3alUVGRko/8DLkyZMnrHXr1rkKnMuizFFWFmO//caYjg7/vyASMTZ6NGPx8VJvqthKfFW+R48esa5du0qGBCo6eb0gE/KfcnJ4T/b16Clvb6mPnsrMzGQbN25kBgYGks7Zzs6O/fPPP7SCnxyTVd82c+ZMpqSkxGrUqMFevHgh1WPLO7m6XiQmflm1k74dJYQUkVz1a1+ZMGECmzJlynf7p06dytzc3GTatqqqKrO3t8+1b+LEiaxZs2b5vqc0jZjKyclhEyZM+CqB5sUAMRs2TNj6NWXRkyeMOTl9+Z7Z0pKxoKCiHevt27dswoQJkhpIysrKbPTo0ezVq1fSDVrBZWVlsTVr1kjq7pYrV475+vrK5Ev5GzcYa9jwy79/o0Z8n7wqzPVCCcWUlZWFtLQ0/PLLLzA1NcWoUaPw4cOH4h6WECILSkqAuztw7x7Qpg2QkgJMnAg4OACPHkmtGVVVVYwfPx5PnjzBwoULUb58eQQHB6NTp05wcHDAlStXpNYWkU+9e/eWbA8fPoSqqioqVaoEd3f3XM8pKh8fH1hYWKBx48ZCh/LF/v1AcjJQuzbg6Ch0NIQQInXbt29H/fr14erqCldXV9SvXx/btm2DkpISPDw8JJu0GRkZwcLCIte+evXqISYmJt/3qKurQ0dHJ9cmj3JycuDq6oqNGzcCEAHYBGAmFi8WYccOQE1N4ADLmBo1gH/+AfbsASpVAsLDgRYtgPHjgYSEwh2rUqVK8Pb2xv3799G9e3fk5ORg69atqFmzJmbPno1Pnz7J5GdQJDdv3kTjxo0xbdo0pKamwsHBAffu3cO4ceOgpFTsVIvEp0+AmxvQrBkQFgZUqABs2gRcvw40aSK1ZgQlYoyxwr5p2bJluHfvHu7du4dHjx6hXLlysLKygrW1NRo2bAgXFxcoKyvLIl65kJiYCF1dXSQkJMjtRYSQ/yQW8x7N05MnqDQ0gMWLgSlTABUVqTb19u1beHl5wdfXFxkZGQCAjh07YvHixWiiKL2pApBm3zZ8+PACvW7Hjh3Fakfeyc31gjHAzg4IDQXWruW/54QQUgRy0699w7GACXeRSITz589Lte2BAwfixYsXuHz5smTflClTcOPGDQQFBRXoGPJ4XrOysjB06FAcOHAAgBKAnVBTG4zt24HBg4WOjrx/D0yfDny+lTI2Bry9gV69AJGo8Me7cuUKPD09Jf9nK1SogGnTpsHd3R3a2tpSjLz0+/DhA+bPnw9fX18wxlCxYkWsWbMGLi4uUk1IMQYcPAhMngy8ecP3DRkCrFkDVK4stWZkplD9WlGGZIlEImZubs4WLVqU51xuRSevQ5gJKZLoaMY6dMg9JjQsTCZNvXjxgo0ZM0ayZC0A1q1bNxYcHCyT9kjhUN8mfXJzTm/c+FKU4v17YWMhhJRqctOvyZGbN28yFRUVtmzZMhYVFcX27t3LtLS02J49ewp8DHk7r+np6axHjx7/3q+pMOAQ09VlLDBQ6MjIt86dY6xmzS+38t27MxYTU7RjicVidvz4cWZhYSG5V9fX12crVqwoM/WkfyQzM5N5e3szPT09yfkZMmQIeyODQmtPnzLWqdOXf9c6dRj7dz2hUkPmNaZatWrFdHR0mEgkYlpaWqxZs2Zs/PjxbPv27SwsLEzhl52UtwsHIcUmFjPm58dYhQq851NRYWz2bMbS0mTS3JMnT5iLi4tkqWEArEePHiyUlnQRlCz6tpii3hkpCLm5Xgwfzn+3hw4VNg5CSKknN/2anPnzzz9Z/fr1mbq6Oqtbt+5/rsr3LXk6r2lpaaxz587/3qOpM+AvVrUqY3fvCh0ZyU9qKmNz5jCmqsov9+XLM7Z+PWNF/VienZ3N9u3bx2rXri25V9fT02MLFy5k78vgF1xisZidPHmS1atXT3I+6tevL5MVDTMzGVu58ktJYDU1xhYvZiw9XepNyVxh+rUiTeX7LCoqCsHBwQgJCUFwcDBCQ0Px6dMnqKuro0GDBrh582ZRDy3X5HGoLSFSERvLa04dOcIf16oFbNkis3o0jx49wpIlS7Bv3z6IxWIAQM+ePbFgwQI0bNhQJm2S/Mmib1NSUkLFihVhbW0tme5tbW2NjIwM+Pj4YNeuXVJpR17JxfXi40egalUgLQ0ICgLs7YWJgxCiEOSiX/vKiBEjCvQ6Pz8/GUdSPPJyXtPS0tCzZ0+cPn0agCaAE6hXrz3++QcwNRUsLFJA4eHAmDHA1av8sa0tv5Vv1Khox8vOzsa+ffuwZMkSPH78GABQrlw5jB07FpMmTYKJiYmUIpdfFy9exNy5cyU1cvX19bFkyRKMGjUKKlIuf3LzJjBqFHD3Ln/s6Ahs3szLg5ZGhenXipWYykt0dDRu376N0NBQLF++XJqHlhvycuEgRGYCAngVxdhY/tjFhU9m1teXSXMPHjzA4sWLceDAAXzuknr06IEFCxbAxsZGJm2S78mib3v27BnCwsIQFhaG0NBQhISE4PXr1wAAHR0dfPz4USrtyCu5uF789hsvTmBlxStmFqXwBCGE/Esu+rWvKCkpwczMDDY2NvjRx5qAgIASjKrw5OG8pqSkoHv37v/W4CoH4G/Y27fBX38BenqChESKQCwGfv+dl5H99Ilf9sePB5YuBXR1i3bMnJwcHD58GF5eXrhz5w4AQFlZGT///DMmT56MZs2aQaRg9xdXr17FwoULcfbsWQB8wYIJEyZg7ty5qFChglTbSkoC5s7lNcIY479vv/4KDBtWum/bZF5jqqyTp6G2hMjMp0+MubkxJhLxcaSVKjHm78+n/clIREQEGzBgABOJRJJhst27d6caVCWkpPq2oKAgVqtWLXb06FGZtiOkjRs3snr16kmGwAt2vRCLGatbl/8Ob9okTAyEEIUib/fB48aNYxUrVmTW1tbst99+K7XTjIQ+r0lJSaxNmzb/3n+VZ8Bl5uzMWEqKIOEQKYiLY2zw4C81iqpUYWzv3uLdyovFYvb3338zBwcHyb06ANa4cWP2+++/l/o6VNnZ2ezw4cOsWbNmkp9NVVWVubm5sZcvX0q9PbGYsSNHGKta9cu/05AhjMXHS70pQci0xtSdO3dYTk5OgV9///59lpWVVdhm5JrQFw5CStTVq4zVr/+lt3RwYCw8XKZNRkREsIEDB+ZKUHXr1o3dvn1bpu2WdSXZt508eZLZ2trKvB2hCX69CAz8UmwiMVGYGAghCkXwfi0P6enpbN++fax9+/ZMS0uL9e3bl/3zzz9MLMMv06RNyPOalJTEWrVq9e89lw4DrrFffmEsI6PEQyEycO4cY7Vrf7mVb9OGsfv3i3/csLAwNnz4cKauri65Xy9fvjwbOXIku3r1aqn6/Xvz5g1bs2YN++mnnyQ/i5qaGnN1dWXR0dEyaTM6mrEuXb78u9Sowdjp0zJpSjAyTUwpKSmx+EKk8LS1tdmTJ08K24xck8cLMiEy9W0VPhUVxmbMYEzG34pERkayQYMG5SqSTgkq2ZFF35aZmZnn/qioKKauri61duSV4NeLX37hv7NjxgjTPiFE4Qjer/2HZ8+esYULF7IaNWowExOTUjOCQ6jzmpiYyFq2bPnvfZYuA26yceOKXjSbyKf0dMaWLct9K+/hwSdIFNebN2/YqlWrchVKB8DMzMzY9OnT2e3bt+UySZWZmcn+/vtv1rt371wrhuvp6bG5c+ey2NhYmbSbns7Y8uVf/i1UVRmbN48XsFc0Mi1+rqSkhNGjR0NLS6tAr/f19UVERARq1KhRmGbkmjzMASdEEM+fA5MmAceP88fVqgFr1wJ9+sh0AvTDhw+xbNky7N27V1IkvVu3bliwYAHs7Oxk1m5ZI4u+TU1NDRYWFrCxsUHDhg1hY2MDY2NjeHt7Izo6GidOnJBKO/JK0OtFfDz/Hc3KAkJCAKrXRgiRAnm/D46JiYG/vz/8/f2RmZmJBw8eoHz58kKH9Z+EOK9JSUlwdnbG1atXAVQAcAazZzfC0qWlu64Nyd/z58CUKbycLABUrgwsX87LySorF+/YjDFcuXIFv//+O44ePYrk5GTJc+bm5nB2doazszMcHR1Rrly54jVWRKmpqTh16hQCAgLw119/5ap12qRJE4wYMQKDBw+WWXynT/N1ph494o8dHIBNm4C6dWXSnOBkWvzcwcGh0IXN9u3bByMjo0K9R57J+wWZEJn780/A3R149ow/btuWV+uzsJBps48ePcLSpUtzJai6d++OBQsWwNbWVqZtlwWy6NuuXLmCO3fu4M6dOwgLC0N4eDjS0tIAAE5OTrCzs4OVlRWsrKxQr149qbQpTwS9XqxcCcycCTRtCly/XrJtE0IUljzeB2dkZODo0aPw8/PDlStX0LVrVwwfPhydOnWCkpKS0OEVSEmf18TERDg7OyMoKAg8KXUWK1bYwdNT5k0TOfDPP3xdlIcP+WNbW75WSsuW0jl+WloaTp48iYMHD+Kvv/6S3PsBvIi4vb09mjdvjhYtWsDe3h4VK1aUTsPfSE9Px40bN3Dx4kVcvHgRQUFBSE9PlzxvYGCAQYMGYcSIEWjQoIFMYgCAmBhg6lTg8GH+2NCQrys1aJBiJ4EFXZWvLJDHCzIhJS4tjX/wXbkSSE/nX7NMmAAsWADI6OLyWV4Jql69emHRokUyvagoupLo28RiMR4+fChZqe9z0io+Ph45OTkyaVNIgl0vxGKgZk0gOhrYsYN/FUoIIVIgb/fBbm5uOHDgAExNTTF8+HAMHjwY+jJaRViWSvK8JiYmolMnZ1y7FgSgIoAz8PGxg5ubTJslciYzE/DxARYuBBIT+b7evYEVK4BataTXTkpKCs6fP4+TJ0/if//7H54/f/7da8zNzWFhYSHZTE1NYWxsDGNjY2hra/9wYEx2djbevHmDly9f4tWrV3j48CHu3buHu3fv4uHDh8jOzs71ejMzM/Tq1Qu9e/dG8+bNoVzcoWI/kJICrFrFt88flyZO5Oe8qCskliaUmJIxebsgEyKo6GjAwwM4dow/rlQJWLIEcHUFVFRk2vTDhw+xePFi7N+/X7JEdL9+/bBo0SLUVdQxsTIkZN/25s0bGBoalmibJUGwc3rpEtCmDb/riY0FNDVLrm1CiEKTt/tgJSUlmJqawsbG5ocfXo8ePVqCURVeSZ3Xb5NSItFZ+PvbYuhQmTVJ5Fx8PDBvHvD77/x7LRUVwM2N76tUSbptMcbw8OFDXL58GUFBQQgKCsKjz/Pa8qGhoYFy5cpBU1MTGhoaUFVVRXp6OtLS0pCeno7ExETJF9V5MTQ0hIODA9q0aYM2bdqgXr16hZ4BVliMAfv3AzNmAK9e8X1t2gAbNgBWVjJtWq5QYuorvXr1QmBgINq1a4fDn8fOAfjrr78wdepUiMVieHp6wtXVtcDHlLcLMiFy4exZXn8qIoI/rl+f15/q0EHmTUdERGDhwoU4dOgQAH6TOmTIECxYsADm5uYyb19RSKtvu3v3LurXr1/g6RPh4eGoU6cOVGScyCxJPj4+8PHxQU5ODh49elTy14uJE4GNG4HhwwE/v5JrlxCi8OTtPtjFxaVAHzJ37NhRAtEUXUmc18TERHTs6Izr13lSSkXlLA4csMXPP8ukOVLKhIfzRMrJk/yxtjb/7tnDA5Dlr/r79+8RHh6O8PBwRERE4MGDB3j16hVev36NhISEAh1DRUUFRkZGqFq1KmrUqIEGDRpINhMTE5knor524QI/j7dv88fVq/Npe717K/a0vbxQYuorFy5cQHJyMnbu3ClJTGVnZ8PCwgIXLlyAjo4ObG1tcePGDejp6RXomPJ2QSZEbmRlAZs38/GpHz7wfZ07A6tXy7z+FADcuXMH8+fPlxTUVlVVxejRozF37lxUqVJF5u2XdtLq25SVlREXFwcDA4MCvV5HRwdhYWEKtUjGZ4JcL8RiXvQ8Nhb4+2/+O0gIIVJC98GyIevz+m1SSlX1LI4ds6VLBPnOuXPA9OlAaCh/rKfHS1aOHw8UcP0zqUlJSUF8fDzS0tIkW1ZWFjQ1NSUjqHR1dWFgYCDTKXkFcf8+4On5JbFXvjw/b1OnAhoagoYmGEpMfSMwMBAbN26UJKaCgoKwevVqBPy7HMGkSZPQrFkzDBgwoEDHowsyIf/hwwdg8WI+cT07G1BSAkaNAhYt4tX+ZOzmzZuYM2cOzp49CwDQ0tLC5MmTMWPGDOiWhQndRSStvo1Wb/1CkOvF1au8eqmuLvDmDaCuXjLtEkLKBHm6D1akEbqyPK+JiYlwcnLGjRs8KaWhcRZ//WWLdu2k2gxRIGIxcPQon8734AHfZ2DAkyxubnw0FeEePOAfew4c4FP4VFSAMWP4uVPAKhWFUph+Ta6Xqbh06RK6desGY2NjiEQiHPtcw+Yrvr6+MDc3h4aGBuzs7HD58uX/PO7r169RtWpVyeNq1arh1efJn4SQ4tPTA9av52OCe/bkV7ctW4CffuLJqa+Wj5WFJk2a4MyZMzh37hyaNGmC1NRULF++HDVq1MDatWtzrcZBpK9169Z4+PAhQkNDC7TZ29tDk2ogSc+/U1rRvTslpQghCs3Gxgbv378v8Ovt7e0RExMjw4jkT2JiIjp0+JKUKlfuLM6epaQU+TElJaBPH+DePcDfHzA3B96+5SOAzMx4IqYQv3oK6dEjYPBgwNKS15NiDPj5Z/7xZ+NGSkoVVrG/LsjKykJcXBxSU1NhYGBQ4OlwBZGSkgJra2sMHz4cP+cx+fngwYOYPHkyfH190aJFC2zZsgXOzs6IiIiAqalpvsfNa5DYj+adZmRkICMjQ/I48fOyBYSQH6tdGwgI4IWYp00Dbt3i0/w2beKr97m6AqqqMmu+bdu2uH79Oo4fP45Zs2bhwYMHmDp1Kry9vbFs2TL88ssvpWYZ6dIkMDBQ6BDKLrEYOHKE/71vX2FjIYQQGWOMYd68eQUeoZuZmSnjiORLYmIi2rd3xq1bPCmlrX0W58/bolEjoSMjpYWKCjBsGDBwIE++LFvGEzILFvCFuV1ceInZ2rWFjrTk3LzJV9k7epQnowD+XeDChYCNjaChlWpF+kSWnJyMLVu2wMHBAbq6uqhevTosLCxgYGAAMzMzjBo1Crdu3Sp2cM7Ozli6dCl69+6d5/Nr167FyJEj4erqinr16mH9+vUwMTHBpk2bfnjcqlWr5hoh9fLlSxgZGeX7ei8vL+jq6ko2ExOTov1AhJRVrVsDN24Af/zBR029ecPHAVtY8H0/WEmjuEQiEXr27Il79+7h999/h7GxMZ49e4ZBgwahadOmuHTpkszaJqTE3bwJvHzJx9iXwMIDhBAiJBqhm7+EhAS0bdtRkpSqUOEsrlyhpBQpGlVVYOhQvsbRgQNAw4ZAairg6wvUrcsTM//8I9NbekHl5AB//gk4OABNm/LvABkDunblRc6PH6ekVHEVusbUunXrsGzZMpibm6Nbt25o0qQJqlatCk1NTXz48AH379/H5cuXERAQgGbNmsHb2xu1atUqfqAiEQICAtCzZ08A/BsPLS0tHDp0CL169ZK8btKkSQgLC8PFixcl+76tMZWdnY169eohMDBQUvz8+vXr0NfXz7PtvEZMmZiYyMXcekJKncxMYNs2PgY4Pp7vs7MDli/nH6RlvFxFamoq1q1bh5UrVyIpKQkAX71z1apVqFmzpkzblnfyVDdEUZT4OZ06la+GOXAgsHev7NsjhJQ5dK2QDWme14SEBDg6dkRo6A0AFVGp0llcvmyLunWlEyshjAGBgfyW46+/vuyvXp2XlR0+HPjBuI9S480bYPt2YOtW4Plzvk9Vld9mTZvGFyEn+ZNpjamgoCBcuHAB169fh6WlJVq0aIEGDRqgZs2aaNKkCUaMGIEdO3bgzZs36N69e64EkTS9e/cOOTk5MPxm8qahoSHi4uIkjzt27Ii+ffvi5MmTqFatGm7dugUVFRX8+uuvcHR0hI2NDaZPn55vUgoA1NXVoaOjk2sjhBSRmhpf1uPxYz7mtXx5IDgY6NgRaNsWuHZNps1raWlhzpw5ePz4McaNGwclJSUEBATAwsIC06dPL/CytITIHcaAf7+AoWl8hBBSNn369Alt2jj9m5TSQ5Uq53D9OiWliHSJRICjIx9FFBkJuLsDFSoAz54Bc+YAJiZAp07Arl1AaauCk5HBK5H06cN/jjlzeFJKT4+vVvj0Ka+7RUkp6SrWqnyampoIDw8vkZWUvh0x9bmAeVBQEOzt7SWvW7ZsGXbv3o0Hn5cPkAH6pogQKXr7FvDy4iv4fa790K0bsHQpYGUl8+bDw8MxdepUnDp1CgBQuXJlLF++HMOHDy9z9aeob5MeHx8f+Pj4ICcnB48ePSqZc3rrFtCkCVCuHP+9KiPTVQghJYuuFbIhjfP68eNHtGnTEffu3QKgDxOTc7hyxRo/KL1LiNSkpvLvx7ZsAYKCvuzX0AC6dOHT/Tp3BipVEi7G/KSnAxcv8vgPHwY+ffryXLNmwLhx/Ds/urUqnBJbla9JkyaIjo4uziGKrFKlSlBWVs41OgoA4uPjvxtFRQiRYwYGfBxwVBQwciRfBuTPP/nk9YED+X4ZsrS0xD///IO///4btWvXRnx8PFxdXdGsWTOp1MojX6Smpub7nFDXElkZP348IiIiSvb/0OfV+Lp2pTsnQggpY96/f48WLdr/m5SqhBo1zuPGDUpKkZKjpcXrUF29ym/fFy3iRdHT03lNpmHD+Ep1rVrxCh5Xr375TrqkMQZERwN+fkCvXjxZ1qkT8PvvPClVtSofHXXnDp/MMXQo3VrJWrESU+7u7pg9ezZevHghrXgKTE1NDXZ2djhz5kyu/WfOnEHz5s1LPB5CSDGZmvKrQUQE0K8fv2Ls3w/Uq8cnq8t4eefOnTvj3r17WLNmDbS1tXHr1i00bdoUo0aNKtRS1CR/urq6mDt3LrKzs797Lq+VV0kh0DQ+Qggps96+fQt7+3aIjAwBYIA6dc7j+nUrhajxQ0qnmjWB+fOBBw94cfC5cwFra14c/coVPj2uZUs+/a99e/748GHgyZMvK91JU0ICT4T5+AADBvApejVq8O/Ejx0DUlIAY2Ng9Gjg/Hk+dW/VqhKZvEH+VaypfJ+nuZQvXx7du3eHg4MDbGxs0KBBA6ipqRU7uOTkZDx+/BgAYGNjg7Vr18LR0RF6enowNTXFwYMHMWTIEGzevBn29vbYunUrtm3bhvDwcJiZmRW7/fzQEGZCSkBYGL+K/f03f6ymBowdC8yezb9ukaHY2Fh4enpi9+7dAAB9fX2sWrUKLi4uCj29T9Z9W926ddG8eXPcuXMHe/fuRd2vCl7Y2NggNDRU6m0KrcSuFyEhfBEBLS0+ja+AS6cTQkhh0X2wbBT1vL558wbNmrXDs2fhAAxhZXUegYEWqFhRdrESUlQxMfzW/vx5Xjz93bvvX6OtDdSqBZib8+RR9ep8RJO+Pt90dQFl5S9bdjavY5WUxP9884a38+IFTzBFRHwpXP41FRWgUSPA2ZkPNrexkfkaTGVOofo1VgzPnj1jx44dY4sWLWK9evViNWrUYEpKSkxVVZU1aNCgOIdmjDF24cIFBuC7bdiwYZLX+Pj4MDMzM6ampsZsbW3ZxYsXi93uf0lISGAAWEJCgszbIqTMu3qVMQcHxvgXKIxpaTE2cyZj79/LvOnLly+z+vXrS/qeFi1asPv378u8XaHIum+zsbFhjDEWEBDATExM2Pr16797TtGU2PVi5kz++9Gnj2zbIYSUefJ2H+zo6Mg+fvyY7/Nv375l5ubmJRcQY2z58uUMAJs0aVKB31OU8/rq1StmYlL33/sUY9a06QOWlFSEgAkRgFjMWHg4Y76+jLm6MmZnx5ia2pdbfmlv1aox5uzM2KJFjJ0/z1hKitBnQPEVpl8r1oipvCQlJSEsLAx3797F+PHjpXlouUHfFBFSwhjjX63Mng3cvMn36erydVonT+Yr+8lIVlYWvL29MX/+fKSkpEBVVRWenp6YM2cONDQ0ZNauEGTdt9na2iIkJAQA/4Z31KhRSE1Nxc6dO9GtWzfJc4qkRK4XjPEiDo8fAwcOAP37y6YdQgiB/N0HKykpIS4uDpUrVwYAZGRkQF1dXfL8mzdvYGxsjJycnBKJ59atW+jXrx90dHTg6OiI9evXF+h9hT2vz58/R9Om7fDmzRMAJnBwOI+TJ2tSHRxSqmVlAY8e8Sl90dF8BbyYGOD9e759+MBHReXk8E0s5qOmtLUBHR3+Z6VKvEKIqSmfslenDtCgAWgUoQAK069JPTFVFsjbBZmQMoMxXhh97lzg3j2+z8CAJ6zGjuXLfsjIixcvMGHCBJw4cQIAULt2bfz+++9o1aqVzNosabLu2/Karrd161YsWbIEaWlpeJfXeO5SrkSuF3fv8sINGhp8Gp8ME7WEECJv98FfJ6bu378PZ2dnDBs2DEuWLIFIJCrRxFRycjJsbW3h6+uLpUuXomHDhvkmpjIyMpCRkSF5nJiYCBMTkwKd16ioKDRr1g4fPrwAUAOdO59DQEB1SKGSCiGESI1MV+WLKWQB4levXhW2CUIIyZtIxNeaDQvjhdFr1uQfxKdM4SNGtm/nE81lwMTEBMeOHcORI0dgZGSER48eoXXr1pgwYQKSkpJk0qaiGTVq1Hf7Ro8ejYsXL6JPnz4CRKQg/vmH/9muHSWlCCFl1uXLl9G6dWs4Ojpi27ZtaNeuHeLj40s0hvHjx6NLly5o3779f77Wy8sLurq6ks3ExKRAbYSHh8POrvW/Sak6+OWXSzhxgpJShJDSrdCJqcaNG2PUqFG4+Xk6TR4SEhKwbds21K9fH0ePHi1WgIQQ8h0lJeCXX3g1w61bgWrVeIVDV1egfn2+rIcMBoOKRCL07t0bERERGDlyJADAx8cHVlZWCAwMlHp7iqZLly5ISEiQPD579izc3d0REBCA3377TcDISrnPq9M6OQkbByGECCQgIACdOnXC3LlzsWvXLoSEhCAzMxMNGzYssevzgQMHEBISAi8vrwK9ftasWUhISJBsBVnl/PbtYDRp4oCkpDgADTBu3EXs21cVysrFDJ4QQgRW6MRUZGQkdHV10alTJxgaGqJLly4YNWoUJk6ciMGDB8PW1haVK1eGv78/Vq9ejYkTJ8oibkIIAVRVgVGjgKgo4Ndf+VIdDx8CffsCTZvyulQyUKFCBfz+++84e/YszMzM8OzZMzg6OsLd3R2pqakyaVMR9O3bFykpKQCA27dvo1+/fjA1NcX9+/cxduxYgaMrpdLSgMuX+d87dBA2FkIIEcikSZOwdetWeHh4AACqVq2KixcvYsCAARg4cKDM23/x4gUmTZqEPXv2FLj+pLq6OnR0dHJtP3LhwiU0b+6I1NR3ABphzpxA+PgY0ipihBCFUOQaU+np6Th58iQuX76MZ8+eIS0tDZUqVYKNjQ06duyI+vXrSztWuSFvc+sJIf9KTOQJql9/Bf5NgKBjR2DlSl6DRwaSkpIwbdo0bN26FQBQt25d7N27F7a2tjJpT5Zk3bdZWVnh7t27AIDJkydDQ0MDK1asAGMMVlZWuPe5bpgCkfn14vRp/n+8alU+apA+oRBCZEze7oOHDx+OgQMHokM+yfljx47h+PHj2LFjh8xiOHbsGHr16gXlr4Yu5eTkQCQSQUlJCRkZGbmey8uPzuuxY3/j55/7QCxOB+CAX389Dg8P4c89IYT8SIkWP4+OjsbFixehrq4OW1tb1KlTpziHKxXk7YJMCPnGmzfA0qXA5s285pRIBAwezPeZmsqkyVOnTmHEiBF4/fo1VFRUsGTJEkyfPv0/b0Tliaz7tgYNGuDu3bsQiUSwtLTEhg0b0K5dO8lzipSY8vHxgY+PD3JycvDo0SPZXS+mTwfWrAGGDwf8/KR/fEII+QbdB38vKSkJz58/z7Vv+PDhqFu3Ljw9PQv0hX1+59XPbz9cXYeCsWyIRN2we/dBDBpES+8RQuSfTIuff239+vWoWbMmJk2ahAkTJsDCwgKNGzfGnTt3inNYQggpHkNDwNsbePAA6N+f15vavZsXSJ85E/iqzpG0dOzYEXfv3sXPP/+M7OxszJo1C05OTnj9+rXU2yqt+vXrhw4dOqBv375QVlaGo6MjAODp06fQ1tYWODrpGj9+PCIiInDr1i3ZNvS5vhRN4yOEkHyFhYXJ9Pja2tqoX79+rq1cuXLQ19cv1iySZcu8MXLkIDCWDRWVQfjf/45QUooQopCKlZhatmwZZs2ahY8fP+L9+/d4/PgxWrVqhebNm+PKlSvSipEQQormp5+AAweAW7cABwcgI4NP66tZE/DxAbKypNqcvr4+Dh06BD8/P5QrVw7nz5+HtbU1Tp48KdV2Sqt58+ZhyZIlGDhwIIKCgqCkxC9B2dnZ2Lhxo8DRlUJv3gCfvwj6d+QZIYQQLiEhAb6+vrCzs0OjRo2EDqdQGGMYN24u5s51B8CgoTEBQUG70LGjqtChEUKITBQrMZWcnAwXFxfJhwtzc3OsXbsWs2bNwtSpU6USICGEFFujRrwQ+p9/AnXrAu/eARMmAFZWgJSTRiKRCMOHD0dwcDAaNmyId+/eoUuXLpg5cyays7Ol2lZpZG9vDxsbG5QvX16yr3bt2nRuiuLsWf6njQ1QubKwsRBCiJw4f/48Bg8eDCMjIyxatAjVq1dHMSuXFElgYCDWr19f6PdlZ2ejR49R2Lx5GQBAT28p7t7dgMaNi/WxjRBC5FqxejgrKytcu3btu/39+/eXFLglhBC5IBIBXbsCd+/y0VKVKvGpfl26AM7OQGSkVJurU6cOrl27JlmZdOXKlWjbti1N7QNfne/rWhynT5/G6NGjBYyolKJpfIQQAgB4+fIlli5dip9++gndu3cHYwyHDx/G69evsWjRIqHDK5QOHQbjzz+3A1BC9erbEBk5B7Vq0cIWhBDFVqzE1K+//oqpU6fi4MGDub6JuHbtGmrVqlXs4AghROpUVQE3N+DxY2DaNP74n3+ABg2AKVOAT5+k1pSGhgY2bNiAQ4cOQVtbG5cvX4aNjQ0CAwOl1kZp5Ofnhz59+uDZs2c4ePAg5s6di1OnTgkdVunCGCWmCCEEQOfOnVGrVi1cu3YNixcvxps3b7B371507twZysrKEJWy1UpDQv4HQB12dkdx/74rDYglhJQJxV6V7+TJkxg7dizS09NhbW2NzMxMhIeHY8+ePejUqZO04pQrtBoJIQokKoonqE6c4I8NDAAvL77KmZL0hs1HRUWhb9++uHPnDpSVlbFixQpMnTpVrm6YS7Jvu3fvHn755RdUrlwZx48fV9i+VGbnNDwcqF8f0NAAPn7kfxJCSAmQt/tgJSUlDBw4EJMnT86zllR4eDisrKyQk5MjQHQF9/m8AhXRvftfOHKkOVRUhI6KEEKKrsRW5QP4txRRUVHw9/eHjY0N1NTUIBKJ0KVLFxgYGKBt27aYPHlycZshhBDZqFULOH4cOHWK1596+xZwdQWaNweCg6XYTC0EBQVh6NChyMnJwfTp09G3b18kJSVJrQ15curUKYjF4lz7GjdujCZNmqBJkyYYOXIkPn36hLdv36J9+/Zo0qSJQJGWUp9HS7VuTUkpQkiZdvXqVWhqaqJt27aoU6cOFi9ejMePHwsdVpFNmnQOx45RUooQUrYUe8RUfl68eIGwsDCEhoYiNDQUAQEBsmhGEPL2TREhREqysgBvb2DhQiApidelGjsWWL4cqFBBKk0wxrBlyxa4u7sjKysLFhYWOHbsmFxMf5Zm36asrIzY2FhU/moOwtd1pfJiZmZWrDblkcyuF1268ML9q1fzEX+EEFJC5PU+ODU1FQcOHICfnx+uXbuGxo0bY9CgQbC0tESHDh1KzYgpeTuvhBBSVIXp12SWmFJkdOEgRMG9fg1Mnw7s28cfGxoCa9cCAwbwZJUUXLt2DT///DNiY2Oho6ODffv2oUuXLlI5dlFJs29TUlJCXFxcrsRUWSST60VGBqCnB6SmAmFhgLW1dI5LCCEFUBrugx8+fIjt27dj9+7dePPmDUQiESWmCCGkhJXoVD5CCFE4xsbA3r3A+fNAnTrAmzfAoEGAkxPw5IlUmrC3t0dwcDCaN2+OxMREdOvWDatXrxZkSWtSyly7xpNShoa8aD8hhJBc6tSpg1WrVuHly5c4evSo4F/8EEII+TFKTBFCSH4cHYE7d4ClS3kdn7NneSJg1So+7a+YjIyMcOHCBYwePRqMMcyYMQPDhg1Denq6FIIX3saNG3Hq1Cm8e/dO6FAUy+f6Uu3bS7VAPyGElFbPnz/H6dOnERsbm2u/srIyevbsic2bNwsUGSGEkIKgO1pCCPkRdXVgzhzg3j2gXTsgLQ3w9ASaNOHTqIpJTU0Nmzdvhre3N5SVlbF79244OTnh48ePxY9dYJs2bYKzszMMDQ1hamqKnj17YsmSJfj777+/+/BACiEwkP/Zvr2gYRBCiDzYv38/atasiU6dOuGnn37C7t27AfBk1YoVK9C0aVOYmpoKHCUhhJAfocQUIYQURM2afKSKvz+v7xMWBjRuDMybx2v+FINIJMKECRPwzz//QEdHB5cvX0arVq3w4sULqYQulPDwcLx8+RInTpzAqFGjIBKJ8Pvvv6Nbt26oVq0ajI2NhQ6x9MnIAG7f5n9v2VLYWAghRA4sWbIEEydOxL1799ChQweMGzcOc+bMwU8//QR/f380adIER48eFTpMQgghP0DFz4uAihMSUsa9eQNMmAAcPswf168P7NoF2NgU+9B3796Fs7MzXr9+jWrVquH06dOoV69esY9bELJele+zDx8+4Pbt2wgLC8OMGTOK1Y68k/r14vp1wN4eqFQJiI+XWjF+QggpKHm7D1ZXV8ejR49gZmaGly9fwtTUFA4ODvDx8Smx66c0yNt5JYSQ4qLi54QQIkuGhsChQ3wzMADu3+dT+5YsKXbtKSsrK1y7dg316tXDy5cv4eDggPDwcCkFXnJ+9J2Hnp4enJycFDop5ePjAwsLCzRu3Fi6Bw4K4n82b05JKUIIAZCVlQVNTU0AQLVq1aCpqYk1a9aUqqQUIYSUdZSYIoSQourTBwgPB3r3BrKzgfnzgVatgMePi3VYU1NTXL58GQ0bNkR8fDwcHBxw9+5dKQVdMv73v/9BV1dX6DAEM378eET8v717j4uqzvsA/hnuiNyRmzBGwpJCeAEvKBpqoGQorE+37UGt7FmeSNfopllbuhlpuV4y2HBbtS3Dto1yXZXQRFt9vA2OsWLmNVBxEVRuwnA7zx+nGRi5zcDMnAE+79fr9zpzzpzzmy+/Ocyc853z+53CQhw/ftywFbdOTBEREQBg27Zt+PHHHwEAFhYWcHV1lTgiIiLSBxNTREQ9MWiQ2KXvs88AZ2fg6FFg5Ejg44+BHvSUdnd3x759+xAeHo6ysjJMnToVp06dMlzcRjZ9+nTY2tpKHUbfIgjA//2f+DgyUtpYiIjMRFRUFN58802EhITAw8MDdXV1WL9+Pb744gsUFhaisbFR6hCJiKgLTEwREfWUTAb85jfADz8ADzwA1NQACxYAjz8OVFR0u1o3Nzfs3bsXY8eORXl5OaZNm9arklNkYEVFwLVrgJUVEBEhdTRERGbh4MGDqKiowI8//oiNGzfihRdeQEFBAZKTkxEaGgoHBweEhYVJHSYREXWCiSkiIkORy4F9+4B33xWTB198IQ6IfuxYt6t0cXFBTk6OVnKqt3XrIwNRd+MbNQoYMEDaWIiIzExQUBAef/xxrF69Gvv27cPNmzdx4cIFfPrpp3j44YelDo+IiDrBxBQRkSFZWgKvvgp8/z1wzz3ApUvAxInA+vXd7tqnTk6NGTMG5eXlmDp1KhQKhWHjJvPHbnxERHoJCAjAI488gnfeeUfqUIiIqBNMTBERGcP48cDJk+IA6Y2NwOLFwKOPApWV3arOxcUF3377rSY5NWXKFBw8eNCwMZN548DnRERERNQHMTFFRGQsLi5id74NGwBra3GQ9IgI4N//7mZ1Lti7dy+io6NRVVWF6dOn45///KdhYybzVFMDKJXiYyamiIiIiKgPYWKKiMiYZDJg4UKxa59cDpw7B4wbB3z+ebeqc3Jywq5duxAfH4+6ujokJiZi586dBg6azM6JE0BTE+DnB/j7Sx0NERG1kpaWhjFjxsDR0RGenp5ISEjA2bNnpQ6LiKjXYGKKiMgUxo0DFAogJga4c0e8i9/vfgc0NOhdlb29Pf7+97/j0UcfRUNDA+bMmYNvv/3WCEGT2VB34+P4UkREZufAgQNISUnBkSNHkJubi8bGRsTGxqKmpkbq0IiIegUmpoiITMXDA9i9G3jtNXF+wwYxUVVaqndV1tbW+PTTT5GYmIj6+nrMnj0b+/fvN3DAZDY4vhQRkdnas2cP5s+fj5CQEIwYMQKbN29GUVERb1RCRKSjfpuYWrt2LUJCQjB8+HAsWrQIQjfvlkVEpBdLS2DlSiA7G3B0BA4cAMLDxa5aerK2tkZWVhZmzpyJuro6/PrXv8alS5eMEDRJShBa7sjHxBQRkdmrqKgAALi5uXW4jkqlQmVlpVYhIuqv+mVi6saNG9i4cSMUCgUKCgqgUChw5MgRqcMiov4kIQE4ehQIDgauXAEmTQKysvSuxsbGBl9++SXGjx+P27dv45FHHoFKpTJ8vCSdc+eA8nLAzg4YOVLqaIiIqBOCICA1NRVRUVEIDQ3tcL20tDQ4Oztrij/HDySifqxfJqYAoLGxEXV1dWhoaEBDQwM8PT2lDomI+pthw8Tk1MyZQF0d8MQTwBtvAM3NelVjZ2eH7du3w83NDQqFAi+++KKRAiZJqLvxRUQANjbSxkJERJ16/vnn8cMPP+DzLm5ysnTpUlRUVGhKcXGxiSIkIjI/vTIxdfDgQcTHx8PX1xcymQxff/11m3XS09MREBAAOzs7hIeH4/vvv9c8N2jQILz00kuQy+Xw9fXFgw8+iKFDh5rwLyAi+oWzM/DNN8DLL4vzb78tJqj0vOpJLpfj008/BQB8+OGH2L59u6EjJakcPSpOOfA5EZFZW7hwIXbs2IH9+/fDz8+v03VtbW3h5OSkVYiI+qtemZiqqanBiBEjsHHjxnaf3759OxYvXoxly5bh5MmTmDRpEuLi4lBUVAQAuHXrFnbu3InLly/j6tWrOHz4MA4ePGjKP4GIqIWlJbB6NbB1K2BtDXzxBTBjBvDLGBW6iouLw2u/DKy+ZMkSNHTjjn9khvLzxWlEhLRxEBFRuwRBwPPPP4+vvvoK3333HQICAqQOiYioV7GSOoDuiIuLQ1xcXIfP//GPf8QzzzyDBQsWAADWrVuHnJwcZGRkIC0tDXv37kVgYKBmQMKZM2fiyJEjmDx5crv1qVQqrTFbODghERnF3LnA4MFAYiKQlyeOO5WTA/j46FzF8uXLcefOHaSmpsLa2tp4sZJpNDYCP/wgPub4UkREZiklJQXbtm3DN998A0dHR1y/fh0A4OzsDHt7e4mjIyIyf73yiqnO1NfXQ6FQIDY2Vmt5bGwsDv8yToe/vz8OHz6Muro6NDU1IS8vD8HBwR3WycEJichkpk0DDh4EvL2BggLggQeAq1d13tzKygpr167l51Rf8dNP4vhjAwcCgYFSR0NERO3IyMhARUUFoqOj4ePjoynsVk9EpJs+l5gqKytDU1MTvLy8tJZ7eXlpfr0YP348HnroIYwaNQphYWEYOnQoZs2a1WGdHJyQiExq5EhxwOshQ8Q7sj3wAMDPnf7p5ElxOmIEYNHnvrKJiPoEQRDaLfPnz5c6NCKiXqFXduXThUwm05oXBEFr2cqVK7Fy5Uqd6rK1tYWtra1B4yMi6lRAAHDgADBlCnDhgpicyssD5HKpIyNTUiem2I2PiIiIiHqqpgYoKwNu3gRu3RKnt28DlZViqagAqqvFUlMjTuvqxBszqVRAfb1Yj0wmFktLwM4OsLdvmQ4YIE6tdE839bnElIeHBywtLTVXR6mVlpa2uYqKiMisDRminZyKjmZyqr9RKsXpqFGShkFEREREZqy2VuxhoS5Xr4rl2jWgpAQoLRXLnTtSR9quPpeYsrGxQXh4OHJzc5GYmKhZnpubi9mzZ0sYGRFRN/j7i8mo6GgmpySSmJiIvLw8TJs2DV9++aXpXlgQWq6YYmKKiIiIqP9qaACKioCLF1vKpUtiuXxZvApKVzY2gJubWFxdxeLk1FIcHcXxTQcOBBwcxCuh7OwAW1txW0A8ThUE8UY9KpWYGKurExNftbViuXULePddnULqlYmp6upqnD9/XjN/6dIlKJVKuLm5QS6XIzU1FUlJSYiIiEBkZCQyMzNRVFSE5ORkCaMmIuomP7+2yal9+8TufmR0ixYtwtNPP42tW7ea9oWLi8XLq62sgJAQ0742EREREZlWVZV4rH93uXhRTEo1NXW+/cCB4o/a/v7inb4HDwZ8fcU7fHt5AZ6eYhk4UOyGZ2yVlX07MXXixAlMmTJFM5+amgoAmDdvHrZs2YLHHnsM5eXlWLFiBUpKShAaGopdu3ZhyJAhUoVMRNQzdyenoqKA3Fxg+HCpI+vzpkyZgry8PNO/sLob3/Dh4i9URERERNR7CYLYnU6dbDp/XjsBVVra+fZ2dsC994o/Tqun6iKXAy4upkk4GUGvTExFR0dDEIRO13nuuefw3HPPmSgiIiIT8PMDvv8eiIkBTp8GJk8Gdu8GxoyROjLJHDx4EO+99x4UCgVKSkqQnZ2NhIQErXXS09Px3nvvoaSkBCEhIVi3bh0mTZokTcD6YDc+IiIiot6lpkbsWnf5ckt3u9bd72pqOt/e3R0YOrSlBAa2PPby6rN3ae6ViSkion7Lx0ccEP2hh4Bjx4Bp04C9e4GxY6WOTBI1NTUYMWIEnnrqKcyZM6fN89u3b8fixYuRnp6OiRMn4qOPPkJcXBwKCwsh/2WcrvDwcKhUqjbbfvvtt/D19dUrHpVKpVVXZWWlnn9RK7wjHxEREZH5EATgxo2WAcaLioCffxZLUZGYjLpxo/M6ZDKxq92992onoNTFxcUUf4nZYWKKiKi3cXcXk1Hx8WKSavp0ccyp0aOljszk4uLiEBcX1+Hzf/zjH/HMM89gwYIFAIB169YhJycHGRkZSEtLAwAoFAqDxZOWlobly5cbpjJeMUVERERkfIIA3L4N/Oc/4h3srl8Xp9eutdzdTl3a+TGzDRcX4J572na7GzpUvOs2h2hog4kpIqLeyNER2LkTmDEDOHRI7N63fz8QFiZ1ZGajvr4eCoUCS5Ys0VoeGxuLw4cPG+U1ly5dqhn3EBCvmPL399e/ops3xV/eAF4xRURERKQrlUpMMlVUiNNbt8Tjqps3gfJysZSVieXGDXFcp9JS8a53upDJAG9vcYgNuVxMNKnLPfeI03561VNPMDFFRNRbDRwI7NoFxMYCR48CcXHigNmDBkkdmVkoKytDU1MTvLy8tJZ7eXnh+vXrOtczffp05Ofno6amBn5+fsjOzsaYDsb1srW1ha0hfgVTD3weEAA4O/e8PiIiIiJDEwTxTnHNzeK0dWlsbDtVl4YGsdTXt0zr68WkkrrU1bWU2lrgzp2WUlMDVFeL06oqsVRWilNdrmjqiJOTOGxG66K+u93gwWIyytcXsLExXBsSACamiIh6NycnYM8eIDIS+PFH4MknxQHRLS2ljsxsyO66O4kgCG2WdSYnJ8fQIXWN3fiIiIhITRC0kzDV1eK0pqalqJM2tbUtRZ3gUSd71AkgdVEniBoatBNHrRNK6sftJZ+6uCGZpJydxeLmJhZXV3Hq4aFdvLwAT0/xh107O6mj7reYmCIi6u1cXIAvvxQHQM/NBd5+G3jzTamjkpyHhwcsLS3bXB1VWlra5ioqQ/vwww/x4YcfoqmpqXsVqK+YYmKKiIio72hsFLuU3bghlrKylu5l5eXic+ruZ7dvt3RJq6wUE0O9jYWF+GOppSVgbS2W1o9tbFqmtrbi1MZGTBCpi60tMGBAS7G3F3sNODi0TJ2cxOLo2PKYP9L2KkxMERH1BSEhwJ/+BMydCyxfDkyYII471Y/Z2NggPDwcubm5SExM1CzPzc3F7NmzjfraKSkpSElJQWVlJZy70xWPd+QjIiLqHQRBTB5duyaW1oNnX78uDqitLuXlPbvKyMJCTMY4OorFwaGlqJM29vbi49aJHTs77eSPrW1LcsjKSjtpZGXVkjxSJ5XUy1snmlqXjpbrcYU69W9MTBER9RVJScD33wObNgELFgDnz4sHFX1YdXU1zp8/r5m/dOkSlEol3NzcIJfLkZqaiqSkJERERCAyMhKZmZkoKipCcnKyhFF3obZW7JYJ8IopIiIiKQmCeFVTcXFLuXJFvDvblSstj2tr9avXzU3sOubhId5t2d29pcuZm5t4NbyrqzhVd0lzchITTkz2UB/ExBQRUV+yYYN42ffrr/f5pBQAnDhxAlOmTNHMq++IN2/ePGzZsgWPPfYYysvLsWLFCpSUlCA0NBS7du3CkCFDpAq5a83NwAcfiMkpX1+poyEiIuq7VCox2VRUBPz8s1iKilpKcbE4RpMuXFzE721fX/GubT4+4vhF3t7iVD2Wkbu7ePUREWnIBMGcRywzT+quGRUVFXBycpI6HCIig+Bnm+GxTYmor+HnmnGwXY2kqqol4dReuX5dt6513t7iHdn8/Vumre/S5usrXs1ERBr6fK4xVUtERGRgPR78nIiIiDrX3CyO29T6Cif1FU/qxNOtW13XY28PDBkiFrm8bRk8WByTiYiMhokpIiIiA+vx4OdERET9WXMzUFqqPZbTlSttx3pqaOi6LlfXlsRTe8XDg+M2EUmMiSkiIiIiIiIyvoYG4MYN8Uon9Z3r1HevU9/V7upVcVljY9f1WViIVzT5+7dc8dR6OmSIePc6IjJrTEwRERERERGRfgQBqK4Gbt4US3m5eAe7sjLx8Y0bYiktbSllZbrXb2Ehju00eHBL8qn1GE9yuTi2EwcSJ+r1+F9MRERERETUQ+np6XjvvfdQUlKCkJAQrFu3DpMmTZI6rLaamoCampZSXd1SqqrEUlkpTisqtMvt2+K4TeqpLlc13c3CQrw7nbd3y93rfHzEJNPgwdp3tmPSiahf4H86ERGRgXHwcyKi/mX79u1YvHgx0tPTMXHiRHz00UeIi4tDYWEh5HK57hXt2AFYW4td3hoaxMRPQwNQX98ybV1UKu1SV6ddamvFcudOy1SlMuwfb2sLuLuLYzkNGiSO2eThIS7z9BTLoEHi1MtLXG5padgYiKhXkwmCLvfHpNYqKirg4uKC4uJi3s6ViPqMyspK+Pv74/bt2xyw20D4fUFEfQ2/K9o3btw4jB49GhkZGZplw4YNQ0JCAtLS0tqsr1KpoGqVIKqoqIBcLkcxAJN+Wzg4iGXgwJbHjo6Ak5O4zNERcHYW552dxeLi0jJ1dRXvasfBw4noLvp8X/CKqW6oqqoCAPj7+0scCRGR4VVVVfFkw0D4fUFEfRW/K1rU19dDoVBgyZIlWstjY2Nx+PDhdrdJS0vD8uXL2yw3+beFujtfaampX5mI+gldvi+YmOoGX19fFBcXw9HREbJffh0YM2YMjh8/rlmn9bz6sTpjaIxfzu9+fUOs39k67T2ny7KO5o3VNvq2iy7b6NsuHS3XZZ8BYDZtw32mZ9vos290tNzY+4wgCKiqqoKvr6/O21Dn2vu+MAfG/D4yNsYuDcYuDXOMnd8VbZWVlaGpqQleXl5ay728vHD9+vV2t1m6dClSU1M1883NzQgPD0d+fn6H3xfGOJaQch/rzjGXoerp6XFtV8/zPTH8eYYu6/E9MVw9hnhPBEFAeHi4Tt8XTEx1g4WFBfz8/LSWWVpaau2krefvfs7JycngO/Tdr2GI9Ttbp73ndFnW1byh20bfdtFlG33bpaPl+uwzgPRtw32mZ9vos290tNwU+wx//Tas9r4vzIkxvo9MhbFLg7FLw9xi53dF++5OKAmC0GGSydbWFra2tm2Wdda2xjiWUJNiH+vOMZeh6unpcW1Xz/M9Mfx5hi7r8T0xXD2Gek9sbGxgYWHRZT1dr0E6SUlJ6XD+7udM8fqGWL+zddp7TpdlXc0bWnfq72obfdulo+XcZ/rPPtPZOr1pnyEiIqK2PDw8YGlp2ebqqNLS0jZXUXWmu8egPTmWkJKhYjH1sZsuz/M9Mfx5hi7r8T0xXD2meE9a4+DnJlRZWQlnZ2dUVFSY1a9e5oBt0zG2TfvYLh1j21BXevM+wtilwdil0Ztj72/GjRuH8PBwpKena5YNHz4cs2fPbnfwc3PBfcz88D0xP3xPjI9d+UzI1tYWb775ZpvLdolt0xm2TfvYLh1j21BXevM+wtilwdil0Ztj729SU1ORlJSEiIgIREZGIjMzE0VFRUhOTpY6tE5xHzM/fE/MD98T4+MVU0RERERERD2Unp6O1atXo6SkBKGhoVi7di0mT54sdVhERGaPiSkiIiIiIiIiIpIEBz8nIiIiIiIiIiJJMDFFRERERERERESSYGKKiIiIiIiIiIgkwcQUERERERERERFJgokpM1VcXIzo6GgMHz4cYWFh+Nvf/iZ1SGYjMTERrq6u+K//+i+pQ5Hczp07ERwcjKCgIPz5z3+WOhyzwv2kLX6u0MGDBxEfHw9fX1/IZDJ8/fXXUoekk7S0NIwZMwaOjo7w9PREQkICzp49K3VYOsnIyEBYWBicnJzg5OSEyMhI7N69W+qwuiUtLQ0ymQyLFy+WOpQuvfXWW5DJZFrF29tb6rB0dvXqVfz3f/833N3dMWDAAIwcORIKhULqsKif47GV9Hjsb374f2EYTEyZKSsrK6xbtw6FhYXYu3cvXnjhBdTU1EgdlllYtGgRPvnkE6nDkFxjYyNSU1Px3XffIT8/H6tWrcLNmzelDstscD9pi58rVFNTgxEjRmDjxo1Sh6KXAwcOICUlBUeOHEFubi4aGxsRGxvbK/ZfPz8/vPvuuzhx4gROnDiBqVOnYvbs2Th9+rTUoenl+PHjyMzMRFhYmNSh6CwkJAQlJSWaUlBQIHVIOrl16xYmTpwIa2tr7N69G4WFhVizZg1cXFykDo36OR5bSYvH/uaJ/xeGYSV1ANQ+Hx8f+Pj4AAA8PT3h5uaGmzdvwsHBQeLIpDdlyhTk5eVJHYbkjh07hpCQEAwePBgA8NBDDyEnJwdPPPGExJGZB+4nbfFzheLi4hAXFyd1GHrbs2eP1vzmzZvh6ekJhUKByZMnSxSVbuLj47XmV65ciYyMDBw5cgQhISESRaWf6upqPPnkk9i0aRPefvttqcPRmZWVVa+6Skpt1apV8Pf3x+bNmzXL7rnnHukCIvoFj62kxWN/88T/C8PgFVPdpEt3iPT0dAQEBMDOzg7h4eH4/vvvu/VaJ06cQHNzM/z9/XsYtfGZsl16u5621bVr1zRfTID4q/zVq1dNEbrRcT9qnyHbpTd9rhDdraKiAgDg5uYmcST6aWpqQlZWFmpqahAZGSl1ODpLSUnBzJkz8eCDD0odil7OnTsHX19fBAQE4PHHH8fFixelDkknO3bsQEREBB555BF4enpi1KhR2LRpk9RhkZnjsZP547G/+eH/jflgYqqbuuoOsX37dixevBjLli3DyZMnMWnSJMTFxaGoqEizTnh4OEJDQ9uUa9euadYpLy/H3LlzkZmZafS/yRBM1S59QU/bShCENtvIZDKjxmwqhtiP+iJDtUtv+1whak0QBKSmpiIqKgqhoaFSh6OTgoICDBw4ELa2tkhOTkZ2djaGDx8udVg6ycrKQn5+PtLS0qQORS/jxo3DJ598gpycHGzatAnXr1/HhAkTUF5eLnVoXbp48SIyMjIQFBSEnJwcJCcns6sIdYnH4OaPx/7mh+ccZkSgHgMgZGdnay0bO3askJycrLXsvvvuE5YsWaJzvXV1dcKkSZOETz75xBBhmpyx2kUQBGH//v3CnDlzehqi2ehOWx06dEhISEjQPLdo0SLhs88+M3qsptaT/aiv7SetdbddevvnChlOe/tQb/Dcc88JQ4YMEYqLi6UORWcqlUo4d+6ccPz4cWHJkiWCh4eHcPr0aanD6lJRUZHg6ekpKJVKzbIHHnhA+N3vfiddUN1UXV0teHl5CWvWrJE6lC5ZW1sLkZGRWssWLlwojB8/XqKIqLfhMbj547G/+eE5h7R4xZQR1NfXQ6FQIDY2Vmt5bGwsDh8+rFMdgiBg/vz5mDp1KpKSkowRpskZol36C13aauzYsfj3v/+Nq1evoqqqCrt27cL06dOlCNekuB+1T5d26YufK9S/LFy4EDt27MD+/fvh5+cndTg6s7GxQWBgICIiIpCWloYRI0Zg/fr1UofVJYVCgdLSUoSHh8PKygpWVlY4cOAANmzYACsrKzQ1NUkdos4cHBxw//3349y5c1KH0iUfH582V9QNGzaMv9BTt/HYyfzx2N/88P/GtDj4uRGUlZWhqakJXl5eWsu9vLxw/fp1neo4dOgQtm/fjrCwME1f17/+9a+4//77DR2uyRiiXQBg+vTpyM/PR01NDfz8/JCdnY0xY8YYOlxJ6dJWVlZWWLNmDaZMmYLm5ma88sorcHd3lyJck9J1P+oP+0lrurRLX/xcof5BEAQsXLgQ2dnZyMvLQ0BAgNQh9YggCFCpVFKH0aVp06a1uZPdU089hfvuuw+vvvoqLC0tJYpMfyqVCmfOnMGkSZOkDqVLEydOxNmzZ7WW/fTTTxgyZIhEEVFvx2Nw88djf/PDcw7TYmLKiO7u8ysIgs79gKOiotDc3GyMsCTXk3YBgJycHEOHZLa6aqtZs2Zh1qxZpg7LLHTVNv1pP2mts3bpy58rpJvq6mqcP39eM3/p0iUolUq4ublBLpdLGFnnUlJSsG3bNnzzzTdwdHTUHBA6OzvD3t5e4ug699prryEuLg7+/v6oqqpCVlYW8vLy2txp0Bw5Ojq2GcfLwcEB7u7uZj++10svvYT4+HjI5XKUlpbi7bffRmVlJebNmyd1aF164YUXMGHCBLzzzjt49NFHcezYMWRmZnJcQOoxHoObPx77mx+ec5gGE1NG4OHhAUtLyza/QJSWlrbJuPYnbBfdsa06xrZpH9uFdHHixAlMmTJFM5+amgoAmDdvHrZs2SJRVF3LyMgAAERHR2st37x5M+bPn2/6gPTwn//8B0lJSSgpKYGzszPCwsKwZ88exMTESB1an3blyhU88cQTKCsrw6BBgzB+/HgcOXKkV1x1NGbMGGRnZ2Pp0qVYsWIFAgICsG7dOjz55JNSh0a9FI8RzB/fI/PD98S0OMaUEdjY2CA8PBy5ublay3NzczFhwgSJopIe20V3bKuOsW3ax3YhXURHR0MQhDbFnJNSANqNWT1mmrn7+OOPcfnyZahUKpSWlmLv3r29OimVl5eHdevWSR1Gl7KysnDt2jXU19fj6tWr+Pvf/95r7oQIAA8//DAKCgpQV1eHM2fO4Nlnn5U6JOrFeIxg/vgemR++J6bFK6a6qavuEKmpqUhKSkJERAQiIyORmZmJoqIiJCcnSxi18bFddMe26hjbpn1sFyIiImoPjxHMH98j88P3xIyY8A6Afcr+/fsFAG3KvHnzNOt8+OGHwpAhQwQbGxth9OjRwoEDB6QL2ETYLrpjW3WMbdM+tgsRERG1h8cI5o/vkfnhe2I+ZIIgCMZLexEREREREREREbWPY0wREREREREREZEkmJgiIiIiIiIiIiJJMDFFRERERERERESSYGKKiIiIiIiIiIgkwcQUERERERERERFJgokpIiIiIiIiIiKSBBNTRERERAawfv16BAQEYMCAAUhISEBFRUW760VHR0Mmk0Emk0GpVJo2SB1FR0dj8eLFem0zf/58zd/19ddfGyUuIiIi6nuYmCIiIiLqoddeew0bN27E1q1b8a9//QsnT57E8uXLO1z/2WefRUlJCUJDQ00Ype6++uor/OEPf9Brm/Xr16OkpMRIEREREVFfxcQUERERUQ8cP34cq1atwvbt2zF58mSMHj0av/3tb7Fz584OtxkwYAC8vb1hZWXV7detr6/v9rZdcXNzg6Ojo17bODs7w9vb20gRERERUV/FxBQRERFRD7z//vuYOnUqRo8erVk2aNAglJWV6VXPnj17EBUVBRcXF7i7u+Phhx/GhQsXNM9HR0fj+eefR2pqKjw8PBATEwMAaG5uxqpVqxAYGAhbW1vI5XKsXLlSa7uFCxdi8eLFcHV1hZeXFzIzM1FTU4OnnnoKjo6OGDp0KHbv3q21TeuufNHR0Vi0aBFeeeUVuLm5wdvbG2+99ZaeLUVERETUFhNTRERERN2kUqnwj3/8A4mJiVrLa2tr4ezsrFddNTU1SE1NxfHjx7Fv3z5YWFggMTERzc3NmnW2bt0KKysrHDp0CB999BEAYOnSpVi1ahXeeOMNFBYWYtu2bfDy8tKqe+vWrfDw8MCxY8ewcOFC/O///i8eeeQRTJgwAfn5+Zg+fTqSkpJw586dDuPbunUrHBwccPToUaxevRorVqxAbm6uXn8jERER0d2YmCIyIXMY8JaD0xIRGU5+fj5qa2vx4osvYuDAgZry8ssvIzg4WK+65syZg1//+tcICgrCyJEj8fHHH6OgoACFhYWadQIDA7F69WoEBwfjvvvuQ1VVFdavX4/Vq1dj3rx5GDp0KKKiorBgwQKtukeMGIHXX38dQUFBWLp0Kezt7eHh4YFnn30WQUFB+P3vf4/y8nL88MMPHcYXFhaGN998E0FBQZg7dy4iIiKwb98+/RqMiIi6zRzOJTry4osvIj4+vsPneQ5CnWFiisjE7h7wVv0hnZyc3Gbd5557DjKZDPPnzzfY63NwWiIiw/npp59gZ2eHgoICKJVKTRk6dCgmTpyoV10XLlzAb37zG9x7771wcnJCQEAAAKCoqEizTkREhNY2Z86cgUqlwrRp0zqtOywsTPPY0tIS7u7uuP/++zXL1FdYlZaW6lQHAPj4+HS6PhERdWzy5MmaRI2NjQ2GDRuGbdu2dbmd1OcSavPnz8eSJUs080qlEiNGjOhwfZ6DUGeYmCIyksbGxnaXtzfgrb+/P7KyslBbW6tZVldXh88//xxyudygcXFwWiIiw6msrISnpycCAwM1xcbGBj/++CPmzJmjV13x8fEoLy/Hpk2bcPToURw9ehSA9iDnDg4OWtvY29vrVLe1tbXWvEwm01omk8kAQKvboC51dLY+ERG1TxAEKJVKvP/++ygpKcHZs2cxY8YMzJ07F5cuXep0W6nPJQDxu+Kf//wnZs+erVl26tSpThNTPAehzjAxRWQAly9fhkwmw5dffonJkyfD1tYW2dnZOm8/evRoyOVyfPXVV5plX331Ffz9/TFq1CitddWD3z7//POaAXJff/11CIKgWaergXCJiMgwPDw8UFlZqfUZvHLlSjz00EMYPny4zvWUl5fjzJkzeP311zFt2jQMGzYMt27d6nK7oKAg2Nvbs0sdEVEvcu7cOVRVVWHGjBnw9vZGQEAAnnnmGTQ1NeHs2bN612fqc4lDhw7BwsIC48aNAwAUFxejvLwcFhYWiImJwYABAxAcHKz5gYWoK0xMERmAuo+3evDZ06dPIzY2Vq86nnrqKWzevFkz/5e//AVPP/10u+uqB789evQoNmzYgLVr1+LPf/6z5nldBsIlIqKemzp1Kurq6vDuu+/i8uXLeOedd7Bjxw5kZGToVY+rqyvc3d2RmZmJ8+fP47vvvkNqamqX29nZ2eHVV1/FK6+8gk8++QQXLlzAkSNH8PHHH3f3TyIiIiNTKBRwdXXV/IBx5coVLFu2DLa2tlrdrPVhynOJHTt2ID4+HhYWYjpBfS70wQcfYOnSpTh16hTkcrlWVz+izlh1vQoRdeXUqVNwcHDA3/72N9xzzz3dqiMpKQlLly7VXH116NAhZGVlIS8vr826/v7+WLt2LWQyGYKDg1FQUIC1a9fi2Wef1QyEu3HjRsybNw8ANIPhEhGRYXl5eWHLli14+eWX8Yc//AFTp07Fv/71L/j7++tVj4WFBbKysrBo0SKEhoYiODgYGzZsQHR0dJfbvvHGG7CyssLvf/97XLt2DT4+Pu2ONUJEROYhPz8fFRUVcHR0RHNzM2pra2Fvb48//elPGDx4cLfqNOW5xI4dO/D+++9r5pVKJVxdXfHFF1/A09MTAJCQkKD3jzTUfzExRWQASqUSs2bN6nZSChC7g8ycORNbt26FIAiYOXMmPDw82l13/PjxmvFAACAyMhJr1qxBU1OTzgPhEhGRYTz22GN47LHHelzPgw8+qHUHPgBaXSvaO7kAxKTWsmXLsGzZsnafb2+7y5cvt1nW2Wu1VwfvqkRE1D0KhQIpKSlYtGgRbt++jZdeegmRkZE9GqTcVOcSZ86cwZUrV/Dggw9qlimVSsyePVuTlAKAixcvIjAwsNt/D/Uv7MpHZACnTp3S6Vftrjz99NPYsmULtm7d2uGlt13RdSBcIiKSTnp6OgYOHIiCggKpQzGY5ORkDBw4UOowiIjM3smTJzFhwgQEBgYiIiIC6enpWL16dZcDn3fFFOcSO3bsQExMjNZ6SqUSkZGRWuudPHkSI0eO7FYM1P8wMUXUQ5WVlbh8+XKbgQW7Y8aMGaivr0d9fT2mT5/e4XpHjhxpMx8UFARLS0sOhEtEZOY+++wzFBYWQqlUIjg4WOpwDGbFihVQKpU4d+4cYmJipA6HiMgsXbx4Ebdv30ZoaKhm2fDhwxEYGIjPP/+8R3Wb4lzim2++waxZszTzVVVVuHTpUptzIaVSycQU6Yxd+Yh66NSpU7CwsOj2QIWtWVpa4syZM5rHHSkuLkZqaip++9vfIj8/Hx988AHWrFkDQHsgXBsbG0ycOBE3btzA6dOn8cwzz/Q4RiIi6pnujh9i7jw9PbW6cRARUVsKhQJWVlb41a9+pbU8JiYG2dnZeO2117pdt7HPJeLj43H8+HGtrtxKpbLNudDPP/+MW7duMTFFOmNiiqiHTp06hfvuuw92dnYGqc/JyanLdebOnYva2lqMHTsWlpaWWLhwIf7nf/5H8zwHwiUiIiIiMj/5+fn41a9+BRsbG63lMTExSE9Px5UrV+Dn59ft+o15LvGPf/wD48aN0/oRor1zoZMnT8LFxaVH4+9S/yITWo90SURGFR0djZEjR2LdunWS1gEAMpkM2dnZSEhI6FE9RERERERkfFKfS8yaNQtRUVF45ZVXuv36PAeh9nCMKSITk3rAWw5OS0RERETUO0l5LhEVFYUnnniiW9vyHIQ6wyumiEzo6tWrqK2tBQDI5fI2l/Dqoqe/lJSWlqKyshIA4OPjAwcHh27VQ0REREREpmMO5xLdxXMQ6gwTU0REREREREREJAl25SMiIiIiIiIiIkkwMUVERERERERERJJgYoqIiIiIiIiIiCTBxBQREREREREREUmCiSkiIiIiIiIiIpIEE1NERERERERERCQJJqaIiIiIiIiIiEgSTEwREREREREREZEkmJgiIiIiIiIiIiJJMDFFRERERERERESS+H9Kd/zsMom+6AAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "f, ax = plt.subplots(1, 3, figsize=(12, 3))\n", "\n", "ax[0].loglog(r_arr, prof_baryons.real(cosmo, r_arr, M_halo, a_sf) * fb, color=\"red\")\n", "ax[0].loglog(r_arr, prof_matter.real(cosmo, r_arr, M_halo, a_sf), color=\"black\")\n", "ax[0].set_xlabel(\"$r$ [Mpc]\")\n", "ax[0].set_ylabel(r\"$\\rho(r)$ [$M_\\odot$ Mpc$^{-3}$]\")\n", "\n", "ax[1].semilogy(theta, T_kSZ_150 * glasz.constants.sr2sqarcmin, color=\"red\")\n", "ax[1].set_xlabel(r\"$\\theta$ [arcmin]\")\n", "ax[1].set_ylabel(r\"$T_{\\rm kSZ}$ [$\\mu$K $\\cdot$ arcmin$^2$]\")\n", "\n", "ax[2].semilogx(R, R * ds_b, label=\"baryons\", color=\"red\")\n", "ax[2].semilogx(R, R * ds_dm, label=\"DM\", color=\"blue\")\n", "ax[2].semilogx(R, R * ds, label=\"total\", color=\"black\")\n", "ax[2].set_xlabel(r\"$R$ [Mpc/$h$]\")\n", "ax[2].set_ylabel(r\"$R \\Delta\\Sigma$ [Mpc $M_\\odot$ pc$^{-2}$]\")\n", "ax[2].legend()\n", "ax[2].set_ylim(0, 14)\n", "ax[2].set_xlim(0.05, 60)\n", "plt.tight_layout()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "glasz", "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.11.11" } }, "nbformat": 4, "nbformat_minor": 2 }