{ "cells": [ { "cell_type": "markdown", "id": "06b263a1-b920-4b97-9bdf-238a87893e92", "metadata": {}, "source": [ "# Surface tension calculations using DFT for Hydrogen\n", "Comparing Feynman-Hibbs corrections to first and second order" ] }, { "cell_type": "code", "execution_count": 17, "id": "7f9082f3-e3e2-421b-b09a-a39160a83183", "metadata": {}, "outputs": [], "source": [ "from feos import *\n", "from feos.si import *\n", "from feos.dft import *\n", "from feos.saftvrqmie import *\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import seaborn as sns\n", "import json\n", "\n", "sns.set_context('talk')\n", "sns.set_palette('Dark2')\n", "sns.set_style('ticks')" ] }, { "cell_type": "markdown", "id": "0c514e92-9e5b-45a5-98b1-81f83eb4c104", "metadata": {}, "source": [ "# Pure fluid surface tension correlations" ] }, { "cell_type": "code", "execution_count": 18, "id": "b334936e-587b-44d0-a92e-6dacd5f08e73", "metadata": {}, "outputs": [], "source": [ "def surftens_mulero2012(fluid, tr):\n", " \"\"\"\n", " Calculate pure fluid surface tension using Mulero 2012 correlation (doi:10.1063/1.4768782)\n", " Args:\n", " fluid (str): Component name\n", " tr (np.ndarray): Reduced temperature\n", " Returns:\n", " sigma (np.ndarray): Surface tension (mN/m)\n", " \"\"\"\n", " ff = open(\"mulero_2012_parameters.json\", \"r\")\n", " complist = json.load(ff)\n", " ff.close()\n", " sigma = np.zeros_like(tr)\n", " for i in range(len(complist[fluid][\"sigma\"])):\n", " sigma[:] += complist[fluid][\"sigma\"][i] * \\\n", " (1-tr[:])**complist[fluid][\"n\"][i]\n", " return sigma * NEWTON / METER / (MILLI * NEWTON/ METER)" ] }, { "cell_type": "markdown", "id": "f445dc40-3c4f-417d-a430-ed4b11d9216f", "metadata": {}, "source": [ "# Hydrogen parameters and functional" ] }, { "cell_type": "code", "execution_count": 19, "id": "4c8f318f-56e0-4a48-b21b-1ecab6bc10a3", "metadata": {}, "outputs": [], "source": [ "# Critical and normal boiling point temperature\n", "Tc = 33.145\n", "Tnb = 20.369\n", "# First order Feynman-Hibbs corrections\n", "parameters_fh1 = SaftVRQMieParameters.from_json([\"hydrogen\"], \"../../parameters/saftvrqmie/aasen2019.json\")\n", "func_fh1 = HelmholtzEnergyFunctional.saftvrqmie(parameters_fh1)\n", "state_fh1 = State(func_fh1, temperature=Tnb * KELVIN, pressure=20.0 * BAR)\n", "model_fh1_tc = State.critical_point(func_fh1).temperature\n", "# Second order Feynman-Hibbs corrections\n", "parameters_fh2 = SaftVRQMieParameters.from_json([\"hydrogen\"], \"../../parameters/saftvrqmie/aasen2019_fh2.json\")\n", "func_fh2 = HelmholtzEnergyFunctional.saftvrqmie(parameters_fh2)\n", "state_fh2 = State(func_fh2, temperature=Tnb * KELVIN, pressure=20.0 * BAR)\n", "model_fh2_tc = State.critical_point(func_fh2).temperature" ] }, { "cell_type": "markdown", "id": "2ee15889-18da-4802-842f-5e042d673ecb", "metadata": {}, "source": [ "# Calculate and plot surface tension" ] }, { "cell_type": "code", "execution_count": 20, "id": "c284f4b3-af36-4d50-8391-a7e44148b050", "metadata": {}, "outputs": [], "source": [ "# First order Feynman-Hibbs corrections\n", "dia_fh1 = PhaseDiagram.pure(func_fh1, Tnb * KELVIN, npoints=50)\n", "sft_dia_fh1 = SurfaceTensionDiagram(dia_fh1.states, n_grid=1024, l_grid=200 * ANGSTROM , critical_temperature=model_fh1_tc)\n", "surf_tens_fh1 = sft_dia_fh1.surface_tension / (MILLI* NEWTON / METER)\n", "# Second order Feynman-Hibbs corrections\n", "dia_fh2 = PhaseDiagram.pure(func_fh2, Tnb * KELVIN, npoints=50)\n", "sft_dia_fh2 = SurfaceTensionDiagram(dia_fh2.states, n_grid=1024, l_grid=200 * ANGSTROM , critical_temperature=model_fh2_tc)\n", "surf_tens_fh2 = sft_dia_fh2.surface_tension / (MILLI* NEWTON / METER)" ] }, { "cell_type": "code", "execution_count": 21, "id": "4e69ad65-f86a-432a-8c1b-b6a2378ff9ec", "metadata": {}, "outputs": [], "source": [ "tr = np.linspace(Tnb/Tc, 1.0, 100)\n", "s_corr = surftens_mulero2012(\"hydrogen\", tr)" ] }, { "cell_type": "code", "execution_count": 22, "id": "0247cea8-f6cd-436f-9e98-7dc14130e5ec", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlsAAAHNCAYAAAApEr6yAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACv6ElEQVR4nOzddXhVx9bA4d/JiXsCIcQIwQIEt+ItVqBIKRJcixW3UooVt+JuxV2LtbhrKe6aIAEixD1H7h8pJ+QmQCBGwnqfJ893s/beM7P7kZyVvWfWKLRarRYhhBBCCJEu9DJ7AEIIIYQQ2ZkkW0IIIYQQ6UiSLSGEEEKIdCTJlhBCCCFEOpJkSwghhBAiHUmyJYQQQgiRjiTZEkIIIYRIR/qZPQABNWrUIDAwECMjI5ydnTN7OEIIIYRIgefPnxMTE4OtrS3Hjh1753mSbH0GAgMDiY6OJjo6mpCQkMwejhBCCCE+QmBg4HuPS7L1GTAyMiI6OhpjY2Py58+f2cMRQgghRAo8evSI6OhojIyM3nueJFufAWdnZ0JCQsifPz87duzI7OEIIYQQIgWaNm3KrVu3PjgFSCbICyGEEEKkI0m2hBBCCCHSkSRbQgghhBDpSJItIYQQQoh0JMmWEEIIIUQ6kmTrC7XzwU4mnJ+Ad4h3Zg9FCCGEyNak9MMX6GHQQ0afHQ3Afu/9LK2zlKI5imbyqIQQQojsSZ5sfYHMDc0x0DMAICQmhK4HunLd/3omj0oIIYTIniTZ+gLlNsvNlGpT0FfEP9gMiwuj28FuXPK9lMkjE0IIIbIfSba+UN/m/ZaZ38zUPeGKVEXy0+GfuPDyQiaPTAghhMheJNn6gtXIU4O5NedipIzf0ylKFUXvI7057XM6k0cmhBBCZB+SbH3hqjpVZX6t+ZjomwAQo46h39F+HHt6LJNHJoQQQmQPkmwJKjpUZFHtRZjqmwIQp4lj0PFBHPQ+mMkjE0IIIbI+SbYEAGXty7L026VYGFgAoNKq+Pnkz+x9vDeTRyaEEJ9m3rx5uLu7Z/YwPivDhg2jZs2aqWpjzJgxdO7cOY1GlLE8PT2ZNm1ahvcryZbQKWlXkmV1l2FlZAWARqth+Knh7HywM5NHJoTIbnbs2IG7uzvu7u78+++/SY5rtVq+/vpr3N3d6dGjRyaMMHOdOHGCefPmZfYwknj27Bnbtm1L9P+Tly9fMn/+fJo3b0758uX56quvaN++PWfPnk22jdDQUEaNGkXFihUpVaoU7du359atW0nO++uvvxgyZAjffvst7u7utG/fPtn2rl+/zrhx42jQoAGlSpXim2++oX///nh5eSU5t1u3bmzYsAF/f/9P/C/waSTZEol45PDgj2//wNbYFgAtWkafHc2mu5syeWRCiOzIyMiIvXuTPkH/559/ePXqFYaGhpkwqsx34sQJ5s+fn9nDSGLNmjU4OTlRsWJFXezIkSMsW7YMV1dXBgwYQK9evYiIiKBz585s37490fUajYbu3buzd+9e2rVrx88//0xgYCDt27fH29s70bkbN27kyJEj5M6dGysrq3eOafny5Rw8eJBKlSoxYsQIPD09+ffff2natCn3799PdG6tWrUwNzdnw4YNqf+P8REk2RJJuNu6s7LuSuxM7HSxiRcmsvrW6kwclRAiO/r666/Zv38/KpUqUXzv3r14eHhgZ2f3jiszXlRUVGYPIVPFxcWxZ88e6tevnyj+1VdfcezYMWbMmEHbtm3p2LEjmzZtIl++fMydOzfRufv37+fKlStMnjyZPn360LZtW9auXYtSqUzyJG/atGlcunSJNWvWkCtXrneOq1OnThw9epSRI0fSokULevXqxfr161GpVCxdujTRuXp6etStW5ddu3ah1WpT+V8k5STZEsnKZ52PlfVWktssty42/d/pLL2+9D1XCSHEx2nQoAHBwcGcOXNGF4uNjeXAgQM0atQoyfkXLlzA3d2dCxcS1wR8/vw57u7u7Nix44N97tq1i6ZNm1KiRAkqVKjAwIEDefnyZaJz2rdvT8OGDbl58yZt27alZMmSzJw5E4DXr18zfPhwKleuTPHixWncuDE7d6Z8usWJEydo06YNpUqVonTp0nTv3p0HDx7ojg8bNoz169cD6F61fmju2eHDh+nevTtVq1alWLFi1K5dmwULFqBWq1M8rg+5dOkSQUFBVK5cOVG8YMGC2NraJooZGhry9ddf8+rVK8LDw3XxAwcOkDNnTr799ltdzNbWlvr163PkyBFiY2N1cQcHB/T0PpymlClTJskT0Lx581KwYEEeP36c5PzKlSvj4+PDnTt3Pth2WpFkS7yTq6Urq+qtwsncSRebd2Uecy/PzdC/CIQQ2ZeTkxOlSpVi3759utjJkycJCwvju+++S/P+Fi1axC+//IKrqyvDhg2jQ4cOnDt3jrZt2xIaGpro3ODgYLp160aRIkUYPnw4X331FdHR0bRv357du3fTqFEjhg4dioWFBcOGDWP16g8//f/zzz/p0aMHpqamDBkyhF69evHw4UPatGnD8+fPAWjZsiVVqlQB4p/uvPl6n507d2Jqakrnzp0ZMWIEHh4ezJ07l+nTp3/if6mkrly5gkKhoGjRlO2l6+/vj4mJCSYmJrrYnTt3KFq0aJIkqnjx4kRFRSU7z+pTaLVaAgICsLGxSXKsWLFiAFy+fDlN+koJ2YhavJeTuROr6q2i28FueId6A7DsxjJi1DEMKTcEhUKRuQMU4gumiYlBGxOT4vP1LCyS/Myqw8IghX88KQwM0HvrgzOtNGrUiBkzZhAdHY2xsTF79uyhfPny2Nvbp2k/Pj4+zJs3jwEDBtCzZ09d/Ntvv+WHH35gw4YNieL+/v6MHTuWVq1a6WKrV6/m0aNH/P777zRu3BiAVq1a0b59e2bPnk2zZs0wNzdPtv+IiAgmTpxIixYtGD9+vC7+ww8/UK9ePZYsWcL48eMpXbo0efPm5cyZM3z//fcpurcZM2ZgbGys+75169aMHj2ajRs3MnDgwDSZ+/b48WOsrKzeeX9ve/LkCYcOHaJevXoolUpd3N/fn3LlyiU5/81rQj8/vzRZQbp79258fX3p169fkmP29vYYGBjw8OHDVPeTUpJsiQ/KbZablfVW0u1gNx4Gx//jXHN7DTHqGIZ/NRw9hTwgFSIzvF66jIAFC1J8fqF/LqC0tEwUe1izFpqwsBRdb9WkCY5TJn/UGFOifv36TJo0iWPHjlGtWjWOHz/OyJEj07yfQ4cOodFoqF+/PoGBgbp4zpw5cXV15cKFC4mSLUNDQ5o2bZqojZMnT2JnZ0fDhg11MQMDA9q3b8+gQYO4ePEiNWrUSLb/s2fPEhoaSoMGDRL1r6enR8mSJZO8Gv0Ybyda4eHhxMbGUq5cOTZv3szjx48pXLjwJ7f9RnBw8Hsnqr8RFRVF//79MTY2ZvDgwYmORUdHJ5v4vYnFfMQfD+/y6NEjxo0bR+nSpfnhhx+SPcfKyoqgoKBU95VSkmyJFMlpkpMVdVfQ41AP7gTGv+fefG8zsepYfqv0G0o95QdaEEKI5Nna2lKpUiX27t1LdHQ0arWaunXrpnk/3t7eaLXaRPOF3qavn/gj0d7ePkli4OPjg6ura5LXYPnz5wfgxYsX7+0foGPHjskeT8kTo3d58OABs2fP5vz584nmSAGEpTCZTokPTSFRq9UMHDiQhw8fsmzZsiRPJ42NjRPNy3rjTczIyChV4/P396dHjx5YWFgwZ86cRE/V3qbVajP0zUyWTLa0Wi1Xrlzh6NGjXLp0icePHxMeHo6FhQVFixalSZMmNGrU6JP/Qx44cIB169Zx9+5d4uLicHV1pXHjxnTo0AEDA4M0vpusw8bYhmXfLqPX4V5cD7gOwM6HO4lRxzCx6kT09bLkPychxGegYcOGjBo1ioCAAKpXr47l/z2Be+Ndv9c1Gs0H+9BoNCgUCpYtW5bsh7CpqWmi799+WpQW3iQq06ZNS3aV5bsSgw8JDQ2lXbt2mJub069fP/LkyYORkRG3bt1i+vTpKfpvkxLW1tZJ5rX9v5EjR3L8+HGmT59OpUqVkhy3s7NLtsaVn58fwHtXHX5IWFgY3bp1IywsjPXr17/3NXRoaGiy87nSS5b8dDx//jydOnXSfe/i4oKTkxM+Pj6cOXOGM2fOsG/fPubNm/fR76mnTp3KihUrAMiTJw8mJiY8ePCAadOmcezYMVasWPHF1n0BsDKyYkmdJfQ+0pvLfvGTC//y+otYdSzTqk/DQPnlJqNCZLQc3bth27FDis/Xs7BIEitw9MhHzdlKL3Xq1OG3337j6tWrzJo1653nvUnC/v9pjY+Pzwf7yJMnD1qtFmdnZ9zc3D5pnE5OTty7dw+NRpPo6dabVW+Ojo7vvNbFxQWAHDlyJFnR9/8+5mHBP//8Q3BwMPPnz6d8+fK6+JsJ92klX7587Nmzh7CwMCyS+bc0depUduzYwfDhwxO9Zn1b4cKFuXTpUpL/ftevX8fExOST//8SExNDz5498fb2ZuXKlRQoUOCd5/r6+hIXF6d7GpkRsuRkmzc/LCNGjODs2bMcPnyYHTt2cOHCBaZOnYqhoSHHjx9nzpw5H9XuoUOHdMnUwoULOXToELt372bPnj04Oztz8eJF3dLfL5m5oTmLai/iK4evdLHDTw/T/1h/olXRmTgyIb4sekZGKC0tU/yV3Ae40sIixdenx+T4N8zMzBgzZgx9+/Z973YyTk5OKJVKLl68mCi+cePGD/bx7bffolQqmT9/fpLXYVqtNkVzeKpXr46/vz9//fWXLqZSqVi7di2mpqaJkp3/V61aNczNzVmyZAlxcXFJjr89j+vNCr4PPUkCdEnL2/cUGxub4sKdT58+5enTpx88r1SpUmi1Wm7evJnk2PLly1mxYgU9e/Z852tSgHr16hEQEMDBgwl77wYGBrJ//35q1KjxSQ8z1Go1AwYM4OrVq8yZM4fSpUu/9/w34//QeWkpSz7ZKlGiBPv370/2lV6TJk149eoVs2bNYtu2bQwePDhFdToAXbXebt26UatWLV08f/78TJgwgU6dOrF+/Xq6d++epKbIl8bUwJT5Necz6PggTvmcAuCUzyn6HOnD3JpzMTUw/UALQgiR2LsmM7/NwsKCevXqsW7dOhQKBS4uLhw/fpzXr19/8No8efIwYMAAZsyYgY+PD7Vr18bMzIznz59z+PBhPD09+fHHH9/bRsuWLdm8eTPDhg3j1q1bODk5ceDAAS5fvszw4cPfO+/K3NycMWPGMHToUJo2bcp3332Hra0tL1684MSJE5QpU4bRo0cD4OHhAcCECROoWrUqSqWSBg0aJNtu6dKlsbKyYtiwYbRv3x6FQvFRRTvfvCk6evToe88rW7Ys1tbWnDt3LtErwkOHDvH777+TN29e8uXLx65duxJdV6VKFXLmzAlA3bp1KVWqFL/++isPHz7ExsaGjRs3olar6du3b6LrLl68qEuqAwMDiYyMZOHChQCUL19el9hOmTKFo0ePUqNGDYKDg5P0//8rOs+ePYujo2OKS1ikhSyZbH1oEmH16tWZNWsWwcHBBAYG6v6f/D7e3t7cvXsXiP9h+n+VKlXC1dWVJ0+ecOTIEVq0aPFpg89GjPWNmVNjDkNPDuXw08MAXHh1gZ6He7Kg1gIsDJM+ZhZCiNQaOXIkKpWKTZs2YWhoSL169Rg6dOg7X129rXv37uTNm5dVq1ax4L+VnLlz56ZKlSop2qDZ2NiYtWvXMn36dHbu3El4eDhubm5Mnjw5ycrF5DRq1IhcuXKxdOlS/vjjD2JjY7G3t6dcuXKJrv/2229p3749+/btY/fu3Wi12ncmWzY2NixevJipU6cye/ZsLC0tady4MZUqVfpg8vgxDA0NadSoEfv372fQoEG6+JvPTm9vb4YOHZrkujVr1ug+h5VKJUuXLmXatGmsXbuWmJgYihcvzuTJk8mXL1+i686fP59ky6I3b6z69OmjS7be9H/s2DGOHTuWpP+3ky2NRsOBAwdo3rx5hk6QV2izYXXKy5cv07p1ayC+4m1KVnj8+eef/PLLL7i4uHD48OFkz/n111/ZsWMHLVq0YMKECWk23qZNm3Lr1i08PDxSVP34c6PSqBh5ZiT7HicUJfTI4cGSOkt0m1oLIYTI+p49e0b9+vVZtmxZshPgP3eHDx9m8ODBHDp0KFWT8d9I6ed3lpyz9SFvKhEXLlw4xUtp3yzJzZMnzzvPeXMsrSrcZhf6evpMrDKRZgWb6WK3Xt+iy4EuvI768KN9IYQQWYOLiwvNmjVLsudgVrFs2TLatm2bJonWx8iSrxHf5+bNm2zatAmIf1ycUiEhIQDvLdj25lhKJixu2rSJLVu2pKjvR48epei8z5lST8noSqMxUhqx4W78pMz7QffptL8Ty79djr1Z2laCFkIIkTnGjh2b2UP4ZJs3b86UfrNVshUQEEDfvn1RqVTUqVPnne+3k/Omau376mi9WSURHf3hFXf+/v7cunUrxf1nB3oKPYZVGIaJvgl/3PwDAO9Q7/iEq+7yRHssCiGEEF+KbJNsvSlm9uLFCzw8PJgyZcpHXf+mam1yy3HfeFPhNiWF7uzs7HSrST7k0aNHKUrgsgKFQkH/Mv0x0Tdh/tX4iY3Pw5/rnnC5Wrpm8giFEEKIjJUtkq2IiAi6du3K7du3KViwIH/88cdHb3vwplDem9eJyXlz7F2Vjd/WqlWrRJuXvs+bCXbZhUKhoEfJHhjrGzP93/gd519FvKLj3x1Z9u0yCtoUzOQRCiGEEBkny0+Qj4qKokePHly9epW8efOycuXKTyrB/6Zq7ZMnT955zpuib3nz5v2ksX5pOnp0ZMRXI3Tfv45+TZcDXbj1OvsklkIIIcSHZOlkKyYmhp9++omLFy/i5OTEqlWrkt1vKiVKliwJxG9v4Ovrm+w5ly5dAuKr6IqUaVW4FeOrjEdPEf9PLTgmmK4HunLF70omj0wIIYTIGFk22YqLi6Nv376cO3cOe3t7Vq9ejYODwye35+bmRqFChYDkVyucO3eOJ0+eYGBgkKi6vPiwJgWaMLX6VPQV8W+tw+PC6XGoB+dfns/kkQkhhBDpL0smW2q1msGDB3PixAns7OxYvXq1boPPD6lZsyY1a9Zk//79SY716dMHiK/D8fa2BY8fP2bkyJEAtGnT5ovfqudT1Mtbj1k1ZmGgF7/aM0oVRe/DvTn5/GQmj0wIIYRIX1lygvzff//NgQMHgPhyDMOHD3/nuaNGjUq0/9GbneEjIyOTnFu3bl06duzI6tWr+emnn8iTJw+mpqY8ePAAtVpN2bJlGTx4cBrfzZfjG5dvWFBrAf2P9SdKFUWsJpb+R/szpfoU6uatm9nDE0IIIdJFlky23pRggPjk6U0ClZywsLCPanv48OGULl2aDRs2cOfOHfz8/MifPz+NGzemU6dO763DJT6skmMlFtdeTK8jvYiIi0ClVTH05FBi1DE0zt84s4cnhBBCpLksmWw1bdo0RRt+JufevXsfPKd+/frUr1//k9oXH1bGvgx/fPsHPQ73ICQmBI1Ww4jTI4iKi6Jl4aSbgAshhBBZWZZMtkTW55HTgxV1V9DtYDcCowMBmHBhAlGqKDoV65S5gxNCZIgdO3bw66+/JnusW7duDBkyhJo1a1KwYEGWLFmS5JwLFy7QoUMH5syZQ7169YD4uot//PEH165d48aNG4SEhDB58uSP+gO9Zs2aujcmCoUCc3NzHBwcKFWqFM2bN9etXn+bu7t7sm3lzJmTwYMHv/M+3+bk5JRovrDIPiTZEpmmkE0hVtdbTdeDXfGNjC+3MePSDCJVkfxU8icUCkUmj1AIkRH69euHs7Nzotib1eEfKygoiAULFuDo6Ii7uzv//PPPJ7VTpEgROnfuDMQncI8fP2b//v1s2bKFTp06JZs8ValShe+//z5RzNjYmKJFizJt2rRE8ZEjR1KiRAk8PT11MTMzs08aq/j8SbL1BdNqtcR6eWGUL1+mjSGvVV5W119N1wNdeR7+HIBF1xYRERfBkHJDJOES4gtQvXp1ihcvniZt5cqVi9OnT2NnZ8eNGzdo3rz5J7Vjb2+fJHEaMmQIgwcPZtWqVbi6utKmTZtEx/PmzZvkmjf+f8X8mDFjcHFxeef5InvJkqUfRNp4vXQZXt83IXj7jkwdh5O5E6vqrcLNyk0XW3N7DePOj0OtUWfiyIQQWY2hoeEnF7f+EGNjY6ZNm4a1tTWLFy9Gq9WmSz8i+5Fk6wsVdvw4/rNno42L4+WIEbyaOAntezbhTm/2ZvasqreKwraFdbFt97cx4swIVBpVpo1LCJH+wsPDCQwMTPT1NpVKleR4YGDgR682TwtmZmbUrl0bX19fHj58mOhYTExMkjG+vXpefLnkNeIXSmlhgdLWFvXr1wAErV1LzP37OM2ehf4n7C2ZFmyNbfmj7h/0OtyLa/7XANj3eB9RcVH8/vXvGCoNM2VcQnyuomLVxKo1mToGQ6UeJobKVLXRqVOnJLG3V46fPn2aSpUqpaqPtFSwYEEgfr/cN/8bYNu2bWzbti3RuR87OV9kT5JsfaFMy5bFbdtWnvfuQ/Tt2wBEXriAd/MWOC+Yj3Hhwh9oIX1YGlqytM5S+h3tx4VXFwA4+uwofY/2ZXaN2Zjom2TKuIT43Izdc4vVZ73RZPKbLD0FdKycl98aeXxyG6NHj8bNze2dx0uWLMmAAQOSxO/evcvUqVM/ud9P9WYie0RERKJ4rVq1aNeuXaJYgQIFMmxc4vMlydYXzMDBAdcN63k5ajShe/YAEOfjg3frNjhOnoTlf0upM5qpgSkLai9gyPEhHH9+HICzL87S81BP5teaj4WhRaaMS4jPyZpzTzI90QLQaOPHkppkq0SJEu+dIG9jY0PlypWTxJXKT3+iFhYWRnR0tO57AwMDrK2tU3TtmyTr/1cP5s6dO9lxCiFztr5wesbGOE6bSq6ffwa9+H8O2qgofAYMxG/WbLSazHlFYaQ0YmaNmdTLm5DwXfa7TNeDXQmODs6UMQnxOelQyRW9z2CxrlJPQYdKrpk9jI82ceJEqlatqvvq27dviq998OABAK6uWe++ReaQJ1sChUJBjh+7YFSoED6DB6MJDQXg9ZIlxNy9i9PMGehlQv0XAz0DplSbgqmBKTsexK+YvP36Np0PdGZpnaXYmabPiiMhsoLfGnkwtG7hbDFnKzN07dqVxo0TtgiztLRM0XUREREcPnwYBwcH8ufPn17DE9mMJFtCx7xaVdy2bOZZ7z7EPnoEgDYuDoWRUaaNSamn5LdKv2Gqb8q6O+sAeBj8kI77O7Ls22U4mTtl2tiEyGwmhkpMyHqJzuegQIECHz2fKjo6mqFDhxIcHMzAgQOlDqBIMUm2RCKGefOSd/MmXgz9hZiHD3GaOQOFfub+M9FT6DG0/FBMDUxZen0pAM/CntHx7/iE6+36XEIIsW7dOkJDQ/Hz8wPg2LFjvHr1CoD27dtjYfHheZ++vr7s2rULgMjISB49esT+/fvx9/enS5cutGrVKv1uQGQ7kmyJJJTm5jjPn4cqIABlCieMpjeFQkHf0n0xMzBj1qVZAPhG+tJpfyeW1FmSqD6XEOLLtmLFCt3ehgAHDx7k4MGDADRu3DhFydadO3cYOnQoCoUCMzMzHBwcqFGjBi1atKBEiRLpNnaRPSm0UgI30zVt2pRbt27h4eHBjh0ZU8194z9PufYsmHYVXSnmZJXi68JPnyHkzz9xGDcWPVPTdBzhu225t4UJ5yegJf6froWBBQtrL6RUrlKZMh4hhBBfppR+fstqxC/QvVdh/LrjBpsuPqPpwrOsPeedom0nYp88wWfQIEL37sW7TVtin/t88Jr04OnuyaRqk1Aq4ueqhMWF0f1Qd86/PJ8p4xFCCCHeR5KtL5CliT4mBvGJSqxaw6hdt+i78QrhMe/fFuf1qlW6lYoxd+/i3aIFERf+SffxJqdhvobM/GYmBnoGAESpouh1uBfHnh7LlPEIIYQQ7yLJ1hfIwcqEFZ3Kk9M8YZXh3usvaTzvNHdehr7zutwjRmDbsaPue3VQEE+7dCFw7bpM2ZC1Zp6aLKi1QFdVPk4Tx8DjA9n3eF+Gj0UIIYR4F0m2vlCV8ufgr/5VqZjPVhd7HBBBkwVn2HLxWbLJk0JfH/tfh+EwZTIKw//2KVSr8Z04kZcjR6LJhA1XKzlWYmmdpVgYxE94VWvV/HrqV7be35rhYxFCCCGSI8nWFyyXhTHru1akb80CvCkXE6PSMHT7dYZsvU5kbPKvFa2bNMF1/Tr07e11sZDtO3jSvj1xvn4ZMfRESuUqxR91/8DGKH4DbS1axp0bx6qbqzJ8LEIIIcT/k2TrC6fUUzD4W3dWda6ArZmhLr798nO+n3+Gh35hyV5nUrw4btu2YlKmjC4Wfe063s2bE3X1anoPO4kiOYqwqt4qcpnm0sVmXJrB3MtzM+UVpxBCCPGGJFsCgK8L2bGvX1XKudroYg/8wmk07wx/Xkl+1aG+nR2uq1Zi7empi6n8/Xnx63C0qvdPtk8P+azzsbreapzNnXWxZTeWMfmfyWi0mbuliRBCiC+XJFtCx8HKhI3dK9Kjej5dLCpOzYDNV/l1xw2i49RJrlEYGuIwbiy5x/wG+voojI0zteq8s4Uza+qvoYB1wjYcG+9uZOTpkag0GZ8ACiGEEJJsiUQMlHr8+l0RlnUoh6VxQsK08Z+nNFt0lievI5K9zqZVK1xXrcTx92kYFymSUcNNlp2pHSvrrqRYjmK62J7Hexh8fDCx6oyfxC+EEOLLJsmWSFadovbs61eNEs4J1eVvvQil4bzT7L/5KtlrTMuVw7JOnSTx2Oc+RN++nW5jTY61sTXL6y6nfO7yutjRZ0fpfaQ3kXGRGToWIYQQXzZJtsQ7udiasrVnJTpWctXFwqJV9Fx3ifF7bxOr+vA8KE1kJM9798a7TVtC9mZs/SszAzMW1lrIN87f6GLnX56n+6HuhMSEZOhYhBBCfLkk2RLvZaSvZOz3xZjXujRmhkpd/I/TXrRaeo4XwVHvvf7lqNHE3LuHNjqaF0OG4Dd9Olp10rlf6cVY35iZNWZS362+LnbN/xpdDnQhICogw8YhhBDiyyXJlkiRRiUd2d23KoVzW+hil58G02DuKY7fe3dtLaumP6BnlfAq8vXyP3jWoyfqkIx7smSgZ8DkqpPxLJSwavJ+0H067e/Ei/AXGTYOIYQQX6bMWTImsqT8dubs7FWF0btusvXScwCCIuPovOoifWoUYEDtQij1FImuMa9SBbdtW3neqzcxDx4AEHH6NF4tPHFZMB+jggUzZOxKPSUjK47EzNCMlTdXAvAk9Akd/u7A0m+Xks8q3wdaEEKktR07dvDrr78me6xbt24MGTKEmjVrUrBgQZYsWZLknAsXLtChQwfmzJlDvXr1ALh+/Tp//vknFy5cwMfHB2tra0qWLMmAAQNwc3NL0bhq1qyJj098yRuFQoG5uTkODg6UKlWK5s2bU7JkySTXuLu7J9tWzpw5GTx48Dvv821OTk4cPXo0RWMUWYskW+KjmBgq+b1FScq72TJ6102i4zRotTDv6EMuPQliTqvS2FkYJbrG0MWFvJs28uLX4YQdPAhA3NOneLdshcPUKclOqk8PCoWCQWUHYWloyZzLcwDwjfSl09+dWFxnMUVzFM2QcQghEuvXrx/Ozs6JYoUKFfqktpYvX87ly5epV68e7u7u+Pv7s379epo2bcrmzZtT3G6RIkXo3LkzABERETx+/Jj9+/ezZcsWOnXqlGzyVKVKFb7//vtEMWNjY4oWLcq0adMSxUeOHEmJEiXwfKtOoZmZ2cfersgiJNkSn8SznAslnK3ote4yjwPiy0GcffSaBnNPMb9NGSq42SY6X8/MDKc5s3m9ZAn+c+aCVosmMhKfvv2I6dWLnH16o9DLmLfaXYt3xcLAgokXJqJFS1BMED8e+JH5teZT1r5shoxBCJGgevXqFC9ePE3a6tSpE9OnT8fQMGFHjO+++45GjRqxdOlSpk+fnqJ27O3tkyROQ4YMYfDgwaxatQpXV1fatGmT6HjevHmTXPOGi4tLou/HjBmDi4vLO88X2YvM2RKfrHBuS3b1qUKDEg66mF9YDK2XnWfxiUdJtslRKBTk7NkT54UL0DM318UDFi7kxZAhGTZugJaFWzK52mSUivhJ/+Fx4fQ41IOTz09m6DiEEGmrTJkyiRItiE+CChYsyOPHj1PVtrGxMdOmTcPa2prFixfLVmAixeTJlkgVC2MD5rcuTXlXGyb+dYc4tRa1RsuUv+/yr3cQM1qUxMrUIPE1NWqQd8tmnvfuQ6yXFwDmNWpk+Ngb5GuAuYE5g08MJkYdQ4w6hv5H+zO52mTqudXL8PEI8bGiVFHEaeIydQwGegaY6Jukqo3w8HACAwMTxWxtE56Oq1SqJMcBwsKS37v1/2m1WgICAiiYBnNEzczMqF27Ntu2bePhw4eJ2oyJiUkyTnNz8yTJn/jySLIlUk2hUNCpihslXazps+EKPv+Vgzh8x5eG80+xsE1Zir9VHBXAKF8+8m7ZzIshP2OYLx9WjRplxtD52uVrFtVeRN+jfYmIi0ClVTH05FBCY0PxdPf8cANCZJKp/0xlw90Nmb7vp55CjzaF2/BLhV8+uY1OnTolid27d0/3v0+fPk2lSpU+uf3du3fj6+tLv379PrmNt71JsJ4+fZoo2dq2bRvbtm1LdO7kyZNp2rRpmvQrsi5JtkSaKZ3Hhr19qzJwy1WO3/MH4FlgFM0WnWV0o6K0/SoPCkXCakWlhQXOixaCJumHhTY2FgwMEp2fXsrnLs8fdf+g56GeBMcEo0XL+PPjCY0NpWvxrunevxCfYuPdjZmeaAFotBo23t2YqmRr9OjR710p+GY14f+7e/cuU6dOfW/bjx49Yty4cZQuXZoffvjhk8f4tjcT2SMiEm9fVqtWLdq1a5coVqBAAYSQZEukKRszQ1Z0LM+iE4+YcfAeGi3EqjWM/PMm/3oHMqlpcUwNE/7ZKfT04P8mxmu1Wl6MGIk2Lg7HiRPQy4AVOh45PFhdbzXdDnXDLzK+bticy3MIiw1jQJkBGZL0CfExWhdu/Vk82VIqlLQu3DpVbZQoUeK9E+RtbGyoXLly0r6VymTOTuDv70+PHj2wsLBgzpw5ic4PCwsjOjpa972BgQHW1tYpGu+bJOv/Vw/mzp072XEKIcmWSHN6egp61yhA6TzW9Nt4hYDw+M2f/7z6glsvQlnUriwFcpm/8/rAlasI3bMHAO/Hj3FeMB/D/1vJkx7yWedjTf01dD/YnadhTwFYcXMFobGhjPxqJEq99/9iFyIj/VLhF/qV6Zct5mylh7CwMLp160ZYWBjr16/H3t4+0fGJEyeyc+dO3fcVKlRg7dq1KWr7wX81A11dXT9wphDxJNkS6aZy/pz81a8afTZe4R+v+EmjD/zCaTz/NFOalaBxScck12g1GiLOn9N9H3P/Pt7NW+A0ayZmGfAXo5O5E6vrr6bHoR7cD7oPwLb72wiLDWNy1ckYKA0+0IIQGcdE3wQTPr9EJ7PFxMTQs2dPvL29WblyZbKv8rp27Urjxo1131taWqao7YiICA4fPoyDgwP58+dPszGL7E1KP4h0lcvSmA1dv6Ln1wm/lCJj1fTbeIXRu24So0q8T6JCTw+XhQux7dJFF1OHhPC0azder1yVIUutc5rkZEXdFZSyK6WLHfA+QN+jfYmMi0z3/oUQn06tVjNgwACuXr3KnDlzKF26dLLnFShQgMqVK+u+ihUr9sG2o6OjGTp0KMHBwfTs2VOmF4gUkydbIt3pK/UYVr8wZV1tGLTlKmHRKgDWnHvCtechLGhTGmcbU935Cn197If+jHGRIrwcORJtTAxoNPhNnUr07ds4jB+HnrFxuo7ZysiKJXWWMOj4IM68OAPAmRdn6H6oOwtqLcDKyOoDLQghMsOUKVM4evQoNWrUIDg4mF27diU6ntIior6+vrprIyMjefToEfv378ff358uXbrQqlWrNB+7yL4k2RIZpk5Re/b1rUavDZe46RMKwLVnwTScd5pZLUtRwz1XovOtGjXEMJ8bz/v2RfXiJQChe/YQ++gRzvPnYeCY9DVkWjI1MGVezXkMPz2c/d7748frf43OBzqzpPYS7Ezt0rV/IcTHu3v3LgDHjh3j2LFjSY6nNNm6c+cOQ4cORaFQYGZmhoODAzVq1KBFixaUKFEiTccssj+FVkrgZrqmTZty69YtPDw82LFjR2YPJ91Fx6kZt/c2Gy48TRTvWzP5zaxVgYH4DBhI5D//6GJKW1ucZs/CrEKFdB+vWqNm4oWJbL2/VRdzNndm6bdLcbFI/4n7QgghPk8p/fyWOVsiwxkbKJn0Q3FmepbExCBhhd+8ow9p/8cFAsJjEp2vb2tLnj+WY9OhvS6mDgwk/MiRDBmvUk/JqIqjEtXceh7+nI5/d+RB0IMMGYMQQoisS5ItkWmalnHmz95VyGeXUKvmzWbW/3on3vJCYWBA7uHDcZg0CYWhIaYVKpArA/dTVCgU9C/Tn8FlB+ti/lH+dNrfiWv+1zJsHEIIIbIeSbZEpnLPbcHuPlVp+NZm1r6hMbRaep7lpx4nWX1o3fQHXNevw2n2LBQGGV+GoVOxToytPBY9RfyPTmhsKN0OduOsz9kMH4sQQoisQZItkenMjfSZ17o0YxoVxUAZP19LpdEyYd8deq2/TFh04qKNJsWLo//WJrVvBG3dSuSlS+k+3qYFmzL96+kY6MUne1GqKHof7a2bRC+EEEK8TZIt8Vl4s5n15h6VcLBKKOvw981XNJ5/hjsvQ997fcT5C7waM5YnHTsRtHFjutfjquNahwW1FugqZ6s0KoaeGMqWe1vStV8hhBBZjyRb4rNS5r/NrKsVzKmLeQVE8MPCM2y79DzZazQxMbz45RdQq0Gl4tXYcbwaPRpNbGy6jrWSYyWWf7tcV3PrzQbWy64vy5Diq0IIIbIGSbbEZyeHuRGrOlegf62CvCnQHB2nYcjWawzbfp3ouMRV5/WMjHCaPQt9u4S6V8Fbt/G0fQfifP3Sdawl7Eqwut5qcpkm1Aibe2Uuv//7e6ZvECyEEOLzIMmW+Cwp9RQMrFOIlZ3KY2OaMBF+08VnNFt0lqevE2+bY1q6NHm3bcOkZEldLOraNbybNyfyypV0HWt+6/ysrb+WvJZ5dbG1t9cy6syoTN8kWAghROaTZEt81r5xz8XeftUo5WKti916EUrDeac4fNs30bkG9rnIs3YN1i2a62Iqf3+edOhI0NatpCdHc0dW1VtFEdsiutjuR7sZdGwQ0arodO1bCCHE502SLfHZc7I2YUuPSnSs5KqLhUar6LrmX6buv4tKnfC6Ts/QkNzjxpF7zG+g/99uVHFxvBo1mpdjxqBNx3lcOUxysKLuCsrZl9PFjj8/Ts/DPQmLDUu3foUQQnzeJNkSWYKhvh5jvy/G3NalMTVMqDq/6Pgj2v/xD/5hCVXnFQoFNq1a4bp6FcqcCRPtgzdtxm/GjHQdp7mhOYvrLKaGSw1d7JLvJboc6EJAVEC69i2EEOLzJMmWyFIal3RkV+8qFMhlroudexxfdf7i/1WdNy1bFrdtWzEuXhwA/dy5ydGtW7qP0UhpxMxvZvJ9/oQNb+8G3qXj3x15Hpb8ikohhBDZlyRbIsspaG/Brt5VaFTSURfzC0u+6rxB7ty4rluLdcuWOM+bh/5bT7rSk76ePuOrjKdj0Y662NOwp3T4uwP3g+5nyBiEEEJ8HiTZElmSmZE+c1uVYmxjD13VefU7qs7rGRnhMHYMJsWLJWkn5vHjdJvHpVAoGFxuMAPKDNDF3uyneMUvfVdICiGE+HxIsiWyLIVCQcfKedncoxKOyVSdv/vq/VXn41684Em79jzp3AWVv3+6jfHH4j8m2k8xLDaMbge7cfL5yXTpUwghxOdFki2R5ZXJY8PeftWSVJ1vsuAMOy6/o+p8dDTP+/RFHRhI1KVLeDVvQdT16+k2xqYFmzLzm5kY6hkCEKOOod/Rfux5tCfd+hRCCPF5kGRLZAu2ZobJVp0ftOUaI3beIEaVuOo8CgVGhQrpvlX5+vKkbTuCt29PtzHWylOLxXUWY2ZgBoBaq2b46eGsubUm3foUQgiR+STZEtnG21Xnrd+qOr/+wlNaLD7Hs8CEqvN6RkY4TJ6E/YgRoIwvJaGNi+PliJG8GjcebVz6VH4vn7s8K+uuxNbYVhf7/d/fmX1ptuynKIQQ2ZQkWyLb+cY9F3v7VqWks5Uudv15CI3mn+bYvYS9EhUKBbbt25Fn5QqUtgnJT9CGDTzp3BlVQPrUxSqSowhr66/FydxJF/vj5h+MPTcWlUaVLn0KIYTIPJJsiWzJ2caULT0r0a5iHl0sODKOLqsuMvPgPdSahKdIZhUq4LZ9G8bFElYrRv17Ca9mzdNtHlceyzysrb+WQjYJrzK3P9jOkBNDiFHHvOdKIYQQWY0kWyLbMtJXMqFJcWZ6lsTYIP6fulYLc48+pNPKf3gdnpDUGDg44LpuLVZNmuhiKl9fnrRrT8i+fekyPjtTO1bWW0mZXGV0sSNPj9DzkGzvI4QQ2UmWTbb8/f35888/mTBhAi1btqREiRK4u7vTvn37T25z3rx5uLu7v/dr48aNaXgXIiM0LePMn72r4JbTTBc79SCARvNOc+VpkC6mZ2ycZB4XWi2GLi7pNjZLQ0uW1FnCN87f6GL/+v5L5/2d8Y9Mn3IUQgghMpZ+Zg/gU+3bt4/JkyenS9s5cuTA1dU12WN2dnbp0qdIX4VzW7K7TxV+3nqd/bdeAfAiJBrPJecY3bAo7Sq6olAodPO4jNwL4TNgILkGD8KkRIl0HZuxvjGzasxizNkx7Hq0C4B7Qfdo/3d7ltZZSh7LPB9oQQghxOcsyyZb5ubmVK5cmeLFi1O8eHFu377NwoUL06Tt6tWrM2XKlDRpS3w+LIwNWNSuDMtPeTFl/13UGi1xai2jdt3i0pMgJjUtjqlh/I+EWYUK5P/7L5RWVkna0Wo0KPTS9qHwm+19bE1sWXlzJQA+4T60/7s9i2ovomiOomnanxBCiIyTZV8jNm/enJUrVzJo0CDq1KlDjhw5MntIIgtQKBR0q56Pjd0qYmdhpIv/efUFTRac4ZF/uC6WXKKliYnhSfsO6VKPS6FQMKjsIIaUG6KLBUYH0uVAFy68vJDm/QkhhMgYWTbZEiI1KrjZsq9fVSq4JZR8uO8bTuN5p/nrxstkr9FqtbwaO46oS5d4OWIkL8eOTZd9FTt6dGRS1UnoK+KfskXERfDT4Z846H0wzfsSQgiR/iTZSsbdu3cZPHgwHTp04KeffmL27Nk8ePAgs4cl0lguC2M2dP2KHl/n08UiYtX0Wn+ZCXtvE6fWJDpf5edP+IkTuu+DN27iSafO6bKvYqP8jZhbcy7Gyvg9H+M0cQw5MYTNdzeneV9CCCHSlyRbybhz5w579+7lwoULHD16lEWLFtGoUSMmTZqEWq3+cAMiy9BX6vFr/SIsblcWC6OEKYzLT3vRZtl5fEOjdTED+1y4bduKcfHiuljU5cvx9biuXk3zsVVzrsayb5dhaWgJgBYtEy5MYOHVhVJtXgghspAsO0E+PeTKlYt+/fpRrVo1nJ2dMTc3x8vLiw0bNrBp0yZWr16Nvr4+Q4cO/WBbmzZtYsuWLSnq99GjR6kdukilesVy457bgp/WXeLuq/gaVxe9g2gw9zTz25SmYr74OYFv6nG9GjuOkB07AFD5+eHdvgO5R43ExtMzTcdVKlcp1tRfQ49DPfCN9AVg0bVFBEQFMOKrESj1lGnanxBCiLQnydZbWrZsmSTm7u7O2LFjcXZ2Zvr06axevZo2bdrg7Oz83rb8/f25detWeg1VpAO3nGbs7FWFEX/eYMdlHwACwmNou/wCQ+u60716PhQKRfy+ihMnYFzMA99Jk0Glgrg4Xo3+jeibt7AfOQI9Q8M0G1d+6/ys+24d3Q91xyvEC4Ct97cSFB3ElOpTMFIafaAFIYQQmUleI6ZQly5dyJUrFyqViqNHj37wfDs7Ozw8PFL0ZWxsnAF3IFLCxFDJjBYlmfhDMQyV8T8eao2WyX/fpcfaS4RGx29QrVAosG3TBtfVq1DmzKm7PnjLFp526Eicr1+y7X+q3Ga5WVNvDSXsEmp+HX56mB6HehAaG5qmfQkhhEhb8mQrhZRKJSVLluTQoUM8efLkg+e3atWKVq1apajtpk2bylOwz4hCoaDtV64Uc7Si1/rL+ARHAXDwti+N551mUbuyFHGIn0dlWrYsbtu38bxvP6L/20cx6upVIs6exfqHJmk6Lmtja5bVWcbPJ3/m5POTAFzyvUSn/Z1YXHsxuUxzpWl/Qggh0oY82foIBgYGAKhUqkweicgIJV2s2du3KtULJewa4P06kh8WnmHH5ee6mIG9ffy+is2bAWDdsmWaJ1pvmBqYMrvGbL7P/70u9iDoAe3/aq97xSiEEOLzIsnWR3hT/iF37tyZPBKRUWzMDFnZqTwDahdEoYiPRcdpGLTlGiN23iBGFb86Vc/QEIfx43GaPRv7EcPTdUwGegaMrzKersW76mIvIl7Q4e8O3PC/ka59CyGE+HiSbKXQ8ePHdclWlSpVMnk0IiMp9RQMqF2IlZ3KY21qoIuvv/AUz8XndK8ZFQoFlvXqJpkcr9VqCVi2jLhXr9JsTAqFgv5l+vNL+V90seCYYH48+COnfU6nWT9CCCFS74tLtlq3bk3NmjVZtWpVoviDBw8YPXo0d+/eTRTXaDTs3buXwYMHA1CjRg1KpPPGxOLz9I17Lvb2rUoJ54RtfK49D6Hh3FOcuP/uwqbBmzfjP2MmXs2aE3nxYpqOqV3RdkyrPg19vfjpl1GqKPoe6cueR3vStB8hhBCfLstOkH/58iVNmjTRfR/737Yply9f5quvvtLFu3btSrdu3XTf+/r64uPjQ1hYWKL2VCoVmzdvZvPmzVhbW+Po6IhSqeTp06eEhIQAUK5cOaZNm5aOdyU+d842pmztWYlxe26z/sJTAIIi4+i08h8G1CpE35oF0NNT6M6PffaMVxMnAaB+/Zonnbtg/8sv2LRri0KhSLaPj1XfrT7WRtYMODaASFUkKq2K4aeH8zrqNR09OqZZP0IIIT5Nln2ypVarCQ4O1n1FRkYC8UnT2/Ho6OgPtBTPycmJAQMG8M0332BhYcGTJ0+4e/cuBgYGVK9enWnTprFmzRosLS3T87ZEFmCkr2TiD8WZ0aIkRvrxP0JaLcw6fJ8uqy8SHJmwX6KhiwsOv41G8d/iClQqfCdO5OWwX9Gk8N9mSlRyrMSKeiuwNU7Y63HGpRn8/u/vaLSa91wphBAivSm0su9HpntT+sHDw4Md/1UlF1nD7Reh/LT+Ek9eR+piTtYmLG5XluJvvW6MunaN5/36o/L11cWMixbFed5cDJyc0mw8T0Of0vNwT56FPdPF6uetz4SqEzBUpl2hVSGEECn//M6yT7aE+BwUdbRkd5+q1Clqr4v5BEfRbNFZNv7zVLeHoUnJkrht34ZJubK686Jv38arWXMizp1Ls/HksczDmvprKJqjqC72t/ff9Drci/DY8DTrRwghRMpJsiVEKlmZGLCkXVl+qVeYN9O1YtUaft1xg6HbrhMdF18eQj9nTlxXrsSmfXvdtergYJ7+2JXXf6xIs82lc5rkZEXdFVR2rKyLXXh1gc4HOuMf+e6J/EIIIdKHJFtCpAE9PQU/fZOfdV2/Iqd5wuu6rZee88PCszx5HQGAwsCA3COG4zh1Cgqj//Y01Gjw+/13grdtS7PxmBmYMb/mfBrka6CL3Q28S/u/2+Md4p1m/QghhPgwSbaESEOV8+dkb99qlHW10cXuvAyl4bzTHL6dMF/L6vvvybtxAwaOjkD8a0ar779P0l5qGCgNmFR1Ep09OutiPuE+dPi7A9f9r6dpX0IIId5Nki0h0lhuK2M2da9I5yp5dbGwaBVd1/zL7wfuotbEvy40LlqUvNu3YdmgAU5z5yQphpoW9BR6DCo3iJ/L/ayLBcUE0fVgV93+ikIIIdKXJFtCpAMDpR6/NfJgXuvSmBoqdfEFxx7RYcUFAsJjANC3scFpxnQM7O2TtBF18xZaTdqUbejg0SFJ8dN+R/ux88HONGlfCCHEu6VZUVMvLy+uXr2Kn58fgYGBxMTEYG1tja2tLfnz56dMmTKYmJikVXdCZAmNSjpSxMGCHmsv8cg/ft7WmYevaTj3NAvblaFMHptkr4u8coUnHTpi/nV1HKdMQWlunuqx1Herj42xDQOODSAiLgK1Vs3os6PxjfSlR4keUvxUCCHSSaqSrStXrrBlyxZOnTrF69ev33uuUqnEw8ODRo0a8f3332NhYZGaroXIMgrksmBXn6r8sv06+66/BOBVaDQtl5xjZIOidKjkmijRifPzw6dff4iLI/zwEbxbeOK8YD5G+fKleiwVHSqysu5Keh3pRUBUAAALri7AL9KP4V8N1z35EkIIkXY+6TXirl27aNSoEW3atGHnzp0EBASg1WoxMTHBwcGBIkWKUKpUKdzc3LC1tUWhUKBSqbh27RoTJ06kevXqjBw5kpcvX6b1/QjxWTI30md+69KMblgU/f/qQ8Sptfy2+xYDNl8lMlalO1fPyAijokV038d6eeHdwpOww4fTZCxFchRhbf215LXMq4ttvb+VgccHEqWKSpM+hBBCJPioCvIXLlxg6tSp3LlzB61Wi5WVFXXr1qVcuXKULFkSV1fXZK+LiIjg5s2bXLt2jaNHj3L16lUAjIyM6NChAz169MA8DV6TZFVSQf7L8q93IL03XMY3NEYXK2RvzqJ2ZclvF/9zoNVoCJg/n4CFixJdm6NHD+z69UWhVJJaQdFB9D3al2v+13SxEnYlmF9zPjbGyb/eFEIIkSCln98flWwVLlwYgKpVq9KqVSu+/vprDN7s+fYRnj17xq5du1i3bh0hISH06dOH3r17f3Q72YUkW18e/7AY+m68zPnHgbqYuZE+vzcvQf3iDrpY2NGjvPh5KJqICF3MrGpVnKb/jtLaOtXjiFJFMfTkUI4/O66L5bXMy6Lai3C2cE51+0IIkZ2ly3Y9VatWZfPmzSxfvpzatWt/UqIF4OLiQp8+fTh27BiDBg3COg0+NITISuwsjFj341f0/Dq/LhYeo+Kn9ZeZuO82cer4VYgWNWuSd+tWDPMnnBdx+jRezVsQffduqsdhom/CrG9m0aJQC13MO9Sb9n+3587rO6luXwghxEcmW8uXL6dkyZJp1rmJiQndunWjbdu2adamEFmFvlKPYfULs6R9WSyMEiamLzvlRdtlF/ALjQbAKJ8beTdvxqJOHd05cc+f492qNeEnU18rS19Pn1EVR9GnVB9dLCAqgE77O3H2xdlUty+EEF86qbMlRCar65Gb3X2rUjh3wgrdf7wDaTDvNBcex6/yVZqb4TR3DnYDB8J/Kxf1zM0xcndPkzEoFAp6lOzBuMrjUCri54NFqiLpfbg3ex7tSZM+hBDiSyXJlhCfAbecZuzsVYWmpZ10Mf+wGNosv8Cyk4/RarUoFApy9uiOy7JlKHPmxHn2rGSLoabGDwV/YF7NeZjox9fEU2lVDD89nOU3lqfZRtlCCPGlkWRLiM+EiaGSGZ4lmdCkGIbK+B9NtUbLxL/u8NO6y4RFxwFgXrUKBQ4dxLRcuSRtaGJjUz2Oas7VWFl3JbbGtrrYnMtzmHhhImqNOtXtCyHElyZNKxiq1Wq8vb0JDQ1FpVK999zy5cunZddCZAsKhYJ2FV0p5mRF7/WX8QmOr3u1/9Yr7vmGsbhdWdxzW6CXzG4MKn9/vFu1Jkf37lh7tkhVRXiPnB6sq7+On478xJPQJwBsvrcZ30hfplWfpnvyJYQQ4sPSJNl69eoVM2fO5NChQ0RHR3/wfIVCwe3bt9OiayGypVIu1uzpW5X+m65w6kF8pXevgAiaLDjD5KbFafLW60YAbWwszwcMJM7Hh1e//UbUjevkHjUKPSOjTx6Di6ULa+qvoe+RvlwPuA7A8WfH6XqgK/NqzUv05EsIIcS7pfo14rNnz2jevDl79uwhKioKrVb7wS9NGm2uK0R2ZmtmyKrOFehXq6AuFhWnZsDmq4zedZMYVcIrvVgfH2KfPtF9H7JtO0/atiPuxYvUjcHYluV1l/ONyze62PWA67T/qz3PQp+lqm0hhPhSpDrZmjlzJgEBAdjY2DBx4kROnDjBrVu3uHv37nu/hBAfptRTMKhOIVZ2Ko+VSUJduzXnntByyXle/Pea0cjNDbdt2zEpXVp3TvTNm3g1a07E+fOpGsObWlwt3VvqYk/DntLu73bcDLiZqraFEOJLkOpk6+zZsygUCmbOnEmzZs2wt7dHmQZbiQghEtQonIu9fatS3MlKF7v6LJgGc09x6oE/AAb2uXBdvQrr1q1056iDgnja5Ude/7EiVasJ9fX0GfHVCPqX6a+LBUYH0uVAF04+T32tLyGEyM5SnWzFxsZibGxMxYoV02I8Qoh3cLE1ZWvPSrSukEcXC4qMo8OKf5h/9AEajRaFoSEOv/2Gw6RJKAwN40/SaPD7/Xd8Bg1KtO3Px1IoFHQt3pVJVSehr4if7hmliqLv0b5su78tVfcmhBDZWaqTLWdnZ6m/I0QGMTZQMrlpcX5vXgIj/fgfX60Wph+8T9c1/xISGV8ewrrpD7hu2IC+41v7LP69H6+WLYnx8krVGBrlb8TC2gsxMzADQKPVMPbcWOZeniu/C4QQIhmpTrbq169PTEwM586dS4vxCCFSoEU5F3b0qkweW1Nd7OhdPxrMO8VNnxAATIp54LZ9O6aVEp46x3o/Qf36dar7r+RYidX1VpPLJJcutuzGMoafHk6cOi7V7QshRHaS6mSrS5cuFC5cmFGjRvHsmaxOEiKjeDhasadvVWoXSUh4ngdF0XTRWTZffAqAvo0NeZYvJ0e3bgDY//JLssVQP4W7rTvrvltHAesCutjex3vpebgnobGhadKHEEJkBwptGjz3Dw4OZtSoUZw+fZq6detSvHhxzMzM3ntNkyZNUtttttG0aVNu3bqFh4cHO3bsyOzhiCxGo9Gy+OQjph+4h+atn2bPcs6M+74YxgbxC1airl/HuHjxVBU7TU5obCgDjw3kn1f/6GIFrAuwqPYicpvlTtO+hBDic5LSz+80KWrq4+NDQEAAUVFR7Nq1i127dr33fIVCIcmWEGlET09Br28KUMrZmr4br/A6In7Lni3/PuemTyiL25UlTw5TTEqUSHKtNjaWF8NHkOPHLhgXKfJJ/VsaWrK49mJGnR3Fvsf7AHgY/JC2+9qysPZC3G3TZrNsIYTIqlL9GvHu3bu0b9+eq1evAmBgYECuXLlwcHB451fu3PLXrhBprXKBnOzrV40yeax1sdsvQ2k47xSHb/sme43v1GmE7t2Ld6vWhHzgj6T3MVAaMLnqZLoV76aL+UX50XF/R876nP3kdoUQIjtI9ZOtefPmERkZiYuLC+PHj6dChQro6cn+1kJkhtxWxmzqXonJf99h5RlvAEKjVXRd8y+9a+RnUB13lHrxrxFjHj4kaONGALQxMbz4ZRhR165jP+yXhLIRH0GhUNCvTD8czB2YeH4iaq2aiLgIeh/pzehKo/mh4A9pdp9CCJGVpDorunLlCgqFglmzZlGxYkVJtITIZIb6evzWyIN5rUtjaphQYHjBsUd0WHGBgPAYAIwKFCDP8mUora115wRt2MCTjp2I8/X75P5bFGrB3JpzdZtVq7QqRp8dzcKrC6U0hBDii5TqzCgqKgoTExOKFSuWFuMRQqSRRiUd2d2nCvntEharnHn4moZzT3PpSSAAZpUr47Z9G8YeHrpzoq5cwatZMyIvXvzkvqs7V2dl3ZXkMM6hiy26toiRZ0ZKaQghxBcn1cmWq6srKpUKtVr94ZOFEBmqQC4LdvWpSsMSCcVNX4VG03LJeVae8UKr1WLg5ITrhvVYNWuqO0cdEMCTTp0JXL36k59GeeT0YH2D9bhZueliux/t5qfDP0lpCCHEFyXVyVaTJk2IjY3l6NGjaTEeIUQaMzfSZ17r0vzWqCj6/83XUmm0jN1zm74brxARo0LPyAiHCRPIPXYsCoP/NrxWq/GdPIUXg4egiYz8pL6dzJ1YW38tZe3L6mIXXl2g498deRn+MtX3JoQQWUGqk6327dtTsWJFRo8ezZUrV9JiTEKINKZQKOhcxY3NPSpib2mki++9/pLvF5zhoV8YCoUCm5aeuK5fh/5bK4ZD//qLsEOHPrlvKyMrltZZyndu3+liD4Mf0vavttx5feeT2xVCiKwi1asRFy1aRKlSpbh9+zZt2rShbNmyKSpq2qdPn9R2LYT4SGVdbdnXrxr9Nl7h7KP4bXse+oXTeP4ZpjQrQeOSjpiUKIHb9m34DBpM5IULWH3fGMvGjVPVr6HSkCnVpuBk7sSyG8sA8I/yp+P+jkz/ejrVnaun+t6EEOJzleoK8oULF9ZVpH7TVEoqVN+5I3/RviEV5EVGU2u0zDx0jwXHHiWKd6qcl+HfFcFQXw+tSkXgunXYtGyJnolJmvW9/f52xp8fj1obP89TqVAy/KvheLp7plkfQgiRETKsgnz58uVT24QQIoMp9RT8XLcwpV1sGLjlKmHRKgBWnfXm2vNgFrYtg4OVCTk6dUr2+qCNG7GoVw99G5uP7rtZoWbkNsvNoOODiFRFotaqGX9+PD7hPvQv0x89hZSPEUJkL2myN6JIHXmyJTLT09eR9Fx3idsvE1YI2poZMrdVaaoWzJnk/OCdf/Ly118xcHTEae5cTIp5JDknJe4F3qPXkV74RSbU9KqXtx4Tqk7ASGn0niuFEOLzkNLP74/+E/LcuXNoNJpUDU4I8fnIk8OUHb0q07Kciy4WGBFL+xUXmH/0AZq3dreOffqUV2PGABD34gVP2rQhePv2T+rX3dad9d+tp5BNIV1sv/d+uh7oSlB00KfdjBBCfIY+Otnq3LkzlSpVYujQoezfv5/IT1wSLoT4fBgbKJnavARTmxXHUD/+14JWC9MP3qfrmn8Jjozf3NrAxQW7fn3hv50itLGxvBwxkpejRqOJifnofnOb5WZ1vdVUcqiki131v0q7v9rxJPRJGtyZEEJkvo9Ottzd3QkJCWH37t0MHDiQihUr0qNHDzZv3oy/v396jFEIkUFals/Djp8qk8fWVBc7etePhvNOc/15MAqFghw//kieFX+gtLXVnRO8dStP2rYj7sWLj+7T3NCcBbUX0KxgM13sadhT2v7Vlku+l1J3Q0II8Rn46GRr165dHD16lBEjRlChQgU0Gg0nTpxgzJgxfP3113h6erJ06VIePXr04caEEJ+dYk5W7OlTldpFculiz4OiaL7oHOsvPEGr1WJWsWL8Nj8lS+jOib55E6+mzYg4e/aj+zTQM+C3Sr/Rv0x/XSwkJoRuB7ux7/G+1N2QEEJkslRPkA8NDeX48eMcPnyY06dPExkZqSv9kCdPHmrVqkXNmjUpW7ZsikpCfIlkgrz4HGk0WhaffMT0A/d4a9oWTcs4MbFJcUwMlWhiY/GdPJngjZsSTtDTw65fP3J074biEzam3++1nxGnRxCridXF+pbuS7fi3eR3iBDis5LSz+80XY0YGxvL+fPnOXz4MMeOHdO9VlQoFNjY2FCjRg1q1qxJ1apVMTKS1UZvSLIlPmdnHwXQb+MVAsITkp/CuS1Y1K4sbjnjixcH7/yTV2PGoH1r3pb9yJHYtmv7SX1e8btCv6P9CI4J1sWaFGjC6EqjMdAz+LQbEUKINJYpydb/u379OocPH+bIkSO614oKhQJjY2MqV65MrVq1qFWrFlZWVuk1hCxBki3xufMNjab3+sv8+yRhlaC5kT7TW5SgXrH4Ta6j79zhed9+xD1/jlHRIuTdsAE9Y+NP7vNp6FN6HemVaKL8Vw5fMfObmVgaWn76zQghRBpJt9IPH6NEiRIMGjSIffv2ceDAAYYOHUrp0qWJiYnhyJEjjBgxgrVr16bnEIQQacDe0piN3SvStaqbLhYeo6LnustM3HebOLUG4yJFcNu+DcuGDXGeMydViRZAHss8rKu/jjK5yuhiF15eoMNfHfAJ90lV20IIkZEyrFSzq6srXbp0YcOGDZw+fZoJEyZQo0YNTNJwGxAhRPoxUOoxsmFRFrUtg7lRwuYTy0550WbZeXxDo1FaWeE0/XcMXVySXB95+Qra2Ngk8fexNrZm6bdLqe9WXxd7FPKINvvacN3/+qffjBBCZKBM2RfD1taW5s2bs3DhQn788cfMGIIQ4hPVL+7A7j5VKJzbQhe76B1Eg7mnOPsoINlrom7e4mmnTjzp0JE4X9+P6s9IacSUalPoVrybLhYYHUiXA1046H3w025CCCEykGxCJoT4aPnszNnZqwpNyzjpYgHhsbRbfoEFxx4mqjqvCgrCp18/tLGxRF29Gl8e4sI/H9WfnkKPfmX6Mb7KePT14p+qxahjGHxiMH/c+APZdUwI8Tn76I2oL168mCYdywbWQmRtJoZKZrQoSfm8tvy2+xaxKg0aLfx+4B6XnwQx07MUVqYG6JmaYlq5EiHb4rf1Ub9+zdMuXcg1aCC2Xbp8VDmHJgWa4GTuRP9j/QmLDQNg9uXZPA17ysiKI2WlohDis/TRyVb79u1TXetGoVBw+/btVLUhhMh8CoWC1hXyUMzRip/WX+J5UBQAR+760WDeKRa1LUtxZyscJ0zApGRJfMdPiJ+3pVbj9/t0oq5ew2HyJJTm5inus3zu8qz/bj29DvfiefhzAHY82IFPmA8za8hKRSHE5+eTXiNqtdpUfwkhso/izlbs61uNWoUTV51vtuisruq8TYsWuK5fj4Gjo+6csEOH8G7egpgHDz6qPzcrN9Y3WE8pu1K62IVXF2j/V3uehz1P9f0IIURa+ug6Wz4+H7/kOjAwkIULF3L8+HG0Wi0KhYI7d+58dDvZldTZEtmFRqNl0YlHzDj4f1XnSzsx4YdimBrqowoK4sXPQ4k4fVp3XGFigsP48Vg1bPBR/cWoYxh1ZhR/e/2ti9ka2zKnxhxK5SqV2tsRQoj3Sunn90e/RnRycvrwSf+JiopixYoVrFy5koiICLRaLfny5WPQoEEf260QIgvQ01PQu0YBSuexTlR1fscVH26+CGFRu7Lkt7PBZcliAhYsJGDhQgC0UVG8GDIE9esAbDt2THF/RkojplabSh6LPCy5vgSIX6n444EfmVB1QqKSEUIIkVnSZTWiWq1m/fr11KlTh/nz5xMeHo69vT0TJkxg79691K5dOz26FUJ8Jirnz8m+ftUon9dGF7vvG07jeafZe/0FCqUSu359cVmyGD3L+DlWSisrLD7hd4NCoaBP6T5MrDpRt1IxVhPL0JNDWXxtsUxbEEJkujRPtvbt20f9+vWZMGECAQEBWFpaMmTIEA4cOEDz5s3R+4SNaYUQWY+9pTEbulWkR/V8ulhErJo+G64w5r/Vi+Zff43bju0YFyuG44wZGHzEk/P/1zh/Y5bWWYqVUcL2XwuuLmD46eHEqj+umKoQQqSlNMt8Tp8+TdOmTRkyZAhPnz7FyMiIrl27cujQIbp27SobTwvxBTJQ6vHrd0VY0r4sFsYJsxZWnfXGc8k5fIKjMHR2Ju+WzZhXrZLk+jhfX7QaTYr7e7NS0dXSVRfb+3gv3Q52Iyg66D1XCiFE+kl1snXjxg06depEt27duH37Nnp6erRo0YIDBw4wZMgQLC1lGbYQX7q6HrnZ27cqRR0Sfh9cfRZMw7mnOHHfH0UyT7xVQUF4t2rN8z59UYeGprgvV0tX1n+3nvK5E2r5Xfa7TJt9bXgc8jh1NyKEEJ/gk5OtJ0+e0L9/fzw9PTl//jxarZY6deqwZ88exo8fj729fVqOUwiRxbnmMGNHr8q0rpCwb2JQZBydVv7DzEP3Ub+1fFGrVvNi8BBUL18SfvQoXs1bEH33bor7sjKyYkntJfxQ4Add7Hn4c9rta8e5F+fS5oaEECKFPjrZ8vf357fffqNBgwYcOHAArVZL+fLl2bx5M/PmzSNfvnwfbkQI8UUyNlAyuWkJprcoibFB/K8frRbmHnlAxxX/EBAeEx+LjgaDhNeOcU+f4t2yFcE7/0xxXwZKA8ZWHsuAMgN0sbC4MH46/BPb7m9Lk/sRQoiU+Ohkq06dOmzZsgWVSkXhwoVZunQpa9eupWTJkukxPiFENtS8rDM7e1XBLaeZLnb6YQAN5p7iX+9A9MzMcFm0iJz9+sJ/O1ZoY2J4+euvvBz9G5qYmBT1o1Ao+LH4j8z6ZhbGSmMA1Fo1Y8+NZfrF6ag16rS/OSGE+D8fXWcrOjoahUKBQqEgNDSUsWPHfnSnCoWCw4cPf/R1Qojso4iDJbv7VOGX7df568YrAHxDY2i59DzD6hWmazU37Hr1wqREyfgaXMHBAARv2UL0rVs4zZmDoXPKVi/Wdq2Ng5kDfY/2xT/KH4DVt1fzJPQJU6pPwczA7AMtCCHEp0vVdj0vXrzAx8fnk76EEMLC2IAFbcrwW6Oi6OvFP8FSa7RM/OsOPdddIiQqDvOqVeLLQ5Qoobsu+tYtvJo1I/zEiRT35ZHTgw0NNuBu466LHX9+nA5/d+BF+Iu0uykhhPg/H/1kq0+fPukxjo/m7+/PmTNnuHnzJjdu3ODOnTvExMRQoUIF1q5dm6q2z58/z8qVK7l27RqRkZE4OjpSr149unfvjqmpaRrdgRAC4p90d67iRkkXa/qsv8yLkGgADtzy5c7L0yxsW4ZiTo64rluL35QpBG3YCIAmJIRnPXqSa+hQcnTpnKK+cpvlZk39Nfxy6heOPzsOwP2g+7Te15q5NedS0k6mQwgh0l6WTbb27dvH5MmT07zdtWvXMnHiRLRaLblz58bBwYGHDx+yaNEiDh48yIYNG7C2tk7zfoX40pXJY8PeftUYsPkqJ+/Hv+p7GhhJ00VnGdPIg9YVXMg9ejQmpUvzcvRvaKOiQE8P46JFP6ofUwNTZn8zmzmX57Dy1kogfoufLvu7ML7KeL7L912a35sQ4suWZcu5m5ubU7lyZXr06MH8+fPp1atXqtu8efMmkyZNAmDcuHEcP36cnTt3cvjwYTw8PHj06BGjRo1KdT9CiOTZmhmyqlN5Btcp9GZePLEqDcN33mDwlmtExqqwatSIvJs3YZg3L3b9+2NW8auP7kepp2RQuUGMqzwu0RY/v5z6hQVXF6DRpryQqhBCfMhHP9n6XDRv3pzmzZvrvvf19U11mwsXLkSj0dCkSRNatmypi9vb2zNz5kzq16/PwYMHuXv3LoULF051f0KIpPT0FPStVZDSeWzov+kKryMSb2a9sG0ZChQqhNv2bShMTJJcr4mMRGFiguJNtvYePxT8AWcLZwYeH0hITAgAi68txivEi/FVxmOin7R9IYT4WGmabKnVary9vQkNDUWlUr333PLly7/3eEaLiIjg1KlTAHh6eiY5njdvXipWrMjZs2fZv3+/JFtCpLOqBeM3s+678TIXveO32rnvG07j+WeY3LQ435dKuhJRq1LxrEdP9O3scBg/Dj2zD68yLJ+7PBu/20jvo73xCvEC4ID3AXzCfJhbcy52pnZpe2NCiC9OmiRbr169YubMmRw6dIjo6OgPnq9QKLh9+3ZadJ1m7ty5Q2xsLIaGhpR4a9XT28qWLcvZs2e5du1aBo9OiC9Tbqv4zaynH7jHkpPxW+1Exqrpv+kq/3gFMqphUYwNlLrz/efMIfLiRQCi793Dee4cjPLn/2A/LpYurPtuHT+f+JmzL84CcPP1TVrta8W8mvMomuPj5oUJIcTbUp1sPXv2jNatW/P69Wu0Wu2HL4AUn5eRvLzi/6J1dHTEwMAg2XPy5MmT6Nz32bRpE1u2bElR348ePUrhKIX48rzZzLpcXlsGb7lKaHT8U/P1F55y7XkwC9uUJU8OU7QqFVE3b+qui330CK8WnjiMG4dVwwYf7MfS0JIFtRYw9Z+pbLq3CQC/SD86/t2RSdUmUce1TvrcoBAi20t1sjVz5kwCAgKwtbVl8ODBVK1alZw5c6JUKj988WckJCR+voaVldU7z3lz7M257+Pv78+tW7fSZnBCCOoUtWdfv2r03nCZ68/jfwZv+oTSYN4pprcoSV2P3ORZtgz/ufN4vXQpANrISF4MGULU5cvkGvYLeoaG7+1DX0+fERVHkM86H1P/mYpaqyZaHc2g44PoXao3PUr0SNFcMCGEeFuqk62zZ8+iUCiYOXMmFStWTIsxZYqY/7b/eNdTLQDD/35Rx6RgqxA7Ozs8PDxS1PejR49S9PpViC+di60pW3tWYsLeO6w9/wSAsGgVPdZeomtVN36pX5hcgwZiUqoUL375BU1YGABBGzYQdfMmzrNnYeDo+MF+WhdujaulK0OODyEsLr6NBVcX8Dj4MeOqjMNY3zj9blIIke2kOtmKjY3F2Ng4SydaAEZGRgDExcW985zY2NhE575Pq1ataNWqVYr6btq0qTwFEyKFjPSVjG9SjPJutvy6/ToRsfH7Gy4/7cWVZ8HMa10ax5o1cNuxnef9+xNz+w4A0dev49W0GY6//455taof7KeyY2XWN1hP36N9eRIan9j97f03z8KeMafmHHKZ5kq/mxRCZCuprrPl7Oz8Wc7B+lgpeUWYkleNQoiM0bikI7v7VsXd3kIXu/QkiAZzT3H8nh+GLi7k3bgR6xYtdMfVwcE8694d/4ULU9SHm5Ub679bz1cOCbW8br6+Set9rbn9+vNa5COE+HylOtmqX78+MTExnDt3Li3Gk2ny5s0LwIsXL975dOvp06eJzhVCZK78dub82bsKzcs662JBkXF0WnmR6QfuodE3wGH8OBwmT0Zh/N+rP60WPZOUb7tlZWTFotqLaOmeUHvvzcT5g94H0+xehBDZV6qTrS5dulC4cGFGjRrFs2fP0mJMmaJIkSIYGBgQGxvL9evXkz3n0qVLAJQqVSoDRyaEeB8TQyXTW5RkWvMSGBsk/Eqbf+wh7f64gF9oNNY/NImvOu/qikWdOth26vhRfRjoGTCy4khGfDUCpSJ+8U+0OprBJwaz6NqibPF0XwiRflI9Z8vY2JiVK1cyatQoGjduTN26dSlevDhmHygm2KRJk9R2nabMzc2pWrUqx44dY8uWLZQtWzbRcW9vb86fPw9AvXr1MmOIQoj38CznQglnK3qtv8xj/wgAzj8O5Lu5p5nbqhSV3d3Ju30bQJIVhVqVCpTKD640bFW4Fa6Wrgw+MZiw2PiJ8wuvLuRR8COpOC+EeKc0KWrq4+NDQEAAUVFR7Nq1i127dr33fIVCkWnJVuvWrfH19aVDhw506tQp0bFevXpx/Phxdu3aRZkyZfD09EShUODn58egQYPQaDTUrl1bqscL8ZkqnNuS3X2qMnzHDXZfewFAQHgM7f64wIDahehdowBKvaQJld/v04l7+RKHiRNQWlgkOf62So6V2PDdBvoe7Yt3qDcQX3H+aehT5tacS26z3Gl+X0KIrC3Vydbdu3dp3749UVFRQHzpBBsbm3Svs/Xy5ctECdublYKXL1/mq68SJrN27dqVbt266b739fXFx8eHsP+WhL+tRIkSDBs2jClTpjB69GgWLVqEjY0NDx8+JDY2Fjc3N8aPH59+NyWESDVzI33mtCrFV/lsGbv7NrFqDRotzDx0n4vegcxqWYqc5gkrikMPHCRw9WoAou/dxXnOHIw/8AdVXqu8uorz517Gz1e9E3iHVntbMbvGbErlKpVu9yeEyHpSnWzNmzePyMhIXFxcGD9+PBUqVEBPL9VTwT5IrVYTHBycJK5SqRLFP7Z+VadOnXB3d2fFihVcv36d169f4+joSL169ejevfsHX48KITKfQqGg7VeulHS2ptf6yzwNjATg1IMAGsw9xbzWZajgZgtA6F9/6a6Le/IU75atyD16FNbNmr23DysjKxbWXsiMf2ew7s46AF5Hv6bLgS78Vuk3vi/wfTrdnRAiq1FoUzmzs3LlygQFBbF161aKFSuWVuP6oryps+Xh4cGOHTsyezhCZCuh0XEM3Xqd/bde6WJKPQWD6hTip6/zo9BqCFi0mIAFC+CtX4dWTZuSe9RI9Ew+PA9rx4MdjD8/HpVGpYt1LNqRgWUHotTLWrtpCCFSLqWf36l+BBUVFYWJiYkkWkKIz5KlsQGL2pVhdMOiGCjj52upNVp+P3CPLqsvEhStxq5Pb1yWL0NpY6O7LmTHDrxbtiImBXuhNi3YlD++/QNbY1tdbPXt1fQ+2ls3kV4I8eVKdbLl6uqKSqVCrVanxXiEECLNKRQKulR1Y2vPyjhZJzypOn7PnwZzT/GvdyDmVargtnMHJqVL647H3L+Pd/MWhO7f/8E+ytiXYWODjbjbuOtiZ3zO0GZfG10FeiHElynVyVaTJk2IjY3l6NGjaTEeIYRIN6VcrNnXryq1i9jrYi9Domm59DyLTzxCmcse1zWrsX1rpbImIgKfAQN5NXES2g/8Uelo7sia+muolaeWLuYd6k3rfa05++Jsmt+PECJrSHWy1b59eypWrMjo0aO5cuVKWoxJCCHSjbWpIcs6lGVkgyLo6yW8Vpzy9126rvmX4Fgt9sN+wWnuHPTMzXXXqYOCIAWLf0wNTJn5zUx6luypi4XFhvHT4Z9Yc2uNFEAV4guU6tWIixYtolSpUty+fZs2bdpQtmzZFBU17dOnT2q7FkKIT6JQKOhaLR+l89jQd8NlXoTEr1o+etePBnNPMb9tGcp8+y3G7u48HzAQbUwMDmPHfLDo6Rt6Cj16l+pNAesCjDw9kmh1NBqtht///Z17QfcYXWk0RsoPb2gvhMgeUr0asXDhwrpfQG+aSskvpDt37qSm22xFViMKkXmCI2MZvOUaR+766WL6egp+qVeYrtXc0MbEoH79GgMnpyTXajUaFB942nXn9R36HevHq4iE1ZAlcpZgdo3Z2Jnapd2NCCEyXEo/v1P9ZKt8+fKpbUIIITJN/GvFciw79ZhpB+6h1mhRabRM/OsOF7xeM71FSayTSbRCdu8meMtWHGfMwMA+1zvbL5KjCJsabGLQ8UFc9rsMwPWA67oCqMXtiqfbvQkhPg+pfrIlUk+ebAnxefjXO5C+G6/wMiShGLKjlTHz2pShrGtCWYiYhw/xauGJNioKZY4cOE3/HbNKld7bdpw6jkn/TGLb/W26mKGeIb9V/o3G+Run/c0IIdJdhtXZEkKI7KJcXlv+6leNmoUTnlS9CImm5ZJzLDnxCI0m/m/TkL170f63RZn69WuedvkR/4UL0Wo072zbQGnAb5V+Y1TFUegr4l8qxGpiGXF6BL9f/D1RQVQhRPYiyZYQQrzFxsyQ5R3KMfy7wrpNq1UaLZP/W60YFBGLXf/+5B4zBoWhYfxFWi0Bc+fxrFt3VIGB723f092Tpd8uxcYo4UnZmttr6H2kNyExIel2X0KIzCPJlhBC/B89PQXdq+dnS4+KOFoZ6+JH7/rx3dxTXH4ahE2rlrhu3ICBi4vueMSZM3j90JTIS5fe23753OXZ2DBxAdSzL87Sel9rHgY9TPsbEkJkqo9KtpYvX/7RGzt/yI0bNzhx4kSatimEEGmhrKst+/pVo9ZbrxVfhkTjuSS+CKpRkaK47diORZ06uuMqX1+edOjI6+XL3/ta0cnciTX11/Ct67e62LOwZ7T9qy1Hnh5JnxsSQmSKj0q2pk+fTu3atVm1ahWhoaGp6vjff/+lR48eeHp6cuPGjVS1JYQQ6cXGzJDlHcsx4rukRVC7rL5IsMIQp7lzsB/+KxgYxF+kVuM3fQbPf+qF+j2/K00NTJn+9XT6lu6Lgvi2I1WRDDg2gAVXF6DRvjtZE0JkHR+VbPXo0YPw8HCmTp1K1apV6devHwcOHOD169cfvDYuLo7r168ze/ZsateuTfv27Tlx4gTFixendu3an3wDQgiR3hQKBd2q52NLz0pJ9lb8bu4p/vEKxLZDB/KuW4u+o4PueJyvb8K8rve03b1Ed+bWnIu5QULF+sXXFtP/WH/CY8PT/oaEEBnqo0s/+Pr6MnPmTPbu3YtardYVMHVwcMDd3R1bW1usrKwwNDQkNDSUkJAQnj17xt27d4mLiwPii5/myZOH/v3706BBg7S/qyxGSj8IkXUER8YyZOt1Dt/x1cX0FDCwdiF61SgAoSG8GPYrkRcv4rZ9G4Z586a47cchj+l/tD/eod66mJuVG3NrzCWvVcrbEUJkjJR+fn9ynS1fX1+2bNnC9u3befUqoTJyctXj33Shr6/P119/TcuWLalWrVqKt77I7iTZEiJr0Wq1rDzjzeS/7xCnTvgVWqVADma1LIWdmSGx3t4Y5cuX9Fq1GoVS+c62w2LD+PXUr5x4njCX1cLAginVp1DduXra3ogQIlXSPdl62/3797l48SLXr1/Hz8+PwMBAYmJisLa2xtbWlgIFClCuXDnKlCmD+Vsbu4p4kmwJkTVdfx5Mnw1XeBoYqYvlNDdkdsvSVC2YM8n50bdv4zPkZxynTMakRIl3tqvRalhwdQFLry/VxRQo6FemHz8W+1H+UBXiM5GhyZZIHUm2hMi6QqPjGLb9On/dePsJP/T+pgADahdEXxk/NVYdFoZXs+bEPX0KBgbY/zwEm/bt35s4HXpyiBGnRxClitLF6rjWYUKVCZgamKbfTQkhUkQqyAshRAawNDZgQZsyTGhSDEP9+F+pWi3MP/aQNssu8DIkPlGK/Ocf4nx84i+Ki8N30mR8+vV772rFOq51WP/depzNnXWxQ08O0favtjwNfZp+NyWESFOSbAkhRCopFAraVXTlz15VyJfTTBf/xzuQ7+ac4sgdXyxq1cJ19Sr0cyXU7Ao7dBivps2Iek/5m4I2BdnUcBOVHSvrYg+DH9JqbytOPj+ZPjckhEhTkmwJIUQaKepoye6+VWlSylEXC4qM48fV/zJuz230S5XBbecOzKpU0R2Pe/4c7zZtCVyzlnfN6rAysmJhrYX8WOxHXSwsLow+R/qw9PpSqcclxGdOki0hhEhD5kb6zGpZimnNS2BikLDqcMUZL5otOosPxrgsW4rdgP6g99+v4Lg4fCdNeu9rRaWekgFlBzD96+mY6MfX+tKiZd6VeQw6PoiIuIh0vzchxKeRZEsIIdKYQqHAs5wLe/pWwd3eQhe/4RNCg7mn2XPjFTl79iTPqpXo29npjr95rRh97947266bty7rv1uPi0XCnoxHnh6hzb42eIV4pc8NCSFSRZItIYRIJwVyWbCrTxXafJVHFwuPUdFv4xV+3XEdvVJlcftzJ2aVE+ZjacLCUFpYJNecTkGbgmxssJGqTlV1scchj2mzrw3Hnh5L+xsRQqSKJFtCCJGOjA2UTPqhOAvalMHCSF8X3/jPM75fcJrHKkNcli/Drn8/UCpxmDIZA0fH97QYz8rIivk159OteDddLDwunH7H+rHw6kKZxyXEZyTVydaECRMICwtLi7EIIUS21aCEA/v6VaOks5Uudt83nMbzT7P53+fk6NmT/Af2Y1GjRpJrNTExybap1FPSr0w/Zn0zC1P9hLpbi64tou/RvoTEhKT9jQghPlqqk61169ZRr149tm/fnhbjEUKIbCtPDlO29qxMt2puulh0nIZhO27QZ+MVonPaJ7lGFRTE4wYNeb1i5TtXK9Z2rc2GBhvIa5lXFzv5/CSt9rbiXuC7538JITJGqpOtYsWK8fr1a0aOHEnLli25efNmWoxLCCGyJUN9PUY0KMqKTuWwMTXQxfddf8l3c05x5WmQLqbVaHg57Ffinj/Hb9o0nv/UC1VQUHLNkt86PxsabOAbl290sefhz2n3Vzv2Pt6bbvcjhPiwVCdbW7duZdy4cVhbW3Pt2jU8PT0ZPXo0Qe/4hSCEEAJqFrbn7/7VqZjPVhd7HhRFi8XnWHT8ERqNFpWvL9F37+qOhx8/jlfTZkRevpxsmxaGFsypMYe+pfuiIH4boGh1NL+e+pUp/0whThOXvjclhEhWqpMthUKBp6cnBw4coHXr1igUCrZu3Uq9evXYsGHDOx97CyHEly63lTHru1ZkcJ1C6P23RaJKo2Xq/rt0XPkPwea28asVv66uu0b18iVP2ncgYOkytJqkk+D1FHp0L9GdRbUXYWloqYuvv7Oerge64h/pn+73JYRILM1WI1paWvLbb7+xfft2SpcuTUhICOPHj6dZs2ZcfsdfYUII8aVT6inoW6sgm3tUwtHKWBc/9SCA+nNOcdZfhcuiReT6+WfQ/281o1qN/8yZPOvWHdXr18m2W8WpCpsbbqawbWFd7LLfZVrubclVv6vpeUtCiP+T5qUfChcuzIYNG5gyZQo5cuTg9u3btG3bll9++YWAgIC07k4IIbKF8nlt+at/Nep6JEySDwiPocOKf5hy4B6WnTqRd93aRGUhIs6c4XGTJkScP59sm84Wzqypv4bG+RvrYv5R/nTe35kNd+TNgxAZJd3qbDVp0oQDBw7QqVMnlEolu3fvpl69eqxatQq1Wp1e3QohRJZlbWrI4nZlGd+kGIb6Cb+el5x4TPPF5/BzKYjbzh1Y1KmtO6b2D+Bp5y68XrUq2TZN9E2YUGUCI74agb4i/smYSqti8j+T+fX0r0TGRabrPQkh0rmoqZmZGcOGDWPHjh3ky5eP8PBwpk6dSsOGDZkyZQq7d+/m4cOH8teVEEL8R6FQ0L6iK7t6V6FALnNd/NqzYBrMPc1er3Cc5s7FfuRIFAYJqxmNCxd5b5utCrdiZb2V2JkkbA+07/E+2v7VFu8Q73S5FyFEPP0Pn5Jy0dHReHl58ejRIx4/fqz7v0+fPiU2NhaFQoFWq8XLywtvb2/ddcbGxri7u7Np06a0HI4QQmRZRRws2d2nCuP23GbTxWdA/FY//Tdd5eT9AMa2aIlr6VL4DBqEVcNGmFX86oNtlspVii2NtjDkxBAu+V4C4GHwQ1rva82EKhOo5VorXe9JiC9VqpOt0NBQpk+fzunTp3n16lWip1T//8TKzs6OfPny4erqSkBAAHfu3OHly5dERUVx7dq11A5FCCGyFVNDfaY0K0HVgjn5dccNwqJVAGy//JxLTwKZ17oMRbfvQM/EOMm1qqAglGZmKAwNE8VzmuRk2bfLmHt5LqturQLit/kZcHwAXYp1oW/pvujrpenf4UJ88VL9EzV8+HCOHDmiS6wMDQ1xdXXFzc0NNzc38uXLR758+XBzc8Pc3DzJ9UFBQdy+fZvbt2+ndihCCJEtNSzhSElna/pvusLlp8EAeL+OpOmiMwytW5gfq7r9V1UrnlalwqdvPzSxsTjNmI6hi0ui9gz0DBhcbjAl7Eow6swoIuIiAFhxcwU3A24yrfo0cpjkyKC7EyL7S3Wydf6/VTBdu3alefPmuLq6olAoPnBVAhsbG6pUqUKVKlVSOxQhhMi2XGxN2dKjEnOPPGDesYdotRCn1jLxrzucfODPDM+S5LKIf8IVsGgxkf/+C4DXD01xGD8Oy/r1k7RZx7UOBawLMPDYQB6FPALgn1f/4LnXkxlfz6BUrlIZdn9CZGepniBvaWmJpaUlQ4YMIW/evB+VaAkhhEg5faUeg751Z0PXiuS2TFyT67s5pzh+zw+tVovqdUKZHU14OD4DB/Fy9G9ooqKStOlm5caGBhuonzchGfOL9KPz/s6sv7NeFjAJkQZSnWwdPXqUo0ePpsVYhBBCpECl/Dn4u3816hR9uyZXLJ1WXmT83jvYjhyF08wZ6JmZ6Y4Hb9mCt6cnMQ8eJGnP1MCUqdWnMqzCsETlIab8M4VfTv4i5SGESKU0Kf1g9tYPtBBCiPRnY2bI0vZlGf+9R6KaXCvOeNFkwVn8ylXDbecOjIsV0x2LefAQrxaeBG3ZkuSJlUKhoG2Rtqyst5JcJrl08b+9/6bVvlY8Cn6U/jclRDaVrnW2hBBCpB+FQkH7SnnZ06cq7vYWuvidl6E0nHearS+1uK5fh23nzrpj2uhoXo3+DZ9Bg1CHhSVps1SuUmxutJmvcieUkvAK8aL1vtbse7wvfW9IiGxKki0hhMji3HNbsKtPFTpVzquLRcdpGLHzJj9tuYFhnwG4LFmM0sZGdzzs7/14e7ZEGxubpL2cJjlZUmcJ3Yp308WiVFEMOzWMCecnEKtOeo0Q4t0k2RJCiGzA2EDJmMYe/NGxHLZmCbW1Dtzypf6cU1x39sDtzz8x/SrhiZV1S88kdbjeUOop6VemHwtqLcDS0FIX33xvMx3+7oBPuE/63YwQ2YwkW0IIkY3UKmLP/v7VqFYwpy72KjSatssvMOtKIA7LlmHXvx/mNWti27HjB9ur7lydLY224JHDQxe79foWnns8Ofn8ZLrcgxDZjSRbQgiRzeSyNGZ15wqMbFAEA2V8OR6tFhYef0TzpRcIa9EB5/nzkpTq0cTGEvHPP0naczJ3Yk39NbR0b6mLhcaG0vtIb+ZcnoNKo0rfGxIii5NkSwghsiE9PQVdq+VjZ68q5LNLWDF+/XkIDeaeYsul50lWJPrPmMHTDh3xnfZ7krlchkpDRlYcyZRqUzDRN9HFl99YTreD3fCL9EvfGxIiC5NkSwghsrFiTlbs7VuV1hUStuyJjFXzy/Yb/LTuMkER8UlV2PHjBK5eA0DgihV4t25DrLd3kvYa5GvAxgYbcbNy08X+9f2XFntacP7l+fS9GSGyKEm2hBAimzM11Gdy0xIsblcWa1MDXXz/rVfUm3OSMw8DMHB0xLBAft2x6Fu38GrajJBdu5K0l986P5sabOI7t+90scDoQLof7M6iq4tQa9Tpe0NCZDGSbAkhxBeiXrHcHBhQnaoFEibP+4bG0O6PC8x6qMZx42asWybMy9JERvLil2H4/DwUdXh4orZMDUyZUm0KoyqOwlAvfkWjFi0Lry2k5+GeBEQFIISIJ8mWEEJ8QewtjVnTJX7yvKEy/iNAq4UlJx/TfOVlInoPxmnOHPQsE8o9hO7Zg9cPTYm6fj1RWwqFAk93T9Z9tw4Xi4TXlOdfnsdzjycXX13MmJsS4jMnyZYQQnxhdJPne1emQC5zXfzWi/jK87ssC+G2cwcm5crqjsU9e4Z3m7YELF2GVqNJ1F6RHEXY3HAzdVzr6GL+Uf50PdiV5TeWo9EmPl+IL40kW0II8YXycLRiT5+qtK/oqotFx2kYtesWPQ76YDpvCTn79AG9/z4qVCr8Z84k4tSpJG1ZGFow4+sZ8ZtZ68VvZq3RaphzeQ69j/QmKDooQ+5JiM+RJFtCCPEFMzFUMr5JMf7oWI4cb1WeP3bPn3rzz3KlRjNc165B39EBAMsGDTCrXj3Ztt5sZr22/loczRx18dM+p2m+pzmXfC+l780I8ZmSZEsIIUR85fkB1alZOJcuFhgRS/e1lxj/xAD7zduwaduW3GN+S1IM9f8Vy1mMLY228I3LN7qYX6QfXQ50Yen1pfJaUXxxJNkSQggBgJ2FEX90LMeEJsUwNkj4eNj4zzMar7rGy469UVpYJLnOf+48ou/dSxSzMrJibo25DCk3BH1FwmvFeVfm0fOQrFYUXxZJtoQQQugoFAraVXRlX79qFHey0sW9X0fSfPE55hx+gEqd8GQqZM8eAhYuxLuFJ4Fr1iaqSq9QKOjo0ZE19dfgZO6ki597eY4We1pw4eWFjLkpITKZJFtCCCGSyG9nzo5elelTowB6/701VGu0zDp8nxZLzuEdEIE6NJRX48YDoI2NxXfSJJ716IEqIPFTq+J2xdnSaAu189TWxQKiAuh2sBsLry6UIqgi25NkSwghRLIMlHoMqevO5h6VcLZJ2A/xytNgvpt7is13gnGcORNlzoQiqREnT/H4+yaEnziRqC1LQ0tmfjOTXyv8ioFefBV7LVoWXVtEt0Pd8I/0z5ibEiITSLIlhBDivcrnteXv/tVoVsZZF4uMVTN85w36PzLCYv0WzL/+WndM/fo1z3r05NWEiWhiYnRxhUJBmyJtkhRBvfjqIs33NOeMz5mMuSEhMpgkW0IIIT7IwtiAGZ4lWdS2DDZv7a949K4f362+wa2+v2E/aiQKIyPdsaB16/Bu3oLoe/cTtVU0R1G2NNxCvbz1dLHA6EB6Hu7JzH9nEqeOS/8bEiIDSbIlhBAixeoXd+DAgOp8426niwVFxtFz3WUmGRbDbsMmjAoV0h2LefAA7xYtCFy7LlE75obmTKs+jd8q/YaRMiFBW3lrJR33d+RZ2LP0vxkhMkiWT7bOnz9Pjx49qFixIiVKlKBevXrMnj2byMjIj25r2LBhuLu7v/fr5MmT6XAXQgiRdeSyNGZlp/JMaFIMEwOlLr7t0nMa73mB3+9LsOnQXhfXxsYS++RJknYUCgXNCzVnQ4MN5LPKp4vfCLiB5x5P9nvtT98bESKDZOlka+3atXTq1Injx49jZGRE/vz58fHxYdGiRTRv3pzg4OBPatfBwYEyZcok+2VlZfXhBoQQIptLKBFRlZIu1rq4T3AUrVddYmWpH7BfvARlzpwYFSxAriGD39lWIZtCbGq4iWYFm+li4XHh/HzyZ8acHUOUKio9b0WIdKef2QP4VDdv3mTSpEkAjBs3Dk9PTxQKBb6+vvz000/cunWLUaNGMW/evI9uu1mzZvTt2zethyyEENlOPjtztvesxMLjj5hz5AFqjRatFpacfMyJ3BZMW7YBd2MVesbGia7TajRoY2LQM4lf5Wiib8KYymOo6FCRsefGEh4XDsD2B9u54neF37/+nUI2hZL0L0RWkGWfbC1cuBCNRsP3339Py5YtddtH2NvbM3PmTPT09Dh48CB3797N5JEKIUT2pq/Uo1+tguz4qTL57Mx08buvwmi27gbLvFSJCqFC/OR5r6bNiLp5K1G8nls9tjTaQvGcxXWxxyGPabOvDVvubUlUNFWIrCJLJlsRERGc+m/XeU9PzyTH8+bNS8WKFQHYv1/e+QshREYo6WLNvr7V6FjJVReLU2uZfvA+zRef47F//NOq6Hv38Zs+g1gvL7xbtSJgyVK06oTCpi4WLqyut5rOHp11sRh1DOPPj2fg8YGExIRk3E0JkQayZLJ1584dYmNjMTQ0pESJEsmeU7ZsWQCuXbv20e1fuHCBfv360aFDB/r06cOiRYvw8fFJ1ZiFEOJLYGKoZOz3xVjf9SscrBJeHV59Fl8IddUZL8IvXkQb9195B5UK/1mzeNKxI7HPE37PGigNGFRuEItrL8bW2FYXP/L0CE13N+Xiq4sZdk9CpFaWTLa8vLwAcHR0xMDAINlz8uTJk+jcj3Hx4kUOHDjAhQsXOHToELNnz6Zu3bosW7bs0wcthBBfkCoFcrJ/QHWalknYEzE6TsOYPbfpH50fg0V/oO/goDsW9e8lvL7/nuA//0z0qrCKUxW2NdpGRYeKuphfpB8/HviRuZfnEqeRmlzi85clJ8iHhMQ/Qn7fysA3x96cmxKurq4MGzaMihUr4uTkhKGhIffu3WPFihXs37+f6dOnY2pqStu2bT/Y1qZNm9iyZUuK+n306FGKxyiEEFmFlYkBMz1LUdcjN8N33OB1RCwAZx6+5odn+owav4TKO5cQtm8fAJqICF4O+5Xw4ydwGPMbSmtrAOxM7VhSZwmrb61m7uW5qLQqtGhZdmMZF15dYEq1KYkq0gvxucmSyVbMf9s/vOupFoChoWGic1Pip59+ShIrWbIkc+bMYezYsWzYsIHZs2fTpEkTzMzMkmkhgb+/P7du3XrvOUII8SWo65Gbsq42jNh5gwO3fAEIi1ExdO99apdoxa+Vv0E1ZRyasLD4Y/v3E3X5Mg6TJ2FepQoAego9OhfrTIXcFfjl1C88CY2v23Xd/zot9rRgVMVRNMjXIHNuUIgPyJKvEY3+2w4iLu7dj49jY2MTnZtagwYNwsDAgNDQUM6fP//B8+3s7PDw8EjRl/H/LYkWQojsJqe5EYvblWVGi5JYGCX8nX/4ji/Nbhly5/eVmJYvr4ur/Px49mNXwk8n3i/RI6cHWxpuoUmBJrpYRFwEw04NY/ip4YTHhqf7vQjxsbLkk62UvCJMyavGj2FhYUHBggW5ffs2T5KphPz/WrVqRatWrVLUdtOmTeUpmBAi21MoFDQr60yl/DkYuu06px8GABAcGceA/d58V7c/g6peI27+bIiLw6R0acwqVUzSjqmBKeOrjKeyY2XGnxtPWFz8E7E9j/dwxe8KU6tPpYRd8ounhMgMWfLJVt68eQF48eLFO59uPX36NNG5aeHNa0uVSpVmbQohxJfG0dqEtT9WYHyTYpgaJmz389fNV7Tyd+Hh9BWYlCqF47SpKJTKd7ZT360+WxtvpZRdKV3sefhzOvzdgcXXFqPSyO9q8XnIkslWkSJFMDAwIDY2luvXryd7zqVLlwAoVapUmvSpUql4/PgxALlz506TNoUQ4kulUChoX9GV/f2rU8EtobRDQHgsfY/5MqvJUCJz2Ce5LmjLlkT7LDqZO7Gy3kp+KvkTeor4jzS1Vs2CqwvocqALz8Oep//NCPEBWTLZMjc3p2rVqgDJrvjz9vbWzauqV69emvS5efNmwsLC0NfX1xVMFUIIkTp5cpiyqVtFRjcsirFBwkfSrqsv+HbWSY7e9dXFIi9d4tWYsTxu8gNBmzbrSkTo6+nTq1QvVtZdiZN5QqmJK35XaL6nOXse7ZHK8yJTZclkC6BXr14oFAp27drF5s0JP3R+fn4MGjQIjUZD7dq1KVy4cKLratasSc2aNZNUlj9z5gy///473t7eieKxsbGsXbuWyZMnA/FzsXLlypV+NyaEEF8YPT0FXaq68Ve/apTJY62L+4XF0GXVvwzZeo2g4HBe/DIMNBq0UVG8GjOGZz17Eufnpzu/jH0ZtjbaSsN8DXWxiLgIhp8eztCTQ6XyvMg0Cm0WTvdXrVrFlClT0Gq1ODg4YGNjw8OHD4mNjcXNzY0NGzZga2ub6Bp3d3cAJk+eTNOmTXXxw4cP07t3bwBy5syJvX3842svLy8iIyMBqFu3LtOnT9eVlUgrbybIe3h4sGPHjjRtWwghshK1RsvyU4+Zceg+saqE/RTtLY0YmScO93ljUb+1OEppZUXusWOw/L+3GH97/Z1o8jxAbrPcTKo6ifK5yyNEWkjp53eWfbIF0KlTJ1auXEn16tWJiori4cOHODo60rNnT7Zv354k0XofDw8PevXqReXKlTEyMsLLy4v79+9jYWHBt99+y6JFi5g7d26aJ1pCCCESKPUU9Pg6P/v6VqWEc8Jqct/QGPre1LCw12zUX9fUxdUhIfgMGIjPkJ8TJWH13eqzvfF2ytmX08VeRbzixwM/MuvSLOLUUnleZJws/WQru5AnW0IIkZRKrWHpqcfMPvSAWHXCUy47CyOG24VQZPEktP+9eQDQz5ULh4kTMa9WVRdTa9SsvLWSBVcWoNImrE4sYluEydUmk986f8bcjMiWvognW0IIIbIvfaUevb4pwL5+VSnpYq2L+4fFMPCxMQt7zyW2XMKCJZWfH8+6dePlmDG6ja6Vekq6Fu/KugbryGuZV3funcA7eO7xZO3ttWi0CYmcEOlBki0hhBCftYL2FmzvWYlf6xfGUD/hY2v3g2B+LNyGmz2Ho3hr+zZVQADoJ67Z7ZHDg80NN+NZyFMXi9XEMu3iNLof7M6riFfpfyPiiyXJlhBCiM+evlKPHl/n569+1Sj91orFgPBYfn5ly5ze84j0KIXSxgaHsWNRKBRJ2jA1MGVUpVEsqLWAHMY5dPELry7QdFdT9j3eJyUiRLqQZEsIIUSWUSCXOdt6VmZkgyIYvfWU62/vcLqW7My1UXNRJrM4KvbZM93/ru5cnZ3f76R2ntq6WFhcGMNODePnkz9LiQiR5iTZEkIIkaUo9RR0rZaPv/tXo5yrjS4eFBnHzyd9+XH1v7wMidLFww4f5lG9+vjNno02NhYAG2MbZn4zk4lVJ2JmYKY794D3AZruasoZn8QbYAuRGpJsCSGEyJLy2ZmzpUclxjQqmmiPxaN3/fh25kk2XHhK7OvXvPxtDKjVvF68BC/PlkTfvQvEbxnUOH9jtjfeTln7srrr/aL86Hm4JxPPTyQyLvL/uxXio0myJYQQIsvS01PQqYobBwZUp1rBnLp4WIyK4Ttv0G7tNV5ZO+jiMXfv4tW8Bf4LFuhWLDqZO/HHt38wuOxgDPQSJtpvureJFntacNXvaobdj8ieJNkSQgiR5bnYmrKmSwWmNS+BpXHCSsR/XkXRvXQX/m7RHzX/TZpXqQiYNx+vli2JvncfiC8R0alYJzY13EQhm0K665+GPaXD3x2YeWkmMeqYDL0nkX1IsiWEECJbUCgUeJZz4fCgr/m2qL0uHq3SMDfOhWGdZvI0f3FdPOb2HbyaNydg8WLdU65CNoXY2GAjXYt3RU8R/xGpRcvKmytptbcVt17fytibEtmCJFtCCCGylVyWxixpX5YFbcqQ0zxhi7WbwWp6l+zExqaDiNH77+lXXBz+s+f895TrHgCGSkP6l+nPmvprEhVCfRj8kHb72rHo6iLiNLLdj0g5SbaEEEJkOwqFggYlHDg08Gt+KO2ki6s0WtZoHOnfehq3CiVsSB1z+w7BW7claqOkXUm2NNpCuyLtEq7Xqlh4bSFt97XlYdDD9L8RkS1IsiWEECLbsjEzZFbLUqzuUgFnGxNd/EmEhiFFW7Loh58JMzDBwNmZXAMHJLneRN+EXyr8wh/f/oGjmaMufifwDp57PVlxcwVqjTojbkVkYZJsCSGEyPa+LmTHwYHV6VrVDb23isvv1trzU9MJ3OjzGwpT00TXaDUaXV2uCg4V2PH9DpoVbKY7HqeJY9alWXTY34HHIY8z5D5E1iTJlhBCiC+CqaE+IxsW5c/eVSjiYKmLv47RMvB8CN3WXEpUDDVo40a8mjUn6sYNAMwMzBhTeQwLai3AzsROd951/+u02N1CnnKJd5JkSwghxBelhLM1u/tU4Zd6hRNt+XP4ji+1Z5xg5Rkvop4+w2/GTGIePMC7ZSv8pk9HEx0NJGz30yhfI921sZrY+Kdcf3fgcbA85RKJSbIlhBDii2Og1OOnb/JzYEB1KuVL2JQ6IlbN2D23abryCvcN/tsKSKPh9fI/8Pq+CZEXLwJgZWTFpGqTmFdzXuKnXAHXabGnBX/c+AOVRpWh9yQ+X5JsCSGE+GLlzWnGhm5fMa15CaxNE6rH34lSMqDGQJZVbkukvhEAsU+e8KR9B16NG4c6PByAb1y+Yef3O2mcv7Hu2lhNLLMvz6b9X+15FPwoY29IfJYk2RJCCPFFe1MM9cigr2laJqFMhAbYkas0PzUew1nHhGKoQRs28rhRY8KOHwfin3JNrDqR+TXnJ3rKdfP1TVrsacHyG8ulLtcXTpItIYQQAshhbsRMz1Js6PoVbjnNdHG//7V35/E1Xfv/x18nyckkghDzENMRQwUxlhqrhBqKGkrVLaXa0ttBq1q3v6tVWr1KubRoaaume81tUaraUqU1VkRJCEEQkpBB5v37wzcncpOQnORk4P18PDwesvdaa+/1Odn7fLL32munmXmn9VO82/V5ItzKAJASHs75Z8dz4dVJpERGAtCpRqcsV7mS05KZe3Auw78dTtC1oMLtkBQbSrZERERu82C9Cmx58SEmdquPs2PG1+Qez9qM6/km6yxdSP2/V/nc+OYbrm/cZC2TfpXr393+TUW3itblQZFBDPt2GB8f/FjvWLwPKdkSERH5H65mR17ubuG7Fx+iTW0v6/KbhgOLG/XmxZ6TCfTywaVRQ7yeHJGlfsfqHVnffz0D6g+wLks1Uln852IGbRrEwcsHC6UfUjwo2RIREclBvYoerBrbllmDmlLutgH0IS5evNrxBeYHvEhUYlqmOmlJSRgpKXg6e/LPB//Jou6LqOaRMRYs9EYoT219ium/TScuOa7Q+iJFR8mWiIjIHZhMJh5vWYMfXunMkJY1Mq1bFxJL13/tYuX+c6SlGQBcnTeP0MFDuBkYCEC7qu1Y13cdTzZ6EhMZ09ev+msVj218jN0XdhdeZ6RIKNkSERHJBa9Szrw/qClrx7fDt3Jp6/Lo+GTeWPcnAxb+ysE9R7j2+VISjh8n9PHBXJ75Pmlxcbib3Xmt1Wt81esr6papa60bHhfO+B3jmfLLFKISooqiW1IIlGyJiIjkgX8tL76Z0IGpjzailLOjdfnhsGgGbQ7jk0aPEufkCmlpRC5bRkifPsT8+CMAft5+rOmzhvF+43EyOVnrbj69mX4b+rE5ZDOGYRR6n8S+lGyJiIjkkZOjA6M71OaHVzrzaNMq1uVpmNhY9yGe6fEG22u0JA0TKRfDOT/+Oc5PfJHky5dxdnTmuWbPsbrPapqUb2KtG5UYxZTdUxi7fSxhN8KKoltiJ0q2REREbFS5jCvzn2jBV6NbZ5qbK8pcitn+Q3m10wsEl7k1OD7m++853as3kV9+hZGSgqWcheW9ljOp5STcnNysdX8L/43HNj3GZ39+pslQ7xFKtkRERPLpofrebP37Q7zS3YKrOeOrNahcLSZ2/jvz/AZyw+xOWlwcl99779YA+mOBODo4MrLxSDb020DH6h2t9RJTE5lzcA5DvxnK0YijRdElKUBKtkRERAqAi5MjE7rVZ8fLnQhoUtm63DCZ+K52O8Y88gbf+LQjFRMJx4+TdDbUWqaqR1Xmd53PrE6zKO+a8WLsk1EnGfHdCGbsm6FpIkowJVsiIiIFqHo5dxaO8Gf56DbU9c64tRhjduPfzQYysfPfOd3pUTx79cpUz2Qy0dOnJxv7b2SQZZB1uYHBihMr6LehHzvO7tAA+hJIyZaIiIgddKhfgS0vduTNXg3xcMl48vB02Wo8X64zL60+TPj1m9blKVevknDyJGVcyvB2u7f5oucX1ClTx7r+cvxlXtr1Es//8DxhMRpAX5Io2RIREbETZycHnulYh52vdGJA82qZ1m04fJGuH/7E3B2nuJmUyuUPPuDMgIFcnvk+qbGxtKjUgv/0+Q/PNXsOs0PG7PW/XPiFxzY+xqdHPiUpNamwuyQ2ULIlIiJiZxU9XZk9pBn/fbYdjap4WpffTE7lox0n6TLjezYdCMNISbk1N1dAANc3bcLsYGa833jW9V1HmyptrPUSUxOZf3g+AzcNZF/4vqLokuSBki0REZFC0tLHi80TOjBjwAOUL+VsXX7pZhrvtxrBqw89z8my1UmNuMrF117n7JNPknDiBD5lfFjcfTEfdPyACm4VrPVCb4Qy5vsxvP7z61y9ebUouiS5oGRLRESkEDk6mBjWuiY/TurM2I51MDtmvC/xePnavNj578xuPoRIl9Lc/OMAZwYM5NK0aaRGRxNQO4BN/TfxhO8TOJgyvsK/O/Mdfdb3YUXQClLTUouiW3IHSrZERESKgKermSm9GvL9S53o3qhSpnXba7VidPfJrLR0I8HkSNSKlZzuGUDUypV4OLrzRps3WNF7RaYZ6GOTY5mxfwZDvx3KoSuHCrs7cgdKtkRERIpQ7QqlWDyyJV+PaUODShkvuE5wcuHLRgGMefh1ttdoSdL1G1z65zQuvvEGAI3LN2Z5r+W81eYtSpsz6p2IPMHILSOZ8ssUIuIjCr0/kpWSLRERkWKgfb0KfDuxA+/0b0I594ynD6+5lWW2/1Amdv47h7zrU+7xx63rHB0cGeI7hE2PbaJv3b6Z2tt8ejN9NvThi8Av9NqfIqZkS0REpJhwcnTgyba12DWpC+M61sHZMeNr+nTZakxpP47njhn8dSkmU73yjp5M7zCdLwO+xNfL17o8LjmOD//4kIGbBvLrxV8LrR+SmZItERGRYqaMm5k3ejXkh1c60devaqZ1u/6KIGDuz0xee5QrNxJIOHmS4G4PE7VyJc28HmBV71VMbTuVMi5lrHXOXD/DuO3jeHnXy1yMvVjY3bnvKdkSEREppmp4ufPxsOZsfL49rWt7WZenGbDq9zA6f7iLmfM3cz3q1niu0/36c/OX3TxueZxv+n/DYMtgTGQ87bj97Hb6bejHgsMLiE+OL4ou3ZeUbImIiBRzfjXKsnpsWxY96U+dChnvW4xPSmWZuy9Pd5/C+roPEXvmLGHjnuXc00/jEnKBqe2msurRVfh5+1nrJKQmsPDIQvps6MPmkM2kGWlF0aX7ipItERGREsBkMvFI48pse6kj0/o1xuu2SVFvuJRi0QP9GPPw63xfsxUxv+0ndOAgLkx6jXo3y/BlwJdM7zCd8q7lrXWuxF9hyu4pDP92OIevHC6CHt0/lGyJiIiUIGZHB0a28+GnSZ2Z2K0+7s6O1nUR7uX4qMUQxnd5hT1VmnB982ZOBwQQ8f4senk9xDePfcPoJqMzvWvx2LVjPLnlSV776TWN57ITJVsiIiIlUGlXMy93t/Dza10Y9aBPppnowzwr8W6bUbzUcSKHy9S89b7Fh7vjcvUGf/f/O5v6b+KRWo9kam9L6Bb6bujLxwc/1niuAqZkS0REpASr4OHC/+vbmJ2vdGZAi2qYMnIu/vKqyeQO45ny4FiC/Ttjrnrrycbqpavzr87/YlnPZTT0amgtn5iayOI/F/Po+kdZd2qdXv1TQJRsiYiI3ANqeLkze3Aztr7YkYcbZn79z6GKFp736sLIz/dz6FyUdbmfkw8rH/mSd9u/i7ebt3V5xM0I3v71bQZuGsiusF0YhlFY3bgnKdkSERG5hzSoXJolT7Vk7fh2maaLAPj5ZASPLfiVUUv3cyQsmvC33+ZMQG86HU5m86PrGdt0LC6OLtbyIddDmLBzAqO2juJoxNHC7so9Q8mWiIjIPci/lherx7Zl+eg2+Ncql2ndrr8i6PfvPbwSW4ugOBPhb00lvN/jPBlShU2919K3bt9M83MdvHKQ4d8N5+VdLxN6PbSQe1LyKdkSERG5R5lMJjrUr8B/n23HF0+3plmNspnW76/ciAldXmJam1EE3Ugj/K2pxA38G6+GNWF1zxW0r9Y+U/ntZ7fTf2N/3v3tXa7evFqIPSnZlGyJiIjc40wmE50s3qx/7kGWjmpF0+plMq3fW6UJE7q8xNR2Yzic6Mqlf07DafALTA9ry6KH5tOofCNr2VQjldV/rabXul7MPzSfG0k3Crs7JY6SLRERkfuEyWSii29FNj7fns+eaknjqp6Z1v9RyZdJDz3Pqx2eY69Rjsvvv4/X0NeYG9KO9zvMpJpHNWvZmyk3+fTop/Rc25NFRxcRlxxX2N0pMZRsiYiI3GdMJhPdGlbimwkdWPSkPw9Uy3ylK7BCHf7x4Bhe6PwSuzxqczPwOL3q9mZz/81Mbj2Zci4ZY8BikmKYd2gePdf2ZOmxpdxMuVnY3Sn2lGyJiIjcp9JfAbTphfZ8Nbo1betkfnrxdNlqzGj9JKNqDmD17+cwDEeGNxzOdwO+49kGT+Nh9rCWjU6MZvaB2QSsDeDroK9JTE0s7O4UW0q2RERE7nMmk4mH6nuzamw71o5vRzffipnWn41N5fW1f9Lh/Z3M33mKpGRn+iwOZOF/yjLC1A43Rzdr2WsJ15i5fya91/VmzV9rSE5NLuzuFDtKtkRERMTKv5YXn41qxZYXH6KvX1UcbpuR/kpMIh9+f5J27+1gRmxlLofF0/e9X5i3MIVB0fVxMWW8c/Fy/GXe+e0deq/vzcoTK0lISSiC3hQPSrZEREQki4ZVPPl4WHN2vtKZYa1r4OyYkTIkphpsqd2OcQ+/xtR2YzjlUIPHFwbx8dyb9D5dFjMZL8cOjwvnvX3v0WNtDz778zNik2KLojtFSsmWiIiI5MinQilmDGjKnsldebFbfcqXcs60/o9KvrzVfizju77KvgqtGfqfaOb+O5GHj4BjWsZlsciESOYcnMMj/32EeYfmEZkQWdhdKTJKtkREROSuvEu78FJ3C3smd+WDgU1pUKl0pvVnPSszt/lgRvT8B2tr9qXHL+WZtzCZgAPgetsrgGKSY1h0dBE9/tuD9/e/z6W4S4XdlUKnZEtERERyzdXsyOBWNdj694dYProNXRp4Z1of6+zOhnodebbbJGY2fQ6LMZiN/bbwzAPPZHp6MSE1geVBywlYF8DUPVM5EXmisLtSaJyKegdERESk5El/FVCH+hUIvhLL0j1nWHfwAjeTU61lAivUIRCYN/sQA1p0Y2HHx/kj6lu+OLiYaOIBSElLYUPwBjYEb6B56UaMaP40XWt1w8nh3klR7p2eiIiISJGoV9GD6Y89wOQAXzYcvsiKfecICs94jc/1m8ks3RPK0j2htPJpxOBT/YiNWM22lqlcK5MxrutQzHEO/fwqFW6a6Rtbj74eD+Jd00LpgABMpoxyhmFk+rm4U7IlIiIiBaK0q5kn29ZiRJuaHD1/nZX7z7HpyEXikzKudv0eGsXvZj/MlZvQ4uAxmpt/JqTuBc5UyWjnqlsyn7sF8VXycTr+4sL4tnVp4NXAuj4pOJgLr71OjU8WYq5UqTC7aJMSn2z99ttvLF26lCNHjhAfH0/VqlXp2bMnY8eOxd3d3aY2t23bxvLlyzlx4gTJycnUqlWLvn37MnLkSMxm890bEBERuY+ZTCb8apTFr0ZZ3uzdkI3/d7Xr+G1Xu5IdHNlX2Q/wwy06habX9pPm/iOnat0g1fHWVatks4kffJP4YfMgWlRsQd+6fenu0x3zzZskBgUR/sYb1FiyBJND8R6CbjIMwyjqnbDVV199xfTp0zEMg8qVK+Pl5UVwcDBJSUnUrVuXFStWULZs2Ty1+f777/P5558DULNmTdzc3AgODiY1NZVWrVrx+eef4+zsfJdW8mbAgAEEBgbSuHFj1q1bV6Bti4iIFAeGYXD0/HX+cyCM7/68RGRcUrblPN1iqO29h8uu+4hzzPqeRWcHZx50tNBq5RGahxhUe/MtvIYPt/fuZyu339/FOxW8g2PHjvHee+8BMG3aNHbt2sX69evZsWMHjRs3JiQkhKlTp+apze3bt1uTqQULFrB9+3Y2bdrE5s2bqV69Or///juzZ8+2R3dERETuaelXu97t/wD7pnRj2d9aMaBFNTxcMt9ku3GzNEfO9eTSqTdJvjQEN8Mn0/qktCR2JR9j1iBHxk50ZPqB9/n9yBaK87WjEptsLViwgLS0NPr168eQIUOsA+UqVarE7NmzcXBw4Pvvv+fEidw/Sjp//nwAnnnmGbp162ZdXrduXd59910Avv76ayIj75+J2ERERAqa2dGBzg0qMntwM/5462EWDm9BQJPKODvdlpYYTiRENefKiXHEnXmepMgHIdUjUzuxbia+9zN4+vBr9FrXi7kH5/Lz+Z+5nni9kHt0ZyVyzFZcXBy//PILAIMHD86y3sfHh7Zt2/Lrr7+ydetWfH1979pmaGioNTEbMmRIlvXt2rWjVq1anD17lh9++IHHH388n70QERERV7MjAQ9UIeCBKsQkJPN94GW2HLvE3pCrxCWlAibSEmqQmFCDxMu9cfQ4hdnzIE6lj2NySLG2cz72PEv+XGL92cfTBz9vP/wq+tG0QlPqla2Ho4NjNntgfyUy2QoKCiIpKQlnZ2eaNm2abRl/f39+/fVXjhw5kqs2Dx8+DECNGjWolMOTDf7+/pw9e5YjR44o2RIRESlgpV3NDPSvzkD/6iSlpHHwXBS/nIrg55NX+fPCdcCR1FhfUmN9wSEBp9LHMJc5hKP7aUymzLcRQ2+EEnojlI0hGwEoZS5FkwpNeLDqgwzzHYabk1uh9atEJltnzpwBoGrVqjk+HVizZs1MZe8mNDQ0U738trlq1SrWrFmTq22HhITkqpyIiMj9wtnJgbZ1ytO2Tnkm9YBrsYnsDr7Kzyev8supCK7EQMr1lqRcb4nJKRonz6M4uZ/Bwe0cDk5xWdqLS45jX/g+9oXvIzohmpdbvlxofSmRydb167fuxZYpUybHMunr0ssWZJs3btzIsUy6iIgIAgMDc7VtERERubPyHi70a1aNfs2qYRgGIRFxBF+JJeiPQE6FGFx0e5SQiHhiEpMxma/h6Hbu1j/3c5jdLpNmZMz1ZVC4g+lLZLKVmJgIcMc5r9KnZ0gvW5BtJiQk3LU9b29vGjdunKtth4SE5KpNERERufVkY72KHtSr6EHPJpWtyw3DICI2kZArcYRExBJ+/Sbt6lSghY87gdcCORpxFEeTI0N9hxbq/pbIZMvF5dbbw5OTk3Msk5SUlKlsQbbp6up61/aGDh3K0KG5+zDT5+kQERER25lMJiqWdqViaVfa1S2faV2ryq1oVblVkexXiZz6ITe3CHNzW/B2np6euW4zvayIiIjI3ZTIZMvHxweAixcv5ngl6ty5c5nK3k3t2rUBOHv2bI5l8tqmiIiISIlMtho2bIjZbCYpKYmjR49mW+bAgQMANGvWLFdt+vn5AXD+/HkuX75cIG2KiIiIlMhky8PDgw4dOgBkO71CaGgov/32GwA9e/bMVZu1a9fGYrEAsHr16izr9+7dy9mzZzGbzZlmlxcRERG5kxKZbAE899xzmEwmNm7cyOrVq63vRLpy5Qovv/wyaWlpPPzww1lmj+/atStdu3Zl69atWdp84YUXAFi8eDE7d+60Lj99+jRvvfUWAE888QReXl726paIiIjcY0rk04gATZs2ZfLkycycOZN//OMfLFy4kHLlyhEcHExSUhK1a9fmnXfeyVLvwoULAMTHx2dZ16NHD5566im++OILxo8fT82aNXF3d+fUqVOkpqbi7+/PK6+8Yve+iYiIyL2jxCZbAKNGjaJBgwZ8/vnnHD16lGvXrlG1alV69uzJ2LFjKVWqVJ7bnDJlCs2bN2fFihUEBQVx5coV6tatS9++fRk1atQd5+ESERER+V8lOtmCWy+IbteuXa7L//XXX3ctExAQQEBAQH52S0RERAQowWO2REREREoCJVsiIiIidqRkS0RERMSOlGyJiIiI2FGJHyB/Lzh//jwAISEhDBgwoIj3RkRERHIjJCQEyPgez4mSrWIgMTERgISEBAIDA4t4b0RERCQv0r/Hc6Jkqxjw8vIiMjISFxcXqlevblMbISEhJCQk4OrqSt26dQt4D+9Piql9KK4FTzG1D8W14N1rMT1//jyJiYl3fbOMyUh/z42UaAMGDCAwMJDGjRuzbt26ot6de4Jiah+Ka8FTTO1DcS1492tMNUBeRERExI6UbImIiIjYkZItERERETtSsiUiIiJiR0q2REREROxIyZaIiIiIHSnZEhEREbEjJVsiIiIidqRkS0RERMSOlGyJiIiI2JHejXiPGDx4MBEREXh7exf1rtwzFFP7UFwLnmJqH4prwbtfY6p3I4qIiIjYkW4jioiIiNiRki0RERERO1KyJSIiImJHSrZERERE7EhPIxYiwzA4dOgQO3fu5MCBA5w+fZrY2FhKly5No0aN6N+/P3369MFkMmVbPy4ujkWLFrFt2zYuXryIu7s7fn5+PP3007Rp08bm/bp27RoLFy7kxx9/5MqVK3h6etKqVSvGjRtHw4YNbW63MNga06ioKHbs2MGvv/5KYGAg4eHhODg4UKVKFTp06MCoUaOoXr16nvfn/PnzdOvW7Y5l/Pz8WLNmTZ7bLkz5+V1t0KDBHduuUKECe/bssWm/7HUMFAZbYzp58mTWr1+fq23s3LmTatWq5Xqf7PlZFaYtW7ZYj+UrV64QHR2N2WzGx8eHTp068dRTT1GuXLls6+q8mj1bYqrzas70NGIh2rt3L6NGjbL+XKNGDTw9Pblw4QLR0dEAdO7cmXnz5uHs7JypbmRkJE888QRnzpzB2dmZevXqERkZyaVLlzCZTEydOpXhw4fneZ/Onj3LE088wdWrV3F3d6d27dpcunSJa9euYTabmTt37l1/yYuSrTF94oknOHDgAAClSpWiZs2aJCQkEBYWRkpKCu7u7sydO5eOHTvmaX9uPym0aNEi2zL169dn2rRpeWq3sOXndzX9C7xJkyZZ1gGULVuWhQsX5nmf7HUMFBZbY/rJJ5/w008/5dhuaGgokZGRVKlShZ07d+LgkPsbFvb6rApbv379OHHiBM7Oznh7e1OuXDkiIyO5ePEiAOXLl+fzzz/H19c3Uz2dV3NmS0x1Xr0DQwrNnj17jK5duxpffPGFcfXq1Uzr1q9fbzRp0sSwWCzGBx98kKXus88+a1gsFuOxxx4zLl26ZBiGYaSlpRmrVq0yLBaL0bBhQ+P48eN52p+0tDSjX79+hsViMUaPHm3cuHHDMAzDSE5ONubMmWNYLBajWbNmxuXLl23ssf3ZGtMRI0YYL7/8srF//34jJSXFujw8PNz429/+ZlgsFqNFixZGREREnvYnLCzMsFgshsVisb1TxUB+flfT+x8WFlag+2SPY6Aw5SemOUlLSzO6detmWCwWY/bs2XneJ3t9VoVt9erVxv79+42kpKRMy0+cOGE8+uijhsViMXr16pWlns6rObMlpjqv5kzJViGKiYnJ8ot7u4ULFxoWi8Vo3bq1kZqaal0eGBhoWCwWw9fX1wgNDc1Sb9KkSYbFYjFeeOGFPO3P9u3bDYvFYvj7+xvR0dFZ1g8fPtywWCzGjBkz8tRuYbI1ppGRkTnWuXHjhtG2bVvDYrEYS5cuzdP+3AsnBcOwPa6GYZ8vcHsdA4UpPzHNyb59+6zxzi4ud3OvJFt3cuTIEWs/g4ODrct1XrVdTjHVeTVnGiBfiDw8PDCbzTmuT7+0Gh0dTWRkpHX5tm3bAGjbti21atXKUm/IkCEA/PTTT8THx+d6f7Zs2QJAz549KVOmTJb1gwcPzlSuOLI1pjmN3wAoXbo0zZo1A+DMmTMFs6MljK1xtRd7HQOFyR4xTR/L1aJFi2zjIlCnTh3r/2/evGn9v86rtssppjqv5kwD5IuRhIQE6/9dXV2t/z98+DAALVu2zLZe06ZNcXZ2JjExkaCgIPz9/XO1vSNHjtyx3fTlly5d4vLly1SqVClX7RYnOcX0bhITEwFwc3Ozedvvvvsup0+fxmQyUa1aNTp06MDDDz+cpzE1xVVu4rpgwQKuXLlCamoqlSpVom3btvTq1SvbsUF3Y69joDjJ6+9qfHw8W7duBWDAgAH52nZBflbFTfoYovSxU+l0XrVdTjG9m/v5vKpkqxj59ttvAfD19cXDw8O6PDQ0FICaNWtmW89sNlOlShXOnj3LmTNncnVSSEpK4sKFC3dst0qVKpjNZpKTkzl9+nSJPCnkFNM7uXz5Mvv37wdyPmHmxldffZXp59WrV9OwYUPmzZtHjRo1bG63OMhNXNeuXZvp5/Xr1/Pxxx8zb948GjdunKft2eMYKG7y+ru6bds24uPjcXNzIyAgIF/bLsjPqjhIS0sjIiKCPXv28OGHHwLw6quvUqpUKWsZnVfzJjcxvZP7/byqZKuYOHbsGKtWrQJg7NixmdZdv34dINtL0unS1924cSNX24uNjSUtLe2O7ZpMJjw9Pbl27Vqu2y1O7hTTO3nnnXdITk6mXr16dOnSJU/bdHJyom/fvvTu3Zt69epRsWJFoqKi+Omnn5gzZw5BQUGMHj2adevW5Tr5K27uFtdu3brRr18/fH19qVy5MnFxcezdu5ePPvqIsLAwnn76aTZs2ECVKlVyvU17HAPFiS2/q+vWrQOge/fuNv8u2eOzKkrLli1jxowZmZY1bdqUmTNnZnkCTufV3MlLTO/kfj+vFu/rbveJq1evMmHCBFJSUujevTu9e/fOtD790uudxnukX+6//VbEnaS3eXvdgmi3uLhbTHOyaNEitm/fjtlsZubMmTg6OuZpu5UrV2bWrFl07tyZ6tWr4+zsTKVKlRg8eDArV67Ew8ODs2fP8uWXX9rSrSKXm7guWLCAHj16UKtWLVxcXPDy8qJ3796sWbOGqlWrEh0dzfz58/O0XXscA8WFLb+rYWFh/P7770D+biHa47MqSpUqVaJFixb4+fnh7e2NyWQiKCiIjRs3ZklsdF7NnbzENCc6ryrZKnIxMTE888wzXLx4kcaNGzNz5swsZVxcXABITk7OsZ2kpCQg9+OS0tu8vW5BtFsc5Cam2Vm/fj2zZ8/GZDIxffp0HnjggQLdr1q1ajFs2DAAtm/fXqBtFwZb45rOy8vLetVmx44dGHmY4s8ex0BxYGtMN2zYgGEYVKtWjbZt2xb4fuXnsypKAQEBrFy5kjVr1rB79242bNiAn58f33zzDSNHjiQ1NdVaVufV3MlLTLOj8+otSraKUFxcHGPGjOH48ePUr1+fzz77LNtLoJ6enkDGZe/spK9LL3s3Hh4e1gGFObVrGIb1L5fctlvUchvT//Xdd9/x5ptvYhgGb7/9Nv369bPL/jVv3hzIGC9SUtga1/+V3v/o6GjrRJ65YY9joKjZGlPDMNiwYQNwa+LJnN44kV+2flbFia+vL59++inlypUjKCjIOi4OdF611Z1i+r90Xs2gZKuI3Lx5k3HjxnH48GF8fHxYunRpjo/N+vj4ALdmJc5OcnKydVbf9LJ34+zsTNWqVQE4d+5ctmXCw8Otf/Xl5YmTopKXmN5u+/btTJo0idTUVF5//XXrX0n2kH7L4m5/DRYntsY1O7ffsslLDOxxDBSl/MR0//79nD9/HpPJlO+nEO/E1s+quPHw8KB169YABAYGWpfrvGq7nGJ6O51XM1OyVQQSExMZP348v//+O9WqVWPZsmV4e3vnWD59bpL0x23/19GjR0lOTsbFxSVP79xKb/ePP/7Idn368sqVK1O5cuVct1sU8hrTdD/99BMvvfQSKSkpTJw4kaefftqu+3nq1CmAYh/PdLbGNSfp/XdxcaFs2bK5rmevY6Ao5Dem6XNrtWzZ0q5PX9n6WRVHKSkpQOYvY51X8ye7mKbTeTUrJVuFLDk5mQkTJrB3714qVarEF198cdcnfXr06AHAvn37sv0rbPXq1cCtSRFz+xju7e1u3bo120ve6S/17NmzZ67bLAq2xBRuvatuwoQJJCcnM27cOJ5//nm77mdcXBwrVqwAoH379nbdVkGwNa45SUlJYenSpcCtiSSdnHL/MLS9joHClt+YxsXFWSfjfOyxx+y1m/n6rIqb6Oho65QDtydNOq/aLqeYgs6rOSq6yevvPykpKcaECRMMi8VitG/f3jh9+nSu644bN876Dq/0d2rd/g4vX19f49ixY1nqLV261OjSpYsxdOjQLOtSU1ONPn36GBaLxRgzZoz1HV4pKSnWd3j5+flZ3xlWHNka04MHDxrNmjUzLBaL8d577+Vpm3eK6VtvvWVs27bNSExMzLQ8ODjYGDp0qPW9aMX99Si2xnXWrFnGunXrjJiYmEzLL168aIwfP96wWCxGo0aNjMOHD2epu2XLFqNLly5Gly5dsm3b1mOguMjP8Z9u7dq11t+h2NjYu5a/U0zz81kVJ/v27TP+/e9/Z3tMHTt2zBgwYIBhsViMhx56KEvMdF7Nnq0x1Xk1ZyX3z5USaMuWLda/Sp2dnZkyZUqOZadOnUqjRo2sP7/33nsMGzaMwMBAunXrRr169YiKiiI8PByTycSUKVOynXwwJibGOsne/3JwcGDu3LkMHz6cn3/+mY4dO2Z5O/2sWbOK9aR7tsZ08uTJxMfH4+joyNGjR3McT9CoUSOmTp2aadmdYnr06FHWrFmD2WymZs2aeHh4EBUVZR2/UaZMGebMmUP16tXz3NfCZGtcT58+zeLFi3nzzTepUaMGZcqUISYmhjNnzmAYBi4uLrz77rv4+fllaSc+Pj7HuILtx0BxkZ/jP136LcQePXrk6mrLnWKan8+qOLlx4wZz585l7ty5eHt7U7FiRRwdHQkPDyciIgK4NX3Bp59+miVmOq9mz9aY6ryaMyVbhej2R4EvXLhwxy+WmJiYTD97eXmxdu1aFi9ezNatWwkODsbd3Z2OHTsyevRomx//rl27Nps2bWLhwoX8+OOPnDx5Ek9PT3r06MGzzz6b7Qm/OLE1pukDVFNTUzl48GCOdfJ6+2TcuHH88ssvHDt2jKtXr3L27FlcXV1p3LgxHTt2ZPjw4fka81RYbI3rsGHDqFChAseOHePKlStcuHABs9lM/fr1adeuHSNGjMhxZu27sdcxUFjyc/xD5rm1CuIWoj0/q8LUvHlz3njjDfbt20dwcDChoaEkJSXh6elJmzZt6Nq1K4MGDcr2SU+dV7Nna0x1Xs2ZyTBKyAQqIiIiIiWQBsiLiIiI2JGSLRERERE7UrIlIiIiYkdKtkRERETsSMmWiIiIiB0p2RIRERGxIyVbIiIiInakZEtERETEjpRsiYiIiNiRki0RERERO1KyJSIiImJHSrZERERE7EjJloiIHX300Uc0aNCARYsWFXjbo0ePpkGDBuzdu7fA2xaRgmMyDMMo6p0QESlMb7zxBuvWrbOpbtOmTfnPf/6Tq7KXLl2iR48euLu788MPP+Du7m5dN2/ePObPnw/AX3/9lWMba9euZerUqaSmptKyZUs++eQTSpcuDcDhw4cZMmQIjRo1Yu3atTg46O9nkeLIqah3QESksJ0/f54KFSpkWR4fH098fDxAtusBWrdunevtfPTRRyQkJDBx4sRMiVZuLVu2jJkzZ2IYBl26dGHOnDm4urpa1zdr1owOHTqwe/duNm3aRP/+/fO8DRGxPyVbInLf+eqrr7JdPnXqVNasWUOVKlXYtWtXvrZx+fJlNm/ejNlsZuDAgXmuP2fOHBYuXAhAnz59mDlzJk5OWU/Zw4YNY/fu3SxZskTJlkgxpWvOIiL/JzAwEIBGjRrlu601a9aQmppKp06dKFu2bK7rGYbBtGnTrInWiBEjmDVrVraJFkDHjh0pW7Ysp06d4sCBA/nebxEpeEq2RESAlJQUTp06BeQ/2TIMg//+978APProo3nah0mTJvH1118D8PzzzzN16lRMJlOOdZydnXnkkUeAWwmeiBQ/SrZERIDg4GCSkpKA/CdbJ0+e5NKlSwC0bNkyV3USExN54YUX2Lx5MyaTiTfffJOJEyfmqm76Nnbv3m3bDouIXSnZEhEh4xYiQOPGjfPV1h9//AFAlSpV8Pb2vmv52NhYxowZw48//oiTkxMzZ85k5MiRud6en58fAFevXiUkJMS2nRYRu1GyJSICBAUFAVC+fHkqVaqUr7aOHDkCgK+vb67Kjxw5kv379+Pi4sK8efPyPNDdx8fH+rTj4cOH81RXROxPyZaICBlXtho2bJjvtq5cuQJAuXLl8rTtQYMG0bVrV5u2mb6t9G2LSPGhZEtE7ntpaWmcOHECyP8tRIDIyEiAXD+F2Lx5cwC+/vprvvzyS5u2WaZMmUzbFpHiQ8mWiNz3zpw5Y53MtCCmfUhMTARuPSmYG0uWLKFFixYATJ8+nWXLluV5my4uLpm2LSLFh5ItEbnvpY/XgoK5spV+S+/69eu5Ku/h4cGSJUvw9/cHYMaMGSxdujRP20zfVm5vXYpI4VGyJSL3vfQxU6VLl6Z69er5bi+vyRZAqVKlWLJkCa1atQJg5syZfPbZZ7mur2RLpPhSsiUi973jx48DtwbH32kC0dyqV68eAGFhYXmq5+7uzqJFi6zvX/zggw9YvHjxXevFxsYSFRUFQN26dfO4tyJib0q2ROS+V5CD4yFjktETJ05YJ0rNrfSEq02bNgB8+OGHLFq06I51jh07RlpaGk5OTtaxXyJSfCjZEpH72vnz54mOjgYKZtoHAH9/f5ycnEhOTs40Hiy33NzcWLRoEe3atQPgX//6F5988kmO5Y8ePQrcGtxfqlQp23ZaROxGyZaI3NfSbyFCwV3Z8vDwoFOnTgDs3LnTpjZcXV355JNPaN++PQAfffQRCxYsyLZs+jby8h5GESk8SrZE5L6Wnmy5ublRp06dAmt3yJAhAGzevBnDMGxqw9XVlYULF9KhQwcA5s6dy/z58zOVCQsL49ChQ7i6uuZ55nkRKRwmw9azgIiI5CgtLY0ePXpw7tw5li9fbn3KsKDNnz+fefPmMWDAAGbMmGGXbYhI/ujKloiIHTg4OPDiiy8C3HWAu63i4+NZvnw5zs7OvPDCC3bZhojkn5ItERE76d27N02bNuXnn3+2DmIvSMuXLycqKoonn3ySatWqFXj7IlIwnIp6B0RE7lUmk4lp06axY8cOu7yz0N3dnQkTJvDUU08VeNsiUnA0ZktERETEjnQbUURERMSOlGyJiIiI2JGSLRERERE7UrIlIiIiYkdKtkRERETsSMmWiIiIiB0p2RIRERGxIyVbIiIiInakZEtERETEjv4/V1wSlszwzkAAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1, 1)\n", "ax.plot(tr * Tc, s_corr, color='tab:red', label=\"Mulero et al. (2012)\", ls=\"--\")\n", "ax.plot(dia_fh1.liquid.temperature/KELVIN, surf_tens_fh1, color='tab:blue', label='FH1-DFT')\n", "ax.plot(dia_fh2.liquid.temperature/KELVIN, surf_tens_fh2, color='tab:green', label='FH2-DFT')\n", "ax.set_xlim([19.5, 34.0])\n", "ax.set_ylim([-0.1, 2.2])\n", "ax.set_xlabel(r'$T$ (K)')\n", "ax.set_ylabel(r'$\\gamma$ (mN/m)')\n", "legend = ax.legend(loc='best', frameon=False, fontsize=12, ncol=1)" ] } ], "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.10.6" } }, "nbformat": 4, "nbformat_minor": 5 }