{ "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": "iVBORw0KGgoAAAANSUhEUgAABOkAAAKpCAYAAADpKXgkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAACMQAAAjEAFnL93WAAD+YklEQVR4nOzdeZxkZXn+/+tU9TIrO8riwogRBVEEUdSoGCSGVUGUmCAS0ShqjAGjURMVNV81wSj8VDQgqCgiiJCImAhCAgRRRAGRRZZBQZCBQWC2Xqrq/P5oZ2So++o59/Sprumez/v14pX49NPPeerUWe46U/1cRVmWpQAAAAAAAAD0TaPfEwAAAAAAAAA2djykAwAAAAAAAPqMh3QAAAAAAABAn/GQrs8WLlyohQsX9nsawDpxrGIm4DjFTMGxipmA4xQzAccpZgqOVVTBQzoAAAAAAACgz3hIBwAAAAAAAPQZD+kAAAAAAACAPuMhHQAAAAAAANBnPKQDAAAAAAAA+oyHdAAAAAAAAECf8ZAOAAAAAAAA6DMe0gEAAAAAAAB9NtDvCWxsWq2Wbr311jX/u9PpqCgK3X///X2cFbBuZVlKEscqNmgcp5gpOFYxE3CcYibgOMVMwbE6+3U6HT344INr/vcf/dEfaWAg99itKFcfKZgWN910k3beeed+TwMAAAAAAAA9cuONN+oZz3hG6nf4c1cAAAAAAACgz3hIBwAAAAAAAPQZa9JNsy222GKt/73d+/5ezfnz+zQbzBr80TqA6VL0ewIANhrUNwCmC/UNatBesUL3fPxf1/zvxz7/qYKHdNOs0Vj7y4vN+fPVXLCgT7PBrEERC2C6UMQCmC7UNwCmC/UNeuCxz38q/U4P5gEAAAAAAAAggYd0AAAAAAAAQJ/x5659VhZSyaNSFTP1zxlm6rx7ZIN/Hzf0+aG/NvA/cyg38PlNuw18f/B+YTLcL2cW3i/MaBv4/Yj75WPM0P3B+zihjv3A4yEAAAAAAACgz3hIBwAAAAAAAPQZD+kAAAAAAACAPuMhHQAAAAAAANBnBEf020A58d9M0MNplnWMXcsYyZUeg+7pxYWn+e2vZfHj9Guc+gqafZn3TNvebNOrBWh7uLBt6lBPHtRh9+RrmfZFfWvYXnrOdVwsZuCxN1MXmt5g9PJ6HYxdy+Z6Wfe45hl4L07Nua65JS5cPQ2l6NXY1DdTszHciyPmYLfnwIYy74x+1GUbSt1DjTOhhmc7fJMOAAAAAAAA6DMe0gEAAAAAAAB9xkM6AAAAAAAAoM94SAcAAAAAAAD0GQ/pAAAAAAAAgD4j3bXPysGOysHOev5yvXNZe+xpjlDpYSJq4V5LOIYZ2G0vSmxzu878IGytI6QnMWfLzbmWsePm1JHX0+MmOfYUt7dBjR3pR6pSYpu5VNXezSM3bi5VMXyNNhHNJKitc1Lr2N5k6ni/onnXMY902lqP5tFLPY2IdNvs0bj9SI6cgXWPrW8y2yxzEwm7Z2uFzLXMDVFHYmsdCaw11EMzsu7pyzk6zdvbwJMqp7vusR9vUtusoe5J3uuSH8tygjFS9U1i3EnbE33D+iY7j7r6T1Uv654pvJb1frbzKHyTDgAAAAAAAOgzHtIBAAAAAAAAfcZDOgAAAAAAAKDPeEgHAAAAAAAA9BkP6QAAAAAAAIA+I921z4qhtorh9uSdepg4lgr1qiNFzLXXklBmutaSABY3h8lqNtAlHqQM4mPs3rDzqN43lcKW3R+me7RPepoQm3jt6WCgzH4yehZGVNe4Pbrk9DItNJUMZscw52giLdQFtmbGqOV1N5I5bL16zxvJ7SX2k00oqyHdNZUyO91JtW7s7NA9Sy1Ozq1HNU4ytNQMkkmHn2ycGsaO0uR7eP+yybGJhFg3v6LTqwuOazd1T1SvZS+RiferjoTYTH3j9LL+Sm9zqjaGuiedTh6NkahvkttM1T3m60H2WhG+lniDfl9XT7vP7o9U+qwbO9gnqfpmsvaw7wZeO2XGdl17msy87nkUo+t4tlMB36QDAAAAAAAA+oyHdAAAAAAAAECf8ZAOAAAAAAAA6DMe0gEAAAAAAAB9RnBEnw0Mt9Ucbq3fLydWIy2TK5eGSyJmx8gsfuwWETbbDMd28+u4RIlE6IOdR/fYhV1d1Iwdza+OhbQzARFSLeEO4SLMpr9bOLqX86tlIedAT8dwerUIs9O7QzI1dh2LMPt5mPM8+ues7OLCQYCCnYf757PMP6u5c9G+xqkv9psKx3CLR0cnQvafE6PXUkOARWHHqD52kVwM2i5+HIyTPjXqWMm9V6vB13AP7GXdk6lv7FySAQBhneQmYu+jUd1jtucGCe/Frv6Km6PutQQz2Foybs4EZtn6xrSH66cn5xe/xup93Tx6GvqQDfTK2MDrnlxQ1TSP4W4mdYQ+2Pt88Ia549fd5xM1XOHmZ2+wPboHJl6LlK2/zJzDOjXR1/TP1CCSqZOy9Y0z3XXPFI6PdYaCVsA36QAAAAAAAIA+4yEdAAAAAAAA0Gc8pAMAAAAAAAD6jId0AAAAAAAAQJ/xkA4AAAAAAADoM9Jd+2zO8LgG5oyv1++mwpZsQlkmITZu73TiZ702dCzYpkt3zcy7dClibuxEyplNIwra7S61CWDBa6khiCiTkur6uxSxwoTW2OSyKA3OjZ1JiK0lZbZ6X7vNbArudKet1aCOVNX02Im00DCB1Y2dHSPxz1mub9EIzvPs9qKE2KbbnhnbDB3vqBoOykyamevv+tr2YNjkGFF/N4Y9fBvdFxeXZtZwY7sEtUQKW2YMp4Zw51r0su6JapbsGB1bywRtrnYKW+OaxdU3NrE1ut9FCfOTjBEdDVHSvTRJ6mt07GVvMonE3FRqfLK+SdUy2bpnqtuTelY71ZLiuhHUPelU+xrGyCTS97Z2Cn7B9o3f3LDGcWO4cyNbD0V9M+9XojaZGCOREJuoZWzNYtuDtmT9ENZOdoywWY2gdpqsf9y3v3VPa3T9nu08Gt+kAwAAAAAAAPqMh3QAAAAAAABAn/GQDgAAAAAAAOgzHtIBAAAAAAAAfcZDOgAAAAAAAKDPSHfts03mjGpwronnW4dOJpk1OUaUXNY2SWQ2GMz0b0dppsmktHa7e2ybqubmnUhKK9tmX0ftJhlMUcqRZFKz3DuWGcOM4NqDeft0V5Ny5lLRorGT6WKZ+WUS1NLprsEBX0tSmlNHgloNakl3dYdvKjXLjeFiohJjJNptX3M5D/vbJFg3dvdr7JgDxB4eiX+ac8eYHTuTcube82bQP2rTJAllzSBV1bxumyIWjO0SWJvB9tzY9mWbeTQTqa+NZJpZ1D97mrttTlWmvpHiYzJT37j+rm/b1SaJxFZbIwX1zUT/6jWLGztKPyzN9lxSYljjmJrAnf829TXsW7nrJMWuGTt4LY1kfVNHveETZRNj1JEy61J6o/7Z2mQ2pdpnLk81pKqGKanJsX1qfGKMTH0j81rcR96gvvn9b1RqmhjDjJA5Jh138CVqOztGpu5J1EM2Nd7UG1F/X99UrzeargZx8zCHQlQnZeobKXfqTqW+GR8dXe/fXbP9KY8AAAAAAAAAYEp4SAcAAAAAAAD0GQ/pAAAAAAAAgD7jIR0AAAAAAADQZwRH9NkWc1ZoeO4fFkLs2JU9q+sEyyK6BZRde7T4cTSuJLXc4seJ4Ijxdrzap1ucuRX0t9trxws/toPFge1iy2ZR+k7mObdZGDh6Cwq3CL4Trs4ad80sLly0kgERZuxGFPrQMmNkFlA2fRvmPU8toFzHYtA1LMLs1LKYth08MXQm9CG5CHO0WLIPdzABCrWEPgTzcH1Neydozy7CHGmYndexaSZmoMTazKn30S5WXX0B5SgIYqK9epBDw4zhFkWO2t3ixwNujGCR4wGzULJrd4ES0YLGTXMhcosfR+22b+Ii0rCJObFM3ePqkKiWSdc9wTwyNZLk66Go3Y5hghyi2smFTNjwiaDd5WLZMIgwjceM4S44UYBFbm388B6YDXGKQiJsDeJqllQoRfW+rj0Kr5p0jB6FT7h5pAIl6gjG6kd9kxg7G3YV1yzV6xtJqcCsjglsiOqTdO0UPGWwd4dMwTGQq2/s+xj+IDt24kJkapaovTGQq3uiEIZmokaSTM3SjC8sLtRqMOjvaidX37h6KKpDMjWSa8/UN85j657RkVVTH3PKIwAAAAAAAACYEh7SAQAAAAAAAH3GQzoAAAAAAACgz3hIBwAAAAAAAPQZD+kAAAAAAACAPpsV6a4jIyP653/+Z51zzjm68847tfnmm2vffffVxz/+cW2//fbh74yPj+uEE07QGWecocWLF2vrrbfWoYceqg9/+MPabLPNpm3uW85ZoblzTCTTOrhEtDDlzMTSRGlhktQKoghbZns2zczE/YwFyazjUfShfOrreKd7n0WJr5I01jLJsc3uebdaJglWJrooyCmyia8uXigXhFeZDaWySWTVU85sQplJP2sE7XZs8x6EKWfZ+YVpa4kEVklF0N+nu5r2KBUtmUpXQxhRTjZwOEoGS6SZTWyz+0XaJDKb+hrNzaQ+2sTWYB7m7unGju62pXlvo0Q0KXepsAnR5qCsIVg8PkaSKWdRkmvDJLk1TOrYQJCK5tLMoiQyKU5sHTJ93RhDwYXIJZFFfSVpwLYHKWw23dW89qDdJ6L16EaVlKl7xm1fV7MENYFNazX1hmmPklxdX1f3RLVTqxHPb7xh5hEl8pn3vFPEY3SCK1FpL8Bxs7IJ9hXH9vft6kn1to7Jpr5G9Ybta867TO1k655EzWJTcKdes0Rp9+61WBtI3eNrlqivqVkS92JXV2QSW10Cq0sLDWs4l14/YOYXHQuudsrUPa6WNPvUXMrqOZwyqfYumTWoWRoumXUgUfe4esONEdQyg2YMV/cMBhciVzsNmAuAS3eNahxXI7lE2SjJtRd1z6o5I+v9u2u2P+UR+uxXv/qVXvSiF+ljH/uY7rrrLu2yyy568MEHdcYZZ+j5z3++rrvuuq7fGRkZ0Z/8yZ/o/e9/v375y19ql1120UMPPaQTTzxRz3ve87R06dI+vBIAAAAAAABsrGb0Q7qyLHX44Yfrpz/9qV796ldr6dKluuaaa/Tggw/qVa96lX7zm9/oVa96lVqttf8p6thjj9UVV1yhHXfcUb/+9a/1k5/8RIsXL9ZOO+2kW2+9VW9961v79IoAAAAAAACwMZrRD+nOPPNM/ehHP9K2226rr3/965ozZ44kaf78+frGN76hJzzhCbrzzjv1zW9+c83v3HHHHTrllFMkSd/4xje03XbbSZK23HJLff3rX5cknXvuubrpppum+dUAAAAAAABgYzWjH9KdffbZkqSjjjpKw8PDa/1szpw5et3rXidJOvXUU9e0n3nmmWq1Wtp111215557rvU7e+yxh3bZZReVZalvf/vbPZ49AAAAAAAAMGFGP6S74YYbJEnPeMYzwp+/8IUvlCTddttta9ouuugiSdKBBx4Y/s4+++wjSbrqqqtqmycAAAAAAAAwmRmd7nrfffdJkubOnRv+fJtttpEk3XPPPRobG9PQ0JCuv/56SdIzn/nM8Hee/vSnS5Luvffeuqcbz3HoEc0fHp20j0tmddpB3M+4ieSJEtFc/zETZ+iSWUdt+ln3OFFq2WTbHG13t4+24r4u3SZOfc2dEmWQ+holWElSGSV6SnFKUTL4LAymseldJuUs6G9TuhIprpLUGI/6mhSmxNhujEzKWaav3abd1y6pNkg5czFT7rCJjrNeJp+5lDObIhr0tQll1ZPLbF+XDJZIOXOJrdEYLpWuNEmk0bnUGTQps/HQcaCZO35d2GLyOIsHN9e4qN1dy6KUScVJc5kUV0kaDJLLXJrZ0EB8wYnSyIab1ftK0pxm94XPpbgOm4unSy6L+g+ag9Ill0X9Xd+mS47t0UXHJ9Jn0l2nXveMmouCS3d1/aMk15H2YNzX1ENR3eP6Ns29cczFHAbMrVhlsJ9S9Y3ia4UPfHU3n4ptmiz1NdE3XfcE1zKX7pqpexJJsJKre3L1V1TL+ITY6nVPqr6ZpH8topolUd9M9A/aknVP9PUZm55qk+q7d1Q21T7apquR3Hvu5h32Ne3hW5D8fOOuT1GdZA8x91LCdNfq9Y1rz6S4SnEtY+seU7MMB2PYusdccKKxXd9hU99ECbFunDrqHlffOFXqnhVDkz/bqbadGexpT3uaJOnmm28Of95uT7wRnU5HS5cu1fLly/XQQw9JkrbaaqvwdzbffHNJEw/2AAAAAAAAgOkwo79J98IXvlA/+9nP9O///u/6+7//+7XWpRsfH9fHPvaxNf9788031+9+97s1/3tdD+mWLFmiTqejRqP6c8yFCxeus09pv8YAAAAAAACA2WDRokUq/NfFJUnLli1b63/P6G/SfeADH9D8+fN11113ad9999VPfvITrVy5Utddd50OPfRQff/735c08fBszpw5Ghz8w58WNJvx93sf/e27de1MAAAAAAAAoA4z+pt02267rc455xwdccQRuvzyy9dKa33KU56i973vffrnf/5n7bLLLpImvj03PDys0dFRLV26NBxz9bftFixYkH5I99gnoJH7779fj3vc41LjAgAAAAAAYOZYvHixtt5669TvzOiHdJK033776dprr9U555yjq6++Wptssome97zn6dBDD9U73/lOSdK+++67pv+2226rO++8Uw888EA43uqHdE984hN7P3kAAAAAAABAs+AhnTTxQO3YY49dq+2hhx7Sd7/7XUnSoYceuqZ9++2315133mnDJm644QZJf0h57bXHDz2shUND6/W7LqGsHfwVs0s5s+1BhM+oiU/06a5x/1Wd7tcbpZZJ0iqTfray6B7Dps8V8RiFjdOMxPPrBElCPqnSxlgFbVP/U2v38mx7EJDjkk9te5DiOtFePeXMjhGk1TWDcV1fKU4u833N8RQls7rkKDO2oqQ0s16lS5kN9XLNS5vSZ6bSDFZTcH1tYmv3GDahzKSINYL+LlW10TLX1OAS4sKgomuCJBMZ5t6v6olybnM+ea/6FtOHU7SAhks5M+2NoL1pEtEGbGJrkMxqUlznDsYXnCjRbN7AWNg3SnGV4uSyuc14DJfu6trnBBdKl3KWaW+YA6Rp4vRcgtpUZeqbif411D1B+0gnrh9cfWPrnnZ3zeKOG5f6urLVPcZAI95epoKwgc/mPYiSEl19Y+sel2wZDhI3RyOkr3tBu61NEimuE+3V2ibGrp7Y6uoem0gf9E+l15v+rkZy9VD0PhZtM4Y/KM0PapCou21iazSGrW+qp75GCfMTY5jrYTC2S553tVP01rj0VPeWx+9XLjE3TLt1Q5ikWnc8pY4mm+4ajJKue7p3dibFVZKGB6NkVpM8PxBfiOYk6p65ru4JxpjbyNU9Lt11OPgQN+jGSNQ9vahvlg3FrzljRq9JN5mjjz5av/vd73TggQdqt912W9P+whe+UJLWPMB7rIsvvliStM8++/R8jgAAAAAAAIA0Cx/StVotvfe979W3v/1tbbLJJvrUpz611s8PP/xwSdLVV1+95ltzq11xxRW6/fbbNTw8rEMOOWTa5gwAAAAAAICN26x5SPfQQw/pnHPO0Yte9CL9y7/8ixYsWKCzzjpLT3va09bqt8cee+iggw5SWZY6/PDDtWTJEkkTgQ5HHnmkJOmYY47RtttuO+2vAQAAAAAAABunGb8m3Yc+9CGddNJJeuSRR9TpTPxN8Z577qkzzzxTT33qU8PfOeWUU/SCF7xAN954o7bffnvttttuuuWWW7Rs2TK9+MUv1ic+8YnpfAkAAAAAAADYyM34h3TLly9Xo9HQM57xDO2zzz7af//9tc8++2hgwL+0xz/+8br66qv1wQ9+UOeff75uuOEGLVq0SEceeaSOO+44DQ7Gi/b2wrYDD2nTwcnfhnawILIkdcwXIaP+bqHkMbuAcvecRkq3gLJZ5DgIiJCkuZ3uxSZXNeIxBhrDpr174ceBYFFlabKFH7v7u4WS2514Xw9Ei7Oavn5hZTO9jGho97LNAsrxor5x1yiAQcotfuwWUG6OVV/Q2PV1gQ3Nse4X7xdbNgsaB6vp2gWUTbs61edhF0qO2u3qvUmZRb0b5niPFlBOLn4cjd0xfRumPVqcudMy11QTKBEtltw2Kxf7dWaD/i5jwi2KHGVxmIWS/XnuNtqjBbnNW1uYMIhGs3vizaBNkgbNoshRSIQLiJhrFlCOFkueb0If3ALK8wdGu8c1CyhHQRCSNGza5xTd7W4Mv4By934aMn0b5oBygU1Tla17olomqmMm2uOTJgqJyNY9LmhiZXCMrGjF9c2gDcEKAphqCO5wt42OqWWicBwXmFOaGqIWYd1jrm+J4Ih83WPag9PRhl2Z9qi/D6qoXsu4GskGZkXtQR0z2RiKwifcfceMHfftYd3jwiRckFZU47gaydYs3e2NKIhLUmfQ3UeDmmUoHsMd7+FnAvu2JGpGc81yAWJR3aN03ZNsz4jqskQwlhSHRDSDMAkpDsaS4pCIeYNxveHCIKJ2V98saHbXN66/q01cPZSpe6I6RvK1TBgc4e65U6h7HjYBHxkz/iHdpz71qa5156rYcsst9bnPfU6f+9znejArAAAAAAAAoLpZsyYdAAAAAAAAMFPxkA4AAAAAAADoMx7SAQAAAAAAAH3GQzoAAAAAAACgz2Z8cMRM97jmMm3e/ENMTTuRkNMx6WdjQeyNSzNz6WdRUppLOVvZiRPKhk17lPo6WMR9XXJZM5Fo1jFpX1G769s2KUxR6mthEnniiCKFSUeleX1FJkHJ7CIXSpNKOUu2R+lnLuXMpp8FSa4ugTVKcZWkIugfppZJKsbjFxO12zSzltkhYbqrGcOln7n+mTGcKLnMJb6alLMiSDQrTUJZMRBfn8pmd3vTJaINxmN0BoNz1MYZmnO0hgjm6NLi0sxcAljZDvonzmdpkvTZzCGSuAy5DWbSz9IpZ1G6ayLFVZIWBsmsUVqr5FPOouSyeYm+kjS/EfeP013jMQYV76co/cwlpTXNAdIIDjTX14nqHlffuBopqmVser2JBhwJahNX96xI1DeSNNiO9nX8vrhU+3qSXIPUR1P3tNrxfmoH13Z3Pk8Sdx30rd5VUpw+max7ouunrWNMAnsm7T6T4iqZusfUN02XSB/VPbavefFBjVMEx/REX1M7RTWLS3G1qa9Bex31jWNqFjtG1N+ku9q6J6hxXH1TmHO0E4zh6h6X+hrVPfZ0drs02E+unMp8rsh8jpnYqGlPcJ/Lwp3iymVTy0R1j6tvhkyq/Zygxpljklld3bMgTLVP1j1Be7a+mZeoe3yKa3xhju67Q6ZGiuob57F1T5S2m8U36QAAAAAAAIA+4yEdAAAAAAAA0Gc8pAMAAAAAAAD6jId0AAAAAAAAQJ/xkA4AAAAAAADoM9Jd+2yzxri2aEyeAOKyRVzKWZTUNW6ex46Y9LMo0WykE6ec+WSV6u3ZNLiITSjrxK+xFbz28U485zGToBSlDjZckpNLbIym3cOUs0zqqwu2sSlnNt01SCgL2iSffhalkbkUV5d+1hgLklkTKa6SpPHuxKAiaJPkk8ui9DOTlFba1NegPZty5kTJZS6d2CSURumuRZDWKkkyKWdh+6C5bZl93QjO/9JcE+zFNh45bC1N6luUaGYuvypdylniHK0lxTUr2KgLwXOXuEZwTR0wKVmDJiltMLinDkUR05LmmvSzuc3uNDKXZrawORK2Lwjas2lmrv/8ont+UfKZ5O/Fc4L2hjlAhsyB1qt/6XWn4piJBuwER5Srb1za/UiwT1eUJq3VJMcNJRLlmpmEbsXJrC2zP1omqXqs0f3aB4M2yZ93rXZwbTeJue78j5MSs4VPNHDcnEl9tddUm/pq2oNDwVyGwhTXifag7hk152Imkb6Ousel17v2oMYpXY3kEtg3kLqncEn1UY3j6ptxc/UMapzC1Nylq8WDNFi3S1PM/nCJrVF7p2lS7d1r7HT3z6a49jTVPjxsqqfXS/FnSZdqP2iuy1Hq65xmru6Jklx7WffYeqgw6a6NIMHW1D1RfSNJg0F10TQHiE21D1vXVpr3L4Nv0gEAAAAAAAB9xkM6AAAAAAAAoM94SAcAAAAAAAD0GQ/pAAAAAAAAgD4jOKLPNm0U2swEDazWNgsX+oWVu38yHi22Kr+w4kiwcvlKmcUqzSrnDTPDZmJ1dheOES0S7RaDdsERY0H7aBGfEgNmAchGZqH0GtZETjFBGpkFVO0CyjZQwvQP2t0izG6B3ChoogjCJCbaEwsoj5mVnE0YRBgSYfqWLTN21J4NjogWXO7lAsouEMUFJYQLKMfnV9EyYw8GYTVu9WPbHjeH2zOvpREsdFw2qx+nE/2DqdmACHPuBq/R9bULImcXVs4IQ3DMwtbumhosrNxMLLYsxQsoD5mLkw+O6G6f1+gOFpisPVoUeWEjXmw5HSgRBkdUD4iQpMHg/YrjoaRBc24061jsP+DqHlfLRO/iHFObuECJwaB/w2zPLSidCcHqmNXW26Z9PLimRgFYkjTajq+10bkxau4x7ryLzlF3PtvAoRoOm/DSUktglrlm2UAJs5+Cdnd/SNU9rWTdMxoFRyTqGymuccbja6etWTJ1j6tlNpC6RyY4Qpm6x7UHr8ftj6ivk/1s0ghCwVzdU7RNzRK8ve5zgrnUxuddOiDCXXB6k6SV3tfBxN3nzgFzIYrCsTL1jSTNCwKz5pngiCggQpI2aazqHsPUMZuYesiFQcwLgyNcQET83kZ1z5Cr/c2Nqkrd0zL3zwy+SQcAAAAAAAD0GQ/pAAAAAAAAgD7jIR0AAAAAAADQZzykAwAAAAAAAPqMh3QAAAAAAABAn5Hu2mdzG4Oa14jTuVZrm/SejoktHA5SYkZNbM6YGbupIIXJPdI1iTztIBloonv3QGMuoayIs+bGi+75DZqElwETJRQl57i0wCh5R5IaQf/CxQvZ9ri5ZzLTSyeluVS0KJUyl3IWtdu+mdTXGlLOyrE4iShMM5NJfbXpriYKK9p/Lh6rDi6lL5FyVpj5laVJOas0sXX0ja5DLZPkZM7/Mkh3jRJfJ/q6c6C7fzY9uY6Us1pSXDOSKWfR9dNdl336Wfc5M9yMz8WoryTNiVLEgjZJGjbtUaLZHJsQWz3FVZLmBffAOebAmWP29ZzgnG6YN2ZQ8T262aPYclf3jJv3qxP0HzHFScMk1TcySX/msucS6aP2MXPdc0n1o0H/VaZGcsf7SJD66s4jd95F56g9DKa5vklf92pItc+0Z1JcpTjJteHSXV17UEOkUlylMMk1lV4vqYzGzqbaT3eNYz7HhCmuktQMPhOYWrd0NXDFNsmX6FFSvUuvL0wNl6q5E4nI7nyxCazROZr8bNJTmVT7xGdJ97kzUw+5z7/u8/JwWFfE5/McU5tEdZJLr8+kuErS/GAuw+awcYmtw8F57tJaG+ZGX6XumeuSoRP4Jh0AAAAAAADQZzykAwAAAAAAAPqMh3QAAAAAAABAn/GQDgAAAAAAAOgzHtIBAAAAAAAAfUa6a58Na1DDxTrSXU1CWcfE2IyXQWqLDSKpPvagSVVyKTGDJkFtMEhnGXJjmFSZZhAP5ObRdGk6wWu0Ka62vbvNprtu6GpId80kTZq3y6ezhUlTJs3MJAOqEyRvuRQx2x5MvGNejEtmDdptimsm9dW97qwouSh5XEeXHJtEZtrLIIWpCNLTJIXvraT4fTQJZW6MohP8e1YN50Y6oSwzxgyVCYN0aZzR9dr1dfeN+P4QHx9D7j4V3UfNfXHI3i/dNrvbXYrroEtsDc6vQVOPuJSzRvCONV0StNEOagubbGeOkPFg/7WTJ2knqIfGTY00FNVZ8u9vdCy448YdZ9ExmTl+J8aufm440x1I30u1pNpn+mfHjlIz7b3ODBLdA90Ypj1MW82kuEqm7nHbM/VQ9Dmkh3VPYVJcU7WMS4J0ybHtKGXWJLC6RPqo1jX7ydbRwTx8fW7ae/W5YobWPTbdNTgU0p9Hg53iPv/6z8vRZ+v4fM58brfPCezzA/e8IZpH9RTXibG7z2lf98RjV6lxhmv4HhzfpAMAAAAAAAD6jId0AAAAAAAAQJ/xkA4AAAAAAADoMx7SAQAAAAAAAH1GcMQM4BYo7JiFi5vBIorjZpHNzFNatwBls65FWwN+4Uyz6G3UN7OA8oa+CvtsWrHZiRZKzqrjmDRjlFG7mXPYV4oXZ3Zj1LE/ssJ5m3OujA/KMmhPH77RItFuwWs7Rm/2nw0nmU3MezvdaRUbynU5Wvx/0v7BfSpqk/w9zW0zde+u4cYxlQWU1yUcwwRVZc5+HzJRfZ/a/e/er2R72HcDWRV92s87nyA0rdPoh57eTzJju/urOR9z0wjm4ca17dN7TLr6ywbEBee5qwPDUBDTv5bjo44x+lGPTrcN/HKTCZTIhgJF97Ve3gfcvc49b4hkK5DoGYnv29/vsvFNOgAAAAAAAKDPeEgHAAAAAAAA9BkP6QAAAAAAAIA+4yEdAAAAAAAA0Gc8pAMAAAAAAAD6jHTXPmuro/Y6kqs6Jv2kY3LOxoPU145J9Rk37e2g2c2zbaJwOmX8DLgdtLfN8+Lxsll5bDdGx86vuz1qm8y0Bx1N9/ayKUemf7hb3dgN84Po7XV9XXpP0F6aMYpmfDwVrWAM17dtkk+bwXHtEsAUnwNxuPPUE9gmNtr9egq7r81rj96Dhvl3IZegFO0nO0biWHD/PJU5bmxfM3Ydpjt1zCVsTfN1yF2X3bU909fdN+K+ZozMvc70dfe6dhGnuI8Hcxk0b0zbppkGaXAuAdDMI6otXBKsE9U4rr5pu5olGMO97viVxPvUvef2/arhWHDbjPvm6p5M32w9NGXuumKT/no0vz5cw/39JGh3fTP3NTtG4l7cNvWGGSOqCUpXP8SnV5y2WkPy7MRGg7onet2Srz2D/mEtJEmuboz6m76ufs2953XU3HFzqvZ3prvu6UOAbeazpK2Hos+0ybon+mxt73VmjOi+Zp8HmGPSP2/o3lEurd0lx44HFcCgOcbGzftSpcZp1/B5jG/SAQAAAAAAAH3GQzoAAAAAAACgz3hIBwAAAAAAAPQZD+kAAAAAAACAPuMhHQAAAAAAANBnpLv22WjZ0uij0spcollk3CQaheknJhFtxCSXjARJLCvL+HAZse2DYft40H+k4/rGqTJR/1bHpOZ14jFaYQJNnNjixi6D/lHbxA+mO6LIcNPIpDClE8qCvs344HNpVVG77Ttg/v0heB+LtknvcsllA93Hr9tNNpMuSIi1qWrtOIswSp8tXTKji45ySV3R9mxCWbz/wlQ0l5QW7FNJKqL2AZN2a9rD/mZfu+MmPvbizdlzI+jvxrDt0Xm0gV9uHHeoRs3umuqv193vuUsLi/pK0min+9gbNfepOcV42B7dAwfLVti3ae7ng3GMczyG4rFdQvycoN5wfV2aWbNH0XsumdUm1Qf9R2yqfbzNkaDeiNom2uNjIdPu6ht3nEXHpDt+M8e7O4/ceRefo2HXDUbmOmn7Zq/XUXt67ODeY+5TRcsMEtwDXa1QdOLrUDTtMGnV9J34QfATU9/Y9syBlqx7Uon0tu4J+mfqG9c/OY/oGMnUN67d9jXlV/SeZz4n2PbkZ5Npl72mBu3uc6ete6JkVvf514wdfbYebpj6JlEPufuiq4e8qL97bpJ5RmLqL5M+HWk85ntvo4nazY8JAAAAAAAAoK94SAcAAAAAAAD0GQ/pAAAAAAAAgD7jIR0AAAAAAADQZwRH9NmyzriGOpMvLuh+6pZKjBZFHjOr0o6b57TxAsrxwo8rO8Nh+wrTvqwzp/IYK9txe2ZR79HEwplukU23CHO02Gen4/qGzT5dYKoKG1sQT6OOBZTN4rHRbi2a8eANEyjRGQxCH8zCwHa94OhNqJ7VMrHNaNHbaLFgmfAESWXQXtgFlM0Eg4VO3f5IixYHdguouvCJRvAaXbiDC5QY7D7PbUBE0FeSyqFgjEH3vpiFeoNjL2qbaHcLK1cPn8gsIm7/qa0fgRLhivKurwmDCK7LbbPIsWsfCy44Y8E9Q5JGzf03use4hZIHi/g+1Qiuwc3sBccYDxZnHjfn6Jwifo0jRfdcBs38zGEdcmuIO5nllV3ow3hwwNvgKdsevOem7nH1TaYecn3dcRbXPfFx7Y736NzInl/ROeqLhbi5V3WPvb659jDQxy2k78KuzFyCa767P0ShVpIUnKL2Pu/ua+G4pt1WjUFwgQ2TcuFOre4zvXCff1woRbTAez/qHlMrRLWMD4ioXg+Vrr6x7UFAjKt7bC0ThE/YGt/cz4P+bow6ArMs+3moBsHQNsPNfj7sbnefOzNhEO7zr/+83B3M4D6HN83VohldtJLajdG4PbrPmwpixIRBDEV1jzk+GmaMKlfaZSaEJ4Nv0gEAAAAAAAB9xkM6AAAAAAAAoM94SAcAAAAAAAD0GQ/pAAAAAAAAgD7jIR0AAAAAAADQZ6S79tkDnaY6Jn1lXTom+WUsePbqUs7GgzQzKU40cylnLoksle5q0mOWJ9pXtc38Eu3jbZPuapK3ojSdKKVn4gdxcyoRMSGbzBqnnJkxEimuklQMdE/GBQC1h9xGo1SvuK87ozpF9/HecAllJq0qTOQKUsskSSaxtWgF551NcTUHQ5hUW09yZChId5u8PUrBzaWqRvvaJ7NWTy7rDMV9fXtwngfHtCR1bHv1vu78Sp2jyZTDXqW+libNzN2/ovSzVtukk5vzP0o/G2nF94EBcyFqBO1RWutkoiQye9+2iaPxvOcH6WdzgsRXSRo06a5DQSqa6+uS3LL7pCq3n6J9KsU1zpi5E/jU1+597eoYV/e4xNZpr3vM8R6dG+48cudddI6698ud/3UIN5m8voXXz2zdY67jUQprYfaHOe3UsQVbN/caw9BSU/cUA2Z748EETTppYWqZMMHeJbO6Wma6axxXH9q6MaoVkqn2Ud1j3hdXD0W1jEtxbbt01+C4dvW5r3uCa4VLd82kvrpzNFsPZdTwGc59Pow+S7rPne5z6kij+5o/0HBp7dXrnqzoHt02b0zmXizFNY6ve7qTaifau69DQybVfir1zYMuuTqBb9IBAAAAAAAAfcZDOgAAAAAAAKDPeEgHAAAAAAAA9BkP6QAAAAAAAIA+4yEdAAAAAAAA0Geku/bZktYCjUVJj4/SMc9SXcpZlKLiUlxdskqUNDdSDoV9R23KWdw/al/ZjvuuMu0rgvYVLbM9074qSD8bNak54y2T+hqkn3VMIppNOasjdiiTcpZIRbJpSyb5tByIk3CiQC4XIlTY/RElb8U93fwaUfKWSRYtTGpWoxVMPGqTSTOTpHYimdWln0Xpri4Jtg6JNDPXv3RJsC65LHgffbqrSR2L0l1dEplLORsKEsoGzfXXtIcpZ+bSn0lPtuF/2XTnzGXIppwFgyRTH6Nrasu8t+7e2QwSuRo1RGa7pEV3Hx0N3uDRZny/nNeIEz2HG3Fy2bIo5cz0dSlnQ0HKmetr011NKtpUZeueqMZxibk21T6oZVzKXLruCWoW13dFKz4WoiTX5eOmr0l3HQ3OGXceRYmDUnyOpuub8FoRd7WiIWq4HtqUepvuavpH+yT5GsO6rDD3S1fbBffXxrirWUxia5DuWkR1jKTS1UNRjeNqlmzafR2iGsfVPVFkrhQmxNp62aa7RjWLec9Ne1j3mNrEJdVHSa5uDHcORO3m8jvJfgrakumutr6pJfU1qFNdqr35fBhdU93nzlGT2FoU3df8bDppJ0ykT9Y9wRs8GiTPStLKZN0zpxjraovqGMkn1Uc1TtN8IJ1K3fNQK66lMvgmHQAAAAAAANBnPKQDAAAAAAAA+oyHdAAAAAAAAECf8ZAOAAAAAAAA6DOCI/rsntYWWjYeL6i4Wjv5LDVazDG78GPUHi2qLMWLY0/e3j1OtCDyZO0j7e6xR9wYZgHlkWCx5NFxE7BhFtMNgyPMgqF2YeUoVCEbJpEIjsgsoOwWeHUZB36R6GjxTbOfzOEerR3vwgIaLbOgcRBs0QkWx5X8wsqdIMihMNsrXOhDNIZZKNmu+7qhBEeY9jD0wSy27NoVjNEx77lb/LgTzcMGR0w99MEGRwSXIXNJnWSB56DNhrvE7faAqmMB5Yhb/zuxsHLLLKA8Zl5LkVgsOVoo2bWPmRXlR037qiAkYnk7Xih5btMslGzDILoXRXaLLfsFlLvbm+YNcwtQu/5T5eqejrnHRP0z9Y0U1yaur6uHXM0S1UO+vqnengmIkOK6Z8ycX+68i87RTH1Tm7DuMfd+tyh9lGXg6h5bl1W/3pT2Pmr6B/fGsunqG3MtC+qThql7XC3TiGpgExzh6p6o3YVP2FpmAwmOsKVuVJ8k657ofXR9XdhVdNnK1Deu3Qdjhc1x3WO3F48R1j11BUQk6h73uawMjmtb3yQCs8YbLsSljhCseB5RWNBoM35j5poQrFXN7kCk5aY2GW7E4QquPaqHMvWNJDWCkAjX16lS9yxrxa85g2/SAQAAAAAAAH3GQzoAAAAAAACgz3hIBwAAAAAAAPQZD+kAAAAAAACAPuMhHQAAAAAAANBnpLv22d1jm2v+2B/S3to2sqZbJuXM9XXJZa0grW7czC3qK/l01yghbyxIa3V9JWk06O/SzMZMMut4kFzmUlzHx6unu5ZtExfk2qP3poYAK58+ZZLBgtQhF06qyQOJuwUJWWXDpbDFQ0TpbCYASIXZ143u0CGbLtYwYxTBNn2aWTy/KKHMzcOnu5r2cIzcARUm0CUOX0lhopk9Jl26a3DJcclgNiktSjkzSX82XSxKmXVpZrY9Spk1fV3qa9DfXCJt+pltjxqzia/R8e7Srs351QmOvXZh0smL+MVHr6Vt5jEepJlJ8b3ApZytDNLMJGlOkNg61IhTxAZM+7BpHwzaB5IpZwON7jcsm+4aadoLX6xXdU+UVCdNUvcE7ePJRF9XD0W1jEtxdXVIWPeYvi6xNRrbpbi2WyZhNzh3bd1TR+qru29EbcnrXvh2peuv6gmgvu4xNUtUb5ggwoZJ3mwEtYUbw9UyjehYMO+hrYei5kxfKXxvaqlvnGTdE76PLiHWHquJWsGlFkc1SzIJPpNq75NZE2O4uieYn/+cYM6vOtLr7TEZDG6uh65ObTem/p2pMphHO1HfSNLoQPebM9KMP2itjD5QSRoK+ru6x6W4unpoMKgtfF+X7tr9Rvai7lkxNlr5dx2+SQcAAAAAAAD0GQ/pAAAAAAAAgD7jIR0AAAAAAADQZzykAwAAAAAAAPqMh3QAAAAAAABAn5Hu2me/XrW55qycO2mfTiL5bKJ/kKZjIopc+lk0RpR8NtkYLlUmSjmzCTSmvRUllAVJq5PNL0o0a7vEQZOE02lFKWcuRiwTURZ3daJwIROI6BOlEilnNpjNJl4FXV0qlUt+ChPKTAKrSz8L2gvzvti01eDFuwCgVLqrS0rLHAs1pAJPbDSxSdc3SrbLpruG6XhxV5tyFvRPnQOm3aeqmTEyyaw2yS2xPZdy5m4nmfQz80ZGx2rpzgGXfhZcuMzpHKaZSVInOKfbzVzK2Vize6sjJt21GaSkStJgY7i7LRhX8omoUQKrJA0kUs5cQlnU3jAXET9GLsm1Klf3uHTXqMaxfU17lMzasqn2pu4x/aPjzNU3PpGvu90dv1HyvGtvm3TXKMVVMjVOMt21CKNPkxGMifuDux5G20zXNy6cPLxe55LF42TWXN0T1TjutI22Z/vbuieR2OrGmKF1T3j8ZdNdU8d19bGz6a5Ruz1+E7WMTZNNzcP0dfVNYl9bicRhn2pv7mthZzOGmUdU49iE8+jiJGm01f2GDZiaZdDUJlE95NJdXe3k657ucVxt4saIapxMiutE/3XXPSOrVqXGDLcz5REAAAAAAAAATAkP6QAAAAAAAIA+4yEdAAAAAAAA0Gc8pAMAAAAAAAD6jOCIPvvNis00NDhvvX7XLZodyS6gHC2h6BYzdmO0zcKZnWCctl0A3GwzWhjcLMhZJubhFkr2i4AG7dkFlKP27EK4waKXpVmsVm7R1mCjhRkjG0pRRMERNiDCjJFY/NgvoBy1VV9Q2o6RWVBWPqwiNUZdiyVXnUZyTe9McIRdsDmzgLJdGDgYxIZPmHnUET7RqzFcQETitUhKLZrtjr0wJMKN2zLX62BVZHcNLzvxBSC6P7QaJvjALC481ui+QDWbZiFiN3awo9xCyU0zRmF2djSOW/zYvQVhcIQbY7ovOIYNC4kCABL1jevv6h43D1f3ROPYusyNEdQ4rq+rnaJzKaxjNEkIVqruiZuj9vQhFm0yeU0NjwZT95hd6uuh6DUOmjFMzRK+vzb0IW6P97ULiKhe9/iF9KsHC2Vr3Q2m7smUcMl7biY4wo4RhHFtMDVLDcFd2brHv4+JAypRz9vttczQwU7pmHPUfU7tBPukbeqsMVNvRDWOrW9cLRPVFcm6x9UhUd2TqW8ma49Mpe4ZW9EdHJbFN+kAAAAAAACAPuMhHQAAAAAAANBnPKQDAAAAAAAA+oyHdAAAAAAAAECf8ZAOAAAAAAAA6DPSXfts6fL5GmjOX6/fNcEvpq9J70pENrrt2bFdIl/Q34ZEuRTM4Bd8AmAiQskmkSUSW7OJnmFqVjJKM5MGlRrDJICZBCWbtppIKHNjhCFsrq8RJ7NWP8bsPLLprjWknGW2V4sa0l2dTIJaTxNi65hHIn3WJ9WaNzIzRnp+UUK0GcMJz1GXsOkO1urXirIdX4jC/Wf2acddzIL+hU1gjYeIksFsWphNVc2MXb2vk00zc9ucqkx9M9E/U8u42mRqfSd+kKh7kmOHtYyts+Ixwv6Z2kRK1U627snUr/ZaFg2SS/SNklkLkxyZfY1xmryZiBs7USvYS0umxqyhZkldQrI1y4ZU41QcI1uzhP1rGCP7mSDsb9NkzSBhzWKSWRO1k61jsvPLvOdu8E5wL3bXocR1uWy7OjBuj8K42yZ62tUyrajd1RWuHoraXH2TTLXP1Bt+jEw9VH17j9VaPvULCt+kAwAAAAAAAPqMh3QAAAAAAABAn/GQDgAAAAAAAOgzHtIBAAAAAAAAfcZDOgAAAAAAAKDPKqe7djodNRrr90xv2bJlWrhw4Xr97my38pFhNTWn9xvKxvT1NJ2peiJqapt1vEabVpVIwkomfUbzTgbshS/dJdiUjcRrSSYl2XcgjM0zY2T3X6JvnFCWO256mlzWqzH6YbqT0uwY1XdgKoUtk8yWHcMMUc883MnhNhoNboaIEsoyiYNS/E+HJmKrtK8luiDGXVOnc3bfZfZpHe9LLefcDL3gpC8M0Rg19M3Mo4ZaIa1XdU8mCVY1JZwnEultTF9UmriPOvYUdXG8iTEyx0Itqao11D3UN9M2Rh11TzZlNtN3ptY9mf3qzq8wddsUPjbtOo5EjbeXqRWS+zocO33sRW3JE72Xdc801TjtR1pTHqPyU7fnP//5uvnmm9Mb+M///E/tvPPO6d8DAAAAAAAANhaVH9Jdc8012n333XXCCSeodP80/igPPPCAXve61+mQQw7RPffcM6VJAgAAAAAAALNZ5Yd0RVFoZGRE733ve/XHf/zHuu2222zfr33ta3rGM56hs88+W2VZamCg8l/VAgAAAAAAABudyg/prrrqKu22224qy1JXXXWVnv3sZ+ukk05aq89dd92l/fffX294wxu0dOlSlWWpPffcUz/+8Y9rnzgAAAAAAAAwW1T+ituee+6pn/zkJzrppJP0oQ99SMuWLdPf/d3f6bzzztOpp56q//qv/9L73vc+rVixQmVZauHChfrYxz6md7zjHSrcgq1QY9mgGuVgv6cxu83UhWmnyC5+bHbIhrKbpn/d8g3llWM2qGNt9xkrWnM4tRq0pHb1zRV1rN69Mb9fmL02gttaeK316VWJ1ulH3YOZbKOoe1KBLb3bIdQ9M0O5fOrPdlJxrY1GQ+9617t000036bWvfa3KstRll12mnXbaSe985zu1fPlylWWpww8/XDfddJP+5m/+hgd0AAAAAAAAwDqkHtKttt122+mss87SV7/6VZVluea/7bffXpdccom+8Y1vaLvttqt7rgAAAAAAAMCstF4P6Tqdjk488cQ1f8q6Ou31nnvu0Wmnnabf/va3tU4SAAAAAAAAmM3SD+muueYaPfe5z9Wxxx6rRx55RM1mU+9973v1J3/yJyrLUl//+tf1tKc9TZ/85CfVarV6MWcAAAAAAABgVqn8kG7VqlU67rjjtNdee+m6665TWZbaddddddVVV+njH/+4Lr74Yn32s5/V/PnztXz5cr3//e/Xbrvtpssvv7yX8wcAAAAAAABmvMrprjvvvLN+/etfqyxLDQwM6H3ve5/+8R//UYODf0iveNvb3qb9999fRx99tC699FLddNNN2nvvvXXkkUfq9NNP78kLmOkGlhUa6PQ+XKOW5J26phmMk55fpr/pW0sgT+a1uPiuXu0P97oz86vrfQlfYzJdLJXkVoPM/Dak82tDVkegXHaMOi5+mW3aBLDEPNz2onbbN/e6w8M9M4/0GNXnZ0/FRHv6KMi8FiN1iUsf18n+gelPlJxefal7Mqd5D+/zdptVx52kPTXvOuqeVF1RZVKPHiNxsci2B2zdU0NNm1JH/ZXeZg1jbOimu8bpZaxqLffA6CRNbM9J1hWpGqKO+3wd86uhJkiGXU+9r2p6Lc4MqnuKFVM/Nyt/k+5Xv/qVyrLUs571LP34xz/W8ccfv9YDutV22GEH/eAHP1jzrbqyLPXVr351yhMFAAAAAAAAZqvKD+kGBwf14Q9/WD/5yU+02267rbP/2972Nl1//fXae++9pzA9AAAAAAAAYPar/OeuP/7xj/XsZz87NfgOO+ygSy65RCeffHJ6YgAAAAAAAMDGovI36bIP6B7tmGOOWe/fBQAAAAAAAGa7yt+ke7QlS5boi1/8oq655hrdeeedkqQLL7xQ2223XZ1zAwAAAAAAADYKqYd0rVZLf/d3f6dTTz1VY2NjkqSyLFUUhUZHR7v6H3vssbrgggt0xhln6PnPf349M55lBpcVGmivZwJID9OWwgCabIpYJvHKfaczk+qVHKNsdEe82DmbsaNtFnafxj+Ix4jjZzKvMXp9kvz3ZzNjZBNso23WMHaRHCPcr7ZvYh6uayohNhc5ZOc3zcoaksjKxMXCbs6OXX2Q0iVtB2Pb123GCPvb7Zmxo9fYibsWpj3z2m0QYSZV1b1GN+/E++VeY9juUtUS+8m+7hpei5Pa171MiN1QkmB7mIKZqXvqSGZN1xuZ+dn7fFD3JOobKb73+Nfi6p7q91E3dtierAnCGiex7yYbO5xLpq/Mfd7MI1MPZeobM4Rl657EBWpG1jdS6sLg6p5wk5n6xgxi66xMPVTHGPaeW31sW1eYHVLaWqF6+qyvFYLaKVMTmG3a7SXG9vvJjBHWTmYe2f0U9a0jwbaOtODE9qoqVk5juqskHXbYYfr85z+v0dFRlWWphQsXTtr//vvv1+23367zzjtvSpMEAAAAAAAAZrPKD+nOOecc/ed//qckab/99tN1112nhx56aNLf+dM//VOVZamrrrpqSpMEAAAAAAAAZrPKf+56+umnS5Je8IIX6D/+4z80MLDuX33CE54gSbr11lvXc3oAAAAAAADA7Ff5m3Q///nPVRSFjj766EoP6CTpcY97nCSt8xt3AAAAAAAAwMas8kO6Bx54QJK02267VR78vvvuS08IAAAAAAAA2NhU/nPXrbfeWr/5zW9099136znPeU6l37n22msl/eEbdeg2tKzUQGs940My6amJlNSJ/lEclOmbTlWd3jF8ezC469t0YwRJaaZvJqbFjuH6p1LfakhgbboxTDJY0N8lkbn2RiLdtWHbu2OHXIpY1FeKd2vT9TWxQybwLjVGPG7vIhg7yTjDKEms4xLAEqmvbh6uPRzDJF5l2stOfELb1xIcImXbjOESudpBWzJ9MpeKlnvPw7GT6bNF8BrdnKO+buxs2m0jGjv7WsJE3xrGkOJ521BFd+Iltuf06pLT02RWd3JkxjDt7j4f1Sw1jNFxtYIdO7iW1VD3yPWNm1UEEyxdXeFEB6ut9xI1i9v/Zn62lmkG9YZNzHVjd4/h65vq7b42SdRfbs6JuidT37h51CVT42QS6bN1TzuoLWywuK1Dqvd1dU9Uh9g6KzFG2XaJ767uCcZwu9/eR838qn/UtTL1Ri/rnui1p+dRQ+3UqzFc/0x9M9HfjJ0Yo4rGqqlfqyo/Nnje854nSbrooosq9V+2bJn+7d/+TUVR6AUveMH6zQ4AAAAAAADYCFR+SPeGN7xBZVnq5JNP1ve+971J+y5dulSvf/3rdc8990iS3vjGN05tlgAAAAAAAMAsVvnPXQ866CDtv//+uvDCC3XQQQfpjW98o/7iL/5izc9/8Ytf6JprrtHPf/5zff7zn9eDDz6ooij0V3/1V3r5y1/ek8kDAAAAAAAAs0Hlh3SS9M1vflOvfOUrdckll+hLX/qSvvSlL6n4/doer3zlK9f0K3//98G77rqrPv3pT9c4XQAAAAAAAGD2ST2kmz9/vi6++GKdd955et/73qdf/vKXYb8FCxboQx/6kN71rnep2UyugL+RGVxeavDRwRHZFSsDdSygHC6ym1zkONPuFy5284v6ujGqt9tFmN0iltH83FqRA24Vy2AIuypqYiXMxOLTkuJFkd1CycGCyBPtcf9muIBy9b6ufcAENrjQh2awzcFmvAprM7GAsguOcIscN4NVUbMLImf6u77ZMIjMGFG7O6475uLSChY6jhZVnmwe7aC91Y5P9LZZqDeah1uEuW3CINrB4scdc764hZXLRrCgdCsZ6GHaw0WKk5ehTBCBW/y4Ebx2GxDRMu3RwsV2e3F7dM13fYu2uU5mFlCuY4HnTMiEzILLdSy23EOZECwfEGHOu2iMZH3jgwsSY5u6J6pPgkuC7TsxdmIebj3uYH4de+AkxvYpE5XbUwERUljjFKZWy9Y9UR3STNRIbgxX90T1jSQNBDVOpr6ZGDtRf5k3MgqJyARVZNVR92RrpKh/25xgtmaJ6o1EfSNJ40GN4+obV8uE9Zetb6qHYLXNRcvWPVHairtv2wuO6Z4JzDL7OjrM7H3b1GthcEQiGEuK6xM7RqLdh11Mve7Jtoc1pg0L6WFgVoVxi5GpD5p6SLfaIYccooMPPlg/+9nPdMstt+iWW27RkiVLtOWWW+pZz3qWXv7yl2vLLbec8uQAAAAAAACAjcF6PaSTpGazqec+97l67nOfW+d8AAAAAAAAgI1O5XRXAAAAAAAAAL3BQzoAAAAAAACgz9b8uetJJ52kE088sScbKYpCt912W0/GBgAAAAAAAGa6NQ/pHnroIS1evFhFUah0iRiPUjwqNWt1/9Vtj/39wiVsQUMrOhpsuRiTCTZFLLFby4ZLzTEpZ9E2zftoU8RsylmQ+mjGsOmHUUKZG2PAvPaov1mlsTTt4TvngoHc2xzsJ5c0k8mKce+tHTxIBnNpZg2TftYwSakDA1HKWTy2S1sdCPq7voMmdSzqP2giG11y2UAQf+T7Vk99bZi+NoUtcTRkk9JS6WfmgA9TVc3FIpPu6sZomQtOlHKW6StJ48E8xlpmHg2THNsMXksrfl/acvup+xgp7YU2mXJWx206k7xlkubCdDGT4mrTVoP+Nt3VJcQGyWU+3bV6eyORiDbp2J3geu0u+ZmktEQimmO356aR+ZuOTFJ9pr6RS3dN1A+TjN0JElF98rxLF436mh3iapngvUnVN4at84PjVJokpTcjGsONm0i1t3WPS2YdSNQ9rlYwtcxQMHamvpHiOiTTV5IGgotfNt01GmM21TcTY1dPZnV1SFj3mItF1FeSxoM6KVPfuP4tk+Lqxg5TX13qbnSBU7LuMRct95YXicInlZ6eTTMNawXT19VDqbrH3GMSCbGZeiiqVybac2NHabw+3dW0m7nE23NjrPtXi9FkIRRYc2veZJNN9IQnPGGdv3D//fdrbGxMZVmqKAptuummGhiYGGZkZETLly+fmFxR6OCDD9aOO+6oBQsWTHmiAAAAAAAAwGy15lH0u971Lv3617+e9L8PfvCDarfb2mKLLXTaaadpxYoVevDBB7VkyRItWbJEjzzyiH7729/qYx/7mObOnasrr7xShxxyiD784Q/38SUCAAAAAAAAG7bKf2RwzTXX6G1ve5vmzZunq666SkcddZTmzJnT1e9xj3uc3v/+9+sHP/iBli1bpsMOO0xLliypddIAAAAAAADAbFL5Id1JJ52kVqulv//7v9eOO+64zv7Pf/7z9Td/8ze677779LnPfW5KkwQAAAAAAABms8oP6S6//HIVRaG999678uAHHnigJOn888/PzgsAAAAAAADYaJhMp2733nuvJKnTqZ5WsTpQ4s4778zNaiPSXNXWQNvFmExIp2BF6ak23dWNEfQ1Y9j0VJegFgT4uITYIpHM2nFHs0mrK1MJT7Ew0Mwl3pj9F6WfpdLupFpSzoow5WzqKa4T7UFCWSLNTJKGg/7DzTjmyI4dxB8NJfpK0mAQ6xOllrm+rn/THGUuuSxqd2PUoW0OHJdcFrWPu4Qy0x4llLlEtNF2fAFoDQTJrKbvSCNuj+bhUunGgpRkSRo3abAZZZD66tI4S5dglbqY5YQjbCgpZ7avOb+C/pkkWCmed6bvZP3rSFALk9Jc8JkbOxesWJmvTRK1jDmmc8nzLmE+l3YfhT67vm7sIrg8dUzSutuB0eUiETw/0T8K0nX71J3nUVKt2Z4VvufV0+sn5lG97smkuEomTd6M4eqQ4YHui45LpJ9j6qGhoN2luLqaKkyINfMYMHGQg0F7pr6Rpr/GydQ3bgyXSD9u012DZFbTd8x88BkL0lbHmqa+Mcmso0E9NGrS6wvzfo2bxNaIub2qjOpDW/e4z1lm8EzZk0j6zKTX2zES9Y1rLxL1jR0jmRCbqXtsuqtLn43qkDrSXTOJr5OMvdaQY5M/26mi8qOALbbYQpL0f//3f5UHv/LKKyVJzebUP5wAAAAAAAAAs1Xlh3QvfOELVZalTjzxRN11113r7H///ffr85//vIqi0DOf+cwpTRIAAAAAAACYzSo/pHvzm98sSVqyZIle8YpX6IorrrB9b7jhBh1wwAFavHixJOmNb3zjFKcJAAAAAAAAzF6V16T70z/9Ux1xxBH62te+pltuuUUvfelLtddee+mFL3yhnvKUp2hoaEh33nmnrr76al100UWSpKIodOihh+qoo47q1fwBAAAAAACAGa/yQzpJ+tKXvqTh4WGddtppKstSV111la666qqufuXvF3Tca6+9dNppp9Uz01mqOdJWs/2olRoTi3fbhZWDMdxCxG7RzGgBZRcQUbRcuINpDxY6dgERblHJcK3UZMJGJ1gFtDD732U+RN9FdeuWutfSo3W3/YKodjXNYH+YxZYbpr3ZrL6AsguImBMslCzFCxpnFkqe6D/ePa5Z/HjYrKw63OgeY6BhXndiAWXXt2kOnIZdDb43Ombh4rZp7wQHoAuOiIIZXP9Rs1DyXLP26ar2YFebWzTb7dNocWa3sHWGW3za7CZ1gouLDQUy56hdTD86/+0NwrUHm0sstuzGsAsDJwIlMgERkhSc5qmFkl1/GxzhFkq2r7F7p9gwDhv6ELSb+TnTHhzhBPVGaYMjqgdKlEHwjOQDJTJhV65vxy1indop8RjRddmVnTbAKjrPk4t3h+u+29dnFomPBre1mnm/gvaGuZ+7+qbp7v9BjeMCIuYOBhccxXWPC7Vy9VA0xtygFpJ86ENUD9n6xgZpBcER5r2dqXVPVLPYuscFaXW6xx7tdNcxE+3xvo4CHkbapr6JEmkUh0G4gIgMF9oX1TeS1AnO0fDclyYpOKYejJXh78Vxe3TK+DAp0x7UEJn6ZqJ/YoxELeNCMHygRGJse79MbNMEkaQO98eM0RyfxuAISRocHNQpp5yiq666Ss9//vNVlmX43/bbb68zzjhDV155pRYuXDjlSQIAAAAAAACzWeqbdKvtueee+uEPf6g777xTV155pe677z6Njo5qyy231LOf/WztscceJLoCAAAAAAAAFa3XQ7rVdthhB+2www41TQUAAAAAAADYOKX+3BUAAAAAAABA/ab0TTpJ+t3vfqdWy6wI+Chbb731VDcFAAAAAAAAzErph3R33XWXPv3pT+t73/uebrvtNnU6607aKYqi0oO8jVFjrK1GuY4EEJvG5wYNUs5MAquLLY3SzwqXvONSXE36kQlRyYn2iYlhcUluRTA9l5qTSTmzca2Z9p5Fvsqmn4UpZ2afunRX178Z9B80iWhDJtIoSkXLpLhKcaKZS3F16WdRQtkcE5VUR7qrSzNrBpGXdSSOSj51NDJexreRdibd1aSLRf3dfnKpr1F6nEuUy3D7qNU0172gre2uqWbs6LyzadLJYNZM31ROWiIJVorTz7JJaVHal7u2Z5LLMimuktQYj1LOXEKsSRdzqa/RawwSX11fSXFKp0kAnErK2TrVkGpfBuedTy11yfPdY5Smxi1NMqtLIjaXScOdeNEblqsPo/u8q+0y5519y5OJrakhIsnrYdRu6xvTPmASW6NU+6hN8omoUZKrS3F1NUtc95g0WZtqn0l3jceI+jfNcTDo4iADvax72uYPzlzqax3priNFkEhv6uWGuV43TE2VESVBu2TWdsOk4AbtdXyucOn16VTwjDo+w2Xas6n2Ue3kahaXzBqmzOYS6cO6J5FSP9k2o3rI1jeZVHvbNx66So3TqCHdNXUWf+tb39Ib3vAGjYyMSJLKWp62AAAAAAAAABu3yg/p7rrrrjUP6FY/nHviE5+obbfdVvPmzevZBAEAAAAAAIDZrvJDug9+8INatWqViqLQ61//ev3TP/2TnvrUp/ZybgAAAAAAAMBGofJDuiuuuEJFUejVr361vvKVr/RyTgAAAAAAAMBGxS2H3+Xee++VJB1xxBE9mwwAAAAAAACwMar8TbonP/nJuvnmm/XEJz6xl/OpxQ033KDvfve7+s53vqObbrpJ73jHO3T88cev1Wd8fFwnnHCCzjjjDC1evFhbb721Dj30UH34wx/WZpttNm1zLVodFS5C6/dKE0sVpZNKCuP+CpcAZpJZw+e3LgXPTcNE05TBxG3yixm8ngTAoM2EoRTJhMKUWtJue9XX7A+b4mqSSIN219eNPRAcJAPmzR1MtLs0M5dcFiW5uqS0OtJdoxRXSWoGr6Xh4qCSOvbi0m3QpFNHyWUuyS1K/5WkZjn11xOlsw2YZLZM+4CZm0thG2937w8TMmnPgfAy7uKnakq8S0klgMbNYYBl8iVm7g+5lNnqSauu3fZ1aWauf5CKZsfIJLaaNFObNNejALFs3RO9Fpv0ZxKYwxcZJL5KklzN4rZZS80SDJ44fl17+vyqI5F+ui9PiRfjk2Bd+qTpH7S5+0OmfcAkwWbaMymuUlz3ZOqbifbqCbFRfeP0su5pm5qgbS5EjaAuiNoma48nZ5rNwRelz7bMPnXtUc09VsSJtHXU/q3Mxd3JprvWkQbbo8+B/pLlaoVqbdl229eFl0Zj1JDiOjGXKDnW1Ddm7DCZ1dUxU6h73GvOqPxJbJ999pEk/ehHP5ryRnvlkUce0ete9zrtuuuu+od/+Af95Cc/0bbbbtv1YHFkZER/8id/ove///365S9/qV122UUPPfSQTjzxRD3vec/T0qVL+/QKAAAAAAAAsDGq/JDuXe96l+bMmaOPfvSjuu+++3o5p/Vy7bXXarfddtNZZ52l5zznOTrllFN033336YYbbtCb3vSmtfoee+yxuuKKK7Tjjjvq17/+tX7yk59o8eLF2mmnnXTrrbfqrW99a59eBQAAAAAAADZGlR/SPeUpT9Hpp5+uBx54QM9//vN13nnnqezRnzlkLV26VAcffLAWL16sP//zP9cPf/hDvelNb9Kmm27a1feOO+7QKaecIkn6xje+oe22206StOWWW+rrX/+6JOncc8/VTTfdNH0vAAAAAAAAABu1ymvSffazn9XDDz+sF7zgBbrssst02GGHaWhoSE960pM0MDD5MEVR6IYbbpjyZJ1jjz1Wd911l/bbbz+deeaZdv01STrzzDPVarW06667as8991zrZ3vssYd22WUX/eIXv9C3v/1tfeADH+jZnAEAAAAAAIDVKj+ke+c737nm4dfq/zs6Oqrbbrtt0t8ry3LSh2ZTtXjxYp155pkaGhrSZz7zmXVu66KLLpIkHXjggeHP99lnH/3iF7/QVVddVftcY6VfsPD33Csqo0WEJSlahN3tF7ModTiGm6ZdnTkWLvBsx3bt0Q96d5w5/VibvWcSu6+Hp7QadnHm7vZM34n2YFFfF4JhVkttBAelDUSw7cFiuomAiIl5BGPUdUAmFjSOghkk836Z19hwB18QPpEV7mv7nueOp7Cvec+jMdxC5LNJT19iDWOngg/S96loe7kxfJiRGSfTNxo7u1Byj/6iIl33ROeXGcXNOA5EiHuXpu5x+zocJns8hduzr6b6ID08RzeCS9wkgRLRPTC3Q6IxsvcvV1tMVVRPTd6euY9Wn3NP6x57IYqbO8EvdEyN1HG1QlRjJvdpptato+auQ+pzxfR/3OupDeY6GYYquL5xc+q1pD7jm3b3HMNuc7rqnqm/qZUf0kkK/7y133/y+rWvfU2tVkuHHHKInvSkJ+nUU0/V//7v/2pkZES77767Dj30UO20005r+l9//fWSpGc+85nheE9/+tMlSffee2/vJw8AAAAAAAAo8ZDugQce0Pz58zU0NNTTb8ZlXXzxxZKkXXfdVXvvvfda6bPf+ta39NGPflSf/exn9cY3vlHLly/XQw89JEnaaqutwvE233xzSdI999yTnsvChQvX2affDzUBAAAAAADQW4sWLVrn87Nly5at9b8rB0dsscUWGh4e3qAe0EnST3/6U0nSCSecoLIs9X//939avny5rrnmGr3iFa/QqlWr9OY3v1nXXnvtWi9+XQ/plixZok6nN18RBwAAAAAAAB4t9eeuG5pVq1Zp+fLlkqRtt91WF154obbccktJ0u67764LL7xQe+yxh6699lq9//3v11e/+tU1v9tsxmsdtdttSVKn00k/kHzsE9DI/fffr8c97nGpcQEAAAAAADBzLF68WFtvvXXqdyp/k25DtPpPV6WJb9KtfkC3WqPR0LHHHitJuvLKK7XVVltpeHhYkrR06dJwzN/97neSpAULFmxw3xoEAAAAAADA7NT1Tbqbb75ZX/jCF9Rut3XMMcdo5513liSdccYZU9rQ61//+in9fmTrrbfWwMCAWq2WnvSkJ4V9Vs//4Ycf1m9/+1ttu+22uvPOO/XAAw+E/Vc/pHviE59Y+3xjxVqRNmXmwaB7xBqM4QJYy6YZJJpH06WqmXk0TLJa1OzGsO3VX2Mv9WOboTqWOsyEHNrAm6nvkI4ZI2rP9J1o7z7eXd920FeSOsGxN25SSG3KWXTyZt/DYNo2rDmZ7halkbn94dJd4/fL9TVjR++XuSi4eUTtbfue546nsK+ZXzRGHefLhs69RBuaV0fCW+Ie4++5UTTr1Odh90d2fmFasBnD3LuLdtBoApXLqK/ZaCoxV8m6x9YhwfmVrFmi98bv/+pj2P7Z4yncXvIEy8zDbrM3fWcqdx0P74HJnR2Nkb1/hfdAmzhavR6y93NTzHSiOqQwn0HK+A+9wpTTRBr9xDyqfz/F1T2u5ov3dbxP3RjRe2Dr0Rpq3Tpq7gx3vqRuG3Ut776BLBM/3ddJf58KfuDiWhP3wMI8D7Bju3t3FD7rnmO04+tC+AWshklxt2nyVeqeqb+pXVfBAw44QHfeeack6b/+67906623SpLe8IY3rPc3y4qi6MlDuoGBAT3pSU/SHXfcod/85jfafffdu/o8+s9a58+fr+2331533nmnbr755nDMG264QdIfUl4BAAAAAACAXut6/FgUhcqyVFmWGh8fX+tnq9vX579eWf0w7cc//nH489UPGZ/4xCdq4cKFeuELXyhJ+u53vxv2X50Wu88++9Q9VQAAAAAAACDU9U26Cy64QJ///OfVarX01re+dU37Jz7xiWmdWFXvfOc7deGFF+qUU07R+9//fs2dO3etn59++umSJr4hKEmHH364/vVf/1VXX321brjhBj3zmc9c0/eKK67Q7bffruHhYR1yyCHT9yIAAAAAAACwUet6SPf0pz9dJ510UlfH97znPdMyoaxXvOIV2muvvXTVVVfpwAMP1Le+9S1tvvnmkqRPfvKT+t73vqcFCxboQx/6kCRpjz320EEHHaTvfOc7Ovzww3XppZfqcY97nO6//34deeSRkqRjjjlG2267bd9eEwAAAAAAADYu8cqcM8y3vvUtvepVr9Ill1yirbfeWrvuuqt+85vf6P7779dmm22mM888U9tss82a/qeccope8IIX6MYbb9T222+v3XbbTbfccouWLVumF7/4xRvstwYBAAAAAAAwO82Kh3Tbb7+9LrvsMh1//PE699xzddNNN2nTTTfV4Ycfro997GN66lOfulb/xz/+8br66qv1wQ9+UOeff75uuOEGLVq0SEceeaSOO+44DQ4OTtvcy8GmykETp7a6TzZ5K0o6cUEpLnElSFbpmAQVEzqkjklniRLXTMiRHTtqd6+ll2OkwltqSHKrRSpByaUw5VKi2p0g8Spok6R2I25vBQdJy6R6jZs3rNHpvuS5BNaMtkkoc4lcg41Wd1sYtSiNm5SzaN7NZIprhk9xrZ5+5tLMMu0jnfj6PBq8t5I0FrRHbZLU6sTzaAXH6li7el8pTt5rm+SojhmjE/W30d3V07Fq06PrYfolRmMk7gOSFB0KhXm/7BjBYWZCxOyZ63MIg5+0XQqbSS6LEmLdGsJm4tHQtYSNp1Pjo8Q2935lahNzjg5UH8P19/WGGSM4JrPHdZxgW72vpFoSYntW97iDL/Eiw+us/HW53YnP3uia7+4P7n4yENQFrSLuO2bO86ZLUMwIpu32k6tlBoN5u76uLmsmk1wzohrHJs8nElvHTb2RqXtcfZNpd3Vxy7Z3z8PVtLaeD9rd5wSbdp+Kk67etTY9uh76a7u7b3S/+J6O4R5dhKnWhrmGuKDSRjSSTZ43zyyCY6QwSbAu7b5K3VO6X06YFQ/pJGnu3Ln6xCc+UflbcFtuuaU+97nP6XOf+1yPZwYAAAAAAABMzv8jLQAAAAAAAIBpwUM6AAAAAAAAoM94SAcAAAAAAAD02axZk26mag811R5az8UFzcKK0aNXu0ikGSMMVRjIBUdkFlA2a59OMkbU141RvT29CHO4r+O+qeCIXoZJmIVVy2ABT7fAq11Y2fRvtbt3yrgJiGiYBZQbNawIGy0CbOfcrL6or18o2bQHYRB2AWXzups1BF5kuAWDO+ZgjfpngyOiRbZHTXCEWxQ5ColY1Y7HcIt3jwT9x03IxLgLlAja3WLL7vyKFrx16/y79ow61hu3g7h7T2YBZXdtD/Zf2THhCW6MaPFju1Pd4seVu6poVV/IWZIaQUhEYYIjCvPao0uLDY5wl5s6DrRIor6RTC2TDMyK6o1MAJYkmZyfcJxMfeP6p0O36hgjU7PUcBFxQ6QysBLXSXsKmOuyu45H1/xx82Iapj269zjuXhyJQi2kScKTgnptoBFfFDL1kKtjZmrdEwdHTD0wy4VdjbarB0eMmXm4Y2y0FYRPmPpm3J4D1WvuzOeK6B4/8YO4uRY1fIar47Nkph5yOSv2Ghe2V6/VpDivyWVUdky4gzvNy+7MPZkcHVvLRPVQaY7fqdQ3bTexBL5JBwAAAAAAAPTZen+TbtWqVbrlllt05513qixL7bfffpozZ06dcwMAAAAAAAA2CumHdL/61a/0mc98RqeddpqWL1++pv22227TokWL1up777336tJLL9WBBx6oTTbZZOqzBQAAAAAAAGah1EO673//+zrssMO0YsUKlY/6O93CrB2y33776ec//7n+/d//XUcfffTUZgoAAAAAAADMUpXXpLv77rt1+OGHa/ny5SqKQocccoj+7d/+bdLfec1rXqOyLHXRRRdNeaIAAAAAAADAbFX5m3Sf/exn9fDDD2t4eFj//d//rZe85CWSpGOPPdb+zq677ipJ+vnPfz7Fac5e7TmNdae72mRW0z/4ZqPvGzdH6Wc+JbV6QqzkUs6q95VMuqtNVYvbw5Qzm6pm2qPXmElEm6R9ylwojdlgnIho9n+Q2CRJJqBQRZBcVrTinRqlKkpxwpNL2IpSQSVpqNGdLjZmUlyjdCxJGgzGGDBRRFHfif7d7S7drWliFaP+bowsl76V6RsllEXpuhN9zfEUnLyZFNeJMbr726S0IM1son+QtuZSzmy6a/c8XCpg25xfnaC/TTlLtheZC1EisdVe98x1Mvpivk9gNe3BKeOPaZP6Gs3DpqGbFLF2sD+CNkkqmi6JzGyzFVxTEymuE/2jvmYemUuLm4fjot8CqbonmbQe1T1RGv1kY9vE1ih5r466x43h5p2Yh2uP9l8qCdaMkU2Cja5ZLs3QXffCVHt7XY6HbpjzvxUlMJvEP/fKo+uWTZONog8V3wMHGnGi51AjHmOw6O4/YOqbQZfYGrTbVPvEsdDLusfdN9qm0I/6t1yKq0vSDU4m19clto4Fqa+urxs7qod83eNquO72lqn93eeKzGcTX3DEzXWIa4W4rwkzjdNds3VPJpHe7r7qNVwQ+DzRP/FaXH1jk+oHugdK1z1Re7buqVDjtN0OSqg8wgUXXKCiKPTud797zQO6ddluu+0kSffcc8/6zQ4AAAAAAADYCFR+SPfrX/9akvTHf/zHlQcfGJh4Cr9y5crktAAAAAAAAICNR+WHdM3f/0nYlltuWXnwW265RZJIdgUAAAAAAAAmUfkh3aJFiyRJP/3pTysP/h//8R+SpB133DE5LQAAAAAAAGDjUfkh3QEHHKCyLPWZz3xGq1atWmf/73//+zrrrLNUFIX+7M/+bEqTBAAAAAAAAGazyumuf/M3f6PPfvazuuWWW3TwwQfr1FNP1ZOf/OSufmNjY/p//+//6ZOf/KSKotD222+v4447rtZJzybt+U21hteR7mq4pLkweauOdNdEksvk7UECTTJdLGp3iWiZxFYT+phKd/WvO5Gyk018DdNqTN9E+plLIuy4iCKTyBNdamyKayeeYJRoNm5SzgYbJm21GSSztuKUs6ivFCe5upQzlzoWjeH6RoloktSsKdGsqnYyhjhKZnVpvC5BLUwGMyejS/SNkstsX5fMGh17JonMjREltqZTzhJpoTbd1anhcIqGyKSZSeZa68K7XLJlNA+3O0yyaJgu5hLLzb6OLgsdkxBpTvNJ2ru36VLO7BiJ+0bhUuKmWR11T6Y2qaO+kVzNEo9hj7Mw3bV6mqwU10Op9HrT3ybpujFM95RU4rB5X6JrqjnGOubFjJtNRqeMu4+2zfzCuqeRS/qMklxt8rypnaJEetfX1jLBG+ZqJ2dDqXtcUn2c7lq9rxTXTm4MW7NECbEumdWlBQf9o5R6yddUUY3Tdp8rzBhhjWPT6+Nme63IHE6Jz9Hu8567P0SXlkx94xTufhSHOIdJ1XYe5n0M654a6hvXv550V9d3/a837dH1e7bzaJUf0j3ucY/TaaedpsMPP1yXXHKJdtppJ+21115rfv6Wt7xF99xzj2677TaNj4+rLEsVRaFTTjlFCxcunPJEAQAAAAAAgNmq8p+7StIhhxyi//7v/9aznvUsjY2N6bLLLlPx+39x+sEPfqCbbrpJY2NjKstS2267rb7xjW/oFa94RU8mDgAAAAAAAMwWqYd0kvSyl71MP/vZz/TVr35V+++/v57ylKeo0WioLEuVZamddtpJ//AP/6Cbb75Zr33ta3sxZwAAAAAAAGBWqfznro91xBFH6IgjjpAkjY+Pa+nSpdpiiy00NDRU2+QAAAAAAACAjcF6P6R7tMHBQW2zzTZ1DLXRGZvXUDkn/YXGCYl1wTOLLUsmECG5CLNdMDgTtuDCIBKLbGYWP7aLLbsxwkWYzUKTbj9lgj6c4Bfc+rp2Hcxgsc/SrPreMUuXluZNKIOFPd1i6+1mPPZ4sJB70yyg3HSLHwdju75ugeJoUWQ3hg+D6G6PFlWWpCIxRj+4xY+jYJBscES0aLZboNge1onQB7dwcbR4dBQEMVl7FAbhtucW5A0XULbBEXFzemHljHABZdPXLaxsrvkRt3Z0dNmyARYuyCHYry5YKFhTfaJ/cN2rY6FkSeH7ZcdILJptF0pO35R6xNUhUY1j+5r26BhJh09Ub8+ETNgxagjdSoW42DHcQunxGOlwrHDwoM2dL2570XluirXSnBsuBKsMzkdX37SbJhggqllcTWADs6rXPW43RSER6boneMNmU30jxTWOGyNT99i+mdrJ1Syu7onm4er2VN2TqG8kldHYru4x+6OwF/3qFyLXNfx4ba+p7riOrkOJ7cVD2LqncNf2YL/amsUEBUYhEa42cS8yE3Zlx64jDGIKdc+4uaZnTH0EAAAAAAAAAFOy5t+Ff/7zn+vaa6/t2YZe//rX92xsAAAAAAAAYCZb85DuvPPO0/HHH9+TjRRFwUM6AAAAAAAAwFhrhRW33sJU9WpcAAAAAAAAYDZY85Du4IMP1hOe8IRJOz/44IM6/vjjtWLFCi1atEiHH364ttlmGw0MTAwzMjKi2267TWeddZYefvhhHX300fqXf/kXbbrppr19FQAAAAAAAMAMtuYh3W677abddtvNdhwbG9Nee+2llStX6m//9m/1qU99So1GnDvxqU99SkceeaS+9KUvaXx8XKeffnrtE58tWvMLae76pYekQkcyiWiufx1JaWacbFJaFHlix8gk2JrEG59U2/0t0XTKbLiva/j2aTblLPpBK+5ZuoTNIM1QiveTS7B0KWdR2ldhEuVsYljQP2qbdIwomdX0dWNE7Sboz3Jj94pLM3OiQ8Gm4CXSz2xQpRsjSH5yqaruS9/RGO4csMllUXumrxQnmtkU12zKWcU2yV9DMknV5jrZic5zc5+yyWVB2qrr6zLuo3Qxn8zq2qMdkhsj1b+HaWYbSNiiP556lHafTSftZd0TJ8RW354dI1uzBAmF2forU2P6i34whE1aNWOEnd19yqUZmvbgpl64hE1XhwQJqr5mCZvD1Nds3RMGWFL3dOlV3ZNOmQ3rntwYUR1i5+GSWcPaKVHfSHEtk0wLte0Z5tgL7yfpWM7oYpa7lkXXYJfi6k6jsK6oI5E+WfdkEunTdU/inJ7K5abVnHqEeeXD6POf/7yuvfZa7b333vr0pz9tH9BJ0ty5c3X22WfrOc95jr761a/qwgsvnPJEAQAAAAAAgNmq8kO6s846S0VR6O1vf3ul/kVR6N3vfrfKstQXv/jF9Z4gAAAAAAAAMNtVfkh32223SZJ23HHHyoM/4xnPkCRde+21uVkBAAAAAAAAG5HKD+lGR0clSTfccEPlwW+55RZJ0v3335+cFgAAAAAAALDxqPyQbvU36M4///zKg3/ve9+TJG2zzTa5WQEAAAAAAAAbEZNr1u1Vr3qVrr/+en3729/Wv/zLv+g973nPpP1PO+00feUrX1FRFNp3332nPNHZamxBoc689UwAqSHd1YkTRxN9k3NJJ7bVksIWJJS57WVS2OpIwc0mXoWJUqavS0VKjGGTldou3TVoc7E5Zl+3E+mubv+FyWC2b9we7WzbNRMNVENqmZ1zDTIBkX6QZLpYLWMHbYlzwI2dPjei/jahLPFGJpO3Mq8xK7qmFubssO9teG03aYEusS2obGwSmX0fp9bX9k+eR6mUMztG8ljI6FXYYh3Xsjpqln7UPZmk1HSqffV7YCqpNl071XDgRBs1Sas29TVzjroES/daEjVmx6atBvGMNprRNEfX5UR946Trno24xom3l7g3Zo5fM4hNqnUhmGE6eTJVNTxHE32lOPg0Xd/EzalDMnudnOIYhbmW2XojUffUkRqf2qc9rHv8e1hD3TNT0l3f9a53aauttpIkve9979PLX/5yXXzxxVqxYsWaPuPj47rqqqv0mte8Rn/9138tSdp88831gQ98YMoTBQAAAAAAAGaryt+k22yzzfStb31L++23n1atWqVLL71Ul156qSRpq6220tDQkO69916Vv3+sX5aliqLQl7/8ZT3pSU/qzewBAAAAAACAWSDzJU295CUv0fXXX68DDjhAZVmu+e/+++/XPffco06ns6Ztp5120ne/+10deOCBvZo7AAAAAAAAMCtU/ibdajvuuKO+853v6KKLLtLXv/51XXnllbrvvvs0OjqqLbfcUs9+9rP1yle+Um9605vUbAZrKQAAAAAAAABYS/oh3Wr77rsvgRA1GF8glfNrHrSGxVlrWEO8lmCLOhZ4rmPtyFSARXZB3jr2dWJhVbvQbBSI4BYdtXOu3t8GRyT2R/Y4DbeYPV96GNiSUsN5XotaLhZu7F71rWHO/VjQP7Mgb/I1Zg4nN3TUbBdVdyuAZ+aRWVzYro5dfXu1LD6dPW2nN3vG6+V5ntHDF1lLvVFD/3SQVqJvfA+c+vbqqHuyh1gYmJUJ3ZFZKD0dZGZqrVRwVNxcy62qhnlMuW8W9c3695V6W+Mk+oY1Th33wGTdU0cAk5UIiHHhHeG1zP2tY+p0TtY9ddQb0xTMsM55ZPTgPG+l/lY1VsMQAAAAAAAAAKaCh3QAAAAAAABAn1X+c9e3v/3tqYHLslS73dbIyIgefvhhnX/++dm5AQAAAAAAABuFyg/pTj75ZBVuDZlJlGW5Xr8HAAAAAAAAbCxSwRGlXaW5W1EUmjdvnrbYYgttuumm6YkBAAAAAAAAG4vKD+n+53/+Z519Op2OFi9erO9973s699xz9YxnPEP/8R//oW233XYqc5zVWgs7KhdE8VLAFKWTCKf3G6/FNG8PQE1qSSKML0QbSKYfuvDOYAMTpTinD9Pe1SHUOMAM1MOkZeqeDVX970C7MfVnO5Uf0r3kJS+p1G/vvffWX/3VX+nss8/WX/7lX+qggw7SlVdeqaGhofWeJAAAAAAAADCb9Szd9bWvfa3++q//Wj/72c/0la98pVebAQAAAAAAAGa8nj2kkyYe1JVlqa9+9au93AwAAAAAAAAwo/X0Id3qP3H95S9/2cvNAAAAAAAAADNaTx/SXXLJJZKkVqvVy80AAAAAAAAAM1rl4Iisb37zm/rIRz6ioii0884792ozM15nQUvFwpofYtYRQVNH0ElmjF7O2aacJraZGdv2jbeXSiNLzMOOm5hfegwnml+ib2ZcqaZ9WkP/fNJcQj1RU1PXwxeZeok19K1lezVcVjJj2zln35fU2HWM4U5e03/K8+jh9jJqeF9q6ZvV04tZQi+vez29b0xzXZape+q4z9dR96RruMQYTh11Tw2Xljreg1rqNeqbtc3E+ibZP32fT/QN914Nr6WuuiIcp47aro751TKP5MVpumucWrZXxxj1n+cdTf3ZTuWHdP/6r/866c/LstTDDz+s3/72t7rkkkv061//WmVZqigKvfOd75zyRAEAAAAAAIDZqvJDuve+970qiuqPK8ty4qnka17zGr3mNa/JzwwAAAAAAADYSKTWpCvLsvJ/ixYt0kknnaRvfOMbvZo7AAAAAAAAMCtU/ibdl7/85XX2GRoa0mabbaaddtpJixYtmsq8AAAAAAAAgI1G5Yd0Rx55ZC/nAQAAAAAAAGy0epbuimqGFoyruXBsvX637GF6ahn0t5szY9v5Be3R9ibdaDS/Tm6MsL/bXmZsF6bTcfOI+uberzCYxozh5pFJiM2MYfu71LfM/LLzqCENLjNGT1NmezWGYw/JxLGaTVWtI00r6usWeqgjtdSMHYYqNswblhijsH3NBN17EM0l+b6E+yT7GlPzSCSXuXnY97G7v12W14xdRPNzb0vyNaaC9+x+mvoFI7FUcUqqvrGDuNoksb5yeuzqA9l5mHt3OHYv655EXWZrlk48eDSPdOJ7Zh6ZU8DUFfZ0cf2j15gdo44UXFcnVd2ek72szMgaJ3mB61Gqajrstoa09rBOyvSVwoMhW39l5uFuSG6bRXDvzu7rsGbJ1DdSLt01M0ay7on2RzYxd7rrHlvfOIn+U6lv2uX4+v/y71V+SPfRj35UkvTOd75Tm266aaXfuf/++3Xttdfq8Y9/vJ71rGet3wwBAAAAAACAWa7yQ7oPfehDKopCRxxxROWHdHfffbde8YpX6NnPfrZ+9rOfrfckAQAAAAAAgNksle6aNTAw8Qxw8eLFvdwMAAAAAAAAMKP19CHdueeeO7GRRk83AwAAAAAAAMxo9s9dzz33XN14441d7SeddJK22GILO2BZlnrkkUd09dVX64orrlBRFHr2s59dz2xnoYXzRzS4YP0eYmYWPzZr99oxovaO6+u22YlfV7T4sevbcQsoB+12fokxynZuIeew3YxR2oUwg3nUsCBveq3wKMCibcZwoRS2f7W2dHtygeeejZFdyDl6g3u39mkt6lm4OBdm0KvFj+tYuNiNkWpv5BY5LpuJRZibZozUgsvJgyyzgHLwWibmEbS7vm4x6Gb3SRouiCwfvBEtRlw04otFw4wdtbtFjt0Yrn8zmp+7xyQuFuaQzC/O3COuZnE1TmaMqL2dDJ+wNUtUU5m+vu7pPljta7FBVVHdY04Ct1ODGsddV6L65ve/UWncSUUbzd7Pg5rF1iBmfjaYoUd1Tx3BXT0N3aLuecwYNdQ9iTApSeG9sadj2PbgWmFrAjePqCZwfU27/Qxs+kd9M/vJ1jdmjKi/Dd0y9/lgDF/3mDokrHtyNUsjqJPcKRD1lSbJzKih7nE1TmaMKsbbI+v9u6vZh3R33323PvzhD6/VVpalTjrppMqDl2WpRqOhD3zgA+s9QQAAAAAAAGC2s1/hOvjgg1WW5Zr/Vnt027r+23333XX++efr5S9/+bS8GAAAAAAAAGAmst+kW7RokS6//HK1222VZamXvexlKopCX/va17T99tvbAYui0FZbbaXtttuucgosAAAAAAAAsDGzD+kk6UUvelFX21577aVFixb1bEIAAAAAAADAxmbSh3SP9oY3vEFFUWjhwoW9nA8AAAAAAACw0an8kO7000/v5Tw2WlvOXaGheeuXHuLSTCM2icxkqLQSKWJRX8mnorWD/m2TZtZqxxE+0Wtvm4Qym6AW9O80zRgusbWVSOZ1Y0QJgC55y2beBBKJY1KcctZoJVNcM+mupm/DtEfzzvSd6J/Y14mUuF4mpYWJaJP0n3YuETWKczJJSe5S1rNkVpd8mkhbtWMk2jvphLJgfm4Md4C4u37imHRh13EqXS6hLEo5i1LLJtpNMljQ36WINc3YzWBsN8aAGzt4jU3b17TbhLLqY0d9JakRvME+EW3DuODYFPcoPdVcFNwYUW1i+ybGcO227rFp98EYpu5pm3ojGsOdoy71NXzp5l7sk+qrJ8S6S0V4L467TlKzdP9GHfWNFNcn6bEz6bMmjTczRh0JsfZSEb1fs6m+Mf1LGy1sxo7qjUSyqBuj03QnWNwc1ieZvorrHpvu6hKpg3l3zAGSPmwSH+Fs/RpdoDIprqa9GDD3czNGVJ80TI3kklnDuse87oFmfNFK1T1ubJf6GvR3faP6xo3hTKXuGRtbud6/u2b7Ux4BAAAAAAAAwJR0/Zv6zTffrC984Qtqt9s65phjtPPOO0uSzjjjjClt6PWvf/2Ufh8AAAAAAACYrboe0h1wwAG68847JUn/9V//pVtvvVXSH9akWx9FUfCQDgAAAAAAADC6HtIVRaHy92sEjI+Pr/Wz0q0dAAAAAAAAAGC9dT2ku+CCC/T5z39erVZLb33rW9e0f+ITn5jWiQEAAAAAAAAbi66HdE9/+tN10kkndXV8z3veMy0T2thsPWeF5sx1EZWTc+linSA6xyWUZdJdWya+p2Wienzqa3f7uElxHTdjRP1bJonMjd0OEmhaLZfCZNLWgrbSxi3FzWFsU/Yvy4MxbMqZTSKbespZo1W9f6bvRHuQ2NjDpLQoCdb190lpiTESiWgT/af3m802zcz1TyRe+XTXIAEwkWYmxclgrm/HpGaFCWUmzawzYNIno/4madUmxAbTywZBuwS19DWn6hg2lc4ltlZPOYuSyKQ40WzAjeESW4PksqEBk2ZmTt7BaAxz0bLzcO3Bxcz3rZ762jB93WucbnXUPZmapWVORlvfmPaxoH6ytYmZ91grqHsaZgyTVB/VOB1T35iAWEWZc/Z25KKg67jeRNu0yaemvZbaJG6PxvF9E8msboxE+mxUT7m+dgyXJptJd7U3MDOPDaXuSSTV+0T6RJq8rVmqj+FSQV29UYR1T7y9qK8kdYIax9ZOcXOoYfZdxx18rj4Mr7W5NN6wPZFeL8U1TsOlu7r7fNDf1UhRbSJJA0F/13fQzCPq3zQnerZmGQguctl01+mqe0ZGRtb7d1cj3RUAAAAAAADoM/Pv+N1Wp7seeuihmj9/fqXfWbFihe655x4tXLhQ22yzzfrNEAAAAAAAAJjlKn+T7g1veIOOOuooLVmypPLgN954o57+9KfrVa961frMDQAAAAAAANgo9PTPXefPn6+yLHXrrbf2cjMAAAAAAADAjFb5z13Xxw9+8ANJ0vj4eC83M6NtNbRc84Yn3z829MG0t4Nnr34ME7YQtLvgiNFoZVD5xZnHgsWSx5rxGG5h5dFGd/9Rs4By0yzgOeZWOg20zELCChZR7ZjFdO3isdHCyslFlcO1Ld3ap2511nBx4bhrNrChERzimYWSJ/pXa5sYwywYGi7CXL2vFC+WbBeUNsdCuK/tYsuJxZnd9rKCwAYnE/pgwwzM9spgkV3f143d3WYuZSoSoQ/msmff8zBQovrh8fvBgyaXO2T2k1t4u471uMvooHQLKJvrYREsXNwwY0QBEZI0ONi9U9zixy4MIlqMeHggvuAMmovFUBgcEY/hAiWihZIlaThod/NwCyhHYzfN0Rcttuxk+kq+Psn0jeoeV7O4wKzxoP+oC8wy7VFAxER79wUjqoXcPKR4EeuWOQeikAlJKoK6x1Wg8aLqUifYZukuqubcjQJsimzhEwUR2LrHXA8TYVKZ+kZyNcvU655s+ERY97jQh17WPWFwRDJ8ItLDusdemlzoQxgiEA+RqWVcXx92FXw2cTWSC4MIapzOQPXzeeIHwTxMV1d2hgF9rn4wg7v6JjrO7NFUQ3BEYWuZIMwgERAhSYNBLWPrG1MPDQftw01T97iaKrhoRbWQ5AMihk2dFNUsgzYYq3otk617qtQ4K4dG19lnXexDuv/5n//RXXfd1dV+7rnn6vGPf7wdsCxLPfLII7r66qv1zW9+U0VRaNddd53yRAEAAAAAAIDZyj6k++EPf6h//Md/XKutLEu9973vrTx4WZYqikLveMc71n+GAAAAAAAAwCxn16R75StfqbIs1/y32qPb1vXftttuq5NOOkmve93rpuXFAAAAAAAAADOR/SbdzjvvrDPOOEPt9sTf/x511FEqikKf+MQnJv1z16IotNVWW2m77bbTM57xDA0NDdU/awAAAAAAAGAWmTQ44i//8i/X/P9HHXWUJOmwww7TokWLejopAAAAAAAAYGNSOd31pS99qYqi0Ny5c3s5n43OZgMrtMAkxq1L28TbRMll4yb60I0R9XcprnOTqa9RCutI2ySURfFCkoogWSVqmxjDpC0mUn1cylkZDeLSDE0aVM/YF2OSnBKJo5kkMtff9jXtmaS0VMqZS3dNpLD5pDTTHiWlJZJgJdUTx5nhjl+XXBb0t2lmbowgdcwlkdmxgzQyl+LqksHCS5lNzI7HiK5PjeQ+DYOgXWJuIsVZkqKwa59y5l5kxTb5lLMieO0uxbVp2qNkVpdy5lLHoiRXl3I2pxnHO0aJrW4Ml/rqUs6i9kFzYc6kvjZdIpo5Glz/qcrUN65/JsVVytU9rj1KcZWk0Xb3ezBmEulH2oNhe5QoN9qqXMZLimuZjrumuron2E/ufI7uAxO/ELW5K465XkeN2eteDan2Pu0+uPck6puJ9uC+4fomahk3j0wt4+ub6mO7tPEZW/dECbE2xbV6LRMl3Uu5RProHi/5pPqwxnHBva49GMKmuNpaIWizKa6uIJru48a0u7onShw17/lAIrF1yDxfsHVPUJ/McTVLoh6K0ugn2qunuLr+ru5xtUnUvxf1zdBAD9NdH+vSSy+d8sYAAAAAAAAAdLPBEQAAAAAAAACmR+578r+3YsUKXX755br77rt13333qdPp/jpgWZZqt9tatWqVlixZoi9/+ctTnSsAAAAAAAAwK6Ue0i1fvlzHH3+8Tj31VD3yyCOpDfGQDgAAAAAAAIhVfkjXarX08pe/XFdffbXKxAKMCxYs0DbbbLNekwMAAAAAAAA2BpUf0n3xi1/Uj3/8Y0nSM5/5TL35zW/WggULdPTRR6soCn3/+9/X4OCg2u22fv7zn+vjH/+4Go2GLrroIu288849ewEz3SbNVVpoUlZWayeXDowSyjqJFFfXPljEiWPjJqGsYRJbG0Fia9OkekVpZo5PIjNpcEEyYKcT76d2kBYoSUUQUxS1STYMMk40qyEJ1gYwuvYo5SyRiDZZe5Q6lk5QC1POXJpZ9TEKmwRrUpiiMWpId3WxoD79LEpKi7vWwR6/NrG1esqZXMpZO5FEOBCfu6lQLx9nGjQlU5yD6bXdpd0m25n2aBp1HAt1jJG+EAXnuUlEs+muQXvTXMMHzf03SkR1Ka4u/SxKIptrxhhuuPY60l3jMaL+TfOmN8zFvWkP1qlxdY+rZdqJVPtxkxof1z1mn5r20UQ6bqa+kaROIpm1bS7Y0bnh+tp6KKxZkue5jT8M1HAdSu1q1zeZHBvVMj4htnot49PrTft49fqraJnjN6rhXJR5Lan2rh6Ku9chPA1sPW8Sh3uV7mpSXN2xVyYS6XOFhblWuLInuHd3zP7IfK6wfZOXodThlPoMZ7raVPvuF9QwNYv7vDwYXNsHzRhR8rwU1zJzBuLaxNUyUVK9q2PcGO7+Oieok7L36KiWcX2dKnVPozmWGjNS+SHdWWedJUnafffdddlll2nevHlaunTpmp8/97nP1aabbipJetnLXqYDDjhAe+21lw444ABdd9112mSTTaY8WQAAAAAAAGA2qvwVrVtuuUVFUehtb3ub5s2bJ0naZJNN1GhMDPHAAw+s1X/HHXfUcccdp1/96lc66aSTapwyAAAAAAAAMLtUfki3YsUKSdJznvOcNW2Dg4N68pOfLEm6+uqru35n7733liSdc845U5kjAAAAAAAAMKtVfki3aNEiSX94WLfa7rvvLkm66KKLun5ncHBiDbM77rhjvScIAAAAAAAAzHaV16TbcccdddNNN+n222/XH//xH69pP+CAA3Tuuefqq1/9qt70pjfpBS94wZqfnX322ZKkrbfeusYpzy5zG2Oa96iFJDuJkIi2WUA5WgDRLaDcKM0isaY90jR9o0WOHbcY9IBZHXQgWPSyaQIs3PyihcTdQp2r/6z7saI1VG3mQ3ah+V6pYwHlZHsulMItOlzDGMGCxpmAiIn+wUbdGLa9ewy7IHLHvMjEeu02fMKIFkUubMiBWQQ4GMMulOxWHW52b9RnLbgd0v0b7pSz7dGC10GohTTZfqo27kR79YWV7YLIvVxA2clcyty+jtrMpO2lNujvFlt295hoYeUBFz5hFmEeDhZhriMgQqpnAeXBontsF2YwFPR1GskwiUzdMxauiB7XGy4Eo+mCSBJ1T1YnOFqjNklqmXooOibbJtzBHe/RuZE5j2z/bBnTo7LH3kcT18N8YJYLP0iM4YIcwrAF09ecotHYmYCIif7RPdAEqmWCtNyNygZmBePWUN9I5pA0l6bCfCYIP/YEdYzvLBVBneSukZmyzNUmqXqjhpqlrs8VU+5bl7BocddfM0TQ7gKzbKBE0O4Cs1x79Nnafg437VHNkqljJF+zRPVTvu6JArNcSEf1muCxdU/LvLaMypXRy172MpVlqfPOO2+t9sMPP1xbb721Op2ODj74YH3kIx/RaaedpqOOOkonnHCCiqJY609kAQAAAAAAAKyt8kO6o48+WgsWLNAll1yiBx98cE373Llz9elPf1plWerBBx/U8ccfrze/+c0644wzVJalms2m3v/+9/dk8gAAAAAAAMBsUPkh3cKFC3XJJZfo3nvv1RZbbLHWz/7iL/5CX/jCFzRnzhyVZbnmv0033VRf+tKXtMcee9Q+cQAAAAAAAGC2qLwmnSQ997nPtT/767/+a732ta/Vd7/7Xf3mN7/RokWLtPfee7MeHQAAAAAAALAOqYd067LZZpvpL//yL+scEgAAAAAAAJj1an1Ih7xmUa6VwtVUkN7lUlJdOlYQCONSxFxzlLbqklJcyqxLiYvaGyZBxbdHY0x/rI9LP5t2dUyjl6mvNcwjTGFLpMnasV1aVSZdLDlGOO9kims4djLlzAkvCyZ9zgRHx9z74iKvwtdoUsRsSFyQSpdMPi0SSdV1HL+zSjLFMXNNzd1jqqdd2u2ZNzfTnpmz5NPFonuju19GKa5umy7F1SW2uhTRjEzd4+Y3FpSzg4r7ts1B2Qja/f7PvY/hsZA8niLu+K3j3MhuM+6cGnpGSt1PkqdLfN9I1iyJe6Cv7aJ7senr0m4TNUuYBOsk6x67n6JEepc8b66HYepropaU4vTZwiR92kDq4PppU4F7ePyGNqQU1w2kBpvuz5Ibyudzd69zzxuixHbf19QsQbu9z0+h7qmjNqqeew8AAAAAAACgJ9b6p8f58+f3ZCNFUWj58uU9GRsAAAAAAACY6dZ6SLdq1SoVRRH+edBUFO5PmQAAAAAAAACs/ZDuyU9+Mg/UAAAAAAAAgGm21kO6xYsX92seAAAAAAAAwEaLdNc+a5eFT2/9vY7J93Cpqu1EHkiU4joxRvec3PZscqzdZnd/Nw/fvmF847PcQOZRi8xL6eXLdkGfQbvb/eawiUNEG/EgLjUrTNkyCWD28AhSf+y3mG2qV5DeZTaXFaWLpUVjmH3d2+MpSGwzr8+3Z7YXN2eO31mlh6FlddwHMtdwd6/LtLs5u3Z7343aXQK7qwmCi5zta9PkXbzg1GTrnmj/ZWqhiTG6+2e2N2l7dCwkj6dIHTVIT+upDSQ9sZdsHRK1J3d1fN/I3afi7jW8Me61uPt8mBZq6ptEDdLTusddQhLvgX0ttuZLjJGoqVxd3MvjN2VD+QzSBxvKZ8np/nzu7nXuvhslyrp5uMTWsC5wl8Mp1D3rerZTBemuAAAAAAAAQJ/xkA4AAAAAAADoMx7SAQAAAAAAAH3GQzoAAAAAAACgzwiO6LPRclAj5dCkfeyC0uYZa7j4sVn9cLxsxu2d7kNjtIwPl5HOYNg+GowhSa1gmy2z8KNt7wRjdMxC04l2t3inb6/WNqk6Flbu1YKrNvjAtLuFaaN1Ou0itmZfBwEKdgwXBhH0Dw4lSVKjE78x0YLmDZlFRON13KUimog5EEx7NESZPviMxKLNZTMRwmDfF9MejO221xmo/p6by17qePJ94/bwnMmeX5nwiX6sQ5w5/Ezf6Fprr79m6Oja7u4Pg+Z4iu49bozxIj6gBoKFiwdM34Y5d127u+SEXaOLheJFmF1d0UwERDSCYJzJZBaadvOLahm3oLSre6JaxtZI2fbgRhO1Sf44i45Jt/B2ph5y71aqHupH3ZNRx+L4dVyvk/eN+N5jwhbcfS04dd390n9zI/hJ27yJ5vwvo8Asc4wVbuxo3B7WPfbSlKhlXH0jV8s0g3PUvDEdV38l5pE5Jn19Xn0MH8aRbJ9q3/XpH0kcfu5QjT9L5j6PRtd29/nX3mOC+qTViO/9/nN79xij7vON446RaCqmr3vu0Qx2dlQLTfSN26vUOKM1XJb4Jh0AAAAAAADQZzykAwAAAAAAAPps2h7SjY2NTdemAAAAAAAAgBml8pp0F198sX7wgx/oyiuv1D333KOlS5dqk0020ZOf/GQ96UlP0tOe9jS9/vWv1w477BD+/oc+9CF95zvf0V/91V/puOOOq2v+AAAAAAAAwIy3zod0X/va1/Sxj31Mt95661rtZVnq4Ycf1l133bWm7SMf+Yj2339/vfvd79aLX/zitfr/5Cc/0Y033qgrr7ySh3QAAAAAAADAo9iHdCtXrtTrX/96nX/++ZK603M233xzjY+Pa/ny5Wva2u22LrjgAl144YX69Kc/rXe84x1rfnbNNdeoKAodeOCBNb+EmW1Fe46K9h/SXV0aScQll7WDdpc45pLBotQxl9bqxnb9x4L2Ve04IXasHY89lkhKG0+knLnEm04nEQ1mo9JqGCMhszkpmcBaQ3smiUzKpR917A7s7u+Ct1wiYpQ6VjZMX/MaiyDpyKaZmTioZIDilPkUUfOD8HiaerqrSzNzizd0BquPYRNig/fRpsy64zoxRjo5NuqbPP9DtSSf5SYSndPu+tux1/bui4i7to+24vtUI7iGZFNLIy7JtNXMpYUON1qV+w4Wccx0M0g0c30dl4o2Va6+caLXHtVCrq9rz9QxkjTaNu1B/xHT1409HtRD7vj16a7R+ZWre+J01zouOEZiiOw0ouZ0feOu48F93oUcFnZfB2NkL8xR0Pp43LVt7sWNoD5xNUtU37h2Gxo9MDPrnjCRPpuqmki1z6Tdm8vKJGNXa5u0vUefK2zfHl6GUp/t7ESqj+Guv+7a3gzOL/f5t2kK1ajGabSnfo9vu89qrh4y8x5vdM97sKxew0lxjRPVQlKuvmk+Znsr2lNfUS4cYXx8XK961at0/vnnqyxLlWWppz/96Tr55JP105/+VL/73e+0dOlSPfLII1qyZImuuuoqfe5zn9PTnvY0lWWpdrutv/3bv9U//dM/SZJuuukmPfTQQyqKQvvvv/+UJw0AAAAAAADMJuFDuuOPP14XX3zxmodzF1xwgW688Ua95S1v0W677aZNN910Td+tttpKz3ve83TMMcfopptu0tlnn60/+qM/UlmW+vjHP65rrrlGl19+uSRp99131+Mf//jpeWUAAAAAAADADNH1kO6uu+7Spz71KRVFoec///n64Q9/WPnbb0VR6LDDDtMPf/hD7b777irLUscdd5wuv/xyFUWhV7/61bW/AAAAAAAAAGCm6/pD3rPOOkujo6PadNNN9a1vfWutb81VtcUWW+jss8/Ws5/9bF1++eVq/H69psMPP3zqMwYAAAAAAABmma5v0l1wwQUqikJveMMbtP3226/3wE95ylN0xBFHqCxLdTod7bnnntphhx2mMlcAAAAAAABgVur6Jt3ixYslSfvtt9+UBz/ooIP0xS9+URLfonMeHJ+vkfHhSfu0XWyhEaWl2HRXl6wS9HfpqaOmvWXao2RWm5RmksvCMVwSrBujFaTBmTQW1x6lopUuCbY3IXieTaVy8WJBolQiCVbyCU/xcZaL6SoT8zNhq6mU2aIdDx4ns5oxbMpZNLlkmlk0jx4mn9l0LJdcFgVeuShd+35FKWdxX5fYGr/nLhGt+tg2Kc2NEaWtZZPSorQ1d3tw7cnQsbizGySKYDZdXZphcG64vu66PO5ilafIpSq2zAUnugeOmRRXlyIapbhK0qpG90Vn0CSRDQR9JakZ3JRcgq1rj8aog6t7XM0StbsxXG0yHpxMrq9PfTV1SJDk6vraWiuoZXzdY8YO+rvzqI5z1KcZRm25CMZoCBc27u8xQVsyed4ltkavxyUOukzlqO5pNKrXcFL8Gt39MkpxlaQyqIdcfeMuCWHtlKhvJu1fg/DwS9Q3kuI6OlHfSO4+X72vFL+/2dopusT5+sa1V59HKrE1/fnG9O8Vcw6462T0WbLt0pPNaxyrNrNJhfdRk7TaKuPaJEpmHWrEB0hUx0jSsGkfDNoHzAcwV7NE6a6ur1Ol7lk5PpoaM9K1J5csWSJJesITnjDlwZ/0pCet+f9JdQUAAAAAAABiXQ/pNt98c0nS7373uykPvnz5cknSnDlztNNOO015PAAAAAAAAGA26npIt/obdNdff/2UB7/uuuskSdtss82UxwIAAAAAAABmq66HdC9/+ctVlqVOP/30KQ9+2mmnqSgK7bHHHlMeCwAAAAAAAJitulbyO+yww/TJT35SP/3pT3XSSSfpne9853oNfPLJJ+vqq69WURQ6+OCDpzzR2eqB8QWaOzZn0j52oeTESphuIWIfHNG98KPr60IfosUjpXgRSruAslkUORp73C0ibsaIFktumcWWbXBEsJhuafoqs7BycnHc6K1xCyinQh9y6wJb0SY7bpFjGwbRPRm7EK4JfYjWIu3Y0AfXXn0B5aivZNa2devdmkCJzDHiXotjwwjCwc0YiQWU7RhRcERiAXDJLKCcDD6J2m1whFuEOVoMOrkIc2qRc7O4ePY9yCiCk9cttu4WBo+un+Z0TinNhcW1R/epcXNPG2zEL3KwGYQ7BAECkjRgxhgywRFRf7f48YC5ADSC9qYLjkhccLKLMLvaIuxrDtR2MEbHHOxRfePm4eqYbN0T1WCuNnHHWSvo36qh7mm3coFZYY1jz3N3D6zjgtPd5O8P5piM7g+u7jHXZXdJDXeJqXsKN7/omm8uiC6IKKpDXEBEx71fwYuxwRGulonGmEX1jeufDmCroe6J6pA6gqpsQJeth4IxsnVPVDslwyf8+2j6h50T7S50x32WDE+OXGhkuD3T7k7dsO5pVA+BlOKQiKapb7J1T1Rb+PrGBEeYUIpwDLMHq9Q4q8YHK2/Hb/8x9thjD7361a9WWZZ697vfrS984QvpQU899VS9853vVFEUWrRokf78z/98yhMFAAAAAAAAZqvwMe2nP/1pbbvttmq1Wnr729+uPffcU+ecc446ncn/ueKCCy7QXnvtpbe85S1qt9tqNBr67Gc/q4EB83gcAAAAAAAAQPefu0oT4RHf+9739PKXv1wPPPCAfvrTn+rP//zP9YQnPEHPe97ztGjRIu2www4aHh7WHXfcodtvv10/+9nPdNttt0mSyrJUURT69Kc/rT/7sz+b1hcEAAAAAAAAzDT2K27PetazdP311+vII4/UxRdfLEm6++67dffdd9vByt+vLbDFFlvoy1/+sg488MCapwsAAAAAAADMPpOuSrjNNtvo+9//vi699FIdfPDBmjt3rsqytP899alP1cc//nEtXryYB3QAAAAAAABARZUWi3vpS1+ql770pWq1WvrJT36i22+/XQ888IBWrlypLbbYQltvvbX23HNPPfGJT+z1fGed+0cXaHhk7pr/nUk5c6Ix3Lgu2S6T7hqlwUhx2prr78ZwyWXR2C6JzCezBq/RzMOlhYYpW8mUszrSXasOK8k+mg/Tz0ySk2NXrYyGMelMJkgo3NeFTWZ1yWVBX/Pe2sSwYOhUiqvM/Owblhs75PrWkZrpxghTzqonokkmzdQdv24eURChTQYz53mULlZDQmw2oayXKWd1HAupa5lNP4tazbXdTSO6B5rtufbxIP2waS5OA0GK60T/6ilnLlXVpZ9F47jEMZdQVgT97RjJxNZesWn3UbKw62sO9mgMW5vUUveY+sslswZjuOM3V/eYc9GkvsZ1T6K+maw9I9pkMkkzqkNc4KjjXnoUzhhf33wdEiVeZuueKMzQHXu+Zonacu95mORqzpdaLjc11D3pj2Rh4rA5vzLHajYhNlOzZOoeU2+4mipVf2Vqp2R6fS13rxrSXX2qfVCz2A9r5v0KttluxvVDu2lSwYP+tjYxY0e1Sbbucf0z6a5RfePGqKO+eewYo6NTz2NIjTAwMKC99tpLe+2115Q3DAAAAAAAAGDCpH/uOlWXXXaZXv3qV+uwww7r5WYAAAAAAACAGW3q38WbxN13363zzjtPhfsTJwAAAAAAAAC9/SYdAAAAAAAAgHXjIR0AAAAAAADQZz39c1es2wMr5mtwcN56/W4mCdamnJlAk6i/G8OlmdltBgk0LlXVpWxFY5QumdUll0XtNnHQJfUE7a6vTc0KYzDjzk6UTGMSm5wwQcmm45j9ZNOjuttcempprkphfzOGTVuN+ifSzCSTOmb7uhOs+ntTS4prHZKHZPgS3amRSX3NppxF/RPHqWvP9JUUp76l01273+B0yqx9D6r3dcdkdLi7czFzibOnkXnxUfKbS+hum8TWRpAuVphEuUaQ4iqZ9FQ7hkkoC1vjNDKXZuZTzqr3dXqV+ppNug8TfRP1jR/D9DXzsAnxUc2SqJGkuGaxtZO7v4a1k9nXmbrH3nOr10P2upJJpHZ9XfpkMBGf1ppLgo/SXd2nLjtGmKrq+rpk4aBvsu7JJHdn6p705WMDqXFSl6ds3RPVG9nU4jBlNjlGpu5JzM/VJq4u6wR1j/2KUbY9cwC6JOLgnHGfnWQSs6MipzQ3sI5rD65PhUn5bps6pIhS411fd/imxjDprolaJlPfTNY/HmP9LzjjK6b+PTi+SQcAAAAAAAD0GQ/pAAAAAAAAgD7jIR0AAAAAAADQZzykAwAAAAAAAPqsawnTt7/97bUN/stf/rK2sWarh1fMVbO5fsERbm3W3CDVF1BOb86OHbQlFhd2Y6fHCAMbEn1d/0xAhGRDB1Iy2RPupQSLuYYLH0vqmEVAMyEMmQWlXXtqkWOZNWKTB3ZmjCK1On5uHj1dQDmjlgWUN5AxMosw2xX9zdjhQs7uwDFjRGMnF3K271cdAQDhOW0WADYTCe9rdmH7eOxwsWr3+sz+awcXP7dQsjKLItvFjM3YbgHlqtubTKK/nd80S9U9NYRP2M0l6hs3kAuOsAFW4b0nWVdE/TN9JXMvTs7D9c9IBGbZ9yW6LtsgiGzdkwhKyARV1RH6YNj5pcZI1nZT7dtL2cO0jvCJGrLkavlMEN0b02NUD7vKzCP9WrL1UEb0+cbVN/YaEqY4xV1N3ROHhZiQCRuOk0j0SYyRrW+cVN1TQ007lbqnvWLq2+96SHfyySer2FCqMQAAAAAAAGAjEIaBl7V8RQsAAAAAAABAFV0P6f7nf/6nD9MAAAAAAAAANl5dD+le8pKX9GMeAAAAAAAAwEaLdFcAAAAAAACgz8I16WaSq6++Wl/5yld03XXX6ZZbbpEk7brrrvqLv/gLHX300eHvjI+P64QTTtAZZ5yhxYsXa+utt9ahhx6qD3/4w9pss82mcfbS2LIhNTXc+w3VkgRbwxjSesQURWNMvW+YQFVHwmYikXZiHokxEkO7JDK7+xP7I0qCnax/mAaXfpGJ7RmpsJ9pTiKrI1xzQ1fHqb+hpK1lt5d6e+tIfUsnfSb6JsK+fOq2SS1tRxczMw8XchXup6nvj+zxm3vPe/h+9XKMDVk/6p5aLnJu7Iptk0jVPTUkffay7smkGTpRwq5Ltc8k97of2KOjlpp26mPU1j9AjVNRDfVGpI6E2Ow9I3zLa6iz0nVPNvU1GsKlJwd1iLtW+DEqtk32g0zdk6lTa3nPazr5Z1Dd0142PuUx1vshnXsA1rWBgQGdfPLJajTq/9LeL37xCz3vec+TJM2dO1fPfOYzdffdd+uSSy7RJZdcoosuukhnnXXWWr8zMjKifffdV1dccYWazaZ22203/fKXv9SJJ56oCy+8UD/84Q+15ZZb1j5XAAAAAAAAwOl6SHfJJZfoAx/4gA4//HC9613vsr94+umnq3D/iv0YRxxxhF784hev9ySdXXbZRXvvvbfe9a53af/999fg4KAk6eSTT9bb3vY2ffOb39Thhx+uQw45ZM3vHHvssbriiiu044476rLLLtN2222npUuX6kUvepFuueUWvfWtb9U555xT+1wBAAAAAAAAp+vrbd/73vf0ox/9qPKDqrIsJ/1Pki677LJ6Z/0ol156qV75yleueUAnScccc4z22WcfSdJ3vvOdNe133HGHTjnlFEnSN77xDW233XaSpC233FJf//rXJUnnnnuubrrppp7NFwAAAAAAAHisrod0V199tYqi0Ctf+cpKA/zv//6v7rzzzvC/Qw45RGVZ6qc//WntE1+X5zznOZKkX/3qV2vazjzzTLVaLe26667ac8891+q/xx57aJdddlFZlvr2t789rXMFAAAAAADAxq3rz13vuOMOSdLTnva0SgM84QlP0JOe9KTwZ69+9at13nnn6Ze//OUUprh+br31VknSk5/85DVtF110kSTpwAMPDH9nn3320S9+8QtdddVVvZ8gAAAAAAAA8HtdD+l+97vfSZK22WabKQ/+xCc+UZJ03333TXmsjJ///Oe6+OKLJUlHHXXUmvbrr79ekvTMZz4z/L2nP/3pkqR77723txN8lOayppo2KhN16GXY2gajhiTCKJOnL2FcG8P7hdlrI0iwy6U4O9N7ohdmextD4iA2PrOq7qkhmdGj7gGmbDbdR2tIPqbugZZP/dlO15+7jo2NSZI222yzKQ++1VZbSZIeeeSRKY9VxT333KMvfvGL2nfffbVixQr9/d//vV7ykpdIkpYvX66HHnporXk91uabb75mHAAAAAAAAGC6dH2Tbuutt9a999675ht1zmc/+1lJ/oGXJI2MjExsZKBrM7U68cQT10qiXbhwoc466ywdfvjha9qWLVu25v9f10O6JUuWqNPpqNHoeoY5qYULF66zz+owDQAAAAAAAMxOixYtUlFM/k3KRz+rkiZ5SLd48WK94AUvsAO97W1vW+eEVn8jbYsttlhn36l44hOfqH322Uf333+/br/9di1btkzvec97tHz5ch199NGStFb6a7MZfwWx3W5Lkjqdzjp3JAAAAAAAAFCXrq+KPfe5z1VZlrrkkkumPPhll10m6Q9r0/XKoYceqosvvljXXXedlixZov/v//v/9Nvf/lZvetOb9MlPflLSxLfnhoeHJUlLly4Nx1n97cEFCxas10O6ZcuWrfO/xYsXr+erBAAAAAAAwEywePHidT4jeqyub9Ltu++++tKXvqRzzz1XJ5xwwnqvTTc+Pq6zzz5bRVHoxS9+8XqNsT7mzZund7zjHZo7d67e9KY36cMf/rDe9KY3acstt9S2226rO++8Uw888ED4u6sf0vX6oeKjDTzS0EA792e16yX7zDHRP714b7BiZX6Mim2TjZ0Zo6fzqL7qaGY/le6wciuGJl6Ln18NYzuZlU5T73nyT9Az867jC7kzdYXXOlb1jl56elHfzEmTmUfy9SVeS+HGzuwP255YSDh7aqTmZ+bRqT6GfQem/bUkxnAy74sdI9G3l2P0Qy3X2updU/fz9BjV23taOyXGqKPeqOO1pN+XmVj3ZI+bquNO2n+K26trHhuKXtU3dfSV4vlla4Wq407av/r2+lL3RPWGYQ/VRH2YG8PMwwyR29dTH8OJa6fccdPTumeaLjnFiqk/2+kaYb/99tPmm2+uRx55RB/5yEfWe+DPfe5z+tWvfiVJOuCAA9Z/huvpta99rZrNpkZGRvTTn/5UkrT99ttLkm6++ebwd2644QZJf0h5BQAAAAAAAKZD10O6hQsX6u/+7u9UlqVOPPFEfe1rX0sP+oMf/EDvfe97VRSF9txzz2n9Jt1q8+bNWxNYsXqtuRe+8IWSpO9+97vh71x88cWSpH322WcaZggAAAAAAABMCL+Ld+yxx2qXXXZRWZY66qij9IEPfGBNUutkOp2OPvOZz+iAAw7Q+Pi4ms2mTjjhhNonXcUFF1yg0dFRSRPr7Elak/Z69dVXr/nW3GpXXHGFbr/9dg0PD+uQQw6Z3skCAAAAAABgoxY+pJs3b57OP/98bbHFFup0OvrEJz6hHXbYQR/4wAf0v//7v1q6dKnKcuKPen/3u9/p//7v//TRj35UT3nKU3TcccdpbGxMkvTxj39cf/zHf9yzyZ9//vn62te+1rXY3oUXXqi3vOUtkqR3vOMd2mqrrSRJe+yxhw466CCVZanDDz9cS5YskSTdf//9OvLIIyVJxxxzjLbddtuezRkAAAAAAAB4rKJc/bQtcPvtt+uVr3ylbrzxxonOj0o8Xf3/P/bXy7LUwMCA/u3f/k3veMc7ejHnNV784hfriiuuUFEUevKTn6xtttlGt9xyy5oAiAMOOEDf+ta3NGfOnDW/c9999+kFL3iBFi9erIGBAe2222665ZZbtGzZMr34xS/WRRddtCYFthfuv/9+Pe5xj1vzv5/y7o9oYP6Cnm1vDYIjpjZGT+dRQ3BEjxZy3uAXULZ9E+0ER/QGwREVxoi7bvDBEZnFiAmOqIbgiKkhOKKamRgcUcf+mE11D8ER/UdwRMX+1be3wQRH1BK2QHDE2mNsfMERrRXLdccJH1zzv5csWaKtt946NcakD+kkaeXKlfrMZz6jT33qU2sefk3moIMO0vHHH6/ddtstNZH1cccdd+jUU0/Vj370I91222164IEHtOOOO2qnnXbSG9/4Ru23337h7y1dulQf/OAHdf755+vBBx/UokWLdOSRR+q4447T4OBgT+f82Id0T//rj2hg3no+pKvhYUHqvMmOYYq2sH8PH0r54jEzhjlNgv52jMR+ctuzryUxhn1fov7ZwtvOLxo7HqNw8w72n+tbmH1dBPOO2iZ+UH1sN0Yt80hw26vD5HeKqH/1ybi+0TYzfSd+0N2/Y4tBM3ZQyNnX13FjBO0uWcyMEba7w9eOUX2bLvnMJqJlxnDvV9A/O4/wtfdwHk7UP1UcT9KeeYg4IwvefvzDYo8eHLlx7P3cCPv3od6oYx7xGHXMo4a6x+47N3b1bdr6xtYyQa2QrMsa4YO+OmqW6n2dbN2zIdU48RjdE8zWLPEYboNu7Brmkag3wr5SXLNka5No3pl7f7Z/H+qNcOzka0w9cMy8xhoe9GXqrOzYTh1jTKW+aa1crpv/vccP6VZbuXKlLr74Yv33f/+3fvnLX2rJkiUaHx/X1ltvre22204vfelL9YpXvEKLFi3KvYqNDA/peEhXbQwe0nX/oHpfHtLVj4d01ebBQ7qKY/CQbm08pFsbD+mq9echXaUxeEhXrS8P6aY6Bg/p1sJDurX78pBubTykm9RA1Y7z5s3TwQcfrIMPPji1AQAAAAAAAACTS/47HgAAAAAAAIC68ZAOAAAAAAAA6DMe0gEAAAAAAAB9VnlNOvTG0LJSg+PrtzJhT0MfMoscZxYGllQ2ggVNs6EPQXunWb2v6++3ZxZcDcboNM37aeYX9rar45qFlaPubp+6+UULF9vXYhYddu3B2A2zUHKjGa8kGvVvNOK+TTNGM1hJ1I1h3nI1g/5u8eNwwWYzD8eNUUfQREYmCEKKwxlcYINdUzbo3+7EB3bbBTakxojbO8HYUZskddpmjGgh57aZsxkjXJzZjhE22+t12Jx8zxPrk6sw84vaC/MaU2Nk5xGFPiT62nbb15znmYWfN5TwiR6a7tCHbBBUVN9IyoVM2Xv31PpKJsCiWb2+mWjvPhhcX1f3RMeZPcQSQQ6+vjFjRP0TdYwkFabeiOohdz/P1T2mrnD1UKJmaSbGtrVJ2Br3d2M4G0rdY8OnAu3kGJmaxdw21An6u3m0Xc0SjJGue4L+LmTC1UNhjePqHlvfVA/NsO+sDQsLxnD3eVNDNILX4/pmgi0ayTHC/ukgjajNXCtsWEhimzUEWzhTudw0RqZ+reKbdAAAAAAAAECf8ZAOAAAAAAAA6DMe0gEAAAAAAAB9xkM6AAAAAAAAoM94SAcAAAAAAAD0GemufTb8SEeDoy7e5PdyAXtxAE0yiSxMAMsmkdlE1CCRK5lQFiWzNpIJZUVw9NvtuTMlCG9puCQnG0FTbdxJmsMImtKkdIUprpI00N1eDFRPLZN8qmqUXDZgxnYJZQPN7tihQbO9ATPGYDCG7WtikaL+DfPODJgxokSzARNzlEk/a9iopJyOjTmO+rrjvXq6a8sml3W3t8zcxtvxyRv1d0lp4649GLtl0szca2m1usdwqWrtdvyed1rd/UubBJ2M484cOjaON9iaS0RtxfNotBJjuKS0WsYIrofZdNegv0tby6afxQlqcd8wwk41pbv2KmyxjkR629fUCnUksGbS5F16aqIe6gxkk1mjMUxfczxF6YwuZdIm0keN2a8NJFLtbVJ9UIe4+qbhahaXJh/WPdXrCkkaCMaI6hhJGkyku7p6w40d9Xf1hnstYbqrOT5mat0T1Ru+7olP0qh2cnXFuBkj6m9rJFv3BPWXGcPVVK2gZnFJsO2GSbUPPkuWhXmvgnu/ZG+B5p5p0uQT90ZfV1RPqs/UN65/pu9Ee3COJmunsDbpZd1j39vqH6TTKa4V+jfW9WynAr5JBwAAAAAAAPQZD+kAAAAAAACAPuMhHQAAAAAAANBnPKQDAAAAAAAA+ozgiD4benhcQ0Pja/63XwQ8kFlY2YU42PZEX7tAsVlQPlxAuXpfKV4suRMEH0hSaY7yoh2NEfe1QQ5Bu1sqsjD7L1pxObF+7e8HD9rcGKY9Wiy5YfZpwy1cPGgWHQ4WPx4yCygPDcQrnQ4F2xxuVu8rSUPBKqpDZlVUO3YwxqBdQDkeuxkcJW6x5aivlFtYuZlcWLmdWPy4bQ6oqL8bY9xcRFpB+2g7PknHzUkzFpzUY2bx45H2YDx2cCEabcXzcGOPBe+vW8h5fNwszhy2mcWW3UrJ9gKVuPeY9zFsTYQquHa3+HFjPG4vgv5REMRkY8djmHmYoI84wCK3ULKbd7gosl1s2YwRzKWOgIjCHXtu6B7VPa6ucLVMdCq5vjawwQQzRXNx9UbH1UNBf3ePtmMH83bHXie+HIaHmXsLXS1TRHWgTfpIhG65A9iEQUR1T9MFRJiaxdU9UQiDrXtMzTIc1EMu1GpOM74gRmNHdYzka5moHhowF/FB0x4GR5j3K1P31FHfONm6Jwph8H3N/T+YX1THSL4eGgvGdmO4miUa29Y3rh5qVA/dkql7OsHFomPuMVGozcQP4uZUMGMifKIw80iFXSXrnmgMWz+4sTO1Uy3hE642Mf2DWsaHTLgaeHrqnmLMvFEJfJMOAAAAAAAA6DMe0gEAAAAAAAB9xkM6AAAAAAAAoM94SAcAAAAAAAD0GQ/pAAAAAAAAgD4j3bXPBpaNaWBwdPJOmeQzmTSydLprd3s5YFIEbRKZS4MNUs4Gc2NEIUU2zcyMHSa/2HQxM3Y0rhmiMI/Ew2DLXDhejk056341DZM46tLMXHLZ8GB3ZNAck+LqUlXnDnQn5cwJ2iRprkk5i9pdytm8xljYHiWrDZuopEGboNY9hktES6Wc2ejOqcukuErSeBCr7FJcM+0jJnJw1KWfBe2r2kNh31Um3XUkSDkbacZ9V7Xi9maje4xifOonuhuhY9LF7DU/ircyFzN7lawh5SxqT6ecBfs1k4hmx0gmxEZjuDSzxrg5z20abDC2mV+UZjYxl+5t+pSzRMR5HeyN1B2/QVvT3HTdORDUGzbd1dU9g9XrJF/3uMTW6mNk0vHcNdyVQ9FbYN8ul7AbHZPxbcDOI7xmuTRZk7rbCJJcsymuLpF+TqLuieobKU5sneNqJFv3dNcyw+ai5drnBBdQV7Nk2gfNxT2b+lqHqMbJprtGNYutb1zNEtVOJgl2ZSeuZaIkV1ff2LonqGVcfTNqPitEdc/oeO7Rw3jwHphd6tNd27nPdvHg1dsz9Y1k0uRtfRO3NzN1j6k9U+murj1T97j0WVf3tIKaxdU3mXRX19epUPd0xuPPkBl8kw4AAAAAAADoMx7SAQAAAAAAAH3GQzoAAAAAAACgz3hIBwAAAAAAAPQZD+kAAAAAAACAPiPdtc8ay0bUMKmpayRTzsLkskYumTXqb9NdTXvH9Q/SzzotN4ZJPxvqbm+7FEEbwhL1N53tvg5GMKlDhQmlitrTOUTRL7gXbtqLIPW1GSSfSdJAkAQrSYMmFS1KNJs3GCffzBuo3r7A9J1rklmj9DOX4jqvGacuzymCtDUToRT1lXIpZ01z4DSDY7XRw+Szjvk3nTFzwEf9XV+X2BolxK7sDMdjBH0laWW7u/9Kk4LnUs5WtLrHGDBpZi6VzrVnlEHKmU1xdamKreT9JCEMiE1c96Q4/cylp7o00+h0zCSiSVJzLEpKq55mNtEeJGYH6WSSVJh01yjNTDLJrKavbEJs0D8YV9IGn+6qIMm1NHWP3LkR1CwuIdbVbVF9I8X1kK17XOprcHkyb20qqT5MSZXUcCm4wbRder09z+s4bBJ1T5ReL8UJ9gOm7smkuErSvMHui45LcXV1z/wgmXX+QFybuHTXqMbx9Y1Lu+/u72oWVw8NBRd3V7NEfSfrX4eoZmmbarwdnQTKpdqPlCZVNTjRXd+ovpHi1NdVDZPu2owTYpc3utsHGvH70mzFY0RcYq5tD85H11cu4dx+Hqp+nXQ9w89wybon2q2u7nE1S3TaRXXMRF9T90QJsW4MW38FtYlLr3f1TaY9qmPkE2LD/i7ddQp1T6M1ss4+6xxjyiMAAAAAAAAAmBIe0gEAAAAAAAB9xkM6AAAAAAAAoM94SAcAAAAAAAD0GcERfVasWKXCBA2sk1sUOWo3ix8XLlBioHtSRdA20deMPWgWlG93H3ZuAfBiyLzGaNHW5ELE0bqjpVlEtNGIB48WZzfr3fr3OVwANfliwgWUTVfzWqKF7ZtmseWBZvwih0z7cLN7BdQ5bpFjs4DyJsFiyW4B5YXNeMHOaAHlBabv/GChZMkFR8RzdsERUXsUJiHFARGuvY5wAsct1Dtm/q2nE5yjbvHjsUb1hZV9cIQJfQjem+XtOWHfZUXcHu3X7L6O9l+rE++7tmlvBYv6N8w13Azh/2kukRthX3rUnlzHOVpA2V1T3cLKUXsmIEKKF0tujpnFj117FBwxFr+YYty8SBs0EbzIthnbhUHUERxhFm2eMlffJIIjbH3j6qFm93WoHDRlsg2OiK9lxVBQU9lAj0TdkxXsPzesC8GKzlGT2+NLmShkJnv7ig4Fd9jYcjk4z00o0KCphzJ1j6tvFppaJqpxFpjQB1f3RDVOpr6R4uAI19eFPkT9G2Yl/SETEDHdNY4NjjDtUUiEDYiwwRHdIQyZ+kaKa5yVQRCEJA2Ya34j8TnE1YdRLePrnupjtF3ojvusZlMfTHtGMHgmGMu1u2AGk9cShkHYgIge1j1RuENjLJ60q3tsCFaru38RtEnytUl0vLu+U6hvivaq9f7d1fgmHQAAAAAAANBnPKQDAAAAAAAA+oyHdAAAAAAAAECf8ZAOAAAAAAAA6DMe0gEAAAAAAAB9Rrprn5UrV6o0aTRrNGxMZ9weJZcFqWUTQ5iUs4Hg0DDprjLpZy65LAptKofXN+L20Rs0aWtm/4UBsWYaZSKpJ5V86PrXEWCVTHeN2hvJlLNBl3I2EKW7xmk/C0z62dxmd3smzUySFja62xc24wQel342v+ieh09xjV/jcHDgDJk4qEHzPkZHew1nkSQpehddxpE7zceDA3CkjPfTuIkXXFF2p5G5fb2yjFNfG8HMmzY5rnqSk0sza5nXEiWauZSz8Xb8Tkapgy552l6IepiOF27OXTvdrg7abQqmTUoLrmUmKc2ln8XJrC7lLJ5I1F4kU84UpbjKJJq5lLNWPEYZppyZMUxybKiTPMZcjRMxtYyChOjCpbhG9Y0U1zjudbvUV5MG14jScV3BYYSXHFv3xGN0gk0W7Xj/+xoukWbYo/BfK3ndy6S7Dpj2KMVVihPs55pU+6i+keIk100HVoZ9o/pGipNZNzE10rzCpL42uucd1UKSNGje9DnBBXvQpddv4HWPuW1oLDjxRk3dM2IikUeC/RrVQpI0aGrMqMbJ1DdSXON0zAeLVnRhkdQKrnFjrq+th7rnnfkcI02S7ppRw2c4e3mKPkuatytK15biVPuoTZos9TVR97hk1tHudpfuausb0x7WOK6+semuwRg21d6lvq67xik7pLsCAAAAAAAAMx4P6QAAAAAAAIA+4yEdAAAAAAAA0Gc8pAMAAAAAAAD6jId0AAAAAAAAQJ+R7tpnnZFRdYpHxc4kUs4Kl3JWBMlbpm+ZSCgrysF4c1FqmXxAafRk2IT6qAheiyQVQfSTSyJzqX6d4KW7tMAozUxS+CLTaWa9CltMpjgWQX93NLr0syETOzQURAy5RLSoryTNC9LP5jXidDHXHiW2uhRXl5QWp7vGr9u1Dwc7dsgc64PmXWgE/Zuur/n3mI7NLuvWNgfquEk/Gg+uCy7tdsREbzWDBCV3jNmXkoh+Gzdpa+PB8TRurqlRmpkkjba7x3Zpga69GSYwm5SzsDW8PUiq6TKUSap2KWdRumsyOTI6RDJJsFJ83yhaLuUsbo8SW9MprqNxMmCUaFaOm74muayMUtFcmqlJSivN/X+q3L1fJpE+TH11Ka4mDa4Y7K5xiuC8lfxhbSu44PVE13DJJ9I3grqnbLrkYzNG8Pa6ELzUeZc8z3tV97jDxr0vcbprPLns9Tq6V/l017g9qmXcfTRKcZXiJFeX4rqJqXuibc4ztdocU3tGnyAGXd1jUotdjRPJ1j1RjdMx17dx8xrHgnpoyNRIg7Y9SMd0KZNGJ0iZbbv9YaKgoxpn3NQ3rp6PklyHGvE1dSxI6Jbiczf6vDKpXqbaJ4a2nw+jMdz116XdB/Wy+/xrxwjaXd3japkwydXVN2OuvonHDmsccz932yyjGsc9x5hCqn2njK+xGXyTDgAAAAAAAOgzHtIBAAAAAAAAfcZDOgAAAAAAAKDPeEgHAAAAAAAA9BnBEX1WtloqC7Nw4mpmAVW7wHO0gLJZFLH6Eqz5hZILs8Bz2exutwtTNk17KwiOGDAhE8FimhPtwRiJxcwn+ge/YBbB7+W6pSluYeWgvWEm7RelN4uwB+0DZqcOmhXeo3bXd04Q7iBJcxrd59occ/4NqXoYhAuImGd2VLRY8pwiXjTXLX48GPRvpM5orxMtoGxWsR10gRLBPhm1ix/HY7Rtkku3seC9laROkEozUsQhOMNmjJEgUGKwYxbNNa9lIFhE3J0DLpglDHdx5+gGc8Ex3LU2CuNJ9JVM+IS7ByYWZ3aLLbuV98P7mlkQuTDtblHkMurv+poFlKMaojQhEzLnbhndX5OLnIc1jgvRMvVQEbwHpXtfTGBWGYxt73VBHTMxhrn/B/VQOWDm1zavMTomM7WJknVPpn0Dv9xkrpOur7suu+t4dM139wdb90ShW+Y+NWTroaDuMWO4UIrhsO6JX8sccw4MB+fXoEl2apoxonqol3XPuKkDB0091Aj6N+3JUb2+GQ/CJCRfv44H4QxjLtTK1EOZmtudA1Htb8/FsNWNYTrXcyhMmS2/7DU1c2037WHoVjI4IqhZClMTuM/tYfCkrW8SARFSXOO4gAhXD0XPTszzlLC+mfhB3P7oLqWpuxL4Jh0AAAAAAADQZzykAwAAAAAAAPqMh3QAAAAAAABAn/GQDgAAAAAAAOgzHtIBAAAAAAAAfUa6a7+VpdYZi2VSfUqTilREiUEmCidMOZFURAmxNuUk2R6mwbkkMtOeSADMpAimE8qmOdHMvcZeTSObEOnSYKP2RiINSopT0dwYTTtGd3+XvPX/t3fn4VFW9///X5OEJEBAgoACgqwKKAoIgguKgFVR3IpLrQsW6ccFly9u4Edxqa2KFqVVrCIfcUNFUVtF60IL4spm3VFEUEQQQRSBhCxzfn/kl7sJOe+QI5PcgTwf1+V1xXvOnDlzz7nv+8zhnvOyUqx87ahnJEpZCWX1PClnISmuJeUr1p1uJUEH8r1isZFaLGM/JT3lrdZZ2zM9n2+hUTrdSFvyfeZmvzES23x9JN3sv2F9MoR1bFQb84RTwxFqgef2sDqqXomVmGl0G3/dVh1mqqr1mp7jLvAa7bztM1JcrTT5gP1n8oxxnJWS7j8desvbn1fAPnXG2cn6vDICPoPAZNYq1yuFHaOpOL5SwXq9Gj7vpeo86zvn29eHgGuPlSxqXZO8456w8ZfvWhx6PU/3jln8/TTDTH2txntLPOc+a/xVaHSRep79lzT2dT2jwxd4jmlrPJpuHP/ez9wcn1v9xtd/A+vwHgNhx1dNJ9XX9DBLUkrGMiHfaUOSu63v4an4jh9Uh7HdHJsY6e7eVPuQMUFVpWBsxJ10AAAAAAAAQMyYpAMAAAAAAABixiQdAAAAAAAAEDMm6QAAAAAAAICYMUkHAAAAAAAAxIx01x2BkX5klw+Ye00LKWulO1Zj0p+VShvykkZZbx01HFoYqqZDFePgSwUNZSaRhtRhdIakZ3uxFdhoRCslPduTRjqW7/UkKc1IP0uFYk/Kma/NJWWN7SmIBvR9BqH9I1lN/xZVbKU+muV3wIPXTPSt2WbUlvOyM65HZvN85a3rZej1NeS6m4prtDmu8Jy3QlPNvPspBcdt8NipdnxeVj+rcbXlODfPnXHEMFZdyDk/9HoSIuQaaF1fi719Mmz/+8YEaVY6qZks7kstrc59Z43LjO2e92OOkcywcM8YMwUHYyrGIKkYn9d2tf4tVuN3Wn8dNf8d37y+estaCexGMqv3fGGVtXZ2zVx7uJMOAAAAAAAAiBmTdAAAAAAAAEDMmKQDAAAAAAAAYsYkHQAAAAAAABAzgiNilshIVyKxjY/BWhTVWFgx4VvosF49f9l0o27P9oS1OKNRhwuo21r40aVXfRFmq6y5CLNns7lgaMj22r7oqMG3DqYzdoi1ZKa1qKxve1HSH3xg1VHoKpb3bZPshZILXcVjrcCoo56xvdCzYHChEfqQbuwp73bPgsiV8S3CnF6Nnc9a/LjQWHC10PN+8o06thgdqtCzmHZBwGcr+T9fq6zVn3zb7b7ub59vuxUKEnIcWcfojqqmQ4GcGQBQ9ToUcO0xr6Pm9dUIiEl6zhf1/AeSuUt9x6N1vSw2zk/V1f+McY85ZvHsv4S17+oZYy5feauOwM/ROw4xxzf+zd4xS0j/VTUuLr6D8p0/Q86/UmXncc8537w+bP+4JxXXwIKE/3pez9PuLdZi6/YIseKmwEyQdM9xVOSsdoTxjXGsgAjf+EaSCj0NLzRer9B4875xj/WZW+NX3+drj4u3f9wTMp636rB6je8YraF1+1MvBd8lQ76nWt9/U/HdOmROwJoPSCSN7cVVn1dIGOdUJ2Pew3dEWtfRZNUDbCq2q0gq2maxSnEnHQAAAAAAABAzJukAAAAAAACAmDFJBwAAAAAAAMSMSToAAAAAAAAgZkzSAQAAAAAAADEj3TVmicxMJRKZ2yhkpI4YySre9DMzocy/PZHh2Z7pT0ox6zYS1FyGJ/nFs02yE2FchieBxkhnMYKLvOXNOgKSAc2yocmx1cUKq0lBylmxkbJT5Enw8SVYSdKWpL+fbUlWjMmxym5KZHm3p3mSutKNZFYj1Mv/TxtG2WIjKa3Y85r1Ev4PJtP4wNI8yWrpKYrd8yXHWrvDTn2tKN9MsPP3hU2ehLLNSf9nu8n5z6P5rmIfyTf6ja+sVX5L0n9+s7YXed6j77iQpGJze9VTzsxk5ppORUtBYrbRPYK2J42EsrQ0/w7xlbeSRV2GsVPrVbz4JIwPwBnbg3afleRmbPemnxYbSYnW9urqUCFtloxkVmOMlGEMfX3brXGMlRDrGztJcp6+EDru8fZh4xiw+rvv2Ag9voKSj2t4fJOK86HvPFuyverjG8l/zg+9bviuPVlp/utUvaSRzJqoOHZKt9IJQ8Y9hmKjkqRnXFHPSGa1xmW+03WqUu1Dxj0F1rjHsznfGusaX042ecYhm4xxz2ZnjIeSFcdDVh32mLtin7TKWmM437FhfU9ISap9HKmv3lRVf9Gw75LGWMH6nuo7t6cH1uH7bm2c36zrlzzlE0aatJmeakyL+FrtzCR4Y8ziu75a6fVJK9V+2x0t4ZKkuwIAAAAAAAA7OibpAAAAAAAAgJgxSQcAAAAAAADEjEk6AAAAAAAAIGZM0gEAAAAAAAAxI901ZonsLCXSsrdRyIiDMRJN5EltSRhJsN40M8mbUOaM1DI7/azqKWfJTH/ZZD0jBcizPelJpal0u+clzZQz460HpZzVNDPNrOopZ2aKq5H2U1js31EF6RX7SIGRRLbF3O5L+vRHAKUZyWDpZlZXRcXGMZP0/NuGVTbTSC7b4klbq2e02dqe7omxSktRtJXvFYuNjm0ldRV69pOVAJZvJD/50lbNlDMr9dWz3Sq7udiqu2I/sxPRjGPAc2wUGmWLrZQzz3Hn2yZVEnKWogRgr4DzYUiiZHC6q+/cbqRdJusZDfQllFmnDyPpK82TXJY0hl0J8zpvpJz6toemmRZVPA+ZKbO1JN3VSpP37g/rfRv7yVe3leJqb7fGPb4xS9XHN9Z2q/9a/d17bKQgPTkoCbay7dvLGt8YxUPOqdZ52TqP+875W9L8Za3ryWZPumA941j0jQksSeMDs8Yyvut/YcKX4S5le8Y3kpTvS7U3xmTpRrJlbRn3FBpfCnxjHGt8UyB/Hb7xiTXu8aW4WnVY42Vruz/d1XgvAdut7wnW9wp/uqu3aCWpr7Xki1lIqr3xvdM61/pOQwnj+2/CuJwnfJ+BdXgZDfQ2L3B8o0L/OcQ3l5EoMt6MZ3wjSc5zTk1YKbNWCnZV0l2TTtq8zWKV4k46AAAAAAAAIGZM0gEAAAAAAAAxY5IOAAAAAAAAiBmTdAAAAAAAAEDMCI6IWSI7W4n0+mU2BCxu6QmIKNlece7VWQslW4ESnvLOWLA5ZKFkq7y5UHKmtYByxfduLaBsBUf4Fxf3Fg1aqDN4EeZUhE/41rA0F/s0NnsXUPY3pKjYCAsw+tmW4oqnmgxjFXZruxUGEcK3WHKhsaivvb3ih+4LOJCkbGNh5Xre4AhrMWgrUMJYLNVbR9jCyr7Fks2Fps2FlT1hIcYBZu0/3/Z8c4HtqgdKWAsl/1zsD/HZVORZyLnYX0d+sfFePNutBZTNhZU9JwvrGHXG9hStse3nW/w4cOF476L05iL4/jeT9JRPGGVDFla2jgGT53qeSDfOb8Z7lLXdF+RkLaDsWShZkhKePpkwypoLJfsWXA4Nk/CNe6zxjTVGSsW4xxuYVfVxjCQlrfGQp57QcY/z9NWQYKyS7Z5F2K1xj3Xc+ZoXul57tQVHGJuN86Hv/GkGRFjn6zT/MZOfVvH4yjDKWkEJIZIBIQdW8EHINdoa31jbMz0nVd9YSLLHLL5xYCrGNxZrPxUb97j493XVg7Ek/xgnNDDLN8axxj0bPeMbSdroCdLKs8Y3RUaQlmfsX2gERFjfK3zHqO/7SskD/s3B230CznGh3wO9oQ+BQYa+7db1wfo6lfCe+1Iw7jGu51aopR126RnjFBtvxviukPCVt8Y9VqCEz1bjnkTx9p/TuZMOAAAAAAAAiBmTdAAAAAAAAEDMmKQDAAAAAAAAYsYkHQAAAAAAABAzJukAAAAAAACAmJHuGjPXoL5cRv3KC4WmnHm2WyldVsqZr7yVlGannxlJPZ56zGRWq46gdFfvZhX76jDK+lLVJHmnuc0AQDPlrGICTEoCGANTzpwn1afYSFsqMvqClX62xfMe0wLSuyQp6Uu2NHa2lchV6Eko9iWwSVJ+ur8z+NK0stKslLMC73Z/ypmR7mql3XpSX62yoYoDUiyt5DJfgpqZcmakMIUk6W72JJFJ0hbPa242klk3GXX4Es02Ffnr2Gxs96WcbTGOFyvlzHc8msdzaLprSIhVSJJjYOqjr+tZSZUJawTjOQz8qWVh7MQ24xro2WyluFrX6ISRFppW5HmTvm2yE1udJ30sJemuqRA47vF+BmY6qZUc5xmbWOMbq26jvC9lLyS93truG8dUWofnmLHGPdZx501gDk1xDkmCtgSk2oecJ81xj9FvtqT5d1R6WsUdm1aNKa6+MZJkjHuMa25+upFUn6w4xrHGLNnGeMiX5Gols1qprz7VOe5JGveyWGMk35jFTsw1xgqezybfGDuZ4x7PQe1La5WkPGs85BnL5HvGMZKUZ6W7FnnSXa1xj5HY6h33GH29OlPtzSp8x7R5nTLq8J1TreRu65rrPakaqfbmydazr41zlj2Gq/hAmpUym2GMWYzv/glP3zFTVK2xjK+8Nb7ZjnGPKyLdFQAAAAAAANjhMUkHAAAAAAAAxIxJOgAAAAAAACBmTNIBAAAAAAAAMWOSDgAAAAAAAIgZ6a4xS+ZkKVkv+5c92UqP8QW8GGmcVoKaL53FSprxlZX8aWZWPVZZc3tAMquVfuYLSzJfz0g5872mnQAYtj2IN+XMiN4xU848Ra2UsyIjPdVIAUpFolmRJyHLSjkrMnaqL/Gqfro/iWxz0p945Usus1LOrNRXX/lUpJylYj9brOS4YuPfenwpZ6FpvL4EOt9nWPn2inX40lorq2OzJ7ks36rDSD/L96ScFVjHkZF+5ks5sz4X6zi3I7n8m7eXed6zkiM95yHfNsk+X5uJXF5Wiqin2oA0M8mfAJo0Er/SjHNtwijvij07xUg5SxhJZN7tVlkzFbiaOk7A+EaSdyxjfi5Wcqzv8zJT7a2xQirGPf7m+cc9xnnZf3qSL9TTHKtZx6hnu1m2Om8FSMG4x3f+tNJdC40U10QKxj3WedyXeFmY4W+Hde3JS684lrHGJlYCaFZaxXGIneLqHw+leVJYzVR7X0S34hnj+ISl2hufV0iqvVHWGrP4xjgFAWUlf1K9L61VstNdCzxjmUJj3GN9r/B9D7HGBHGMe7xp14Hf97wfjfUeraRUz7XY1x///8Lezd7v/uYYzro2Vqw7WeRvR5qRhm6OWTzjIausOZbxjJ2qY3yTLPSf20JwJx0AAAAAAAAQMybpAAAAAAAAgJgxSQcAAAAAAADEjEk6AAAAAAAAIGYER8SsqGE9pWX6F6iPVH0tU0mS8yy4bC7eGxIcYS4MHBgc4anHXGzZDGyoenCEvd1Th7HYsrFGrH8BZXNfG9t9bz3wMw9irIPpfIvSG20uNhYoL0gYH5iHtUivuYCyZ8f6FlWWpPw0/weZmV5x8eNNaVbogz+wIctTR4ax+LG1KHJGWsWFS0MXSk73LMKcVl2r/8sO6Sg2Orzvc7QWULa2F3kSWwqN17MWRfYtpl1gJMFYYRC+IAdrke4tRuiDLySiwFiE2VpAudiz3RVbQTDWAur+zQnP5xW8Rrd3AeWqL1Bcst2zzbr2WOcyz3tJWsdGQBhE0KLPknynkDTj80paoQ/G2sMJz/naXGzZeuu+17TKmjvbv3m7met/BwR9WOOKkM88cGxilfedckIDG1Ix7vEFWISMb6TQhdKNDhJwi4C9qLfnvRjHgHU+9J0/i2UFRBjtCGAGMBljGd9261qXbYRg+YIBfOMYSco0xj31POONDGPsZI17fNtDxjdSzY9xQsY3kj9IyxqnmuMez3ZzvGH0Bd8Yp8CowxoP+V7TCo7wBURI/pAIKzDLCmxJevaf7/tKyQP+zSlhHf++7aEBgr7vkta53WiGryFm0IoVzBRwbreuX2mecYU1jkkzgrHMsYynnuDgCM8ONMc3lioULyooCKvTgzvpAAAAAAAAgJgxSQcAAAAAAADEjEk6AAAAAAAAIGZM0gEAAAAAAAAxY5IOAAAAAAAAiBnprjErapihRNYv/BjMhLKAsilIeAlJ6ZNSk+7qr8MoayTQ+MqbaTpW+llIEqGZcujZGJoi5qvaigCyUs485ZNGEqES2z+/X2ymnBnprp6Epy1p/g+mXro/SijTk0Zmlc0w0sV8CbFWWStZqZ6nHVZqmVVHmvGa1SVpxhb7+dJ4rc/cl+Iq+dPWrBQxXyJaSd0V21FovJ4vxbWkvCdhz0pmNVLHfHVbaWbFRVbKmSfN0Eo5s7qHGW1plPcx+qTzfF5mImLI9aGKzapMIvg6VfE9WtcpI+TQm35mJZFZqa/WYe6rxzwlWAFvKaijxgWMe6yxSSrqMBPvzOu/Lzk2sA7PKcccOwUktgYnxPraF5JeL6M7hcZJ+4pb5zcr9dWXam8cBNZ52eJNmTaSPovSjWtPRsXy9Yr97duSboyHPCcoX8K8ZCfVZ3rGSdaYxUp9Tfd8vtbYyVJbxj1WuqtvzGKNN6w6fGMnK5nVV1byjzesFFdr3OMdOwWMbyT/GMdMr7fSXT3XRju9PizVPmzcY1QRkGov63ztKW6mxvs3+9+ief0yzocBKeRmYqsvgdVKqTc+LzPd1ZfMao2RAhLpq2PcU7Rl+6fYuJMOAAAAAAAAiBmTdAAAAAAAAEDMmKQDAAAAAAAAYsYkHQAAAAAAABAzJukAAAAAAACAmJHuGrPCBmlSdpm50tBUTw9vYJCZ7mps9yXChKTJBm63E1GN7Z6EHDMRLWB7aDt86TuhdaTiM/dWa3y4vhRXSZInfcdK/02a0YzWa3o+rwx/Q6z0s0JPopmZUGYkpeV7yvsSxyqrO92z3UooSwTUbaW4WqxkNW/ZwES0kCRXX5pZSR2ez9xKRLOSfj3t8CWOSVKxtd1Tt1XWSi7z9UmzHVZiqy9F0Hw94zjypAs6K4HZ2m51m1Skd3pe0k7BNFJOrYuSr46AdliB1NZ2335NWCmYxQFJZMEprsZ2z36yEtHMdFffvjbL+reHvJ4p4Bpodg/fWwkoKxkJrKka9/jGGwHjG6t8yPjG2p6KsZN1PFfruMeb7mqUNY47/+drnMOtZljXNd+4x0yv9x/ovrTwNGNsUi+96mMWa9zjKyv5xyzpgan2vu1m2Vo+vrHrrljeGitY/caX2GqNWcxxSMC4x0qk99Vtjc9DkuqtOnwprpIxxrHGN1ZaaMC4IljAuEdWordnm5niaoWW+64PVqqqdY3xJrNWfXwj+T/H4PGN9Sa9yaxG+8yLt6+sv2iQreooTN/+++C4kw4AAAAAAACIGZN0AAAAAAAAQMyYpAMAAAAAAABixiQdAAAAAAAAEDMm6QAAAAAAAICYke4as8IGCan+L0udCQqrCU05C0mrCUmIteoJTVALSUoLSZkNeD2zHYEpZ96AstAu4f3AwpL+vPFCnmSmkqr9lSSN6FjnSbkxEz2N10zzJJcVpRnJYOb2qqecWcms6Z667bJh6WfedhjbQ9Ngt5f1eVl8rbPqMNPxvAllYXX4Es3MlFmjbl8amVnWSDnztc9KM7Piu4JSzoz3GJKaZQq5nljXDev05D1/WnFm/s3et2icw9OsyEZfMquVRBaQ2GoetmbKmZVc5tlofeYhib6Bp5XqOg0FXwNTkGofkhBrjW9Cgs+Dxxu+8mZ6sn+79z0GJtJ7j1HrOA8ZNwaPezxVWOdOq25fqr35gkaCpTHe8J/zjTqMZFbfGMdKdy2wxj2eg9Qam4SMe6wxiFWHr3zoR15bxj0hrTDHugGpr8HjHl/fs1JVzXFP1fuvuZ88dViJ9M6o2zvGsQJ9Q9PrAy40Qdck8xYoqyGefW2VtBJbQ1LtrXN+ChLp/edlox2hifTeMVXgmDbg4N2e002hNUYIwJ10AAAAAAAAQMyYpAMAAAAAAABixiQdAAAAAAAAEDMm6QAAAAAAAICYERwRs+IGCSViDI5ISR2hgRIBi4vbizYHvF5IoIRZtuphEOb+CNkeuFqlr7iVG2EurFzljTZnvHnnWRw0Ye5T//Ziz6KyVh3m4sKe7VbIhG8R1pLtVQ+OSEUdlu1fkjRM6Pqp1oLGIWV9fdgMmQhYWNnMVAlYyNmuo+ppBnZZ/2ZvIkJoHSEfZOjxH7IYfMA/ETrr+DLeu2+xZGvhYmut6qDQh4BFjoPKKnBR5Kqf9iotH1RHDavWcU8KwidCxkOpGDuFjjeCArOMa6OvbrOOVIw9U3HOss6HAVdS65wvY4HwpG9RdGOfWtcv72Lw1vnQDIOoeh32eMg3ZvEWDRrL1IVxT8j4JlV1ePtTQECXZIxPAsc93u0BYyRJ/uCIwO83wdfGkMIBAQH2y/nSFoyiVihjSNhVKkIfAs619hhk+0MpLHGPe4qtE2QA7qQDAAAAAAAAYsYkHQAAAAAAABAzJukAAAAAAACAmDFJBwAAAAAAAMSMSToAAAAAAAAgZqS7xqwoW1L9X/jkgOCQoCQto+7gOkLqTkVCbHDaWsXYluB2+FLOQpLjrPKp2NeBKWcJz4ua6W6+tKVK6vbFgJmJuWaiUdXTxZJmuquvaVUva76mlbbmr8KoNzBGqKbjFlNwAjATylLQDjN5LyAhNii5LLAOb3JZSFmrfIpSzlJygvf1SeMgNV8u3VevUdY6h3hTzowUMTNhz1eH0Y6AvpeKxDGznjhOIdV1GkpBd4xl3JOK63wqUu3NOnyJo2F1eNNdjaL2e0xBx/HG3YYlQTvvBSIwfdIcy3jKh457fKmq1rjMF2stY2wRMr4x2mExq0jFZ17Lxz1BqfYBrxkyvjHbEXjt975m6LgiZNwTMKYyxzepSLW3pOLcHvIdM+QcLv/4xOqPQemu1ZhqH6rGxz3b0e6iFLxn7qQDAAAAAAAAYsYkHQAAAAAAABAzJukAAAAAAACAmDFJBwAAAAAAAMSMSToAAAAAAAAgZqS7xqw4U0pkbaNQqlJVPYKCi6o15SwsBiUlSWmpqMO3MTQprbrSqswXNF7Pl4ho1RGYKOfdJ1Y6k7k/PIlXgZ+XN9E3Ff068DM0E9S2V6r6UsqinLeqNrR5Ie1IRWpmSGJrKlJSQ1OzQhLbzPSusJTTgCq8h4ZZ1jpP+toR+s+JvvRJK/Ux5DMwPxerjopvPvgQTVEarE+1BSWG1ltN58OUnMZCrzGGkHFP0HUt9FocUodRRXW1I/Tz8u0+bxK3VVj+VHv7vVhRmgEnxOLAcY9nsBA6lvRuTUW/ZtxTtWpTMe5JxXUjNJE+FWMnH2NsYn7fCBp/Wdur8Yu073tFWmCavG+f+JLuJTvR1zPuSRidLyQtOHzMEjDuSdV4yKNag6CrUHdx0fa/DHfSAQAAAAAAADFjkg4AAAAAAACIGZN0AAAAAAAAQMyYpAMAAAAAAABiRnBEzJJZTsVZv3B1w2pcBzMVdccRSuFtRypeM3RRZG8d1sLA21mv5F/s0ywasDhrQJslKWGtDOxbFDUl+zSgrFSNKxdbrxewkHNw3amoJEB1LsJqvmZ1rShffeXNxY9D6q3GhZyrdTFdn7TQhYurvuiw+VYCFpquekxNeB0hK4an5HOJ5Ritpnpr+vyWotes6XFP8O4PCcwyBIV/1Y5Lrj3uSQYco9b4IQUhDKnpN4Ht216hgW8pec1UVBKgWhefr+E3E8e4J2jMErA9cN+lIrggLEfPGvsblfhCIlLR5hSclsODT3aicc921J0s2v6GcScdAAAAAAAAEDMm6QAAAAAAAICYMUkHAAAAAAAAxIxJOgAAAAAAACBmTNIBAAAAAAAAMSPdNWbJDClRL+5WVNHOlCZrSUXUTCqS0lIhFWFrRuOC2xywXxPVmVBW08lgNZ0mu7OprtSm2pIGlYKDP/iUVY3v3fd2zPaZiYgVnxCeLpaKC8r2VxFSb0peLo50151JLR/j+MQy7qkt47WQ9MmAhphtCz2+QsY91TlWqLZ0V8Y326Wmz9e1JUkzVcmsQa+Zgjqsqn3jHqtwQFR98Dky4AnmsKyG93Xwy+1AY5xkCmbYuJMOAAAAAAAAiBmTdAAAAAAAAEDMmKQDAAAAAAAAYsYkHQAAAAAAABAzJukAAAAAAACAmNXpdNfCwkLdcccdeuSRR7Rs2TI1b95cJ598sm644QY1adKkRtrgMpxcxg4UV1JddtCQqGpLZq1OqQiwDU1srCU7qpY0A7VUSpKtqlNtb59HatInjc3VGFFG2iqqXS2/Hu1U18uQJNjQ8Y35mjW7A3eqzwspx/im5gQfi76E2GpMk7Yw7kmdVMzt1Nk76fLz8zVw4EBdc801+vzzz7XPPvvoxx9/1MSJE3XggQdq3bp1cTcRAAAAAAAAdUSdnaQbPXq03njjDXXs2FFff/21FixYoGXLlmnvvffWkiVLdP7558fdRAAAAAAAANQRdXKS7ssvv9TkyZMlSY8//rhatWolSdp111312GOPSZJmzJihTz/9NLY2AgAAAAAAoO6ok5N006ZNU1FRkbp3764+ffqUe+yAAw7QPvvsI+ecnnnmmZhaCAAAAAAAgLqkTgZHvPrqq5Kk4447zvv4oEGD9PHHH+udd96p9ra4tJL/gJ1CLV8wtNYvnAugymr9Qum1vX0Aqq6Wjx8Y3wA7j1o/vrHsqO1OsVTM7dTJ6aEPPvhAkrTvvvt6H+/SpYskadWqVTXWJgAAAAAAANRdde5Ouo0bN+rHH3+UJDVr1sxbJjc3V5L07bffBtXdqFGjbZZxjn/qAgAAAAAA2Jm1b99eiUTltxn+/PPP5f6/zt1JV3YHbGuSbs2aNUomkzXSLgAAAAAAANRdde5Ounr16kV/p6ene8sUFxdLkpLJ5DZnPcvaegbU57vvvtPuu+/+39fatKnK9QMm1gAAUFO4IRxATWF8A6CmML5BCmw9v7N06VLttttuQXXUuUm6Zs2aKSsrS1u2bNG6deu8ZdavXy9JysnJCZqkq4offvih3P9/e+vtKa0fAAAAAAAA8frhhx+CJ+nq3M9dJally5aSpLVr13ofL52ka9OmTY21CQAAAAAAAHVXnZyka926tSRp8eLF3sc/+ugjSf9NeQUAAAAAAACqU537uaskHXzwwXrzzTc1c+ZMjRs3rsLjr732miRp0KBBKX/tzp0765NPPon+v3fv3kokElq2bFnKXwtIpfbt20sSfRW1Gv0UOwr6KnYE9FPsCOin2FHQV3d+yWSy3BJnnTt3Dq4j4Zyrc0skLly4MJoc++CDD7TvvvtGj73xxhvq37+/srKytGzZsuinsdWlUaNGkqoWOgHEib6KHQH9FDsK+ip2BPRT7Ajop9hR0FdRFXXy564HHHCAhg4dKuecTjvtNK1Zs0aS9P333+vss8+WJF1wwQXVPkEHAAAAAAAASHV0kk6SJk+erPbt2+uTTz5R69at1adPH3Xs2FHLli1T//79deutt8bdRAAAAAAAANQRdXaSbrfddtP8+fN14YUXqkWLFvroo4+0xx576JZbbtGsWbOUlZUVdxMBAAAAAABQR9TJNelqE36Xjh0FfRU7AvopdhT0VewI6KfYEdBPsaOgr6Iq6uyddAAAAAAAAEBtwSQdAAAAAAAAEDMm6QAAAAAAAICYsSYdAAAAAAAAEDPupAMAAAAAAABixiQdAAAAAAAAEDMm6QAAAAAAAICYMUkHAAAAAAAAxIxJOgAAAAAAACBmTNIBAAAAAAAAMWOSDgAAAAAAAIgZk3QAAAAAAABAzJikAwAAAAAAAGLGJB0AAAAAAAAQMybpAAAAAAAAgJgxSRejwsJC3XLLLerWrZvq16+vtm3b6rLLLtOPP/4Yd9NQh8yfP1+jRo1S//791aJFC7Vo0UKDBg3SlClTzOfQd1EbfPTRR7rtttt06KGHatddd9X1119foQx9FXHIz8/Xddddpy5duig7O1stW7bU2WefrZUrV5rPoa+iNti0aZPGjx+voUOHqn379mratKkGDx6s8ePHq7i42Hwe/Rc7Cvoqatp3332nUaNGqUePHsrJyVG3bt105pln6osvvjCfQz+t4xxikZeX5w499FAnyaWnp7sDDjjANWrUyElynTt3dmvXro27iagDPvroIyfJSXL169d3ffr0cS1btoy2nXbaaRWeQ99F3H766Sd3+umnR/00KyvL7bPPPm7y5MnlytFXEYfly5e7Xr16OUmuQYMGrlevXi4zM9NJcq1bt3b/+c9/KjyHvorq9v3337uLL77YZWRkOElu+fLlFcrMnz/f7b777k6Sy83Ndf369XNt2rSJzrWHHXaY27BhQ4Xn0X+RSlXpq845l5+f7xKJRNQ/t/4vJyenwnPoq0iVqvbTd955xzVu3NhJcs2aNXN9+/Z1TZo0ib57vfzyyxWeQz8Fk3QxueCCC5wk17FjR7dy5UrnnHNr1651e++9t5Pkhg0bFnMLUVcMGDDAPffcc66goCDaNmnSpGiQ88wzz5QrT99FnN577z3Xvn17J8n17NnTTZ482f3444/esvRV1LRkMun69u3rJLlf//rXLi8vzznn3MaNG92JJ57oJLl27dq5wsLCcs+jr6K65OXludtuu83tsssu5SYwfF8o16xZ43Jzc92kSZPKjQlef/1117x5cyfJXXfddRWeR/9FKoT0VeecW7x4sZPkWrZs6QYPHlzhv+OOO67Cc+ir2F4h/XTjxo3RmHXs2LGuqKjIOedcQUGB+/3vf+8kuTZt2kRjhVL0UzBJF4OlS5dGs+7z5s0r99iCBQucJJdIJNwnn3wSUwsB5wYNGuQkuXPPPTfaRt9FnNauXRvd1XH66ae7/Px8syx9FXF49NFHoy+NW/fPvLw8t8ceezhJ7tFHH42201dRnd577z0nyTVv3tzdfffd25z4WLdunXd76T/e5eTkuOLi4mg7/RepEtpXZ86c6SS5MWPGVKl++ipSIaSfvvrqq06S69Kli0smk+Ue27x5s8vNzXWS3Ntvvx1tp5/COedYky4G06ZNU1FRkbp3764+ffqUe+yAAw7QPvvsI+ecnnnmmZhaCEg9e/aUJH311VfRNvou4jR69GitWLFCxxxzjKZNm6asrCyzLH0VcZg+fbokafjw4RX6Z3Z2tn7zm99Ikh544IFoO30V1Wm//fbTTTfdpKVLl+qiiy7aZvmmTZt6tx999NGSpI0bN2r58uXRdvovUiW0r3755ZeSpM6dO1epfvoqUiGkn77//vuSpH333VeJRKLcY/Xr11enTp0kSStWrIi2008hERwRi1dffVWSdNxxx3kfHzRokCTpnXfeqbE2AVtbsmSJJGnPPfeMttF3EZdly5Zp2rRpyszM1F133VVhsLM1+iri8NFHH0mSunbt6n384IMPlqRyi0XTV1Gd0tLSdN1116lRo0bbVU9OTk70908//RT9Tf9FqoT21aVLl0qq+iQdfRWpENJPd911V0nSxx9/LOdcuceKioq0bNkySeXHDPRTSEzSxeKDDz6QVDKr7tOlSxdJ0qpVq2qsTUBZH374oV577TVJJXeElKLvIi6PPvqoioqKdOyxx6pt27Z64IEHdNZZZ+mUU07RLbfcos8++6xcefoq4vDdd99JKvkXcp/dd99dkvTtt9+qoKBAEn0VO4ZPPvlEkpRIJLT33ntH2+m/iEvonXT0VdS0E044QY0aNdKnn36qG264IZqoSyaTuu6667R27VoNGjSoXJ+kn0Jikq7Gbdy4MYpObtasmbdMbm6upJJBPFCTvv32W91333068sgjtWnTJl155ZU67LDDJNF3Ea/SSePu3btrwIABGjlypB599FE9/fTTuuaaa9SzZ0/93//9nyT6KuKz1157SZIWL17sfby4uFhSyQB93bp19FXsMGbNmiVJ2nvvvdWgQQNJnGsRr7J30mVnZ6tDhw465ZRTovFCWfRVxCE3N1ePPfaYMjMzddNNN6lfv36655571L9/f916663q1auXHn/88ag8/RSlmKSrYT///HP097YOvjVr1iiZTNZIu1C3TZw4UYlEQq1bt9b555+vzZs364knntD48eOjMvRdxGnRokWSpDvuuEPOOb355pvauHGjFi5cqKOOOkp5eXkaOXKk/vOf/9BXEZvSn7Pef//92rJlS7nHCgsLdfPNN0f/n5ubS1/FDmHdunWaOHGiJOnSSy+NttN/Eafs7Gwdcsgh2m+//bTrrrtq+fLlevrpp3XkkUfqjDPOiP5RRKKvIj5Dhw7VnDlzlJWVpXnz5mnUqFF66623tNtuu2n69Olq3rx5VJZ+ilJM0tWwevXqRX+np6d7y5T9l/ZtrbsEpEKbNm00aNAg7bfffmrYsKF+/vlnXXXVVZoyZUpUhr6LuOTl5Wnjxo2SpJYtW+rFF1/UwQcfrIYNG6pXr1568cUX1aNHDyWTSV1zzTX0VcTmf//3f9WwYUOtWLFCRx55pBYsWKDNmzfr/fff18knn6xXXnlFktSoUSNlZ2fTV7FDuPbaa7VhwwbtscceOvfcc6Pt9F/EacGCBXrjjTf05ptvauXKlVqxYoXOO+88SdLjjz+uG2+8MSpLX0VcXnnlFZ144okqLCzUySefrCuvvFK77767vvvuO+2777567LHHorL0U5Rikq6GNWvWLEp8W7dunbfM+vXrJZUs0svBh5pw8skn67XXXtP777+vNWvW6K9//atWr16t8847T7fddpsk+i7iU3rrv1RyJ13pQryl0tLSNHr0aEnSW2+9RV9FbFq2bKmnnnpKTZs21dy5c9WnTx81bNhQPXr00CeffKKxY8dKkvbZZx9JnFdR+02bNk1/+9vflEgkNHny5HKpxfRf1CatW7fW5MmTdfbZZ0uS/vKXvygvL08SfRXxePfdd3X88cfrp59+0vPPP68ZM2Zo/Pjx+vLLL3X++ecrPz9fZ599tubOnSuJfor/YpIuBi1btpQkrV271vt46cHXpk2bGmsTUKpBgwYaNWqUJk2aJEm64YYbogsFfRdxaN68uTIyMiRJbdu29Zbp1q2bpJLUwdWrV9NXEZtjjjlG//nPf/TnP/9Zp59+un7/+9/rgQce0IIFC/TVV19Jko488sioPH0VtdWHH36okSNHSiq5S/Too4+uUIb+i9qm9B/tfvrpJ33++efRdvoqatrYsWO1ZcsWXXbZZRoyZEi0vX79+po0aZIGDhyoZDKpP/zhD9Fj9FNITNLFonXr1pLshaU/+ugjSf9NbwHicOqppyo9PV35+fnRemD0XcQhIyMjmpxbuXKlt0zZnwU0bNiQvopYtWnTRqNHj9bjjz+u++67TyNGjFAikdDMmTMlldy9XIq+itpoxYoVGjJkiDZv3qwhQ4bohhtu8Jaj/6K26dixY/R3fn5+9Dd9FTWpsLAwukPu2GOPrfB4IpHQqaeeKklauHBhtJ1+ColJuliULixdOljfWmkq0aBBg2qsTcDWGjRoEN29VLr+AX0XcSkdjMybN8/7+JIlSySVTI40atSIvopaZ8SIEVq/fr2OO+449ejRI9pOX0Vts3btWv3qV7/SN998o4MOOkhPPfWUuT4S/Re1Tendc+np6erevXu0nb6KmvTzzz+rqKhIksyAh8aNG0sq+elqKfopJCbpYnHaaadJkubPnx/Nhpd64403tHTpUmVlZemkk06Ko3mAJOmFF16I0gl79+4tib6L+FxyySWSpMmTJ0drzJT14IMPSvrvv1bSV1FbFBUV6eqrr9Yzzzyjxo0b689//nO5x+mrqE02btyoIUOGaPHixerWrZteeOEFNWjQwCxP/0Vtc/fdd0uS+vbtW67v0ldRk5o2bRr9CuTZZ5/1lpk1a5Yk6cADD4y20U8hSXKIxdChQ50k161bN/fdd98555xbs2aNa9++vZPkLrvssphbiLrg2WefdY888ojbsGFDue0zZ850u+22m5PkRo0aVe4x+i7i0q9fPyfJDRw40P3www/R9ltvvdVJcjk5OW7VqlXRdvoq4rR+/Xo3ffp0d+CBB0b988UXX/SWpa+ipkhyktzy5csrPLZlyxY3ePBgJ8l17tzZffvtt1Wqk/6L6lBZX33qqafcP//5T1dcXBxt27hxo7v88sudJJeVleU+/PDDCs+jryLVKuunf/rTn5wkl56e7v72t79F25PJpLv33ntdIpFwGRkZ7q233ir3PPopmKSLyerVq6MDLSMjw/Xu3ds1atTISXL9+/d3+fn5cTcRdcChhx7qJLlEIuHatWvn+vXr53Jzc6MLzrHHHuvy8vLKPYe+i7h88803rnfv3tGAp0ePHq558+ZOkmvSpEmFCRD6KuIwbtw416RJE5eWlhadS/v06eOWLFliPoe+ippS2RfKM888M3q8e/fubvDgwd7/xo8fX+559F9Uh8r66v777x9Nxu27776uV69ern79+k6Sa9q0qZs+fbq3TvoqUq2yflpcXOxOPfXUqEyzZs1cv379XIsWLZwkl5mZ6SZNmlThefRTMEkXo7Vr17oLL7zQtWrVymVnZ7uuXbu6W265xRUUFMTdNNQRS5cudWPHjnUDBw50bdu2dQ0aNHDdu3d3w4YNM+/4cI6+i/hs3rzZXX311a5Tp04uKyvLtWjRwp122mnmBAh9FTVt9OjRrmnTpm6fffZxl1xyifvnP//pCgsLt/k8+ipqQmVfKNu1axc9Xtl/w4cPr/Bc+i9SrbK++u6777qRI0e63r17u2bNmrkmTZq4Qw891I0ePdqtWbOm0nrpq0ilyvppqZkzZ7rjjz/ederUKfqudc4557jPP//cfA79tG5LOOfcL/uhLAAAAAAAAIBUIDgCAAAAAAAAiBmTdAAAAAAAAEDMmKQDAAAAAAAAYsYkHQAAAAAAABAzJukAAAAAAACAmDFJBwAAAAAAAMSMSToAAAAAAAAgZkzSAQAAAAAAADFjkg4AAAAAAACIGZN0AAAAAAAAQMyYpAMAAAAAAABixiQdAAAAAAAAEDMm6QAAAAAAAICYMUkHAAAAAAAAxIxJOgAAAAAAACBmTNIBAAAAAAAAMcuIuwEAAABATfn3v/+t119/XZI0fPhw7bnnnjG3qOY553TTTTdJktq1a6dzzjkn5hYBAABJSjjnXNyNAAAA2NG1b99ey5cvV3Z2tvLy8uJuDgzXXXedbr75ZknS7Nmzdfjhh8fcoppXXFysjIySf6sfMGCA/v3vf8fcIgAAIPFzVwAAUI02bNigJ598Uuedd5569Oih1q1bKzMzUzk5OWrTpo0OP/xwXXHFFZo1a1aF586ZM0eJREKJREIXXHDBdrdlwYIFGjNmjHr37q3dd99dWVlZatOmjQ477DCNHz9eX3/99Xa/Rl3y+uuvq1WrVurRo4dWr14dd3NqPfYXAADYFn7uCgAAUm7Tpk268847NWHCBK1fv77C44WFhdq0aZO++eYbvf7665o7d67efffdamnLihUrNHr0aD399NMVHvvmm2/0zTffaO7cuRo3bpwuv/xyXX/99crMzKyWtuxMbr/9dq1atUqrVq3StGnTNHr06LibVKuxvwAAwLYwSQcAAFJq6dKlOvHEE/XRRx9Jkho3bqyTTz5ZgwYNUvv27ZWbm6sff/xRq1ev1ueff66///3vGjZsWLW05YMPPtDRRx+tVatWSZJ69uyp3/3ud9pvv/2Um5urb775Ri+99JIeeughbdiwQX/60580b948Pfvss8rJyamWNu0sioqKor+Li4u9ZT7++GP9/e9/V79+/TRw4MCaalqtVJX9BQAA6jYm6QAAQMosXbpUffv21bp16yRJl1xyia6//no1bdrUfM6YMWOqpS2rVq3S4MGD9f333ys9PV233367LrvsMiUSiahM9+7ddcwxx+jKK6/UsGHDNG/ePL322ms666yz9Oyzz1ZLu3YWo0eP1oIFC9SqVSudeeaZFR5/9NFHddZZZ0mSpk6dWsOtq322tb8AAACYpAMAACmxadMmnXDCCVq3bp0SiYQefPDBWFMjzznnHH3//feSpIkTJ+qiiy4yy7Zp00Yvv/yyDjroIC1evFjPPfecHn74YZ199tk11dwdzpFHHhntXx/uFitvW/sLAACA4AgAAJASEydO1McffyxJuuqqq2KdoJszZ45effVVSVKvXr0qnaAr1aRJE91xxx3R/998881yzlVbGwEAAICymKQDAADbbdOmTdEEV/PmzTVu3LhY23PnnXdGf1955ZVVft6QIUPUpUsXSdKSJUuiib5UWr16tcaNG6c+ffqoadOmql+/vtq3b6+zzjpLb7zxhvm80jsUE4mE/vKXv0gq+XnxJZdcom7duqlhw4bq1KmThg8fXuWk2vfff18jRozQnnvuqezsbLVv317Dhg2L2nHxxRcrkUgoOzu7wnMvuOCCqD1z5syJth9xxBFKJBIaPnx4tG348OFR2UQioYceeih67MYbb4y2P/nkk2ZbH3rooajcbbfdZpZbuXKlrrjiCnXt2lUNGjRQixYtdPjhh2vKlClKJpNV2i9Syb69/PLL1a1bN+Xk5Khp06bq2bOn/vCHP+iHH36ocj2lrP0lSQsXLowe+8c//iGp5LM599xztddee6lBgwbq2rWrLr744uin5Nvy1FNP6ZhjjlGLFi3UoEEDde7cWRdeeKE+//zzKrc5mUzqkUce0THHHKNWrVopOztb7dq107Bhw/TKK69UKD9z5szoffTt29fc3zNmzIjKnXDCCVVuDwAAOz0HAACwnZ577jknyUlyY8aMSUmds2fPjuo8//zzq/y8goIC16BBAyfJ7bHHHq6wsDDode+9997oda+66qoqP69du3ZOksvOzjbLTJ48OWqb77+0tDR3xRVXuPz8/ArPzcvLi8rddNNN7oknnjDrql+/vps+fXql7b3xxhtdWlqa9/mJRMLdcccdbtSoUU6S22233So8//zzz4/Kz549O9o+YMAA8/2V/jd16tSo/A033BBtf+KJJ8z2Tp06NSp36623estMnz7dNWzY0HzdI4880l155ZXedpc1adIkl52dbdbTqlUrt3Dhwkr3b1X3l3POffrpp9FjDz/8sLvzzjtdenq697Vzc3Pd66+/br7OTz/95I488kiz7Q0aNCh3vA4YMMBbz6pVq9yBBx5Y6ed4wQUXuGQyWe5555xzTvT4XXfdVaHeDRs2uNatWztJbtddd3WrVq0K2o8AAOzMuJMOAABst5dffjn6+1e/+lWMLZEWLFigzZs3S5IOPPBAZWSELcHbt2/f6O+33norZe2aOnWqRo4cqc2bNysnJ0c333yz3n77bb333nuaMmWKOnfurGQyqTvuuEPHHntshednZmZGf5eGWzRu3Fh33323FixYoGeffVa//vWvJUl5eXk6++yz9dVXX3nbMnHiRF1//fVKJpPKycnRDTfcoFmzZunNN9/UH//4RzVq1EhXXnlldMdX2dfelscee0xffPFFubvdbrvtNn3xxRfRfyeffHKV66uq119/Xaeffro2bdokSTr33HP15JNPatGiRXr00UfVo0cPvfrqq5o0aVKl9TzyyCO68MILlZ+fr379+mnKlClasGCBnnvuOY0aNUqJRELffvutjj32WG3YsCElbS+7fx966CH9v//3/9S+fXs9+OCDWrRokZ544gkdccQRkqT169dr2LBh+vnnn711/frXv47uAO3UqZNuv/12vfXWW5o1a5auuOIKFRUV6be//W2l7SkoKNDAgQM1b948ZWZmasyYMXrppZf07rvv6p577lHXrl0lSffee2+5n4hL0l133aVWrVpJkq699toKd3Ved911WrlypSTpnnvu0e67717V3QQAwM4v7llCAACw4xs0aFB098yaNWtSUucvvZPukUceiZ53+eWXB7/u+vXro+d36NChys+r7E661atXu9zc3OjuoXfffbdCme+++8517949eu3HHnus3OPJZLLcXUxt2rRxK1eurFDmrLPOispceumlFV5n1apVrnHjxtEdYR988EGFMh9//LHbZZddonr23HPPCmUquzPMufJ3vpW9c25rqbiTrqioyO2///5OksvIyHAPP/xwhecXFRW5k046qdw+3LrdmzZtci1btnSS3BlnnOG9C7Ns/5owYYLZ3q1Vtr++/PLLcu3q0aOH27BhQ7kyW7ZscQMHDozK3HnnnRVe4+mnn44eP/jgg93atWsrlHnrrbfK3SXou5NuwoQJ0R2Z8+fPr/D4Tz/95Hr16uUkuRYtWri8vLxyj//jH/+I6j/22GOj7YsWLYruEBw2bFil+wsAgLqIO+kAAMB2W7NmjSQpkUioadOmZrl77rmn3NpkZf8799xzU9KWtWvXRn+3a9cu+PlNmjRRbm6uJFV5/a9tGT16tNavXy9JmjBhgg488MAKZVq0aKFnnnlGiURCknTFFVcoLy8verx0e6kJEyZEdyyVLXPddddF///SSy9VeJ2xY8dGd4D98Y9/VPfu3SuU6datW7l6tn7t2uaBBx7Q+++/L0k688wzddZZZ1Uok56erqlTpyonJ8esZ8KECVq1apV23XVX3XPPPd67MM8880z16tVLkqK1AbfX1vt38uTJatSoUbltmZmZGjt2bPT/W3+2BQUFuvzyyyWVvNeHH35Yu+66a4XXOuigg3TVVVeZbfnpp5/0xz/+UVJJX+ndu3eFMo0bN9b1118vqeTY33otwaFDh0afwcyZM/Xss8/KOaeLLrpIxcXFatGixTbvaAQAoC5ikg4AAGy3wsLC6O+4J3RKf+oqSVlZWb+ojtL307hx4+1uT15enmbMmCFJatWqlX7zm9+YZTt16qQhQ4ZIklatWmUGSWRmZuqkk07yPta5c2e1adNGkvTFF1+UW7y/oKBA06dPlyS1adNGZ555ptmWM844o5J3VbtMmzYt+nvMmDFmucaNG+u4444zH3/88ccllQSINGnSxCx36KGHSpKWL1/+i0Iktlb2mOnQoYN3YkySBgwYEJXdOgDizTffjH7efMopp6hjx47m61X22b7yyivR5HRlP4vt379/9PfChQsrPD5x4sTop6yXXnqp7r77br399tuSSn4m27x5c7NuAADqKibpAADAdiv9wu2ci+4Y8znzzDPLrU02e/bslLel7J18X375ZfDzv//+e23cuFGStOeee253exYtWqQtW7ZIKpncqVevXqXlBw8eHP29YMECb5m2bdsqPT3drKNly5aSStI5y34e7733XjSJ2atXr0rX62vZsuU221obOOf0zjvvSCr57Pfee+9Ky1troBUXF+uLL76QVLIunXXHZ9l0Xalkoi6VOnToYD6WkZERHWtl7xiVyq+fePDBB1f6GpWtA/fZZ59Ff3fs2NHcB2WPM98+yM3N1X333SdJWrFihS655BJJJRN/1bEmIQAAOwMm6QAAwHYr+7PLTz75xCy3yy67qGPHjtF/v+TnqNtS9g6dpUuXBj+/7MRe27Ztt7s9q1evjv5u3779NsuXLVP6M+KtbSsMo0GDBtHfRUVF3rZUdqdVqcomAmuL9evXq6CgQFLlE1ylrH23fPnyqJ4Qzrng51Smqp9t2c9Vkr777rvo723th8peo+wkXVVZ++D444/XaaedFv1/o0aN9Ne//jW4fgAA6oqwuDMAAACPww8/PFqX6l//+le5n8LVtD59+kR/v/fee3LOBf0Et+xP9/baa6/tbk/ZyZSqTOiUnSj6pT/XrUrd9evXT2ndqfBLJrxK71KUwlJot1Z2QvScc84ptyZfZVIxkZsK1bEfFi5cqF122WWbz7HW+cvPz9cHH3wQ/f/PP/+suXPn6vjjj//F7QMAYGfGJB0AANhuRx11VPT3lClT9L//+7/bvCOourRt21adOnWKflL7/PPPV3lSIJlMauLEiZKktLQ0bwBBqBYtWkR/b72OmE/ZMpX9LPGXKDvpuK2faRYXF5dba7C6lJ1Q+umnn8xyZdfWK2tbP7vcWn5+vnd72ZCFNWvWVOlOw9qkbPu3tR+sfbB1PUVFRdu1H8aNG6dPP/1UWVlZysnJ0bp163ThhRdqwIABKVnvEQCAnQ0/dwUAANutQ4cOGjp0qKSS9adSlXr5S40YMSL6++abb67y85555plokuyYY46p0s8nt6Vnz57Rz0bnz5+/zYmvuXPnRn8fdNBB2/36Ze21115KSysZ/i1atKjSsl999ZWKi4t/8WuVvXuxsjvkygY0lP057taWLVvm3Z6VlaUuXbpIKgnbWLlyZaXtKl13bmsdOnSI7gj76KOPUv4z1uq2//77R3/Pnz+/0rLWPpCk/fbbL/r7ww8//MXteffddzVhwgRJ0uWXX67x48dLklauXFlpuiwAAHUZk3QAACAlbrzxxujuuTFjxui1116LrS0XXnhhdIfV/Pnzdf/992/zOevWrdM111wjqeQuuiuuuCIlbWnSpEl0p+HKlSs1ZcoUs+zbb7+tl19+WZK0xx57lPvpbirUr18/Sg799NNP9eKLL5plS5NOf6lmzZpFf2/YsMEsVzboYdasWWa50v3ic8IJJ0gqmQwsnRjyWbFihV599VXvYxkZGVHy64oVK/Tggw+a9dRGgwYNUsOGDSWVBF9Y6xlKqrQPHn/88dFE7m233VbpXXeW/Px8nXvuuSouLlbr1q11zTXX6Nxzz1W/fv0kSffff7/mzJkTXC8AADs7JukAAEBK9OzZM5ogKSws1JAhQ3TnnXf+osX4t1fjxo01derU6G6uCy+8UPfcc49Z/ssvv9RRRx2lJUuWSJImTJigAQMGpKw9EydOVHZ2tiTp2muv1dtvv12hzPLly8vdAXjTTTdFkyWpdMstt0R/X3rppd4E3EWLFkV3Pv1SZcNEKruzq1+/ftHk0htvvFEupbTU7bffrnnz5kX/v/VdbmPGjIl+pnnvvfdq5syZFerIz8/Xb37zm0rvDrz11lujYIYrr7yy3PqEtV2zZs00ZswYSSXvdfjw4dq0aVOFcn//+98rnaTr1q2bLrjgAknSkiVLdOmll5Zb764qSn/mKknjx49Xw4YNlUgkNGnSJKWnp8s5p5EjRyovLy+oXgAAdnasSQcAAFLm4osvVn5+vsaOHavCwkKNHj1aEydO1KmnnqrDDjtMrVq1UnZ2ttauXasVK1ZoxowZVar3888/12OPPVZpmT322EOHH3549P9Dhw7VxIkTddlll6m4uFijRo3SQw89pN/97nfaf//9lZOTo6+//lovvPCCHnnkEW3atEmJREJjxozRpZdeul37YWudOnXSuHHjdM0112jdunU64ogjdNlll+mII45QgwYNNHfuXN15551au3atJOm3v/2thg8fntI2lBo4cKBOPPFEPffcc/riiy900EEH6ZJLLtGhhx6qzMxMvfzyy5owYYJyc3O1efNmFRUV/aIggq5du6ply5ZatWqVnnzySfXp00eDBg3Spk2b1LVrVzVq1EhSSVrp//zP/2jChAlyzumYY47RTTfdpIMPPljr16/X1KlT9fjjj+vwww+P7r7aeqKtSZMmuvnmm3XBBRcoLy9PJ5xwgi666CINHDhQbdu21Xvvvadbb71VS5Ys0YgRI8xJqj333FPjxo3TmDFj9MMPP6hv374aOXKkBg8erM6dO2vz5s369ttvtWjRIj3//POaP3/+doU0pNrll1+uKVOmaPny5XrppZd0yCGHaMSIEerXr582bNigp556Sg888IA6d+6s5s2b64033vDWc/PNN2vGjBlavXq17r//fs2ePVuXXnqp9t13X+Xm5mr16tX66quvNHPmTB1yyCHl7jp95513osn6QYMG6Ywzzoge69mzpy655BLdeeedWrJkiW644Qbddttt1btTAADYkTgAAIAUmzNnjuvTp4+TVOX/RowYUa6O2bNnBz3/xBNP9LZl5syZrn379tt8fv/+/d28efN+8Xtu166dk+Sys7PNMg8++KDbZZddzDYkEgl3xRVXuIKCAu/zS8t16dKl0rYMGDAgKrt69eoKj+fn57sRI0aY7WjZsqVbuHChy8rKcpJcnz59KtRx/vnnR+Vnz57tbcdf//pXb/2LFy8uV27z5s2uf//+ZnsGDx7s1q1bF/3/tdde6329+++/32VnZ5v1XH311e6DDz7YZrsfeugh16RJk232mblz51b6OVR1fy1fvjx67Oijj660nm31s5UrV1a6Lzt06OA+/vhjd/LJJztJbsCAAd56vv32W3fcccdtcx8MGjQoek5eXp7r0qWLk+QyMzPdZ599VqHejRs3urZt2zpJLj093S1YsGBbuw4AgDqDO+kAAEDKHXbYYZo3b57mzJmjmTNnau7cuVqxYoV++OEHFRUVqXHjxtptt9201157qVevXjrssMN0yCGHVEtbhgwZosWLF+uJJ57QP/7xDy1cuFCrV68ut9ZWWlqa2rRpo0WLFqlhw4bq2rVrueCDVBk+fLgGDRqk22+/Xa+++qq+/vprNWrUSPvuu6/2339/nXfeeeratWvKX3drWVlZeuCBBzR48GDdd999ev/991VQUKB27drppJNO0sUXX6zMzMzoZ45lEz9DjBo1Srm5ubr77rv14Ycfql69eurUqVN0F12p+vXra/bs2br33nv10EMPafHixdHncN555+m0005Tenq6cnJytHHjRv3444/e1xs5cqT69u2rW2+9VXPmzNG6devUsmVL9enTR6NGjdJhhx1mBlCUdfbZZ2vw4MH685//rNdee01ff/21Nm7cqJycHLVp00Y9evTQUUcdpV69ev2i/VKdWrVqpX//+9+666679Pjjj+uzzz5TRkaGOnTooFNPPVXnn3++dtlllygkw9KyZUs9//zzevLJJzVt2jQtXLhQ33//vdLT09W0aVPtvffeOvTQQ3XKKadEzxk3bpwWL14sSbr66qvLpQmXatiwoe655x4NHTpUxcXFGjFihBYsWBBbGjQAALVJwrkdLLoKAAAgBUqDE5588skKa6Y1adJEBxxwgHr37q0bbrghWk+uLnn77bd18MEHS5LGjh2rP/3pTzG3CAAAYOfGJB0AAKjzvvzySz355JOaMWOGFi1aFAUTdO3aVZ988knMrYvHlVdeqTvuuEOS9NJLL+noo4+OuUUAAAA7NybpAAAAyli1apVefPFFvfDCCzrooIN01VVXxd2klPv555+Vk5Nj/qT3X//6l371q1+puLhYe+21lxYvXlwtP/8FAADAfzFJBwAAUMeMGjVKL774os466yzts88+ateunerVq6cvv/xSzz33nB577DE551SvXj3NmjVL/fv3j7vJAAAAOz0m6QAAAOqYzp0764svvqi0zG677aapU6fyM1cAAIAawiQdAABAHZOXl6fp06drxowZWr58uVasWKFkMqmOHTuqY8eOGjRokH77299WSGEFAABA9WGSDgAAAAAAAIhZWtwNAAAAAAAAAOo6JukAAAAAAACAmDFJBwAAAAAAAMTs/wObcCguMOimGQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": { "image/png": { "height": 340, "width": 628 } }, "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.566370614359176\n", "Variance of Ylm = 1.0000000000000002\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.566370614359178\n", "Variance of Ylm = 1.0000000000000004\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 (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.11.5" } }, "nbformat": 4, "nbformat_minor": 4 }