{ "cells": [ { "cell_type": "markdown", "id": "fbc5e807", "metadata": {}, "source": [ "## SAFT Equations of State" ] }, { "cell_type": "markdown", "id": "85bef0a4", "metadata": {}, "source": [ "### PC-SAFT" ] }, { "cell_type": "code", "execution_count": 1, "id": "bf5e0da8", "metadata": {}, "outputs": [], "source": [ "import yaeos\n", "\n", "m = [1.0582, 3.3004]\n", "epsilon_k = [145.5257, 224.0780]\n", "sigma = [3.6316, 3.8639]\n", "model = yaeos.PCSAFT(m, sigma, epsilon_k)" ] }, { "cell_type": "code", "execution_count": 5, "id": "16e01088", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAGyCAYAAAACgQXWAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvc2/+5QAAAAlwSFlzAAAPYQAAD2EBqD+naQAAZxpJREFUeJzt3Xd4VFX+BvB3JpNM+qT3Sk0ghRJKQAUlCMJSBHQpKig/XV1wAdcCFrCgCLoKNlBRQBfEVWmCoPQuPXRCC6SQQkiZtMm08/tjwshQk5Dkzkzez/PMk+TeO8P3cmHmzTnnniMTQggQERER2Sm51AUQERERNSSGHSIiIrJrDDtERERk1xh2iIiIyK4x7BAREZFdY9ghIiIiu8awQ0RERHaNYYeIiIjsGsMOERER2TWF1AVYA6PRiEuXLsHDwwMymUzqcoiIiKgGhBAoLS1FSEgI5PLbtN8ICW3dulX87W9/E8HBwQKAWL58uXmfVqsVL7/8soiLixOurq4iODhYPP744yI7O9viNa5cuSJGjhwpPDw8hEqlEk899ZQoLS2tVR2ZmZkCAB988MEHH3zwYYOPzMzM237OS9qyU15ejsTERDz11FMYMmSIxb6KigocPHgQb7zxBhITE1FUVIQJEyZg4MCB2L9/v/m4UaNGIScnB+vXr4dOp8OTTz6JZ555BkuWLKlxHR4eHgCAzMxMeHp61s/JERERUYNSq9UIDw83f47fikwI61gIVCaTYfny5Rg8ePAtj9m3bx86d+6MixcvIiIiAidPnkSbNm2wb98+JCUlAQDWrVuHfv36ISsrCyEhITX6s9VqNVQqFUpKShh2iIiIbERNP79taoBySUkJZDIZvLy8AAC7d++Gl5eXOegAQEpKCuRyOfbs2XPL16mqqoJarbZ4EBERkX2ymbCj0WjwyiuvYMSIEeb0lpubi4CAAIvjFAoFfHx8kJube8vXmjFjBlQqlfkRHh7eoLUTERGRdGwi7Oh0Ojz66KMQQmDu3Ll3/XpTpkxBSUmJ+ZGZmVkPVRIREZE1svpbz68GnYsXL2LTpk0WfXJBQUHIz8+3OF6v16OwsBBBQUG3fE2lUgmlUtlgNRMREZH1sOqWnatB58yZM9iwYQN8fX0t9icnJ6O4uBgHDhwwb9u0aROMRiO6dOnS2OUSERGRFZK0ZaesrAxnz541/5yeno7U1FT4+PggODgYw4YNw8GDB7F69WoYDAbzOBwfHx84OTkhNjYWffv2xdNPP4158+ZBp9Nh/PjxGD58eI3vxCIiIiL7Jumt51u2bMH9999/w/bRo0fjzTffRHR09E2ft3nzZvTs2RMAUFhYiPHjx+PXX3+FXC7H0KFD8cknn8Dd3b3GdfDWcyIiIttT089vq5lnR0oMO0RERLbHLufZISIiIqothh0iIiKyaww7REREZNcYdoiIiMiuWf2kgkRERAAghECV3ogKrQFavRFymWkRabkMUDjI4a5UwEEuk7pMskIMO0REJDmDUSCrqALnLpfh/OVy5JRokKvWIF9t+lpSoUOF1gC98fY3EHs4K6BycYTKxRFBns4I93FFmLcLwrxdEe3nhmb+bnB0YKdGU8OwQ0REjapKb8CpnFIcySpGamYJjl8qwfmCcmj1xhq/hkIug4ApJF2rVKNHqUaPrKJKHL+kvuF5Tgo5Wgd6oG2IJ9qGqtA5ygetAt0hk7FFyJ4x7BARUYPS6o1IzSzGzrMF2H3uClIzi6E13BhsnBRyNPNzQ3N/d4R6uyDQ0xlBns4I9FTCy9UJ7koFXJUOcHOy7K4SQkBnEFBrdCiprH5U6HCppBKZhZXIKqpAZlElzuWXoaxKj6PZJTiaXQLsMy0C7efuhK7NfNGtuR/ubemHcB/XRvu7ocbBSQXBSQWJiOrb5dIqbDyZh/Un8rDr3BVU6gwW+71cHZEQ5oV2YSrEh3mhdaAHQr1dGnTMjdEokFFYgRM5ahy/VILDmSXYf7EQGp1l8IoPVaFvXBAeigtCM/+az8ZPjY8zKNcCww4R0d3LLdHg18OXsPZYDg5lFuPaTxcfNyckN/dF9+Z+6NrMB9F+blbRdVSlN+BwZgl2n7uCnWcLsP9iIa7tGYsJ8sCwjmEY2iEM3m5O0hVKN8WwUwsMO0REdaPW6LDuaC5WpGZj9/krFgEnIUyF3rGBeCA2ALFBnpDbwJ1SV8qq8MeJPKw9lotdZwvMA6KdHOToExeEEZ3DkdzM1yqCGjHs1ArDDhFRzQkhcDirBEv2XMSvh3MsuqiSIr0xsF0IHmwThCCVs4RV3r2SCh1+PXIJS/dl4Fj2X4OdWwa449kezTGwXQjv7JIYw04tMOwQEd2ZRmfAytRsLNp1ESdy/vrwb+7vhiEdwjAwMcRuB/cezSrBD/sysPJQNsq1pnAX6uWCZ+5rhkeTwuHi5CBxhU0Tw04tMOwQEd3albIqfP/nRXy/+yKulGsBmO6c6h8fjJFdIpAU6d1kunXUGh3+++dFfLsjHQVlpr8LXzcn/KtXS4zsEsGWnkbGsFMLDDtERDfKKanE3C3nsHRfpnkOnBCVM8Z0j8IjHcOb9IBdjc6Anw5k4att55BZWAkAiPZzwyt9W6NP26AmE/6kxrBTCww7RER/MYecvZnm+XASwlT4v3uboV9cEBRsvTDTGYxYui8TczacNrf0dIjwwtQBbdEu3Eva4poAhp1aYNghIjLdOj53y1n8cE3I6Rztg4m9WiK5Oe9Aup2yKj2+2noOX29PR6XOAJkMGJ0chX8/2Aoezo5Sl2e3GHZqgWGHiJoytUaHLzafw7c7083dVZ2jfDCxd0veZl1LeWoN3l97CssPZQMAgjyd8dagtujTNkjiyuwTw04tMOwQUVOkr+6C+Xj9afPA405R3piU0ootOXdp+5nLeG35MWQUVgAA+rQNxHsPx8PXXSlxZfaFYacWGHaIqKnZevoy3l1zAqfzygAAzfzd8Fq/WDwQE8CQU08qtQZ8sukMvt52HnqjgL+HEh89moh7W/pLXZrdYNipBYYdImoqMgsrMHXlMWxOuwzAtEbVxF4tMaprJG+bbiDHL5Vg4tJUnMk3Bcv/uycaL/VtDaWCc/PcLYadWmDYISJ7pzMYMX97OuZsPA2NzghHBxlGJ0fh+QdaQuXKAbQNrVJrwLu/ncB//8wAALQJ9sRnI9tzodG7xLBTCww7RGTPDlwswqvLjiItrxQAkNzMF9MfjkNzftA2uvUn8vDyz4dRVKGDh1KB2cPboVdsoNRl2SyGnVpg2CEie6TW6DBz7Sks3mNqTfB2dcTr/dtgSIdQjsuRUJ5ag/FLDmLfhSLIZMALKa0w7v4WNrFQqrVh2KkFhh0isjc7zxbgpZ8O41KJBgDwSMcwTOkXC58mPOuxNdHqjXh79XFzt1bftkH48NFEuCsVEldmWxh2aoFhh4jsRaXWgPfXnsSi3RcBABE+rpg5NAHJzX0lroxuZuneDExdeRxagxFxoZ74dkwnBHjY9mrxjammn98cek9EZCcOZhSh/yfbzUHnsa4RWDvhXgYdKza8cwR+eKYrfN2ccCxbjSFf7MK5y2VSl2V3GHaIiGyc3mDER3+kYdjcXThfUI4gT2cseqozpg+Ohxu7Raxex0hvLPtnN0T5uiKrqBJD5+7CgYtFUpdlVxh2iIhsWE5JJUZ8/Sc+2XQWRgE83D4Uv0+8Dz1aceI6WxLp64ZfnuuGxHAvFFfoMPLrP7HpVJ7UZdkNhh0iIhu14UQeHpqzHfsuFMFdqcCnI9rj47+347w5NsrXXYkfnu6CXjEBqNIb8Y/vD+D347lSl2UXGHaIiGyMVm/EO6tP4P++24/iCh3iQ1VY8697MCAxROrS6C65Oikw7/GO6J8QDJ1BYNzig/jtaI7UZdk8duYSEdmQ7OJK/PO/B3A4qwQA8FT3aLzyEJcesCeODnLM+Xs7KOQyrEy9hOd/OAS9UWAgw2ydMewQEdmIXecKMH7JIRSWa6FyccSHjySidxvOvmuPFA5yfPRoOyjkcvxyMAsTlx6CDGDrXR0x7BARWTkhBL7ZkY4Za0/BYBRoG+KJeY91RLiPq9SlUQNykMvwwbAEKOQy/Lg/E5N+TIWHswI9WwdIXZrN4ZgdIiIrVqk1YMLSVExfcxIGo8CQDqH45bluDDpNhFwuw3tD4jEgMQR6o8Cz/z2A/RcKpS7L5jDsEBFZqeziSgyZuwurDl+CQi7DmwPa4D+PJMLZkeNzmhIHuQz/eSQRPVv7Q6Mz4smF+3DiklrqsmwKww4RkRU6lFGEQZ/txMkcNfzcnbDk6a4Y0z2aC3g2UU4KOeaO6oikSG+UavR44tu9yLhSIXVZNoNhh4jIyqw+cgnDv/oTBWVViA32xKrx96BztI/UZZHEXJwc8M2YTogN9kRBWRWeWrQPJZU6qcuyCQw7RERWQgiBTzeewfglh1ClN6JXTAB+ejYZIV4uUpdGVkLl4oiFT3ZCkKczzuaXYfySg9AZjFKXZfUYdoiIrIBWb8S//3cY/1l/GgAw9p5ofPVEEty5thVdJ9DTGfNHJ8HF0QHbzxTgzVXHIYSQuiyrxrBDRCSxUo0OTy3ch2WHsuEgl+Hdh+Pwxt/awEHO8Tl0c3GhKswZ3g4yGbB4TwYW7LwgdUlWjWGHiEhC+aUa/P3LP7HjbAFcnRzw7ZhOGNUlUuqyyAY82DYIUx6KAQBMX3MC289clrgi68WwQ0QkkfOXyzB07i6cqL7jaukzXblaOdXK0/c2w6NJYTAK4F8/HEJ2caXUJVklhh0iIgmkZhZj2LzdyCysRKSvK355rhsSwrykLotsjEwmw9uD4hAX6omiCh3GLT4IrZ4Dlq/HsENE1Mh2nSvAyK//RGG5FglhKvzyXDdE+rpJXRbZKGdHB8wd1RGezgqkZhbj3TUnpC7J6jDsEBE1os2n8vHkgn2o0BpwTws//PB0V/i5K6Uui2xcuI8rZg9vBwBYtPsiVqZmS1uQlWHYISJqJL8dzcEz3+9Hld6IlNhAzB+dBDfeWk715IGYQDz/QAsAwORfjuLc5TKJK7Iekoadbdu2YcCAAQgJCYFMJsOKFSss9gshMHXqVAQHB8PFxQUpKSk4c+aMxTGFhYUYNWoUPD094eXlhbFjx6KsjBeYiKzLLweyqieAExiYGIK5j3XgGldU7yamtEK35r6o1Bkw6cdUTjhYTdKwU15ejsTERHz++ec33T9r1ix88sknmDdvHvbs2QM3Nzf06dMHGo3GfMyoUaNw/PhxrF+/HqtXr8a2bdvwzDPPNNYpEBHd0fd/XsS/fzoMowD+nhSOj//eDo4ObFin+ucgl+E/jyZC5eKII1klmLPhzJ2f1ATIhJVMuyiTybB8+XIMHjwYgKlVJyQkBP/+97/x4osvAgBKSkoQGBiIhQsXYvjw4Th58iTatGmDffv2ISkpCQCwbt069OvXD1lZWQgJCanRn61Wq6FSqVBSUgJPT88GOT8iapq+2nYO7/12CgAwplsUpv6tDeScLJAa2JojORi35CDkMuDHfySjU5R9rq1W089vq/3VIj09Hbm5uUhJSTFvU6lU6NKlC3bv3g0A2L17N7y8vMxBBwBSUlIgl8uxZ8+eW752VVUV1Gq1xYOIqL7N3fJX0Bl/fwtMG8CgQ42jf0IwhnYwzb8z6cdUlGqa9oKhVht2cnNzAQCBgYEW2wMDA837cnNzERAQYLFfoVDAx8fHfMzNzJgxAyqVyvwIDw+v5+qJqKn7ett5zFxnCjr/7t0KL/ZpDZmMQYcaz5sD2yDcxwVZRZWYtuq41OVIymrDTkOaMmUKSkpKzI/MzEypSyIiO/LNjnS8+9tJAMCklFZ4vldLiSuipsjD2REfP9oOchmw7GA2Np3Kk7okyVht2AkKCgIA5OVZXpy8vDzzvqCgIOTn51vs1+v1KCwsNB9zM0qlEp6enhYPIqL6sHBnOt5ZbZrU7V+9WmJCCoMOSScpygdj74kGALy+/BjKqvQSVyQNqw070dHRCAoKwsaNG83b1Go19uzZg+TkZABAcnIyiouLceDAAfMxmzZtgtFoRJcuXRq9ZiJq2r7ffQFv/moKOuPub45JDDpkBV7o3RoRPq64VKLBrOqu1aZG0rBTVlaG1NRUpKamAjANSk5NTUVGRgZkMhkmTpyI6dOnY9WqVTh69CieeOIJhISEmO/Yio2NRd++ffH0009j79692LlzJ8aPH4/hw4fX+E4sIqL6sGRPBt5YaRoX8Y8ezfDigxyjQ9bBxckBM4bEAzBNg7DvQqHEFTU+ScPO/v370b59e7Rv3x4A8MILL6B9+/aYOnUqAODll1/G888/j2eeeQadOnVCWVkZ1q1bB2dnZ/NrLF68GDExMejVqxf69euHe+65B1999ZUk50NETdPyQ1l4dflRAMDT90Zjct8YBh2yKt1b+OHRpDAIAbzyyxFodAapS2pUVjPPjpQ4zw4R1dX6E3l49r8HYDAKjE6OxJsD2zLokFUqqdAh5eOtuFxahQm9WmJS71ZSl3TXbH6eHSIia7f73BWMW3IQBqPAkA6hmDaAQYesl8rVEdMGtAEAzN16DpmFFRJX1HgYdoiI6uBIVjH+b9E+aPVG9G4TiFlDEzhhIFm9/vHB6NbcF1q90XzXYFPAsENEVEtn8kox+tu9KNcakNzMF5+OaA8F17oiGyCTyfDWwLZQyGX440QetqTl3/lJdoD/O4mIaiGzsAKPf7MXRRU6JIap8PXoJK5eTjalZaAHRneLAgC89esJVOntf7Ayww4RUQ1dKavCE9/uRa5ag5YB7lj4ZGe4KxVSl0VUaxNSWsLPXYn0gnIs2HlB6nIaHMMOEVENVGoNGLtoP9ILyhHq5YLvx3aBt5uT1GUR1YmnsyMmPxQDAPh801kUlmslrqhhMewQEd2B3mDE8z8cRGpmMbxcHbHoqc4IUjnf+YlEVmxI+1DEBnuitEqPzzadlbqcBsWwQ0R0G0IITF11HBtO5kOpkGP+E0loEeAudVlEd00ul2FKdevO939esOtb0Rl2iIhu4/PNZ7FkTwZkMmDO8PZIivKRuiSienNfK3/c29IPOoPAB7+nSV1Og2HYISK6hZ8PZOHDP04DAN4c0BZ944Ikroio/r3SNwYyGbDq8CUcySqWupwGwbBDRHQT289cxuRfjgAwLex59VZdInsTF6rCw+1CAQCz1tln6w7DDhHRdc7ml+Kfiw9CbxQY1C4Er/SJkbokogY1qXcrODrIsONsgV2uis6wQ0R0jcJyLZ5auB+lGj2SIr0xaxiXgSD7F+7jimEdwwEAczackbia+sewQ0RUrUpvwLPfH0BGYQXCfVzw5eMdoVRwdmRqGsbd3xwKual1Z7+dte4w7BARwXSL+avLjmHvhUJ4KBX4dnQn+LorpS6LqNGEebvikaQwAMCcjfbVusOwQ0QEYO7Wc/jlYBYc5DJ8NqoDWgZ6SF0SUaP7Z88WUMhl2H7Gvlp3GHaIqMlbdyzHfBfKtAFt0KOVv8QVEUkj3Mc+W3cYdoioSTuaVYKJP6YCAEYnR+KJ5ChJ6yGS2j97toBDdevO8UslUpdTLxh2iKjJyldr8H/f7YNGZ8R9rfzxxt/aSF0SkeTCfVzRLz4YADB/e7rE1dQPhh0iapKq9AY8+98DyFNXoUWAOz4b2R4KB74lEgHA0/dGAwB+PXwJl4orJa7m7vF/NhE1SW+uOoGDGcXwcFbg6yeS4OnsKHVJRFYjIcwLXZv5QG8UWLjrgtTl3DWF1AUQETW2xXsu4oe9psU9PxnRHtF+blKXRGR1/nFfc/x5vhBL9lxEO8dclBTkIjg4GPfeey8cHGxr/im27BBRk7L/QiHeXHUcAPDig61xf+sAiSsisk49Wvkj0MWIsioDRrz+KUaOHIn7778fUVFRWLZsmdTl1QrDDhE1GbklGjy3+CB0BoF+8UH4Z8/mUpdEZLVWrFiOE8s+BQB4dBwAyEyRITs7G8OGDbOpwMOwQ0RNwtUByZdLqxAT5IEPhiVCJuOaV0Q3YzAYMGHCBJQd3wJDRQkUngFwaZ4EwDTbOABMnDgRBoNByjJrjGGHiOyeEAJvrDiG1MxiqFwc8eXjHeGm5JBFolvZvn07srKyAIMOZUc3AAA82vUz7xdCIDMzE9u3b5eqxFph2CEiu7dkbwb+tz8Lchnw6Yj2iPTlgGSi28nJyTF/X5a6DgDg3KwDFKrAWx5nzRh2iMiuHckqxlurTgAAXuoTg/u4FATRHQUHB5u/1xfnoDL9IGQyOdzb9b3lcdaMYYeI7FZxhRbP/fcgtAYjercJxLM9mkldEpFNuPfeexEWFmYe11Z6aC0AwD2+N+CggEwmQ3h4OO69914py6wxhh0isktGo8CkH1ORXVyJSF9XfPgIByQT1ZSDgwPmzJkDAJDJZKg8uwf60gI4uHnBrXV3AMDs2bNtZr4dhh0isktzt57D5rTLUCrk+GJUB6hcOEMyUW0MGTIEP//8M0JDQwFhRNnh3wEAvp0G4ueff8aQIUMkrrDmGHaIyO7sPFuA//yRBgB4Z1Ac2oaoJK6IyDYNGTIEFy5cwObNm/He0wMBALKg1uh8f987PNO6MOwQkV3JLdHgXz8cglEAjyaF4dFO4VKXRGTTHBwc0LNnT/xz9HAkN/MFAKw4lC1xVbXDsENEdkNnMGL8koO4Uq5FbLAn3h4UJ3VJRHZlaMcwAMAvB7PNkwvaAoYdIrIbH/6Rhv0Xi+ChVGDuqA5wdrSNwZNEtuKhuCC4OjkgvaAcBzOKpS6nxhh2iMgubEnLx5dbzwMAPngkAVFcyZyo3rkpFegbFwQA+OVglsTV1BzDDhHZvDy1Bv/+32EAwONdI9E3zjYmOiOyRUPam7qy1h7Ngd5glLiammHYISKbZqieT+dKuRYxQR54rX+s1CUR2bWuzXzg4+aEogoddp+/InU5NcKwQ0Q27YvNZ7Hr3BW4ODrgs5Ecp0PU0BQOcnNX1m9HuTYWEVGD2pteiI83nAYAvDM4Di0C3CWuiKhp6B9v6ipedywXOhvoymLYISKbVFSuxYSlpvl0hrQPxbDqW2KJqOF1ifaBb3VX1p820JXFsENENkcIgZd+PoKcEg2i/dzw9mDOp0PUmBQOcvSp7spac8T6u7IYdojI5izYeQEbTubByUGOT0e0h7tSIXVJRE3O1a6s34/nWv1dWQw7RGRTjl8qwftrTwEAXu0Xg7hQrntFJIUu0T7wcnVEUYXO6icYZNghIptRqTVgwtJUaA1GpMQGYnS3KKlLImqyFA5y9GjlDwDYdCpf4mpuj2GHiGzGjLUncTa/DP4eSswcGg+ZTCZ1SURN2gMxAQCATafyJK7k9hh2iMgmbD6Vj+92XwQAfPhIInzdlRJXREQ9WvlDLgNO55Uhq6hC6nJuyarDjsFgwBtvvIHo6Gi4uLigefPmeOeddyxWWhVCYOrUqQgODoaLiwtSUlJw5swZCasmovpWUFaFl342LQfxZPcoc9M5EUnLy9UJSZE+AEy/kFgrqw47M2fOxNy5c/HZZ5/h5MmTmDlzJmbNmoVPP/3UfMysWbPwySefYN68edizZw/c3NzQp08faDQaCSsnovoihMArPx9BQZkWrQM98ErfGKlLIqJr3F/dlbWRYadudu3ahUGDBqF///6IiorCsGHD8OCDD2Lv3r0ATG+Cs2fPxuuvv45BgwYhISEB3333HS5duoQVK1ZIWzwR1YvFezKw8VQ+nBzkmD28HZeDILIyV8ft7Dp3BRqdQeJqbs6qw063bt2wceNGnD5tmg7+8OHD2LFjBx566CEAQHp6OnJzc5GSkmJ+jkqlQpcuXbB79+5bvm5VVRXUarXFg4isz9n8MkxfcwIA8MpDMYgN9pS4IiK6XqtAd/h7KKHVG5GaWSx1OTdl1WFn8uTJGD58OGJiYuDo6Ij27dtj4sSJGDVqFAAgNzcXABAYGGjxvMDAQPO+m5kxYwZUKpX5ER4e3nAnQUR1otUbMfHHQ9DojLi3pR+e5G3mRFZJJpOhazNfALDapSOsOuz873//w+LFi7FkyRIcPHgQixYtwocffohFixbd1etOmTIFJSUl5kdmZmY9VUxE9WX2htM4lq2Gt6sjPnwkEXI5bzMnslZdok2DlPecL5S4kpuz6jnWX3rpJXPrDgDEx8fj4sWLmDFjBkaPHo2gINO6HHl5eQgODjY/Ly8vD+3atbvl6yqVSiiVvG2VyFoduFiEeVvPAQBmDIlHoKezxBUR0e1cbdk5mFEEjc5gdWPrrLplp6KiAnK5ZYkODg4wGk1rcERHRyMoKAgbN24071er1dizZw+Sk5MbtVYiqh8VWj1e/OmweTXzvnHBd34SEUmqub8b/NyVqNIbcdgKx+1YddgZMGAA3n33XaxZswYXLlzA8uXL8dFHH+Hhhx8GYOonnDhxIqZPn45Vq1bh6NGjeOKJJxASEoLBgwdLWzwR1cmsdWlILyhHsMoZ0wa2lbocIqoBmUyGLs2qu7LSra8ry6q7sT799FO88cYb+Oc//4n8/HyEhITgH//4B6ZOnWo+5uWXX0Z5eTmeeeYZFBcX45577sG6devg7MxmbyJbs/NsARbuugAAmDk0ASoXR2kLIqIa69rMF2uO5ODP81fwr14tpS7HgkxcOx1xE6VWq6FSqVBSUgJPT97aSiQFtUaHvh9vw6USDR7rGoHpg+OlLomIauFsfilSPtoGZ0c5Dk97EEpFw4/bqennt1V3YxFR0/HOrydwqUSDSF9XTHkoVupyiKiWmvu7w8/dCRqdEUeySqQuxwLDDhFJbsOJPPx0IAsymWmRTzelVfewE9FNyGQydIk23ZW118rG7TDsEJGkCsu1mLzsKADgmXuboVOUj8QVEVFdtY/wAgCrm0mZYYeIJCOEwOsrjqKgrAqtAt0xqXcrqUsioruQEOYFADjKbiwiIpNfj+Tgt6O5UMhl+OhRLvJJZOviQj0hlwG5ag3y1RqpyzFj2CEiSVwurcLUlccAAOMfaIG4UJXEFRHR3XJ1UqBlgAcA4LAVte4w7BCRJKatOobiCh3aBHti3P0tpC6HiOpJfJjpF5cjWcXSFnINhh0ianS/Hf2r+2rWsAQ4OvCtiMheJFaHnaPZbNkhoiaqsFxr7r56rmdzdl8R2ZnYYNPkfmm5pRJX8heGHSJqVG/9ehwFZVq0CnTH+AfYfUVkb1oGmsbs5JRoUFKpk7gaE4YdImo060/kYWXqJchlwAfDEhtlOnkialwqF0eEqEzrU57Os47WnRpNU/rCCy/U+oVff/11+PhwcjAiMimp0OG15abJA5++rxkSw72kLYiIGkzrIA9cKtHgVG6pVUwUWqOwM3v2bCQnJ8PJyalGL7pjxw6MHz+eYYeIzN5ZcwL5pVVo5u+GSSmcPJDInrUK8sDmtMs4bSXjdmq8AM3y5csREBBQo2M9PDzqXBAR2Z/Nafn4uXrtqw+GJXDyQCI717p63I61dGPVaMzOggULoFLV/I6JL7/8EoGBgXUuiojsh1qjw6vVa1891T0aHSPZ4ktk76L83AAAGYUVEldiUqOwM3r0aCiVShgMBmzbtg3FxcW3PX7kyJFwc3Orj/qIyMbN+O0Ucko0iPR1xYsPtpa6HCJqBFG+pgyQU6KBRmeQuJpa3o3l4OCABx98EEVFRQ1VDxHZkd3nruCHvRkAgJlDE+DixO4roqbA29URHkrTSJlMK2jdqfWt53FxcTh//nxD1EJEdkSjM+DV6ruvRnWJQNdmvhJXRESNRSaTIdLPFQBw4YoNhp3p06fjxRdfxOrVq5GTkwO1Wm3xICICgE83nUF6QTkCPZV45aEYqcshokYW6WPqyrp4pVziSmpxN9ZV/fr1AwAMHDgQMpnMvF0IAZlMBoNB+r45IpLWyRw1vtxqagF+e1AcPJ0dJa6IiBpbpK+pZeeiFbTs1DrsbN68uSHqICI7YTAKTF52FHqjQN+2QejTNkjqkohIAlfDzgVbbNnp0aNHQ9RBRHZi0a4LOJxZDA9nBd4a1FbqcohIIuHeprCTU6KRuJI6hJ2rKioqkJGRAa1Wa7E9ISHhrosiItuUVVSBD/9IAwBMfigGgZ7OEldERFIJ8FQCAPLUNhh2Ll++jCeffBJr16696X6O2SFqmoQQeGPFMVRoDegc5YMRnSKkLomIJBRQ/ctOqUaPSq1B0qknan031sSJE1FcXIw9e/bAxcUF69atw6JFi9CyZUusWrWqIWokIhvw65EcbE67DCcHOd4bEg+5XHbnJxGR3fJQKuBSvTRMfqm0rTu1btnZtGkTVq5ciaSkJMjlckRGRqJ3797w9PTEjBkz0L9//4aok4isWHGFFm//ehwAMO7+FmgR4C5xRUQkNZlMhgBPJS5eqUCeugqRvtKtrFDrlp3y8nLzgqDe3t64fPkyACA+Ph4HDx6s3+qIyCa8u+YkCsq0aBngjud6Npe6HCKyEoEepq4sqcft1DrstG7dGmlppgGIiYmJ+PLLL5GdnY158+YhODi43gskIuu282wBfqpe0fz9oQlwUtT6bYWI7NTVQcr5pVWS1lHrbqwJEyYgJycHADBt2jT07dsXixcvhpOTExYuXFjf9RGRFbt2SYjHu0aiY6S3xBURkTUJqG7ZsbkxO4899pj5+44dO+LixYs4deoUIiIi4OfnV6/FEZF1+2zTWVy8UoEgT2e81IcrmhORJU8XU8wo1eglraPO8+wApltNXVxc0KFDh/qqh4hsxNn8Uny57RwA4M2BbeHBJSGI6Dru1Sufl0kcdurUuf7NN98gLi4Ozs7OcHZ2RlxcHObPn1/ftRGRlRJC4LXlx6AzCPSKCUCftoFSl0REVujqunhlVTbWsjN16lR89NFHeP7555GcnAwA2L17NyZNmoSMjAy8/fbb9V4kEVmXXw5mY096IVwcHfDWoLYWiwITEV3l7ny1G0snaR21Djtz587F119/jREjRpi3DRw4EAkJCXj++ecZdojsXFG5Fu/9dhIAMCGlJcKq178hIrre1W4sqcfs1LobS6fTISkp6YbtHTt2hF4v7ckQUcN7f+0pFJZr0TrQA2PviZa6HCKyYldbdqTuxqp12Hn88ccxd+7cG7Z/9dVXGDVqVL0URUTWad+FQvy4PxMA8O7DcXB04Jw6RHRrns7W0bJTo26sF154wfy9TCbD/Pnz8ccff6Br164AgD179iAjIwNPPPFEw1RJRJLTGYx4rXpOneGdwpEU5SNxRURk7dyVfw1QFkJINr6vRmHn0KFDFj937NgRAHDunOm2Uz8/P/j5+eH48eP1XB4RWYv529NxOq8MPm5OmPxQjNTlEJENuNqNZTAKaHRGyVY+r1HY2bx5c0PXQURWLLOwAnM2ngYAvNYvFl6uThJXRES2wNXxr3BTodVLFnbY4U5EtyWEwLRVx6HRGdG1mQ+GdAiVuiQishFy+V/dVkYhYR01OWjIkCFQq9U1ftFRo0YhPz+/zkURkfX4/XguNp3Kh6ODDNMHx3NOHSKqlat5Rwjp0k6NurFWrlyJy5cv1+gFhRD49ddf8c477yAgIOCuiiMiaZVV6fHmqhMAgGd7NEeLAHeJKyIiWyOXyWAUAhI27NQs7Agh0KpVq4auhYiszEd/nEauWoNIX1eMu7+F1OUQkQ2Sy2QABIzW3rJTlwHKoaHs1yeyZceyS7BwVzoA4O1BcXB2lGZgIRHZuOpuLCnH7NQo7PTo0aOh6yAiK2I0Cry+4hiMAvhbQjB6tPKXuiQislFXx+wYJUw7vBuLiG7w04FMpGYWw12pwBt/ayN1OURkw+TVNzVI2IvFsENElorKtXh/7SkAwMSUlgj0dJa4IiKyZeawI+EQZasPO9nZ2Xjsscfg6+sLFxcXxMfHY//+/eb9QghMnToVwcHBcHFxQUpKCs6cOSNhxUS2bdbvaSiq0KF1oAdGd4uSuhwisnFXJ6uw+nl2pFJUVITu3bvD0dERa9euxYkTJ/Cf//wH3t7e5mNmzZqFTz75BPPmzcOePXvg5uaGPn36QKPRSFg5kW1KzSzG0n0ZAIB3BnOhTyK6e1UGIwBAqZDu/aRGA5Svp9frsWXLFpw7dw4jR46Eh4cHLl26BE9PT7i71988HDNnzkR4eDgWLFhg3hYdHW3+XgiB2bNn4/XXX8egQYMAAN999x0CAwOxYsUKDB8+vN5qIbJ3BqPAGyuOQQhgSPtQdI7mQp9EdHcMRgGt3hR2XCS8o7PWMevixYuIj4/HoEGDMG7cOPNkgzNnzsSLL75Yr8WtWrUKSUlJeOSRRxAQEID27dvj66+/Nu9PT09Hbm4uUlJSzNtUKhW6dOmC3bt33/J1q6qqoFarLR5ETd0PezNwNLsEHs4KTOkXK3U5RGQHKnUG8/dSrYsF1CHsTJgwAUlJSSgqKoKLi4t5+8MPP4yNGzfWa3Hnz5/H3Llz0bJlS/z+++947rnn8K9//QuLFi0CAOTm5gIAAgMDLZ4XGBho3nczM2bMgEqlMj/Cw8PrtW4iW3OlrAof/J4GAHjxwdbw91BKXBER2YNK7V9hx6a6sbZv345du3bBycly1eOoqChkZ2fXW2EAYDQakZSUhPfeew8A0L59exw7dgzz5s3D6NGj6/y6U6ZMwQsvvGD+Wa1WM/BQk/b+2lMoqdShTbAnRnWJkLocIrITV8OOi6ODpOvq1TpmGY1GGAyGG7ZnZWXBw8OjXoq6Kjg4GG3aWM7xERsbi4wM0wDKoKAgAEBeXp7FMXl5eeZ9N6NUKuHp6WnxIGqqDlwsxE8HsgCYBiUrOCiZiOrJ1W4sVwm7sIA6hJ0HH3wQs2fPNv8sk8lQVlaGadOmoV+/fvVZG7p37460tDSLbadPn0ZkZCQA02DloKAgi+4ztVqNPXv2IDk5uV5rIbJHeoMRr684DgB4NCkMHSO97/AMIqKauxp2pF5uptbdWB9++CH69u2LNm3aQKPRYOTIkThz5gz8/Pzwww8/1GtxkyZNQrdu3fDee+/h0Ucfxd69e/HVV1/hq6++AmAKWhMnTsT06dPRsmVLREdH44033kBISAgGDx5cr7UQ2aPv/7yIkzlqqFwc8UrfGKnLISI7U6HVA5B2cDJQh7ATHh6Ow4cP48cff8Thw4dRVlaGsWPHYtSoURYDlutDp06dsHz5ckyZMgVvv/02oqOjMXv2bIwaNcp8zMsvv4zy8nI888wzKC4uxj333IN169bB2ZmzvhLdTn6pBh/9cRoA8HLf1vB156BkIqpfGt1fY3akJBOi5qtV6HQ6xMTEYPXq1YiNtZ9bU9VqNVQqFUpKSjh+h5qMST+mYvmhbCSEqbD8n93hIJdu8CAR2adfDmTh3z8dxr0t/fD92C71/vo1/fyu1ZgdR0dHzkxMZAf2nL+C5YeyIZMB7wyKY9AhogZxpbwKAODr5nSHIxtWrQcojxs3DjNnzoRer2+IeoiogekMRryx8hgAYETnCCSGe0lbEBHZrStlWgCQvJu81mN29u3bh40bN+KPP/5AfHw83NzcLPYvW7as3oojovq3cOcFnM4rg4+bE17u01rqcojIjhWYw460LTu1DjteXl4YOnRoQ9RCRA0st0SD2RtMg5In942Bl6u0b0BEZN+udmP5udlYy861i3ISkW2ZvuYEyrUGdIjwwrCOYVKXQ0R27oqVtOxwqlSiJmLX2QKsPpIDuQx4e1Ac5ByUTEQN7EpZ9QBlWxuzEx0dfdv1Lc6fP39XBRFR/dPqjZi6yjRT8uNdIxEXqpK4IiKyd0IIFJRXt+xIfDdWrcPOxIkTLX7W6XQ4dOgQ1q1bh5deeqm+6iKierRgZzrO5pfB180JLzzIQclE1PBKq/TQ6o0ApO/GqnXYmTBhwk23f/7559i/f/9dF0RE9SunpBJzNp4BAEx+KAYqF0eJKyKipiCrsBIA4O3qCFenWseNelVvY3Yeeugh/PLLL/X1ckRUT95dcxIVWgM6RnpjaAcOSiaixpFRWAEAiPBxlbiSegw7P//8M3x8fOrr5YioHuy0GJTcloOSiajRZBWZwk64FYSdWrcrtW/f3mKAshACubm5uHz5Mr744ot6LY6I6k6rN2Jq9UzJj3eNRNsQDkomosZztWXHJsPO4MGDLX6Wy+Xw9/dHz549ERMTU191EdFd+nZnOs5dLuegZCKSRKYVdWPVOuxMmzatIeogonqUU1KJT6oHJU/pF8tByUTU6Gx6zM7Bgwdx9OhR888rV67E4MGD8eqrr0Kr1dZrcURUN9OrByUnRXpjSPtQqcshoibGaBTILDLdjRXubYNh5x//+AdOnzatrXP+/Hn8/e9/h6urK3766Se8/PLL9V4gEdXOjjMFWMOZkolIQpfLqqDVG+EglyHYy1nqcmofdk6fPo127doBAH766Sf06NEDS5YswcKFC3nrOZHETDMlmwYlP5EchTYhnhJXRERN0fnL5QCAUC8XODpIvzJVrSsQQsBoNM2IuGHDBvTr1w8AEB4ejoKCgvqtjohq5Zsd6Th/uRx+7k6Y1LuV1OUQUROVlqsGALQK9JC4EpNah52kpCRMnz4d33//PbZu3Yr+/fsDANLT0xEYGFjvBRJRzVwqrsSnm6oHJT/EQclEJJ20vFIAQEyQjYad2bNn4+DBgxg/fjxee+01tGjRAoBpUsFu3brVe4FEVDNXZ0ruFOWNIR04KJmIpHMyxxR2WltJ2Kn1recJCQkWd2Nd9cEHH8DBwaFeiiKi2tl+5jLWHDUNSn5rYJzFxJ9ERI3JaBQ4Xd2yExtsHWGn1i07mZmZyMrKMv+8d+9eTJw4Ed999x0cHdlsTtTYtHojpq06DoCDkolIellFlajQGuCkkCPK103qcgDUIeyMHDkSmzdvBgDk5uaid+/e2Lt3L1577TW8/fbb9V4gEd3eX4OSlRyUTESSO1U9OLmFvzsUVnAnFlCHsHPs2DF07twZAPC///0PcXFx2LVrFxYvXoyFCxfWd31EdBuXiv+aKfnVfjEclExEkkvLta7ByUAdwo5Op4NSqQRguvV84MCBAICYmBjk5OTUb3VEdFvT15xApc40KPlhzpRMRFbgRI6pZcdaBicDdQg7bdu2xbx587B9+3asX78effv2BQBcunQJvr6+9V4gEd3c9jOX8dvRXDjIZXh7EAclE5F1SM0sBgAkhntJWse1ah12Zs6ciS+//BI9e/bEiBEjkJiYCABYtWqVuXuLiBpWld6AaSuvDkqORGwwByUTkfRySiqRU6KBg1yGhDCV1OWY1frW8549e6KgoABqtRre3t7m7c888wxcXaVf7IuoKfhmRzrOF3BQMhFZl0MZxQBM43VcnWodMRpMnYZJCyFw4MABfPnllygtNQ1EcnJyYtghagTZxZX4dONZAKZByZ7OHJRMRNbh4MUiAED7CC9pC7lOrWPXxYsX0bdvX2RkZKCqqgq9e/eGh4cHZs6ciaqqKsybN68h6iSiatNXc1AyEVmngxmmsNMhwvsORzauWrfsTJgwAUlJSSgqKoKLi4t5+8MPP4yNGzfWa3FEZGnb6ctYe4yDkonI+lTpDTh2yXQnlrWFnVq37Gzfvh27du2Ck5OTxfaoqChkZ2fXW2FEZKlKb8Cb1TMlj06O4qBkIrIqJy6podUb4ePmhEhf6xrWUuuWHaPRCIPBcMP2rKwseHhYzz31RPZm/va/BiVP7N1S6nKIiCwcuDpeJ9zL6lqdax12HnzwQcyePdv8s0wmQ1lZGaZNm4Z+/frVZ21EVC27uBKfbjLNlPxafw5KJiLrs/NsAQCgazPrm3Ov1t1YH374Ifr27Ys2bdpAo9Fg5MiROHPmDPz8/PDDDz80RI1ETd701Seg0RnROcoHg9txUDIRWRet3og96YUAgO4t/CSu5ka1Djvh4eE4fPgwfvzxRxw+fBhlZWUYO3YsRo0aZTFgmYjqx9ZrByUPbmt1zcNERKmZxajQGuDr5mRVa2JdVauwo9PpEBMTg9WrV2PUqFEYNWpUQ9VFRLhxUHJMEAclE5H12VHdhdWthR/kcuv7haxWY3YcHR2h0WgaqhYius787elILyiHvwcHJROR9bo6XueeFtY3XgeowwDlcePGYebMmdDr9Q1RDxFVyyqq+GtQcr9YDkomIqtUqtGZF/+0xvE6QB3G7Ozbtw8bN27EH3/8gfj4eLi5uVnsX7ZsWb0VR9SUTV990jQoOdoHg9qFSF0OEdFN7TlfCINRIMrXFWHe1jW/zlW1DjteXl4YOnRoQ9RCRNW2pOVj3XHToOR3OFMyEVmxzWn5AKy3VQeoQ9hZsGBBQ9RBRNWuHZQ8plsUWlvhnQ1ERABgNApsOJkHAEhpEyhxNbdW4zE7RqMRM2fORPfu3dGpUydMnjwZlZWVDVkbUZP09bbzuHClwjQoOYWDkonIeh3JLkGeugpuTg7o1tw6BycDtQg77777Ll599VW4u7sjNDQUc+bMwbhx4xqyNqImJ6uoAp9tPgsAeL1/LDw4KJmIrNj6E7kAgJ6tA6BUOEhcza3VOOx89913+OKLL/D7779jxYoV+PXXX7F48WIYjcaGrI+oSXmneqbkLtE+GJjIQclEZN3+OG7qwnqwrfV2YQG1CDsZGRkWa1+lpKRAJpPh0qVLDVIYUVOzOS0fvx/PM82UzEHJRGTl0gvKcSa/DAq5DD1bB0hdzm3VOOzo9Xo4OztbbHN0dIROp6v3ooiamiq9AW9VD0p+koOSicgGXO3C6trMFyoX6+5yr/HdWEIIjBkzBkql0rxNo9Hg2WeftZhrh/PsENXe1UHJAR5KTOCgZCKyAb9Xd2H1tuK7sK6qccvO6NGjERAQAJVKZX489thjCAkJsdjWkN5//33IZDJMnDjRvE2j0WDcuHHw9fWFu7s7hg4diry8vAatg6g+ZRb+NSj5NQ5KJiIbkFlYgQMXiyCTAX3jgqQu545q3LIj9fw6+/btw5dffomEhASL7ZMmTcKaNWvw008/QaVSYfz48RgyZAh27twpUaVEtcNByURka1amZgMAujX3RaCn8x2Oll6t18aSQllZGUaNGoWvv/4a3t7e5u0lJSX45ptv8NFHH+GBBx5Ax44dsWDBAuzatQt//vmnhBUT1cymU3n444RpUPI7gzkomYisnxACK1JNNycNahcqcTU1YxNhZ9y4cejfvz9SUlIsth84cAA6nc5ie0xMDCIiIrB79+5bvl5VVRXUarXFg6ixVWj1eGOFaVDy2Hui0SqQg5KJyPodv6TG2fwyOCnkNtGFBdRhuYjGtnTpUhw8eBD79u27YV9ubi6cnJzg5eVlsT0wMBC5ubm3fM0ZM2bgrbfequ9SiWplzoYzyC6uRKiXC2dKJiKbcbULKyU2AJ42MsbQqlt2MjMzMWHCBCxevPiG297vxpQpU1BSUmJ+ZGZm1ttrE9XEiUtqzN+RDgB4e1BbuDpZ/e8dREQwGAVWHbatLizAysPOgQMHkJ+fjw4dOkChUEChUGDr1q345JNPoFAoEBgYCK1Wi+LiYovn5eXlISjo1k1rSqUSnp6eFg+ixmIwCry6/CgMRoGH4oLQK9b6b9skIgKAP89fQZ66Cp7OCvRs7S91OTVm1b9O9urVC0ePHrXY9uSTTyImJgavvPIKwsPD4ejoiI0bN2Lo0KEAgLS0NGRkZCA5OVmKkonuaMneDKRmFsNdqcC0AW2lLoeIqMaW7M0AAAxIDLHqtbCuZ9Vhx8PDA3FxcRbb3Nzc4Ovra94+duxYvPDCC/Dx8YGnpyeef/55JCcno2vXrlKUTHRb+WoNZq09BQB4qU9rBKms/5ZNIiIAuFxahT+Om8bDjuwSIXE1tWPVYacmPv74Y8jlcgwdOhRVVVXo06cPvvjiC6nLIrqpt1afQGmVHolhKjzWNVLqcoiIauynA5nQGQQSw73QNqRhJxGubzYXdrZs2WLxs7OzMz7//HN8/vnn0hREVEOb0/Kx5kgO5DLg3Yfj4SDnnDpEZBuMRoGle00384zqbFutOoCVD1AmsheVWgPeWHEMAPBU92jEhdrWb0VE1LTtOFuAjMIKeDgr8LfEYKnLqTWGHaJGMGfjGWQVVSJE5YxJvVtJXQ4RUa0s2WMamDykfahNTpXBsEPUwE7lqjF/+3kAwNuD4uCmtL03CiJqunJKKrHhpGmB7ZFdbHOsIcMOUQMyGgVeXXYUeqNAn7aBSGnDOXWIyLYs3HkBeqNAl2gftA6yzWVtGHaIGtD3f17EwYxiuDk54M2BnFOHiGxLWZXePLfO0/c2k7iaumPYIWogWUUVmLXONKfO5IdiEKxykbgiIqLa+XFfJko1ejTzd8MDMQFSl1NnDDtEDUAIgdeWH0O51oBOUd4YZaP93ETUdOkNRnxbvYbf/93TDHIbni6DYYeoAaxIzcbW05fhpJDj/aEJNv0mQURN09pjucguroSvmxOGdLCdRT9vhmGHqJ4VlFXhrV9PAAAm9GqJ5v7uEldERFQ7QgjzXaSPdY2Es6PtrIN1Mww7RPXsrV9PoLhCh9hgTzxzn+0O6COipmv7mQIcziqBUiHH48m23w3PsENUjzacyMOvhy9BLgNmDU2AowP/ixGRbRFCYPaG0wBMrTp+7kqJK7p7fCcmqidqjQ6vVy8J8fR9zRAfxiUhiMj2bD9TgIMZxVAq5PhHD/tonWbYIaonM9eeQq5agyhfV0xK4ZIQRGR7rm3VGdUlEgEezhJXVD8YdojqwZ7zV7C4eu2YGUMSbH4wHxE1TTvO/tWq86ydtOoADDtEd02jM2DysqMAgBGdI5Dc3FfiioiIas/UqnMGADCySwQCPO2jVQdg2CG6a//5Iw3pBeUI9FRiSr8YqcshIqqTTafyceBiEZwUcjzXo7nU5dQrhh2iu7DvQiHmV88wOmNIPDydHSWuiIio9vQGI2asNS1v82T3KLtq1QEYdojqrEKrx0s/HYYQwCMdw/BADFc0JyLb9L/9WTibXwYvV0f8s2cLqcupdww7RHU0a10aLlypQLDKGW8MaCN1OUREdVJepcdH6013YD3/QEuoXOyvhZphh6gOdp+7goW7LgAAZg5NYPcVEdmsr7efR0FZFSJ8XPF4V9ufLflmGHaIaqm8So+Xfj4MwHT31X2t/CWuiIiobvJLNfhqm2kNrJf7toaTwj5jgX2eFVEDeu+3k8gqqkSolwte6x8rdTlERHX2/tpTqNAakBjuhf7xwVKX02AYdohqYfuZy+bJAz8YlgB3pULiioiI6mZveiGWHcyGTAa8NbAtZDKZ1CU1GIYdohpSa3R45ecjAIAnkiPRrYWfxBUREdWN3mDE1JWmtfyGdwpHu3AvaQtqYAw7RDU0ffUJXCrRIMLHFZMf4uSBRGS7Fu2+iFO5pfBydcRLfez//Yxhh6gG1h3Lxf/2Z0EmM3VfuTqx+4qIbFO+WoOPq281f7lPDHzcnCSuqOEx7BDdQZ5agynLTN1Xz9zXDF2ace0rIrJd09ecRFmVHonhXhjeKVzqchoFww7RbRiNAi/+dBhFFTq0DfHEv3u3lrokIqI623AiD6sOX4JcBrwzqC3kcvsdlHwthh2i21i46wK2nymAUiHHnOHt7HYOCiKyfyWVOry24igA4Ol7myEhzEvaghoR37mJbiEttxTvrzMtjPda/1i0CPCQuCIiorp7d80J5KmrEO3nhkm9W0ldTqNi2CG6iSq9AROWHoJWb8T9rf3tdgp1Imoatp2+bL7JYtawBDg7OkhdUqNi2CG6iQ/WpeFUbil83Zwwa1iiXU+2RUT2raxKjynLTN1Xo5Oj0CnKR+KKGh/DDtF1dpwpwPwd6QBMvwH5eyglroiIqO6mrTyO7OJKhHm74KU+TfMmC4YdomsUlmvx759SAQCjukSgV2ygtAUREd2FXw9fwi8HsyCXAR892g5uTXSJG4YdompXbzPPU1ehmb8bXu/fRuqSiIjqLLu4Eq8uN3Vfjbu/BTpHN73uq6sYdoiqzd9xHptO5cNJIcdnIzrAxalpDeAjIvthMApMWpqKUo0e7cK98K9eLaUuSVIMO0QADmYUYda6NADAtAFt0CbEU+KKiIjqbu6Ws9h7oRBuTg6YM7wdHB2a9sd90z57IgAlFTo8v+QQ9EaB/gnBGNk5QuqSiIjqbPe5K/ioeu2rtwfFIdLXTeKKpMewQ02aEAIv/XwY2cWViPR1xftD4nmbORHZrDy1Bs//cBBGAQxpH4ohHUKlLskqMOxQk7Zw1wX8cSIPTg6mcToezo5Sl0REVCc6gxHjFh9EQZkWMUEeePdh/vJ2FcMONVlHsorx3m8nAQCv9otBfJhK4oqIiOpuxm+nsP9iETyUCsx7rCNvsrgGww41SSWVOoxfcgg6g0CftoEY3S1K6pKIiOrs18OX8O1O02So/3k0EVF+HKdzLYYdanKMRoFJP6Yio7ACoV4umDWUy0EQke06klWMl34+DAB4tkdzPNg2SOKKrA/DDjU5n2w6g02n8qFUyPHl4x2hcuU4HSKyTbklGjz93X5odKZFi5vqchB3wrBDTcrmU/mYs/EMAODdh+MRF8pxOkRkmyq1Bvzfd/uQp65CywB3fDKiPRzkbKW+GYYdajIuXinHhKWHIATwWNcIDOsYJnVJRER1YjQKvPC/VBzLVsPHzQnfjunEu0lvg2GHmoRKrQH/+P4A1Bo92kd4Yerf2kpdEhFRnc1cdwprj+XC0UGGLx/viHAfV6lLsmpWHXZmzJiBTp06wcPDAwEBARg8eDDS0tIsjtFoNBg3bhx8fX3h7u6OoUOHIi8vT6KKyRoJITB52RGcyi2Fn7sT5o7qCCeFVf/TJyK6pfnbz+PLbecBAO8PSUCnqKa7wGdNWfU7/tatWzFu3Dj8+eefWL9+PXQ6HR588EGUl5ebj5k0aRJ+/fVX/PTTT9i6dSsuXbqEIUOGSFg1WZsFOy9gZeolOMhl+HxkBwSpnKUuiYioTlamZmP6GtP8YC/3bY2h7I6vEZkQQkhdRE1dvnwZAQEB2Lp1K+677z6UlJTA398fS5YswbBhwwAAp06dQmxsLHbv3o2uXbvW6HXVajVUKhVKSkrg6ckFIO3JlrR8PLVwH4wCeONvbTD2nmipSyIiqpNtpy9j7KJ90BkEnuwehal/a9Pkp82o6ee3VbfsXK+kpAQA4ONjarI7cOAAdDodUlJSzMfExMQgIiICu3fvvuXrVFVVQa1WWzzI/pzNL8XzSw7BKIBHOobhqe5RUpdERFQnhzKK8Ox/D0BnEBiQGII3+jPo1IbNhB2j0YiJEyeie/fuiIuLAwDk5ubCyckJXl5eFscGBgYiNzf3lq81Y8YMqFQq8yM8PLwhSycJFJVr8dTC/Sit0qNzlA+mPxzHNwYisklHs0rwxLd7UaE14J4WfvjwkQTIeYt5rdhM2Bk3bhyOHTuGpUuX3vVrTZkyBSUlJeZHZmZmPVRI1kKrN+LZ/x5ARmEFwn1cMPexDlAquEYMEdmeE5fUePzbPSjV6NEpyhtfPt6R72d1oJC6gJoYP348Vq9ejW3btiEs7K/BWEFBQdBqtSguLrZo3cnLy0NQ0K2ny1YqlVAqlQ1ZMklECIGpK49hT3oh3JUKfDO6E3zdea2JyPacySvFY9/sQXGFDu3CvfDtmE5wU9rEx7bVseqWHSEExo8fj+XLl2PTpk2IjrYcXNqxY0c4Ojpi48aN5m1paWnIyMhAcnJyY5dLVuCbHelYui8Tchnw6cj2aBXoIXVJRES1diavFCPn70FhuRbxoSoseqozJw28C1YdEceNG4clS5Zg5cqV8PDwMI/DUalUcHFxgUqlwtixY/HCCy/Ax8cHnp6eeP7555GcnFzjO7HIfqw7loN3fzPdkvlqv1jc3zpA4oqIiGrvWLZpjE5huRaxwZ74fmxnqFwYdO6GVYeduXPnAgB69uxpsX3BggUYM2YMAODjjz+GXC7H0KFDUVVVhT59+uCLL75o5EpJavsuFOJfS1MhBDCqSwRvMScim3TgYhHGLNiLUo0eCWEqLHqyM7xcnaQuy+bZ1Dw7DYXz7Ni2s/mlGDp3N0oqdejdJhDzHuvIxfCIyObsOleA/1u0HxVaAzpFeeObMZ3gya6r26rp57dVt+wQ3UmeWoPR3+5DSaUO7SO88MlwrvpLRLbn9+O5+NcPh1ClN+Leln748vGOcHXiR3R94d8k2axSjQ5jFuxDdnElmvm54ZvRneDixFsyici2fL/7AqatOg6jAFJiA/HZyPZwduR7WX1i2CGbVKU34Nn/HsDJHDX83JVY9FRn+LixX5uIbIcQArN+T8PcLecAACM6R+CdQW2hcLDqG6VtEsMO2RydwYjxSw5h59krcHNywMInOyHcx1XqsoiIakyrN2LyL0ew7FA2AODfvVth/AMtONN7A2HYIZtiNAq89NNhrD+RByeFHF8/kYS4UJXUZRER1VhBWRWe++8B7LtQBAe5DO8PiccjSVy2qCEx7JDNEELg9ZXHsCL1EhRyGeaO6oBuLfykLouIqMZO5qjxf4v2I7u4Eh5KBT4b1QE9WvlLXZbdY9ghmyCEwHu/ncSSPRmQy4CP/94OvWIDpS6LiKjG1h3LxQv/S0WF1oAoX1fMH90JLQLcpS6rSWDYIZswe8MZfL09HQDw/pAEDEgMkbgiIqKaMRgF5mw4jU82nQUA3NPCD5+NbM/JAhsRww5ZNSEEZm84gzkbzwAApv6tDR7txL5tIrINV8qqMGFpKnacLQAAjOkWhdf7x/KOq0bGsENWSwiBj9afxqfVvw1NfigGT3EZCCKyEQcuFmHc4oPIVWvg4uiAGUPiMbh9qNRlNUkMO2SVrp9/4vX+sfi/e5tJXBUR0Z0ZjQLzd5zHrHVp0BsFmvm7Yd5jHdEq0EPq0poshh2yOkIIvL/2FL7cdh6AqeuKLTpEZAvy1Br8+3+Hzd1W/ROCMXNoAtyV/LiVEv/2yaoYjQJvrz6BhbsuAADeHtQWTyRHSVoTEVFNrD+Rh5d/PoyiCh1cHB0wdUAbDO8UzokCrQDDDlkNncGIF386jJWplwAA0wfH4bGukRJXRUR0e6UaHd777SR+2JsJAGgb4ok5w9vztnIrwrBDVqFSa8Bziw9gS9plKOQy/OfRRAxqx4F8RGTdtp+5jMm/HEV2cSUA4Ol7o/Fin9ZQKriQpzVh2CHJlVTo8NSifThwsQjOjnLMfawj7m8dIHVZRES3VFalN090CgDhPi6YNTQRyc19Ja6MboZhhySVXVyJpxbsQ1peKTydFVjwZCd0jPSRuiwiolv643gupq06jpwSDQDgieRIvNI3Bm4chGy1eGVIMkeyijF20X5cLq1CgIcS343tjJggT6nLIiK6qUvFlXhz1XH8cSIPABDh44r3h8ajW3Ou0WftGHZIEr8fz8WEpYeg0RkRE+SBb8Z0QqiXi9RlERHdQKs34rvdF/Dx+tMo1xqgkMvwjx7N8PwDLeHsyLE5toBhhxqVEALf7EjHu7+dhBBAj1b++Gxke3g4O0pdGhHRDbak5ePt1Sdw/nI5ACAp0hvvDYnnBIE2hmGHGo1GZ8DrK47h5wNZAIDHukbgzQFtuUYMEVmd9IJyvLP6BDadygcA+Lk74eU+MRjWMQxyOefNsTUMO9Qosooq8Nx/D+JodgnkMuDVfrEYe080J9siIqtSUFaFTzaewZI9GdAbBRRyGZ7sHoXne7WEJ1ugbRbDDjW4XWcLMP6HQygs18Lb1RGfjeyA7i04oI+IrEeFVo/529Px5dZzKNcaAAA9W/vj9f5tODmgHWDYoQYjhMDX28/j/bWnYBRAfKgKcx/rgDBvV6lLIyICYOpeX7InA19sOYeCsioApveqKQ/FoBt/KbMbDDvUIK6UVeGln4+Y+7uHdQzD9MFxvHOBiKxCld6ApXsz8cWWs8hTm0JOhI8rXurTGv3jgzkux84w7FC923m2AJN+TEV+aRWcFHK88bc2eKxLBMfnEJHkKrR6LN2bia+3nzdPChjq5YLxD7TA0A5hcFLwhgl7xLBD9UZnMOKj9acxb+s5CAG0DHDHpyPbc6JAIpJcSYUOi3ZfwIKd6Siq0AEAgjydMe6BFng0KYxrWdk5hh2qF2m5pXjxp8M4ml0CABjZJQJv9G8DFye+gRCRdC5eKceCnRfw0/5M88DjSF9X/OO+5hjSIZRd600Eww7dFb3BiC+3ncecDWegNRjh6azA+0MT0C8+WOrSiKiJEkJgb3ohvtmRjvUn8yCEaXtMkAf+eX8L9IsL4vxeTQzDDtXZ9a05vWIC8N6QeAR6OktcGRE1RSUVOiw7lIUf9mbgdF6ZeXuPVv4Ye0807m3px7GDTRTDDtVahVaPzzadxdfbz0NnEPB0VuDNgW3xcPtQvpEQUaMSQmD/xSL8sCcDa47moEpvBAA4O8rxcPswjL0nCi0CuLRDU8ewQzUmhMD6E3l469cTyC6uBACkxAbivYfjEMDWHCJqRMUVWiw7mI0f9mbgTP5frTgxQR4Y1SUCg9qHcsZjMmPYoRpJLyjH9NUnsLF63pxQLxdMG9AGvdsEsjWHiBqFRmfAplP5WJV6CZvS8qGtbsVxcXTAgMRgjOgcgXbhXnxPohsw7NBtXaleJ2Zx9Toxjg4yPHNfM4y/vyXvtCKiBqczGLHjbAF+Tb2EP07koaxKb94XG+yJkV0iMKhdCFtx6LYYduimKrUGfLszHXO3nDO/uTwQE4BX+8VynRgialB6gxEHLhbh1yOX8NvRXBSWa837Qr1c8LfEYAxMDEGbYE+24lCNMOyQBY3OgB/2ZmDe1nPmKdTjQj3xar9YdGvOdWKIqGGUVemx/fRlrD+Zh82n8s0T/wGAr5sT+ieYAk6HCG8u5UC1xrBDAEwhZ/GeDHy59RzyS00hJ9TLBS/1aY2BiSF8cyGiepdTUokNJ/Ox4UQedp+7Aq3BaN6ncnFESmwgBrYLQffmvpwXh+4Kw04TV1SuxZK9GVi46wIuXxNynuvZHI9wCnUiqkflVXrsvVCIHWcKsPNsAU7lllrsj/R1Re/YQKS0CURSpDcDDtUbhp0m6mx+Gb7dmY5lB7Og0Zl+mwr1csG4+1tgWEcuhkdEd09vMOJIdgl2nCnAjrMFOJRRBJ1BmPfLZECHCG+kxAaid5sANPd35xgcahAMO02IzmDExpP5+GFvBraevmzeHhfqibH3RKN/fAhDDhHVWXmVHqmZxdh3oRD7LxThUEaReT2qq8K8XXBvSz90b+GHbs394OPmJFG11JQw7DQB5y+X4cf9mfjlQBYKykx3NchkQO/YQIy9Jxqdo3342xQR1YoQApmFlTicVYyDGUXYf6EIJ3LUMBiFxXEqF0d0b+GL7i38cE8LP0T6uklUMTVlDDt2KrdEgzVHc/Dr4UtIzSw2b/f3UGJYxzAM7xTONx0iqhEhBPLUVTicVYyjWSWmr9klKL7mjqmrQr1c0CnKGx2jfNApyhutAjx4gwNJjmHHjlwqrsTGk3n49UgO9l0oNK/0K5cB97cOwN87heP+mAA4ctAfEd2CRmfA2fwynMotxakctelrrtrcKnwtJwc5YoM9kBjuhaQoHyRFeiPEy0WCqoluj2HHhhmMAqmZRdh4Mh+bTuXfcGdDUqQ3BiSG4KH4IAR4cO0qIvpLcYUW5y6X4/zlMpwvMH09d7kc6QXlN3RFAYCDXIaWAe5IDPNCfJgKiWFeaBXkzjs2ySYw7NgQo1EgLa8Uf56/gj/PX8Ge9EKLZmS5DGgf4Y2+bYPQPyGYv2ERNXEllTpkFVUgq6gS6dWB5vzlcpwvKLeYlfh6Xq6OiA3yREywh/lrywAPLhFDNothx4qVVOhwJLsYR7JKcDizGHsvFN7QR+7prECP1gF4IMYfPVoF8M4GoibCaBQorNAit0SDS8WVyCqqRGZ1sDE9KlCq0d/2NYJVzmjm74Zmfu6mr/7uiAnyQICHkjctkF1h2LECRqNAdnElTueV4kx+GU5cUuNIVjEuXKm44VhXJwckRfmgazMfdG3mi/hQFcfgENkRjc6AwnItCsu1uFKuRUFpFXLVGuSrNchTVyGvVIO8Eg3yS6ugv0l30/X83J0Q6u2KCB9XNPNzQzN/NzT3d0e0nxvclPwIoKbBbv6lf/755/jggw+Qm5uLxMREfPrpp+jcubPUZZkZjAJ5ao3pt6/CCmQWVSDjSgXO5JfhbH4ZKnWGmz4vwscVCdX94x2jvBluiGyEVm9EqUYHtUZv+lpZ/VWjQ3GFzhxm/vpahcIy7Q3z0tyOTGZaNyrUywVh3q4I83apfpi+D/V2gauT3bzNE9WZXfwv+PHHH/HCCy9g3rx56NKlC2bPno0+ffogLS0NAQEBktX1+oqjOJdfjpySSmQXV1rMHHo9Jwc5mvm7oWWgB1oFuCMh3AsJoSp4s1uKqEEJIaA1GKHRGVGpNaBSZ0CFVo9KrQEV1Q+N7ur31dt1BtOxWgPKtfprAs1f4ebqzOR1oZDL4OPmBB83J/i5KxHo6YwglelrgIczAj1N3/t7KPnLD1ENyIQQd24HtXJdunRBp06d8NlnnwEAjEYjwsPD8fzzz2Py5Ml3fL5arYZKpUJJSQk8PT3rra6Uj7bibH6Z+WeFXIYQLxeE+7gg3NsV4T6uaO7vjlaB7ojwceU6MFTviiu0KKvSQwiYHhDVXwGjENXTE5i+Gq/dL0z7gb+eZxSmYCCubrvu+6vPx3WvbzpGWPz5eqOA0SigNwoYrn2I6m0GIwwCMBiNMBgtv+qrjzMYqr8aBXQGAa3eCK3BCK3ecM33Rmiv7tMb/tpW/dAZhMXikw3BXamAh7MCns6O8HRRwMPZESoXR3OY8b361d0JPm5K+Lg5wdNZwTEzRDVQ089vm2/Z0Wq1OHDgAKZMmWLeJpfLkZKSgt27d9/0OVVVVaiqqjL/rFarG6S2Cb1awmAUCFI5I8zbBUGezgw01KhmrkvDD3szpC7DZjjIZXB1dICLU/XD0QGuTg5wdVLA2fy9aZ+reb8Cni6OFoHG09kRns6OcHdWwIET6hFJzubDTkFBAQwGAwIDAy22BwYG4tSpUzd9zowZM/DWW281eG0DEkMa/M8guh0nBxmUCjnkMhlkMkAGQHaL768eA9y4zeJYGSC77hhcc4z8mv24dts12x3kMijkMjiYH3Io5DLIZZbbFXIZ5Dcce912mQwKBzmcFHI4Xf2qkEOpkMPRwXKbxTHXbVdWb2OLCpH9sfmwUxdTpkzBCy+8YP5ZrVYjPDxcwoqIGsZbg+Lw1qA4qcsgIpKUzYcdPz8/ODg4IC8vz2J7Xl4egoKCbvocpVIJpVLZGOURERGRxGx+AImTkxM6duyIjRs3mrcZjUZs3LgRycnJElZGRERE1sDmW3YA4IUXXsDo0aORlJSEzp07Y/bs2SgvL8eTTz4pdWlEREQkMbsIO3//+99x+fJlTJ06Fbm5uWjXrh3WrVt3w6BlIiIianrsYp6du9VQ8+wQERFRw6np57fNj9khIiIiuh2GHSIiIrJrDDtERERk1xh2iIiIyK4x7BAREZFdY9ghIiIiu8awQ0RERHaNYYeIiIjsGsMOERER2TW7WC7ibl2dRFqtVktcCREREdXU1c/tOy0GwbADoLS0FAAQHh4ucSVERERUW6WlpVCpVLfcz7WxABiNRly6dAkeHh6QyWS1eq5arUZ4eDgyMzOb3LpaTfXcm+p5Azx3nnvTOvemet6A7Zy7EAKlpaUICQmBXH7rkTls2QEgl8sRFhZ2V6/h6elp1f8gGlJTPfemet4Az53n3rQ01fMGbOPcb9eicxUHKBMREZFdY9ghIiIiu8awc5eUSiWmTZsGpVIpdSmNrqmee1M9b4DnznNvWufeVM8bsL9z5wBlIiIismts2SEiIiK7xrBDREREdo1hh4iIiOwaww4RERHZNYadm9i2bRsGDBiAkJAQyGQyrFixwmK/EAJTp05FcHAwXFxckJKSgjNnzlgcU1hYiFGjRsHT0xNeXl4YO3YsysrKGvEs6uZO5z5mzBjIZDKLR9++fS2OscVznzFjBjp16gQPDw8EBARg8ODBSEtLszhGo9Fg3Lhx8PX1hbu7O4YOHYq8vDyLYzIyMtC/f3+4uroiICAAL730EvR6fWOeSq3V5Nx79ux5w3V/9tlnLY6xxXOfO3cuEhISzBOnJScnY+3ateb99nrN73Te9nq9b+b999+HTCbDxIkTzdvs9bpf62bnbdfXXdANfvvtN/Haa6+JZcuWCQBi+fLlFvvff/99oVKpxIoVK8Thw4fFwIEDRXR0tKisrDQf07dvX5GYmCj+/PNPsX37dtGiRQsxYsSIRj6T2rvTuY8ePVr07dtX5OTkmB+FhYUWx9jiuffp00csWLBAHDt2TKSmpop+/fqJiIgIUVZWZj7m2WefFeHh4WLjxo1i//79omvXrqJbt27m/Xq9XsTFxYmUlBRx6NAh8dtvvwk/Pz8xZcoUKU6pxmpy7j169BBPP/20xXUvKSkx77fVc1+1apVYs2aNOH36tEhLSxOvvvqqcHR0FMeOHRNC2O81v9N52+v1vt7evXtFVFSUSEhIEBMmTDBvt9frftWtztuerzvDzh1c/4FvNBpFUFCQ+OCDD8zbiouLhVKpFD/88IMQQogTJ04IAGLfvn3mY9auXStkMpnIzs5utNrv1q3CzqBBg275HHs59/z8fAFAbN26VQhhusaOjo7ip59+Mh9z8uRJAUDs3r1bCGEKinK5XOTm5pqPmTt3rvD09BRVVVWNewJ34fpzF8L0Jnjtm+L17OXchRDC29tbzJ8/v0ldcyH+Om8hmsb1Li0tFS1bthTr16+3OF97v+63Om8h7Pu6sxurltLT05Gbm4uUlBTzNpVKhS5dumD37t0AgN27d8PLywtJSUnmY1JSUiCXy7Fnz55Gr7m+bdmyBQEBAWjdujWee+45XLlyxbzPXs69pKQEAODj4wMAOHDgAHQ6ncV1j4mJQUREhMV1j4+PR2BgoPmYPn36QK1W4/jx441Y/d25/tyvWrx4Mfz8/BAXF4cpU6agoqLCvM8ezt1gMGDp0qUoLy9HcnJyk7nm15/3VfZ+vceNG4f+/ftbXF/A/v+v3+q8r7LX686FQGspNzcXACwu9tWfr+7Lzc1FQECAxX6FQgEfHx/zMbaqb9++GDJkCKKjo3Hu3Dm8+uqreOihh7B79244ODjYxbkbjUZMnDgR3bt3R1xcHADTNXVycoKXl5fFsddf95v9u7i6zxbc7NwBYOTIkYiMjERISAiOHDmCV155BWlpaVi2bBkA2z73o0ePIjk5GRqNBu7u7li+fDnatGmD1NRUu77mtzpvwL6vNwAsXboUBw8exL59+27YZ8//12933oB9X3eGHaqV4cOHm7+Pj49HQkICmjdvji1btqBXr14SVlZ/xo0bh2PHjmHHjh1Sl9LobnXuzzzzjPn7+Ph4BAcHo1evXjh37hyaN2/e2GXWq9atWyM1NRUlJSX4+eefMXr0aGzdulXqshrcrc67TZs2dn29MzMzMWHCBKxfvx7Ozs5Sl9NoanLe9nzd2Y1VS0FBQQBww8j8vLw8876goCDk5+db7Nfr9SgsLDQfYy+aNWsGPz8/nD17FoDtn/v48eOxevVqbN68GWFhYebtQUFB0Gq1KC4utjj++ut+s38XV/dZu1ud+8106dIFACyuu62eu5OTE1q0aIGOHTtixowZSExMxJw5c+z+mt/qvG/Gnq73gQMHkJ+fjw4dOkChUEChUGDr1q345JNPoFAoEBgYaJfX/U7nbTAYbniOPV13hp1aio6ORlBQEDZu3GjeplarsWfPHnN/d3JyMoqLi3HgwAHzMZs2bYLRaDT/47EXWVlZuHLlCoKDgwHY7rkLITB+/HgsX74cmzZtQnR0tMX+jh07wtHR0eK6p6WlISMjw+K6Hz161CLsrV+/Hp6enubuAWt0p3O/mdTUVACwuO62eO43YzQaUVVVZdfX/GaunvfN2NP17tWrF44ePYrU1FTzIykpCaNGjTJ/b4/X/U7n7eDgcMNz7Om6826smygtLRWHDh0Shw4dEgDERx99JA4dOiQuXrwohDDdeu7l5SVWrlwpjhw5IgYNGnTTW8/bt28v9uzZI3bs2CFatmxp9bdfC3H7cy8tLRUvvvii2L17t0hPTxcbNmwQHTp0EC1bthQajcb8GrZ47s8995xQqVRiy5YtFrddVlRUmI959tlnRUREhNi0aZPYv3+/SE5OFsnJyeb9V2/LfPDBB0VqaqpYt26d8Pf3t/rbMu907mfPnhVvv/222L9/v0hPTxcrV64UzZo1E/fdd5/5NWz13CdPniy2bt0q0tPTxZEjR8TkyZOFTCYTf/zxhxDCfq/57c7bnq/3rVx/F5K9XvfrXXve9n7dGXZuYvPmzQLADY/Ro0cLIUy3n7/xxhsiMDBQKJVK0atXL5GWlmbxGleuXBEjRowQ7u7uwtPTUzz55JOitLRUgrOpndude0VFhXjwwQeFv7+/cHR0FJGRkeLpp5+2uA1RCNs895udMwCxYMEC8zGVlZXin//8p/D29haurq7i4YcfFjk5ORavc+HCBfHQQw8JFxcX4efnJ/79738LnU7XyGdTO3c694yMDHHfffcJHx8foVQqRYsWLcRLL71kMf+GELZ57k899ZSIjIwUTk5Owt/fX/Tq1cscdISw32t+u/O25+t9K9eHHXu97te79rzt/brLhBCi8dqRiIiIiBoXx+wQERGRXWPYISIiIrvGsENERER2jWGHiIiI7BrDDhEREdk1hh0iIiKyaww7REREZNcYdoiIiMiuMewQETWyLVu2QCaTQSaTYfDgwbV67pgxY8zPXbFiRYPUR2RvGHaImrirH5y3erz55ptSl1jvoqKiMHv2bKnLQFpaGhYuXGj+ecyYMTeEn59//hnOzs74z3/+AwCYM2cOcnJyGrFKItunkLoAIpLWtR+cP/74I6ZOnYq0tDTzNnd3dynKqjUhBAwGAxSKxntb02q1cHJyqvPzAwIC4OXldcv98+fPx7hx4zBv3jw8+eSTAACVSgWVSlXnP5OoKWLLDlETFxQUZH6oVCrIZDKLbUuXLkVsbCycnZ0RExODL774wvzcCxcuQCaT4X//+x/uvfdeuLi4oFOnTjh9+jT27duHpKQkuLu746GHHsLly5fNz7vagvHWW2/B398fnp6eePbZZ6HVas3HGI1GzJgxA9HR0XBxcUFiYiJ+/vln8/6rXUFr165Fx44doVQqsWPHDpw7dw6DBg1CYGAg3N3d0alTJ2zYsMH8vJ49e+LixYuYNGmSufUKAN588020a9fO4u9m9uzZiIqKuqHud999FyEhIWjdujUAIDMzE48++ii8vLzg4+ODQYMG4cKFC3d1XWbNmoXnn38eS5cuNQcdIqobtuwQ0S0tXrwYU6dOxWeffYb27dvj0KFDePrpp+Hm5obRo0ebj5s2bRpmz56NiIgIPPXUUxg5ciQ8PDwwZ84cuLq64tFHH8XUqVMxd+5c83M2btwIZ2dnbNmyBRcuXMCTTz4JX19fvPvuuwCAGTNm4L///S/mzZuHli1bYtu2bXjsscfg7++PHj16mF9n8uTJ+PDDD9GsWTN4e3sjMzMT/fr1w7vvvgulUonvvvsOAwYMQFpaGiIiIrBs2TIkJibimWeewdNPP13rv5ONGzfC09MT69evBwDodDr06dMHycnJ2L59OxQKBaZPn46+ffviyJEjdWr5eeWVV/DFF19g9erV6NWrV62fT0TXkXjVdSKyIgsWLBAqlcr8c/PmzcWSJUssjnnnnXdEcnKyEEKI9PR0AUDMnz/fvP+HH34QAMTGjRvN22bMmCFat25t/nn06NHCx8dHlJeXm7fNnTtXuLu7C4PBIDQajXB1dRW7du2y+LPHjh0rRowYIYQQYvPmzQKAWLFixR3Pq23btuLTTz81/xwZGSk+/vhji2OmTZsmEhMTLbZ9/PHHIjIy0qLuwMBAUVVVZd72/fffi9atWwuj0WjeVlVVJVxcXMTvv/9+03qu1l5UVGSxffTo0cLJyemGv7+bASCWL19+22OIyIQtO0R0U+Xl5Th37hzGjh1r0QKi1+tvGDOSkJBg/j4wMBAAEB8fb7EtPz/f4jmJiYlwdXU1/5ycnIyysjJkZmairKwMFRUV6N27t8VztFot2rdvb7EtKSnJ4ueysjK8+eabWLNmDXJycqDX61FZWYmMjIzanP4txcfHW7TWHD58GGfPnoWHh4fFcRqNBufOnav16yckJKCgoADTpk1D586dbWbMFJE1Y9ghopsqKysDAHz99dfo0qWLxT4HBweLnx0dHc3fXx0Dc/02o9FY6z97zZo1CA0NtdinVCotfnZzc7P4+cUXX8T69evx4YcfokWLFnBxccGwYcMsxgPdjFwuhxDCYptOp7vhuOv/vLKyMnTs2BGLFy++4Vh/f//b/pk3Exoaip9//hn3338/+vbti7Vr194QpIiodhh2iOimAgMDERISgvPnz2PUqFH1/vqHDx9GZWUlXFxcAAB//vkn3N3dER4eDh8fHyiVSmRkZFiMz6mJnTt3YsyYMXj44YcBmMLI9YOFnZycYDAYLLb5+/sjNzcXQghzYEtNTb3jn9ehQwf8+OOPCAgIgKenZ61qvZXIyEhs3brVHHjWrVvHwEN0F3g3FhHd0ltvvYUZM2bgk08+wenTp3H06FEsWLAAH3300V2/tlarxdixY3HixAn89ttvmDZtGsaPHw+5XA4PDw+8+OKLmDRpEhYtWoRz587h4MGD+PTTT7Fo0aLbvm7Lli2xbNkypKam4vDhwxg5cuQNrUpRUVHYtm0bsrOzUVBQAMB0l9bly5cxa9YsnDt3Dp9//jnWrl17x/MYNWoU/Pz8MGjQIGzfvh3p6enYsmUL/vWvfyErK6vOfz/h4eHYsmUL8vPz0adPH6jV6jq/FlFTx7BDRLf0f//3f5g/fz4WLFiA+Ph49OjRAwsXLkR0dPRdv3avXr3QsmVL3Hffffj73/+OgQMHWkxg+M477+CNN97AjBkzEBsbi759+2LNmjV3/LM/+ugjeHt7o1u3bhgwYAD69OmDDh06WBzz9ttv48KFC2jevLm5qyk2NhZffPEFPv/8cyQmJmLv3r148cUX73gerq6u2LZtGyIiIjBkyBDExsZi7Nix0Gg0d93SExYWhi1btqCgoICBh+guyMT1ndRERA1szJgxKC4ubrLLHWzZsgX3338/ioqKbjup4O3IZDIsX7681stNEDVFbNkhIpJIWFgYRowYUavnPPvss7xDi6iW2LJDRI2uqbfsVFZWIjs7G4BpOY6goKAaPzc/P9/cnRUcHHzD3WFEdCOGHSIiIrJr7MYiIiIiu8awQ0RERHaNYYeIiIjsGsMOERER2TWGHSIiIrJrDDtERERk1xh2iIiIyK4x7BAREZFd+38jNydEJe6rBQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "z = [0.5, 0.5]\n", "\n", "env = model.phase_envelope_pt(z, kind=\"dew\", p0=0.01)\n", "\n", "env.plot()" ] }, { "cell_type": "code", "execution_count": null, "id": "2a996103", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "thermo", "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.13.11" } }, "nbformat": 4, "nbformat_minor": 5 }