{ "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": "iVBORw0KGgoAAAANSUhEUgAAAjkAAAGxCAYAAACA4KdFAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAo3lJREFUeJztfQmYJFWV7smllq7eu6uazRYElGVkERh5bA99jxFHxRl5fo8HCsogjgujwnMDEQRUHBQGR8EeUYQ36sCMouM38OEoigyKw8jiuAHKYoPd9AIMvVRlZmzvOzczsiKzIiLj3nvOjYjM+/MlVZWdGRGVlRn3j///zzmVIAgCsLCwsLCwsLAYMlTzPgALCwsLCwsLCw5YkmNhYWFhYWExlLAkx8LCwsLCwmIoYUmOhYWFhYWFxVDCkhwLCwsLCwuLoYQlORYWFhYWFhZDCUtyLCwsLCwsLIYSdRhR+L4PGzZsgKVLl0KlUsn7cCwsLCwsLCwyANv7bd++HXbffXeoVtO1mpElOUhw1q5dm/dhWFhYWFhYWCjgySefhBe84AWpjxlZkoMKTvgiLVu2LO/DsbCwsLCwsMiAbdu2CZEiXMfTMLIkJ7SokOBYkmNhYWFhYVEuZIma2OCxhYWFhYWFxVDCkhwLCwsLCwuLoYQlORYWFhYWFhZDCUtyLCwsLCwsLIYSluRYWFhYWFhYDCUsybGwsLCwsLAYSliSY2FhYWFhYTGUsCTHwsLCwsLCYihhSY6FhYWFhYXFUMKSHAsLCwsLC4uhRCFIzl133QUnnXSSmCiKbZq//e1vD3zOnXfeCYcddhhMTEzAvvvuCzfccIORY7WwsLCwsLAoBwpBcnbu3AmHHHIIXHPNNZke//jjj8NrX/taeOUrXwkPPvggvO9974O3ve1t8N3vfpf9WC0sLCwsLCzKgUIM6PzTP/1TccuKdevWwYte9CK48sorxc8HHHAA3H333fA3f/M3cOKJJ8Iw4plnnoEdO3aUYiCaPYbROIa892+PwR6DPYaFWLFiBdTrhVjaC4FSvhL33HMPnHDCCT33IblBRScJzWZT3KKj2suC2267DV73utdBEAR5H4qFhYWFRYGxzz77wC9/+UuYnJzM+1AKgVKSnKeffhp22WWXnvvwZyQuc3NzsGjRogXPufzyy+GSSy6BMgIzSkhwkJ0XiaEXiXQV6VgQ9njKcSwIezzlOh6LdDz66KPw29/+Fg466KC8D6UQKM6KyYzzzz8fzjvvvO7PSIjWrl0LZcC9994rvt58881w8skn021420aAZbvRbQ8ANs9uhjVTa4AazqZNMNZHbKmw47kGLFnJc9WzbesWWDY9A1xw/6sB9RW8V2z/9V//JSRwTszN/QEWLdoDTOD3c03Yc9GEkX1t3taAZYvGYHKsBsbw7GMAq/YG0/j9tt/Dnsv2hDzQfPxxmHjRi3IjgK05F1oND5asNPO+SgJmW3/xi1/AU089ZUlOkYLHsth1111h06ZNPffhz8uWLYtVcRBYhYX/Hr2VBU888YT4euCBB9JuuElv2W1j2CbCe/55tpNUc9YFLjR38uaofMZjD18fVEe54brPQxB4YALPOmb2g5hteeJmFDu2mN1fB8/MPQN5wd2yJbf8Dd7clg9uy+v+nNdtt113Fcf1hz/8IZfXo4goJck56qij4I477ui573vf+564f9jg+343P7Ry5Uq6DbstgNZOoMZ2ZztwwGfKUHmuDw7jItTYwfN6mCI5rutCq9UCbnj+HHgeP5lCPOPwvmZRzDlIcsztT2BnTiSnkR/J8bZuhTzhND2h5OSNXde0VWNLcgpGcrBqCEvB8RaWiOP369ev71pNZ5xxRvfx73jHO+Cxxx6DD37wg/DQQw/BtddeC//4j/8I5557Lgwbtm/f3pVEly9fTrdhtwHgzgexC6/kMJEcvALzWj5wYY6d5Dis23ccR9y44XsN8PwG+37ws/SsQZKDKs6caSVndiv+okZ3KV7XuWchL7hbt+aaHWo1XHAahslsDNZMT4uvluQUjOT87Gc/g5e97GXihsDsDH5/0UUXiZ83btzYJTwILB+/9dZbhXqDHiSWkn/pS18ayvJxzEMgxsYnaNPySHCQ6BBjW4vLruJTclyHj+SUXclBgoNqDjdQxfENKDkNP4CdHt/fe8H+HE+oOUbR3AHgmFHFQjS8Buxw8mlxgeTG27EDggY/SU5TcvCWN2Y6ar8lOQULHr/iFa9IZeFx3YzxOQ888AAMO57vZFEWLyXOEDEpObPOLLi+C/Uq7VvL37kTAs+DSo02wCm8dKZFyPc9aM7OAhcC14eAWSVAgmNCyTFlV835PswaJDmo5NSqhvumOLPt2/iUsV3OuXPQ9Jrg+fj7GgxZ4+cAyY3ngz87C9WETCY3nEYx7KrplW2135Kcgik5FoOVnKnFS2k37LUAXPpFBU90Ldw2MYJmg+VKzXN8ceOA13LAbdETyR6Sw6hCmbWrcIE0QHI8Hxqeb8zaQBXHqF3lNPDFBGiZVVXw4kZ8dflIfRKQ3ES/jrKSs6pTUIPVVRZtWJJTEpKzaAmxkoNyNoOS0/JbQrqmht9sihs1UMXhsqucVhMC3wfP5SEJSHACPxBkp+wkBwmOb0jJ8SGAOd8QyWm5ZqurOmQDWmYXfFRyomTHJPxO9R+qvXlmcvCWN1Z3SA52yG/kaN8VCZbklMSumlxCrOQwZHL8wBcqDrWSE2B1D159s5AcLP3kIQlupyrJYThuRGhVcao5puwqJDgm7KrQqkJFxwRQxTFaXRVWTBomG6GCk4uSs7MYSg6eR3xD5DkJk/UqjI+Nie83bNiQ67EUBZbklETJmZhaUvhMDlpV0a9U8DtkwWeyq7gyOSHJCb9SI1RwAsZgaxg85rR3RHDUR7uK/8ozJDezvhmSM2s6eNxVcnbko+TkYlft7CE7eUBYVUGQq2WFqjFeUM2ssuHjKCzJKYmSMzG1lIHk0C4qoYLTIN5uqOBwKTlcmZwwj8NGcjrHzankIMlBEsKp5vh+C8/QxqqrTCk5nh9A0/HN2lWhkmPYrupmcpwcMzk52VXYZwttY/F9jpZVS1wEBjCzapX42ZKcNizJKZGS06LMXgi7ikfJIberOgoOC8npnKA4iI7XIQZc4WNTJCf6lQN+J3BsQsmZ9dqEw0SFVajgYBm5sR4uIckwTDZCJSf8OkrBY6ysCpFnhZXTaL/2oZJjw8dtWJJTEiUHg8dNl/AD5HaqMDyn+HZVM7Srmix9chAuQ3h3GJScsEcOZ6+cMItjJnjcUXIM2FVhVRXyG2OWVStnuyoHJSfIW8mJWFRRwmMazbn26zCz2io5UViSUxIlZ3LxUmhSLsSh2kLYNKzp8pCcoEMWwq+UCEPHqOjwZXJ4lRyf0Q4xo+S0FRwjfXLCTI4JJSfydzFmWTn52FW5KjlhdVVeSk5z/gKgFfne+HF0XgebyemFJTllKSGnJjkhuSG0rJo+E8np2FVcwePoVxaSw1Vd1VH2uEvIo1854HVaDuCATr/zHuJCqOCYyORE1RtjvXLCz3VIdkahT04ncIxVmIGBSsB+RNWbPJWcVteuskpOFJbklMmuopS8Q3JDGBIOszhcdlXQ+UqJsLKKo4w8VHAch9euAq5mhp4nbqYyOe19NoamuipaOm5MyQmDx57L0gcrdpdeC7zOBHnTdhVWFIVKTl6WVavHrsoxeGyVnFhYkjOqSk5IbiiVHKZMDnY7jn4lLV122xkNjjJyh7uEPLSruKrDIjkcXiUnQnIYux57QQCNnJQcY71yQpLT/z0josTG8R1wCHN+gxDgwh4JdedhWRUteDzdITnYJ8c31CqhyLAkp0TNAFkyOQxKDnV1VdjpmNquEhZV5wTJ0fXY6yg4/NVVPCfWKLFhJTmRqirO8HGU2JjO5Bixq3BBi36eDZGc/hyOScsqquLkRnKiSk4zf7tqesUKqFQq4jO7ZcsWGHVYklNgoNLQo+RQLWZuq11ZxaTk0PfJ4bGrosSGNZPDXV3F1LHZFMnxIxYVp10VragyUl3Vo+QYWPyQbERL1Q1ZRwtIjkHLqt+eysOuigaPsYDBNzgANs6uqtfrsGp5e7zDH6xlZUlOkYGzR1qdBXLRYiwhJ/rwREkIISHp9snB5m4lsKuixMYtGckRwzk75dBcwWNzSk4kU8FoV0XVG2FdMS9GUWKDnY/Z0V9RZaiMvF+5MVlh1a/c5K3k5GVZ+T4OCJ2/YJ224eMuLMkpgVVVqVZhfNEUHcmJ2kkufQk52lU4x4rarsJqIsrqiWgOh7qE3BWdgn02uypKbMJBnZyZHM4+Ob1KDqNd1afecKs52ASwuy8TmZz+iipDZeS52lX9JMewkuN5Czum52FZOZ1uxyFs+HgeluSUoUfO1BKoVqt0zQCjvXEY7Kr+73URtanCSqui21VRYoMVIEh6KNHfAJBDzclDyeHsejzn9RLBWZNKjgm7qp/UjIJdlbOSE1cynsc08tCqCmFJzjwsySlD+fji9twqnINDr+TQB4/7v9dF1KaitKw47ap+i4pazVlAchgW0X6SwzWaIKrkmAoecys5+FpFlRwjJMfJya7q269ZJWcuVyUnTrXJo1fOQpJj7aoQluSUpNsxgieTw6PkUIWPAw8tqvkrI8r5VdHeOOR21QKSQ5xT6st4cCg5UYsKF20Oywqb/4W2nmm7ilPJwdBxlBPi3DmXO5DaKohdNUJKTpxqkwvJafT+3mEZ+VN2fpUlOWVpBIggs6sYgseu73YbglEqOf2khnJ+FaeSE5aP85GcfiWH166K+5kC/dVUnH1y+kkNZ6+cuJJx9vBxP7kwRDb6lRuzSk5fddVcQ9jDuSo5uWVy5rHGKjldWJJTCiVnCbGS0yRXcvpJDVUmJwwdc9hV0eCx7/rgE4Z3o5UORuwqp6Qkp4/U+EOi5MTZUw1uy2qBkrMjFyUnz+oqlM9MqjmxmZwc5leFwzlD2EzOPCzJKTDme+Qs60reRVVy+kkNFcnpV3JI7ao+YkBpWfWTGnYlh9mu4iI50TwOwhP2FU/2p2EwkxPN4xjL5fQrN9gPC8c7MAIV3P4LHFN2lbCyYwoRwllWJhCn2rTmclByEoLHzz//POzMaTp7UWBJTinsqnYmx/MDGqLTr+QQLCpsJKdPhuWyq8TPhEShv5qq7MHjuJ8psKCaCq/EmSqs+pUb00oOO8mJy+AwD+qMs6bQtqYe7ZLY7Tjm3NVvYZm3q/Kvrlo8NQWLJifF96Ou5liSU6LgMVkuJ6re4EmCwLLiIjn9JeNBi1PJ8fmCx8STyMMJ5GW3q+LsKa7wcTi3yoSSE+12PH8f8+IXR2iYRzvMRdtRGFZzkhQbo0pOTPAYL54oL5hkRjpEYS2rNizJKUUJeduuIsvl9JMaAsuKTcnpIzWU86sWKDmERKFfuaGeRN4fNKYmOXEl4zzB44UnZw4lp+X74PT9PqzVVaaVHDF1POY9xlxhlZS/MRE+TlJsTCo50QnkeVVY+Z4XqxSHJOepEa+wsiSnRMFjOpLTt4gQKDlcweN+u4o2k9N7IqKcRM5eQu7ykpy4cnGeEvKFhMbz6BfIuEoq1uoq05mcJFsqB7vKmJKTEDA21SsHu4wnqb8mGwL2W1UhbK+cNizJKVEJOYJkSGeJlJwFdhURycETlGfSrmKvrqJdQONUGyOZHCa7ai6mcg6VHVR4OBBHaFgnkScpNtx2VYKSY6LCKkgiOYaqq0QeJyHPaLKMvL9HTghrV7VhSU6BEZ1AzqvklMmuotmuG/M6UvbKMV5dxWBXZbmPJZPDYFclqTZcao55JSeJ5AyxXZWgYBglOQkwOaTTKjnpsCSnDCSnU11V5ExOv11F1gyQya7yYlQbykyOt6C6qsUygdykXWWiGSBXr5ykkPEsk5ITN5BzlnNIZyLJ4e2Vk2RLmQke78zVrkojOSYrrPobAYawSk4bluSUwK6ajAaPdW0Jz8FVhDyT0+hbrKjGOvTbVT5RM8A41YYqk+P7HnguXwl5XE8cJD2UvXLM2VVxSs6cQSWHvicPBrbnYkh0XO8cMiQpNs4wB4/zVXLScjcmg8f9jQBDWJLThiU5BYXnebB9+3Z6uyqOfHAoOX6LpKnbgg7Hng8BwWIbR2ioMjlxqg3lJPKkEQ5lJDn9zQCT7mNTcjz6xQg/o37Me9/xArrRLAs2vjOXTE4Rg8d4jqCswiy6XdXfCDDEdMeu2rhxI0vRQFlgSU5BsW3btu73pNVVcaoNQ58cJDiOr78oxtlT/aMeVBBnTVHZVUnWFJWak0RmKOdXxREa8TclJDq+IMKekeBxUrl4XCBZF2kBY7bw8QgGj9MUGxOWVVHsqrgeOYhVK5ZDrVoF3/dh06ZNMKqwJKfgeZyxiUmoj41379e+EoxTbRIaeskgLoPTb2FR2FVxOZ2i2VVJZIYql5NUSUVZYZV05UdJcuLyOFx2VRLJ4eiVkzaIMy6QTIIk5QQ/70y5Iz/woZlwgYQKD9d4DgSquUHK58mEZZVmSRkNHieQHCQ4q1eugFG3rCzJKXwjwHmrCtF0yqHkUISP8SQZdyKjCB97hu2qtPtlkRQy5rarEJSyt59AZjjsqv5ux5zVVWlqDVuFVZJig0SDyTpCtSaAIDmXxKjmDCIxeSs5OPCXspAhDa2U12La5nIsySnTSIciZ3LiSI5uGbkgM3GzaQhITpySg63YKa4+2e0qJx+7Ku1+SiUnycZiUXIYVI5c7Ko0NZbJshpEYjjDxwNJjgElZ1DDPxPTyLHAob/IIYqZlbaM3JKcwpePz1dWIVocdpUmyUHFJo4cJEnZWZGk2JDYVa34CiXfDcqr5DjlIjlp4xuoczlJ2RsWJcfJQclJ62zM1PV4ULiYVclJCNsWRckR/27AskrqkRNixo52sCSndHaVtpITs9BqkpEkxUZXyUlq/McVPKbK5eSm5JQukzNnhOQgAW+YzOSk9MNh6ZWDn2mcXZUEpoaAA5UcxgqrQUM4uYd0ui1vQa+qPHI5ST1yQsxYu8qSnNGzq+bIlRwukpM0cZwik5PU3Zii67HrJASPqUZSJKh5VEoOVmMYITkpAeOkvI4KGqjQJWRHOCaRp/XDYbGrBik1OdlVuWZymO2qLAQmbkI5+XEk9MgJMbPa2lWW5BR+blUvyfH8ABydq8841UY0CFTfZlLAWJvkJFylUPTAiAseU4WPk5Qcqknk3HZVGpEhtatSAsaUSk4akWkioSMuI0+zpFjsqkFKDZddNSBzw6rkDJg0zm1XZZlNlTSh3ERlVYgZq+RYklM2JUdbzUlSbTSuupJKxbXtqpjycUSQcD+NksNpV5WD5KRVUNEqOSkkh3B+1SBLilrNSc3kcJSQDyITXHaVM7rB4ywkpxiZnHklh7Okv8iwJKfwmZze4LH2aIek/I1GLodNyUm0q3j65CAoyj6T++TwZnJ8ogXUnJKTYldRKjmDSA5xLie9uorBwmjlY1cNVHJyJDnYeoKiM7qOFTWo+sqkkrNz586eBrOjBEtyCj+BfL7bMa+S0yheJifJrtLMtviiiso3bleVRclJIzKUfXJSlRxDdhVHGXk6yaFpUyCl5Dg5ZXIImowmIcig1HBaVmVRciYnJmDp4qmRrrCyJKfodlVfCbk+yWmSk5xEJUezaiuJzAQJVVdZkabWFN2uiptAPr+DlH8rmZJDaVeZVHKwI3laxgdnWmkXD8jaUcPYJ2fA4s5tWWUKHpvI5GR4HaYjltUowpKcwgeP40hOsewqNiUnieRo2j5pREbXrsKrdM4S8rSuxigQUHQ9LkImx/dmzSk5hCSnkUEJJA8fD6yumo1tqqmDLB2NG25DjH6gBl78BI5beCXH93xRas4JZ4BdhVgz4uFjS3LKaFc5xVJyksiM9liHBJKjW12VZknplpB7ggQktLr3/dTupFkwyJKisKzSiExaebkMfN+BwHeluyGzKDmEdtVshsWXvFfOICUHiQaxdYTFBoMIDI584Cgjz6rQcCo5WQdwcvbKcZ30bsdx4eNRhCU5ZQweq16ti1LxhA+nw2BXaVdXJTwfMzUa1k+qXaWZyRlkSelaVnmTnCz/rtvtmHpI56BJ45RKTpY+OOS9crKUahOXc2clLxxl5FnyOPwkJ9vfkHMa+aAeOSGmrZJjUTSgFJxeQq54kkxTazSUHJSl4+AFHji++oKYlr3RaQiYptYk9c/JvO0BlpSjGZoe1NWYouvxIKWGguQMUmpMVldRkpwsVhS5XZUlc0Ocy8lMchhyOZmVHCa7yvOyD9/kVHKyWFWIUe+VY0lOATE7O9tdaPqbAWrZVWm5G4YS8kH/Nghp2Rud+VXpwWOr5JhRctJP0EHgg+/TlNwPIjGU1VVp3Y5ZSI6YMj5nnORkVWiG0a6SqZrirLDKEjqO2lW2usqicFZVtVqD8cl2+R+JXcWk5KTZUkkqTxaklYrrlJGnBY/17ar049INHw8LyclSIk6VyxmUuaGsrspCYOYy5HYyAz9fWcK9Q2RX5U5yJKqmOHvlZCY5q20mx6JgiFpVlUqF0K5iUnJ8eiVHZG5SshRcdhVWRGAfHfVtOzkrOfpXjoNIDEXwOEuJOEUux8OSbYPVVWndjlmUnKzdjK1dlY+Sw1hGPqgRYL9dtXnzZmgR9eoqEyzJKdHcKl4lR21BwQoLBwPNCVANHw8iMToVVoP8dJ2yz2FQckxkcrJkbihyOVlUmiYG2YlKrLOEiklJTtZGfznZVTxKTrb3hT/XEBWN1GhJhImLYFetWLoUxup18f3GjRth1GBJTgGRFjrWGuvAoOSgUoOlouQkZwCJ0ZlfNajhn06vHPZMzgAVT5fkIMHBMnF2uyqLkkNgV2UpDxelzkRqThYlZy4PJScvuypHJQfzShyWVWHsqka2361SqcD0ytENH1uSU+jycZNKjtqCMojEqJKcpOGcFPOrBuVudMLHg0iMU3AlJwuBoVFyGtrh5CzIakVRhY9nR0XJyUhe8szkiMfunM3druIajOlkVHJGvcLKkpwSKjktFpLTLBTJGdTVuLx2VbGrq7LkbcwFjwnsqozkhSqXk6W6CjN1OrmvUc7kIGHwM/aH4crlyCg5OGaFYh5eHFoS58DpDskZxQorS3KK3O04ZqQDAmfjuConZTdlgfW9drPAopAcVrsqPyWHP5PjlULJyRIqJiE5XjYyQdH12PH8TBcgeGE/SzQxPrMN5eRDctDOdlM6WyudFyTOfXnbVSqPz3qe8b3sr+uaEa6wsiSnwHbVZMxIBy3LapAlpWBZDaqeUh3SOahEXMeuGqTkFDqTMyx2VYYeOBTVVVmzNhRKTpY8Dvloh6wKDS6IaRc5MruUJC6Uak6WwZw9j5+lV3JkczYcuZyslVUhrF1lUdC5VfFKjjrJGbCwKBASNiVngFKj2idHlIh75bSrAm/wlPHAC8TjVJGFwJDMrspUXaUfPM6ataEIHssEirMM8swEmcxLawfJLmUb/M0Rzs2SVWao7SoV+4mjwiprZVWIUZ5EbklOCUvIlXvlDFJqFE5GA5UcZZLTUB75oGtFadlVA/vkqNtVWVUaHTUnC4HxPE/cVOH7WMHlDJ+SI0FysgzyJM3kIIhCwLIkh1TJkQwSU9tVDv6NJYPEHPOrss6tCmGVnALgmmuugb322gsmJyfhyCOPhHvvvTf18VdffTXst99+sGjRIli7di2ce+650NCcTl08JSeF5DjDreQMUmp8RbsqixWlalf5njfQJ9eZRG6C5GS1onQsq6xVUyR9crIqOQRBYJmqKbIKK5msDVH4WLZiipTkSNpP5CRHQZXhmF+lY1cFTNVeRUUhSM7NN98M5513Hlx88cVw//33wyGHHAInnnii6NAYh69//evw4Q9/WDz+N7/5DXz5y18W27jgggtguDI5ywqfyeGzqwZlclpKH1ZOJSerSqM6pHNYSI6XcSYVPk73hJxdyfGMZnJIeuUggZO5MCEiOdJKDmEZed52lUqImMOukikfR4R9cprNJjzzzDMwSigEybnqqqvg7LPPhjPPPBMOPPBAWLduHUxNTcH1118f+/if/OQncMwxx8Bpp50m1J9XvepVcOqpp6aqP/jH3bZtW8+t+ErOEmK7qjxKzsCxDUEAgUK+JYtKo5rJyZq3Uc3lZK2c0qmwMqLkZFVosFxY07LKSnIolBwpu4qC5KCKI0MCiciGrDJDOaQzkFzcZR/PESLmCB43Je2q8bExWLFs2UhaVrmTHJylcd9998EJJ5zQva9arYqf77nnntjnHH300eI5Ial57LHH4LbbboPXvOY1ifu5/PLLYfny5d0bWlxlLSFXUnLQRhlUEcGg5GAVhofl6ZLIEixWmV81qNuxnpLDTXKKkcnRVnIkbCjdMnKTfXKMV1fJZuiGQcmRVGawQ7hOTy0SJaeZf/B4lHM5uZOcrVu3ihDjLrvs0nM//vz000/HPgcVnEsvvRSOPfZYGBsbg3322Qde8YpXpNpV559/vrCBwtuTTz4JxberCDM5WQgMQwm5qpqTJViscvJyGTM5WbsZq4aPA3c47Co/w0gHitEOLd8HN6PS0cDKNU1rTIa4kNhVsqSFKpOTo5KjkrGhtKxUSQ51DsZROPfNWJJTHtx5553wyU9+Eq699lqR4bnlllvg1ltvhcsuuyzxORMTE7Bs2bKeWxGBi8fOzoeStLoqixXFYFdlfYxsx+P2YxTsqgzln0iEVE5KXsaFX1nJaQ0HyZEhLjp2lUxZuI/zqzQtqyzdjmntKskFn6q6SlJBog0eK5AcwvCxSr4Gy86p1Zysc6uiGFWS0x5NmiOmp6ehVqvBpk2beu7Hn3fdddfY53z0ox+F008/Hd72treJnw866CBBDN7+9rfDRz7yEWF3lRXRrFCqkuNyKDnyC0qWZn+yJCdw3YGDKMXjmJQczDl4rg/1sZrctrmVnMx2ldoJFQdzZi0NN1FdpWtXyc6jQlI0VVM/d8gQFxlri6R8PEclh8quwspEnCxeNiUnfN74JM1y6zTVpqtPd3rljNpoh9zZwPj4OBx++OFwxx139Jxs8eejjjoq9jmzs7MLiAwSJUTZy+PCPM745BTUanVCJaeRm5KTxdKKws/YAyfr41SsKJV5M+yZnIx/c1UlB/M4WT8/Og0BZZScLNPKdUc6UA3plCEuOP4Bx0BoQXZUQ06ZHMd3wFEYGRMbIlY4v1MqOaohYsoKK5U8DsIqOTkCy8ff8pa3wBFHHAEvf/nLRQ8cVGaw2gpxxhlnwB577CHCw4iTTjpJVGS97GUvEz11fve73wl1B+8Pyc4wNwJUy+RksavkFhQRKg4Gf3gbkrmKLFaVzONkg8equZyyB49l1Bm9EvLsJ2mdXjkNBSVHFZ4fSH8mUflZvqhqLniM5wAsQEi5eBq4Cd+VvmgJ1Z/lteWgA1WyQjmJXLXykrLCSrZHTog1I9r1uBAk55RTToEtW7bARRddJMLGhx56KNx+++3dMPL69et7lJsLL7wQKpWK+Ip/sJmZGUFwPvGJT8CwTyDntavkSEPWk53sSTFr1RSbXSVBhpT65HDbVYonYlMkR2Zcg07XY9mKKZ0KKxX7CTM8yxeNKe9TSZlB9UeDbKiGiNGyWj6RF8nZSUZwfEl1MIRVckac5CDOOecccUsKGkdRr9dFI0C8DRuydDvmCx7LkYasWRvZTE7Wqqny2lXMJEelUaSkBWVOydGxqySVHA27SqUkXDt8rEJyMMczaZ7kUFRYKZMcIrtKp3NxizB4LNsIsD+T8+yzz8Lc3JyYFjAKyD2TYyFfPs5XQt6S8rwzkxxJhShr1ZSSXZWRvKj0ynGd4thVKtm0Qio5GnaVLGnRUXJUBm5q98pRCfRqhoBVQ8QUFVayE8ipSY5OhZRTALtq6eIpmBgfF99v2LABRgWW5JSwESDCxWm4MiflLEQj8LHGOvMms9pQ0tVVmZUc+at8zy1nJgfVmUETyLuPxYcpqDkyxEV1SGcQeOD7rULaVTqZHJWBm9q9clQIh+Ykch27SheqthOVXaUzaJNyfpXscM4QGPEILatRqrCyJKeoweMBSo50LierFSURZswaKJa2qzJncuS263l+Zk9dLZPDZ1fJWlAq4WPZiikVNUe2uZ9vUsnRsKtUCIuWXYUBYlReucvOiRQZEiVnVu29gMpwoKE8UuRqSDM5ikrOqOZyLMkpafBYmuRkVWgkrCU2JafJY1dlaQSomslBeygryVGZRC5LWnwFkiNLWlRIjmxzP993wB80joQqk+OViOTIlo/rPq8ISo6G7USh5mjZVS0vsxLLFTxGzIxghZUlOSUtIZcOH2dVciTCx1zB46DJEzyWydnI2lXtPE72k5isZSVLclSUHFnSotIrR2VMg2rXY9kOxqarq+YULC5tRaaVE8khUXLUj50il6NlOQU0XY/xYkplpEMIS3IsSlNdJR0+zqrQSCg5bNVVWe0qlKElArYyFpSsXSVLWhzJ4aKyXYxVyshNKDkqGRsVYoTvC5xHZa66yrSSo0pycgoe563kEJAcXZJCQXJEt2PMTipixtpVFsUhOcvyyeRIKDl8weOMjw8CqUnkMuqMrF0lm7NhV3LcYmZyfENKDqo4OI9KBp4CMdJScnRIjqoik1PwWLeEPPC8zDZ2Ee0qqoaAOioOYtqSHIvilJAvGRq7Clu6yyguMlkbGZIjY1fJlpDLkhZZUmQieGxGyVGoiFMIH6vma1TVHBXCokVyVJURTUVFlaxgZ3SVQb09SozGyB4aJcfNXclpav4eMyM4v8qSnJKWkEspOb7XrsbISckJIJA6wWW1q2RzOTJKjsdsV0krOa3hIDkq1VIqxEiVrMwaVHKwDYR0U09dJQerJxVfGz/woSHZMJTKslKtrKJScrAyU6WtRBStOQolR+91mOkoORs3bhQzIkcBluQUNnicgeRk/dDJnJgkMjkyM6lkSI5MaXjWkLKsOoOl5nhi47OrZDM5siRHbvHEnjeyJz1TmRy/4EqO7wdiRIPS/lTVHFXCgGqI4nOR4OAFiyp0wsc6oeP28/UUEIoScIpeOao9ckKsXrEcqpWKsKY3b94MowBLcgoEtHTmS8gJ7SqZjsMMSo7MY7G8WqanhZRdJbmgSCk/LYdXyZG84pdVflQIi5qSozCKQ4XkGFRyGq6n7KQoh491AsQ5dS3WU3Jm8yU5JKFhCqKkp+TU63VYuWL5SOVyLMkpEHbs2NG9miYNHjMpOTLqTFbVR4a0SNtVkrkWmfCxrDIjO6STO3hsiuQoVVf5BpUchefpVEkpP1en341i+HhOduo5IUkKdEmOpl1FoeTQkBy94PEoVlhZklNAq6paq8PYxCShkiPxwZA4kcmQnKxKjuyoBim7SrZiSsKCcDLOrSpqnxxzSo4Zu0o1W6PyPJ0AsbpdpUE4WiOo5Mw1hEqsipZm6JiquqpFEKCeWTlavXIsySlk6HipmDNCl8mht6tkK6YyKzmSi79MSFk2OChXct4qdSZHpbGfibEOqs9RrpJSeJ5K6Fh7SGcOdpVuGbheJkdzcQ8CrfAyhQqDF1mY38rTrkKM2vwqS3JKOIG8CHaVbDlo5kyObJM8phJybruKv09OINVGXo2wyA3pxCZmftAyYlc1FBcT00rOrApBchoAiqMudOwqXSVHhyTpVlfpWlYks6c0ux6jEqXbJ2cUe+VYklPSbsd5B49lSU7Wx8vbVTIkR1LdYFVyWiwTyHueJ2W3qQ0wlFGARIBYIaFr0q5SUXJ0MjlKBEm3e3ArJyUnR7uqvY2duSo5YjsalhV2O5YZHZOENautXWVRgvJxthJyz8nUR0Oa5GQkWrJdTbn65LQfn/3E5jHaVSrdi8XzJH5fVZIj8zxfod9NVwHym0aCx0rVVVp2VQ4kJ6fRDHpKDgHJ2TmbO8nRKSNvapaPh7DBY4tSTCAPm4m5WU7KMkqOeHwjNyVHJkgs83hPQQ3hVHJkJpGrNPYTz5Ox2xQyObIkRyVbo/pc5UyO4eqqOZVMjub8KeXqKoJMjkyOLwS2lJDN6lErORShYV3bS2f6+CgP6bQkp8R2VeZcjmyX0gyPl+mRI2VXMWVyZPM44jkZSYLnuuBjV2np7bd4SY5bNCVH/STtebNSM6iaiiTHCQJoST53zriSs7OUdhUSHJVtUKg4OtvBiyPZykyOKi2HiORMd5Scbdu2wfbt22HYYUlOEe0qcpJTJiWHx65Sacme1a6SVXFknydbKdV9nsQCaoLkqDT1U7G6VK0q1ecrqTEdNBxfXt3QVXIUbCdVgkIRXiYjOYrBY2FVaczNIlNyCELHiMWLFsHUosmRUXMsySmiXbVEhuR4uSg5hbGrUMrOcOUt2+1YPCcjMfIke+TI5nJMKDlG7CrFTI4sQVK1qkLMGlRyfCQPss/XzeSg6ijZZwc/vzi7ShdlVHKo8jhiWxq2F0X5eIg1I2RZWZJTSCUnW/A4c/hYWslp5lddJWlX4RVWFsuK065yZI9Z8nkmMjlG7CqdTI6E1aVaWRVizgvkFA5NK0PaslIdztmzDbkFn0LFUQ0vF0LJIYJO8NghCh4jpi3JsShLJqeVKXhcJiVHXhXJUnauYle5hbGreKurcKE2oeToZHJkCJIuyZmV6P2DdjGqMTqQLiPXVXIUwse6PXK0lByNqqgoAsVMC1XoWJcwUSo5MyNUYWVJTiFLyJcSKzmyJKdJHjzGhRS7JFPbVeI5GciCipLjYUO9DAuYbCNA2eepkxwvM1FRqXqR75NTDrtqTqIKT6eyqrs/JwclR5Io6c6t0lJyiBQMHHIr24eLWslB29xXJOFU1VUIS3IsSlFCnimTg/675+au5GQd7SBtV+HJK5OSo3CiCoJMCpCr2kgvq5IjOYFclhypqjjymRyN6iqJ5+oGj2WUIJ08Tnd/uSg5OdlVOQaPVS0rSpKjalmJbseKtvioj3awJKeQs6skMjmDwqWyeRwmJScLMUI1QaUfRhZipKLkZH2eupJTjBJy1TyOyUyOzHNllJjY58uQHI3KKqX5Vai44ViHktpVuZMchW2RjHTQJE1tq4qmwmvUeuVYklPE2VVTS+iUHFmrCpFBmm4obHcQMRIER2V8QRaSoxgOzfI89UxOMYLHpkiOVnWVSSVHwu7SDR23t+HJfTYJqpyk7SoiJUfF9gqGTslxc7WqEDMjNNrBkpyCoNVqwVznjSyl5AxaAFVIThYlx2+R21Wywzml7CrF0QhZwseqJMcpSCZHx67C5/oZSIHKaIbcMjkSJEl5injPNmRIDtGCL5nr0R3poKXkEC7wKkoOZfBYVRkiJzmr2nbVpk2btD7/Q09y8CruySefhIcffhieffZZuqMaYRVHXsnhIDkNkhCxrJKjEgoUz8tQkaXSJwfhldmuCrI9V0fJyfp80cxPowoJCRISpWHM5EgpORShY4XtUCk5qADL9NtBKzpw6BZhlUotaiVHZXsOYWUVYuWyZVCr1cQFytNPPw3DDGmSg22gv/CFL8Dxxx8Py5Ytg7322gsOOOAAmJmZgT333BPOPvts+I//+A+eox2BPM7E1GKo1mqEdhV9JgezNYGCPzzI4lKdT5OlIks5k8NqV7XYJpB3n59hETZBcnQqqwSCIHPXY+0+OVJ2FUXw2DWv5EhuhyqTg+cNGcJEmcdRsavw4kjn80el5FAN5wxRrVZhesWKkbCspEjOVVddJUjNV77yFTjhhBPg29/+Njz44IPwyCOPwD333AMXX3yxkL5e9apXwatf/Wr47W9/y3fkQwaVRoB5KTkqlVWZMjmKdhVXnxxuuyqLAqQ6gVzm+bpydWYlRxNZLauGpl2Fc69w+K2x6iqZbeiOdMhZyZG1vijzOCqkSad5H2kmh1jJic6wGvYKq7rMg1Ghueuuu+CP/uiPYv/95S9/OfzFX/wFrFu3ThChf/u3f4MXv/jFVMc61FApH8+WyaFXclQqq7KQI1W7alADQbTXVDM52UrIFRWoziTyWn0s+TGK5Kz7/AxKlK6Sk4Uk6cytmt9GhsGxSFAI5gyhmrO0WjPSJwc/w54fQK1a4R/OGQLtZrcFUB8fvEvPAdens4xkVCFyJUdye9RWlSpxohrOOYq9cqRIzj/8wz9ketzExAS84x3vUD2mES8flyQ5HHYVntDwJFgbI1VyBpaQK3Q7bj8vffHzsamfouQ8yOYSREWDJKAKxEpyMpA7I3aVRo8cmY7JulZVdDtL64NJToNAyQkVoSUTdXNKTkiYMpAcKqtKRcnJ265yNKaGU26TOng8ShVWysHjv/3bv4UNGzbQHs0IQ2UCOcLxAnEVSGpXDXgeH8lRDR43ScYzqDy3reKoKweDrC7VCeQh/AxKg5HgsW4mJ6MapFtZJbsdCiWnvR3XbCZHgjBRWlXSSg7x4o65Pxzqm1ePnFAdllWWOeyqmRFRcpRJzvve9z447rjjRHVVfyn0fffdR3FsIwVVu2qgmqNMcprGMzkq3Y4RQWMQyVFf+AYFj1XzOPPPH0D8dJUcpxiZHJ0eOVIkR2K4pq4ihJ+71AsMCWQOMFNVV0lsK1clR3GoJtU2OewqWfLke55yBWcaRqUhoFYJOYaPscoqSnSee+45kc2xUJ1bJRc8HpjLUbGrxPMa5JmcgX1yBpCVxOfh7KWUoYqqoeMsz9UlOYNatWsHj4tSQu6ZsavIlJwMJIeiskpaEaJUcjLme6jmVnW3J1VdRa9gyFhgHMHj9nbdXK0qhFVyBqBSqcBll10Gb3rTmxYQHdVhf6MMlQnkmSqsSqTkBBpXK2mVWardjrM8V/cKa7BdxUtyPM8TN/4ScjPBY90eOTJKDkVlVXd/o6Tk5Bg8lt0mm5IjsV0Oq6p/ftUwr9lSweM4INFBwoNE50c/+hGMj4+Lny0URzqQ21WKi3DK1RtXJkfVrgqfW52aiv031cqq9nO9fO0qzbEBg5Qgim6n2UrI9eX2LOFlKpKTRRGiyuOI/WXZFg7bVf08FymTk2PwWN6ucvMnOcQ9cvrtKuy0jxfZK1e2Sc+wQVnJiTK/Sy+9FE4//XRBdH7/+99THdtIoUxKjqpdhWWoaaWoqnbVYCVHfTHyMdidsnDqk5wWywTy7vMH/O66VpXR6iqDweNZ43aVa9aqkrGrcgoe4xrjzw2nklMEu2pifByWLVky9JaVspLziU98AhYvXtz9+ZJLLhFfTzrpJJojGzGolpDnkclRVXJCglSv1untqpQeO65mrsVr+VBbVCXtkVMUu4qC5GRRgygyOVnCy2Ql5L5ZuyrTtijLx2XsKmJyhecAvNhJOg/0fKaJ/p4qSg5e3Ojk+aiCx1x2VWhZbduxQ5Ccl770pTCMUFZyzj///B6SExKdv/qrv4KlS+UX6lEHi12F8rbCjClOkpMWPtazq1qpJEUHadVZZa+uoiI5aUM6RTNGhYGuKkRJZoJ46naKGDymagSYs12VVc3hsKpktstRPq4SaOZoBDhK4WMpkrN+/fqBj7nwwgu7qsQwv3B8dpVCdVWSUqHj3zMEj9OsLjG3SuOqLa3Hjk4JuXh+ygJU+ODxgNlXFJkcJDFp2xF5HImhjMnbccAf0HmXLJPjmc3kZNoWuZKzI5fgcdaKrdxJDpNV1d62WwglZ9qSnF788R//MfzlX/5l6gBOVCSuu+46IX1985vfpDjGESshJ1RyVPM4A56rmslJI0i+piKSNhLC07QVOJUch1nJEdtIsesolJxB28lS+p0Vg7ZFRXKafgD+gIoTqm7HiLk8Mjl4IYNqb9pDfFfr815IJSejXVUUJYfXrlo19POrpDI5v/71r0UW50/+5E9gcnISDj/8cNh9993F99gfB//9V7/6FRx22GFwxRVXwGte8xq+Ix8ioNSvZVc55VFykiaRqw7nnH9+i03JSfPlOYPHgac3gby7HVQJxmu5kRyK8vH5bTWgXl+aqCg1iJrziWnZng+LU0Y7UCo52Lm85fowXq+aKR+PbnNymVGrKneSM9cQ41gq1fRr/BZTZZU4Bred96mNVXMLHiOsXdWH1atXi0nkGzduhM9//vNi+ObWrVu708axZw52O8aJ5JbgZMf27du71WpKzQAT7SodJSf+g+X5ntawvkS7SnE4Z9ntKi+N5BCFHvNWcrL0t8kKz0te+OZ8pCZ0/T4G5Xsog8eZtket5GQgTmwkx8mP5ABWbWXYNqddlZVEea4LnmZxw6h3PZaurkLvfdGiRfDGN75R3Cz0Eao4tbExGBufkH4+XgGaUnJ0VJy04HFacDgL/JTyc90KiVS7SpMk+L6XOImcjOSkBK8pMjmD7So6kpO2LSqrKmv4OJPFJAEMMi9fNGYuk4MYQDa4SE6W7bKRHGFZzUKtUz6dh10Vbn/RkvysKoRVcmLwwhe+EP76r/+6G5S1yDd0bDqTo0tykoPHunZVMzkUy2pX6TdnS7KsyEiOk7OSQ5jJSbO+qCqrspAcx/OFxUS6v0Gkibq6KoOSQ10+LrPdgJPkzO7MX8nJQKI4K6uiSs6WLVugqRkZGBqSc+6558K1114La9euhfe+973w2GOP8RzZCEF1AjlvdVVLyLrUJCfp+bp2lZ9gV4lux5oty5NIEiow6O3zkRyak2zadqhITmp1FaVdlUKYqJWctMaClHmczNtsDY9dVQQlJ2+S42RoCMjV7TjE8qVLYHysrR5u2LABhhHSJOcDH/iAIDZf/OIX4ac//Snst99+wrb693//d54jHKUJ5AqVVWyZHCz5jVFddCstEqurNO2qpG7JOnOruttIIAm6oeP57TRzU3JM2FWUweM0wtQwqORQVlblmskZoA4NY/A4q5Ij05VYBa0MJIrbrqpUKjC9csVQW1ZKzQBrtRqceuqpgtj84Ac/EAP+jjnmGHH71re+NdTDvniVHDW7CjM5flxVie6cmxiSxKbkpASHsyDJ7qLoWJpElKhITtIkct0J5N3tJLwG+Dk1U0JuRsmh6nacRRkyruSgsqra2FNDHeLokSO2O4CwoUKKVVBcGESgsKqR4gJJN/PDWVk1Kr1ylDsehzjuuOMEsXnkkUdESflb3/pWeMlLXkJzdCOXyRmQQpNVc3SUHPH8Jr2Sk0C8dLodIwLHhSBGldDN46QRJYo8Tp6ZHFRxqC5ITFVXpXU9Jg8epyhDlN2O57fpms3j5GhXOb4DTgppC3BxZ7xYHtQrx8G/L/PFehaliFvJQawZ8gor6eqqiy++WCgPcTdcrGdnZ21OR9muUlNywvDxov5eKLpKTkxX0rSxDHqZHH3CgOHjSr2eyWqisauckthV8cdPpeIUJXhMNZwzC2mac1yzSk6GDsEsdhXXfjsq0fLacuNWVZbtc1dWFUnJmbEkpxeXXXaZaP6Hig02/Vu+fDksW7ZM3MLv8auFueCxSSVHu7rK56muCtWgat88NQq7ynMDIV9XqhXDSo7HruRQoQh2lckS8jkGKyOV5HCEjrNUVzHZVWLbziwsn8iJ5AxScphDx1n3Ye2qHEjOHXfcAVdeeSVcf/318H/+z/+B97///UM7vdS8XaVBchwOktMgt6vw+X7gQ7VSJbWrksrIKewqlK1xkvlYn1JGFzzOx66iVHJC6wuDjLx2VUqfHKJux1m2N7DcW2V/qUoOk12FSg3aMjF/N9FBWvccokig8lZyuEPH4hg8XzQarSd0I0c4BuyqmQ7JGdbRDtKZnFe+8pXwL//yL/Dzn/8cJiYm4Mgjj4RXv/rVgvxYqEFnpENqrxzt4DG9kpO0DQq7Km5+le4E8rTtlL26ipLkJIWYPY9mOOf8fvz2Nk2UkONYjYRcBnW344Hb5FJy8PdLCAFjHoeyg7RM+Jib5OAw4LR5dyaUnEG9crDRKLap4MbMkNtVysFjLB3/u7/7O3jiiSfgv/23/yZGOrzsZS+Dr33ta6LaykJByVEsIeezq+iVnKRtUNhVsUoOUYVSXC7HdWjsKoeb5OC4g5jXgdKuStoe5XDOQdukzuSI+VUJag5H8Njzg+TSdC4lJ8Wy4rSqBoWa2UnOAMvKFMlJm0bO3SOnX8nBPjnDWBmtXV01MzMDH/vYx+Chhx6Ck08+Gd7znvfA3nvvLb2da665Bvbaay+R90F16N577x1IDN797nfDbrvtJhQlrOi67bbbYBRLyGOVHDzh65acMik5/RJ44HmiOoqF5BAtRnHlpGWxq5K2RankJG2PskdO2ja9IIAmMclJy+VwKDlif0nvVy4lR2x7p9HKqmx2Fb9Nk0akTASPByk5JqyqaCan1WqJWZQw6pmc//W//ldsZRWe4EIWKDvy4eabb4bzzjsP1q1bJwjO1VdfDSeeeCI8/PDDsGbNmgWPxz8GTkLHf/vGN74Be+yxB/z+97+HFSvaTY1KW11Fmcmh8NIZ+uTEKTm63Y7T5leJjscEiMv2cJIcqgnk3e3h8S8yT3IoQ8dp26TukTNIHeLok9PergurFo8v/Aem8Qpp2+Ya6ZBl+yOj5KSQHBOhY8RYvQ4rly+D557fJiwrFC5GmuRMTU3B7rvvLghF2k0GONn87LPPhjPPPFP8jGTn1ltvFeHmD3/4wwsej/c/++yz8JOf/ATGOi2pUQVKA87liM7m2LZtGwy1XVVgktO/DZ+ILMQ1FKRq6BVrVxFVV8VNIqdUcdrb83JSchpGlBzqPE4aeUJbKXEoriYS7aoBVVDDp+Qw/r5ZlBxjdlX+JCe0rEKSc+ihh8JIk5y///u/Jz0AVGXuu+8+OP/887v3VatVOOGEE+Cee+6Jfc53vvMdOOqoo4Rd9c///M+CeZ522mnwoQ99SHRjjsPll18Ol1xyCYxM8Fg3dMyYyeknOUnDNWURV6FFUUIutsNoV8VNIqcnOfyZnHglh8GuitkmdR6nu90YksNRWTW/7aRMztzQZXLyV3LSSA5/ddWgKi4TjQBDzKxcBY88/vuhrLDSzuToAj1ADCrvsssuPffjz08//XTsc7DZINpU+DzM4Xz0ox8VZe0f//jHE/eDJCpqrz355JNQBDQaja7CtEirGWDfydijJznYoRTLv8lJDpFdFVehRdEMsL0dPpITty0TJKesmZy4MnIuuypuu1x5nHSSY96u4lZykrYvMnqas+x07CrM8fnEE+aLruRMD3GvHGklpwjwfV/kcXBIKCo3OE4C/zif/vSnRUfmOGA4GW9FVXGwv8jEot5GdlqZHIfCrmqSW1VxXZN1h3MmVWhhRgwb+VGgXxHC2TqU5Z1IciamFpM3AkzbnplMDs17ZqBdxaXkxGyXo7Iqddv4Wfa9obOrvMATRQiT9cmFKo6BKp8ktSgtDGw0k9MwU12FWLN6eMvIcyc509PTgqhs2rSp5378edddd419DlZUYRYnak0dcMABQvlB+2t8PCa4V1CEeZyJqSXCpqOzq+gzOVQkZ0HwWHM4Z4j+vhdCfSE6WfYrQpQqTnt7fQRtSOwqnuqq2VwzOcaVHOYAcKJdxb3fDpFaSHLMKBhJSo6pPM4gu8ohUrhlysiHkeTkblchIUElJtpMEJUa/BlzN3HAaee/+93vxONC4IBQJD9lIjg95eMaVlV88JjCrmq1S9EJ8zicmZx+u4oqjxNnVyX1tiGzq1xekoOfHSMkh6VPToO92/H8duMyOZwkxzUbOs7RrkrK/ZgIHbf3M5s7ycEKSjEMNHe7apX4akkOE7B8/LrrroMbb7wRfvOb38A73/lO2LlzZ7fa6owzzugJJuO/Y3XVe9/7XkFusBLrk5/8pAgil7d8XH0COVt1Vd92qJScBdVVTHYVVWVV3LaolRynaVbJobaqkraZNoahrNVVDU67yslDycmR5MT8biZCx4ig1YIg5j3rGBjp0Lu/GCu51QTfM3ccM0Os5CjbVUg8cMTDf//v/x322WcfrYM45ZRTYMuWLXDRRRcJywlL2G6//fZuGHn9+vU9Vs7atWvhu9/9Lpx77rlw8MEHiz45SHiwuqq8c6v0lBwsafX9AKrhEEkKJUdspwEwPsVKcqjsqsDFpoIOVDptBby4URcjalf5fcfPQXLilCEOJccz2ScntrrKcCaHW8nx3XbuZ2zeNsKsDEWRgYqSExgiOaFlVetreWJSyQktq6llvQ6EY1DFiZKc5557DmZnZ0WrGBh1koO2EJZln3XWWYJkHH/88fCKV7xCfH3xi18svb1zzjlH3OJw5513LrgPrayf/vSnUHbMdzvWU3JCNWdROOyNTMmZX3zJSI7LY1eFqlCtQ3IolZx+64uqR05ydRV18Ljv+ImtKkTYEDQ6pJNDycEwM86wqkSGvBoNHjNmcnDbCwadGsjGCCIVITkmVJy8lZxwX3mTnDglx6RVhVgyNQWTExPQaDaFmqOyhg+dXfWlL31JWEVYin3FFVfAkiVLRBn3/vvvDy94wQtoj3KIMd8IUE/JWRA+JsrPlMmu6leFSCaQh9v1sVJrfnsesRJiuoScQ8nBxTlKnkIywrCjBbkcLrsKx0U0+rbNWV2FOfkFJIpzpEPCbCxjJMfNn+TkWV2V1JPHZI8cBJLqYbWstDM5K1euhNWrV4uv2Om4Xq8PXVvoojcCjM3lUDUPi5AcquAxbic6CI7KrhLbilQkUAaP+0kTvZLTp25Rd9R1e8dEcJCc/u1ydDue33bv+5tLyUHM9m2bU8kR++snUZzDORMsMVMkJ24/RklOTIWVebsqTskx9xoM+zRyZZJzwQUXwNFHHy0IDo5ewKZ2+BUzNQ888ADtUQ4x5jM5BCQnuqgXOHiM051bfovNrqJuBNjdXmTxoc/k9Ck5hFab2F7QS5yMkByGPE4cgcLBnKi4cKFfJeLM5Ij9LSA5BghHn1pkonw8aT9pnYipEUeoTHU77u4vjuQYVnIQw6rkKGdyPvWpTwnFBpvv4fRxnAJukV8JOZ9dRZ/JCbc1UZtgtauolZzo9sirq5iDx91tdjJbHJmcfpLjM/TICeH5s+xWVVyoGcP9C3pSUe+vn+RwB49ztKtyV3L69uV5Pvl5Q6VXjulMDsKSnD6gWvOjH/1IhIIxi4NB5DB8jDdLesxNII+1qwqs5EStL7StsJSTClFViDKTI7YXUVc4g8fUE8i723VMKzn03Y7jCBSnVdW//XYwmHV3C3vlmCA5ffvgnlsV3U80aI3VkXFl3absqrQOxFyIs8dMNgLs75UzbPOrlEnOIYccIm7vec97xM8///nP4W/+5m9ErxpsNIZzpSwM21XhFSaelD2HXMmhyuSEJapdUkK4aviRhoD0mZz597TDWELOoeK0t+sZJTmsSk7ErppjnjMUVXK48zhiH1ElBz/LVBcsBczkIMHBfU2NTRlXceL2ZzqPE+6zv6LOKjkFIDn4R0E1B5UcvN19992wbds20bcGFR2LHILH4eLoERIHbiWH+Iol2hCQXMlhtKvCWVg4iZyP5ESOn8muim6XY6RD3HRzdiUnSnKY8zgL7CokGwbmOPWXqZvK5Ih9ubP5kZwCKDmi63HTg/HJei5zq4Z9fpUyyVm1ahXs2LFDqDlIas4++2w47rjjRIWVhUoJOaFdRXnlx5jJoc7jiO1FSFM0KEwBj9Guam+zZYzkGFFyYpr2cSg5XI0Au9vvs6u4MRvdh4ny8RyVHLEvDFYvgnxIzlxDXGBUOs1m81Bywv2GJMdpto8pLyVn48aNwomJzoYcSZLz1a9+VZCaZcv0A7OjDNoSco+223GksgOVO4fKAouQnP5RDJSZnGhfGwq4jH1yopPIqRsB5mVXxXUmpkK0cos7eBzdPndlldhfNJNjonw8Z5ITzf+YJjmi59LsHNSWLBY/tgxXVvUoSMvzs6oQK5cvh2qlIgjO5s2bxSzIkSY5r33ta4UKgaFjnDeFOPDAA0UH5OXLO38ti1TgmwktPoqxDtxKDpISLP0mJznUdlWH5IjmfeTB4/YC54rOvvQLa6gOsSk5USWq5NVV0W03mK96Z/O0q0wpOXgBg7famLiYcXwnJ5JjfoEXox06JCcPu6q/wiqP8nFEvVaDVStWwNbnnhPh42EhOcp9cn72s5+JmVUYNsZhmXjD7/G++++/n/YohxQhwSHvk0Op5HQIE6VVxWtXNReoLlQISROHVdXebienxHDs0e0iwcHigFI3A/TN2VW91VWuWZJjMBsDrR1GK6t67CrDE8ijiO4zN7sqQq7yUnKGNXysTHJwOObrX/96eOKJJ+CWW24Rt8cffxxe97rXwfve9z7aoxxyq6o+PgH18d4BbSpoenxKDmVlVY+SQ9jtOLq9aH6GOnhMHTrun0TOncnhsqpMNgOMKjn9HYl5lRz+rAQO23XDfZooH+/ueNa4VZW7XdW3z7xITiuyX9PDOYed5NR1lJzrrrtOjHHobqxehw9+8INwxBFHUB3fUINybhWiGWYuKJUcnFDsuWxKDmW3Y4Tf2R51t+Meu4qJ5LDbVZ3XhJPk9MyuYlRyfN8Rt2p1jD2T4wbYoduH8Wp1YQ8bxvDxslrVrJLT2ZdpJSd3khOpsIprzGcCTgHsKsTMEFZYKSs5GDhev379gvtxYOfSpfrWy0h1O57Sn0DOlsnpbI9cyekQsZCUkAEb6bVa5OXjvUoOs13FrORw5XEQaIMhifJ9nE/Ge1UcVm9xk5yomtMwUF0l9hNaVkaVnB35KDkRIhfkqORgji/a8NMkogpSHnOrhlnJUSY5p5xyiggZ33zzzYLY4O2mm26Ct73tbXDqqafSHuWwdzsmUnIczxdt50mVHITbLI2SExInjtbsvtt+ffmUnJDkMFVXdU7gnEpOSKI48zjRPjxY9ddg6A6dlMsxUV3Vsx+jmZyOXWViVlaSkpODVRMqOQ6+5iZ6Eg0Y0pmrkrNq+JQcZbvqM5/5jOjQeMYZZ3SvDMfGxuCd73ynmGtlIaHkEISOEfj5bHk+TJKTnDlyktNtBshAcnCbbks/4xQHJE+uU1K7yg/EyAhukoPbr1bnjJAczONQVv2lKTmCUBmaa9QlOaaqq8S+duZiV6Gq62O1YgtHOri5KTl5VVaFVjh+PivVSiGCx08N0WgHZZKDs6o++9nPwuWXXw6PPvqouA8rq6am2p0rLcw2AoxWWE2S21X0Sk6jc6UfHcNABdym5/I0ssKsD7+Sw7eQ4rZNkJx6nV/JQbuqwTzSIQRaYkhwfENX+t0qLlN9ciL7Mm1XIUnFfU7M8nyuMis5OYWOu12PWx6MTdRymVsVYnoI7SplkhMCSc1BBx1EczQjO5yTrqGiaAjIkMmhJjl45Ya9OFiUnBYqOZPAAfTsy05yODM5IcmZmDBjV3GPdIgqOSa6HfcoOZ6Lb4qhr64KcznjOZGcoKOc5ElyxP4bHlTAZenBJWtX4TQDbHEyDM1+pUjOeeedl/mxV111lcrxjKhdRRM87oaPSxA8RuA2qTsehw0GXWcpn13FFjxusk0gN63kcM6tipaomwgdI+b8wFhlVZfkmMzjRO0q0/vtWGTL8iI5rifGweRVWRUC9x8EPOeWrJianIQlU1OwY3ZWqDkjR3JwIGcU2PQPrwr3228/8fMjjzwi5l0cfvjhtEc5pKAuIZ9XcoofPA4DjmzBY6aGepx2lYMkhznzgaFmIySHsUdOtFeOOSXHg4bBTKrorGyysipHuyokVn5OJCe0rJxmPqHjnmnkBj43WSyrkOQccMABMFIk54c//GGPUoOl4jfeeCOsXNn28Z577jk488wzxUwrC7Nzq0KIYCSHkuMzKDmNnTDGoFqI4LHvlc6uwqF87hzvlZwpu4qzR06PXWVSyTFtVzmGsxkuXhy0WFTbLEpOriRndhacBk+xgkyFVWBAAc0SPn7iqT8MTS5HuYQcZ1Zh6DgkOAj8/uMf/7j4NwsJJYeQ5IiGgISDNAUc+kwOojG3HTiA0jNHCTkCt+sxVVch3FlmktMyZFcxDufMw64SmRxD5eNif2iNmbaNggDm5p4xUq0Wr+TM5ktycs/kuLmWjw9rhZVy8BhDSVu2bFlwP963fTvP4jW0mRxCu6rVmqPv9YDB4wr94tuabTcfo0bQbIFb47Kr+JQcsX1uJcc1Q3I4h3OGQLXIlF2FZGrWMbf4i6aDJsvHO5ibe9b4PsV+3TnwmQn+ILuq1VgEeSs5vuHy/VHolaOs5LzhDW8Q1hTOrELGh7dvfvObokHgySefTHuUQ19dRRc8dpsMHxKmTE5rjovkYPCYaQBl0wHPdUpLcvyWy25XiWaAhpQc7uGcIXA/3S7EBuB4AbQa5i8WZ3MiOW27Kj8Vw9852x3bkhdQScqzR86wdj1WVnLWrVsH73//++G0007rXhni7CokOZ/+9Kcpj3EE7Co6JcdlkDs9ZxbcMfqFsTXHE6x05xrgV3gWv+Ys7+LtzrWgBhNs2281+HMPppQcz6CS0wp82NE0W33jzO0A0ymRucZzkAeEXbUzPyWntW1nu1t8ztVVtUoRSM4q8XXkSQ72x7n22msFoYk2A1y8eDHl8Q0tsHsqR/DYRbuKGE3MBozRN9dz2EiOA8DUk7LF0Lywn+RMMJIcxxDJMaHkQODDLKOq1o/nDZaQh0qn6bPpXHMb5IFZZyf4c/kFj1vb58D4ix1T1GCiKnEQrJLTByQ1Bx98MM3RjBDm5ua6Chhlx2MOktNy0FZaXhqS4zFaPq053sXbRxKibCJnJDk8zaC78LyWqNKpVhl/kQ52YruECh8pjGJbw4UKmIPT4LFz0zDb/C+AMf6/Wz/cuVnwvQpUK+b3jWjtaOROcgLfg+bOWajVTb7LkpWcTZs2iTUKxzWNLMm54447xG3z5s1i+nAU119/ve6xDTVCFadSrcLEIrpPl4dlp8TvyaYoZV3OcmLjWHGdBt/VvdNoQpXxHOShJcKYf2xhaJo9X9kCz/fYSY4bADR9HCFhhuRsb7mwrMbMECPwGjuBUdSLxVxzO8AY/Wd9EGoNBxy/AhM1w79wlOTsArnCc3CwsAe1urb2oIUVy5ZCvVYD1/Ng48aN8MIXvhDKDOWz0CWXXAKvetWrBMnZunWr6JETvVlkDR0vFYNOqeBz2FW+gyO4ybfrNnmkWc+rQODxhAidZpN1ULHPnPtwWybsnRb4BgLBzaAGQWDGQnJdH5qGQs7dfTZ35mAb5VMZ2yY55qzH+EZ8+Y1TQHiiT1G+uSAEXpysXrliaCwrreDxDTfcAKeffjrtEY0IqCeQR68GqNEMvDbJqdJeYXgiJE1XWdbdblCBwHWhwnDVjT1ygrEAKjUeOcdjJiEOzkPCgCWnHAWYyeGvVGkEVQgMLYxIchxDwzkRVd/p9GMyaxXMtcxbZFElJy+4fhUCx4HKRD5KUlfJYerULos1q1fBpq3PDAXJqerI3kcffTTt0YwQOBoBCrgN8mZeLUFy6BetNsmhh+dXUc5h2bbvOaxVGAFzGavjY1fVYCiUnIZfA9+QkuO4PrgGq2/G/DlwDCtHiLkc5laFJMfNU8nxqhAw9r/KrOTk8DePw0ynye9Ik5y3ve1t8PWvf532aEbRriIMHSNqGPokPhl3lRxi+E2eEG+o5LCRHEaS4DOTHNd3WQeAttFckNFjU3IM2lUmlZwxD0lOYLziU8ytMvh7hqg1HXCoO7UrKDl5AvOUfgHsKsT0EJWRK/sPjUYDvvjFL8L3v/99UV3Vn8C2U8gz2lVTxCTHb4HrBVAnzHw2wWchOdiZ2F/kk1dUeD5fJoeV5ATticiccNGuEleLNfbgsZFMjim7ysmD5Ji9qp8LPAgC/Kw7ALXxHOyqfCqrikNyimNXzQxRGbkyyfnP//xPOPTQQ8X3v/zlL3v+jTJIO6zgmECOqPtNeiUHCQ4ujsSotlzwJj16koNX+FxKjot2Fc+JqOJXwGciZyFcI3YVEkH+k/VcYNaucozaVQ1wjZOczmvpjSbJqeZNctxG4UjOU0Mwv0qZ5EQnklvIY74RIG3wto5KDjXJYVJyqo4nrvjHqrThSj+ooL8AHPA8h40kIMlBywCrPLC1AJddhX9OXphRchr+8NtV+F/FUHee2S7JwWzK4pEhOXi6xAujipHKw3Qlpyh21cwQ2VX5UecRB5eSUwtQyfEZgse0i0nF9aDi+eAyLFJCei6hXYUkp70Pt+RKjqkScrPVVa5BklP350T5ALUqm4a58DOeQzYmz0wOni8QRbCr8HPDn5kbjJnV8yQHL7xGRsk577zz4LLLLhNdjvH7NNhMjvm5VfNKzhhD8NgjV3EQHFf8PnfwmOmKc57keFBjqhxGJadC/P7ITckxaVc5niA5ppSVMa8dyneI83WZ7SqT8ANBctycrrmxsipvkoPqLXYKR2CvnPp4vpGP6U6fHMzeYt+7VR1lZ+hJzgMPPNAdRYDfJ8FmcnKyq4IAaj6WYtYLX10VkhwOJQeDx8AQ4MUrmraSM85OcjiAxMMPAqgOiZLTrq7y2tVAzOcct2MjYC5nnLXH0HwJudif58Mihrlxg+0qc6ihTYSWEfjgBR7UKua6SvcqOa1cy8fDqjYsI6/nbLJMjI/D8qVL4PntO4SaMzIkJ5rDsZmc4tlVtQA/pAG5xM1hV2HoGMFxxS/8dQbLJ8DXAIkOk5xcDaqsJAetqvYODJAcIyXkuBgGIpdTqfCqU06HNKOaYyKSO6/k+ObtKsP9ajCPEwK7HtcMjs5AuH57f4HjiguZPC7SxTie8PvChI9XCZKD4eODDjoIygplunjGGWfAV77yFXjsscdoj2hEwDGBHK0qBHnwmJPkECs5+Ku37SoG8tSZeM2fyWEaSdH5G3JmclBZCcAFP/DFjRPNzhW4HzhGMjkIU+HjMb/dlM9kr5zZ8HX0ciQ5OeRyQiUnJDq5KTkdFCV8PD0kZeTKJGd8fBwuv/xy2HfffWHt2rXw5je/Gb70pS/Bb3/7W9ojHFJwdDyu+e0Pikd4YnSCtsXBZlcRKzmi27HYMEPzws4JOGC6uuYOHptRcuYlf07LqhlUwOtkY4TCxgi0D0L1zkgZeRBAPQ8lBy9m8rCr+pScfElOPpZVdBxPYZSclSNOcpDQPPLII/Dkk0/CFVdcAUuWLIErr7wS9t9/f3jBC15Ae5RDCF4lh+5D0gyvxokXklrXrnLJux0jAo+P5HDZVdxKjigf74Qc+TC/SHCGj5sdiwHBHT4OVRxTSg5erFQ6o1lM9srJq7oqSnLc3ElOThVeBbSr1kQqrMoM7XTTypUrYfXq1eLrihUroF6vw8zMDM3RDSlc14UdO9qD8BYRjnXAbscIykyOsKoQbMFjaiWnQ3JctE2YSE5J7SrR7VjsoN1dudxKTmRhYrarsBGgSZKDjQC7+zNoV+VVXYWVVXkqOdH+PEFOvXKsXVVAknPBBReIAZ1IcD784Q+LUjP8+vTTT6dWXlnMqzjUdhV2O6bO5IjQMQIVHcITfJVNyYm8pYnJQteu8jHsCiVUcua3y6fmmFFysNtxiIC5XL1HyTFgV2EjwO7+DCk5eDHjhaotfmXs1dSP2lyU5LgjaldFlJyiDOlcNRxKjnKt8ac+9Smh2Fx88cVw8sknw0te8hLaIxsBkjM+sQhq9THi6iomJSfM5RA1cKl2An7USo7odtwB9sqpEFZqhCSn/b0PNeIGJpXAjF0lgApBrbxKDpaPG1NynPnfw0RDwLB83KSS07WqQoj5VbStKIoYPMY/Z6j+5mlXFTKTs2o4Rjson6VRrfnIRz4C9957LxxzzDGwxx57wGmnnSaGdmJWx8L8BPIwk8NCchCEV1nzzQCJlZzoVRlx+Nh3owt4UG4lh23xnF8kbCanPEpOt0dOCIPh4zztKtEdHQpActxGz2eTK/enQnKeeeYZ4dSMHMk55JBD4D3veQ/ccsstsGXLFrjttttExdW73/1uOOCAA2iPclgnkBNaVT3VVWinEIUuunYVgnDRCu0q8kxORMnhsqvE99QnoSCq5PAs2k50u2wn0RyUHOaF0bhdFcnkiM+yAWLVzeN0d+yMRHVV1KrKi+Sg3YqDf6MoQi5n2ZIlMD7WVu43bNgAZYWyHokfPFRz7rzzTnG7++67Ydu2bXDwwQfD8ccfT3uUQwa+CeTtBSaceVMn6Mza6CE5dCeAaqu9XR9nKRG2yudUcnA4J5eSE6o45VdyIgFKxiqudiNAM0qOk6OSE4RdlmsVs3aVKSXH97sXPHlUV/WTHMiB5ERDx/P3+VAfz7frcaVSEWrOHzZtFrmcvffeG0aK5GCbZ6wQQkUHSc3ZZ58Nxx13nKiwsshYPj5FO4E8rK4Kw8cUkRE2JaeTyWkTMg/q1Tq5kkPdENA3RHLak8g9qFRrfJkcA0oOp13Vm8lhtqt6Mjlij/gXM5LJCS2r8VrVsF3lGFdxENiTC9+nVOeDQXAKoOREy8eLl8tZ1SU5ZYXyO+mrX/2qIDXLltGqEaMAPiUnEl4jWsS6fXIQhPmZsE9OOL+qrv5W7EE0REhdIRKVlKknBUdJTndIJyPJMZHJ4S0hrxkLHkftKiQ4qOaMMbb+DxsBdvfvDa9d1U9yQsvKFMlZYFehPYgFC3Uz++8PHYfgalOhmsspM8lRvjx47WtfawmONskhzuRETvZUJ8ZeJYeINPgBVCILB+UVf7SEnDx43FddxU1yKIHqUE936ZIrOXNRWzLwxc2EXSV+Zr7IjlNyzFdXtXIlObnZVTmoOUl2VZF65TxV4gqrfE2/EQVHt2NEvS+wWNTqqtCqCkE5ibynHLSkdlV7+7TH3j8+g2s0hbngca/KxWlZ9So5mMvxjWVyTJGc3OyqSGXVyJIcp7gkZ2YIeuVYkpPr3Cqe4DHlaAcWkhOxqjiVHMrqKszIRJvOUVdXhRPI+UhO/5U6vZIjFJWoXcU0pBMPvdX3enFVWOHfub85G2vhS+D3fI5N9copkl3l5k1yWq3cGgF27ytAdRWMul1lQaHkEAaPgwBqkZMDlZLDYVeFPXJCUE4i728GyFFZVU4lp/e14MnkLFwcONQcHOnQf/RcFVb9Kg53GXlbxQnMKzk5VVfF2lXeqCk5C0mOXxAlZ80QzK+yJGdIgsftbsfzJ0eq0Q69JeQeeeiYUsnBX5mL5EStKjMkh3bRXjDtnWWhXrgwcuRy+q0qTrsq2u3YRBl5tEeOKZKD9pvT368K/27M4zIQNpNTjkzOhg0bWFtCcEIqQn7eeedlfuxVV12lcjwjAc4J5JRKDoZVHTHNkVfJoWoIGO2R076D7iS9oFmXgeoqtkaAJVdyouXjIQLfoJLDSnLmFh4Ds121wKoKgcSeuMKvSCTH9Ss9F0V5DemMU3KwygsvpKrM/ZEGYXrFCtEvx3Ec2Lp1K6xZswaGmuRkHbyJL4oKrrnmGvj0pz8thnxi/53Pfe5z8PKXv3zg82666SY49dRT4c/+7M/g29/+NpQnk0M5gbz3aoBCyWkBVq0EBjI5Ln23Y2YlJyh58Li9kwCAoGGkaSUnWj4ewmcqI3diwuuUA3AHlY+bUHIWWFUhBNmYzCF47Oam4phWcrAZapJqi1mwKuHsPRXU63VYuWwZPPv886LCauhJzg9/+EO2A7n55puFUrRu3To48sgj4eqrr4YTTzwRHn744dQX9oknnoD3v//9omfPKNtVHEpOT+hYbJSH5MQuwARKTuB5ZG3bFthVnUnkVO1SwpEOxoLHHaJWYSY5HBJ3tHyc265yYwKgppUc7lESyUoOfy6nPtfKLZNTBJITp+L0jHYYh0KEj599/nmRyznssMOgbChMJgftLeyafOaZZ8KBBx4oyM7U1BRcf/31ic/xPA/e9KY3wSWXXFKaltOojHDYVdFuxwiX4OqvJ3RcguBxv5IjQKTm9JMcZDiUDQFzU3JKaVfFZHKYLA6n773KTnL6ysfDCxaPcZ8Lyse7O+Zd7Cuu19MvK4QbOGSz94pPcprJ/1awXM4fSho+1m7r+Otf/xrWr18Prb6yu9e//vWZt4HPve++++D888/v3letVuGEE06Ae+65J/F5l156qVB5zjrrLPi3f/u31H00m01xC4FztvLAzp07BTmjtqui3Y6plJye0DECy4Hxyrxapc3kkCk5Mf46UffSBSSnEz6m8sxNV1d1lRzSvZiyqxa+/4amuiomeCz26flQq9fMNAI0RHLi8jgI5HP4fh2rtodDGic5qAD7PlQ0z3Oy08eLSnLWlLzCSvns/9hjj8Eb3vAG+MUvfiEyOGF2I8zjhAt5FmCgCR+/yy679NyPPz/00EOxz8GBoF/+8pfhwQcfzLSPyy+/XCg+eSNUcarVGoxPLiKuriLO5MQFgvGEWB2nzeSQKTkxNgYRWegPHov7qBa7yARyruqquJxDMExKjkGS44pzHc/8qjglR+wTM2D14bKrkkhOGD7Oi+SEvXIqk7x5pMFKTlF65awqNclRpqrvfe974UUvehFs3rxZ2Eq/+tWv4K677oIjjjhCTCXnxPbt2+H000+H6667DqanpzM9B1UiJBjh7cknn4S8RzqoBrRzy+QQWVYLOh4zKjlUdpXn8S3g/SoOi5ITtz1yAmKohDw2k8MVPF74GqGVwrX+JJEczvBxXnZVXOjYZC4nleQYsqzSMznFUHKmR9WuQhvpBz/4gSAZaC3h7dhjjxWKyXve857MlVgI3EatVoNNmzb13I8/77rrrgse/+ijj4rA8UknnbQg4IhpcAwr77PPPj3PmZiYELe8wTXSob+6CkkO/lfRuNrkIjm1Fk8mx49TcrgyOYS9cuJIDvUkcnfIlRwffz8GcSVOyQlzOXWGIZ1xwWNukpNsV+Wr5HDD8QpAcmJ65Mz/m1+orsdPlXR+lbKSg/bS0qVLuyQFmwUh9txzT0EyZDA+Pg6HH3443HHHHT2kBX8+6qijFjx+//33FzYZWlXhDTNAr3zlK8X3a9euhVEe6YAICNScVlxLfgolp8+uQvXfI+iVExc8JrOr4kiOz0dy2vukU0Fi1TLy/iv5ZXKEvsJgWSWSHKZcTnT+nKleOYl2FTPRSCM5JkY7FEPJsXZVYZWcl770pfDzn/9cWFZY8n3FFVcIsvLFL35RqdIJy8ff8pa3CLsLe+NgCTmGdLHaCnHGGWfAHnvsIZSiyclJsf8oVqxY0T2ukZxAHjM1GHM5dY3sXCPu5EdiV8X1HnGFmldUuyqO5FD1ykkjOTWCWILnx8+QGiYlB4Ekp0IYXBGT2xMWmnYuhxaVwO0ZzTL0dlWqkuOOBMlx0+wqA+M8ZJQcdCFwTV68eDGUCcpnhAsvvFD8wmGV0+te9zrRq2b16tWiOZ8sTjnlFNiyZQtcdNFFohngoYceCrfffns3jIwVXGiJlR1sE8j7gsdQVCUnCBZkchA0Sk61nEpOwKvkxFlVAoTqQP9wzhBIrpAsUOXPkHM4Sa9X4ECVsHldO48T/xpxlJGPxTQC7O6PScnxgyC+wED8I1qAPkClOnR2FX50484XRbKritL1ePGiRbBochLmGg2h5rzkJS+BkSA52KgvxL777iuqoJ599llYuXKl8gntnHPOEbc4DAoz33DDDVAGcCk5da9JLnHHZ3I8fRUn5rAobI34EnL97eIiHVftRJXJqSZcUdKRnPjt0Co5yYsC/m3rtTpbt2Ou0Q5uzNwqTrsqKXTc3p/PZlX1dDVf4CM7APUJ88FjZpKTpuKYGu2AF07BgCpKzOXk3fW4UqkINWf9ho3DT3LQUrrsssuEXDVojpWdXTVAyZkitqvilBzNq83EEnJiqwrhEuQpYq/MCOwqoeLEvJacweOyKTlxVlWP7F7jm1sVgjqTk5THYVNyEkLHYn9M1kWiVRWCk+SkKTnMVlla6NiUkpOWx+k+xvVhbCJfkoMISU4Zw8fSs6twUFf4fRIoS6OHDRwjHZIyOR6LkuOSho7ZlRwCohBnVYn7ia6uk0mOWyIlJ5nkUIaP5/zkEz51Q8C48nHOTE6aksMVPE6srDJQYZWnXTVQyXFauVpV84+x4ePcZldxzrEaBZIzuXgJ3UaDICF47BeP5MTkcVLVhozAC924/r0UJeSJJKckSk7/BPL5HQSEZddNIwHKVLvKpJLDYVclVFblquQwkQ282Kmk/E44tFe3BUYa3BSyjAicto3HecGe1iOnaOHj6RL3yil/krekdhWlklMLHKjEBF2KGDyu9vXIoQoeJ4YICdSQspOcNAJJp+aYUXLiy8d55lc5aZkcBiWnnqLk4J9J96JFqnw8BJNtlKbiIAJmNWeQkiOOIeGCjOwYnPIoOWtGRckZlMOJwmZyzPXJiVNxdEc7eIEPTqmUnHiiQBE8jhvpILZdmuqqlO2IvAyF5+8YuRqdC8zZVeaVnLn04/H0WkIUya4aRHLCXM645giZJLgJ54sFltX4WC5zq+YfU6wy8j8MO8npz+Hcf//94Lou7LfffuLnRx55RPQ6wcZ+FoNKyJewDeekUHKacSoORXUVk5IT1+24O2xPU3ZOUnKQ5OCtUq2UWsmplEjJiRvpwGVXFSmTI47H82FyrGY+eGy4sioEr5Iz+HXkDh9nCR4XZbTDzKiQnGgOB5Ua7Hh84403irJxxHPPPSea92G/HAtzwWMOJSe1d4YGaonBY00lJ0ENaf+jh/M+yEmO+Dc/gFqJSQ5dhVVKdRWhzWLSrkpTcrArEFYv1ggzG2mZHK5eOUW1q7i7HjtZ7CrmMvIyBY+nO3bVxo0bhbCB45PKAmXx88orrxTdh0OCg8DvP/7xj4t/sxiQySFsBsij5PCQnES7SjeTk3aFrxk+TiU5ugtPzARyU9VVpjI5lHZVUrdjFrsqJZPDkcvJouRQo9B2Ve6ZnFb+Sk5BgserViyHWrUqLlj6Z0wOLcnZtm2b6FDcD7wPp4RbLESr1YLZ2VnyjsdxPXIQOkHFRCUHT+waRCfRrmJUcnRJjsdIcpJUHCPVVYaUHNIS8tTqKmw0GRixqzhyOUlzq7r7Y1jw8qquymRXMalIosdhpkwOt101OJODFyG6rUAogAQHiU4ZLStlkvOGN7xBWFO33HKLaBCEt29+85tw1llnwcknn0x7lEOm4tBncuIXGJ0PR6KSg9AhOQlKjnZ1VdqVmSZZSAoei3/z+UhOexK5/sI2TEpOml2FshiVmtOeW5V+3JS5nJrfhOqAzwA1qcLfMZNdxZA/ylPJQRUnSxKNk+R4bguCjOS/MLmc1eWssFI21tatWwfvf//74bTTTus2CESfDknOpz/9acpjHDqSM7FoCmpEbe7DEyR1JqcJKR8sjSvzxI7HBVZy/BTJnlPJaW/fE1dROgtZqpJSMiWnOaiJW4DnIv2KmPZgzsCYXVVPmVvVPSZiJacRtEP5qQiVW4pJsQUiOVnAS3KaEo8tSNfjlUhyHh0dkjM1NQXXXnutIDSPPvqouG+fffYp3YTSspePpyo5WP2j2FCrmUY6tOwqN3FQIA5zrCoOA0yTn3W7HqdlcnSVkMEkx4XamPoCgwQD3wNJoFCKcDhnagm5TzOksxlUwBvwXhbzqwjWAzdD6wFKZWVQ+bjYH7FtMdCqCoHvf3KSMzjvwjWJPCvJAU6SkyGPU7Tw8UxJK6y05QQkNQcffDDN0Qw5uCaQJ1VXBR2iU1eo/klXcnTsquTFAy0rZZKTZmNoKzkumxWTRclh65GDIFk4nVQiRTWks5mh7JfKrhqUx6FWctImkHf3R6zkDLSqQgglc4pux1iV1hy8b6vkFMyuWtUmOWWbX6VNcn7961/D+vXrRag2ite//vW6mx46sE0gT7CroEtyCIPHTEpOaFmNVcfolZwCV1dVUzMm+iRn0NUwTSZn8FU5xZDOtOGcvXaVPgblccRjAtNKDjHJyfo5Jg4Ai3NAhvcdEmMddVenfLzbB8t1ocJQLp0ldNx9bMFIzh9GRcl57LHHRPj4F7/4hZChQ283lKQ9osqQYQKXklNPqK4KczkT5MFjtb9txfWgknJy0wkfpyk5OnaVj3ZPyu9rIpOjg4FZJ88QySGwxdLKx3vsKkMkh9SuGlA+HnY8ppznlN2uahnP40TVnInaRD52VUfN4SE55bOrpks62kGZIr/3ve+FF73oRbB582aRz/nVr34Fd911FxxxxBFw55130h7l0GVyzNhVOr1yUkmO4pVdUvk4RfjYZwoe+276CZ6zuoqiV84gu8qUkkMRPs6i5JDZVQN65IjHGFZyAuJp5JntKmLbSJbkUEOW5OQ10qFoSs6aSHXVwMD6MJCce+65By699FKYnp6GarUqbscee6xoEPie97yH9iiHbgK5ObtKtcIq3a7yyK0qHZIzsO+FxvyqNKtqKJQcPHxtojN4IaAoI8+SyaEa7eAWMJNDbVnN5mRXWZIjp+QUpSHgTMeu2rlzp+iTN/QkB+0oHOuAQKKzYcMG8f2ee+4JDz/8MN0RDhE4JpAPVHIUF2EsL02EIhlJ6pHT3ayiXYUqTlrfi0BDDSk9ycnwfH01p1kYJYdqtEOW4LFOiwYVJYe6wkoueJwPyXG9nElOX9Y0n0xOwNGqSBqTExOwdPFU6SwrZZLz0pe+FH7+85+L74888ki44oor4Mc//rFQd/bee2/KYxwacAWP00iOatdjjuBxWmWVjpKTOrdKcxL5QJKjmTVJGulgKngsoL1wZgweM3Y77u7HoJKDxfnY+sBUnxzdLubqJMcx3u04hFVy2lJ1UdSc6U4up0wVVsok58ILL+z6ckhsHn/8cTGY87bbboO//du/pTzG4QseT9GRnKrvAOoYRjM5imQkaTinbvA4tdsxQieTM+AEjy1idNa53EvISZSc4mRyyOyqDJkcSssqS/CYWskph11FPI/Mr6Tm9/oRDFCflbsdi95SMs8pgJQD5aywUoqNY4djVG6w6zFi3333hYceegieffZZMaRTt+nXsIJjAnlaHkdVUncCP/0KtWxKjkYzukEkBxkOkoRKrVLOTA7CK0cmp5WluiowZ1eJxwVq1Ys9CIKBc6u6+yO8oi+DXUWt5MioOFx2lYxV1dsrpwhdj1eWjuQoKTljY2Pwn//5nwvuX7VqlSU4mUrIlxixqlSVnFSrijF4rK7kZHjPKao5XsrcKu1cTsoEcrrqKndolJy5LJOjAyS0BPO+MpIcilwOXqikqbEcJAeVWi/r64SPI1RUSkVyGOwqmUaA3ecUYEhnWedXKdtVb37zm+HLX/4y7dEMOViUnJQeOaon4dTQMULxxDO4uko1eJxh8VNURAYqORokZ5CKY8quAu2Fszh9cigsKyQ4WUtkKeyqrKFjsT+ixS5zI0AGy0oqk+MNIclRUHKKUkY+Myp2FcJ1Xbj++uvh+9//Phx++OELZlZdddVVFMc3VODok5M0nJNVyQmDKJKq3SC7ymOyq3R65WQiOb6alJyF5ISTyCuKQzrdDEqQieoqXbsKD7GVgcx2K6yq4+wqDhnJyZjHoRzSmdmqilpWY4v0d+zjSAeJ6qq8lRzP0/r8JWVyZFGc0Q7lCx4rk5xf/vKXcNhhh4nvH3nkkZ5/s5ZV/EIY9hagJDlJwzl1qjFSQ8cak4kHKjlcwWOxcbeUSo7OJHJsiZ/JktBQB9qKx+DXVndIZzOoZjR08Pd2tNILWfM44rG+aSXHN9vtOASRoiIIjsRLhu9ftLFrlVouJKfb9XiCruuy17JKTilIzg9/+EPaIxly7NixoyvZTxLaVQMzOQoL2ECSg1AhOXkqOSW1q3QmkWfpkaOv5OBwzsEnYBxHgO//Wq3GVj5OZldlrKwSjyVQcrKWj+drV7WM53GiuRzV9w0JycHwMSHJcSW6HRe1umrz5s1iXuX4uLpiagq0k88sBoaOa/UxGBufKHR11UC7CqGQnxnUDBCv2gZNs459Xhbbh1HJUSUJg0LH88egRtCyVqsFWgtn9sVPJ3yMSk5W6M6vKrRd5dP05pnLS8lRITmEuRzHU1NycuuRUzAlZ/nSpVDvEM6NGzdCGWBJTg55HEo7j6O6qpnF4lBYSAb1yVEtI08bzqlNctxi2FUqyNxjxBDJ0cnlNDKMdOjuR1PJKbJdRVXRJW1XEWVjZPI4HBVWqnYVJXyF6io8xxSh63G1Wu3mcspiWVmSU/YJ5AOUnEAhl5NNyZFfSAbZVapX+5lKyBWIApYi+xl+T9UhndWMJ1x1JceEXWVGycnSCJCyuioraKqr5OwLilxOmewqd4hIDubSXIXqqiJ1PZ4pWS7HkpySj3QYVEKuouYMLCFHyJ4ksUIhyzyggEvJkV9gRY+aDIuY6slHJpOjgsyqGL4/goIrOTKZHM1F0clAxikzOTJ2FRnJKZNdRURy8G2e5VzBSXJ8rKxSfM8UJZcz3SE5ZamwsiTHuF1lbjinqrzNoeRksao4lRwVuypLHqfIdpVM3yF1NUeC5Gj0ypHJ5PgGlRwkOSo5Mj2Sk4NdlSvJcXNTcRBBy8m1R07RcjkzVsmxKJJdpaLkZK6ukkC1lW3BlVVyRDV7lgCvAlEQV11ZHucXleRIvJYGSI6eXSVTXaWp5EguJrq5HNlMzqjZVVTBY5XQMbWS4yrkcYraK+cPluRY5D2BXFXJyUZyPPI8jtis5HaxfDyAkio53NVVEjZXoLxwmgoeF7O6SjeXU8EeMBk+w1G4mkoOzqbLpNZGgZ9Lgq7VeQaPlZUcSrtKobKqaHbVjFVyLMwqORkyOV7+dtWgRoCqDQGzThTmJDnKJeQFsquGSckRdlVgjuTo5HJkeuRQKTnSeRxCNSfPTI4rUaHXT3KyjvkYeAzDYFettkqOhaGRDlnGOoBCdVUzQ3M3aZIzoEeOagl5pm7H4oEKWR8JkqNiWRWmhFyrV44ZJUcmkyNaDyou5J6HFXWBMbtK1qoiITmqSpcu2fD9zBc7RVJyEEHG89cgeAqNAAtHclbNKzlU5I8TluQYJjmUSk7VdzJNL5bJ5OCbNk8lR3YSeZZuxwicP4M3DiVHiSRkmECuX101PEqOTMdjhCrJkc3jiOdonOjHlJScwGzoOIRmNkZFxemWXROEj92MFxWxx+C0cmsEqGuLU2N6ZZvkNJtNeOaZZ6DosCTHsF1FOoE8g4ojm8lpQcYJzNJKjpevkqNQRi5DcmSv/rOqOOYyOaonUInXSFGBcIOKuMkA51cp7cs0yclDycnJrlIlOVRqjqpdRZnL8TSCx0VRcsbHxmDF0qWlsawsySlx8DhrYFFGyckUOmasruJSctoPdsm7HXcf6/GRnHASOauSY8KuQjVNISwjZ1XphY9l5laR2FWS5eMUk8iV7aqclBwqkuNo2VVOIZScorhD0yUKH1uSU+LgcZZGgLIkJ3PVhWyfnIyetuyQzkzdjhXDx1JKDiPJaW9fjvyhIiZDKFTC01iqnWU4Z/fx+GiFBXpO4Qq8NHaVgpKDbzVPY59ls6volBwNkoNDOjWBHdR1lBxkOEVRc2YsybEw0QwwS+hY1q7KruT4UiWlmaurpEvIJewqr5x2VftY5BYm2ddRTclRGDSokMuRGemg2xBQxa4yXV2la1nlZlcplI9TjXbAP5HMBRGHkuNpdDsuWq+cNSWqsLIkx3QmZ/ESo+XjbHYVQkJ1yW5X8Sk5wKrkyJ18soaOlZUcycer2GEyeRydrsdNydCxzmgHpUyO4eoq3V45eVVXaSk5mioSqjhZ+mklH4Cbq1VVtF450yVqCGhJjgE0Gg1xow4eZ83kyJSQt7JMIFchOZntKk/KZvFllBxOkuMXz66SgtKiKX9lr2JXqSg5ynaVY7q6as64kjM7gsFjHasK4RNUV+mMdOhuo2B21VMlmF9lSY5BFQcxsWix8eqqQig5GaurAsnwsUzwWNqukggey1YnZZ1AborkqDU0bBmyq2oFt6vaeSOTSk4+dpUzsiQHSOwqAiXHTiKXhiU5JkPHU0ugWlMvY1RWcrwi2FUuy0LIZVdhBgaDgtkfX6xMjiNLJvDwpYmOISVHYYFSnV+lEjzGF081l6PSJ0enV44fBHKf8Z4nYyfpICeS4+ZKckgyOQRKjl8Qu2rG2lUW8eXjtBPIsyo5Qh3JuIDJkZyMjw0CqEr0qJFTcniCxzJW1VDYVUpqTnGVHOUSckU7QCWXU/VdqKraaopX9KjiKHepFeldddum3mjllsnRKR8PPxsqo2HoMzl+oUjOs88+C3NzamqkKViSY7R8nC50LFNCLpPLkRrcl3EhEVZVwLNAy5WQM5KcopWQq3RJllYHipvJ8Q1mcsTzFIjDmD8LqlDtlaNsVYVQJBsV14OKxgLtinYFQX52FYGaozPSoWgkZ+niKZgYHxffb9iwAYoMS3JKWj6OkJlezKPkZCQ5kvNqsl7t47qSdUCnvF3FTHK4q6sUFBN5JcdUCblinxxJ4oFqnEr1lyrJUS0f17GrZnXHIyhmY3SsKgS+vDqjHQpBcgiUnKLYVZVKpTS5HEtySj2BPPuHxi0RyXEzXm1i6FimLDSQUDdETwsJyC6O/HaVQu5CmuQolJArKBBNxQVKVs1RVXEQKhfYqqFjXbtKC4p2lS7J0Q0fDw3J8XAGXzGIzkxJKqwsySnpSAdETeJDn6uSk7GySvZqX7a5F6ddhdVVMhfz3CRHJagpP7/KTCZHZayDSvhYNY+DcGRaL2iGjrWUnJzsKp1GgBS5HBKSo9H1GPtQeZol+CFsrxw5WJJT0gnkw6zkZG0IKBM6bj9BrrpKFpnDxxITyFWPR0naN5HJkVS8mkEFPMUmbrLhY0dygGvPcwOzSo5M7yuSRoC6JGfElRxRPk40eMqWkcvBkhyj3Y6JlRyJ4HFWJUeuGaDHouS4XEoOVkhkXBxklRyx/YwkQVbFMZfJ8Y0Ej2UCpE2NydG+SSVHwUJQbQSIwN2pEJ0y21Wqox1cvyKX3UtAoNH1mKJ8vGjh4zVWybHgLCGv+g5UJAhJlmoM7KHhMCg5NSYlR6bbcYisZaC+wsk8a95EheTITCLHx6rYQjJKjuxwzu7z8FkSi7NKZVV3X5ILusoE8u5zlaqr9EpvVSyr3OwqEiXHzU3F0bWrKBoBFi18PG2VHAvWCeQSVlVWJUe6SVjmTA6P1SLT7Xj+SR6bkpO1wkqF5LS3n+3YVatQ5AKN6id8mfDxnEJllbpdpZPJUVFy9K7uVcLHcyWtrtKxq8hIjo5dRRA6LmqvnKds8NiCI3gsUz6eNZPDRnIyDueUDx4zKjkSIx26z/G5SY7L0+04hJQyoHFVK3F8KuXjqtVVpu2quq9HclR65eRnV7VGm+QMoZIz01FyNm7cqNx6wQQsyTGZyZkiVHIk8jh8Sk7WTI5sCXnWieUVRruKUclRzAeMopLT8g3aVVrN6vC1C4xlclTsKrQxSewqBdWq1nRzq64iIzmep5Bb6xwDZSanIMHj1SuWi345ruvC5s2boaiwJMdkdVUOE8hllBypbscIzARlIDrySk7WieXWriLrdixdQm5GydGzqxxjdhVmjWQvsPUzOXLH2wg89ZEOPR043ZG0q3TUnGEMHtfrdVi1fFnhczmFIjnXXHMN7LXXXjA5OQlHHnkk3HvvvYmPve666+C4446DlStXitsJJ5yQ+vhidDxekhvJYVFyEBlOeHxKDmfwmM+ukp1ALq/kqA5gLJ6SoxM8NmlXyVpWNb8pVThAoeRoW1UhJD8b2EKiQqA+4MWPymgHUpKjGD4exkwOYmZ18SusCkNybr75ZjjvvPPg4osvhvvvvx8OOeQQOPHEExNlsDvvvBNOPfVU+OEPfwj33HMPrF27Fl71qlcV7sVGr3L79u3k1VXywePBH4ymQrVMJpIjqeSI6qAMREdFycnSEBCnj6v0yeEsIZfL5KgvZtnVHA2SY6y6SlLJ0aiuEs+XUEl0rSqxP0niQEdyWsZVHESgaFk5XgGUHMJMDn5GZQcCc2FmpSU5mXHVVVfB2WefDWeeeSYceOCBsG7dOpiamoLrr78+9vFf+9rX4F3vehcceuihsP/++8OXvvQlcfK84447oEjYtm1bVyKmra4qgF2FyLCg1hT6S2SxNFSUnCwNAVVUHJnFm7+6Sr2pXXY1x4xdpdMnRyaTg4uGbtZBiuRoho7F/iSPV3tuVQjJzwcVyVG1rPK2qzDH40uOiSlb+PipAldYFYLktFotuO+++4TlFKJarYqfUaXJgtnZWXAcB1Z1ytr60Ww2BeGI3kyGjutj4zA2PpFLI8D87Sr57WaZX8Wl5KhUVg1DJkcufOwU3q5CRS7I+J7WtarENgwrOVm7mJMrOZJEg2KkQ5lJDsX08aJaVtMl6JVTCJKzdetW8DwPdtlll5778eenn3460zY+9KEPwe67795DlKK4/PLLYfny5d0b2ltlbQSoYlfh+dAbcBLmIDkVJDgK0iqXkhNwKjklr64SyExAmoUvIZfplUNBcmQyOXnYVdqVVTnbVSpdj9vnvZxJDmEep2gkZ8aSHDP41Kc+BTfddBN861vfEqHlOJx//vlCVQlvTz75pOFGgHShY5XgcRY1R43keKSh46xKjijyYCoh95TtqqD0dpWRTE7GxRlfzpbmApXVstKprOpuQ0LJqWtWViFcHAorEcTVbgRYQruKUsVBBC0FJYmwskr2gspUQ8Aik5w6FADT09NQq9Vg06ZNPffjz7vuumvqcz/zmc8IkvP9738fDj744MTHTUxMiFt+lVXUSk5LqXnYeK1qVMmRHemQ9WpfqdsxIkvwmFvJKWjH46JlctCq0j2V4/yqWiGVHP2FL+gQnbFapdDVVaR2leS+KUPH6nbV8Co5a2x1VTaMj4/D4Ycf3hMaDkPERx11VOLzrrjiCrjsssvg9ttvhyOOOAJGaQI5lqBSKzkcwWPZyqqs86tUuh1ntatUlRyUlwaqOQoTyItZXaWRyckY0Na1qqSUHM3KKulMDoGSI2tZDYNdlbuS4ygo6ENsV0137CrMuIZVxEVDIUgOAsvHsffNjTfeCL/5zW/gne98J+zcuVNUWyHOOOMMYTmF+Ou//mv46Ec/KqqvsLcOZnfwtmPHDihkt2PCkQ7KSg6HXTVg4VW1q7iUHM7gsXjuAJKgquIUqboKiUMAGpZYxgGiTYIshdFMjuHqKrFPCdvC2lVAMolctqEiZSPA7jYLUl21eNEimFo0WWg1pzAk55RTThHW00UXXSTKwh988EGh0IRh5PXr14sZGSG+8IUviKqsN77xjbDbbrt1b7iNIoHLruLI5LRUmpMxKTmDGgKqKjmotmB7do7gMTfJyTKJHMmDTlfbbEqOfjlsFjWnoVE+3t1PRvXCvF01q70/sU+pYadU1VWt0lRXuQTvoTiik7dd5RdEySlDLqcQmZwQ55xzjrglNf+L4oknnoAygCt4XA+axVByfC4lx+XJ5HRm0FRqNRaSM6gEW4fkhGpOrVplsaraOzBEcnBxrvOVjxc9eEyRyZEhOWhFu5odlrtAgorvs2qG5SMIiJUcyS7Wmp+3RMtqfCxnJadIJGcl/P4PGwpLcgqj5AwrOErIq76r1BI+TclxAh88FiVHsbpqkF2lc/IaUGHlK05abj+Xm+S4fFaVQSUni11FksnJeOXvEmRysNJpUJuG9gMD6RYQScDgsVGrKkTGC4EqDuYk7M6L7xtf4jxFbVephI85Mjl4MVWUCqvplcUuI7ckp4TBY5XQMcJNsQiUQseZSsh5gse+TpO4AbkcvUxO+glYNXQ8v32PrRFgewe+OSXHAMkxaVdlVXPqIo8TGFVyyELHkiSH0qpSsawcFrsq+/59JGUaynAZppGvKXiFlSU5poLHlCMdJLsdh/BSmL+SVcVYQj5QydGyqwYQKI2rXt+AXTUySg7BVXjW+VVOhkB6pu1kUC0oGgHK7I80jxMio9pJaVWpkBwWJUdiSCeHVVXU0Q5/sCRnNDEfPC6CkhPQKzn4vJSrVy4lRzl4nMGu8twi21WDSI7mYoaHHxRDyTFVXYVBbZdowcii5FCVj4v9Zbyaz8uuypPk4J9Cy9YmUHI4rKqi5XJmOsHjos6vsiTHWAk55QRyRSXHZ1ByROthlzyTw9YMMItdpVNdNeDKuqp5ZcmdyUEEAxfO8mRyhF0VZLGqaEhOll45pEpObnZVfkqOm/HziSpOADmTHIa5VUUjOdNWyRltdDM5hH1yVMrHByk5TZ3KCw6SE/ipAUOdK7Q0u0oQHI0S7NIrOWInQUEyORSnp2BghRVVHiezXUXUI0fsL2vwmJrk+MVXcjisqmIpOcWyqzZt2gRuhrE5pmFJDiNQBmfJ5CjaVV5KqLSpcxJMuSpXtasQXoq6pDN0L21+lW5IMMiZ5GiXkGfK5Rjqk0Og5CAGkRyK8vHutgwrOajO+hn2OZuXXdXUf6+ovse5SA5IkRzOTE4xlJyVy5eLthb4mc46UNskLMlhRKPREA0LqZsB1suk5GiQnDRVQq+E3GOprMoUPGavriII0A5UIxx2u8oNKuJGAZxflbov0ySHMJMj9plBFcsveOyOtpLD0AhwftvFIDm1arXQZeSW5BiwqiqVCowvmiLbbk21uoojeIxIICMV14OKRpkjm5KTQgSU51ZlLSEvgV01WMnRP3EPep1orKps4WOKuVVZG24i6kSNALv7zGBZ5VZCnqNd5XCRHD9IVYNHza4qei7Hkhwj3Y6XQjWlS62pTA6ef5OalSkHj8WGXXIVJ+1qH38PXyt4zGdX5Z/JMaHk8AePKUY6hLBKTo7VVRx9ciSCx1zIqua4nHZVQfrkRHM5RaywsiSnhBPIdbqlJqk5eiTHYyE5bsLVp3ZZaApR0M7k4FVeQD+BPHt1Fa+Sozucc347OLHdN6PkDHhvlzmTI/Y5gFjjOActpTYO+D4blKvyAxaS4+ZsV8mQHE67qkhdj2eskjOa4Agd62RyEG4C+9cjOQ5pZdWgq30dq4pbyWlvI2BRcdrb9tL7vZAoOWmLF12QNE3NoQodZxntYL66ipjkDCAb5HmcEAPOQ4LgMKzBWHmZZmUXheTguWRQ49FhyeXMFHhIpyU5JZtbpdMM0LiSo0lyklQJHasKgSeegJPkJCw6FCQnbRI5CcEZmMlxjMjttErOAPWLMJPjQ3q1UyVwoSY5SVvXriKvrAox4LPCoeLI5HJYSU6Grsc6TUXLlsuZtkrOaIJtArmOkpNAclpFzOQkHJNWt2MEvgQJiohudRW3ktPevsfXI0fsIMhfySHN5JjrkzPIsqKaPt6zvwGWBZuSM4jkzDGSnAwXI3krOZzl40XL5awp8PwqS3KMjHQgVnIUq6tSlRygLyFXnVsVwkvYrk6340GWFYWSE7CTHJeV5KQrOXQkJy2TQzHSIev8KspMzkCSQ9gIcJAFzU9yMthVOSk5rl/RVnzTELSKQXKKY1et7AaPUW0uEizJKWHwWLW6KknJwTdlEZUcl0vJgRQlh8SuSiA5VH1fEo7doeiRI3aQdpJqGrKrKDM5hpWclNevThw6zqLk5GZXMZSPZw0fc6o4WZUcl7F8vGh21UwnkzM3N9dd94oCS3KMzK2iIzno6Vc1CEmckuNgkkCHfSeRnBFVcspuV5nK5KQHj6tG7CokONRXnmnzq6hDx1kyObnZVY08lZz8SY7PWFlVNCVnYnwcli5eXEjLypKckk0g18njINwYi0ArdJxaQu4yKTmW5LBmcsROAn67ypSSI4Z0BkZUnMGZnDnjdhV5I8Cs1VWcJMcrPsnh7JFTtNEORS4jtyTHUDPAIlhVCC9mAdYnOUkl5B6TkqP/tk3qeqzb8TjNrtKdQD44k0PXCyVZzTEVPKY9NSU1BKTO4wyyqzgyOfinirt4YWsE2N3xCCs5npdY5Wg2k1MMu6rI4WNLckpWQq7TCDApk6NPcnwMPpRKyYEYJQdPWhR9LZIUCm4lh2I4ZxGUHBRCmoRKTnubLnv5eF5KzqDRDrkFjxutoSU5WdQczkaA8/sojpIzbZWc0QNHM0BtJYeD5CBiFlhdJcf3PQhiOtqQKDkxJIfCqsrXrhoOJacFFZ1aP6nwMYdd5eZActIaAuaWyWnyNcIbROiNkJwBvXI451bNH0QQq87nGT5+qmCjHSzJMaLkLC20ktPSmUAeImbBqmkqOaKdTcx2STI5MUTByzhZeeC2c6quoszkBDkqOZQ9copnVzGRnISFDhsTNqhHOnQ3npx1QnVXt/BAR8lxvHyVHNHtmPCiowy5nBmr5IweOErIdYPHHkfwOFHJ0T/Jxc2vIul/4XrDp+RQlZAjEsOs/H1yKEPHA+0qw8HjOkMmJy183MBZY1x9S3C7CZ8ZzjzO/AgTt7B2lYnQcdEsqxlLckYLruvCjh07yJsB6jQCTFZyCBbH/pOdH0CF4MMXq+SQBI9dlm7Hw1JdFafk4KDLAOj24Qe+uHGWjw+yqxzNXk5F6Hgs9pnwnmOrrArh50NyBqk5eZMcI1ZVwcLH0wWdX2VJDhO2bdvW/Z5yrIOuXcWXyendBpVU3a/k4OFTKDkBp5ITN4mcYAK50RLy2EWTfuGKs6w4lJwku4pDyfFwvlhMlqzqO1BlIh1JvXLYKqtCJFi8eZIcPEdQXAhpkRwDoeOi2VVrOiRny5Yt0Gya+/0HwZIc5tDx2MQk1MfGCxM8bp8AAna7SreyKknJIel2jGAMHiPD6VdCqFScpBJytCHTBkPSZHLoT1xxSl2T4So8ya7iyOSI7ca8flx5nDSSw67k5GRXpZEcEyrOoNEOJsrHu/sqCMlZvnQJjI+Nie83btwIRYElOSVqBEiRyYlTc3hIDo0N4PWdpCm6HSfaVVQkJ6ZXDiXJiZtETqriJCo5rdIqOSarq5IsK67KqjS7iq2yapCSwzi3atBoBxOhY4S1q3pRqVRgeuWKwlVYWZJTokaAiJqmXWWK5OgO50xWcohIjov5EkaS4/GRnPb2PVaSE6/ktIwoORyZnKTRDsNDchKUHHa7Kv4zU5/j65EzqOuxMSUn1a4aPSWnqL1yLMkpUSNAKiWnvzsqTQk5j5LTv3iTeu19RKFcJMdl65GTfyaHw65aeOye5yd2p9ZFXMCfo9txd39JmRxrV7GSnKTKNZNKTtp4lLx65ViSM1Ll43ShY4rqKkR/8ygWu4pKyek7NpLy8YSGgFTVVWJbfUSSKnTc3X4fQSPtdtzp/pybkuObsau48jh52VVxYWf+TE78e6LK2AhwMMmhf/8kIUjIHpoMHqNdxdUlYBjKyC3JYe92TK3kNMmvMgtdXcUVPI4hORRzq7rbNm1XUfbIETvodGMclkwO4LTxvtfMMMmpMwaP240zA/PVVbmWkCfNcKP9rKUhcFq5KznIcJLaVpiGVXJGCBzdjhG1AZ0+syB6MsSKHGcEg8fsdhVj8NhEJidezTFUXcVU/tuv5nDMrUqtrmLqkdPdZ8xCl5ddVWecWzU4k2NSyYmxQd2WsW7H8/sshmU1Y5Wc0QFH8LgSuCR9NqJKDomKw2hXmVRyODM5VBPIjWVyEAsWTQYlp49I4VuzxUVy+j47xu0qRiVH7LNPFcO8SB4kp+J6JI1AB8ENnFiLzjGUyUkkOQbLx4vWK2emgPOrLMlhLyGnnEBOs8hElRySbsclVXKiJAcXhHIFj00oOYHx4DGGjrmE9/6GgO4QZXLEPvteS7yAoeydFAssWuj73JiwqsSug/j3vangsUCskmO+EV7RlJwNGzbwjRORhCU57CRnSWEaAcZVV9EpOR5LM8CFSg6PXSWsDMIP5TDYVSaUnH67iiOPk6TkcJKcuOoqrrlVSXYVe+g4RE4kJy58LMZpGczk+DGTyPNQcorSK2e6Q3JarRZs3boVigBLcrjtKsIScorQcb+SQ0Zy8IousmDVWlxKDo9d5RFWVsUqOczVVRx21UIlhz94zJXHicvkOJyZnH7CHATGlRx2qyqhwspEI8CkXA6qOAGYIzkQczFnNHQc7rMgZeRj9TqsWLasULkcS3JK1PGYTslhIDmIyCJClckRV2aRY6S8SovOr6K0qsS2+0463EqOE9PBWRsRohYEPgQcdlXfkE6O8vE87Cq8sI7mRbCJZ4XNiOvss5/kcFdWhfCLo+QYtarEr94qBMkpSiYHsWZ1scLHluRwl5ATVleRKTmRxaspaoUZSA7hdOeoFUPbDHB+u35Cv4/y2FXcSg5ftUxUzeFoBJhH8BgpTvRigrMR4CjbVW7OJCcuk+Ma7HZcNLsKMbOyWGXkluSUKHhcS5ikLAuXI3iMCMlIEJBlchBlVHJ67CrCCeTz2+8NTfNncvhITjSXw5rJ6S8hZ776jVpW3FZVoeyqEVJy4qqr/DzsqgIpOdOdXE5RKqwsyWEALjocJeR0mRyG4DGis1gJFYfwwiJcwPEcTum3RzM5cZO9tbbtz3chpVZx+pWcuF4z9H1yTCk5nHaVuUzOApLDXD4ep+QYs6sKFDw2WT7e/Zz3taJw8yghx47XBRFzZgrWK8eSHAbMzs6C23njF9Gu4lZyqPI4/UqOT21l9NhVfOXRHCQnOomceqSDaSUn2iunwbhIRedXoZ3Y36Nn2JSc3OyqnIPHphFVc/BzmUcJebvrcTHUnJmCdT22JIfRqqpWazA+OVW44DFynLB/BkfwmDKPE1UqSLsd99tVLkflUMBGctrb9/gaARrM5PTaVWaqq7hVHLEPw5kcrJqM9iYZDbvKLRTJEeeRnCSVouRyZqySM/yIWlWVSqVwzQCjak6DheTQnlzn7Srat2vgeV1XjUXJ8blJTvt1YcnjmFRyIlegTW67KjCTxxH7MKzkBH3EatZWVxknObmoON19F0TJWW2VnKEH2wRyIrsq2iuHx64iVnICHiVHoGMrUg7nXKDkcBy3VXIUgEpH++/tRFQ8M0oOP8kR+4wQxjyUHLSqKwZtE893e1oQ5EJyWk6ueZyilZHPdJSc5557TkQ38oYlOazl49QTyCmVHJ/RriqHkoMIQ4M8mRwzdhVLjxwEHn53oS5/8Fjsq5PLMaHkOBGOWGcezhnC7bzncOiuG1n8WYHnkg5RNaniIII+JTMfJaeVa7fjotlVS6amYHJiojBqjiU5JWkEiKgFdAtN2CunRXki7JxsakzBYw4lBy0rhE/c8Vhsm92uCpUcviv2oGtZ8Ss5blARN06ESg7nBPIQTuSzZVrJMWZV9ak5VYOh437LyvUr4DO/f+IQRC7q8mgE2N13QZScSqVSqFyOJTkM4CgfpwweRzM5LCXk1HZVV8nhs6s4lRzqCeQLMzmM1otBJWeO1arqDR/zNgJsI3phbSKTEyU5xqyqEJ3PT92wktNLcvJZzoLI/CqbyentlWNJzrArOQW2qzCTg3K2x6DkkNtVXCXk7HYVXwn5MCo5nCMdcrGrQoIY+KSf3dR9dv5exsrHQ3Q+P6btqmjX49xITjR4nKOS0z8vL0/MrLQkZ6jBETyuBB5UCU9cqOSQho4RHWWBTclhsqsCv30rX3UVb/DYmJLTzYcZUHI673lz1VXhYM7ArJJj2q7qkLhawwyZK5SS01NdlWcmpzhKzpoCVVhZksMZPKYc6UB8JYhKDmn5OKOS083kcJzEXJelsspkdRVb8LgzaBSHcwK4rCQH+7twh44RQWdBNGFX4YBODB+b6JGzIJMzQkpO2BAwN5IjLpT8Qig5Rel6PN1pCFiE0Q6W5JQkeEzV7ThaXUWu5OD2xNwqj6e6ikPJcT0Wq8pkdRXXWIf58LTTM1GbA/g78JaP9wWPDV31omJqKnQs9td5z+WVyak13ZFTckI1By8G8szk4Lm3KGrOTIEyOfW8D2AYwTGBnEPJaVKXmOJlhO+Rj3XA7szYC4MjeCyuwhgqq8zaVYwLS4fkcAPzS2YyOZ3gsYHqKrGfIDBWPl6E6qpclJyCkBy/GuTW7binV85Y/trFjCU5o5LJobOrqIOLeIXZtiGI4bvkSk5oWXkcV/quy6bkiNAuwwTyKMkR5I+xH0o7eMyfs0Alx0gmx3cE+fQ8cyTHpJIThp1zU3JyzOQ4Xo4kp9UCr8rfYLIsvXJmOiRn48aN4Hke1Gr8FzBJyJ/yDXUmhy54XA9oZVDsk0NuV4Ukh1jJCdUKFiWHkeQIJcfj69uBJeRuZBo5z07wpNk0ouTMBWaqq0xZVSHpMFU+Hiq0XhDkk8kJgpzsKrcQSk6ejQBDFMWuWrViBVQrFUFwNm/enOuxWJJTkhLyGoOSQ9ojp4NKqwmVnnEANGhhABZ47CoukiOka8ZzPoZ1WwyDRRcqOc7wKDmBa5TkYIWVyeBxaFkZr67yWlBFgsPw2R+4a7+taOZOcvLM43RgSqEchHqtJohOESwrS3JY7SrK4DF1JsdnITnVJs+iiySHA5xKjoDL24HVYSY57SGdLTOZHBPVVYEHDoOdmmpXGVRyEA3X41Fp0+C7UJvLb5Fveg6PnV0yJccviF0VtazyrrAqFMm55pprYK+99oLJyUk48sgj4d577019/D/90z/B/vvvLx5/0EEHwW233QZ5w3Gc7lAy2uAxsV3FpORUm81SkRxAJYcpeIzgXmtazCWr7eoqM5mchqEr8VZk1tCw2VWI7Yzv5zTU5vIbxtjI2abBIZ15lo8Xza4qUvi4MCTn5ptvhvPOOw8uvvhiuP/+++GQQw6BE088MdHP+8lPfgKnnnoqnHXWWfDAAw/An//5n4vbL3/5SyhCHgcxObWkuEpOgMMQOUgOzwmW67OLdpUXmaJMDuYrqxb3gmaK5HioLNaMXYiYgqiuMmxXbed8P6egttMsmYtiLm+SI+yq/JWcIpGc6U6vHEtyOrjqqqvg7LPPhjPPPBMOPPBAWLduHUxNTcH1118f+/jPfvaz8OpXvxo+8IEPwAEHHACXXXYZHHbYYfD5z38eikByxienoFYfK2wmB9Fg8O257Kpui3wG+C3Gk1PZ7Sqh5vBfoTa8Cpg6PZskOe1MjtnFf0duSk5+i3wrZ5umbVflr+QU0a76Q84kpxAl5K1WC+677z44//zzu/dVq1U44YQT4J577ol9Dt6Pyk8UqPx8+9vfjn18s9kUtxDbtm0DzjzO1JKlsGic7sp0vF6FytgioIRX2QmTddptTgZVqC6aBGr4tXEYq/Bw8moVoD4xybPt2hhUcAdM8CsBjI3Rkek41PD4a7Tvk344lUmYqpm55poNfFg0ZkY1GqsEUK/VAfBmCHO1JvnnOgsmPWD57GeBUx+DsTyv2Ss+VKoB23lEBrWxCaiYH8a+ALuuWSO+WpIDAFu3bhWlZrvsskvP/fjzQw89FPucp59+OvbxeH8cLr/8crjkkkuAG3Nzc7B06VLYY5dpeMfx+xBumXJbFhYyeK2RvfwJDCteZnyPxxvfo4VFL3YsXQW1z39BrO0w6iTHBFAliio/qOSsXbuWfD/HHnus2HY4dNDCwsLCwmLUcPzxxwv3JM9GgIUhOdPT0+KF2LRpU8/9+POuu+4a+xy8X+bxExMT4mYKaLdZWFhYWFiMImo5k5sQhViJx8fH4fDDD4c77rijex8qIfjzUUcdFfscvD/6eMT3vve9xMdbWFhYWFhYjBYKoeQg0Ep6y1veAkcccQS8/OUvh6uvvhp27twpqq0QZ5xxBuyxxx4iW4N473vfK+SwK6+8El772tfCTTfdBD/72c/gi1/8Ys6/iYWFhYWFhUURUBiSc8opp8CWLVvgoosuEuHhQw89FG6//fZuuHj9+vU9FtDRRx8NX//61+HCCy+ECy64AF784heLyqqXvvSlOf4WFhYWFhYWFkVBJcABOCMIDAcvX75c9LVZtoxuxpSFhYWFhYVFMdbvQmRyLCwsLCwsLCyoYUmOhYWFhYWFxVDCkhwLCwsLCwuLoYQlORYWFhYWFhZDCUtyLCwsLCwsLIYSluRYWFhYWFhYDCUsybGwsLCwsLAYSliSY2FhYWFhYTGUsCTHwsLCwsLCYihRmLEOphE2esbOiRYWFhYWFhblQLhuZxnYMLIkZ/v27eLr2rVr8z4UCwsLCwsLC4V1HMc7pGFkZ1f5vg8bNmyApUuXQqVSIWeZSJ6efPJJOxeLEPZ15YN9bflgX1se2Nd1dF/bIAgEwdl99917BnfHYWSVHHxhXvCCF7DuA98cRXyDlB32deWDfW35YF9bHtjXdTRf2+UDFJwQNnhsYWFhYWFhMZSwJMfCwsLCwsJiKGFJDgMmJibg4osvFl8t6GBfVz7Y15YP9rXlgX1d+TAxRK/tyAaPLSwsLCwsLIYbVsmxsLCwsLCwGEpYkmNhYWFhYWExlLAkx8LCwsLCwmIoYUmOhYWFhYWFxVDCkhxFXHPNNbDXXnvB5OQkHHnkkXDvvfemPv6f/umfYP/99xePP+igg+C2224zdqzD+rped911cNxxx8HKlSvF7YQTThj4dxhlyL5nQ9x0002iK/if//mfsx/jqLy2//Vf/wXvfve7YbfddhMVLC95yUvsOYHgdb366qthv/32g0WLFomOveeeey40Gg1jx1sG3HXXXXDSSSeJbsH4uf72t7898Dl33nknHHbYYeK9uu+++8INN9wApQFWV1nI4aabbgrGx8eD66+/PvjVr34VnH322cGKFSuCTZs2xT7+xz/+cVCr1YIrrrgi+PWvfx1ceOGFwdjYWPCLX/zC+LEP0+t62mmnBddcc03wwAMPBL/5zW+Ct771rcHy5cuDp556yvixD9trG+Lxxx8P9thjj+C4444L/uzP/szY8Q7za9tsNoMjjjgieM1rXhPcfffd4jW+8847gwcffND4sQ/T6/q1r30tmJiYEF/xNf3ud78b7LbbbsG5555r/NiLjNtuuy34yEc+Etxyyy1YWR1861vfSn38Y489FkxNTQXnnXeeWL8+97nPifXs9ttvD8oAS3IU8PKXvzx497vf3f3Z87xg9913Dy6//PLYx//v//2/g9e+9rU99x155JHBX/7lX7If6zC/rv1wXTdYunRpcOONNzIe5ei8tvh6Hn300cGXvvSl4C1veYslOUSv7Re+8IVg7733DlqtlsGjHP7XFR/7P/7H/+i5DxfmY445hv1YywrIQHI++MEPBn/0R3/Uc98pp5wSnHjiiUEZYO0qSbRaLbjvvvuENRKdg4U/33PPPbHPwfujj0eceOKJiY8fRai8rv2YnZ0Fx3Fg1apVjEc6Oq/tpZdeCmvWrIGzzjrL0JGOxmv7ne98B4466ihhV+2yyy7w0pe+FD75yU+C53kGj3z4Xtejjz5aPCe0tB577DFhAb7mNa8xdtzDiHtKvn6N7IBOVWzdulWcjPDkFAX+/NBDD8U+5+mnn459PN5vof669uNDH/qQ8Jn7P5CjDpXX9u6774Yvf/nL8OCDDxo6ytF5bXHx/cEPfgBvetObxCL8u9/9Dt71rncJgo5dZi3UXtfTTjtNPO/YY48VU6pd14V3vOMdcMEFFxg66uHE0wnrF04qn5ubE/mnIsMqORZDgU996lMiIPutb31LhBQt1LF9+3Y4/fTTRbB7eno678MZOvi+LxSyL37xi3D44YfDKaecAh/5yEdg3bp1eR9aqYHhWFTErr32Wrj//vvhlltugVtvvRUuu+yyvA/NIkdYJUcSeNKv1WqwadOmnvvx51133TX2OXi/zONHESqva4jPfOYzguR8//vfh4MPPpj5SIf/tX300UfhiSeeEBUY0YUZUa/X4eGHH4Z99tnHwJEP5/sWK6rGxsbE80IccMAB4ooZbZrx8XEYdai8rh/96EcFOX/b294mfsYq1p07d8Lb3/52QSLR7rKQR9L6tWzZssKrOAj7V5cEnoDw6uuOO+7oWQDwZ/TZ44D3Rx+P+N73vpf4+FGEyuuKuOKKK8SV2u233w5HHHGEoaMd7tcWWx384he/EFZVeHv9618Pr3zlK8X3WJprof6+PeaYY4RFFRJHxCOPPCLIjyU46q8rZvL6iUxIJO2IRnWUfv3KO/lc1tJGLFW84YYbREnd29/+dlHa+PTTT4t/P/3004MPf/jDPSXk9Xo9+MxnPiNKnS+++GJbQk7wun7qU58SJabf+MY3go0bN3Zv27dvz/G3GI7Xth+2uorutV2/fr2oAjznnHOChx9+OPiXf/mXYM2aNcHHP/7xHH+L8r+ueF7F1/Uf/uEfRNnzv/7rvwb77LOPqG61mAeeH7HtBt6QAlx11VXi+9///vfi3/E1xde2v4T8Ax/4gFi/sG2HLSEfAWCvgBe+8IVikcVSx5/+9Kfdfzv++OPFohDFP/7jPwYveclLxOOxHO/WW2/N4aiH63Xdc889xYe0/4YnOwv992wUluTQvrY/+clPRBsJXMSxnPwTn/iEKNm3UH9dHccJPvaxjwliMzk5GaxduzZ417veFTz33HM5HX0x8cMf/jD2vBm+lvgVX9v+5xx66KHi74Dv16985StBWVDB/+WtJllYWFhYWFhYUMNmciwsLCwsLCyGEpbkWFhYWFhYWAwlLMmxsLCwsLCwGEpYkmNhYWFhYWExlLAkx8LCwsLCwmIoYUmOhYWFhYWFxVDCkhwLCwsLCwuLoYQlORYWFhYWFhZDCUtyLCwsLCwsLIYSluRYWFiUAq94xSvgfe97H9ljo4/Bxu84rXrVqlVQqVTEIFILC4vyo573AVhYWFjkgVtuuQXGxsbE9zjF/oYbboA777wT9t57b5ienhYk6NBDD4Wrr74670O1sLBQhCU5FhYWxtBqtWB8fByKAFRtQjz66KOw2267wdFHH53rMVlYWNDC2lUWFhZsQDXknHPOEbYQqiMnnngi+L4Pl19+ObzoRS+CRYsWwSGHHALf+MY3ep63c+dOOOOMM2DJkiWCfFx55ZULto3POeigg8Q2Vq9eDSeccIJ4Xgjczwc/+EFBZnbddVf42Mc+tuDY8Lje+ta3wl/91V/B+vXrhVW11157ift+9KMfwWc/+1lxH96eeOKJ2N/xk5/8ZPcx0ZtVgCws8oclORYWFqy48cYbhXrz4x//GNatWycIzv/7f/9PfP+rX/0Kzj33XHjzm98sSEWID3zgA+Lnf/7nf4Z//dd/FTbS/fff3/33jRs3wqmnngp/8Rd/Ab/5zW/Ev5988skiWxPd7+LFi+Hf//3f4YorroBLL70Uvve97y04PiQy+G8veMELxHb/4z/+Q9x31FFHwdlnny3uw9vatWtjfz8kSOFj8IbP2XPPPeGNb3wj+WtpYWEhB2tXWVhYsOLFL36xIBmIZrMplI/vf//7gkQgMANz9913w9/93d/B8ccfDzt27IAvf/nL8NWvfhX+5//8n13CgiQkBJIJ13UFsUFCgUBVJ4qDDz4YLr744u4xfP7zn4c77rgD/uRP/qTnccuXL4elS5dCrVYTik8IJGZTU1M998UBn4s3xEc/+tEuKYser4WFRT6wJMfCwoIVhx9+ePf73/3udzA7O7uAaGBW52Uve1k3H4M/H3nkkd1/R8tpv/326/6MFhcSICQ2aIG96lWvEsrJypUre0hOFGh7bd68Gbhw0UUXwd///d8LgoOWl4WFRf6wJMfCwoIVaBmFQJUGceutt8Iee+zR87iJiYnM20TVBa2nn/zkJ0I5+dznPgcf+chHhDWFWR9EWDkVAnMymNPhACpGaMFZgmNhUSzYTI6FhYUxHHjggYLMYMh333337bmFmZd99tlHEBQkLCGee+45eOSRRxaQlmOOOQYuueQSeOCBB4S99K1vfYvsWHF7nudlIjhop1mCY2FRPFglx8LCwhgwu/L+979fhI1RVTn22GPh+eefF6HkZcuWwVve8hZRUXXWWWeJ8DFWTa1Zs0aoNNXq/DUZEiDM16BNhf+OP2/ZsgUOOOAAsmNFwoLbxaoqPCa0zKLHgPj4xz8OX/jCF+A73/kOTE5OwtNPPy3uR9tMRpmysLDggSU5FhYWRnHZZZfBzMyMqLJ67LHHYMWKFXDYYYfBBRdc0H3Mpz/9aWFtnXTSSYIY/d//+38FGQqBhOiuu+4SZdrbtm0T4WMsM//TP/1TsuNEMoakC9Wnubk5ePzxx3uUGqzkwuPE/Ych6hD33nsv/PEf/zHZsVhYWKihEkRrLi0sLCwsLCwshgQ2k2NhYWFhYWExlLAkx8LCwsLCwmIoYUmOhYWFhYWFxVDCkhwLCwsLCwuLoYQlORYWFhYWFhZDCUtyLCwsLCwsLIYSluRYWFhYWFhYDCUsybGwsLCwsLAYSliSY2FhYWFhYTGUsCTHwsLCwsLCYihhSY6FhYWFhYUFDCP+Pzg2T/mQ8st3AAAAAElFTkSuQmCC", "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": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.13.2" } }, "nbformat": 4, "nbformat_minor": 4 }