{ "nbformat": 4, "nbformat_minor": 0, "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.6" }, "colab": { "name": "buckley_leverett_1d_notebook.ipynb", "provenance": [], "collapsed_sections": [], "include_colab_link": true } }, "cells": [ { "cell_type": "markdown", "metadata": { "id": "view-in-github", "colab_type": "text" }, "source": [ "\"Open" ] }, { "cell_type": "markdown", "metadata": { "id": "6t7dTDp6Iwst" }, "source": [ "# Buckley-Leverett 1D Simulation using Finite Difference (Forward-Backward Scheme)" ] }, { "cell_type": "code", "metadata": { "colab": { "base_uri": "https://localhost:8080/" }, "id": "Iwri8UT2Fmau", "outputId": "f5163ae1-8242-4990-9b8d-ee3e4e56ab00" }, "source": [ "!git clone https://github.com/yohanesnuwara/pyreservoir" ], "execution_count": 1, "outputs": [ { "output_type": "stream", "text": [ "Cloning into 'pyreservoir'...\n", "remote: Enumerating objects: 12, done.\u001b[K\n", "remote: Counting objects: 100% (12/12), done.\u001b[K\n", "remote: Compressing objects: 100% (12/12), done.\u001b[K\n", "remote: Total 927 (delta 6), reused 0 (delta 0), pack-reused 915\u001b[K\n", "Receiving objects: 100% (927/927), 24.37 MiB | 33.10 MiB/s, done.\n", "Resolving deltas: 100% (495/495), done.\n" ], "name": "stdout" } ] }, { "cell_type": "code", "metadata": { "id": "x5uREYoJcfFq" }, "source": [ "import numpy\n", "from matplotlib import pyplot\n", "import pandas as pd\n", "from scipy import interpolate\n", "\n", "import sys\n", "sys.path.append(\"/content/pyreservoir/fluid_flow\")\n", "from relperm import interpolate_relperm\n", "from twophase import buckley_leverett1d\n", "\n", "pyplot.style.use('seaborn')" ], "execution_count": 2, "outputs": [] }, { "cell_type": "markdown", "metadata": { "id": "UYO7nN_LcfFs" }, "source": [ "Relative permeability data" ] }, { "cell_type": "code", "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 620 }, "id": "eq5bX_j4cfFs", "outputId": "d50c06b6-eaf3-48d9-a279-b3ca30c9d46c" }, "source": [ "Sw = numpy.arange(0.2, 0.9, 0.05)\n", "krw = numpy.array([0, .002, .02, .04, .07, .11, .15, .22, .3, .4, .5, .6, .7, .8])\n", "kro = numpy.array([.6, .5, .4, .3, .23, .17, .12, .08, .05, .03, .02, .01, .005, 0])\n", "\n", "df = pd.DataFrame({\"Sw\": Sw, \"krw\": krw, \"kro\": kro})\n", "print(df)\n", "\n", "pyplot.plot(Sw, krw, '.-', label='krw')\n", "pyplot.plot(Sw, kro, '.-', label='kro')\n", "pyplot.xlim(0.2, 0.85)\n", "pyplot.xlabel('Sw'); pyplot.ylabel('Relative Permeability')\n", "\n", "pyplot.legend()\n", "pyplot.show()" ], "execution_count": 3, "outputs": [ { "output_type": "stream", "text": [ " Sw krw kro\n", "0 0.20 0.000 0.600\n", "1 0.25 0.002 0.500\n", "2 0.30 0.020 0.400\n", "3 0.35 0.040 0.300\n", "4 0.40 0.070 0.230\n", "5 0.45 0.110 0.170\n", "6 0.50 0.150 0.120\n", "7 0.55 0.220 0.080\n", "8 0.60 0.300 0.050\n", "9 0.65 0.400 0.030\n", "10 0.70 0.500 0.020\n", "11 0.75 0.600 0.010\n", "12 0.80 0.700 0.005\n", "13 0.85 0.800 0.000\n" ], "name": "stdout" }, { "output_type": "display_data", "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAe8AAAFYCAYAAAB6RnQAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdeVyU9f7//8c1DDuI7AqKIIosgoBb5Zpbmpplm5lpx9I8WWrLseJ8y9/pc/KctOWktpgtZpniQqVmalaWpYYLIIKIIgKyg8giOzO/PyiKEMeF2eB1v9265TBwzZOXg0+ua655X4pWq9UihBBCCLOhMnYAIYQQQlwbKW8hhBDCzEh5CyGEEGZGylsIIYQwM1LeQgghhJmR8hZCCCHMjNrYAa5WfX0DJSWVxo5hUpyd7WQmfyEzaU7m0ZLMpCWZSXOmNA93d8fLftxs9rzVagtjRzA5MpOWZCbNyTxakpm0JDNpzhzmYTblLYQQQohGUt5CCCGEmZHyFkIIIcyMlLcQQghhZqS8hRBCCDMj5S2EEEKYGSlvIYQQwsxIed+gnTu3s2rV/4wdQwghRAei1xXWli5dSkJCAoqiEBUVRVhYWNN969evZ9u2bahUKvr27cs///lPfUYRQggh2g29lXdsbCwZGRlER0eTlpZGVFQU0dHRAFRUVPDhhx+yZ88e1Go1s2fPJj4+nvDwcH3FaVJdW0920SW83eyxsWrbb/+991ZhY2PD2bNpFBUVYmFhwbx5TxIS0penn36SAQMGMX36Q3z66ce4ubkzYcKkNn18IYQQ7Ud1bT2nMi7Qp4dLi/v0Vt4HDx5kzJgxAPj7+1NaWkpFRQUODg5YWlpiaWlJZWUldnZ2VFVV4eTkdEOPt+n7MxxOKbji52i1Wkov1dKg0WKhUnCyt0JRlFY/f2CgB/eN6nVVj//993spKMhnwIBBHDz4C++99xHbt39JUlIigYFBWFioSElJBiAxMYGnn37+6r85IYQQHUp5ZS1LPorlYkUt21+f0uJ+vZV3UVERISEhTbddXFwoLCzEwcEBa2tr5s+fz5gxY7C2tmbixIn4+fnp3GZrC7QD2NpZYWHRehED1NVradBoAWjQaNEC6it8ja2d1RUfE8DR0Ybs7Ax+/nkfO3fu5OuvvyYyMhwPj06MGjWMlStXUlKSS79+oZw6dQo3NwdKS0sIDe19xe1eLV35OiKZSXMyj5ZkJi3JTJoz5jxq6xpYtiGOixW1rX6Owa4qptVqm/5cUVHB6tWr2bVrFw4ODsyaNYuUlBQCAwOvuI3CwvJW75t8kw+Tb/K54tdX19bzf58cIbe4kq6udrw4a4DOQ+dXekyA8vJqMjOz8PHxY/PmL2loaKCuTkthYTkODm5kZmbx448H8PcPpKCgmO3bd+Hr669zu1fD3d2xTbbTnshMmpN5tCQzaUlm0pwx51Fb18CqmERSMkqwtlRRU6e57Ofprbw9PDwoKipqul1QUIC7uzsAaWlpdO/eHReXxuP4AwYM4MSJEzrL+0bZWKl5cdaANn/N++abh/Lgg7N4/PFHmDjxjmb3eXh0Yf/+fbz88n8oKytj06YNjB8/sU0eVwghRPtRW9fAyq3HSTpXQpi/K49MDKKOyx8d1ttbxYYMGcLu3bsBSEpKwsPDAwcHBwC8vb1JS0ujuroagBMnTuDr66uvKM3YWKnx93Jq85PVnJ2deeSRx1i/fl2zj0dERFJYWECnTk6EhIRy5EgsERH92/SxhRBCmLeaugbe2tJY3OG93Jh/VyiOdlaXPVkNQNH++Xh2G3vttdc4cuQIiqKwZMkSkpOTcXR0ZOzYsWzcuJGYmBgsLCyIiIhg8eLFOrcnh3Wak0NdLclMmpN5tCQzaUlm0pyh51FT28BbWxJIybxIRG83/n5nX9QWqqYsl6PX8m5r8uRqTn7gWpKZNCfzaElm0pLMpDlDzqO6tp63Nh/nVNZFIgPcmTclpKm4f89yOQY7YU0IIYQQf6iured/mxJIPV9K/z7uPHZH8+K+EilvIYQQwsCqaup5c3MCZ86XMiDQg7mTg6+6uEHKWwghhDCoqpp63tgUT1p2GYOCPJgzORgL1bWdPy7lLYQQQhhIZXU9b26KJy2njJuCPXlkUtA1FzdIeQshhBAGUVldx+vRCaTnlnFziCePTAxGpbryyqCtkUuC3iC5JKgQQghdLlXX8drGeNJzyxjSt8sNFTfInrcQQgihVxVVdby+MZ6M/HKGhnbl4QmBN1Tc0AHLu7q+htxL+XS198RGbd2m2/7rJUH/9a+l/PDDd3z33R4Ahg0bwYwZD7fpYwohhDBdFVV1vLYhjsyCCob368rM8YGornA1y6vVbso75swO4goSr/g5Wq2W0toyNFoNKkWFk1WnK14SNMIjlKm9ru6a25e7JGhubg7ffLOdNWsal0ydO3cWt946Bm/vblf/jQkhhDBL5ZW1vLYxnqyCCkaEe/HQbX3apLihg73m3aBtQKNtvEKLRquhQdvQJttNTz/Lu++u5Lnn/h8AQUHBKIrC6dOnCAkJRa1Wo1arCQ3tx5kzqW3ymEIIIUxXWWUtyzfEkVVQwa0R3m1a3NCO9ryn9pqkcy+5ur6GZUdWkl9ZgKedB4sHPNkmh87z8nLw8+vJvn3fAaBWW/52j9LsUqh1dXUoSof6fUkIITqcsku1LN8YR3bhJUZFevPg2IArHuW9Hu2mvK+GjdqaxQOebPPXvFu7JGhAQB8++uh96uvrAUhOTmLmzNlt8phCCCFMT+mlxj3unKJLjO7fjeljerd5cUMHK29oLHA/J5823+7vlwR9/fVXmwq8a1cv7rjjLp58ci4ajZbJk6fQpUvXNn9sIYQQxnexooblG+LILa5k7IDuTBvdSy/FDXJVMbMmVwJqSWbSnMyjJZlJSzKT5q5nHiXlNSzbEEf+hUpuG9Sd+25tm+KWq4oJIYQQelBSXsOyz4+RX1LF+ME+3DvSX2973L+T8hZCCCGu04WyapZtiKOgpIrbb+rB3SN66r24QcpbCCGEuC7FpdUs23CMwovVTLqlB3cNM0xxg5S3EEIIcc2KSqtY9nkcRaXV3DHElylD/QxW3CDlLYQQQlyTootVLNvQWNxThvoxZaifwTNIeQshhBBXqeBiFcs/P0ZxWQ13DvPjjiGGL26Q8hZCCCGuSkFJJcs2xHGhrIapw3sy6RZfo2WR8hZCCCF0yL/QWNwl5TXcPaInE2/2NWoeKW8hhBDiCvIuVLLs82NcrKjl3lv9mTC4h7EjSXkLIYQQrcktvsSyDXGUVtRy/6he3Dao7ZfXvh5S3kIIIcRl5BQ1FnfZpVqmje7NuIHdjR2piV7Le+nSpSQkJKAoClFRUYSFhQGQn5/Ps88+2/R5WVlZPPPMM0yePFmfcYQQQoirkl1YwfINcZRV1jF9TG/GDDCd4gY9lndsbCwZGRlER0eTlpZGVFQU0dHRAHh6evLpp58CUF9fz0MPPcSoUaP0FUUIIYS4KtW19fx4LIv3Yo5TUVXPg2MDGN2/m7FjtaC38j548CBjxowBwN/fn9LSUioqKnBwcGj2eV988QW33XYb9vb2+ooihBBC6FRdW89LH8ZSVFoNwLTRvU2yuEGP5V1UVERISEjTbRcXFwoLC1uU9+bNm/noo4+uaputXRqtI5OZtCQzaU7m0ZLMpCWZCXx/JKupuAEGhHQx2bkY7IS1y102PC4ujp49e7Yo9NbI9Wabk2vwtiQzaU7m0ZLMpCWZCWTklbM6JqHpdldXO+zUitHnYvDreXt4eFBUVNR0u6CgAHd392afs2/fPm6++WZ9RRBCCCF0OpdXxusb46mqaWDGuADCAz2xUyvYWJnuG7JU+trwkCFD2L17NwBJSUl4eHi02MNOTEwkMDBQXxGEEEKIK0rPLeO1DfFUVtcze2IQoyK70aeHi0kXN+hxzzsyMpKQkBCmTZuGoigsWbKEmJgYHB0dGTt2LACFhYW4urrqK4IQQgjRqrScUt6ITqC6tp5HJwVzc98uxo501fT6q8Wf38sNtNjL3r59uz4fXgghhListOxS3tgUT3VtA3MmBXNTiPkUN8gKa0IIITqYM+cbi7u2TsPcySEMDvY0dqRrJuUthBCiw0jNusibmxOoq9Pw2JQQBgZ6GDvSdZHyFkII0SGcyizhf5uPU9+gYd6UEAaYaXGDlLcQQogO4FRmCW9uTqChQcu8KX3p38dd9xeZMClvIYQQ7drJjBLe2tJY3I/f2ZeIAPMubpDyFkII0Y4ln7vAii3H0Wi1zL8rlPDebsaO1CakvIUQQrRLSekXWLH1OFqtliemhhLm3z6KG6S8hRBCtEMnzhazYmsiAE/eHUZoz/a1IJiUtxBCiHbleFoxq2ISURR4cmoofdtZcYOUtxBCiHYk4UwRb3+RiKIoLLgnjBBfF2NH0gspbyGEEO1C/OnG4rZQKSy8J4ygdlrcIOUthBCiHYhLLeSdL09gYaGw6J5+BPZwNnYkvZLyFkIIYdaOnirgva+SUFuoWHRvGH182ndxg5S3EEIIM3YkpbG4LS1VPHVvPwK6dzZ2JIOQ8hZCCGGWYk/m8/62ZKwsVTx1Xz96d+sYxQ1S3kIIIczQr8n5rNneWNxP3x9OL28nY0cyKClvIYQQZuVQUh5rdiRjY2XB0/eH4+/VsYobpLyFEEKYkYMn8vjg62RsrNQ8c384Pb06GTuSUUh5CyGEMAu/JOby0dcnsbVW88y0cPy6dsziBilvIYQQZmB/Qg5rv0nBzkbNs9Mi6NHF0diRjErKWwghhEn76bfitrdR848HIvDx7NjFDVLeQgghTNi++GzW7TqFg60lz04Ll+L+jZS3EEIIk/TDsfN8uicVRztL/jEtgm4eDsaOZDKkvIUQQpic746eZ/23qXSys+QfD0Tg7S7F/WdS3kIIIUzK3iNZfL73NJ3srRqL283e2JFMjl7Le+nSpSQkJKAoClFRUYSFhTXdl5uby9NPP01dXR3BwcG8/PLL+owihBDCDOw5nMXG707jZG/F4ukRdHWV4r4clb42HBsbS0ZGBtHR0bzyyiu88sorze7/73//y+zZs9myZQsWFhbk5ORccXuni9Oprq/RV1whhBBGVF1bz/o9qY3F7SDFrYveyvvgwYOMGTMGAH9/f0pLS6moqABAo9Fw9OhRRo0aBcCSJUvw8vK64vb+uXcZy46skAIXQoh2prq2nufeO8h3x86jUik8dW8/KW4d9HbYvKioiJCQkKbbLi4uFBYW4uDgwIULF7C3t+c///kPSUlJDBgwgGeeeUbnNvMrC6m0LKW7m7++Ypsdd3d528RfyUyak3m0JDNpyZgzeWdrAuWVdQBoNFocHG2M/ndk7MfXxWAnrGm12mZ/zs/PZ+bMmXh7ezN37lz27dvHyJEjdW5nZ9I+Oge5oVL0dtDAbLi7O1JYWG7sGCZFZtKczKMlmUlLxpzJ9gPn+ObAOSxUCg0aLV1d7bBTK0b9OzKl50hrv0Torbw9PDwoKipqul1QUIC7uzsAzs7OeHl54ePjA8DNN9/M6dOnr1jeUcOfZO3RzRzMO4JKpWJan6lS4EIIYca2/ZLOl/vTce1kw8J7w6ipa8DbzR4bK3kjlC56a78hQ4awe/duAJKSkvDw8MDBofF9emq1mu7du3Pu3Lmm+/38/K64vfCuwSyKnEd3By9+yYllQ8pWNFqNvuILIYTQo69+bixuNycbnpseQTd3B/y9nKS4r5LephQZGUlISAjTpk1DURSWLFlCTEwMjo6OjB07lqioKJ5//nm0Wi0BAQFNJ69dib2lHU9GzGVV/BoO5B5Go9XyYNA9sgcuhBBmQqvV8uX+dLYfOId7ZxsWPxCJq5ONsWOZHUX75xejTdzvr0FU1lWyKv5DMsqzGNylPzOC7u2QBW5Kr8uYCplJczKPlmQmLRlqJlqtlpifzvL1wQw8OtuyeHoELp1Mr7hN6TnS2mveZtl4dpZ2PBH+KL6dfPg17yjrkjfJIXQhhDBhWq2WrT82Frens+kWt7kwy/IGsLO05YnwR/Dr5MPh/GN8kryRBk2DsWMJIYT4C61Wy+Z9aew8lIGnix2Lp0dKcd8gsy1vAFu1LfPDH6WnUw+O5MdLgQshhInRarVEf3+GXb9m0sXFjuemR+DsaG3sWGbPrMsbwFZtw/x+j+Dv5MvRggQ+Tt4gBS6EECZAq9Wy4bvT7DmcRVfXxuLu7CDF3RbMvrwBbNQ2PN7vEXp19iOu4DgfJX0uBS6EEEak1Wr5fO9p9h45j7ebPYunR+Ikxd1m2kV5A9iorXm83yP07tyT+MJEPjzxGfWaemPHEkKIDker1bL+21S+O3qebu72/OOBCJzsrYwdq11pN+UNYG1hxeP9ZhPg3IuEoiQ+kAIXQgiD0mi1fLonle+PZdPN3YFnH4igkxR3m2tX5Q1gZWHF38MeJtC5N4lFyaxJ/JQ6KXAhhNA7jVbLul2n2BeXjY+HA/94IJxOdlLc+tDuyhsaC/yxsIcJcgngRPFJ1iSuo66hztixhBCi3dJotXzyTQo/JeTg49m4x+0oxa037bK8AawsLHksdBbBLn1IKk7hfSlwIYTQC41Gy9qdKew/nkuPLo7844EIHGwtjR2rXWu35Q1gaWHJ3NCZhLgGknzhFKsTP6FWClwIIdqMRqPlo50n+TkxF7+ujvxjWjj2NlLc+tauyxsaC3xO6Ez6ugZx8kIqq4+vpbah1tixhBDC7Gk0Wj78OpkDJ/Lo6dWJZ+4Px06K2yB0lvdTTz3FgQMHDJFFbyxVah4NfYhQt2BSSk7z7vG11EiBCyHEdWvQaPhgRzIHk/Lx9+rE0/dJcRuSzvIeO3YsGzduZOLEibz99tvk5eUZIlebs1SpebTvDPq59yW15AzvJnwkBS6EENehQaNhzfZkDiXn08vbiafvD8fORq7DbUg6y/v2229nxYoVbN68GV9fX55++mnmzJljlnvjapWaR0IeJNw9lNMXz/JOwodU19cYO5YQQpiN+gYNq7clE3uygN7dnHjqvn7YWktxG9pVveZdVVXFnj172Lx5MxqNhltvvZV169bx5ptv6jtfm7NQWTA7ZDoRHmGcuZj+W4FXGzuWEEKYvMbiTuJISgEB3TtLcRuRzvJ+4YUXuO2220hISOD5559n48aNTJ8+nXfffZf9+/cbImObs1BZ8LfgB+jv0Y+00nO8nfAhVVLgQgjRqvoGDe99lcTRU4UE+nTmqXv7YWMlxW0sOicfFBTEiy++iJ2dXdPH4uPjCQ8P59///rdew+mThcqCWcHTUBSFI/nxvB3/IfPDZ2OrtjV2NCGEMCn1DRre/fIEcaeLCOrhzIJ7wrC2tDB2rA6t1T3vsrIyMjMz+frrrykuLiYrK4usrCzOnj3Lc889B0BwcLDBguqDhcqCmUH3M9AzgvSyDFbFf0hVfZWxYwkhhMmoq9fwdkwicaeLCPaV4jYVre55x8XF8cknn3Dy5ElmzZrV9HGVSsXQoUMNEs4QLFQWzAy+H5Wi4te8o6yM+4Anwh/FzlL2wIUQHVtdfQNvf3GC42nFhPi58OTUUKykuE1Cq+U9YsQIRowYwYYNG3jggQcMmcngVIqKGUH3oqBwKO8IK+Pf58nwOdhZ2un+YiGEaIfq6htYGZPIibMX6Nuzsbgt1VLcpqLV8t66dSt33303+fn5vPXWWy3uX7hwoV6DGZpKUfFg0D2oFIUDuYdZEb+GJ8PnYC8FLoToYGrrGos7Kf0CYf6uzL+rrxS3iWn1NW+VqvEutVqNhYVFi//aI5Wi4oHAuxniNYis8mxWxL1PRd0lY8cSQgiDqalrYMXW4ySlX6Cfvyvz75I9blPU6p73lClT0Gg0PP7444bMY3QqRcW0PlNRFBU/Zx9iRdz7LAifi4OVvbGjCSGEXtXUNbBiy3FOZpQQ3suNv9/ZF0t1u78EhllqtbyDg4NRFKXFx7VaLYqicPLkSb0GMyaVomJawF2oUPgp+yBvxa1mQcRcHK0cjB1NCCH0oqa2gbe2JJCSeZHIAHfmTQlBbSHFbapaLe+UlJQb3vjSpUtJSEhAURSioqIICwtrum/UqFF06dKl6RD8a6+9hqen5w0/ZltRFIX7Au5EUVT8eP6XpgLvZOVo7GhCCNGmqmrqeXNzAqlZF+nfx53H7pDiNnU6T1i73MlqoPuEtdjYWDIyMoiOjiYtLY2oqCiio6Obfc6aNWuwtzfdw9GKonBv7ztQofDD+Z9569hqHgv7G5fqLtHV3hMbtbWxIwohxA0pKa/h5U+OcC63jAGBHsydHCzFbQZaLe/fT1i73pPTDh48yJgxYwDw9/entLSUiooKHBzM69Czoijc3XsyiqLwfdZ+/v3razRoG/C082DxgCelwIUQZqukvIbnVx+krl6DjZUFs8YHSHGbiVbL+6677gLgiSeeoKysjHPnzqEoCn5+fldVwEVFRYSEhDTddnFxobCwsNnXLlmyhOzsbPr3788zzzxz2dfY/8zd3XiHrB9zf4A6VQ37M2IByK8soNqqnO6ubkbLBMadiamSmTQn82hJZgKXqur419rD1NVrAKiubaBGo+ArswFM/zmic23ztWvX8u677+Lr64tGo+H8+fM88cQTPPjgg9f0QFqtttntBQsWMGzYMJycnJg/fz67d+9m/PjxV9xGYWH5NT1mW7vTZxLHsk9wqb4SC8WCmgothRrjZXJ3dzT6TEyNzKQ5mUdLMhOorK7njU3xZOSVY2tlQVVtA11d7bBTKx1+NmBaz5HWfonQWd5bt25l7969ODo2bqC0tJQZM2boLG8PDw+KioqabhcUFODu7t50+84772z68/Dhw0lNTdVZ3sZmY2nDyzc/z6bTX/Fr3lHeSfiIRRGP4WzT2djRhBDiqlRW1/F6dDzpueXc0rcL08f0ploDdmpFrhJmRnS+uOHp6dlU3ABOTk74+Pjo3PCQIUPYvXs3AElJSXh4eDQdMi8vL+eRRx6htrYWgMOHD9O7d+/r+gYMzcbShoeC7mOC72iKqor537H3uFBdYuxYQgih06XqOpZvbCzuIaFdmH17EHY2lvTp4SLFbWZa/dvasmULAF5eXsybN49bbrkFlUrFoUOHruotXZGRkYSEhDBtWuNlN5csWUJMTAyOjo6MHTuW4cOHc//992NtbU1wcLDJ73X/maIoTOp5G4qiYmf6t/zv2GoWRjyGq62zsaMJIcRlVVTV8frGeDLyyxkW1pVZEwJR6TjPSJguRfvXF6N/88ILL1zxC//zn//oJdCVmMprEH+2M/1bvk7/FhcbZxZGPIabrYvBHtuUXpcxFTKT5mQeLXXEmVRU1fHahjgyCyoY3s+LmeP7NCvujjiTKzGleVzza95XKud169bdeKJ24na/sagUFdvP7uZ/x95jUeRjuNm6GjuWEEIAUFZZy2sb4jlfWMHICG9mjAuQPe52QOeLHCdPnuS9996jpKTxdd3a2lry8vKYOXOm3sOZi/G+o1Gh4quz3/DmsfdYFDEPdzspcCGEcZVdqmX5xjiyCy9xa6Q3M8YG6HxLrjAPOk9Y+9e//sW4ceMoLS1l9uzZ+Pr6smzZMkNkMyvjfG/lTv/buVhTyv/i3qOgstDYkYQQHVjppVqWb2gs7tH9u0lxtzM6y9vGxoaJEyfi6OjIyJEjeeWVV/jwww8Nkc3sjO0xkrt6TWws8GOryb9UYOxIQogOqLSihmWfHyO76BJjBnRj+pjeUtztjM7yrqmpITU1FWtra2JjYyktLSU7O9sQ2czSGJ8R3N1rEqW1ZfwvbjV5UuBCCAO6WFHDsg1x5BZXMm5gdx4YLcXdHuks72effZbMzEwWLFjAiy++yLhx45g8ebIhspmtUT7Duaf3HZTVlvO/uPfIu5Rv7EhCiA6gpLyGVz9vLO7xg3y4f1QvKe52SucJa/3792/68++Lrgjdbu0+FJWiYlPql/zvWOPlRL0cuhg7lhCinSopbzxUnl9SxYSbfLhnhL8Udzumc887JSWFqVOnNi2i8s4775CQkKD3YO3BiG63cH/AXZTXVfBW3GqyK3KNHUkI0Q5dKKvm1fWNxT3x5h5S3B2AzvJ++eWXWbp0adO65BMmTDDKAi3mani3m3mgz1Qq6i6xIu59zpfnGDuSEKIdKS6t5tXPj1FwsYrJt/gydXhPKe4OQGd5q9VqAgMDm277+fmhVssauNdiqPdNTA+8m0t1layIe58sKXAhRBsouljFq58fo/BiNXcM8eUuKe4O46rKOysrq+kJ8eOPP7a4vKfQbYjXYKYH3kNlfRUr4laTWX7e2JGEEGas8GIVr34eR1FpNXcO9ePOYT2NHUkYkM5d6Oeee47HH3+c9PR0+vfvj7e3N6+++qohsrU7t3gNRKUofHZyMyvi1rAgfA4+nboZO5YQwswUXKxi+efHKC6r4a7hPZl8i6+xIwkD01neffr0Yfv27Vy4cAErK6umy3qK63NT1wEoKHx6chMr4t/nyfA59OjU3dixhBBmIr+kkmWfx1FSXsPdI3oy8WZfY0cSRqCzvE+dOkVMTAzl5eXNDpfLSWvXb3DX/qgUFZ8kb2RF3BqeCH8UPyfd10gXQnRs+RcqWbahsbjvHenPhJt6GDuSMBKd5b1w4UImTZpEr169DJGnwxjYJQJFUVibtIFV8WuYH/4oPZ3kB1EIcXm5xZdYtiGO0opa7ru1F+MHyy/8HZnO8vb29uaJJ54wRJYOZ4BnOAoKa5M38Hb8B8wPf4SeTr7GjiWEMDG5xZdY9nkcpZdqmTaqF+MGSXF3dDrLe8qUKbzzzjtEREQ0e4vYwIED9Rqso+jv2Q9FUfg46XNWxX/A4/0eoVdnP2PHEkKYiOyiSyzfEEfZpVoeGN2bsQPlHBlxFeW9bds20tPT+fnnn5s+pigK69ev12uwjiTSIwyVouLDE5/xdsKHPB42m97O8rYPITq684UVvLYhjrLKOh4cG8Do/vLuFNFIZ3lfuHCB7777zhBZOrRw97482vchPjzxGe8kfMjf+80mwNnf2LGEEEZyvqCC5RvjKK+s46FxAdwaKcUt/qBzkZaBAweSmZlpiCwdXj/3EOaEPkSDVsM7CR9x6sIZY776ouEAACAASURBVEcSQhhBZn45yzY0FvfM2/pIcYsWdJb3L7/8wu23387QoUMZOXIkI0aMYOTIkQaI1jGFugUzN3QmWq2Gd49/RMqF08aOJIQwoMz8cpZviKOiqo5Z4/swMsLb2JGECdJ52Hz16tWGyCH+pK9bEHPDZvF+4jreO/4xj4U+TJBrgLFjCSH0LCOvnNc2xlFZXc/fJgQyrJ+XsSMJE6Vzz3v58uV4e3u3+E/oV4hrII+FzkILvJe4lqTiU8aOJITQo/TcMpZv+K24bw+S4hZXpLO8u3XrxpYtW0hLSyMrK6vpP6F/wa59mBf2MArw/vG1nCg6aexIQgg9SM8t47WN8VTV1DN7YhBDw7oaO5IwcToPm+/cubPFxxRFkTPQDSTIJYB5YX/jveNrWZO4jkdDHyLULdjYsYQQbSQtp5Q3ouOprm3g0cnB3BzSxdiRhBnQWd7ff//9dW986dKlJCQkoCgKUVFRhIWFtfic119/nfj4eD799NPrfpz2LtClN4/3+xvvJHzMmsRPeaTvDPq5hxg7lhDiBp3JbizumroG5kwO5qZgKW5xdXQeNs/OzmbBggU89NBDAGzevJlz587p3HBsbCwZGRlER0fzyiuv8Morr7T4nDNnznD48OFrT90BBTj3Yn6/2VgoKj448SmH8+I4XZxOdX2NsaMJIa5RdW09++KzeW1jHLV1Gh67I0SKW1wTneX94osvMmXKlKYrivn6+vLiiy/q3PDBgwcZM2YMAP7+/pSWllJRUdHsc/773//y1FNPXU/uDqm3sz/zwx/FQrFgbfIG/rl3GcuOrJQCF8KMVNfW8+IHv7Ju1ylq6zTMnhjIoCBPY8cSZkbnYfO6ujpGjx7N2rVrgatf07yoqIiQkD8O7bq4uFBYWNh0PfCYmBgGDRp0TWeuu7s7XvXntlfu7mFUKPex5mjj8rT5lQVUW5XT3dXNyMlMhzxPmpN5tGTMmew8kE5x2R+/cAf2dDOJvyNTyGBKTH0eOssboKysDEVRADh9+jQ1Nde+p/fna4FfvHiRmJgYPv74Y/Lz8696G4WF5df8uO1RoH0gLjaduVB9EYCU8+l01kh5Q+MPnDxP/iDzaMmYM0nJKOGDrxKbbnd1tcNOrRj970ieJ82Z0jxa+yVCZ3nPnz+f++67j8LCQiZPnkxJSQnLly/X+YAeHh4UFRU13S4oKMDd3R2AQ4cOceHCBR588EFqa2vJzMxk6dKlREVFXe3306HZqK3556BnSK1MYe2xLXyWsgULlZqBXSKMHU0I0YqT5y7w1pbjaDTw2JRg3Jxs8Xazx8bqqvahhGhG57Pmpptu4ssvvyQ1NRUrKyv8/PywtrbWueEhQ4awcuVKpk2bRlJSEh4eHk2HzMePH8/48eMBOH/+PC+88IIU9zWyUVsz2n8onXBmVfwHfJK8ES1aBnWJNHY0IcRfJJ27wIotx9FqtcyfGkp4LzlSJm5Mq+Xd0NDAxo0bSU9PZ8CAAU1le7UiIyMJCQlh2rRpKIrCkiVLiImJwdHRkbFjx95wcNHIt5MPT4bPYWX8B6xLjkar1TK4a39jxxJC/OZEejErtyai1cITU0MJ85fiFjdO0f75xeg/efnllykpKaF///788MMPDBs2jIcfftjA8ZozldcgTMWfX5fJLD/Pyrg1VNVX82DgPdzsdXUnFrY3pvRalSmQebRkyJkknm0sboAFd4fSt6erQR73WsnzpDlTmkdrr3m3+lax5ORk3nzzTWbMmME777xzQ4u1CP3zcezGgoi52KltWZ+yhQM5scaOJESHdjytiJVbj6MosOAe0y1uYZ5aLW9LS8umP1tbW9PKDrowId0dvRsL3LKxwH/OPmTsSEJ0SPFnilgVk4hKUVhwTxh9/aS4Rdtqtbx/f2tYa7eFaerm6MXCiMdwsLRnw6kY9mcfNHYkITqUuNRC3v6tuBfeE0aIr4uxI4l2qNUT1pKTk3nwwQebbp86darZ7fXr1+s3mbhu3g5dWRjxGG/FrWbjqS/QaLWM6HaLsWMJ0e4dPVXIe1+dwMJC4al7+9HHx9nYkUQ71Wp5v/POO4bMIdqYl0MXFkXO46241WxK/RKtVsvI7kOMHUuIdutISgGrtyWhtlCx6N4wKW6hV62W96BBgwyZQ+hBV3tPFkU8xltx77P59Fdo0DCq+zBjxxKi3TmcUsDqr5KwtFTx1L39COje2diRRDun88Ikwrx1+a3Anawc2Xp6O99l/mTsSEK0K7En81n9VRJWliqeuS9cilsYhJR3B+Bp78HCyHk4WXUi5swOvs3YZ+xIQrQLh5LyWL0tCWsrFc/cH06vbk7GjiQ6iKsq79TUVPbu3Qs0XqREmB9PO3cWRc6js7UTX6btZE/GD8aOJIRZO3gijzU7krGxUvPM/RH4e0txC8PRWd5r164lKiqKFStWAI0nssnJbObJw86NRRHzcLbuzFdp37DrnCy8I8T1+CUxlw92JGNrpebZaeH09Opk7Eiig9FZ3jt27GDTpk04OTX+Vrl48WL27dun71xCT9ztXFkU+RjO1p3ZfnYX36TvNXYkIczKz8dz+ejrk9jZqHn2gXD8ukpxC8PTWd729vaoVH98mkqlanZbmB83W1eeipyHq40zO9L38HX6t8aOJIRZ2J+Qw8c7fyvuaRH4dpHiFsahs4V9fHxYtWoVZWVl7Nmzh0WLFuHv72+IbEKPXG1dWBgxD1cbF3amf8uOs3tkCVwhruDH+Gw+/iYFe1tL/vFABD26XP6CEUIYgs7yfumll7C1tcXT05Nt27bRr18/lixZYohsQs9cbZ1ZFPkYbjYufHNuLzvO7pYCF+IyfojL5pNdp3D4rbh9PKW4hXG1ukjL71asWMGUKVN45JFHDJFHGJiLjXPTSmy7Mr5Hg5Y7eo6XteyF+M33x87z2Z5UHO0ai7ubu4OxIwmhe8/bzs6Op556iqlTp7J27VqKiooMkUsYkLNNZxZFzsPDzo09GT/wZdpO2QMXAvjuaGNxd7KzZLEUtzAhOsv773//O9u3b2f58uWUl5czd+5c5syZY4hswoA6WzuxKGIennbu7M38kS/OfC0FLjq0bw9nsf7bVJzsrVg8PRJvKW5hQq76tHFra2tsbW2xtbWlqqpKn5mEkThZd2JhxDw87Tz4Lusntp7ZLgUuOqQ9sZls+O40Tg5WLJ4egZebvbEjCdGMzte8V69eze7du6mrq2PSpEm8+uqrdOvWzRDZhBE4WTuyKLLxYiY/ZP2MVqvlnt53yGvgosPY9Wsmm344Q2eHxj3uLi52xo4kRAs6y7u0tJSlS5cSGBhoiDzCBHSycmRRxGOsiHuffed/QaPVcl/AFClw0e7tPJTBln1pODtas/iBCDyluIWJarW8t27dyt13342VlRW7d+9m9+7dze5fuHCh3sMJ43G0cmBBxFxWxq/hp+wDaGkscJUiC/SI9mnHgXPE/HQWl06Nxe3hLMUtTFer/xL/voqaWq3GwsKixX+i/XO0cmBB+Fy8HbqyP/sgG099gUarMXYsIdrc9l/SifnpLK6drFk8PVKKW5i8Vve877rrLgAcHBx4+OGHm933+0VKRPvnYGXPgoi5rIpbwy85v6LVankgcKrsgYt246uf0/nq53RcO9nw3PQI3DrbGjuSEDq1Wt6HDh3i0KFDbNu2jdLS0qaP19fXExMTw4IFCwwSUBifg6U9T0bMZVX8Gg7kxqJBw4OB90iBC7Om1Wr56ud0tv1yDjcnGxZPj8DNSYpbmIdW//Xt2bNn0xrmfz5cbmNjwxtvvGGwgMI02Fva8WT4HHo4dudQ7hE+O7lZDqELs6XVavlif2Nxu3e24bnpkVLcwqy0uuft4eHB5MmTiYiIaPHWsHXr1jF48GCdG1+6dCkJCQkoikJUVBRhYWFN923atIktW7agUqkIDAxkyZIlcjazibOztOOJ8EdZlfABv+YdRYuWh4Lukz1wYVa0Wi0xP53l64MZeHS2ZfH0CFw62Rg7lhDXROdbxcrLy1m4cCElJSUA1NbWkpeXx8yZM6/4dbGxsWRkZBAdHU1aWhpRUVFER0cDUFVVxddff8369euxtLRk5syZxMXFERkZ2QbfktAnO0tbngx/lLfjPyQ27xh1DXXc2n0o3g5e2KitjR1PiCuqqqnjzc+P8cOx83g627J4eiTOjvK8FeZH5y7Tv/71L8aNG0dpaSmzZ8/G19eXZcuW6dzwwYMHGTNmDAD+/v6UlpZSUVEBgK2tLZ988gmWlpZUVVVRUVGBu7v7DX4rwlBs1bbMD38U307diStM5I1j7/LqkRVU19cYO5oQraqsrmPxuwf54dh5LFQKC+/pJ8UtzJbOPW8bGxsmTpzIhg0bGDlyJMOGDePxxx9n0KBBV/y6oqIiQkJCmm67uLhQWFiIg8Mf6wO///77rFu3jpkzZ9K9e3edYd3d5TJ8f2W8mTgyU3U3L+/7HwAFlYWcrjrFqJ5DjJTnD/I8aU7mARWVtby6IY5L1fUANGi0WNlaymz+RGbRnKnPQ2d519TUkJqairW1NbGxsfTq1Yvs7OxrfqDLrZE9d+5cZs6cyZw5c+jfvz/9+/e/4jYKC8uv+XHbM3d3R6POxEnjioedOwWVhQCsPryejMJcJviOwUJlnLUAjD0TUyPzgKyCClbFHKfwYjXWlipq6jR0dbXDTq10+Nn8Tp4nzZnSPFr7JUJneT/77LNkZmayYMECFi9eTHFxMY8++qjOB/Tw8Gh2+dCCgoKmQ+MXL17k9OnTDBw4EBsbG4YPH86xY8d0lrcwLTZqa54bsIDcS/lU1lWy4VQM35z7jqTiFGYFT6OLvaexI4oO7lByHmu/SaG2TsOkW3owfpAP1RqwUyvYWOn8508Ik6Xz2fvnQv3rEqlXMmTIEFauXMm0adNISkrCw8Oj6ZB5fX09zz//PNu2bcPe3p7ExETuuOOO64gvjM1GbY2fkw8A/+z8FJtTt/Fr3lH+e/gtpvjfzohut8jZ6MLgGjQaNv+Qxp7DWdhYWfDE1FAiAxp3HnqY0F6VENer1fKePn36Fd+6tX79+ituODIykpCQEKZNm4aiKCxZsoSYmBgcHR0ZO3Ys8+fPZ+bMmajVavr06cPo0aOv/7sQJsFWbcvM4PsJcw9hQ8pWtpzexomik8wIuhdnm87Gjic6iLJLtbz31QlSMi/S1dWOJ6aG0tVVLukp2hdF28oFm2NjY6/4hbpOWNMH+W25OVN6XeavSmvK+TxlMyeKU7BV23J/wJ0M8AzX+3v5TXkmxtDR5nE2p4y3v0ikpLyG/gHuzJ4YhK11832UjjaTqyEzac6U5nHNr3n/uZz37dvH+fPnmTFjBpmZmVd1Zrjo2JysHZkX9jcO5MSy5cx21iZv4HhREtP6TMXeUi76INreTwk5fLbnFA0NWu4e0ZPbb+ohCz+Jdkvna97Lly8nIyODnJwcZsyYwfbt27lw4QIvvviiIfIJM6YoCkO8BxPg3It1JzdyrOA4aRfPMSPoXoJd+xg7nmgn6uo1fL43lR/jc7C3UfPYPSH09XM1diwh9ErnmUSHDx9m1apV2Ns3vmY0f/58kpKS9B5MtB/udq48Ffl37ug5noq6S7yd8CHRp76ktqHW2NGEmbtQVs2rnx/jx/gcfDwceOnhgVLcokPQuedtbd24AtHvh58aGhpoaGjQbyrR7qgUFbf5jiLYNZBPkjfwU/YBUkpSmRU8Dd9OPsaOJ8zQqcwS3v3yBGWVddwc4snM8YFYWxpnfQEhDE3nnndkZCQvvPACBQUFfPzxx8yYMcMoJ6uJ9qG7oxfPDVjAqO7DKKws5vWj77Dj7B4aNPILobg6Wq2Wb49k8drGeC5V1zN9TG8enRQsxS06lFbPNv+zXbt28euvv2JlZUX//v0ZN26cIbK1YCpn/5kKUzoj8nqklqSxLjmakpqL+Dh2+21hF48b2qa5z6Sttbd51NQ1sG5XCgeT8ulkZ8nf7+xLHx/na9pGe5tJW5CZNGdK82jtbHOd5V1SUoKzc/MfjvPnz7e4TKghmMowTYUpPcGuV1V9VdPCLpYqNXf6T2R4t5uve2GX9jCTttSe5lF4sYq3YxLJLKjA36sTj98Vel0XFmlPM2krMpPmTGkerZV3q/9CHjlyhGHDhnHbbbcxfvx4MjMzAfjss8+YPn26flKKDuf3hV0e7fsQVhZWbD79FW/Hf8jFmlJjRxMm5MTZYl5ee5jMggpGhnvJpTxFh9fqCWtvvvkma9euxd/fn++++44XX3wRjUaDk5MTmzdvNmRG0QFEeITS08mX9SmbSSpO4d+/vsG0gDsZ0CXC2NGEEWm1WnYeyiDmx7NYWCg8PCGQ4f28jB1LCKNrdc9bpVLh7+8PwOjRo8nOzmbmzJmsWrUKT0+54IRoe07Wjvw97G880GcqDdoGPk7ewEcn1nOprtLY0YQRVNXU8/YXJ9j641k6O1rz/IP9pbiF+E2re95/XZmoa9eujB07Vu+BRMemKApDvW9qXNglOZqjBQmcuZjOQ0H3EeQaYOx4wkByiy+xKiaR3OJKAn06M29KXzrZWxk7lhAm46rPCpJlBoUhedi58VTkPCb3HE95XQWrEj6QhV06iGOphfzfJ0fILa5k3MDuPDMtXIpbiL9odc87Li6OkSNHNt0uLi5m5MiRaLVaFEVh3759BognOjILlQXjfUcR7BrAJ8nRTQu7PBz8AD06yfr67Y1Go+WL/Wf5+mAGVmoVc+8I5qbgLsaOJYRJarW8d+3aZcgcQrTKx7Ebzw9YwLazu/g+az+vHX2b8T1GMd53NBYqWZijPaioquP97UmcOHsB9842PDE1jO4eDsaOJYTJarW8vb29DZlDiCuytLDk7t6TCXULYl3yJnae20tS8SlmBd+P5w0u7CKMKzO/nLe/SKTwYjWhPV2Ze0cw9jaWxo4lhEm7vpUwhDCSAOdeRA16ikFdIskoz+I/h99i3/lfuIqFAoUJOpSUx9JPj1J4sZrJt/iy8J4wKW4hroLOC5MIYWrsLG2ZFTyNULdgNqbEsDn1K04UneTegClcLC7Cpt4RG7Us4GHK6hs0bP4hjW+PZGFjZcGTU0OJCHA3diwhzIaUtzBbkR5h+Dv58lnKZpKLT/F/h15DixZPO3cWD1ggBW6iSi/V8t6XJziVdZGurnY8MTWUrq72xo4lhFmRw+bCrDlZd+LxsNmM9RmJlsZD5/mVhRwrSDByMnE5aTmlvLz2MKeyLtK/jzv/b+YAKW4hroPseQuzpygK431HE1+YSGFVMQDrU7aQXprBHf4TcLSSs5aNrbq2nh0HzrE7NhONFu4e0ZPbb+oh60cIcZ2kvEW7YKO25vmBi6i2Kie7oIgv03ZyIPcwcYUnmNzzNoZ6DZa3lRlJeWUtz68+RFVNPYoC8+/qS2SAvENAiBsh5S3aDRu1Nd1d3eiscSPQpTf7sw+xI303m1K/5JecX7kv4E56dfYzdswOJSOv8W1gVTX1AGi14OQg5yIIcaOkvEW7ZKGyYGT3IfT37MeXaTs5lHuEN4+9y0DPCO7sdTudrZ2MHbFda9Bo2Hkok20/p9Og0WJnraaypp6urnZ4u8lr3ELcKClv0a45WjnwUNB9DPUazKbULzmcH8fxoiRu9xvLyG5DUKvkR6Ct5ZdU8sGOZNKyy+jsYMXsiUH08nYiu+gS3m722FjJzIW4UfJTJDoEP6ce/GPAkxzMOcxXZ7/hizNfcyDnMPcFTCHQpbex47ULWq2WH+NziP7+DDV1DQwK8mDGuD442DYuuuLvJUc7hGgrei3vpUuXkpCQgKIoREVFERYW1nTfoUOHeOONN1CpVPj5+fHKK6+gUsk714T+qBQVQ7wHE+4Ryo6zu9mffYiV8WsIdw/l7t6TcLFxNnZEs1VaUcPH36RwPK0YO2u1XFRECD3TW3nHxsaSkZFBdHQ0aWlpREVFER0d3XT/Sy+9xLp16+jSpQsLFixg//79jBgxQl9xhGhib2nH/X3u4havQWxK/ZL4wkSSilO4rccoxvgMx9JClue8FkdSCli3+xQVVXWE+Drzt9uDcOlkY+xYQrRreivvgwcPMmbMGAD8/f0pLS2loqICB4fG99zGxMQ0/dnFxYWSkhJ9RRHisro7evN05OPE5h3ji7Sv2ZG+m0O5h7kn4A5C3YKNHc/kVVbX8/neVA6cyMNSrWL6mN6M6t8Nlbx3Wwi901t5FxUVERIS0nTbxcWFwsLCpsL+/f8FBQX88ssvLFy4UOc23d0d9RPWjMlMWrrWmUzyGMmowMFsSfqanad/4L3ja4ns2peHI+6li6P5vx9ZH8+RxDNFvLnxGIUlVfTq3pmnH4iku6f5PBfl56YlmUlzpj4Pg52wdrmrPhUXFzNv3jyWLFmCs7Pu1xsLC8v1Ec1subs7ykz+4kZmMqHbbfTr3I/Np7dxLPcEx/NOMsZnBON8R2FtYdXGSQ2jrZ8jdfUNbP3xLHsOZ6FSFO4Y4sukW3xRq8zn51N+blqSmTRnSvNo7ZcIvZW3h4cHRUVFTbcLCgpwd//jqkEVFRXMmTOHRYsWMXToUH3FEOKaeDl0YUH4HOIKE9l6eju7Mr7n17xjTO09iQj30A69nGdGXjkf7Egmu+gSni52zJkUTE+vTsaOJUSHpLfTu4cMGcLu3bsBSEpKwsPDo+lQOcB///tfZs2axfDhw/UVQYjroigKkR5hvHTTP7itxyjKa8v58MRnrIxfQ+6lfGPHMziNRsvXB8/x73VHyC66xKhIb/6/vw2U4hbCiBTt5Y5nt5HXXnuNI0eOoCgKS5YsITk5GUdHR4YOHcrAgQOJiIho+txJkyZx//33X3F7pnIYw1SY0qEdU6GPmRRUFrL59DaSi0+hUlTc2m0oE/zGYKs2/TOqb3QeBSWVfLDjJGeyS3FysOKR24Po29O1DRManvzctCQzac6U5tHaYXO9lndbM5VhmgpTeoKZCn3NRKvVcqL4JJtTt1FcfYFOVo7c1WsiAz0jTPpQ+vXOQ6vV8lNCDhu/a1xwZWCgBw/d9seCK+ZMfm5akpk0Z0rzMPhr3kK0J4qiEOoWTKBzb/Zm/sjujO/5JHkjP2cf4t6AO+nu6GXsiG2mtKKGtd+kkPD7giuTgxkc7GnSv6QI0dFIeQtxDSwtLJngN4ZBXSKJObOD+MITvHr4LYZ538yknuOwt7QzdsQbcvRUAZ/salxwJaiHM49MlAVXhDBFUt5CXAdXWxfmhM7kZHEqm09/xU/ZBzhWkMAd/uOJdA8jr7KQrvae2KjN4/KXldX1bNibyi+/LbjywJjejJYFV4QwWVLeQtyAINcAopyf4oesn9l5bi+fp2xlU+pX1Gvq8bTzYPGAJ02+wE9llvDBjmSKy2ro0cWROZOC8ZLLdgph0qS8hbhBapWasT1GMrBLBJ+e3ETKhdMA5FcWcDg/jmHeNxk54eXV1TcQ89NZ9sRmofx5wRULuUCQEKZOfkqFaCOdrZ2Y03cmLtadmz628VQMbx1bzemSNCMmaykzv5yXPznC7tgsPJxteeGhSO4c1lOKWwgzIXveQrQhG7U1/xz8DLmX8qlpqGFv5o+cvJBKalwavTv35Ha/MfTu7G+0M7c1Gi3f/JrBl/vTadBouTXCm/tu7YW1lYVR8gghro+UtxBtzEZtjZ+TDwCBLr1JL81g57m9JBef4q249/F38uV2v7H0ce5l0BIvuFjFBzuSOXO+ccGV2bcHEWrmC64I0VFJeQuhZ35OPZjf7xEyyrLYmb6XE8UnWRm/hp5OPZjgO4YglwC9lrhWq2X/8Vw27D1NTV0DAwI9mNlOFlwRoqOSFdbMmCmtAmQqzGEmmWXn+ebcdxwvSgLAt5MPE3xHE+Ia2KYlXl1bT35ZDTHfnyHxbDG21mpmjAvgpg6+4Io5PEcMTWbSnCnNQ1ZYE8JE+HTqxmNhs8gqz+Gbc3tJKDzBu8c/podjdyb4jaava9ANl2tJeQ3/WnuYsku1APTp3pk5k4NlwRUh2gkpbyGMpLujF3NDZ5Jdkcs36XuJK0zkveNr6e7ozQTfMYS5BV9TiWu1Ws7mlPFjfA6HkvOob/jjoNrdI3tKcQvRjkh5C2Fk3g5deTT0IXIq8th17juOFRzn/cRP8Hboyu2+YwhzD0GltP4WrkvVdRw4kcdPCTlkF14CwLWTNTV1Giqq6ujqakc3d4dWv14IYX6kvIUwEV4OXZjd90EmXBrDrnPfcTQ/gTUnPsXLvgsT/MYQ7t63qcS1Wi2pWRf5KSGHI6cKqavXYKFSGBDowYh+XgT5OlNb10BlvRY7tYKNlfyoC9GeyAlrZsyUTqowFe1pJnmXCth17nuO5MehRUtXe09Gdh1JeY4r+4/nkXehEgBPZ1uGh3sxpG9XOtlbNdtGe5pHW5GZtCQzac6U5iEnrAlhZrrYe/BwyDTG+45ic9IuUsqT2HAmGk2VPRpVLwYHhzGiXzf6+HTu0GePC9ERSXkLYaJKK2r4OTGXnxJyKLzYDcXamU5+WdR1ykDVM4E8u2xKrUaj0YZjocgKaUJ0JFLeQpgQjUbLifQL/JSQQ8KZIho0WqzUKob07cKIcG/8vTtRXF3C7nPfcyjvCOtORvPNub2M9x3NQM8ILFRS4kJ0BFLeQpiAC2XV/Hw8l/3HcyguqwGgu4cDI8K9uCnYEzubP1ZDc7N14cGgexjvO5o9Gd9zMPcIn57cxDfpe7nNdzSDu0RKiQvRzskJa2bMlE6qMBXmNJMGjYbjacX8FJ/D8bPFaLVgbWXB4CBPRoR74dvF8apey75QXcK3Gfs4kBNLvbYBVxtnbusxinD3vtTZVGFT62jy1xQ3JHN6jhiKzKQ5U5pHayesSXmbMVN6gpkKc5hJ0cUqfjqey8/Hc7hY0bgCml9XR4b382JQkCe21td3QKyk+iLfZu7jl5xY6jX1qBQVGq0GD1s3nhu4UAr8SIM51gAAHB9JREFUN+bwHDE0mUlzpjQPOdtcCCOqb9AQf7qIHxNySE6/gBawtbbg1khvRvTzwsfz8j+g18LZpjP3BdzJuB63sjV1B8cKEwAoqCpibdIGJviNxsexm5yZLkQ7IOUthB7lX6jkp4QcfknMpayyDoBe3k6MCPdiQKAH1pZt/9p0Z2snHgy6h6yK8xRWFaNSVCQWJ5NYnIy3Q1eGeA1moGcEdpa2bf7YQgjDkMPmZsyUDu2YClOYSXllLT/EnSc5vYTU86UA2NuouaVvV4b364q3gZYqra6vodqqHOsae86WZXIg51eOFyWj0WqwVFkS6RHGEK/B9HTq0aH2xk3hOWJqZCbNmdI85LC5EHqWU3SJ74+d54e4bH7/lTigmxMjI7zp38cdS7VhzwC3UVvT3dWNwsJyQlz7EOLah9Kacn7NPcIvubH8mneUX/OO0sXekyFegxjUJRIHS3uDZhRCXB+9lvfSpUtJSEhAURSioqIICwtruq+mpoaXXnqJ06dPExMTo88YQuhNTV0DR1IK+DEhhzO/7WX/2b2jeuHv5WSEZJfnZO3ION9bGdNjBKklaRzIiSW+8ARbT2/nq7RvCHfvyxCvwfTu3LND7Y0LYW70Vt6xsbFkZGQQHR1NWloaUVFRREdHN92/bNkygoKCOH36tL4iCKE3mfnl/JSQw8GkfKpq6gEI8XXm5r5d2HEgg7wLlXR1tcPbzTT3ZFWK6v9v796D2yzvBY9/dbeti2PJlu/Oxbk4cRKHYAjBuRwgWbLdsMvOYZJAIVwyZZlJyZbS9tB0WzMspTlTYGZPu9OWne5MS9tZLhNY2sOWHigBGnIncWzHgcSJHd8t2bEl2brr3T/kKFEcO6HEeqX495nxWNL7Wvr5N6/00/O8z/O8VNnnUWWfhzfk40DvET7tPsjhvmMc7juGMzuf20tu5bbiWqxGuSKZEOlmyor3vn37WLduHQCVlZUMDw/j8/mwWOIfBE899RRDQ0O88847UxWCENdVIBThYEs/Hx3r5myPB4Bci5G7bp7J6qUlFMyIDwBbPr+ALvcIpfnmjLial9VoYV3FWu4qX0PrcBt7uw9wtP84b7e+yx/PvMfS/EXUlaxggX3upJcmFUKkzpR9srjdbqqrqxP37XY7LpcrUbwtFgtDQ0Nf6jknOnE/nUlOxrueOVEUhdOdQ7y3v52Pj3biD0bRaqB2YSF33zaTWxYWotONL2jlpXnXLYav6svkw+lcysp5S/GFRvik7SAfnNnLUVcjR12NOM0O7pxTxz/MXok9e8YURjz15H0znuQkWbrnI2XNgusxqD1dRv+li3QaEZkurldORgMRDpzo5aNj3Zzr9wFgt5n4d7dUsHppMXZbFgCDgyNf+bWm0lfJR21eLTcvv5k2zzn2dh/kSN8x/k/jO7ze9CeqHVWsKlnBIseCjGuNy/tmPMlJsnTKR8pHmzudTtxud+J+f38/BQUFU/VyQnxliqLQ2u3h42PdHDzZRygcQ6vRcNO8fNYuK2XxbDta7fQaxKXRaJidO5PZuTP5x3n3cLjvKHu7D9LoPkGj+wQzTLmsLL6F20tuwZ6VPr0NQtzopqx419XV8bOf/YwtW7bQ3NyM0+lMdJkLkU58/jD7mnv5+Fg3Xe54Szo/N4u1y0qoW1LMDIssKwqQrc9idelKVpeu5Jynk73dBzjcd4z/1/Y+f277gIWO+awqWcFix0K5MIoQU2zKivfy5cuprq5my5YtaDQa6uvr2b17N1arlfXr17Njxw56e3s5e/YsDz30EJs2beKee+6ZqnCESKIoCl90DPFRQzeHT7qIRGPotBpuqXKyZlkJC2fmoZWpUhOqsJVRYSvjP8/dyGf9x9nbfYATA59zYuBzbEYrtxXXUldyK/nZDrVDFeKGJCusZbB0Oi+TLq6WE89oiE8be/m4oZvewVEACu05rK0p4fYlRdhyjKkKNSVSeYx0+XrY232Qg72f4Y/4AajKm8ftJbdSlTeXfv8AxeZC1S+QIu+b8SQnydIpH7LCmpi2YopCS/t5Pj7WzWdfuIjGFPQ6LbdVF7K2poT55TNkQZLroNRSzKb5/4l7K7/G0f7j7O0+yMnzpzh5/hQaNCgozDDl8r3aJ8k12dQOV4iMJsVb3LCGfEH2NvbwcUM3rqEAAKX5ZtYsK2FldRGWbIPKEd6YjDoDK4pvZkXxzfSO9PHu2fc50h+/wtlQcJgffvoT5udVsjh/IUscC3Fk21WOWIjMI8Vb3DACoQgtbYO0dZ5nf3Mfx065iSkKRoOWVUuKWbOshMoSm7SyU6jIXMgDVffR6eumb9SFWZ+DzWSlZfALWga/4A3+L8XmQhY7FrI4fyGzbRUy2E2IayDFW2Q8RVE40+Phf7xxHJ8/nHi8otDC2poSViwqIidLDnW1ZOlNfK92Bz0jfYlz3oOB8zQPnKTJ3cLn50/zb+f28G/n9mDW57DIsYDF+QtZZJ9PjiFH7fCFSEvyiSYyUjgS5eS5IRpOu2k4PcCAJ5C0/bH/UMWqJSUqRScul6U3MTu3InHfnpWXmHYWiob4/PxpmtwtNA2c5FDfUQ71HUWr0VKZO4vF+QtZ7FhIYU6B9JoIMUaKt8gY571BGs8M0HDaTXPbIKFwDIAck57aqgJOdw4z5AtR7MihdoFT5WjFtTLqjCzJX8SS/EUoikKnr2eskLdweugsp4bO8NbpfyU/28GSse71uTNmo9fKx5eYvuToF2krpii093oTrev2votTN4odOdTMzaem0sHcslx0Wi2BUITRiEKOXpMRFwQR42k0GsqtJZRbS/j3s+/CE/LSPPA5Te4WWgY/58POv/Fh59/I0pmoss8fa5VXyZXPxLQjn3AirQRCEZrPnqeh1c3x1gE8IyEAdFoN1bPyWDpWsJ1548+FZhn1lJemz/xM8dXZjFZWFteysriWSCzC6aGzNLlbaHSf4JirkWOuRjRomGkrTwx6K7MUS/e6uOFJ8Raqcw35463r1gE+P3eeSDS+bpAtx8CqJcXUzHWwaJadbJMcrtOZXqtPXIP8H+fdQ9+oi6aBFprcLbQOt9HmOcefzr7HDFMuix1VLM5fyIK8uRh1N9bCO0KAFG+hgmgsRmuXJ1Gwu90Xr8xVUWihpjKfmrn5zCq2yhKl4oo0Gg1FZidFZifrKtYyGh7lxOAXNLlbODHwOX/rPsDfug9g0OpZkDc3MegtW5/NqQE3WRGr6iu9CfFVSPEWKTESCNN4ZoDjpwdoPDPASCACgFGvZdncfJbOdbB0jiNxqU0hvowcQw61hcuoLVxGNBblrOdcYtBb08BJmgZOAm+h1+iIKFFyjTaeWPowZdbSjLukqRAga5tntHRaf/dyiqLQMzBKQ2t8sNnpzmFiY4ea3WYaa107qKrIw2i4fotypHNO1CD5ALd/kKaBFg71HqXNcy5pm1FnpMxSQoW1lAprGeXWUgpzCqbdQjFynCRLp3zI2uZiyoUjMb7oGJt73epOLEmqAeaU2hLd4WUFZhlQJFImP9vOP5TVcVtRLf98+F/oH3VhNVhYkDeX7pFezg63c2a4LbG/QWugzFJCubU0XtRtZRTlOKddQRfpTYq3+LspikKne4RDLX109o/Qcu48wVAUgGyTjtoqJzWVDpZUOm64q3WJzJOlN/FPtTsIGL1khS6e8w5FQ3T6eujwdnHO20mHt4t2bwdnPe2JvzVo9ZRYisda5yVUWMsoNhfKXHOhGjnyxDWJKQqu837aer2093lp7/XS1uPBP1asAfJnZLG2poSaSgfzymeg18m5RJFesvQmyh35SV2iRp2RObkzmZM7M/FYOBqma2SsoHu66PB20untpt3TkdhHr9FRYimi3FpGhbWUcmspJZZiDFLQRQrIUSbGiSkKfYOj8QLdGy/U5/q9+IPRpP3sNlNS8f4v/7GaypLcVIcrxHVn0BmYZatglq0CSuOPhWMReny9idb5OW8X3b4eznm72Dv2dzqNjhJzIeVj588rbKWUmosx6OQKduL6kuI9zcViCj0DI7T3XVqofYnub4ifsy5y5FAz18qsQiszi6xUFFrRaOC//+YwPQOjFDtyKM03q/ePCDHFDFo9FbYyKmxliccisQg9I/10eDs55+2iw9tFl6+bDl839MT30Wq0FJsLx86hx4t6QbYdt/984kItQnxZUrynkWgsRo97NFGk2/viLeoLa4QDaDRQ4jAzs8jKzEShtky43OgPH66lyz1Cab5ZliQV045eq08s53r72GPRWJTe0f6xYh5vpXd6u+ny9bC/53DS3xu1BhY7FpKXPYMZRhs2k41co3Xst00Ku5iQfNreoCLRGN3ukaRz1B39PsKRi4Vaq9FQkm9mVlG8SM8sslJeYMFkvPZRtVlGvXSVC3EJnVZHqaWYUksxK4trAYgpMXpH+unwdtHkPslnrgYAQrEwn7mOT/hcJp2RXKONXJMNm9FKrumS22OP55qsZOmyZAbHNCPFO0MFQhE+bx8kR69Bp9XS5fbFW9Nj56k7Xb7EMqMQXxu8ND/eoo4XaxtlBebrOsdaCHFlWo2WEksRJZYiagoW0zXSQ99oP4U5BXxjyVYCkSDDIQ+eoIfhkJfhoIfhkIfhoAdP0IvLP4DCxEtyGLSGeCFPtNrHCr3Rhs10sdDn6LPHFflAJCirzmUgKd4ZZjQQ4WzPMP/rTy14RkLodVpisRixS97Xep2G0gJLvEiPdX2XFVgw6GX0txBqy9Kb+F7tk/SM9F3zOe9oLIon5MWTVNi9eMYK/HDIiyfo4cxw+6RFXq/Vxwu8Md5itxjMHHM14QuPkGfKZdviB8nPdmA25MjKc2lOVlhLU9FYjN5BP10uHx39Pjr7fXS6RhjwBMbtW+zIoaoiL3GeurTAPG2naaXTykjpQPIx3o2ck5gSwxvyJbXah5MKfPy+J+QlpsQmfB6dRofNaMVmsiadi7+8+95iMN+QRT6djhFZYS1NKYqCZyREh8tHZ/8Ina54oe4eGCUSTX5z5ZqNVM+2U2TP5vBJF8MjIYrsOfzw4VoZLCaEQKvRJgorV/7MB+JFfiQ8Sv+om//d/HuGgsNYDGZq8qvxRUbxBD0MBT3xue1Kx4TPo9VosRosiXPvuZcV+gvd9laDRVaou86k5Z1CwXCUbvcInf0+Olw+ulwjdPT78PnDSfsZ9FpK882UOS2UFVgoLzBT6rQkrVIWCEUYjSjk6DVSuC+RTt+Y04HkYzzJSbJAJDhu1bkLFEVhJDI6QSs+ufs+HItM+BoaNFiNlkta72PF/pLbWbosvCEf5ZYSsgzqXqAonY4RaXmnUExRcA/56XRdLNSdrhH6B0fHnY1yzshmXlku5WOFusxpwTkjG6128pGjWUY95aXpc4AJITLTlVadu0Cj0WAxmLEYzJRaiid8DkVR8Ef8FwfbBeNd84liP1bg+0Zd8TnwV4tJZyJLn4VJZyJLZ8KkN2HSGcffnnC7KbHdqDPckF37U1q8X3jhBRoaGtBoNOzcuZOlS5cmtn366ae8/PLL6HQ61qxZw/bt26cylCnj84fpGivOHf0+Osda1MFw8mpk5iw988tnjBVoM2UFFkoLZG60ECLzaTQacgw55BhyKDYXTrifoigEosGkVvtwyEOHp4vD/ccS+1mNlnirPzzCYOA84Vh4wue8amxoMOoMiWJ+efE36Uxk6S/cNqLT6tD2KVgUG1ajZWybMWm/dPgyMGWV4+DBg7S3t/Paa6/R2trKzp07ee211xLbn3/+eX79619TWFjIgw8+yN13383cuXMnfL4L06K+SrGLxRSisRiRqEI0NvYTjSVuR6IxooltV77t84c51TnE8EiYnoERznuDSa+h02ooduQkurzjP2byrCaZhymEmNY0Gg3Z+iyy9VkUmp2JxwORIB2+7rHpc06+V/tkUhd+NBYlGA0RjAYJRoMEokECkWDisfjti9uCkdDF29EgwUj8tj8aYCjkIRQNfaX/w6C98GXAiElvSvpicOHxcb0AE3xpyNKZJhwPcGEa3zzH7HHbpqx479u3j3Xr1gFQWVnJ8PAwPp8Pi8VCR0cHubm5FBfHu2HWrl3Lvn37Ji3e3/mXTzAZdMwptYHCFQpsciGORi8t0vH9rvfJ/RkWI0vmOBIt6fICC0WOnGk70lsIIf4eV5s+p9PqyNFmk2PIvi6vF1NiF78MRC4U+RDnvJ28dfpfE/vdXnwL2frsi18CrvClwecfJRgNTjpF72r0Gl3yl4Cx7v52TweBaJDXN/9i/N/83a92FW63m+rq6sR9u92Oy+XCYrHgcrmw2+1J2zo6Jh7ReEEwHKWl7Xzivk6rQafTotfFFyrR6+L3TUYtOq0hcV+v06DXadHrtMl/o9Oi12rRjW1P/NYm39eP7T8wFODtj1sTr//fHlvBgpn2K4WaMhMNZpjOJCfJJB/jSU7GUz8nVsrJVzWCm8NVHOw7Qpe3l1JrEU+sfOCaBs8pikIwGooP/gsH8EeCBCIB/OH470AkiD889ntsn/jtQHx7OP64PxJgOOwhMBogOslUPkjhgLXrMai9MC+bZx68GXOWHp1Wk/Ju6EAowoHmnsSFOHL0GlUHjKXTiMh0ITlJJvkYT3IynuTkom/ftD0x+t47FMbLlznfrkFHNhaysQAYxn6+pPj5/lFeOvI/6fe7r7jPlBVvp9OJ233xRfv7+ykoKLjitr6+PpxO57jnuNSLO1arPi0qy6iXC3EIIcQNbLLR96mi0WiwGM380y3/lYDxynFM2cnZuro63nvvPQCam5txOp1YLBYAysrK8Pl8dHZ2EolE+PDDD6mrq5v0+RbMtKdFsbxwIY50iEUIIcSNK0tvuuJgNZjClvfy5cuprq5my5YtaDQa6uvr2b17N1arlfXr1/Pss8/y9NNPA/C1r32N2bOvHKAQQgghkskKaxlMzlONJzlJJvkYT3IynuQkWTrlY6KBhDKnSQghhMgwUryFEEKIDCPFWwghhMgwUryFEEKIDCPFWwghhMgwUryFEEKIDCPFWwghhMgwUryFEEKIDJNRi7QIIYQQQlreQgghRMaR4i2EEEJkGCneQgghRIaR4i2EEEJkGCneQgghRIaR4i2EEEJkGL3aAVzJCy+8QENDAxqNhp07d7J06dLEtv379/Pyyy+j1WqZPXs2P/7xj9Fqb+zvIJPl4/XXX+fNN99Eq9VSVVVFfX09Go1GxWhTY7KcXPDSSy9x7NgxXn31VRUiTL3JcnLnnXdSVFSETqcD4MUXX6SwsFCtUFNisnz09PTw7W9/m3A4zKJFi3juuedUjDR1JspJX18f3/nOdxL7dXR08PTTT3PPPfeoFWrKTHac/P73v+edd95Bq9WyePFifvCDH6gY6WWUNHPgwAHl8ccfVxRFUU6fPq1s2rQpafv69euVnp4eRVEU5cknn1T27NmT8hhTabJ8jI6OKlu3blVCoZCiKIry0EMPKUeOHFElzlS62jGiKIpy6tQpZfPmzcqDDz6Y6vBUcbWc3HHHHYrP51MjNFVcLR87duxQ/vKXvyiKoijPPvus0tXVlfIYU+1a3jeKoijhcFjZsmXLtDheJsuJ1+tV7rjjDiUcDiuKoiiPPvqocvToUVXivJK0a7Lu27ePdevWAVBZWcnw8DA+ny+xfffu3RQVFQFgt9s5f/68KnGmymT5yM7O5je/+Q0GgwG/34/P56OgoEDNcFPiascIwK5du3jqqafUCE8V15KT6WSyfMRiMY4cOcKdd94JQH19PSUlJarFmirXeoy89dZb3H333ZjN5lSHmHKT5cRgMGAwGBgdHSUSieD3+8nNzVUz3CRpV7zdbjd5eXmJ+3a7HZfLlbhvsVgA6O/vZ+/evaxduzblMabS1fIB8Morr7B+/Xo2bNhAeXl5qkNMuavlZPfu3dx6662UlpaqEZ4qruU4qa+v5/777+fFF19EucEXVpwsH4ODg5jNZn7yk59w//3389JLL6kVZkpdyzEC8MYbb3DfffelMjTVTJYTk8nE9u3bWbduHXfccQc1NTXMnj1brVDHSbvifbkrfcgMDAzwxBNPUF9fn5T46eBK+Xj88cd5//33+eSTTzhy5IgKUanr0pwMDQ2xe/duHn30URUjUt/lx8mOHTv4/ve/z6uvvsqpU6d47733VIpMHZfmQ1EU+vr62Lp1K7/73e84ceIEe/bsUS84lVzps+To0aPMmTMn0Uiabi7Nic/n41e/+hV//vOf+eCDD2hoaODkyZMqRpcs7Yq30+nE7XYn7vf39yd1Bft8Pr7xjW/wrW99i1WrVqkRYkpNlo+hoSEOHToEQFZWFmvWrOGzzz5TJc5Umiwn+/fvZ3BwkK9//et885vfpLm5mRdeeEGtUFPmau+be++9F4fDgV6vZ82aNXzxxRdqhJkyk+UjLy+PkpISKioq0Ol0rFy5klOnTqkVaspc7RgB2LNnDytXrkx1aKqZLCetra2Ul5djt9sxGo3U1tbS1NSkVqjjpF3xrqurS7QKmpubcTqdSd8Cd+3axcMPP8yaNWvUCjGlJstHJBLhmWeeYWRkBIDGxsa06taZKpPlZMOGDbz77ru8/vrr/PznP6e6upqdO3eqGW5KTJYTr9fLtm3bCIVCABw6dIh58+apFmsqTJYPvV5PeXk5bW1tie3T/X1zQWNjI1VVVWqEp4rJclJaWkprayuBQACApqYmZs2apVao46TdVLHly5dTXV3Nli1b0Gg01NfXs3v3bqxWK6tWreLtt9+mvb2dN998E4CNGzeyefNmlaOeOpPlY/369Wzfvp2tW7ei1+tZsGABd911l9ohT7mr5WQ6ulpO1qxZw+bNmzGZTCxatIgNGzaoHfKUulo+du7cyTPPPIOiKMyfPz8xeO1Gdi3vG5fLhcPhUDnS1LlaTrZt28bWrVvR6XTcdNNN1NbWqh1yglwSVAghhMgwaddtLoQQQojJSfEWQgghMowUbyGEECLDSPEWQgghMowUbyGEECLDpN1UMSGEOj766CNeeeUVtFotfr+fsrIynnvuOWw2m9qhCSEuI1PFhBCEQiFWr17NH//4R5xOJwA//elPcTgcPPbYYypHJ4S4nHSbCyEIBoOMjo7i9/sTj333u9/F5/Px29/+FoAPP/yQmpqaxEpt27Zto7GxUZV4hZjupHgLIbBarTz55JPce++9PPLII/ziF7/gzJkz3H777Rw+fBiAAwcOUFtbS0NDA+FwmLNnz1JdXa1y5EJMT1K8hRBA/Op0f/3rX7nvvvvo7u5m06ZNtLS0JK6kdPz4cR544AEOHjxIY2Mjy5YtQ6uVjxAh1CAD1oQQAPj9fvLy8ti4cSMbN25kw4YN7Nq1i8rKSo4ePUpOTg4rVqzgD3/4A1qtdlpc1U+IdCVfm4UQfPLJJ2zevBmfz5d4rKOjg5kzZ1JXV8cvf/lLamtrsVgs+P1+Dhw4QF1dnYoRCzG9SctbCMHq1atpa2vjkUceITs7G0VRcDgc/OhHP8Ln8/H888/z+OOPA7B48WL2799PYWGhylELMX3JVDEhhBAiw0i3uRBCCJFhpHgLIYQQGUaKtxBCCJFhpHgLIYQQGUaKtxBCCJFhpHgLIYQQGUaKtxBCCJFhpHgLIYQQGeb/Ay4uCeDf9AaQAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "tags": [] } } ] }, { "cell_type": "markdown", "metadata": { "id": "ZjKuTkgzcfFt" }, "source": [ "Interpolate relative permeability data" ] }, { "cell_type": "code", "metadata": { "colab": { "base_uri": "https://localhost:8080/" }, "id": "tGoa5kUBcfFu", "outputId": "17c46b1a-ecb8-406b-99c5-610cc3d701c6" }, "source": [ "# Test intepolated relperm for Sw=0.575\n", "Sw_new = .575\n", "krw_new, kro_new = interpolate_relperm(Sw, krw, kro, Sw_new)\n", "\n", "print('At Sw={}, krw={} and kro={}'.format(Sw_new, krw_new, kro_new))" ], "execution_count": 4, "outputs": [ { "output_type": "stream", "text": [ "At Sw=0.575, krw=0.25827076696551754 and kro=0.06379291182871477\n" ], "name": "stdout" } ] }, { "cell_type": "markdown", "metadata": { "id": "0FkD9mmXGS1z" }, "source": [ "Initial condition" ] }, { "cell_type": "code", "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 382 }, "id": "HQhjO5Y-cfFv", "outputId": "1788dea6-8486-4f4e-a253-b1a105cc9f23" }, "source": [ "# set parameters for initial condition\n", "L = 4 \n", "nx = 41\n", "\n", "x = numpy.linspace(0.0, L, nx) \n", "Sw0 = numpy.full(nx, 0.2)\n", "Sw0[:15] = 1\n", "\n", "pyplot.plot(x, Sw0)\n", "pyplot.xlim(min(x), max(x))\n", "pyplot.xlabel('x'); pyplot.ylabel('Sw')\n", "pyplot.title('Initial Condition', size=20)\n", "\n", "pyplot.show()" ], "execution_count": 5, "outputs": [ { "output_type": "display_data", "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfYAAAFtCAYAAAD1Skg8AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3de3xU5Z3H8e/cuWSkDM6gAtoUKJFQVrlZDJcSAwWXfa3LWggqiRWhvgS1FlzZ2DZr26TqIirYXRXFKoKE2pStdXdjuy/tFkwJsEghrqWgZLkJM4CRQC4kmf0DcyCQTGYGjjln+Lz/ac6cMyfPj1P88jzPOedxRKPRqAAAQEpwdnYDAADAxUOwAwCQQgh2AABSCMEOAEAKIdgBAEghBDsAACmEYAdMkp2drezs7It2XFv27dunQYMGadGiRUl936pmzZqlQYMGGdsbN27UoEGDtGzZspjfi/c4IJW5O7sBgBXt27dPN910k8aMGaOXXnopqXMUFhae99m//du/6YorrtANN9wQ8zgzHTx4UK+++qrWr1+vAwcOqLGxUcFgUIMHD9Ytt9yS9D8yzDRw4EA988wzGjBggPHZiRMntGLFCuXn5+uyyy5r9zjgUkOwAyYZP378eZ89/fTTmjZtWqtgb+s4s7z11lv6x3/8RzkcDk2dOlV5eXlyOp36+OOPtW7dOpWVlWnq1Kn66U9/Kq/X+4W1qyOBQECTJ09u9dn27dv17LPP6u/+7u+MYG/rOOBSQ7ADX5AjR47owIEDnfb7N2/erIceekh9+/bVSy+9pH79+rXaP3/+fD3wwAP6zW9+o759++rBBx/spJbGZ/v27Z3dBMCSmGMHEjBo0CDdeeedOnTokB544AHdcMMNGjJkiKZNm6b169e3OvbsufNly5bpxhtvlCQ9++yzGjRokEpLS887rsWhQ4f0k5/8RBMmTNCQIUM0evRozZo1Sxs2bEi67cXFxYpGo1q6dOl5oS5JXbp00eLFizVu3Dj17du31b7du3fre9/7nrKyspSZmakbb7xR999/v/785z+3Om7RokUaNGiQDh48qBdeeEE5OTkaMmSIxo4dqyVLlqipqanV8R988IHy8/N13XXXacSIEZo7d64++uij89p27tx5dna2Fi9eLEm66aabjPn49ubYzWo/YEX02IEE1dbWKi8vTyNGjNDDDz+s/fv3a8WKFZo/f77efvtthUKh874zZcoUORwOLVu2TJMnT9aUKVP0ta99rc3znzx5UnfccYc++eQT3XXXXRowYICOHTumNWvW6K677tLPfvYz5eTkJNTmXbt2qbKyUllZWcrIyGj3OL/fr+XLl7f6bOfOnZo5c6bcbrdyc3OVnp6uffv2adWqVcrNzdXq1at17bXXtvrOU089pd27d+vOO++U2+3WmjVr9Pzzz6tXr17Kz8+XJB04cEB5eXlqamrSrFmzlJ6erj//+c+6++671b1795j1FBYW6sUXX1RFRYUKCwsVCATaPdas9gNWRbADCXr//ff10EMP6e677zY+czqdWrp0qf77v/9bt95663nfGTBggEaOHGn8HGseuKqqSv369dNtt92mb3/728bnY8aM0ZQpU/Taa68lHOzbtm2TJI0aNSqh70nSE088oZqaGq1Zs0bXX3+98fm4ceP0rW99S08++aRefPHFVt/ZuXOn1q5da8zTjx07VjfddJPefvttIxhfeeUVHT9+XEVFRa3+zK699lo9/PDDMds0fvx4/cd//IfRjnNHGL6I9gNWxVA8kCC32628vLxWn7X0vsPh8AWf/9prr9WKFSuMUK+trdVnn32myy+/XG63W/v370/4nEeOHJEk9e7dO6HvnTx5Uhs2bNCgQYNahaIkDR06VF/96ldVXl6u+vr6Vvtuv/32Vjff9enTR7169Wr151NeXi6n06mbb7651XenTp2qtLS0hNrZGe0HrIpgBxLUu3fv8+4Y9/l8kqTGxsaL8jvKy8uVl5en4cOH67rrrtPIkSM1cuRINTY2JjXP63A4JEnNzc0Jfa+qqkrNzc0aOHBgm/vT09PV2Nioffv2tfr86quvPu9Yn8/X6s9n7969CgaD6tatW6vj3G63rrnmmoTa2R4z2w9YFUPxQILMfgxs/fr1mjNnjvx+v7797W9r8ODBxpzz7Nmzkzpny7x/onflnzhxQpLOC98WLf+gOXnyZJufx1JXV9fm/QjS6Rv5LgYz2w9YFcEOWMzPf/5zNTc365lnntHo0aONz+vq6pK+K3vYsGGSpA0bNui+++6LeezRo0eNm9Fa/kFxbvC1qK2tbXVcIrp06XLeEHiL9n5fosxsP2BVDMUDFrNv3z45nU59/etfb/X5li1bEh5Kb9GvXz+NGDFCW7duVXl5ebvHNTQ0aObMmbr99tvV2NioL3/5y3K5XNq5c2ebx+/atUterzfmzWvtueqqqxQOh88L94aGBlVVVSV8vraY2X7Aqgh24Avicrkkqd1eaovLL79czc3NrYbNq6ur9cwzz6hr166qq6tL6vcvWrRIHo9HCxcuVGVl5Xn7T548qfvuu0979uzR2LFj5Xa71bVrV40fP147d+7U5s2bWx1fUVGhjz/+WBMmTEhqeqLlnoHf/va3rT5/88034+qxO52n//MV68/TzPYDVsVQPPAF6dOnjxwOh95880317NlTgwYN0pgxY847bsqUKdq0aZO++93vaubMmfrss8+0evVqzZgxQ06nU1u3bjVenpJIIH3ta1/T0qVLtWDBAk2fPl3f/OY3NWrUKPl8Pn300UcqLS1VJBLRPffco+985zvG9/7hH/5Bmzdv1rx58zRr1iz17dtXe/bs0erVq9WzZ0899NBDSf155Ofn61e/+pUKCwu1a9cuXX311frwww/129/+VkOGDNGOHTtifr+ll/3P//zPGjlypP72b/+2zePMaj9gVQQ78AW58sordc8992jlypV69tlndd9997UZ7Lm5uTp27JhKS0v16KOP6pprrtHcuXN16623auDAgSooKNC//uu/6qqrrtJ1112XUBuys7P19ttv69VXX9Xvf/97/f73v1d9fb0uv/xyjRs3Tvn5+ee9wCY9PV1r167V0qVLtXr1alVXV6tnz57Kzs7WvHnz2nyLXTzS09O1YsUKPfnkk3r55Zflcrl0/fXXa/ny5Vq6dGmHwZ6bm6v169dr/fr12rFjhyZNmtTu7zGj/YBVOaLRaLSzGwEAAC4O5tgBAEghBDsAACmEYAcAIIUQ7AAApBCCHQCAFGKbx90aG5t07NjFec2kFfXs2Y36bCyV60vl2iTqs7tUry8Y9Cf8Hdv02N1uV2c3wVTUZ2+pXF8q1yZRn92len3JsE2wAwCAjhHsAACkEIIdAIAUQrADAJBCCHYAAFIIwQ4AQAoh2AEASCEEOwAAKcTUYN+5c6dycnL02muvnbfvvffe06233qoZM2boZz/7mZnNAADgkmFasJ88eVI//vGPNXr06Db3/+QnP9GyZcv0+uuva8OGDdq1a5dZTQEA4JJh2rvivV6vli9fruXLl5+3b+/everRo4euvPJKSdL48eNVXl6uAQMGtHu+0nf+ouM19WY1t9N17+7TiRPUZxVdfW6NHXql3C5mqwDYi2nB7na75Xa3ffpwOKxAIGBsBwIB7d27N+b5Xv7NBxe1fUBHvtKvp0Zc2zvu45NZrMEuUrk2ifrsLtXrS5RtVnd7dM5oHfs0dVfw6dGjq6qrazu7GaaxU307Pj6i323epwOffKbw5d3i+k4w6Fc4fNzklnWOVK5Noj67uxTqS1SnBHsoFFIkEjG2Dx06pFAoFPM7wzJCKX/xqM8aTtadkiTVNzZ1cksAIHGdMoHYt29f1dTUaN++fWpsbNQ777yjrKyszmgKcB6f5/QykA0NBDsA+zGtx75jxw49/vjj2r9/v9xut8rKypSdna2+fftq4sSJ+qd/+ictWLBAknTzzTcrPT3drKYACfF+Huz1jc2d3BIASJxpwT5kyBCtXLmy3f0jR45USUmJWb8eSJrRYz9Fjx2A/fAsD3AOr+f0X4t6gh2ADRHswDnosQOwM4IdOIfXCHbm2AHYD8EOnMPHUDwAGyPYgXN4GYoHYGMEO3AOt8spl9OheobiAdgQwQ60wetx0mMHYEsEO9AGr8fFHDsAWyLYgTb43C418OY5ADZEsANt8Hpcqudd8QBsiGAH2uDzONXA6m4AbIhgB9rg9bjU2BRVUzPD8QDshWAH2uDj7XMAbIpgB9rAQjAA7IpgB9rA2+cA2BXBDrShZSiet88BsBuCHWhDy1A8PXYAdkOwA20402Mn2AHYC8EOtMHr5q54APZEsANt8HnpsQOwJ4IdaIPXzRw7AHsi2IE2MMcOwK4IdqANxnPsrPAGwGYIdqANvpY3z7HCGwCbIdiBNngZigdgU6YGe3FxsWbMmKHc3Fz96U9/arXvd7/7nf7+7/9eM2fO1GuvvWZmM4CE+RiKB2BTpgV7RUWFqqqqVFJSoqKiIhUVFRn7mpub9eMf/1jLly/XqlWr9M477+iTTz4xqylAwrwMxQOwKdOCvby8XDk5OZKk/v37q7q6WjU1NZKkY8eO6bLLLlMgEJDT6dTXv/51vffee2Y1BUjYmR47wQ7AXtxmnTgSiSgzM9PYDgQCCofDSktLUyAQ0IkTJ7Rnzx716dNHGzdu1KhRozo8ZzDoN6u5lkB91uFvaJQkRR2OuNttp/oSlcq1SdRnd6leX6JMC/ZzRaNR42eHw6HHHntMBQUF8vv96tu3b1znCIePm9W8ThcM+qnPQlr+/1pzoiGudtutvkSkcm0S9dndpVBfokwL9lAopEgkYmwfPnxYwWDQ2B41apRWr14tSXryySfVp08fs5oCJMzhcMjrcXJXPADbMW2OPSsrS2VlZZKkyspKhUIhpaWlGfvvvvtuHTlyRCdPntQ777yj0aNHm9UUICk+j4tXygKwHdN67MOGDVNmZqZyc3PlcDhUWFio0tJS+f1+TZw4UdOnT9ddd90lh8OhuXPnKhAImNUUICleN8EOwH5MnWNfuHBhq+2MjAzj50mTJmnSpElm/nrggvi8Ln12oqGzmwEACeHNc0A7vG4nPXYAtkOwA+3weVxqaGxW81lPdACA1RHsQDta3hd/6hSvlQVgHwQ70A5jhTfePgfARgh2oB3Gmuy8Lx6AjRDsQDta3hdfzwpvAGyEYAfa0bLCG3fGA7ATgh1oh7HCG8EOwEYIdqAdxlA8wQ7ARgh2oB3GzXM87gbARgh2oB0tc+z02AHYCcEOtIM5dgB2RLAD7fAac+wMxQOwD4IdaIfPzeNuAOyHYAfa4fVyVzwA+yHYgXb43NwVD8B+CHagHfTYAdgRwQ60w5hjZ3U3ADZCsAPtMO6KZ3U3ADZCsAPtMJ5jZ3U3ADZCsAPtcDodcruczLEDsBWCHYjB53HyHDsAWyHYgRi8Hhc9dgC2QrADMXg9Lp5jB2ArBDsQg8/DHDsAe3GbefLi4mJt27ZNDodDBQUFGjp0qLFv1apV+vWvfy2n06khQ4bokUceMbMpQFJahuKj0agcDkdnNwcAOmRaj72iokJVVVUqKSlRUVGRioqKjH01NTV66aWXtGrVKr3++uvavXu33n//fbOaAiTN53EpGpUam6Kd3RQAiItpwV5eXq6cnBxJUv/+/VVdXa2amhpJksfjkcfj0cmTJ9XY2Kja2lr16NHDrKYASfN+/vY5huMB2IVpQ/GRSESZmZnGdiAQUDgcVlpamnw+n+bNm6ecnBz5fD799V//tdLT0zs8ZzDoN6u5lkB91tPD30WS5L+sqy7/UteYx9qxvnilcm0S9dldqteXKFPn2M8WjZ4ZyqypqdHzzz+v//zP/1RaWpry8/P14YcfKiMjI+Y5wuHjZjez0wSDfuqzoOam03fEH/ikWtFTje0eZ9f64pHKtUnUZ3eXQn2JMm0oPhQKKRKJGNuHDx9WMBiUJO3evVv9+vVTIBCQ1+vViBEjtGPHDrOaAiTNeK0sj7wBsAnTgj0rK0tlZWWSpMrKSoVCIaWlpUmS+vTpo927d6uurk6StGPHDn35y182qylA0rwe5tgB2ItpQ/HDhg1TZmamcnNz5XA4VFhYqNLSUvn9fk2cOFGzZ89WXl6eXC6Xrr/+eo0YMcKspgBJO7MQDMEOwB5MnWNfuHBhq+2z59Bzc3OVm5tr5q8HLtiZpVsZigdgD7x5DojB9/lQPD12AHZBsAMxGD125tgB2ATBDsTAXfEA7IZgB2LgrngAdkOwAzGc6bET7ADsgWAHYvC6mWMHYC8EOxCDz8scOwB7IdiBGFpWd2MoHoBdEOxADC09dobiAdgFwQ7E0DLHTo8dgF0Q7EAMbpdDTodD9Y3MsQOwB4IdiMHhcMjrcaqhgR47AHsg2IEO+DwueuwAbINgBzrg9TiZYwdgGwQ70AGfx0WwA7ANgh3ogNfj4nE3ALZBsAMd8HlcamyKqqmZeXYA1kewAx048/Y5gh2A9RHsQAfOvC+e4XgA1kewAx1ghTcAdkKwAx04syY7Q/EArI9gBzrg9Z7+a0KPHYAdEOxAB3wsBAPARgh2oANeT8scO0PxAKyPYAc64PN8/rhbIz12ANbnNvPkxcXF2rZtmxwOhwoKCjR06FBJ0qFDh7Rw4ULjuL1792rBggX6m7/5GzObAyTF6LGzwhsAGzAt2CsqKlRVVaWSkhLt3r1bBQUFKikpkST17t1bK1eulCQ1NjZq1qxZys7ONqspwAXxeXjcDYB9mDYUX15erpycHElS//79VV1drZqamvOO+9WvfqVvfvOb6t69u1lNAS5IS4+9gaVbAdiAaT32SCSizMxMYzsQCCgcDistLa3Vcb/4xS+0YsWKuM4ZDPovahuthvqsKfRZvSTJ5XHFrMGu9cUjlWuTqM/uUr2+RJk6x362aDR63mdbt27VV77ylfPCvj3h8PGL3SzLCAb91GdRtSdOB/un1bXt1mDn+jqSyrVJ1Gd3l0J9iTJtKD4UCikSiRjbhw8fVjAYbHXMu+++q9GjR5vVBOCi4HE3AHZiWrBnZWWprKxMklRZWalQKHRez3z79u3KyMgwqwnARWE87sbNcwBswLSh+GHDhikzM1O5ublyOBwqLCxUaWmp/H6/Jk6cKEkKh8Pq1auXWU0ALgovd8UDsBFT59jPflZd0nm98zfffNPMXw9cFCwCA8BOePMc0AGPm0VgANgHwQ50wOlwyOtxMscOwBYIdiAOXreLHjsAWyDYgTj4PC7m2AHYAsEOxMHrcdJjB2ALBDsQB5/HxbKtAGyBYAfi4P18KL65jVcjA4CVEOxAHFqeZT/FCm8ALI5gB+Lg9fAsOwB7INiBOJx5+xzBDsDaCHYgDqzwBsAuCHYgDqzwBsAuCHYgDl43Q/EA7IFgB+Lg8zIUD8AeCHYgDl43Q/EA7IFgB+LgM26eI9gBWBvBDsTBy+NuAGyCYAfi4ONxNwA2QbADcfDyuBsAmyDYgTgYPXZWeANgcQQ7EAdjjr2BoXgA1kawA3FoefMcPXYAVhdXsFdWVprdDsDSuCsegF244znoqaee0p49e/RXf/VXGjt2rLKyshQMBs1uG2AZZ1Z3YygegLXFFewvvviiGhsbtX37dm3cuFEPP/ywjh49qnXr1pndPsASWI8dgF3EFeynTp3Stm3bVFFRoS1btqi+vl433HBDh98rLi7Wtm3b5HA4VFBQoKFDhxr7Dh48qO9973s6deqUBg8erB/96EfJVwGYzOV0yu1yMBQPwPLiCvaRI0dq5MiRuu2225Sfn6/u3bt3+J2KigpVVVWppKREu3fvVkFBgUpKSoz9jz32mO666y5NnDhRjz76qA4cOKCrrroq+UoAk3ndLnrsACwvrpvnnnrqKQ0YMEAvvviiHnzwQb3wwgt6//33Y36nvLxcOTk5kqT+/fururpaNTU1kqTm5mZt2bJF2dnZkqTCwkJCHZbn8xLsAKwvrh77hAkTNGHCBEnS+++/r+eee05PP/20Pvjgg3a/E4lElJmZaWwHAgGFw2GlpaXp6NGj6t69u37605+qsrJSI0aM0IIFCzpsRzDoj6e5tkV91tbV51ZtfWO7ddi9vlhSuTaJ+uwu1etLVFzB/tprr2nTpk368MMPNWDAAI0dO1YPP/xwQr8oGo22+vnQoUPKy8tTnz59NHfuXL377rv6xje+EfMc4fDxhH6nnQSDfuqzOLfTodr6xjbrSIX62pPKtUnUZ3eXQn2JijkUX1NTo5///Oc6duyY8vLyNGvWLO3du1cbN26U3x/7l4VCIUUiEWP78OHDxiNyPXv21FVXXaWrr75aLpdLo0eP1l/+8peEGw98kbwelxpONbf6RyoAWE3MYP/hD3+oI0eO6L777lMgENDSpUu1aNEiZWVlqaioKOaJs7KyVFZWJun0C25CoZDS0tIkSW63W/369dOePXuM/enp6RehHMA8Po9TzdGoGpsIdgDWFXMofu/evVqyZIkkqaysTJMnT9aNN94oSfrNb34T88TDhg1TZmamcnNz5XA4VFhYqNLSUvn9fk2cOFEFBQVatGiRotGovvrVrxo30gFWZbx9rrFJHjdvYwZgTTGDvVu3bsbPFRUVuvXWW41th8PR4ckXLlzYajsjI8P4+ZprrtHrr78ed0OBzmas8NbQpO5dPJ3cGgBoW8xuR1NTk44cOaL/+7//09atW5WVlSVJOnHihGpra7+QBgJWcabHzmtlAVhXzB77nDlzdPPNN6uurk7z589Xjx49VFdXp9tuu03Tp0//otoIWILxWtkGnmUHYF0xg338+PFav3696uvrjRvfunTpooceekhjxoz5QhoIWIXvrDl2ALCqDp9j93g88nhazycS6rgUtQzF8/Y5AFbGrb1AnFi6FYAdEOxAnFi6FYAdEOxAnM702Al2ANZFsANx8rpb5tgZigdgXQQ7ECef9/RfF3rsAKyMYAfidKbHTrADsC6CHYgTd8UDsAOCHYgTd8UDsAOCHYgTb54DYAcEOxAn71mruwGAVRHsQJx8rO4GwAYIdiBObpdDDgdz7ACsjWAH4uRwOOTzuHiOHYClEexAArweF2+eA2BpBDuQAJ/HSY8dgKUR7EACvAzFA7A4gh1IgI+heAAWR7ADCfC6nWpsalZzc7SzmwIAbSLYgQS0PMvOI28ArIpgBxLgNRaCIdgBWBPBDiTA6LHz9jkAFuU28+TFxcXatm2bHA6HCgoKNHToUGNfdna2rrjiCrlcp/9DuXjxYvXu3dvM5gAXrGWFtwbeFw/AokwL9oqKClVVVamkpES7d+9WQUGBSkpKWh2zfPlyde/e3awmABfdmR47wQ7Amkwbii8vL1dOTo4kqX///qqurlZNTY1Zvw74Qhhz7PTYAViUaT32SCSizMxMYzsQCCgcDistLc34rLCwUPv379fw4cO1YMECORyOmOcMBv1mNdcSqM/6Al/qJknq0s13Xj2pUF97Urk2ifrsLtXrS5Spc+xni0ZbP/d7//33a+zYserRo4fmzZunsrIyTZ48OeY5wuHjZjaxUwWDfuqzgVMNpyRJ4SM1repJlfraksq1SdRnd5dCfYkybSg+FAopEokY24cPH1YwGDS2b7nlFvXq1Utut1vjxo3Tzp07zWoKcNEYc+wMxQOwKNOCPSsrS2VlZZKkyspKhUIhYxj++PHjmj17thoaGiRJmzZt0sCBA81qCnDRGHPsPO4GwKJMG4ofNmyYMjMzlZubK4fDocLCQpWWlsrv92vixIkaN26cZsyYIZ/Pp8GDB3c4DA9Yge/zx9148xwAqzJ1jn3hwoWttjMyMoyf8/PzlZ+fb+avBy46r5s3zwGwNt48ByTA5+Vd8QCsjWAHEuB1f/7mOZZuBWBRBDuQAFZ3A2B1BDuQAFZ3A2B1BDuQgDM9dobiAVgTwQ4kwNOyuhs9dgAWRbADCXA6HPK6ncyxA7Asgh1IkNfj4s1zACyLYAcS5PM4eVc8AMsi2IEEne6xE+wArIlgBxLk9biYYwdgWQQ7kCCfx6WGU82KRqOd3RQAOA/BDiTI2/LIGzfQAbAggh1IkI+3zwGwMIIdSFDL0q3MswOwIoIdSFDL0q2s8AbAigh2IEEtS7fSYwdgRQQ7kCDm2AFYGcEOJKjlrnhWeANgRQQ7kCB67ACsjGAHEuT1cFc8AOsi2IEEGT12XlADwIIIdiBBxhw7K7wBsCCCHUjQmR47wQ7Aegh2IEHMsQOwMlODvbi4WDNmzFBubq7+9Kc/tXnMk08+qVmzZpnZDOCiOnNXPHPsAKzHtGCvqKhQVVWVSkpKVFRUpKKiovOO2bVrlzZt2mRWEwBTnHmOnR47AOsxLdjLy8uVk5MjSerfv7+qq6tVU1PT6pjHHntMDz74oFlNAEzBc+wArMxt1okjkYgyMzON7UAgoHA4rLS0NElSaWmpRo0apT59+sR9zmDQf9HbaSXUZw/d0rqc/sHhbFVTqtTXllSuTaI+u0v1+hJlWrCfKxqNGj9/+umnKi0t1csvv6xDhw7FfY5w+LgZTbOEYNBPfTbR2HR6bv34iXqjplSq71ypXJtEfXZ3KdSXKNOG4kOhkCKRiLF9+PBhBYNBSdIf//hHHT16VLfffrvmz5+vyspKFRcXm9UU4KJyu5xyOR0MxQOwJNOCPSsrS2VlZZKkyspKhUIhYxh+8uTJ+vd//3etXbtWzz77rDIzM1VQUGBWU4CLzudxcfMcAEsybSh+2LBhyszMVG5urhwOhwoLC1VaWiq/36+JEyea9WuBL4TX4+RxNwCWZOoc+8KFC1ttZ2RknHdM3759tXLlSjObAVx0Po9LdbxSFoAF8eY5IAlej4tXygKwJIIdSILP41J9Q3Orpz0AwAoIdiAJXo9TzdGompoJdgDWQrADSfCxEAwAiyLYgSR4WQgGgEUR7EASfCwEA8CiCHYgCV43C8EAsCaCHUiCz8scOwBrItiBJHjdp//qMMcOwGoIdiAJ3BUPwKoIdiAJZ+6KJ9gBWAvBDiSBHjsAqyLYgSR4PcyxA7Amgh1IAj12AFZFsANJMObYWeENgMUQ7EASjB57A0PxAKyFYAeSYMyx05b5BbUAAAvMSURBVGMHYDEEO5AE5tgBWBXBDiSB1d0AWBXBDiSB1d0AWBXBDiTB7XLKId48B8B6CHYgCQ6HQ16vix47AMsh2IEk+dxO5tgBWA7BDiTJ66HHDsB6CHYgST6vizl2AJbjNvPkxcXF2rZtmxwOhwoKCjR06FBj39q1a/XGG2/I6XQqIyNDhYWFcjgcZjYHuKi8bpfqGYoHYDGm9dgrKipUVVWlkpISFRUVqaioyNhXW1urt956S6tWrdKaNWv00UcfaevWrWY1BTCFz+NUY1Ozmpujnd0UADCYFuzl5eXKycmRJPXv31/V1dWqqamRJHXt2lWvvPKKPB6PamtrVVNTo2AwaFZTAFN4efscAAsybSg+EokoMzPT2A4EAgqHw0pLSzM+e+GFF/Tqq68qLy9P/fr16/CcwaDflLZaBfXZy2VpPkmS/7KuklKvvrOlcm0S9dldqteXKFPn2M8WjZ4/XDl37lzl5eVpzpw5Gj58uIYPHx7zHOHwcbOa1+mCQT/12Uy0+fT8+oFDn6nnZV1Srr4WqXjtzkZ99nYp1Jco04biQ6GQIpGIsX348GFjuP3TTz/Vpk2bJEldunTRuHHj9D//8z9mNQUwRctCMA0NDMUDsA7Tgj0rK0tlZWWSpMrKSoVCIWMYvrGxUYsWLdKJEyckSdu3b1d6erpZTQFMYcyxs3QrAAsxbSh+2LBhyszMVG5urhwOhwoLC1VaWiq/36+JEydq3rx5ysvLk9vt1qBBg3TTTTeZ1RTAFPTYAViRqXPsCxcubLWdkZFh/Dxt2jRNmzbNzF8PmMrbssJbI8+yA7AO3jwHJMnosfO4GwALIdiBJHndPMcOwHoIdiBJPm9Lj52heADWQbADSfK6T//1YSgegJUQ7ECSfLxSFoAFEexAkrwehuIBWA/BDiTJ1/K4Gz12ABZCsANJ8vK4GwALItiBJDHHDsCKCHYgSS1vnmvgzXMALIRgB5JkLALDu+IBWAjBDiTJ6XDI43aqgdXdAFgIwQ5cAJ/HpXoedwNgIQQ7cAG8Hid3xQOwFIIduACne+wEOwDrINiBC+B1u3jzHABLIdiBC+D7fCg+Go12dlMAQBLBDlwQr8elqHiWHYB1EOzABWh5+1xdfWMntwQATiPYgQvgZSEYABZDsAMXwMfb5wBYDMEOXABeKwvAagh24AK0BHtdA3PsAKyBYAcugI85dgAWQ7ADF+BMj51gB2ANbjNPXlxcrG3btsnhcKigoEBDhw419v3xj3/UkiVL5HQ6lZ6erqKiIjmd/DsD9sLNcwCsxrQkraioUFVVlUpKSlRUVKSioqJW+3/4wx9q6dKlWrNmjU6cOKE//OEPZjUFMI3xuBtz7AAswrQee3l5uXJyciRJ/fv3V3V1tWpqapSWliZJKi0tNX4OBAI6duyYWU0BTONzn+6x/0f5HlVUftK5jTGJz+dWfQq/gIf67C2V6/O4nCq464aEv2dasEciEWVmZhrbgUBA4XDYCPOW/z18+LA2bNigBx54oMNzBoN+cxprEdRnP4Ojktvl1McHPtPHBz7r7OYASCEOR3LfM3WO/WxtLZJx5MgR3XPPPSosLFTPnj07PEc4fNyMpllCMOinPhvyOaSlD4yR/7JuOnKkprObY4pevdJStjaJ+uwuletzu5JLdtOCPRQKKRKJGNuHDx9WMBg0tmtqajRnzhx997vf1ZgxY8xqBmC6Ll63vuT36VRdQ2c3xRSpXJtEfXaX6vUlw7Sb57KyslRWViZJqqysVCgUMobfJemxxx5Tfn6+xo0bZ1YTAAC45JjWYx82bJgyMzOVm5srh8OhwsJClZaWyu/3a8yYMVq3bp2qqqr0xhtvSJKmTp2qGTNmmNUcAAAuCabOsS9cuLDVdkZGhvHzjh07zPzVAABckngjDAAAKYRgBwAghRDsAACkEIIdAIAUQrADAJBCCHYAAFIIwQ4AQAoh2AEASCGOaFurswAAAFuixw4AQAoh2AEASCEEOwAAKYRgBwAghRDsAACkEIIdAIAUYup67MkqLi7Wtm3b5HA4VFBQoKFDhxr73nvvPS1ZskQul0vjxo3TvHnzOrGliYtVW3Z2tq644gq5XC5J0uLFi9W7d+/OamrSdu7cqXvvvVd33nmn7rjjjlb77H79YtWWCtfviSee0JYtW9TY2KjvfOc7mjRpkrHP7tdOil2fna9fbW2tFi1apCNHjqi+vl733nuvJkyYYOy3+7XrqD47X7uz1dXVaerUqbr33ns1bdo04/OEr1/UYjZu3BidO3duNBqNRnft2hWdPn16q/1TpkyJHjhwINrU1BSdOXNm9C9/+UtnNDMpHdU2YcKEaE1NTWc07aI5ceJE9I477oh+//vfj65cufK8/Xa+fh3VZvfrV15eHr377ruj0Wg0evTo0ej48eNb7bfztYtGO67Pztfvrbfeir7wwgvRaDQa3bdvX3TSpEmt9tv92nVUn52v3dmWLFkSnTZtWvSXv/xlq88TvX6WG4ovLy9XTk6OJKl///6qrq5WTU2NJGnv3r3q0aOHrrzySjmdTo0fP17l5eWd2dyExKotVXi9Xi1fvlyhUOi8fXa/frFqSwUjR47UM888I0m67LLLVFtbq6amJkn2v3ZS7Prs7uabb9acOXMkSQcPHmzVW02FaxervlSxe/du7dq1S9/4xjdafZ7M9bPcUHwkElFmZqaxHQgEFA6HlZaWpnA4rEAg0Grf3r17O6OZSYlVW4vCwkLt379fw4cP14IFC+RwODqjqUlzu91yu9v+v5Xdr1+s2lrY+fq5XC5169ZNkvTGG29o3LhxxtCm3a+dFLu+Fna+fpKUm5urTz75RM8995zxWSpcuxZt1dfC7tfu8ccf1w9+8AOtW7eu1efJXD/LBfu5oin8xttza7v//vs1duxY9ejRQ/PmzVNZWZkmT57cSa1DolLl+v3ud7/TG2+8oRUrVnR2U0zRXn2pcP3WrFmj//3f/9VDDz2kX//617YLt460V5/dr926det03XXXqV+/fhflfJYbig+FQopEIsb24cOHFQwG29x36NAhWw2LxqpNkm655Rb16tVLbrdb48aN086dOzujmaax+/XrSCpcvz/84Q967rnntHz5cvn9fuPzVLl27dUn2fv67dixQwcPHpQkXXvttWpqatLRo0clpca1i1WfZO9rJ0nvvvuu/uu//kvTp0/XL37xC/3Lv/yL3nvvPUnJXT/LBXtWVpbKysokSZWVlQqFQsZQdd++fVVTU6N9+/apsbFR77zzjrKysjqzuQmJVdvx48c1e/ZsNTQ0SJI2bdqkgQMHdlpbzWD36xdLKly/48eP64knntDzzz+vL33pS632pcK1i1Wf3a/f5s2bjRGISCSikydPqmfPnpJS49rFqs/u106Snn76af3yl7/U2rVr9a1vfUv33nuvbrzxRknJXT9Lru62ePFibd68WQ6HQ4WFhfrggw/k9/s1ceJEbdq0SYsXL5YkTZo0SbNnz+7k1iYmVm2vvPKK1q1bJ5/Pp8GDB+sHP/iB7YbSduzYoccff1z79++X2+1W7969lZ2drb59+9r++nVUm92vX0lJiZYtW6b09HTjsxtuuEGDBg2y/bWTOq7Pztevrq5OjzzyiA4ePKi6ujrNnz9fn376acr8d7Oj+ux87c61bNky9enTR5KSvn6WDHYAAJAcyw3FAwCA5BHsAACkEIIdAIAUQrADAJBCCHYAAFIIwQ4AQAoh2AEASCEEO4A2vfzyy/r+978vSfroo480efLklFuNEEhFBDuANuXn5+vjjz/Wli1b9Oijj+pHP/pRq5UIAVgTb54D0K6qqirdcccdmjx5sh555JHObg6AONBjB9Cu6upqdevWzVhZC4D1EewA2lRfX6/CwkI999xz8ng8WrduXWc3CUAcGIoH0KYnnnhC3bt317x58xSJRDRjxgytWrVKV1xxRWc3DUAMBDsAACmEoXgAAFIIwQ4AQAoh2AEASCEEOwAAKYRgBwAghRDsAACkEIIdAIAUQrADAJBC/h8ELBHZbaNXXgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "tags": [] } } ] }, { "cell_type": "markdown", "metadata": { "id": "FCZb5n30Ie2B" }, "source": [ "Run simulation" ] }, { "cell_type": "code", "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 555 }, "id": "h33gOIgUcfFw", "outputId": "775b750a-81e6-4483-b338-16c8482bcadd" }, "source": [ "# Set parameters for simulation\n", "nt = 70\n", "L = 4\n", "sigma = 0.1\n", "bc_value = Sw0[0]\n", "u_max = 1\n", "\n", "muw = 0.5E-3\n", "muo = 1E-3\n", "q = 200 # m3/hr\n", "A = 30 # m2\n", "poro = 0.24\n", "\n", "# Simulation\n", "nt = [10, 50, 70, 90]\n", "\n", "pyplot.figure(figsize=(16,9))\n", "for i in range(4):\n", " pyplot.subplot(2,2,i+1)\n", " buckley_leverett1d(nt[i], Sw0, L, nx, sigma, bc_value, muw, muo, q, A, poro, Sw, krw, kro)" ], "execution_count": 7, "outputs": [ { "output_type": "display_data", "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA7UAAAIaCAYAAAD/ZGq6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdeXxU9bk/8M+ZfZKZJDPJZE8gCZBAkE1AIQiCgNbtWn9WqLZ6r1u9andbvdQr3mqotdVrq71ttba2tbX0Z7lq3ejPurGEnRAIS0gCIRvJTPaZzGTW3x9JBiLZmTNn5szn/Xr1ZSZzcs7zdMI8eeb7Pd+vEAgEAiAiIiIiIiKKQgqpAyAiIiIiIiKaLDa1REREREREFLXY1BIREREREVHUYlNLREREREREUYtNLREREREREUUtNrVEREREREQUtURtaquqqrB69Wq89tprFzy3c+dO3HLLLVi3bh1+8YtfiBkGERERDWBtJiIiuRGtqe3t7cWTTz6JJUuWDPv8U089hRdeeAGvv/46duzYgerqarFCISIiIrA2ExGRPInW1Go0Grz88stITU294Ln6+nokJiYiIyMDCoUCK1asQFlZmVihEBEREVibiYhInkRralUqFXQ63bDPWa1WmM3m4GOz2Qyr1SpWKERERATWZiIikieV1AGM13++/Rs43Z5hnzMbdchMNQAAnC4vaho7RzxPQVYS9Lr+tJta7WjvcQ17nF6rQkF2UvDxkRrbiOfMshhgSuj/I6Gj24VGq33EY2cXpAS/rmnohLPPO+xxzOnic2pwnoJf04Prcm7GnUvXjHgtIiKanMfefhku9wjv+Ql6ZFriAQC9Li9qR3nPn5adBJ22/z2/sdWOjhHe8+N0auRnJQIAAgGgsnZ8dayty4Vm2/B1TICA4oLk4OPq+k7mJGJOjT1N8MXZMNe8AD9Yc++I1yIimghJmtrU1FTYbOfe4FpaWoadCnW+J2+8B1Zrz/guMG2cgWSP8zgA103g2LFYLMb+XCZyzgjNKZjLJK4vdk4bP3gVNhxFY0fLmL87F+QRxZhLZJJLLnLJA+jPhc6ZTG1+6sZ7Ja3N1+eE7pyRUptDkVOk1ebzc3r2zW2oxd9xrOMomlo6oFaM/Keo3N5vmEvkkUsucskDmHxtlmRLn+zsbNjtdjQ0NMDr9eLjjz9GSUmJFKGQjKXF939K7RYcEkdCRBT5WJspHApTc6DxJsEdcOFY2wmpwyEimRBtpPbIkSP48Y9/jMbGRqhUKmzduhWrVq1CdnY21qxZgyeeeALf/e53AQDXXnst8vLyxAqFYtQVMwtQWbENAXWv1KEQEUUE1maS2vVLp0Jzegneqn0f+1rKMcdSLHVIRCQDojW1s2fPxh//+McRn1+0aBE2b94s1uWJYNaZAAAdrpHvDyIiiiWszRQJLk2bi7dq30eF7Shc3j7oVFqpQyKiKCfJ9GOicDCqEwAAba4OBAIBiaMhIiIir88Ptd+AKYZcePweHLYdlTokIpIBNrUkW4eruhHwquDxe+DwcAoyERGR1HYeOYtvvbAd/o5MAMC+loMSR0REcsCmlmRLr1Uh4O7f+qDd1SFxNERERKQf2BJI78yGQlDgaHsV7B4u6EhEF4dNLcmWXqNEoE8PgE0tERFRJNBrlAAAr0uNQtM0+AN+HGw9LHFURBTt2NSSbOl1KgTcA01tHxeLIpJCc3MT7r77q8M+19ZmwzPPlAIAyssPoKOjHQDw6KPfGfF8Bw7sw2OPfT/0gRJRWAyO1Pb2+bAwbR4AYH9LuZQhEcUcOdZmNrUkW3qNCoE+Tj8milTJySn4/vd/AAB49923g4Xz6aefkzIsIhKRbqCpdfZ5MdcyGyqFCtWdp7hTAVGEiNbaLNqWPkRS67+ndnD6MYslkZRKS59ASooFJ04cQ0vLWTz++FNISEjAY489gvvvfxDbtn2CU6dq8dRTz+Duu2/Hu+/+E3v37sZvfvMrqNVqGI1G/PCHT0udBhFdpLjBptbthV6lw+zkmSi3Hsb+1kNYnbtC4uiIYoucajObWpItvVYF/+A9tU6O1BLd9fRHIz53xzWFuHJeFgDgk/JG/OGDEyMe+9tHV03q+m63G8899yLefPMNfPDBu7j11i8DABYtuhzTps3Ad77zfaSnpweP7+npwcaNTyEzMwtPPvk4du8uQ1xc3KSuTUSRQa/tv6fW2ecFACxKm4dy62HsaylnU0sxibU5NNjUkmypVQrcuWoeNrfsQnsfm1oiqc2dOx8AYLGk4ejRyjGPT0pKwo9//BR8Ph+amhpx6aWLIqJwEtHkadVK/PtNs4PNbXFyEXRKHep7GtHiaEVafKrEERLFFrnUZja1JGvLZk3F31qVcHh60edzQ6vUSB0SkWTG+ynulfOygp8Mf57FYoTV2jOp6yuVyuDXgUBgzON/9KMn8ZOfPI+pU/Pw3HM/ntQ1iSiyCIKARUXnGle1Uo15ltnYdXYf9rWU47r8tRJGRxR+rM2hwYWiSNYUggJJuiQAQAcXiyKKWAqFAj6fb8j3HA470tLS0dPTgwMH9sPj8UgUHRGJaXAV5H2t5eP6o5qIwiOaajObWpK1HYeb4XNqAQBtXCyKKGLNm7cAjz32CGpra4Lfu/nmL+Hf//1uPPNMKW6//Q689tqraGuzSRglEYXCtoomvPFJDdq7XQCAGaYCGNUGtPbaUN/TKHF0RDQommqzEIiij8QmO6weaS5mikCkifRcnvnzAdQot0FlacT6wptxRdblwx4X6XlMBHOJTHLJRS55AP250MWT0+8DcwmPH722HycbuvDIbfNRmGsCAPy16k182rATV+Usx83TrwcQ+XlMBHOJTHLJRS55AJOvzRypJVnTa1UIDK6AzOnHREREktMH96o9N61xYVr/YjX7Ww/BH/BLEhcRRS82tSRrOo0KAbcOALixOxERUQQ419R6g9/LS8hFss6Ezr4u1HSekio0IopSbGpJ1uI4UktERBRR9JqBvWrd55paQRBw6cCCUXtbyiWJi4iiF5takjWdVomAe7Cp5UgtERGR1IYbqQXOrYJc3noYXr/3gp8jIhoJm1qStTjtuenHnX1d8Pl9Y/wEERERiUk3zD21AJBlyEBmfDoc3l4ca6+SIjQiilJsaknWjHEamAx6aBGHAALo7OuSOiQiIqKYlhCnhsmohUZ14Z+hg1OQ93EKMhFNAJtakrVlczLw7IMlyEywAOAUZKJIUF5+AB0d7ZP62Y8++hBr1lyB2trq4Pf27t2Ne++9A1/72r/h1Vd/E6owiUgkK+Zl4dkHS3DjsrwLnhucglxhrYTL2xfu0Ihi1mRr8yuv/Brr138RDz10Hx566D68886bAMJfm9nUUkww65IAcLEookjw7rtvT6pwHjy4H7t27UBBwfQh3//Zz36Kp556Br/85SvYs2cXTp2qDVWoRBRmKXoz8hJy4fZ7sL+pQupwiGLGZGszAHzpS+vx4osv4cUXX8L1198EIPy1WSXq2YkihFnXv7k7R2qJwue99/6OiopydHZ24MyZOtx221eRlpaObds+walTtXjqqWeQnp4OANi5czv+/Oc/DPn5G2+8GWvXXhN8XFhYhPnzL8VDD90X/F5jYwOMxgSkpfWfZ8mSEuzfvwd5eflhyJCIxLAwbT5OdZ/B9rq9mFFUJHU4RLIS6to8HClqM5takrWz7b14+k8HoM3oAswcqaXY9T+HfovKtuMhPWdxchEemHvXqMfU1FTjV7/6LRoa6rFx4wa8+uqfMW3aDHznO98PFk0AWLp0GZYuXTbqueLi4i/4Xnt7G5KSTMHHJpMJjY2NE8yEiMKp0ebAT14/CEuiDj+4Y+EFzy9Im4M3Tr6N8rNH4SjoRbw6ToIoicQnh9oMAB9//E9s2/YpNBoNvvWt70lSmzn9mGRNrVSg2+FGn0MDgE0tUbjNnj0HSqUSFksqHA676NcLBES/BBFdJLVSQLfDjS6He9jnEzRGFJqmwef3ocJ2NMzREclfKGvzkiUluOee+/H88/+DtWu/gOef/8kFx4SjNnOklmRtcC88t0MLAUBHH6cfU2wa61Pb8bJYjLBae8Z9vFKpDH4dGKWqTXaKU0qKBe3tbcHHVmsrUlJSxh0fEYXfSPvUnm92ykwc7ziJ4+1VWJJx4WgukRzIoTbPmjU7+PWyZSvwy1++IEltZlNLsqbT9P+jddk10KP/ntpAIABBEKQNjCiGKRQK+HxD96cc7xSnz8vIyITD4UBzcxMsllTs3Lkdjz/+ZKhCJSIR6M/bp3akmjzT3L8g3In2atZtojCYbG1+/vmfYuXKqzB37nwcPLgP+fkFktRmNrUkawqFAJ1GCZcb0Ct1cPpcsHscMGoMUodGFLPmzVuAxx57BD/60bPIzy8Y98+9886b+OCD91BdXYVNm36IKVOm4j//84d4+OFH8cQTPwAArFq1Brm5U8QKnYhCQKVUQK1SwOP1w+31Q6tWXnBMWlwqTPpEdDi70OQ4iyxDhgSREsWOydbmG264CT/5ySaoVCoIgoBHHnkMAMJem4XAaGPOEWYiw+qRbKJTBCJZNOTy3V/sQEdPH/JWlOOs8yy+v/DrmJKQM+SYaMhjvJhLZJJLLnLJA+jPhS6enH4fmEv4fOuF7eh2uPHcQyVIMmiHPWZz7d/w2enduHna9bgqd3mYIwytaHhNxou5RB655AFMvjZzoSiSvcFpTkZVAgCgjYtFERERSUo/cHvQaPfVzkmbCQA43n4yLDERUfTi9GOSvTULs+Hs88EW14mTPUAH96olIiKS1NrFuehz+xCvV494zCVp/XvUVnfWwuP3Qq3gn61ENDy+O5DsrZiXBQD48MwpoIXb+hAREUlt5fysMY8x6RORGZ+OJsdZnOqqwwzT+O/zI6LYwunHFDPMuv5NoNs5UktERBQVCs3TAAAnOAWZiEbBppZkr6HVjl1Hz8Ln0gHgSC0REZHU6gdqc3ObY9Tjikz9W/sc62BTS0QjY1NLsrfjSDNeevsoGpv6995iU0tERCStbRVNeOntozhc0zbqcdOS8qEUlDjT3YBeT2+YoiOiaMOmlmRPr+m/ddzXp4ZaoUKv1wmX1yVxVERERLFrsDb3jrL6MQDoVFrkJeYigACqOmvDERoRRSE2tSR7g1v6uNx+mHRJAHhfLRERkZTO1WbfmMcOTkHm1j5ENBJRVz/etGkTDh06BEEQsGHDBsyZMyf43J/+9Ce8/fbbUCgUmD17Nn7wgx+IGQrFMJ323F545hQTWnttaHd1INOQLnFkREThx9pMkUA/UJvHGqkFgCLzdLxz6h9cLIqIRiTaSO2ePXtQV1eHzZs3o7S0FKWlpcHn7HY7XnnlFfzpT3/C66+/jpqaGpSXl4sVCsW4uIFPg519Xpg5UktEMYy1mSJFcKR2HE1trjEbepUOrU4b2pxcF4OILiRaU1tWVobVq1cDAAoKCtDV1QW73Q4AUKvVUKvV6O3thdfrhdPpRGJiolihUIzTD2lqB7f1YVEkotjD2kyR4vzaPBalQokZSf171J7gKshENAzRph/bbDYUFxcHH5vNZlitVhgMBmi1Wjz44INYvXo1tFotrrvuOuTl5Y15TovFKFa4YcdcwifD2V8wPf4AplgygFNAL+wXxB3peUwEc4lMcslFLnnEItbm0TGX8MmwewD01+bRYh187tLc2Thkq8Qpxyn8i+WqsMQYapH+mkwEc4k8csljskS9p/Z8gUAg+LXdbsevf/1rfPDBBzAYDLjzzjtx/PhxFBUVjXoOq7VH7DDDwmIxMpcwMmoUeOFbV0CnUaK26xQAoLnLOiTuaMhjvJhLZJJLLnLJA+AfAABr8/nk9rsd6bkk6ZXB2jxSrOfnka3JBQBUnD2OltYuKIToWus0Gl6T8WIukUcueQCTr82ivSOkpqbCZrMFH7e2tsJisQAAampqkJOTA7PZDI1Gg4ULF+LIkSNihUIxTqVUIF6nhlKhOG/6Me+pJaLYw9pMkeL82jweqfoUmLRJsHscaLQ3ixwdEUUb0ZrakpISbN26FQBQWVmJ1NRUGAwGAEBWVhZqamrgcvXvFXrkyBFMnTpVrFCIgpK0iRAgoKuvG17/2PfxEBHJCWszRStBEFBk5tY+RDQ80aYfL1iwAMXFxVi/fj0EQcDGjRuxZcsWGI1GrFmzBnfffTfuuOMOKJVKzJ8/HwsXLhQrFCI8/38PoafXg0dvX4BEbQI6+7rQ2deFFH2y1KEREYUNazNFikAggOf+egi9Li82fHXBuEZsi0zTUNa8F8fbT2LNlCvFD5KIooao99Q+/PDDQx6ff1/O+vXrsX79ejEvTxRU29QNu9MTXAG5s68L7a5ONrVEFHNYmykSCIKA2qZuOPu8cLl9iNeN3dQWDozU1nSdgsfngVqpFjtMIooS0XWXPdEkDW7yPnSvWm7rQ0REJJVgbXaN73Ygo8aALEMGPH4varvqxAyNiKIMm1qKCcH98Nzcq5aIiCgSnKvNvnH/TJFp4L5a7ldLROdhU0sxIW6wcLrOH6nlCshERERSCTa1feNfuJGLRRHRcNjUUkzQac59Gjw4UtvBppaIiEgyes3Em9ppSXlQCUrU9zTC7nGIFRoRRRk2tRQTzv80mNOPiYiIpHf+ehfjpVFqkJ84FQEEUNVRI1ZoRBRlRF39mChSzMhJhCAAKYk6mHXxAID2vk74A34oBH62Q0REFG6FOUlQKRUwJ+gm9nPm6ajqrMHx9pNYkDpHpOiIKJqwqaWYsGJeFlbMywo+jlfHweHpRY/bgUStUcLIiIiIYtPKBdlYuWDiPzfTPB1/r/0AJ3hfLREN4BAVxSSzltv6EBERRaMcYxbiVHrYXO2wOdukDoeIIgCbWooJzj4vGqx2tHb0AsC5xaL6uFgUERGRFIK1udM5oZ9TCArMME0DwFWQiagfm1qKCYeqbXj8lT3Y8lktAHCxKCIiIokdqLLi8Vf24K1ttRP+2SIzm1oiOodNLcWEc6sf92/wfm6vWja1REREUoj7XG2eiCLTDABAVUcN/AF/SOMioujDppZiwuc3eOdILRERkbR0A7W5dwJb+gxK0ZuRrDPB4e1FQ09TqEMjoijDppZiQrCpdX++qeU9tURERFIYHKl1TaKpFQQBRebpADgFmYjY1FKM+PwG76bg9GM2tURERFLQDdTmyYzUAkChaaCp7WBTSxTr2NRSTPj89GODOh5qhRpOrxNOr0vK0IiIiGLS52vzRBWapkGAgJqu03D7PKEMjYiiDJtaigl6zeAUJx/8gQAEQeB9tURERBIK1ma3D4FAYMI/b9DEI9uYCa/fi5quU6EOj4iiiErqAIjCQaEQ8Mht86HTnPuVN+uS0NLbyqaWiIhIAmqVAo/cNj84YjsZM80zUN/TiP+tfhe587MRr44LYYREFC04UksxozDXhCnpRigEAQAXiyIiIpJaYa4JuWlGCAO1eaJW5VyBtLhUNNqb8WL5y3B6nSGOkIiiAZtailmDe9V2sKklIiKKSkaNAd+Yfy9SdGac6WnEL8p/CxfXyiCKOWxqKWZ8uK8ev3vvGJrbHAC4Vy0REZHU/rG3vza3dPRO+hxJ2kR8Y/7XYNIm4VR3HX5V8SrcPncIoySiSMemlmLGoZo2bKtohrWz/xNcNrVERETSKj9pxbaKZti6Lm50NVlvwjfm34dETQJOdtbipcN/gIcrIhPFDDa1FDMGF6Jwufu3DjAH96plU0tERCSFYG2e5LY+50uNS8E35t8Ho9qAY+1V+M2R1+D1X/x5iSjysamlmBH3uU3eEzUJUAgKdLl7+GkuERGRBOIGmtreEDS1AJAen4qvz78X8ao4HGk7ht9Vvg6f3xeScxNR5GJTSzFDd95etQCgVCiRqEkAALQ5uVgUERFRuOm0Q2tzKGQZMvDQvHugV+lQbj2MPxzbDH/AH7LzE1HkYVNLMWO4T4MH76u1OdokiYmIiCiWDU4/doZopHZQbkI2Hpx7N7RKDfa1lOPPx//GxpZIxtjUUszQDVM4B5taq6NdkpiIiIhimf5ztwaFUl7iFPz7nLugVqhR1rwXf616C4FAIOTXISLpsamlmGE2apGbZkCSQRP8XvLAYlHWXja1RERE4ZacoENumgGJ59XmUJpuysf9c/4VKoUK2xrL8OvDv0ezo0WUaxGRdFRSB0AULguLUrGwKHXI985NP2ZTS0REFG6LZ6Zh8cw0Ua9RZJ6Oe2d/Fa8ceQ2HbUdxxHYMi9MX4Lq8NUjWm0W9NhGFB0dqKaYFm1qO1BIREcnW7JSZeGLJI1ietQSCIGD32f34r10/wV+r3kRXX4/U4RHRRWJTSzElEAjA7Tm3wqKJ04+JiIgk9fnaLJZEbQLWFX4RGy//HhanL4A/4MenDTvxRNnTeKvmffR6ekWPgYjEwaaWYkZLey/ueeZjbPztnuD3zANNbVtvB1dFJCIiCrMmmwP3PPMx/uvVvWG7Zoo+GXfOWo8Ni7+NuSnFcPs9+Efdx3i87Gl8cPoj9PncYYuFiEKD99RSzNBqlAgEhq5+rFFqYFDHw+5xoNvdgyRtooQREhERxRbdMLU5XDIN6bhvzp043X0Gf6/ZiuMdJ/H32g/wScN2rMq5AksyFsGoMYQ9LiKaODa1FDOCe+G5h05xMutMsHscaHd1sqklIiIKo5FqczhNTcjF1+ffixPt1Xir9n3UddfjrZr38U7tPzDPMhvLsi7H9KR8CIIgWYxENDo2tRQzNCoFlAoBHq8fXp8fKmX/7HuzzoQzPQ1ocbQiP3GKxFESERHFDq1GCQFAn9sHvz8AhUK6xrHQPA3fMz2Eyrbj2Na4C5Vtx7G/9RD2tx5CWpwFJZmX4bKMS2FQx0sWIxENj00txQxBEKDTKOFwedHb50VCXP+eeIWmApRbD+Oj+m24LONSKATeak5ERBQOCkGATquEs88Hp9uLeJ1a0ngEQcDslJmYnTIT7a4O7Gzai51Ne9DSa8WW6nfwdu0HmG+5BMuyLkdB4lSO3hJFCFGb2k2bNuHQoUMQBAEbNmzAnDlzgs81NzfjO9/5DjweD2bNmoUf/vCHYoZCBKB/mpPD5YXrvKZ2SeZi/LPhMzQ5zmJ/yyEsSp8vcZREROJhbaZIo9Oo+pvaPumb2vOZdSZcn78WX5h6FY60HcP2xt041l6FvS0HsbflINLj01CSsQiXps1DojZB6nCJYppoQ1J79uxBXV0dNm/ejNLSUpSWlg55/umnn8Zdd92FN954A0qlEk1NTWKFQhQUvHen79y9O2qFCl8qvg4A8M6pf8Dnl+6+HiIiMbE2UySKG6Y2RxKlQom5ltl4cN7deGLJI7h6yioYNQacdbTgb9Xv4Ac7SvHCwZdR1rwPTq9L6nCJYpJoI7VlZWVYvXo1AKCgoABdXV2w2+0wGAzw+/3Yv38/nnvuOQDAxo0bxQqDaIgbS6bC5fbBnKAd8v3lUy/DlsoP0NJrRVnzXizLulyiCImIxMPaTJHoX5blwe31wWTUjn2wxFL0ZtxYcA2uy1uDCttR7Dl7AJVtx3G84ySOd5zE5hNbMDtlFhalzUdxciFUCt7pRxQOov1Ls9lsKC4uDj42m82wWq0wGAxob29HfHw8fvSjH6GyshILFy7Ed7/7XbFCIQq6tDB12O8rFUpcl7cWv638E94//U8sTr8UGmXkTIEiIgoF1maKRAuLhq/NkUypUGJ+6iWYn3oJHJ5eHGytwN6Wg6juPIWDrRU42FqBOJUe81PnYFHafCSnXCJ1yESyFraPjwKBwJCvW1pacMcddyArKwv33XcfPvnkE1x55ZWjnsNiMYocZfgwl8iztngpPmr8FKc7G3Cw6wCuL1wtdUiTJpfXBGAukUgueRBr8+cxl8gTDXlYYMTUzDX4ItbA5mjH9jN7sb1uL850NWJH027saNqN5OMmLM6eh8uy56MopQAKRXQvShkNr8t4ySUXueQxWaI1tampqbDZbMHHra2tsFgsAACTyYTMzEzk5uYCAJYsWYKTJ0+OWTit1h6xwg0ri8XIXCRS09iF2uZuTMtKRF7GuUUdLBYj2mwOfCF3DX7Z+Tv8rfJ9zEmYC71KJ2G0kxNtr8lomEvkkUseQGz+AcDaPDK5/W5HUy7VDV041dyN6TmJmJo+tDZHUx791ChJWYqSlKVotDdj79mD2NdSjjZnB94/+THeP/kxjGoD5liKMd9yCWaYCqBUKKUOekKi83UZnlxykUsewORrs2gfE5WUlGDr1q0AgMrKSqSmpsJgMAAAVCoVcnJycPr06eDzeXl5YoVCFFRebcPrH57EkVPtwz5fnFyE/MSpcHh68XH9tjBHR0QkLtZmikQHqqx4/Z8ncex0h9ShhFSWIQM3TbsWP1z6KEpXfx+rc1cgRWdGj8eOHU278eKh3+CR7T/EH45uRoW1Eh6fR+qQiaKWaCO1CxYsQHFxMdavXw9BELBx40Zs2bIFRqMRa9aswYYNG/Doo48iEAhgxowZWLVqlVihEAXpNP2fhjr7vMM+LwgCbsy/Bs8f/BX+eeYzLM9eyk3WiUg2WJspEum0/bW5d4TaHO0UggLTk/OQNC0FNxVci0Z7M8qth3HQegRnHS3YfXY/dp/dD61Sg+LkIlySMguzkgv59wfRBIh6T+3DDz885HFRUVHw6ylTpuD1118X8/JEFxjcNsA1SuGcbsrHTPMMHGuvwod1n+KmadeGKzwiItGxNlOk0Qdrc2Ru6RNKgiAg25iJbGMmrs+/GmcdrSi3HkG59TDqexpxoLUCB1orIEBAfuIUzE6ZiUtSZiE9LhWCIEgdPlHE4jrjFFMGC+dYnwbfkH81jrVX4ZOGHbgypwRJ2sRwhEdERBRz4sZZm+UoPT4V18SvwjVTV8HmbMdh21EcsR3Dyc5a1HSdRk3XabxV8z6SdeaBBncmpiXlQ82tgoiG4L8Iiim6cW7wPiUhB/Mss1FuPYKtpz/CusIvhiM8IiKimKPTDNbm2Gtqz5eiN2NlzjKszFkGp9eFY+1VOGI7hsq242hztePThh34tGEHtEoNZppnoDi5CLOSC/nBOxHY1FKMGfw02Okeu3Ben381Dlkrsb1pN67KXY4UfbLY4REREcWcuOAocvgAACAASURBVIF7al3jqM2xQq/SYUHqHCxInQN/wI/T3fU4YjuGI23HBu7JPYJy6xEA/QtSFScXYZa5EPmJU6JuNWWiUGBTSzFFr1VBMc57UjLi07A4fQF2n92P9059iDtmrRM5OiIiotij142/NscihaBAfuIU5CdOwY0F16DN2YHKtuM42n4cJ9qr0WhvRqO9Gf+o+xh6lQ5FpunBUdxEbcLYFyCSATa1FFOyLfF4+ftXjnuxhWvz1mBfSzn2nD2A1bkrkGlIFzlCIiKi2DIlzTih2hzrkvUmLM9eguXZS+DxeVDddaq/yW07gZZeKw5aD+Og9TAAINuQiUsGFpvKMWZBIYi2myeRpNjUUkyZaMFM0ZtRkrkYnzWW4d1T/8C9l9whUmRERESxic3s5KmVasw0z8BM8wxgOmBztqGy7QSOth3HiY4aNNib0GBvwvun/4lETUKwwS00TYNaqZY6fKKQYVNLNIZrpl6FsuZ9KLceQV13PaYk5EgdEhEREdEFUvTJWJG9FCuyl8Lj86CqswaHbcdw2HYUnX1d2N60G9ubdkMzsNjUJSmzMDu5CEaNQerQiS4Km1qKOaV/3Icuuxv/ddfi4BY/o0nUJmBF9lJ8eOZT/L12Kx6ad08YoiQiIoodT/5+H3p63XjqnsugUXOho1BQK9UoTi5CcXIR1s24CQ32JlTYjuKw7SjqexpxyHoEh6xHIEBAXuIUXJo2F4vT5iNOHSd16EQTxqaWYo6ty4Uuuxsut29cTS0ArJlyJbY37sKx9iqc7KjBdFOByFESERHFDluXEz29HjjdPja1IhAEATnGLOQYs3Bd3hp0uDqDI7hVHdWo7TqN2q7T+N/qdzHfcgmWZi7CtKR83oNLUYNNLcWcOK0KXXY3evu8MBm14/oZgzoeq3KX471T/w8fnP6ITS0REVEI6TWq/qa2z4vEeI3U4cieSZcUXGzK5XXhSNtxlDXtxfGOk9jbchB7Ww4iRZ+MJRmLcHnGpdwLlyIem1qKOYObvLsmuMn7yuwSfHjmUxzvOIn6nkbkGLPECI+IiCjmDM6cck6wNtPF06l0WJg2DwvT5qHN2Y6y5n0oa94Lm7MNf6/9AO/UbsXslCIsyViM2clFUodLNKxxzSmorKwUOw6isBnc5H2ihTNOHYdlmZcBAD4882nI4yIimgjWZpIT/SRrM4VWst6M6/PX4sml/4EH5t6NeZZLIAgCDtuO4aXDv8djOzfhrWP/gM/vkzpUoiHGNVL73//93zh9+jTmzp2LK664AiUlJbBYLGLHRiQK3cCnwb2TKJwrc5bhk4YdONBagRvyr0GK3hzq8IiIxoW1meSEI7WRRSEoUJxciOLkQvS47dh9dj92Nu1FS28r/lTxv/jEsAu3F93CHSEoYoyrqf3Nb34Dr9eLw4cPY/fu3XjkkUfQ3t6ON998U+z4iEJusHC63BP/lNGsM2Fh2jzsOXsAH9Vvw60z/iXU4RERjQtrM8nJ4K1Bzj6OAEYao8aA1bkrcFXOchxtr8Ib1W+i0d6Mn+x7EVfmlOD6vKuhU41vjRIisYxr+rHH40F5eTnKysqwd+9e9PX14bLLLhM7NiJRXJKfjGsW5yIrJX5SP786dwUAYGfTHtjdjlCGRkQ0bqzNJCdzCpJxzWW5yJxkbSbxCYKA4uRCPHvN48G/hT6u347SPc+hsu24xNFRrBvXSO2iRYuwaNEi3HbbbbjzzjsRH883HIpei4pSsagoddI/n2XIwKzkQhxtO4HPGnfi2rw1IYyOiGh8WJtJTi6blYbLZqVJHQaNg1alwRenXYdL0+biz8feQL29Cf9z6LdYmDYPt0y/EUaNQeoQKQYpn3jiiSfGOqiwsBBOpxPvv/8+PvroIzQ3N0OpVCI9PT0MIZ7T2+sO6/XEEh+vZS4RZqJ5JGkTsPvsfjQ7WrAieymUisjZU08urwnAXCKRXPIA+nOJZqzNoSW332055CKXPAB55pKoTcCSjEXQqXSo6TyNBnsTypr3IkFjRJYhA4IgSB3qmOTyusglD2DytXlcI7UrV67EypUrAQDl5eX41a9+heeffx5Hjx6d1EWJpNTT60aD1YE4rQpT0o2TOsf0pALkGrNxpqcBu5r3YXn20hBHSUQ0OtZmkpPuXjcaW+2I16uRmza52kzhp1QosTp3BeZZZuP141twvOMk/nBsM/acPYAvF92MFH2y1CFSjBjXPbWvvfYavvnNb+Lqq6/Gyy+/jCuuuALvvvuu2LERieLEmU785PWD+PvO05M+hyAIWDPlSgDAP898xqXtiSjsWJtJTo6ebsdP/lKOd8vqpA6FJiFFn4yH5t2DO2auQ7wqDsc7TuKZfS+gpdcqdWgUI0Ztau12O1599VV0dHTgjjvuwFe/+lXU19dj9+7dMBr5KRpFp1BtGzDPMhspOjNsrnaUW4+EIjQiojGxNpMcxQ3WZje39IlWgiDgsoxL8Z+XP4wi03Q4PL34n/JX0OO2Sx0axYBRm9rHH38cbW1t+PrXvw6z2Yyf//znePTRR1FSUoLS0tJwxUgUUqFqahWCAlcNrP734ZlPEAgELjo2IqKxsDaTHHGfWvkwagy495I7kGvMgs3Vjl9W/A5unzzu96TINWpTW19fj+9+97sAgK1bt+Kaa67B0qVLsW7dOthstrAESBRqem3/ok6hKJyXZyyEQR2PMz2NqOqouejzERGNhbWZ5EjPfWplRafS4v45dyFZZ0Jddz1+V/k6/AG/1GGRjI3a1MbFxQW/3rNnDy6//PLg42hY0YxoOMFPg90XXzg1SjWuzC4BAPy/M59c9PmIiMbC2kxyxJFa+UnUGvHA3LsQp9KjwlaJN06+zVltJJpRm1qfz4e2tjacOXMGBw8eRElJ/x/vDocDTqczLAEShVqoC+cV2UugUahxrL0KDT1NITknEdFIWJtJjtjUylN6fBruu+ROqAQlPm3YiX/WfyZ1SCRToza19957L6699lrccMMNeOCBB5CYmAiXy4XbbrsNN910U7hiJAopjUoBhSDA4/XD67v4qTAGdTyWZi4GAHx45tOLPh8R0WhYm0mOdJr+W4Ncbh/8HM2TlemmfHx11joAwP9Wv4sDrRUSR0RyNOo+tStWrMD27dvR19cHg8EAANDpdPje976HZcuWhSVAolATBAFP3rMYOo0KSkVopuqtyrkCnzWWYX/rIdyQfw2S9aaQnJeI6PNYm0mOFAoBpfdeBp1GBU6il5+FafPQ4erEmzXv4fdH/4IEjRHTkvKkDotkZMx9atVqdbBoDmLRpGiXkRwPk1EbsvvPkvVmLEidA3/Aj48btoXknEREI2FtJjkKdW2myLI6dwWWZy2B1+/FSxW/R4ujVeqQSEbGbGqJaHxW514JANjRtAcOT6+0wRARERFFEEEQcMv0G3FJykw4vL34xaHfotvdI3VYJBNsaikm/X3nafz8jQqcaQndm2mOMRMzzTPg9rmxrbEsZOclIiKKBW9vP4Wfv1GBhla71KGQSJQKJf6t+HbkGrPR5mrHrw69ij7uYUshwKaWYlJtYxfKq21o63aF9Lyrc1cAAD6p3wG3zxPScxMREclZ9UBtbu8JbW2myKJVavDvc/8NyToz6nrq8bvKP3OrH7pobGopJg1uHeAK8SbvhaZpyDFmocdjx+6z+0J6biIiIjk7t61PaGszRZ4EjREPzr0LepUeh21HUdtVJ3VIFOXY1FJMGiycvSHeD08QBKzKuQIAsPfswZCem4iISM702v5tfbhXbWxIi09FycCWiHtaDkgcDUU7NrUUk3QiFs7ZyTMhQMCp7jNweftCfn4iIiI5OjdSy6Y2VixOXwAAONByCF4/X3eaPDa1FJPiBgunO/RvoHFqPXITsuEP+FHdWRvy8xMREcmRXiNebabIlGXIQGZ8Onq9TlS2nZA6HIpibGopJol9306RaToA4ERHtSjnJyIikptgbXbxntpYMjhau/cspyDT5LGppZiUmqRHcZ4ZGclxopy/0DQNAJtaIiKi8Uo16TE7z4x0kWozRaaFafMgQMDhtmPo9TilDoeilErMk2/atAmHDh2CIAjYsGED5syZc8Exzz77LMrLy/HHP/5RzFCIhpidn4zZ+cminT8/cQrUChUa7c3ocdth1BhEuxYR0USwNlOkmjstBXOnpUgdBoWZSZeE6aYCVHVUo9x6GEsHFo8imgjRRmr37NmDuro6bN68GaWlpSgtLb3gmOrqauzdu1esEIgko1aqUZCYB4CjtUQUOVibiSgSLU6bDwDYwynINEmiNbVlZWVYvXo1AKCgoABdXV2w2+1Djnn66afx7W9/W6wQiEbkDwRgd3rQ3i3eBu+F5oEpyO1saokoMrA2UyTzBwLo6XWjo4c7B8SaeamzoVaocLKzFu2uDqnDoSgk2vRjm82G4uLi4GOz2Qyr1QqDoX8a5pYtW7B48WJkZWWN+5wWizHkcUqFuUjrbJsD3/jZx0g16fHKY2sBhD6PJcq5eKvmfZzsqkZKigGCIIT0/KOJxtdkJMwl8sglj1jE2jw65iKtRqsd3/z5dmSkxOOl/+j/8CUa8xgJcxmNEQuz5qKsfj+O2Y/hppyrQ3z+kcnldZFLHpMl6j215wsEAsGvOzs7sWXLFvzud79DS0vLuM9htfaIEVrYWSxG5iIxp9MDALD3emC19oiSR3wgCXEqPay97Th2pg6WOPHu4T1ftL4mw2EukUcueQD8AwBgbT6f3H63ozEXp8MNALD3ukWrzVJhLmOba7oEZfX78XFNGZYmLwnLYIBcXhe55AFMvjaLNv04NTUVNpst+Li1tRUWiwUAsGvXLrS3t+P222/HQw89hMrKSmzatEmsUIguoNMoAfTvhXf+H3WhpBAUmBFcBfmkKNcgIpoI1maKZPrB2tzHfWpj0UzzDMSr49DsaEGDvVnqcCjKiNbUlpSUYOvWrQCAyspKpKamBqc3XXPNNXjvvffw17/+FS+++CKKi4uxYcMGsUIhuoBKqYBGrUAgAPR5xNsPb3Brn+NcLIqIIgBrM0UytUoBpUKA1xeAx+uXOhwKM5VChUtT5wLgnrU0caJNP16wYAGKi4uxfv16CIKAjRs3YsuWLTAajVizZo1YlyUaN71GBbfHDWefeE1t0cBiUVUd1fAH/FAI3BqaiKTD2kyRTBAE6LUq2J0ejtbGqEXpC/BZYxn2tRzETdOu5d9NNG6i3lP78MMPD3lcVFR0wTHZ2dncB48kodeq0OVwi1o4LfoUmLRJ6OjrRKO9GTnG8S++QkQkBtZmimR6rbK/qXWzqY1FeQm5SNGZYXO1o6qjBkXm6VKHRFGCH39QzNJr+z/TEbOpFQQhuLXP8XbeV0tERDSacNRmilyCIGBR+gIA3LOWJoZNLcWsL11ZgG/fOhcZyXGiXqfI1P8p4wneV0tERDSqdaum49u3zkVqkri1mSLXovT5AIBy62G4fW6Jo6FoEbYtfYgiTdEUU1iuM7gCcnXnKXj8XqgV/GdHREQ0nJlhqs0UudLiLJiSkIO67npU2I5iYdo8qUOiKMCRWiKRJWqNyIxPh8fvwemuOqnDISIiIopoi9P6pyBzFWQaLza1FLOOnW7Hls9qcKyuQ/RrBe+r5RRkIiKiEVUO1OYTZ8SvzRS5Lk2bC4WgwNH2KvS47VKHQ1GATS3FrBP1nXhnZ11YCufgfrUnuFgUERHRiI7XdeCdnXWoqu+UOhSSkFFjwCzzDPgDfuxvPSR1OBQF2NRSzDq3wqJ4+9QOmp6UD4WgQF1PA5xep+jXIyIiikbhrM0U2QZXQd579qDEkVA0YFNLMSuc2wboVDpMTciBP+DHyY5a0a9HREQUjfQaJQCgl1v6xLw5KbOgVWpwuvsMWnutUodDEY5NLcWsYFMbpg3eC7m1DxER0agGa7MrTLWZIpdGqcE8yyUAOFpLY2NTSzFLr+3/NDhcG7wXmfubWi4WRURENLzBppYjtQQAiwemIO9pOYhAICBxNBTJ2NRSzNJrwjf9GACmJuRAo1DjrKMFXX3dYbkmERFRNAnnrUEU+WaYCpCoMcLmbMPp7jNSh0MRjE0txaw4nQqJBg0Mek1YrqdSqDDNlA+AU5CJiIiGE69TIcmggUGnljoUigAKQYGFafMBAHu4Zy2Ngk0txayM5Hj890PL8O1b54btmoNb+xzn1j5EREQXyLIY8NxDy/DNL4WvNlNkG1wFeX/rIXj9HMGn4bGpJQqjovMWi+K9IURERESjyzZkICM+DQ5PL47YjkkdDkUoNrVEYZRpSIdBHY/Ovi4uT09EREQ0BkEQsDRzMQDgs8YyiaOhSMWmlmLaD17ehfuf/QTdDndYrqcQFMEpyLyvloiI6EKP/roM9z/7CXpdHqlDoQhxefpCqBVqnOioxllHq9ThUARiU0sxze3xw+3xh7VwFpoH7qtlU0tERHQBt8cHt8cPh5P3T1K/OLUeiwYWjNrG0VoaBptaimnB/fBc4SuchQP31VZ11MAf8IftukRERNHg3F61HKmlc5ZnLwUA7Grejz5feGbYUfRgU0sxTa9VAgAcYRypTdGbkaIzw+l1or6nMWzXJSIiigbBppYjtXSeHGMm8hKmwOVzYS+396HPYVNLMS24yXsYR2oBoNDcP1rLrX2IiIiG4kgtjWR59hIA/QtGcRcJOh+bWoppcQOFM5wjtQC4WBQREdEIOFJLI5mfOgcGdTwa7c2o7aqTOhyKIGxqKabpgoVTmqa2pus03D5+Ek1ERDRIrwn/rUEUHdQK1Xnb++yUOBqKJGxqKaYtmJGCL60sQNFUc1iva9DEI9uQCa/fi9qu02G9NhERUSS7tDAVX1pZgMIpJqlDoQi0LPNyCBBwsPUwut09UodDEYJNLcW02XnJ+MJlU1CQnRT2aw9u7cMpyEREROfMKeivzXmZiVKHQhEoWW/C7JSZ8AV82Nm0V+pwKEKwqSWSSNHA1j77Ww7B5e2TOBoiIiKi6LA8q3/BqO2Nu+Dz+ySOhiIBm1qKaR09fdhzrAUV1dawX3uGqQAZ8Wloc7XjjZNvh/36REREkai924U9x1pwuNomdSgUoYrM02HRJ6OjrxNH2o5LHQ5FADa1FNNOn+3Gr96qxJuf1oT92iqFCncV3w61QoWy5r3Y31Ie9hiIiIgiTW1Tf21+e1v4azNFB4WgwBUDo7WfNXDBKGJTSzFOrxlY/TjM+9QOyjSk4+ZpNwAA/nx8C2zOdkniICIiihTBLX0kqs0UHZZkLIRaocbxjpNo6Q3/jDuKLGxqKaadK5zSbRtwRdblmJtSDJfPhVcrX+e9IUREFNMioTZT5ItTx2Fh2jwAwLbGMomjIamxqaWYptf1F06HhJ8GC4KA22begiRtIk511+G90x9KFgsREZHU9NrBfWo5UkujW57dPwV5V/M+9PncEkdDUmJTSzFtcIN3p8SfBhvU8bhz1noIELD19Eeo6uB9REREFJsGR2qdbGppDLnGbExNyIXT68K+loNSh0MSYlNLMW2wcDpcXgQCAUljmWEqwNVTVyGAAH5/9C+wexySxkNERCSFc7WZ049pbMuDC0aVSf63HEmHTS3FNJVSAbVKAb8/ALfHL3U4uHbqauQlTEFnXxf+dOwNvjkTEVHM0agUUAgCPF4/PF7pazNFtgWpc2BQx6PB3oRT3WekDockwqaWYt6zD5Zgy49vgHZgKrKUlAol/q34y9CrdKiwVXLhAyIiijmCIOD5byzDlh/fALWKf6rS6NRKNZZkLALQP1pLsYnvFBTzDHp1RBXNZL0ZXy68GQDwt+p30GhvljgiIiKi8Iq02kyRbVnW5RAg4GDrIfS47VKHQxIQ9d1i06ZNWLduHdavX4+Kioohz+3atQu33nor1q9fj//4j/+A38/pJUSDLk2bhyUZi+D1e/G7yj/DzRX9iChEWJuJSG5S9GYUJxfBG/ChrGmv1OGQBERravfs2YO6ujps3rwZpaWlKC0tHfL8448/jp///Of4y1/+AofDgW3btokVCtGo3vikBg//7DNUN3RJHcoQX5rxL0iLs6DZ0YK/Vb8jdThEJAOszRQt/vpRNR7+2WeobeqWOhSKEoPb+3zWWAavnytnxxqVWCcuKyvD6tWrAQAFBQXo6uqC3W6HwWAAAGzZsiX4tdlsRkdHh1ihEI2qpb0XJ8504P3ddUg/GRf8/vScJMyblgIAaO3oxaflTSOe4+rLcpEQpwEA7Dp6FvUtw099sSTpceX8LACAx+vHm9tqRzzn4plp+Lfi2/DTfS9ie+MuNDcCOhjghx9++BAY+C8UAcycmgif3wtvwIczu3rg6B1+ZNdk1CHVpAfQv6l9fevIKyznZSRAq+6/z7ipzY5ux/CrUMZpVchNMwIAfP4ATjZ0jnjOjOQ4JMZrAQDt3S60djqHPU4hCJiRkwR9kwbOXjdqm7tGXMjLZNQizdT/ujlcHtS3jjztSMqctBoV8tKNwcfRmpNarYTH4wNw7nWK1pyeXXf/iMfLFWszRYuzg7V5V12wbgFAYa4JcwqSAfTX788OjVybv3D5FBj0agBA2ZGzaLAO/76TatJjxbz+2tzn8eHt7adGPOflxenISe3/N3KsrgNHatuGPU6tUuCmK/KDj9/46CRabcNfP1pzem9XHRzO4d+fpcgpEFBCByM6+jrxzU82QK1QQavUQhFQw+dVQAn15/6ngk6jRm6qAQIEQABON/ectziZACF4dgHmBB1Sk/p/F8eqY/mZCdCo+utYo82Ont4R6phOhdzU/jrm9ftHHWDJTIkP/q3Z1u2CdYTarFQImJ59rjbXNHWNuOBaJOY02dosWlNrs9lQXFwcfGw2m2G1WoPFcvC/ra2t2LFjB775zW+OeU6LxTjmMdGCuUSODIsBqLLi4EnbkO9rdWqsWZIHADjb1Yf3d4+8ot6NV04L/v9wvP4EPjvYOOxxswuS8aW1RQD6m8rRzjlrmgVX5hfhds8X8fvyN1CDXcMf6AeqR+6Nh6jvAdAzzmPHeytvH3BinIsN1rcBGL5WDhUA6sZ7zmjJyQ1UyyGn82uYHF8nmWNtHh1ziRwZFgPKq23YX2Ud8v24OA2uunwqAKCxwzlqHb35qhmwJMcDAI6eOY4dFcM3VnOnp+CWNf212d7rHvWccwpTg//fflLRPOKx8Xo17r15bvDxB2W70NLeO+yx0ZiTxWLEtormiMtJmZwH9ZRjEFQeePxeeMYasXUDtQ2jHzKovhvAOCcO1I/cww/lmkAds43y5Pm12Q+cHu85IzmnCRKtqf284bYmaWtrw/3334+NGzfCZDKNeQ6rdZx/EUU4i8XIXCLI6gVZSE+JR2fX0BG2vPRzuakRwC1XFox4Do/LEzx2br4ZqYm6YY8zG7XB4zxe/6jnTNKrYLX2YJFpEY4lN+FURyMUggIKKM/9T1BCJShQkGmGUlBCpVDhbLsLPY6+Ea+fOoHRMo26/w6FJptjxE/k9FrlkNGy0T6Ry0iOQ0J8/ydy7d19Qz5lPJ9C6B8pj4/XwuHoQ21T94ifMpqM2vNGn70Rm5NOq8KUNEPwcbTmpNYo4XEPjtT2v07RnlMsY20+Rw71bJAccll7aRYyLQZ0dQ/9t5p/Xm3WChi1jvb19sE6cF/4/GnJyDDrhz0uOUF3Xm32jXrOBK0yeGxOctyIx6qUiiGvwf9ZOQ2tbcPPjoq2nAZ/v9YuyoGzb/imUbqcCgBcgXnTkpFi0sDl68PJxjZUNdvgC3jghSf4X2/AA71ehalpBmDgvfBkYxc8Xh/6Hw19fzQbtbAM1DGH04MG68iz3fIzEqAeqGONVgfsI4xox2lVyBn422CsOpaZEgfj4Khmlwu2LlfwuaG1WcD0nMTgczVN3fCONFIbwTlNlBAQaSPMF154ARaLBevXrwcAXHXVVXjrrbeCnwLb7Xbccccd+Na3voXly5eP65zR/gY9SA7FZpBccpFLHgBziVRyyUUueQDRP5I1GazNI5Pb77YccpFLHgBziVRyyUUueQCTr82iLRRVUlKCrVu3AgAqKyuRmpoaLJoA8PTTT+POO+8cd9EkIiKii8PaTEREciTa9OMFCxaguLgY69evhyAI2LhxI7Zs2QKj0Yhly5bhzTffRF1dHd544w0AwPXXX49169aJFQ4REVHMY20mIiI5EvWe2ocffnjI46KiouDXR44cEfPSRERENAzWZiIikhvRph8TERERERERiY1NLREREREREUUtNrVEREREREQUtdjUEhERERERUdRiU0tERERERERRi00tERERERERRS02tURERERERBS12NQSERERERFR1GJTS0RERERERFGLTS0RERERERFFLTa1REREREREFLXY1BIREREREVHUYlNLREREREREUYtNLREREREREUUtNrVEREREREQUtdjUEhERERERUdRiU0tERERERERRi00tERERERERRS02tURERERERBS12NQSERERERFR1GJTS0RERERERFGLTS0RERERERFFLTa1REREREREFLXY1BIREREREVHUYlNLREREREREUYtNLREREREREUUtNrVEREREREQUtdjUEhERERERUdRiU0tERERERERRi00tERERERERRS02tURERERERBS12NQSERERERFR1GJTS0RERERERFGLTS0RERERERFFLTa1REREREREFLXY1BIREREREVHUErWp3bRpE9atW4f169ejoqJiyHM7d+7ELbfcgnXr1uEXv/iFmGEQERHRANZmIiKSG9Ga2j179qCurg6bN29GaWkpSktLhzz/1FNP4YUXXsDrr7+OHTt2oLq6WqxQiIiICKzNREQkT6I1tWVlZVi9ejUAoKCgAF1dXbDb7QCA+vp6JCYmIiMjAwqFAitWrEBZWZlYoRARERFYm4mISJ5UYp3YZrOhuLg4+NhsNsNqtcJgMMBqtcJsNg95rr6+fsxzWixGMAnw1AAAIABJREFUUWKVAnOJPHLJA2AukUouucglj1jE2jw65hJ55JIHwFwilVxykUsekxW2haICgUC4LkVERETjwNpMRERyIFpTm5qaCpvNFnzc2toKi8Uy7HMtLS1ITU0VKxQiIiICazMREcmTaE1tSUkJtm7dCgCorKxEamoqDAYDACA7Oxt2ux0NDQ3wer34+OOPUVJSIlYoREREBNZmIiKSJyEg4tyjn/70p9i3bx8EQcDGjRtx9OhRGI1GrFmzBnv37sVPf/pTAMDatWtx9913ixUGERERDWBtJiIiuRG1qSUiIiIiIiISU9gWiiIiIiIiIiIKNTa1REREREREFLVE26f2YmzatAmHDh2CIAjYsGED5syZE3xu586deO6556BUKrF8+XI8+OCDEkY6utHyWLVqFdLT06FUKgH03+OUlpYmVahjqqqqwgMPPIB//dd/xVe+8pUhz0XTawKMnku0vS7PPPMM9u/fD6/Xi6997WtYu3Zt8Lloel1GyyOaXhOn04lHH30UbW1t6OvrwwMPPICVK1cGn4+m12SsXKLpdQEAl8uF66+/Hg888ABuvvnm4Pej6TWRGmtz5GFtjkyszZGFtTkyXxcgxLU5EGF2794duO+++wKBQCBQXV0duPXWW4c8/4UvfCHQ1NQU8Pl8gS9/+cuBkydPShHmmMbKY+XKlQG73S5FaBPmcDgCX/nKVwKPPfZY4I9//OMFz0fLaxIIjJ1LNL0uZWVlgXvuuScQCAQC7e3tgRUrVgx5Plpel7HyiKbX5N133w289NJLgUAgEGhoaAisXbt2yPPR8poEAmPnEk2vSyAQCDz33HOBm2++OfC3v/1tyPej6TWREmtz5GFtjkyszZGHtTlyhbI2R9z047KyMqxevRoAUFBQgK6uLtjtdgBAfX09EhMTkZGRAYVCgRUrVqCsrEzKcEc0Wh7RRqPR4OWXXx52v8Joek2A0XOJNosWLcLPfvYzAEBCQgKcTid8Ph+A6HpdRssj2lx77bW49957/z97dx4fd13tf/w1a2aSmeyTfW3TNS3dKaVsIiBXBJErUOUHKojXC+i9AiriAi5wwYt67wURRUUF5BYVK16VqsjefW/TNWmzr5N1Jvsyvz/STFvokqQzme39fDx4tJOZ+c45fNOcnO/3swDQ0NBwwtXRSDoncPpcIk1FRQXl5eVccsklJ3w90s5JKKk2hx/V5vCk2hx+VJvDU6Brc9gNP3a73ZSWlvofp6am0tLSgsPhoKWlhdTU1BOeq6mpCUWYZ3S6PMY88MAD1NXVsWTJEu655x4MBkMoQj0js9mM2Xzyb5VIOidw+lzGRMp5MZlMxMfHA/Db3/6Wiy66yD/cJJLOy+nyGBMp52TMqlWraGxs5KmnnvJ/LZLOyfFOlsuYSDkvjz76KF//+tdZs2bNCV+P1HMSCqrN4Ue1OTzPi2pz+FJtDi+Brs1h19S+my9Kdhx6dx6f//znufDCC0lKSuLOO+9k7dq1XHnllSGKTsZE4nn5+9//zm9/+1t+/vOfhzqUs3KqPCLxnPzv//4v+/bt44tf/CIvv/xy2BaU8ThVLpFyXtasWcPChQvJz88PdShRRbVZplIknhfV5vCj2hw+glGbw274cUZGBm632/+4ubkZl8t10ueamprCdqjK6fIAuPbaa0lLS8NsNnPRRRdx8ODBUIR51iLpnIxHpJ2Xt956i6eeeoqnn34ap9Pp/3qknZdT5QGRdU727NlDQ0MDAHPmzGF4eJi2tjYg8s7J6XKByDkvr7/+Oq+++io33HADv/nNb3jyySdZt24dEHnnJJRUmyNLJJ2T8Yi086LaHF5Um8NPMGpz2DW1K1euZO3atQCUlZWRkZHhHxaUl5eH1+ultraWoaEhXnvtNVauXBnKcE/pdHl4PB5uu+02BgYGANi8eTMzZswIWaxnI5LOyZlE2nnxeDx897vf5cc//jHJycknPBdJ5+V0eUTaOdmyZYv/arbb7aanp4eUlBQgss4JnD6XSDov//Vf/8Xvfvc7XnzxRa6//nruuOMOzj//fCDyzkkoqTZHlkg6J2cSaedFtTn8qDaHn2DUZoMvDMcQPfbYY2zZsgWDwcADDzzA3r17cTqdXH755WzevJnHHnsMgCuuuILbbrstxNGe2uny+OUvf8maNWuIi4tj7ty5fP3rXw/bYRB79uzh0Ucfpa6uDrPZTGZmJpdeeil5eXkRd07OlEsknZfVq1fz+OOPU1xc7P/a8uXLmTVrVkSdlzPlEUnnpK+vj69+9as0NDTQ19fHXXfdRUdHR0T+/DpTLpF0XsY8/vjj5ObmAkTkOQk11ebwotocnudFtTn8qDaH53kZE6jaHJZNrYiIiIiIiMh4hN3wYxEREREREZHxUlMrIiIiIiIiEUtNrYiIiIiIiEQsNbUiIiIiIiISsdTUioiIiIiISMRSUysiIiIiIiIRS02tiIiIiIiIRCw1tSJR4itf+QqPP/44AJWVlXzgAx+grKwsxFGJiIjELtVmkalh8Pl8vlAHISJnr6mpieuuu46nn36ae++9l29961ssXbo01GGJiIjELNVmkalhDnUAIhIYmZmZXHvttdx00038z//8j4qmiIhIiKk2i0wNDT8WiRKtra28+eabxMfHk5OTE+pwREREYp5qs8jU0PBjkSjQ1dXFJz/5ST772c/S2trKG2+8wVNPPRXqsERERGKWarPI1NGdWpEI19vby7/8y7/wsY99jCuuuILrr7+eI0eOsGHDhlCHJiIiEpNUm0Wmlu7UioiIiIiISMTSnVoRERERERGJWGpqRUREREREJGKpqRUREREREZGIpaZWREREREREIpaaWhEREREREYlYampFREREREQkYqmpFRERERERkYilplZEREREREQilppaERERERERiVhqakVERERERCRiBbWpPXjwIJdddhnPPffce55bt24dH/3oR7nxxhv54Q9/GMwwRERE5CjVZhERiTZBa2p7enr49re/zYoVK076/He+8x0ef/xxXnjhBd555x3Ky8uDFYqIiIig2iwiItEpaE2t1Wrl6aefJiMj4z3P1dTUkJSURHZ2NkajkYsvvpj169cHKxQRERFBtVlERKKTOWgHNpsxm09++JaWFlJTU/2PU1NTqampOe3xblj9rwGNT2LHwvRFfOXS2zEYDKEORUQkpFSbJZwUJuVyYdFyLihcRqo9OdThiEgEC1pTG2iXZPwT7V09J30uyRFHVmo8AH39Q1Q1eU55nMJMJ7a40bQb23ro9Paf9HU2q4nCrET/4wPV7ac8ZmZqPMmOOAA6vP00tZ08ToBZBSk4HHF4vf1UNXbRNzAc0TnZ4swUZjr9j8Mpp4N1bfQk72WHezvPbfkDVxa9/5Sf53I5aWk5dTyRRLmEp2jJJVrygNFc5OxcnPEBOrr6TvpcsiOOzKM/83v7h6g+zc/8oiwncdajP/Nbu+nsHjjp6+xxZgqO1hyfz8fBmo5THjMrNZ6ksTrm6aep/eR1zGCAmfnHanNlQxf9gyevY5GSk81moTDD4X8cTjkdrO3EM9xKQnYLVZ11VO18ied3/p5ZKSUsz17CAtc84kxWIPp+3iiX8BMtuURLHjD52hySpjYjIwO32+1/3NTUdNKhUMe7433XjP9kTR9nIHnjfB1w8QReeyb+b7yJHDNMc3rPP6IwyunHL5exucJM3Mzt/PHwWtLtaSzNXDiBAEVEYsdkavOd77s2pLX5kvzxv/ZMxurZhGpjmOb07tocTjk9+fvdbDncwvULZhPvamNj4zb2uPexv/0Q+9sPYTVZWeiax/KsJaSlqWaLyPiEpKnNy8vD6/VSW1tLVlYWr732Go899lgoQpEoNq84FYd9CZaMNN5o/jvP7nuRVFsy05KKQh2aiEjYUW2WqTB/WtroyK0UB8WuHBa45tE92MO25p1satzG4c4qNjVuY1PjNl48lM7nFnyGFJuGJovI6QWtqd2zZw+PPvoodXV1mM1m1q5dy6WXXkpeXh6XX345Dz74IPfccw8AH/zgBykuLg5WKBKjVs7PZuX8bHy+Gfgs3bxZt54f7/ol9y65C1d8WqjDExGZcqrNEmoXLsh5z9cSLPFcmLuCC3NX0NzjZnPjNtY3bKGp283b9Ru5etoHQhCpiEQSg8/n84U6iPGKprHiymVqDY8M89TuX7C39QCZ8S7uXXIn8ZZ4//ORksd4KJfwFC25REseoDm1gRJN3w/KJXzsce/jR7ueocCZx5eXfT7U4Zy1aDgnY5RL+ImWPGDytTloW/qIhJq3d5Dy2k7q3N2YjCZuLb2JnIQsmnpa+MnuXzE0MhTqEEVERGKKp2eA8tpO6t3dp33dzJTpWIxmqj21eAa8UxSdiEQqNbUStXZVuHn4ua38aX0lAHazjTsW3EqS1cmhjsP8ev/viKCBCiIiIhFvx6HR2vyXjVWnfZ3VZGVuxgwA9rUdnIrQRCSCqamVqGU/uiVQb9+xO7IptmQ+e86nsBotbGzcytqqf4QqPBERkZjjr839J99i6HgLs0oB2Nt6IKgxiUjkU1MrUctuHSucJw4zLkjM45OlH8eAgT8eXsuWxu2hCE9ERCTmHGtqzzwFaGH2aFO7r+0gI76RoMYlIpFNTa1ELX/hHHjv1eAFrlKuK7kKgGf3/4YD7oopjU0kVjQ01HPbbTef9LnWVjff/e5DAOzYsY329jYA7rvv7lMeb9u2LXzta18KfKAiMiUm0tTmODNJtaXgHeymxlMX7NBEYkY01mY1tRK17HEm4NSF8335F3JR7gqGRob47ttP0d7XMZXhicS8tLR0vvSlrwLwpz+97C+cjzzy/VCGJSJBdKbafDyDwcDctFmAhiCLTJVIrc1qaiVqnelqsMFg4KMzrmFmSgmefi+bmzQMWSRYHnroQX784x9y9913cdNNH+XAgf3+K8WbN2/grbde5+GHv0VjYyNXXfV+ADZv3si//MunuOuuz/CVr9zD4OBgiLMQkbN1ulFUJzM3dbSpLVNTKxJw0VSbzaEOQCRYjl+MwufzYTAY3vMak9HEiuylHGwv53Bn5RRHKDK1bn3k1Auj3XLlLC5ZmAvA6zvq+NUrp/4F8uf3XTqpzx8YGOD733+CNWt+yyuv/IkbbvgYAMuWnUdJyUzuvvtLZGVl+V/v8Xh44IHvkJOTy7e//Q02blxPfHz8qQ4vIhFgIsOPAWalTMdkMFHZVU33YA8JFv0MkOii2hwYamolaplNRh781DJ/AT2VaUlFABzurDpl8ysiZ2/BgkUAuFyZ7N1bdsbXJycn8+ij32F4eJj6+jqWLFkWFoVTRCbPah5fbR5jM9uYnlTEwY4K9rcdZEnmwiBHKBJboqU2q6mVqFaQ6Tzja9JsKaTYkmjv66Spp4WshIwpiExk6o33Ku4lC3P9V4bfzeVy0tLimdTnm0wm/9/Hs0f0f/zHt/nP//wvioqK+f73H53UZ4pIeDEYDOOqzcebmzaLgx0V7G1VUyvRR7U5MDSnVmKewWBgVvp0ACo6j4Q4GpHYZDQaGR4+cY5dd7eXzMwsPB4P27ZtDZt5OyIytfyLRbUdGNcv3SISGJFUm9XUSlT747pKnvz9bhpau0/7ulnp0wA43FE1FWGJyLssXLiYr33tyxw+fGx7reuuu55//dfb+O53H+Kmm27hued+QWurO4RRikggvPz2EZ78/W6a2nrG9fqchCyS45LoGvBQ620IcnQiMiaSarPBF0GXvCZ7Wz3cnM0QgXAT7rn85wvb2VfVzj2rFlJalHrK13WZ2vjK3x4hw57OAysiew/McD8nE6Fcwk+05AGjucjZi6bvB+UyNR55fhsHazr40scWMbsw5ZSvOz6P5/f9hnUNm/nwtH/iiqL3TVWoARPu52QilEv4iZY8YPK1WXdqJar5V1nsO/0qi4XJeViNFpp73XgGvFMRmoiISEyKn+AKyABzjg5BLmvbH5SYRCSyqamVqGa3jm+Td7PRRFFiAYC29hEREQkiW9xobe6ZQFM7O2UGRoORw51V9A71Bis0EYlQamolqk1kk/dpyUUAVHRUBjEiERGR2Ga3jtbmvnHU5jHxFjvFiQWM+EY40F5x5jeISExRUytRbSKbvB/br7YyiBGJiIjEtrHaPJE7tXDcKsitGoIsIidSUytRbWJNbQEGDFR76hgYDo/lyUVERKKNPW58U4Pe7VhTe1Bb+4jICdTUSlRzJduZlZ9MepLtjK+1m+3kOLIY9g1T7amdguhERERiT0ZK/Lhr8/HyHDk4LQ7a+zto6G4KUnQiEonMoQ5AJJiWzHKxZJZr3K+fllREnbeBwx2VlCQXBzEykdi1Y8c2CguLSEk59TZbJ/ODH3yXiopyAPr7+3A4nPzgBz/kr3/9Cy+++AIGg4EPf/gjfOhD1wYjbBEJkGWzM1g2O2PC7zMajMxNm8XGxq3sbTtAjiMrCNGJxKbJ1ubWVjcPPfRN+vv7SElJ4f77HyQ+Pp7Nmzfyk5/8EKPRxIoVK/nkJz8dpMhH6U6tyHGmJRUCUKF5tSJB86c/vUx7e9uE3/eFL3yJJ574CU888RNWrLiAq6/+CL29vTzzzNP81389yRNP/JjVq39NV1dnEKIWkXAwN3UmAPtaD4Y4EpHoMtna/Oyzv+DCCy/mhz98mgsuuJjf/vZ/Afjv/36M73znu/zoRz9j06YNHDlyONAhn0B3aiWq+Xw++gaGGRweITHeesbXTz9usagR3whGg677iEzWn//8R3bt2kFHRzvV1VV8/OM3k5mZxVtvvc6RI4f5zne+S1bW6J2Wdeve5te//tUJ77/mmuu44oor33Pcrq4utm7dzKc+dTvbtm1hzpxSHA4HAPPnL2DXrp1ccMFFwU9QRCZlorX5eLNTZ2LAQHnHYfqG+rGZ44IUpUh0CnRtrq2t5sorrwJg+fIVfP3r9/H+91+B05lIZubocVasWMnWrZsoLp4WtLzU1EpUq2/t4es/3Uh2WjwP3X7eGV+fakshyZpI50AXTT0tZCdkTkGUIsH35M6fUxbgFUNL02Zzx4JbT/uaiopynnrq59TW1vDAA/fzi1/8mpKSmdx995f8RRPg/PMv4PzzLxjX5/7xj7/ngx+8GoPBQGtrK8nJyf7nUlJSaW11Ty4hEZkSNc1eHnxmM3muBL512/IJvddhTaAwMZ/KrmoOdVQwP31ukKIUCb5oqM3TppWwfv3bzJ49hw0b1tHR0U5bWyvJySn+16SkpFBXV3d2iZ2BbkNJVLNbJ7bCosFgYPrR/WoPa79akbM2b945mEwmXK4Muru9ATnm3/62lssuu+Kkz2lFVJHwd2xngvHvU3u8sSHIe1sPBCwmkVgSyNp8882forLyCHfd9Rna2lpPWoenojTrTq1EtckUzmlJRWxr3kVFZyUrcyd2BVkkXJ3pqu14uVxOWlo84369yWTy//10Ded4hx/X1FSTlJRMXNzoqqnp6em0trb6n3e7WygtnT/u+ERk6k1ku72TmZs2mz9X/p2y1gP4fD4MBkMgwxOZMtFQm51OJ9/85sMAVFdXsnXrFtLTXbS1HavNLS3NpKenjzu+yVBTK1EtzmrCAPQPDjM8MoLJeObBCcfPqxWRwDMajQwPn3ihabzDj/ft20tJyQz/49LSeTz66HfweDyYTCZ27drJ5z9/T8BjFpHAsY2NohoYmlRTWpiYR4Ilnta+Nlp63WTEj3+XAxE5ucnW5pdf/j0jI8Nce+1H+dOf/sjKlReSnZ1Dd3c3DQ31uFwZrFv3Nt/4xreDGb6GH0t0MxoM2I5eEe4bGN/d2lxHNlaTlZbeVroGxn/VS0TGZ+HCxXzta1/m8OGKCb+3tdVNSsqxeTpxcTY++9m7uPvuu/j3f7+DW2+93b9olIiEJ7PJiNVixOcbveg8UUaDkTlHhyCXaQiySEBMtjZfeOHF/O1va/nMZz5Jc3MTH/7wdQDce+99PPjgV7nzztu59NLLKSgoDEbYfgZfBE1Amsht9XA20SEC4SwScvnik+/Q2tXPdz+7gvRk+0lf8+48/nv7TzjYXs7t825mYUZkDWWMhHMyXsol/ERLHjCai5y9aPp+UC5T5wtPvE2nd4Dv3bmSFOfJVzA+XR4bG7byq32rmZs2izsX3BbMUAMiEs7JeCmX8BMtecDka7Pu1ErUG7tT2zOBuTtjQ5C1X62IiEjg2a0Tr83Hm5M2eqf2UHsFA8ODAYtLRCKT5tRK1LvhfSUMDY2QnmQb93uOzautClJUIiIisWvV+0sYHvaR4pjcPrOJVif5zlxqPHWUdxxmbtqsAEcoIpFETa1EvfnT0ib8nqKkAgwYqPHUMTA8iNVkCUJkIiIisemc6We/Eurc1FnUeOrY23ZATa1IjNPwY5GTsJtt5DiyGPYNU9VVE+pwRERE5F3GGlntVysiamol6u2rbOPld45QXtc5ofdpXq2IiEhw7D1amyvqJ1abj1ecWIDdbKOppwV3b1sAoxORSKOmVqLe7iNtrHnrCAdrOib0Pu1XKyIiEhy7KlpZ89YRDtVMvqk1GU3MShndt7qsdX+gQhORCKSmVqKefWyT9wmusDgtuQgYXSxqxDcS6LBERERilm2StfndFrhKAVhXv4kI2qVSRAJMTa1EPfvRLX0mWjhTbSkkxyXRO9RLY3dzMEITERGJSfFjtXng7JraRa75OCwJ1HrrOdKlHQtEYlVQm9qHH36YG2+8kVWrVrFr164Tnnv++ee58cYb+djHPsZDDz0UzDAkxk22qQUNQRaR6KPaLOHAdha1+XgWk4Xzc84F4I3adWcdl4hEpqA1tZs2baKqqorVq1fz0EMPnVAcvV4vP/vZz3j++ed54YUXqKioYMeOHcEKRWLcsaZ2eMLvnabFokQkiqg2S7gYu1PbN4na/G4X5JyHAQPbm3fTNeA56+OJSOQJWlO7fv16LrvsMgCmT59OZ2cnXq8XAIvFgsVioaenh6GhIXp7e0lKSgpWKBLjzuZO7bTkQgAOd1QGMiQRkZBQbZZwYYsbnVPbc5Z3agHS7CnMS5/DsG+YdfWbzvp4IhJ5zME6sNvtprS01P84NTWVlpYWHA4HcXFx3HnnnVx22WXExcVx1VVXUVxcfMZjulzOYIU75ZTL1OnqH8YZb8HpiDttrCd7LjVtJrbtcbj72rA4Rki2R8YveOF+TiZCuYSfaMkjFqk2n55ymTq5vUM44y0kTqI2n8w1pe9n9xt7WdewiY8vuRqT0RSoUAMm3M/JRCiX8BMteUxW0Jradzt+RTqv18uPf/xjXnnlFRwOB5/4xCfYv38/s2fPPu0xWlqiY0iJy+VULlMoMc7Ef3/+QuDU30Ony6PIWcD+9kNsOlzGooz5QYszUCLhnIyXcgk/0ZIH6BcAUG0+XrR9b4d7Lil281nV5nfLMuaSYU+nudfNP/ZvYqFrXsBiDYRIOCfjpVzCT7TkAZOvzUEbfpyRkYHb7fY/bm5uxuVyAVBRUUF+fj6pqalYrVaWLl3Knj17ghWKyFmZlnR0CLLm1YpIhFNtlmhlNBi5MG8FAG9qwSiRmBO0pnblypWsXbsWgLKyMjIyMnA4HADk5uZSUVFBX18fAHv27KGoqChYoYiclbH9ais0r1ZEIpxqs0Sz87KWYjVaONBerq34RGJM0IYfL168mNLSUlatWoXBYOCBBx7gpZdewul0cvnll3Pbbbdxyy23YDKZWLRoEUuXLg1WKCJ89ekNdHUP8NidK4mzTGyeTXFiAQYM1HjrGBgewGqyBilKEZHgUm2WcOHz+fjKTzbQ3TvIDz53AWbT2d9nibfYWZa1iHfqN/Fm3XpumPnhAEQqIpEgqHNq77333hMeHz8vZ9WqVaxatSqYHy/i1907SHffEH39QxNuam1mG3mObGq89VR21TAzZXqQohQRCT7VZgkHBoPhWG0eGMZhD8zgwYtyz+ed+k1sbNjKNdOuxGaOC8hxRSS8BW34sUg4GdvkfbJbB4wNQda8WhERkcCwWc+uNp9MnjOHaUlF9A33sblpW8COKyLhTU2txISxvWr7Bia3yfu0pCIAKtTUioiIBIS/NgewqQW4OHdswaj1J6zwLSLRS02txIT4s7xTO/1oU3uks4oR30igwhIREYlZ8XGj04F6A9zULsyYj9PqoL67kfKOIwE9toiEJzW1EhNs1qOFs29yhTPFlkxKXDK9Q300dDcFMjQREZGYdLZTg07FbDSzMmc5AG/WaXsfkVigplZiwtid2t6ByRfO6Ufn1R7qOByIkERERGJavH/48eSmBp3OBTnLMRqM7GjZQ2d/V8CPLyLhRU2txITFM1185MJiCjOdkz7G3NRZALxS+Srewe5AhSYiIhKTlswarc15GY6AHzvFlsw56XMZ8Y3wdv3GgB9fRMKLmlqJCYtmurh6ZTEFZ9HULstaRElyMZ4BLy8eWBPA6ERERGLPklkZXL2ymPwgNLUwur0PwDt1GxgeCfzdYBEJH2pqRcbJaDBy85wbsJqsbG3eydamHaEOSURERE5hZsp0suIz6BzwsNNdFupwRCSI1NRKTGj39LPtYAvldZ1ndZx0exrXlXwIgNUH1miejoiIyCS1dfWx7WALFWdZm0/FYDBwYd7Y9j5aMEokmqmplZhwsKaDJ17azV8315z1sS7IWc6c1Jl0D/Xw6/2/0x54IiIik3CgerQ2/31rbdA+Y3nWEuJMVg51HKbe2xi0zxGR0FJTKzEhkBu8GwwGbpr9UexmG3ta97GhYctZH1NERCTWjNXmQO9Te8JnmG2cm7UEgDfr1gftc0QktNTUSkyID3DhTLElc8PMawH47aGXae1tD8hxRUREYoU97uge8kFsagEuyh0dgrypcSu9Q31B/SwRCQ01tRITbGOFcyBwqx8uy1zEAtc8+ob7eW7/bxjxjQTs2CIiItHu2J3a4K5MnOPIYkbyNPqHB3irVndrLLKGAAAgAElEQVRrRaKRmlqJCYG+Uwujw5A/Nus6HJYEDraX81bdhoAdW0REJNrZpmD48Zj35V8AwB8O/4W/Vr6m9TBEooyaWokJNuto4ewJcOF0Wh18bNZ1AKwp/xPNPS0BPb6IiEi0CsYF51M5J72U60o+hAEDfzj8F369/3fau1YkiqiplZhgizNhAPoHhhkZCezV2YUZ81mWuYiBkUGe3feihiGLiIiMg806NjVoKOh3Tg0GA+8vuIhPz/t/WIxm1jVs4ke7ntEcW5EooaZWYoLRYOCxO1fy5N0XYTAE/vg3zPwwSdZEDndW8Wr1m4H/ABERkShjNhn53p0refILF0/ZZy7MmM+/LfosDksC+9oO8v2tT9Le1zFlny8iwaGmVmJGijMOm9WMIQhdbbwlnpvmXA/A/x1eq73wRERExiHFGUec1RSU2nwqxUkFfHHpXWTGZ1Df3ch/bnmcak/w9soVkeBTUysSIKVps1iZs5wh3zC/2vu/mqsjIiISptLtady75A5mJE+jc8DDD7Y9xR73vlCHJSKTpKZWYsZvXivn0ee3UdXoCdpnXFdyFWm2FGq89Ty77zeaqyMiInIaq/9xiEee30Zts3fKPzveEs+dCz/NsszFDAwP8NSuX/Bm7bopj0NEzp6aWokZNS1eDtR00OHtD9pn2Mw2bpm7CrPBxOambXxn4/coa90ftM8TERGJZNVNXg7WdNDRHbzafDoWo5lPzL2Rfyq6DB8+Vh9cw+8O/VGLPopEGDW1EjP8WwcMBHfrgJLkYr687N8odObT0d/Jkzt/zq/2rqZ7sCeonysiIhJpxmpzX3/opuwYDAY+NO0Kbp5zA0aDkX/UvOUfjqzmViQyqKmVmGH374cX/MKZ48jiniV38JGSq7AYzWxs3Mp3Nn6PnS17gv7ZIiIikWKsNgd6H/nJOC97KXct+DQJ5ngOd1byo13P8NCmH7C+fjODI6GPT0ROTU2txAy7deo2eQcwGU1cVnAxXzn3C0xPKqJrwMNPdv+Kn+95Hs/A1M8dEhERCTe2uKN71YZBUwswK7WEb57/Za6d/kGSrIk0djfx3P7f8MC6R/hb1ev0DvWGOkQROQlzqAMQmSr2EBXOzHgX/774s7xZu54/VPyZrc07OdBezvUzP8ySjAVTuo2BiIhIOPFPDQqTphbAbrZzeeElvC//ArY07eDV6jep725kTcWfeaXyVVbmLOd9+ReQYksOdagicpSaWokZthAWTqPByCX5K5mXPpvn9/+Og+3lPFP2a7Y27eT6mdeQakuZ8phERERCzWaduqlBE2U2mjkveynLs5awt+0Af696g4MdFbxa8yav1b7N0syFXJR7PoWJeRgNGvwoEkpqaiVm5KQnsHR2BgWZzpDFkG5P4/MLb2dd/SZeKv8/drnL2O3ey5y0mazIXsb89LlYjPpnKSIisSHXNVqb8zMcoQ7llAwGA6VpsylNm011Vy1/r36Dbc272NS4jU2N23BaHcxLm8P89DnMSpmBzRwX6pBFYo5+e5aYUVqUSmlRaqjDwGAwsDJ3OXPTZvH78j+xo2UPe1sPsLf1AAnmeJZmLWJF9jLynTmhDlVERCSo5k9LY/60tFCHMW4FiXncOu8mrun9J16vfZsdzXto7+9gfcNm1jdsxmwwMSNlOvPT5zIvbQ5pdo3EEpkKampFQiTFlsyt827CO9DN5qbtrG/YTJ23gTdq3+GN2nfIc+SwInsZS7MW4rAkhDpcEREROSrdnspHZ1zDP5dcTX13I7vd+9jj3kdlVzX72g6yr+0gL7KGnIQs5qXPYSWLSBpJ12gskSDRvyyJGUPDI7R19TE84iM7LXyaRIc1gfflX8AleSup8daxoWELmxu3U+ut5zeH/sDvy/+P+a5Szs1cxOzUGVhN1lCHLCIiEhDhWpvHy2AwkOvIJteRzZVFl+IZ8FLWup897n3saztIfXcj9d2N/LXqNcxGM0WJ+ZQkT6MkuZjixEINVRYJEDW1EjNaO/v4yk82kJFs55HPrgh1OO9hMBgocOZR4MzjI9OvYpd7LxsatrCv7SDbm3exvXkXFqOZWSkzmJ8+h3npc0iOSwp12CIiIpPW3N7L1366kazUeB7+zHmhDuesOa0OzsteynnZSxkaGaK84wh73Puo8ByhurOO8o4jlHccAUYXkcx35FKSXExJcjHTk4tJsMSHOAORyKSmVmKGLYw2eD8Ti8nCkswFLMlcQHtfB5sat7GzpYwqTw17Wvexp3UfHIACZy7z0ucyP30O+Y5cbQ8kIiIRxR6GW/oEitloZnbqDGanzsDlclJZ30RFxxHKO0cb2xpPHVWeGqo8Nbxa8yYA2QmZFCcWUJRYQFFSAdkJmVpZWWQc1NRKzIg/uk9t30BkFc4UWzIfKLqUDxRdSmd/F3vc+9jdupf9beVUe+qo9tTx5yN/IzkuiXlps1k5uIQMQxY2sy3UoYuIiJyWfw/5CKvNk5FgieccVynnuEoB6Bvq50hX1Wij23GEI13VNHQ30dDdxLqGzQBYTVYKnXmjTW5iPkVJBRqlJXISamolZljMJswmA0PDPgaHhrGYTaEOacKS4hJZmbuclbnLGRge5ED7If/iFB39nbxdv5G36zdiNBgpTizwXyEudOZjMkZeviIiEt3iLCYMBhgYHGFoeASzKXbuStrMccxJncmc1JkADI4MUeOpo7KrmsrOaiq7amjta+NQx2EOdRz2vy85LomixHz/lKX8xFwtKCkxT02txBSb1Yy3d5Ce/mGSIrCpPZ7VZGF++lzmp8/F5/NR46ljt3svhzwVlLdWUtE5+t+fjvwNm8nGjJRpo01uygwy410aqiwiIiFnMBiwW8309A/RNzCMwx47Te27WYxmpiUVMi2pEPJHv+YZ8I42uV01/ka3o7+THS2d7GjZ439vmi2FfGcehUeb3AJnnubnSkwJalP78MMPs3PnTgwGA/fffz/nnHOO/7mGhgbuvvtuBgcHmTt3Lt/61reCGYoIAPFxo01tX/8QSQnRs4qwwWCgIDGPgsQ8XC4nVfXNHOqoYH/bIfa3H6K5x81u9152u/cCo1d5Z6WUMDNlOjOSp5FmD/3+vSIyNVSbJdzY40z09A/R2z+Ew24JdThhxWl1+C9gA4z4RmjuaaGyq4ZqTx01nlpqPPW09rXT2tfOjpbd/vem2VIpcOaS58wl35lDniOXpDhnqFIRCaqgNbWbNm2iqqqK1atXU1FRwf3338/q1av9zz/yyCPceuutXH755Xzzm9+kvr6enJycYIUjAhxbkCISFos6G/EWOwtc81jgmgdAW187+9vKOdB+iP1th+jo72Rj41Y2Nm4FINWWwozkacwYa3JtKbqTKxKFVJslHI3W5v6oXCwq0IwGI1kJmWQlZHJe9lIAhkeGaeppodpTO/pfVx213jpa+9po7Wtj+3GNbqLVSZ4zh3xHrv/PdHuqar5EvKA1tevXr+eyyy4DYPr06XR2duL1enE4HIyMjLB161a+//3vA/DAAw8EKwyRE9z8gVmMjPjIToutITmpthTOz1nG+TnLGPGNUO9t5EB7OYc6DlPecYS2vvYTmtyUuGRmpkynJHkaM5KnqeCJRAnVZglHn/in2fh8kJkSW7U5UExGEzmOLHIcWSc0uo09zVR76qj11FHjqafOW0/XgIe9rQfY23rA/36byUaeM5tcRw55R/fczU7IwmrSXXOJHEFrat1uN6Wlpf7HqamptLS04HA4aGtrIyEhgf/4j/+grKyMpUuXcs8995zxmC5X9AyZUC6hcbpYIymPMzlTLpkZSSxiFgAjIyNUd9ZR1nyQspZD7Gs5RHt/x4lNri2JWa7pzE6fzuz0EgqTc6ds4alYOi+RIlryiEWqzaenXEJDtTk4skhmITP9j0d8IzR3t1LZXsOR9hoqO0b/7OjrOmH/XBid1pTjyKQwOZfC5Lyj/+WSak8OSS7BFC25REsekzVlC0X5fL4T/t7U1MQtt9xCbm4un/nMZ3j99de55JJLTnuMlhZPkKOcGi6XU7mEmWjJAyaXSwLJnJt6LuemnsvIzNE7uWOrLZZ3HKa9r5MNNdvYULMNgDiTleLEQqYlFzE9qYiixAJs5riwyCVcRUsu0ZIH6BcAUG0+XrR9b0dDLtGSB4RPLiZsTLfNYHr2DMge/Vpnv4dabx113gbqvA3Uehto7mmhztNInaeRdTVb/e9PsMRTlJJHujWdnITRu8PZCZnYzfYQZXR2wuW8nK1oyQMmX5uD1tRmZGTgdrv9j5ubm3G5XACkpKSQk5NDQUEBACtWrODQoUNnLJwiZ2tnuZs9R9pYUJLGvOK0UIcTlowGI3nOHPKcObwv/4LRX3R7WqjoPEJFx+iKyu7eVva3jy5C5X+PI5vipEKKEgsoTizUkGWRMKTaLOFo+6EW9h5pZ9HMdOYWaeHCqZYU5yQpbjalabP9XxscHqShp4k6z1ijW0+dt4HuwR7Kmg8CB084Rkpc8ugQaH+jm0VWvAuLhjDLFAlaU7ty5Uoef/xxVq1aRVlZGRkZGTgcjtEPNZvJz8+nsrKSoqIiysrKuOqqq4IViohfRX0nr26txRlvUVM7TgaDgayEDLISMliZsxwYvap7uLPS3+jWeuup9tRR7anjDdYB4LAkjDa4SQUUJRZQmJiP3WwLZSoiMU+1WcJReW0nr26rJdlpVVMbJiwmi38f3DE+n4+O/k66zZ3sqztCfXcjDd5GGnqaae/voL2/g7LW/f7XGzDgsqeRlZBJdkImWQkZZCdkkhnvwmqKnh0oJDwEraldvHgxpaWlrFq1CoPBwAMPPMBLL72E0+nk8ssv5/777+e+++7D5/Mxc+ZMLr300mCFIuI3tvqxVlg8O0lxThZlzGdRxnwA+ob6qeqqobKrmiNd1RzprMI72M2e1n3sad0HjBa3rIQMio82uIWJ+eQkZE3Z3FwRUW2W8HSsNg+HOBI5HYPBQIotmZmufPLMhf6vD48M4+5tpb67iXpvw+if3Q209LTS3OumudfNLnfZseNgIM2WckKzmxmfQWa8i3hLZA5jltAL6pzae++994THs2cfG9ZQWFjICy+8EMyPF3kPu1WFMxhs5jhmpZYwK7UEGL2a29rXTmVn1WiT21VNraeehu4mGrqbWNewGQCz0UyeI4fCxNGrwYWJ+WTGuzAajKFMRySqqTZLuPE3tQO64ByJTEYTmQkZZCZk+C92w+gQ5uZeNw3dTTR2N9HQ3UxjdxPNvW7cfW24+9r8F77HJFqdZMa7yEzIIOtoo5sZn0GKLUm/G8hpTdlCUSLhQHdqp4bBYCDdnkq6PZWlWYuA0eJW462nsrOKKk8t1V21NPe6qeyqprKr2v/eOJOVfGcuhc585vWWkORLw2VPUzETEYlS9rjRETuqzdHFYrKQe3SLoOMNjQzR0ts6eqHb20hTTwuNPc0097TQNeCha8DDoY7DJx7LaCEz3kVGfPrRP4/9PVIXqZLAUlMrMUWFM3QsJgvTkgqZlnRsyFLPYM/oXNyuWqo8NVR11dLe3+HfWuDVmjcBsJniRjeJd+aS78ilIDFPd3RFRKKEfxRVn2pzLDAbzWQfHXpMxjn+r4/4Rmjv66Cxp4Wmnmaaupv9Da9nwEutt55ab/17jue0ON7T7GbEu0i3p2ExqtWJFTrTElM0xCm8xFvimZ06g9mpM/xf6xrwjDa5XTU0DTRR0VpNR3/ne/bQsxot5DpGG90CZy75zlyyEzI1R1dEJMIcq82aGhTLjAYjafZU0uyplKbNOuG5nsFemnqaae5x09zTQlPv6J/NPW48g148nV4qOitPeI8BA6m2ZFz29OOa3XRc9nTSbCn6fSHKqKmVmOKwW8hOiyctUavwhqtEq5N56XOYlz7Hv+9a14CHGk8dNUdXWK7x1NHW186RriqOdFX532s2mMhxZI/e0T3a7OYkZGlLARGRMOaIH6vNgd/vXKJDvMVOcVIhxceN9oLRu7ud/V00HW1wRxveFlp63LT2tfv/G9uCcIzRYCTdnkpecjbJpmR/s5sZ7yIpLlEjwSKQmlqJKdlpCTx0+3mhDkMmKNHqpDTtxD30vAPd1HjrqOmqG/3TU0dLbyvVnlqqPbX+1xkNRrITMsl35JKfODpXN8+RrUZXRCRM5Lkcqs0yKUaDkRRbMim25BNGfcHo3N3Wvnaae0ab3KZeNy09bpp73LT3dxxtgt3vOabFaMFlT/MPY86wp5OZMLpgVYIlfqpSkwlSUysiEclhTWBO6kzmpM70f613qJdaT/3RO7r11HjraOpups47unn8hsYtAJgMJnIcWaNbCznzKUzMIzshU1dmRUREooTZaD66erLrPc8NDA/S0uum39JNeUPN6NZDxw1nru9upL678T3vc1gS/CsyZya4jq7QnEGaPUW/Q4SYmlqJScMjIxgNBgwGQ6hDkQCym+3MSJnOjJTp/q/1Dw9Q522g2lNLTVcdlZ4amrqb/cOZ32YDAFaTlXxHLoWJeRQlFlCSPI2kOGeoUhERiTmqzTJVrEdXZna5nEyLKznhud6hXv9d3LFmt+nof97Bbryd3e+Zv2s2mMiId5GVkEGuI5uchCxyHdmk2lL0/TxF1NRKzPnC42/T2T3AE/9+EfE2/ROIdnEm63tWXe4b6qPGU0dlVw1VntFFqdr62qnoPEJF57HFqDLjM5iRXDzaKCdPIykuMRQpiIhEvc//91t4ewf50d0XE2fVAj4SOnazfXQkV2L+CV/3+Xx09Hf6V2Ru6j66SnNPCx39nf67u9uad/nfYzPZyHFkkuPIJjcha/RPR5a2IQoC/UYvMcdoHL1i1jcwpKY2RtnMtvfc0fUMeKk62uQe7qjkcGfl0WLVzNv1GwHIjHcxI3mamlwRkQAbq829A0NqaiUsGQyGU87f7Rvqo6mnhYbuJuq8DdR7G6nrbsAz4OVwZxWHO6tOeH26LZVpyUVHL7oXaQpUAOg3eok59jgz7Z5+evqHSA11MBI2nFaHf9VlgOGRYao8tRxqr+BQx2EqOiv9w4+Ob3IXueazLGsRWQmZoQxfRCSi2a0murpH95FPdmgVZIksNrPtpHd3uwY81Hsbqfc2UOdtpL67gYbuJtx9bbgb29jUuG30/SYbxUkFTE8qojipkKLEAmxm/TuYCDW1EnPsR68A9/Zrr1o5NZPR5B+2/AEuZXhkmGpPLYfaD3Owo8Lf5L5S9Q9eqfoHeY4clmUtYmnmQpLjkkIdvohIRBnbq7ZHtVmiSKLVSWKq84Q7u8Mjw9R3Nx69g1tJRUcl7f0d7Gs7yL62g8DoHrt5jmxmpExnccY5FCUWaG7uGYyrqS0rK6O0tDTYsYhMCf8m7/3a5F3Gz2Q0+ffIu4L3MTwyTHnHEbY0bWd7y25qvfXUltezpvzPzEiexrKsRSx0zSfeonkzEhyqzRJNxmpzn2qzRDmT0US+M5d8Zy4X550PQHtfh7/JPdxZRa23npqj//2j5i3SbCksyVzI0syF5CRkqcE9iXE1tT/4wQ+orKxkwYIFXHjhhaxcuRKX673LY4tEgmNNra4Gy+SZjCZmpZYwK7WEG2ZeS1nrfjY37WBP6z4OdlRwsKOC1Qd+T2n6HJZlLuLStOWhDlmijGqzRBPVZollKbZkltiSWZK5ABjduaGqq5rd7n1sbdpJa187f616jb9WvUZWfAZLMxeyJHMBGSfZrihWjaup/elPf8rQ0BC7d+9m48aNfPnLX6atrY01a9YEOz6RgLPHHR1+PKDCKYFhMVlYmDGfhRnz6RnsZUfLHrY0bedgewU7W/aws2UPa2te5aPTP0xJcnGow5Uoodos0URTg0SOiTNZmZlSwsyUEj5SchUVHUfY0rSD7S27aexp5v+O/JX/O/JXCpy5LMlcyBUJK4n1WaXjyn5wcJCdO3eyadMmtm7dSn9/P8uX666DRKbz5mZRlJXIzLzkUIciUSjeYuf8nGWcn7OMjv5OtjXt5I3addR01vODbT/ivKylXFvyQZxWR6hDlQin2izR5Px5WUzPS6IkT2sSiBzPaDD6d2y4Yea17G8/xJamHexqKaPaU0e1p44/VPyFS/JWclXx5djMtlCHHBIGn8/nO9OLFi5cyLJly/j4xz/OueeeS0JCwlTE9h4tLZ6QfG6guVxO5RJmoiUPUC7haGB4kHfc77Bm71qGfMPEm+18ePo/cX7OuRG3hH+0nBMYzSWSqTYHVrR9b0dDLtGSByiXcBXJuQwMD1LWup8tTdvZ2VKGDx9J1kQ+OvMaFrnmR+y828nWZtODDz744JleNGvWLHp7e/nLX/7CP/7xDxoaGjCZTGRlZU3qQyerp2dgSj8vWBIS4pRLmImWPEC5hCOT0cS5xfOZ7ZhNU3cLDT1N7Gndx762gxQ4cyNqv9toOScwmkskU20OrGj73o6GXKIlD1Au4SqSczEZTWQnZLIkcwEXlCyhvKWKpt4Wtjfv4nBnFUVJBTgsobnYeTYmW5vHdaf2eDt27OCpp57izTffZO/evZP60MmK1Csp7xbJV4XeLRJzaenoZW9lGylOG+dMTwMiM49TUS7haSwXn8/HtuZd/O7Qy3QOeDBg4KK887l62hXYzeG/UnK0nZNoodp89qLtezvScmk+WpvTE23Mm6baHM6US/hxuZw0NXfyTv0mXq74Cz1DvZgNJi4rvIQPFF6K1WQJdYjjNtnaPK45tc899xybN29m//79lJSUcOGFF/LlL395Uh8oEmrVTV5++coBFs1I9ze1IlPFYDCwJHMBc9Nm8acjf+X1mnd4o/Ydtjfv4p9LPsSSzIURO2RIppZqs0STyoYufvXKAZbOcvmbWhEZP6PByIW557HQNY81FX9mQ8MWXql8lc2N27l+5jXMT58b6hCD6rSTubxeL7/4xS9ob2/nlltu4eabb6ampoaNGzfidEbPFW6JLf7Vj7XCooSQ3WzjozOu4cvL/o3ixEK6Bjw8s/cFfrVvNSO+kVCHJ2FMtVmikX9LnwHtUytyNpxWBzfPuYG7F99BTkIWrX1tPLXrF/x41y9p7W0PdXhBc9qm9hvf+Aatra187nOfIzU1lf/5n//hvvvuY+XKlTz00ENTFaNIQB3bC0+FU0Iv35nD3Uv+lZtmf5Q4k5VNjdv43aE/MsGZIRJDVJslGmmfWpHAmp5cxH3L/o1/LvkQcSYru9xlfHvjY+xx7wt1aEFx2qa2pqaGe+65B4C1a9dy5ZVXcv7553PjjTfidrunJECRQFPhlHBjNBg5P+dcPjP/E5gNJl6vfYe1Va+FOiwJU6rNEo20T61I4JmMJi4tuIhvnPdFFmWcw+DIID/d8xyHOytDHVrAnbapjY+P9/9906ZNnHfeef7HmvMlkerYECcVTgkvs1Nn8InSj2HAwB8Pv8I7dRtDHZKEIdVmiUa64CwSPMlxSdxWehMrspcxODLIj3Y+Q723MdRhBdRpm9rh4WFaW1uprq5m+/btrFy5EoDu7m56e3unJECRQIvXnFoJY4szzuGGmdcC8MKBl9jRvDvEEUm4UW2WaKSpQSLBZTAY+Nis65ifPpeeoV5+uPNntPd1hDqsgDltU3v77bfzwQ9+kKuvvpo77riDpKQk+vr6+PjHP8611147VTGKBJTZZMRkNGA0GBgc0oI8En4uylvBVcWX48PHM2W/5mB7eahDkjCi2izRKM5qwmAAHz5GRrSmgEgwmIwmbi29ielJRXT0d/LEjp/iHewOdVgBccZ9agcHB+nv78fhcPi/9vbbb3PBBRcEPbh3i4Z9pCB69sSCyM1laHgEs+nYNZ1IzeNklEt4mmguPp+PFw/+gTfr1mEzxfHviz9LvjM3iBGOT7Sdk0il2hx40fa9HYm5qDZHBuUSfiaaR89gD9/f9iMaupsoTizgc4s+Q5zJGsQIx2+ytfm0d2oBLBbLCUUTCEnRFAmk44umSDgyGAxcP/MaFmecQ99wPz/c8TOae7QIkIxSbZZopNosMjXiLfHctfDTpMQlc6Srmp/veY7hkcge+q+fHiIiYcpoMHLL3FXMTpmBZ9DLEzt+Smd/V6jDEhERkQiXHJfEXQs/TYI5nj2t+3l+/28jejtBNbUSk5796wG++vQGyms7Qx2KyGlZjGZun38zBc48Wvva+OHOn9EzqMWARCT6/PKV/Xz16Q0crtfFO5GpkJWQwb8uuBWr0cLGxq38oeIvoQ5p0tTUSkxq7+qnobWHrp6BUIcickY2s407FtxKRnw6dd4Gntr1CwaGB0MdlohIQLV7VJtFplpxUgGfnn8LRoORv1W/zj+q3wx1SJOiplZikl3b+kiEcVod3LXgdpKsiVR0HuHX+38X6pBERALKZlVtFgmF0rRZ/L/Z1wPwu/L/Y1PjthBHNHFqaiUmaZN3iURp9hTuWvhpjAYjW5q209kf+Ss2ioiMGavNfarNIlNuefYSPlJyFQDP7nuRQ+0VIY5oYtTUSkxSUyuRKseRRWnabHz42Nq8I9ThiIgEzFht7lFtFgmJywou5n35FzDiG2Ft1WuhDmdC1NRKTPI3tQORvXy5xKZzsxYDsDkChweJiJyK/ejw4z7VZpGQubLo/RgNRg60l+MZ8IY6nHFTUysxSXdqJZLNS5uDzWSj2lNHY3dTqMMREQkI3akVCT2HJYG5qTMZ8Y2wrXlXqMMZt6A2tQ8//DA33ngjq1atYteuk/9P+d73vsfNN98czDBE3iM/w8Gli3OZU5gS6lBEJsxqsrA4Yz4Amxq3hzgaiTSqzRKuCjKdXLo4l9kFqs0iobQkcyEAW5oiZ5qTOVgH3rRpE1VVVaxevZqKigruv/9+Vq9efcJrysvL2bx5MxaLJVhhiJxUSW4SJblJoQ5DZNKWZS1mXcNmNjdt50PTrsBo0MAbOTPVZglnM/OTmZmfHOowRGLeOemlWIwWDndW0trbTpo9/C80Be23oPXr13PZZfbsGDAAABsCSURBVJcBMH36dDo7O/F6TxyX/cgjj/CFL3whWCGIiEStkuRikuOSaOtr53BnVajDkQih2iwiImdiM8dxTvpcgIhZlDJoTa3b7SYl5VhXn5qaSktLi//xSy+9xLnnnktubm6wQhA5pcGhYSrqOjlU2xHqUEQmxWgwsixzEUBE7icnoaHaLOFsYHCY8rpOyms7Qx2KSMxbGmFDkIM2/PjdfD6f/+8dHR289NJLPPPMMzQ1jX+RE5fLGYzQQkK5hFZzWw8PPbuV9GQ7z3z9CiAy8zgV5RKeAp3LBywX8Lfq19nRsos7zr8Ji2lqhotG0zmJdarNJ1IuoVXv9vLws1vJSovn6fsvByIzj1NRLuEpWnIJdB4XpS7huf0vUudtoM/qIT8pJ6DHD7SgNbUZGRm43W7/4+bmZlwuFwAbNmygra2Nm266iYGBAaqrq3n44Ye5//77T3vMlhZPsMKdUi6XU7mEWG/fIADdvQO0tHgiNo+TUS7hKRi52Ekkz5FDrbee1w9sYaFrXkCPfzLRdk5ijWrzqUXb93Yk5tLbMwCAt2dQtTmMKZfwE6w8FqTPZ13DJv627x2unn5lwI9/MpOtzUEbfrxy5UrWrl0LQFlZGRkZGTgcDgCuvPJK/vznP/Piiy/yxBNPUFpaesaiKRJINuvo9Zy+/mFGjrtTIRJplmWNDkHWnrUyHqrNEs7s1mPb7flUm0VC7vghyOH+bzJod2oXL15MaWkpq1atwmAw8MADD/DSSy/hdDq5/PLLg/WxIuNiNBqIs5roHximX5u8SwRbmrmQNeV/Zo97Hz2DPcRb4kMdkoQx1WYJZxazEbPJyNDwCINDI6EORyTmzUiZRpLVibuvjcquGoqTCkId0ikFdU7tvffee8Lj2bNnv+c1eXl5PPvss8EMQ+Sk4uPM9A8M06tN3iWCJcclMSulhP3th9jWvIsLcs8LdUgS5lSbJZzFx5no6hlRbRYJA0aDkcWZC3it5m22Nu0I66ZWGxtKzLJZTQAqnBLxxoYgb2rcHuJIRETOji3u6BBkjaISCQtjOy1sbd7JiC98R1CoqZWYFT9WOPtVOCWyLXTNw2K0UNF5hNbetlCHIyIyafa4Y/NqRST0Cpx5pNvT6BrwcLC9ItThnNKUbekjEm4+/aG5YIBUpy3UoYicFZvZxgJXKVuadrC5aQdXFl0a6pBERCbls9eUqjaLhBGDwcCyzIX8pfJVtjTtYHbqjFCHdFK6UysxKzM1nsyUeCxm/TOQyDc2PGhT47awX6FQRORUVJtFws/YKsg7WnYzOBKeoyj0E0NEJArMSZ2Jw5JAU08zNZ66UIcjIiIiUSIrIZM8Rw69Q33sbd0f6nBOSk2txKzN+5t56g972LK/OdShiJw1k9HEkqNXUjc1ac9aEYlMG/c28dQf9rDtYEuoQxGR4xy/Z204UlMrMave3c2mfc3UNHtDHYpIQCzPWgyMFpzhES2AJiKRp87tZdO+ZmpVm0XCypLMBQDsdu+lb6gvxNG8l5paiVlaYVGiTYEzj4z4dDwDXg60l4c6HBGRCRurzT2qzSJhJdWWwvSkIgZHhtjl3hvqcN5DTa3ELLv2qZUoYzAYODdz9G6t9qwVkUhkt+qCs0i4Wnp0UcrNTeH3O4aaWolZdm3wLlFoWdZowdnZspu+of4QRyMiMjGqzSLha1HGfIwGI/vbDuEZCK8pAmpqJWbZbboaLNEn3Z7GtKRCBkYG2eUuC3U4IiIToqlBIuHLaXUwO3UGI74RtjfvDnU4J1BTKzFrbIiT5u1ItDk3a2wIslZBFpHIYo/T1CCRcLY0IzxXQVZTKzErMd7CzLwkirOcoQ5FJKAWZZyDyWBif9shOvs9oQ5HRGTcEhOszMxPpjBTtVkkHC1wlWIxmqnoPEJbX3uow/FTUysxKz3Zzn3/bwm3XDk71KGIBJTDkkBp2mx8+NjaHF5XUkVETiczJZ77blrMzR+YFepQROQkbGYb89LnArC1aWeIozlGTa2ISBQaWzBKQ5BFREQkkJZlht8QZDW1EtP6B4fp8PYzMuILdSgiATU/bQ52s50aTx1VXTWhDkdEZNz6B0Zrs8+n2iwSjuamzcZutlHrraexuynU4QBqaiXGffHJdf+/vbsPjqq+9zj+2ezmkWxiFpIQSPAiKirVWhHbEkwUIdaO9w5lFIKl1RlEvdDbdiidKa2UmU5hqpfrWFMdlGvv9HbuVFqkaNvppdXyVAiPRZFA5UEIeSLZJBCyeWQ35/6RZCU0m3CVzTm/3ffrH909Z875fv3G/eZ7zsn+tOxnu9Ta3m13KMA1lehO1PS8aZKk7dW7bY4GAK7espf/qmU/26W2Tr4sCnCixASP7sy+XZJz7tYy1CKupST1fstiW+clmyMBrr2i/OlyyaWD9e/pYjdfGAXADCn9qxN00JsBp5qW2/tnTrtq9+lSj/0XoBhqEdfS+tbDa+dqMGLQmFSfbh9zm4JWSLtq9todDgBclXBvZlkfwLFuzpqk8el5utjdqn3nDtodDkMt4ltqeKjlajBiU3H+dEnSzppyhXpCNkcDAMPr781t3KkFHMvlcmnWhGJJ0jtnt6vH6rE1HoZaxLVU7tQixk3OulFjR+WqpbtVh/wf2B0OAAwrJbn3T4M6uFMLONrUnM/Kl5KlhvZGfdB41NZYGGoR11L7Gid3ahGrXC6X7ssvlCRtq9plczQAMLw07tQCRnAnuDWz4F5J0p8rt9n6jeUMtYhrKeHGydVgxK57xt6lVE+qTl+sZHkfAI4X/qIoLjgDjjd93D0a5UnT6YtndarljG1xMNQirhXdMU5Lv3K7pt2Wa3coQNQku5NY3geAMYrv7O3NU2+hNwNOl+xOUlHf93f8uXKrbXEw1CKuXT/Wq6mTszV29Ci7QwGiiuV9AJhiYl6Gpk7OVo4vze5QAFyF4vzpSkxI1JGmv6s2cM6WGBhqASAOsLwPAACIBm9Sur7Y90TYO2e32xIDQy3i2rnmdv1u9xlt+1u13aEAUcfyPgBMUNvYpt/tPqMdh+jNgCkemHCvXHJpf/0hne+8MOLnZ6hFXDvX3K7f7vhI2xlqEQdY3geACeqaenvzjkM1docC4CqNSR2tu3LuUI/Vo79U7Rzx8zPUIq6lJvUu6cOyAYgHLO8DwAT9y+218e3HgFFmX3+fJGlX7V61X2of0XMz1CKupfYt6cMC74gXLO8DwOn6e3N7J70ZMEmBd7xuybpJXaFu7ajZM6LnZqhFXAsv8M7VYMQJlvcB4HRpyaxTC5iq/27ttqq/qjs0cv8PM9QirqX0N04eP0YcYXkfAE7W35vbOrhTC5hmctaNKvCOV+ulgPaeOzhi52WoRVzr/7ud9q6gLMuyORpgZIxJ9ekzY25leR8AjpTW15s7ui7RmwHDuFwuzZ5QLEl69+x29Vg9I3JehlrENXdCgrK8ycq+LlXdwZH5nw5wgv4vjNpZU65gD3dDADhHosetjFFJ8mWkKBhiqAVMc2f27RqT4pO/o0nv+Y+MyDkZahH3/mNpoV5/tkTJiW67QwFGzOXL+4xUwwGAq/Xiv83Q68+WKNHDr6qAadwJbj0woUiS9E7l9hF54oJPCgCIQ73L+0yXxPI+AADg2vpC3t1KTxylytYqnbhwKurni+pQu2bNGs2fP1+lpaU6fPjwgG179uzRvHnzVFpaqhUrVqinh0c/AWAk3TN2qlI9KSzvE2fozQCAaEtyJ4X/1OnPldujfr6oDbX79u1TZWWlNmzYoNWrV2v16tUDtv/whz/USy+9pDfeeENtbW3auXNntEIBhvSfvz+qx1b+URWnm+0OBRhRvcv73CNJ+kvVTnUGu0bsCx1gD3ozTPHa2xV6bOUfdazyvN2hAPiE7s3/opISEnW0+UNVt9ZG9VyeaB24vLxcs2bNkiRNmjRJLS0tCgQCSk9PlyRt2rQp/O8+n0/nz/OhBXt0Xwqptb1bv935kbYdqgm//9kbx2jGHXmSpGp/QG/tPB3xGI/NvllZ3mRJ0p/2ndWJ6pZB9xs3ZpS+UnRD+Lzrf3c04jFnTyvQzQXXSZIOnfBr9wfnBt0v0ZOgp/5lSvj1q5sOq84fGHRfU3P6nz8d14VAl9E5pY9K1uMP3uy4nLoTfJJXOlD/ng7UvydJcrvcSnInKSkhUVZPgrq6pATLI5flljvBI6vHkuSSy+XS9bleuVyS5FJ1Q0DBUI9kufrO6gqfPzM9SWMyUyRJHV1B1Ta2R8ypIDddSX1/R1ff3K5AhGU9UpLdGj9mlCSpp8fS6brIyxPlZKXKm5YoSboQ6FJTS5de/9qKiPvHKnozTNHV15s37Til60Ylh9+/6+ZsffEzYyVJVQ0Bvf3XyJ+PCx+crMxRSZKk/917VqdqBv98HJ89SnPu7f187OwO6vXfH4t4zAfvmaAb8zMlSQc/bNCeivpB90tKdGvxP98Wfv3KxvfV0NQ26L6m5vTfWz5Ua1u30TlleJP1tdkf92ZTc0pK9qi7q7dXOqlOHitZ6R2T1Jz8d/373lfl6UlR7+8GLrmUoCxvikYlJ8rlSlCg/ZJaApf0XwtXRoxhKFEbahsbGzVlyse/lPp8Pvn9/nCz7P9nQ0ODdu3apW9961vDHjM72xudYG1ALs5x/bhMHfjQr49qLw54vyAvI5zbuZYuHTzuj3iMJ79ye3jf2vMdEfftCvWE92vvvDTkMe+bNiG8b1tFfcR9U5PdA2pw6HiDavyDN04Tc8rO9urvVedjIqflC6eGXzspJ0/BRHnGVEsJPXK5QwpZIXUEO9Shjt4drvwOtcteH2m67P0ERXz+p7Vbqr48tKSIKenolXNUhH1bLcl/lcdsbZN0+X/uIfaNZfTmoZGLc1w/LlOHTjTqVM3A3vxP4zPDudUM8ZknSc888lllj+698FXT1B5x36BlhY8ZaO8e8pizvnB9eN/Ww3UR9x2VmjigBn/7sEH1zYNfzDMxp+xsr45Vno+JnJY99nFvjpWcLv/Zszun2mO5Sr79hILuDgUTOgbs29Eh6fK3EiOeflguK0pfR7Vy5UoVFxeHrwgvWLBAa9as0cSJE8P7NDU1afHixVq2bJlmzJgx7DH9/shX4U2Sne0lFwfpvhRSzYVONV1xBTXXl6aCnN5f8C62d+v42QsRjzFlok+pfYvFf1R7Uc0XOwfdz5uWqMkTsiRJwVCP3jvRGPGYE/MyNLrvzta55nZVNwx+99Wd4NLnbs4Ov65q6lB9hJqYllP/z9eRj5rU2R0yOqesrDRNyk0Pv3ZqTpZlqWBsmrzpbnWHulV7vlU1jS0KWkEFraCS0xLU3tZ7h9nlkm4Yl6H+JnL23EV1B3tk6R/bSmZ68hV3agcf6CWpICddSX3fRt57p/bSoPulJrs1bkzvf6dQT4/ODHuntneS7b1T26l/LXko4v6xit4cWSz0s36xkEvXpZBqz3eqqXngZ8VYX5ry+z4fW9q6daIq8ufjZ27wKSWp9/PxVE2LzrcO/nRMxqik8BM3w33m3zAuQ76M3s+yuqa2iBcn3W6XPnfTx735bFO7GiI8RWVaTv0/Xx981KSuCH3MlJx8WWm64bLebGpOGZmputjSOx1e+bPnhJw6Q+1q7wmox7Jkqff3BMuylDcmVelpHvVYlhpb2uVv6dBTs2dHjGEoURtqy8rKlJ2drdLSUknSAw88oLfeeit8FTgQCOjrX/+6vv3tb6uoqOiqjmn6B3S/WGg2/WIll1jJQyIXp4qVXGIlD8n8O1mfBL05slj72Y6FXGIlD4lcnCpWcomVPKRP3puj9kVRhYWF2rJliySpoqJCOTk54aYpST/5yU/0+OOPX3XTBAAAnw69GQAQi6L2N7V33XWXpkyZotLSUrlcLq1atUqbNm2S1+vVjBkztHnzZlVWVmrjxo2SpIcffljz58+PVjgAAMQ9ejMAIBZFbaiVpOXLlw94fcstt4T//ciRI9E8NQAAGAS9GQAQa6L2+DEAAAAAANHGUAsAAAAAMBZDLQAAAADAWAy1AAAAAABjMdQCAAAAAIzFUAsAAAAAMBZDLQAAAADAWAy1AAAAAABjMdQCAAAAAIzFUAsAAAAAMBZDLQAAAADAWAy1AAAAAABjMdQCAAAAAIzFUAsAAAAAMBZDLQAAAADAWAy1AAAAAABjMdQCAAAAAIzFUAsAAAAAMBZDLQAAAADAWAy1AAAAAABjMdQCAAAAAIzFUAsAAAAAMBZDLQAAAADAWAy1AAAAAABjMdQCAAAAAIzFUAsAAAAAMBZDLQAAAADAWAy1AAAAAABjMdQCAAAAAIzFUAsAAAAAMBZDLQAAAADAWAy1AAAAAABjMdQCAAAAAIzFUAsAAAAAMBZDLQAAAADAWFEdatesWaP58+ertLRUhw8fHrBt9+7deuSRRzR//ny9/PLL0QwDAAD0oTcDAGJN1Ibaffv2qbKyUhs2bNDq1au1evXqAdt//OMfq6ysTL/61a+0a9cunTx5MlqhAAAA0ZsBALEpakNteXm5Zs2aJUmaNGmSWlpaFAgEJElVVVXKzMxUXl6eEhISVFxcrPLy8miFAgAARG8GAMSmqA21jY2NysrKCr/2+Xzy+/2SJL/fL5/PN+g2AAAQHfRmAEAs8ozUiSzL+tTHyM72XoNInIFcnCdW8pDIxaliJZdYyQP05iuRi/PESh4SuThVrOQSK3l8UlG7U5uTk6PGxsbw64aGBmVnZw+6rb6+Xjk5OdEKBQAAiN4MAIhNURtqCwsLtWXLFklSRUWFcnJylJ6eLknKz89XIBBQdXW1gsGgtm7dqsLCwmiFAgAARG8GAMQml3Utnj2KYO3atTpw4IBcLpdWrVqlo0ePyuv1avbs2dq/f7/Wrl0rSSopKdGiRYuiFQYAAOhDbwYAxJqoDrUAAAAAAERT1B4/BgAAAAAg2hhqAQAAAADGGrElff4/1qxZo/fff18ul0vf//73dccdd4S37d69Wy+88ILcbreKioq0dOlSGyMd2lB5zJw5U2PHjpXb7ZbU+zdOubm5doU6rOPHj2vJkiV64okntHDhwgHbTKqJNHQuptXl+eef18GDBxUMBvX000+rpKQkvM2kugyVh0k16ejo0Pe+9z01NTWpq6tLS5Ys0f333x/eblJNhsvFpLpIUmdnpx5++GEtWbJEc+fODb9vUk3sRm92HnqzM9GbnYXe7My6SNe4N1sOs3fvXuupp56yLMuyTp48ac2bN2/A9oceesiqra21QqGQtWDBAuvEiRN2hDms4fK4//77rUAgYEdo/29tbW3WwoULrWeffdb65S9/+Q/bTamJZQ2fi0l1KS8vt5588knLsiyrubnZKi4uHrDdlLoMl4dJNfnDH/5gvfbaa5ZlWVZ1dbVVUlIyYLspNbGs4XMxqS6WZVkvvPCCNXfuXOvNN98c8L5JNbETvdl56M3ORG92Hnqzc13L3uy4x4/Ly8s1a9YsSdKkSZPU0tKiQCAgSaqqqlJmZqby8vKUkJCg4uJilZeX2xluREPlYZqkpCStX79+0PUKTaqJNHQuppk2bZp++tOfSpIyMjLU0dGhUCgkyay6DJWHab785S9r8eLFkqS6uroBV0dNqok0dC6mOXXqlE6ePKn77rtvwPum1cRO9GbnoTc7E73ZeejNznSte7PjHj9ubGzUlClTwq99Pp/8fr/S09Pl9/vl8/kGbKuqqrIjzGENlUe/VatWqaamRlOnTtV3vvMduVwuO0Idlsfjkccz+I+KSTWRhs6lnyl1cbvdSktLkyRt3LhRRUVF4cdNTKrLUHn0M6Um/UpLS3Xu3DmtW7cu/J5JNbncYLn0M6Uuzz33nFauXKnNmzcPeN/UmtiB3uw89GZn1oXe7Fz0Zme51r3ZcUPtlawYWXHoyjy++c1v6t5771VmZqaWLl2qLVu26Etf+pJN0aGfiXV55513tHHjRv385z+3O5RPJVIeJtbkjTfe0LFjx/Td735Xb7/9tmMbytWIlIspddm8ebPuvPNOFRQU2B1KTKE3YySZWBd6s/PQm50jGr3ZcY8f5+TkqLGxMfy6oaFB2dnZg26rr6937KMqQ+UhSXPmzNHo0aPl8XhUVFSk48eP2xHmp2ZSTa6GaXXZuXOn1q1bp/Xr18vr9YbfN60ukfKQzKrJkSNHVFdXJ0m69dZbFQqF1NzcLMm8mgyVi2ROXbZt26Z3331X8+bN029+8xu98sor2r17tyTzamInerNZTKrJ1TCtLvRmZ6E3O080erPjhtrCwkJt2bJFklRRUaGcnJzwY0H5+fkKBAKqrq5WMBjU1q1bVVhYaGe4EQ2VR2trqxYtWqTu7m5J0v79+3XTTTfZFuunYVJNhmNaXVpbW/X888/r1Vdf1XXXXTdgm0l1GSoP02py4MCB8NXsxsZGtbe3KysrS5JZNZGGzsWkurz44ot688039etf/1qPPvqolixZounTp0syryZ2ojebxaSaDMe0utCbnYfe7DzR6M0uy4HPEK1du1YHDhyQy+XSqlWrdPToUXm9Xs2ePVv79+/X2rVrJUklJSVatGiRzdFGNlQev/jFL7R582YlJyfrtttu08qVKx37GMSRI0f03HPPqaamRh6PR7m5uZo5c6by8/ONq8lwuZhUlw0bNqisrEwTJ04Mv/f5z39ekydPNqouw+VhUk06Ozv1gx/8QHV1ders7NQ3vvENXbhwwcjPr+FyMaku/crKyjR+/HhJMrImdqM3Owu92Zl1oTc7D73ZmXXpd616syOHWgAAAAAArobjHj8GAAAAAOBqMdQCAAAAAIzFUAsAAAAAMBZDLQAAAADAWAy1AAAAAABjMdQCAAAAAIzFUAsAAAAAMBZDLRAjVqxYobKyMknSmTNn9OCDD6qiosLmqAAAiF/0ZmBkuCzLsuwOAsCnV19fr7lz52r9+vVavny5fvSjH+nuu++2OywAAOIWvRkYGR67AwBwbeTm5mrOnDn66le/qpdeeommCQCAzejNwMjg8WMgRjQ1NWnHjh1KS0vTuHHj7A4HAIC4R28GRgaPHwMx4OLFi3riiSf0zDPPqKmpSdu3b9e6devsDgsAgLhFbwZGDndqAcN1dHTo6aef1oIFC1RSUqJHH31Up0+f1p49e+wODQCAuERvBkYWd2oBAAAAAMbiTi0AAAAAwFgMtQAAAAAAYzHUAgAAAACMxVALAAAAADAWQy0AAAAAwFgMtQAAAAAAYzHUAgAAAACMxVALAAAAADDW/wHZWu0TAleWZAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "tags": [] } } ] } ] }