{ "cells": [ { "cell_type": "markdown", "id": "7337e7c3-fef0-41be-92e0-3680f66544b7", "metadata": {}, "source": [ "# This notebook shows how to visualize how CO2 density relates to pressure and temperature using an EOS" ] }, { "cell_type": "code", "execution_count": 1, "id": "291862ad-de8c-42de-9ff1-3d8ef5f225d6", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'0.0.91'" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import numpy as np\n", "import pandas as pd\n", "import matplotlib.pyplot as plt\n", "import DiadFit as pf\n", "import CoolProp.CoolProp as CP\n", "pf.__version__" ] }, { "cell_type": "markdown", "id": "dbbddb92", "metadata": {}, "source": [ "### Install DiadFit if you havent already! You might also have to install CoolProp if you want to use Span and Wagner EOS - the error message will give you instructions, else reach out" ] }, { "cell_type": "code", "execution_count": 3, "id": "fad7fd3d", "metadata": {}, "outputs": [], "source": [ "#!pip install --upgrade DiadFit" ] }, { "cell_type": "markdown", "id": "5012a43f-ae4e-4c30-994c-0a2d06a10887", "metadata": {}, "source": [ "### Calculations of CO2 density as a function of pressure and temperature" ] }, { "cell_type": "code", "execution_count": 2, "id": "bedf37a8-ca90-4dc2-a4bc-a6d6d6f4e90b", "metadata": {}, "outputs": [], "source": [ "# Make a linearly spaced vector of pressures, \n", "\n", "# we do 2 here to get a smoother curve and combine\n", "P=np.linspace(100, 2000*10**6, 100)\n", "Psmall=np.linspace(0, 100, 100)\n", "P_combo=np.concatenate((Psmall, P))\n", "\n", "# calculating for different temperatures 1300-1000 celcius\n", "CalcD_P_1300=CP.PropsSI('D', 'P', P_combo, 'T', 1300+273.15, 'CO2')\n", "CalcD_P_1200=CP.PropsSI('D', 'P', P_combo, 'T', 1200+273.15, 'CO2')\n", "CalcD_P_1100=CP.PropsSI('D', 'P', P_combo, 'T', 1100+273.15, 'CO2')\n", "CalcD_P_1000=CP.PropsSI('D', 'P', P_combo, 'T', 1000+273.15, 'CO2')" ] }, { "cell_type": "markdown", "id": "70656156-f07e-42bb-9709-052d0c31c7c1", "metadata": {}, "source": [ "### Some things to allow us to plot in pressure vs. depth space" ] }, { "cell_type": "code", "execution_count": 3, "id": "cccf8549-b4fa-467b-a687-1a734d8d5b76", "metadata": {}, "outputs": [], "source": [ "rho=2700 # Crustal density\n", "P_axis_min=0 # Min pressure to show\n", "P_axis_max=20 # Max pressure to show\n", "# Calculating how they relate to each other\n", "D_axis_min=100000*P_axis_min/(9.8*rho)\n", "D_axis_max=100000*P_axis_max/(9.8*rho)" ] }, { "cell_type": "markdown", "id": "f0036eeb-1147-4f2a-a39e-3c05598e3c36", "metadata": {}, "source": [ "### Now lets make a pretty plot showing how CO2 density, pressure and temperature are related" ] }, { "cell_type": "code", "execution_count": 4, "id": "6d440d64-f696-4180-935e-b52c359e9b97", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAi8AAAHoCAYAAACB0Q/VAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAC3NklEQVR4nOzddVjVZxvA8S8poIJiAXZjd8+ardhdM2YsnD2nblO3d+pmzZo5u9gsFEXF7pgoM2eiAoJYIEhznvePHx48Uioogvfnus6l59fPQeHmifs2UkophBBCCCHSCeO0fgAhhBBCiDchwYsQQggh0hUJXoQQQgiRrkjwIoQQQoh0RYIXIYQQQqQrErwIIYQQIl2R4EUIIYQQ6YoEL0IIIYRIVyR4EUIIIUS6IsGLEOKDdefOHYyMjPD09Hzjcw8cOICjoyM6nQ6ASZMmUbFixdR9QKBatWps2bIl1a8rhEicBC9CiAT17dsXIyMjjIyMMDMzI0+ePDRp0oTly5frA4LUvl+7du1S7Xpjxozh+++/x9j43X6b+/HHHxk7duw7+UyEEAmT4EUIkajmzZvj5+fHnTt32LVrFw0bNmTYsGE4OTkRHR2d1o+XqBMnTnDjxg06d+78zu/VqlUrgoKC2LNnzzu/lxBCI8GLECJRmTJlws7Ojrx581K5cmXGjx/Ptm3b2LVrFytXrtQfFxQUxKBBg8idOzfW1tZ8+umn/Pvvv/r9L4ZsFi9eTP78+bGysqJz584EBgbq969atYpt27bpe3sOHTqkP//27ds0bNgQKysrKlSowMmTJ5N8bmdnZ5o2bYqFhUWix3h5eVGsWDG+/PJLdDodK1euJFu2bOzYsYOSJUtiZWVFp06deP78OatWraJQoUJkz56db775hpiYGP11TExMaNmyJRs2bHizD1cI8dYkeBFCvJFPP/2UChUq6Od5KKVo1aoV/v7+uLm54eHhQeXKlWnUqBFPnjzRn3fz5k3+/vtvXF1d2b17N56ennz99dcAjB49mi5duuh7evz8/Khdu7b+3O+//57Ro0fj6elJiRIl6N69e5I9P0eOHKFq1aqJ7r906RJ16tShc+fOLFy4UD+0FBoayty5c3F2dmb37t0cOnSIDh064ObmhpubG2vWrGHJkiVs2rTJ4HrVq1fn6NGjb/5hCiHeimlaP4AQIv1xdHTkwoULABw8eJCLFy8SEBBApkyZAJgxYwYuLi5s2rSJQYMGARAeHs6qVavIly8fAPPmzaNVq1bMnDkTOzs7LC0tiYiIwM7OLt79Ro8eTatWrQD46aefKFOmDDdv3sTR0THB57tz5w4ODg4J7jt58iROTk6MGzeO0aNHG+yLiopi4cKFFC1aFIBOnTqxZs0aHjx4QJYsWShdujQNGzbk4MGDdO3aVX9e3rx5uXfvHjqd7p3PsRFCSM+LEOItKKUwMjICwMPDg5CQEHLkyEGWLFn0Ly8vL27duqU/p0CBAvrABaBWrVrodDquXbuW7P3Kly+v/7u9vT0AAQEBiR4fFhaW4JDRvXv3aNy4MT/88EO8wAXAyspKH7gA5MmTh0KFCpElSxaDba/e29LSEp1OR0RERLJtEUKknPS8CCHe2NWrVylcuDAAOp0Oe3t7gzkqL2TLli3Ra7wIfl78mRQzM7N45yW1uidnzpw8ffo03vZcuXLh4OCAs7Mzn3/+OdbW1one58W9Etr26r2fPHmClZUVlpaWybZFCJFy0vMihHgjBw4c4OLFi3Ts2BGAypUr4+/vj6mpKcWKFTN45cyZU3/evXv3uH//vv79yZMnMTY2pkSJEgCYm5sbTIRNiUqVKnHlypV42y0tLdmxYwcWFhY0a9aM4ODgVLnfpUuXqFy5cqpcSwiRPAlehBCJioiIwN/fH19fX86dO8eUKVNo27YtTk5OfPbZZwA0btyYWrVq0a5dO/bs2cOdO3c4ceIEP/zwA2fPntVfy8LCgj59+vDvv/9y9OhRhg4dSpcuXfRzXAoVKsSFCxe4du0ajx49Iioq6q2fu1mzZhw7dizBfZkzZ2bnzp2YmprSokULQkJC3vo+Lxw9epSmTZum+DpCiNcjwYsQIlG7d+/G3t6eQoUK0bx5cw4ePMjcuXPZtm0bJiYmgDaM4ubmRr169ejfvz8lSpSgW7du3Llzhzx58uivVaxYMTp06EDLli1p2rQpZcuWZcGCBfr9AwcOpGTJklStWpVcuXJx/Pjxt37uXr16ceXKlUTn02TJkoVdu3ahlKJly5Y8f/78re/l6+vLiRMn6Nev31tfQwjxZoyUUiqtH0IIkbFNmjQJFxeXt0rz/7bGjBlDUFAQixcvfqf3+fbbbwkKCmLJkiXv9D5CiDjS8yKEyJC+//57ChYsmGrzaBKTO3du/ve//73TewghDMlqIyFEhmRjY8P48ePf+X2+/fbbd34PIYQhGTYSQgghRLoiw0ZCCCGESFckeBFCCCFEuiLBixBCCCHSFQleXrJgwQIKFy6MhYUFVapUyXBVYo8cOULr1q1xcHDAyMgIFxcXg/1KKSZNmoSDgwOWlpY0aNCAy5cvp83DppKpU6dSrVo1smbNSu7cuWnXrl283B8Zsd0LFy6kfPnyWFtbY21tTa1atdi1a5d+f0Zs86umTp2KkZERw4cP12/LiO2eNGkSRkZGBq+Xi1tmxDaDll+nV69e5MiRAysrKypWrIiHh4d+f0Zsd6FCheJ9rY2MjPTV2TNimxOlhFJKKWdnZ2VmZqaWLl2qrly5ooYNG6YyZ86s7t69m9aPlmrc3NzU999/rzZv3qwAtXXrVoP9v/76q8qaNavavHmzunjxouratauyt7dXz549S5sHTgXNmjVTK1asUJcuXVKenp6qVatWqkCBAiokJER/TEZs9/bt29XOnTvVtWvX1LVr19T48eOVmZmZunTpklIqY7b5ZWfOnFGFChVS5cuXV8OGDdNvz4jtnjhxoipTpozy8/PTvwICAvT7M2Kbnzx5ogoWLKj69u2rTp8+rby8vNS+ffvUzZs39cdkxHYHBAQYfJ337t2rAHXw4EGlVMZsc2IkeIlVvXp19cUXXxhsc3R0VGPHjk2jJ3q3Xg1edDqdsrOzU7/++qt+W3h4uLKxsVGLFi1Kgyd8NwICAhSgDh8+rJT6eNqtlFLZs2dXf/75Z4Zvc3BwsCpevLjau3evql+/vj54yajtnjhxoqpQoUKC+zJqm7/77jv1ySefJLo/o7b7VcOGDVNFixZVOp3uo2nzCzJsBERGRuLh4RGvNknTpk05ceJEGj3V++Xl5YW/v7/BZ5ApUybq16+foT6DoKAgAGxtbYGPo90xMTE4Ozvz/PlzatWqleHb/PXXX9OqVSsaN25ssD0jt/vGjRs4ODhQuHBhunXrxu3bt4GM2+bt27dTtWpVOnfuTO7cualUqRJLly7V78+o7X5ZZGQka9eupX///hgZGX0UbX6ZBC/Ao0ePiImJMajDApAnTx78/f3T6KnerxftzMifgVKKkSNH8sknn1C2bFkgY7f74sWLZMmShUyZMvHFF1+wdetWSpcunaHb7OzszLlz55g6dWq8fRm13TVq1GD16tXs2bOHpUuX4u/vT+3atXn8+HGGbfPt27dZuHAhxYsXZ8+ePXzxxRcMHTqU1atXAxn3a/0yFxcXAgMD6du3L/BxtPllkmH3JUZGRgbvlVLxtmV0GfkzGDJkCBcuXEiw2nBGbHfJkiXx9PQkMDCQzZs306dPHw4fPqzfn9Ha7O3tzbBhw3B3d8fCwiLR4zJau1u0aKH/e7ly5ahVqxZFixZl1apV1KxZE8h4bdbpdFStWpUpU6YAUKlSJS5fvszChQv11c4h47X7ZcuWLaNFixY4ODgYbM/IbX6Z9LwAOXPmxMTEJF50GhAQEC+KzaherE7IqJ/BN998w/bt2zl48CD58uXTb8/I7TY3N6dYsWJUrVqVqVOnUqFCBebMmZNh2+zh4UFAQABVqlTB1NQUU1NTDh8+zNy5czE1NdW3LaO1+1WZM2emXLly3LhxI8N+re3t7SldurTBtlKlSnHv3j0gY/+/Brh79y779u1jwIAB+m0Zvc2vkuAF7Zt8lSpV2Lt3r8H2vXv3Urt27TR6qvercOHC2NnZGXwGkZGRHD58OF1/BkophgwZwpYtWzhw4ACFCxc22J9R250QpRQREREZts2NGjXi4sWLeHp66l9Vq1alZ8+eeHp6UqRIkQzZ7ldFRERw9epV7O3tM+zXuk6dOvFSHly/fp2CBQsCGf//9YoVK8idOzetWrXSb8vobY4njSYKf3BeLJVetmyZunLliho+fLjKnDmzunPnTlo/WqoJDg5W58+fV+fPn1eAmjVrljp//rx+Ofivv/6qbGxs1JYtW9TFixdV9+7d0/0yuy+//FLZ2NioQ4cOGSwxDA0N1R+TEds9btw4deTIEeXl5aUuXLigxo8fr4yNjZW7u7tSKmO2OSEvrzZSKmO2e9SoUerQoUPq9u3b6tSpU8rJyUllzZpV/70rI7b5zJkzytTUVE2ePFnduHFDrVu3TllZWam1a9fqj8mI7VZKqZiYGFWgQAH13XffxduXUducEAleXvLHH3+oggULKnNzc1W5cmX9ctqM4uDBgwqI9+rTp49SSlteOHHiRGVnZ6cyZcqk6tWrpy5evJi2D51CCbUXUCtWrNAfkxHb3b9/f/2/5Vy5cqlGjRrpAxelMmabE/Jq8JIR2/0il4eZmZlycHBQHTp0UJcvX9bvz4htVkopV1dXVbZsWZUpUybl6OiolixZYrA/o7Z7z549ClDXrl2Lty+jtjkhUlVaCCGEEOmKzHkRQgghRLoiwYsQQggh0hUJXoQQQgiRrkjwIoQQQoh0RYIXIYQQQqQrErwIIYQQIl2R4OUVERERTJo0iYiIiLR+lPfqY2z3x9hm+Djb/TG2GT7Odn+MbYaPr92S5+UVz549w8bGhqCgIKytrdP6cd6bj7HdH2Ob4eNs98fYZvg42/0xthk+vnZLz4sQQggh0hUJXoQQQgiRrpim9QN8aKKjowHw9vbGxsYmjZ/m/QkODgbA19eXZ8+epfHTvB8fY5vh42z3x9hm+Djb/TG2GT6cdut0Oh48eEClSpUwNX13IYbMeXnFsWPHqFu3blo/hhBCCJFunTlzhmrVqr2z60vPyysKFCgAaB+8vb19Gj+NEEIIkX74+flRvXp18uTJ807vI8HLK4yNtWlA9vb25MuXL42fRgghhEh/XvwsfWfXf6dXF0IIIYRIZRK8CCGEECJdkeBFCCGEEOmKzHl5C0opoqOjiYmJSetH+WiZmJhgamqKkZFRWj+KEEKI90yClzcUGRmJn58foaGhaf0oHz0rKyvs7e0xNzdP60cRQgjxHknw8gZ0Oh1eXl6YmJjg4OCAubm5/OafBpRSREZG8vDhQ7y8vChevPg7n9kuhBDiwyHByxuIjIxEp9ORP39+rKys0vpxPmqWlpaYmZlx9+5dIiMjsbCwSOtHEkII8Z7Ir6tvQX7L/zDI10EIIT5O8t1fCCGEEOmKBC9CCCGESFckeBFCCCFEuiLBy0fgyJEjtG7dGgcHB4yMjHBxcYl3jFKKSZMm4eDggKWlJQ0aNODy5csGx0RERPDNN9+QM2dOMmfOTJs2bfDx8TE45unTp/Tu3RsbGxtsbGzo3bs3gYGByT7jzZs36devH/ny5SNTpkwULlyY7t27c/bs2ZQ0XQghRAYkwctH4Pnz51SoUIH58+cnesy0adOYNWsW8+fP559//sHOzo4mTZoQHBysP2b48OFs3boVZ2dnjh07RkhICE5OTgbJ+nr06IGnpye7d+9m9+7deHp60rt37ySf7+zZs1SpUoXr16+zePFirly5wtatW3F0dGTUqFEp/wCEEEK8cyEhIRw6dOi93MtIKaXey53SCR8fH/Lnz4+3t3e8qtLh4eF4eXlRuHBh/dJcpVSaJKyzsrJ6qxwzRkZGbN26lXbt2um3KaVwcHBg+PDhfPfdd4DWy5InTx5+++03Bg8eTFBQELly5WLNmjV07doVgPv375M/f37c3Nxo1qwZV69epXTp0pw6dYoaNWoAcOrUKWrVqsV///1HyZIl4z2PUopy5cphYWHBmTNn4q0gCgwMJFu2bAm2JaGvhxBCiPcjODiY48ePc+jQIQ4dOsTZs2f1v8wm9DM0NUmelxQKDQ0lS5Ys7/2+ISEhZM6cOVWu5eXlhb+/P02bNtVvy5QpE/Xr1+fEiRMMHjwYDw8PoqKiDI5xcHCgbNmynDhxgmbNmnHy5ElsbGz0gQtAzZo1sbGx4cSJEwkGL56enly+fJn169cnuPQ5scBFCCHE+xUaGsqJEyc4ePAgBw4c4J9//olXJidv3rz4+vq+82eR4EXg7+8PQJ48eQy258mTh7t37+qPMTc3J3v27PGOeXG+v78/uXPnjnf93Llz64951Y0bNwBwdHRMWSOEEEKkqoiICE6dOsXBgwc5ePAgp06dIjIy0uCYQoUK0aBBAxo0aED9+vUxNTUlf/787/zZJHhJISsrK0JCQtLkvqnt1WEopVSyQ1OvHpPQ8Uld58WopZRZEEKItBUVFcXZs2c5cOAABw8e5Pjx44SHhxsckzdvXj799FMaNmxIw4YNKVSokMH+VxdxvCsSvKSQkZFRqg3fpBU7OztA6zmxt7fXbw8ICND3xtjZ2REZGcnTp08Nel8CAgKoXbu2/pgHDx7Eu/7Dhw/j9eq8UKJECQCuXr1KxYoVU6U9QgghkhcTE8P58+f1w0AvFmK8LHfu3DRs2FAfsBQrVuyD+GXzg15t9DpLfBNz/PhxTE1N5QfiayhcuDB2dnbs3btXvy0yMpLDhw/rA5MqVapgZmZmcIyfnx+XLl3SH1OrVi2CgoI4c+aM/pjTp08TFBSkP+ZVFStWpHTp0sycOROdThdv/+sssxZCCJE8pRSXLl1i9uzZtG3blhw5clCtWjXGjBnD7t27CQkJwdbWlg4dOjB//nwuX76Mv78/zs7ODBo0iOLFiycYuESHh3Pv0CGOTZiAS/v276UtH3TPy4slvv369aNjx46vfV5QUBCfffYZjRo1SrAn4GMTEhLCzZs39e+9vLzw9PTE1taWAgUKYGRkxPDhw5kyZQrFixenePHiTJkyBSsrK3r06AGAjY0Nn3/+OaNGjSJHjhzY2toyevRoypUrR+PGjQEoVaoUzZs3Z+DAgSxevBiAQYMG4eTklOBkXdB6rlasWEHjxo2pV68e48ePx9HRkZCQEFxdXXF3d+fw4cPv+BMSQoiMydvbm/3797Nv3z727dsX72eitbU19evX1/eulCtXLtm6cbroaPzPnuXe/v3cO3AA3+PHiYmIACDwXTXkFR908NKiRQtatGjxxucNHjyYHj16YGJi8ka9NRnV2bNnadiwof79yJEjAejTpw8rV64EYMyYMYSFhfHVV1/x9OlTatSogbu7O1mzZtWf9/vvv2NqakqXLl0ICwujUaNGrFy5EhMTE/0x69atY+jQofpVSW3atEkyvwxA9erVOXv2LJMnT2bgwIE8evQIe3t7ateuzezZs1PpUxBCiIwvMDCQQ4cO6YOVa9euGey3tLSkbt26NGrUiE8//ZRKlSoZfA9PiFKKR5cucW//fu7u34/P4cNEvpQDDCCzvT0FGjbEvEIFJsem3HiX0k2el4TykyRkxYoVLFiwgJMnT/LLL7/g4uKCp6dnosdHREQQERsxAvj6+lK6dOnXzvMi0o58PYQQH7uIiAhOnjypD1b++ecfgyF4Y2NjqlevTuPGjWnUqBG1atUiU6ZMyV430MtL61mJ7V0JDQgw2G9ha0v+hg0p8OmnFPj0U2xLlsTIyCjJXGmp6YPueXlTN27cYOzYsRw9ehRT09dr2tSpU/npp5/e8ZMJIYQQKafT6bhw4YI+WDly5AhhYWEGx5QsWZLGjRvTuHFjGjRo8Fr5sp4HBOB94AB3YwOWIC8vg/2mVlbkq1uXAo0aUbBRI3JXrIhRMsNL71KGCV5iYmLo0aMHP/30k34Fy+sYN26cfhgF4npehBBCiA+Bt7c3e/bsYd++fezfv59Hjx4Z7M+TJ48+WGnUqNFr5VmJDA7G+/Bh/VDQo4sXDfYbm5piX6MGBRo1okCjRtjXqIFpUj02MTHg6QmbNr1NE99YhglegoODOXv2LOfPn2fIkCGAFqEqpTA1NcXd3Z1PP/003nmZMmUy6EJ79uzZe3tmIYQQ4lWhoaEcOXKEPXv2sGfPHq5evWqwP3PmzDRo0EAfsJQpUybZ5cu66Gj8Tp/mzt693N27F7/Tp1GvZMfNVaGCvmclX926mL805zEepeDmTdi3D/bvhwMH4OnTt27zm8owwYu1tTUXX4kcFyxYwIEDB9i0aROFCxdOoycTQgghEqeU4vLly/pg5ciRIwZzMY2NjalRowZNmjShSZMm1KhRAzMzs2SvG3jrFnfc3bnj7s69AweIfOWX82xFi+p7Vgo0bIhVrlxJX/DBAy1QeRGw3LtnuN/aGqpX1/a/Yx908JLcEt9x48bh6+vL6tWrMTY2pmzZsgbn586dGwsLi3jbhRBCiLT0+PFj9u7dy549e3B3d+f+/fsG+/Pnz0+zZs1o1qwZjRo1ileaJSHhgYHcO3CAu7EBy6vzVixsbSnYpAmFmjShQOPG2BQsmPQFg4PhyBEtGNm3Dy5dMtxvbg61a0PjxtCoEVStCv7+8LGXB0huia+fnx/3Xo38hBBCiA9MdHQ0p06d0veunD17lpcX+1paWtKgQQN9wFIydvVOUmKiovA/c0bfu+J/5gzq5ZVGZmY41K5NoaZNKdS0KbkrVcI4qWXR0dFw9iy4u2uv06e1bS8YGUHFilqw0rgxfPIJvINSNa/jgw5eGjRoQFIruV/kKEnMpEmTmDRpUuo+lBBCCPEa7t69y+7du9mzZw/79++PN6eybNmy+mClbt26yaZ8UErph4Lu7t2b4FCQraMjhZo2pWDTpuSvXx/zLFmSfsg7d+KClf374dWs5kWLar0qjRtDw4aQM2fi1woPAO+tSd8vlXzQwYsQQgiRXkRFRXH8+HF27tyJm5sbV65cMdhva2tLkyZNaNasGU2bNiVv3rzJXjMyOJi7+/dzZ/fupIeCmjalYJMmWCc3ZBMcDAcPxgUsN24Y7s+WTQtUmjbV/kxqvmh0GDw8Cv57wW8vBP4Lj5NtUqqQ4EUIIYR4S35+fuzatQs3Nzf27t1r0LtiYmJCzZo19b0rVapUee1stl67duG1eze+x46hi4rS73/joaCYGPDwiAtWTp40HAoyMYFatbRgpWlTbd5KYtdTOnj6rxas+O+FgKOgizA8Jltp4EqCp6cmCV6EEEKI1xQTE8OZM2dwc3Nj586dnD9/3mB/rly5aNGiBS1btqRp06avNdE24tkz7u7bh9euXdzZvZtgHx+D/dmKFqVwixYUataM/A0aJD8UdO9eXLCyb1/8JczFisUFKw0agI1N4tcK9dF6Vfz3gv8+iHhouN8qH9g1iX01gkeRwEc+YVekjiNHjjB9+nQ8PDzw8/NLsMzCli1bWLx4MR4eHjx+/Jjz58/Hq8gdERHB6NGj2bBhg7620YIFCwxSQD99+pShQ4eyfft2QKttNG/evGQzPN68eZPJkyezd+9eHj58iIODAzVr1mTUqFFUrVo1NT4GIYR4K48fP2bPnj3s3LmT3bt38+TJE/0+IyMjqlWrRsuWLWnZsiVVqlRJtrChUoqHFy7oe1fuHz+O7qXeEFMLC/I3bEjhFi0o3KIF2YsVS/oBw8Lg8GHYtUsLWP77z3C/tbU2b+VFwFKkSOLXig6FgMNwfzf4u8OzV65lmgVyNwD7JmDXFKxLahN59QwDr3dFgpePwOtU537+/Dl16tShc+fODBw4MMFjhg8fjqurK87OzuTIkYNRo0bh5OSEh4eHviu0R48e+Pj4sHv3bkCrKt27d29cXV0Tfb6zZ8/SqFEjypYty+LFi3F0dCQ4OJht27YxatQoqSothHivlFKcP38eNzc33NzcOH36tEG9oGzZstGsWTNatmxJ8+bNyZ07d7LXDA8MNOhdCXllaXT2EiW0YKV5c/LVr4+ZpWXSF7x5UwtWdu2CQ4e0AOYFY2OoUSMuWKleHRIrmaMUPLuqBSt+uyHgiOFQkJEx2FYD+6Za70qOGmBinmx73zUJXlJIKUVUaOh7v6+ZlVWyy+heeJ3q3L179wbgzp07Ce4PCgpi2bJlrFmzhsaNGwOwdu1a8ufPz759+2jWrBlXr15l9+7dnDp1iho1agCwdOlSatWqxbVr1yhZsmS86yql6Nu3L8WLF+fo0aMGv7FUrFiRYcOGvVYbhRAiJUJDQ9m3bx/bt2/Hzc0NPz8/g/3ly5enZcuWtGrVipo1ayZbP08pxcN//+W2mxteu3Zx/+RJg4y2ppaWFPj0U33vSrakekPAsHdl1674E23z5oUWLaB5c62XJane7shA8N+vBSt+u7WhoZdZFQCH5lrPit2nYJ780Nf7JsFLCkWFhjI3ufHHd2BoSAjmmTO/t/t5eHgQFRVF06ZN9dscHBwoW7YsJ06coFmzZpw8eRIbGxt94AJQs2ZNbGxsOHHiRILBi6enJ5cvX2b9+vUJdrW+TkExIYR4G35+fuzYsYPt27ezb98+wsPD9fsyZ85MkyZNaNmyJS1atHitCslRYWHc27+f2zt2cGvHDkJ8fQ322zo66oOVfHXrYprM0miD3pWDB+Gl58PUVMuz0qKF9ipb9pXhm5coHTw5FxesPDoF6qXSACYWkLs+2DfXXvGGgj48EryI1+Lv74+5uXm8yWd58uTB399ff0xC3ae5c+fWH/OqG7G/PTg6OqbyEwshhCGlFBcvXmT79u1s376df/75x2B/oUKFaNOmDU5OTtSrV8+g7l1ign18uL1zJ7d27ODe/v1EvzR8Y2plRcFGjfQBi02hQklfLCxMGwJ6EbC8lGEeiOtdadlS612xtk7iWv7g5w5+e7S5KxGGxRyxdowLVnLXA9Nkhqk+MBK8pJCZlRVDQ0LS5L4fAqWUwfBVQkNZrx7z6r7EzhNCiJSKjIzk8OHDbN++HVdXV+7evWuwv0aNGrRp04bWrVtTtmzZZL8XKZ0O/7NnueXqyu0dOwjw9DTYn7VAAYo6OVHEyYkCDRsm37ty+zbs3Jny3hVdNDw6ETd35anhKihMs4JdY204yL4ZZE6mNMAHToKXFDIyMnqvwzdpxc7OjsjISJ4+fWrQ+xIQEEDt2rX1xzx48CDeuQ8fPiRPnjwJXrdEiRIAXL16Nd7qJiGEeBtPnjzBzc2N7du3s3v3boKDg/X7LC0tadKkCa1bt8bJyQk7O7tkrxcZHMydvXu5vWMHt3fuJDQgIG6nkREOtWrpA5acyQVA0dFarhVXV9ixA16pGE2+fHHBSnK9K+EP4f4uuO+m9bBEBRruz145NlhpDjlrgnHyxRzTCwlexGupUqUKZmZm7N27ly5dugDaePGlS5eYNm0aALVq1SIoKIgzZ85QvXp1AE6fPk1QUJA+wHlVxYoVKV26NDNnzqRr167x5r0EBgbKvBchRLJu3LihHw46fvw4MS9NjrWzs8PJyYk2bdrQqFEjrF6j5zrw9m1u79zJ7R078D50iJjISP0+c2trCjVrRlEnJwq3aJF8NeanT2H3bi1Y2bXLMO+KiQnUrasNBbVoAWXKJD135el58N2pBSyPzwAvldDJlAPsmoFDC21lkGXCvzRmBBK8fASSq84N2m8q9+7d01c2vXbtGqD9p7ezs8PGxobPP/+cUaNGkSNHDmxtbRk9ejTlypXTrz4qVaoUzZs3Z+DAgSxevBjQlko7OTklOFkXtJ6rFStW0LhxY+rVq8f48eNxdHQkJCQEV1dX3N3dZam0ECKeF8uZt2zZwtatW+Ol4i9fvjytW7emTZs2VK1aNfncKzodfmfOcNPFhVuurjx+5XrZihWjaOvWFHVyIu8nn2BinsRyYaXg+vW43pVjx7RMty/Y2mrBipMTNGuW9MqgqGdacrgXAUv4K/MHs1cEh1baK0d1ME46g2+GoYQBb29vBShvb+94+8LCwtSVK1dUWFhYGjzZ2zt48KBCC88NXn369NEfs2LFigSPmThxov6YsLAwNWTIEGVra6ssLS2Vk5OTunfvnsG9Hj9+rHr27KmyZs2qsmbNqnr27KmePn2a7DNeu3ZNffbZZ8rBwUGZm5urggULqu7du6tz584lek56/XoIId5OdHS0Onz4sBo+fLgqWLCgwfcqU1NT1aRJEzV37lzl5eX1WteLCg9Xt9zc1J5Bg9QCOzs1HfSvGSYmyrlBA3Vmxgz1+L//kr9YRIRS+/YpNXy4UsWKKaWFMHGvMmWU+u47pY4eVSo6OvHr6HRKBf2n1JWZSu37VKkNZkqtI+71V2alDrdT6sZSpZ77vN4H9x4l9TM0NRkplUTZ5o+Qj48P+fPnx9vbO97SuPDwcLy8vChcuHCy1T/FuydfDyEyvoiICA4cOMDWrVtxcXHh4cO49PRWVla0aNGC9u3b06pVq9caYg4PDMTLzY2b27Zx282NqJcWXJhbW1OkZUuKtW1LoebNsUjueo8egZub1ruyZw+8XOHZ3FxLve/kpL2SKnAYEw4PDsP9ndor5Lbh/qzFtZ6VvK0gV10wSX4VVFpJ6mdoapJhIyGEEB+UkJAQdu/ezZYtW9i5c6dBscNs2bLRpk0bOnToQJMmTV5r/kqwjw83t23j5rZteB88aJCKP4uDA0XbtqV4u3bkb9Ag6eEg0IaDXFxg2zZt4u3Lv//nzg2tWkHr1lpF5qxZE79OeIA2FOS7XVvSHPNSslNjcy3vikMrcGgJ1sWTbePHRoIXIYQQae7Jkye4urqyZcsW3N3dDRLG2dvb065dOzp06ED9+vUxM0t61YxSikeXL3Nr2zZuuLjw4OxZg/05SpemWLt2FGvbFruqVTFKaj6MTgf//BMXsLy6OqhiRS1YcXLSKjIndi2ltDpBvtvBZzs8OonBZFtLh7hgxa4xmL3/5KfpiQQvQggh0sT9+/dxcXFhy5YtHDp0yGCFUNGiRWnfvj0dOnSgRo0ayU641cXEcP/kSW66uHDTxYXAW7fidhoZkbd2bX3Akr14Mj0ZkZFazpUXAcvLpQLMzODTT6FtWy1oSWpoRBcND4/HBSwhrySds60CedtA3tbaxFvJd/XaJHgRQgjx3vj4+LB582b+/vtvTpw4YbCvfPnydOjQgfbt21OuXLlkE8bFREVx78ABrm/axM1t2wh7aT6MSaZMFGzShGJt21K0dWsyJ5JrSi8oSFvG7OKizWN5KTcMWbNqq4PatdOWM9vYJH6dqGdazhWf7dr8lciXlkUbm0OeRpCvDeR1Aqt3Nycko5PgRQghxDvl4+PDpk2b2LhxY7yApVatWvqApWjRosleKzoigrt793J90yZubd9O+Es5UyyyZ6dIq1YUa9eOQs2aYZ5c3bn792H7di1gOXAAoqLi9tnba70rbdtCw4aQVKmA597g6wo+2yDgIOheuk6mHODgpAUsdk3ALIl5MOK1SfAihBAi1SUVsNSpU4cuXbrQsWNH8ubNm+y1osLCuLN7N9c3b+aWqyuRL03gtcqTh+Lt21OiY0fy1a+PSVLzYZSC//7TghUXFzhzxnC/o6PWu9KuHVSrlvT8lafntd4V3+3xU/FnLQ752mpDQjlrfzy5V94jCV6EEEKkihcBy99//83JkycN9n3yySd07tz5tQOWyJAQvHbt4vqmTdzeuZOo58/1+7LkzUvxDh0o0akTeevUwdgkieBAKTh/HjZtgs2btdVCL6tVSwtW2raFRJJpAqCLgYfHwGcreG+F0Htx+4yMtSAlbxuth8U6ieuIVCHBixBCiLfm6+vLxo0b4wUsRkZG1KlT540Clohnz7i9YwfXN23Ca9cuol9acZS1QAFKdOpEyU6dsK9RI+kVQkppvSqbNmmvO3fi9pmbazWD2rXTJtza2yd+nZgI8N+vBSw+2yAibk4NJlZagcN8bbUVQhbJlAgQqUqCFyGEEG/k0aNHbNq0CWdnZ44cOWJQHf7lISEHB4dkrxX+9Ck3t2/n+qZN3HV3N6ghlK1oUUp06kTxjh21Jc1JTeDV6eDEibgeFh+fuH2WltqE206dtD+TKnYYFQJ+u7TeFd8dEP3SxF3z7FrvSv4O2vwVU8tk2yfeDQlehBBCJOvZs2e4uLjg7OzM3r17iX4p0VudOnXo2rXrGwUsN1xcuPb339zbt88gaZxtyZKU6NyZEh07kqtCheQrNB85ogUrW7aA/0t1f7Jk0XpWOnaE5s0hc+bErxPxWJtw671FSxini4jbZ2kP+dprAUvuehmqMnN6JsHLR+DIkSNMnz4dDw8P/Pz82Lp1K+3atdPvj4qK4ocffsDNzY3bt29jY2ND48aN+fXXXw2+EUVERDB69Gg2bNhAWFgYjRo1YsGCBQYpoJ8+fcrQoUPZvn07AG3atGHevHnJpu2+efMmkydPZu/evTx8+BAHBwdq1qzJqFGjqFq1aqp+HkKI1xMWFoabmxsbNmxg586dBonjKlWqRPfu3enatau+wGtSIoODubl9O/85O3Nnzx50L63syVmuHCU6dqREp07kKF066YAlKkpbGbR5M2zdqqXof8HGRpu70qkTNGkCSZUNCfUFHxctYAk4DOqlwolZikH+2IAlR3VtTov4oEjw8hF4/vw5FSpUoF+/fnTs2DHe/tDQUM6dO8ePP/5IhQoVePr0KcOHD6dNmzacfSkz5fDhw3F1dcXZ2ZkcOXIwatQonJyc8PDwwCR2wlyPHj3w8fFh9+7dgFZVunfv3ri6uib6fGfPnqVRo0aULVuWxYsX4+joSHBwMNu2bWPUqFFSVVqI9ygqKop9+/axYcMGXFxcCH4p30nJkiXp3r073bp1S7RSvMG1QkO5vXMn//31F147dxrMYclZtiwlu3alZOfO2CZ3rYgI2LtXGxLavh1eWh5Njhza/JWOHbW5LEml9392I3bC7RZ4fNpwX7YKWrCSvz3YlJWEcR84Kcz4ijcuzKgUhIYmcKV3zMrqrf5zGRkZxet5Scg///xD9erVuXv3LgUKFCAoKIhcuXKxZs0aunbtCmjZMfPnz4+bmxvNmjXj6tWrlC5dmlOnTlGjRg0ATp06Ra1atfjvv/8S/GanlKJcuXJYWFhw5syZeFk0AwMDE+21kcKMQqQOnU7H8ePHWbduHZs2beLx48f6fQUKFKBbt250796dCskN46DlYbmzezf//fUXt7ZvN1gllL1ECUp27Ypj167kLFMm6YeKiAB3d/jrLy1geTlpXJ480L691sNSvz6YJvF7+LMb4L0R7v4Ngf++tMMIctaKC1iyFEn6ecRrkcKM6UVoqDa2+r6FhCQ9hptCQUFBGBkZ6QMHDw8PoqKiaNq0qf4YBwcHypYty4kTJ2jWrBknT57ExsZGH7gA1KxZExsbG06cOJFg8OLp6cnly5dZv359gum/X6dKrBDi7fz333+sXbuWdevWceelFTm5c+emS5cudO/enZo1ayabmj8mKop7+/fzn7MzN11ciAgK0u+zLlgQx27dKNm1K7krVny9ISFnZ21I6KXrkDcvdOigBSx16kBSy6NfBCz3NsJTz7jtRiaQ51MtYMnXVpvPItIlCV5EPOHh4YwdO5YePXpgHTsr39/fH3Nzc7Jnz25wbJ48efCPnSTn7+9P7ty5410vd+7c+mNedePGDQAcHR1TswlCiEQ8ePAAZ2dn1q5dazAsnDVrVjp27EiPHj1o2LAhpkn1ZgBKp8P3+HGurl/P9Y0bCXuptyaLg4O+h8WuevWkA5aYGDh8WOth2bwZXroODg7QuTN06QI1ayaeNA6SDljsGkOBzpCvnZbxVqR7EryklJWV1guSFvd9B6KioujWrRs6nY4FCxYke7xSyuAbU0LfpF495tV9iZ0nhEgdoaGhuLi4sHbtWtzd3fUFEE1NTWnevDm9evWiTZs2WFomv/T34aVLXF23jv82bODZ3bv67Va5c1OiUyccu3Ujb506yVdqPnFCC1g2bTJcJZQrl9a70q0bfPJJ0gFL8E0tWLn3dwI9LI2gYBcJWDIoCV5SysjonQ7fvE9RUVF06dIFLy8vDhw4oO91AbCzsyMyMpKnT58a9L4EBARQu3Zt/TEPHjyId92HDx+SJ5GiaCVKlADg6tWrVKxYMRVbI8THLSYmhoMHD7JmzRq2bNlCyEu/ZNWoUYNevXrRtWtXcuVKPrnaM29v/tuwgavr1vHwwgX9dvOsWSneoQOlevakQMOGGCfVW6MU/POPFrD8/bdhHpbs2bUJt127QoMGSc9h0QcsGw3T8r8IWF70sFjkTLZdIv2S4EUAcYHLjRs3OHjwIDlyGP6mUqVKFczMzNi7dy9dunQBwM/Pj0uXLjFt2jRAK7AWFBTEmTNnqF69OgCnT58mKChIH+C8qmLFipQuXZqZM2fStWvXN5qwK4SI7+LFi6xevZr169dz//59/fYiRYrQq1cvevbsqf+lISnhT59yfdMmrqxbh8+RI1rwARibmVGkZUsce/SgaOvWmCXXW3P1KqxbBxs2wO3bcduzZtUm3XbtCo0bJ71KKPjWSz0srwYsn0KBLhKwfGQkePkIhISEcPPmTf17Ly8vPD09sbW1pUCBAkRHR9OpUyfOnTvHjh07iImJ0c9RsbW1xdzcHBsbGz7//HNGjRpFjhw5sLW1ZfTo0ZQrV47GjRsDUKpUKZo3b87AgQNZvHgxoC2VdnJySnRZpZGREStWrKBx48bUq1eP8ePH4+joSEhICK6urri7u8tSaSGS8fjxYzZs2MDKlSvx8PDQb8+ePTtdu3ald+/e1KpVK/mVQuHh3Nqxg6vr1uHl5maQ7TZfvXqU6tmTEp06YWlrm/QD+fpqwcq6deDpGbfdygratNEClubNk87D8twb7v0Fd53hSVyb4gKWzlryOAlYPk5KGPD29laA8vb2jrcvLCxMXblyRYWFhaXBk729gwcPKiDeq0+fPkoppby8vBLcD6iDBw/qrxMWFqaGDBmibG1tlaWlpXJyclL37t0zuNfjx49Vz549VdasWVXWrFlVz5491dOnT5N9xmvXrqnPPvtMOTg4KHNzc1WwYEHVvXt3de7cuUTPSa9fDyFSQ1RUlNqxY4fq1KmTMjc31/+fNTMzU+3bt1dbt25VERERyV5Hp9Mp7yNH1O4BA9Qca2s1HfSvFeXKqVO//qqC7t5N/oGePFFq6VKlGjRQyshIKa2vRilTU6WcnJTasEGpkJCkrxH2UKnrC5Vyr6vUOuJe602U2t9YqRtLtGPEByupn6GpSfK8vOKN87yINCNfD/ExunLlCitXrmTNmjUGq/gqVapE37596dGjBzlzJt8bEXj7NpdXr+bK6tUEeXnpt2fNn59SPXpQqmdPcpUrl/RFwsNhxw6th8XNDV7qqeGTT6BnT221UI4kJsxGBWuZbu9sAP+9oOJKBZCrLhTqDvk7SeHDdELyvAghhAC0shvOzs6sXLmSM2fO6LfnzJmTXr160bdvXypUqJDsdSKCgri2cSNXVq/G5+hR/XbzrFkp0bkzZT77jHx16ya9UigmBg4e1AKWLVvg2bO4fWXLagFL9+5QsGAS1wiH+7vgznq4v0N7/0L2ylrAUqArZM6fbJvEx0mCFyGE+ADFxMSwb98+VqxYgYuLCxERWrFAU1NTWrVqRd++fWnZsiXmSU10BXQxMdzdu5fLq1Zx08UlLkW/kRGFmjShTJ8+FGvXDrPk0i9cvAirV2tBi59f3Pb8+aFHDy1oSaqnRhcNDw7A3Q1aev6ol4Ie65JQsDsU7Kb9XYhkSPAihBAfkLt377J8+XKWL1+Oz0vLicuVK0e/fv3o2bNngskgX/Xw0iWurF7NlbVref5SsGFbqhRl+vShdK9eZM2bN+mLPHgA69drQcvLE29tbbXhoJ49tWy3ifXUKKVNtr2zVgtawgPi9lnl14KVgt0he0WpJSTeiAQvQgiRxqKioti+fTtLly7F3d1dn7zR1taWnj170q9fPyoml1ofCHvyhKvr1nF55UoenDun326ZIweO3btTpk8f8lSpkvR1wsK0WkKrV8OePdowEYCZGbRuDZ99Bi1aJL20OeQO3FmnBS3P/ovbnimntkqoYHfIVUeqNYu3JsGLEEKkkevXr7Ns2TJWrlxJQEBcr0SjRo0YMGAA7du3J1OmTEleQ+l03N2/n0vLl3Nj61ZiYoeXjM3MKNKqFWX69KFIy5aYJBVsKAXHj2sBy99/G9YUqlFDC1i6dk164m3kUy0Xi9daeBg3nwYTCy0HS6HeYN8EjM2SbI8Qr0OCFyGEeI/Cw8PZvHkzS5cuNchhZGdnR79+/fj8888pWrRostd5du8el1as4NKKFQZp+nNVqEC5/v1x7NEDq+RWHd26BWvWaK+XE8gVKAC9e2uvRHI0ARATAffdtB4W3x2ge7HayEjLxVK4t1ax2cw68WsI8RY+6ODlyJEjTJ8+HQ8PD/z8/Ni6dSvt2rVL9PgtW7awcOFCPD09iYiIoEyZMkyaNIlmzZq9v4cWQogEXLx4kaVLl7J27VqePn0KgLGxMS1atGDgwIG0bNkSM7OkeyWiIyK46eLCxWXLuLtvnz7rbSYbG0r17EnZ/v3JU7ly0sNCwcGwcSOsWAHHjsVtz5JFqyn02WdQv37S81geHtcClnt/az0uL2QrD4V6QaEeYJXMfBohUuCDDl6eP39OhQoV6NevHx07dkz2+CNHjtCkSROmTJlCtmzZWLFiBa1bt+b06dNUqlTpPTyxEELECQ0NxdnZmSVLlnD69Gn99gIFCvD555/Tr18/8udPfjnwwwsXuLhsGVfWriX8yRP99vwNG1Lu888p3qFD0mn6XwwLLV+uDQs9f65tNzKCJk20gKVdu6TrtIV4gddquL0KnsflhcHSAQr11IKW7OWTbYsQqeGDDl5atGhBixYtXvv42bNnG7yfMmUK27Ztw9XVVYIXIcR7c+3aNRYtWsTKlSsJDAwEtCXObdu2ZcCAATRp0gQTE5MkrxERFMTV9eu5uGwZD15K+Z81Xz7K9O1L2X79yFakSNIPcv++No9l+XK4cSNue/Hi0L+/NiyU1IqjqBDw3gy3V0LAobjtplkgf0dtWCh3AzBOui1CpLYPOnhJKZ1OR3BwMLZJ1OGIiIjQ508ACA4Ofh+PJoTIYKKjo9m+fTsLFixg//79+u2FCxdm8ODB9O3bN9Hq6i8opfA/c4Z/Fy/mP2dnosPCAG3ybbG2bSn3+ecUbNIE46QCn8hILevt8uWwaxfodNr2zJmhSxctaKlTJ/GlyUoHAUfBa6U2ATc6tpcGI7BrBIX7avNYTJPJCyPEu/ROiw+kIkBt3br1jc6ZNm2asrW1VQ8ePEj0mIkTJyZY0ycj1TY6fPiwcnJyUvb29ol+jhMnTlQlS5ZUVlZWKlu2bKpRo0bq1KlTBseEh4erIUOGqBw5cigrKyvVunXreJ/TkydPVK9evZS1tbWytrZWvXr1eq3aRjdu3FB9+/ZVefPmVebm5qpQoUKqW7du6p9//kn0nPT69RAZi6+vr5o0aZJycHDQf/8wMjJSrVu3Vm5ubiomJibZa4QHBanzCxaolRUqGNQWWl6mjPpn1iz1PCAg+Qe5eFGpESOUypkzrq4QKPXJJ0otX65UcHDS5wffVurCJKVcChvWFdpeXKlLk5UKuZf0+UKo91fbKMP2vGzYsIFJkyaxbdu2JBM6jRs3jpEjR+rf+/r6Urp06ffxiO/N68wdKlGiBPPnz6dIkSKEhYXx+++/07RpU27evEmuXFpNkeHDh+Pq6oqzszM5cuRg1KhRODk54eHhoe8C79GjBz4+PuzevRvQqkr37t0bV1fXRJ/v7NmzNGrUiLJly7J48WIcHR0JDg5m27ZtjBo1SqpKiw+OUooDBw6wcOFCXFxciInNhZI7d24GDBjAoEGDKJhUevxY/mfPar0sGzYQFTsPxSRTJkp26UKFwYNxqF076cm3QUFa9ebly+Gff+K229tDnz7Qt2/Sq4USHRbKqiWQK9IXctaSBHLiw/NOQ6NUxBv0vDg7OytLS0u1Y8eON77PG1eV1umUigp5/y+d7o3bptTrf45BQUEKUPv27VNKKRUYGKjMzMyUs7Oz/hhfX19lbGysdu/erZRS6sqVKwow6LE5efKkAtR///2X4H10Op0qU6aMqlKlSoK/oSbVayM9L+J9e/Lkifr9999VyZIlDXpq69atqzZs2PBaVZwjnj1TnosXq9WVKxv0sixzdFRnf/9dhT5+nPQFdDqlTp1Sqn9/paysDKs3d+ig1I4dSkVFJX3+g6NKneyr1F+ZX+plMVJqfxOlvNYpFfX8DT8ZITTS8/KWNmzYQP/+/dmwYQOtWrV69zeMCYW/s7z7+7yqSwiYJrEyIAUiIyNZsmQJNjY2+mJvHh4eREVF0bRpU/1xDg4OlC1blhMnTtCsWTNOnjyJjY0NNWrU0B9Ts2ZNbGxsOHHiBCUT+A3Q09OTy5cvs379eowTWJqZLVu21G+gEG/o4sWLzJ07l3Xr1hEWOw8la9as9O7dmy+//JKyZcsme40H589zYfFirqxbR1RICAAm5uaU6NSJ8oMHawURk+tlWbcOFi+GCxfitpcqBQMHaqn6kyobEPYgdrXQMnh2LW571uJaD0uh3lIIUaQbH3TwEhISws2bN/Xvvby88PT0xNbWlgIFCjBu3Dh8fX1ZvXo1oAUun332GXPmzKFmzZr6cvGWlpbY2NikSRvSkx07dtCtWzdCQ0Oxt7dn79695IxNcuXv74+5uTnZs2c3OCdPnjz6z9nf3z/BIbrcuXPrj3nVjdgVEI6OjqnZFCFSLCYmBldXV+bOncvBgwf128uVK8dXX31Fz549yZo1a5LXiAoL4z9nZ/5duBD/l4Z1shcvTvlBgyjTt2/SieSUgtOnYckScHbWUvcDWFhok28HDYLatRMf1tHFgL873PoTfLaDita2m2bWqjYX7Q85kzhfiA/UBx28nD17loYNG+rfv5ib0qdPH1auXImfnx/37t3T71+8eDHR0dF8/fXXfP311/rtL45/J0ystF6Q980k9Wf6N2zYEE9PTx49esTSpUvp0qULp0+fTnLOkFLK4LfFhH5zfPWYV/cldp4QaSEwMJDly5czf/58vLy0fCYmJiZ07NiRb775hjp16iT77zXQy4t/Fy7k4rJl+rwsxmZmFG/fngqDB5O/YcOkrxEYCGvXakHLxYtx20uXhsGDoVcvrThiYkLuwO3lcHsFhMYVdyRHDSg6AAp2BbOkAy8hPmQfdPDSoEED/Q+3hLwakBw6dOjdPlBCjIze2fDN+5Y5c2aKFStGsWLFqFmzJsWLF2fZsmWMGzcOOzs7IiMjefr0qUHvS0BAALVr1wa09OYPHjyId92HDx8mukS0RIkSAFy9epWKFSumfqOEeE3Xrl1j7ty5rFq1iuexk2dtbW0ZNGgQX331VbLJ5JROx919+zg/fz63duzQZ7+1LliQCl98Qdn+/cmc1LCOUnDqlBaw/PWXYS9L165aL0utJCbPxkSAjwvcWgb++9Cm4wDmtlD4Myj6OWRLfnhLiPTggw5eRNpSSulz4FSpUgUzMzP27t1Lly5dAPDz8+PSpUtMmzYNgFq1ahEUFMSZM2eoXr06AKdPnyYoKEgf4LyqYsWKlC5dmpkzZ9K1a9d4814CAwNl3ot4Z3Q6HXv27GHu3Ln6FXIAZcuWZdiwYfTs2RPLpDLXoiWTu7xqFef/+IOn16/rtxds0oRKQ4ZQpFWrpPOyPHum9bIsXAiXLsVtL1tWC1h69YJXhmsNBF7SApY7ayDicdx2uyZawJKvHZgkXdxRiPRGgpePQHJzh54/f87kyZNp06YN9vb2PH78mAULFuDj40Pnzp0BsLGx4fPPP2fUqFHkyJEDW1tbRo8eTbly5WjcuDEApUqVonnz5gwcOJDFixcD2lJpJyenBCfrgjZctGLFCho3bky9evUYP348jo6OhISE4Orqiru7uyyVFqkuJCSEVatWMW/ePK5d0yavGhkZ0aZNG4YOHUrD5IZ1gIeXLuH5xx9cWbNGv8zZPGtWyvbrR8WvvsI2qSXKAJcvw4IFWgbc2Am8WFrG9bLUrJl4L0t0GHhvghuL4NGJuO1W+aBIP+2VpfBrfRZCpEvvdC1TOvTGS6XTgYMHDyaYiK9Pnz5KKa1d7du3Vw4ODsrc3FzZ29urNm3aqDNnzhhcJywsTA0ZMkTZ2toqS0tL5eTkpO7dM0xc9fjxY9WzZ0+VNWtWlTVrVtWzZ8/XSlJ37do19dlnn+mfoWDBgqp79+7q3LlziZ6TXr8eIu3cvXtXjRw5UllbW+v/H1hbW6sRI0aoW7duJXt+dGSk+m/jRuXcoIFhMrnSpdX5BQtUxLNnSV8gMlKpv/5Sqn59w0Ryjo5KzZmjVHL/V4L+U8pjpFIbbeOWOK83VepIB6V83ZSKiX7tz0KId+F9LZU2UiqJSSUfIR8fH/Lnz4+3tzf58uUz2BceHo6XlxeFCxfGwsIijZ5QvCBfD/G6zp8/z4wZM/jrr7/0CeVKlCjB0KFD6dOnD1myJJ3uIPTRIy4sXoznwoWE+PoCYGRiQrF27ag0ZAj569dPuqfG11eby7JkCbxYeWdiAm3bwtdfQ8OGScxlidTmstxcBA/iVj2RuSAUGwRF+oOl3et+FEK8U0n9DE1NMmwkhMiQlFLs2bOHGTNmGNQaatSoESNHjqR58+YJ5hZ62aMrVzg3ezZX1qwhOjwcAMtcuagwaBDlBw/GOqlJvErBwYPa0JCLC8QGTeTJow0LDRoESX1zD/GCm0u1vCzhAdo2I2NwaAXFvgD7ZlIQUXy0JHgRQmQokZGRbNiwgRkzZnApdgKsiYkJXbt2ZfTo0clWmFdKccfdHY/ff+fOnj367XkqV6by8OGU7NIF00xJTIANCtLmsSxYAP/9F7e9Xj346ito3x7MzRM+VxcN93dqc1n89qBfMWRpD0UHasucJZGcEBK8CCEyhsDAQJYsWcKcOXO4f/8+AFmyZGHgwIEMHz6cAgUKJHl+VFgYV9euxWP2bB5fuaJtNDKieLt2VBkxgryffJL00NC1azB3LqxaBbETeMmcGT77DL78EsqVS/zc0Ptwa6nW0xLmG7fdrikU/wLyOoGx2et8DEJ8FCR4EUKka/fu3WP27NksXbqUkNhVOw4ODgwbNoxBgwYlu9Q+xM8PzwUL+HfRIsIePQLALEsWyn3+OZWHDiVbkSKJn6wUuLvDnDmwa1fc9lKltLksvXuDtXXi5z48Btfng/eWuOy3mXJpmW+LDoSsRV/3YxDioyLBy1uQOc4fBvk6fNw8PT2ZPn26wSTcsmXLMnr0aLp37455YkMzsR6cP4/H77/zn7MzuqgoQEsoV3noUMp9/jmZkiop8vw5rFmjBS0vhoaMjKB1axg2LOkJuNGhcGe9FrQE/hu3PVddKP4V5G8veVmESIYEL2/AzEzrtg0NDU02cZV490JDQ4G4r4v4OBw7dozJkycbJJVr1KgRo0ePplmzZkkO7Sidjls7dnB21ix8XsoflLdOHSoPH07xdu0wNk3i2+K9e/DHH9qqocBAbVvWrNC/P3zzDRRNoqck5DZcX6AllIuKPdfEEgr1ghJfQ/YKr9F6IQRI8PJGTExMyJYtGwEB2sx/KysrqcmTBpRShIaGEhAQQLZs2TBJKnupyBCUUri7uzN58mSOHj0KgLGxsX4SbuXKlZM8Pzoigqvr1vHP9Ok8ie0pMTIxoWSXLlQZPhz72IzQidwcjh/Xelm2bo1bNVS0qBaw9OuXxNCQDvzctV6W+27oJ+BmKQLFv4ai/cA8iey5QogESfDyhuzstHwKLwIYkXayZcum/3qIjEmn0+Hi4sKUKVPw8PAAwNzcnL59+zJmzBiKJtXTAUQ8e8a/ixdzbvZsQmIn8ZpbW1Phiy+o/M03ZE1qqXJkpFZjaM4ciL03AI0aaUNDLVtquVoSPDcIbq+EG39A8I247fbNocQQ7U9Z5izEW5Pg5Q0ZGRlhb29P7ty5iYodJxfvn5mZmfS4ZGDR0dFs2LCBqVOncvXqVUDr6Rw8eDCjRo0ib968SZ4f4ufHuTlz8Fy4kMhnzwDI4uBAlREjKD9oEJkS6ykBePwYFi2C+fPjEspZWGg1hoYOTXrVUNAVuDZPqzMUHbviyMxaSyRX/CuwLv7an4EQInESvLwlExMT+eEpRCoLDw9n5cqVTJs2DS8vL0CrqzVkyBCGDRtGrly5kjz/ybVr/DNjBldWryYmMhIAW0dHqo0ZQ+mePTFJahLv7dvw+++wfDnEzqfCwUFbNTRoEOTMmfB5SmlDQ9d+j83NEsumDJT4Bgr1BLOkM/gKId6MBC9CiDQXEhLC4sWLmTlzJn5+fgDkypWLESNG8NVXX2GT1MofwO/0ac789hs3XFy0YAJwqF2b6t99R1EnJ4ySyqT7zz8wfTps3gw6nbatYkX49lvo3BkSmxAeHab1sPw3G55djd1oBPnaQslhkLt+4iuOhBApIsGLECLNBAUFMXfuXGbPns2TJ08AyJcvH99++y0DBgzAysoq0XOVUnjt2sWZ337D58gR/faibdpQfcwY8tapk/iNdTrYuRNmzICXzqV5cxg9Gj79NPHAI8wPrv+h1RqKeKxtM80CRT+HkkO1ybhCiHdKghchxHsXFBTEnDlz+P333wmMXXJcrFgxxo4dS+/evZPM0aKLieH6xo2cmjKFRxcvAmBsZkbpXr2oOno0OUuXTvzG4eGwdi3MnBmXn8XMDHr0gFGjkp7P8uQc/Pc73PsLdLHz3TIX0gKWIv3BPOneISFE6pHgRQjx3rzoaZk1a5Y+aCldujQ//PADXbp0SXIeWUxUFFfXr+f0lCk8vX4dAPOsWSk/eDBVhg1LeuXQkyewcCHMmwcPHmjbrK3hiy+0SbiJTQDWxYCvqzafJeClHppcn4DjCMjbBozl26gQ75v8rxNCvHPPnj3TBy1Pnz4FoFSpUkycOJFOnTolGbRER0RwedUqzvz6K0Gxk3gtbG2pMnw4lb75Bouk0v/fuQOzZsGyZXGTcPPnh+HDYcCAxPOzRAXDreVwfa6WXA7AyBQKdNGClhxV3+wDEEKkKglehBDvTEqClqiwMC4uW8Y/v/1GsI8PAFa5c1N11Cgqfvkl5lmzJn7jq1dh6lRYvz4uqVzFitp8li5dEp+EG+YP1+bCjYVxWXDNbaHYYC0LrlXSS7SFEO+HBC9CiFQXEhLCvHnzmD59ukHQMmHCBDp37pxk0BL5/Dn/LlrE2RkzeB6bZyWLgwPVxoyh/MCBmCUxiZezZ7WgZetW/aojmjSBMWO05HKJTcJ9dg2uzgSvVaDTlliTtYTWy1L4MzBN4p5CiPdOghchRKoJDw9n8eLFTJkyRZ+F2tHRkYkTJyYbtEQ8e4bnH39wdtYsfXVn64IFqT52LGX79sXUwiLhE5XSVgxNmaJVeH6hQwcYNw6qJjHE8+gUXJkGPi7oU/fnrAWlxkC+NmCUxBJrIUSakeBFCJFiUVFRrFq1ip9++gmf2CGeokWL8tNPP9GtW7ckg5bwp0/xmDOHc3PmEBE7iTdb0aLUGD+e0r17Y5LYEI9S4OamBS0nTmjbTEy0lUNjx0Jiq46UDnx3wtXp8PBo3Pa8baD0GMiVxBJrIcQHQYIXIcRb0+l0ODs7M3HiRG7evAloeVomTJhA3759k6z4Hf70KWd//51zs2cTGRwMgG2pUtT8/nscu3ZNvLpzTIyWUG7KFPj3X21bpkxaZedvv4XChRM5LwLurNeClhdJ5YzNoFBvKDUabEq91WcghHj/JHgRQrwxpRTbt2/nhx9+4NKlS4CWEXf8+PF88cUXWCQ2xANEBAXhMWcOHrNmEREUpJ1bvjw1f/iBEh07Jp4NNzJSy9Hy669wI7bYYZYs8OWXMGIE2Nsncl4Q3FwC12ZDmFacETNrKP4llBgKVg5v8xEIIdKQBC9CiDdy+PBhxowZw5kzZwCt9tC3337LsGHDyJIl8Ro+kcHBnJ8/n3+mTyc8dhJvzrJlqf3TTxRv1y7xoCU8HP78E377DWKHpMieXavs/M03YGub8HlhD7T8LDcWQpRWnBFLB20SbrFBWgAjhEiXJHgRQryWixcvMnbsWNzc3ADInDkzw4YNY/To0WTPnj3R8yKfP8dzwQL+mTZNPxHXtlQpak+aRMlOnZIPWqZOhfuxPSb29lom3EGDILGl0qE+cGU63FoCMeHaNpvSUOpbKNgDTJIoziiESBckeBFCJOnevXtMmDCB1atXo5TC1NSUQYMGMWHCBPLkyZPoeVFhYfy7aBFnfv2V0NiVR9mLF6fWxIk4duuGcWKTeBMKWvLn11YO9esHiQ1JBd+CK7/GLneOTd+fowaUGQ95nWTlkBAZiAQvQogEPXnyhKlTpzJv3jwiIiIA6Ny5M5MnT6Z48eKJnhcTFcWlFSs4+fPPhPj6AmBTpAi1JkygdM+eiU/EDQ/XMuFOnQqx55EvH4wfr03GzZQp4fMCL8OVqXB3g7aSCCB3Ayj7A+RJosCiECLdkuBFCGEgLCyMefPmMXXqVH39ofr16zNt2jSqV6+e6HlKp+Pa339zfMIEnsZOqM2aPz+1JkygTJ8+iS95jojQgpYpU94saHlyDi5PBu8tcdvsW0DZ72W5sxAZnAQvQghAW/a8bt06xo8fr8/VUrZsWX777TdatGiBUSI9GEopvHbv5tj48QR4egJgmSsXNb//ngpffIFpYsHHi6Bl6tS4ibh582pBy+efJx60BBzTgha/3XHb8nfUhodsK79N04UQ6YwEL0IIjh49ysiRIzl79iwA+fPn53//+x+9evVKMsGc7/HjHB03Dp+jWrI386xZqTp6NFVHjEi89lBEBCxfrvW0vG7QohT479OCloDD2jYjY20Cbplx2oRcIcRHQ4IXIT5it2/fZsyYMWzevBmArFmz8v333zN06FAsLS0TPe/hhQscHT+e2zt3AmCSKROVhgyh+tixWOXMmfBJUVFa0DJ5Mnh7a9scHOKCloQm4ioFfu5wcRI8PqVtMzaDwn2h9HeQtehbtlwIkZ5J8CLERygwMJDJkyczd+5cIiMjMTY2ZuDAgfz0009JriB65u3N8R9/5PLq1aAURiYmlOvfn1oTJpA1X76ET4qJgQ0bYOJEuH1b2+bgoK0eGjAg8aDlwX64MAEendS2mVhA0UFaNtzM+VP4CQgh0jMJXoT4iERHR7NkyRImTpzIo9icK02bNmXmzJmULVs20fPCAwM58+uveMyeTUzsyqOSXbpQ53//w7ZEiYRPUgq2bYMffoDLl7VtuXPD999reVoSW/L84KAWtDw8pr03sYBiX2p1hyzt3qrdQoiMRYIXIT4Su3fvZuTIkVy9qtX1KVWqFDNnzqR58+aJTsaNiYzEc+FCTv3vf4Q9fgxAvnr1qD99OvZJrDxi3z4tSInNwku2bPDdd1pG3MyZEz7nwWG4ODFuTotxJig2GMqMBctEUv8LIT5KErwIkcHdunWLkSNHsn37dgBy5MjBzz//zMCBAxMtnKiU4trff3N0/HiCYod6bEuVov60aRRp1SrRYIdTp7Sg5cAB7b2VlVZ3aPRoLYBJSMAxLWh5EHuOsTkUHahNxLXK+7bNFkJkYBK8CJFBPX/+nKlTpzJjxgwiIiIwNTVl6NCh/Pjjj2RLLJAAvI8c4fDo0fj/8w8Ame3sqPPzz5Tt1y/xBHOXLmlBS2yAhLk5fPGFNhk3sTk0D09qQYv/Xu29sRkUHQClx8mcFiFEkiR4ESKDUUrx999/M3r0aH2+liZNmjBnzhxKlSqV6HmBt25x+NtvubF1KwBmmTNTbcwYqo4ahXliQz0+PjBhAqxaBTodGBtD377atoIFEz7n0WktaPHbo703MoWi/bU8LZkTOUcIIV4iwYsQGcjFixcZOnQohw4dAqBQoULMmjWLdu3aJTrUExEUxKnJk/GYPRtdVBRGJiaUHziQ2pMmkTmxXpPAQPj1V5gzR0vrD9CxI/zyCzg6JnzOU0/49we4ry2vxsgEivSFMt9DlsJv22QhxEfog65UduTIEVq3bo2DgwNGRka4uLgke87hw4epUqUKFhYWFClShEWLFr37BxUijT19+pShQ4dSqVIlDh06hIWFBT/99BNXrlyhffv2CQYuupgY/l28mD+LF+ef6dPRRUVRqGlT+vz7L00WLkw4cImIgFmzoGhR+O03LXCpWxdOnoRNmxIOXIJvwvHusKuSFri8CFqcrkGNPyVwEUK8sQ+65+X58+dUqFCBfv360bFjx2SP9/LyomXLlgwcOJC1a9dy/PhxvvrqK3LlyvVa5wuR3iilWL16Nd9++y0PHz4EoGPHjsycOZOCiQ3bAHf37+fgiBE8ungRANuSJWkwaxaFEysDoNPBunXw449w9662rXRprffFySnh4oehvnDpZ7i1DFSMtq1gNyj3E1gnsrxaCCFewwcdvLRo0YIWLVq89vGLFi2iQIECzJ49G9CWgp49e5YZM2ZI8CIynCtXrvDVV19x+LC2tLhUqVLMnTuXxo0bJ3rO0xs3ODR6NLdiJ9ZaZM9O7UmTqPDll4kXTty7F779Fv79V3vv4AA//wx9+kBCE3gjHsOVX+H6fIiJHVJyaAkVJkP2im/bXCGE0Pugg5c3dfLkSZo2bWqwrVmzZixbtoyoqKgEl4VGREQQEZt0CyA4OPidP6cQKREaGsovv/zC9OnTiY6OxsrKiokTJzJixIhElz5HPHvGyf/9j3Nz5ujntVT6+mtqTZyIpa1twje6elVb4uzmpr23toaxY2HYMG0J9KuiQuC/3+G/GRD1TNuW6xOoMAVy102FlgshhCZDBS/+/v7xUpvnyZOH6OhoHj16hL19/ERXU6dO5aeffnpfjyhEiuzcuZMhQ4Zw584dAFq3bs28efMSHSJSSnF1/XoOf/stz/38ACjcsiUNZswgR2Irjx49gkmTYNEiLbW/qSl89ZU2ZJRQ3aKYSLi5RBsiitCGrshWQQtaHFokPKQkhBApkKGCFyDeeL1SKsHtL4wbN46RI0fq3/v6+lK6tFSoFR8WHx8fhg0bxpYtWwCt6vO8efNo27Ztouc8vHCB/UOG6Cs+ZytWjE/nzKFIy5YJnxARAfPmaSuGgoK0bW3bwrRpkFAJAKXAezN4joOQm9q2LMWg/P+gYBet6rMQQrwDGSp4sbOzw9/f32BbQEAApqam5MiRI8FzMmXKRKZMmfTvnz179k6fUYg3ER0dzdy5c5k4cSIhISGYmJgwcuRIJkyYQJYsWRI8JzwwkOMTJuD5xx8onQ5TS0tq/vADVUeNwvSlf+t6SsGWLTBmTFzhxIoVYeZM+PTThB8s4Cic/xYen9beW+SGcpO0JHPGicydEUKIVJKhgpdatWrh6upqsM3d3Z2qVasmOhdAiA/V+fPn+fzzzzl//jwAtWvXZtGiRZQrVy7B45VOx6WVKzkydixhsSuPSnTuTIMZM7AuUCDhm5w9CyNHQmzvDHZ2MGUKfPYZmJjEPz7oKniOBd/YTLomVlDqWyg1Csyypqi9Qgjxuj7o4CUkJISbN2/q33t5eeHp6YmtrS0FChRg3Lhx+Pr6snr1agC++OIL5s+fz8iRIxk4cCAnT55k2bJlbNiwIa2aIMQbCwsL4+eff2b69OnExMSQPXt2pk+fTr9+/TA2Tngoxt/Dg/1ff43faa0nxLZUKRrNm0fBRo0SvsmDB1rq/uXLtfeWltrk3DFjIKEenTA/uDgJbv0JSqflaik6AMpNlKKJQoj3T33ADh48qIB4rz59+iillOrTp4+qX7++wTmHDh1SlSpVUubm5qpQoUJq4cKFb3RPb29vBShvb+9UaoUQr+/w4cOqePHi+n/rnTt3Vv7+/okeHx4YqPZ+/bWabmSkpoOanSWLOjNjhoqOiEj4hMhIpWbNUsraWiltwEipnj2Vuncv4eOjQpT6d4JSzlZKrUN7HW6nVODVVGitECKjeV8/Q42Uip3RKgBtYmT+/Pnx9vYmX758af044iPx7NkzvvvuO31GaAcHBxYsWJDohFylFNc2buTg8OH6VUSlevSg/vTpZHFwSPgme/dqy5yvXtXeV6kCc+dC7doJ3EAHXmvg3/EQdl/blqMmVJoOuT9JUVuFEBnX+/oZ+kEPGwnxMdixYwdffvmlvojiwIEDmTZtWqKVnwNv32bf119zZ/duALIXL07jhQsTHyLy8oJRoyC24CI5c8LUqdCvX8LzWgKOgMcIeHpOe5+5EFSaBvk7ybJnIcQHQYIXIdJIQEAAw4YNw9nZGYCiRYuydOlSGjZsmODxMZGRnJ05k5M//0x0eDgm5uZUHzeOGmPHYmphEf+E0FAtff+0adoyaBMTGDJEy+GSUGAUfAs8x4C3thwb06xQ9gcoORRMEri+EEKkEQlehEgDmzZt4ssvv+TRo0cYGxszatQoJk2ahFVCmWsBn2PH2PvFFzy+fBmA/A0b0mThQmxLlox/8IulzyNGgLe3tq1RI60CdJky8Y+PDILLv8C1uaCL1PKzFB0I5X/WlkALIcQHRoIXId6jJ0+eMGTIEP0KuPLly7Ns2TKqVq2a4PFhT55wZMwYLi5bBoBlzpw0mDWL0r16JZx48dYtrXcldkiJQoW0fC3t28cf8tHFwK2lcOFHiHikbbNrApVnQraEl2MLIcSHIFWCl6ioKPz9/QkNDSVXrlzYJlYrRYiP2M6dOxkwYAD+/v6YmJgwbtw4fvzxR8zNzRM8/tqmTez/+mtCAwIAKDdgAPV++y3hWkTh4drw0JQp2hCRublWh2jsWG0Z9KsCjsLZbyAwttiidUmoNFMroCjzWoQQH7i3Dl5CQkJYt24dGzZs4MyZMwbFDfPly0fTpk0ZNGgQ1apVS5UHFSK9evbsGSNHjmRZbO+Jo6Mjq1atonr16gke/9zfn31ff82N2FIAto6ONF26lHyfJLLKZ+9e+PpruHFDe9+4MfzxR8Ip/UN94fwYuLtee2+WDcr/BMW/lMy4Qoh0462Cl99//53JkydTqFAh2rRpw9ixY8mbNy+WlpY8efKES5cucfToUZo0aULNmjWZN28exYsXT+1nF+KDd+DAAfr168e9e/cwMjJixIgR/PLLL1gm0BuilOLy6tUcGjGC8KdPMTY1pfrYsdT84YeE0/rfv6/Na/n7b+29vT3Mng2dO8fvPYmJ0Co+X/4Fop8DRlBsIJT/BSxypXq7hRDiXXqr4OXEiRMcPHgw0TTl1atXp3///ixatIhly5Zx+PBhCV7ER+X58+eMHTuW+fPnA1C4cGFWrlxJvXr1Ejw+6O5d9g4ezJ09ewDIU7kyzZYtI3fFivEPjo6G+fNhwgQIDgZjYxg6FH76Cayt4x/v6wbnhkNwbM9MzlpQdR7YVkmFlgohxPv3VsHLxo0bX+u4TJky8dVXX73NLYRIt06fPk2vXr30pS2++OILpk+fnmAhRaXT4blwIUfGjiUqJASTTJmo/dNPVBs1CmPTBP57enjAgAHg6am9r1kTFi7UCim+Kvimlq/l/g7tvYWdlq+lUE+p+CyESNdktZEQqSQmJoZff/2ViRMnEhMTQ968eVm2bBnNmjVL8Pgn16/jPmAAPrFFEfN+8gnN/vwz4eXPz5/DxInw+++g00H27PDbb/D551rPy8uiw+DyFLg6LXbpsyk4DoeyP4JZAj0zQgiRzrxx8PL06VOUUtja2vLw4UOOHDlCyZIlKVu27Lt4PiHSBW9vb3r16sWRI0cA6NKlC4sWLSJ79uzxjlU6HefmzePo2LFEh4djliUL9X79lYpffolRQoUX3d3hiy+0TLkA3btrc1tyJ5CD5f4uODsEQm5r7+2bQeXZYOOYOg0VQogPwBv1Hf/5559UrVqVKlWqsHDhQtq3b8/+/fvp1q0bS5YseVfPKMQHbdOmTVSoUIEjR46QOXNmVqxYgbOzc4KBS9Ddu/zduDEHhw8nOjycgk2a0PfSJSp9/XX8wOXRI/jsM2jWTAtcChSAnTth/fr4gUuoDxztBIdaaoGLVT6ouxka7JLARQiR4bxRz8u8efO4fPkyoaGhFChQAC8vL3LlysWzZ8+oV68egwYNelfPKcQH5/nz5wwbNky/BLpatWqsW7cuwcnpSikurVjBweHDiQwOxtTKigYzZ1Jh8OD4yeaU0gKU4cO1AMbISJuQ+8sv8Oq8GV2Ulhn34kRtFZGRCZQcDuUmgVn8OTZCCJERvFHwYmJigoWFBRYWFhQrVoxcubQlltbW1gln+xQig/Lw8KBHjx5cv34dIyMjxo4dy08//YSZWfxcKc/9/dkzcCC3d2gTZ/PWqUOLVavIVrRo/AvfuQNffhmXIbdcOVi6FGrUiH/swxPwzxcQeFF7n7M2VFsI2cunUiuFEOLD9EbBi6mpKeHh4VhYWHD48GH99uDg4FR/MCE+RDqdjpkzZ/L9998TFRVF3rx5Wbt2LQ0aNEjw+GsbN7Lvyy8Je/wYE3Nz6vzyC1VHjsT41WrOOp22aui777TJuZkyaUuhv/0WXg2IIh6D53dwS+vxwdxWW0VUpJ+sIhJCfBTeKHg5cOAAmWKTZdnY2Oi3h4WF6bvOhcioAgIC6NWrF3v37gWgQ4cOLF26NMFyGGFPnrB/yBD+i61hlLtSJVqsXk2uhCa2e3lB//5w6JD2vl49WLIEXl11pBTcWa/lbHlRi6hIf6j4G1jkTKVWCiHEh++NgpeE8lQA5M6dm9wJrXwQIoM4duwYXbt25f79+1haWjJ37lw+//zzBIdL7+7bx64+fQi5fx8jExNqjB9PrR9+wOTVGkY6HSxerPWuPH8OVlba8uevvoq//Pn5XTjzBfjFDifZlIFqiyB3IiUDhBAiA0uVPC/h4eFcuHCBgIAAdDqdwb42bdqkxi2ESBNKKWbMmMG4ceOIiYmhVKlSbNy4kTJlysQ7NiYykqPff8/ZGTMAsC1ZkharV2OfUA2jO3e0HC0HDmjv69WD5cvh1Xkwuhi4Ph8ufK9NyDU2h7IToNS3YJJwQUchhMjoUhy87N69m88++4xHjx7F22dkZERMTExKbyFEmnj69Cl9+/Zl+/btAPTo0YPFixcn2AP55Pp1dvbowQMPDwAqfPklDWbMwMzKyvBApbQhodGjISREq/j8668wZEj83pbAi3B6ADw+o73PVReqL5Glz0KIj16KZ/cNGTKEzp074+fnh06nM3hJ4CLSq7Nnz1K5cmW2b9+Oubk5ixYtYu3atfECF6UUF1esYE3lyjzw8MDC1pZ2Li40WbAgfuBy9y40baolnAsJgU8+gQsXtGXQLwcuMeHw74+wq7IWuJhZa0NEjQ9J4CKEEKRCz0tAQAAjR44kT548qfE8QqQppRQLFy5kxIgRREZGUqRIETZu3EjlypXjHRseGMjewYO5FlvVOX/DhrRcs4asefO+elFtSGjECK2QoqUlTJkSP2gBCDgKZwbCs2va+3ztoOp8sHrlmkII8RFLcfDSqVMnDh06RNGEclYIkY4EBwczaNAgnJ2dAWjfvj3Lly8nW7Zs8Y71PX6cHT16EHzvHsamptT53/+o9u238ZdAP3oEAweCi4v2vnZtWLECSpQwPC76OXiO1ea3gFZEsdofkL9D6jZSCCEygBQHL/Pnz6dz584cPXqUcuXKxUvSNXTo0JTeQoh37urVq7Rv355r165hamrKb7/9xogRI+KtJtJFR3Nq8mRO/vwzSqfDpkgRnDZsSHhS7u7d0K8f+PtruVp++QVGjYJXA5yAI3CqX1w9oqIDoNJ0MM/2bhorhBDpXIqDl/Xr17Nnzx4sLS05dOiQwTd7IyMjCV7EB8/V1ZWePXsSHBxMvnz5+Ouvv6hdu3a840Lu32dHt276KtCle/em0fz5ZLJ+pVJzWBiMGQPzY3tRSpXS0v1XrGh4XPRz8BwP1+dq763yQ41lYN8klVsohBAZS4qDlx9++IGff/6ZsWPHYpxQRVwhPlA6nY4pU6YwYcIElFLUq1ePjRs3Jpiz6N6BA+zo3p3QgADMs2al8cKFlO7ZM/5Fz5+Hnj3h6lXt/TffaLlbLC0Njws4Gtvbckt7X3QAVJ6pTc4VQgiRpBQHL5GRkXTt2lUCF5GuhISE0KdPH7Zs2QLA119/ze+//x5v2FPpdJyaMoUTEyeidDpylS9P640bsX11zkpMDMyYAT/+CFFRYGenzW1p3tzwuOhQ+Pd7uDYHUFr15+p/gkOzd9haIYTIWFIccfTp04e//vorNZ5FiPfi9u3b1KpViy1btmBmZsbSpUuZP39+vMAl9NEjNrdqxfEff0TpdJTt148eJ0/GD1zu3oVPP4WxY7XApV07uHgxfuDy8ATsqgjXZgNKS+3f8pIELkII8YZS3PMSExPDtGnT2LNnD+XLl4/3A2DWrFkpvYUQqWbfvn106dKFp0+fYmdnx5YtW6hVq1a84+6fOoVrly4Ee3tjamFBowULKNevX/wLbt6sZcoNCoLMmWHuXG2S7ssTfWMi4MKPcHUGoMDSAaovhbwt311DhRAiA0tx8HLx4kUqVaoEwKVLl1L8QEK8C0opfv/9d7799lt0Oh3Vq1dny5Yt5H0lJ4tSinNz53J49Gh00dFkL16cNps2kat8ecMLhodrK4cWLNDe16gB69bFT+8feBlO9ITAf7X3hftAld/BPPs7aqkQQmR8KQ5eDh48mBrPIcQ7ExYWxqBBg1i7di0Affv2ZeHChVhYWBgcF/HsGbv79+fG5s0AlOjcmWZ//hl/NdH169C1K3h6au/HjNGWQb/c66h0Ws6W82NAFwGZckKNPyFf23fVTCGE+GikOHiZOnUqefLkoX///gbbly9fzsOHD/nuu+9Segsh3lpAQABt27bl1KlTmJiYMGvWLL755pt4+VseXb7MtvbteXrjBsZmZjSYMYNKCRzH+vUweLCW3j9nTli9Glq0MDwm9D6c7g9+e7T39s2h5gqwtHuHLRVCiI9HiifsLl68GEfH+PVWypQpw6JFi1J6eSHe2rVr16hVqxanTp0ie/bsuLu7M3To0HgByc1t21hXsyZPb9wga/78dD96lMqvHhcaCgMGaMugQ0K0KtCenvEDF+8tsKu8FriYWGip/Ru4SeAihBCpKMU9L/7+/tjb28fbnitXLvz8/FJ6eSHeypEjR2jXrh1Pnz6lSJEiuLm5UbJkSYNjlE7HyV9+4cTEiQDkb9CA1hs3YpUzp+HFLl/WhokuX9Ym4v74o/Yyfem/T1QweAyH28u199krQe11YFPqHbZSCCE+TikOXvLnz8/x48cpXLiwwfbjx4/j4OCQ0ssL8cbWrVtH//79iYyMpGbNmmzfvp1cuXIZHBMZEsKuPn24EZvnpdI339Bg5kxMDOatKC1Xy5AhWtZcOzttUu6nnxre8NEpONErNuGcEZQeA+V+BhPzd9xSIYT4OKU4eBkwYADDhw8nKiqKT2O/qe/fv58xY8YwatSoFD+gEK9LKcXkyZP58ccfAejYsSNr1qzB8pXstoG3b+PSti2PLl3CxNycxgsXUu6VOVuEhcFXX8HKldr7Jk1gzRp4uXq60sGV37Rl0CoGrApArdWQp/47bKUQQogUBy9jxozhyZMnfPXVV0RGRgJgYWHBd999x7hx41L8gEK8jqioKAYPHsyKFSsA+Pbbb/n111/jZX6+u38/rl26EP7kCZnt7Gi7ZQsOr+Z5uXMHOnTQUv0bG2srib77Tvv7C2EP4ORn4O+uvS/YDaotlGKKQgjxHhgppVRqXCgkJISrV69iaWlJ8eLFyZQpU2pc9r3z8fEhf/78eHt7ky9fvrR+HPEaAgMD6dSpE/v378fY2Jj58+fz5ZdfGhyjlOLcnDkcGj0aFRODXbVqtN26layv5HnB3R26d4cnT7TVRH/9FX+YyP+Alrsl3B9MLLVJuUVeSUwnhBAfoff1M/Ste17Gjx9Pu3btqF69OgBZsmShWrVqqfZgQryOu3fv0qpVKy5fvkzmzJn5+++/adnSMHNtdHg4e7/4gsurVgFQ+rPPaLp4MaYv53nR6WDqVG0irlJQrZqWPTd//peOiYFLP8Ol/wEKbEpDnb8hW5n30FIhhBAvvHXw4ufnh5OTEyYmJrRu3Zq2bdvSuHHjdNvjItKfCxcu0Lx5c/z8/HBwcGDHjh36bM8vhD56hEvbttw/cQIjY2MazJxJ5WHDDJdBBwVBnz6wbZv2fuBALc3/y8FNqK/W2xJwWHtf9HOoMhdMrd5xK4UQQrzqrYOXFStWoJTi2LFjuLq6MmrUKHx9fWnSpAlt2rTBycmJnK8uORUilRw9epTWrVsTFBRE2bJlcXNzI//LvSTA0xs32NyyJYE3b5LJxobWGzdSqEkTwwtdvgzt28ONG2BuDn/8oeVzedn9Xdr8lohHYJoFqi+GQj3ecQuFEEIkJkVJ6oyMjKhbty7Tpk3jv//+48yZM9SsWZOlS5eSN29e6tWrx4wZM/D19X3reyxYsIDChQtjYWFBlSpVOHr0aJLHr1u3jgoVKmBlZYW9vT39+vXj8ePHb31/8eHZvn07TZs2JSgoiE8++YSjR4/GC1x8jh1jfa1aBN68iXWhQvQ4eTJ+4PL331pNohs3tOGhY8cMAxddFJz/Dg611AKX7BWh+TkJXIQQIq2pd+TBgwfqzz//VG3atFHTp09/q2s4OzsrMzMztXTpUnXlyhU1bNgwlTlzZnX37t0Ejz969KgyNjZWc+bMUbdv31ZHjx5VZcqUUe3atXvte3p7eytAeXt7v9Uzi3dr+fLlysTERAGqdevWKjQ0NN4xV52d1SxzczUd1Jpq1VSIv7/hAdHRSo0erZQ2u0WpRo2UCggwPCb0vlLunyi1Du31zxClosPeYcuEECL9e18/Q99Z8JIaqlevrr744guDbY6Ojmrs2LEJHj99+nRVpEgRg21z585V+fLle+17SvDy4Zo+fboCFKD69u2roqKiDPbrdDp1aupUNR3UdFBb27VTkc+fG14kKEipVq3iApcxY5R65Toq4JhSW+y1oOVva6XubX7HLRNCiIzhff0MTXGel5EjRya43cjICAsLC4oXL06bNm2wtbV9o+tGRkbi4eHB2LFjDbY3bdqUEydOJHhO7dq1+f7773Fzc6NFixYEBASwadMmWrVqleh9IiIiiIiI0L8PDg5+o+cU755SiokTJ/K///0P0HK4/PbbbwaTbmOiotj/9ddcWLoUgCrDh1N/xgyMTUziLnTnDrRuDZcuaZNxV67U0v7H3UirBH1uJKhosCkDdbeAdYn30EohhBCvLaXRT4MGDZS1tbXKnDmzqly5sqpUqZLKkiWLsrGxUTVq1FDZsmVT2bNnV5cvX36j6/r6+ipAHT9+3GD75MmTVYkSJRI9b+PGjSpLlizK1NRUAapNmzYqMjIy0eMnTpyo/23+5Zf0vHwYdDqdGj58uP7rMmXKlHjHhAcFqY3NmqnpoGYYGyuPuXPjX+jYMaVy5dJ6W+ztlTpzxnB/1HOljveMGyY62lWpyOB31CohhMiY3lfPS4qrSr9YIn3//n08PDw4d+6cftVR9+7d8fX1pV69eowYMeKtrv9qBWClVLxtL1y5coWhQ4cyYcIEPDw82L17N15eXnzxxReJXn/cuHEEBQXpX1euXHmr5xSpLyYmhkGDBjF79mwA5s2bFy9rc7CPD85163Jnzx5Mraxou3Urlb/5xvBCq1drieYePoRKleDMGS2Pi/4it8C9FtxZB0YmUPl3qLMBzLK84xYKIYR4KymNfhwcHBLsVbl06ZJycHBQSinl4eGhcuTI8UbXjYiIUCYmJmrLli0G24cOHarq1auX4Dm9evVSnTp1Mth29OhRBaj79++/1n1lzsuHITIyUnXv3l0BytjYWK1YsSLeMQ8vXVIL8+ZV00H9kSeP8vvnH8MDYmKUGjs2bn5Lhw5KhYQYHuOzQ6m/s2m9LZvzKPXg8LtrlBBCZHDppuclKCiIgICAeNsfPnzIs2fPAMiWLZu+7tHrMjc3p0qVKuzdu9dg+969e6ldu3aC54SGhsarZWMSO+dBpU4VBPEeREZG0q1bNzZs2ICpqSnOzs707dvX4Jj7J0/iXLcuIb6+2JYqRc/Tp7GrWjXugJAQ6NgRfv1Ve//997BxI2TOrL1XOrgwCQ47QVQg5KwFzT0gd7330EIhhBApkeIJu23btqV///7MnDmTatWqYWRkxJkzZxg9ejTt2rUD4MyZM5Qo8eaTHkeOHEnv3r2pWrUqtWrVYsmSJdy7d08/DDRu3Dh8fX1ZvXo1AK1bt2bgwIEsXLiQZs2a4efnx/Dhw6levToODg4pbap4DyIjI+natSsuLi6Ym5uzZcuWeBOub7u5sb1TJ6LDwnCoVYv2O3Zg+fKE8Hv3oE0b+PdfyJQJli2Dnj3j9kc9g+M94P5O7X3xr6HyLDAxfw8tFEIIkWIp7boJDg5WAwYMUObm5srY2FgZGxsrc3NzNXDgQBUS20V//vx5df78+be6/h9//KEKFiyozM3NVeXKldXhw3Hd+n369FH169c3OH7u3LmqdOnSytLSUtnb26uePXsqHx+f176fDBulnfDwcNWmTRsFqEyZMqldu3bFO+bS6tVqhomJmg5qU4sWKuLVYaB//lEqTx5tmCh3bqVOnDDcH3xLqR2ltWEiZwulbq9+hy0SQoiPy/v6GZqqVaVv376NUoqiRYuSJUv6nOwoVaXTRkREBB07dmTnzp1YWFiwbds2mjZtanCMx5w5HBw+HIDSvXrRbPlyTMzM4g7Ys0cbKnr+HMqXh+3boWDBuP0BR+FoBy1brqUD1N8OtlXeQ+uEEOLj8L5+hr7VnJd79+7F25YlSxbKly9PhQoV4gUuKSkPIDK+iIgIOnTowM6dO7G0tMTV1dUgcFFKcfJ//9MHLlWGD6fFqlWGgcuaNeDkpAUujRvD0aOGgcut5XCgkRa42FaFZv9I4CKEEOnUWwUv1apVY+DAgZw5cybRY4KCgli6dClly5Zly5Ytb/2AImOLioqiS5cuuLm5YWlpyY4dO2jcuLF+v1KKI999x/EJEwCo8/PPNJg1C6MXE7OVgt9+g88+g+ho6NEDdu4Ea2ttvy4Gzo2G059rtYoKdIHGh8FK5kAJIUR69VYTdq9evcqUKVNo3rw5ZmZmVK1aFQcHBywsLHj69ClXrlzh8uXLVK1alenTp9OiRYvUfm6RAURHR9OzZ0+2b9+OhYUFrq6ufPrpp/r9Sqdj39df8++iRQA0/P13qsT2vgAQEwMjRsC8edr70aO1QOZFYPPqxNxyk6DsBEgkT5AQQoj0IUVzXsLDw3Fzc+Po0aPcuXOHsLAwcubMSaVKlWjWrBlly5ZNzWd9L2TOy/sRExNDnz59WLduHebm5mzbto3mzZvr9+uio9ndvz9X1qwBIyOaLllC+ZcrPoeHQ+/esGmT9v733+HlwCbECw63hqDLYGIBNVdBwS7vp3FCCPGRel8/Q1O0VNrCwoIOHTrQoUOH1Hoe8RHQ6XQMGjSIdevWYWpqyt9//20QuMRERbGzZ0+ub9yIsakpLdeswbFbt7gLBAZCu3Zw+DCYmWkZdF/ebzAx1x7qbYMcL2XUFUIIka6lOM+LEG9CKcXQoUNZvnw5xsbGrF+/nrZt2+r3x0RGsqN7d25s2YKxmRltNm6k2Ev78fWF5s214opZs4KLi5b6/4U76+FUX21+i20VLXCxyvve2ieEEOLdk+BFvFcTJkzgjz/+wMjIiFWrVtG5c2f9vpjISFy7dOHmtm2YmJvTdutWirRsGXfy1avQrBl4e4O9PezaBRUqaPuUgqszwHOM9j5/J6i1Ckyt3mPrhBBCvA8SvIj3Zvbs2fzyyy8ALFiwgF69eun3RUdE4Nq5M7dcXTHJlIl2Li4UfmkoCQ8PLXB5/BhKloTdu6FQIW2fLgbOjYTrc7X3JUdA5RlglOLqF0IIIT5AEryI92LVqlX6yuKTJ082qPQdHRHB9o4dub1zJ6YWFrTbto1CLyeoO3oUWrWC4GCtGrSbG+TMqe2LCYcTvcE7duJupZlQauT7apYQQog0IMGLeOe2bdvG559/Dmj1qsaNG6ffFxMVxY6uXbXAxdKS9q6uFGzUKO7k3buhQwcIC4P69cHVVZvrAhD5FI60g4AjYGyurSgq9NLEXSGEEBlSivvV+/bty5EjR1LjWUQGdPjwYbp27UpMTAx9+/ZlxowZGMXmWdFFR7OzRw9ubtuGqYVF/MBl82atwGJYGLRsqc1xeRG4PPeGvZ9ogYuZNTTcLYGLEEJ8JFIcvAQHB9O0aVOKFy/OlClTpBSA0Lt06RJt27YlIiKCdu3asXTp0rjAJSaGXX36cH3TJv3kXIPAZd066NIFoqKgc2fYuhUsLbV9gRfBvRYEXdFqFDU5BnkapkELhRBCpIUUBy+bN2/G19eXIUOGsHHjRgoVKkSLFi3YtGkTUVFRqfGMIh3y9fWlRYsWBAUF8cknn7BhwwZMTbVRSqXTsWfAAK6uX4+xqSltNm0ynJy7cqWWgE6ng/79YcMGMDfX9j04pPW4hPmCTWloehKylXvv7RNCCJF2UmU5Ro4cORg2bBjnz5/nzJkzFCtWjN69e+Pg4MCIESO4ceNGatxGpBNBQUG0aNECHx8fHB0d2bZtGxYWFoCW52Xf119zeeVKjExMcHJ2pmjr1nEn//mnFrAoBV98AUuXgomJtu/eRjjYTEv7n6uu1uOSuUAatFAIIURaStW1pH5+fri7u+Pu7o6JiQktW7bk8uXLlC5dmt9//z01byU+UJGRkXTo0IGLFy9iZ2fHrl27sLW11e8/On68VqvIyIiWq1dTomPHuJMXLYKBA7XAZcgQWLAgrk7RzT/hWFfQRWo5XD51B/Ps77l1QgghPgQpDl6ioqLYvHkzTk5OFCxYkI0bNzJixAj8/PxYtWoV7u7urFmzhp9//jk1nld8wJRS9O/fnwMHDpAlSxbc3Nwo9CIXC3Bm2jTO/PorAE0XL6ZUjx5xJ8+fD19+qf19xAiYOzeugOLVmXBmIKCg2GCo46zVKxJCCPFRSvFSaXt7e3Q6Hd27d+fMmTNUrFgx3jHNmjUjW7ZsKb2V+MD99NNPrFu3DhMTEzZt2kSlSpX0+y4sXcqR774DoN60aZQfODDuxD/+gG++0f4+Zgz8+qsWuCgFFybAZS2xHaXGQMVfpSq0EEJ85FIcvAwbNoxRo0ZhZWWYhl0phbe3NwUKFCB79ux4eXml9FbiA7ZhwwZ++uknABYtWkSzZs30+/77+2/cBw8GoPrYsVT/9tu4E5cs0YaIAMaOhSlTYgMXHXgMh+vztH0VpkKZse+jKUIIIT5wKR42mjRpEiEhIfG2P3nyhMKFC6f08iIdOHXqFP369QNg9OjRDBgwQL/vjrs7br16gVJU+OIL6k6ZEnfiihUQG9QwalRc4KKLgdMD4gKXqn9I4CKEEEIvxT0vSqkEt4eEhOhXmIiM6+7du/pcLm3atOHX2DktAP4eHmzr0AFdVBSO3brRaP58fZ4X1qyB2Ky7DBsG06fHBi7RcKof3Fmr1SaquRIK937/DRNCCPHBeuvgZeRIrX6MkZEREyZMMBg2iomJ4fTp0wnOfxEZx7Nnz3ByciIgIIAKFSro57sABN6+zZaWLYl6/pwCjRrRYtUqjF8seXZ2hr59tTktX34Jv/8eG7hEwYlecO9vMDKB2uuhYJe0a6AQQogP0lsHL+fPnwe0npeLFy9i/iKJGGBubk6FChUYPXp0yp9QfJB0Oh09e/bk0qVL2NnZ4erqSpYsWQAIffSIzc2bExoQQK4KFWi7ZQsmL/59uLpCr15aAroBA7RVRkZGEBMJx7uBz1YwNoM6f0P+dmnXQCGEEB+stw5eDh48CEC/fv2YM2cO1tbWqfZQ4sM3ceJEduzYQaZMmdi2bRv58+cHICo0lK1OTjy9cQPrggXp6OZGphf/Ng4d0lL9x8RoAczixVoel5gIONoJ7u/QCizW3QJ5W6Vd44QQQnzQUjznZcWKFanxHCId2bJlC7/8oi1fXrp0KdWrVwe0Qos7unXD7/RpLGxt6bh7N1kcHLSTzpyB1q0hIgLattUm6xobaz0uxzprgYuJBdTbBvZN06ppQggh0oG3Cl5GjhzJ//73PzJnzqyf+5KYWbNmvdWDiQ/TlStX6NOnD6Atk+/dO24y7aFRo7jl6qpViN6+nRyOjtqOS5egRQsICYFGjbQ5L6amsUNFXcDXVQtc6ruCXeO0aJYQQoh05K2Cl/Pnz+uLLr6Y+5IQI0kmlqEEBgbSrl07QkJCaNCgAdOnT9fvO//HH5ybOxeAlmvXkrdOHW3H7dvQtCk8eQI1aoCLC1hYaJNzj3cDn21gnEnrcZHARQghxGt4q+DlxXyXV/8uMi6dTkevXr24ceMGBQoU4O+//8bMzAwAr927OTB0KAB1p06Nq1cUEKAFLn5+UK4cuLlBliyxgUuP2Mm55lDPRYaKhBBCvLYUJ6kLCwsjNDRU//7u3bvMnj0bd3f3lF5afEB+/vlndu7ciYWFBVu3biVXrlwAPLx4EdcuXVA6HWX79aN6bAkAQkKgZUu4dQsKFYI9e8DWVktAd/Iz8N4UOzl3Kzg0T7uGCSGESHdSHLy0bduW1atXA9qwQvXq1Zk5cyZt27Zl4cKFKX5Akfb27t2rL6y5ZMkSKleuDMBzf3+2ODkRGRxM/gYNaLJokTZUGBkJHTuChwfkzKkFLvb2Wsr/fwbDXWdtOXTdzZC3ZVo2TQghRDqU4uDl3Llz1K1bF4BNmzZhZ2fH3bt3Wb16NXNj50CI9MvX15eePXuilGLw4MH6CbrR4eG4tGtH8L17ZC9enDabN2u5XHQ6LXOuuztYWcHOnVCihJaQ7txIuLVMy5xbez3kdUrj1gkhhEiPUhy8hIaGkjVrVgDc3d3p0KEDxsbG1KxZk7t376b4AUXaiY6Oplu3bjx8+JCKFSsye/ZsQEtMuO+rr7Ql0dmz02HnTixtbbWTxo6FtWu11USbN0PsMmouTIBrc7S/11gOBTq9/wYJIYTIEFIcvBQrVgwXFxe8vb3Zs2cPTZtqEy8DAgIkcV0698MPP3Ds2DGsra3ZuHGjvlbV+fnzubRiBUbGxjj99RfZixfXTvjjD61GEcCyZdA8di7LlWlwWcsLQ9X5UKTPe26JEEKIjCTFwcuECRMYPXo0hQoVokaNGtSqVQvQemEqVaqU4gcUaWPHjh389ttvACxfvpxixYoBcO/gQQ6OGAFAvWnTKNSkiXaCmxvErjhiyhT47DPt7zeXgGfsJN4KU6HE1++tDUIIITKmFGfY7dSpE5988gl+fn5UqFBBv71Ro0a0b98+pZcXaeDevXt8Fht8DBs2jI6xS5+D7tzBtXNnVEwMpXv1ouqLBIWentC1a9x8l7Fjte3eW+CfL7W/lx4HZca+55YIIYTIiFIcvADY2dlhZ2dnsO1FyniRvsTExNCrVy+ePn1K9erVmTZtGgCRz5/j0q4dYY8fk6dKFZosWaKtLPLxASenuOy5CxdqhRYfHILj3bUVRsUGQYXJadswIYQQGUaqBC/79+9n//79BAQEoNPpDPYtX748NW4h3pOpU6dy9OhRsmTJwvr16zE3N0cpxd5Bg3j4779Y5c5N261bMbO0hOBgrV6Rry+UKgWbNoGZGTz1hCNtQRcJ+dpD1QVaQCOEEEKkghQHLz/99BM///wzVatWxd7eXkoCpGOnTp1i0qRJACxYsICiRYsC4LlwIVfXr8fIxITWGzdinT+/Vhm6Rw9tyCh3bm1JdLZsEHwLDjaHqGeQuz7UWQ/GJmnWJiGEEBlPioOXRYsWsXLlSoMCfSL9efbsGT179iQmJobu3bvTq1cvAPzOnOHg8OEA1PvtN/LXq6ed8P33sGOHVqdo+3YoXBjCA+BgMwh/ANkqaPWKTCzSqEVCCCEyqhSvNoqMjKR27dqp8SwiDX3zzTfcvn2bggULsnDhQoyMjAh7/JjtnTuji4qiePv2cRN016+H2JVILF+uFVyMDoXDbSDkFmQuDA13gblN2jVICCFEhpXi4GXAgAGsX78+NZ5FpJENGzawevVqjI2NWbduHTY2NiidDrfevQm+d49sxYrRfMUKbUjwn3+0FUWgrSrq3l2rV3SiJzw+Dea2WuBiaZ+2jRJCCJFhpXjYKDw8nCVLlrBv3z7Kly+vrzT8wqxZs1J6C/EO+fj48OWX2nLmH3/8kTp16gBwasoUvHbtwtTCgjabNpHJxkarDt2uHYSHayuMfolNPHd+FPi4xFWIti6ZJm0RQgjxcUhxz8uFCxeoWLEixsbGXLp0ifPnz+tfnp6eKX7ABQsWULhwYSwsLKhSpQpHjx5N8viIiAi+//57ChYsSKZMmShatKiseEqEUooBAwYQFBRE9erV+eGHHwAtEd3xCRMAaLxwIbkrVNAClvbt4f59bWXRunVgYgL/zYlL+19rNeSum1bNEUII8ZFIcc/LwYMHU+M5EvTXX38xfPhwFixYQJ06dVi8eDEtWrTgypUrFChQIMFzunTpwoMHD1i2bBnFihUjICCA6Ojod/aM6dmff/7Jnj17sLCwYNWqVZiamhL66BFuvXqBUpTt35+yfftqRRW//hpOn4bs2bUJutbW4O0C57Rsu1T8FQp2TcvmCCGE+EgYKaVUWj9EYmrUqEHlypVZuHChflupUqVo164dU6dOjXf87t276datG7dv38b2RaHAN+Tj40P+/Pnx9vYmX758b/3sH7q7d+9StmxZQkJCmDlzJiNHjkQphUvbttxydcXW0ZFeZ89injkzLF0KgwaBsTHs2QONG8OTc7D3E4gJg2KDodpCyeUihBAfuff1MzTFw0YAR48epVevXtSqVQtfX18A1qxZw7Fjx976mpGRkXh4eOgLPb7QtGlTTpw4keA527dvp2rVqkybNo28efNSokQJRo8eTVhYWKL3iYiI4NmzZ/pXcHDwWz9zeqHT6ejfvz8hISHUqVOHYcOGAXD+jz+45eqKibk5Ts7OWuDyzz8wZIh24uTJWuAS5qcloYsJA7umWrFFCVyEEEK8JykOXjZv3kyzZs2wtLTk/PnzREREABAcHMyUKVPe+rqPHj0iJiaGPHnyGGzPkycP/v7+CZ5z+/Ztjh07xqVLl9i6dSuzZ89m06ZNfP114sUAp06dio2Njf5VunTpt37m9GLRokUcOHAAS0tLVq5ciYmJCQH//svh0aMBqD9jhjbP5dEj6NQJIiOhbVv47juICYcj7SHUR5uY+8lfYJwqiZqFEEKI15Li4OWXX35h0aJFLF261GClUe3atTl37lxKLx8vY69SKtEsvjqdDiMjI9atW0f16tVp2bIls2bNYuXKlYn2vowbN46goCD968qVKyl+5g/ZrVu3+PbbbwH47bffKFasGJHPn7OjWzdiIiIo4uREpSFDtAy6PXvCvXtQrBisWqVd4NTnsUuis0M9VzDPlnaNEUII8VFKcfBy7do16r3IuvoSa2trAgMD3/q6OXPmxMTEJF4vS0BAQLzemBfs7e3JmzcvNjZxydFKlSqFUgofH58Ez8mUKRPW1tb6V9asWd/6mT90L1YXhYaG0qBBA32P1MHhw3ny339ktrePy+cyaRK4u4OlJWzZAjY2cGUq3F0PRqbwySawLp62DRJCCPFRSnHwYm9vz82bN+NtP3bsGEWKFHnr65qbm1OlShX27t1rsH3v3r2JZvStU6cO9+/fJyQkRL/t+vXrGBsbZ+jJt69r5cqVHDp0CEtLS5YtW4axsTG3XF25+OefYGREq7VrscqZE3btisvhsnQplCsHPtvg3++1bVXngd2nadcQIYQQH7UUBy+DBw9m2LBhnD59GiMjI+7fv8+6desYPXo0X331VYquPXLkSP7888//t3fncVWW+f/HX4cdEVBxAQMJF9xXNMUlUxOXFq0mLR2Xyin7uoRaaVq51Py0abJt0jaXaoqsURlLTMndNDNFszDLXMDCTCVESQTO9fvjABOyiAicc+T9fDx4yLnvz3Wfz8U1Z86n614uFi9ezIEDB5g0aRJJSUmMHTsWsJ3yGTlyZH78sGHDCAgI4L777iMxMZEtW7bw2GOPcf/99+Pt7X1VuTi7kydPMmXKFMC2mGbDhg3J+O031o4ZA0DHKVNo0Lu3bYXovL/p//2f7dRR2vewPXftqvDx0GSsPbogIiIClMNzXh5//HHS0tLo1asXFy5c4MYbb8TT05NHH32U8Xl3qZTR0KFDOX36NHPmzCElJYVWrVoRFxdHaGgoACkpKSQlJeXHV69enfj4eCZMmEDHjh0JCAhgyJAhPJs3i1CFTZ48mdTUVNq1a8ekSZMwxhD/8MNknDxJQMuWdH/mGdt1Ln/9q+1C3Xbt4IUXbKtDb70DstNtq0R30BOTRUTEvsrtOS8ZGRkkJiZitVpp0aIF1atXL4/DVrpr8Tkva9eupX///ri4uPDll1/SqVMnEt9/n7i//hUXNzeGf/UV9dq3h9mzbde6+PjAnj3QpDFs/QscXwne10H/3eBd9PVGIiIilfUdelUzL1arlaVLl7JixQqOHj2KxWIhLCyMv/zlL4wYMaLYu4Kk8mRkZOSvXTRhwgQ6depE+vHjrM+9WDdy5kxb4bJ5M8yZY2v0+usQHg7f/T9b4eLiAT1WqHARERGHUOZrXowx3H777YwZM4aff/6Z1q1b07JlS44dO8bo0aO54447yjNPKaPZs2dz5MgRQkJCeOaZZzDGsPaBB8hMSyPwhhvoPG0a/PYbDBsGVivcd5/t1NEvn8E+21pHdHwNat9g346IiIjkKvPMy9KlS9myZQvr16+nV69eBfZt2LCBwYMH8+677xa4oFYq1759+3jhhRcAeO211/D19WXvwoUcXbcONy8vBr77Li6urjB69P8WXHz1VTh3FLYPAww0+hs0HmPPboiIiBRQ5pmXmJgYpk+fXqhwAejduzfTpk3j/fffv6rkpOyMMUyYMIGcnBzuuusubrvtNs4mJbH58ccBuPG556jVtCksWABxceDpCcuWgZc7bBsCF1Mh4AbbbdEiIiIOpMzFyzfffEP//v2L3T9gwAD27dtX1sPLVYqJiWHr1q1Uq1aNF198Mf/uoqxz57iuWzfbU3S//x5ylwTgH/+wPc8l4TE4s8v2BN3uH4Orp307IiIicokyFy9nzpwp9km3YFuDKDU1tayHl6tw7ty5/CUApk+fTkhICN/HxHAkLg5XDw+i3n4bS3a27dqWCxegb1/b4otJy+GHV2wHiXwXfBrYsRciIiJFK/M1Lzk5Obi5Fd/c1dWV7Ozssh5ersKzzz7LL7/8QsOGDZkyZQoZv/3GhokTAejy1FMENGsGTz4Ju3dDzZqwZAlkHIWd99sO0PwxuO5W+3VARESkBGUuXowxjB49Gk/Pok8r5K0uLZXrhx9+YP5824PkXnrpJby8vFj9wAP8cfo0tVu35obHH4ft22HuXFuDN96AwNoQ3832QLraXaHt3+3YAxERkZKVuXgZNWrUZWN0p1HlMsYQHR1NVlYWAwYM4NZbb+Wn1as58MEHWFxc6LdoEa6ZmTBihO226BEj4O67Ydd4OLMbPAOg+zJwcb/8m4mIiNhJmYuXJUuWlGceUg4+/fRT1qxZg7u7Oy+99BJZ587xee4D6jpERxPUqRM8+CAcPgyhobbboo+vgh9fsx0g8j2odm08VVhERK5dV722kTiGzMxMoqOjAds6RuHh4WycPJn05GT8w8LoNmcOxMfbVokGeOcd8MiAnQ/YXjebAvUH2Cd5ERGRK6Di5RqxYMECDh8+TFBQEE8++SS/7d/Pnldsdw7dvGABHlYr/O1vtuDx4+HGHrBxAGSegprtdJ2LiIg4DRUv14DU1FSeeeYZAObMmYOPjw+f/N//YXJyaHLnnYT17w/jxsGxYxAWZrtY9+ArcGIduHpD1w/0PBcREXEaZX7OiziOefPmkZqaSsuWLRk9ejSJ773Hz9u24VatGr1eesm26OKCBbbgt96CrJ9g71Tb6w7zwb+53XIXERG5Upp5cXLHjh3j5ZdfBuC5554j+9w5Nuc+oC7y6afxCwiAvCUcHnwQenaFtR3BehGuux0aP2Sv1EVERMpExYuTe+qpp8jMzKRXr14MHDiQDRMnknHyJLWaNaPjpEkwdSr89BMEB8Pzz9tmXNISwSsQOr8NFou9uyAiInJFLMYYY+8kHMnx48cJCQkhOTmZ4GDHvm04ISGBiIgIjDHs2rWLEFdX/t2xI8ZqZcj69TSoVg26dgVjYM0aaO8J63vbGt/0GdTvZ98OiIjINaWyvkM18+KkjDE8/vjjGGO49957iejQgZju3TFWK83uuYcGPXpARIStcBk1Cvp0g7g2tsaNH1LhIiIiTksX7DqpdevW8fnnn+Ph4cHf//53vv/wQ37ZsQP36tXp+c9/wssvw/79EBAA//ynbbXo80fB53po/7y90xcRESkzFS9OyGq1Mm3aNADGjx9PcGAgW3Jfd5k+Hd+cHJg50xb8/POQtQcOvWF73WUxuPvaI20REZFyodNGTig2Npa9e/fi6+vL9OnT2T1/PunJyfg2aECH6Gi4917IyIAePeDewfBZ7umi8PFQr5c9UxcREblqmnlxMlarlVmzZgEQHR2NV1YWO3NXiL5x3jzc162D//4X3Nxg4ULYOwUyjkP1RtBunh0zFxERKR+aeXEyy5cvZ//+/fj7+zNp0iS2Pf44WefPE9S5M81uvRVatrQFPvYY1EyCfUsAC3RZCm4+9kxdRESkXKh4cSI5OTn5sy6TJk0iKymJ/YsWAXDT/PlYZs+G5GS4/nqYFg0bb7A1bPoI1O1ul5xFRETKm4oXJ/LRRx+RmJhIjRo1eOSRR1j3l7+AMTQdMoTr/P3hpZdsga+9Boeeg/PHwCcU2jxj17xFRETKk655cRI5OTnMnj0bgClTpnB62zaS1q/H1cODHnPnwqRJkJMDgwZBl3pw8CVbw04Lwb26/RIXEREpZ5p5cRIxMTEcPHiQWrVqMX7cOFZ2t50Gipg0iRrffgvx8eDhAf98DnbeC8YKofdC/QF2zlxERKR8qXhxAtnZ2fmzLo899hg/f/oppxMT8apZkxsmTYJu3WyBU6ZA1ieQmgAeNaHDi3bMWkREpGKoeHEC77//PocOHaJ27do8/OCDfNSxIwCdpk7F6513bAsvBgXBI/fCps62Ru3/Cd717Ji1iIhIxVDx4uCsVivPPfccAI8++ihHli0j7cgRfAIDaX/33dC2rS1w3lxIfAxy/oC6N0HD++yXtIiISAVS8eLgPv30Uw4cOICfnx9jRo3iow4dAOjy5JN4PPssnDsHN9wAPavBF2vBxRNueAMsFjtnLiIiUjF0t5EDM8Ywb57tqbgPP/wwh//9b86npOB3/fW0ad8eliyxBb74HCRMsf3eYir4hdspYxERkYqnmRcHtm3bNnbs2IGHhwf/d//9rIqMBKDr00/jOnWqLWjECPBZDxnJtme6tJhqx4xFREQqnooXB5Z3rcvo0aM5HhPDhTNnqNWsGS1q1oRt28DbG55+EL7uY2vQ4UVwq2bHjEVERCqeihcH9e2337J69WosFgsT7r+fDX37AtBt1ixcZsywBUVHw/F5YL0IgVEQPNhu+YqIiFQWXfPioP7xj38AcNddd3EmNpaL6enU69CB8HPnIDERataEUW3gl9Xg4g4dX9FFuiIiUiVo5sUBHTt2jJiYGAAmP/wwuwYPBqDrtGlYJk+2Bc14HA7mzsA0nQR+Te2QqYiISOXTzIsDmj9/PtnZ2fTp0wezbRsX09Op06YNDY8cgePHISQE+lyAc4fBuz60etLeKYuIiFQazbw4mDNnzvD2228D8OiECey5/34AIidOxPLYY7agZ6bAwem239s/D+6+9khVRETELjTz4mDefvttMjIyaNu2LbUOHuTCmTPUDA+nycGDkJoKrVpB028hJwMCOtsWXxQREalCHL54WbBgAWFhYXh5eREREcHWrVtL1e6LL77Azc2Ndu3aVWyC5SgnJ4cFCxYAMGHsWL6ePx+Abg89hOXVV21Bz46FI4ttv3d4QRfpiohIlePQxcuyZcuIjo5mxowZJCQk0KNHDwYMGEBSUlKJ7dLS0hg5ciR9+vSppEzLx6effsqxY8cICAig9R9/kPHrr/iFhhJ+8CBcuADdu4PfajBWCLkL6nSzd8oiIiKVzqGLl/nz5/PAAw8wZswYmjdvzksvvURISAgLFy4ssd1DDz3EsGHDiMx9Iq2zeDV3dmXMffeR8OKLAHQdMwaXvGUAnhwMKWvA4gbt5tkpSxEREfty2OLl4sWL7N69m6ioqALbo6Ki2L59e7HtlixZwk8//cTMmTNL9T6ZmZmcPXs2/yc9Pf2q8i6rxMRE1q9fj4uLC/3r1iU9ORmfwEBaJCdDVhb0vgly3rMFh48D38Z2yVNERMTeHLZ4OXXqFDk5OdSrV6/A9nr16nHixIki2/z4449MmzaN999/Hze30t1INXfuXPz9/fN/WrRocdW5l8W//vUvAAbdfjs/vfUWcMmsy5Qu8Ps+cPeHVk/ZJUcRERFH4LDFSx7LJRekGmMKbQPbxa7Dhg1j9uzZhIeXflXlJ554grS0tPyfxMTEq875SqWlpfHuu+8C8Nd27Uj98Ue8atWi1a+/2mZd+vSAzNxZl5YzwDOg0nMUERFxFA77nJfatWvj6upaaJbl5MmThWZjANLT0/n6669JSEhg/PjxAFitVowxuLm5sW7dOnr37l2onaenJ56envmvz549W849ubylS5dy/vx5WrZsyYX4eABuGD4c1zfesAU80gzSt9pWjW46odLzExERcSQOO/Pi4eFBREQE8blf5nni4+Pp2rVroXg/Pz/279/P3r1783/Gjh1L06ZN2bt3L507d66s1K+I1WrNP2X08K238ssXX+Di7k679HS4eBFu7gaZ/7EFt3kGXL3smK2IiIj9OezMC8DkyZMZMWIEHTt2JDIykjfffJOkpCTGjh0L2E75/Pzzz7z77ru4uLjQqlWrAu3r1q2Ll5dXoe2OZO3atRw6dAh/f3+uO3yYQ0C7227D44MPbAETGsO5L8CvOYQOs2uuIiIijsChi5ehQ4dy+vRp5syZQ0pKCq1atSIuLo7Q0FAAUlJSLvvMF0eXN+sy5u67Obx0KQBd3N1tsy59I+HCCltgmzng4mqnLEVERByHxRhj7J2EIzl+/DghISEkJycTHBxcoe+VnJxMaGgoxhhiH3yQH998kyaRkQzavdtWvMTeC+djoEZbGLAHLA57lk9ERKTSvkP1bWhHS5cuxRhDn+7dSfroIwBuDArKvdalE1z4ry2wzTMqXERERHLpG9FOrFYrixfb1ii6u2FDMn//nToNG1Jj/XpbwIOBuYsv3gDX3WrHTEVERByLihc72bRpE0ePHsXP1xeX3CcG92nTBktaGnRqBDnrbIFtntHiiyIiIn+i4sVOFi1aBMCYbt1IO3QIb39/rtu1y7ZzfAhYM6FODwjsa8csRUREHI+KFztITU1l+fLlADQ/dQqAPt26Yfn5Z2heB9y/sAW2fVazLiIiIpdQ8WIHMTExZGZm0i08nN+//hqLxUKTn36y7XykMZgsqNcb6t5o30RFREQckIoXO8g7ZXRn/foARHbqhOvBgxBUDfwTbEEtZ9grPREREYem4qWS7d27lz179lDN3R33ffsAaH/hgm3nlFZgvQABnaFeLztmKSIi4rhUvFSyvFmXke3bk5maSli9enh/8w34uULwAVtQyyd0rYuIiEgxVLxUogsXLvD+++8D0Pr8eQB61qlj2zmpNeSkg39LuO42e6UoIiLi8FS8VKLY2FhSU1NpGxhIxnff4eviQsCBA+ABND9mC2rxhJ6mKyIiUgJ9S1ai9957D4C7Q0IAuKlJEyw5OTCmEeSkgk8YhA61Z4oiIiIOz6FXlb6WnD59mnXr1uEOVD94ECvQ+NdfwRXo9rstqMVUcNGQiIiIlEQzL5XkP//5D9nZ2dwWEkLW2bO0rVMH199/h1tqAKfBOwgajrJzliIiIo5PxUsliYmJAaBL7l1Enby9wQLckTvT0mwKuHrZKTsRERHnoeKlEvz8889s2bKFIMCSlESgqyu+SUnQ3hU8ToG7HzR+0N5pioiIOAUVL5Xgo48+whjDoKAgAHrkXrDLiNzbpBuNAXdfO2UnIiLiXFS8VIKYmBhcgUbp6XgBIb/8AtcBdU/YbosOn2DnDEVERJyHipcKdujQIXbt2kUziwXOnSPC1xeXixdheIAtIPgOqH69XXMUERFxJipeKtiHH34IQP86dbAA7VxcoDrQ5qwtoGm0vVITERFxSipeKpAxhpiYGLyBwDNnCAO809LgFi+wZEGtCKjTzd5pioiIOBUVLxVo//79JCYm0sHVFbKz6eTra3soXf/c26ObRmsBRhERkSuk4qUC5T3bpbe/Pz5A8Pnz0BnwOAdegdBgiF3zExERcUYqXiqIMYYPP/yQAMDvzBlaWixYrFb4i48tIHwcuHrYNUcRERFnpOKlgnz11VccPXqULh62AqW9tzeEA/XO256k2/gh+yYoIiLipFS8VJCVK1cC0MXdnfqAb0YGDHC17bx+OHjVsV9yIiIiTkzFSwWJjY3lesDr/HnaurmBL9DR2HY2GWfHzERERJybipcK8P3333Pw4EE6urjgDjQFuBFwsUKtTlCrvX0TFBERcWIqXipAbGwsbkAHFxeaAm452dDP3bazia51ERERuRoqXirAf//7X5oC7tnZtPXwgBZAQJZt9ejQe+ydnoiIiFNT8VLOUlJS+PLLL2kD1AKCLl6EPrkPorv+r+DmY8/0REREnJ6Kl3K2atUq3IDWLi60AvADOuXubPyg3fISERG5Vqh4KWexsbE0BTytVlq5uEBPwMVAQGeo2dbe6YmIiDg9FS/l6OzZs6xfv57WwPVANWOFm3P/xHoonYiISLlQ8VKO1qxZgzUri9YuLjQDaAnUtoK7P4QOtXN2IiIi1wYVL+UoNjaWcMDHaqWJxQJ9cneEjQC3avZMTURE5Jqh4qWcZGZmsnr1atoCDQF3PwMRuTt1ykhERKTcqHgpJ5s2bSIjPZ2WFgvNAboBrkDtSKjRyr7JiYiIXENUvJST2NhYmgA1jCEMbMsBADQcbbecRERErkUqXsqBMYZVq1bRFmgCuIYCIYCLJzQYYt/kRERErjEOX7wsWLCAsLAwvLy8iIiIYOvWrcXGrlixgr59+1KnTh38/PyIjIxk7dq1FZ7jN998w6+//EJLsJ0y6pG7I3gQeNSo8PcXERGpShy6eFm2bBnR0dHMmDGDhIQEevTowYABA0hKSioyfsuWLfTt25e4uDh2795Nr169uO2220hISKjQPNeuXUtjoDYQ4gp0zd0RNrJC31dERKQqshhjjL2TKE7nzp3p0KEDCxcuzN/WvHlzBg8ezNy5c0t1jJYtWzJ06FCefvrpUsUfP36ckJAQkpOTCQ4OLlWb3r17U3vjRh4GerUHHgW86sLg4+DiXqpjiIiIOLuyfIeWhcPOvFy8eJHdu3cTFRVVYHtUVBTbt28v1TGsVivp6enUqlWr2JjMzEzOnj2b/5Oenn5Feaanp7N961Zagu3BdN1zd4QOV+EiIiJSARy2eDl16hQ5OTnUq1evwPZ69epx4sSJUh3jhRde4Pz58wwZUvxFs3PnzsXf3z//p0WLFleU58aNG6mfnc11QFA1MB1ydzTUKSMREZGK4LDFSx6LxVLgtTGm0LaixMTEMGvWLJYtW0bdunWLjXviiSdIS0vL/0lMTLyi/NauXUsLcmdduoDFA6jRGmpoEUYREZGK4GbvBIpTu3ZtXF1dC82ynDx5stBszKWWLVvGAw88wMcff8zNN99cYqynpyeenp75r8+ePVvqHI0xrFmzhr9wyV1GYSOhFAWWiIiIXDmHnXnx8PAgIiKC+Pj4Atvj4+Pp2rVrMa1sMy6jR4/mgw8+4JZbbqnQHA8dOkT6kSM0BwLqAeEALnD98Ap9XxERkarMYWdeACZPnsyIESPo2LEjkZGRvPnmmyQlJTF27FjAdsrn559/5t133wVshcvIkSN5+eWX6dKlS/6sjbe3N/7+/uWe32effUZzoDH870LdoCjwDir39xIREREbhy5ehg4dyunTp5kzZw4pKSm0atWKuLg4QkNDAUhJSSnwzJc33niD7Oxsxo0bx7hx4/K3jxo1iqVLl5Z7fnnXuzSB/xUveraLiIhIhXLo57zYQ2nvUb9w4QLX1arFs3/8wcONgDmAazW46zdwq1Zp+YqIiDiKKv+cF0e3bds2Qv74g6YAXXI3Bg9S4SIiIlLBVLyUUf4pIwtwQ+5GLcIoIiJS4VS8lNHaNWuIAILzFjVyrQ71+9s5KxERkWufipcyOH78OOe++45WgKVz7saQweDqZcesREREqgYVL2VQ4JRRXvGiU0YiIiKVQsVLGaxbt472QGgToBbg6mt7vouIiIhUOBUvV8gYw97PP6cz4Jp3l1GDO8HVs6RmIiIiUk5UvFyhxMRE6p45o7uMRERE7ETFyxXatGkTbYCGTYGagIsfBJa8+KOIiIiUHxUvV2jT559zE+Ced8ro+rvA1cOOGYmIiFQtKl6ugNVq5eiGDTR34U+njIbaMyUREZEqR8XLFfjuu+8IPHuWRs0Af8DiB4G97Z2WiIhIlaLi5Qps3LiRXoB3/imjv4CLuz1TEhERqXJUvFyBLz77jM4WoEPuhlDdZSQiIlLZVLyUktVqJWXLFsIbYrvLyHhBvZvsnJWIiEjVo+KllL755huanj9P7YjcDXX76sF0IiIidqDipZQ2btxIP8CSV7w0ucee6YiIiFRZKl5K6atVq2hXDwgGY1yg/kB7pyQiIlIlqXgphZycHFJ37CA0d9bF4t0ePGrYNScREZGqSsVLKezbt4+emZm4550yav5Xu+YjIiJSlal4KYUN8fH0rg6E525ocIc90xEREanSVLyUwt6VK2neAXABYxqAT6i9UxIREamyVLxcRnZ2Nh579uCX+2A6S6O77ZuQiIhIFafi5TISEhLoTxa0zt0QPsyu+YiIiFR1Kl4uY8uaNXRrBXiBuegLNdvbOyUREZEqTcXLZRz65BMC826RrtsfLBb7JiQiIlLFqXgpgTGGBt9+g2veQoztHrBrPiIiIqLipURHjhzh5uCL4A8m0xXq9bJ3SiIiIlWeipcSbPvsM5q3s/1utbQFVw+75iMiIiIqXkr043/+g08b2++ubUbYNxkREREBVLyUqOH3O7E0yn3RfKhdcxEREREbFS/FOHv2LN3CMgDIOl0dvIPsnJGIiIiAipdifb11K6Ftbb9bavS0bzIiIiKST8VLMY58/BGeude7uHV72L7JiIiISD4VL8UI++Ur8APrH0BoX3unIyIiIrlUvBSjRcg5AP74tbZukRYREXEgKl6KUaeF7V/XhoPtmoeIiIgUpOKlGK5htn+9oh6zbyIiIiJSgIqX4ljgYooL1Aq3dyYiIiLyJypeSpD2m57tIiIi4mhUvJTAvcN99k5BRERELuHwxcuCBQsICwvDy8uLiIgItm7dWmL85s2biYiIwMvLi4YNG/L666+X6X3NRahxxxNlaisiIiIVx6GLl2XLlhEdHc2MGTNISEigR48eDBgwgKSkpCLjjxw5wsCBA+nRowcJCQlMnz6diRMnsnz58it+7/NHLOBe7Wq7ICIiIuXMYowx9k6iOJ07d6ZDhw4sXLgwf1vz5s0ZPHgwc+fOLRQ/depUVq1axYEDB/K3jR07ln379rFjx45Svefx48cJCQlh64gGdH/32NV3QkREpIrI+w5NTk4mODi4wt7HrcKOfJUuXrzI7t27mTZtWoHtUVFRbN++vcg2O3bsICoqqsC2fv36sWjRIrKysnB3dy/UJjMzk8zMzPzXaWlpAJxufAfHjx+/2m6IiIhUGSkpKQBYrdYKfR+HLV5OnTpFTk4O9erVK7C9Xr16nDhxosg2J06cKDI+OzubU6dOERRU+O6huXPnMnv27ELbB898GWa+fBU9EBERqZqSk5Np0KBBhR3fYYuXPBaLpcBrY0yhbZeLL2p7nieeeILJkyfnvz5z5gxhYWF8++23+Pv7lzVth5eenk6LFi1ITEzE19fX3ulUGPXz2lNV+qp+XluqSj/T0tJo1aoVzZs3r9D3cdjipXbt2ri6uhaaZTl58mSh2ZU8gYGBRca7ubkREBBQZBtPT088PT0LbQ8JCcHPz6+M2Tu+s2fPAnDdddepn9eAqtJPqDp9VT+vLVWln3l9c3Or2PLCYe828vDwICIigvj4+ALb4+Pj6dq1a5FtIiMjC8WvW7eOjh07Fnm9i4iIiDgfhy1eACZPnszbb7/N4sWLOXDgAJMmTSIpKYmxY8cCtlM+I0eOzI8fO3Ysx44dY/LkyRw4cIDFixezaNEiHn30UXt1QURERMqZw542Ahg6dCinT59mzpw5pKSk0KpVK+Li4ggNDQVsVzX/+ZkvYWFhxMXFMWnSJF577TXq16/PK6+8wl133VXq9/T09GTmzJlFnkq6lqif15aq0k+oOn1VP68t6mf5cujnvIiIiIhcyqFPG4mIiIhcSsWLiIiIOBUVLyIiIuJUVLyIiIiIU6mSxcuCBQsICwvDy8uLiIgItm7dWmL85s2biYiIwMvLi4YNG/L6669XUqZlM3fuXDp16oSvry9169Zl8ODBHDx4sMQ2mzZtwmKxFPr5/vvvKynrKzdr1qxC+QYGBpbYxtnGEuD6668vcmzGjRtXZLwzjeWWLVu47bbbqF+/PhaLhdjY2AL7jTHMmjWL+vXr4+3tzU033cR333132eMuX76cFi1a4OnpSYsWLVi5cmUF9aB0SupnVlYWU6dOpXXr1vj4+FC/fn1GjhzJL7/8UuIxly5dWuQ4X7hwoYJ7U7zLjefo0aML5dulS5fLHteZxhMoclwsFgvPP/98scd0xPEszXeJvT6jVa54WbZsGdHR0cyYMYOEhAR69OjBgAEDCtxy/WdHjhxh4MCB9OjRg4SEBKZPn87EiRNZvnx5JWdeeps3b2bcuHF8+eWXxMfHk52dTVRUFOfPn79s24MHD5KSkpL/06RJk0rIuOxatmxZIN/9+/cXG+uMYwmwa9euAn3MexDj3XffXWI7ZxjL8+fP07ZtW/71r38Vuf8f//gH8+fP51//+he7du0iMDCQvn37kp6eXuwxd+zYwdChQxkxYgT79u1jxIgRDBkyhJ07d1ZUNy6rpH5mZGSwZ88ennrqKfbs2cOKFSv44YcfuP322y97XD8/vwJjnJKSgpeXV0V0oVQuN54A/fv3L5BvXFxcicd0tvEECo3J4sWLsVgsl31sh6ONZ2m+S+z2GTVVzA033GDGjh1bYFuzZs3MtGnTiox//PHHTbNmzQpse+ihh0yXLl0qLMfydvLkSQOYzZs3FxuzceNGA5jU1NTKS+wqzZw507Rt27bU8dfCWBpjzCOPPGIaNWpkrFZrkfudcSyNMQYwK1euzH9ttVpNYGCgmTdvXv62CxcuGH9/f/P6668Xe5whQ4aY/v37F9jWr18/c88995R7zmVxaT+L8tVXXxnAHDt2rNiYJUuWGH9///JNrhwV1c9Ro0aZQYMGXdFxroXxHDRokOndu3eJMY4+nsYU/i6x52e0Ss28XLx4kd27dxMVFVVge1RUFNu3by+yzY4dOwrF9+vXj6+//pqsrKwKy7U8paWlAVCrVq3LxrZv356goCD69OnDxo0bKzq1q/bjjz9Sv359wsLCuOeeezh8+HCxsdfCWF68eJF///vf3H///SUuUArON5aXOnLkCCdOnCgwZp6envTs2bPYzysUP84ltXE0aWlpWCwWatSoUWLcuXPnCA0NJTg4mFtvvZWEhITKSfAqbNq0ibp16xIeHs7f/vY3Tp48WWK8s4/nr7/+yurVq3nggQcuG+vo43npd4k9P6NVqng5deoUOTk5hRZ2rFevXqEFHfOcOHGiyPjs7GxOnTpVYbmWF2MMkydPpnv37rRq1arYuKCgIN58802WL1/OihUraNq0KX369GHLli2VmO2V6dy5M++++y5r167lrbfe4sSJE3Tt2pXTp08XGe/sYwkQGxvL77//zujRo4uNccaxLEreZ/JKPq957a60jSO5cOEC06ZNY9iwYSUu4NesWTOWLl3KqlWriImJwcvLi27duvHjjz9WYrZXZsCAAbz//vts2LCBF154gV27dtG7d28yMzOLbePs4/nOO+/g6+vLnXfeWWKco49nUd8l9vyMOvTyABXl0v9iNcaU+F+xRcUXtd0RjR8/nm+++YZt27aVGNe0aVOaNm2a/zoyMpLk5GT++c9/cuONN1Z0mmUyYMCA/N9bt25NZGQkjRo14p133mHy5MlFtnHmsQRYtGgRAwYMoH79+sXGOONYluRKP69lbeMIsrKyuOeee7BarSxYsKDE2C5duhS42LVbt2506NCBV199lVdeeaWiUy2ToUOH5v/eqlUrOnbsSGhoKKtXry7xy91ZxxNg8eLFDB8+/LLXrjj6eJb0XWKPz2iVmnmpXbs2rq6uhaq7kydPFqoC8wQGBhYZ7+bmRkBAQIXlWh4mTJjAqlWr2LhxI8HBwVfcvkuXLg5T9ZeGj48PrVu3LjZnZx5LgGPHjvH5558zZsyYK27rbGMJ5N85diWf17x2V9rGEWRlZTFkyBCOHDlCfHx8ibMuRXFxcaFTp05ONc5BQUGEhoaWmLOzjifA1q1bOXjwYJk+s440nsV9l9jzM1qlihcPDw8iIiLy79bIEx8fT9euXYtsExkZWSh+3bp1dOzYEXd39wrL9WoYYxg/fjwrVqxgw4YNhIWFlek4CQkJBAUFlXN2FSczM5MDBw4Um7MzjuWfLVmyhLp163LLLbdccVtnG0uwLbQaGBhYYMwuXrzI5s2bi/28QvHjXFIbe8srXH788Uc+//zzMhXTxhj27t3rVON8+vRpkpOTS8zZGcczz6JFi4iIiKBt27ZX3NYRxvNy3yV2/YyW+tLea8SHH35o3N3dzaJFi0xiYqKJjo42Pj4+5ujRo8YYY6ZNm2ZGjBiRH3/48GFTrVo1M2nSJJOYmGgWLVpk3N3dzX/+8x97deGyHn74YePv7282bdpkUlJS8n8yMjLyYy7t54svvmhWrlxpfvjhB/Ptt9+aadOmGcAsX77cHl0olSlTpphNmzaZw4cPmy+//NLceuutxtfX95oayzw5OTmmQYMGZurUqYX2OfNYpqenm4SEBJOQkGAAM3/+fJOQkJB/l828efOMv7+/WbFihdm/f7+59957TVBQkDl79mz+MUaMGFHgbsEvvvjCuLq6mnnz5pkDBw6YefPmGTc3N/Pll19Wev/ylNTPrKwsc/vtt5vg4GCzd+/eAp/ZzMzM/GNc2s9Zs2aZzz77zPz0008mISHB3HfffcbNzc3s3LnTHl00xpTcz/T0dDNlyhSzfft2c+TIEbNx40YTGRlprrvuumtqPPOkpaWZatWqmYULFxZ5DGcYz9J8l9jrM1rlihdjjHnttddMaGio8fDwMB06dChwC/GoUaNMz549C8Rv2rTJtG/f3nh4eJjrr7++2P8xOgqgyJ8lS5bkx1zaz+eee840atTIeHl5mZo1a5ru3bub1atXV37yV2Do0KEmKCjIuLu7m/r165s777zTfPfdd/n7r4WxzLN27VoDmIMHDxba58xjmXdb96U/o0aNMsbYbsWcOXOmCQwMNJ6enubGG280+/fvL3CMnj175sfn+fjjj03Tpk2Nu7u7adasmd0Lt5L6eeTIkWI/sxs3bsw/xqX9jI6ONg0aNDAeHh6mTp06Jioqymzfvr3yO/cnJfUzIyPDREVFmTp16hh3d3fToEEDM2rUKJOUlFTgGM4+nnneeOMN4+3tbX7//fcij+EM41ma7xJ7fUYtuQmKiIiIOIUqdc2LiIiIOD8VLyIiIuJUVLyIiIiIU1HxIiIiIk5FxYuIiIg4FRUvIiIi4lRUvIiIiIhTUfEiIiIiTkXFi4hUCYsWLSIqKir/9ejRoxk8eHCx8bNmzaJdu3YVksv+/fsJDg7m/PnzFXJ8kWudiheRKmD06NFYLBYsFgvu7u40bNiQRx99tMp8eWZmZvL000/z1FNP2TsVAFq3bs0NN9zAiy++aO9URJySiheRKqJ///6kpKRw+PBhnn32WRYsWMCjjz5aZGxWVlYlZ3d5V5PT8uXLqV69Oj169CjHjMomrx/33XcfCxcuJCcnx84ZiTgfFS8iVYSnpyeBgYGEhIQwbNgwhg8fTmxsLPC/UySLFy+mYcOGeHp6YowhLS2NBx98kLp16+Ln50fv3r3Zt29f/jH37dtHr1698PX1xc/Pj4iICL7++msAjh07xm233UbNmjXx8fGhZcuWxMXFAbB06VJq1KhRIL/Y2FgsFkv+67LmVJQPP/yQ22+/vcSY3bt3U7duXf7+978X2P7GG28QEhJCtWrVuPvuu/n999/z9+3atYu+fftSu3Zt/P396dmzJ3v27CnQ3mKx8PrrrzNo0CB8fHx49tlnAejXrx+nT59m8+bNJeYlIoWpeBGpory9vQvMZhw6dIiPPvqI5cuXs3fvXgBuueUWTpw4QVxcHLt376ZDhw706dOHM2fOADB8+HCCg4PZtWsXu3fvZtq0abi7uwMwbtw4MjMz2bJlC/v37+e5556jevXqV5RjWXIqytatW+nYsWOx+zdt2kSfPn2YPXs2M2bMKPT+n3zyCZ999hl79+5l3Lhx+fvT09MZNWoUW7du5csvv6RJkyYMHDiQ9PT0AsefOXMmgwYNYv/+/dx///0AeHh40LZtW7Zu3XpFfxMRAcqyTLaIOJdRo0aZQYMG5b/euXOnCQgIMEOGDDHGGDNz5kzj7u5uTp48mR+zfv164+fnZy5cuFDgWI0aNTJvvPGGMcYYX19fs3Tp0iLfs3Xr1mbWrFlF7luyZInx9/cvsG3lypXmz/+XVNacLpWammoAs2XLlgLb8/4msbGxxtfX13zwwQcF9s+cOdO4urqa5OTk/G1r1qwxLi4uJiUlpcj3ys7ONr6+vuaTTz7J3waY6OjoIuPvuOMOM3r06CL3iUjx3OxcO4lIJfn000+pXr062dnZZGVlMWjQIF599dX8/aGhodSpUyf/9e7duzl37hwBAQEFjvPHH3/w008/ATB58mTGjBnDe++9x80338zdd99No0aNAJg4cSIPP/ww69at4+abb+auu+6iTZs2V5RzWXK61B9//AGAl5dXoX07d+7k008/5eOPP+aOO+4otL9BgwYEBwfnv46MjMRqtXLw4EECAwM5efIkTz/9NBs2bODXX38lJyeHjIwMkpKSChynuFkfb29vMjIyium9iBRHxYtIFdGrVy8WLlyIu7s79evXzz+9k8fHx6fAa6vVSlBQEJs2bSp0rLzrVWbNmsWwYcNYvXo1a9asYebMmXz44YfccccdjBkzhn79+rF69WrWrVvH3LlzeeGFF5gwYQIuLi4YYwocs6gLcsuS06UCAgKwWCykpqYW2teoUSMCAgJYvHgxt9xyCx4eHkUeI0/eNTl5/44ePZrffvuNl156idDQUDw9PYmMjOTixYsl9iPPmTNn8os9ESk9XfMiUkX4+PjQuHFjQkNDCxUuRenQoQMnTpzAzc2Nxo0bF/ipXbt2flx4eDiTJk1i3bp13HnnnSxZsiR/X0hICGPHjmXFihVMmTKFt956C4A6deqQnp5e4FbtvGtayiOnP/Pw8KBFixYkJiYW2le7dm02bNjATz/9xNChQwsVUElJSfzyyy/5r3fs2IGLiwvh4eGA7VqaiRMnMnDgQFq2bImnpyenTp26bD/yfPvtt7Rv377U8SJio+JFRIp08803ExkZyeDBg1m7di1Hjx5l+/btPPnkk3z99df88ccfjB8/nk2bNnHs2DG++OILdu3aRfPmzQGIjo5m7dq1HDlyhD179rBhw4b8fZ07d6ZatWpMnz6dQ4cO8cEHH7B06dKrzqk4/fr1Y9u2bUXuq1u3Lhs2bOD777/n3nvvJTs7O3+fl5cXo0aNYt++ffmFypAhQwgMDASgcePGvPfeexw4cICdO3cyfPhwvL29S/X3PXr0KD///DM333xzqeJF5H9UvIhIkSwWC3Fxcdx4443cf//9hIeHc88993D06FHq1auHq6srp0+fZuTIkYSHhzNkyBAGDBjA7NmzAcjJyWHcuHE0b96c/v3707RpUxYsWABArVq1+Pe//01cXBytW7cmJiaGWbNmXXVOxfnb3/5GXFwcaWlpRe4PDAxkw4YN7N+/n+HDh+c/e6Vx48bceeedDBw4kKioKFq1apXfB4DFixeTmppK+/btGTFiBBMnTqRu3bql+vvGxMQQFRVFaGhoqeJF5H8s5tITzyIi16AhQ4bQvn17nnjiCXunQmZmJk2aNCEmJoZu3brZOx0Rp6OZFxGpEp5//vkrfs5MRTl27BgzZsxQ4SJSRpp5EREREaeimRcRERFxKipeRERExKmoeBERERGnouJFREREnIqKFxEREXEqKl5ERETEqah4EREREaei4kVERESciooXERERcSr/H3Xb12RtQPHHAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, (ax4) = plt.subplots(1, 1, figsize=(6,5))\n", "\n", "\n", "\n", "ax5 = ax4.twiny()\n", "ax4.plot(P_combo*10**(-8), CalcD_P_1000/1000, '-', color='black', label='1000 C')\n", "ax4.plot(P_combo*10**(-8), CalcD_P_1100/1000, '-', color='darkred', label='1100 C')\n", "ax4.plot(P_combo*10**(-8), CalcD_P_1200/1000, '-', color='red', label='1200 C')\n", "ax4.plot(P_combo*10**(-8), CalcD_P_1300/1000, '-', color='orange', label='1300 C')\n", "\n", "\n", "\n", "ax4.set_xlim([P_axis_min, P_axis_max])\n", "ax4.set_ylim([0, 1.5])\n", "\n", "ax5.set_xlim([D_axis_min, D_axis_max])\n", "ax4.set_ylabel('Density (g/cm$^{3}$)')\n", "ax5.set_xlabel('Depth (km)')\n", "ax4.set_xlabel('Pressure (kbar)')\n", "\n", "ax4.legend()" ] }, { "cell_type": "code", "execution_count": 5, "id": "126ee50e-3637-4f05-b07a-f11675f65a7f", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAn8AAAHoCAYAAAA4xYmdAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAABg1UlEQVR4nO3deVxU9f7H8dcICK64IuKeu6KmkFtpaYqSWmq5lksuXUtzy0rzVubtZrd98bqVS5apLWobqZR7aipiLlhXc0EJIjfADRDO74/zY3QETNnOMPN+Ph7zYOZ7zsx85tsI777nfL/HZhiGgYiIiIi4hSJWFyAiIiIiBUfhT0RERMSNKPyJiIiIuBGFPxERERE3ovAnIiIi4kYU/kRERETciMKfiIiIiBtR+BMRERFxIwp/IiIiIm5E4U9EJAeOHTuGzWZjz549t/zcdevW0aBBA9LT0wGYNm0at99+e94WCNxxxx2sWLEiz19XRAo3hT8RKXSGDh2KzWbDZrPh5eVFpUqV6Ny5MwsWLLAHqrx+v549e+bZ6z3zzDNMnTqVIkXy91fw888/z+TJk/OlT0Sk8FL4E5FCqWvXrsTGxnLs2DG+//57OnTowLhx4+jevTtXrlyxurxsbd26lUOHDtGnT598f69u3bqRkJDAmjVr8v29RKTwUPgTkULJ29sbf39/qlSpQosWLXjuuef46quv+P7771m0aJF9v4SEBB577DH8/PwoXbo0HTt25JdffrFvzzjkOnfuXKpVq0bx4sXp06cP586ds2//6KOP+Oqrr+yjjRs2bLA//8iRI3To0IHixYvTrFkztm3bdsO6ly1bRkhICD4+Ptnuc/ToUerUqcPjjz9Oeno6ixYtokyZMnz77bfUr1+f4sWL89BDD3HhwgU++ugjatasSdmyZXnyySdJS0uzv46Hhwf33XcfS5cuvbXOFRGXpvAnIi6jY8eONGvWzH6em2EYdOvWjbi4OMLCwoiIiKBFixbce++9nDlzxv68w4cP89lnn/HNN9+wevVq9uzZw+jRowGYNGkSffv2tY80xsbG0rZtW/tzp06dyqRJk9izZw/16tVjwIABNxx53LRpE8HBwdlu379/P3feeSd9+vRh9uzZ9kPDFy9e5L333mPZsmWsXr2aDRs20Lt3b8LCwggLC+Pjjz9m3rx5fPHFFw6v17JlSzZv3nzrnSkiLsvT6gJERPJSgwYN2Lt3LwDr169n3759xMfH4+3tDcAbb7zBqlWr+OKLL3jssccAuHz5Mh999BFVq1YF4P3336dbt268+eab+Pv7U6xYMZKTk/H398/0fpMmTaJbt24AvPTSSzRu3JjDhw/ToEGDLOs7duwYAQEBWW7btm0b3bt3Z8qUKUyaNMlhW2pqKrNnz6Z27doAPPTQQ3z88cf8+eeflCxZkkaNGtGhQwfWr19Pv3797M+rUqUK0dHRpKen5/s5hiJSOOg3gYi4FMMwsNlsAERERHD+/HnKly9PyZIl7bejR4/y+++/259TvXp1e/ADaNOmDenp6fz2229/+35Nmza1369cuTIA8fHx2e5/6dKlLA/5RkdH06lTJ/75z39mCn4AxYsXtwc/gEqVKlGzZk1Klizp0Hb9excrVoz09HSSk5P/9rOIiHvQyJ+IuJSDBw9Sq1YtANLT06lcubLDOXoZypQpk+1rZITHjJ834uXllel5N5pdW6FCBc6ePZupvWLFigQEBLBs2TKGDx9O6dKls32fjPfKqu369z5z5gzFixenWLFif/tZRMQ9aORPRFzGunXr2LdvHw8++CAALVq0IC4uDk9PT+rUqeNwq1Chgv150dHR/PHHH/bH27Zto0iRItSrVw+AokWLOkykyI3mzZsTFRWVqb1YsWJ8++23+Pj40KVLF5KSkvLk/fbv30+LFi3y5LVExDUo/IlIoZScnExcXBwxMTHs3r2bV155hQceeIDu3bszePBgADp16kSbNm3o2bMna9as4dixY2zdupV//vOf7Nq1y/5aPj4+DBkyhF9++YXNmzczduxY+vbtaz/Hr2bNmuzdu5fffvuNU6dOkZqamuO6u3TpwpYtW7LcVqJECb777js8PT0JDQ3l/PnzOX6fDJs3byYkJCTXryMirkPhT0QKpdWrV1O5cmVq1qxJ165dWb9+Pe+99x5fffUVHh4egHkYNCwsjPbt2zNs2DDq1atH//79OXbsGJUqVbK/Vp06dejduzf33XcfISEhBAYGMmvWLPv2kSNHUr9+fYKDg6lYsSI//fRTjut+5JFHiIqKyvZ8wpIlS/L9999jGAb33XcfFy5cyPF7xcTEsHXrVh599NEcv4aIuB6bYRiG1UWIiFhl2rRprFq1KkeXacupZ555hoSEBObOnZuv7/P000+TkJDAvHnz8vV9RKRw0cifiEgBmzp1KjVq1Miz8wiz4+fnx7/+9a98fQ8RKXw021dEpID5+vry3HPP5fv7PP300/n+HiJS+Oiwr4iIiIgb0WFfERERETei8CciIiLiRhT+RERERNyIwt81Zs2aRa1atfDx8SEoKIjNmzdbXZIlNm3aRI8ePQgICMBms7Fq1SqH7YZhMG3aNAICAihWrBj33HMPBw4csKZYC8yYMYM77riDUqVK4efnR8+ePTOt2ebOfTR79myaNm1K6dKlKV26NG3atOH777+3b3fnvrnejBkzsNlsjB8/3t7mzv0zbdo0bDabwy1joW1w777JEBMTwyOPPEL58uUpXrw4t99+OxEREfbt7txHNWvWzPT9sdlsjB49GnDvvrmewt//W758OePHj2fq1KlERkbSrl07QkNDiY6Otrq0AnfhwgWaNWvGzJkzs9z+2muv8dZbbzFz5kx27tyJv78/nTt3zrPLUTm7jRs3Mnr0aLZv3054eDhXrlwhJCTEYTFed+6jqlWr8uqrr7Jr1y527dpFx44deeCBB+y/ZN25b661c+dO5s2bR9OmTR3a3b1/GjduTGxsrP22b98++zZ375uzZ89y55134uXlxffff09UVBRvvvmmw3Wq3bmPdu7c6fDdCQ8PB6BPnz6Ae/dNJoYYhmEYLVu2NEaNGuXQ1qBBA2Py5MkWVeQcAGPlypX2x+np6Ya/v7/x6quv2tsuX75s+Pr6GnPmzLGgQuvFx8cbgLFx40bDMNRHWSlbtqzx4Ycfqm/+X1JSklG3bl0jPDzcuPvuu41x48YZhqHvzosvvmg0a9Ysy23u3jeGYRjPPvuscdddd2W7XX3kaNy4cUbt2rWN9PR09c11NPIHpKSkEBERken6lyEhIWzdutWiqpzT0aNHiYuLc+grb29v7r77brftq4SEBADKlSsHqI+ulZaWxrJly7hw4QJt2rRR3/y/0aNH061bNzp16uTQrv6BQ4cOERAQQK1atejfvz9HjhwB1DcAX3/9NcHBwfTp0wc/Pz+aN2/OBx98YN+uProqJSWFTz75hGHDhmGz2dQ311H4A06dOkVaWprDtT4BKlWqRFxcnEVVOaeM/lBfmQzDYOLEidx1110EBgYC6iOAffv2UbJkSby9vRk1ahQrV66kUaNG6htg2bJl7N69mxkzZmTa5u7906pVKxYvXsyaNWv44IMPiIuLo23btpw+fdrt+wbgyJEjzJ49m7p167JmzRpGjRrF2LFjWbx4MaDvz7VWrVrFuXPnGDp0KKC+uZ6u8HENm83m8NgwjExtYlJfmcaMGcPevXvZsmVLpm3u3Ef169dnz549nDt3ji+//JIhQ4awceNG+3Z37ZsTJ04wbtw41q5di4+PT7b7uWv/hIaG2u83adKENm3aULt2bT766CNat24NuG/fAKSnpxMcHMwrr7wCQPPmzTlw4ACzZ89m8ODB9v3cuY8yzJ8/n9DQUAICAhza1TcmjfwBFSpUwMPDI1P6j4+Pz/R/Ce4uY+ad+gqefPJJvv76a9avX0/VqlXt7eojKFq0KHXq1CE4OJgZM2bQrFkz3n33Xbfvm4iICOLj4wkKCsLT0xNPT082btzIe++9h6enp70P3LV/rleiRAmaNGnCoUOH3P67A1C5cmUaNWrk0NawYUP7xET1ken48eP88MMPjBgxwt6mvnGk8If5hyooKMg+MyhDeHg4bdu2tagq51SrVi38/f0d+iolJYWNGze6TV8ZhsGYMWNYsWIF69ato1atWg7b1UeZGYZBcnKy2/fNvffey759+9izZ4/9FhwczMMPP8yePXu47bbb3Lp/rpecnMzBgwepXLmy2393AO68885My0r973//o0aNGoB+92RYuHAhfn5+dOvWzd6mvrmORRNNnM6yZcsMLy8vY/78+UZUVJQxfvx4o0SJEsaxY8esLq3AJSUlGZGRkUZkZKQBGG+99ZYRGRlpHD9+3DAMw3j11VcNX19fY8WKFca+ffuMAQMGGJUrVzYSExMtrrxgPP7444avr6+xYcMGIzY21n67ePGifR937qMpU6YYmzZtMo4ePWrs3bvXeO6554wiRYoYa9euNQzDvfsmK9fO9jUM9+6fp556ytiwYYNx5MgRY/v27Ub37t2NUqVK2X8Pu3PfGIZh7Nixw/D09DT+/e9/G4cOHTKWLFliFC9e3Pjkk0/s+7h7H6WlpRnVq1c3nn322Uzb3L1vrqXwd43//ve/Ro0aNYyiRYsaLVq0sC/d4W7Wr19vAJluQ4YMMQzDXE7gxRdfNPz9/Q1vb2+jffv2xr59+6wtugBl1TeAsXDhQvs+7txHw4YNs/87qlixonHvvffag59huHffZOX68OfO/dOvXz+jcuXKhpeXlxEQEGD07t3bOHDggH27O/dNhm+++cYIDAw0vL29jQYNGhjz5s1z2O7ufbRmzRoDMH777bdM29y9b65lMwzDsGTIUUREREQKnM75ExEREXEjCn8iIiIibkThT0RERMSNKPyJiIiIuBGFPxERERE3ovAnIiIi4kYU/q6TnJzMtGnTSE5OtroUp6T+uTH1T/bUNzem/rkx9c+NqX+yp77JTOv8XScxMRFfX18SEhIoXbq01eU4HfXPjal/sqe+uTH1z42pf25M/ZM99U1mGvkTERERcSMKfyIiIiJuxNPqApzNlStXADhx4gS+vr4WV+N8kpKSAIiJiSExMdHiapyP+id76psbU//cmPrnxtQ/2SsMfZOens6ff/5J8+bN8fTM/2imc/6us2XLFtq1a2d1GSIiIuJmduzYwR133JHv76ORv+tUr14dMP8DVK5c2eJqRERExNXFxsbSsmVLKlWqVCDvp/B3nSJFzNMgK1euTNWqVS2uRkRERNxFRgbJ9/cpkHcREREREaeg8CciIiLiRhT+RERERNyIwp+IiIiIG1H4ExEREXEjCn8iIiIibkThT0RERMSNKPyJiIiIuBGnDn+bNm2iR48eBAQEYLPZWLVq1U0/96effsLT05Pbb7893+oTERERKWycOvxduHCBZs2aMXPmzFt6XkJCAoMHD+bee+/Np8pERERECienvrxbaGgooaGht/y8f/zjHwwcOBAPD49bGi0UERERcXVOPfKXEwsXLuT333/nxRdfvKn9k5OTSUxMtN+SkpLyuUIRERER67hU+Dt06BCTJ09myZIleHre3KDmjBkz8PX1td8aNWqUz1WKiIiIWMdlwl9aWhoDBw7kpZdeol69ejf9vClTppCQkGC/RUVF5WOVIiIi4vauXIE9e2DOHBg6FO6+u0Df3qnP+bsVSUlJ7Nq1i8jISMaMGQNAeno6hmHg6enJ2rVr6dixY6bneXt74+3tbX+cmJhYYDWLiIiIGzh1CrZtg61bYft22LkTLlywrByXCX+lS5dm3759Dm2zZs1i3bp1fPHFF9SqVcuiykRERMRtpKVBVNTVsLd1Kxw6lHm/0qWhVSto3Rrq1IEhQwqsRKcOf+fPn+fw4cP2x0ePHmXPnj2UK1eO6tWrM2XKFGJiYli8eDFFihQhMDDQ4fl+fn74+PhkahcRERHJEwkJ8PPPV4Pezz9DVkcRGzWCNm2u3ho0gCL/f/bdyZMFWrJTh79du3bRoUMH++OJEycCMGTIEBYtWkRsbCzR0dFWlSciIiLuxDDMUbytW6+O7B04YLZfq2RJc1SvbVsz6LVuDWXLWlNzFmyGcX3F7u3kyZNUq1aNEydOULVqVavLEREREatcvGien5cxqrdtG5w+nXm/2rWvBr22bSEwEDw8bvptCjp7OPXIn4iIiEiB+eMP+Okn2LLF/Llnj3kO37V8fCA42Ax5bduao3qVKllSbk4p/ImIiIj7MQz49Vcz6GXcjhzJvF/Vqo6jerffDkWLFni5eUnhT0RERFxfSgpERFwNej/9lPkQbpEi0KwZ3HUX3HmnGfaqVbOm3nyk8CciIiKu59w58xy9jLC3Ywdcvuy4T7Fi5mHbu+4yb61bm0uwuDiFPxERESn8Tp50PIS7d2/mWbgVK14NenfdBc2bg5eXNfVaSOFPRERECpeM8/U2bYLNm82wd/x45v3q1nUMe3Xrgs1W8PU6GYU/ERERcW7p6bB/P2zcaAa+TZsgPt5xHw8PcyQvI+jdeSf4+1tTr5NT+BMRERHncuUK/PKLGfY2bjRH986eddzHx8ecgdu+PbRrZy6qXLKkNfUWMgp/IiIiYq3UVNi16+rI3pYtkJTkuE+JEuZo3t13m7fgYPD2tqbeQk7hT0RERArW5cvm7NuMkb1t28yraVzL19c8fJsR9tx0ckZ+UPgTERGR/HXxohnwMsLezz9DcrLjPuXKmYdwM8Je06a3dIk0uXkKfyIiIpK3kpNh+3ZYvx7WrTPvp6Y67lOpkhnyMgJfo0bmIsuS7xT+REREJHcyztnLCHs//ZR5QeUqVeCee66GvXr1tOyKRRT+RERE5NakpcGePVfD3ubNcP684z5+ftCxI3ToYP6sXVthz0ko/ImIiMiNpafDgQNXw97Gjebl065VtqwZ9DLCXsOGCntOSuFPREREHBkGHDpkBr31683bX3857lOqlHn4NiPsNW2qc/YKCYU/ERERMa+N+8MP8OOPZuj74w/H7cWKmYspZ4S9Fi3AUzGiMNJ/NREREXeUmGgevg0PN2+//uq4vWhRaNv2athr2dJsk0JP4U9ERMQdpKaaCyuHh5sjfNu3mxM3MhQpYl41o1MnuPde89JpxYpZV6/kG4U/ERERV2QY8NtvV0f2NmzIfMm0OnXMsNe5sznCV7asJaVKwVL4ExERcRV//mmO6mXcTp503F6+vDmq17mzGfpq1rSkTLGWwp+IiEhhdfEibNp09VDu3r2O2729zUkaGaN7t9+uGbmi8CciIlJoGAb88gusXg1r1sDWrZCS4rhP8+ZXR/buukvn7UkmCn8iIiLO7PRpc2QvI/DFxTlur17dDHudO5uzcitWtKZOKTQU/kRERJxJWhrs3GmGvdWrzRm6hnF1e4kSZsjr0gVCQsxJG7qShtwChT8RERGrxcaao3qrV8PatXD2rOP2Jk2ga1fzdued5rl8Ijmk8CciIlLQUlLgp5+uHsr95RfH7WXKmIdxu3Y1R/iqVLGkTHFNCn8iIiIF4ejRq4dy162D8+evbrPZ4I47ro7u3XGHLp0m+UbfLBERkfxw+bJ5+bSwMDPw/e9/jtv9/K6Gvc6doUIFa+oUt6PwJyIikldiYuC778zbDz+Y6/Bl8PQ0r5WbEfiaNdOae2IJhT8REZGcSkszZ+NmBL49exy3BwRAt24QGmrO0PX1taRMkWsp/ImIiNyKs2fNSRrffWcezj116uo2mw1atTIDX/fu5uielmERJ6PwJyIiciOGAVFRZtj79lvzqhppaVe3lyljzsjt1s08nKtFlsXJKfyJiIhcLyXFnKzx9dfwzTdw/Ljj9kaNro7utW2rmblSqOjbKiIiAnDmDHz/vRn4vv8ekpKubvP2hg4dzLB3331Qq5Z1dYrkksKfiIi4r99/N8Pe11/D5s2Oh3P9/aFHD/PWsaN5WTURF6DwJyIi7iNjdm5G4IuKctzepAncf795Cw7WUizikhT+RETEtV24YK659/XX5oSN+Pir2zw94e67zbDXo4cO54pbUPgTERHX89df5kSNlSvN4Hf58tVtvr7meXv332/Ozi1TxrIyRayg8CciIq4hOtoMeytXmufvpadf3VazJjzwgBn42rUDLy/LyhSxmlOfzLBp0yZ69OhBQEAANpuNVatW3XD/FStW0LlzZypWrEjp0qVp06YNa9asKZhiRUSkYGWsv/fvf0NQENSoAePHm0u0pKdD8+YwfTrs3QtHjsA775gTNxT8xM059cjfhQsXaNasGY8++igPPvjg3+6/adMmOnfuzCuvvEKZMmVYuHAhPXr04Oeff6Z58+YFULGIiOSr9HTYufPqCN///nd1W5EicNdd0KsX9OxpjvaJSCZOHf5CQ0MJDQ296f3feecdh8evvPIKX331Fd9880224S85OZnk5GT746Rr13USERHrpabCpk1m2Fu1CmJirm4rWhQ6dTID3/33g5+fZWWKFBZOHf5yKz09naSkJMqVK5ftPjNmzOCll14qwKpERORvXb4Ma9fCihXmLN2zZ69uK1nSvLpGr14QGgqlS1tXp0gh5NLh78033+TChQv07ds3232mTJnCxIkT7Y9jYmJo1KhRQZQnIiLXunQJ1qyBzz83Z+peeySmQgVzwkavXnDvveDjY12dIoWcy4a/pUuXMm3aNL766iv8bnAYwNvbG29vb/vjxMTEgihPREQALl6E1avNwPftt3D+/NVtVarAgw+atzvvBA8P6+oUcSEuGf6WL1/O8OHD+fzzz+nUqZPV5YiIyLUuXICwMPjiC/juO/NxhurV4aGHzFurVrrChkg+cLnwt3TpUoYNG8bSpUvp1q2b1eWIiAiYI3rffWeO8IWFmYd4M9SoAX36mIGvZUuw2ayrU8QNOHX4O3/+PIcPH7Y/Pnr0KHv27KFcuXJUr16dKVOmEBMTw+LFiwEz+A0ePJh3332X1q1bExcXB0CxYsXw9fW15DOIiLitpCTzUO7nn8P33zteZaNWLTPw9eljrtGnwCdSYJw6/O3atYsOHTrYH2dMzBgyZAiLFi0iNjaW6Oho+/a5c+dy5coVRo8ezejRo+3tGfuLiEg+u3jRHOFbtsz8ec1SWtSufTXwNW+uwCdiEacOf/fccw+GYWS7/fpAt2HDhvwtSEREMktJMZdlWbYMvvrKcdJG3bpXA1+zZgp8Ik7AqcOfiIg4qbQ02LABli411+K7dh2+GjWgf3/zpsAn4nQU/kRE5Oakp8O2beYI3+efw59/Xt3m7w/9+pmBr1UrBT4RJ6bwJyIi2TMMiIw0A9/y5XDNedaUK2fO0O3fH9q31zp8IoWEwp+IiGT266/w6adm6Dt06Gp7qVLmVTb69zevqevlZV2NIpIjCn8iImKKizPD3iefQETE1XYfH+jRwwx8992nS6uJFHIKfyIi7uzCBVi1ygx8a9ea5/UBeHpCly4wcKAZ/EqVsrRMEck7Cn8iIu7myhX48Ucz8K1c6Xh5tdat4ZFHoG9fqFjRuhpFJN8o/ImIuIOMiRuffGIuz/L/V0ACzMWXBw2Chx+GOnWsq1FECoTCn4iIKzt2zJy48ckncPDg1fby5c1z+B55REuziLgZhT8REVeTkACffQYffwybN19t9/GBBx4wA1+XLpqpK+KmFP5ERFxBWhqsWweLFplX3Lh82Wy32aBjRzPw9e4NpUtbWqaIWE/hT0SkMDt0CD76CBYvhhMnrrY3bgyDB5uzdatWta4+EXE6Cn8iIoVNYqJ5ebVFi2DLlqvtZcuaYW/oUAgK0nl8IpIlhT8RkcIgPR02bDAD35dfwsWLZnuRIub5e48+aq7HpwWYReRvKPyJiDizI0fMw7offQTHj19tb9DADHyPPAIBAdbVJyKFjsKfiIizOX8evvgCFi6ETZuutvv6woAB5mHdli11WFdEckThT0TEGRiGeT3dDz4wF2FOSjLbbTYICTED3wMPQLFilpYpIoWfwp+IiJUSEmDJEjP07dlztb1uXfOw7qBBmq0rInlK4U9EpKAZBvz0E3z4obkY86VLZru3Nzz0EIwcCe3b67CuiOQLhT8RkYJy6pS5Ht+HHzpeai0w0Ax8jzwC5cpZV5+IuAWFPxGR/JSeDuvXm4d1V66ElBSzvXhx89q6I0fq2roiUqAU/kRE8kNsrDlbd/58c7mWDEFBZuAbMECXWhMRSyj8iYjklfR0+PFHmDULvvnGvN4umCHv4YfN0Ne8ubU1iojbU/gTEcmts2fNK2/Mnm1eazdD27Zm4OvTB0qUsKw8EZFrKfyJiORURIQ5yrd06dUZu6VKwZAhMGoUNG5sbX0iIllQ+BMRuRWXLpnLs8yaBTt2XG1v2hSeeMI8vFuypHX1iYj8DYU/EZGb8fvvMGcOLFgAZ86YbUWLmuvyPfGEeYhXM3ZFpBBQ+BMRyU5aGoSFmaN8q1dfba9RwzysO2wY+PlZV5+ISA4o/ImIXC8+3lyiZc4ciI4222w26NrVHOULDQUPD2trFBHJIYU/EZEMu3bBe+/B8uVXF2MuVw6GD4d//ANq17a2PhGRPKDwJyLuLTUVVqwwQ9/WrVfbW7UyR/n69IFixayrT0Qkjyn8iYh7OnUK5s0zz+eLiTHbvLygXz8YOxbuuMPa+kRE8onCn4i4l7174d13YckSSE422/z84PHHzUO7lStbW5+ISD5T+BMR15eWBl9/bYa+jRuvtrdoAePGmaN93t7W1SciUoAU/kTEdZ0/DwsXwjvvwJEjZpuHBzz4oHloV2vziYgbUvgTEdcTEwPvvw9z58K5c2Zb2bLmYd0nnoBq1SwtT0TESgp/IuI69uyBN9+EZcvgyhWzrU4dmDDBvN5uiRKWlici4gwU/kSkcEtPh++/N0Pf+vVX29u3h4kToUcPKFLEuvpERJyMwp+IFE6XL8PixfD22/Drr2abhwf07WuGvuBga+sTEXFSTv2/w5s2baJHjx4EBARgs9lYtWrV3z5n48aNBAUF4ePjw2233cacOXPyv1ARKTjnzsErr5jX1/3HP8zgV7o0TJoER4/Cp58q+ImI3IBTh78LFy7QrFkzZs6ceVP7Hz16lPvuu4927doRGRnJc889x9ixY/nyyy/zuVIRyXcxMWbAq1YNpk41r79bvTq89RacOAGvv66JHCIiN8GpD/uGhoYSGhp60/vPmTOH6tWr88477wDQsGFDdu3axRtvvMGDDz6YT1WKSL46eNAMdp98Yl6KDaBJE3jmGXN9Pi8va+sTESlknDr83apt27YREhLi0NalSxfmz59PamoqXln8kUhOTiY5Y5V/ICkpKd/rFJGbsHUr/Oc/5uLMGe6+2wx9oaFan09EJIec+rDvrYqLi6NSpUoObZUqVeLKlSucOnUqy+fMmDEDX19f+61Ro0YFUaqIZMUw4NtvoV07uPNOM/jZbNCrF2zfDhs2wH33KfiJiOSCS4U/ANt1fxQMw8iyPcOUKVNISEiw36KiovK9RhG5TloafPYZ3H67uTTLli3m4dzhw83DvitWQKtWVlcpIuISXOqwr7+/P3FxcQ5t8fHxeHp6Ur58+Syf4+3tjfc11/RMTEzM1xpF5Bqpqebs3Bkz4LffzLaSJeHxx2H8eAgIsLQ8ERFX5FLhr02bNnzzzTcObWvXriU4ODjL8/1ExCKXL8OiReY5fceOmW1ly8K4cfDkk1CunJXViYi4NKcOf+fPn+fw4cP2x0ePHmXPnj2UK1eO6tWrM2XKFGJiYli8eDEAo0aNYubMmUycOJGRI0eybds25s+fz9KlS636CCJyrQsXYN48c/ZubKzZ5ucHTz1ljvaVKmVtfSIibsCpw9+uXbvo0KGD/fHEiRMBGDJkCIsWLSI2Npbo6Gj79lq1ahEWFsaECRP473//S0BAAO+9956WeRGxWkIC/Pe/5tU4MiZfVa1qztwdPhyKF7e2PhERN2IzMmZECAAnT56kWrVqnDhxgqpVq1pdjkjhlpAA77xjhr6EBLPttttg8mQYPBiuOd9WRMRdFXT2cOqRPxEppBIS4N13zdB37pzZ1rAhPPcc9O8PnvrVIyJiFf0GFpG8k5gI771nXnLt7FmzrWFDePFF6NMHirjc6lIiIoWOwp+I5F5iIrz/Prz55tXQ16DB1dDn4WFtfSIiYqfwJyI5l5R0NfSdOWO2NWgAL7wAffsq9ImIOCGFPxG5dRcumKHvjTfg9GmzrX59M/T166fQJyLixBT+ROTmJSeb6/T9+9/w559mW926ZugbMEChT0SkEFD4E5G/d+UKfPIJTJsGx4+bbbfdZj4eMECzd0VEChH9xhaR7BkGfPklPP88/Pqr2Va5sjnSN3w46LKJIiKFjsKfiGRmGLB2LUydChERZlu5cjBlCoweDcWKWVufiIjkmMKfiDj66SdzMeZNm8zHJUvCxInmzdfX2tpERCTXFP5ExLR/v3nZte++Mx97e5ujfJMnQ8WK1tYmIiJ5RuFPxN3FxJiLMS9cCOnp5ozdYcPM8/yqVbO6OhERyWMKfyLuKjERXnvNvBTbpUtm20MPmcu41KtnbW0iIpJvFP5E3E1qqrlW30svwV9/mW1t25oLNrdpY21tIiKS7xT+RNyFYcDKleY5fIcOmW316sGrr0LPnmCzWVqeiIgUDIU/EXewdSs8/bT5E8wJHNOmwciRWqtPRMTN5En4S01NJS4ujosXL1KxYkXKlSuXFy8rIrl15Ag884y5UDOY6/M99ZQZBEuXtrY2ERGxRJGcPvH8+fPMnTuXe+65B19fX2rWrEmjRo2oWLEiNWrUYOTIkezcuTMvaxWRm5WUZC7I3LChGfyKFIERI+DwYfjXvxT8RETcWI7C39tvv03NmjX54IMP6NixIytWrGDPnj389ttvbNu2jRdffJErV67QuXNnunbtyqGM84tEJH+lp8OiRVfP5UtJgU6d4Jdf4IMPICDA6gpFRMRiOTrsu3XrVtavX0+TJk2y3N6yZUuGDRvGnDlzmD9/Phs3bqRu3bq5KlRE/sbWrTBuHOzaZT6uUwfefBN69NBkDhERsctR+Pv8889vaj9vb2+eeOKJnLyFiNysEyfg2Wdh6VLzcalS5gLNY8eaV+kQERG5hmb7ihRWFy+aa/O9+qq5SLPNZl6Z49//hkqVrK5ORESc1C2Hv7Nnz2IYBuXKleOvv/5i06ZN1K9fn8DAwPyoT0SuZxjw2WfmjN0TJ8y2du3gnXegRQtLSxMREed3SxM+PvzwQ4KDgwkKCmL27Nn06tWLH3/8kf79+zNv3rz8qlFEMhw8aE7g6N/fDH7Vq8Py5bBxo4KfiIjclFsa+Xv//fc5cOAAFy9epHr16hw9epSKFSuSmJhI+/bteeyxx/KrThH3duECvPyyOYEjNRV8fMylXJ5+2ly7T0RE5CbdUvjz8PDAx8cHHx8f6tSpQ8WKFQEoXbo0Ns0mFMl7hgFffWXO4o2ONtu6d4f33oNataytTURECqVbOuzr6enJ5cuXAdi4caO9PSkpKW+rEhHz6hzdu0OvXmbwq1HDDILffKPgJyIiOXZL4W/dunV4///SEb6+vvb2S5cuMX/+/LytTMRdXb4M06dDo0YQFmZee/e55yAqCu6/3+rqRESkkLulw74lS5bMst3Pzw8/P788KUjEra1eDWPGwO+/m487dYKZM6F+fWvrEhERl5En6/xdvnyZvXv3Eh8fT3p6usO2+zVSIfL3/vjDXJT5yy/NxwEB8NZb0Levrs4hIiJ5Ktfhb/Xq1QwePJhTp05l2maz2UhLS8vtW4i4rvR0mD/fnLWbkAAeHubkjmnTzCt1iIiI5LFbOucvK2PGjKFPnz7ExsaSnp7ucFPwE7mB//0POnaExx4zg98dd0BEhLmci4KfiIjkk1yHv/j4eCZOnEglXU5K5OakpsKMGdC0qbk4c/Hi5iHebdugWTOrqxMREReX68O+Dz30EBs2bKB27dp5UY+Ia9u1C4YPh717zcchITBnjpZuERGRApPr8Ddz5kz69OnD5s2badKkCV5eXg7bx44dm9u3ECn8LlyAF1+Et982z/MrV868Fu8jj2hCh4iIFKhch79PP/2UNWvWUKxYMTZs2OBwpQ+bzabwJxIeDv/4Bxw9aj4eONAMgVoeSURELJDr8PfPf/6T6dOnM3nyZIoUyfUphCKu49w5mDABFi0yH1erZh7ive8+K6sSERE3l+u0lpKSQr9+/RT8RK61di0EBprBz2Yz1/A7cEDBT0RELJfrxDZkyBCWL1+eF7WIFH7nz8MTT0CXLhATA3XqwJYt8O67Wr5FREScQq4P+6alpfHaa6+xZs0amjZtmmnCx1tvvZWr1581axavv/46sbGxNG7cmHfeeYd27dplu/+SJUt47bXXOHToEL6+vnTt2pU33niD8uXL56oOkb+1eTMMHQpHjpiPx4yBV1+FEiUsLUtERORauR7527dvH82bN6dIkSLs37+fyMhIh1tuLF++nPHjxzN16lQiIyNp164doaGhREdHZ7n/li1bGDx4MMOHD+fAgQN8/vnn7Ny5kxEjRuSqDpEbunwZJk2Cu+82g1+1avDDD/D++wp+IiLidGyGYRhWF5GdVq1a0aJFC2bPnm1va9iwIT179mTGjBmZ9n/jjTeYPXs2v//+u73t/fff57XXXuPEiRM39Z4nT56kWrVqnDhxgqpVq+b+Q4hr27ULBg+GgwfNx48+as7k9fW1ti4RESk0Cjp75Hrkb8aMGSxYsCBT+4IFC/jPf/6T49dNSUkhIiKCkJAQh/aQkBC2bt2a5XPatm3LyZMnCQsLwzAM/vzzT7744gu6deuW7fskJyeTmJhovyUlJeW4ZnEjKSnmun2tW5vBr1Il+PprWLBAwU9ERJxarsPf3LlzadCgQab2xo0bM2fOnBy/7qlTp0hLS8t02bhKlSoRFxeX5XPatm3LkiVL6NevH0WLFsXf358yZcrw/vvvZ/s+M2bMwNfX135r1KhRjmsWN7F/vxn6pk+HtDTo29ecydujh9WViYiI/K1ch7+4uDgqV66cqb1ixYrExsbm9uUdFo0GMAwjU1uGqKgoxo4dywsvvEBERASrV6/m6NGjjBo1KtvXnzJlCgkJCfZbVFRUrmsWF2UYMHMmBAdDZKR5lY5ly2D5ctCEIhERKSRyPdu3WrVq/PTTT9S67tqkP/30EwEBATl+3QoVKuDh4ZFplC8+Pj7TaGCGGTNmcOedd/L0008D0LRpU0qUKEG7du14+eWXswyp3t7eeHt72x8nJibmuGZxYadOwbBh8M035uP77oMPP4QsvlMiIiLOLNcjfyNGjGD8+PEsXLiQ48ePc/z4cRYsWMCECRMYOXJkjl+3aNGiBAUFER4e7tAeHh5O27Zts3zOxYsXMy027eHhAZgjhiI5sn49NGtmBr+iRc01+779VsFPREQKpVyP/D3zzDOcOXOGJ554gpSUFAB8fHx49tlnmTJlSq5ee+LEiQwaNIjg4GDatGnDvHnziI6Oth/GnTJlCjExMSxevBiAHj16MHLkSGbPnk2XLl2IjY1l/PjxtGzZMlejkOKmUlPhpZfglVfMQ77165uHeW+/3erKREREcizX4c9ms/Gf//yH559/noMHD1KsWDHq1q3rcCg1p/r168fp06eZPn06sbGxBAYGEhYWRo0aNQCIjY11WPNv6NChJCUlMXPmTJ566inKlClDx44dczXrWNzUsWMwcCBs22Y+Hj7cHPHTun0iIlLI5Xidv+eee46ePXvSsmXLvK7JUlrnT1i+HB57DBIToXRpmDcP+vWzuioREXFRhWadv9jYWLp3707lypV57LHH+O6770hOTs7L2kQK1oUL5ghf//5m8GvdGn75RcFPRERcSo7D38KFC/nzzz/57LPPKFOmDE899RQVKlSgd+/eLFq0iFOnTuVlnSL5a88eCAoyF2m22WDqVNi0CWrWtLoyERGRPJWr2b42m4127drx2muv8euvv7Jjxw5at27NBx98QJUqVWjfvj1vvPEGMTExeVWvSN4yDPOwbuvW8NtvEBAAP/4IL78MXl5WVyciIpLncj3h41oNGzakYcOGPPPMM8THx/PNN9/w9ddfAzBp0qS8fCuR3LtwAR5/HD7+2HzcvTssXAgVKlhbl4iISD7K0/B3LT8/P4YPH87w4cPz6y1Ecu633+DBB83LshUpAjNmwKRJ5n0REREXluvwN3HixCzbbTYbPj4+1K1bl/vvv59y5crl9q1E8sby5TBiBJw/D/7+5tp9d99tdVUiIiIFItfhLzIykt27d5OWlkb9+vUxDINDhw7h4eFBgwYNmDVrFhMnTmTLli00atQoL2oWyZnUVHjmGXjnHfPxPffA0qVmABQREXETuT7G9cADD9CpUyf++OMPIiIi2L17NzExMXTu3JkBAwYQExND+/btmTBhQl7UK5Iz8fHQufPV4Dd5MoSHK/iJiIjbyfEizxmqVKlCeHh4plG9AwcOEBISQkxMDLt37yYkJKRQLP+iRZ5d0I4d5vl9J09CqVLw0UfQq5fVVYmIiACFaJHnDAkJCcTHx2dq/+uvv0hMTASgTJky9uv+ihSo+fOhXTsz+NWvbwZBBT8REXFjeXLYd9iwYaxcuZKTJ08SExPDypUrGT58OD179gRgx44d1KtXL7dvJXLzUlLMZVxGjDDv9+xpBr8GDayuTERExFK5nvAxd+5cJkyYQP/+/bly5Yr5op6eDBkyhLfffhuABg0a8OGHH+b2rURuzqlT8NBDsHGjebWOf/0LpkzRMi4iIiLkwTl/Gc6fP8+RI0cwDIPatWtTsmTJvHjZAqdz/gq5ffvg/vvh2DHz/L6lS6FbN6urEhERyVahOOcvOjo6U1vJkiVp2rQpzZo1yxT8dHk3KRBffQVt25rBr3Zt2L5dwU9EROQ6OQp/d9xxByNHjmTHjh3Z7pOQkMAHH3xAYGAgK1asyHGBIn/LMMwrdPTqZS7c3LEj/PwzaF1JERGRTHJ0zt/Bgwd55ZVX6Nq1K15eXgQHBxMQEICPjw9nz54lKiqKAwcOEBwczOuvv05oaGhe1y1iSk6GkSOvXp939Gh4+23w8rK2LhERESeVq3P+Ll++TFhYGJs3b+bYsWNcunSJChUq0Lx5c7p06UJgYGBe1logdM5fIXLmjDnat2kTeHjAzJkwapTVVYmIiNySgs4euZrt6+PjQ+/evendu3de1SNycw4fhvvug0OHoHRp+OIL8woeIiIickO5XupFpMBt2WKu23f6NNSoAd9+C4VwlFlERMQKWvhMCpelS+Hee83gd8cd5oxeBT8REZGbpvAnhYNhwBtvwMCB5hU7eveGDRvA39/qykRERAoVhT9xfunpMHEiPP20+XjCBPj8cyhe3Nq6RERECiGd8yfOLTkZBg+Gzz4zH7/5phkERUREJEdyPfI3dOhQNm3alBe1iDhKSICuXc3g5+UFn36q4CciIpJLuQ5/SUlJhISEULduXV555RVdyk3yxh9/QLt25nl9pUrB99/DgAFWVyUiIlLo5Tr8ffnll8TExDBmzBg+//xzatasSWhoKF988QWpqal5UaO4m99/h7vugn37zAkdmzaZM3xFREQk1/Jkwkf58uUZN24ckZGR7Nixgzp16jBo0CACAgKYMGEChw4dyou3EXewf78Z/I4ehdq1YetWuP12q6sSERFxGXk62zc2Npa1a9eydu1aPDw8uO+++zhw4ACNGjXi7bffzsu3Elf088/Qvj3ExUGTJuZizrVqWV2ViIiIS8l1+EtNTeXLL7+ke/fu1KhRg88//5wJEyYQGxvLRx99xNq1a/n444+ZPn16XtQrruqHH8xDu2fPQps2sHGj1vATERHJB7le6qVy5cqkp6czYMAAduzYwe1ZHKLr0qULZcqUye1biatauRL69zcXb+7c2XxcooTVVYmIiLikXIe/cePG8dRTT1H8ugV3DcPgxIkTVK9enbJly3L06NHcvpW4ouXL4eGHIS0NHnwQliwBb2+rqxIREXFZuT7sO23aNM6fP5+p/cyZM9TS+VpyI598Yl6uLS3NXMh52TIFPxERkXyW6/BnGEaW7efPn8fHxye3Ly+uauFCM/Clp8OIEeZjT11wRkREJL/l+K/txP+/0oLNZuOFF15wOOyblpbGzz//nOX5fyLMmwf/+Id5//HHYeZMKKLLTIuIiBSEHIe/yMhIwBz527dvH0WLFrVvK1q0KM2aNWPSpEm5r1Bcy8yZ8OST5v2xY+Gdd8Bms7QkERERd5Lj8Ld+/XoAHn30Ud59911Kly6dZ0WJi7o2+D31FLz+uoKfiIhIAcv1SVYLFy7MizrE1c2dezX4TZ4Mr7yi4CciImKBHIW/iRMn8q9//YsSJUrYz/3LzltvvZWjwsSFLFgAo0aZ9ydNUvATERGxUI7CX2RkJKmpqfb72bHpD7x88ok5mxfMc/xee03BT0RExEI5Cn8Z5/tdfz8/zJo1i9dff53Y2FgaN27MO++8Q7t27bLdPzk5menTp/PJJ58QFxdH1apVmTp1KsOGDcvXOiULn30GQ4aAYZizejW5Q0RExHK5Pufv0qVLGIZhX+rl+PHjrFy5kkaNGhESEpKr116+fDnjx49n1qxZ3HnnncydO5fQ0FCioqKoXr16ls/p27cvf/75J/Pnz6dOnTrEx8dz5cqVXNUhObBypbmAc3o6DB9uTvZQ8BMREbGczchuleabFBISQu/evRk1ahTnzp2jfv36FC1alFOnTvHWW2/x+OOP5/i1W7VqRYsWLZg9e7a9rWHDhvTs2ZMZM2Zk2n/16tX079+fI0eOUK5cuRy958mTJ6lWrRonTpygatWqOa7drYWHQ/fu5rV6Bw82F3DWOn4iIiJZKujskeu/yLt377Yfhv3iiy/w9/fn+PHjLF68mPfeey/Hr5uSkkJERESm0cOQkBC2bt2a5XO+/vprgoODee2116hSpQr16tVj0qRJXLp0Kdv3SU5OJjEx0X5LSkrKcc0C/Pwz9OplBr+HHjIneyj4iYiIOI1cH/a9ePEipUqVAmDt2rX07t2bIkWK0Lp1a44fP57j1z116hRpaWlUqlTJob1SpUrExcVl+ZwjR46wZcsWfHx8WLlyJadOneKJJ57gzJkzLFiwIMvnzJgxg5deeinHdco19u+H0FC4cAFCQszJHh4eVlclIiIi18j1kEydOnVYtWoVJ06cYM2aNfaRuvj4+DxZ+Pn6GcOGYWQ7izg9PR2bzcaSJUto2bIl9913H2+99RaLFi3KdvRvypQpJCQk2G9RUVG5rtktHT1qBr6zZ6F1a1ixAry9ra5KRERErpPr8PfCCy8wadIkatasSatWrWjTpg1gjgI2b948x69boUIFPDw8Mo3yxcfHZxoNzFC5cmWqVKmCr6+vva1hw4YYhsHJkyezfI63tzelS5e23zJGMeUWxMVB584QGwuBgfDdd1CihNVViYiISBZyHf4eeughoqOj2bVrF6tXr7a333vvvbz99ts5ft2iRYsSFBREeHi4Q3t4eDht27bN8jl33nknf/zxB+fPn7e3/e9//6NIkSKavJFfzp2DLl3g99+hVi1YuxZyONlGRERE8l+enInv7+9P8+bNKXLNif0tW7akQYMGuXrdiRMn8uGHH7JgwQIOHjzIhAkTiI6OZtT/Xy1iypQpDB482L7/wIEDKV++PI8++ihRUVFs2rSJp59+mmHDhlGsWLFc1SJZSE42J3fs3Qv+/uYs38qVra5KREREbiDXEz4AfvzxR3788Ufi4+NJT0932JbdRIub0a9fP06fPs306dOJjY0lMDCQsLAwatSoAUBsbCzR0dH2/UuWLEl4eDhPPvkkwcHBlC9fnr59+/Lyyy/nuAbJRno6PPoobNgApUrB6tVQu7bVVYmIiMjfyPU6fy+99BLTp08nODiYypUrZ5qMsXLlylwVWNC0zt9NmjwZ/vMf8PSEsDDznD8RERG5ZQWdPXI98jdnzhwWLVrEoEGD8qIeKQxmzzaDH8CHHyr4iYiIFCK5PucvJSUl2wkY4oK+/hrGjDHvT59uXrtXRERECo1ch78RI0bw6aef5kUt4ux27ID+/c3z/UaMgH/+0+qKRERE5Bbl+rDv5cuXmTdvHj/88ANNmzbFy8vLYftbb72V27cQZxAdDT16wKVL5lU8Zs2CbBbbFhEREeeV6/C3d+9ebr/9dgD279/vsC27K3FIIXP+PNx/P8THQ7NmsHw5XBfyRUREpHDIdfhbv359XtQhzio9HQYNgl9+AT8/85w/XQVFRESk0MqTRZ7FhT3/PKxaBUWLmj+rV7e6IhEREcmFPAl/mzdv5pFHHqFNmzbExMQA8PHHH7Nly5a8eHmxypIl8Mor5v358+H/r9ssIiIihVeuw9+XX35Jly5dKFasGJGRkSQnJwOQlJTEKxnBQQqf7dth+HDz/uTJ8Mgj1tYjIiIieSLX4e/ll19mzpw5fPDBBw4zfdu2bcvu3btz+/JihZMnoWdP89q9DzwA//631RWJiIhIHsl1+Pvtt99o3759pvbSpUtz7ty53L68FLTkZHjoIfjzT2jaFD75BIro1FARERFXkeu/6pUrV+bw4cOZ2rds2cJtt92W25eXgjZ+PPz8M5QtCytXQsmSVlckIiIieSjX4e8f//gH48aN4+eff8Zms/HHH3+wZMkSJk2axBNPPJEXNUpBWbQI5swxF29esgQU3kVERFxOrtf5e+aZZ0hISKBDhw5cvnyZ9u3b4+3tzaRJkxiTcQ1YcX67d8OoUeb9adPMq3iIiIiIy7EZhmHkxQtdvHiRqKgo0tPTadSoESUL6eHCkydPUq1aNU6cOEHVqlWtLqdgnD4NwcFw7Bh07w5ffaXz/ERERApIQWePXI38paens2jRIlasWMGxY8ew2WzUqlWLhx56iEGDBunyboVBWho8/LAZ/GrXho8/VvATERFxYTn+K28YBvfffz8jRowgJiaGJk2a0LhxY44fP87QoUPp1atXXtYp+WX6dFizBooVgxUroEwZqysSERGRfJTjkb9FixaxadMmfvzxRzp06OCwbd26dfTs2ZPFixczePDgXBcp+eTHH+Ff/zLvf/CBubSLiIiIuLQcj/wtXbqU5557LlPwA+jYsSOTJ09myZIluSpO8lF8vHnVDsOAkSPNQ78iIiLi8nIc/vbu3UvXrl2z3R4aGsovv/yS05eX/JSeDoMHQ1wcNG4M77xjdUUiIiJSQHIc/s6cOUOlSpWy3V6pUiXOnj2b05eX/PTGG1fP81u+HIoXt7oiERERKSA5Dn9paWl4emZ/yqCHhwdXrlzJ6ctLftm+HaZONe+/95458iciIiJuI8cTPgzDYOjQoXh7e2e5PTk5OcdFST45dw4GDIArV6BfPxg+3OqKREREpIDlOPwNGTLkb/fRTF8nYhgwYoS5nt9tt8HcueZl3ERERMSt5Dj8LVy4MC/rkPy2cCF8+SV4ecGyZeDra3VFIiIiYgFdysEdHD0K48aZ919+Ge64w9p6RERExDIKf64uLQ2GDoXz5+Guu+Cpp6yuSERERCyk8Ofq3nkHNm2CEiXgo4/Aw8PqikRERMRCCn+ubP9+eO458/7bb5sTPURERMStKfy5qpQUGDTI/NmtmznTV0RERNyewp+rmj4d9uyB8uXhww+1rIuIiIgACn+uaft2mDHDvD9nDvj7W1uPiIiIOA2FP1dz+bI5uzc9HR55BB56yOqKRERExIko/Lmaf/8bfvvNHO177z2rqxEREREno/DnSvbtg1dfNe/PnAlly1pbj4iIiDgdhT9XkZZmzui9cgV69oTeva2uSERERJyQwp+rmDkTduyA0qXN+5rdKyIiIllQ+HMFx47B1Knm/ddegypVLC1HREREnJfCX2FnGDBqFFy4AO3awciRVlckIiIiTkzhr7D79FNYswa8veGDD6CI/pOKiIhI9pw+KcyaNYtatWrh4+NDUFAQmzdvvqnn/fTTT3h6enL77bfnb4FWOn0axo8377/wAtSvb2k5IiIi4vycOvwtX76c8ePHM3XqVCIjI2nXrh2hoaFER0ff8HkJCQkMHjyYe++9t4AqtcjUqXDqFAQGwtNPW12NiIiIFAJOHf7eeusthg8fzogRI2jYsCHvvPMO1apVY/bs2Td83j/+8Q8GDhxImzZtCqhSC0REwLx55v1Zs8DLy9p6REREpFBw2vCXkpJCREQEISEhDu0hISFs3bo12+ctXLiQ33//nRdffPGm3ic5OZnExET7LSkpKVd1F4j0dBgzxpzsMXCgOdFDRERE5CY4bfg7deoUaWlpVKpUyaG9UqVKxMXFZfmcQ4cOMXnyZJYsWYKnp+dNvc+MGTPw9fW13xo1apTr2vPdxx/D9u1QsiS8/rrV1YiIiEgh4rThL4PtusWKDcPI1AaQlpbGwIEDeemll6hXr95Nv/6UKVNISEiw36KionJdc75KSIBnnjHvv/ACBARYW4+IiIgUKjc3PGaBChUq4OHhkWmULz4+PtNoIEBSUhK7du0iMjKSMWPGAJCeno5hGHh6erJ27Vo6duyY6Xne3t54e3vbHycmJubxJ8lj06ZBfLw5s3fcOKurERERkULGaUf+ihYtSlBQEOHh4Q7t4eHhtG3bNtP+pUuXZt++fezZs8d+GzVqFPXr12fPnj20atWqoErPP/v3w/vvm/ffew+KFrW2HhERESl0nHbkD2DixIkMGjSI4OBg2rRpw7x584iOjmbUqFGAecg2JiaGxYsXU6RIEQIDAx2e7+fnh4+PT6b2QskwYOxYSEuD3r3huokwIiIiIjfDqcNfv379OH36NNOnTyc2NpbAwEDCwsKoUaMGALGxsX+75p/L+PxzWL8efHzgzTetrkZEREQKKZthGIbVRTiTkydPUq1aNU6cOEHVqlWtLseUnGye43f8uHnO300uYyMiIiLOr6Czh9Oe8yfX+O9/zeBXpYqu5CEiIiK5ovDn7M6ehZdfNu//619QvLi19YiIiEihpvDn7GbMMANgYCAMHmx1NSIiIlLIKfw5s+PHzSVdAP7zH/DwsLYeERERKfQU/pzZ88+bkz06dIDQUKurEREREReg8OesIiPhk0/M+6+9Bllc0k5ERETkVin8OatnnzUXdh4wAIKDra5GREREXITCnzNauxbCw8HLC/79b6urERERERei8Ods0tPhmWfM+2PGQK1a1tYjIiIiLkXhz9l8/jn88gv4+sLUqVZXIyIiIi5G4c+ZpKXB9Onm/YkToXx5a+sRERERl6Pw50y++AKioqBMGRg71upqRERExAUp/DmLa0f9JkwwA6CIiIhIHlP4cxbXjvqNG2d1NSIiIuKiFP6cQVoavPSSeX/iRHOyh4iIiEg+UPhzBp9/DgcP6lw/ERERyXcKf1a7foavRv1EREQkHyn8WS1j1K9sWY36iYiISL5T+LOSRv1ERESkgCn8WenaUb8nn7S6GhEREXEDCn9WMQx45RXz/oQJGvUTERGRAqHwZ5W1a2HfPihRAsaMsboaERERcRMKf1Z5/XXz58iR5mFfERERkQKg8GeF3bvhxx/BwwPGj7e6GhEREXEjCn9WeOMN82e/flCjhrW1iIiIiFtR+Ctox4/DZ5+Z959+2tpaRERExO0o/BW0t9821/fr1Aluv93qakRERMTNKPwVpLNn4cMPzfsa9RMRERELKPwVpNmz4cIFaNYMOne2uhoRERFxQwp/BeXyZXjvPfP+00+DzWZtPSIiIuKWFP4KyiefwJ9/QrVq0Lev1dWIiIiIm1L4Kwjp6fDmm+b9CRPAy8vaekRERMRtKfwVhLVr4ddfzev3jhhhdTUiIiLixhT+CsKsWebPRx+FUqWsrUVERETcmsJffjt2DL791rz/+OOWliIiIiKi8Jff5swBwzCXdqlXz+pqRERExM0p/OWny5evLuo8erS1tYiIiIig8Je/Pv8cTp82l3fp1s3qakREREQU/vLVf/9r/hw1Cjw9ra1FREREhEIQ/mbNmkWtWrXw8fEhKCiIzZs3Z7vvihUr6Ny5MxUrVqR06dK0adOGNWvWFGC114iIgJ9/Ntf00/IuIiIi4iScOvwtX76c8ePHM3XqVCIjI2nXrh2hoaFER0dnuf+mTZvo3LkzYWFhRERE0KFDB3r06EFkZGQBV87V5V369AE/v4J/fxEREZEs2AzDMKwuIjutWrWiRYsWzJ49297WsGFDevbsyYwZM27qNRo3bky/fv144YUXbmr/kydPUq1aNU6cOEHVqlVzVDdnzkCVKuaEj59+grZtc/Y6IiIi4vLyJHvcAqcd+UtJSSEiIoKQkBCH9pCQELZu3XpTr5Genk5SUhLlypXLdp/k5GQSExPtt6SkpFzVDcDChWbwa9YM2rTJ/euJiIiI5BGnDX+nTp0iLS2NSpUqObRXqlSJuLi4m3qNN998kwsXLtC3b99s95kxYwa+vr72W6NGjXJVN+npkDFSOXo02Gy5ez0RERGRPOS04S+D7brwZBhGprasLF26lGnTprF8+XL8bnDO3ZQpU0hISLDfoqKiclfw2rXw++/mdXwHDszda4mIiIjkMaddf6RChQp4eHhkGuWLj4/PNBp4veXLlzN8+HA+//xzOnXqdMN9vb298fb2tj9OTEzMedEAc+eaP4cMgRIlcvdaIiIiInnMaUf+ihYtSlBQEOHh4Q7t4eHhtL3BBIqlS5cydOhQPv30U7oV9MLK8fFXr+M7cmTBvreIiIjITXDakT+AiRMnMmjQIIKDg2nTpg3z5s0jOjqaUaNGAeYh25iYGBYvXgyYwW/w4MG8++67tG7d2j5qWKxYMXx9ffO/4CVL4MoVuOMOCAzM//cTERERuUVOHf769evH6dOnmT59OrGxsQQGBhIWFkaNGjUAiI2NdVjzb+7cuVy5coXRo0cz+ppr6Q4ZMoRFixblb7GGYc7yBXj00fx9LxEREZEccup1/qyQ47V2IiIgOBi8vSE2FsqWzb8iRURExGVonb/CKmPUr1cvBT8RERFxWgp/eeHyZfj0U/O+DvmKiIiIE1P4ywtffw1nz0LVqnDvvVZXIyIiIpIthb+8kHHId8gQ8PCwthYRERGRG1D4y62YGPOqHgBDh1paioiIiMjfUfjLrcWLzev5tmsHdepYXY2IiIjIDSn85YbW9hMREZFCRuEvN7ZuhUOHzGv49uljdTUiIiIif0vhLzcyRv369IGSJa2tRUREROQmKPzl1MWLsHy5eV+HfEVERKSQUPjLqW+/hfPnoVYtc7KHiIiISCGg8JdTS5eaP/v3B5vN2lpEREREbpLCX04kJEBYmHl/wABraxERERG5BQp/ObFyJaSkQKNGEBhodTUiIiIiN03hLyeWLTN/DhigQ74iIiJSqCj83aq//oIffjDv9+tnbS0iIiIit0jh71Z98QWkpUFQENSta3U1IiIiIrdE4e9WXXvIV0RERKSQUfi7FSdPwubN5v2+fa2tRURERCQHFP5uxWefgWGYizpXq2Z1NSIiIiK3TOHvVly7sLOIiIhIIaTwd7MOH4Zdu8DDAx56yOpqRERERHJE4e9mZUz0uPde8POzthYRERGRHFL4u1kZ4U+HfEVERKQQU/i7Gfv2wYEDULQo9OpldTUiIiIiOabwdzOWLzd/hoZCmTKWliIiIiKSGwp/N2PlSvNnnz7W1iEiIiKSSwp/f+d//4OoKPD0hG7drK5GREREJFcU/v7OqlXmzw4ddMhXRERECj2Fv7+TcchXEz1ERETEBSj83UhsLGzfbt6//35raxERERHJAwp/N/LVV+bPVq2gShVraxERERHJAwp/N5Jxvl/PnlZWISIiIpJnFP6yk5gI69aZ93W+n4iIiLgIhb/srFsHqanQoAHUr291NSIiIiJ5QuEvO6tXmz816iciIiIuROEvO+vXmz8V/kRERMSFKPxl5+JFc4ZvUJDVlYiIiIjkGYW/G+nZE4qoi0RERMR1OH2ymTVrFrVq1cLHx4egoCA2b958w/03btxIUFAQPj4+3HbbbcyZMyfnb64lXkRERMTFOHX4W758OePHj2fq1KlERkbSrl07QkNDiY6OznL/o0ePct9999GuXTsiIyN57rnnGDt2LF9++eWtv3np0nD33bn8BCIiIiLOxWYYhmF1Edlp1aoVLVq0YPbs2fa2hg0b0rNnT2bMmJFp/2effZavv/6agwcP2ttGjRrFL7/8wrZt227qPU+ePEm1atU40asXVVesyP2HEBEREbkBe/Y4cYKqVavm+/s57chfSkoKERERhISEOLSHhISwdevWLJ+zbdu2TPt36dKFXbt2kZqamuVzkpOTSUxMtN+SkpLMDaGhuf8QIiIiIk7GacPfqVOnSEtLo1KlSg7tlSpVIi4uLsvnxMXFZbn/lStXOHXqVJbPmTFjBr6+vvZbo0aNzA3t2+f+Q4iIiIg4GacNfxlsNpvDY8MwMrX93f5ZtWeYMmUKCQkJ9ltUVJS5oUSJXFQtIiIi4pw8rS4gOxUqVMDDwyPTKF98fHym0b0M/v7+We7v6elJ+fLls3yOt7c33t7e9seJiYm5rFxERETEeTntyF/RokUJCgoiPDzcoT08PJy2bdtm+Zw2bdpk2n/t2rUEBwfj5eWVb7WKiIiIFBZOG/4AJk6cyIcffsiCBQs4ePAgEyZMIDo6mlGjRgHmIdvBgwfb9x81ahTHjx9n4sSJHDx4kAULFjB//nwmTZpk1UcQERERcSpOe9gXoF+/fpw+fZrp06cTGxtLYGAgYWFh1KhRA4DY2FiHNf9q1apFWFgYEyZM4L///S8BAQG89957PPjgg1Z9BBERERGn4tTr/FmhoNfaEREREfemdf5EREREJN8o/ImIiIi4EYU/ERERETei8CciIiLiRhT+RERERNyIwp+IiIiIG1H4ExEREXEjCn8iIiIibkThT0RERMSNKPyJiIiIuBGnvravFdLT0wHzusEiIiIi+S0jc2RkkPym8HedEydOANCyZUuLKxERERF3cuLECapXr57v72MzDMPI93cpRM6cOUP58uXZv38/vr6+VpdjuaSkJBo1akRUVBSlSpWyuhzLqT+uUl84Un84Un84Un84Un84SkhIIDAwkNOnT1OuXLl8fz+N/F3H09PskmrVqlG6dGmLq7FeYmIiAFWqVFF/oP64lvrCkfrDkfrDkfrDkfrDUUYfZGSQ/KYJHyIiIiJuROFPRERExI0o/F3H29ubF198EW9vb6tLcQrqD0fqj6vUF47UH47UH47UH47UH44Kuj804UNERETEjWjkT0RERMSNKPyJiIiIuBGFPxERERE3ovAnIiIi4kYU/q4xa9YsatWqhY+PD0FBQWzevNnqkvLFjBkzuOOOOyhVqhR+fn707NmT3377zWGfoUOHYrPZHG6tW7d22Cc5OZknn3ySChUqUKJECe6//35OnjxZkB8l16ZNm5bpc/r7+9u3G4bBtGnTCAgIoFixYtxzzz0cOHDA4TVcoR8y1KxZM1N/2Gw2Ro8eDbj+92LTpk306NGDgIAAbDYbq1atctieV9+Hs2fPMmjQIHx9ffH19WXQoEGcO3cunz/drbtRf6SmpvLss8/SpEkTSpQoQUBAAIMHD+aPP/5weI177rkn03emf//+Dvu4Qn9A3v37cJX+yOp3ic1m4/XXX7fv4yrfj5v5u+pMvz8U/v7f8uXLGT9+PFOnTiUyMpJ27doRGhpKdHS01aXluY0bNzJ69Gi2b99OeHg4V65cISQkhAsXLjjs17VrV2JjY+23sLAwh+3jx49n5cqVLFu2jC1btnD+/Hm6d+9OWlpaQX6cXGvcuLHD59y3b59922uvvcZbb73FzJkz2blzJ/7+/nTu3JmkpCT7Pq7SDwA7d+506Ivw8HAA+vTpY9/Hlb8XFy5coFmzZsycOTPL7Xn1fRg4cCB79uxh9erVrF69mj179jBo0KB8/3y36kb9cfHiRXbv3s3zzz/P7t27WbFiBf/73/+4//77M+07cuRIh+/M3LlzHba7Qn9kyIt/H67SH9f2Q2xsLAsWLMBms/Hggw867OcK34+b+bvqVL8/DDEMwzBatmxpjBo1yqGtQYMGxuTJky2qqODEx8cbgLFx40Z725AhQ4wHHngg2+ecO3fO8PLyMpYtW2Zvi4mJMYoUKWKsXr06P8vNUy+++KLRrFmzLLelp6cb/v7+xquvvmpvu3z5suHr62vMmTPHMAzX6YfsjBs3zqhdu7aRnp5uGIb7fC8MwzAAY+XKlfbHefV9iIqKMgBj+/bt9n22bdtmAMavv/6az58q567vj6zs2LHDAIzjx4/b2+6++25j3Lhx2T7HlfojL/59uFJ/XO+BBx4wOnbs6NDmqt+P6/+uOtvvD438ASkpKURERBASEuLQHhISwtatWy2qquAkJCQAZLqY9IYNG/Dz86NevXqMHDmS+Ph4+7aIiAhSU1Md+iwgIIDAwMBC12eHDh0iICCAWrVq0b9/f44cOQLA0aNHiYuLc/iM3t7e3H333fbP6Er9cL2UlBQ++eQThg0bhs1ms7e7y/fienn1fdi2bRu+vr60atXKvk/r1q3x9fUt9H2UkJCAzWajTJkyDu1LliyhQoUKNG7cmEmTJjmMdLhaf+T234er9UeGP//8k++++47hw4dn2uaK34/r/6462++PgrmCsJM7deoUaWlpVKpUyaG9UqVKxMXFWVRVwTAMg4kTJ3LXXXcRGBhobw8NDaVPnz7UqFGDo0eP8vzzz9OxY0ciIiLw9vYmLi6OokWLUrZsWYfXK2x91qpVKxYvXky9evX4888/efnll2nbti0HDhywf46svhfHjx8HcJl+yMqqVas4d+4cQ4cOtbe5y/ciK3n1fYiLi8PPzy/T6/v5+RXqPrp8+TKTJ09m4MCB9ovUAzz88MPUqlULf39/9u/fz5QpU/jll1/spxS4Un/kxb8PV+qPa3300UeUKlWK3r17O7S74vcjq7+rzvb7Q+HvGteOboD5H/D6NlczZswY9u7dy5YtWxza+/XrZ78fGBhIcHAwNWrU4Lvvvsv0j/daha3PQkND7febNGlCmzZtqF27Nh999JH9RO2cfC8KWz9kZf78+YSGhhIQEGBvc5fvxY3kxfchq/0Lcx+lpqbSv39/0tPTmTVrlsO2kSNH2u8HBgZSt25dgoOD2b17Ny1atABcpz/y6t+Hq/THtRYsWMDDDz+Mj4+PQ7srfj+y+7sKzvP7Q4d9gQoVKuDh4ZEpNcfHx2dK6a7kySef5Ouvv2b9+vVUrVr1hvtWrlyZGjVqcOjQIQD8/f1JSUnh7NmzDvsV9j4rUaIETZo04dChQ/ZZvzf6XrhqPxw/fpwffviBESNG3HA/d/leAHn2ffD39+fPP//M9Pp//fVXoeyj1NRU+vbty9GjRwkPD3cY9ctKixYt8PLycvjOuFJ/XCsn/z5csT82b97Mb7/99re/T6Dwfz+y+7vqbL8/FP6AokWLEhQUZB9mzhAeHk7btm0tqir/GIbBmDFjWLFiBevWraNWrVp/+5zTp09z4sQJKleuDEBQUBBeXl4OfRYbG8v+/fsLdZ8lJydz8OBBKleubD8Uce1nTElJYePGjfbP6Kr9sHDhQvz8/OjWrdsN93OX7wWQZ9+HNm3akJCQwI4dO+z7/PzzzyQkJBS6PsoIfocOHeKHH36gfPnyf/ucAwcOkJqaav/OuFJ/XC8n/z5csT/mz59PUFAQzZo1+9t9C+v34+/+rjrd74+bn7vi2pYtW2Z4eXkZ8+fPN6Kioozx48cbJUqUMI4dO2Z1aXnu8ccfN3x9fY0NGzYYsbGx9tvFixcNwzCMpKQk46mnnjK2bt1qHD161Fi/fr3Rpk0bo0qVKkZiYqL9dUaNGmVUrVrV+OGHH4zdu3cbHTt2NJo1a2ZcuXLFqo92y5566iljw4YNxpEjR4zt27cb3bt3N0qVKmX/7/7qq68avr6+xooVK4x9+/YZAwYMMCpXruxy/XCttLQ0o3r16sazzz7r0O4O34ukpCQjMjLSiIyMNADjrbfeMiIjI+2zV/Pq+9C1a1ejadOmxrZt24xt27YZTZo0Mbp3717gn/fv3Kg/UlNTjfvvv9+oWrWqsWfPHoffJcnJyYZhGMbhw4eNl156ydi5c6dx9OhR47vvvjMaNGhgNG/e3OX6Iy//fbhCf2RISEgwihcvbsyePTvT813p+/F3f1cNw7l+fyj8XeO///2vUaNGDaNo0aJGixYtHJY+cSVAlreFCxcahmEYFy9eNEJCQoyKFSsaXl5eRvXq1Y0hQ4YY0dHRDq9z6dIlY8yYMUa5cuWMYsWKGd27d8+0j7Pr16+fUblyZcPLy8sICAgwevfubRw4cMC+PT093XjxxRcNf39/w9vb22jfvr2xb98+h9dwhX641po1awzA+O233xza3eF7sX79+iz/bQwZMsQwjLz7Ppw+fdp4+OGHjVKlShmlSpUyHn74YePs2bMF9Clv3o364+jRo9n+Llm/fr1hGIYRHR1ttG/f3ihXrpxRtGhRo3bt2sbYsWON06dPO7yPK/RHXv77cIX+yDB37lyjWLFixrlz5zI935W+H3/3d9UwnOv3h+3/ixYRERERN6Bz/kRERETciMKfiIiIiBtR+BMRERFxIwp/IiIiIm5E4U9ERETEjSj8iYiIiLgRhT8RERERN6LwJyIiIuJGFP5ERJzA/PnzCQkJsez94+PjqVixIjExMZbVICIFQ+FPRCw3dOhQbDYbNpsNLy8vbrvtNiZNmsSFCxesLq1AJCcn88ILL/D888/b26ZNm4bNZqNr166Z9n/ttdew2Wzcc889mfa32Wx4eHhQrVo1RowYwV9//XVTNfj5+TFo0CBefPHFXH8eEXFuCn8i4hS6du1KbGwsR44c4eWXX2bWrFlMmjQpy31TU1MLuLq/l5uavvzyS0qWLEm7du0c2itXrsz69es5efKkQ/vChQupXr16ptdp3LgxsbGxREdHM3v2bL755hsGDx5803U8+uijLFmyhLNnz+bsg4hIoaDwJyJOwdvbG39/f6pVq8bAgQN5+OGHWbVqFWCOat1+++0sWLCA2267DW9vbwzDICEhgcceeww/Pz9Kly5Nx44d+eWXX+yv+csvv9ChQwdKlSpF6dKlCQoKYteuXQAcP36cHj16ULZsWUqUKEHjxo0JCwsDYNGiRZQpU8ahvlWrVmGz2eyPc1pTVpYtW8b999+fqd3Pz4+QkBA++ugje9vWrVs5deoU3bp1y7S/p6cn/v7+VKlShe7duzN27FjWrl3LpUuXWL16NXfddRdlypShfPnydO/end9//93h+U2aNMHf35+VK1fesF4RKdwU/kTEKRUrVsxhNO3w4cN89tlnfPnll+zZsweAbt26ERcXR1hYGBEREbRo0YJ7772XM2fOAPDwww9TtWpVdu7cSUREBJMnT8bLywuA0aNHk5yczKZNm9i3bx//+c9/KFmy5C3VmJOasrJ582aCg4Oz3DZs2DAWLVpkf7xgwQIefvhhihYt+rf1FStWjPT0dK5cucKFCxeYOHEiO3fu5Mcff6RIkSL06tWL9PR0h+e0bNmSzZs3//2HF5FCy9PqAkRErrdjxw4+/fRT7r33XntbSkoKH3/8MRUrVgRg3bp17Nu3j/j4eLy9vQF44403WLVqFV988QWPPfYY0dHRPP300zRo0ACAunXr2l8vOjqaBx98kCZNmgBw22233XKdOanpeufOnePcuXMEBARk+R7du3dn1KhRbNq0iaCgID777DO2bNnCggULbljbr7/+yuzZs2nZsiWlSpXiwQcfdNg+f/58/Pz8iIqKIjAw0N5epUoVIiMjb6kfRKRwUfgTEafw7bffUrJkSa5cuUJqaioPPPAA77//vn17jRo17CELICIigvPnz1O+fHmH17l06ZL9cObEiRMZMWIEH3/8MZ06daJPnz7Url0bgLFjx/L444+zdu1aOnXqxIMPPkjTpk1vqeac1HS9S5cuAeDj45Pldi8vLx555BEWLlzIkSNHqFevXrZ17tu3j5IlS5KWlkZycjL33HMP8+bNA+D333/n+eefZ/v27Zw6dco+4hcdHe0Q/ooVK8bFixdvsgdEpDBS+BMRp9ChQwdmz56Nl5cXAQEB9sOzGUqUKOHwOD09ncqVK7Nhw4ZMr5Vxvt60adMYOHAg3333Hd9//z0vvvgiy5Yto1evXowYMYIuXbrw3XffsXbtWmbMmMGbb77Jk08+SZEiRTAMw+E1s5rQkZOarle+fHlsNtsNJ1kMGzaMVq1asX//foYNG5btfvXr1+frr7/Gw8ODgIAA++gjQI8ePahWrRoffPABAQEBpKenExgYSEpKisNrnDlzxiHQiojrUfgTEadQokQJ6tSpc9P7t2jRgri4ODw9PalZs2a2+9WrV4969eoxYcIEBgwYwMKFC+nVqxcA1apVY9SoUYwaNYopU6bwwQcf8OSTT1KxYkWSkpK4cOGCPeBlnNOXFzVdq2jRojRq1IioqKhs1/lr3LgxjRs3Zu/evQwcOPCGr5VVH54+fZqDBw8yd+5c+4ziLVu2ZPka+/fvd1hCRkRcjyZ8iEih1KlTJ9q0aUPPnj1Zs2YNx44dY+vWrfzzn/9k165dXLp0iTFjxrBhwwaOHz/OTz/9xM6dO2nYsCEA48ePZ82aNRw9epTdu3ezbt06+7ZWrVpRvHhxnnvuOQ4fPsynn37qMOkipzVlp0uXLtmGsQzr1q0jNjY22xHEGylbtizly5dn3rx5HD58mHXr1jFx4sRM+128eJGIiAhLF5sWkfyn8CcihZLNZiMsLIz27dszbNgw6tWrR//+/Tl27BiVKlXCw8OD06dPM3jwYOrVq0ffvn0JDQ3lpZdeAiAtLY3Ro0fTsGFDunbtSv369Zk1axYA5cqV45NPPiEsLIwmTZqwdOlSpk2bluuasjNy5EjCwsJISEjIdp8SJUrkKPgBFClShGXLlhEREUFgYCATJkzg9ddfz7TfV199RfXq1TOtNygirsVmXH9ii4iIFLi+ffvSvHlzpkyZYlkNLVu2ZPz48Tc8tCwihZ9G/kREnMDrr79+y+sM5qX4+HgeeughBgwYYFkNIlIwNPInIiIi4kY08iciIiLiRhT+RERERNyIwp+IiIiIG1H4ExEREXEjCn8iIiIibkThT0RERMSNKPyJiIiIuBGFPxERERE3ovAnIiIi4kb+D2LrSBAeGiGBAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Lets do it for just one temperature (e.g. perhaps for a conference talk)\n", "fig, (ax1) = plt.subplots(1, 1, figsize=(7,5))\n", "ax1.plot(100*P_combo*10**(-8), CalcD_P_1200/1000, '-', color='red')\n", "ax5 = ax1.twiny()\n", "ax1.set_xlabel('Pressure (MPa)')\n", "ax1.set_ylabel('Density (g/cm$^{3}$)')\n", "ax5.set_xlim([D_axis_min, D_axis_max])\n", "ax1.set_xlim([P_axis_min*100, P_axis_max*100])\n", "ax5.set_xlabel('Depth (km)')\n", "fig.savefig('Pressure_Depth.png')" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.5" } }, "nbformat": 4, "nbformat_minor": 5 }