{ "cells": [ { "cell_type": "markdown", "id": "4c4ed7f1-9e71-4d8c-bc51-9b972bf5a8cc", "metadata": {}, "source": [ "# Phase diagram of a pure substance" ] }, { "cell_type": "code", "execution_count": 1, "id": "b6b2b5bb-4c0c-49c9-bea2-1d34055a57dd", "metadata": {}, "outputs": [], "source": [ "from feos.si import *\n", "from feos.eos import EquationOfState, PhaseDiagram\n", "from feos.pcsaft import PcSaftParameters\n", "import matplotlib.pyplot as plt\n", "import seaborn as sns\n", "import pandas as pd\n", "\n", "sns.set_context('talk')\n", "sns.set_palette('Dark2')\n", "sns.set_style('ticks')\n", "colors = sns.palettes.color_palette('Dark2', 1)" ] }, { "cell_type": "code", "execution_count": 2, "id": "a19bbe80-cc85-4911-a169-d9501c962983", "metadata": {}, "outputs": [], "source": [ "parameters = PcSaftParameters.from_json(\n", " substances=['methanol'], \n", " pure_path='../parameters/pcsaft/gross2002.json'\n", ")\n", "eos = EquationOfState.pcsaft(parameters)\n", "phase_diagram = PhaseDiagram.pure(eos, 150.0 * KELVIN, 201)" ] }, { "cell_type": "code", "execution_count": 3, "id": "588d342b-720f-4f59-9046-e468c31608fa", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA6UAAAGvCAYAAACjJXMJAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAC1c0lEQVR4nOzdd1gU5xbH8e8uVZqABbsiCio2FHtBxd5jj11j7CYaTaImMYkxscRYEmOLvRt7F7vGLvYC2BAQFRUs9Dr3D2SvKzYUGMr5PA/PldnZ3d8uuTt7Zt73vBpFURSEEEIIIYQQQggVaNUOIIQQQgghhBAi+5KiVAghhBBCCCGEaqQoFUIIIYQQQgihGilKhRBCCCGEEEKoRopSIYQQQgghhBCqkaJUCCGEEEIIIYRqpCgVQgghhBBCCKEaKUqFEEIIIYQQQqhGilKR6YSFhRESEqLKc8fExBAUFKT7fePGjTg5OXHq1ClV8ojU8+zZMwYNGkTFihWpUqUK165dUzVPQECA7t93797FycmJv/76S8VEb5bR8wkhMo8GDRrQo0cP3e89evSgQYMGKibKmBRF4ffff6datWpUrFiRlStXqponODiYiIgI3e+v/h0zmoyeLzsyVDuAEClx5coVBg0axNSpU6lWrVq6PndgYCB9+/ZlwIABtGvXDoAqVaowZcoUHBwc0jWLSH1z587lwIED9O7dm+LFi1OsWDHVsowbNw5fX1+WL1+uWgYhRNYTHh6Oq6srCQkJ77X/6dOnyZkzZxqneruBAwcSGRmpaoaM6NChQyxYsIB69erRsGFDKleurFqWw4cPM2rUKDZt2oSZmZlqOUTmJkWpyFSuX7/Ow4cPVXnuu3fvcufOHb1thQsXpnDhwqrkEanLx8cHa2trxowZo3YUjh49SsGCBdWOIYTIYuLj45k0aZLettWrV3P+/Hm+/fZbcuXKpdtubGysekEKUKtWLbUjZEg+Pj4AfPXVVzg5Oama5dKlSzx//lzVDCLzk6JUCCGA2NhYzM3N1Y4hhBBpxsrKijZt2uhtW7p0KSYmJvTs2RNDQ/lamFnExsYCyHFLZBkyp1SkuWfPnjF69Gjq1atH2bJladiwIX/88QfR0dF6++3evZvu3btTuXJlypYtS4MGDZgyZQoxMTEA/PXXX7qrWD179tTNMRk9evRrzxK+un306NE0bdqUlStXUqVKFapUqcKRI0cAOHHiBP369aNatWo4OztTp04dxo0bpzvzt3HjRnr27AnAmDFjdI/7ujmlkZGR/PHHHzRo0ED3OqZOnao3/Cjpft7e3owcOZIqVarg4uLCkCFDuHv37lvfz6T7Xrp0STcHsnbt2kycOJGoqKhk+3l4eNCgQQMqVKigm/MXHR3N9OnTdRnd3d2ZOXOm7r1O4uHhQfv27XFxcaFy5cr06dOHs2fP6u1z+vRpunXrhqurKy4uLnTp0oUDBw7obn/TfMNXtyf9vmTJEj799FPKli1L79699V5P27ZtKVeuHNWrV2f06NHvfdV83bp1tGnTRnffkSNH6t7npOc9ffo0gYGBODk5MXr06Nc+TtK+W7duZfLkydSsWRMXFxcGDx5MSEgIly5dokuXLlSoUIEmTZqwc+dOvfsnJCSwaNEimjZtStmyZalTpw4TJkwgLCxMt4+TkxOBgYGcPn0aJycnNm7cqLstLi6O6dOnU7duXcqXL0+nTp04ffq03nOEhYXxxx9/0LRpU8qVK4eLiwudOnVi//79yV7H5s2bdY9Xrlw5OnbsyMmTJ1P0/gkhMq/Y2Fhu3LiBk5PTBxWkISEhjB49mmrVqlG6dGlq167NrFmzUBTlve6/c+dO2rRpQ/ny5WnZsuVrP39enVOqKAqrV6+mQ4cOuLi4UK5cOZo2bcr8+fOTPe/hw4fp2LEjFStWxN3dnRUrVvDdd9/pPV6DBg34/vvvGTt2LOXKlaNu3bqEhIS89/M0aNCA8ePHs27dOpo0aUL58uVp3749ly5d4tGjR3z55Ze4uLhQp04dpk+f/l5Dpz09PenduzcuLi64uLjQs2dPzpw5o/ecs2bNAsDd3f2tc24/Nt/58+fp06ePLkvfvn25dOmS7vbRo0frZXl1nua2bdto0aIFZcuWpUmTJqxevVrv9pS8z+PGjWPLli20aNGCcuXK0bhx49fOpX3X+ycyJjklJtLc8OHDuXbtGj179iRv3rycP3+e+fPn8/TpU3755Rcg8Uvv999/T4MGDRg1ahSxsbHs3buXhQsXYmZmxtChQ2nUqBGPHj1i7dq1DBw4kHLlyqU4y/379/n7778ZOnQoDx8+pEKFChw9epTPP/+cSpUq8cUXX6DRaDh27Bhr164lNjaWiRMnUqVKFQYOHMjcuXPp3LnzG+duxMTE0KdPHy5cuEC7du0oW7Ysly5d4p9//uHs2bMsW7YMIyMj3f6DBg3CwcGBESNGEBAQwNKlSwkKCmL9+vXvfC1ffvklefPmZeTIkXh5ebFkyRJu3rzJwoUL9fYbM2YMPXr0wNLSkooVKxIfH8+AAQM4d+4cnTp1wsHBgStXrjB37ly8vLyYM2cOGo2G06dPM2LECOrWrUvHjh2JjIxkxYoV9OnThx07dlC4cGFu377NgAEDKF26NCNGjADg33//ZfDgwaxYsQJXV9cU/41mzpyJm5sbrVq1wsTEBIBZs2bx119/0aRJEzp16kRQUBArVqzg9OnTrF+/Hltb2zc+3uTJk1m0aBE1atTgm2++4eHDh6xYsYLjx4+zbt06bG1tmTJlCnPnzuXJkyeMGTOGIkWKvDXj1KlTyZMnD0OHDuXmzZusXLmSJ0+ecPv2bdq1a0erVq1YtmwZ33zzDc7OzhQtWhSA7777js2bN/PJJ5/Qu3dvbt26xerVqzl37hyrV6/GxMSEKVOmMHHiRGxsbBg4cCCVKlXSPe/SpUspUKAAn332GVFRUSxcuJD+/fvj4eGBnZ0diqIwYMAArl27Rvfu3SlSpAgPHjxgzZo1DBs2DA8PD73h5jNnziRHjhz07duX2NhYFi1axIABAzh06BA2Njbv9f4VKlQoxX9jIUTGcPPmTWJiYihdunSK73v9+nX69OmDpaUlvXv3xsTEhH///Ze//voLa2trunfv/tb7b9y4kTFjxuDi4sLXX3+Nn58fAwcOJCEh4a3TF2bMmMHcuXP55JNP6NSpE+Hh4WzevJk//viDPHny8MknnwBw8OBBhgwZgqOjIyNGjCAoKIjJkydjZmaW7Orijh07sLe357vvvuPx48fY2toyffr093oegP3797Nnzx569eqFoijMmTOHYcOGYWlpScmSJRk9ejR79uxh7ty5FCtWTO++r9q/fz9Dhw6lSJEiDBo0CEj8jtS7d2/+/PNP3N3dGTt2LJs3b2bv3r2MGTPmnZ/DH5rv2LFjDBgwgFKlSvHll18SExPDxo0b6datG4sXL8bV1ZXOnTsTFhamy1KyZEnd816+fJnr16/TvXt3bG1tWbNmDT/99BN58uShYcOGKfp7Avz333+6Cxi5c+dm7dq1jB8/nkKFCuHm5vbe75/IoBQh0tDjx48VR0dHZcGCBXrbR48erfTq1Uv3e9OmTZXOnTsrCQkJum2xsbFK3bp1lZYtW+q2bdiwQXF0dFROnjyp2/btt98qjo6OyZ771e1Jv2/YsEFvv88++0ypX7++Eh0drbe9U6dOiouLi+73kydPJrv/q3lWrVqlODo6KosXL9Z7rH/++UdxdHRUVq5cqXe/oUOH6u03btw4xdHRUfH19U32el59zvbt2+tlnjZtmuLo6KgcOXJEb79vv/32tfdP2i/JmjVrFEdHR2Xv3r2KoijKjz/+qLi4uOj9Tby9vZXGjRsru3btUhRFUebPn684OjoqwcHBun1CQkKUxo0bK8uWLVMURVECAgIUR0dH5c8//9R7vle3J/3eqFEjJSYmRrefv7+/UqpUKWXq1Kl69/fx8VGcnZ2VX3/99Y3v1c2bNxUnJydlyJAheq/jwoULipOTk/Lll1/qtnXv3l2pX7/+Gx/r5Yx169ZVIiMjddvbt2+vODo6KqtWrdJtO3bsmOLo6KisXbtWUZT///ezevVqvcf877//FEdHR2XJkiW6bfXr11e6d++e7Hnr1KmjhIaG6rZv2rRJcXR0VNatW6d7Xa97jiNHjiiOjo7KokWL9B7Pzc1NCQ8P1+23Y8cOvczv+/696W8shMjY1q9frzg6Oipr1qxJ0f2ioqIUd3d3pVu3bkpERIRue1BQkFKmTBmlTZs2b71/XFycUqNGDaV9+/Z6n/dJx6eXP/9e/myOiYlRKlWqpIwYMULv8UJDQ5WyZcsqAwYM0G1r2LCh0rhxY73P6r179yqOjo56n/X169dXSpUqpfj5+em2peR56tevrzg5OSne3t66bZMnT1YcHR2V4cOH67aFh4crzs7OyldfffXG9yXpe4+bm5veZ/2zZ8+UOnXqKHXq1NG9X3/++afi6OioBAQEvPHxPiZffHy84u7urnTp0kWJi4vT269Ro0Z6f+PXZUl63itXrui23b17V3FyclK+/vprRVE+7H328vLSbXv48KHi5OSky5yS9+/V46xQnwzfFWnK0tISMzMzVq1ahYeHh65d+MSJE1myZIluv61btzJ//nw0Go1uW3BwMFZWVnotxlND7dq19X6fN28eGzZswNjYWLftyZMnWFhYpPi5Dxw4gIWFBd26ddPb3rNnTywsLPSGUAI0a9ZM7/eks9WPHz9+53P17dtXL3OfPn10GV726uvds2cPtra2ODs7ExISovtxc3PDwMCAQ4cOAZAvXz7Cw8OZMGECt27dAtANB27atKluH4BffvmFK1euAGBjY4OHh8cHt1qvXr263tXkvXv3kpCQQIMGDfTy5s6dm9KlS+vyvs6BAwdQFIX+/fvr/bdVoUIFatWqxaFDh4iLi0txxjp16mBqaqr7PalTb6NGjXTbks5cJw0x3rNnDxqNBjc3N73XUaZMGfLkyfPW15GkXr16WFhY6H5PGi2Q9N9LhQoVOHPmjK47NCQ2NkkajhUeHq73eG5ubnqdEkuVKgXAo0ePgLR7/4QQGUPS0lcpvVK6fPlyAgMD+e2338iRI4due968eSlRogQPHjx46/2vXr1KcHAw7dq10/u8b9OmzVubKxkZGXH8+HHGjx+vt/3VY7a3tzf+/v506dJF77O6YcOGr+2WX6RIEb0RMu/7PC/f/+XpQvb29oD+McHMzIxcuXLpPl9f59q1azx48IBu3brpfdZbWVnRvXt3goKCdMfalPiQfNeuXSMgIICGDRvy7Nkz3TErKiqK+vXr4+Xl9c6/c7FixXB2dtb9XrBgQWxtbXXHrJS+z/b29rrjFECePHnInTu37vHS6v0T6UOG74o0ZWxszPjx4/nhhx/44osvMDY2pmrVqjRu3Ji2bdvqhmcaGRlx5swZtm/fzu3bt/H39yc4OBgg1buQvtxdEMDAwICAgABmzpzJzZs38ff311uLNCXu3r1L4cKF9Q6ykPg+FC5cmMDAQL3tSUMkX94PEguJd3n1wGptbY21tXWy53j19fr7+xMSEkKNGjVe+7j3798HoHv37hw9epQVK1awYsUKChUqRP369enQoYPuoNC0aVP27t3Lzp072blzJ3ny5MHNzY1PPvnkg4buAsmG4vr7+wPQpUuX1+7/6nv9sqR5j0kH4Jc5ODhw9OhRnjx5Qp48eVKU8dX3NGku1svbtdrEc37Kizkx/v7+KIpCvXr1XvuY79Os4tXnTfr/T1LDi6Qsa9as4fTp0/j5+eHv76+ba6y8Mt/q1fc66b+/pCL2fd8/IUTmdO3aNQwNDZP1ZYiKiiI0NFRv28vDXnfv3k3FihVfO9UhKipKV1i+6XGSjlOv3t/AwEA33eFNjIyMOHToEPv378fX1xc/Pz+ePXsG/P8zzs/PD+C1j2Vvb4+Xl5fetlc/W9/3ed50fwMDAyD5Z6yBgcFb59u+7TO3ePHiANy7dw8XF5c3PsbrfEi+pGPvlClTmDJlymsf9/79+7qT0+/zvACmpqZ6x6yUvM+vm6pjbGz8Xsesj3n/RPqQolSkuVatWlGnTh327dvH4cOHOX78OEePHmXVqlWsW7cOY2Nj/vjjD+bPn0+ZMmWoWLEibdq0wcXFhV9++UVXJKXUmwq7pA/jJGvWrOHHH3/E3t4eV1dXGjduTIUKFVi+fDnbtm1L0XO+7WCTkJCQrIBKKlw+xOuKsfj4+GSP+erv8fHxFCtWjB9//PG1j2tlZQWAhYUFK1as4MKFC+zbt48jR46wfPlyVq5cyZQpU2jVqhVGRkb8+eef+Pj4sHfvXo4cOcLGjRtZv349I0eOpH///m/M/6ZmD6/+fZL2mzNnjt4Z7/fxrr8HvL2ofZM3NQR5+Wri657P3Nxc1xDiVUkF5tu867+X58+f06VLFwICAqhVqxYNGjSgVKlSFCxYkI4dO6b48d73/Xu1aZkQIuNLSEjA29ub4sWLJ/v82blzZ7LlsYYOHcqwYcN49uwZV65cee1omKdPn3L37l3dPMA3PU7SvMPXfXa8rRGQoih8/fXXbN++ncqVK+Pi4kLnzp2pUqUKvXr10u2XNILj5dFESV73Wfvqced9nyfJhxwT3vT63nVbeh6zILF/RcWKFV+7T1Kh9ybvc4xJyfv8Mcesj3n/RPqQolSkqfDwcLy8vChZsiQdOnSgQ4cOxMTE8Pvvv7Ns2TKOHj2Kk5MT8+fPp02bNsnOxr3PMNakD6mYmBi9A9D73Dc6OppJkyZRrVo1Fi1apPfBPXPmzPd9mToFCxbkwoULxMbG6n3wxcTEcPfu3Q++evg6AQEBegeEkJAQQkNDdUNJ36RQoUJcuXKF6tWr633AJzWXSjrr6evrS2hoKBUrVqRixYqMGjWKmzdv6hoctGrVinv37nHv3j1cXV1xcnJi6NChPHjwgF69euma8CQd7F/t7Pu2IUwvS7pSnj9//mRDzA4fPqw3ROd1rxXg9u3bVKhQQe82X19fzMzM0m0dvoIFC3L06FHKli2rK/yTeHh4YG1t/dHPsWzZMm7dusWSJUv0roSfO3fugx7vfd+/l7sHCyEyhzt37hAREaE3vDJJ7dq1Wbx4sd62pCZp165dQ1EUvaH/SbZv305cXBzNmzd/6+MkdbZ/de1vRVEIDAzUa5bzMk9PT7Zv387gwYP58ssvddvj4uJ4+vSpLmPS/965cyfZFJZXn/Njnie1JR3vbt++new2X19fgLdemUyLLGZmZtSsWVPvtkuXLvHs2bMUnyh+VWq/zxnp/RMpJ3NKRZq6ceMG3bp10+sma2xsTJkyZYDEs5NJwzRKlCihd9/Dhw9z584dvTlrSUXUy2dSk4Zevjwc58GDB5w/f/6d+aKiooiMjKRYsWJ6BamXl5duqY2k508qrt52FrdBgwaEhYUla1G+atUqwsPD3zh080OsWLFC76xgUtfdl+eIvCnj06dPk7VlX7NmDSNGjODEiRMATJgwgcGDB+vNQyxevDhWVla6v8PcuXPp3bu33nDnfPnyYWdnp9vH2toaQ0PDZMOldu3a9V6vs379+kDi3N+XX6+XlxeDBg1i6dKl77zvP//8o3ffq1evcvz4cdzc3FJ8JvtDJbXsnzNnjt72AwcO8MUXX+hdlddqte+1bMCrnj59Cuj/f0lRFFasWAGQ4vmfGen9E0KkrrfNJ82bNy81a9bU+0kqEJI+y0+ePKn3uXD79m1mzJhBjRo1qF69+lsfp0yZMhQsWJDVq1frLZe2Y8eOt04JeN1nHCR2fY+MjNR9xpUtW5b8+fOzfv16vROiFy5c0L3ut3nf50ltzs7O5MmTh9WrV+ud7AsLC2PVqlXkyZOHsmXLpslzv6ps2bLkyZOH5cuX630PCAsLY/jw4YwZM0b3vejV6SrvK7Xf54z0/omUkyulIk1VqFABV1dXpk+fzv3793FycuL+/fusWLGC4sWL667mFChQgLlz5xIdHU2+fPm4dOkSmzZtwsTERO/DMGk+werVq3n8+DGtWrWiWbNmzJs3jxEjRtC7d2+io6NZuXIldnZ27zwjmjNnTipUqMDGjRuxsLDA3t6eGzdusG7dOt2HbHh4ODlz5tTN/9y6dSuKory2pXvHjh3ZtGkTkyZN4vr165QtW5YrV66wceNGKlSo8NohlB/q1KlTfP7559SvX5+LFy+yZcsW2rZt+8blal7N+Msvv3D16lXKly/P9evXWbt2Lc7OzromOX369OHzzz+nW7duuvm/+/btw9/fn8mTJwPQrVs3tmzZQrdu3ejcuTM5c+bk5MmTnDp1ii+++AKAHDly4O7ujoeHB2PHjsXFxYWTJ09y/vz59xpG4+joSI8ePVi+fDlPnz6lYcOGPH36lBUrVmBubq53dvVVJUuW1N23T58+NGzYkEePHrF8+XKsrKwYOXLk+77dH83NzQ13d3cWLVrE3bt3qVmzJoGBgaxcuVK3zEsSW1tbvL29WbVqFVWrVn3vs9F169Zl+fLlDBgwgA4dOhAbG8uuXbu4cuUKWq02WaOjd8lI758QInUlFWevu1L6NlevXiVXrlwEBgYyaNAg6tSpw927d1mzZg3W1tZMmTLlncMsNRoNP/zwA0OGDKFz5860b9+eoKAgVq5c+dZRIy4uLlhYWDBx4kTu3buHlZUVp06dYufOnXrfF7RaLaNHj2b48OF06dKFNm3aEBISwrJly147pPdDnye1GRkZ8cMPPzB8+HDat29Phw4dAFi/fj0PHz7kzz///KhpPx+apV27dnTo0AETExPWrVvHvXv3mDp1qu5kftJ3swULFlC3bt33XnYltd/njPT+iZSTv4xIUxqNhr///psuXbpw8OBBxo8fz7///kvjxo11BwdjY2Pmz5+Pi4sLy5YtY/LkyVy9epWxY8cyatQowsLCdN3SatSoQbNmzTh8+DC//PIL0dHRlCpVihkzZmBubs6UKVP4999/+fzzz+nUqdN7ZZw5cyYNGjRgw4YN/Pbbbxw/fpz+/fszdepUAN1i3g4ODvTo0YMrV67w22+/ce/evWSPZWxszJIlS+jTpw/Hjx/nt99+4/Tp0wwYMCDZGqUf67fffkNRFCZPnszZs2cZOXIkEydOfOf9Xs548uRJJkyYwKFDh/j0009ZuHChrpNi7dq1mTNnDjly5GD27NlMmjSJZ8+eMW3aNNq2bQskduNdvHgxRYsWZdGiRfzyyy/cuHGDH374gcGDB+ue8+eff+aTTz5h3759/Pbbb0RGRrJ8+fL3fj++++47fvzxR0JCQpg8eTKrVq3C1dWVVatWvbaT4qv3HTduHMHBwUyaNIn169fTqFEjNm7cmGZDsF5Ho9Ewc+ZMhg8fzvXr1/n111/ZunWrbvHv3Llz6/YdNmwYOXPm5LfffmPv3r3v/Rx169ZlwoQJREZGMmnSJBYsWIC1tTVr166ldOnSnDp1KsW5M8r7J4RIXV5eXmg0Gr1upu97v9KlSzN79mzdZ/LWrVtp1aoV69evJ2/evO/1OPXr12fevHmYmpoybdo09u3bx6+//vrWeYq5c+dm/vz5FC5cmNmzZzNt2jTu3bvHtGnT6Nq1Kzdv3tRN3WnatCnTp08nLi6O33//ne3btzNmzBjKli37zsI0Jc+T2po0acKiRYvImzcvf//9N/PmzaNQoUIsW7ZMt7ZneknKYmdnx+zZs5k5cybm5ubMmTOHli1b6vZr0aIFNWvWZOPGjbrvTu8jLd7njPT+iZTRKCm91i6EUFXSguPLli2jWrVqascRQgiRTURGRlK5cmV69+7NN998o3acN4qPj+fZs2ev7dbaqlUrrKyskk2zEUKoS4bvCiGEEJlUXFwclSpVStbF1MzMTDevvlGjRrrlHV524sQJ3Zf2y5cvM2XKFK5cuYK5uTnt2rVj2LBh0qlS6PH29iY+Ph5HR0e1o7xVfHw8devWpV27dnprYF6/fp0bN27o1vUWQmQcUpQKIYQQmZSvry/R0dFMnjxZr/P2y3PiAwICGDlyJFWrVtW7b1IXaD8/P3r37o2LiwszZszg1q1bTJ8+nbCwMMaNG5dur0VkfEnzUFM65De9GRsb07RpU9avX49Go6Fs2bI8fPiQ1atXY2NjI0WpEBmQFKVCCCFEJuXt7Y1Wq6VJkya6+eAv8/HxQVEU3N3d3zj/ev78+VhaWjJ79myMjY1xc3PD1NSUCRMmMGDAAOzs7NL6ZYhMwsvLC0NDw3euT5kRTJgwAXt7e7Zu3cqmTZuwtLSkRo0aDB8+/L3nvQoh0o/MKf1ASUslHDx4UOUkQgghsqspU6awf/9+PDw8Xnv7ypUrmTx5MufPn9ct3/CqevXqUb9+fX788UfdtpCQEGrUqMFvv/1G+/btU5RJjo9CCCFSSq6UfqDQ0FC1IwghhMjmfHx8MDY25rPPPuPcuXMYGhrSrFkzvvnmGywsLPDx8cHa2pqvvvqKY8eOER8fT7169Rg7dix58uQhMjKS+/fvY29vr/e4tra2WFhY6Bacf5mrq+tbM4WGhmJpaZmqr1MIIUTWJkvCCCGEEJmUt7c3/v7+uLm5MX/+fAYPHsz27dsZNGgQiqLg7e3N48ePKVmyJHPnzmXMmDGcOXOGnj17EhUVpTvBamFhkeyxzc3N9RagF0IIIdKKXCkVQgghMqnp06eTM2dOnJycAKhSpQq5cuXi66+/5vjx43z//fcoikKFChWAxKucDg4OdO3ala1bt+Lm5gYkrqP7KkVRXrvQvKen51szvetKqhBCCPEquVIqhBBCZFJVq1bVFaRJ6tWrByReRS1fvryuIE1SuXJlLC0t8fb21l0hfd0V0YiICBmGK4QQIl1IUSqEEEJkQsHBwaxbt46AgAC97VFRUQCYmpqyYcMGvL299W5XFIXY2FhsbGwwNzfHzs4OPz+/ZI8dFhaWbK6pEEIIkRakKBVCCCEyIY1Gw7hx41ixYoXe9p07d2JgYECNGjWYPHkys2bN0rt9//79REVF6dYtrVWrFgcPHiQmJka3j4eHBwYGBsnWNhVCCCHSgswpFUIIITIhW1tbunXrxvLly7GwsMDV1ZWzZ88yd+5cunXrRvHixRk0aBCTJk1iwoQJNGjQgOvXr/PXX3/h7u5OtWrVAOjXrx87duygf//+9OrVizt37jBt2jQ6depEgQIFVH6VQgghsgNZp/QDJTVyeFfDByGEECKtxMbGsmTJEjZs2EBgYCB2dnZ06tSJfv366ZoUrVu3jmXLluHv70/OnDlp1aoVw4YNw9TUVPc4np6eTJkyBS8vL2xsbGjbti3Dhg3DyMgoxZnk+CiEECKlpCj9QHLQFUIIIZKT46MQQoiUkjmlQgghhBBCCCFUI0WpEEIIIYQQQgjVSFEqhBBCCCGEEEI1UpQKIYQQQgghhFCNFKVCCCGEEEIIIVQj65RmAuNObuVg4HV+rtaKBoWc1I4jhBBCiCwuPiGBpzERPI4MJzgqjOCocIKjwnkcFUZwZNLvif97N+wp0fGx5DK1IE8OC3KZWmBraoaNiTm2pmbYmpiTy9QcW1Pzl7abY2IgX0OFEInk0yATuBQciO/zx/Tau4QvKtTnq4oNMdDKRW4hhBBCfJgn0RFceBTAleB7PIh4TkhSwfmi0HwSHUFCClcNfBwVxuOosPfe38LIBFsTc2xeFK62pmbkMjWngLk1pW3yUcY2Pzam5il9aUKITEiK0kxgdr2uDDq4krOP/Jl58QDnHvkzy60LuUwt1I4mhBBCiAwuOj6OayH3Of/In/OPAjj/KIA7ocEpegwbk8SCMZepBblMzcmdwwJbU3Nym1pgZWxKeGw0OQyNeRIdTkhUBMFR4YREhet+D4lOXuiGxUYTFhuNf1jIG583v1lOytjm1/spZplLTs4LkcVoFCWFp8EEkP6Lg8fEx/Gr5y4WXjsGQD4zK+bW64arXdF0eX4hhBDifaT38VG8nl9oMAcCfDhw14fjD24RHR+XbB8DjRYnGzuKWNgmKzSTCtDcOcyxMTHDUGvw0ZniExJ4HhNJSHQEIS+K1uAXReuTqHBCohOHCN95Hozv82AUXv8V1cLIhMp5ilDVrhhV89lTMXdhchgafXQ+IYR6pCj9QGoddLf7XmLk0fWEx8VgqNHyfZXmfFamFhqNJl1zCCGEEK8jRak6FEXh3CN/tvle4sBdH24/f5xsn4Lm1rjkKfzipwjlchUgh6GxCmnfLSI2Bq8nD/AKuc+1J/e5FpL4ExEXk2xfI60B5XMVpFo+e+oUKEFVO3uZrypEJiNF6QdS86B769kj+h9Ygc/TIABaFCvH1FrtsTQ2TfcsQgghxMukKE1fPk+C2Hz7AltuX0w2DDaXqTn1CjrSoFApquezx87MSqWUqSNBScAvNIRzjwI4E3SH00G+XH/6MNl+pgZG1MxfHLeCjtQr6Ehxq9xy8l6IDE6K0g+k9kE3IjaGMSc2seHWeQDsrXIzv353StvmUyWPEEIIAeofH7ODyLgYtvpeYpn3SS4+vqt3W2mbfDQrWpYGhZwon7sgWk3WnnsZEhWO50M/TgXd4cT9W1wKDky2TzHLXDQp6kzTImWolKeIzEcVIgOSovQDZYSDrqIorLx+mnEntxKTEI+pgRETa7SlY8nKqmUSQgiRvWWE42NWFRAawsJrx1h38yzPYqJ024ta2tLGvgJtilfEycZOxYTqC44K40jgTQ4HXufIvRs8jAzVuz1PDgsaFS5DG/vyVM9XXApUITIIKUo/UEY66F56fJcBB1cSEPYEgK6OVRlfrRWmMulfCCFEOstIx8es4kpwIHOuHGG772XilQQAjLUGNC9Wjh5O1ahqV0yGp76GoihcDbmHh/81dvtdxevJA73b85lZ0bZ4Rdo5uFDGNr9KKYUQIEXpB8toB90n0RGM+O9f9gV4A1DWtgDzG3SniKWtysmEEEJkJxnt+JiZXQ2+x9Tze9kb4KXblt8sJ71L16CLo6ssDZdCfqHB7PG/xjbfy5x75K93WymbfHxSvCKfFK9IAQtrdQIKkY1JUfqBMuJBN0FJYPblw0w5t4cERSGncQ5muXWhfiEntaMJIYTIJjLi8TGz8X3+mMlnPdh+57Jum6N1XgaVdaNN8QoYS2fZj+b7/DGbb19g460L+L7UqViDhur57OlcsjIti5WXUWdCpBMpSj9QRj7oHrt/i8GHVhEcFY4GDaNcGjKsQv0s3+xACCGE+jLy8TGjexodwcyLB1jidYLYhHggsRgd6dKIZkWd5TieBhRF4cLju2y8dZ6tvhcJjgrX3WZjYkankq70cKpGMatcKqYUIuuTovQDZfSD7r3wZ/Q/sIILjwMAaFS4NDPrdsZKlo0RQgiRhjL68TEjUhSFzbcv8tPpbbqiqJCFNd9UakIb+wrSjCedxCbE89+9m6y7cZZdfleIezF/F8CtoCO9SlXHvVAp+XsIkQakKP1AmeGgGx0fx4+ntrHC5xSQuGzMPw26U8pGlo0RQgiRNjLD8TEjuRv2hDHHN3Mw0AcAc0NjhlVoQL8ytWToqIoeRoSy5sYZVvic4l74M932AuY56VmqOt2dqmFtYqZiQiGyFilKP1BmOuiuvn6G709uITo+DjNDY/6o3YFW9uXVjiWEECILykzHRzXFJySw2Os4U87tISIuBoAmRcowoXob8pvnVDmdSBKXEM+Buz4s9T7J4cDruu3mhsZ0dapKvzK1KSiNkYT4aFKUfqDMdtC9+Pgu/Q+sIDD8KQD9nWsz1rUZhloDdYMJIYTIUjLb8VEN98KfMezwak4F3QEgbw5LJlRvQ/NiZdUNJt7qzvNglnmfYNX1M4TFRgNgqNHSpngFBpStK8vKCPERpCj9QJnxoBsSFc7gQ6s5ev8mADXyFWdOva7kziEt5YUQQqSOzHh8TE97/a8x4uh6nkZHAIlri3/n2oycJjlUTibe17PoSFb4nGLRtWMERYbqttcr6MjwCu642hVVMZ0QmZMUpR8osx504xLimXJuD7MvHwYSF46e36A7lfIUUTmZEEKIrCCzHh/TWnR8HBM9d7Hg2jEgsbPr9DodaVi4tMrJxIeKjo9j063zzL1yhJvPHum2uxUoyVcujaicV75bCfG+pCj9QJn9oLvjzmW++m8d4XExGGsNGF+9Nd0cq6LRaNSOJoQQIhPL7MfHtPAwIpR+B5Zz7pE/ANXz2fNX3S4ydzSLSFAS2OvvxZ+XDnLx8V3d9voFnRjp0pCKeQqrmE6IzEGK0g+UFQ66N54+pN+B5dx6cXavc0lXfq3eRrr9CSGE+GBZ4fiYmq4G36PP/qXcC3+GBg3DKzZgeAV3WVYkC1IUhX0BXky7sJ/LwYG67Q0Ll+LbSk0pbSurHwjxJlKUfqCsctANjYlixH/r2O1/FYDyuQoyv0F3ClnYqJxMCCFEZpRVjo+pYZffFb44spbIuFjMDY2Z5daFRkXKqB1LpDFFUdjjf41pF/ZxNeQ+ABo0dCxRiVEujSgg3XqFSEaK0g+UlQ66iqIw+/JhJp/zIEFRsDExY269rtQqUELtaEIIId4iLi6OSpUqER0drbfdzMyM8+fPA3D06FGmT5/OzZs3yZUrF927d6dv3756+1++fJkpU6Zw5coVzM3NadeuHcOGDcPIKOUjZ7LS8fFjLPE6wfcntwBQ2MKGRe695EpZNpOgJLDb7xpTznno5pyaGBjSr0xtBpdzk+ZWQrzEUO0AQn0ajYYh5etRPndBBh9azZPoCLruWcQPVZrzWZlaMs9UCCEyKF9fX6Kjo5k8eTLFihXTbde+GBp67tw5Bg4cSLNmzfjyyy85e/YsU6ZMQVEUPvvsMwD8/Pzo3bs3Li4uzJgxg1u3bjF9+nTCwsIYN26cGi8rU1MUhT8vHuD383sBcM1blAUNekin+2xIq9HSvFhZGhcpzZobnkw7v4+HkaH8ffkQK6+f5quK7vQsVV2W5xMCuVL6wbLqmeCA0BA+O7Ccay+Gm3RwqMTEmp+QQ+aZCiFEhrNt2za++eYbzp07R44cya+69O7dm4iICP7991/dtt9//51///2XY8eOYWxszHfffcexY8fYs2cPxsbGAKxatYoJEyZw8OBB7OzsUpQpqx4f34eiKPxyZgfzrx4FEhvdzG/QjRyGxionExlBeGw086/+x9zLRwiPiwHAydqO8dVayeg0ke1lu1n2//33H+3bt6dly5YMHDiQsLAwtSNlKIUtbdnSYhBt7CsAsP7WOTrsmse9sKfqBhNCCJGMl5cXRYoUeW1BGh0djaenJ40bN9bb3qRJE54/f865c+cAOHbsGPXr19cVpABNmzYlPj6eo0ePpu0LyEIUReH7k1t1BWlr+/IsdO8hBanQMTcyYUTFhhzt8DVdHauiQYPP0yA6eyxgwMGV3A17onZEIVSTrYrSu3fv8vXXX/P777+zfft28uXLx7x589SOleHkeNGM4TvXZmg1Gi4+vkvzbbM4HXRH7WhCCCFe4uPjg7GxMZ999hkuLi5UqVKFcePGERYWRkBAALGxsdjb2+vdp2jRokDi0N/IyEju37+fbB9bW1ssLCzw9fVN9pyurq5v/QkNDU27F5xBKYrChDM7Wep9AoCujlX5q24XjA1klpRILk8OS6bUasfOVkNxzZv4/8cddy7jtvEPZlzYT0x8nMoJhUh/2aoo3bt3L40aNaJ48eIAfPXVV3z66acqp8qYNBoNg8q5saxRH3Iam/I4KoxOu+az3Puk2tGEEEK84O3tjb+/P25ubsyfP5/Bgwezfft2Bg0apCsOLSz05zKam5sDEBYW9sZ9kvaT0UTvZ+r5vcy7+h+QuLzapJptZckX8U7lchdkU/OB/Fm3M3Y5LImOj2Pq+b003fonnkF+ascTIl1lq1N4fn5+GBoa8sUXX+Dn50eZMmUYM2aM2rEytHoFHdneaij99i/H52kQY05s5krwPcZXb42JnAEWQghVTZ8+nZw5c+Lk5ARAlSpVyJUrF19//TXHjh0DeGOzOq1WS1JbidftoyiKrmHSy941VzRpTml28c/V/5h58QAAbYtXZErNdmg1UpCK96PRaGjn4EKTImWYdn4f/1w7yvWnD2m7cw49S1VndOWmWBmbqh1TiDSXrT414+LiOHToEGPGjGHTpk2Ym5szadIktWNlePZWudnScjBNizgDsPL6aTrv/oegiOcqJxNCiOytatWquoI0Sb169fR+f/VqZ9LvlpaWuiukr7siGhERgaWlZSqmzXq237nM+NM7AWhaxJkZdTrKFVLxQcyNTPihagt2tBxKuVwFAVjmfZL6G/9g550rSF9SkdVlq0/OPHnyUKVKFfLnz49Wq6VVq1ZcvHhR7ViZgoWRCfMbdONrl0YAeD70o/m2WZx75K9yMiGEyJ6Cg4NZt24dAQEBetujoqIAyJUrFwYGBvj7639OJ/1ub2+Pubk5dnZ2+PnpDxUMDg4mLCws2VxT8X+ng+7w5ZG1KCi45i3KX25dZGkP8dHK5S7ItpaD+bFqC3IYGhEUGUr/gyvod2A598KfqR1PiDSTrYrSBg0acPLkSYKCggDYv38/zs7OKqfKPLQaLV9WdGexe08sjEwIinhOh53zWHsj+7X9F0IItWk0GsaNG8eKFSv0tu/cuRMDAwNq1qyJq6sre/bs0bvK4uHhgaWlJWXLlgWgVq1aHDx4kJiYGL19DAwMqFq1avq8mEwmIDSEz/YvIzo+juJWuVns3lOWThOpxlBrwOfOdTjQdgT1CyaOhPDwv0bDzdPZcOu8XDUVWVKmLEq9vLxwdnbmwYMHyW7bvn07LVq0oHz58jRr1ozNmzfrbitXrhzffPMN/fr1o1mzZly/fp3Ro0enY/KsoVGRMmxvOYTiVrmJSYhn5NH1fH9yC7EJ8WpHE0KIbMPW1pZu3bqxfPly/vrrL06cOMGsWbOYOnUq3bp1o2jRogwaNIhz584xYsQIDh8+zIwZM1i4cCEDBgzQLSPTr18/Hj16RP/+/Tl48CCLFy9m4sSJdOrUiQIFCqj8KjOeyLgY+h1YzpPoCGxMzFjeuA82puZqxxJZUGFLW5Y16s1st0/JZWrO85govjyylgEHVxISFa52PCFSlUbJZKdbbt++Te/evQkKCuLw4cPky5dPd9uuXbsYMWIEPXv2pE6dOuzbt481a9Ywc+ZMmjZtmqLneVejhtDQUCwtLbPl4uBJnsdEMezwGvbf9Qagmp098+p3I3eO5F0chRBCpL7Y2FiWLFnChg0bCAwMxM7Ojk6dOtGvXz9dk6K9e/fy559/4uvri52dHd26daNv3756j+Pp6cmUKVPw8vLCxsaGtm3bMmzYMIyMUn71L+n4mRWPj4qiMPTIGrbcvoiBRsuqxn2pVaCE2rFENvA4Moxvj2/Ew/8aAHlyWPB7rfY0LFxa5WRCpI5MU5TGxcWxdu1a/vjjD4yMjHj69GmyorRRo0aULVuW6dOn67YNHz4cHx8fdu3alaLnk6L0/SQoCfxxfp+u82AB85wsaNCD8rkLqZxMCCGEGrJyUbrU6wTfndwCwE9VW9LPubbKiUR2oigK626eZdypbYTFRgPwqWMVfqzaEgsjE5XTCfFxMs3w3bNnzzJ16lT69u3LqFGjkt0eEBCAv78/jRs31tvepEkTbt++nawRxLt4enq+9Uc6EibSarR8Xakx8+p3w8zQmHvhz/hk51w23bqgdjQhhBAi1Xg/ecD4MzsAaFO8Ap+VqaVyIpHdaDQaOpV0ZV/b4dTIVxyA1dfP0HjzTGk8KTK9TFOUOjg4sG/fPoYOHYqBQfLudrdv3wZI1imwaNGiAPj6+qZ9yGysRbFybG05mKKWtkTHxzHsyBomnd1NgpKgdjQhhBDio0TGxTLk0Gqi4+MoYmHLpBqfvHH9VyHSWiELG9Y27cdPVVtiYmCIf1gI7XbMZd6VI/K9S2RamaYozZ07N7ly5Xrj7aGhoQC6NdeSmJsnNh943RpsInWVssnH9lZDqZ0/cX7NrEuH+Gz/ckJjolROJoQQQny4Xz134vM0CAONlr/cumBpbKp2JJHNaTVa+jnXZmerYTha5yVOSeCXMzvps2+ZNEESmVKmKUrfJWlq7KtnLpO2a2Ux63SR1ImwT+maAOwN8KLtjjn4hQarnEwIIYRIuf0B3izxOgHASJeGVM5bROVEQvyfk40dO1oNpXPJxLnc++9603jLTE49kBGCInPJMpVa0hzPV6+IhoeH690u0p6R1oBfqrdmcs12GGkN8HkaRIttf3Ps/i21owkhhBDvLTQmitHHNwGJHeaHlKunbiAhXiOHoTF/1O7AzLqdMTM05kHEczrt/oe/Lh6U4bwi08gyRWnSXFJ/f/2J3n5+fnq3i/TTzakqa5r0w9bEnKfREXT1WMjSF2ebhRBCiIxu8jkP7kc8w9TAiD9qd8BARl2JDKy9gwu7Ww+jjG1+4pUEJp/z4LP9y3ku06hEJpBlPl2LFi1KoUKF2L17t972PXv2UKxYMVkAXCXV8tmzs/VQStvkI15J4LuTWxhzfBMx8XFqRxNCCCHe6OxDP5Z6nQTgK5eGFLN6c18LITKK4jnzsLXFYLo7VQMSp1G13DaLG08fqpxMiLfLMkUpwJAhQ9i+fTvjx4/nyJEj/PTTT+zatYsvv/xS7WjZWiELGza3GESzos4ALPc5Rdc9C2UivhBCiAwpJj6Ob45tREHB2TY//WU9UpGJmBoaManmJ/xeqz3GWgNuP39My22z2OV3Re1oQrxRlipK27Vrx88//8zRo0cZMmQIp0+fZvLkyTRv3lztaNmeuZEJ8+p3Y0RFdwBOPvClxbZZXAu5r3IyIYQQQt+ia8fxeRqEVqPh91rtMdQmX4pOiIzuU8cqbGg+kHxmVoTHxfD5gRX8fm4P8Qkyz1RkPBolqT2tSBFX18QuZ56enionyXy237nMiP/+JTIuFjNDY/6s25mmL66iCiGEyNwy+/ExOCqM2ut/JzQ2mt6lazChehu1IwnxUR5FhjLw4CpOBSV25HUvVIq/632KhZGJysmE+L8sdaVUZA4ti5Vjc/NBFDS3JiIuhn4HljPzwn7k/IgQQgi1TTu/j9DYaHIamzKyYkO14wjx0fLksGRN03665fr23/Xmkx1zCAx7qm4wIV4iRalQhXOuAuxoNZQqeYsC8Pv5vQw+tJrIuBiVkwkhhMiurj8NYoXPaQCGV3THxtRc5URCpI6k5fqm1GqHoUaL15MHtNw+iwuPAtSOJgQgRalQUe4cFqxt+jldXiz4vO3OJdrtnMc9OXMnhBBCBRPO7CReScDeKje9StVQO44Qqa6rY1VWNO6LlbEpjyLD6LBrPtvvXFY7lhBSlAp1GRsY8nut9oyv1goDjZbLwYG02D6Lsw/91I4mhBAiGzn70I8Dd30AGOvaFGMDQ5UTCZE2ahcowdYWgylqaUtUfCwDD65k1qWDMo1KqEqKUqE6jUZD3zK1WNG4DzmNc/AoMoyOu+bz743M2SRDCCFE5vPH+X0AlM9VkKZFpPmeyNpKWOdlW8shVLUrBsCksx58f3KLdOYVqpGiVGQYdQqUZFvLIZTImYeYhHi+Oro+cSiVfEAKIYRIQ2eC7nDk3g0AvnJpiEajUTmREGnP1tSc1U360aZ4BQCWep9k0KFVRMXFqpxMZEdSlIoMpXjO3GxtOYT6BZ0AmHvlCP0OLCcsNlrlZEIIIbKqaRcSr5JWyF0I90KlVE4jRPoxMTDkr7qd+dy5NgA7/a7Qfe8insdEqZxMZDdSlIoMx8rYlCUNe+k+IPcGeNF2xxwCQkNUTiaEECKrOf8ogP/u3QTgq4pylVRkP1qNlh+rtuR71+YAnHzgS/udc3kQ8VzlZCI7kaJUZEgG2sQPyKTW5d5PHtBy+9+cCbqjdjQhhBBZyPwr/wHgbJufBoWcVE4jhHoGlqvLjDqddEvGfLJjDneeB6sdS2QTUpSKDK2rY1VWNfkMaxMzgqPC6bz7H9bfPKt2LCGEEFlAQGgIO/wSl8PoX7auXCUV2V6HEpVY3LAXOQyNCAh7Qvtd87jx9KHasUQ2IEWpyPBq5ndgW8vBugZIw/9bx0TP3SQo0gBJCCHEh1tw7RgJikI+Myta25dXO44QGUL9Qk6sbtwPK2NTgiKe02HXPK4G31M7lsjipCgVmYK9VW62tBiMW0FHAP6+fIjPD6wgXBogCSGE+ADPY6JYc/0MAJ+VqYWR1kDlREJkHK52RVnb5HNsXoxU67R7Puce+asdS2RhUpSKTCOnSQ6WNuxF39I1AfDwv8YnO+cSGPZU3WBCCCEynU23zhMeF0MOQyO6OlZVO44QGU653AVZ32wAdjkseRYTxae7F3DiwW21Y4ksSopSkakYag0YX701E2u0xUCj5VrIfVpun8XZh3L2TgghxPtRFIUVPqcAaGNfgZwmOVROJETG5GRjx/rmAyhobk14XAy99i7h1ANftWOJLEiKUpEp9ShVnZWN+5LTOAePIsPotHs+m25dUDuWEEKITODcowC8njwAoLtTNZXTCJGx2VvlZmPzgRSxsCUiLoaeexfjGeSndiyRxUhRKjKt2gVKsK3lYIpb5SY6Po5hR9Yw5ayHNEASQmRbQ4cOpVGjRnrbGjVqhJOTU7KfkJD/r/18+fJlevTogYuLC7Vr12batGnExsamd/x0s/LFVdKytgWokLuQymmEyPgKWljzb7PPKWSReMW0+95FMsdUpCopSkWmVjxnHra2HEzt/CUA+PPSQQYcXElEbIzKyYQQIn1t2bKFvXv36m0LDw8nICCAkSNHsnbtWr0fKysrAPz8/OjduzcmJibMmDGDvn37snjxYiZOnKjGy0hz4bHRbLtzCYBuTlVlGRgh3lMhCxvWNv2c/GY5CYuNppvHQi48ClA7lsgipCgVmZ61iRnLG/ehV6nqAOzyu0q7nXO5F/5M5WRCCJE+goKC+PXXX8mXL5/edh8fHxRFwd3dnYoVK+r9GBoaAjB//nwsLS2ZPXs2bm5u9O3blzFjxrBmzRqCgoLUeDlpapffVSLjYjExMKS1fQW14wiRqRS1zMW/zT7HzsyK0Nhouu1ZyOXHgWrHElmAFKUiSzDSGvBrjbZMqN4GA42WKyH3aLltFuflDJ4QIhv4/vvvqVWrFjVq1NDb7uXlhYmJCcWKFXvjfY8dO0b9+vUxNjbWbWvatCnx8fEcPXo0rSKrZtOt8wA0LFxaGhwJ8QHsrXKzrunn5H3RlbeLxwKuhcg6puLjSFEqspTepWuwrFFvrIxNeRgZSodd89hy+6LasYQQIs2sW7eOq1ev8sMPPyS7zcfHB2tra7766itcXV1xcXFhxIgRPHr0CIDIyEju37+Pvb293v1sbW2xsLDA1zd5l01XV9e3/oSGhqbNC00FQRHP+e/+TQDaFa+obhghMrHiOfPwb9PPyZPDgmcxkXTbs4g7z4PVjiUyMSlKRZbjVtCRrS0GU8wyF9HxcQw5vJqp5/dKAyQhRJYTGBjIxIkT+fHHH7G1tU12u7e3N48fP6ZkyZLMnTuXMWPGcObMGXr27ElUVJSugLSwsEh2X3Nzc8LCwtL8NaSnbb6XSFAUrE3MqF/ISe04QmRqJazzsrpJP91KCF09FhIU8VztWCKTMlQ7gBBpoYR1Xra1GsKAAys4/uA2My7s5+bTh0yv05EchsbvfgAhhMjgFEVh7NixuLm50aRJk9fu8/3336MoChUqJM6ddHV1xcHBga5du7J161bc3NwAXtvsR1EUtNrk5649PT3fmsvV1TWlLyXd7PS7AkDzomUxNpCvQEJ8rFI2+VjWqDddPBbgHxZCtz2LWN+sP9YmZmpHE5mMXCkVWZaNiRkrm3xGV8eqAGy/c5n2O+fxQM7iCSGygJUrV+Lj48PYsWOJi4sjLi4ORVEAdP8uX768riBNUrlyZSwtLfH29tZdIX3dFdGIiAgsLS3T/oWkk0eRoZx5sbZi06LOKqcRIuuonLco8+t3x1CjxfvJA/rsW0pknKyCIFJGilKRpRlpDZhc8xN+qtoSrUbDpeBAWm6bxZVg6RQnhMjcPDw8ePLkCbVr18bZ2RlnZ2c2b96Mv78/zs7OrFq1ig0bNuDt7a13P0VRiI2NxcbGBnNzc+zs7PDz89PbJzg4mLCwsGRzTTOzPf5eKChYGJlQK7+D2nGEyFLqF3JiRt1OaNBw5qEfAw6uJDYhXu1YIhORolRkeRqNhn7OtVnasDcWRiY8iHjOJzvnssf/mtrRhBDig/3888+sX79e76d+/frky5eP9evX07x5cyZPnsysWbP07rd//36ioqKoWjVxFEmtWrU4ePAgMTH/v7Lh4eGBgYGBbp+sYLffVQDcC5XCRIbuCpHq2havyITqrQE4cNeHkUfX60ZvCPEuUpSKbKN+ISc2txhEIQtrIuNi+Wz/cuZdOSIfmEKITKl48eKUK1dO78fa2hpjY2PKlSuHjY0NgwYNYu/evUyYMIHjx4+zZMkSvv32W9zd3alWrRoA/fr149GjR/Tv35+DBw+yePFiJk6cSKdOnShQoIDKrzJ1hMZEcexF110ZuitE2ulVugZfVWwIwMZb55l6fq/KiURmIacKRbZSyiYf21oO4bP9yzn3yJ9fzuzk1rPH/FqjDUZaA7XjCSFEqurTpw8WFhYsW7aMdevWkTNnTrp06cKwYcN0+zg4OLBo0SKmTJnCF198gY2NDX369NHbJ7M7FHidmIR4TAwMpeuuEGlsREV3giKes/L6aWZePEBhCxu6OFZRO5bI4DSKXCb6IEndBd/VhVBkTJFxsYw6up4tvolrmNbOX4K59btKtzghhPhIGfH4OPLoetbe8MStoCMrG/dVO44QWV5cQjy99i3lcOB1DDValjfuQ50CJdWOJTIwGb4rsqUchkbMcuvC8IruABy9f5M2O+bg+/yxysmEEEKkJkVROBJ4A4B6BeVLsRDpwVBrwNx6XSllk484JYH+B1bg/eSB2rFEBiZFqci2NBoNo1wa8WfdzhhrDbj17BGtts/m5IPbakcTQgiRSm4+e8T9iGcAcqVGiHRkaWzKsoa9sTOzIjQ2ml57lxAky/KJN5CiVGR77RxcWNv0c2xNzHkaHcGnHgtZf/Os2rGEEEKkgsOB1wGwy2GJk7WdymmEyF4KWFiztGEvzAyNCQx/KmuYijeSolQIoIpdMba1GkzJnHmJTYhn+H/rmHzWgwQlQe1oQgghPsKRe4lDd+sWLIlGo1E5jRDZT9lcBZlTr6tuvfhRRzfIygciGSlKhXihqGUuNrcYRN0Xw7v+unSQQYdWyxk9IYTIpGIT4jn5wBeAugUcVU4jRPblXrgUP1RpDsAW34vMvnxY5UQio5GiVIiX5DTJwdJGvenhlLh+3447l+mwaz4PI0JVTiaEECKlrgTfI+LFicUa+YurnEaI7K1fmdp0LFEJgElnPdgf4K1yIpGRSFEqxCuMtAb8VqMtP1VtiVaj4eLju7TcPotrIffVjiaEECIFzgTdAaCopS35zKzUDSNENqfRaJhY4xNc8hRGQWHo4dXcfPpQ7Vgig5CiVIjX0Gg09HOuzWL3XpgbGnMv/Bmf7JgjZ/WEECITOf2iKK2St5iqOYQQiUwNjfinQQ/sclgSGhtN3/3LeBYdqXYskQFIUSrEW7gXLsWmFoMoYJ6T8LgY+uxfyoKrR2WCvhBCZHCKonDm4R0gsZmdECJjyGdmxT/uPTAxMOT288cMPbxGGksKKUqFeJcytvnZ1nIIFXIXIkFR+On0dsae2ExsQrza0YQQQryB7/PHBEeFA1DFrqjKaYQQL6uUpwiTa34CwMFAH2ZePKByIqE2KUqFeA92Zlasb9afFsXKAbDc5xS99i6RISdCCJFBnXnoB4C1iRklcuZROY0Q4lUdSlSmV6nqAEw7v1+3prDInqQoFeI95TA0Zk69TxlWvj6QuPZd2x1z8AsNVjmZEEKIV116fBeAirkLodXI1x0hMqJxVVtSMXdS46M13At7qnYkoRL5lBYiBbQaLd9WbsL02h0x0hpw49lDWm2brevwKIQQImO4+DgQgAq5C6mcRAjxJiYGhsyr3w0bEzOeREcw4NBKYuLj1I4lVCBFqRAfoGPJyqxu8hnWJmaERIfTefc/bLx1Xu1YQgghgNiEeLyeJC7jVS5XQZXTCCHepqCFNX+5dUGDhvOPAvjlzE61IwkVZNuidOnSpbRr107tGCITq56vONtaDqa4VW5iEuL54shapp7fK515hRBCZdefBBH94mpLeblSKkSGV6+gI8MrNgBgsddxtvleUjmRSG/Zsii9du0aCxYsUDuGyALsrXKzteVgauV3AGDGhf0MPbyGqLhYlZMJIUT2dSk4cehublML8ptZqZxGCPE+hldwx61ASQC+ObYB/9AQlROJ9JTtitLw8HDGjRvHV199pXYUkUVYm5ixonFfPnWsAsAW34t02v0PjyJDVU4mhBDZ05UXRWm5XAXRaDQqpxFCvA8DrZaZdTuTJ4cFobHRDDm8Wpbfy0ayXVH6008/0bdvXwoUKKB2FJGFGGkNmFKzHT9UaY4GDece+dNq+994P3mgdjQhhMh2vJ8EAYnrTAshMo/cOSyYWaczAOcfBTD13F6VE4n0kq2K0o0bN2JkZETz5s3VjiKyII1Gw4CydVnQoDs5DI24G/aUT3bM4UjgDbWjCSFEtqEoCtefJhaljjZ2KqcRQqRU3YIlGVzODYDZlw/z3z35HpUdZKuidNu2bVy8eJE2bdrw/fffc+vWLXr37q12LJHFNCnqzMZmA7EzsyI0NpoeexezwueU2rGEECJbeBwVxpPoCACcrPOqnEYI8SG+rtQYlzyJ65d+eeRfHkeGqR1JpLFsVZQuXryYHTt2sGXLFiZMmICDgwNLlixRO5bIgsrlLsi2lkMoY5ufeCWB0cc38cuZnSQoCWpHE0KILO36i6G7Wo0Gh5xSlAqRGRlpDfjb7VMsjUx4GBnKV0fXyeoGWVymLEq9vLxwdnbmwYPk8/W2b99OixYtKF++PM2aNWPz5s3pH1AIoIB5TjY2H4h7oVIAzLtyhP4HVhIZF6NyMiGEyLq8XwzdLWJhSw5DI5XTCCE+VBFLWybVTFy+8cBdH1ZeP61yIpGWDNUOkFK3b99mwIABxMXFJbtt165djBo1ip49e1KnTh327dvHt99+i6mpKU2bNtXbt1q1amzcuPGNz+Pq6vrWHKGhoVhaWn7YixDZhoWRCQvde/Dz6R0s9jrObv+rdNg1n8XuvchrJv/9CCFEaku6Uuok80mFyPTaFK/A3gAvNt++wPjTO6idvwTFrHKpHUukgUxzpTQuLo6VK1fSoUMHoqOjX7vPtGnTaNasGWPHjqVOnTr8/PPPNGvWjJkzZ6ZzWiH+z1BrwC/VW/NztVZoNRouPr4rnXmFECKN6JocWUtRKkRWMKF6a+zMrIiIi2H4f/8SnyBTobKiTFOUnj17lqlTp9K3b19GjRqV7PaAgAD8/f1p3Lix3vYmTZpw+/ZtAgICUvR8np6eb/2Rq6QipT4rU4tF7j0xMzQmMDyxM+/hwOtqxxJCiCzl9vPHADjkzKNyEiFEarA2MWNa7Q4AeD70Y+6VIyonEmkh0xSlDg4O7Nu3j6FDh2JgYJDs9tu3bwNgb2+vt71o0aIA+Pr6pn1IId6hYeHSbGw+QNeZt+feJazwls68QgiRGsJiowmOCgfAXob4CZFluBV0pFepGgBMPb+XayH3VE4kUlumKUpz585NrlxvPsCEhoYCYGFhobfd3NwcgLAwaSUtMoayuQqyveUQnJM6857YxC+nd0hnXiHERxs6dCiNGjXS23b06FHat29PhQoVaNCgAYsWLUp2v8uXL9OjRw9cXFyoXbs206ZNIzY2Nr1ipxq/58G6fxe1lKJUiKzkO9dm2FvlJjYhni+O/Et0fPL+MiLzyjRF6bsktYnWaDSv3a7VZpmXKrKA/C868zYs/KIz79X/pDOvEOKjbNmyhb179+ptO3fuHAMHDqR48eL89ddftGrViilTprBw4ULdPn5+fvTu3RsTExNmzJhB3759Wbx4MRMnTkzvl/DR7oQmFqXmhsbkMjVXOY0QIjWZGRkzo04ntBoN3k8eMPXc3nffSWQaWaZSS5rj+eoV0fDwcL3bhcgozI1MWNigJ31L1wTQdeYNiniucjIhRGYTFBTEr7/+Sr58+fS2//nnn5QpU4bff/+dunXrMmLECD777DPmzp1LTEziSbD58+djaWnJ7NmzcXNzo2/fvowZM4Y1a9YQFBSkxsv5YH6hIQAUtcqV7CS1ECLzq5y3CEPL1QNg3tUjXHiUsp4xIuPKMkVp0lxSf39/ve1+fn56twuRkRhotYyv3ppfqrXW68zrFSKdeYUQ7+/777+nVq1a1KhRQ7ctOjoaT0/P1zYAfP78OefOnQPg2LFj1K9fH2NjY90+TZs2JT4+nqNHj6bPC0gld14M3y1qaatyEiFEWhle0R0nazsSFIVRxzYQI8N4s4QsU5QWLVqUQoUKsXv3br3te/bsoVixYhQoUEClZEK8W58yNVns3gszQ2PuhT/jk51zOHjXR+1YQohMYN26dVy9epUffvhBb3tAQACxsbFvbQAYGRnJ/fv3k+1ja2uLhYXFa5sEurq6vvUnqceDGvxCk4pSmU8qRFZlbGDI1NoddMN4Z18+rHYkkQqyTFEKMGTIELZv38748eM5cuQIP/30E7t27eLLL79UO5oQ7+ReuBSbmg8kn5kVYbHR9N63lGXeJ9WOJYTIwAIDA5k4cSI//vgjtrb6VwffpwHgm/ZJ2i+zNQkMCH0CyJVSIbI6lzyF+bxMbQBmXjygW59YZF6GagdITe3atSMmJoZFixaxbt06ChcuzOTJk2nevLna0YR4L865CrCt5RD67FvKlZB7jD2xmTvPH/Oda3MMpFmXEOIliqIwduxY3NzcaNKkyWtvh+QNAJNotdq37qMoymubBHp6er41l6ur6zuzp4UEJYEHEc8AKGBurUoGIUT6GVWpEbv9r+EXGsyooxvY1HygfFfKxDLlX65du3b4+Pgka+gA0KVLF/bs2cPly5fZuXMnbdu2Tf+AQnyE/OY52dB8AI0KlwZg/tWj9D+4gohY6cwrhPi/lStX4uPjw9ixY4mLiyMuLk5XZMbFxb2xAWDS75aWlrorpK+7IhoREZGpmgSGREUQkxAPJH6OCiGythyGxkyp1Q6Ac4/8Wex1XOVE4mNkyqJUiKzO3MiEBQ168FmZWgB4+F+jw655PJDOvEKIFzw8PHjy5Am1a9fG2dkZZ2dnNm/ejL+/P87Oznh6emJgYJCsAWDS7/b29pibm2NnZ6drCpgkODiYsLCwTNUk8H74M92/pSgVInuold+Bbo5VAZh8zgP/Fx24ReYjRakQGZSBVsvP1VoxoXobtBoNl4IDabXtb66F3Fc7mhAiA/j5559Zv3693k/9+vXJly8f69evp2nTpri6urJnzx7dFVRILGYtLS0pW7YsALVq1eLgwYO6JWKS9jEwMKBq1arp/ro+VNLQXVMDI6yNc6icRgiRXr6r0hw7Mysi42IZfXyT3uedyDykKBUig+tdugaL3XthbmjM/YhnfLJjDgekM68Q2V7x4sUpV66c3o+1tTXGxsaUK1cOGxsbBg0axLlz5xgxYgSHDx9mxowZLFy4kAEDBpAjR2Lh1q9fPx49ekT//v05ePAgixcvZuLEiXTq1ClTda6//2IkSX7znLJGqRDZiJWxKRNrtAXgyL0bbPO9pG4g8UGkKBUiE3AvXIpNLQaS3ywn4XEx9N63RDrzCiHeqUaNGvz111/cunWLIUOGsG3bNr755hs+//xz3T4ODg4sWrSIiIgIvvjiCxYvXkyfPn347rvvVEyecknDd/ObWamcRAiR3hoXKUPjF704fj69ndCYKJUTiZTSKHKN+4MkdRd8VxdCIVLTg4jn9Nm3lMvBgQB87lyb76UzrxAiA1Hr+Dj8yL+sv3WOdg4u/Fm3c7o+txBCfXfDnlB/0zQi42LpW7om46u3VjuSSAH5JitEJpLPzIoNzQbozgb+c/Uonx9YLp15hRDZ3v2IpCul0uRIiOyokIUNwyu4A7DE+wSXHweqnEikhBSlQmQyZkbG/NOgB587Jy4avSfAi/bSmVcIkc09jAgFwM4s8yxjI4RIXZ8718bROi8JisKYE5uJT0hQO5J4T1KUCpEJGWi1/Fi1Jb/VaItWo+GyrjPvPbWjCSGEKh5HJa61mieHFKVCZFfGBob89qLp0YXHAay6flrdQOK9SVEqRCbWs1R1ljTs/VJn3rnsD/BWO5YQQqSruIR4nkZHAmBraq5yGiGEmqrnK07HEpUAmHR2N48jw1ROJN6HFKVCZHINCjmxqcUgCpgndubts38pS7xOqB1LCCHSzZPoCBQS+zbmkqJUiGzv+yrNyWmcg2cxUUw+56F2HPEepCgVIgsoY5ufbS2HUD5XQRIUhe9PbuHHU9tkLoUQIlsIjgrX/VuKUiFELlMLvq7UGIA11z259PiuyonEu0hRKkQWYWdmxfpmA2hSpAwAC68dk868Qohs4eWi1MbETMUkQoiMortTVUrZ5ENB4cdT25BVMDM2KUqFyELMjIyZX787/V/qzNth1zyCpDOvECILC34xZ8zGxAxDrYHKaYQQGYGh1oCfq7YE4MxDP7b4XlQ5kXgbKUqFyGIMtFrGVW3Jr9XboNVouBQcSOvts/F+8kDtaEIIkSaSrpTK0F0hxMtqFShBs6LOAEw4s1NGj2VgUpQKkUX1Kl2Dxe69MDM0JjD8KZ/smMN/926oHUsIIVJdcLQUpUKI1/uhSgtMDAx5EPGcvy8fUjuOeAMpSoXIwtwLl2Jj8wHYmVkRGhtNjz2LWXP9jNqxhBAiVT2JigDAWuaTCiFeUcTSloFl6wIw98oRAkJDVE4kXkeKUiGyuLK5CrKt5RBK2+QjTklg1LENTD7rQYIinXmFEFnD85jENUqtTXKonEQIkRENKVeP/GY5iY6P4zfP3WrHEa8hRakQ2UAB85xsbD4Qt4KOAPx16SDDDq8lKi5W5WRCCPHxQmOjALA0MlU5iRAiIzIzMma0a1MAtt25xLlH/ionEq+SolSIbMLS2JQlDXvRzbEqAFt8L9J1z0KevLSUghBCZEbPYxKLUitjuVIqhHi9T4pXoFyugkBi0yNZIiZjkaJUiGzESGvApJqf8J1rMwBOB92h9Y45+D5/rHIyIYT4cM+iE4fvWhnLlVIhxOtpNVq+f+n7j4f/NZUTiZdJUSpENqPRaBhUzo259bpiYmCI7/PHtN4+mzNBd9SOJoQQHyRp+K4UpUKIt6lVoAQNCjkB8JvnLmIT4lVOJJKkSVEaEiJdrYTI6Fral+ffpp9ja2LOk+gIungsYOttWVhaCJH5yPBdIcT7+s61OVqNhtvPH7PK57TaccQLhh9yp82bN7Nnzx4iIiJISPh/B8/4+HjCw8O5efMmV65cSbWQQoi0UTlvUba2HEzPvYu5/fwxgw+vxj/sCUPKuaHRaNSOJ4QQ7xSfkEBYbDSQOHdeCCHexsnGjs4lXVl9/QzTLuyjnYOLfHZkACm+UvrPP/8wZswYjh49yvXr1zlz5gwBAQFcvnyZs2fPcufOHXr06JEWWYUQaaCYVS62tBxMNbtiAEw6u5tvjm+UIS1CiEwhaeguQE75YimEeA+jXBqRw9CI4Khw5lw+rHYcwQcUpRs3bqRUqVIcP36ctWvXoigKy5Ytw9PTk3HjxhEdHU2FChXSIqsQIo3YmJixqkk/2havCMDq62fovXeJbkicEEJkVKEvfU7J1Q4hxPuwM7NiQNm6AMy/epT74c9UTiRSXJQGBgbSpk0bLCwsKFy4MDlz5sTT0xMDAwO6du1K8+bNWbp0aVpkFUKkIRMDQ/6q25kvKzQA4PC9G3yyYw6BYU/VDSaEEG/x8pVSK1mnVAjxngaWrUtuUwui4mOZdmGf2nGyvRQXpYaGhpibm+t+L1q0KD4+Prrfq1Wrxp07d1IlnBAifWk0Gr6u1JhptTtgqNHi8zSI1tv/5vLjQLWjCSHEa4XHxuj+bWZkomISIURmYmFkwoiK7gD8e+Mst5/J8nhqSnFR6uDgwPnz53W/29vb6zU1ev78OTExMa+7qxAik+hU0pUVjftiZWxKUGQo7XfNY1+Al9qxhMh0Nm/enKL979+/z2effZY2YbKoyLhYAAw0Woy1BiqnEUJkJp86VqGwhQ3xSoJcLVVZiovSdu3asXHjRkaNGkVERAQNGjTA09OTWbNmsXPnTpYsWUKpUqXSIqsQIh3VLlCCzS0GUcjCmoi4GPruX8YSrxNqxxIiUxk7dixr1qx5r303btxI69atOX78eBqnyloi4hJPhOcwNJKu4UKIFDE2MNRdLd1y+yJeIQ9UTpR9pbgo/fTTTxkwYACHDh3C0NCQxo0b06JFC2bNmsVXX31FVFQUo0aNSousQoh05mhtx9YWQ6iQuxAJisL3J7fw8+ntxL+0FJQQ4s2KFCnCzz//zOLFi9+4z+PHjxk4cCDfffcdUVFRDBkyJB0TZn5JV0rNDI1VTiKEyIzaObhQImceFBSmnt+jdpxsK8VFKcCIESM4efIkxsbGaDQa/vjjD1asWMGsWbPw8PDAxcUltXMKIVSS18yS9c3606RIGQD+uXqU/gdXEBknw/SFeJdVq1ZRunRppkyZwqxZs5LdvmPHDlq2bMmhQ4coX748mzZtYujQoSokzbxevlIqhBApZag1YJRLIwA8/K9x/lGAyomyp/cuSmNjY/Hy8uLy5ctERkZiaGiod7urqysNGzYkV65cqR5SCKGuHIbGzK/fnX5lagGJH9odds3nUWSoysmEyNhsbW1Zvnw5rq6u/P3330yZMgWAJ0+e8OWXXzJq1Ciio6MZPXo0a9asoUSJEil6fEVRWLJkCU2aNKF8+fK0bt2abdu26e3TqFEjnJyckv2EhITo9rl8+TI9evTAxcWF2rVrM23aNGJjYz/+DUgHkbqiVK6UCiE+TPNiZXG2zQ/AlHMeKqfJngzfvQssWbKEv//+m7CwMACMjY3p2rUrI0eOTFacCiGyJgOtlp+qtaKoZS5+PL2Ni4/v0nr7bJY26o2jtZ3a8YTIsMzNzVm4cCHDhw9n8eLFBAQEcP78eR4/fkzNmjUZP348hQoV+qDHnjdvHn/++SfDhg2jYsWKHDlyhFGjRmFgYEDz5s0JDw8nICCAkSNHUrVqVb37WllZAeDn50fv3r1xcXFhxowZ3Lp1i+nTpxMWFsa4ceM++vWntaThu3KlVAjxobQaLd9UakKvfUv4795Njt+/Rc38DmrHylbeWVFu3ryZSZMmUbBgQdq0aYNWq+XUqVMsWbKE+Ph4xo4dmx45hRAZRJ8yNSlsacOgQ6sICHtC2x1z+Kd+d2oVSNkVHiGyE2NjY2bNmsV3333Hpk2b0Gq1TJgwgQ4dOnzwY8bGxrJo0SI+/fRTBg0aBECNGjW4cuUKK1asoHnz5vj4+KAoCu7u7jg4vP4L1vz587G0tGT27NkYGxvj5uaGqakpEyZMYMCAAdjZZeyTTklXSmVOqRDiYzQo5IRr3qJ4PvRjyrk9bGo+UJqnpaN3Dt9dtWoVFStWZOfOnXz//feMHTuWzZs306RJE9auXSvLvwiRDTUsXJqNzQZil8OS5zFRdNuziHU3zqodS4gMTavVMnHiRPr06UNCQgLHjx8nLi7ugx/PwMCA5cuX079/f73tRkZGREdHA+Dl5YWJiQnFihV74+McO3aM+vXrY2z8/6KuadOmxMfHc/To0Q/Ol14ipCgVQqQCjUbDN5UaA+D50I/D926onCh7eeeV0lu3bvHVV19hYvL/Bak1Gg29e/fGw8OD27dvyxIwQmRD5XIXZGvLIfTcuxifp0GMOLqOO6HBjHJpJGcWhXjh3r17ybb16NGDiIgI/v33X8LDw/nhhx/QavXPERcoUOCdj63VanFycgIS55YGBwezceNGjh8/zvjx4wHw8fHB2tqar776imPHjhEfH0+9evUYO3YsefLkITIykvv372Nvb6/32La2tlhYWODr65vseV1dXd+aKzQ0FEtLy3fmTy0yfFcIkVpq5negVn4Hjt2/xfTz+3ArUFK+06STdxalkZGRrz24FCpUCEVReP78eZoEE0JkfAUtrNnUYhADD67kyL0bzLx4gICwJ/xeqz0mBjLfXIgGDRq88QuNoigcPnyYI0eO6G3XaDRcu3YtRc+zZ88evvjiCwDq1atH69atAfD29ubx48eULFmSHj16cPv2bf7880969uzJpk2bCA1NbFZmYWGR7DHNzc11vSQysghpdCSESEXDK7pz7P4tzj7y59j9W9SW6Unp4p3fGhMSEl57QDUwMAAgPj4+9VMJITINK2NTljbqzdgTm1l9/Qwbb50nMOwpC9x7YGNipnY8IVTVtm3bdDnLXqZMGVasWIGPjw8zZ86kf//+LFu2jO+//x5FUahQoQKQeJXTwcGBrl27snXrVtzc3ABem1FRlGRXcAE8PT3fmuVdV1JTW9KVUlMDuVIqhPh4NfIVp5qdPaeCfJl+YR+18jvI1dJ0IJcyhBAfzUhrwJSa7ShqmYtJZ3dzKsiXtjvmsLRhb4pZyTJRIvuaNGlSujxP4cKFKVy4MFWqVMHCwoJvv/2W8+fPv3bd8MqVK2NpaYm3tzctWrQAeO0V0YiIiHQdhvuhYhMST46byugMIUQqGVHRnS4eCzgVdIcTD25LJ9508F7rlD59+pR79+7p/Tx48ACAkJCQZLe9bg6NECJr02g0DC1fj9lun2JiYMitZ49ovX02Zx/6qR1NiCzp6dOnbN68maCgIL3tZcqUAeDOnTts2LABb29vvdsVRSE2NhYbGxvMzc2xs7PDz0///6fBwcGEhYUlm2uaEcUkJDaLMnoxgksIIT5WrfwOVMlbFIAZF/arnCZ7eK/Tir/99hu//fbba28bNWpUsm0fMh9GCJE1tC5egfzmOem7fxkh0eF02v0PM+t2pmWxcmpHEyJLSUhIYPTo0QwePFg3nxQSu+kCVKhQgS5dulC1alVmzZqlu33//v1ERUXp1i2tVasWBw8e5JtvvtF14PXw8MDAwCDZ2qYZUcyLaURGWilKhRCpQ6PRMLyiO932LOL4g9uceuBLtXwZ/yRdZvbOovSTTz5JjxxCiCykil0xtrYcTM+9S/B9/piBB1fynWszBpatK/MyhEgltra2dO3alfnz52Nqakq5cuU4e/Ys8+bNo2PHjhQvXpxBgwYxadIkJkyYQIMGDbh+/Tp//fUX7u7uVKtWDYB+/fqxY8cO+vfvT69evbhz5w7Tpk2jU6dO79UFWG1Jw3eluZoQIjXVLVASlzyFOf8ogBkX97M6Xz+1I2VpGkVRFLVDpKeFCxeyceNGtFothQoV4tdff8XW1jbFj5PUyOFdDR+EyM6eRIXz2YHlnA66A0B3p2pMqN4aQ7miIUSqiI2NZcmSJaxfv5579+6RL18+OnbsSL9+/XRNitatW8eyZcvw9/cnZ86ctGrVimHDhmFqaqp7HE9PT6ZMmYKXlxc2Nja0bduWYcOGYWSU8uZB6X18bLrlT66E3OOHKs0ZULZuujynECJ7OHDXh557FwOwufkgXO2Kqpwo68pWRemZM2f46aef+PfffzE3N+ePP/4gJCSEX3/9NcWPJUWpEO8nKi6WkcfWs+X2RQDqF3RiTv2uWBiZvOOeQojMKL2Pj+6bpuPzNIjx1VrRt0ytdHlOIUT2oCgKLbf/zcXHd6lX0JEVjfuqHSnLeq9GR1mFra0tP/74I+bm5kBiMwhpyiRE2jI1NOKvup35onx9AA4G+tBh5zweRMgaxyLra9SoET/++CMeHh48e/ZM7ThZUsyL4bvGMnxXCJHKNBoNwys0AOBQ4HUuPApQOVHWla2KUgcHB13ThrCwMGbPnk2jRo1UTiVE1qfVaPmmchOm1GqHgUbLlZB7tN7+N95PHqgdTYg09dlnn/H06VN+/PFHatSoQfv27fnjjz84ceIEMTExasfLEmLiE7vvGsu0ACFEGmhYuDTOtvkB+PvyIVWzZGXZ8rRiUFAQgwYNwsXFha5du6odR4hso6tjVQqYWzPgwAruhT/jkx1z+KdBD2oXKKF2NCHSRJcuXejSpQuKonD16lWOHz/OiRMnWLZsGQCVKlWiRo0a1KxZk7Jly6qcNnNKanRkpM2WX2mEEGlMo9EwpFw9Bh9ezW6/a9x69giHnHnUjpXlZKsrpQDe3t507twZd3d3xo8fr3YcIbKdegUd2dh8IHZmVoTGRtN9zyLW3Tirdiwh0pRGo6Fs2bL079+fxYsXc+bMGWbPno2zszO7d++mU6dOVKtWjS+++IK1a9eqHTdT+f/wXblSKoRIG82LlaWoZS4UFOZcPqx2nCzpnUVpVpoP8+jRI/r06cPIkSMZMmSI2nGEyLaccxVgW4vBlLLJR5ySwIij65h2fh/ZqO+ayOaMjY2pVasWo0aNYuPGjRw/fpyffvoJa2trFixYoHa8TEWG7woh0pqh1oBBL7p7b7h1nvvhmbsmyojeWZRmxPkwXl5eODs78+BB8vlo27dvp0WLFpQvX55mzZqxefNm3W0LFiwgIiKCBQsW0KZNG9q0acPw4cPTL7gQQqeAhTUbmw+kboGSAEy7sI+RR9frvmAKkZ1YW1vTrFkzxo8fz969e9WOk6nohu9KoyMhRBrqUKISeXJYEJsQz4KrR9WOk+W895Iwr86HOXfuHJD+82Fu375N7969CQoK4vDhw+TLl093265duxgxYgQ9e/akTp067Nu3jzVr1jBz5kyaNm2aoudJamn/JqGhoVhaWsqSMEJ8pNiEeEYf38jaF0N4a+cvwfwG3bEyNn3HPYUQGVF6LgmjKAqFl4wB4N+mn1Mzv0OaP6cQIvuadekQk87uxtzQmJOdRmNjYqZ2pCzjveeUqj0fJi4ujpUrV9KhQweio6Nfu8+0adNo1qwZY8eOpU6dOvz88880a9aMmTNnpnoeIUTqMNIaMLVWB0a5JHbCPnr/Jp/smENg2FN1gwmRSkJDQ9WOkGUlzScFMJZGR0KINNazVHUsjUwIj4thmdcJteNkKR/c6Ci958OcPXuWqVOn0rdvX0aNGpXs9oCAAPz9/WncuLHe9iZNmnD79m0CAlK2rpCnp+dbfywtLT/q9Qgh/k+j0TC8ojsz6nTCSGuAz9MgWm//myvBgWpHE+KjVa1aFQ8PD7VjZElxLxWlhtps17tRCJHOrIxN6VGqOgALrx0nMk6W9kotqfYJntbzYRwcHNi3bx9Dhw7F4DUd9m7fvg2Avb293vaiRYsC4Ovrm+qZhBCpq0OJSqxo1AcrY1OCIkNpv3MeB+76qB1LiI+iKIqsSZpG4hISdP+WolQIkR4+K1MLEwNDQqLDWXNdpvGllkzzCZ47d25y5cr1xtuThkdZWFjobTc3NwcgLCws7cIJIVJNrQIl2NR8EAXNrQmPi6HPvqWs8DmldiwhRAaUoPy/KNVqNComEUJkF3ZmVnQsURmAeVeP6JqtiY+TaYrSd0nq16R55aCUtF0rZ1CFyDScbOzY2nIw5XIVJF5JYPTxTUz03K33BVQIIRL4f69GrUaO80KI9DGgbB20Gg13w56y884VteNkCVnmEzxpjuerV0TDw8P1bhdCZA52Zlasb9afBoWcAPj78iGGHV5LtCwZIzKh4OBgEhLkpEpqi0/4f1FqIEWpECKd2FvlpkmRMgDMv/qfrLOeCrLMJ3jSXFJ/f3+97X5+fnq3CyEyD3MjExa596SHUzUAtvhepKvHAp5ER6icTIiUmTx5MpUqVaJz58789NNPrFu3jitXrshc048U/9LoCQMZviuESEf9nesAcPHxXc489FM5Teb3Qf3THz9+zK5duwgMDMTMzIwyZcpQq1YtcuTIkdr53lvRokUpVKgQu3fvplGjRrrte/bsoVixYhQoUEC1bEKID2eoNeC3Gm0pbGnLb567OBV0h7Y75rCsUW+KWr55nrkQGUnFihWJjo7Gy8uLixcvAonTTQwMDHBwcKBMmTK6n1KlSun6IYi3S1Bk+K4QQh2ueYvikqcw5x8FMP/Kf1S1K6Z2pEwtxUWpp6cnn3/+OVFRUXqXqq2trRk8eDA9e/ZM1YApMWTIEMaMGUPOnDmpV68eBw4cYNeuXUyfPl21TEKIj6fRaBhczo1C5tYM/+9fbj17RJvtc1jcsBcueQqrHU+Id+ratSutWrUiLi6Omzdvcu3aNd2Pt7c3Pj4+bNq0CY1Gg1ar5erVq2pHzhTkSqkQQi0ajYb+znUYdGgVHv7X8H3+GHur3GrHyrRSXJROnjwZgN9++40aNWqQkJDAhQsXWLJkCb/99hsXL17kjz/+SPWg76Ndu3bExMSwaNEi1q1bR+HChZk8eTLNmzdXJY8QInW1Ll6BfOY56bt/GY+jwui4az5/u3WhSVFntaMJ8V4MDQ0pVaoUpUqVol27dkBiQ747d+5w7do1rl69ire3t8opMw/9olSulAoh0lezos4UNLcmMPwpC68dY0L1NmpHyrQ0Sgpn5lasWJE+ffrw5ZdfJrtt3bp1jBs3ju+++47u3bunWsiMyNXVFUi8ciyESF+3nj2i597F+IWGoEHDz9Va0rdMLbVjCfFapUqV4vfff6dVq1ZqR0kX6Xl8vPXsEW4bE0+Ee3YeSz4zqzR/TiGEeNn8K/8x/swOchgacabTGKxNzNSOlCml+LSiubn5G+dnduzYkebNm7NmzZqPDiaEEG/ikDMPW1oMxiVPYRQUxp3axs+nt8uSMSJDkiXJ0o4M3xVCqO1TxypYGJkQGRfLSp/TasfJtFJ8pKxWrRqHDx9+4+3Vq1dP1gFXCCFSW+4cFvzb9HOaFkkcuvvP1aMMOLiSyLhYlZMJoc/T05N69eqpHSNLkiVhhBBqszQ2patjFQAWeR0nRpau+yAp/gTv1KkT58+fZ+nSpa+9/e7du+TNm/ejgwkhxLvkMDRmXv1ufPZi6O4uv6t03v0PwVFh77inEOnHzMxM1spOIy+PjtDKlVIhhEr6lqmFVqMhKOI52+9cVjtOppTiorR37948e/aMSZMm0aNHD7Zs2cKtW7e4c+cO69atY9myZfTp0yctsgohRDIGWi0/V2vFT1VbokHDuUf+tN4+h9vPHqkdTQiRxqTRkRAiIyhkYUOLouWAxDmmKWzZI/iA7ruDBw/G29sbb29vzpw5w5kzZ9C8dHaydOnSWFlZcf36dYoXL46h4QcthSqEECnSz7k2hSxsGHp4DX6hwbTeMYfF7j2pIuuGCZFlxSsyfFcIkTH0L1uHbXcucSXkHice3KZmfge1I2UqKa4Yv/jiC92/nz9/rltjzcvLi2vXrnH9+nW+/vprNBoNhoaG2Nvb4+TkxO+//56qwYUQ4lVNizrzb7PP6bNvKcFR4XTxWMDMOp1oaV9e7WhCiDQgw3eFEBmFS57CuOYtiudDPxZdOy5FaQp91GVMKysrqlevTvXq1XXbYmJiuHHjBl5eXrpC9eDBgx8dVAgh3kelPEXY0mIwPfcu5vbzxww8tIrvwp4wsGxdvVEdQojMT+9KqXQ5FkKorG/pmng+9GNPwDXuhj2hkIWN2pEyjVQfW2tsbIyzszPOzrKYvRBCHcWscrGlxSA+O7Cc00F3+NVzFwFhTxhfrRWGWgO14wkhUoksCSOEyEiaFSuL3RkrgiKes9TrJN9VaaZ2pExDTisKIbIkG1NzVjX+jNYvhu4u8z7JZ/uXEx4brXIyIVKPoigsWbKEJk2aUL58eVq3bs22bdv09jl69Cjt27enQoUKNGjQgEWLFiV7nMuXL9OjRw9cXFyoXbs206ZNIzY24y+vlJDw8vBd+UojhFCXkdaAnk7VAFh944wsU5cC8gkuhMiyTA2NmOXWhcHl3ADYf9ebDrvm8zAiVOVkQqSOefPmMWXKFNq2bcu8efOoVasWo0aNYufOnQCcO3eOgQMHUrx4cf766y9atWrFlClTWLhwoe4x/Pz86N27NyYmJsyYMYO+ffuyePFiJk6cqNbLem8JJA7flfmkQoiMoptTNYy1BjyNjmDz7Qtqx8k0pDWuECJL02q0jHVtRhELW8ae3Mzl4EBa7/ibZY364Ghtp3Y8IT5YbGwsixYt4tNPP2XQoEEA1KhRgytXrrBixQqaN2/On3/+SZkyZXTNBuvWrUtcXBxz586lR48eGBsbM3/+fCwtLZk9ezbGxsa4ublhamrKhAkTGDBgAHZ2Gff/J0lzSqXzrhAio8idw4JW9uXZcOs8i64do0tJV+lp8R7kU1wIkS10L1WNxe69MDM05m7YU9rumMOx+7fUjiXEBzMwMGD58uX0799fb7uRkRHR0dFER0fj6elJ48aN9W5v0qQJz58/59y5cwAcO3aM+vXrY2xsrNunadOmxMfHc/To0bR/IR8hQZErpUKIjKdP6ZoAeD15wKkgX5XTZA5ypVQIkW24Fy7FhmYD6LVvCQ8jQ+m+ZxFTa3egvYOL2tGESDGtVouTkxOQOLc0ODiYjRs3cvz4ccaPH09AQACxsbHY29vr3a9o0aIA+Pr6UqFCBe7fv59sH1tbWywsLPD1Tf5lytXV9a25QkNDsbS0/JiXJoQQmVrFPIWplKcI5x75s9jrBNXzFVc7UoYnV0qFENlKudwF2dZyCE7WdsQmxPPlkbXMuLAf5aWlJYTIbPbs2UOtWrX4448/cHNzo3Xr1oSGJs6dtrCw0NvX3NwcgLCwsDfuk7RfWFhYGicXQoisKelq6W6/q9wLe6pumExAilIhRLZT0MKajc0HUvPFmcup5/cy6tgGYhPiVU4mxIcpU6YMK1as4IcffuDcuXP0799fd6LlTXOZtFrtW/dRFAXta9b+9PT0fOuPXCUVQghoUawseXNYEq8ksMznpNpxMjwpSoUQ2VJOkxysaNxXN3R37Q1Peu9dQmhMlMrJhEi5woULU6VKFbp37853333HqVOndAXnq1c7k363tLTUXSF93RXRiIiIDF9gyggHIURGZWxgSDenqgCs8jlDlCwP81ZSlAohsi1jA0Nm1OnE8IruABy+d4N2O+dyL/yZysmEeLenT5+yefNmgoKC9LaXKVMGgLt372JgYIC/v7/e7Um/29vbY25ujp2dHX5+fnr7BAcHExYWlmyuaUalQRodCSEynu5O1TDSGhASHc5W34tqx8nQpCgVQmRrGo2GUS6N+KN2Bww1WryePKD19r+5FnJf7WhCvFVCQgKjR49m7dq1etuPHTsGQLly5XB1dWXPnj16VxQ9PDywtLSkbNmyANSqVYuDBw8SExOjt4+BgQFVq1ZNh1cihBBZk52ZFS2KlQNg0bXjMrrjLaQoFUIIoHNJV5Y16oOFkQkPIp7Tfudcjt67qXYsId7I1taWrl27Mn/+fObPn8+JEyeYNWsW06ZNo2PHjhQvXpxBgwZx7tw5RowYweHDh5kxYwYLFy5kwIAB5MiRA4B+/frx6NEj+vfvz8GDB1m8eDETJ06kU6dOFChQQOVXKYQQmVtSw6MrIfc498j/HXtnXxpFSvYPktQS39PTU+UkQojUdC3kPj33LuZBxHMMNVqm1m5PhxKV1Y4lxGvFxsayZMkS1q9fz71798iXLx8dO3akX79+uiZFe/fu5c8//8TX1xc7Ozu6detG37599R7H09OTKVOm4OXlhY2NDW3btmXYsGEYGRmlOFN6Hh/3B3jTa98STA2MuNnzlzR/PiGESClFUWi+bRaXgwNp7+DCzLqd1Y6UIUlR+oGkKBUi67oX/oweexbh8zRxrt7XLo34okKDN3YxFUL8nxpFaQ5DI270kKJUCJExrbp+mm+ObcRYa4Bn57HYmpqrHSnDkeG7QgjxigLmOdnUYhC18jsA8Pv5vXxzfCNxsmSMEEIIIVKorX1FrIxNiUmIZ+0NuaD1OlKUCiHEa1gZm7K8UR/dkjGrr5+hz75lhMdGq5xMCCGEEJmJmZExHRwqAbDC5zQJSoLKiTIeKUqFEOINkpaM+aJ8fQAOBvrQYdd8HkaEqpxMCCGEEJlJj1LVAfALDeaINFJMRopSIYR4C41GwzeVmzCp5icYaLRcDg6k9Y6/ufH0odrRhMj2FKQthhAicyhpnZca+YoDsMzrhMppMh4pSoUQ4j10d6rGIveemBkaczfsKW13zObkg9tqxxJCABqkCZkQIuPr9eJq6b673gSGPVU3TAYjRakQQrwn98KlWN+sP3lyWPAsJoquHgvZevui2rGEEEIIkQk0KepM3hyWJCgKq66fVjtOhiJFqRBCpED53IXY0mIwJXLmISYhnsGHVzP38hFkdS0hhBBCvI2R1oBPHasAicvExMTHqZwo45CiVAghUqiIpS2bWgyiml0xACZ47uT7k1uJT5BuekKkJzkZJITIbLo5VkWr0fAoMgwP/2tqx8kwpCgVQogPYGNixsrGn9GqWHkAlnqfoP/BFUTGxaicTIjsR2aUCiEyiwIW1jQqXBqAZd4nVU6TcUhRKoQQH8jU0Ii/63VhQNm6AHj4X6PT7n8IjgpTOZkQQgghMqqk5WFOPLjN9adBKqfJGKQoFUKIj6DVaPmhSnN+qdYaDRrOPwqg9fY53H72WO1oQgghhMiA6hYoQVHLXAAs9z6lcpqMQYpSIYRIBX3K1OSfBt0xMTDELzSYtjvmcPahn9qxhBBCCJHBaDVaejhVA2D9zbNExMrUHylKhRAilTQt6sy/Tftja2JOSHQ4nXb/wy6/K2rHEiLLkjZHQojMqnPJypgYGBIaG81m3wtqx1GdFKVCCJGKKuctwpaWgyhmmYvo+Dj6H1jJomvH1I4lRJam0UirIyFE5mJjak7LYuUAWOZ1Mtt3E5eiVAghUpm9VW62tBxEpTxFUFAYd2ob409vJ0GRJWOEEEIIkajni4ZHV0LucfHxXZXTqEuKUiGESAO5TC1Y27QfTYqUAWD+1aMMPrSaqLhYlZMJIYQQIiOolKcIZWzzA7DcJ3s3PJKiVAgh0kgOQ2Pm1+9O79I1ANh+5zJd9yzkSXSEysmEyBqy+3A3IUTmptFodA2Ptty+yNNs/P1AilIhhEhDBlotv1RrzQ9VmgNwOugObXfMwT80ROVkQmQdMqNUCJFZfeLggrmhMVHxsWy4dV7tOKqRolQIIdKYRqNhQNm6zKnXFWOtAbeePaLNjtnZfv6IEEIIkd1ZGJnwiYMLACu8T2XbESBSlAohRDppZV+e1U36kdM4B48iw+iwax77A7zVjiWEEEIIFXV3qgrAjWcPORXkq3IadUhRKoQQ6ahaPns2txhEIQtrIuNi6bN/KSuyeXMDIT5U9ryeIITIasrmKohLnsJA9m14lO2K0l27dtG8eXNatWrFL7/8QlxcnNqRhBDZTEnrvGxtMYRyuQqSoCiMPr6JyWc9su2QHSE+lqxTKoTI7JIaHu28c4XHkWEqp0l/2aooffToEb/++iuLFy9my5YtPHr0iI0bN6odSwiRDeU1s2R9s/40KOQEwF+XDvLFkbXExMuJMiGEECK7aWVfnpzGpsQmxPPvzbNqx0l32aooPXbsGJUqVcLOzg6tVkv79u3Zvn272rGEENmUuZEJi9x70tUxcS7JptsX6LF3Mc9jolROJoQQQoj0lMPQmA4lKgOw0ucUCUqCyonSV7YqSh8+fIidnZ3u97x58xIUFKRiIiFEdmeoNWByzU/4tlITAI7dv8UnO+ZwL+ypusGEEEIIka6ShvD6hYbw372bKqdJX9mqKE1ISNCbd6IoClpttnoLhBAZkEajYViF+syo0wlDjRafp0G02jGbayH31I4mRIamSKsjIUQWUsI6LzXyFQdguXf2aniUrSqyfPny8fDhQ93vjx490rtyKoQQaupQohIrGvfF0siEoIjntNs5jyOBN9SOJTKwhIQEVq9eTatWrXBxcaFhw4ZMnDiRsLD/N8lo1KgRTk5OyX5CQkJ0+1y+fJkePXrg4uJC7dq1mTZtGrGxsWq8pA8ibY6EEFlF0tXSvQFe3A9/pnKa9JOtitJatWrh6elJYGAgCQkJbNiwgXr16qkdSwghdGoXKMHG5oPIb5aTsNhoeu5dzL83PNWOJTKoBQsW8Msvv1CvXj3+/vtv+vTpw+bNm/nyyy8BCA8PJyAggJEjR7J27Vq9HysrKwD8/Pzo3bs3JiYmzJgxg759+7J48WImTpyo5ksTQohsqWlRZ3KZmhOvJLD6+hm146QbQ7UDfAgvLy86dOjA/v37yZcvn95t27dvZ86cOQQEBFCwYEEGDBhA27ZtAciTJw8//vgjAwcOJCYmhsqVK9OtW7fXPoerq+tbM4SGhmJpaZkqr0cIIV5W2jYfW1oOptfexXg9ecBXR9cTGP6U4RXcZekLoaMoCgsWLKBz586MHDkSgJo1a2JjY8OIESPw8vIiMjISRVFwd3fHwcHhtY8zf/58LC0tmT17NsbGxri5uWFqasqECRMYMGCAjCgSQoh0ZGxgSJeSVfj78iFWXT/NFxXqY6g1UDtWmst0V0pv377NgAEDXru+6K5duxg1ahS1atXi77//pmrVqnz77bfs3r1bt0+jRo3Ytm0bHh4e/PbbbxgZGaVnfCGEeC8FzHOysflA6hQoAcAf5/fx9bENxCbEq5xMZBTh4eG0bt2ali1b6m0vXjxxPpK/vz9eXl6YmJhQrFixNz7OsWPHqF+/PsbGxrptTZs2JT4+nqNHj6ZJ9tQiS/sKIbKirk5V0KDhQcRz9gd4qx0nXWSaK6VxcXGsXbuWP/74442F5LRp02jWrBljx44FoE6dOjx79oyZM2fStGnTFD2fp+fbh8u960qqEEJ8LEtjU5Y27M03xzay/tY51tzw5H7Ec+bV74aFkYna8YTKLCws+P7775Nt37dvHwAlSpTg2LFjWFtb89VXX3Hs2DHi4+OpV68eY8eOJU+ePERGRnL//n3s7e31HsPW1hYLCwt8fX2TPX7GHEkkIwiEEFlHUctcuBUsyaHA6yz3OUWTos5qR0pzmeZK6dmzZ5k6dSp9+/Zl1KhRyW4PCAjA39+fxo0b621v0qQJt2/fJiAgIL2iCiFEqjE2MGR6nY4Mr+gOwOHA67TfOZcHEc9VTiYyoosXLzJ//nwaNmyIg4MD3t7ePH78mJIlSzJ37lzGjBnDmTNn6NmzJ1FRUYSGhgKJBe6rzM3N9RomCSGESD9JDY8OB97ALzRY5TRpL9NcKXVwcGDfvn3kypWLjRs3Jrv99u3bAMnO9hYtWhQAX19fChcunPZBhRAilWk0Gka5NKKguTWjj2/iash92myfzYrGfSlpnVfteCKDOHv2LAMHDqRQoUJMmDABgO+//x5FUahQoQKQeJXTwcGBrl27snXrVtzc3ABeO1f5TcumyUgiIYRIe+6FS5HPzIoHEc9Z5XOGMa4pG/WZ2WSaK6W5c+cmV65cb7z9TWd7zc3NAeRsrxAi0/vUsQpLGvbCzNCYwPCntN0xh9NBd9SOJTKAnTt30qdPH/Lnz8+SJUuwsbEBoHz58rqCNEnlypWxtLTE29tbd8x83TEyIiIiwzf0k3VKhRBZlaHWgK6OVQFYc+MMMfHJ++lkJZmmKH0X5UW3g1fP9iZtf93ZXiGEyGzqF3JifbP+5MlhwbOYSD71WMD2O5fVjiVUtHjxYr766isqVqzIypUryZs38ep5REQEGzZswNtbv0mGoijExsZiY2ODubk5dnZ2+Pn56e0THBxMWFhYstFHGZV0pRZCZEWfOlbBQKMlOCqc3X5X1Y6TprJMpZZ0NvfVs73h4eF6twshRGZXPnchtrQYjEPOPETHxzHo4CoWXM3YXVJF2li3bh2TJk2iWbNmLFiwQO9YZ2JiwuTJk5k1a5beffbv309UVBRVqyaega9VqxYHDx4kJiZGt4+HhwcGBga6fYQQQqS//OY5aVS4NADLfU6pnCZtZZmiNOlsrr+/v972pLO/meVsrxBCvI8ilrZsbj4Q17xFUVD46fR2xp/eToKSoHY0kU6Cg4P59ddfKViwIN26dePatWtcuHBB9/Ps2TMGDRrE3r17mTBhAsePH2fJkiV8++23uLu7U61aYhONfv368ejRI/r378/BgwdZvHgxEydOpFOnThQoUEDlVymEENlbj1KJn9UnHtzm5tOHKqdJO5mm0dG7FC1alEKFCrF7924aNWqk275nzx6KFSsmB1YhRJZjY2rO6ib9GHZ4Dbv9rzL/6lEeRDxnep1OmBhkmY938Qb//fcfkZGRBAYG0q1bt2S3T5kyhT59+mBhYcGyZctYt24dOXPmpEuXLgwbNky3n4ODA4sWLWLKlCl88cUX2NjY0KdPH719hBBCqKNOgRIUtbTFLzSE5T6n+LlaK7UjpYks9a1lyJAhjBkzhpw5c1KvXj0OHDjArl27mD59utrRhBAiTeQwNGJe/W78eHobS7xOsNX3Eg8jQ1nQoAfWJmZqxxNpqG3btrRt2/ad+3Xs2JGOHTu+dR9XV1f+/fffVEqWfqTNkRAiq9NqtHRzqsZvnrtYf/Msoys3IYehsdqxUl2WGb4L0K5dO37++WeOHj3KkCFDOH36NJMnT6Z58+ZqRxNCiDRjoNXyS7XWfOfaDICTD3xpv3Me98KeqhtMiHQibY6EEFlZ55KVMdIa8Cwmim2+l9SOkyYyZVHarl07fHx8yJcvX7LbunTpwp49e7h8+TI7d+58r7PIQgiR2Wk0GgaVc+Ovul0w0hrg8zSIVjtmcy3kvtrRhBBCCPERcpla0LxYWSDrNjzKlEWpEEKI1/vEoSIrGvfF0siEoIjntNs5l6P3bqodSwghhBAfoYdTYsOj848CuBIcqHKa1CdFqRBCZDG18juwsfkg8plZERYbTY+9i9l467zasYRIfYrMKhVCZA/V7OwpmTNxHeoVPqdVTpP6pCgVQogsqLRtPra2GIyTtR2xCfF8cWQtsy4dQpEv8SIL0sisUiFEFqfRaOj+YnmYTbfOExYbrXKi1CVFqRBCZFEFLKzZ2HwgNfIVB2DS2d18f3IL8QmylqkQQgiR2XRwqISpgRHhcTFsymIjoKQoFUKILCynSQ5WNO5LG/sKACz1Pkn/gyuIjItVOZkQQgghUiKnSQ7aFE88ni/3OZWlRj9JUSr+196dx9WU/38Af91WqovKviVRVFpIxpoSwjQMYxiEjH03zIiHdYbJPs0YISlfy9iHMURhshWjlBlDxUyosZOo0Pr5/WG6v64W7ae6r+fjcR/f6ZzPOff9eX+vTu/7OZ/PIaIqTltdA+sdhmCiZTcAQGDcDQw9sQUJb1IkjoyIiIiKInvBoxsJDxDxJF7iaEoPi1IiIhWgJlPDgvZ98XUHV8ggw5Unceh/bCPuJj2TOjSiYqs6YwRERIVjXbsx2hg2AgDsiL4kcTSlh0UpEZEKGWPeGZsdh0NbXQO3Xz5F/6Mb8cfTf6UOi6hEZFzniIhUhEwmw4j/Rkt/vfMnnqe+kjii0sGilIhIxfRtZok9vceilrYOnr5JxuDjPvjt3xipwyIiIqJCGNDcGnqa2kjNzMCBv69IHU6pYFFKRKSC2tdrhsN9J6KJnj5eZaTB/dT/sPtmmNRhERER0XvoampjkIktAGBHdNVY8IhFKRGRimpRqy5+6TcZlgYNkSmy8GXIQayNPFklLm6kGvhZJSJVNfy/W3hjXz5F6MNYiaMpORalREQqrK6OHAf6ToBDI1MAwHdXT2NOyEGkZ2VKHBlR4cnASaVEpFrMDRrArq4RAGBn9O8SR1NyLEqJiFScnqY2tjmPwpCW7QAAe2+Fw/3U/5CSnipxZERERJSf7AWPjt/9C09eJ0kcTcmwKCUiImiqqWNN508w06YHAODMvZsYfNyn0l/kiIiIqqoPm7VBLW0dZIgs7L0VLnU4JcKilIiIALxdZn6ObU+s7DQQ6jI1/PnsHvof3Yh/XjyROjSiPHFGKRGpsmoamhjS4u1dTrtiLiMzK0viiIqPRSkRESkZbmYPvx4jUV1DE3HJCRhwbCOuPL4rdVhE+eJzSolIVQ03swcAxCc/x9n7tySOpvhYlBIRUS49mrTCfpfxMKymi+epr/DpiS04cfe61GERERFRDs1r1kGXBi0AADuiL0kcTfGxKCUiojzZ1GmCX/pNhnGN2kjNzMC433ZiW9RFqcMiIiKiHEa0ervg0el/o3EvOVHaYIqJRSkREeWrWQ1DHO43EbZ1mkBAYMGlX/Bt+HFkico7b4WIiKgq6d3UHHWry5ElBH66eVnqcIqFRSkRERXIsJoe9rmMQ68mrQEA3tfOYsa5fUjLzJA4MlJ1gksdERFBU00dQ1vaAQD23AyrlM8aZ1FKRETvVV1DCz5OI+D23zPRDsVehdtJf7xMeyNxZESADFzpiIhU2zAze8ggw6PXSTgZFyV1OEXGopSIiApFQ00d33YcAI92LgCAkAf/YGDAJtxPeSFxZERERKqtsZ4+nBqbAQB2xvwucTRFx6KUiIgKTSaTYapVd3h1/RQaMjVEP3+I/ke9Ef38odShERERqTS3/xY8Onf/Fm6/fCpxNEXDopSIiIrskxZtsaOXO/Q0tfHg1QsMDNiE0Af/SB0WqRjBKaVERAqOjczQSLcWAGBXTOVa8IhFKRERFUvXhi1xsM8E1Ksux8u0NxgR5IdfYv+QOiyVkpWVhd27d8PV1RW2trZwdnaGp6cnkpOTFW0uXLiAQYMGwdraGk5OTvDz88t1nmvXrsHNzQ22trbo0qUL1q1bh/T09PLsSolwRikREaCupobhZvYAgL23wpFaiRYkZFFKRETFZmHYEEc+nIKWNesiLSsTU87uxua/zkFwCKtc+Pr64ptvvkH37t2xYcMGuLu74/Dhw5gxYwYAICIiAhMnTkTz5s2xfv16uLq6YtWqVdi6daviHHfv3sXo0aOhra0NLy8vjBkzBv7+/vD09JSqW0REVExDWtpBQ6aG56mvcOzONanDKTQNqQMgIqLKrZFeLRzqNxGfn96B3x/dxjdhAbif8gKL2veDuhq/+ywrQgj4+vpiyJAhmD17NgCgU6dO0NfXx6xZsxAVFYUffvgB5ubmWL16NQCgW7duyMjIwKZNm+Dm5gYtLS34+PhALpfD29sbWlpacHBwQLVq1bBs2TJMmDAB9erVk7KbRERUBPV0aqC3kQWO3bkGvxuhGGhiK3VIhcK/FoiIqMRqaetgV68xcG1mBQDYeiMEk8/uxpuMynMLaGWTkpKCjz76CB9++KHS9ubNmwMAbt26hfDwcPTq1Utpf+/evfHy5UtEREQAAEJCQuDo6AgtLS1FGxcXF2RmZuLChQtl3AsiIiptY1p3AgBcfRqPK4/jJI6mcDhSSkREpaKahiY2dB+K+mE1sOX6BRy7cw1PXydha4+RqKWtI3V4VY6enh4WLFiQa/upU6cAAObm5khPT4exsbHSfiMjIwDA7du3YW1tjQcPHuRqY2BgAD09Pdy+fTvX+e3s7AqMKykpCXK5vEh9KS4B3iZORPQu+3rNYGnQEH8l3IffjRC0q9tU6pDeiyOlRERUatRkalhs/yEWte8HAPj90Z23zzJNTpQ2MBXxxx9/wMfHB87OzkhKSgLwtnjNSVdXFwCQnJycb5vsdjkXTKrQZFzqiIgom0wmw+fmnQEAx+5cw4NK8DxxjpQSEVGpG2/ZFfV0amDW+X24mfgYHx3zxo6eY9DaoL7UoVVZV65cwcSJE9G4cWMsW7ZMMcopy6dgU1NTUyxIlVcbIQTU8pgTHB4eXmAc7xtJJSKisvdRc2ssDz+Op2+SsT36Eua26y11SAXiSCkREZWJ/s2tsaPXGMg1tfHw1UsMDNjIZ5mWkYCAALi7u6NBgwbYtm0b9PX1FbfQvjvamf2zXC5XjJDmNSL66tWrcrsNl4iISpe2ugZGtOoAANgZ8zteV/A1HliUEhFRmencwAQH+05EPZ0aSEpPxYggPxzhs0xLlb+/P7744gvY2Nhg165dqFu3LgCgadOmUFdXR1yc8iIX2T8bGxtDV1cX9erVw927d5XaPHv2DMnJybnmmlY0nFFKRJS/kWYfQFNNHc9TX+FQbKTU4RSIRSkREZUpc4MGONJvMkxrvX2W6eSzu7Hl+nmpw6oS9u/fjxUrVqBPnz7w9fVVGtnU1taGnZ0dgoKClJ4bGxgYCLlcDktLSwBA586dERwcjLS0NKU26urqsLe3L7/OlABnlBIR5VZXR46PjP9bFf96SIV+hjiLUiIiKnON9GrhYN+JsK/XDACw9PIxfH35KLJElrSBVWLPnj3D8uXL0ahRIwwfPhw3btzA1atXFa+EhARMmjQJERERmDVrFs6ePQsvLy9s3boVEyZMQPXq1QEAY8eOxZMnTzB+/HgEBwfD398fnp6e+PTTT9GwYUOJe0lERCWRveBRTOKjCj2FhgsdERFRudDX1sFPvT7HtHN7cPzudfhcv4BHr5KwrutgaKvzclRU58+fx+vXr3Hv3j0MHz481/5Vq1ahf//+WL9+PX744QdMmTIF9erVw1dffYUxY8Yo2pmYmMDPzw+rVq3C9OnToa+vD3d3d0ybNq08u0NERGXAqnZjtK9rhLDHd7H1Rgg6N2whdUh5komKPI5bgWWvLvi+VQiJiEhZZlYWFl/+FduiLgIAOtVvDt8eI1FDq5rEkVFpKM/r44G/IzDz/D7U16mB8CHzy/z9iIgqo6O3/8TEMz9BBhnOD5qDZjUMpQ4pF96+S0RE5UpdTQ3fdPgI8+36AABCH8ZiUMCmSvEcNaqYZJxVSkSULxcjCzTUrQkBAf+oUKnDyROLUiIiKncymQyT2zjAq+un0JCpIer5Q/Q/5o2biY+kDo2IiKhK0VBTx6hWHQEAe2+FIyntjcQR5cailIiIJPNJi7bY3tMduhpauJ/yAh8f24jLj+5IHRYREVGVMszMHtXUNZGcnoo9t8KkDicXFqVERCSpbo1a4mDfCahbXY4XaW/wWaAvAu78JXVYREREVYa+tg6GtGwHAPC9HoL0rEyJI1LGopSIiCRnadgIh/tNQvMatZGamYEJwbvgf6NiznuhioRrNRIRFdZYiy6QQYZ7KYk4duea1OEoUbmidOvWrejXrx9cXV0xadIkJCQkSB0SEREBaCo3wOF+k9C2TlMICCz8/Qg8w09U6Id9U8Ug4zpHRETvZVyjNvoYWQAANl07V6GurypVlIaFheHnn3/Gvn378Ouvv6JFixZYu3at1GEREdF/DKrpYq/LWPRq0hoAsOHaGcw8vw9pmRnSBkZERFQFTLDsBgD4K+E+Qh/8I3E0/0+lilIDAwMsXrwYurq6AABzc3Pcv39f4qiIiCin6hpa8HEagRFmHQAAB/+JxOhT/0NyeqrEkREREVVu7eo2hX29ZgCAjX+dkzaYHKpcUfrrr7/Czs4u1ysxMREmJiawt7cHACQnJ8Pb2xs9e/aUOGIiInqXhpo6PDsOwJe2b39Hn7t/C58EbMbjV0kSR0YVSQW684yIqNKY+N9o6Zl7NxH9/KHE0bxV5YpSV1dXhIeH53rVqlVL0ebRo0cYOXIkbG1tMWzYMOmCJSKifMlkMsyw6YG1XT6BukwNfyXcR/9j3vjnxROpQ6MKRgZOKiUiKiznJq1gUrMOAGBzBRktrXJF6ftER0djyJAh6NGjB77++mupwyEiovcY0tIO/s6jUF1DE/HJzzHg2EZceRwndVhERESVkppMDeMtugIADsf+gQcpLySOSMWK0idPnsDd3R2zZ8/GlClTpA6HiIgKyamxGfa7jIdhNV08T32FISe2ICjuhtRhERERVUqDTGxRu5oe0rMy4VcBHsFWYYvSqKgoWFhY4OHD3Pc5Hz16FP369YOVlRX69OmDw4cPF+qcvr6+ePXqFXx9fdG/f3/0798fM2fOzLNtXvNSc76SkjiviYioPNnUaYJf+k2GkdwQbzLTMfa3HdgZ/bvUYREREVU61TQ04d66IwBgZ8wlvEx7I2k8GpK+ez5iY2MxYcIEZGTkfgTA8ePHMWfOHIwcORJdu3bFqVOnMHfuXFSrVg0uLi4FnnfevHmYN29eWYVNRERlrFkNQ/zSbxJGndqGP57+C4+Lh/Dg1QvMse0JGR9WqXIEuNIREVFxjWz1ATZcO4uk9FRsj76IqVaOksVSoYrSjIwM7N27F2vXroWmpmaebdatW4c+ffpg/vz5AICuXbvixYsX+P77799blBZFeHh4gfvt7OxK7b2IiKjwalfXwz6XcZgY/BOC78Xg+z9+w8NXL7Gi08fQVFOXOjySAL+PICIqOv1qunAz64DN189jy/UL+Ny8M6praEkSS4W6fffKlStYs2YNxowZgzlz5uTaHx8fj7i4OPTq1Utpe+/evREbG4v4+PjyCpWIiCSkq6kNP+eRGNLy7ReEe2+FY8yp7Ujhs0yJiIgKbZxlV2ipqePZmxT8dDNMsjgqVFFqYmKCU6dOYerUqVBXz/1td2xsLADA2NhYabuRkREA4Pbt22UfJBERVQiaaupY03kQZtr0AAAE34vBpye24OnrZIkjIyIiqhzq69RQfMG76do5pGXmnj5ZHipUUVq7dm0YGhrmuz97cSE9PT2l7bq6ugCA5GT+IUJEpEpkMhnm2PbEik4fQ00mwx9P/0X/Yxtx++VTqUOjcsAZpUREJTepjQPUZWp48OoFDv4TKUkMFaoofR8h3l5+3l3MInu7mlql6g4REZWSEWYd4Ovkhmrqmrib9AwDjm3E1Sec0qEqZOCkUiKi4moqN8CA5tYAgA1/nkFmVla5x1Cpqji5XA4g94hoSkqK0n4iIlI9vZqaY6/LWOhr6+DZmxQMPuGD0/HRUodFRERU4U21coQMMtxJeoajd66V+/tXqqI0ey5pXFyc0va7d+8q7SciItXUrq4RDvebhCZ6+nidkY4xp7djj4QLNxAREVUGLWvVhYuROQBg/Z/BijtRy0ulKkqNjIzQuHFjnDhxQml7UFAQmjVrhoYNG0oUGRERVRQmNevgl36TYWnQEJkiC3NCDsLr6ulyv8BS2eP/p0REpWfaf88pjX7+EKfio8r1vStVUQoAU6ZMwdGjR/H111/j3LlzWLJkCY4fP44ZM2ZIHRoREVUQdXXkONB3Aro1bAkAWBN5Eh6hh5CRlSlxZFQWOKeUiKjkrGo3hkMjUwDAD+U8WlrpitKBAwdi6dKluHDhAqZMmYLLly9j5cqV6Nu3r9ShERFRBaKnqY1tzqMwyMQWALDr5mWM+20nXmekSRwZERFRxTT9v9HSyCfxOHf/Vrm9b4UtSgcOHIiYmBjUr18/176hQ4ciKCgI165dQ0BAAAYMGFD+ARIRUYWnpa4Br66fYqpVdwDAyfgoDDnhi4Q3KdIGRkREVAF1qG+MTvWbAwDWRp4qt9HSCluUEhERlQaZTAaPdi74psNHkEGGiCdxGHBsI+KSEqQOjYiIqML5wtYZABDxJA5ny2m0lEUpERGpBHfzTtjsOBza6hqIffkU/Y9549rTe1KHVWqioqJgYWGBhw8fKm3v2bMnzMzMcr0SEv6/KL927Rrc3Nxga2uLLl26YN26dUhPTy/vLhQZlzkiIip9H9Rvjs4NTAAAayNPlstoqUaZvwMREVEF0beZJWpXHwv3U//Dk9fJ+OT4Zvg4jVAs7FBZxcbGYsKECcjIyFDanpKSgvj4eMyePRv29vZK+2rUqAHg7WPVRo8eDVtbW3h5eeGff/7Bd999h+TkZCxatKjc+lASMq5zRERUqr6wcUbIg38Q+SQeZ+7dhGNjszJ9P46UEhGRSrGv1wyH+k1EQ92aSMlIw6iT23Dg7wipwyqWjIwM7Nq1C5988glSU1Nz7Y+JiYEQAj169ICNjY3SS0Pj7ffSPj4+kMvl8Pb2hoODA8aMGYN58+Zhz549ePToUXl3iYiIKoAO9Y3RpUELAMC6q2U/t5RFKRERqRzTWvXwS7/JaK1fHxkiCzPP78PGa2elDqvIrly5gjVr1mDMmDGYM2dOrv1RUVHQ1tZGs2bN8j1HSEgIHB0doaWlpdjm4uKCzMxMXLhwoSzCJiKiSiB7bmnkk3gE37tZpu/F23eJiEglNdCtiYN9J2Ls6e0IfRiLb8NP4NOW7WBYTU/q0ArNxMQEp06dgqGhIX7++edc+2NiYlCrVi188cUXCAkJQWZmJrp374758+ejTp06eP36NR48eABjY2Ol4wwMDKCnp4fbt2/nOqednV2BMSUlJUEul5esY4Wkqaau9L9ERFR67Os1Q9eGLXD+/t/4JfYqnMrwFl4WpUREpLJqaFXDjl5jsC7yFFIyUmGgrSt1SEVSu3btAvdHR0fj6dOnaNmyJdzc3BAbG4sffvgBI0eOxKFDh5CUlAQA0NPLXYjr6uoiOTm5TOIuLT2atEL/5tZwaWohdShERFXStx0HwDP8BAa3aFem78OilIiIVJq2ugbm2blIHUaZWLBgAYQQsLa2BvB2lNPExATDhg3DkSNH4ODgAODtY3PeJYSAmlruWT7h4eEFvuf7RlJLk762DjY4fFZu70dEpGqMa9SGj9OIMn8fziklIiKqoqysrBQFabZ27dpBLpcjOjpaMUKa14joq1evyu02XCIiUm0sSomIiKqgV69e4eDBg4iOjlbaLoRAeno69PX1oauri3r16uHu3btKbZ49e4bk5ORcc02JiIjKAotSIiKiKkhbWxsrV67Ejz/+qLT99OnTePPmjeK5pZ07d0ZwcDDS0tIUbQIDA6Gurp7r2aZERERlgUUpERFRFaSuro5Jkybh5MmTWLZsGUJDQ7Ft2zbMnTsXPXr0QIcOHQAAY8eOxZMnTzB+/HgEBwfD398fnp6e+PTTT9GwYUOJe0FERKqACx0RERFVUe7u7tDT08P27duxf/9+1KxZE0OHDsW0adMUbUxMTODn54dVq1Zh+vTp0NfXh7u7u1IbIiKisiQTQgipg6iMslcXfN8qhERERKqE10ciIioq3r5LREREREREkmFRSkRERERERJJhUUpERERERESSYVFKREREREREkmFRSkRERERERJJhUUpERERERESS4XNKi0kul0sdAhERUYXD6yMRERUVn1NKREREREREkuHtu0RERERERCQZFqVEREREREQkGRalREREREREJBkWpURERERERCQZFqVEREREREQkGRalREREREREJBkWpURERERERCQZDakDqOocHR2RlJQkdRhERFQO5HI5goODpQ6jUuD1kYhItRR0jeRIaRWRlJTEi3sBmJ+CMT/vxxwVjPkhqfCzV3jMVeExV4XHXBUec5U/mRBCSB0ElZydnR0AIDw8XOJIKibmp2DMz/sxRwVjfkgq/OwVHnNVeMxV4TFXhcdc5Y8jpURERERERCQZFqVEREREREQkGRalREREREREJBkWpURERERERCQZFqVEREREREQkGRalREREREREJBkWpURERERERCQZPqeUiIiIiIiIJMORUiIiIiIiIpIMi1IiIiIiIiKSDItSIiIiIiIikgyLUiIiIiIiIpIMi9IK4OjRo+jXrx+srKzQp08fHD58uMD2KSkpWLp0KTp37gxbW1uMGzcOd+7cUWqTkZEBLy8vODg4wNraGsOGDcOff/6p1CYtLQ2bNm2Ci4sLbGxs4OLiAm9vb6SlpZVyD0tGyvx4eXnByckJ1tbWGDJkCC5cuFDKvSs5qfLzbvtBgwZh9OjRJe9QGZAqRw8fPoSZmVmu14cffljKPSwZKT9DJ0+exMcffwwrKys4Ojri+++/R0ZGRin2jqqyon52q6KsrCzs3r0brq6usLW1hbOzMzw9PZGcnKxoc+HCBQwaNAjW1tZwcnKCn59frvNcu3YNbm5usLW1RZcuXbBu3Tqkp6eXZ1fK1dSpU9GzZ0+lbcyTsrCwMHz22WewtrZGly5d8M033yAlJUWxn/n6f7t370afPn1gY2MDV1dXHDlyRGk/c1UIgiQVEBAgzMzMxPLly8W5c+fEokWLhKmpqTh+/Hi+x4wbN0588MEH4ueffxaBgYHC1dVVdO3aVbx8+VLRZsmSJcLa2lrs2LFDnD59WowYMULY2tqKuLg4RZuFCxcKa2trsXnzZhEaGio2b94srKysxMKFC8u0z0UhZX4WLFggbGxsxI4dO0RISIj48ssvhbm5ubhy5UqZ9rkopMxPThs2bBCmpqZi1KhRpd3FEpMyR8HBwcLU1FScP39eREZGKl5RUVFl2ueikDI/J06cEGZmZmLRokUiNDRU+Pj4CEtLS7F69eoy7TNVDcX57FZFmzdvFq1btxZr1qwRISEhYufOncLe3l6MGTNGCCHElStXhIWFhZgzZ444e/asWLdunTAzMxO+vr6Kc9y5c0e0bdtWfP755+LMmTNi69atwtLSUixdulSqbpWpw4cPC1NTU+Hs7KzYxjwpi4yMFBYWFmLatGkiJCRE/PTTT8LOzk7MnDlTCMF85bRnzx5hamoqVqxYIUJCQsS3334rTE1NRUBAgBCCuSosFqUSc3Z2VvwDzzZjxgzh4uKSZ/uwsDBhamoqzp49q9j27NkzYWNjIzZv3iyEECI+Pl60bt1a/PTTT4o2qamponv37mLRokVCCCESExOFmZmZ2LJli9L5N2/eLExNTUViYmKp9K+kpMrP06dPRevWrcXWrVsVbTIzM4WTk5OYM2dOqfWvpKTKT05RUVHCyspKdO7cuUIWpVLmaOPGjaJTp06l2Z1SJ1V+srKyhKOjY65/T6tXrxYDBw4slb5R1VbUz25VlJWVJdq3by+WLFmitP3YsWPC1NRU3LhxQ4waNUoMHjxYaf+qVauEnZ2dSE1NFUIIMX/+fOHg4KD4WQghdu3aJVq3bi0ePnxY9h0pRw8fPhTt27cX3bp1UypKmSdlw4cPF8OHDxdZWVmKbTt37hQ9evQQr169Yr5yGDJkiHBzc1PaNmzYMDFixAghBD9bhcXbdyUUHx+PuLg49OrVS2l77969ERsbi/j4+FzHhISEQFdXF507d1ZsMzAwQPv27XHu3DkAwKVLl5CZmYnevXsr2mhpaaF79+6KNsnJyRg6dCicnJyUzt+8eXNFbFKTMj81atTA/v37MWjQIEUbNTU1aGhoIDU1tVT7WVxS5idbeno65s6dCzc3NxgbG5dm90qF1DmKioqCmZlZaXer1EiZn+vXr+PevXsYPny40vnnzJmDgwcPllofqWoqzme3KkpJScFHH32Ua0pA9rX81q1bCA8PzzNPL1++REREBIC3/64dHR2hpaWlaOPi4oLMzMwKOW2lJBYsWIDOnTujY8eOim2pqanMUw4JCQkIDw/HZ599BplMptg+fPhwnDp1CmpqasxXDqmpqdDV1VXaVqtWLSQmJvKzVQQsSiUUGxsLALn+mDcyMgIA3L59O89jjIyMoK6urrS9adOmivaxsbGoWbMmDAwMcp33/v37ePPmDRo1aoQlS5YoLlzZTp8+DU1NTUUMUpIyP5qamrCwsEDNmjWRlZWFBw8eYPny5YiLi8PgwYNLrY8lIWV+sv34449IT0/H9OnTS96hMiB1jqKjo/HmzRt89tlnaNOmDTp16oS1a9dWmDkiUuYnJiYGAKChoYHPP/8clpaW+OCDD+Dl5YWsrKzS6SBVWcX57FZFenp6WLBgAdq1a6e0/dSpUwAAc3NzpKenF5in169f48GDB7naGBgYQE9Pr0rlcv/+/bh+/ToWLlyotD0+Pp55yuHmzZsQQqBmzZqYOXMmbGxs0K5dOyxevBhv3rxhvt4xcuRInD9/HsePH0dycjJOnDiBM2fOoH///sxVEWhIHYAqS0pKAvD2opJT9rctORcpyJacnJyrffYx2e0LagO8/Wa1WrVqufafPHkShw4dwsiRIyGXy4vYm9JXUfKzbds2rFy5EgAwZMgQdOrUqTjdKXVS5+fPP/+En58fdu3apfTNXkUiZY6EEIiLi8OLFy/w5ZdfYtasWbh06RJ8fHzw+PFjxWdKSlLmJyEhAQAwZcoUDBgwAOPGjcPly5exadMmqKurY9q0aSXoGVV1xfnsqoo//vgDPj4+cHZ2LlSe8muT3a6q5PLevXvw9PSEp6dnri/MmCdl2b+fPTw80LNnT2zcuBExMTHw8vJCamoqhgwZAoD5ytavXz9cunQJM2fOVGz7+OOPMXbsWERGRgJgrgqDRamEhBAAoHRrRM7tamq5B7Kz9+Ulu31+bfJ7PwAICgrC7Nmz0a5dO8yePbsQ0Ze9ipKfbt26wdLSEhEREdiwYQPS09Ph6elZyF6UHSnzk5qaCg8PD4waNQpWVlZFD76cSJkjdXV1+Pn5oVGjRmjatCkAwN7eHpqamvDy8sKkSZPQrFmzonWolEmZn+zR4r59+2LWrFkAgA8++AAvX77Eli1bMG7cuDy/PCMCivfZVQVXrlzBxIkT0bhxYyxbtkwxwpLXdR94m6eC/jYQQlSJXAohMH/+fDg4OChNK8i5H2CesmX/fm7bti0WL14MAOjYsSOEEFi5ciU+/fRTAMxXtkmTJiEyMhLz5s2Dubk5/vjjD3h7e0NPTw99+/YFwFwVhmr0soLKHo189xuQ7OW28xqt1NPTU1qOO+cx2d+wFNQme39O27Ztw4wZM9C2bVts3rwZ2traxehN6aso+WnRogXs7e0xceJETJgwAYcOHcKjR4+K0aPSJWV+sm+xnDx5MjIyMpCRkQHxduE0xX9XBFLmSEtLCx07dlQUpNm6d+8O4O2tvVKTMj/Z3xJ369ZNqU2XLl2Qmpqa6xEzRDkV57Nb1QUEBMDd3R0NGjTAtm3boK+vn2+esn+Wy+WKf7d5jca8evWqSuRy165diImJwfz585WuWcDbx1cxT8oK+v0shMC1a9cAMF8AEBERgQsXLmDBggUYPXo07O3tMW7cOHh4eGDHjh3Q0dEBwFwVBotSCWXfOx4XF6e0/e7du0r73z0mPj4+1x/9d+/eVbRv3rw5EhMT8eLFi1xtGjdurHSr5bfffgtPT0/07dsXW7ZsyfPWAalImZ9Hjx7h4MGDuf6wtrCwgBACjx8/LlnnSoGU+QkMDMTt27dha2sLCwsLWFhYICwsDJcuXYKFhQUuX75cav0sCSlzFB8fj7179ypug8qWPd9UX1+/BD0rHVLmJ3uU+N3nIleU+bZUsRXns1uV+fv744svvoCNjQ127dqFunXrAng711tdXT1XnrJ/NjY2hq6uLurVq6fIXbZnz54hOTm5SuQyMDAQz58/R5cuXRTXrMOHDyMuLg4WFhYIDw9nnnJ43+/nxo0bM1//uX//PoC3o8o52dnZAXi74CFzVTgsSiVkZGSExo0b48SJE0rbg4KC0KxZMzRs2DDXMV26dMHLly8RGhqq2Ja9Slr2XMfs/w0MDFS0SUtLw9mzZ5XmQ3p5eeF///sf3N3dsWbNmgo3L1DK/Dx58gTz58/P9d4hISHQ0tKqEL8gpMzPxo0bceDAAaWXhYUFrKysFP9dEUiZo5cvX2LRokU4evSo0vkDAgKgp6cHc3Pz0ulkCUiZHzs7O1SvXh0BAQFK5w8ODkatWrVgYmJSOp2kKqk4n92qav/+/VixYgX69OkDX19fpVEVbW1t2NnZISgoSOmLpMDAQMjlclhaWgIAOnfujODgYKUiJDAwEOrq6rC3ty+/zpSRpUuX5rpmOTo6on79+jhw4ABcXFyYpxxMTEzQqFGjPH8/a2howNbWlvn6T/bfg2FhYUrbr169CuDtl7TMVSGV4uNlqBgOHjwoTE1NxdKlS8XZs2fF4sWLhampqTh27JgQ4u3z/yIjI0VSUpLimBEjRgh7e3uxb98+ERQUpHhwfc5ni86dO1e0adNG+Pv7i99++024ubkJW1tbcefOHSGEENHR0aJVq1Zi4MCBIjIyMtcr5/tJSar8ZGVliQkTJgg7Ozuxc+dOERISIpYvXy5atWolvL29yzcJBZAqP3kZMWJEhXxOqZQ5mjhxorCxsRH+/v7iwoULYvny5cLMzEz4+/uXW//fR8r8+Pj4CFNTU7F8+XIRGhoqVq1aJczMzMT27dvLLwFUab3vs6sKnj59KqytrYWjo6MICwvLdS1/9uyZCA0NFWZmZmLGjBnizJkz4rvvvhNmZmbCx8dHcZ6///5btGnTRowaNUr89ttvws/PT1haWorFixdL17kyNnfuXKXnlDJPyo4dOybMzMzE7NmzRUhIiNi8ebOwsLAQnp6eQgjmK6fJkycrrvWXLl0Svr6+wtbWVowdO1YIwVwVFovSCmD37t2iZ8+ewtLSUvTp00ccOnRIsS/7onvp0iXFtsTEROHh4SHs7OxE27Ztxbhx48Q///yjdM7U1FSxfPly0bFjR2FtbS2GDRsmrl69qti/fv16YWpqmu8rLCyszPtdWFLkRwghUlJSxOrVq4Wjo6OwsLAQ/fr1E/v27SvTvhaHVPl5V0UtSoWQ/jPk5OQkLC0tRd++ffkZesfevXuFi4uLsLCwEL169RK7d+8us35S1VPQZ1cVHDp0qMBr+eHDh4UQQgQFBYkPP/xQWFhYCCcnJ7F169Zc5woLCxODBw8WlpaWomvXrmLt2rUiLS2tvLtUbt4tSoVgnt518uRJMWDAAGFpaSkcHByEt7e3yMzMVOxnvt5KTU0Va9euFQ4ODsLS0lL07t1brF+/XqSmpiraMFfvJxOigqxIQkRERERERCqHc0qJiIiIiIhIMixKiYiIiIiISDIsSomIiIiIiEgyLEqJiIiIiIhIMixKiYiIiIiISDIsSomIiIiIiEgyLEqJiIiIiIhIMhpSB0BERERERNJYsmQJgoODkZSUBF1dXbi4uODLL7+ElpaW1KGRCuFIKZGK2bVrFz755BO0adMGbm5uUodTZJU9fiIioopkxIgROH78OCIiIvDLL78gOjoaW7ZskTosUjEsSolUTJ06dTB+/HiMHj1a6lCKpbLHT0RElcOePXvg7Oyc7343Nzc4OTmVY0Rlo0WLFtDR0VH8LJPJEBcXJ2FEpIpYlBKpmF69eqFXr14wNDQs9jned6H28PCAmZkZ/v3331z7EhIS4OLiglatWuHQoUNFfu/SiJ+IiOh9QkND0alTJ6nDKBc+Pj6wtbVFx44dERMTgxEjRkgdEqkYFqVEVGTFvVAnJydj7NixuH37NhYuXIiPP/64DKIjIiIqmaysLPz+++/o2LGj1KEUS1paGlJSUvJ9ZWZmKrUfP348IiMjERAQgCFDhqBOnToSRU6qigsdEVGRZF+olyxZUqTj0tLSMGnSJFy/fh1ffvklhg8frrQvPT0932OrVasGdXX14oZMRERUJDdu3MCLFy/wwQcfSB1KscyfPx+//vprvvu3b9+ODh065NpuYmKCVq1awcPDA9u2bSvDCImUsSgloiIpzoU6MzMTM2fOxOXLlzFlyhSMHTtWaX9xL55EREQFSUhIwKpVqxAcHIyXL1/C0NAQQ4cOxZQpUyCTyfI9LjQ0FObm5tDX1y/0e6WkpGD06NH4+++/4evri3bt2gEAzp49ix9//BG3bt2CoaEh3N3dERUVhYsXL+K3337L93xOTk7o3r07WrduDV9fXzx48AAtW7bE4sWL0aBBAyxbtgznzp2Dnp4eBg4ciBkzZkBN7e1NkGvWrMGaNWsKHXtOWVlZuHv3brGOJSouFqVEVCTFuVAvWLAAp0+fhru7O6ZPn55rf0kunkRERHm5efMm3N3dIZfLMXr0aGhra2Pfvn1Yv349atWqVeC8ydDQ0CLdupuWloapU6ciJiYGmzdvVhSkwcHBmDJlCkxNTTFr1iw8evQIK1euhI6ODnR1dd973tOnTyMoKAijRo2CEAIbN27EtGnTIJfL0bJlS3h4eCAoKAibNm1Cs2bNijwtJiUlBSdOnEDPnj0hl8tx8+ZNeHt7o0uXLkU6D1FJsSglUjEZGRnIzMxEZmYmsrKykJqaCplMVujnkRX1Qr1y5Ur8/PPPcHZ2hoeHR3HDVihp/EREVPWlpqZi8uTJMDY2xpYtW1C9enUAwIcffghHR0ccOHAg36I0NTUVERERGDduXKHeKysrC3PmzEFYWBg2bNigdI389ttv0aRJE+zZswfVqlUDALRt2xZTpkwpVFH66NEj/PLLLzAzMwMAJCYmYuvWrWjbti2+++47AICrqyvs7e1x4cKFIhelMpkMR48excqVK5Geng4DAwP06tUrzy+QicoSi1IiFbNx40b8+OOPip+trKxgb2+PHTt2vPfYol6ofXx8sHfvXshkMkRGRiIhIQEGBgbFjh0oWfxERKQaduzYgXv37sHPz09RkAJA3bp10aJFCzx8+DDfY8PDwyGEUIx2vs/ixYsRGBiIFStWwMHBQbE9OjoacXFx8PDwUBSkAODs7AwTExO8efPmvedu2rSpoiAFAGNjYwBAz549Fdt0dHRgaGiIJ0+eFCrenHR0dODv71/k44hKG4tSIhUzbdo0TJs2rVjHFvVCvXfvXri4uKBLly5YsGABlixZgh9++KFY752tJPETEZFqOHHiBGxsbNC0adNc+968eYOaNWvme2xoaCjatm2rVEjm5969e9i/fz8AICIiQmmkMnteppGRUa7jjI2NERUV9d7zv/v4s+xF/979glddXR1CiPeej6ii4iNhiKjQinKhBoCuXbti9erVGDx4MDp16oTAwMACFzQiIiIqqRcvXuCvv/6CpaVlrn2JiYn4999/0b59+3yPv3jxYqEfeyaTybBkyRL0798f+/fvx5UrVxT7MjIyACDP6SXa2tqFOr+GRt7jRwUt0kRUGbEoJaJCK8qFGnh7S1P2xfjrr7+Gjo4Oli1bhsePH5dViEREpOJu3LgBIQR0dHRy7Tt69CgyMjLQt2/fPI99/vw5oqKiCr12QsOGDTF06FB89dVX0NPTw+LFixWPOGvSpAkA4M6dO7mOy2sbkSpjUUqkQszMzAr1yktRL9SA8je5TZo0wYwZM5CYmIiFCxeWuC9ERER5yb4t9tKlS0q3tMbGxsLLywsdO3bM97FmFy9ehFwuz3OUtSC1a9fG9OnTcevWLfj5+QEALC0t0aBBAxw4cABpaWmKtlevXsWNGzeK2i2iKo1zSolUSExMTLGPLe6FOqeRI0ciICAAZ86cwYEDB/DJJ58U+1xERER5uX79OgwNDXHv3j1MmjQJXbt2xb///os9e/agVq1aWLVqleJ5nu+6ePEiOnTokO/+ggwfPhwHDx6Et7c3+vbtiyZNmsDDwwMzZ87E0KFD0b9/fyQkJGD79u1cMZ7oHRwpJVIhAQEBsLW1ha2tLaytrWFmZqb42dbWVrFYQ15KcqHOpqamhuXLl0NTUxOenp64f/9+sc9FRESUl6ioKLRu3Rre3t5ISEjAypUrceTIEbi6uuLAgQOoW7duvscW9bFnOamrq2Px4sVITU3F0qVLAQAuLi747rvvkJGRgdWrV+Po0aOYN28eLC0tWZgS5SATXKqLSCWdOHECixcvxu+//16o9j169MDnn3+OYcOGvbeth4cHDh06hNOnT6Nx48a59v/4449Yv349OnXqBD8/Py7YQEREpeL169do164dRo8eja+++krSWDIzM/HixYs8H4Xm6uqKGjVqYNeuXRJERlTxcKSUSEVFRUXBwsKi0O1Pnz5dqIIUAFasWIGYmJg8C1IAmDp1KmJiYuDv78+ClIiISk10dDQyMzNhamoqdSjIzMxEt27dsGjRIqXtN2/exK1bt2BlZSVRZEQVD+eUEqmo69evo3Xr1lKHQUREVGqyFxBq1aqVxJG8fRSMi4sLDhw4AJlMBktLSzx+/Bi7d++Gvr4+3N3dpQ6RqMJgUUqkoqKiojBw4ECpwyAiIio1UVFR0NDQQPPmzaUOBQCwbNkyGBsb48iRIzh06BDkcjk6duyImTNnFji3lUjVcE4pkQp69OgRunXrhsDAQDRr1kzqcIiIiIhIhXFOKZEKunHjBnR1dWFkZCR1KERERESk4liUEqmg7PmkXGSIiIiIiKTGopRIBUVFRcHc3FzqMIiIiIiIOKeUiIiIiIiIpMORUiIiIiIiIpIMi1IiIiIiIiKSDItSIiIiIiIikgyLUiIiIiIiIpIMi1IiIiIiIiKSDItSIiIiIiIikgyLUiIiIiIiIpIMi1IiIiIiIiKSzP8BVhErJ5ug6PoAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1, 2, figsize=(15, 6))\n", "ax[0].set_title(\n", " \"saturation pressure of {}\".format(\n", " parameters.pure_records[0].identifier.name\n", " )\n", ")\n", "sns.lineplot(\n", " y=phase_diagram.vapor.pressure / PASCAL, \n", " x=1.0/phase_diagram.vapor.temperature * KELVIN, \n", " ax=ax[0]\n", ")\n", "ax[0].set_yscale('log')\n", "ax[0].set_xlabel(r'$\\frac{1}{T}$ / K$^{-1}$');\n", "ax[0].set_ylabel(r'$p$ / Pa');\n", "#ax[0].set_xlim()\n", "#ax[0].set_ylim()\n", "\n", "ax[1].set_title(\n", " r\"$T$-$\\rho$-diagram of {}\".format(\n", " parameters.pure_records[0].identifier.name\n", " )\n", ")\n", "sns.lineplot(\n", " y=phase_diagram.vapor.temperature / KELVIN, \n", " x=phase_diagram.vapor.mass_density / KILOGRAM * METER**3, \n", " ax=ax[1], \n", " color=colors[0]\n", ")\n", "sns.lineplot(\n", " y=phase_diagram.liquid.temperature / KELVIN, \n", " x=phase_diagram.liquid.mass_density / KILOGRAM * METER**3, \n", " ax=ax[1], \n", " color=colors[0]\n", ")\n", "ax[1].set_ylabel(r'$T$ / K');\n", "ax[1].set_xlabel(r'$\\rho$ / kg m$^{-3}$');\n", "#ax[1].set_ylim()\n", "#ax[1].set_xlim()\n", "\n", "sns.despine(offset=10)" ] } ], "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.9.7" } }, "nbformat": 4, "nbformat_minor": 5 }