{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Making line plots" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Import packages" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "# import plons scripts\n", "import plons\n", "import plons.SmoothingKernelScript as sk\n", "import plons.PhysicalQuantities as pq\n", "import plons.ConversionFactors_cgs as cgs" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Setting information about data" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "prefix = \"wind\"\n", "loc = \"/STER/matse/Papers/Esseldeurs+2023/Phantom/High/binary6Lucy/\"\n", "outputloc = \".\"\n", "dump = loc+\"wind_00600\"" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Loading setup and dump" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "setup = plons.LoadSetup(loc, prefix)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "dumpData = plons.LoadFullDump(dump, setup)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Making x and z lines on which to smooth" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "n = 200\n", "rx = np.array((np.linspace(-150, 150, n), np.zeros(n), np.zeros(n))).transpose()*cgs.au\n", "theta = pq.getPolarAngleCompanion(dumpData['posComp'][0], dumpData['posComp'][1])\n", "rz = np.array((np.zeros(n), np.zeros(n), np.linspace(-150, 150, n))).transpose()*cgs.au" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Smoothing the data on the lines" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "smoothx = sk.getSmoothingKernelledPix(100, dumpData, ['rho'], sk.rotatePixCoordAroundZ(theta, rx))\n", "smoothz = sk.getSmoothingKernelledPix(100, dumpData, ['rho'], rz)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Plotting the lines" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmYAAAJMCAYAAABQG2raAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAACa4klEQVR4nOzddXjd5d3H8feR5MQ9adMm1dRdkBbX4u4wAwbjgQls8IxtMNjG2DNnGzaswBjurgVKhbbU3Zu0TdO4J0efP37nnLgfSU4+r+vKdU6O/M5dJP3kvr/39zZ5PB4PIiIiIhJ25nAPQEREREQMCmYiIiIi/YSCmYiIiEg/oWAmIiIi0k8omImIiIj0EwpmIiIiIv2EgpmIiIhIP6FgJiIiItJPKJiJiIiI9BMKZiIiIiL9RFCD2X333cf8+fOJi4sjJSWlzfOlpaWcccYZDBs2DJvNRm5uLrfccgtVVVXBHJaIiIhIvxTUYGa327n00ku56aab2v9ws5nzzz+ft956i+3bt7Nw4UI++eQTfvCDHwRzWCIiIiL9kikUh5gvXLiQn/zkJ1RUVHT52n/84x/86U9/oqCgINjDEhEREelXrOEeQHMHDx7ktdde44QTTujwNY2NjTQ2Nvq/d7vdlJWVkZ6ejslkCsUwRUREZBDzeDxUV1czbNgwzObALj72i2B25ZVX8uabb1JfX8+5557L448/3uFr77//fu69994Qjk5ERESkrYKCAnJycgJ6zR4vZd5zzz1dBqOVK1cyd+5c//ddLWUeOnSIiooKtm3bxi9+8QtOOOEEHnrooXZf23rGrLKykhEjRlBQUEBSUlJP/igiIi3kl9Vy1gNfEW+zcM1RI3n0y91ceWQuvzx7criHJiL9SFVVFbm5uVRUVJCcnBzQa/d4xuyWW27hiiuu6PQ1o0aN6tE1hw4dytChQ5k4cSLp6ekcd9xx3HXXXWRnZ7d5rc1mw2aztXk8KSlJwUxE+iSu0YzZFoc1xkpcQiJmWxxRsfH62SIi7QpGCVWPg1lGRgYZGRkBH4iPbwKv+ayYiEgouL3rBxazCbP3B67LHcYBicigE9Qas/z8fMrKysjPz8flcrF27VoA8vLySEhI4L333qOoqIgjjjiChIQENm/ezB133MExxxzT41k3EZG+8v1iaDaZsJhbPiYiEgpBDWZ33303Tz/9tP/7WbNmAbBo0SJOPPFEYmNjeeyxx7j11ltpbGwkNzeXiy66iJ///OfBHJaISLtc/mDWtEThciuYiUjoBDWYLVy4kIULF3b4/EknncTSpUuDOQQRkW5ze5ctjRkzI5gpl4lIKOmsTBERL3ezpUxvLvM/JiISCgpmIiJe7mZLmb7ifwUzEQklBTMRES/fsqW5xa5MBTMRCR0FMxERL3eLXZmaMROR0FMwExHxcrubLWX6gpn6mIlICCmYiYh4tVzKNO67NGMmIiGkYCYi4tViKdNbY6YGsyISSgpmIiJeLZYyVfwvImGgYCYi4uVfyjSZmmrMlMtEJIQUzEREvNRgVkTCTcFMRMTLf1amGbXLEJGwUDATEfHyFfpbTCYdYi4iYaFgJiLi5etZZmq2K1O5TERCScFMRMTL1eysTIv3p6NbyUxEQkjBTETEy7+UaW62lKkaMxEJIQUzEREv3+SYljJFJFwUzEREvFwtzso0HtNSpoiEkoKZiIiXu9lSptmkdhkiEnoKZiIiXp7mnf/VLkNEwkDBTETEyxfCTCaTv8GsJsxEJJQUzEREvPxLmSbwTphpV6aIhJSCmYiIV/OlTItqzEQkDBTMRES8fLNjzZcytStTREJJwUxExKtpVyb+BrPKZSISSgpmIiJe7uZLmWbtyhSR0FMwExHxcvsbzJrw5jLVmIlISCmYiYh4+UKYWQ1mRSRMFMxERLyaljJp1mA2jAMSkUFHwUxExKv5UmZTg1nNmIlI6CiYiYh4+Zcym9WYqcGsiISSgpmIiFeLpUz1MRORMFAwExHxaj5jZlEfMxEJAwUzEREvf42ZdmWKSJgomImIeLVcyjTuq8GsiISSgpmIiJerRfG/ZsxEJPQUzEREvDz+szKbHWKuXCYiIaRgJiLi5ZsdM5mML9BSpoiEloKZiIiXr8t/812ZoCazIhI6CmYiIl7NlzLNzYKZZs1EJFQUzEREvJovZfoazBqPh2tEIjLYKJiJiHi1WMpsEcyUzEQkNBTMRES8fAHM0uyszOaPi4gEm4KZiIiXx9/HDNWYiUhYKJiJiHi5/DVmLYv/lctEJFQUzEREvHwBrHmDWWg6Q1NEJNgUzEREvFouZTY97lKNmYiEiIKZiIiXr5bMZDJ5v4zHVfwvIqGiYCYi4tV8KROaNgC43eEakYgMNgpmIiJe7mZLmYD/WCbNmIlIqCiYiYh4+Yr8fTNlZu9PSLXLEJFQUTATEfHy5S9/MPPeasJMREJFwUxExKujpUztyhSRUFEwExHx8gczbzLTrkwRCTUFMxERL3ezQ8yhaXemGsyKSKgomImIeDUtZbZql6FcJiIhomAmIuLVusbMt6SpXZkiEioKZiIiXv5dmf4Gs77HFcxEJDQUzEREvFovZarBrIiEmoKZiIiXy62lTBEJLwUzEREvT0dnZSqXiUiIKJiJiHj5lixNrdtlaClTREJEwUxExKv1Uqa/waymzEQkRBTMRES8/EuZrYr/dSSTiISKgpmIiFfrpUwdYi4ioaZgJiLi5VKDWREJMwUzEREvd5tdmcb3WsoUkVBRMBMR8fK0bjBrNrV4XEQk2BTMRES8fEuWvt2YvoDmcodrRCIy2CiYiYh4dbSUqT5mIhIqCmYiIl4dLWWqj5mIhIqCmYiIV+ulTJOOZBKREFMwExHx8i1ZqsGsiISLgpmIiJcvf/n6l5nNvscVzEQkNBTMRES82jSYNanBrIiEloKZiIiXu1Xxv1k1ZiISYgpmIiJebm+/Mu3KFJFwUTATEfHqaMZMxf8iEioKZiIiXv5g5v3JqAazIhJqCmYiIl6+FUstZYpIuCiYiYh4+QKYiv9FJFwUzEREvPwNZn1LmWa1yxCR0FIwExHx8uUvk0mHmItIeCiYiYh4tV7KtPiXMhXMRCQ0FMxERLxan5WpQ8xFJNQUzEREvJqWMo1bX62ZasxEJFQUzEREvPxnZZpb7srUIeYiEioKZiIiXp5WS5lNuzLDNiQRGWQUzEREvJoazBq3Fh3JJCIhpmAmIuLlqyVr3S5DS5kiEioKZiIitAxfvqOY1GBWREJNwUxEhJYtMXwzZTqSSURCTcFMRISWs2Km1oeYaylTREJEwUxEhJbhyxfIfP3M3JoyE5EQUTATEQE87SxlalemiISagpmICC3Dl9nUusFsWIYkIoOQgpmICC2XMs1tGswqmYlIaCiYiYgAnmbd/bWUKSLhEtRgdt999zF//nzi4uJISUnp9LWlpaXk5ORgMpmoqKgI5rBERNpofynT+F4NZkUkVIIazOx2O5deeik33XRTl6+97rrrmD59ejCHIyLSoRZLmWowKyJhEtRgdu+993Lrrbcybdq0Tl/38MMPU1FRwc9+9rNgDkdEpEO+YOabJTPuq8GsiISWNdwD2Lx5M7/5zW/4+uuv2b17d5evb2xspLGx0f99VVVVMIcnIoOE21tj5gtjABaz7zklMxEJjbAW/zc2NnLllVfypz/9iREjRnTrPffffz/Jycn+r9zc3CCPUkQizZNf7eE/y/e1eMw/Y9ZsyqxpxkzBTERCo8fB7J577sFkMnX6tWrVqm5d684772TSpElcc8013f78O++8k8rKSv9XQUFBT/8IIjKI1dmd/Pbdzdz95kbq7S7/450tZbqUy0QkRHq8lHnLLbdwxRVXdPqaUaNGdetan332GRs2bOCVV14BmnY+ZWRk8Mtf/pJ77723zXtsNhs2m61ngxYR8XK4PHg84AFKahrJTYsD2l/K9IU0zZiJSKj0OJhlZGSQkZERkA9/9dVXqa+v93+/cuVKrr32WhYvXszYsWMD8hkiIh0pq7U3BTNv+LK0qDHzLmWqxkxEQiSoxf/5+fmUlZWRn5+Py+Vi7dq1AOTl5ZGQkNAmfJWUlAAwadKkLvueiYj0VVmt3X/fF8xMzZcyzaoxE5HQCmowu/vuu3n66af938+aNQuARYsWceKJJwbzo0VEulTaTjBrr/jf5UZEJCSCuitz4cKFeDyeNl8dhbITTzwRj8ej2TIRCYmy2qbWO77VyhZLmdqVKSIhprMyRSRiNTpdfLm9mAaHq93n25sxMzULZiYV/4tIiCmYiUjEem55Pt9+cgWPfdl+8+qymqZg5jt2qXm7DIuOZBKREFMwE5GItau4BjDaYrSnefG/b1LM0k6NmSbMRCRUFMxEJGI1D17tabf4v3kfM82YiUiIKZiJSMQqrekqmDXNpPnCV4t2GaoxE5EQUzATkYhVUtv+EqZP8xozdztLmdqVKSKhpmAmIhGrqxmzWrvLv2PT08lSplYyRSRUFMxEJCLZnW4q6x1tn2gVsnx1aO0vZarGTERCS8FMRCJSeV3ns2U+vmDWboNZ709Ij5YyRSREFMxEJCJ11CKjNd/OzPaWMn3NZl0KZiISIgpmIhKRuqov8/Edy+Rq5xBzi87KFJEQC+oh5iIi4VLayY5MM27mmzdxqeULFry7Eap/hHvod4GOGsxqxkxEQkPBTEQiUoczZm4Xr0ffzQyz95gmJ7DuRdxDvgO03pVp3Kr4X0RCRUuZIhKRSjoIZubyXcww78busfCy83jjwbLdeBzGDJu5nV2Z6mMmIqGiYCYiEams2VJm81hlKdoAwAbPGG533kidOR48LmKr9gBNvcugaVlTE2YiEioKZiISkTpayjQXbwJgi3sEYCLfMgKAuKqdxvOmtjVmmjETkVBRMBORiFTSwQHmlsPeYOYZCcBOTw4AcZU7gNZLmcatasxEJFQUzEQkIpV20MfMcngjAJvdRjDb7MwGIKFqF9Byxsy3lKkJMxEJFQUzEYlI7S5l1hRjri3C7TGx1ZMLwEa7L5h1vJSpGTMRCRUFMxGJOHV2J/Xew8lb8Bb+7/UMoZ4YAHa4jaXM+Jp9WHH6W2SAasxEJPQUzEQk4nTYw+yQsYy5xWMU/KfERVFIGq6oeMweJyNNRe32MVMwE5FQUTATkYjT4TmZRd5g5q0vS4uPBkzUJeUBMN60v2WNmZYyRSTEFMxEJOJ0NWO22bsjMz0+GoDy+NEAjDMdaLEr01f873QpmIlIaCiYiUjEafecTGcjlGwDWs+YQZFtFADjzPtbnJUZbzNOratzuHRepoiEhIKZiEScdo9jKt4KbifumBQKSQMgPcEGQIHVCGp5pgOYTG2DmcvtodHpDvKoRUQUzEQkApW111zWu4zpypwCGOHLt5S5E2Nn5hhTIVaadnPGRVn892sanUEarYhIEwUzEYk4vuaycdFNwcpX+O/Kmup/yLeUuacxFaclBpvJyRBXof95s9nkv0atgpmIhICCmYhEnFLvjFl6QnTTgwdWA+DOmuJ/yLeUWVRjpyJ+DADZjvwW1/ItZ9Y2ttMXTUQkwBTMRCTi+GrM0uON4JXasB8KlgMmHCOP878uNzUWgILyesrjxwIwu/qLFteK982Y2TVjJiLBp2AmIhHHt5SZ4Z0xm1v6lvHEuNPwJOX4XzciLQ6A4upG1g65GLfHxBHVn8D2j/yv8c2YqcZMREJBwUxEIorb7fEX/6fFRxOFk1ll7xpPzvlui9cmx0aRFGMEr68dY3jSdYbxxDs/gYYqoFnLDC1likgIKJiJSESpanDg9HbqT4u3cbp5FQnOCkjMhnEL2rx+ZHo8AHtLavmz8zJKooZB1QH45NcAJPhrzDRjJiLBp2AmIhHFV1+WGGMl2mrmSsunxhOzvgUWa5vX+5Yz95bW0oCNV4ffYTyx6kko3+vflamlTBEJBQUzEYkoh6sbAMhIsJFSX8Cxlk24McHsb7V5rclkYkS6Ecx8gW5P4hzImmy8oGy3ZsxEJKQUzEQkouwqrgVgdEY840o/AWBn4lGQMgKA1icr+WbMfEwmE8RnGt/UljS1y7CrxkxEgk/BTEQiys6iagDyshJIajgEwP74KR2+vnUws5hpP5hpxkxEQkDBTEQiys7iGsAIZvGOUgCqrWkdvr51MDObTBCfYXxTW9zUx0zBTERCQMFMRCLKzsPNgpndCGa1nQSz7OQYrOamg8vbBDP1MROREFIwE5GIUdXgoKjKaC6bl5VAnG/GLKrjYGa1mMnxngAAvmDWtJTpK/6vU42ZiISAgpmIRAzfbNmQJBtJNivxjjIAaqzpnb4vt9lyptlEs2BWrHYZIhJSCmYiEjGaL2PSWIXVbcye1XQyYwYt68zM5mYzZnUlapchIiGlYCYiEcMfzDIToOYwAFWeWBzm6E7fNzK9+YyZCeK8M2zNdmVqKVNEQkHBTEQihj+YDUmEmiIAij0pXb5vREdLmfYaEswOQEuZIhIaCmYiEjHamzErIbnL941Ii/ffN5tMYEsEiw2ARHc5YCxlelp3pxURCTAFMxGJCA0OFwXldQCMG9IUzLozY5ab1mxXptkEzXZmxjkrAHC6PTQ63YEdtIhIKwpmIhIRdhXX4PFASlwU6fHRzZYyu54xS4yJIi3eqEPztzTz9jKLtZf7X6c6MxEJNgUzEYkIzZcxTSZTsxmzjoOZqdl9X52Z2eR91BvMLHUlxEQZPyq1M1NEgk3BTEQigi+YjRuSYDzgmzEjpVvv9+3MjLJ4fyw262WWoO7/IhIi1nAPQEQkEHzBbGymN5jVNs2Y5TR7XUfl+9cdOxqX28M507ONB1ody1RSY6fOrmAmIsGlYCYiEaFFc1lo2pXZKph1ZHpOCv+6anbTA/4ms6XER/tmzFRjJiLBpaVMERnwXG4Pe0trAe+Mmdvdo12Z7YprPmNmHMukGjMRCTYFMxEZ8IqqGnC4PFjNJoalxEJ9GXiM2a1Sknp30WY1ZvGqMROREFEwE5EBr6DM6F82LCUWi7lpR2adNQVnbys2/DVmzY5lUjATkSBTMBORAa+gvB5o1ijWuyOzLrrzw8s75Z8xKyHe1y5DfcxEJMgUzERkwNvv7fifm+o989I3YxZlHEbeq5OUfDNmrkZSo+zGZTVjJiJBpmAmIgNeQZkxY5aT2nLGrDY6vfcXjYqFaGOHZ5apyriegpmIBJmCmYgMeL4zMnPTfDNmvqXMPgQz8M+apfmDmZYyRSS4FMxEZMA7UO6bMfMGs9pioGkps9e8dWap7grjspoxE5EgUzATkQHN4XJTWNl+8X+fljLBH8ySPJXG9dT5X0SCTMFMRAa0gxX1uD1gs5rJTLAZD7Yq/u+IydTp0/6lzERXBaAZMxEJPgUzERnQ9pc3Ff6bfEkrUDNm3u7/8Y5y43qqMRORIFMwE5EBzddc1l/473JAXRkAdVFt+5h5etI7w7uUGWs3rqd2GSISbApmIjKg+XZk+ltl1JYAHjBZaIhK7tvFvcHM5g1mqjETkWBTMBORAc23lNnUXNZYxiQ+E4/J0reLe2vMohq9M3BayhSRIFMwE5EBrc1SZl2pcRvXx/oy8AczS71xTbvLjd3p7vt1RUQ6oGAmIgNaQXmrrv/1RqE+cWmY6GrbZRdijRo1U305YNSmaWemiASTgpmIDFgNDhfF1Y1As6XMhgrjNqaP9WUAsSkAmNwOUqwOQHVmIhJcCmYiMmD56ssSbFZS4qKMB30zZrGp/tf15gxzAKLiwGxcd2i0EQDVMkNEgknBTEQGrOY7Mv09zOorjFvvbFefmEz+6wyJNkKgWmaISDApmInIgLW/zBfM4poe9Aez1LZv6I2YFAAyrUYwU42ZiASTgpmIDFj+Vhm+MzKhWY1ZSmA+xDtjlmkxPqtONWYiEkQKZiIyYDUtZTafMWtbY9Yn3oCXajE+q0Y1ZiISRApmIjJg7Ss1wtLItPaWMlO6fL+py1PMm66TZq4FtJQpIsGlYCYiA5LH4yHfF8zSgz9jlmzyzZgpmIlI8CiYiciAVFZrp9obknKbz5h1UWPW49YZ3hmzJIwZM9WYiUgwKZiJyIC0z7sjMzs5hpgo75mYzkZwGI8HesYs0VMDqI+ZiASXgpmIDEi+ZcwR7dWXYQJbUmA+yDtjFu+uBrSUKSLBpWAmIgPS3lJjabH9+rIUMAfox5t3xizWZcyY1ds1YyYiwaNgJiIDUlPhf3zTg63qy7qz6bJL3hkzm7MKgEanOwAXFRFpn4KZiAxIvhqzoO7IBH/I8wUzu0vBTESCR8FMRAakfb6lzLRmM2Yd9DDz9PoU86ZrRTmqAA+NDi1likjwKJiJyIBT0+ikpMYOwIgQzZhZ3A5isGvGTESCSsFMRAYcX31ZalwUybFRTU8E+pxMgOh4MFsBSKYWu2rMRCSIFMxEZMDJLzOWMUc0L/yH4MyYmUzNuv8rmIlIcCmYiciAs9c7Yzaq+TImNKsxC2AwA3+dWTK12pUpIkGlYCYiA067h5dDyz5mgaQZMxEJEQUzERlwfEuZI1svZfpqzAI9YxaTDHhrzFT8LyJBpGAmIgPO3pJ2ephB04xZIIv/oWkpUzNmIhJkCmYiMqDYnW4KK+uBVq0yoFs1Zr3qadZsKbPRqT5mIhI8CmYiMqDsL6/D7YG4aAuZCbamJzye4NWYea+XRC0Olwe3uy8da0VEOqZgJiIDiq/wf0RaHKbmh2Haa8Djnc0KeI1ZCmDMmIGOZRKR4FEwE5EBpaC8KZi14Jsts9ggKhaAQJxhDrRolwE6yFxEgkfBTEQGlMLKBgCGpcS2fCJYPcyg7YyZgpmIBImCmYgMKEXeYDYkKablE8GqL2t2zRQtZYpIkCmYiciAcqjKCGbZya2CWac9zPpYrO+bMUMzZiISXApmIjKgHOpqxizQPcygaVemyVdjppYZIhIcCmYiMmB4PB7/jNnQ1jNmIagxs+HAhl0zZiISNEENZvfddx/z588nLi6OlJSUdl9jMpnafD3yyCPBHJaIDFDVjU7q7MZs1dBQ1pjZEsFkAbzHMimYiUiQWIN5cbvdzqWXXsq8efN44oknOnzdU089xRlnnOH/Pjk5OZjDEpEBylf4nxRjJTba0vLJHp6TaepJLw2TyTgvs75MxzKJSFAFNZjde++9ACxcuLDT16WkpDB06NBgDkVEIkCHy5gQ3BozMGbi6stIpkZ9zEQkaPpFjdktt9xCRkYGRxxxBI888ghud8c/9BobG6mqqmrxJSKDQ4eF/xDcGjNodV6mgpmIBEdQZ8y647e//S2nnHIKsbGxfPrpp/z0pz+lpKSEX/3qV+2+/v777/fPxInI4FLkmzFrN5h1r8bM09vWGc26/6uPmYgES49nzO655552C/abf61atarb1/vVr37FvHnzmDlzJj/96U/5zW9+w5/+9KcOX3/nnXdSWVnp/yooKOjpH0FEBqhOlzJ7WGPWY81mzFRjJiLB0uMZs1tuuYUrrrii09eMGjWqt+Ph6KOPpqqqiqKiIoYMGdLmeZvNhs1m6/X1RWTgOlTZCHSxlBnMGjN8S5nqYyYiwdHjYJaRkUFGRkYwxgLAmjVriImJ6bC9hogMXh0uZToaoNFbbxrf9POpRzsvu9Ks+79mzEQkWIJaY5afn09ZWRn5+fm4XC7Wrl0LQF5eHgkJCbz99tscOnSIefPmERsby6JFi/jlL3/JDTfcoFkxEWmjw6XM2sPGrcVmtLUIhmYzZsUKZiISJEENZnfffTdPP/20//tZs2YBsGjRIk488USioqJ46KGHuO2223C73YwZM4bf/OY33HzzzcEclogMQA6Xm5KaDpYya7zBLGFIgKfJmonPBCCDSg4omIlIkAQ1mC1cuLDTHmZnnHFGi8ayIiIdOVzdiMcDURYT6fHRLZ+sKTJuE7Lafa+nj2eYA5Bo9FrMMlWoXYaIBE2/6GMmItIVXw+zrMQYzOZWs2L+YNZ2w1DAJGYDMMRUrnYZIhI0CmYiMiD4Cv+HJLVTf+pfymx/xiwgvDNmqaYaXPaG4H2OiAxqCmYiMiD4Zsza7WEWihmzmBScJmMJNbqhOHifIyKDmoKZiAwITa0yYts+2YsZsx5vETCZqLMZGwBiGw739N0iIt2iYCYiA0JTq4z2ljJDMGMG1McYwSyuUTNmIhIcCmYiMiB0eoB5iIJZQ4wxI5dgVzATkeBQMBORAaHDrv8eT2iK/4HGWG8wc5QG9XNEZPBSMBORfs/j8XTc9b+xCpzeXZLdCWZ96Gnm8AazJEdJ7y8iItIJBTMR6fcq6x00OIzeYR12/bclQ1TLjQGmAJ8C4Io3lkqTXZoxE5HgUDATkX6vqMo4iik5NoqYKEvLJ7vo+h9I7nhvLzMFMxEJEgUzEen3Oqwvg5AV/gO4ExXMRCS4FMxEpN/zBbOscHX99zJ5g1kiteCoD/rnicjgo2AmIv3e4WpjKbO3rTICcog5YIlNpt7jPUC9+lBgLioi0oyCmYj0e4fDfU6mV3SUhSJPqvGNgpmIBIGCmYj0e77i/3A2lwWItpgpwhfMCoP+eSIy+CiYiUi/V1TtrTFLDG8ws0VZOOxJMb7RjJmIBIGCmYj0e4f9M2aBW8rsTY+zaIvZv5TpUTATkSBQMBORfs3t9nC4uoNzMt0uqPWeWxmKpUyr2T9j5q7SUqaIBJ6CmYj0a+V1dhwuY1tlZmKrGbO6UvC4wWSG+Iygj8VmbT5jpmAmIoGnYCYi/Zqv8D89PpooS6sfWb76srgMMLc6ESAIoi1mDqNdmSISPApmItKv+Qv/A7Qjsy8tzcxmE6WmNON+jYKZiASegpmI9Gv9pYeZT7klHQCzvRrstSH7XBEZHBTMRKRf8+/IDHOrDB+nNZ5ajzckajlTRAJMwUxE+rWi6v41Y2azqvu/iASPgpmI9Gu+4v9A1Zj1VbS1+QYA7cwUkcBSMBORfq2pxqydYFa+17hNHt7pNTx9KvlvKdpq5oAno+Xni4gEiIKZiPRrRR11/fd4oGSncT99XMjGE20xs9udbXxTuitknysig4OCmYj0Wy63h+KaDg4wry2GxkrABGljQjYmW5SZvZ6hxjdlCmYiElgKZiLSb5XWNuJyezCbjAazLZTsMG5TRkBUO8ucQRJtMbPHF8xKd4bsc0VkcFAwE5F+y9cqIyPBhrV11/9SbzDLCN0yJhg1Zv4Zs7pSqCsL6eeLSGRTMBORfquos8L/0r7Vl5l6OSab1UwdMdTFeFt0lO3u5ZVERNpSMBORfqvDwn9oVvg/NoQjMvqYAVTGjjAe0AYAEQkgBTMRCbvDVQ2c+cBinlqyp+Xj3uayme11/Q/jUiZARexI7zhUZyYigaNgJiJht2RXCVsKq/jP8n0tHu9wxszlaOohFsJWGWAU/wOUxeQaDyiYiUgAKZiJSNiV1tgB2Fdah93p9j/eYXPZ8r3gdkJUPCQN6/C6pnYKyTx97DVrizJ+bBZHK5iJSOApmIlI2JXVGsHM6fawr7TW/3iH52T6WmWkj20/fQWRb8asKNp72kDZ7r6nPRERLwUzEQk734wZwI7DNQC43R72ldQBMCwlttUbfMEsLyTja85XY3bYkg0mM9hrms7sFBHpIwUzEQm70tpmwazICGZ7S2upbnRis5rJy0xo9Qbv8mGIC/+hKZg1uC2Qog0AIhJYCmYiEnZltY3++zsOVwOw4UAlAJOHJbVtLtvDMzIDudLoa5fR6HQ3zdgpmIlIgCiYiUjYlTWbMdvpXcrc6A1m04Ynt32Dv1VG+JYy7U53Uw819TITkQBRMBORsGu+lLm7pBany836/R0Es/oK4wBzgLTQNpeF1sHMN2OmYCYigaFgJiJh1eh0Ud3gBMBsMgLPvrI6Nh2sAmBaTqtg5ls2TBgKMUmhHCoANu+yqt3VfMZMS5kiEhgKZiISVuW1DgAsZhMThhpB65PNRdQ0OomJaqfwvyQ8Hf99fH3MGp2uphmz8j3gdoVlPCISWRTMRCSsSr2F/6lx0UwYYoSw11YfAGBydjuF/0UbjdvMiX363N62P/P1MbM73ZCUA9ZYcNk1ayYiAaFgJiJh5Sv8z0iIZtyQRAC2FRk7M9st/D+03rjNnh6S8bXWosbMbIbhc4wn8peHZTwiElkUzEQkrHzNZdPio8nLarlsObV1MPN44NAG4/7QaaEYXhu+YNboOzpqxNHGrYKZiASAgpmIhJVvR2ZafDTjWgWz6TkpLV9cuR/qy8FshazJIRphS74+Zv4zPUfMM27zl4VlPCISWRTMRCSsfM1l0+OjGZEW56/hiokyMzYzvuWLfcuYmRPB2ur8zHaYaFtI5qFv3WbbzJjlHmkczVS+B6oP9enaIiIKZiISVmX+GTMbVouZ0RlGGGu38N+/jBme+jJoVvzv8gazmCQYMsW4r1kzEekjBTMRCStfjVl6QjQAed6dme0W/hd6Z8zCVF8GrYr/ffzLmaozE5G+UTATkbDy1ZilxxvB7NpjRnH0mDS+PX9U2xeHeUcmgM3arI+Zj28DwL6lYRiRiEQSa7gHICKDW1mz4n+AOSPTeOGGeW1fWFcGlQXG/TDOmNk6mzEr2ggNVWE5kUBEIoNmzEQkrEprvMX/3qXMDvnqy1JHQUw7y5yd6Fu5f0u+pUy3B5y+OrOkYZAyEjxu2L8ygJ8mIoONgpmIhI3d6abKe05menwXuyx9y5hhLPyHpmAGzXZmgurMRCQgFMxEJGzK64xlTIvZRHJsVOcvLuwnwazZTtGWy5m+RrPamSkivadgJiJh49uRmRoXhdncxeGV/aDwH8BqMWPxjtXfMgOajmbyneUpItILCmYiEjatC/875KiHku3G/QDNmLXXfLa7Whxk7pOeZ9zWl0NtaV+GJiKDmIKZiIRNqbfrf5fBrHirUVgflw6JQ0Mwss7ZotppmREdB8kjjPu+ECki0kMKZiISNr4Zs/SELgr/i7cZt1mTwdT7ma5A8c2YtSj+B8gYZ9wqmIlILymYiUjY+Lv+d2fGDCBzQpBH1D3tdv8HyBhv3CqYiUgvKZiJSNiUdrfGzDdjljmxR9dvb3LNE4CmZh0HM2+dWcmOvn+IiAxKCmYiEjZl3hqzLmfMDm8xbvvJjJnNagHaW8rUjJmI9I2CmYiETdOuzE5qzBz1UL7XuN/DGbNg6XIps2IfOBtDPCoRiQQKZiISNv4as86OYyrZAXggNhXiM0MzsC7YfO0yXK2CWcIQsCUZO0jLdodhZCIy0CmYiUjY+GrMOl3K9NeXTeoXOzKhkxkzk0k7M0WkTxTMRCQs7E43lfUOoIvi/wDsyAxEwX9zCTYrAMXV7SxXqs5MRPpAwUxEwmL5bqM7fkZCNKlx3Qlm/aO+DGDWiBQAvt7TTof/dO3MFJHeUzATkbB4b0MhAAumDO38nMx+1sMMYP7YDAC+3l2Gs3WdmX/GTMFMRHpOwUxEQs7hcvPhpkMAnD0tu+MXOhubiuj70YzZ5GFJJMVYqW50svFgVcsnmwezQK+hikjEUzATkZBbvruU8joHafHRHDk6reMXlu40djjGJAf+jMw+7COwmE0cPSYdgKW7Slo+mTYaTBawV0P1oT4MUEQGIwUzEQm59zYYgWXBlKFYLZ38GGpeX9ZPdmT6HJNnLGcu29Wqzsxqg9RRxn1tABCRHlIwE5GQcnZ3GROatcroP/VlPvPHGjNmK/eW0eh0tXxSLTNEpJcUzEQkpL7eU0ZZrZ3UuCiOHtPJMiY0O4qp/9SX+eRlJZCRYKPB4WZtfkXLJ311Zoc3h3xcIjKwKZiJSEi9s75pN2any5geDxzaYNzv5YxZMBc/TSaTf9ZsaevlzJwjjNv85UEcgYhEIgUzEQmZvSW1vPrNfgDOmzGs8xcfWg/le8Biawo6ARDIfZK+YNamzmzkfOP28GaobafXmYhIBxTMRCRkfvvOZuwuN8eNy2CeN9R0aN2Lxu2EM41dmf2Qr5/ZmoJy6uzOpifiMyDDO8uXvywMIxORgUrBTERC4tMtRXy69TBWs4lfnzsFU2e7LF1O2PCycX/GFaEZYC/kpsUyJMmGw+VhS2GrfmajjjFu9y0N/cBEZMBSMBORoGtwuPjNO0Yh/HXHjiYvK6HzN+z5HGoPQ2wa5J0a/AH2kslkYsowYzZvU+tGsyN9wWxJiEclIgOZgpmIBN1/v85nX2kdWYk2fnjKuK7f4FvGnHoxWKL6/PmegFaWtTQ5OwmATQdaBbMR84zbQ+uhodVzIiIdUDATkaDyeDy8sDIfgB+eMo4Em7XzNzTWwNZ3jPv9eBnTZ8owbzArrGz5RPJwo9Gsxw0FK0I/MBEZkBTMRCSo1u+vZHtRDTarmfNndrETE4xQ5qiDtLEwfE7wB9hHvqXM7YdqcLQ+0FzLmSLSQwpmIhJUL39TAMAZU4eSFNONZclt7xm30y7td8cwtSc3LZbEGCt2l5sdRTUtn/S1zdAGABHpJgUzEQmaBoeLt9YeBOCSOTldv8Hjgb3e2aWxJwVxZIFrPmsymZrqzA62Ws70BbMD34CjPkCfKCKRTMFMRILm481FVDU4GZYc4+/51anibVBXAtZYGDY7+AMMkA53ZqaOhsRscDtg/8owjExEBhoFMxEJmpe9Xf4vnpODxdyNOaq9i43b3CPBGh3EkQWWbwPA5tbBzGSCUcca9/d8GeJRichApGAmIkFRVNXAVzuKgW4uY0JTkbwvzPRRqErUpgz3BrPCKtzuVq05Rp9g3O7+IjSDEZEBTcFMRIJiTX45bo/R52tkenzXb/B4YO9Xxv0ABbNQGZuZQLTVTE2jk/yyupZPjvEGswPfqJ+ZiHRJwUxEgmK7d4fixOzE7r2hZAfUFoM1JqhtMjyewDebjbKYmTjU+HO2qTNLGWHUmnlcapshIl1SMBORoNheVA3A+CHdDGa++rKcI8BqC9KogsffaLb1zkxomjXTcqaIdEHBTESCYudhY8Zs/JAuzsX0CXB9WahN7mhnJjTVme1RMBORzimYiUjAOV1udhfXAjAuqxszZs37l/m65Q8w04cbwWz1vvK2JwCMPt64PbwZag6HeGQiMpAomIlIwO0trcPuchMbZWF4SmzXbyjdBTWHwBINOXMDP6DgnWHuN3V4MqlxUVQ3OlmTX9HyyfgMGDLNuK+2GSLSCQUzEQm4Hd76snFDEjB3p3/Z7kXGbc4RENWNINcPWcwmjhuXCcAX29uZFfPXmX0eukGJyICjYCYiAbfDW1/WrWVMgJ2fGrd5pwRpRKFxwngjmH2+rbjtk83rzIKwM1REIoOCmYgEXNOOzG4U/jsbm5b38k4N4qiC73hvMNt0sIrD1Q0tnxw531iqrciHku1hGJ2IDAQKZiIScDu8PczGdSeY7VsKjlpIGAJDpwd5ZE2CcSpAZqKNad5NAF9uL2n5pC2hadZs6zuB/3ARiQhBDWb33Xcf8+fPJy4ujpSUlA5ft3DhQqZPn05MTAxDhw7llltuCeawRCSIHC43u0t6sJS58xPjNu/U0J2hFES+5cwvtreznDnxbON267shHJGIDCRBDWZ2u51LL72Um266qcPX/PWvf+WXv/wlP//5z9m0aROffvopCxYsCOawRCSI9pXW4XB5iIvu5o7MHR8btwN8GdPnxAlGMFu8oxhX63MzJ5wFmIzjmaoKQz84Een3rMG8+L333gsYM2LtKS8v51e/+hVvv/02p5zSVPQ7ZcqUYA5LRILIvyMzqxs7MivyoWQbmMww9qSAj8VE6GfgZuamkBRjpaLOwbr9Fcwekdr0ZOIQY+fp/hWw7T044rqQj09E+rew1ph9/PHHuN1uDhw4wKRJk8jJyeGyyy6joKAgnMMSkT7wnZGZ15NlzJwjITa189cOEFaL2d82Y9HWdtpmTDzLuNVypoi0I6zBbPfu3bjdbn7/+9/z97//nVdeeYWysjJOO+007HZ7u+9pbGykqqqqxZeI9B/bD/dgR+YObzAbFxnLmD6nTR4CwFvrDrY9NH3iOcbtni+hoZ1zNUVkUOtxMLvnnnswmUydfq1atapb13K73TgcDv7xj3+wYMECjj76aJ5//nl27NjBokWL2n3P/fffT3Jysv8rNze3p38EEQminUW+MzK7mDFz2pvOjsw7LcijahKKFmKnTxlCXLSFfaV1rCmoaPlkxjjIGA9uR1N9nYiIV4+D2S233MKWLVs6/Zo6dWq3rpWdnQ3A5MmT/Y9lZmaSkZFBfn5+u++58847qays9H9p2VOk/3C7PewpMc7IzMvqYsYsfxnYayA+M6RtMkIhLtrKgilDAXhjzYG2L+hnuzNX7i3jlW/2h3sYIkIviv8zMjLIyMgIyIcfc4xxWPG2bdvIyckBoKysjJKSEkaOHNnue2w2GzabLSCfLyKBVVzTiN3lxmyCockxnb94Z7PdmObIa6l4wazhvL7mAG+vO8hd50wmytLszzjhbPjqb8aJBy4HWKLCNk6328NN//mGkho7M3NTug7UIhJUQf1pmJ+fz9q1a8nPz8flcrF27VrWrl1LTY13qWP8eM4//3x+/OMfs3TpUjZu3Mh3vvMdJk6cyEknBX6HlogE1/7yegCGJsW0DCLt2dGsf1mQheMApGPGppORYKO8zsGXrXuaDZ8NcRnQWAn5y8Mwuia7S2opqTFqegvK68I6FhEJcjC7++67mTVrFr/+9a+pqalh1qxZzJo1q0UN2jPPPMNRRx3F2WefzQknnEBUVBQffPABUVHh+w1SRHrnQIURzIandtG/rHI/FG/xtsk4OQQjCz2rxcx5M4YB8Hrr5UyzBcadbtzf/kGIR9bS6n3l/vvFVY1hHImIQJCD2cKFC/F4PG2+TjzxRP9rkpKSeOKJJygvL6e0tJTXXntNBf0iA9RBXzDrqrGsr03G8LkQlxbkUYXPRbOHA/Dx5iKqGxwtn5xwhnEb5mD2TbNgVlTV0MkrRSQUIq+wQ0TC5kB5N2fMfLsRx4VuN2Y4TBmWRF5WAo1ON2+sPdjyyTEngTkKSndCyc7wDBD4Jr8pmB2u1oyZSLgpmIlIwPiXMlPiOn6R0w67fW0ywte/LBSnAphMJq4+agQAC5fsadnTLCYJRhkboNjxYdDH0p6KOjs7D9f4v9eMmUj4KZiJSMB0a8as4GuwVxttMrJnhmZgYXTJnBwSbFZ2FdeyeEdJyyfHe5czt70f+oEBa/IrWnyvGTOR8FMwE5GA8Hg8zWbMOglmvjYZY0+JyDYZrSXGRHHJHKMd0MKle1s+OX6BcZu/DOorQjougNXeZcyxmfEAHNaMmUjYRf5PRREJiap6JzWNTqCTYObxwJa3jfshWMY0hf4M83Z9d/4oTCb4bOthfwNeANLGeE8BcDZtiAghX+H/mVONZt/FNY243eFoLiIiPgpmIhIQ+yuMHljp8dHERlvaf9HexVC2G6ITmw7zHgRGZcRz0oQsAJ5uPWs26VzjdtWTIR2T0+Vmrfe4qDOmGqcUOFweyuvaP6dYREJDwUxEAqJb9WXfLDRup18G0fHBH1Q/8t35owB4eVUBxc1rueZeB2Yr7FsCB1aHbDxbD1VTZ3eRaLMyOTuJtPhoQHVmIuGmYCYiAdFlfVltCWx+y7g/57uhGVQ/cty4DKbnJFNrd/HXj7c3PZE8HKZebNxf9mDIxuOrL5s1MhWz2URWonHUnXZmioSXgpmIBIR/xqyjYLb2v+B2wLBZkB1Zh5Z3h8lk4ldnTwbgxZX5bCmsanpy3i3G7abXoaIgJONZv78SgFm5KQBkJRlnm2rGTCS8FMxEJCA6PY7J42laxhyEs2U+R45O4+xp2bg98Lt3Nzf1NcueDqOPB48Lvn4kJGPxzYzlphk954Z4Z8y0M1MkvBTMRCQgOl3K3PsVlO2C6ISmZbtB6udnTiTaYmbJzlI+3XK46QnfrNk3T0NDVftvDiDfweWZ3kCWleQNZpoxEwkrBTMRCQjfUuaw9oLZmmeN26kXgy0xhKMytOi4H2a5aXFcd9xoAO59ZxP1dpfxRN5pkD7OaL676fWgj8O3ASEjwSj6H+JdylSNmUh4KZiJSJ/V212U1hozMDmtlzIbKmHzm8b92d8J8cj6p5tPymNYcgwFZfX8/RPvRgCzGWZdY9xf+1xQP9/l9lBWawQz/4xZombMRPoDBTMR6TPfMmZ8tIXk2KiWT258FZwNkDkRhs8Ow+j6nwSbld9eMBWAx7/aw8YDRiE+0y8Hk9k4tiqIB5uX1dpxe4wGvGlxxoyZv/i/SsFMJJwUzESkz5oX/ptat9tf4539mXVN/2nFT/iHcsqkIZw9PRuX28PPX1uP0+WGpOymExGCOGvmW8ZMj4/GajH+GvDNmBVXN/arpV+RwUbBTET6rMNWGYe3woFVYLIYs0HSwq/PnUxSjJWNB6p47ut848GZVxm3614Atyson1tS46svs/kf8y1p2l1uKuocQflcEemagpmI9NkB73FMbVplrP2PcTv+DEjICvGo+r+sxBjuOGMiAP/8bAe1jU6YcBbEpED1Qdi9KCif65sx84UxAJvVQmqcsQxdVK0NACLhomAmIl1qdLqoauh4FqVpxiyu6UGXw5j1gaaidmnj8iNyGZkeR0mNnaeW7AGrDaZdajy59r9B+UzfjFlmsxkzaNqZqTozkfBRMBORLv3g2W+Y9/tPKSira/f5vaXG4yPTmwWzHR9DbTHEZ8K400IxzAEpymLmttPGA/Dol7upqLPDrKuNJ7e8bRxlFWD+VhmJLYNZpo5lEgk7BTMR6VRBWR2LthVTa3fx0eaiNs97PB52F9cAMCaz2cHka7zLmDOuAEtUm/dJk3OnD2Pi0ESqG5w88sVuyJ5pHF3lssPqZwL+ecVdzZipZYZI2CiYiUinPth4yH//qx3FbZ4vrbVT1eDEZIJR6d5gVnMYdnxo3J+pZcyumM0mfnb6BAAWLt3DlkPVcMT3jSdXPdnlJoD95XX+5cnu8Bf/J0a3eDxLxzKJhJ2CmYh06t0Nhf77X+8pw+50t3h+d3EtAMOSY4mJshgPrn8R3E4YPheyJoZsrN3RXztBnDIpi/lj02lwuLnm8a/ZNeQ0iE2DygLY/kGH73t9zX6O/+Mi5v7uE47/4yJufXFth0vOPv7i/4SYFo9rxkwk/BTMRKRDByrqWVtQgckEiTFW6uwu1uSXt3hNm2VMj6dpGVNF/91mMpl4+Jo5TB2eRGmtnSufWk/lpCuMJ1c81u573t9QyE9fWofbGzbzy+p4fc0B7nt3S6ef5Tsns6MZM9WYiYSPgpmIdMi3jHnEqDROmmC0u/hqZ8ti9N0lxozZ2MwE44EDq6F4K1hjYepFoRtsBEiOjeLZa49i4tBEDlc38p0N0/BgMtpmlOxo8dpF2w7zoxfW4PbApXNyWHf36Tx4lXGywkebD7G/vP1ZM4fLTZn3+KzWNWZZ/vMyNWMmEi4KZiKDnMfj4dqFKzn7H4t5bfV+HK6mpcr3vcuYZ00dyrHjMgBYvKNVMGs9Y+brXTb5PIhJDvLou6efrl62KzU+mmevO4rxQxJYW53MF3iPsfr6UcD49/Xssr18/+lVOFwezpmezR8unk5yXBRnT8/mmLx03B54dvm+dq/vC2UWs4nUuJYzZmO9/w4PVNRT2oOaNREJHAUzkUGuutHJZ1sPs+lgFbe9tI4T//Q5D3yyg8+2FrFqn7Fseea0bI7NM4LZ+v0VVDbrDO+bMRuTkWAcWL7hFeOJmVeH9g8SQTITbbx04zzmjkzlUfsZADhWPc2bX63h9lfWc9ebm3C6PZw7Yxh/u3wmFnPT+VLfnT8agBdWFFBvb7tpoPlxTGZzy3OpUuKiycsyZj6/2Vfe5r0iEnwKZiKDXF2j8Ze3yQQZCdEcqKjnb59s59qFqwCYOzKVIUkxDEuJZUxmPG4PLNtdChjLYvneHmZjMuNh1VPQWGUcWD7quPD8gSJESlw0/7n+KOLHn8Qadx5RHjuFH/yFV77Zj9kEd545kX9cMZMoS8sf4ydPzCI3LZbKegdvrD3Q5rr+HmatljF95o5MBRTMRMJFwUxkkKu1OwFItFn56n9P5o+XTOesaUNJjLECcNVRI/yvPc47a/bVTqNtRkFZHU63h9goC0PjgOUPGS885sdg7t8/XvrPceodi4my8Oi35+I+9mcAfCfqE44dbuGZa4/ixhPGtj0wHmOJ8jvzRgHw1JI9bQ4k9/cwS2w/mM3xBrNVCmYiYdG/f3KKSNDVNhrBLMFmJSbKwmVzc3no6jmsues0VvziFC6aneN/7bHjMgH4cnuJt7GssYw5OiMe84YXoaYIkobD1EtC/weJUBaziTmnXQFDphHrqec/U9f46/06cuncXOKiLWwvqmnRhw66MWM2Kg2ADfsraXAE5xB1EemYgpnIIFfrXcqMs1lbPG61mP279HzmjU0nNspCflkd3+wrZ3eJUfg/NiMGlvzD+6KbwdqyqFz6yGSC439q3P/6YWio6vTlybFRXHesUWt215sbKfcW/EOzczI7mDEblR5HRkI0dpebjQcqAzB4EekJBTORQc43YxbfKpi1J8Fm5dwZ2QD8d0W+f8bsdPMqKNsFMSkw+ztBG2tPtbfUN2BNOg8yxns3WLzU5ctvOTmPcVkJlNTYueftTf7H/c1lOwhmJpNJy5kiYaRgJjLI+WrM4qMt3Xr9lUcaNWfvri9kbUEF0Tg48YDRyoEjbwBbQlDGOeiZLTDzKuP+1ne7fLnNauHPl87AbII31x7kw03Gkqb/OKaEjmc15440ljNX7VUwEwk1BTORQc63lNmdGTOAmbkpTByaSKPTzdZD1VxveZfEmj0Qn2UsY0rwTDzHuN3zJdRXdPnyGbkp3HjCWAB+9cZG6uzOLmfMAOaMMmbMVueXt9k8ICLBpWAmMsjV9XDGzGQy+WfNckyH+ZH1deOJBfdBbEowhig+GeOM5Uy3E3Z+0q23/PiUceSmxVJc3ch/lu9rdk5mx8Fs6rBkbFYzZbV2f586EQkNBTORQa7GW2PWuvi/MxfMGo7NauIe69PEmBxGz7JplwZriAHlGVDnALTDN2u25e1uvTwmysKPTh4HwMOf76Kqwfj33dmMWbTVzIycFAC+0XKmSEgpmIkMcnXe7vAJPQhmybFR/HTMAU61rMGJFc7+i7FzUILPF8x2fgKO7h02fuGs4YzOiKfce2JDlMVEcmxUp++ZO8q3AaCs92MVkR5TMBMZ5Hy7MuO6uZTpc43lYwD2j70SMicEfFzSgWGzIDEb7DVGrVk3WC1mfnzKOP/3GQm2LneszhphBLN1BWqZIRJKCmYig1zzBrPdVlNM3L5PARi1oP8X/EdU/brZDBPOMu5vfafbbzt3xjD/IeUdNZdtbkaucQD99sPV/uVuEQk+BTORQa7Wu5QZF92DYLbhJaMAfdhsyJoUpJFJhyaebdxuew9cjs5f62Uxm7jjjImYTDA9J7nL12clxjA8JRaPBzWaFQkhBTORQa6pwWw3lzI9HljznHF/1tVBGpV0atRxEJcOtcWw7F/dftuCKUP58vaTuPvcyd16vW/WbG1BRW9GKSK9oGAmMsj5ZsziuztjVrgWDm8Ciw2mXhy8gQXZgD4VwBoNp//OuP/5H6B0V7ffmpsWh83avRDu25m5TsFMJGQUzEQGuZ4cyQQ0zZZNOgdiU4M0KunSjCth7MngbIC3fgRud+A/IjcFUDATCSUFM5FBrq4nS5n2WtjwsnF/ppYxw8pkgnP+DlFxsO8rWP10wD9i2vBkzCY4WNnA4aruteYQkb5RMBMZ5Gp6MmP22e+goQJSRsKYE4M6rkAYwIuV3ZM6Ek6+y7j/yT1QF9ieY/E2K+OHJAKqMxMJFQUzkUHM4/H4G8x2WWNWsAKWP2zcP+evxqHaEn5H3QhZU4zA/MX/Bfzy/jqz/RUBv7aItKVgJjKI2V1unG6jyVdcZ0uZjgZ482bAYyxh5p0amgFK18wW45xSgBWPQfH2gF6+qc5MLTNEQkHBTGQQq210+e93OmP25R+hZDskDGkKAdJ/jD0Jxp8JHhd89KuAXtrXMmNdQQVudyR16hXpnxTMRAYx347MmCgzFnMHFVmF6+Crvxv3z/7LgN+JGVGnADR3+u/AbIUdH8LW9wJ22QlDEomJMlPd6GR3SW3Arisi7VMwExnEau1dHMfkchhLmB4XTLkQJp0bwtFJj2TkwZE3Gvdf+R7s+CQgl7VazEwb3jRrJiLBpWAmMoj5ljI7PI5pyd/h0AaITYMz/xS6gUnvnPpr4xxNZwO8cGXAZs6mDDOC2ebCqoBcT0Q6pmAmMoh12ly2eBt88Ufj/pn/BwmZIRxZYEXq6mUbVhtc+jRMPh9cdnjpW1C4vs+XnTwsCYDNBxXMRIJNwUxkEKvzLmXGR7ezI/Oju4y/3MefAdMuDfHIpNes0XDxkzBugXHQ/Jd/7PMlJ2d7g1lhFZ6ILdIT6R8UzEQGsRrvUmabGbODa40icpMZFvze6DIvA4fFCqf9xri/5W04vKVPlxs/JJEoi4nKegcHKuoDMEAR6YiCmcgg5p8xa93D7EtvPdnUSyB9bIhHFRoRHzWzJsKk84z7i//Sp0tFW83kZRknAGg5UyS4FMxEBjH/cUzNi/+LNsHWdwATHP+z8AxMAsP372/jq1C6q0+Xar6cKSLBo2AmMojVtbeU6ZtdmXw+ZE4Iw6gCZ9CvwGbPMGrNPG746q99utQU7waATZoxEwkqBTORQczXxyzOV/xfsgM2vmbcP/72MI1KAsr373HdC1CR3+vLaGemSGgomIkMYm3aZSz+C+CBCWfD0KnhG5gETu4RMOZEY4fmkgd6fZlJ3qXMAxX1VNY5AjQ4EWlNwUxkEKu1e5cyoy1QtgfWv2Q8odqyyOKbNVv9LFQV9uoSybFR5KTGArCpUAeaiwSLgpnIINZixuyrvxpHL+WdCsNnh3lkElAjj4ER88DVCMv+1evLTNFypkjQKZiJBJjb7eGNNQf4xesbKK1pDPdwOuUr/k93HYa1zxsPHn9HGEckQWFqtsN21ZNQW9Kry0zO1tFMIsHWwQF5ItIba/LLufftzaz1HvYcY7Vw97mTwzuoTvjaZUzY+QS4HTD6eBhxVJhHFVyDtm/92FNg2Cw4uAaWPWicq9lD2gAgEnyaMRMJkHUFFVz88FLWFlRgNRt9Gt5adxCnyx3mkXXM0ljO/1n/zfAdzxkPaLYscplMTbVmKx6D+vIeX8K3lLnjcA0NDlcgRyciXgpmIgHyyZYi3B6YOzKVL+44idS4KEpqGlmyqzTcQ2vfrkUsrLuFy62fG9/P/yGMOjasQwoWne/oNf5MyJoM9mojnPVQdnIMafHRuNweLWeKBImCmUiA+JZ3zpmezfCUWM6dMQyAN9YcCOew2tdQBa9eTzqV7HAP59DFb8Lpv1NH1khnNsNxPzXuL38IGmt69HaTycTM3BQA1uZXBHZsIgIomIkEjK8j+pThRoH0BbOGA/DBxkP+3Y/9xtJ/Ql0Ju91DOcd+H9bR88I9IgmVKRdC2lhjKfObp3r89lm+YOatoxSRwFIwEwmA0ppGDlU1ADBxqHHY86zcFEalx1HvcPHx5qJwDq+l6kP+lgl/cl5OI9Ek2AbhPqDBOjlotsCxtxr3l/6zx7NmM0ekALCmoOc1aiLSNQUzkQDYUlgNwMj0OBJjogBj2cc3a/Z6f1rO/PwP4KjDMXQ277uPxGwCm1U/CgaV6ZdDci7UFMEfR8MTpxs1Z92oxZuRm4LJBAVl9ZT083YwIgORfhqLBMBmbyd03641nwtmGsFs8Y5iymrtIR9XG8XbYPUzxt15vwJMxNusmCK0tiwy/1QBYI2G8/4BidngskPB1/Dez2DHR12+NSkmirGZCYDqzESCQcFMJAB89WWTs1sGs1EZ8QxPicXtgb2lteEYWpO6MnjhaqO7//gzKcuYC0B89CBcxhQYezLctgV+tMaYQQP4+pFuvXWm6sxEgkbBTCQAfDsyJ7eaMQPISrIBcNhbgxYWjnp4/goo3QFJOXDO3/zNZeNtlvCNS8LLZIK0MXDSL8Fkhl2fGbOqXZjlrTNTMBMJPAUzkT5qcLjYVWwUUE8Zltzm+SGJMQAcrg5TPY7bDa/dYCxXxSTDNa9AUjZ19mbnZMrgljoSJpxl3P/60S5f7psxW1dQgdutHnEigaRgJtJHWw9V4/ZAenw0WYm2Ns83zZiFKZit+y9seQss0XDFfyFrEgA13nMytZQpABz1A+N23fNdngowYUgisVEWqhud/l9KRCQwFMxE+qj5MmZ7RfS+sHa4OgxLmQ1V8Mm9xv2T72rR2b9OS5nS3KhjIWsKOOqMWTNnx5tVrBYz03KM2eE12gAgElAKZiJ95NuR2V59GUCWdymzKBwzZov/DLWHjYaivhkRr1q7d8ZskC1l6nimDphMcNSNxv3P74ffZ8Mjx0L+8nZf7ms0u0Z1ZiIBpWAm0kcd7cj08S9lhrrGrHQXLHvIuL/g90aLhGZ8pxHEaSlTfKZfDhPPAVsyuJ1waAO8ej042s72+jYAfLWzWHVmIgGkYCbSBy63h63e5rKte5j5+GbMikO5lFm+F968GdwOGHsKjF/Q5iW1vuL/aC1lildUDFzxHPx8H/x4PSQOg8oCWPVEm5cePz6TpBgrBWX1LNp2OAyDFYlMCmYifbDtUDX1Dhdx0RZGZyS0+xrfjFlprR2Hyx3cATVUwvs/h3/OhfxlYI0xZsvaqX2rbRw8uzI1n9NDJpOxU/OkO43vv/yT8d9WM3HRVq44cgQATy3ZG9LhaYZOIpmCmUgfLNlZAsCRo9OwmNvvM58WF43VbMLjIbhH2Hg88Mp18PXD3pmyk+H6TyBrYrsvr/PtylTxv3RkxlWQMcHYpbnkgTZPf3veSMwm+GpnCduLqkMypC2FVRxx3yfc9+7mkHyeSKgpmIn0wVfeYHZsXkaHrzGbTWQkhKBlxo6PYOfHYI6Cq1+Fb70OQ6d1+PKaQTRj1h4d19QNFiuc+mvj/rKHoKqwxdM5qXEsmDIUgKeW7Onx5Q9W1Pt3NXeHx+Ph7jc3Ulpr56PNRT3+PJGBQMFMpJcanS5W7CkD4NhxHQczgCHB3gDgtMOHvzDuH30TjDu1y7fU2dXHTLphwlmQexQ46+H929s8/b1jRgPw2uoDlPfgPFiHy82ljyzjnH8uZuXesm695611B1m51+ixVljZoB22EpEUzER6aU1+BfUOFxkJ0UwYktjpazP93f+DtAFgxb+hdCfEZ8Lxbf/ybI9vWTUlLio4Y+oHIvVw9pAymeDsv4DZClvehk1vtHj6iFGpTBmWRKPTzSvf7O/2ZRdtPcyBinrcHvjl6xu6rL+sszu5/72t/u/tTjdlPQiCIgOFgplIL321w1jGPCYvo8sAENTu/zXF8MX/GfdP+TXEtL87tDmPx0N+WR0AI9LiAj8miSxDp8Gxtxr33/sZ1DXNcJlMJi6dkwPAF9uLu33Jl1YV+O9vL6rhia86Xwp9aNEuDlU1kJMaS6r3l4nCyjCePysSJApmIr3kqy87ppP6Mp+gdv//6JfQWAXZM2Dm1d16S0mNnTq7C7PJqBMS6dLxtxsbAWqL4f07jDNYvXz/D6zcW0aj09XlpQ5XNbBomxHibjpxLAAPfLKD/eV17b6+uLqRfy/eDcCvzp5MrveXiUMKZhKBFMxEeqGy3sH6/RVA54X/Pr5eZgGfMdvxCax/EUxmOPtvYO7e/9L5ZbUAZCfHEm3VjwHpBqsNzv8XYIINL8Oz50OlsXSZl5VARoKNRqe7W0c0vbJ6Py63hzkjU7ljwQSOGp1GvcPFb99pf6fl00v3Yne6mZmbwoIpQxiaZPz/VFhZH6g/nUi/oZ/IIr2wfHcpbg+MyYxnWEpsl68PSvF/Yw28411eOuomyJnT7bfuK9UypvRC7pFw/oMQFQd7voSH58O2DzCZTMwfmw7A0l2lnV7C4/Hw0kpjGfPyubmYTCZ+d8FUTCb4cFMRBWUtZ81qG508s2wvAD84YSwmk4nsZF8w04yZRB4FM5FeWNKNNhnNZQWj+H/RfVCZDykj4ORf9uitvmA2Ml3BTHpo1tVw42IYPsdoOvvqdVCR7w9my3aVdPr2r/eUsbe0jvhoC2dPzwZg3JBEjhuXCcDzK/JbvP6FlQVUNTgZkxHPaZOHAJDt/WVIwUwikYKZSA81OFx8uOkQ0L36Mmgq/i+ubsQViK7lG1+D5Q8b98/5G0TH9+jtvlmJEQpm0hsZeXDtRzBiHthr4O2fMH+MEczW5FdQ5z3uqz1vrj0IwLkzhrXooXfVkbkAvLRqv3+HpsPl5glvbdn1x43xN3FumjHTUqZEHgUzkR76z/J9FFU1Miw5hhPGZ3brPenx0ZhM4PZAaW0flzN3fw6v3QB44MgbIK/rnmWt7RvEOzLV+SpALFY4759gscGuT8kteJ3hKbE43R5/r7H2+GozT5yQ1eLxUyYNITPRRklNI594m8e+ve4gBysbyEiI5qLZw/2vzU42ZsxU/C+RSMFMpAdqG508/PkuAH50yjhiorp3nJHVYg5M9/+Da+GFq40jlyafD2f8oVeX8S9lpvVspk2khYxx/vM0TR/+ggUjjNi7tIPlTIfLzY6iGgAmZ7ds6xJlMXPZXKPtxn9X5PPF9mLufnMTAN+dP6rF/2vNa8zUZFYijYKZSA88tWQPpbV2RqXHcbG3d1N3+VpmFPd2A0DpLnjuEmPpaNRxcNFjYO75OZd1dqe/ueygWcrU393BM++HkD0TGiq5qewPWHGyrIMNALuKa7C73CTarOSktt00c8URIzCZYPGOEq5duJKaRidHj0nj2mNHt3idrzSg0emmvM4R8D+SSDgpmIl0U2Wdg0e/NOpdbj1tPFGWnv3v4wtmRVW9WH6pPgTPXmj0kBo6Ha74r9G+oBd8jWWTY6NIjo3crv8SIhYrXPgoRCeQWbKCX1qfY+OBSirbCUxbCo1zMSdmJ2I2t23KnJsW598E4HJ7uGjWcJ659ijiWh0bZrNa/DPQqjOTSKNgJtJNf/14G9UNTiYOTeTc6cN6/P6mnZk9nDFrqIT/XAIV+yB1NFzzare6+3dEOzINOq4pgLImwoWPAPA964dcZP6CT7e2PWR8S2E10HYZs7kfnpzH0KQYbjttPH+5bEaHffb8y5kVqjOTyKJgJtINi7Yd5ull+wD4xVmT2v1tvyv+Y5l60jLD5YSXvwtFGyA+C771GiRkdfm2zhQM4sJ/CaJJ58IJ/wvAfdYn2fPl821esvmgMWM2qZNgdsSoNJb/4hR+dMq4TsPzUF8w680MtEg/pmAm0oXi6kZuf3kdYBQhH9/NnZitZSX1ovv/h7+AXZ/hiYqDq1+GtDG9+uzmBlNzWU2KhdgJP6c+72xsJgc/rfgdpe/fB97ifI/H41/KnDys9zO+PsP8M2ZaypTIomAm0gmPx8Mdr6yjpMbOhCGJ/PzMib2+VtN5md0MZquehBWPAvDDhh9w53Kzv2i/L3ytMgb7UqYEgdlM7JXP8HHihQCkf/1HeOsW8Hg4XN1Iaa0dswnGD0ns80cNVcsMiVAKZiKdeHnVfhZtKybaauYfV87qdnuM9vhqYnYerqGizt7ua77cXsz9723h3489iPOdnwHwJ8dlvOOYy/MrCjjpT5/zxFd7+tQioGkpU60yJAgsVjxn/IFfOK7DiRnW/Ac2vMxm72zZ2MyEPv1/5DMsRccySWQKajC77777mD9/PnFxcaSkpLR5fuHChZhMpna/Dh8+HMyhiXSptKaR37+/BYCfnjaeCUP79lv+5OwkJgxJpKbRyV8/3t7iOY/Hw98+3s63n1zBnq9e5Lv778KKizfdx1I66xYeuWY2U4cnUd3o5LfvbObxxXt6NQaX28P+cnX9l+A6eWIWn8Sdxd8dFxsPvH8He/cZ/812Vl/WEzrIXCJVUIOZ3W7n0ksv5aabbmr3+csvv5zCwsIWXwsWLOCEE04gK6tvBc4iffX797ZSUedg4tDENn2UesNqMfPrcycDxukBWw8ZMwiNThc/fWkdD3y6gwXmFTwc/Q+iTS7255zFvNte5A+XzOCMqdm8dfOx3L5gAgB/+GArX+/u/LDo9hysqMfh8hBtMfv/YhMJNKvFzKVzc3jEdS57o8ZCfTkz1/8OCFww83X/V5NZiTRBDWb33nsvt956K9OmTWv3+djYWIYOHer/slgsfPbZZ1x33XXBHJZIl5buKuHV1fsxmeD3F03rcc+yjszPy+CMKUNxe+DetzbzxpoDnPXAYl5bc4BJlgIeinkQCy6Ydhk533uWrJQE/3vNZhP/c+JYLpg5DJfbwy3Pr+nxoei+ZcyctFj/uYMiwXD53BG4TVb+p+Y63CYrs2q+5FuWj5icndD1m7thSLKazEpk6lc1Zs888wxxcXFccsklHb6msbGRqqqqFl8igeTxePxHwVx15Ahmj0gN6PV/efYkbFYzy3aX8pMX17KruJasWA8vZz6Jxe2AcQuMnlAWa5v3mkwmfn/RNMYPSaC4upHvP/MNB7u5K62izs7b6wuBwbEjU8JrRHocd545ic2eUfzLcQ4Av41ayPxPL4Fdi/p8faPJbDSg5UyJLG1/8ofRk08+yVVXXUVsbNujOnzuv/9+7r333hCOSgabjQeq2Hm4hrhoC3ec0ftdmB3JTYvjf07M42+fbCc1LorrjxvD9+ufIHrFNojPhPMf7PSopbhoKw9fM4cLHlzCuoIKzvj7l/z2gqnkpsWxel85u4pribKYiImy4HZ7qGl0UlJjZ/GOYhqdbgCmBKBdwUClVa/Quf640RyoqOeBpRdj90Txg6h3SDi8Hp69ABb8Hubd3KfrD02OoaTGzqHKBqYMSw7MoEXCrMfB7J577ukyGK1cuZK5c+f26LrLli1j8+bNPPPMM52+7s477+S2227zf19VVUVubm6PPkukMx9vMTqWHz8uM2hHFv3olDyOG5/B+CGJJOxfDM8+bDxx3r8goes+aWMzE3jrlmP5yYtrWVdQwY9fWNutz52UncTVR43gsrn6f0aCz2Qycdc5kzlU2cC/Nl3I/rGX8/es94xWMB/9CjLGw7jTen397ORYNh6o4qB2ZkoE6XEwu+WWW7jiiis6fc2oUaN6PJDHH3+cmTNnMmfOnE5fZ7PZsNl6d0agSHd86g1mp04eErTPMJlMzM5JgqX/hM9+B3hgzndhwhndvsbojHhe+cE8/vnZTh75YhcJNiuzR6T4j7tp8M6OJdqsJMZYmZGbwszclEF5FJFHp5iHjcVs4u9XzOSNNQc4Ji8DUk8FtxNWPwOvXAvXfwKZE3p1bV8LmkNaypQI0uNglpGRQUZGRkAHUVNTw0svvcT9998f0OuK9NTBino2HazCZIKTJvSuw3+31JbCS9+GfV8Z3088x1ja6aEoi5nbThvPj07Ow2I2DcrQJf1fTJSFK44c0fTAWX+Bkp2QvxSevwJu/BJsPW9Hkx5v/JJeVqvif4kcQS3+z8/PZ+3ateTn5+NyuVi7di1r166lpqamxetefPFFnE4nV199dTCHI9Il32zZnBGppCcEcWb2g/81QllUPJz3T7j8PxDd+4avVotZoayH9I8rjKzRcPmzkJQDZbvh8z/06jJp8UapQXlt+w2bRQaioAazu+++m1mzZvHrX/+ampoaZs2axaxZs1i1alWL1z3xxBNcdNFFpKYGdvebSE99ssVobBzMZUz2fwMbXgZM8J23Yfa3lRJk8InPgHP+Ztxf/jAUberxJdL8M2ahD2bqnSbBEtRdmQsXLmThwoVdvm7p0qXBHIZIt9Q0Olm2y2jaeuqkIDU49niMg8kBZlwBOZ3XVErfKO72c+NPh0nnwpa34Z3b4Hvvg7n78wWp3hmzsg6OOAuGynoH97+3hTfXHmRYSgwzclI4emw6F8/OUW9ACYh+1cdMJJy+2lGM3eVmVHocYzMD0wSzjS1vQcFysMbCyXcF5zNEBpIz/mAs6Rcsh1VP9KifSVq80ccsVEuZH206xGl//YIXVhZQ73Cxq7iW19Yc4I5X1nPlY8u73VNQpDP9qo+ZSDj5lzEnDQlsvdaKx2DHR8b9g2uM2/k/hOThgfsMkYEqOQdO/Dl8fBe89zNY/TQccT2MOg5SR3Xa0y8tzhvM6uy43R7MQZyx2lJYxY3/+QaPB8ZkxPPr86bgdntYnV/Ok1/tYcWeMs58YDF/vWwGp0wKYimERDwFMxGv1fnlABw7LoC7jre+Z/xl01zCUDjmx4H7DJGB7uiboHK/EcoObYC3vf9/WGwwdCqcfh+MnNfmbaneGTO3x1hi9H0fDGvyK/B4YPaIFP77/aOJiTIC40kTs7h4dg4/fmEN6/ZXcvN/V7P056f4Z/NEekpLmSJAnd3JnpJagMB1EK/cD2/+j3F/+uVGR//zH4Rr3wdbkJZKRQYiSxSc9Ue4bQuc/jvIngHWGHA1woFvYOFZsOh+cDlbvC3KYiYxxphfCHadWUG5cc7stOHJ/lDmMyojnldums/k7CQaHG5eWJkf1LFIZFMwk7DasL+S/d4feOG07VA1Hg9kJNjITAxAmwyXE169HurLIXum0dF/1jXGV9qYvl9fJBLFpRnL/Dd+Cb84CD9aAzOuBI8bvvgDPHM+OFp2+Q9Vndn+cqN+LLeDc2ajLGauPXY0AM8u24fD5Q7qeCRyKZhJ2GzYX8n5D37FpY8so9HpCutYthRWAzA5EGdIejzw4Z2QvwyiE+HSp4y+TSLSfWaL8UvMhY/AxU+ALcno/ff+HS1eluqtMysNejAzfoHMSe34LOdzZ2STkRBNYWUDH246FNTxSORSMJOwefiLnbg9UFjZwJtrDoZ1LJsLKwGYlN3z7uMtuF3w9o9gxb+N7897QDNkIn017RK47GnAZNShrf2v/6n0EM+Y5aS2P2MGYLNauMp7wsFTS/YGdTwSuRTMJCz2lNTy/sam3ygf/2p3WBs2+mfMsnsxY1a8Ddb8xziY+aVvG2cAmsxw/kMw9eIAj1T6To1BB6SxJ8OJdxr337nVaNRM0waAYNaYNThcFFc3Ap3PmAFcc/RIoiwmvtlXzvr9FUEbk0QuBTMJi8cW78bjgaPHpJFgs7K9qIYvtheHZSxut4cthVVAL4LZvqXw0Dx482bjL4ut74A5Ci55CmbpiLH+Qk3aI8Txt8PYU8DZAI+fDP8+kQXVr2LDHtQZM99sWYLNSnJsVKevzUqK4exp2QA8s2xf0MYkkUvBTELucHUDr3yzH4DbTpvAFUfkAkZYC4f8sjrq7C5sVjOjM3pwXmVdGbz6ffC4YMhUmHA2TL0EvvU6TLkgaOMVGbTMZrj4cZhwFpgscHANp+U/wF+iHqa0pjFoH9u8vqw7PQ6vOmokAB9uPBT2+lkZeNTHTEJu4ZK92J1uZo1I4YhRqQxPjeWppXtZsrOUTQcrA9euops2e2fLJgxNxGrp5u8qHo9RS1a136ghu/YDsPWxPk1CSofnDFBxaXDl81BTDJtew/3BLzjH8jXFRW8Cs4Lykd2pL2tu7shUhiTZKKpq5KsdJWo4Kz2iGTMJuTfWHADgxuPHYDKZGJ4Sy1neqf+nl+4N+Xh8y5iThvZgGXPVE8b5fuYouORJhTKRUEvIhKNuZOe0nwJwVdmDcHhLUD6qKZh1Xl/mYzabOHOq8TPt3Q2FQRmTRC4FMwkpl9vDoSqjD9HsEan+x307mT7YeAi7M7T9fzYf9NaXdbdVxs5P4f3/Ne6feg8MC85v6RIAgTxaS/ql6tk38IVrOjbs8PL3wBH48yq70yqjNd8vmx9vLtJypvSIgpmEVFmtHbfH+Puy+ZElR45OIyMhmqoGJ8t2l4Z0TP4Zs+4U/heuM3Zeup1GPdnR/xPk0YlIZ1LjY/ip4yZKPMlQvAU++13AP6Ogh0uZYCxnZiXaqG5wsmRnScDHJJFLwUxCqsRboJsWF92instiNrFgylAA3g/h1H9FnZ2DlcYM3sTOepi5XbBrETx3KdhrYPTxcMFDRjGyiIRNeryNEpK53XGD8cCyByF/eUA/40AvZsyM5UzjZ9q769VsVrpPf6tISPl6AbV37JFv6v/DTYdwhug4E1/hf25aLEkx7WyDb6yBj++Gv06GZy+AmiJjB+bl/wFrAI5uEpE+SYyxYjGbWOSeRf2UKwAPvHET2ANz1Fu93UVJjdGKo6PjmDri+5n20WbtzpTuUzCTkOosmB01Oo20+GjK6xws313WretV1Nl5eVUBd762gbMeWMz1T6/014x1h++17Rb+l+6CJ06DJQ9AzSGISYE53zPaYcSEdueoiLTPbDaRGmf8UrVv7q8gaTiU7TZ+oQpAA7sDFUbAS4zpuodZa3NHpWk5U3pMwUxCqti7lJmZ0DaYWS1mFkwxtpW/t7Hr5cyVe8s47W9fcvsr63l+RT6bC6v4ZMthzvnnYu58bQOVdY4ur7HxgHEU07ThrYLW9o/g3yfB4c2QMAQuewZ+tgPO/TskZHV5XREJHf95ma5YOO8fxoMrH4OPfgXuvs2+F5T1vL7Mx2I2+WfN3gjzsXMycCiYSUj5Zswy2pkxA/xbzD/ceAiXu/3fdj0eD08t2cOV/15OcXUjo9LjuPH4MTx41WzOmZ6N2wPPr8jnV29u7HI8673BbGpOs2C29V14/gporITco+DGL2Hy+TqIXKSf8m0kKqu1Q96pcPp9xhPL/mUsa7q6/iWtI73ZkdncBbOGA8ZyZk2js9fjkMFDwUxCqqSTGTOAeWPTSYmLorTW3uHU/+OL93Dv25txuj2cO2MY7/34OO48axJnT8/mX1fN5rFvzwVg0dbDODqpVatucLCnpBZoNmO26zN4+btGN/9pl8F33oHEob3804pIKPiCWbnvvMz5t8AFjxinA6x/AV76Tq/DWU97mLU2IyeZMRnxNDjcfLBRmwCkawpmElKd1ZgBRFnMnDdjGAB//2R7m4PNi6sb+fsn2wH42enj+ccVM4mLbnmAxSkTs0iNi6Km0dnpIcKbDlbh8cCw5BgyEmzGTq7nrwKXHSadBxc8rFmyCKRzMyOP7yDz0ppm52XOvNI4IcBig23vwqvXg6vnM1a+YJbbi6VMAJPJ5J818zXXFumMgpmEVFfBDODmk/KIiTKzOr+CjzYXtXjub59sp9buYnpOMv9zYl6759aZzSbmjU0HYMnOjnui+erLpg5PNs69fOk74KyHvNPg4ifAohPLIoXCWGRLi2s1Y+YzfgFc8RxYomHzG8ayprtnuyP7upQJcKE3mC3ZVcIhb3sekY4omElI+Yv/OwlmQ5JiuO7Y0QD88YOt/tYZ24uqeWFFPgC/OnsyZnPHXd2PycsA4KtOdkJt8Aaz6cOT4J2fGDsvM8Ybhf6aKRMZMFrUmLU27jS49GkwW2HDS0aD6B600uhNc9nWctPiOGJUKh4PvLlWs2bSOQUzCRm7002Fd6dkRgc1Zj43njCW1LgodhXX8tzX+ewpqeW372zG7YEFU4Zw5Oi0Tt9/zFgjmK3JL6fO3v7yxYb9RjA7xfEFbH7T+MF90b8huvc/gGXgaG+2VQamNjVmrU08Cy55yljW3PoOPH2OcQi6owFqSzucUq2sd/jDXk5a72fMoGkTwOtazpQuKJhJyJTWGrNlVrOJlC76ASXFRHHLyeMA+PVbmzjpz5+zeEcJVrOJn585qcvPGpkex/CUWBwuDyv2tO2JZhT+V3O8eR0TVt9rPHjCz3XuZYRR9Boc2q0xa23yefDtNyE2FQ58A3/Og/uGwJ/GwL9PgIr8Nm/ZdND45S0ntYMG1D1wzrRhRFvMbD1UzYebtAlAOqZgJiHjb5WRYOt0GdLnmqNHMHGocUxSgs3KmIx47jpnMqMz4rt8r8lk4pg8o85s6a5WdWYNVZR99Gc+j76NZ6L/D7O9GnKOgGNv7eGfSET6gw5rzFobOQ+u/xTSx7V8vHAdPHYyFKxo8bCvAfXUYX1vKJ0cF8W13hKNX76+gVJvWYdIa6pulpDpTuF/czarhfd+dBx1DhcJtp7/p3pMXgYvrdrPVzu8dWb1FfD1I7D8YUY2VIAZ6szxxM29Bk64Q8X+IgNUWoI3mNU68Hg8nS9Tp4+Fm1cYNaXRCdBQCS9cCYc2wMJz4MJHYOpFQPMNQu2cDNILt542jkVbD7OtqJpfvL6BR66ZoyV1aUMzZhIyPQ1mYOyw7E0oA5jvrTPbXFhFWWUVPLkAPr8fGioois7lDsf3eXb+h3DWHyE+o1efISLh55sxs7vc3WviajZD0jCISYKUXPjeBzDxHHA1wivXwjcLAdjonTGbEoAZMzB+2fzLZTOwmk18uKmI11ZHRr1Zg8PFu+sL2998IT2mYCYh07SUGZodj5mJNiYMMZZCi9//PyjeCvGZcMmTXBX1AC+5TmLiCDWPFRnoYqMtxEQZf52V1/aikawtwdiNPfdawANv/xj7F39hV3E1AFMCNGMGRnueH59iLKX+4vUNA/4MzeLqRq58bDk3/3c1d3XjtBXpmoKZhExJN1plBNqJEzMZZSpkzLZHjQfO/CNVeeexq9ToJdTmjEwRGZDS442fK8W9rd0yW+Dsv/prTaMX/Yb3o37O9+KXkhUb2L8qbzpxLKdOGkKj0811T69k+e6O+y32Z5sOVnL+v75iTX4FAJ9uKaLe3rM+cdKWgpmETGcHmAfLxbOG8xvrQqI8DuyjToQpF/L+BuOA9FHpcf5t9iIysI0bkgDA59sO9/4iJhOceg8suB+HJZaJ5gJ+7foXPHQ0lO0JzEABq8XMg1fP4qQJmTQ43Fy7cCWr88sDdv1QOFzVwBWPLudgZQNjMuIZkmSjweHmyx3F4R7agKdgJiHTVGMWE7LPHH/wDY63bKDRE8Wbw27D7YHHFhs/YK85emTIxiEiwXXx7BwAXv1mPy53H496mPc//G78K/zBcQW1UWlQtguePAOKNgVgpAab1cLD18zhuHEZ1Nld/PC/a6hu6P1h66H2zLJ9VDc6mZSdxOv/cwxnTcsGUCuQAFAwk5DpTfF/r3k8sPgv8NYPAXjQeT5PbTHz+fbD7DxcQ6LNyuVH5AZ/HCISEqdNHkJybBQHKxsCUre1qsjDI67zWHH6m5A1xdjF+dSZsG9pAEZriIkywlluWiwHKur5zdubA3btYGpwuHju630A/OjkPJLjolgwxajX/XTLYRze01qkdwZVMPtiezHPfb3PHxAktIJS/F99CLa+C1vehs1vwYZX4JunjQOLP/0NAI1zvs+/uYjNhVXc6/3Bd+VRI0jsY8NIGZh0bGZkiomycMHMYQC8tKqgT9dqdLrYXmQU/o/Ly4PvvQs5RxqtNRaeDZ/cC87A/D2SYLPyl0tnYjLBy9/s56MBMOP02uoDlNc5yEmN5XRvIDtiVBpp8dFU1jvabeot3TdoGjdVNzj4/tOrsLvc3PXGRo7Jy+C208Yza0RquIfWr2w+WMXPX1uP1WwiKzGGockxjEqPY0xmAjNyU0juomN/R+rsTmq9RaEBmTGr3A9f/Q1WPwOuDrZom8xw5h+xHfl9TqlazbsbCtlXWofVbOK780f1fQwi0q9cOjeXp5ft46PNRVTU2UmJ690vgTuKanC4PKTERTE8JRZMcfDtN+DtnxjnbX71V9j2Pkw6F9LzYPgcyMjr9biPHJ3GDceP4dEvdnPnaxuYOSKFrBCWfPSEx+PhySVGOch354/C4m0WbjGbOHVSFi+t2s+Hmw75zyuWnhs0wWxXcS12lxuzCdweWLyjhKp6B2/ecmy4h9Zv1Ntd3PL8anYX17b7fHJsFC/eeDQTh/Z863hJtRGeYqLMve5LBoDLYSxRLv5LUyDLnGT0I8IEliijaaQtEWZ/C0YfD8Alc3J411v0f+6MYQxL6du5dzKweDRPNihMGZbEpOwkthRW8da6g1xz1EiqGhxsL6phS2EVxdWNzBqRwtFj0onv5OeQr7HslGFJTQ1go+Ph4sdg0jnwzq1QvMX4AsAEp9wFx95mbCDohdtOG88X24rZeqiaG5/9hue/fzQxUZZeXSuYvtxRws7DNSS0Uw6yYMpQXlq1n482FXHPuVO6dcKLtDVogtnu4hrA+M3k52dO4oIHl7C5sIpGpwubtf/9xx8Ov39vC7uLaxmSZOOucyZTWmPnYEU9u0tq2XigksLKBr731EreuPkYhiT17Le54hqjPUVmoq33na6Lt8PrN8DBNcb3I4+FE/8XRh3X5Q/D48ZlkJMaS2FlA9cfN7p3ny8RRX9lRB6TycRlc3O49+3N3Pv2Zn791qZ2zyePspg4cnQaC6YMZcGUoW1+nq33dfxvr7Hs5PNhxHxj5qx4GxzeDPtXGqUTh7fAef+EqJ7/4ufbDHDBg0tYk1/BL17bwF8um9HvTgZ4yjtbduncnDblIMfkZRAfbeFQVQPrD1QyMzclDCMc+AZNMNvlDWZjMxOYkZNMalwU5XUOthRW6z8eYNG2wzy73Cjm/POlMzhuXGaL5yvq7Fz08FJ2F9fyvadW8tIP5vVo5stf+N+bVhnORljyD1j8Z3A2QEwKnP0XmHZJty9htZh56cZ5VNQ5mDwscM0ipX/rZ3+nSQhcMHM4D3y6g4q6ph2O2ckxTMpOIj0+muV7Sikoq2fJzlKW7Czl129t4rhxmfzPiWOZPSKVv368nedXGAead/h3Q0ImzLu56fuVj8P7/wsbXjbaalz9MsSl9XjsozPiefCq2XznqRW8tuYA44YkctOJY3t8nWAprm7ky+1GO4zvzBvV5vmYKAvHjsvgw01FLN9dqr9be2nwBLPDxvLcmMwETCYT03NS+GJ7Mev3Vwz6/3g2Hqjk9pfXA0bNQOtQBpASF83T3zuSCx8yZhrvfmMjf718Zrc/o/kB5t3mdsOuT+GDO6F0h/HY2FPg/H8Zx6n00LCUWC1hikS41PhovvjZSZTUNpIYYyUpJqrFkqDH42FvaR0fbz7E+xsPsSa/gi+3F/Pl9mJS4qL8ge7KI0f4dxp26YjrIWMCvPQtOLDKaK3xrdcheXiPx3/suAzuOXcyd725iT9+uJW8rAROmzykx9cJhg82HcLtgRm5KYzKiG/3NdNzUvhwUxFbCqtCPLrIMWh2Ze4u8c2YGf8xzcgxpqjX768M25j6g5dXFXDxw0spqWlk4tBEfn7mxA5fm5sWxz+umAXAx5uLcPegV1BxjVEP1q3C/6qDsPRf8K+58NwlRiiLz4KLn4BrXu1VKBORwSM5LoqxmQlkJca0qdMymUyMzojnhuPH8vr/HMPiO07iW0ePJNpqpqLOQUpcFI9cM4f7L5rWsxqp0cfBtR9C0nAo2WaczVu4rlfj/9a8UXzr6JF4PPCTF9aw9VD/CDnvrj8IwDnenmXtmZxtrEgomPXeoJgxc7rc7C2pA4ylTIBpOSkArN9fEaZRhd8f3t/KI1/sAuDkiVn87bKZXRabHjk6jZgoM9WNTnaX1JKXldCtz+pyxqzqoLHLctdnULqz6fHoRJj9bTjhDohN6dZniYh0V25aHL+9YCo/PCWPz7Yc5uSJWWT1sIbWL3OCEc6evdD4hfLR42HyBXDinZDV8S+97bn73MnsLqlhyc5Srn96FW/efAzp7fz89Hg8bC6sorzWQWKMlbT4aHLT4no3/k4crm7ga28bjDOndTyT6CsV2VVcS4PD1S83MPR3gyKY7S+vx+5yY7Oaja3PNM2Y7TxcQ22js9MdOpFo2a5Sfyi79dTx/PDkvG79dmi1mJkyLJlv9pWzfn9Ft4PZocp6AIYmt/MDb9Mb8PaPoaHC+N5khmGzYdY1MO1S44BhEZEgykqM4YojR/T9Qim5cO0H8N7tsOk12PwGbHkL5v8QTvolWLtXzhFlMfPgVbO54MEl7C2t4+x/fMUPT8nj0jm5WMwmCsrq+HTrYV5aWcA2b881nyuPzOX3F04L6MaBDzYewuOBWSNSyEntOPhlJdpIi4+mrNbO9qJqpnsnQaT7BkUa8RX+j86I94ePrKQYhibFcKiqgY0HKjlqTHo4hxhSdqebX72xAYCrjxrBj08d16P3T8/xBbNKLvIeg9KVwkpjV2Z282BWtAmWPADrXzS+z54Jx98Oo47V7JiIDFzxGXDpU3D8z2DR72HrO8bPuh2fwIUPQ/aMbl0mJS6ax79zBN95cgUHKur55esb+etH26m1O2lwNHXXt1nNjEyPo6bBSWFVA8+vKGBSdhLfbqdAv7feWW+0Gzq7k2VMMJaKJ2cn8dXOErYUVimY9cKgCGa+vlxjW83uTM9J5tDmBjYMsmD22OLd7CquJSMhmjsW9Gx6HWBGL5aBm4JZrNGYcfFfjC3mAJjg2FuN6X6rDhUXkQgxZApc8ZxxOslbP4LDm4zlzTEnwpzvwoSzu/yZl5eVwKc/PYHnV+Tz4KJdlNQYZSHRVjNThiVx0ewczpsxzN/8+7Evd3Pfe1v47TubmTIsmTkj+95EvaiqgZV7jWXMs7oIZmAsZ361s4TNB1Vn1huDIpj5W2W02kUyIzeFjzYXsW4QbQDIL63jH58aOxx/dfZkkuN63sl/uncZeNPBKhwuN1GWzveQ1NmdVNY7AA8jNz0Ii+83njBbYeLZMO8WyD2yx+MQERkQJp5tHOn0/u1G6cbuz42v+EyYeTXM+Q6kjenw7TFRFr53zGguPyKXdQWVDE2OYURanL/rfnPXHzeaNQXlvLfhEP/z3De8dcuxPe472dq76wvxeGDOyNRu7Wz3bQDYrA0AvTIodmX6g1k7M2YwuDYA/P3T7TQ63cwfm875M3u3u3FUejyJMVYanW7/eXKdKaxsIBoHD9geJcYXyo68AW7bApc9o1AmIpEvIRMuXQg/XgfH/QwShkBtMSz5O/xjFrz8PaMHWifioq3MG5vO6Iz4dkMZGEuJf7xkBmMz4ymqauTyR5dxoKK+18NucLh4bPFugG7/nTHJvzOzuke798UwKIKZfykzs1UwG54CwL7SOirqOjhvMYJUNTh4z3ss0c8WTOh1YajZbGLa8O63Gyksr+PvUQ9yvulLMFng7L/CWX+ChKxefb6IyICVOtI4vunWTXD5f4zejGBsFPjXEcamgX3LjOPneinBZuWp7x5JTmose0vruOyRZewtaf+ova48uWQPhZUNDE+J5bK5uV2/ARiTGU+01UxNo5P95b0PhYNVxAez8lo7pbVG6BrdaikzOS6KkenG7pINByJ/OfOddYU0ONyMy0pgVh+b6k7vQZ1Z2oo/cZZlBQ6i4KqX4Ijr+vTZIn3R3hE9IiFniTIOQf/Wa3DjYhhzErgdsOLf8NQZ8H+j4Y3/gdrSXl1+RHocL904jzEZ8RyoqOeyR5exoxsrHM2V1dp5eJGxe/+np4/vduuLKIuZCUMSAdhcGPl/twZaxAczX2PZ7OSYdlti+ALGuoKKEI4qPF5aVQDAZXNz+7yNutsNete/zOSd/wbg1eG3w7hT+/S5Ir2hMCb9WvZ0+PYbcM1rMPViiE0DezWsfQ4eOgq2vteryw5LieXFG+cxYUgih6sbufzfy9l0sPtB6Z+f7aC60cnk7CQumNmzUwwmZXuDmTYA9FjEB7NdHSxj+vhmjtbkV4RoROGxvaiatQUVWM0mLpzd82NCWpvu/ee27VA1DQ5X2xdUFcIn98CbxnlyDznPo3DUhX3+XJFA0Tma0u/knQKXPAm374LvvgeZE406tBeuhLd/As6el9xkJtp44YajmTY8mbJaO1f+ezlr8su7fN/iHcX8x3t+8i/OmtSzUxBovgGgZ7N0MiiCmTFjNiaz/XO9Znu3Eq/OL8cTwb9Wv7TSmC07ZVJWz86r7MCw5BjS46Nxuj1NO2/cbtj9Bbx2I/x9mtHJ39XIyrjj+JPzspY9zERCwITSlwxAZjOMOgZu+AKO+TFggm+egmfOh9qSHl8uNT6a575/FHNGplLV4OTahSvZX17X7msbHC5++85mvvXEChwuD6dMzOLYcRk9/szJw4xVFR3N1HMR3y7Dd3h5RzNmk7OTsFnNlNc52FNSy5gOXhdoS3eW8OSSvfz8zInd7p7fW3anm9fXHADodvFmV4yD4JP5YlsRh9Z/BltWG9vAq/Y3vWjEPJj/Q+56Lx4PtWTrAHERke6LioHTfgMjj4VXr4P8pUYftGGzvM/HQspISB0FeadCUsc9xpJionj2uiO5/NHlbDhQyc3PrealH8zDZm2qG1u6s4Rfv7WJHYeNCY1rjh7BL8+a3KuhT/QuZR6oqKe81k5qvHpUdlfEB7P8MiOYjcpof8Ys2mpm2vBkVu0rZ3V+RUiC2cq9ZVz79EoaHG4q6+28dOO8gB6d0dpHmw9RWmsnM9HGCeMzA3PR+gq+736ZP9ueJ31Vs6nqmGSYcpFxnFLOXAAOvvAhYMyyiYhID40/Ha7/BJ6/Asp2Q9WBtq+JSYZrP+r0TM64aCsPXT2bc//1Fev2V3Lv25u59dTx7Ciq5tnl+3h/4yEA0uOj+eMl0zll0pBeDzkpJorRGfHsKall3f4KTpygXfjdFdHBzOPxUFBmbNXNTe14tmbOyFRvMCvnkjndO2Kot7YUVnHtwpX+4zRW7i3n/Y2HutVNuTc8Hg8PenfVXHXkCKxdNIPtUn0FfP0ILHuI+Y2VYIJqUwKJ0842miiOW2D8ludV2+ikqsEJdHBOpoiIdC1zAtzwubERwOltQdFYDeX7YM+XxqHpz11iBLjEjg8Zz02L44ErZvHdp1bw36/z+e/X+f7nzCb41tEjufW08aTE9X2Ga9aIFPaU1LI6X8GsJyI6mJXV2qn3FqYP7ySYzRrhrTPb13VBZF8crmrgW0+soLrByRGjUpk9IpVHv9zN/e9v4ZRJWS2mlANl0bbDbCmsIi7awveOGdX7CzULZDQau3rsaRP5yaHTWWQ+ig3nn9Vu6PMdxZRos5IY0/NTBkRExCsmGWZe2fbx2lJ44jQo2wXPXQrfew9siR1e5oTxmfzs9An86cNtmEwwMi2OKcOTueWkPH9z2ECYPSKV11Yf6NZmA2kS0cGswNvYbkiSrdPQM3tkCgDbiqqpbnAELUA8s2wfJTWNjB+SwOPfOYIoi4nX1xygoKyep5fu5Ybjxwb08zweD//6bCcA1xw9sve/Ae38FF7/AdQeNr7PnAQn/i/Wiefx+b0fU293sbe0lrystj8ICiuNfweaLRMRCZL4dLjmFXj8NDi0Ht76oXHKQCduPimPC2cNJy0+utv9yXpq1ogUANbmV+B2e3q8s3Owiuhdmb5dJzmpcZ2+Lisxhty0WDweWFcQnGZ4Tpebl78xdkb++JTxJMdGERdt5fYFEwD456c7A376wPLdZazOryDaaub6Y0f3/AJOO3x8N/znIiOUpY8z/me/aSlMuRCzxcKEob4mgu1viS6s8B5ersJ/EZHgSRsDV75gnK6y6XXY/GaXbxmWEhu0UAYwYUgicdEWqhud7PR2SJCuRXgw67q+zGe2dznzmyAtZy7aVkxRVSPp8dGcNrmpoPLi2TlMGJJIdaOT9zYcCuhnPrjImC27fG4uWT09xNZRD/+9DJY8YHw/9zr4wWKYcqGxldvLN+29tYMt0b6lzOw+HqIrIiJdyD0Cjr3VuP/uT3t9akCgWC1mZnibuAe7VCiSRHQwKyjr3owZNAWz1UFaC39hhVFgefGcHKKtTf/YzWYTF3kbvr697mDAPm/JzhK+2lmC1WzihuPH9OzNjnpj98/uRRAVD5c9C+f81dia3cok74zZ1kMdzJh5lzKzUxTMRESC7oQ7jHKT2mL44H/DPRr/cmaw/m6NRBEdzPwzZmndnzFbk1+O2x3YRrOFlfUs2mbUZ11+RNs+YmdPN3ZkLt9TyuGqhj5/ntPl5jdvbwaM2rLctK6DqZ+9Fp6/EnZ/boSya16Fyed1+PKJ3hmzjpoI+mbMhiVrKVNEJOisNrjgQTCZYcPL8Or3oXh72IbTNOlREbYxDDQRHcwKulljBkYzvNgoC1UNTv9pAYHy8qr9uD1w5Oi0dhvd5qTGMWdkKh4PvLuhsM+f9/zKArYVVZMSF8VPTh3X/Tce+MZoXuibKbvmFRg5r9O3+GrMCisb2q2RU/G/9DceIveEDxEAhs+Bk35h3N/wEjx4JPz3cvjiT7D9Q6PNRoj4Zsx2Hq6hst4Rss8dyCI2mHk8Hg74a8y6DmZRFrP/P6DluwO3Lu92e3jRexzSFe3Mlvmc65016+tyZmWdg79+tA2AW0/19qJxu2HHJ7DxVdj1GRzaCC5n05tqiuGz++CJ06F0JyRmw7deg5Hzu/y8pJgocrw1fO0tZ/qK/4dpKVPCSFFMBp3jbzf6nk04C/DA9g9g0e+M2uF/zjX+TgiB9AQbI9ONv4PXFlSE5DMHuohtl1Fc3Uij043Z1P36pvlj01m6q5Slu0r51rxRARnHVztLOFBRT1KMtdMmsmdNz+Y372xmdX4FBWV1PVt+bOZvn2ynvM7BuKwErj5qBLhdxkHi655v+UJbMow+zjjJedv74PYGtSkXwtl/hbi0bn/mxKFJ7C+vZ0thFUePSfc/Xt3goLrR11xWS5nS32jrvkS4YbPgyueNX8Z3L4JDG2DvEuPovOcuhrnXwum/g+j2T8YJlNkjUtlXWsfqfeWBO30mgkXsjJmvh1l2cixR3ex2P2+scVDrst2lAasze2GlUfR/4azhnW5LzkqM8Yea3i5nfrm9mIVL9wJw1zmTseKG1280QpnJAiPmQ9ZksCUZTWK3vgNb3jZC2fA5cMlTxlcPQhnAZO+ZaFtbtcw45GsuG2MlwRaxvwNIPxbEk85EBo6hU2H+D+Gif8MPV8FRPzAeX/UkLDzHWDUJotne1aiVe8uC+jmRImL/tvT1MOus439r03OSiY+2UFHnYMuhKqYMS+7TGEpqGvl4cxEAVxw5osvXnzN9GEt3lfLGmgPcePyYHp2fWVzdyG0vrQPg6qNGcHxeGrz2fWP50myFS56EyecbL3a74OAa4zcoRwNMvQiGTOn5H9DLtwFg66GWGwBU+C8i0s9ExcKZ/wfjz4BXroWDq41TA655FdID2+TcZ97YDEwmWLqrlPc3FHKmd/WoqKqBTQcrmTcmg9jo4PVTG2gidsZsfw/qy3yiLGaOHG3MFi3b1fc6s1e/2Y/D5WFGbkq3jrk4a9pQ4qItbD1UzVs9qDVzuz389OV1lNQ0MmFIInedPQk+/IU3lEXBZc80hTIAs8U4YPz42+GUu/oUyqCpl9m2ompczWYafb8d+eoLRESknxh7Elz3MaSMhPI9RjjbtSgoH5WXlcCN3pNtfv7aBgor61mxp4wz/v4l1y5cxRH3fcIdr6xTDZpXBAcz347Mns3WzPcuZy7tYzDzeJqK/q/spOi/uZS4aG46wfiP948fbKPBe85nZ5wuN3e/tZEvtxdjs5r551WziPnmUeNcS4CLHjUOFw+iEWlxxEZZaHC42VNSCxhh8dVv9gNwzoxhQf18ERHphYw8I5xlz4S6Unj2Qvj0Ny03hwXIbaeNZ3pOMpX1Dr71xAqufnw55XUObFYzNY1OXlq1nwseXMKPX1jj380/WEVsMCso8/Uw69lszbyxRp3X17tLcbjcvf78FXvK2F1SS1y0pUfB5PvHj2FYcgwHKup54qs9nb62st7B9xau5D/LjTq2v5yexvjt/4YPf2m84LTfwtSLe/1n6C6L2cS04cay7+trjDC2dFcpBysbSIqxcnqzkw5ERKQfSRwC33sf5nwX8MDiv8ADM+Dp84yNY+/cBm//BN67A/Z8CZ7e1V9HW808cMUs4qIt7Dxcg8Pl4axpQ1lz92m8eMPRXDhrOCYTvLn2ICf/+Qt+8sIaFi7Zw8YDwTkmsT+L+Bqzns6YTc5OIjk2isp6BxsOVPqb4/XUf72d/s+bMaxHhe8xURbuOGMiP3lxLQ8t2smlc3PISmy7q3Tn4WpufGYVMaWbuD16Nd9K3UTSZ1ubXnDE941izxD5/vFjWLG3jCe/2st354/2nwt63sxhQT2LTURE+ig6Ds59AEafAG//2Ni1WbW/7etWPArpeTD72zDhbGPGrQdGZ8TzfxdP5563NnH1USP4yanjMZtNHDUmnaPGpHPdsaO59+1NrNxbzhtrD/LGWqOk55qjR3DveVOxDJJD0CMymLndHg5UGDNmPQ1mZrOJeWPS+WDTIZbtKu1VMNtVXOPvR3b1USN7/P7zZgzjqaV7WVdQwbULV/LnS2cwcWhTjdpb6w7y+qv/5XH+zWibsbmASoxOz7lHG8X8c68N6Za0UydlMTM3hbUFFfzh/a18sNE49/OSOd1bxhURkTCbehGMPRmKNkJFgRHOXE6jLrlyv1G3XLoTPr7b+EobA+MWwPjTYeQxxqkDXTh3xjDOmZ7d7ua2qcOTeenGeSzbXcrKPeWszi/nyx3F/Gd5PuV1Dv562Qxs1sj/RT8ig1lRdQMOlwer2cTQXhyePT/PCGZLd5Vw80k9+40A4O+f7MDtgVMnDWFaTs93dprNJn53/lSuenw5Gw9Ucc4/vuKao0diMZvYW1TG0Xse5CnrewB4rLGY8k4xmgiOPwPi07u4enCYTCbuWDCBqx7/mldXG79pjctKYEYv/vwiIhImsSkw6tj2n1twn3HM06Y3YN9SKNsNXz9sfEUnwOzvwAm3Q2znExqddRwwmUzMH5vhr/d+Z/1Bbn1xLe+uL6Sq3sGT3z2i2y2wBqqI/NP5dmRmp8Rg7cW/QN9/ECv2lPW4CHFLYZV/tuy208b3+LN9puUk88ltJ3D65CE43R5eWLqNxmX/5q591/J9byhzz/4Optt3wBXPwayrwxbKfObnZXBMXtMYLp2b06OWHyIi0o/ZEo3VmO+8BXfshsuehVnXQMIQsNfA8gfhgZmw7CFwtj2irzfOmT6MJ797BHHRFhbvKPFvKgu3msbAb5DwicgZs4Iyo76sJ60ymsvLSuDI0Wms2FPGY1/u4e5zJ3f7vX/72Dgs9uzp2Uwe1nWLjDY8HijfC/tXMaRwLY/GF1E64gCxpRuJdxl9whwxGURd+C/ME87s+fWD7GenT2DJzqVYzSYumDU83MMREZFgiEmCyecZX2437PrUWN48vBk+vBNWPgan/QYmntPnsprjxmVy22nj+d27W3jw851cPCcnLLNmpTWNfLS5iI82HeLLTQVB+5yIDGZ7S3tX+N/cLSfl8e09K/jvin3cfNJY0hO6Xjv/Zl85H20uwmwyzqnsNo/HOCpjw0uw8TWoOuB/ygRk+L5JGQnzbiZq5tVga3sYen8wa0Qqj1wzh9hoS7ubFkREJMKYzTDuNBhzEqz9j3H2ctluePEaSM41atQ8boiKN5Y5E7LgiOtg9PHd/oirjxrJI1/soqCsnjfXHuSSOTlB/AO1taOomoseXkp1gzFT5u5D14auRGQwW5NfDsC0nJReX+O4cRlMG57MhgOVPLVkLz9bMKHT1+8uruHGZ1cBcOGsHPKyuhGcPB7Y+Qks+r3RfdnHHAXZM4xjklJGQHwGJA03DhU39//CxzOmDg33EERaaL7Dv5e7/UWkKxar0XZj6sWw5AFY+k+o7GBmafMbMOtbcPpvu6xJA4iNtnD9cWP4w/tbeXDRTi6cNTxkuzTtTjc/eXEt1Q1OxmTGc/HsHObnxjL778H5vIgLZk6Xm9X7jGA2d2TvWl2AUYB480lj+cF/VvP0sr3ccMIYkmKi2n1tQVkdVz/+NSU1diZnJ3H3Od1Y+tz/DXzwc9i/wvjeYoMJZ8C0SyHvVOPYDBEJGpU/igSJLRFO/pVxJmfJDmNCwWQ26tDqymDPF/DNQljzLOz4CC56DMac0OVlv3X0SB79Yhd7Smp5Z/1Bzp8ZmnKZBz7dzqaDVaTERfHC948mKymGqqqqrt/YSxFX/L/1UDW1dheJMVbGD0ns07VOnzyUvKwEqhuc/OmDbS2OGwKju/9Hmw5xxb+XU1jZQF5WAs9edyTJce0HOMA4m/LjX8MTpxqhzBoD826BWzcZRydNOlehTCRAlL1Ewig+A0bOg9wjjWMAx5xotOQ49wGjqW36OKgpgmcvgC//ZNSqdXY5m5Xrjh0NwD8/24nbHfzp71V7y3j4810A3H/hNLJ60emhpyIumK3yns84e0Rqn6c5zWYTPzl1HADPLt/H1Y8v51BlA3tKannlm/2c+6+vuOHZbzhQUc/I9Dieu/6ojmvRXA7Y8Ao8ejws+bux3j7tMvjxOmMLckJmn8YqIiIyYIycDz9YbOzq9Ljhs9/BM+fBnsWd1ht8e/4oEm1Wdh6u4YvtxUEdosfj4c7XNuD2wEWzhvsPXw+2iFvKXBWAZczmzpk+DIfLzS9f38jy3WXM+8OnLf6biYu28L1jRnHDcWPbzpS5XXBwrVFHtvrppqL++Cw49+9BP8NSRESk34qKhfMfhBHz4N2fwt7Fxtew2UYLqLxTIXVUi7ckxURxxZG5PLZ4D49/tZuTJmYFbXir9pWz43ANcdEWfn3elKB9TmsRFcw8Hg+r9nqD2ai0gF33wlk5zMhJ4Zb/rmFzYRXRVjNThyVxTF4G350/qu0sWfF2WPYvo7ixodk5X/GZcOQNcMT1EBe48YmIiAxYs64xTg5Y9i9Y8x9jM5xvQ1zqaBg6FTLGGycNJAzl+nFJPLfEwZKdpWw+WNW71lTd8NJKY+PC2dOySY7tpEQpwCIqmB2oqOdQVQNWs4mZuSkBvfaYzATe/uGx7C+vIzs5lmhrq1Vgj8fohLz0n7D9/abHY5KNLcETzoIpF0GUWkiIiIi0kDYazv4LnHgnrH4GdnwMBV9D+R7jq5khwDcx8bxqP5oPPq5h8rcvC/huntpGJ+9uKATgsiNCe7RgRAUz32zZlOHJxEYHvq2ExWxiZHp80wMej3F+2L4l8PWjzVpemIwgdvRNA6bFhYiISNjFZ8BxtxlfDZWwf6Wxs7NkO5Tvg5rDULWf2PpyrrF+Cns+xf7EQqJPvwdGHBWwYby7oZA6u4vRGfEBK43qrsgKZvuMwv+g/UN0NkLBCuM/lAPfGLc1RU3PW2ww8yqYdzNkjAvOGERERAaDmGSjzizv1JaPu92wdzFfvfRXjqj/Ctv+pfDk6ZB3mtFHLe/UPq9OvbzKWMa8ZE7ojxaMrGC2N7CF/4AxK7bnC+Pg1i1vt6wZAzBbYcgUGH+mUTum3ZUiIiLBYzbDmBNoOG8iJz7zIT+Jfp3LLF9g2vkx7PwYbMkw5QKY/yPIyOvx5XcX17BybzlmE1w8O7QnDEAEBbPKegfbiqoBmDMqAMHM7Yatbxu9VQ5taHo8YYixgyRnLgyfa3Toj+7dmZwiIiLSO6dMyuLfo8bxv3uvZ/eU67gza5lxrGH1QaMTwppnYfIFRt1aZvePSXzFe1D6CeMzGZoc+rrwiAlmL68qwOOBMZnxfT+jceen8NFdcHiT8X10Aky7xOjKP2K+kdZFREQkbEwmE784exIXPLiEf282ce7JdzD1tN8add/L/gXbP4BNr8HWd+D42+HYW8HS+e5Kl9vDq6uNYHbp3NAW/ftETMJ4bPFuAG4+sefTln6lu+C5y+A/FxmhzJYEx98BP9lgdCoedaxCmYiISD8xMzeF82YMw+OB+97dgsdkgtHHwVUvwg++gnELwGWHRffBv0+E3V902sD2yx3FFFU1khoXxSmTgtcjrTMRkzJqG11Mz0nmwlm9ODvL44FVT8LDx8COD426saNuMrryn/xL9RwTGfCafhDrEHORyHL7gglEW80s213Kx5ubbcgbOs0IaBc/AbFpULTROF3gqTNh2wdQU9zmB4Kv6P+CWcOxWcPTUSFiljIB7j5nMuaeHsNUXQTv3Arb3jW+H32C0UtFuypFIprO0RSJDLlpcVx37Gge/nwX97y1iWPyMoi3eeONyWSUIo0+Ab74P6P2LH+Z8QUQmwoZEyBjHHXJeRRtMWNiNJfOCc8yJkRQMDtjypCedfuvL4clDxj9xxx1YI6CU38NR9+s5UqRCBHiXe4iEiY/Onkc76w/SEFZPX/5aDt3nzu55QsSMuHsPxv90ZY8ANveh4p8IwsULIeC5cQBr1qhPCqF1FXnwQn/C8naldlrt57WxY6LhirY/bnRSfjAaihcawQygOFz4Jy/GTssRUREZECJjbbwuwum8Z0nV7Bw6R4unDWcaTnJbV+YNAzO/D/jy14HZbugeBuUbGfp0i+YZl9HKhXG6QNb34WLH4exJ4f0zxIxwWy4pwhWf2x02U8dZbS1KNkOheuMHRr7loHb0fJNWZPh5F8ZXfr1q7WIiMiAdcL4TM6bMYy31h3kf19dz4s3Hk1iTCe7MKPjjDq0odNYW1DBVR/OIN7q5uurbCR8+Vs4tB6evQhO+iUc/7OQ5YSICWY8ehzYuviHljYWxpxozJANnw2ZExXIREREIsRd50zmi+3FbC6s4rJHl7Pwe0cwJKnzFloOl5tfvGb0K10wLZeEyTNh3HHw/h3GzNmi30FULMy/JQR/gkgKZiarcU6WNRrK90JVoXESffaM/2/v3oOaPPc8gH+JTUIESYEgISqR9VQ5axSUeoEdi1WBsoDttkPFmZ2hPaed0W1q6fAPbW0DWKbWcmhnqg7TltW2OxZnjtrL8TJeDmpZqKUxTAPaFsolXkCPDhAummD47R+s75xAgIC54u8z8/yR93ny5sl33gk/3utwEfaHDUD4Am/PkjHGGGNuEjFLii//vBJ/2l+Hyx1m/Mee/8Vfno/HypgwzBjj4sA9Vc241GFG6Ewx3vj3Pw4vFMuAjR8D4Y8Bp94ebhGxwGMbHK7DlaZPYfZ6AxCh8vYsGGOMMeZFS+c+iiP/9W/I3fcjWv7Rj82f/oDwIAnWLpqN9X+cjTWPKYRDnI3Xe7D7780AgKKnNYiYJbVfWdKrw6dFGb4E/von4OUzbr9rw/QpzKTB3p4BY4wxxnzAvLCZOLQlCTuOXsKpSzdwu9+KQxev4tDFqxDPCMAi5SzYhoDr3Xdwb4iQrlEia2nU6BUFBAzfQutW0/DVm//zLPCfRwCJ+24+O30KM8YYY4yx/xcaJEHZ8/EYtA3hp7Yu/P2XGzhz+SZabvWj4ZpZGKcIlmLHMxoEjHXO+SNSYNOXQEUq0NUKVKQAWf/ttnm7tTArKSnB0aNHUV9fD4lEgu7u7lFj6urqUFBQAL1ej4CAAKxYsQK7du1CfHy8O6fGGGOMsYeAeIYIiQvCkbggHG9l/Cta/tGH5pt9CBTPQKB4Bv4wOxhhQZLxVxI8G/jzSeDA88B1A3Bgk9vm69Y7qVqtVmRnZ2Pr1q0O+3t7e5GWlobo6GhcuHAB1dXVCAkJQVpaGgYHBx2+hzHGGGNsqv4lIhipi5V4YmEEVsaETVyU3Rc8G8j9G/BYKmC767b5uXWPWVFREQBg//79Dvt//fVXdHV1obi4GPPmDT/+QKfTYenSpTCZTFiwgK+iZIwxxpiPkAYDOV8BxwoBlLjlI7x6jtmiRYugUChQUVGBN998EzabDRUVFVi8eDHUarXD91gsFlgsFuF1T08PAMBsNjsczxh7eA309WLIMgDLQJ/wG9HXa8aQZQA28T3+3WCMTYl52RYAJaARD0F3CfKAffv2kVwud9jX0NBACxYsIJFIRCKRiGJjY6m9vX3Mdel0OgLAjRs3bty4cePm1fb777+7vGYKIJpcuVdYWCgcohxLXV0dHn/8ceH1/v37kZeXN+rk/zt37mDt2rWIjY2FVquFzWZDaWkpfvnlF9TV1UEmk41a98g9Zt3d3VCr1TCZTJDL5ZP5KmySzGYz5s2bhytXriAkJMTb05m2OGfP4Jw9h7P2DM7Zc3p6ehAdHY2uri48+uijLl33pA9larVa5OTkjDtm/vz5Tq3rwIEDaGtrQ21tLUQikbAsNDQU33zzjcPPkUqlkEqlo5bL5XLeED0kJCSEs/YAztkzOGfP4aw9g3P2nPu1iytNujBTKBRQKBQu+fCBgQGIRCK7e4fcfz00NOSSz2CMMcYY8xduvV2GyWRCfX09TCYTbDYb6uvrUV9fj76+PgBASkoKurq68Morr+Dy5ctobGzEiy++iEceeQRPPvmkO6fGGGOMMeZz3HpV5jvvvIPPP/9ceL1s2TIAQFVVlXBu2XfffYeioiIkJiZCJBJh2bJlOHHiBKKiHDwawQGpVAqdTufw8CZzLc7aMzhnz+CcPYez9gzO2XPcmfWkT/5njDHGGGPu4dZDmYwxxhhjzHlcmDHGGGOM+QguzBhjjDHGfAQXZowxxhhjPsJvCrOSkhIkJSVh5syZY95lNyAgYFQrLy+3G2M0GpGcnAyZTIY5c+aguLjYPc+68mPOZG0ymZCVlYWgoCAoFAps27YNVqvVbgxnPXnz588ftQ0XFBTYjXEmezaxvXv3IiYmBoGBgUhISMD333/v7Sn5tcLCwlHbrlKpFPqJCIWFhVCpVJDJZFi7di0aGxu9OGP/cP78eWRlZUGlUiEgIABff/21Xb8zuVosFrz66qtQKBQICgrCxo0bcfXqVQ9+C/8wUdYvvPDCqG189erVdmNckbXfFGZWqxXZ2dnYunXruOP27duHjo4OoeXm5gp9ZrMZKSkpUKlUqKurw8cff4zS0lKUlZW5e/p+ZaKsbTYbMjIy0N/fj+rqalRWVuLQoUPIz88XxnDWU1dcXGy3DW/fvl3ocyZ7NrGDBw8iLy8Pb731FgwGA9asWYP09HSYTCZvT82vLV682G7bNRqNQt+uXbtQVlaG3bt3o66uDkqlEikpKejt7fXijH1ff38/4uLisHv3bof9zuSal5eHI0eOoLKyEtXV1ejr60NmZiZsNpunvoZfmChrAHjqqafstvFjx47Z9bska5c/fdPNxnsgOgA6cuTImO/du3cvyeVyunv3rrDsvffeI5VKRUNDQy6eqf8bK+tjx46RSCSia9euCcu++uorkkql1NPTQ0Sc9VSp1Wr68MMPx+x3Jns2sZUrV9KWLVvslsXGxlJBQYGXZuT/dDodxcXFOewbGhoipVJJO3fuFJbdvXuX5HI5lZeXe2iG/m/k3zhncu3u7iaxWEyVlZXCmGvXrpFIJKITJ054bO7+xlE9kZubS08//fSY73FV1n6zx8xZWq0WCoUCK1asQHl5ud2jnWpra5GcnGx3Q7i0tDRcv34dbW1tXpitf6qtrYVGo4FKpRKWpaWlwWKxQK/XC2M466l5//33ER4ejvj4eJSUlNgdpnQmezY+q9UKvV6P1NRUu+Wpqamoqanx0qymh6amJqhUKsTExCAnJwctLS0AgNbWVnR2dtplLpVKkZyczJk/AGdy1ev1GBwctBujUqmg0Wg4+yk4e/YsZs+ejYULF+Lll1/GzZs3hT5XZe3WO/972o4dO7B+/XrIZDKcOXMG+fn5uHXrlnAoqLOzc9QD1iMjI4W+mJgYT0/ZL3V2dgq53RcaGgqJRILOzk5hDGc9ea+99hqWL1+O0NBQ/Pjjj3jjjTfQ2tqKzz77DIBz2bPx3bp1CzabbVSOkZGRnOEDWLVqFb744gssXLgQN27cwLvvvoukpCQ0NjYKuTrKvL293RvTnRacybWzsxMSiQShoaGjxvD2Pjnp6enIzs6GWq1Ga2sr3n77baxbtw56vR5SqdRlWXt1j5mjk0VHtp9++snp9W3fvh2JiYmIj49Hfn4+iouL8cEHH9iN+ecHpgMQTkYfuXy6cXXWjvIiIrvlD2vWI00m+9dffx3JyclYunQpXnrpJZSXl6OiogK3b98W1udM9mxijrZPznDq0tPT8dxzz2HJkiXYsGEDjh49CgB2j+XjzN1jKrly9pO3adMmZGRkQKPRICsrC8ePH8dvv/0mbOtjmWzWXt1jptVqkZOTM+6YkXtdJmP16tUwm824ceMGIiMjoVQqR1Wt93dDjvyPY7pxZdZKpRIXLlywW9bV1YXBwUEhx4c565EeJPv7V/w0NzcjPDzcqezZ+BQKBWbMmOFw++QMXScoKAhLlixBU1MTnnnmGQDDe2/++TnInPmDuX/V63i5KpVKWK1WdHV12e3JuXnzJpKSkjw74WkmKioKarUaTU1NAFyXtVf3mCkUCsTGxo7bAgMDp7x+g8GAwMBA4ZYPiYmJOH/+vN05OydPnoRKpXqgAtAfuDLrxMRENDQ0oKOjQ1h28uRJSKVSJCQkCGMe1qxHepDsDQYDAAg/us5kz8YnkUiQkJCAU6dO2S0/deoU/6FyIYvFgsuXLyMqKgoxMTFQKpV2mVutVpw7d44zfwDO5JqQkACxWGw3pqOjAw0NDZz9A7p9+zauXLki/D67LGunLxPwsvb2djIYDFRUVETBwcFkMBjIYDBQb28vERF9++239Mknn5DRaKTm5mb69NNPKSQkhLZt2yaso7u7myIjI2nz5s1kNBrp8OHDFBISQqWlpd76Wj5poqzv3btHGo2G1q9fTxcvXqTTp0/T3LlzSavVCuvgrCevpqaGysrKyGAwUEtLCx08eJBUKhVt3LhRGONM9mxilZWVJBaLqaKigi5dukR5eXkUFBREbW1t3p6a38rPz6ezZ89SS0sL/fDDD5SZmUmzZs0SMt25cyfJ5XI6fPgwGY1G2rx5M0VFRZHZbPbyzH1bb2+v8BsMQPiNaG9vJyLnct2yZQvNnTuXTp8+TRcvXqR169ZRXFwc3bt3z1tfyyeNl3Vvby/l5+dTTU0Ntba2UlVVFSUmJtKcOXNcnrXfFGa5ubkEYFSrqqoiIqLjx49TfHw8BQcH08yZM0mj0dBHH31Eg4ODduv5+eefac2aNSSVSkmpVFJhYSHfvmGEibImGi7eMjIySCaTUVhYGGm1WrtbYxBx1pOl1+tp1apVJJfLKTAwkBYtWkQ6nY76+/vtxjmTPZvYnj17SK1Wk0QioeXLl9O5c+e8PSW/tmnTJoqKiiKxWEwqlYqeffZZamxsFPqHhoZIp9ORUqkkqVRKTzzxBBmNRi/O2D9UVVU5/D3Ozc0lIudyvXPnDmm1WgoLCyOZTEaZmZlkMpm88G1823hZDwwMUGpqKkVERJBYLKbo6GjKzc0dlaMrsg4g4luxM8YYY4z5gml3HzPGGGOMMX/FhRljjDHGmI/gwowxxhhjzEdwYcYYY4wx5iO4MGOMMcYY8xFcmDHGGGOM+QguzBhjjDHGfAQXZowxxhhjPoILM8YYY4wxH8GFGWOMMcaYj+DCjDHGGGPMR3BhxhhjjDHmI/4P3ZpBI+4KZ1MAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1, figsize=(7, 7))\n", "\n", "ax.plot(rx[:,0]/cgs.au, np.log10(smoothx[\"rho\"]+1e-99))\n", "ax.plot(rz[:,2]/cgs.au, np.log10(smoothz[\"rho\"]+1e-99))\n", "ax.set_xlim(-150, 150)\n", "ax.set_ylim(-18, -13)\n", "\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "matse", "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.10.12" }, "orig_nbformat": 4 }, "nbformat": 4, "nbformat_minor": 2 }