{ "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": [ "from __future__ import print_function # only necessary if using Python 2.x\n", "%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 using a retina display, uncomment 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": "iVBORw0KGgoAAAANSUhEUgAABOkAAAKbCAYAAACgm96CAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAjEAAAIxABZy/d1gAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3XmUrNtd3vdnV3X1OedeXY1IQkIMAmRsBwESBNssZMBDDMsBDGLUwgYWYJbtxJiAYEGMFkmYDA4mGCI7TLFFADPYyAwmBBaRmcNkEjAzAgwySAyS7nBOd1fVzh9VR7fvqefpU7+u9+3q0+f7Weuuvnef3fvd9dY7/Oq9dfbTeu/Cva219rAk9d4f2vdcgKFxfOOq4xjH/YjjHlcZxzeuOo5x3I8u6rifjDk4AAAAAAAAgLvjIR0AAAAAAACwZzykAwAAAAAAAPbsYN8TwCCetO8JjOWqr3fA67u38frubby+exuv7553ZWsX6eq/f7y+e9tVfn1X+bVJvL57Ha/vnnchtQvfpDun1tpLWmuvbq29vrV2s7X2C621z22tXd/33AAAAAAAAHBv4SHdObTWXi7pNZI+RNJS0q9L+i8kfaGkH2utPWWP0wMAAAAAAMA9hod0Ra21D5X0pZK6pJdJepve+wsl/XlJj0p6kaSv3d8MAQAAAAAAcK9hTbqC1tpM0les//Nze+/ffPvPeu8/2Vr7BEnfJukjWmsv7r3/7EjzmEp6gWl/5hjb27MmXdnXJvH6LmqMsVzmuQ2B13dvz+OyvL6x8PqugCv8+q76+8frG/f3x3bZ57eLq/zapMvz+qhdzofXdwW01p4l6emnmn6t974YbPze+1BjXXmttb8u6bsl3ZL0rN77w3f8+UTS6yQ9W9L/3Hv/zJHm8acl/dIYYwMAAAAAAGArf6b3/stDDcZfd6354PXP77nzAZ0k9d6XWj3Ek6T3ubBZAQAAAAAA4J7GQ7qaF69//vQZff6/9c+3GXkuAAAAAAAAuCJYk67mHdY//+CMPn+0/vnc1tpk/e26rbXWNr6hJ+lJZ/3Ocz/n5Zo++GBlM8Dj+BvvAMbQ9j0BAFcWtQuAsVC/4C4Wjz6q133xl51u+qXWnnDgPHLn7/TeH9p2fB7S1Tx5/XObh3QHkp55l76DmD74oKZPOvM5HpBR6AIYA0UugLFQuwAYC/UL9oyHdDVzSdckTc/oc3qfPlrdwDZPWNdpKa+vjg0AAAAAAIDBPKv3/oahBmNNuprfXf989hl9nrH++ebe+8bXHAEAAAAAAIA78U26mt+V9C6SnnVGn7da/3zd+NNZ6U3q9/nj1nYv/rWHe3HOI7nU799lnhv26xL/dYh+ied24S7xvuB9gnOp74kS98VTLvV7dZnnhv25xPcd7ol3uMf2x/30/o39Wu/zRztl/2n9873P6POi9c/fGXkuAAAAAAAAuCJ4SFfz79c/P7C1trF2XGvtQNIHrf/z+y9sVgAAAAAAALin8ZCu5jslPaxVeMTfNX/+cZKeJulI0rdc4LwAAAAAAABwD+MhXUHv/U8kfcH6P7+otfaRrbWJJLXW3lfSK9d/9uW999/bxxwBAAAAAABw7yE4ou4fa7Xu3MdI+lZJv99a+0NJ77r+8++Q9IoLndFBX/1zLxhpmn2IcQcZo7CKZOhaWoT4gt/2nRdIrv7+jqtyDrKg80WfWvfIqXwpjbWI60jjlg/vwgFtuxa3d+ELAO+4vTH3p//93X59L2PfR4s6D2rM67IZe5DNjVnTmOZ78X574fVW4SI1aiDFWGNTv5zPPXavlQaqD8xBHo/7yzLnisL2xtqf9TH2/Pv7Gvs8Rn72wjfpinrvS0kvk/TJkn5a0pMlPU/Sj0r6BEkf2Xuf722CAAAAAAAAuOfwTbpz6L13SV+3/gcAAAAAAADYCd+kAwAAAAAAAPaMh3QAAAAAAADAnvGQDgAAAAAAANgz1qS7AvpsqT5bnvOXh53L4+NecATLiAldzb2WOEYhpi20210X9qdtHSLcpzhn33dzdsOM65tLR9xIx8ulT7e76OS1q5KkNVbfquYiFX1Xfx1J4/o/qLyUsfZnHNfNuZzuumVb0NMJf1mOF7u9i47PHGncC0/dvsQ1Tehva5fq9vr2E4ldK7VAmvJY9+HCnEvbG6DWGa2mqc7johNwK65KTXPRyaUD9C2dqnF7O9Y0Ui31vjJuRWXO1XvwjnVKnttIF6OrVNPc5bWc+9nLlvgmHQAAAAAAALBnPKQDAAAAAAAA9oyHdAAAAAAAAMCe8ZAOAAAAAAAA2DOCI66AdrhQu7Y4u9OIix4X1hWuzaOy8O5Y46r6+rYfNy7qbNeh9IN0s6pl3BNxHtv33TkMorrOutkXo4ZPFI6tYRan3n6I0dZGHWLckS4v5cvWjovb1re3ufPy4samqbxA8pZtSdrepLAE9IiLAvdJYXuFfVEKeCjtz+IFrbKI9BALNe940UjH5yAqcxupfind2+Mg6aCrjDHAuC6zasR7VClQy0wkza0tx7qZpPZwf3d1WOVULe/P3cbIddGOARhnzWPHvqNm2Ow69pj3uUGu7TtuL+x833f7bUn+vhFfs/mqULw2xH2x+Qd5H28OHqufwr4fJOAqfG3K1i8DHBel2uOi66Lq2K7raCEvj8+hHd3l2cuO+CYdAAAAAAAAsGc8pAMAAAAAAAD2jId0AAAAAAAAwJ7xkA4AAAAAAADYMx7SAQAAAAAAAHtGuusVcHBtoem1+fl+uRAz1Ct9B9ne9oP3lAhmtpfHTWmrboOFZNbQt4eJNBsd5DdXmltFJdHUptHW0spsklron9Lfdp5b2N4gaWzBrmOU09HGTFNzdjwUSymn1TEK4+Z5mOtL+l9flTRSl3Ka5pG2V/lfcOn8s69v90Sw0vtaSTyr/m9H91oGeH0tjrF9Slsr9fXtbnul02GI+MWxIhwHuM+NWdPY23tMkK8kpYbtuftiKnbsvTLVNGF7bpB4vzVjhL7pLdk1QTWmxg5wf7f1S7H2sPMrvb7t+8Z5XHDtMcil4RLXNKPWHruOkW4aOyazxv7pnuiO20oNFbbX0tzcIGPe59w8inVfKRXebq/QN/Qv1RiV13HW2Ntur/Drq+0Nd6Fr10h3BQAAAAAAAK40HtIBAAAAAAAAe8ZDOgAAAAAAAGDPeEgHAAAAAAAA7BkP6QAAAAAAAIA9I931Crh+7UQH10/O9buVTJuUhFZKSHOBYEv/rDgGk7nE1pRAZtNdUzpa2J4ZO6bJuvaUulVJNyukpsV0tCECbcw8clqZaQtBODF5zY2R9oVN1g19KyltQ/S9zGmyIxoimaw0bilBdcdx0xiF//WV+rZJuEZVtmcStvo0bc+328Ml7qQBLjAuWa6SFJb6xnYzbGGM1De2277+guYSzyZp3LA/XXsM9bN9t39PB8gUH8RYNU2lHknty1inmLZUF/lp7FynpPrHJsiHMdJR4JLsbeLr6g98e+VmUkjLLdU6hfqlVKek7aW+hXErqa9D1EXl/VwYY+dxBzBe2up4Y9ga4cLrotDZ9g3p9qF+sWOk88H1DcMmtX1fSKGvpLCOWae4uqhSY5TrlM22SaEuSipzjmPc5c/nR+d79rItvkkHAAAAAAAA7BkP6QAAAAAAAIA94yEdAAAAAAAAsGc8pAMAAAAAAAD2jOCIK+DJ1480u5FW1DzbshL6UBgjLaa8MIshx3WJw8LJi8L23ELNi4UfNy7q7OYc+5q5LcI+ju2mLS26asME0jtVGSOM4BYbjospu75+DnEM0z5MUIVvHy2oIhzkdn8WFnqOhli8eQA7L7JcWTQ59E/BCt2uVhv6FtpLfcNlO4dBbDbF0IepWbg+HBjxsCj8bzx3bMVxK4ssp/d6avq6Np2xkPF082RL+9MtZJzGTQEPU7O9uECy+/3Qd1pYqHlSWEw59rWtXhpjV5XaRfLHYhrDhj4Ug7MWrvYohEHE+qdQv6TgCDd2DM4K27P9w/3W1TppcfgYKGH7bt31jAI2jG1ey6RSvxSDHHatX6rBEbY9BpSZxeGLdUop+OpeDI7YsdbJdUpod58Jhgi4qtQ0hfolvg43hqld1r23b05zGyCAzR50hf0ZD9pQv7j2Vqh1YuBUqCdcf1e7pL4pmGEa2l0NlD7yuhoobe8i6peTo6PCaHV8kw4AAAAAAADYMx7SAQAAAAAAAHvGQzoAAAAAAABgz3hIBwAAAAAAAOwZD+kAAAAAAACAPSPd9Qp4+vVHde3G44knyxjbs51lyD5J6WauPaWYubHnKcWskO56svBRPi5hbR77pjS1zXSXRUj5sqlpLkVS0rLyjDykfLm3pIXtRW6QSoJqSgSbb45bTVuduHTXeRjD9S0kt622V0gxGylhzSWpxTGKqVSDpOHZgQvDVpJZi0lqLvEsp61uvsBKMqvkk8ls6lrqGxLIlimZzG7P93Um6doeI5PdJHxXG5pWTcJzryUkk/nEs5CUGpLQbNpqGMP1TYlnKcXswI0REtYOTLtrO2sMl0w2DRcj1zclm9m+xYvRJF5YN1Vqmkr9UqppwhzSGK6eSLWOa0/1yDykrbq6KCXBuvZl6pvOVZfYGuOxXWMYOKXMbh9m6VPTCzWN5JNcY/1iapJ67bF9X1/rhLphpITY3Hf7G0SsRyp9k4uuXwrj+vTTUDcMUL9UEluXJlm1mkJv65TwtMEeRvG4CDv0YPtUeDtEfFO3r4t6vMBsn0Kf2icHJtE09TW1wLRQ/0i+njiY+guJS5afpb4p9bVQ67g6o1L/pPZq/XLb0fWb5/q9bfFNOgAAAAAAAGDPeEgHAAAAAAAA7BkP6QAAAAAAAIA94yEdAAAAAAAAsGcER1wBz7j+qG5cD6vD3oVbkDkuphxW4nQLFs/Dyudzs724mHJYrfTYBD+chO25QImTpd9XKVDieL7ZvpiGOc9NyITCqqthhWQbKJEWNt1+7e2SuDauW0A4rCztFix2iyNLOQxi4hZkTgsvm31fWXhZqi3eXAmZaKZv6p8XuDZ/UFwM+5xro55fZUFmt9hwcTFlt/PcQshpjLwQcrgm2jAIv5OX5m6bxk135m7e19KCzEEMmzH7c8dcovW4qX37RZZdSMTELSAtaRIWMj4wCzKnxZTdYsguCEKSDsP23BiH4aLoFkNOfQ9iu1mQOQZHbLbPSiETI92MilLIhKtrTmJfU6eEvjkMYvPicBzqFBcSkfqmkCxXF80nfm4nE1PTpICWcENats0xluGq0+3F1m9O1eAru0EzbLw3b1+/xDplx4ArKdU64XpWqYtiDVSoX0y7q38k1ULHUkhW5ePMJalp7OWhEnAVzrNKTTJEGIQLJIjBESn46sBsL90eXF1UrWnMNSPtT3PZGuYQqoRhFQIiJGli+k8PCjVNCFY4TGOYOmUWxnA1zSxcdFJddGBO+BQc4eqXVP+koAoXEpFDJs6ua27eunXmn++Kb9IBAAAAAAAAe8ZDOgAAAAAAAGDPeEgHAAAAAAAA7BkP6QAAAAAAAIA94yEdAAAAAAAAsGeku14Bb334Zj147ejMPimZ1VmEiKCTEOXjU9NCApmJOEzJrEeh3Y3hks1S36OFP+yP5r7dJeO4xNeV7U+pHlJfXeJVd8mekk9CK4aj2VCbmPK1OXglEayS4ipJkxPXN6Q2mTHiuGEMn+46QN+wPZuWG95rm8aWYqliUN/2CbGDMMdiSmy1fWMyayWxNfQ1p2pMgo0pZtuP4RLrekgjTefUcrY5dsqesrs5Jf2lBDk3vcrxkhLrYoSxaQupky6FrpLiKkkzk26WEs8ODzYvJimt7NrUX3hc/+tTc5GTT3K9Fi5oKd3M9Z+F6ESXblbpO01JsCNdYHLafCXddfea5shdBOTTXVNfl+R6azHzfUOt4+qa1Hdq7kfHLvbwDO5I7GEf2folJXenVEbbPcVnbtmms1JfC30rtcdJIbG1UC+ltNU0D1dPpDHcPGKdEutBU78UEmJTEuxo9Us6tCr1SxrDJrOGzimx1SWoxhR6l9i6fWK929aq3c/NvddpDDtuaI/Bz7aODtcXdxym+idszh4bMd3VbC8m1vt2l+Saahpbp6SaJtUvZoxY05gLTBrX9ZWka6Z+iQmxO9Y0qX+qX+y4p46MRw/PfvayK75JBwAAAAAAAOwZ36S7B7TWHr6zaS8TAQAAAAAAwG2vbe3xr/D13h/aZTC+SQcAAAAAAADsGd+kuwfc+SS2tfZMSa/f03QAAAAAAAAgPb/3/oahBuObdAAAAAAAAMCe8U26K+DZh2/SQ4eH5/pdl2K2CM9uU7qZaz8JsT9HJlIxp7v6MW4uN19rSmy9aRLSHmt+X8V0urY5Rouxmo6f2zKlD5mUpx4SDm16YoxA2l4MX3SbC4mRLpwnBPbYFNdVu0kgGyAJdhoT1kwSUyEhts3DMZRS01zSVEyCdemu2yfBRik1bQjmWMyJZ+a6k/rGxNbNMWKKmUkbm6S+JlVVkiZzc/30oYw29S5dA3KsmPuD7dPm0uZKCc2+a+0wSv970CWhhWvfxLRPQzraQUxsNQmqJtlMkm7MNi8wKfHsgYNj2+6SXF2ymSTdmG6OkdJdU/t1c1FMSWiuPaembR4Y05DJlxLWduVqFynXL0uT+lqqaULfW0t/wrv6JdY0i82aJKX+ptTXx+abYxxM/PYqFUIMEDf7P6XQu/ol1jQp5dIO7Jtt+GLhGpfac1q8aUs1Rqx13LjbJ7ammiYmzpv+lXT61DfVQDYNNr1/i80xcpL9SPVLsY62ia1pDFfjh/TTVOu4FHlX/0ghCTYky7u+5u2QlBNU7WUgvk9mjJiU6pvdNFJYtTuOykeQvcCk61mlpvE72iW5uhRXSbo2c8ms/j5+/cBfjK6buibVNDdcTRPqohsTP4ZNoU+J9eYD2SxcmCu1znnrl4cP/WsaCt+kAwAAAAAAAPaMh3QAAAAAAADAnvGQDgAAAAAAANgzHtIBAAAAAAAAe0ZwxBXwnIM36imzs9/KhVk0WZKW5jlt6psWTj62iyz7+dzqm4seH4WFlx8zARGSdGO5uVDlzYkf42ByzbT5xSQPzMLLUlo40vd1iykvln5/HqSFYk3/tMhyeKtq7CKvoa9bWzMu/mvaCiEMq3bTFhZenh67kIkQHGH6ruZnFqM/Dgshm4WTJyepb2g3iyynhZe1NH1TQERapNe1h8W+SyqLfU/8QWsXWQ5BDmmBZDf2MvSduJCJcE4u5+H6aQIl0mLKC9Oe16TdPgwiBXE0l8MRFlMOa+b6wccMGrFzDteGqQkvMG2SNAsLJ7uQCBcQIUk3zCLLaTHlB03og+QXWX7w4Mj2fcAssuyCICTpWmi/3jbb0xg+OMJfsA9N30k4iFIg064qNY1Uq1NcrZMCIlxNI/m6Jo3xmDkuHp1v1i6SNIsBV+Z+NEBoR7o9LE2dkoJwXHsPIVKDsDVNuFBWQnPCnN1pUgm4knzwQwqfqPTN7a722D4MYlKoU9IYSuETtk4pXkd2rWtSTZPCIFxIVqhfXJ2SaprUPjFBW8tZuleazyWHflx7jKcaP74lhXrQXKNSWIaraSRJrq6JNc2WbVVhbi4kwoVeST4gQpKmJlDChV5JPiTigZmvR1L94tpd7SJJT5pu1i+pb6o9XK1TqWkqdcqqvwmOSPdVcyCd7vumEOAxFL5JBwAAAAAAAOwZD+kAAAAAAACAPeMhHQAAAAAAALBnPKQDAAAAAAAA9oyHdAAAAAAAAMCeke56BTxr+rCeNn083mZRSNZZmoS0YxuVk9NdXUKaS1KTfBLaY0ufYnYttLvU11nzfV262bSYeLY0qWCuLbUvQjpTSn1tJsknxhq5ZKTw+lolcamQ5pQSnmy6a0pHC+0uIc0lm6X2SUhxTSmsLsm1hb4u3ayd+BcS213i2TzsDJvuGnZ+SuBM/StjOCnxzKV0hb7NJZ6ZBDNJagf++tKnm+3TlJo22+y7nIXtxYhD03+AyOUUROhSz1JSWF+YQQpJhlJIn61cPguXnLRBl44m+dftUtCkM5LQXLqrSXGVfOLZQyGZNSW2uiQ0l2wmSQ8U+j448dvz6a5+jJk291FOR9vcb9NwYEzCwZX6O66mcbVL6ivV6pQTUwPdCmnzKd31UVO/pMT62cIl64YUepNYL+2e5JpqmkVony8299EiJHfbc7g638q1pLC5mFBZSHedmJTScpK9qV9SrePqmklIoZ+mxHnTHtPpXf0S0l2bOZZX/TfbY/3iklxTPZLuzUPUL06oSewYqa85T2JNU6hfmjknJZ9wn2oam/oarrXxFHa7M+xjN3T1c0Llc8kQSa72M1a8wJimUKekWs7VL4chsf66qV+uh7TVlO76JNP+oKlHpFDThL6V+uWBQk1TqVNW7WZ/mvpHyvXLbQfhfRgK36QDAAAAAAAA9oyHdAAAAAAAAMCe8ZAOAAAAAAAA2DMe0gEAAAAAAAB7RnDEFfDUyYmePjl78cK09KFbZDktFHwSnuneMosvp8WUby032/Oij9u3VxahTuICycvN1zdPC04vN+d2HBZzTYucT9yismlxeLtAq+2aVdZ0dn0Li9G7BZals4IjzALJhYWXKwERkl98eXK8fRhECojQiV/AtLl2t2iy5AMlwiLNPQZKFBZkrkgLL7vAExcmIR8c0UwQhCQpLLJs22fhNmf288Sc65LUQ3u8sFomGKOwmLIkuctOH2Ax5cr5PogUbuMWWQ5DTMz1My3iOwvX2pm5bx66tBpJN8ziyzemYeHlsHDyQ9Nbpu9mm1RbTDkFRzzYNufnFl6W/H31ergwT8yBcRgOrrH+T3A69Y7DybM0R5KrXSQfknXL7EtJerSHMAizaPVhYSHraSXkRz74YR72xdwE3hxP/L6YhXZ3rs0XYYF5E8SRbhkp+Kpe2Nw5cBh1pJCsVJKHy4ttj8FXpk6ZHoXzrxBmVQq+SjVNCr6q1C+uBooBEcXwrG1VahpJzYRElOqXk3ClDPVLM7V0T/W1C8ka4t4e9oW77KSaZjk1YVjhdfRlqB3dIVA536v7ovC5y312SwER6TOha5+FWscFSlyfbl/TSD4kYqyaJrU/0EJwxMQEY4SaJtUvM1M9pEBJ92zh9KHs6tAh8U06AAAAAAAAYM94SAcAAAAAAADsGQ/pAAAAAAAAgD3jIR0AAAAAAACwZzykAwAAAAAAAPaMdNcr4CmTpqeaZMTTFiG+xuWSHIe0pJPQ7hJUboXIwcdk0vRC30nIb5sWIhVdeq1LeZN8opvk012PQ9rjUds8pQ5Simsl4XDHYLOykHRbSUaqpKOlJDQXzpOSYF3aVUqCbSH1tZR4dmwSkyoprqF/n4e+rr2a7upS08ZMd3WJrSnR1CWhHYRks3lIWJuZVOkUY+bai0FN7rVMTFqZJPXp9sdnuBTJhEerxcSzzbFT35huVkl+rkjXMzN4iynYm32nxdQ0l4R2GC5GPt3Vp4o9MPFJoK49JZ49NNlMSEt9Y+qrTXf1r8+1z8L75LLbZ+G8nu6ayhmkmibVKe6duh5qD5f6Okt1StieS4WrpNAvQxziIrSfmOtnSqE/WmxeV925sOqbEge3Ty2053BIhhzicCklOJbql/D6XN+QUJmSK929II7h+s63r2kkqR25WqdQv8Saxl8TbU2Sah1T1/RUp7iaZvULvn1bhZpGkmTSXRXTXTfPv5ZqnfA63P5Ife24KY3UtE8m29c0ktQWm/1T3e5utymwN31+KNUp9rPNWDH2tc9z6TOh+wx5EHaGS6dPdUqsX0xq/QMh3dUluT55ctOPG+qUJ5taJyW2PmDTXVOKq9+frq45DG/KxNyQTtc083DPGwrfpAMAAAAAAAD2jId0AAAAAAAAwJ7xkA4AAAAAAADYMx7SAQAAAAAAAHtGcMQVcGMy0wOTsEDp2iIsKLo0K6VfC4tXHoXVPI/N2FOFBWHdY+GwGOgiLFa6NIMcpwWS2+YS1yfNz20WFp88MCubpjAIt0B5Wgx0EsZorn9aBfXCAyXMFCoLMsdFmtOCzG7x+7RY7fYLL8d2EyiRFl52CydXAiIkqR+bxVHDYso2UCIGR4Q5u32XVukdglsc3C2wLNlFlluYW+9hkeWtJxb6hmuO5mFRWXMO9xAc4QIl4sLLMQxiu7bYns6/YvtoCossu+tkCohI1+sDc82/NvXnn+t73SxifFb7NdOeFlO+bkImKgERkvSAudddDwfMdbOfr4fF/SfmTZnJ34OnI6UepZrmJNzHl6b/rVB8TEzA1aS6uLjZdS7IKrUfh2tcCrg6Mv1vmvpH8sf4LRMmIdVqHVu7KJzDF1y7lK9xrtYpXGtL12UVg69MSMQkBUekdlMjlOqXFBBRCL7qaXsuOCKFYY1Vv6Rgk1QjuJCIadj3pg7rqa71W7Pt6VB2AVcpwKuZ+qxcR9u6PczN1im1gKtSUMxYKmFYxc+E7jNkpdZxn2Gl/Jn3mq0b/Ll63YVThfonBV+5kAgXECFJD5p5XAv7PoVBXDPndgq4mpgb+ema5sZdQjt3xTfpAAAAAAAAgD3jIR0AAAAAAACwZzykAwAAAAAAAPaMh3QAAAAAAADAnvGQDgAAAAAAANgz0l2vgGua6Vq7S7prSDFbmgickx6SIWMa1+bYblxJmpkkppQwMzMJa6v+m+kuh2kM03caYobSPKYmWSclvbkUnpjuWkgtvNQGSHetpDaFtyn0TalUIXXLJQYuUzqaaU8JZDFt1bSnvjbxbPu+sX9ISSyJEZzbj11KKwvt3aQ2tZCwZt/X+P6FmZgx2jL8v6+RzofKGPfapeUslXDI0vU63bvMhSeP64+jQ3c/Stsz97/DeE9M29tsdymuqzFMYmtIOJyZesOloK3a/QanKT3RWJi6Id1XU0rbidl3i8LJsww11EloPzR1VKpp3DHgjhUpH1vuWEw1jesb65RCTOJFh82PabRr7QD1kk24D3VKvHe5e10aw7THtNWUTu+SXGOd4rYXap2U7rprXRPTT/1nLbe1eD64sdMHgkVKmTW1Tkj87K5OCfvH1bUtzCEmJpu3qnQ+DHHuXGIx3TUcAqXPlWZnuM+w0lmfed3ndH9eu8/eqW/8rO/GCPfVmdlHlRSVx7rNAAAgAElEQVTX1dib57CraSRfv5yuXa6N/F03vkkHAAAAAAAA7BkP6QAAAAAAAIA94yEdAAAAAAAAsGc8pAMAAAAAAAD2jId0AAAAAAAAwJ6R7nqfSElqS5NANg1JKSchQafypNelzEyHSJcMfCpOSINKY1SS0C5zfOJVil5zXOJZ1a7HYvj9nsY1c859XRpb2N4Q+6Ii7jcz5+4PxG7ay4esS3pLiXV+EtUtbi2lqV0Z7n294OvhZbn+psRW2zfcj1x7unel7ZXuzTveIIZIcU3sGCHVsXJ3z0mwLvU3jbF9wm/lvU4qx9aYLvxcs/HfV7uoGfWeURnb3UNTqmppCmEObuyxUlyLUm1lEzvDee1et03sDX2lAY6NIfbbRdeZF+0SX14qn0ErCd2rMbZPUx+Cu6elRFqnWmGkZxy+736/y8Y36QAAAAAAAIA94yEdAAAAAAAAsGc8pAMAAAAAAAD2jId0AAAAAAAAwJ4RHHEFLLTU4i4L6C7DwpFLs8zyiQmTkKRlWGj0xLQvwpqPbp6LsDrnsvtnyAvTvgjPm0/6dPtxwxhLM79l2N+p3ffduuswLnp7lUVXQ1+7O9O4E/MH6X9DuL6S5BYUDYuMdjNGm/oNtnkYw/Rvi9R381hOi/82mb6S/Km9+wLQCourNrefU1+3nyfhDUyLubp9FMfY/r2Ox1HleLF9w7hDuOhFj91Cvxd8zYnX5cLOSH3T/cH3DWO4e1fhPufuZ5K0aP6efWLmMQtvysKGJYTFqd11J8wh1Q0paMKOYebhahdJWoTtudfn2lbtm9y+XPVNQVub71XlvY59C/utUtMkqW+l1hmEe6vi4uIjze2Cr9f2niFd+L3L3kPdvVaSFqGeMGPYe76kbu7vLWwuhmTtGmyR6pT0ul2tE/ra151qx/SemP6uJl0Psl2bNFAdvdkULxcDfE4YzYXXL5W+21+XqzWN+4wcaw8zRrp3xc/e5lhMzzTcPdsFPUk5ZOnE3OFn4dhygZmna5fFEJ+fzsA36QAAAAAAAIA94yEdAAAAAAAAsGc8pAMAAAAAAAD2jId0AAAAAAAAwJ7xkA4AAAAAAADYM9Jdr4CjPtfRqUSzlHrmnJgEpJR45lJcJemWab4VUlwe65uH3C3TtmqfhXmYMZap72YiTeo7X4aE2OXmGPOYXrMZEZPG7SG9xrZfdJJaUklbLfTNKWam79Qfhy7ZKqVdxfYD816F968tTNpRSjY78Me4m0XMq3MJsSm5dOGTFl36bA/ntY2aSmlegU8xC4lnrj0lqaX96doPQtKtaw990352x0s+5kxjOh/C2+raY9/CZeSyXF6cdHi65nRNTUloc3NtT6liru/R0h+HR+Eec72dbLSl+9yszzfapiGxcBYS2Z2pNseVfILq9VBLuL4prXU6UiRfqlNiCr3pfysm1m+23Qrpdrl9831N77VrT2l66dhyx6I7ZiV/jKe+6dxx51q6d6Vz+DKoXBNj38p1OaatVsYw+97VLpLaPAxi7nWpFmjLzetAOqtT2qrtn+o+V7+EmqZ0cKX4TLM/Y6pqqrlMreLqLUm2frG1S+gb5xHqJVunhONlkDraTWOI1PsBPmtcOHedrHz2k/8MGWsal8waru3ps6n7jHxtslm7pL6uzpFqtU7m+qZnHiEB3iXEptoqpDzfdqtQe50H36QDAAAAAAAA9oxv0t0DWmsP39m0l4kAAAAAAADgtte21t7yVb3e+0O7DMY36QAAAAAAAIA945t094A7n8S21p4p6fV7mg4AAAAAAACk5/fe3zDUYDykuwIeXp7ocHn24oXpT91SiW7RZEk6DqvYnpgvZFYWU35sec32fTS0P7y8vvUYjy0226uLfR+54IhCyERaiDwucr50C43arnml5l21GF+wOYXKwstpIeSwXr9b77RN/QYnJlBiOQuLKYcFhG1WQgxW8M12e2kRYrOwsA1QkNTdwsRpMeUUYGEWR037oiQFSrhFV1PfiXl9KcghBUrMNs9tGxCR+h76a0OfpffELOgbjjnXvpylhZe3D5+oLC4evzt/0YESletZvE5uvphFuC6n9mNzgTmO94fNcy3dS1I40axt3o8m4Vo7rVxgghOzgPNJWAj5ett8fbdaWEzZzC0cylE4K63K0sypfjkxB3kKZ3D1SzXgytUvlVon9U3HljsW0/HpjnF3Lki1c8qdk5JCMeC7jlXTxGtZarfX2nRd3px0rnXC9cydQCm0ypyW6T6e7l123NDuRm4hQCHWOi6oYh4CrtznmRRIERZ8jyER26rUNJKt5WJIViHgKoZBuPrFtK3aN8dYpprG1ikhZCLW7Zv7Llxe7BiVMKyz2q342WZHKeDKtLvPeKu+vt19hqyEQbjPsKv2dC/ZDGdwn6UlaWpe+DTUDRWLyZFvd/fxUCHcCteGQzO/WTguJmaM03vzEROqMyT+uisAAAAAAACwZzykAwAAAAAAAPaMh3QAAAAAAADAnvGQDgAAAAAAANgzHtIBAAAAAAAAe0a66xXwh8uplik65y6WJjXmODy7TUloJyb1LCWeufaUVlZKdw3JM4+YdtcmSTcXYc6m3bVJ0snCpLuGhK6UEGsTfkppiKFvQUxRclOO6ZKmrZDiKkntYHMiKThoceg2mJK//KTdNJbNXyYnLsUsJLfF5C6XbhYSW9vczCOmuJaipnzfIbgEuJAKZ9PUYtJtGMMlyMVkVtc3JLMe+jFc+9Idh5K6OZaXpm3Vbptt/3RO+XRC37WSfDha4qukblLP3D1K8glp80VK3d4+jfvW3F/bD8yFZxIuRimx1XFpZVK4NxdS0yXpQZOQdt0kvkrSzKS7HobUNNfXpbxJtX1RkY6LtD9d/XIcMmZd30qKq+TrmpTYumtNk9pjTWOO8ZQWmM4dd66l1EL3XrlzfSj20CgmQ5ZqHbPr0rU9pbA2sz/MabYaO13I3dxSmLppT8ms7cBs7yRMLtQ6zdQqMZ3e7aNUp4xVv6Q6JaXX2now1QJmH6XE+rA/u3lPUq3j6pRU6yxcums4lm3NHfrHWselu6Yk2HTYD1HrVFQ+d5n2lOKaPhO6z5Dus6Yk3ZpsXtsPJiGlPdYvu51T6R68CG9K5X7r6pdc02ym1K7aXa1zvlruj10S9YD4Jh0AAAAAAACwZzykAwAAAAAAAPaMh3QAAAAAAADAnvGQDgAAAAAAANgzHtIBAAAAAAAAe0a66xXw+vmTdOySH09ZhuexLoUlJ7D4bdjUtJhCd7jRdhTSXR9bbvZN7Y8tfN+bpv3R0PfRedieab8ZEgCPTOLOydzvi0VIIlya9piENkRUUSUJrZKiZFOb/MD9wCfouOCuFDzU7L5ICV2+2c1vkhK6THsL6VqTeZi0aY+JZ4tK4lkh3TUlwQ6hkHjm+vaUsObS5uTfv5zualLMUrprSiZzSWiHqa+51pq2M7dnLsGVxOSYjlZJOKxccmLi2fbXs3Ttc9fPeTi20v1xapK7JgPEY6f0RXevPApRvkfTzXvMAxOf7Hlt4tPNHnZJaKGvS0I7DPGSrm9Mdw2pabuq1DSSr19SnWIT60OdklLoXF0TaxpTk6S+j863T6d/5CT0NfXLUThH0rnjkghTTWPP4cI1YNXumy23uQGufTGF3qa7hr6xlvPNtqurw1rY96k+M/fQyYk/V5upa1tId22uTpHUXa2T6hdXk1ST7HdVSXGVfFxuqF9snRLTXVNNYj4npL6uTgm1h0uhjymusX7Zrm21PdOWPicUUl9L5/sgia9hzi6xPlwn0/XTfYY8ComtrW1e26sJ60uzQ0o1TXhWcGSSZyXpMVPXpJrmejveaMt1SqF+CR8sXV1zuqZ549wnyA6Fb9IBAAAAAAAAe8ZDOgAAAAAAAGDPeEgHAAAAAAAA7BkP6QAAAAAAAIA9IzjiCnjd/Ol6+MQvyHjbovA81i0EKdUWjkxjuMWX48LZsX1zDLdocmq/tfDj3kpjmEWWb4XFlI9ONttPzKK70hnBEWah0Rgc4YIVqmESlYVUCwsyu7VD0zrBebFot+BpWKDVzCNlFLjQAEmazDe3F0MtzGK6aeHlZQhyaGZ7LYU+mPYWFlOO68RehuCIsPCyXUzZLcZ8RrvMGMsU/GEWSF7GYJPtwyAqoQ8xOCJczl17XgDatKWFl0O7PZCGWGQ5MYdzuva5xZfnIaTnOJwQrbCgsl1MOeyM47DC/JFpv2kCIiTpkcXmYso3pmEx5RgGsblwclqQ2fVNCy9PzRuVFqd2fYeQappluJe4/qlO8QEf/n2q1DqpTnG1Tuqb6hTXHgOuTP2SaprjcE65cy0tiG7P4XEOixVb04RAg3TNd9ei8InJH3O1+2q398rQ19z/+jS8vnQ/MrXHJAQEuDplEupaG3AlX9ekWseGT6Q65ZIER7hDINWZLmQi1jrh/XP9U/CVu0Tl0Ifta5ocfOW2Z7uG7fm+6fwrBVwNEBzhPmP1cCzb4IhCGJYknUxccMsQAVd+Hi4U6Gjqd/4NU7/cnPrQo0dC7XFtshm+4NokX+vkgAjfPjEhEamvc7qmeXjuX9NQ+CYdAAAAAAAAsGc8pAMAAAAAAAD2jId0AAAAAAAAwJ7xkA4AAAAAAADYMx7SAQAAAAAAAHtGuusV8LvHT9ODx48nwS1S3Kbh0l2qqWku3Wwe0u1OzNxS35Tu6pLzjkNiq+t7FPq6xLPV2Ca9NiSeuSTXk5NaumtfmP3s2qQQKeW7VsSwVZMgFtMX3TzODiHeZJK0+iSltLk2P7cQHKRm9vPEBxXZBLJJeJ9a2p5NPPN9XYqZTUHTWemuod2OsX1nm0wn+TTglKTlEs8K6XarMTabUjKZT8jzw6bUV3eJSsecTTxLaYFpzqXUNNOW0l1TCrJpj0dFJSEtHeMu2TqcU0tzzC1aSEdr/oW717II17MTk3iWkrtTEtpjJvXsekhsPZxspo0dmDZJuhbaZ26MQhLawcS/UZV012QaL3Sbdq1pJF/XuBQ7KdQ0KQm2kOSbah1Xp6QU13jMmbrmKPR1ia1p3JSYvJhv7rtlOFftOVxIrI/S/cG1Fa5xUrhWlg7x7VNAV/MwtVWhfknhhJOQwDkxtUMaw9UpE/P+S4rvn6110v6s9A3tg9QvduA0hmkL75+va8O4MfV1sy2m09vU+7A9lwRbSKxP7TFluFLTFOYc046HSKe3F5gwsLn2pfd0Mdn9e1PdfaZP97l0Lzkw6d9T/yHmMfMB6TD0dTWN5JNcU60zM3VD7pvSXTffwJRCf7e65rGjozP/fFd8kw4AAAAAAADYMx7SAQAAAAAAAHvGQzoAAAAAAABgz3hIBwAAAAAAAOwZwRFXwO/cfJquP3bjzD7LARZeXoYVU93iy2kMt/hyWrw5LXbpFlmOC2Oa9nlaIDkEObj5xcWUzSKhy7C95TwtsuxWv6+skOy7Ju6tCuuv+4Xkw2Kubh5xTegYVGG6psVx3QK0KVghBTyYdUbDmqRq5j2JQQ5pMWXTvRQcURg3GiBopBIaEBfuLYRMxOAIN0ZcGLyyeHNotwsWF7aXFkJOCzK77VX6xoWXwyLL7nVXFl4Ob2A6Prs7ntNi9OYiFU5Vu5iyJC3NObyYbr/I8vHUb/FWCI6YmiCG2eSa6SnNzNgpbCEFPBwUFll2CySnRZMn5qIR+xYCIipSTVOpX2Jf055CH+ZhHq5uSAEYNnCqWBedmLohLQzuQqtSkNUi1DouJMLWLpI/h0NN0+JFf7fF/fM1PF2MzPFS2V7KNorXYLO9cG33oQ/b1zRSqF/CC7Tbi3VK2t5uYRDFXJrd65piwIANjkgfuyrBEbEGMm2pNnYfKUohDKlvuH4W6hQbVFGYW2qP+97tt2qYROVzlw3DCveuuD0zRtieq19iQFK4GLkQxYNQ68xM7eHqHCkHR7j+uabZHCPVHmmMSv1if//Uxe/WrZtb/9558E06AAAAAAAAYM94SAcAAAAAAADsGQ/pAAAAAAAAgD3jIR0AAAAAAACwZ6WHdK21P3XeDbXWPvS8vwsAAAAAAABcZdV0159urX1q7/2bt/2F1toNSV8h6ZMlpRxI7OD3Hn2qDmcPnOt3U+qdU0lCSzkpLpksjbsI6V9Lm5q2fV+X6CfldLNu+rtxJZ94llLMekjjKiWhufZqqpVJtekmfUqSPYOXYYPNjFFJjZWk5tJdU5qlCQ6y+0dnpKLaMUJfk1YW0ywrKawxHa0SoVrY3ohKqVmFhLyYvFZJ9bMpXymSzzfbFLIBkpFLcy6MEVNc0515x3TXUoprGnserssm3sxdq1ftfoPuXjCfhJRSkxR2PAkprtOQKmbGTqliLvFsGubWCmOk7bk9lxNbTVLjRV9cgpjk69I6Q1/3SnKdko5Pkxycao9CXVSpX1JfV7/EcyelK7t6qVTT+K6pvXR42YPZd43XPhsxGvanGTvWOun1zUxbSGa172tMZvXtrn8LkZGVJNhYf5rjuZLumlyamqZSnhXuq5V01zjGpJAmW0hKHWuMSorrqt3UAqmv20fVg8hdGtL109Y0YdiYWL65QftZU9LS7ItFqKGOQz3h6hdXu6zaTZ2S6oaQturqmkpdlE69Sv2S3K2uOX702tZjnUf1r7s+SdI3ttb+WWvtrjNrrb2HpJ/R6gEdAAAAAAAAAKP6kO6mVg8tP0XST7TW3jl1bK19hqSfkPQu69/5w/NOEgAAAAAAALjKqg/pXqTVN+OapHeT9LOttY8+3aG19pzW2vdL+lJJh+u+3y/p3XefLgAAAAAAAHD1lB7S9d5/VdJfkPSFWv2t7CdJ+qbW2v/aWru2Dof4fyX9Za0ezh1J+geSPqj3/vuDzhwAAAAAAAC4IqrfpFPvfd57/zxJ7yvpV7R6GPepkn5d0r+W9Ix1209Jeq/e+1d2t8IzAAAAAAAAAEn1dNe36L3/xDoY4gskfYak52r1cO4xSZ8n6St77ylXCAP6o0ce1MH0wXP9buXxaUpNqyTEuu3FcVNaTiFN1o4REwcL80iv2YXXpFTO1F5IDrLbK0VqqpSIue3vr8YwCYAhcSknqG6fYmbHiEmpvr3St7n9XEwrK6W7FvqWjPm/UHZMd00qCWulNLZKkloau5LSVk1Nc0loIXXLjV16HXGMkBBbea8LCcY9HqDbXxv6wl947GsJr2/pLl4pbTWmsLq2QmJr7Gub7RiVvkmt79ZdS6r/67dWp5jjMNUNsS5yjYWapjKuQp0SayjTWK1TKrWHSxIt1D9Jvm65QbZP8pV8OmsL6diV15drne3GlQaoG9L2YpTo9uMOkiw/Vt+KIa5blVqgmO5aql92rPFj30q9FGsak8xarItsnTJAKm7kBl+m1PTt7yXxer3Yvk7xodt+x6U6xSbcp7rBfc7zXXOt48Yo1jql7Q1Yv8wfGanAWSt/k+4Oz9PqG3W3X3HXKlzidTygAwAAAAAAALZz7od0rbVPkfQfJL23Vg9Of1mrNeiertU6df+mtfbWg8wSAAAAAAAAuMLKD+laa09rrX2HpH8m6UGtHtD9U62SX99bq4d1TdKHSPrF1tpHDjddAAAAAAAA4OopPaRrrb2/Vumtf0OrB3FvlPThvfdP670f995/QdJ7Sfq69Z8/TdK3tNZe1Vp7yqAzBwAAAAAAAK6IanDED65/Nkk/Lulje++/c7pD7/2mpE9prf2gpH8u6SFJL5P0FyW9/W7ThfPYm69pquvjb6iyMviYC79WFuyvbG+I12cX2C0ukFwJCDBjlxbo9UPEhTX7pLDvC4vVxj1vk0bCGEMEK5QWSN7+eBlt0eQxx7hoIy3UXAo0KJ48lQWgdw61SGOErrsuLL3qv/0Cwn7gMGxhQfsY8mIXgA7jxtXT3cXPdy2dwnF7W7YllXHjGIW+1XlcZtVApY3fH6D/SDVGeewBtlcKTiqEFFQDCfwg5tdTwE5aGdyVHulrDfYymRI+CmMMUCPWao8Bahrqlwv5/V3rmnLtUei7c11U3F6ttirsi+2HiMGDrqiJoTk5cWFz2AFqjHwr2bHuq9ST1TF27Vudx10s3jwfbCyn+tdd2/qffyTp/e58QHda7/1bJL1Y0s+sf+d5550kAAAAAAAAcJVVH9K9QdJf671/Tu/9ro8Pe++/Iel9JH3FeSYHAAAAAAAA3A+qf9313Xvvv1/5hd77iaT/rrX2A8VtAQAAAAAAAPeF0jfpqg/o7vjd7z3v7wIAAAAAAABXWfWbdG/RWrsu6RMlvb+kF2i17tyH9N7/0zBTAwAAAAAAAO4P53pI11r7O5I+X9Jb3W7SKijkxh39ppL+iaSXSPqw3vtvnXeiyCYPzzTps31P4+q6F5OmdhRTzMzOuCy75+IDBy/LK8e9atdAxntWKYkw7KTF9ptrQ0Sa3q/vFa6eK37rKqVLhp1xGXYRNQ3uRVe+rqkkMY9UOFDTXA79kXGfvVSDI9Ra+1xJXyXpmVq9xY+mvr33haT3kvRukv76OecIAAAAAAAAXGmlh3SttRdK+h/W//lzkv6ipKfc5dderdXDvJeUZwdJUmvt4dP/SHrtvucEAAAAAABwn3vtHc9rdlL9666fLGkq6fckfVDv/fWS1NqZ35n8zfXPtynPDgAAAAAAALgPVB/S/RWt/jb2l91+QLeF/7z++dzitrDWe3/o9H+31p4padv9DwAAAAAAgOE9v/f+hqEGqz6ku/1tuJ8q/M7tNeueU9wWtnTwcNPBcvwVIAdZDHSkMUpzq87BbW+IsUNfv+hxWml9y9+vzmOkuZXmENp7ZTXlwriDqK70vOs8rvrCr0OsY10ZY4iLXGl7vrlV5pG259pj3+23Fw/xwvZqY4w0t9BeOgKq2zNKl4zq+bDj+XPxC9dfrAuvaSqn9Yj3cbu9yrihvVaH1eqGUhjEIHVR4eKw476XQl0zZl1rxxigtiptb4AxLrN7sX6pzKF0n0sn9m7j5jH89uwhPuZ9fNe5FedRqmlGqj0GeR0DzCMZsq5pj457EasGR9x+qHdS+J13Wv88Lm4LAAAAAAAAuC9UH9K9bv3zXQu/86L1zz8obgsAAAAAAAC4L1Qf0v2wVt+c/KvbdF6vnfb3tPqC4o8VtwUAAAAAAADcF6oP6b5p/fNjWmv/1VkdW2sPSvoaSU9eN31bcVsAAAAAAADAfaH0kK73/oOSvk+rb9O9urX26a21O5NHn9dae6mkn5f0wVp9i+5neu/fO9CcAQAAAAAAgCulmu4qSR8n6UclvYukfyzpi0/92c9IeuDUfzdJfyLpY887wTGtv+338Vr99d0XSHpHrdbd+0VJX917/37zO4eSPkrSh65/550l/dH6d14l6Vt67xeaiTZ7uOlgcc6EkYtOWx0p3aynx82V5NLCGH3i32I75zBu2l6z+9NP2o3RQnTNEK/Pvpa4PyvpaIXtpbkVkmdbYYy0P31f3zXOw7WVEtZql5o4vwtUujqGi0AvXBzi5swYcW6hvbtU7Thn0xhSueM87PZCXzePpe/aQrsb275mhTDE4v5sbuw050qabOV1p+Q12zeMG9+/Qt+R0tRGS5O9LEmwF522WkkjLcwj/n6qJ3atdYp1g689Ql9bF6WaZvt7ZQyMtNsrpsm6/rGGKtQ6lXmkfTFareO7urHLYbmlG4T7/eIGRzJE/eLHDedDYexK/RJrqMI9P5+AhTEK9+BS3RB2Rnf34Pj60vZc3+3rhnpdZBoLddEwr8/3tTVGqieDXWuSfabJtsfGvTCVH9L13v+4tfYXJP1vkj5C0uGpP37wju6/LOnDe++/cf4pjmP9sO03JD173fQGSb8q6U9rlUj7Ia21L+u9f9Ydv/rjkl68/vc3rX/nT0n6oPU/f6O19jEX/aAOAAAAAAAA967qmnSSpN77G3vvHyXpPSV9laRfkPSwpBNJ/1nSd0v6BEnv1nv/5WGmOqze+7Gkn5D0tVo9ZHt27/09JD1d0tevu728tfZ+d/zqj0r6TknvLunpvfcXS3qqpC9Y//lHSfqbI08fAAAAAAAAV8h5/rrrW/Tef07Szw00l3342N77zdMNvffHWmt/R6u/Avu2Wn1b8DWnuny2+Z15a+0Vkv6KpD8v6SMl/ctRZw4AAAAAAIAr41zfpLsq7nzYdqr9WNJPrf/znbb8nS7px9zvAAAAAAAAAGe5rx/S3cUL1j9/e+TfAQAAAAAAwH3uCX/dtbX2dpLeYaCx23r865KuSbree/+mgcYeVWvtAyW9cP2f//uWv/Mirf66qyR9wwjTig4f7jqYnzOnYqR0s5RG6pNEtx839a8ll24/bt5emJzrO03jhvQh17/wOD1tL/YvpcIVEs9c32lKcQ3JZKZ/Sitz7ZNi4tnEjTHxUUXuEI99bas0Nf1Tums65Oz2CnFJaR8NYbljutkypaoWUl/THFx7GncZ0rFce+rbl5snSnwdIR2rL8wYKblrYdqKIW2V1LTKzSSmjdnktTCGeX1pXNc3jV1Jup2EcSspuqUk2ErSreSTyWIirYvkK25v2zkMoRiqVrrPuYv7AOmusU5xdUM1hd60L1MtYOe2fYK85OuMmArv+saphXmk2sEOsn2dkpPsC7VOpU6ZhhrBpuWm+mVzDFe7VNtz7VGorQrtqaa5DPVLpXaRaonzrq5Jv78wdYPkz59l6ms6p76ufkk1RtpHrn8aoy9conuqacIYrjneK83ctv+4GlUS5Cs1zar/5kxS31LdMMCcK+muQ4zh013TB4UwjwHrl8nNcTNC71yT7hMlvWKkbXVJl/YhXWttJun5kl4m6eXr5s/vvf/kGb8z1Wrdug+T9N9LuqFV6MS3jztbAAAAAAAAXCUuOKL4/yfvba21p0r6kzuaf0fS3+29f88Zv/dGSU851fQnkv6WpP+j9/Q9iK3m8/A23c47PgAAAAAAAAbx2naXrwL33h/adrA7H9L9B0n/Yovf+6uSnrv+95/Vag22Y60eHj1N0ntIenV5DG8AACAASURBVKZW3577J5J+eNsJ7cFcq/TWmaTnSHp7SW8n6Qtba5Pe+3eF3/sRrR7SPVurvyL8NK2+hThrrX3DOkgCAAAAAAAAuKsnPKTrvb9a0qvP+oXW2mdK+nhJvyfpo3vvP2b6HK77/C+SPk3SD6/HvnR6749Iev/b/91aextJf1/SZ0n6t621v917/xrze//1qd95uqRPkvTFkr5O0jtK+ofnnM9dn7C21p4p6fXnGR8AAAAAAACDeH7v/Q1DDeb+umvUWnuhVg+ibkp6v977b7p+vfdjSV/TWvt1ST8g6etba3+29/4Hu054bL3335P02a21I0mfJ+nLW2vf3nu/86/Env6dP5b0Za21P9LqId3ntNZe1Xv/lYuY8+yRrtnp4Igd/zJsJSBi1d8sxJkW4x0tyCH0nbq5pb7bt6e+dqHm9JefzdxWg5u2g8L+jG9gYRXNwuLUMQzCLaacFk0OY0xN/7Qgs+vr2iTpIAQ8uOCHadjebLq5cuu0sPDyamyzAHTqa1ZSrS6aXOnv+lYXU3YGCXIIF6O5WQw5LbzstrcI25sv/Am/MAv6ujlIfqHmhQmCWLWH123OE7vwsqQ+MSET8+Ji2KYthj645kqggfylKC2QPDGvOy68PA/tboHkuD3TGK7tKVCiLcw1sbLI8hALQBfek8qCzCPmz5QME3BlzrNikINtH6mmkXztYS4BsW8OuArtbvH7NDd7wGw/7uoPTFvlva4EREi+fgl1mKtrUk2TwqWmpn+qX9wYqaZJ9cvBjvWLq13OmsfEvIFDBFVUpDEqdc0QfRfmpIp9B6hfTkz94moXydcpqabJ9cv2AVcLV6fEgIhwMbL328Lnq2oYltnP6fC09+BQh8X6xYZk+b6u9qiMW93ezjVNaB8kqOIC6pd2a9zCp5AZKUn6VK1yml6ZHtCd1nv/IUnfKOmpWoVS3Etetf75JEnvueXv/CtJJ1rt15eMMSkAAAAAAABcPdWHdB+g1TPIHyr8zr/R6nH0hxW3tW9vOvXvKcj+Tje1Wpuv8jsAAAAAAAC4z1Uf0r3t+mflr63+7vrnOxS3tW8ff+rf/58tf+elkh5c//tPDjsdAAAAAAAAXFXVh3S3/+bvCwu/82fXP59c3NaoWmsvaa29qrX259ZBF7fbD1trny7pf1o3/aPb69G11t6htfbq1tpfaq3dOPU7B621j5N0O2DiX0n6+Qt6KQAAAAAAALjHlYIjJP2WVg/o/pKkb9jyd95n/fP3i9sa24dL+rj1P/PW2q+t21+gx/fLq/TElNYPlfQh6396a+03JD0m6V0kXVv3+T8l/e3e40qUAAAAAAAAwBNUH9J9n6R3k/Sy1tq39N6/56zOrbX3lfRJWq1j95rzTXE0L5f045I+QtI7S3pHrTJiflbSL0v6p733n77jd75S0q9o9WDvBZLeSdKzJP2ipF+V9HW99x+4kNmfcvjoUrN5ikJZSQldlSTYPkmJOyYJLW5vcwybiHrGGG4eOW3VzK2Q4ipJy4Ptt+fOqB7OsviOuWSy1Nnso5RSU3lq7N7TOHhID3OJZ5OQjjYxSWOSdHCwfWKrS1s9KPSVpJlJJst9zfZS4lmIUXL9D0JUkU88C4luKTWtcBRU0tRKiWfhouOSyebh4lBJd01jzM0FxqWgpb6p/0lIQjueb/adT0Jq7NSPMZ9vvieLsPzp0lxheryoFpLQdg/4zemuNuVr+yS0lOIa01ZN/5juavuG86ySmhbTZAupaWl7SzNGYd/nFF4X7en7JjG9zU2j8nc+Cin0aWfYtNWY7hrqIpegGtNWt68xekgNbTaxNUza1Snh/SjXL0Zz+94cm9IZ6bwVbow0bjq2Cun0E5fMerB9TSP5tFSXwCpJh2ZsV7tIuX5xdUal70G4yFXSXdMYrn6p1C7S5ahfluGEd9tLNUZKVp2bC0Tqe2JqoFTruPol1kUh3dWmyYa+7hq8dBcz+ZpGCnVNup6520CxqLGHVjrcCsmllWT5WOuUappwLynVKdv3dfXIqn27cSXZJN6c7lqbh9+e+/1T/35UuRPWVR/SvVLSP5A0k/SNrbXPlvS1vW/e3ltrHybpq7UKUOiSvnbHuQ6q9z6X9K3rf7b9na7Vg8rvG2teAAAAAAAAuP+UHtL13n+7tfYKSV+i1Rpzr5T0Ga2112j17bObWn3D7CWSXqzV/7vqkr659/4jQ04cAAAAAAAAuCqq36RT7/1LW2tvJekz103vvP7nTre/P/pDkj75fNMDAAAAAAAArr5quqskqff+WVqFKPyiVg/j3D9vkvR5kv5a7/3WILMFAAAAAAAArqDyN+lu671/V2vtuyW9h1Zpr8+T9ICkP5T0U5J+oPf+yCCzxJmmNxc6WKRVFldKi/GmhZBjcIQbI8zDjOGCGVZ9wxhmTdEUPtEKoQ/LdDaYhSp7aaFZz62lLGkVX3Jn37Dv3eLLpUW2pZ0XWW5xIWuz+G8hIGLVbhZILiymfC30vTb1q666sQ/dSqySDgt9Z2FlU7dwcqXvNBxdadFk157GGMLCHDBpkWbXfhJOVrdosuQXSE4LMh8tNk/4+UEIfTB9JenWZLPdzUHyi2Efh9CVExMyUdVNoIRd8F9STwvpukMjXri2F0fYcZHlSkCEFBZZjn3N9awwruSDJnLIxG59V+2mrbB4s7v3SWFB5jTuSDn38RZcqVPCsWyDI2JARDinTGhDqlNcfkwM1AphEM1copYhqMntjPBWx7uDu1KGjBc1tz/TeZ1ed5iHH8S1hRFi8JULjvB9XUhEqmliEJUZw9UYknTtYPMC44KsJOl6qHUOTXsKrXL1Uuqb5uHCs2ZhhXkfkhX2/Uj1i6tdpFr9ksZwYVYnMTgi1Dqm/3H4EHNsghyOp76vC304MnXOqt3PrZn36iSd2Ea4farHNB3zOShdjGxNs9W0Tg1uhkjXMzOPHKwQ2s0OqdQeLQRElMYohE+U6xQXHJFCLVydMURwRCFM4vS4y+Ozn73s6twP6aS3BCn83PofAAAAAAAAAOdwrr/uCgAAAAAAAGA4PKQDAAAAAAAA9uzcf921tfZcSR8s6b+U9ExJD93lV3rv/S+fd3sAAAAAAADAVVV+SNdam0r6Ekl/v/D7TcV1XwEAAAAAAID7xXm+SfeNkj5K9TwUjGR6a6Hp4lQcSyF9r5R4FlLF3JGQktBckmub17bXTWKZS3GVfGpMTnGtJLaG5Bmz71Iwa/rL5i6UKKX+jPbkO83ZxeWECJ1mUtMmIUltOt0+Cc2luErSdZN4llJcK4ln16cntu81E2F0LcQlXZv4MQ4mm687JZ659tR3Gg6YSTqQRrI0CWSLkGK2NAddSndNCaqu/1E44W9MN/veXMxs35Rk5/ZnSlhL6XQVLkEuplWbC0lMqAznpbsX9HS+2xtBuJCkzRVS0+wYhSRYKaWmbZ+EFk7rPIZLpE0pZi7dNfb182jzzR2SU+jMfTVFfoaUNjvuRae7Jqae6DHd1Z07YR4hEdqlvqY6xZ3Dqe8ypdCVdsjmGO76K+Vy0u6PmKZn2tLLSMeLfX0hWddeSMK4Kd3V1i/+5HH1yzT0dSmukk9yvTELtYepU1KyfKp13Bg3Qq3jkllTrRPrF5tkH9Jdzft6mWsaydc1sX4x7THJfhnS6ZebtcrR0u9Pl8J6a+H327GJiXZprWe1V3RzXrvaRZKW4Vz15/v2Nc2Y7LUvHLIpFdW153t+IYW+UL/EMSp1SiH1PqXQu7FT33g/cvVLJcn+lOnJuNef0kO61tr7SfporV76QtK/lvStkn5b0h+Kb8sBAAAAAAAAZdVv0n32+udc0gf23n9o4PkAAAAAAAAA951quuu7afVtuS/nAR0AAAAAAAAwjOpDuievf/77oScCAAAAAAAA3K+qD+l+c/3zjUNPBAAAAAAAALhfVdeke42kF0r6c5J+fPjp4DwmxwtNeoh6uS2lmLnmlAAYUlhd/5Qi2ExqzzKluIYUpZj+ta2U7hNiXNxraeHxtkvcSalwMWbFtQ/RdwgueTYkLrk0y5TumpIvp6b/LKSmHZoIJJeYtmpPSWibcUcp8cylm6W+KcXsuolXGiLdNSWeTU385Vipo8lJ97edRSXd1SSQpf5pH7nUV5cqd1Z7hdtH82m4xoUxFu76Gfa9O9dCaFpOcCwEobm+5Rw1N0glMTIml6Z2kxSWEtZcumtIOY0pbS417SQloZm+YXuurxRen0l8TX1zOppJja2eIpUb+a6J9ZK6OddycqmpaVKK69Lvz+7SWVNKey3CLTS7hMPt6750H3c126p9uzYpvNXx4rL9cVEKtC1e+1x7rF9M+0GoPVxifWp3iaiST3JNKa6pJnHtKYXe1Tr1dNfN/jGd3hwDsxQNGexa16T76iJ8t8Wlvg6R7nqrhcR5UwdPQmLrJNRL20rJzy6ZVZIWk8194dqkYT4nuHT6cvp3xVif0VK7u9amdNdC30r9EvtWappCCn2l1onprqHdptbHvqbx1O9PTmrXparqN+m+SqtU15e31p58t84AAAAAAAAA7q70kK73/muSXi7pOZJ+sLX2vFFmBQAAAAAAANxHal+yb+0hSd8m6e0kfbqkX2+t/VtJvyTp6G6/33v/ovNMEgAAAAAAALjKqn9R/XRgRJd0KOmlhd/nIR0AAAAAAABwh+pDOrcU47bLM465nP19rc2XamnF3rUeVse1AQhpoeC0CLENfgh/k9otvu17ahlWWe5m0nFBSjN4fXHxQl+zIGWLK1n75pIhxhhg8WXf1+yLGBARgg5Me+rrxj4IB8ZBeANnpt21SX7h5EpAxGqMcYIjXECEJE3Na5mk1cwLlilJxZiFkBu3cLJbQFrygSKSNO27vZa0KPRBSH9x7bnv5tzc4s+SdLLwi0i7PJ50TtnLdVpMe4DwkJLK5gpTrr68yrXdzaMSSJHaS33TYsppDLMgcxyjEgbhwhIK98QhlGoa+dcSFxe3gS7hBYZACZnLXMy/GCQEZaTQlQHOvyHO99EUXmAOmTCLp6e+YRruXpDuD75vqHUK7SkMwrWnmqZWv2wfPuFql7PsWtekmmYR7u8L038S6hHXnvpG7hIcDjoXajEP+9O1pzr6uPk6pVK3u/Z55SKejPW5Jhnx85y/1qZaYLu2anvs6w6N1DeEVrmQiFIYVghMSbWOTTJKdYrbF6f6ptc0lOpDuldKOpF0vP55Ih6+AQAAAAAAADspPaTrvf+9sSYCAAAAAAAA3K9K6a4AAAAAAAAAhsdDOgAAAAAAAGDPeEgHAAAAAAAA7Jldk6619gGSvkirXIvP7b2/Zt3+il021nv/H3f5fSQ9J5OspfCa7hLBQnJijLZyyXJpDNccI9bCNArBLH57qfMQET/bu+hAxdEUdltMtxvAxCas+Z1caZ+EWCPXNyWQTVJKqWmPiaYmZqiS4rqahxljiAOxkEyWElTtvg+vb5IOOpMQW5H3Z3hPCseck46LNEZKcr0qRnt5A4xbSikt3Y/S9ip9w/FSCcOr3Cvvkni21Rg7KtU0kk/rDKO4GefkUv8H3dQ1aQw7xABJqfmYLdyIRzonr/ilLF6rY+1R2NG+9tj+HiXle92ucr1kEk1L98rafHeua1JNEy88m03L0HlpaqBlqgVSPVjYn5V9P0Qdvavy54SL/eg2mktzTbTpp6nvZlP5dVRqD9ce0u3z9uzJatl76BPaxn3TUnDEN0p6a60O/VdJert1++fvOCMe0gEAAAAAAAB3SA/pHtHjz6aP7/iz8z6zvizPiO85rbWH72zay0QAAAAAAABw22vbqa9R994f2mWw9JDupZL+oaQ3SHrlqfZP3GVjAAAAAAAAADbZh3S991+Q9DGm/V+MPiNsuPNJbGvtmZJev6fpAAAAAAAAQHp+7/0NQw2WvkmHe0p7wkqbvbLqplvDPfx+ynfoUzNImsPUtKfpTgrzSGPYvrXXN5aL3l60619Er6ynHtcC3X1nLM0Yrq3avuw+6MD1XaS+4Zg7MUEHceFld7JW3zszjbTmamWhZrcQsuT3RwqOqO377beXFm9280hzW4TjxbWnY8uJC0uHMYY4Ty4z9/LiOt2VXVG6P6TtuT9ISQ5DbM90Lc1NPiwh3sdN0MHC95XJZ+mpb9hgJYijVNO4GiPMw73mVV83B9817/vtx6gcF6XMh8LchtneOH3vRelaHWuPwo72tcf29yipdv9zc4u1Tmo318plXK3d1Tr+42oKyaqEWaX6xUmv29VysZ4w+9P9/mpu6X01tU485ip9d6+jK9x5Us4ausDPMGO6NJ9B3X0jpUEU7nMtfKa3Y6d7s8txcM8gJGnhrwHNvb4QdmmDqJ7wvGX7a8d5jDs6AAAAAAAAgLviIR0AAAAAAACwZzykAwAAAAAAAPaMh3QAAAAAAADAnvGQDgAAAAAAANgz0l2vgD6bqs98MtBb+lRSt1JSSgpbcYktIZllacYOYUlahjFcIlsIRrJjp+3Z1zHAGKlvJTVtiOS1QVRSkGxqUy0larE0CWSmTZIWk832eTgw5iH568S8WZOlv0ymFNaKhUkGSulhs8l8sy3EL56EJDQ352khxbXKpZul9DeXbpYSzyrtt5Yz2/fIvK/H4b1O7fPl5vbm4fg8XmzfNybymaSpZRhj6VKpckRlaPfNOxvp2ld+eW6MeG3fbDNv/2pqbt+nMUIV5sLG0pma/2+r+ZNFSmnbbO8h0c0ms4Z0tBQKt3MgXzkV3iW6pfdp+yTY1L482H4M27cwN8kfi5VjeYhzZ5A6ZayaJh1whRdor6ny1+DF0p+t6dru7gXuniFJB+a+P2++73E4AafpxNxWqtvDPnK1yizM2fVN9da0kOJakVPod09sPTEX/Wqt4+oX15baXa0rSXPTnuroVKtW6nZX+8cU+9IH2e27DmLEa5//DJruD+Y+Xv7MWxjDHRoxwTow1wybwCpp4kaJyfJ+0q58aSEJ1iXZP+HSGT/gD4Nv0gEAAAAAAAB7xkM6AAAAAAAAYM92+uuurbVnSHpvSS/Q6oHf1/fe3zzExAAAAAAAAID7xbke0rXW3lbSyyV9kqTrp/7oeyW9+VS/JulFkl4i6Wt774+ef6oAAAAAAADA1VR+SNdaey9J/07S0/XEZQ7TMo3fJemtJf2WpFdXtwcAAAAAAABcdaWHdK21hyR9u6RnSJpLepWk71m3bei999bad0j6byR9oHhIN4rF4VSLw7PTXSObeOa7xiQYM0ZMgjnYPt21koSWEvLcGLnv9u3DpMkW2ivJQWMmvppH8T0lnlXS0VLi2WJzZ5yYFFdJmpgktMkAEU8pzcvNeT6tJXT5xDMfVTQzia2pb3rd0wESaStc+tcyHKCubzXxzCXkHYV0V5dullJcby78GC5971boe2LiF09Cet88tLuEtHROuQQr13ZWe8WuYYFxkJDyVUpmTddrs+/60r8Qex+IadVpZ7gUs627qs23T2OTpIlJcm0h3bW51x2DYAvRs0McXHYSYcelY8D1LyTWx3qkkkIfag83Rq5//Bg7p97H11EYo1Sn7H5cDJEcXLkmpr7uGpzSLNO1/cS8mEl4gekeY+dWKAhz8qxJJg912MHEXwgqtY6rUy5zTSP5/ZzTXXdPsne1ytFi+3TX4xBN7o6to7kfN9UvJ+bYd7W8NMznBPsZZMx01x2TrYf4TFipdVIAcrye2fbt67DwqCC2L00Kazqt+9zMLNyjbJ0iX+v0cL2+W/2ySBsfSPWbdJ8k6e20Ksc+tPf+76Qclbv20+uf71SeHQAAAAAAAHAfqKa7vlSr59P//PYDui38/vrnc4vbAgAAAAAAAO4L1Yd0t78NV/lrq3+w/slDOgAAAAAAAMCoPqR7xvrnH5zZ64metf55o7gtAAAAAAAA4L5QXZPuDZKeI+ntJf38lr/zZ9Y//6i4LWxpcX1y9+CIGPpgGgsLda/6m75pezaEIS2G7TfnF1nevm9ceDku6mz6pjmbMSoLL0uyj87Tvo/vya7SWplmgyk4wrUvw+KxYT10NbNwcpv7HeoWcU8LzabFlF3wwOHEL258bEIi3AK9kjQLYxyY1VFz3832tLD0NKzi7vqnMSrSfq70dQskp9AOF/og+QWuU1+38LJ7/1NfyS+onBZkdiETOTjCz8MtRr5ICzKbvulcVaG9VS466dgaIAzCDRGvtWlxeHOa5GPZLDYcr8vbBzy0Rehr2ts0LYQcxpib7YVgDHfNj2uymx1avoykeThpxWkj1ynbB1xVahoXZJXGzuETrq+fWgqqcJeoOIaZcyU4K7bHIA7TVgmZkEoHmLtGpcXT07XP1i/h/rAwt+xJCHOZp+AWswB5esXuGhWDKtxK6/L3uoOJD6Q4nGyOMWu+70GoX2YuDCJcYFygRLVO2bWuibVjaF+Ywj0GopmTx4VJrPqGGsj0z7XH9nWKGzcFR7iaRpJOTE2Saqu5qefT54TKZ418c/bNu6p8RjM5Cau+lTCIGHjjPgelOiyM4e6VKQxix1ALydcvMeDqYHOQSk0jhVtJCplwzae2twgBOkOpjv4f1z9fvE3n1tpE0t/Ualf9QnFbAAAAAAAAwH2h+pDu1Vo9e/1vW2vP3qL/p0h6z/W/f19xWwAAAAAAAMB9ofqQ7hskvU7SUyV9Z2vtbV2n1tqstfYlkr563fRmSV9/7lkCAAAAAAAAV1hpTbre+2OttU+U9D2S3lvSf2ytfdepLp/XWnuOpPeQ9DStvnXXJX1a7/3NA80ZAAAAAAAAuFLKK9713v8vSS+V9CZJD0r6aD2+PN/LJH2AHn9AJ0mf03v/l7tPFQAAAAAAALiaqumukqTe+3e11t5V0mdL+luSnnJnF0n/t6RX9N5/ZKcZ4q4WD/7/7N17kGzrXd735+2enr3P0TkSIE4kcdEFEZzYAXGTMYaADCQxBhcGoQCxXWA7NmCEicuhTByXDa4gxw6JwSZQuFwxMQVRjABzLUQkJVGMMSBFBgMBBZCCuAgdiYuO9j4z05c3f3RvndHu55m9ftOrZ2bP/n6qTo30nrfftXr1Wu/6zTo97zPV4sY90l0Dm0JXTRKtpLtWkmBiu0mvKSSQVVPTKomtlYS1yvt2qYCS/GdVTXy1STehrwvjSolLJp1wlWKNQpKPm6JcepEkrVbbO5cSz+ahfTYxaavTkLa62E43i31DiplLQkupZG6M1Delpk1HSHKtWBaSQF0ya0rhjalp5nN1SWqpb0xYS+eRSTfLfc32QjpaSmwtJaGZ68+1rQcpTBojnEJpCDc9lJLCUshXmq/dPqRDYe47KYEs/efPlTnOIQxRK9Oe0lZzu9nnkITmxojThUuCTQlyFywl69o6pZJGWqhpymOY8zPVNPHcsrVHIU02pd4X6pdSrZOOm2+uqQxSqF/SubUyb2YeNpcuE3evXIZ9c3XNfBKSPcOJ5JJcY7K8qYtc2nzqK4Vk+fBBpYRY56rUNC6JPqe7FvoWUl9jWnyhr6uNF8UUeluHmdpFkpbu94RQQ5Xql5hMPrDtLIXfu9zvbuk+EOdEd3/wXa2W7js++NkmU8caynwmsaYp1C8psd73HSPdNfU9++RYHp/v2ctQ53pIJ0m999+U9FWttb+q9Z+3PlfSQ5LeKelNvfd3jrOLAAAAAAAAwPV27od0d/TeV5L+780/AAAAAAAAAIrKa9IBAAAAAAAAGBcP6QAAAAAAAIBL9j5/7tpa+ypJXzXS2G0z/k1JNyTd6L3fGGlsnHLy8ET95jmft1bWC68syJwWwCz0jYtolhYXN4v/VhYiT9uL4ROFfYsLNZuFKtMxKiyGHZkXpLV47RqaYZHQblaBX4XlTns4oN0sCOoWX5ek5XR77HlY2H0aFlmeugWSzbipb1rEOAU8uDFyGMTwhZdbcT8uUlog2QWCVIMj3ILaKcjBnsqF0AfJL3CcFpZ2YRApICKFQbjtuYV7pbDIclx42Tfb9jFOocKC/QqhOWlRfCflYrigipRt012QQziebhH/dX/TtxLkUFhMWZIPeEhj7LqYcvnGsycx+Gq3kKxYp5TCJyp9h4dMxDFGCNSq1Et5DLd4uu9bDr6yGzRt6RpJ23ML2odCrJvrIQVc9UL9spyG+5HpG2uPEOQwc2OEvm6PU0BEGmPX+uUq1C7S/uqXSk2T+sfaw4ZB+HFd37gP6Vw2Y+eaZnjAVQ9j2LomHItmJ/fapOO658Ap83r3O956FNtaCrhy+xYOW0vztTmeMZzK/I6VAyKG159xe4Wa5rxhEO87xtnnxjzM0WO5u5R8P0nP173L8NN7PfRsuxozLAAAAAAAAHDFpOCIIY+V+6bfWX27pB+TdFTcLwAAAAAAAOCB8T7f0+u9f13vfXLWP5I+RtJvav1w7p9I+ihJz5D0kKSHJX2wpM+T9OObPs+W9BW998+9uLcFAAAAAAAA3D9Kf0zbWnu6pO+X9BxJX957/0u995/rvT/Rez/uvR/13n+r9/4vJH2KpG+V9CJJr2wt/VU0AAAAAAAA8GCrPjj785KeK+l1vfdvO6tjX6+k+pWSfl7SJ0nim3QAAAAAAACAkdakSz5f63Xm/umQzr33VWvtH0r6NklfLOl7itvDAIunNemh80VilUJtRkhNq6SRVtK/KqlpldRYqZb0ZtPN4vZSauFu+xajWStKSWjhA1xsN/WUtBkSDt0xSmmWLgktJYK1cOxtqljo69qrqaquPY3hE89s1yiNvS8p4c5xp0BMyAvtLvEsvWObxhZSolyq6noM1zel+pnEs9g3pWCZ95f6uuskprhWktDCGK69cB+QwnyWkrTttRPeR0o3c4mt6b5jqqWYQFZorySTVVJcY/89JZ5dkfDFXNPsmmS/xzpl55omtI+yvVKabLqPm7bq5+TaC3ORncsk9XRN2c7pfuQSDv3OrcJNu7n071SnuFT4WNPYZpv6mscwc60fNtdchVqnUtfc3aW4WAAAIABJREFUbzVNGqNS06T2OIa5x8Taw42b6qJ0PZjaI41h22MKfaG9klhfZc45e8+Qil+FSjdndzzDCC6xPszh6dKxdUMhsXWMOiVe1rum0EulBx/3ml4W0/2m2le/SffCzc83F17zs5uf/0FxWwAAAAAAAMADofqQ7v03P59TeM0zz/EaAAAAAAAA4IFRfUj39s3PFxde8/Gbn+8ubgsAAAAAAAB4IFQf0v241ksRfFlr7Z7fjGutvZ+kL9H6L4Z/rrx3AAAAAAAAwAOg+pDuOzY/nynpla21D0wdW2sPb/o/f9P03eW9AwAAAAAAAB4ApXTX3vuPttZeK+nTJX2ypF9prf19Sf+7pF+UdCTpwyV9iqS/IelZWn+L7i2S/qcR9xunnDzStHr4nAkjI6S7OpXkrlLCbHUMm37qu9aS0EKKmRu7mNK2cwJuNfHTpmqGvi5FKSZYmcaYRFlJhQsbNH2XxXRXd+xiepjt67umiCB7elbSykZINov7PIJKYKQfoJbuWtqcTUcLXQsJnHmfTWNl3NQ/jWFf75srCaP1CdsMkVIEzRURP9PCvBzTxmxia0q6NY3peMbPr9B3jGmgNEYhybdiX+GLY8xblZThCx4jvr6S9F5KrB9+n1v3L/Qt1UUjnDBugyFqM6a+FlIEXRplfB+F2nEV01ZNXGOhTpF8DVStXwqb83XNFa5fdq5dpNIkUK5pCvWLGySm1Lq+sT5Iaat76jvC7yU++Tn0TSrJ3QNff9YYzc1dqfYo1DRjpMKXjueeapr8+RXr63P03Xe6a+kh3cYXSvoJrR/GPSLp75zRt0l6QtJLe+8n59gWAAAAAAAAcO1V/9xVvfd3SfrDkr5P64dwZ/3zS5I+uff+M2PtMAAAAAAAAHDdnOebdOq9/56kl7bWXizpz0j6NEkfIulhSe+U9NOSXiXpf+m9L0faVwAAAAAAAOBaOtdDujt67z+t9QM5AAAAAAAAAOe000M6XA3zR6T+tBEHHGEx1xHWFh8l1GKMBaDtIqGFYSvjSmGfi4sCl7iFOMNCo3Zh2hSK4MaI+5tWNzb7MMKxqJwXeU3Syiqow7uO8pnGsfe1invBKJNDGntffUfY58pC5EFpcf/K9grvr3IKpWHT1uxi62ll8Mp+VBYhjkk4w7cXF2Te2/k5vOvepoB9XtcVe3qDYwRcjdG3FJJV2F517e2d97l4H68cfxuGVQnYUVg8vRRQFmqoEeqGUW5H91v9chVqF2l/81z17V1k/TLCvqWAllHuc4WAK3sajXFqFQJhUmhHOsVtoETp8i3WNLvWE/s8l4fuQ9U5r5FFedG4mj0PDwAAAAAAAOBeSt+ka619SqW7pJmkm5JuSLrRe/+uyvYAAAAAAACAB0H1z13/D53/i4ldEg/pAAAAAAAAgLucZ026K7LwCAAAAAAAAHA9VB/S/c8D+hxK+oOSPmrz/79J0o9IendxWwAAAAAAAMADofSQrvf+54b2ba19gqT/VdLLJb2m9/5TxX3DQItHV+qPuCgq4JxK6YQX++Xaxpd5gfvDKOmEITWtvDPYLz4RXBE2XbI6yH7qDOoX4D6xt3Rlapqr53xHfznZ77OXvaW79t5/UtKnSXpS0ne21h7b17YAAAAAAACA+9neHtJJUu/9V7X+c9dnSPryfW7rOmutPXH6H0lvuex9AgAAAAAAeMC95a7nNTvZ60O6jddtfn7WBWwLAAAAAAAAuO+cJ921arn5+WEXsK1rqff+6On/v/nT4Xdc0u4AAAAAAABAekHv/fGxBruIh3SfcIHbeiCtHlmoPboYb8DdV8UcZwXMyhj73Ge7CHFhe4VxY3vYXmkx5NL7K4xR2bfqeeH2rdC3Mq400vHcsW99geuCMa6TXe3xDZbe3gh9d95e8eMo5baY9ri/lc+kciwK+5bHSBdrGLvS1x6jPW6vYoTPZOe+FXuduAr2NceNcO/KfS+45qrUNNX6ZWjfak1TqCfGqT2292+UfSvtQ7HdbW9fn191e7u6CrWLtLc3WH57e5rzS/fxwvbiUdtXHVa4j+e6qNhujFLT7PiZlGuaXS+16jWy8/Z2fL107ut6pRGfvRh7/XPX1tpHSPrrWl+fv7zPbQEAAAAAAAD3q1G/3dZaa5Iek/Shkl4q6S9KeqbWD+n+2ZjbAgAAAAAAAK6L0kO61try3r2sN0v61nO+FgAAAAAAALjWqn/u2s7xz09K+o977/v9w10AAAAAAADgPlX9c9df072XEL0l6V2S3iTph3vv/9t5dgwAAAAAAAB4UJQe0vXen7+n/cAODh+Za/royble2/eUQNZDX7u51LeQSpW2V0kV66vhY1T6qtI3tLdV6GraW9peYENtwhh2PwoJa+l9xDHs9kbYt8J+jJGOVhljbwmz+xwjsR9V4fyspqrumrqV+qbvnBdSt+z2wrhxOpuYD6swRot9wwZdEprbh9C3/PlV3p/rW02vtZ/J7n3T4XT9W2Hf4rhhjFJwnp2wd58c4j7vqFS7xEFS7VGodQpjj1LThHuzHXtfNY3kj12835q+K9857YdtLdVQtffn66IR+qb3t+v2inVKrIGGbi/ZZ+JuYXs7i1PA7vWLU61pdk2Wr2yvVP+k/uEk8n0L46b+4cbjxmjhPl6rEcOJaLZXqqFS+wh1StqePR6lOqW4vUICrv+oq2myw/vfq35Z9nlt20V7TXcFAAAAAAAAcG88pAMAAAAAAAAuWekhXWvtda2117bWPqTwmhe11r6xtfby+u4BAAAAAAAA1181OOIlWq8A8HDhNc+U9Fck/bKkby5uDwAAAAAAALj2qg/pzuOJzc/nXMC2HkiPPu1Is0fO95fLlQWSwzq/dow07sr1Tdtb+ffkFkhOfVdmkd60MLHbt9Q/jdGXboHk2gLQMmP0uIjm8MWpxwgvKC1uvDRtadFk03fdf1hbtW9lQebKGPG4VQIz4rF3q32HcYMR1oEv2X1x4+GBBnmMQt8YUrCfMcrbm5jrPY0x3f6w47jTMEZhQebSyVhZZNm8j/V+mPbUN7zvNt2+MNMi0i50Iy1Y3Cb+gp+YsV1bGrvSV5Kmpj3uc2FyMKdhffHmPcm1x25jpHGXhfAJV49U+8YxTA0UA7zM6ZlrmnDxuAPq6h+FEJs4iYcPKow9eIOVEAal+iX13d5eufbYV61TqEn2Gqhlx0i/VIQxjPutplmPYS+I2hjuXhkDi4a9Po5RrFN8rTM8yKGHeiTWDa5/qnVsSE/om4aoHCNXk6TPKdYvLlAr3MfNGLmm8e0TVzekvram8RNJKh1d/3gqj1DThFOxNMYd8+XR8MHO4SLWpPuczc/jC9gWAAAAAAAAcN8585t0rbVPl/TB5l99bmvtt8546VTS+0v6VEmfrfWz6589704CAAAAAAAA19m9/tz10yR9zV1tTdIrCttoWj+k+weF1wAAAAAAAAAPjHv9ueurtX7IduefO1rhn8clfUXv/QdH3XMAAAAAAADgmrjXN+l+QtLXnfr/f1vrb8V9i6R3nvG6laS3S3qbpNf33m/vspMAAAAAAADAdXbmQ7re+1ynHtK11v725n/+o977m/e5YxjumQ/d0uHD54s4SommTkxNMzksi5jMOrxvSk1bmv7LkEy2WG7H/qT3vAwpZi5NbZX6TrfbVynxbFHIbUljmGSdmNAV83KMSqJpSGadLEziWUpxraS7pu259vA+bN/Qf7IMyUFu7GqCnOtfSdatJMGe0f9ClRLP/A6naauSoFpKZk3pp6W01cK4oX1VSjEzbzClpqUTw1UIhfS+FGCdE+uGp5i5JDSXbLZuD2ljpn9KJpuavtMwbhrjwLRPw/ub2r5h32KKmdnnlDzr0mTDh22TZ69IumtMaXcJquFEdGO4uuOs7bn6JY1Rqmlikr0ZI9QpS1NPuNdL4ZqUT32N5aS53+YU+uEJsWlqsCntYddyTbJb/ZLqsFR72DFG2F4LscaVMWwyayF5Vhqhfrka08soia29cILuK4V1NU1Rm66v75q25/rnGsq0hXOoh31emZOjdLpUIzTN55fmSfu5phTXVL8cbB+QSejrao9JrFNC7WH6p/v7wXR7IqnUNJKvX1ytJPnaI/Wt1C/Jveqak5P9fgftItJdAQAAAAAAAJzhXn/u+j56j99JAAAAAAAAAHBOPHQDAAAAAAAALll8SNda+wOttX/3IncGAAAAAAAAeBDZP3dtrX2jpK+88797739t879/dYdt9d77C3d4PQAAAAAAAHAtpTXp/tzmZ5P0n0v6a5v//3ydL2OnnfN1GOCxm7d086EUV3k2l0C2CksPphSzSrrrwsT+LML2currdvvcpLhK0tyMkfouQhKa678M6TWLhUttCsfTtkrdfcE1JioOT5SKzBgxCc2mn/relbSyySJsz/Sv9G0hmXWUhLVdk2BTe4i9s0m3haTNdf+LnYZt4lnq666pmI6WxiikrbrLLKSHpTFWJmGrkti6OgjbS2lq5o4dt+fSVsOw8Ti7F1Tnl8L2bHuYa10yq0tBk85IYTXtB2kMkyDmks0k6fDAt7sUs1kaw0xSKR0tpZsdmMkr9x2e6DYxfVPC7EVLqfCurkk1jatJYk0TLkDXP6W7npi6KNUp6f2dLExtNQljmBR6W7tIWoX6xd/2Q0KsGzpFP48xv9jthc2l9p1rD9+3VusUklkrqfcK9UShfimnyVbSXeONyo17RWoal/hZSWwN41bql5TY6vqmdFA3nbVY04Ta39Y6aQzT13eNJubYrdL9yNV9MZa6kNiahnD1SyHFVZImLt013cdN31T/pNrjwPRPfWdmP1LfabiwXU3i6pF13+E1TUp3tUn2YXv3qmuOjo7O/Pe7Sn/u+qN66pT7sbv+XTvHPwAAAAAAAACC9E26L5L0TyU93nt/451G0l0BAAAAAACA8dmHdL33ldbfpgMAAAAAAACwZ3wzDgAAAAAAALhk6c9drdba39r8z2/uvf/OwNc8R9InSPqd3vvri/uHAT7w8D16+Mb8zD4x9MG0L8Oz2zyGCWeIYRDbq4QehxVFU6DEiVlQ+WTqx3CLLx9PfN/jsMjy1Cz8eZJWUjUWYaFghUVXV2bh3fiH5m7x5eIqkDsv6FtZeLkQ2CBJE3NaVxZTTos058WbzYKicZ8LfdNiym4B6NDXHeccHFFYvDltr8IENiSV0IcYaBC2112QQ+zrxvXbM9OWJKmZ4IcU5OCmufRZp0AJd11WwiDSdZY+P7co9xjrdPe4qLNpD59JM4sbT0LIhAuIkKTZbPuAxCAHEwaRFiy+ceAnmJmZIA5jcMT2GC5MYr0fvv2GaXf7IPlFltO4U3PWpZCJpNI/1R6Vvq6ucfWI5MOw5qHvcWh3Y7uAiHX79uTg6pyz9sMtcL0I570LmWihpkmVpVtsfRW2190+h2vVhtVIapXCxmUQxZqmEHxVqXXCgcs1iQnCKYRP5L7Dw7PGqFNy0JbpHz8Tcyyq951d65pUN8TQI1MLxDABM26hTkn9XZDVegzXN4xr+rbwpGB1ULiGY/DgtlROxrA9czzTterql3RuxTNox+CIFuuU0O6CFULIxMzUKSnIKgVf3TDtN6ahpjF9Xe0i5VrHhUTcCGO4mmQWPuwUBuFqD1fTpL6n224fHtvXjaX0kE7S12p93r5S0qCHdJI+RtL3SnqTpI8rbg8AAAAAAAC49i7iz11/a/PzBRewLQAAAAAAAOC+cxEP6T5m8/PwArYFAAAAAAAA3HfO/HPX1tq/I+mm+Vcf1Fo7OuOlU0nvL+lT9dSfyP7COfcRAAAAAAAAuNbutSbd10n6S3e1NUmvLWyjaf2Q7lWF1wAAAAAAAAAPjHs9pPsxSV9q2ov5kfrnkv674msw0Psd3NIjIU3uXpYmFsclm0nSPMQMuTFSX5fk+lBMTRuewnq0DClmJpaohSif1H5i2lNwlA1KDalGPQ3ikn9MctRexdRQk/yUEs/c2yiklaX+sa9prySppfZSumsheVYKKWZm3Ni3kAQraZxozgp33oZFFrpLR0uJZ2kMk0zm2tLYPaSVuRRXyZ/7Ydqy1076+NJcNCkcTxv8nEJj0/niEoVTKpxrjJGKYXtun0MSWjPvO6W4TkO7S2dNSWgumSyluKYktJvT7cjHlNjqxkipaSkJzbXPwmTkUl9d4qskTc2JX+k7Bld3SLl+KdU6piap1DSp3aW4StLxcvvYn4S0+aPlzLa7FLrjxfBsuFSnrNL8afr3cIzcNezm+3XnsIN2LgnzsmsszHGSn9tj/VJKsg/z2Y71S+wbtmfTZFNtValTYv0ybFzJp4pf6ZpGqiW2utqjUKes+5vPJNQp7rJM93E7RaWY2vR7kAvyTYfNtIfAZXvPl/y5leaziz9fTFuqaUK9NDWf9UEhsfUw1CkpbdXVHjdDTXNo+/poa5c2v27fHiMn1g+vaVLt4fqft345PLjcdNfXSnq9nroUP3Xzv98g6fYZr1tJerukt0l6Ve/9DTvuJwAAAAAAAHBtnfmQrvf+bkkvufP/W3vvI8U/23t/8x73CwAAAAAAAHhgXES6KwAAAAAAAIAzDF+wYu2PbX7+2tg7AgAAAAAAADyoSg/peu//5752BOf39OmTejQsAHnHsvClybRA8ios1Oz6pzFmbXvR43lYIHkSwiAmJgxiGhdaH75IaAx4MO3LsBD5arV9jJZmcXJJamFBWNee1kD16Qy7h0zE9d7d5ioLLxf6Sn7R49LizYWAiHX78L7NhkwM7yvtHhyRkgfswsuhf+ESKbPnbWUx5dBXaZHlpQlnSH0Ptq/V8prCPi3B961cq2G6Xrr2uKB2aDdGOQfGGMMeo+GT0SQsyByDI0z7NMzXM3OPdWELUl442S2+nEIfHjJj3JikBZnHCI4oLMhsPuxJmMSn8QTdTappUp2yNNdlqlPmpsbINU04nqb9OBwjt2h1pXaRpJULpglz39L0dddC6iuFWqdUOKT3V6hfRphzSoc59XWHLtY0w9tz6INrC33TGPPhtVVbmKCYSp0iSa6WS/VLqdbxm9tVrLlT/WLu5anGHyU4woVEhHOuV8IgSoVDmBtcvlW4N6/M+6v+nmB/1yhMReVTyL3tQkhWDsMK9wfTnn7nnZl5fBbGTaFVrk65eeBrD1enVAOu3BjpvnrT1ECVe7Dka5XU1zld00ymJ4Nfdx4X8ueube0ZF7EtAAAAAAAA4H5T/XPX92qtfbKkL5b0YZI+WNKh6TaT9HRJj2j9sPrc2wMAAAAAAACuq/JDs9ba8yR9t6SPO928+Xn39y93/7s7AAAAAAAA4JorPaRrrT0k6TVaf3vu9AO4Ow/n7n4o1yU9LumNkn79nPsIAAAAAAAAXGvVNem+XNILN//7NVqnvX6snno4d6P3PpH0kKTPlfRbm//93/Tev3T33QUAAAAAAACun+qfu37e5udPS/rs3vtJa+3GqX//dEnv6r0fS/r+1tovSnqDpFe21j6y9/7u3XcZd3tocqKHT6XFrArPXpcmCS2lnKR0s0k3yU+mLZmGvi6tLEmJbgcmxeUgJNpMQ8qs27+UAOhSeCYTv28p2NG2x+SgC/6L8l2T0NLr09srJcS69LDh48YxUjqaSyuLCWthgza9NiQ/mTFistkqbK8QtBgTYg2XbCZJzaaRhr5mjJh4FsbQdHuDaTZc2YOx+7UaE+tM8qw9Puud82PYczwcz0LiWfqoR0lCswPv3tdOk+ENxnBs0z+lprl7SUpHO0gJsab/DZOkJvkk10qKq7R7Etqs+XFd8uhh6JtMCpNRpaY5sVGGvp5wKbWSNDXpe6lOGcPKnKGuTZIWodZx5+LSJLBK/hwf49qJl/UI1/uu4r1yb0n2KWG9MEZKW63UOuGydGO7FFepVuu4OmW9f8NrHXtDion1vnnX+iWehmEqaqbOj7/CmDoldW6hBnJzYiwnbMKo7+vqiUq9XB1jjN8Tdu47hjhRurk2dA3tLhnX/a4p+d9N0++rqd39juzqn9S+zzrF1UXVdFfXnlLop/HEXZuH9zSW6jfpPlzrU/87eu8nkrR5IHfn4dv7ne7ce/8lSd8q6UMl/YXddhUAAAAAAAC4nqoP6e48hHvTXe2/vfn5IvOa12x+fkFxWwAAAAAAAMADofqQ7k74w91/T/C2zc9PNK95fPPzw4vbAgAAAAAAAB4I1Yd0b9n8fOFd7a/T+i+yv6S19sy7/t1Hb35e8OJZAAAAAAAAwP2h+pDu57R+2Hb3N+Zetfn5AZL+cWvtEUlqrX2kpK/Teh27f7vDfgIAAAAAAADXVjXd9dskfZWkP91a+5re++9IUu/9za2175H0Ukl/StK7Wmu/K+kxrR/qdUmvHm+3cdq09fdJ65rKpHylmCHXHEOUho+R0lZdqopLmJV8glxqn4QEFtdeGXefUprahdt1N8ZIXBrjULg0tkJymxQSqGKCqkltqiaQFcaw+1xMcbVjF1LQkvg1aROuFEKivRSslGKw7PsLqWn2UBSOfWhvhVTqSrLgWe3XRuHQVebP2r0knAOV7YUP1rXHvoV9Swlk/v4XkmdNDGTanktyTWmtKS23olLTpJTZE1PuzuT7Ls2JOAknZz72wz+/0nlRuFmmc7a0b4XPr1TTXPO/rSnP4ZVDZ+876d41PEG1tG9x3DCGSwKt7FtKgk0KdY09RqHGiMnyZv5zia/rQcz2Qq3jkmfXY2+PEQOozVwZE1sLdfS+6vbyuBf5O8weXfTvhFfld2x3T0vJrC6RPfcNNYlpj/dx2/epfThIEdwjKT2k673/YmvteyXdkvSopN859a+/XNLHS3qepJmeekAnrb+B9w077y0AAAAAAABwDVW/Safe++eH9ne21l4s6b+X9J9KuqH1w7xXSvrq3vt8lx19kLXWnri76VJ2BAAAAAAAAHe8pZ36SmTv/dFdBis/pDtL7/2dkr64tfbnJT1T0jt7j1/ABQAAAAAAAKCRH9Ld0XtfSnrHPsZ+EN39JLa19pg4vgAAAAAAAJfpBb33x8cabC8P6XCxlr3lYIiNVQjydaENy2LorwuJcAsvp+3FQIq4ve3+KajCtbvXX4Z+RfZjZ5W3sc+3bMZOhzi2m9MoZRRosv0v0veG3SK/kiSzCHE8Lcyiqy3tXNieW/R4jI8kLW5c4sYwx3jdd/fNDd4H5ffn2kuXdfX8vCZTRrSndZPHmPMr83W6p7n22Nfe54bfV9f93YTmFzq29/0wodm+6fDs6Y8pKjWN5I9dpdZJNUZle+nzK50XhclvjBpjb/XSFVmsfV/Kc3glb8jWOukeFTZn20f4UAr1Ug4SM6EIxRpj17M2bi9NGYVawI5d3J4do1AvhWkrnFu+7/7qsD33v6Kuyu+EF/07trun5UDJ7Xoi7VsKg/C1Ttg5G6751Lj3evayq62HdK215+5rY733X9vX2AAAAAAAAMD9yn2T7q3az3/j6mF7AAAAAAAAwAMtPTS7Gt+5BAAAAAAAAB4A7iHd63XtV4sAAAAAAAAAro6th3S995dcwn4AAAAAAAAADyzWiLsGjvtMR/3wzD4xFc6knMQUs/BX0PM+3W5b+VPruG+3H61mvm8YY2G2twj77NoXq+3Xr9vD+zbtrk3y6TwpsccEWJ3Z7jsX+ib7Smd1fWPyl292H2tOpTLHPqWcxu25MXxfdxpNVn57KYlwIpM+5MMXpWbGCNtL7W6IXjrhgkLyWp8WElRDWpn7nNLYaXurg+GftZlyYv+4b4VzOV4nlWvqKqS0JSPMcaW5Ngzt5vF0H5iZ8yjed8IY87Z9Ih2EBLID03cSrtXU7qaXZGUmB5ekJvm6YVpMcZ2YtOqkkhaXEltdnZJqHdc31Smub2qPfc3NxLVJ+dxy52JKgnVjxJrGttauvzGu972ppLDuaV6WirWOrVNCrZPuXeZydfdEKZVnYeeW4QM013sPc0AzB7+lcYOd65piWq6rVXItYNpTnTIN16UrB1NtVdq3YW3SGYnClVrHtY9Rp+yrb1I43aq/+1Xm2srvq/FeYmqPxcTf3+3v2GHSOU6/rzjpfHG7kaai8MFOzYFOtY6ra07XLsd7vl8Nz58HAAAAAAAAsBc8pAMAAAAAAAAuGQ/pAAAAAAAAgEvGQzoAAAAAAADgkhEccQ3cWt5UWz4VHJEWS3TcwsnLwmLKkl+cOC167MIg0rgpOOLEtD+59Is6nyy3xz4J+5b2eb7jQpyrVWVlYvkFSGNftzq871pR2twYoQ+VMQoLIVcWYpWklT14YTFl0+wWX5fyosd9st3frNm6bjeLrsbFlMMbL6zVPgq/+HY4uUohDMODI9Jiym57q1lavLkQPpHOT7tvqe/wMarXlO1bnKKsURZfHp524a7hNNeuYijQ9qSRFlM+XmzfdyZhsq2EIiQuLGExHR5SIEk3JovBfWdtO7FmGhZTdn2TtCDzrlLoQ+Led6XWGaNOOV6GQC3T9yj0deNK0tzUOu6clVJwRO3acddaDo6opNgUFIaoznGuuVanhHkrhUuZ5hY+ExeKkEJCIpfTNPddl+Z+Owm1R6pJbP0SjoXNoDm4wjWNZOuaeA4UAsriGDuGZIVpJIyb+g5vv/DfE8aoaZJ9/Y6Wfi8x80CqU6bmmnK/w677+g/Q1S+T5e738WV6f67WSeFbk+19nplASinXZ5Va5171y3vMfXdMfJMOAAAAAAAAuGQ8pAMAAAAAAAAuGQ/pAAAAAAAAgEu200O61tozW2uPFvo/ssv2AAAAAAAAgOuo9JCutfYnWmvf1lr75dbaiaR3SPq91tpxa+2NrbVvaq193BlDfENr7Vdaa1+7y04DAAAAAAAA18mgdNfW2qdI+h8kfczp5lP/eybpozf/vLy19lpJL++9v/muoT5a0gskvfjce4wtvzN/mo7mN87ssyw8j3VJK9IZ6a4umSX0dQmqxyFhZhHaXTprSjxz6WYp3dUdWsXSAAAgAElEQVQlwUrSiRtj4fsulyY1zbRJOTWtu0Sv/QTkZTHBykWQhfSpShJsCMjx5+LwOK9e2DdJMmGrpX1uSz9wSjFzIYmxrzsHqimuLmFtj+lo9jinZFZ7vgxPZk3947lVSUqNqWnDxk19QyiVTWNbj+32zY9RSViLt4ddw6Nj4llKJTZNKeHQpTmHvmkOnqco5R2lpMWFmWDifc4kuaYkUZfiKklPTrYnmFlIKzswfafhxmPT38JEksbYVappUv3i2tMY7jOZh4snfX7us8q1x/C+KYXe1jqxpjHvL/RN14671mwqZ+ibEw5T+/BYRpvRnl6+p2T5cBuP78MlEaYMZVfXTCZhg6kGcunm4d7lklx7sdZx00CsdVxzpe8I4ulWqF8qxz4n1vvd2LnWKfSNSbCxfnFtw99frMPi7yWuLV0PYYx9KdQpMYXeXH8tvL+T4XsW2Xtl2LdF3649UjLr4cSfMK5OuWHaJGlm2g9C2nyqSVy6a+rrnK5pbs8HPUY7t3s+uWmtfbWk12r9gK6d+mep9Tfp3qX1PfH0v/t0SW9srf3JU+NMJX3Upu8Pj/ouAAAAAAAAgPvYmQ/pWmt/RdLf2/Rrkt4m6W9K+jhJN3rvz+69PybpUOtvx/1NSb++6fs0Sa9qrf2xzXAfIenm5n//yMjvAwAAAAAAALhvxYd0rbVP1PpPXLukhaS/JemFvfdX9N7f1PtTX/DuvS9772/svb9C0odp/bBuofWfwX5Ha+2mnvpT2V/qvb91L+8GAAAAAAAAuA+d9ce0f1/rh3gnkj6/9/5DQwbsvS8kvaK19jOSvk/ScyT9ZUnP3nThW3QAAAAAAADAKfabdK21T5L0SVp/i+5vDH1Ad1rv/Ycl/dda/+nrX5f0qZvxeEgHAAAAAAAAnJK+SXcn8OFXJH3jDuN/o6Qvl/Q8SY9JekLS/7XDeHvVWntM69CLz5D0EknPkvSy3vuPjvmasb1z/ogeOrl5Zp+YeFaIukmpYj7ddXjyWkpmTQkxLmUmJqGZxLI07jwlAJoxUuLZwqSmxXTXkI7VXf+QBmRjzIppV+7USElolcQztx/VnD935FaFxNaUeBaTrcxnEkKGtLLJrL5vC5+fSzeLfQufdSukZybpvdhhh4dHx3Qtl1gXx43nZyGxtZBuV0n4jeeWS3dNSbAp0c2mpoW+lXTCdJ249z1COloL9yObDhlS/dw8GabUkp72rZB4Ng/ts8n2G5xNQ4qZSfw8MK+XpMOQ7ur6pxSzA3PBT8IkMHXprsUbTyVNLdUvtm84QZdmjFVKbDXtaR9SPeHqmlh7uDTZkLaazq2Fq3VGqGmWi5Q4uN1uaxfJ3/jjfW6MCWa7Kd8HwnnoEjFTcKmZg9OtK95Wzf2vpX1zc3uY/NK9xNUZLsVVklamb651htcecQx3oIt17VWoX3Lf4fsQU19tQmwYwyW2VlLhC0mwkq9TKnVRrGnGSLJ3x7465VR+73Ip2Ol3wnhBVE7QbXHXwr9wdc18EpLCzYmUUlynoX5xdUqqaVzd4GqX1FfyCbGJq2tOj/vkfDZ4rPNID+leovXn+l2n156r6r3PW2uvlPQ1m/Fe03ufn3e8fWmtPU3S35H0cq1DMCTpNyS9WdLvjvUaAAAAAAAAwEkP6T548/MnR9jGvzz1v6/cn7q21j5I6/16kaTfl/QNkr5b0s/07v+72XleAwAAAAAAACTpId1jm5+/OcI2fuvU//7xEcYbTWvtYUnfr/XDtp+X9Lm99/937NcAAAAAAAAAZ0kP6Z6Q9AGSHh1hG3f+YHch6ZdGGG9MXyvp4yW9TdIf7b2/e0+vAQAAAAAAAKL0kO7tWj+ke57e989Vz+N5m5+PX6U/Bd38yepXbv7vVwx52Hae11yEx48f0Y2jh977/yuLLDsxZKKwoHYlOCItvu0Wek790xhu4eQ4blo42bSv0sKfZj9yQET4nAqLLI8RHDF0WEl+weLCwstJXkzZjeu7unVN0zFuYd1Qt+jxKq3lasaOixWngAe3IHPsW1lMuRA+kbi+YyzMn8awC/qGzpVgk7SYshsjrXseAxfMosCxr2krLISc2scItSgFdIxwDsTz1rWHua/bazjM4WlzLkwgbM+1z8NC69OwyPKBCYmYFhZZdoENUg6UcGOkxZTdAskt9XXBEaXJZX8q9UsKCXHhE2ncWHu4OqVU04SgihAo4cZI53KtpgnXnwmUyDWNm9x91zHql8q8Fec+N38W9i2VbGkNeDefuRpD8gvrV2oayQdipXPOXdrVkCwbHBEDrobXRSVpjMI9rfTrVSXgqnp+VsInKnVDpaaJgRKmsVBblesi8ztIfH++uWbH4IgchpXuR654TOFb5nfsqb9Yl9MQLGT6x9rD9S3UI2ns1LcSHHER9cvx0eUER7xJ0h+S9FmSvnPHbfzxzc+37jjO2L5I0k1JPyfph1pr/5Gkz9P6oeJvSPo3kr6z9/57O74GAAAAAAAAOFN6SPfDkv6MpM9prX1o7/1t5xm8tfZsSS/T+hnza8+3i3vzWZufPyXphyT9CdPnv2qtfUnv/TU7vKaktfbEkG7nGRsAAAAAAACjeUtLX+Hb6L0PXkoufZH2hyQ9rvW3xr61tfQF7ay11iT9I0mPbJq+pzrGnn385ucXSnqxpL8o6dmbf16m9ZpzHyzp+1prH7LDawAAAAAAAIAz2Ydvvfdbkr5e629sfaakb2+tHQ4ddPNQ75slvVTrb9H9SO/9Z3ff3XFsElrvPMk8kvRHeu//pPf+25t/XiXpD2sdoPGIpL97ntecZ99674/e6x9JLzj3mwcAAAAAAMAYXjDgGc5gZ31D7lskvU7rB3V/WtIbW2ufca8BW2sv1jps4ss2TU9I+quVnboAzzj1v/9u7/1X7+7Qe3+7pP9x839fcs7XAAAAAAAAAPeU1qRT733RWnuZpNdrHSLxByW9urX2K5J+VNIbJP22pBNJz5T0UZL+Ez31J6FN61C1L+69//Le3sH5vEPr/T6U9P+c0e/fbH5+yKZ/6TWttfe7iBCJd956mmazh8/12koSbExCc6GThSTYlHgWt2fSa1yq6nqMYa+XpB7GcP1dgs66s+lbSTyTJNe/kA5ai5+Sj80yCU9JTFyyCTrh2Mekqe22lCrWXeJZSlstJJNVEltj4llMfnJ9hyeeJTunuI5l19S0lHhWSX2tJMGm/S2cn5VzuZSqqmq6ayHxLF7Dw9pSezoP0ynurr/KdBYvnfDGXSpcSuNemsTWSUggayHxemKSXGMCmRkjbs+2+sSytD3Xnm4D91hy5Z77MIZqir1N8o3nSyEJNtUvrm+lxijUP5KvSWJdZE6jVNPEWsddJ4W6qJICKvm5pJIUHi+SMPe5HcmJrcPrhrh4UKF+ce3VtNWlrT2Gb6+azmvrmvABXon6pVhG75r6GmuaQjr9KEmwrr3SV7V9s7VHTIINH3Zhn2179R7l0ofT7yVubJOMve4cUl/NjWoVbl4rNxeF5O5lqFOaS4VPfc1nneqUPIZJdy3UKTnF1TaPWr/Mb+83IiA+pJOk3vvvttb+iKRv1/pPVyXphZK+4oyX3dnjJyS9rPf+Y7vu5Nh678vNw8Z/X9Jzz+i62PxcSrolqfqaox13FQAAAAAAAA+AewZC9N5v9d5fpvXadD+l9UO4s/6ZS/o2SX/gKj6gO+UXNj8/9Yw+f2jz85d670fnfA0AAAAAAABwpjO/SXda7/3VWv+56/O1/rPWF0l6TOu12t4l6e2SflzSq3vvT4y+p+P7dq2/Hfiy1trX9t5/8fS/3ARF3FlX77t3eA0AAAAAAABwpsEP6e7ovb9V62/K3e9+WOu19f64pNe11l4q6V/33ntr7QMlfZek50h6q6S/t8NrAAAAAAAAgDOVH9JVtdZeqM2fgPbef2Df2xtq82DtP5P0Q5L+qKR/Jenx1trbtd7fidYBES/rvT953tdchN+/9ZCm0/MFR6SFtmuDDF/0uLS5GChh2iqLEKdxK2Psq29oj4ssp1CDikr2hOsbFl52CySvwsKhORjDtVU+6+Hjpv5x3dDCyVwZo5VWyh/e9Vz992HHBZal3RdpHm0MuwhxYXuFhZfX23MnTOq7+/b8AuwjnETpGjYnaAtv0N674iL3YUFmGx4yfHHqZVgFPq0BLvP5xUWMXRBH/JzCIsuurfL5FfrGfbtgpZpmhOCsuDlbF6W+w7eXznE7dhxjx76pfyXgKoY6jXAiFcKw4mfi3l4Mg3AfYOgb3p+91CpBVCPUOklle3mM4eeLdRVqF2mU+sUp3YPP6j9wjFKNXwy48mMMvweXahpVQ8eG9y1J16oZPB2L+KGYm2sPNU3l2OcgnMIvQq5OKYRMnDn2sM2dUUPtPmncq65Z3trvxHTPNelG8Kck/QtJ33sB2yrpvf+u1uvLfamkf6l1cutHSHqzpFdIenHv/ed3fQ0AAAAAAABwlr1/k+6UK/LfWd9X730h6R9v/tnbawAAAAAAAIDkIr5JBwAAAAAAAOAMPKQDAAAAAAAALhkP6QAAAAAAAIBLdpFr0mFPTp441FQ39r+hUZJgxxhjx+UNR0jEjAmcY6RVFVLaxkjdckOndB67G4VjkZJgS0m3lTe4z2TWHcc9V/+7jBG0eZXteqlL2lsamzRCUtgI00gpYW2EdK3S+xseCLbu7uLGQvxiW7qJKwycIrrsMdr9WOwtoLl6we96/VzJlYRHdNE1zSgTWhp7YFtQrml2rT0qNU117ELCYeLSdUOYc0jWjSPbVns0Rqg99la/jHDtUL8McNFpspVxx7jPjZC26uuwQvppcXs2MTnUGG5uiInLpeM8Qk1TmfJH+azHmDR2H2LMumb5xHy8wQy+SQcAAAAAAABcMh7SAQAAAAAAAJeMh3QAAAAAAADAJbNr0rXWnjviNj5gxLEAAAAAAACAaycFR7xV4yypCwAAAAAAAOAezkp3ve65XtfG9ImppjE2E7vaZyDblbBz2lFIgj3Xzuzgun9OuH6u+38KKyY1hrzHcfZloJb24rp/VnhgXJuaZoR05Wz7gqemAQa4LvfKMdKVC637Qk2zR+/Z77OX9JDu9bo+lxkAAAAAAABwpdmHdL33l1zwfuAMrbUn7m66lB0BAAAAAADAHW9p7anvKPbeH91lMNJdAQAAAAAAgEt21pp0uCLufhLbWntM0jsuaXcAAAAAAAAgvaD3/vhYg/GQ7ho4ePdEB8sL+FJk5Y9sC33Li/ya1S5LY6S+lUWIC2PExR0rY8S+ZvQRFlPu6XQqbM+2x30rrLpa+qyLS2uWPuvC2Hu6dvIY9+GSomOs9u3edmnx3+I+VLZXugCHt7c0bmnfUvv22PHUqlwOlTHC+2ur4fsQj/yO50v5WNjPb/j24rEYYYHrkvtwetl5Xq1ODfZeMvzA5Xv+8Pa91kWVvq6eqNzztft7KdVQaT9Kx744Ge1a1+xr3Nh3hO2NsR9Xxa71S/Ut73rfGOE+fiVqmtReqWlcLXEGe3rus4YqXO6VY1Grw9IGK+MOP1/2WtOMOL20W/t99rLTQ7rW2usGdv3t3vsX7bItAAAAAAAA4LqyD+laa8+T9P2SXinpm3rvT4bXv0T3fibZJPXW2n/Ze/+N8+4oAAAAAAAAcF2l7+l9pqSPkvTVkk4GjNPO+OeOzzj/bgIAAAAAAADXV3pI9xKtvyH3A7335YBxPl3SC8w/Hybp3276fMJOewoAAAAAAABcU2lNug/a/PzXA8f5jd77/+f+RWvtn0v6yFNjAgAAAAAAADglPaR71ubn20fYxi9sfj57hLFgHL5bOlic88UjpH/tLXUrfM+zN/OCEZJZU6KpTyYb3rdPwrKNI2xPk+1/EbdXSmEdvs+l91cYd91u+oftNdce3rPtK8meWmGfbXvsG/bD9K/13T2mKG1vDL2we70wkaS+bnuVvmkyWsW0MXP9haQwux+rsG+hXWbs3Hd4oluLY5gk7dC3cDnE9+1S1vIYw15/drv7TEJftx/F7VX61hLkdh9jb2lq+wpqLKfCD+9aqhtK99VwvbtLtRgaZ/vHGqrQN9wrfa3jx3Dnlq3jzhrD7EdhWo51Q+l9xzpseJ2StmdrklinmFognS/hwp7sWL/kuijsRmGCqfUd3LWkUrvkMVLtYeqGsL0x6he/vUINleqDtM/2vjq8Tumm7lj/i+F1UaWmiSdRJW21ULPF4OdUNxTeXyn1vlK/7LPGqNR9o6TMDh/jXn0ntwtjnUOa1u88UBvjId3jm5/PGWEsAAAAAAAA4NpJ36SbbX6mVNc7TrR+BnrWs8bf2/x81hl9AAAAAAAAgAdWekj3uKQPkfTMs17ce785YBvP2Pw8LuwXAAAAAAAA8MBIf+56509UnzfCNu78mevvjjAWAAAAAAAAcO2kb9L9rKSPlfQSSd++4zb+w83P39xxHASHT3TN5udb5XTX0Ic4RiWcoRzk4MISCmOEvqupb3f9U1+/vbBAaxhjNTWfZdqea4wr6YYFoCuLSLt9Swsk2/cRFhtO7WbsSdjeZLq9+mjsO/Erpk7NGNO08LIZI3zUmobtuQWS7YLOZ+yH3bdK2MUeVcIgXDhDCmyI68+a/suVP5mXbsHisL08xnb7Kizo69pXSz9ufN9Ls89hDLt4s3n9egzbbOfm+IkWPuu4cLJbkDnsm2tv4f2Vxkj75vqmhZfT9gphF36x6HBdVxaGriyGPUbfC1YJcqiFPoTtFcIZXO0iKQQyFbdnaoRS3/T+psPrF1sfhL6ppomL0bvGQpBD2rf4VYVK/WK210wtsW4fHuTgahrJ1zWp1km1h2tP9cHUjJ3GjbVHoW9qt+NekZomhksZy0qtU6hJUt7CyvR1+yBJS1NPuNev28OxsGMMD8lydc5654a3p4+jmTMxhXaUap1CCEOqDyaV+qWwvUmq78IYtn/l/cW+YW7YNeyiGlQxRvjExuRov/NPuj29ZvPzM1trQ/6k1WqtTSX9Sa0Pyb867zgAAAAAAADAdXbWQ7qFpA+U9Bd2GP8LJT1/879fv8M4AAAAAAAAwLVlH9L13t8h6Z9p/W3Pr2+tfVh14NbaB0n6hs3//Q1JP3LenQQAAAAAAACus/RNOkl6haS5pEcl/WBr7blDB22tPVvSD0h6ltZ/6voPeu+LXXYUAAAAAAAAuK7iQ7re+69K+i+0/jbdvyfpp1trf7a1uCq92toXSHqDpI/R+gHdT0j6h6PuNQAAAAAAAHCNpHRXSVLv/Vtbax8p6cu0Xp/u2yV9XWvtByX9lKR3bLo+S9KLJX2WpBds2prWf+b6Bb3H7DiM4Ma7V5odp4i4jUKKaynxTCFttZCEVklxlXxKVyXFLCWzTgopZi1cOXZ76SpLYWMu+SnG1wwfN2bQmPianhJbXftBSOw5MOlhIdnMpapKPt3swIwr+bSxg6mfemZhewdmjFkYw/YNMUquryRNzKdyEMZwiWcHIUapko42ifGSw63SBWj7hkQwczKnvouUNmb2YxH2bb7cvlhT35TuOjftblxJWpjEs/Q+Fgs/hkteWy79Z71abPft8b+xpc9veIqZFWN4w9ZcgurC7/PEfDe/kuI6xhiTRZj7KumuMenNvb6Qjpa2lz4/E3E3SrrrvsLPKmnzoX/uu2OCfGhPfV1NktLfS7XOQSWZNexbqnXMeWQTpZVSJ0MSrN/c2X/3c7dCYn1KbFWhfpmYvqmmSe0HB9sXfKobDswYqU6ZFdJdUz3hxs61x/DaKqa7mrOgUtOctR9DVWqadf/tky7VE7Zv+MXE1UXr/qb2CGPYvqlOsX1DnRLGcPXSwtQjkk+CXU5CYn36nbCZ/uFv91ySaz5VQtpq4f7n64bdU+hd7TJG33W7uf4K+1ZPvTe1R6WmiUmww38ZLk0vp/pO7vXsZUdnPqSTpN77X26t/Zqkr9f6jH2epJef8ZI7Z98bJX1O7/03d95LAAAAAAAA4Bob9J8Keu//raSPlfRqrR/CnfXPr2v9zbtP5AEdAAAAAAAAcG/3/CbdHb33n5H0mZsAic+W9CJJz9mM8XZJb9H6Id4bendfhAcAAAAAAADgDH5Id0fv/dckfcse9gUAAAAAAAB4INVWxgQAAAAAAAAwuvI36XD1HP7+XIeH8/f+/5zgZ1QSz1Kyjkt3LSS2xhSzkLbqk9CG902JZ6uQUtrNVdKWaQw3gO1qU4YkH3zYwrF3sWnFUKpaEpppz4lnJiEoJZDNQiqqSTE7NClo6/btCKPDsL0bUx935PofhmikQxN3FMcNY8xMVFFKd52aMyOlv7m+Ui0hbVpIR1sWUsyW4eRyfVO66zxMGgvTfrz0t7m52eeTEGV4ElLMjpaz7XFDwtrxYnvsNO5J+FxdItt8HhLWbFtITStNRrtHhccRCumnlRSzydy3N9M/JbbaJNi0vdTuUtNi8trwxLOYMuvSzWJqmtleSA6200gxxbWlc84NvWNNsx7DtKVEd3e/LabQuzojpaa7/UipqqtU65j+7h6c+qa6KJ1zq+2pLwY/u48v1Skt1YN2Limk3qd7X6hfXF0zTcnypiZJNU1KYXV1TaxfTK2TkuVvTv3kV6l1XJ2Sap2DMGHPTHtMdzXt1ZpmjPrFSTWJu7em9Hbf15/4rk6RfK2Sap0TM3al1onjpvrF1TqT4an3CjXNKkwOK3MvSUnTdsoo3F7SGDmNdHvwStq8FOqUUNO4MWJ9UKhfYl3k0l1D3ZATYl3tEfqWkmBTXbtbkv37OAkHcSR8kw4AAAAAAAC4ZDykAwAAAAAAAC4ZD+kAAAAAAACAS8ZDOgAAAAAAAOCSERxxDRw8caKD2fHZnQoLL6eFkGvBEaHvwfZz4Rj6EBYytossz4aPERdkTu1m7LggZUzdMOOGdvdRtbRotVtHtbiAd4lZTLmZcAdJmpjF79NiyikM4sZse1HOm2bRZMkvZPzQgV9d9WZof8gssuzaJL/I8sOTE9s3Lep8w6zcOouLN5vFqcPqqpVFllPfMbgFknMYxPYFmAIiKu1HboVzScdu4eUwCTy5PAzt22MfhUWWj6bbfZ9c+H2bTvwYbb7bxZ1evQqLLNu5Pa2wayauSkCEVFtk2bXHhZDTIsvmeJYWZA6fRyV8Io3hFlOezMN1nQIe3ALJYd/sYspLvz0fHFFMSNpVqmlCuw+OCDdWd96neiQFR9g6JQS3lGqaFAYxfAx77YSaJs3XrrnykbQU2pFqKxcYFkNCzBgpqCKEeUxMSIQLiJB8XeOCrCTppqlpJF/XxPrF1CQ3Q5BDql8emm7XKjdC7eHab4aJMtUklfplZibyGBCxp/qlEnCV+lfqlHmoPY5dcp18QNXtla9TXEiEq11S+1GoU1L9cmxq/1TTHM+HP4aYp7nIzQ0pOCKE/pUUgpPsXJtqmkL9kuqUaaWmifWLaxte61RqmjRGDK1abJ9b6d5VCo5IfZ1TNU2f3+PZy474Jh0AAAAAAABwyXhIBwAAAAAAAFwyHtIBAAAAAAAAl4yHdAAAAAAAAMAl4yEdAAAAAAAAcMlId70GJk8caWJSU99HJXYrprsOTyaLfV26a9j3VWjvJiFttUhjmMSzw5DOFNKAfKhUSggyndOxLyS2tpiGaNp818y9ICY4unTXkLpl0tEOQhLsLKSmucSzh2c+QfXhg+121yZJj4T2h0w6a0pHc0muD0990s/NFlLaTOxS6uvT0Xwc1DScMFNzfk72mO66Mif5SUg8q/RNia0uIfb26oYfw/Vd+r63wzngktBuLfwYBybxLCXWpfaKbpLQYoprSlpcFCIcC+L04uazNPeZKSOEE8ZE011T06Ynadzh7TGx1aWYhb4u8Uzy6aypr0yaWlulvqb9iqe7yiS59lCnyFwPqU5JCbGuJnO1i+RrnVjTpNRXMyWGgLxSCn1Kc564BNyYoGra0nU9xulSqGkq6fQHpqaRfJJrSnF9eBbSVk2Sa6pfnmaSWZ924GuPceoXl2Tv+6aaxNU6hyHm0tUklb5jcPWIJC1Dhb00J7+rR9btJoW++5om1Tquf6xfTOrrk5OQ2Drd7vueiU+NPZj4z2S68P0dl5abEnRX4fqz/VOKuZsHijWNnV4qv6OFvuFw2rrG1SOSr1+qdYpNiE1j2HTXUI9Uap1U07j2UKekhFjbP6W7uvrlVNtkceRfNxK+SQcAAAAAAABcMh7SAQAAAAAAAJeMh3QAAAAAAADAJeMhHQAAAAAAAHDJCI64BtqtJ9X8+ur35hZOjosp+/Zm+vcDv0PNtKcFmdssLDC/3D5t08Lg7dCMnVY3LixYnNZd7mYB0snED5wWa3fr48bP1+5IceVluwpq6GreS1rkfmoWZD6Y+pVRD0P7jalZkDkthGwWWX56WEw5LbL86HR7EVC3wLIkPWL6Pi0sppyDI7bHjn1NuwuTkHxARGofI6QgcQv6noT/NrQy12VaTPlkEgIl3GLKMTjChD6Ez/o9y5u2/Ym23T5GGERaOHmx2j5GS9MmSQuzuP8kzO1hCP+f8QprLMe3XGhPY7hFlsPlEAMlKgsyu8WX02LK05OwQLJpj8ERJ9tvps3DG0yLLM/NG1z6MVzIRFqQuRQckcbYVapTCsERrnaRfBhEm/o5p89CGW2DI0JddLjdnha9buliTXXNUOG4xXLJvJW08LldP784N5RuU+6tFEItJF+3TU2YhCTNTK1TqWkkX788WqhfHgmhD66mWfffrX5JwRGpfnHBD6nvxKysfxgCIvZVv6R7cAyOMO0uIELytUcOjvAhDLvWL7dTGISrG4q/U7hjF+sUW9OkcL8Q5mFCdlbp967C7zsl4XyxwRFhnsz1iwl1SjWNDacaXtNIvq6p1DQp9GFy4nfa1TUx4Grh+oYDV6lfUt971C9t+eSZ/35XfJMOAAAAAAAAuGQ8pAMAAAAAAAAuGQ/pAAAAAAAAgEvGQzoAAAAAAADgksTjYCwAACAASURBVPGQDgAAAAAAALhkpLteA/32bfWQZPNekxjXud0WUlwV0s1cQlo7SIlnZoyQjpbSzVzIU79x3njb0xsMSW/m2FUSz3oIh4kJP+5tVxLPxgi7KqS7ujZJmpgkNJeCtm4PSWgHLt3VJwQ9YtLRHpr6tKtK4tmjE9/30el2qk9KR3ta8/vh0s1uhtimG+aEOXTRUZJm4fNzp+0IV47CqWzz2NJlPTcn3VH36W/zcAHe6tuJZSlB7nbfTn2dhAS5aUyWG55c6RLPFuF9uMSz1D5f+k/QJRGmpOkYnbjH5F+7OZeElg6xTU0LiZgxTc0koYWkcJ+aNjzxbN2+vSOuTZKaSUKL6a4uxVUh9SwloS22x+guBU2SVmaMkBobrQrnVqpfnFCnyCRCt5RY7+qXkFgf37era0JS3MQl44ZkyMSnFqaaZrttFTbXlimF3uxzSGUsXdf7UpzjKumuB6Y9pbimdPqHTHuqX1yS6zMObtu+qX5x6axPD3XRw227782Jfx+p1pmZD/xmmJhnpoidFmoaafe6plLTSJIL0DwJ9/djU9cc2Qhk6SgcT1frzELt6OqXXWsXSVqFXxQWZjJZhPnsxPUN9c883I8qv5eEt1JTSKGv9I11innbKUnbJdbHdPuY+mrOl1TTuGTWY79zKd3V1S82mV7y9YupXSSpx3RXV7+Evu4X+FO1S1+R7goAAAAAAABcazykAwAAAAAAAC4ZD+kAAAAAAACAS8ZDOgAAAAAAAOCS8ZAOAAAAAAAAuGSku14Dq6NjrdqpyJpCElpzSWgtJHSF1LReSDFrfWbaQgqPbfVPlkMYkJp5Ly3ERKVkQJf2twpXjk3nCYlnMbG1knq2r/DFQqpjC33du07paIchqujQxBKl1DTX9+GQjvbwZHh7TGw17SlJLae7br/vlHh2wxzQw3CtzkLq1sT0n6a+5kpbxWwzb2lO0HmIO56beSAl3R6FiK6pSYxM55Z9K8VIuLlJZJuHc2tu5s+UeHa89BOMSxF0bZI0NelmKd013THc6TXKlFNIQqvMk2nuTO3u1BgjCbYtUhLadntKbLXtKcX12Cct2sTWeehr0s16SE2z6Wjhnt/D/X1X7t4uSTJp85J86mtKoTfvu822axdJauFade86Vmbmvbi5WvJp85I0MXVNn6a0Y7O9NE0WrqlYu1RS6PdU06TTJX0mPt3V71xlXk73I5/u6q9VV6eke6VLcZV8kqtLcZWkp5u6Jm8vpNqaOtFfUdLMfFizkFSc6hfH1TSSr2tc7SJJqzCfzc37Owm1zqFpn4W+sx6SMkN/Z2VSZpfpWJi+rnaRpHmoX1yN7lJcJelwsj1/npgkbumsOqUwaewrsb6yC5V5Uqol2ZsaONYpaQzTnmoaV6dUUlwlqZ2YeS6k0Nv6JdUpYXvd1S/pOcQ9UutX3df7Y+Eh3X2gtfbE3U2XsiMAAAAAAAC44y3t1FPj3vujuwzGn7sCAAAAAAAAl4xv0t0H7n4S21p7TNI7Lml3AAAAAAAAIL2g9/74WIPxTToAAAAAAADgkvFNumugLxbqLSwGfUdYdNUuAB0WCU0LK1YWyKssptzCAtB9ut0eF7WcmoW6FyE44iC0m4U4WwiDcGuSxsXM04LaZlH8fa11Wmbedly/2+x0deFXN8ZBOKAzs+K7azur/aYJeLg58dfWTXPNHSqNO7z94XCQ3GLKN5u/VtMCyTPTfzLCEpertMiyWfF2FvrOzbE4josj+zGWadV/48R8rquQQHPU/BLXN8wYRyZMQpJmK7PAbngfB2FxcXfupzAWd03FkJcrM8EYKTjCzbWFvlJY/D7d5yqLN4d2twp/une5hZNbWEw5LZzcXf/U1y2yHBZN7iZkQuFa7eb+eVZ/y9UvKSAr1DrNHPueAq5MGFZP4/q9UDN1Sk+hXKbW6Qdh35ZhP9xi34Xao1LTpPY4jVzh6aUyJ6a+bg5OdUqa2929INYvLlAr1CmHsfbY7p9rHbO9OK4/RjfNuX8jXFMzk+A0DddOqnV2rWtSTTMPNd7M1DqT0Hdqxx5eu0jS3ARKuM9UkuYunCGFVplap1pHu3Pf1fJSuM5sz7PGcI1hkAtmdznOqYXaoxSoVQyOMDVJc/f80DcGR8b6Zbs9Bly5+iUFRKRax9U1qdZxx+5U7dL7PZ697Ihv0gEAAAAAAACXjId0AAAAAAAAwCXjIR0AAAAAAABwyXhIBwAAAAAAAFwyHtIBAAAAAAAAl4x01+ugd90zPsskAUlSNylKLaUMhXQll5TSUkKsTUpJ6W8pIW/4GDaZtZgA6A7tKClmF5x4lvZ5X7tRSYxMqU2ufRJS02zf8O7SGFM7RuhbSGNL+zEzl1RKMZuZJLRKiuu6//bY05T8XBCudi1NUrHCMVqZvmnPUvuh+VznoffUpEumzzqeL6a/T26TpjbxrJBwWZSuqb1x23Of/z7tc15O9yMjJne7j7ty/wsJa/leaa61mLa63d7jvpmk1JCOVjlukalfeji3wtRn++fPybS7YylJPcxG7rM6KHzW1VT4oeNKtevyqiS2uu1d8By3zznV3QvcPUPy9xh3L0rjSv5el8fY3p67167H9Vz7NERwuhroIFQZY9QvVkifTrXV3Jwas3C+VFLvT8I17GrNaehrP+tYc7u6aHh9Hsco1PjJRafQX3jo/Qh1Si1Ndnh7Soj1967i7/SVMUx7rD1CyqxLco0p9Om+/94X7vck4Zt0AAAAAAAAwCXjIR0AAAAAAABwyXhIBwAAAAAAAFwyHtIBAAAAAAAAl4zgiAdFWIze9y0+u50U+k/cQvJ7XFzcjF1ey9z0j2Nc8DrpFRe9hvtFc8EDVTbooPL6cAKsQvvSrQ0fVnldmfZVWLw5bW8SIx52swyLLLt9XoaFVtP7Lu2Hed+V82K1x/9utUwLzF8XNiTkgvfhisxxPdzTbGu6/7l2d/+sjpHseh+OdUMx7MKO7c6tEa6nMY5Paq/UOoX9SOfWhbsK13ac2y96xff92Oc9o3Kvc/fQZTwPhx/7dM+fmLkhBgyE2mNfgRIu9CG1r1KtY9pdLbgeN9SOO16Au9a60jg191V2pd9e2De7z8X34ce42N/TY63jpGcQKVDCzg2FwMw9h0Wcds1/awAAAAAAAACuPh7SAQAAAAAAAJeMh3QAAAAAAADAJeMhHQAAAAAAAHDJeEgHAAAAAAAAXDLSXa+BdjBVa/f4KFPSkUlQaSnFZTbzQ0/N2K5NUnMpLKFvD+22f0h36dPhCWu2r0Ka2hjJOvdhQqyTgm66ORgpEyelRLn2xconlLq+8+77pnaXeDbv/to6MWPM0vZCItjcpLNOw1Gy7SHZLHFpatM9nnAuxWweUpTm5r0chZPrOJxIc5OGdxL+W5T7XN1nmvqu27f7x3PLnJ+rkN4X2wvpta7dXZP3q9JbGWG+7jats7Y9mXtMTIJ197SUYjYNqc0rMz/M/MVj9yJN7m6flynFdU/nXKhpbD0i2WPX0nGbmes99U3thc/PngOhHqmcc/acTX1HqEeu0fRi58rKXJtSOStzftpe5b5TqV9i37Z9z56F93GcUhJtXRPmDJvk67tOw/y56Gk/hkkp9Cnd1dUv87DTc9uWPutUv2x/fpX6JaX77lrTSL5GL9UptmeuXy4wbHMcI/xOWJmvU42Rfud19/fUd4zf6dvK3JuXw8do4RrpCs8s3BWY7pWrs3/van0hLfxLx8A36QAAAAAAAIBLxkM6AAAAAAAA4JLxkA4AAAAAAAC4ZDykAwAAAAAAAC4ZwRHXQDs8VGuH9+gUFkW0iykXF6eebLe3g9D30CzkWFm8WVI/2N4/1yb5RSb7QVhEMywc6dZMjX1Ne1rgM46xp0WdR+HW0Kwsphz6LsPCnwuzoGhaSPd4tX1uHa/8ip6uryTdaje22iZhoeCpCX1I6yDH/xxi+i/NIs3r9u3Os+ZXzD0MS+9OzKLO0xFW+3aBFJI/HGlBZreY8lFcONsf0FtmgeTbq+3PdN13e8486v68OArni+uf+h6vtvfNtUnSonA9LE3but0FR9iuV2dB5h2Dd8JhK7WvwgLJk8n2wUh90z20H5gDOvP3v2YOfg8fSOn2kBaRNu0xWGFp5ijXJu3vJEqBG5WAh7SQ9YG5Ll2blOsU1x7qom7OgUpNI4VzMZz3ru8Y107pRLzg2mWMuc/Nqev27YPh5mopz+2V+4O7x9yY+PvObOWvy1nbro2mKYjKNRe/4rE0g6xC3TAzoQ+23pJkpuV1/x1PsEpNI0knLiQr7NuROQeOQzjDrVCT3DJ1ze0eap3Vdq3jXi+lOtqfh6mOdvVZuh5c7V8JmZDCNXwVahfVfp/LvxOaWiD+Dmrapv5gxDHM78g9fH72PhX6thBMY8MZwiMNmykTQ6RCTeLuoSn4yoVvnbrWW18RHAEAAAAAAABcZzykAwAAAAAAAC4ZD+kAAAAAAACAS8ZDOgAAAAAAAOCS8ZAOAAAAAAAAuGSku14D7eYNtcnNe3QKcTIuFSUkvjSTBCvJp56lFDPXHtPRhiehrQ5939XMJAeZNklahdRX174KAXI+Wcf3rUXyXbBCilJKTXNJTCmJcr70B+lkun1unIS0smPTnlKpbpu0K0mamASxaczz2rYM18gq/PcQ1//QJJtJ0rFJY5uFxLPUPjWRV5MRYrDSEVqakzkl+c7NMUoprkchJcqlraYUM5f6WukrSbeXJmEtnFs+NS2c9+F6mJv+y5SEZq411yadEYQ2QvKvtafE1nJCpUvuDomtq5lpTylm6YIwk+XEJZtJWpnyrMX7eLhnu/ZKounCx5a5lNl2RdJdU1q8PRYpsdUdozCuTXEN7bmmMddqSHeN9YutdVJivbmPF2qa1L63JNgxpJomdK/Mn24OdnO1lOf248l2e0rPvG0SB2fh+nP3/GQVPkBXp7h7uyTNm8tpl26a+uUo1S+mopiGJPv0/natayo1jSTNzQVUqV9O5M+LmE5v6xdfe7gxUp3i2nO6q28/Me2uTfK1f/o9Iae7urZC3OpFK859tk6Jyd1m2PC7bQu37OaOf7qczM7Fb39V6pe5rz3cc4i2CG8k1S9m/myhDpNLvD6d7rrq0m3/0jHwTToAAAAAAADgkvGQDgAAAAAAALhkPKQDAAAAAAAALhkP6QAAAAAAAIBLxkM6AAAAAAAA4JKR7noNtJs31aYPnWoopNe4JNeQUNlDKpztH/p2k5BWSTxL/WPi2eHwxLPYbpJxYhLajik8qX9MTXO7XA0vcqE2McnHNKUkptX2jiyWIdkznC/Hy+0p6iAkgrl2l9ZalRLP5iahy7Wt2/2H7dJIb4Z0tJlNd02Jbin1NaQg2TGGp6OlxDN37HI6mkkEKxy31H4UE/Jc4llKd/VJaE8st1O1by1CGttye4yjZXgfod0loaVkZJc46K5JSeqhfYTgX68SvJb62iTtlGYZElRN/xb6uoS0dDmlOcMPHJLXpmY+C+9Pqd2lqackNJd4Fs7DZvqeEfPt2yupr+4YhRT6WP+YOqVU08TE+kqdMnyMSk2zHsNc7ymx3uyGuxaks2odUxel096lNl94umtoDnOfmytTkradlyf+Hnw08dfUgemfEk0rVqU00pBIa9rTPTjVL679MEygrtaJKa6FJPsk1SROOkZLm04fUtpNrRNrmlC/uHTXSv2Sapr3mPrlPSbFXpKeTPXLwiTZm1pekubm94f0e0KuXypppAPbzlKYzyq/z6XfCZutdcL23Nye0l3Dr0fNvpkRappwz7aJreE5hA7MnLEMbyRcO831dzWNlOuX94515r/eGd+kAwAAAAAAAC4ZD+kAAAAAAACAS8ZDOgAAAAAAAOCS8ZAOAAAAAAAAuGQER1wD/eGH1A8eOrtTZZHlsOhjWpTbL8icxjB948LLYfFQM0YOgxgjOGK7bVno6xZ0lhQfkdtFReMiy9uLWo6y1nthkeUeVkxdmsVfF2Gh7rT4/bF5f5PCAsKrsG85DMIspmzCTiS/APTR1E+paUHfG5PtxZRvthPb1y2yHIMjwoqwExMokfpWLAsL5adwDbd4c+qbFlN2n19akPm2WQz5OGzvtgl9kKRbZoy0mPKthVm82bRJeZHlY3OdpEWW3fWXFkkvBUcUJpi4mHJh0fjKgsxx4eVU6ZhT3y+aXBODfty9MmUXuIX5w321hUCCycK8QdcmHwbRlyFEY4zgiF0VgyPcsU+BG7ZvqFNWqX5xn18cw4Q+hICISv0S6xTT7mqXM9sLIVmlQJjqnGE7D2xTbU50c6okLcz5cjzx1+R04g/oZMeQiBQQkWogV9ek++rR1ARcrVLAla9JbppaxwVESD70IfVNdq1rUk2zCoW7618L4gi1QPpMTK3iapr1GNt9cxjE9n642kWSjkKd8qQLjliE4AhX04QwunT9ud9B9hmGZYdI16+5H8VyOc2f7v4QA5lcERXu43FSNXVKeH92c+HePEkBFgcmtCr8/t/M+dJCnRLDIFz/VL+49lPHvod6aix8kw4AAAAAAAC4ZDykAwAAAAAAAC4ZD+kAAAAAAACAS8ZDOgAAAAAAAOCS8ZAOAAAAAAAAuGSku14Dq0duaDW7eb4Xu+SZFIoTkjldylpKd3EpNamvSzxLY6S+NjWtkMwq+YS0EALptxcSB9P2XPJPTgv07SU2CS2cBC7dNYTbrFy668IfjHlIDto18WwR0rVSEtrCHFCXjCVJD02308pur3wKlks2k3wSmkt8TX1dCtq67/AktF2P8Vlcstwy/Lchl4RWSeGVfDpd+vxce0pSS4mtbozbJtlMko7MGCnF9SgkoZ2Y6yclI7sktJT0565rSSG6y3cdg537wvzp5p04F6WE7ZToZbkk9DBsTKQdnti6Wphk65Bu10xfSepLc0BiYut2u2tb75zpmw5l6RgXpBTXGFhXqFNce/qcQl1k65SUJluqaWxzSGwNyetmigpTX67DXLprulYLSbCj2LGmkcK9K1x/c5Pk2kaoadJ87dIv5wf+4Kd7zJPT7Vol1R4uCfTGxNcYlVpnEhJYfa0TEuv3VL/Ee2XgkuhdYv267/ZnlWqPMWodV7+cFPrGFPpQp7h015NQp8xNTZN+T3C/U0jhvn/BNU3ld7TUN83t9tRPSanmfpvOw3TTtnVKen/m/tcP/LirRUp93R4j1R6u1qnUKZLUTF103vplNfdp1mPhm3QAAAD/f3t3HiZbUR98/Pu7l8uOKCqILAKKiju4Lyhq4r5EBQVxi9G4xEQxIvF9jcFdVFxiXOKKC+7gStx4DYqQiBiIIIIQZFFUwBXZ7731/lGnM4e5VT1dM93TM3O/n+fpp/ucrqpT5/Tp0zU1deonSZIkTZmddJIkSZIkSdKU2UknSZIkSZIkTZmddJIkSZIkSdKUGThiBVi71RpWbVqe2PN/NcyBmmoTMte6dFsmZC5OINwYOKJQRnVC5lLa1gmZC+urZRTmfq0FmahOslw6zrVjX6pG23y3bQpzaKbKZK6FuY1ZV5m0/LqoHIxSuZUJfUvrS4EgoDzxMsA1qzb8ADddXZ4g+cpVpaAP5bSbVcrYpDBBcmnSZIBNVm042WnrZMqrCxM1r5pgJIBSgI51lc+k9PnVJk2urV9biNJyfWV7pYmTa5NsX1eJ/lIKBlEL5FAq+9pK2lKAiLx+wzJqkyyvK6xP6yrXydoky4XTKyrfv6b5u6uTLG9YSH1y/8K62jW1UrdU2Jf1te9Dab7phrpB+bekcslgVeGzWl8L+lCZuzgK1+bqhMyF1VHZXukQlSasrqUdi+q84LVJsgtpa+2Gls+6oe1RDRLS0E6pnePFoFUNbZpS8Iq8vlxGSzCI8uTplROjYehAfbLvUoCryiTplWtf6Vq5jvLBr51yLYqBKirtlNL6WiCAzQsBrqAcIKDWTtm0cJFaUwn6sEmhXZTTj97WKbVfSm0XmFz7pRZcrKX9UguSVQz8UWvTVNYX2xOVc6DUfrmuoa1TaxfVAkeUgkSUAkRAua1TC9CyvvJ9KP4NMsEAV01/d7UEAqy1X0p/g1Y2V6pINbhKw9/69TZNYXuN7ZRVDcEgSmU0B44otmsbTphe0rXXXTd6vnlwJJ0kSZIkSZI0ZY6kWwYi4orZq6ZSEUmSJEmSJA38LGJm6GJKaZuFFOZIOkmSJEmSJGnKHEm3DMzuiY2ImwOXTqk6kiRJkiRJgt1TSpeNqzBH0kmSJEmSJElT5ki6FWDtVpsQm83zoyxGMRs9LZSjqbVF+WqLkLfQ6K7V6GjVyDqFqGktUXhq0dEa9q8WCa1YRuuMhaWia6GDCpHQUiXt+lIkyVj4/wXWVaJLrivUrRbF9dpV5Q9lzeoNQwdtWotWVki7SSUCWS1CbCl9LRLTmkI9apHNamWsqtRvUtZXLyYbKkXirX3WpSiuUI7IVoo0lre34fra+XJ9ZXulSK7XV8ooRmatRDGrRYgtRT1bt7YWCa3wXa1sr/p9L4a5rKQtqZyHqXKRKgYba/kdGLFaw0RTNNny/tV+j0qXkloUz1LEslLE15y2uLpYRvUSUIrY2hAdbaLR9FrU2ikNbY+WdlFL+6XabiiU0d5OKaUdPWJrPWrs6OtrZRS/w7XPqVJEU/joUtLKb0ktAmDpWrm+8uWpXYOL5dai0xd+N9aurvzubLJh2jXrynW7dnWlrVO4GJUiyEM5Cv2mhfYP1Nskpaivqyufaa0dVdzeEmnTlKK71iLEltoTpfxQbhdBOTprLW2pPVGLWF9KW20XNbRfahFbS9Hpq9FdW6LTN0Ssb/7tavp9KBReuy7XgpSWosJXqlb8NlR/oyrXvoZo46Voq5U/mepR6AufVUvbo9r+aYg439Iu6lt77WS70RxJJ0mSJEmSJE2ZnXSSJEmSJEnSlNlJJ0mSJEmSJE2ZnXSSJEmSJEnSlNlJJ0mSJEmSJE2Z0V1XgOu3XAWb9/pbW6N7zlILglWP7lpYN5aoaaOvr0UVK6etRNeqldEQIbYcmbWcthr1taGMhX7WNVE5CYqRXCsRe0pRf9dXQzWOHrVp/SblMDyl6GjXVyKeVaOYFaKpXVNJW4pMVit3dTVq2obroxbxrFBGLYprTS3yWjFtQ9S0liiutYhnxehotUh4tQi/hXrUIpOtK6yvRw4ePbpZ6Tys1aMWxawe3ay0vcp3pxBxMFWio1FbX4ySWE7apCEC5/pKBNVV1R+qQhkN9agFoC6trx3PqEXEXDd6tNUolF2P4lpZXzhGtahpxYhntWNcTFtOWtUUJbih2Ib2S1PaWoS8pmiylbQNv/m19ktTu2gcUe+bytjww55om6blulWLGFlcXblel/JXo7iOvn7t+vIXuxQVfFWljbFm9ehtklq01VLaaluncjEqtVXqUegLaRt/eBa7/VIutxadvvDb3BjdtdQmqbZ1CmXX2jSlc6tWbq2t0xaFvtCmqX0nW9ovle9Z9TdtoRp+B6hdJytFl4qoBa8tt1MqaWu/JYX0UWi7QLntUfv8WtovpYivuXKltJW6TbL90ilF2h4nR9JJkiRJkiRJU2YnnSRJkiRJkjRldtJJkiRJkiRJU2YnnSRJkiRJkjRlBo5YAa7fMmCL+U2G2TSH5hgmZC5OollNO/oEyS2TdrYEmaiWMY6gFtX1o0+yXJwfufVUKEb+GH1y8epsp6WJ6yvlri9GpIC0ujCpbG1y/8L2VlUmTV67qjJhcWF9bULm0iTLtaAPqyvbK6WvBZloCRJROwVaA00sVO2zKimeWrWAIg0BJdbVAis0TKZcDWBRCmxSnWS5kLYSIKK6fw0TJBcnWa4GiKhN9Fu6NpSLKBfQuL50OGqXomJAidEDFkHl0FWuy6tKky+3BnJoCQbRcK2tT5xcWFn7rFsm2284ByZ1yWn+nSulb/k6tH11yu2XhrhJLUEmqulrQatKZdT2oykoV+XDbqhb8zWjpDi5eFvgj1I0iPqpXJgov/KbX7+2F8poaL/U2inX1do6DW2Pcjtl9HJrZVTTFteWLZU2TUstyu2UtrZAqV1TTVsKytXSTqkFO6kFuCpsrxQEDsrtlFQpt9p+KZ22Lb9djT8mTcmLu1KrXOU4l1LWgkG0BMOqXdtLgRxagkFUdq8eOKKQoaGMlgARQ9eXtjdH2utLB3yMHEknSZIkSZIkTZmddJIkSZIkSdKU2UknSZIkSZIkTZmddJIkSZIkSdKU2UknSZIkSZIkTZnRXVeAdVsGMcXorgstoxbFrGl7Y4ia1hQhtlZGMW0lyldDhNhqd3ox0m1btKtS8lpw11KEtOrWGqqRKgcjFaL+RPV4brh+XSVKVK2MYqSwyvEsRYKtBfqpRX0trW8po1ZuzWTjEG2opXa1yGQtaUvnbTWaXkN0tGqw40KEtKa6VSKeVSO2ltbXoisX0zZER4PGD7AhaUsEx4Z/Jaba96yy36WoZ7UIZKVLSVNkVqhEnVx42qboZrVre0t011HzT8HE2jSNUUcX2vYYR7uorZ1SKbe6vvCBt0S9H0d7chzXp9r3r+HXsvj7UIruC6yvRGKmcDxrUTWLERxr175qxNbRyyi3dWptmuLqprZKU9qRU47HONo0pfXVNsYYyiieR7V2UUNk1mrk9VLbqqWtU0tbjU6/4arab34xbetvVylD5fs+YhWGv1NsF9XaiIXsjVHo29oehZXN7aLRo96PpU0yxnbNOqO7SpIkSZIkSSubnXSSJEmSJEnSlNlJJ0mSJEmSJE2ZnXSSJEmSJEnSlBk4YgVYuzmwxTwzN8x5uOhBJhrKbpqQuXUC6IbgDOW0DeXSNml10wTXLRomWY7KBosTS9cmfq1NElqYlLMaiKM4uWrb5Mbri4EjalVrSVteXzqPmr5mrTPeLvbs7gv8wlcnTR5DHYqTLDdO3lxKX59kuXQxKietTpxcnKR3DEEDmsoYwwWmdh6Wvu+1za0u5a+kYWFJkwAAIABJREFUrV0zipMsVyZrbwqwU6nHIgZyGEcZY7lcTOqSM4bTcCkEzqqWMZZ2yujbq/2utgTlqqYtrWz4TWxW3MG2oDKpeHFvmPy+2k6plFFIX2+rjh5QqxgdhzG0X8YQtKq5/VIseOm2aZqCYTVub8Htl4ZJ/Gu/c5Nrp7QFqih+h8cRDKtmoX9XNgYFKv/NW7ueFdJWKtcUOGJCwbBaTaxdVDNH2rUTvvw4kk6SJEmSJEmaMjvpJEmSJEmSpCmzk06SJEmSJEmaMjvpJEmSJEmSpCmzk06SJEmSJEmaMqO7rgDrNoXYbI5E44j4WTGO6GYtacuR0EYPsdIa3XXBkXzaAhUVu85bonyNRXWDhe1VIvlEqYzWCHmlfyPUokQVj0UlMlbT59cQTbZmDJ9fNULsQk0qmt4YVKOKjaMeY4gS1RSxdRyRUpuig7ZEk61tb8MyWk6X2u5Vry6lS0YtYmSpHq3/dixda2tRIFvOl6aIraNHk61qrMeoJho4saXsCV37JhVtfuj6UevR+rszqajwLe2XSbWhhqUvFdESdbtynItR66t1a4gEW9uPUuT7cUTBrp1HxbS1MirrG7ZXTDrBv0sWfPGaUJsGGts1jZFOR07bUu4k22GlCOut+9z02zyhz7X2PVs1ehuqFG0eKEeyr5VRaNNE5YSrnoct0V2L+RvbNGOMtjrS9haqV+66tRPaRseRdJIkSZIkSdKU2UknSZIkSZIkTZmddJIkSZIkSdKU2UknSZIkSZIkTZmddJIkSZIkSdKUGd11BVi/WWLdZvMMYzKxiJELL2Kxo8ZW6zGh7bXtXy3K1wLLhXIkn2rShoiRpTrXAp7VwnyVoh2N5Xi2pJ1kCLLS9hqisTWXPY5CGkwyOmRxe5MKAzmZtNUoZi1ljyVi1hgijI7Dqg03WI9A1hA1rba9huM5euzotjJaQv2N5fNY9O/khMpd7GvZUm7TVNI3HfpJRamFprbApNSvDbW2R0uU0tEjwY7leDaV2xhldqEaLlK2afrlLvKOTKhNA42R5Vu2N4aIrU1R2hs2V69BoU1TS12K4louoqnOtSZGyxnXFpG47YAuuF0zxSj069dOtkHlSDpJkiRJkiRpyuykkyRJkiRJkqbMTjpJkiRJkiRpyuykkyRJkiRJkqbMwBErwPpNINZMuxYjWimBKmrGMbP3pCYWbjGGiUZLlWuub8PxjFLtxvKZjqGMpu0t9gZXkIlN6rxEyl3gF7758jSh/a7tRrF+DZPRN85XvPAL6CJPWDyWzS124IiVYpKX5UnFu1nsNs1SaYeNo8ot7ZeW71RLm2aSbYGJtcVtv8zLcgzoM5YyFhj0oXl7YyijVGxLfJbGCFdN18SGxC0x/5q0BuKYYNmTtn7CvWiOpJMkSZIkSZKmzE46SZIkSZIkacrspJMkSZIkSZKmzE46SZIkSZIkacrspJMkSZIkSZKmbKOP7hoRdwFeDjwUuAlwMXAs8MaU0u+H5NsN+AfgUcAOwK+A44DXp5R+Mdla31DaJJE2WWIhTxbbMgwoNbHIrJO0wNOsGjWoGjVt+gdpCVRBS9RYoo1NylKuW8WCv2stUWOHvjEao61qYpb4786K+V1siQTb2n4pbm9xD9yK+Zw0VrZdFkfT96+1/bLgxGULLmEFfX5zmXTfy4oZSRcRm0TEX0XExRGRIuKoEfIcBJwKHAxsAZwN7E7utDs9Inap5Hsw8CPgecBNgR+TO+peAPwoIu42hl2SJEmSJEnSRmJFdNJFxP7AmcAHgZ1HzHMP4GPAGuBlwPYppbsBewG/BG4FHBMRq2bl2xX4IrAN8Hbg5imlfbr0ZwLbAV+OiC3HsGuSJEmSJEnaCKyITjrgjcDtgE8Cx4+Y5x3k233fl1I6MqV0HUBK6VzgL7o09wQePyvfG4AbAV8DXpZSuqrL9+su7Z+AXcij7CRJkiRJkqQ5rZROujcBe6eUDgbmnA8uIu4I3L9bPHL2+ymlU4BTusWDe/m2A548yJdSWj8r3/nAv83OJ0mSJEmSJA2zIgJHpJQ+1Jjlsd3zD1NK51XSfAW4F3C/3rqHkW+PvRT4zpB8TwbuHhGbpZSubaxbs7QqP6RlawlPNLqkJ9iVNKclPXn6Uq6bpLkt4TaC7RdpeVvS7ZeS5VbfBZh038vG2rWzT/d86pA0Z3TPO0TEmln5TksprZ0jH8At51k/SZIkSZIkbURWxEi6edite/71kDS/6Z5XAbcALm7MBzmIxc9aKhYRVxRWb91ShiRJkiRJkibu0ogbDCX80+wEKaVtRi1sYx1Jd6PuedTOtsGIuPnmkyRJkiRJkqo21pF0g1tVVw9J0z82Vy4w38hG6WGNiO3pdRSuu7J5M9KMjWj+AEmLyPmQJE2KbRdJk2L7RXMo9L/skFK6dFzlb6yddD8H7gjsMCTNTXuvBxFjf949t+Ybt+36C5e86S0T2owkSZIkSZIqtiMHFx2LjfV210Fn2/ZD0tyse74G+P088gFc0l41SZIkSZIkbWw21k66i7vnew1Js3f3fFFKaTDodZBvn4iojUIc5LsGuGz+VZQkSZIkSdLGYmO93fW73fNdI2LPlNK5hTSP7Z6/2Vt3Yvd8U+DBwLeG5Ds+pbR+wTUtOxfYq7f8k+552Ai/5WoQHXf3qdZicty/xSljUpZy3cbB/Vsck6rHUtm/SXH/lrfBbSErse0CK//zc/8mm3/Slnr9FmIl7xssnf2z7TI/7t/yNmi77MANpyAr9SfN28bcSXcBsBvwUuAF/TcjYj/gzt3iUYP1KaXzI+JEYF/gkIj4f/2OuIjYC/jzbvEjk6k6pJTWAWf3tjtYv+JG7kVEgpW5b+D+LVYZk7KU6zYO7t/yrsdS2b9Jcf+Wt5XcdoGN4vNz/yaYf9KWev0WYiXvGyyd/bPtMj/u3/LWa7tcyhjnoJtto7zdtevkOqxbfH5EvDgi1gBExO2Bz3bvfSql9MNZ2Q8D1gOPBI6IiC26fDsCXwQ2B04GvjDZvZAkSZIkSdJKsVF20gGklD4LvLlbfAfw64j4L+BM4ObAScBfF/L9B/B35ODMLwMui4hTyUM7b0se4bZ/bx47SZIkSZIkaaiN9XZXAFJKh0XEfwCHkAM+7An8N/Ax4D0ppesr+d4dEWeSR9XdB7gDcB55BN6RKaUrF6P+PX9a5O0tmpTSNtOuwyS5f8ub+7e8uX/Lm/u37K3Ytgus/M/P/VveVvL+reR9A/dvuXP/lr1FabuEA76Wv4i4AjaKL4U2Qp7fWuk8x7Ux8rzXSub5rZXOc1wbo8U67zfa210lSZIkSZKkpcJOOkmSJEmSJGnK7KSTJEmSJEmSpsw56SRJkiRJkqQpcySdJEmSJEmSNGV20kmSJEmSJElTZiedJEmSJEmSNGV20kmSJEmSJElTZiedJEmSJEmSNGV20kmSJEmSJElTZiedJEmSJEmSNGV20kmSJEmSJElTZiedJEmSJEmSNGV20kmSJEmSJElTZiedJEmSJEmSNGV20i1jEbFvRHwpIi6NiKsj4syI+D8Rsfm06yYNExFbRcQLI+IL3Xl7VUSc153PD6vk2TQinhYRn4uI0yPiTxFxYUT8W0QcFBGx2PshtYiIm0fEgRHxwe58vyIiHjHuPNJSFhF3iYhPRMQvI+KaiDg3Io6IiBvPkW+3iHhfRFwUEdd21//3RMROi1V3qWa+bZSIWB0Rj4uIoyPiBxHxh4i4JCKOj4gXRcTqxd4XaRQRsW1EPD4i/iUizurO3edNKp+0EJE9PCI+HBH/GRF/jIjfRcSJEXHQHHkXvd0SKaXWfdQSEBGHAkcAAfwauAy4U/f2acCDU0p/mFL1pKqI2BS4CNihW3UZcAlwe2Czbt1bUkovn5Xvh8A+3eIfgPOB2wJbdes+CxyYvKhpiYmIrYDXAC8CNu1W/4J87X5hSun748gjLaaI2AR4JnA4sDPw0ZTSs+bIcxDwUWAN+Tp+Abntshq4ENg3pXRxId+DgS8B2wBXAecAdyD/ZvwWeGhK6fQx7JYEtJ/f822jRMQxwBO7xauAs4E9gMEffyeRz+9rF7ZH0g3N5xrey3cY8A/A1t3qS8lt+VenlL44znzSQkXELsBRwEO6VT8Dfg/sBQwGN30ceHZKae2svFNptziSbhmKiMcDbwYS8FRgp5TSnYH7AFcCewMfnF4NpbqU0nXAf5LP0dsCO6SU7gZsB3y4S3ZoRDxoVtaTgC8CdwW2SyntQ27Evq57/8nA0ydcfalJRNySfO6+FLgaeAP5Gr1LSunulQ665jzSYoqI/YEzydfxnUfMcw/gY+SG7suA7btr/17AL4FbAcdExKpZ+XYlX/u3Ad4O3Ly7/t+qq8N2wJcjYssx7Jo0r/Ob+bdRTgJOBO4HbJtSujv5nH5B9/79yR0b0tjM8xwnIrYBvkE+rzcF/pl8ju6YUtp7SAfdvPJJY/IH4J7ACcCdUkp7dNfobYAXA2vJ1+cbjOicarslpeRjGT26k+Rn5A66wwrv79+9l4B9pl1fHz5KD2CLyvrBKLsEvGvEPAH8R5fnK9PeNx8+Bg9gS+AH3bl5JrDnJPL48LHYD+Dc7hw9GvhW9/qoOfJ8r0v33sJ79+q1XZ4w671PdOv/DVg16709gCu69w+Z9nHxsTIe8zy/59VGAbagu7Op8N6nu3xnTPuY+FhZj3me46uBL3dpLwbuNeK25pXPh49xPsj/CFldee/d3fl5af96PM12iyPplp+HAbsB1wDvKbx/LPl2KMij7KQlJ6V0dWX9deQOCoBbj5gnASeX8khTdjhwD3Kj9H4ppXMnlEdabG8C9k4pHUy+DXuoiLgjedQEwJGz308pnQKc0i0e3Mu3HXkEEsCRKaX1s/KdT24E3yCftEBN5zfMv42SUrq6S1NyUi2ftEDN5zh5lNFjgT8B9++u25PMJ41NSunklNK6ytuf7Z5vTjeydNrtFjvplp/Hds/HpZSumP1mdyJ8tVu836LVShqfPbvnCyecR5qY7pbVv+0W/yal9MdJ5JGmIaX0odQ2B9yg7fLDlNJ5lTRf6Z77bZeHke8guBT4zhz57h4Rm1XSSCObx/k9l/m2UWzbaCJaz/Hutrx/7BZfkVK6aJL5pEX2m97r7bvnqbZbNhlaXS1Fg0lpTx2S5ozu2YhnWla6qJV37haPGjHP3sCfdYsfmUC1pPk4iDwZ7ZnAVyPiz8mTg9+K/F/r04GjU0q/X2AeaTloabvsEBFrUkrX9/KdlmZN5lzIB3BL8pQg0pIw3zZKN6fRYDSGbRtN2yOBW5Dn9npfRNwHOBC4DbmD43TgMymlS8aUT1pMd+y9Pqd7nmq7xU665We37vnXQ9IMeoNvGRGrZg+zlJaSiFgD7E6+PfvQbvXhacjk+BGxGtgFeALwf8lzunwY+PxkayuN7NHd8ynk0c2PKqR5RUQ8K6V0/ALySMvBbt3zKG2XVeQ/6i5uzAf5NhU76TRV822jdJOP7wg8gjz6aAfg6+RJ9qVpGrRPTgLeD/zlrPefAbwqIl6YUvrUGPJJi2lwnp6TUvpT93q37nkq7RY76ZafG3XPo3zwm5DvrR6WVpqKiLgx8LtZqy8CXphSOm5Ivt8D2/ZW/Y78I3+0HdJaQu7RPR9Ijrr9XGaGt+8LvI38R9wXImKvlNLP55lHWg5a2i6Q/7N88TzzSVMz3zZKRJxOjgw7cC3wEvKE5ddNoq5Sg0H7ZD/yhPcvAz4DXEWeQP8twJ2AT0TEOSml/1pgPmlRdBFcB5G339J7a6rtFjvplp+1wGbkSDk1/c/1yslWR5q3teR79deQ/3N8K2BX4PXdCNCvVPJ9j9wA3oH834qbAK8C1kTER4ZMwCwtim4Olm26xWuA+3QTxQ58PiK+B/y0S/fGiHheax5mGhXSUje45aO17TLffNK0zLeNcirwR/I/13cjT31wKLB5RLx1yITn0mK4Rfe8CfCAlNIPeu99PSK+C5xFbsu/k/yPxYXkkyYuIrYCPtAtng18tPf2VNstBo5YfgYjJ3YYkuam3fMfe0M2pSUlpfSnlNJ+KaX7p5T2IHfQvZn8n+QvR8RzK/kek1LaN6V0W/KP/8vJt8t+CHjtIlVfGqY/iuKNszrbAEgp/Yoc8h3yf5jnk0daLlraLjATbXC++aSpmG8bJaX0nJTSA1NKe5E76p4L3JgchfOjtXzSIhm0UT4wq6MNgJTSVeR/HgLctzcZ/nzzSRPVTUtwNHA38sjOp86aQ26q7RY76ZafwQe//ZA0N+uenYRTy0ZK6RcppcOYacS+LSJuMkee36aU3gL8dbfqFRFxu0nWUxrBpcDg9qSfDEk3iKy2c5e+KU93y7i0HLS0Xa4BBsFRWvKB7R4tIfNto3T/xPwg8Lhu1cER8dBJ1VMawSAq6yjtk9XA7ReYT5qYiAjgrcDju1XPSCmdNivZVNstdtItPxd3z/cakmbv7tkw11qOPt49bw3cfcQ8nwGuJ1/THCqvqepuS/qfbnHXIUkH/7FbRx7u3prnmvnWUVpkTW2X3i2Bg3z7RERtipZBvmuAy+ZfRWli5ttG+Xdm/oB78LgrJTU4t3sepX0CMLiTa775pInoOuheS57zE+DvU0rHFJJOtd1iJ93y893u+RERsc3sN7uT4ZHd4jcXrVbS+Pyh93rY/fx9VzMzCmnUPNIkndU9P2hImkHI93NSStfMM4+0HAzaLneNiD0raR7bPffbLid2zzel3kkxyHe8wYO0RM2rjdL90TdoE9m20TT9b/uk6+QoGbRPrmQmWuV880mTcjg56jbAoSmlt1XSTbXdYifd8vNF4Apy8IgXFt5/GnmS2muBTy9ivaRxeWbv9Skj5nkSsFX3+vvjrY40L0d1zwdExAa3b3TBJZ7fLX5uAXmk5eC7wAXd65fOfjMi9gPu3C0eNVjfzc04aPAeEhGrZuXbC/jzbvEjY6utNF7zaqNExL2BvVrzSRMwmBfx3sAGt1531+bByKTP9zoe5ptPGruI+EdyIB/IHXRvHZJ8qu0WO+mWmZTS74DXdYtviIgDBh9+RDwAeG/33ttSSk6grCUnIvaNiI9HxL0jYtPe+k0j4hBm5qQ7ojvfiYjdIuJLEfGQiNiil2eTiHgaM5F5PgP89yLtijTMccDXyb+z346I+w7+ixwRNyP/w2VHcgPgiAXkkZa87hbww7rF50fEiyNiDUDXIf3Z7r1PpZR+OCv7YcB68l0CRwx+AyJiR/J3YnPgZOALk90LqWy+bZSI2DoivhYRj+vfHRMRqyLikcCx3aqTyL8P0lSklM4A/rVbPDYiHt37+3Mrcmfc3uQIxYctNJ80TpEdDrwGSMAL5+igm3q7JeqRwLVUdRe3o4EDu1W/Ai4H7tQtHwMcOCtCibQkRMTbmfmv2Vpm5qvYk5mQ1B8Hnj04hyPixcA7uvcSee6uq4DbkUeVAnwDeHJK6Y8T3QFpRF3gk68C9+tWXUa+Xt+R3BH3E+CAlNKPF5JHmqaIOIo8AvqjKaVnzZH2CHK0S4DfkTuc70K+le8k4BGlqPQR8TfAu4Ag3xJ1NrnNs1n3+iEppV8ufG+kGxrl/J5vGyUiHk/+g23gQuC3Xb4tu3WnAE9KKf0caQJGvYZHxObAJ4EndKv+QL499Y7AGnKUyoNSSieOI580Lt0Iutd0i78DfjQk+ddTSm/q5Z1Ku8VOumWqG13xbPKtT3cgz3XxY/J/6z6W/GC1RHXzJj4R2B+4DbAHeRL888gXrXellE6dlSeAh5Nv594TuDX5h/084KfAh1JKxy/WPkij6s73ZwNPJw+L35zcOD0WeENK6cpx5JGmpaWTrkv/F8Ah5NETQb6Gfwx4T0rp+iH5HkT+7/R9yN+J88n/yT7S74QmZcROunm1Ubp89yVf7/fq8m1D7uQ7lzzy7nO26TVJjf9oCeAA4Dnk4G7bkKNZHge8OqV0+TjzSeMQEScwfL7nvg2+B9Not9hJJ0mSJEmSJE2Zc9JJkiRJkiRJU2YnnSRJkiRJkjRldtJJkiRJkiRJU2YnnSRJkiRJkjRldtJJkiRJkiRJU2YnnSRJkiRJkjRldtJJkiRJkiRJU2YnnSRJkiRJkjRldtJJkiRJkiRJU2YnnSRJkiRJkjRldtJJkiRJkiRJU2YnnSRJkiRJkjRldtJJkiRJkiRJU2YnnSRJkiRJkjRldtJJkiRJkiRJU2YnnSRJkiRJkjRlm0y7ApIkSdJSERG7Ac/qFk9PKX1xapVZgiLiJcCNAVJKh0+3NpIkrSyRUpp2HSRJkpa1iDgBeFC3uHtK6YLp1UYLERH7Af/eLX40pfSs6dVm6YmIC4BbAaSUYrq1kSRpZfF2V0mStCARsUNEHBwRH4mIH0bERRFxdUT8MSJ+GhEnRsRnIuIF3SilUhknRETqHsU086jX1hFxYER8PCLOjIjLI+LaiPh5RHw/Io6IiAdEhB0NYxARO0fEqRHxm4g4eNr1Wc4i4qDufP2viNhl2vWRJEmLw9tdJUnSvETErYFXAM+k3KbYHNgG2LNbfjJwBnCXCddrDfAC4JXAzQtJduoe9wJeDpwUES9PKZ08yXptBP4WuHv3+m3A0VOsy3L3DuCm3eNFwGHTrY4kSVoMdtJJkqQm3ciz5wPvBNZ0qy8HvgacAFwI/AbYArhF97gP8DjgMxOu23bAZ4GHdquuAb4MfLWr1xXALYF7A08Fbg3cHzgxIl4K/HNyLpD5uqz3en0pQURsCzwEeDhweErpV4tRsWXoMmD77nXxWEqSpJXHTjpJkjSyroPu3eSRapA7vV4JfDCldNWQrO/tRrhtOsG6bQ18B7hTt+rrwPNTShfOSnoacFxEvAZ4CfAGcmfjO8ij/46YVB1XuKOBpwM7A39fSfMlZubue9NiVGqZej3wHuB88vdNkiRtBOykkyRJLQ5jpoPuHOCxKaVzR8mYUroeuH4Sleo6Dz/CTAfdJ4FnppTWDqnPWuCtEXEOcAy5o+4NEfH9lNIJk6jnSpZS+iVw12nXYyVIKX0K+NS06yFJkhaXgSMkSdJIIuKu5FFnkG9vfeSoHXSL4GHA/t3rnwB/NayDri+l9BXgtd3iKuBtBpOQJEnSYrOTTpIkjepNwKDz6tCU0s+mWZlZXt17fUhK6ZrG/G8CLu5e7w08eiy1KoiIu0fEOyLivyPitxFxXUT8IiK+GREviYgbzZH/rC4K7mndckTEfhHx6S6a7tURcW5EfLQL7jFqve4VER+IiP+JiGsi4pcRcXJXp227NIMovOcV8u/Xi9B7VGk9M7e6Avyslz5FxOG9PLv11p8wQt0v6G1jWLrNIuKvI+Ib3f5dGxEXRsRxXUTV1XMeqBuWt0tEHB4RP4iIy7ryLoiIr3SRhefV1o6Io3r7v1/h/UmdA/eLiA9GxHldGZd3EXv/MSJu0bgPayLimRHx5cgRn6+NiEsj4pSIeGVE7FDI82cRsb7bt1+X0vTSHtw7RmdExJYt9ZMkaSmyk06SJM0pIm4FPKJbvBD4+BSrcwMRsTs5EATA/wDfai2juxX3X3urHj+Gqt1ARNwkIj4HnAq8mBzl9ibk22xvCfw58Hbggog4eEhRg1uGt4mILYCjgH8HnkKOpLs5cBvgGcDZEVGbH25Qr80j4kPA94HnAHsAm5EDfty3q9NpEXHHXrbWTtCpi4h9gLPIn/PDyPu3KbAr8CjyLdLfBm48QlnRHddzgH8C7gHcrCvvVsBjyLernhARpQjDCzXuc2DLiDgaOAn4K3JAlc3J0WXvDrwG+HFE/NkolYuIvcmRnI8CHgvsQj42NwfuSR65elZEPKqfL6V0PPDGbnF74CMRG45q7a5H7+kWrwQOmGNOTEmSlgXnpJMkSaN4TO/1V1JK66ZWkw09vPf66yml+UbD/Brwuu71Q4clbBUR25M7QG7TrboA+GfgB8BV5I6xA4EnkTvuPhERO6WU3lwo7rrueVvg/cDTyMEw3km+1ffW5E7Gp5Dbem+NiHNTSl8u1GsT4CvAoPPlSuBd5A67S4C9yJF879Ol+0OXrqWT7pSuHICPkTtp6Lb5i166yxvKbBIR9yRHHh6MtjoD+AS5w/R68lx6fwc8kJnPaJg3AP/QvT4b+CDwI+C35OP/bPJ5uS/w1Yi435i/M+M8BzYDvkmOcgz5fHwf+Rw4vyvrAPK5eSwzo2mLumP9bWBr8nnynq6sc4HtyJ3Rh3Svj42IB6SUTu0V8U/kEZf3Bx5J/lze2St/Nfk8Gow4fV5K6exhdZIkadlIKfnw4cOHDx8+fAx9kDsCUvd4+gTKP6FX/m6Ned/by/usBdRhNbC2K2ctsGoc9Sd3anyn9/5XgZtUynkCudMoAeuBBxfSnNwrK5E7iFYXtvmqXpqzgCiU9bpemjOB2xbSBHl0U3+b3yuk26/3/lEL+ZyB3XrpThjh+F8wSF94byvgol557wM2L6RbQ45Qm4btB3nU3eD9D1TKWjXr2O7feC4e1cu734TPgSNnpblDpU6P652btWO9NbljL5FHte5ZKevewO8G34fC+7uSOzwTcC1wl957h/bq8P7W77kPHz58+PCxlB/e7ipJkkaxY+/1pXMljjynWao8Thhz3fq3E/56voWkPNLpN93iavKtfuPwRPIILcidSU9MKf2uUocvMHO7XwBHFJL1I+SeRR5JdINRWimlRJ5n75Ju1V7A7ftpujnGXtotrgeemlL6aaFOCfg/5I7GgZGCciwRLyTfbgl5VN8LU2HOwpRveX4meWRcUXfr5eDzOR34m0pZ64HDgcHcfYfNt/IV4zoHdgZe1C2uB56SUjqrtMGUR+EdPke9ngvs3r0+KFUCy6SUvk++jRrg0RFxp1nvXwQ8q1vclDyydLOIuD0zQV7OIN82LknSimEnnSRJGsVSniJj297r66qpRtMPHLDQsgae1Xv9Tymlucp9GzO3ld4zImbfftkPjvDe2Z0z/5vdur9cAAAKjElEQVQob+cLvVX3nJXkycAW3euPpZR+VKtQ1+Hzztr7S9zTe68PTUNuh045IvB7au8D+5DnEgT412GfZVfW17vFe0TE1iPWdxTjPAc27V5/NKV0xhzb/ZdZ257t2d3z91NKp8xR1nG91w+a/WbXKTg45+5M7hz9MHm+xD+R56G7eo5tSJK0rCzlBrckSVo6ftV7XY242HNXbjh31U7A8WOt0Yzf9F7Pe/RbNz/bdt3iWuCKhVSqKzOAh/RWfW2uPCml30fEd8i3F0Kem2uDaKqdn8xR3IW917MDGPTn3TthrnqRgyQsK13Qhjt3i9eR5wWcy4VD3tuv9/q9EfHehurchjz6btwWcg70z81vz7WhlNIfIuL35HkTbyAibgYMRsTdO+aItDvLnpX1hwEPIAevOKS3/nkppWV3PkqSNBdH0kmSpFH8rPd69micDaSUzkkpnT14kOenmpR+B+JdqqnmdjdmOhbPHDbiqsGNmQlW8AdGD47Q74DYsZoqz9c1TH9728x6b5fe61Em3l+O0TN36r0+tzbibJZhIx13XUBdJtXuXsg5sHPv9ajBF2rHZ5fK+lEUj01K6VpyQJX+7b1fSCl9cgHbkiRpyXIknSRJGsU3gFd3rx8dES8ZscNjMZzEzNxUDwdeOc9y+lFiv7egGs3YbJ75Nu+9Xkjn2LDRTP0Omz8uYBuTtnruJFVb9l6PemvksOilW/Rev4T8vRjV+Q1px2nYOTDO49M/Nt8hRwUe1W+GvPdAclCPgYdGxG4ppQsaypckaVmwk06SJI3iVODn5JE3u5PnsvrUVGs04xvkkTZryHN/3bubmH5kEbGGG3YqfGxMdbucPCH/KvLceTuRj+Nc7tx7PUr6+bicfAsm3fNct02OK5DGKK7svZ59i2bJppX1/VFke0ZEdPPrDbPBrZyV8rbtRokuZ5czc6vpnsCPhyWOiFXk0aG1sgZ2HcexiYi9gHd1i4Prz43IgST26+b9kyRpxfB2V0mSNKdu1NzrequOjIhbTqs+fSmlP5InlB94Zze/XItDmbn179vkTslx1G0t8J+9VfvOlScibgTs3S2uA747jroU9PdxvxHS33UM2+zfQjxshFw/+u1OXedQUUTchPotwef3ytqWG3Z+1tx7yHs/HDHdcvGD3usHVlPNuCv1DtHzgd93r3ePiO0XUrGI2AL4DHm039XAg5kJNnF/csRhSZJWFDvpJEnSqD7MTCfFjsBXI2KnIekX06uZCfRwb+Bdwzp2+iLi0czcynsV8KIRRlu1+Hjv9ZsiYqs50r+emYi1x6WURp3HrtUxvdcviIjq3HcRsZobTtw/X/35A6sj1rrOzUG02a2Bhw0p8zn9hS5YR7+cfnTT1w+rXBf84Fn9VbOSfIuZ8+xR3bmznH2+9/oFEbFzNWX2sv5C4Vh/sff2v/Tfn4e3M9Op+n9TSueRR7sOjv+rIuK+CyhfkqQlx046SZI0kpTS9cDjgUu6VXsDP4yIp81j5NpYpZR+CRxAHnkG+Y/5z87R8bRpRBxK7ljYhBzR9eCU0ly3fbb6EDO3ku5KvlWvFB1zVUS8CPibbtW15BF+k/IdZkbpbQF8sTQ6suug+yAzkTsX4uLe6wfNkfa43ut/iojZQQ+IiEcBr521es2s5TcyE+zgMRHx2m6fZpd1426bN+qtvsF53Y3aPKK36sMRMXQ/pv3dGCaldCLw/7rFzYFjImKDABCRHQY8ddZbs/ftdcwEsjgAODwiqvMyduVucHwi4gDged3ifwD/3NX358BLu/WrgaO7kaeSJK0IS7bRIEmSlp6U0i8i4j7k29DuC+xAHin22oj4CnAi8EtyJNMbk+cTuwN5DrtR/XVE/HaEdF9KKZ3bq9s3uj/uP0oOivAk4OER8Tng68CF5LnrdgQeQI4auVuX/VfA01NKxzfUcyQppesj4nHAyeTj8RfAvSLi3eRbTv9Ejix7IDO3w64DnpJS+um469OrV4qIZ3R1uBlwL+DUiDiS3DGylhwt9yXAHYHjgT/rso8aZGC2zwMv714fHhG/B04jfyanpZQu6aV9J/Ai8md5H+C7EfFe4HTycTwAeDr5fDuRfPwgdzb9bwTSlNJ5EfF8Zm6JfiXw4Ig4ljwyNLry/747Dh8gn6/bcsMAHgNHkM/9RwPbAydExDHkc+xs4Jpuf25Jvo34gRGx6xIKtDLbM8m3ve5IPgdOi4j3k8+Li4HbkzvM7k/uoP8OcFCXd3N6kVdTSv8TEX8JHE0+rq8CntqVdxbwC/IxvgV5lNz+wGH0RvRFxO7kzwDy3ITPmHXsPgQ8hvwPg92Bd5PPA0mSlj076SRJUpOU0sXd6KG/I98CuRO5s+tvu8dcLp3j/VeMWJXzgHP7K1JKX4iIH5NvlXsU+VbJv+weJZeTO13ek1JaSBTVobqOonuQOy8eQO7Aqd16eR7wzJTSyZOqT69eF0bEA4AvAbcjd9S8tZD04+QRaWd1y78rpBllez+IiM+SO8G2Jo/QG3gScGwv7a8jYn/gs+QOs7sB/zqryHPInTXP6K3bYORWSukjEXFdl38rcofT/QtV/AD5vB50KpfKWhsRTySfN39HvjPlSd2j5i7kzsglp+t4vx/5tuC7kQOElL6DPyF3jD63t24zZm4/HZT3qa7z9SPkTvzbAG8eUoV96TrpImJT4NPM3O790u421375KSKeQ+5Q3BF4WkR8LaX0yRF2V5KkJc1OOkmS1Ky79fXIiHgXuTPsgcD9yJ1PNyH/8f5r8iinX5FH4HwfOD6ldHGx0PHV7afAoyPibuSOk4eSOxFvzoZtnwvIHQlPiIjTgJ9OKmJkSumiiHggeX61g8iddTuSR32d1T2+BRy7mFErU0rnRMRdyB2ZTyGPcLoRedTTqcD7U0rHR8RDe9ku2bCkkT2V3GG1P3mU5XXkz2GD0ZMppW9GxG3Io98eRh459RtyR+bRwCdSSldHxB972Yq3V6aUjo6I48mj8x4D7EE+H35OHuX4wZTSSQC9qdRKI+lIKV0HHBIR/wL8Ffkc24M8evS35ONzCXAm+RbaM+Y4JlOVUrqg60R+Kvnc3AfYjtyhfg75WH86pXTVrGnmasfna92IuKeTRxzuQx5Bt5aZY3MBeXTm13tZX0fufAP4N2ZG1M0u//KIeCbwzW7VeyPi5JTSBaPvtSRJS0+Md15kSZKkpambh+x25BFUjyLfurn1rGTXkAMWnA4c2s1BJiAiDgHe1i3un1I6Zlh6SZIktbGTTpIkbZS6Ce33JXfYPYrcgTdwAbDHmKO8LlsRsTl5JORdyPPR7ZJS+s10ayVJkrSyeLurJEnaKKWUriXfbnc88NKI2IPcWfdo4NSNpYMuIu4AnFMLbNB10L2b3EEH+bZQO+gkSZLGzJF0kiRJG6muA+6X5CianyYHB7iAPK/azuSori8iBwcB+G9g35TSFRsUJkmSpAWxk06SJGkjFRH3BL7NhnPzlRwDvDClNFd0XkmSJM2DnXSSJEkbsYjYGjgQeCywS/fYBrgIuBD4MTm650ZzC7AkSdI02EknSZIkSZIkTdmqaVdAkiRJkiRJ2tjZSSdJkiRJkiRNmZ10kiRJkiRJ0pT9f0ohvBkY5MMZAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "image/png": { "height": 333, "width": 628 } }, "output_type": "display_data" } ], "source": [ "grid = coeffs.expand(grid='GLQ')\n", "fig, ax = grid.plot()" ] }, { "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." ] }, { "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**2) * grid.weights[np.newaxis,:].T).sum() * (2. * np.pi / grid.nlon)\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:" ] }, { "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.nlat)\n", "\n", "N = ((grid_dh.data**2) * weights[np.newaxis,:].T).sum() * 2. * np.sqrt(2.) * np.pi / grid_dh.nlon\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.1" } }, "nbformat": 4, "nbformat_minor": 1 }