{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# More about spherical harmonic normalizations and Parseval's theorem" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## The variance of a single spherical harmonic\n", "\n", "We will here demonstrate the relatioship between a function expressed in spherical harmonics and its variance. To make things simple, we will consider only a single harmonic, and note that the results are easily extended to more complicated functions given that the spherical harmonics are orthogonal. \n", "\n", "We start by initializing a new coefficient class to zero and setting a single coefficient to 1." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import pyshtools as pysh" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "pysh.utils.figstyle(rel_width=0.75)\n", "%config InlineBackend.figure_format = 'retina' # if you are not using a retina display, comment this line" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "lmax = 100\n", "coeffs = pysh.SHCoeffs.from_zeros(lmax)\n", "coeffs.set_coeffs(values=[1], ls=[5], ms=[2])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Given that we will perform some numerical integrations with this function below, we expand it onto a grid appropriate for integration by Gauss-Legendre quadrature:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAABOkAAAKpCAYAAADpKXgkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAACMQAAAjEAFnL93WAADjfklEQVR4nOz9eZzs213X+79XVVfvvc+Q6SRBAhHCEMYgEML4A5UwKIPiZZBJJIAgXEQUELw/Bb2KKCij4IAoIAhBDEQJgwKiV4FIZJQwXQgGCFMSkpyz9+nuGtb9o2onO7ve79716f5WVw+v5+OxH33O2qtWrfrWd/jUd1evd+u9CwAAAAAAAMDujHY9AQAAAAAAAOCq4yYdAAAAAAAAsGPcpAMAAAAAAAB2jJt0O9Zae7i19vCu5wHcC/sqLgL2U1wU7Ku4CNhPcRGwn+KiYF/FJrhJBwAAAAAAAOwYN+kAAAAAAACAHeMmHQAAAAAAALBj3KQDAAAAAAAAdoybdAAAAAAAAMCOcZMOAAAAAAAA2DFu0gEAAAAAAAA7xk06AAAAAAAAYMf2dj2Bq6a1Npb0lnc2rdqftJsZARtjX8VFwH6Ki4J9FRcB+ykuAvZTXBTsq5dfk/SEO/7/V3rv89IAvfdhp4RjtdbeWtIv7HoeAAAAAAAA2Jq36b3/YuUB/LorAAAAAAAAsGPcpAMAAAAAAAB2jDXpzt4r7/yfp/z1z9P4/vt3NRdcFvzWOoCz0nY9AQBXBvUNgLNCfYMBzG/e1Mu+5MvubHpl6ptwk+7svV65Mb7/fo0feGBXc8FlQREL4KxQxAI4K9Q3AM4K9Q22o3wl49ddAQAAAAAAgB3jJh0AAAAAAACwY/y66471JnVulapd1F9nuKjz3pJz/z6e9/lht875rzn0cz6/M3fOtwfvF47D9fJi4f3ChXbOr0dcL+9yQbcH7+PSENuB20MAAAAAAADAjnGTDgAAAAAAANgxbtIBAAAAAAAAO8ZNOgAAAAAAAGDHCI7Ytb2+/HMRbHGafYixBxmjuNKj6V5eXPiM3/5BFj8uv8bTr6C5k3lftOe7bLa1AO0WF7Yt7erFndp2L76WM1/Ud4DnK895iJPFBdz3LupC0+fGNs/XZuxBnm6bdU9qvoDX4tKch5pb4cS11VCKbY1NfXM6V+Fa7ISdPR4D52XeFbuoy85L3UONszTAvR2+SQcAAAAAAADsGDfpAAAAAAAAgB3jJh0AAAAAAACwY9ykAwAAAAAAAHaMm3QAAAAAAADAjpHuumN9slCfLE744GHn8vpjn3GEyhYTUVt6LXaMMHB6PpfYljZd+AvbOkRIT2HOUZrzIGP75tKet9X9pjj2KZ/vXI3t7CJVqfCctVTV7c2jNm4tVdG+xpiIFhLU7jmpezzfcYZ4v9y8h5hHOW1tS/PYpq1GRKbn3NK4u0iOvIB1T6xvKs/ZaxOx3au1QuVcloYYIrF1iATWAeqhC1n37OQYPePnO+dJlWdd98SPN6XnHKDuKV7rih/LaswYpfqmMO6x7YW+tr6pzmOo/qe1zbrnFK/lxPd27sA36QAAAAAAAIAd4yYdAAAAAAAAsGPcpAMAAAAAAAB2jJt0AAAAAAAAwI5xkw4AAAAAAADYMdJdd6ztz9WuzY/vtMXEsVKo1xApYql9kISy0HWQBDDfbJPVYqCLH6Sb+Ji4NeI8Nu9bSmGrbo/Q3W2TrSbEFl57ORiosp2CrYURDTXulk4520wLLSWDxTHCMVpIC02BrZUxBnndo2IO27be81Hx+QrbKSaUDZDuWkqZPeuk2jR2deitpRYX57alGqcYWhoGqaTDHzfOAGO7NPktXr9icmwhITbNry22dcJJ7aHucfVa9RRZeL+GSIit1DfJNuuv8nOe1lWoe8rp5G6MQn1TfM5S3RO+HhTPFfa1+CfM23rztPvq9iilz6axzTYp1TfHtdu+57x2qoydum41mfne82iH97i3swG+SQcAAAAAAADsGDfpAAAAAAAAgB3jJh0AAAAAAACwY9ykAwAAAAAAAHaM4Igd27s21/ja7GQPLqxG2osrl9olEatjVBY/TosIh+e0Y6f5LVKiRCH0Ic5jfewWVxcNY7v5DbGQdiUgQhok3MEuwhz6p4Wjtzm/QRZyNrY6RrKtRZiT7e2SpbGHWIQ5zyMc5+6fs6qLC5sAhTiP9M9nlX9WS8difI2nX+y3FI6RFo92B0L1nxPdaxkgwKLFMTYfuxUXg46LH5txyofGECu5b2s1+AGugduseyr1TZxLMQDA1klpIvE66uqe8HxpEHstTvWXb3bdBwlmiLWkb64EZsX6JrTb9dOL8/OvcfO+aR5bDX2oBnpVnPO6pxZUdcZjpIvJEKEP8Tpv3rC0/6brfKGGa2l+8QK7pWtg4bVI1forzNnWqYW+oX+lBpFCnVStb5KzrntOsX/cMxR0A3yTDgAAAAAAANgxbtIBAAAAAAAAO8ZNOgAAAAAAAGDHuEkHAAAAAAAA7Bg36QAAAAAAAIAdI911x65fm2rv+vREjy2FLcWEskpCrG9fLPy93hg6Zp4zpbtW5t1Tilgau5ByFtOITHvcpDEBzLyWAYKIKimpqX9KEWshtCYml7k0uDR2JSF2kJTZzfvG56ym4J512toAhkhVLY9dSAu1Caxp7OoYhX/OSn3byBzn1edzCbHj9Hxh7DC031AD7JSVNLPUP/WN7WbY4hiufxoj7r6j9ZNLSjMbpbFTglohha0yRjJAuPMgtln3uJqlOsYi1jKmLdVOttXXLKm+iYmt7nrnEuaPGcPtDS7pXjom9dXte9WLTCExt5QaX6xvSrVMte457fNJW6udBklxvQJ1TznVfoAxKon0262dzANiX//m2honjZGOjWo95PpW3q9CbbIco5AQW6hlYs0S201bsX6wtVMcwzZrZGqn4/r7vrute2aHJ7u3cye+SQcAAAAAAADsGDfpAAAAAAAAgB3jJh0AAAAAAACwY9ykAwAAAAAAAHaMm3QAAAAAAADAjpHuumOPuX6oyY0Qz3cPi0oya3EMl1w2D0lkMRgs9J+7NNNiUtp8vj52TFVL8y4kpfV52NauPSSDyaUcSSE1K71jlTHCCKndzDunu4aUs5SK5sYupotV5ldJUCunu5odfpCktGSIBLUBDJLumnbfUmpWGiPFRBXGKLTHvuF0bvvHJNg09vprXIQdJO4ehX+aS/tYHLuScpbe87Hp79p0TELZ2KSqhtcdU8TM2CmBdWyeL40dX3aYx7iQ+joqppm5/tXDPD3naVXqG8nvk5X6JvVPfeepNikktsYaydQ3y/6b1yxpbJd+2MPzpaREW+OEmiAd/zH11fbduOsxxW4Y27yWUbG+GaLeyImyhTGGSJlNKb2uf7U2uUyp9pXT0wCpqjYltTh2To0vjFGpbxReS/rIa+qb1SM2alqOEUao7JNJ2vkKtV0co1L3FOqhmBof6g3XP9c3m9cb41SDpHmEXcHVSZX6Rqoduqepb6aHhyd+7Guf/9QjAAAAAAAAADgVbtIBAAAAAAAAO8ZNOgAAAAAAAGDHuEkHAAAAAAAA7BjBETv2hOs3de3G6xZCXMSVPTe3MMsipgWUU7tb/NiNK0mztPhxIThiOverfabFmWemf3y+uV/4cW4WB46LLYdF6ReV+9xhYWD3FrS0CH5iV2f1XSuLC7dZMSAijD1yoQ+zMEZlAeXQdxTe89ICykMsBj3AIszJIItpx8ELQ1dCH4qLMLvFknO4QwhQGCT0wcwj9Q3tC9NeXYTZGYWNt4hpJmGgwtrMpfcxLla9+QLKLghi2b55kMMojJEWRXbtafHjvTSGWeR4LyyUnNpToIRb0HgcTkRp8WPXHvsWTiKjmJjjVeqeVIe4WqZc95h5VGokKddDrj2OEYIcXO2UQiZi+IRpT7lYMQzCpvGEMdIJxwVY1NbGt9fAaoiTC4mINUiqWUqhFJv3Te0uvOrYMbYUPpHmUQqUGCIYaxf1TWHsatiVr1k2r28klQKzFiGwwdUn5drJ3GWIV4dKwbFXq2/i+2j/ojp24UQUahbXPtqr1T0uhGFcqJGkULOM/YklhVpNTP9UO6X6JtVDrg6p1EipvVLfJHfXPYcHj55+zFOPAAAAAAAAAOBUuEkHAAAAAAAA7Bg36QAAAAAAAIAd4yYdAAAAAAAAsGPcpAMAAAAAAAB27FKku7bWrkn6/0v6KElvIulVkn5Q0hf03n8rPGYi6XMkfYKkp0n6PUnPk/S3e++v2v6slx66flM3rodIpntIiWg25SzE0ri0MEmamSjCWXi+mGYW4n6OTDLr1EUfKqe+Thfr28wlvkrS0Swkx47X5z2bhSRYhegik1MUE19TvFAtCG9jMZQqJpFtnnIWE8pC+tnItMexw3tgU86q87Npa4UEVknN9M/prqHdpaIVU+kGCCOqqQYOu2SwQprZ8jnXX2RMIoupr25uIfUxJraaeYSrZxrbXW17eG9dIppUO1XEhOiwUw4QLO73kWLKmUtyHYUkt1FIHdszqWgpzcwlkUk+sXU/9E1j7JsTUUoic30laS+2mxS2mO4aXrtpz4loW7pQFVXqnmnsm2oWUxPEtNZQb4R2l+Sa+qa6x9VOs5Gf33QU5uES+cJ7vmh+jIU5E/V4AvbNqibYbzh2vm5vnlQf65hq6qurN2LfcNxVaqdY9xRqlpiCe/qaxaXdp9cSnZO6J9csrm+oWQrX4lRXVBJbUwJrSgu1NVxKr98L83P7QqqdKnVPqiXDNg2nsmF2p0qqfUpmNTXLKCWz7hXqnlRvpDFMLTMJY6S6Z2JORKl22gsngJTu6mqcVCOlRFmX5LqNuufR6wcnfuxrn//UI+xYa+1NJP13SX9T0lMlvVjSEyR9vKQXttbe0TzmuqQfkvQlkt5I0k9Lerykz1495qEzmDoAAAAAAAAg6YLfpGvLrwo8V9IztfwW3BN778/U8ibd87W8AffdrbW778//I0nvLelXJL1Z7/09tPw23S9JerqkrzqbVwAAAAAAAABc8Jt0Wn5b7t0k/bakj+29PypJvfebkj5a0m9p+euvH337Aa21N5P0qav//dje+ytWj3nFajxJ+ojW2hPO5BUAAAAAAADgyrvoN+k+cvXzG3vvh3f+Re/9QNK3r/73k+/4q4/R8jfif6r3/qK7HvMiST8v6Zqkj9vKjAEAAAAAAIC7XPSbdG+/+vkL4e//++rnW9zR9gGrny8Ij/mhu/oBAAAAAAAAW3XR012ftPqZIjR+Z/XzKa21/d77kaRnrNp+LjzmF1c/33CA+d3TH9p/je6/dnhsn5TMmsxN3M80RPK4RLTU/yjEGaZk1sOYfrY+jkstO+45D+fr7Ycz3zel2/jU19oh0U3qq0uwkqTuEj0ln1JUDD6zwTQxvSuknJn+MaWrkOIqSaOp6xtSmApjpzEqKWeVvvE547ZOSbUm5SzFTKXdxu1n20w+SylnMUXU9I0JZZsnl8W+KRmskHKWElvdGCmVrockUncsLSYhZdYP7QPN0v6bwhaL+5kfPJzjXHs6l7mUSfmkuUqKqyRNTHJZSjPb3/MnHJdGdm28eV9Juj5eP/GlFNdr4eSZkstc/0nYKVNymeuf+o5TcuyWTjo5kb6S7nr6uucwnBRSumvq75JcD+YT3zfUQ67uSX3H4dp4lGIOjXApVjfbqVTfyJ8rcuBruvhs2KbjUl8Lfct1jzmXpXTXSt1TSIKVUt1Tq79cLZMTYjeve0r1zTH9B+FqlkJ9s+xv2op1j/v6TExPjUn16xuqmmrvnjPVSOk9T/O2fUO7fQuKn2/S+cnVSXEXSy/FprtuXt+k9kqKq+RrmVj3hJrlmhkj1j3hhOPGTn2vhfrGJcSmcYaoe1J9k2xS99zcP/7ezmbPc7H98urnW4e/v73VR5Ieaq09oGWKqyS9PDzmD1Y/z+QmHQAAAAAAAHDRv0n33yW9s6S/0Fr7stU6dJKkVaLrX7+j7x/odTfoJOkVYczbN+neoLU26r1vfHu1tfbwJt02HQ8AAAAAAAAX0ktai78zJUnqvT945/9f9G/S/T1Jj0h6qqT/2Fp7l9bafa21PyLp30n60FW/h1c38O783ZP0nf/b7SNt94vVAAAAAAAAgKQL/k263vvvtNY+UtK3SnpvST9xx1//mqS/K+kLJb141fYKSUeS9iU9FIa9/W27R3qPq/ek+Tx4rz6ttSdJ+r3KuAAAAAAAALhQntZ7//3KAy70TTpJ6r1/f2vtHSV9pKRnSXqNpBdK+i5JX7Hq9oOrvr219jJJbyrpiWHI2zfpfnNLUwYAAAAAAABez4W/SSdJvfffkPTld7a11u6T9EGr/33eHX91+ybdO0j6NjPc269+/sKws/TeYP/VenB//0SPTQllc/NbzCnlLLabCJ/DEJ+Y0119/0cX66/XpZZJ0qMh/exWWx8jps81P8Y9fjX8Ln5+C5MklJMqY4yVaTv90oXp5cV2E5CTkk9ju0lxXbZvnnIWxzBpdWMzbuor+eSy3DfsTy6ZNSVHhbHlktLCF3dTyqxV+/JvTUzpC1MZm9UUUt+Y2Lo+RkwoCyliI9M/paqOZuGcak4hKQzKnRMkhcUT0vu1eaJcerqcvLf5M5Z3J7eARko5C+0j0z4OiWh7MbHVJLOGFNcbE3/CcYlm9+0d2b4uxVXyyWU3xn6MlO6a2q+bE2VKOau0j8IOMg5xeilB7bQq9c2y/wB1j2k/WPj6IdU3se6Zr9csab9Jqa+3Zutj7I3881UqiBj4HN4Dl5SY6ptY96RkSzuIb3YjlM97pj3WJoUU12X7Zm3LsTdPbE11T0ykN/1L6fWhf6qRUj3k3sc2D2PknTL8xQAKdXdMbHVjxPpm89RXlzC/HCOcD83YKXk+1U7urUnpqekt9+9XLTHXpt2mIcJCVml/Ku1NMd3VjFKue9Y3diXFVZKuTVwya0ie3/MnouuFuudGqnvMGDdGtbonpbteMx/iJmmMQt2zjfrm4X3/misu+pp0x/lKSU+S9P2995+8o/1HVz8/aO0RS89e/fzhLc0LAAAAAAAAeD2X8iZda+1TJf0FSTclfc5df/3c1c93aK29612Pe09Jb6HlunXfte15AgAAAAAAANIluknXlp7ZWvsOSf9M0qOSPqH3/uI7+/XeXyTpe1b/+62ttYdWj39I0jet2v9p7/1lZzR1AAAAAAAAXHEXfk261tpflPTJkt5Ey19vlaSfkfSxd9+gu8OnSPpxLb8196uttRdLejtJj9Hy12E/f6uTBgAAAAAAAO5w4W/SSXpzSc+Q9Fta/orqC7Rchy6u2Nd7/93W2rtI+juSPkzSO0n6dUlfKunLjnvs0N5w71V67OT4t2FuFkSWpEX4IqTrnxZKPooLKK/P6aCnBZTDIscmIEKSbizWF5t8dOTH2BtdC+3rCz/umUWVpeMWflzvnxZKni/8tt5zi7OGvnlh5TC9Cjd0etlhAWW/qK/v6gIYpNrix2kB5fHR5gsap74psGF8tP7i82LLYUFjs5puXEA5tGux+TziQsmuPa7eW1RZ1HsU9ne3gHJx8WM39iL0HYV2tzjzYhbOqSFQwi2WPA8rF+d1Zk3/lDGRFkV2WRxhoeR8nKcn3dKC3OGtbSEMYjRen/jYtEnSJCyK7EIiUkDEjbCAslss+f4Q+pAWUL5/73B93LCAsguCkKRrof16W29PY+QFlNe3037oOwo7VApsOq1q3eNqGVfHLNv9QeNCIqp1TwqauGX2kZszX99MYgiWCWAaILgjXTYWoZZx4TgpMKeHGmIQtu4J57dCcES97gnt5nCMYVeh3fXPQRWb1zKpRoqBWa7d1DHHjSEXPpGuO2Fs33eLdU8Kk0hBWq7GSTVSrFnW20cuiEvSYpKuo6Zm2fdjpP3dfiaIb0uhZgznrBQg5uoeleueYnuFq8sKwViSD4kYmzAJyQdjST4k4r6JrzdSGIRrT/XNA+P1+ib1T7VJqocqdY+rY6Rcy9jgiHTNPUXd8+oQ8FFx4W/S9d4/T9LnneBxr5D0Gas/AAAAAAAAwM5cmjXpAAAAAAAAgIuKm3QAAAAAAADAjnGTDgAAAAAAANgxbtIBAAAAAAAAO3bhgyMuuiePH9bjx6+LqZkXEnIWIf3syMTepDSzlH7mktJSytmthU8ouxbaXerrpPm+KblsXEg0W4S0L9ee+s5DCpNLfW0hkcdHFMkmHfXw+lolQSlsohRKU0o5K7a79LOUchbTz0ySa0pgdSmuktRMf5taJqlN/Ytx7THNbBY2iE13DWOk9LPUvzJG4pLLUuJrSDlrJtGsh4SytufPT3283j5OiWgTP8ZiYo7RGGcYjtEBIpjdqSWlmaUEsD43/QvHs3RM+mxlFymchtITVtLPyilnLt21kOIqSQ+aZFaX1irllDOXXHZfoa8k3T/y/X26qx9jIr+dXPpZSkobhx1kZHa01DdxdU+qb1KN5GqZmF4fogEPTG2S6p6bhfpGkiZzt639+5JS7YdJcjWpj6Humc39dpqbc3s6no+JuzZ9N+8qyadPFused/6MdUxIYK+k3VdSXKVQ94T6ZpwS6V3dE/uGF29qnGb26WXfUDu5miWluMbUV9M+RH2ThJoljuH6h3TXWPeYGifVNy0cowszRqp7Uuqrq3vi4Zw2qdlOqZyqfK6ofI5ZPmloL0ify+xGSeVyqGVc3ZPqm/2Qan/d1DjXQzJrqnsesKn2xbrHtFfrm/sKdU9OcfUnZnfd3Q81kqtvkrvrHpe2W8U36QAAAAAAAIAd4yYdAAAAAAAAsGPcpAMAAAAAAAB2jJt0AAAAAAAAwI5xkw4AAAAAAADYMdJdd+xxo6meMDo+ASRli6SUM5fUNQ33Yw9C+plLNDtY+JSznKyyeXs1Dc6JCWUL/xpn5rVPF37ORyFByaUOjlKSU0psdNPeYspZJfU1BdvElLOY7moSykyblNPPXBpZSnFN6WejI5PMWkhxlSRN1xODmmmTlJPLXPpZSErrMfXVtFdTzhKXXJbSiUNCqUt3bSatVZIUUs5s+yRctsK2Hpnjv4dzQjzZ+pFtaw+pby7RLJx+1VPKWeEYHSTFtco8aQrBS6e4kTmn7oWUrElISpuYa+q+i5iWdCOkn90Yr6eRpTSzB8cHtv0B015NM0v972/r83PJZ1K+Fl837aOwg+yHHW1b/9KbDsWjEA24MHtUqm9S2v2B2aY3e0hrDclx+4VEuXEloVs+mXUWtscsJFUfjdZf+8S0Sfm4m83NuT0k5qbj3yclVgsfN7BvrqS+xnNqTH0N7WZXCKchm+K6bDd1z2E4FiuJ9EPUPSm9PrWbGqenGiklsJ+TuqelpHpX46T6ZhrOnqbGaaHm7qkWN2mwaZOWhO2REltd+2IcUu3Ta1ys96+muG411d7uNpun10v+s2RKtZ+E87JLfb0+rtU9Lsl1m3VPrIdaSHcdmQTbUPe4+kaSJqa6GIcdJKba29bX18P7V8E36QAAAAAAAIAd4yYdAAAAAAAAsGPcpAMAAAAAAAB2jJt0AAAAAAAAwI4RHLFjjx01PS4EDdw2DwsX5oWV1/9m6hZbVV5Y8cCsXH5LYbHKsMr5KMxwXFidPYVjuEWi02LQKTjiyLQfNn9I7IUFIEeVhdIHWBO5JARpVBZQjQsox0CJ0N+0p0WY0wK5LmiimTCJZXthAeWjsJJzCIOwIRGhb5+FsV17NTjCLbi8zQWUUyBKCkqwCyj746vNwtgTE1aTVj+O7b7ZPl94LSOz0HEfb76fLvubqcWAiHDsmteY+sYFkasLK1fYEJywsHU6p5qFlceFxZYlv4Dyfjg55eCI9fb7RuvBAse1u0WRHxz5xZbLgRI2OGLzgAhJmpj3y8dDSZNwbIyHWOzfSHVPqmXcu3g91CYpUGJi+o/C86UFpSshWIuw2vo8tE/NOdUFYEnS4dyfa92xcRiuMem4c8doOp5j4NAAu409tQwSmBXOWTFQImwn056uD6W6Z1asew5dcEShvpF8jTP1585Ys1TqnlTLnJO6RyE4QpW6J7Wb15O2h+ubVD+bjEwoWKp72jzULObtTZ8TwqnWH3flgIh0wtlOklZ5W5uJp8+de+FE5MKxKvWNJN1nArPuC8ERLiBCkh4zenR9jFDHPCbUQykM4j4bHJECIvx76+qe/VT7hwvVJnXPLFw/K/gmHQAAAAAAALBj3KQDAAAAAAAAdoybdAAAAAAAAMCOcZMOAAAAAAAA2DFu0gEAAAAAAAA7Rrrrjt0YTXTfyKdz3TYP6T2LEFt4zaTEHIbYnKMw9lgmhSnd0g2JPHOTDLTsvj7QUUooaz5rbtrW5zcJCS97IUrIJeektECXvCNJI9O/pXih2O6bt6YyvXJSWkpFc6mUtZQz1x77VlJfB0g560c+icimmSmkvsZ01xCF5bZfiscaQkrpK6SctTC/3kPK2UYTu0dfdx6ahSSncPx3k+7qEl+XfdMxsN6/mp48RMrZICmuFcWUM3f+TOflnH62fsxcG/tj0fWVpOsuRcy0SdK10O4Sza7HhNjNU1wl6T5zDbwedpzrYVtfN8f0KLwxE/lr9HhLseWp7pmG92th+h+E4mQUkupHlaS/cNpLifSu/Sic91JS/aHp/2iokdL+fmBSX9NxlI47d4zG3eCM65vyeW+AVPtKeyXFVfJJrqOU7praTQ1RSnGVbJJrKb1eUndjV1Ptz7rGCZ9jbIqrJI3NZ4JQ6/ZUA2/YJuUS3SXVp/T6Fmq4Us1dSEROx0tMYHXHaPGzyVZVUu0LnyXT585KPZQ+/6bPy9dsXeGP5+uhNnF1Ukqvr6S4StL9Zi7Xwm6TEluvmeM8pbWOwoV+k7rnRkqGLuCbdAAAAAAAAMCOcZMOAAAAAAAA2DFu0gEAAAAAAAA7xk06AAAAAAAAYMe4SQcAAAAAAADsGOmuO3ZNE11r90h3DQllixBjM+0mtSUGkWw+9iSkKqWUmElIUJuYdJb9NEZIlRmbeKA0j3FK0zGvMaa4xvb1tpjuet4NkO5aSZoMb1dOZ7NJUyHNLCQDamGSt1KKWGw3E1+EF5OSWU17THGtpL6m113lkouK+7U75cQkstDeTQpTM+lpkux7K8m/jyGhLI3RFubfswY4NsoJZZUxLqhKGGRK43Tn69Q3XTf89cHvH/vpOuWuo+G6uB+vl+k519tTiuskJbaa42sS6pGUcjYy79g4JUEHc1NbxGS7sIdMzfabFw/ShamHpqFG2nd1lvL76/aFtN+k/cztk5X9dzn25sdGctaB9Ns0SKp9pX91bJeaGa91YRB3DUxjhHabtlpJcZVC3ZOeL9RD7nPIFuueFlJcS7VMSoJMybFzlzIbElhTIr2rdcN2inW0mUeuz0P7tj5XXNC6J6a7ml2h/HnUbJT0+Td/Xnafrf3xXPncHu8TxPsH6X6Dm8fmKa7LsdeP6Vz3+LE3qXGuDfA9OL5JBwAAAAAAAOwYN+kAAAAAAACAHeMmHQAAAAAAALBj3KQDAAAAAAAAdozgiAsgLVC4CAsXj80iitOwyGblLm1agHI81KKtRl44Myx66/pWFlA+76uwX6YVmxO3UHLVEPtkGKO79jBn21fyizOnMYbYHlV23uGY636n7Ka9vPu6RaLTgtdxjO1svxhOcpmE9/as0yrOy3nZLf5/bH9znXJtUr6mpecsXbsHuHCcZgHle7FjhKCqytGfQyY236Zx+6f3q9hu+56TVdHP/LjLCUJnOo1d2Or1pDJ2ur6G47E2DTOPNG5sP9t9MtVfMSDOHOepDrShIKH/IPvHEGPsoh49a+f8dFMJlKiGArnr2javA+lal+43ONUKxN0jyX13+102vkkHAAAAAAAA7Bg36QAAAAAAAIAd4yYdAAAAAAAAsGPcpAMAAAAAAAB2jJt0AAAAAAAAwI6R7rpjcy00v0dy1SKknyxCztnUpL4uQqrPNLTPTXOa5zxE4Sy6vwc8N+3zcL942scbj53GWMT5rbe7tuOcedDRWT9fNeUo9LebNY09Cn/h3t7UN6X3mPYexmhjvz+1mRkj9Z2H5NOx2a9TApj8MeDDnU+fwLZ80vXX0+K2Dq/dvQej8O9CKUHJbac4RmFfSP88VdlvYt8w9hDOOnUsJWyd8XkonZfTub3SN103fN8wRuVaF/qma928+RT3qZnLJLwx85hmatLgUgJgmIerLVISbOJqnFTfzFPNYsZIr9u/Er9N03se368B9oX0nL5vre6p9K3WQ6eWzisx6W9L89vBOTxfT0x76lu5rsUxCtfieag3whiuJuipfvCHl09bHSB5dvmkpu5xr1vKtafpb2shSUp1o+sf+qb6tfaeD1Fz++ZS7Z+cdd2zgwDbymfJWA+5z7TFusd9to7XujCGu67F+wFhn8z3G9Y3VEprT8mxU1MBTMI+Ng3vyyY1znyAz2N8kw4AAAAAAADYMW7SAQAAAAAAADvGTToAAAAAAABgx7hJBwAAAAAAAOwYN+kAAAAAAACAHSPddccO+0yHd6SVpUQzZxoSjWz6SUhEOwjJJQcmieVW97vLQWyf2Pap6X+wSH19qozrP1uE1LyFH2NmE2h8Yksau5v+rm35F2cdURSkaVRSmMoJZabv2O98Ka3Ktce+e+HfH8z72OYhvSsll+2t779pM8VMOpMQG1PV5j6L0KXP9pTMmKKjUlKXe76YUOa3n01FS0lpZptKUnPteyHtNrTb/mFbp/3G73v+6eKxYfqnMWK7O47O+ekmSbuqa07n1Hy+Xn/PU1qY6ytJh4v1fe8wXKeut6ltd9fASZ/ZvuNwPZ/4GGc/hvzYKSH+uqk3Ut+UZjbeUvReSmaNSfWm/0FMtffPeWDqDde2bPf7QqU91TdpP3P7ZNp/K/t7Oo7SceePUdv13KicJ2Pf6vnatZfHNteecJ1qszCIuQamWqEt/HnITdsmrYa+y78wfxPqm9he2dGKdU8pkT7WPaZ/pb5J/YvzcPtIpb5J7bFvKL/ce175nBDbi59Nzlz1nGra0+fOWPe4ZNb0+TeM7T5bXxuF+qZQD6XrYqqHMtc/3Tep3CMJ9VdIn3ZGd33v7bBQu+UxAQAAAAAAAOwUN+kAAAAAAACAHeMmHQAAAAAAALBj3KQDAAAAAAAAdozgiB17eDHV/uL4xQXT36alEt2iyEdhVdppuE/rF1D2Cz/eWlyz7TdD+8OL6xuPcWvu2yuLeh8WFs5Mi2ymRZjdYp+LReprm3O6wGm1GFvgpzHEAsph8Vi3WdvYDz4KgRKLiQl9CAsDx/WC3ZuweVbL8jndordusWCF8ARJ3bS3uIBymKBZ6DRtjzK3OHBaQDWFT4zMa0zhDilQYrJ+nMeACNNXkvq+GWOS3pewUK/Z91zbsj0trLx5+ERlEfH4T227CJSwK8qnviEMwpyX52GR49R+ZE44R+aaIUmH4frrrjFpoeRJ89epkTkHj6snnGBqFmeehmP0evOv8aCtz2US5hd2ayutIZ5UlldOoQ9Ts8PH4KnYbt7zUPek+qZSD6W+aT/zdY/fr9P+7o6N6vHljtFcLPjmbdU98fyW2m2gT1pIP4VdhbmYc366PrhQK0kyh2i8zqfrmh03tMeq0QQXxDCpFO40Wz/SW/r8k0Ip3ALvu6h7Qq3gapkcELF5PdRTfRPbTUBMqntiLWPCJ2KNH67npn8aY4jArCh+HhqAGTpmuMXPh+vt6XNnJQwiff7Nn5fXgxnS5/BxOFuM3UmraD469O3uOh8qiIMQBrHv6p6wf4zCGJucaR8OITwVfJMOAAAAAAAA2DFu0gEAAAAAAAA7xk06AAAAAAAAYMe4SQcAAAAAAADsGDfpAAAAAAAAgB0j3XXHXr4YaxHSV+5lEZJfjsy915RyNjVpZpJPNEspZymJrJTuGtJjHim0PzoP8yu0T+ch3TUkb7k0HZfSs/wL31xKRCyoJrP6lLMwRiHFVZLa3vpkUgDQfD89qUv18n3TEbVo6/v7KCWUhbQqm8hlUsskSSGxtc3McRdTXMPOYJNqh0mOtEy62/HtLgW3lqrqtnVOZt08uWyx7/vmdnOcm31akhaxffO+6fgqHaPFlMNtpb72kGaWrl8u/Ww2D+nk4fh36WcHM38d2AsnopFpd2mtx3FJZPG6HRNH/bzvN+ln103iqyRNQrrrvklFS31Tklt1m2wqbSe3TSVf4xyFK0FOfV3f1qmOSXVPSmw987on7O/u2EjHUTru3DGa3q90/A/BPmXx/GbPn9W6J5zHXQprC9sjHHZaxIJtXXqNNrQ01D1tLzzf1EwwpJO2UMvYBPuUzJpqmbOucVJ9GOtGVysUU+1d3RPel1QPuVompbjOU7qr2a9TfZ7rHnOuSOmuldTXdIxW66GKAT7Dpc+H7rNk+tyZPqcejNbP+XujlNa+ed1T5a7R8/DGVK7Fkq9xct2znlS7bF8/D+2HVPvT1DevTMnVBXyTDgAAAAAAANgxbtIBAAAAAAAAO8ZNOgAAAAAAAGDHuEkHAAAAAAAA7Bg36QAAAAAAAIAdI911x35v9oCOXNLjHRbhXmpKOXMpKinFNSWruKS5g75v+x7GlDPf37Xfmvu+j4b2m6b95iw8X2h/1KSfHYbUnOkspL6a9LNFSESLKWdDxA5VUs4KqUgxbSkkn/Y9n4TjArlSiFCL28Mlb/meaX4jl7wVkkVbSM0azczEXZtCmpkkzQvJrCn9zKW7piTYIRTSzFL/npJgU3KZeR9zumtIHXPprimJLKWc7ZuEskk4/4Z2m3IWTv2V9OQY/ldNd66chmLKmRmkmProzqmz8N6ma+fYJHKNBojMTkmL6Tp6aN7gw7G/Xt438ome10Y+uexhl3IW+qaUs32Tcpb6xnTXkIp2WtW6x9U4KTE3ptqbWialzJXrHlOzpL43Z35fcEmuj0xD35DuemiOmXQcucRByR+j5frGnit818gNMcD5MKbUx3TX0N9tk+JrtHVZC9fLVNuZ6+tommqWkNhq0l2bq2Mk9VQPuRon1SzVtPshuBon1T0uMleyCbGxXo7prq5mCe95aLd1T6hNUlK9S3JNY6RjwLWH0+8x28m0FdNdY30zSOqrqVNTqn34fOjOqelz52FIbG1t/ZxfTSdd2ET6Yt1j3uBDkzwrSbeKdc/1drTW5uoYKSfVuxpnHD6QnqbuedXM11IVfJMOAAAAAAAA2DFu0gEAAAAAAAA7xk06AAAAAAAAYMe4SQcAAAAAAADsGMERO/ay2RP08NQvqHjbvHgv1S3mWF340bW7RZUlvzj28e3r47gFkY9rP5ivj32QxggLKB+YxZIPpyFgIyyma4MjwoKhcWFlF6pQDZMoBEdUFlBOC7ymjIO8SLRbfDNsp7C7u7XjU1jAaBYWNDbBFguzOK6UF1ZemCCHFp6vpdAHN0ZYKDmu+3pegiNCuw19CIstp3aZMRbhPU+LHy/cPGJwxOlDH2JwhDkNhVPqMQs8m7YY7uLb4w41xALKTlr/u7Cw8iwsoHwUXksrLJbsFkpO7UdhRfnD0P6oCYl4ZO4XSr4xDgslxzCI9UWR02LLeQHl9fZxeMPSAtSp/2mlumcRrjGuf6W+kXxtkvqmeijVLK4eyvXN5u2VgAjJ1z1H4fhKx507Riv1zWBs3ROu/WlRepdlkOqeWJdtfr7p8Toa+ptrYx+n+iacy0x9Mgp1T6plRq4GDsERqe5x7Sl8ItYy5yQ4Ipa6rj4p1j3ufUx9U9iVO21V6pvUnoOxbLOve+Lz+TFs3TNUQESh7kmfy7rZr2N9UwjMmo5SiMsQIVh+Hi4s6HDs35gbIQTr0fF6INIjoTa5NvLhCqnd1UOV+kaSRiYkIvVNNql7Hp7511zBN+kAAAAAAACAHeMmHQAAAAAAALBj3KQDAAAAAAAAdoybdAAAAAAAAMCOcZMOAAAAAAAA2DHSXXfsN48er/uPXpf2No+RNesqKWepb0oum5m0ummYm+sr5XRXl5B3ZNJaU19JOjT9U5rZUUhmnZrkspTiOp1unu7a5yEuKLW792aAAKucPhWSwUzqUAon1fGBxOtMQlYfpRQ2P4RLZwsBQGphW4/WQ4diutgojNHMc+Y0Mz8/l1CW5pHTXUO7HaO2Q9kEusLuK8kmmsV9MqW7mlNOSgaLSWku5Swk/cV0MZcym9LMYrtLmQ19U+qr6R9OkTH9LLa7xmriq9vfU9p1OL4WZt+bt5BO3vyLd69lHuYxNWlmkr8WpJSzWybNTJKum8TW/ZFPEdsL7ddC+8S07xVTzvZG629YNd3VGccTn7etuscl1UnH1D2mfVpM9E31kKtlUoprqkNs3RP6psRWN3ZKcZ3PQsKuOXZj3TNE6mu6bri24nnPvl3l+mvzBNBc94SaxdUbIYhwFJI3R6a2SGOkWmbk9oXwHsZ6yDVX+kr2vRmkvkmKdY99H1NCbNxXC7VCSi12NUsxCb6Sap+TWQtjpLrHzC9/TgjH1xDp9XGfNIOH82GqU+ej039nqpt5zAv1jSQd7q2/OQdj/0HrlvtAJWnf9E91T0pxTfXQxNQWuW9Kd11/I7dR99w8Otz4sQnfpAMAAAAAAAB2jJt0AAAAAAAAwI5xkw4AAAAAAADYMW7SAQAAAAAAADvGTToAAAAAAABgx0h33bGXPvp4Xb9149g+i0Ly2bK/SdMJEUUp/cyN4ZLPjhsjpcq4lLOYQBPaZy6hzCStHjc/l2g2T4mDIQlnMXMpZylGrBJR5rsmLlwoBCLmRKlCylkMZouJV6ZrSqVKyU82oSwksKb0M9PewvsS01bNi08BQKV015SUVtkXBkgFXj5p4SlTX5dsV013tel4vmtMOTP9S8dAaM+pamGMSjJrTHIrPF9KOUuXk0r6WXgj3b7a0zGQ0s/MiSsczjbNTJIW5piej2spZ0fj9Wc9COmuY5OSKkmT0bX1NjOulBNRXQKrJO0VUs5SQplrH4WTSB6jluS6qVT3pHRXV+PEvqHdJbPOYqp9qHtCf7efpfomJ/Ktt6f91yXPp/Z5SHd1Ka5SqHGK6a7NRp8WIxgL14d0PnTPWa5vUji5PV/XksV9Mmut7nE1Tjps3fPF/rHuKSS2pjEuaN1j979qumtpv9587Gq6q2uP+2+hlolpsqV5hL6pvils66iQOJxT7cN1zXYOY4R5uBonJpy7k5Okw9n6G7YXapZJqE1cPZTSXVPtlOue9XFSbZLGcDVOJcV12f/edc/Bo4+WxrTPc+oRAAAAAAAAAJwKN+kAAAAAAACAHeMmHQAAAAAAALBj3KQDAAAAAAAAdozgiB37rZuP0/7kvhM9Ni2a7VQXUHZLKKbFjNMY87Bw5sKMM48LgIfndAuDhwU5e2EeaaHkvAioaa8uoOzaqwvhmkUve1isVmnRVvOkLYxRDaVoLjgiBkSEMQqLH+cFlF3b5gtKxzEqC8oqh1WUxhhqseRNp1Fc07sSHBEXbK4soBwXBjaDxPCJMI8hwie2NUYKiCi8FkmlRbPTvmdDItK4s3C+Nqsip3N4X/gTgLs+zEYh+CAsLnw0Wj9BjcdhIeI0ttlQaaHkcRijhY3txkmLH6e3wAZHpDHO+oQTxLAQFwBQqG9S/1T3pHmkuseNE+uyNIapcVLfVDu5Y8nWMTomBKtU9/hm117exdxTFs+pdm8IdU/YpLkecq9xEsYINYt9f2Pog2/32zoFRGxe9+SF9DcPFqrWuuem7qmUcMVrbiU4Io5hwrjOTc0yQHBXte7J72NhhyrU8/H5ZmFos1EW4RhNn1MXZpvMQ511FOoNV+PE+ibVMq6uKNY9qQ5xdU+lvjmu3TlN3XN0cz04rIpv0gEAAAAAAAA7xk06AAAAAAAAYMe4SQcAAAAAAADsGDfpAAAAAAAAgB3jJh0AAAAAAACwY6S77tgrHrlfe+P7T/TYEPwS+ob0rkJkY3q+OHZK5DP9Y0hUSsE0D8gJgIUIpZhEVkhsrSZ62tSsYpRmJQ2qNEZIAAsJSjFttZBQlsawIWypb+CTWTffx+I8qumuA6ScVZ5vEAOkuyaVBLWtJsQOMY9C+mxOqg1vZGWM8vxcQnQYI7HHaErYTDvr5ueKPvcnIrv9wjZdpJOZ6d9iAqsfwiWDxbSwmKpaGXvzvkk1zSw952lV6ptl/0otk2qT0/Vd/kWh7imObWuZWGf5MWz/Sm0ilWqnWPdU6td4LnOD1BJ9XTJrC8mR1dfo0+TDRNLYhVohnloqNeYANUvpFFKtWc5TjbPhGNWaxfYfYIzqZwLbP6bJhkFszRKSWQu1U6xjqvOrvOdp8IW5FqfzUOG83OepDvTtLox7HqKnUy0zc+2prkj1kGtL9U0x1b5Sb+QxKvXQ5s93t9kjpz+h8E06AAAAAAAAYMe4SQcAAAAAAADsGDfpAAAAAAAAgB3jJh0AAAAAAACwY9ykAwAAAAAAAHZs43TX1tqo917MUnztYx/svT98ksdedrdec01jXd/+E1Vj+raazrR5ImrpOYd4jTGtqpCEVUz6dPMuBuzZl54SbPqo8FqKSUnxHbCxeWGM6vYr9PUJZbX9ZqvJZdsaYxfOOiktjrH5BiylsFWS2apjhCGGmUc6ONKTusHDEC6hrJI4KPl/OgwRWz2+FndC9F1Lh3N121W26RDvyyDH3AU94ZRPDG6MAfpW5jFArVC2rbqnkgSrgRLOC4n0MabPlSbp6wvxEE1xvIUxKvvCIKmqA9Q91DdnNsYQdU81ZbbS96LWPZXtmo4vm7odCp+Ydu0jUf3zVWqF4ra2Y5f3PddWPNC3WfecUY0zf83s1GNUvkn3wtbaW1efoLX2pyS9uPo4AAAAAAAA4Kqo3KR7pqSfbK19bmvpn6Rep7X2xNbat0n6LklPOekEAQAAAAAAgMuucpOuS7ou6R9I+m+ttbdIHVtrHy/pFyR9lJZfRDzRr8kCAAAAAAAAV0HlJt27S/opLW+6vbukn2mtfdadHVprT22tfa+kb5L00KrvT0l6j2GmCwAAAAAAAFw+GwdH9N5/orX2LEmfKenvSHpQ0le01v6MpE+R9CckfYmk+7W8OfewpC+U9DUnDZy4CkYPTzTqk11P43K7qAvTnlJc/DhskPOymc5+3fLz8spxGQyxtvuF5dYcLq0GLWm++dO1IVbvvsrvFy6vK3BZs+fanF5VaD171D24yK5E3VMKbNneBqHuuRj6I6e/t1P5Jp1674ve+1dLeitJ/0bLt/l9JP2SpK+W9MCq7bmS3qb3/lXcoAMAAAAAAACOV7pJd1vv/Xd67x8v6cNXTW3157ckvW/v/WN67y8baI4AAAAAAADApXaim3Rt6dMlff3tptXPp0j6pNbaHxpicgAAAAAAAMBVUL5J11p7Z0k/IekfS3qClivI/ANJP6zlzbqPk/TLrbXPb61tvOYdAAAAAAAAcFVtfJOutXajtfYPJf24pHfS8obcz0l69977X++9v5+WoRI3tVyb7u9J+unW2nsPP20AAAAAAADg8qh80+3Fkv6wljfnZlomuf7d3vv0dofe+9e11r5X0jdI+uOS3kbSj7TWvrn3/pzhpn157D3ctLfYfszKIMk7Q03TjFOeX6V/6DtIIE/ltaT4rm1tj/S6K/Mb6n2xr7GYLlZKchtAZX7n6fg6z4YIlKuOMcTJr/KcMQGsMI/0fK499q29bru7V+ZRHmPz+cVDsdBe3gsqryUoneLK+3Wxv3H2iZJnayd1T+Uw3+J1Pj7npuMe016a9xB1T6mu2GRSd45ROFlU241Y9wxQ05YMUX+Vn3OAMc67s65xthmrOsg10B2khedLinVFqYYY4jo/xPwGqAmKYden76uBXktygeqedvP0x2bl113fRMv3+2clvWvv/YvuvEF3W+/913vvz9brvlXXJH3CqWcKAAAAAAAAXFKVm3RTSX9L0rv03n/6Xp17718n6R0k/chJJgYAAAAAAABcFZVfd33X3vvPVAbvvf+6pPddJcECAAAAAAAAMDb+Jl31Bt1dj/0nJ30sAAAAAAAAcNlVvkn3Wq21J0v6VEnPlPS0VfMH9d5fNtTEAAAAAAAAgKuidJOutTaR9BWSPlnS/u1mLfM2rpn+XyHpgyX9ud77C0831ctp8nDT3vyECSBbTFuyATTVFLFK4lX6Tmcl1as4Rh+tR7zEOYex3XO2uE39X/gxfPxM5TW61ycpf3+2MkY1wdY95wBjt+IYdrvGvoV5pK6lhNha5FCc3xnrAySR9cLJIj5dHHvzQXpK2jZjx9cdxrD94/OFsd1rXPiuLbRXXnsMIqykqqbXmOZdeL/Sa7TtKVWtsJ3i6x7gtSSlbb3NhNjzkgS7xRTMSt0zRDJrud6ozC9e503dU6hvJH/tya8l1T2bX0fT2La9WBPYGqew7Y4b286l0lfhOh/mUamHKvVNGCKKdU/hBHUh6xupdGJIdY99ykp9EwaJdValHhpijHjN3XzsWFeEDdJjrbB5+myuFUztVKkJwnPG5yuMnbdTGMPWTmEe1e3k+g6RYDtEWnDh+TbVbp1tuqskfaekT9fyhlyT9PA9+j9Z0ptL+jP1qQEAAAAAAABXw8Y36VprHyXpQ1f/+wJJz+i9P+4eD/seLW/mvfuJZgcAAAAAAABcAZVfd33O6uePSvozvffZBo/55dXPtyzNCgAAAAAAALhCKr/u+gwtfzv3Gza8QSdJj6x+Pq4yKQAAAAAAAOAqqdyke2j182cKj3mjQl8AAAAAAADgSqr8uuvvSXpjSU+V9FMbPuaPrH7+bmVSV8n+w117sxPGh1TSUwspqcv+Lg4q9C2nqp7tGLndDJ76jtMYJikt9K3cEo9jpP6l1LcBEljHaYyQDGb6pySy1D4qpLuOYvt67FBKEXN9Jb9Zx6lviB0KgXelMfy424tgXBTjDF2S2CIlgBVSX9M8UrsdIyReVdr7wh/Q8bWYXaTPwxgpkWtu2orpk7VUtNp7bscups828xrTnF3fNHY17Xbkxq6+FpvoO8AYkp93DFVMB17h+ZJtnXK2msyaDo7KGKE9XeddzTLAGItUK8SxzblsgLpHqa9vVjMT7KmuSNzOGuu9Qs2Stn+YX6xlxqbeiIm5aez1MXJ9s3l7rk0K9Veac6HuqdQ3aR5DqdQ4lUT6at0zN7VFDBaPdcjmfVPd4+qQWGcVxujzlPie6h4zRtr88Toa5rf5R92oUm9ss+5xr708jwFqp22NkfpX6ptl/zB2YYxNjB49/bmq8k26H1/9/IBNOrfWHpT0OVq+xB8rzgsAAAAAAAC4Mio36f6VljeY/2Jr7UOO69hae0jSv5b0lFXTN5xsegAAAAAAAMDlt/Gvu/bev7+19t2SPkzSd7fWvlHSt9zR5e1aa8/UMmDiMyQ9Qctv0f2L3vsPDzVhAAAAAAAA4LKprEknSR8r6XmS/oSk56z+3P6l2+ff0e/2L1L/tJa/8goAAAAAAAAgKN2k670fSPqg1a+7/n1Jbxu6PiLpiyR9de89LXEISZNHuiZ3BkdUV6w0hlhA2S6yW1zkuNKeFy5O83N90xibt8dFmNMilm5+aa3IvbSKpRkiropaWAmzsPi0JL8ocloo2SyIvGz3/cd2AeXN+6b2vRDYkEIfxuY5J2N/ihoXFlBOwRFpkeOxWRW1uiBypX/qWw2DqIzh2tN+vQgnl5lZ6NgtqnzcPOamfTb3B/o8LNTr5pEWYZ6HMIi5Wfx4EY6XtLByH5kFpWfFQI/QbhcpLp6GKkEEafHjkXntMSBiFtrdwsXx+Xy7O+envm0ezpOVBZSHWOC5EjKhsODyEIstb1ElBCsHRITjzo1RrG9ycEFh7FD3uPrEnBJi3+XYhXmk9bjN/BZxxymMnVMmNm4vBURItsZpoVar1j2uDhkXaqQ0Rqp7XH0jSXumxqnUN8uxC/VXeCNdSEQlqKJqiLqnWiO5/vNwgMWaxdUbhfpGkqamxkn1TaplbP0V65vNQ7Dm4aQV6x6XthLvLKQTTuheCcwK29rtZvG6Heo1GxxRCMaSfH0Sxyi057CL09c91XZbY8awkC0GZm0wbjs4/aDVb9It59H797TWvlfSO0l6q9WfN5D0Skk/I+k/9d5feerZAQAAAAAAAFfAiW7SSVLvfSHpf67+AAAAAAAAADihSrorAAAAAAAAgC3gJh0AAAAAAACwY6/9ddfW2mdJ+ssDjt20vAl4XdKDvff7BxwbAAAAAAAAuDTuXJPucZKeppyBIq1HnmyaxXROcsHOn/2bC01mKcZkKaaIFUKH+iil5oSUM/ecISktpojFlDOT+hjGiOmHLqEsjbEXXrvrH1Zp7KHdvnMpGCi9zWY7paSZyoGU3ts4uEkGS2lmo5B+NgpJqXt7LuXMj53SVvdM/9R3ElLHXP9JiGxMyWV7Jv4o99089XUU+sYUtsLeUE1KK6WfhR3epqqGk0Ul3TWNMQsnHJdyVukrSVMzj6NZmMcoJMeOzWuZ+fdlrrSd1veRHk+0xZSzAZLFS8lbIWnOpouFFNeYtmr6x3TXlBBrkstyuuvm7aNCItqxYy/M+Tqd8itJaYVEtCQ+X5pG5Xc6Kkn1lfpGKd21UD8cM/bCJKLm5PmULur6hg2Sahnz3pTqm6ClJF2zn0rHpPRWuDHSuIVU+1j3pGTWvULdk2qFUMvsm7Er9Y3k65BKX0naMye/arqrG+My1TfLsTdPZk11iK17wsnC9ZWkqamTKvVN6j8LKa5pbJv6mlJ33QlOxbonnLTSW94KhU8pPb2aZmprhdA31UOluidcYwoJsZV6yNUry/ba2C6NN6e7hvYwF/98aYx7P7QdFgsh485L80zSwQaP2dfyMnj7z2z1R5Im0ms/XXRJPyzppZIeOfVMAQAAAAAAgEvqtbeie+9/r/d+33F/JH2ipENJNyX9NUlv1Hvfv+Pv9yU9XdKXSppKeitJX9t7H/LXaAEAAAAAAIBLZeNfMmitvb2kb9byG3Lv3Xv/h733376zT1/6f3vvXyDpT0h6sqTvbq09YchJAwAAAAAAAJdJZSWQz9bym3Jf3nv/mXt17r3/iKSvlfTGkj79JJMDAAAAAAAAroLKTbo/ruW36P5T4THPX/38iMJjAAAAAAAAgCslZDpZT1n9rMRVHK1+Pq3wmCtl/Ohce/MUY7JUTsFy6akx3TWNYfqGMWJ6akpQMwE+KSG2FZJZF2lvDml1vZTw5NlAs5R4E7afSz8rpd1Jg6ScNZtydvoU12W7SSgrpJlJ0jXT/9rYxxzFsU380X6hryRNTKyPSy1LfVP/cdjLUnKZa09jDGEedpyUXObapymhLLS7hLKUiHY49yeA2Z5JZg19D0a+3c0jpdIdmZRkSZqGNNiKblJfUxpnTwlWpZNZjR3hvKScxb7h+DL9K0mwkp93pe9x/YdIULNJaSn4LI1dC1bcWK5NCrVM2KdryfMpYb6Wdu9Cn1PfNHYzp6dFSFpPG9CdLgrB88v+Lkg3bdN0nLuk2vB8kX3PN0+vX85j87qnkuIqhTT5MEaqQ67trZ90UiL99VAP7Zv2lOKaaiqbEBvmsRfiICemvVLfSGdf41TqmzRGSqSfxnRXk8wa+h6FDz5HJm31aBzqm5DMemjqocOQXt/C+zUNia1OuLyqu/ow1j3pc1YYvFL2FJI+K+n1cYxCfZPaW6G+iWMUE2IrdU9Md03ps64OGSLdtZL4eszYrzfk0fH3djZRuRXwqtXP/1/hMe+5+rm9syoAAAAAAABwwVVu0v2Ylvec/1Jr7Y3u1bm19pCWa9F1SS8+2fQAAAAAAACAy69yk+6fr36+oaQfaK29R+rYWntbSS+Q9Oarpn91sukBAAAAAAAAl9/Ga9L13r+/tfbtkj5a0ttI+m+ttR+T9KOSXqLl+nNvKuldJL2/lt+665L+fe/9GwaeNwAAAAAAAHBpVIIjJOkTJR1K+gQtb8K9x+rP3W6vmPgTkv78SSd3FYwP5hrP71ipsbB4d1xY2YyRFiJOi2a6BZRTQESbpXCH0G4WOk4BEWlRSbtWajFhY2FWAW1h+6fMB/dd1LRuaXotW1p3Oy+IGlfTNNsjLLY8Cu3j8eYLKKeAiOtmoWTJL2hcWSh52X+6Pm5Y/PhaWFn12mh9jL1ReN2FBZRT33HYcUZxNfjtWISFi+ehfWF2wBQc4YIZUv/DsFDyjbEf49H5ZK0tLZqdtqlbnDktbF2RFp8Om0kLc3KJoUDhGI2L6bvjP14gUrt5usJiy2mMuDBwIVCiEhAhSeYwLy2UnPrH4Ii0UHJ8jesbJYZxxNAH0x7ml5x5cERi6o0egyM2D5ToJnhGyoESlbCr1HeRFrEubRQ/hjsvp7IzBli547y4eLdd9z2+vrBIvBs81mrh/TLto3A9T/XNOF3/TY2TAiJuTMwJR77uSaFWqR5yY9wwtZCUQx9cPRTrmxikZYIjwnt7UeseV7PEuicFaS3Wxz5crNcxy3a/rV3Aw8E81DcukUY+DCIFRFSk0D5X30jSwhyj9tiXjik4Th+MVZGvxb7dHTI5TCq0mxqiUt8s+xfGKNQyKQQjB0oUxo7Xy8JzhiCS0u5+1xjj6dkGR6j3ftR7f46WN+ZeqOXl0P15maRP7L2/W+/9NaeeJQAAAAAAAHCJVb9JJ0nqvf8PSe/RWntTLW/YvYGk65JeIelnJb2o9376W4gAAAAAAADAFXCim3S39d5/XdKvDzITAAAAAAAA4Ioq/borAAAAAAAAgOGd6pt0ktRae/wm4/Tef/+0zwUAAAAAAABcRuWbdK21Pyzpr0r6E5LeXJt9G6+f5LmugtHRXKN7Ld8X0/jSoCblLCSwpthSl37WUvJOSnEN6UchRKXGbZMQw5KS3JqZXkrNqaScxbjWSvvWIl8V089sylnYpindNfUfm/6TkIi2HyKNXCpaJcVV8olmKcU1pZ+5hLLrISppiHTXlGY2NpGXQySOSjl11Jl2f2qfV9JdQ7qY65+2U0p9delxKVGuIm2j2Tic90zbPJ1Tw9juuItp0sVg1krfUk5aIQlW8uln1aQ0l/aVzu2V5LJKiqskjaYu5SwlxIZ0sZT66l6jSXxNfSX5lM6QAHialLN7GiDVvpvjLqeWpuT59TH6wm+PHpJZUxJxOE0G6cBzb1itPnTX+VTbVY67+JYXE1tLQzjF86Frj/VNaN8Lia0u1d61STkR1SW5phTXVLP4uiekycZU+0q6qx/D9R+H/WCS4iCNbdY98/ARN6W+DpHuetBMIn2ol0fhfD0KNVWFS4JOyazzUUjBNe1DfK5I6fXlVPCKIT7DVdqrqfaudko1S0pmtSmztUR6W/cUUuqPe05XD8X6ppJqH/v6oTepcUYDpLvWyobWPkrSN0q6drvp1DMAAAAAAAAArriNb9KtvkH3r7S8QXf75txvSnqZpFvDTw0AAAAAAAC4GirfpPu/Jd3Q8kuZ3yLp7/bef3krswIAAAAAAACukMpNuvfS8gbd83rvn7Cl+QAAAAAAAABXziahD7e94ernt2xjIgAAAAAAAMBVVfkm3f+W9NaSXrqluQymtfYMSR8s6UMkvY2kf9x7/6K7+kwkfY6kT5D0NEm/J+l5kv527/1VZzbX2UItRWit9BBL5dJJJdm4v5YSwEIyq71/m1Lw0jRCNE03E4/JL2HwYRIATVtIbGnFhMKSQdJut9U3bI+Y4hqSSE176pvG3jM7yV54cyeF9pRmlpLLXJJrSkobIt3VpbhK0ti8llGKgypaxJPLuklIp3bJZSnJzaX/StK4n/71uHS2vZDMVmnfC3NLKWzT+fr2CCGT8Riwp/EUPzVQ4l1JKQHUN9sAy+JLrFwfaimzmyetpvbYN6WZpf4mFS2OUUlsDWmmMWlukLj2ddW6x72WmPQXEpjtizSJr5KkVLOk5xykZjGDF/bf1F4+voZIpD/r01PhxeQk2JQ+GfqbtnR9qLTvhSTYSnslxVXydU+lvlm2b54Q6+qbZJt1zzzUBPNwIhqZusC1HdfuJxeaw87n0mdnYZumdldzHzWfSDtE7T+rnNyTapzlEPGXW/ocmE9ZqVbYrK3aHvum8FI3xgAprsu5uOTYUN+EsW0ya6pjTlH3pNdcUfkm3Q+tfr7bqZ91S1prj22tPVfSz0r6EknP0vLm22/c1e+6lq/nSyS9kaSflvR4SZ8t6YWttYfObtYAAAAAAAC46io36b5S0qOS/kZr7cnbmc7JtdbeUdJPSfooST8j6dMk/aHe+9v03v/FXd3/kaT3lvQrkt6s9/4eWn6b7pckPV3SV53VvAEAAAAAAICNb9L13n9N0idJepKW3zb70y39DuUZa609UdJ/0PJG23MlvVvv/Z/33v/A9H0zSZ+6+t+P7b2/QpJWPz9+1f4RrbUnbH/mAAAAAAAAQGFNutbaZ0p6rKQfk/Q+Wq7fdtha+w3F1cpeq/fe3/7Es7y3L5f0xpK+T9LH9H7sLwt/jJav+6d67y+68y967y9qrf28pLeT9HGSvmZL8wUAAAAAAABeqxIc8dV63dKGt39el/QW93hc0xaXiF19M+5jJR1J+ux73KCTpA9Y/XxB+Psf0vIm3QfoTG7S9bxg4Ur6umJ3iwhLkluEPX3pMSxKbcdI04yrM3t2gec4dmp3f3H2X+zcxdrsW1PYfNv8Du0oLs683l7pu2w3i/qmEIywWurI7JQxECG2m8V0CwERy3mYMYbaIQsLGrtgBim8X+E1jtLOZ8Inquy2ju95bX+yfcN77sZIC5FfJlt9iQOMXQo+KF+n3PPVxshhRmGcSl83dnWh5C0FR5TrHnd8hVHSjH0ggu/dQ92TtrUdpro/2eeLr2bzQbZ4jF6BU9wxgRLuGljbIG6M6vUr1Ran5eqp49sr19HN57zVuieeiHzzwjxgEWqkRaoVXI1Z3KaVWneImnsIpc8V5+L3+IZzbs6TNlQh9fXNpddS+owf2tN9jPicZ1X3nP5Nrdykk/xhsetD5eMkjSX9e0kvba19iqQ/quUNxJ+U9Lze+y/d0f8Zq58/F8b7xdXPN9zCXAEAAAAAAIA1lZt0T5R0U9LRBt9WO0vvv/r5c5J+RK+fPvsRkv5ma+0ze+//srX2gJYprpL08jDe7XXsyjfpWmsPb9KtOi4AAAAAAAAulJe0e/z6TO/9wTv/vxIc8cre++E5u0EnSe+0+vm5Wt4Aey9JD0h6pqTvl3RD0te31v6IpDtf/CvCeLdv0r1Ba62SfgsAAAAAAACcSPXXXc+V1toNLW/ISdJvS/qg22mtkn6ytfZBkv6nljfyvkTSJ9zx8LTY0e32kYq/UHz3HdAw5ydJ+r3KuAAAAAAAALhQntZ7//3KAy76N8Ued8d/f+4dN+gkLSNlJX3F6n/fQ8tvzx2t/v+hMObtX4d95Bx+axAAAAAAAACX0No36Vprby3pL2r5jbKv673/wqr9z53miXrv//o0jw9+X9Jcy7m+NPR58ern4yQ9WdLLJL2plmvsObdv0v3mIDO8p/Z6kTa9Em+TbrGaMVIAax+HQdw8xilVLcxjFJLVKvEjsX3z17hNu3hOa4jbyZWQwxh4c/oNsghjuPZK32X7+v6e+s5NX0lamH1vGlJIY8qZO3ir76GZdgxrLqa7uTSytD1Suqt/v1LfMLZ7v8JJIc3Dtc/je17bn2zfMD83xhDHy3mXXmIMzRsi4a1wjcnXXBfNevp5xO1RnZ9NCw5jhGt3m5vG8DsG3fUNT1pKzFWx7ol1iDm+ijWLe2/y9t98jNi/uj/Z5yseYAPEv1WO0StwiovncXsNLG5sN0b1+mWvgTFxdPN6KF7PQzGzcHVIWlWo+1/0simnhTT65Tw2/35KqntSzee3td+maQz3HsR6dIBad4iauyIdL6XLxlBfoTknX8U56/Nkvk6Zv0jLqBWugS3cD4hjp2u3C59N9zHm/rzQ3GschRT3mCa/Sd1z+jfVnQVfoOVNLEn6QElPX/33N+nku3OXNPhNut77rLX2UklPk/RGWqa53m12x38/qtfdpHsHSd9m+r/96ucvDDdTAAAAAAAAIHO3H7uWt/+a1m/itVP82ZZfXP18Vvj7t1z9/M3e+2sk/ejq/z8o9H/26ucPDzA3AAAAAAAA4J7cN+k+VNKna/lLD//kjvYvOJMZ1X2NpD8p6S+01v5e7/3grr9/zurn965+PlfLJNh3aK29a+/9f9zu2Fp7T0lvoeW6dd+13WkDAAAAAAAAS2s36VZr0H2Waf/SM5lRUe/9+1prL5T0bpL+fWvtI3vvr5ak1tpf0/Ibczcl/a1V/xe11r5H0odI+tbW2rv33l/RWntIy1/plaR/2nt/2Vm/FgAAAAAAAFxNfmXOi+cjtfzm2/tLenlr7WclPVXSkyS9WtLH995/+47+nyLpx7X81tyvttZeLOntJD1Gy1+H/fwznDsAAAAAAACuuEtxk673/huttffW8ttyH67lDbdXS/oOSX+j9/4rd/X/3dbau0j6O5I+TNI7Sfp1SV8q6ct670dnNvfJWH0S4tRu96kmb7mkkxSUkhJXTLLKIiSohNAhLUI6i0tcCyFHcWzXnl7LNscorbY4QJLbIEoJSimFqZYSNV+YxCvTJknzkW+fmZ1kFlK9puENGy3WT3kpgbViHhLKUiLXZDRbb7NRi9I0pJy5eY+LKa4VOcV18/SzlGZWaT9YTGzfQ/PeStKRaXdtkjRb+HnMzL56NN+8r+ST9+YhOWoRxli4/jG6e/N0rMFs6XxYfolujMJ1QJLcrtDC+xXHMLtZCBGLR27OITR/M08pbCG5zCXEpoi9MHE39CBh42n/iKnxLrEtvV+V2iQco3ubj5H653ojjGH2yep+7RNsN+8raZCE2K3VPWnnK7xIe55VPi/PF/7odef8dH1I15M9UxfMmu97FI7zcUpQrDDTTtsp1TITM+/UN9Vl42KSa4WrcWLyfCGxdRrqjUrdk+qbSnuqi2exfX0eqaaN9bxpT58TYtp9KU56866D2dL5MJ/b03Vj/cVvdYx068KmWgfhHGITWCWN3EgxeT7cszD7SAtJsCntfpO6p6cHF1yKm3SS1Ht/VMtvwG30Lbje+yskfcbqDwAAAAAAALAz+R9pAQAAAAAAAJwJbtIBAAAAAAAAO8ZNOgAAAAAAAGDHLs2adBfVfH+s+f4JFxcMCyu6W69xkcgwhg1V2KsFR1QWUA5rnx4zhuubxti8vbwIs93Wvm8pOGKbYRJhYdVuFvBMC7zGhZVD/9l8faNMQ0DEKCygPBpgRVi3CHCc83jzRX3zQsmh3YRBxAWUw+seDxB4UZEWDF6EndX1rwZHuEW2D0NwRFoU2YVEPDr3Y6TFuw9M/2kImZimQAnTnhZbTseXW/A2rfOf2iuGWG88DpKuPZUFlNO53Wy/vgjhCWkMt/hx3Khp8eONu6rNNl/IWZJGJiSiheCIFl67O7XE4Ih0uhliR3MK9Y0UapliYJarNyoBWJIUcn7sOJX6JvUvh24NMUalZhngJJKGKGVgFc6T8RAI5+V0Hnfn/Gl4MaPQ7q49SboWOy7UQjomPMnUa3sjf1Ko1EOpjrmodY8Pjjh9YFYKuzqcbx4ccRTmkfaxw5kJnwj1zTQeA5vX3JXPFe4av/wL3zyIAT7DDfFZslIPpZyVeI6z7ZvXapLPa0oZlYsQ7pAO876euaeQoxNrGVcP9bD/nqa+maeJFfBNOgAAAAAAAGDHTvxNutbaDUlPl/Q0LW+zfl/v/WCoiQEAAAAAAABXRfkmXWvtTSR9tqRPkvTAHX/1FpJeclffp0j6Y5K+p/f+mhPPEgAAAAAAALjESjfpWmsfKOnfSrpfr/9LyumXdn9A0ttK+lRJ33CSCQIAAAAAAACX3cZr0rXWnirp27X89lyX9O8k/dV7POzbtbyZ9/4nnSAAAAAAAABw2VW+SfeXJD1W0oGkD+i9/zdJaq19+TGP+fHVz2ecbHqX3/z66N7prjGZNfQ3KWe5r2926Wc5JXXzhFgppZxt3lcK6a4xVc2325SzmKoW2t1rrCSiHdN+aun7reEJfSJi2P4msUmSQkChmkkuazO/UV2qouQTnlLClksFlaT90Xq62FFIcXXpWJI0MWPshSgi13fZf709pbuNQ6yi65/GqErpW5W+LqHMpesu+4b9yRy8lRTX5Rjr/WNSmkkzW/Y3aWsp5Symu67PI6UCzsPxtTD9Y8pZsb1VTkSFxNZ43gvnSRfSmRNYQ7s5ZPI+HVJf3TxiGnpIEZub7WHaJKmNUxJZeM6ZOacWUlyX/V3fMI/KqSXNI0nRb0ap7ikmrbu6x6XRHzd2TGx1yXtD1D1pjDTvwjxSu9t+pSTYMEY1Cdads1KaYTrv2VT7eF72Q4/C8T9zCcwh8S+9cnfeimmyLvpQ/hq4N/KJnvsjP8akrfffC/XNJCW2mvaYal/YF7ZZ96TrxjwU+q7/LKW4piRdczClvimx9cikvqa+aWxXD+W6J9Vw6+2zUPunzxWVzya54PDNQ/C1gu8bwkx9umu17qkk0sfNt3kNZwKfl/0LryXVNzGpfm99oHLd49qrdc8GNc48baCCyggfrOVL/oe3b9Bt4LdXP59SmhUAAAAAAABwhVRu0j119XPTG3SSdPufXu4rPAYAAAAAAAC4Uio36W5/B/oPCo9569XPVxceAwAAAAAAAFwplZt0v7b6+c6Fx/zpux4LAAAAAAAA4C6Vm3T/QculBv9Ka+3+e3VurX2ApI/Wch277z3Z9AAAAAAAAIDLr5Lu+o+1THh9S0n/obX2yb33l9zdqbW2L+n/kvT5Wt6g+w1JxyXAXmnz+8eaXbtHumuQkuZs8tYQ6a6FJJfj200CTTFdzLWnRLRKYmsIfSylu+bXXUjZqSa+2rSa0LeQfpaSCBcpoigk8rhTTUxxXfgJukSzaUg5m4xC2urYJLPOfMqZ6yv5JNeUcpZSx9wYqa9LRJOk8UCJZpuaF2OIXTJrSuNNCWo2GSwcjCnR1yWXxb4pmdXteyGJLI3hElvLKWeFtNCY7poMsDu5ISppZlI416bwrpRs6eaRNkdIFrXpYimxPGxrd1pYhITIcJgf077+nCnlLI5RuG60lBJ3xoaoeyq1yRD1jZRqFj9G3M9suuvmabKSr4dK6fWhf0zSTWOE7iWlxOHwvrhzatjHFuHFTMNTukMmXUfnYX627hnVkj5dkmtMng+1k0ukT31jLWPesFQ7Jeel7klJ9T7ddfO+kq+d0hixZnEJsSmZNaUFm/4upV7KNZWrcebpc0UYw9Y4Mb3eN8dzRWV3KnyOTp/30vXBnVoq9U3S0vXIhzjbpOo4j/A+2rpngPom9R8m3TX1Pfn5Zn54sns7d9r4Jl3v/RWttT8v6XmS/qikX2yt/dgdXf5Za+0pkt5Cy8CIpuXL/pTe+yOnnikAAAAAAABwSVV+3VW99xdIerakn9TyRtz76HX3H58t6W0k7Wt5g+5lkv5s7/0HB5stAAAAAAAAcAmVbtJJUu/9v/XenyXpYyS9QNKvavnty7b680uS/r6kt+69f+eAcwUAAAAAAAAupcqadK+n9/5cSc+VpNbaRNJDkl7Zez8aaG4AAAAAAADAlXDim3R36r1PJf3OEGNdNUf3jdSvl7/QuFRYF7yy2LIUAhGKizDHBYMrYQspDKKwyGZl8eO42HIawy7CHBaaTNupEvSRmAek9XXjOphmsc8eVn1fhKVLe3gTulnYMy22Ph/7sadmIfdxWEB5nBY/NmOnvmmBYrcochojh0Gst7tFlSWpFcbYhbT4sQsGqQZHuEWz0wLFcbcuhD6khYvd4tEuCOK4dhcGkZ4vLchrF1COwRG+ubywcoVdQDn0TQsrh3O+k9aOdqetGGCRghzMdk3BQmZN9WV/c94bYqFkSfb9imMUFs2OCyWXL0pbkuoQV+PEvqHd7SPl8InN2yshE3GMAUK3SiEucYy0ULofoxyOZQc3bel4Sc/njvNQrPVwbKQQrG6Ox1TfzMchGMDVLKkmiIFZm9c9aTO5kIhy3WPesMtU30i+xkljVOqe2LdSO6WaJdU9bh6pbi/VPYX6RlJ3Y6e6J2yPFk/6m5+IUlf78TqeU9N+7c5DhefzQ8S6p6Vzu9musWYJQYEuJCLVJulFVsKu4thDhEGcou6ZhnN6xelHAAAAAAAAAHAqr/134dbaMyS944BjNy1vAl6X9Jje+5cOODYAAAAAAABwadz5yxt/RtIXbel5uiRu0gEAAAAAAADG3SusnJNFRwAAAAAAAICr486bdD8o6eAe/e+X9FdXP18h6XmSXippuvr765KeruW38q5L+teS/pakVw82YwAAAAAAAOCSee1Nut77j0r60dSxtTaW9F8k3SfpmyV9Wu/9MPR9nKR/I+nPSXqk9/6ZA875Upnd36QbJ/sCYyl0pJKIlvoPkZQWxqkmpbnIkzhGJcE2JN7kpNr1lJhyyqzd1gMkW1VTztxfzHzPnhI2TZqh5LdTSrBMKWcu7auFRLmYGGb6u7Zjx3DJrKFvGsO1h6C/KI29LSnNLHG7QkzBK6SfxaDKNIZJfkqpqin0yY2RjoGYXObaK30ln2gWU1yrKWcbtkn5HFJJqg7nyYU7zsN1KiaXmbTV1Ddl3Lt0sZzMmtrdBqmNUeq/xTSzcxK2mPenLaXdV9NJt1n3+ITYzZ8vjlGtWUxCYbX+qtSY+aRvhohJq2EM2zldp1KaYWg3F/WWEjZTHWISVHPNYptt6mu17rEBltQ9a7ZV95RTZm3dUxvD1SFxHimZ1dZOhfpG8rVMMS00tleEfc9eT8qxnO5kVjuXuXNwSnFNh5GtK4ZIpC/WPZVE+nLdUzimT3O6mY1P/8upld3o0yS9p6QX9t4/Md2gk6Te+6skfaikn5f06a21Z59qlgAAAAAAAMAlVrlJ93Fa3q/8yk06997nkv6Blv8A8xnlmQEAAAAAAABXROUm3dNXP3+p8Jj/tfr5zMJjAAAAAAAAgCulcpPuxurn2xYe85arn08uPAYAAAAAAAC4Uio36X5t9fNPFR7zJ1c/f7fwGAAAAAAAAOBKCblm1vMlvb2kj2qt/UTv/cuP69xa+0RJz9FyHbsfPPEML7mjB5oW950wAWSAdNfEJ44W+hbnUk5sGySFzSSUpeerpLANkYJbTbyyiVKhb0pFKowRk5XmKd3VtKXYnLCt54V017T9bDJY7Ovb3caOXSvRQAOklsU5D6ASEJkHKaaLDTK2aSscA2ns8rHh+seEssIbWUzeqrzGKndObeHoiO+tPbeHtMCU2GYqm5hEFt/H0/WN/YvHUSnlLI5R3BcqthW2OMS5bIiaZRd1TyUptZxqv/k1sJRUW66dBthx3JOGpNWY+lo5RlOCZXothRpzEdNWTTxjjGYMze68XKhvknLdc4VrHP98hWtjZf8Ng8Sk2hSCadPJi6mq9hgt9JV88Gm5vvHNpV2yep485RgtnMtivVGoe4ZIjS9t0y3WPfk9HKDuuUDprl8u6ZWr//6y1toPtNber7V2/+0OrbVJa+3dWmvPlfT1q+ZXS/riU88UAAAAAAAAuKQ2/iZd7/0PWmsfIen7JF2T9H6rP2qtvVzSkaQ31OtuXTYt70E+p/f+6wPOGQAAAAAAALhUKt+kU+/9RyS9o5Y36todf54k6Smr8W63/YqkP9V7f/5w0wUAAAAAAAAun8qadJKk3vsvSfrg1toHSPpYSe8h6Q0kXZf0Ckk/K+nfS/r63vtswLkCAAAAAAAAl1L5Jt1tvff/KOk/DjiXK2n6gNTvv3e/kgEWZx1gDfFBgi2GWOB5iLUjSwEW1QV5h9jWhYVV40KzLhAhLToa57x5/xgcUdge1f3UPmP1eNliYEvJAMf5IAY5WaSxt9V3gDnvYkH/yoK8xddY2Z3S0K45LqqeVgCvzKOyuHBcHXvz5xtk8enqYXu22TPZNo/zii2+yEHqjQH6l4O0Cn39NfD0zzdE3VPdxWxgViV0R2Gh9HKQWai1SsFRvnmQS9UA8zh13yrqm5P3lbZb4xT62hpniGtgse4ZIoApKgTEpPAOey5Lv+tYOpyLdc8Q9cYZBTPccx4VWzjOZ6XfVfUGGAIAAAAAAADAaXCTDgAAAAAAANixjX/dtbX2tcWxm6SxlmvVPbb3/mHFxwMAAAAAAABXQmVNuk/XyX57uJ3wcQAAAAAAAMCVUA2OqKys1yXdkvRKSa8uPg8AAAAAAABwZVRu0v2xDfqMJD1N0p+U9OGSXizpw3rvv12f2tUwe3Ch/oCLlwJOqZxEeLbpfe2Mnw/AQAZJIvQnIr52f17xzuCccSnO5d10e3UINQ5wAW0xaZm657wa/h2Yj05/b2fjm3S99/+6YdcfkfSvWmsfJelbJT2/tfZevffpCeYHAAAAAAAAXHpbS3ftvX+HpH8u6ZmS/vy2ngcAAAAAAAC46LZ2k27lO7T84ugnbPl5AAAAAAAAgAtr2zfpjlY/n77l5wEAAAAAAAAurG3fpHv26udky88DAAAAAAAAXFiVdNeS1trHSPqbWkZmvHhbz3PRLR6YqT04G3bQISJohgg6qYyxzTnHlNPCc1bGjn3985XSyArziOMW5lceI3HzK/StjCsNtE0H6F9PmisYJmrq9Lb4IksvcYC+gzzfAKeVythxztX3pTT2EGOkgzf0P/U8tvh8FQO8L4P0rdrqyaxgm+e9rV43zrguq9Q9Q1znh6h7yjVcYYxkiLpngFPLEO/BIPUa9c3ru4j1TbF/+Tpf6Gu33gCvZai6wo4zRG03xPwGmUfx5HTWNc4gzzfEGMMf5wud/t7OxjfpWmufd68ukh4r6Q9p+Q26p67auqSvPukEAQAAAAAAgMuu8k26f6Dav8fcvrf5b3vv/7bwOAAAAAAAAOBKqa5J1wp/XiLpsyR9zFCTBQAAAAAAAC6jyjfpPnGDPkeSXiXpl3rvLznJhAAAAAAAAICrZuObdL33b97mRAAAAAAAAICramvprtjM/gNTjR88OtFj+xbTU7vpH58ujB3nZ9rd8x37pG5+i9oYtn96vsrYKUxnkebh+tbeLxtME8ZI86gkxFbGiP1T6ltlftV5DJAGVxljqymz2xojibtkYV+tpqoOkabl+qaFHoZILQ1j21DFUXjDCmO02DdMML0Hbi7F98Vuk+prLM2jkFyW5hHfx/X+aZOmsZubX3pbiq+xFLwXt9PpTxhx3qdUqm/iIKk22XzSO6l7wrXbjr3NuqdQl8WaZeEHd/MoJ75X5lE5BEJdEQ+X1N+9xuoYQ6Tgpjpp0+dLqqeVC1njFE9wW0pVLYfdDpDWbuukSl/J7gzV+qsyj3RBSs/ZzLW7uq1tzVKpb6RaumtljGLd47ZHNTH3rOueWN8khf6nqW/mfXryB69U0l3/5uo/v7r3/uoNH/MkSe8o6Xd77z9bnx4AAAAAAABw+VW+Sfe3tfw3jG+RtNFNOklPlfQDkn5G0jvVpgYAAAAAAABcDdV016rZ6ufTtvw8AAAAAAAAwIW17Zt0H7H6WVgRAQAAAAAAALha4q+7ttY+XNLbmr/6rNbaK48Zs0l6rKRnSXovLX9F9mdOM8nL7MH7DzR54GT3SiuLH4e1e+MYrn2R+qbnXPjX5RY/Tn0XaQFl0x7nVxijz2sLOdv2MEaPC2GaeQywIG95rXAXYDEPY6RQith/s7Zye3GB562NUV3I2b3B21v7dBDDLFxcCzPY1uLHQyxcnMYotY9qixz3cWER5nEYo7TgcnEnqyygbF7Lch6mPfVNi0GP1w9SuyCycvCGW4y4jfzJYhTGdu1pkeM0Ruo/dvNL15jCySLskvXFmbck1SypxqmM4drnxfCJWLO4mir0zXXP+s4aX0sMqnJ1TzgI0kY1NU46r7j6ZvWIjcY9lnvS6vXc1CyxBgnzi8EMW6p7hgju2mroFnXPXWMMUPcUwqQk2WvjVseI7eZcEWuCNA9XE6S+oT1+Bg79Xd/Kdor1TRjD9Y+hW+E6b8bIdU+oQ2zdU6tZRqZOSoeA6ysdk5kxQN2TapzKGJuYzg9O/NjbjluT7g0l/a272pqkzyqM37S8TH1xbVoAAAAAAADA1XHcV7i+W8ubbLf/3NYKf35S0of13n9w6IkDAAAAAAAAl0X8Jl3v/Tdba+8laaLlDbf/rOUXRj9e0m8dM2aX9HJJv9V7f82AcwUAAAAAAAAupeN+3VW99x+7/d/tdb9Q/OO995dsc1IAAAAAAADAVXLsTbq7fJOW35J7eEtzAQAAAAAAAK6kjW/S9d6fs82JXFUP3bip/ftOlh6S0kydmEQWMlRmhRQx11fKqWhz038e0sxmcx/h4177PCSUxQQ1038xDmOkxNZZIZk3jeESAFPyVsy8MQqJY5JPORvNiimulXTX0HcU2t28K32X/QvbupASt82kNJuIdkz/M5cSUV2cU0hKSqeyrSWzpuTTQtpqHKPQvignlJn5pTHSDpKu+oV9MoVd+1S6WkKZSzlzqWXL9pAMZvqnFLFxGHtsxk5j7KWxzWscx76hPSaUbT626ytJI/MG50S083HCiSnuLj01nBTSGK42iX0LY6T2WPfEtHszRqh75qHecGOkYzSlvtqXHq7FOal+84TYdKqw12Lf9ZiaZf0RQ9Q3kq9PymNX0mdDGm9ljCESYuOpwr1fl6m+Cf17jBYOY7t6o5AsmsZYjNMB5pttfVLpK1/3xHTXlEht5r0IO0h5tyl8hIv1qztBVVJcQ3vbC9fzMIarT0ahRkrJrLbuCa97b+xPWqW6J42dUl9N/9TX1TdpjOQ0dc/R0a0TP/a1z3/qEQAAAAAAAACcytq/qbfW3lrSX9Ty3+e/rvf+C6v2P3eaJ+q9/+vTPB4AAAAAAAC4rNwvvrxA0puu/vsDJT199d+316Q7iS6Jm3QAAAAAAACA4W7Sdb3ut6nv/vvCglgAAAAAAAAANuFu0n2opE/X8tdd/8kd7V9wJjMCAAAAAAAArpi1m3SrNeg+y7R/6ZnM6Ip50vWbun4jRVQeL6WLLUx0Tkooq6S7zkJ8zyxE9eTU1/X2aUhxnYYxXP9ZSCJLY89NAs1sllKYQtqaaesxbsk329im6ndWzRgx5SwmkZ0+5Ww027x/pe+y3SQ2bjEpzSXBpv45Ka0wRiERbdn/bOPPYppZ6l9IvMrpriYBsJBmJvlksNR3EVKzbEJZSDNb7IX0Sdc/JK3GhFgzvWoQdEpQG+R78m6MmEqXEls3TzlzSWSSTzTbS2OkxFaTXLa/F9LMwsE7cWOEk1acR2o3J7Pcd/PU11Hom17jWRui7qnULLNwMMb6JrQfmfop1iZh3kczU/eMwhghqd7VOItQ34SAWLnMuXg5SlHQQ5xv3HPG5NPQPkht4tvdOLlvIZk1jVFIn3X1VOobx0hpspV013gBC/M4L3VPIak+J9IX0uRjzbL5GCkVNNUbzdY9/vlcX0lamBon1k6+2RqFbbdIO1+qD+25tpbGa9sL6fWSr3FGKd01XedN/1QjudpEkvZM/9R3Eubh+o/DgV6tWfbMSa6a7npWdc/BwcGJH3sb6a4AAAAAAADAjoV/x193R7rr83rvNzd8zP2SniLp4d7775xgfgAAAAAAAMClV/km3TdJ+kZJTy485u0l/aKk7y48BgAAAAAAALhStv3rro9o+dvab7nl5wEAAAAAAAAurI1/3fWEnr36Odny81xYT9x/RPddmx7bJ4Y+hPa5ufeaxwhhC6Y9BUccupVBlRdnPjKLJR+N/RhpYeXD0Xr/w7CA8jgs4HmUVjo1ZmEhYZlFVBdhMd24eKxbWLm4qLJd2zKtfZpWZ7WLC/uu1cCGkdnFKwslL/tv1rYcIywYahdh3ryv5BdLjgtKh33Bbuu42HJhceb0fFUmsCGphD7EMIPwfN0sspv7prHX28KpTK0Q+hBOe/E9t4ESm+8eq8FNU8odCtspLbw9xHrc3e2UaQHlcD5sZuHiURjDBURI0mSyvlHS4scpDMItRnxtz59wJuFksW+DI/wYKVDCLZQsSddMe5pHWkDZjT0Oe59bbDmp9JVyfVLp6+qeVLOkwKyp6X+YArNCuwuIWLavnzBcLZTmIflFrGfhGHAhE5LUTN2TKlC/qLq0MM/Z00k1HLsuwKZVCx8XRBDrnnA+LIRJVeobKdUsp697quETtu5JoQ/brHtscEQxfMLZYt0TT00p9MGGCPghKrVM6pvDrsxnk1QjpTAIU+Ms9jY/npd/YeYRuqay0wb0pfohDJ7qG7efxb1pgOCIFmsZE2ZQCIiQpImpZWJ9E+qha6b92jjUPammMictVwtJOSDiWqiTXM0yicFYm9cy1bpnkxrn1v7hPfvcS7xJ11r7Y5Keav7qw1trv3vMmE3SYyU9S9Kf1XJ//7mTTxEAAAAAAAC43I77Jt27SvqSu9qapH9QGL9peZPuHxfnBQAAAAAAAFwZx61J911a3mS7/ee2Vvjz25I+q/f+bYPPHAAAAAAAALgk4jfpeu+/0lr7WEn7q6Zv1PJbcV8g6bhfd+2SXi7ptyT9Yu/9aJipAgAAAAAAAJfTscERvfdvv/3frbVvXP3nd/beX7LNSQEAAAAAAABXSSXd9b9o+S25R7c0lyvpcXs39UBIjLuXeYi3ccll0xB9mMZw/VOK641i6qtLYT2Yh4QyFy8kqZlkFde2HCOkLRZSfVLKWXeDpDTDkAa1NfHFhCSnQuJoJYks9Y99Q3slKa2UcpbSXQspbDkpLbS7pLRCEqykYeI4K9L+m5LLTP+YZpbGMKljKYksjm3SyFKKa0oGs6eymJjtx3Dnp1Fxm9og6JSYW0hxliQXdp1TztKL3LBNOeWsmdeeUlzHod0ls6aUs5Q65pJcU8rZ9bGPd3SJrWmMlPqaUs5c+yScmCupr+OUiBb2htT/tCr1TepfSXGVanVPancprpJ0OF9/D45CIv3BfGLbXaLc4axSxvtaZpHOqanuMdspHc/uOrB8gGtLZ5xwvnaN1fPeAKn2Oe3eXHsK9c2y3Vw3Ut9CLZPmUallcn2z+dgpbfzC1j0uITamuG5ey7ike6mWSO+u8VJOqrc1TgruTe1miJjiGmsF0xZTXFNBdNb7TWhPdY9LHA3v+V4hsXU/3F+IdY+pT66nmqVQD7k0+mX75imuqX+qe1Jt4vpvo77Z39tiuuvdeu9//NTPBgAAAAAAAGDNccERAAAAAAAAAM5A7XvyK621+yS9j6Q3lvQG8jf7mqSxpBuSntx7/8QTzhEAAAAAAAC41Eo36Vpr90v6IkmfIumxxef6xGJ/AAAAAAAA4ErY+CZda20s6T9JejflJRGdQ0m/VZwXAAAAAAAAcGVUvkn3aZLeffXfvyjpX0p6RNLXaZn38gGSpqsxnynp81b//+ze+y8ONeHL5jHjR/VgSFm5bV5cOtAllC0KKa6pfdJ84tg0JJSNQmLryCS2jkOql0szS3ISWUiDM8mAi4XfTnOTFihJzcQUuTYphkH6RLMBkmBjAGNqdylnhUS049pd6lg5Qc2mnKU0s83HaDEJNqQwuTEGSHdNsaA5/cwlpfmuQ4j7b0xs3TzlTCnlbF5IItzzx24p1CvHmZqmYoqzmd48ndpjsl1od9MYYl8YYozyicgc5yERLaa7mvZxOIdPwvXXJaKmFNeUfuaSyG6EMa6NUvsQ6a5+DNd/HN70UTi5j+POejqp7km1zLyQaj8NqfG+7gnbNLQfFtJxK/WNJC0KyazzcMJ2x0bqG+shW7MUj/PKv/UPcB4qberUt5gc62qZnBC7eS2T0+tD+3Tz+qvNwv7rargUZT5Iqn2qh3z3IdjDINbzIXF4W+muIcU17Xu9kEhfKyzCuSKVPebavQjbo/K5IvYtnoZKu1PpM1zoGlPt11/QKNQs6fPyxJzbJ2EMlzwv+Vrm+p6vTVIt45LqUx2TxkjX1+umTqpeo10tk/omm9Q9o/FRaUyncpPuY1Y/f1rSe/feb7bWHtLyJp0kvaj3/urVf/9wa+35kn5U0ve21t6x9/6aU88WAAAAAAAAuIQqX9F6Ky1vOn9N7/3mqu01et2N6Cfe2bn3/suSvlzSm0r6S6ebJgAAAAAAAHB5VW7S3b/6+dO3G3rvU0kvXf3vu5rH/OfVz48szwwAAAAAAAC4Iio36V6y+nnfXe0/tfr5/uYxh6ufb1aZFAAAAAAAAHCVVNak+1VJbyPpzSX99zvaXyDpwyT9udbaN/Te7/y7j1r9fPlpJnmZ3Rgd6b47FpJcFO6bzsMCym4BxLSA8qiHRWJDuzMOfd0ix0laDHovrA66Zxa9HIcAizQ/t5B4WqhzNPLzc2uoxsyH6kLz2zLEAsrF9looRVp0eIAxzILGlYCIZX/zpGmM2L4+RlwQeRFeZGG99hg+EbhFkVsMOQiLAJsx4kLJadXh8fqT5qyFtEHWH5EOudjuFrw2oRbScdtps3GX7ZsvrBwXRN7mAspJ5VSWtrVrC5OOp1rTPy22nK4xbmHlvRQ+ERZhvmYWYR4iIEIaZgHlSVsfO4UZ7Ju+yagYJlGpe47siui+3kghGOMURFKoe6oWZm91bZI0C/WQ2yfnIdwh7e/u2KgcR7F/tYzZUtkTr6OF82E9MCuFHxTGSEEONmwh9A2HqBu7EhCx7O+ugSFQrRKklS5UMTDLjDtAfSOFXTKcmlr4TGA/9pg6JneWmqmT0jmyUpal2qRUbwxQswz1ueLUfYdii5Z0/g1DmPYUmBUDJUx7CsxK7e6zdfwcHtpdzVKpY6Rcs7j6qV73uMCsFNKxeU1wd90zC6+tovJNuh/Rclf803e1f7ukV6zG+u7W2he21j6ptfavtEx47brjV2QBAAAAAAAAvL7KTbp/IekRSe/XWnvc7cZViMTnaHkD7wmSvkjS10v6hFXbTNIXDzRfAAAAAAAA4NLZ+CZd7/1hSe8r6Q1776+66+++WdJnarkGXbvjz6skfXLv/X8ONF8AAAAAAADg0qmsSafe+4uO+buva619m6QPlvQULYMm/nPvnfXoAAAAAAAAgGOUbtLdS+/9DyR9y5BjAgAAAAAAAJfdoDfpUDdu/fVSuMYy6V0pJTWlY5lAmJQilppd2mpKSkkpsyklzrWPQoJKbndjnH2sT0o/O3NDTGObqa8DzMOmsBXSZOPYKa2qki5WHMPOu5jiascuppwl9rQQ0udCcLSX3pcUeWVfY0gRiyFxJpWumHzaCknVQ+y/l0oxxbFyTq1dYzZPu4zPF97cSntlzlJOF3PXxnS9dCmu6TlTimtKbE0pohWVuifN78iUsxP5vvOwU45Me97+tffR7gvF/clJ++8Qx0b1OX3n0tAXUul6Ujxc/HWjWLMUroG5tnPX4tA3pd0WahabBJsU6564nVwifUqeD+dDm/paqCUlnz7bQtJnDKQ258+YCrzF/dc6Tymu56QGO+vPkufl83m61qX7DS6xPfcNNYtpj9f5U9Q9Q9RGleAIAAAAAAAAAFvwev/02Fq7uaXn6b33B7Y0NgAAAAAAAHCh3f37ATe0/PLn0F9QPydfKAUAAAAAAADOn7tv0v1vcUMNAAAAAAAAOFOvd5Ou9/60XU0EAAAAAAAAuKpId92xeW85vXVlEfI9UqrqvJAH4lJcl2Oszyk9X0yOjc+53j/NI7efj8iwfk7mMYjKS9nmy05Bn6Y9bf6w2/gQ0ZEfJKVm2ZStkAAWdw+T+tNSwmlM9TLpXeHpqly6WJkbI2zr7e5PJrEtvL7cXnk+31zZfy+VLX43f4jrQOUcnq51lfY059Qer7uuPSWwp5rAnORi35gmn+IFT6da97jtV6mFlmOs968837Htbl8o7k/OEDXIVuupK/C7ObEOce3FTe2vG7XrlO8+wBuTXku6ztu00FDfFGqQrdY96RRSeA/ia4k1X2GMQk2V6uJt7r8l5+UzyA6cl8+SZ/35PF3r0nXXJcqmeaTEVlsXpNPhKeqee93b2QTprgAAAAAAAMCOcZMOAAAAAAAA2DFu0gEAAAAAAAA7xk06AAAAAAAAYMcIjtixwz7RQd8/tk9cUDrcY7WLH4fVD6d97NsX67vGYfe7y8FiYtsPzRiSNDPPOQsLP8b2hRljERaaLrSnxTtz+2ZtxxpiYeVtLbgagw9Ce1qY1q3TGRexDdvaBCjEMVIYhOlvdiVJ0mjh3xi3oPlIYRFRv4671NxEwo4Q2t0QvbzzBYVFm/u4EMIQ35fQbsZOz7fY2/w9D6e90v6U+/p2e8xUj69K+MQu1iGu7H6hrzvXxvNvGNqd29P1YRL2J3ftSWNMm9+h9szCxXuh7ygcu6k9nXJsV3eykF+EOdUV40JAxMgE4xynstB0mp+rZdKC0qnucbVMrJGq7eZC49qkvJ+5fTItvF2ph9K7VaqHdlH3VAyxOP4Q5+vidcNfe0LYQrqumUM3XS/zNzfM38zDmxiO/+4Cs8I+1tLYbtwt1j3x1FSoZVJ9o1TLjM0xGt6YRaq/CvOo7JO5Pt98jBzGUWw/bd+T9HcKu1/aVf1nydrnUXduT59/4zXG1Cezkb/258/t62Mcps83SdpH3FRC33TfY2w2tquFln19+yY1zuEApyW+SQcAAAAAAADsGDfpAAAAAAAAgB07s5t0rTX/O5EAAAAAAADAFbfxmnSttWdLerak95L0hpKeKOk1kl4q6X9L+hVJ/7r3/pIwxP/dWvsQSV/fe//qU80aAAAAAAAAuETueZOutfaxkv6mpKff/VeSHivpqVreuJOkL2ytfa+kL+u9/z939X8XSW8r6V1PNWMAAAAAAADgkok36VprNyR9s6T/43bTXV3+QNK+pPvvfJikD5b0J1trf6X3/o/v+Ltnrn6+4FQzvmRuzq+rzV+X7prSSJyUXDY37SlxLCWDudSxlNaaxk79j0z7o3P/29BHcz/2USEpbVpIOUuJN4tFIRosRqUNMEZB5emkYgLrAO2VJDKpln60iBtwvX8K3kqJiC51rI9C3/Aam0k6imlmIQ6qGKB4ajlFNPyF3Z9On+6a0szS4g2LyeZjxIRY8z7GlNm0XxfGKCfHur7F498aJPmsNhF3TKfz7yKe29dPIuncfjjz16mROYdUU0udlGQ6G9fSQq+NZhv3nTQfMz02iWapb5JS0U4r1TeJe+2uFkp9U3uljpGkw3loN/0PQt809tTUQ2n/zemu7viq1T0+3XWIE05QGKI6Dddcrm/Sedxc51PIYYvb2oxRPTG7oPWp7zoP1+KRqU9SzeLqm9QeQ6P3LmbdYxPpq6mqhVT7Stp9OK0cM/Zmbce2b+lzRey7xdNQ6bNdnMjmY6Tzbzq3j83xlT7/jkOh6mqc0fz01/h5+qyW6qEw7+lofd6TvnkNJ/kax9VCUq2+Gd/1fDfnp19Rzo6wWj/uu7S8QddWf35R0qdLemdJj+u9P9R7f1DSkyW9u6T/U8tfeW2SxpK+qrX2xavx3kbS4yTNJX3fqWcNAAAAAAAAXCLpNt8XSvoALW+4/ZKkD+29v23v/Z/13n+69/6a2x177y/vvf+P3vs/kfTWkj5Kr7tZ9/mttWdJeu9V9x/rvb9qS68FAAAAAAAAuJDWbtK11p4q6XO1/FLmCyW9e+99o19R7UvfKek9JP2kljfqvkLLm3Rd0vMGmjcAAAAAAABwabhv0v1ZSde0TG79yN77q6uD9t5fqeU36m5pecPuo7W8SfdvTz5VAAAAAAAA4HJyN+k+RMsbat/ce//Nkw7ce/81Sd+i5bfpRlr+quvLTjoeAAAAAAAAcFm5SIynrX4OEfDwHyR92uq/v2OA8S6dV07v18H02rF95nHpQM+lpcR015SsYvqn9NTD0D4L7S6ZNSalheQyO0ZKgk1jzEwaXEhjSe0uFa2nJNjthOBlMZUqxYuZRKlCEqyUE578flaL6eqF+YWw1VLKbJv7wX0yaxgjppy5yRXTzNw8tph8FtOxUnKZC7xKUbrx/XIpZ75vSmz173lKRNt87JiUlsZwaWvVpDSXtpYuD6m9GDrmO6dBXARz6JrSDM2xkfqm8/I0xSqfUkpVnIUTjrsGHoUU15Qi6lJcJenR0fpJZxKSyPZMX0kam4tSSrBN7W6MIaS6J9Usrj2NkWqTqTmYUt+c+hrqEJPkmvrGWsvUMrnuCWOb/uk4GuIYzWmGrq0WweiGSGHj+Rpj2orJ8ymx1b2elDiYMpVd3TMabV7DSf41puulS3GVpG7qoVTfpFOCrZ0K9c2x/Qdgd79CfSPJ19GF+kZK1/nN+0r+/a3WTu4Ul+ub1L75PEqJreXPN6H/toRjIJ0n3WfJeUpPDq/xaLOZHcteR0PS6qz72sQls+6P/A7i6hhJuhbaJ6Z9L3wASzWLS3dNfZNN6p5b08PSmI7b8m+w+vkbpx5deukd//39A4wHAAAAAAAAXDruJt0frH4+foDx71/9POi9/8oA4wEAAAAAAACXjrtJd3sduncYYPzbY/zOAGMBAAAAAAAAl5K7SfeDWv729nMGGP+TtFw24n8OMBYAAAAAAABwKbmV/P6tpM+X9M6ttb/ce/+qkwzcWvs0Se+m5U265598ipfby6cP6MbR9WP7xIWSCythpoWIc3DE+v3b1DeFPrjFIyW/CGVcQDksiuzGnqZFxMMYbrHkWVhsOQZHmMV0e+irysLKxcVx3VuTFlAuhT7U1gWO3FMu0iLHMQxifTJxIdwQ+uDWIl3E0IfUvvkCyq6vFNa2TevdhkCJyj6SXksSwwjs4GGMwgLKcQwXHFFYAFwKCygXg09cewyOSIswu8Wgi4swlxY5D4uLV9+DimYO3rTYeloY3J0/w+Fc0sOJJbW769Q0XNMmI/8iJ2MT7mACBCRpL4yxH4IjXP+0+PFeOAGMTPs4BUcUTjjVRZhTbWH7hh11bsZYhJ3d1TdpHqmOqdY9rgZLtUnaz2am/2yAumc+qwVm2RonHufpGjjECWe9KV8fwj7prg+p7gnn5XRKtZsk1D0tzc+d88MJMQURuTokBUQs0vtlXkwMjki1jBvjEtU3qX85gG2AusfVIUMEVcWArlgPmTGqdY+rnYrhE/l9DP1t50J7Ct1JnyXtwVELjbTPF9rToWvrntHmIZCSD4kYh/qmWve42iLXNyE4IoRS2DHCFtykxnl0Otn4efLz36X3/pOSvlPLXfrLWmufUR20tfZJkr5Gy/3j1yQ995TzBAAAAAAAAC6tdJv2r0j6bS2/afc1rbUXtdY+qrXjb+u21j64tfajkr5+9diFpM/sPeT0AgAAAAAAALC/7qre+2+11j5Q0g9LeqKkd5L0bZL+YWvthZJ+XdJLJB1IejNJbyHpHSW95WqIpuW36D679/4DW5w/AAAAAAAAcOGF3wCXeu//q7X2DEnfLOn9V81vJOn/OGa827+E/UpJf773/oJBZgkAAAAAAABcYsf++mrv/Xd77x8o6Y9pGf5wS8sbcenP/yvpr0t6GjfoAAAAAAAAgM3Eb9Ldqff+XyX919banqRnavnrrU+UdJ+W35r7fUk/0Xv/jW1N9LL6/cMHdO3gxmv/v5Jylrgx0rgp2a6S7urSYCSftpb6pzFScpkbOyWR5WRW8xrDPFJaqE3ZKqacDZHuuumwkuKteZt+FpKckhiy5YYJ6UwhSMhu6xaTWVNymekb3tuYGGaGLqW4KswvvmG1sa3Ud4jUzDSGTTnbPBFNCmmmaf9N83BBhDEZLBznLl1sgITYakLZNlPOhtgXSueymH7mWsO5PU3DXQPD86X2qUk/HIeT055JcV323zzlLKWqpvQzN05KHEsJZc30j2MUE1u3Jabdu2Th1Dfs7G6MWJsMUveE+isls5ox0v5bq3vCsRhSX33dU6hvjmuvcE9ZTNJ0dUgKHE3SS3erePvzW65DXOJlte5xYYZp38s1i2urvec2yTUcL4Ocbgaoe8ofyWzicDi+KvtqNSG2UrNU6p5Qb6SaqlR/VWqnYnr9IFevAdJdc6q9qVnih7XwfpnnnI99/TAfh1Rw0z/WJmFsV5tU657Uv5Lu6uqbNMYQ9c3dYxwebnSL7VilEVYBEC9c/QEAAAAAAAAwgGN/3fW0Wmvv01r7d62179zm8wAAAAAAAAAX2em/i3e8N5b0ZzTQN00BAAAAAACAy2ir36QDAAAAAAAAcG/cpAMAAAAAAAB2bNu/7op7ePnN+zWZ3Heix1aSYGPKWfhFZNc/jZHSzOJzmgSalKqaUrbcGD0ls6bkMtceEwdTUo9pT31japaNwfSdE5dMExKbEpugFNNxwnaK6VHrbSk9tYezku0fxohpq65/Ic1MCqljsW86wDZ/bwZJcR1CcZe0LzEdGpXU12rKmetf2E9Te6WvJJ/6Vk53XX+Dyymz8T3YvG/aJ93uno7FyikuHkbhxbvkt5TQPQ+JrSOTLtZCotzIpLhKIT01jhESymyrTyNLaWY55Wzzvsm2Ul+rSfc20bdQ3+QxQt8wj5gQ72qWQo0k+Zol1k7p+mprp7CtK3VPvOZuXg/F80olkTr1TemTZiI5rbWWBO/SXdOnrjiGTVVNfVOysOlbrHsqyd2Vuqd8+jgnNU7p9FSte1y9UU0ttimzxTEqdU9hfqk2SXXZwtQ98StG1fbKDpiSiM0xkz47KSRmuyKnhwvYIrWb81MLKd/zUIc0lxqf+qbdtzRGSHct1DKV+ua4/n6Mk59wpjdP/z04vkkHAAAAAAAA7Bg36QAAAAAAAIAd4yYdAAAAAAAAsGPcpAMAAAAAAAB2bG0J09ba1w44/tMHHOtSevXNGxqPTxYckdZmrQ2y+QLK5aeLY5u2wuLCaezyGDawodA39a8EREgxdKCkkj2RXopZzNUufCxpERYBrYQwVBaUTu2lRY4V1ogt7tiVMVppdfzaPLa6gHLFIAson5MxKoswxxX9w9h2Iee044Qx3NjFhZzj+zVEAIA9psMCwGEi9roWF7b3Y9vFqtPrC9tvbk5+aaFkVRZFjosZh7HTAsqbPt9xCv3j/M5Yqe4ZIHwiPl2hvkkDpeCIGGBlrz3FusL1r/SVwrW4OI/Uv6IQmBXfF3dejkEQ1bqnEJRQCaoaIvQhiPMrjVGs7U7bd5uqu+kQ4RMDZMkN8pnAXRvLY2wedlWZR/m1VOuhCvf5JtU38RxiU5x811D3+LCQEDIRw3EKiT6FMar1TVKqewaoaU9T98xvnv75Xc7Qp+v8nB4BAAAAAACASy+EgQ9ybxkAAAAAAADABtxNuj921pMAAAAAAAAArrK1m3S99/+6i4kAAAAAAAAAVxXprgAAAAAAAMCOpTXpLozW2rMk/XlJf0TSW62af07St/be/2V4zETS50j6BElPk/R7kp4n6W/33l+17Tnf6ejhfY11bftPNEgS7ABjSCeIKXJjnL6vTaAaImGzkEi7nEdhjMLQKYksbv7C9nBJsMf1t2lw5RdZeL6gFPZzxklkQ4RrnndDHPrnJW2t+nylt3eI1Ldy0mehbyHsK6duh9TSuTuZhXmk6C27nU6/Par7b+093+L7tc0xzrNd1D2DnOTS2Bu2HaNU9wyQ9LnNuqeSZpi4hN2Ual9J7k1/EfeOQWra048xWH+DGmdDA9QbzhAJsdVrhn3LB6izynVPNfXVDZHSk00dks4VeYwN2477i0rdU6lTB3nPBzr4L1DdM394euoxTnyTrrX2DRt2nUn69N5TqPnJtdbeTtL/WP3vo5L+l6Q3lvS+kt63tfb+vfePuesx1yX9R0nvLek1kn5a0ttJ+mxJH9Rae8/e+yuGnisAAAAAAACQrN2ka629r6QvlvTc3vtXHvPY52jzf1v5N5L+S3l299B7//nW2o9I+kpJ39t7n0pSa+3TJX2dpI9urT239/7ddzzsH2l5g+5XJL1H7/0VrbWHJP13Lb+J91WSPn7ouQIAAAAAAACJ+1L3n5T0bpI+csMx2j3+SNL7nG6aWe/9j/fen3/7Bt2q7Z9I+qHV//6p1060tTeT9Kmr//3Y29+YW/28fWPuI1prT9jWfAEAAAAAAIC7uZt0z9LyG3LP33CMPyrpTcOf79LyRt07n2KOJ/VTq59vckfbx2j57cGf6r2/6M7Oq///eUnXJH3cmcwQAAAAAAAAkL9J92arn7+84Ri/2Xt/qfsj6d+t+rzlqWdad/s5//cdbR+w+vmC8JgfuqsfAAAAAAAAsHUuOOLxq5+/M8D4v7H6+QYDjLWx1tozJL3f6n+/6Y6/esbq58+Fh/7i6ucbbmNezvjhscYxKhND2GbY2rkxQBKhW2JyJ2FcV+H9wuV1BRLsainOydke6C0831VIHMTVc6nqngGSGTPqHuDULtN1dIDkY+oe6JHT39tx36TbX/181alHl16++vngAGPdU2vtKa21T5P0nyTdL+nLe+//ZfV3D+h1NyBfHob4g9XPM7tJBwAAAAAAALhv0v2+ljepHm/+7k6fufqZbnhJ0vXVz1lxXiWttb+sZcLrbY9I+rje+7+5o+3OG4WvCEPdvkn3Bq21Ue99UZzHw5t0q4wJAAAAAACAC+clrR3/Hcbe++t9qe24m3RPk/Rjxwz0dRtM6Cmrn39wbK/T+w0t15N7kqQ3l/SApC9ure313r951Wd6R//0HcTb7SNdri/vAgAAAAAA4BxzN+leJOmPSHpfSf/G/H3F+6x+/saxvU6p9/48Sc+TpNbafZKeI+nLJX1Ta+2h3vtXaPntuSMtf533oTDU7W8PPtJ7L9+ku/sOqNNae5Kk36uODQAAAAAAgAvjab333688wN2k+0+SPlnSh7fWPqf3/uqTzKS1NpH0Z7X8Rtr/c5IxTqL3fkvS17bWDiT9C0l/t7X2Tb33V7bWXibpTSU9MTz89k2639z+TJf2XjPS3twtDTiw6i/ZFvqXF+813/asj7Fh23FjV8bY6jw2X3W0sp162q3St20LryXPb4Cxk8pKp6X3vHg/vjLvIX65/aKu8DrEqt7upZcX9a0cNJV5FF9f4bW0NHZle8T2wkLC1UOjNL8wD7fIROG0Evtv9bUUxkgq70sco9B3m2PswiDn2s27lq7n5TE2b99q7VQYY4h6Y4jXUn5fLmLdU91vNh332P6nfL6h5nFebKu+GaKv5OdXrRU2HffY/ps/307qnsKiVnFXLdSHtTHCPMIQtW19+jESXzvV9put1j1ndMppN09/b8eN8H1a/nrqYyR94SnG/kxJf3j13997inFO6jskzSXdJ+mdV20vW/18h/CYt1/9/IUtzgsAAAAAAAB4PWs36XrvD0v6Ci1v1n52a+3jqoO21p4t6e9reb/yJ28nrJ6xW3pdYMXtbwz+6OrnB4XHPHv184e3NSkAAAAAAADgbum7eF8u6ee1vFH3Ta21v9taux76vlZrbdxa+yuSXiBpImkh6fOGmmzRh0i6tvrvF61+Pnf18x1aa+96Z+fW2ntKegst1637rjOZIQAAAAAAAKBwk261rtuHSXrlqs9fl/TrrbUvbq390dbaQ621Jkmttce31t6rtfaFkn5V0j/UMpxBkv5G7/1HtjX51tqfbq19fGvtwbvaP0jSP1v979f23l++el0vkvQ9q/Zvba09tOr/kKRvWrX/0977ywQAAAAAAACcERccIUnqvf9qa+3dJD1f0ttKerKkL1j9kaR++z7dXQ9tWv6a6ef23r968Bm/vs+R9N6rufxvSb8j6a30ugCIF0j63Lse8ymSflzLb839amvtxZLeTss1+H5U0udvec4AAAAAAADA64k36aTX3qh7V0mfreUNscff8dcpquM/SPpbvfefGmSGx/uU1Z9nSXpzLQMhflXSD0n6l73377v7Ab33322tvYukv6PltwXfSdKvS/pSSV/Wez86g3m/1v5rpL3ZvftZA6RMlgJXqmOEX6buzTxggGTWlC6WU8cqY4Q4GNO/8nzLMdb/Ij7fEKmllddYTWxL/e3YfoxWeO2pr9vFlu3r/V3b8i82HzuNMcg8CtLzDaEXp9cLJ5fU1z1npe/yL9b7L2KKWBjbJIDF17dIY5j2kCxm+6ax0+4bx/APcM9ZPDTs/FJ6Wh7D9E1jxPbNt/UQ80hc/1Kq2jHtlfTZC5mUtotE+lJaaGoPx38lETWw/WOdFQaxNUs4JxRqqrSP2XovjZHmUTldF2q15TwK6alx7M2fM9Y3sZYxtUKxLhvZhNghapbN+ybVuuc81Th+DFPPF2sWP0Z6wjT2APMo1BulmiXWSPEkYsbwXat1j92hqtfiAeoNW1sUX6Mdo/J81TEK26lSZ1XHToYY4zT1zejWyR9727E36aTX/urr32utfaWk95P0gZKeruU36yaSfl/L1NT/IukHeu8vOf20NtN7/2VJf+0Ej3uFpM9Y/QEAAAAAAAB26p436W5b3az796s/AAAAAAAAAAZS/LI9AAAAAAAAgKFxkw4AAAAAAADYMW7SAQAAAAAAADu28Zp02I79h7sm05PFh2w1mbWSRJbGjmlaLs20Noa7vbwYb9439c/PF9KIzBiLcXg/w/xs7xhhVUg/S9s0zc+li8XXEpLBUrsZexTSzEZjH/fj+o9Gvu84jDE2cT9pjPCWa2z6p4Qym6oW5pGkMYZIg62opLVKPkE1parG4CfTf77wO/Y8JYaVxvDtCzO2a5OkxTyM4dLW5mHOYQyboBbHsM3xfG2bi+95IURQLczPtbfwGktjVOfhklkLfWN77BuO80o623lJiN2i7SazmjGK6amuvpFUSoKP7YWaxfWVQsrsePP6Ztm+vjOkvqnucftZ3MUKaau5vgljuP6FOkaSWqg3XD2Urue1uifUFakeKtQs48LYsTaxrb5/GiM5L3VPTIg35sUxKjVLCi1dmP5pHvNUs5gxynWP6Z+SYFM9ZGucVPfE+mbzZNv4zqbBK4n0oYYYmdeT+uYawo1bG8P2r9Q3oT3WN2GMUirtAOmzyWlON6OD05+r+CYdAAAAAAAAsGPcpAMAAAAAAAB2jJt0AAAAAAAAwI5xkw4AAAAAAADYMW7SAQAAAAAAADtGuuuOXXvNQpPDFG+yUgvY8wE0xSQymwBWTSKLiagmkauYUOaSWUfFhLJm9v74fOlIMeEto5TkFCNoNhv3mGYbQdNDSpdNcZWkvfX2trd5apmUU1VdctleGDsllO2N12OHJuH59sIYEzNG7BtikVz/UXhn9sIYLtFsL8QcVdLPRjEqqWYRY45d37S/b57uOovJZevtszC36dwfvK5/SkqbpnYz9iykmaXXMputj5FS1eZz/54vZuv9e0yCLsZxV3adGMdrni0los78PEazwhgpKW2QMcz5sJruavqntLVq+plPUPN9bYSdBkp33VbY4hCJ9LFvqBWGSGCtpMmn9NRCPbTYqyazujFC37A/uXTGlDIZE+ldY/VrA4VU+5hUb+qQVN+MUs2S0uRt3bN5XSFJe2YMV8dI0qSQ7prqjTS265/qjfRabLpr2D8uat3j6o1c9/iD1NVOqa6YhjFc/1gjxbrH1F9hjFRTzUzNkpJg56OQam8+S/YW3itz7ZfiJTBcM0OafOHamOuKzZPqK/VN6l/pu2w3x2ixdrK1yTbrnvjebv5BupziukH/0b3u7WyAb9IBAAAAAAAAO8ZNOgAAAAAAAGDHuEkHAAAAAAAA7Bg36QAAAAAAAIAdIzhix/ZfPdX+/vS1/58XATcqCyunEIfYXugbFygOC8rbBZQ37yv5xZIXJvhAknrYy9vcjeH7xiAH056Wimxh+7kVlwvr164GN21pjNDuFksehW06SgsXT8Kiw2bx4/2wgPL+nl/pdN8857Xx5n0lad+sorofVkWNY5sxJnEBZT/22OwlabFl11eqLaw8Li6sPC8sfjwPO5Trn8aYhpPIzLQfzv1BOg0HzZE5qI/C4scH84kf25yIDmd+HmnsI/P+poWcp9OwOLNtC4stp5WS4wmqcO0J76NtLYQqpPa0+PFo6tub6e+CII4b248R5hGCPnyARW2h5DRvuyhyXGw5jGHmMkRAREv7Xhp6S3VPqitSLeMOpdQ3BjaEYCY3l1RvLFI9ZPqna3Qc28w77XsLfzq0u1l6C1Mt01wdGJM+CqFbaQcOYRCu7hmngIhQs6S6x4UwxLon1CzXTD2UQq2uj/0J0Y3t6hgp1zKuHtoLJ/FJaLfBEeH9qtQ9Q9Q3SbXucSEMuW+4/pv5uTpGyvXQkRk7jZFqFjd2rG9SPTTaPHRLoe5ZmJPFIlxjXKjN8i98cymYsRA+0cI8SmFXxbrHjRHrhzR2pXYaJHwi1Sahv6llcshEqoHPpu5pR+GNKuCbdAAAAAAAAMCOcZMOAAAAAAAA2DFu0gEAAAAAAAA7xk06AAAAAAAAYMe4SQcAAAAAAADsGOmuO7b38JH2JofHd6oknymkkZXTXdfb+15IEYxJZCkN1qScTWpjuJCimGYWxrbJLzFdLIztxg1DtHBL3AZb1sLxamLK2fqrGYXE0ZRmlpLLrk3WI4OuhxTXlKp6Y289Kee6aZOkGyHlzLWnlLP7Rke23SWrXQtRSZOYoLY+RkpEK6WcxejO06ukuErS1MQqpxTXSvtBiBw8TOlnpv3R+b7t+2hIdz0wKWcHY9/30ZlvH4/Wx2jT0x/oaYRFSBeL53wXbxVOZvEsOUDKmWsvp5yZ7VpJRItjFBNi3RgpzWw0Dcd5TIM1Y4f5uTSz5VzWnzOnnBUizocQL6Rp/zVt43DRTceAqTdiumuqeyab10m57kmJrZuPUUnHS+fwVA65tyC+XSlh1+2T/jIQ52HPWSlNNqTujkySazXFNSXSXy/UPa6+kXxi6/VUI8W6Z72WuRZOWqn9ujmBppql0j4JJ/dq6usQXI1TTXd1NUusb1LN4mqnkAR7a+FrGZfkmuqbWPeYWibVN4fhs4Krew6ntVsPU/MehE2a013ntc92fvDN2yv1jRTS5GN949vHlbon1J6ldNfUXql7UvpsqntmpmZJ9U0l3TX1TTaoexZT/xmygm/SAQAAAAAAADvGTToAAAAAAABgx7hJBwAAAAAAAOwYN+kAAAAAAACAHeMmHQAAAAAAALBjpLvu2OjhA41CauprFVPObHLZqJbM6vrHdNfQvkj9TfrZYpbGCOln++vt85QiGENYXP/QOW5rM0JIHWohlMq1l3OI3APSCw/tzaS+jk3ymSTtmSRYSZqEVDSXaHbfxCff3Le3efsDoe+NkMzq0s9Siut9Y5+6fL2ZtLUQoeT6SrWUs3HYccZmXx1tMflsEf5N5yjs8K5/6psSW11C7K3FNT+G6StJt+br/W+FFLyUcnZztj7GXkgzS6l0qb2im5SzmOKaUhVnxetJgQ2ILZz3JJ9+ltJTU5qpOxwriWiSND5ySWmbp5kt201itkknk6QW0l1dmpkUkllDX8WEWNPfjCvp3Ke7yiS59lD3KB0bpmZJCbGpbnP1jeTroVj3pNRXc3oKb20pqd6mpEoapRRcM+2UXh+P8yF2m0Ld49LrJZ9gvxfqnkqKqyTdN1k/6aQU11T33G+SWe/f87VJSnd1NU6ub1La/Xr/VLOkemjfnNxTzeL6Htd/CK5mmYdqfO4OAtVS7Q96SFU1B3rq6+obyae+PjoK6a5jnxD7yGi9fW/k35fxzI/hpMTc2G6Ox9RXKeE8fh7a/DyZetrPcMW6x23WVPekmsUddq6OWfYNdY9LiE1jxPrL1CYpvT7VN5V2V8coJ8Ta/ind9RR1z2h2cM8+9xzj1CMAAAAAAAAAOBVu0gEAAAAAAAA7xk06AAAAAAAAYMe4SQcAAAAAAADsGMERO9ZuPqoWggbuKS2K7NrD4sctBUrsrU+qmbZl3zD2JCwoP1/f7dIC4G0/vEa3aGtxIWK37mgPi4iORn5wtzh7WO82v892AdTii7ELKIeu4bW4he3HYbHlvbF/kfuh/dp4fQXU62mR47CA8mPMYslpAeUHx37BTreA8gOh7/1moWQpBUf4OafgCNfuwiQkHxCR2ocIJ0jSQr1H4d96FuYYTYsfH402X1g5B0eE0Afz3jwyv277Ptx8u9uu1W3ttt9s4bfdPLTPzKL+o3AOD0Pkf5or5EbEl+7ai+s4uwWU0zk1Lazs2isBEZJfLHl8FBY/Tu0uOOLIv5g2DS8yBk2YFzkPY6cwiCGCI8KizaeW6ptCcESsb1I9NF4/D/VJKJNjcIQ/l7V9U1PFQI9C3VNltl8aNoVguWM05PbkUsaFzFQvX25XSLtNLJfNcR5CgSahHqrUPam+eTDUMq7GeSCEPqS6x9U4lfpG8sERqW8KfXD9R2El/f0QEHHWNU4MjgjtLiQiBkTE4Ij1EIZKfSP5GueWCYKQpL1wzh8VPoek+tDVMrnu2XyMeQrdSZ/VYupDaK8wg1eCsVJ7CmYIeS02DCIGRGyx7nHhDqMjP+lU98QQrNl6/2baJOXaxO3vqe8p6ps2f/TEj72Nb9IBAAAAAAAAO8ZNOgAAAAAAAGDHuEkHAAAAAAAA7Bg36QAAAAAAAIAd4yYdAAAAAAAAsGOku+5Yv3VLPaTRvNYoxnT6dpdcZlLLlkOElLM9s2uEdFeF9LOUXOZCm/q1k0bc3vmEIW0tbD8bEBum0QtJPaXkw9R/iACrYrqrax8VU84mKeVsz6W7+rSfB0L62Y3xenslzUySHhyttz849gk8Kf3s/rY+j5zi6l/jNbPj7Ic4qEl4H93ePsBRJEly72LKOEqH+dTsgAfdb6dpiBe82dfTyNK2vtV96uvIzHwck+M2T3JKaWaz8FpcollKOZvO/TvpUgdT8nQ8EW0xHc8+XTp3pk1t2mMKZkxKM+eykJSW0s98MmtKOfMTce2tmHIml+KqkGiWUs5mfoxuU87CGCE51loU97FU4zihlpFJiG4pxdXVN5KvcdLrTqmvIQ1u5NJxU8ER2FNOrHv8GAvzlG3ut3+u4QpphlsK/42K571KuuteaHcprpJPsL8RUu1dfSP5JNfH7t2yfV19I/lk1seEGum+FlJfR+vzdrWQJE3Cm37dnLAnKb3+nNc94bKhI3PgHYa65yBEIh+Y7epqIUmahBrT1TiV+kbyNc4ifLCYuROLpJk5xx2lvrEeWp935XOMdEy6a8UAn+Hi6cl9lgxvl0vXlnyqvWuTjkt9LdQ9KZn1cL09pbvG+ia02xon1Tcx3dWMEVPtU+rrvWucviDdFQAAAAAAALjwuEkHAAAAAAAA7Bg36QAAAAAAAIAd4yYdAAAAAAAAsGPcpAMAAAAAAAB2jHTXHVscHGrR7oidKaSctZRy1kzyVujbCwllrU/807nUMuWAUndnOIT6qJnXIknNRD+lJLKU6rcwLz2lBbo0M0n2RZbTzLYVtlhMcWymf9obU/rZfogd2jcRQykRzfWVpPtM+tl9I58ultpdYmtKcU1JaT7d1b/u1H7NbNj9sK9PwrswMv3HqW/495hFzC5bNw876jSkH03NeSGl3R6E6K2xSVBK+1h8KYXot2lIW5ua/WkazqkuzUySDufrY6e0wNQ+tgnMIeXMttrLg6SBTkOVpOqUcubSXYvJkW4XqSTBSv660WYp5cy3u8TWcorroU8GdIlmfRr6huSy7lLRUpppSErr4fp/Wunar5BIb1NfU4prSINrk/Uap5njVsq7dazgzOtx53ApJ9KPTN3Txyn5OIxh3t4Uglc67orH+bbqnrTbpPfFp7v6yVXP1+5aldNdfburZdJ11KW4Sj7JNaW4PibUPe457wu12vVQe7pPEJNU94TU4lTjONW6x9U4i3B+m4bXeGTqof1QI01iu0nHTCmTwcKkzM7T9ghR0K7GmYb6JtXzLsl1f+TPqUcmoVvyx677vHKsbabaF4aOnw/dGOn8m9LuTb2cPv/GMUx7qntSLWOTXFN9c5TqGz+2rXHC9Tw9Z3c1TrqPcYpU+0X359gKvkkHAAAAAAAA7Bg36QAAAAAAAIAd4yYdAAAAAAAAsGPcpAMAAAAAAAB2jOCIHeuzmXoLCyfeFhZQjQs8uwWUw6KImy/BWl8ouYUFnvt4vT0uTDkO7TMTHLEXQibMYprLdjNGYTHzZX/zgLAI/jbXLS1JCyub9lGYdF6UPizCbtr3wkadhBXeXXvqe92EO0jS9dH6sXY9HH/72jwMIgVE3Bc2lFss+Xrzi+amxY8npv+odERnC7eAcljFdpICJcw2OYyLH/sx5jHJZd2ReW8laWFSaQ6aD8G5FsY4MIESk0VYNDe8lj2ziHg6BlIwiw13ScfouTnhBOlc68J4Cn2lED6RroGFxZnTYstp5X17XQsLIrfQnhZF7q5/6hsWUHY1RA8hEwrHbnfX1+Ii57bGSSFaoR5q5j3o6X0JgVndjB2vdaaOWY4Rrv+mHup7YX7z8BrdPlmpTVSseyrt5/x0UzlPpr7pvJzO4+6cn64Pse5xoVvhOrUf6yFT94QxUijFNVv3+NdyPRwD18zxNQnJTuMwhquHtln3TEMdOAn10Mj0H8eDY/P6ZmrCJKRcv05NOMNRCrUK9VCl5k7HgKv947FoW9MYofMwu8KpxfIrnlMr5/bQbkO3isERpmZpoSZIn9tt8GSsbwoBEZKvcVJARKqH3L2TcD/F1jfLv/Dtd3bpoe4q4Jt0AAAAAAAAwI5xkw4AAAAAAADYMW7SAQAAAAAAADvGTToAAAAAAABgx7hJBwAAAAAAAOwY6a671rvuGYsVUn16SEVqLjEoROHYlBNJzSXExpSTYrtNg0tJZKG9kABYSREsJ5SdcaJZeo3bmkY1ITKlwbr2USENSvKpaGmMcRxjvX9K3kopVm4ek5AolRLKJiblrJLiuuy/PvY4JUEXuWech9Rihe20MP3T7FL7vnl/p6H3OKQtufc87jchsc3tI+O4/xaTLQvSsbE18YRzxhFqxXN7bYzNB0mJmWG38WOnMWKqanpOc9wVr9Hdzi+kuKY0+cL2i0yN01NKuj8d2v75/Sps0x7OTun92iu8B8Vk1o3HlWrH6BDH1xDS853xeW+b59lcsxSuPSlZNF2TbN1Tq7/ctbh6PR/bmsXvp3sx9XWL3y0x575Uf03DLjIx228RtvUk7PBH5phO9eg4HP/2PY/1edpvXM1dHMPW/rXj66yT6s+6zJI0SC1T+UxbSe5On8OH+IxfGiO0x9okpLvbVPtKTbCpAWojvkkHAAAAAAAA7Bg36QAAAAAAAIAd4yYdAAAAAAAAsGPcpAMAAAAAAAB2jJt0AAAAAAAAwI6R7noRhPSj3L9w73VU6ZvSHbeY9JdSaStPGfraMc44tLDqrEMVd8GlglbFJNLKGGFnWJj2eQpsDNFKC9O+COlY7vkkaRTSz4YwNylnbs7LvqF9gGhA9x5U94/Flv4tap5SHy+TmOh7ttM4L+flHq5HcXquf7peVq+vlevuENfoWFeY81Y11cxupwGOr3LtdD7er7SfnbnzcpzHc/4uYhi3Y5vXk8o1MF1f53afrG1/VxOMUjppTBZ3qaXb3HapLgvt5vXEGimGhZsac4CDcYi6eIj6/Lw79y9xi59p/Rhn/xk/Xl9t35TAHpJZ7fki9U0b+2yuPVfgUwYAAAAAAABwvnGTDgAAAAAAANgxbtIBAAAAAAAAO8ZNOgAAAAAAAGDHCI7YsbY3Vmv3eBvSoqhhYcXmFjqcTHzfcRjbtLe0OGMYoxfGTgs/9vHmizCnvnERZtMcFwyttJ/3RUcDtw5mDxskLZmZFpV17bOFDz5IY0z7en/XJuWFkqd9/Vg7CmNMQvvULBg8DaEP47ClbLtZEPk4bhHm8RZ3vrT48TQsuDo1r+cgjHEYdqipWUz7qPDeSv79TX3T/uTa877u5+faUyhI5ThKx+hFddahQD0GAGw+hgrXnngdjdfXEBCzMOeLiT+Q4iZ1x2O6Xs7D+Wlb+1+oe2LNYrZfS9tuEmou1z+NUXwfbR0S6xvfbGuWyv6rLS4ufkG582fl/Csddx435/x4fTh93TPENfCo+ev5xMz7MC22nivE9aZiJsjYHEeznuZR42qcFBDh6htJmpqJT8PzTcOLd3VPes9T/ere31wXn77uqdTzaYy017hj9IzW7R/eAJ8lK59T0+ffIT5bV+4JpPsBbRHa55vfV2jhnNoV7nu4IzJdRxebB9isz2smze7Z7Vh8kw4AAAAAAADYMW7SAQAAAAAAADvGTToAAAAAAABgx7hJBwAAAAAAAOwYN+kAAAAAAACAHSPddcfa/r5a279Hp5A6EpJVbPpZTCjz7W3PtO/7pJQ4dkhQ63sm+cW0STkRpu+ZBJqQzhKCi2z/OEYhGTD2rSbHbksKqxkg5WweUnZmJsHHJVhJ0uHC72eHi/WYnNT3Zrtm20cmqWsckllDqJf/p43Qdx6S0ubmOSfNvzH74Q0bmWS18UCxey45Nm2OnPq67iAm2Pl94aZJKLu18O/tze7Powd9fR85CPuN65v6Hy78+S21z8xrdMeFJM1j++YpZzGZ+axT0QZIzA67R6l9ERLKRiO/QVz/lCza98JGnaxffFp4A3poL22+lOQW2m366TwkJab2be1QlTlLIZk11Eh7ofR17amOSQmxrnaS1M2+UK177D4cjoG0v7tjo3p8lZKPz7i+GeJ86M6zy/bN6xvJn/Or1w137bk28tepySIks7b12mmc0gkrdU8wD4MsTF0xCcmsqS5zp+uhUu0rdc9RqntM80GqdcOHk5umDrkZ6p5bPdRDi/V6KI2Ra+71fTL1TTWcOzbS54RBUu13kfpqU1V919pnyVArpM+p7tw+Lo7hPluH81u6fsn0byFNOqanhtsibtY9JsGHmsVdX1N6/SKl2t97R2t9QborAAAAAAAAcNFxkw4AAAAAAADYMW7SAQAAAAAAADvGTToAAAAAAABgx7hJBwAAAAAAAOwY6a471q5fUxtdv0enEAcTEk1kUltaSIK1aWaSTSjrIbUsp59tnnK22Pd9F5OQAmTaFyaV5th285Qx5Sy89FLK2VmLaWabp5zFFNeQ9jOd+w11NF7fR45CEtlhbHdJnz4CaBSSwcYxq2vdPBwzC/NvG6nvfkguOzRpa5Mw59Q+NjFWo4GirdwzzsOOnZK6pmY7pQSwg5D85NJWY8pZSn017anvrXkae30/y4lo4Rgwx8Y09J2nlDNz3Lk26ZiQs4ESgK3C+bCSKFlOd3Xn9pB2uZiECbqEsnT6CElfI5NctghlV4vX+ZBy6tqraaaz9fNQTJk9J+muKU3ebo/0usN2cmOnFNfcnuoeV7NsXt+k9rT/pv3dHhsDpCeXkmCPaz+tVN+E7pVzajovp/O4O+cfjnzfdD25ZdIFJ+FYdDVBsghvWKpl3PV/2lyGu3Td1DeSdOBS7UNNNg7Jluel7pmGDwWuxkn1zZH8GK4+SXWPS3FNY6R6ObX7dNfwWgrt6XNC+lzh011t12NSX8/JB7NKqn343JnOte401MLn3xYu5829B+nwChO00yvWN5r6c4i7l9Fm4cWY+kaSujmntpQym1KwN0l3XXTp1j27HYtv0gEAAAAAAAA7xk06AAAAAAAAYMe4SQcAAAAAAADsGDfpAAAAAAAAgB0jOGLH2vXrauMbdzQUFrc0ARHL9vV7rz0tlJwCJUz/HhZsriyUnPrHhZL30wLK6689LaCcgiP84uK2a2mhzvIizEOET7g1LONin6HZLqDsJzKbh7CAsJ8dztdPNXthFfbUnsIgKtxiydOwqG9uX3/TXcCBJF0PCytPbHBEWgw6BUqExVLtGLWFld1iyXGh6biwsgkLCQdY2n6u/SAusL15oERaKPnhuQ/xuTkzCznP/RgH8/BaTHtaQDkurGxOFukY7aF9oDW2Pbf4cXHheLsofVwE37+YhenfQt/KwsrpGIjM9byNw/ktvEaldhfklBZQNgslS1Iz+2QLfeNCyW7B5WqYhKt7Un2TaqQh6h4bmLV5HSNJi1QPmXGqdU83+2olGGvZbhZhT3VPOu7c9KrrtW8tOCI0h/OhO3/GgIh0vh75Y+ZgtH587YW+KSihYlEIOUjBB5VrdKpvUvu+Oam6WkjKNYurA4eob5K0nebhOy5+W28ejCX5GqcamOVqnFT3PGLqG0l6xARpPZrqm1kI0jK1/zQERKTPFe4YdZ9Xln/hm8vtTuEcV/0caEMfikGGrj1dH9LHqWbPfQPUPeF6nkItc9ilqXHm4cWEzwrN9U91TwqUcO6qe9r89Od0vkkHAAAAAAAA7Bg36QAAAAAAAIAd4yYdAAAAAAAAsGPcpAMAAAAAAAB2jJt0AAAAAAAAwI6R7rpj/b4b6ns3ju9UTTkz7SmlK6Wcuf4pKS2nn4WkHjNOTGZNY5TSXW2z5m6M0Nelqkmyt7ljAGBMOVtPgBkkgLGYctZNqs88pC3Nwr6Q0s8OzWscFdK7JGnhki3Dxk6JXFOTUOwS2CTpYOx3BpemdW2UUs6ObLtPOQvprint1qS+pr5V80KKZUoucwlqMeUspDBVknRvmSQySTo0z3krJLPeDGO4RLObMz/GrdDuUs4Ow/GSUs7c8RiP52q6ayXEqpLkWEx9dLteSqpsqYIxh4FPLavJiW3hGmiaU4pruka3kBY6mpkX6dqUE1u7SR8bJN11CMW6x74HMZ00JceZ2iTVN2ns0N+l7FXS61O7q2OOHcMcM6nuScedTWCupjhXkqCTQqp95TwZ656w3xyO/IYaj9Y37GiLKa6uRpJC3ROuuQfjkFS/WK9xUs1yPdRDLsk1JbOm1Fdnm3XPInyXJdVIrmbJibmhVjDvzUGonWLdYw5ql9YqSY+mesjUMgemjpGkR1O668yku6a6JyS22ron7OvbTLWPQ7hjOl6nwhjunJqSu9M1155UQ6p9PNmabR3OWbmGW/+LUUqZ3Qs1S/js38y+E1NUUy3j+qf65hR1T5+R7goAAAAAAABceNykAwAAAAAAAHaMm3QAAAAAAADAjnGTDgAAAAAAANgxbtIBAAAAAAAAO0a6644tHrimxeT6yR6c0mNcwEtI40wJai6dJSXNuL6STzNL46S+sb2QzJrSz1xYUny+kHLmnjMnANbaS2zKWYjeiSlnpmtKOZuF9NSQAjREotnMJGSllLNZ2Kgu8erG2CeR3Vr4xCuXXJZSzlLqq+s/RMrZENs5Sclx8/BvPS7lrJrG6xLo3Ht4fPv6GC6t9bgxbpnksoM0Rkg/OzApZ0fpOArpZy7lLL0v6TjPkVy++bTieS8lR5rzkGuT8vk6JnJZKUXUDFtIM5N8AugiJH6Nwrm2hf59bjZKSDlrIYnMtqe+MRV4SztOob6RZGuZ+L6k5Fj3fsVU+1QrDFH3+On5uiecl/3pSS7UM9Zq6Rg17bHvNr8KMEDd486fKd11GlJc2wB1TzqPu8TL6Z6fR7r2PDper2VSbZISQK+N1uuQnOLq66GRSWGNqfYuolu7qXGcWqp9eL8qqfahb6pZXI1zVOgr+aR6l9Yq5XTXI1PLTEPdkz5XuM8hqSbYRd1j066Ln/fsW5NeY0pKNdditz+uOttm+9k/1nDp2rg+9v/X3p3HWVaUBx//PTPsiKKIIEpEMSjBHVzQ1903UTRxCzEo7rsx7tGYxd1ogkqIGjXGhSi+IYs77qAmbgQVd0VQorKohCWyD8w87x91rn2mp6qnD3O6T8/07/v59Ofee07dutXn1j23+uk69Wy4ut6ONY1s6M0xS2U81CrbHMtUxk5LMb7ZcFX93DaEM+kkSZIkSZKkiRmkkyRJkiRJkiZmkE6SJEmSJEmamEE6SZIkSZIkaWImjpjY1btuz5od6gvU/9ri1zIFICsLLjcX7x2SOKK5MPDAxBGVepqLLTcTNiw+cUR7e6WOxmLLjTVi6wsoN491Y3vtVx/4ng/SWAcza4vSN9q8vrFA+bpovGEVrUV6mwsoVw5sbVFlgCvW1N/IHdZuuvjxpWtaSR/qCRt2rNSxXWPx49aiyNut2XTh0qELJa+tLMK8ZqlW/6edpGN9o8PX3sfWAsqt7VdXMrZc1Xi91qLItcW01zUywbSSQdQSObQW6b6ykfShliRiXWMR5tYCyusr23N9KxFMawH1+uaovF+D1+iuLqC8+AWKy/bKttZ3T+tcVvldNrQ+GwOSQQxa9BmonULWNN6vDa2kD421h6Nyvm4uttz61Wuv2SrbPNj1zVusuf73gEQfrXHFkPd84NikVb52yhmasGGMcU8tgcWQ8Q0MXSi90UEGTBFoL+pd+V0an4HW+bB2/lxPK0FEox0DNBMwNcYyte2t77qdGkmwaokBauMYgB0a457tK+ON7Rpjp9a4p7Z9yPgGln+MM2R8A/VEWq1xanPcU9neHG80+kJtjLOuUUdrPFR7zVbiiFqCCKgniWglzGolbNlQOX61v1fKjvrmUbQ+/7XtQxMI1v6WbJ3bG82oNaSZaKWVmGnAub31/bWmMq5ojWPWNBJjNccylXoGJ46oHMDm+KZlEcWvXrduWJ0VzqSTJEmSJEmSJmaQTpIkSZIkSZqYQTpJkiRJkiRpYgbpJEmSJEmSpIkZpJMkSZIkSZImZnbXiV2963bEjtfwbWhmKBtQdoQML0Oy9ME42V3rdTTKNjLQ1Mo3s+m0sp8NyUTYzHJY2Tg0i1it6lYKoFaWs0r5DY1MhMSWx/fXN7OcNbK7VjI8Xbmm/sZsv7aeSmiHSjayVtntGtnFahliW2VbmZW2r7SjlbWsVceaxmsulQ3NtMV1tWy8rfe8lsUV6tnWWlnEahnRSt2btuOqxuvVsriW8pUMe63MrI2sY7W6W9nM1l/dynJWyWbYynLW6h7N1JaN8jWNPpmV96uZEXHI98Mim7WQGPw9tenv2PqeaiQ5rGY/a2Uia2V9bX3Ma/U0TwmtBG8j1LHsBox7WmOTMepoZrxrfv/XMscOrKNyymmOnQZkbB2cIbbWviHZ62l0p6HppGvFW+e3VtbXWlb7xoegdV5uqWaZbmT6vHpt47tnu03Lb7++3r4r1zbGQ5UTVC3DPLQz1e9QGSe1xiytrK9rK+9va+zUslLGPa3srrUxS2u80aqjNnZqZWatlYX6eKOVxbU17qmOnQaMb6A+xmlmr29ld618N7az1w/Laj9s3NOoYkBWe1rn60rxZtb4+ub6r9j8/mqcDwdkIW9mbK1lYG1lqW+8X83srrXMrK0x0oCM9Esx7rn6yi0PsTmTTpIkSZIkSZqYQTpJkiRJkiRpYgbpJEmSJEmSpIkZpJMkSZIkSZImZpBOkiRJkiRJmpjZXSd21S5rYKderHRoVs+KasKgZnbXxvZaRpgh2WQHbm9nRG1sr2TIaWZEG7B9aDtq2XeG1jHGe16ttvHm1rK4AlDJvtPK/ruhmZqx9ZqV92u7ekNa2c+uqmQ0a2Yoa2RKu6JSvpZxbKG611a2tzKUxYC6W1lcW1qZ1aplB2ZEG5LJtZbNrNRRec9bGdFamX4r7ahlHANY39peqbtVtpW5rNYnm+1oZWytZRFsvl7jc1TJLpitDMyt7a1uM0b2zspLtrNgNrKctr6UanUMaEcrIXVre+24RisL5voBmcgGZ3FtbK8cp1ZGtGZ219qxbpatbx/yek0DvgOb3aP2qwwoC40MrGONe2rjjQHjm1b5IeOb1vYxxk6tz/OSjnuq2V0bZRufu/r72ziHt5rR+l6rjXua2evrH/RatvA1jbHJ9msXP2ZpjXtqZaE+Zlk7MKt9bXuz7Aof37Tr3rR8a6zQ6je1jK2tMUtzHDJg3NPKSF+ruzU+H5KpvlVHLYsrNMY4rfFNK1vogHHFYAPGPbQyele2NbO4tpKW174fWllVW98x1cysix/fQP19HDy+af2S1cysjfY1v7xrZetFB5lXx1Vrt3wenDPpJEmSJEmSpIkZpJMkSZIkSZImZpBOkiRJkiRJmphBOkmSJEmSJGliBukkSZIkSZKkiZnddWJX7RKw8zXLOjMoWc3QLGdDstUMyRDbqmdoBrUhmdKGZJkd8HrNdgzMclZNUDa0S1TfsGGZ/qrphSqZmUrV9Uo2NFLHZiXLTTOjZ+M111Qyl129ppEZrLl98VnOWplZ11bqbpcdlv2s2o7G9qHZYLdU6/1qqbWuVUczO141Q9mwOmoZzZpZZht117KRNcs2spzV2tfKZtZK3zUoy1njdxySNatpyPdJ63ujdXqqnj9b6czqm6u/YuMcvqaVsrGWmbWViWxAxtbmx7aZ5ayVuayysfWeD8noO/C0slSnocHfgSNktR+SIbY1vhmS+HzweKNWvpk9ub69+jsOzEhf/Yy2PudDxo2Dxz2VKlrnzlbdtaz2zRdsZLBsjDfq5/xGHY3MrLUxTiu767rWuKfyIW2NTYaMe1pjkFYdtfJD3/KVMu4Z0ormWHdA1tfB455a32tlVW2Oexbff5vHqVJHKyN9NuqujnFaCX2HZq8f8EUz6DupOQWq1ZDKsW6VbGVsHZLVvnXOHyEjff283GjH0Iz01THVwDHtgA/vlpxurmqNEQZwJp0kSZIkSZI0MYN0kiRJkiRJ0sQM0kmSJEmSJEkTM0gnSZIkSZIkTczEERNbv0sQEyaOGKWOoQklBiwu3l60ecDrDUko0Sy7+GQQzeMxZPvA1SprxVt5I5oLKy96Y1s2fvmsLA4azWNa376+sqhsq47m4sKV7a0kE7VFWMv2xSeOGKOOli1fknSYoeunthY0HlK21oebSSYGLKzczKkyYCHndh2Lz2bQLlvfXM2IMLSOIW/k0M//kMXgB/yLMFufr8bvXlssubVwcWut6kFJHwYscjyoLAMXRV78aW/B8oPqWGZLOu4ZIfnEkPHQGGOnoeONQQmzGt+NtbqbdYwx9hzjnNU6Hw74Jm2d82ksEL6htih645i2vr+qi8G3zofNZBCLr6M9HqqNWapFB41lVsO4Z8j4Zqw6qv1pQIIuaIxPBo57qtsHjJGAeuKIgX/fDP5uHFJ4QIKA9svVsi00iraSMg5JdjVG0ocB59r2GGTLk1K0TD3uWd86QQ7gTDpJkiRJkiRpYgbpJEmSJEmSpIkZpJMkSZIkSZImZpBOkiRJkiRJmphBOkmSJEmSJGliZned2NU7ATtfwycPSBwyKJNWo+7BdQype4wMsYOzrW2atmVwO2pZzoZkjmuVH+NYD8xyFpUXbWZ3q2VbWqDuWhqwZsbcZkajxWcX29DM7lpr2uLLNl+zlW2tXkWj3oFphJY73eIIJ4BmhrIR2tHMvDcgQ+ygzGUD66hmLhtStlV+pCxno5zga32y8SFtvtzaWr2Nsq1zSDXLWSOLWDPDXq2ORjsG9L0xMo4165niFLJUp6ERuuMk454xvufHyGrfrKOWcXRYHdXsro2i7d9xhI5TTXc7LBN0Vr8gBmafbI5lKuWHjntqWVVb47JaWmsaY4sh45tGO1qaVYzxnq/wcc+grPYDXnPI+KbZjoHf/dXXHDquGDLuGTCmao5vxshq3zLGuX3I35hDzuHUxyet/jgou+sSZrUfatnHPVvQ7qtH+J2dSSdJkiRJkiRNzCCdJEmSJEmSNDGDdJIkSZIkSdLEDNJJkiRJkiRJEzNIJ0mSJEmSJE3M7K4TW78DxI6bKTRWVtWKQYmLljTL2bA0KKNkShujjtrGoZnSlipbVfMFG69Xy4jYqmNgRrnqMWllZ2oej0rGq4HvVzWj7xj9euB72MygtqXG6kujpXKeV+3Q5g1pxxhZM4dkbB0jS+rQrFlDMrY1s3cNy3I6oIrqR6NZtnWerLVj6L8Ta9knW1kfh7wHzfelVcemv/zgj+hI2WBrlixR4tB6l+h8OMppbOh3TMOQcc+g77Wh38VD6mhUsVTtGPp+1Q5fNRN3qzD1rPbt36WVSnPACXH9wHFPZbAwdCxZ3TpGv3bcs7hqxxj3jPG9MTQj/Rhjp5rG2KT598ag8Vdr+xL+IV37u2LNwGzytWNSy3QP7Yy+lXFPNDrfkGzBw8csA8Y9Y42HKpY0EfQi6l5/9Za/jDPpJEmSJEmSpIkZpJMkSZIkSZImZpBOkiRJkiRJmphBOkmSJEmSJGliJo6Y2IYdk/U7XsPVDZdwHcwx6p4iKUW1HWO85tBFkat1tBYG3sJ6ob7YZ7PogMVZB7QZIForA9cWRR3lmA4oC0u4cnHr9QYs5Dy47jEqGWApF2FtvuZSrSi/dOWbix8PqXcJF3Je0sV0a9YMXbh48YsON3+VAQtNLz5NzfA6hqwYPsr7MslndInqXe7z20ivudzjnsGHf0jCrIZByb9Wxldue9yzYcBntDV+GCEJwzj9ZmD7ttTQhG+jvOYYlQywpIvPL/MvM8W4Z9CYZcD2gcdujMQFw/Lotcb+jUpqSSLGaPMIp+XhiU+2oXHPFtS94eotb5gz6SRJkiRJkqSJGaSTJEmSJEmSJmaQTpIkSZIkSZqYQTpJkiRJkiRpYgbpJEmSJEmSpImZ3XViG7aD2H7qVizStpRNtmWMVDNjZEobwxjJ1hqNG9zmAcc1ljJD2XJnBlvubLLbmqXK2rRSskGN8OEffMpawt+99us029fMiLjpE4ZnFxvjC2XLqxhS7ygvN0V2123JCh/j1Ewy7lkp47Uh2ScHNKTZtqGfryHjnqUcKyxZdlfHN1tkuc/XKyWT5liZWQe95gh1tKqujXtahQekqh98jhzwhOawbJmP9eCX24rGOBtGiLA5k06SJEmSJEmamEE6SZIkSZIkaWIG6SRJkiRJkqSJGaSTJEmSJEmSJmaQTpIkSZIkSZrYqs7uGhHbA88HHgPcFPgl8H7g5Zl50XK0IbdLcrutKF3JUtlKk0QtWWbWpTRGAtuhGRtXyIFaIc3QCjVKZqultNLbVzFO9snG5iVMUWa2VS25Ff59tE19Xw7JBDt0fNN8zeU9gNvU+6XROb5ZPoM/i7UMsUuYTbrFcc94xojtrNqZdBGxE3Ai8BrgRsA3gOsCzwFOjog9JmucJEmSJEmSVpVVG6QDXg/cHTgduFlmHkqZTXcacABwzIRtkyRJkiRJ0iqyKoN0EXEz4Cndw0dm5vkA3e2R3fbfj4jrTdE+SZIkSZIkrS6rMkgHHEFZj+/UzPxqf0f3+LvAjsCjJmibJEmSJEmSVpnVmjjit7vbExr7TwQO6sq9cSkbkmvKj7RNWOELhq74hXMlLdqKXyh9pbdP0uKt8PGD4xtp27HixzctW2u7RzZGbGe1hodu3d1+u7H/B93tDZehLZIkSZIkSVrlVt1Muoi4FiWLK8D/NIpd2N0OCtJFxMWLKTakTkmSJEmSJG11zoxYeL5zZu7Wf7waZ9L1D8D5jTKzIN1eEbEaj5EkSZIkSZKW0aqbSQdc1bu/tlFmtn0NA1ahmB8BrYmIGwC/mD1ef+mli61eanN+pqTl4tpHkpaL4xtJy8XxjUZQie/cLDN/OaSO1RikOx9YB+wA7NEoM7sc9pLMHPvjer3+g3Nee9TI1UuSJEmSJGli1wMGBelW3aWcXdDtnO7h9RvFZkG6s5a+RZIkSZIkSVrtVl2QrjML0t2msf9W3e33l6EtkiRJkiRJWuVW4+WuAF8C7gocBry4sv++3e1JS/DapwMH9h5/tbu96RK8ljSmM7tb+6pWMvupthb2VW0N7KfaGthPtbWwr277go2XODt9cAXjL7m28kXEIcAp3cM7Z+Z/9fbdFfgiZd26m2bmOZUqxmzLxbC4pBPSlOyr2hrYT7W1sK9qa2A/1dbAfqqthX1Vi7EqL3fNzK8CH+0eHhcRewB0t8d229+61AE6SZIkSZIkCVbv5a4ATwK+Atwc+FFEfA84CLg25XLYF03YNkmSJEmSJK0iq3ImHUBm/gI4BHgLcBlwe0pCib8A7p2ZV0zYPEmSJEmSJK0iq3JNupXE69K1tbCvamtgP9XWwr6qrYH9VFsD+6m2FvZVLcaqnUknSZIkSZIkrRQG6SRJkiRJkqSJGaSTJEmSJEmSJuaadJIkSZIkSdLEnEknSZIkSZIkTcwgnSRJkiRJkjQxg3SSJEmSJEnSxAzSSZIkSZIkSRMzSCdJkiRJkiRNzCCdJEmSJEmSNDGDdJIkSZIkSdLEDNJJkiRJkiRJEzNIJ0mSJEmSJE3MIJ0kSZIkSZI0MYN0kiRJkiRJ0sQM0k0oIraPiD+NiO9FxOUR8ZOIODoidp+6bVo9IuKOEfGmiPjPiPhl93NiRDxhgefYdzW5iLh11w+/EBHnR8TLK2Xsq1p2EbFjRLwiIn7Q9btzI+I9EXGjBZ5jX9XkImLXiHhhRHwkIs6MiAsi4jPdtrULPM/+q62CfVXLLSL2iog3R8Q3I+LSru+9JyL2X+A59tNVLDJz6jasShGxE/Ap4O7Ar4DvAQcBuwE/BO6amedP10KtBhFxEPCd7uHl3f0bAzfstv1zZh4x7zn2XU0qIq4D/APwB92mdcCPgddn5j/2ytlXtewi4ibAvwMHA5cBPwBuBewAnA08KDO/Me859lUtqYi4PvAS4OnAdsB+mfmTeWUOAT4C7A1cROm7NwL27Yp8HvjdzLx43vPsvxrNYvpqV25Hytg1GlVdkpm7zXuOfVWjGNBP7wJ8Erg28D/Aj4ADgOtS+u+DM/PT855jP13lnEk3nddTPninAzfLzEOBmwKnUT64x0zYNq0Smfld4HPAQ4DrZOadMnMf4BldkT+MiIfMe5p9V5OJiNsBp1ICdN8EngrsnZkH9gN0HfuqllVEBHA8JUD3fuD6mXkwcD3gQ5SAxwcjYrt5T7WvaklExE4R8ULgDOCPKX9MtvwE2BH4I2CvzDw0M38DuAflj8t7Ai+oPM/+qy02sK8C7EcJ0J0LfKby87nKc+yr2iJD+mlE7Ar8P0qA7jWU8epdgBsAbwN2Bt7RBeX67KernDPpJhARN6N8yLYD7piZX+3tOwQ4BbgS2CczL5imlVrtIuIzwH2Bd2XmE7pt9l1Npvuv5amU2Z7HA4/NzCsbZe2rWnYR8Wjgnyh/NN603z+7QfgZlEDdozPzvd12+6qWTO8fG+cBLwfe1O1qzfq4Xq2fRcQzgDcDFwO7Z+aGbrv9V6O4Bn31MOAE4LWZ+eJF1G9f1RYb0k8j4n7Apymzkn8re4GXbkxwDmVG3aGZ+ZVuu/1UzqSbyBGUD96p/Q8eQPf4u5T/ZD5qgrZJM6d2tzfpbbPvakpvoAToPg4c0QrQdeyrmsLh3e275/fPzLwC+Ofu4RN7u+yrWkrfolyStX9mvnlzhRf4o+/j3e1ulBlMM/ZfjWVQXwVu1t2evsj67asaw5B+etvu9jv9AB38ekxwRvdw394u+6kM0k3kt7vbExr7T5xXTprCb3a3/f8K2Xc1ie4/i4+krD/3nPmDnQr7qqZwq+72+439X+xub97bZl/VksnMDZn5yvnryF0D/efv3rtv/9UorkFfnQXpzliw1Bz7qrbYwH46WzfuoG45jF/rlr24afewP2awn8og3URu3d1+u7H/B93tDRv7pSUVEbcG7tc9PLa3y76rqTwKWEsZtPw0Ip7UZcb614h4cUTcYl55+6qmsGd3e0Vj/8+7230iYofuvn1VW4ODutukXIo1Y//VVGaZMRc7k86+quX2Ico/OA4EXj4L1EXEGuCVwPWBEzPzO73n2E9lkG65RcS1KNeeQ1mEt+bC7tYPn5ZVROwTEU+lrJ+wK/CGzPx8t8++qyn93+7225TFoN8OHAn8PvBXwKkRMVs70b6qqfywu71lY/+G7nYNsId9VVuR+3S3P8zMS8FzrSY3m0n33Yi4IiLOiIjjI+K+8wvaVzWFzLyQ8k/mdcBfAidHxDOB/wT+FPg65fJWwH6qOQbpll8/FXgrdfLsw7dXF2mXllREPDsiEjgbeCslQPeozHx+r5h9V1O6fXf7Ako2t7sB16Jk0fwEJUPW2yPitthXNZ3Z5axPnp+trbu0pb+4+YXYV7UViIjrAc/uHvazCtp/NaX1wJcolwpeQAna/QHwmYh4X0Ss7ZW1r2oSmfkR4F6UQN0dgTcCdwV+AfxBZp7XK24/FWCQbgpX9e6vbZSZbV9DuaxAWmo/o6xx8C3gUkrw49UR8ZheGfuuJhERO1P6JJSsmYdl5pcy89LM/DpwGCXRyRpKinv7qqbyV8AllEWgPxURh0TELl3w+N+B3+3KXdwtGm1f1dbglcB1gLOAd/a22381mcy8XWberfvZh3LenfXPIyiL+8/YVzWJiLg/8EFKMogPAEdRlr7YC/hORDyyV9x+KsAg3RTOp0TSAfZolJlNc71kEYujS1ssM9+fmffLzNsCNwCeCewDHBsRz+2K2Xc1ld1791+QmRv9d7Hra0d3Dw/FvqqJZObPKRleLwDuDpxC+cfHNyhJJV7VFf1ed2tf1YoWEX8IPIPyx+CT52Uttv9qxcjMszPzicB7uk3P6s1otq9q2UXEnSkBuusAD8nMh2XmCymzPt8G7AS8JyLu3j3FfirAIN2y6z5M53QPr98oNvvwnbX0LZI2lpmXdSnFn9FtelVEXM++qwmdR7msBeCnjTKzoMfulECzfVWTyMxPALcDng/8M/APwBOBQ4CbdMU+05X1vKoVKyJuBfxj9/DVXd/+NfuvVqg3dLe7AweAfVWTeQ2wI/C33WWvAGTm5cDTgc9S4jF/0W23nwowSDeV2YfvNo39t+puv9/YLy2Hf6EERnYB7tBts+9q2WXm1cwF527UKHZ17/7l2Fc1ocz8WWa+ITOPyMynZuY7gSspl2YDvL9X3L6qFSci9gU+Rlmj9mPAyxpF7b9aac7o3d+ld9++qmUTEdtTZtQDfHT+/i4gd3z38JDeLvupDNJN5Evd7WGN/bOsRCctQ1uklsuYC3xs193adzWVWcr5Ozb2/2Z3e1Zm/gr7qlaevwX2BD7RraU4Y1/VihIRewCfpKzx9WXg8Mxc3yhu/9VKc4vudj1lreUZ+6qW027M/f3UWl/uV93tJb1t9lMZpJvILGp+m4i4U39HRNwVuDnlevQPLHfDpJ4HUaZoA3y1u7Xvaipv7G43yZrZeXx3+7Hu1r6qFSMingI8mbI+3fPn7bavasWIiF2BE4ADKcsIPCgzL1vgKfZfrTTP7G5Pmdd37ataNpl5ASUxH8BDGsVmAbdTetvspzJIN4XM/Cpz016P6/5jOfvP5bHd9rdm5jm150tjiYgHR8SREbHbvO2HURY0BXhzZv4P2Hc1ncz8OHAysDfw4Yi4zmxfRLyQ8h/HS+kuybKvampRHBwR/0I5n14OPCYzv9cvZ1/VShERO1Auxb4zcDpwv+4PzSb7r6YQEb8fEQ+IiLW9bbtGxOuAx1GCGE/tP8e+qgm8pbt9ZvfPOuDX44OnAU+gXLU0S35mPxUAYVKQaUTEXsBXgP2A/6X8t/Ig4NqUaa73zcwrJmugVoWI+A/KegkJ/ISSEvwWzC1KegLw+/2+aN/VVLo1kj4AHEwZ1HyLcjnWnpS+eGRmfrRX3r6qZdcNvJ9ISRKxZ7f5m8Aj5wfoes+xr2pZRMRs4L9fZv5k3r5/Ah7dPfw28ItGNZ/KzKN6z7P/anSb6avfAG5LWevzdEpQ7kBgZ0p27Wdk5vHMY1/V2DbTT9dQEkgd3m06j7Jm4v6UJGdXAc/tEvb1n2c/XeUM0k2oi4i/kjIF9rrAfwPvBY7KzHXtZ0rjiIgDgCdR1vnan5Lu+0fAacA7u9lLtefZdzWJiNiZMlvu4cCNKYOXzwF/kZmnV8rbV7WsIuIo4I+BsylZXE+grEO3YH+zr2o5bOYPyjMpfxRuzrsz8/H9DfZfjW0zffVOlPHr7Sl9djvgO8B/Aa/NzPMWqNe+qtEs1E97ZQ6jzOz8LWAfyt9aX6dkzd5k7No9x366ihmkkyRJkiRJkibmmnSSJEmSJEnSxAzSSZIkSZIkSRMzSCdJkiRJkiRNzCCdJEmSJEmSNDGDdJIkSZIkSdLEDNJJkiRJkiRJEzNIJ0mSJEmSJE3MIJ0kSZIkSZI0MYN0kiRJkiRJ0sQM0kmSJEmSJEkTM0gnSZIkSZIkTcwgnSRJkiRJkjQxg3SSJEmSJEnSxAzSSZIkSZIkSRMzSCdJkiRJkiRNzCCdJEmSJEmSNLHtpm6AJEmStFwi4oHAId3Dd2fmT6ZszxQiYlfgBd3D/87MY6dsjyRJKiIzp26DJEnSVi8izgT2A67IzJ0nbo4aIuItwNO6h/fKzM9P2Z4pRMRewM+7h5/LzHtP2R5JklR4uaskSVoyEXHtiHhERPxjRHwjIs6OiHURcUlE/CwiPh8Rr4uI+1aee8+IyO7nLSO05eCIeG1EfDUifh4RV0bET7s2/ElE/MaWvsZq0r0/53bv695Tt2el83hJkqTN8XJXSZI0uu5yuucCzwOuWymyPbArcGPgHsDdgTsvUVtuDLwBOLyye9/u5x7AKyPidcArMnPdUrRlG/MCYO/u55GUY6w2j5ckSVqQQTpJkjSqiNgf+CBwq27Tr4D3AycCZwIXArtTghUHAA8G/m2J2nIr4JPAPt2mU4F3At/q2rEv8ADgMcC1gT8H7hgRD8vMS5eiTduQtY37vxYRBwAPBU7JzJOWpVUr12aPlyRJWt0M0kmSpNF0AbqTgT26TccAL8/MCxd42muXqC17AycBewLrgRcCR+fGC/J+G/hYRBwF/CtwJ+C3gX8CHr4U7dqGvJ6SgOEc4L3zd0bEEcD7uoePW75mrVgLHi9JkiSDdJIkaRTdJa4fogToEnhcZv7ThE06lhKgA3huZr6xVTAzfxoRvwN8BbgF8LCIODIzDaY0ZOaJwA0WKLLDcrVla7CI4yVJklY5E0dIkqSxPBs4qLv/11MG6CLiHpQZcQCnLhSgm8nMiyjrhs38ZUTEEjRPkiRJ2oRBOkmStMW6WXSzANd5wCsnbA6UpBUzRw143gnAD7r7BwD3G61FnYjYOyJeERGnRMQFEXF5RJwZEe+JiP+zwPP26GW7fVa3bf+I+LuI+F5EXBoRZ0TEsRFxk0W25bYR8Y6I+EnXjh9HxL9GxKHd/jd2r3dF5blv6bXnnr3tn42IBN7dK/7uXtmMiMf2yr+0t/0RC7T1sb1yL1qg3L4RcXREnNb9Tj+PiBMj4sjFHJNePTeJiL+OiO9ExMURcWGXmfWlEbHH5mvYpL7q8er2Hdzb93vdtttGxLsi4ocRcVlEfD8i3hQR11/Ea0VE/GFEfDIizuuef1pEHBMR+w1o86yej0bJzHxF11fe3808nV/+gb3f4+SIqP6tEREP75X70GLbI0nSts7LXSVJ0hjux1wW13dk5mVTNSQitgf+b/fwbMpac4uSmRkRxwBv6TbdF/j0iG17MvC3wC7zdu3X/TwyIo4G/jwzr5xXpp/I4jpdQOud8+rav/s5PCIen5nHL9CWlwAvZeN/2t60+3l4RDynt/2ihX6vlSIi/hB4Bxsfk52AvYD7RMRDgfMXUc8TgDey6fu0O3Bb4KkR8aDM/PoY7WbT9/Y5wOvYOMHELbufR0XEgzPzPxptvw4lUct95u06oPt5HPCUzTWoCwZ+AJgfOP6N7uehEfE24OmzdR4z84SIOBZ4LGV9x2cCfzev3mtRPgNQ3ounbq4tkiStFs6kkyRJY+jPqvnUZK0oDgZ27e6fnJlXD3z+yb37dxunSRARjwf+gRL4uRT4C+BQ4PbAE4HTKWOz5wMfq1Sxrnf/fsB7gP+lBEIOAR5GCc4A7Awc25o1FRHPBl7evd4lwMsoAcm7UTLcXgwczdzvv27TWpoeBdwc6M92e1G3bfbz/srztkhE3JuSqGIXypqI/wgcTukPT6DMkHwY8MjN1POI7rm7AKcAT6Yc3wdTAncJ3BA4ISKuPVLz+8f3sZRj/2NKQO0OwBHAZ7v9uwP/FhG7Ner6AHMBuh9S+tNdKJd/v4nye71jocZ0ge4TKQG6dcDfAIcBdwaeAXy3K/pU4E/mPf05lOQYAK+OiH3n7X8FcOPu/h9l5s8XaoskSauJM+kkSdIYDujd/85krShu3rv/39fg+T/u3b/hljWliIi9KNk9AS4ADsvMfjDwGxHxEeAzwG0os74elZnH9cr0s9LeA/gpcJfMPLfb9rWI+CAlM+2RwI7A84BnzWvL3pRACcC5wP0z81u9Il+KiA8AX6YEEAE2LPZ3zcxzutf5RW/zLzLzR4utY6iIWEvJJBzA1cBjMvP/9Yp8PSKOAz4B3HuBenaizGAL4Hjg0Zl5Vbf7a8CHI+IrwHHA3pTg6tEj/Ar99/a+wNeBe2bmJd22UyPi/b3270kJHr5hXvsPZ+73+wLwe/MyK386Ik4APrqZ9jyV0g+v6NrxX719/9Udy5MoAdDnR8Qxs5mfmXlRRDwN+DBwLUpg8MFd+27LXH/8t4VmekqStBo5k06SJI1hlrUyKUGoqoj4o3lrk/V/3jVSW/prdp059MmZ+b/ALLCx2fW/Fun1zF0O/Nx5AbrZ655Hmek1C4i9LiJ27u3PeU95bi9A1y/zqt6m+1fa8lpgNgPsz+YF6Gb1fJ+5QB5sHERaiZ4E3Lq7f+y8AB0AmbmOErxcaGblsymzvC6gXMZ51fwCmfk+SsAO4I+3pNH9auc9fkovQDd73XWU925mo/c2InZgbv3FWaCyH6Cb1fMJ4G2thnQz9F7SPXztvADdrI5fUWZfQvnsP2Le/o9QZnoC/F5EPLC7/ybKJby/pMzIkyRJPQbpJEnSGLbv3Z86oNNfR2yThAeLtEN3+6stbAtdoO3h3cNzgX9ule1mm80udd0buHuj6JVAdcH9zDwNOKt7uH9/8f4ukHN49/As4L0LNH2TQNcK1r+E9W9ahbpZfp9t7Qce3d2eUAtw9Xyxu71pRFx3gXKL1f/M/Cgzv9Yo99le2d+ct+9uwCxhyL9m5kIB6uMW2Pc7lJl6MBdoq/li7/7Blf3PBmaXsh4TEU9ibn27p3dBaUmS1GOQTpIkjWH2B3cwN2Os5r1svDbZvZagLf2ZfDcb+uSI2JO5Ne1+OkJ77kBJXgDwn92MqIWc2Lt/SKPMzzJz/QJ1zNYEW8PG78ftmQtifm2h9fq6WXqbzCRbaSIiKGuuAZyfmT/czFOq72kXzJxdKv3oBWZ8JhtfQnzTLfoFNtUMrnUz+2aftfmzPA/t3f/SZl5joX59i979Hy1wDPqfs02OQRfknCWF2J+yHiPAcZk5+pqEkiRtCwzSSZKkMZzTu/9brUKZ+b+Z+aPZD9dszbjN6c/QuXmzVFs/sPeTLWwLlBlxM4u5/LZfZs9Gmc0lw+hn1+2vQdxvS3/tvZaFAoErxe7MzXxczPFtBR73o6zjN1Rcg+csZLHv7fy1pffq3d/ccVgo+HqLBfa1VI9BZn6YsrbfrMzFjHeJsCRJ2xwTR0iSpDF8nrl1qe4L/OeEbTmld/92ERGV9dwW0r907/QR2tMfby0moLND7/6VI7x+q+5LR657DNck4LVT7/6QLLTz9YNcxwKvXOTzxphtOYalOA4HUzIIb84ltY1dIo7b9DbtRrmE+8PXuHWSJG3DDNJJkqQxfLJ3/wkR8eraovvLITN/GhFnUGbR/SbwIOAji3lud8njc2ZVsfCaXIvVn9k3fx2xmn6m3J83S10z/UtB91uoYJcxdfuFyoykH1C6zgLlWleA9C+73G8Rr7dTY3u/nhssZTbaJXJ+7/5+mynbOgaw8XHYbguPwyuAAynB5kuAPYC/j4jPdcknJElSj5e7SpKkLZaZP2YuELYvG6/ZNYV39O7/+YDnPYy5QNrHRwrUfJ25y0bvGBGb+ydpP1nEl0d4/b4fMpc9trbYf99NKJk4r6n+7MWFZshd1Lu/d6sQjbXfMvNK4LTu4Q0jYp/NtOuAxvYfMzcj7FbdWndbk36W3tZahjOtYwDwzd79WzdLbUZE3Bl4Xvfw9cALu/s3YoHkHpIkrWYG6SRJ0lhextx6Wq+JiPtN2Ja/Z25G0J27zJILiog9gL/qHm4AjhqjIZl5EfCp7uGNgccv0Ia7ULJrAvyMjS/dHaMtlwNf7R4eGBGHLVD8iC18uf/p3V9ohtxpvfsL9ZnfWWDf7PLJAJ7bKhQRtwTuWtvXzfz8aPdwXxZ4n1aoE5lbr+7RXQKUlicusO/DzAWVX9RdsjpI95x3UYK8Z1M+V+8CvtIVeUpE3HNovZIkbesM0kmSpFFk5teBF3QPtwc+FhHPi4gdFnjaUrXlV8BjmZvN9ZaIeHqrfETcjHLJ7mwW3fMy83MjNulZwBXd/VdHxJ0qbdiPjWcAvjQzN8wvN4IX9+4f0/3u89tye+ZmPl1TZ/fu33GBcl9hbn28/xMRd6u050+A/jGbP8vttcwFZZ8REfev1HEt4Lj5m+c9flGvLUdFxOZmG64YmXkeczPUdgbeFRG7zC8XEY8E/qC/aV4936MEuaF8Ho6JiKEJNWaXuQK8MDMv7daFfAYlABjA2yNi54H1SpK0TXNNOkmSNJrMPKYLDLyKEqh7PfDsiDge+A9KFtgrgOtTZis9fJFVHxARj9pMmbMy8/O9tnw0Ip4FHEMZ8/x9RDwWeCfl0sBLgd+grFn3aGBXSlDvbzLzmEW2a1Ey84yIeAVlRtGewOcj4mjgs8DllEtcn8tcNtfjgHeP2YZeW06KiA8BD6as2/fliPg74AuU9eF+G3g+cCGwC+XYXZNEBD8AzgVuCDwiIk6hzPa6FvDdzLy4a89lEfE25i6N/FhEvBT4InBd4HGUWX2fB2azrza6DDczL4iIlwBv6tr8kYh4M3AScBbl0t4/o6zV9h7K+w3z/mHdrWf4MsosyusBJ0fE24FPU5KIBLAP5f16MHCHzNySJA1jOwp4AqVfPxD4QkS8AziZciyPoASvv0X53W9N/Z/2f0n5bO4DPAW4V0QcA3ybssbidShZnB8CfCEzfz3rtJsNOnsvT8zM9832ZeapEfEm4NmUAODLKIFRSZIExLBkZ5IkSZsXEfehBAzuMOBp78jMX1+W2l0O97kBz/9gZj600pYHUII3m8wYm+fLwHMz8+QBr9l/nTMpQaArMrM6QygiHgf8Le3LPxN4A/DiWuKNiJgN3H6QmQfO398r91ngXt3DvTPzF/P270Q5Jq3LHs+lBC+/BOwInJKZG83+i4i3AE/rHt6rHyDtlXkm8MZK/bfMzNN65XamzGS8e6UslMuFj2AuOcKrMvMvK6/3NMrxrc38SuAllGDxrK0PzMyPVep5FPB3lEDdQu6emV/YTJlZnc3jFRE3Af67e/iJzHzAAvUs2M8i4sbA8TQu6wW+D/wu8Hbg3lTe266evYG3UoKRCzkxM+/XPWcn4FTglpTA7q0zs5+shIjYrWvDjSiz6u6cmV/bzGtIkrQqeLmrJEkaXWaelJkHA/cBXkeZyXMOZdbYVZRLE79PWf/qZV25p1Ur2/K2fJxy6d1jgH+nBEOumFdsA3AmcLuI+K2lShqQme+mzF56I2Wm2WXALykzvt4AHJSZL1jqzLiZeUUXED2CEgi9kDKz8LuUWZC3oyRSmAW7zt+0lkW9zpuAI5m7pPUiypp4F88rdzklqPhMyjp8F1MChSd1bTwsMy9gLrHD7o3XeytwKCVIdS4lq+iZwPuAQzPzVb06AK7dqOc44CBK3/1m1+6ru9tvU2bjHUlJCrKiZOZZlBmHL6IEzC6lvL+nUGa43alLiDI7Dq1j8PPMfAhwOPAhyuXL6yif4bMp780rmMuGTPf4lt39v54foOvqvZgykw7KjMh3LCKZiiRJq4Iz6SRJ0qoUEbcA7k9Zn+tQNl6b6yLga5SA0ku7DKKrSkQcSplJB/CazPyzKdsjSZK0rTNIJ0mSVr3uEsHDgYdRAnazNc++l5kHTdawCUXE6yhr0wHcPzM/OWV7JEmStnUG6SRJknoi4nrAAygL75+SmUdP3KTRdeuCXZKNgWC3puCnKMHK04ADW2UlSZI0DoN0kiRJq0yXYfOBwHuB71DW6VtHSa7xUOCRlMt/1wH3ycwvTtNSSZKk1cMgnSRJ0ioTEacDN99MsXOAx2Xmp5ehSZIkSaueQTpJkqRVJiJ2piTMeDiwH7AvZebcj4AzgM8A78vMS6dqoyRJ0mpjkE6SJEmSJEma2JqpGyBJkiRJkiStdgbpJEmSJEmSpIkZpJMkSZIkSZIm9v8B/BHSlNavoRUAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "image/png": { "height": 340, "width": 628 }, "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "grid = coeffs.expand(grid='GLQ')\n", "fig, ax = grid.plot(show=False) # show=False is used to avoid a warning when plotting in inline mode" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Next, we would like to calculate the variance of this single spherical harmonic. Since each spherical harmonic has a zero mean, the variance is equal to the integral of the function squared (i.e., its norm $N_{lm}$) divided by the surface area of the sphere ($4\\pi$): \n", "\n", "$$N_{lm} = \\int_\\Omega Y^2_{lm}(\\mathbf{\\theta, \\phi})~d\\Omega$$\n", "$$Var(Y_{lm}) = \\frac{N_{lm}}{4 \\pi}$$\n", "\n", "When the spherical harmonics are $4\\pi$ normalized, $N_{lm}$ is equal to $4\\pi$ for all values of `l` and `m`. Thus, by definition, the variance of each harmonic is 1 for $4\\pi$-nomalized harmonics." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can verify the mathematical value of $N_{lm}$ by doing the integration manually. For this, we will perform a Gauss-Legendre Quadrature, making use of the latitudinal weighting function that is stored in the `SHGrid` class instance. Note that it is necessary to ignore the redundant longitudinal band at 360 E." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "N = 12.56637061435917\n", "Variance of Ylm = 0.9999999999999999\n" ] } ], "source": [ "N = ((grid.data[:, :grid.nlon-grid.extend]**2) * grid.weights[np.newaxis,:].T).sum() * (2. * np.pi / (grid.nlon-grid.extend))\n", "\n", "print('N = ', N)\n", "print('Variance of Ylm = ', N / (4. * np.pi))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Alternatively, we could have done the integration with a 'DH' grid instead. In this case, it is necessary to ignore the redundant longitudinal band at 360 E and the latitudinal band at 90 S." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "N = 12.566370614359172\n", "Variance of Ylm = 1.0\n" ] } ], "source": [ "grid_dh = coeffs.expand(grid='DH')\n", "weights = pysh.utils.DHaj(grid_dh.n)\n", "\n", "N = ((grid_dh.data[:grid_dh.nlat-grid_dh.extend, :grid_dh.nlon-grid_dh.extend]**2) * weights[np.newaxis,:].T).sum() * 2. * np.sqrt(2.) * np.pi / (grid_dh.nlon-grid_dh.extend)\n", "\n", "print('N = ', N)\n", "print('Variance of Ylm = ', N / (4. * np.pi))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Parseval's theorem\n", "\n", "We have seen in the previous section that a single $4\\pi$-normalized spherical harmonic has unit variance. In spectral analysis, the word *power* is often used to mean the value of the function squared divided by the area it spans, and if the function has zero mean, power is equivalent to variance. Since the spherical harmonics are orthogonal functions on the sphere, there exists a simple relationship between the power of the function and its spherical harmonic coefficients:\n", "\n", "$$\\frac{1}{4 \\pi} \\int_{\\Omega} f^2(\\mathbf{\\theta, \\phi})~d\\Omega = \\sum_{lm} C_{lm}^2 \\frac{N_{lm}}{4 \\pi}$$\n", "\n", "This is Parseval's theorem for data on the sphere. For $4\\pi$ normalized harmonics, the last fraction on the right hand side is unity, and the total variance (power) of the function is the sum of the coefficients squared. Knowning this, we can confirm the result of the previous section by showing that the total power of the `l=5`, `m=2` harmonic is unity:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Total power is 1.0\n" ] } ], "source": [ "power = coeffs.spectrum()\n", "print('Total power is ', power.sum())" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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.9.1" } }, "nbformat": 4, "nbformat_minor": 1 }