{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "# Matter shell definition\n", "\n", "This example demonstrates how matter shells are defined, and their angular power\n", "spectra computed.\n", "\n", "The the angular power spectra are saved here, so that they can be reused in\n", "other examples without recomputing.\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Compute\n", "Here we define the shells for these examples, and use CAMB to compute the\n", "angular matter power spectra for the shell definitions." ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np\n", "\n", "import camb\n", "from cosmology import Cosmology\n", "\n", "# almost all GLASS functionality is available from the `glass` namespace\n", "import glass\n", "import glass.ext.camb\n", "\n", "# creating a numpy random number generator for sampling\n", "rng = np.random.default_rng(seed=42)\n", "\n", "# cosmology for the simulation\n", "h = 0.7\n", "Oc = 0.25\n", "Ob = 0.05\n", "\n", "# basic parameters of the simulation\n", "lmax = 1000\n", "\n", "# set up CAMB parameters for matter angular power spectrum\n", "pars = camb.set_params(\n", " H0=100 * h,\n", " omch2=Oc * h**2,\n", " ombh2=Ob * h**2,\n", " NonLinear=camb.model.NonLinear_both,\n", ")\n", "\n", "# get the cosmology from CAMB\n", "cosmo = Cosmology.from_camb(pars)\n", "\n", "# shells of 200 Mpc in comoving distance spacing\n", "zgrid = glass.distance_grid(cosmo, 0.0, 1.0, dx=200.0)\n", "\n", "# triangular radial windows, equivalent to linear interpolation of n(z)\n", "shells = glass.linear_windows(zgrid)\n", "\n", "# compute angular matter power spectra with CAMB\n", "cls = glass.ext.camb.matter_cls(pars, lmax, shells)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Plot\n", "We can do a quick preview of the shell definition." ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjkAAAGxCAYAAACA4KdFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAACxvElEQVR4nO29e7hdVXk1Pvbe55aT6wm5cDHlptwKBAklX8D8gl+jtFps5fMpH1hAithaqEI+LyBICqihYCJW0RQUpa0ttBZpn8JDq6nRoilUSKzKrQrIJRcSgZwkZ++91ppz/v7YZ+3bWbd5eedae585nicP4Zy95lpn5+w1xxpjvO9bEkIIODg4ODg4ODj0Gcp5X4CDg4ODg4ODAwUcyXFwcHBwcHDoSziS4+Dg4ODg4NCXcCTHwcHBwcHBoS/hSI6Dg4ODg4NDX8KRHAcHBwcHB4e+hCM5Dg4ODg4ODn2JgbwvIC9wzrF9+3bMnj0bpVIp78txcHBwcHBwyAAhBPbt24dDDz0U5XKyVjNtSc727duxZMmSvC/DwcHBwcHBQQEvvvgi3vCGNyS+ZtqSnNmzZwNovElz5szJ+WocHBwcHBwcsmB8fBxLlixp7uNJmLYkJ7So5syZ40iOg4ODg4NDjyFL1MQFjx0cHBwcHBz6Eo7kODg4ODg4OPQlHMlxcHBwcHBw6Es4kuPg4ODg4ODQl3Akx8HBwcHBwaEv4UiOg4ODg4ODQ1/CkRwHBwcHBweHvoQjOQ4ODg4ODg59CUdyHBwcHBwcHPoSjuQ4ODg4ODg49CUKQXK+//3v45xzzsGhhx6KUqmE+++/P/WYzZs349RTT8Xw8DDe+MY34utf/zr5dTo4ODg4ODj0DgpBcg4cOIClS5fi9ttvz/T65557Du985zvx1re+Fdu2bcOVV16J97///fjXf/1X4it1cHBwcHBw6BUUYkDnb//2b+O3f/u3M79+48aNOPLII7F+/XoAwPHHH4+HH34Yn/vc53D22WdTXWau+NWvfoX9+/fnfRmZBqK5a5ge15D3+d01uGtw1zAV8+bNw8BAIbb2QqAn34ktW7Zg9erVHV87++yzceWVV8YeU6/XUa/Xm/8/Pj5OdXnG8eCDD+J3fud3IITI+1IcHBwcHAqMo48+Gj/96U8xMjKS96UUAj1Jcnbu3InFixd3fG3x4sUYHx9HtVrFjBkzphyzbt063HDDDbYu0Sjuv/9+CCEwMDBQKIZeJNJVpGsB3PUkoUjXArjrSUPRrschGb/4xS/wP//zPzjppJPyvpRCoDg7JjGuueYarFmzpvn/4+PjWLJkSY5XlB2PPvooAODee+/Fueeea27h8R3AnEPMrQfglYlXsGh0kdE1AcDftQuDXcTWFPa/VsOsMZqnnvE9uzFnwUKStQEgeL2GgXm0T2yvv/465s2bR3qOavVlzJhxGOk5QvyyWsfhM4atnOuV8RrmzBjEyGDFyvkAAK8+C8w/yt75JvHL8V/i8DmHWz8vANSfew7DRx5p/bwhAfSqAbwaw6wxO79XcVi6dCl+8pOf4KWXXnIkZxKFCB7L4uCDD8auXbs6vrZr1y7MmTMnUsUBgOHhYcyZM6fjT6/g+eefBwCccMIJZheum7fsxgnWBAC2dy/JukII1CcCkrUBoH6ANkfFCa8daLw/1WqV9BwAEAR7IQQjPw8AvOrbOQ8ATHgME5698wEA9u+2e75J/Kr6q1zOCwDB7nx+5lKphFKphMDjCDzW/P+8/hxy8MEAgJdffjmX96OI6EmSs2LFCmzatKnja9/+9rexYsWKnK6IDpzzZn5obGzM3MKBB3gHzK03iX3+PuNrAgAnylCxgMMn3IRq+2nejxDUJCcIAnieR3oOAGC8CsboyRQA/Mqnfc/aUfUZJjx75wMAHMiJ5NTyIzlsz57czg0Afp3Bq1kmsxE4eFFDNXYkp4VCkJz9+/dj27Zt2LZtG4BGifi2bdvwwgsvAGhYTRdddFHz9X/8x3+MZ599Fh/72Mfw1FNP4Utf+hL+/u//HldddVUel0+Kffv2NSXRuXPnmls4qAFBPf11kiBTcohITuBxMI+TrA0AVXKS45Ou7/s+fJ/2HADAWQ2M18jPI4TAqxZJzoTHULWt5EzsASznaIQQeLX6qtVztiPYsyfX7JBXC+DXLJPZCCxasACAIzntKATJ+dGPfoQ3v/nNePOb3wwAWLNmDd785jfj+uuvBwDs2LGjSXgA4Mgjj8QDDzyAb3/721i6dCnWr1+Pr3zlK31ZPv76668DAAaHhs2m5YN6g+gYxrhHZVfRKTmBT0dyel3J8X0fQUB/82asCm5ByalxgQOM7t97yvl8hqpFewwAUN8P+HZUsRA1VsN+P58WF0IIsP37IWr0JDkOfp3Br+ev5CycVPsdyWmhEMHjs846K5GFR3UzPuuss7B161bCqyoG9k5mUWbONpwhIlJyJvwJBDzAQNnsrxY/cACCMZQqZgOcgccREG1CnDPUJyZI1gYAEXAIYpUgCAIrSo4tu6rKOSYskpwJj6FSttw3xZ9o/BkatXbKalBFndXBOEOlbDFkDTTIDePgExMox2QyqeHXimFXLRhrqP2O5LRQCCXHIR6hkjM6c7bZhZkHBOY3lTqrw2PmMxyiXiN5UmM+ByNScpjnI/DME8kQIuAQhCoUYNOuqoNxCySHcdQYt2ZtVH3LdpVfAzgDPLuqyoTfIPMTAR2pjwOffJDghA8UaSiKkjN/sqDmpZdeyvlKigNHcgqOkOTMmGVYyfGrJEqOxz3UmHkywut18Lr56w18RmZX+V4dgnOwgIYkCJ9DcAER0BEdWySHcTt2VZVzcAhUuSWS4wV2q6smyQY8uxt+dfKBKSQ7NsEnq//4AfOFFFnh1QJ4BcjkHDRJcn71q1+hlqN9VyQ4klNwhHbVyCzDSg5BJocLDo95xpUc4XkA4xAkJIcjIAoeB5NVST7BdQNoWlWUao4tu4ozO3ZVaFVVLVlWVc9ydVVYMWmZbIQKTi5KzoFiKDmBx8Etkec4jAyUMTQ4CADYvn17rtdSFDiSU3CESs7w6CyzCxNkcuqs3vFfU+CTZIET2VVUmZyQ5AREJdihgiMIg61h8JjS3hFCgPG6leqqkNxMcDskZ8J28Lip5Ni1q5pKTi52VYPYhWQnD/h1BgiRq2UlOIdfr2PhfBc+bocjOQVHqOQMj5pWcmrGlZxQwakZXjdUcKiUHKpMTpjHISM5Pu/4LwV834cQglTN4dwDBLdWXQXYUXIYF6j73K5dFSo5lu2qZiYnD7sqzOTkZFf5HoOY/L3Ks4zcq9UACCycPx+AIzkhHMkpONqVHM9k9iKokyk5xu2qSQWHhORM3qAoiA6bJAZU4WNbJKf9vxTgk4FjG0rOBGOT/6UnOaGCU/OZvR4uIcmwTDZCJadKUMyQhryDx35bVVWeFVZ+rfHeh0qOCx834EhOwREqOTNmzUE9MPgBCiarMJi5zYvMrqqHdpV5ssAmiWNAEN7tByUn7JFD2SsnzOLYCR5PKjkW7KqwqkoI2LOsvJztqhyUHJG3ktNmUfk5kpx6tfE+LDzIKTntcCSn4AiVnJGZs1E3uRGHaovBpmH1gIbkiEmyIAgUkTB0HBBYCq1MDq2SwwntEDtKTkPBsdInJ8zk2FBy2v5drFlWfj52Va5KTlhdlZeSU289AHj1/Owqv9qp5DiS04AjOQVHs4TcNMkJyY1By6rOiUjOpF1FFTxu/69JNEkOVXXVpLJHXULe/l8KsMmWA0IwcE7XVwhoKTg2Mjnt6o21Xjnh59q3q2rk2idnMnAsPA/CQiVgN9rVmzyVHK9pVzklpx2O5BQcHXaVSck7JDcGQ8JhFofKrhJ187ZPWFlFUUYeKji+T2tXgaqZIWNgkxkWG5mcxjlpczk2q6vaS8etKTlh8JgFJH2wIk/JPLDJCfK27SrBeVPJAfKxrLwOuyrH4LFTciLhSE7BQabkhOTGpJJDlMkR9VrHf42tKwRY0MhoUJSR+9Ql5KFdRVUd1pbDoVVy2kgOYddjJgRqOSk51nrleAei/06IdmLjcx++wZxfGkS12jGMNA/LqmjB4wWTJGf79u3gllolFBmO5BQc7c0ASTI5BEqO6eqqsNOxabuK+bx5g6Toesx8O5kcqj457cSGlOS0VVVRho/biY3tTI4Vu4rzzs+zJZLTncOxaVm1qzhATiSnXcnJsU9OaFctmDcPpVIJvu9j9+7duV1PUeBIToEhhOhUckxtZoHXqKwCSJQc831yaOyqdmJDmsmhrq4i6thsi+TwNouK0q5qr6iyUl3lWw4eB52qhq0y8ikkx6Jl1W1P5WFXtQePA4+BWxwA247QrhoYGMD8uY3xDs6yciSn0KjVavAmN8gZM+eYU3LaSYhBQtLsk8MN98khsqvaiQ2FkkNJckTAmw3IqILH9pSctkwFoV3Vrt4wIVAj3ozaic2EjRLy7ooqS2Xk3cqNzQqrbuUmbyUHyMey4px1jI9Z4MLHTTiSU2CEVlWpXMbQjFFzJKfdTjJ4QwpLyD3mgQtzG0hoV4mAGa2eaM/hmC4hD3wfQoTl6ebtqnZiEw7qNI32TA5ln5xOJYfQrupSb6jVnFpHdZWFTE53RZWlMvJc7apukmNZyWFsasf0PCwrf7LbcQgXPm7BkZwCo9kjZ3QWyuWyuWaA7b1xCOyq7r/rot2m4gYtK0q7qp3YCM4RGFZCuhsAUqg5eSg5lF2Pq6yTCFLnciZs98npJjXTwa7KWcmJKhnPYxq515VNciSnBUdyCoxm+fjMxtyquqmNuEPJMR887v67LtptKpOWFaVd1W1RmVZzppAcgk20m+RQjSZoV3JsBY8BWiVHCNGh5FghOd3kwpZd1XVeu0pOV/DYspITpdrk0StnKslxdlUIR3IKjPZuxwCIMjk0So6p8LFgDMJvPRmZnF/V3hvHuF01heQYzil1ZTwolJx2i0oIQWJZcV5v2nqAXbuKUsmp+qwjA+wFHAF1ILW7miovu2oaKTlRqk0uJKfW+XMvcPOrmnAkp8BobwQIwJxdRRA8DnjQbAgGmFNyukmNyflVlEoO86lJTreSQ2tXRf2/CXRXU1H2yekmNZS9cqJKxsnDx93kwhLZ6FZu7Co5XdVV1RqExd4wkUpObpmcFhY5JacJR3IKjJaSMwuASSWnHv13DXSTGlOZHN5FckzaVe3BYx5wcIPhXb/rusntKoLqMCskp4vUkNpVFpWcKHuqRm1ZTVFy7NhV3UpOntVVEMKqmhOZyclhflU4nDOEy+S04EhOgdHqkdNQcrwCl5B3kxpTJKdbyTFqV3URA5OWVTepIVdyiO0qgIbk8ClKTp0s+9NdMk6ZyalFqDbkuZxu5SbwGuMdCBHwYMoDji27SjAW2TsrnGVlA1GqjVfNQcmJCR7v3bsXB3Kazl4UOJJTYLTsqkYmh3Fhhuh0KzkGNhUyktMlw1LZVQDADBKF7mqqXg8eR/2/CUypphKiOZXcNLqVG9tKDjnJicrgEA/qjLKmmGDGR7tEgXeNdGh+fcLeph5tV+VfXTVzdBQzRkYAODXHkZwCozt4DBjK5bSrN0IYsayoSE53ybgwSBamKjkGSU538NjwJHLR9XvQq3ZVlD1FFT6uWeyTU41Qcqo+8eYXRWiIRztU/eh/KxtqTpxiY1XJiQgeM58bfWDKgnCkQzucZdWAIzkFRquEfE7za0ZyOd2kxoBlRabkeN3BY5oS8qj/10G3cmN6Enl30Ng0yYkqGacJHk+9OVMoOR7n8IW9PjmRwWNKJYcFDXuqG8QVVnH5Gxvh4zjFxqaS48WEjG1WWHHGIpXiha7CCoAjOYVGd/AYMEVyujYRA0oOVfC4264ym8lhif+vtTZ1CXlAS3KiysVpSsinEhrGzG+QUZVUpNVVtjM5cbZUDnYVYEnJiQkY2+qVI7iIVX9tNgTstqpCuF45DTiSU2B0l5ADMDOks4eUnCl2lSGSI7iYotyQ2lXk1VVmN9Ao1cZKJgc0dlU1onLOFwIekWUVRWhIJ5HHKTbUdlWMkmOjwkrEkRxL1VV+ncXmGW2WkXf3yAnh7KoGHMkpMNonkIegUXKKS3Km2lVm1g0i3keTvXKsV1cR2FVZvqaLyEwOgV0Vp9pQqTn2lZw4ktPHdlWMgmGV5MTA5pBOp+Qkw5GcAqNJcmaZJjnmlZxuu8pYM0Aiu4pFqDYmMzlsSnWVwTEXwdSBnDbsKhvNAAGaXjlxIeMJIiUnaiDnBOWQzliSQ9srJ86WshM8jsnkWLKrkkiOzQqr7kaAIZyS04AjOQVGaFeNzDRoVzEf4F1rGMjk1Lo2K1NjHbrtKm6oGWCUamMqk8M5AwvoSsijeuIILoz2yrFnV0UpOQQkJ1bJMd+TRwiBagSJjuqdYwxxig0x2cg3eJyvkpOUu7EZPO5uBBjCkZwGHMkpKBhj2LdvHwDDdlUU+aBQcrhnpKnblA7HjEMY2GyjCI2pTE6UamNyEnncCIdeJDndzQDjvqaLWCWHmd+M6gEHj/jd95kwN5plyuIx6gVxJqeIwWMwbrQKMw5Fsau6GwGGWDBpV+3YsYOkaKBX4EhOQTE+Pt78u9HqqijVhqBPjhACPtffFKPsqe5RDyqIsqZM2VVx1pQpNSeOzJicXxVFaIQQRokO5x6EmLoZUASP48rFowLJukgKGJOFj6dh8DhJsbFhWRXFrorqkQMA8+fNRaVcBuccu3btsnY9RYMjOQVFmMcZHB7BwOBQ8+vaT4JRqk1MQy8ZRGVwui0sFXTbVcDUnI4KKO2qODJjKpcTV0llssIq7snPJMmJyuMANHZVHMmh6JWTNIgzKpBsBHHKSVADiHJHXHDUYx6QJoIJsvEcACB8HyLh82TDskqypKwGj2NITqVcxkFj8wBMb8vKkZyCotUIcHbH1+u6aoMlJQfQDx8LISJvZCbCx8yyXZX0dVnEhYyp7SrAbK8cHkNmKOyq7m7HISiqq5LUGrIKqzjFRgiyXE41qEIgmsgIIUjVnDQSk7eSwwNutJAhCV7Ce7HA5XIcySkqokY6AMXN5ESRHN0yclGPnqtlwq6KUnJYwI08fZLbVXEkh9iuSvq6CuKUnDgbSwexSg6BypGLXZWkxhJZVmkkhjJ8nEpyLCg5aQ3/bEwjZ4E/pcihHQvHXBm5IzkFRat8fE7H1z0Ku0qT5HgsOmQcJ2VnRZxiY8SuiiAEggvwgJLkECs5Bp8cbZCcpPENpnM5cdkbEiUnwZIiU3KSOhsTdT1OCxeTKjkxYdvm93NWcgA7FVZxPXJCuNEOjuQUFrF2lbaSE7HRapKROMVGV8mJa/xHFTwGzORyclNyei6TE3+DNklyhBCo2czkJPTDIemVE3iN2VVxIGoImKrkEFZYpQ3hpB7SGXhsSq+qbtjI5cT1yAnhysgdySks6OyqiBuTppJDRXLiJo6byOTEdTc20fU48GOCx6ZGUsSoeaaUHM65HZKTEDCOy+uooMYFeEx2hGISeVI/HBK7Kk2pycmuyjWTQ2xXZSEwURPKjV9HTI+cEAsPcnaVIzkFRWtuVSfJYVzA13n6jFJtmK9VgREXMNYmOTFPKSZ6YEQFjwEz4eM4JcfUJHJquyqJyBi1qxICxiaVnCQiU+ccgeEy8iRLisSuSlNqqOyqlMwNqZKTMmmc2q7KMpsqbkK5ScRVVoVwSo4jOYVFnJIDaKo5caqNxlNXXKm4tl0VUT4OACLm6zKIV3Io7areIDlJFVRmlZwEkmNwflWaJWVazUnM5FCUkKeRCSq7KqX1RD8Hj7OQnGJkclpKDmVJf5HhSE5B0crkzJnyPa3RDnH5G41cDpmSE2tX0fTJAcw0BIzvk0ObyeGGNlB7Sk6CXWVSyUkjOYZzOcnVVQQWRpodRWRXpSo5OZIc4XlGOqPHIYsVlVZ9ZQJZlZwDBw50NJidTnAkp6BoTSCfNeV7NEqOOnEgy+TE2VWa2RbOBXjMe0hpV/WKkpNEZEz2yUlUcizZVYD5MvJkkmOmTUEH0pQcIrsqNZNjoMloHEQGpYbSsuoVJWdkeBizZ44CmL4VVo7kFBRNu2pWhJKjRXLilBx1khOr5GhWbcWRGRFTdZUVSWpN0e2qqAnkrRMkfE8CRVByTNpVNpWcesASMz5cCP3igW6k2VH92CcnZXMHaC2rTMFjG5mcDO/DgvnTO3zsSE5B0QoeR5GcYtlVZEpOHMnRtH2SiIyuXSWEIC0hT+pqLISZrsdFyORwZm6DSlVyDJKcWgYl0Hj4OLW6aiKyqaYOsnQ0rgU1cGG+eo3X6xB+uqKYt5LDGUdA1RcpvI4UuwoAFk3z8LEjOQVFol2lsxETKDlxZEZ7rEMMydGtrkqypHRLyJnvA3Gt7jlP7E6aBWmWlAnLKonIJJWXy4BzH4LHrxPXDVkFqUqOQbtqIsPma7xXTpqSI7iR+XTtqLF0AiNAM9ohq0JDqeRkHcBJ2Ssn8JO7HYdY6JQchyIiMXis+rTOfCBuY/EJ7Crt6qqY47kA17B+Eu0qzUxOmiWla1nlTXKyfD8LkrodA2aHdKZNGjep5GTpg2O8V06WUm3D5dxZyQtFGXmWPA5ATXKy/RtSTiNP65ETYrrPr3Ikp4AQQqSUkCveJJPUGg0lpxZzLBMMPlffEJOyNzoNAZPUmrj+OZnXTrGkfM3QdFpXYxNdj9OUGhMkJ02psVldZZLkZLGijNtVWTI3hnM5mUkOQS4ns5JDZFcxln34JqWSk8WqAlyvHEdyCoiJiYnmRtPdDBDQsKuScjcEJeRp30tDUvZGZ35VcvDYKTl2lJzkG7QQHJybKblPIzEmq6uSuh03z2eS5AiRzYoyTHKyKjT9aFfJVE1RVlhlCR0DLbvKVVc5FAahVVUuVzA0Mjrl+8p2FZGSk2RLxak8WZBUKq5TRp4UPNa3q5KvSzd83C8kJ0uJuKlcTlrmxmR1VRYCU82Q28mMoNbI3KShj+yq3EmORNUUZa+czCRnmo92cCSngGi3qkql0pTvq9tVREoON6/kcM8DErIUVHYVZxxcoww7SCEA9EqO/pNjGokxETzOUiJuIpfDhEDdYnVVUrfj5vlMKjlZuxk7u8oYpJQcwjLytEaAIUK76pVXXoFnqFdXL8GRnAIibm5VCBolR21D4YLDZ/Gbomr4OI3E6FRYpfnpOmWf/aDk2MjkZMncmMjlZFFp6lyAGyqxzhIqNkpysjb6y8muolFysv1e8GoNgmAAqycRJi6CXTVv9mwMDgwAAHbs2EF2PUWFIzkFRFLoGNAY60Cg5HjMg4gpmQY0SE4KidGZX5XW8E+nVw55JidFxdMlOUEQgKdsDEbsqixKjgG7Kkt5uIAwZlllUXKMVldlVXLysqtyVHIgBIllVRi7qpbtZyuVSlgwNn3Dx47kFBCt8nGbSo7ahpJGYlRJTtxwzhA686vScjc64eM0EuMXXMnJQmDMKDnp/35p4eQsyGpFmQofW6+uykvJyUhe8szkAAA/QEByJO0qqsGYfkYlB5jeFVaO5BQQaUqOR0Jy1DZfKpKT1tW4d+2qYldXZcnb2AseG7CrMpIXU7mcLNVV9YBp5b46MM0yOUII8Iz9YQCaXI6MkiO4MDIPLwqexD0w7JUzHSusHMkpIJrdjiNGOgBAwAUClZtykLDBctZoFigJMpJDalflp+TQZ3L0VAJbSk6WULERksOykQkTXY99xjM9gAgBTBiaGJ/ZhjI8pDMryfGYhyChs7UsRK0GSNz78rarVF6fBYHngbPs7+uiaVxh5UhOARHaVSMRIx1CKFlWaZaUgmWVVj2lOqQzrURcx65KU3IKncnpF7sqQw8cE9VVWbM2JpScLHmc5vlMjXbIqtCwIPkhR+aUksTFpJqTZTBnx+snzCs5sjkbilxO1sqqEM6ucigUWnOropUcQJXkpGwsCoSETMlJUWpU++RwxsFTNrSi2lWCpU8ZF0xAaGzYWQiMkdlVmaqr9IPHWbM2JoLHMoHiLIM8M0Em8+LtN3JK2QZ/VYNzs2SVGdN2lYr9RFFhlbWyKsR0nkTuSE4BkVZCDij2yklTahRuRqlKjjLJSbGrEkY+JCGLFaVlV6X2yVG3q7KqNDpqThYCwxgDY+o3bs4D8AzjPnpOyZEgOVkGeWZC1kwOYKzCSpbkGFVyJIPEpu0q32PSE90p5ldlnVsVwik5BcDtt9+OI444AiMjI1i+fDkeffTRxNffdtttOPbYYzFjxgwsWbIEV111FWqa06mLgpaSk0ByVDayHlJy0pQarmhXZbGiVO0qzliqT64zidwGyclqRelYVlmrpoz0ycmq5BgIAstUTRmrsJLJ2hgKH8tWTBklOZL2k3GSo6DKUMyv0rGrqKq9iopCkJx7770Xa9aswdq1a/H4449j6dKlOPvss/HKK69Evv5v//ZvcfXVV2Pt2rV48skn8dWvfhX33nsvPvGJT1i+chq0Mjmm7SrzmRw6uyotk+MpfVgplZysKo3qkM5+ITks40wqxuvaN+TsSo7+RiSTyTHSK4dzuQcTQyRHWskxWEaet12lEiKmsKtkyscBNPvk1Ot1/OpXvzJ+PUVGIUjOhg0bcNlll+GSSy7BCSecgI0bN2J0dBR33XVX5Ot/+MMf4swzz8QFF1yAI444Am9/+9tx/vnnJ6o/9Xod4+PjHX+KipaSkxQ8VrGrekfJSR3bIASEQr4li0qjmsnJmrdRzeVkrZzSqbCyouRkVWiE0O6Vk5XkmFBypOwqEyTHPyBnnRgiG7LKjMkhnUJyc5d9fRpUQsQUweO6pF01NDiIeXMaD83TzbLKneR4nofHHnsMq1evbn6tXC5j9erV2LJlS+QxZ5xxBh577LEmqXn22Wfx4IMP4h3veEfsedatW4e5c+c2/yxZssTsD2IQaSXkgIKSwwIgrSKCQMkJeADG5W/oWYLFKvOr0rodN16jquRQk5xiZHIATSVHwobSLSO32SfHenWVbIauH5QcSWVGBEyrp1Y3lJQcghJy2eAxMH1zObmTnD179oAxhsWLF3d8ffHixdi5c2fkMRdccAFuvPFGvOUtb8Hg4CCOPvponHXWWYl21TXXXIO9e/c2/7z44otGfw6TaNlVBjM5WQgMQQk5oKbmZAkWq9y8shAY1UxO1m7GquFjkZHYFt2u4hlGOoTQGe3gcY4go9JRY1zbGpMhLkbsKlnSYiqTk6OSo5KxMWlZqZIc0zkYX+He50hOD2Hz5s34zGc+gy996Ut4/PHHcd999+GBBx7ATTfdFHvM8PAw5syZ0/GniPB9HwcmP5RGq6uyWFEEdlXW13QjreNx4zUKdlWG8s/AV9vwWMaNX1nJyVi6WnSSI0NcdOwqmbJwDqFtWWXpdhzCjF0lueGbqq6SVJDMBo8VSI7B8LFKvkZwYVzNyTq3qh3TleQM5H0BCxYsQKVSwa5duzq+vmvXLhx88MGRx3zyk5/EhRdeiPe///0AgJNOOgkHDhzABz7wAVx77bUol3uSuwFAR1YoUcmRtasyKTnyG0qWZn+yJEcEQeogSiC9K3IUMllRQoAFHAODFbm1qZWczHaV2g2Vc565NNxGdRWgZ1fJzqOqMo7Rivq9Q4a4yFhbsZApHwdyU3JM2VWCc/Cq/Gc+byUnPG5oxMx269fVpquHvXKm22iH3NnA0NAQli1bhk2bNjW/xjnHpk2bsGLFishjJiYmphCZSqWxIfV6eVyYxxkaGUWlEv+hkFdyspAcGiUni6XVDp6xB07W17UjqxWlMm+GPJOT8d9cVckJgiDz50enIaCMkpNlWnkcso50CKE7pFOGuHgBh6+bA5Id1ZBTJsfnPnyFkTHdENWqdI8awKySoxoiNllhpZLHAZySkyvWrFmDiy++GKeddhpOP/103HbbbThw4AAuueQSAMBFF12Eww47DOvWrQMAnHPOOdiwYQPe/OY3Y/ny5fj5z3+OT37ykzjnnHOaZKdXkaURIKCSycliV8ltKAEPwET6h7cmmavIYlXJvK4dWYLHgFoup9eDxzLqjF4JefabtE6vnJqCkqMKxoX0Z3LCY5g7Q+M5UzZ4HNQbBQgJD0+pS/BA+qEFaKg/cytzlc8LqJMVk5PIVSsvTVZYyfbICbFomnY9LgTJOe+887B7925cf/312LlzJ0455RQ89NBDzTDyCy+80KHcXHfddSiVSrjuuuvw8ssvY+HChTjnnHPw6U9/Oq8fwRjSJpCHoLGr5EhD1pud7E0xa9UUmV2F7GSo45isfXKo7SrFG7EtkiMzrkGn67FsxZROhZWK/VTzGebOGFQ+p5Iy4x8ANMiGaoh4wp/A3OG8SI4ZBSvwGLikOhjCKTn5oRAkBwCuuOIKXHHFFZHf27x5c8f/DwwMYO3atVi7dq2FK7OLLN2OAargsRxpyJq1kc3kZK2a6l27ipjkqDSKhJwFZU/J0bGrJJUcDbtKpSRcO3ysQnK8CWDEPskxUWGlTHIM2VU6nYs9g8Fj2UaAIcJMzquvvopqtYoZM2YYu6YiI/dMjkMnspSPA1Ql5J6U552Z5EgqRFmrppTsqozkRaVXTuAXx65SyaYVUsnRsKtkSYuOkqMycFO7V45KoFczBKwaIjZRYSU7gbx5nCGSo1Mh5RfArpo9cxTDQ0MAgO3btxu7nqLDkZyCIUsjQAAIuEAgc1POQjQEBySspaw2lHR1VWYlR/4pn2VUwIqWyRFB+gTy5msFAAU1R4a4qA7pFIKB8+w/v027SieTozJwU7tXjgrh0JxErmNX6ULVdjJlV+kM2jQ5v0p2OGeIUqnUtKymU4WVIzkFQzN4nKLkAJK5nKxWlESYMWugWNquypzJkVuXMZ7ZU1fL5NDZVbIWlEr4WLZiSkXNkW3upxM8llZyNOwqFcKiZVexoKG8ykK27LwLqoqMESVnQu13QXgehIbyGEInV2M0k6Oo5ADTM5fjSE7BkDV4DEiSnKwKjYS1RKbk1GnsqiyNAEPIZnKEEJlJjsokclnSwhVIjixpUSE5ss39OPfB08aRxEA6k6MTPLZNcmTLx3WPm0SuSo6G7WRCzdGyqzyWWYlNg2rwGAAWTsMKK0dyCoasJeSAZPg4q5IjET6mCh6LOk3wWCZnI2tXNfI42W9ispaVLMlRUXJkSYtKrxyVMQ2qXY9lOxjbrq6qKlhcTagqMpq9cpRJjhElR/3aTeRytCwnYabrsRBCaaRDCEdyHHJH1uoqQDJ8nFWhkVByyKqrstpVnicVsJWxoGTtKlnS4ksOF5XtYqxSRm5DyVHJ2KgQIyEEalarq2wrOaokJ6fgcd5KjgGSo0tSTJAcv16DEOq/p86ucsgdLZKTPluLJJMjoeTQBY8zvl4IqUnkMuqMrF0lm7MhV3IUgsc2MjkqJeEqSk6VC3AJZQ0AmAIxap5PRcnRITmqikxOwWPdEnLBWGYbOwp521WAmYaAOioOACxwJMchb7RKyGelvrZX7Cqf+VKKi0zWRobkyNhVsiXksqRFlhTZCB7bUXIUKuIUwseq+RpVNUeFsGiRHFVlRFNRUSUrTDClQb0h+MSE0kiHjuM1oVNd1TheX8mpa/4cC6fh/CpHcgqGrCXkgISSw1mjGiMLCJQcASF1g8tqVwFyuRwZJYcR21XSSo5sELqgJEelWkqFGKmSFdVcjoqSE3Ah39QzhKqS41cBxfeGC46aZMPQduhYVqqVVc3jNZUcxrhSW4l2eFUTSo7e+xDaVTt27ADXnNXWK3Akp2BoBY8zkJysHzqZG5NEJkdmJpUMyZEpDc8aUgbk1BnOBJjEhidvV8lmcmRJjtzmyRiTvunZyuSoECObSg7nAjXFqeLKao4qYRBC+dhaUIOQtADboRM+1gkdN47XU0BMlICb6JWj2iMnxEHz5qJcKiEIArzyyiva19MLcCSnQBBCtJWQG7SrZDoOEyg5Mq8VnEv1tJCyqyQ3FCnlx5OsTJJVciSf+GWVHxXCoqbkKIziUCE5FpWcWsCUnRTl8LFOgDinrsV6So7eubVJjpHQsAmipKfkDAwMYGxeY6zHdMnlOJJTIOzfv7/5NG00eEyk5MioM1lVHxnSAkjaVZK5FpnwsawyIzukkzp4bIvkKFVXqQSPVZUcheN0qqSUj9Xpd6MYPq7KTj3vgg5JErokR9OuMqHkmCE5esFjYPpVWDmSUyCEVlW5MoDB4ZHU12dXciQ+GBI3MhmSk1XJkR3VIGVXyVZMSVgQfsa5Va1rKVafHHtKjh27SjVbo3KcToBY3a7SIByKKlBPKznVGoRGBsXTDB0DZqqrPAMB6oVj06tXjiM5BUIrdDwbpVIp9fXZMznm7SrZiqnMSo7k5i8TUpYNDsqVnBNXVxFnclQa+9kY66B6jHKVlMJxKqHjEMpDOnOwq3TLwPUyOZqbuxBa4WUTKkzgcXDNrse6dhWAaTe/ypGcAiHrBPIQedpVsuWgmTM5sk3yiErIAVq7ir5PjpBqI69GWOSGdArBwYV8rxMVu6qmuJnYVnImVAiSXwMUR10AULardJUcHZKkW10F6FlWRmZPaXY9Fpxr98kBpl+vHEdyCgSZbsdAvsFjWZKT9fXydpUMyZFUN0iVnOyvl5lA3nGclN2mNsBQRgFirKrU68SmXaWi5OhkcpQIkm73YEUVSFvJydGuaqyhQXIMKDkA4GtYVn69BpnRMXFYdJCzqxxygkz5OEBUQs78TH00pElORqIl29WUqk9O4/XZb2yM0K5S6V4MyKk/qiRH5jiu0O8GmFSAuNzvm2rwWKm6SsuuyoHk5DSaQU/JMUByDqivYYrk6JSR1zXLx0O44LFDbpCZQA40mokFWW7KMkoOkIkUUSk5MkFimdczBTWEUsmRmUSu0tgPkCsjV8nkAHIkRyVbo3qscibHcnVVVSWTozl/Srm6ykAmRybHF0L4vnRWLwo6So6J0DCgZ3vpTB9vx3Qb0ulIToEga1cBGXM5sl1KM7xepkcOIGFXEWVyZPM4QPZMDgsCcK7Q8TbjjVuZ5EgoQHaUHPWbNGPZN3YmBOqKJMcXAp7ksXrBYxUlR3MOU052lRBCaQ0TKo7OOoIL6crMOOhUafmGSE6YyRkfH8e+ffuMrFlkOJJTIDTtKuMkp5eUHBq7SqUle1a7SlbFkT1OtlKqeZzEBmqD5Kg09QshY3WpWlWqxyupMZOo+Vxe3dBVchRsJ1WC0g2V8LIxkqMYPPbrTGtuVsdaOkqOgdAxAMycMQOjMxotSqaDmuNIToHQtKtmyZCcDB8aAiWnMHaV72fqfyHb7RjIrv4wyR45zfUz5nJsKDlW7CrFTA4gR5BUraoQExaVHC6E/PG6mRzOpPvs1FkdXOirGb2o5JjK4wB6wWMT5eMhFk0jy8qRnAKhpeRkCx4DGcPH0kpO+uvJqqsk7SoIkcmyorSrfNlrljzORibHil2lk8mRsLpUK6tCVFn2p3YhBKqaVoa0ZaU6nLNjDbkN34SKA6iFlwuh5BiCTvDYNxQ8BoAFjuQ45AGVTI6XKXjcS0qOvCqSpexcxa7KWnJOb1fRVlcJIawoOTqZHBmCpEtyJiR6/9QDDq5pZUiXkesqOYB0+Fi3R04IJSVHoyqqHUIx02IqdAzoESaTSs50qrByJKdAaJWQS9hVmZQcWZKTTkhkg8dCCPgsfUOUtauAbF2SVZQcFohMeQnZRoCyx6mTnGw3VN+X617dDrk+Ob1hV1UlqvB0Kqua55O1q0woOZJESXduVQglJceQgiECJt2HCzCr5AQeA1ck4aaqqwBHchxygmwJOZAhk8MZwCSfRAiUHCDbaAdpuwqAyKTkKNyohMikAAWqjfSyKjmSE8ibx2UkR6oqDiCbydGorpI4Vjd4LKME6eRxmufLRcnJya7KMXgMqFlWJkkOoGZZCc6VbfEoTKfRDo7kFAit2VUSmZy0cKlsHifjMbJKDpBOjIQQSv0wshAjFSUn63HqSk4xSshV8ziyx+pkcmSOlVFiIo+XITkalVUhpOZXCdEY66CLnOyq3EmOwlpGRjq0r6dAmhpWlZkKL2B69cpxJKdAaM6uGp2V+ZhUJUfWqgIyVV7UFNZNI0bC8wCV8QVZSI5iODTLceqZnGIEj22RHK3qKptKjoTdpRs6bqwhsen5VcBAlZO0XWVIyVGxvUTfKTnyxNikVQUAC6fRaAdHcgoCz/NQnfxFllJy0jZAFZKTRcnh8ht7ml0lO5yzeVwWu0pxNEKW8LEqyfELksnRsauCIADPQApURjO0w2omR4IkKU8R71hDhuQY2vAlcz26Ix2a66goOQY3eBUlx2TwGFBThoyTnEm7ateuXVqf/16AFsnxfR8vvvginn76abz66qumrmlaIlRxAFklh4LkpJCRjCHibqQpOSqhQADgGSqyVPrkANmqsgprV4lsx+ooOVmP57ym1VCN8zpERgWj1zI5UkqOidCxwjqmlJxaUJPqt8PrdQjf3CasUqllWslRWc83WFkFAGNz5qBSqYBzjp07dxpdu2iQJjn79u3Dl7/8ZaxatQpz5szBEUccgeOPPx4LFy7E4Ycfjssuuwz/9V//RXGtfY0wjzM8OhPlSiXzcel2lflMTp3VIRT84TSLS3U+TZaKLOVMDqldlX6c6gTy5vEZNmEbJEensgoAIETmrsfafXKk7CoTwWOJTdyUkiO5jqlMjoBc52STeRxA3q4KPKb1+YuCipJjajhniHK5jAXz5gHof8tKiuRs2LABRxxxBL72ta9h9erVuP/++7Ft2zY888wz2LJlC9auXYsgCPD2t78dv/Vbv4X/+Z//obruvoNKI0AgHyVHpbIKyJDJUbSrqPrkALR2VRYFSHUCuczxunJ1ZiVHE1ktq5qmXVXnHEHGjc1IdZXMGrojHZrr5KPkAHLWl8k8DiBPmnSa98WvqZDJMazkAK0ZVv1eYTUg8+L/+q//wve//338+q//euT3Tz/9dPzhH/4hNm7ciK997Wv4j//4D7zpTW8ycqH9DpXycSBLJse8kqNSWQWkkyNVuyqtgaAQQjmTk62EXFGBmpxEXhkYjH+NIjlrHp9BidJVcrKQJJ25Va01MgyO5RyBgTlDVc4xu5yuqJrok1P3ORgXqJRL6S/WHc4ZgvlA4AEDQ+mnZD4Cbs4yklGFjCs5kuuZtqoANeJkajhnO6ZLrxwpkvN3f/d3mV43PDyMP/7jP1a6oOmKVvm4JMmhsKt40LgJVqI3X1UlJ7WEXKHbceO45M2PB0JZck6zuQTnYBokIfA8WpKTgdxZsas0euSEyNIxWdeqal9n9kA6yakZUHKAhiI0azjD7diUkgM0CFMGkmPKqmquJ6Hk5G1X+RpTw02uaTp4DEyfCivl4PFf/MVfYPv27SavZVpDZQI5APhMgCVt4Cp2VcpxdCRHNXicvG7W8QwqxzZUHHXlIM3qUp1AHoJnUBqsBI91MznIpgbpVlbJrmNCyWmsk3HjM5XJATITJpNWFSCp5Bje3IXnQUj8vpvukQM01GFZZZnCrpouSo4yybnyyiuxcuVKvPjiix1f9zwPjz32mPaFTTeo2lVAipqjTHLiiQNVJkel2zEAiFoayVHf+NKCx6p5nNbxKcRPV8nJYrdZyOTo9MhprpGF5EgM10xCFkWoHrDkBwwJZA4wm6quklgrVyVHcaimqTUp7CpAjjxxxpQrOJMwXRoCapWQr169GqtWreogOq+99hpOP/107QubbmjNrZILHgMpuRwVuwpIJEeqmZzUPjkpZCX2ON+HSBiqqBo6znKsLslJa9WuHTwuSgm5gUxOFrvKmJKTgeSYqKwKkVkRMqnkZMz3mJpb1VxPqrrKvIIhY4FRBI8b62Z/sKCwqgCn5KSiVCrhpptuwnvf+94pREd12N90hsoE8hCJFVY9pOQIjaeVpMos1W7HWY7VfcJKt6toSQ5jDExi6nYUspWQ2wke6/bICZFFyTFRWdU833RScnIMHsuuSabkSKxLYVUBnfOr+nnPlgoeR+Gmm25CqVTCqlWr8L3vfQ9DQ0MolTJUCTh0oDnSwbhdpbgJJzy9UWVyVO2q8Njy6Gjk91QrqxrHpmRyqO0qzbEBaUqQiW6n2UrI9eX2LOFlUyQniyJkKo8DZFSFOFP/PEchr0xOjsFjQNauoukGLEVyDPfICRHaVdVqFa+//jrGxsZIzpM3lJWcduZ344034sILL8SqVavwy1/+0siFTTf0kpKjalcFPEgsRVW1q4A0JUd9M+JMgCVsnPokJ0XJUZxA3jw+5WfXtaqyrmGkuspi8DiTkmPUrsqwmZq0qoDsdlVOwWMhBDjBBl8EJacIdtXw0BDmzGp01+9ny0pZyfn0pz+NmTNnNv//hhtuAACcc845+lc1DaFaQg7Yz+SoKjlAgyANlKN/7bTsqoQeO4FmroV5HJUZ0c8Dqj1ymsfnbFeZIDlZ1CATmZws4WVjJeQZyJJJuyrTWibLx4HsdpVhcuUxDwEPYu8DIUStBhj692xHViWHMa6V50uCTPCYyq4CGpbV+P79ePnll3HiiSeSnSdPKCs511xzTQfJARpE50//9E8xe7b8Rj3dQWJXcdbod6MCIpKTFD7Ws6viyQLTtHySqrN6vbrKFMlJGtIphABTGOjajSxESWaCeOI6RQwem2oEGCInuwrIpuZQWFUy61KUj4eQCTRTNAIMMR3Cx1Ik54UXXkh9zXXXXddUJfr5jTONll2lUF0Vp1To+PcEwWMg3uoSnqf11JbUY0enhBxItrsKHzxOmX1lIpMjhEhch/M6IDGUMX4dHzyl866xTE6GdUxmcjKtZVzJ2Z/pZaaDx0C2iq3cSQ6RVdVYW8KuIlRyFjiS04nf+I3fwB/90R8lDuDcu3cv7rzzTpx44on4x3/8R+0LnC5olZAbVHJU8zgpx6pmcoB4gsQ1FZGkkRBM01agVHJ8YiUHSA4fm1By0tbJUvqdFWlrmSI5dS7AUypOTHU7BoBqHpmcoN5Qe5NewgOtz3scclVyMtpVRVFyaO2qRvi4n+dXSWVynnjiCXz605/G2972NoyMjGDZsmU49NBDMTIygtdeew1PPPEEfvazn+HUU0/FLbfcgne84x1U191X4Jzr2VVxG2EBlZy4SeSqwzlbx8ffiHWVnCRfnjJ4LJjeBPLmOh4DhqJHFNggOSbKx1tr1TAwEP0ZEUKgZqg5n4BAlXHMTBjtYFLJ8ZmAF3AMDSQ8d5osH29fcyRePaawqoCcSU61BsE5SuXkZ3yPqLIKAHjQyPtUBtN1BqrgMeDsqik46KCDsGHDBuzYsQNf/OIX8aY3vQl79uxpTht/73vfi8ceewxbtmxxBEcC+/bta1arKTUDjLWrdJSc6A8W40xrWF+sXaU4nLN5fI/aVSyJ5BgKPeat5GTpb5MVjMVvfFUuIDRGbHQjLd9jMnicaT3TSg6QSpzISE6Gn4WK5ECITGtT2lVANhLFggBMs7ghCdOh67F0dVUQBJgxYwbe85734D3veQ/FNU07hCpOZXAQg0PD0sd7FjM5OioOEB88TgoOZwFPKD/XrZBItKs0SQLnLHYSuTGSkxC8NpHJAdLsKnMkJ2ktU1ZViLTwcSaLSQJVj2HujPhhrcYzOUAqcaIiOVnWJSM5APiBCVQmy6fjQGlXhevPSL4EUqsKcEpOJH7t134Nf/7nf94MyjroQyd0DNjN5OiSnPjgsa5dFX28EILYrtJvzhZnWRkjOQnrWFFyDGZykqwvU5VVzfUSSI7POHxDc7Ka50sjTaarq4BUJcd0+bjMuoKS5Eykv5fkSk4GEkVZWQW0lJzdu3ejrhkZKCqkSc5VV12FL33pS1iyZAk+/OEP49lnn6W4rmkF1QnkIWiqqzwgInipS3Lijte1q3iMXcUCHvlzyCCOJLHAhzCwscaTHDM32aR1TJGcxOoqk3ZVAmEyreQkNRY0mcfJvCaFkpOTXVUEJScN1CTHz9AQkKrbcYi5s2dhaLChHm7fvp30XHlBmuR89KMfxbPPPos77rgD//mf/4ljjz0W73nPe/DII49QXN+0QHMCuUJlFUCUyREciFBddCstYqurNO2quG7JOnOrmmvEkATd0HFrnRjiZ0HJsWFXmQweJxGmmkUlx2RlVYhcMjkp6lA/Bo8ba6crOTJdiVXgZSBR1HZVqVTCgrF5APrXslJqBlipVHD++efjkUcewb//+7+DMYYzzzwTZ555Jr71rW/19bAvCrSUHDW7ygs4eFRVie6cmwiSRKbkJASHsyDO7jLRsTSOKJkiOXGTyHUnkDfXiXkPhBCWSsjtKDmmuh2HSFKGrCs5gafe2DMJKeoQRY8cIN2uEpyDV8393nQjjUAJLow8ICUhS+aHsrIqRL/3ylHueBxi5cqV+Na3voVnnnkGy5Ytw/ve9z4cc8wxJq5t2qCVyUlJoSUgUs3RUXKASJKkreTEEC+dbscAIPwAIkKV0M3jAPFEyUQep7FOPpmcIAiMPZDYqq5K6npsPHicoAyZ7HbcWjNBOaDI4wC52VU+9+EnkDZRrWrbzElI65Xje4z0/EA2pYhayQGARX1eYSVdXbV27Vrs3bs38s/rr7+OiYkJl9ORRMuuUlNygEb4eEZ3LxRdJSeiK2nSWIYsiM/k6BMGUa+jNND5Kx1nNckg3q4ylGcht6uir9+UipO2lq3gsanhnM31EkhT1TdvZSQqORk6BCshza6iOi8aKtHcytzI71FaVVnWp66synoOG0pOv5eRS5Ocm266CSMjI3jf+96HU089FXPnzsWcOXMwZ86c5t/nzo3+xXWIhm7wGLCn5GhXV8XMMNKtrgIaalC5a56aCbuKBQKCC5TKpY6v0ys5poLH8UqOKRTBrrJZQl4lsDISSQ5F6BhIr64isquAhmU1dzgnkpOm5BCHjrOew9lV+pAmOZs2bcL69etx11134f/+3/+Lj3zkI307vdQWWnaVBsmJ2si0Sc7U43XtKo954IKjXOp0SnXtKiC6jNyEXQUhEAQcg11KmbngcT52lUklJ7S+SqXSlO+ZtasS+uQY6nacZb3Ucm+V8yUqOUR2lT9pC0X8uwkhYjuUm0ASgcpbyaEOHQMAZxyBxzAQ040cAHwLdlXYK6dfRztIZ3Le+ta34l/+5V/w4x//GMPDw1i+fDl+67d+C5s2baK4vmkBnZEOISJ75WgHj80rOXFrmLCrouZX6U4gT1qn16urTJKcuBAzY2aGc7bOwxtrRsB4CTnjsZkl092OU9ekUnKEiK3aqgZVox2ku5EUPqYmOWA8cd6dDSUHSO6VE/g+WEAQNu9Cv9tVysHjY489Fn/5l3+J559/Hv/rf/0vvPe978Wb3/xmfOMb3wBjdn5B+gVNJUexhBygsqvMKzlxa5iwqyKVHEMVSlG5nMA3Y1fFDek0RnK4iKzUMmlXxa1ncjhn2pqmMzkCIlbNoQgeMy7iS9OplBwg1rKitKqA5FAzOclBsmVli+QkTSOn7pETIlRytm/f3peV0drVVQsXLsSf/dmf4amnnsK5556LD33oQzjqqKOk17n99ttxxBFHYGRkBMuXL8ejjz6a+PrXX38dl19+OQ455BAMDw/jmGOOwYMPPqj6Y+QK3RJyIELJ4Vy/5JRIyemWwAVjEAaCnJEkx9BmFFVO2it2VdxaJpWcuPVM9shJWpMJgbphkgPE53IolBwgIZdDpeQAsSSHqrIqRLJdRW/TJBEpG8FjIFnJsWFVAa1Mjud52LNnj5Vz2oR0Juf//J//E1lZ5ft+kwXKjny49957sWbNGmzcuBHLly/HbbfdhrPPPhtPP/00Fi1aNOX1nufhbW97GxYtWoRvfvObOOyww/DLX/4S8+bNk/1xCoFmdZXJTI4JL52gTw4wVcnR7XYcImp+FTOm5NglOaYmkDfX8zkwo/NrNkiOydBx0pqme+SEiFOHKPrkNNYNMH/m0NRvEI1XSFqbaqRDlvWnjZKTQHJshI4BYHBgAGNz5+C1veN4+eWXsXDhQivntQVpkjM6OopDDz0U8+bNS/wjgw0bNuCyyy7DJZdcAgDYuHEjHnjgAdx11124+uqrp7z+rrvuwquvvoof/vCHGJxsSX3EEUcknqNer3fM5hgfH5e6RkqQ2FUFJjnda3BDZCGqoaCphl6RdpWh6qqoSeQmVZzGelOv346SY57kRCk5pvM4IaLIE+MifiiuJmLtqpQqKC0UUskh/Hmb50hQcqzZVfmTHKBhWYUk55RTTrF2XhuQJjl//dd/bfQCPM/DY489hmuuuab5tXK5jNWrV2PLli2Rx/zzP/8zVqxYgcsvvxz/9E//hIULF+KCCy7Axz/+cVQq0Un1devW4YYbbjB67aZAEjzWDR0DZJmcbpITN1xTFlEVWiZKyAHa4HHUJHLzJIc+kxOt5BDYVRFrms7jNNeNIDkUlVWtteMyOYQbXk6ZnPyVnCSSQ19dBSRXcdloBBhi4dh8PPPcL/uywko7k6OLPXv2gDGGxYsXd3x98eLF2LlzZ+Qxzz77LL75zW+CMYYHH3wQn/zkJ7F+/Xp86lOfij3PNddc02Gvvfjii0Z/DlXUarWmwjRDqxlg183YgOLSTXJ85oMbqJSZQnIM2VVRFVommgE21qEjOVFr2SA5vZrJiSojp7KrotalyuMASSTHvl1FreTErS8Yg9CcZZcFcXZV4DFwwxPm41AUJaefe+VIKzlFAOccixYtwh133IFKpYJly5bh5Zdfxq233oq1a9dGHjM8PIzh4WHLV5qOUMUplUoYnjEz5dXxmJLJ8U3YVZ2kwYRVBUztmqw7nDNEd4WWEAIsMHOz6laEBOdGyzsDz8PwaOvf31QjwKT17GRyzPzOtCPSrqJSciLWpaisSlzbrwGc0D7Jya5igqEW1DAyMNLxdT4xQT5SoXmeCCSFgU0jMZNTs1NdBQCLDurfMvLcSc6CBQtQqVSwa9eujq/v2rULBx98cOQxhxxyCAYHBzusqeOPPx47d+6E53kYGooI7hUUYR5neHQWymV1YW2qXWU+k2OK5EwJHmsO5wzR3fci8Lmxm2W3ImRSxWms10XQ+sSuoqmumnrzt5nJsa7kEAeAY+0q6vOiQaSmkhw7CkackmMrjwMk21W+IYU7Cxb2sZKTu101NDSEZcuWdTQT5Jxj06ZNWLFiReQxZ555Jn7+85+Dtz1lPfPMMzjkkEN6iuAAbeXjGlYVEBU8NmFXeY1S9EmYyOMAdJmcbrvKVB4HmGpXxfW2UV6/264yHGztJjmcczskh6RPztSbv+lux611ozI5lCQn4t+EMnQM5GZXAdG5Hxuh48Z5on9umyRHcNEYBhoBu3ZV/yo5uZMcAFizZg3uvPNO3H333XjyySfxwQ9+EAcOHGhWW1100UUdweQPfvCDePXVV/HhD38YzzzzDB544AF85jOfweWXX57Xj6CMVvm4+gRygKi6qmsdU0rOlOoqIrvKVGVV1FqmlRy/blfJMW1Vxa2ZNIZBFXlXV9Uo7aoolYhcycmR5ET8bDZCxwAgPA8i4nfWtzDSofN8EVayVwdn9q6jn5UcZbvqoosuwlvf+lb8f//f/4ejjz5a6yLOO+887N69G9dffz127tyJU045BQ899FAzjPzCCy90WDlLlizBv/7rv+Kqq67CySefjMMOOwwf/vCH8fGPf1zrOvJAa26VnpLjBRycC5TDIZImlBygQXKGRgHQkRxTdpUIGITvozTZVoBFjbpQRK/bVbzr+ilITpQyRKHkMJt9ciKrqyxncqiVHB40cj+DLduoFtSMFBmkIUrJEZZIDtCwrCpdLU9sKjlAw7IandPpQPgWVRygRXJee+01TExMYHR01Or5KaFMcoaGhrBu3TpceumlOOyww7Bq1SqcddZZWLVqFd70pjdJr3fFFVfgiiuuiPze5s2bp3xtxYoV+M///E/p8xQNrW7HekoO0FBzZoTD3owpOa3N1xjJCWjsKqChClUmSY5JJafb+jLVI6e1Xnd1lengcdf1G7aqADQbgrYP6aRQcjivQwiOUtuQV6vBY8JMTtVnUwedWsjGwDvQQXJsqDhAvkpOeK68SU6UkmPTqgKAWaOjGBkeRq1ex8svv6y0hxcVynbVV77yFTzzzDN48cUXccstt2DWrFlYv349jjvuOLzhDW8weY19jVYjQD0lB+gKHxvKz/SSXQV0qkJGJpCH63LR0T2ZGVZCbJeQUyg5QogO8hSSEYITTcnlUNlVTAjUutamrK4SIoJEUY50CNE1G8sayYnM5NglOd2wWV0FRPfksdkjB2hU9/arZaWdyRkbG8NBBx2EsbExzJs3DwMDA33XFpoSJhoBhujI5ZhqHtZGckwFjz3mdQyCM2VXAZ09d0wGj4FO0mReyelSt0x31A06x0RQkJzudSm6HbfW7vz9plJyAGCia21KJQeIsMMoh3OG8PIhOVHnsUpyIiqs7NtVUUqOvfcgRL9OI1cmOZ/4xCdwxhln4KCDDsLVV1+NWq2Gq6++Gjt37sTWrVtNXmNfo5XJMUBy2jf1AgePBQQ83iJMpu2qEKYaATbXa9t8zGdyupQcg1Yb0FAI2omTFZJDkMdprt1GoOqcgxH2VelWiSgzOUCEUkTZ7ThEl1pko3w87jxJnYhNI4pQ2ep23DxfFMmxrOQA/Rs+Vs7k3HzzzVi4cCHWrl2Lc889F8ccc4zJ65o2MFVCDlDZVeYzOeFaw5VGc0Yqu8q0ksM6lBzD1VXEwePmmpOZLYpMDtBJcjhBj5wQjLc2JyqrKkR7qJlzMbUnlenzdZMc6uAxkJtdlbuS03Uuxrjx+0Yaonrl2M7kAI7kTMHWrVvxve99D5s3b8b69esxNDTUDB+fddZZjvRkhIkJ5CE67KoCKzlAy/oSQkAYJAztqpDJTA7QGWSmDB6bnkDeXNe3reSY73Ycop1AUVpV3es3gsGkp5vaK8cGyek6B/XcqvbztAethe9HlnVToduuSupATIUoe8xmI8AQYa+cfptfpUxyli5diqVLl+JDH/oQAODHP/4xPve5z+Hyyy8H5xyM2f9l6UUYtavCJ0zOAWboRtGm5JjK5ACNElVgkpQY3DV4W0NA85mc1u+0T1hCTqHiNNZtu34LJIdUyWmzq6rEc4balRzqPA7QZVdxbu6BJQk5ZXKEEKgGVYwONkqWbao4UeeznccJz9ldUeeUHHNQJjlCCGzduhWbN2/G5s2b8fDDD2N8fBwnn3wyVq1aZfIa+xpGg8fh5sgMEgdqJcfwE0t7Q0DjSg6hXRXOwqoMDBKSnLbrJ7Kr2telGOkQon26ObmS005yiPM4QJddFVStzHHqLlO3lckBGmpObiSnAEqO4AJ+nWFopLUd25xbFaJf51cpk5z58+dj//79WLp0KVatWoXLLrsMK1euxLyungMOyWiVkBu0q0w++RFmcgCzeRygc35VYHhDYoR2VWNNzxrJsaLkRDTtM4V2JYeqEWBzfW5XyZloP4eN8nEgNyUHAKp+FZjR+Lt1klOtQXCO0mSz2TyUnPC8Icnx641rso1QydmxYwcYYx2zIXsZyiTnb/7mb7By5UrMmaMfmJ3OMFtCPvkBNdXtGGhWdggh4JuywNAiOd2jGHTRnslhhsuwA8I+OUBrErnpRoAhbNtVUZ2JTaG9cos6eNy+PnVlFQBU2zM5NsrHgVxJTnv+xzbJgRDgE1VUZs0EAHiWK6tC+DUGzG38PQ+rCgDG5s5FuVQCYwyvvPIKDjnkkFyuwzSUSc473/lOvP7661i/fj2efPJJAMAJJ5yASy+9FHPnzjV2gf0MxhjGx8cB6I91AGiVnDqrQ8CcbN4kOabtqkmSI7gwn8mZ3OAC3ydpcheqQ2RKjkdvV9nK5LSvXSN+6p3I066ypeQwv/GnMgif+fC5vfBvJ8mxv8HzAweaJCcPuwrorLDKo3wcAAYqFcyfNw97XnsNL730Ut+QHOU+OT/60Y9w9NFH43Of+xxeffVVvPrqq/jc5z6Ho48+Go8//rjJa+xbhAQHMNwnx6SSM0mYTFpV7euZt6sa6wamm+mhFWSmsKoa607mlAiuvX3dIAjAiYiBtWaA3J5d1VldRf+k30FyLGZj4O1vnN9SZVWIalsfIFsTyNvRfs7c7Ko2cpWXkgP0Z/hYmeRcddVVeNe73oXnn38e9913H+677z4899xz+J3f+R1ceeWVBi+xfxFaVQNDwxgYGkp5dTrqjE7JMVlZBbQpOQa7Hbevxww30wNawWPToeMQ4SRy6kwOlVXVvTZlM8B2Jae7I7FpdCo59FkJL+AIwnPaKB9vnrhBbmxaVUDOdlXXOfMiOV7beW0P52xHP5IcZbvqRz/6Ee68804MDLSWGBgYwMc+9jGcdtppRi6u32FybhUA1H2CTA4PABaQKTkmux0DAJ9cz3S3Y6DNriIiOeR21eR7QklyOmZXESo5nPvg3Ee5PEieyQmEgMc5hsrlqT1siDDhM8yplO0qOZPnsq3k5E5y2iqsohrz2YBfALsKABb2YYWVspIzZ84cvPDCC1O+/uKLL2L2bH3rZTqg2e14VH8COUCUyZlcz7iSE4R2lWHrh3EIzzNePg60KznEdhWxkkOVxwEAzjl83wfnHoSgfSoOq7eoSQ7QUnNqFqqrAKAWWlZWlZyGXWVdyWkjciJHJUdw0dHw0ybaFaQ85laF6EclR5nknHfeebj00ktx77334sUXX8SLL76Ie+65B+9///tx/vnnm7zGvkWz27EhJcdnHJwLs0oOAAT1nlFygAZxomjNzoPG+0un5IQkh6i6yqO3q4AGiaLM44RgrAohBGoE3aG7EeZybFRXdZzHaiZn0q6yMSurDR1KTg5WTajk+B6z05MoAu1DOnNVcvpwSKeyXfXZz34WpVIJF110UfPJcHBwEB/84Adx8803G7vAfkZTyTEQOgYan0+PcYwYJzlV4ySn2QyQgOSIeh2Bp59xigLzOQK/R+0qLiAYJyc5vu+jXKa/UTNWxQTnRqv+4jDBeINQWZpr1CQ5tqqrgKZqZNuuqgd1cMEBz4ewEOzuRqjk5FVZBTSscMEFSuVSIYLH/TTaQZnkDA0N4fOf/zzWrVuHX/ziFwCAo48+GqOjo8Yurt9hshFgiLrPMWLcrjKv5NQmn/TbxzCYAq/VwQKaRlaBzywoOXQbqfDtkJyBAXolh/MaasQjHUJUGUfN5+CWnvSbVVy2+uS0ncu2XSXQGO0wPEHzuUpDU8nJKXQMTHY99hgGhyu5zK0KsaAP7SplkhNidHQUJ510kolrmXZoDec011CxHjCSTI5pksMFh899GiXHqyPwRoyvCzSGdPY6yaHM5AANkjM8bMeuoh7pEGKCcSvdjpvn8xjAAiCwuPHnVF0FNHI5QzmRHDGpnORJcoCGklRCQNKDKytCu2r//v0YHx/vi2a/UiRnzZo1mV+7YcMG6YuZbmjZVWaCx8Bk+LgHgsdAw7Iy3fEYaDQYDHya8DvzOWHwuE42gTyELSWHcm5VCMarVkLHAFDlwlplFTBJcmzmcYCWXWX7vGhYZHPyIjkBA6/VcqusCuHVAghBc2/JitGREcwaHcX+iQm8/PLL04/kbN26teP/H3/8cQRBgGOPPRYA8Mwzz6BSqWDZsmXmrrCPYbqEHAiVnOIHj4FGwJEseEzUUI/SrvK9OqmKAzRCzVZIDmGPnBDcqpLDULOYSa16zG5lFZCbXQU0iBXPieQADcvKr+cTOg7h1xmEhc9NGhbMH2uSnOOPPz7vy9GGFMn57ne/2/z7hg0bMHv2bNx9990YG2v4eK+99houueQSrFy50uxV9ilMzq0KUfOJlBxOoOTUDmCQQLUQ9ToCTiM9U9pVgnMEVdonOVt2FWWPnBCMWVZybNtVvuVsRlAHCzwS1TYNE0HOJGdiAn6NplghK7wag7CggKZh4fwxPP/Sy32Ty1EuIV+/fj3WrVvXJDgAMDY2hk996lNYv369kYvrdzSVHIMkp+6zxgwak/DNZ3IAoFbdZ3xNoDGJnKKEHGjYVYyougoAgglikuNZsqsIh3OGsGlXTTBuZW5V83xeYN+uEgLV6q+sVKt1o6Hk5Ncfhk9MFCCTE+RaPh6i3yqslIPH4+Pj2L1795Sv7969G/v20Wxe/YZmJsegXeV5VfO9HoIa6iXzm683sd/4mgAg6h6CCpVdRafkAKBXcgI7JIdyOGcIzmrW7Koq45jw7W3+NZ/ZLR+fRLX6qvVzAg2LjBMT/CTwAwfg1Wbkdn6goeRwy+X7Uei3XjnKSs673/1uXHLJJbjvvvvw0ksv4aWXXsI//uM/4tJLL8W5555r8hr7Fq3qKnPB46BO8CEhyuR4VSqSUyPpeAwAft0HCwjHIhCTHO4F5HZVEATWlBzq4ZwhJhhvdSG2AJ8JeDX7D4sTOZGchl2Vn4rBD0w0x7bkBb/Ocu2RE6Lfuh4rKzkbN27ERz7yEVxwwQXNJ8OBgQFceumluPXWW41dYD+jZVeZU3ICArmT+RMIBs1vjF6VJlgZVGvgJZrNrz5Bu3kHVQ8VDJOt79Xocw+2lBxmUcnxBMf+ut3qG7+6H7ZTItXaa5bP2MCEPwF+ID8lxxs/0OgWnyO8WoBKqQgkp7+UHGWSMzo6ii996Uu49dZbO5oBzpw509jF9TOEECTB48Az/yGp+xPAoPnmej4ZyfEBop6UHkHzwnYEVQ/DhCTHt0RybCg5EBwThKpaN/ZaLCEHGkqn7btptT5u+YwNTPgHwKv5BY+9fVVYf7O7EHjcSlViGpyS04WZM2fi5JNPNnEt0wrVarWpgJnseExBcjx/P4C5xtelIjmM0PLxqrSbN695GiZyOvyaB9A0g26CMQ8s8FAuE/4gkzgQ1IESHSlsx3gtQMnKmRrwazR2bhIm6q8Dg/T/bt0IqhPgrIRyyf65AcDbX8ud5AjOUD8wgcqAzd+yqQiVnF27dsH3fQwODuZ6PbrQIjmbNm3Cpk2b8Morr4B3ycZ33XWX1oX1O0IVp1QuY3iGuU8X82uA4d/Jul8DBckJqhOg2HH9Gt3TvV+ro0x4D2L1ACDMP3qeR7r+5FnAOCMnOYEA6tzHwIAdkrPPCzCnQswQ28BqB0Ao6kWiWt8HDJr/rKehUvPh8xKGK5Z/4El4+2vA4lxO3QTz62A+Q2VAW3vQwrw5szFQqSBgDDt27MCv/dqv5Xo9ulC+C91www14+9vfjk2bNmHPnj147bXXOv44JKMVOp6NUsncrskp7CruA9y8VB/UaaRZxkoQjCZE6NfrpIOKOXHuI/Bs2DseuIVAcF1UIIQdCykIOOqWQs7Nc9YtNwMEMOHnUxnbIDn2rMduNBrx5TdOAQBYUAcL8s0FAUC5XMZBY/MA9IdlpRU8/vrXv44LL7zQ5PVMG5ieQB6C+eatmrpgDZJTNvuEwWpVAOYqy5rrihJEEKBE8NTNfA9iUKBUoZFzGDEJ8VkAcAFSOQo+GFEzxnbURBnC0sYYBBy+peGcAFDm/mQ/JrtWQdWzb5EBLSUnLwS8DOH7KA3noyQBk0oOUad2WSw6aD527flVX5AcZSXH8zycccYZJq9lWoGiESAAIKgZb+blCQYQbFqMqPEV42WASMnhzCetwhDEZaw+ZxDkk7vtKDk1XgG3pOT4AUdgsfpmkFfhW1aOAKCaw9wqoEFygjyVHFaGIOx/lQUsqIPl8G8ehYVj/RM+ViY573//+/G3f/u3Jq9lWqFpVxkMHQNAhXlghm/GTSXHMHidJsQbKjkU4MwHJyQJnJjkBDwgHQDaQH1KRo8CNVG2alfZVHIGWRU+ORnthBCiMbfK4s8ZolL34Zvu1C6BUMnJE8yvgRfArgKABX1URq7sP9RqNdxxxx34zne+g5NPPnlKAttNIU9G064aNUxyuIeACQwYzHzWwUlIjqh74DO48YoKxukyOaQkRzQmIlMiYAHAOGhLrDwrdlVdVOzZVX4eJMfuU31VMAjBAe4DFbsdehp2VT6VVUBRSE5x7Kp+KiNXJjn//d//jVNOOQUA8NOf/rTjeyaDtP0KignkADDA6+aVHB4AzDzJKXsB2AgzT3JEmU7JCXwylaLES+BE5CxEYMWu8q3YVVVh167yrdpVNQTWSc7ke8mmJ8kp501yglrhSE4/zK9SJjntE8kd5NFqBGg2eDvAPePZASolp+wzMM4wWDYbruSiBBCRHMZ8MpJQ4iUIISA4R4mo/DrgAUB+H7Wj5NR4/9tVAgIlS915Jpokx4PtpjF5khwuGg9GJSuVh/Fgfr0wdlU/dT3OjzpPc1ApORVRBzOsNHgEmZxSwFBiHAHBJhXwck/aVaXJ6hJOoJqFsKPk2Coht1tdFVgkOQO8CgEYV2WTUOVtSo5l5JnJCSbJVRHsKs64hcxcOhYe1CI5IoeMlklIKTlr1qzBTTfdhJkzZ2LNmjWJr3WZnGRQzK0CQiXHrDJSJ6iuKvuN9Sie+Dl18JjoibNFchgqRJXDAQ9QMvz7MRWWlBybdpXPEAh7ysoga4TyfcP5uiR02FU2wQUqdR9BTs/cPsuf5AjOwVijuosFAgND+UY+Fkz2yanVanjttdcwf1LZ6UVIkZytW7c2RxFs3bo19nUuk5MOErtKCFS4j4Cb7WdDUV0VkhwKJYfxEkAQ4BVCTCo5NHmFdpJDAcYZuBAo94mS06iuYo1qIOJ7TjBpI/hcYIi0x1ADg5MzjHzGMYNgblwUOu0qe6h4PiAABg4mGCole12lgXYlJ78SchbUm1VtjHEM5GyyDA8NYe7sWdi7bz9efvnl6UNy2nM4LpOjBwq7qiI8AMK4xE1hV5Unhx1SPPEzUUaJwPIRPACEIOuTUxaNGxsVyQnC95pcDvcslZBXAAgIEaBUolWn/CAk5cLKZPCWkmMviNq0qyz3q6m0jWHxuY+KxdEZABDwxvmEH0AIkctDOvNb7TSKEz6ej7379uOll17CSSedlPflKEOZLl500UX42te+hmeffdbk9UwbUEwgH+CNJxHjwWNKkmNYyeEitKsIyNPkxGv6TA7RSIrJf0PKTI4QDAIBuODggvZmXZ98AueCflMOJjceW+HjQd5oymezV85E+D5atqs6SE4OuZygzX4Wvt1J8yFY0OpUX5Tw8YI+KSNXJjlDQ0NYt24d3vjGN2LJkiX4gz/4A3zlK1/B//zP/5i8vr4FRcfjCm98UJjBG6MvGhYHmV1lWMlh4Q2LIJPDJ2/Agujpmjp4bEfJaUn+lJZVXZTAJrMxgqDyrx2M8aZ6Z6WMXAgM5KHkiMnfD9t2VZeSYxudJCcfy6p9HE9hlJw+6XqsTHK+8pWv4JlnnsGLL76IW265BbNmzcL69etx3HHH4Q1veIPJa+xL0Co55j4k9fBp3PBGUmnaVWbXZWJy4yMgCiHJobKrqJWcIFRySK2k1iZBGT6u85alQR0+Dto2HRtKToXXUZoczWKzV05e1VXtJCeP0Q6dJCenCq8C2lWLDuqPMnLtdNPY2BgOOuggjI2NYd68eRgYGMDChQtNXFvfIggC7N/fGIQ3w+BYh8okyTGZyamHT3dkwWPTSs4kyQmY4QlebSSnR+2qICR+HDD+5jRhS8lp25iI7SrfMskZ5K0Nz6ZdlVd1VaWer5LT3p9H5NQrx9lVdFAmOZ/4xCdwxhln4KCDDsLVV1+NWq2Gq6++Gjt37kysvHJoqTiAWbtqYNKuMpnJ8UISIrjRmTZlMiWn7VfaMFlo2lVckIz3oVdyWuvSqTl2lJyqaCk5grhcvUPJsWBXDbLW4FpbdlVdMLBQtRWcpMN5HCrVdpJjPxNTDLuqTckpypDOPmkIqFxrfPPNN2PhwoVYu3Ytzj33XBxzzDEmr6uvEZKcoeEZqAyYqwppVFcRKTlhLsdQA5fyZMDPtJLDRasyQgQBSgYrNXjbEy5nHBXDDUxKwo5dBQBggmh8lR0lp2ZTyfFbP4eNhoBh+ThgT8mpdpML7gMVs60o4pBn8FiIlvoL5GdXFTKT0yejHZTv0lu3bsW1116LRx99FGeeeSYOO+wwXHDBBbjjjjvwzDPPmLzGvgPVBPIBSrsKMGpZtZoBGlZy2p/KDIePedC+gZvffKwqOWSbZ2uTcJkcNeSh5Ex0v4cWw8d52lUBL0OgACQnaCk5gtG1qZBBSHJ+9atfoVarpby6uFAmOUuXLsWHPvQh3Hfffdi9ezcefPBBDA0N4fLLL8fxxx9v8hr7Ds0J5AatKqCtuoo3OrOagNdBcsxtWqFdZTyT06bkUNlVAEH4WLQrOTSbtt++LtlNNAclh3hjtG5XtWVyGBdW2upXp5Ace5t9ntVVQVf38jxIjuAMPOg8bxFyOXNmzcLQYEO53759e85Xow5lPVIIga1bt2Lz5s3YvHkzHn74YYyPj+Pkk0/GqlWrTF5j34FuAnljgwln3gwY6Mxa6yA55m4AZa+xLufMaKt8SiWHddhVZm9CpTbJvLeVnLYAJWEVV03YU3KsB4/blByByS7LFdoGdVPsKltKDufNBx7AfnVVN8lBDiSnPXTc+hrHwFC+XY9LpRIWzh/Dy7tewcsvv4yjjjoq1+tRhTLJmT9/Pvbv34+lS5di1apVuOyyy7By5UrMmzfP4OX1J5rl46NmJ5CH1VVAI3xsIjJCpuRMZnIahIxhoGzG/2cdmRxCJYeQ5DQmkTOUymZDMx2ZHAtKDqVd1ZnJIbarOjI5QOO3lo50tGdygIZlNVSh3fCm2lV2Nvt2FQcAuBAIeGDsfpCG7snneSg57eXjIYqTy5nfJDm9CuXfpL/5m7/BypUrMWeOWTViOoBOyWkLrxnaxOrtXWsN5mcq7U9vIsCA+q9iB9pDhKYrRNolZdOTgttJDjA5pJOQ5NjI5NCWkLdVVxEHj4OODUfAFwKDhK3/w0aAzfNbCB/nZVd1kxygYVnZIjlT7CouGgULA3bOD3SGjkNQtamQxcI+KCNXfjx45zvf6QiOIlokx3Amp+1mb+rG6FEEj7lAqW3jMPnE315Cbjx4zOg28CiSYxJCiM7u0j2u5FTbbUnBIQhHSPhdT9U+8UN2lJJDjanVVXbsqjiSYwtT7CrYV3Pi7KoiYEEfVFjla/pNU1B0OwaAga7AoglQVFeVu+bDmJxE3lEO2qN2VWN9s9fePT6DajSFveBxp8pFaVkF3SSHeCZXeyYHsENycrOr6o7kRCk5RSE5/dArx5GcHNCaW0UTPAbMjXYgITle5zpUSo7J6irBWUfTOdPVVWXR+VE0T3K6n9TNKzlCcIh2u4poSCcXgCe6bQaajYlzMaU5G2nhi+Adn2PATq+cItlVNsPHkSTHs9sQkEVmcpxdZQqO5OSAlpJjMHgsBCptNwdTSg6FXRX2yAlhchJ5dzNAU2BdN/3eU3I63wuaTM7UzYFCzamL8pQGCVQVVt0qDkBbRt5QcTrXt6Lk5FRdFWlXWSxfL4aSM5Xk8IIoOf0wv8qRnBxAETxudDtu3RxNjXaoEVRXVYiUHC7oSA63TnLMbtpTpr2TbNRTN0aKXE63VQXQ2VV+RACHsoy8vUdO83zEJMcXHH53vyrOjFZTxsFlcnojk7N9+3bSlhCUkIqQr1mzJvNrN2zYIH0x0wWUE8hDmFByhBDwYb66qlvJMdUQkHXfsAyqIVOadVmorjIJn/WPklMTERsT0cyjSCWHlORUp3yNurpqilUVgvmA4Qq/buRJcgJe6ngoCmF7SGeUkiO4AGcCZeL+SGlYMG8eSqUSfN/Hnj17sGjRolyvRwVSJCfr4M2SYnnl7bffjltvvRU7d+7E0qVL8YUvfAGnn3566nH33HMPzj//fPzu7/4u7r//fqVz20Qrk2NyAnnn04AJJccD7+y2SpbJMbMu67phUSo5pkmC7eBx4yQCMNAwsgU7Sk49QsnhRGXkfkR43eQA3G50l48D9ErOFKsqBPcBjJCeOzp4bGdIZ5SKA9hVcjhnsaotYxxlg7P3VDAwMICxOXPw6t69eOmll/qf5Hz3u9+lug7ce++9WLNmDTZu3Ijly5fjtttuw9lnn42nn3468Y19/vnn8ZGPfAQrV64kuzbToLCrKJScerfCYshC6SY5kRuwArqVHMGYsbZtU+yqyUnkptqllITdTA7QIGolYpJDIXFXoywGskzO1M+RbSWHepREvJJDn8sZqE49h61MThFITpSKE4IHAhiydimxWDh/DK/u3YuXX34Zp556at6XI43CZHI2bNiAyy67DJdccglOOOEEbNy4EaOjo7jrrrtij2GM4b3vfS9uuOGGnmk5LYQgsasqXSQnMPD0503x6YsdPO5WcgAAhtScbpIDIYw2BMxNyTEKW3ZVRCaHyOLw/anvGynJYVNJDuMCjPCcU8rHmyem3exLAevolxUiEL6x2XtJKAbJmZrHaX6vYLmcXg0fa7d1fOKJJ/DCCy/A6yq7e9e73pV5Dc/z8Nhjj+Gaa65pfq1cLmP16tXYsmVL7HE33ngjFi1ahEsvvRT/8R//kXiOer2Oer31CzU+Pp75+kziwIEDYJMbmEm7aqDLrjKh5NS6SY7gAOdAWY8bT8nkGFNyIvx1Q91Lp5AcwKhnbru6CphUcoyexZZdNfX3r2+qqyKCx0DDsqoM0FgXUxoBhiAmOVF5HAAQovH7OlgeJD1/LMlhDIJzlDTvc1nQPn186veKQXJ6vcJK+e7/7LPP4t3vfjd+8pOfoFQqNbMbYR6HSdyk9+zZA8YYFi9e3PH1xYsX46mnnoo85uGHH8ZXv/pVbNu2LdM51q1bhxtuuCHzNVEhVHHK5QqGRmYYW7dRXdWCkUxOVCCYB0BZT0OdkskxpuRE2BiGyEJ38BgwGD4WUXaV2U07KudgejRFrkqORZITCAGq+VVRSg4wGT4mmjSQl10VR3KARvg4L5IDNHrllEZo80hAmpJTlF45vU1ylKnqhz/8YRx55JF45ZVXMDo6ip/97Gf4/ve/j9NOOw2bN282eIlTsW/fPlx44YW48847sWDBgkzHXHPNNdi7d2/zz4svvkh6jXFoH+mgGtCOgpVMDmDEsprS8ZhQyTFlV7GIG76pDbxbxWmsbVjJiVrPOAGxVEIemcmhCh5PfY8EBFlDwDiSQxk+zsuuigodh7CRy0kkOZYsq+RMTjGUnGlrV23ZsgX//u//jgULFqBcLqNcLuMtb3kL1q1bhw996EOZK7EAYMGCBahUKti1a1fH13ft2oWDDz54yut/8Ytf4Pnnn8c555zT/FoYcBwYGMDTTz+No48+uuOY4eFhDA8Py/yIJKAa6dBdXcW4gIBASeNpk4rkVDyaTA6PUnKoMjkw1ysniuSYnkQeaVf1kZLDeUAirkQpOUAjlzNAMKQzKngM0JKceLsqXyWHGj4rAMmJ6JHT+l4xSM7CHp9fpazkMMYwe3Zjo16wYAG2b98OADj88MPx9NNPS601NDSEZcuWYdOmTc2vcc6xadMmrFixYsrrjzvuOPzkJz/Btm3bmn/e9a534a1vfSu2bduGJUuWqP5Y5LAx0gFo3O911RwvqiW/CSWny64SAmAGeuVEBY+N2VVRJMcQSYgiOY1zmlNBItUy4/1X8svkAILEsoolOUS5nIGYTA5lr5xYu4qYaCSRHBujHYqh5Di7ihrKSs6JJ56IH//4xzjyyCOxfPly3HLLLRgaGsIdd9yhVOm0Zs0aXHzxxTjttNNw+umn47bbbsOBAwdwySWXAAAuuugiHHbYYVi3bh1GRkZw4okndhw/b9685nUVGWQTyCOmBgdcYEAjO1eLuvkZsauieo8EqGj2hKC0q6JIjqleOUkkp2IglsB49AypflJygEYup2QwuCKEiCwhB8JcjlmURNAxmqUdfWlXJSo59L1yikBygiS7ysI4jywIlZy9e/fiwIEDmDlzZs5XJAflO8J1112HAwcOAGhUOf3O7/wOVq5ciYMOOgj33HOP9HrnnXcedu/ejeuvvx47d+7EKaecgoceeqgZRn7hhRdQtpB2pwbZBHIR9RRdQCVHiCmZHMCUkkMYPKZUcqJK32FOyYmyqgAYVXK6h3OG4KLRUNJU/iwQgB/3fgkfZYPN6xp5nOj3iKKMfDCiEWDzfERKDhciusAAaHzWBQdKNPfdPO0qLqLvFyGKYFcVpevxzBkzMGNkBNVaDS+//DKOOeaYXK9HFsok5+yzz27+/Y1vfCOeeuopvPrqqxgbG1O+oV1xxRW44oorIr+XFmb++te/rnRO26BScgbY1A+LrsQdncnR23jLPovcN0zYGtEl5PrrCiEiq51MZXLKMU+U5khO9DpmlZz4TYFxhoGKGYUlqttxCNOjHYKIuVUhKOyquNBx43w0T/VVEXR2NW+HEA01Z4Amy5gYPCYmOUkqDmBntANnPkRKFSUL8u96XCqVsHD+GF7YvqP/Sc6aNWtw0003YebMmalzrNzsqmg0lZxRw3ZVlJKj+bQZW0KugSirCgACA3mKyCczA3YVZ37jhj/l672RybGh5ERZVSE444Ch+3TU3KoQpjM5cXkcgEjJiQkdA3R2VaxVFYKS5CQpOcRWWVLoGLCj5CTlcZqvCTgGh/MlOQCaJKcXw8fSs6v8yX/8pOopk6XR/QaKkQ5AdCaHkSg5miTHiz6eTMkxQBSirCrA3MiCeJJjZtO2o+TEkxyT4eMqj7/hm24IGFU+HoIik5Ok5FAFj2Mrq0IQVljlaVelKjk+/UiLJKuq9RoXPtaF8uwqyjlW/YyQ5IzMnGVuUSFigsd6mzAJyYnI4wAJakNGMA5E9e81UUIeS3J6RMnpnkDeOoEwWHYdf8M2GaBMtKtsKjkUdlVMZRWQo5JDRDbKXoBSws/EeKDdAiMJQQJZBgDhB0azZFFI6pEToijh417uldP7Sd4eQ2hXmVRyKsJHKSLoUsTgcdmL3rh1g8exIUIDakivk5wkAmlOzbGj5ESXjzdgen6Vn5TJIVByBhKUHC70H1qiEFs+HoLINkpScYAG96ZUc9KUHKBBdCgRZLKriqHkLJouSk5aDqcdLpMTDYo+OVEqDqA32oEJDr+nlJxoomAieBw10gEwRxDoq6sS1mEcMBJsjN+QTD6NVhOUHNN2lX0lJ57kAA3LSqclRBTysqvSSA7QyOUMaY6QiUMQc79oh/A9YIhutETS3KrWa4qh5CzsYSVHOpPTjscffxxBEODYY48FADzzzDOoVCpYtmyZuSvsM7RKyM3ZVd3DOUPoKDn1KBUH0K+uIlJyorodA5PD9jRl5zglR/DGJPJSWU/SzlvJMSPI21FyokY6hDBtVxUpkwM0LKuRQbMh1EzBYwIkVVaFoFVy0t9H6vBxluBxUUY7TBuS057D2bBhA2bPno27774bY2ONN+C1117DJZdcgpUrV5q9yj4CRfCYQslJ7J2hgUps8FhTyYlRQxrfZIDGJPI4kgM0euVUepjkmKuwSqiuMmiz2LSrkpQcDgEmBCoGMxtJmRyApldOUe0qgLbrsZ/FriIuI++l4PGCSbtqx44dCIIAAxr3U9tQFj/Xr1+PdevWNQkOAIyNjeFTn/oU1q9fb+Ti+hHNTI7BZoA0Sg4NyYm1q3QzOUlP+Jrh40SSo7vxREwgb61NW10F2MnkmLSr4rodAwR2VUImBzCfy8mi5JhGoe2q3DM5tBVWmZScggSP58+bi0q5DM75lBmTRYcyyRkfH8fu3bunfH337t3Yt2+f1kX1KzzPw8TEBACzHY+jeuQAekHFWCVHCC2iE2tXESo5uiSHEZKcOBWnsTZxdRVgRckxWkKeWF3FIvsZqSLJrgLM53Li5lY1z0ew4eVVXZXJriJSkYSIz/B1vI7crkrP5AgutFuBmEClXMb8eXMB9J5lpUxy3v3ud+OSSy7Bfffdh5deegkvvfQS/vEf/xGXXnopzj33XJPX2DcIVRzAdCYneoPR+XDEKjmAHsmJUXK0q6uSnsw0yUJc8BjQH+2QRHIak8j1N7Z+UnKS7CpAGFNzGnOrkq/bZC6nwusop3wGTJMqIUQ2u4ogf5SnkhPwcmS7iW5QkhwWeBAZyX9hcjkH9WaFlbKxtnHjRnzkIx/BBRdc0GwQODAwgEsvvRS33nqrsQvsJ4QkZ3jGKCqG2twDjRtkFHQyOXUkfLA0nsxjOx4XWMnhCZI9pZLTWJ+hojGzTQiRrKT0mJJTT2viJnwA+hUxjcGcye+NSbtqIGFuVYjAsJJTEyx+pEOIULk1MSm2DXmTnCygJTnpVlXrtQXpejw2H8Avpg/JGR0dxZe+9CXceuut+MUvfgEAOProo3tuQqlNUJSPAwlKDhfKDbXqSaRDy66KPpYLAS44yorDAJPkZ92ux0mZHF0lJJ3kBKgMqm8wjDOIhM3ahFIkBEdiCTk3M6SzLkpgKb/LggdGRkgEGVoPmFRW0srHAfPB41SrKgTzCUhOet6FahJ5VpIDSpKTIY/TfG1Bwse9WmGlLSfMnDkTJ598solr6XtQTSCPq64SaBCdAYXqn2QlR8euit88mGDqJCfJxtBWcuKP17Visig5OkjskQMYUnL8RCIFmBnSWc9Q9mvKrkrL4wBmlZykCeTN8xlWclKtqhDMAzBq7sRCoFJPP7dTchoojF01SXJ6bX6VNsl54okn8MILL8DzOjfad73rXbpL9x3IJpDH2FVASHLk14wNHgMkSg7QsKwGy2pPjIlKToGrq8qJGRN9kpP2NGwmk5P+VG5iSGfScM4QDbtKH2l5HMBsJiebkmOY5GT9HBsOAJe9oNHCOe20nGmpu3HIUj4OTPbBCgKUCMqls4SOm68tGMmZNkrOs88+i3e/+934yU9+glKp1PR2Q0maGaoM6SdQKTkDMdVVQCOXozJDODl4rPZvWwoYSgk3N53wcZKSo2NXcc4SA4I2Mjk6SM06GVFyMpAcA7ZYUvl4CGHI4shCcozaVSnl40Cj47HJeU7Z7SqzpdRZ8jghfO5juGJ2CnpmuwoNNYeG5PSeXbWgR0c7KFPkD3/4wzjyyCPxyiuvYHR0FD/72c/w/e9/H6eddho2b95s8BL7B61Mjh27ClDvlZNIchSf7OLKx0PohI85UfCYB8k3eMrqKkC/V06aXWVLyTERPs6i5Bizq1J65ACG7aoMSo6A2Wnkme0qw7aRLMkxDVmSQ4EsIx1ary2GkrOorboqNbBeICiTnC1btuDGG2/EggULUC6XUS6X8Za3vAXr1q3Dhz70IZPX2DdoTSC3Z1epVlgl21VqG1aSVQWok5zUvhca86uSrKrG93tcyRHIZB0kI30jMFFGniWTY2q0QyYlx3ImBzBrWU3kZFc5kiOn5BSlIWBoVx04cADj4+M5X012KJMcxhhmz25s1gsWLMD27dsBAIcffjiefvppM1fXZ6CYQA6kKDmKm3CNIJMT1yOnuayiXcVFKbHvhdBQQ3qe5GQ4Xl/NydCe3pKSY2q0Q5bgsU6Lhm5kUXIAsxVWcsFjc5AhOQFBQ0ApkuPRdD2Wy+QIilZF0hgZHsbsmY0Aei9ZVsok58QTT8SPf/xjAMDy5ctxyy234Ac/+AFuvPFGHHXUUcYusJ9AFTxOIjmqXY8pgsdJlVWAupKTOLcKepPIU0mOZtYkbqRD6/y0wWMABnI5GYPHmkjqdtw8j0Ulh0GAG9p9svTJAfS6mHcjO8kxrORk6HYcwik5AIQojJoT5nJ6qcJKmeRcd911TV/uxhtvxHPPPYeVK1fiwQcfxF/8xV8Yu8B+QjN4PGqO5JS5j1JC+S5JJkeRjMQN5wyhGjxO7HYMaJWQp5EcwfUawuZeQg4TSk5xMjnG7KoMmRzAnGWVJXgMmFVyesOuMjyPjJcS83vdECnqswpY4E32lpI5pgBSDnqzwkopNu77Pm655RZs3LgRAPDGN74RTz31FF599VWMjY1pN/3qV1BMIE/K4wBqkrovePITaq8pORrN6NJIDoSA4AKlitrvfO6ZHMCAkmMnk+Nlqa4yVEKexa4CGiRHu/ZHiNS5Vc3zGXyi7wW7yrSSI6PiADR2lYxVFaLRK6cIXY97j+QoKTmDg4P47//+7ylfnz9/viM4CWiVkJubW5VkVQFqSk6iVQWQBY/VlZwMv3OKag5LmFsVQjmXkzCBvLW2bnVV+vG9ouRUs0yOFlz6KTkKWewqwEwuZ4DXE9XYdpgiOXXBwLK+T4Jr9cbqRk+RHAK7SqYRYPOYAgzpBHpzfpWyXfUHf/AH+OpXv2ryWvoeJEpOQo8cQO0mnBg6BpRLStOrq1SDxxk2P0VFJFXJgTrJSVNxGmvT21XQ3jiL0ycH0LesgoBnLpE1YVdlDR0D5uyqzI0AQxi0rKQyOYatskKQHAUlpyhl5NPGrgKAIAhw11134Tvf+Q6WLVs2ZWbVhg0btC+u30DRJyduOGcIEiUnDKJIqnZpdhUjsqsA9V45mUgOV5OSs5CccBJ5SXFIZ5BBCbJRXaVrV3EBeBnILDBZYVUeUj5XVhUHMERyMuZxAHNDOjNbVSGYBwzO0D8xF1IkJ8hbyWFM6/MXBZbSeysKxRnt0HvBY2WS89Of/hSnnnoqAOCZZ57p+J6zrKaCc97sLWCS5MQN5wyhUo2RGDoGlCcTpyo5VMFjQNmuylvJaayvNomcC57NktBQBxqKR/p7qzuksy7KGQ0dgAtfK72QNY8DmOl6LKfkmNnsMnc7DmFIUanU/bTh7p2nFRxMMFRKZvIosiQHmOx6PGyu6zLznJJjE8ok57vf/a7J6+h77N+/vynZjxi0q1IzOQobWCrJAdRITp5KTo/aVY311SaRZ+mRA+gqOT5E0jDX8BwQ4JyjUlHbrLKUjzfPpWtXZaysAszMr8paPg7kaVeZCeDK5HFC+NxX/r3phhLJ8TzAIMkJJLodhyhaddUrr7wCz/MwNKSumNqC2clnDrEIQ8eVgUEMDpn7wFBUV6XaVYBS+DitGSATPHWadeRxWWwfQiVHlSSkhY5b16BG0LJWqwmtjTP75qcTPq5ntKoA/flVhbaruJnePPJ2lSElR4XkGMzl+ExNyTEJqR454TEFUXLmzp6NgUnCuWPHjpyvJhscybGE9jyOSTuPorqqnsXiUNhI0vrkAGpl5EnDOUMokxzC6ioZu0oFmXuMWCI5OrmcWoaRDs3zaCo5RbarADMVXdJ2laFsjEweJ4TJCitVu8okuEJ1FWfF6HpcLpebuZxesawcybEEsgnkKUqOgHwuJ5uSI7+RpNlVgNrTfqYScgWiIAQHz/Bzqg7pLGe84aorOTbsKjtKTpZGgCFMVFdlhZnqKjn7wkQup5fsKpPh47xJjhACgUJ1VZG6HvdaLseRHEugGumQVkIOyKs5qSXkgDzJ4RylLPOAFDaobEqO/AbLWZCpnbHqzUcmk6OCzKoYF1Jh0E5YUnJkMjmam6KfgYyHMJHJkbGrAEMkp5fsKkMkh4ts94pumCQ5PPCUW6QXJZezYJLk9EqFlSM5ltCyq+wN5wwhK29TKDlZrCqATslRsauy5HEaryumXSXTd0hdzZEgORq9cmQyObp2lYySEwihlCNrhzzJycGuypXkGBrVoaDiAIDwzJEclR45zWMLkstxSo5DJPKyqwB5JSdzdZUEyl62DVdWyREC2WbRKBAFnrGfhapdRU9yJN5LCyRHz66Sqa7SVHIkNxPdXI5sJme62VWmgscqoWPArJITKORxQhStV44jOQ4dyGMCeQhZJScbyZHbsLLkcQCAS67LRAkCParkUFdXSdhcQnnjtBU8LmZ1FaCXyykJlukz3I5AU8nxBc+m1raDM8BA1+o8g8fKSo5Ju0qhsipEUewqp+Q4RIJOycmQyZG8KVLYVWmNAEPINgTMOlGYkuQol5AXyK7qJyWHi0AjYyRPcnRyOTI9ckLoKjnSeZwQBtScPDM5gUSFXjuE72ce85F6Df1gV/XY/CpHciyBYqQDkD7WAZCvrqpnaO4mTXJSeuSEkC0hz9TtGFCyq5gEyVGxrApTQg6dXjl2lByZTA4glCusGOPS/5Y6dpWsVQUYIDmqSpcu2eA888NOO/JWcgBAZLx/pYEpNAJsHVsQktOm5Jgif5RwJMcSQpJjUskpcz/T9GKZTI4QIlclR3YSeZZuxwAgOIeQJHtZlRxAgSRkmEDeug7V6qr+UXJkOh4D6mXksnkcQM+uGlRScvQ2FunQcQjNbIyKigNMll0bCB8HGR8qIq/BN5NJUmkEGELVFjeNBWMNklOv1/GrX/0q56tJhyM5lhDaVUYnkGdQcQC5TI6HjBOYpZWcjMFjKiUH8mXkMiRH9uk/q4rTuA4bmRzVG6jEe6SoQASihCAjIWyeSzF8LGtVAZokJw8lJye7SpXkAGbUHFW7CjCXy2EaweOiKDlDg4OYN7vxsN4LlpUjOZZAETzOGliUUXIyhY4BsuoqKiWn8WK5a87S7bj5WkmSIENywknkspBScmzYVVxtbIecVdWAavhYZm5VCC27SrJ8HNCfRK5sV+Wk5ABmSI6vZVcZIjmaSk5R3KEFPRQ+diTHEiiCx1kaAQJyJCdz1YVsn5yMnrbskM5M3Y4nIRs+llJyCElOY3058hfwQIpQqISnhcg2nLP5eqh1ba0qPIH3jF2loOQwATCNc/aaXQWYUnI0SI6nb1cJwbWUHAhRGDWnlyqsHMmxBIpmgFlCx4CcXZVdyeFSJaWZq6ukS8gl7CpJolAUu6pxLXIbk+z7qKbkKAwaVMjlyIx0CKHaEFDFrrJdXQXoWVa52VUK5eMhdEc7CCH3QDTleANKDtPodhyiKL1yFvVQhZUjOZbQzOTMnGVszSzl4wCRXQVIqTnZ7So6JQekSo7czSdr6Li1vqSSI/l6FTtMJo8TQqXrcV0ydAyoj3ZQyuRYrq4C9Hrl5FVdpaXkaKpIAS9n6qcVfwH6wWcdq6q5RkF65SzooYaAjuRYQK1WQ63WeGIzGTzOmsmRKSH3skwgDyFDcjLbVUzKZuEySg4lySFXcuTtKikobZryT/YqdpWKkqNsV6lkcrSqq9RIjo6So25X9W7wWMeqAgBuoLpKZ6RDc42CKDkLe2h+lSM5FhCqOAAwPGOmsXWzVlcVQsnJWF0lIBc+lgkeS9tVEsFj2eqkrBPIm9dCTHLUGhrK3/jV7Cp5JceuXQXl+VWqSk4+dlXvZnJ0SQ6M2FUGlBw3iVwajuRYQDN0PDoL5Yp6GWM3Mis5EhswnV2V/bUyGyGVXcVZACGhatEHj+U2Jl+WTAgo9MqxpOQobFCq86tUgseAUM7lqPTJAdR75XAh5D7jHQcHWpkSPZKjOapDk+QYyeQYUHJ4QeyqXppf5UiOBbTKx81OIM+q5AhkV3PkSE7G1wqBskSPGjklhyZ4LGNVAX1gV0FFzSmukqNcQq5oB6jkcso8QFnVVlN8oq+KQL1LrRBaltVATf1Y3UyOTvk40PhsqIyGaYeZTE5RlJwGyXn11VdRraqpkbbgSI4FtMrHzYWOgewl5ED2XI7U4L6MG0nZZ1KzhGQ2aLkSckKSU7QScpUuydLqQHEzOap2lUomB1DL5QzyCaVzAeq9cpStqhCKZKMUMJQ0NuhA+MqWIGDAroK+mqMz0qG1RjFIzuyZoxgeGgIAbN++PeerSYYjORZAUT4OZLerAColJyPJkZxXk/VpX4jsAzoBSNpVxCSHurpKQTGRV3JslZAr9smRJB6cC6XqL0CN5KiWjwPqdtWE7ngExWyMjlUFNP4pdUY7FILkGFByimJXlUqlnsnlOJJjAXQTyLN/aLL2yikCyQkyPm0yUZIqCxUS6gYL5FQK2c2R3q5SyF1IkxyFEnIFBaKuuEHJqjmqKg4AqDxgq4aOAT27SguKdpUuyQH0wsd9Q3IYVywSMI9eqbByJMcCKEY6AEBF4kOfq5KTsbIqRNanfdnmXpR2lZBsuU5NclSCmvLzq+xkclTGOgDy4WPVPA4A+DKtFyahGjoGNJScnOwqnUaAIXRyOUZIjkbXY8E5mGYJfgjXK0cOjuRYAMUEcqB/lZysDQFlQseNA+Sqq2SROXwsMYFc9XqUpH0bmRxJxasuSmCKTdxkw8e+5ADXjmMV9h0dJUem91U7lBsBhlAlOdNcyWFBXbvbcXMtV0YuBUdyLKDV7diwkiMRPM6q5Mg1A8y2KcgqOVmtFmklh2cfdCmr5ADZlRBZFQewlcmRvXmqBY9lAqR1jcnRspPItZQcBQtBtREg0HAWVYhOL9tVqqMdAl6Sy+7FQGh0PTZRPt5cqyDh40VOyXEIQVFCXuY+ShKEJEs1BhcCPoGSUyFScmS6HYfIWgbKFW7mWfMmKiRHZhK5EELJFpJRcmSHczaPg1y4V6WyqnkuyQ1dZQJ581il6iq90lsVyyo3u8qIkqPYFsCAigPo2VUmGgGGKEr4uFcmkTuSYwEkE8glrCogm5Ij3SQscyaHxmqR6XbcOijbz6ii5GStsFIhOY31s127ahWKXKBR/YYvEz6uKlRWhZC3q3QyOSpKjt7TvUr4WNuuyqm6ClC3q4yRHB27ykDouLlWQZScsFeOCx47kASPZcrHgWyZHDKSk3E4Z4jswWNCJUdipEPzmIwkQZ3kZLt26W7HIaSUAY2nWonrUykfDyFbXWXbrhrgeiRHpVdOfnaVfui2p0lOHyo5YSZnx44dyq0XbMCRHAtoZnJGDSo5EnkcgErJyZrJkS0hzzqxXMH2yWxXESo5ivmA6ajkeBoblLRdpdWsTkCq4yX0MjmAvF0lhDBjVymoVpW6/hRv1eoqYySHMYXc2uQ1mMzkFCR4fNC8uSiVSgiCAK+88krelxMLR3IsoFldlcME8hBZlBypbscAIHgmoiOv5GSdWO7sqigodTuGbAm5HSVHz66S+zfUsasEBGQfsPUzOXLXWxNMfaRDCCGkZtaF6Ae7ClBXc/oxeDwwMID5cxt7WpFzOYUiObfffjuOOOIIjIyMYPny5Xj00UdjX3vnnXdi5cqVGBsbw9jYGFavXp34+jzR6nhsbqyDLMkhUXKATDc8OiWHMnhMZ1fJTiBvrp9ZyVEdwFg8JUcneGzTrgLkLKsKr0sVDkSeT1LJ0baqQkh+NspegJIB9YHxQGm0g1GSoxg+7sdMDgAsPKj4FVaFITn33nsv1qxZg7Vr1+Lxxx/H0qVLcfbZZ8fKYJs3b8b555+P7373u9iyZQuWLFmCt7/97YV7sznn2LdvHwCz1VXyweP0D0ZdoVomE8mRVHKEEJmGdKooOVkaAgrBlfrkUJaQAzKZHPXNLLuao0FyrFVXSSo5GtVVgFz4WNeqAuSVHHMkR+7f3oSKAzTMQBXLymcFUHIMZnIEE9IDgamwcMyRnMzYsGEDLrvsMlxyySU44YQTsHHjRoyOjuKuu+6KfP03vvEN/Mmf/AlOOeUUHHfccfjKV74Czjk2bdpk+cqTMT4+3pSIzVZXFcCuAjKRnIpCf4ksloaKkpOlIaCKigNk37zpq6vUm9plV3Ps2FU6fXJkMjmcC+2sgxTJ0QwdA/IkR3tuVQjJz4cpkgOoWVZ521WCc3DJMTFpKFr4uMgVVoUgOZ7n4bHHHsPq1aubXyuXy1i9ejW2bNmSaY2JiQn4vo/5k2Vt3ajX6xgfH+/4YwNh6HhgcAiDQ8PG1pVpBAjkbVfJr5tlfhWVkqNSWQX0fiYHkAkfq29ctuwqIThExt9pXasKkOuVY0LJydrFPIQxJUeSaJgY6RCiF0mOienjU9cshmXVC71yCkFy9uzZA8YYFi9e3PH1xYsXY+fOnZnW+PjHP45DDz20gyi1Y926dZg7d27zz5IlS7SvOwsoGgEC8nYVFwBLuQlTkJySzxQGP9IpOVmGdCorOT1eXQUAyExA1OV3WyXkQPZeOSZIjkwmJw+7SruyKkROdhUg3/W4cd/LmeQYzOM01ywIyemF0Q6FIDm6uPnmm3HPPffgW9/6FkZGRiJfc80112Dv3r3NPy+++KKVa2s1AjQXOgbkg8dAupqjRnKSj5ENHYdIU3KEgFKr9izBY6ZsVxVEydGwq6xkcjJuzlwAnuYGldWy0qmsaq4hoeQMaFZWAUDAhFQQV7sRYIgesqtMqjgAIDwFJclgZVWIrA9U1FjYA6MdBvK+AABYsGABKpUKdu3a1fH1Xbt24eCDD0489rOf/SxuvvlmfOc738HJJ58c+7rh4WEMD5uzi7KiVVllWsmR32QCxjFUif/QUyg5siMdQqQ97St1OwaALMFjaiWnoB2PGycpTianJsoKtTSd4MJHFi3IvpKjv/EJNIjOYCXb71Ne1VVG7SrJc5sMHQOqdlX/KjmLXHVVNgwNDWHZsmUdoeEwRLxixYrY42655RbcdNNNeOihh3DaaafZuFRpUE0gr0jaVUC6kkMRPJatrAqRNr9KpdsxkM2uUlVyIDJUPShMIA9RrOoqjUxOxoC2rlUFSCg5mpVVgGQmx4CSA8hZVv1gV+Wu5PgKCnof21VhJmd8fLxZRVw0FILkAMCaNWtw55134u6778aTTz6JD37wgzhw4AAuueQSAMBFF12Ea665pvn6P//zP8cnP/lJ3HXXXTjiiCOwc+dO7Ny5E/v378/rR4hEs9uxwZEOgKKSQ2FXpWy8qnYVlZJDGTwG0tUcVRWnsXYxqquECCCgYYllHCBaN5ClsJrJsVxdBcj1ynF2lT6EH0g3VDTZCLC5ZkGqq2bOmIHRGY2ISFHVnMKQnPPOOw+f/exncf311+OUU07Btm3b8NBDDzXDyC+88AJ27NjRfP2Xv/xleJ6H97znPTjkkEOafz772c/m9SNEgsquosjkeCrNyYiUnLSGgKpKDoSASCELqsHjxrF0JCfLJHLG9braZlNy9Mths6g5NY3y8eZ5MqoX9u2qCe3zAXJKjrnqKkklJ8fqqsDA71A3hOSDG4VdxQui5ADFz+UUIpMT4oorrsAVV1wR+b3Nmzd3/P/zzz9Pf0EGQBU8HhDyHxwSJSeN5CgrOSl2lWomB40ZNKVK/M1Ph+SklWDrkBygoeZUyvEET8eqapzAEslhPPXuo1M+HqKowWMTmRwgO8nxBEOg2WG5Cc4bn/tyhu1DCMNKjmQXa83PWxSE7wFDg5lfT6PkFInkjOGXL28vLMkpjJLTr6AoIS/zQKklfJKS4wsORqLkKFZXpdlVOjevlAorrjhpuXEsNclJvnYtqwr2lJwsdpWRTE7GJ//AQCZHQKS2aWi8UEi3gIhDkNGuMmZVhcj4IFCuB0otJGJPyxm4xH3KtF0FyIePKTI5govCVFgtGCt2GbkjOcSgCB6rhI4BIEiwCJRCx0CGEnKa4DHXaRKXksvRy+Qk34BVQ8et9ZOvXacRYOMEWTYQQ0pOCkyQHJt2FZBNzRngNchOLI89X0Ylx1joOERGkmPSqgohY1n5JHZV9vNzzrSU4SQUZRp50SusHMkhRjN4bHKkg2S34xAsgfkrWVUAWQl5qpKjZVelECiNp9606ioTdlUS+krJMfAUnnV+lZ8hkJ5pnQyqhYlGgDLnAwzmcUJkVDtNWlUhZEgOiZIjMaSTwqoKUbTRDo7kTFO0gsdFUHLiPxTKSo5gjc58MaBScpSDx0CqXcU05szQ21VpJEdzMxPIIDLYUXJsVVcJIRAY2jCyKDmmyseB7EpOXnZVniRHCE1bO25dCSWHwqpqrl2QXE4YPC7q/CpHcojRKiE3OYFcUclJIDnKSo4QiWqOaiaHrBkgMthVOtVVKU/WZc0nS+pMDgCI1I2zdzI5XASppK1hVZkhOVl65RhVcnKzq/JTcoKMn8+AlyGQM8khmFvVWrsYJKfo86scySFGM5NjsE+OSvk4kKzk1HUqLyhIjuCJAUOdJ7Qku4ozP1GZSkPPKzlAhqCorUyOiduTSK2wMpXHATLaVYZ65ADZ++QYt6syqil5KjkUVhVQJCWnWHbVrl27EGQYm2MbjuQQQghBk8lRtKtYQqi0rnMTTHgqV7WrAIAlqEs6Q/eS5lfphgTTMi3UJEe7hBxZcjmW+uQYUHKA9DJyE+XjzbUsKzmMC/AM55zIy66q6/+udCPr7zgVyYEUyaHM5BRDyRmbOxeVchmc88wDtW3CkRxC1Go1eJMhNZPNAFXtqlyUHA2Sk6RK6JWQx1+TTmUVkCF4TF5dZSBAm6pG6D+dp9lVgSgh0HyvQvCU8LFRJcdyJgfIZlnlFzw2/2TfU0oOQSPA1trFIDmVcrnQZeSO5BAitKpKpRKGZowaW7eiWl1FETwGYklOKWAoaZQ5kik5CURAeW7VJFJLyHvArkpXcvRv3GnvkxmrqoG08LGJuVUh0hpuAsCAoUaAzXNmsKxyKyHP0a7yqUgOF4lqcDumg10FFDuX40gOIVrdjmejnNClVhaqmRwuENusTDl4DMSSHB0VB4h/2ucC4FrBYzq7Kv9Mjg0lhz54bGKkQwin5ORYXUXRJ0cieEyFrGpOQGlXFaRPDtDK5RSxwsqRHEJQTSDX6ZYap+bokZzoY3VJThDz9KldFppAFLQzOVzE55Y1JpCHSK+uolVydIdzttYRibkco0pOyu92L2dygPTwcSC4nlIbBR6kN47kgoTkBDnbVUB2kkNpVxWp63GRe+U4kkMIitAxoJ7JAYAghv3rkZzoD7xqZVWIuKd9HasKoFVyGmtE33h0VZzG2vH/TkIIQ0pO0uZlLkiapOaYCh0D6aMd7FdXGSY5KWTDeB4nRMp9qFL3TVXmd4AJnmhlh8ib5HDmpzYe1UVRcjlFHtLpSA4hKOZWAerNAAHLSo4myYlTJXSsKqBRQh537zVCcmI2HRMkJ2kSuRGCg7RMjrkn8yS53aySk6J+GczkcCRXO5VEgIrkJO00pNlVxiurQqR8VihUnBBZcjmkJCdD12OdpqJZUZRcjsvkTFOQTSDXUXJiSA5F8Fg7kxNzTVrdjoHG02WMIqJbXQXQKjmN9aOv3UiPHCAlk2NJyTGaybHXJwdItqxMTR/vOF+KZUGm5KSRnCohycnwMJK3kkNZPh6iKLmcIs+vciSHEK2RDoaVHMXqKiBByYH5EnLVuVUhWMy6Ot2OQ8RZViaUnDglxBzJib52UyQnWckxR3KSMjkmRjqESJtfZTKTA6SQHIONAEPEWdAh6EhOBruKCGlKTsBL2opvEoRXDJJTHLuqFTwWGs1UKeBIDiGogseq1VVAtJIjhCikkhNQKTlAvJJjxK6KITmm+r7EXLtvokcOkKLkmAtSJttVJjM5lpWchPdvwHDoGEhXcnKzqwjKx0OkhY8pVRwgm5ITEJaPhyiKXRVmcqrVanPfKwocySFEa26VOZJTEgHKGoQkSsnxwfXYdxzJmaZKTq/bVbYyOcnBY3O3piS7Kgg0f/ej1kxUcihITl5KTn4kJ13JyZ/kcMLKqhBFUXKGh4Ywe+ZMAMWzrBzJIQTFBHKdPA4ABBEWgVboGEgoIdcMHscqOY7kRMFYJgdIUHMM2lW2lBwRxM4jM63iAGmZHPMkJ82uMt4IMERadRUlyUn5nBaB5FD2yAlRlNEOQHHLyB3JIUR7M0BT0LGqAIBFbMD6JCeuhFy3GWCckqP/axvX9Vi34zEQb1fpTiBvrh+byTHXCyVezbEVPDZ7a4prCGg6jwMk21UUmRwmoh9eQhhvBNg88TRWchiLrXIMYSeTUwy7Cihu+NiRHEJQlJDrNAIEojM5+iSHAxGzr4qs5CBCyRGcG+lrEadQUCs5JoZztk6Sn5IjBFA3qOQ01owhhgbLx0PYVnKA5NEOuQWPa3Ql1HmTHCBdzaFsBNg6R3GUnKKWkTuSQwiKZoDaSg4FyQEiczm6Sg7nDCKio40RJSeC5Jiwqhrr5GVX9YeS46GkU+sXibjwMYVdlZjJISI5SQ0Bc8vk1Oka4aUReiskJ6VXDuXcqtZFiEh1Pg+E4eOijXZwJIcQLSXHZCbHvJLj6UwgDxGxYVU0lRyB6I3QSCYngiiwjJOVU9fOqbrKZCYn7mewoeSY7JHTPFdh7CoikhOz0XEhUDM90qG5eHzWCZxrFx4kIU3J8Vm+Sg5nPoTBh44kFCWX4zI50xAUJeS6wWNGETwGYpQc/Ztc1PwqI/0vgqk/c88rOaZKyAEgNsxK3yfHZOg4RKxdZTl4PECQyQHiw8c1EdD1LREiVs2hzOM0Ti0SSX3edpWN0HGIolhWjuRMMwRBgP379wMw2wxQpxEgEKfkGNgcu292XKBk4MMXqeQYCR5HkCcD3Y6B/qiuilJyhGAQMHcOLjh4hIposnw8RJxd5Wv2copc03LHYyBeySGrrAoRo6hQkxwgWc3Jm+RYsarCcxUkfLygoPOrHMkhwvj4ePPvJsc66NpVdJmczjVMSdXdSg4XZpQcQankRE0iNzCBvLm+jRLyyE3T/MYVZVlRKDlxdhWFksOEiMySlbmPMhHpiOuVQ1ZZFSLG4s2T5HBh5kEoDYkkx0LoOERR7KpFkyRn9+7dqNft/fxpcCSHCGHoeHB4BAODQ8bW1Q0eN24AnTdgCrtKt7IqRLeSY6TbMRBZXWWK5ECIKUqIKRUHiC4hZ5wnDoaURXQmx/yNK0qpqxM8hcfZVRSZHCA6l0OVxwHiSQ65kpOTXQXEkxwbKg6QPNrBRvl481wFITlzZ8/C0OAgAGDHjh05X00LjuQQgaIRIKCfyQGmqjk0JMeMDcC6btImuh0DMXaVKZKDqb1yTJKcqEnkRlUcIEbJMV8SbEvJsVldBURbVlSVVUC8XUVWWRUiTskhnFsVIm60g43QMeDsqm6USiUsGJsHoFgVVo7kEIGiESAAVDTtKsAOydEdzhliqpJjiOQEbIqhYJTkMDqS01i/830xTXKilRzzJCdKyaHI5MSNdugfkhOj5JDbVdGfmYEqXY+cEHFdj60pOYl21fRTcoBi9spxJIcIFI0AATNKTnd3VDMl5DRKTvfmbdRr7yIKvUVyOt8Xkz1yAOScyaGwq6ZeO2M8tju1LqIC/hTdjpvni8vkOLuKDML3YyvXbCo5SeNRbGNhAcPHjuQQoVU+bi50DOhXVwFTRzuQ2FWmlJyuazNSPj6J7oaApqqrgKnl0aZCx831uwia0W7HQEzLektKDkGfnCi7iiqPA+RjV0WFnekzOdG/E2XCRoAh4kmO+d+fOIiY7KHN4DELIgodckIRy8gdySFCq9uxaSVH/8PT/ZRZ6OoqquAxppIcE3OrmmvbtqtM9sgBAA5M3TN7OJMDDtH1e05lVQHRJGeAMHjcaJw59Zzk1VW5lpDHzXAz+1lLgvCjPxM2lRwIEdu2wjackjONQNHtGAAqKZ0+s6D9ZsiFgD8Ng8eNxQjtKsLgMUCfyQGi1BxL1VVE5b/dag7F3KoQkdVVRD1ymueM2OjysqsGCOdWhYjP5NhUciJs0MCz1u24dc5iWFZOyZlGoAgel0RgpM9Gu5JjRMUByOwqm0oOZSbH1ATy1vrEmRwgIpdDoOR0ESkuAI+K5HR9dqzbVYRKDjA1fCyEyIXklAJmpBFoGgLhR1p0vqVMDhBDciyWj4coSq+cIs6vciSHCK0ScpMTyM1sMu1KjpFux0BPKjntJEcI0WPBYxtKTvcGQh88rolyxLZl6Fxd4WPbdhVlJgeYSnLqghntnRQJwacQHRtWFdCYKhH1e28reAwAiFRy7DfCK5qSs337drpxIpJwJIcILZJjLnis2wgwRHt1lTklpyuTY6gZ4FQlh8auEknDBhXQD3aVDSWn266iyOOE6FZyKElOVHUV1dyqEN12FXnoOEROJAeYGj4WwvA9IgU8YhJ5HkpOUXrlhCXknudhz549OV9NA47kEKFpVxksITcROgY6lRxjJEfwDqJT8aiUHBq7ihmsrAIilBzi6ioKu2qqkkMfPKbK4wBTMzk+ZSanmzALYV3JIbeqQnRVWNloBBiiO5cT8DIE7JEcRDzMWQ0dh+csSBn54MAA5s1p7HlFyeU4kkMEio7H5pQcApIDdFhWpjI5QnSWkZt8SmufX2XSqgIA0XXToVZy/IgOztpoI2pCcAgKu6prSCdF+XjrXPbsqkCgIy9S4XWUyIy4yXN2kxzqyqoQvDhKjlWrCgCPqK7Kg+QUJZMDAIsOKlb42JEcIjRLyA1WVxlTcto2rzoMfjjaSY7B6c7tVozZZoCtdXlMvw9V2LerqJUcumqZdjWHohFgCJvBY0B0PExQNgIMMR3tqu7RDrZJTlQmJ7DY7ThEUewqAFg4VqwyckdyiEARPK7ETFKWRUARPAZaJEcIY5kcoDeVnA67yuAE8tb6naFp+kwOHclpz+WQZnK6S8iJn37bLStqqwookF01jZScqOoqnoddVSAlJ8zlFKXCypEcAgghSErIzWVyCILHQDOTU/ZZRCM5dYQbOOMw6re3Z3KiJntrrc1bXUhNqzhAp5IT1WvGBDr75NhScijtKnuZHKCL5BCXjwNTlRxrdlWBgsc2y8eByc95VyuKII8Scua6HsfBkRwCTExMIJj8xS+iXUWt5JjK44QIlRxu2srosKvoyqMpSE77JHLTIx2asKTktPfKqRFuUu3zqzgXU3r0mEbeSk5udlXOwWPbaFdzhBC5lJA3uh4XQ80pWtdjR3IIEFpV5XIFQyOjxtY1FTzmAs3+GRTBY5N5HKClVBjtdowuuyqgqBxqvMcUJKexfuP6SRoBwl4mp9OuslNdRa3iAJ1dj21kchgXHb1Jpodd1d0UM1+SwwPPaCsKGRQll+OUnGmAdquqVDK3wZlqBgi01JwaCckxe3Nt2VVmf10FY01XjUTJ4dQkp/G+kORxAHtKDmu3T4ntqskfiTqPAwCBZSVHoJNYTbjqKitoJzm5qDjNcxdEyTnIKTl9D7IJ5IbsKqDVK4fGrjKs5AgaJQcAMGkrmhzOGaKp5FBcN5ySIw/RrLDyA5r3rB2dSg49yQE6Las8lJyyF6Bk0TZhPOhoQZALyfFa94488jghilJGHio5r732GiYmJnK+GkdySNAqHzc9gdykktP4QNDYVb2h5ACt8DFNJseOXUXSIwdoSAPNjbr3g8dAq1eODSXHb+OIA8TDOUMEk79zvmAIhKVNj7cagdpUcYDGr2i7kpmPktP6bOTR7bh57oLYVbNGRzEyPAygGGqOIzkEoGgECAAVYW6jCXvleCZvhJM3mwpR8JhCyRGTRIEb7ngMtJQQ+kwO3RO7YPQkJ1RyAlFCQKR6hQiVHMoJ5CH8ts+WbSXHmlUVYlLNKVsMHYcILauAl8CJf3+iINoe6vJoBNg8d0GUnFKpVKhcjiM5BKAoHwfMBY+BViaHpITctF3VVHLo7CpKJcf0BPLW+mEmh9B6sajkVEmtqgbC8DFtI8AG2h+sbWRygBbJsWZVhZj8/AxYVnKAdpKTz3Ym2uZXuUxOAwscyelvNJWcAttVjAsEgoMRKDnG7SqqEnJQ21V0JeSN9ftLyaEc6RDCql0VEkTBjX52E885+e9lrXw8xOTnx7ZdBbS6HudGctqDxzkqOd3z8vLEwjFHcvoaFMHjkmAoG7xxBVyYDR0Dzb4zZEoOkV0leOOPadBXV9EGjxsnsaDkNPNhFpScyd95e9VV4WBOOxtQU8mxbVdNkrhKzQ6Za0fuSk5HdVWemZziKDmLClRh5UgOAZrBY5MjHQw/CTIuzJaPA2RKTjOTQ3ETCwKSyirAXnUVWfAYjUGjQnAAdOfgnEMIQR46BgAxuSHasKsEBHxhp0dOiGYmZxopOWFDwNxIDmPNxpx5KzlF6Xq8YLIhYBFGOziSQwCK4LGpbschAs7NKzmCTc6tMrtuQKnkBIzEqgLsVVdRjXUAwvC03zFRmwKMM+Ly8QaawWNLT70BF9ZCx0CruiqvTE6lbvm8yF/JARpqjhA810wOhCiMmlOk4PFA3hfQj6CYQE6h5NRNl5gKAXBmfKwDFwJccJLgsWAMgqCyCrBpVxFuLJMkhxqccUuZnMngsYXqKqAx2sFW+TiQf3VVLkpOQUgOL4vcuh2H4AEHBvPXLhzJ6XO0Mjnm7CrTwcWAi0kbwjB4YFzJARqWFaN40g8CMiVHMEEygTwEZ6xB/gj7oTSCx/Q5C8aZnUwO98G5ALPUsM4XdpWcMOycm5KTYybHZzmSHM8DK9M3mExDUXrlhCRnx44dYIyhUqF/gIlD/pSvD9HK5JgLHg8IszIoYwTBY6BBcgwrOUBDrSBRcghJDucCYHR9OzgLEDDiGysXAOgleM44qhYyOVz41qwqoEE6bJWPAw2FlgmRTyZHiJzsqrCVQr5KTp6NAEMUxa6aP28eyqUSGGN45ZVXcr0WR3IIQFFCbtquCrgw2yNnEiWvjhI3/zThMQ4BGruKiuRACMq8LoQQ8AgGi3acg9mxq6wpOSKwSnICIawGj4GGZWW9uop5KNeDtmo8i6fmDUUzd5KTZx5nErYUyjQMVCqYP28egPwtK0dyCNCyq0wGj01ncjgJySnXaTZdj+jDS6nkAAAC2g6sPjHJgSW7ijNup7pKMPgEdmocfGFXyQGAWsBoVNok8ACVan6bfJ35NHZ2RhRFyeEFsauAlmWVd4VVoUjO7bffjiOOOAIjIyNYvnw5Hn300cTX/8M//AOOO+44jIyM4KSTTsKDDz5o6Urj4ft+cyiZ2eCxYbuKSMkp12ludFQkB4yRjHQIQb3XeMQlq43qKjuZnJqlJ3HPt5cbsW1XAcA+wt/nJFSq+Q1jrOVs0wjPz7V8PERR7CqgOOHjwpCce++9F2vWrMHatWvx+OOPY+nSpTj77LNj/bwf/vCHOP/883HppZdi69at+L3f+z383u/9Hn76059avvJOhHkcABgZNZjJMa3kCJjvkwO62TVUn13BGBgj3PSIn6w86g3NFslhHHULSg7QeBCxBV8IDFi2q/ZR/j4noHLALplrRzVvkuP7uTYCDFEkkhP2ynEkZxIbNmzAZZddhksuuQQnnHACNm7ciNHRUdx1112Rr//85z+P3/qt38JHP/pRHH/88bjppptw6qmn4otf/KLlK+9ESHKGRkZRGRg0tq7pTA4A1Ah8eyq7yif0+rlHeHPqdbsKgDCsIkahxkqwdXu2SXICy9VVALA/NyUnv03ey9mmadhV+Ss5RbSr8iY5hSgh9zwPjz32GK655prm18rlMlavXo0tW7ZEHrNlyxasWbOm42tnn3027r///sjX1+t11NuslPHxcf0Lj0CYxxmdNRszhsw9mQ4NlFEanGFsPQBgpQMYGTC75ogoozxjxOiaAMArQxgs0XDychkYGDZ/zQBQrgyiVKZ7luAlgcFBc2Q6CpVyGaiY/T3phl8awWjFzjPXhOCYMWhHNRosCQxUBoCKvVtttVI3/rnOghEGks9+FvgDgxjM85m9xFEqC7L7iAwqg8Mo2R/GPgUHL1oEwJEcAMCePXvAGMPixYs7vr548WI89dRTkcfs3Lkz8vU7d+6MfP26detwww03mLngBFSrVcyePRuHLV6AP151tMGVTa7l4CCDd1o5y9usnCUPvNn6GVdZP6ODQyf2z56Pyhe/DEbd5iIFhSA5NnDNNdd0KD/j4+NYsmSJ8fO85S1vwfj4eHPooIODg4ODw3TDqlWrUK/Xc20ECBSE5CxYsACVSgW7du3q+PquXbtw8MEHRx5z8MEHS71+eHgYw8PDZi44A8qEFoWDg4ODg0ORkTe5CVGInXhoaAjLli3Dpk2bml/jnGPTpk1YsWJF5DErVqzoeD0AfPvb3459vYODg4ODg8P0QiGUHABYs2YNLr74Ypx22mk4/fTTcdttt+HAgQO45JJLAAAXXXQRDjvsMKxbtw4A8OEPfxirVq3C+vXr8c53vhP33HMPfvSjH+GOO+7I88dwcHBwcHBwKAgKQ3LOO+887N69G9dffz127tyJU045BQ899FAzXPzCCy90WEBnnHEG/vZv/xbXXXcdPvGJT+BNb3oT7r//fpx44ol5/QgODg4ODg4OBUJJiJxnw+eE8fFxzJ07F3v37sWcOeZmTDk4ODg4ODjQQWb/LkQmx8HBwcHBwcHBNBzJcXBwcHBwcOhLOJLj4ODg4ODg0JdwJMfBwcHBwcGhL+FIjoODg4ODg0NfwpEcBwcHBwcHh76EIzkODg4ODg4OfQlHchwcHBwcHBz6Eo7kODg4ODg4OPQlCjPWwTbCRs/j4+M5X4mDg4ODg4NDVoT7dpaBDdOW5Ozbtw8AsGTJkpyvxMHBwcHBwUEW+/btw9y5cxNfM21nV3HOsX37dsyePRulUsno2uPj41iyZAlefPFFNxfLINz7Sgf33tLBvbc0cO8rHYr+3gohsG/fPhx66KEdg7ujMG2VnHK5jDe84Q2k55gzZ04hf0F6He59pYN7b+ng3lsauPeVDkV+b9MUnBAueOzg4ODg4ODQl3Akx8HBwcHBwaEv4UgOAYaHh7F27VoMDw/nfSl9Bfe+0sG9t3Rw7y0N3PtKh356b6dt8NjBwcHBwcGhv+GUHAcHBwcHB4e+hCM5Dg4ODg4ODn0JR3IcHBwcHBwc+hKO5Dg4ODg4ODj0JRzJUcTtt9+OI444AiMjI1i+fDkeffTRxNf/wz/8A4477jiMjIzgpJNOwoMPPmjpSnsLMu/rnXfeiZUrV2JsbAxjY2NYvXp16r/DdIbs72yIe+65B6VSCb/3e79He4E9DNn39vXXX8fll1+OQw45BMPDwzjmmGPcPSECsu/rbbfdhmOPPRYzZszAkiVLcNVVV6FWq1m62t7A97//fZxzzjk49NBDUSqVcP/996ces3nzZpx66qkYHh7GG9/4Rnz9618nv05jEA7SuOeee8TQ0JC46667xM9+9jNx2WWXiXnz5oldu3ZFvv4HP/iBqFQq4pZbbhFPPPGEuO6668Tg4KD4yU9+YvnKiw3Z9/WCCy4Qt99+u9i6dat48sknxfve9z4xd+5c8dJLL1m+8uJD9r0N8dxzz4nDDjtMrFy5Uvzu7/6unYvtMci+t/V6XZx22mniHe94h3j44YfFc889JzZv3iy2bdtm+cqLDdn39Rvf+IYYHh4W3/jGN8Rzzz0n/vVf/1Uccsgh4qqrrrJ85cXGgw8+KK699lpx3333CQDiW9/6VuLrn332WTE6OirWrFkjnnjiCfGFL3xBVCoV8dBDD9m5YE04kqOA008/XVx++eXN/2eMiUMPPVSsW7cu8vW///u/L975znd2fG358uXij/7oj0ivs9cg+752IwgCMXv2bHH33XdTXWLPQuW9DYJAnHHGGeIrX/mKuPjiix3JiYHse/vlL39ZHHXUUcLzPFuX2JOQfV8vv/xy8b//9//u+NqaNWvEmWeeSXqdvYwsJOdjH/uY+PVf//WOr5133nni7LPPJrwyc3B2lSQ8z8Njjz2G1atXN79WLpexevVqbNmyJfKYLVu2dLweAM4+++zY109HqLyv3ZiYmIDv+5g/fz7VZfYkVN/bG2+8EYsWLcKll15q4zJ7Eirv7T//8z9jxYoVuPzyy7F48WKceOKJ+MxnPgPGmK3LLjxU3tczzjgDjz32WNPSevbZZ/Hggw/iHe94h5Vr7lf0+v41bQd0qmLPnj1gjGHx4sUdX1+8eDGeeuqpyGN27twZ+fqdO3eSXWevQeV97cbHP/5xHHrooVM+kNMdKu/tww8/jK9+9avYtm2bhSvsXai8t88++yz+/d//He9973vx4IMP4uc//zn+5E/+BL7vY+3atTYuu/BQeV8vuOAC7NmzB295y1sghEAQBPjjP/5jfOITn7BxyX2LuP1rfHwc1WoVM2bMyOnKssEpOQ59gZtvvhn33HMPvvWtb2FkZCTvy+lp7Nu3DxdeeCHuvPNOLFiwIO/L6TtwzrFo0SLccccdWLZsGc477zxce+212LhxY96X1tPYvHkzPvOZz+BLX/oSHn/8cdx333144IEHcNNNN+V9aQ45wik5kliwYAEqlQp27drV8fVdu3bh4IMPjjzm4IMPlnr9dITK+xris5/9LG6++WZ85zvfwcknn0x5mT0J2ff2F7/4BZ5//nmcc845za9xzgEAAwMDePrpp3H00UfTXnSPQOX39pBDDsHg4CAqlUrza8cffzx27twJz/MwNDREes29AJX39ZOf/CQuvPBCvP/97wcAnHTSSThw4AA+8IEP4Nprr0W57J7pVRC3f82ZM6fwKg7glBxpDA0NYdmyZdi0aVPza5xzbNq0CStWrIg8ZsWKFR2vB4Bvf/vbsa+fjlB5XwHglltuwU033YSHHnoIp512mo1L7TnIvrfHHXccfvKTn2Dbtm3NP+9617vw1re+Fdu2bcOSJUtsXn6hofJ7e+aZZ+LnP/95kzgCwDPPPINDDjnEEZxJqLyvExMTU4hMSCSFG9GojJ7fv/JOPvci7rnnHjE8PCy+/vWviyeeeEJ84AMfEPPmzRM7d+4UQghx4YUXiquvvrr5+h/84AdiYGBAfPaznxVPPvmkWLt2rSshj4Ds+3rzzTeLoaEh8c1vflPs2LGj+Wffvn15/QiFhex72w1XXRUP2ff2hRdeELNnzxZXXHGFePrpp8W//Mu/iEWLFolPfepTef0IhYTs+7p27Voxe/Zs8Xd/93fi2WefFf/2b/8mjj76aPH7v//7ef0IhcS+ffvE1q1bxdatWwUAsWHDBrF161bxy1/+UgghxNVXXy0uvPDC5uvDEvKPfvSj4sknnxS33367KyGfDvjCF74gfu3Xfk0MDQ2J008/Xfznf/5n83urVq0SF198ccfr//7v/14cc8wxYmhoSPz6r/+6eOCBByxfcW9A5n09/PDDBYApf9auXWv/wnsAsr+z7XAkJxmy7+0Pf/hDsXz5cjE8PCyOOuoo8elPf1oEQWD5qosPmffV933xZ3/2Z+Loo48WIyMjYsmSJeJP/uRPxGuvvWb/wguM7373u5H3zfC9vPjii8WqVaumHHPKKaeIoaEhcdRRR4mvfe1r1q9bFSUhnI7n4ODg4ODg0H9wmRwHBwcHBweHvoQjOQ4ODg4ODg59CUdyHBwcHBwcHPoSjuQ4ODg4ODg49CUcyXFwcHBwcHDoSziS4+Dg4ODg4NCXcCTHwcHBwcHBoS/hSI6Dg4ODg4NDX8KRHAcHBwcHB4e+hCM5Dg4OPYGzzjoLV155pbHXtr9GCIEPfOADmD9/PkqlErZt26Z1rQ4ODsXAQN4X4ODg4JAH7rvvPgwODgIAHnroIXz961/H5s2bcdRRR2HBggU466yzcMopp+C2227L90IdHByU4UiOg4ODNXieh6GhobwvAwAwf/785t9/8Ytf4JBDDsEZZ5yR4xU5ODiYhrOrHBwcyHDWWWfhiiuuwJVXXokFCxbg7LPPBucc69atw5FHHokZM2Zg6dKl+OY3v9lx3IEDB3DRRRdh1qxZOOSQQ7B+/fopa3/zm9/ESSedhBkzZuCggw7C6tWrceDAgeb3Oef42Mc+hvnz5+Pggw/Gn/3Zn025tiuvvBLve9/78Kd/+qd44YUXUCqVcMQRR+B973sfvve97+Hzn/88SqUSSqUSnn/++cif8TOf+UzzNe1/nALk4JA/HMlxcHAgxd13342hoSH84Ac/wMaNG7Fu3Tr81V/9FTZu3Iif/exnuOqqq/AHf/AH+N73vtc85qMf/Si+973v4Z/+6Z/wb//2b9i8eTMef/zx5vd37NiB888/H3/4h3+IJ598Eps3b8a5554LIUTHeWfOnIlHHnkEt9xyC2688UZ8+9vfnnJ9n//853HjjTfiDW94A3bs2IH/+q//wuc//3msWLECl112GXbs2IEdO3ZgyZIlkT/fn/7pnzZfs2PHDlx22WU4/PDD8Z73vMfgu+jg4KACZ1c5ODiQ4k1vehNuueUWAEC9XsdnPvMZfOc738GKFSsAAEcddRQefvhh/OVf/iVWrVqF/fv346tf/Sr+5m/+Br/5m78JoEFY3vCGNzTX3LFjB4IgwLnnnovDDz8cAHDSSSd1nPfkk0/G2rVrm9fwxS9+EZs2bcLb3va2jtfNnTsXs2fPRqVSwcEHH9z8+tDQEEZHRzu+FoXZs2dj9uzZAIBPfvKTTVLWfr0ODg75wJEcBwcHUixbtqz595///OeYmJiYQjQ8z8Ob3/xmAI18jOd5WL58efP78+fPx7HHHtv8/6VLl+I3f/M3cdJJJ+Hss8/G29/+drznPe/B2NhY8zUnn3xyxzkOOeQQvPLKK0Z/tnZcf/31+Ou//mts3rwZRxxxBNl5HBwcssORHAcHB1LMnDmz+ff9+/cDAB544AEcdthhHa8bHh7OvGalUsG3v/1t/PCHP8S//du/4Qtf+AKuvfZaPPLIIzjyyCMBoFk5FaJUKoFzrvpjJGLt2rX4q7/6K0dwHBwKBpfJcXBwsIYTTjgBw8PDeOGFF/DGN76x40+YeTn66KMxODiIRx55pHnca6+9hmeeeaZjrVKphDPPPBM33HADtm7diqGhIXzrW98ydq1DQ0NgjKW+bu3atbj77rsdwXFwKCCckuPg4GANs2fPxkc+8hFcddVV4JzjLW95C/bu3Ysf/OAHmDNnDi6++GLMmjULl156KT760Y/ioIMOwqJFi3DttdeiXG49kz3yyCPYtGkT3v72t2PRokV45JFHsHv3bhx//PHGrvWII47AI488gueffx6zZs3C/PnzO64BAD71qU/hy1/+Mv75n/8ZIyMj2LlzJwBgbGxMSplycHCggSM5Dg4OVnHTTTdh4cKFWLduHZ599lnMmzcPp556Kj7xiU80X3Prrbdi//79OOecczB79mz8v//3/7B3797m9+fMmYPvf//7uO222zA+Po7DDz8c69evx2//9m8bu86PfOQjuPjii3HCCSegWq3iueee61BqhBC49dZbMT4+3gxRh3j00UfxG7/xG8auxcHBQQ0l0V5z6eDg4ODg4ODQJ3CZHAcHBwcHB4e+hCM5Dg4ODg4ODn0JR3IcHBwcHBwc+hKO5Dg4ODg4ODj0JRzJcXBwcHBwcOhLOJLj4ODg4ODg0JdwJMfBwcHBwcGhL+FIjoODg4ODg0NfwpEcBwcHBwcHh76EIzkODg4ODg4OfQlHchwcHBwcHBz6Ev8/ODZP+T/ayacAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# plot the individual shells\n", "for shell in shells:\n", " plt.fill_between(shell.za, np.zeros_like(shell.wa), shell.wa, alpha=0.5)\n", "\n", "# plot the total of all shells\n", "all_za = np.linspace(0.0, zgrid[-1], 1000)\n", "all_wa = sum(np.interp(all_za, shell.za, shell.wa) for shell in shells)\n", "plt.plot(all_za, all_wa, \"-k\")\n", "\n", "plt.xlabel(\"redshift $z$\")\n", "plt.ylabel(\"radial window $W(z)$\")\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": ".venv", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.12.7" } }, "nbformat": 4, "nbformat_minor": 0 }