{ "cells": [ { "cell_type": "markdown", "id": "1ec37f01", "metadata": {}, "source": [ "# VLLE @ constant pressure\n", "\n", "Following the approach described in Bell et al.: https://doi.org/10.1021/acs.iecr.1c04703, but slightly different because the pressure is fixed rather than the temperature, but the same basic principles hold\n", "\n", "for the mixture of nitrogen + ethane, with the default thermodynamic model in teqp, which is the GERG-2008 mixing parameters (no departure function).\n", "\n", "Two traces are made, and the intersection is obtained, this gives you the VLLE solution." ] }, { "cell_type": "code", "execution_count": 1, "id": "29a2031a", "metadata": { "execution": { "iopub.execute_input": "2025-10-15T23:08:17.007635Z", "iopub.status.busy": "2025-10-15T23:08:17.007489Z", "iopub.status.idle": "2025-10-15T23:08:18.549381Z", "shell.execute_reply": "2025-10-15T23:08:18.548846Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "rhovec / mol/m^3 | T / K\n", "[4921.97976373 9.6755684 ] 125.14729018874179\n", "[ 6008.68040253 15630.22353351] 125.14729018874179\n", "[18948.39537895 1540.60935171] 125.14729018874179\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjwAAAHKCAYAAAAHLJEDAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjcsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvTLEjVAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAe4dJREFUeJzt3Xd4FPXaxvHvpveEBNKA0DsISA0gICAdAfFgQUTFHixwjth7QTm+do9YARXECkIEFKkivSm9d0gogVRSd94/hixEQgkkmc3m/lzXXLtTduaZWWAfftVmGIaBiIiIiAtzszoAERERkZKmhEdERERcnhIeERERcXlKeERERMTlKeERERERl6eER0RERFyeEh4RERFxeUp4RERExOUp4RERERGXp4RHxElUr16dO+64w+owSk1aWhrh4eFMmjSpyJ/NycmhatWq/O9//yuByErehAkTsNlsrFq1yupQLqo8f0/iWpTwiDipTZs28cILL7Bnzx6rQykR7777LoGBgdx8882ObYcPH+aJJ57g2muvJTAwEJvNxoIFC875rKenJ6NGjeLVV18lMzOzFKMumv/9739MmDDB6jCuSGHf09y5c7nrrruoW7cufn5+1KxZk7vvvpvDhw8X+GxZ+Z6kfFDCI+Iktm7dyqeffupY37RpEy+++KJLJjw5OTm8++673H333bi7uzu2b926lTfeeIODBw/SpEmTC57jzjvv5NixY0yePLmkw71sZT3hOd/39Pjjj7NgwQIGDhzIe++9x80338x3331H8+bNSUhIKHCOsvA9SfmghEfESXh7e+Pp6XlZnzUMg1OnThVzRCUnPj6eo0ePMnjw4ALbW7RowfHjx9m2bRujRo264DlCQkLo3r37ZScUL7zwAtWrV7+sz5YX5/ue3nrrLXbs2MEbb7zB3XffzWuvvUZ8fDyJiYl88MEHBY690u9JpLgo4REpQS+88AI2m40dO3Zwxx13EBISQnBwMHfeeScZGRkFjj27Dc+ECRP417/+BcC1116LzWYrUL1TvXp1+vbty6+//krLli3x9fXl448/BmDXrl3861//IjQ0FD8/P9q2bcsvv/xyTmx79+7l+uuvx9/fn/DwcEaOHMmvv/5aaDXS8uXL6dmzJ8HBwfj5+dGpUyf+/PPPy77XadOmUb16dWrVqlVge2BgIKGhoZf8fK+77joWL15MUlLSJX+mONjtdt555x0aNWqEj48PERER3HfffZw4ccJxTPXq1dm4cSMLFy50fH+dO3cucJ6srCxGjRpFpUqV8Pf3Z+DAgRw9erTAMT///DN9+vQhOjoab29vatWqxcsvv0xeXl6B4zp37kzjxo3ZtGkT1157LX5+flSuXJmxY8eeE39WVhbPP/88tWvXxtvbm6pVqzJ69GiysrIKHHe+76ljx464ubmdsy00NJTNmzefcz2rvieRs3lYHYBIeTB48GBq1KjBmDFjWLNmDZ999hnh4eG88cYbhR7fsWNHHn74Yd577z2eeuopGjRoAOB4BbP655ZbbuG+++7jnnvuoV69eiQmJtKuXTsyMjJ4+OGHCQsLY+LEiVx//fX88MMPDBw4EID09HS6dOnC4cOHeeSRR4iMjGTy5MnMnz//nFjmzZtHr169aNGiBc8//zxubm6MHz+eLl268Mcff9C6desi3+uSJUu4+uqrr/i5tmjRAsMwWLJkCX379r3i812q++67jwkTJnDnnXfy8MMPs3v3bj744APWrl3Ln3/+iaenJ++88w4PPfQQAQEBPP300wBEREQUOM9DDz1EhQoVeP7559mzZw/vvPMOI0aM4Ntvv3UcM2HCBAICAhg1ahQBAQHMmzeP5557jpSUFP773/8WON+JEyfo2bMnN9xwA4MHD+aHH37g8ccfp0mTJvTq1Qswk7Xrr7+exYsXc++999KgQQPWr1/P22+/zbZt25g2bZrjfEX5ntLS0khLS6NixYrn7LPqexIpwBCREvP8888bgHHXXXcV2D5w4EAjLCyswLZq1aoZw4YNc6x///33BmDMnz//nPNWq1bNAIzZs2cX2P7oo48agPHHH384tqWmpho1atQwqlevbuTl5RmGYRj/93//ZwDGtGnTHMedOnXKqF+/foFr2u12o06dOkaPHj0Mu93uODYjI8OoUaOGcd111xX5XnNycgybzWb8+9//LuyRXdL95zt06JABGG+88cYFz1WY559/3qhWrVqRP/fHH38YgDFp0qQC22fPnn3O9kaNGhmdOnU65xzjx483AKNbt24FnuvIkSMNd3d34+TJk45tGRkZ53z+vvvuM/z8/IzMzEzHtk6dOhmA8eWXXzq2ZWVlGZGRkcagQYMc27766ivDzc2twJ8RwzCMcePGGYDx559/GoZx6d9TvpdfftkAjLlz556z70q+J5HioiotkVJw//33F1i/5pprOH78OCkpKZd9zho1atCjR48C22bOnEnr1q3p0KGDY1tAQAD33nsve/bsYdOmTQDMnj2bypUrc/311zuO8/Hx4Z577ilwvnXr1rF9+3ZuvfVWjh8/zrFjxzh27Bjp6el07dqVRYsWYbfbi3SvSUlJGIZBhQoVLvve8+Wf49ixYxc9Nj/2/CUjIwO73X7O9n9W6/zT999/T3BwMNddd12Bz7Vo0YKAgIBCS8nO595778VmsznWr7nmGvLy8ti7d69jm6+vr+N9amoqx44d45prriEjI4MtW7YUOF9AQAC33XabY93Ly4vWrVuza9euAvE3aNCA+vXrF4i/S5cuAI74i/I9LVq0iBdffJHBgwc7znO2onxPIiVFVVoipSAmJqbAev4PwIkTJwgKCrqsc9aoUeOcbXv37qVNmzbnbM+vCtu7dy+NGzdm79691KpVq8CPLUDt2rULrG/fvh2AYcOGnTeO5OTkAj+Kl3qvhmGc95yXKv8c/7yPwlSqVOmSto8fP/6C4yFt376d5ORkwsPDC91/5MiRi8aS70LPKt/GjRt55plnmDdv3jkJcnJycoH1KlWqnPMsKlSowN9//10g/s2bN5/3efwz/ot9T1u2bGHgwIE0btyYzz77rNBjivI9iZQUJTwipeDsLr1nu5If/bP/519S8ktv/vvf/9KsWbNCjwkICCiwfrF7DQ0NxWazFfhRv1z55yis3cg/zZkzp8D6l19+yW+//cbXX39dYHujRo0ueB673X7BgfjOl0gU5mLP6uTJk3Tq1ImgoCBeeuklatWqhY+PD2vWrOHxxx8/p3TtUv6c2e12mjRpwltvvVXosVWrVgUu7Xvav38/3bt3Jzg4mJkzZxIYGFjocUX5nkRKihIeESd1Of8brlatGlu3bj1ne37VR7Vq1RyvmzZtwjCMAtfZsWNHgc/l984JCgqiW7duRY6nMB4eHtSqVYvdu3df8bnyz3F2Y+7z+Wf8ixcvxsfHp8j3VatWLX7//Xfat29/0aTzSks0FixYwPHjx/npp5/o2LGjY/uVPLtatWrx119/0bVr1wvGd7Hv6fjx43Tv3p2srCzmzp1LVFTUec9VlO9JpKSoDY+Ik/L39wfM/+Vfqt69e7NixQqWLl3q2Jaens4nn3xC9erVadiwIQA9evTg4MGDTJ8+3XFcZmZmgYEPwexdU6tWLd58803S0tLOud4/u1BfqtjY2GKZVmH16tXYbDZiY2Ov+FyXavDgweTl5fHyyy+fsy83N7fA9+Xv71+k7++f8ktszi6hyc7OvqKpGgYPHszBgwfP+a4BTp06RXp6umP9fN9Teno6vXv35uDBg8ycOZM6depc8JpWfE8i/6QSHhEn1axZM9zd3XnjjTdITk7G29ubLl26nLftCMATTzzBN998Q69evXj44YcJDQ1l4sSJ7N69mx9//NExdsp9993HBx98wC233MIjjzxCVFQUkyZNwsfHBzhTMuHm5sZnn31Gr169aNSoEXfeeSeVK1fm4MGDzJ8/n6CgIGbMmFHke+vfvz9fffUV27Zto27dugX2vfLKK4DZdgXgq6++YvHixQA888wzBY6dM2cO7du3JywsrMgxXK5OnTpx3333MWbMGNatW0f37t3x9PRk+/btfP/997z77rvceOONgJkwfvTRR7zyyivUrl2b8PDwQhv1nk+7du2oUKECw4YN4+GHH8Zms/HVV19dUVXo0KFD+e6777j//vuZP38+7du3Jy8vjy1btvDdd985xnaC839PQ4YMYcWKFdx1111s3ry5wNg7AQEBDBgwoMA1rfieRM5hSd8wkXIiv6v20aNHC2zP75a8e/dux7Z/dks3DMP49NNPjZo1axru7u4FumhXq1bN6NOnT6HX3Llzp3HjjTcaISEhho+Pj9G6dWsjPj7+nON27dpl9OnTx/D19TUqVapk/Pvf/zZ+/PFHAzCWLVtW4Ni1a9caN9xwgxEWFmZ4e3sb1apVMwYPHlygC3JR7jUrK8uoWLGi8fLLL58TF3De5WwnT540vLy8jM8++6zQ53Axl9stPd8nn3xitGjRwvD19TUCAwONJk2aGKNHjzYOHTrkOCYhIcHo06ePERgYaACOLur5z2TlypUFzjl//vxzuuL/+eefRtu2bQ1fX18jOjraGD16tPHrr7+ec1ynTp2MRo0anRPnsGHDzrnP7Oxs44033jAaNWpkeHt7GxUqVDBatGhhvPjii0ZycrLjuPN9T/nDIhS2/PNaV/o9iRQXm2EUQ1cJEXEJ77zzDiNHjuTAgQNUrly5RK/18ssvM378eLZv337exrYX8s477zB27Fh27txZKg24yyt9T+Iq1IZHpJz659xbmZmZfPzxx9SpU6fEkx2AkSNHkpaWxpQpU4r82ZycHN566y2eeeYZ/YiWMH1P4ipUwiNSTvXq1YuYmBiaNWtGcnIyX3/9NRs3bmTSpEnceuutVocnIlKs1GhZpJzq0aMHn332GZMmTSIvL4+GDRsyZcoUbrrpJqtDExEpdpZXaY0ZM4ZWrVoRGBhIeHg4AwYMOGcckYSEBIYOHUpkZCT+/v5cffXV/PjjjwWOSUpKYsiQIQQFBRESEsLw4cML7UYrIqZHH32UDRs2kJaWxqlTp1i9erWSHRFxWZYnPAsXLiQuLo5ly5YxZ84ccnJy6N69e4GxIG6//Xa2bt3K9OnTWb9+vWMm4LVr1zqOGTJkCBs3bmTOnDnEx8ezaNEi7r33XituSURERJyM07XhOXr0KOHh4SxcuNAxsmhAQAAfffQRQ4cOdRwXFhbGG2+8wd13383mzZtp2LAhK1eudIwfMXv2bHr37s2BAweIjo625F5ERETEOThdG578yfBCQ0Md29q1a8e3335Lnz59CAkJ4bvvviMzM5POnTsDsHTpUkJCQhzJDpjDyLu5ubF8+XIGDhx4znWysrIKzIpst9tJSkoiLCxME9yJiIiUEYZhkJqaSnR0tGNw1cI4VcJjt9t59NFHad++PY0bN3Zs/+6777jpppsICwvDw8MDPz8/pk6d6pjZOSEh4ZzRZz08PAgNDSUhIaHQa40ZM4YXX3yx5G5GRERESs3+/fupUqXKefc7VcITFxfHhg0bHMPI53v22Wc5efIkv//+OxUrVmTatGkMHjyYP/74gyZNmlzWtZ588klGjRrlWE9OTiYmJob9+/cTFBR0RffxT3a7Qc93F3HoZCavDmhM/+aFjHFiGJC4EbbNhq2z4NjZDbdtULkl1OsBdXpCaPVijU9ERKSsSklJoWrVqgQGBl7wOKdpwzNixAh+/vlnFi1aRI0aNRzbd+7cSe3atdmwYQONGjVybO/WrRu1a9dm3LhxfPHFF/z73//mxIkTjv25ubn4+Pjw/fffF1ql9U8pKSkEBweTnJxc7AkPwIfzd/DfX7fSsloFfnig3cU/cHwnbJ4Bm6fDwdUF90U0gQZ9oX5fiGgEqoITEZFy6lJ/vy3vpWUYBiNGjGDq1KnMmzevQLIDkJGRAXBOvZy7uzt2ux0wZ/Q9efIkq1efSQzmzZuH3W6nTZs2JXwHl+ZfLarg7mZj1d4T7DiSevEPhNWCDo/CPfNg5Cbo9V+o0RFs7pC4HhaMgXHt4b3m8NszsG85nH4eIiIiUpDlJTwPPvggkydP5ueff6ZevXqO7cHBwfj6+pKTk0PDhg2JiorizTffJCwsjGnTpvHYY48RHx9P7969AXPU2MTERMaNG0dOTg533nknLVu2ZPLkyZcUR0mX8ADcPXEVv29O5O4ONXimb8PLO0lGklnltSUedsyFvDMNrwmIgHq9zdKf6h3Bw6t4AhcREXFSl/r7bXnCc74eUePHj+eOO+4AYPv27TzxxBMsXryYtLQ0ateuzX/+858C3dSTkpIYMWIEM2bMwM3NjUGDBvHee+8REBBwSXGURsLz+6ZE7v5yFaH+Xix7siteHldYwJaVBjt+N6u+tv8GWSln9nkHQ93uZrVX7W7gfWnPQUREpCwpMwmPsyiNhCc3z0671+dxJDWLcbe1oGfjyGI8eTbsXgRbZsCWmZB+5Mw+d2+o1cUs+anbC/zDiu+6IiIiFlLCU0SlkfAAjJm1mY8X7qJbg3A+G9aqZC5iz4MDK82Sny3xcGLPmX02N4hpd7rRcx8IiSmZGEREREqBEp4iKq2EZ8eRNLq9tRB3NxvLn+pKxQDvErsWYHZ3P7IJNsebpT8J6wvuj7wKGvQzk5/whurxJSIiZYoSniIqrYQHoP+Hf/LX/pM827chwzvUuPgHitOJPWaV15Z42LcUjLN6dlWoYSY+DfpBlVbg5l66sYmIiBSREp4iKs2E56ule3j25400jApi5iPXlOi1Lij92OkeX7/AznkFe3z5V4J6vaB+P7M7vKePdXGKiIichxKeIirNhOdEejatX/udnDyDXx/tSL3IC48OWSqy0mDnXLPqa9uvkJV8Zp9XANS5zuzxVec68Am2Lk4REZGzKOEpotJMeADu+XIVczYlcn+nWjzRq36JX69I8nJgz2Kz2mvLL5B6+Mw+N0+ocY2Z/NTrDUFR1sUpIiLlnhKeIirthGfW+sM8MGkNUcE+/Pl4F9zcnLSxsN0Oh9eebvT8yz/m+MKc46t+HzMBqlTXmhhFRKTcUsJTRKWd8GTm5NHqld9Jzcrlu/tiaV0jtMSvWSyObTe7u2+daXZ9P1tYnTPJT+UW4Gb5zCUiIuLilPAUUWknPAD/+f4vflh9gCFtYnh14OXN+m6p1AQz8dnyC+xaCPacM/vyp7mo39esAvMo4e73IiJSLinhKSIrEp4/th9l6OcrqODnyYqnu+HpXoZLRDJTYMccM/nZPqfgNBdegVCnmxo9i4hIsbvU32+PUoxJ/iG2ZhgVA7w4lpbNkp3H6VS3ktUhXT6fIGg8yFxys2HPH2bys+UXSEuAjVPNxdHouc/pRs/RVkcuIiLlgEp4TrOihAfgmWnr+XrZPga3rMLYG5uW2nVLjd0Oh9ae6fH1z0bP0Vef1ei5nkZ6FhGRIlGVVhFZlfAs23Wcmz9ZRpCPB6ueue7KZ1B3dsd2wNbTJT/7VwBn/fELrQX1T7f70UjPIiJyCZTwFJFVCU+e3aDtmLkcTc3iizta0qV+RKld23KpibDt9EjPuxZAXvaZfY6RnvtCjU4a6VlERAqlhKeIrEp4AJ7/eQMTl+7lXy2q8N9/uWC11qXISoUdc83k558jPXv6Q+2uZvJTtzv4VrAuThERcSpqtFyG9GwcxcSle5mzOZGcPHvZ7q11ubwDodEAc8kf6Tm/y3vKQdg83Vxs7lC9/ZmRnkOqWh25iIiUASrhOc3KEp48u0HrV3/neHo2k+5uQ/vaFUv1+k7NMODwujM9vo5sKrg/qinU62M2fI5opEbPIiLljKq0isjKhAfgiR//ZsrK/QxtW42XBzQu9euXGUm7YMvpkp/9y8Cwn9kXUs0s+anfB6q2AXcVYIqIuDolPEVkdcIzf+sR7hy/kvBAb5Y92dV559ZyJunHYNtsM/nZOQ9yM8/s8w093ei5D9S8Frz8rItTRERKjNrwlDHtaoUR6O3BkdQs1uw7QcvqZWRuLSv5V4Tmt5lLdrqZ9GyZafb8OpUE6yaZi4fv6UbPfaBuT/DTsxURKW+U8DgJbw93ujQI5+d1h5izOVEJT1F5+UODfuaSlwv7lp5p95O87/TAh/Fgc4Nq7c+M9FyhmtWRi4hIKVCV1mlWV2kBzPjrEA99s5ZalfyZ++/OlsTgcgwDEjfA5tMjPSeuL7g/ssmZdj8RjdXoWUSkjFEbniJyhoQnJTOHq1+aQ67dYP5/OlOjor8lcbi0E3vONHret+QfjZ5jzmr03FaNnkVEygAlPEXkDAkPwJDPlvHnjuM806cBd19T07I4yoX042c1ep6rRs8iImWQEp4icpaEZ/yfu3lxxiba1Ajl2/tiLYuj3MlOh53zT4/0PAtOnTizT42eRUSclhKeInKWhGd/UgbXjJ2Pu5uN1c90I8TPy7JYyq3CGj3ns7lDtXZnqr400rOIiKWU8BSRsyQ8AD3eXsTWxFTevbkZ/ZtVtjSWcs8wIGH96eQn3mwAfbaoplC/HzToC5Xqq9GziEgp0zg8ZVjn+pXYmpjKwq1HlfBYzWaDqKvM5donIWm3OcfX5nizFOjwX+Yy/xUIrWUmPvX7QeUW4FYO50QTEXFSKuE5zZlKeJbuPM4tny4jzN+LlU9306jLzirt6JkJTnfNh7zsM/sCIqF+b7Pqq/o14KGqSRGRkqAqrSJypoQnJ89O85fmkJaVy89x7WlaNcTSeOQSZKXC9jlmtde23yA79cw+72Co28Ms/andzRwkUUREisWl/n5bXuY+ZswYWrVqRWBgIOHh4QwYMICtW7eec9zSpUvp0qUL/v7+BAUF0bFjR06dOuXYn5SUxJAhQwgKCiIkJIThw4eTlpZWmrdSbDzd3WhfOwyABVuPWhyNXBLvQGh8A9z4BYzeCUN+gKuHgX8lyEqG9d/Bd7fD2Jow+WZY+7XZLV5EREqF5QnPwoULiYuLY9myZcyZM4ecnBy6d+9Oenq645ilS5fSs2dPunfvzooVK1i5ciUjRozA7aw2EkOGDGHjxo3MmTOH+Ph4Fi1axL333mvFLRWLzvXCAViw7YjFkUiReXhDnevg+vfg31vhrl8hdgRUqG6O9bNtFvwcB2/Whgl9Ydk4OLnf6qhFRFya01VpHT16lPDwcBYuXEjHjh0BaNu2Lddddx0vv/xyoZ/ZvHkzDRs2ZOXKlbRs2RKA2bNn07t3bw4cOEB0dPRFr+tMVVoAh5NPETtmHjYbrHnmOir4qw1ImWcYkLjxdI+vGWbvr7NFNz89H1h/qFjbmhhFRMqYMlOl9U/JyckAhIaag7sdOXKE5cuXEx4eTrt27YiIiKBTp04sXrzY8ZmlS5cSEhLiSHYAunXrhpubG8uXLy/0OllZWaSkpBRYnElUsC91IwIwDFi6S1UfLsFmg8jG0PlxuH8xPPIX9HgNYtoBNji0Fua+BB+0gA/bwrxX4fDfZqIkIiJXxKkSHrvdzqOPPkr79u1p3LgxALt27QLghRde4J577mH27NlcffXVdO3ale3btwOQkJBAeHh4gXN5eHgQGhpKQkJCodcaM2YMwcHBjqVqVecbQK597YoA/LH9mMWRSImoUB1i4+CuWfCfbdD3HajVFdw84OhmWDQWPr4G3msGvz0D+1eA3X6Rk4qISGGcKuGJi4tjw4YNTJkyxbHNfvof+Pvuu48777yT5s2b8/bbb1OvXj2++OKLy77Wk08+SXJysmPZv9/52lBcU8dMeBbvUMNllxcQDi3vhKE/wWM7YOAnZpd2D19zwtMl78Pn18FbDeCXf8OuBeaI0CIickmcZuDBESNGOBobV6lSxbE9KioKgIYNGxY4vkGDBuzbZw75HxkZyZEjBRv35ubmkpSURGRkZKHX8/b2xtvbuzhvodi1qRGGh5uN/Umn2Hc8g5gwTWBZLvhWgKY3mUt2Ouz4HTbPgG2/QloCrPzMXHwrQL0+Zrufmp3B08fqyEVEnJblJTyGYTBixAimTp3KvHnzqFGjRoH91atXJzo6+pyu6tu2baNatWoAxMbGcvLkSVavXu3YP2/ePOx2O23atCn5mygh/t4eXB1TAYA/VMpTPnn5Q8P+MOgzs+RnyA/QfCj4hZkTnK77Gr65Cf5bC76/Ezb8ZI4JJCIiBVhewhMXF8fkyZP5+eefCQwMdLS5CQ4OxtfXF5vNxmOPPcbzzz9P06ZNadasGRMnTmTLli388MMPgFna07NnT+655x7GjRtHTk4OI0aM4Oabb76kHlrOrEOdiqzYk8Ti7ccY0qaa1eGIlfK7u9e57swEp5tnmEvqIdj4k7m4e0OtLtDweqjXG3xDrI5cRMRylndLt51nssXx48dzxx13ONZff/11PvzwQ5KSkmjatCljx46lQ4cOjv1JSUmMGDGCGTNm4ObmxqBBg3jvvfcICAi4pDicrVt6vjX7TnDD/5YQ7OvJ2mev0zQTci67HQ6tgc3TYdN0OLH7zD43T7O6q2F/c3Z3v1DLwhQRKQmaWqKInDXhycmz0+zF30jPzmPmw9fQMNp5YhMnZBhwZJOZ+Gz62eztlc/NA2p0PJ389AX/itbFKSJSTMrsODxSkKe7Gy2rm/8rX6bxeORibDaIaGTO7B63DOJWwLVPQ0RjsOfCznkw4xF4sy5MvB5Wfg5pGs1bRFyfEp4yILaWOa+WBiCUIqtUDzqNhgf+hBGroetzEHkVGHmweyH8Mgr+r545xcWKTyG18HGrRETKOlVpneasVVoA6/afZMCHf6odjxSfpF1nqr0OrTlrhw1i2kLDAWZ39+DKVkUoInJJ1IaniJw54cnNs9P0dDueXx7uQKPoYKtDEldyYu/pBs8/w4GVBfdVaQ2NBkCD6yHE+UYjFxFRGx4X4uHuRqsa+e14kiyORlxOhWrQ7iG4+3cYuRF6jIGqbc19B1bAr0/BO43h0y7w57vmyM8iImWMEp4yom3N0+14dqodj5Sg4CoQ+yAM/xVGbYZeY6Fae8AGB1fDnOfg3abwybXmdBfJB6yOWETkklg+8KBcmjanS3hW703CbjfUjkdKXlA0tLnPXFITYcsMs9prz2Kz3c+hNeakplXbQuNBZnf3wAiroxYRKZTa8JzmzG14ALJz7Vz14q9k5tj5fVRHaocHWh2SlFdpR8zEZ8NPsG/Jme02N6jeARrdYCY/GuRQREqB2vC4GC8PN5pWCQFg5Z4T1gYj5VtAOLS+B+6aBSM3QY/XoHJLMOywexHEPwpv1oGvB8G6yZCZbHXEIiJKeMqSVqcHIFylhEecRXBliI2De+bCI39Btxcgsok5yOGO32HaA/Df2vDNrbD+B8hKszpiESmn1IanDGlR3Zw5fdVe9dQSJ1ShOnQYaS7HtptVXht+hGNbYesv5uLhC3V7mG1+6lwHnr5WRy0i5YQSnjLk6pgK2Gyw93gGR1IzCQ/0sTokkcJVrAOdHzdHeT6y6Uzyc2I3bJpmLl4B5oSmjW4wZ3f38LI6ahFxYUp4ypBgX0/qRQSyJSGV1XtO0KtJlNUhiVxY/txeEY2gyzNweJ2Z+GycBsn74e9vzcUnBBr0NUt+qncEd/3TJCLFS214ypiWp6u11HBZyhybDaKbQ/dX4JG/4a7foM39EBABmSdh7dfw1UBzbq/4UbBvuTn7u4hIMVDCU8a0rHZ6PJ59SnikDHNzg5g20OsNc4DDYfHQ4k7wC4OMY7Dqc/iiO7zXDOa/Bsd3Wh2xiJRxGofnNGcfhyff3uPpdPrvArzc3Vj/Yne8PdytDkmk+OTlmLO4r//BnNw0J/3Mvsot4aqboPEN4F/RuhhFxKloHB4XFRPqR6i/F9l5djYdSrE6HJHi5e4JtbvBwHHw2Ha44TOofR3Y3OHgKpj1GLxZFyYNNpOi7AyrIxaRMkItA8sYm81G0yrBzN96lHX7T9I8poLVIYmUDC9/uOpf5pJ2xGzs/Pe3cGgtbP/VXLwCoeH1cNVgqH4NuKnEU0QKpxKeMqhZVTPJWbf/pLWBiJSWgHBo+wDcuwDiVkLHxyAkBrJTYd0k+LI/vN0YfnsWEjZYHa2IOCElPGVQs5gQQAmPlFOV6ppd3B/+C+6cbTZ29gmB1EOw5D0Y1x7+1w4WvwPJB62OVkSchBotn1ZWGi0DJGfk0PSl3wBY/Uw3wgK8LY5IxGK5WbD9N7PKa9uvkJd9eofNnNC06c3Q4Hrwce6/2yJSdGq07MKC/TypWdEfgL8OnLQ2GBFn4OENDfrBTV/Df7ZBv3chph1gwJ4/4Oc4c0LT7++ArbPM3mAiUq4o4SmjmlUNAWDdvpOWxiHidHwrQIs7zNncH/kbujwLFetCbiZsnArf3Gz29PrlP7B/pQY3FCknlPCUUfnteNaqHY/I+VWoBh3/A3ErzAbPbR8E/3A4lQQrP4XPu8H7V8Mf/wcph62OVkRKkBKeMiq/hGf9wWTUDEvkIvKnteg5xhzZ+bYfoclg8PSDpF0w9yV4uxFMvhm2zIS8XKsjFpFipnF4yqi6EYF4uNk4mZHDwZOnqFLBz+qQRMoGdw9zcMPa3SArDTb9DGu+hP3LYNsscwmIhGa3QvPbIKyW1RGLSDFQCU8Z5ePpTt2IQAA2HNSIyyKXxTsAmg+B4b+a1V6xI8z5vNISYPFbZnXXhL7w93eQc8rqaEXkCijhKcMaVza73208lGxxJCIuoFI96PEqjNoCg780S4Cwmb28frrHnMX9l//A4b+tjlRELoMSnjKsceVgADYcVMIjUmw8vKBhf7Odz6ProfNTEFwVMpPNhs4fXwMfd4KVn5nbRKRMUMJThjWKPp3waBJRkZIRUhU6Pw6P/AW3/QQNB4CbJxxeB7/8G96sB1Pvh71L1L1dxMlZnvCMGTOGVq1aERgYSHh4OAMGDGDr1q2FHmsYBr169cJmszFt2rQC+/bt20efPn3w8/MjPDycxx57jNxc1+5p0TAqCDcbHE3N4khKptXhiLguN3eo3RUGT4R/b4Eer0Gl+pB7Cv76Bsb3gg9amtNZpB2xOloRKYTlCc/ChQuJi4tj2bJlzJkzh5ycHLp37056evo5x77zzjvYbLZztufl5dGnTx+ys7NZsmQJEydOZMKECTz33HOlcQuW8fVyp3Z4AAAb1I5HpHT4V4TYOHhwGQz/HZoPBU9/OL4Dfn8e3moAU4acnuLCtf/TJVKWON1cWkePHiU8PJyFCxfSsWNHx/Z169bRt29fVq1aRVRUFFOnTmXAgAEAzJo1i759+3Lo0CEiIiIAGDduHI8//jhHjx7Fy8vrotctS3NpnW3Ut+v4ae1BRl1Xl4e71rE6HJHyKSvVHMV5zZdwYOWZ7YHRZi+w5rdBheqWhSfiysrsXFrJyWZJRWhoqGNbRkYGt956Kx9++CGRkZHnfGbp0qU0adLEkewA9OjRg5SUFDZu3FjodbKyskhJSSmwlEWN1HBZxHregXD17XD37/DAUnNEZ99Qcwb3Rf+Fd5vCxOth/Q+Qo+pnESs4VcJjt9t59NFHad++PY0bN3ZsHzlyJO3ataN///6Ffi4hIaFAsgM41hMSEgr9zJgxYwgODnYsVatWLaa7KF2No81sVgmPiJOIaGiO6PzvLXDjeKh5LWCD3Qvhx+HwVn2Y9QQkFv6fMREpGU410nJcXBwbNmxg8eLFjm3Tp09n3rx5rF27tliv9eSTTzJq1CjHekpKSplMehqeTngOJWdyPC2LsABviyMSEcCcwb3xDeZyYi+smwRrv4aUg7D8I3Op3AJaDofGg8DTx+qIRVya05TwjBgxgvj4eObPn0+VKlUc2+fNm8fOnTsJCQnBw8MDDw8zRxs0aBCdO3cGIDIyksTExALny18vrAoMwNvbm6CgoAJLWRTo40mNiv4AbDpcNqvlRFxehWpw7VPmuD5DfoQG14ObBxxcDT8/CG83NOfzSj5odaQiLsvyhMcwDEaMGMHUqVOZN28eNWrUKLD/iSee4O+//2bdunWOBeDtt99m/PjxAMTGxrJ+/XqOHDnTHXTOnDkEBQXRsGHDUrsXqzSIMqeY2JqQanEkInJBbu5Qpxvc9JU5onO3F8xBDTOOmzO2v9MEvhumcX1ESoDlVVpxcXFMnjyZn3/+mcDAQEebm+DgYHx9fYmMjCy0lCYmJsaRHHXv3p2GDRsydOhQxo4dS0JCAs888wxxcXF4e7t+FU+9iCBmrk9g82ElPCJlRkAl6DASYh8yJyxd/rE5jcWmaeYS0QTa3AdNbgRPX6ujFSnzLC/h+eijj0hOTqZz585ERUU5lm+//faSz+Hu7k58fDzu7u7ExsZy2223cfvtt/PSSy+VYOTOo/7pEp4tCarSEilz3D2gQT+4Ix7u/xOuHgYevpC4HqaPgLcawu8vwMn9VkcqUqY53Tg8Vimr4/AA7DueQcf/zsfL3Y1NL/XAw93yPFZErkRGEqz9ClZ8Bsn7zG02N6jfB9rcD9XaQyGDsIqUR2V2HB4puioVfPH3cic7z87uY+eOUC0iZYxfKLR/BB5ZBzdNghodwbDD5hkwoQ981B5WT4TsDKsjFSkzlPC4ADc3G/UizWqtzWq4LOI63NyhQV8YNsMc0LDFneDpB0c2woyHzWksfnvW7PYuIhekhMdF1I8yi/G2qGu6iGuKaAj93oFRm6D7KxBSDTJPwpL34L1m5vxduxepd5fIeSjhcRENIvMbLquER8Sl+VaAdg/Bw2vh5m+gZmezumtLPEzsBx+1g1VfQLaqt0XOpoTHRdSLVAmPSLni5g71e8PtP8ODy80Rmz394cgmiB9pVnf9+rR6d4mcpoTHReS34TmUnElyRo7F0YhIqQqvD33fMqu7erxmzsyemQxLPzCru6beD0c2Wx2liKWU8LiIYF9PooLNuXh2HFW1lki55BsCsXHw0Fq49Tuzd5c9F/76Bv7XFr65BfavsDpKEUso4XEhtcMDANiWmGZxJCJiKTc3qNvD7N11zzxzYENssHUmfH4djO8D239XA2cpV5TwuJC6EWa11rZElfCIyGmVW8BNX0PcCmh+G7h5wt7FMGkQjLsG1v8AeblWRylS4pTwuJC6EWYJz44jKuERkX+oVBf6fwiP/AWxI8wGzonr4cfh8EELs2dXTqbVUYqUGCU8LqSOSnhE5GKCK0OPV2HkBuj8FPiGwok9Zs+ud5rA4rfNBs8iLkYJjwvJb8OTmJJF8in11BKRC/ALhc6Pm4lPzzcgqAqkHzEnKn27sfmammh1lCLFRgmPCwnyOaun1hGV8ojIJfDyh7b3m/N2DfgIKtaDrBSzpOedJmbJT9Iuq6MUuWJKeFzMmWotteMRkSJw94Rmt8KDy+DmyVC5JeRlmW173m8BPwyHhPVWRyly2ZTwuJi6jq7pKuERkcvg5gb1+8Ddv8OweKjV1Zy6YsMPMK4DfH2jxvKRMkkJj4upo55aIlIcbDaocQ0M/QnuWwSNbgCbG+yYY47lM2kwHP7b6ihFLpkSHheT33B5pxIeESkuUU3hX+NhxCpzLB+bO2z/FT6+Br4bBke3WR2hyEUp4XExtSqZCc+h5EzSszSYmIgUo7Ba5lg+cSug8SBz26Zp8L82MPUBs3u7iJNSwuNiQvy8qBjgBcDuY+kWRyMiLqlibbjxC7j/T6jX22zj89dkeL8lxI+ClMNWRyhyDiU8Lqjm6VKenUdVrSUiJSiyMdzyDdw9D2peC/YcWPW5OUP7r09D+nGrIxRxUMLjgvKrtdSOR0RKRZUWcPs0s1dX1TaQmwlLP4B3r4J5r2rkZnEKSnhcUK1K/gDsPKoqLREpRTWugbt+hVu/h8irIDsNFo2Fd66CP96CbP2bJNZRwuOCaoWra7qIWMRmg7rd4d6F8K+J5sjNmSdh7ovwbjNYNg5ys6yOUsohJTwuqPbpKq3dx9LJsxsWRyMi5ZKbGzQaAA8uhYEfQ4Xq5lxdsx+H966G1RMhTz1JpfQo4XFB0SG+eHu4kZ1n58CJDKvDEZHyzM0dmt5sjuHT920IjIaUAzDjYfiwFfz9PdjtVkcp5YASHhfk7majRsX8djyq1hIRJ+DuCS3vgofXQI/XwK+iOSnpT3fDZ11h33KrIxQXp4THReX31Nqlhssi4kw8fSE2Dh75C7o8A16BcGgNfNEdfrwbkg9YHaG4KCU8Liq/hGeXBh8UEWfkHQAdHzNLfJoPBWyw/ntz8MIFr0O2quOleCnhcVE1T3dN360SHhFxZgHh0P8DuHcBxMRC7ilYMAY+aAXrfwBDHS+keCjhcVH5JTyaXkJEyoToZnDnLLhxPARXNRs2/zgcvugJh9ZaHZ24AMsTnjFjxtCqVSsCAwMJDw9nwIABbN261bE/KSmJhx56iHr16uHr60tMTAwPP/wwyckFR+7ct28fffr0wc/Pj/DwcB577DFyc8tvl8f8hCchRZOIikgZYbNB4xtgxEq49mnw9IP9y+CTa2FaHKQmWh2hlGGWJzwLFy4kLi6OZcuWMWfOHHJycujevTvp6WbJxKFDhzh06BBvvvkmGzZsYMKECcyePZvhw4c7zpGXl0efPn3Izs5myZIlTJw4kQkTJvDcc89ZdVuWC/HzItTfnER0z3GV8ohIGeLpC51Gm13Zr7oJMGDd1/D+1bD4bQ1cKJfFZhjOVUF69OhRwsPDWbhwIR07diz0mO+//57bbruN9PR0PDw8mDVrFn379uXQoUNEREQAMG7cOB5//HGOHj2Kl5fXRa+bkpJCcHAwycnJBAUFFes9WeXGj5awau8J3r+lOf2aRlsdjojI5dm/0hyw8OBqc71Cdej+CtTva5YKSbl2qb/flpfw/FN+VVVoaOgFjwkKCsLDwwOApUuX0qRJE0eyA9CjRw9SUlLYuHFjoefIysoiJSWlwOJq1I5HRFxC1VYw/HdzxObAKDixB769Db7sD8d3Wh2dlBFOlfDY7XYeffRR2rdvT+PGjQs95tixY7z88svce++9jm0JCQkFkh3AsZ6QkFDoecaMGUNwcLBjqVq1ajHdhfOoUUkJj4i4CDe3MyM2X/MfcPeG3Qvho3bwx/9BXo7VEYqTc6qEJy4ujg0bNjBlypRC96ekpNCnTx8aNmzICy+8cEXXevLJJ0lOTnYs+/fvv6LzOaPqYUp4RMTFeAdA12chbjnUvBZyM2HuS/BJZziw2uroxIk5TcIzYsQI4uPjmT9/PlWqVDlnf2pqKj179iQwMJCpU6fi6enp2BcZGUliYsHW+/nrkZGRhV7P29uboKCgAouryU949qrRsoi4mtAaMHSqWc3lGwqJG8wpKmY9AVmaUkfOZXnCYxgGI0aMYOrUqcybN48aNWqcc0xKSgrdu3fHy8uL6dOn4+PjU2B/bGws69ev58iRI45tc+bMISgoiIYNG5b4PTir6hX9ADiRkUNyhop7RcTF2Gynq7lWwlU3AwYs/wj+1xa2/WZ1dOJkLE944uLi+Prrr5k8eTKBgYEkJCSQkJDAqVOngDPJTnp6Op9//jkpKSmOY/Ly8gDo3r07DRs2ZOjQofz111/8+uuvPPPMM8TFxeHt7W3l7VnKz8uDiCDz/nerlEdEXJV/RbjhY7jtJwipBsn7YfK/4Ie7IO3IxT8v5YLlCc9HH31EcnIynTt3JioqyrF8++23AKxZs4bly5ezfv16ateuXeCY/HY37u7uxMfH4+7uTmxsLLfddhu33347L730kpW35hTyq7X2qB2PiLi62l3hwaXQ7iGwucGGH80pKtZ8pSkqxPnG4bGKK47DA/D4D3/z7ar9PNK1DiOvq2t1OCIipePQWpj+MCT8ba5Xvwb6vQthtayNS4pdmR2HR4pX9dNj8Wi0ZREpV6Kbwz3z4bqXwMMX9vyhLuzlnBIeF1fjdMPlPcczLI5ERKSUuXtA+0fMaq5/dmE//LfV0UkpU8Lj4hwlPGrDIyLlVWFd2D/tAn++B3a71dFJKVHC4+JiQs0SnuRTOZzMyLY4GhERi5zdhb1eH7DnwJxn4cvrIfmA1dFJKVDC4+LO7pq+V9VaIlLe+VeEmyeZDZg9/c607dnwo9WRSQlTwlMOVAtTw2UREQebDVrcAfcvhuirITPZHLPnp/sg0/UmkhaTEp5yoHqYWa2lEh4RkbOE1YLhv0HH0ea4PX9PgXHtYe9SqyOTEqCEpxxQCY+IyHm4e0KXp+HOWeYozSf3wYTeMPdldV93MUp4yoEzk4iqhEdEpFAxbc0qrqa3gGGHP96Ez6+DYzusjkyKiRKecqDa6SotdU0XEbkAnyAYOA5uHA8+IeZozR9fA+u+sToyKQZKeMqB/ITneHo2aVm5FkcjIuLkGt8ADyyBGh0hJwOm3Q/xoyA3y+rI5Aoo4SkHAn08CfP3AmCv2vGIiFxccGUY+jN0egKwwarPYXxvSD5odWRymZTwlBMx6qklIlI0bm5w7ZNw63fgEwwHV8HHHWHXQqsjk8ughKecqBaqhEdE5LLU7Q73LoTIJpBxDL4aAIvfAcOwOjIpAiU85UQ1R08tVWmJiBRZaA0YPgea3mr24vr9efj2Ng1UWIYo4SknHD21lPCIiFweT18Y8D/o+za4ecKWePj0Wjiy2erI5BIo4Skn8kt49qlKS0Tk8tls0PIuuGs2BFWG4zvg066ai6sMUMJTTuSX8BxOySQzJ8/iaEREyrgqLeG+Rae7rqebc3H99izY7VZHJuehhKecCPP3IsDbA8OAAydUyiMicsX8K8JtU6HDSHN9yXvwwx2Qc8rSsKRwSnjKCZvNRox6aomIFC93D+j2AtzwKbh7waafYeL1kH7M6sjkH5TwlCNnGi4r4RERKVZXDYahU83xeg6sMOfhOr7T6qjkLEp4ypH8wQf3qaeWiEjxq97B7LoeEgNJu+CzbrBvmdVRyWlKeMoRx6zpSSrhEREpEZXqwd1zIfpqOJVkVm9tnGp1VIISnnIlf7RldU0XESlBAeFwRzzU6wN5WfD9HfDnuxqZ2WJKeMqRahXNEp79JzLIs+svnohIifHyh5u+gjb3m+tznoNfRkFerrVxlWNKeMqRyCAfvNzdyMkzOHRS3SZFREqUmzv0egN6vo454/oXMOVWyFYpuxWU8JQj7m42qoT6AuqaLiJSato+YJb2ePjC9l9h8mDISrM6qnJHCU85c6bhsnpqiYiUmgb94PZp4BUIe/6Ar2+AzGSroypXlPCUMzFquCwiYo2YtnD7z+ZYPfuXw5f9ISPJ6qjKDSU85Ux1zZouImKdKi1gWDz4hcGhtRqVuRRZnvCMGTOGVq1aERgYSHh4OAMGDGDr1q0FjsnMzCQuLo6wsDACAgIYNGgQiYmJBY7Zt28fffr0wc/Pj/DwcB577DFyc9Ua/p/yZ01XGx4REYtEXQV3/AIBEZC4Hsb3htQEq6NyeZYnPAsXLiQuLo5ly5YxZ84ccnJy6N69O+npZ0ogRo4cyYwZM/j+++9ZuHAhhw4d4oYbbnDsz8vLo0+fPmRnZ7NkyRImTpzIhAkTeO6556y4JaeWP73E3uMZGBoTQkTEGuEN4I6ZEFQZjm2F8b0g+YDVUbk0m+Fkv3pHjx4lPDychQsX0rFjR5KTk6lUqRKTJ0/mxhtvBGDLli00aNCApUuX0rZtW2bNmkXfvn05dOgQERERAIwbN47HH3+co0eP4uXlddHrpqSkEBwcTHJyMkFBQSV6j1bKzrVT/9lZ2A1Y8VRXwoN8rA5JRKT8OrEHJvaDk/vMKSnunAXBVayOqky51N9vy0t4/ik52Wy1HhoaCsDq1avJycmhW7dujmPq169PTEwMS5cuBWDp0qU0adLEkewA9OjRg5SUFDZu3FjodbKyskhJSSmwlAdeHm5Eh5hd0zWJqIiIxSpUN5Oc0Fpm0vPlALXpKSFOlfDY7XYeffRR2rdvT+PGjQFISEjAy8uLkJCQAsdGRESQkJDgOObsZCd/f/6+wowZM4bg4GDHUrVq1WK+G+eV3zVdDZdFRJxAcBWz91ZQFTi+Hb4aqC7rJcCpEp64uDg2bNjAlClTSvxaTz75JMnJyY5l//79JX5NZ3GmHY8SHhERpxBS1Ux6/CpCwt8w+WaNyFzMnCbhGTFiBPHx8cyfP58qVc7UX0ZGRpKdnc3JkycLHJ+YmEhkZKTjmH/22spfzz/mn7y9vQkKCiqwlBc1KuaX8Ogvk4iI06hYG4ZOBe9g2LcEvrsdcrOtjsplWJ7wGIbBiBEjmDp1KvPmzaNGjRoF9rdo0QJPT0/mzp3r2LZ161b27dtHbGwsALGxsaxfv54jR444jpkzZw5BQUE0bNiwdG6kDDnTNV0lPCIiTiXqKhjynTkNxY45MPU+sOdZHZVL8LA6gLi4OCZPnszPP/9MYGCgo81NcHAwvr6+BAcHM3z4cEaNGkVoaChBQUE89NBDxMbG0rZtWwC6d+9Ow4YNGTp0KGPHjiUhIYFnnnmGuLg4vL29rbw9p5Q/+ODeY2bXdJvNZnFEIiLiENMWbv7arNba+BN4B0K/d0H/Vl8Ry0t4PvroI5KTk+ncuTNRUVGO5dtvv3Uc8/bbb9O3b18GDRpEx44diYyM5KeffnLsd3d3Jz4+Hnd3d2JjY7ntttu4/fbbeemll6y4JadXNdQPmw1Ss3JJSldxqYiI06ndDQZ9CjY3WDMR5r5odURlntONw2OV8jIOT752Y+ZyKDmTHx+IpUW1UKvDERGRwqz5CqaPMN9f/z5cfbu18Tghy8bhSUvTlPdlQfXTDZd3H1PDZRERp3X1UOj0hPk+fhTsWWxtPGVYkRKet99++4L7U1NT6dGjxxUFJKUjP+HZc0wNl0VEnFqnx6HRQLDnwLe3QdIuqyMqk4qU8Dz11FN8+eWXhe5LT0+nZ8+eHD9+vFgCk5JVIyy/hEcJj4iIU3NzgwEfQfTVcOoETL4JTp20Oqoyp0gJz1dffcV9993H9OnTC2xPT0+nR48eHD16lPnz5xdrgFIyalRUwiMiUmZ4+sIt35yebHQb/HAn5OVaHVWZUqSE58Ybb+T999/nlltuYcGCBcCZkp3ExEQWLFhAVFRUScQpxcxRpXU8XbOmi4iUBYGRZtLj6Qc758GvT1odUZlS5EbLd999N88//zz9+/dnwYIF9OrVi0OHDjF//nyio6NLIkYpATGhfrjZICM7j6OpWVaHIyIilyKqKdzwqfl+xSew4lNr4ylDLquX1ujRo3nggQfo2rUrBw8eZMGCBQWmgxDn5+XhRuUK5qzpu1StJSJSdjToC12fN9/PfgL2r7Q2njKiSCMt33DDDQXWPT09qVixIo888kiB7WcPCijOq0bFAPYnnWLPsXTa1gyzOhwREblUHUaak4xunArf3wH3/wF+GlPtQoqU8AQHBxdYv+WWW4o1GCldNSv6s2jbUTVcFhEpa2w26PceHP4bknaac27d8q3Zo0sKVaSEZ/z48SUVh1igZiWz4fLOo0p4RETKHJ8gGDwRPusG23+DJe+aJT9SKKWC5diZrukaHVtEpEyKbAK93jDfz30Z9i6xNh4npoSnHKtZKQCAfUkZ5ObZLY5GREQuy9XD4KqbwMiDH+6CtKNWR+SUlPCUY1FBPvh4upGTZ7D/xCmrwxERkcths0Gft6BiXUg9DDMeAY2vdg4lPOWYm5uNGhXNUp5dR1WtJSJSZnkHwI3jwc0Ttv4Cf31jdUROp0gJz3PPPcfq1atLKhaxQK3TDZd3qeGyiEjZFtkYrj09+vKsx+HkfmvjcTJFSngOHDhAr169qFKlCg888ACzZs0iOzu7pGKTUlDrdDuenSrhEREp+9o9AlVaQVYK/BwHdrXPzFekhOeLL74gISGBb775hsDAQB599FEqVqzIoEGD+PLLL0lKSiqpOKWEnOmaroRHRKTMc/eAgR+Dhy/sXggrP7M6IqdR5DY8bm5uXHPNNYwdO5atW7eyfPly2rRpw8cff0x0dDQdO3bkzTff5ODBgyURrxSzMyU8qtISEXEJYbXgupfM93Oeg2M7rI3HSVxxo+UGDRowevRo/vzzT/bv38+wYcP4448/+OYbNZgqC/JLeJLSs0lKV/WkiIhLaHU31OgEuadg+ghVbVHMvbQqVarE8OHD+fnnn/nPf/5TnKeWEuLn5UHlEHMS0R1HVK0lIuIS3Nyg/4fg6Qf7lsJfk62OyHLqli7UiTCrtbYfSbU4EhERKTYhVaHzE+b7356FjPLdzrZICc/vv/+OocGMXE7t0+14VMIjIuJi2j4I4Q3hVJLZnqccK1LC06NHD44e1ZDVria/hEcJj4iIi3H3hL5vm+/XfgX7llkbj4WKlPCodMc11Q4PBGB7ohIeERGXE9MWmt9mvo8fBXk51sZjEbXhEWqHmyU8CSmZpGaWz78IIiIurdtL4FsBjmyE5R9bHY0lipzwfPTRR8ydO5cTJ06URDxigWBfTyKCvAHYrmotERHX4x8G3V403y8aWy4bMBc54fnggw+47rrrqFixItWrV+eGG27glVdeYebMmSQkJJREjFIK6kaY1VrbEtRTS0TEJTW/DSIaQ2YyLPqv1dGUuiInPBs3buTAgQNMnz6d4cOHYxgGn376KX379qVy5cpUrly5JOKUElbvdMKzNVEJj4iIS3Jzh+4vm+9XfApJu6yNp5R5FOVgm80GQHR0NNHR0fTp08ex7/jx46xevZp169YVa4BSOupGni7hUcIjIuK6anWB2t1gx+/w+4sweKLVEZWaYuulFRYWRvfu3Rk9evQVByWlz1HCoyotERHXdt1LYHODTdNg/wqroyk1RUp4Zs+eTXBwcLEGsGjRIvr160d0dDQ2m41p06YV2J+WlsaIESOoUqUKvr6+NGzYkHHjxhU4JjMzk7i4OMLCwggICGDQoEEkJiYWa5yurk5EADYbHEvL5lhaltXhiIhISYloBM2GmO9/f9HaWEpRkRKe7t274+3tXawBpKen07RpUz788MNC948aNYrZs2fz9ddfs3nzZh599FFGjBjB9OnTHceMHDmSGTNm8P3337Nw4UIOHTrEDTfcUKxxujo/Lw+qhfoBKuUREXF5nZ8Ady/Yuxh2/2F1NKXC8nF4evXqxSuvvMLAgQML3b9kyRKGDRtG586dqV69Ovfeey9NmzZlxQqzGC45OZnPP/+ct956iy5dutCiRQvGjx/PkiVLWLas/I4oeTnqRwYBsPlwisWRiIhIiQquAlffbr5f8Lq1sZQSyxOei2nXrh3Tp0/n4MGDGIbB/Pnz2bZtG927dwdg9erV5OTk0K1bN8dn6tevT0xMDEuXLj3vebOyskhJSSmwlHf1o8x2PJsPq4RHRMTldRhVrkp5nD7hef/992nYsCFVqlTBy8uLnj178uGHH9KxY0cAEhIS8PLyIiQkpMDnIiIiLjgu0JgxYwgODnYsVatWLcnbKBMaRJklPFsSlPyJiLi84Mpw9TDz/YIx4OLTR5WJhGfZsmVMnz6d1atX83//93/ExcXx+++/X9F5n3zySZKTkx3L/v37iynisqvB6Sqt7Ylp5OTZLY5GRERKXIeRp0t5/oQ9i62OpkQVaRye0nbq1Cmeeuoppk6d6hjz56qrrmLdunW8+eabdOvWjcjISLKzszl58mSBUp7ExEQiIyPPe25vb+9ib4Bd1lWp4EugtwepWbnsPJrmaNMjIiIuKrgyNB8Kqz6HJe9BjWusjqjEOHUJT05ODjk5Obi5FQzT3d0du90sgWjRogWenp7MnTvXsX/r1q3s27eP2NjYUo23rHNzszmqtTYeVLWWiEi5EBsH2GD7b3Bki9XRlBjLS3jS0tLYsWOHY3337t2sW7eO0NBQYmJi6NSpE4899hi+vr5Uq1aNhQsX8uWXX/LWW28BEBwczPDhwxk1ahShoaEEBQXx0EMPERsbS9u2ba26rTKrYXQQK/YksfFQCoNaWB2NiIiUuLBa0KAvbJ4BSz+A/h9YHVGJsDzhWbVqFddee61jfdSoUQAMGzaMCRMmMGXKFJ588kmGDBlCUlIS1apV49VXX+X+++93fObtt9/Gzc2NQYMGkZWVRY8ePfjf//5X6vfiChpFny7hOZRscSQiIlJq2j1sJjx/fwtdnoXACKsjKnY240LzRZQjKSkpBAcHk5ycTFBQ+W27sulQCr3f+4NAHw/+fr67Y/40ERFxcZ93h/3L4Zp/Q9fnrI7mkl3q77dTt+GR0lcnIgAvDzdSM3PZezzD6nBERKS0xI4wX1eNh1zXm2JICY8U4Onu5mi4vP6gqrVERMqNer0hqDKcSoJN0y9+fBmjhEfOcVVlc4LYvw+ctDYQEREpPe4eZwYiXPWFtbGUACU8co4mVfITHpXwiIiUK1cPBZs77FsCRzZbHU2xUsIj52haJQSADQeTsdvVpl1EpNwIioZ6vcz3q8ZbG0sxU8Ij56hVyR9fT3fSs/PYdSzN6nBERKQ0tbzLfP3rG8g5ZW0sxUgJj5zDw92NJqfb8azdd9LaYEREpHTVvBaCYyArBbbOtDqaYqOERwrVPCYEgLX7T1oah4iIlDI3N7hqsPn+rynWxlKMlPBIofITnjV7T1gbiIiIlL6mN5uvO+ZC2hFrYykmSnikUM1jKgCwLTGVtKxci6MREZFSVbEOVG4JRh6s/8HqaIqFEh4pVESQD5VDfLEb8LeqtUREyp/8Up6/vrE2jmKihEfOy1GttU/VWiIi5U6jG8DNAxL+huM7rY7miinhkfNqVT0UgBV7lPCIiJQ7/mFQ/Rrz/aafrY2lGCjhkfNqWd1sx7Nm7wnyNAChiEj507C/+aqER1xZ/cggArw9SMvKZUtCitXhiIhIaavfF2xucHgdnNhrdTRXRAmPnJe7m42rq5mlPKtUrSUiUv4EVIJq7c33m8v2DOpKeOSCWp1OeJbvPm5xJCIiYokG15uvW36xNo4rpIRHLqhtrTAAlu9KwjDUjkdEpNyp19N83b8cMpKsjeUKKOGRC2paJQRfT3eOp2ezLVETiYqIlDshMVCpARh22DnP6mgumxIeuSAvDzdHb62lO49ZHI2IiFiibnfzdduv1sZxBZTwyEW1rWlWay3dpXY8IiLlUp0e5uuO38GeZ20sl0kJj1xU7Ol2PMt2JWk8HhGR8qhqG/AJhlNJcGit1dFcFiU8clFXVQ4m0MeD5FM5rD+YbHU4IiJS2tw9zoy6vGuBpaFcLiU8clEe7m60r1URgEXbjlocjYiIWKJmZ/N190JLw7hcSnjkklxT10x4/tiuhEdEpFyq0cl83bccck5ZG8tlUMIjl6RjnUoArNl3ktTMHIujERGRUlexDgRGQV6WOSZPGaOERy5J1VA/alT0J89u8OcOdU8XESl3bDao0dF8v/sPa2O5DEp45JJ1rmeW8szbcsTiSERExBLV2pmvKuERV9a1fgQA87Ycxa7u6SIi5U9MrPl6YBXkla3mDUp45JK1rhGKv5c7x9Ky1D1dRKQ8qlgXfEMh9xQc/tvqaIrE8oRn0aJF9OvXj+joaGw2G9OmTTvnmM2bN3P99dcTHByMv78/rVq1Yt++fY79mZmZxMXFERYWRkBAAIMGDSIxMbEU76J88PJw45rTjZfnqlpLRKT8sdkgpq35ft8Sa2MpIssTnvT0dJo2bcqHH35Y6P6dO3fSoUMH6tevz4IFC/j777959tln8fHxcRwzcuRIZsyYwffff8/ChQs5dOgQN9xwQ2ndQrnStUE4AL9tTLA4EhERsUTVNubrgZXWxlFEHlYH0KtXL3r16nXe/U8//TS9e/dm7Nixjm21atVyvE9OTubzzz9n8uTJdOnSBYDx48fToEEDli1bRtu2bUsu+HKoW4MI3N1sbElIZc+xdKpX9Lc6JBERKU2VrzZfD5atKSYsL+G5ELvdzi+//ELdunXp0aMH4eHhtGnTpkC11+rVq8nJyaFbt26ObfXr1ycmJoalS5ee99xZWVmkpKQUWOTiKvh70bZmKACzVcojIlL+RDUDbJC8D9LKzmC0Tp3wHDlyhLS0NF5//XV69uzJb7/9xsCBA7nhhhtYuNAc2johIQEvLy9CQkIKfDYiIoKEhPP/II8ZM4bg4GDHUrVq1ZK8FZfSs3EUALM3KOERESl3fILMxssAh9ZYG0sROHXCY7fbAejfvz8jR46kWbNmPPHEE/Tt25dx48Zd0bmffPJJkpOTHcv+/fuLI+RyoUfDCGw2WLf/JIdOlr3hxUVE5Ao5qrVWWxtHETh1wlOxYkU8PDxo2LBhge0NGjRw9NKKjIwkOzubkydPFjgmMTGRyMjI857b29uboKCgAotcmvAgH1pWqwDAL38ftjgaEREpddGnE57Df1kbRxE4dcLj5eVFq1at2Lp1a4Ht27Zto1q1agC0aNECT09P5s6d69i/detW9u3bR2xsbKnGW55c36wyANPWHbQ4EhERKXWRjc3XhA3WxlEElvfSSktLY8eOHY713bt3s27dOkJDQ4mJieGxxx7jpptuomPHjlx77bXMnj2bGTNmsGDBAgCCg4MZPnw4o0aNIjQ0lKCgIB566CFiY2PVQ6sE9WkSxYvTN7LxUAo7jqRSOzzQ6pBERKS0RDQyX1MOQEYS+IVaG88lsLyEZ9WqVTRv3pzmzZsDMGrUKJo3b85zzz0HwMCBAxk3bhxjx46lSZMmfPbZZ/z444906NDBcY63336bvn37MmjQIDp27EhkZCQ//fSTJfdTXoT6e9GxrjkI4c/rDlkcjYiIlCqfYAiOMd8nbrQ2lktkMwxDkyIBKSkpBAcHk5ycrPY8l+jndQd5ZMo6YkL9WPhYZ2w2m9UhiYhIafnmFtg6E3q+Dm0fsCyMS/39tryER8qu6xpG4Oflzr6kDFbtPWF1OCIiUpoiTrfjSSwb7XiU8Mhl8/PyoO9V5pg83yzfd5GjRUTEpYTXN1+Pbbc2jkukhEeuyK1tzN5y8esPczIj2+JoRESk1OQPPnh0K5SB1jFKeOSKNK0STIOoILJz7Uxdqy7qIiLlRmgtwAaZJyHjuNXRXJQSHrkiNpuNW1ub03J8s2IfagMvIlJOePlByOlpmY5tszaWS6CER65Y/+aV8fF0Y1tiGmv2qfGyiEi5kV+tpYRHyoMgH0/6XRUNwOTlmpNMRKTccCQ8zt9wWQmPFItb2pgDUM34+xDH0rIsjkZEREpFWG3zVSU8Ul40rxpC0yrBZOfa+WrpXqvDERGR0hBa03w94fz/7ivhkWJhs9m4p6P5B/+rZXvJzMmzOCIRESlxFaqbryf2gN1uZSQXpYRHik3PRpFUqeBLUno2P645YHU4IiJS0oKrgs0d8rIgLcHqaC5ICY8UGw93N4Z3qAHAZ3/sxm5XF3UREZfm7nGma/qJPZaGcjFKeKRYDW5ZlSAfD3YfS+fXjc6d7YuISDE4u1rLiSnhkWLl7+3BHe2qA/Du3O0q5RERcXXBp0t4kp27KYMSHil2wzvUJNDbgy0JqcxWKY+IiGtzJDzOPQ6bEh4pdsF+ntx5ui3PO79vUymPiIgrC65svqqER8qj4R1qEOjjwbbENH5Zf9jqcEREpKQEVzFflfBIeRTs6+nosfV/v20lO9e5x2cQEZHLFJSf8BwEJ55AWgmPlJi7r6lJxQBv9hzP4Otlzj8Kp4iIXIbASPM1Jx2yUq2N5QKU8EiJCfD24N/dzYnl3p27nZMZ2RZHJCIixc47ALyDzPepzttRRQmPlKjBLatSLyKQ5FM5vD9vh9XhiIhIScgv5Ul13jabSnikRLm72Xi6TwMAvly6hx1H0iyOSEREil1glPmqhEfKs451K9Glfjg5eQZPT12P4cSN2kRE5DIERJivaUesjeMClPBIqXjx+kb4eLqxfHcSP6x27q6LIiJSRAHh5mu6Eh4p56qG+jGym9mA+bWZm0lKVwNmERGX4V/RfE07am0cF6CER0rNXR1qUD8ykBMZObwSv8nqcEREpLj4q4RHxMHT3Y0xNzTBZoOf1h5k9gbnbdwmIiJF4KjSUgmPCADNYypwX8daADz503qOpGRaHJGIiFwxv1DzNeOEtXFcgBIeKXWjrqtLw6ggTmTkMPrHv9VrS0SkrPMLM18zjlsbxwUo4ZFS5+Xhxjs3N8PLw40FW49q2gkRkbIuP+HJPQXZGdbGch6WJzyLFi2iX79+REdHY7PZmDZt2nmPvf/++7HZbLzzzjsFticlJTFkyBCCgoIICQlh+PDhpKVpgDtnVjcikCd61gfg5V82s/5AssURiYjIZfMKAHcv872TlvJYnvCkp6fTtGlTPvzwwwseN3XqVJYtW0Z0dPQ5+4YMGcLGjRuZM2cO8fHxLFq0iHvvvbekQpZicke76nRrEE52rp37v16truoiImWVzQY+Ieb7U87ZjsfyhKdXr1688sorDBw48LzHHDx4kIceeohJkybh6elZYN/mzZuZPXs2n332GW3atKFDhw68//77TJkyhUOHDpV0+HIF3NxsvHVTM2pU9OfgyVM89M0acvPsVoclIiKXw7eC+Zp50tIwzsfyhOdi7HY7Q4cO5bHHHqNRo0bn7F+6dCkhISG0bNnSsa1bt264ubmxfPny8543KyuLlJSUAouUviAfTz4e2gI/L3f+3HGc//621eqQRETkcviGmK+nTloZxXk5fcLzxhtv4OHhwcMPP1zo/oSEBMLDwwts8/DwIDQ0lISE809TP2bMGIKDgx1L1apVizVuuXR1IwL5741NAfh44S5+WqOpJ0REypz8Eh5VaRXd6tWreffdd5kwYQI2m61Yz/3kk0+SnJzsWPbv31+s55ei6XNVFA90NsfnGf3D3yzY6ryjdYqISCF8gs3XTOfshOLUCc8ff/zBkSNHiImJwcPDAw8PD/bu3cu///1vqlevDkBkZCRHjhT8cczNzSUpKYnIyMjzntvb25ugoKACi1jrse71GNAsmly7wQNfr2Hd/pNWhyQiIpfK+/TvaJZzNhFx6oRn6NCh/P3336xbt86xREdH89hjj/Hrr78CEBsby8mTJ1m9erXjc/PmzcNut9OmTRurQpfL4OZmY+yNTbmmTkVO5eRx14SV7Dyq4QVERMoEn9MJT6ZzJjweVgeQlpbGjh07HOu7d+9m3bp1hIaGEhMTQ1hYWIHjPT09iYyMpF69egA0aNCAnj17cs899zBu3DhycnIYMWIEN998c6Fd2MW5eXm4Me62Ftzy6TL+PpDM7Z+v4IcHYokK9rU6NBERuRCV8FzYqlWraN68Oc2bNwdg1KhRNG/enOeee+6SzzFp0iTq169P165d6d27Nx06dOCTTz4pqZClhPl7e/DFHa2oHubHwZOnGPzxUvYdd86RO0VE5DQnL+GxGZrICICUlBSCg4NJTk5Wex4nceBEBkM+W87e4xmEB3oz6e421IkItDosEREpzN/fw093Q41OMGx6qV32Un+/LS/hETmfKhX8+P6+WOpFBHIkNYvBHy/VFBQiIs7Ky998zXbOtpdKeMSphQf5MOXetjStEsyJjBxu/XQZK3YnWR2WiIj8k3eA+Zqdbm0c56GER5xeBX8vvr67Da1rhJKalcvtXyxnvsbpERFxLvklPFkq4RG5bIE+nky8szWd61UiM8fO8AkrGbdwJ2qCJiLiJLxOt7HMTrU2jvNQwiNlhq+XO58MbcngllWwG/D6rC2M+GYt6Vm5VocmIiKep4cPyTllbRznoYRHyhQvDzfeGHQVrwxojKe7jV/+PswN/1vCnmPOWWcsIlJu5Fdp5WVDnvP9R1QJj5Q5NpuN29pW45t72lIp0Jutialc/8FitesREbGSp9+Z9znON3aaEh4ps1pWDyX+oQ5cHRNCSmYud45fyQvTN5KR7Xz/sxARcXke3sDpib6V8IgUr4ggH765ty23x1YDYMKSPfR+9w9W71XXdRGRUmWzOXU7HiU8UuZ5e7jzUv/GfHlXa6KCfdhzPIN/jVvKmFmbyczJszo8EZHyw8PHfM3NtDaOQijhEZfRsW4lZj/akUFXm724Pl64i+s/WKyBCkVESotKeERKR7CvJ/83uCmf3t6SigHebEtMY/DHS3nom7UcOul8fwFFRFyKSnhEStd1DSOYM7Ijt7aJwWaDGX8dosv/LeC9udtVzSUiUlIcCU+WtXEUQgmPuKwK/l68NrAJ8Q91oHX1UDJz7Lw1Zxtd/28hP605QG6e3eoQRURci4eX+ZqXbW0chVDCIy6vUXQw397XlvdvaU50sA8HT55i1Hd/0f3tRUxde4A8u6anEBEpFu7e5qtKeESsYbPZ6Nc0mrn/7szonvWo4OfJrmPpjPz2L657eyHT1h5U4iMicqXcPc1XlfCIWMvXy50HO9fmj8e78FiPeoT4ebLraDqPfruOLv+3gM8X7yY1M8fqMEVEyiZ3VWmJOJUAbw/irq3N4rMSn73HM3g5fhNtX5vLC9M3slvzc4mIFE1+wmN3vhHvPawOQMRK+YnPne2rM3XtQcb/uYcdR9KYsGQPE5bsoXO9StzUsipdGoTj7eFudbgiIs7N7fS/k3nOV1KuhEcE8PPyYEibatzaOobFO44x/s89zNtyhAVbj7Jg61FC/Dzpd1U0g1pUoWmVYGw2m9Uhi4g4n/w2PE5YwmMzDEMtNYGUlBSCg4NJTk4mKCjI6nDECew+ls53q/bz05oDJKac6XFQOzyAAc2i6dk4ktrhgQAs/fEblnw3mfaDh9B20M2O9XaDbyV20C1W3YKIlBMpc/eRMmcvQddVI6hrzDnrJW3lL7tZMWM3bWr+RcuMF1gZ+Qkr1lWidb8atOpTo0Svfam/30p4TlPCI+eTZzf4c8cxflxzgF83JpCZc2b8npqV/OmZsx73tb86tsU0bsq+DX851tsNHqKkR6SYGIZBjj2H7LxssvKyCrzPtmeTk5dDjj2HnLwcco1cc904vW7PJceeQ56Rh92wk2vPxW7YyTPyHO/zF8B8jx3DMMyFgj+XttMzg7vZ3LDZbLjZ3HCzueFuc8dms+Fuc8fd5o6Hmwcebh6O955unni6eeLl7mW+d/fE293bsfi4++Dl7oWPhw9+nn54unle8JnkJzf5vGuHkLXjpGO9pJOe/GQnX0TuBhI9GjvWSzrpUcJTREp45FKkZuYwa30CszYc5s8dx8nOszNi90fkV3Aln8okOSOTmLCQsz5l49/fzrAgWpHSl5OXQ3pOOhm5GWTkZHAq95Rjycj9x/pZ+zPzMsnKzeJU3imycrPIysviVO4psvKyCmzPtjtf75+S5uHmga+HL74evvh5+BHgGYC/pz/+nv4EeAVwf3wPR/K1fP9fxIREExVYqcA5qrx+TYnF9+H98wpuMAxz5vSzxI3rUmLXv9Tfb7XhESmCQB9PBreqyuBWVUnNzGHB1qMs+6kLoVvmkZGdw8sz5gLwTN8uhPiZk+hV6TaQzJw8fDzV6FmcU549j/TcdNKy00jLSXO8puekk5qdSnpOumM9f8nIyTDf5571PiedHHvpNlbNLynxdvd2vPdy83KUqni6eRb63t3mjrubu6MUJv+9m80NG7YzpTaYrzZsUFjTPaNgSVB+6VH+kmvPdZQg5S/5JU75pVPZ9mxHKVV+speVl0WeYU6Dk2vPJTU7ldTs1EKfgX/FHG4/1o/7pj7LzG0L6VyjDV8N/q9j/5eVZjBrylOE+oQS5htGmE8YFX0rUsmvEpV8KxHhF0GEfwQRfhH45E8NUQSt+9UoUMLzz2SnzfUlW6V1qZTwiFymQB9P+jWNpl/TUXz70jHenvi1Y192bh4GsN+nCu/vjOC5F3+jWdUQWlSr4FhC/LysC15chmEYZOVlkZqdSkp2SoHX1OxU0nLSSMlOMZOY7DRScsz3qdmp5mtOKqdyi39iXS83L/w8/QqUTPh6+hZc9ziz7uPhg4+7D94e3mfeu3vj6+FrVvV4eDu2ebl7OaqD3GyuO7pKTl6Oo1TM8XpWcpmfkKblpLEpfi8zty0EYOuxXQAYGKz138I3FWdBFpzMOsmu5F0XvGaoTyhR/lFEB0RTJbAKVQKqUC2oGtWCqhHhF1Foh41WfWpwaPtJDmw5cc6+KvUr0LK3Eh4Rl7D0x29YuGA+y3ftB+DBa2MJDwoAICbzAJ0z1rHArxnLdyexfHeS43M1K/rTpEowTSqbS6PKwQR4669keWQ37GZikpVCSvbp5az3qdmp56yfndgUV6mKt7s3/p7+BHoFmtUlngHm4hXgWPfz9HNUp/h5+J1Z9zi9zdPvktqdyMV5unsS7B5MsHfwBY9LmbuPB38aAUDD8NpMHzoOMNsYXZ3egDmR35Pa2p2kzCSSMpM4duoYR08d5VjGMY5kHOHIqSMkpCdwKveU45iNxzeecx1fD1+qB1WnVkgtaofUpm6FutQPrc+e+WmFJjsAB7acYNXM3U6R9OhfV5ErtHDyV/ywaj0AbWvGULNSaIH9TRKX8cz7j7B6bxKr955g1d4T7Dqazq5j5vLzukOOY6uF+VEvIpD6UUHUjwykbkQA1cL88XR33f/FugrDMBylKfnJSXJWcoEEJjk7uUDikr8/LTvtnAaxReVmcyPQK5BAz0ACvQIJ8goy170CCfAKcOwr7H1+YuPpriSlLNr84wq+XT8TgJevexRvj4Klx7nzj1Knx4Xb8BiGQUp2CofTD3Mo7RAH0w5yIPUA+1P3sz91PwdSD3Aq9xSbkzazOWlzgc/et/QdRxui0ycrUK21fLoSHhGXsPpUHkdT0wn08ab3VfWJadKMfevXOfa3G3wrtcMDqB0ewE2tzJ4SJ9Kz+ftgMusPnOTvA8lsOJjMoeRM9h7PYO/xDH7blOj4vKe7jZhQP2pWCqBmJX9qVQygeUwIdSICS/tWy6Xjp47zx8E/SM5KPpPAnF0SczpxSc1OdbS5uFw+7j4EeQUR5B1kvp5OWvLXz05k/vnq7+mv8aHKqYmJs8i159GhWgtaV7mq0F5aF2Oz2Qj2NkuT6ofWP7Mj5RDsnEeOkceBrJPsyjrOzsxjbM88yrbMo+zJSmJl1Zm03t/H8ZEKJ7ZwIrSBY71Vr6rFcp9XSr20TlMvLbkcM2fOpE8f8y/63de0ZvgjIy97HJ7jaVlsTUhlS0IqWxJS2JKQyo4jaWRkn/sjen+nWjzRq34hZ5HilJCewM3xN3M88/glf8bb3duRrAR7BxdMYM5KZBz7ztru5a52XVI0WVlZREZGcvLkSb7615tcf9+/inccnp3z4KuB5919ymZjq5cnq9IGk3n8BmrujqfG3tnsrtaT3dX7UmOPuW7z9sY9OJioMa8R0L795cdTCHVLLyIlPFJUiYmJXHXVVRw5coRHHnmEd955p9ivYbcbHE7JZOeRNHYdTWPXsXR2H0vnltYx9G4SVezXk4Lm7J3D6IWj8fbwplOVTgWSlLOTmWCvYEfScjm9XEQu17Rp0xg4cCCVK1dm7969uLsXc2/Qw3/BvFcu+fBtb/xFXrYZgy0wAONUJuSeGXU5ZuJE/Nu0LtYQL/n327DYwoULjb59+xpRUVEGYEydOtWxLzs72xg9erTRuHFjw8/Pz4iKijKGDh1qHDx4sMA5jh8/btx6661GYGCgERwcbNx1111GampqkeJITk42ACM5Obk4bktcnN1uN3r16mUARpMmTYxTp05ZHZKUkBWHVxgHUg9YHYZIoYYOHWoAxqhRo6wOpVB2u93ITU01svYfMDLWbzByU9OK/RqX+vtteUvI9PR0mjZtyocffnjOvoyMDNasWcOzzz7LmjVr+Omnn9i6dSvXX399geOGDBnCxo0bmTNnDvHx8SxatIh77723tG5ByqG33nqLWbNm4e3tzeTJk/Hx0f/qXVWryFZUDqhsdRgi5zAMg/nz5wPQq1cvi6MpnM1mwz0gAK8qlfFt3Aj3AH/rYjEM56nSstlsTJ06lQEDBpz3mJUrV9K6dWv27t1LTEwMmzdvpmHDhqxcuZKWLVsCMHv2bHr37s2BAweIjo6+pGurSksu1eLFi+ncuTN5eXl8+OGHPPjgg1aHJCLl0M6dO6lduzaenp6cPHkSPz8/q0OyxKX+fltewlNUycnJ2Gw2QkJCAFi6dCkhISGOZAegW7duuLm5sXz58vOeJysri5SUlAKLyMUkJiYyePBg8vLyuOWWW3jggQesDklEyql169YB0KxZs3Kb7BRFmUp4MjMzefzxx7nlllscWVxCQgLh4eEFjvPw8CA0NJSEhITznmvMmDEEBwc7lqpVnaPbnDiv3NxcbrnlFg4fPkyDBg345JNP1A1YRCyzebM5Hk7Dhg0tjqRsKDMJT05ODoMHD8YwDD766KMrPt+TTz5JcnKyY9m/f38xRCmu7Omnn2b+/PkEBATw008/ERAQYHVIIlKObdmyBYD69TVExaUoEwMP5ic7e/fuZd68eQXq6CIjIzly5EiB43Nzc0lKSiIyMvK85/T29sbb27vEYhbXMmnSJMaOHQvA559/rn9gRMRy+bUYqqG4NE5fwpOf7Gzfvp3ff/+dsLCwAvtjY2M5efIkq1evdmybN28edrudNm3alHa44oJWrVrF3XffDZglg4MHD7Y4IhERSEoy5+YLDQ29yJECTlDCk5aWxo4dOxzru3fvZt26dYSGhhIVFcWNN97ImjVriI+PJy8vz5HRhoaG4uXlRYMGDejZsyf33HMP48aNIycnhxEjRnDzzTdfcg8tkfM5fPgwAwYMIDMzk759+/LKK5c+AJeISEk6ccKcsFMJz6WxvFv6ggULuPbaa8/ZPmzYMF544QVq1Ch8wrH58+fTuXNnwMxyR4wYwYwZM3Bzc2PQoEG89957RWpjoW7p8k9ZWVlce+21LF26lAYNGrBs2TL92RARp1GxYkWOHz/Ohg0baNSokdXhWOZSf78tL+Hp3LkzF8q5LiUfCw0NZfLkycUZlpRzhmHwwAMPOIY9+Pnnn5XsiIhT8fAwf8Lz8q5s0trywunb8IhY4c0332T8+PG4ubnx7bffUqdOHatDEhEpwNPTEzDbusrFKeER+Ycff/yR0aNHA+YUEt27d7c4IhGRcynhKRolPCJnWbFiBbfddhsAI0aM4OGHH7Y4IhGRwuUnPNnZ2RZHUjYo4RE5bc+ePfTr14/MzEx69+7N22+/rZGURcRp5U+xlN89XS5MCY8I5hxtffr04ciRIzRt2pQpU6Y4GgSKiDijqKgowBw+Qy5OCY+Uezk5OfzrX/9i06ZNREdHEx8fT2BgoNVhiYhckBKeolHCI+WaYRjExcUxZ84c/P39mTFjBlWqVLE6LBGRi1LCUzRKeKRce/PNN/n0009xc3Pjm2++4eqrr7Y6JBGRS5I/m8CBAwcsjqRsUMIj5dbZ3c/ffvtt+vXrZ3FEIiKXrm7dusCZWdPlwpTwSLm0dOlSdT8XkTKtQYMGAOzdu5eMjAyLo3F+Snik3Nm2bZuj+3nfvn15++23rQ5JRKTIKlWqRFhYGIZhsG3bNqvDcXpKeKRcSUxMpGfPnhw/fpxWrVqp+7mIlGn169cHYNOmTRZH4vyU8Ei5kZ6eTt++fdm9ezc1a9YkPj4ef39/q8MSEblsTZo0AWDdunXWBlIGKOGRciE3N5fBgwezatUqwsLCmD17NuHh4VaHJSJyRVq2bAnAypUrLY7E+SnhEZdnGAYPPvggM2fOxMfHh/j4eM1+LiIuoVWrVgCsXr0au91ucTTOTQmPuLxXX33VMdbOlClTaNu2rdUhiYgUi4YNG+Lr60tqaqoaLl+EEh5xaRMmTODZZ58F4L333qN///4WRyQiUnw8PDxo3rw5AMuWLbM4GuemhEdc1m+//cY999wDwOOPP05cXJzFEYmIFL9rrrkGgIULF1ociXNTwiMuad26dQwaNIjc3FxuvfVWXnvtNatDEhEpEZ07dwaU8FyMEh5xOXv37qV3796kpaVx7bXX8sUXX+Dmpj/qIuKa2rdvj7u7O7t372bv3r1Wh+O09CsgLiUpKYlevXpx+PBhGjduzE8//YS3t7fVYYmIlJjAwEBH93SV8pyfEh5xGZmZmQwYMIDNmzdTuXJlZs2aRUhIiNVhiYiUuPxqrQULFlgahzNTwiMuwW63c/vtt/PHH38QFBTErFmzqFKlitVhiYiUCiU8F6eER1zCY489xvfff4+npydTp051DLcuIlIenN2OZ/fu3VaH45SU8EiZ98477/DWW28B5rg7Xbp0sTgiEZHSFRgYSPv27QGYNWuWxdE4JyU8Uqb98MMPjBo1CoDXX3+dW2+91eKIRESs0bt3bwBmzpxpcSTOSQmPlFl//PEHt912m2OurNGjR1sdkoiIZXr16gXAvHnzyMzMtDga56OER8qkzZs3079/f7Kysujfvz/vvfceNpvN6rBERCzTpEkTKleuzKlTp9Q9vRBKeKTMSUhIoHfv3pw4cYI2bdowefJk3N3drQ5LRMRSNptN1VoXoIRHypT09HT69u3Lnj17qF27NjNmzMDPz8/qsEREnIISnvOzPOFZtGgR/fr1Izo6GpvNxrRp0wrsNwyD5557jqioKHx9fenWrRvbt28vcExSUhJDhgwhKCiIkJAQhg8fTlpaWinehZSG3Nxcbr75ZlavXk3FihWZNWsWlSpVsjosERGn0bVrVzw9PdmxYwfbtm2zOhynYnnCk56eTtOmTfnwww8L3T927Fjee+89xo0bx/Lly/H396dHjx4FGmQNGTKEjRs3MmfOHOLj41m0aBH33ntvad2ClALDMHj44YeJj4/Hx8eH6dOnU7t2bavDEhFxKoGBgXTq1AmA6dOnWxyNkzGcCGBMnTrVsW63243IyEjjv//9r2PbyZMnDW9vb+Obb74xDMMwNm3aZADGypUrHcfMmjXLsNlsxsGDBy/52snJyQZgJCcnX/mNSLF74403DMCw2WzGjz/+aHU4IiJO64MPPjAAo127dlaHUiou9ffbw8pk62J2795NQkIC3bp1c2wLDg6mTZs2LF26lJtvvpmlS5cSEhLimDgNoFu3bri5ubF8+XIGDhxY6LmzsrLIyspyrCcnJwOQkpJSQncjl+uHH37g8ccfB2DMmDF069ZN35OIyHl07doVgCVLlrBt2zYiIyMtjqhk5f8eGIZxweOcOuFJSEgAICIiosD2iIgIx76EhATCw8ML7Pfw8CA0NNRxTGHGjBnDiy++eM72qlWrXmnYUoKeeOIJnnjiCavDEBEpE+rVq2d1CKUmNTWV4ODg8+536oSnJD355JOOEXrBnHwyKSmJsLCwYh3PJSUlhapVq7J//36CgoKK7bxSkJ5zydMzLnl6xiVPz7jklfYzNgyD1NRUoqOjL3icUyc8+cVwiYmJREVFObYnJibSrFkzxzFHjhwp8Lnc3FySkpIuWIzn7e2Nt7d3gW0hISHFE3ghgoKC9JerFOg5lzw945KnZ1zy9IxLXmk+4wuV7OSzvJfWhdSoUYPIyEjmzp3r2JaSksLy5cuJjY0FIDY2lpMnT7J69WrHMfPmzcNut9OmTZtSj1lEREScj+UlPGlpaezYscOxvnv3btatW0doaCgxMTE8+uijvPLKK9SpU4caNWrw7LPPEh0dzYABAwBo0KABPXv25J577mHcuHHk5OQwYsQIbr755osWb4mIiEj5YHnCs2rVKq699lrHen67mmHDhjFhwgRGjx5Neno69957LydPnqRDhw7Mnj0bHx8fx2cmTZrEiBEj6Nq1K25ubgwaNIj33nuv1O+lMN7e3jz//PPnVJ9J8dJzLnl6xiVPz7jk6RmXPGd9xjbjYv24RERERMo4p27DIyIiIlIclPCIiIiIy1PCIyIiIi5PCY+IiIi4PCU8IiIi4vKU8BSDDz/8kOrVq+Pj40ObNm1YsWLFBY///vvvqV+/Pj4+PjRp0oSZM2eWUqRlV1Ge8aeffso111xDhQoVqFChAt26dbvodyKmov5ZzjdlyhRsNptjfCw5v6I+45MnTxIXF0dUVBTe3t7UrVtX/2ZcRFGf8TvvvEO9evXw9fWlatWqjBw5kszMzFKKtuxZtGgR/fr1Izo6GpvNxrRp0y76mQULFnD11Vfj7e1N7dq1mTBhQonHeY5SmLndpU2ZMsXw8vIyvvjiC2Pjxo3GPffcY4SEhBiJiYmFHv/nn38a7u7uxtixY41NmzYZzzzzjOHp6WmsX7++lCMvO4r6jG+99Vbjww8/NNauXWts3rzZuOOOO4zg4GDjwIEDpRx52VLU55xv9+7dRuXKlY1rrrnG6N+/f+kEW0YV9RlnZWUZLVu2NHr37m0sXrzY2L17t7FgwQJj3bp1pRx52VHUZzxp0iTD29vbmDRpkrF7927j119/NaKiooyRI0eWcuRlx8yZM42nn37a+OmnnwzAmDp16gWP37Vrl+Hn52eMGjXK2LRpk/H+++8b7u7uxuzZs0sn4NOU8Fyh1q1bG3FxcY71vLw8Izo62hgzZkyhxw8ePNjo06dPgW1t2rQx7rvvvhKNsywr6jP+p9zcXCMwMNCYOHFiSYXoEi7nOefm5hrt2rUzPvvsM2PYsGFKeC6iqM/4o48+MmrWrGlkZ2eXVohlXlGfcVxcnNGlS5cC20aNGmW0b9++RON0FZeS8IwePdpo1KhRgW033XST0aNHjxKM7Fyq0roC2dnZrF69mm7dujm2ubm50a1bN5YuXVroZ5YuXVrgeIAePXqc9/jy7nKe8T9lZGSQk5NDaGhoSYVZ5l3uc37ppZcIDw9n+PDhpRFmmXY5z3j69OnExsYSFxdHREQEjRs35rXXXiMvL6+0wi5TLucZt2vXjtWrVzuqvXbt2sXMmTPp3bt3qcRcHjjL757lU0uUZceOHSMvL4+IiIgC2yMiItiyZUuhn0lISCj0+ISEhBKLsyy7nGf8T48//jjR0dHn/IWTMy7nOS9evJjPP/+cdevWlUKEZd/lPONdu3Yxb948hgwZwsyZM9mxYwcPPvggOTk5PP/886URdplyOc/41ltv5dixY3To0AHDMMjNzeX+++/nqaeeKo2Qy4Xz/e6lpKRw6tQpfH19SyUOlfCIS3v99deZMmUKU6dOLTD/mlyZ1NRUhg4dyqeffkrFihWtDsdl2e12wsPD+eSTT2jRogU33XQTTz/9NOPGjbM6NJexYMECXnvtNf73v/+xZs0afvrpJ3755Rdefvllq0OTYqYSnitQsWJF3N3dSUxMLLA9MTGRyMjIQj8TGRlZpOPLu8t5xvnefPNNXn/9dX7//XeuuuqqkgyzzCvqc965cyd79uyhX79+jm12ux0ADw8Ptm7dSq1atUo26DLmcv4sR0VF4enpibu7u2NbgwYNSEhIIDs7Gy8vrxKNuay5nGf87LPPMnToUO6++24AmjRp4piw+umnn8bNTeUCV+p8v3tBQUGlVroDKuG5Il5eXrRo0YK5c+c6ttntdubOnUtsbGyhn4mNjS1wPMCcOXPOe3x5dznPGGDs2LG8/PLLzJ49m5YtW5ZGqGVaUZ9z/fr1Wb9+PevWrXMs119/Pddeey3r1q2jatWqpRl+mXA5f5bbt2/Pjh07HMkkwLZt24iKilKyU4jLecYZGRnnJDX5CaahubWLhdP87pVqE2kXNGXKFMPb29uYMGGCsWnTJuPee+81QkJCjISEBMMwDGPo0KHGE0884Tj+zz//NDw8PIw333zT2Lx5s/H888+rW/pFFPUZv/7664aXl5fxww8/GIcPH3YsqampVt1CmVDU5/xP6qV1cUV9xvv27TMCAwONESNGGFu3bjXi4+ON8PBw45VXXrHqFpxeUZ/x888/bwQGBhrffPONsWvXLuO3334zatWqZQwePNiqW3B6qampxtq1a421a9cagPHWW28Za9euNfbu3WsYhmE88cQTxtChQx3H53dLf+yxx4zNmzcbH374obqll1Xvv/++ERMTY3h5eRmtW7c2li1b5tjXqVMnY9iwYQWO/+6774y6desaXl5eRqNGjYxffvmllCMue4ryjKtVq2YA5yzPP/986QdexhT1z/LZlPBcmqI+4yVLlhht2rQxvL29jZo1axqvvvqqkZubW8pRly1FecY5OTnGCy+8YNSqVcvw8fExqlatajz44IPGiRMnSj/wMmL+/PmF/hub/1yHDRtmdOrU6ZzPNGvWzPDy8jJq1qxpjB8/vtTjthmGyuxERETEtakNj4iIiLg8JTwiIiLi8pTwiIiIiMtTwiMiIiIuTwmPiIiIuDwlPCIiIuLylPCIiIiIy1PCIyIiIi5PCY+IiIi4PCU8IlKudO7cmUcfffSKzmEYBvfeey+hoaHYbDbWrVtXLLGJSMlRwiMilrnzzjt55plnrA6jyGbPns2ECROIj4/n8OHDNG7c2OqQROQiPKwOQETKp7y8POLj4/nll1+sDqXIdu7cSVRUFO3atTvvMdnZ2Xh5eZViVCJyISrhEZFL8s033+Dr68vhw4cd2+68806uuuoqkpOTi3y+JUuW4OnpSatWrQrd37lzZx566CEeffRRKlSoQEREBJ9++inp6enceeedBAYGUrt2bWbNmlXgc1lZWTz88MOEh4fj4+NDhw4dWLly5XnjsNvtjBkzhho1auDr60vTpk354Ycfznv8HXfcwUMPPcS+ffuw2WxUr17dEe+IESN49NFHqVixIj169ADM0qAOHToQEhJCWFgYffv2ZefOnQWuP3bsWGrXro23tzcxMTG8+uqrl/oYReQSKeERkUty8803U7duXV577TUAnn/+eX7//XdmzZpFcHBwkc83ffp0+vXrh81mO+8xEydOpGLFiqxYsYKHHnqIBx54gH/961+0a9eONWvW0L17d4YOHUpGRobjM6NHj+bHH39k4sSJrFmzhtq1a9OjRw+SkpIKvcaYMWP48ssvGTduHBs3bmTkyJHcdtttLFy4sNDj3333XV566SWqVKnC4cOHCyRTEydOxMvLiz///JNx48YBkJ6ezqhRo1i1ahVz587Fzc2NgQMHYrfbAXjyySd5/fXXefbZZ9m0aROTJ08mIiKiyM9TRC7CEBG5RDNmzDC8vb2NV155xahQoYKxYcMGx74BAwYYISEhxqBBgy7pXHXq1DHi4+PPu79Tp05Ghw4dHOu5ubmGv7+/MXToUMe2w4cPG4CxdOlSwzAMIy0tzfD09DQmTZrkOCY7O9uIjo42xo4d6zjvI488YhiGYWRmZhp+fn7GkiVLClx7+PDhxi233HLe2N5++22jWrVq58TbvHnzC9+0YRhHjx41AGP9+vVGSkqK4e3tbXz66acX/ZyIXBm14RGRS9a3b18aNmzISy+9xG+//UajRo0c+x555BHuuusuJk6ceNHzbN68mUOHDtG1a9cLHnfVVVc53ru7uxMWFkaTJk0c2/JLQo4cOQKYbWtycnJo37694xhPT09at27N5s2bzzn/jh07yMjI4LrrriuwPTs7m+bNm1/0Pv6pRYsW52zbvn07zz33HMuXL+fYsWOOkp19+/aRkZFBVlbWRZ+DiFw5JTwicslmz57Nli1byMvLO6fapXPnzixYsOCSzjN9+nSuu+46fHx8Lnicp6dngXWbzVZgW351WH4SUVRpaWkA/PLLL1SuXLnAPm9v7yKfz9/f/5xt/fr1o1q1anz66adER0djt9tp3Lgx2dnZ+Pr6XlbcIlJ0asMjIpdkzZo1DB48mM8//5yuXbvy7LPPXva5fv75Z/r371+M0Zlq1arlaEOTLycnh5UrV9KwYcNzjm/YsCHe3t7s27eP2rVrF1iqVq16xfEcP36crVu38swzz9C1a1caNGjAiRMnHPvr1KmDr68vc+fOveJriciFqYRHRC5qz5499OnTh6eeeopbbrmFmjVrEhsby5o1a7j66quLdK4jR46watUqpk+fXuxx+vv788ADD/DYY48RGhpKTEwMY8eOJSMjg+HDh59zfGBgIP/5z38YOXIkdrudDh06kJyczJ9//klQUBDDhg27ongqVKhAWFgYn3zyCVFRUezbt48nnnjCsd/Hx4fHH3+c0aNH4+XlRfv27Tl69CgbN250xPvBBx8wdepUJUUiV0gJj4hcUFJSEj179qR///6OH+s2bdrQq1cvnnrqKWbPnl2k882YMYPWrVtTsWLFkgiX119/HbvdztChQ0lNTaVly5b8+uuvVKhQodDjX375ZSpVqsSYMWPYtWsXISEhXH311Tz11FNXHIubmxtTpkzh4YcfpnHjxtSrV4/33nuPzp07O4559tln8fDw4LnnnuPQoUNERUVx//33O/YfO3asQDd2Ebk8NsMwDKuDEBHXsGDBAj744IMLjmNz/fXX06FDB0aPHl2KkYlIeacSHhEpFt26deOvv/4iPT2dKlWq8P333xMbG3vOcR06dOCWW26xIEIRKc9UwiMiIiIuT720RERExOUp4RERERGXp4RHREREXJ4SHhEREXF5SnhERETE5SnhEREREZenhEdERERcnhIeERERcXlKeERERMTlKeERERERl6eER0RERFze/wOIRLFjFD0llwAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import teqp, numpy as np, matplotlib.pyplot as plt, pandas\n", "import CoolProp.CoolProp as CP \n", "\n", "names = ['Nitrogen', 'Ethane']\n", "model = teqp.build_multifluid_model(names, teqp.get_datapath())\n", "pures = [teqp.build_multifluid_model([name], teqp.get_datapath()) for name in names]\n", "p = 29e5 # Pa\n", "\n", "# Trace from both pure fluid endpoints\n", "traces = []\n", "for ipure in [1,0]:\n", " # Init at the pure fluid endpoint\n", " anc = pures[ipure].build_ancillaries()\n", " rhoLpure, rhoVpure = [CP.PropsSI('Dmolar','P',p,'Q',Q,names[ipure]) for Q in [0,1]]\n", " T = CP.PropsSI('T','P',p,'Q',0,names[ipure])\n", "\n", " rhovecL = np.array([0.0, 0.0])\n", " rhovecV = np.array([0.0, 0.0])\n", " rhovecL[ipure] = rhoLpure\n", " rhovecV[ipure] = rhoVpure\n", " j = model.trace_VLE_isobar_binary(p, T, rhovecL, rhovecV)\n", " df = pandas.DataFrame(j)\n", " plt.plot(df['xL_0 / mole frac.'], df['T / K'])\n", " plt.plot(df['xV_0 / mole frac.'], df['T / K'])\n", " traces.append(j)\n", " \n", "# Do the VLLE solving\n", "for soln in model.find_VLLE_p_binary(traces):\n", " T = soln['polished'][-1]\n", " print('rhovec / mol/m^3 | T / K')\n", " for rhovec in soln['polished'][0:3]:\n", " rhovec = np.array(rhovec)\n", " rhotot = sum(rhovec)\n", " x = rhovec/rhotot\n", " p = rhotot*model.get_R(x)*T*(1+model.get_Ar01(T, rhotot, x))\n", " plt.plot(x[0], T, 'X')\n", " print(rhovec, T)\n", " \n", " # And also carry out the LLE trace for the two liquid phases\n", " opt = teqp.PVLEOptions()\n", " opt.integration_order = 5\n", " opt.init_dt = 1e-10\n", " # Or could be 1 depending on the initial integration direction, do not know the direction \n", " # a priori because not starting at a pure fluid endpoint\n", " for init_dt in [-1]: \n", " opt.init_c = init_dt \n", " rhovecV, rhovecL1, rhovecL2, T = soln['polished']\n", " j = model.trace_VLE_isobar_binary(p, T, np.array(rhovecL1), np.array(rhovecL2), opt)\n", " df = pandas.DataFrame(j)\n", " plt.plot(df['xL_0 / mole frac.'], df['T / K'], 'k')\n", " plt.plot(df['xV_0 / mole frac.'], df['T / K'], 'k')\n", "\n", "# Plotting niceties\n", "plt.ylim(top=280, bottom=100)\n", "plt.gca().set(xlabel='$x_1$ / mole frac.', ylabel='$T$ / K', title='nitrogen(1) + ethane(2)')\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.14" } }, "nbformat": 4, "nbformat_minor": 5 }