{ "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 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": "iVBORw0KGgoAAAANSUhEUgAABOkAAAKaCAYAAABrxw0nAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAjEAAAIxABZy/d1gAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdebSk213e92dXnTrdfa9mEEJCTAaB41ggCQI2Cxk8w0oAM4MWjiGAvYxXjB0sWCiGRRKDGRzsxNgyYVCwCMbMMoMJgcWSsRnC6ETYzAHMJMQgdO/t7nNq2PmjqtHRrec5/f5OvXXq9OnvZ627jrR7n/3ueusdfvV29X5a710AAAAAAAAADmdy6AkAAAAAAAAADzse0gEAAAAAAAAHxkM6AAAAAAAA4MB4SHcNtNYea609duh5APvA8Y3rjmMcDyOOe1xnHN+47jjG8TC6rOOeh3QAAAAAAADAgfGQDgAAAAAAADgwHtJdUGvtpa2117TWfru1dqe19rrW2itaazcPPTcAAAAAAAA8WHhIdwGttZdLeq2kD5W0kvQLkv5zSZ8v6Qdba08/4PQAAAAAAADwgOEhXVFr7cMkfbGkLullkt6u9/5CSX9C0hOSXizpKw83QwAAAAAAADxoeEhX0FqbSfpHm//7it77v+i9LyWp9/4jkj5x82cf1Vp7yQGmCAAAAAAAgAcQD+lq/oKkd5J0V9I/NX/+LZJev/nfL7ukOQEAAAAAAOABx0O6mg/Z/PzO3vtjT/7D3vtK0nds/u/7XdqsAAAAAAAA8EA7OvQEHjD3/gnrj53T5//d/Hy7fU2itTaV9IKzTZv2Z+9rm8ABcXzjuuMYx8OI4x7XGcc3rjuOcTyM7h33byPpWWfaf/7eMmhj4CFdzTttfr7+nD6/u/n5vNbaZPPtusFaa1vf0HPdJD1q2n+7si3gAcPxjeuOYxwPI457XGcc37juOMbxMHry86AnWmv9vF/ovT916OD8c9eap21+DnlIdySJv1kAAAAAAADAffFNupqFpBuSpuf0ObtPn6huYMgT1tbaH5X0H6tjAwAAAAAAYDTv3Xv/mbEG4yFdza9JendJzzmnz1ttfr6p9/74nubxe2f/z3Nf8XJNH3X/+vXh0Q49gYs49wuxl+gKzOP8Lwcf2FWeGw7rCl94+lWY21WYg3R15mFweYFzhQ/Ztatw4F6FOYj6BQ+gK3yBuRK1i3R19tFVmcdAD9MlZ/nEE/rNL/iSs02/l/peBA/pau49pHubc/q89ebnb+xxHm9xDkwffVTTpz5lj5u7+q50kZRclTlfgXlc6ffvKs8Nh3WFi6crUehehTlIV2cexpV4n3DlXOl7onQ17otXYQ664u/VVZ4bDucK33euzD2ReVzIlXn/DmPUKy5r0tX8p83P9zmnz4s3P391z3MBAAAAAADANcFDupp/s/n5Qa21rbXjWmtHkj5483+/59JmBQAAAAAAgAcaD+lqvk3SY1qHR3ya+fNPkPRMSSeSvv4S5wUAAAAAAIAHGGvSFfTef7+19vckfZGkL2it/ZKkb+69r1pr7y/plZuuX9p7//VLm9hRX//3INjTNPsY444yRuEf44eupfVNLvlt33ntlerv77i4wShrxVz2qfWAnMpX0r7WwtjTuOXDu3BA267F7V362iI7bm+f+9P//m6/fpCxH+71Yi5un9dlM/Yom9tnTWOaH8T77aXXW4WL1F7XutvX2NQvF/OA3WulkeoDc5DH4/6qzLmisL197c/6GAf+/UONfRF7fvbCQ7q6f6D1unMfJ+kbJP1Wa+13JP3xzZ9/s6TPPdDcAAAAAAAA8ADin7sW9d5Xkl4m6VMk/Zikp0l6vqR/J+kTJX10731xsAkCAAAAAADggcM36S6g994lfdXmPwAAAAAAAGAnfJMOAAAAAAAAODAe0gEAAAAAAAAHxj93vQb6bKU+W13wl8edy5vHveQIlj0mdDX3WuIYhZi20G53XdiftnWMcJ/inH3f7dmNM65vLh1xezperny63WUnr12XJK199a1qLlLRd/XXkTSu/4PKS9nX/ozjujmX010HtgU9nfBX5Xix27vs+Mw9jXvpqdtXuKYJ/W3tUt1eHz6R2LVSC6Qp7+s+XJhzaXsj1Dp7q2mq87jsBNyK61LTXHZy6Qh9S6dq3N6ONY1US72vjFtRmXP1HrxjnZLntqeL0XWqae7zWi787GUgvkkHAAAAAAAAHBgP6QAAAAAAAIAD4yEdAAAAAAAAcGA8pAMAAAAAAAAOjId0AAAAAAAAwIGR7noNtOOl2o3l+Z32mExWCP+qzaOSjrWvcVV9fcPHjclrNizGD9JN9EzcE3Eew/vunNhaDUM0+2KvCbGFY2ucBLnhQ+wtwGiMcfd0eSlftnZMoKpvb3vn5QQy01ROMRvYlqTtTQo5bXtM7uqTwvYK+6KUwlran8ULWiXpbYw0tR0vGun4HEVlbnuqX0r39jhIOugqY4wwrguW3+M9qpR6byaS5tZW+7qZpPZwf3d1WOVULe/P3cbIddGOKbXnzWPHvnsNmt517H3e50a5tu+4vbDzfd/h25L8fSO+ZvNVoXhtiPti+w/yPt4ePFY/hX0/Sgp9+NqUrV9GOC5Ktcdl10XVsV3XvSUxv3kO7eQ+z152xDfpAAAAAAAAgAPjIR0AAAAAAABwYDykAwAAAAAAAA6Mh3QAAAAAAADAgREccQ0c3VhqemNxsV8urGDaK31H2d7wwXtabNhsL4+bghzcBguhD6FvDxNpdlVSv7nS3CoqYQk26KK2ELJdpDn0TwtL7zy3sL1RFnoOdh2jvPDyPhdqdnY8FEsBCtUxCuPmeZjrS/qrr0rQgQtQSPNI26v8FVw6/+zr232x4dL7WllMufrXju61jPD6Whxj+ALQrdTXt7vtlU6HMVZ239fq8CPc5/ZZ09jbewynqoQwhO25+2Iqduy9MtU0YXtukHi/NWOEvukt2TWcIQZSjHB/t/VLsfaw8yu9vuF94zwuufYY5dJwhWuavdYeu46Rbho7hj7E/ume6I7bSg0VttfS3GyC1x7vc24exbqvFDhlt1foG/qXaozK6zhv7KHbK/z6envjXejuG9q5I75JBwAAAAAAABwYD+kAAAAAAACAA+MhHQAAAAAAAHBgPKQDAAAAAAAADoyHdAAAAAAAAMCBke56Ddy8MdfRzfmFfreSaZOS0EoJaS4QbOWfFcdgMpfYmhLIbLprSkcL2zNjxzRZ155StyrpZoXUtJiONkagjZlHTiszbSEIJyavuTHSvrDJuqFvJaVtjL5XOU12j8ZIJiuNW0pQ3XHcNEbhr75S3zYJ16jK9kzCVp+m7fl2e7jEnTTCBcYly1WSwlLf2G6GLYyR+sZ229df0Fzi2SSNG/ana4+hfrbv8Pd0hEzxUeyrpqnUI6l9FesU05bqIj+NneuUVP/YBPkwRjoKXJK9TXxd/4Fvr9xMCmm5pVqnUL+U6pS0vdS3MG4l9XWMuqi8nwtj7DzuCPaXtrq/MWyNcOl1Uehs+4Z0+1C/2DHS+eD6hmGT2r4vpNBXUlj3Wae4uqhSY5TrlO22SaEuSipzjmPc588XJxd79jIU36QDAAAAAAAADoyHdAAAAAAAAMCB8ZAOAAAAAAAAODAe0gEAAAAAAAAHRnDENfC0myea3Uorap5vVQl9KIyRFlNemsWQ47rEYeHkZWF7bqHm5dKPGxd1dnOOfc3clmEfx3bTlhZdtWEC6Z2qjBFGcIsNx8WUXV8/hziGaR8nqMK37y2oIhzkdn8WFnqOxli8eQQ7L7JcWTQ59E/BCt2uVhv6FtpLfcNlO4dBbDfF0IepWbg+HBjxsCj8NZ47tuK4lUWW03s9NX1dm85ZyHi6fbKl/ekWMk7jpoCHqdleXCDZ/X7oOy0s1DwpLKYc+9pWL42xq0rtIvljMY1hQx+KwVlLV3sUwiBi/VOoX1JwhBs7BmeF7dn+4X7rap20OHwMlLB9B3c9p4ANY5vXMqnUL8Ugh13rl2pwhG2PAWVmcfhinVIKvnoQgyN2rHVynRLa3WeCMQKuKjVNoX6Jr8ONYWqXTe/hzWluIwSw2YOusD/jQRvqF9feCrVODJwK9YTr72qX1DcFM0xDu6uB0kdeVwOl7V1G/TI/OSmMVsc36QAAAAAAAIAD4yEdAAAAAAAAcGA8pAMAAAAAAAAOjId0AAAAAAAAwIHxkA4AAAAAAAA4MNJdr4Fn3XxCN269OfFkFWN7hlmF7JOUbubaU4qZG3uRUswK6a7zpY/ycQlri9g3paltp7ssQ8qXTU1zKZKSVpVn5CHly70lLWwvcoNUElRTIthie9xq2urEpbsuwhiubyG5bb29QorZnhLWXJJaHKOYSjVKGp4duDBsJZm1mKTmEs9y2ur2C6wks0o+mcymrqW+IYFslZLJ7PZ8X2eSru0xMtlNwne1oWnVJDz3WkIymU88C0mpIQnNpq2GMVzflHiWUsyO3BghYe3ItLu288ZwyWTTcDFyfVOyme1bvBhN4oV1W6WmqdQvpZomzCGN4eqJVOu49lSPLELaqquLUhKsa1+lvulcdYmtMR7bNYaBU8rs8DBLn5peqGkkn+Qa6xdTk9Rrj+F9fa0T6oY9JcTmvsNvELEeqfRNLrt+KYzr009D3TBC/VJJbF2ZZNVqCr2tU8LTBnsYxeMi7NCj4anwdoj4pg6vi3q8wAxPoU/tkyOTaJr6mlpgWqh/JF9PHE39hcQly89S35T6Wqh1XJ1RqX9Se7V+uefk5p0L/d5QfJMOAAAAAAAAODAe0gEAAAAAAAAHxkM6AAAAAAAA4MB4SAcAAAAAAAAcGA/pAAAAAAAAgAMj3fUaeKubT+jWzRDhdB8uNS0mnoW4HJcqtgjxhAuzvZh4FiKFTk066zxsz6W+zld+X6XU19PFdvtyGua8MEmwCtFIIcbMpr6m9KHhAXklMcDKpXyF+DeXKuYSzKSc2DpxqWkpHc3s+0o6mlRLWKskwTbTN/XPKXTmD4qJdRcMMLq4SmqaSwQrJp65nefSytIYOa0sXBNtYqvfyStzt03jpjtzN+9rKTUtiInQZn/uGB6+GTe1D09Cc0muE5fyJmkS0saOTGpaSjxziWUurVWSjsP23BjH4aLoEstS36PYblLTYrrrdvuslAS7p5tRUUqCdXXNPPY1dUromxNbty8Op6FOcUmuqW9Ksnd10WLi5zafmJompSiHG9KqbY+xCledbi+2fnOqptPbDZph4715eP0S65QdU+ilVOuE61mlLoo1UKF+Me2u/pFUqklikn3l48wVqWns5aGSQh/Os0pNMkZiq0sNjemuKZ3+yGwv3R5cXVStacw1I+1Pc9ka5xCqJNYXUlwlaWL6T48KNU1IPz1OY5g6ZRbGcDXNLFx0Ul10ZE74lO7q6pdU/6Q0WZfkmpNgz69r7ty9e+6f74pv0gEAAAAAAAAHxkM6AAAAAAAA4MB4SAcAAAAAAAAcGA/pAAAAAAAAgAMjOOIaeNvjN+nRGyfn9kmhD84yrD46D6uE+gWZw+LGZvX0FPpwEtrdGG7R5NT3ZOkP+5OFb3eLbrowibXhp1QPgRJuMd3uQgMkv8hycd1lu15mXEB4e/DKYsOVgAhJmsxd37AgrBkjjhvG8MERI/QN27NBHOG9tgs9pxVv4xrgw8MnRmGOxRQGYfvG0IdKGEToa07VGDIRF0gePoZbDLuHoIN0Tq1m22OnZW3tbk6LiKfFqd30KsdLWgw7pqOYtrCgvVvguhIQIUkzs3ByWkz5+Gj7YpIWQr4x9Rce1//m1Fzk5EMiboQLWlo42fWfhVXZ3cLJlb7TFDKxpwtMDrKqBEfsXtOcuIuAfHBE6utCIu4uZ75vqHVcXZP6Ts396NStqH4OdyT2sI9s/ZJCgdKC77Z7Wpl/YJvOC5Qo9K3UHvNCGEShXkpBDmkerp5IY7h5xDol1oMu2Gt4+EQKmdhb/ZIOrUr9ksawoQ+hcwqDcOEMMeDKhUEMD8Ny21q3+7m59zqNYccN7TFTxtbR4frijsNU/4TN2WMjBkeY7cUwLN/uQiJSTWPrlFTTpPrFjBFrGnOBSeO6vpJ0w9QvMXxix5om9U/1ix33zJHxxPH5z152xTfpAAAAAAAAgAPjIR0AAAAAAABwYDykAwAAAAAAAA6Mh3QAAAAAAADAgfGQDgAAAAAAADgw0l2vgecc/4Geenx8od91KWbL8Ow2pZu59nmI/TkxkYo53dWPcWe1/VpTYusdk5B2u/l9FdPp2vYYLcZqOn5uq5Q+ZFKeekg4tOmJMQJpuBi+6DYXEiNdOE8I7LEprut2k0A2QhLsNCasmSSmQkJsW4RjKKWmuaSpmATr0l2HJ8FGKTVtDOZYzIln5rqT+sbE1u0xYoqZSRubpL4mVVWSJgtz/fShjDb1Ll0DcqyY+4PhaXNpc6WEZt+1dhilvx50SWjh2jcx7dOQjnYUE1tNgqpJNpOkW7PtC0xKPHvk6NS2uyRXl2wmSbem22OkdNfUftNcFFMSmmvPqWnbB8Y0ZPKlhLVdudpFyvXLyqS+lmqa0Pfuyp/wrn6JNc1yuyZJqb8p9fX2YnuMo4nfXqVCiAHiZv+nFHpXv8SaJqVc2oF9sw1fLFzjUntOizdtqcaItY4bd3hia6ppYuK86V9Jp099Uw1k02DT+7fcHiMn2e+pfinW0TaxNY3havyQfppqHZci7+ofKSTBhmR519e8HZJygqq9DMT3yYwRk1J9s5tGCqt2x1H5CLIXmHQ9q9Q0fke7JFeX4ipJN2YumdXfx28e+YvRTVPXpJrmlqtpQl10a+LHsCn0KbHefCCbhQtzpda5aP3y2LF/TWPhm3QAAAAAAADAgfGQDgAAAAAAADgwHtIBAAAAAAAAB8ZDOgAAAAAAAODACI64Bp579EY9fXb+W7k0iyZL0so8p01908LJp3aRZT+fu3170eOTsPDybRMQIUm3VtsLVd6Z+DGOJjdMm19M8sgsvCylhSN9X7eY8nLl9+dRWijW9E+LLIe3qsYu8hr6urU14+K/pq0QwrBuN21h4eXpqQuZCMERpu96fmYx+tOwELJZOHkyT31Du1lkOS28rJXpmwIi0iK9rj0s9l1SWex74g9au8hyCHJICyS7sVeh78SFTIRzcrUI108TKJEWU16a9rwm7fAwiBTE0VwOR1hMOayZ6wffZ9CInXO4NkxNeIFpk6RZWDjZhUS4gAhJumUWWU6LKT9qQh8kv8jyo0cntu8jZpFlFwQhSTdC+8223Z7G8MER/oJ9bPpOwkGUApl2ValppFqd4mqdFBDhahrJ1zVpjNvmuHhisV27SNIsBlyZ+9EIoR3p9rAydUoKwnHtPYRIjcLWNOFCWQnNCXN2p0kl4ErywQ8pfKLSN7e72mN4GMSkUKekMZTCJ2ydUryO7FrXpJomhUG4kKxQv7g6JdU0qX1igrZWs3SvNJ9Ljv249hhPNX58Swr1oLlGpbAMV9NIklxdE2uagW1VYW4uJMKFXkk+IEKSpiZQwoVeST4k4pGZr0dS/eLaXe0iSU+ZbtcvqW+qPVytU6lpKnXKur8Jjkj3VXMgne37ByHAYyx8kw4AAAAAAAA4MB7SAQAAAAAAAAfGQzoAAAAAAADgwHhIBwAAAAAAABwYD+kAAAAAAACAAyPd9Rp4m+ljeub0zfE2y0KyzsokpJ3aqJyc7uoS0lySmuST0G6vfIrZjdDuUl9nzfd16WbTYuLZyqSCubbUvgzpTCn1tZkknxhr5JKRwutrlcSlQppTSniy6a4pHS20u4Q0l2yW2ichxTWlsLok1xb6unSzNvcvJLa7xLNF2Bk23TXs/JTAmfpXxnBS4plL6Qp9m0s8MwlmktSO/PWlT7fbpyk1bbbddzUL24sRh6b/CJHLKYjQpZ6lpLC+NIMUkgylkD5buXwWLjlpgy4dTfKv26WgSeckobl0V5PiKvnEs6eGZNaU2OqS0FyymSQ9Uuj76MRvz6e7+jFm2t5HOR1te79Nw4ExCQdX6u+4msbVLqmvVKtT5qYGuhvS5lO66xOmfkmJ9bOlS9YNKfQmsV7aPck11TTL0L5Ybu+jZUjutudwdb6Va0lhczGhspDuOjEppeUke1O/pFrH1TWTkEI/TYnzpj2m07v6JaS7NnMsr/tvt8f6xSW5pnok3ZvHqF+cUJPYMVJfc57EmqZQvzRzTko+4T7VNDb1NVxr4ynsdmfYx27o6ueEyueSMZJc7WeseIExTaFOSbWcq1+OQ2L9TVO/3Axpqynd9Smm/VFTj0ihpgl9K/XLI4WaplKnrNvN/jT1j5Trl3uOwvswFr5JBwAAAAAAABwYD+kAAAAAAACAA+MhHQAAAAAAAHBgPKQDAAAAAAAADoyHdAAAAAAAAMCBke56DTxjMtezJucnjKR8EpeEltK85uGZ7l2TkJYSz+6utttzMsvw9kpSXBJTzFbbr2+RUuFW23M7DYlLKYlw4pKfUoKjTVGyXbNK8JrrW0iMdClo0nnpribFrJCOVklxlXxC2uR0eGJrSnHV3KcMNdfuks0kn/oaktR6TH0tpKZVpHQ0l0rsEl/l012bSWuVJIUkNNs+C7c5s58n5lyXpB7a44XVMum1hcQzSXKXnT5C4lnlfB9FSqB2SWhhiIm5fqakrVm41s7MffPYRUpLumUS0m5NQzpaSDd76vSu6bvdJtUSz1K666Nte34uHU3y99Wb4cI8MQfGcTi49vU3wenUOw0nz8ocSa52kXyS/V2zLyXpiR4SW02y3HEhbW5aSeKWT2ddhH2xMKnUpxO/L2ah3Z1ri2VIgTRpuemWkdLp64XNkwcOo+4pyT6V5OHyYttjOr2pU6Yn4fwrJM6X0ulTTZPS6Sv1i6uBYoprMeF+qEpNI6mZJNdS/TIPV8pQvzRTS/dUX7sk+zHu7WFfuMtOqmlWU5NYH15HX4Xa0R0ClfO9ui8Kn7vcZ7eU4po+E7r2Wah1XOrrzenwmkbySa77qmlS+yMtpLtOTHptqGlS/TIz1cM03Ajcs4Wzh7KrQ8fEN+kAAAAAAACAA+MhHQAAAAAAAHBgPKQDAAAAAAAADoyHdAAAAAAAAMCBERxxDTx90vQMs+j6WcuwMqZb8vA0LMQ6D+1ucca7YTXz2zILdYe+k7A09LSwWrsLxnALSEt+sWjJB0echoXkT9r2KXWUAiIqi6fvuGZyWQjRqCy6Wll4OS2y7Nb9TCETbiHdFDLRQqBEaTHlU7MYayUgIvTvi9DXtVeDI9yCzPsMjnBhECkswS2yfBQWTV6ExZtnJrAmrZDs2otrwLrXMjELIUtSnw4/PsOlSCaXRi0uprw9duobF06uhMpUpOuZGbzFgJ3tvtPigsxukeXjcDHywRF+weJHJj5kwLWnxZSfOtlefDn1jYESNjjCvz7XPgvvk4uFmoXzerrrgv9BqmlSneLeqZuh9nCBErNUp4TtuQWnKwFXq7DS+jK0z831MwVcnSy3r6vuXFj3TYuZD18Q3Z7DYdH5MQ6X0uLwpfolvD7XNyx+nxbFd/eCOIbruxhe00hSO3G1TqF+iTWNvybamiTVOqau6alOSUFb+wqOSJ+1THCEYnDE9vnXUq0TXofbH6mvHTcFHZj2SQj7cjWNJLXldv8YDmfaUxZI+vxQC79zr2VfCVm1z3PpM6H7DHkUdoYLvkp1SqxfTCDWIyE4woVEPG1yx48b6pSnmVonhUE8YoMjUkCE35+urjkOb8rE3JDO1jSLcM8bC9+kAwAAAAAAAA6Mh3QAAAAAAADAgfGQDgAAAAAAADgwHtIBAAAAAAAAB8ZDOgAAAAAAAODASHe9Bm5NZnpkElKENpYh9Wdl4gxvhISZkxC5c2rGniqkNrnHwiGxZxkShVZmkNOUYta2c+jmzc9tFhJijkz8UEpsdSmCKbFnEsZorn+KKrr01FczhUpqWkxSS6lpLqEyJUoNT0eL7Sb1NaWjuXSzSoqrJPVTk2AUEs9s6mtMdw1zdvsuRWmNwSX4uRQ0ySahtTC33kMS2uCJhb7hmqNFSH4y53AP6a4u9TWmo8XE1mFtsT2df8X2vSkkobnrZEpxTdfrI3PNvzH155/re9MkjZ3XfsO0p8SzmyYJtpLiKkmPmHvdzXDA3DT7+WZI4JyYN2Umfw+e7imaPNU083AfX5n+d0PxMTEp9JNqAqDZdS5tPrWfhmtcSqE/Mf3vmPpH8sf4XZP4KtVqHVu7KJzDl1y7lK9xrtYpXGtL12UV0+lNkuskpbumdlMjlOqXlOJaSKfvaXsu3TUl1u+rfknpw6lGcEmu07DvTR3WU13rt2bb06HsUuhdmyQ1U5+V62hbt4e52TqllkJfSnPel0piffEzofsMWal13GdYKX/mvWHrBn+u3nQJ8qH+Sen0LsnVpbhK0qNmHjfCvk+JrTfMuZ1S6CfmRn62prmV0p5HwjfpAAAAAAAAgAPjIR0AAAAAAABwYDykAwAAAAAAAA6Mh3QAAAAAAADAgREccQ3c0Ew32n2CI8ICySuzuua8h0Xn40K/22O7cSVpZhZ5TYtXzszizev+2wtHHqcxTN9pWME0zWNqFu1Mi0i7BT5jcERhQfQrbYTgiMqCsOFtCn3TgrdhQV+3GPkqLbxs2tPixjHIwbSnvnYx5eF9Y/+wAHtJXN1/+NilhZBDezcLwraweLN9X+P7F2Zixmir8HdfezofKmM8aJeW81TWnS9dr9O9y1x48rj+ODp296O0PXP/O473xLS97XYXELEew4RBhMXTZ6becAssr9v9BqdpYXZjaeqGdF9NC0DPzb5bFk4eF7K1Hjfse1NHpZrGHQPuWJHyseWOxVTTuL6xTimswH7ZOVb7tLdr7b6CtkKdEu9d7l6XxjDtMcghBV+5kIhYp7jthVonBUfsWtfEYAX/WcttLZ4Pbuz0gWCZAixMrRPCBLqrU8L+cXVtC3OIYSzmrRojYK7cfkXF4IhwCJQ+V5qd4T7DSud95nWf0/157T57p77xs74bI9xXZ2YfVQIi1mNvn8OuppF8/XK2drmx5++68U06AAAAAAAA4MB4SAcAAAAAAAAcGA/pAAAAAAAAgAPjIR0AAAAAAABwYDykA+3uWj0AACAASURBVAAAAAAAAA6MdNeHREpSW5kEsmlISpmHBJ3Kk16XMjMdI10y8Kk4IQ0qjVFJQrvK8YnXKXrNcYlnVbsei+H3exrXzDn3dWlsYXtj7IuKuN/MnLs/ELtpLx+yLuktJdb5SVS3OFhKU7s23Pt6ydfDq3L9TYmttm+4H7n2dO9K2yvdm3e8QYyR4prYMUKqY+XunpNgXepvGmN4wm/lvU4qx9Y+Xfq5ZuO/r3dRs9d7RmVsdw9NqaqlKYQ5uLH3leJalGorm9gZzmv3um1ib+grjXBsjLHfLrvOvGxX+PJS+QxaSehejzE8TX0M7p6WEmmdaoWRnnH4vof9LhvfpAMAAAAAAAAOjId0AAAAAAAAwIHxkA4AAAAAAAA4MB7SAQAAAAAAAAfGQzoAAAAAAADgwEh3vQaWWml5n5SrVUh3WZkstLlJfJWkVUgDmpv2ZQhmcfNchgidVffPkJemfRmeN8/7dPi4YYyVmd8q7O/U7vsO7jqOy95eJRkp9LW7M407MX+Q/hrC9ZUkl/oTkoC6GaNN/QbbIoxh+rdl6rt9LKeEribTV5I/tXdPaVNIQGpuP6e+bj9PwhuYEpfcPopjDH+v43FUOV5s3zDuGC47mcylcV3yNSdelws7I/VN9wffN4zh7l2F+5y7n0nSsvl79tzMYxbelKVNNA0Jcu66E+aQ6oaUBmvHMPNwtYskLcP23Otzbev2bW5frvuGhFjzXlXe69i3sN8qNU2S+lZqnVG4tyomAO5pbpd8vbb3DOnS7132HurutZK0DPWEGcPe8yV1c39vYXMxyX7X9NlUp6TX7Wqd0Ne+7lQ7pvfE9Hc16WaQYW3SSHX0dlO8XIzwOWFvLr1+qfQdfl2u1jTuM3KsPcwY6d4VP3ubYzE903D3bJfGLuUk9Lm5w8/CsTU3Q5ytXZZjfH46B9+kAwAAAAAAAA6Mh3QAAAAAAADAgfGQDgAAAAAAADgwHtIBAAAAAAAAB0ZwxDVw0hc6ObNYclpQ2ZmbxVXTYsouIEKS7prmu2GByNt9+5C7a9rW7bMwDzPGKvXdXuwy9V2sQvjEanuMRVwYc3v1yTRuDwtj2vbLXqQ5qQQ5FPrmBZJN36k/Dt2iuWkh3dh+ZN6r8P61pVlINS2afOSPcTeLuBS2C59IoQhLv4i7C7bo4by2q9imhYIDv0ByWEzZtadFmtP+dO1HIUTDtYe+aT+74yUfc6YxnQ/hbXXtsW/hMnJVLi9OOjxdc7qmpkWWF+banhYsdn1PVv44PAn3mJttvtWW7nOzvthqm4bF0Gch7MmZantcyYcz3Ay1hOubgiCme1rtO9UpMeDK9L8bw7C22+6GhbNz+/b7mt5r154W6k7HljsW3TEr+WM89U3njjvX0r0rncNXQeWaGPtWrssxyKEyhtn3rnaR1BZhEHOvS7VAW21fB9JZnYIcbP9U97n6JdQ0pYMrrczvQsAqIRqSrVVSkJirX2ztEvrGeYR6ydYp4XgZpY520xgjUGuEzxqXzl0nK5/95D9DxprGhT6Ea3v6bOo+I9+YbNcuqa+rc6RarZO5vumZRwiXcuETqbZKIXUbdwu110XwTToAAAAAAADgwHhIBwAAAAAAABwYD+kAAAAAAACAA+MhHQAAAAAAAHBgPKQDAAAAAAAADox012vgsdVcx6vzE0bSn7o8E5dsJkmnIWpqbp71VhLPbq9u2L5PhPbHVjcHj3F7ud1eTeQ7cemuhSTYlBYYkwhXLg3Ids1xartqMWN0ewqVdLSUVhZCNV0oUZv6DU5M6utqFhLPQsqXDTSN6ae+2W4vJYWZ9C+bciqpu/SwlHiWUmZNglHaFyUp9dUlI6W+E/P6UtpqSn2dbZ/bNsU19T3214Y+S++JSd0Kx5xrX81SOtrwhNhKAmD8a7nLTn2tXM/idXL7xSzDdTm1n5oLzGm8P2yfa+lekhLEZ237fjQJ19pp5QITzE3K2jykld1s26/vbguJZ2Zu4VCOwllpVfLTUv0yNwd5SlB19Us1hd7VL5VaJ/VNx5Y7FtPx6Y5xdy5ItXPKnZOSQjHgu+6rponXstRur7Xpurw96VzrhOuZO4FSsrw5LdN9PN277Lih3Y3cQspprHVcmuwipNC7zzMpNTakMsYk16EqNY1ka7mYZF9IoY+Jra5+MW3r9u0xVqmmsXVKSIKNdfv2vguXFztGJbH+vHYrfrbZUUqhN+3uM966r293nyEria3uM+y6Pd1LthNU3WdpSZqaFz4NdUPFcnLi2919PFQId8O14djMbxaOi4kZ4+zefNwkX4+Jb9IBAAAAAAAAB8ZDOgAAAAAAAODAeEgHAAAAAAAAHBgP6QAAAAAAAIADIzjiGvid1VSrtCrnfazMgpSn4dltWmR5bhZUTospu/a0EHIpOCIsavm4aXdtknRnGeZs2l2bJM2XJjgiLP6bwifs4qGlhdZD34K4QKubcly43rQVAiIkqR1tTyStSbo8dhtMiwr7SbtprJq/TE7cAslhUei4KLBbODmEQbSFmUcMiCitYuv7jsEtLh0WnLYLNccQjTCGW5w6hj64viH04diP4dpX7jiU1M2xvDJt63bbbPunc8ovfO67VhZV31uYhKRuFlR29yjJL768WKZAn+FBP3cX/tp+ZC48k3AxSmEQjlsIWQr35kIgkyQ9ahZfvmnCJCRpZoIjjsOCzK6vW0Baqu2LinRcpP3p6pfTEF/h+lYCIiRf16QwiF1rmtQeaxpzjKeFyNO54861tCC6e6/cuT4We2gUF50v1Tpm16Vrewp4aGZ/mNNsPXa6kLu5pZwm055CH9qR2d48TC7UOs3UKjH4yu2jVKfsq35JdUoKxrD1YKoFzD5KYVhhf3bznqRax9UpqdZZuuCIcCzbmjv0j7WOC45IIRPpsB+j1qmofO4y7SkgIn0mdJ8h3WdNSbo72b62H01CAFSsX3Y7p9I9eJnCJwv3W1e/5JpmOwBj3e5qnYvVcr93n9DOXfFNOgAAAAAAAODAeEgHAAAAAAAAHBgP6QAAAAAAAIAD4yEdAAAAAAAAcGA8pAMAAAAAAAAOjHTXa+C3F0/RqUt+PGMVnse6FJacwOK3YVPTYgrd8VbbSUh3vb3a7pvaby993zum/YnQ94lF2J5pvxMSAE9M4s584ffFMiQRrkx7TEIbI6qokoRWSVGyqU1+4H7kE3RccFcKHmp2X6SELt/s5jdJCV2mvYV0rckiTNq0x8SzZSXxrJDumpJgx1BIPHN9e0pYc2lz8u9fTnc1KWYp3TUlk7kktOPU11xrTdu52zOX4EpickxHqyQcVi45MfFs+PUsXfvc9XMRjq10f5ya5K7JCPHYKX3R3StPQpTvyXT7HvPIxCd73pj4dLPHXBJa6OuS0I5DvKTrG9NdQ2rario1jeTrl1Sn2MT6UKekFDpX18SaxtQkqe8Ti+Hp9I/PQ19Tv5yEcySdOy6JMNU09hwuXAPW7b7Zcpsb4doXU+htumvoG2s532y7ujqshX2f6jNzD53M/bnaTF3bQrprc3WKpO5qnVS/uJqkmmS/q0qKq+TjckP9YuuUmO6aahLzOSH1dXVKqD1cCn1McY31y7C29fZMW/qcUEh9LZ3voyS+hjm7xPpwnUzXT/cZ8iQktra2fW2vJqyvzA4p1TThWcGJSZ6VpNumrkk1zc12utWW65RC/RI+WLq65mxN88aFT5AdC9+kAwAAAAAAAA6Mh3QAAAAAAADAgfGQDgAAAAAAADgwHtIBAAAAAAAAB8ZDOgAAAAAAAODASHe9Bn5j8Sw9NvepKfcsC89jXVqLVEt3SWO4hLSYbhfbt8dwyWap/e7Sj3s3jWGS0O6GxLOT+Xb73CRjSeeku5o0oJju6tJPq4mvlbSjQmqaC/hJYV450c2lEoUUJTOPFCTqkj0labLY3l5MnjWJVykdbRXSVpvZXkvJrKa9hcSzGOZ0FdJdQzqaTTxziWnntMuMsUrpvCbFbBXTh4cntlaSWWO6a7icu/ac0mbaUjpaaLcH0hhJaIk5nNO1zyWkLUKS9mk4IVoh9cwmnoWdcRpiIE9M+x2T4ipJjy+3E89uTUPiWUxs3U43S6lprm9KR5uaNyolyLm+Y0g1zSrcS1z/VKf4FF7/PlVqnVSnuFon9U11imuPKfSmfkk1zWk4p9y5llIL7Tm8n8NizdY0IXU0XfPdtSh8YvLHXO2+2u29MvQ1978+Da8v3Y9M7TEJKZ6uTpmEutam0MvXNanWsQmxqU65Iumu7hBIdaZLgo21Tnj/XP+UTu8uUTmZdXhNk9Pp3fZs17A93zedf6UU+hHSXd1nrB6OZZvuWkisl6T5xKUrj5FC7+fhkrtPpn7n3zL1y52pTyZ/PNQeNybbCamuTfK1Tk5x9e0Tk+Sa+jpna5rHFv41jYVv0gEAAAAAAAAHxkM6AAAAAAAA4MB4SAcAAAAAAAAcGA/pAAAAAAAAgAMjOOIa+LXTZ+rR0zcvMr1MK/kbbuHI6oLMbuHkRVg4e27mlvqm4Ai3KPdpCINwfU9CX7eY8npsE4wRFlN2IRHzeS04oi/NfnZtUlit1netiDkOZnHiuLC7m8f5+SbbzCK9fZIWgHZtfm5hTVI1s58nfg1Uu7jxJLxPLW3PLqbs+7oFku0CyzovOCK02zGGd7aLXks+aCQt0usWUy4snL0eY7spLXrsF9/2w6ZACXeJSsecXUw5LUSe5lxakNm0peCIFLBi2uNRUVl8OR3jLjQnnFMrc8wtW1h4ufkX7l7LMlzP5mYx5RQKlBZZvm0WVL4ZwiCOJ9sLGR+ZNkm6EdpnbozCIstHE/9GVYIjkmm80G3btaaRfF3jFsiWQk2TQiYKISGp1nF1SgqIiMecqWtOQl8XBpHGTWEsy8X2vluFc9Wew4UwrCjdH1xb4RonhWtl6RAfHjCwnoeprQr1S1r3fBIW95+Y2iGN4eqUiXn/JcX3z9Y6aX9W+ob2UeoXO3Aaw7SF98/XtWHcGCix3RaDr2ygVtieC5kohGGl9hhgUqlpCnOOQSpjBF/ZC0wY2Fz70nu6nOz+vanuPtOn+1y6lxyZYKGp/xBz23xAOg59XU0j+ZCIVOvMTN2Q+6bgiO03MAVc3a+uuX1ycu6f74pv0gEAAAAAAAAHxkM6AAAAAAAA4MB4SAcAAAAAAAAcGA/pAAAAAAAAgAPjIR0AAAAAAABwYKS7XgO/eueZunn71rl9ViOko61CrJFLSEtjuIS0lLCWEmlcElpMrzHti5RiFtJW3fxi4plJ8lmF7a0WKQnNRVRWYsx818S9VSEk0ac9hsQlN48Y3BbTZE3XlGDlUqJS+mlKYTVhQCE4SM28JzFtNSWeme6ldNfCuNEIacCVZM+YrlVIgo3prm6MmN5XSVgL7TZVrLC9lFaWUtPc9ip9YzpaSEJzr7uSjhbewHR8dnc8p8RIc5EKp6pNPJOklTmHl9PhSWinU7/FuyHddWrSUmeTG6anNDNjp0TUlMJ6VEhCcylmKdlsYi4asW8hxbUi1TSV+iX2Ne0pmXUR5uHqhpRSa1Phi3XR3NQNKb3PJcuntPllqHVckqutXSR/DoeapsWL/m4JnPkani5G5nipbC8FkMdrsNleuLb7ZNbhNY0U6pfwAu32Yp2StrdbYmsxPHr3uqaYAmrTXdPHrkq6a6yBTFuqjd1HilJSauobrp+FOsWmyRbmltrjvnf7rZr4WvncZRPrw70rbs+MEbbn6peYYh4uRieL7TfrKNQ6M1N7uDpHyumurn+uabbHSLVHGqNSv9jfP3Pxu3v3zuDfuwi+SQcAAAAAAAAcGA/pAAAAAAAAgAPjIR0AAAAAAABwYDykAwAAAAAAAA6sFBzRWnu33vvPXWRDrbUP672/5iK/i/P9+hPP0PHskQv9blpQ26ksspyWYHSLHqdxl2Fh4ZVdkHl4X7dYuJQXTu6mvxtX8osppwWSe1jot7TIsmuvLphrFszsZmFbSZJb5DVssJkxKoEUktRccERaKN+sSWr3j84JXLBjhL5mIeS4UH4l4CEuvFxJZyhsb49KC/IWFt+OizpXFgy3Cwin1b59s13geITQldKcC2PEgIgU/rJjcEQpICKNvQjXZbNysrtWr9v9Bt29YDEJAQhmEeLTSQiImIYFi83YacFit5jyNMytFcZI23N7LodBmEXgL/viEsSQEBcEEPq6V5LrlHR8mlCSVHsU6qJK/ZL6uvolnjspuMXVS6WaxndN7aXDyx7Mvmu89tn0grA/zdix1kmvb2baQuiDfV9j6INvd/1bWI2+EjIR609zPFeCI5IrU9NUyrPCfbUSHBHHMEFbpTCsMQK1xqhpCsFXsa/bR9WDyAa+pYAr07gIw8YwpO0N2s+aklZmXyxDDXUa6glXv7jaZd1u6pRUN4QgB1fXVOqidOpV6pfkfnXN6RM++Gss1W/S/Vhr7eMrv9Bau9Va+3JJ31LcFgAAAAAAAPBQqD6ke4qkr22t/bPW2n0fH7bWXiTpxyV9ykUmBwAAAAAAADwMqg/p7mj9zcJPlfTDrbV3TR1ba58h6Yclvfvmd37nopMEAAAAAAAArrPqQ7oXa/3NuCbpPST9RGvtY892aK09t7X2PZK+WNLxpu/3SHrP3acLAAAAAAAAXD+lh3Sb0Ig/KenztV468SmSvq619k9bazdaax8m6f+R9Ge1fjh3IulvSfrg3vtvjTpzAAAAAAAA4JoopbtKUu99IelzWmvfKelVWv9z1r8m6UMkPU9vDtr4UUn/Te/9p0eaK4LfffxRHU0fvdDvhjCn0Dek1xQiHN324rgpLaeQJmvHiImDhXmk1+zCa1IqZ2ovJAfZ7ZUiNVVKxBz6++sxTAJgSFzKCarDU8zsGDEp1bdX+ja3n4tpZaV010Lfkn2mo+2Y7ppUEtZKaWyVJLU0diWlrZqa5pLQQuqWG7v0OuIYISG28l4XEox7PECHXxv60l947GsJr2/lLl4pbTWmsLq2QmJr7Gub7RiVvkmt7+CuJZXaZd2/UqeY4zDVDbEuco2FmqYyrkKdEmso01itUyq1h0sSLdQ/Sb5uuUGGJ/lKPp21hXTsyuvLtc6wcaUR6oa0vRglOnzcUZLl99W3YozrVqUWKKa7luqXHWv82LdSL8WaxiSzFusiW6eMkIobucFXKTV9+L0kXq+Xw+sUH7rtd1yqU2zCfaob3Oc83zXXOm6MYq1T2t6I9cvi8T0VOBvVf+76h3rvPyzpRZL+Z63fk3sP6G5L+gxJ78cDOgAAAAAAAOD+LvyQbuP5kt5fb/67jK51uMRv9N6XO44NAAAAAAAAPBQu/JCutfapkn5K0vto/Q26n9F6Dbpnab1O3be21t52lFkCAAAAAAAA11j5IV1r7ZmttW+W9M8kPar1A7p/rHXy6/to/bCuSfpQST/dWvvo8aYLAAAAAAAAXD+lh3SttQ/UOr31L2n9IO6Nkj6i9/7pvffT3vvrJL23pK/a/PkzJX19a+3VrbWnjzpzAAAAAAAA4Jqoprt+3+Znk/RDkj6+9/6rZzv03u9I+tTW2vdJ+nJJT5X0Mkl/StI77jZdOLffdENT3dz/hirxfftMZ6qkala2N8brsylYxRSzSoqnGbuUouWHiOk3fVLY94VEqbjnbRxwGGOM9NNSitnw42VvyWb7HOOy7SlNrZQ6Wjx5KiltOyfPpjFC113T39b9h6d8+YHDsIXUyZjEbFPawrgx4tBd/HzX0ikctzewLamMG8co9K3O4yqrpp5v/f4I/fdUY5THHmF7pXTzQpJoNTXUD2J+PaVgp/g+V3qkrzXYy2SK4S2MMUKNWKs9RqhpqF8u5fd3rWvKtUeh7851UXF7tdqqsC+GD+HTtSVb1MRk6xyLuj3sCDVGvpXsWPdV6snqGLv2rc7jPpZvWow2llP9565t898XSfqAJz+gO6v3/vWSXiLpxze/8/yLThIAAAAAAAC4zqoP6d4g6S/23j+7937fx4e991+U9H6S/tFFJgcAAAAAAAA8DKr/3PU9e++/VfmF3vtc0n/XWvve4rYAAAAAAACAh0Lpm3TVB3RP+t3vuujvAgAAAAAAANdZ9Zt0f6i1dlPSJ0n6QEkv0HrduQ/tvf+ncaaGoSaPzTTps0NP4/p6EBex3VFcINnsjKuyey5/LfOr8srxoNp1rfcHVmmR87CTlsM318ZIS3hY3ytcP9f81lVauD7sjKuwi6hp8CC69nVNJeRlT4UDNc3V0B/f77OXCz2ka639dUmfJ+mt7zVpfdjeelK/qaR/KOmlkj689/7LF50oAAAAAAAAcF1VgyPUWnuFpC+T9GytH849kfr23peS3lvSe0j6Ly84RwAAAAAAAOBaKz2ka629UNL/sPm/PynpT0l6+n1+7TVaP8x7aXl2AAAAAAAAwEOg+s9dP0XSVNKvS/rg3vtvS1Jr5/7D5l/a/Hy78uwAAAAAAACAh0D1n7v+Oa3XnvuSew/oBvjNzc/nFbcFAAAAAAAAPBSq36S79224Hy38zr01655b3BYGOnqs6Wi1/5iWURJ79jRGaW7VObjtjTF26OuTyVIc4sDfr85jT3MrzSG090rkWWHcUVTj2Hadx3VPZxojbK4yxhgXudL2fHOrzCNtz7XHvsO3Fw/xwvZqY+xpbqG9dARUt2eULhnV82HH8+fy0yUv16XXNJXTeo/3cbu9yrihvVaH1eqGUmLrKHVR4eKw476XQl2zz7rWjjFCbVXa3ghjXGUPYv1SmUPpPpdO7N3GzWP47dlDfJ/38V3nVpxHqabZU+0xyusYYR7JmHVNe2K/F7HqN+nuPdSbF37nXTY/T4vbAgAAAAAAAB4K1Yd0v7H5+ccLv/Pizc/XF7cFAAAAAAAAPBSqD+l+QOtvTv75IZ1ba8+W9De0/oLiDxa3BQAAAAAAADwUqg/pvm7z8+Naa3/hvI6ttUclfYWkp22avrG4LQAAAAAAAOChUAqO6L1/X2vtuyV9kKTXtNZeIekrz/ZprT1f0vtK+iJJ76z1t+h+vPf+XeNMGU82e6zpaHnBxQsvO8hhTwsn9/S4uRKKUBijT/zKk3bOYdy0vWb3p5+0G6OFVTHHeH32tcT9WVl4ubC9NLdCqEUrjJH2p+/ru8Z5uLbS4s21FVDj/C5RH2Eh5F64OMTNmTHi3EJ7d4E9cc6mMQT+xHnY7YW+bh4r37WFdje2fc0K66wX92dzY6c5V4IqKq87Leps+4Zx4/tX6LunhZr3FlRxVUImLjvIoRJ0UJhH/P1UT+xa6xTrBl97hL62Lko1zfB7ZVyL3m6vGFTh+scaqlDrVOaR9sXeah3f1Y1dzuEo3SDc7xc3uCdj1C9+3HA+FMau1C+xhirc8/MJWBijcA8u1Q1hZ3R3D46vL23P9R1eN9TrItNYqIvGeX2+r60xUj0Z7FqTHDKoot2+WsERkvQJkn5W0g1J/0DSG8782Y9L+hVJ36D1A7om6Y2SPn63ae5Ha+3R1tqntda+tbX2utba7dbaL7TWXpO+KdhaO26tfUJr7Rtbaz/VWnu8tfYrrbXvaq19fGtX5VYCAAAAAACAB0Xpm3SS1Hv/vdban5T0v0n6KEnHZ/740Sd1/xlJH9F7/8WLT3E/WmvHkn5R0nM2TW+Q9HOS/qjWibQf2lr7kt77Zz7pV39I0ks2//sPNr/zbpI+ePPfX2qtfVzvpb9vAQAAAAAAwEPsIt+kU+/9jb33j5H0XpK+TNLrJD0maS7pNyV9h6RPlPQevfefGWeq4+q9n0r6Ya3/ue67SXpO7/1Fkp4l6as33V7eWvuAJ/3qv5P0bZLeU9Kzeu8vkfQMSX9v8+cfI+kv73n6AAAAAAAAuEbK36Q7q/f+k5J+cqS5HMLH997vnG3ovd9urf11rRNs317rbwu+9kyXzzK/s2itfa6kPyfpT0j6aEn/fK8zBwAAAAAAwLVxoW/SXRdPfth2pv1U0o9u/u+7DPydLukH3e8AAAAAAAAA59npm3TX3As2P39lz7+zs+PHuo4WF1wCb0/pZimN1CeJDh839a8llw4fN28vTM71naZxQ/qQ6194nJ62F/uXUuEKiWeu7zSluIZkMtM/pZW59kkx8Wzixpj4qCJ3iMe+tlWamv4p3TUdcnZ7hbiktI/GsNox3WyVUlULqa9pDq49jbsK6ViuPfXtq+0TJb6OkI7Vl2aMlNy1NG3FkLZKalrlZhLTxmzyWhjDvL40ruubxq4k3U7CuJUU3VISbCXpVvLJZDGR1kXyFbc3dA5jKEZ0le5z7uI+QrprrFNc3VBNoTftq1QL2LkNT5CXfJ0RU+Fd3zi1MI9UO9hBhtcpOcm+UOtU6pRpqBFsWm6qX7bHcLVLtT3XHoXaqtCeapqrUL9Uahepljjv6pr0+0tTN0j+/FmlvqZz6uvql1RjpH3k+qcx+tIluqeaJozhmuO90sxt+MfVqJIgX6lp1v23Z5L6luqGEeZcSXcdYwyf7po+KIR5jFi/TO7sN37gLR7StdbeQdI7jTR224x/U+sk2Ju9968baey9aq19kKQXbv7v/z7wd16s9T93laRX7WFaAAAAAAAAuKae/E26T5L0uXvaVpd0ZR/StdZmkt5Z0sskvXzT/Hm99x8553emWq9b9+GS/ntJt7QOnfimHebx2JBuFx0fAAAAAAAAo/j/2n2+Ctx7f+rQwdw/d32oHgC11p4h6fef1Pyrkj6t9/6d5/zeGyU9/UzT70v6ryX9H72nf6wEAAAAAAAAbHvyQ7qfkvQ1A37vz0t63uZ//4TWa7Cdav2A75mSXiTp2Vp/e+4fSvqBMSa7Jwut01tnkp4r6R0lvYOkz2+tTXrv3x5+799q/ZDuOVr/E+Fnav0txFlr7VWbIImyIU9YW2vPlvTbFxkfAAAAAAAAo3jn3vsbxhrsLR7S9d5fI+k15/1Ca+3vSPorkn5d0sf23n/Q9Dne9PlfJH26BBqo8AAAIABJREFUpB/YjH3l9N4fl/SB9/5/a+3tJP1NSZ8p6V+11v5q7/0rzO/9V2d+51mSPlnS35f0VZL+iKS/u9+ZAwAAAAAA4Loopbu21l6o9YOoO5I+oPf+S65f7/1U0le01n5B0vdK+urW2h/rvb9+1wnvW+/91yV9VmvtRNLnSPrS1to39d6f/E9iz/7O70n6ktba72r9kO6zW2uv7r3/7GXMefZ41+xsuuuO/2C5kuK67m/SclJi1t7SVkPfqZtb6ju8PfW1aWrpHz+bua0HN21Hhf0Z38BC1E0hQS4mtrrEs5RsFsaYmv4pNc31dW2SdBRSWF066zRsbzbdjleaFtLR1mOblLbU18QdVZPNKv1d32rimTNK2mq4GC1MYllKR3PbW4btLZb+hF+a1C03B8mnqS1NWuu6Pbxuc57YdDRJfWKSYBfFxDrTFpNZXXMldVT+UpRSzCbmdcd0tEVodylmcXumMVzbU+prW5prYiUJbYyUtsJ7UklN22NIdMk4KfTmPCumrdr2PdU0kq89zCUg9s0p9KHdJVSmudkDZvi46z8wbZX3upLiKvn6JdRhrq5JNU1KgJ+a/ql+cWOkmibVL0c71i+udjlvHhPzBo6RJluRxqjUNWP0XZqTKvYdoX6Zm/rF1S6Sr1NSTZPrl+Ep9EtXp8QU13AxsvfbwueramK92c/p8LT34FCHxfrFJtn7vq72qIxb3d7ONU1oHyVN9hLql3Z3v4VP+iie/DWtw9RfmR7QndV7/35JXyvpGVqHUjxIXr35+RRJ7zXwd/6lpLnW+/Wl+5gUAAAAAAAArp/qQ7o/rfUzyO8v/M63av04+sOL2zq0Pzjzv8PfMW65o/XafJXfAQAAAAAAwEOu+pDu7Tc/K/9s9dc2P9+puK1D+ytn/vf/PfB3PlLSo5v//SPjTgcAAAAAAADXVfUh3b1/+fvCwu/8sc3PpxW3tVettZe21l7dWnvfTdDFvfbj1trflvQ/bZq+6N56dK21d2qtvaa19mdaa7fO/M5Ra+0TJN0LmPiXkv79Jb0UAAAAAAAAPOBKwRGSflnrB3R/RtKrBv7O+21+/lZxW/v2EZI+YfPforX285v2F+jN++XVesuU1g+T9KGb/3pr7Rcl3Zb07pJubPr8n5L+au9xJcrRHT+x0myRVllcS4v/VkIm+iQt5mkWWY7b2x7Dhi2cM4abRw5yMHMrBERI0upo+PbcGdXDWRbfMbfoceps9lFaALNyQLr3NA4eFiZ2iylPwsLLE7OIsSQdHQ0Pg3BBDkeFvpI0M4se575me2kx5bBCq+t/FFZB9Ysph8Wi04LMhaOgslBzaTHlcNFxix4vwsWhEhyRxliYC4xbYDn1Tf3nYZHl08V238UkBFJM/RiLxfZ7sgwrK6zMFabHi2phkeXds0NycIRdQHj4IsspICIGOZj+MTjC9g3nWWVB5hhUUViQOW1vZcYo7Psc8OFSA3zfJC4M7aZR+evkQsBV2hk2yCEGR4S6yIUzxCCH4TVGD4EEzYZBhEm7OiW8H+X6xWhu35tjUzon+KPCjZHGTcdWIfhq4kIfjobXNJIPYnDhDpJ0bMZ2tYuU6xdXZ1T6HoWLXCU4Io3h6pdK7SJdjfplFU54t71UY6TQhoW5QKS+c1MDpVrH1S+xLgrBETaoIvR11+CVu5jJ1zRSqGvS9czdBopFjT200uFWCEWohFbFWqdU04R7SalOGd7X1SPr9mHjSrIhHzk4ojYPvz33+2f+90nlTlhXfUj33ZLeQ9LLWmtf33v/zvM6t9beX9Ina/0yX3uxKe7NyyX9kKSPkvSukv6I1hkxPyHpZyT94977jz3pd/5XST+r9YO9F0h6F0lvI+mnJf2cpK/qvX/vpcweAAAAAAAA10b1Id0rJf0tSTNJX9ta+yxJX9n79t/BtdY+XNI/0TpAoUv6yh3nOqre+0LSN2z+G/o7XesHld+9r3kBAAAAAADg4VN6SNd7/5XW2udK+kKt15h7paTPaK29Vutvn93R+htmL5X0Eq2/YN4l/Yve+78dc+IAAAAAAADAdVH9Jp1671/cWntrSX9n0/Sum/+e7N4/8v5+SZ9ysekBAAAAAAAA11813VWS1Hv/TK1DFH5a64dx7r8/kPQ5kv5i7/3uKLMFAAAAAAAArqHyN+nu6b1/e2vtOyS9SOu01+dLekTS70j6UUnf23t/fJRZ4lzTO0sdLVMUylopMSullcV0VzdGmIcZw6WnrvuGMUzwT0qIbYVk1lU6G0yaTC+lQXku8EzSOr7kyX3DvncJaaUkPGnnJLQW0+ZMQlchxXXdblLMColnN0LfG1MfjeTGPnZxSZKOC31nIX7IpZtV+k7D0ZWSzVx7GmMMS3PApCQ11z4PJ6tLNpN8illKTTtZbp/wi6OQzGr6StLdyXa7m4PkE+tOQzLy3CTBVnWT+mpTOSX1lHblDo144RoujrBjElolxVUKSWixr7meFcaVfBpsToLdre+63bQVEtbcvU8KqWlp3D3l3MdbcKVOCceyTXeNKa7hnDLJqqlOcSHPMfU+JLY2c4lahTR1tzPCWx3vDu5KGYKY1dz+TOd1et1hHn4Q1xZGiOn0Lt3V93VJrqmmiWnxZgxXY0jSjaPtC4xLm5ekm6HWOTbtKVne1Uupb5qHS7ifhRhIn2Qf9v2e6hdXu0i1+iWN4RLn5zHdNdQ6pv9p+BBzatJWT6e+r0tmPTF1zrrdz62Z92qeTmwj3D7VY+S1+RyULka2phk0rTODmyHS9czMI6efhnazQyq1RwsprqUxCgmx5TrFpbum5FlXZ4yR7lpIfD077ur0/Gcvu7rwQzrpD4MUfnLzHwAAAAAAAIALuNA/dwUAAAAAAAAwHh7SAQAAAAAAAAd24X/u2lp7nqQPkfRfSHq2pKfe51d67/3PXnR7AAAAAAAAwHVVfkjXWptK+kJJf7Pw+03FdV8x3PTuUtPlmZUeCwt7lxZTDgsWu0U30yLLLiSiLWrb62YxZBcQIfkFKXNARCUMIixqafZdynxI32N1652mBUX3dlKlObuVOMPqnM0syDwJizRPp8MXWXYBEZJ00yymnAIiKosp35zObd8bZnXUG2El1hsTP8bRZPt1p8WUXXvqOw0HzCQdSHuyMosbL8MCyStz0KXgiBTO4PqfhBP+1nS7753lzPZNi2S7/ZkWb04LX1e4xaljEI65kMTF78N56e4FPZ3v9kYQLiRpc4UFme0YhZAJKS3IPHyR5XBa5zFc2EVaINkFR8S+fh5tsb1D8gLX5r6a0gTCAtB23MsOjkhMPdFjcIQ7d8I8QtiMC5RIdYo7h1PfVVrgurRDtsdw118pl5N2f8SFuk1behnpeLGvL4R22AtJGDcFR9j6xZ88rn6Zhr4uIELyIRG3ZqH2MHVKCq1KtY4b41aodVzoQ6p1Yv1iQ7JCcIR5X69yTSP5uibWL6Y9hmStQvDVartWOVn5/ekCHu4u/X47NQk0LgjivPYKF8LnahdJWoVz1Z/vw2uafbLXvnDIpsAF157v+YWAq0L9Eseo1CmFQK0UcOXGTn3j/cjVL5WQrDOm8/1efy7yTbqvlfQxquehAAAAAAAAADBKD+laax8g6WO1fj65lPQtkr5B0q9I+h3xbTkAAAAAAACgrPpNus/a/FxI+qDe+/ePPB8AAAAAAADgoVNNd30Prb8t96U8oAMAAAAAAADGUX1I97TNz38z9kQAAAAAAACAh1X1n7v+kqQXSnrjHuaCC5qcLjXpIerlnpRi5ppTAmBIYXX9U4pgM6k9q5TiGlKUYvrXUCndJ8S4uNfSwuNtl7iTUuHiCo6ufYy+Y3DJsyFxyaVZpnTXlHw5Nf1nITXt2EQgucS0dXtKQtuOO0qJZy7dLPVNKWY3TbzSGOmuKfFsauIv95U6msy7v+0sK+muJoEs9U/7yKW+ulS589or3D5aTMM1LoyxdNfPsO/duRZC03KCYyEeyvUtp0u5QSqJkTG5NLWbpLCUsObSXUPKaUxpc6lp85SEZvqG7bm+Unh9JvE19c3paCY1tnqKVG7kuybWS+rmXMvJpaamSSmuK78/u0tnTSnttdWhQ7NLOBxe96X7uKvZ1u3D2qTwVseLy/DjohRoW7z2ufZYv5j2o1B7uMT61O4SUSWf5JpSXFNN4tpTCr2rderprtv9Yzq9OQZmKRoy2LWuSffVZfhui0t9HSPd9W4LifOmDp6ExNZJqJeGSsnPLplVkpaT7X3h2qRxPie4dPpy+nfFvj6jpXZ3rU3proW+lfol9q3UNIUU+kqtE9NdQ7tNrY99TeOZ35/Ma9elquo36V67+fm+Y08EAAAAAAAAeFhVH9J9mdapri9vrT3tfp0BAAAAAAAA3F/pIV3v/eclvVzScyV9X2vt+XuZFQAAAAAAAPAQqa2E0dpTJX2jpHeQ9Lcl/UJr7V9J+o+STu73+733L7jIJAEAAAAAAIDrrLqa5NnAiC7pWNJHFn6fh3QAAAAAAADAk1Qf0rm8lKEZKvvMnHyotcVKLcVqbfQQYWVTSlOaV0oKs+ms4V9Su4Q831OrEIXWzaRjaowZvJ4AWOhrUmNajJvzzSVjjDFCQprva/ZFTHENaaSmPfV1Yx+FA+MovIEz0+7aJJ9uVklxXY+xn3RXl+IqSVPzWiYpcrBgleKOjVlIonbpZi7lTfKpv5I07bu9lpTcdhQiml177rs9N5fQJknzpU96c6HZ6Zyyl+uUeDdCwm9JZXOFKVdfXuXa7uZRSY1N7aW+KfEsjWFS0+IYlcRWl2hauCeOoVTTyL+WmABoU5fDCwyprzKXuRhSO0pS8Z6SkUc4/8Y43/em8AJzEqxJOEx9wzTcvSDdH3zfUOsU2lNiq2tPNU2tfhmeEOtql/PsWtekmmYZ7u9L038S6hHXnvpG7hIcDjqXPLsI+9O1pzr6tPk6pVK3u/ZF5SKe7OtzTbLHz3P+WptqgWFt1fbY1x0aqW9IlndJrqXE+pBqnGodGzee6hS3L870Ta9pLNWHdK+UNJd0uvk5Fw/fAAAAAAAAgJ2UHtL13v/GviYCAAAAAAAAPKxK6a4AAAAAAAAAxsdDOgAAAAAAAODAqmvS4UrqedHDjbQuZneLDYdF2eOquW7R6jSGa46rN4dpFNZ89NtLncdYPXS4y16rfW8Kuy0unD2CiV282e/kSvskrJjq+qbFjScpAMG0x7AEs4JpJSBiPQ8zxhgHYmHR4xTOYPd9eH2TdNCZ8ImKvD/De1I45px0XKQxUkjEdbG3lzfCuKUAhNL9KG2v0jccL5V1tiv3yvsspjxojB2VahrJBwGEUdyMcyiC/4Nu6po0hh1ihBCGfMwWbsR7Oiev+aUsXqtj7VHY0b72GH6PkvK9ble5XjJhCaV7ZW2+O9c1qaaJF57tplXovDI10CrVAqkeLOzPyr4fo47eVflzwuV+dNubK3NNtMEKqe92U/l1VGoP1x6Cs/L27Mlq2XvoW7Tt902zD+laa39a0hdoPe1X9N5fu2n/3F021nv/H3f5fQAAAAAAAOA6St+k+1pJb6v18+lXS3qHTfvnabfHhjykAwAAAAAAAJ4kPaR7XG/+Aunpk/7sol8svSpf5AQAAAAAAACulPSQ7iMl/V1Jb5D0yjPtn7T3GQEAAAAAAAAPGfuQrvf+OkkfZ9q/Zu8zAgAAAAAAAB4ypLteC+0t4nB6JRrHBS2G308hrH1qBklzmJr2NN1JYR5pDNu39vr25bK3F+36D9EroYcxsGf3nbEyY7i2avuq+zRS13eZ+oZjbm7SSGM6mjtZq++dmUYKRqqkqbm0Msnvj5TuWtv3w7eXEtbcPNLcluF4ce3p2HJi+lsYY4zz5CpzLy+G6VV2Ren+kLbn/iDFrY6xPdO1NDf5RNN4HzdppEvfVyZEuae+YYOVtNxSTeNqjDAP95rXfd0cfNe874ePUTkuSsGshbmNs7399H0QpWt1rD0KO9rXHsPvUVLt/ufmFmud1G6ulasYqehqHf9xNSXZVxLnU/3ipNftarlYT5j96X5/Pbf0vppaJx5zlb6719EV7jwpB4Jf4meYfboyn0HdfSNFthbucy18prdjp3uzC1t1zyAkaemvAc29vklIaXdp8W/xvGX4teMi9js6AAAAAAAAgPviIR0AAAAAAABwYDykAwAAAAAAAA6Mh3QAAAAAAADAgREccQ302VR95hcd/cM+lQV90yKMaR1HtxhkWPRxZcYO67BqFcZwiz2HNVft2Gl79nWMMEbqW1mQeYxFnUdRWWDVLghbW4B2uTKLG5s2SVpOttsX4cBYhEWF5+bNmqz8ZTIFPFQszaKjaWHi2WSx3RZWdp+HRZbdnKeFgIgqt3ByWljaLZycFlOutN9dzWzfE/O+nob3OrUvVtvbW4Tj83Q5vG9c7NssYrsKY6zcgrd59fvQ7pt3tqdrX/nluTHitX27zbz966m5fZ/GCFWYW8c4nan5b1vNnyzTAtDb7T0sFm1DH8LCy2m96Z3X+i4HTrnFotP7NDxkIrWvjoaPYfsW5ib5Y7FyLI9x7oxSp+yrpkkHXOEF2muq/DV4ufJna7q2u3uBu2dI0pG57y+a73saTsBpOjGHSnV72EeuVpmFObu+qd6aFgIiKnLA1e5hEHNz0a/WOq5+cW2p3dW6krQw7amOTrVqpW53tX8MyCp9kB3edRR7vPb5z6Dp/mDu4+XPvIUx3KERw3ECc82w4Q4KgXYxtMpP2pUvLYRMuJCst7h0xg/44+CbdAAAAAAAAMCB8ZAOAAAAAAAAODAe0gEAAAAAAAAHttOadK21t5L0PpJeoPUDv6/uvb9pjIkBAAAAAAAAD4sLPaRrrb29pJdL+mRJN8/80XdJetOZfk3SiyW9VNJX9t6fuPhUAQAAAAAAgOup/JCutfbekv61pGfpLbNIUpbKt0t6W0m/LOk11e3h/pbHUy2Pz093jWzime8ak2DMGDEJ5mh4umslCS0l5Lkxct/h7eOkyRbaK8lB+0x8NWd5T4lnlXS0lHi23N4Zc5PiKkkTk4Q2GSHiKaV5uTkvprWELp945qOKZiaxNfVNr3s6QiJthUv/WoUD1PWtJp65hLyTkO7q0s1SiuudpR/Dpe/dDX3nJn5xHtL7FqHdJaSlc8olWLm289ordg0LjIOElK9SMmu6Xpt911f+hdj7QEyrTjvDpZgN7qq2GJ7GJkkTk+TaQrprc687BsEWomfHOLjsJMKOS8eA619IrI/1SCWFPtQeboxc//gxdk69j6+jMEapTtn9uBgjObhyTUx93TU4pVmma/vcvJhJeIHpHmPnVigIc/KsSSYPddjRxF8IKrWOq1Ouck0j+f2c0113T7J3tcrJcni662mIJnfH1snCj5vql7k59l0tL43zOcF+BtlnuuuOydZjfCas1DopADlez2z78DosPCqI7SuTwppO674wMwv3KFunyNc6PVyv71e/LNPGR1Jak6619lRJ3yTprbQOvX2VpI9K/XvvXdI3a/3uftDFpwkAAAAAAABcX9Vv0n2ypHfQ+u9MP6z3/q8lqcW/CpYk/djm57uUZwcAAAAAAAA8BKrprh+p9ZdIv/zeA7oBfmvz83nFbQEAAAAAAAAPhepDunvfhqusLff6zU8e0gEAAAAAAABG9SHdW21+vv7cXm/pbTY/bxW3BQAAAAAAADwUqmvSvUHScyW9o6R/P/B3/rPNz98tbgsDLW9O7p/uGpNZTWMhTW/d3/RN27NJqSmxzm/OJ6EN7xvT0WLymumb5mzGqKSjSbKPztO+j+/JrlKgjdlgSnd17auQ8BRCC9VMullb+B3qkhZTGlRKPHPpoMcTn0B2apJcXYqWJM3CGEcmwij33W5P6W/TELXo+qcxKtJ+rvR1KWYpWdcls0o+hS71delo7v1PfSWfepZS01wSbE539fNwiYHLlJpm+qZzVYX2VrnopGNrhMRWN0S81qYER3Oa5GPZJILF6/LwFNa2DH1Ne5umtLIwxsJsL6TXumt+DE40O7R8GUnzcFIsnJHrlOEp9JWaxqXNp7FzQqzr66eW0mTdJSqOYeZcSbeP7TEt17RVkmCl0gHmrlEp4TBd+2z9Eu4PS3PLnoTE5UVKVzYpgekVu2tUTJN1cYjy97qjiU+NPZ5sjzFrvu9RqF9mLrE1XGBc6mu1Ttm1rom1Y2hfmsI99V2Yk8clvq77hhrI9M+1x/A6xY2b0l1dTSNJc1OTpNpqYer59Dmh8lkj35x9864qn9FMmOm6byWxNaZSu89BqQ4LY7h7ZUps3TF5VvL1S0yhP9oepFLTSOFWkpJgXfOZ7S1DyvVYqqP/h83Plwzp3FqbSPrLWu+q1xW3BQAAAAAAADwUqg/pXqP1s9f/trX2nAH9P1XSe23+93cXtwUAAAAAAAA8FKoP6V4l6TckPUPSt7XW3t51aq3NWmtfKOmfbJreJOmrLzxLAAAAAAAA4BorrUnXe7/dWvskSd8p6X0k/YfW2ref6fI5rbXnSnqRpGdq/a27LunTe+9vGmnOAAAAAAAAwLVSDY5Q7/3/aq19pKSv0fobdR+rNy/P97IzXe+tBPjZvfd/vtMsca7lo1MtbtwnOCKwC1xXQwoqwRGVRSZju1kYs7C4cXVB5koYRGXx5srrdguOS/LvVTVMwi6iGfq6dX7TYq5m4fNVWjE1LBLqLlFuYVRJWq22J5cWU56H9tnEBDlMQ5DDYnvh5Ng3LJDsFllOCx67MVLftCDzdISQiIplIWTAhT6kgI+4ILN5X90izalvXLw5HUdm4eTc12wvLLycwiBKiyyb88+1rQcpXDRGOITSEO7yUFqEOK0fnK7Xbg5pV5j7TlrcOP0bhZXZz2Gdda1MewpyyO1mzmGRZTdGvFy4kIm0OPUlS6Edtk6pBB0UapryGOb4TDVNPLZs7VEIqkiBWoX6pVTrpP3mm2sqgxTql3RsrcyLmYfNpdPE3SuXYW6urplPQmhAOJBcSEQMrTJ1kQuySn2lEFoV3qgUPuFclZrGhVzl4IhC30KgRAyiKvR1tfGiGHBl67AQ+LZ0nxNCDVWqX2Lo0cC28xQ+d7nPbuk+EK+J7v7gu1ot3Xf+f/buPUi29S7v+/N2T8/e+1wkQChIIHRBBCe2kRAgYy5BMuRiLg4GoQJiu4Tt2AYjTCjHZRK7bHAZOXaoGAyBwkUCMQVRjABzLUQkylGMMSBFgDEBRYCMuAgkcdHR3mdm+vLmj+4tjU4/z+z1m149M3v291O1a5/znrfftXr1Wu/6zTo97+MzZWzoTayhzGcSa5pC/ZLCsHzfMYIjUt+zT47l8fmevQx1rliK3vsPSPqjkr5e619lbU/4I0n/UtIn997/0e67CQAAAAAAAFxf5W/S3dV7/01JX9Za+3Ktf731mZJuSXqHpDf23t8xzi4CAAAAAAAA19u5H9Ld1XtfSfp/Nn8AAAAAAAAAFJ3r110BAAAAAAAAjIeHdAAAAAAAAMAle59fd22tfZmkLxtp7LYZ/6akG5Ju9N5vjDQ2Tjl5aKJ+85zPWyuhfpXUtJRSU+gbk25KCYAmoauSFpi2FxNiC/sW09RMmkw6RoXEusi8IAVm2aCbkOTTTVTjKmQS9XBAu0ntcQmJkrScbo89D+mL05CENnUpZmbc1DcljaUUVjdGTmwdno7WivtxkVKKmUvtraa7utS7lLZqT+VCMqvkU8hS+ptLbE0primx1W3PpWtJIQktpqP5Zts+xilUSNVUSLZOyZVOCq91abIpgLq7tNVwPF3S5rq/6VtJWy0knknyKaxpjF0Tz8o3nj2J6fS7JdnHOqWUEFvpOzwJNo4xQup9pV7KY7iEQ9+3nE5vN2ja0jWStudSJ0Mh1s31kFLoe6F+WU7D/cj0jbVHSFuduTFCX7fHKcU1jbFr/XIVahdpf/VLpaZJ/WPtYRNb/biub9yHdC6bsXNNMzyFvocxbF0TjkWzk3tt0nHdcyq8eb37GW89im0tpdC7fQuHraX52hzPmCBvfsbKKa7D68+4vUJNc97E1vcd4+xzYx7m6LE8sZR8P0nP1r3L8NN7PfRsuxozLAAAAAAAAHDFpOCIIY+V+6bfWX27pB+VdFTcLwAAAAAAAOCB8T7f0+u9f1XvfXLWH0kvkPSbWj+c+xZJz5P0ZEm3JD0k6UMkfY6kH9/0eZqkL+m9f/bFvS0AAAAAAADg/lH6ZdrW2pMkfZ+kp0v64t77X+69/3zv/bHe+3Hv/aj3/lu9938h6ZMlfZOk50t6ZWvpt6IBAAAAAACAB1v6ddfkL0h6pqTX9t6/+ayOvffeWvtSrR/WfaKkz5b03efaS5xp8XCTbp1vtd3SepkjLMhcCTqoLCxcWZC5Ekgh1RaRtgsnx+2lBdF327eY+lBRWmQ5fICL7aaeFvEPi6e7Y5QWyneLLKfFhls49nbB4tDXtVcDG1x7GsMvpmy7RmnsfUmLZzvuFIiLb4d2t5hyesd2oeewAK0LbFiP4fqmBcPNYsqxb1pg17y/1NddJzEgorLIchjDtRfuA1KYz1JIj712wvtICye7MIh03zHVUlzcuNBeWfS4EhAR++9pMeUrsq57rml2DcnaY52yc00T2kfZXimoIt3HTVv1c3LthbnIzmWSerqmbOd0P3KLp/udW4WbdnPBQqlOcYFTsaaxzTZQIo9h5lo/bK65CrVOpa6532qaNEalpkntcQxzj4m1hxs31UXpejC1RxrDtseAq0J7JQyrypxzMVix9PWkdHN2xzOM4MKwwhyeLh1bNxTCIMaoU+JlvWvAlVR68HGv6WVRCDA7j+q32z5X68PxrUM6995Xkv6J1vP5y4rbAgAAAAAAAB4I1Yd0z938/abCa35u8/cfLW4LAAAAAAAAeCBUH9K9/+bvpxde85RzvAYAAAAAAAB4YFQf0r1t8/cLC6/52M3f7ypuCwAAAAAAAHggVB/S/bhm6R3ZAAAgAElEQVTW68t9UWvtnt+Ma629n6Qv1Hodu58v7x0AAAAAAADwAKimu367pC/Q+ldYX9lae0nv/R2uY2vtoU3/Z2v9kO67dthPnOHkkabVQ+dMGBkh3dWpJHeVEmarY9j0U9+1loQWUszc2MWUtp0TcKuJnzZVM/R1KUoxwco0xiTKSipc2KDpuyymu7pjF9PDbF/fNUUE2dOzklY2QrJZ3OcRVAIj/QC1dNfS5mw6WuhaSODM+2waK+Om/mkM+3rfXEkYrU/YZoiUImiuiPiZFublmDZmE1tT0q1pTMczfn6FvmNMA6UxCkm+FfsKXxxj3qqkDF/wGPH1laT3UmL98Pvcun+hb6kuGuGEcRsMUZsx9bWQIujSKOP7KNSOq5i2auIaC3WK5Gugav1S2Jyva65w/bJz7SKVJoFyTVOoX9wgMaXW9Y31QUpb3VPfEX4u8cnPoW9SSe4e+Pqzxmhu7kq1R6GmGSMVvnQ891TT5M+vWF+fo++VSnftvf+IpNdq/c4/SdIvt9b+VmvtE1prH9Bae6i19rzW2sslvVnSp2v9Fn9V0v868r4DAAAAAAAA10L1m3SS9PmSfkLSh0t6RNLfO6Nvk/SYpJf03k/OsS0AAAAAAADg2quuSafe+zsl/TFJ36v1Q7iz/vySpE/qvf/sWDsMAAAAAAAAXDfn+Sadeu+/L+klrbUXSvqzkj5F0jMkPSTpHZJ+WtKrJP3vvfflSPsKAAAAAAAAXEvnekh3V+/9p7V+IAcAAAAAAADgnHZ6SIerYf6I1B8eccAREpdGCAAcJXl2jJQ2m+RTGLYyrhT2uZjcVeLSckIakE2PSsmlboy4vymCzOzDCMeicl7k4KBKVNHwrqN8pnHsfUUtFowyOaSx99V3hH2upAUGpQTOyvYK769yCqVh09ZsImKK76vsRyUpLMZVD99eTE3b2/k5vOvepoB9XtcVe3qDY6TQj9G3lGRf2F41IG/nfS7exyvH3ybWV1KwFRIOK8cizFuxfhk6rka6Hd1v9ctVqF2k/c1z1bd3kfXLCPuWUpRHuc8VUujtaTTGqVVIbU7JuukUt6mvpcu3WNPsWk/s81weug9V57xGFuVF42r2PDwAAAAAAACAeyl9k6619smV7pJmkm5KuiHpRu/9OyvbAwAAAAAAAB4E1V93/Zc6/xcTuyQe0gEAAAAAAABPcJ416a7IwiMAAAAAAADA9VB9SPe/DehzKOkPS3re5t+/TtIPS3pXcVsYaPHoSv0Rt8otcE6lhc8v9rl94/8TAPeHURY+Dwsyl3cG+8UngivCLlxfHWQ/dQb1C3Cf2FtwCzXN1XO+o7+c7PfZS+khXe/9zw/t21r7OEn/h6SXS3pN7/2nivsGAAAAAAAAPBD2lu7ae/9JSZ8i6XFJ39Fae+q+tgUAAAAAAADcz/b2kE6Seu+/ovWvuz5Z0hfvc1sAAAAAAADA/WqvD+k2fmzz92dcwLYAAAAAAACA+85FPKRbbv7+sAvYFgAAAAAAAHDfqaa7nsfHXeC2HkirRxZqjy7GG3D36JpxYmoqY+xzn21SWGF7hXFje9heKbGs9P4KY1T2rXpeuH0r9K2MK410PHfsW0+hKxjjOtnVHt9g6e2N0Hfn7RU/jlK4smmP+1v5TCrHorBveYx0sYaxK33tMdrj9ipG+Ex27lux14mrYF9z3Aj3rtz3gmuuSk1TrV+G9q3WNIV6YpzaY3v/Rtm30j4U29329vX5Vbe3q6tQu0h7e4Plt7enOb90Hy9sLx61fdVhhft4rouK7cYoNc2On0m5ptn1UqteIztvb8fXS+e+rlca8dmLsddv0rXWPkLS39T6+nzzPrcFAAAAAAAA3K9G/XZba61JeqqkD5X0Ekl/SdJTtH5I98/G3BYAAAAAAABwXZQe0rXWlvfuZb1J0jed87UAAAAAAADAtVb9ddd2jj8/Kek/773v9xd3AQAAAAAAgPtU9dddf033XkL0tqR3SnqjpB/qvf+f59kxDHf4yFzTR0/O9dq+p8WNe+hrN5f6Fha8TdurLFjcV8PHqPRVpW9ob6vQ1bS3tL3ArpcZxrD7UVi8Ob2POIbd3gj7VtiPMRZeroyxt/CKfY6R2I+qcH5WAxt2XdA39U3/O6uwoK/dXhg3TmcT82EVxmixb9igW2TZ7UPoW/78Ku/P9a0GY9jPZPe+6XC6/q2wb3HcMEZpTW47Ye8+OcR93lGpdomDpNqjUOsUxh6lpgn3Zjv2vmoayR+7eL81fVe+c9oP21qqoWrvz9dFI/RN72/X7RXrlFgDDd1ess8wj8L2dhangN3rF6da0+waWlXZXqn+Sf3DSeT7FsZN/cONx43Rwn28ViOGE9Fsr1RDpfYR6pS0PXs8SnVKcXuFcA3/UVeDKob3v1f9suzz2raLSg/peu/P3tN+AAAAAAAAAA+svaa7AgAAAAAAALg3HtIBAAAAAAAAl6z0kK619mOttde21p5ReM3zW2tf21p7eX33AAAAAAAAgOuvGhzxYq2X6Xyo8JqnSPprkt4s6RuK2wMAAAAAAACuvepDuvN4bPP30y9gWw+kRx8+0uyR8/3mciXFLIRx2THSuCvXN21v5d+TSzFLfVcmSSulh7l9S/3TGH3pUsxqKW0yY/SYdDM8QW6MhNFSAtnStKVkM9N33X9YW7VvJTWtMkY8bpVU23jsXSRfGDcYIayxZPcEsuGpo3mMQt+YJLqfMcrbm5jrPY0x3f6w47jTMEYhNa10MlaS0Mz7WO+HaU99w/tu0+0LMyW9uWTclCrWJv6Cn5ixXVsau9JXkqamPe5zYXIwp2E9YW1Pcu2x2xhp3GUhIdbVI9W+cQxTA6V9dqnwuaYJF487oK7+UUiajpN4+KDC2IM3WElKVapfUt/t7ZVrj33VOoWaZK+p93aM9ENFGMO432qa9Rj2gqiN4e6VMVV82OvjGMU6xdc6w9NWe6hHYt3g+qdaxyZph75piMoxcjVJ+pxi/eJS78N93IyRaxrfPnF1Q+praxo/kaTS0fWPp/IINU04FUtj3DVfHg0f7BwuYk26z9r8fXwB2wIAAAAAAADuO2d+k6619qmSPsT8p89urf3WGS+dSnp/SS+S9JlaP7v+ufPuJAAAAAAAAHCd3evXXT9F0lc8oa1JekVhG03rh3T/uPAaAAAAAAAA4IFxr193fbXWD9nu/rmrFf68XdKX9N5/YNQ9BwAAAAAAAK6Je32T7ickfdWpf/+7Wn8r7hslveOM160kvU3SWyW9rvd+Z5edBAAAAAAAAK6zMx/S9d7nOvWQrrX2dzf/+PW99zftc8cw3FNu3dbhQ+eLOEqJpk5MTTM5LIuYzDq8b0pNW5r+y5BMtlhux/6k97wMKWYuTW2V+k6321cp8WxRyG1JY5hknZjQFfNyjEqiaUhmnSxM4llKca2ku6btufbwPmzf0H+yDMlBbuxqgpzrX0nWrSTBntH/QpUSz/wOp2mrkqBaSmZN6aeltNXCuKF9VUoxM28wpaalE8NVCIX0vhRgnRPrhqeYuSQ0l2y2bg9pY6Z/Siabmr7TMG4a48C0T8P7m9q+Yd9iipnZ55Q869Jkw4dtk2evSLprTGl3CarhRHRjuLrjrO25+iWNUappYpK9GSPUKUtTT7jXS+GalE99jeWkud/mFPrhCbFparAp7WHXck2yW/2S6rBUe9gxRtheC7HGlTFsMmsheVYaoX65GtPLKImtvXCC7iuFdTVNUZuur++atuf65xrKtIVzqId9XpmTo3S6VCM0zeeX5kn7uaYU11S/HGwfkEno62qPSaxTQu1h+qf7+8F0eyKp1DSSr19crST52iP1rdQvyb3qmpOT/X4H7SLSXQEAAAAAAACc4V6/7vo+eo/fSQAAAAAAAABwTjx0AwAAAAAAAC5ZfEjXWvtDrbX/8CJ3BgAAAAAAAHgQ2V93ba19raQvvfvPvfe/vvnnX9lhW733/twdXo/gqTdv6+attBL+2dzixqvwW81pgeRKcMTCrCi6CNvLgRLb7XMTECFJczNG6rsIiyy7/suwMOZi4RaEDcfTtkrdPTuPi7UPX6w2MmPERZZtsILvXVkIebII2zP9K31bCH0YZfHmXUMmUntYUduGaBQW8V/3v9jVl+1iyqmvu6biwstpjEKQg7vMwsLEaYyVWby3EgaxOgjbSws1mzt23J4LcgjDxuPsXlCdXwrbs+1hrnWhD26BZemMgAfTfpDGMIsTu0WTJenwwLe7BZJnaQwzSaWFl9PCyQdm8sp9hy8WPTF9U3jFRUuBU66uSTWNq0liTRMuQNc/BUecmLoo1Snp/Z0sTG01CWOYgCtbu0hahfrF3/ZTYJgd2PYdZX6x2wubS+071x6+b63WKYQ+VAK1lEK5CrVHNaiiEhwRb1Ru3CtS07gwgUoYRBi3Ur+kMAjXNwUPuOmsxZom1P621kljmL6+azQxx26V7keu7ouJN4UwiDSEq18KARGSNHHBEek+bvqm+ifVHgemf+o7M/uR+k7Dhe1qElePrPsOr2lScIQNyQrbu1ddc3R0dOZ/31Vak+7Pb/5ukv5rSX998+/P1vkydto5XwcAAAAAAABce+nXXX9E730u/KNP+G/tHH8AAAAAAAAABOmbdF8g6Vslvb33/oa7jaS7AgAAAAAAAOOzD+l67yutv00HAAAAAAAAYM/4ZhwAAAAAAABwydKvu1qttb+z+cdv6L3/7sDXPF3Sx0n63d7764r7hwE+8PDdeujG/Mw+MZnVtC/Ds9s8hklQjYmt21E+xyH2J6W+npjUs5OpH8MlpB1PfN/jkIQ2Nek8JynuyFiENC+FZKSVSceKv2juEtKKq0DunLpVSUcrpKpK0sSc1pXEs5SklhPWTOpP3OdC35R45lLaQl93nHO6ayFhLW2vwqSqJpVk1pg6GrbXXdpq7OvG9dsz05YkqZl01pS26qa59Fmn1Fd3XVYSW9N1lj4/l5w3Rphej8lrpj18Js0kkE1CEqxLcZWk2Wz7gMS0VZPYmlLFbhz4CWZmJojDmO66PYZLfF3vh2+/YdrdPkg+CS2NOzVnXUqCTSr9U+1R6evqGlePSD6xfh76Hod2N7ZLcV23b08Ors45az9cCt0inPcuCbaFmiZVli4RcRW2190+h2vVJkpLapXCxgWFx5qmkE5fqXXCgcs1iUmrLiTE5r7DE+7HqFPSftj+8TMxx6J639m1rkl1Q0wmN7VATPw04xbqlNTfpc2vx3B9w7imbwtPClYHhWs4/fxo2lI5mWodV7ela9XVL+ncimfQjumuLdYpod2ln4Yk2JmpU1LafEqnv2Hab0xDTWP6utpFyrWOS3K9EcZwNcksfNgpsdXVHq6mSX1Pt905PLavG0vpIZ2kr9T6vH2lpEEP6SS9QNL3SHqjpI8pbg8AAAAAAAC49i7i111/a/P3cy5gWwAAAAAAAMB95yIe0r1g8/fhBWwLAAAAAAAAuO+c+euurbX/QNJN858+uLV2dMZLp5LeX9KL9N5fkf2Fc+4jAAAAAAAAcK3da026r5L0l5/Q1iS9trCNpvVDulcVXoOC9zu4rUfCQtX3sjQrbrpFkyVpHlYwdWOkvi4k4lZckHl4wMPRMiyQbFY8bWGV0NR+YtrTmrQ2gyEsmNrTIG5RUbMo7V7FQAKzqGxaTNm9jcJCyKl/7GvaK4s0p/ZScEQh1EIKCySbcWPfQsiEpHFW/a9w5234/nZ3Cy+nxZTTGGbRY9eWxu5hIWQXECH5cz9MW/baSR9fmosmheNpM2VSHkU6X1xYSVpw2jXG1drD9tw+h0WWm3nfKSBiGtpd8ENaZNktepwCItIiyzen26vJpzAIN0ZakDktsuzaZ2EycoESLkxCkqbmxK/0HYOrO6Rcv5RqHVOTVGqa1O4CIiTpeLl97E9CkNXRcmbb3QLXx4vhy06nOmWV5k/Tv4dj5K5hN9+vO4cdtHNJmJddY2GOk/zcHuuXUkhWmM92rF9i37A9G1SRaqtKnRLrl2HjSj6w6ErXNFItDMLVHoU6Zd3ffCahTnGXZbqP2ykqpn2FZpcRkg6baQ9ZLvaeL/lzK81nF3++mLZU04R6aWo+64NCGMRhqFNSkIOrPW6GmubQ9vWpOS7Iat2+PUYOwxpe06Taw/U/b/1yeLDf4Ih7/brrj2p9it39c1cr/JGkfy7pfxxtrwEAAAAAAIBr5F7/m+21kl6n9z4vf9Hmn18v6c4Zr1tJepukt0p6Ve/99TvuJwAAAAAAAHBtnfmQrvf+Lkkvvvvvrb3ne39/rvf+pj3uFwAAAAAAAPDAuIh0VwAAAAAAAABnGL6q7Nqf2Pz9a2PvCAAAAAAAAPCgKj2k673/X/vaEZzfk6aP69GQ0nLXsvClyZRitgppaq5/GmPWtpPJ5iHFbBISWycmsXUa0xCHJ/nEFFbTvgxpgavV9jFamgRBSWohtcm1p6AiH6G6exJsDGV0m6ukoxX6Sj6ZrJSwVkhxXbcP79tsEuzwvtLu6a4pHtSmo4X+hUukzJ63lcSz0FcpCW1pElRT34Pta7Uc/OUjTX3fyrUapuula4+pd6HdGOUcGGMMe4yGT0aTkJoW011N+zTM1zNzj3WJqFJON3MJaSmZ9ZYZ48YkpaaNke5aSE0zH/YkTOLTeILuJtU0qU5Zmusy1SlzU2PkmiYcT9N+HI6RS5ar1C6StHLp0WHuW5q+7lpIfaVQ65QKh/T+CvXLCHNO6TCnvu7QxZpmeHtOZnVtoW8aYz68tmoLk+ZcqVMkydVyqX4p1Tp+c7uKNXeqX8y9PNX4o6S7uiTXcM71SmJrqXAIc4MLoQ/35pV5f9WfE+zPGoWpqHwKubddSLLPifXh/mDa08+8MzOPz8K4KVne1Sk3D3zt4eqUagq9GyPdV2+aGqhyD5Z8rZL6Oqdrmsn0ZPDrzuNCft21rT35IrYFAAAAAAAA3G+qv+76Hq21T5L0MkkfJulDJB2abjNJT5L0iNYPq8+9PQAAAAAAAOC6Kj80a609S9J3SfqY082bv5/4/cvdf+8OAAAAAAAAuOZKD+laa7ckvUbrb8+dfgB39+HcEx/KdUlvl/QGSb9+zn0EAAAAAAAArrXqmnRfLOm5m39+jdZprx+t9z6cu9F7n0i6JemzJf3W5p//fu/9r+y+uwAAAAAAAMD1U/1118/Z/P3Tkj6z937SWrtx6r8/SdI7e+/Hkr6vtfaLkl4v6ZWttY/svb9r913GE92anOihU2kxq8Kz16VJQkspJyndbNJN8pNpS6ahr0srS1Ki24FJcTkIiTbTkDLr9i8lALoUnsnE71sKdrTtMTnogn+jfNcktPT69PZKCbEuPWz4uHGMlI7m0spiwlrYoE2vDclPZoyYbLYK2ysELcaEWMMlm0lSs2mkoa8ZIyaehTE03d5gmg1X9mDsfq3GxDqTPGuPz3rn/Bj2HA/Hs5B4lj7qUZLQ7MC797XTZHiDMRzb9E+pae5ektLRDlJCrOl/wySpST7JtZLiKu2ehDZrflyXPHoY+iaTwmRUqWlObJShrydcSq0kTU36XqpTxrAyZ6hrk6RFqHXcubg0CaySP8fHuHbiZT3C9b6reK/cW5J9SlgvjJHSViu1Trgs3dguxVWq1TquTlnv3/Bax96QYmK9b961fomnYZiKmqnz448wpk5JnVuogdycGMsJmzDq+7p6olIvV8cY4+eEnfuOIU6Ubq4NXUO7S8Z1P2tK/mfT9PNqanc/I7v6J7Xvs05xdVE13dW1pxT6aTxx1+bhPY2l+k26D9f61P/23vuJJG0eyN19+PZ+pzv33n9J0jdJ+lBJf3G3XQUAAAAAAACup+pDursP4d74hPbf3vz9fPOa12z+/rzitgAAAAAAAIAHQvUh3d3whyf+PsFbN39/vHnN2zd/f3hxWwAAAAAAAMADofqQ7lc3fz/3Ce0/pvVvZH9ha+0pT/hvH7X5+4IXzwIAAAAAAADuD9XgiJ+X9Klaf2PuW0+1v0rS35f0AZL+aWvtZb33d7fWPlLSV2m9jt2/HWF/YUxbf5+FgKcyCwinFUxdc1yfdfgYKcjBLdjowiskvzh1ap+ExR1de2XcfUoLNV+4XXdjjMVcxzgUbqHnwqLQUljcNoYzmAVhq4sbF8aw+1wMiLBjFxZYTuL/gTHrtob8GS+t2ZpW2LXvLyzIbA9F4diH9lYIvKksWn5W+7VROHSV+bN2LwnnQGV74YN17bFvYd/S4sb+/hdCLcwK82l7LiQiBUGkII6KSk2TAixOTLk7k++7NCfiJJyc+dgP//xK50XhZpnO2dK+FT6/Uk1zzf+3fXkOrxw6e99J967h4QylfYvjhjFcyEBl31LIRFKoa+wxCjVGDK0y858Lk1gPYrYXap0cyrU9Rsy2MXNlDIMo1NH7qtvL417kzzB7dNE/E16Vn7HdPS2FPriwp9w31CSmPd7Hbd/37sNBSvcZSfWbdN+8+fvPtNY+4G5j7/1Nkr5b69vun5b0ztba2yT9jKRnbLq9esd9BQAAAAAAAK6l0jfpeu+/2Fr7Hkm3JT0q6XdP/ecvlvSxkp4laSbpqXrv/yv7eUlfs/PeAgAAAAAAANdQ9Zt06r1/bu/9Zb33f/+E9ndIeqGkb5d0rPUDutuS/hdJL+q9z0fYXwAAAAAAAODaqa5Jd6bNg7qXtdb+gqSnSHpH7/G35AEAAAAAAABo5Id0d/Xel5J+Zx9jAwAAAAAAANfNXh7S4WIte8vprRur8JvNLll1WfwtaJfk6tLR0vZiamzc3nb/lCbr2t3rL0O/Ivuxs8rb2OdbNmOnQxzbzWmUgkQ12f4P6XvDLolLkmSSwuJpYZKRWtq5sD2XTDbGR5ISyErcGOYYr/vuvrnB+6D8/lx76bKunp/XZMqI9hRuNsacX5mv0z3Ntce+9j43/L667u8mNJ9GZu/7YUKzfdPh2dMvU1RqGskfu0qtk2qMyvbS51c6LwqT3xg1xt7qpSuSqLgv5Tm8Egpua510jwqbs+0jfCiFeiluzqSfVmuMXc/auL00ZRRqATt2cXt2jEK9FKatcG75vvurw/bc/4q6Kj8TXvTP2O6elu6rLk027VtKbPW1Ttg591ZO1TT3evayq62HdK21Z+5rY733X9vX2AAAAAAAAMD9yn2T7i3az//j6mF7AAAAAAAAwAMtPTS7Gt+5BAAAAAAAAB4A7iHd63TtV4sAAAAAAAAAro6th3S99xdfwn5gB8d9pqN+eGafuOC0WUAxLpAcvmA579PttpX/kuZx324/Ws183zDGwmxvEfbZtS9W269ft4f3bdpdm+QX/kyLgZq1cc9s950LfZN9BT+4vnFRYd/sPta84K059ilAIW7PjeH7utNosvLbS4ucT2QWNvXrukvNjBG2l9rdEL10wgWFRZ37tBDOEBZCdp9TGjttb3Uw/LM2U07sH/etcC7H66RyTV2FBaCTEea40lwbhnbzeLoPzMx5FO87YYx52z6RDsLixgem7yRcq6ndTS/JykwObpFmydcN02JAxMQE4SSVhahTGISrU1Kt4/qmOsX1Te2xr7mZuDYpn1vuXEwhE26MWNPY1tr1N8b1vjeVgIc9zctSsdaxdUqoddK9y1yu7p4opfIs7NwyfIDmeu9hDmjm4Lc0brBzXVMM4rBBYpX6JdUp03BdunIw1ValfRvWJp0RVlKpdVz7GHXKRQfaFU636s9+lbm28vNqvJeY2mMx8fd3+zN2mHSO088rTjpf3G6kqSh8sFNzoFOt4+qa07XL8Z7vV7UYTwAAAAAAAACj4yEdAAAAAAAAcMl4SAcAAAAAAABcMh7SAQAAAAAAAJeMh3QAAAAAAADAJfPxmbiv3F7eVFu+N901JZo4Lt1sWUg8k3yCWEomc4mtadyU7npi2h9f+uS1k+X22Cdh39I+z3dMy1mtKvFh8ilBsa+LcPRdK0qbGyOZtTJGIa2skpYkSSt78ELimWl2CYlSTibrk+3+Jlhp3W6SkWLiWXjjhUDFUfiEvHBylZJSh6e7psQzt73VLCWsFRJi0/lp9y31HT5G9ZqyfYtTlDVKQtrwSFp3Dae5dhWTu7cnjZR4drzYvu9MwmRbSS5NXKLpYjo8SVSSbkwWg/vO2nas9DQknrm+SUpN21VKZk3c+67UOmPUKcfLkHpv+h6Fvm5cSZqbWseds1JKd61dO+5ay+mulajpgsIQ1TnONdfqlDBvpQR409zCZ+KSS1OSb+TC1Oe+69Lcbyeh9kg1ia1fwrGwQdEHV7imkWxdE88Bm7YatpfG2DHJPkwjYdzUd3j7hf+cMEZNk+zrZ7T0c4mZB1KdMjXXlPsZdt3Xf4Cufpksd7+PL9P7c7VO2Of5ZHufZ92fzKk+q9Q696pf3m3uu2Pim3QAAAAAAADAJeMhHQAAAAAAAHDJeEgHAAAAAAAAXLKdHtK11p7SWnu00P+RXbYHAAAAAAAAXEelh3SttU9vrX1za+3NrbUTSb8j6fdba8ettTe01r6utfYxZwzxNa21X26tfeUuOw0AAAAAAABcJ4PSXVtrnyzpf5L0gtPNp/55JumjNn9e3lp7raSX997f9IShPkrScyS98Nx7jC2/O39YR/MbZ/ZZFp7HuqQV6Yx0V5fMEvq6BNXjkDCzCO0unTUlnrl0s5Tu6pJgJenEjbHwfZdLk5pm2qScmtZdotd+AvKymGDlIshC+lQlCTYE5PhzcXicVy/smySZsNXSPrelHzilmLmQxNjXnQPVFFeXsLbHdDR7nFMyqz1fhiezpv7x3KokpcbUtGHjpr4hlMqmsa3Hdvvmx6gkrMXbw67h0THxLKUSm6aUcOjSnEPfNAfPU5TyjlLS4sJMMPE+Z5JcU5KoS3GVpMcn2xPMLKSVHZi+03DjselvYSJJY+wq1TSpfnHtaQz3mczDxZM+P/dZ5dpjeN+UQm9rnVjTmPcX+qZrx11rNpUz9GIV6bQAACAASURBVM0Jh6l9eCyjzWhPL99Tsny4jcf34ZIIU4ayq2smk7DBVAO5dPNw73JJrr1Y67hpINY6rrnSdwTxdCvUL5VjnxPr/W7sXOsU+sYk2Fi/uLbh7y/WYfHnEteWrocwxr4U6pSYQm+uvxbe38nwPYvsvTLs26Jv1x4pmfVw4k8YV6fcMG2SNDPtByFtPtUkLt019XVO1zR35oMeo53bPZ/ctNb+hqTXav2Arp36s9T6m3Tv1PqeePq/faqkN7TW/tSpcaaSnrfp+0OjvgsAAAAAAADgPnbmQ7rW2l+T9A83/Zqkt0r625I+RtKN3vvTeu9PlXSo9bfj/rakX9/0fVjSq1prf2Iz3EdIurn55x8e+X0AAAAAAAAA9634kK619vFa/4prl7SQ9HckPbf3/ore+xt7f+8XvHvvy977G3rvr5D0YVo/rFto/Wuw395au6n3/qrsL/Xe37KXdwMAAAAAAADch876Zdp/pPVDvBNJn9t7/8EhA/beF5Je0Vr7WUnfK+npkv6qpKdtuvAtOgAAAAAAAOAU+5CutfaJkj5R62/R/fdDH9Cd1nv/odba39L612X/pqS3bMbjId3I3jF/RLdObp7ZJy6mXFhFMy1Y7IMjhi/qnEIf0uKTbgHLuMiyWQw5jTtPi4ubMdJiyguzIHMMjggL73bXPyw0aldILi6k606NtMhyZTFltx/VJcTdkVsVwiDSYspx0VzzmYT1S7WyoQ++bwufn1s4OfYtfNatsDB/kt6LHbaSE57OLfO5xnHj+VkIgygsnF0JD4nnlguOSCETabFouyBz6FtZ+DxdJ+59j7Dwcgv3I7vwfFgw3M2TYUot6WnfCospz0P7bLL9BmfTsECyCRM4MK+XpMMQHOH6pwWSD8wFPwmTwNQFRxRvPJWFmlP9YvuGE3RpxlilMAjTnvYh1ROurom1hwuqCEEO6dxauFpnhJpmuUiLmW+329pF8jf+eJ8bY4LZbsr3gXAeusX2UyaCmYPTrSveVs39r6V9c3N7mPzSvcTVGS4gQpJWpm+udYbXHnEMd6CLde1VqF9y3+H7EAMlSiFnZu6rBE4VQiYkX6dU6qJY04wRkuWOfXXKqfzc5QJ20s+E8YKonKDbCjlyknxdM58MD3JMARHTUL+4OiXVNK5ucLVL6iv58InE1TWnx318Phs81nmkT/5u4MMvS/raHcb/Wq0fzn2g1mvWvVvS/73DeHvVWntqa+3zW2vf0lp7c2vtsdbanxz7NQAAAAAAAMBp6dddX6z1w9fvPL32XFXvfd5ae6Wkr9iM95re+/y84+1La+1hSX9P0su1DsGQpN+Q9CZJvzfWawAAAAAAAAAnPaT7kM3fPznCNv7VqX++cr/q2lr7YK336/mS/kDS10j6Lkk/27v/cvt5XgMAAAAAAAAk6SHdUzd//+YI2/itU//84yOMN5rW2kOSvk/rh23/TtJn997/v7FfAwAAAAAAAJwlPaR7TNIHSHp0hG3cXVVvIemXRhhvTF8p6WMlvVXSJ/Te37Wn1wAAAAAAAABRekj3Nq0f0j1L7/vrqufxrM3fb79Kvwq6+ZXVL93865cMedh2ntdchLcfP6IbR7fe8++VJDQnJsEWUu8q6a4pIc+lsaX+aQyXbhbHTelmpn2V0nnMfuQU1/A5FZLQxkh3HTqsJJ8qVkhHS3LimRvXd3XhQ+kYtxDu45LJVilwyYwdE8VSCqtLTYt9K4lnhYTYxPUdIz0zjWFTt0LnSvpwSjxzY6RwwpiKapK7Yl/TVkgrS+1jJM+WUnRHOAfieevaw9zX7TUc5vC0OZf4Gbbn2uchDXEaktAOTJLrtJCE5lJVpZz66sZIiWcuxaylvi7dtTS57E+lfklJvi4hNo0baw9Xp5RqmpAmG1Jf3RjpXK7VNOH6M6mvuaZxk7vvOkb9Upm34tzn5s/CvqWSLQU1uvnM1RiST7+s1DSST61P55y7tKtJ9jbdNabQD6+LStIYhXta6cerSgp99fysJMRW6oZKTRNTX01jobYq10XmZ5D4/nxzzY7prjmxPt2PXPEY5mWzveXUX6zLaUj/Nv1j7eH6FuqRNHbqW0l3vYj65fhov+mu6SHdGyX9EUmfIek7dtzG3aTTt+w4zti+QNJNST8v6Qdba/+ZpM/R+qHib0j6GUnf0Xv//R1fAwAAAAAAAJwpPaT7IUl/VtJntdY+tPf+1vMM3lp7mqSXav2M+bXn28W9+YzN3z8l6Qclfbrp89+11r6w9/6aHV5T0lp7bEi384wNAAAAAACA0fxqS1/h2+i9D15KLn2R9gclvV3rb419U2vpC9pZa61J+npJj2yavrs6xp597Obvz5f0Qkl/SdLTNn9eqvWacx8i6Xtba8/Y4TUAAAAAAADAmezDt977bUlfrfU3tj5N0re11g6HDrp5qPcNkl6i9bfofrj3/nO77+44Ngmtd59kHkn64733b+m9//bmz6sk/TGtAzQekfQPzvOa8+xb7/3Re/2R9Jxzv3kAAAAAAACM4TkDnuEMln7dVZK+UdJ/KelTJP0ZSS9orX35vX6Ns7X2QklfJ+njNk2PSfryyk5dgCef+ud/0Hv/lSd26L2/rbX2P0v6CkkvPudrLsQ7bj+s2eyhc722EjIRF1l269kXQibSYspxe2ZhTBfYsB5j2OslqYcxXH+3OOe6s+lbWUxZklz/QvBAbWVb+RV5zeKxSVzM1S7OGY59XMR2uy0tWNzdYsopyKGw6HElDCIuphwXlXV9hy+mnOwcEDGWXRdkTospVwIlKiETaX8L52flXC4FNqgaHFFYTDlew8PaUns6D9Mp7q6/ynQWL53wxt2C0ynoZ2nCICZhceMWwnQmJiQiLm5sxojbs61+MeS0PdeebgP3+G2Oe+7DGKoBWTYkJJ4vhZCJVL+4vpUao1D/SL4miXWROY1STRNrHXedFOqiSsCA5OeSSghRvEjC3Od2JIdBDK8b4u8lFeoX114Nclja2mP49qrBH7auCR/glahfimX0roESsaYpBF+NEjLh2it9Vds3W3vEkInwYRf22bZX71Eu2CT9XOLGNqE7684hUMIG2vm+KzcXpVDEUKc0FziV+prPOtUpeQwTHFGoU3JAhG0etX6Z39nv6mPxIV3vfdFae6mk12kdIvGHJb26tfbLkn5E0usl/bakE0lPkfQ8Sf+F3vsroU3rULWX9d7fvLd3cD6/o/V+H0r6f8/o9zObv5+x6V96TWvt/QiRAAAAAAAAwL2c9U069d5/r7X2xyV9m9a/uipJz5X0JWe87O5jxcckvbT3/qO77uTYeu/LzcPG/1jSM8/outj8vZR0W1L1NUc77ioAAAAAAAAeAPcMhOi93+69v1Trtel+SuuHcGf9mUv6Zkl/6Co+oDvlFzZ/v+iMPn9k8/cv9d6PzvkaAAAAAAAA4ExnfpPutN77q7X+dddna/1rrc+X9FSt12p7p6S3SfpxSa/uvT82+p6O79u0/nbgS1trX9l7/8XT/3ETFPFFm3/9rh1eAwAAAAAAAJxp8EO6u3rvb9H6m3L3ux/Sem29Pynpx1prL5H0b3rvvbX2gZK+U9LTJb1F0j/c4TUAAAAAAADAmcoP6apaa8/V5ldAe+/fv+/tDbV5sPZfSfpBSZ8g6V9Lentr7W1a7+9E64CIl/beHz/vay7CH9y+pen0fOmuKQ2vNsjwZLLS5mLqq2mrJIWlcStj7KtvaI9JaCl5tKISEOv6hnQ0l2K2Cuk+Ob3WtVU+6+Hjpv4x3KdwMlfGaKU4y+Fdz9V/H3ZMQZN2T1IbbQybFFbYXiEdbb09d8Kkvrtvz6ckjnASpWvYnKAtvEF774pJlCE1zSb8Dk+QW4aoxhTUJ/P5xaQxl5YbP6eQhObaKp9foW/ctwtWqmlGSLePm7N1Ueo7fHvpHLdjxzF27Jv6V1LoY/L6CCdSIbE+fibu7cXEVvcBhr7h/dlLrZIWP0Ktk1S2l8cYfr5YV6F2kUapX5zSPfis/gPHKNX4xRR6P8bwe3CpplF4L5XjNsa9K12rZvB0LOKHYm6uPdQ0lWOf06oLPwi5OqWQBHvm2MM2d0YNtfukca+6Znl7vxPTPdekG8GflvQvJH3PBWyrpPf+e1qvL/dXJP0rrZNbP0LSmyS9QtILe+//btfXAAAAAAAAAGfZ+zfpTrki/5/1ffXeF5L+6ebP3l4DAAAAAAAAJBfxTToAAAAAAAAAZ+AhHQAAAAAAAHDJeEgHAAAAAAAAXLKLXJMOe3Ly2KGmurH/DY2SBDvGGDsubzhCImZM4BwjraqQ0jZG6pYbOqXz2N0oHIuUBFtKuq28wX0ms+447rn6P8EYQZtX2a6XuqS9pbFJIySFjTCNlBLWRkjXKr2/4YFg6+4ubizEL7alm7jCwCmiyx6j3Y/F3gKaqxf8rtfPlVxJeEQXXdOMMqGlsQe2BeWaZtfao1LTVMcuJBwmLl03hDmHZN04sm21R2OE2mNv9csI1w71ywAXnSZbGXeM+9wIaau+Diuknxa3ZxOTQ43h5oaYuFw6ziPUNJUpf5TPeoxJY/chxqxrlo/NxxvM4Jt0AAAAAAAAwCXjIR0AAAAAAABwyXhIBwAAAAAAAFwyuyZda+2ZI27jA0YcCwAAAAAAALh2UnDEWzTOkrq4ANPHpprGFfmxq32u9Xwl7LyQagiZONfO7OC6f064fq77Xba4CHxYSn6cfRmopb247p8VHhjXpqYZIbgl277gqWmAAa7LvXKM4JZC675Q0+zRu/f77OWsdFduDwAAAAAAAMAFSA/pXqfr8ywcAAAAAAAAuNLsQ7re+4sveD8AAAAAAACABxbprgAAAAAAAMAl4yEdAAAAAAAAcMnOCo7AfeLgXRMdLC/geWslSqTQt5zEZSJpSmOkvpWksMIYcXHHyhixrxl9hMSznk6nwvZse9y3QjRS6bMuLq1Z+qwLY+/p2slj3IdLio4Ryefedimhq7gPle2VLsDh7S2NW9q31L49djy1KpdDZYzw/tpq+D7EI7/j+VI+FvbzG769eCxGSKEruQ+nl53n1erUYO8lww9cvucPb99rXVTp6+qJyj1fu7+XUg2V9qN07IuT0a51zb7GjX1H2N4Y+3FV7Fq/VN/yrveNEe7jV6KmSe2VmsbVEmewp+c+a6jC5V45FrU6LG2wMu7w82WvNc2I00u7vd9nLzs9pGut/djArr/de/+CXbYFAAAAAAAAXFf2IV1r7VmSvk/SKyV9Xe/98fD6F+vezySbpN5a+297779x3h0FAAAAAAAArqv0Pb1Pk/Q8SX9D0smAcdoZf+76T8+/mwAAAAAAAMD1lR7SvVjrb8h9f+99OWCcT5X0HPPnwyT9202fj9tpTwEAAAAAAIBrKq1J98Gbv//NwHF+o/f+791/aK39c0kfeWpMjOzwXdLB4pwvHmFh4b0t6BseIfdmXjBC6EMKS/CLHg/v2yfhN8JH2J4m2/8hbq8U8DB8n0vvrzDuut30D9trrj28Z9tXkj21wj7b9tg37IfpX+u7+wqoaXtj6IXd64WJJPV126v0TZPRKi5kbK6/sAix3Y9V2LfQLjN27jt8segWxzAhPaFv4XKI79st4JzHGPb6s9vdZxL6uv0obq/St7Y49e5j7G2h5n2tAV8OnBretVQ3lO6r4Xp3l2pxPWrbP9ZQhb7hXulrHT+GO7dsHXfWGGY/CtNyrBtK7zvWYcPrlLQ9W5PEOsXUApWwL0mTHeuXXBeF3ShMMLW+g7uWVGqXPEaqPUzdELY3Rv3it1eooVJ9kPbZ3leH1ynd1B3r/zC8LqrUNPEkqgQ5FGq2mCmT6obC+ysFalXql33WGJW6b5/BXufoO7lTGOsc0rT+QZu/3zbCNn5h8/fTRhgLAAAAAAAAuHbSQ7q7D9TGeEj39s3fTx9hLAAAAAAAAODaSb/uOtv8nVJd7zrR+ouKZ30h8Pc3f3/QGX0AAAAAAACAB1Z6SPd2Sc+Q9JSzXtx7vzlgG0/e/H1c2C8AAAAAAADggZF+3fXur6g+a4Rt3P01198bYSwAAAAAAADg2knfpPs5SR8t6cWSvm3Hbfwnm79/c8dxEBw+1jWbny+KaNdk1jhGJUG1nLbqEk0LY4S+q6lvd/1TX7+9kKIUxlhNzWeZtucaY9xVSGmrJL25fUspZvZ9hESw1G7GnoTtTabbEUGx78THGk3NGNOUjmbGCB+1pmF7LsXMpq6dsR923yqJtHtUSWx1CaopVTWGRJn+y5U/mZcuVSxsL4+x3b4KqVuufbX048b3vTT7HMawCWvm9esxbLOdm+MnWvisY7qZS00L++baW3h/pTHSvrm+KR0tba+QSOsT3cJ1XUlvqyTWjdH3glXSVmvJrGF7hQRVV7tICqnpxe2ZGqHUN72/6fD6xdYHoW+qaWJipGsspK2mfYtfVajUL2Z7zdQS6/bhaauuppF8XZNqnVR7uPZUH0zN2GncWHsU+qZ2O+4VqWliAryxrNQ6hZokhaKuTF+3D5K0NPWEe/26PRwLO8bwJHtX56x3bnh7+jiaORNTsm6p1ikkpab6YFKpXwrbm6T6Loxh+1feX+wb5oZdE2mrabJjJMRuTI72O/+k29NrNn9/WmttyK+0Wq21qaQ/pfUh+dfnHQcAAAAAAAC4zs56SLeQ9IGS/uIO43++pGdv/vl1O4wDAAAAAAAAXFv2IV3v/Xck/TOtv+351a21D6sO3Fr7YElfs/nX35D0w+fdSQAAAAAAAOA6S9+kk6RXSJpLelTSD7TWnjl00Nba0yR9v6QP0vpXXf9x732xy44CAAAAAAAA11V8SNd7/xVJ/43W36b7jyT9dGvtz7UWV6VXW/s8Sa+X9AKtH9D9hKR/MupeAwAAAAAAANdISneVJPXev6m19pGSvkjr9em+TdJXtdZ+QNJPSfqdTdcPkvRCSZ8h6Tmbtqb1r7l+Xu8xOw4juPGulWbHKSJuo5DiWko8U0hbLSShVVJcJZ/SVUkxS8msk0KKWQtXjt1euspS2JhLforxNcPHjRk0Jr6mp8RW134QEnsOTHpYSDZzqaqSTzc7MONKPm3sYOqnnlnY3oEZYxbGsH1DjJLrK0kT86kchDFc4tlBiFGqpKNNYrzkcKt0Adq+IRHMnMyp7yKljZn9WIR9my+3L9bUN6W7zk27G1eSFibxLL2PxcKP4ZLXlkv/Wa8W2317/H9s6fMbnmJmxRjesDWXoLrw+zwx382vpLiOMcZkEea+SrprTHpzry+ko6Xtpc/PRNyNku66r/CzStp86J/77pggH9pTX1eTpPT3Uq1zUElmDfuWah1zHtlEaaXUyZAE6zd39u/9PFEhsT4ltqpQv0xM31TTpPaDg+0LPtUNB2aMVKfMCumuqZ5wY+faY3htFdNdzVlQqWnO2o+hKjXNuv/2SZfqCds3/GDi6qJ1f1N7hDFs31Sn2L6hTgljuHppYeoRySfBLichsT79TNhM//C7ey7JNZ8qIW21cP/zdcPuKfSudhmj77rdXH+Ffaun3pvao1LTxCTY4T8Ml6aXU30n93r2sqMzH9JJUu/9r7bWfk3SV2t9xj5L0svPeMnds+8Nkj6r9/6bO+8lAAAAAAAAcI0N+l8Fvff/QdJHS3q11g/hzvrz61p/8+7jeUAHAAAAAAAA3Ns9v0l3V+/9ZyV92iZA4jMlPV/S0zdjvE3Sr2r9EO/1vbsvwgMAAAAAAABwBj+ku6v3/muSvnEP+wIAAAAAAAA8kMoP6XD1HP7BXIeH8/f8e14c3KgsppwW7XTBEYUwiLhAcghy8IssD++bFlNehQCEbq6StkxjuAFsV7uAqeTXVG/h2LsVmYvr3dYWWTbteTFls/hoWtx4FgIXzALJh2aB5XX79uqoh2F7N6Z+JVXX/zCsunpoVlKN44YxZmYV1BQcMTVnRlpY2vWVaosvTwsLLy8LCyQvw8nl+qbgiHmYNBam/Xjpb3Nzs88nYZX0k7BA8tFytj1uWLz5eLE9dhr3JHyubrHn+Tws3mzbwoLMpclo9xSiOEIhWKGyQPJk7tub6Z/CIGzIRNpeancLMsdFnYcvphwDLNzCyXFBZrO9EEpip5FiQERL55wbeseaZj2GaUthUe5+Wwy4cnVGCmRy+5ECG1ap1jH93T049U11UTrnVttTX8yUcR9fqlNaqgftXFII1Er3vlC/uLpmmkKrTE2SapoU8ODqmli/mFonhVbdnPrJr1LruDol1ToHYcKemfYYHGHaqzXNGPWLk2oSd29NwVC+rz/xXZ0i+Vol1TonZuxKrRPHTfWLq3UmwwO1FGqaVZgcVuZekkJs7JRRuL2kMXLQwfbglSArKdQpoaZxY8T6oFC/xLrIBUeEuiGHT7jaI/QthUykuna3kKz3cRIO4kiqP84DAAAAAAAAGBkP6QAAAAAAAIBLxkM6AAAAAAAA4JLxkA4AAAAAAAC4ZDykAwAAAAAAAC4Z6a7XwMFjJzqYHZ/dqZCOltLKaumuoe/B9nPhmMwa0sZsEtps+BgxNS21m7FjakyMxjXjhnb3UbWULOfCjoopeyUm8ayZBFZJmpiEypR4lhJbb8y2k3NummQzyaeN3TrwEUg3Q/stk4Tm2iSfhPbQ5MT2TclrN0y80iwmrJkEuRCBVElCS33H4FLMcmLr9gWYUlwr7UcuhlDSsUtHC5PA48vD0L499lFIQjuabvd9fOH3bTrxY7T5bhd3evUqJKHZuT3FYJmJq5LiKtWS0Fx7TCtLSWjmeJZS08LnUUmITWO4xLPJPFzXKYXVpZiFfbOJZ0u/PZ/uWowx31WqaUK7T3cNN1Z33qd6JKW72jolpCuXapqU2Dp8DHvthJomzdeuufKRtJSsm2orM+WncsvOUSlNNiTuTkySq0txlXxd49LmJemmqWkkX9fE+sXUJDdD2mqqX25Nt2uVG6H2cO03w0SZapJK/TIzE3lMcd1T/VJJoU/9K3XKPNQex6YuWvffHuPOytcpLsnV1S6p/SjUKal+OTa1f6ppjufDH0PM01zk5oaU7rqsRrm6DQ5sU5hrU01TqF9SnTKt1DSxfnFtw2udSk2TxojJ8ovtcyvdu0rprqmvc6qm6fN7PHvZEd+kAwAAAAAAAC4ZD+kAAAAAAACAS8ZDOgAAAAAAAOCS8ZAOAAAAAAAAuGQER1wDk8eONDGBDO+jsqJvDI4Yvuhx7OuCI8K+r0J7N4svrxZpDLOY8mFY+DUsNOrXq02Lj5rO6dgXwiBaXGjdtPmumXtBXBzeBUeEBX3NwssHIWRiFhZkdospPzTz4QwPHWy3uzZJeiS03zLBD2nhZRcS8dDULyJ6s4UFoM2KrqmvX3jZrzQ7DSfM1Jyfkz0GR6zMSX4SFlOu9E1hEC584s7qhh/D9V36vnfCOeAWWb698GMcmMWU02LYqb2im0WWY0BEWsR9UVgdviBOL24+S3OfmTLCuucxLGHXBZmnJ2nc4e0xDMItkBz6usWUJR/8kPrKLNTcVqmvab/iwREyIRE91Cky10OqU1L4hKvJXO0i+Von1jQpUMJMiWHt7VLAVQqKmdjAMD+GC76K1/UYp0uhpqkEXx2YmkbyIREpIOKhWQhyMCERqX552IQ+PHzga49x6hcXkuX7pprE1TqHYQV9V5NU+o7B1SOStAwV9tKc/K4eWbebgKvua5pU67j+sX4xgRKPT0IYxHS777snPpDiIASiTRe+v+OCOFI4xypcf7Z/Ckhy80CxprHTS+VntNA3HE5b17h6RPL1S7VOseETaQwbHBHqkUqtk2oa1x7qlBQ+YfvHMEjTfqptsjjyrxsJ36QDAAAAAAAALhkP6QAAAAAAAIBLxkM6AAAAAAAA4JLxkA4AAAAAAAC4ZDykAwAAAAAAAC4Z6a7XQLv9uJoPQbw3l24WE898ezP9+4HfoWbaU2pam4UUyOX2aZvS+9qhGTtFkBVSxVI4WjcpQZOJHzglKroQq/j52h0pxqPZqKLQ1byXlEQ5NalpB1MfX3QY2m9MTWpaSiszSWhPColnKQnt0el2Uo9LQZOkR0zfh0PiWU533R479jXtLvFV8imuqX2MJNHEpW6dhP83tDLXZUo8O5mE1FeXeBbTXU0ya/is3728adsfa9vtYyS2pnSzxWr7GC1NmyQtTALnJMztYQj/v/EKQWjxLRfa0xguCS1cDjH1tZKa5hLSUuLZ9CSkmJn2mO56sv1m2jy8wZSENjdvcOnHcEmwKTWtlO6axthVqlMK6a6udpF8Ymub+jmnz0IZbdNdQ110uN2ekulaulhTXTNUOG6xXDJvJaUT2pDL4txQuk25t1JInpV83TY1ia+SNDO1TqWmkXz98mihfnkkJLO6mmbdf7f6JaW7pvrFpbOmvhMTf3kYUlz3Vb+ke3BMdzXtLsVV8rVHTnf1Sam71i93UmKrqxuKP1O4YxfrFFvThGMcxliaJOxV+rmr8PNOSThfbLprmCdz/WKS11NNYxPkh9c0kq9rKjVNSmadnPiddnVNTKFfuL7hwFXql9T3HvVLWz5+5n/fFd+kAwAAAAAAAC4ZD+kAAAAAAACAS8ZDOgAAAAAAAOCS8ZAOAAAAAAAAuGQ8pAMAAAAAAAAuGemu10C/c0c9JNm8xyTGdW63hRRXhXQzl5DWDlLimRkjpKOldDMX8tRvnDfe9vQGQ9KbOXaVxLMewmFiwo9725XEszHCrgrprq5NkiYmCc2loK3bQxLagUt39QlBj5h0tFtTn3ZVSTx7dOL7PjrdTvVJ6WgPN78fLt3sZohtumFOmEMXHSVpFj4/d9qOcOUonMo2jy1d1nNz0h11n/42Dxfg7b6dWJYS5O707dTXSUiQm8ZkueHJlS7xbBHeh0s8S+3zpf8EXRJhSpqO0Yl7TP61m3NJaOkQ29S0kIgZ09RMElpICvepacMTz9bt2zvi2iSpmSS0mO7qUlwVUs9SEtpilmij8QAAIABJREFUe4zuUtAkaWXGCKmx0apwbqX6xQl1ikwidEuJ9a5+CYn18X27uiYkxU1cMm5Ihkx8amGqabbbVmFzbZlS6M0+h1TG0nW9L8U5rpLuemDaU4prSqe/ZdpT/eKSXJ98cMf2TfWLS2d9UqiLHmrbfW9O/PtItc7MfOA3w8Q8M0XstFDTSLvXNZWaRpJcgOZJuL8fm7rmyEYgS0fheLpaZxZqR1e/7Fq7SNIq/KCwMJPJIsxnJ65vqH/m4X5U+bkkvJWaQgp9pW+sU8zbTknaLrE+ptvH1FdzvqSaxiWzHvudS+murn6xyfSSr19M7SJJPaa7uvol9HU/wJ+qXfqKdFcAAAAAAADgWuMhHQAAAAAAAHDJeEgHAAAAAAAAXDIe0gEAAAAAAACXjOCIa2B1dKxVO7UaZmGR5eYWWW5h8d+wIHMvLJDc+sy0hQU+bat/shzWGVUz76WFFWjTouNuIfFVuHLswp9hMeUYBlFZUHlf67oXFoxvoa9712nh5cOwCuqhWfE0Lcjs+j4UFl5+aDK8PYZBmPa0SHMOjth+32kx5RvmgB6Ga3UWFvSdmP7T1Ndcaau4bLK3NCfoPCSpzM08kEI0jsLqv1OzGH06t+xbKa42PTeLPc/DuTU382daTPl46ScYt0C5a5OkqVk4OQVHpDuGO71GmXIKiyxX5sk0d6Z2d2qMETLRFmmR5e32FAZh21NAxLFfxN2GQcxDX7Nwcg8LMtuFl8M9v4f7+67cvV2SZIKsJPlAiRRwZd53m23XLpLUwrXq3nWszMx7cXO15IOsJGli6po+TUEqZntpmixcU7F2qQRc7ammSadL+kx8cITfucq8nO5HPjjCX6uuTkn3ShcQIfmQCBcQIUlPMnVN3l4IzDB1or+ipJn5sGYhBCXVL46raSRf17jaRZJWYT6bm/d3EmqdQ9M+C31nPSzCnxLpjJUJsFimY2H6utpFkuahfnE1uguIkKTDyfb8eWJCfqSz6pTCpLGvMKzKLlTmSakWkmVq4FinxKAtE8QRahpXp1QCIiSpnZh5LgRc2fol1Slhe93VL+k5xD0CsVbd1/tj4Zt0AAAAAAAAwCXjIR0AAAAAAABwyXhIBwAAAAAAAFwyHtIBAAAAAAAAl4yHdAAAAAAAAMAlI931GuiLhXoLiW13hWQkm9IWknxS+snwbKVa4lkLKW19ut0ek2emJk1vEdJdD0K7SctpIbHVBQfFxMGUemeSK/cVSFRm3nYM2TM7XU1ncmMchAM6M7GMru2s9psmhfXmxF9bN801d6g07vD2h8JBcolnN5u/VlOK2cz0n5SuYG+VktBMLNUs9J2bY3EcE8z8GMsUzWmcmM91FWKij5rPobthxjgyia+SNFuZFKzwPg5CAqA791NisrumYhLzlZlgjJTu6ubaQl8pJFSm+1wlYS20u6jMdO9y6WYtJJ6ldLPu+qe+LgktJJt1kwSrcK12c/88q7/l6peUYh9qnWaOfU8p9Caxvqdx/V6omTqlh5ulq3X6Qdi3ZdgPl8hXqD0qNU1qj9PIFZ5eKnNi6uvm4FSnpLnd3Qti/eJS70Odchhrj+3+udYx24vj+mN005z7N8I1NTMx69Nw7aRaZ9e6JtU081DjzUytMwl9p3bs4bWLJM1N6qv7TCVp7hJUU7K8qXWqdbQ7910tL4XrzPY8awzXGAa5YHaX45xaqD1KqffFdFdTkzR3zw99UyR4rl+222MKvatfUoprqnVcXZNqHXfsTtUuvd/j2cuO+CYdAAAAAAAAcMl4SAcAAAAAAABcMh7SAQAAAAAAAJeMh3QAAAAAAADAJSM44jroXfdcmdcsMipJ3SzQ2tICpmHhVrcIY0vhE3YRxrSwdFp8e/gYNvShuLi4O7SjLJB8wYspp33e125UFqNPC8K69klYkNn2De8ujTG1Y4S+hYWe037MzCWVFkiemUWWKwER6/7bY09TqExBuNq1NCEoCsdoZfqmPUvth+ZznYfeU7Nwffqs4/li+vtFoaWpXUy5sHh+Ubqm9sZtz33++7TPeTndj4wYCuQ+7sr9LyzenO+V5lqLQQ7b7T3umwlhCAsvV45bZOqXHs6tMPXZ/vlzMu3uWEpST6Fc5rM6KHzW1cCpoeNKtevyqoRB2NSxi92Jfc6p7l7g7hmSv8e4e1EaV/L3ujzG9vbcvXY9rufap2F1f1cDHYQqY4z6xQrBNqm2mptTYxbOl0qg1km4hl2tOQ197Wcda25XFw2vz+MYhRo/ueiAqwvP0xqhTqkFVQxvT+ET/t5V/Jm+MoZpj7VHCLBwIREx4Crd99/zwv2eJHyTDgAAAAAAALhkPKQDAAAAAAAALhkP6QAAAAAAAIBLxkM6AAAAAAAA4JLxkA4AAAAAAAC4ZKS7PihCYqTvW3x2Oyn0n7i0xz0mAJqxy4GDpn8c44LDDCsuOmjxorl00CqbRlp5fTgBVqF96QIcQxTTyrSvQsJa2t4k5rDuZhmS0Nw+L0MaUnrfpf0w77tyXqz2+P+tlikF8rqwSb4XvA9XZI7r4Z5mW9P9z7W7+2d1jGTX+3CsG4qJtHZsd26NcD2NcXxSe6XWKexHOrcu3FW4tuPcftGxjPuxz3tG5V7n7qHLeB4OP/bpnj8xc0NMAQ21x75SX10ya2pfpVrHtLtacD1uqB13vAB3rXWlcWruq+xKv72wb3afi+/Dj3GxP6fHWsdJzyBS6qudG1Jf9wDg4u4v1/ynBgAAAAAAAODq4yEdAAAAAAAAcMl4SAcAAAAAAABcMh7SAQAAAAAAAJeMh3QAAAAAAADAJSPd9RpoB1O1do+PMiUdmQSVllJcZjM/9NSM7dokNZfCEvr20G77h3SXPh2esGb7KqSpjZGscx8mxDop6Kabg5EycVJKlGtfrHxCqes7775vaneJZ/Pur60TM8YsbS8kgs1NOus0HCXbHpLNEpemNt3jCedSzOYhRWlu3stROLmOw4k0N2l4J+H/RbnP1X2mqe+6fbt/PLfM+bkK6X2xvZBe69rdNXm/Kr2VEebrbtM6a9uTucfEJFh3T0spZtOQ2rwy88PMXzx2L9Lk7vZ5mVJc93TOhZrG1iOSPXYtHbeZud5T39Re+PzsORDqkco5Z8/Z1HeEeuQaTS92rqzMtSmVszLnp+1V7juV+iX2bdv37Fl4H8cpJdHWNWHOsEm+vus0zJ+LnvZjmJRCn9JdXf0yDzs9t23ps071y/bnV6lfUrrvrjWN5Gv0Up1ie+b65QLDNscxws+Elfk61RjpZ153f099x/iZvq3MvXk5fIwWrpGu8MzCXYHpXrk6++eu1hfSwr90DHyTDgAAAAAAALhkPKQDAAAAAAAALhkP6QAAAAAAAIBLxkM6AAAAAAAA4JIRHHENtMNDtXZ4j05hUUS7mHJxcerJdns7CH0PzUKOlcWbJfWD7f1zbZJfZLIfhEU0w8KRbs3U2Ne0pwU+4xh7WtR5FG4NzcpiyqHvMiz8uTALiqaFdI9X2+fW8cqv6On6StLtdmOrbRIWCp6a0Ie0DnL83yGm/9Is0rxu3+48a37F3MOw9O7ELOo8HWG1bxdIIfnDkRZkdospH8WFs/0BvW0WSL6z2v5M132358yj7s+Lo3C+uP6p7/Fqe99cmyQtCtfD0rSt211whO16dRZk3jF4Jxy2UvsqLJA8mWwfjNQ33UP7gTmgM3//a+bg9/CBlG4PaRFp0x6DFZZmjnJt0v5OohS4UQl4SAtZH5jr0rVJuU5x7aEu6uYcqNQ0UjgXw3nv+o5x7ZROxAuuXcaY+9ycum7fPhhurpby3F65P7h7zI2Jv+/MVv66nLXt2miagqhcc/ErHkszyCrUDTMT+mDrLUlmWl733/EEq9Q0knTiQrLCvh2Zc+A4hDPcDjXJbVPX3Omh1llt1zru9VKqo/15mOpoV5+l68HV/pWQCSlcw1ehdlHt57n8M6GpBeLPoKZt6g9GHMP8jNzD52fvU6FvC8E0NpwhPNKwmTIxRCrUJO4emoKvXPjWqWu99RXBEQAAAAAAAMB1xkM6AAAAAAAA4JLxkA4AAAAAAAC4ZDykAwAAAAAAAC4ZD+kAAAAAAACAS0a66zXQbt5Qm9y8R6cQJ+NSUULiSzNJsJJ86llKMXPtMR1teBLa6tD3Xc1McpBpk6RVSH117asQIOeTdXzfWiTfBSukKKXUNJfElJIo50t/kE6m2+fGSUgrOzbtKZXqjkm7kqSJSRCbxjyvbctwjazC/w9x/Q9NspkkHZs0tllIPEvtUxN5NRkhBisdoaU5mVOS79wco5TiehRSolzaakoxc6mvlb6SdGdpEtbCueVT08J5H66Huem/TElo5lpzbdIZQWgjJP9ae0psLSdUuuTukNi6mpn2lGKWLggzWU5cspmklSnPWryPh3u2a68kmi58bJlLmW1XJN01pcXbY5ESW90xCuPaFNfQnmsac62GdNdYv9haJyXWm/t4oaZJ7XtLgh1DqmlC98r86eZgN1dLeW4/nmy3p/TMOyZxcBauP3fPT1bhA3R1iru3S9K8uZx26aapX45S/WIqimlIsk/vb9e6plLTSNLcXECV+uVE/ryI6fS2fvG1hxsj1SmuPae7+vYT0+7aJF/7p58TcrqrayvErV604txn65SY3G2GDT/btnDLbu74p8vJ7Fz89lelfpn72sM9h2iL8EZS/WLmzxbqMLnE69Pprqsu3fEvHQPfpAMAAAAAAAAuGQ/pAAAAAAAAgEvGQzoAAAAAAADgkvGQDgAAAAAAALhkBEdcA+3mTbXprVMNhYUxXUhEWPy+hwWnbf/Qt5vFlyuLKaf+cTHlw+GLKcd2s+hmXGR5xwU+U/+4ILPb5eq6qG69zLhIqGlKi7yutndksQyhAeF8OV5uT1EHYbFh1+6CIKrSYspzs/iva1u3+w/bBR3cDAsvz2xwRFosOgVKhAVW7RjDF15Oiym7Y5cXXjaLDReOW2o/iotvu8WUU3CEX2T5seV2YM/tRVjoebk9xtEyvI/Q7hZZTqErbjFzd01KUg/tI2SKeJU1nVNfG9KTFsoP4Qymfwt93eLL6XJKc4YfOCzqPDXzWXh/Su0uqCktsuwWUw7nYTN9z0gQ8u2VQAl3jELAVax/TJ1SqmliGFalThk+RqWmWY9hrvcUhmV2w10L0lm1jqmL0mnvAmEuPDgiNIe5z82VKaTHzssTfw8+mvhr6sD0T2EJFatS0EEIuzDt6R6c6hfXfhgmUFfrxICIQkhWkmoSJx2jpQ2+CgFQptaJNU2oX1xwRKV+STXNu0398m4TkCVJj6f6ZWFCskwtL0lz8/ND+jkh1y+VoIOBbWcpzGeVn+fSz4Qu37HyM2i6D6Qfj5p9MyPUNIVQyhxUaeaMZXgj4dpprr+raaRcv7xnrDP/8874Jh0AAAAAAABwyXhIBwAAAAAAAFwyHtIBAAAAAAAAl4yHdAAAAAAAAMAl4yEdAAAAAAAAcMlId70G+kO31A9und2pkoQWkllScp5PTUtjmL4xHS0k/JgxcmLrGOmu223LQl+XuiYpPiK3yT8xCW07eWaUQMZCEloPsUZLk9C0CGl6KaHy2Ly/SSHlaxX2LSe2msQzk0gs+ZS2o6mfUlPq1o3JduLZzXZi+7oktJjuGmKbJib1NfWtWBbSLFMCrktYS31T4pn7/FJq2h2TWHYctnfHJLNK0m0zRko8u70wCWumTcpJaMfmOklJaO76S0mGpXTXwgQTE88KyY6V1LSYjpYqHXPq+2SzmpjG7e6VKWDUpWeG+2oLqaGThXmDrk0+sbUvQ9LtGOmuuyqmu7pjn1Jxbd9Qp6xS/eI+vziGSWYNKa6V+iXWKabd1S5ntheS7EupzdU5w3Ye2KbanOjmVElamPPleOKvyenEH9DJjkmuKcU11UCurkn31aOpSaFfpRR6X5PcNLWOS3GVfDJr6pvsWtekmmYVCnfXv5aWG2qB9JmYWsXVNOsxtvvmxNbt/XC1iyQdhTrlcZfuugjprq6mcWmtytef+xlkn4n1doh0/Zr7USyX0/zp7g8xNd0VUeE+HidVU6eE92c3F+7Nk5Qye2CS5cPP/82cLy3UKTGx1fVP9YtrP3Xse6inxsI36QAAAAAAAIBLxkM6AAAAAAAA4JLxkA4AAAAAAAC4ZDykAwAAAAAAAC4ZD+kAAAAAAACAS0a66zWweuSGVrOb53uxS55JoTghmdOlrKV0F5dSk/q6xLM0RuprU9MKyayST0gLIZB+eyFxMG3PJf/ktEDfXmKT0MJJ4NJdQ7jNyqW7LvzBmIfkoF0TzxYhXSsloS3MAXXJWJJ0a7qdVnZn5VOwXLKZ5JPQXOJr6utS0NZ9hyeh7XqMz+KS5Zbh/w25JLRKCq/k0+nS5+faU5JaSmx1Y9wxyWaSdGTGSCmuRyEJ7cRcPykZ2SWhpaQ/d11LCtFdvusY7NwX5k8378S5KCVsp0QvyyWhh2FjIu3wxNbVwiRbh3S7ZvpKUl+aAxITW7fbXdt650zfdChLx7ggpbjGwLpCneLa0+cU6iJbp6Q02VJNY5tDYmtIXjdTVJj6ch3m0l3TtVpIgh3FjjWNFO5d4fqbmyTXNkJNk+Zrl345P/AHP91jHp9u1yqp9nBJoDcmvsao1DqTkMDqa52QWL+n+iXeKwOXRO8S69d9tz+rVHuMUeu4+uWk0Dem0Ic6xaW7noQ6ZW5qmvRzgvuZQgr3/QuuaSo/o6W+aW63p35KSjX323Qeppu2rVPS+zP3v37gx10tUurr9hip9nC1TqVOkaRm6qLz1i+ruU+zHgvfpAMAAAAAAAAuGQ/pAAAAAAAAgEvGQzoAAAAAAADgkvGQDgAAAAAAALhkBEdcA4uHZ5oc+oU936OwBmpPCzKnR7qVBZntAsLF4AgzRlyQ2fWtLshs2uMYZu3XFDIRF1l2xzkde7cbtfVua8wamj0s5mrWNtYyLFp+0sLBcOOGBX1duwuCkPzCy5J0NNn+AA+nfoHk2xMX+uD73ghjHJgFkt2iyZJ0MNle7LS6mPLULNQ82WMSgAvoWIbPxH1+adHk1L4wKS3zsD23cHJaZPskpL+4MIgU5ODGPg59XUDEun17jLTI8tK092WYJ9Miy+b0auH6K63fHRdZ3h4kL+5v2tKcGvatm/eySteDW2+6sG+Sv5eEKUMT81mtUuhDWLu4mbk5LshsmlvYnjtEbsHq1HcUcV3wtEi26ZvqhspnXag9YkhIoU5J57gNrSrUNC68Yt3ux6iEQfjF08OJUfjqQF7s2wVchUXSw9zn5sr/v707j5OtKA8+/nvu5bIjigoiooCi4g7uO2rivkQFBXGL0bjExCUi8X1NxF3cjXGJKy64gytxCa9BERIRAxFEEIIsioq4K+u993n/qNNOM7dqps9M9/TM3N/38+lPd5+uqq7uPn265pk69Wyg/ua3drk+qokqGuOU2vZWIoCtKwmuoJ4goDVO2bJykFrXSPqwRWVcVMqPPtapjV9qYxeY3PillVysz/illSSrmvijNaZpbK+OJxr7QG38clWPsU5rXNRKHFFLElFLEAH1sU4rQcvGxveh+jfIBBNc9fq7q08iwNb4pfY3aOPpah1pJlfp8bd+e0xTeb6e45Q1PZJB1NronTiiOq7tscMMFV1/1VWj11sAZ9JJkiRJkiRJU2aQTpIkSZIkSZoyg3SSJEmSJEnSlBmkkyRJkiRJkqbMIJ0kSZIkSZI0ZWZ3XQXWb7cFsdUCP8pqFrPRy0I9m1q/LF/9MuQtNrtrMztaM7NOJWtanyw8rexoPV5fKxNatY2+mcZqTbdSB1UyoWWj7MZaJslY/P8FNjSyS26o9K2VxfXKNfUPZd3aTVMHbdnKVlYpu0UjA1krQ2ytfCsT07pKP1qZzVptrGn0b1I2Ng8mm6pl4m191rUsrlDPyFbLNFaeb9Ptrf3l6sbz1TK5Xt1oo5qZtZHFrJUhtpb1bMP6Via0yne18XzN73s1zWWjbE1jP8zGQaqabKzP78CI3ZpL9MomW399rd+j2qGklcWzlrGslvG1lK1urrbRPATUMrb2yI420Wx6fbTGKT3GHn3GRX3GL81xQ6WN/uOUWtnRM7a2s8aOvr3VRvU73PqcGk30Sh9dK9r4LWllAKwdKzc2vjytY3C13VZ2+srvxvq1jd+dLTYtu25DvW9Xrm2MdSoHo1oGeahnod+yMv6B9piklvV1beMzbY2jqs+3TMY0teyurQyxtfFErT7Ux0VQz87aKlsbT7Qy1tfKNsdFPcYvrYyttez0zeyufbLT98hY3/u3q9fvQ6Xx1nG5laS0lhW+0bXqt6H5G9U49vXINl7Lttr4k6mdhb7yWfUZezTHPz0yzvcZFw1bf+Vkw2jOpJMkSZIkSZKmzCCdJEmSJEmSNGUG6SRJkiRJkqQpM0gnSZIkSZIkTZmJI1aBq7ddA1sPxVv7Jg6YpbW+bjtxRGXbWBZkHn17a8HietnGwr2tNnokn6gnfaiXbSaU6NHGYj/rlmjsBNUkEY3FQGsJRTY2V4EffUHYjVvUV/isLbx8dWMx5eYCyZWFmq9olK0tetxqd21zQeZNt0drMeVKG60EES2tRZ2rZXssyNwnQURrMeXqwsutRbZbyUMq/Wgteryhsr2dlGT0hZNr+2GrH60FktsLJ9eer/HdqSxmno2Fl2ltry7AXi/aS4/F/Tc2kjOsaf5QVdro0Y9Wbpva9tb7Ga3F9jeMnsghKm23E0Q0tlfeo9aCzNXFlFvvcbVsvWhTrwQkPZrtMX7pVba1+HavRBWNsj1+81vjl17jonEk1OrVxqYf9kTHNH2OW63F6KubG8frWv1mgojRt6/fWP9i1xIOrWmMMdatHX1M0krkUCvbHOs0Dka1sUo7wVWlbM8fnqUev9TbbSW+qvw290wcURuTNMc6lbZbY5ravtVqtzXW6ZfgqjKmaX0n+4xfGt+z5m/aYvX4HaB1nGw0XWuilRejPk5plG39llTKR2XsAvWxR+vz6zN+qSWTKJ2rlW30bZLjl04tic84OZNOkiRJkiRJmjKDdJIkSZIkSdKUGaSTJEmSJEmSpswgnSRJkiRJkjRlBukkSZIkSZKkKTO76ypw9bYB2ywsY02vRDdjyJpWzXTTLDt6FrM+mXX6ZIJttjGOzLPN7aNnQqsmMeu7K1TT846eAbCZkqiWXbLR7sZq2ljItZXMT60MnJXnW9PIbLZ+TSOrWGV7K2taLRNaKzPr2sbz1cq3MsH2yeTa2gX6ZoNdrNZnVVPdtVpZf3tkfd3Qyn7aI+NZM8tsLftwMxNapWwji2vz9fXIYlbNhNbM4trKxlU7NtSbqDfQc3vt7WgdiqpZX0fPKg6Nt65xXF5Ty5DWN9tqn4ytPY617exmlY2tz7pPRswe+8CkDjm9f+dq5ft8Hfp9derjlx7Jzftkgm2Wb2WWr7XReh09nq/+nWz0o2921z6fdzUDYL/svLWUre1duZLNsvGb3z62V9roMX5pjVOuao11eow96uOU0dtttdEsW91at1zGNH16UR+n9BsL1MY1zbKV7b3GKa2MxK0s9JXnyx7jlGy02xy/1HbbPr9dPX9MehWvvpRW5xrvc61kK2Nrn4z1rWN7Ldtqn4ytjZfXzu5aqdCjjT5ZXOfcXnu+ecpeXXvDx8iZdJIkSZIkSdKUGaSTJEmSJEmSpswgnSRJkiRJkjRlBukkSZIkSZKkKTNIJ0mSJEmSJE2Z2V1XgQ3bBjHF7K6LbaOVxazX840ha1qvDLGtNqplG1m+emSIbYbTq5lu+2W7qhVvJXetZUhrPluPbmTjzchK1p9ovp+bbt/QyBLVaqOaKazxftYywbYS/bSyvta292mj1W7LZPMQbapP71qZyfqUre23zWx6PbKjNZMdVzKk9epbI+NZM2NrbXsru3K1bI/saNDzA+xRtE8Gxx7/SszW96zxumtZz1oZyGqHkl6ZWaGRdXLxZXtlN2sd2/tkdx21/hRMbEzTM+voYsce4xgX9RunNNptbq984H2y3o9jPDmO41Pr+9fj17L6+1DL7gtsbGRipvJ+trJqVjM4to59zYyto7dRH+u0xjTVzb3GKr3KjlxyPMYxpqltb44xxtBGdT9qjYt6ZGZtZl6vja36jHVaZZvZ6Tfd1PrNr5bt+9tVq9D4vo/YhbkfqY6LWmPESvWeWej7jT0qG3uPi0bPej+WMckYxzUbzO4qSZIkSZIkrW4G6SRJkiRJkqQpM0gnSZIkSZIkTZlBOkmSJEmSJGnKTByxCqzfGthmgZV7rHm45EkmerTda0HmvgtA90jOUC/bo136LVrda4HrPnosshyNJ6wuLN1a+LW1SGhlUc5mIo7q4qr9FjfeWE0c0epan7L17bX9qNfXrO+Kt0u9uvsiv/DNRZPH0IfqIss9F2+ulW8vslw7GNWLNhdOri7SO4akAb3aGMMBprUf1r7vradbW6vfKNs6ZlQXWW4s1t4rwU6jH0uYyGEcbYzlcDGpQ84YdsPlkDir2cZYximjP1/rd7VPUq5m2drGHr+JvVVfYL+kMlk9uPdY/L45Tmm0USnfHquOnlCrmh2HMYxfxpC0qvf4pdrw8h3T9EqG1fP5Fj1+6bGIf+t3bnLjlH6JKqrf4XEkw2pZ7N+VPZMC1f/mbR3PKmUbneuVOGJCybD6mti4qGWesusnfPhxJp0kSZIkSZI0ZQbpJEmSJEmSpCkzSCdJkiRJkiRNmUE6SZIkSZIkacoM0kmSJEmSJElTZnbXVWDDlhBbzVNoHBk/G8aR3axP2XomtNFTrPTN7rroTD79EhVVQ+d9snyNRfMJK8/XyOQTtTb6Zsir/RuhlSWq+l40MmP1+vx6ZJNtGcPn18wQu1iTyqY3Bs2sYuPoxxiyRPXK2DqOTKm9soP2ySbber5N2+izu7ReXvPoUjtktDJG1vrR99+OtWNtKwtkn/1BW4qpAAAgAElEQVSlV8bW0bPJNvXsx6gmmjixT9sTOvZNKtv8nNtH7Uff351JZYXvM36Z1BhqrvK1Jvpk3W68z9Ws9c2+9cgE23odtcz348iC3dqPqmVbbTS293i+atEJ/l2y6IPXhMY00HNc0zPT6chl+7Q7yXFYLcN639fc67d5Qp9r63u2ZvQxVC3bPFDPZN9qozKmicYO19wP+2R3rdbvOaYZY7bVkZ5vsYba3bB+Qs/RcSadJEmSJEmSNGUG6SRJkiRJkqQpM0gnSZIkSZIkTZlBOkmSJEmSJGnKTByxCmzcKtmw1QJXSJzYYvSLb2KpE1I0+zGh5+v3+loLCC+yXagvEtos2mMx+lqfW2spt1YQri2kOpb3s0/ZSa5uXHu+Hgs99257HI30MMmF56vPN6kV5idTtrlAcp+2x7IY7xiSF4zDmk2fsL24cY8FmVvP1+P9HD0tTb82+qwiPpbPY8m/kxNqd6mPZct5TNMo3+utn1QCDOg1FpiU9rGhNfbokwBh9CQTY3k/e7XbM4HFYvVJ4jaW5xtHIz1MbDH6JX4hExrTQM+kVX2ebwzJIHolgOrxdO0eVMY0rdK1BBH1Jnr1uTXE6LPH9Ut20u8NXfS4ZooJrjaun+yAypl0kiRJkiRJ0pQZpJMkSZIkSZKmzCCdJEmSJEmSNGUG6SRJkiRJkqQpM0gnSZIkSZIkTZnZXVeBjVtArJt2L0a0WrLJtowj/d6ksn/1MYZsQLXO9e5vj/czar0by2c6hjZ6Pd9SP+EqMrHMa8uk3UV+4Xsfnib0ulsvo9q/HhkjeyYVW/wBdImzio3l6ZY6u+tqMcnD8qSSUi/1mGa5jMPG0eU+45c+36k+Y5pJjgUmNhZ3/LIgKzHr9ljaWGRm1t7PN4Y2as32SaLcMw19r2Nij8Ktt2Kps632Kr7Mxi8bJxxFcyadJEmSJEmSNGUG6SRJkiRJkqQpM0gnSZIkSZIkTZlBOkmSJEmSJGnKDNJJkiRJkiRJU7bZZ3eNiNsBLwYeAFwHuAg4BnhtZv5mjnp7AP8APBTYBfgZcCzw6sz8yWR7fU25RZJbLLOUJ0ttBSaUmlhm1kla5G7WzBrUzJo2/TdpGXRBy9RYso1NynLuW8Oiv2t9ssbO+cBozLaqiVnmvzur5nexTybYvuOX6vMt7Ru3aj4njZVjl6XR6/vXd/yy6MJ1i25hFX1+85l07GXVzKSLiC0i4q8i4qKIyIg4coQ6BwOnAIcA2wBnAXtSgnanRcTujXr3A74HPBO4LvB9SqDu2cD3IuIOY3hJkiRJkiRJ2kysiiBdRBwAnAG8D7jRiHXuBHwYWAe8CNg5M+8A7AP8FLgJcHRErJlV78bA54AdgLcA18/M/bryZwA7AV+IiG3H8NIkSZIkSZK0GVgVQTrgtcAtgI8Bx41Y562U033fnZlvysyrADLzHOAvujJ3Bh41q95rgGsBXwZelJmXdfV+3pX9A7A7ZZadJEmSJEmSNK/VEqR7HbBvZh4CzLseXETcGrhnd/dNsx/PzJOBk7u7hwzV2wl43KBeZm6cVe884N9m15MkSZIkSZLmsioSR2Tm+3tWeUR3/d3MPLdR5ovAXYB7DG17IOX02EuAb8xR73HAHSNiq8y8smffess15SKtWMt4odFlvcCupHkt68XTl3PfJM1vGY8RHL9IK9uyHr/UrLT+LsKkYy+rIki3APt116fMUeb07nqXiFiXmVcP1Ts1M9fPUw/ghsCPFt7Npmt8BTb88Y8TeAptNjajA6qkJeQfiJImxbGLpElx/KJ5VOIvY/1V2lyDdHt01z+fo8wvu+s1wA2Ai3rWg5LEoleQLiJ+P0qx4TsXv+4NfZ5CkiRJkiRJi3dexNzzlzNzh1Eb21xPkrxWdz1qsO2Gi6wnSZIkSZIkNW2uM+kGp6qunaPM8HszmM+40HojGyXCGhG3BH7Qt21JkiRJkiSNzZ0y86xxNba5Bul+DNwa2GWOMtcduj3IGPvj7rpvvXE7B9hn6P5gbb09J/R80jQNThl3/9Zq5T6uzZH7vVYz92+tdu7j2hwN9vu9gJ2Gtp8zzifZnIN0ADvPUeZ63fUVwG8WUA/g4v5dm19mbgD+FKkdnP+cmb+YxPNJ0+T+rdXOfVybI/d7rWbu31rt3Me1ORra7y8BLpnU82yua9Jd1F3fZY4y+3bXF2bmYBHAQb39IqIV4BzUuwLwoCVJkiRJkqR5ba5Bum9217ePiL0bZR7RXX9taNsJ3fV1gfvNU++4zNy48C5KkiRJkiRpc7E5B+nO726/cPaDEbE/cNvu7pGD7Zl5HjOBuhdExJpZ9fYB/ry7+8Gx9VaSJEmSJEmr2mYZpOvWdDusu/usiHheRKyDP2VO/VT32Mcz87uzqh8GbAQeAhwREdt09XYFPgdsDZwEfHayr0KSJEmSJEmrxWYZpAPIzE8Br+/uvhX4eUT8N3AGcH3gROCvK/X+E/g7IIEXAb+IiFMomT5uTknocMDQOnaSJEmSJEnSnDbbIB1AZh4GPJpy+usWwN7A/wDPB+6XmX9o1HsHZU26LwNXAbcCzgVeBtwpM386+d5LkiRJkiRptQgnfK18EfF7gMzcYdp9kcbN/Vurnfu4Nkfu91rN3L+12rmPa3O0VPu9QTpJkiRJkiRpyjbr010lSZIkSZKk5cAgnSRJkiRJkjRlBukkSZIkSZKkKTNIJ0mSJEmSJE2ZQTpJkiRJkiRpygzSSZIkSZIkSVNmkE6SJEmSJEmaMoN0kiRJkiRJ0pQZpJMkSZIkSZKmzCCdJEmSJEmSNGUG6SRJkiRJkqQpM0i3gkXEvSPi8xFxSURcHhFnRMT/iYitp903aS4RsV1EPCciPtvtt5dFxLnd/vzARp0tI+KJEfHpiDgtIv4QERdExL9FxMEREUv9OqQ+IuL6EXFQRLyv299/HxEPHncdaTmLiNtFxEcj4qcRcUVEnBMRR0TEteept0dEvDsiLoyIK7vj/zsjYrel6rvUstAxSkSsjYhHRsRREfGdiPhtRFwcEcdFxHMjYu1SvxZpFBGxY0Q8KiL+JSLO7PbdZ06qnrQYUTwoIj4QEf8VEb+LiF9HxAkRcfA8dZd83BKZ2fc1ahmIiEOBI4AAfg78ArhN9/CpwP0y87dT6p7UFBFbAhcCu3SbfgFcDNwS2Krb9obMfPGset8F9uvu/hY4D7g5sF237VPAQelBTctMRGwHvAJ4LrBlt/knlGP3czLz2+OoIy2liNgCeApwOHAj4EOZ+dR56hwMfAhYRzmOn08Zu6wFLgDunZkXVerdD/g8sANwGXA2cCvKb8avgAdk5mljeFkS0H//XugYJSKOBh7T3b0MOAvYCxj88XciZf++cnGvSLqmhRzDh+odBvwDsH23+RLKWP7lmfm5cdaTFisidgeOBO7fbfoR8BtgH2AwuekjwNMyc/2sulMZtziTbgWKiEcBrwcSeAKwW2beFrgb8EdgX+B90+uh1JaZVwH/RdlHbw7skpl3AHYCPtAVOzQi7jur6onA54DbAztl5n6UQeyruscfBzxpwt2XeomIG1L23RcClwOvoRyjd8/MOzYCdL3rSEspIg4AzqAcx280Yp07AR+mDHRfBOzcHfv3AX4K3AQ4OiLWzKp3Y8qxfwfgLcD1u+P/Tbo+7AR8ISK2HcNLkxa0f7PwMcqJwAnAPYAdM/OOlH362d3j96QENqSxWeA+TkTsAHyVsl9vCfwzZR/dNTP3nSNAt6B60pj8FrgzcDxwm8zcqztG7wA8D1hPOT5fY0bnVMctmellBV26neRHlADdYZXHD+geS2C/affXi5faBdimsX0wyy6Bt49YJ4D/7Op8cdqvzYuXwQXYFvhOt2+eAew9iTpevCz1BTin20ePAv69u33kPHW+1ZV7V+WxuwyNXR4967GPdtv/DVgz67G9gN93j79g2u+Ll9VxWeD+vaAxCrAN3ZlNlcc+0dU7fdrviZfVdVngPr4W+EJX9iLgLiM+14LqefEyzgvlHyFrG4+9o9s/Lxk+Hk9z3OJMupXngcAewBXAOyuPH0M5HQrKLDtp2cnMyxvbr6IEKABuOmKdBE6q1ZGm7HDgTpRB6T0y85wJ1ZGW2uuAfTPzEMpp2HOKiFtTZk0AvGn245l5MnByd/eQoXo7UWYgAbwpMzfOqnceZRB8jXrSIvXav2HhY5TMvLwrU3Niq560SL33ccoso0cAfwDu2R23J1lPGpvMPCkzNzQe/lR3fX26maXTHrcYpFt5HtFdH5uZv5/9YLcjfKm7e48l65U0Pnt31xdMuI40Md0pq3/b3f2bzPzdJOpI05CZ789+a8ANxi7fzcxzG2W+2F0Pj10eSDmD4BLgG/PUu2NEbNUoI41sAfv3fBY6RnFso4nou493p+X9Y3f3JZl54STrSUvsl0O3d+6upzpu2WLO7mo5GixKe8ocZU7vrs14phWly1p52+7ukSPW2Rf4s+7uByfQLWkhDqYsRnsG8KWI+HPK4uA3ofzX+jTgqMz8zSLrSCtBn7HLLhGxLjOvHqp3as5azLlSD+CGlCVBpGVhoWOUbk2jwWwMxzaatocAN6Cs7fXuiLgbcBBwM0qA4zTgk5l58ZjqSUvp1kO3z+6upzpuMUi38uzRXf98jjKDaPANI2LN7GmW0nISEeuAPSmnZx/abT4851gcPyLWArsDjwb+L2VNlw8An5lsb6WRPay7Ppkyu/mhlTIviYinZuZxi6gjrQR7dNejjF3WUP6ou6hnPSinqRik01QtdIzSLT6+K/BgyuyjXYCvUBbZl6ZpMD45EXgP8JezHn8y8E8R8ZzM/PgY6klLabCfnp2Zf+hu79FdT2XcYpBu5blWdz3KB78F5dzqucpKUxER1wZ+PWvzhcBzMvPYOer9BthxaNOvKT/yRxmQ1jJyp+76IErW7WcwM7393sCbKX/EfTYi9snMHy+wjrQS9Bm7QPnP8kULrCdNzULHKBFxGiUz7MCVwPMpC5ZfNYm+Sj0Mxif7Uxa8fxHwSeAyygL6bwBuA3w0Is7OzP9eZD1pSXQZXAeZt98w9NBUxy0G6Vae9cBWlEw5LcOf6x8n2x1pwdZTztVfR/nP8U2AGwOv7maAfrFR71uUAfAulP9WXAf4J2BdRHxwjgWYpSXRrcGyQ3f3CuBu3UKxA5+JiG8BP+zKvTYintm3DjODCmm5G5zy0XfsstB60rQsdIxyCvA7yj/X96AsfXAosHVEvHGOBc+lpXCD7noL4F6Z+Z2hx74SEd8EzqSM5d9G+cfiYupJExcR2wHv7e6eBXxo6OGpjltMHLHyDGZO7DJHmet2178bmrIpLSuZ+YfM3D8z75mZe1ECdK+n/Cf5CxHxjEa9h2fmvTPz5pQf/xdTTpd9P/DKJeq+NJfhWRSvnRVsAyAzf0ZJ+Q7lP8wLqSOtFH3GLjCTbXCh9aSpWOgYJTOfnpn3ycx9KIG6ZwDXpmTh/FCrnrREBmOU984KtAGQmZdR/nkIcPehxfAXWk+aqG5ZgqOAO1Bmdj5h1hpyUx23GKRbeQYf/M5zlLled+0inFoxMvMnmXkYM4PYN0fEdeap86vMfAPw192ml0TELSbZT2kElwCD05N+MEe5QWa1G3Xle9XpThmXVoI+Y5crgEFylD71wHGPlpGFjlG6f2K+D3hkt+mQiHjApPopjWCQlXWU8cla4JaLrCdNTEQE8EbgUd2mJ2fmqbOKTXXcYpBu5bmou77LHGX27a5Nc62V6CPd9fbAHUes80ngasoxzanymqrutKT/7e7eeI6ig//YbaBMd+9b54qF9lFaYr3GLkOnBA7q7RcRrSVaBvWuAH6x8C5KE7PQMcp/MPMH3P3G3Smph3O661HGJwCDM7kWWk+aiC5A90rKmp8Af5+ZR1eKTnXcYpBu5flmd/3giNhh9oPdzvCQ7u7XlqxX0vj8duj2XOfzD7ucmVlIo9aRJunM7vq+c5QZpHw/OzOvWGAdaSUYjF1uHxF7N8o8orseHruc0F1fl3aQYlDvOJMHaZla0Bil+6NvMCZybKNp+tP4pAty1AzGJ39kJlvlQutJk3I4Jes2wKGZ+eZGuamOWwzSrTyfA35PSR7xnMrjT6QsUnsl8Ikl7Jc0Lk8Zun3yiHUeC2zX3f72eLsjLciR3fWBEbHJ6RtdcolndXc/vYg60krwTeD87vYLZz8YEfsDt+3uHjnY3q3NOBjwviAi1syqtw/w593dD46tt9J4LWiMEhF3BfbpW0+agMG6iHcFNjn1ujs2D2YmfWYo8LDQetLYRcQ/UhL5QAnQvXGO4lMdtxikW2Ey89fAq7q7r4mIAwcffkTcC3hX99ibM9MFlLXsRMS9I+IjEXHXiNhyaPuWEfECZtakO6Lb34mIPSLi8xFx/4jYZqjOFhHxRGYy83wS+J8leinSXI4FvkL5nf16RNx98F/kiLge5R8uu1IGAEcsoo607HWngB/W3X1WRDwvItYBdAHpT3WPfTwzvzur+mHARspZAkcMfgMiYlfKd2Jr4CTgs5N9FVLdQscoEbF9RHw5Ih45fHZMRKyJiIcAx3SbTqT8PkhTkZmnA//a3T0mIh429PfndpRg3L6UDMWHLbaeNE5RHA68AkjgOfME6KY+bol2JnAtV93B7SjgoG7Tz4BLgdt0948GDpqVoURaFiLiLcz812w9M+tV7M1MSuqPAE8b7MMR8Tzgrd1jSVm76zLgFpRZpQBfBR6Xmb+b6AuQRtQlPvkScI9u0y8ox+tbUwJxPwAOzMzvL6aONE0RcSRlBvSHMvOp85Q9gpLtEuDXlIDz7Sin8p0IPLiWlT4i/gZ4OxCUU6LOoox5tupu3z8zf7r4VyNd0yj790LHKBHxKMofbAMXAL/q6m3bbTsZeGxm/hhpAkY9hkfE1sDHgEd3m35LOT311sA6SpbKgzPzhHHUk8alm0H3iu7ur4HvzVH8K5n5uqG6Uxm3GKRbobrZFU+jnPp0K8paF9+n/Lfuw+kHq2WqWzfxMcABwM2AvSiL4J9LOWi9PTNPmVUngAdRTufeG7gp5Yf9XOCHwPsz87ileg3SqLr9/WnAkyjT4remDE6PAV6TmX8cRx1pWvoE6bryfwG8gDJ7IijH8A8D78zMq+eod1/Kf6fvRvlOnEf5T/ab/E5oUkYM0i1ojNLVuzvleL9PV28HSpDvHMrMu087ptck9fxHSwAHAk+nJHfbgZLN8ljg5Zl56TjrSeMQEccz93rPwzb5Hkxj3GKQTpIkSZIkSZoy16STJEmSJEmSpswgnSRJkiRJkjRlBukkSZIkSZKkKTNIJ0mSJEmSJE2ZQTpJkiRJkiRpygzSSZIkSZIkSVNmkE6SJEmSJEmaMoN0kiRJkiRJ0pQZpJMkSZIkSZKmzCCdJEmSJEmSNGUG6SRJkiRJkqQpM0gnSZIkSZIkTZlBOkmSJEmSJGnKDNJJkiRJkiRJU2aQTpIkSZIkSZoyg3SSJEmSJEnSlG0x7Q5IkiRJy0VE7AE8tbt7WmZ+bmqdWYYi4vnAtQEy8/Dp9kaSpNUlMnPafZAkSVrRIuJ44L7d3T0z8/zp9UaLERH7A//R3f1QZj51er1ZfiLifOAmAJkZ0+2NJEmri6e7SpKkRYmIXSLikIj4YER8NyIujIjLI+J3EfHDiDghIj4ZEc/uZinV2jg+IrK7VMssoF/bR8RBEfGRiDgjIi6NiCsj4scR8e2IOCIi7hURBhrGICJuFBGnRMQvI+KQafdnJYuIg7v99b8jYvdp90eSJC0NT3eVJEkLEhE3BV4CPIX6mGJrYAdg7+7+44DTgdtNuF/rgGcDLwWuXymyW3e5C/Bi4MSIeHFmnjTJfm0G/ha4Y3f7zcBRU+zLSvdW4Lrd5bnAYdPtjiRJWgoG6SRJUi/dzLNnAW8D1nWbLwW+DBwPXAD8EtgGuEF3uRvwSOCTE+7bTsCngAd0m64AvgB8qevX74EbAncFngDcFLgncEJEvBD453QtkIX6xdDtjbUCEbEjcH/gQcDhmfmzpejYCvQLYOfudvW9lCRJq49BOkmSNLIuQPcOykw1KEGvlwLvy8zL5qj6rm6G25YT7Nv2wDeA23SbvgI8KzMvmFX0VODYiHgF8HzgNZRg41sps/+OmFQfV7mjgCcBNwL+vlHm88ys3fe6pejUCvVq4J3AeZTvmyRJ2gwYpJMkSX0cxkyA7mzgEZl5zigVM/Nq4OpJdKoLHn6QmQDdx4CnZOb6OfqzHnhjRJwNHE0J1L0mIr6dmcdPop+rWWb+FLj9tPuxGmTmx4GPT7sfkiRpaZk4QpIkjSQibk+ZdQbl9NaHjBqgWwIPBA7obv8A+Ku5AnTDMvOLwCu7u2uAN5tMQpIkSUvNIJ0kSRrV64BB8OrQzPzRNDszy8uHbr8gM6/oWf91wEXd7X2Bh42lVxURcceIeGtE/E9E/CoiroqIn0TE1yLi+RFxrXnqn9llwT21ux8RsX9EfKLLpnt5RJwTER/qknuM2q+7RMR7I+J/I+KKiPhpRJzU9WnHrswgC++5lfr7D2XoPbK2nZlTXQF+NFQ+I+LwoTp7DG0/foS+nz/0HHOV2yoi/joivtq9visj4oKIOLbLqLp23jfqmu3tHhGHR8R3IuIXXXvnR8QXu8zCCxprR8SRQ69//8rjk9oH7hER74uIc7s2Lu0y9v5jRNyg52tYFxFPiYgvRMn4fGVEXBIRJ0fESyNil0qdP4uIjd1r+3mtzFDZQ4beo9MjYts+/ZMkaTkySCdJkuYVETcBHtzdvQD4yBS7cw0RsSclEQTA/wL/3reN7lTcfx3a9KgxdO0aIuI6EfFp4BTgeZQst9ehnGZ7Q+DPgbcA50fEIXM0NThleIeI2AY4EvgP4PGUTLpbAzcDngycFRGt9eEG/do6It4PfBt4OrAXsBUl4cfduz6dGhG3HqrWNwg6dRGxH3Am5XN+IOX1bQncGHgo5RTprwPXHqGt6N7Xs4GXAXcCrte1dxPg4ZTTVY+PiFqG4cUa9z6wbUQcBZwI/BUlocrWlOyydwReAXw/Iv5slM5FxL6UTM5HAo8Adqe8N9cH7kyZuXpmRDx0uF5mHge8tru7M/DBiE1ntXbHo3d2d/8IHDjPmpiSJK0IrkknSZJG8fCh21/MzA1T68mmHjR0+yuZudBsmF8GXtXdfsBcBfuKiJ0pAZCbdZvOB/4Z+A5wGSUwdhDwWErg7qMRsVtmvr7S3FXd9Y7Ae4AnUpJhvI1yqu9NKUHGx1PGem+MiHMy8wuVfm0BfBEYBF/+CLydErC7GNiHksn3bl2533bl+gTpTu7aAfgwJUhD95w/GSp3aY82e4mIO1MyDw9mW50OfJQSML2aspbe3wH3YeYzmstrgH/obp8FvA/4HvAryvv/NMp+eW/gSxFxjzF/Z8a5D2wFfI2S5RjK/vhuyj5wXtfWgZR98xhmZtNWde/114HtKfvJO7u2zgF2ogSjX9DdPiYi7pWZpww18TLKjMt7Ag+hfC5vG2p/LWU/Gsw4fWZmnjVXnyRJWjEy04sXL168ePHiZc4LJRCQ3eVJE2j/+KH29+hZ911DdZ+6iD6sBdZ37awH1oyj/5SgxjeGHv8ScJ1GO4+mBI0S2Ajcr1LmpKG2khIgWlt5zn8aKnMmEJW2XjVU5gzg5pUyQZndNPyc36qU23/o8SMX8zkDewyVO36E9//8QfnKY9sBFw61925g60q5dZQMtTnX66DMuhs8/t5GW2tmvbcH9NwXjxyqu/+E94E3zSpzq0afHjm0b7be6+0pgb2kzGrdu9HWXYFfD74PlcdvTAl4JnAlcLuhxw4d6sN7+n7PvXjx4sWLl+V88XRXSZI0il2Hbl8yX+Eoa5pl43L8mPs2fDrhzxfaSJaZTr/s7q6lnOo3Do+hzNCCEkx6TGb+utGHzzJzul8AR1SKDWfIPZMyk+gas7QyMynr7F3cbdoHuOVwmW6NsRd2dzcCT8jMH1b6lMD/oQQaB0ZKyrFMPIdyuiWUWX3PycqahVlOeX4KZWZcVXfq5eDzOQ34m0ZbG4HDgcHafYcttPMN49oHbgQ8t7u7EXh8Zp5Ze8Iss/AOn6dfzwD27G4fnI3EMpn5bcpp1AAPi4jbzHr8QuCp3d0tKTNLt4qIWzKT5OV0ymnjkiStGgbpJEnSKJbzEhk7Dt2+qllqNMOJAxbb1sBTh26/LDPna/fNzJxWeueImH365XByhHfNDs78qVB5ns8ObbrzrCKPA7bpbn84M7/X6lAX8Hlb6/Fl7klDtw/NOU6HzpIR+J2tx4H9KGsJAvzrXJ9l19ZXurt3iojtR+zvKMa5D2zZ3f5QZp4+z/P+y6znnu1p3fW3M/Pkedo6duj2fWc/2AUFB/vcbSnB0Q9Q1kv8A2UdusvneQ5JklaU5TzgliRJy8fPhm43My4OuT3XXLtqN+C4sfZoxi+Hbi949lu3PttO3d31wO8X06muzQDuP7Tpy/PVyczfRMQ3KKcXQlmba5Nsqp0fzNPcBUO3ZycwGF537/j5+kVJkrCidEkbbtvdvYqyLuB8Lpjjsf2Hbr8rIt7Vozs3o8y+G7fF7APD++bX53uizPxtRPyGsm7iNUTE9YDBjLi7xjyZdmfZu7H9MOBelOQVLxja/szMXHH7oyRJ83EmnSRJGsWPhm7Pno2zicw8OzPPGlwo61NNynAA8XbNUvO7AzOBxTPmmnHVw7WZSVbwW0ZPjjAcgNi1Waqs1zWX4efbYdZjuw/dHmXh/ZWYPXO3odvntGaczTLXTMcbL6Ivkxp3L2YfuNHQ7VGTL7Ten90b20dRfW8y80pKQpXh03s/m5kfW8RzSZK0bDmTTpIkjeKrwMu72w+LiOePGPBYCicyszbVg4CXLrCd4Syx31pUj2ZstcB6Ww/dXkxwbK7ZTMMBm98t4jkmbe38RZq2Hbo96qmRc2Uv3Wbo9vMp34tRndej7DjNtXR2iTIAAAhhSURBVA+M8/0Zfm++QckKPKpfzvHYfShJPQYeEBF7ZOb5PdqXJGlFMEgnSZJGcQrwY8rMmz0pa1l9fKo9mvFVykybdZS1v+7aLUw/sohYxzWDCh8eU98upSzIv4aydt5ulPdxPrcduj1K+YW4lHIKJt31fKdNjiuRxij+OHR79imaNVs2tg/PIts7IqJbX28um5zK2Whvx26W6Ep2KTOnmu4NfH+uwhGxhjI7tNXWwI3H8d5ExD7A27u7g+PPtSiJJPbv1v2TJGnV8HRXSZI0r27W3KuGNr0pIm44rf4My8zfURaUH3hbt75cH4cyc+rf1ylByXH0bT3wX0Ob7j1fnYi4FrBvd3cD8M1x9KVi+DXuP0L524/hOYdPIZ5rhtxw9tvduuBQVURch/YpwecNtbUj1wx+ttx1jse+O2K5leI7Q7fv0yw14/a0A6LnAb/pbu8ZETsvpmMRsQ3wScpsv8uB+zGTbOKelIzDkiStKgbpJEnSqD7ATJBiV+BLEbHbHOWX0suZSfRwV+DtcwV2hkXEw5g5lfcy4LkjzLbq4yNDt18XEdvNU/7VzGSsPTYzR13Hrq+jh24/OyKaa99FxFquuXD/Qg2vH9icsdYFNwfZZrcHHjhHm08fvtMl6xhuZzi76avn6lyX/OCpw5tmFfl3Zvazh3b7zkr2maHbz46IGzVLFi8avlN5rz839PC/DD++AG9hJqj6fzPzXMps18H7/08RcfdFtC9J0rJjkE6SJI0kM68GHgVc3G3aF/huRDxxATPXxiozfwocSJl5BuWP+U/NE3jaMiIOpQQWtqBkdD0kM+c77bOv9zNzKumNKafq1bJjromI5wJ/0226kjLDb1K+wcwsvW2Az9VmR3YBuvcxk7lzMS4aun3fecoeO3T7ZRExO+kBEfFQ4JWzNq+bdf+1zCQ7eHhEvLJ7TbPbunb3nNca2nyN/bqbtXnE0KYPRMScr2Pa3425ZOYJwP/r7m4NHB0RmySAiOIw4AmzHpr92l7FTCKLA4HDI6K5LmPX7ibvT0QcCDyzu/ufwD93/f0x8MJu+1rgqG7mqSRJq8KyHTRIkqTlJzN/EhF3o5yGdndgF8pMsVdGxBeBE4CfUjKZXpuyntitKGvYjeqvI+JXI5T7fGaeM9S3r3Z/3H+IkhThscCDIuLTwFeACyhr1+0K3IuSNXKPrvrPgCdl5nE9+jmSzLw6Ih4JnER5P/4CuEtEvINyyukfKJllD2LmdNgNwOMz84fj7s9QvzIintz14XrAXYBTIuJNlMDIekq23OcDtwaOA/6sqz5qkoHZPgO8uLt9eET8BjiV8pmcmpkXD5V9G/Bcymd5N+CbEfEu4DTK+3gg8CTK/nYC5f2DEmz6UwbSzDw3Ip7FzCnRLwXuFxHHUGaGRtf+33fvw3sp++uOXDOBx8ARlH3/YcDOwPERcTRlHzsLuKJ7PTeknEZ8n4i48TJKtDLbUyinve5K2QdOjYj3UPaLi4BbUgJm96QE6L8BHNzV3ZqhzKuZ+b8R8ZfAUZT39Z+AJ3TtnQn8hPIe34AyS+4A4DCGZvRFxJ6UzwDK2oRPnvXevR94OOUfBnsC76DsB5IkrXgG6SRJUi+ZeVE3e+jvKKdA7kYJdv1td5nPJfM8/pIRu3IucM7whsz8bER8n3Kq3EMpp0r+ZXepuZQSdHlnZi4mi+qcukDRnSjBi3tRAjitUy/PBZ6SmSdNqj9D/bogIu4FfB64BSVQ88ZK0Y9QZqSd2d3/daXMKM/3nYj4FCUItj1lht7AY4Fjhsr+PCIOAD5FCZjdAfjXWU2eTQnWPHlo2yYztzLzgxFxVVd/O0rA6Z6VLr6Xsl8Pgsq1ttZHxGMo+83fUc5MeWx3abkdJRi57HSB93tQTgu+AyVBSO07+ANKYPQZQ9u2Yub000F7H++Crx+kBPFvBrx+ji7cmy5IFxFbAp9g5nTvF3anuQ63nxHxdEpAcVfgiRHx5cz82AgvV5KkZc0gnSRJ6q079fVNEfF2SjDsPsA9KMGn61D+eP85ZZbTzygzcL4NHJeZF1UbHV/ffgg8LCLuQAmcPIASRLw+m459zqcEEh4dEacCP5xUxsjMvDAi7kNZX+1gSrBuV8qsrzO7y78Dxyxl1srMPDsibkcJZD6eMsPpWpRZT6cA78nM4yLiAUPVLt60pZE9gRKwOoAyy/IqyuewyezJzPxaRNyMMvvtgZSZU7+kBDKPAj6amZdHxO+GqlVPr8zMoyLiOMrsvIcDe1H2hx9TZjm+LzNPBBhaSq02k47MvAp4QUT8C/BXlH1sL8rs0V9R3p+LgTMop9CePs97MlWZeX4XRH4CZd/cD9iJElA/m/JefyIzL5u1zFzr/flyNyPuSZQZh/tRZtCtZ+a9OZ8yO/MrQ1VfRQm+AfwbMzPqZrd/aUQ8Bfhat+ldEXFSZp4/+quWJGn5ifGuiyxJkrQ8deuQ3YIyg+qhlFM3t59V7ApKwoLTgEO7NcgERMQLgDd3dw/IzKPnKi9JkqR+DNJJkqTNUreg/b0pAbuHUgJ4A+cDe405y+uKFRFbU2ZC3o6yHt3umfnL6fZKkiRpdfF0V0mStFnKzCspp9sdB7wwIvaiBOseBpyyuQToIuJWwNmtxAZdgO4dlAAdlNNCDdBJkiSNmTPpJEmSNlNdAO6nlCyan6AkBzifsq7ajShZXZ9LSQ4C8D/AvTPz95s0JkmSpEUxSCdJkrSZiog7A19n07X5ao4GnpOZ82XnlSRJ0gIYpJMkSdqMRcT2wEHAI4Ddu8sOwIXABcD3Kdk9N5tTgCVJkqbBIJ0kSZIkSZI0ZWum3QFJkiRJkiRpc2eQTpIkSZIkSZoyg3SSJEmSJEnSlP1/yKc/qXXnvCQAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "image/png": { "height": 333, "width": 628 }, "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "grid = coeffs.expand(grid='GLQ')\n", "fig, ax = grid.plot(show=False) # show=False is used to avoid a warning when plotting in inline mode" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Next, we would like to calculate the variance of this single spherical harmonic. Since each spherical harmonic has a zero mean, the variance is equal to the integral of the function squared (i.e., its norm $N_{lm}$) divided by the surface area of the sphere ($4\\pi$): \n", "\n", "$$N_{lm} = \\int_\\Omega Y^2_{lm}(\\mathbf{\\theta, \\phi})~d\\Omega$$\n", "$$Var(Y_{lm}) = \\frac{N_{lm}}{4 \\pi}$$\n", "\n", "When the spherical harmonics are $4\\pi$ normalized, $N_{lm}$ is equal to $4\\pi$ for all values of `l` and `m`. Thus, by definition, the variance of each harmonic is 1 for $4\\pi$-nomalized harmonics." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can verify the mathematical value of $N_{lm}$ by doing the integration manually. For this, we will perform a Gauss-Legendre Quadrature, making use of the latitudinal weighting function that is stored in the `SHGrid` class instance." ] }, { "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.4" } }, "nbformat": 4, "nbformat_minor": 1 }