{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# 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" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "pyshtools.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 = pyshtools.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": "iVBORw0KGgoAAAANSUhEUgAABOcAAAKZCAYAAADzmk86AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAjEAAAIxABZy/d1gAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdf5Bta17X98+zd+/uPvfODPPTERh+DL8NYaT4EX4IJkBEg0CJimQEDGCVEQIBsShGi6TKBIzRApmgEMHgYDAEAzJgBEw5ImplJJgwCIJg4hQDMwjMAHLn3HO6e+/95I/d594zd38+fda319q9u/u8X1W3euY5Tz9r7bXXetZ3r7PP82m9dwEAAAAAAAC4erN97wAAAAAAAADwuOLhHAAAAAAAALAnPJwDAAAAAAAA9oSHczdQa+2p1tpT+94PYBc4v/E44rzHbcb5jccR5z1uM85vPI52fd7zcA4AAAAAAADYEx7OAQAAAAAAAHvCw7mi1tpha+01rbWfba3db629pbX2ja21F+173wAAAAAAAHCz8HCuoLV2R9IbJP13kj5I0r+U9CJJXy7px1prL93j7gEAAAAAAOCG4eFczddL+gRJ/6+k9+69f6SkV0r6eUkfKOm1e9w3AAAAAAAA3DA8nBuotfb+kv7E+f99de/9rZLUe3+7pM87b//DrbWX7GP/AAAAAAAAcPPwcG64V0uaS3pT7/2fP/wHvfcfl/Qzkg4l/dE97BsAAAAAAABuIB7ODfdJ5z//XvjzN5z//L1XsC8AAAAAAAC4BQ72vQM3yCvOf/5U+PN/df7z3afecGttrs2ads80nbe/bOptAdcA5zceR5z3uM04v/E44rzHbcb5jcfRg/P+t0l68UPt/7r3vho7OA/nBmitNT37cO7todtvnP8sPZxrrT01pJukJ037r1a2BdwwnN94HHHe4zbj/MbjiPMetxnnNx5Hv/Kc/3+3tdYv+oXe+/MfNSj/rHWYF0p64vx/vyP0efBw7uXn33QDAAAAAAAALsQ354Z5+DilB28P2k8lrYcOPOQJamvtQyT97NAxAQAAAAAAsHMf1Xv/V4/udjEezg3zDklnkhaSXhL6vOj85y/33i/8SuMl/PrD/+fd/+xXaf6k+1euj4+27x24jKnPisu4Dvsg6eIv/e7Zdd437M81nnT6ddi367AP0vXZD4OpBc41PmWvz0l7DfaDugU3zjWeXK5F3SJdj2N0HfahiClHWt29q1/+83/p4aZfT30reDg3QO993Vr7ZUnvLemlodszD+d2sQsP/5/5k09q/vzn7WAzN8e1LpKS67DP12EfdM3fv+u8b9ifa1w8XYsi9zrsg3R99sO4Fu8Trh3uhwNcg/3gfcKNc43vOdfmfngd9uM67EPRtXn/rpdJZmLWnBvubec/XxX+/N8///nWK9gXAAAAAAAA3AI8nBvuh85/flr48085//n9V7AvAAAAAAAAuAV4ODfct2sT9PBhrbWPefgPWmu/S9L7S/pNSX9nD/sGAAAAAACAG4g15wbqvf9Sa+37JX2WpO9srf2Hvfe3tdZeKuk7zrt9e+/93s535qBv/rsJdrSbk0RujB2j+g/uTffSGiZX/JZPsr5K6fWNX8Dgyvd5CjfkUr52drXexQ7X0Sid4oWT2XYtvo4rXz9k5PZ2dSzzGOOHuPKxWRPmcnY1J4dxJ9ncruqZ1Dx6eyN/v6i8v1dYH+50LbsrPpfxCDfsPitNVBuYkzye99dlnysK2xu9b9elnqFu2dmzGB7O1XyJpI+Q9AGSfrG19hOSPkjS8yW9UdLX7HHfAAAAAAAAcMPwz1oLeu//VtJHSfofJf2qNiEQb5P0X0n6pCv51hwAAAAAAABuDb45V9R7f7ukLz7/DwAAAAAAALg0vjkHAAAAAAAA7AkP5wAAAAAAAIA94Z+13kB9sVZfrC/5y9Puy7PjXnG8yo5STlt6HXaMYvSaaY+HzfxBPMJTBPe4McrpZmafC8cij7vdVD7bdna+VHdk3PauxbjJbUnH2lXfqubinX1X+/piEpr/g8pL2dXxjOO6fZ5iH0r7Fg7odTlf7PauOhJzR+NeeYL2zatnYu0ydIxeO8i2e9iHUpr0ru6/sRaZIFV1ZP10I+uZ255qf42TKK+6nildqnF7I+uZ4r3MlkQ7TCgdvc8T1Ch533Y0Gd2WembE67j0s5hH4JtzAAAAAAAAwJ7wcA4AAAAAAADYEx7OAQAAAAAAAHvCwzkAAAAAAABgT3g4BwAAAAAAAOwJaa03UDtcqR2tLu60o7SxYqBXbT8qiVdjxw3tpddXTCK1aWox6GV7kB4iZSoprrW0sULyWhzDd7Ujh2MxTZLswHGlWupsYXtXngKbjB17hwlNV51CVtueuSaLCVsugLU0xhSveVbIXtvlez0rbK9wLGwy2SRJaIXJrJI6W96PwriVYXeWTljct+tSu9hB3AVc+f3CuKl/Cokfux/VFNhC4qtrbutdTi6uLdQzqQarTIel4zm8b22MHaXOJtQzjzS6ninuWyUx1Petba9Uz5ivAqU5y96Lwg0qH+PtwWPlU6ovwiCVusMdiynqi+Q610SVcV3XXSbGDtiPdvKIZzGXxDfnAAAAAAAAgD3h4RwAAAAAAACwJzycAwAAAAAAAPaEh3MAAAAAAADAnhAIcQMdHK00P1pe7pcLK5T2St9Jtjd84J4WEzbbiwsl24CGtMFCmEN4zd3sSIurjpq2yr5VVUIQJghusIswx75Xu2+jgzSCnY4xctxJ7OY0LI1bvhRKAQ1mbkl/vVVYbNcGI6T9SNur/DVbmrfs6xu/4G8p8CK8DrtYcuU1p9cxwetrNtiiFvzQCgs2+8Wy/fZKazBPsWr7rlZ+n+Aet6t6JtcXhX0o1Ci59nEbDOO6McK4PQ1i5pEYHmH6pq6ThC6411IJVwj9Y3iEPRZhe4V9m6R+miBQa2zdMcm0cI3rmSmCrHY2Rloxf2SYQ+yf7qnunC3UT3G+SPvmBtnlPc7tRwzlqgQ0FGqf0r6FmqFSXxTqp1Jww1XXMyPOi0eGc14S35wDAAAAAAAA9oSHcwAAAAAAAMCe8HAOAAAAAAAA2BMezgEAAAAAAAB7wsM5AAAAAAAAYE9Ia72Bjo/OdHB8dqnfLYUu2WSyWmqaCxBbr/0zYRs2lpJPU7JYYZ+7SxBL4xbSzWIaYiGxrJKEFlOMxobVVJJP5RPEWgiysYllaXtu3JiUG8aoJK+N7Ju2d23SYXdkirSx0riVxNBK2lhhjJjWGrj+bRau68r2TGpWn/uuOd3MNaaDP8GJWEiu9cmuhQTW8Jpj2phpT31dezxlZ34icUlms7Q9czxTMpkbN/cd/p5OEMw8iV3VMymBtTLG2tzj4ripJrJ9x9corvaJafCpzjFngUukl0Ltks63yo2kkH5bTng3tYutW1SsUSp9g7E1UaWGS2OU02+HjhvcxHqmlD5fSo4P2yslhu5mez3UM+4eHJPqXe2S6pZ0/brt+a5+3PKxGJmqWqjtUu1TqlFSTVSpL0o1it/ezNRElWTXfdczy5PLPYt5FL45BwAAAAAAAOwJD+cAAAAAAACAPeHhHAAAAAAAALAnPJwDAAAAAAAA9oRAiBvoBccnWtwJq30/wrqwyqnNCIgLKPv2lVnoOK47bPquittzizCvVv4ZtF2wOSzM7PuGfVuFY+zaQ2CC3KKqMSQg/UFlDPPbacHfGPLg+oYFlN1iy2ncQkBDZd92GkBhTvLSIswTLKp8ExdQHrsocg5+qCxSHMYY2VfyCx3nkIftphjmMN9+fetwYsTTovBXde7ciuNWFlBO7/Xc9HVtCosfz0MQQ3jNdpHisPixC26Yh+25cSX/sueh79y9vnCxzwqLO7u+ad+SNMZYY+uWNEasIwp1xyrUAa72ifsW6g5Xu7haRPL1SAzfcuFUoU5KC67b2iXc7921ngKuktFBA8WaYWZeS+pbCXmoBWcVxij0zWFfYQ4o1Chjw7CSa1PPjKxzKuFUMVxhitCqSj0TPm5WQiXsGKZuOe89qGkzRhihch466YQrHM84RqWeMe2VWmTTf/sCTn1d7VIJp5KkuWlPtY87xVPt47a377rl7OTk0r97Eb45BwAAAAAAAOwJD+cAAAAAAACAPeHhHAAAAAAAALAnPJwDAAAAAAAA9oSHcwAAAAAAAMCekNZ6A734+K6O7jybZrKOcTzDrEOuiUssSylmLpksjb1MyWSFtNazlY/occlpy9jXbG/lU1tWJrkrJqGFZMh15Vm4Se9KKUgtJVE6MUrJjFtJ/5LUlsPTzdwYs9R3adqq6Wam/yy812OTXVN77GtiparpsH7c4X1rKVaFvgqnXDoNK31NzFNOTw3JpSMTy1KaWkw3M+3rQt/KVD9L83qMDXY74bvaILRqup17LeF9cklmMYHV9I3pqaHd9U9juGSyg9Q3pJAdmHbXlsZIaWNzM5GkvpX2WWHCmMXJzKvUM66+SDVKqZ4J++D6p9rH1TmV2keSliZBNdVELtk1JtW7vuk6TQmsNu4xjOEmgXBaxDrHtRUSQ9Np6FJZpZAob2qRNHZOgR3eN9c5Zj4s1ETl2s6OMfzmEGuRSt9kR7XLFOnzPs00JRq7tuF1i6RSAuvaJKVWUllT/x6eKthTK54X5hgdFOqWMER+U4cnx/fKpBMSWF377GB4PZOST+dhe7aeSbXIfHsicQnxkrQwfSVfE1Vqn1RfVGqf2F6aMNy477oPJ/fvjRovbmcnowIAAAAAAAB4JB7OAQAAAAAAAHvCwzkAAAAAAABgT3g4BwAAAAAAAOwJD+cAAAAAAACAPSGt9QZ6yfFd3TkO0UyP4JLQYrqZiatJSWHLEDm4NNtLiWVLE/1zGpJWz8L2XIrr2dofK5fieroMya5z8zqWIdlVIfLI5BXFBFd3nGuhdyX2XY3Jp8PTzSoJrLOUhObGDce+knoW921ssqukZvrntFY3biEJbYJk10lUktBSKphJMovJZObAuQSyi8awCWlhDDfl9JCOtQ53Vjt26mve11ISWhDTnc3xHBkEfj5uah+ebuaSWWchvW1mEsQOQhJaSmB1KWQpgfXQ9E0pZodh0nEpZKnvgWlPya4HNq01vObQbtNaiwmsu1KpZ85s31SLVBJY/WR2atpTKqvrK/l6JtVEy9n22GezUM+49L2Qbrdufoy1mXV6nGhNWyVlPknBiTY9dXjdIoUaZYJEeVfnuHoh9U1j53qmUIvERFuzf5VjbxLpN/vmx7CuST1jT/FKony4ziq1yBQJrC4FtJIyL0nrA7O9dGswtUupnkkpt+F4umlrklOokj4f6hmXzDpLSasHhXom1AGHboxQoyzMGKmeWYRJx9VEB+Fit0n1YSJytY9LhpVyKuvU9cy94/uX/t2L8M05AAAAAAAAYE94OAcAAAAAAADsCQ/nAAAAAAAAgD3h4RwAAAAAAACwJwRC3EC//fC39OTRyYV9XJhDsgqri56ZVUDzYstpUeTtUyyFOZyYdvf7Ul4U2fU/WfkxTpbb7WlBTR8UUbt8ugmKSAvldhcIEBdyH74PYe1MuwJrW4cFlAuLCbuFkiW/0PHsLPU1CxoXxk1jVBZQrvRN24sBG+a9TuPa9y8tzBzOrZ0trJwWRXbnberrFjQ2iypv2l1b6JsWHi4soOymohRAkRa67ibEIF1P68X22GnZWjs1pHM2rdleOLf8wGEuS5OO22e3UL384tUu+EHyiyUvzILIkl8oedO+PZG4RY4l6Wg+vO/x3E9yLvzhKExmblHk1HdhTkS3IHLqm/rPU3jEjiaXVM+42qVSo1TrmRMzCaRACNc3BT/cXy1su6tzUj3j+s5DcNJpCHlwwi1V3RwjW7dIdoKKC7nHesb8QTrdCsFJMRBidD0zPOShXM+4GiWGR7h6Jhz7NIarUSaoZ+z7dJ3rltA/TCOlesZ9VcYFLqRxJWltQx6GB1yl7aWAK/depzHsuKHdHvp4voVj5M7DQl5NDrIK7aZ2cXVLanfBD5KvZ1x9smkfHtBwFMaw9UyYXFKd4/ofhQ9RLlQibW+KeibVLs6Qeubu4cXPYi6Lb84BAAAAAAAAe8LDOQAAAAAAAGBPeDgHAAAAAAAA7AkP5wAAAAAAAIA94eEcAAAAAAAAsCektd5ALz/8d3r+4eGlftelk63CM1qXWJZSzM5CnM+JiUnMaa3bfe+t/etMiWX3TOrZ082P4VJbZs2nprUYc+r4fVubVKGcRGljrPzmcrzZYG7ouLmQbubCeFLKqUtmraSbVZJdJWluxk59XWJZ7htSC006mUu22ozhkl1TWmshCS1JaWhjhfPQJZn1efh7Idc3JrBujxGTyUKC2Mz0dympkjRbmrnTTxc5gdUli8W3w/3B8AS5EGJWS1v2XWunUPorQJduFtJaZ6Z9HpLQDkxiWUoxS4lldxbbE4xLMZOkJw5Ot9pSKmtKLLsz3x4jJbC69uMwIbrEspRiltNat0+MeYjOS8lpY6VUVVe7rAvp87GeCe3319sXvKtbUvu9la9F0vniUlyfXvoxDmbb26tUBjH4NBx7l4Zo65bUnlIrE1ejhK52LivMe1KoO1J6qqsvUo3iap9QX6TtuXomJsdX+qY6x9UzofaxdU6sJbfHiFPIFdctiU1gTWOYc7ynBFbTN6XBu9pH8nWOS4hPfc3bISknotoyNb5PZoyYfGq2leqZEDztzqPSGRQnlzBKqZ7ZPtAulVXyyaxHi5S0GlLiD7YnneNCPXMn1TNhjDuz4fWMS2s9CrHRCzNGtZ5xtcuYuuWpw+3XOgW+OQcAAAAAAADsCQ/nAAAAAAAAgD3h4RwAAAAAAACwJzycAwAAAAAAAPaEQIgb6N0PflPvtrj4rVuFRZHX5nls6usWRT6NCyj7/bnf3QLKfhX1p034w521X4jy3syPcTA7Mm1+YcgDs7ByXhhyu29aKHm19sfzwC3WGvq6BZTD21QTF9s1bWFx2Ly4r2lLCyibhYfzYsvbbfPT2oLGrr8NYpA0PzWLFMfFlsOiyG6h47SAsmtfp3HdKrfhDUntIWxisOpC3jMT3JAWUHaLLYfFj92469B3Ftrd4szrZZg7TVBEWih5Fdr99DI85CFMOWpml9NCyfm6dhvc0SLckv2rwRZCHmZzE0pg2iRpYRZFrgQ/SNIds4CyWyhZkp40YQ5pAeUnD05s+xNmAeUU8nBk2o+b7+vGyAso+2N0aPrPwknkQpamUKlnKjVKJfhBqtUzboynw3lxd7ldt0jSwoZWhWtkZBhHui2sQ43iwm1s4I2kHoKhRov1jEvI8V1jUESpnjFtIRDChlOFMCzXN/WPgVqmRkm1Twy4qtQorm+on1zAVRo32lU9k2oU056CqGztE2uR7fZZCM5aL9J90oQ8HPox7G0g1fbxLSnUgmZ+SmFfrp5RuZ4Z2FaVsszMZzYXZCX58Ie5CYmQfJhVCn54YuFrFFe7pHrG1S7Pm/u6JdU5ru5wNY5Uq2dcjeLqk03fEAjh7qkj6pl/F2rLsfjmHAAAAAAAALAnPJwDAAAAAAAA9oSHcwAAAAAAAMCe8HAOAAAAAAAA2BMezgEAAAAAAAB7QlrrDfTb5k/pRfNno2tWhcSctUk9Ow0xOC7JLKWypoQ0l2729Nonkx2ZdpfgKkmL5sdwiWXzQorZOsQhuvbUdxXSmFyKawsJPT5+0b+OHtpbJUnJDBFTzCrpZiGlzbWntFaXWBZTzFKKq0ksc6msktRMX5tWJqmd+Rfo2lMSmpZmjJjWatpToqbrm1RSOVOKWSH1rJkUM0nqJp2sHfi5pc+32+cpCW3hx1gvzPZibKEZe4II5ZTA6pLMUvpXX02QTjg23awY4us26BLPJP+6S+lmKa3VpLJKPsns+SFp1SWwpnSzlFj2hOmf+j452+6b01q3x1jIz1k59Wz72M3DiTEzJ1fqm7h6xtUtqW+lRjkLtc/9UHe4euZuqGdc7bJY1ZLlXPr82FRWydcuqzARLVf+GK3MHJ6u31S7+L7Du8Z0yUI9E+dJV6Ok1FGX1hoSUV2qarmeMbXLPCXHu3om9g0Fm6ldWjiXXT1j6xbJ1znVlPkpahcnJKXaMVJfd42kesbULqluaeGadGn1qZ6xKa5hno2Xrzuc4Ri7oSufD6qfR8Yms6bPVfFgmJedPt+5esbVLZJ0aJJZj0PdchzSU10987yYPl+oZwp1jqtbNn2H1zOuRkkp8+meemjqH1e3JM+tZ1Jy7lh8cw4AAAAAAADYEx7OAQAAAAAAAHvCwzkAAAAAAABgT3g4BwAAAAAAAOwJgRA30AtnZ3rx7OJFCNPyhm4B5bQQ8Jl5dnu/EPwgSffX2+154Wm32KPvW11k2nGve7n2r2/pFpNe+307DYu1ugXMZ2FhfLewclo4flcLKMe+hYVZ4wLKNhAiLH5s2t2iypt2v3Mu/MEtqrxpNwsaF4IfJEln24uUNtMmyS+K7EIiJMkswtzTYss9tY+8dtKiyi7ERLJBESkQopmQB4UFlG37ItzSQsDGzFzvPcwBcVL1I9vWbo5dypRwU21PCyibfSsFP0ijF1COUmCNW0A5DDEzc+dBWIh3Yfouwv3yMKTQ3DELK9+Z+wWU3WLJz5/fD319u1ss2S2UnPo+2fy+uYWV0z31OLTPzIlxGE6uXf1tb7r0Ts3Fsw5nkatdXOiVJN0Px/NuNyEPYUHqQ9Me65lCeE8KolqaY7F0ITaSTmfbr3th2qR8nS1XZvH5ELDhrvW44Ho5WcZwgRAT1DNpAXs3vaSAK9deCX6QpPmJC3kohFNNUc+kGsW1h/CI7u7LMfjhetQzzYVWubpF8jXKWZglTe3SQv3cU11tAiTS4SwJx8LVLqmeWc9NwFV6fWtTMxaDH3YWcBVPF/M5JQTkuM+CKeBqYeZfFxIhScfz4fWMC36QdlfPxECIZgIhZiHwwtQzqW5ZhKrBBUTGgCvb+q56CnUciW/OAQAAAAAAAHvCwzkAAAAAAABgT3g4BwAAAAAAAOwJD+cAAAAAAACAPeHhHAAAAAAAALAnpLXeQO82a3phSDt8YBXSR1yuyGlIQToz7SkZ5X6IEXxaJmkm9J2ZvZvXIhJtGm1Kb3NJbSmt9dS0nzR/+RyE9JaZSYmJQVETBJaVmAS4aqpjJTHSpZulJDTXnlKsUuJrMymulcSydhqi10ICq01mDX370rS7NqmW1hoSSneWbpYSWF1CaUw3M4llyzDuwiREp2iy2O6b7fbM65iZBDJJ6vPh52cIjJQLg24mxWzTbuaW0DcmlhUSDkti6pnZ5zR3mtSzeSEJLaWbHYYUV5/W6hPEnphtJ3u6Niknlj1/tp16FtPNSmmt268vppuF98nlsC/CHDCfImnTSPWMq1H8uyQdm7ojpc8vUo1itpeS3iqJ8usQcbgy7Wdh7nSJ8icrX6O46+Ek3EfSdeauyXT92rTHCU6VUo1STmt1c2roa+qRlETp7gHlemZZqGdOXFproW6RfO1y5q80W49U6plUn1yTekYmrVWVesa0SVIzryMdC9c3qXzGmM1q9Uxbmc8NhUTjFMBrr7NKKqtkP9PsLpI+pM+nU8vsdPrceGAORiVlPrU/EdLnnzBprSmV9QWze34MU6O8wNQ4kk9gfSKmtW6fRIvwnqZ65tDV8eFmNKSeWYb741h8cw4AAAAAAADYEx7OAQAAAAAAAHvCwzkAAAAAAABgT3g4BwAAAAAAAOwJD+cAAAAAAACAPSGt9Qa6M1voiVlIBzq3Cmk+axNPeBTibk5MlM5pGHeukMbkHv+GhJ6VSQpah+fHpyFl7aRtZ8udNb9vC5P8chDS+1ySjksFlHwSjyTNTP+WooZc+5UnuPrmaRLSXBJaSIRyfUO6WaXdJbhu2s05MEG6WT8NOYImycwmuEohrTXEY6WE0hSRNZZL5JNK6WbN7FvvId1s8I5d0Nelky19b3f99pDWWklxzQmsw9pie/H6nSSZtaKQbubmyTT/urn6IETIHc39deb6H6cEMdN+FPq6FLPNGC7x1fd1yaxPhHvcsTkxjsMxPg7X78y8KQv563e+o5jxVM+cmfdpHfreN4XHzKTJb9oLF0OY9lxyvGuTpNMwx7lE+ZPQ956pfdL5fd+kuKa0wHSduWsyvv1XXLvY8mmH6fOlvi7ZtZDKuulvasnQt5n6oJTKKtlk1lijuHomjVtJn7/quiUkl9raZR6OvanBeqpzB7ZJOXPUJcq7Nklqpi4r19WVRGNbowxPlK9evztTSZ8vfBZMnxvd/BuTXcPnV/dZ9yjWDNvtxykNPtQ5LmnepbJKPpn1ybBvR+bYu/TVTV9/XbsE1lm4iQ+pZ+6kVOeR+OYcAAAAAAAAsCc8nAMAAAAAAAD2hIdzAAAAAAAAwJ7wcA4AAAAAAADYEwIhbqAjLXTUHhEIEVIX1mb1zLMeFpS3ayEOH1eSFmYRV7c4pSQtzOLMi7Aw5GEaw/SfhxVK3X7M0wKe5vXF4IfYvt0WAyGuswkCISoLvrq3OvdNC9qaBZTDwuFam75pkeLYbnZ6Ha4z1zeEPNjwh0pfSUqve6i4av/wcSsLHce+ZsHXFhZmdu+pJP/+hXPInhfr8PdbU4Sm7Oh6uokq68mX5up03zKTTuo7M/eXw3DfmqftmXvfYQgrWNjt+fPbhT8swvW7CAsoL0ytkRZQnpl3ap4WXA9WpmaIi2Wb7Z2F47YqXHwuOGsztjn2oX5y72l6/9P54s6tac5Zc40UV1a/6nyqXZkk4GrkXF0JzpLk72fpvuXucel+GNptSEMIhLDhD7FGMSFLsW5JQRG7qWdaCK2q1Ch27BQ0sTLzZFhsvoVAgO5qlHB8bE3s9kHT1OZjPx9cefDDBGIghD0txs/V6fOr/6ybPhdvX7+Vz9tpe/lzv3tGYLva8IcU/JBCq3w94zc4pHY52tF33PjmHAAAAAAAALAnPJwDAAAAAAAA9oSHcwAAAAAAAMCe8HAOAAAAAAAA2BMezgEAAAAAAAB7QlrrLZVSRtYmWWwe0orOTPBL9WmuS4+Zj01XuoBNtgmpNPb3K+lm1z0O8bbEqTkpxaxiivMwjNFde9hn3zelqW337VMci4gwMfoAACAASURBVIp43FKa2vaJ2E2bVDxlXXpbOm5xjN0cu5gEfFu49++K58PrMv+mBE7bN9yLXHu6b7ntVe7LLuG0aky62aPYMUJSY+Vq98muKcU3jTG8vnDvaXr/k8q5tStXfp3leO4r3Y2rtNP7RWXsdP9MSamDdyHsgxt3V6msRamusgmc4bp2rzsl8Nq+O6xRS666xrxq13hqqXz2rCRup/vWrub7dC9LCbNOqc4JzzV83+v3PbXrt0cAAAAAAADAY4KHcwAAAAAAAMCe8HAOAAAAAAAA2BMezgEAAAAAAAB7QiDEDbTSWqtHLI67Dosvrs0SymcmJEKS1mYh0bOwuOgqrOno9nMVVt9c9+1nxSvTthnDt5/1+aBx0xjruG/b7a7tIle+pupVb69yOEzfeDhd+yx0Tn/d4PqnBUNNew/ba3O/wbY0Y6S+K9d3+zze/IFZNFi+b7isVVtG3QiLp7b0npj+LR37mRk7LdbqjpH7/c0Gh7enc8j1DeP2uL0w9lhXvaCxW8T3iuebNP+mObzSN91ffF9zj0v3nEK7u5dJ0qptX9hn4XUszJuyiiEIod3d880+SL5mSOERiatdXN0iSSuzvfT6XHuaItPxdO91fJ8q9UysUYYfu0o941Rqn51K84hdOHyH+3aFc3W+X1ztfSveP929dhXqCDNGut93Vxuk0icV0CPDKmI9k2owV+eEvvZ1pzqw0DfVo7X3unAOxe1tN5Xq+OSq65krr10qfYd/9qzUM+lzsbufpXoo3bfsZ/lwHvpnBP4AuQCnFDRxFu7uC7MbLvxSGla7rMZ+ngr45hwAAAAAAACwJzycAwAAAAAAAPaEh3MAAAAAAADAnvBwDgAAAAAAANgTHs4BAAAAAAAAe0Ja6w100pc6eSilLCWZOWcm2Sgmo5gktPsh1eR+SG15um+fYvdN26Z9YfYh9F1v9930306aSX2Xa5OQt/apS0vz+lJijhtXkrrp79o2f3DVcUVGNXXJtceX51KlQt/59kmX0qoq7f0g/N2Eef/aKiR3pcSyg+3zNh02m0Fn0l4l+TS1lU8lSkmy3aUvpviolNLltpdSwUySWUxCc+3mWEpSc+0HIbk2tNv+IbHOnS/5fPObs+lmoa9rj30LqWnXYWpJQhi4vUbS3Onm5WWY11Pquet/svbn4Ym5vxy3M9vX3eMkadGXW23zkEK4yDHM22Noe9yU5H4c6gjXP6WYzXcUs5dqFJsoH/ret+nzfnv3QwKra0/vqWtPfVPiqzu30nnoztnK+Z3qmXSd+WvSdr0WqvOhnVMnmKtr94Bw7M29qC3DIOYeZ2sASW3t5wC3Fyk91e5xqg1c7RLqmdLJVahnSsnxUqhn0rE3dWCoZ2ydU0nPVahRQp1ra+JqXe12o5JgP8FniStPdk0Kn+9ce/rcaOuZlLSaPr+asdPn4qPZdu2S+lbqHFfjZKmvm59Cknt8rrHdfxESm53Zc77TdlKoxyr45hwAAAAAAACwJzycAwAAAAAAAPaEh3MAAAAAAADAnvBwDgAAAAAAANgTHs4BAAAAAAAAe0Ja6w301PpMh+uLE0LSn7pck5RYdmoSYc7C89ycbrad2vL0+sj2vWvan1of275pjKdX2+2VlL2TQtpNSsZJCWkuoWe9Tn1do+06jWYzQ23XSdLNzKELh1Ntvj3wzCS4StJ64TfYTHpXDPSyaaa+b2ITwEKil0su7SGNq9l0s7BzIe3RHYuSlOCaEo9c/1l4fS49NSW7Lrav65jKavpKUj80YyxCEpp5/9L5lttNopc5v6ViWqs9xr5vTlAO/cdKp5ud40ICq5l/VyHdzLWfhsnlNN4btq+zdB9xSWaL5u9PMzvPSvPqBPMcZyE17cxck8fNVwf3W0iHNftmTuMoXJFRJfvM1S5n4QS3Se6xbhmeKO/qFsnXKJW+kj+3cj2z3Z7Ob3c9VK4nyV+TuTgY2DYRuxvpnC3MqX3md9rO1WFed/cAlxAvSeGStPfwdN+yvx/a01vSTGpoTDl16bDLkCjvPsekFNhQz8RCbqhqPePquFSvFRLl3Rg91S2x3aQwp3rG1CipbgnTpNbmHE91vBtjivT5KNxrR0uJ8u5jQ/x8t92ePje6z5kp2TV9fvWfdX0iqvsMPQ8vep4mqILV7GS7Ld3DTXVwP8wLh6meMefFLIwxZEZ9KiRcj8U35wAAAAAAAIA94eEcAAAAAAAAsCc8nAMAAAAAAAD2hIdzAAAAAAAAwJ4QCHEDvX091zqtuvkIa7Po5Gl4RusWUD4rLJSc2t0ix1IxEMIsWilJ7zTtrk2S7q3Mvpm21H62CoEQYbFOt4inWxh08wcD2y5qL6iEOeQFlE1bXEh2u60d+GPh1vVcHaadSIsGb/dPV9C6bZ/js7T4cVj02S76GxZFlgl5aMswNbvwB5seckH7jhYwlVk0OranRZhtOEYY1y08HcMcCosiH/q+rn0dzsMezuW1aQ9rttu+6Xpy12RlsWVJdpXwnYVESOpmsWR3f5L8wsrLVbhvFcJ77i/9fH9gJp1ZWGA4hTw4aaFje18uhCw9aRZVlqRjExSxCIEQhyGKwfVPC0RXjkVFOi/c8XR1iySdmhk/9U31jKtRUj3jQh6uvJ4J57e7Htx1I+XrzF2T6X1y1/oUppjLKnVOrGfMXJ2Cl5o5FuGS1DrunNm38PrcrTaFObSDsL0zs4Mh2KCZGsUGWUk+zCHVJ1ddt6SarxD2ZQOqUsCVq2fC+5HqHFejuBpHklau9gl1S6q3fT0T7nEuECKFR7jNFT53SBPULtXPW6Y9fb5znwXT50b3OfP+LNQtsxSCMLyeqUj1zMqFSRbuta5uSe2L5oMtcp3jjsXl65Zff0Q452XxzTkAAAAAAABgT3g4BwAAAAAAAOwJD+cAAAAAAACAPeHhHAAAAAAAALAnPJwDAAAAAAAA9oS01hvoV5fP02lKczy3Ds9dXbqKS1aRfDJrTEKLSSyHW20nMd1su69rk6SnV779nmm/G/reXZrtmTZJumdSz05CWuvZMqS4mtSzdUhCs+lmU0QnVtLNqulINo0ppLcdbKfjpDAuFyrU4rFIyVvD921mkrdSYmgLqVmzpdlp16aQZLYqJK2GVLhyiutYMdHWHLvQt7vktJRYZt6/nNYa0sYKiWW272HqG1KsTHvcnpniK+nHMeivklpYmXJiilnaoHn/QqqjmzuX4T1198Z5SOOaTRB17RIV033yJETznsy37y9PzHwq59FsO7HsqZRuZvqmdLPDkG7m+se0VpOENoVKPZMS5V2NEtPnQ43ikuVK9UyoRVKdc3e5fQ64VFZJeueZ6RvSWk/MNZJqSpcsKPlrMqay2jj4QlJ9MsFcluZJmygf01pN33gsfLPtGpMozbFPfc39c3aWapGQwGrSWluoUbqpc1oq7lwt4hLpU98pVFJZJR9/GxJfbY0S01pdLZJSYAv1TKhFXKJ8TGUNY7jzPqXPu6k2fj4o1DOlOmeK0OhCInX6fOfmzvS58cQksLbm5/VK6ug6JscX6plw/zwxabJPF+qZ43Zq+7oaJaWypjpnbj5QjqlnfnPptzMW35wDAAAAAAAA9oSHcwAAAAAAAMCe8HAOAAAAAAAA2BMezgEAAAAAAAB7QiDEDfS25Yv11JlfEPKBVeG5a1rssbIwZGp3CyvHRbFNe1psOS2K7Nrvr8Kiz6ZvWkD5vlks+eTMj3sWFtW1gRBh0WC7sHIKTKgERRQCIaoLKLu1QdM6wH6B6LSYqVlUN+xDWCPeBgLMlmFBYxdWERbKTQsrr01IQwvbay7QIYQ8NLNYclwD9roEQpj2uBCwWWzZtUmSzBjrcAK4xY83/c32YiDEuDAHKQRChKnctefFnU1bWHc6tdsTaYoFlBNz6aQF5d3CysuwgPKpeR2tsFDyZtfM+xQOxqlZOf7ErSYv6Z4JfpCkd662F0u+M6+EPPhFkd1iy3kBZd8+N29UWnja9Z1CqmfW5j6S+roaJQd3+PfJ9U/hEa4WSbVPqmdcjeLapBBaFUIeXD1zGq6ndJ25azIGQuzmtLigntk+P+M9J+UPuHom1lrD55fu7pOx1koL5pvXl+5Fpu6YhXom1SgzV9OGQAhXz9gaRyFU4hrXLVIoXVPhWahn3PuX+rrgh81+uL6pRhlez7i6ZdPfbc92DdvzfW09Uwy4GhsIkT5X9XAu20CIQsDV2SyFsUwRWrW9Hyno52S+ffDvhLrl3twHGb3T1B1HMx+c4NpdjSP5GiXVLTOXJHhBf2dIPfPU0u/rWHxzDgAAAAAAANgTHs4BAAAAAAAAe8LDOQAAAAAAAGBPeDgHAAAAAAAA7AkP5wAAAAAAAIA9Ia31Bvql0xfpydNn091WMcZmWyXdzPVN6WbLkE53ZvYt9XVJZi4JT5JOQwKrTc4LfV2S2WlIWj0ziWUplfXsbHhaa1+F+CDXnpLCJgix8glUw1OJpBAwenGo8LtK6Vgzl7zmh0iJbC4oqIVjPzMBRDZVTNIsjNHc9lJimQsECslkbj9yWmtod+MWktBs2pyUE+dce0osc+dhSms101ZKG8upd2aMcA65ZLF0vqUUMptulva5lIRm2tI1khKNTXs8KyoprinwyswjaT5cm3Nu1UIqZ9t+4el1rMJcdmaSzNJ879LNng4pZschgfVwtp0gdmDaJOnItC9C34NCutnBzL9RlbRW+/shNS3ZVT3j0uliPZNS6QvJvK7OSfVMSmB151ysZ0zflMDqxk2prKtlqA/NtRrrmUL6fOTuDaGrncvCaZXmyVpdZY5FTGAdniQbAg7lLuFZSNScuZohBBamGmXmzoHw/tm01nQsK31D+yS1ix04jWHawvtnk+rj/dclx/u+qUZx+xHrAFf7FNPnbU1UGCPWM2bfcs0fPqeMTZqP5+Hwz2yp7lzNxn0vqod9WIUEVnsfOQhJ3vPtSedp96FI0qHpK/l6JqW1ujpnEWoG3zeltQ5PlB9Tz9w9PRn8uxV8cw4AAAAAAADYEx7OAQAAAAAAAHvCwzkAAAAAAABgT3g4BwAAAAAAAOwJD+cAAAAAAACAPSGt9QZ6y70X6fjpOxf2WY9MPJOktYkrcolnF43hUs/SGC5pJqWbxVQa074MKXtLk56a9s0lma1SsmDY3nrp0s1ShFglmiy0u64pVKqQbpZSk9x+xEA2m3gVupoEqphiFVJVXTJrSixzwYctpDqmFFcXNpQCgUpprYVxo7HpvsUULHvOVZLQUlqr7eu7xnSzked9Omdzmpppi4lsI/tW083c666812FySednd+d9mFO7maDC5WuTzNbh+l3Nh6ebnc79Fu+btNZ5SD5dzI58uxk7pZy6VNWDQrpZSiaL7WbCyGNUIziHSfWMqztc3RL7prolJbDa9PlQz5i+MeG9UBOdhZrBje0S4lP7KqS1ulRWKdQuhbTWFuNMJ0jUdPN6mPfS9uyZXNjlWSFRPqVhuqRVqVbPuNolXaZxe7ZGSdsbl8B6E+uZ+HGrktZaGDelw7r+lbTWXDOEa6RQo9h02Cn2LR17NweU6pliu02fD/ctO254T832Ut2S7kVnZjI6Wfo36sDUIouU5B7aXVpr6uvrmeEJrCll3tUtaYxkSD1z/969weNV8M05AAAAAAAAYE94OAcAAAAAAADsCQ/nAAAAAAAAgD3h4RwAAAAAAACwJwRC3EBvvftCHS6euNTvusWyk8oCymmJRbegcRpjZRbUXMfFlsMCpW57aTFws1hnD33tuIXFjyWpu/6VBZTDuKUFccNimN0sXKu02H3YYDNjuKAJyS/i2lIghJml1oXFjzftpi2O4dpqi0mXghvcosjpvR74+xdub0dKi+1W1gJP41YWUI6LMLuVvH1fu3hxdXtuoePCGJVFkWPwQwp0GRkIUQp+SGMv/cHoZlXkNFf39fYG0z1gOQuLBptFhk9nvlyaz7f7ztK44SC5xZLnYYxmxkiLLbvtpbe0EhTh9mEfbPhHqg1Me1zbO9YopmYo1DMpyCrusxsjLTJeqJ/ctWPrE10QWlWqZwa2qXjfiifzdlOc9+IN1BzPcChsoFZ6fQvTFmqRNG+5Y+eCrFLfFgOnhtczecF8F/6R+oZ2tw/XpZ6plGaFe2otaCJtrxBAMbIWmWKMHEDh2mr1jK8lK+l5vjldI3Z7yzC0ORjrcE26z5nrcCxWoX46NbWEq1skX7u4ekiS5qlmKNQzrr5I9Yw7xNWAKzvuiMnl9K4P+RqLb84BAAAAAAAAe8LDOQAAAAAAAGBPeDgHAAAAAAAA7AkP5wAAAAAAAIA94eEcAAAAAAAAsCektd5A73jnkzqYP3mp3w2BMKGvSf4pRTL67aUxbIJYMR3WJvjFFEGzvZiIatpTCmFM2Bq+bzYRKG2v8p7ElKfhQ+QxTKpfSmC1iajDk8lseph0wfEM/Qf2bekYVxJYd9S3bFepZxOktTqV1LRywlol8bWQJBv3wyWWVdLNQuJVadzSvqV05zCG7eyb3fXe48k5fD7sq+1JJx+3kJDmJq6UnurmvXB8UiqYbY99K+MO75tU+qfXPdbYuqXaN23P10Rp8EI9U9iPWKPY+insm60vCkmrUqkmsvf24n3IHqJ0XZv5IvZMCawuEbHw+lLNYdurybWFmsEnrdbqGTv2FCnxu+pbMcWcNUH6/Nik+jRGpbavp8OaxljPmMT1SjpssZ6pJOWWBl6He62bcwrzb18Nr1FygLb/A1ejpKR6d4zc74eum3aX8F5In6/UEaWaKo4xfHvPtXznbooevjkHAAAAAAAA7AkP5wAAAAAAAIA94eEcAAAAAAAAsCeDH861llaPGvS7L7js7wIAAAAAAAC3VSUQ4sdaa5/fe//ZygZaa58p6ZslvaK0Z4ie/q0jzXW8+w1VVv3e2SKwxYX4K33Hvr64IG5hnysL/4dxKwvwpl1zi2f2WfHYFxajtSPHVbjN708RmFBa/Li26OfOFkXexe/vw44WYS6+TaU3alcLNpfCKkJXO0Y1rMKGEqQNuoHDsGGheXe5x+AWu7hzGNe+jjTx+ebaLapw3EqhKcMXbM5jFPpW9+O6Kk8CbowJ+k5RP01wLxo+ru9aqmcK4RGT3MML4VRxDkhlh5tz4iXpkjt830nqmcoY9vcnqGeuczjVLo29/KYIw4pjmzp+knuRb75p9UwpO68QJLj5AxNskPra11epZ/ywlXrGjnvR2JW+lcloV/XMFdUtq99a7mTcyj9r/QhJ/09r7ataS/lGz2qtvay19r9K+j5J737ZHQQAAAAAAABuq+qac0eS/oKkf9pa+8DUqbX2+ZJ+VtJna/Os8yb+/QgAAAAAAACwU5WHcx8r6U3aPGz7GEk/2Vr78tae/S5ma+29W2s/JOl1kl583vdNkj5usj0GAAAAAAAAbonBD+d67/+XpI+W9Kck3ZV0LOkbJP1Ia+0DWmtfKumnJX2qNg/l3inpKyV9dO/9x6fecQAAAAAAAOCmK/2z1t77qvf+WkkfLOm7tHkI94mSfk7SayU977ztb0v6Hb33b+y9r6bdZQAAAAAAAOB2qKS1PqP3/suSPre19r2SvkebNeWapLdJ+mO993843S7iuWZPLTTri33vxu30mK6O6FNV/cG4Dofo6gMEr8Orxk03RdDijVNJRU5xXIW/4mtTRJQ+ju8Tbp9bftuqJVFSzzzrOrxq3HS3up4p1S3SrooG6pnrrb9zN89iqoEQkqTW2qy19iWSvk3PPpiTNqmsX9haI50VAAAAAAAAeITyw7nW2kdK+nFJ3yTpRZLWkv6ipH+kzUO6Pyrp51trr2mt8fUuAAAAAAAAIBj8cK619kRr7esl/TNJH67Ng7iflvSxvffXSPoUSV8m6WlJT0r6Oklvaq397sn3GgAAAAAAALgFKt+c+xlJXyFprs235b5W0kf13v9vSeobf1XSqyT9qDYP7z5EmzTX10250wAAAAAAAMBtUAmEeO/znz8l6Qt67z/hOvXe39xa+2RJXyLpL2jzLbrPl/QFI/YTDzl4qulgvfsVHidZ7HPsGOH3S/s2Qd/K2uWVsfOCxmaLuzwWY/ftorGH9o2vr7B4cWmB6AlMsW+l7U0wxnU2dp3q6u+PneTK29tuapV9SNtz7bHv8O3F07uwvdoYV7xvg7eWx6icA6WF33fVN7j6Remv1pXXM5XLekf38Li9Ccat7fME9UypZhiyUw/GKEwOUxz7wrGIRtfVxYt9R3X8rTHF3Fma73d0QCe5xxU+RE1xjwvbs6f4Lu/hhXqmtG9Dt6ULLrMd1RKlY1Fxg+qZdnc312Llm3NLSX9Om2/L2QdzD5hv0QEAAAAAAAB4jso35/6D3vubKoM/9C26L67tFgAAAAAAAHD7DX44V30w99DvdUnffJnfBQAAAAAAAG6zyjfnntFae7mkPyHpIyW98rz503rvb51qxwAAAAAAAIDbrvRwrrW2kPSNkv64pMWDZm2W7zsy/V8r6dMkffZlv3kHAAAAAAAA3FbVb859j6RP17OhIE9Jev4F/V8m6f0k/SFJPJybyOKppoPVJRNCrjo9dUeJZT1FmVQSvdwYse929Evc37BvbnstHs/tP0ivuYVYmrGvL8bFVMaoJL6m7Y0dV1IrjGGPZ+w7fD9i11Li6/C+cd+uWB+ZQtYLE0PcVBjD7lsYpLuE7Mq4kmTGqPQtJZatfdcW2t3Y9jUrhC8Wk7uaGzvtcyWB1YyRX3N4fa5/JY228jouGHvwGNWk3JH7MEk64RSuOj21kmpe2IdKLVHat3gPN/VMoW6R/P0lvw5Xz4R7dRjCjZ23V0iHDfth65FKsmsat9C3VF8UtlcZNwwR+frp8atb8rihZiiMG/fN3cML6anpfh9fX2UMe08dPq6tF6R4MLqtAyrb810rNUPlfh+3Vxh3mtc3sO2iMVzfKdJox6bcFrc3RHt6z2mtrbU/Iukzzv/vD0r6nZJe+Ihf+2Ft5vWPudTe7Uhr7UNaa/9Da+2fttZ+7fy/N7TW/nhr27eGan8AAAAAAABgiMo3577o/OcbJX1W7/1Mkh7xbOqXzn++R33XdqO19u9J+pfn//e+pJ+W9ApJn3z+3+9prb36PMii3B8AAAAAAAAYavA35yR9mDZf/vv2Bw/mBviV85/vWdqr3fpZST8q6bMkvVvv/aO1eXj4X5z/+edI+gMj+gMAAAAAAACDVB7OveT8Z2XtuJef/9wKi9iX82+4fVLv/fW999MHbb33b5b0I+fdPvOy/QEAAAAAAIChKg/nfvX853sVfud3nv/8lQt7XbEL/gnqT5z/fJ8x/QEAAAAAAIAhKmvO/ZikPyjpUyV9/6M6t9ZeIOlPa/NPYd94qb27eh9w/vMXdtR/EodPdR0sL7nE3djEskK6aOpfSVqtJoiNHSP3NQOnvvM0hklIC30rj83jGK5vKemtkGIm+cSxeRrDpI2Fvi6FLCWTzQpprbM0xmw7giid3q6v5C+Teepr9tmdbkkp7VX5GI21niCxbF1IG3PtaR9Sux0jJF659tS3r/1F4rbnUsUkqa+2x4hJaCvTVgxes0lf8e+lhr/XMUGskKrazOtL4/q+ad98u0vvmrljLNVeR0zmHTlGJRUuvadXnfhasbOk1eF1Sx4j7IebAlLNMMEYa1cHxHGHp8FX6hmlvnbXwtyZagYnnZy2thtei6QxUj1ja5R5uN+bcdM9vIWawdUuuZ4x+xa2l2oD1x73uVDPVGqX61C3SLX0+Eo9szI1Q7xdxPpieF9Xu6T6ItZPhTH6yqWzp3rGjJsOe7xPDh+jlEZcqCV2Vc/E1zzFvlWS6nc1RqxF3AWV+vp2a8TUMru3m3mp8s25b9fmHP6TrbXPuKhja+2lkv5nPRsE8V2X272r01p7laT/+Pz/fsfU/QEAAAAAAIDnGvzNud77D7XWXq9N+MHrW2t/Q9LfeqjLh7bWPlKb4IgvlvRibZ5H/uPe+9+dcJ8n0zZRs+8h6TMk/TlJT0j6y9oEQIzuP3AfnhrS7bLjAwAAAAAAYCfe3B7xleDe+/MfNUjln7VK0udK+gFJnyLpC8//e7ATr3+o34OHST8n6Y8Ut3ElWmtfoc2DtQfuSvp8SX/LrTFX7Q8AAAAAAAA8SunhXO/9XmvtUyV9qaSvkvSK0PVU0v8k6TW99yHfDNuHX5T0DyW9TNL7S3pS0tdqs2qG+2eq1f6DDHmC2lp7mZ4N5AAAAAAAAMD+vbL3/mtjB6l+c+5Bcuk3tda+RdLHS/o4bf6p5x1J75D0k5L+Qe/97WN3bpd6798r6XslqbX2pDbfAvwGSa9rrb2k9/4NY/rv0uKdXYuHAyFG/qPXKRZQjovt7iygIbTPxy10HBc/Nu12AWYpLtYps29xIcoDt5K375oXqLWrgfuuhYWnY8iDaY+LIpu+87iA8vC+qf3ALKycwhzmZnuLuV8Nfl5YQDkFQti+YaXVyqLIU/StLpY8dIxSQEM48ZdmoWO3qPJF21uZ9uXKX9grs2Cv2wcpL8K8MiEPK7dQsqS1uUbsosqS+swsJr0sLnRt2mKYg2suhBJIfiqyix9LmpnXnfq2pWkrLMy82Z5pDGO4vm0VFkuvLKA8xeLOhfCISlDEjtZmLxsdWhVeSCWgIYZT7aieSXWHmQJiX1v7pH1Ip4Xbt3hyDR83ziOV99rVo6lGTSEPpgar1DOxvijUPmkMV8+4ukWSDkztUqlbNmNvb8/tgyTNzBtYCaCYIvhhinpmir4rc1HFvqZmqNQtknRmahdXt0i+Rkn1jKtbUnsKhFi5GiXVMy41JQUyxUnAtRUml3CM0+lp78uhBrOBEKmeMeOmcKo4xsjtTVHP7Co8YqcBVwPGbfd3UwyVH8490HtfSvrH5//daL33u5L+SmvtvqRvk/S1rbXv6L2/Y4r+AAAAAAAAgFNJa30c/G/a/H31HUkfsYP+AAAAAAAAXYNLiwAAIABJREFUwDN4OPeu7kp68I9ihnyrsNofAAAAAAAAeMYzD5Raa18u6csnGredj30k6VjSUe/9aKKxd+nTJR2e/+8f30F/AAAAAAAA4BkPf9vrhZLeV3k5Vml71cShqytek+WDpdbaZ0l6iaS/3Xv/rfO2JunTJP21825/9UGgRbU/AAAAAAAAMNTDD+eWkk4G/M6Db4q18/9Wevafdi707D+V7ZJ+RNJbJL1z9J5O5yslfYKkb2ut/YKkfyvpg7V5OClJPyjpq0b037nDu2stlikSdCMmiBUC/PrMJJ+mdLO4veFpYz4JLaT5hDFswmElgfWgsL3wD5l7aLfvWAoacp3DMU7pM5Un4vZ9jcmuIbnHJI7NXOqspJlJEDs4GJ7AmtJTD0Lqmeu/CGljvm/YXkpTM1FIua85bjHdzKS0pb6FM6CSkFZNcHVpqyltbGkmh0paq/t9SVqGCcqlm1X6noV0s9Nl2I+ZSVOb+zGWy+33ZCU/7trMLj3GLxbSzcaH9daSt0LSm00bM6msUkhPDX1jWqtLfA2JZT6ttbg9M3YlYbatw5zspvVUOlRSXAs3l7i9tBuVxVYKifLuYJQvEVcTFeqZtUtsv2gMm4Ce+pqxU41i3pNS3RK0dOzN+TlBEHien1x7MX3e1jOh7/ygUM+45NNQzxyacSVfu8SayFyAlb6SdGAmuUpaq/t9ydculbpF2l3tkuoOV7usw0TitpfqC1vPhIkhpaqemfrH1S2bvttjpL7LkNZq02FDXzf/rsNkVqpnwgTl3upWKGjiaZXaKwml5nKICaymFnH1SRp309/cRwo1Sjkd1sz36VjYsUMCq09rHb4PUaFG3epzUixwBnrmbO+9f13v/c5F/0n6Am0e4N2T9BpJr5B0+NCfH2rz4OovafPA7oO1+VbZf7mTvb+cz5X0dZJ+VNJc0qsk/ZKk75X0+3vvv7/3fm9EfwAAAAAAAGCQwSEGrbUPk/Q3tXno9om99594bp/e+1rSv5b01a21H5b0w5Je31r78N77Oyba51F672+R9DW76g8AAAAAAAAMVfkHBF+uzT9b/cvuwdxz9d5/RNI3S3pPSV98ud0DAAAAAAAAbq/Kw7lP0uZf5v4fhd95/fnPP1z4HQAAAAAAAOCxMPiftUp6j/OfldXvTs9/vrLwO3iE+b2VDlZhJcZzpcV2U+iCDYRIY4T9MGPE0AUXCBGDH3x7M2PHMdzZHxai7KWFZD27VnJaUNO9J2GBy7ELaMf2tK5rXEDZLO6bFh42iyUfpMWPzRhpoeSjsL2j+fYKqmlR5EOz2uphoa8kLcxKomlR5ErfuTm7cnhEWLy6NIUPtwonjFsUOS3MfGYu1rQoslv8OC22fLLyt7rlwXb/09D3/my73e2DlEM6Tk2YylkIj6joJijCLuQvqaeFct1pkRZ4L4gjXIcFlFOohFlAOY0xNjwi9a/0nWJh5rgIs2suBFBMIdcdhRrFnMsxECIGUbnAKd/XTQ0mDyb2TWO3FNzgwpfCwXBvdTotUnnhsltaOp7uWk+vOWzPiu+1GSUGWQ1vd8EPkq9nUn2xMGOk+uLoINQXZnI4NjXOZuzt9hT84Oqk1D+GZJnJaBEmKFejXOe6JbWnMVxA1VkMhBje99R+eJFOTUDD6dz3dWEOJ6bG2bT7i7WZ9+osXdhGuKWqu5ov1jMphMa0VcqZYkiAC7PK92XTVqgvcj0T6o6RoRKuHpKKNUoKhHBhFYUatRwIUQiKGFLPrE8vfhZzWZWP9b95/vMTCr/z8ec/dzOrAgAAAAAAADdY5eHcG7V57vxlrbVXPKpza+2l2qw11yX9zOV2DwAAAAAAALi9Kg/n/vr5z98u6e+31j4+dWytfaik/13S+503fefldg8AAAAAAAC4vQavOdd7/8HW2ndL+hxJHyLpn7TW/pmk/1PSm7VZX+59JX2UpN/z0K/+pKRvnWqHAQAAAAAAgNuiEgghSV9w/jt/6Pz/f+z5f8/1YGXEn5L0ab333ayYBwAAAAAAANxgpYdzvfcTSZ/dWnu1pK+W9KrQ9R2SvknSX+y93x+3i3iu+f2V5quHolYKiXo2gCj8vk0hC5tK6WYumbUth2+vuwQy+VRWySfChGCjUqTt2kT3tHDcwqGw/4g8BQ2517GjILzzDbq2FIMT3hOThjYLCWnz+fB0M5fMehxSzFLamEsycylmm75n2+OGWKKjEHl0NNse42Dmc3FckllKN3Pt8xCDNIuxjLuxDsliK9O+DhOJS2tNiaiu70m42O/M/Rj3VouttpRk545nSk1LiXMVLhUuJk+biSQmToZr0qZZpmvdvX9pPi1MIzlJ1LQVkl0rqaybdtO2fUnHMVKKWdqeTWuNfU3b0h8MmwoXk1bDGxWSZ/0Yg7uWlNLnQ3pqt2mtw5PjJambdGeX4CrVkuNTPbN271XpYPh9c/NvKiNjGnzh+nXnRTrd8uszSbmxRjFtKa01tM/M/drVLZI0N31dKqvkk1nvLPzkkuoZlxKf0lrdGHdMjSP5pFXJ1zmxRjGTX0xrNe/pTaxnXC2S2lPf5Xp73JP1dn2yaffH06Wq3l/543ZqYp9d+upF7RXdXNeubpGktbkm87U+vJ7ZJX+v9X1dfZCTT01boW6RfO1Sqn1S+nyqUQoJ9q4eiePa+6HtGrfnbjyl0/s5vz8/2813z6rfnJMk9d6/S9J3tdZeKenjJL1c0rE2D+X+haR/3ntPKckAAAAAAAAAdMmHcw/03t+szXpzAAAAAAAAAIoqaa0AAAAAAAAAJjTqm3Nts+jWi4aM03v/1THbAgAAAAAAAG6b8sO51tr7SPpKSb9P0vsrRgS8i36ZbQEAAAAAAAC3WemBWWvtcyS9TtLhg6apdwiPNjtdadYfkRCSElhtilXo61JVY7pZSk81CT0pTc2kIMVEr4qU2mMiWuLrMP8APARQldLNYgTr2L5TSIctpZuZ45nSWl3feei7MElohyHayCWhbdpduplPLHNJZimVNaWeuXSy4xD3ODatNaWYzUOc5dgkUZciepGzvn2bWVXSWk2qWOqbjltKcXVpca6tKh2j5dzMcWGMlZs7w7juOotJ0CmVsfC2ur7lYsANkpK3CklopYTSmKpq2lJimUteS0loZyndzCW+Du8bX59JcU19c+qZGWNEutkjjUyf7+YaS8PGtFaTyipJfb19LHpIWnXTr5kKH8FdaJW05dDVzBeuVtu0hzHc6yslsNbOi9Jtxx224nzo5tRUzxyYuiOlz7t2l3Aq+VRWySezplrE1zMhHTamz1fSWof3nZtzYJGitYNd1TOrsOqTS3GdIq31fjPJ8ab2laRZSGCdhVppKJfiLPmkVUlazbaPhWuTxn8+SCnzxXJ0uOrnrbGf2Srp84W+kq8lYqJ8ITk+1jOFRHm3vZRGa2uXavq86R+nkAG1y2zfaa3n35j7G9o8mHtwObxV0tskPT39rgEAAAAAAAC3W+Ux+38j6VibZ77/i6T/tvf+czvZKwAAAAAAAOAxUHk497u0eTD3fb33z9vR/gAAAAAAAACPjbQ6lvPu5z+/cxc7AgAAAAAAADxuKt+c+wVJHyzpLTvaFwzUlmu1tCLkuR5WuXXBBmnl8OYWGA5hDvE5r1tYO4ywNitidrvDeeFLG9xQWji80DcsFtnSCqVj15mfJBxjR30lu6pmCyttzs1Ct64ttadxD8KJcWDewEV4s117Wig5LXTswh/SIsxjAyFS8MM8vL5ZWnl2oHW4JpOFCa9JiyK7BaJTUMi8j3sdkl/0+SAkurj23Nfvm1vg+Wzlj4Vbqz6d93a6jwvH7ypBJiiFB/hmt8uVl1eZ19N+5DEKAQ2F9pYCIVzfsNiyG8MFPEgXnBYmBCEuXj1JgtO2Sj2TXoe9LYfwiPgCXVBEqkXc9ornoa9RhtcXtXom9K20T3Ct70zxBbr3L82/dq4Ou+HuAWnh/9y+fdK5ttReCX6QfD1TqVFSyIOtZx7x2ea5dlXPrMK9fWX6z8L93rWnvn7nQnMIsnFhFctwPF17qp9PW6jXRtbxy1RLVmqUXX6mcaaYt0r1jKsvQt9dfdZNmQdpDFOPTBJw5cKpwrgxzMG1j6hnUu01VuVT1hvOf37MLnYEAAAAAAAAeNxUHs59o6T7kr6mtfbyHe0PAAAAAAAA8NgY/HCu9/7/SfpCSS+T9GOttT/QWvHfNwEAAAAAAAB4xuA151prXybp3SS9UdInSvpeSaettV+UXVnsXfTe+4deei8BAAAAAACAW6gSCPFaPbuM4YOfR5Le/xG/13T1S78CAAAAAAAA117l4Zzks06myD9BSc9JJOfSm9Jd0ldIQ7RxVSFFJY5hE72GnzIxeS29fLu91PlqT92rDkncmcJhCyF7o81iatrw9tzXJA2l1NmYiGpSR8PfUdiE0pTAatqrqazzsSdiMSXVJaLGY2/2eZZOuJD4WuGPZ3qvh59DiTsv0hgpGfC22NnLm2DcUupoJbVyguRLm95WCbeL+xb+oJBU/qi65LJK9Uy6nswosTIoHKMe6hk3Rqluuajdbs91LtyAdzjd3PKpzM7V+R43/GDkusPUDKnvyDTTxNVJqb1yn6ymr+6snomTznbTOnRem9pnHd5/WzPGYzy+dq3VxDua1686afWauBbzYbyHD2xT8XVU6o5Yi1TqsquqZ3bzZlYezr1c0qmkM0mnvfdH/VNWAAAAAAAAABcY/HCu9/5ru9wRAAAAAAAA4HFD2ioAAAAAAACwJzycAwAAAAAAAPZk65+1ttZeJek1kn5d0jf33n/mvP2PjdlQ7/1vjvl9PKy9y0qavbKqpnscG37frXPc5+F5btqHuVuw2XfVzCzYnPpW2guvb5euenvWFGtXTrAmZ1pEe6h1+P1Ke+5rFvENfVemryStzTl3FgIM7CK+6e9NKu9fOMRuTdXKIsxukWMpHwsXCFE79sO3lxZmdvuQ2lfxvR5+DiVu/9IYY6+R6869vLgG99iFo4v3EX9PTasiT7E907WybzFUxPy+uydLaiu/PZlpq6e+ZoOlcA0V6xlbX4TrqVCLpPPNH/vCGFPUM0Fl36aIdqtck7d8KrNzdbzHFQ50rju22+N9y93jwv0w3j8L92UXeLBO9UUzY3S/2lIKR6gEVKXaxUn1jKvjcn0xvA50xz7WlxPUrrWaePwF7K6R0q3hij+77NK1+OyZ7rOulijeJ5v5LB/TI9x9OX1udM8fVv76b+n1mfBKGyy1GWS7aeuk3c2b6WbB10t6n/P//fskfcD5/36dLn9qd0k8nAMAAAAAAAAe4h7OdT37KPC5f37L//4LAAAAAAAAuDru4dwflPRntPlnrd/yUPufuZI9AgAAAAAAAB4TWw/neu8/Kek/Ne3//ZXsEQAAAAAAAPCYIK0VAAAAAAAA2BMfi4NrrS/m6guf+PNMn0qSVkxgNX1dCovkE1ckrc3YIYBIazNGSpYLgUd27LQ991py32G/f9EYpRUbJ0hTG60a/2LTmIanP63WIfHKtK9mvu8ynBhLk+h1Ft6o2Xp7WnSJqlUrl0wmn8i1mC1t34WJVDwL6WZpn+eFZNaKmABnXl9KLHPtlb731wvb98S8p5J0atpdmyQt19vbW4Zz9nQVzkN3LqeUPZMgtQ7bW7u0qRw5Gdp982hj573QXnl5sW9hvjdv/2bXzLGP97hQcZkAsXiVhlnLd1655DX/RveU+Ooi9dwOh6EnCdkrJbynlDb3PhWSXUP7+mCCvrGWMGNMUftMMF2MrlF2Wc+4k674At2cmubf1Xr7ak3zursHpPvFQYhQXrbt/qfh+p2nlMQKNx+GhENXoyzM/qa+sW4ppLJWVOqWTf/hCaxnZsKv1DOpbqm0pzp3adpT/ZxSY21tnmqUwueD2gfY4V0nMUGytv9s7bv6z6nhHr6jMdLn7Th32saQEm/u17NU/bjpMHyuSknAzaW7ptpuQD3T0y+PxDfnAAAAAAAAgD3h4RwAAAAAAACwJzycAwAAAAAAAPaEh3MAAAAAAADAnvBwDgAAAAAAANgT0lpvoNXhXKvDSyaEuCSzSsJLSEKLKTEHw9NaK+lmKfXOj5H6DmuLfQtJf6m9lOw6RUpQRUi76SF1x6UxpUQv13e58gfjzCSzzkK62WyC6CaX0uX2V5KW8+HJWy6ZLLUvQgKrTTcLr3k+QcJsRUr0WpsTNPWtpLW61LuTkNaaEstcMuu9lR/DJerdD33PQqTimRljGc5ll3qWrieXTJXSqlJ7xRQBgHaQdH+ppJu5uTodt3VILHNjxORpN0ZKJrPNtntbpnvt9vZmLpVVUjPtLbzmNHXatNY0tUxxctmdSAmsZhdS30L6fEprLSXKmync/f6F2zNjVBLsK7VWrH0Kia+5Rhl/XoxNAq7Oh/a0D/OIm6vTvH5mXsgsHJ90f3HcfTbJSbLhXmtqsIOZnwR8WquvfVyNchPrmZzWOi59PiXHn6yGp7WehvfUnVsnSz+uq1sk6cye98PTWiufD9I9/MpT5iuJ8hN8Fqykp6ZAYzuXxePmPvf7nuH2advXIVXVXdZ9GfbMvG5bnyjXOd2lCY+oW1YhiXqsSz+ca609IemDJL1Sm3fzh3rv96baMQAAAAAAAOC2Kz+ca629r6SvkPRFkp586I8+UNK/eU7f95T0H0n6gd77U5fdSQAAAAAAAOA2Kj2ca639XknfI+kJvet3H9N3Av++pN8h6fMkfddldhAAAAAAAAC4rQYHQrTW3kvSd+vZb8t9n6Q//Yhf+25tHuJ9+qX2DgAAAAAAALjFKt+c+zJJL5B0IulTe+//RJJaa19/we/8i/Of73e53YOzOp49OhAiLXRsF/Gt9PWbywsru7bhoRJ5AWW/H65/DISwCzaHvm4B5ULQhBQW/CwstpwW5ZyE++5r2GBeXN0c+7A4rFvrvIVFkdty+4DmxdlDuzlxXaCAJB3OthcvPg3BD24BXklamDEOwoLGvq9fQNktHD0Pq7OnRaZT+1DpGFf6p0WRXRhHCnNwi1dXgh82Y2z3j4swm8WS02LLLjxCSoEQfp/dIuOrtNiy6RsXUC60t8p7nc6rCUIe3BBxnnULx4eFkvO5vD1I6uoCCFxogyS1VTgWpr3N00LHZnvLwqLI6W0qhDyUppC88rSXVpk2SvWMO7eK9YwLqCoFXBUCpzb9xwVcpUCtsWFYkmoLoO8oPCLNT/Z6L9Qtkp9TV/62rJm53pcppMUsIJ5ecZqfbABFWEXd3eMOZj5o4nDmx1i07f4Hpm7Z9DUhD6n2cQFXxfpkV/VMPPamaI+BYS7gKoVuuNon1Rex7hheo7ixUyBErmfMeRjq6qWp49PnA3dN5kCIqw2KyHXAdlvIQKiFPNiwr0o4lfxcXajLTB7Mpm8hlMvVLVIIrTrwB6NUz6T3v1LPDKhdVungjFQZ9dO1OQxf/+DB3ABvO//5HqW9AgAAAAAAAB4DlYdz73X+c+iDOUk6O//58sLvAAAAAAAAAI+FysO5B993/o3C73zw+U+SWgEAAAAAAIDnqDyce/P5z48o/M5nnv/8N4XfAQAAAAAAAB4LlYdzP6DNsoJ/qrX2vEd1bq19qqRXa7NO3Q9dbvcAAAAAAACA26uS1vpXJH2ppA+Q9Hdba1/Ue3/zczu11o4k/VlJX33edE/St4zdUTxr9eRcy6NHpLUGNqGlEoJTTWt1qTSVBJuUmlZIFqskoVUSWCupaVL1WJiUmJhuFtrtwIX2lN6XUs9M4uA6xRXZJDN/QF0y0Xrtd86lmEnSmWlfzEKC2Nykpy59upnrK/lk1pRu5tLGUrKr65uS0OYjU8yqVoVkT5e0KvlU3ZiEZhPrQhJaPC9Memrq65JWY9/hY6QE1lK6mbn2XNtmkMKEMcEplIZwU0Ml/SsNbOfqtA/pUJj7jksV2/Q12wvHOEwBWpv2lJ7q2mMSmkkbi+MW7g0tpcJdsbH1TL7/TlHPmGsy1iK+3Z5bse4opMOaMUop86F/6utexyRnUCk1eHjdIvlzax2+z3Bm2tIl4u6Tq7BvsZ6ZbR/8lMrpklldQvymr58cXHp86mtrlPBGpZrIuS71jEuUz2mthb6uFgkXlKsjUv/U19XEyzRuqlHcGKZukaSV+3wQzm97TcaUed9sT7nKKVT8vOU+s6X7iJ0P070h7IbT0n3LhDC7hOm0H+swR8Z6xtYovq+rXXLfKdJaXd/Lzy2rk8s9i3mUwQ/neu9vb639Z5K+T9LvlvRzrbU3PtTlr7XW3kObh3cH2pzCXdKX9t5/ZcJ9BgAAAAAAAG6Fyj9rVe/970n6TyS9VZsHcJ+oZ59DfrKkD5G00ObB3Imk/7z3/rqpdhYAAAAAAAC4TUoP5ySp9/4GbR7C/UlJ/0jSfW0exj347+clfYOkD+q9//XJ9hQAAAAAAAC4ZSprzj2j9/60pG+V9K2ttSbphZLuSHpH7/1kwv0DAAAAAAAAbq1LPZx7WO+9S/qN8/9wBU6fmKkfl7/0uFFZC7yy2HJaNLiwCLNdJLMQHiH5RZGnWNDYLpJZXUDZLsI8fJHxUnBHEjq7xTPjGplhEdBuVnhfh+VMuzmg3S32Kb+4+mruxz0Li7bPzQLK87T4sRk79U2LFLtFkdMYPuQhjGtWM22FfdgHtwCyC/mQaoEQbrHsFNAQT2XTP4U5uMWL06LRKeTBtaeQB7e9tDCvXUA5BkL4Zts+xSlUuWfERYqHT3Ju3eiUS5PauwtoCMfTLdBv1k/f9A1zXC3kwTSmRZHdYsuVUCCFxZJLN50dcvVFSvmwfX3XWKOUQiUqfcN1PbJGiSEPI8fNY4Trt1DPlKRz2V0jaXvpujZFWA9Fkbuf9RTyYOqL1Tws/B/qHFd3zGLA1fB6Jh0iF/5QqmfCG+Vql+tct0ghoCwcOTdGpZ6JfSs1UapFXD0T6qcUcFSrZ4aHVnU3Rqpn4meawoQ/fNgLQqTMGOnznatzU89CMGOsc9xcHY6nrUXC56qQy+drl0rAVaGeKdUtyYh65izM3WPtZlQAAAAAAAAAj/TM3/m21l4l6cMnGredj30k6VjSce/9z080NgAAAAAAAHArPPwPMj5L0n+9w23xcA4AAAAAAAB4yHNXS9nVQiLXYxEBAAAAAAAA4Bp5+OHcP5B0/xH9n5T0lZKekPTrkv6OpLf8/+zde7Bk613e9+ft3r1n5lwEAoEAgx0cLGIcbDBB4la2hUSsQC44AoIxYCCkILYJkGAHqsLlH9kmLgTGgIMJdwrHXIwEMQLZQQiXKSqIEDAlI8UWFiAsI3FA55yZs/fuy5s/ukcaTT/PnvXbvXr3zJ7vp+rUnPPOu9+1evVa7/r1Or3fR9J88/fXJT1P62/hXZP0A5K+TtI7RttjAAAAAAAA4Ip458O53vu/kPQvUsfW2pGk10q6ofVDty/qvT8T+r73ps9nS3q69/5Xx9zph93i0SbduNiXHEuhJDsmoeUxfNdKolcpxbWSplbYXkoxy9vb/gJpTEIr7FuO3yuopJulv1hsN/WUnmlSC1PSm0uoTOlmKemrmbFjyqnp69rOHaOQQubGyOPa5tK+7UtKYHVCaKVPvUvJrq5v2F4cw6Q/uZTU9RjDfl7K571NLAtjuNjR2Nclb8VU1kq6WRjDtRfuAVKYz8J8uLLXSHgd5tCn9NSUbuZy7G1KamjPfdPBGD7Gzn0rKWaSfaPuk0BFfw7tmkh/zhg71ygj1DOl2meEeianuLp7eOhbeZ9SuzvnYlKfmzvDuIm9F6XUwu32VbhZN5NEuUz1RUpEtfWM7WpTXKv1jBuaeuZd9lXPVBJjpVTPDB8j1RdxP0zdkcaw7TFR3tUzI6TPV4TzLd4zSvGa7sacjptpK6SySv6lxJph1+R4qVSj2MNcqDtjPVO4fneZWhbT/fzCaeV0+iJJHyvpF3vvn5sezElS7/33JP0Xkl4v6Ytbay/ebTcBAAAAAACAq6fycO6ztH5u+Y1DOvfeF5K+Xuv/6fJX6rsGAAAAAAAAXG2Vh3PP2/z5hsLP/Nrmz48q/AwAAAAAAADwUKg8nLux+fPDCj/zxzZ/vm/hZwAAAAAAAICHQuXh3Js2f/6XhZ95yebPf1/4GQAAAAAAAOChYDLJoh+X9B9L+vTW2i/23r/hvM6ttb8s6fO1Xqfun118F3G3s8eaVo9cMCFkx7TWJCeJFvoW9qGyvXHS1EwyWdqHMdLU3HFLj9IriVcxVcr0LSQLpjFiwtKykvRmBg59l4W01nTcbCpY7Ovb3YGOXStRQSMklsV93lE1BNIPMjzdrLS5mHpm2irnfWFcSf56qPRN+2Z/3jeXruvSZB2GTcmA5oqI76mdD8O47nWHSiemkNm5LIyxa1+pdDKX0s3sz9fm9ZJ9BSqOMWftWouMMUa1nqmktpfS54ff4yp1R2nfwvVbkjZo4jPTeR/natce0iXtawmvz+3GKqanVuIXQ9dCsmulvijVMw9h3ZK3V6hnRjhnY+psqfYpJKWO0XeEzyM+2TP0darz4Y5jtBT569JTi/XMrgnvORW7MEZQqWf89kaoZx7wtNaXS3pi8+//a2vt1a21T2qtPXa7Q2vtuLX2Ma21H5L0nZvm5eZnAQAAAAAAANxh8Dfneu9PtNY+XdKrJM0kvWjzj1prvyfpTNL73/Ejtx8nfnnv/fXj7C4AAAAAAABwdVS+Oafe+2skfYKkN2r98O32P8/R+sHcnW3vkPTZvfdvHXOHAQAAAAAAgKuisuacJKn3/rrW2p/QOhjiL0j6WEnPlXRd0u9J+lWtv1333b33J0fcVwAAAAAAAOBKKT+ck6Te+0rSKzb/4JLNH5P6oyMOOMJirSOsGz5KWMXOizOPsVb2GOEYhQV/S9JCm2Yh0bjobApdcIuRxn0e3tcGQhSPReW88GuOFq+RPQWv1PZhX6uzF40yObhx99V3hP0dZaHcQnhEJSRxZNm+AAAgAElEQVSg8Poqp1AaNm3NL0q/+wRs560078V0m+Hb23nh6RG2l+xtCtjXNV21pxc4RsjDzn3TfoxRJxX6jhMMVgigGCPgqhCcExdGrxx7M2/ZuiXuhG8e5VY0wn7s3LeCuuVifaXd93mE+9YogUMj1DO7BidFhaCXFMbh9i0GTZQu30I9M0YdcUmhC/fcj4qRr+tF6fdPh9vTsAAAAAAAAADuhYdzAAAAAAAAwIEM/rXW1tq3FcZtko61XofuuqTrvfdPKe4bAAAAAAAAcKVV1pz7Yl3st4bbBX8OAAAAAAAAuNKqgRCVlfS6pGckPSHpHcXtAAAAAAAAAFde5eHcCwf0mUj6YEkvkfRpkl4v6VN7779zgX1DsHh8pf5YipgCLqCUOHi5SX3tkrcH4AJGSRwMSWjlncF+8Y7gPpBSmEun537qC+oW4AHhklZHuXypZ+4/4x795WQ/z2IGP5zrvb92YNfXSPqu1tpnSPpBSa9srX187/3sIjsIAAAAAAAAXFV7S2vtvf+QpG+X9Kcl/eV9bQcAAAAAAAB4UO3t4dzGD2v9hdHP3fN2AAAAAAAAgAfOvh/O3f5V1ufteTsAAAAAAADAA6ea1lr1os2fsz1v56Gyemyh9vhivAHHWPlyjDUWK2Psa5/jAsOF7aXXUdiee321RY5r27NjF/atPIYd1zfbrY1wLErHc4xjH5Tf18H7cJ8sSL2nF1h6eSP0HWV7u04jhXHj/lbej8qxKL5mP0ZlEihsL+7DnrZXNcJ7snPfir1NWkX7muP2dW+47FqrUs9U6pbi9kr1TKG9PIbtO0I94/atsAtjHPtxjsUI29sVdcsd+1AcvFJfVO7hhe2NkiM3Rj1TGmNg2zl2rmfGqBnH2F7FvuqWvY4x7nW90ojPYu6wt2/OtdY+S9JXa32tvn5f2wEAAAAAAAAeVIO/Odda+xv36iLpPSS9n6RPlPRBm7Yu6ZsvuoMAAAAAAADAVVX5tda/rdoXpG9/AfFHNsmtAAAAAAAAAO5Q/bXWVvjnzZK+TNJfHGtnAQAAAAAAgKuk8s25zx/Q50zSH0h6Q+/9TRfbJQAAAAAAAODhMPjhXO/9e/e5Ixju+LG5po+fXehn+56SxXroazeX+hbSsdL2fEJp2N5qeNRQpa9c39Q/heus3D6kvsPfpxhUY8Zw+7DekeFjV8bIfUfYN9d/jOTaoDLGXlNj9zWGE0/DQsTSrglbQSlsOX2fvJQCGsYwY8epbGLeqLBvbowW+4YNumPv9iH0rR+LwutzfSvpbZV9SP1DX3s4Y9/h20tvk3t9pSC8uA+7Twxxn3dUqlviINs7F+uINMTAcddjDx3gnP0w9+V4LPZVzxRqrViLrLY7233QOedyqX4a/vpyTbRj3/T63BiptqsksFZqn2CMWmSM+qmyvZ3tqW5JKvVMObh2x/qpVPuk/uEk8n0r4/q+6abjxmjhvlyrD8OJ6Gq7EeonW0uMUM+kY1Gpny67nonbGzhu7LpD3bLs84v/8Dn2ltYKAAAAAAAA4HyVtNav2fzrN/fe/2DgzzxX0p+S9Nbe+69eYP8AAAAAAACAK6uy5tzXaf3F4h/Qel25IT5U0qsk/ZykF5b2DAAAAAAAALji9v1rrU9q/evDz9vzdgAAAAAAAIAHTuWbcxfx0s2fN/a8nYfK44+eaPbYxZ6rVhZANmv45uCH0L5yizCn7a22X1Na/Nj1Xbeb7YWFee2+hb6uvS+HL+Ic28MY3S6SWVh4WruHElQWP5aktjRtaVFk2zeMaxdbHt5XUmlxZztGccHmSgCFXVg7van7WQ91FKXFi2N4gFtJdvgY1fAIt2hwZVHk1LfSHrc3MXNOGndqFlxP407DGHaR4nRACydXZQFl8zrW+2HaU1+3EPTUX8BpUWQXppEWI26T7bEnYdzU7sau9J2mfXP3keLEYE7D8hj74uoOV7dUfv689mUhVMLWIoW+cYxQ+7ixU4iUr2fChJEOqKld0vzra5fh40YxSMPsQ5paTC0ihboj7FupZhihnikFXBW2N0rIA/XMHWPsWM8UAqAkhVCCwhij1DPDAxp6qEVszRD7hnb7+TX0dT9eOMZSqF3SGK5vrJPMPTzUPrmeMfVFrGeG1yITU/tIvmyMfd24I9Qz4TQcPO5Q8+XJhX/2PPHhXGvt0yR9mPmr/6G19sQ5YzZJz5L0fEkfp/UlwnpzAAAAAAAAwF3O++bc+0v62rvamqQvKYzftH4497LifgEAAAAAAABX3nm/G/kKrR+u3f7ntlb455clfWrv/Z+OvucAAAAAAADAAy5+c673/luttU+QNNP6QdvPaP0tuM+R9JZzxuyS3i7pLb33d4y4rwAAAAAAAMCVcm4gRO/952//e3vXCn+/0Ht/0z53CgAAAAAAAHgYVNJav0/rb8U9tad9wUDvfeOmjh+5WLqISyhNbBJaiJ9ZFBLEUl+XhLZMfUO62WK5HemTXvPSpJPF1DTXdxoSY1MC66KQsOuS0FIST0reilFBRiVtLKSbTRbb24tJaJW0VtN3EsZN+2z7x74mEaiQAiuF9NtCQloMD9pTEtpe2VTVFOtnUkf3lLQax0hppqX01OHtqe/KtcdkMvMCU7pZOjFcNVA4D1NIdT72hcQyk04WE8tMMmvqmxLEpqb/NCS+ujGO0rjh9U1Nf9cm+WTWlG7mxkh9J+HNtkmy90laq01cT4mo5kSMtUFKg3c1SmGMcj1j+qekelfPLEMt4sZI9UVKcbUvO9yX3S2qpaTcMF+43YsJrK4t1iKFRPlCsmuqUUrjVpJkQ6pupdYaI/G1kuxqa5f7Y2oZJYG120Iw9HW1SDUx1LSvpumCcn1r23P9c/1k2lKatNnnVTgxSqdL4SNYmlzSPGnf15Qo7+qZo1BfmL6pbpnEGqVQz5jXfTT1E1GpngnH09VKKT3V9a3ULcku9czZ2a0L/+x5Bj+c671/3l72AAAAAAAAAHhIVZ4jAwAAAAAAABjR1jfnWmt/UtJXSnpC0rf13l+/af/cXTbUe/++XX4eAAAAAAAAuGrcr7W+QtIf2fz7SyR9yObfv0cXXw2ga71mHQAAAAAAAICNtAR0C39fWGEeAAAAAAAAwHncw7n/WtJXaf1rrX//jvavupQ9wj29z/Wbun4jxVWezyWLrULkoEsmq6a1LkyczyJsz42xDH3nJpVVkuZmjNR3YVLIUt+lSaVZLFJqTziepq3H2CXXlhIufbMfNySTubaYZrp7utlksZ++zSStSj45bYwkWZfsmvrnhLTh6bA7J6HtUUxgdX0LiVc5rbWQnhpTVYePsTKpWZVUVklaHZntpYQ0c3eO23PpqWHYmBbnfmCM/x1XSXFNidSFdDOXQpZSzI7SGC6BNSSWHR9tt6dkslkaw0xQKa3VJZYdhcnM9x2e0rZuH568dtl2rWcqtci6v6lnYgLrdvtZiEOMdYfZ57NFqGcm2+3LkCjvapdVqFvC7V5uyep4y3FptGPMLXF7ZnPl9PnhfW2ivPn53Hd40moco5LsWqhb4hgpHbZSo8SblBv3PqlndkyUVxi3kj6fElhdf5f2KflaosV6JtT8pv8qxEy67RXefk3CcVule5E7njH+dniqbmx3tUtKlDd1xySltbp7eKH2kXzdcVToOwu1SKpnpubCTgn2rh45CpNZJa31suqZk5OTC//sebYuo977r0j6TNP+9XvZAwAAAAAAAOAhRVorAAAAAAAAcCDhC6jb7khr/ce996cH/sxjkj5A0lO99393gf0DAAAAAAAArqzKN+e+R9J3S3rfws88X9K/kvS9hZ8BAAAAAAAAHgqDvzl3QW/XevnED93zdh4qzzl+Wo9cm5/bxy1+nNqX4RmtDYRIAQ1xYeXtVUBPw4qhbnHms7BQ8tnUj+EWVj6d+L6nZgHlaViQ/CytlGoswkLAMouqrsKiuvZwukWVpdLCynHdS7cealqttbCwclyk2PUNp7RbLLmyiHNqr4VHDO8r+cWS4+LO7hwoBUIUFmaWpHDODWaCGM7jpiIX5iDJB0KEvt0FNMS+Yd/MdRbWbFcrhDmkRZHde+1CItaDm3F9T3vc0vmW3j+34PYYa3D3uGCzaQ//u7C5kIcwV7vwh9nMH4wY0GBCHtKCxteOtieXWZgYjmMgxPYYLiRivR/b7ddCX7cfaQHltAjz1Jx1abFlp9JXyrXL0L6pnnG1SAq4modJ4NQFXIW+LvzhLEwMqc5x+5EWr16Y8z6FRzRTz6SKMi2ivjLb62nyNNeqDaCR1CoFTQxDMm2hfsr35eF9Xe2Sa5Hd65lKqIStZ1KYw77qmRgIYcat3nP2VM/EacgEE8S+NuRheI2S+rpwqvUY5jNGIeShhbpldRS2Z1Puwj6btlQG2mFTbRDmEVe7pHPLNo8QCNFijbLd7oIfJB/+MDP1ieTrFsmHWV0Ltci1qalnCnWL5OucFER1zYyRapGZGcMFPKzb/bGv1DNDapdbx6f37HMR+15z7oWbP5+95+0AAAAAAAAAD5z4zbnW2gslfZD5q5e21v79OWM2Sc/S+lda/xutH0r/2i47CQAAAAAAAFxF5/1a6wskveyutibpbxfGb1o/nPuW4n4BAAAAAAAAV955v9b6Y1o/XLv9z22t8M9bJX1p7/0HR99zAAAAAAAA4AEXvznXe39Da+0vSTreNH231t+C+ypJ5/1aa9c6COItkn69976f1fIAAAAAAACAB9y5aa299394+99ba9+9+dcf6b2/aa97hXO959FNPWYS4oZYmsibmFhm4oPcz6e+kk9mvRGT0Lb7ukRVSTpZ+oSWMxM31ELiims/C30LAVQx3ay7QUKaj0uE2iubAhrSzQpJopUUslLfQirrut2kI1XSzVJaaxjDvZackGZSnlLfQrLrKFGbFemcdYlloa9LJ0spXS5tzLWlcSWpmxQyl8oq+fM+pbLGxDJ3jcSkKDNGOBYuiDC9HSkl0Z1HKaTap5uluSxsz+1zmA+bed0ulVWSpqY9Ja2mdDOXNuZSWSWfbnZ96rMvUwKrGyMlobl0M9cmSTMzEaUk2ZTiOnUJaYW+Y0h1h6tdKn1TKmulnknp8y6Z9XTpj/1ZrHNmW20pQe50cW45/25cjbJKc2eqZ8wxStevne/jvJDuWyZdMvS098RUt4yQPl+pGVztUqlbUv9JoUZJ9VOpRkmvz91H0riuRnkQ65lUX7h6plCjuER6KdcobtpK93A7baXY2ZRy6j4XpcNm2kOAsr3f51TWtM+Xfb6YtlTPmDluGt7rI5ccH1Pmh6enuppDkq67WiT29XWOS49PNYpPnx9ez6Saw/WVfO2yS91yfLSf758Nv5tLP6f1JfrMXvYEAAAAAAAAeMgMfjjXe/9ze9wPAAAAAAAA4KFzXiAEAAAAAAAAgD2q/FrrO7XWHpX0ZyR9oKTnyj/kO5L0mKTHJT3We/+LF91JAAAAAAAA4CoqPZxrrT0m6WslfaGkZxW3xcO5kTxr+oweDwtC3rYsfCkyLX68Mgsrp76pfda2FzSeh8WPJybkYWICHiRpGhdRH74IqFtINC0uujSLjK9W/hgvw+LjzayI6tqksKhqXHB99/AIG+aQNldYWLnSNy1oXFqYOS6g7AIhho/RYnhEWvS5EPJQCIRwiQJ2UeXQd93fd99VWpfXrQRcWUBZaQHlZWFR8yN/rZbWDPYpCL5v5VoNU/XStcfFskO724Ux3v8xxihNOub6DYst20CIFB4R5uqZub+mIAW3KLJbVFnKCx3fMGNcm6TFlncNhBjeV5Km5s2ehIl9Gk/Q3aR6xtUoy0LA1TzUF7meMcczHLfTQpBGpW5ZpbAZM/ctQ193PaS+sc7Z8fo9J87B23HOKc97rn8hTKcSHpHDHMIYrp5JY8xdmIMfty38C3Rjxxpl54Cr+7duWfcv1PG7BkKE4Id0HvZKyEOpaAhzgytnwn155T7/jPFZojDllE6hamBNKeDK3BtCLeI+685CPTNLwVdm0kk1yvWj7brD1SfrcYfXOWkMd/+8Hmqfyv031SipvzOknplMzwaPVzH44Vxr7UjSP5X0fNXuqmeS3lLcLwAAAAAAAODKq3xz7oskvWDz72+Q9F2SnpT097V+IP3nJc03Y36UpL+++e8X9d7/1Vg7DAAAAAAAAFwVlYdzn7n581clfULv/enW2jWtH85J0i/13n9/8+//V2vtlZJ+XtJPttY+ovf+jnF2GQAAAAAAALgaKmmtH6r1N+T+Xu/9aUnqvZ9Kuv1A7jl3du69v0HSyyX9EUlfsvuuAgAAAAAAAFdL5eHcY5s/f/mu9rdu/vxo8zOv2fz56ZWdAgAAAAAAAB4GlV9r/Q1J/5GkR+5qf+Om/ZMk/eBdf3e6+fOPXmjvYN2YnOmRO1JgVoVnrEuTbpbSS1xi2aSHNKfQ7kxD35RC5vv613xkElqOQmzW1KTGpn1zqX4pXWcy8fvmAo9i0KpNQts9lbVkjHSzNIZ7eaXkppT+NcIYJm2sksq67m82GNNoXZpaSPy1CXLhRReCE2Piq+HSyiSppWnIRHq1MIZLMusuEkySptsbjLsQD8b2T6TLzF6/8RwqHKN4zrrUuzCuO7/DW5re6p3TzZLKtJWOvWsLL7DSN6V+u/uISzyTpCNzH0jJrtdiiqtLax2ehDZGutms+e25JNHj0Nf+fDHBtVLPnJl4wlRHuNTZaUr8LdQzFatwgqf2halz0nm4NKmq6fx210MsRQrX2RjX+hjsLleCZFVNn3epo4VxU3pqoc5Jl6Qbu5LKuu5fqFHcfSuMa29GMX3eN+9au8TTMExDzdT38aOLqVFS52ZqnzQXFkqtWJe5WqJ0fhfH2PXzQbSnFN+o8JmtUkumlFv3OTOlzFfa0+di+xk6vKn7qlFSUn2tnknp8+Z4xpN22931zCLs664q35z7N5s//8O72v+J1qfr57TWPuGuv/uMzZ9vv8C+AQAAAAAAAFda5eHcz2r9EO6/uqv9/5D0xObvXtFa+9rW2n/bWvseSV+h9XPt/3f3XQUAAAAAAACulsrDuf9d0k1JL26tPft24yYc4n/U+uHcsyV9jaR/IOlzNm0rSS8ba4cBAAAAAACAq2Lww7ne+5OSPlHS+/fef/+uv/s+rRNZz7R+IHf7nyclfUHv/XWj7TEAAAAAAABwRVQCIdR7/8Vz/u5bW2v/UNKnSPoArQMkXtN7f9tuuwgAAAAAAABcTaWHc/fSe39C0vePOSa2TVt/txSuqbZTSZYpPsg1x3AkF/0T+ob0VJeY4hJjJZ8K59rW7SFtyrTnMS435ielnl2qMXZhjCSlXfejmrxWSIe1Y6e0qkKyWEy8MmPEU8Uls1ZTzArpZnbc9Bchcc6EPmfutaTIK5v0FpLQ4qEYfuztOVRImF5vcNi457VfGYVDV5k7a/eR4WmWcXvmTXVt57YX9tkli1Xuh5VUVskns6YE1pQOWlGpZ9y+nYWydqbtvstwEk5Cuz/2hfe0el4UbpTunC3tW/G9K9Uzlxw0f9n8vSF0rhy2eG9w963h9/vyvtl7bejr6pzKvqVk16RQz9hjlJLjU0q8mftcgut6kOEpvjZJNiR4xjBpM0/GBNY9nbNRJa31MvfhAC77s+Cun61zfbH7vc89I3DJ6rlvqEVijWJe3w71zBg1j1NZcw4AAAAAAADAiHg4BwAAAAAAABzIu33/v7V2a0/b6b33R/c0NgAAAAAAAPBAuntxjuta/1b22CtE3Ce/6Q0AAAAAAADcP+5+OPeb4kHafW/ZWw582FiF31h2YQzLwm83p+CHtLCy254Nmojb833Tfrj2NMZl6/fJfuys8jL29ZJTRkChPZxCfn3gSVo0OIzhFvINCwzbfY4L1ZvOcdHgEI5gW4dzCxeXpTHccd7nZWP2I70+116+pN15WDyXr4w9VRpjzPeVudrdz9I9Lrab7aXXYe+pcTIz4Qrpfh8mM9s/HZ64UvluKvVMPG5j1DmF7dn3tHpeFCa/XeuLvdZJV/xTha0v0uEsHOZ8bxh+L/K3sxHekPQ63D28EDRRrS/2Vs+k6aJw7O3Yhe3lfRse+JduDfs6Z0vuh88SB3A/fBasfLbOn7d3q5Mkf09NARRuP1zAg3RenWPadqhn7vUs5qLe7eFc7/0/2MtWAAAAAAAAAGwhEAIAAAAAAAA4EB7OAQAAAAAAAAfCwzkAAAAAAADgQHg4BwAAAAAAABzI3WmteACc9plO+vG5fSqJZZUE1nmf2r7zlT+VTvt2+8lq5vuaMRZhe4uwz659sQpjrExyrWlL7SlxJ7cPa4vGSDzbZzqSTR0NfV0yVQrXKaRY9ZRcareXxthuC6eQJiu/PZcuOFFI/tkOVJRaOBhue2Ef0hC9dNK5gWsnRp8OT0R1KWTxfXLjmjZJWh0Nf6/DlLPzORTbKwlphb6XnryWjDDHuTk1zrOmLc3r7h4gSTNzHsV7jhlj3vxJdBSSxY5M/0m4Tm17ISR1FSaGlJDmaoZpIZV1EpKnk0oCnNu3VKO4Oif1TTWK6x9rItc33EhSuzu30nno0vDS+W3rGduzWOdcdj1TUU3F3nH+jdsr1T7D65x43zKXarof5m9rmL9ZhjfQXO89pc+bg9/SuMG+6pl4XlRqFNeeapSpuSZTGRjGcNur1Cg52XX3MezJNUaNctmJr4XTLZ2a/rPg8Hm28jlVCveRUKMsJtsTRv687cc4DZ9JLDd0Ki/cNBTe1Gmqn0ydk+qZIbXL6Z7uY3xzDgAAAAAAADgQHs4BAAAAAAAAB8LDOQAAAAAAAOBALu3hXGvt/EXSAAAAAAAAgIfM4ECI1tqLJb1I0sdLen9Jz5H0pKQ3b/75/yR9f+/9N8IQL2utvUTSd/bev2mnvX7I3VxeV1u+61lnWhDRcYsiL8Nij25BY7fosJQXNHYhD2kBZdf3LARNPLP0CzafLbfHPisswjwvLLSZFlterQorD8dVmAsrvI+wIGVpc5XghjH6mrcvrUeeF2DdfjGreOC2+4Y1dePi6m5R4z4Jfc3raynkwS2WHF50cR32neWFtc1fFBbArgRCpIWS0/ZWs+FjuEW04yLccdFn13f4GKWgiSBe15e+sPLwScdev2GeXdmgHz9hpDn8dLF935mE+aIaeOC4EITFdHjQwLXJYnDfWXMJNNI0BEKk/k4KldhVCq1yUn3h6pxKmINUq1FOl9vt7ucl6cT0TWPPTY0j+XM2B0K466lWz/hAiDEml6AwRGXX0rC1GsXM1SksyjS3eIxT4FThYLjb79x3XYZ77cTUHSm4wdUuqZ6xddzRg1fPpPu9r2fC9goBV5XgqzDlhO2lvsPbL/3zwb6mnMpns9ReGCPNs24On4brKX1+nZrCM9Utk+Xu9/Cle33huNlArYk/4WYmYDLVZZU6p1K3TO/a3s3lfr7jds+Hc621z5L01ZKed/dfSXqWpA/U+oGdJH1Na+1Vkv5O7/3n7ur/sZI+TNJH7rTHAAAAAAAAwBURH/m11h5prf2IpO/X+sFcu+ufP5D0zF1tE0mfLOlnWmtfcsdYM0kftfnPV4z/MgAAAAAAAIAHj304t1kf7sck/QW968HbGyT9Fa0fsr1X7/29eu+PSXo/SR8n6Ysl/XO96yHdN7XW/tZmyI+QdE3SiaSf2turAQAAAAAAAB4g6ddav1rSJ23+/Y2S/idJP9nNAgi999+V9LuSfkHSP2itvVDSt0v6EEl/o7X2Ckkv2HR/de/9mRH3HwAAAAAAAHhgbT2ca619kKSv0Hopw9dJeknv/feHDth7f01r7WMlvVrrb8y9XNJvbsbjV1oBAAAAAACADffNuc/U+ldQn5L0aZUHc7f13n+vtfYZkn5F0sdo/c25laSf2GFfsfHE/FGdzK+d22eZlxPcklJUbFprSlxJqWcmJeY0pKcuTHtKWo0JaSaxLI7hkl3Nz6/bt/suQ0pLandpaD0lu+4n9C5zu5HisVz6pkJYUSGBNZ1blTjaXti3EJ7qk6JCWlVb+u3ZxLIQeujTzXxfl8waU8xS4uueUs/i2+cSy2IS8I7JruF9Sgms/r1OSWjDx40JaW4Mk7C2HtvtWxi3kJoWbw2FsDG/E+lNTRHKpimlFprrIfV18+/cRSKPxCUnLsLk4u5xknRmkllTsqdLZn1m4ieXmZlIjkLfabjpuFS3lPSWxthVqmfcPSPdR9wY6f2Yh4vH9U/vk6s7zmIqa6F+SjWKrWfCuKZvqlvGuCbtRFJNQywMYcuAYtJ1JSXe3mrD63DpgikPOdUzk4mLfB1+j0v3LZfKKknd1DkpgdVNATF93jXfx3VL7F849jl9fre+kn9fK30rdUvqn+qZUrKrPca+b/6cEvrvS2E+tInyKf3YvL6z2p75fUj3SbNvi+7T4FMa+PFk+8RINco10z4LfY/Mh6hUi6S01tTfGVLP3JqfDh6vwh3Z/1zre9739t5/66ID997/jaQf0LvWrPu53vsTFx0PAAAAAAAAuGrcw7kP3vz5qhHG//E7/v3HRhgPAAAAAAAAuDLcw7nnbv688Lfm7nDnGK8eYTwAAAAAAADgynAP526vMffsEcZ/dPPnae/9jSOMBwAAAAAAAFwZ7uHcWzZ//qkRxv+Tmz/fOsJYAAAAAAAAwJXiMlj+maSPlPT5kr5lx/G/QOtwiV/acZy9aq01SR8u6VO0DsT445K+pff+NWP0H9vb54/pxtn1c/ukJBaXLJe4pLCc1jo8TS0lrbrkF5ccsx5jeApZSpSZu1Q/8/OSTzJbhCS0mNZqEq966CuX8hPTzUK76xrefhc2VUkxS/tRye5LIZIrs3PpddgUM4V02JC06oKCVjFpNbUPTzezfQvvdSukYSYxHdYNOzwIejO4GaOQbpamLJvWWkzks+lmlfS+lISWUs9cQlpKYLVJaKFvIWG4p2ukcOwrWrhYbfJhuEbcPBku32oX0RsAACAASURBVJKe9s20p3vR3LTPJv6Cmk1DqqpJ8TwKYxybtNbU1yWTHYWLfRLapy6ttTC5VNLRpPNSu01fc4IuU+1jTvBK3SL5WqJSz7iaSsp1hzu3FqHvYsd6Zrmopc/b2iXdR+w9bozJxTf7lMxwHqYkb1fPhPnXHaF4SzX3vpb2Ld2LzOSX7iOuvkiprKuQLummhpjW6mqUVCe5g1xMZd1b7ZLOLVe7xL7D92GMeqZUX1TqmXSNmOuhUhPlOmn4uPE9tXVn6GsHLrabayd9vlvZC6paWG+Lu2b+ItYzE5P6HU4il8oqSVNTj1TqmVQzuNolprWGxFcn1TNDapdn5rPB26lw784Pbf78iNbal1104NbaF0l6/uY/X3nRcfattfbBkn5a0q9I+puSPlrS7yqsuVftDwAAAAAAACRbjz1777/UWvtRSS+V9Hdaa/Pe+7dWBm2tfYHW37rrkv6tpH80wr6OrrX2UZJ+StJzJP2qpG+T9MO99yfG6A8AAAAAAACcJ32P8ssl/TtJU0nf3Fp7dWvt4+41WGvtI1trPy7pOzY/2yX9td77fKwdHktr7X20/kbfc7R+ePiC3vu3n/NgrtQfAAAAAAAAuBf7C8O9999urb1E0s9Iem9JL5L0otbab2q9fty/lfQbkk4k/VFJH6L1Gmwfuhni9i9ff3nv/VV72/vdvFzSH9L6m3Cf1btdAWeX/gAAAAAAAMC5wnKMUu/9X7bWPlzS90l68ab5D2/+SW4/lPt9SZ/Xe/+JUfZyZK21D5H0WZLmkr70Xg/aqv337W2nj+nayY13/ndlAWUnhkeY9rSAdmVh5bQQpVvIOfYN7W5R5LRAtF0UOYY5bLevwj6koAG3iG9lAeUxAiESO0RcjLa2sLJjX3b6cbMWqVmzVFI4xpKaWRs0LWi8cmu1pnHT+2cXRQ5j2L6VkIfh40ap7xiL7rsxxlhs2S76XdgHyZ7jcQFlt2BzcZHi0kLHru8Y4xYWUB7j/S/NW+Ea6XZt33BvcD+f7nFhe659HhZRn5rJ6CgEP0wLCyi7IAbJL6zsfl7yCxqnxY9b2J4doxjysC+uvkj1jDsHUkBWGsPVHSlwqlbPhLrKBDqkMdw5W6tnwnELQRG+nkmT9cC2qjQ/FRbid/WF5AMhEvey0/rubi5L9UVaML9Sz7h10dP5li5rHwgx/L3OoVWFMKyKEeqZ0seqSmhVoZ6p1D5SsWao1DMxKMI0xvCt4eP611EIstII08sIgRA54Mrdi4YXrz1ce8uprwOWUxMWFPq6uuMo9Q11R6WecX0rgRCVuuW89qHu/vnT0/gYbSfnjtp7f6uk/7S19mclfZmkT5L0yDk/8iZJ3ynpW3vvT462l+P7S1pPIa+U9NuttS+U9Gcl3ZD0/0j6x733X9+hPwAAAAAAAHBPgx759d5fK+m1rbWZpI/S+tdYn6P1g7onJL1N0ut672/e146O7PY3Ad8s6ZclPe+Ov3uppK9urf213vt3XrD/YK21p4Z0q44LAAAAAACAvfqNlr7Ot9F7f/xeg6QveqcB5733X+i9/0Dv/Zt673+z9/6/9d5/9AF6MCdJf3rz538v6e2SPkHS45L+E0k/Lem6pO9orX3EBfsDAAAAAAAA97SfX5bd2Pw67JdK6r33l+5zW0O11h7Vu341919KelHv/WTz37/UWvtkrUMvPkLS32qtfVqlv6T/rLI/Q56gbpJif7cyLgAAAAAAAPbqg3vvb9t1kNI35y7gAyV96uaf+8V73PHvL7vjQZskaRP28I2b//xYSe9Z6d+aW9EcAAAAAAAA2LbXb87dp96mdVjkRNJvhT6v3/z5Hpt+lf7PlfTWUfY0ePvNRzWbnZfLkVWSXW26WQpdKiSkpfRUu72YphcSxMz+xRQyM0buOzyZLCWG2v6pbyHtsxQrlX4V3iQ3JTFhyab6hWNUSNp0SWE9pZjF9FvXt5DAWkgxk8Jhjn1N58J7Wg4f2lfQ4hhJaC5lL/3/jkoKbCXFNaabFcattKd9LqW1br+plUQ3yR+LmCzn0gnjvcG3u+uvMpXFAEDzwlPSW0rWXpoE1klIJmtm7ElIZY3JYnaMsD3382Fct72cbmabY38/xn4ml13rlvUYw/vmMcw5G/bD1SixviikCaekPre9Hu45bow0bqxRCjWRvdemeWF4sGcprTWlsqYd8Qmsw2uGlNbqPnXFOqLQnuqZpetbrGcqabuVemaURPldFb9OsWuKa6xnXH1RrGcqKcW2vVrPVLZXSnY1b3Z6zZX2ygmXztk0p7qxQ9K1K156+HC9Mu2rNA+FdO6lqS9aSnh3fdMpWxojpLUWahTXfui6ZX5zP99x2/c35+47vfe5pN/c/OcfCt0Wd/z7k8X+z1x87wAAAAAAAPAweegezm38+ubPjw5//8c2f76l9/6OC/QHAAAAAAAA7ulhfTj39zZ//nettRvm7z9/8+dPXrA/AAAAAAAAcE8P68O5V0n6v7VeH+7HW2vvKUlt7Su1Tly9JenrLtgfAAAAAAAAuKeHMRBCvffeWvscST8r6cWS3t5a+1Wt02XfR+t15j679/47F+m/b++4eUPT6cUCIdIi2sMHqC1oXNqcW2w5LQ4cwxGGjRvHiKs779g39S8slp4CBUoqmRKFheolvwDyKizE7gMvwriVY19Z6LiyIHXxuqmM0Uqr4Bd2Yl8LKFftvIDy7tsaZYxKmEPaXmnRZ3eRpL7DtnXe9nzAxggnUZoPzQnawgu094G4gH0hHCMG5JhhwwrvdrHkMO/FRYrd4vPxfXLHLXStvH+FvvdLJn2pnhkhVMJuLtZEQwfI23PneHzNboy4wP+OfaVaaJUbO84LBYWAq3jc0suzIQ+VeqYQglAJlorbG943qWwvj1Go13bc1l6NUDM4tTCswj5UxyiEOVSCV1L4Uql+KtRJY3ymKUnXqhk8HQufqJUCBit1YCXcpvABKIUuxDpn4LhBqZ4ZoUbdpZ5Z3tzPhPWwfnNOvfc3SvpwSd8g6c2S/sTmr35Y0kf33n9il/4AAAAAAADAvWx9c6619m0jjv+8EccaXe/99yR9xeaf0fsDAAAAAAAA53G/1vrFun++WAwAAAAAAABcWWnNuftkRREAAAAAAADg6nIP51546XsBAAAAAAAAPIS2Hs713l97iB3BcGdPHWuqa/vf0Bi/3LyndNi97UMlUbP62iopoGZ7YyRppcPpUnfioS+kt6Vk11JybekFFtvd9vaVLDbC9TRGeOb9bOfLfYR0s2SU9K9dp5Exkt5KCZ7Du+brN3S3Kcwh/Wvp4sbCwC56Kya67X4s9hW2vLf3aZ9j3M8uu54Zo3ap7MO+EjVH2N6+6plqaqHtG94nF86ck3IHN/pdHqGW3FvdcpH+d6FuGeCy02Er445xj9uxfirXMzum38b045SqWklQLh2LEeqZSt05ynu964Sx24+PNsYAy6fmexn3oU1rBQAAAAAAAA6Nh3MAAAAAAADAgaRAiHtqrX3XwK5zSV/ce/xCOAAAAAAAAPBQ2no411r7eEnfIOkf9d6/8Zyf/TwNX4ngOyS9rrx3AAAAAAAAwBXmvjn3yZKeL2kl6byHc7fda9m9LumTxMO50UyfmmoaV9rHLva1hvN9Y4RFzdMz+Uv9auxVf59wNV3l748XF3gPy7CPsy8DtbQXV/l9wkPjStUzOy7wnlHPABdyVe6TY4QGFlr3hXrmAJ7ez7MY93DuBVqffq8cOMafk/Tm8HffKOlTJX1kec8AAAAAAACAK849nHvu5s83Dhzjt3vv9uFca+1HtX449/4X2DcAAAAAAADgSnNpre+3+fOtI4z/W3eNCQAAAAAAAGDDPZx71ubPPxhh/Ldv/uSbcwAAAAAAAMBd3K+1vl3rb7o9+x4/+yV39E+ub/5kuVMAAAAAAADgLu7h3Nu0fjj3wZJ+Pv1g7/1bB4z/AZs/f7++a0iOnpzoaOm+9DiyyiPVQt9SwlaImamNMbw9juvaU/JpZT8q6akpcmeEBNbuTqfi9mqvz4xdGTepxBKVXkcx7mhP147/+QcwimmMmD33squHonSRFLa367ihvaVxS/uW2rfHjqdW5TKrjBFeX1sN3wc7wqW/juIYdtwRktfGmBoetOlljP8VPEY9U3ijRqk7dr3njFDP2DpC8sdiX68jjXGV6pnCvpXGjf133N5Y+3E/2LV22ef87fZthHvflalnXB0RlGufXeunPdYzdj/2WhMNP1/2Vs9c0tTSbu7nWYwb9XVanw+fOML4f2bz52+d2wsAAAAAAAB4CLmHc6/e/PnS1tp7XnTg1tqxpM/Q+vnlP7/oOAAAAAAAAMBV5R7O/ZTWYRCPS/rqHcb+q5L+8Obff3KHcQAAAAAAAIAraevhXO/9SUnfqPWvtn5Za+2zq4O21l4s6eu1/tbcL/fef3bH/QQAAAAAAACunLSS3cslvV7rB3Tf21p7WWvtxr0Ga61NW2tfLumfaB020SX99bF2FgAAAAAAALhKXFqreu83W2ufKukXJL2XpK+U9IWtte/Uek26X5P0hNYP354t6cO0DpD4AkkfpHcFjPwvvffX7PUVPISOn5SOFhf84R2TsMpBRZUxzKPi3kLnUkpX6Ou2F9PGhvftkxATs+P2NPEvJG6vkrDlInNGeH0xisf1jeNut7fKaw7946ll9tm1rf/Ct7ux0xi1vrtHEKXXvate2LVemEhSX7e9St/1X2z3X8UEse32HtK/4utbuTFSXzNu7Ds8pa3FMbZ/IG3PhhOmY5zGMK8vjzHs5+O4hWMc96OwvSTu845papVUuPIUcj8kpF12cvwI9YWbaOPUUgh8i31t/VToG+6puc7Zbkvnlqvj0rjp3m6DKOPrG6GesTXY8Bolbc++vjBurAMKtdakMFnXaqLhfZNa38FdSyp1Sx4j1R3mvE/XSKFvOvH9GIX6Kd0n0z67/oUapZuaY/0XblzftVLPxJOokp5aqCXixwb3WgqvL9Yche1Vxqgm15bqp0rtM/Tni2PsUrdMbl38Z89jH85JUu/9X7fWPkbSKyX9cUnvI+l/3vxzniZpKekreu9/d6wdBQAAAAAAAK6ac/9/Xe/9X0t6vtbBEH+g9YO3e/3zf0p6AQ/mAAAAAAAAgPPFb87d1nu/KellrbW/K+lFkv68pOdJel9JM0lvk/Q7kl4r6ad772/a3+4CAAAAAAAAV8c9H87d1nt/WutfcX3l/nYHAAAAAAAAeHgMfjiH+8fxU12z+cVWMCwFOlQW4K0sdFwKaEghCIXthb6r6W5988LMYZ/NGKtpeB/d9sLm8kq5ZoH3ygLRad/Sosiufxijmfa0EPTEtE+mfnVR13fdvt1/GsaYmhVG3c9L8a3W1PRPix+7BZvdPiR2wedztrcvlZCHFLrg2mNGgOm7XPmLchkW8XX7nMfYbl+FcWP70oyRFmxemmNhfl4KCzObn1+PYZvtvBzf0cJ7XclSaWHfXHsLr8/3Le6bGyMtaFzpmxZhtgs2h+vaHbfKAtG1bJtxQiX2pBQAVQpzCGMUQhds7VKpW9L2XN1S7etqu+nwumXdbuqL0NfVM3GR+TBEJaDB1i6pXos12PAapZlawtU4kr9f1+sZUzOEGqVSi0zD9twYse4wbalvarfj3if1TAyMMpaFwKlKLZJyFFamv9sHSVq6WiRtr1LPpFrLBUKEe7itXVI9EzMDhwds2CFiqFcYw93DQ90xMa8l9a2EYU0KY1T6VusZV7uUAij2GFbh7DK1TE72My8VMqIAAAAAAAAAjImHcwAAAAAAAMCB8HAOAAAAAAAAOBAezgEAAAAAAAAHwsM5AAAAAAAA4EBIa30AXXtypdlpijPZKKSyVhJYY3pqKYE19TXpOiHxKo5hE1F934lJJ0tpY81cKanvKl1VLkAsJULZWJrh48bmEEvTXUpXSmU9CildR8MTy1xSakosOzLjpmSyo6mPIJqZsY/CGDMzRuwbIo9c/0l4o47MGCnF7MjEGFUSz9b97zF33MMqxhSn/iaxLJzMru8iJpZtty/Cvs2X/mJ1/VNC2ty0x3FDqqp7LYuFH8OlqS2X/r1eLbb79pjinN4/079yqsRY3bA1l3K68Ps8WQz7+dTufr4+Rpj3dkx2lXxyWkxr3TVNLUTWjZLWuq9AxRESWH3f3dPgK31TLTJGAuvqaHg9Y+ukULf0cG659MWUImmT41PXyu2lMsWlVFZTX0i+dpmEvq6eSWnwR0fD64ujMIarUWaFtFZXR6RxU/9UR9jaJ6W1mrPgQaxnUt3h6xl/UbqaKNU+8zCG6x9rFNs31C1hDFcrLUwtIvlk1+UkpM+7z4ItvE/h3u5uc/lUMempxVRzXzPsnihfqX0qdU4LtaStRVKdVEqwH6Gese9pLfa7NL0M6Du517OYC+KbcwAAAAAAAMCB8HAOAAAAAAAAOBAezgEAAAAAAAAHwsM5AAAAAAAA4EAIhHgAHb9jruPj+Tv/Oy/8bbiQh7iorlmUMwVClEIeQl8T0JAXUA4LytuFjkNfE2zQwxXRlmbfCsEPkl+gNC0l2dxxDgtfltazrSygHMZNIQ8TczwnKaBhZhZFDosfH5sFlI+P/Mqnx2F716bb/VPfY7Oq6nEIfnDjpjFmcQHl7bGn4cxwiy2nvpWFlaeFRZWXhcWPJWlpTqTU17XPw4SxMO2nS39RzsM+n5mL+CwsfnyynG2PGyao04XfDzf2WVjI2y3kPJ+HhZltW1hsOQQC2NOocm8J72kcwS0EPELIw2S+3dZS3xDyYBdhjmOYtrTYclywebt/uiTdPsdFke1iyynYohAUUVhUuaXzLdhXPeNqhlTPpHuf65/qCxeyFOuWUEusXE0U+rr7bxzX7HM631bb09663bSlt85Nvy3VdrEgLYRkuZM21C0xtMoFUZlaRPL1TApXsPVMqltCneOCqK5PzcQXxnb1iZRrFFfnHIXJembaYyCEaa/WM2PULk6lnknBDb6vP/FdjeLqEynXOWdm7DSGq0XiuKEmOjN1ztmkEJIV6pmVmRxW4T7igmnWf2HaCreWGCgQw5fMnFoJeSjUM7H2SaFVrkYpjBE+/hRDJSqBWpXwiFT7XE490878vLsrvjkHAAAAAAAAHAgP5wAAAAAAAIAD4eEcAAAAAAAAcCA8nAMAAAAAAAAOhIdzAAAAAAAAwIGQ1voAOnrqTEez0/M7FRLPcmJZJa01tB9tP//NiWUu2TX0nQ0fIyeWDR/XJsLEVLGwPTduGKKZx+Yp5baSNFMS081CmpZJnXQpZpJPLLs28/FB101iWUpJvXHkk3Oum/YbId3Mtad0s0cmZ7bdpaldC3FMM5vsOjwJrZxuFjOCd5PSQX0Ca0pVNQml4cR37SchWvA0pZ6Z9meWx7bvMyat9SSkm51M/X48s9hun078GG2+24WdfnoV0s1ssnZK/DUTVyWVVaqlm9lEr0q6WTiWrq8UEljTGC7dLO7b8KTUyTxc1zbZNaW0uYjwtA9+ez7drBBJPoZ4o3TnrO/ap2Z+inXL8GvEJapK0mpWqH1iPeMSWIePEa8ncw6kdMp0PN1bEt8ml5SbUvZSiqtN5g1juHRYk54rSROTyir5ZNZcz2xf8NcL9UysW0KNct3UP7me2a5RroUJKrVfNxNlpUbJfbe3dz/XLal/pUaZp1rE1EQpDf7WytcoLpnV1S2p/cTUJ5KvWyTp1NT8qZ45nQ9/3DA3xz59/olprcva57PtgWvttkZJ86+tL3xf1z6t1jOmfymtNaXaV+qZlCTr6plFqEXMPaOc1pr6OwPqmdXcfwbcFd+cAwAAAAAAAA6Eh3MAAAAAAADAgfBwDgAAAAAAADgQHs4BAAAAAAAAB8LDOQAAAAAAAOBASGt9AE2eOtHEpKC+m0K6WUos02R42pjrK4W01rDvK9fXJJ5J0mqRxjBpaschdcmlBcZwFjfG8CRDSfZReEogaiasxrVJ56QkVjq7Fx4ORgsprlOTenYUkl1nJgnNpZhJ0iOz7TScR458Qk5qf8y03whJqy71LKWyPjL1qcnX2/YYLvEs9a2km03DiTEN5+dkT6lnq/D/es7MSV7pmxJYXeLrrdU12/ckpMPeWm73vxVS71y62c2F396RSTGTfBJdSqer6CbdLKayhjm8LwrxiwUx8LUyx5nLIaWNuYTSShKa5NPQpmfDE8tyOmw4L0w6WQt9XZJZTFp1qWcmHU2S2iocfDf2fZzWKpfKKqm7GiVdC6FGcYmvqRZztYurcaRz6hmTwBqmQ/+2FhLlU/LpJCXa2kRUP7a91sc4VQr1TCVlXpKOTD3jUlkln8z6yCykp5pk1lS3PGqSViXp0aPtuiOltbraJdctKZV+u7+rRSRf5xyHesbVIpW+Y0i1yDKcXEtz4lfS50+6v4BdnZP6urpF8imuz0xCAut0u+/TE58CezTx78l04fs7Lv02JeKuzLWX+iqkgPrPNMPnwzi1xPT54X3d4Uz1jKtFYt1SqFFi4qsZI6e1hrrDtLsaRwr1TOgrU6O4tNfUd91u+u9Qz0wWJ/fscxF8cw4AAAAAAAA4EB7OAQAAAAAAAAfCwzkAAAAAAADgQHg4BwAAAAAAABwIgRAPoHbzGbUQJHBPblHkEObgFlZuMfjB71Az7Wmx5TYzC8cv/SnqFv2WpHbsUhfC6yssSOzWIu1hcdHJxA/sFmIPa9/69zcu7lx4IXGVU9MUXkdawH5qFlw+mvoXeGzar03DAsMuoCEsoPwss1Cy5BdQfnzqF/J0Cyg/Fvo+ahZKllIghN9n29e0ST4oIgU/xECIUVbi3pYW7D0z/w9oFa5JtwDy2cTPLa5vDoTwiyLfNO/J08vrtu9Tbbs9HcvKMU7HbbHaPkZL0yZJC7No/yTM1WEI/7/qCnkQ8SUX2tMYbgHlNHe6hZUriy1LfmFlt1Dyuq859qZNOicQ4mz7xbR5eIE2PCK8wKUZN4RHxAWUK4EQaYxdpRqlEAjhahcX8CBJbernnD4z9UgMhNgeox2HOimGdBTqmQpz3GKZFGpNd02GdfH9HFBYxz1K81MhrCLVa1MTFDELoRKVesbVLo8X6hZJeswEOqR6xtUulbpF8oEQqa8LdEh9J2bF/OMQ/HDZdUsMhDDtLvhB8nVHDoTYDleo1C2Sr11upZAHVzNUPkvIH7tYo9h6JoX1mb4pNCd93nJD755vFT+H2UCI9PnOBkKEe4CrZ2Lg1OXWMym4YXK2vdOpnrFjLFJf016pW1L/HeqWtnzmwj97Hr45BwAAAAAAABwID+cAAAAAAACAA+HhHAAAAAAAAHAgPJwDAAAAAAAADoSHcwAAAAAAAMCBkNb6AOq3bqmHhJp3moRYGhdZFRLLZBLLUlprOwqnkktxdYlnCollIUSlX7toXO2dGzTpbeG4uSSzlGLWwz67hJ5KwmE5DbGikNYaU1wL6WYzl252lNJat9sfC2mtN6a+3SWZpQTWxyfb7Y9PfSJPSj17tG3vR0osu27imK6FmKdjEwk1C5d6+j8vu145IYAqXapyl/U8xGad9O1jNA8xgjf7dgpZOsa3uk9xnZi9nsa0uOGJTjGB1bwWl2KW2udL/+65ZMGUQhgnkj2l4SU23SwdYpuEFuYnm+wa5qyQkOZSz3LSqks381dJam+FdDOZZFabYib51LOFn2d7TDczY5gU2GhVPK9S7eK4VNWQ7txc+nylbpH86w71jEuAm6SU21RMuK4xcd3VM77rymyuLUMKYbh2ZJIW0/VbmDrH4eaylKxdSGs9Mm2ST2Z1KfOSdMO0p7rFpbJK0nsc3dpqc3WL5JNWn5WS6ltIcZ1s77OrcSRpZt7s66GemZnidXqf1zMuEPMsXGinpp45CZHGJ+Z4uhpHkmYuwlO+dhmjblmFem1hJpJFmMvOXN9Q+8zNvaj6eSTOk0NVy6TCZzZbo4S3yaVip/T52G5OWpfKKoV6JiWtnqY6x+xISJS3SfOxntnu22NaaxjDps+n2ufetUtfkdYKAAAAAAAAXCk8nAMAAAAAAAAOhIdzAAAAAAAAwIHwcA4AAAAAAAA4EAIhHkCrk1Ot2h2rXRYWUG5uAeUWFgI2fXth8WNJan1m2sICnqYtPT0O64iqmdfSwgqzbqHjtED4yrzssMatXShZkl8wtLJQ8j7Xay8sBt9CX/eq3aLKknRsVjk9DquZusWWU99HwsLKj0y2212b5EMeUvBDWoTZB0L4E8a1Xwun0LE5v2dhsd5JuK6npv8kXGmruCzytmU4QedmsdV5mANcOMZJWNl3ahZrdeeVpLy6c2E16blZyHkezqG5m2flF0s+Xfo51S0+nhYkn5pFkdNC53E9efMXo0w5lSCbtIByITzCtafTIs3hlXtDW7gFlMP9MC2s7NrTAsqnZqH5FPIwN31D8EMPY9iFlcP9vqfAgx25+7okyQVUhWtPLvwhvOY2265bJKmZazW9YrvH4XWkudoFVE1CPdOn23uSQh7c9ZDW1S6FPIxwrY/BHc407+VAiO32NP+69nQv8oEQPjwi1SjuPumCHyQf/pCCH54V6hm/vRDgZepDfzVJM1fPuNA6+bolqdQzqW5Zhblsbl7fWVhQ/ti0z0LfWd8+XyZpofpgZYIplulYmL6pbpmHkAdXm7vgB0k6nmzPnWcpvMe1VQOr9hVwVRg2fr5zY8QgQTOvh6CCWKO4MVKglqlnUt1igx8kH1p15uc4F/5g6xbJ369DndRTIIS5rmNf565jv+p+Lt0V35wDAAAAAAAADoSHcwAAAAAAAMCB8HAOAAAAAAAAOBAezgEAAAAAAAAHwsM5AAAAAAAA4EBIa30A9cVCvYU0k9tC4pFNXkvpZibVZHhe0mYI0xbTAk3yWp/61+ESZSSpTU3SzCKktR6ZZNeQgtNMAmsKA4rpZi79KSRR7itoqCS8UTE4z+x0ToY0aYjhRR+ZHkwEWQAAIABJREFUAzoLMYuV9usmUVWSrk+2r63r4Xo71vAE1pTW+og5SC7FbD3G9rWakslmpu+6f/UqfnerlG4W4qZmpv88HItTm07mt7eMccnbzsx7KkkrE/t80ny23DUzxolJcJWk2SqkW5nXchRS/dx5n9KP3fUUU5Xvi8klSHOqS7qu9K3MyaF/TDdz7Sm9PNy3XGJZM23rvtvpZL3SN6Sb2dpAUnfpriFFsLv7ZzFx0NYuKZHe9G0pSda0t5A+30P9ZFMEQ43SXXK8S5eV1I/C+bI0ry+dh+6cjed3oZ6ptN/HU0t1PnTtaf51c3Wa1909INYtKcHe3IuOYxq8qWfC/dClskrSNVvPpMT17XPrWrieZiYufZoSjUOds696Zh5qu5mpcyah79SOPbxumZsEVynXo3OXiBqSVk9NnVOtq915n+p4W6PYnuGzRPwwEdovWWU+dPNyrlFMW0zQLqS1ps/Q5n4f65ZY55jzJdQoNpk1Jceb2iWlzMcEVtNu65b1X/j2O7v0sK874ptzAAAAAAAAwIHwcA4AAAAAAAA4EB7OAQAAAAAAAAfCwzkAAAAAAADgQHg4BwAAAAAAABwIaa0Pot51z1iskPLTTTpSS+lBJh4nJaC0lPhq09tSMopLvSv0lU+rGSNtbJRksktOMnP7vM9dqKRAujSmlPI0KSRCuSS0NMY0jmH6hnFTipVPZLNdbTrZLKSbucSyairrNCU5DxSudC1D8rDMMVqFvm7P0t4em/d0HnpPQ+qSe6/duSJJ08J5MQ1jpLF3la6HvbGTyyXHpqWXXDkUcYzhg9hEzPQ2V+59LiU19Q3pwJX7b4/7tr0fMQmtcNwiU7v0lGpuJqPYt1JfpOPZzfyS3qejQu1TTBO2Cmnw/ueL7fuStnfJc9y+5lRfi6R7Tri/FO5bvp5J9ye/PXevTfdl1z6Ntch2+1GoMHatW6JQG6S6am4O0SwcN5dg79LrJenMXL+pvpyGecG+15W6OtbPI9TmhevpshPlLz3AfscaJe7vGGn3lc/slc/nlTHCuLbuCImxMX3ebm94gvL2gHu6T+xlVAAAAAAAAAD3xMM5AAAAAAAA4EB4OAcAAAAAAAAcCA/nAAAAAAAAgAMhEOKqMgut5r6FZ7ST4vPcidmPyr5VuRCLyuZCXzvGJa9/XnXZ67NfphQoUBEDDCpjhJNgZdqXYd3QpVnFdRVWdl2ZBXjdtiRpEqMbdrMMCyinfV6aBVPday7vh3nd1fNitaf/P7V0C8dfJe44X/Z8c5/Mb93cc+KupXufa3f3ztS3dL8f4cDFmsHMDdUFk+3rG+F6GuMYXeb7JH9uXbr75bq2c/tlr+S+H/u8X1Tucen+ubTn4fBjn+73EzM3xOCAUHfsKyjChTmk9lWY42ztEzNhTM04wsU3Rp07Rr19P7uvX57Zt7i/lfyfOMaeDkYaN90/bV8XyJQCLdO8YPqnfdtT2MMQV/wTBAAAAAAAAHD/4uEcAAAAAAAAcCA8nAMAAAAAAAAOhIdzAAAAAAAAwIHwcA4AAAAAAAA4ENJaH0DtaKrW7vHWpaQSk4zSUlLJbLbddxrGDe3NpauEvt21p+2F1Ng+HZ6Q5vrGdLQxEnMewMRXJwXYdHNAUtaNS39KiVCL1XbqaOo77z6h1LWnFLN53762zsK4s7Q9k+g1N0mrkjQ1R8m1SYqJZU5KSJvu6aRzyWSSNDfpSPPwOk7MGKfhUMxNwt1Z4T2V/Pua+rpzKJ1v6fxcmX12bZJPb8vjmmvvvo4gq9k1cbs6V3ebMl4Yw92HlO8v9j6ZktGn5pxbhXlhtn3xxEOZJna3z8uwvX2dc6GesfVIOG7NHbdZqKNc39Se3ifTHuuLcL5UaobKOTtG+vxVmV7SPFmpUfxcHeZ1ew/YvZ4Z4x531nzy4czs86lLPZTkq74wXxRCYKfh2ln0tB/DpLolpbW62mUednpu29J7vX2M0/uf6lH3vuY6d/d6plKb2xrF9gyfJR7EYOYRPgtW5up0f9n1c7GtT6TSZ/m2CvflpblPpucJ5hrp2n5OIUnNXn3y6bCri3/ean0hLe7ZrYxvzgEAAAAAAAAHwsM5AAAAAAAA4EB4OAcAAAAAAAAcCA/nAAAAAAAAgAMhEOIB1I6P1drxPTqFFSPNwo455MEtfuwXDG1HYQHlY7NYY1ps2SzO3I/CgpOp3byWfhQWvjQLQ4b1UH1ft7Ck8kLJdow9Lti8s0Lwg1RbQHlpFvZchAVD3UK5pyu/COjpyq/M6frfbNds34lZCHgawhzSOsf2f3uEvkuzCPMybG/Wtt+U4/BGTcKCzdMdV/JOQRPpULgFl8NSrTox++bef0m6aRY/vrXy7+nN7ufLk759XpyEc6vS93Tlb62ufRFen7seluEaWa6GL6Ccrt9LX3C5Mse5kIdw23Ltqe8qLMQ/mWwfjNTX3T/7UTiYs3D/NAe/hzekdGswdUAMpUj35aWZR1ybtL+TqLLPMczB1D5HoQRO7a5GSaESpibq4f2v1DPpPHT3nNS3co2k9l2v331yp2F13nNzapp/3Vyd5nV3D0j3i3R/uTbZbp+tQphD266Jpmmh80o9EyzNIKtQM8xMmEOqtcyUvO6/48lVrWfOXD0T9u3E1a7hQ8ZNU1/cDPXMrR7qnNV2nZPGcDVxrlv8eehqs1THu5p/lICr+6Fukf/MVvss6F+I/+wZxp0Wxkifi837l+5PSiEPLnAmhS6Y0jxd0d0EQ7VpqEXScw0XZpUCtQbUM62vCIQAAAAAAAAArhIezgEAAAAAAAAHwsM5AAAAAAAA4EB4OAcAAAAAAAAcCA/nAAAAAAAAgAMhrfUB1K5fU5tcv0enkHdi0k4UUkebSXaNKWYhrbW79pBu5pLMUrrZ6ji0z0wikGmTpJVJq3FtkrQym8uJOb69FrN3yQrpSCnAxiUspXSz+XL7IJ1N/XlxZlLITkMyWUqbumVSrCYhFWwac7q2Ld01Imll/r9H6ntsEstOTcKaJM3MPrs2SZrGFNfd4q1iKms4mV1K1zz8fyGX/nXikp/k01NjullKcTXtqe+tpevrU2BzmrA57821IElz03eZ0s3MdebapHPCzXZM8Y0uOYHV9k0p3CHNcjUz7SmZzF0QYZKchMSylSnFWryHb7fHvpWE0oWfc1xqbLtP0lpdMms8Fu51V5Lq5ZNZU1qrr2dSLRLaK/WMOWfT+e2uh2paa+WavPTax9UzoWuaJ117mn/dXJ3m9dPJdnu6X9wKKYIzc/2l+72zCm9UqlHcvX3efOb6dVO7nKT0eVNNTENq5f1Sz8zNxZMS5V3tciZ/Xri6I9UzLpU1jZFqFNdeSZmXpDPT7tokX/Onzwc+rdV2PSfF9T74cFWpfVKNYlO4w7Dh82szt+sWjr29nMLOxW93uXtwui/Pt+cL++xBUluYF5LqljB3NleDpfTqIWmtqy7dume3Mr45BwAAAAAAABwID+cAAAAAAACAA+HhHAAAAAAAAHAgPJwDAAAAAAAADoSHcwAAAAAAAMCBkNb6AGrXr6tNb9zRUEilccmsIRmluySz0DelnnWTepYSWF2SWUxrTYllxy7dbHgiX0prtelmhXQdqZj46pLQCsk/UQqfsQk9oWtMN9vekcUypHKa8+V06aejI5P05dqknMBa4ZLM5iExNLdvv9kuXVSSrpvUs1lMa3UpbSEJzUU0BZWkt5RiFhPgXNJbOG5nhePm2k9i6t3wFNeUbvbUcjsh++YipKkt/RgnS7PPpk3y6WauTfIpgu56lKQe2ncMvcsqYWqpr01gTXP19gtZhb7N9JV86lm6nNJ57wdO+7F9DU/CPsu1h7R0uXSzlGIWzsPm+ufI7uF9E3eMQqK87ztCPRPT5ys1iql9Cqmskq9nekqUd+nz4bRw10M10die9iNc66OwgXwp6TqljJv7VkrWnmxfIycTfz0dmb4pobRiVUoX9a/D3X8lf691dUtqPw6Tp6tzYiprMZXeSbWLk47R0nynJfV1dU6sZ0ztUk2fd7VLqmeeNrXL0yaRXpKeCfeGk8V2e6rj5+ZzQ/p84K7J9LljjM80VjEEtpJe7eblll5e4bNn+vzqLp0WP1COUM+Y+3VMYHXtR6HYWpoXEmr+5vpKvv5xdUtyVz3TlvspnPnmHAAAAAAAAHAgPJwDAAAAAAAADoSHcwAAAAAAAMCB8HAOAAAAAAAAOBACIR5A/ZEb6kc3zu9UWUA5LOpoFwKOiy2nMdyCxmEMsyjyKgRN5JCHXQMhbFctC33Tgs3uUXhcS9wuoOwXnhxlOcrCAso9LCS6NIu7LsL75xZWPg2vb2J2Li0OvAr75kMewiK+JsQkLe58MvUngVuw99okLaB8ttWWF1A2gRApHCMERaT+Qy0rC+DLL4qcFma2CyiHBV8roRu3wkLHp2Z7t0KYw00zRloo+ebCj3HLtKcFlE/NNZIWUHbXXrx+K4EQlXVyiwsou1OgsthyWuy+ucMZTvm8KPJwfiHocI8Ll44Lf0j31GaCBiaL8AJNuw14kNTDwsY7B0KMoVDPpGPvgjRi3xTc4AIhUu1j+qYFu13wgxRqlFDP2Bol1T7mGkn1TLrOdg2tKs8XtvPw9up8aOuZcL6cTrYP0nTiD+hkT+EPqfax9Uy4p55MQ2jVyoVW+RrluqlzUsCVC3NIfZN91TOr8N0V178SsHHSfW1wat6TkxCcFesZcxHnkIft/Uh1y0moUZ5xgRCLEAjh6pkQ8mDrmXB+7yvgKv54un7tvSiM4ebOdG+wIUvpNYd7uO0fahTz+nJd5v9i4gK1jkIQlfnM3kLwjg1fCPWMUlCDO0Y71DN9QSAEAAAAAAAAcKXwcA4AAAAAAAA4EB7OAQAAAAAAAAfCwzkAAAAAAADgQHg4BwAAAAAAABwIaa0PoNVj17SaXb/YD7tEmZTE4pI2QzpLSm1xCTSpr0syiwk2hfZKYplLPFvvR2EfQrqZ215K8/EJgL5vSUw3M68lppv5IVYu3WwRElFNItAYKWaLkJrl0s0W4YC6xKsbU5+0emvl0618YplPN3MprqnvGOlmYxxnJ6XFLc3/A0rpZpVUXZc4596789u3x0gJrG6MWyatTJJO0hgm9ewkpJudmWvHJZ5JPt0svR/purZzwJ7CN6Uwx6VkSBcYmuYhNy+nNK7IpZr7nvawpfthSPZcmcSvSUjmbaZvXw5PLGshgSy1u8SyOIWUj/NAKVHehtANr1HS++SSXSWfHh9rFJcOW65nTFusZ7bbwzQkF9aZ9iFek6Y99t3X1wBGqGfifctcf3OTyipJbcd6Ju1DSrOcH5m0+5SoOd2uUVJyfEr2vDbZrjFcjSP52mUSElVt+nyI1r7suiWppc+b96mSPh/6pnrG1S5nhb4uTV7KCawurfUs1ChzU8+kzwfus0S63192PTPGZzb7lqTX55JPw302nYfuhh1rFFuXhWTXI78fq8X22JOQgO7qDlfjpL4pUbWFmshOIzvULau5/6y2K745BwAAAAAAABwID+cAAAAAAACAA+HhHAAAAAAAAHAgPJwDAAAAAAAADoRAiAfQ4tGZJsd+4c53Kqxx2tNiy+7RbTUQwi0aXAmESOEKYcFm17+02HK4InzQhO/rwiPSvsWFkl17ek9r69kOF9bI7GGhcrd28dIsSC5JZy28sW5cs7BrXEA5HFC3sPLJxL+Bx9PtxY9vTlKYgw9juGbGOAohD25R5KOJX8zULZacFkqehkWYJ3ta5d+FbkjS0rwn6f1ziyKnQIiFSV6Zh/c/LYrsFtE+C4kuLuQhBTSkxblPTX8X/LBu3x4jLaC8NO19mQJdUlDEdlML71Npbe64gPKuixSHcc2+9fA6VulaqIQ8mH1L95EwXWhi3qtVmDvdNNJSeEQhzKGF7blDlBak3luASFzzuxDc4QIaUt9CoEcOTDDvaQyLutx6xgZ1FeoWqbYAurvWq18N8At5pznOzC0p4CrMk0ttv/BQKpe4e98yBD+kdnc/ux5Cq9zC/64+kaTjMEHNTC1xFGoiV8+kgCtXu9zPdYsU3r9wMru6M9YzLjwi1RHhwna1y1kYw/WN20uhVaaeccEPkq9zXOiKJK3McUufOy77nhPbK/Oh+yyYPnsWdiKGprggyEKwRbrnTAo1yqQQ8hCmi1ogRAjYiLWLM6Dr4uxs+HgFfHMOAAAAAAAAOBAezgEAAAAAAAAHwsM5AAAAAAAA4EB4OAcAAAAAAAAcCA/nAAAAAAAAgAMhrfUBtHj0SO3aBd86m0JX6FtJdg3tYyTyVdJa8ximb0pCq6TrpNSzSuKgTTL0fUtprSl9xiXbpHSzkIKzcqlnbffn/0ubbhbSWkO62elk+02ZTX0k0LFJIUt9j0IkkEt8TX1dwtIsJKG5xLKU0DRJcUV7soon6LaUquvea5fKKvmUNZcett5eSnw1aWphey6ZdZ7S9EK62cIkjqXEV5dktlykdDOTxpXSzdJp4ZIPK0lo4TzsYYKyt5JCutkYZ3cr3Yv863P3onD5xtQzl0LmElzXfYf9fOobA2pHGOPSFeoZW3fEtNbhY6S+tmYopNrHMVLd4ZJkCwmsOal++Bjxf/dXppZKFHTq6uaycH7H9Hlz4qf5146bEqLNPWMxDfecI7+92XJ7306n/o1ytURKg0+J8sem/knpqS7FdRre01QTOfdLPePSWlPiq6slUlK9q4lS0mqqn1wtkdLnXd9UP8/DNeLGSAmsLmk+prWae19OmR+ePl+rZ3xzKZE6fU51oaMp4d20xSsh1jPm82QhOTylp6Y6xybKh/fJprWmpFVX+xST48euZxan+3mMxjfnAAAAAAAAgAPh4RwAAAAAAABwIDycAwAAAAAAAA6Eh3MAAAAAAADAgRAI8QCaPzKRrt/xXLUSCGCE9UlDIEQYIyxEWVuwudA3Bim4voWQhzRupW9aINpcbaUxdnyfz9PMG5uCHxQWAXVhIau4wvvwBV9XR9s74hZVlqS5WShZ8gsgH4VFmE9M37igcVhYeeq2F1Y5bWbsNG4Kf7B9C6ucVhZbrgQ/SH6x5LQosl1sOfRdmv1ICxovU7sNGxm+KHI6D+N+uJCHuCiy215Y8NcsVN5DoIBSuztdxlj4v7Cw8iqELkziTcqMUdiHlFfj2tPxbGZeXy0L4QqSmhm7tijy8MWW40LJ8ebu+vquQ3/+XIX7nN3l9DIqfdPC2rvWM8WawdUulZqoUs+EdeOLY4TQhX3VMzEQwrSF6zdPLWauTptz961Uz5j2xcpf7C5ASJImpj6YTYfXIqmecX0lX49MCwFXObTK9C1OGPuqXVLIgx93eEBZCgpxIQ+pFon1RaGecedWGjfVObXQKlPPpGvStae6JQUNFGqGkkpQRApdMG3pzHQvL9UtPQU0uPtIDHMYHtCQ3j8fWhW2515gIbSqUres+/v2we76+fl0P99x45tzAAAAAAAAwIHwcA4AAAAAAAA4EB7OAQAAAAAAAAfCwzkAAAAAAADgQHg4BwAAAAAAABwIaa0PoPkjTbpxsSSaUoDNrkloqqV0uYS0GK5U2F4p2XWfSbI2TW14ulkMJqu8p/GNGp7q5+MQJbnESDeupJWJgu0h9cYlYaVEqElILFtMTCqYaVu3D083c0mrkjQ1Y+e+29urpLKmt78yxhhSYpkTTy2X4ltIdl2mNNOU+OoSxNL2zNgpxSwm9Zl0s/j6XBJWen2VdLOY7Dk8NcsPUGx3hy5sz8+TwxPCw2GLaWoTl2RWSU8tpJhJIUEsbs8llvm+7r3OfYvtxr6mnJ3rltBeSapP/WNSvbuHV+okhZoh1UQ2/Tj0da+jnCRr3uz/v737DpOmKhM2fj+AJEFAkCQKCCIqGAiuAQOCnxGVT1gXFMVVMeEiZldRxISyqKgoqyuCi+6aWNHFwCcoKiKCohiIEkSCCqIEJb7P98epZuodqnq6Zrqne2bu33XN1dVVp06drq7uOfPMqfN0aFvn74vGRrRU0fCBb32v27LPN65tyVrZ8Pu+/Xu9IcNlh34LNPdRbm0r2/ChbMvK2qU/09a/aKqjtWzj2maT0p/p0ormPkrzNdR0vXTuzzRmDR68j9Kl3wIt/bUOfZRsqbex79L2d0eX31td+qhd/8RufCltjWs4920lG76fomv2+abfI0PIHN/+/du0rq1zN9j+0JKZdcz9ltvafv/PkSPnJEmSJEmSpDExOCdJkiRJkiSNicE5SZIkSZIkaUwMzkmSJEmSJEljYnBOkiRJkiRJGhOztS5Ad6wexJiytXbOYNOY3WwIx2vLIDbX43XI1trahpasWV0yvjbW3fo6OmT2bEuY05gxp2O2qg4Zb7LhhWdLNp9oyoTWco7vaMn+1FRHa/avpmxjLcdrzZrUUEdbJrSmOtrLzj2L66h0SXjUlm2sS9nGBMMdMq22lW9JMEw2ZD3r0rZSR4cMrI1lm+ttrKMtRWnXrGdzLNspI2aHfxdm22ek4XW3ZjFrORdNXyOdMq12yGLWVn5UGcs6Z2ttbMPgZUdpZP2ZIWR87ZapfkTHG0b2+ZbffY3ZaLv0y7r+gprr91PbZ6/Db8q273UaMvYta8tO2HA+234/depftGZgHWx/6NfPaerPDN62Nou9P9Opf9EhU31bHY3XUYds922ZVluzqDdlRW7td3ToozRla+2QmbmtfKffW22FO2TnbD9cU+NaijYcrzXTapf1XTKttr2QTn2iDhnsh9EXmaf+zB1tX4Rz5Mg5SZIkSZIkaUwMzkmSJEmSJEljYnBOkiRJkiRJGhODc5IkSZIkSdKYmBBiAbp9VWC1We7cYe7CuU623LmOLvUOI6lEp+QRTUkJBq+3rXyntnWdvLqLDhMoR8sBGyeObprYta3ulok1m+ptTxQyeNKFZa0JIZr2H7xs2/HaZh3t9DHrMnPpfM/aPucPe8ukyENoQ+tE3h2SSjSWbZ2MuMPEyl0mYe5SdggTKA/jPW29Dhs+JK2HW7Fp/5ayTd8XrRMlt0xS3HhdtNXRsLLjJMVdkkoMvP+w6uhiVF85Q7gM570/M6oEFENJHtHwRnVI5tBWvvXt7/D7sJPWF9iQwKA1wVXbbPcdJrBv7KN0KdtctLXP0PT+tWS9aewzdOm39GlHY9FB29DVBPdnOiW46nC8Lv2W1na01tHheCPro3RI5jDKBFdNuv4O6JLgqvFD0nK8hs96W+KHtuuwU0KIESW46mLe+zNzaO/tI/pacuScJEmSJEmSNCYG5yRJkiRJkqQxMTgnSZIkSZIkjYnBOUmSJEmSJGlMDM5JkiRJkiRJY2K21gXojpUhVpmh0DAyeDbonLxvZNnNBk+R0im7TpcMPR0zrzW2uEuGtFFmq2o8YMvx2jIfNtXR5Xy2/augKUtT67loycjW6f3rkE2tzRzfv9asacMw1+toGBk826rulGGpS+rpjusHPV7Xeufa5tYMWx3a1pph6651dLlU2l5a2ytu/Lpo+z5sakeXfy221tuSDbFL9tSmzHId3/85Zxab5GSIXeudlL5Lky6/R7q0oeXkz7XfMZS+T0sVo+o/dXmf2q7ZxizardnSB//iytYDDt73acxg3yGjddvhWq+hxnqbD9etv25/ZsaqR9W/aK1jwHpby3Y9Xodz19DvaOy3tLVjlH2tLto+ZysM3n9qzB7flJEeGl9fW8buaLngGld3ydbapkt/ZkTZU0eaxHmAuu+4fTSHduScJEmSJEmSNCYG5yRJkiRJkqQxMTgnSZIkSZIkjYnBOUmSJEmSJGlMTAixAC1bJbljlVnOgjiquU+HUG+n+TtHVLbT5Mcdy3Z7fU2zjA+h3rZJQBuLdpwhukObo2mG4LaJ6rtMdN2m07kf5ezFTcdrSEAxlHqHUUkHo5yYtfF4o5o5fkRl6TMB8qB1jyqxBSOeWLfJCg3XfYfJnVsnWx54Zfv6ptWt71zThM2tZTtMoj6M92PeP5Mjqne+v8sWYH+m06mfY2IL6Ji0avy/Utv7M8u6JDbokPVmGIlCWhsy8MoR9vk7JGUbyvGGUUkHI0uQM88vZL77M8M4XpfEFi3mmpSga4Krpspbv3Oakj8MoW/X1r3o9Cdip/dvkfRn5lDvsttH0yhHzkmSJEmSJEljYnBOkiRJkiRJGhODc5IkSZIkSdKYGJyTJEmSJEmSxsTgnCRJkiRJkjQmZmtdgJatBHG3cbdiQIslO2xjG4aQpWWUGb26GDzRap867rrHnDPUthUdZWayec8MON8HXCQmMHPTUOue4we+89fTiF53a9azLilRmzIaz+O5LHXMvYoudQ/lcPOdrXWxmOB+S5t578/M8XhD6c90aHJ7tuXmLY2ru5yiLv2ZUfYBRnYt22+ZlYWYQXsodQx+vUxyZs9OWVw7pJTv9H3Y8cuzsas1z+e40+EWUL9l2YiiaI6ckyRJkiRJksbE4JwkSZIkSZI0JgbnJEmSJEmSpDExOCdJkiRJkiSNicE5SZIkSZIkaUyWbLbWiFgZeC3wQmBz4I/A8cA7M/O6ln3WBN4K/COwMXA58Dng0My8eT7aDZArJbnSAkpnMgoLMFHUSLOtjsIwktG21dGYCW0yTtCENEMTZijZrUZlktvWYu7ZJVtWjyg1mdlTNTIT/Dtn0fw+7JqMtks26cbjze+JWzTvk4bKfsv86PT5ayjb6X2yPzMxRhWLWVQj5yLiXhHx0Yi4LSIyIjZrKbcacDLwPmAr4NfAOsABwBkRsV7DPusBZwJvAu4N/BK4D/AO4OSIWHX4r0iSJEmSJEmL2aIIzkXEahHxJuAiYH9mHhF4OLBTVf6+mbk9ZfTcBcD9gSMa9vks8ADgNGCjzNwReBBwLfBo4KAhvBRJkiRJkiQtIYsiOAdsDRwK3EIJzrWKiC2A/aqne2XmFQCZeQ3w/Gr9HhGxbm2fxwFPBf4G7J2Zf672uRh4eVXsxRGxZG+Rv2LrAAAgAElEQVQTliRJkiRJUneLJTh3DuX20i0z88gZyu4FrAj8PDPPqm/IzDOB3wArA3vXNr2gejwpM383rb6vAtcDG1ACeJIkSZIkSdJAFsVIr8y8AzhkwOI7V48ntmw/mXK76pOBj860T2beHhGnArtV+3x9wHbMWq5QfqQFa4InEp3oCXQl9TXRE6NPctsk9TfBfQP7LdLCNdH9ljYLsc1DNqpYzKIIznW0SfX4y5bt51WPGwFERAy4z269fUZguY/AHTfdNKLDaNHzy1TSKPjHoaRRse8iaRTsu2iWGuIxQ/lNtaSCc9MCbde0FLuueuwF2taj3ObaZZ8ubbphkGL1J1ceeljXw0iSJEmSJGm4Lo7oP445M9ecqZKldnPk2sDq1fK1LWV6gbYNImJFpoJ5g+wzqpFzkiRJkiRJWoSW1Mg5ln+9K7aU6a2/FVjWcZ+buzZokAhqRGwNnNu1bkmSJEmSJI3MDpl53szF+ltqwblrgduAuwHrtpRZp3q8KjMzIq6sbVuX5tFzd+4zlFbe1YXAA2vPe1lmNx/R8aRxuqR69PrWUuJ1r8XM61tLkde9FjOvby1Fvev+fsA9a+svHEblSyo4l5nLIuIq4L6UueSaTA+0/YEyXWRU+1wwwD5DVWWjvTMS27ufOTP/NIrjSePk9a2lyOtei5nXt5Yir3stZl7fWopq1/0fgT8Ou/6lNuccQG8k3ENatm9TPV4BkJm3UwJ0A+8jSZIkSZIkDWIpBue+WT0+rWX7LtXjCYPsExGrAI9t2EeSJEmSJEnqaykG546mJHrYNiL+ob4hIh4DbAH8BTi+tumT1eOuEbHZtPr2AFaj3Hb6oxG0V5IkSZIkSYvUkgvOZebvmRrhdlxEbAwQEesBx1brj87Mv9d2OwM4mxKEOy4i1qn2uR/w0arMxzMzR91+SZIkSZIkLR5LLjhXeSVwGbAlcHlEnAVcTBk1dzrwtnrhKuj2fOA64DHA1RHxE+DXlGQQxwNHzlvrJUmSJEmStCgsyeBcZl4N7AAcRcmysQ0lUcRBwM7TRs319vkNsB3wOeB6YFvgfGB/YM/MXDY/rZckSZIkSdJiEd6JufBExA0AmbnmuNsiDZvXt5Yir3stZl7fWoq87rWYeX1rKRr1dW9wTpIkSZIkSRqTJXlbqyRJkiRJkjQJDM5JkiRJkiRJY2JwTpIkSZIkSRoTg3OSJEmSJEnSmBickyRJkiRJksbE4JwkSZIkSZI0JgbnJEmSJEmSpDExOCdJkiRJkiSNicE5SZIkSZIkaUwMzkmSJEmSJEljYnBOkiRJkiRJGhODcwtIRKwcEW+OiHMj4uaI+F1EfDgi1hl326SZRMTWEfGRiPhhRPyp+jk5Il4cETHX8tJCEMVDIuItEXFaRPw5Ig4Zd7ukYZtNnyUi1oyIQyPi4mqfCyPi4IhYdT7bLs2G/RYtFfZltFBFxBoR8aaI+N+IuDQi/h4Rv4mIt0bESi37zFsMJjJz2HVqBCJiNeAkYCdgGfBzYCtgDeBC4NGZec34Wii1i4gHAb+unt4M/ArYBNiwWvcFYK+svpC6lpcWgojYHPh34EnVqtuAi4HDM/NTVZmDgXcMUN1lmbnZCJoptYqIewFvB14OrARsnpmXNpTr3GeJiPWAHwIPAG4FzgG2BVYBfgTskpk3j+SFSX0Mct0Ps98SEZ8H9qqeNn7GpGEZ9Hu9Vn7GvkxV7mDsz2iCRMSOwNeBDYC/AucB96Z8VwN8H9gtM6+v7TPrGExEbAa8B9gbIDNn/CeNI+cWjsMpF8VFwH0zc3tgc+AC4P7AEWNsmzSTc4FTgd2BtTJzR2Bj4FXV9ucCz55DeWmiRcT2wE8ondlzKJ3gDTNz63pnltLBPbnPzyVVuT/MU9MlImK1iHgTpQ+yP+UPuH5m02f5LCUwdxqwUfW9/yDgWuDRwEFDeCnSwDpe90Ppt0TEXkwF5qSRmcX3epe+DNif0eS5jPIPv/2BDTLzkZl5H+DxlL7G44DXT9unc38mItaOiPdTgn97d2mgI+cWgIjYAjgfWBHYMTPPqm3bkfIleSuwcWZeO55WSv1FRDT9xzgiTgF2Bo7JzBfNtrw0qar/Sp9N+e/cF4B9ZzMCqLot6ifADsA+mXncUBsqtYiIhwM/A/4EvBP4WLWpaQRR5z5LRDyOEtj4G/DAzPxdbZ89gC9R/oDbJDNvH8VrlKbrct1X5efUb4mITYBfUv4I3KHfsaS5msX1PZS+TFWX/RmNRUSs2xQviYhXUT4DNwJrZ+Yds43BRMTuwPHApcAngPeDI+cWk70oF8XP6xcFQGaeCfwGWJmOkVlpPvW5lePs6nHTuZSXJtgHKZ3ZbwF7z+HWvJ0pHdnLKR1jab6cQ7k9acvMPHKGsrPps7ygejypHpirfBW4nnIbylNn13xpVrpc93Pqt0TECsAxwNrAa7s1U5qVTtc3w+vLgP0ZjUmfgUzfrB7XADarlmcbg/kB8AbggcAXu7TP4NzCsHP1eGLL9pOrxyfPQ1ukYduyerxsROWlsYmILSm/tG8DDsjMZXOo7s3V44cy87Y5N04aUGbekZmH1Odh6WM2fZbWfaqRcqc27CONVMfrvp9B+i0HALsAX8zMH8zxeNKMulzfQ+7LgP0ZTZ4bastrV4+zisFk5jWZ+W+zCWDPeG+5JkJvksJftmw/r3rcaB7aIg1NRDwE2LV6euywy0sT4HmUf4SdAPw+Il5CmdtiNcrtJMdn5nl99gcgIrajzPHyV+A/Rtdcac469Vmq25sG2Wc37OdogRmk3xIR2wDvA64EXjlPTZO6GEpfBuzPaGI9uLZ8fvU47zEYg3MTblqntS0b63XVo51WTbzqmt6Y8ofWO4HVgQ8xNTJiTuWlCdP7o+wyyq1NW9W2PQc4KCL2z8xPz1DPm6rHT2TmDX1LSmMyyz7LepTbQrrsI02sLv2WiFgFOI7yGdjHuaM1oYbVlwH7M5pMT6weL8jMG8cVg/G21sm3NuWXOpQsIk16F8YGEbHi6JskzU5EvIaShvr3lAky7w7sA7yuZRLlTuWlCbRd9fgKyi/3nYA1KXOtfBtYFfhURDysrYLqdpI9KJPOfmSkrZXmZjZ9lk1q22bax+CcJtos+i2HAA8FPpCZp8xbQ6Vu5tyXAfszmkwRsS5lagGYyr46lhiMwbnJVx/d2Pam99bfSukQSJPqcuAUyvDgv1E6re9majLwuZaXJkZE3J2pX+y/BHbJzNMy88bM/CnwNODnQFBuaWrzesrv6//MzKtG2WZpjmbTZ+myz1wmIJfmw8D9lipL8RuAs4C3z2MbpYENsS8D9mc0md4F3AO4Aji6WjeWGIzBucl3LWXyTYB1W8qsUz1e5WgiTbLM/Epm7pKZDwHWB15NufXjmIi4S3ayruWlCbNWbfk90yeGrSZU/lD19FHVEPrlRMSGwL7V08NH0UhpiGbTZ7mytm3GfebWPGm0Bu23RMRawGcpAby9M/PWsTRYmtmc+zJgf0aTKSL2oowIBXhp7foeSwzG4NyEq77wep3R9VqK2WnVgpOZN2Xmx5ia/Pjd1bDioZSXJsCfmPpP2uUtZX5TPa4FbNCw/QBgFeDrmXnucJsnDdcs+yx/AHqdWvs5WjRm6Le8FdiU8sfff0XEWfWfWjVfr9btNo9Nl+qG0ZcB+zOaMBGxLVNJSd6Tmd/sbRtXDMbg3MLQ+6/yQ1q2b1M9XjEPbZGG7UuUX/qrMTWnxTDLS2ORmbcBv6ue3rul2O215b/XN1QjK3p/1B023NZJI9Opz5KZt1MCdAPvIy0wTf2W9avHtYHtG356tqme+89IjcVc+zJgf0aTJyLuC3yDcsv2N4GDG4rNewzG4NzC0IviPq1l+y7V4wnz0BZp2G5i6pf6IBmku5aXxqmXZn3Hlu33rx6vyMy/Ttv2csocGGcAPxxB26RRmE2fpXWfKpvlYxv2kRaKu/RbMnPfzIy2n9q+m1frjpnvRks1c+nLgP0ZTZCIWI+SyGQT4HRgz+ofhdPNewzG4NzCcDTlP27bRsQ/1DdExGOALYC/AMePoW3SXD0DWLlaPnME5aVx+mj1+NKIWK1h+4uqx2/UV0bEqsBrqqcfcD5RLSCz6bN8snrcNSI2m1bfHpQRR+cBPxpBe6VRs9+ihW5WfRmwP6PJEhFrACcCWwPnAs/IzJtais97DMbg3AKQmb9nKiJ7XERsDHdGfY+t1h+dmXcZRixNgojYPSJeEhH3qK2LiHg68O/VqiMz85rZlJcm2DeBn1DmYPlaRKwNd17PbwaeSpkQ/OBp+70A2BC4CEcLaQGZZZ/lDOBsShDuuIhYp9rnfkz9Ufhx/6jTpLLfokVutn0ZsD+jCVGNxD8eeARwIbBrZv65rfw4YjBhP2dhqDLc/JgyeewySid2K2BNynDMXQzOaVJFxA+AnaqnlwFXAw+gzLUC5T9te/Su4a7lpUkWEVsB3wM2Au4AzqEMpb8XcD3w/Mz8eq38ipRRQlsCr8jMo+a7zVI/EdHrPG6emZc2bO/cZ4mIB1Fud1oHuBX4BbAtsCqlM71nNUGzNBb9rvth91tm+oxJwzbA93qnvky1j/0ZTYyI+E/g+dXTXzE13+10J2XmB6p95hSDqe4GuARg2pQFzeUNzi0cVZT2XcCzKZ3XS4HjgMMy85YxNk3qq5p0cz9Kx3ULysTGvwXOp/zHYfotfZ3KS5OuytD3FmB3Smf2r5RO7tsy84JpZfcEvkjJkLapQWhNmkECB7Pps1Sd2HcDTwbWoHznfwr4hIE5jdsMwbmh9lsMzmm+Dfi9PnBfpipvf0YTIyIupQTZZnJsZu5b22/WMRiDc5IkSZIkSdIC4ZxzkiRJkiRJ0pgYnJMkSZIkSZLGxOCcJEmSJEmSNCYG5yRJkiRJkqQxMTgnSZIkSZIkjYnBOUmSJEmSJGlMDM5JkiRJkiRJY2JwTpIkSZIkSRoTg3OSJEmSJEnSmBickyRJkiRJksbE4JwkSZIkSZI0JgbnJEmSJEmSpDExOCdJkiRJkiSNicE5SZIkSZIkaUwMzkmSJEmSJEljYnBOkiRJkiRJGpOVxt0ASZIkaRwi4hnADtXTYzLz0jE2ZywiYg3g9dXTSzPzmDE2R5KkJSkyc9xtkCRJWlQi4lJgU+CWzFx1zM1Ri4g4CnhZ9XTnzPzeGJszFhGxIXBV9fTUzHzCGJsjSdKS5G2tkiRpXkTEWhHx3Ij4dET8IiKujIhbI+KaiPhVRJwcEe+NiJ0a9n1CRGT1c9QQ2rJDRHwgIn4WEX+o2nF5RHw/It4YEfed6zGWkoh4UkRcXb2vnrsZeL4kSVKdt7VKkqSRqm6bOxB4LbB2Q5F1q58HA08EdgUeMaK23Af4ILBHw+ZNqp/HAodExL8B78rMW0bRlkXm3cAG1c8rgTePtzkTz/MlSZLuZHBOkiSNTERsCXyVEngDuAE4HjgZuAS4jhKw2xDYCngW8OURtWVb4NvARtWqXwBHV4/XAfcBngq8EFgDeCvwiIj4v5l54yjatIj8uWX5ThGxNbA78JPMPHleWjW5ZjxfkiRp6TA4J0mSRqIKzJ0B3LNadQRwSGb2C0a8LyJiBG3ZCDgFWA9YBrwR+FBmLqsVOwc4MSIOowQIdwCeBHw2Ip6TTtTbz+GU83UVcNz0jRGxN/C56umL5rFdk6rv+ZIkSUuLwTlJkjR01a2sJzAVmNs3M48dZN9hB8GqYN+xlMAcwIGZ+ZE+x78sIp4MnE4Zzbc78DwMorTKzO8A9+pTZOX5astCMMD5kiRJS4gJISRJ0igcADyoWn7/oIG5EXkcZQQclNFzH51ph2p03+tqqw6KCPtNkiRJGjo7mZIkaagiYk3g9dXTa4B3jbE5UJJR9BzcYWTeN4Dzq+WtKIkqhioiNoqId0XEWRFxXUTcHBGXRsRxEfGkPvutX8tee0C1bsuI+GhEnBsRf4uI30bEZyNiswHb8tCIODoifle145KI+HJEPLra/rHqeLc37HtUrT1PqK3/XkQk8Jla8c/UymZE7Fsrf3Bt/T/1aeu+tXKtyRQi4r4R8eGIuKB6TVdXWYH36XL7dERsVmX3/XVE3BgRf6kyrR4cEevNXMNd6ms8X9W2R9S2Pata97CIOCYiLoyIv0fEeRFxZETMOPouIlaIiL0i4qSI+FN1bZwfEUdExOYd2tyr58QomZZvqa6V/6lGmk4v/4za6zgzIlZsqfc5tXJfG8Vt7ZIkTTpva5UkScO2C1NZWT+dmTeNqyERsTJTo+Z+B/xw0H0zc1lEHAF8vFq1C3DSENu2H/BhYLVpmzatfp4XEV8D9snM66eVubW2vFZEPJcSAKvXdb/qZ8+IeFFm/ndLOwJ4e/VT/8ftZtXPcyLiNbX1N8/86sYvIvYCPs3y56SXIfWJlNuVrx2gnhdTRltOf58eUv3sFxG7ZeZPh9Fu7vreHgj8G8u/Nw+ofp4XEc/KzFNb2r42JQHLztM2bVX9vAjYb6YGVUHA/wEeM23TfaqfZ0fEUcCrevM4Zub/RsRngRdQ5tfbnzLvZL3eNWvr/gzs59yOkqSlyJFzkiRp2OqjaIYWzJql7YHVq+WfzuIP/zNqy9MDE7NWBVz+nRLwuYkSGHs0sB3wEqZG7D0T+HYVZKy7rba8K2U+vL8CrwZ2BJ5DCaYArAoc22eU1EHAwZR+4U3AIZRA5E7A2ygZdj/E1Ou/pcNLfR5wf6A+uu3N1brez/Ed6htIRDwT+DxTAbX/APakBIleDFxACc7tPUM9L6j2XQ04ixLI2hF4NnBkVWwjSiKRtYbU/Pp7+0Lgg8DFwD9TrufnAb1g3FrAlyPiHg1tvxvwTaYCcxdSRrQ+ivIZ/Thwd0oAs1VErE65HfwxVdsOA54OPBJ4FXBuVfTlwBum7X4gcHW1/J6IuO+07YcA966WX5WZVyNJ0hLkyDlJkjRsW9WWfzW2VhRb1pYvncX+l9SWN5pbU4qIeBhlJBTAdcDTMvPHtSJnR8QXKIGVnShBkDdTAhk99dFVjwUuBx6ZmVdW686KiK8C/0kJQK0MvJYSvKu3ZTvgHdXTq4GnZOYvakVOq+o5HXhYtW7gkXOZeUV1nD/UVv8hMy8atI6uqkDV0dXTO4AXZObna0V+GhGfA74NPL5PPfcGPlY9/SJlBGPvvJ8FnBARp1MCoxtQgn4fHMJLqL+3TwTOBh6fmTdU634WEV8BvgU8gZLo5KWUDLB1r6dcOwCnAc+clin5pIg4Efj6DO15F7ANJSj7+MysB6zPqM7lKZTA8usi4ojMvBnK3I0R8TJKcpi7U87nM+HOz8G/VPV8BfjCDO2QJGnRcuScJEkatvVry39uKxQR+0+be6z+c8yQ2lKfD+yS1lItMvM64C8Ndc3FQUz1wQ6cFpjrHfdGyi2HvZF+B1YjmHqmz/t2YC0w16tjGcvP9/eUhra8udaWf50WmOvV82uWDwzeZc65CfNiYN1q+dhpgTkAMvMW4PmU4F2b/YE1KQHUV9QCc3WfB35WLb96SPOlTT+/L60F5oA7239obdVy721ErMrUXIu9AOVdPouZ+Q3gk20NiYh7Aq+onr5vWmCuV8dfgXdWT+8FPHfa9q8xlel4t2ouuqAE6lYA/kQ5v97OKklasgzOSZKkYbtbbXncf3DXA1qd50qrggi9W0qnz/vWWRVge3r19Grgv9rKVqPLvlE9XZsyQq63rX5eb6WMTGqq4zzgiurpFvVJ+SNiFeAZ1dPfMxVAaVIPcI37PZ3J7rXlD7QVyszfA9/tU08vIcWJTYGtqo6kjEqDMj/fOoM3s1X9/F7cZy67etvvP23boymBMoAvZebFfY7X733fjalbg/uVq8/luH3D9gOA3ujJIyi3bvduk355Zv6pT92SJC16BuckSdKw1f/Q7hesOI7l5x6bPmn9MNSDKvebxf7rMxXgu2zuzWE7YJVq+Qcto7HqTq4t79BS5vLM7DearTeiLlj+/Xg4U4GXn2bmbbS7iskfMUdErMDUrZx/zszz+5WnJAlpqmdDSrAN4Pl9Rngmy98qPHD20wG1jvasrp3eZ236qM5H1ZZ/NMMxGs9BQz0X9TkH9cQadzkHVXDzZdXT+zE1Wu/zmTn0OQclSVpoDM5JkqRhq99e+aC2Qpn5l8y8qPfD7OaEm0k9ULhla6l29YDeMIJzG9aWB7nNtl7mXi1lZgqa/a22XJ9vuN6WfiOreiPE+t0COinWZmrk5iDnty0gufEsjz+M21rrBn1vp88jvUFteabz0C8oO5vz0HgOMvMEytx9PTcybQ5ESZKWKhNCSJKkYTsV+Mdq+YnA98fYlrNqy9tGRHSc2+ohteULh9Cern2vepbWLllSu9b9t9ZS4zObQNeqteWZRiX2s1pt+bMsP3dfP/1Goc2nUZyHHSgZgWdyY9PKiFgNeGht1RqUhCdfm3XrJElaJAzOSZKkYft2bfmfI+I9M9wyOTKZeWlEXEwZAfcAym2lbXN4NXlerypK5tO5qo/kmz5PWJN65turh3D8ugtqy5v1KxgRK7H8XIKjUg8krdWnXNvdH/XbmDcb4Hirtqyv36a5/iizy45Ivf2bzVC27RxMr2elOZ6HQyifwVspAbx7Ap+IiFOrpBKSJC1Z3tYqSZKGKjN/C/xv9fQ+wL+MsTkAn64tv7O11DQRsTNTSRi+Vb2uuTobWFYt7xgRMwW8dqotzzR3WFf14Nz2M2Qa3ZS59RvroxX7HecvteUNW0u1zO2WmTcz9bo2ioh7z9CurVrWXwzcVC1vM6QsrPPpnNpy21yFPW3nAKCevXfb2TYmIh4JvLZ6ejjwxmp5Y/ok7ZAkaakwOCdJkkbhYKbmKDs0Ip40xrYcCVxXLT89Il4y0w4RsR5Tk9YncNgwGlJNjH9q9XQT4EV92vAo4MnV08uBM4fRhlpb/sZUEGdr4Gl9iu81x8NdU1u+R59y9QQOuzYVqAJlT66vmlbkm7XlA9sOFBEPYvmEB3fWUyVb6I0A3QT45/YmT6STmZqvbp+IWL9P2RfXlqefyxNry2+qbk3tpNrnGMrfHVcC7wU+A5xRFdkvIp7QtV5JkhYTg3OSJGnoMvOnwOurpysB34iI10bEKn12G1Vb/grsW1t1VES8om00VERsQQnM9BJIvDYzvzvEJtXnL3tPRPxDQxs2Z/kRf+/IzGXTyw3BEfXl6rVPb8t2TI10mq0rasuP6FPux0zNf/eYiNipoczrgR1rz1eetv1jTAWGXxkRT51eQUSsSckWXDd9upd6QPawiNi+T7snSmb+kanXtyrwmYi4+/RyEfF8puaHhGnnIDPPYSrYuSXlGun6Ge7dzgrwxsy8sbqWX8nUKNJPRcTqjXtLkrQEOOecJEkalSOA1YF3U/ochwMHRMQXKUkirgRuBtaj3P76nAHr3aoKKvTz+8z8Xu9JZn4tIg4APgysCHwc2DcijqaMHrupasNuwD5VuwEOy8wPD9iugWTmdyPiMOANlNf+vYj4MPBdyvnYiTLia71ql89RRh6NwjHA3sAuwBbA6RFxBHAaZW6w/wO8jjLycHXKuZtNgoHzKHPmbQg8NyLOpLze9YAzM/N6gMy8KSI+Cbym2u8bEfF2yi2961CCrP9EuX4eV5VZbs60zLyo2uc9lIQGX4+IjwGnUIKE2wP/SrlV9zigdy2tMq2eH0fEe6uy6wA/qdr2HUpykDWAjar6ngVsl5nDTtoxF2+kjDDciDIq8gfV9X4G5fXsBbyQcv2vAGzDtHNQ2Y8yT+P6wEuBx1fXyK8otyFvSDmXuwGnZeb7eztWoz97t7OeAny+ty0zf1a9L/9CCfwdzNyDwJIkLUjRLWGZJElSNxGxC2UU0sM77PbpzLzz9tPqtrcuo9dOyMxnN7Tl6ZSRVZvNsP+PKSPmTu9wzPpxLqUELG7JzLtMuF+N2nsHcBD972T4IPCW6jbL6XX0OnHnZ+bWfdryPeDx1dONMvPqadvvAZwAPKGliquBZ1ACZCsDZ2VmfeQaEXEU8LLq6c71wGitzKuBjzTU/8DMPK9W7u7At1h+vr2671ACdL1bZT+cmcvdvlqd30PpH+x5OyXI12vrnpn55YZ63kp5r2b6p/bjMvMHM5Tp1dt6viJiM+CS6um3M/Mpfeq5lP7X2QMpI982baniPEpQ7VOU9//XmblNQz33B46nBPD6OSUzd6n2WY0yx+IDgNuAbTOzftty79o7lzL33DLgEdWoW0mSlhRHzkmSpJHKzJOrWwJ3pozgeSxllNo6lPncrqIEgK6i/DH//+iWUbVLW06MiO8Az6WMdtqeMvKnPmIogd8Bj4iIO4CfNwXH5tiOBA6OiC8ArwKeRDknN1JGJP0c+I/M/M0wj9vSluurAOo+lDnwHkoJwl0G/A8loHYrU7ePXttUzwA+RhmBtz8lyHMdJXnDDdPac1MVjH0l8ALKfHg3UoI4nwK+lJm3R8RNwN0po+Omv6akzJH2NcrIrMdRsoNeCZwOfKQaGVe/VfUuwa2qnndHxH8Dr6DMg7dpdcyrqvouoAQ3z+p+SkYrM8+NiAcDBwB7UpI/3ApcBPwX8KnMvCEibqx2aczcmpkXRsTDKaPt9qAkmbgX5b27kjIi8fvAV2q71W9nff/0wFxV7/XViNYvUYLUR0fEDuPK7ixJ0rg4ck6SJC1p1eioB1Bu4dydEsipj2a7hRIsPB14W5URdEmJiEdTbnUFODQz3zLO9kiSJC0mBuckSZJqqkytuwHPpgTseqOJfpOZDx5bw8YoIg5nau6wp2bmt8bZHkmSpMXE4JwkSVKLav6zp1ACdWdn5gfH3KShq+b96mXQbNq+KyV77QqUWzgfOKLMsZIkSUuSwTlJkqQlLCKOBJ5OyVz6K+BSyrxk96Pc5rt3VfQ24ImZ+cMxNFOSJGnRMjgnSZK0RFXz7V0IbDFD0auBF2bmSaNvlSRJ0tJicE6SJGkJix5juQYAAABQSURBVIjVgX8EngNsRskaG8BvKVk9vwN8PjNvbKtDkiRJs2dwTpIkSZIkSRqTFcbdAEmSJEmSJGmpMjgnSZIkSZIkjYnBOUmSJEmSJGlM/j/MhWYNJHJTiAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "image/png": { "height": 332, "width": 627 }, "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "grid = coeffs.expand(grid='GLQ')\n", "fig, ax = grid.plot(show=False) # show=False is used to avoid a warning when plotting in inline mode" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Next, we would like to calculate the variance of this single spherical harmonic. Since each spherical harmonic has a zero mean, the variance is equal to the integral of the function squared (i.e., its norm $N_{lm}$) divided by the surface area of the sphere ($4\\pi$): \n", "\n", "$$N_{lm} = \\int_\\Omega Y^2_{lm}(\\mathbf{\\theta, \\phi})~d\\Omega$$\n", "$$Var(Y_{lm}) = \\frac{N_{lm}}{4 \\pi}$$\n", "\n", "When the spherical harmonics are $4\\pi$ normalized, $N_{lm}$ is equal to $4\\pi$ for all values of `l` and `m`. Thus, by definition, the variance of each harmonic is 1 for $4\\pi$-nomalized harmonics." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can verify the mathematical value of $N_{lm}$ by doing the integration manually. For this, we will perform a Gauss-Legendre Quadrature, making use of the latitudinal weighting function that is stored in the `SHGrid` class instance. Note that it is necessary to ignore the redundant longitudinal band at 360 E." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "N = 12.56637061435917\n", "Variance of Ylm = 0.9999999999999999\n" ] } ], "source": [ "N = ((grid.data[:, :grid.nlon-grid.extend]**2) * grid.weights[np.newaxis,:].T).sum() * (2. * np.pi / (grid.nlon-grid.extend))\n", "\n", "print('N = ', N)\n", "print('Variance of Ylm = ', N / (4. * np.pi))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Alternatively, we could have done the integration with a 'DH' grid instead. In this case, it is necessary to ignore the redundant longitudinal band at 360 E and the latitudinal band at 90 S." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "N = 12.566370614359174\n", "Variance of Ylm = 1.0000000000000002\n" ] } ], "source": [ "grid_dh = coeffs.expand(grid='DH')\n", "weights = pyshtools.utils.DHaj(grid_dh.n)\n", "\n", "N = ((grid_dh.data[:grid_dh.nlat-grid_dh.extend, :grid_dh.nlon-grid_dh.extend]**2) * weights[np.newaxis,:].T).sum() * 2. * np.sqrt(2.) * np.pi / (grid_dh.nlon-grid_dh.extend)\n", "\n", "print('N = ', N)\n", "print('Variance of Ylm = ', N / (4. * np.pi))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Parseval's theorem\n", "\n", "We have seen in the previous section that a single $4\\pi$-normalized spherical harmonic has unit variance. In spectral analysis, the word *power* is often used to mean the value of the function squared divided by the area it spans, and if the function has zero mean, power is equivalent to variance. Since the spherical harmonics are orthogonal functions on the sphere, there exists a simple relationship between the power of the function and its spherical harmonic coefficients:\n", "\n", "$$\\frac{1}{4 \\pi} \\int_{\\Omega} f^2(\\mathbf{\\theta, \\phi})~d\\Omega = \\sum_{lm} C_{lm}^2 \\frac{N_{lm}}{4 \\pi}$$\n", "\n", "This is Parseval's theorem for data on the sphere. For $4\\pi$ normalized harmonics, the last fraction on the right hand side is unity, and the total variance (power) of the function is the sum of the coefficients squared. Knowning this, we can confirm the result of the previous section by showing that the total power of the `l=5`, `m=2` harmonic is unity:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Total power is 1.0\n" ] } ], "source": [ "power = coeffs.spectrum()\n", "print('Total power is ', power.sum())" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.5" } }, "nbformat": 4, "nbformat_minor": 1 }