{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Predicted ionospheric range delay (for different SAR bands and range bandwidths)" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Go to directory: /Users/yunjunz/Papers/2022_GeolocRg/figs_src/iono\n" ] } ], "source": [ "%matplotlib inline\n", "import os\n", "import numpy as np\n", "from matplotlib import pyplot as plt, ticker\n", "from mintpy.simulation import iono\n", "plt.rcParams.update({'font.size': 12})\n", "\n", "work_dir = os.path.expanduser('~/Papers/2022_GeolocRg/figs_src/iono')\n", "os.chdir(work_dir)\n", "print('Go to directory:', work_dir)" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "center incidence angle on the ground : 42.0 deg\n", "center incidence angle on the ionosphere: 38.7 deg\n", "--------------------------------------------------\n", "20 TECU in L-band --> refraction angle of 6 deg --> 27% overestimation\n", "20 TECU in S-band --> refraction angle of 20 deg --> 20% overestimation\n", "20 TECU in C-band --> refraction angle of 29 deg --> 12% overestimation\n", "20 TECU in X-band --> refraction angle of 35 deg --> 5% overestimation\n" ] } ], "source": [ "# constants\n", "K = 40.31\n", "Re = 6371e3 # earth radius in meters\n", "speed_of_light = 299792458 # speed of light in m / s\n", "\n", "# SAR parameters\n", "bnames = ['L', 'S', 'C', 'X']\n", "freqs = [1.257e9, 3.2e9, 5.405e9, 9.65e9] #carrier frequency for NISAR-L/-S, S1 and TSX\n", "rg_bandwidths = [24e6, 44e6, 80e6] #range bandwidth (sampling frequency) in Hz\n", "inc_angle = 42 #LOS incidence angle on the ground in deg\n", "height = 747e3 #NISAR altitude in meters\n", "\n", "inc_angle_iono = iono.incidence_angle_ground2iono(inc_angle)\n", "print('center incidence angle on the ground : {:.1f} deg'.format(inc_angle))\n", "print('center incidence angle on the ionosphere: {:.1f} deg'.format(inc_angle_iono))\n", "\n", "# Numbers for 20 TECU\n", "VTEC = 20\n", "\n", "# refraction angle\n", "print('-'*50)\n", "for bname, freq in zip(bnames, freqs):\n", " ref_angle_iono = iono.iono_incidence2refraction_angle(inc_angle_iono, VTEC, freq)\n", " rg_delay_inc = (VTEC / np.cos(np.deg2rad(inc_angle_iono)) * 1e16 * K / (freq**2)).astype(np.float32)\n", " rg_delay_ref = (VTEC / np.cos(np.deg2rad(ref_angle_iono)) * 1e16 * K / (freq**2)).astype(np.float32)\n", " perc = (rg_delay_inc - rg_delay_ref) / rg_delay_ref * 100\n", " print(f'20 TECU in {bname}-band --> refraction angle of {ref_angle_iono:3.0f} deg --> {perc:3.0f}% overestimation')\n" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "save figure to file /Users/yunjunz/Papers/2022_GeolocRg/figs_src/iono/iono_refr_impact.pdf\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAqwAAAFHCAYAAABkqR7sAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy89olMNAAAACXBIWXMAAA9hAAAPYQGoP6dpAADp20lEQVR4nOydd1gU19eA3y30KqAgKs0CNkSw995N+5KYmGJiqib5pffElm4SY4opJjEm0VRjTIy9YDRWLIgVFBSwIEWkly3z/TGysALKLgu7i/fl2YedOzN3zpwpe+bOKQpJkiQEAoFAIBAIBAIbRWltAQQCgUAgEAgEgqshDFaBQCAQCAQCgU0jDFaBQCAQCAQCgU0jDFaBQCAQCAQCgU0jDFaBQCAQCAQCgU0jDFaBQCAQCAQCgU0jDFaBQCAQCAQCgU0jDFaBQCAQCAQCgU0jDFaBQCAQCAQCgU0jDFaBQCAQCAQCgU0jDFaBQCAQCAQCgU0jDFaBQCAQCAQCgU0jDFaBQCAQCAQCgU0jDFaBQCAQCAQCgU0jDFaBQCAQCAQCgU0jDFaBQCAQCAQCgU0jDFaB2SxevBiFQsHevXsbZb2GRqFQMGvWLGuL0aTYvXs3N998M0FBQTg5OeHv70/fvn159tlnr7nukCFD6NKlSyNIWXcqzt3Tp09bW5QmS0JCAvfffz+hoaE4Ozvj7u5OdHQ0c+fO5eLFi1ddV5wz1wdTp07FycmJQ4cOVZv37rvvolAoWLlyZa3rh4SEMGHChIYU0WRmzZqFQqGwthg2jTBYBQJBg7Bq1Sr69etHfn4+c+fOZf369Xz88cf079+fX3/91driCWyQr7/+mpiYGOLi4nj++edZu3Ytf/75J7fddhtffvklDzzwgLVFFNgA8+fPJyAggClTpqDRaAzthw4dYubMmdx3331MnDjRihIKGgK1tQUQCARNk7lz5xIaGsq6detQqytvNXfccQdz5861omQCW2Tnzp1MmzaNkSNHsmLFCpycnAzzRo4cybPPPsvatWutKKHAVvD09OTbb79l1KhRvPnmm8yePRuNRsM999yDv78/8+fPt7aIggZAjLAKrEZubi73338/Pj4+uLm5MXHiRFJSUoyW2bBhAzfeeCOtW7fG2dmZdu3a8cgjj5CdnW20XMXrlCNHjnDnnXfi5eWFv78/U6dOJS8vz2jZ/Px8HnroIXx9fXF3d2fMmDEkJSU1+P5eb+Tk5ODn52dkrFagVNb91rNt2zb69OmDi4sLrVq14vXXX0en0xktM3v2bHr37o2Pjw+enp5ER0fz7bffIkmS0XIVrwLXrl1LdHQ0Li4uREREsGjRomrb3bVrF/3798fZ2ZnAwEBefvllo9EcgWV5++23USgULFy40MhYrcDR0ZEbbrihTn2Jc6bpM2LECB599FHefvtt9u3bx6xZszh48CDffvstXl5ederjzz//JDIyEmdnZ8LCwvjkk0+M5peWlvLss88SFRWFl5cXPj4+9O3bl7/++qtaXwqFgscff5wff/yRjh074urqSrdu3fjnn3+qLbtq1SqioqJwcnIiNDSUDz74wDwlXGeIEVYrMfHT/8gqKLO2GAA093Bi5RMDGn27DzzwACNHjuSnn34iPT2d1157jSFDhpCQkIC3tzcAycnJ9O3blwcffBAvLy9Onz7NvHnzGDBgAIcOHcLBwcGoz//7v/9j0qRJPPDAAxw6dIiXX34ZwPDjIkkSN910Ezt27GDGjBn07NmT7du3M3bs2Ebd9zrx1WAozLS2FDLuLeCRf01apW/fvnzzzTf873//46677iI6Orra8boWGRkZ3HHHHbz00kvMmTOHVatW8eabb5Kbm8tnn31mWO706dM88sgjBAUFAbLh8MQTT3D27FlmzJhh1OfBgwd59tlneemll/D39+ebb77hgQceoF27dgwaNAiAo0ePMnz4cEJCQli8eDGurq58/vnn/PTTTybJ3xhM+mcS2SXZ116wkfBz8ePXCaa5fOh0OjZv3kxMTAxt2rSp1/bFOXN1Tv3frWizbeN8Ufv5EfrHMrPXf//991m3bh233nor6enpPProo4wcObJO68bHx/PUU08xa9YsAgICWLp0KU8++STl5eU899xzAJSVlXHx4kWee+45WrVqRXl5ORs3buSWW27hu+++49577zXqc9WqVcTFxTFnzhzc3d2ZO3cuN998M4mJiYSFhQGwadMmbrzxRvr27csvv/yCTqdj7ty5XLhwwWw9XC8Ig9VKZBWUkZFfam0xrEqPHj349ttvDdOdO3emf//+LFiwgFdffRWARx991DBfkiT69evHkCFDCA4OZs2aNdVGXB544AGef/55QH4CP3nyJIsWLeLbb79FoVCwbt06YmNj+fjjj/nf//4HyK8bHR0dDdu0GQozoeCctaUwm3fffZfjx4/z6aef8umnn+Lg4EDPnj2ZOHEijz/+OO7u7tfsIycnh7/++stwnEeNGkVJSQlffPEFL7zwgsHY+O677wzr6PV6hgwZgiRJfPzxx7z++utGwQzZ2dls377dsO6gQYPYtGkTP/30k8H4mDNnDpIksXnzZvz9/QEYP368zQX0AGSXZJNZbCMPNmaSnZ1NcXExoaGh9e5LnDNXR5udjbaJGEdubm68+eabTJ48mYCAAN5///06r3vu3DkOHDhAt27dABg7diyZmZm88cYbTJ8+HVdXV7y8vIzOE51Ox/Dhw8nNzWX+/PnVDNaSkhI2btyIh4cHANHR0QQGBvLbb7/x0ksvAfDqq6/i7+/Phg0bcHZ2BmD06NGEhITURxXXBcJgtRLNPaq/8rIWDSWLVqs1mlapVEY/AnfddZfR/H79+hEcHExsbKzBeMzMzGTGjBmsWrWKc+fOodfrDcsfO3asmsF65XRkZCSlpaVkZmbi7+9PbGxsjduePHmy7Rms7i2sLUElZsji6+vLtm3b2Lt3L5s2bWLv3r1s2bKFl19+ma+++oq4uDh8fHyMjqlCoUClUhmmPTw8qh3TyZMn8/XXX7N161buvvtuADZv3szbb79NXFwc+fn5RstXHPsKoqKiDIYHgLOzMx06dCA1NdXQFhsby/Dhw43WU6lUTJo0idmzZ5usi4bEz8XP2iIY0dDy6PV6cc7UA7Wf7Zwv9ZVFr9fz6aefolQqyczM5ODBg/Tv398w72rnSefOnQ3GagWTJ09mw4YN7N+/nwED5LeOv//+O/Pnz+fgwYMUFRUZlq0wNqsydOhQg7EK4O/vT4sWLQznSVFREXFxcUyfPt1ofQ8PDyZOnMj3339fH3U0eYTBaiWs8Qq+MTl9+nS10ZLY2FiGDBlimA4ICKi2XkBAADk5OYB8wxk1ahTnzp3j9ddfp2vXrri5uaHX6+nTpw8lJSXV1vf19TWarvCFq1g2JycHtVpdbbmaZLE6Jr6Ct1V69OhBjx49ANBoNLz44ot89NFHzJ07F1dXV6Mf8+DgYKP0P1V//CuoOFYV58mePXsYNWoUQ4YM4euvv6Z169Y4OjqyYsUK3nrrrWrnyZXHHuTzpOpyOTk5tZ6ftoapr99tET8/P1xdXTl16tQ1l50zZ444Z+pBfV7B2xoffPABO3fu5JdffmHmzJlMnTqV+Ph4XFxcmDp1qpEBOHjwYLZs2WKYvtqxqjhPli9fzu23385tt93G888/T0BAAGq1mi+++KJGH+ZrnSe5ubno9Xq7OE9sEWGwChqEwMBA4uLijNrCw8ONpjMyMqqtl5GRQbt27QA4fPgwBw8eZPHixUyZMsWwzMmTJ82Wy9fXF61WS05OjtHNpSZZBJbHwcGBmTNn8tFHH3H48GG++eYbo3yIVwbb1OTXVXGsKo7fL7/8goODA//884/RqMWKFSvMltPX17fW81NgeVQqFcOHD2fNmjWcOXOG1q1b17rsww8/LM4ZAUePHmXGjBnce++9TJo0ieDgYPr378+rr77KvHnzmDVrFo8//rhh+aojn1D77w9UnidLliwhNDSUX3/91ejtYFmZefEnzZo1Q6FQiPPETESWAEGD4OjoaBhZq/hcecNYunSp0fSOHTtITU01jMJW3CCu/EH66quvzJZr6NChNW7b1gIjmgLnz5+vsf3YsWOA/FATGBhodI507drVaNmCggL+/vtvo7affvoJpVJp8B1UKBSo1Wqj130lJSX8+OOPZss+dOhQNm3aZGT86HQ6kT+2AXn55ZeRJImHHnqI8vLyavM1Gg0rV64U54wArVbLlClT8PPz4+OPPwagT58+PPPMM3z88cds376dkJAQo/PkygGTI0eOcPDgQaO2n376CQ8PD6KjowH5PHF0dDQyVjMyMmrMElAX3Nzc6NWrF8uXL6e0tDKGpaCg4KqFDgQyYoRVUG82b95cYxWXcePG4erqWut6e/fu5cEHH+S2224jPT2dV199lVatWjF9+nQAIiIiaNu2LS+99BKSJOHj48PKlSvZsGGD2bKOGjWKQYMG8cILL1BUVESPHj3Yvn17vX6oBDUzevRoWrduzcSJE4mIiECv1xMfH8+HH36Iu7s7Tz755DX78PX1Zdq0aaSlpdGhQwdWr17N119/zbRp0ww+hePHj2fevHlMnjyZhx9+mJycHD744IMaUyPVlddee42///6bYcOGMWPGDFxdXVmwYIGRD5vAsvTt25cvvviC6dOnExMTw7Rp0+jcuTMajYYDBw6wcOFCunTpcs2E8OKcafq888477N27lzVr1hgyygC88cYbrFy50sg1oDYCAwO54YYbmDVrFi1btmTJkiVs2LCB9957z/C7NWHCBJYvX8706dMNmQjeeOMNWrZsyYkTJ8yS/Y033mDMmDGG3MI6nY733nsPNze3a1Zyu+6RBAIz+e677ySg1s+pU6euut769eule+65R/L29pZcXFykcePGSSdOnDBa9ujRo9LIkSMlDw8PqVmzZtJtt90mpaWlSYA0c+ZMw3IzZ86UACkrK6vGbVWV5dKlS9LUqVMlb29vydXVVRo5cqR0/Pjxan0K6sevv/4qTZ48WWrfvr3k7u4uOTg4SEFBQdI999wjHT169JrrDx48WOrcubO0ZcsWqUePHpKTk5PUsmVL6ZVXXpE0Go3RsosWLZLCw8MlJycnKSwsTHrnnXekb7/9ttqxDw4OlsaPH1/jtgYPHmzUtn37dqlPnz6Sk5OTFBAQID3//PPSwoULr3puC+pPfHy8NGXKFCkoKEhydHSU3NzcpO7du0szZsyQMjMzr7quOGeaPvHx8ZKDg4P00EMP1Th/586dklKplJ5++ula+6g4psuWLZM6d+4sOTo6SiEhIdK8efOqLfvuu+9KISEhkpOTk9SxY0fp66+/NvzeVAWQHnvssRq3NWXKFKO2v//+W4qMjJQcHR2loKAg6d13362xT4ExCkm6IkuyQCAQCAQCgUBgQwgfVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTSMMVoFAIBAIBAKBTWNzBus333yDQqHA3d292rz9+/czYsQI3N3d8fb25pZbbiElJcUKUgoEAoFAIBAIGgubMljPnj3Lc889R2BgYLV5x48fZ8iQIZSXl/Pbb7+xaNEikpKSGDhwIFlZWVaQViAQCAQCgUDQGCgkSZKsLUQFEydORKFQ4OPjw7JlyygsLDTMu/3224mNjSU5ORlPT08AUlNTad++PU8//TTvvfdejX3q9XrOnTuHh4cHCoWiUfZDIBA0TSRJoqCggMDAQJRK85/3xX1JIBBYCkvdl2weyUb48ccfJQ8PDyk9PV2aMmWK5ObmZpin0WgkFxcX6ZFHHqm23qhRo6T27dvX2m96eroEiI/4iI/4WOyTnp5er/uduC+Jj/iIj6U/9b0v2TpqbIDMzEyeeuop3n33XVq3bl1tfnJyMiUlJURGRlabFxkZyYYNGygtLcXZ2bnafA8PDwDS09MNI7NXQ6PRsH79ekaNGoWDg4MZeyOoitCnZRH6tBzm6DI/P582bdoY7ivmIu5L1kXo07IIfVoWU/VpqfuSrWMTBuv06dMJDw9n2rRpNc7PyckBwMfHp9o8Hx8fJEkiNzeXli1bVptf8bpNq9Wi0WgAUCqVuLq6olar0ev1hmWVSiV6vR5XV1dcXFwMJ4pKpUKpVKLVapGqeFBUtFf0W4FarTZssy7tDg4O6PV6dDqdkdwV8tXUrtPpqsmuUqlqba9N9obep9LSUoM+HR0dm8Q+WfM4qdVqXF1dcXd3x8nJqUns05XtjbVPFbr09PQ0XPtXyl5cXExJSYmhvWLd+r7Gr1jf09OzzgZrhazCIKg/Qp+WRejTspirz6buXmR1g/WPP/5g5cqVHDhw4JrKvtr8a60bGhpqND1z5kxuvPFG0tLSDG3h4eG0b98eZ2dnNmzYYGiPiooiODiYrVu3UlBQYGjv27cvLVq0YP369UY/mkOHDsXFxYXVq1cbbXPcuHGUlJQQGxtraFOr1YwfP57s7Gx27txpaPfw8GDYsGGkp6cTHx9vaG/evDn9+vXjxIkTJCYmGtqDgoLo3r07CQkJ1fYpIiKCPXv2GAWnNfY+bdiwocntE1jvOF28eJGWLVs2qX2yxnFydnZGqVTWuk9PP/00CxcuxNoolUqCgoKatn9aIyL0aVmEPi2L0GfNWDXoqrCwkHbt2nH33Xfz2muvGdqnT5/O33//zZkzZ3BwcODMmTNERESwYMECpk+fbtTH888/z4cffkhxcXGNLgH5+fl4eXlx6tQpw3D51UZYbXlEqGp7UxnlEvsk9smW9+nKEdbCwkJCQkLIy8ur08hobVTcl+rbj0AgEFwv9xOrGqynT5+uNvJ5JTfeeCPLli3D09OTKVOm8MUXXxjNHzNmDCkpKSQlJdW4vqkHUqfTkZCQQGRkJCqVqu47I6gRoU/LIvRpOczRpaV+GMR9yboIfVoWoU/LYqo+rxeD1arjzQEBAcTGxlb7jB49GmdnZ2JjY3nzzTdRq9VMnDiR5cuXG71CTEtLIzY2lltuucViMun1etLS0oxGWwTmI/RpWYQ+LYc96dKeZLUHhD4ti9CnZRH6rBmr+rA6OzszZMiQau2LFy9GpVIZzZs9ezY9e/ZkwoQJvPTSS5SWljJjxgz8/Px49tlnG09ogUAgEAgEAkGjYvWgq7oSERHBli1bePHFF7n11ltRq9UMGzaMDz74gObNm1tsOxlr/yI/MwVJGmuxPgUCgUAgENgvkiShlyr/6yUJSQKJK6arLCchtyMht1FlnoTR+hXLS5JEuUbL+WIo0+oRSRcqsUmDdfHixSxevLhae0xMDBs3bmyw7Zbl55H34uv00EPS98vxGTAUt759cO3dG4eAgAbbblNGqVQSHh4uoh0thNCn5bAnXdqTrPbA9axPvV6iXKdHo9Oj1UlodHo0egnt5TaNTpLb9Xp0esmwnE4vob28nFZvPK3R6jivDyBrZxp6QKcHnV5eTn95OV3FR6r8rjd8p/K7JK9jNP+yQVexniQhL1exrCShv9yH/gqDsqIfvb6yvdLwNF626vpXGqaNj5qhQ0oJd3GyxsZtEps0WK3F34sepctllxFVdh55K1aQt2IFAI7Bwbj26YNb71649u6N2tfXeoLaESqVioiICGuL0WQQ+rQctqDLnj17olKpeOyxx3jsscdqXc4WZG1KWFOfer1EiUZHcbmOUo2OEo2OkvLL/zU6Sst1lGp1lGr0lGp0lGn1lGn0lGp1lGn0lGl1lGv1cvvl7+U6eZlynd4wrdHp0Whlg7O8ijGq0zek9ZXdgH1ff1zrgWrBggUsWLDAKPNJU8aqWQIaA1Oi5/IP/suK+Q/idN6BiDMSjtral3Vs1xa3Xr1x7d0b1149UTdrZmHJmwZarZY9e/bQq1cvQxoigfkIfVoOc3RprSwB4rhbFnP0WabVkVeiIb9EK/8v1ZBfoqGwTEtBqZaCUg2FpVoKy3QUlWkpKtdSVKaluFxH4eX/RWVayrQikMaSKBSgUihQKhTyd6X8XamQ09ZVnVZenlZc/l7xX2n4L7cpqrRduawCUCov/69pvoIq84ynlUpQIDcY+rq8PQUYlkeSyM7O4q07+tLKx/2aOrhesgSIO18VPLsN5q5RvViUsZX7PTxpd15B19N6epx1ISS9DKo8xZSfTKb8ZDK5P/0EgFOHDrj26oVrr5649hQGbAWSJJGVlUUTfy5qNIQ+LYc96dKeZLUH9Ho9qeez8Mks5FKpjuzCMnIKy8kpKudiURm5RRouFpWTW1xOXomGS8UaSjS2P4qlUICjSomjWomTWomDSv44qpWolQrD/4o2eb4CtUqJg/Lyf5UCtVKJWqVArTSep1LKbarLfVRMS5KeI4cSiOkehZOjAyqlApVCgepyHyqFAmWVdSuMyGrfLxt1qirrVPxXVjFGrzRKmxoajYbVq1fTwkO4A1RFGKxXIPWZzkPf/0P7cg0vtWnOsSAVv1GOH568630foScLKNq9h9LDh6FKyomypCTKkpLIXbIEAKf27XHt2bPSgBUuBAKBQNDgSJJETlE55y6VcO5SCWcvlXLuUgkZeaVcyC/lQkEpF/LLKNeqIW57g8vj7KDEzVGNq5MKV4fL/x1VuDiocXFU4eqgwsVRhbODCmcHpfxfffn/5TYnBxVOl9scVUqcHZQ4qlQ4OSgNBmqFMWoNA06j0eCacZBxkS1FaVZBg3HdGKx19RXTBvYk16UtQ0qS+ensOf7XrhuppVlkU8BDlz7jsRGP8dAzPyMVFVOybx9Fu/dQvHs3pceOGRuwJ05QduKEYQTWsW1bXHv0kI3Ynj1x8G/R4PssEAgsw/XmK2brSJJEdmE5yVmFpGQVkZJVSNrFYtIuFpN+sZiicsscJ0eVEm9XB7xdHfByccDLxRFPFzVeLg54ODvg6azG09kBD2c1Hs4OuDurcXdS4eaklj+OalTKpjcCKBBYA+HDWoV9qRe57cudjFHs4nPHTwAoaDuMV1u1ITa9sl75kNZDeGvgW3g6VvanKyigeN8+ivfEUbxnD6VHjxoZsFfiEBwkG7AxPXDt2QOH1q2b5KsNvV5Peno6bdq0uS4jci2N0KflMEeX1vJhvZ6Pe1GZluMZBRzPyOfY+XyOny8g6UIB+aVXCTK4Ct6uDrTwcMLTQaK1nxctPJ3xc3fE180JH3dHfN0caebqiI+bI66OqiZ5X7Y01/P52RCYqs/rxYdVGKxVOJlZyIh5/6JCR5zH8/hoMgDQT9vJtxe28+mBT5EzpUFr99Z8OORDOvl2qrEvXWEhJfv2URwXR1FcHKWHjxj5wF6J2t8f15gYXHv2wCUmBqd27VCIC18gsCmsZbBeL+j1EkmZBRxIu0R82iUOnrlE0oUC6hrY7qBS0LqZK218XGndzIVW3i4EejvT0suFQC8XWng64ewgSocKmhbXy/1EGKxVuFRcTtScDQC84f8v9+R9Jc+Ivhdu+JTtZ7fz4rYXySvLA8BB6cBLvV7itg63XfMpXF9URPGBeIr3xlEct5fShAQkjabW5VVeXrhER+PaIwaX6GhcOndG4ehowp7bBlqtlq1btzJo0CAR3WwBhD4thzm6tGaWgKZ43CVJIvFCAbuSc9iZksPuUxe5VFz7fbGCll7OtGvhTtvm7rRt7kZYc3eCfV1p6eVSp1fwTVWf1kLo07KYqs/rxWAVZ1YVvFwccFAp0Ogk/pSGco/Tz1CWDwd/hWEz6N+qP79N+I3n/n2OQ9mH0Og1vLHrDfZe2MvMvjNxc3CrtW+lmxvuA/rjPqA/APrSUkoSEijeu5eSvXspjj+IVFxsWF6Xl0dhbCyFsbIrgsLJCZfISFxionGNicElKgqVh0fDKsQCSJJEQUGBiG62EEKflsOedGlPsl6LUo2OnSk5bDp2gc3HMjmXV1rrsiqlgnB/DyJbe9GxpScRAR5EBHji5Vq/wJ6mpE9bQOjTsgh91owwWKugUCjwcXPkQn4ZacVq6Hkv7PwMdGWw+wsYPoNA90C+H/M98/bNY8kxOSPAmlNrOJpzlPcHvU9H34512pbS2Rm3Xr1w69ULAEmjofTYMYr37qN43z5K9u1Dd+mSYXmprIziuDiK4+LIkYWVU2nFROPSPRrXmGgcAgMtrBGBQCCoP1qdnm0nsll+4Cybjl2guJagKE9nNb3DfOkV4kNUkDddAr1wcRSv8AUCQR0NVnOGmBUKBVu3bqVbt24mr2tN/Nxlg/ViUTm63tNQ7f4K9BrYvRD6Pg6uPjioHHix14t0b9GdGTtmUKQpIjU/lbtW38WzPZ5lcsRkkx31FQ4O8ghqZCS+U+9H0uspP3WK4r37KNm/j+K9+9CcPVu5giRRlphIWWIiuT/9DIA6IACX7lG4do/GJToa54hwFOL1jEAgsBLHM/L5Le4Mfx88S3ZhebX5jiolfdr6Mqi9H33b+hIR4Cmi6gUCQY3UyYdVqVQybtw4mjdvXqdO9Xo9S5YsIS4ujujo6HoLWR9M9e2YsmgP/yZlAbD3tRH4xb4I+76TZw5+CYa+bLR8en46z299niM5RwxtQ9oM4Y1+b+Dt7G2x/QDQZGRQsn8/xfv2U7x/P2WJiVfNRKBwcZGN4O5RuHbvLrsReHlZVKZrodfryc7Oxs/PT0SPWgChT8thji6tmSXAXo67Xi/x74ksvt12iv9OVi/V6e3qwMiO/gzv6M/A9n64OTX+Q7U96dMeEPq0LKbq83rxYa2zwbpr1y56XX59fS20Wi2Ojo7s3bvX7gzW534/yLJ9ZwBY99Qgwp0uwqfRoNeCkxc8fQicjY0+jU7Dx/s/5vuj3xvaWri24K0Bb9GnZR/L7lAVdIWFlMQflI3YA/spOZhg5AdbE45t28oGbFQULt274xgaKrIRCAR1RGQJqB2tTs/y/WdZuC2Fk5mFRvMc1UpGdGzBzd1bM7hDcxzV4p4jEFiKpng/qYk6Pdo++eSTtGzZss6dqlQqnnzySfz9/c0WzFr4uFaqJLuwjPCAYOh2BxxYAmV5sGchDHreaB0HlQPP9XyOXi178dp/r5FblktmcSYPrX+IKZ2m8L/o/+GosnyEv8rd3SiQS9JqKU1MpORAvGzExh9Ae+680TrlycmUJyeTt+wPAJSenvIobFSU/OkWadFgLo1Gw/r16xk1apSogGIBhD4thz3p0pZllSSJ2MRM3l59vJqhGuLryv39Q7mpeyu8XGxHblvWpz0i9GlZhD5rpk4G60cffWRSpwqFwuR1bAU/98ravdmFZfKXgc9C/M8g6WDnAuj9KDhVN+oGtR7EshuW8dp/r7Hz/E4Avj/6PbvO7+K9Qe/R1rttg8quUKtx6dwZl86d4e67gMtuBPHxlBw4QPH+A3JFLm1lwm19fj5F//1H0X//Xe5EgWPbMFy6dbtswHbDqW1bFCrzAx+0WvMSfAtqRujTctiTLm1R1sNn83h79TF2JOcYtfcK9eHBAaEM7+hvsz6ptqhPe0bos47o9VCcA4UXrvhkyf+LMlEXXGBE7lkYeAAcRFXMCq6biJy6lmb1dascCc0quGyw+oRB19sg4RcoyYW4b2HAUzWu38K1BV+O/JIlR5cwf/98NHoNibmJ3L7ydp6KeYq7Ot6FUtF4r8McAgJwGDMGzzFjADmdVunhw5TEx1N8IJ6S+Hh0OVV+bCSJ8pPJlJ9MJu+P5YCcksu5a1fZiO3WDZdukah9fRttHwQCa9NQpVnrel+yNcq0OuZtSGLh1hSqOpV1D/Lm1XEd6RHiYz3hBAJrIEmyfZB/DgrOX/5kVPl/+VOUKbsYXgUF4ARoirLAs3aD9XorGW2ywbp169Za5ymVSry9vYmIiLC55MFxcXF18u3wc680WI2iWgc9Bwm/AhLs+BR6PQyOrjX2oVQoubfzvfRu2ZsXt75Icl4y5fpy5sbNJTY9ljf7v0mgu3VSUCmdneWSsD164Iv8Ok9z9qzsRhAfT8nBg5QeP248CltURPGuXRTv2mVoc2jVCpdukThHRuIS2Q3nTh1ROjtbYY8EgoanwqCs8BWzFHW9L9kSSRcKePKXeI6dzze0Bfm48tLYCMZ2CRClTAVND0mComzIPwN5Z2WjNP/M5f/nIP+sbIxqa88pbNLm1C4UK91x1Fw9JqWh7ku2ismVrpRK5TVvSO7u7jzzzDPMnDmzXsJZAlOdkY+dz2Psx/Lr8VtjWvPBbVXSci2bCodl309GvQn9nrhmf6XaUj7e/7EhZyuAm4MbL/Z8kZva3WSTN3d9aSmlR45cNmATKImPR5uZefWV1GqcO3TAuVskLl0jcekWiWNoKCgUFBQU4OHhYZP7am9UJJQW+qw/5ujSWkFXtnDc9XqJ73ee5p01xynXytlJHFQKnh7ZgQcGhOKktp98qbagz6aE3etTWy4bnZfSIC8dLqVD3hnIS5P/55+zgDGqALfm4OEP7gGV/939wb1Flf8tkBzcKCgsrLM+r5egK5MN1lWrVvHEE08QHh7OnXfeib+/PxkZGfz8888kJiYyZ84ctm3bxqJFi/joo4944olrG3UNiakHMquglJ5vbQJgSHhzFt9fJTNC5jH4vC8ggUszePJgtYwBtbH7/G5e3/4654sqg6AGtR7EjD4z8Hez/eA02Rf2ICUJCZQkHKT08BGk0qtfwLIrQRccO3XGrVs3XCK7og4QIzD1QZIktFotarVa6LGemKNLaxqs1jzuZVodz/52kH8SKu9fHfzd+WhSFJ0D7W9kx9r6bGrYvD71evnVfO5p+XMpFXJT5f+X0mSDlHpUlXL2Ao9A8GwJnoHg0bLKJ0D+uLUAVd3ePJuqT2Gw1sITTzxBQUEBixcvrjZvypQpODs789VXX/Hwww+za9cuEhISLCWrWZh6IEvLyuk4cz0SCrq08uSfJwYaL7D84cuuAcDA52D463WWpaC8gPf2vMdfyX8Z2jwcPHix14vc0PYG27zQa0HSaik7cUIegb1sxJYnp8A1TidVcz9cunTFuWsXXLpG4tylM+pmzRpJavtHo9GwevVqxo0bJ6JH64k5urSWwWrN415QquHRJfvYfrLS1/3+/iG8OCYCZwf7GVWtiriOLItN6FOnkY3PiylXfE7JhqmueuGKOuHoAV6twasVeLaSv3u2qpz2DATH2suym4Op+rxeDFaTHU1//fVXfvzxxxrnTZ48mbvvvpuvvvqKCRMm1LqcLaNSKnB3gAIN5NRQmYWhr8KRP+WTf+cC6PWQ/PRUBzwcPXhzwJsMDxrOnF1zyC7JpkBTwGvbX2Pd6XXM6DuDALe69WVtFGo1zh074tyxI83umATIeWFLDx+mJOEQpYcSKEk4hPbCBaP1dFnZFMbGUhgba2hzaNMGl65dcO7cBeeuXXDu1BmVu2VvAAKBwHSyCsq4f/EeDp+V/VWdHZR8dmc0IzrZ/lshQRNEr4eCc5B9AnJOXv4kw8VkecRUMiP4yNUXvIPBuw14B4FXkPzdq41snLp4W3w3BOZhssFaVFREVlZWjfMuXLhA8eXE9R4eHjYXeFVXPKoYrJIkGY98NguGHg/A7i9AWwL/vgcTTEvhNTRoKNH+0by35z1WpqwEYNvZbdz01008Hf00t4Xf1qiZBCyFyt0dtz59cOtTWSyh5MwZtv/wA12cnCg/cpSSw4fR5+cbradJT0eTnk7+6jVyg0KBY2gozl0649KlC85duuAcEYHSteYgN4FAYHnScoq5Z9FuUnPke7qXiwOL7utJTLB4IyJoYLTlsjGanSgbp1mJ8vecZLhGIFI11C7QLET++ITKxmmz4MtGahA4uTfEHggaAJMtyv79+/P666/Ts2dPwsPDDe3Hjx9nxowZDBgwAICUlBRat25tOUkbEQ8HCVBQrtOTX6LFy/WKIflBz8mFBMoLYN/30Ocx8Gtn0ja8nLx4e+DbjAoZxZydc8gqyaJIU8Sbu99k9anVzOo3i1CvUMvtlJVQ+/tT1LkzvpdfbUiShCY1lZJDhyk9fEj+f/SosT+sJFGekkJ5Sgr5f8sGPUolTm3DcO7UGefOnXHu0lkYsQJBA5FdWMbkb3ZxJrcEgJZezvwwtRft/S1XVEQgQKeVDdPMo3KMSNYx2TjNSTZttNTBDXzbyikoK/77hEGzUPkNqB252wlqx2Qf1qNHjzJo0CAuXbpEly5d8Pf358KFCxw+fJhmzZqxdetWOnbsyBtvvIGDgwMvvfRSQ8leJ8wJbnjqlwP8dVAOLtj4zGDatajhCezfuRD7lvy9041w+w9my5hXlse8ffNYfmK5oc1R6cjDkQ8ztctUHFT262NVF+dxSaulLDnlsgF7iNIjRyk7fhxJo7l650oljmGhuHTuIhuxnTvJRqxb03UnsPngBjtCBF3VTJlWx11f72Zvai4A7Vq488PUXgR6uzTodhsTcR1ZljrpszATMg7BhSNw4bD8Pzup7r6lCqVsgPq1B992lf9928kR9k3oOIqgq5ox2WAFyMjIYN68eWzdupWcnBx8fX0ZPHgwTz31lEklXBsDc34YZq44yA+7zwLwy8N96BNWQ5L8skL4pLucBBjgwc3QOqZesu4+v5tZO2ZxpvCMoS3MK4wZfWcQ41+/vq2FuelO9OXllCWdkH1iD182Yk+cMMoPWyMV7gSdO+PcqdPlT0eLlpu1JnafPsaGEGmtat7O88sSWLZPvge18HDi78cHEODVtHIsi+vIshjpU5Ig9xScj5cN1PMJ8v+ia6RGrEDlBH4doHk4NI+A5h3kaZ8wUDtde/0mgKnnpzBYmwjmROM++/Ua/k6To18/m9ydCZG1JPnf8zWsfk7+HtQP7l9d76e8Em0Jn8d/zo9Hf0RX5ZXI/7X/P56OeRovJ/tKIWPJ6FF9WRllSUmyEXvkSN2NWMAhKKjSgO3YEefOnVD72F81HpuIxm0iiCwB1Vm4NZm3Vx8HwEmt5LdH+tKtjXeDbc9aiOvIQkgSXExBmx7H6e1/EuaSj/LCYSjLv/a6CpU8StqiI7ToVPnfO7jO6Z+aKiJLQM2YfVbk5eWxa9cusrOzGTduHM1sPDWRKSUQPaqcH9kV5VlrIuY+2PWFHKGYtgOOLIcu/1cvOV3ULjzb41nGhY5jzs45HM45DMAfJ/5gc9pmno55mhvb3WiXQVn1RenkhEvXrrh07UrF2WYwYo8cpfTIEUqPHqU0KQmucCfQpKWhSUujYO1aQ5va3x/njh1x6hghG7GdOuHQqpUYcRFU43oogbjp2AXeWXPcMP3h7d2apLEqqAeFWXB2L5zZC2f3wbkDUHoJNXDVKA4XHwjoKn/8O4N/F3kE9ToZMRVYBrMM1jfeeIN3332XkpISFAoFcXFxNGvWjOHDhzNy5Eir+63WhCklEI0M1ppSW1WgcoAx78BPt8vT62dAhzEWycnW0bcjS8Yt4dfEX/nkwCcUaYrILctlxo4Z/HnyT17t/SrhPuHX7qiJU9WIrUAqL6fs5EnZeD16lJIjRyhLTKpW6EB74QKFFy5QuGVLZX+enjhHRODcMQKny2m7nMLCUIhRmOuapl4C8UJ+KU/9Em9Io/zk8Pa1v1kSXB/otJB5BNJ2w5k9cCZOTrp/LTxbQcsoaNkNWkbKRqpnqyblYyqwDiYbrJ9//jmzZ89m+vTpjB07lvHjxxvmTZgwgeXLl9ukwWoK3i6VybCzC68ywgrQYTS0HwUn1su1hf+bD8NetYgcKqWKyR0nMzxoOO/vfZ91p9cBcCDzAJP+mcSdEXcyPWo6Ho627Z/Z2OnNFI6Ohtf/FUhaLeWnTlF67Jg8GnvsGKXHjqEvKDBaV5+fT/GePRTv2VPZn4MDTu3byyOxER1lYzY83Gp+sfaaLs4WsSddNqSsb646RkGZ7FozrmsATw5v32DbshXs6dg3CuVFslGauhPSdsojqOWFV1/HrQW0ikEX0I19GRLdx92Hg7d40LEE4vysjsk+rJ06dWLChAnMnTsXnU6Hg4MDe/fuJTo6mn/++YcHH3yQjIyMhpLXZMzx7biQX0rvt+XyrCM6+vPNlB5XXyH7JHzeB/Qa2WH88T1yzjcLs+PcDt7e/Tap+amGNh9nH56Kfuq6dROoD5IkoTlzhtKjxyg9dpSyY8cpPXYMbWbdggMcWrfGKSJcNmIjwnGKiBAuBU0cS/uwdujQoc6uSg3FthNZ3POt/IDWzNWBzc8OoZmbo1VkETQiZQXy6OnprXB6uxwkpb9KPIDaGQK7Q+se0KqH/F+MnFqVqq5KSUlJVvNhvXjxolnreXt7o1TW3W4x2WB1dnZm1apVDB8+vJrB+u+//zJ69GhKr1FjvjEx9QdGr9eTcSGTfh/vAyCqjTcrHut/7Q2tfx12fCJ/j5gAdyytj9i1Uq4rZ/GRxXyd8DWluko9R/pF8nLvl+ni16VBtmsuer2e7Oxs/Pz8TDoxrYk2J4fSY8cpO36M0mPHKT1+nPJTp+QqK9dA6e6OU3g4zuEdcAqPkA3Z9u0tli/WHvVpq5ijS2sFXTXUcS/T6hgzfxunsosAmPt/kdzes43F+rdVrsvrSFsG6bshZQuc2gpn918916lHIAT1hjZ9oE1PCIiU3eBq4LrUZwNiqj6tHXSlVCpNHqhRKBTs2bOH6OjoOq9j8pizl5cXF64ot1nB6dOnadGihald2hQ6nY64Pbtp5upCbrHm2i4BFQx6HhJ+hcILcPwfSN4MbYdZXD5HlZyfdULYBD7Y+wEbUjcAkJCdwJ2r7uSGtjfwZPSTtHC1jeOg0+nYuXMn48aNs5sbmdrXF/cB/XEfUPmgoi8poezEicsG7DHKjh6j9MQJpGLjqiv6wkJK9u2jZN++ykaFAsegIJzCw+UR2fBwnMLDcQgMRGGiTuxRn7aKPemyoWT96t8Ug7HaI7gZt8bYZ7EXU7GnY282kiQn5D+5UTZSU3fK1Rlrw68DBPeTM94E95WrQNWR60KfjYg96nPq1Kl1Lhal0+l46623TN6GyQbr8OHDmTt3LjfeeCPOznJuPoVCgVar5YsvvmD06NEmC2GL+Lo5GgzWauVZa8LZE0bMhhWPytNrXoRH/2uwKMhA90DmDZnHrvO7eHf3uyTnJQPwd/LfbEjdwINdH+TeTvfirG5a+ROthdLFBZfISFwiIw1tkl6PJi2N0uOJlCYel10KEhPRnj9vvLIkUZ6aSnlqKgXr11f26eaGU4cOOIV3kI3YDh1w6tChyeSMFdg2qTlFfBZ7EgCVUsEbN3VBqRSvd+2aklxIjoWTmyB5ExScr31Zv3AIGQChAyG4P7jbxiCHwD556KGH6NWrV52W1el0vPnmmyZvw2SDdc6cOfTs2ZNOnTpx8803o1Ao+Oyzzzhw4ABpaWn89ttvJgthi/i5O3Iyq4hSjZ6ich3uTnVQVeQk2Put7LienQT/fQRDGjYArU/LPvx+w+/8evxXPj/4OQXlBZRoS/j0wKf8kfQHT8U8xZiQMcKvsgFQKJU4hoTgGBKC55jKBzVdXh6liYmUVRiyiUmUnTiBVGY8Wq8vKqLkwAFKDhwwalcHtsS5fQeDAevUoQNOoSEoHIVfocAySJLEzL+PUK6V3Vym9g+hY8umm7+xySJJcinTE+sgaR2k7ar9Nb9nKwgbAmFDIXQQePg3qqiCpstHH31ESEhInZdXqVR89NFHBAXVfRQfzDBY27Vrx/bt23nmmWf4/PPPkSSJH374gaFDh7J06VKTBbA1FAoFHh4e+BVWjoxmF5TVzWBVKmHCR/DVYPmmse1D6HyznG+uAXFQOnB3p7uZEDaBBfEL+D3pd3SSjnNF53hh6wv8ePRHnuvxHNH+dfcVsRQV+ryeDGaVlxduvXrhVuVpU9JqKU9NlXPGJiZSlphEaeJxtOeqj4Boz52n8Nx5Cv/9t7JRrcYpNBSH9u0IUKspcnPDNSLCLLcCgYw9nZuWlnVHcg5bErMACPB05skRHSzSr71gT8e+Gnqd7It6fJX8yT1V83JqF3n0tN0I2Uj1a99gAVJ2rU8bxN70+eSTTzbKOvWqdFVWVkZOTg7NmjXDxcU260yb64w8e+URvtt+GoBlj/alR4gJVZE2zpJHVwGC+sJ9q2VjtpE4mXuS9/e+z45zO4zaRwSN4MnoJwnxCmk0WQRXR1dQQNmJE5QlJlJ6PJGypCTKkpLQFxXVaX2lqyuO7dvh1L49zu3by+m32rdH5ednNzc7e8JaQVeW5sHv49h4TM6G8fEdUdwY1arRZRCYgLZcDpQ6ugISV0NxTs3L+baD9qOh/Uj5t8dBuIRdD1j7ftJY1CvRl5OTE4GBTSvnml6vJz09HZ8qaV3qHHhVweAX4cgK+ck3bSfsXww9plpUzqvRrlk7vhr5FdvPbueDvR9w8pLsp7YxbSOx6bHc2uFWHu32KH4ufg0uS4U+27RpYzfO442JysMD1+hoXKtESkqShPbcOUqTkgzuBGVJSZSdOlWtDK2+uJjSgwmUHkwgr2q/3t6Xjdd2OLZrh3P79ji2a4faxivSNSb2dG5aUtbUnCI2HZeN1ZZezozv2tISItoVdnHsdRrZH/XoCjmQtzSv+jIKlRwoFT5WNlT9rlpvqsGwC33aEU1Fn3l5eTz11FOsXr0aSZIYO3YsH330ET5mlkWvk8E6Z86cOneoUCh4/fXXzRKmIalraVadTkd8fDy+/pXBNVlXq3ZVEw4uMHE+/HCjPL1hJnQYC56N+8PQv1V/erfszV8n/+Kz+M/ILslGJ+n4NfFX/k7+mymdp3Bf5/twc6h/Za7aqNBnYGCgXV94jYlCocChVSscWrXCY+hQQ7tUXk7RiRPsWbaMTu7uaJJTKEtKQnPmTLU+dJcuURwXR3FcnFG7ys8Pp3btKj/t5f+qJli96VqYcm5auzSrJa+j73ekGipa3dM3GLXq+rsubfa+pNdB6g44vAyO/iUHUV2Jgyu0Gy6nT2w/ClzN+/G3JDarTzulqejz8ccfJzMzk2+++YaCggLeeOMNpk+fzi+//GJWf3UyWGfNmmU0rVAouNKToOrrR1s0WE0pzQrg617pw5pj6ggryM7tUXdB/FIoy4c1z8OkJab3U0/USjX/1+H/GBs6lu+Pfs/iw4sp1hZToi3hy4Nf8uvxX3mw64NMipiEk0rUdbZlFI6OOHXoQEFUFL7jxuFwuVysvqiIsuTkypHYEycpO3ECbVZWtT502dkUZ2dTvGuXUbu6eXN5NLZthTHb9ro1ZGuiqZRmLSzT8vvedACc1Eru7GnfMQdNhsxjEP8THPq95sh+Rw8IHwOdboS2w8HRMnmdBQJLkJiYSHh49ViddevWcfLkSYPt5e3tzV133WX2dupksOqrJEw/ceIEY8eO5YEHHmDy5MkEBASQkZHB0qVLWbRoEWvWrDFbGFvCrz4uARWMelOO3CzOhmMr4fAf0OX/LCShabg6uDKt2zRu63AbXx38imVJy9BKWnLLcnl/7/v8eOxHpnWbxg1tb0CtFCXh7Amlm1u1lFsgj7KWnTwpf5JOyP9PnECXW33URpuVhTYri6IdO43aVc39cGrbDqe2bXFq1xbHtrIhqzbzlY7Auizbm24owXpz91aiopU1KcqBQ7/BwZ/h/MHq8x1cIXyc/JvRdpjwRxXYLN27d+ell17i5ZdfNgykAHh6epKcnEz37t0BSE5Oxtvb2+ztmGyZPPnkk9x77728/PLLhrbg4GBeeeUVNBoN//vf/+zaaFUoFDRv3pzmHlWzBJjoElCBqw+MfQ/+eECeXvWsnJS5kV0DquLn4serfV7lnk738Fn8Z6w5JR+rjKIMZu6YyXeHv+PRbo8yJmQMKqWq3tur0KcIALIMpuhT5e2Na48euPYwLi2szcmRR2FPnqQs+STll7/rLl2q1ocuK5virOojsipvbxzbtcUp7LIhG9YWp7ZhqAMC7OZY29O5aQlZ9XqJ73dWlnW+r3+IBSSzT6x27PV6OLUF9v8Ax/6Ry3lXRamGdiOh662yX6pjw7lrWRJ7upbsAXvTZ2xsLI888gi//PILCxcuZMCAAQA888wzDB48mNGjR1NUVMTGjRt5//33zd6OyVkCPDw8WLFiBcOHD682b+PGjdx8880UFBSYLZClMTd6rlSjI+L1tYBcAWbZtH7mCSBJsOx+OPKnPN1uJNz1u83UX068mMgnBz5h65mtRu3tvNsxPWo6w4OGo1TYrw+NoO5oc3IoO5lM2ckTlCcny9+Tk9Hl1BKRXANKV1ccw8JwahtmMGIdw9ri2KY1CoeayzraE/acJSD2eCb3L5Z9mvuG+fLzw30aZbsCoDAT9n8P+36AvLTq8wO7Q7c75dFUt4YPhhU0LWwhS4BOp+P999/nzTffZPLkycydOxdvb282bdrE2rVrkSSJ0aNHM3LkSLO3YfIIq5OTE3v37q3RYN27dy+OJiY3j4+P59VXX+XQoUNkZWXh4uJCeHg4jz32GHfffbfRsvv37+eFF15g165dqNVqhg0bxgcffEBYWJipu1ErOp2OEydO0L59ezyc1RSUas13CQDZMB0/T3akL7wAJzfAvsXQ436LyVwfwn3CWTB8Afsv7OfTA5+y98JeAE5eOskzW54hwieCR7s9yrA2w8x62quqT5Wq/iO21zsNqU+1ry9qX1/cehtXK9Hm5hoZsLJBm4I2M7NaH/riYkoPH6b08GHjGQ4OcnnasLBKgzY0DMfQUFTu1hlFsqdz0xKyfrfjtOH7/dfx6Co00rGXJDmRf9zXcPTv6qOprn4QdSdE3Q0tIhpGhkbCnq4le8Ae9alSqXjppZe4/fbbmTZtGh07dmT+/PlMmjSpRnvRHEw2WG+++WZmz56Nu7s7kydPplmzZuTm5rJ06VLmzJljskPtpUuXaNOmDXfeeSetWrWiqKiIpUuXcs8993D69Glee+01AI4fP86QIUOIiorit99+o7S0lBkzZjBw4EDi4+Np3ry5qbtSI3q9nsTERNq2bUtzd6fLBquZLgEVuPrADZ/BT7fJ0+tehbDB4GM5Q7u+RPtHs2j0InZn7OazA59xMEv2qTp+8ThPxT5FeLNwpnWbxtCgoSaNuFbVp71ceLaMNfSpbtYMdQ2uBbr8fMpTUmQjNjlFNmpTUuSsBVe+uNFoKE9Opjw5uXr/AQE4hYXiGBIqG7NhoTiGhja4e4EtnJt1zV5SX1lPZhayNUkOwmvj48Lwjtd3laMGPfaaUjnKf9cXcOGKBzcUcoR/9L1y5hh10/AhtoVrqSlRV31aO3tJTYSFhbFu3TqWLFnC//73P77//nu++OILgoOD6923yQbrvHnzSE5O5oknnuB///sfarUarVaLJEkMGjSIefPmmdTfkCFDGDJkiFHbhAkTOHXqFAsXLjQYrDNmzMDJyYl//vnHMOQdExND+/bt+eCDD3jvvfdM3ZVr4ufuREp2EYVlWko1Opwd6nEhdhgFMffJo6uaIvhzGty/GizgJ2opFAoFfVr2oXdAb7ad3caC+AUczTkKQGJuIk9tkQ3XhyMfZkTwCOEqcJ2j8vTEJSoKl6goo3Z9aSnlp09TlpxMecopylKSKU9Oofz0aaTy6g9/2owMtBkZ1QK+FK6uOIYE4xQahmNYKE6hsiHrGBKC0kYLlZiKqdlLzGX5/srUZ1P6hqBS2oZLUpOiMBPivpXLcxddkaHD1Vc2UmPug2Yh1pBO0ASxpewl27dvZ9WqVZSWltK/f3/uvvtuxo0bx3PPPUeXLl2YNWsWTz/9dL3SdJlssHp4eLB582bWrl3Lli1byMnJwdfXl6FDhzJq1CiLjYj4+fmRefmVo1ar5Z9//uHee+81urkHBwczdOhQ/vzzzwYxWH3dK59+swrKaONTz1Qio96ClC2QexrSd8mlWwe/UL8+GwCFQsGg1oMY2Gog285u4/P4zzmScwSQDddn/32WMK8wHop8iDEhY0RWAYERSmdnnCMicI4wfs0p6XRozp6lLCWF8uQUyk6lUJ5yivLkZHR51ZOiS8XFlB09RtnRY9XmqVu2NDJgHUNDcQoNQd2ypShVWwOxiZUGlKhqZWFykmHHp3JaKt0V7mOtekCvh6HzTaAWaQMFTZOvv/6aadOmMXz4cNzd3Vm4cCHr1q1j4cKFLFq0iHvvvZdHH32UJUuW8PXXX9Pjird1dcVsS2PMmDGMGTPG3NWrodfr0ev15Obm8vvvv7Nu3To+++wzQE6FUFJSQuQVaXsAIiMj2bBhA6WlpTg71z/th1KpJCgoCKVSiV/VXKxF5fU3WJ3c4aYvYfE4kPSw5R25SknIgHpK3TBcabh+Ef8Fh3PkV1wpeSm8vO1lvoj/gqldpjKx7UQcVdVfb1XVp6D+2LM+FSoVjkFBOAYFwRVvVbQXL1J+6pRszJ46LbsaVLgXVEmrZ1j+/Hm0589TtMO4/LDCyQnH4OAqhmwITiEhOIaEoLoinYo96bI+smbklXLsfD4Aka29jDKgXK9Y5NifPwj/zZcrUUlVzlGFSs6X2mc6tOlZX1HtAnu6luwBe9Pnu+++y4cffsiTTz4JyKOtgwYN4t1338XHx4chQ4aQkJDAm2++yeDBgymqY+nxK7GZobHp06fz1VdfAeDo6Mgnn3zCI488AkDO5Sjlmsp5+fj4IEkSubm5tGxZe7qoixcvotHITu9KpRJXV1fUarVRjlmlUolKpaJr167o9XqauVaqJ7ugzOD6UIFKpUKpVBr6rUCtltfTXlFGU61WQ1Af9ANfQLX1XZD0SMseQDFtO3oXHyM/FIVCYZCvpnadTlej7LW11yZ7XfdpYKuBDAgcwLb0bXxz5Bvis+IBSCtIY9bOWXwe/zl3d7ybW9regquDKw4ODoaHkC5duqDX65Ekyab26arHqYb2in2y9nHq0qWL4U1GU9knPDxw7t4d15gYo32SysvRnztHecopSpJPojl9mvLTp9GcOo2+hmwkUlmZXDwhKanaPJW3Nw4hITgEBeEQEoxjcDAdQkJQlJejc3SsUfbi4mJKSkoM7YWFhdX6bQxUKpUhl6GpxCZWBscNDW9hKZHsmvrok7P74d/3IGmtcbujB/S4D3o9At5t6i2jPVEvfQqqYW/6zM3NNSoc0KFDByRJ4tKlSwa7zdHRkTlz5jB58mSzt1MngzUyMpKffvqJLl261KlTvV5PVFQUv/76Kx07dqzTOq+88goPPvggmZmZrFy5kscff5yioiKee+45wzJXcze4litCaGio0fTMmTO58cYbSUurTDESHh5O+/bt2bRpE6WlpVy4oABkH9MLBaVs3brVKGVX3759adGiBevXrzcyBIYOHYqLiwurV6822ua4ceMoKSkhNj+Cfu6daF54FEVhBvz5KNkjP2Pnrt2GZT08PBg2bBjp6enEx8cb2ps3b06/fv04ceIEiYmJhvagoCC6d+9OQkJCtX2KiIhgz549ZFWpfBQVFUVwcLDJ+3Qp/hL/J/0f0W7RbCnbQoo2BYDMkkzm7Z/HFwe+oJ9zP2beMBNNvoadOyv9Em11n2o9TrGxhja1Ws348ePJzs62iX3q1asXLVu2bFL7dNXjFBJCbGkJtGoF/fuDJDGqTx8KExM5vH4DjtlZOGZl4Zidg2NuLlxhyINcSEEXH09plX0CULq7o//xB5KqGLkV+/T000+zcOHCan01NjqdjoSEBCIjI00Oaok9XsVgjRAGK5ipz7P7YMt7cGKdcbtbc+gzDXo8AC7eFpfVHqjP+Smojr3pc/To0TzzzDPk5+fj6urKZ599RocOHWrM4BQRYX5GjDrlYVUqlezevZuePev2ekOn0+Hg4MDevXuJjo42S7Bp06bxzTffcO7cOS5evEhERAQLFixg+vTpRss9//zzfPjhhxQXF9foElDhjHzq1Ck8PDwM+1PbCKter2f16tWMHDmS3an53P/9PgAeGRzG8yPbW27kriAD9TdDUBRnA6AfMQdd78p9s7UR1tr26WjuUb49/C2x6bFG7c4qZ25qdxOT2k7iyI4jjBw5EkdHR7vYJ1sejdRoNGzYsIExY8bg5OTUJPbpyvb67pNaoaAsPZ3SlBQ0qanyiGxqGprUVLQZGVyJU8eOBC/7vc4jrCEhIY2eh1Wj0bB69WrGVSnJWxfKtDq6z9lAcbkOXzdH4l4dgVIEXJmmz8xjsOkNSFxl3O7ZGgY8Bd3vBoemEQRoLuaen4KaMVWf1s7DeunSJZ555hlWr15NaWkp/fr1Y/78+XTo0MGi26mzS8BNN92Ek1PdfZ/qG3zVq1cvvvzyS1JSUoiJicHFxYVDhw5VW+7QoUO0a9fumv6rPj4+NR7IK59eKn60HBwcaOdfuXxqdrHhB/JKajuhrtru0wZu+QqWyKValZvnoAzuV83nSalU1ujHolKpanzyqq29Ntnru0/dWnTjk2GfcDL3JIsOL2L1qdXoJB2lulJ+SfyF35J+o7O6M8H5wUQFRNnFPl2tvTbZG3ufKrbVlPbpajKa2u4cGorzFW9VQM4TW56WRvnp05QkJ3Nqxw5CLqeWqkl2V1dXXF0rfdft7cc47lQuxeXyQ8fg8ObCWDWF3FTY8q5cOpUq4zqerWHQsxB1lwikEggAb29vFi1a1ODbqZPBOmXKFLM69/Mzv2JHbGwsSqWSsLAw1Go1EydOZPny5cydO9cwUpqWlkZsbCxPP/202du5Gi29nHFQKdDoJE7nmOckfFXajYABT8N/H4FeC7/dCw9vAQ/7y5HYrlk73h74No93f5wfj/7IHyf+oERbgl7Sc0hziHvW3UN0i2imdJ7CkDZDREosgVVQuroaMhi4aDTsat2amHHjrC1Wg7G5ijvAMOEOUDdKcuHf9+WE/7oqadg8AmHQc/KIqjBUBYJGp04G63fffddgAjz88MN4enrSq1cv/P39yc7O5vfff+fXX3/l+eefNxQEmD17Nj179mTChAm89NJLhsIBfn5+PPvssxaTR6lUEh4efvmVoJI2zVxJyS4iNacYSZIsn8h86KuQthvSdkDBOfj9PpjyN6jsaySngkD3QF7s9SKPdnuUX47/wtJjS8ktywVgf+Z+9mfuJ9gzmLs63sWNbW/E1aGemReuM6qen4L6YU+6NFfWLZcDrlRKBQPbW6a4SlOgRn3qNLD3Ozl7S8nFynZnbxj4LPR66Lp/9V8b9nQt2QOW1OfmzZtZsmQJO3bsID09HW9vb3r06MGMGTOIiYkxWtbcaqLz5s3jrrvuwt+/7oNt8+bN45577jGp6FOdfFgbku+++47vvvuOY8eOcenSJdzd3enWrRsPPvhgtdKs+/bt48UXX2Tnzp1Gymzbtm2t/dfXt+P+7/YYchjufmU4/p71T51VjcJM+GqwbLCCHGU6bq7lt2MFynRl/JP8D98f/Z5TeaeM5nk4enBr+1uZ3HEyAW4BVpJQIKg7lvIVawyfs9PZRQz5YAsAvUJ9+O2Rvg2ynSbByY2w9hXIrgwmRO0CfadDv/9dt8FUgsZDkiQyijI4nHOYw9mHSc1P5aMhH9VpkOxq95PbbruNnJwcbrvtNjp16kRWVhYffvghe/fuZd26dQwbNgyQq4n26tWLqKgoo0HB3Nzca1YTValU7Ny5k169etW6TFV0Oh2Ojo7ExcWZFOdk9bRW999/P/fff3+dlo2JiWHjxo0NKo9Wq2XPnj306tULtVpNsK8bIBusp7OLGsZgdW8Bk36E78bKr6D2fAWB3eU603aOSlLRMrsly8YvY9eFXXx/5Hv2ZOwBoKC8gO+OfMcPR39geNBwJnecTHSL6AYtx2nvXHl+CszHnnRpjqwinVXtGPQZ3gr1xlfh2ErjBSInwfAZ4NXaOgLaGfZ0LdkKeWV5HMk+QkJ2AoezZSM1pzTHaJn0vHSCvIPqtZ0FCxbQooXx9T9mzBjatWvH22+/bTBY61NNVJIkvv32W9auXVvrMlWpGuBqCuLMugJJksjKyjJEL4f4Vr6yTs0ppneYb8NsuHUPGPcBrPyfPP3PU9AiQjZc7ZgKfSqQixAMaj2IxIuJ/Hj0R1afWo1Gr0En6Vifup71qeuJ8IlgcsRkxoaOxVndAA8Hds6V56fAfOxJl+bIWrW6lfBfNUbSluF1eDGq2JWgKa6c0boXjHkXWsfUvrKgGvZ0LVkDjV7DidwTJGQlkJCVwKHsQ5zOP33N9Y5fPF5vg/VKYxXA3d2dTp06kZ6eDlimmujXX39dLznrgjBYr0Gwn5vhe4MEXlUlZgqcOwD7vgNtKfw8GR7aBJ6BDbvdRibcJ5w3B7zJUzFP8Vvib/yW+JvhyfL4xePM2DGDD/d9yM3tbub28Ntp43F9JeEWCOpLcbmWXSnyNRXo5UwHf3crS2RDpO9B/dfjdK76+t+tOYx6Ux5ZFW94BPUkuySbg1kH5U/mQY7kHKHsyrK9V+Dp6Eln38508etChHcEFw5eYGiboQ0iX15eHvv37zeMrta3mqi5I6amIgzWaxDiW2mwpuYUX2VJCzH2Pcg8Cum7ZZ/WnybB/Wvksq5NDD8XP6ZHTefBrg+yPnU9S48uNZR+zSvLY/GRxXx/5HsGtBrAHRF30D+wPyql7SdRFgjqSl0r8FVQNUft1fLZ7jiZQ7lW7mNwBz+0Wq3I0VtehMPWd5F2f4nicpoqSaFEipmKcvjr6J080VVZ3i72Cds4ThX7ULFuU9inK9tr26ey8jJS8lMMBmpCTgLpBelcDQelA+HNwunq15XIFpF09ulMK9dWBnc4rVbL+kPrjfazquz1rcD32GOPUVRUxKuvvgpYpppoYyAM1itQqVRERUUZcjK2buaCSqlAp2+g1FZXonaCSUvhm+FwKRUyEmD5QzBpCdihsXalPmvCUeXIhLAJTAibQEJWAj8f/5l1p9eh0WuQkNh2dhvbzm6jlXsrbu1wKze1uwk/F/NTptkzddGnoG7Ygi7rWoGvQ4cOeHh4sGHDBkP71SqGbUmq9F91z09l9erTTbcKWh32qXn+YaLSv8OhPIuK8dNc11AOtrkPvWskw1y8SU9Ntat9Ats7Trm5uQQEBDSpfbryOGkkDWd0Z1AFqTiSe4R95/dRSilXo5myGX2C+xDmHIYuXUeAKgC1To1HsQfDwoaRmprKmn/XGO1TVFQUp06dsngFvtdff52lS5fy6aefVssSUJ9qoo2BWVkCsrKymDdvHlu2bCEnJ4c///yTzp0789VXX9GrVy+bqoFbET3XoUMHVCoVjz32GI899phJfQx+P5bUnGLcndQcmjWqcQ5c5nH4dhSU5cnTfR+H0W81/HZthJySHP48+Se/Jf7G+aLzRvPUSjXDg4ZzW4fb6BnQU+R0FTQ4CxYsYMGCBeh0OpKSkiyWJaCuFfhMHRG66fMdJJzJQ6GAA68Ow81J3aRHuWrdp5J8lJtmodpXmdRcUjujH/wS+l6PglJtf/vUFI+TDe9Tsa6Yg9kH2Xt+L/sz93Pk4hG0emN5q+KkcqKTbye6+nYl0i+SSL9I/Fz8LLpP5lbgmz17NrNmzeKtt97ilVdeMbQnJiaaXU20MTHZYD116hT9+/cnLy+Pbt26sXv3bkNqgscee4zi4uIGzdtqKqamj9FqtWzdupVBgwYZLqZ7F+1ha5L8JLb3tRH4uTdS0ujkWLkSlnT5RB4/D3o+0DjbthA16dOk9fVatp3Zxm9Jv7H97HYkjE/XII8gbu1wKze0vQFflwYKiLMh6qtPQSXm6NJaaa1MkVWr09N55jrKtHpC/dyIfW6I2XLaNWf3wfKHIedkZVtwf7jhU7ReweI6siBN6b6UV5bH3gt72Zuxl30X9pGYm4heqt1H08fZh+4tutO9RXeiWkTRyacTDvXMo26qPutyP6kwVmfNmsXMmTOrbc/T05MpU6bwxRdfGM0bM2YMKSkpRiO91sLkM+uFF17A29ubvXv30qJFCxwdHQ3zBgwYUE0R9oYkSRQUFBg9YYX4urL18vfUnKLGM1jbDoUJ82Dlk/L06ufkFFgdJzbO9i1ATfo0BbVSzdCgoQwNGkp6QTrLkpax4uQKLpbKib3TCtKYt28enxz4hKFthnJL+1vo27Jvk/V1ra8+BZXYky5NkfV0ThFll/1XO7b0aGjRbA+dVq4e+O+7cgVBkHOqjnoDejwASiWSRmM3x94esKdr6UryyvLYd2EfcRlx7MnYQ1Lu1Q2zYM9gYvxj6N6iO9Etomnj0cbib10trc833niDWbNm8dprr9Voo1mrmqipmGywbtq0iS+++ILAwMBqDsEtW7bk3LlzFhPOVgiuEnh1OruYmODqjskNRsx9cDEFtn8Mkh6WPQD3LIeQAY0ng43QxqMNT8c8zeNRj7M5fTO/J/3O7vO7AXkkdkPqBjakbiDALYCb293Mje1upJV7KytLLRA0LkfPV/rfdQxomKIENkv+Ofkembajsi2wO9zyNfi1t55cApuhWFPM/sz97D6/m93nd3P84vFqb+4qUKCgfbP29PDvQYx/DNH+0XYXP/Hhhx8yY8YMxowZw/jx49m1a5fR/D59+gCNV020PphssJaWltYYSQZQVFTUJEuzGedibYTAqysZPgsKLkDCL6Arg5/vhPtXQ0DXxpfFBnBQOTA6ZDSjQ0aTlp/GHyf+4K+TfxlSY2UUZfDFwS/44uAX9G7Zm5vb3czwoOEir6vguuDY+XzD944tryOD9eQm2QWgOFueVihh4HMw+AW7LXUtqD8avYbD2YfZdW4Xu87vIiErAa1Usw+qAgURPhH0COhhMFK9nLwaWWLLsnKlXBRj7dq1NSb2rxjFjYiIYMuWLbz44ovceuutRtVETSmfWsHx48f5999/yc7O5oEHHiAgIIBz587RrFkzXFzMK3FsssEaHh7Oxo0bGTlyZLV5W7dupUuXLmYJYiuoVCr69u1rFDlsNMLaGKmtrkSphBs/k+tbn1gPZfmyb+vUdeATeu31rUhN+rQkQZ5B8qhr98fZemYrf574k21ntxl8jiqeoj0cPBgdOpob295It+bdbCLi0RwaWp/XE/akS1NkNTJYA68Dg1WnhS3vwLYPoWKkzLM13PotBPWpcRV7Ovb2gC3pU5IkTuefZse5Hew8t5O4jDiKtbX/boc3C6dnQE96BfQiJiAGT0frXzOW1OeWLVvqvKwlqonqdDoefvhhFi9ejCRJKBQKxo4dS0BAAI888gjdu3dnzpw5ZvVtssH60EMP8cwzzxAYGMhdd90FQHl5OcuWLePzzz/ns88+M0sQW0GpVFarDNHGxwWFAiTJSiOsII8Q3LYYfrgRzsRB4QX48WY5R6undXOjXY2a9NkQOCgdGB40nOFBw7lQdIGVKSv588SfpBXI6UwKNAUsS1rGsqRlBHsGc0PbG5gYNpGW7raru5poLH1eD9iTLk2R9eg52WD1dFYT6NXE3yoU5cCy++DU1sq29qPg5q/AtXbXLXs69vaAtfWZV5bH7vO72XFuBzvO7aiWWaYqIZ4h9G7Zm14BvegZ0JNmzs0aUdK6YW191oe33nqLn376iffff58xY8YYDWKOHTuWxYsXm22wmpXW6uGHH+abb75BqVSi1+tRKpVIksRDDz3El19+aZYgDYWp0bgajYb169czatQoHBwqXyP1f3czZy+V4OXiwMGZoxpS5KtTfBEWjYGKKi1+HeC+1eBu+pB9Y1CbPhsDSZLYn7mfP0/8yfrU9ZRoS6ot0zOgJxPDJjIyeCTujrZfnMGa+mxqmKNLa2UJqKusOYVlxLwpj5D0DvXh10f6mi2jzZNxCH6ZDJcu59hUqGD4DOj3P/mt1FUQ15FlaWx96iU9x3KOse3sNraf3U5CdkKtkfw+zj70admHPi370DewLwFuAQ0uX30xVZ+Wui9ZgrCwMB588EFeeeUVdDodDg4O7N27l+joaNasWcO9995rlP/WFMzKP7Fw4UKmTp3KqlWruHDhAn5+fkyYMIF+/fqZJYStcWVOOIAQP1fOXiohr0RDblE5zdwca1izEXD1gXv+hO/GyoUFspPkUdf7/rnqiII1qUmfjYFCoSDGP4YY/xhe6f0KG9M28tfJv9iTscewTFxGHHEZcby1+y2GtRnG+LDx9GvVDwel7f6IWUufTRF70mVdZD1WNeCqKfuvHvkTVkwHzeVXvW4t4PYfILjuBro9HXt7oKH1mVeWx85zO9l2dhv/nf3PkCnmShyVjkT7R9MvsB/9AvvRvll7u8zVba/n59mzZ+nbt+br0NnZ2agog6mYnTCtT58+huiy64FgXze2n5SDek7nFFnPYAXwagVTVsJ34yD/DGQekY3WKSvBxdt6ctkwrg6u3ND2Bm5oewNnC8+yKmUVK5NXcjr/NABlujLWnF7DmtNraObUjFEho5gQNsGu/V0F1x9V/Vc7NUWDVa+HLW/D1vcr2wKj5UqAXiIjSFNCkiRO5Z3i3zP/8u+Zf4nPjEcn6WpcNtQrlP6B/enfqj8x/jG4qM0L6hHUnxYtWpCSksLQoUOrzUtMTKR169Zm923fGX4bEeNMAcV0D7Ky30uzYJjyt2y0FmbIJVyX3CKPvjrbd1RjQ9PKvRUPRz7MQ10f4lD2IVYmr2Tt6bVcKrsEQG5ZLr8m/sqvib8S6BbI2NCxjA0dS4dmHYTxKrBpmnSGAG0ZrJgGh/+obOt2J0yYDw5N3Ff3OkGj17D/wn62pG9hS/oWzhSeqXE5F7ULfVr2YUCrAQxoNYBA98BGlVNQO+PGjeOtt95izJgxBATI7hcKhYK8vDw++eQTJk40P498nXxYQ0ND6/xDrVAoSE5ONlsgS2NqadaKhL0eHh5G+7z+SAYP/7gPgKdGtOepER0aXPY6kZUEi8dB0WWfkMDucPdym3EPqE2ftoZGp2HHuR2sSllFbHospbrqtaHberVlTOgYxoSMIcQrpPGFxH70aQ+YosuGKs1a137qKuuY+Vs5nlGASqngyOzRODtYP2rbIpTkwi93Qep2eVqhhFFvQp/pYMZ1IK4jy1IffRaWF/Lf2f/YnL6Z/878R4Gm5lfGwZ7BDGw1kEGtBxHjH4OjyopvORsYU/VpSz6sFy5coGfPnuTl5TF06FBWrlzJqFGjOHz4sMGftbbUqNeiTgbrfffdZ9JJaM+lWSVJQqvVolarjfY56UIBoz6SI1Fv7t6KjyZFNZTIpnPhCHw/EYpllwX8u8K9K8DN+gmOa9OnLVOkKWJT2iZWn1rNrnO7anwN1dGnI2NCxzAqeBStPcx/xWEq9qhPW8UcXVor6KouspZr9XSeuRaNTqJ9C3c2PDPYbPlsitxUWHpbZaCp2gVuXQQR48zuUlxHlsVUfWaXZLM5bTOb0zez+/xutPrq/ppqhZpo/2gGtx7MoNaDrDZIYA1M1actGawgG60zZ86sFuc0Z84cw6irOZiVJcCeMCcad/Xq1YwbN84oOq9UoyPidTnpbvcgb/6c3r/BZDaLzGOyH2vhBXm6eQTc+zd4+FtVrNr0aS/klOSwIXUDq0+t5kDmgRqX6ezbmdEhoxkVMqrBK2vZuz5tCXN0aWmDta5vfuoi69Fz+Yz7ZBsAN3QL5JM7u5stn82QcVhO31eUKU+7NYfJv0KrmHp1K64jy1IXfZ4tPMum1E1sTNtIfGZ8jdWlPBw8GNB6AEPbDKV/q/42kRPVGtT1/LT0mx9bR/iw1hFnBxUtvZw5n1dKqjWKB1yLFh3l9FbfT4SCc5B1XM4kMOVv8Gq8EcCmhq+LL3dE3MEdEXeQUZTButPrWHtqLYdzDhuWOZJzhCM5R5i3bx6dfTszMngko4JH0cazjRUlF9gDcXFxFvuBaXL+q+lxsPT/oDRPnvZtD3f9bvPFUgSVpOensz51PetT13M052iNy/i7+jMsaBjDgoYR4x9j0xlabI2KB92KB2BbYNiwYXz++edERERUm5eUlMSjjz7K5s2bzerbZIN169attc5TKpV4e3sTERGBWt30bOFgX1fO55VysaicvBINXi42dmH5tZNLtn5/A+SlwcVk+HaUHIjVPNza0tk9AW4BTOk8hSmdp5Cen8661HWsP72eYxePGZapMF7n759PR5+OjAgewYigEYR5h1lRcsH1gLHB6mFFSSzAqa3w0x2guVyopVUP2Vi1Ed98Qe3Udm+sSphXGMODhjMieAQdfToKt4wmxJYtW8jPz69xXkFBAf/++6/ZfZtsVQ4ZMuSaJ5e7uzvPPPMMM2fONFswWyTE141dKXLut7ScYrq2to0nGiN8QmWj9Ycb4GIK5J+FRaPhrmXQuoe1pWsytPFsw4NdH+TBrg+Slp/GutPr2JC6wegGfeziMY5dPManBz4l1CuUEUEjGB48nE4+ncQNWmBxjmU0kZRWiWvht3tBVyZPhw6CO34GJ9sv7HG9cq7oHNtKt7F07dJajdSOPh0ZGTyS4UHDxQP8dcr58+dxdXW99oK1YLIP66pVq3jiiScIDw/nzjvvxN/fn4yMDH7++WcSExOZM2cO27ZtY9GiRXz00Uc88cQTZgtnCSwZ3PDFlmTeW3scgE/v7M7EbjacSqMwS36ddv6gPO3gCpN+hHYjGlWM6y24IT0/nQ1pG9hweoOR20BVAtwCGNam8hWYWln358brTZ8NSVMKupIkiZg3N3KxqBxfN0f2vjbCPs+Po3/DsvuhIginwxi47XuLp60S11H9ySnJYX3qetacWnNV//5RIaMYGTySNh7CRaqu2FvQ1V9//cVff/0FwOLFixk/fjzNmxtX3ywpKWHLli20b9/+qm/qr4bJI6xr165l0KBBLF682Kh9ypQpTJkyhf/++4+FCxcC8PXXX1vdYDWHkpISPDyqv1IzzsVa1JgimY57c5jyj1y68PQ2uSLMT5PgxgXQ7Y5GFaU2fTZF2ni2YWqXqUztMpVzhefYlLaJjakbOZB5wBBkkFGUwU/Hf+Kn4z/h6ejJoNaDGBY0jP6B/XF1uPbT5/Wkz4bGnnR5NVkzC8q4WFQOyP6rdmmEJa4xNlY73wK3LARVw7he2dOxtxWKNcVsStvEqlOrrppBZXTIaEaHjG7UDCpNDXs6P48ePcrvv/8OyKlNN2/ejPKK8shOTk507dqVjz/+2OztmFyv7Ndff+XOO++scd7kyZNZvnw5ABMmTODEiRNmC2YttFotsbGxNZZFC23uZvieeKGwMcUyD2dP2RWg4w3ytF4Lfz4CW96DRkoOcTV9NnUC3QO5p9M9fD/2ezbfvpnX+7xO/1b9jUZU88vz+SflH57Z8gwDfxnI9I3T+S3xNy4UXaixz+tZn5bGnnR5LVmP2rv/6omNshtAhbHa7U74v28azFi1p2NvbbR6Lf+d/Y+Xtr3EkN+G8Mp/r7D97HYjYzXMK4zhzsNZMWEFv038jQe6PiCM1Xpgb+fnyy+/TEFBAQUFBUiSRGxsrGG64pOdnU1sbCyRkZFmb8fkEdaioiKysrJqnHfhwgWKi+UIeg8PjyYXeNW2uTtOaiVlWj2Hzlyytjh1w8EZblsMq5+Hvd/KbVvehtzTMPFjUDfd5Mu2hJ+LH7eH387t4bdTUF7Af2f/Y1PaJv47+x9FlwNLyvXlbDu7jW1nt/EGb9DRpyND2wxlUJtBdPTpaJf1sAWNg11nCEjZIr8J0skjxHS9TX4TpGwiRQ/slKTcJP4++TerTq0iuyS72vyqVQBD3UNZs2YNQZ5BVpBUYEvo9foG69tki7J///68/vrr9OzZk/Dwysjz48ePM2PGDAYMGABASkpKvWrG2iIOKiWdAj05kHaJ0znF5BVr8HK1sUwBNaFUwfgPoVkIbHhdbjv4E+Sly36tLlYuM3ud4eHoYbjRl+vKicuIIzY9ltj0WDKLMw3LVQRtfX7wc5q7NGdQ60H0b9mfMqnMitILbJFj5yurA9mVwZq6Q84GUBFg1elGuOlLYaxaidzSXFafWs1fJ/+qMXjK09GT0SGjmRA2gagWUYaHaI1G09iiCq5DTDZY58+fz6BBg+jcuTNdunTB39+fCxcucPjwYZo1a8b8+fMBOHfuHFOmTLG0vGbTs2fPOiXoBq46MtyttTcH0i4BcOhsHgPaW7+aVJ1QKKD//6BZMCx/GLSlsm/rNyPgzl/Ar32DbbqpjbRbEkeVI/1b9ad/q/682vtVjl48yr/p/7IlfYvRD0ZWSRZ/nPiDP078gQoVGzZvYHAbuQKMGNUwn7qem1UTdFuLq8maeDlDgINKQdvmdhJNf+GIbKxqS+Tp8PHwf9+CqnHuF+K+JKPVa9lxbgcrTq4gNj22WtUptVLN4NaDmdh2IgNbDay1JKrQp2WxZ30uWbKE+fPnc+zYMUpLq5c5N/c+alalq4yMDObNm8fWrVvJycnB19eXwYMH89RTT9GyZUuzBGkoLB0998e+Mzz7uxx5//zocB4b2q7efTY66XHw8x1QfPk1j5OnXOqw/UjryiUwIqMog61ntrL1zFZ2nd9Fma7mkdUgjyAGtBrAgFYD6BnQE2e1ZSOqBZVYK0vAteg8Yy1F5TpC/dyIfW5IvftrcC6lyzmiC87J0+1GwB0/gdrJunJdR5wpOMOfJ/9kxckVRm92Kujq15WJbScyNmQs3s7ejS+goM5YO0tAVf7++2/+7//+j/vuu49vv/2WqVOnUlpayt9//01gYCB33nmn2SlPRWnWK9Dr9WRnZ+Pn51ctyg3gZGYBI+bJKRnGdA7gy3vqVyLQauSehp8nQ+aRyw0KGDkb+v1PHo21ENfSp6BulGhLiMuIY0v6Fv5N+5fM0uo/MABOKid6BPSgf6A8ahvqGWqfEeONgDnnprUM1qvJWlCqoeus9QD0CfPhl4f7mi1Xo1B8ERaNgexEeTowGu77Bxzdrr6eBble70sanYZN6ZtYnrScned3Vpvv6+zLxLYTuandTbT1blvnfq9XfTYUpurTlgzWfv36MXDgQN5++20cHBzYu3cv0dHRZGRkMHDgQF5++WWmTp1qVt/izLoCnU7Hzp07ax2yDvNzx81R9q9KsJfAq5poFgIPrIeICZcbJNgwA5Y/BOWWS9l1LX0K6oaL2oVBrQfxco+XecLpCX4f9ztPxzxNjH8MKkWlv1+ZroztZ7czN24uN664kTF/jGH2ztlsTN1IfnnN1UeuV+zp3LyarBfyK1+5BXja+Oi6pgR+vrPSWPVpK1ewakRjFezr2FuC9IJ05u+bz4hlI3j+3+eNjFWVQsWQNkP4ZOgnbLhtA8/2eNYkYxWuP302NPasz8TEREaMqMwDXZHpICAggNdee4158+aZ3bdZThJLlizhp59+IjU1lZKSEqN5CoWC5ORkswWydZRKBV1aebH71EXO5ZWSVVBGcw87fY3l5A63/whb58KWd+S2Q79DxmE5GKsB/VoF5qNQKGjr3ZaI5hFM7TKVgvICdp3fxX9n/+O/s/8Zvd47V3SOZUnLWJa0DJVCRVe/rvQL7EffwL508etiUtECgW2SkVfpKhLg5WJFSa6BXi8/EKfvkqfdWsDdf4CbncQB2Bk6vY6tZ7bya9Kv7Di7w5AHuoLW7q35vw7/xw1tb6CFawsrSSloauh0OhwdHVEqlbi5uZGRkWGYFxQUREpKitl9m/xr9d577/Hyyy/TqVMnunXrhpOTnRpr9aBbG292n5JLtB46e4lhEf5WlqgeKJUw5CVo0QlWTIPyQsg6BguHwk0L5KhdgU3j4ejByOCRjAweiSRJnLx0ku1nt/Pfuf/Yf2E/Gr0cwauTdMRnxROfFc/nBz/H3cGdngE96RvYl74t+xLsGSzcBxoZU4JBayPDaITVhu/HW96BYyvl744ecPcyuZS0wKJcLL3I8hPL+S3xN84XnTeap1aoGRY0jFs73Ervlr1Fqjw7xxaCQa8kNDSUc+dk3/Ru3brx888/c8MNci74ZcuW1SvOyWSDdeHChTz22GN8+umnZm/UllEoFHh4eFz1hzuytZfh+8H0PPs2WCvodAO06Ai/3g1Zx6G8QE7k3ecxGDHL7HytddGnoO5cS58KhYL2zdrTvll77utyH8WaYvZe2MvOczvZcW4HKXmVT7eFmkJDOi0Af1d/+rTsQ++WvenTsg/NXZvXuI2mgi2cm3FxcXXyObuarEYuAV426hJweLn8JgdAoZRzQ7fsZjVxbOHYW5ojOUf46dhPrDm1xvCQWkEr91bc2uFWbmp3E34ulh/Rbor6tCZ11WfFg26FD6stMHz4cDZu3Midd97Jk08+yaRJk4iLi8PR0ZHExETeffdds/s2OejKzc2NlStXMmzYMLM32pg0hDNyWk4xg96Xf+SHRbRg0X09LdKvTVBWCP88JbsGVNAqRk43I0ZD7J7zhefZdX4XO8/tZNf5XeSW5da6bJhXGL1b9qZ3QG96BPTAy8k2bojWxBazBLy24hBLdqUB8Of0fnQPsrG8yufi5SCrivRVo9+GvuaNJguM0eg1bErdxNJjS4nPijeap0BB/1b9uTPiTvoH9kclcts2WWwp6KqsrIyysjKDHMuXL2fp0qUoFAomTJjAfffdZ3bfJo+wxsTEkJycbDcGq6no9XrS09Np06ZNrdF5bXxc8HZ14FKxhoQzl5Akqek8WTq5wy1fQ5vesPZl0Gvg7D74ahBMnA9d/s+k7uqiT0Hdqa8+W7q35Ob2N3Nz+5vRS3qOXzzOrvO72H1+N/su7DNKnZWSl0JKXgo/H/8ZBQoifCLoFdCLXi17Ed0iGndHO8n3WQv2dG5eTVZjH1YbG2EtzJSrWFUYq1F3QZ/p1pUJ+zr2NZFXlscfJ/5g6bGl1VJSeTp6ckv7W7i9w+208WzTKPLYuz5tDXvWp5OTk5Gr6C233MItt9xikb5NNljnzZvH3XffTXR0NDExdprS6SrodDri4+MJDAys9URRKBR0beXFthPZZBeWcy6vlFbeNhzsYCoKBfR6CFr3gGVT4WIKlOXL31O2wJh36xzVWxd9CuqOJfWpVCjp5NuJTr6dmNplKmW6Mg5mHpQN2IzdHMk+YqgXLiEZKm99f/R7VAoVnXw70SOgBz39exLtH42bQ+NGetcXezo3ryZrhUuAUgHN3W3Ih1VbLrsY5Z+Vp1v3ggkfWTRtnrnY07GvSlp+GkuOLWHFyRWUaI0Dntt5t2Nyx8lMCJuAi7pxf4/sVZ+2itBnzZhssN5///3k5OTQq1cvAgIC8PX1NZqvUCg4ePCgxQS0Vbq19mbbCTnx/qEzl5qWwVpBYHd4ZCv88wwc+k1u2/8DnN4OtyyUDVpBk8FJ5USvlvIIKkBheSH7M/ez6/wu4jLiOH7xuGFZnaTjUPYhDmUf4rvD36FSqOjo05EeAT3o4d+D7v7d8XS0oxKhdkxF0FVzDyfUKhv6cds0G9J3y989W8GkJaIwgJkkZCWw+MhiNqZuNIr2V6BgcOvB3N3pbnoF9Go6b/oEds2KFStYunQpqamp1Spd1cdGNNlg9fX1xc/P/tKQWCIatypGgVdn8hjTxbYqfFkMJw/ZOG07FFY9C5piuJgsV6kZ9DwMeg5UDtaWUtAAuDu6M6j1IAa1HgTApdJL7L2wl93nd7P3wl5OXjppWFYn6Ticc5jDOYdZfGQxSoWS8GbhxPjHEOMfQ/cW3fF18a1tUzaPLUbjAmh0erILZZcAm8rBenwV7PxM/q5ylNPkeTSB4NRGRC/p2XpmK98d/o79mfuN5rmoXbih7Q3c3fFuQrxCrCOgQFAD77//Pi+++CLNmzenXbt2uLlZ7s2byQbrli1bLLbxxsSUaNzmzZtf80k1srW34btdFxCoCwoFRE2WX+n9+bDs0yrp4N934eQGuPmrWnO21lWfgrphTX16O3szIngEI4JHAJBTksO+C/vYk7GHfRf2GRmweklvcCFYcmwJACGeIcT4xxDtH010i2haubey6nlhii6tHY1bm6xZBWVUhM3624rBmntaTpFXwai35MBNG8KW70tavZY1p9aw6PAio2sKwM/Fj7s63sVtHW6zqSBIW9anPWLP+vz888+ZOnUqX331FSqVZQP9RGnWetDrrY1kFpTh4awmYeYouzy5TEanhW0fwr/vyUYrgMoJhr0KfR8HEYl63XKx9CL7Luxjb8Ze9l7Yy4ncE9WSlVelhWsLoltEE9UiiugW0XRo1sHmI5ltLUvA/rRcbvl8BwD39g1mzo1dzO7LImjLYdFoOHd5RLDTjXDb9zbht2rrlOnK+PPEnyw+spizhWeN5oV6hXJf5/uYEDYBR5V5KQYFTRdbyhLg6enJihUrGiQw3+wyN3l5eSQlJVWrdAUwaNCgegllTXQ6HSdOnKB9+/bXfDqIbO3NxmMXKCjVcjqnmFA/+wo6MQuVGoa8CO1GyKOtOSdBVyaXdT36N9z0OTQPNyxuij4F18aW9enj7GMoYACQX55PfGY8ey/sZd+FfRzNPopW0hqWzyzOZO3ptaw9vRYAV7Urkc0j6d6iO1HNo4hsHtmgmQhsWZdXUpusGXmV/mE2McK6YUalsdosFG741CaNVVs69sWaYn5P+p3FRxaTXZJtNK9b82480OUBBrcZbNNJ/m1Jn00Be9Zn//79OXbsmG0YrFqtlkcffZQffvihVn8uW/PzMgW9Xk9iYiJt27atg8HqxcZjFwDZLeC6MFgraB0Dj2yD2Ldg5wJAgrN74cuBMPgF6P8kqBxM0qfg2tiTPj0dPY18YEu0JRzKOsS+zH3sv7CfhKwEirXFhuWLtcXsOr+LXefl0p0K5CIIUc2jiGoRRbfm3Wjj0cZibzLsSZe1yVrVYLW6D2vSOtj9hfxd5SgXB3C2ndfWVbGFY1+kKeKX47/ww9EfuFh60Whe/8D+PND1AXr497CLN3e2oM+mhD3rc/78+dx88820adOGMWPG4OhouTcCJhusH330EStXrmTRokXce++9LFiwAAcHB77++mvy8vL45JNPTOpv8+bNLFmyhB07dpCeno63tzc9evRgxowZ1dJm7d+/nxdeeIFdu3ahVqsZNmwYH3zwAWFhYabuhkWoGniVcCaPG6NaWUUOq+HoCqPfgo43wF/TK0dbN78Bh/+AiZ9AQJS1pRTYCC5qF6MsBFq9lqTcJA5kHmD/hf3EZ8aTWVKZU1JCIik3iaTcJH5LkrNU+Dj7ENk8km7Nu9GteTc6+3bG1cHVKvtjC9hMlauiHPjr8crp0W9DYJTVxLFlijXF/HT8J74/8j2Xyi4ZzRseNJyHIh+is29n6wgnENSTdu3aMWLECG6++WYUCgWursb3Z4VCQV5enll9m2yw/vjjj7z66qvceeed3HvvvfTu3Zvo6GgefPBBRo8eTWxsLKNGjapzf1988QU5OTk8+eSTdOrUiaysLD788EP69OnDunXrDMPKx48fZ8iQIURFRfHbb79RWlrKjBkzGDhwIPHx8TRv3vhlJKsGXh1Iq71iUJMnqDc8+l/laKukh8yj8O1IlDFTUet6WVtCgQ2iVqoNeWDv6ngXkiRxrugcBzIPEJ8Zz8GsgyTlJqGX9IZ1LpZeZEv6FrakbwFApVDRvll7Iv0iiWweSdfmXQnxDLHp16eWJMMWDFZJglVPQ9Hlh432o6Hng9aRxYYp1hTzS+IvLD682KjCnAIFY0LG8FDkQ7RvVnPwqkBgL7zwwgt89tlnREVF0bFjR+uOsKakpNCtWzdDMtuqObYeffRRnnzySd55550697dgwQJatGhh1DZmzBjatWvH22+/bTBYZ8yYgZOTE//884/BqTgmJob27dvzwQcf8N5775m6KzWiVCoJCgqqU7JeHzdH2rVw52RmIfHpl8gr1uDlep2meHJwgVFvypWwVj4J5w8CEqp93zLKaQWqtkqIvM0m/dnsCVPOT3tDoVDQyr0VrdxbMSFsAiC/Nj2UfYiDmQeJz4onISuB/PJ8wzo6Scfxi8c5fvG4YRTWw9GDLr5d6Nq8K5F+kXTx61JjSi170mVtstqES8ChZXD0L/m7i4/N+q1WpTGPfZmujN8Tf+frQ18bvfpXKpSMCx3Hw5EPE+pl32Wv7elasgfsWZ+LFy/mxRdfNMkOrCsmG6xubm6Ul5ejUCjw8fEhNTWVfv36AeDi4kJOTo5J/V1prAK4u7vTqVMn0tPTAdlv9p9//uHee+81ioALDg5m6NCh/PnnnxYzWFUqFd27d6/z8oM7NOdkZiF6CbYnZzOuaxPNx1pXArvDg5th95fyiKumGIeyHPjzITjwA4z7AFpEWFtKu8XU89PecXNwo0/LPvRp2QeQ02Wdzj/NwcyDHMw6yKHsQ5y8dNJoFLagvICd53ey8/xOQ1ugWyBd/LrQ1a8rnf0608m3E24Obnajy9qOe4VLgIeTGjcns2NozSfvLKx+tnJ6wjy7yLfaGNeRRq/h75N/82XCl2QUZRjaFSgYGzqWR7o9QpiXddzZLM31dl9qaOxZnzqdjpEjRzZI3ybf4SIiIjh16hQA/fr1Y968eQwcOBBHR0fmzp1LeHj4NXq4Nnl5eezfv98wupqcnExJSQmRkZHVlo2MjGTDhg2Ulpbi7Fz/EQadTkdCQgKRkZF1cnYe3KE53/4n6+PfxCxhsIKcSaDf49BxItKaF1EkrZHbT2+DL/tDn2kw6AVwFpWQTMXU87OpoVQoCfMKI8wrjJvb3wzIo7BHso8YDNhD2YeqRVufKzrHuaJzrE9dD8hGQ5hXGIHKQAa0G8DtEbejVja+wVfXgiY1HXdJkgwuAf7WcAeQJPjrMSi97I/W9TbofHPjy2EGDXkdSZLE+tT1fHrgU1LzU43mjQoexfSo6bT1bmvRbVqb6/2+ZGnqqk9bLGgyatQodu3aZRtZAiZNmkRSUhIAs2fPZtCgQQQHBwPg4ODA8uXL6y3UY489RlFREa+++iqAYdTWx8en2rI+Pj5IkkRubi4tW9ZuLF68eBGNRgPIw+2urq6o1Wr0+sqRGaVSiV6vJy0tjfDwcBwc5Nf7KpUKpVKJVqulatpalUpFr1AfnB2UlGr0bEnKpLy83LCeVqulKmq1usZ2Bwc5mr7qSadQKAzy1dSu0+mqya5SqWptr0n22vZJqVQadHUt2a+6T15tKL1pEQnL5tIr5w+UeWmg18KOT5EO/opuyCtIkXeiUKntZ5+sfJw0Gg1paWl06tQJlUrVJPbpynZT98kRR7r7dae7X3dD+7nCcxzMPMiRnCMczjnM0ZyjlOoqX59LSCTnJZNMMocSDnF7h9vR6Cq3WyF7cXGxUeq+wsJCLEldC5pU3Je6dOli+AHLL9FSqpH1bRV3gH2LISVW/u7REsa93/gymElN+rQEe87v4aN9H3E457BR+6DWg3g86nE6+na02LZsiYbS5/VKXfVp7YImNfH6668zadIk3NzcGD9+fK12mzmYbLBOnz7d8L179+4cPXqUFStWoFAoGDlyZL1HWF9//XWWLl3Kp59+Wi1LwNXSe1wr9UdoqLGP0MyZM7nxxhtJS0sztIWHh9O2rfzku2HDBkN7VFQUwcHBbN26lYKCAkN73759adGiBWHueo7mwoX8Mr79Yw13jR+Ki4sLq1evNtrmuHHjKCkpITY21tCmVqsZP3482dnZ7NxZ+QrTw8ODYcOGkZ6eTnx8vKG9efPm9OvXjxMnTpCYmGhoDwoKonv37iQkJFTbp4iICPbs2UNWVlad92n9+vVGBsLQofXYJ6/urPLoTJfcDYSeXQG6MhRFmahXPcWlzfNJ7/IYXcc/bF/7ZOXjlJOTQ2BgYJPaJ0seJ4cSBzSHNXS4/Hd3y7sJ6h7E1qSt7Dy1k7O6s2ToMtCjp4tvF06ePFnjPj399NMsXLgQW+R8fqUh3eg5WAsuwIaZldM3fgYuzRpXBhsiKTeJefvmsf3sdqP2GP8Ynop+iqgWUdYRTCBoZLp16wbAM888wzPPPFPjMuaOCNtUpavZs2cza9Ys3nrrLV555RVDe2JiIhERESxYsMDIYAZ4/vnn+fDDDykuLq7RJaDiyePUqVN4eHgA1x5hXb16NSNHjqzTCKtSqeSbrSd5c7X8Y/fi6A48OqQdcH2NctW2T6WlpWzYsIGRI0fi6OiIuuAM0vrXURz722h5Ok5EN3QG+maVDxa2uk/WHmHdsGEDY8aMwcnJqUns05XtjbVPhaWFLF2/lH79+tGtRbcaZa9phDUkJKTRK11pNBpWr17NuHHjDPelLYmZ3PddHACPDW3L86Mb0Td82QNweJn8vdtkuPmLxtu2BahJn+aQVZzFgvgF/HnyTyM/6vbN2vNU9FMMbDXQLvKo1hdL6VMgY6o+banS1axZs655zs+cOfOq82vDCl76NVNhrM6aNcvIWAVo27YtLi4uHDp0qNp6hw4dol27dtf0X/Xx8anxQNY03B4eHo6Tk1O1eRU/kFcyrGOAwWD9LzmHaUPl1CS1nWg1tSuVyhojAmtrV6lUNcpeW3ttstfWborstbUrlUqcnJyM9dksBMWkH+HUNlj7Ely4/Ors2EpUiWtQ9ZgKg18ENz+b3SdrHielUkl4eLhhuinsU11kNLW9LvvkqfRkROcRtG/RvlbZXV1djfIIWuvHuOK4V90n4xysLo0nzMmNlcaqSzM5O4idUZM+TaFEW8LiI4v57vB3lGgrH2haurXkie5PMC50nM2XGbYk9dWnwBh71uesWbMarG+b0MYbb7zBrFmzeO2112q0vNVqNRMnTmT58uVGrxDT0tKIjY3llltusZgsKpWKiIgIk/xwQv3caOMj/2DEncqlqEx7jTWuH2rVZ+hAeGQrTPwY3C9HFeu1sGchfBwF/86FsoJq/V3vmHN+CmrGnnRZk6wZeWWG743mw1peDP9Uec036k1wq54yzNYx99hLksTqlNVM/HMin8d/bjBW3R3ceTrmaVbevJKJbSdeV8Yq2Ne1ZA8IfdaM1UdYP/zwQ2bMmMGYMWMYP348u3btMprfp4+czmb27Nn07NmTCRMm8NJLLxkKB/j5+fHss8/W1LVZaLVa9uzZQ69evWodAboShULB4A7NWbIrjXKdnl0pOQzvaPupXRqDq+pTqYKY+6DLrbDzM9j+CWiKoLxATom1+ysY9Bz0mApqJ6vIb2uYc34KasaedFmTrEZFAxrLYN36Ply6HPke3B+i7mqc7VoYc4794ezDvLfnPeKz4g1taoWa28JvY1q3aTRzvn59eO3pWrIH7E2fP/zwA+PHj8fX15cffvjhmsvfe++9Zm3H6ppYuXIlAGvXrmXt2rXV5lf4uEVERLBlyxZefPFFbr31VqPSrJasciVJEllZWZjq2ju4QwuW7JIDTv5NyhIG62XqpE8ndxjykmy8bnkH9v8Ikg6Ks2W3gZ0LYNDzEDUZVNe3f5S556egOvaky5pkreoS4O/VCA90F47Cjsult5UOMGG+zRcIqA1Tjv3F0ovM3zefP0/+adQ+sNVAnuv5XJPJpVof7OlasgfsTZ/33Xcfu3btwtfXl/vuu++qyyoUCvs1WLds2VLnZWNiYti4cWPDCVMP+rb1xUGlQKOT+Dcp69orCKrjESC7CPR9Ara8DYf/kNvz0mHl/2DbhzD4BYi8Q871KhBcx1RUuVIrFfi5NbDBKkmw5gXZbQdgwNPQvEPDbtPKaPVafkv8jc/iP6OgvNI9KdQrlBd6vsCAVgOsKJ1AYDucOnXKkFa0Ik9/QyB+9S2Eu5OaHsE+7EzJITWnmNPZRYT4uVlbLPvErx3cugj6PwWb34ATcrJ3LqXKicq3fiCPuEbeft2PuAquXypGWFt4OKFUNvBI5/FVcuEPgGYhMNBybli2SHxmPG/uepPE3MpUZ+4O7kyPms4dEXfgoBT3HYGggopc/Fd+tzRmBV0tWbKEcePG0blzZ8LCwow+FXlM7RWVSkVUVJRZzs6DwytdE8Qoq0x99EnLSLjrd3hgA4QNrWzPPQV/TYdPoiHuW9CU1t5HE6Ne+hQYYU+6vFLWMq2OnKJyoBGqXGnLYcPrldMj3wAHKxQqsCC1HftLpZeYtWMW96y5x8hYvbHtjay8eSX3dLpHGKs1YE/Xkj1gz/oMCwvj4MGDNc47fPgwYWHmu9CYPML63nvv8fLLL9OpUye6deuGk1PTCoZRKpVmPyEM7tCcd9ccB2SDdUq/EAtKZp/UR58G2vSCe1dA2i6IfRtO/Su356XBqmfkQJB+T0D0FNkftgljEX0KAPvS5ZWyZuZXZgho2dAG656FcDFF/h48ADpObNjtNQJX6lOSJP5K/ot5e+eRW5ZraO/o05FXer8iEv9fA3u6luwBe9bn6dOnKSsrq3FeaWkpqampNc6rCyYbrAsXLuSxxx7j008/NXujtoxWq2Xr1q0MGjTI5Oi8iAAPWng4kVlQxvaT2RSUavBwvr6fxuujz2oE9YEpf0PqTtj2gZwPEqDgPKx7RU6F1esh6P2oUR7XpoRF9XmdY0+6vFLWqhkCGrTKVVGOfF0BoIDRb9ltoFVVquozvSid2Ttns+/CPsN8dwd3nuj+BJPCJ113KarMwZ6uJXvA3vVZW+GAlJQUQwEnczBZExkZGdx8881mb9Ba9OzZE5VKZai9WxuSJFFQUGBWdJ5CoWB05wB+3JVKmVbPmsMZ3N6jTX3Etnvqo89aCe4LwX/AuQOyP+vxf+T20kvyaOuOT6H73dBnOvjat4vKlTSIPq9TTNHlggULWLBggdklBevLlbJWBFxBA6e0+vddKMuTv0dNhsCohttWIyJJErn5uXx96Gu+OfwN5fpyw7yxIWN5vufzNHe1XPaZpo64L1kWe9Pn999/z/fff2+YnjZtWrVCTSUlJRw8eJDBgwebvR2TfVhjYmJITk42e4PWIi4ujqNHj17VWLUEt0S3Mnxfvv9Mg27ruiewO9yxFKbvkstDKi8/f2lLIe4b+DQGfrkLUnfIUc4CgZk89thjHD16lLi4OIv227NnTzp16sSCBQtMWs+4ylUDGaxZibKPOICDKwx7/erL2xGHsw/zRcEXfJ7wucFYbe3emq9GfMXcwXOFsSqwCxYsWECnTp3o2bOnVeUoLi4mKyuLrKwsFAoFly5dMkxXfDQaDZMmTeKrr74yezsmj7DOmzePu+++m+joaGJiYszecFMlqo03YX5upGQXsSvlImdyi2ndzPXaKwrMp0VHuZb5sFdh5+ewb7FcgABJHn09/o9s3PaeBp1vEkUIBDZDXFycWbW/q46wNphLwIYZcj5kkDN2eLZsmO00IqXaUhbEL+CHoz+gl/QAKBVKpnSawrSoabioG7HErUBQTyreGOfn5+Pl5WU1OaZNm8a0adMACA0N5Y8//qBbt24W345CMnHMuWvXrmRkZHDx4kUCAgLw9TUuy6dQKGqNELMGFQcyLy+vTj8Mer2e7Oxs/Pz8zK7j+9nmE3ywPgmA50Z14PFh7c3qpylgCX2aTEmubLTu/kr2b62KWwu5claPqeBhf8UdrKLPJoo5ujT1fmKpfq6U9fGf9vNPgnxub3luiOVT6KXvgW9Hyt89AuGJfeBo3w/e8ZnxvL79dU7nnza0RfhEMLvfbDr5drKeYE0AcV+yLKbq01L3JVvH5BFWX19f/PyaZkALyNF5LVq0qFcfN3VvZTBYl+8/y2ND29XqhNzUsYQ+TcalmZzYvM9jcHSF7NOakSDPK8qU/fK2fQidboCeD0JQX7sJJLGKPpso9qTLK2VtcJeA2Lcqvw950a6N1VJtKZ8c+IQlR5cgIY/POCgdmB41nfs634daaX9BLbaGPV1L9oBCkvAuKbFL4z8hIYFLly4xaNAgAAoLC3nhhRfYv38/o0aNYvbs2WbbQyZfqaZUprJHNBoN69evZ9SoUTg4mBfh37qZK33CfNiVcpGU7CLi0y/RPej6rDNtCX2ajdpRLi7Q9TY5JdbuL+HYSvk1p14jV9I6/Ae06CSPuEZOAmfbfjq1qj6bGPakyytlPX/ZJcDb1QFnBwtHsZ/eDilb5O/ewRB1l2X7b0QOZx/m5W0vG42qRvpFMqP3DE7uPokUIZmZjVxQFXu6lmwRTUYGJQkJlCYkUHIwgZLDh9GXl9Nuz24c3eyrANEzzzxDdHS0wWB99dVX+frrr+natSvvvPMOzZs354knnjCrb/FoWQNarbbefdwS3ZpdKRcBeZT1ejVYwTL6rBcKxeXMAn0h74wcSLL/eyjOkednHoXVz8GGmdD1/yDmPgiMttlRV6vrswlhT7qskFWSJEMeVotnCJAkOddxBYNftMtqchq9hoUJC/k64Wt0l/1wHZWOPNH9Ce7pdA96nZ7j2uNWlrJpYU/XkjXRl5RQeuQIJQcPUhJ/kJKEBLQXLlRbTgGUHT+Oo53FCh0+fJjHH38ckO9VS5cuZfbs2bzyyiu89tprLFq0qHEN1osXL/LRRx+xadMmcnJy8PPzY8SIETz11FM0a3b9GmZVGdslgBl/HaZUo2dlwjlem9ARJ7XI52d1vFrDiJkw5CU4+pecTSB9tzxPUwT7f5A/AV3lQgRdb5VdDAQCGyG/REu5Tg4Yau5h4QDCU/9C6n/yd9928lsHOyPlUgov//cyR3OOGto6+3bm7YFvE+YlV9nRX9afQNCQSJKEJjVVNk4vG6iliYlwjfR46pYtyfXzAzvMwXrp0iWD2+jBgwfJzc3l9ttvB2D48OH1yuFvsjbOnj1L//79SUtLo2PHjgQFBXHu3DneeOMNfvjhB7Zv305gYKDZAjUVPJwdGN05gL/iz3GpWEPs8UzGdLH/KNsmg9pJdheIvB0yDsHeRZDwO5QXyPMzDsmjrutelSv7dL8bQgeDHfoUCZoWBWUaw3dPFwuOflYbXX0JVPbzgylJEr8n/c77ce9TqpNdJtQKNQ93e5gHuz4oSqoKGhx9URElhw5REh8vj54ePIguN/eq6yhdXXHu2hWXbt1w6RaJc9eu0KwZR1evpnsn+wsG9PX1JT09HYDY2Fj8/f1p164dAOXl5fXKLWvy3eiVV16hpKSE3bt3G+X+iouLY+LEibzyyissXrzYbIGsjVqtZujQoRapLnFLdGv+ij8HwB/7z16XBqsl9dlgBHSFCR/JNdKP/ClnGDi7V56nK4PDy+SPVxB0mwTd7rRaQQK70KedYE+6rCprYVmJod3DyYKyn9xU+baheQR0ucVyfTcwF0svMnP7TLac2WJoC/MK4+2Bb9PZt3O15e3p2NsD16M+JUlCc+YMJQcOUBIfT/GBeMoSE0F/ldF7hQKndu1w7haJa1QUzpGROLVti0Jl/PZVkiS71efAgQOZNWsW2dnZfPTRR4wfP94w78SJE7RpY34xJZO1sXbtWt58881qiWp79uzJnDlzeP11+08u7eJimVx8/dv6Gkq1xh7PJKugzPKv8OwAS+mzwXFyh+h75E/GYYhfCgd/gRLZF5m8NLmS1tb3oXUv6HYHdL4ZXH0aVUy70acdYE+6rJC1oLTSV9DdUgarJF2RGeAlsJOSpDvO7eDV/14luyTb0HZnxJ08E/MMzurafXzt6djbA01dn/ryctn39EA8JQcOUBx/AF1W9lXXUXl7yyOnUd1wiYrCuWtXVO7uddqeJfVZUFDAG2+8QXx8PAcOHCA7O5uZM2cya9asasvu37+fF154gV27dqFWqxk2bBgffPABYWFhddrWO++8w9ixY3nyySdp27YtM2bMMMz7/fff6dOnj9n7YfL7zby8PEJCQmqcFxoaSl5entnCNCR1rSij1WpZvXq1RRzI1SolN1+ufKXVSyzecarefdobltRnoxLQBca8A88eh9t/gHYjQVHlcjmzB1Y9Ax90gJ8nyyOzmpLa+7MQdqtPG8QUXVq7okxVWQurGqzOFjJYU7bAuf3yd/8u0PFGy/TbgGj0Gubvm88jGx4xGKs+zj4sGL6AV3q/clVjVVxHlqUp6lObm0vB5lgyP/yQ03fdTVKPnqTeOZnMuXMp2LChurGqUOAUHo73pEm0fOcdwtaspv3OHbT56kv8pk3DrW/fOhurltZnTk4OCxcupKysjJtuuqnW5Y4fP86QIUMoLy/nt99+Y9GiRSQlJTFw4ECysrLqtK3Q0FCOHz9OdnZ2tRHVzz77jHfffdfs/TD5bhcaGsqqVasYOXJktXlr1qwhNDTUbGEaEnMrytSX+/qFsOi/U2h00v+3d95hUR1dA/9toyysFEEFEVTsUWNDxV5jj2800ajJmxjT1PTET5O89jSN6d0UTWI0VU1UNKJiL1hixNiJiIqKdJa+u/f748rCSnEXFtjF+T3Pfbg7M8w9c+7u7Nm5Z87h+30XmNqvmf1WRQRVj9oV2oyWj8yrcOxn+HuVHFkA5PBYpzfIh2sdaDUS2o6Fpn2dcne1oHQcJaMMQGZe0ZeYzs1O77G9Hxad937R4X21r+iv8H87/4+j14+ay3o27MlrPV/Dz732xgkXVA2Fj/ezDx8m5/ARso8cIf8WKeiVnp64d+iAe8cO8t8777TaIK1uQkJCSE1NRaFQkJSUxFdffVVquzlz5uDq6sr69evN9lLnzp1p3rw5S5YsYdGiRVZf09e35JPHdu3aVWwAN7DZcpo8eTKzZs3CZDLx0EMPERAQwJUrV1ixYgUfffRRpazn2kiAlzujOzTk18OXyMg18GN0PI/2tm5pXeBg6BpAz2egx9NyIoJjP0PMr6C/KtfnZcDfK+XD3Vc2ctuOgZCeTvN4VeD4FF9htYsP65VjELtNPvdpLL9vHZio+Cj+t+d/ZORnAPLGquc6P8eDbR5EqXBsQ1vgGEhGI3lnzpB96DDZRw6Tc+gwhlusIGoaNULbqSPuHTvi3rETrs1K+p46KtYE6jcYDKxfv57//ve/Fot7ISEh9O/fnzVr1lhtsJ46dYr58+ezfft2kpOT2b9/P506dWL+/Pn06dOH/v37V2gcNs92M2bMIDY2lo8//tji8bokSTz++OO89NJLFRKkNvNEn6b8evgSAF/tOs9/wxvjohYTq9OiUEDAnfIxeAGc3wkxv8hJCfLkL1FyUuDwMvnw8JcjDbT5j2y8OtHOa0HVkpKSQkGBvOtfqVSi1WpRq9WYim3cKJ7tpqCggPTsPPNrTzc1BoPBYuetSqVCqVSa+y2kcAPHzY8Z1Xs/ovDrzNh1mhzyyWhCo9FgMpkwFgvBo1AozPKVVm40GkvIrlKpyiwvS/bSyk2YeO/Qe3x38jtzeaBnIIt7L6aNTxuMBiNGjOWOtXBMhbopKCio0THZdJ9uMaaavE+FYyj8X0cbkzE3V/Y/PXyY3CN/kfvXX5j0espErca1VSvcOnbEo3MnPDp3RvL2tmyjVCJJUpWMqfB/jTeFvyocU3Z2Njk5RS5o+vLGYiWxsbHk5OTQvn37EnXt27cnMjKS3Nxc3NzKj/189OhRevfujU6no1+/fvz8888Wcn7++efVZ7AqFAq++OILXnjhBaKiokhOTqZu3boMGDCAFi1aVEgIR0KtVjN8+HC77s5rXl/HoNb12HIykasZufx+9DL3dan4Tjlnoir06VAoVRDaXz5GvAvntsjZs05vBMONCSXruhw269A3oPWDVsNlH8EmfeRsXDZQ6/VZjTiCLm92oZo7dy6jR48mPj7eXNayZUtatmyJn58fkZGRHI1XUrj9wNNVzc6dO8nMzDS3Dw8Pp169emzevNniS7N///64u7sTERFhLnPPT2Lwid8AyFPriLzqgzEiArVazYgRI0hKSmLfvn3m9jqdjgEDBnDx4kWOHj1qLvf396dHjx6cPXuW06dPm8uDg4Pp2LEjx44dKzGmVq1aER0dbeEb16FDB0JCQkqMqUXHFiw6sYjD1w6by+7Q3ME7fd+hnlc9izEBDB8+nJycHKKiosxlpY0pMjKyxsZky32yZUxQc/cpLS2N+vXr1/iYFPn5hOTkEJiWRmLUdjh7FmU5/qAmFxfU7dri27MnpxUKUvz8kFzkuTm8Uye8/P3ZsGFDtY5p+PDhnDt3jjNnzpjLC+/T888/z9KlS8scT0VITpYT6ZT2KN/X1xdJkkhNTSUgoPxoR7NmzTIbuC4uLvz000/muq5du/Lbb79VWEaFVJmgWE5Aoc9Zenq6VT6skiSRmZmJTqercL7b0jgUl8K9n8tv1Gb1PNn8XB+USsfMpGRPqkqfDk+eHs5sghNr4WwkGHJLtnH1ghZDoNUIaDZIjlJwC25bfVYBFdGlrfPJrfo5f/48Op0OKH+FValUkpaWhoeHB69FnOa7/bJRsXpaD9oH6iq8cqeMfBVV9BcAGHv/H6Y+/2euc4SVO4BD1w7x8p6XSc6Vv1DVCjUvdHqB8S3Gm9OA2rrKZTAY0Ov1eHp6olQqxQprJcckSRJ6vR4vLy9UKlW1j0mRm4v+4EGyow+Sc/gweSdOQDkGqsrPD/dOHXHr1Am3jh1xbdECtaurw9wngJycHLRarYUs5a2wNm7c2Kp5KSkpCX9//xJRAvbu3UvPnj358ccfGT/eMmHIm2++ySuvvMKVK1do0KBBuf3rdDpWrFjB6NGjMRqNaDQaDh06RKdOndi5cydDhw4lOzu73D7KQizT3ITBYCAqKorhw4fbNSdyl8a+dAnx4dCFVM4l6tl6KpHBberbrX9Hpar06fC4espZstrdKxuvZ/+Ef9bKK7AFNz6seekQ87N8qFzljVqtRkCLobK/bCnctvqsAhxBl76+vqV+wahu8o0rKChg586dDB8+nOyCIoNC56ouc4W4rDGZy3NS4a8V8rnaHVX3J1Hd9D+FxvLNlFWuUqlKyF5eeVmyq9VqJEni+xPf8+7hd83pVetr67Ok7xI61OtQ+phuorRypVKJQqFg165dFivs1TEma2W0tbwm7xPI789CfapUqiofk6TXy/6n0dFkHzxI7smT5cY/1QQFoe3SBW1YF7SdO6MJCSnzR6oj3KeCgoJy5yatVotWq72lDLZQt25doGiltTgpKSkoFAq8b3aLKAVJknBxKf3JYWpqKq6uFQ/taZXBqlKp2LdvH127djV/2MtCoVDUqtAW9uTJvqE8+p0ckP7zHbG3hcEqQDZe246Vj/xsiN0q+7ue3iQbrSAnKDi7WT4AAjtCy+Gy8dqgnew3KxBw06arykQJOPi1nI4Y5ExuHnUrKZn9yDHkMG/vPCLOFz1uDQ8I560+b+HrVr1xjwU1jzEjo8hAjY6WDdRyHg67hIbKBmqXzmi7dEFzi8fYAggNDcXd3Z2YmJgSdTExMTRr1uyW/qsg+7uuWbOGYcOGlajbtGkTnTt3rrCMVhmsc+bMISgoyHwuHkVWjAGt6tG8nidnE/UcvpDKvthkwkMd50tCUA24aOUNWK1HgSFfztt+agOcioDMhKJ2CX/JR9TroAuE5oNl94FGPWpOdoFDoC8W1qrCcVgLcuGA7AqAQgnh0+0gmX24rL/Mc1HPcSrllLnssXaPMb3DdFQi2sZtgVGfRc7hQ2QdiCb7wIFbrqC6tmiBNizsxtEFdV3xvWorarWaUaNGsXr1ahYvXmx2VYqPjycqKornn3/eqn6effZZJk6ciIeHBw8++KC5j23btvHNN9/w66+/VlxGaxrNnTvXfF5aZoTaRlVtwlAqFTzZN5QXf/kbgNc2nOCPp3qhquW+rGKDUBmoXSB0gHwMXyIHbj+9Cc5shKvFfuVmJsCRb+HIt6hVLvT0aI5y/3loOUROoSl+QFYYZ3pvFspaGIdVoQCtpoIG3D+rIStRPm8zGnwdI3529JVoXtzxIml5aQC4q915vdfrDA4pGfe7sjjTvXcGKqNPU04OOX/9Rdb+A2QfOEDO8eNwk09ncVxbtULb9YaB2qULah+fCl/bUbH3+3Pjxo1kZWWZN8mdOHHCbDwOHz4crVbL/PnzCQsLY+TIkcyaNYvc3FzmzJmDn58fL774olXXGT9+PLGxscybN48PP5TjO48dOxa1Ws38+fMZNWpUhcdg86arBQsW8OijjxIYGFii7sqVK3z55ZcWqbhqGnttkrAXRpPEyI92c/KKHP5o8dj2jAu7PSIGCGwg7aK8aevMnxC3q/RNWwB1gqDZAAgdKEcdqOY0sbcb9t50VZF+Br6zndjrWehc1cTMH1IxAb4cAJdv7LifEgmNulasHzvy8+mfefPAmxgk2SAP1gXzQf8PaObTrIYlE9gbqaCAnJgYsvbtI3v/AXKOHkW6aVNTcVxbtkTbtatspNZSA7UyWDOfNG7cmAsXLpRad/78eXMG08OHDzNz5kz27dtnkZo1NDTUKlny8/NRq9UkJCTw559/cu3aNfz8/BgyZAghISEVGl8hNhusxf1Zb+bw4cN07dq1xG63mqTwRrZo0QKVSmXOWFMWJpOJpKQk/Pz8SnWOtgf7YpOZ8OV+APw8XYl6qa/9MtY4GNWhz1pPfrZstJ75E+nsnyjSL5XeTqGEwE5yiK2m/SCoq81hs24nbHlvfvLJJ3zyyScYjUbOnDlT7QZrcVm7v7mNxMw8Arzc2PfyQNsvnnAUlvaVzxu0gyd21egqvcFk4O2Db7Py1EpzWc+GPVncZzF1XKpmkUHMS/blVvqUTCbyzpwha99+svbtJfvQYaRydoq7hIbi0a0b2m7d0HYNu+0MVFvfn46yMJebm4uHhwe//vor99xzj937t3nNuTz7Vq/XO+zOZWtTsxqNRvbt28fw4cOrbCILD63L0DsasOmfqyTp8/h0eywzh7aqkmvVNNWhz1qPi1b2X20xBEN+PjvXfE2/hgWozkdB3B55wxaAZILLh+Rj59ug0cqJCpr2hSZ95Rzx4h6YseW9WdOpWYvLWujDWuEUz4e+LjrvMqVGjdWM/Axm7JjB3oS95rL/tvkvL3R+oUr9VcW8ZF9K02fB5ctk7dtH1t59ZO3fjzElpcz/1zRqhEf3bmi7dkPbrSuaevWqS3SHxFnfn25ubtStWxcPD48q6d+qGe/YsWMWQW0jIiI4deqURZucnBx++OEHq5eNb3deGd6abacSyTea+HrXeSaEBRNcV3vrfxTc3igU6N0CMXUbjqrXM/Lqa/xeiI2S02smnihqW5AN5yLlA8DdBxr3ll0HGvcG/5bC/9XJMJoksvPlJ1i6imy4yk2X0wkDuOig3X12lM42LmVeYvrW6fyb/i8AaqWaOd3ncE9z+6/MCKoeZU4O+q1byT1wgKy9eym4EF9mW5WfHx7du+MR3h1tt+64BDWsRkkFVcmoUaNYs2YNd911l937tmrGW7NmDfPnzwfksFULFiwotZ27uzvLli2zn3S1mOC6Wh7p1YTPd8SSbzTx5saTfPZAxcM9CG5TXLRy0oFmg+TXGVfg3+1Fh/5qUducVDj5h3yAnDK2cS/5COklDFgnIMsiQkAFnmb9/WNRHOA777cqWUVVEHM9hqe2PUVKrrzq5uPqw3v936NzfTEHOgtSQQE5x46RtWcPmbv3EBoTw9UynsAqPTzQdu2KR3g4HuHdcWnWTEQbqqXcf//9TJkyhUceeYQxY8YQEBBQ4l536tSpQn1bZbA+/vjjjBw5EkmS6Nq1K8uWLaNt27YWbVxdXc1xvJwZhUJRbVmEpvcP5dfDl0jS57Hx+FX2xibRI9Svyq9bnVSnPm8HbqnPOgHQYYJ8SBJcPw3nd8L5HbIfbG56Udus6/DPGvkA0NaF4HDZjSAkHOq3A1Xt3UntTO/NQlmz8ov2B+hsdQmQJDn2aiFhU+wknW1svbCVWbtmkWuUNxI2rtOYTwd9SiNd9W0+daZ770jkx8ej372brD17yd6/H1NWlrnOQpNqNe4d7sSjRw88wsNxb9cOhYjKYDXO/P4cMkTeCLp8+XK+/fZbizpJklAoFBXe52TzpqsdO3bQuXNnPD1r5pe5rTiKM3JZ/HQwnpm/ySGMGvm6s/HZPhX3TRMIysNkhCt/Q9xu+YjfB3kZZbd38ZR3jweHQ3B3aNgZXKrGN8lZqOkoAaevZjLk/Z0AjO/SiEX3trf+oud3wbcj5fPgHvDIRltEtgvfn/ietw++jYT8tdOlfhfe7/8+Xq7V7xcsuDVGfRbZB/bLRuruPRRcvFhmW9fmzfDo0QNteDgeYWEoq8iPUVASR7JzbjZSS+Ohhx6qUN82W0Zt2rQhISGBFi1alKg7c+YMvr6++Pk57yqhyWTi4sWLNGrUqFqcne/t3IhfD1/iYFwqF1NyWLjuhG1fQg5OdeuztlMpfSpV0LCTfPR8BowGuHpMNl4v7JV9YYuvwObrZb/Y2G3ya4UKAtpDo+4Q3E2OQuDlvL5nzvTeLJQ13VS0UGBz0oBDNbe6apJMvHf4PZb/s9xcNqLpCBb0WICLqvojWTjTva9OJJOJvFOn0O/aTdbu3WT/9ReUkblS5esrr6D27Il7925czc/HX+jTLjjz+7Oixqg12GywTp8+HS8vL7788ssSde+88w4ZGRmsWrXKLsLVBEajkaNHjxIYGFgtbxSVUsE793Vg2Ac7yco38tOhiwxqU7/WpG2tbn3WduyqT5Xa0oA1meRNWxf2yKuvF/ZZ+sBKxqIMXAc+k8vqNISgMHklNigMGrQHza3T9zkCjvDeDAsLsyrcXqGsHs3CzGU2PYnJvCanAwbQ+smZ1qqJAlMBc/bMYf2/681lj7d/nKc6PFVjjzwd4d47Csa0NPR79pC1azf63bsxJiWV3lCjQdupEx49e+LZqyeurVqhKIwIUFDA0YgIoU87Ye37s3i4PUfk9OnTJCUl0aFDB7tEDrDZYN2zZw8ffPBBqXVDhgzh2WefrbRQtxvBdbXMGdXG7Brw8upjdAruQ11P1xqWTHBboVRCg7by0e0J2ecxNU42XuP3w8VouH7S8n8yLsOJy3Bi7Y0+NHJsz6AwCOoiuxH4NhWbucrA2nB7hejzivmw2rLCevQHMN1YKev0X1BXz9ySXZDNC9tfYE/CHgCUCiWvdnuVcS3HVcv1BSWRTCZyT5xEv3MHWTt3kXPsWJlpTzUhwXj26o1Hr554dO0qHvM7GDUdbq8svvvuO1555RWuXLkCyPNcp06dGDduHIMHD+axxx6rUL82G6xJSUnULSNPr4+PD9evX6+QILc747o0IvLENbacTCRJn8/Lq2P44sHOTul0LaglKBRyyk7fJtBholyWkwoXD8LFA3ApGi4fkV0HCjEVyClmE45A9I1c9W7esuHasBMEdpSTG9QJqPbh1Ab0xaIEWG2wSpIcHaCQTv+1s1Slk5qbyrQt0ziefBwAF6ULi/ssZmBIBZIdCCqFMSODrL170e/YiX7XrjJXURXu7nh064ZH71549u6NS3BwNUsqcHZ++eUXHn74YUaOHMmwYcMsnhx16tSJn3/+ufoM1vr16xMTE0P//v1L1MXExJRpzDoLCoUCf3//ajcUFQoFb45pz1/v7yQ5K5/NJ67xy6FLTp+2tab0WVupcX26+0CLu+QD5I1ciSfk1dfLh+HSQUg6Y/k/uWkQu1U+CtEFQECHGwZsB/lcV71uMDWuSxsolPVKsRVWT1crw1pdOQpJp+Xz4B7yD5Aq5lrWNZ6IfILY9FgAdC46Ph7wMZ3qVyycjb1xpntfESRJIj82Fv2OHei37yD7yBEo47GxS7NQPHv3wbNPb9w7d0bpYrtPcW3XZ3XjzPp88803mTx5Ml9//TVGo9HCYG3dujUfffRRhfu22WAdOnQor7/+OkOHDrXYeHX27FnefPNNRo4cWWFhqhJrfcXUajU9evSoRsmK8Ne58saYdjzxvZzj+3+/H6dFAx0dGnnXiDz2oCb1WRtxOH0qVbILQIN2RRt5ctJk4/XykRt/D8khtIqTeUU+zhTbqa4LhIA7LY86gVXmTmCLLmvaV6xQ1v2RRT8GrN50VXx19c777SxZSS5mXuSxzY9xWX8ZAH93f74Y/AXNfZpX+bWtxeE+R3bAlJdHdnQ0+qjt6HfsoODy5VLbKdzd8ejeHc++ffDs3RtNw8pvnKyN+qxJnFmfJ0+eZNGiRaXW+fr6kpycXOG+bTZY582bx/r162nfvj39+/cnKCiIS5cuERUVhZ+fnznBgKNhS2rWs2fP0rx5c1SqqksNWBZD7mjApG7B/HAgnnyDice/O8S6p3tRv45zbGS5mZrWZ23DKfTp7g3NBsoHyI+k0y/Km7UuH7mxceso5KVb/l9mgnwUN2K1deWNXA3ayQZsg3ZQt5lsKFcSW3RZ075ihbJm5hRLHGDNpitjQVFmK5UrtBldRRLKnE09yxORT3A9R/6B0tCzIV/e9WW1xli1Bqf4HFmB4fp19Dt2kLl9O1l79yFlZ5faThMSjGffvnj26Yu2a1iFVlHLo7bo01FwZn1qtVrS09NLrbt8+TI+Pj4V7ttmgzUwMJBDhw4xe/ZsNm3axNatW/H39+eBBx5gwYIFBAYGVlgYR8BkMnH69GlCQ0Nr7I0yd9QdnE3UE30+hcTMPB7/7hA/PRGOm8a53rjgGPqsTTilPhUK8A6Wj0KDyWSC1POy8XrlqGzAXvm7ZFzY7GT4N0o+ClG7Q73WsvFa/8YmsXptZEPZBpxJl4WyZuQV+RRa5cN6bitk3/BXbDXcZh3ZwonkEzwe+TjpN36INPNuxheDv6Ce1vHywjvTvS+OJEnknTmDfts2MqO2k3vsWOkNNRq0XTrLRmrfvrg2qVo3EGfVp6PizPrs2bMnH3/8MWPHji1Rt3z5cvr161fhvisUoT4wMJCvv/761g0FFcJFreSzSZ24++M9XE7L4e9L6by8OoZ3x93plD4tAkEJlEqoGyof7e6Vy0wmSIuTDVfzcazI4CrEkFO0sas4XsFQvw3Uv0M2YOvfIa/GqiqQwtRB0efauMJ6rJg7QPuqcwc4dv0YT0Y+SWZBJgBt67bls0Gf4e3mXWXXvF2Q8vPJPnSIzG1R6LdtoyAhodR2Kl9f2UDt1w+Pnj1QOUlyH0HtYs6cOfTq1YuuXbsyceJEFAoFq1evZu7cuezcuZPo6OgK9y1SKjkodT1d+eqhLoz9bC/Z+UbW/HWZ5vU9mdavWU2LJhBUDUqlHALLtynccY9cJkmQeVVOcHDlGFyLgavHIeVf4KYkfenx8nFmU1GZygX8Wsgrsv6tZEO2Xivwblxdo7IrNoW1ykmDUxHyudavyEXDzhxNPMqTW54kq0BO09mpXic+GfgJni7CYKooxowM9Dt3od+2Ff3OXZj0+lLbubZsiWf/fuj698etXTtzXFSBoKbo0qULGzduZNq0abz44osAvPHGGzRv3pyIiAjatm1b4b4rZLCmpKSwcuVKTp48SU5OjkWdQqGwafU1MzOThQsXcvToUf766y+SkpKYO3cu8+bNK9H2yJEj/N///R/79+9HrVYzYMAAlixZQtOmTSsyjFJRKpUEBwc7RPDj1gF1eHdcB55cIW/CWrzpNF7uGiZ1C6lhyazHkfRZG7jt9KlQyCGw6gRAiyFF5Xl6OTrBteOyAXvtH/nIz7T8f2O+3ObacctytTtq/5b0UtdHtf9f6PmsbDA7KIX3PetSUe52D5dbTN8n1oIxTz5vd1+VrDQfunqI6Vunk22QfSe7NujKRwM+QqvR2v1a9sQRP0cFV6+SuXUr+q3byIqOLj3DlEaDR1gYngMGoOvfzy4bpuyBI+rTmXF2ffbv35+TJ08SGxvLtWvX8PPzKzU7qq3YbLDGx8cTFhZGdnY22dnZ+Pn5kZKSgtFoxMfHx+YNCcnJySxdupQ777yT//znP3z11Veltjt16hT9+vWjQ4cO/Pzzz+Tm5jJnzhx69+7N0aNH8ff3t3UopaJSqejYsaNd+rIHQ9s2YMaQlrz9pxyW5n9rj+OmVjG2c1ANS2YdjqZPZ0fo8waunnJ2rUZdi8okCdLiiwzZaycg8SQkny0Kml+IIQfFlaPUBUg7Br2fr07pbabwvuujdgCyO4BSeQv3oL9/Kjq/c7zdZTp49SDTt04nxyAvWvQI7MH7/d/HXe1u92vZG0f5HOXFxpK5ZSuZW7aQGxNTahtlnTp49u2LbuAAPHr1cshH/Y6iz9qCM+vz+PHj5lXU0NBQQkND7da3zQbrrFmzuOOOO1i/fj2enp5s3LiRtm3b8uWXX/LGG2+wYcMGm/oLCQkhNTUVhUJBUlJSmQbrnDlzcHV1Zf369ebd/p07d6Z58+YsWbKkzDAKtmI0Gjl27Bjt27d3GGfnaf1C0ecZ+Gx7LJIEM379G1eNkpHtHX+DmyPq05kR+iwHhQJ8QuSj5bCickM+JJ+TDdnEk3D9FCSeREo9j0IyIfm3wtE9wwvvuz63ALDCfzU1DuL3yuf+reQ4t3bk8LXDFsZq74a9ea//e7iqnCM7X019jiRJIvf4cTIjt5AZGUn++fOlttMEBuI5cCC6gQPRdu6EQuPYfthiXrIvzqzP9u3b06VLFyZPnsyECRPw9va2W982rzfv27ePqVOn4uYmh1mSJAkXFxemT5/OlClTmDFjhk39KRSKW24kMhgMrF+/nrFjx1qEpgoJCaF///6sWbPG1mGUiclkIj4+HlMZqepqAoVCwf8NacnDPRoDYJLguR+PEnniWs0KZgWOqE9nRuizAqhd5M1Y7e6FgbPh/h/gmSMYZlxge8sFGPq9WtMS3pLC+555Y9PVLWOwHvu56Lz9eLvGsj2aeJRpW6aZjdW+QX15v//7TmOsQvV+jiSjkazoaK6+/gbnBgwk7r5xJC9dWsJYdW3dGr+nnqLJmtWEbt1Cg1dfwaN7N4c3VkHMS/bGmfX5ySefoFQqmT59OoGBgUycOJHIyEi79G3zCuu1a9cICAhAqVSiUqnIyCgKQ9O3b18+/PBDuwhWnNjYWHJycmjfvn2Juvbt2xMZGUlubq7ZiK6NKBQK5o5qQ57ByKroixhMElNXHGbJfXfyn46O4cckEDgVGnfStY3tvvpYVZgkyMqXN13dcoX1+Oqi8/bj7CZDzPUYntzypNlntWfDnrzb711cVPaN6+nsSAUFZO0/QObmzWRu3YoxJaVkI4UC986d0A0ahG7QIFyCnMPNSyAoj6lTpzJ16lROnjzJsmXL+OGHH/jxxx8JCgri4Ycf5qGHHqqwm0CFUrOm3PjwNW7cmEOHDpnTtMbFxaFW2z/wQGFmBF9f3xJ1vr6+SJJEamoqAQFl5ydPSUmhoEB+nKZUKtFqtajVaotfMMUdnAvbguxPolQqMRgMSJJUorx4W8CsA8NNTvNllWs0Gkwmk0UWHYVCYZavePm8ka3JLTCx5q/LGEwSz/10lKvp2TzWuykqlQqj0VhiTCqVqkzZq3pMhf0UFBSUOabC8rJkd7QxWXOfqmpMhWMo/N/aMKaby6trTMWvU5bs2dnZFhtL9WXs1q4OigUIKD9CQNJZuH5SPm/UDbzsYwidSD7BE5FPmKMBdA/ozvv93hfG6g1M+flk7dlD5p+bydy2DVNGRslGGg0e4d3RDR6MbsAA1E6eylwgKIvWrVuzePFi3nrrLTZu3Mi3337L4sWLeeONN0rM2dZis3XZvXt3/vrrL+6++27GjBnDggULyMvLw8XFhbfffpsBAwZUSBBrKM914FZuBU1uCpw8d+5cRo8eTXx8vLmsZcuWNG/eHK1Wa7GE3aFDB0JCQti5cyeZmUW7kMPDw6lXrx6bN2+2uAH9+/fH3d2diIgIi2sOHz6cnJwcoqKKgqCr1WpGjBhBUlIS+/btM5frdDoGDBjAxYsXOXr0qLnc39+fJfeFk5+dyYbT8oT41qYznPj3Mu8/3Jdjx46VGFOrVq2Ijo7m+vWi9JjVPabIyMhyx9SjRw/Onj3L6dOnzeXBwcF07NjRYcd0q/tUlWNKSUkhICCgVo2pJu6TVqtFqVSWOabnn3+epUuXUlVYmzJaqVTSsHEoHLwgj6c8g/XkH0Xnre+2i5yxabE8EfmEOc5q1wZd+XDAh7ipnfOpllKppGXLlpXehW3KyyNr924y/vwT/baoUsNPKdzc8OzdG91dd+HZry8qna5S13RE7KVPgYy1+qzplNHWoFQqadasGU2bNsXb29viu8BWFFLxJQ4rOHz4MHFxcYwdO5asrCwmTJjAhg0bkCSJPn36sGrVqnJXOssjKSkJf3//EmGtTp8+TatWrfjkk0+YNm2axf/MmDGDd955h+zs7FJdAgpTKZ4/fx7djYmivBVWR14RKl5uMBj4OCqW97eeM9fdfWcgb91zBxpVkfHuTGOqLSt3Yky1Z0ylrbA2btyY9PR0q1I9l0XhvGRLP2euZXLXezsBGNcliMX33ll6wy/6ytnDAJ79G3waV1hOgEuZl3ho40Mk5iQCcpzVzwZ95vChq6oKs5G6cRP6qChMWVkl2ig9PPDs3x/dXYPx7N0bpbvjR04QOC8VmU+qmszMTH766SeWLVvG/v37cXV15Z577mHy5MkMGjSoQn3avMLauXNnOnfuDICHhwd//PEHGRkZKBQKs0Fob0JDQ3F3dyemlLAfMTExNGvW7Jb+q76+vqXeyJt34BkMBqKjo+natWsJ94ay3B00ZTjF21KuVCpL/TVVVrlarea5wS1p4OXOK2tiMEnwx98JnE/K4tNJnWjkqy3RvjSqekwmk6mEPssak0qlKnVHZFnlNTUmW+6TvcdU/P2pVCprxZiskdHWcmvGdPNnvTTZtVotWm3RZ6ksGaoag8HA3oN/mV97upYhR+qFImM14M5KG6vXs6/z2ObHzMZqa9/WfDzwY6c3Vsub50vDlJ9P1u49ZGzaiH7rttKNVJ0O3YAB6IYOwaNHD5SuzrMJrbLYqk9B+TizPqOioli2bBmrV68mOzubsLAwPvnkEyZMmGBz2NObsUkTOTk5NGvWjM8//5xRo0aZy6vaoler1YwaNYrVq1ezePFis2EcHx9PVFQUzz9vvxiKkiRx/fp1bFx4rjHu7xpMXU9Xnl51hNwCEzGX0xn18W4+vL8jfVrYJzZtZXA2fTo6Qp/2w5l0KUkSV5PSANmoLjNKwMl1ReeVdAdIy03j8cjHuaS/BEBTr6Z8MfgLdC7O/0jbmnsvb5zaT8aGCDK3bsWUmVmijbJOHXQDB1Jn6BA8wsNRuNye/rzO9FlyBpxZnwMHDqRevXo8+eSTPPLII7Rp08ZufdtksLq7u5OTk4OHh4fdBADYuHEjWVlZZn+2EydO8OuvvwKyn5pWq2X+/PmEhYUxcuRIZs2aZU4c4OfnZ07/dbsyuE191kzryZMrDnMhOZu07AIeWhbNS3e1ZGrf0FsHGBcIBA5PbvFNV2VFCSjuv9pmdIWvlV2QzfSt0zmXJrscNfRsyNLBS/Fx86lwn86AZDSSffAQGRERZG7ejDEtrUQbZZ066AYNko3U7t1vWyNVICiNNWvWMHLkyCqJH2vzWvPAgQPZsmWLXTdXTZ06lQsXLphf//LLL/zyyy8AnD9/nsaNG9OqVSu2b9/OzJkzuffeey1Ss9ory5Uz0zqgDn881YsXfz7KlpOJSBK8/edp9pxL4u377qSht/ChEgicmdxbRQnIuAIXD8jn/q3Br3mFrlNgKuClHS9xLOkYAH7ufnw5+Evqe9SvUH+OjiRJ5P79N+kREWRu3IShlE0hSk9PdAMHohs2FM8ePYSRKhCUwejRo5EkiaSkJBQKBb6+vrfcFG8tNhusr7zyCmPHjsXNzY0xY8YQEBBQQpjSwk+VR1xcnFXtOnfuzJYtW2zq21ZUKhUdOnRwuuwSAF7uGpY+2IVPos7x7pYzSBLsjU1m6Hs7mT/6Du7p2NBubxxrcWZ9OiJCn/bDmXSpUqnwbdAQ/r0KlOEScGp90XmbirkDSJLE/L3z2XV5FwA6jY7PB31OozqNKtSfo6JSqWjv7U3KRx+TGRFBwaVLJdoo3N3R9e9PnRHD8ejV67bySbUVZ/osOQPOqs99+/bx1ltvsW3bNrKz5VjNWq2WgQMH8vLLL9OtW7dK9W9zlIDiGxnKMn4cKcSCI+6eqw72nkvipV/+JiE911w29I4GvHZPW/w8xcQrEFQEe80nFennvcgzfLD1LADLJ4fRr2U9ywbLR0KcbGjy5G5o0M5muT448gFfxcjpsTVKDV8M/oKwBmE29+OoFCQkkL5hAxnrN5BXLIxZIQqNBo++ffAaMQLPvn1Rap17c5ng9sAR7JxPP/2UZ599FpAXF5s0aYIkScTFxXH48GEAPvjggxKRnmzB5hXWOXPmVPsqXXViMBjYuXMnffr0cbrdecXp0cyPjc/1Yf66f1h95DIAm/65yt7YJP5vaCsmdg2uFt/W2qJPR0Ho0344ky4NBgMnz/1rfl3CJSArCS7skc99mkD9tjZfY+XJlWZjVYGCt3q/VSuMVWN6Ohmb/iRj3TqyDx0q2UCpxCM8nDojRqAbPKhWxkmtapzps+QMOJs+9+/fzzPPPMPw4cP59NNPCbopa9ulS5eYOnUqzz77LF26dKFr164Vuo5Vmvjjjz/o27cvXl5eFvFRayOSJJGZmemUu/Nuxstdw7vjOjC4dX1eWRNDanYBGbkG/rf2OL8cvsTr/2lL24aVCzNxK2qTPh0BoU/74Uy6lCSJjJwCQH7CVSKs1akNIN2II9vmbrBxUWFb/Dbein7L/HpW11nc1fiuyohco5jy8tDv2EHGunXot+9Auim2L0BOo0Y0mjQJn5EjUPv51YCUtQdn+iw5A86mz3feeYdu3bqxdu3aUsMJBgUF8fvvv9O7d2/efvtt8x4lW7EqLcU999xjzgLTtGlT/v777wpdTFAzDGsXQOQLfRnTqaG57O+Ladz98W5eXRPD9cy8GpROIBBYQ26x/AglfFgtwlnZFh3gn+R/mLVrFhLyl+Oj7R5lYuuJFRWzxpAkiezDh7kyZy5ne/fh8jPPkhm5xcJYdWnaFP9nnyEkYgMXn5qO96SJwlgVCCrJ7t27mT59ermZuZRKJdOmTWP37t0Vvo5VK6zu7u5mB9q4uDjy8pzPwLE2BWJtxc/TlXfHdWBcl0b8b+1xziXqMUnww4F41v51mcf7hPJYnyZoXRz/8YNAUN04QgrE4lECPIuHtcrPhvNyBix0gdCwk9V9XtFf4amtT5FjkLN5DWsyjGc6PmMPcauN/Ph40tf+Tvq6dRRcvFiiXuXvh9fwEdS5exRubdqgUCjkbGqlJKIRCAS2k5KSQnBw8C3bhYSEkJKSUuHrWGWdtG7dmldffZV77rkHgJUrV5ZpJSsUCrsG8rcXBw8etMoZWaVSER4e7nS786yle9O6RDzTm692/8sn286RlW8kK9/Ie1vOsOLABZ4d2Jz7ugThqrbP+Gu7PqsboU/7YYsuC3/oFm5uqG5UKhVqrQ7S5QxLFgZr3C4w3lhEaD7YancAfb6eaVunkZSTBMgpVxf2XOgUexSMmZlkbNxI+trfyTlypES9QqtFN2ggXnePxqN7NxQ3+QGKz5F9Efq0L86mz7p163LhwgV69epVbrv4+Hjq1q1b4etYFSVg69atjB8/npSUFBQKRbl+FQqFQkQJcBKuZ+bx4dazrIyOx2gquqcBXm482TeU8WGNcNM4xwdGIKgOajJKwOB3d3A2UY+Hi4p/FgwtqtjwEhz8Uj4fvwJajyq9g2IYTAae2voUexLkjVrBumBWDF/h0IkBJKORrL37SF+zhsytW5FuftKnVOLRvTte/xmNbuBAlHZOcCMQOCo1beeMGzeOy5cvs2vXrjLdAkwmE7169aJRo0b89NNPFbqOVT6sAwcOJCkpiYsXLyJJEmvWrOH8+fOlHv/++++tO3RgCgoK2LBhg/zIqJbjr3Nl4X/asvn5Pgy9o4G5/Ep6LnP/+Ic+i6P4ate/ZOZWXBe3kz6rA6FP++EIugwLC6NNmzZ88skn5bYrKCjgepoeuMl/VZLgXKR8rlRDk75WXXfJoSVmY9XL1YtPB33qsMZq3r/nSXz3Pc4NGMjFxx4jIyLCwlh1bd6MejNeollUFMHffI3X3Xff0lh1hHtfmxD6tC/W6vOTTz6hTZs2hIXVbDSPF154gQMHDjBmzBiuXLlSoj4hIYExY8Zw8ODBSmUmtclhsWHDhsydO5ewsDACAwMrfFFHx2Aw3LpRLSLU35PPH+zMsUtpfLj1HFtOXgMgMTOP1zac5IMtZxkf1oiHezYmyMf2uIS3mz6rGqFP+1HTurTWVQkg1ygBCkt3gORYSI2Tz4PDwe3Wff125jd+OPkDAGqlmvf7vU9InRAbJa9ajPosMjdtJO231eT89VeJepW3N3VGjsTrP//B7Y42FXJjqOl7X9sQ+rQv1uizpl2VCunevTvvvfcezz//PBEREXTp0oUmTZoAcrbSQ4cOYTKZeP/99ysc0goqEId17ty55vPTp0+TlJREhw4d8BCPX5ye9kHefPVQF45fTufjbefY9I+cVSczz8BXu8+zbG8cQ+9owAPdQ+je1H7p1gQCQfmYTBJ5NzytPN2KhbQqXF0FaDbolv0cuXaE1w68Zn49u/tsujToYi8xK4UkSeQcOkTab6vJ+PNPpJwcywYqFZ59++I95h48+/QR6VEFAgfi6aefplOnTrz55pts376d/fv3A3KmqyFDhvDyyy/To0ePSl2jQlvCv/vuO1555RXz0u/Bgwfp1KkT48aNY/DgwTz22GOVEkpQs7Rt6MXnD3bmzLVMvtl9ntV/XSbfYMJoktgQc4UNMVdo6u/BxK7B3Ns5CG+t+OIQCKqSrHwjEvIPRF3xFdazxQzW5oPL7SNBn8Dz25/HYJJXbia1nsSY5mPsLqutFCQmyrv8f/uN/AsXStS7Nm+O15gxeI0aKUJQCQQOTM+ePVm/fj0mk4mkJHkzp5+fX7nhrmzB5tSsv/zyC+PHj2fkyJEMGzaM6dOnc+jQITp16sRbb73F1q1biYyMvHVH1YStzsiFAXt1Op1YQbxBkj6PH/bH8/3+OJL0+RZ1Lmolg9vUZ2ynhvRp7o9aZfnGFPq0L0Kf9qMiuqypTVcJadn0eCsKgGFtG/DZA53lcFaLGssRAnSB8MKJMiMEZBdk89+N/+V0qhxPu3tAdz4b9BlqZc2EsZMMBvQ7d5H266/od+yAmzbqKuvUoc6I4XiPGYtb2zvs/l4XnyP7IvRpX2zVZ01vuqoubJ6t3nzzTSZPnszXX3+N0Wi0iGnaunVrPvroI7sKWBO4u7vXtAgOhZ+nK88Oas6T/Zry5z/XWHngAvv/lWOp5RtMbDh2hQ3HruDn6cLoDg0Z3SGQdg29zB80oU/7IvRpP5xFl/piWQPMPqxxu4uFsxpUprEqSRLz9s4zG6vBumCW9F1SI8Zq/qVLpP36K+mr12BITCxRr+3eHe+xY9ENHoTSza1KZXGWe+8sCH3aF6HPkti8Tnvy5Enuv//+Uut8fX1JTk6utFA1icFgICIiQjiQl4KrWsXddwby4+PhbHmhL1N6NcHXo8gdIEmfz9e7z3P3x3vo83YUb208xd/xKWzYIPRpL8T70344ky7Tsot2xZujBFj4r5btDrDi5Ao2xm0EwEPjwUcDPsLLtfo2aEgFBWT8uZn4KY8SO2gwyZ9/YWGsquvXp+7UJwmN3EzI8mV4jRpZ5caqM917Z0Do074IfZaOzT+xtVot6enppdZdvnwZHx/HDI0isC/N6nkye2QbZg1rxY7T1/ntyCW2nkwk3yjnM7+YksPnO2L5fEcsdV1V/K08zZC2AXQJ8SnhNiAQCMpHn1f0xWX2YT1bLJxV036l/t/ha4d559A75tev93ydpt5Nq0pMC/IvXiTt519IW7MG4w1/NjMqFZ79++F977149u6NwkkCpAsEgprDZoO1Z8+efPzxx4wdO7ZE3fLly+nXr5895LI7t3tq1qpCo1IyqE19BrWpT1p2PhuPX2XDsSvsjU2iMBdBcp6CZXsvsGzvBXy0Gvq3rEe/VvXo09xPbNgSOAU1nZrVwiXATX0jnNV5uaBR91LDWSVmJ/Li9hcxSrLMU9pOYWDIwCqVUyooIDMqirSffiZrz54S9ZqgILzvvRevMfegqVevSmURCAS1C5sN1jlz5tCrVy+6du3KxIkTUSgUrF69mrlz57Jz506io6OrQs5KY0u8Q0HF8Na6MKFrMBO6BpOkz+PPf66y4e8E9v+bjOnGDufU7AJW/3WZ1X9dRqmADo286deyHj2b+XFnkJdYfRU4JDUd77D4CqunqwbObSyqbF4ynFWBsYAXt79Icq7sotU9oDtPd3y6yuQruHpVXk395RcM169bVqrV6AYNwmfcfWi7d0dhpx3DAoHg9sLmKAEAUVFRTJs2jdOnT5vLmjdvzhdffOFwK6wViRJgMBhQq9Vit6MdkCSJ5Mwcdp1LYcvJRLafTiQrv/RVKp2rmm5NfekR6kf3pnVp1UCHUinuQXHE+9N+VESXNRUl4Ktd//LahpMAfDShI6NininyYX1yDzRoa9F+UfQiVpxcAUCARwA/jfzJ7pmsJJOJrL37SP1xFfqo7SV2+muCgvAeNw7vMfc4XDgq8TmyL0Kf9sVWfYooAeXQv39/Tp48SWxsLNeuXcPPz48WLVoAsqKd/Q2bk5ODTqeraTFqDS4Y+E/HhtzTKYg8g5FDcalsP53I9tPXOZuoN7fLzDOw5WQiW07KGzK83DWENfalWxNfujT24Y5AL1zUYnVGvD/th7PosrhLQB2NEeJ2yS90gVD/Dou2Wy9sNRurGqWG9/q9Z1dj1ZieTtqaNaSt+rFk3NQbvqk+4+/Ho2cPh15NdZZ77ywIfdoXoc+SVCquSWhoKKGhoebXK1euZMGCBZw6darSgtUUBoOBqKgohg8fjkajufU/CMrlZn26qlX0bOZHz2Z+vDoCLqVms+dcErvPJbP3XBLJWUVxXtNzCthy8po5VayrWsmdQd50buxDx0bedAj2pp6uancTOxri/Wk/nEmXGTlFn4v6GTFgyJVfhPa3CGd1KfMSs/fONr+eGTaTO/wsDdqKknviBCkrV5KxfgNSbq5FnbpePbzvuw/v++5F06CBXa5XlTjTvXcGhD7ti9Bn6VhtsKanp7N27VquXbtGixYtuPvuu83ZC1avXs2cOXM4ceIEISGOlZNa4NgE+WgZHxbM+LBgTCaJU1cz2f9vMgfOJxN9PoXU7AJz2zyDiei4FKLjUsxlgV5udAj2pn2QN+0aetE20AsvrfiAC2oXxX1Y/a4X2yfQuLf5tMBYwIwdM8jMzwRgSOMhjGs5rlLXlfLzydgcSeoPP5Dz118l6rXh3fGZMAFd//4oxBerQCCoQqwyWM+dO0fv3r1JTEw0P/Lv27cva9euZcKECWzatAlvb28WL17M009XnWN/dbNgwQJ+/PFHjh8/blNqsT59+tCpUyfef//9qhOuFqJUKmgTWIc2gXV4pFcTTCaJs4l6ouNSOByXwuH4VC6mWOYXT0jPJSHmKhExV81lIXW1tG3oRZsAua82AXWop3N1elcVQe3D2uglFmGtru4vqmjcy3z67uF3OZ58HIBGukbMC59X4fd8QWIiaT/9TOpPP5UISaX09MTrnnvwmXA/rk2rJ0SWQOBMFLcd4uPjadKkCW+//TYvvfRSuf/39ddfM2fOHP755x+rrlPT0UuqG6sM1tmzZ5ORkcG8efPo0qUL//77L6+//jo9evTgxIkTPProoyxevBhvb+8qFrd6UKvVJCQksHjxYpYvX25zHtyFCxcyePBgpk6dSsuWLatISudBra6Y54lSqaBlAx0tG+h4sLu8cp+YkcvhC6kcvZTG0fg0Yi6nk33TJq4LydlcSM5mw7Er5rK6Hi60bKCjRX0drW702by+rihrkBNRUX0KSlLTurQ2eknmDR9WV/JxuXpELvRpDN6NANgab+m3uqTvEjxdPG2WJ+fvv0lZ8QMZmzZBQYFFnWvz5vhMmiQH9vfwsLlvR6Om731tQ+hTpjK2w0MPPcSiRYtYsmQJ4eHht2xf09FLqhurogQEBgby9NNP8/LLL5vL/vzzT4YNG8aTTz7Jp59+WqVCVoaK7p6bOXMmP/zwA/Hx8Ta/6QDatWtHeHg4S5cutfl/BdZjNEmcTczk2KV0jl9OJ+ZyOicSMsgzmKz6/0AvN5rV19GinifN6nkSWs+TUH9PiwxeAkEhNRUlYMynezgSn0Z35Ql+dHlNLuz4AIz+hKtZVxn7x1gy8jMAeKXbK0xoNcFqWQozUaV89x25x45ZVqpU6AYOxOeBSWjDwsRTCoHgFtxsO8TFxVm9wgrwzjvvsHDhQhISEtBqtVZdU0QJKMb169fp2bOnRVmvXvKjqPHjx9tfqhrEZDKRkJDA119/zSOPPGJhrM6fP5+IiAjOnj2LwWCgWbNmTJ8+nUceeaTERP7ggw+yYMEC3nnnndt6p5/JZCIpKQk/P78KGf63QqVU0KpBHVo1qMO4LvJqk8Fo4nxSFieuZHAiIYMTVzI4eSWTJH1eif9PSM8lIT2XnWcsY0f6aDU08fOgiZ8nTfy0NPHzpLGflpC6HjW6KlvV+rydcCZdFkYJ6K0ptqG1cW9MkolXd79qNlYHBQ/i/palp86+GUNqqvzYf+VKi1SpACpvb7zHjcNnwv1oAgLsMwgHwpnuvTMg9CmTn59fqu0Aso5ef/11vvjiCxITE2nbti2LFi1i4EDLZB6TJk1ixowZLF26lGeeeea21ufNWPXNazQacbspt3Ph69pmjBmNRr777juSk5Pp37+/RV1cXBxPPPEEwcHBAOzfv5+nn36ay5cvM2fOHIu2/fr1Y+bMmWzfvp1Ro0ZVm/yOhtFoZN++fQwfPrzaPnhqlZLm9eVH/qM7NDSXJ+vzOH01k1NXMzlzLZOziXrOXMs0P24tTmp2AanxaRyJTytR5+fpQkhdD0J8tTS6cQTfOOrpXKs0dmxN6LO24ky61KgVaJQS3ZUnofCZWEhPlv+znOir8iasetp6zOtxa7/VvHPnSPn2O9L/+AMpz/JHnGvLlvj+90HqjBiB0q32RuBwpnvvDAh9yhw4cKBU2wHg448/JiQkhPfffx+TycTixYsZNmwYO3bssHj836BBA1q2bMlPP/3E9OnTb2t93ozVS0WnT5+28FEpdPItLYRVp06d7CCafbElNWthQoSbx7Fs2TLzuclkol+/fkiSxAcffMDs2bMtvig6duyIQqFgz549t7XB6kjU9XSlRzNXejQrCmIuSRKJmXmcvabn3yQ9sYl6Yq9ncS5Rz9WM3FL7SdLnk6TP5/CF1BJ1Liolgd5uBPloCfJxp6G3O4E3jiAfd+rXcROxZJ2Qmt7csHZqOJvWr6XT8bNgBHwa848xg4/++ggABQre7PUmXq6l+7FJkkTWnr2kLF9O1u7dlpUKBZ4DB+D73/+Kx/4CQSXYt28fULoNZDQaiYyMNC/2DRkyhMaNGzNnzhwiIyMt2nbs2JGNGzeW6ON2x2qD9eGHHy61/MEHHzSfF0YQcMQda7akZk1JSUGhUOB3U3aWbdu28cYbb3Dw4EEyMjIs6hITE6lfv775tUajwdvbm8uXL1deeEGVoVAoqF/Hjfp13OjV3PJ+Z+cbiEvKJi45i/NJWcQlZXEhJZsLyVlcyyjpXgCQbzQRl5xNXHJ2GdcDP09XArzcbhzuNLhx3qCOGw28ZFncNCq7j1VQcRxhc4Nv1jkURjkea3ZIOLN2zsJgkp8OTG47ma4BXUv8jyk/n4x160hZvpy8s+cs6pQeHnjfey8+Dz6AS1BQ1Q9AIKjlJCQklGo7AIwZM8biSbVOp2PUqFGsWrUKo9GISlU059erV4/09HQMBoOIw1oMqwzW4iuLtR2FQoEkSWg0Gos3UHR0NHfddRf9+vXjyy+/JCgoCBcXF9auXcvrr79OTk5Oib7c3NxKLb+dUCgU6HQ6p1y10bqozWG2biYn30h8SjYXU7K5mJptPr+UmsOl1ByLMETFkSS4npnH9cw8jl1KL/PaddzUZkO6ns4V/zqu1NO54eehJqFAS1xyNvW9tehcRSrEiuJM702FQkFg/r/m12+rs4lLjgOgTd02PNXhKYv2htRU0n78kZQfVpYIS6UJCsL3wQfwGjsWlaftkQRqA850750BoU+ZnJycErZDIQ1KSajRoEED8vPz0ev1Fj+E3dzckCSJvLw83N3dq1RmZ8Iqg/Whhx6qajkcBrVazZ133smGDRvIysrC40b4lh9//BGNRsP69estfiWtXbu2zL5SU1NL/aV1O6FWqxkwYEBNi2F33F1U5pBbNyNJEhk5Bi6mZpOQlsPltJxif3O5kp5DYmYe5cXnyMg1kJGrt0hdW5xFh+THuq5qJX6ervjpXPH3dJHPPV3x9XCh7o3Xvh4u1PVwwcfDBY1KuCMU4kzvTbVaTROF/LRmp7sbvybLQfzd1e681fstNCp5FSb/4kVSln9L2urVSDf9WHbv2BHfyQ+jGzgQRSlfqLcTznTvnQGhTxk/Pz/y8/MtbIdCrl69WqL91atXcXFxwfOmH45paWm4urrWmlCh9kIETrsJk8lkNjJjY2Np3749IP+CVKvVFr+ccnJy+P7770vtJyEhgdzcXNq0aVP1QjswJpOJixcv0qhRo9vGeVyhUOCl1eCl9aJtw9IfHxcYTVzLyOVaRi5X0nO5euO4kpHL9Yw8rmXKr28VnivPYOLyDWPYGuq4qal7w4j10brg66HBx8MFX6382lsrv/bRavByl1/XViPXmd6bprwsFBejSVcqmVuvnrl8RtgMmng1IScmhuSvvyFz82YwFXvPKJXoBg+m7uSHce/QofoFd1Cc6d47A0KfMq1atQIsbYdCVq9ezdtvv21e8MrMzGTdunX07t27xIpsbGwszZo1w2Qy2UWfer2e//3vf/z888+kpKTQqlUrZs2axf33WxdRxFEQButNGI1Gc+SD/fv3m990I0aM4N1332XixIk8/vjjJCcns2TJElxdXUvtZ/9+ORtNabsFbyeMRiNHjx4lMDDwtp7IbkajUt7YmFV2nL3CldrEzFwSM/NIzMzlSlo2B4+dRuffkOTsfK5n5pGkzyc1O7/cFdtC5JVbA+eTsqyWVeeqxkurwVurwdvdRTbG3YsO7xt/6xT+dZP/erqpUVVhxITK4kzvTVP8AdSmAt7wr0vSDVF7BfZk6NX6XHj9YbIPHLBor3B3x3vsWHwffkj4p5aCM917Z0DoU6Zfv36Ape1QiEqlYvDgwbzwwguYTCYWLVpERkYG8+fPt2hnMpk4ePAgffv2xWg02kWfY8aM4eDBg7z11lu0aNGClStXMmHCBEwmExMnTqx0/9WFMFhLwd/fn169evH777/z+OOPAzBgwAC++eYbFi1axKhRo2jYsCGPPfYY9erVY8qUKSX6WLt2Le3ataNdu3bVLb6gllC0UquheX35R1RBQQENM04yfHg7C2d8g9FESlY+1/V5pGTlk6zPJ0mfR3JWPsmFZVn5pGTlk6LPJ7MMH9vSyMwzkJln4FKq7f7YOlc1ddw16NzU1HG78ffGa/mQzz1di+o9b7zWuWrwcFWhrqUrvLaguLCHP7XubPT0QGmSGHDWjam/XebSmScs2qnq1sX3wQfwuf9+VOJxokBQrTRq1IjevXtb2A6FPPXUU+Tm5vLMM8+QmJjIHXfcwYYNG0rEuN++fTvp6en07dvXLjJFREQQGRlpNlJBXki7cOECM2bMYPz48aX63DoiwmAtg6eeeopJkyZx+fJlGjaUY3lOnjyZyZMnl2j7yCOPWLzOyMhg9erVvPfee9Uiq0CgVimpV8eNenWsi52ZbzCRlp1PSrZswKZmF5CanS+XZRWQlpNPenYBaTlyeeG50WTFMm4xCo3dyuCuUeHppkbnqsbDVY2HqwpPVw2erio8XGXj1qPwcFGZ23i4yGVac5karUZVpXFyq4qUCztY5OXD4CMm7j5gon5aFkaKdv27hITgO+URvEaPRlnGUx+BQFD1PPvss4wfP95sOzRu3JjiCUVvjtl+M8uWLSM8PJzQ0FC7yLNmzRo8PT257777LMonT57MxIkTOXDgAD169LDLtaoaYbDehEKhwN/fn7CwMMLCwnjzzTf5+OOPberjvffeIzg4uFTj9najUJ+3++5Re2EvfbqobTNw4UYsz3wjadn5pGUXkJFTQHqObMim5xS9Ljwycw1k5haQkWsgPaeAfCvT5d5MToGRnAIj1zNLDyVmK24aJR4usiGLwYUN6UdZ+t8wu/RdFRjTkli/J57Xj6rwzrbUoVv79tR9dIrYSGUjYl6yL0KfRYwZM6bCtkNsbCw//fQTmzdvxsXFxS76PH78OK1bt7aIow+YXRaOHz8uDFZnRa1Wm2/el19+yR9//GGz43OdOnVYvnx5iTfI7UhxfQoqT03qU6FQ4HljRTPIx/b/zzMYycw1kGE2ZmWDNjNXXoXNzC1An2tAnyfXZeQWkJVnICvPiD6v6LWNi7wlyC0wkVuQT/INN163pNJj5lYHKSkpFBQUAKBUKtFqtajVakzFNk7Fp1/mzkMq3POL/k/bswd1H3sMj27dMBqNGEwm82YrlUqFUqk091tI4XxkMBisKtdoNJhMJou42oWbT8sqNxqNFrIrlUpUKlWZ5QaDwWL1qVD2ssrtNSZJkggLC0OSJAwGQ60YU03fp7CwMPP3ZG0ZU/FyW8a0dOlS1q5dS15enlkn1ozp33//5YMPPqBXr17mMRW/bqHs2dnZFiEz9frSo8kAJCcn07Rp0xLlvr6+5npnQVhUN2E0Gjl79izNmzenbdu2tG3b1uY+nn/++SqQzDkprk9n8ZNxZJxZn65qFa6eKvw8K/7IWpIkcgtM6PNkwzar8Mg3oM8zml9n5xvN5dl5Rvlvvmz45uQbzeX6vAI8XGtOj02aNLF4PXfuXEaPHk18fLy5rGXLlniOu5eCH34loW1jjP1HkdewIR0CAvBUKNi5cyeZmZnm9uHh4dSrV4/NmzdbfJn2798fd3d3IiIiLK45fPhwcnJyiIqKMpep1WpGjBhBUlKSOXsPyMHOBwwYwMWLFzl69Ki53N/fnx49enD27FlzpkCA4OBgOnbsyLFjx0qMqVWrVkRHR3P9+nVzeYcOHQgJCRFjctIxde3alYCAgFo1porcp5CQENq3b8+mTZtsHlNgYCDR0dH4+vpiMpk4e/ZsiTE9//zzLF26FGspb6XWmVbFFZJkzd5i56UwM02LFi2sSs1aUFBAREQEw4cPFxkm7IDQp30R+rQfhbocNmwYLi4u5bYtnpr1zJkzpKenW505rzQK56Xz58+bo5KUtcKqVCopSEpiy7p1DJg40Xzfb5eVu6oYU25uLpGRkQwePBgXF5daMaaavE8FBQVERkYydOhQXF1da8WYbi6vzjEZDAY2b97MkCFDLBYmylthbdy4canzUnh4OEajkejoaIvyf/75h7Zt2/LFF1+U2CDmqNw2K6y2pGYVCAS3D9asMFRValZfX99S56WbV89Nvr4Y/f3RaDQlfqiU5XpU1g8aW8qVSmWp7lBllatUqlJX/ssqL0v2qh5TYblGozFfy9nH5Aj3qfjjb2tlL6vcUcZUnoy2lldkTKX1o9Vq0WqLQiKWt3jRrl07Vq1aZXZ9KSQmJgagQk+RawoRL0YgEAgEAoGgFnLPPfeg1+v57bffLMq//fZbAgMD6datWw1JZju3zQqrtSiVSoKDg2/r4Mf2ROjTvgh92g9n0qUzyeoMCH3aF6FP+2JPfQ4bNozBgwczdepUMjIyaNasGatWrWLTpk2sWLHCqfZCiHfXTRgMBn7//fcS/ieCiiH0aV+EPu2HM+nSmWR1BoQ+7YvQp32xtz5Xr17Ngw8+yJw5cxg6dCgHDhxg1apVTJo0yS79Vxe3zaYrazdJpKam4uvrS0pKCj4+FYjdI7BA6NO+CH3aj4ro0tb5xF79iPtuX4Q+7YvQp32xVZ/2mpccHbHCehOFuwmL7yq0B5988old+3OWPoU+hT4dtb+q0mVVcDvf96roU+hT6NOR+3Smuak6EQZrNeEMH5Kq6rMqcJaxC306bn8C57jvVdVnVeAsYxf6dPw+BSWp9ZuuCj0eMjIyrGpfGDA4MzPTrnEujUaj1TLUpj6FPoU+HbW/iuiy8PqV9aQS81LN9in0KfTpyH3aqk97zUuOTq33Yb106RKNGjWqaTEEAkEt4uLFiwQFBVX4/8W8JBAI7E1l5yVHp9YbrCaTiYSEBHQ6nVOlIBMIBI6HJElkZmYSGBhYqZAzYl4SCAT2wl7zkqNT6w1WgUAgEAgEAoFzU3tNcYFAIBAIBAJBrUAYrAKBQCAQCAQCh0YYrDfQ6/U899xzBAYG4ubmRocOHfjxxx9rWiyHZtu2bTzyyCO0atUKDw8PGjZsyOjRozl8+HCJtkeOHGHQoEF4enri7e3NmDFj+Pfff2tAaufiq6++QqFQ4OnpWaJO6NQ6du/ezfDhw/Hx8cHd3Z3mzZuzcOFCizaOqksxL1UMMTdVLWJeqjzOPC/VFMJgvcGYMWP49ttvmTt3Lhs3biQsLIwJEyawcuXKmhbNYfnss8+Ii4vj2WefJSIigg8++IDExES6d+/Otm3bzO1OnTpFv379yM/P5+eff+abb77hzJkz9O7dm+vXr9fgCByby5cv89JLLxEYGFiiTujUOlauXEnfvn3x8vLiu+++IyIigpkzZ1qEf3FkXYp5qWKIuanqEPNS5XH2eanGkATShg0bJEBauXKlRfngwYOlwMBAyWAw1JBkjs21a9dKlGVmZkr169eXBg4caC677777JD8/Pyk9Pd1cFhcXJ2k0Gun//u//qkVWZ2TkyJHSqFGjpIceekjy8PCwqBM6vTWXLl2SPDw8pKlTp5bbzlF1KealiiPmpqpDzEuVw9nnpZpEGKySJD366KOSp6enVFBQYFG+cuVKCZD27NlTQ5I5J/3795datGghSZIkFRQUSO7u7tITTzxRot1dd90lNW/evLrFcwq+//57SafTSRcvXizxxSB0ah3z5s2TACkuLq7MNo6sSzEv2R8xN1UOMS9VHmefl2oS4RIAHD9+nNatW6NWWyb+at++vbleYB3p6ekcOXKEO+64A4DY2FhycnLMuixO+/btOXfuHLm5udUtpkOTmJjIc889x1tvvVVqEGihU+vYuXMnvr6+nDp1ig4dOqBWq6lXrx5PPvmkOTOMI+tSzEv2RcxNlUPMS/bB2eelmkQYrEBycjK+vr4lygvLkpOTq1skp2X69OlkZWXx6quvAkW6K0u/kiSRmpparTI6OtOmTaNly5ZMnTq11HqhU+u4fPky2dnZ3HfffYwfP54tW7YwY8YMvvvuO4YPH44kSQ6tSzEv2RcxN1UOMS/ZB2efl2oS9a2b3B6Ul21GZKKxjtmzZ/PDDz/w0Ucf0blzZ4s6oV/r+O2331i3bh1//fXXLfUidFo+JpOJ3Nxc5s6dy6xZswDo168fLi4uPPfcc2zduhWtVgs4ri4dVS5nQ8xNlUPMS/ajNsxLNYVYYQXq1q1b6mpFSkoKUPqvHIEl8+fP57XXXuP111/nqaeeMpfXrVsXKH01KCUlBYVCgbe3d3WJ6dDo9XqmT5/O008/TWBgIGlpaaSlpZGfnw9AWloaWVlZQqdWUqinIUOGWJQPGzYMkEPGOLIuxbxkH8TcVDnEvGRfnH1eqkmEwQq0a9eOkydPYjAYLMpjYmIAaNu2bU2I5TTMnz+fefPmMW/ePF555RWLutDQUNzd3c26LE5MTAzNmjXDzc2tukR1aJKSkrh27RrvvPMOPj4+5mPVqlVkZWXh4+PDpEmThE6tpDT/L8AcOkapVDq0LsW8VHnE3FR5xLxkX5x9XqpRamq3lyMREREhAdKPP/5oUT506FARPuYWLFiwQAKk//3vf2W2GTdunFSvXj0pIyPDXHbhwgXJxcVFmjlzZnWI6RTk5ORIUVFRJY4hQ4ZIbm5uUlRUlBQTEyNJktCpNfz5558SIL3++usW5e+++64ESLt27ZIkyXF1KealyiHmJvsg5iX74uzzUk0iDNYbDB48WPLx8ZGWLl0qbdu2TXrsscckQFqxYkVNi+awLFmyRAKkoUOHSvv27StxFHLy5EnJ09NT6tOnjxQRESGtXr1aatu2rRQYGCglJibW4Aicg9LiHQqdWseoUaMkV1dXaeHChVJkZKT05ptvSm5ubtLIkSPNbRxZl2Jeqhhibqp6xLxUcZx9XqophMF6g8zMTOmZZ56RGjRoILm4uEjt27eXVq1aVdNiOTR9+/aVgDKP4hw6dEgaOHCgpNVqpTp16kj/+c9/pHPnztWQ5M5FaV8MkiR0ag3Z2dnSzJkzpUaNGklqtVoKDg6WXn75ZSk3N9einaPqUsxLFUPMTVWPmJcqjrPPSzWFQpKK5QITCAQCgUAgEAgcDLHpSiAQCAQCgUDg0AiDVSAQCAQCgUDg0AiDVSAQCAQCgUDg0AiDVSAQCAQCgUDg0AiDVSAQCAQCgUDg0AiDVSAQCAQCgUDg0AiDVSAQCAQCgUDg0AiDVSAQCAQCgUDg0AiDtRZwzz334O7uTlpaWpltJk2ahEaj4dq1a5W+XkJCAvPmzePo0aMl6ubNm4dCoaj0Ncpi+fLlKBQK4uLiymyjUCisOrZv305cXFy5bebNm1ei/3Xr1jFq1Cjq16+Pi4sLvr6+DBw4kB9++IGCgoJy5S/e95IlSwBo3LixVfIuX778luN7+OGHS1xz165djBs3joYNG+Li4oKXlxc9evTgs88+Iysry0K2p556qlS5f/31V7POCnnuuefM1/X09Cx33ILbDzEvWSLmJUvEvCSwFXVNCyCoPFOmTGHt2rWsXLmSadOmlahPT09nzZo1jBw5kvr161f6egkJCcyfP5/GjRvToUMHi7pHH32UoUOHVvoalWHfvn0WrxcuXEhUVBTbtm2zKG/Tpg0pKSkAPP3000ycOLFEX0FBQeZzSZJ45JFHWL58OcOHD+fdd9+lUaNGpKenExUVxbRp00hKSuLZZ58tV74pU6bw6KOPEhISAsCaNWvIy8sz13/11Vd8/fXXbNq0CS8vL3N5aGio+fzee+/lxRdfLNG3v7+/xeu5c+eyYMECevTowcKFCwkNDSU7O5u9e/cyb948zpw5w3vvvVeuvGXx/PPPc//997Nw4UJ27NhRoT4EtRcxL1ki5qUixLwkqBA1mxlWYA8MBoMUGBgode7cudT6zz77TAKkdevWVfo6ubm50sGDByVAWrZsWaX6qwjLli2TAOn8+fNW/09ZOa8lSZLOnz8vAdLbb799y34WLVokAdL8+fNLrb9y5Yq0a9eucvsApLlz55bbZu7cuRIgXb9+vcw+pk+ffkt5f/75ZwmQpkyZIplMphL1GRkZ0p9//mlVv7/88osESFFRUSXqytOv4PZFzEvlI+YlMS8JbEO4BNQCVCoVDz30EIcPHyYmJqZE/bJlywgICGDYsGEAXL16lSeeeIKgoCBcXFxo0qQJ8+fPx2AwmP+n8JHU4sWLee2112jSpAmurq5ERUURFhYGwOTJk0s8oirr0dvKlSsJDw/H09MTT09POnTowNdff22uj4yMZPTo0QQFBeHm5kazZs144oknSEpKsqeqKkxBQQGLFi2iVatWzJ49u9Q2DRo0oFevXtUsWdksWLAAHx8fPvzww1LviU6n46677qoByQS3A2JeqnrEvCS4nRAuAbWERx55hLfeeotvvvnG4lHKiRMniI6OZtasWahUKq5evUrXrl1RKpXMmTOH0NBQ9u3bx2uvvUZcXBzLli2z6PfDDz+kRYsWLFmyhDp16lC/fn2WLVvG5MmT+d///seIESMAy0dUNzNnzhwWLlzImDFjePHFF/Hy8uL48eNcuHDB3CY2Npbw8HAeffRRvLy8iIuL491336VXr17ExMSg0WjsrDFLTCaTxRdjIWq1/BE5dOgQKSkpPPbYY1XqC2ctkiSVKq9KpUKhUHDlyhWOHz/O+PHj0Wq1NSChQCDmpcoi5iWBoAhhsNYSmjVrRp8+fVixYgWLFy82T6TffPMNIH9xgLzSkJqayj///ENwcDAAAwcOxN3dnZdeeokZM2bQpk0bc79ubm78+eefFhNzfn4+IPsude/evVy5zp8/zxtvvMGkSZNYsWKFuXzw4MEW7Z588knzuSRJ9OjRg379+hESEsLGjRu5++67bdaJLcycOZOZM2eWKN+1axe9evUiPj4egCZNmlSpHNby6aef8umnn5Yo//7773nggQccTl7B7YmYlyqHmJcEgiKES0AtYsqUKSQlJfHHH38AYDAYWLFiBb1796Z58+YArF+/nv79+xMYGIjBYDAfhY/lbnZSv/vuuyu1ihAZGYnRaGT69OnltktMTOTJJ5+kUaNGqNVqNBqN2fn/5MmTFb6+tTz77LMcPHiwxHHz5g1HYdy4caXKO3z48JoWTSCwQMxLFUfMSwJBEWKFtRZx77338vTTT7Ns2TLGjh1LREQE165dY9GiReY2165dY926dWVO9jf7ZgUEBFRKpuvXrwPlP5ozmUzcddddJCQkMHv2bNq1a4eHhwcmk4nu3buTk5NTKRmsISgoiC5dupRZX7jqc/78+SqXxRr8/f3tLq9KpcJoNJZaV/iYr6ofgQpqH2JeqjhiXhLzkqAIYbDWItzd3ZkwYQJffvklV65c4ZtvvkGn03HfffeZ2/j5+dG+fXtef/31UvsIDAy0eF1Zv6jCcCaXLl2iUaNGpbY5fvw4f//9N8uXL+ehhx4yl587d65S17YnXbp0wdfXl99//50333zTIfzFyiMgIIB27dqxefNmsrOzrfIXq1+/PpcvXy61rrDcHuGHBLcXYl6qOsS8JOal2wnhElDLmDJlCkajkbfffpuIiAjuv/9+i0lh5MiRHD9+nNDQULp06VLiuPmLoTRcXV0BrFphuOuuu1CpVHz22WdltimcZAv7LeSLL764Zf/VhUajYebMmZw6dYqFCxeW2iYxMZE9e/ZUs2RlM3v2bFJTU3nmmWeQJKlEvV6vZ/PmzebXgwYNIioqyrz6VIgkSfzyyy80btyYZs2aVbncgtqHmJeqBjEviXnpdkKssNYyunTpQvv27Xn//feRJIkpU6ZY1C9YsIDIyEh69OjBM888Q8uWLcnNzSUuLo6IiAg+//zzch+Tgbypwd3dnR9++IHWrVvj6elJYGBgqV8qjRs35pVXXmHhwoXk5OQwYcIEvLy8OHHiBElJScyfP59WrVoRGhrKrFmzkCQJX19f1q1bR2RkpF11Ux7x8fHs37+/RLm/v785MPaMGTM4efIkc+fOJTo6mokTJ5oDdO/cuZOlS5cyf/58evbsWeXyXrt2rVR569SpY96cct999zF79mwWLlzIqVOnmDJlijlA94EDB/jiiy8YP368OYTMnDlzWLduHd26dWPWrFk0b96cq1ev8uWXX3Lw4EF+/vnnKh+XoHYi5qWKIeYlMS8JilEz4V8FVckHH3wgAVKbNm1Krb9+/br0zDPPSE2aNJE0Go3k6+srde7cWXr11VclvV4vSdKtA1evWrVKatWqlaTRaCyCThcGl76Z7777TgoLC5Pc3NwkT09PqWPHjhYBvk+cOCENHjxY0ul0ko+Pj3TfffdJ8fHxJQJaV1WA7rKOSZMmlfif33//XRoxYoTk7+8vqdVqycfHR+rfv7/0+eefS3l5eeXKcvN4SsOaAN1lHT179izRfseOHdK9994rBQQESBqNRqpTp44UHh4uvf3221JGRoZF27Nnz0oPPPCAFBAQIKnVasnb21u66667pK1bt5YprwjQLbAGMS9ZIuYlMS8JbEMhSaWsyQsEgipBoVAwe/Zs5syZY45N6KyYTCZMJhNTpkzht99+Q6/X17RIAoGgAoh5SeAMCB9WgaCaWbhwIRqNhnfeeaemRakUL7zwAhqNhu+++66mRREIBJVEzEsCR0essAoE1cihQ4fM540aNXLq3a0XL17k2rVrgBx6pmPHjjUskUAgqAhiXhI4A8JgFQgEAoFAIBA4NMIlQCAQCAQCgUDg0AiDVSAQCAQCgUDg0AiDVSAQCAQCgUDg0AiDVSAQCAQCgUDg0AiDVSAQCAQCgUDg0AiDVSAQCAQCgUDg0AiDVSAQCAQCgUDg0AiDVSAQCAQCgUDg0Pw/kwhGOnAp9fUAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# setup\n", "vtec = np.array([0.1] + np.linspace(1, 60, 60).tolist()) # TECU\n", "inc_angle = 42 # deg\n", "inc_angle_iono = iono.incidence_angle_ground2iono(inc_angle)\n", "\n", "# plot\n", "fig, axs = plt.subplots(nrows=1, ncols=2, figsize=[7, 3], sharex=True)\n", "markers = ['-']*4 # ['s-', 'o-', 'd-', '^-']\n", "for i, (bname, freq, marker) in enumerate(zip(bnames, freqs, markers)):\n", " # calculate\n", " ref_angle_iono = iono.iono_incidence2refraction_angle(inc_angle_iono, vtec, freq)\n", " rg_delay_inc = (vtec / np.cos(np.deg2rad(inc_angle_iono)) * 1e16 * K / (freq**2)).astype(np.float32)\n", " rg_delay_ref = (vtec / np.cos(np.deg2rad(ref_angle_iono)) * 1e16 * K / (freq**2)).astype(np.float32)\n", " perc = (rg_delay_inc - rg_delay_ref) / rg_delay_ref\n", " # plot\n", " axs[0].plot(vtec, ref_angle_iono, marker, lw=2, label=bname+'-band')\n", " axs[1].plot(vtec, perc * 100., marker, lw=2, label=bname+'-band')\n", "# axis format\n", "axs[1].yaxis.tick_right()\n", "axs[1].yaxis.set_label_position(\"right\")\n", "for ax, ystep, num, xloc in zip(axs, [10, 10], ['(a)', '(b)'], [0.04, 0.88]):\n", " ax.tick_params(which='both', direction='in', left=True, right=True, top=False, bottom=True)\n", " ax.yaxis.set_major_locator(ticker.MultipleLocator(ystep))\n", " ax.yaxis.set_minor_locator(ticker.AutoMinorLocator())\n", " ax.xaxis.set_minor_locator(ticker.AutoMinorLocator())\n", " ax.grid('on', linestyle='--')\n", " ax.set_xlabel('Vertical TEC [TECU]')\n", " ax.annotate(num, xy=(xloc, 0.06), xycoords='axes fraction')\n", "axs[0].set_ylabel('Refraction angle [deg]')\n", "#axs[1].set_ylabel('Slant range delay\\noverestimate [%]')\n", "axs[1].set_ylabel('Overestimate [%]')\n", "fig.tight_layout()\n", "\n", "axs[1].legend(ncol=4, handlelength=2, frameon=False, bbox_to_anchor=(1.01, 1.25))\n", "\n", "# output\n", "out_fig = os.path.join(work_dir, 'iono_refr_impact.pdf')\n", "print('save figure to file', out_fig)\n", "#plt.savefig(out_fig, bbox_inches='tight', transparent=True, dpi=300)\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.13" } }, "nbformat": 4, "nbformat_minor": 4 }