{ "cells": [ { "cell_type": "markdown", "id": "56b72ffe", "metadata": {}, "source": [ "# Quantum PR\n", "\n", "The quantum-corrected Peng-Robinson model of Aasen *et al.* (https://doi.org/10.1063/1.5111364) can be used to account for quantum effects by empirical fits to the Feynman-Hibbs corrections.\n", "\n", "The conventional Peng-Robinson approach is used, with an adjusted covolume $b_i$ given by\n", "\n", "$$\n", "b_i = b_{i, PR}\\beta_i(T)\n", "$$\n", "with \n", "$$\n", "\\beta_i(T) = \\left(\\frac{1+A_i/(T+B_i)}{1+A_i/(T_{ci} + B_i)}\\right)^3\n", "$$\n", "and Twu alpha functions are used to correct the attractive part." ] }, { "cell_type": "code", "execution_count": 1, "id": "b64ec2f0", "metadata": { "execution": { "iopub.execute_input": "2025-10-15T23:11:42.158390Z", "iopub.status.busy": "2025-10-15T23:11:42.158274Z", "iopub.status.idle": "2025-10-15T23:11:43.296668Z", "shell.execute_reply": "2025-10-15T23:11:43.296175Z" } }, "outputs": [ { "data": { "text/plain": [ "'0.23.1'" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import numpy as np, matplotlib.pyplot as plt, pandas\n", "import CoolProp.CoolProp as CP\n", "\n", "import teqp\n", "teqp.__version__" ] }, { "cell_type": "code", "execution_count": 2, "id": "a7069986", "metadata": { "execution": { "iopub.execute_input": "2025-10-15T23:11:43.298427Z", "iopub.status.busy": "2025-10-15T23:11:43.298199Z", "iopub.status.idle": "2025-10-15T23:11:44.728541Z", "shell.execute_reply": "2025-10-15T23:11:44.728101Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAj8AAAG5CAYAAABhrVVvAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjcsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvTLEjVAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAjpdJREFUeJzs3Xd8VFX6x/HP9EnvpJBC6L1IbwKKIKiIZW2ra3ddsWJby9pWZVfUtfwsa9e1N1QsYKMpSO8ltASSkN4zmT7398dNhgQCJCHJTDLPe3dec+e2PGFi5ptzzz1HoyiKghBCCCFEgND6ugAhhBBCiPYk4UcIIYQQAUXCjxBCCCECioQfIYQQQgQUCT9CCCGECCgSfoQQQggRUCT8CCGEECKgSPgRQgghRECR8COEEEKIgCLhRwghhBABxS/CzyuvvMLgwYMJDw8nPDycsWPH8sMPP3i322w25syZQ0xMDKGhoVxwwQUUFBT4sGIhhBBCdFQaf5jba+HCheh0Onr16oWiKLz77rvMnz+fjRs3MmDAAP72t7/x3Xff8c477xAREcHNN9+MVqvl999/93XpQgghhOhg/CL8NCY6Opr58+dz4YUXEhcXx4cffsiFF14IwK5du+jXrx+rVq1izJgxPq5UCCGEEB2J3tcFHMntdvPZZ59hsVgYO3Ys69evx+l0MnXqVO8+ffv2JTU19bjhx263Y7fbva89Hg+lpaXExMSg0Wja/PsQQgghxMlTFIWqqiqSkpLQalunt47fhJ+tW7cyduxYbDYboaGhLFiwgP79+7Np0yaMRiORkZEN9o+Pjyc/P/+Y55s3bx6PPvpoG1cthBBCiPaQnZ1NcnJyq5zLb8JPnz592LRpExUVFXz++edceeWVLFu2rMXnu++++5g7d673dUVFBampqWRnZxMeHt4aJQshhBCijVVWVpKSkkJYWFirndNvwo/RaKRnz54ADB8+nLVr1/L8889z8cUX43A4KC8vb9D6U1BQQEJCwjHPZzKZMJlMR62vu6NMCCGEEB1Ha3ZZ8Ytb3Rvj8Xiw2+0MHz4cg8HAL7/84t2WkZHBwYMHGTt2rA8rFEIIIURH5BctP/fddx8zZswgNTWVqqoqPvzwQ5YuXcrixYuJiIjg2muvZe7cuURHRxMeHs4tt9zC2LFj5U4vIYQQQjSbX4SfwsJC/vKXv5CXl0dERASDBw9m8eLFnHHGGQD85z//QavVcsEFF2C325k+fTovv/yyj6sWQgghREfkt+P8tLbKykoiIiKoqKiQPj9CCCFEB9EWn99+2+dHCCGEEKItSPgRQgghRECR8COEEEKIgCLhRwghhBABRcKPEEIIIQKKhB8hhBBCBBQJP0IIIYQIKBJ+hBBCCBFQJPwIIYQQIqBI+BFCCCFEQJHwI4QQQoiAIuFHCCGEEAFFwo8QQgghAoqEHyGEEEIEFAk/QgghhAgoEn6EEEIIEVAk/AghhBAioEj4EUIIIURAkfAjhBBCiIAi4UcIIYQQAUXCjxBCCCECioQfIYQQQgQUCT9CCCGECCgSfoQQQggRUCT8CCGEECKgSPgRQgghRECR8COEEEKIgCLhRwghhBABRcKPEEIIIQKKhB8hhBBCBBQJP0IIIYQIKBJ+hBBCCBFQJPwIIYQQIqBI+BFCCCFEQJHwI4QQQoiAIuFHCCGEEAFFwo8QQgghAoqEHyGEEEIEFL2vCxDN9+yls1A8HjRaHXM/+rreay1zP/rG1+UJIYQQfk1afjqYuqADoHjcPHPx2fVee3j20lm+LE8IIYTwexJ+Opi6oNPS7UIIIUSgk/DTwWi0x3/LNFpdO1UihBBCdEwSfjqQvL0ZdBtyynH3mfvR1+1UjRBCCNExSYfnDsBWXc2y999k25KfTrjvs5eeKwFICCGEOA4JP34ua8tGFr30LJbysibtr3jcbVyREEII0bHJZS8/5fG4WfnZB3zx5ENYysuITkrmkkefakKfH3lLhRBCiOORT0o/ZK+xsOBfj7Lq849AURh8+plc/u/n6dq3P3M/+sYbcDRaHUOmneU9Tsb5EUIIIU5MLnv5mcriQhb861GKsw+gN5k44/qb6T9xSoN96gKO2+XktZuuBuDcu/9BzxGj271eIYQQoqOR8ONHCvbvZcFTj2EpKyUkKprz7n2Y+PQex9x/79o/qKkoJyQqmvShw9uxUiGEEKLjkvDjJw5u28JXTz2G024jNrUb5937MOGxccfcX1EUNi76FoBBp01Dp5e3UgghhGgK+cT0A1mbN/D1/MdxOR2kDhrKrLn3YQoOOe4xB7ZsJHfXdnR6PYNPP7OdKhVCCCE6Pr/o8Dxv3jxGjhxJWFgYXbp0Yfbs2WRkZDTYZ/LkyWg0mgaPG2+80UcVt579G9by1VOP4XI66H7KSM6756ETBh/F42HFh+8CMGTaWYTFxLZHqUIIIUSn4BfhZ9myZcyZM4c//viDn376CafTybRp07BYLA32u/7668nLy/M+nnrqKR9V3Dr2rF3F108/gdvloteoccy68370RuMJj8tYtYLCrH0Yg4IYfd5F7VCpEEII0Xn4xWWvRYsWNXj9zjvv0KVLF9avX8+pp57qXR8cHExCQkJ7l9cm9m9cy7f/+Rcet5s+Yycy4+Y7m9Rvx2mzseKj9wAYcc75BIdHtHWpQgghRKfiFy0/R6qoqAAgOjq6wfoPPviA2NhYBg4cyH333UdNTc0xz2G326msrGzw8Bc5u7az8Jl5eNxu+o6fxMxb7mpyh+XfP/uAyqICwmLjGH7W7LYtVAghhOiE/KLlpz6Px8Ptt9/O+PHjGThwoHf9ZZddRlpaGklJSWzZsoV7772XjIwMvvzyy0bPM2/ePB599NH2KrvJCrP289W/D/fxOfOmO9DqmjYTe/6+PWz4Tp2364zr5mA0B7VlqUIIIUSnpFEURfF1EfX97W9/44cffuC3334jOTn5mPv9+uuvnH766ezdu5cePY4eC8dut2O3272vKysrSUlJoaKigvDw8Dap/UTK8nL5+OF7qakop2vf/lxw/2MYTOYmHet2ufjg/jsoOpBJ3/GTOOvWu9u4WiGEEML3KisriYiIaNXPb79q+bn55pv59ttvWb58+XGDD8Do0epoxscKPyaTCZPJ1CZ1tkR1WSmfP/EPairKievWndn3PNTk4AOwbuGXFB3IxBwaxpQrr2/DSoUQQojOzS/Cj6Io3HLLLSxYsIClS5eSnp5+wmM2bdoEQGJiYhtXd/KcNhtfPfUYlUWFRCYkcsF9j2IOCW3y8fl7d7Pysw8AmPyX6wiOiGyjSoUQQojOzy/Cz5w5c/jwww/5+uuvCQsLIz8/H4CIiAiCgoLYt28fH374ITNnziQmJoYtW7Zwxx13cOqppzJ48GAfV398isfDDy89S8H+vZjDwrngvscIiYxq8vH2mhq+e2E+Hreb3mMm0P/U09qwWiGEEKLz84vw88orrwDqQIb1vf3221x11VUYjUZ+/vlnnnvuOSwWCykpKVxwwQU8+OCDPqi2eVZ8/B571qxEp9dz7l0PEJnQvJaqX956hfKCPMLjunDGDTej0WjaqFIhhBAiMPhF+DlRn+uUlBSWLVvWTtW0nq2//sjarz8HYPqNt5Hcd0Czjt++7Bd2rliCRqtl5i13N+tSmRBCCCEa55fj/HQG2du38PMbLwEw9sJL6TdxSrOOLz2Uyy9vqi1i4y68jK59+rV6jUIIIUQgkvDTBiqLi1hYO3pz3/GTGHvhZc063l5Tw9dPP47TbiO5/0BGnfenNqpUCCGECDwSflqZy+Fg4bNPYq2qpEt6D6bdeGuz+unUdZAuzc0mNDqGs2+7F622aYMgCiGEEOLEJPy0sl/ffpX8fXswh4Yxa+79GIzNG2to1Rcfs2/dH+j0embdeX+z7gwTQgghxIlJ+GlFW35ZxNZff0Sj0XLWbfcQ0SW+WcfvXbeaVZ9/CMDU6+aQ2LNPW5QphBBCBDQJP60kb08Gv771KgDjL7mCboOHNev4kpxsfvi/pwEYduY5DJxyRqvXKIQQQggJP63CWl3Fwv/8C7fLRc+RYxl17oXNOr6msoKvnnoMh9VKcv+BTLri2jaqVAghhBASfk6Soij8+OrzVJUUEZmQyJk33dGsDs5Ou42v/v0Y5QV5RHSJ55zb/45O7xfDLwkhhBCdkoSfk7Rp8bfsXat2UD77tnsxBQc3+ViPx833Lz5D3t4MzKFhnH/fozJvlxBCCNHGJPychILMfSz735sAnHr5NcR379ms45e99yZ7167yTn0RnXT8meyFEEIIcfLk+koLOaw1fPf8v3G7XPQYMZphZ57TrOPXf/c1G374BoAz58wlud/Atigz4L1046/H3T7nVZkoVgghAo20/LTQL2+9SlneIUJjYpl+423N6uez+4/fWPq/NwCYeNlV9B13aluVGdBOFHyauo8QQojORcJPC2Ss+o0dy39Vx/O59W6CwsKbfGzW5g1898LToCgMOWMGI2dd0IaVCiGEEOJIEn6aqbq0xDth6ejz/tSsmdpzM3by9TNP4HG76D1mAqddc2OzWoxE09mqnb4uQQghhJ+SPj/NoCgKi199Hlt1FfHdezLmgkubfGxh1n4W/OsRXHY73YYOZ+Ytd8qcXa3E7fZQmFVFzq5SCrMqKc6pprrM3qRjh05NwWZxYg4xtHGVQggh/IWEn2bY/OP3ZG3egN5gZMacO5s8Hk9ZXi5fPPkQ9hoLSX36M2vufej08mF7MiqKajiwrYTsnWXk7i7DaXO36Dybfs5mz9oCzrltKDFJoa1cpRBCCH8k4aeJSg/lsOz9twCY+OeriUlOadJxlcVFfPb4g9RUlBPXrTvn3fsQBpO5LUvttKrLbOxdX8ietQUUHqhqsM0cYiC5bxSJPSOJSwklpmsor9+xvEnntVQ4+PbFzVz60GiMQfKfhBBCdHbym74JPB43i155DpfDTtrgYQybflaTjqsqKeazx+6nqriIqMSuXHj/Y5hDpHWhOazVDvatL2TPukIO7S0HRV2v0WpI6hVJav9oUvpFE5scikbbsv5T4XFBVBZZ2bIkmxEz01uveCGEEH5Jwk8TbFr0LXm7d2EMCmLaX29Foz1xP/Gq0mI+few+77QVFz74uIze3ESKRyFndxk7Vhxi/6YiPG7Fuy2xZwS9RsTT45QuBIcbT/przXn1NHauPMSv7+1i7/oiCT9CCBEAJPycQHlBPis+fg+AU/98DeGxcSc8prq0hM8ee4Dy/DzC4+K56OF5TTou0DmsLnauzGPr0hwqiqze9XGpYfQaEU/PEV0Ii276JcM5r57WYByfxl4DdO0dBUB5QQ2KosgdeEII0clJ+DkORVH48b8v4LLbSRkwmMGnTz/hMdVlpXz6zwcoy8slPK4LFz30JOGxXdqh2o6rstjKll9z2LHykLfjssGso8+oBPpPTCIuJazF5z5yBOfGRnSu6+fjdnnwuBR0Bgk/QgjRmUn4OY6tvywme/sW9EYT02645YSXu6rLSvnssfspO5RDWGwcFz30JBFd4tup2o6n6GAV6xcdYP/GQpTaK1tRCcEMPi2FPqMTMJjaZyiAujGBDCYdOoMMfSWEEJ2dhJ9jqCwuYtn76qSlEy75C5EJicffv6iQz/75AOUFeYTFxHHRQ/OI6JLQHqV2OHn7Klj/QxYHtpV416X2j2bw6Smk9otuccfllirLtwAQ2oxLakIIITouCT+NUBSFn994CYfVSmLvvgybcfZx9y89lMvnjz9IVUlRbefmJ4iMl+BzpEN7y1mzcD+5GeUAaDTQc0Q8w89MI6ar7+6Cy9tbAUBijwif1SCEEKL9SPhpRMbK5WRuXIdOr2f6X2877kjMRQez+Lx2HJ/opGQu/MfjhEXHtmO1/q/oYBV/fL2fg9vVlh6tTkPfMQkMm55GZJdgn9amKApZW4sBSOoV6dNahBBCtA8JP0ewWapZ+p464/ro8y8+7mCGeXsz+PLJh7FZqonr1p0L739Mbmevpyzfwupv9rNvQxGgjs3Tb3wiI2Z0a9ZdW22pOLuasvwadAYt6YMltAohRCCQ8HOE3z7+H5byMqKSkhk568Jj7pe9YysL/v0YTpt6aez8vz8iAxjWslU7WfNdJtuX5eLxKKCB3iPjGXl2us9beo604/dDAHQbFCujOwshRICQ3/b15O3NYPNP3wMw9dqb0Bsan38rc9N6vnn6CVxOB6kDB3Pu3f/AaA5qz1L9ktvlYduyXNZ+l4m9xgVAt0ExjJndw6d9eo7FWuVg18o8AAZO6urjaoQQQrQXCT+1PG43P7/+MigK/SdOIXXg4Eb32/n7Mha99B88bhfdTxnJOXfch9548iMNd3RZW4r57fM9VBSqgxPGdA1l/J96ktI32seVHdvWpTm4nB7iUsPo2jvS1+UIIYRoJxJ+am1a/C2FWfswhYQw6YprG91n3cIvvZOb9hk7kRk3N31m986qqtTGik92k7lZ7TQcFG5kzKzu9B2XiLadb1lvDpvFyZYlOQAMm5YqozoLIUQACexP7lrVZaX8/un7AJx62dVHdVpWPB6Wvf8W67/7CoBTZp7L5CuubdIcX52V2+1h88/ZrP0uE5fDg1arYegZKQyf0Q2j2f9/rNZ9n4W9xkVM1xB6nCIjcAshRCDx/0+pdvDbR+/isFpJ6NmbQadNa7DN5XSy+JXn2PX7MgBOvfwaRpx9XkC3FBzaU87SDzMoy1MHB0zqFcmkS/sQnRTi48qapqLIytalaqvPuPN7+nULlRBCiNYX8OHn0O5dbF/2CwCnXf3XBq059poavnnmcQ5u24JWp2P6326n/8QpvirV5xw2F38s2MfWZbkABIUZGHdBT/qMTuhQYfD3z/fgcSuk9IsidUCMr8sRQgjRzgI6/CgeD7++/V8ABkyeSmLPPt5tlcVFfPXUYxQdyMRgDmLW3PvoNuQUX5Xqc9m7Slnyv11UldgA6D8hibHn9cAc0vgdcf5q/6YiMjcXo9VqGH9hL1+XI4QQwgcCOvxsW/YzBfv3YAwKZuKlV3rXF+zfy4KnHsNSVkpwRCTn//0R4rv39GGlvuOwulj55V62r1DHwwmLNjPlL339+i6uY3FYXSz/eDcAQ6el+uXt90IIIdpewIYfm6WaFR++C8DYCy8lJDIKgD1rV/H9i0/jstuJSU7lvHsfDtiZ2XN3l/HzOzuoLrUD6lg4Y8/r0SE6NDfmj6/2YSm3Ex4XxMiZ3XxdjhBCCB/pmJ9ireCPLz7CWllBdFIyw848G0VRWPftApZ/8DYoCt2GnMLZt9+LKbhjdOJtTW63h7ULM1m/+AAoEB5rZsoV/UjuE+Xr0lose0ept6/S5Ev7oDcee742IYQQnVtAhp+S3Gw2LvoWgClX3QBo+Om1F9n6648ADJl2FqdddQNaXeB9QJYX1vDTWzsozKoEoN+4RCZc1KvDtvaAOqbPL+/uANTWq5T+He+SnRBCiNbTcT/RTsLyD97G43bTffgoEnr25st5D3Nw22bQaJjyl+sYNmNWh7p7qTUoisKuVfms+GQ3TrsbU7CeyX/uS8/hHXsMHEVRWPpBBpYKB5HxwYy7IDD7bgkhhDgs4MJP9o5t7F+/Bo1Wy9BpM/nwwbsoO5SDwWTmrNvupsfw0b4usd057W6WfriL3asLAOjaO5LTr+rvNzOvn4xdq/LYt6EQrVbDGdf0xyCXu4QQIuAFXPj5/ZP3AOg29BS+fe4pHNYaQmNiOe+eh+jSrbuPq2t/ZfkWFr22jdJDFjRaDaNnpTNsWlqnGPivOKeKZR+pd3eNPCedLmnhPq5ICCGEPwi48FN0IJOQoCAyN64HRaFr3/6cc8d93ru9AsmedQUs+d8unHY3weFGpl8/gKRenePfwW51sei/23A7PaQOiGH49DRflySEEMJPBFz4AXC7XAAMOWMmU666Hp2+Yw3Ud7LcLg8rv9jrndiza+9Izrh2ACERJh9X1joURWHJezupKLISGmXijKv7o+kELVlCCCFaR0CGH61Oz+nX3MjgqWf6upR2V1PpYNFrW8nbWwHAKdPTGD0rHa2u80zSuumnbPZtLEKr0zD9hoGYQwMr3AohhDi+gAs/hqBgLvrHE3Tt29/XpbS74pxqvn95C1WlNoxmHVOv7k/6kDhfl9WqsrYWs3LBXgDGX9iLhPQIH1ckhBDC3wRc+Lnsn0/TtXsPX5fR7vZvKuKnt3fgsruJiAti5k2DiU7sXAM4luZZ+OnN7aCoc48NmtzV1yUJIYTwQwEXfsJjO1dLx4koisL6Hw6w+pv9ACT3jWL69QM73ISkJ2KzOPnu5S04bG6SekVy6iW9A26sJiGEEE0TcOEnkLidHn793052r1HH7xk0OZnxf+qJrhP17wG1A/ei17ZRWWQlLMbMmTcMRKfvXN+jEEKI1iPhp5OyW1388OoWcjPK0Wo1TLykNwNP7XyXgRRF4df/7SQ3owyDScfMvw0mKMzo67KEEEL4MQk/nVB1mZ1v/28TJbkWDCYdZ/51IKn9Y3xdVpv44+v97F5dgFar4cwbBhKbHOrrkoQQQvg5v7g2MG/ePEaOHElYWBhdunRh9uzZZGRkNNjHZrMxZ84cYmJiCA0N5YILLqCgoMBHFfuvkkPVfPHUOkpyLQSHGznvrlM6bfDZtiyHDYsOADD58r6kDuic36cQQojW5RfhZ9myZcyZM4c//viDn376CafTybRp07BYLN597rjjDhYuXMhnn33GsmXLOHToEOeff74Pq/Y/ubvL+HL+BqrL7EQlBHPBPcOJSwnzdVltYv+mIpZ/rE5dMeqcdPqNS/RxRUIIIToKjaIoiq+LOFJRURFdunRh2bJlnHrqqVRUVBAXF8eHH37IhRdeCMCuXbvo168fq1atYsyYMSc8Z2VlJREREVRUVBAe3vnmeMraUsyi17bhdnlI7BHBzJsGd7o7uurkZpSx8P8243Z66D8hicl/7iN3dgkhRCfVFp/fftHyc6SKCnX04ejoaADWr1+P0+lk6tSp3n369u1Lamoqq1atavQcdrudysrKBo/Oas+6An54dStul4f0IbHMum1opw0+BZmVfPfyFtxOD90GxzLpUrmlXQghRPP4XfjxeDzcfvvtjB8/noEDBwKQn5+P0WgkMjKywb7x8fHk5+c3ep558+YRERHhfaSkpLR16T6xc+UhfnpzOx6PQu9R8Uy/YSB6o87XZbWJktxqFr64CafdXTte0YBONS2HEEKI9uF3nxxz5sxh27ZtfPzxxyd1nvvuu4+KigrvIzs7u5Uq9B9blmTz63u7UBToPzGJqVf173Rj+NQpL6jh6+c3Ya9xkdA9nBk3DkJv6JwhTwghRNvyq1vdb775Zr799luWL19OcnKyd31CQgIOh4Py8vIGrT8FBQUkJCQ0ei6TyYTJ1DlmKW/MhsUHWLVgHwBDp6Yw7oKenfbyT2WJla+f34i10kFMcihnzRmC0exXP7pCCCE6EL9oJlAUhZtvvpkFCxbw66+/kp6e3mD78OHDMRgM/PLLL951GRkZHDx4kLFjx7Z3uT5XP/iMPDu9cwefYitfPbOR6lI7kfHBzLq18/ZnEkII0T784s/nOXPm8OGHH/L1118TFhbm7ccTERFBUFAQERERXHvttcydO5fo6GjCw8O55ZZbGDt2bJPu9OpMNv500Bt8Rs/qzoiZ3XxbUBuqLLby1bMbqSq1EdEliHNvH0ZwuIzeLIQQ4uT4xa3ux2q1ePvtt7nqqqsAdZDDO++8k48++gi73c706dN5+eWXj3nZ60id4Vb3zb9k89tnewB1bJuRZ6Wf4IiOq7LYyoJnN3hbfGbfMYyQyM57GVMIIUTj2uLz2y/CT3vo6OFny5JsVnyiBp8RM7sxelZ3H1fUdiqKrHz1Hwk+Qggh2ubz2y8ue4nj27Y81xt8hp+ZxqhzOm+LT3lBDV8/t5HqstrgM3cYIRESfIQQQrQeCT9+bs+6ApZ9pM5zNmxaKqPP7d5pOzcX51TzzQubsFY6JPgIIYRoMxJ+/NiB7SX8/PYOUGDAqV0Ze16PTht88vdX8O3/bcZe4yI2JZRzbhkqnZuFEEK0CQk/fipvXwWLXt2Kx63Qa0QXTr2k807jkL2rlO9f2YrL7iahewRn3zwYU7Dczi6EEKJtSPjxQ8U51Xz30mZcTg+pA6I5/ar+aLWdM/js31TE4je24XEppPSLYsaNgzGYZORmIYQQbUfCj5+pLLay8IW6aRwiOPOvg9Dp/WIsyla3c2UeS97fheJR6D4sjmnXDEBn6JzfqxBCCP8h4ceP2CxOvv2/zdRUOojpGsJZcwZj6ISTlCqKwrrvs1izMBOAvmMSmHJFX5mkVAghRLuQ8OMn3E4PP7y6lbL8GkIiTZx985BOOY2Dx+1h2Ue72fHbIQBOmZ7GmNmd9w42IYQQ/kfCjx9QFIVf/7eTQ3vKMZh1nH3zEEKjzL4uq9U57W5+fGMbWVtL0Ghg4sW9GTQ5+cQHCiGEEK1Iwo8fWP3NfnavKUCj1XDmDQOJTQ71dUmtzlrl4NuXtlCYVYnOoGXatQPoPjTO12UJIYQIQBJ+fGznykOs/+EAAJP/3IfU/jE+rqj1leVb+PalLVQWWTGF6DnrpiEk9ojwdVlCCCEClIQfH8rbW87SD9TRm0fM7Eb/8Uk+rqj1Ze8oZdHr23BYXYTFmDnnliFEJYT4uiwhhBABTMKPj1SV2vjhv+oghj2GxTHq7M43X9e2ZTks/2QPikchsUcEM24cRFCYjNoshBDCtyT8+IDD5uK7l7dgrXISmxLK6Vf1R9OJBjH0uD389vleti7JAaDP6ASmXN5XxvARQgjhFyT8tDPFo/DLuzspyakmKMzAzL91rhGN7VYXP76+jYM7SgEYM7s7p0xPk1vZhRBC+A0JP+1s3Q9Z7N9YhFanYcZfBxEW3XluaS/Lt3jHKtIbtEy9pj89hnXxdVlCCCFEAxJ+2tGBbSWs+VYd1XjSZX1I7Bnp24JaUebmIn5+ewcOm5uQSBMz/zaILmnhvi5LCCGEOIqEn3ZSUWTlp7e2gwIDTu3aae7sUjwKa7/PYm1tqEvsGcH06wcSEmHycWVCCCFE4yT8tAOXw82i17Zir3ERnx7OxD/18nVJrcJudfHz2zvI2lIMwKDJyYy/sGennYhVCCFE5yDhp40pisKyDzMozlY7OJ95w8BOcddTaZ7av6e8oAadXsuky/rQb1yir8sSQgghTkjCTxvb+Xseu/7IR6OBadcO6BRzdu1em8/S9zNw2t2ERpmYcaP07xFCCNFxSPhpQ8U51Sz/ZDcAY2b3ILlvtI8rOjkup5vfPtvL9uW5AHTtHcm06wYSHC4DFwohhOg4JPy0EYfNxeLXt+F2ekgdEMOwM1J9XdJJqSiysvj1bRQdrALU6ThGntUNra7jX8ITQggRWCT8tAFFUVj+0W7KC2oIiTQx9ep+HXoE5/0bi/jlvZ04rC7MIQamXtOftAGdbwJWIYQQgUHCTxvYtSqPjNX5aLQapl07gKDQjnlZyO3ysGrBPjb/kg1AQvdwpl03sFMNzCiEECLwSPhpZeUFNSz/WO3nM+qcdJJ6Rfq2oBYqy7fw01s7vJe5hk5NYcx5PdDJZS4hhBAdnISfVuR2efjxze24HB669oli+PQ0X5fUbIqisHNlHis+2Y3L4cEUouf0v/QjfUicr0sTQgghWoWEn1a0ZuF+ig5WYQrWM7UDztRur3Gy9IMM9q4vBKBrn0imXjWA0CgZrVkIIUTnIeGnleRklLHhx4MATLmib4cLDHl7y/nprR1UldrQajWMmpXOsGlpaDtYgBNCCCFORMJPK7DXOPn57R2gQP/xiR1qJnO328O677JYv+gAikchPNbMtGsHEp8ugxYKIYTonCT8tILfPt2DpdxORFwQ4zvQvF2lhyz8/M7hTs29R8cz6ZI+GIPkx0IIIUTnJZ9yJylzSzG7/sgHDZx+VX+MZv//J1U8Cpt/zeaPr/bjdqmdmidd2odeI+J9XZoQQgjR5vz/k9qP2aqdLH1/FwBDp6aS2CPCxxWdWGWJlV/f3Unu7nIAUgfEcNpf+hIS0bH6KAkhhBAtJeHnJCz/ZDc1lQ6iEoIZfU66r8s5rrpb2H//bA8Omxu9SceEC3vSf0ISGo10ahZCCBE4WjRindPp5PTTT2fPnj2tXU+HsW9jIXvWFqDRwOlX9kdv1Pm6pGOqKrXx7YubWfK/XThsbhK6R3DJgyMZMLGrBB8hhBABp0UtPwaDgS1btrR2LR2GtcrBsg8zABg2Pc1v74xSFIXtKw6x8su9OG1udHoto2alM3RqqtzCLoQQImC1+LLX5Zdfzptvvsm//vWv1qynQ1j20W6sVU6ik0IYdZZ/Xu6qKLKy5P1d5GaUAZDQPYLT/tKXqIQQH1cmRODZO/UMnDk5J9yv366d7VCNEKLF4cflcvHWW2/x888/M3z4cEJCGn6oPvvssyddnD/as66AfRsK0Wg1TL2qPzqDf8115fEobF2awx9f7cPl8KA3aBkzuweDpiRLa48QPtDU4AOws28/0Onot31bG1clRGBrcfjZtm0bp5xyCgC7d+9usK2z9iOxWZys+ET9XofPSCMuNczHFTVUnFPNkvd3UZhVCUBSr0hO+0tfIuKCfVyZEJ1b0csvU/zCi61zMrebnQMGSgASog21OPwsWbKkNevoEFZ+uRdrlZOoxBBGzOjm63K8XA43a7/LYtNPB/F4FIxmHWPP78mACUkdbn4xITqaVg0+ddzu1j2fEKIBudW9iQ7tKWPn73kATP5zH3R6/7jclb2zlKUfZlBZZAWgx7A4Jl7cm5BIGbdHiPZQ/OL/tf5Jdf5796gQncFJh58dO3Zw8OBBHA5Hg/WzZs062VP7DbfTw5L31bu7+k9MIqlnpG8LAmoqHaz8Yi8Zq/MBCIk0ceolvek+NM7HlQkRGNzVFhyZmYSdcQZVP/7YqueWS15CtK0Wh5/9+/dz3nnnsXXrVjQaDYqiAIf7+7g7UbPt+sUHKC+oISjcyNjZPXxai8ejsH15Ln98vR+H1QUaGDQpmTHndpc5uYRoZYrLhTM3F3tmJo7MLByZmTiy1GdXUVGbfd2dffvJnV9CtKEWf1redtttpKen88svv5Cens6aNWsoKSnhzjvv5Omnn27NGn2qLN/C+kVZAEy8qBfmEIPPainIqmTZhxneiUjjUsOYdGkfvx1nSIiOQFEU3GVlDYKNvS7oZGeD03nMY3WxsWiNRpyHDrVjxUKIk9Xi8LNq1Sp+/fVXYmNj0Wq1aLVaJkyYwLx587j11lvZuHFja9bpE4pHYekHGXhcCmkDY+g5vItP6rBZnKz+ej/bVuSCAsYgPWPO7c6AU7vK7etCNJHHbsdx4MBRLTj2rCw8FRXHPE5jNmNMS8OYno4xvRum9HSM3bph7NYNXXg4O/v1b8fvQgjRGlocftxuN2Fh6q3esbGxHDp0iD59+pCWlkZGRkarFehLGavzObSnHL1Ry6mX9G73W/g9bg87fjvE6m8ysVnUvz57j45n3Pk9ZSJSIRqhKAqugoLa1puGl6qcublQe3n+KBoNhsRENeB069Yg6OgTEtBoj32DQ/R111L6+hut+n3IJS8h2laLw8/AgQPZvHkz6enpjB49mqeeegqj0chrr71G9+7dW7NGn7BbXaxcsA+AkWelEx4b1K5fPyejjN8+3U1JrgWA6KQQTr24N137RLVrHUL4o7rOxo6s2oCTVXupKisLxWo95nHasDCM6emY0rs1DDppaWjN5mbVYNu5k9L3/kfld9+d7LdzlL1nnEHPn35q9fMKIVQtDj8PPvggFov6wfzYY49x9tlnM3HiRGJiYvjkk09arUBfWfttJtZKB5HxwQw5PaXdvm5lsZXfv9jL/o1qZ0pTsJ7Rs7ozYGISWp1/3F4vRHtocWdjvR5jcnJt603DS1W6mJiTasFVPB6qf/2VknfewbpuvXe9edAgov9yBWHTp5MxeEiLz1/Hmd20EaGFEC3T4vAzffp073LPnj3ZtWsXpaWlREVFdfgRnksPWdi6RP3lM/GiXu0ypo/D5mLDogNs+jkbt8uDRqth4KldGXV2OuZQ33WyFqIteTsb14aa5nY2NnXrdtSlKmNyMhpD6/43o7hcVP6wiJLX/ot9z151pV5P+PTpRF9xOUFDh3r3NSQnN3k6i2MxJCef1PFCiONrlXuj625zj46Obo3T+ZSiKKz4dDcej0L6kFhSB8S07dfzKOxek8/KBfuoqVDHSkruG8WEP/Uipmtom35tIdqLp6ZG7WyclaW25GRl4chSX3sqK4953Ik6G7c1xeGg/OuvKXn9DZwHDwKgDQ0l6rLLiPrzZRji4486pufPPx2ez0una/pozbX7GpKT6fmzXPISoi2dVPh58803+c9//sOePXsA6NWrF7fffjvXXXddqxTnC/s2FJGzqwydXsv4C3u16dcqyKxkxae7KchUf/mHx5oZf2Ev0ofEdvjWMxF4FKdTvUyVlVUbbrJq++Nk4SooOPaBdZ2Nva03Te9s3FYUj4fKb7+l6PkX1I7SgC4ykuirriTqsstOGLzqh5emTGwqgUeI9tXi8PPQQw/x7LPPcssttzB27FhAvf39jjvu4ODBgzz22GNNPtfy5cuZP38+69evJy8vjwULFjB79mzv9quuuop33323wTHTp09n0aJFLS2/UU6Hm9+/UIPcsOmpRMS1TSfniiIrq7/ex551hQAYTDpGzOzGkNNS/G6WeCHqUxQFV2HR4XBT/5GdDS7XMY/VRUV5W23US1Vp6nJqarM7G7cVRVGwrFhB4bP/wb5rFwD6uDiir72GqIsuQhvc/EmCJdQI4X9aHH5eeeUVXn/9dS699FLvulmzZjF48GBuueWWZoUfi8XCkCFDuOaaazj//PMb3efMM8/k7bff9r42mVr/Vu8Niw5QXWonNNrEKdPTWv381moH677PYtuyXDxuBTTQd3QCY87rIbeuC7+hKArukhIcBw/iOHAQx8EDOA8c9IYcT03NMY/VmM2HA05tuDHV3k2li4xsv2+iBWw7d1Iw71/UrFkDqHeGxVx/PdFXXI42qH3v9hRCtK0Whx+n08mIESOOWj98+HBcx/nrrzEzZsxgxowZx93HZDKRkJDQrPM2R0WRlY0/qtf0J1zYC4Ox9SYWdDrcbP4lm42LD+Cwqdf/U/pHM/a8HsSlhLXa1xGiqRSPB1dh4eFwczBbDTsHD+I8cOC4AQetFkNysnppyht01Ic+Pt4nl6lOhruykqIXXqTsww/B40FjNBJ1+eXEXH8d+igZWkKIzqjF4eeKK67glVde4dlnn22w/rXXXuPPf/7zSRd2pKVLl9KlSxeioqI47bTTePzxx4mJOXZnZLvdjt1u976uPE6nSoBVX+7F7fKQ3DeK7sNaZ3JQj0dh16o81nyzH0ttZ+bYlFDGnd+TlH4dv3O48G+K240zLw9nbahRg85BnAcP4DiYjVLvv4+j1PbDMaSlYkxNw5iacviuquRkNEZj+30jbURRFCq+/prC+U/jLikBIHzmDLrcdReGpCQfVyeEaEvNCj9z5871Lms0Gt544w1+/PFHxowZA8Dq1as5ePAgf/nLX1q1yDPPPJPzzz+f9PR09u3bx/3338+MGTNYtWoVOl3jLTTz5s3j0UcfbdL58/aWs29jERoNTPhTr5PubKwoCge2lrDqq32UHlLHQgqLNjP63O70HhmPRqakEK2krpNx/XDjvUyVm3vcW8XR6TAkd60NN6kY01IxpKRgTEvDkJyMthMEnGNxHDxI3gMPUrN2LQDG7t1J+MeDhNT2XxRCdG4aRTnWeO9HmzJlStNOqtHw66+/tqwgjeaoDs9H2r9/Pz169ODnn3/m9NNPb3Sfxlp+UlJSqKioILzenRqKovDFU+spyKyk/4Qkplzet0V118ndXcbqb/aTt1edK8gUrGfEzG4MmpQsnZlFi3jsdpzZ2Q374NRepnIeOnTcW6k1BgOG1FQ13KSm1C6nqUEnMbHVx8Pxd4rHQ9nHH1M4/2kUqxVNUBBxc24i+i9/6RStWUJ0RpWVlURERBz1+X0ymtXys2TJklb5oiere/fuxMbGsnfv3mOGH5PJ1KRO0XvXF1KQWYnepGPUOektrqkgs5LV3+wje2cZADqDlsFTkjlleppPZ4IXHYPHYsFRG3CcDVpxDuLKzz/2nFTUdjKua7mpF26MqalqH5xjtI4GGkdOLnkPPkjNH38AEDxqFIlPPoFRBhQUIuC0yiCH7S0nJ4eSkhISExNP6jxup4c/vlLn7zplWmqL7rgqzqlm9Tf7ydpSDIBWp6H/+CSGz+hGaJTcwSUO81itaqDJOqAO+Hcgq/b5AO6i4uMeqw0NxZiaWq8PzuGwo4+Lk3GhTqBi4bfkP/IIHosFjdlMlzvvJOrPl3W4ztlCiNbhF+GnurqavXv3el9nZmayadMmoqOjiY6O5tFHH+WCCy4gISGBffv2cc8999CzZ88GU2y0xNZlOVQW2wiJMDJ0amqzji3Lt7Dm20z21o7Vo9FAn7GJjJzZrd0nQRX+w2O317bc1AYcb9A5cPyB/lAH0WsQbuouU6WloesE08b4gsdmo+CJJyn/7DMAgk45haR5T2JMa/2hLIQQHYdfhJ9169Y16E9U17H6yiuv5JVXXmHLli28++67lJeXk5SUxLRp0/jnP/95UmP92KqdrPs+C4DR53bHYGrapYHSPAvrf8hiz9oC75WIniO6MOrsdKISQlpcj+g4FIcDR05Ooy04rrzjX6LSRkSo0zXUf3RTw44uIqIdv4vOz74/k9w77sCekQEaDbF/+xuxc26Sy4BCCP8IP5MnT+Z4/a4XL17c6l9z3fdZ2GtcxHQNpc+YE18+K8mtZt0PWexdXwi1pXYbHMvoWenEJstYPZ1N3UB/9n37se/bi2PffnWQvwMH1E7GHs8xj9WGhjYMNrXLhrQ0GTemnVR+/z15D/4DT00NupgYus5/ipBx43xdlhDCTzQ7/Dz00EOce+65DB8+vC3qaReVxVa2LlPn2hl3QQ+0x7n1vDininXfZbFvY5F3XfehcYyY2Y24VAk9HZ2iKLjy87Hv3Ydj/z7se/dh378fx969uCsqjnmcJji48RactDR00dFyicpHFI+Houeep+S11wC1U3PS0/MxdOni48qEEP6k2eEnJyeHGTNmYDQaOeecc5g1axann346xg50m+jabzPxuBWS+0aR2r/xgRILD1Sy7vssMjfXdkTVQI9hXRgxsxuxyTLbekfkrq7GnpGBbdcu7LsysGVk4Ni799ijGWs0GFJSMHXvjrFHd+80DYa0NOlk7Ifc1RYO3XMP1bXDbMRcdy1xd9whl7mEEEdpdvh566238Hg8/P777yxcuJDbb7+dvLw8zjjjDM4991zOPvtsoqP9d/Ti0jwLGavzARgzu8dR2wsyK1n7fSYHtqojvqKBXiPiGT4jjZgkCT0dgaIoOHMPYc/YhW3nLvV5VwbO7OzGD9DrMaalqSGnZw9MPXpi6tEdY3q630y4KY7PkZNDzt9uwr5nDxqjkcTH/0nErFm+LksI4aeaNcjhsezcuZOFCxfy9ddfs379ekaNGsWsWbO49NJL6dq1a2vUedLqBkn6/D8ryd9lpfvQOGbcOAhQPyxzM8rYsPiAd5wejQZ6j0pg+Iw06cjs51zFxVi3bMG6ZQu2LVuwbtuO5xjTmejj4zH17YO5bz/Mfftg6tULY2qqDHDXgVnWrCH31ttwl5ejj4sj+f9eJGjIEF+XJYRoJT4f5PBY+vXrR79+/bjnnnsoKirim2++4ZtvvgHgrrvuao0v0WoytxQTbAph9KzueDwKmZuK2LD4AIUHqgDQaDX0GR3P8DO7ERkf7ONqxZE8Viu27duxbtlaG3g24zqUd/SOBgOmHj0w9+1bG3b6YurTRzocdzJlH39C/uOPg8uFeeBAkl/6Pwzx8b4uSwjh51ql5acjqEuO86/+hiETutO1dxQbfjxARaEVAL1BS7/xSQydmiLj9PgRd1UVNevXU7NmLTVr1mDbufPo6Rw0Gkw9e2AePJigQYMJGjwIU8+e0prTiSlOJwXz/qXOxA6Ez5xJ4pNPyGVKITohv2356Ug0GsjeWcruNeqAc6ZgPYMmJzN4SjJBYfJh6WtHhZ0dO466rVzfpQtBQwZjHjSYoMGDMQ8ciC5ULk0GCnd5OTm33+GdpiLu9tuJ+esN0gFdCNFkARd+FAWsVU5Co0wMOT2F/hOSMJoD7p/BbyhuN9YtW6hetgzLb783GnaMaWkEjxqlPkYMx3CS05qIjsu+dy/ZN83BefAgmuBgus5/irBjzO8nhBDH0qJP/QMHDrBlyxbi4+MZNWpUa9fUpiLigphwbj96j4pHp5d5fXzBXVmJ5bffqF62jOrlK3CXlTXY3iDsjBopfTgEANXLlpE79048FguGrl1JfvllzH16+7osIUQH1Ozw89FHH3HVVVfhdDrRaDQMGzaMH374gbi4uLaor9VddN9IIiJlGoH2Zs/MpPrXX6leuoyaDRsa9NvRhoUROnECIaeeSsjYsRJ2RAOKolD61lsUPv0MKArBI0bQ9YXn0fvxkBpCCP/W7A7Pffv2ZezYsdx3331kZ2dz7733MmTIEN588822qrFVtEWHKXF8zkOHqPzhByq++w77jp0Nthl79iB00iRCJ00ieNgwNAaDj6oU/sxjs5H30ENUfrMQgMiLLiLhwQekM7sQAaQtPr+bHX6MRiO7d++mW7duAOzatYvhw4djsVhapaC2IuGnfbiKi6lcvJjK777HumHD4Q06HSFjxhA6ZQqhkydhTE72XZGiQ3Dm5ZFz8y3Ytm8HnY74++4j6s+XScdmIQKMX9zt5XK5CA4+PP5N37598Xg85Ofnk5CQ0CpFiY7FY7dT9eOPVCz4CssffxzusKzREDxiBOFnnUXY9Gkyxo5ospr168m59TbcJSXoIiPp+txzhIwZ7euyhBCdRIs6PL/77ruMHz+ewYMHExoail6vp+ZY8yOJTstx4ABln3xKxZdf4i4v9643DxpE+FkzCZ8xQ/rviGYr+/gT8p94ApxOTH36kPzSSxiT/WOkeCFE59Ds8DNx4kQef/xxqqqq0Gq1pKenY7PZePPNN5k6dSojRowgLExmO++sFKeTqiVLKP/4EywrV3rX6xMSiLzwQiLOORtjWpoPKxQdleJwkP/kk5R//AkAYWeeSdKTT6ANlpHWhRCtq8UjPO/Zs4f169ezYcMG76O8vBytVkuvXr3YuXPniU/SjqTPz8lxlZRQ9sGHlH/2Ga6iInWlRkPIxAlEXXIpoadORKOX8ZJEy7iKi8m57Xas69eDRqMOXHjD9dK/RwjhH31+6vTq1YtevXpxySWXeNdlZmaybt06Nm7c2CrFCd9z5uZS8tbblH/+OYrdDoAuJobICy4g8qI/ScdlcdJq1q4ld+6duIqK0IaGkvT0fMImT/Z1WUKITizg5vaSlp+mse/bR8nrb1Dx7bfgcgFgHjyYmKuuJGzqVLnVWJw0xeOh5I03KXruOfB4MPbsQfILL2Dq3t3XpQkh/IhftfyIzsm6dRslr71G1c8/q3OBACHjxhJzww0Ejx4tlyFEq3CXl3Po3r9TvWwZABHnziLh4Yelf48Qol1I+BEA2PdnUvjsM1T//It3XdgZU4m54QaCBg3yYWWis7Fu2ULO7bfjOpSHxmgk/h8PEnnhhRKshRDtRsJPgHOVlFD80kuUffKpOuWEVkvEOWcTc911mHr18nV5ohNRFIWy9z+g4KmnwOnEkJZK8nPPYe7Xz9elCSECjISfAOWxWil9911KXn8DT+3o3KFTptDlrjsx9ejh4+pEZ+OuribvwX9QtWgRAGHTppH4xOPoZFgMIYQPSPgJMIrHQ8WCryh64QVcBQUAmAcOpMvddxMyepSPqxOdkW3XLnJuuw3ngYNgMBB/991EXXG5XOYSQviMhJ8AYtu9m/yHHsa6aRMAhq5dibvjDsJnzkCj1fq2ONHpKIpC2f/ep/Dpp1EcDvRJiST/5z8EDRni69KEEAFOwk8A8NjtFL/yCiVvvAkuF9qQEGJvuomoKy5HK7esizbgKiri0P0PYFmxAoDQyZNJnPekzO8mhPALEn46Ocsfq8l/+GEcBw4AEHr66ST840EMMgmtaCNVv/5K3gMP4i4rQ2MyEf/3e4m85BK5zCWE8BsSfjopV1kZhfOfpuLLLwHQd+lC/D8eJPyMM3xcmeisPFYrBf/+t3duLlPfvnR9ej6mnj19XJkQQjQk4acTsqxcyaF7/67OwaXREHXpJcTdcYfcWSPaTM2GjeTdd5+3hTH66quJu+N2uawqhPBLEn46EcXhoPD55yl98y0AjN27k/j44wSfMszHlYnOymO3U/TCC5S+/Q54POi7dCFx3pOEjh/v69KEEOKYJPx0EvbMTA7ddTe27dsBiLzkYuLvvRdtUJCPKxOdlXXrVg79/T4c+/YBEHHuucTffx+6iAgfVyaEEMcn4aeDUxSFii8XkP/EEyg1NegiIkh84nHCpk71dWmik/I4HBS/9DIlb7wBbje62FgSH3uUsNNO83VpQgjRJBJ+OjCPxULePx6i8vvvAQgePZqkp/6NIT7ex5WJzsq6dSt5DzyIffduAMLPOov4Bx+QW9iFEB2KhJ8OypGTQ85Nc9QPIb2euFtvJebaa9DodL4uTXRC7moLRc8/T9n774OioIuKIuGRRwifPs3XpQkhRLNJ+OmALH/8Qe7td+AuL0cXG0vyCy9Ip2bRZqp+XUL+Y4/hys8HIPycc4i/7+/oo6N9XJkQQrSMhJ8OxDsr9r/+BW435oEDSf6/F2XAQtEmnIWFFDzxJFWLFwNgSE4m4ZFHCJ0gd3IJITo2CT8dhMfhIP+xx6j4/AsAwmedQ+Jjj6E1m31cmehsFI+H8k8/pfCZZ/FUVYFOR8zVVxE7Z47cPSiE6BQk/HQA7vJysm+ag3XDBtBq6XLnnURfc7VMFyBanXX7dgoefwLrxo0AmAcNIvGxRzH36+fjyoQQovVI+PFzzvx8Dl53HY69+9CGhdH12WcInTjR12WJTsZVVkbRc89T/umnoChog4OJu/12ov58mXSiF0J0OhJ+/Jh93z4OXnc9rrw89PHxpL7xOqZevXxdluhEFJeLsk8+oeiFF/FUVADq7etd7r5L+pIJITotCT9+yrppE9l/vRF3RQXG9HRS33wDQ1KSr8sSnUjN2rXkP/4E9owMAEx9+pDw4AMEjxzp48qEEKJtSfjxQ9XLl5Nz2+0oVivmwYNJ+e+rMoicaDXO/HwKn5rvHRxTGxFB3G23EnXRRWj08itBCNH5yW86P1P5/ffk3nMvuFyETJhA8vPPoQ0J8XVZohPw1NRQ8s47lLz+BorVChoNkRddRNztt0m4FkIEFAk/fqTy++/Jvetu8HgIP+sskuY9icZo9HVZooNTXC7Kv/iSov97EXdRMQBBw4YR/+ADBA0Y4OPqhBCi/Un48ROVixaRe/c94PEQcd55JD7xOBqt1tdliQ5MURSqlyyh8OlncOzfD6gDFcbdcTvhM2fKUAlCiIAl4ccPVC7+kdw77wK3m4jZs0l8/J8SfMRJsW7eTMH8+VjXrQdAFxlJ7E1/I/KSS9BKa6IQIsBJ+PGxyp9+IvfOO9Xgc+4stcVHxlURLeQ4cIDC/zxH1aJFAGhMJqKvvJKY669DFxbm4+qEEMI/SPjxoapffyX3jrngchF+zjkkPvmkBB/RIs7CQkr++xpln3wCLhdoNEScdx5xt94i4/UIIcQRJPz4SM3ateTefocafOo6N0vwEc3kKiqi5I03KPv4ExS7HYCQSafSZe6dmPv09nF1QgjhnyT8+IAtYzfZN81BcTgIPf10kv79LxlfRTSLq7iYktffoOzjj72hJ2jYMOJuu5WQMWN8XJ0QQvg3+cRtZ87cXLKvvx5PVRVBw4fT9ZmnJfiIJnOVlFDyxpuUffQRis0GQNDQocTecjMh48bJHVxCCNEEfnFL0fLlyznnnHNISkpCo9Hw1VdfNdiuKAoPPfQQiYmJBAUFMXXqVPbs2eObYk+Cq6xMnaursBBTr56kvPwSWrPZ12WJDsBVUkLBU/PZO/UMSt9+G8VmwzxkMCmvv07aRx8SOn68BB8hhGgivwg/FouFIUOG8NJLLzW6/amnnuKFF17g1VdfZfXq1YSEhDB9+nRstX/5dgQeq5WcG/+GIzMTfWIiKa+/ji4iwtdlCT/nKi2l8Omn1dDz1luHpzx5/TW6ffwxoRMnSOgRQohm8ovrLTNmzGDGjBmNblMUheeee44HH3yQc889F4D33nuP+Ph4vvrqKy655JL2LLVFFI+H3Lvvxrp5M9qICFLfeF3uwBHH5cjJofSttyn/8kvv5S3zwIHE3XIzIaeeKoFHCCFOgl+En+PJzMwkPz+fqVOnetdFREQwevRoVq1a1SHCT9HzL1D98y9oDAZSXnkZU48evi5J+Cnr9u2UvvkWlYsWgccDqKEnds5NhE6eLKFHCCFagd+Hn/z8fADi4+MbrI+Pj/dua4zdbsdeexcMQGVlZdsUeAIVC7+l5L//BSDx8X8SfMopPqlD+C9FUbCsXEnpm29iWbnKuz5kwgRirruW4NGjJfQIIUQr8vvw01Lz5s3j0Ucf9WkN1i1byHvgAQBirr+OiNrLdkKAOuFo5eLFlLz5JvYdO9WVOh3hM2cSc+01mPv29W2BQgjRSfl9+Emo7RtTUFBAYmKid31BQQFDhw495nH33Xcfc+fO9b6urKwkJSWlzeo8kjM/n+w5tWP5TJlC3O23t9vXFv7NY7VS/sWXlL79Ns7cXAA0QUFE/ulCYq68EkPXrj6uUAghOje/Dz/p6ekkJCTwyy+/eMNOZWUlq1ev5m9/+9sxjzOZTJhMpnaqsiGP3U7OnJtxFxVj6tWLpPnzZfRmgTMvj7JPPqH8409wl5cDoIuKIuqKy4m69FL0UVG+LVAIIQKEX4Sf6upq9u7d632dmZnJpk2biI6OJjU1ldtvv53HH3+cXr16kZ6ezj/+8Q+SkpKYPXu274o+joLHH8e2fTu6yEiSX3kZXWiIr0sSPqIoCjVr1lL2wQdU/fILuN0AGFJSiLnmaiJmz0YbFOTjKoUQIrD4RfhZt24dU6ZM8b6uu1x15ZVX8s4773DPPfdgsVi44YYbKC8vZ8KECSxatAizHw4QWP7lAso/+xw0GpKeeRpjcrKvSxI+4LFYqFi4kLIPPsC+53CwDx49mqjLLiNs6unSGiiEED6iURRF8XUR7aGyspKIiAgqKioIDw9vk69hy8gg66KLUex2Ym+9hbibbmqTryP8lyMri7KPPqL8ywV4qqoA0AQHE3HuLKIvuwxTr14+rlAIITqWtvj89ouWn87AXVVFzq23otjthEycSOyNN/q6JNFOFLeb6hUrKPvgQywrVnjXG9PSiPrzZUTMno2ujQK3EEKI5pPw0woURSHvgQdxHjiIPjGRpKf+jUbrFzOHiDbkKimh4quvKfv4Y5zZ2epKjYbQU08l6vLLCRk/Tn4OhBDCD0n4aQVlH31E1Y8/gsFA8vPPyV07nZjicmH5/XfKP/+CqiVLwOUCQBseTuQFFxB16SUYU1N9XKUQQojjkfBzkux79lD476cAiL/rToIGD/ZxRaItOA4epPyLL6lYsABXYaF3vXnQICL/dCER55wjd20JIUQHIeHnJHjsdnLvutvbzyfqL3/xdUmiFXmsVqp+/JHyz7+gZu1a73pdZCQR584i4vwLMPfp7cMKhRBCtISEn5NQ+Mwz2DMy0EVHkzTvSZl/qRNQFAXbtm2Uf/4Fld99h6e6Wt2g0RAyYQKRF5xP6GmnoTUafVuoEEKIFpPw00LVK1ZQ9t7/AEia9yT62FgfVyROhqukhMrvvqP88y+w797tXW9ITibi/POIPO88DPWmVxFCCNFxSfhpAVdZGYfuux+AqMsvJ3TSJB9XJFrCXW2h+pefqfj2OywrV3pHX9aYTIRNm0bkBecTPGqU3LElhBCdjISfFij45+O4i4sx9uhBl7vv8nU5ohkUh4Pq336jYuFCqpcsRbHZvNvMAwcScf55RJx1FrqICB9WKYQQoi1J+GmmykWLqfz+e9DpSPrXPLQ+mjxVNJ3i8VCzdh2V335L5Y8/4qmo8G4zdutG+DlnE3HWWRi7dfNdkUIIIdqNhJ9mcJWUkP/oowDEXH8dQYMG+bgicSyKomDbsYPKb7+j8vvvcRUUeLfpu3QhfOZMws8+G/OA/tJRXQghAoyEnyZSFIX8Rx7FXVaGqU8fmbfLT9n376fyhx+o/PY7HJmZ3vXa8HDCp08j/KyzCR45QiYVFUKIACbhp4mqFi+m6qefQK8n6V/z0Mitzn5BURTsO3dS+dNPVP30E469+7zbNCYToadNIeLsswmZOFFuTxdCCAFI+GkSd3k5+f98HIDYG27A3K+fjysKbIrHg3XTZqpqA48zJ+fwRoOBkLFjiDjrLEJPn4ouNMR3hQohhPBLEn6aoGD+fNwlJRh79CDmxr/6upyApLhc1KxdWxt4fsZVVOTdpjGbCZ04kbBpZxA6eTK6sDAfViqEEMLfSfg5Acsff1DxxZcAJP7zMbl00o48djuWlSup+ulnqn/5BXe9u7S0oaGETplC2BlTCZ04UebVEkII0WQSfo7DY7OR9/DDAERddinBp5zi44o6P1dZGZbly6lauhTL8hV4LBbvNl1UFGFTTyfsjDMIHjNGgqgQQogWkfBzHCWvvYbzwEH08fHEzZ3r63I6JUVRsGdkUL10KdVLl2HdvBkUxbtdHx9P2NSphE2bRvDwU9Do5UdWCCHEyZFPkmOwZ2ZS8vobAMQ/cD+60FAfV9R5eKxWLKv+UAPP8uW48vMbbDf17UvopEmETZmMefBgmV5CCCFEq5Lw0whFUSj45z9RnE5CTp1I2Bln+LqkDs+Rk0v1sqVUL1tGzR+rURwO7zaN2UzI2LGETppE6KRTZQJRIYQQbUrCTyOqfvgBy8pVaIxGEh58UEYAbgGPw4F1w0Ysv62getky7Hv2NthuSEoidPIkQidPJnjUKLRms48qFUIIEWgk/BzBY7FQ8K9/AxDz1xswpqb6uKKOQVEUHHv3Uv3771hWrqRm7ToUq/XwDlotQcOGqYFn0iRMvXpJqBRCCOETEn6OUPzGG7gKCzEkJxNz3XW+LsevuYqLsaxaheX3lVhWrsRVWNhguy429vDlrAnj0UVG+qZQIYQQoh4JP/U4cnIpffMtALrce4/M2H4Ej81Gzfr1WFauxPL7Suy7djXYrjGZCB4xgpBx4wiZMB5T797SuiOEEMLvSPipp/Dpp1EcDoLHjCFs6lRfl+NzisuFbft2LGvWULPqD2rWr0ex2xvsY+rXj9Dx4wgZN46g4cMlMAohhPB7En5qWdasoWrRItBqib/vvoBssVBcLmw7dlCzZg2WNWuwrluPp6amwT76+Hi1ZWfcOELGjUUfE+OjaoUQQoiWkfADKG43BU/OAyDy4osw9+nt44rah+JyYdu5i5o1qw+HnXojKgNoIyIIHjmCkFGjCBk3DmOPHgEZDIUQQnQeEn6A8i++wL5rF9rwcOJuvdXX5bQZb9hZu5aa1aupWb8eT3V1g3204eEEjxxJyKiRBI8erfbbkUEGhRBCdCIBH37c1dUUPfc8AHE334w+KsrHFbUej8WCdcsWatZvwLphPTWbNqMccRlLGxZG8MiRBI8aScioUZj69EGj0/moYiGEEKLtBXz4KX3rbdylpRi7dSPq0kt8Xc5JcRUVUbNhoxp01m/AtnMnuN0N9tGGhxN8yikEjx5N8KiRmPv2lbAjhBAioAR0+HEVFVHyzjsAxM29A43B4NuCmkHxeHBkZmLduJGa9Ruo2bAe54GDR+2nT0okePgIgoefQtApp2Dq2VMuYwkhhAhoAR1+il5+GaWmBvOQwX4/f5e7ogLrlq1YN23Cunkz1i1b8FRWNtxJo8HUu3dt0BlO8CnDMCQl+aZgIYQQwk8FbPixZ2ZS/ulnAHS5806/uoNJcbux792HdfMmrJs2Y928Gce+fUftpzGbCRo4kKDhw9XAM3QouvDwZn89j+Khwl5BjasGu9uO0+3E7rbjcDvQarTotfoGj2B9MKGGUIINwWg10ookhPBfNc4aqhxVxIfE+7oU4UcCNvwUPfc8uN2ETppEyKhRPqtDURScuYewbduGbfs2rFu3Ydu69ahbzgEMaakEDRlC0NChBA0Zgrl37yZdqnN6nORU5ZBVkUVWZRYHKg+QW51Lqa2UEmsJ5fZy3Ir7hOc5klajJcQQQpghjDBjGDFBMcQGxRITFEOMWV2ODYolPjiehJAEzHqZvFQI0X7+t+N/PLX2KQBWXLyCSHOkbwsSfiMgw49182aqFi8GjYa4uXPb9Ws7CwuxbduObdtWrNu2Ydu6DXdZ2VH7aYODMQ8eXBt2hhA0ZAj66OgmfY1iazGbCjexuWgzm4s2s6NkB3a3/YTHmXVmDDoDJp0Jo9aIUWfEo3hweVzqQ3Hh9Dixuqy4PC48iocqRxVVjiqwAEd/Gw1Em6NJDEkkMSSRhJAEEkMS6RraleSwZJLDkgkxhDTp+xNCiONRFIUrF13JxsKN3nWhxlAfViT8TcCFH0VRKHzmWQAiZs9u0wENXWVlatCpa9HZtg1XQcHROxoMmHv3xjxwIEGDBmIeOFCd9byJd2EpisKO0h0sy17G0uyl7CzdedQ+QfoguoV3Iy08jW4R3UgJSyHWHEt0UDTR5miiTFEYdE3r8K0oCna3XQ0+ziqqHdVU2CsotZVSbC2mxFaiPlvV5zxLHlaXlVJbKaW2UraXbG/0vNHmaJJDk71hKDk0mZSwFNLC04gNivWrS5NCCP9Uaitl0ieTGqxbctES9NqA+7gTxxFwPw01a9ZQs2YNGoOBuFtubrXzuquqsO3Y2aBFx5mTc/SOWi2mHj0wDxqEeeAAggYNwtS7d9PnxCo7AJs+hJ6nsz88lgV7FvB95vcU1hyeUV2Dhp5RPRkSN8T7SAtPa7X+ORqNBrPejFlvJo64E+6vKAqVjkryLfnkWfI4VH2IfEs+hyyHyK3KJac6h3J7uTccbSnectQ5gvRBpISlkBqWSmp4aoPnuOA46XskhOC7/d/x9xV/977WaXSsu3ydBB8AhwV+eQz6nwtp43xdjc9pFEVRfF1Ee6isrCQiIoIt552PfscOoq64goQH7m/2eRRFwVVYhH3XTmw7d6qBZ+dOnNnZje5vTEtrEHTMffuiDWn55R3bK+NZbMnky6gYNtT77zlIH8S4pHFMTpnMxK4TiQnqWHNuVTmqyKnKIac6h+yqbHKq1OfsqmzyLHl4FM8xjzXrzCSHJXvDUUpYCinhKaSEpZAYkii/+ITo5Iqtxdy25Da2FB3+w+mSPpfwwJgHfFiVn1AUyPgePr7s8Lq/Z4O5+TfH+Erd53dFRQXhLbippzEBF37W9OxFWEgIPX/6EX3c8VstFI8H58GDDUKObedO3CUlje6vT0okaOCgw5evBgxo0d1Xjalx1vDpxld4d+ubFOvVy2FajZZTu57K7F6zmdB1AiZd55xR3el2kludy8GqgxysPKg+Vx0kuzKb3Orc43bW1mv0JIUmqWEoNIXU8FRvSOoa1rXT/psJEQjcHjefZHzCvDXzGqz/7xn/ZVyStG6Qvw0W3weZyxuuH3EtnP2sb2pqgbYIPwH5J3HUpZceFXzc1Rbsu3dj352Bbdcu7Bm7sWdkHDWrOQBaLcb0dMz9+qmP/v0w9e3bJlNjWJwW/rfjf7y/830q7BWg15HocvEni4NZ164kft278Nb56s5THoRJd8Oyp2DJkzDlfph0T6vX1N4MOgPdIrrRLaLbUducHid51XlkV2Wrgagqm+xKdTmnKgeHx+ENS42JDYqla2hXkkKT6Bra1bucHJpMQkgCRp2xjb87IURLrMxdyfx189lbvte7TqvRsviCxSSEJPiwMj9QXQRLHocN74HiAZ0Jxs6B5JHw8aWw7k0YdGFAX/4KvJafQYMZ8P7/cOXlY8/IwLY7A/uujMb75wAaoxFTnz61Qacv5n79MPXpgzYoqE3r9SgeFu5byHMbnqPYWgxAqtPNdeXlnF1twQCQPgkylzU88Mh1Ux7oFAGoJTyKh8KaQjUYVR70BqScqhwOVh3E4jx6OIH6NGiIC44jOTTZG44SQxKJD4knPjie+JB4wgxh0hFbiHa0o2QHz294npWHVjZY/6fef+Lvo/4e2H+wuOyw+lVYNh8cVeq6AefB1EchKk19/c0taiiK6Qk3/g4G/x+CRC57nQRv+Onbj9BjfMv6Ll0w9e2DuU8fTL37YOrTG1P37mj07dtAtrVoK/PWzGNr8VYAUsJSuDluLNN+eboFTXUaeKS8lSvs+BRFocJeQW51LrnVuRyqPkROdQ6Hqg95X9vcthOeJ0gf5A1C8cH1HvUCUpQpSgKSECdpc9Fm/rv5v6zIXXHUtnkT53F297N9UJWfUBTY9S38+CCUZanrEofCmfOObt2xlsNLo6E6HybeCac/1M7FNp9c9moNLhea4GBMPXuqLTp9Dwed9p7R/dXNr/LyppeZM3QOfx3yV/5v4//x3y3/9W4P1gdzw+AbuKL/FRg/bOGkq1Okw19jNBoNkeZIIs2RDIgdcNR2RVEotZV6g1BdSMq35FNQU0BBTQEV9gqsLitZlergkcdi1BqJC47zDv4YbY4mJujwc4y59hEUQ7gxXIKSaFVH/p6pe33T0Ju4cciNvi7vuF7Z9Aovb36ZrqFdya3OPWp7v+h+zJ80n7TwNO86RVFwedQxyVyKOkaZ2+PGrbi945O5lMPr3Iobl9uFzW2jxlmDw+3A4XHgdDtxKk4cbgdOjxOn26meS3E3ONatuPEoHhRFQUFRl1Go366g0Wio+5/6/9r/aTToNDp0Gh0arQYdOvRavbpOq8OgMaDXqSPrG7QGDFp1HDaDzoBJa8KsN2MqP4j5j/9izl2PSQF9WAKGSfdgGHQxBoMJvaI0/J0SFAlnPQ2fXA6/P6+2DCUMasu30S8FXMtP4aZNxA4Y0O6tOUd6dfOrvLTpJe/rAcYYtjsOd6TuY4jilX3biDv17zD4Ynh+CKBA/AAoaHycnKN0nwx/+frwa29foAc6Zd+g9mZ1WSmsKaTAUuANREcul9ga7xx/LHqNnmhzNNFB0d5AFG2OJtwYrj5M4UQYIwg3hXvXhRnD0GmbNiaUCCxH/p4ZnTia1Xmrva/nDJ3T7ACkKAoOj0MNCXUPj8M7NU7dssNdGyA8zgb7Oz3OBs912+oCi9PtpNpZzbbibZTZjz1yqkFrIMwQhlNxegdidXvceDj2naGBTIsWrebwQ+d2oHc71WAVEt9ggNu6oUyC9EEE64MJMYQQaggl1BhKmDGMCFMEkcZIQk2hBOuDvfsE6YPa5HeRXPY6CW3xj3cyBr87GIV6//SKAvXSuUZR2JJVe/t86jg4uBKiuh1u0myqBp2gnzi8XvoGtQun20mhtZDCmkJKraWU2ErUh7XEO71IqU1dX1V3jb4FQg2h3nBUPyjVhaMQQwjB+mCCDOovsyB9EMGG4KOWZVgA/1QXOOrCRV2w8IaKeuvsbrv39f2/nXg4jxndZjQIL/WPr39uh7t2H4+zHb5j36prldGgQavRNnh95Pr6LTl1y3XnABr+nq+loKD+Xzlqua4FqbGWpMP7qi1N+GErsVajRafRYdAaMOqMmHVqkArWBxNiVKdDCjeFE2mKJMoc5Z0KKT44Xv1DzxR+1O8hCT8nwd/Cz5F/kR3p5rJy/lpeecztTVfb5+eRSGjkP8Kj9hM+43A7vEHoyHBU6aik0l6pPtc97JXUuBq5G/EkGLVGb0AK1gcTbFDDUZA+CKPOiF6r9059UvfLzaA1YNCpTfLH2mbUGr3Peq2+wV+gWo32qL9K6y4H1H3A6LQ67weOVqM9/AFR/0MBjvrgqNPo/gpHXf6om87leMv1L3u4PC71cojS8NnpcR5uyTjRsrthC0hd60f9YOPyuFr1fW5Nde+591Hvdd3PQf33v/76CnsFu0p3HfeycVOEGEIatISGGcMIN4YTYYog3BhOqDHUG/5D9CHen2uzzoxJb/J+QJt0Jv//A+DQJnhNHcFaAVxaPY7EQThj++CI6409uhv2qG5Y9HpqnDVYXVYsTgs2l40aVw01zhpqXDXq64Kt2HLWYNPqsaWNwaa4sbvt3pY77+U/j7PBz39b06DxXuoz6owYHAaWXLlE+vx0Bqcmn8rrW17H4XE03KAojLHZjgg+Ghh1PQTHwNJ641l0nwz7lzY8/sh1U+4//Fy/5edI0jfI54w6IwkhCc26TdfpcVLlqKLCXnF0QKpdrnJUUeNSfwnW/eKr+6VY46rB6rTiUtQPV4fHgcPuUIdVEH6rLoTWDxR1IaOuT0jd8pbiLZTaSo86R3JoMhf2vlC91NFIcPGuP2Jd3QdS3XJL+6gt3LewSS1Tjekf3Z9/Tvgn0eZoIowRTZ6ap1OI6wPjboGs39AU7sTgsmHI3Qi5GxvuF5YIYQlqWGq0q8MD6hWHb26GgRdCjylN+vJ1rZBWp9X7e8XitFDpqKTMVka5vZwKewUV9grvFEgWpwWL00KNswab24bdZfeGepfHdVTrmILi/QOhxlWD29r6gUvCjw98t/87HvjtgcYTtEbDH0FB/DcyXA1AIbFw1Q8QVzsHmUbbSL+d2lBzvHF+Jt0DWb8dfWs8qIFp0t1t8r2KtmXQGtQ+QuamTXp7LE63s8FfhUcGJZvbdlQ/jbq/CI/sw1G/5aKxbXWdTo960Mi6eo/GLh8cqcHliUY6l9Y91+2r0+rQa/TotGorU91yXSfURpfrrav76/S4zzoDes0Rz/W2e8NFXZipaxlprAWldrmp07m8uvlVluYsbXRbTnUOLo+Lawdd2+Sfk9ak0zTsGzKx60TOTD+TLUVb+CTjk+Meu6N0B0sOLuGvQ/7aliX6J0MQTHtcXfa4oXS/2g/U+9gG5QegKk99gDrmz+pXoKa2D+KSJyBrBVy5EM499hWIxmg0Gkw6EyadiUgiW+Vbcnvc1LhqKLeXU1ijdhMoqinyTnlUWFrITo6es/JkyGWvdvbOtnd4Zv0zR284Xp+fulBzMo7s83OkKQ+SPWgOhVV2esSFEBkcwGNlCL+kKIo3JB0r0IiGjupbeAQNGrZcefRcer6wvWQ7/17z7wYzsR+PP9Xud/Ytgf/NPuFuRSPvRDvpXmJC/Xuk+7b4/JbZINuJR/Ewf+18b/AZGje0wfYxtoZjytxUXu+yw/FCS1MtefIE25/gH19v44JXVjL0sZ/4Yn0ONmfbX9sVoqk0GrWlpq7FpK61RoLPsd009KYGr8ckjjnudl8osZbw8MqHufTbS9lYuJEgfRC3DLuFGwc3vAvNH2v3W028IzhmzTMs2Hj0EAKBQMJPO3B6nNz/2/28t+M9AO4cfif/m/k/5gydgwa1c/Pr+UXMKStHoyjMKSvnxvp9fqa07Lp4A0eeo/vko7ab9Yeboe/8bDOjnviZh77exrZc6f8hREd045Aba3/PaLh56M28Pu117+uW3ObemhxuB+9uf5ezF5zNl3u+REHhrO5nsXD2Qm4YfANzhs3x29r93oGVJ94H+K/uUkyGwBwmo0Nc9nrkkUd49NFHG6zr06cPu3btavI5fHXZy+lxcs+ye/j54M/oNXoeG/8Y5/Q4R71W+/MjsPIFdUedCdz2o0/Qmregn2Ccn8xiC1OeXgqAUa/F4To8XsaApHAuGZnCrKFdiQgKoM6FQohW5fK4WLhvIa9sfoU8i9onZUDMAP4+6u8M7TLUt8V1Fo9EnHifqHS4bVObl9IaAvZW90ceeYTPP/+cn3/+2btOr9cTGxvb5HP4Ivy4PC7+vuLvLM5ajEFr4Lkpz3Fq8qlgr4YvroPdP6g7jrsVBp4Pr01ueIJTroRZL7RLraD2qbjizTVU2px8+bdxrNpfwidrs/lxewEOtxqETHotMwYmcPHIVMZ0j5ZLDkKIJvEoHn488CMvbXzJe2t7l+Au3Dz0Zs7teW6TO3GLEyjPhucGNm3f1uhP2g4CenoLvV5PQkLHmanX7XHzwG8PsDhrMXqt/nDwqcyDjy6GvM2gN8PsV9Tg09g12tAu7VqzRqPh/etGo9QOhz6xVxwTe8VRZnGwYGMun6zNJqOgiq82HeKrTYdIjQ5m9rCuzB6aRPe40HatVQjRMSiKworcFby48UV2laqt9ZGmSK4bdB0X97kYs97/J9bsUBbe2vR9lzzRIcJPW+gw4WfPnj0kJSVhNpsZO3Ys8+bNIzU19Zj72+127PbDl5EqK1tjwMCmURSFR1Y9wveZ36PX6Hlm0jNq8CnarfbAr8yF4Fi49GNIGVl7VCMtKCHtG37qHNmaExVi5JoJ6Vw9vhubcyr4ZG02Czcf4mBpDS/8socXftnDkJRIZg9N4pwhScT6+Z0DQoi251E8LMtexpvb3mRz0WZAHYzwygFXckW/Kwg1yh9Mrc5pg32/Hnt7cCzUFB9+3Rr9STuoDnHZ64cffqC6upo+ffqQl5fHo48+Sm5uLtu2bSMsLKzRYxrrJwS0y2WvFze+yGtbXkOn0TF/0nzOSDsDCnbAe7PAUgSxveGyTyE6/fBBhbvg5dENT3Th22qrkB+qcbj4aUcBCzbmsmJPMW6P+mOk02qY2CuW84Z15Yz+8QQbO0y+FkK0AqfbyXeZ3/HOtnfYV7EPAJPOxGV9L+OagdcQaY70bYGdWMWSF4lY9qD3taIzodEZwFGtTl761xWwfH6Hm9MxYPv8HKm8vJy0tDSeffZZrr228QG6Gmv5SUlJafPw82nGp/zzj38C8MjYR7ig9wXqJa73ZoO1FBIGq5ONBh8xKF3RbnhpZMN1V30H3Sa0Wa2tpajKzrdb1Ethm7PLveuDjTqmD0hg9rCujO8Rg14n1/SF6KwsTguf7/6c/+34HwU1BYA659zFfS7m8v6XExvU9D6aovmslmqC5ndtfGNcX7j+VzCGtG9RrSSg+/zUFxkZSe/evdm7d+8x9zGZTJhM7Xv5ZVn2Mp5YrY7Jc+OQG9Xgk78N3j0HbBWQdApc8SUERR19cGMdh3102au54sJMXD0+navHp7O/qFrtE7Qxl4OlNSzYmMuCjblEhxg5o188Zw5KYHyPWIx6CUJCdAbF1mI+2vURH+/6mEqH2r0gLiiOy/tfzp96/4kwY+Ot86L1ON0eVrx2G9Ma22gIhj+922GDT1vpkOGnurqaffv2ccUVV/i6FK+M0gzuXn43HsXD7J6zuWnITWqv+w8uVINP8ii4/HMwH+MWxMbudAiNa9ui20D3uFDmntGbO6b2YsPBcr7elMu3W/IotTj4ZF02n6zLJsysZ2q/eKYPSGBynzjMATrOhBAdlaIorCtYx6cZn/LzwZ+9E692C+/GVQOu4pwe52DUySjx7cHtUXjow2XMq/i88R3Oeha69G3fojqADhF+7rrrLs455xzS0tI4dOgQDz/8MDqdjksvvdTXpQFQbivntiW3YXVZGZM4hofGPoTGWgbvn6/OrRLXD/786bGDz7F04GvjGo2G4WlRDE+L4qGz+7Mms5QftuWzaHs+RVV2b4tQkEHHlL5xnDkwkdP6diHU1CF+JIUISJWOShbuW8hnGZ95+/MADI4dzFUDr+K0lNPQaeWPmfbi8Sjc/8Um/r1vduM7DLschvrH56S/6RCfNDk5OVx66aWUlJQQFxfHhAkT+OOPP4iL833LiMvj4q7ld5FbnUtyaDJPT3oag9sFH10CxbshvCtc/kXjl7rqa6zlp5OMoaPXaRnXM5ZxPWN5dNYANhwsU4PQtnxyy618vzWf77fmY9RrObVXbG2LUBfiwuSuMSH8wfaS7Xya8Sk/ZP6A1WUFIEgfxMz0mVzc52L6xfTzcYWBR1EUHl24nR6b5zf+Sd6lP8yY3+51dRQdssNzS7TVIIfz187nvR3vEaQP4v2Z79M7vDt8+hfI+E5t6blmMXRpwi+GsgPw/GB1OTgGZv0f9J3ZanX6I0VR2Jpb4Q1CmcUW7zaNBgYnR3Jany6c1rcLA5LC0Wo7RxgUoiOodFTyU9ZPfL77c7aVbPOu7xnZk4v6XMTZ3c+W/jw+4nJ7eGDBNiwbPuX/jC8e3pA4RL3BRquHm/6A2F6+K7IVSYdnP/PTgZ+883U9MeEJekf1hsUPqMFHZ1LH8WlK8IGGrTzTHu/0wQfUS2ODkyMZnBzJPdP7kFFQxQ9b8/llVwHbcivZnF3O5uxy/vPzbuLCTEzpE8ek3l0Y1yOGqBDpTyBEa7O77SzPWc53+79jec5ynB4nAAatgTPSzuDiPhczrMswGdndh6wON7d9vJGqXb/yUf3gM+A82P6VunzefztN8GkrEn5aKKcqh4d/fxiAqwdcrY7ls/VzWPV/6g7n/xfSxjX9hPUvewVgR0GNRkPfhHD6JoRzxxm9Kai0sTSjkF92FvLb3mKKqux8ui6HT9floNHAwKQIxveMZWKvWIanRUmnaSFayO1xs65gHd/t/46fD/xMlbPKu61nZE/O6XEOs3vOJtocfZyziPaQW27lhvfWYcxbz/vGpw9v0BpqBzdUYPjVMOhCn9XYUUj4aQGn28k9y++hylnFkLgh3HLKLZC/Fb6+Wd1hwlw1hTdLvb+kZLh34sPNXDwylYtHpmJ3uVmTWcqSXUX8treI3QXVbM2tYGtuBa8u24dJr2VUejTje8YyoWcs/RPlEpkQx6MoCrtKd/Hd/u/4IfMHCq2F3m3xwfHM7D6Ts9LPok90Hx9WKepbk1nK395fT1zNXt41PUUI9SbC9jjVu4oTBsGZ//JdkR2IhJ8WeH7D82wt3kq4MZynTn0Kg6MGPrkcXFbocTqc9uCJT3Kk+i0/En4aMOl13nnGAAorbfy2t1h97CmmsMrOij3FrNijDtseZtYzIi2KUekxjEqPYlDXSBlXSAQ8RVHYWbqTpdlLWZy1mP0V+73bwoxhTEubxlndz2J4/HCZZNSPuD0K/12+j2d/3E0fZT8fmv9NuGI5ekdjmDqej0E+P5pCwk8zrc5bzbs73gXgn+P/SVJIInx+DZRlQUQqXPAGtORWzwbhR+5yOp4u4WbOPyWZ809JRlEU9hZWs2JPMb/vLeaP/SVU2VwsyShiSUYRoM5EPzQlklHp0YzsFs0paVFyS70ICA63gzX5a1iavZSl2Uu9Iy+DOuXEpORJnNX9LCZ0nSDj8vih3HIrcz/ZxOrMUkZqdvFe0DMEeSzqFEmlmWqLT51zX4SYHr4rtoORT4BmqHRU8uDvaqvORb0v4rTU02DDe7D9S9Do4MK3jp62oqnqdyCU8NNkGo2GXvFh9IoP45oJ6bjcHnbmVbE6s4S1WaWsyyqjxOJgdWYpqzNLAdBqoHd8GEOSIxmcEsGQ5Ej6JIRhkOk3RCdQZitjRe4KlmYv5ffc36lx1Xi3BemDGJ80nimpUzgt5TSZXNRPeTwKH645yL8X7aLK5mK6cQsv6Z9D77FB2nhIHQsr6vX5GXl9C7paBDYJP83wr9X/It+ST0pYCneOuBOK98D3tRPDnf6PejO0t0CAd3huLXqdlkHJEQxKjuC6id1RFIV9RRbWZpWyNrOUNVml5JRZ2ZVfxa78Kj5Zlw2orUMDksIZnBzJ0JRIhqRE0i0mWO5qER1CVkUWS7OXsiR7CZuKNuFRPN5tXYK6MCllElNSpjAqcRQmnfxx5c+25VbwwFfbvPMk3hW7mjk1L6HxuKDXdJj9Mrw06vABiUNh+hM+qbUjk/DTRD8f+JmF+xei1Wh5csKTBGuNsOCvaj+f7pNh3G0n9wXksleb0Gg09OwSSs8uoVw6KhWAgkqbeht9TjlbcirYnF1Opc3FhoPlbDhY7j023KxnUHIE/RLC6ZcYTt/EMHp2CcWklzvLhG9V2CtYm7+W1Xmr+SPvD7Iqsxps7xPVh8kpk5mSMoV+Mf2kD08HUFxt58Vf9vC/Pw7gUSDSBB+lLaTfwY/UHQZeCOe9CitfhJoSdZ0pAv70jnxmtICEnyaosFfw+B+PA+pt7UO7DIVl8yF3vfrDd+7LoG3FXy7S8tOm4sPNTBuQwLQBCYDaxJxVYmFLTgWbakPR9kOVVNpc/L63hN/3lniP1Wk19IgLoW+9QNQvIZz4cJO0Eok2U+OsYX3Betbkr2F13mp2le5C4fD4tHqNnpEJI5mcMpnJKZNJCk3yYbWiOSptTt5Yvp83fsukxuEG4OIBwTzmeBrTwd/UnaY8AKfeDY5qNfzUmf0SRKf7oOqOT8JPEzy97mlKbCWkR6Rz09CbIG8LLKu9nXDmfIjoevJfpMFlL8PJn080mVaroXtcKN3jQpk9TH0vHS4Puwuq2JZbwa78KnbmVbIzTw1Euwuq2V1QzTebD3nPERlsoHeXMNJjQ+geF1J7vhBSo4OlL5FoNofbweaizazOW82a/DVsLdqKS3E12Kd7RHdGJ45mVMIoRieOltGWO5hSi4P/rTrAOyszKatROy4PSY7gnyPsDP7jNig/AMZQdcDCfmerB615Daxq30XG3AT9zvFR9R2fhJ8TWHloJV/t/QoNGh4d9yhGjQ4W3gYeF/Q9GwZf1DpfqN41emn58T2jXsvArhEM7Hp4MlpFUcivtNUGIbXP0K68SvYXWyivcbImS+1TVJ9eqyE1OrhhKIpVn2NDjdJaJAB1jsCdJTtZnb+a1Xmr2Vi4Ebvb3mCfrqFdvWFnVMIo4oJ9P7ehaL59RdW8+VsmX6zPwe5Sf+/3iAvhnjN6MK3sIzSL/wWKGyLT1FkC4vurB9oqD7f6dB0OUx/10XfQOUj4OQ6ry8pjqx4D4JK+lzCsyzBY/Roc2gCmcDjrmdabfLT+eQxBrXNO0ao0Gg2JEUEkRgRxWt9473qb083ewmr2FVWzv8jC/mIL+4uqySy2UONwq6+LLfyyq+H5wsx6useGkBwdTHJUEMmRQXSNCiI5KpiukUGEyO34nVZhTSFbi7ayuXgzW4q2sKNkh3fC0Dox5hhGJY5iTOIYRiWMIjks2UfVipPldHv4dVchn67N5teMQupm1BzUNYLrJqZzVrID/dfXQfYf6oYB58PZzzacEHvXd2AtA3NkbT8f+SP5ZMhv1+N4fcvr5FbnEh8cz22n3AaVefCLGoY4/SEIS2i9LxYUBafeo44RZI448f7Cb5gNuqNaieBwS1FmkYV9tYFIDUfV5JRZqbK52JxTweacikbPGxVsoGtUEF0jDwei5KjDASkiSC6PdgQ2l42dpTvZUrRFfRRvId+Sf9R+YcYwRsaP9Aae7hHdpWWwg9uZV8ln63L4elMuJRaHd/3Ufl24bmJ3RneLRLPpA3jtfnBUqX9Uz3xavaJw5HufNla92jDmbxCZ2s7fSecjs7ofQ2ZFJud/cz4uj4vnJj/H6Wmnq7O17/gauo6Aa39s2WCGQqC2Fh0srWF/kYWcshpyy63kllnJKbOSW26lwuo84TmCjTriwkx0CTPRJcysLoery13qLUcFG+RDtJ0oikJ2VTabi9QWna3FW8kozTiqv45Wo6VnZE8Gxw1mcOxgBscNJj0iXe7K6uAURWFHXiWLt+WzeHsBGQWH50mLCzNx/rCuXDQyhR5xoWrf0e/vguzV6g6pY9X+PVFpPqref8ms7u1EURSeWP0ELo+LiTVWTstcD3aXGnwAkoZK8BEnxWzQ0Ts+jN7xjXdSrbI5yS23klOqhqEjA1KJxUGNw82BkhoOlNQ0eo46Bp2GuFATceG1oSjMREyIkchgI1EhBvU52EhUsLocbtZLWGoCRVHIs+Sxu2w3O0t3srVoK1uLt1JuLz9q3xhzjBp0asPOgNgBhBhC2r9o0XLLnoIlT6p3Xk262/vacerfWdn1GpbvLubHHfnklKmXL2/RfckPps9ZGHM1oWfcx+T8d9AtuxCMd4DDAmtfV/t6GkJg8t9h7Bz5XGlH0vLTiB+zfuTOZXdi8nhYkJtHissNOhPU74A45QGYdE8bVy1E46wONwWVNgqr7BRW2SistHuXi6rsta9t3rtImkOn1RAZZCAyuC4YHfl8eDnMrCfUpCe09tmk13bK4GR1WdlbtpfdZbvJKMtgd9ludpftpspRddS+Bq2BfjH9GBw7mCFxQxgUN4ikkKRO+e8SMJY9BUsODyRYmTie8Lzfva+fcV7Ii+7zATAbtPw7dhHnlr1z+Pj0SZC57OjzDjgfpj3eOncMd2Jt0fIj4ecINpeNc786l0OWQ/ytrIKbyhvvjwEaeKS8TWoVorU4XB6Kq2uDkTcs2Sm12CmrcVJR46SsxkF57XPdOCMtpddqCDGpQSjMrPcu1z1CvEFJR6jJQKhZT4hRh9mgw2zQ1j6rj6C6dXodWm37BAdFUci35HtDTkapGnQOVh1sMGqy9/vV6EmPTKd3VG8GxQ5iUOwg+kb3lXmyOgmn20NWsYWer6SgqTeukqI07JLjQcPfBy7jtL7xTOodh/nJmIb7A/V/ghXA+ecFuNNOxe1R8Cjg8njwKLUbwTuO05Gf0PU/shXUE9ffR6MBrUaDXqtBr9Wi16nLWo360GjocEFcLnu1g7e3v80hyyESdMFcXZF97B2nPNB+RQnRQka9lqTIIJIim3YHoc3ppsKqBqEyi5PyGgdl3oCkLpfX1K13UG13YbG7qbarfVpcHoUKq7NJfZaa+32Y9VqC6oKSXofZqMOs13qDk0GnxahTnw16TcPXOvVDQH2twaDXgsZJmfMgxY4sCmyZ5FszOVSznxr30a05AOGGSNLCe5Ie1ov08J70iOxNWmg3gowmdBoN2toPGIsNajQONGhAo34YaVA/cNRn0KDxfnjWf123X90HnMLhDzYFpcGHXGPr1f0PH+up/WBVP2DVh9uj7u/2KLgVBUVRcHsa7tPgmNr9PB7U4484RqldV/+8h78Wh8/pUXArdTXVnfPw1zp8HmrPX/e1jq7f0+A89ZY9jdRfW3td3fW/Xt2yw+3G5vBgdbqxOt3YnG41iNS6RXcBdxo+974+Mjv8x3Uhn63L4dN1OY3vf8TP0nPOC3j+TSuw+Fg/8j6hqbegqXtd7+dW613WoNWqr3W1yzqNtvZZg07b8GGoDWHqf5vqf5em2mejXovJoMWk12Gq/e/ZVPtHT5BBi9moQ3FYj1Fxy0n4qSffks9bW98C4M7xjxLkfrnxpsruk9VrvkJ0MnWtLvHh5mYd5/EoWBx1QchJtd1Ntc1Vb9mJxeGmqnadxa4uW+wuahwubE6P90NHfXhwuA+3tDhcHhwuD5U213GqaLQyNPpKtKYitMYitMbi2uciNIYKNJqjG74VRYvH3gWPPQG3LRGPPRGPLZEqdyi5aFjp3TO39iE6uxfd5zNau5Px2u0Ngo+iwG+egbzoqj+pqMIOJY1qxUyoxtbgPHX7P+++oH0Kbyal3sLhZaWxPdqVx378fo0tIeGnnhc3vojNbeOULqcw/cDmxoMPwP6l6vQWEoCEANRRssPMBsLMBqB5welY3B7lcBhyebA61GW7y43V4ald78bqcFNhr6LIlkORPZsSey6ljhwqXIeodOXhxnHMr6FXwjB5kjF4uqJzJaFzdUXnjMej6PF4FFwouA0KHn1tS4TncCuHp95yXWuEP9Nq1P5cmtq/1tVldV3dJRFd7V/z2trtWo363upq12lr1x0+5sjj6+2jqf1a2mN8Xe9yM75uvRa2o46pv0/tfnX7HP7ata0W9Zbf/+MAv+wsbPTfLD02hNToIM6v+pAJZduP2q7RwETdNr7s/Tv7e1xB77yFdM/6iNCq/Y2er27/dWM3UjHydnRaDfq6lhJdbd319gW1VVDh6NYmTe3lLk1te03ddkUBt+LB41Ev27k8HpweBbdb/Tl1uWtfexRcta/dCupz/XW1yw6XB2ftMS63B6dbweH2eNd7t9euV/fx4HIrON2K+vXdincfV21Nrno1eP97qtdyV/esKEqbRC7p81NrZ8lOLv72YhQUPpz5IYNemsjxU670+RGiPTk9TnKrcsmqzOJA5QEyKzLJqswiqyKLElvJMY/Ta/WkhKXQLbyb+ohQn9PC04g2R7da/wel9hJN3aWnuktWDS5L1f5SP3If6u1X1z+k/odf3fUIzeFF72W0w+sbfgiCGjp0HbSfR3uqqHESZtYfu2/ZI5GcsNXDGKaO1dMk8vnRHNLnp40oisIz659BQWFG+gwGxQ1Se+Fv/+LwTt0nqy0+dabc395lCtGpKYpCia2E3OpcDlUf8j7XLedU5Rw1Xk59sUGxDcJNekQ6aeFpdA3til7b9r/qNBpNveAhQaMjiQg+wYChU+5vcLdXo3dvOaogtjeMugGq8mHF04e3yeeH35HwA6w6tIrVeasxaA3cOuxWcDvVKSzqTHmwwbgOTLlfbnMXopkURaHUVqqGGUttwKnK9S4fqj501HxWRwrSB5EWnuZtuekW0Y30cDXkhBpD2+k7EQFn0j1qE93SeZA8Eorqz1WjgT4zYdT1asipS8B6U6PjAsnnh38I+MteHsXDJd9ews7SnVze73LuHXUvrH9Hnbw0JA5u3QQm+aUqxIkoikK5vZxD1YfIqc45qvXmkOXQUfNXHUmDhviQeJJCkkgOSyYpNImkkCSSQpNIC0+jS3AXGQVZtK/CXbDtc9j6GZRlHV4fFAWn/AVGXCujMrcxuezVBn488CM7S3cSYgjh+sHXg9OmJnSACXMl+AiB2t+muKaYgpoCiqxFFNYUUlBTQGFNIYU1hRTVFFFQU9CkcBMXHEdyaG2wCU2ia2hX73NCcAIGncxZJnys7ABsXwBbP4eCrYfXG4LVVp5Bf4Iep8nkoh1YQIcfl8fFSxtfAuDKAVcSbY6G1f+FylwIS4IR1/i4QiHalqIoVDoqvUGmLsTUhZq6R6mtlKbecxEXFNcg0NQ9dw3tSkJIggwAKPyPywEHV8Hen2DPTw0va2n10HOqGnj6zACjTEvSGQR0+Fm4byFZlVlEmiK5ot8V4KiB5bWd1CbdDYbWuWVXiPbmUTxU2CsotZVSYi2hxFZCUU2RN8x4w4616IT9bOrotXq6BHUhLjiOLsFdiA+Op0uw+rpuOSEkAZPO1MbfnRCtoCL3cNjZvxQc1Ye3abSQOg4GXQj9z4XgaJ+VKdpGwIYfp9vJf7f8F4BrB16rdpb8/XmwFEJkGgy93McVCtGQw+1oEGbqlkttpepra6l3fZmtDLfS9KkqIk2RRwWZuoATF6SGnShzlPS3ER2XrRJy1qp3ae35GQqPGLcnJA56ngG9pkL3KRJ4OrmADT9f7/ua3OpcYswxXNz3YnWW3d+eUzdO/rtcyxVtzqN4qHJUUWYrOyrMNBZsqpxNHUPksAhTBNHmaKLN0d4Qc2SrTZfgLtJaIzqfykPqpayDf6jPBdvVWdS9NOqdW73OUB8JQ0Ar4T5QBGT4cbqdvL7ldQCuHXQtQfogWPUyWEshKh0GXeTjCkVHY3fbKbeVU24vp8JeQYWjwrtcf325/fBypaOyWa0zoF56ijZHE2OOITpIfY4xx6jrgmK8QScmKIYoU5R0HhaBweNR++kcXAXZq9Xn8oNH7xfVDVLHqn14epwmrTsBLCDDz8L9CzlkOURsUCx/6v0ncNlh5Qvqxgm3gy4g/1kEaif4KkcVlY7Kw0HmiNBy5HKFveKEdzkdT6ghtGFwqRds6oeZaHM04cZwGalXBDaPG0r2Qv5WyNusPh/aCLbyhvtptJAwSA07qWMgZQyEJ/qkZOF/Au5T3uVxeVt9rh5wNWa9WR3XpypPvcNryKW+LVCcFI/iweK0UOmopMpRpQYZe+Xh1071tXebo7LBvjWulk+gp9PoiDBFEGGKINIU6X0+0bLc/STEMThqoHAH5G+BvC1q0CnYDo39sWEIVi9jpY6F1NHqsims/WsWHULAhZ+fsn4ipzqHaHM0f+rzJ3C7Dvf1GXeLOiqn8AlFUbC6rFQ7q9WHQ322OC1UO9TnKmcV1Y7qBoGlfoipdlQ3+Zbs4wnWBzceWsyHlyOM9baZIwgzhEmrjBAt4XaqAwgW74Hi3WrIyd8KJXuO6KdTyxAM8QPVlp3EwZAwWF2Wy7yiiQIu/Ly7410Aruh/hdrXZ8tnUJYJQdEw/EofV9cxuTwualw11DhrsDgtalipDS/eZWc1Fsfh5fqvveHGZcHT2C+6FjDpTIQZwwgzhhFuDG+wfOTrI59DjaHtMheUEAFFUaCm5HDAKdkDxXvV57Is8Bxj3raQuMPhpi7oRHcHra5dyxedS8D9hs+qzCIyPJKL+1ysdpL77Vl1w5ibAmLwKkVRsLlt1DhrsLqsWJwW77PFaaHGVXP0OmdNw3Djql1Xu76p48Q0lU6jI8QQQqghlBBj7bMhhDBDGMGG4KMCS5gxjHBTw3Vy95IQPqAoat+b8mw10NQPOMV7ju6XU58hGGJ6QEwviB8AiUPUwBOW0E7Fi0AScOEH4JK+lxBmDIOMRer1ZGMYjLrO12V5KYqC3W3H6rJic9mwuqxY3VasTis2t/q6LrzUuNRnq/Pw8pHb6l7XPVrjslBjDFoDwYZgQg2h3sASagz1Bpn667yv64WbUEMoocZQzDqzXD4Swh8pCliKoeKgejdVeTZUZDdctlce5wQaiEiB2J5qyIntBTE91eewJLnVXLSbgAs/Jp2JP/f7s/pi1f+pzyOuUiepawKXx4XdbcfmsmFz27C77FjdVuwudZ3VrQaWuu11Aeao1/WW6weTunVtFVDqM+vMBBuCCTGEEKyvfTYEH3M5xBBCiD6EIENQg2PqluW2aiE6OHs1VBdAVb46Ts5RISe78c7GRwqOhcjU2nDT63DYiekBhqC2/z6EOIGACz+9o3rz5Z4vsVUcxF65FVtMNHZDDfZl96hhpjbYNAg4brs35LiOdV26jRi1Rsx6M0H6IO/DrDcTrA8m2BDsXResV5ePty5YH0yQQX02680yWq8QgaDuUlRVAVTnN/Jc+6guaDjFwzFpICwRIlPUVpzI1NrlVHU5IhmMwW39XQlxUgIu/Gwt3soOyw71RUS4+nzwpxady6wzY9KbMOvMmPVmTDqTN5wE6dRns96MWWf2rq973VigCdIHEaQLIsgQ5N1HOt4KIY7idqmDstaUqJehakpqH6XqFD11Yabu2WVr+rkNIWo/m7DEesGmXsgJT5YR8EWHF3CfrNPSphERZMC49QvMHjemoZdjjk7HpDNh0pkahJO6dUH6oKO2GXVGaTkRQpw8RQFbxeHwUlM/zJSApaTh65pidf/mMkdAaAKExavBJjReDTne59ptMjaOCAABF34eG/8Y4Wufh5ISdcTPSU/4uiQhREfncqiXlmwV6sNa3vD1UduOWN+iy+kata9iSCwEx9Q+otX+NmGJapCpCzSh8dLXRoh6Ai784KiBdW+py2Nv8m0tQgjfcrvAUaV29HVUg71KfTiq1XX2qobbvaHliBBzEtObeBlDa8NLjBpg6gJNSEy9cFNvW1CkjHUjRAsFXvjZ8bV6rTwyFfqe7etqhBBN5XGDwwLOmsPPTmu9dTXgtDQMMt5AU7eust5ydeuElvpMEerlpaAIMEeqy97nCDWw1C3X3xYcAwZz69YihDimwAs/699Rn0deJ381CdEaFEWdnsBlBaftxM8NwktdaKk5OsQ4asNN3XIrD6bZgM6otryYQsEUfnjZGKr2gTGFqcvm8GOHGVO4/E4RooMIvPBTuB2CzTDsCl9XIkTrUBRw2dVw4HLUPtvB7Tjiuf72Y+zntKp3Bh3z+RihppWmJWkajToauyFY7cdSt2wMVu9U8oaWUHUAU1NY40HGuz1U5vQTIsAEXvgBGHiBem1diOaoa+FwO8DjVPuLeGpfnyhUNDV8HHO/4+zvcfr6X6YhfZB6CedYz4bgIwJLcCPrQo7YVm+d3gwyArgQ4iQEZvgZ6T9TWXR6Ho96J0uDh7uRdbUPt1N9eOqeXeqHvHf5GOHDu3zk8UeuczXz+HpfX3H7+l+zabQGtSVDZzzi2aSOz9Loc+1+hiA1XBzz+TihRh+knkeCiRDCzwVe+EkYDMnD2+bcHo/6Aelx1z7XftArnhOsq12vuBueo7F19c9xVIg40esTBI8mneOI1+4THN8O03T4lFavhobGAobuyBBynNDR5LBygv10RpkfSQghTiDwwo/TBu/OOkYgqR80Glt3gpAimk6rr/fQ1Vs2qKFBZ6hdrg0XdctaQ23YaGxfQzOON9QGl9rXdcuNnbPB+fUNj5dWDiGE6HACL/yU7IbqPT74wprDH/Ianbqs0R57nUZ3OBgcc129AKHTc+xAcax1je1zkuc4qo5G9tFoJTQIIYTwmcALP2c/B2Fh6qWBBmFCV29dS0OK7hjn1cmHvRBCCOEnOlT4eemll5g/fz75+fkMGTKEF198kVGjRjXvJIMuhPDwtilQCCGEEH6vw/SM/OSTT5g7dy4PP/wwGzZsYMiQIUyfPp3CwkJflyaEEEKIDqTDhJ9nn32W66+/nquvvpr+/fvz6quvEhwczFtvveXr0oQQQgjRgXSIy14Oh4P169dz3333eddptVqmTp3KqlWrGj3Gbrdjtx8eDr+iogKAysrKti1WCCGEEK2m7nNbUVpv6JQOEX6Ki4txu93Ex8c3WB8fH8+uXbsaPWbevHk8+uijR61PSUlpkxqFEEII0XZKSkqIiIholXN1iPDTEvfddx9z5871vi4vLyctLY2DBw+22j+eaJnKykpSUlLIzs4mXDqf+5S8F/5F3g//Ie+F/6ioqCA1NZXo6NablqpDhJ/Y2Fh0Oh0FBQUN1hcUFJCQkNDoMSaTCZPp6MkKIyIi5AfZT4SHh8t74SfkvfAv8n74D3kv/Ie2FUev7xAdno1GI8OHD+eXX37xrvN4PPzyyy+MHTvWh5UJIYQQoqPpEC0/AHPnzuXKK69kxIgRjBo1iueeew6LxcLVV1/t69KEEEII0YF0mPBz8cUXU1RUxEMPPUR+fj5Dhw5l0aJFR3WCPhaTycTDDz/c6KUw0b7kvfAf8l74F3k//Ie8F/6jLd4LjdKa944JIYQQQvi5DtHnRwghhBCitUj4EUIIIURAkfAjhBBCiIAi4UcIIYQQAaVThZ+XXnqJbt26YTabGT16NGvWrDnu/p999hl9+/bFbDYzaNAgvv/++3aqtPNrznvx+uuvM3HiRKKiooiKimLq1KknfO9E0zX3v4s6H3/8MRqNhtmzZ7dtgQGkue9FeXk5c+bMITExEZPJRO/eveX3VCtq7vvx3HPP0adPH4KCgkhJSeGOO+7AZrO1U7Wd0/LlyznnnHNISkpCo9Hw1VdfnfCYpUuXcsopp2AymejZsyfvvPNO87+w0kl8/PHHitFoVN566y1l+/btyvXXX69ERkYqBQUFje7/+++/KzqdTnnqqaeUHTt2KA8++KBiMBiUrVu3tnPlnU9z34vLLrtMeemll5SNGzcqO3fuVK666iolIiJCycnJaefKO5/mvhd1MjMzla5duyoTJ05Uzj333PYptpNr7ntht9uVESNGKDNnzlR+++03JTMzU1m6dKmyadOmdq68c2ru+/HBBx8oJpNJ+eCDD5TMzExl8eLFSmJionLHHXe0c+Wdy/fff6888MADypdffqkAyoIFC467//79+5Xg4GBl7ty5yo4dO5QXX3xR0el0yqJFi5r1dTtN+Bk1apQyZ84c72u3260kJSUp8+bNa3T/iy66SDnrrLMarBs9erTy17/+tU3rDATNfS+O5HK5lLCwMOXdd99tqxIDRkveC5fLpYwbN0554403lCuvvFLCTytp7nvxyiuvKN27d1ccDkd7lRhQmvt+zJkzRznttNMarJs7d64yfvz4Nq0zkDQl/Nxzzz3KgAEDGqy7+OKLlenTpzfra3WKy14Oh4P169czdepU7zqtVsvUqVNZtWpVo8esWrWqwf4A06dPP+b+omla8l4cqaamBqfT2aqT2AWilr4Xjz32GF26dOHaa69tjzIDQkvei2+++Yax/9/e3YZE0bVxAP93q6NimdDmS6GFhr2YtmEZlmBhBUWJECmpmxYmon6slEwULJOQEKMMpYwKshd6ETeiMv2QUqCtIGWJiqXlaqLVUmK1e+4P99Py7K3duNu6azv/H8yHzpwzc41X61ycObNGRiIrKws+Pj5YuXIliouLodfrbRW2w7IkH+vXr0dra6vx0VhPTw/u3buH7du32yRm+oe17t1/zDc8/5fh4WHo9foJ3/bs4+ODV69eTTpGq9VO2l+r1U5bnHJgSS7+LScnBwsWLJjwH5zMY0kunjx5gvPnz6Otrc0GEcqHJbno6enB48ePkZSUhHv37qGrqwuZmZn4/v07CgoKbBG2w7IkH4mJiRgeHkZUVBSEEPjx4wcyMjJw5MgRW4RM//Ore/fnz58xNjYGd3f3KR3HIWZ+yHGUlJSgpqYGt2/fhpubm73DkRWdTgeVSoWqqiooFAp7hyN7BoMB3t7eqKysRHh4OBISEpCXl4dz587ZOzRZamxsRHFxMc6ePYvnz5/j1q1bUKvVKCoqsndoZAGHmPlRKBRwcnLC4OCgSfvg4CB8fX0nHePr62tWf5oaS3LxU2lpKUpKSvDo0SOEhYVNZ5iyYG4uuru70dvbi507dxrbDAYDAMDZ2RmvX79GUFDQ9AbtoCz5XPj5+cHFxQVOTk7GtuXLl0Or1eLbt2+QJGlaY3ZkluQjPz8fKpUKaWlpAIDQ0FB8+fIF6enpyMvLw19/cS7BFn517/b09JzyrA/gIDM/kiQhPDwc9fX1xjaDwYD6+npERkZOOiYyMtKkPwA8fPjwl/1paizJBQCcPHkSRUVFuH//PtasWWOLUB2eublYtmwZ2tvb0dbWZtxiY2OxadMmtLW1wd/f35bhOxRLPhcbNmxAV1eXsQAFgM7OTvj5+bHw+U2W5OPr168TCpyfhangn8i0Gavdu81biz1z1dTUCFdXV3Hx4kXx8uVLkZ6eLry8vIRWqxVCCKFSqURubq6xf1NTk3B2dhalpaWio6NDFBQU8FV3KzE3FyUlJUKSJHHz5k0xMDBg3HQ6nb0uwWGYm4t/49te1mNuLt6+fSvmzJkjsrOzxevXr0VdXZ3w9vYWx44ds9clOBRz81FQUCDmzJkjrl69Knp6esSDBw9EUFCQiI+Pt9clOASdTic0Go3QaDQCgDh16pTQaDTizZs3QgghcnNzhUqlMvb/+ar7oUOHREdHhzhz5oy8X3UXQojTp0+LgIAAIUmSiIiIEE+fPjXui46OFikpKSb9r1+/LoKDg4UkSSIkJESo1WobR+y4zMnFokWLBIAJW0FBge0Dd0Dmfi7+H4sf6zI3F83NzWLdunXC1dVVBAYGiuPHj4sfP37YOGrHZU4+vn//LgoLC0VQUJBwc3MT/v7+IjMzU4yOjto+cAfS0NAw6e//nz/7lJQUER0dPWGMUqkUkiSJwMBAUV1dbfZ5ZwnB+ToiIiKSD4dY80NEREQ0VSx+iIiISFZY/BAREZGssPghIiIiWWHxQ0RERLLC4oeIiIhkhcUPERERyQqLHyIiIpIVFj9EREQkKyx+iGhG27dvH44ePWrR2NTUVMTFxU1ob2xsxKxZs/Dx48ffC46I/kjO9g6AiOhX9Ho96urqoFar7R0KETkQzvwQkc19+PABvr6+KC4uNrY1NzdDkiTU19ebtLm4uKCjowPz5s3D+Pi4yXHi4uKgUqlsFjcROQYWP0Rkc/Pnz8eFCxdQWFiIlpYW6HQ6qFQqZGdnIyYmxtivtrYWO3fuRHx8PPR6PWpra437hoaGoFarsX//fntcAhH9wVj8EJFdbN++HQcOHEBSUhIyMjLg4eGBEydOmPS5e/cuYmNj4e7ujsTERFRXVxv3XblyBQEBAdi4ceN/nqeurg6zZ8822bZt22bSp6+vDxs3bsSKFSsQFhaGGzduWO06iWjmmSWEEPYOgojkaWxsDCtXrkRfXx9aW1sRGhpq3NfR0YG1a9dieHgYbm5u0Gg0WLt2Ld68eYOFCxciLCwMu3fvRn5+/i+Pn5qainfv3qGiosKk/dmzZ0hOTsbo6Ci8vLwwMDCAwcFBKJVKaLVahIeHo7OzEx4eHtN27URkP1zwTER2093djffv38NgMKC3t9ek+KmtrcWWLVvg5uYGAFi9ejVWrVqFS5cuYevWrXjx4sWUFkJ7eHhgyZIlJm39/f0m//bz84Ofnx8AwNfXFwqFAiMjIyx+iBwUix8isotv374hOTkZCQkJWLp0KdLS0tDe3g5vb28A/zzySk9PNxmTlpaGsrIyvHv3Dps3b4a/v7/V42ptbYVer5+WYxPRzMA1P0RkF3l5efj06RPKy8uRk5OD4OBg4+LloaEhtLS0YMeOHSZjEhMT0d/fj6qqqmlZ6DwyMoK9e/eisrLS6scmopmDMz9EZHONjY0oKytDQ0MDPD09AQCXL1/GqlWrUFFRAUmSEBERAYVCYTJu7ty52LVrF9Rq9aRfXvg7xsfHERcXh9zcXKxfv96qxyaimYULnoloxomNjUVUVBQOHz48YV9MTAxCQkJQXl5utfMJIZCYmIilS5eisLDQasclopmJj72IaMaJiorCnj17TNpGR0dx+/ZtNDY2Iisry6rna2pqwrVr13Dnzh0olUoolUq0t7db9RxENHNw5oeI/giLFy/G6Ogo8vPzcfDgQXuHQ0R/MBY/REREJCt87EVERESywuKHiIiIZIXFDxEREckKix8iIiKSFRY/REREJCssfoiIiEhWWPwQERGRrLD4ISIiIllh8UNERESywuKHiIiIZIXFDxEREcnK34jJlHD0tFPaAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "kij_library = {\n", " ('H2','Ne'): 0.18,\n", " ('He','H2'): 0.17\n", "}\n", "lij_library = {\n", " ('H2','Ne'): 0.0,\n", " ('He','H2'): -0.16\n", "}\n", "\n", "def get_model(names, c_factor=0):\n", " param_library = {\n", " 'H2': {\n", " \"Ls\": [156.21],\n", " \"Ms\": [-0.0062072],\n", " \"Ns\": [5.047],\n", " \"As\": [3.0696],\n", " \"Bs\": [12.682],\n", " \"cs / m^3/mol\": [c_factor*-3.8139e-6],\n", " \"Tcrit / K\": [33.19],\n", " \"pcrit / Pa\": [12.964e5]\n", " },\n", " 'Ne': {\n", " \"Ls\": [0.40453],\n", " \"Ms\": [0.95861],\n", " \"Ns\": [0.8396],\n", " \"As\": [0.4673],\n", " \"Bs\": [2.4634],\n", " \"cs / m^3/mol\": [c_factor*-2.4665e-6],\n", " \"Tcrit / K\": [44.492],\n", " \"pcrit / Pa\": [26.79e5]\n", " },\n", " 'He': {\n", " \"Ls\": [0.48558],\n", " \"Ms\": [1.7173],\n", " \"Ns\": [0.30271],\n", " \"As\": [1.4912],\n", " \"Bs\": [3.2634],\n", " \"cs / m^3/mol\": [c_factor*-3.1791e-6],\n", " \"Tcrit / K\": [5.1953],\n", " \"pcrit / Pa\": [2.276e5]\n", " }\n", " }\n", " params = [param_library[name] for name in names]\n", " model = {k: [param[k][0] for param in params] for k in ['Ls','Ms','Ns','As','Bs','cs / m^3/mol','Tcrit / K','pcrit / Pa']}\n", "\n", " if len(names) == 1:\n", " model['kmat'] = [[0]]\n", " model['lmat'] = [[0]]\n", " else:\n", " kij = kij_library[names]\n", " model['kmat'] = [[0,kij],[kij,0]]\n", " lij = lij_library[names]\n", " model['lmat'] = [[0,lij],[lij,0]]\n", " \n", " j = {\n", " \"kind\": \"QCPRAasen\",\n", " \"model\": model\n", " }\n", " return teqp.make_model(j), j\n", "\n", "model = get_model(('H2','Ne'))[0]\n", "modelH2 = get_model(('H2',))[0]\n", "modelNe = get_model(('Ne',))[0]\n", "\n", "def get_traces(T, ipures):\n", " traces = []\n", " for ipure in ipures:\n", " rhovecL0 = np.array([0.0, 0.0])\n", " rhovecV0 = np.array([0.0, 0.0])\n", " if ipure == 1: \n", " rhoL, rhoV = modelNe.superanc_rhoLV(T)\n", " else:\n", " rhoL, rhoV = modelH2.superanc_rhoLV(T)\n", " rhovecL0[ipure] = rhoL\n", " rhovecV0[ipure] = rhoV\n", "\n", " opt = teqp.TVLEOptions(); \n", "# opt.polish=True; \n", "# opt.integration_order=5; opt.rel_err=1e-10; \n", "# opt.calc_criticality = True; \n", " opt.crit_termination=1e-10\n", " trace = model.trace_VLE_isotherm_binary(T, rhovecL0, rhovecV0, opt)\n", " traces.append(trace)\n", " return traces\n", "\n", "for T in [24.59, 28.0, 34.66, 39.57, 42.50]:\n", " if T < 26.0:\n", " traces = get_traces(T, [0, 1])\n", " else:\n", " traces = get_traces(T, [1])\n", "\n", " for trace in traces:\n", " df = pandas.DataFrame(trace)\n", " \n", " # Plot the VLE solution\n", " line, = plt.plot(df['xL_0 / mole frac.'], df['pL / Pa']/1e5)\n", " plt.plot(df['xV_0 / mole frac.'], df['pL / Pa']/1e5, color=line.get_color())\n", "\n", " # Plot the VLLE solution if found\n", " for soln in model.find_VLLE_T_binary(traces):\n", " for rhovec in soln['polished']:\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], p/1e5, 'X', color=line.get_color())\n", " # print(T, rhovec, x[0], p/1e5, 'bar')\n", "\n", "plt.gca().set(xlabel='x/y H$_2$', ylabel='$P$ / bar', xlim=(0,1), ylim=(0,30));" ] } ], "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 }