{ "cells": [ { "cell_type": "markdown", "id": "a086c311-a05a-4077-85d2-717bcbbc215a", "metadata": {}, "source": [ "### Small scales on the PySM2 beta synchrotron template" ] }, { "cell_type": "markdown", "id": "56cb133d-1d33-422c-82b2-c2c4ccfab7d0", "metadata": {}, "source": [ "- we subtract the offset term from the $\\beta_s$ map \n", "- we use the fitted power law from the study by Nicoletta Krachmalnicoff to inject small scales on to the spectral index map\n", "- we rescale the large scale $\\beta_s$ map by a factor `~1.572` which comes from the SPASS latest results \n", "- we modulate the small scales on the $\\beta_s$ map by the same modulation map adopted for intensity " ] }, { "cell_type": "code", "execution_count": 1, "id": "989c7685-3708-45eb-914f-719da0319837", "metadata": {}, "outputs": [], "source": [ "import os\n", "\n", "os.environ[\n", " \"OMP_NUM_THREADS\"\n", "] = \"64\" # for jupyter.nersc.gov otherwise the notebook only uses 2 cores" ] }, { "cell_type": "code", "execution_count": 2, "id": "4f48bbb4-a285-4585-813f-93de83f91f7f", "metadata": {}, "outputs": [], "source": [ "from pathlib import Path\n", "import healpy as hp\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "\n", "# import pymaster as nmt\n", "from astropy.io import fits\n", "\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 3, "id": "0d4ad523-661b-4981-937c-2d882020b616", "metadata": {}, "outputs": [], "source": [ "plt.style.use(\"seaborn-talk\")" ] }, { "cell_type": "code", "execution_count": 4, "id": "b0c70a96-13c7-4963-8c71-b5c16f6a4898", "metadata": {}, "outputs": [], "source": [ "import pysm3 as pysm\n", "import pysm3.units as u" ] }, { "cell_type": "code", "execution_count": 5, "id": "e09fa591-c024-4325-8feb-439ebac13f3e", "metadata": {}, "outputs": [], "source": [ "nside = 512\n", "lmax = int(1.5*nside)\n", "ell = np.arange(lmax+1)" ] }, { "cell_type": "code", "execution_count": 6, "id": "306e6dc4-66c6-4639-8c37-1b0af5f4a2ad", "metadata": {}, "outputs": [], "source": [ "datadir = Path(\"data\")" ] }, { "cell_type": "code", "execution_count": 7, "id": "a1f4fb6a-03a5-46dc-81bf-cd45194f736e", "metadata": {}, "outputs": [], "source": [ "output_dir = Path(\"production-data/synch\")" ] }, { "cell_type": "code", "execution_count": 8, "id": "af7d24e1-c47e-401c-8c06-342594effcc1", "metadata": { "tags": [ "parameters" ] }, "outputs": [], "source": [ "output_nside = 2048\n", "# output_nside = 8192" ] }, { "cell_type": "code", "execution_count": 9, "id": "abea969d-d50a-45d3-ba66-3d4fc6fe0d21", "metadata": {}, "outputs": [], "source": [ "output_lmax = int(min( 2.5 * output_nside, 8192*2))\n", "output_ell = np.arange(output_lmax + 1)\n", "output_cl_norm = output_ell * (output_ell + 1) / np.pi / 2\n", "output_cl_norm[0] = 1" ] }, { "cell_type": "markdown", "id": "bdc9fc22-f9fd-4146-8251-20d0774753b5", "metadata": {}, "source": [ "The modulation used here is different from the one used for the Synchrotron templates, these are based on a heavy smoothing of the input temperature maps.\n", "The actual algorithm used is available in the documentation of PySM 3.4.0b5 " ] }, { "cell_type": "code", "execution_count": 10, "id": "8739464f-3c93-48dd-a59e-15f7f85a0c76", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "-2.9999962\n", "7.461009e-08\n" ] } ], "source": [ "template_dir = \"/global/project/projectdirs/cmb/www/pysm-data/pysm_2/\"\n", "beta_s_pysm = hp.read_map(template_dir + \"synch_beta.fits\")\n", "print(beta_s_pysm.mean())\n", "beta_s_pysm_nomono = hp.remove_monopole(beta_s_pysm)\n", "print(beta_s_pysm_nomono.mean())" ] }, { "cell_type": "code", "execution_count": 11, "id": "56bae391-7618-44cc-945e-fbf27023e9fc", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/tmp/ipykernel_64107/2230800712.py:1: RuntimeWarning: divide by zero encountered in power\n", " cl_spass_fit = 10 ** (-0.82655497) * output_ell ** (-2.60860208)\n" ] } ], "source": [ "cl_spass_fit = 10 ** (-0.82655497) * output_ell ** (-2.60860208)\n", "cl_spass_fit[0] = 0" ] }, { "cell_type": "code", "execution_count": 12, "id": "736969cf-d2c3-4361-972f-41b5d87a99e7", "metadata": {}, "outputs": [], "source": [ "from pysm3.utils import sigmoid" ] }, { "cell_type": "code", "execution_count": 13, "id": "4e0b6047-4011-4fe5-b1f3-08f10951266b", "metadata": {}, "outputs": [], "source": [ "ell_fit_low = {\"TT\": 10, \"EE\": 10, \"BB\": 10}\n", "ell_fit_high = {\"TT\": 36, \"EE\": 36, \"BB\": 36}\n", "\n", "sig_func = sigmoid(output_ell, x0=ell_fit_high[\"TT\"] - 15, width=ell_fit_high[\"TT\"] / 10)" ] }, { "cell_type": "markdown", "id": "d7a65d0b-eb10-452e-ac0a-d4a0490f14ad", "metadata": {}, "source": [ "## Large scales" ] }, { "cell_type": "code", "execution_count": 14, "id": "f1245e83-2f58-4c04-b9d5-1f40665aa8b9", "metadata": {}, "outputs": [], "source": [ "np.random.seed(444)\n", "\n", "cl_spass_fit *= np.sqrt(sig_func)\n", "\n", "bs_ss_alm = hp.synalm(cl_spass_fit, lmax=output_lmax)" ] }, { "cell_type": "code", "execution_count": 15, "id": "0f2aa5d1-a6b7-42e3-ae47-5d2b9f481cf1", "metadata": {}, "outputs": [], "source": [ "## Includes rescaling from S-pass\n", "\n", "alm_bs = hp.map2alm(beta_s_pysm_nomono * 1.5722501, lmax=lmax, use_pixel_weights=True)\n", "\n", "bs_LS_alm = np.empty_like(alm_bs)\n", "\n", "bs_LS_alm = hp.almxfl(alm_bs, np.sqrt(1.0 - sig_func[:lmax+1]))" ] }, { "cell_type": "code", "execution_count": 16, "id": "059293d5-95bd-43ae-b8df-1226a68d893c", "metadata": {}, "outputs": [], "source": [ "if output_nside > 4096: # only write at 8192\n", "\n", " hp.write_cl(\n", " output_dir / \"raw\" / f\"synch_small_scales_beta_cl_lmax{output_lmax}.fits\",\n", " cl_spass_fit,\n", " dtype=np.complex128,\n", " overwrite=True,\n", " )\n", "\n", " hp.write_alm(\n", " output_dir / \"raw\" / f\"synch_largescale_beta_alm_nside{nside}_lmax{lmax:d}.fits\",\n", " bs_LS_alm,\n", " lmax=lmax,\n", " out_dtype=np.complex64,\n", " overwrite=True,\n", " )" ] }, { "cell_type": "code", "execution_count": 17, "id": "683f6e56-5a65-43dc-9e71-44981bb52046", "metadata": {}, "outputs": [], "source": [ "if output_nside < 4096:\n", " bs_ss = hp.alm2map(bs_ss_alm, nside=output_nside)\n", " modulate_amp = hp.alm2map(hp.read_alm(output_dir / \"raw\" / f\"synch_amplitude_modulation_alms_lmax768.fits.gz\").astype(complex), output_nside)\n", " bs_ss *= modulate_amp\n", " bs_ls = hp.alm2map(bs_LS_alm, nside=output_nside)\n", "\n", " bs_out = bs_ss + bs_ls - 3.1\n", " # output_map_alm = hp.map2alm(bs_out, lmax=output_lmax)\n", " clout = hp.anafast(bs_out, lmax=lmax)\n", " clin = hp.anafast(beta_s_pysm, lmax=lmax)" ] }, { "cell_type": "code", "execution_count": 18, "id": "213f8a3d-0119-447f-ba2c-9f5ff26f1529", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAApkAAAHKCAYAAABIXVCWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAACyRklEQVR4nOzdd1zV9f7A8df3DPbeIoqAe+Pe4sqyzFW5Glaa1q2sn1a32x52u3UbN7NlppVmZmY7TS1HqbjALQ4EZCh7z3PO9/fHF0GUPQT0/Xw8iO/5zs8BoTef8X4rqqoihBBCCCFEfdI1dgOEEEIIIcS1R4JMIYQQQghR7yTIFEIIIYQQ9U6CTCGEEEIIUe8kyBRCCCGEEPVOgkwhhBBCCFHvrqsgU1GU3oqi/K0oynZFUf5QFCWwsdskhBBCCHEtUq6nPJmKorQAslVVzVIUZRwwXVXVuxq7XUIIIYQQ1xpDYzfgalJVNeGSl4WAqbHaIoQQQghxLWtyw+WKokxTFGWHoiiZiqJcEQQqiqJXFOVNRVGSFEXJUhRlnaIoHjV8hj3wOvDf+mq3EEIIIYQo1eSGyxVFGQu4AbbAJ6qqGi47/gxwD3AjkAJ8BtipqnqToigOwNZybvuTqqovFV9vBXxffO/vG+htCCGEEEJc15pckHmRoighwOZygsxo4GVVVZcVvw4CTgMBqqpGVXFPPbAG2KCq6qcN0GwhhBBCCEEzm5OpKIoz0BrYf3GfqqpnFEXJBLoDUVXc4na0HlAPRVHuBA6rqvpIDZ7vDrgXv0xRVTWlBs0XQgghhLhuNKsgE3Aq/pxx2f70S45VSFXVr4Gv6/D8R4AXAKysrGjbtm0dbiWEEEII0bydOXMmWVVVz/KONbcgM6v4s/Nl+12AzKvw/MXAVwAdOnSI2Llz51V4pLje5ebmsmvXLgYOHIidnV1jN0eI6478DApRMUdHx+iKjjWrIFNV1XRFUWKAXkA4QHFCdSfg0FV4fgraYiP69OmDg4NDQz9SiBJ2dnbyb06IRiQ/g0LUTFNMYaRXFMUGsCp+bVP8oRSf8gnwlKIoAYqiOAH/ATZWtehHCCGEEEJcPU0uyATuAvKAjYC+eDsP8C8+/jrwE7AXiCs+586r30whhBBCCFGRJjdcrqrqCmBFJcfNwMLiDyGEEOIKJpMJk6l+iroVFBRgMBhKPgtxvTAYDHX6Ny8/LUIIIa4pWVlZqKqK0WikdKZV7dnY2NC3b19sbGzqoXVCNA+qqpKbm4uiKDg6OtbqHhJkCiGEuGZYLBZMJhOurq71dk+z2UxRURE2Njbo9fp6u68QTZ2trS1paWlYLBZ0uprPsGyKczKFEEKIWikqKsLa2rqxmyHENcPa2pqioqJaXStBphBCiGuGqqq16nERQpRPURRqW4JcfhKFEEIIIUS56jKvWYJMIYQQQghR7yTIFEIIIYQQ9U6CTCGEEEIIUe8kyBRCCCEaUUFBAQsWLMDPzw83NzemTJlCVlZWYzdLiDqTIFMIIYRoRM888wwHDx4kPDyc6OhoIiIi+Prrr6t9fUhICNbW1jg4OODs7ExwcDDr1q2r9vWRkZHcfvvt+Pj44ODgQKtWrZg0aRKFhYUl91cUhW+++abMdaGhoSiKQps2bUr2PfXUU3Tp0gUnJyd8fX2ZM2cOqamp1W6LuLZIkCmEEEI0kqKiIpYtW8bbb7+Nh4cHjo6O+Pv7k5ubW6P7PPfcc2RnZ5OSksL06dOZOnUqJ0+erNa148aNo0WLFkRERJCVlcWuXbsYO3ZsmbQ1nTp1YunSpWWuW7p0KZ06dSqzT6/Xs3LlSlJSUjh48CCxsbHce++9NXov4tohQaYQQohrltmikp5bWOePjLyiap9rtlQ/p2BYWBiKotC9e3cAcnJy+OuvvwgODq7V+zUYDDz00EOYzWYOHz7Mhx9+SI8ePcqcc+bMGQwGA9HR0aSkpBAREcG8efNwdnZGURT8/PyYN29emaT2kydPJiwsjMjISEAr3blu3borAsjXXnuN4OBgjEYjnp6ePPzww2zdurVW70U0f1JWUgghxDUrK7+I5X9H1ekeqmqhoKAQa2srFKXqvpl7B7fBxc6qWvcOCwsrCSijoqJ49NFH6dSpE0OHDq1VWwsLC1myZAlGo5EePXrg5eXFE088wd69e+nbty8Ay5YtY/To0fj7+wPQpUsXZs+ezbx58+jTpw+dOnW6IjeijY0NM2fOZNmyZSxatIjVq1czfPhwWrRoUWl7tmzZUhJAi+uP9GQKIYQQjSQ8PJzg4GAWLVpEQEAAv/zyCwsWLKhxAuxFixbh4uKCn58fP/zwA+vWraNt27Y4OTkxbdo0li1bBmh12D///HPmzJlTcu3WrVsJCQnh3XffpWfPnnh7e/PKK69cUeVlzpw5LF++HJPJxCeffFLmHuVZt24dS5cu5X//+1+N3ou4dkhPphBCiGuWo42Rewe3qdM9zGYz2dnZODg4oNfrq/XM6goLC+ORRx5h5syZ/OMf/+Dbb79l2rRpdO/enQ4dOlT7Ps888wzPPvtsucfmzp3L6NGjefvtt9myZQsmk4lbb7215LiHhwevvfYar732Grm5uXzzzTfMmTOHli1bct9995Wc17VrV/z9/XnllVe4cOECN954I6tXry73mWvXrmXu3Ln8+OOP9OrVq9rvQ1xbJMgUQghxzdLrlGoPXVfEbDajNxtxtLOqVpBZXRaLhcOHD5cMl7u4uDB16lTmzp1bUi/6scce4+DBg+Tn5zNv3jxmzZpV4+f07duXoKAg1q5dy/r165k1axZGY/mBsJ2dHbNmzWLx4sWEh4dfcfyBBx7g/vvv5/nnn6/wa7F8+XIWLFjATz/9xODBg2vcXnHtkCBTCCGEaAQRERHk5uayZs2aktXh//d//8eoUaNo3749x48f5+TJkyULZywWS62f9cADD/DWW29x4sQJ3nzzzZL9aWlpvPHGG8ycOZMOHTqgKAo//PADR44c4Z///OcV95k+fTqtWrWid+/e5T7nvffe46WXXmLjxo0lc0DF9UvmZAohhBCNICwsjKCgINavX4+Liwtdu3bF2tq6JB9l+/btGTJkCFOmTOGVV16p07NmzpzJ2bNnGTx4MO3atSvZb2VlRWJiIpMnT8bNzQ1PT09effVVFi9ezO23337FfWxsbBg9ejSurq7lPmf+/PlkZmYyYsQIHBwcSj7E9Um5fGKvqJ4+ffqo+/bta+xmiOtAdnY2W7ZsYdSoUfLLWogq5OfnA1owVF/MZjNZWVk4OjrW63D5E088QW5uLkuWLCn3eE5ODvb29gD07t2bH3/8kZYtW9bqWaqqEhgYyKJFi5gxY0at2yyuP1X9TCmKsl9V1T7lHZPhciGEEKIRhIeHM3ny5AqPz507l8jISEwmExMnTqx1gAmwatUqCgsLue2222p9DyFqSoJMIYQQohGEh4fz/PPPV3h85cqV9fIcT09PDAYDy5Ytw8qqboughKgJCTKFEEKIRpCUlHRNPUeIy8nCHyGEEEIIUe8kyBRCCCGEEPVOgkwhhBBCCFHvJMgUdZJ8LoLo8K0gqbCEEEIIcQkJMkWtXYg6js2y4fh/P4GwxTMozMtp7CYJIYQQoomQIFPUisVkIv2r+3EgD4Dg1F+JfWsIidEnGrllQgghhGgKJMgUtbL/65fpUHgUgN32o7CoCoGmSKyXj+LYtrWN3DohhBBCNDYJMkWNRR0Npcep9wHY5TyOAU98R9iwT8hQ7XEmm45/zGHP8iewmM2N3FIhhBBCNBYJMkWNFObnoX73AFaKmTi86XyvVnO396g7yLxnC6f0QegUlX7Rn3DkvzeSmZLYyC0WQgghRGOQIFPUSPgXCwkwR2FRFVJv+B/OLm4lx1oFdsJvwQ5CnccB0D1vD9nvD+HMob8bq7lCCNHkFRQUsGDBAvz8/HBzc2PKlClkZWU1drOuSyEhIbz66qu1ujY2NhZFUYiKiqrfRjVjEmSKaosI/Y0+casA2OUzk26DbrriHFs7e/o99hW7u75AoWrAV71Ay3UT2LN+8dVurhBCNAvPPPMMBw8eJDw8nOjoaCIiIvj666+rfX1ISAjW1tY4ODjg7OxMcHAw69atq/b1kZGR3H777fj4+ODg4ECrVq2YNGkShYWFJfdXFIVvvvmmzHWhoaEoikKbNm1K9j311FN06dIFJycnfH19mTNnDqmpqdVuy7WuOl/r6nwvX331VRRF4YsvvqjxM6o6Xp8kyBTVkpOZitOGR9EpKmd0beg9680Kz1UUhQG3/R9RE77jPJ7YKEX0O/gsoe/dTX5e7lVstRBCNG1FRUUsW7aMt99+Gw8PDxwdHfH39yc3t2a/K5977jmys7NJSUlh+vTpTJ06lZMnT1br2nHjxtGiRQsiIiLIyspi165djB07FvWS/MedOnVi6dKlZa5bunQpnTp1KrNPr9ezcuVKUlJSOHjwILGxsdx7773Vfh9FRUXVPrc5qs7XuqrvpcViYdmyZbi5ufHxxx/X+BnVaUN9kSBTVMux5Q/TQk2kUNVjmfgxNrZ2VV7TvtdwbP6xg8M2vQDon/oDMf8dRnx09X7xCSHEtS4sLAxFUejevTsAOTk5/PXXXwQHB9fqfgaDgYceegiz2czhw4f58MMP6dGjR5lzzpw5g8FgIDo6mpSUFCIiIpg3bx7Ozs4oioKfnx/z5s3D2tq65JrJkycTFhZGZGQkAFlZWaxbt+6KAPK1114jODgYo9GIp6cnDz/8MFu3bq2wvW3atOHll19mxIgR2Nvbl/Ta5ebmsnDhQgICAnBzc+PGG2/k9OnTJde99957BAQE4OjoSMuWLfnXv/5Vciw7O5uFCxcSGBiIo6MjXbp04a+//gLgf//7Hx07dsTR0ZHWrVvz9NNPY65gkWpVbTh//jy33norzs7OtG/fng0bNlT2ran21/qiy7+XF23cuJHY2Fi++OILdu7cyZEjR6r9jJq2oa4kyBRVOrhpFX3TfgFgX+BDtOs+oNrXuni2oPPCTYT63QdAe/MpbJaPJHzr+gZpqxBCXKrIbCE6JaeOH7mcS8sjOiW3WucXmS3Vbl9YWFhJQBkVFcX06dPp1KkTQ4cOrdX7LSwsZMmSJRiNRnr06MHMmTM5c+YMe/fuLTln2bJljB49Gn9/f9zd3enSpQuzZ8/miy++4NixY+X2aNnY2DBz5kyWLVsGwOrVqxk+fDgtWrSotD1btmwpCaArsnTpUt5++22ys7OZMGECALNnz+bEiRPs3r2b8+fP079/f2655RaKioo4efIk//znP/n555/Jysri6NGj3HrrrSX3u//++wkNDWXLli1kZmby/fff4+PjA4Cfnx+//fYbmZmZ/PDDD3z22Wd8+umn5barsjYAzJw5E71eT0xMDNu3b2fFihWVvs/qfq0vuvx7edHHH3/MTTfdxM0330yPHj345JNPqv2MmrahrpSGvPm1LLCVt7rnWBQejraN3ZQGlXohFj4ciBuZHDN2of2T2zAYjbW616E/VtNm+//hRC5mVWGn/0MMuudV9Hr5W6cy2dnZbNmyhVGjRuHg4NDYzRGiScvPzwe0oAggOiWH4W9uvapt2PZECP7u9tU698EHH8Te3h5XV1eeffZZdDodX3/9Nbfffnu1nxcSEkJoaCjW1tZYWVnRtm1bnn76acaPHw9owZLBYOCjjz7CbDbTunVr3nvvPaZMmQJAcnIyb7/9Nhs2bODIkSO4uLjwyCOP8Oyzz6IoCiEhIYwePZqJEydyww03EBMTw4ABA3jppZdIS0vj2WefLXexy7p165g1axbbtm2jV69e5ba9TZs23HfffTz//PMl+5KTk/H09CQ6OprWrVsD2hCxq6srv/zyC76+vnTp0oXPP/+ccePGlfm9mJiYiLe3N0eOHKFLly5Vfu0WLlxITExMyXzTi+913rx5lbYhICAAPz8/Tp8+TVBQEACbNm3ihhtu4OzZs2XmqV6qOl/ryr6X8fHx+Pv7s3btWiZOnMh7773HCy+8QHx8PLa2ttV6RlXHL3f5z9TlFEXZr6pqn/KOyf/da8nNnMzxt28mLCKysZvSYFSLhZjP5+BGJjmqDc7TP611gAnQfeR0cu7ezFl9G/SKytCYJYT992bSUpPrsdVCCNF8XOzJfOaZZ0hLS+Pjjz9m2rRpRERE1Og+zzzzDOnp6SQmJrJz586SoARg7ty5rF69mtzcXH799VdMJlOZnj8PDw9ee+01Dhw4QHp6Om+88QYvv/wyy5cvL/OMrl274u/vzyuvvMKFCxe48cYbK2zP2rVrmTNnDj/++GOFAeZFlwdkZ8+eBaB79+64uLjg4uKCm5sbRUVFnDt3jsDAQFatWsXSpUvx9fVlyJAh/P777wAlwW779u3Lfdbq1avp27cv7u7uODs7s2TJEpKSkq44r6o2xMbGAuDv719yTUBAQKXvE6r3ta7se3lxLuYtt9wCwJ133kleXh5r1qyp9jOq+/2uD4Z6v+N1ZKi6n9ivxvB9v3eZMO6Wcv8CaM72/rCEfrk7ATjS/Z/0D+xc53u2COxC/oK/CPvkfoLTN9Inbycxi4dwYeLndOzRv873F0KIS/m62LLtiZA63cNstpCTk429vUO1Rl58Xao3wmWxWDh8+HDJcLmLiwtTp05l7ty5KIqCqqo89thjHDx4kPz8fObNm8esWbNq3P6+ffsSFBTE2rVrWb9+PbNmzcJYQYeBnZ0ds2bNYvHixYSHh19x/IEHHuD+++/n+eefR6/Xl3uP5cuXs2DBAn766ScGDx5cZft0urJf04uB26lTp/D09Cz3msmTJzN58mQKCwv56KOPmDBhAikpKSUB66lTp+jcuez/s86dO8edd97Jd999x0033YSVlRULFy5k3759V9y/qjZcDDKjo6NLejIvBqbVVdXX+nIWi4VPP/2U9PR0/Pz8SvabzWY++eSTcv9tVPWMmrahpqQns5aKHP0oxIifksy4Pffw9ZLnyMyr/+X/jSU+KoLO4YsACLcdSL9J8+vt3jZ2jgTPX8OBrs9QqOpprSbQ+rvx7PjuwwadGyKEuP4Y9Tr83e3r+GFHK1db/N3tqnW+sZpTgCIiIsjNzWXNmjWYTCbS09OZP38+o0aNon379pw4cYKTJ0+ydetWdu/ezd13313rr8MDDzzAW2+9xa+//srs2bNL9qelpfH0009z5MgRioqKMJlMrFu3jiNHjpQ7L3T69On8/vvvzJ9f/v8T3nvvPRYuXMjGjRurFWCWx8vLixkzZvDQQw8RFxcHQHp6OuvXryc7O5uIiAg2bNhAbm4uRqOxZAGLTqfDy8uL2267jYceeoioqChUVeX06dOcPn2a7OxsLBYLnp6eGI1Gdu/ezZdfflmrNvj5+RESEsKTTz5JZmYmFy5c4JVXXqn0fdX0a325DRs2EBsby86dOwkPDy/5+OWXX9i1axeHDx+u8hl1bUNNSZBZS0ZHT0z3biTJ6IuVYmZ68mL2/XcSJ6LiG7tpdWY2m0lbdT8OSh5pOOF3z1IUXT3/U1EUet32JLETvyVJccNOKWDooX+y4737yZU0R0KI60BYWBhBQUGsX78eFxcXunbtirW1dcn8wPbt2zNkyBCmTJlSZQBTlZkzZ3L27FkGDx5Mu3btSvZbWVmRmJjI5MmTcXNzw9PTk1dffZXFixeXOy/UxsaG0aNH4+rqWu5z5s+fT2ZmJiNGjMDBwaHko6aWLl1Khw4dCAkJwdHRkW7durF27VoURaGwsJCXXnqJFi1a4OLiwnvvvce6detK5gx+9tln9OzZk+HDh+Po6MiECRM4f/48nTp14qWXXmLChAm4uLjw+uuvM3369Fq1AeCrr76ioKCAVq1aMXTo0Cr/CKjp1/pyH3/8MRMnTqR37974+PiUfNxwww0MHDiQjz/+uMpn1LUNNSULf2qpT58+6r59+1Dz0ji3/F5aJ/4JwFm1BSeHLWHsqFGN3MLa2/Xliww88w4ABwa+T6+xdzXo8zKT44j/dDod8w8CcETfCce7VuLfpm2DPre5kIU/QlRfVYsUasNsNpOVlYWjo2OFQ8S18cQTT5Cbm8uSJUvKPZ6Tk4O9vbaAqHfv3vz444+0bNmyVs9SVZXAwEAWLVrEjBkzat1mcf2RhT+NSLF1pfWD64nr9ywm9AQoCQzbPo2vl75OXmH5ubeasjNH9tD7tFadZ5/z2AYPMAGcPFrSfuEWwlppfwV2NR/HfvlI9m77ucGfLYQQjSU8PJyuXbtWeHzu3LkMGjSIfv36MXHixFoHmACrVq2isLCQ2267rdb3EKKmZOFPfVAUWo57gsygAZi+mYWbOZlpcf9mw3/30fG+j2jj49HYLayWgoI8+O4BrBQT5xVPOtz74VV7ts5gJPj+xZz4oz+tti/AQ8nA7o/7iGq5lTZt677gSAghmprw8PAyqXsut3Llynp5jqenJwaDgWXLlmFlZVUv9xSiOqQnsx45dRiKy2O7iXbRkpXfWLiJ/A9Hsn3XrkZuWfXsX/EUQRZtdVzamHdxdHGv8poCk5mcAlO9taHjyDvJvWczaThhpxSQsvYxVEv1ExsLIURzkZSU1CCLLcp7TkJCAuPGjWvwZwlxKQky65nO0RP/R38luvtjWFDoqEQTvGES3365pEZVIK62I7t/p3/8FwDs8ZlGp0G3VHlNXqGZlbtjWLojktDIlHpbGe4Z0I3Yfs8A0LsglJ2/fF4v9xVCCCHE1SNBZkPQ6fGf/BLpk9eQoTjjqORx25l/seGtezmfmtnYrbtCVmYaLhsfQa+oROta0WPW21Veo6oqm49fIDOvCFWFnWdS+PFgPPlF9TMPtdtNczlu0xOAoP2vkJqWWi/3FUIIIcTVIUFmA3LrPha7R3dyzkGr2zo+93sS3xvFnvBDjdyyso4sfxQ/9TxFqh7zpE+wtqm6HNqxhExOJ2YD4OOsrTiLTMrhq9AYEjPz694oRcHtjvcpVA34kMKRlU/V/Z5CCCGEuGokyGxgRlc/Wj3+B2fb3wdAd07Sdv1N/LD2cyyWxk8ftW/T1wxM+xGA8KC5BHYbVOU1GXlFbI3QynAFeNgzrW8rRnXyQq9TyMgrYs3ecxyJy6hz27wDu3E4QPu6DU5ey5F9O+p8TyGEEEJcHRJkXg16IwEz3uHCuM/IVuxxU7IZf2Q+P737D1Kz8hqtWUkX4vD/W+shPGnsSK8ZL1d5jcWisvHoeQpNFmyt9Izp7I2iKHT3c+GOPq1wtDFgsqhsOnaBTccuYKrjPNQe018iTtcCvaKi//X/KCwsqtP9hBBCCHF1SJB5FXn3m4LxwR3E2rRHp6hMyFxF5NtjOBxx8qq3RbVYiPp8Hp6kk6ta4zh9GXpD+bVsLxV2Lo24NC0wHt3JC3vr0ixYPs42zOzvj7+7HQBH4jJYs+8cGbm1DwwN1nbkjXkDgE6Wk+z65s1a30sIIYQQV48EmVeZtVcQfgt2cNr/DgD6qIfx+moMv/387VWt273r+w/pm7sdgBPdn6RFYMUJgS9Kyirg79MpAHT2daKtl+MV59ha6ZnYsyUDAt1RFEjMLGDVnmgik7Jr3da2A2/loMtoAHqdeo/YmKha30sIIYQQV4cEmY3BaEPbe5cSO/I98rDBW0lnzN45/LDkSbLzCxv88efOnqTrwVcBOGLbl+BJ/1flNSazhQ1Hz2O2qDjZGhne3rPCc3U6hYFB7kzo2RIbo56CIgs/hMez80xyreehBsz8H1nY4ajkEfv1/KsakAshhBCi5iTIbER+w+7BPPsP4o3+GBQLE5M/4fCb4zgdHdNgzzSZTKR9NRsnJZd0HPC9ZxmKrup/BrsiU0jOKkBR4IbO3tgYq67fG+Bhz4z+rfF20lafh0am8n14XK3KbTp5+hHZYyEAA3K3Err52xrfQwghhBBXjwSZjczBrwstFu7klM/NAAw078X2s5H8ueW3BnneztX/pnvRQQBiB72Km49/ldfEpuWyPzoNgF6tXWnlZlft5znbGrmjjx/dWjoDEJ2Sy6rQaBIyar7gqfuExzht1RGAln8/S2ZW08s5KoQQQgiNBJm1lJFXRKGpfir4KNYOtJu7iqiBiyjASEsliUHb7+THpS+SX1h/JRtPHdlHv9P/AyDcZQxdb7i3ymsKTGY2Hr2AqoKHozWDgqouNXk5g17H6M7e3NDFG4NOISvfxNp9sRw8l16jYW9Fp8du0nuYVB2tOM+BVc/VuC1CCNHUFBQUsGDBAvz8/HBzc2PKlClkZWU1drOEqDMJMmspv8jM13tjSMkuqJ8bKgptxj5M7l2/cUHfAmvFxK1x7xD630nEnk+s8+3z8/OxfDcXG6WIRMWddrM+qtZ12yKSyMwrQq9TGNvFG4O+9v9kuvg6M7VfK1zsjJgtKn+cSCxJh1Rdvp36c8hvOgCDEr7kxOF9tW6PEEI0Bc888wwHDx4kPDyc6OhoIiIi+Prrr6t9fUhICNbW1jg4OODs7ExwcDDr1q2r9vWRkZHcfvvt+Pj44ODgQKtWrZg0aRKFhYUl91cUhW+++abMdaGhoSiKQps2bcq8l4CAAJycnPDy8uK2224jJqbhpoCJpk2CzFpSUEjJLmT1nhhOnK+/YVvXoL54LNjFabcQAIYXbqfoo+Hs2lW3ROShnz9NB8tpANLHvIu9i0eV15xOzOZovPbeBgW54+VoU6c2AHg52jC9X2sCPbWqQscTslizN4a0nOoveOoy898kKu5YKWYKf3gMk6l+SlkKIcTVVlRUxLJly3j77bfx8PDA0dERf39/cnNza3Sf5557juzsbFJSUpg+fTpTp07l5MnqpccbN24cLVq0ICIigqysLHbt2sXYsWPLjDR16tSJpUuXlrlu6dKldOrUqcy+u+66i/DwcDIzM4mKiqJ169ZMmzatRu9FXDskyKwlV3sjDtYGiswqvx0+zx8n6p54/CK9nSttH/me08FPU4SeAOLpuWEyP3/xdq2ecXD3FgbHrwDggM/ttB90a5XX5BSY2Hz8AgAtXW3p1dq1xs+tiI1Rz609fBnSzgNFgeTsQr7aE8PpxOoND1nbOZM6bBEA3U2H+fu7JfXWNiHENcZcBKmRdfw4iy49GlLPVu98c/VzA4eFhWkFLbpr5YdzcnL466+/CA4OrtXbNRgMPPTQQ5jNZg4fPsyHH35Ijx49ypxz5swZDAYD0dHRpKSkEBERwbx583B2dkZRFPz8/Jg3bx7W1tYl10yePJmwsDAiIyMByMrKYt26ddx7b9lpVx07dsTZWZuDr6oqOp2OiIiIWr0X0fwZqj5FlMeo1zGjf2t+O3Kec6m5HDyXwfmMAm7u3gJn26qTmldJUWg74Z8ktRuI8u29eFhSuCXyJbb8N5Rusz/Cy716QV9GRgYuGx/GoFg4p/Oj6z3vVnmNqqpsPn6BvEIzVgYdYzv7oNMpdXxDZSmKQt82bvg42fDr4QRyC838dDCB3v75DGnrUeXzOoZM4+j+L+iSvZNuR9/g/IXb8PH2rdc2CiGuARnn4L3aBWwX6QGnmlzwaBi4BVbr1LCwsJKAMioqikcffZROnToxdOjQGrcToLCwkCVLlmA0GunRowdeXl488cQT7N27l759+wKwbNkyRo8ejb+/tvCzS5cuzJ49m3nz5tGnTx86deqEopT9HWxjY8PMmTNZtmwZixYtYvXq1QwfPpwWLVpc0YavvvqKBx98kMzMTAwGA2+//Xat3oto/qQnsw7srQ1MDm5J/wA3AC5k5vNVaEydEo9fzrPzcJzm7+a0Yz8ARuVtIG3xCMLC9lfr+oPLH8Vfjcek6lAnfYSVrUOV1xyNzyQyKQeA4e09cbarh6C5Aq3c7JjRvzW+LtpQ/P7oNL49EEtOQRULnhQF3+mLycMaNyWL0ysXNFgbhRCioYSHhxMcHMyiRYsICAjgl19+YcGCBVcEeVVZtGgRLi4u+Pn58cMPP7Bu3Tratm2Lk5MT06ZNY9myZQCYzWY+//xz5syZU3Lt1q1bCQkJ4d1336Vnz554e3vzyiuvXLEwc86cOSxfvhyTycQnn3xS5h6XmjFjBhkZGSQkJPDiiy/SrVu3Gn5VxLVCkaTWtdOnTx91377SRSdnk3PYcOQ8+UXa/MB+AW4MDHSvvx5Ai5kTa5+n/bEl6BSVLNWWv7u+zA1T5lT4jN2/r2XAztkAhAXMJfieN6p8THpuIatCYyg0WQjycmB89xY1/mVXG2aLyvZTSYTHpANgb63n5u6+tHSxrfS68K9fpueJtwDYN/Ir+gy7uaGbWm+SUlLIz8uhlV/rSs/Lzs5my5YtjBo1CgeHqv9IEOJ6lp+fD2g9b4A2dJ1xrk73NJst5OTkYG9vj746ix+dW4G+en+cDxgwgEceeYSZM2eSnp7Ot99+y9y5czl27BgdOnSo1j1CQkIYPXo0zz77bLnH9+7dy+jRo0lISGDLli3Mnj2b2NhYjMYr25ibm8s333zDnDlz+Pjjj7nvvvvK3H/gwIHccMMNfPbZZ0RFRbF69WqeffZZoqKiyn12YmIigYGBxMTE4ObmVq33I5qWK36mLqMoyn5VVfuUd0x6MuvJxcTjPs7aN2HP2VS+C4urukeuunR6Ok5dRMKtX5GuOOOo5HHj0SfY/O59ZGTlXHH6hQsJBO58AoAzxvb0mPFqlY+wWNSS1d52VnpGd/K6KgEmgF6nMKKDF+O6tcDKoCOnwMz3YXFk5FU+t6n7lH8SbWgDgNufT5GbV/P8m1dTbPQZdqz+DwdeG4XTe+3xXtqT/X9+19jNEuLapTdqQ9d1+gjA4uIPbgHVO7+aAabFYuHw4cMlw+UuLi5MnToV0KYUqarK/PnzCQkJYcCAAaxYsaJWX4K+ffsSFBTE2rVrWbZsGbNmzSo3wASws7Nj1qxZdO/enfDw8CuOP/DAA7zyyivcf//96PVVF+UwmUzk5OQQHx9fq7aL5k2CzHrkbGvk9t5+9GzlAsC51Fy+Co0hNq1mqwQr07L3OKwf/otIW2344YbM74h5ZwTHTxwrOUdVVc58/iBepJGPEccZy9AZraq8977oNOLTtb9YRnf2xs7q6k/Z7eDjyLS+rbC10lNosrA1ovL0TTqjFcp4LfdnoHqO0FUvX41mVptqsXDqcCh/L3uSiFf64Le8F0MjXqNX4T6sFRNWipmW254gPS25sZsqhLjKIiIiyM3NZc2aNZhMJtLT05k/fz6jRo2iffv2nDhxgpMnT7J161Z2797N3XffXetnPfDAA7z11lv8+uuvzJ49u2R/WloaTz/9NEeOHKGoqAiTycS6des4cuRIufNCp0+fzu+//878+fOvOGaxWHj//fdJTNR+b8fGxvKPf/yDNm3a0LFjx1q3XTRfEmTWM4Nex4iOXtzUzQcrg47sAhPr9sexLyq13upt27q3JmDBHxwLmAVAN0sEPqtv4I+fV6GqKjvWf8yg3D8BONVtIV4B3au8Z2JmPrvOpGj3a+lMkGfjDcu6O1gztJ2WYikyKYfTiZXPcW3dI4Rwr4kADDj3KWdOHm3oJlbKXFTIsb9/YvcHc4h/pSPt1t3A4HMf08F8CoBcrDnkNJzwLk9ToBrxIZkTn1/5C1sIcW0LCwsjKCiI9evX4+LiQteuXbG2ti7JR9m+fXuGDBnClClTeOWVV+r0rJkzZ3L27FkGDx5Mu3btSvZbWVmRmJjI5MmTcXNzw9PTk1dffZXFixdz++23X3EfGxsbRo8ejatr+YtPf/31V7p27Yq9vT39+/fHzs6OzZs3YzDIOuPrkczJrKXL52SWJyW7gF8OJ5CSreWADPJyqHbd7+qK+nst7psexZFcLKrCFpfb6Jv+Ky5KDifsetFh4WYUXeXPM5ktfLUnhpTsQpxtjdw5wB8rQ+P+/aGqKt/ujyU2LQ9HGwN3DfTH2lDx+8jLSCH/nZ64ksk+q770eup3dHVIHF8biTEniF3/Am3TduBE2SkMSbgR5TEch+7jadf/JgzWWmnOPV8+R78z7wFwMOQzeoRMueK+MidTiOqrav5YbZjNZrKysnB0dKzWEHF1PfHEE+Tm5rJkSflp2C7OAwXo3bs3P/74Iy1btqzVs1RVJTAwkEWLFjFjxoxat1lcf2ROZhPl7mDNtL6t6ejjCMCZxGy+Co0hMTO/3p7RZvDtqA9sJ9qqHTpFZUzGWlyUHLKwo8Xdn1UZYAL8fSaFlOxCFAXGdvVp9AATtPlIozp5oy8uQ7k7MrXS822d3UkY8DwAfQr38tfPK65CKzWq2cSBNa/hsGwYvdI2lASYkTp/dvndx5mJP+Hx/Gn6PryCTsOmlASYAL2nP89JQ3sAfLY+QVZ6ylVrtxCicYWHh9O1a9cKj8+dO5dBgwbRr18/Jk6cWOsAE2DVqlUUFhZy22231foeQtSU9F83MCuDjhu7+tDS1ZatEUlk5BWxZu85Qjp40bWlU70srHHybYfDwh0cXfEwXeK/BSBu4Ct09Amo8tqYlFwORKcB0LeNW5Wrua8mN3sr+rRxJTQylbCYNDr5OOLlVHHvROexs4k4uIoOeWF0CHuF5CET8HCvea31mkg4FUbW2ofoVXgMFEjFiWNBs/EffAeBgZ2oKlOe3mDE+raPKVh9A95KCvu+eJQ+j65q0DYLIZqG8PBwnn/++QqPr1y5sl6e4+npicFgYNmyZVhZVT0/X4j60vhdVtcBRVHo7ufC1L6tcLI1YrJoyc43Hr1AUT1VCdJZ2dLlgWXET/iGmLHL6XjD/VVek19k5vdj5wHwcrJmQGDDBmS10beNGy52RlQVtpxIxGKpZHqHouA+9X0KVQPepHJ01VMN1i5LUSFhK/+F+8rRtC/UFl3tchiD+aFQhtz1Aq0CO1Vxh1L+HXuxP/BBAPqk/syx7dWvOSyEaL6SkpJqnXS9ps9JSEhg3LhxDf4sIS4lQeZV5O1kw8z+l9btzuTrPTWr210V3+CxtB44GarRQ7o1IpGsfBMGncLYLj7o67mqT30w6nWM7OgFwPmMfA7HZVR6vkebrhwL0lZODkn5lvDQrfXepvhju4j5Tz+CTy/BSjGRgAehgz5mwIK1eHrVrupQ3xnPc8Kg5cTz+PMJcjMrnx4ghBBCNHUSZF5lF+t2D25btm73yQvVq9tdX05eyOJ4gvbMwe088HCwruKKxuPvbl8yr/Wv08lkV5F7tNu0F4nT+aJXVKw3LiS/oH6CeHNBLgeXz8drzTjamM4CsN1lItaPhNL/hml1mvpgNBoxTP6IAtWIl5pCxOeP1EubhRBCiMYiQWYjUBSFfgFuTOnlh11xPshfDiXwZ0Qi5sqGg+tJdoGJLce1PGat3OwILs7r2ZQNa++JtVFHocnC9pNJlZ6rt7Kl8MY3AehkOcX+d25j05evE/rnT5yLi618yL0CseGbufBGb3pEr8CgWIimBXtCVjLssc9xc/eo1Xu6XNvOvQhtMw+A4JSfOfW3JGkXQgjRfMnCn0bUys2OmQP8+fVwAnFpeYTHpHMhI59x3VvgZNMw9cJVVWXTMa38pbVRxw1dvK9aVZ+6sLc2MKStB1uOJxJxPovOLZxo42Ff4fkB/W7h0O6xdE/dyOD8bXBmG5wBtkGy6ky8dRtynNqi9+6Ei383/Nr3ws7F64r7mHLTOb5yId3i12qvVR3bPabR6+7X8Xd2rvf32X/m8xz7z+90Nkfgsnkh+d1DQJGJ+kIIIZofyZNZS9XJk1ldFovK32eS2RelrfK2tdJzYxefSoOo2jp4Lp0/Tmi9mDd29aFTC6d6f0ZDUVWVNXvPkZCRj7OtkbsG+mOsJBemOT+L42tfwiphPx55kbip6ZXeP1VxJskmgAKXdhhbdMHK1gGnXf/B06L1nJ5U2pAx5h36DhpZn2/rCscP7yPw2xuxVoo45HkLgfd8LHkyhaim5pQnU4jmoC55MiXIrKX6DDIvOp2Yze/HzlNQZEFRoH+AO/0D3NDV04Kc1JxCvgqNpsis0t7bkXHdfJpFL+alErPyWR16Douq0j/AjUFtqz9UXZiZTPypMFKjDmO6cAy7jFN4F0ThSXql1xWoBrb63MvAu1/Gyd6u0nPryx/LnmHkufcBOBHyCREZNhJkClENEmQKUb/qEmTKcHkT0tbLAU8Hf34+HE9iZgG7I1NIyMjjxq4+da4jbraobDx6niKzioO1gZEdvZpdgAng5WhDcGsX9kensS86jQ4+jrhXc9GSlZMHbXqPoU3vMSX7VFUlKfE8safCyIg5AokncMw6jZ8pGm8lnXBdZ8zj3mFsnwEN9ZbKNejOFzj6n010sUTgvu1plC6vXtXnCyGEEHUlC3+aGGc7I1P7tKJbS22+X3RKLl+FxhCfnlen++45m8r5DO2vkTGdvbG1ar5/jQ8IdMfRxoDZorLlRGKdasIrioKndwuCh4wjZMaThDz2Gb2f247zs5GcmRdJ53/9Te+rHGAC2FhbYb51CQWqEU81BY9Tq696G4QQ4qKoqCgURSE2NrbBn5WamsrYsWNxdnamd+/exMTE4ODgQHx8fIM/W9QvCTKbIINex+jO3ozt4oNRr5VVXLsvlgMxabUKqM5n5LPnrJZ3sWcrlwaZ63k1WRl0jCjOnRmXlsexhMx6f4aNUU+Qj3ujltjs3rMvW1s+AMCggu1c2P9Do7VFCFH/QkJCsLa2xsHBAWdnZ4KDg1m3TooxfPTRR2RnZ5OSksL+/ftp3bo12dnZ+PpqeYhXrFhB27ZtG7mVojokyGzCOvs6MbVva1ztjFhUlW0RSfxyOIECk7na9ygyW9hwJAGLquJqZ2RIu/pJt9PYgjwdCPLS5ifuOJVMXmH1vybNyZC7XuCwoiVp9/3zcVJP/NXILRJC1KfnnnuuJKCaPn06U6dO5eTJk43apqKiokZ9fmRkJJ06dcJgkBl9zZ0EmU2cp6M10/u3pr23loz81IVsVofGkJRVUK3r/zqVTFpuETpF4cauLSpdjd3chHTwxMqgI6/QzI5TlefObK7sba3JvXkJCaobthRgXDOVrJjDtbpX2A+L2f+/aSTHnannVgoh6spgMPDQQw9hNps5fFj7Gf/www/p0KEDzs7ODBgwgB07dgDaQiRXV1d27doFaEGZoii88MILJffr1KkTa9dqqddyc3NZuHAhAQEBuLm5ceONN3L69OmSc0NCQnjssceYOHEiTk5OvPXWW1W29+DBgwwfPhwPDw9cXV256aabOHNG+92SnJyMXq8vGd7esmULiqKwfPlyAEwmE05OTuzdu/eK+44fP57PP/+czz//HAcHB1544YUyQ/W7du1i3rx5REZG4uDggIODA1u3bq3pl1tcJddOxHENszboGdfNh5AOnugUhbTcIr7eE8PR+MpLLEYl5xB+Lh2AfgFu+DjX32rLpsDJxlhSb/1ofCaxabmN3KKaMZkt5BdV3QPbpWMnfmj5JKmqA45qNoUrJpCfeLb6D1JV9i9fSHDYs/RO+43CZTeTmRhTh5YL0QydOVP5R8Zlv0/Pni1zXHfZa9LSyp4fFaXtr6XCwkKWLFmC0WikR48erF69mueee44vvviClJQU5syZw4033kh0dDR6vZ6QkBA2bdoEwKZNm2jbtm3J6/j4eE6ePMnIkVq6tdmzZ3PixAl2797N+fPn6d+/P7fcckuZHsvPPvuMRx99lIyMDB599NEq26soCi+++CJxcXFERUXh4ODAnXfeCYCHhwfdu3dn8+bNAGzevLlM+0JDQzEYDPTu3fuK+/7000/MnDmTe+65h+zsbF566aUyxwcOHMhHH31EYGAg2dnZZGdnExISUsOvtrhaJMhsJhRFIbi1K7f38cPRxoDJovL70QtsOnaBIrPlivPzCs1sOnYBAB9nG/oFuF3tJl8Vwa1c8HTUVpf/ceLqVEyqi/wiM8cTMvnpYDwfbTvD0u2RnEutOjhu4e1LaP8l5KjWuFtSyPhkHEUZ56t+oKoS/umD9I5eWrLL15JA9sc3kp8aV5e3IkTz0rZt5R8rV5Y9v2fPkmP6Dh1w6tULfYcOped/+GHZ84cO1fbX0KJFi3BxccHPz48ffviBdevW0bZtW5YvX87cuXPp378/BoOB+++/n+7du/PVV18BMHr06DJB3D//+U+OHz9ORkYGmzZtomfPnri7u5OcnMzq1av54IMP8Pb2xsrKihdeeIGEhARCQ0NL2nHbbbcxcuRIFEXBzq7qVG3du3dnxIgRWFtb4+zszAsvvMDu3bvJyckpt32vvvoqW7ZsQVVVNm/ezIgRI9DpJAS51sl3uJnxdbFlZn9//N21XwJH4jJYs/cc6bml9blVVeWPE4lkF5gw6hXGdvFBX0+5NpsanU5hdCdvFAVSsgvZH51W9UVXWXaBiYPn0vnuQCwfb4tkw5HznE7MpsisYrKobDp2gULTlX8oXG7osNFsDX6XAtWAtymeCx/cjCU3vcLzVbOJwx/eQ884bWX6HtuhbO35DoWqHl9zHGkf3oipOoGqEKLBPPPMM6Snp5OYmMjOnTsZP348AOfOnSMwMLDMuUFBQZw7dw7Qgrjdu3eTlZXF1q1buemmmxg0aBB//vknmzdvZvTo0QCcPauNenTv3h0XFxdcXFxwc3OjqKio5F4Abdq0qVG7z5w5w+TJk2nZsiVOTk4MHjwY0IbKL7Zvy5YtpKWlcfLkSSZPnoyHhwcHDx4s0z5xbZNZtc2QrZWeiT1bsicqld2RKSRlFbAqNIaxXbxp6+XIifNZnLyQBcDQdp642V/bZQl9nG3o7ufMwXMZhEam0N7bARe7xn3PaTmFnE7K5kxiNgnFqaMu0usUWrnZ4udqx87TKWTkFbHzTDIhHa4sa3m5myfO4LucNCacfAa/gtNEfzCB1o/+hmJVtudBNRVw4sMZdEvRehL+chhLn0e+xMbaml9NKmMOP0GLohjOf3Aj3o9sQnHwrL83L0RTdMkcxHJ5XLYoMjwcLNoff2azmZycHOzt7UuTsbtdNjq0YweY628BYqtWrUoCxIsiIyNLgtAOHTrg4+PDu+++i7e3N76+vowePZpNmzbxxx9/sGLFCgD8/f0BOHXqFJ6eFf+c17RXcd68efj6+nLo0CHc3d05cuQI3bp1K8mAMmzYMFJSUnj//fcZOnQoRqOR0aNHs379ekJDQ/nss89q9Ly6tFU0HvlONVM6ncKAQHcmBbfE1kpPocnCTwcT+OPEBf6M0MpGtvGwo7tf/dfXbooGBXlgb63HZFH5M6JuuTNrK7/IzM4zyXyxK4oVO6P461RySYBpZdAVV1lqwQPDApkU7EffNm70DXAFIPxcOnHVzIU6acZDrG+5EAD/7HCiPp4KZlPJcbUwl1OLJ9KpOMD8w3kS/eZ/hY21Nq3gpin3sT7oFUyqDp+CsyR9cCPkptbb10GIJikoqPIP58t+VwYElDluuew1rq5lz2/TRttfT2bNmsXHH3/Mnj17MJlMrFixgvDwcKZPn15yzqhRo/jvf//LmDFjSl6vXLmSlJQUhgwZAoCXlxczZszgoYceIi5OmyKTnp7O+vXryc7OrnX7MjMzsbe3x8XFheTkZJ5//vkyx21tbRk4cOAV7Xv33Xdp0aIF7dq1q/WzfXx8SExMJDOz/tPXifp13QWZiqK4K4qSpijKnY3dlvrg727PjP6t8XXRFvUcPJdBQZEFG6OeMZ2bX9nI2rIx6hneXusJjErO5VRi7X951kZkUjZf7oomNDKVlGxt6oKdlZ6uLZ2Z0NOXucMCubl7Czr4OGJjLE2E3z/AHQ8HK1QVNh+7gKmc+bWXUxSFiff/i+/cZgMQkLKdyM/uBYsFS14mZ9+7mfYZOwHY4DaToY8sw8poKHP9bXc+xKqWz2JWFbxyT5P84U2Q1/SmGghxvZoxYwYvvPACd955J+7u7nzwwQf8+uuvZYa1x4wZQ2ZmZkkQ161bN2xtbRk8eDC2trYl5y1dupQOHToQEhKCo6Mj3bp1Y+3atXX6/8M777zDjh07cHJyYujQodxyyy1XnHN5+0JCQsjNza3zUPnIkSMZM2YMAQEBuLi4sG3btjrdTzSc6652uaIo7wBtgTWqqq6s6vyKNETt8rowW1R2nEoiLCYdgJu7tyhJe3S9UFWV9WFxRKfkYm+t5+6BbcoEdA0hr9DMtpOJHE/QpicYdArd/Jxp6+WAr7NtterOn8/I5+u9Magq9G3jdkUu0+zsbLZs2XJF7fL8QhO//+8Bbs3RkjdHBc5Adz6c1rnHAPjJcw43zXsDQwVpqwpNFpZ9+B/mJv8HnaKS5tIN13m/gE1xj46pEPLTteAzLx0UHbTsBbrmWy1KXPukdrkQ9Utql1eToihBgDuwv7HbUt/0OoWQDl6083bEZLbg7968q/rUhqIojOzoxZe7oskpMLPrTEpJZaCGcOpCFn+cSCS3OBF8SxdbxnT2xrWGc2B9nG3o7e/Kvqg09ken0c7bAW+nqv8HaWNlIOThj/j9nQxuKNxMm8ivSo591+JxJsx5odIFX1YGHffMfYLFiwt5JPMdXNMPk/9OL/QGI/qCDHSmK1e9q44tULpPhZ4zwLNDjd6nEEKI60uTGi5XFGWaoig7FEXJVBTFVM5xvaIobyqKkqQoSpaiKOsURalJCZtXgJeqPKsZa+lie10GmBe52FmVpGs6GJteUq+9PuUWmvj5UDw/H0ogt9CMUa8Q0sGT2/v41TjAvGhAoHtJZaffj12odiomJ1srgv/xBX/p+wFgVhW+afUME6sIMC+yszJw97ynedvmHwDYFCRjzEkoN8AEULIS4O93YUk/0t8bRvZfH8kwuxBCiHI1tZ7MNOADwBb4pJzj/wQmAP2BFOAz4EvgJkVRHICt5Vzzk6qqLymKMhBIUVX1zPUyT/F61aeNGxEXskjJLmTLiQtM79u6WsPWVVFVlYgLWWyNSCopY9nKzY4xnbxxtjPW6d5GvY4xXXxYu+8cyVkF7DmbysAg92pd6+lsT/7cNXy15m2s/Xpw+8TbajTXytXeiukPPssjH/vhm3mQDOzJUO2LPzuUvPZS0pii38Fk/Q58lDRcUg/C5oMUbH6WE9630HX2J+iN13YmAyGEENXXJOdkKooSAmxWVdVw2f5o4GVVVZcVvw4CTgMBqqpGVXHPR4EpQB7anMwcYJ6qqrtq0C53tOF2unXrFrFz587qXiqusrj0fL4/qOWADPK0w8/FFnd7Ix4OVrUqrZlTYGLbqVTOpmg9fFYGHQMDXOnSwqFeF1ftOJ3Cobgs9DqF24Jb4OFgRW5uLrt27WLgwIHVSpJcW2aLyoXMAhQFdIqCTtGmIFx8nZFXxKG4LA7FpqGL+osBWRu5QbcXG0WrGnIgYB7tpzzXYO0TojoKCgqwsbGp1zmZFoulJIWRpM8R15v8/Hzy8/OxLs5QcjlHR8cK52Q2myBTURRnIB0IVlU1/JL9GcBdqqr+WIP7vwicrunCn+LrXgBwdXUtqcMqmqZDqQpxOZcHgCp2BnAygqNRxdEKHI1gq4fyYkVVhbhcOJ6uYLJoJ3jYqHR1VbFtgHEAkwX+uqCQZ1JwtlIZ4KXSVPPoF5jhfGYenWO+ZLTlL4pUPZvavkSRU+vGbpq4jhkMBvr27YuVlfSqC1EfCgsL2bt3LybTFbMYAZg4ceI1EWS2AmKAQFVVz16yPxp4pi4rxWvQLunJbEaKzBbCYzNJyCggOaewZIi7PNYGHe72Vrg7GPGwt8Ld3gprg44dZ1KJSc0rOWdIkBsdvO0bNDXUubQ8fjyklQQdGOhKR3fjVenJrK3zySlYLx9NKyWRBJsgHB/cDHr5H7xoHCaTCUVRymRiqCvpyRTXs6ys4uwphvJ7VirryWxqczIrk1X8+fLs4i7AVcnIqqpqCtpcUPr06VOvv8REwxjh7FSynVNgIimrgOTsgpLPqTlFWFQVFUjOs5CcV0BEUkGZe1hZWRHoac+oTt44WDf8j0wnBwdisywcicvgYHwuAe7a2jY7O7sm+W+urYMD63u+SsvwubTIP0PS1rfxnPBqYzdLXKdUVSU9PZ28vDz0en29/EFoNpspLCzEaDRKCiNx3VBVFbPZjKqquLi41OpnqdkEmaqqpiuKEgP0AsIBFEUJBJyAQ43YNNFM2FsbsLc20MajdPW9yWwhNaeQpJLAs5CkrALyi7ReTxujnhEdPeng7XhVE9sPbedBdEoOWfkm/jyZjFPTG3AoY9z42/j26E/cYfoJt7AlqL0nofj1buxmieuQoii4uLhQVFSExVJ1cYPqyM/PZ+/evQwYMAB7++s3e4e4viiKgpWVFfb2tR+9a1JBpqIoesAIWBW/vjhzu0DVxvU/AZ5SFOVPtB7F/wAbq1r0I0RFDHodXk42eF2Sl1JVVbILTKTnFuHpaN3gCd3LY2PUM7KjFz+Ex5OQUUDa1S1gVGPWBj2eExZxZu0+gnQJZK2ZjeOju8BYf4svhKiui/9zrC8mkwmTyYS1tXW9LigS4lrX1CaX3IW2+nsjoC/ezgP8i4+/DvwE7AXiis+5JspDiqZDURQcbYy0crNrlADzokBPBzq10Ko2nczQVnc3ZSO6+fOV7z8xqwqOWZEUbn6lsZskhBCiETWpIFNV1RWqqirlfEQVHzerqrpQVVUPVVUdVVWdrKpqciM3W4gGM7y9F7ZWesyqwm/HkkqG8ZuqO2+7nU8ttwJgDF0CMbvLHFfNJg4ePcrGvUernXBeCCFE89SkhsuFEGXZWukZ1cGd4xGnSMku5MfweCb1almrXJ8xKbnEpecR3NqlwXpoAzzsyR64kBO7D9BRd46itbMxBg2jKDWavKQo7PIS6IGZrqrCb7vu4YYH38LKKL+GhBDiWtSkejKFEFfyd7Ojq5vW6xeXnsevhxNq3At4ND6D78Ji2R2Zwo/h8RSZ62dBRHnmjerMv60epUjVY8w6B+GrMMb8hVNeLAa0nli9onJL6gqOvD2evCwpSymEENciCTKFaAb87LWcmQCRSTlsOnaB6ua4PRKXUXy+9vpioGppoOFqe2sDU8bfwr9M97Pf0o7fzH351HQTiyz3sLzVIo5P+JkjXuMB6JW3k9R3h5Add6xB2iKEEKLxyDiVEM1Er1bOKAYr9kWlcTwhE1srPcPaeVSaWuLguXT+OJEIgI+zDf5udoSeTSUyKYfNxy8wprN3g6RmGt+9Bd8duJ0pESF08HZkRv/WPBzcEmfb4hrvPYfw9+rX6RfxJi3NseR+OpKsCR/j2HMCZJ2HhIMQH659tveAcW+CofySZkIIIZomCTKFaEaGtPUgr9DM0fhMDkSnYWvU0y/ArdxzD8SksS0iCQBfFxsm9GyJtUFHodlCWEw6R+Mzsbc2MLitR723U1EUPrmrDxcy8/Fztb0ykFUUBs94ms0bOtBz13w8yITv76ZggyfW+UlX3M/i0wNdv/vrvZ1CCCEajgyXC9GMKIrC6E7eBHlplX/+Pp3M4diMK87bF5VaEmC2dLVlYnBLbIxa9ZPh7T1LUiPtOZvKgZiGmRNpZdDRys2u0p7S0TdOJuymHzhoCQIoCTDNqkKExY9TlpYA5G19G8zl180VQgjRNElPphDNjE6ncFNXH9aHxRGXlseWExewtdLR1ksLHHdHprDrTAoArdzsuLWHL1aG0r8nFUVhTGcf8orMRCXnsi0iCTsrPR19nMp9XkMbM6AX2xx/4JVv38dUmMcZY1uS7Nth5+CEQ+Ypvsyfj31uLObD36LvOa1R2iiEEKLmJMgUohky6nXc2sOXb/fHkpRVwK+HzzMpWM+51FxCz6YC0MbDjlu6+5ab7kivU7i5my/fHYglISOfjUcuYGPQlym5eTUN79KKwR3/jUWlTEB8JK4rmz76gjH6/WRv+S/O3e8AnQzACCFEcyC/rYVopmyMeiYFt8TFzojZovLdgbiSADPQ057xFQSYF1kZdEzo2RI3eyssqsovhxOIS88jv8hMfpGZvEIzuYUmcgpMZBeYGjwRvEGvKxNgAnRt6cyB1vcC4Jx1ClPEb1deeP4IvNURvpwM5qZdFUkIIa4n0pMpRDNmb21gcrAfa/bFkFOgBYFtvRwY160Fel3Vq8ZtrfRM6tWSb/aeIyvfxDd7z1V6vpOtEV9nG3ycbWjhbIuno3W1nlMXE8dPZOcHHzBId4z0ja/j0XEcXJznmZ0Eq6dDVoL2sf1NGPGvBm2PEEKI6pGeTCGaOWc7I5OC/WjhbEPPVi7VDjAvcrIxMim4JXZWVVcByswr4sT5LLZGJLF6Twwf/Hmab/ae4+C59Grn7aypDj6OHGh9HwAe6YcoPLNdO2AqhG/uhoyY0pO3/xfO7W2QdgghhKgZ6ckU4hrg6WjNtH6ta329u4M19w4OICWnAAUFRQEFiv8DCgr5RWYuZOYTn5HP+Yw8cgrMmCwqcel5xKXnYWPU08HHsV7ez+VuHD+NQ0s+oLsukuTfXsf34WHw60KI2amdcOti2LkYkk/Cd3Ng3l9g7dAgbRFCCFE9EmQKIQBtjmYLZ9tKz2nlZgeAqqpk5ps4n5FPWEwaCRn57DiVRICH/RXzKutDW29Hlre6l+5xz+GbshPTj49hCPtcOzjkceh1N/h0h09HQdpZ2PgvuPW9em+HEEKI6pPhciFEjSmKgrOtkQ4+jozt4oNep5CVb2JvVGqDPXPExHs5rfoCYAhboe1sfyOMfE7b9u0JIU9r2wc+hxO/NlhbhBBCVE2CTCFEnbjaWxHc2gWA/dFppOcWNshz2ng6st/vnpLXRa7tYPJS0F0yl3TI49Cqv7b94yOQFt0gbRFCCFE1CTKFEHXWL8ANB2sDZovKtpNXloWsL4MmPsh+SwfOWTyZkPoIP5zIKnuCTg+TPgYrB8hNhg8GwI63wFRQek5sLLz7LowcqX08/DB89FGDtVkIIa5XMidTCFFn1gY9Q9t78Nvh80Qm5RCZlE2gZ/0vvGnl6czpmT8zZ004qblFzP86nG0nk3h5QlccrIt/nbkFwG2fwfp5kJcKW16GfSshcxhs2g9//132pn/+CYGBMG9e6T6TCQzy61EIIepCfosKIepFB29HDsVmEJeWx7aTSbR2s8NQSTL42hrRwYvfHhvGgm8O8tfpZL47EMf+6DQeCgnC2dYKB2sDDrb98Zn1Nz773oJ9n8HXx2BPeOlNnJxgwgRwdYVjx6D1ZSvzH38cwsNhwQLtvErqrwshhCifBJlCiHqhKAohHTz5KjSG9NwiDsSk0y/A7Yrz0nMLyS4w0dLFFqWWwZu3kw1f3NePpTsieXNjBNEpuTy17vAV5y2evoDxve4Gu0fg1F/gpYdbBsGLP4KDc/k3z86GL76AzEz46y/o3RtefhluukmCTSGEqAGZkymEqDdejjZ099OCtz1nU8jM18o8mi0qEeez+HZ/LMv/jmLtvli2RiTVKYG7Tqcwd3gQ3z00iD7+rng4WGF9Wfqkr/fGQIsesHAbfLoAptqCQzh8dzcU5pR/Y1tbWL0aRo3SXu/fDzffDIMGwfbttW6vEEJcb6QnUwhRrwYFeXDyQjZ5hWb+PJGIm70Vx+IzyS0sW/s8/Fw6BSYLN3T2RleH0pTd/Vz49sFBJa+LzBZ+3hdN3OP/5LOCyaTn9sLFzgomvQHuXvDnqxC5Fb6cBDO+AVuXsjfU62HcOO1j2zZ47jnYsQN274bhw7Xh86VLwdOz1m0WQojrgfRkCiHqlY1Rz6AgdwAik3LYF5VGbqEZRYE2HnaM79GCri213s7jCZn8fDgBk9lSb8836nXcsvQ1Ht71DV9+9S+27j+rHVAUGP4EjP239vpcKPzHHxa1gDfbwnvB8MercGnv6vDhWqD5++/Qq5e27/BhbU6nEEKISklPphCi3nX1deZYfCYJGfnYW+vp4utMV19nnO2MAAR5OmBt0LE/Oo0zidn8EB7P+B6+9VMtaMMGjMs+BWCnf3f2RmYyceglxwc+xIUCAx5/PoleUaEoV/vISYLtb4KigxH/Kj1fUWDMGG34/KuvtMVC1talx3fsgD59tGF2IYQQJSTIFELUO51OYVKvlqRkF+LtZIP+suFwRVEY2s4Da4OOnWdSiEnN5bsDsYzs6EWByUJ+kZm8IjO5hWacbY10alHNnsOMDJgzB4DkrsH8O+RerE4nk1doxtaqNGn728kD2FHwP4J08dzdy50xQfYQ8Suc+Bm2/QfcgqDH1MvfFNx5Z9l9cXEwdqw2dP7aazB9unaeEEIICTKFEA3D2qDH16Xi3j1FUegf6I6VQcfWiCQSMvJZFRpT7rlWBh1B1cm7+cQTWrJ1a2v0y5fDuljyiyxsO5nEjV19AEjLKeT78DgK8CDe4sGho0a23zwC5263wRcTIWYn/PgwuLQG/4GVP2/rVigshJgYLQB95x146y1tmF0IIa5z8ie3EKJRBbd25YYu3hgu6e006hWcbI0lvY+7zqRUvRJ982ZtQQ7Aiy/i2qcHffy1FEq/HztfctqafecoMFmwNeqxMerIyCviw61nwGANU1eCawCYC+HrGZBypvJnzpwJR49qi4FAW4keEgITJ8LJkzX5MgghxDVHgkwhRKPr4uvMnGGB3D80gIdHtuXhke24f0gA43v4ApCUVcCZpOyKb5CVBbNna9t9+sDChQDc0MUbgC3HEykyWzBbVL7cpdUzn9yrJbOHBAKw/O+zxKfngb07zFwLNs5ataClIyF8ddnFQABZ56GguKRlhw7w/fda5aDevbV9P/wAXbrA66/X/YsjhBDNlASZQogmwcaox8nGiPGSKkEtXWzxd7cDqujNTEoCZ2cwGmH58pKSkGO7aEPkGXlF7DmbyubjF4hLzwPgnkFteGB4IK52RgpMFt7dXNzz6NEOpq4CayfIT4fv58Gq2yDhIIR+Ap+Ohrc6wOLecOFoaRtCQmDPHvjyS/Dz00pT+vvX69dICCGaEwkyhRBN2sDidEjJ2YWcvFBBb2ZgIOzdqw2Zd+1asruVmx2dixcNbTx6ns93RgEwKMid9t6OONkYeWRkOwC+3R/LyQvFvZMBQ+Gh3dBurPb69Gb4eBj89gTE7tX2ZV+AFTdD3IHSdlxcHBQRAR9/DNOmlR47fx5+/PHKXlEhhLhGSZAphGjSWjjbEuhpD8DuyBQslgqCNCsrGDbsit0XezPXh8Wx80wKoPViXjRzQGtaudliUeHfvx4v7S11bgkz1sCUZWCnBboY7aHbHXDr+2DrCnlp8PmtEL2r7EPt7OCBB8qWoXz+eW3u5ujRWl10IYS4xkmQKYRo8gYGakFeak4hJ85nlR74/Xc4fbrSa8d21eZlZuWbAG0IfnQn75Lj1gY9C2/oAMCfEUmsD4srvVhRoNtt8PA+uHMdPHEKpiyFXnfBrF/B3gsKs+Dz8fBeL20o/atpcHZH2UaYTBAZqW3/8YeW2H32bEhIqM2XQwghmgUJMoUQTZ6Xkw1BXloKo9CzKZgtKqSlaau7O3fWao1XoIO3Y8m8ToC7B/pfkbfz1h6+DG+vlYl84YejxKbllr2JnRu0HQ1W9qX7vDvDfRvAyQ8sRZB6RhtKP/kbfDcHivJLzzUYYNMmWL8e2rbVhsyXLYN27WDRIsjLq+VXRgghmi4JMoUQzcLF3sz03CKOJ2TCSy9BcrIWwA0eXOF1iqJwQ2et59LGqGNq31blnvPmbd1xtTOSVWDi/745qAWyVXEPggf/gslL4cbXYegCUPSQlQBhX17+EC210dGjWj5NFxfIyYFnn9VWqJ89W90vhRBCNAsSZAohmgVPR2vaezsCEPFHKOr772sHnnoKWreu9No5QwMJ6eDJy7d2xcXOqtxzvJxs+PfkbgDsOZvK0h2R1WuYrSt0vwMGPAijnte2Af56B0wFV55vZQWPPaYN8z/6qBYke3rKSnQhxDVHgkwhRLMxINANBZU+776MYjZjatUannyyyuu8nGxYcW8/7iinF/NSN3Ztwe29/QB46/cIjsZn1LyRw57Q6p9nxl3Zm3kpd3f43//gyBFt6PzScpQrVkBUVM2fLYQQTYgEmUKIZsPdwZrJsQfwP/A3ABtmLWBHbHb1hrar6YVbu9DKzZYis8rT3x2uutLQFY0Mgm63a9s7KujNvFSHDtCzZ+nrY8e0RUEdO8LTT0NmZs2eL4QQTYQEmUKI5iM/n9avPAPA+V4DODVkLPui0liz9xyRSdnEp+eRnF1AZn4RRWZLrR7hYG3gP5O7A3AoNoPNxxNrfpOS3sxYCF9Vs2vPnQMPDygo0CoGtW2r5dw0mWreDiGEaEQSZAohmo9167QFMjodHp9+SK82Wm3yC5n5/BAez5q95/hyVzTLdpzlw61nOJucU6vHDGrrUbLQ6N3NJ2vem+nRDrrepm1v/y/k16A3cuxYOHUK/vUvsLbWqhnNmwfBwVrKJiGEaCYkyBRCNB8zZmhVc155BUNwT4a392RScEtc7IxXnGq2qByKTa/1ox4f0x6Ao/GZ/H7sQs1vMPwp0FtpczM3Pl2zax0dtdRGERHaewZt7ubYsfD44zVvixBCNAIJMoUQzYeiwPjxWi9fsTYe9tw7OID5o9rxYEgQ9w8NYFhxzsvolFwKTOZaPapfgBuD22q9mf/bfKoWvZltYeSz2nbYSojYcOU5Fgv89S6crKCH0t8fVq2CXbtg4EBt380316wdQgjRSCTIFEI0fbm5VZ6i0ynYGPU42Rjp4uuETlEwW1SiU8q/VlXVKgPHx0drvZnHEjLZePTK3kxVVflg62nu+GgXMeU9Z+DD0Kq/tv3To5CbWvb4kW9h8wvw9QzIiK24IQMGwN9/w59/amUpLzp1Cv77X23+phBCNDESZAohmjazWatJPm2atiimGmyMevxcbQE4nZh9xfGMvCI+2R7J2v2xFddCB/q0cWNoOw9Am5t56bmFJgsL1h7kjQ0R7IlK5cvdUVfeQKeHiR+C0Q6yL8CvC8seP/q99tlSBDsXV/6mFAVCQsrue+IJ7aNLF/juO62SkBBCNBESZAohmrbly2H/flizBs6cqfZlF8tQnk3OwXTZSvMDMWnkFpqJS8vjaHzli3IeG90OgBPns7jrs1C+2XeOuPQ87luxl+8OlNY5Dz2bWv4N3INgzMva9pF1cP6wtl2QBac3l563/3PITtK2C3O0+ueWS4b6U87A8nFw4AvttckE3t5afs0zZ2DKFBg+HPbtq/T9CCHE1SJBphCi6UpPL51/edttV/bkVSLIU6szXmiycC6ttDZ4gcnMsUsCy92RKRSaKk531NvfjZu6+gDw9+kUnvz2EINf/4O/TicDlMzbPBKXQVZ+Ufk36XMfuLbRtkM/1j6f+h3MBVoZSmsnMOXB7iWQmQBLR8Lnt8D2N0vvseVliP4bNjytBagGg5baKDwcxozRztmxA/r2hXvugdhKht+FEOIqkCBTCNF0vfyylsLHxkabe1gDjjZGfJxtADhzyZD58YQsCk0W9DoFvU4hu8BEWExapfd6b3owi6cHc0Nnb6z02q9NnQKvTuzKBzN6oyhgUWFfVAX30emh3wPa9uG12tzMYz9orwOGQr852vaeT2H5jZB0Qnu9833IS9Pmax7/SdtXmK3d46Ju3WDjRvjlFy2BO8AXX0D79jXq+RVCiPomQaYQomk6cgQWF89TfPLJWtX2DvLUhszPJGVjsWgLfQ6eSwegvbcj3fycAdgXnUZeYcWr0I16HeN7+PLJ3X3Y99xoFk8PZt2Dg7hzgD/OdkY6t3ACtF7RCvWcCUZ7MOVD6EdwapO2v/MEGPCQNm+zMAvSokBvDQYb7XXox7D3U1Avad++z8rOv1QUGDcODh2C99/XSlYOHQqBgTX+mgkhRH2RIFMI0fSoKjz0kDbv0N8fnnqqVrdpWzwvM7fQTEJmPjGpuaTmFAIQ3NqF/gFuWBl0FJoshJ6tJEC8hJONkfE9fAlu7Vqyr3+ANmReaZBp6wI9p2vb2/8LRbmAAh1vAXsP6H2vdsxoDzPXQp/7tde7P4D9K7TtNkO1z+cPQ9yBK59hNMI//gGnT8PSpVrwedHzz8O2bdV6j0IIUR8kyBRCND1ffqnNLwStN9POrla3cbO3ws3eCtCGzMOLezFbONvg7WSDnZWBPv5asHgoNoOM3ArmVFZhQKBWeehIfGbF8zKhdMj8Yq+k/2Bw8NK2Rz4LN7wKszdD4HAY/KjWo5mfoQ2Z6www6WPw6KCdv++zip/j4gKtW5e+/vtveOUVbU7r5MlaECqEEA1MgkwhRNPTqxcMGQK33qolX6+Di72ZxxMyS8pM9mztUnI8uLUr9tZ6zBaVnWeSa/WMfgFuKIpWZWhfdCXzOz07QOCI0tedby3dtrKDQY+Ad2fttaMP9L7nknMngHNLbRERaCvV8yqfS1p6byvtawqwfj107gwLFkBaNa8XQohakCBTCNH0dO2qDe1++WWdb3VxXmZuoRlVBXtrPe28HEuOWxl0DCiuU37ifBZnkq7Mq1kVFzsrOvpUY14mQP952mdFD52qCKAHzweDLaDAgH9o+3pM1faZ8mDjM9o+iwV2vAW/P6ttX65vX9i7F1asAF9fKCqCt9+Gdu20OZxFtevBFUKIykiQKYRomnQ6cHKq8228naxxtDGUvO7W0gW9TilzTldfZ3xdtJXom49dILfQVOPnXBwyD42sIF/mRe3Hwq3vw7RV4ORb+bnOftrw+X0bwK+3ts/WFYY/qW2Hr4IDX8JPj2gpjnYuhsg/y7+XTqelNjp5El54AWxtISUFHnkE5s6twTsVQojqkSBTCNE0mM1w//1aj1s9UhSlpDdTr1NKVpRfSqdTGNvFByuDjtxCM5uOXahxrfKLi38Ox2WQXVBJkKoo0Osu6HBT9W7s0xVaDyi7b/Bj0La4vOSPD2u10S+K3lm6XZBVNqE7gL09vPiiFmzefbcWfD78cPXaIoQQNSBBphCiafjoI/jsM61O96lT9Xrr4NYueDlZMzDIHQdrQ7nnuNhZMby9JwCRSTkciau8EtDl+gdoPZlmi8ovh+L596/H6fXKJu5bsbfSZO+1otPBpE/A8ZKeUFvt+SVBZkwovO4PX029cgj9yHew40lY8l+IjCydrwmwZw888ABcuLJWuxBC1IQEmUKIxnf+PDxTPL9w+nRtrmA9crGzYmZ/f/q2cav0vC6+TiXlKLedTCStON1RdbjaW9HRR5vr+dS6w3y8PZLUnEL+OJHIa78er33jK2LvDlNXQoseMOxJGFdcHShuHxTlw/7l2ir205tg/yUr0Yvy4afHIOIX2Lu0bP5RVYX/+z8t/VHbtvDvf0N+fv23XQhxXZAgUwjR+B5+GDIytDmYNazsU58URWF0Jy/srfUUmVU2H6/ZsPnAIPeSbRc7Y0nJyRU7o1gf1gBlHv16w9ztMPIZLR0SgLkQzu2GiF9Lz9v0ImQU11k/+RsUZGjb8WFl73chAjrbgpMjZGdrJT07dtTqxtdw+oAQQkiQKYRoXOvWaR8Ab74JPj6N2hw7KwMjO3oDEJuWR+wldc+r8vCItszo35oXx3dm5z9H8vm9/RhYvHL96e8Ol6mZXu+cWoBrgLa97Q0tvyaAtbNWOeiX/9MCxYNrSq+JDysNHrMTYdUk8N0D/5usJcPX6yE6GqZNg8GDITS04dovhLjmSJAphGg8qalahRqAESNgzpzGbU+xIE/7krrnuyJTqt2b6e5gzWuTujFrcAB2VgYMeh3vzwjG19mG/CIL81bur3XC92q52JsZ/bf2uUVPGPeGtn1yg7YC/fSm0vNzkiAzHsxFsHYWZMVr+9P2wZIlWpnKm4oXKO3aBQMHSj10IUS1SZAphGg8jz+uLTCxtb2yDGIjUhSlJHdmXFoe51Kr35t5OXcHaz68szdWeh0xqbnMXxOGxdJAQ8/+g8q+7jQeuk+FLpO013+9DRaTViddV7wAKj5MK3N5MTAFyDgH6TFa0vZff4UNG6BLF5gxA4KCGqbtQohrjgSZQojGYTaDtbW2vWhRkwte2rjb0aK4N3N3DXozy9OjlQsvTegCwNaIJN7dUr+r50tcEWTeqgXut74PXp0v2T8evDpp27F7S0tU9rqnOPk7ZVMhjR0L4eFa7+ZFqqoFnZ9+qn0vhRDiMhJkCiEah14Pn3yiDcM++mhjt+YKZXoz0/OISc0t97yj8Rn8ciiBnMpyYwLT+7VmWt9WALy35RRbjjdAiiDXNqVpjTw6gGd7bdvaQVuJblOcI7TXPeAbrG3vXwE5idr24PnQqq+2fWmQCWAwgPMlOUZ/+glWr9amOPTqBVu21P/7EUI0axJkCiEa14ABWsDZBPlX0pupqio7TiXx+9ELnLyQRfi59Crv9+KtXehenAz+sTXhJbXU642ilA6N955V9ph7EDy4E2ZvgTaDS4PM/HTts19f7ZzWxb2hlweZl+veHe64Q9s+dAhGj9ZqzUdE1Mc7EUJcAyTIFEJcXTt2aDXJm0FKHEVRStISxafns+HIec4m51BosvDbkfPsi0orOTeuGqvQbYx6PryzN272VmTlm3jgi30l1YGiU3J4/bcTHInLqFujR78ID+2GAQ9eeczZD/z6aNsXg8yLuk/VPl8cck85pa04r0ibNlpqo7//hn79tH0//aTVnX/0Ua1kpRDiuiZBphDi6snKgnvv1coZLlzY2K2pltZudrRyswPgxPksvg+L48OtZ4g4nwWAp6M2r/R8Zn61Kvu0dLFlyYxe6HUKpxKz+b814Xyz7xzj/reDj7ad4b4Ve8kvqsMcR4OVNt+yqkVUXp1Bb6Vt64zQdYq27de3dFFQzK6qnzdokDblYdUqaNUKTCZYvBgWLKj9exBCXBMkyBRCXD2PPaalwDEYYObMxm5NtSiKwq09fBnewbNk6NxS3As7INCdyb1aAlo5yfMZ1auOMzDInWdv1hbe/H7sAk9+e4icQi2wTMwqYO2+c/X9Nq5ksAZvbTES7W4Au+JqSFZ24FtcZjJiQ/XupdNpi4BOnIBXXwV3d3j++bLnNIOeayFE/ZIgUwhxdXz3nVabHODll8vWy27irAw6erV2ZVq/1tw3JIDhHTy5tacvA4PcsbMy4O6g9QjGppW/OKg8swa1KQlQAUZ08OSGzloS+A+3nqn/euflGfkctB2tDbFf6mKv5uG1WqWgv96FFbdAamTl97Oz08qDnjsHgYGl+3/7DUaNgrCwiq8VQlxzDI3dACHEdSA+vjTR+rBh8OSTjdueOnC2NdKrtWuZfX6utqRkF9aoOpCiKLw2qRuBHvb4utgyKbglZ5Jy2HT8AvEZ+Xx3IJZp/VrXd/PLajtK+7hcr7tg238gLxW+uVurhw6w9l64f5M2JH+5uAOQfBK63aHlPb3IZNKGzo8fh969YdYsrbfT17dB3pIQoumQnkwhRMOyWLTAIjVVq03+xRdNdjV5bfm5anM2z2fmU2Sufg+kjVHPwyPbMbmXH4qi0NbLgZu7tQDgg61nMNXgXvXKyh76z9O2LwaYAAnh8OerV55/4Av4dDSsnwu7Fpc9ptfDf/4D7dtrQ+bLl2vbr7wCudXv+RVCND8SZAohGtbixbCpuJThBx+Av3/jtqcBtHTReu7MFpWE9OrNy6zIwyPbAhCTmst3YXF1blut9ZsDRntt284DeszQtv9+D6KLFwTlpcMvC+DHR0AtXqy0/S3IuWRluaLA+PFw+DC8+y64ukJOjjZns0MHWLlS+0NECHHNkSBTCNGwnJy0uXrTpzebxT41ZW99ybzM9Lr1znX0ceKmrj4AvLHhRMPWOq+MnRvc8DJ4doQ7Pofx74J3V0CFnx+DyG3wfl/Y+6l2fuuBYOUIBRmw7fUr72dlBfPnw+nT2gIwgwFiY7VsA1FRV+1tCSGuHgkyhRAN6957tQUfH3zQ2C1pUBd7M2syL7Mi/xrXCRujjuTsQt7YeKJkf1JWQcPVPS9P39nwj1BoM0RbjT7+f4ACSSfgi1u1SkEGG20B0d0/wpDHtOv2fab1cpbHzQ3eeQeOHtWSt//jH2UXCUmJSiGuGRJkCiEaXvv24OLS2K1oUCXzMjNqNi+zPK3c7Jg/SisJ+dWeGH46GM+DK/fTd9FmFq49WOe21ppfH+h7f+lr1zYw7y8YtlBbDNR3NqCAxQQxuyu/V/v28MMP8NZbpfvMZhg4UFsYllHHpPRCiEYnQaYQov6tWwcPPggFBY3dkqvGz7V0XmZMai77o9P4cnc0m45dKFOOsrpmDw2gvbcDqgqPrA7jtyPnAfjtyPk6B7F1Mup5aDsGOt4C928Gj3alx2xdoEV3bTv6r+rd79JFYF9+CXv3wptvQtu28OGH2up0IUSzJEGmEKJ+nTypDZF/9BE8/nhjt+aqsbc24Gavzcv8MTye7SeTSM4q4EhcBrsjU2t8P6Nex6sTu5W8vnjvvCIzh2IbsZfPxhnu/BamrQIHzyuP+w/RPkf9pa0mr8minttug+eeAxsbSE6Ghx6CHj1gQzWTwgshmhQJMoUQ9ScnB6ZM0cpHenlpibmvIxd7MwH0OqWk5GTo2RTOJufU+H79Atz4cGYvnr+lM9ufHFEy73PP2ZoHrVdNm8Ha54SDsOp2eCMAUs6UHrdYKg48HRy0RP0nT5YuEjt2DG66Sfs4erRh2y6EqFcSZAoh6oeqakPkR45oZQa//hpatqz6umtI3wA3grwc6Bfgxv1DApjWtxXeTjaoKvx2JIH03MIa3/Ombi24b0gADtYG+gdopR9Dz6ZUcVUjaj0QUEC1wOlNkJ8Oh7/VjuUkw+Jg+KA/5GdWfI9WrbTURqGhMLg4aN2wAZYsaejWCyHqkQSZQoj68ckn2pw60Cq6jBjRuO1pBE42Rm7t4cvgth7YWxsw6HXc3L0FtlZ6Coos/HQooU7zKfsVB5n7otIaL1F7VezcilMdXeL8Ie3zlpchLUqrDLRz8RWXXqFfP9ixA775Bnr2hJdeKj1msVxXc36FaI4kyBRC1N2+ffDoo9r2+PHw1FON254mxNnWyE1dfVAUSM4qYMvx2i0EAugf6A5AdoGJYwmV9AQ2trYjy76O2w/xYVploIt2vQ9ZF7Ttojz4/FZ4Lxi+/0fZ4XVFgdtvhwMHwPOSOaArV0KnTvDtt1ovuhCiyZEgUwhRN2azNn+usBACAuDzz7XhclHC392ewW09ADiekEX4ufRa3aeNux1exfM8m/S8zGFPwE1vwh3FQWVWAvz4KKCCWyDYuEBRLqycAvHhcGgNnN0GqZEQvhI+GQFHvy8bPCpK6bbJpC0QOntWC0CHDdNWpQshmhT5P4EQom70em2YvH17LXWRq2tjt6hJ6uPvSlsvBwC2n0wmLr3mSdsVRSnpzazNivWrxtoR+j8A7caCzqDtuzhkPvJZGFM87H3hMCwbA9veKL3Wzl2rGrT2Hvh8POSW8z4NBvj1Vxg7Vnv911/a0Ppdd8G5cw33voQQNSJBphCi7vr101YBBwc3dkuaLEVRuKGLN272VlhUlV8OxZNdUPMckBfnZe6NSr261X9qw2hTdn6mnQd0HA+9Z8E9P4NzazAXQmZxjfb7N8ED26DNUO111A7449Xy792li7YY6LffoHNnbd/KlVo99Oefh+zsBntbQojqkSBTCFE7W7dqQ+QXXZpUW5TL2qDnlu4tsDLoyCkw8/PB+JIFPLmFJg7EpFVZq3xAcZCZkVfEifNZDd7mOmvZu3S75wytMhBAwFC4+3uwLe759g0Gv77g0gru+QlCntb271sGS/rDp6Nhw9OQl1b2/jfeCAcPamVLPTwgLw/efluCTCGaAAkyhRA1t2kTjB6tDVemNuFh2ybI3cGasV28AUjIyOePE4mcTc5h5e5otkUkse5ALOZKeijbejng4aDNy/zpUPxVaXOdXBpk9rqn7DH3ILjzO+h6G9zybum8S0WBwfPBsYX2OukExO6F3R/A9v9e+QyDQUufdfo0PPEEvPAC+PiUHpd/o0I0CgkyhRA1c+yYttjCbIbEROnBrIW2Xo70D9R6JI/GZ/J9WBw5BWbgYg9lxSvHFUXhjj5+AHy9J4b8InPDN7guukyCzhNg5HPg0fbK4y17wW3LwLdn2f1GWxi7CBQdeHaCdjdo+8NXQVF++c9ydoY33tACzYvy8rRpHJMmwalT9fKWhBDVI0GmEKL6kpLgllsgI0NLJ/Pzz9r/2EWNDQx0J9DTvuS1j7MNrdzsANh7tvL5lncO8EevU0jLLeLH8Cbem2llp60yH7aw5td2nQJPRcNDu+CWd7SAMy8Njv9U8TXn9sBX07SKQ6DN04yJge+/1+ZuPv44pKVVfL0Qot5cV0Gmoig9FEX5TVGUPxVFWdnY7RGiWcnPh4kTtbQx1tbwww9ayiJRK4qiMLaLDz1buTC0nQd39GnFsHZamqO03CJOJlY839LXxbZkyH3Fzqha591sFmyctOFzZ7/S3sw/F8Hvz8GqOyDhUOm5hbnaavWTv8HmF7V9998PX3yhVZ8ymeDdd6FtW3jvPSiqfP6rEKJuahxkKorST1GU/1MUZaSiKPpL9rsoiuJQv82rP4qiWAH/BaapqjpCVdU7G7tNQjQbqgr33Qc7d2qvly+HgQMbt03XABujnhEdvejTxg29TsHLyaakd3Pv2dRKg8dZg7QA/1hCJnujrpOeuSGPg84IaWdh53twaqMWVEZu1Y7vfK/03Iv7dDottdHJk1pddDs7bY7m/PnQtauWCkkI0SBq05N5F9ABCAa+UBTlOUVRnAAF+LQ+G1fPBgBZaG3eqijKrY3dICGajddeg9Wrte2XXoLp0xu3PdewiymKkrMLOZVY8Qrpvm1c6dzCCYBFvx7nQmYF8xSvJa0HwMxvwMoRDLbg4A2mfNjxllY1aNcHpecqOjBdUnbSzk5L4H7qFMyapfWOnjwJe/Zc9bchxPWiNkHmU0AEMBbwAUYBx4GXanm/EoqiTFMUZYeiKJmKolyRQE5RFL2iKG8qipKkKEqWoijrFEXxqObtWwK9gXuAicC/FUWRrNFCVMdNN2npYe68U/sftWgwLZxt8XfX5mZuP5lEoan8GuWKovBgSBAAB8+lM/bd7Ww8ev6qtbPRBI2Exw/DwghtMRFo8y+Pfq8lcb/IYtJqpZ/dUfZ6X1+tJ37fPrjjjrKLhAoK4Px18DUU4iqpMihUFGW0oii7FUVxAVBVNVdV1bdVVb0BGAc8DNwJhAJ1LbWQBnwAPFbB8X8CE4D+gF/xvi+L2+mgKMq+cj5eKD4vFdilqmq6qqrpQDjQro7tFeL60KsX7N4Nn35atryfaBDD23uiUxSy8k1lykdevhhofA9f/jetJ442BtJzi3hw5X6ON+Wa5vXF1hVsnKFFD+11fgZs+4+23f5GrZcTtProX9wKpsIr79GrF6xZA/ali694/31tvuaiRdqqdCFEnRiqcc584I/iwOxyOrRh6B9VVf1TUZSv69IYVVU3AiiKElLBKQ8AL6uqGll83pPAaUVR2qiqGgX0qeT2u4FFiqIY0Yb2uwBna9I+RVHcAXeAbt26kS3JfsVVkJubW+bz1aL76y8sfftqi3wAvL21hRKyWKLBWQNdvG3YH5PBrlPn8bSFM8k5HEvIpquvI0OC3ErOHdXWmU6ze3P3F+Gczyxg1a5I/nlDOamCrkV2ftjrrVDMhdo8TSCvyzSsMuLRXyheEKRayD13CItnx7LXWkyg6Ev/aDKZsFu8GF1ODjz7LJaPPqLwpZcw3X47ucUB59X+GRSiuatOkBkMPF/eAVVV8xRFGQV0BBaqqtpgCdsURXEGWgP7L3n+GUVRMoHuQFRl16uqmqEoyn+BPwEjsERV1aQaNuMR4AWA2NhYtmzZUsPLhai9Xbt2XbVnee/bR79//5vkrl3Z8/TTmG1srtqzhcZsgfMXFPJMCu9Elv49HHVWJe+siu6yDuWeTgobMvV8fyCWXko0+uskd8gw65a45hYHmEZXNp1VGZaZicsl5xzd9h3xrv1LXntkHWPgmTeJch/B4VZ3l+w3vvIKHdesoc1vv6GLjcXm/vtJ/c9/OHHffdCx41X9GRTiWlCdINMdSK7k+CrgtfppTqWcij9nXLY//ZJjlVJV9WugLr2ti4GvAPz8/CJGjRpVh1sJUT25ubns2rWLgQMHYmdn1+DP02/fjs2bb6KYzXjk5xPSv782H1Ncde1TcvnlSCIABr2CyawNl/fo1QIvR+uy56blsWHJHrJNClb+PQhpf318z6zNG+GQFmTqe9/FyCE3YHSLgz+eLTkn2HyA7i5eFPW4CwzW2H34BDrVTGDyZrzueBvs3EtvOGkSeRERWD/7LIYNG3A7eZJh//wnsUOGYP3hh9i2aXOV36EQzVd1gszzQBAVz7c8ArSqtxZV7GLSuMszP7sAV2USkqqqKUAKQJ8+fXBwaLIZm8Q1yM7OruH/zYWGwtSp2gKIgAB0W7bg0LJlwz5TVKirgwNmnRVZ+SZ6tnbhq9BocgrMpBXqCLzs30InBwf6tXFjT1Qqvx5P4ZZebRqn0Vdb6z5waBUAVv3uw8rBAQbNBSsDHFoD8WEYzu3CcG4X1mGfgXdXyLlQcrnD2Q3Qf27Ze/buDb/9Bps3w4IFcOgQHocPY3Jzk9/7QtRAdQZUNgMPVXLcGm2OY4MqnhMaA/S6uE9RlEC0XsxDFVwmhKiuQ4e0VeTZ2doK3M2btQTWolH1aOXCkHYeOFgb8HPVerLj0stflDK5l/b92nwskfTccha7XIu6ToGgUTDiGXArLg5gsIIBD0Kn8WXPTY+GiF/K7tv1PqScKbtv7zLY/iaMGgUHDpC/eDFH7r8fnC4ZNDt6VEvuLoSoUHWCzDeBWxVFebmC4yOAMxUcq5HiFEU2gFXxa5vij4tB7CfAU4qiBBTn5vwPsLF40Y8QorYOHoQxY7Ryex4esGkTBAY2dqvEZVq6aKum49Lyyk3UPq57C6wNOgrNFoa+8SfD3viT5X/XaH1j82PjDHd9B8OfvPLY5Yt9WvSElr3BxgV63qkt/EmPgS8ngrk4YMw6D7/8H/zxKhxZB3o9plmziBs2rPQ+aWkwbJhWE/333xvojQnR/FUZZKqqehKYCTypKMpBRVFmK4rSU1GUQEVR7kObj7mintpzF5AHbAT0xdt5gH/x8deBn4C9QFzxOVK5R4i6MJlg8mRITNTqkG/cqNV4Fk1OS1ctyMwvMpOcfWVPpZONkQk9fQHIyjcRk5rLSz8d4/frIX9mebw6lW5P+gTmboM5f8A/o2HiErjzW+1YegxE/qltJ0WUXnP8x/Lv+8svWtWgI0dg7Fi4+WY4frxh3oMQzVi11h+qqroOGIw2L/ITtBXep9Aq/GxDy21ZZ6qqrlBVVSnnI6r4uFlV1YWqqnqoquqoqupkVVUrW5QkhKiKwaBV8wkMhD//1PIHiibJ3d4KWyutmm9FQ+YvT+jK5/f143/TetKjlQsAT647xNnkHLZGJLL9ZE2TajRjboHaMPqgR6DbbVceDxoJfv207VW3wbEfIeVU6fEzW7WezdyUstfdeSfs2gUDBmivf/0VunWDhx+GZPlfkhAXKZXVxi33AkVpizYv0hY4rKrqgYZoWFPXp08fdd++fY3dDHEdyM7OZsuWLYwaNap+Fx2oatnE6kVFYDTW3/1Fg/jpYDynE7Np5+3ALd19Kz33XGou497bQVZ+2bmDX97fj6HtPBuymc3HnqXw60JtW9FBx1uu6MFUrZ3Z0G4RQ8dOxMHR8ZIDqpbQ/amnICZG2+fsDG+8AQ88cJXegBCNS1GU/aqqlpunvMaZ1FRVPa2q6jeqqn5+vQaYQjR7K1bAxIlQeMmQqwSYzcLFIfOK5mVeqpWbHf+Z0v2K/ev2xzZI25ql7lNLezNVS7lD5EpBBi3Td2P32TD4aCgUFdeJVxSYNg1OnIDXXgMHB8jIkIIFQhS7TtL1CiFK/O9/cO+98OOP8HJF6/lEU+VXHGTmFppJzal6Bfm4bi34bFYf3rq9B4+N1irpbjp2gbzCBqud0bzYOMHsTdB6UNn9urJ/dHWPXYkuLRLOH4Izf5Q919YWnn4aTp+GZ54p24uZmQn79yPE9UiCTCGuF6oKL70Ejz2mvR40CBYubNQmiZrzdLDGxqjNyzyVWL3StiM7ejOltx8z+/ujUyCn0MyfEYkN2czmp3X/sq/7VTLcXdGCIG9vePXVsqMCr70GffrAPfdAXFzd2ylEMyJBphDXA5MJHn0UXnxRez1mjJZ6xcWlMVslakFRFDq10OYFHo7NwGyp/rx6T0drBgVplYB+DI9vkPY1W60Hlm4rOi3PpltQ+eee+BXy0qq+p8kEGzZo2198Ae3aaT+DOTl1bq4QzYEEmUJc61JTtSTr77+vvZ48GX76CeztG7ddotZ6+LkAkF1g4nQ1ezMvGt+jBQB/RCSSlS9zB0u0Hgj2XmDtBDO+AZdWMHc7TF1VcorZvT0YbKAgA1ZPh09GwK4lFd/TYIC9e+G998DNDfLytNGE9u3h88/BYrkKb0yIxiNBphDXMlWFCRO06j0A8+drq2GtrSu/TjRprvZWBHhofySEn6tGj9olxnbxwahXKDRZ2Hz8QtUXXC9snODRMPi/Y9BujLbP2gEChqHauGBWjBSM/whGPa8di9kF8Qdg47/AUsn8VqMRHnlEm6/5+OPa6/h4mDUL+vaF6OgGf2tCNBYJMoW4limKlk7F0RE++wzefVfrXRHNXs/iHJjx6flcyMyv9nUudlYMCHQHYOfplCrOvs5YO4C1Y9l9Nk7k3rWBPzsuwuLRATrdeuV139wNaVGw6QVt9Xl6zJXnuLrC229r5SgnTdL2ZWSAj0+9vw0hmgr5v40Q1xpVhdzc0uHwgQMhKkobrhPXDH93O1ztjKTlFvHXqWRu7t6iZEFQVfoHuLHjVDJ7olIbuJXXBtW5FTk2xcGgU0sw2ILpkmT4J36G+HDILE4N9cM/YPJSsLK/Mmht1w6++w62btWGyy8dVdi4Efr3l7nS4pohPZlCXEtyc2H6dC0HpumSBNwSYF5zFEWht7/2fY1JzWVVaEyFVYAu1y9A68mMTsnlfEb1e0EFoNOBezkLgjIvyT16dju81QH+7QdvtoWj3195fkgIjBxZ+vrcOa2Hs21bWLKk7M+vEM2UBJlCXCtiYmDIEG3O5ebNWqlIcU3r2tKJYe090esUMvOK+G5/LGnVyJ3Z3c8ZK4P26z/0rAyZ15iDV/XPzUmCtffA6S2Vn3exglxKilaesnt3rVxlDavyCdGUSJApxLVgwwZtEUFYmDYP8z//0eori2ua1pvpyrR+rbCz0mOyqByNz6zyOhujvmRO556zMmReY8ol/+u8+0d4/GjV16ycDEfWVXx80iQ4eRLuukt7ffw43HwzjB0Lhw/Xrb1CNBIJMoVoznJz4R//0FIUJSaCk5OWnujJJ8vWJRfXNC9HG7r4OgNw4nwmlmrkzuwfoA21S5BZC22Glm4HDANnP5j8KQSNrDyJ+7f3QeRWiDtQfvoiPz8tn+bevdqoBMCmTdCzZ2mOWyGaEQkyhWiujh2D4GD44APtdb9+2pDbzTc3brtEo7iYoD0r30RsWtVzM/sVB5mnErNJyS5o0LZdcwY8CAP+AVNXlv4x1/12uGs9DHoUjMWL7uzLGVb/YgIsHQG7K8mv2acPbN8O334LgYFaQNqpU/2/DyEamASZQjRX7u6QlgZ6vZbg+e+/tZWr4rrk7mCNt5MNAMcSqh4y7+3vikGnBUh7o2qWa/O6Z7CGG1+DTuOvPObSCubtgPmH4B+hFd/j92e1+ZYpZ2DjM7DqDii6ZBGWosCUKdofk59/DnfcUXosIUGbey3zNUUTJ0GmEM3JpStOvb1h1SrYuROef17yX4qS3szTiVkUmCpJEA7YWRno2lIbYt91JrnB23ZdcQ8CV3+wda38vI3/gsW9YNf7cGojHPjiysTu1tZw991lp788+yxMmwaDB0NoJYGsEI1MgkwhmgNVxfDZZ9Cli1Ym8qIxY7RhciGAjj5O6BSFIrPKqt0xrNkbQ8T5rArPH9ZOq2P+Z0QSqvSK1b9LA8P2N8LY18oe3/1B2de/PaGlPTrzZ8X3NJm0+dcAu3bBgAEwc6aWXUKIJkaCTCGaOCUhgf6LFmEzf762+vRf/2rsJokmytZKT1svBwAy8oqIT8/ntyMJnEvNLff8kI7anMGY1Fwik3OuWjuvK/f8BMF3wvj/Qf8HYea3cNMbFZ9flKslc6+IwaAt7tuwQfujE+Crr6BDB62HM6viPyqEuNokyBSiqcrPh9dewy44GJ+LOfQmTYJXX23cdokmbVQnL27o4s2w9h64O1ihqrDx6HnyCs0kZuWTmFU676+HnwuudkYA/jyR2FhNvrYFDIMJS8DRR0vk3m4M9J8L/eZWfE1mXPn7VRW2vQFhK7XURuHh8NFH4Omp/b5YtEibl33mTIO8FSFqSoJMIZoaVYV167TVpM88g5KTQ6GDA/kffqjt9/Bo7BaKJszGqKeLrzO9/d0Y160FRr1CVr6JT3dEsmp3DF+FxhBfXBlIr1MY3t4TgK0RSY3Z7OvPDa/C4McqPl54Wc9yQRb89Cj8uUjr6cw6r/Vqzp0Lp07BU0+BlRW0agUBAQ3adCGqS4JMIZqa77+H227T6o3r9RTOncvmDz7AdOedkvtS1IiHgzXD22tD4qbi3JmqSpmE7SOKh8xDz6aQXSClDK8agxWMfrHi41tfh6I8+ONVOPg1fHajtjDoovNHSredneH11+HECfjsM63H9KJPPoGIiHpvvhDVIctRhWhqbr1VKynn7Q3vvEOhvz9FW6ooSSdEBbq2dEJRIL/ITFaBifCYdE4lZjGyoxd6ncKwdp7oFCgyq3wfFsf0fq3R6+SPmavi0j8avbpAyD9h/3I484e24rwoD/YuLf/a356AoP1lA8rLezCPHIEHH9TOeegheOEFcHOr//chRAWkJ1OIxlRYCG+/Xbaah14Pf/4JGzeWTuwXopYURaFrS2f6tHGjt7+WUqegyEJ0ijYc62pvRXBrbf+z3x9hwL+3sPO0pDS6aqau0laeT18NnW+F21eAzgiqpeIAEyA1El52hQNfVnxOYiK0bKmtSH/vPWjbFt59V/u9I8RVIEFmbZ0+rZX9upgEu6iosVskmhOLBdavh27dYMECbcL+pUNabm4yNC7qnZONkZYutgCcvFC6CnnhDR0I9NCq1CRlFfDYmnAy8+V32lXR6RaYsUbLqwlg4wxtBlf/+h8fhtxUiN0HcfvBfMn3beRIbQj9lVfA3l4r3vD449C1K/zwgyRzFw1OgszaysrSgssXX9SCTXd3bZhz8WLth1p+eEV58vO1OVKdOsHkyVpKIkWBe+4B1yoSNwtRD9r7aAnbzyTlUGTW6mcPDHLnj4Uh/PjwYKwNOhKzCnhzg8zjazRdp5RuuwVCl8mVn7+4N3w6CpaO1OZwXsrOTkttdOoU3H+/9vvm1CmYOBEefrjemy7EpSTIrK02bbRVfUFB2uusLC132aOPagHEu+82ZutEU1NUpKUe8vfX/t2cPKntHzUK9u6FTz8Fr3LqHAtRz9p5OaAoUGiyEHVZbszufi48OkorTboyNJoDMVJuslH0mA7ubbXtQY+U9nJe5D8YfINLX+ddUqDh73fBVM5weIsW2u+ZAwdgxAht35QpV54nRD2SILO2XF21/GSnT0NkJHz8sbYi+GJv1JAhpedmZUHv3tqw6IYNkCNJj687BoM2PJWYqM25nD4d9u+HzZu1fxtCXCX21gZ8i4fMY8pJ0j5naCDtvBxQVVi6PfJqN08A6I1wz88wdSX0mgWDHgXXNqXHdQZ4YKu2WKg8y8ZoVYPyyvkjoWdP2LJFqxY0cmTp/ogIbVg9t/zE/ULUhgSZ9SEgAB54ANauhaQkrWeqV6/S49u2aX89vv023HSTNt9uxAh47TXYswfMldcYFs3Q3r1le7MVBZ55BubP1/4w+eqrsv9GhLiK/Fy1IDOuOF/mpawMOu4boq1S/jMikRxJa9Q4nFpAp/HaynA7N5h/sPRYix7a54kfgGfHK69NCIcvJ8K62eXfW1G0cpSXeuIJeP55rXLQypXavHEh6kiCzPqm10OfPtrniwIDYeFC7S9I0Fb2bd2qBR39+2vVGhISGqO1oj6ZTNqUieHDtXriCxZovdwXTZyoBZ5t2jRSA4XQ+LnYAZCSXUhu4ZVB5NguPuh1CvlFFv6QSkBNx5Rl0OP/27vv8LiqO43j3zOj3rubXHDFDbANmGoMJllggVAMJqGFFCAkJJDQwhJwKAuEJCwlZIHF4FBCIIYQSmKqjROKscG927g3SVbvmrn7x5lRHcmyNdadkd7P88wj6d4zo2PwWK9O+Z3vwJSb7df9j4Iffw6374Thp7dtv+H9zu0PaGiwS79iYmD7drj8chtC//WvsHZfeh+FzO4wZgw89BB89RXs2WNHsb73PXsyA0BKCvTt29T+7bftyOgrr0BRkTt9ls5paLBTT9dcY/8fnnsufPyxvTdqFOze7W7/RELom56AJ1C9YEdx29HMrOQ4ThiWDcDby/QLcMQYPx3O/yMkpLW8HpcMl82BmaUw6aqW916/1n5sqIO9q0OHzpgYePhhWLnS/hsGdjbm5JPh4otb/rIscgAUMrtbXp5dj/fMM7Bli10H8+yzLcvVzJkDTz8NM2bYUc6jj4bbbrNhpqam/deW7lVVZX9ROP10u2M8+AvBlCl2RHPFCjjhBHf7KBJCXIyHPmnxAGwPMWUO8J/j+wF2ylwnAUWRgce2/HrZy/DSJTD/AXjiOFjw2/afO3KkXTv+wQdwZGBK/tVX7WbW9esPXZ+lx1LIdJMx9k09bVrL60cfbYNKTIz9rXPxYnjwQRtmMjPh2mvd6W9v5vPZJQ7BUUqwpUFGjbKfT5gA999v11vOnw9nn93yJA6RCDMguC4zxEgmNE2Z1zb4+WD1nu7smnRF/jFtr637Byz4nf38w3uhfj+DFaedZn/uPPOMnaE57TQYMSL8fZUeTz8FI9F119mgUlxsp85vuKHp5JeaGoiPb9n+f/7H/mOwdWt397Rn8/thwQJbS27AALtZ6667Wrb53e9sOaIvv7SjzcGSViIRLliUvbCilpr6tpsPM5PjOHlEDgAvfqZ/W6JGzgj41h/glNtCB06wZY4ASrba6fQtn7Rt4/XaZV3r19vSR83dequtlCKyHzq7PJKlpMBZZ9kH2M1B779vT2sI8vng7rttIAX72+bJJ9sRtpEj7dfDhkFCQvf3PxqtXg1//zt8/jl88oldQ9tcSQnU1jYFfZUfkijVPyMRY+xkyc6SaobmprRpc+XxQ5i3toCFm/exYkcp4waku9BTOWATLrMfT7oR7uvT9v68+2H3cljzlv166Z/tek4InBhkwBuIBykp9hH08cfwm9/YxxlnwG9/q+NvpV0ayYwm/frZXX8TmhXhLS626/6S7ZFwrF8Ps2bZ3zTPP98G0ueea/k6f/wjvPOOndpt6KVrrXw+u2byxRdbXl+wwI5Ivv56U8AcN84G+dWr7eat1iPJIlEoIdZLbmpgXWY7U+anjMxtPG5y1r+/7ra+SZjEJsCIb4a+FwyYQY4DtRXw6ET44/GhC7oDpKU1rTX/5z/t2s3rrrPl+0Ra0UhmtMvJgbfesmWRPv8c3nsPliyxYXPjRnvSzMiRTe2Li+0/CEExMbbE0ogRTSOfZ51lT6bpSXbvtv99go8vvrBF8sGuN+pnNzlw3HG27unkyfbxzW/a6gAiPdCgrCT2ltWyuaiSKeS2ue/xGL574hDufGMlby7dydj+6Rw/NJsx/dNCvJpEpPjUps+/eR/MfxBqy9q2qyyETfOgNLA0Yu/KlqcKBR11lC1t9Ne/wi23wObNduDixRft8ZU//al+EZdGCpk9RVycnSY/+eSmaw0Ndp1m8/JIe/faELVli11z2NBg1xSuW2fXf4KdLg6GzJoauPRSG0CDIXTkSLvrvfmOeDf5fPa36N277WPy5KaTlxYvtmeEt7deNS3NlucIhswjjlC5Duk1hmQns2hzMUUVdZRW15OeGNumzYUT83lo7lrKaxq4561VGAPv3jCFEX1SQ7yiRJypv4QdX8LEy+GEn0De4bD1M1j/LuxqVuD9t8NbPq9qH+0yBi66CM45Bx55BO67D8rKbOhcsQJmzz40fxaJOgqZPVlwlLK5UaNsiKqttR/Xr28KmcHPm+8i3LgRXnut7WunpTWFzgcegEGD7PWiIntaRFyc/W02+LH55xMmQEaGbd/QYIv/Nr8fFwexsfb4zd277aabRLtJgR074Fe/std37bIf9+5teTrFRx/B1Kn284EDmwKm12tDZHCUcvJk+99Du8Cll+qfkUhcjKfxHPMjB2a0aZMcH8Oj357Anz7ZzCcbi6ht8POvDYUKmdEiZwT8bEnT18NPt49Bx8ELHZxdXr7Lrs/0tv3Fo1FCgl2addVVdlPkM8/Az38etq5L9FPI7K3i423ts9GjO24XGwtXXtkUQgsL7fWyMli0yD5uvbUpZO7YYXfDd2T+fFuiCWDnTjuy2pEFC5rOgnccW1e0Pca0LGCflwdPPGHXVU6aZMsOiQgAXo9hcHYS6/dUsLkodMgEOHVUHqeOyuOa5xcxd+UevtpawlUndm9fJcyGnmpPD1r6Uuj7b/zYPs54ECZfY/9t9dXD6jdh8AmQ2myGLC/PTpnfeWfTrBDAZ5/BH/5gy7vl5x/aP49EJIVM6djIkS03DhUXN414Bj8Gp6bBjhaOGmVHSuvq7Mfg53WBheTN1+vU1u6/D81PzcnLg+nT7T9kffu2fPTrZ6fxY1r9tf7Rjw74jy3SWxyWk8z6PRVs21dFvc9PrLf9kf0JgzJtyNxW3I09lEPC47WnB21fCEUb7LUzfwP/uKVlu3/eas9CP+dRW2tz/gMwbBpcHmKGq3nAdBw7qvnpp/aAkZtustPpKW2rGEjPpZApByYz057Lfeyxoe+PHQtr1oS+5zh2I1Lzc90HD7btW4fS2lpITbXhsflvwHFx9gQKEQmLIdl293i9z2FzYWWH0+ATAiOd2/ZVU1hRS06KNnhEvcEnNoXMyddAdQnM+++WbZb+2T6CNn6w/9f1+201lOAM2D332Hqb991nZ8e0TKlX0P9l6T7G2JDYPGTGxdmRzyOOsCcdnXii3e195pl2inz4cNX4FDmEkuNj6JNm32NvLdvF859tYU9Z6BNhxuen4/XYDX9LtpZQWduA3x/iLGyJHtPugnEXwkXP2a9PaTWSOfrc0M+bmQ5zfgi+dsrgeb12Fmn9ejuKGRtr19F/73v23/p588L1J5AIppApItLLHXtYFgmx9pe/wvJa3l62i3qfv027pLgYRgVGOv8wbwMT7nmP6178slv7KmGWnA3TZ8HY8+3XrauGXNTBTvHlr8Cqv7W85m91elRGBjz0kK0zfGFgo9FXX9kT1HQeeo+nkCki0ssNz0vhmilDOfeo/niMobS6noVfhy5hM2FQBgBfbS2hrsHPP1fuZsPeim7srRxy/Sfaj4mZdlr7m/e233bO9+G1q+Gj++HevvCbw6BgXdt2w4bZ2prz59tNmJdfrvPQewGFTBERweMxDMtNYeLgDAAWbS6msKLtxrwJgzLbXHt5oc4271Eung2Tr4XvBmonH/djmP4s3LgSzn28bftlf7EbghqqoaYUPg3RJmjKFFi40Fb9CHIcW8/40Uftun3pMRQyRUSk0eTDsklLjMXvOCza3HY08/hh2cTFeIiP8fAfY+252HO+3E5tg69NW4lSGYPgzAehT+BMco8Hxl0A6fm2qPv35nb8/L2roWwXLHsF6qttiFz1d3teevD1mu8y/9vf7FG+P/uZLTf35pv2ORL1tLtcREQaxcV4mDgog3lrC9hYUEmDz09Ms7JGAzISeeenJxEf47Wn/6zaQ3FVPe+u3MM5R/Z3sefSbQYdBzNL7edv3QiLZrW8v30hPHIk+AIj4YNPgi3/sp//qrBtgfdjj7U7zmfPtmXxzj0Xpk2D3/3Ono0uUUsjmSIi0sLwvBSMwZ4EVFQV4n4qA7OSyM9MYsoIe+b535fu7O5uSiQ48Wehr/uaLbUIBkyAnUvath0wwNZjXrSo6aCODz6wp8P98IctayVLVFHIFBGRFlITYumfYY9yXb+nvMO2Z4yzJ798sqGQuoa2O9Klh8scAt9/D0b8B0y7E5KyO24/+xz45LHQ9yZNsqWNXnvNbhRyHFtb89Zbw91r6SYKmSIi0sbIQKmiTYWVIcsZBZ0y0o5kVtb5Qq7hlF5g4LFw6Stw8i9g5Bkdt22ohnfvgEeOgpJtbe8bA+efDytX2unyvn3h179uuu84Wq8ZRRQyRUSkjeZT5luKKttt1z8jkZF97CaO+esKuqt7EqlOucWObgYdfjbcvBH6tVpbWfw1vHEd7FoK9SGK/8fH22Mpt2yBIc1e75134IQT7HGVEvEUMkVEpI2U+BgGBKbMl+8o7bDt1FF5AMxbq5DZ62UOgZ8ttRuDZpbCjBcgOQcu/1vbtl9/DE9OgTfbWdcJ9lS4IJ/Pnh702Wc2aH772zaESsRSyBQRkZCODJxVvrmwiqIQNTODpgamzNfuKWdnSXV3dE2iRfAEoaQsuKMAZrwIN62HmMSmNstehsWzoWBtx6/l9dpamuPH269fftkeS3z77VBWdmj6L12ikCkiIiENz00hNcFWuluyraTddkcPySIpzh5L+cnGou7omkSjmDgYfTak5MHQqS3vvflT+N+ToWo/63q/8Q17LOWTT0JeHtTWwv3329ODnnrKjnZKxFDIFBGRkDwew1GB0czVu8pYu7uc8pq2J7LExXgYNyAdgBX7mVoXAeCCp9pe89XC70fD+zNhw/uw8Gk7wulraNnO64Wrr7Znn//yl3b95t698JOfaPo8wihkiohIu8YNSCfWa6j3ObyzfBfP/nszpVVtg+a4/jZkrtqpaUvphIQ0+P77ba831MC/HoYXLoR3brIjnC9dBMWbYdUbsPYfTW3T0uC//xvWrIEZM+CnP4WhQ5vu64hK1ylkiohIuxJivUwb3Ye8tHg8xuDzO2wOsdt8bP80AFbuLMXvV4kZ6YSBx8Cos/bfbuOH9gShV66AP18CpTta3h8yxK7P/M1vmq41NNiThK6/Hoq0hMMtCpkiItKh0f3SuHTyYIbkJAGwq7RtyZngdHllnY8t+9qeEiQS0pkPwjfvhTHf6vxzSreD3wc1rUbNPc0izQsvwJIl8PjjMHw4/P73UFcXli5L5ylkiohIp/RNSwBgT1nbkDksN5n4GPsjResypdMyBsEJ18P5T8EPPuzcc549Ex4cAg+Pg51fhW4zY4adSk9JgZIS+MUvYMwYeP11FXPvRgqZIiLSKf3SbdmZfZV11NS33MUb4/VweF97StBKrcuUAxWbAPmT4KLZTdf6HmFHOVtzfFBbBrWlMOcHsPafUFvRsk1iot0UtH69Pf/c44GNG+GCC+DUU+HLLw/tn0cAhUwREemkvLT4xrKHu0trcFqNCI0NTJmv3KmRTDlIY8+DO4ttIfdrFwRGOZ9sv33RBvjzDHjt6tD3+/a1pY2++gqmTbPX5s+HWbPC3nVpSyFTREQ6JSHWS3ayPYHly63FPPnxJj5YvafxftPmn7I2AVSk0zytosmRl0D+MR0/Z+3bULge/P7Q9484At57D958EyZPhpkzm+75fFCldcSHgkKmiIh0Wt/AlPmWoiqq63ws31FKabUtFTNhYCZgp9P/sWK3a32UHug7r+y/zeNHw0f3Qvke2LkEalqNqBsDZ59tj6XMyWm6/qc/wciR9mN7IVUOikKmiIh0WnDzT5DjwMrARp8x/dOYEjhi8oF/rKGuQT+wJUySsuBHn8LpM+GXO+DHX4Rut+B38H+nw1OnwAODoHw/v+z4fPDrX8OOHXDllbbs0YIFYe9+b6WQKSIinZafmYgx4DGGgVm2pNHKnWWNtTF/eebhGANb91Xx9IJNbnZVepo+Y+CkGyE+BXJHws/XhG5XurXp803zOn5Nr9dOo593nv168WKYMgWmT7cbhaRLFDJFRKTTMpPjmD4pn0uOHchph+cBUFHbwNeBAu2j+6Vx0aR8AB6au5Y/L9za7muJdElKn/23eetGuLcvzEyHN34SunzRiBG2tNFHH8GECfbanDm25NHNN0OpNrIdLIVMERE5IPmZSfRJSyArOY4BmXaN5ppd5Y33f3X2GCYMygDgl68tZ+m2Ehd6KT1e6w1CodRXQUO1/fyr52HLJ+23nToVvvgCnn0W+vWzxduffBJq2taFlc5RyBQRkYM2so+tjbllX2XjlHlqQix/+t6xDMiwAXTB+gLX+ic93K1bIHsEjJtuyx5d+EzH7Z87y54W1B6vF777XVi3Du68E+6+G/o0GzHds0fF3A+AQqaIiBy0Idl2XWZtvZ9dzU4CSk2IZdJgu9t8dbNRTpGwSsyA6xfB9EC4HD8dfrbMbg469Y7Qz7k7C77ez+aelBS7IeiGG5quVVXB0UfDGWfAihXh6H2Pp5ApIiIHLSMpjoykWAC2FFa2uDe6n62buWqXTgCSbpQ52G4OOuVmO7p5/lNwym0t28w+267T/L/T4eVLwVcP1SUdj1K+/DJs3w7vvgtHHgnXXgt79x7SP0q0U8gUEZEuGZKTDMDmopYFrUf3Sw1cr6SqrqHb+yUCwJEz4NRfhq61uf0LWPMW3JMDDw6GX2dAUTu7yq+6Cl59FQ47zNbTfPJJGD4cHnxQ6zbboZApIiJdMiTbhsw9ZTVU1jaFyTGBkUzHgTW7NWUuLhv5H3DpnP23e2wiNNS2vW6MLW20ejU89BCkpUF5Odx2G4webXeoSwsKmSIi0iX5mYnEeOyh5tuKm0Yzc1PjG4+hXLVTU+YSAYadCkd+Z//t7s2DbQth3gNtRzbj4+Gmm2DDBrjuOrtZaPNmWL78kHQ5milkiohIl8R6PeSlxQOwt6xpBMgYw5jAeeartS5TIoHHC+f/0a7V/MW6jts+8w2Ydz/87brQ93Nz4Q9/gGXL4IorbPAMqq2FraoRq5ApIiJdlpdqj5vcW95ymjG4+UchUyJOah8bNn9V1HG7bZ/B3tXw5fPwyhVtz0QfMwZmz4akpKZrjzwCo0bBHXfYKfVeSiFTRES6LDc1MJJZXoPTbIducPPP6l3l1Pt0lrlEIG8MTLqq4zZPHAd//wmsegM+vLfjtj4fPP203Qx0330wciQ884y93ssoZIqISJflBUJmbb2fspqmzT9HD84CoLrex7LtJW50TWT//vP3cOMq+MliiE/ruO3Cp+CrF+BP59kRzta8Xli0CG65BeLiYPdu+MEPYNIk+PDDQ9L9SKWQKSIiXZadEo83sPmnoLypnMvArCQGZdlpxH+t38+0pIhbPB5IHwA5w+GWTTZwDjut/fZv/Bg2fWRHOENJT7eljVavhosusteWLoVp0+Dcc2HXrvD/GSKQQqaIiHSZ12PITrE7yZtv/gE4cXg2AP/eWNjt/RI5YN5YGzgv72RJosXPwZp3Qt8bOhReeQUWLLCnBYE9Hz01NSxdjXQKmSIiEha5KcF1ma1DZg4AX20tVlF2iS4XPbf/Nm/+DF7+tj1B6NV21naedBJ8/jk8/7zdFJSS0nTvq6+gri4s3Y00CpkiIhIWeWnBHeYtTz85YZgNmfU+h4Vf78Pnd1i8ZR+1Db1vI4REmbHnw40r4c5iuxN9ZmnH7Ve+1nb3eZDHA5ddBhdf3HRt3z47hT5uHLzxRsfHWkYhhUwREQmLvoGQWVnrY3Ozc8yzkuMaT//588Kt3PiXJVz4x0/5n/fXu9JPkQOSnm8DYtA1H8M37oHMIaHbv3gRVHZy/fHcuVBSAuvXw3nn2cC5ZEnX+htBFDJFRCQs+qTF0z/DBs0FGwrx+5tGZS4/fjAAc1fu4e9LdwLw0ucqVi1RqN+RcOJPQ5+FDrDtc3hoKOxaZr8u3gLLXoGGEFPi3/42LF4MU6farz/6CCZOhO9/v0dsDlLIFBGRsDDGcPKIXAAKy2tZsbNp2vCSYwZy8oicVu27tXsi4ZU7Cu4qgROuD33/yZNhzdvw+DHw2g/hk0dCt5swwZY2ev11GD7cTpnPmgUjRsBjjx2y7neHXhcyjTGPG2M+M8Z8YYy5zO3+iIj0JP0zEhnZx+6c/WhNARv22tNOjDH8ZvoR9AkcPwlQUlVPcWXP3PAgvYQxkDOy/fsvfwd8gY1wn/1vx69z3nmwciU8/DBkZEBlJcTEhLO33a5XhUxjzDhgjOM4xwFTgbvd7ZGISM9z6uG5ZKfE4Xcc3l62m8IK+0O2X3oiH99yKp/fPq2x7abCCre6KRIeR8yAwSfCuAs7bldVaHegP3E87FkVuk1cHNxwA2zYAPfcAz/8YdO90lL45JOwdbs79KqQCewEqowxsUAaUOJud0REep6kuBgunJhPcrwXv+O02AQUH+OlT1oCGUmxAGwsqGzvZUSiQ0w8XPUOTJ8FF80G4+24/d5V8NrVHbfJzrbnnjcfybzvPjjxRJgxA77+uuv97gYRFzKNMZcYYxYYY8qMMW0KqhljvMaYh4wxBcaYcmPMHGNMTqjXCqEYWAusA5YC94Sx6yIiEpAcH0P/jEQAClrVzQQYmpMMwCaFTOlJxp4Hvyrcf6mjPcvh0Ykw+1y7E92/n3JePh/Mn28/f+UVGD0abrsNysrC0u1DJeJCJjYIPgHc0M7924BvAZOB/MC15wGMMSnGmEUhHncF2n0DGB54jALuMcYkHqI/h4hIrxYszl5QESJk5tpi1JsKNF0uPUyw3NE5j8DgkyD38NDt9m2Er+fbneh3Z0FD2/dJI6/XTpX/3/9Bnz5QW2uPrRw+HJ58Ehoi85CDiFtR6jjOXABjzNR2mlwN3O04zqZAu1uADcaYIY7jbAaO7uDlDbDPcRyfMaYciOUA/hsYY7KBbIDx48dTUaF/HOXQq6qqavFRJFokeX3U1dWxp7iO4tIyYr1N4xr5aXa6fMPe8oj/t1TvQTkoo6bbB4Dfh3f753jXvUXcktmh29+bR9V33sDfv4MYM2MGnHUWcQ8/TOxjj2EKCuDaa/E9+ig1r7yCM3hw+P8cXWCcCK0uHwiZ7zuOE9PsWjp2HeUEx3GWNLteClzuOM7f9/OaHuAZYAQQD7zkOM7DB9CnmcBdAJmZmTz77LOdfaqISK9T3QDzdtlgeUIfP+lxTfeWFhlmrfPiNQ4PTfbhVTkj6UWmrbqZlNo97d5fPPgadmRMBmNwMGDaTjwnFhQw+vnnGfjxx1T078+Hjz6K48Ju9PPOO2+x4zghk3HEjWTuR1rgY+vFDiXN7rXLcRw/0M7Bop3yGPASQH5+/tpp06btp7lI11VVVfHpp59y/PHHk5SU5HZ3RDrNcRz2fLqNmno/h4/MZky/1MZ7QwormbVuET7HMHj8ZEb1Sengldyl96CEW2zGj+CjmfjyxuLdu7LN/UlbnmTSlicBcIyHqmu/xEnObftCl1xC1cKFeOrqOO2kkxove995B98xx0BuiOd0o2gLmeWBj+mtrmcAh3z1q+M4RUARwNFHH01KSuT+oyg9T1JSkv7OSdTpn5XK9uJqKv3eFn9/xyUl0y89gV2lNby7rphJw/q62MvO0XtQwuaUG+HY7+JNzIRZZ8LW9ksTGcdP8rs/h4R0mHYnZA1t2eC001p+vXUrXHklxMfbHerXX28/d0Ekbvxpl+M4JcBWYGLwmjFmKHYUc5lL3RIRkXbkpgY2/7TaYe71GC6YOACA177cQYPP3+19E3FVYqb9eMUbcM6jHbfd+CGsfB0enbD/1126FGJjbV3Nm2+GMWNgzhx7klA3i7iQGShRlADEBb5OCDyCK3aeAm41xhxmjEkDHgTmBjb9iIhIBGkeMpufZQ5w0aSBjffmryvo9r6JRISYOJh0JVz77861n/OD0OegB51zDqxfD9dcY3e6b9oE06fDKafAokXh6XMnRVzIBC4HqoG5gDfweTUQ3DL1APAm8AWwI9BGx0OKiESgfum2Slxdg79NKaMhOckcOyQLgLeX7er2volElL7jbH3NOwrsx5mlcP5TbdstfxUeGg615VBdEnqEsk8f+N//hSVL4BvfsNcWLIBjjoHbbz+Uf4oWIi5kOo7znOM4JsRjc+C+z3GcmxzHyXEcJ9VxnAscxyl0udsiIhJCZlIsqQl2+f/WfW1LAJ0wPBuANbvL29wT6ZVimpVhyBkRuk1tKdyfDw8OhrkdhMbx42HuXHj7bTg8UK9z4sT224dZxIVMERHpOYwxDMyyO7K3FrUNmSPy7I7zjQUV+PyRWVJPxDX9joKxF9iPtFPn67Mn7JnoM9Nhy6dt7xsDZ50Fy5bBX/4CFzY7Y33HDnjuOfAfmjXRCpkiInJIDQqEzJ0l1W02+AzPs7u1axv8bC9WsXORFjweuOhZuGY+/NdumPFix+2fPQO2fRH6XmwsXHyxDZ1Bd9wBV11lp9GDx1aGkUKmiIgcUsGRzAa/w67SGhzHobSqHsdxGJKThNdjf+it3xPZJ/+IuCo2AUafDZfN6bjdM6fDPXnw4kVQ28F7yueD4GlbX34JU6fCBRfAhg1h67JCpoiIHFIp8TFkp9h1Zqt2lTFvXQGz/v01C7/eR3yMl8HZNoRu0DnmIvs3/PTQG4Ka89XC+nfhvTvbD5peL7z6qh3BDK7TfP11W/Lo5z+H4uIud1UhU0REDrmx/e2hbKt2lrFkawkAS7fbj8Nz7ZS5RjJFOunIGYGd6Hs7brfoGbh/ACya1X6bKVPgiy9g9mzo3x/q6+Hhh2H4cFsKqQsUMkVE5JCbMDCTfukJLa4F9/mMCBwpuWGvdpiLHJCYTp7k89aNMOsMqK8Ofd/jgSuugHXr4K67IDERRoywQbMLFDJFROSQ83gM/zG2L/GxTT92qut81NT7GneYb9hbgePCqSQiUe30X8PgE+HmTXBj23PQG239FO7rC/9+FHZ+FbpNcjLMnGlHMGfNarlJ6LHH7A71A6CQKSIi3SIzOY6rTjiMK44f3HituKqucYd5ZZ2PpdtL3eqeSHQ66Qa46h1Izob0fPjFOrhxFRx1aej27/0KnpoKjx8D6+ZCQ23bNgMG2LWZQcuWwQ03wIQJcPXVsHt3p7qmkCkiIt0mMc5LVnIccTH2x09xZT2j+6UxLDcZgEfeX+dm90SiX2ofSB8A5z0BWUPbb1e4Dl66GN67a/+vWVoKQ4bYeppPP22n0u+/H6rbmX4PUMgUEZFuZYwhK9nuNi+uqsPrMfzs9JEAfLS2gC+3dn1Xq4gAV/1z/20+/yPMe6DjNiefDKtWwW9/C+nptvTR7bfD6NEdPk0hU0REul1mUixgQybA2eP7MTKwAWjO4u2u9UukR0ntA3eVwMk32ZODfvBh6Hbz7rcnBn29oP3Xio+HX/zC1tH88Y9tCaQtWzr89gqZIiLS7TKTAiOZlTZkejyGUw/PA2DFzjLX+iXS4xgD035lTw7KnwTff7/9trPPhienwMYPwXHA72vbJicHHn8cli+3x1V2IKaLXRcRETlgmYHp8pKqevx+B4/HMK5/OgBrdpXR4PMT49U4iEjYDTwGzn0cakrg3Tva3t+1FJ4/337ujYfL/gqHTWnbbvRoePvtljvQW9E7WEREul1wJLPB71BSXQ80FWyvbfCzsaDStb6J9HgTL4cTrrcF3fOPab+drxZmnwM7vjyob6OQKSIi3S4rOY7EOC8Aa3fbIuxDspNJDlxbuVOljES6xXdegUv/Cle+1X6bp0+Fly9tv5h7OxQyRUSk23k9htH9AkdN7irDceyU+ZjAaOaKHVqXKdItkrJgxDfgsJNtfc32rHnLFnN/elqnX1ohU0REXBGcHi+rrmfbvurANbsuUyOZIi5IHwA3rLCnCJ0WYr0mwI5Fdid61b79vpxCpoiIuCInJZ4+afY887V77JR5MHiu2lmG368jJkW6XcZAe4rQlJvh1g5KFP3mMCjc0OFLKWSKiIhrhmQnAVBQbo+2C45kltc2sK24yrV+iQiQmAH/tRuu+Tj0/ccndfh0hUwREXFNVkrTyT+O4zCiTwpxgdJFWpcpEgFiE6HfkXYn+sQrD+ipCpkiIuKa4PGSdQ1+ymoaiPV6GNU3FdC6TJGIc+6jcMmfO91cIVNERFyTmRTXWMt5X+D0n3ED7LrM5TtK+eeK3ewtr3GreyLS2qgz4ZRbO9VUIVNERFwT6/WQnmjPMQ+GzDGBdZkL1hdy7QuLuf215a71T0RaMQZOvd1On8/seLZBIVNERFwVnDJvHMkM7DAPen/1XhxHO81Foo1CpoiIuKopZNod5of3TWvTJrj7XESih0KmiIi4KhgyiyrtDvPEOC+njMxt0SZYR1NEoodCpoiIuCo7OR6A2no/RYEp8ycvn8Qnt51GfmYi0HS+uYhED4VMERFxVW5qPKkJMQAs/NoeVZcQ66V/RiKj+thyRus0kikSdRQyRUTEVV6P4bih2YAdsWy+/nJk32DIrHClbyJy8BQyRUTEdaP7pTWWMlrRrAh785HMmnqfK30TkYOjkCkiIq7zegzD8lIAKChrNpIZCJlVdT6Ou/8Dlm0vcaN7InIQFDJFRCQi5KXaDUAFFbX4/bYu5uF9U7lgwgA8Bkqq6nljyU43uygiB0AhU0REIkJuIGTWNfgpqa4HwOMx/H7GUXxn8iBAG4BEoolCpoiIRISspDhivfYg89bnlWuXuUj0UcgUEZGI4PEYclICU+atTvgJrs3cU1ZLaVV9t/dNRA6cQqaIiESMvDQbMveWhQ6ZAOv2ajRTJBooZIqISMTITUkAYG95LY7jNF7PTI5r3Bik039EooNCpoiIRIzgSGZNvY/y2oYW90b11bpMkWiikCkiIhEjOzkOjwls/mk1ZT4iz4ZMjWSKRAeFTBERiRgxXg9ZKXFA2x3mw/KSAdhSVNXt/RKRA6eQKSIiEaWxKHt5LTtLqqltsMdJDspKAmB3WY2OmBSJAgqZIiISUYIhc1NBJX/5YhsfrN4LwMDMpMY2O0qqXembiHSeQqaIiESUvLSEFl8H12D2z0jEY5drsnWfpsxFIp1CpoiIRJScwJrM5qrqGoiL8dAvPRGAbQqZIhFPIVNERCJKfIyXARmJLa4VVdQBMDBLIVMkWihkiohIxDnnyP5cOnkQqQkxABRV2pAZ3Pyj6XKRyBfjdgdERERaS4zzkhjnJTsljvKaBooqbM3M4Oafrfu08Uck0mkkU0REIlZ2st1pHpwuH5RtQ+b2fVUtjp0UkcijkCkiIhErO7AJqLDSnmWeHxjJLK9tYNa/NytoikQwhUwREYlYOSl2JLO23s9zn2ymf3oCw3LtyT/3vLWKuSv3uNk9EemAQqaIiESs3JT4xs0+JVX1bCio4LUfncj4AekAvLtqt5vdE5EOKGSKiEjE8ngMF07KZ0z/NMCuzUxPiuU/j+gHwKcbizRlLhKhFDJFRCTiBY+aLAzsMj9+aDYAu0pr2FKkckYikUghU0REIl5wbWZ5TQO1DT7G9k9rrKH5ycYiN7smIu1QyBQRkYiX3eyoyX2VdcR4PUw+LAuAzzYpZIpEIoVMERGJeElxMSTGeYGmmpkTBmUCsGJnqWv9EpH2KWSKiEhUyE4O1MwMrMsMbgb6urCSqroG1/olIqEpZIqISFQIrssMjmSO7WdDpuPAmt3lrvVLREJTyBQRkaiQFRjJLK6yITM3NZ6cwFrN1bvKXOuXiISmkCkiIlEhLTEWgIraBnx+B2MMowOjmat2KmSKRBqFTBERiQrpgZDpOFBeUw/AmGDI1EimSMRRyBQRkaiQlhCDMfbz0upAyAxs/lm9q4wGn9+trolICAqZIiISFWK8HlLibQH2YMg8Ij8DgJp6P2v3aPOPSCRRyBQRkagRXJcZDJlDspPISLLXlmwrcatbIhKCQqaIiESNtISWIdMYw5GB0cyvtpa41CsRCUUhU0REokZw809ZdVPx9aMGZgAayRSJNAqZIiISNdJbTZcDTBiUAcDGggrKaupDPU1EXKCQKSIiUSM9sP6ypt7H2sApP+MHpAO2tNF6bf4RiRgKmSIiEjVyU+JJTbA7zN9ZvostRZVkJceRHOcFYHtxtZvdE5FmFDJFRCRqxMV4uPiYgWQHjpNct6cCYwz5mUkA7ChRyBSJFAqZIiISVdISYhnZJxWA7cVVAAzITAx8rZApEikUMkVEJOoMzLIjlyVV9ZRW1zMgw4bMHQqZIhFDIVNERKJO37QEYr32jMntxVXkB0YyNV0uEjkUMkVEJOp4PaZxinzbvupm0+VVOI7jZtdEJEAhU0REolK/dBssC8prGjf+1NT72VdZ52a3RCRAIVNERKJSbmo8APsq6+mbFt94XVPmIpFBIVNERKJSMGT6A9Pj8TH2R9qGvRWu9UlEmihkiohIVEqNjyEh1hZhL6yoY0z/NADuemMly7eXutk1EUEhU0REopQxpnE0s6C8lvsvGE9OSjzltQ08/P46l3snIgqZIiIStRpDZkUth/dN47YzDwfgXxsKqaxtcLNrIr2eQqaIiESt3JSmkUy/3+G0w/PwGKhr8LNgfaHLvRPp3RQyRUQkavUJ7Cqva/Czr6qOrOQ4jh6cBcB7q/a42TWRXk8hU0REolZWchxxgV3lu0trADj18DwAFm4ucq1fIqKQKSIiUcwYQ9+0BKApZI7ulwrA9uJqaup9rvVNpLfrsSHTGPO+MabQGHNHq+uXGWM+DTymudU/EREJj77pNmTuKrMhc1huCgCOA18XVrrWL5HerseGTOC7wE3NLxhjMoCfA1OBs4GHjTHe7u6YiIiETzBkFlXUUtfgZ0BGYmNh9o0FKswu4pYeGzIdx9ke4vJkYIHjOLWO4xQBm4Fh3doxEREJqz6B6XLHgaLKWjwew9DAaKZO/xFxj2sh0xhziTFmgTGmzBjTppiZMcZrjHnIGFNgjCk3xswxxuR08dtmA8XNvi4OXBMRkSiVHOcl1msAKK+xP06G59mQubFA0+UibnFzJLMYeAK4oZ37twHfwo4+5geuPQ9gjEkxxiwK8bhrP9+zCMhs9nVm4JqIiEQpYwxpibEAlFXXAzAsNxmAjRrJFHFNjFvf2HGcuQDGmKntNLkauNtxnE2BdrcAG4wxQxzH2QwcfRDf9nPgfmNMPJAMHAZsPIjXERGRCJKaEENRRR1lNcGQaUcyNxVW4Pc7eDzGze6J9EquhcyOGGPSgUHA4uA1x3E2GmPKgCOwayn39xqzsKOg8caYyY7jnOM4Tokx5n+AeYFmP3ccp9P1LYwx2QSm18ePH09FhX5DlkOvqqqqxUcRaSsWH3V1dRSUVFBRkUxuog2VNfV+tu4tJicl7qBfW+9BkYNjHMdxtwN2JPN9x3Fiml0bCGwFhjqO83Wz61uA/3Ic54Xu7mfg+88EglPytcAKN/rRCRlASQ/4fuF63a68zsE890Ce05m2XqAPsAdQ0T8rg+79O36gMoj+92C4XrMrr3Mwzz3Q53Smvd6DbWUQue/BDKL//RfO1+3K63TmuYMdx8kNecdxHFcf2HJCDa2uZQAOcFSr66XAuS72NRsYGXj81e3/dh3086We8P3C9bpdeZ2Dee6BPKczbQN/3xxgZHf+f43kR3f/HY/0/h2K79cb3n+dba/34KH7+9ET+qafge0/InK63LHT2luBicASAGPMUCANWOZiv4oIbBQyxvzZrX50wpwe8v3C9bpdeZ2Dee6BPKe7/1/1FJH+360nvAd7w/vvYL+HRPZ/t57w/gvn63b3e7CRa9PlgSLoscAU4B0gJXCr1nEcxxjzX8AVwBnYYPcMkOo4zhlu9FfELcaYkcBaYJTjOOvc7o9Ib6P3oMjBcbOE0eVANTAXu96lOvAYHLj/APAm8AWwI9Dmsu7vpojrioBfo3JbIm7Re1DkILi+8UdEREREep4ee6ykiIiIiLhHIVNEREREwk4hU0RERETCTiFTRERERMJOIVNEREREwk4hU0RERETCTiFTRERERMJOIVMkyhlj3jfGFBpj7nC7LyK9iTFmkjHm38aYj40xHwaOPxaRABVjF4lyxph84HQg33Gce93uj0hvYYzpB1Q4jlNujDkL+LbjOJe73S+RSBHjdgdEpGscx9lujHG7GyK9juM4u5p9WQc0uNUXkUik6XIRlxljLjHGLDDGlBlj2vyQMsZ4jTEPGWMKjDHlxpg5xpgcN/oq0tOE4/1njEkGHgB+2139FokGCpki7isGngBuaOf+bcC3gMlAfuDa84e+WyK9Qpfef8aYOOBV4F7HcVYeum6KRB9Nl4u4zHGcuQDGmKntNLkauNtxnE2BdrcAG4wxQxzH2dwdfRTpqbry/gO2AS8BrzmO87dD3VeRaKOQKRLBjDHpwCBgcfCa4zgbjTFlwBHAZmPMLOwoS7wxZrLjOOe401uRnqUT77/jgDOAHGPMZcByx3Gud6WzIhFIIVMksqUFPpa2ul4SvOc4zve6s0MivUiH7z/HcV4AXu7WHolEEa3JFIls5YGP6a2uZwBl3dsVkV5H7z+RLlDIFIlgjuOUAFuBicFrgYLPacAyl7ol0ivo/SfSNQqZIi4LlEhJAOICXycEHsHil08BtxpjDjPGpAEPAnO16Uek6/T+Ezl0FDJF3Hc5UA3MBbyBz6uBwYH7DwBvAl8AOwJtLuv+bor0SHr/iRwiOlZSRERERMJOI5kiIiIiEnYKmSIiIiISdgqZIiIiIhJ2CpkiIiIiEnYKmSIiIiISdgqZIiIiIhJ2CpkiIiIiEnYKmSIiIiISdgqZIiJRyBiTa4x5wRgzy+2+iIiEopApIhKdHgK2YI9FFBGJOAqZIiJRxhiTBFwAvAVsd7k7IiIhKWSKiESfqcBS4GhggbtdEREJTSFTRCT6TAE+Bc4F/upyX0REQlLIFBGJPscClcAI4B2X+yIiEpJCpohI9BmMHcV8yHGcBrc7IyISinEcx+0+iIjIATDGVAI7gXGO49S63R8RkVA0kikiEn0agJsUMEUkkilkiohEEWPMlUAaUG+MGWuMucztPomIhBLjdgdERKRzjDGJwHTgMuB3wC7g+652SkSkHVqTKSIiIiJhp+lyEREREQk7hUwRERERCTuFTBEREREJO4VMEREREQk7hUwRERERCTuFTBEREREJO4VMEREREQk7hUwRERERCTuFTBEREREJu/8HtKUuOfkG8joAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "if output_nside < 4096:\n", "\n", " for ii in range(1):\n", " plt.loglog(ell, clin, alpha=0.5, color=\"C%d\" % ii, label=r\"$ \\beta_s$ PySM2\")\n", " plt.loglog(\n", " ell, clin * 1.5722501 ** 2, alpha=1, label=r\"$ \\beta_s$ PySM2 rescaled SPASS\"\n", " )\n", " plt.loglog(\n", " ell,\n", " clout,\n", " label=r\"$ \\beta_s$ PySM3 \",\n", " )\n", "\n", " plt.grid(True)\n", " plt.plot(output_ell, cl_spass_fit, \"--\", label=\"Power law fit \", color=\"red\")\n", "\n", " plt.legend()\n", "\n", " plt.ylabel(\"$ C_\\ell $\")\n", " plt.xlabel((\"$\\ell$\"))\n", " # plt.xlim(10, 55 )\n", " plt.xlim(2, lmax)\n", " plt.ylim(1e-10, 1e-2)" ] }, { "cell_type": "code", "execution_count": null, "id": "7298c15e-7c96-4dc8-9ca2-5c5e6f50441a", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "condanamaster2", "language": "python", "name": "condanamaster2" }, "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.8.0" } }, "nbformat": 4, "nbformat_minor": 5 }