diff --git a/doc/source/models/MultifluidPlusAssociation.ipynb b/doc/source/models/MultifluidPlusAssociation.ipynb new file mode 100644 index 0000000..2eb5d72 --- /dev/null +++ b/doc/source/models/MultifluidPlusAssociation.ipynb @@ -0,0 +1,140 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "id": "03058ac5", + "metadata": {}, + "source": [ + "# Multifluid+Association\n", + "\n", + "Version 0.22 of teqp adds a new sort of hybrid model -- a combination of multifluid model plus an association model. Conceptually the pure fluid EOS consider the self-association implicitly, so you may want to consider the association model to only allow cross-assocation, which be enforced by supplying the \"self_assocation_mask\" logical mask in the \"options\" of the model. See below" + ] + }, + { + "cell_type": "code", + "execution_count": 69, + "id": "ae2f6b60", + "metadata": {}, + "outputs": [], + "source": [ + "import teqp, json, numpy as np\n", + "import matplotlib.pyplot as plt\n", + "import pandas \n", + "\n", + "BIP = [{\n", + " \"Name1\": \"Water\",\n", + " \"Name2\": \"Ammonia\",\n", + " \"betaT\": 1.0,\n", + " \"gammaT\": 1.0,\n", + " \"betaV\": 1.0,\n", + " \"gammaV\": 1.0,\n", + " \"F\": 0.0\n", + "}]\n", + "jmf = {\n", + " \"components\": [\"Water\", \"Ammonia\"],\n", + " \"root\": teqp.get_datapath(),\n", + " \"BIP\": BIP, \n", + "}\n", + "\n", + "jassoc = {\n", + " \"kind\": \"Dufal\",\n", + " \"model\": {\n", + " \"sigma / m\": [3.0555e-10, 3.3309e-10],\n", + " \"epsilon / J/mol\": [3475.445374388054, 323.70*8.3124462618],\n", + " \"lambda_r\": [35.823, 36.832],\n", + " \n", + " # Note the scaling factors of 0.2 on the bonding energy to yield more reasonable behavior\n", + " \"epsilon_HB / J/mol\": [0.2*13303.140189045183, 0.2*1105.0*8.314462618], \n", + " \n", + " \"K_HB / m^3\": [496.66e-30, 560.73e-30],\n", + " \"kmat\": [[0.0,0.0],[0.0, 0.0]],\n", + " \"Delta_rule\": \"Dufal\",\n", + " \"molecule_sites\": [[\"e\",\"e\",\"H\",\"H\"],[\"e\",\"H\",\"H\",\"H\"]],\n", + " \"options\": {\"self_association_mask\": [False, False]}\n", + " }\n", + "}\n", + " \n", + "j = {\n", + " 'kind': 'multifluid-association',\n", + " 'model': {\n", + " 'multifluid': jmf,\n", + " 'association': jassoc\n", + " }\n", + "}\n", + "\n", + "model = teqp.make_model(j)" + ] + }, + { + "cell_type": "code", + "execution_count": 70, + "id": "469bdad6", + "metadata": {}, + "outputs": [], + "source": [ + "# model.get_assoc_calcs(300, 300, np.array([0.5, 0.5]))" + ] + }, + { + "cell_type": "code", + "execution_count": 71, + "id": "3d9a26b9", + "metadata": {}, + "outputs": [], + "source": [ + "T = 293.15 # K\n", + "pure = teqp.build_multifluid_model([\"Ammonia\"], teqp.get_datapath())\n", + "anc = pure.build_ancillaries()\n", + "j = model.trace_VLE_isotherm_binary(T, np.array([0, anc.rhoL(T)]), np.array([0, anc.rhoV(T)]))\n", + "df = pandas.DataFrame(j)" + ] + }, + { + "cell_type": "code", + "execution_count": 74, + "id": "e9760392", + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjoAAAHJCAYAAACMppPqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABzh0lEQVR4nO3dd3hTZf8G8Ds7bbonHbS0QKGlzLL3RhQQ90RQUVFEAcerIiq+8uJ4fzgriiKIAvIKgsgQAdlDtgJlltGWlpa2dO/m+f1x2pTQAm1pctLk/lzXuZKcnJx8G6C5ec4zFEIIASIiIiI7pJS7ACIiIiJLYdAhIiIiu8WgQ0RERHaLQYeIiIjsFoMOERER2S0GHSIiIrJbDDpERERktxh0iIiIyG4x6BAREZHdYtAhm7Bs2TIoFAosXbq02nPt27eHQqHA+vXrqz3XvHlzdOrUqU7v9eWXX2LBggX1LbVW4uPjodPpsHv3btO+Y8eO4bnnnkOPHj1gMBigUCiwZcuWaq+9cuUKPDw8sHLlSovWeDNHjhyBQqGARqNBSkqKrLVY07hx49CsWbNbOsegQYMwYcIE0+MDBw5g4sSJaNu2LVxdXeHv74/Bgwfjzz//rPbaMWPGYPTo0bf0/rcqJSUFb775Jnr06AEfHx+4ubkhJiYGc+fORXl5ebXj9+7di2HDhsHV1RUuLi4YMGAAdu7cWe24zz77DN27d4ePjw90Oh1CQkLw4IMP4tixY7Wqa/Xq1XjsscfQtm1baDQaKBSKGo87f/48FApFjdtPP/100/dZsGABFAoF9u/fb7Y/PT0dnTt3houLCzZs2FCrmkl+DDpkE/r37w+FQoHNmzeb7c/MzMSRI0dgMBiqPZeUlISzZ89iwIABdXovawSdl19+GUOGDEGPHj1M+/bv34+VK1fCy8sLgwYNuu5rPT09MWXKFLzyyisoKSmxaJ038u233wIAysrKsHDhQtnqsLbp06djxYoV9X79r7/+ip07d2L69OmmfUuWLMHevXvxxBNP4Ndff8W3334LnU6HQYMGVfts33nnHaxZs6bGEHQzp0+fxpQpU9CmTRs4OzvDyckJrVq1wtSpU3Hq1Klan+fAgQNYuHChqb7ly5ejX79+ePbZZ/HUU0+ZHbtv3z707dsXhYWF+OGHH/DDDz+gqKgIgwYNMgv6AJCRkYHhw4fj22+/xR9//IEZM2bg0KFD6NatG06ePHnTulasWIE9e/YgKioK7du3v+nxkyZNwu7du822IUOG1PpzuFpSUhL69OmDs2fPYuPGjfU+D8lAENmItm3bilatWpnt++WXX4RGoxEvvPCC6Nq1q9lzCxcuFADEb7/9Vqf3adOmjejXr9+tlmvGaDSKgoICIYQQcXFxAoD4/fffzY4pLy833f/5558FALF58+Yaz3fp0iWhVqvFokWL6lXP2LFjb+lnLCoqEt7e3qJ9+/YiKChIRERE1PtcjqZr167iwQcfNNuXmppa7biysjLRrl070bx582rPjRgxQgwZMqTW71leXi7eeOMNoVKpRJ8+fcQXX3wh1q1bJ7Zu3Srmz58vRo8eLfR6vfjggw9qdb7MzExRUlJSbf/EiRMFAJGQkGDaN2zYMOHv7y/y8/NN+3JycoSPj4/o2bPnTd+r8t/L9OnTa/VzXltLTc6dOycAiI8++uim56zJ/PnzBQCxb98+IYQQp06dEiEhISIgIED8888/9TonyYctOmQzBgwYgJMnT5pdJtmyZQu6dOmC22+/HQcOHEBubq7ZcyqVCn369AEAzJgxA926dYOXlxfc3NzQqVMnzJs3D+KqdWubNWuGY8eOYevWraam7KsvU+Tk5ODll19GWFgYtFotgoKCMHnyZOTn55vVqlAo8Pzzz+Orr75CZGQkdDodvv/+ewDAnDlz0KRJk2r/41Mqa//Pzd/fH0OGDMFXX31V69c0pJUrVyIjIwPjx4/H2LFjcerUKezYsaPacc2aNcOIESOwevVqdOzYEU5OToiMjMTq1asBSJcAIiMjYTAY0LVr12qXAsaNGwcXFxecOHECw4YNg8FgQEBAAN5//30AwJ49e9C7d28YDAZERESYPuOrHT16FHfeeSc8PT2h1+vRoUOHasdt2bIFCoUCS5YswbRp0xAYGAg3NzcMHjy4WktCTZeuYmNj0bdvX/j5+cFgMKBt27b48MMPUVpaanbcoUOHsHfvXowZM8Zsv5+fX7W6VSoVYmJikJiYWO25MWPGYOPGjYiPj6/2XE2eeuopLFiwAJs3b8a2bdswceJE3Hbbbejbty/GjRuHFStWYNOmTZg9ezbee++9m57P09MTGo2m2v6uXbsCkFo3Ku3cuRP9+/eHs7OzaZ+rqyv69u2LXbt23fSyp6+vLwBArVbftK66/BtqKIcPH0bv3r2hVquxY8cOtG3b1uo10C2SO2kRVVqxYoUAIBYvXmza17ZtW/H666+L3NxcoVarxZo1a0zPhYWFiS5dupgejxs3TsybN09s2LBBbNiwQfz73/8WTk5OYsaMGaZjDh48KMLDw0XHjh3F7t27xe7du8XBgweFEELk5+eLDh06CB8fHzF79myxceNG8emnnwp3d3cxcOBAYTQaTecBIIKCgkS7du3E4sWLxZ9//imOHj0qhBAiPDxc3H///Tf8WW/WoiOEEB988IFQKpXiypUrtfr8rnarLTpDhgwROp1OZGZmijNnzgiFQiHGjRtX7bjQ0FARHBwsoqOjxZIlS8TatWtFt27dhEajEW+99Zbo1auX+OWXX8SKFStERESE8Pf3N7V8Vdap1WpFZGSk+PTTT8WGDRvE448/LgCI119/XURERIh58+aJ9evXixEjRggAYv/+/abXnzhxQri6uormzZuLhQsXijVr1oiHHnpIADBrvdi8ebMAIJo1ayYeeeQRsWbNGrFkyRIREhIiWrZsKcrKysxqCg0NNfs5p0yZIubMmSN+//138eeff4qPP/5Y+Pj4iMcff9zsuHfffVeoVCqRm5t708+4tLRUtGjRQnTs2LHac6mpqQKA+Oyzz256nkWLFglPT09x7ty56x5TVlYmjEajOHLkiDAYDGLv3r03PW9Nxo4dK9RqtUhPTzft02q14rHHHqt2bOWfw/r162usp6ioSBw/flzceeedws/Pz6yVqDZq06Lj7e0tNBqNcHJyEr169RK//vprrc5d2aLz8ccfC3d3dxEdHS2Sk5PrVB/ZDgYdshmZmZlCqVSKp59+WgghRHp6ulAoFKZLQF27dhUvv/yyEEKIhIQEAUC8+uqrNZ6rvLxclJaWinfffVd4e3ubhZTrXbqaNWuWUCqVpubqSsuWLRMAxNq1a037AAh3d3eRmZlpdmzlF9T7779/w5+1NkFnw4YNAoBYt27dDc8lhPSlefX22GOPib59+1bbf/XncD3nz58XSqXS7PJLv379hMFgEDk5OWbHhoaGCicnJ5GUlGTad/jwYQFABAQEmF3OWLlypQAgVq1aZdo3duxYAUAsX77c7Gfx9fUVAEwhVAghMjIyhEqlElOnTjXte/DBB4VOp6v2JTl8+HDh7OwssrKyhBBVQef22283O+5///ufACB2795tVtO1QedqlX+3Fi5cKFQqldnfgeHDh4vWrVtf97VXmzZtmgAgVq5cWePzQUFB4oEHHrjpeVq2bCm+/vpr0+Njx46JPn36CJ1OJ5o1ayZ++OEHERoaavq79vLLL4tHH320VjVebf369UKpVIopU6aY7e/QoYOIiIgwu6xUWloqwsPDq/3HpZJOpxMABAAREREh4uLi6lzPjYJOcnKyeOqpp8T//vc/sX37drFo0SLRvXt3AUB88803Nz13ZdCp/HeelpZW5/rIdvDSFdkMT09PtG/f3jQSaevWrVCpVOjVqxcAoF+/fqYOyZW3V3dE/vPPPzF48GC4u7tDpVJBo9HgrbfeQkZGBtLS0m76/qtXr0Z0dDQ6dOiAsrIy0zZs2LAaR0gNHDgQnp6eZvuSk5MB1Hypoq4qz3Hx4sUbHnf+/HloNBqzbeHChdi2bVu1/Vu3br3p+86fPx9GoxFPPPGEad8TTzyB/Pz8GkfFdejQAUFBQabHkZGRAFDtckbl/gsXLpi9XqFQ4Pbbbzc9VqvVaNGiBQICAtCxY0fTfi8vL/j5+Zm9/s8//8SgQYPQtGlTs3OOGzcOBQUF1TrDjho1yuxxu3btaqzpWocOHcKoUaPg7e1t+rv12GOPoby83KyTb3Jycq3+7L/99lvMnDkTL730Eu68884aj/Hz87vpn/2xY8dw8eJFPPLIIwCAwsJCDBs2DE5OTli1ahXef/99vP3226a/lwAwcuTIOnd0PnjwIO6//350794ds2bNMntu0qRJOHXqFJ5//nlcvHgRiYmJmDBhgukzrely065du7B79278+OOPcHV1xYABA2o98qo2AgICMHfuXNx3333o3bs3Hn74YWzbtg0dO3bEa6+9hrKyslqdZ9SoUcjOzsbkyZNrHG1GjQODDtmUAQMG4NSpU0hOTsbmzZsRExMDFxcXAFLQOXToELKzs7F582ao1Wr07t0bgDS8dejQoQCAb775Bjt37sS+ffswbdo0ANIXwM2kpqbin3/+qRYOXF1dIYRAenq62fEBAQHVzlH5Pnq9vv4fQoXKc9ys9sDAQOzbt89sGzFiBDp16lRtf0xMzA3PZTQasWDBAgQGBiImJgZZWVnIysrC4MGDYTAYMG/evGqv8fLyMnus1WpvuL+oqMhsv7Ozc7XPS6vVVnt95f6rX5+RkVHjn0NgYKDp+at5e3ubPdbpdABu/BknJCSgT58+uHjxIj799FNs374d+/btQ2xsbLXXFhYW3vTPfv78+XjmmWfw9NNP46OPPrrucXq9/qZ/9qdOnULz5s1hMBgASGE9Ly8Py5cvx9ChQ/HAAw9g/vz5Zn2J/P39cfny5Rue92qHDh3CkCFD0LJlS6xdu9b0mVV64okn8P777+OHH35AcHAwQkJCEBcXh5dffhkAzEJwpU6dOqF79+545JFHsHnzZggh8MYbb9S6pvrQaDR44IEHkJGRgdOnT9fqNdOnT8dbb72FxYsX49FHH2XYaaRu3vuLyIoGDBiA2bNnY8uWLdiyZYvZ//QrQ822bdtMnZQrQ9BPP/0EjUaD1atXm33R1GUuGh8fHzg5OeG777677vNXq2kOj8pjMjMza/2+11N5jmvf91parRadO3c22+ft7Y3c3Nxq+29m48aNpv+JXxsKAKlzcFxcHKKioup0Xkvx9vausbNrZQvGzT672li5ciXy8/Pxyy+/IDQ01LT/8OHD1Y718fG54Z/9/PnzTR28v/rqq+vOAwNIf/43m8+ntLTU7O/7uXPnEBERYfp3AQBdunQxe01SUlKtP5dDhw5h8ODBCA0NxR9//AF3d/caj/vXv/6FyZMn4/Tp03B1dUVoaCieeeYZGAyGm4ZrV1dXtG7duk7D3+tLVAxMqEun5hkzZkChUGDGjBkwGo1YtGhRrTpOk+3gnxbZlL59+0KlUmHZsmU4duwYPvzwQ9Nz7u7uphE158+fx8MPP2x6TqFQQK1WQ6VSmfZVzutxLZ1OV+P/lEeMGIH//Oc/8Pb2RlhYWL3qDw0NhZOTU61Hy9zI2bNnAcCqoWLevHlQKpX45Zdfqn2pJSUlYcyYMfjuu+/w3//+12o13cigQYOwYsUKJCcnm1pxAGDhwoVwdnZG9+7db/k9KsPI1S0ZQgh888031Y5t3br1dcP1ggULMH78eDz66KP49ttvbxhyysrKkJiYaBb0axISEoL4+HiUl5dDpVLB398fiYmJpseAFH6urWPYsGE3PC8gBbnBgwcjODgYGzZsqHaZ9lo6nQ7R0dEApFawpUuX4qmnnoKTk9MNX5eeno4jR46YLlFbSmlpKZYuXQofHx+0aNGiTq995513oFQq8fbbb0MIgcWLFzPsNCL8kyKbUjksfOXKlVAqldV++fXr1w+ffPIJAPP+OXfccQdmz56Nhx9+GE8//TQyMjLw3//+t1ozOwC0bdsWP/30E5YuXYrw8HDo9Xq0bdsWkydPxvLly9G3b19MmTIF7dq1g9FoREJCAv744w+89NJL6Nat2w3r12q16NGjB/bs2VPtuYKCAqxduxYATM9v3boV6enpMBgMGD58uNnxe/bsgbe3t9WGs2ZkZODXX3/FsGHDrttv5OOPP8bChQsxa9asGocfW9vbb7+N1atXY8CAAXjrrbfg5eWFRYsWYc2aNfjwww+v2wJRF0OGDIFWq8VDDz2EV199FUVFRZgzZw6uXLlS7dj+/fvju+++w6lTpxAREWHa//PPP+PJJ59Ehw4d8Mwzz2Dv3r1mr+vYsaPZ39V//vkHBQUFN50Ms7LFbsOGDbjtttswfPhwTJo0CS+99BLeeOMN5ObmYuLEiQCkvl4vv/wy1q1bV22Y/7VOnjyJwYMHAwBmzpyJ06dPm13uad68uWlY+NGjR7F8+XJ07twZOp0Of//9N95//320bNkS//73v02vyc7OxpAhQ/Dwww+jZcuWcHJywqlTp/Dpp5+iuLgYb7/9tlkNarUa/fr1w6ZNm0z7Lly4gH379gGA6T8Ty5YtAyBNdVD5eUydOhWlpaXo1asXmjRpgsTERHz++ec4fPgw5s+fb/Yfotp66623oFQqMX36dAghsGTJEoadxkLOntBENXn11VcFANG5c+dqz1WO3NFqtWYjeoQQ4rvvvhOtWrUSOp1OhIeHi1mzZol58+YJAGZDb8+fPy+GDh0qXF1dBQCzETZ5eXnizTffFK1atRJarVa4u7uLtm3biilTpohLly6ZjgMgJk6cWGP98+bNEyqVqtpw1MohrzVt147yMRqNIjQ0VEyaNKmWn5q5+gwv/+STT244CkgIIb766iuzUVKhoaHijjvuqHZcTZ9PTZO4jR07VhgMhmqv79evn2jTpk21/TW935EjR8TIkSOFu7u70Gq1on379mL+/Plmx1SOuvr5559rrOnq42sadfXbb7+J9u3bC71eL4KCgsQrr7wi1q1bV23kXHZ2tnBxcREffvih2esrR5ddb7t2aPj06dOFj4+PKCoqqvYZXOudd94R0dHRIi8vz1Srp6enACBUKpV47bXXRGhoqFAqlWLo0KHi5MmTNz3n1aOOatqu/rxOnjwp+vbtK7y8vIRWqxUtWrQQb775pqmeSkVFRWL8+PEiMjJSuLi4CLVaLYKDg8Wjjz4qjh07Vq0GANX+Dt+orrFjx5qOmzdvnujatavw8vISarVaeHp6imHDhtU41P1GP/+1IzCFEGLmzJkCgLj77rtrnFSRbI9CiKtmUyOiW1ZUVISQkBC89NJL+Ne//lWvc2zatAlDhw7FsWPH0Lp16waukCxp0qRJ2LRpE44dO3bDy1PXU15ejhYtWuDhhx/GzJkzb3p8YWEh+vfvD1dXV/zvf/+Dl5cXysrKcPr0afj7+8PLywunTp2Cv79/g7RwETU2HHVF1MD0ej1mzJiB2bNnV5tRubbee+89PPHEEww5jdCbb76JixcvYvny5fV6/Y8//oi8vDy88sortTreyckJa9asgdFoRGRkJP7zn/8gLi4OwcHB0Gg0+Pvvv/Hzzz8jJibG7DIQkaPgBUYiC3j66aeRlZWFs2fP1rmPzZUrV9CvXz8899xzFqqOLMnf3x+LFi2qsQ9PbVSO7PHw8Kj1a3x8fLBx40b8+OOP+Pzzz/Hmm2+aRhhVTsMwe/bsGy4mS2SveOmKiMjOZGdnm4bYh4aGmk3cSORoGHSIiIjIbrGPDhEREdktBh0iIiKyWw7dGdloNCI5ORmurq71GgZKRERE1ieEQG5uLgIDA2+6pIdDB53k5ORqqx4TERFR45CYmIjg4OAbHuPQQcfV1RWA9EG5ubnJXA0RERHVRk5ODpo2bWr6Hr+RRh90cnNzMXDgQJSWlqK8vBwvvPACnnrqqVq9tvJylZubG4MOERFRI1ObbieNPug4Oztj69atcHZ2RkFBAaKjo3H33XfD29tb7tKIiIhIZo1+1JVKpTJNhlVUVITy8nJwaiAiIiICbCDobNu2DSNHjkRgYCAUCgVWrlxZ7Zgvv/wSYWFh0Ov1iImJwfbt282ez8rKQvv27REcHIxXX30VPj4+VqqeiIiIbJnsQSc/Px/t27fHF198UePzS5cuxeTJkzFt2jQcOnQIffr0wfDhw5GQkGA6xsPDA3///TfOnTuHxYsXIzU11VrlExERkQ2zqSUgFAoFVqxYgdGjR5v2devWDZ06dcKcOXNM+yIjIzF69GjMmjWr2jmeffZZDBw4EPfdd1+154qLi1FcXGx6XNlrOzs7m52RiYiIGomcnBy4u7vX6vtb9hadGykpKcGBAwcwdOhQs/1Dhw7Frl27AACpqanIyckBIP3g27ZtQ6tWrWo836xZs+Du7m7aOIcOERGRfbPpoJOeno7y8nL4+/ub7ff398elS5cAAElJSejbty/at2+P3r174/nnn0e7du1qPN/rr7+O7Oxs05aYmGjxn4GIiIjk0yiGl187Tl4IYdoXExODw4cP1+o8Op0OOp2uocsjIiIiG2XTLTo+Pj5QqVSm1ptKaWlp1Vp5iIiIiK5l00FHq9UiJiYGGzZsMNu/YcMG9OzZU6aqiIiIqLGQ/dJVXl4ezpw5Y3p87tw5HD58GF5eXggJCcHUqVMxZswYdO7cGT169MDcuXORkJCACRMmyFg1ERERNQayB539+/djwIABpsdTp04FAIwdOxYLFizAAw88gIyMDLz77rtISUlBdHQ01q5di9DQ0Hq/Z2xsLGJjY1FeXn7L9RMREZHtsql5dKytLuPw68pYWoK8kjK4GZwb9LxERESOri7f37K36NgjYSxH6qxoqMuLEKP8Gk193BDq7YxQL2eEeBtM931ddbVaeZWIiIjqh0HHAooL8xFgTAUUgHNRCg4nGnE4MavacXqNEiFezgjxckZTL2c09ay49XJCU09nGHT84yEiIroV/Ca1BKXKdHfO/a2RqA7D+YwCJGTm40JGAS5kFCAluxBFpUacSs3DqdS8Gk/j6awxBaBgTycEezohyNMJwZ7OCPJwYhAiIiK6CX5TWli4uwLR4QHV9peUGZGcVYiEzAIkZBYgMbMAiVcKkJhZiMQrBcgqKMWVglJcKcjGP0nZNZ7bw1kjhR8PJwR6SLdBHlIYCvRwgrdBy0tjRETk0Bh0LExRWljjfq1aiWY+BjTzMdT4fE5RKZIqQs/FK4VIulKIpCsFSLpSiItZhcguLEVWgbQdvZhz3fcIdNcjsCIIBbrrEeDhhICKfQHuerjqNQ32sxIREdkahww6Vh1eXlZUr5e56TWICtQgKrDm3uS5RaW4mFWIixXBp/J+csX9tNxilJQZcT6jAOczCq77Pq46NQI89Ahwd0Kghx5N3JwqHktbE3cnuPASGRERNVIcXm6B4eVFhfnQfxAIACi8az6c2t/dYOeurZIyI1JzipCcVYjk7EIkZxXhYlYhUrIKkZIt7c8pKqvVuVz1alPoCXDTo4kpBFXcd3OCm5Oal8mIiMgqOLzchlzv0pWladXKihFc15/HJ7+4DCkVISglWwpAl7KLkJJd9Ti3qKxiu36naQBw0qgQ4K6Hv1tVCKp67IQm7np4G7RQKhmGiIjIehh0LK28WO4KrsugU6OFnyta+Lle95i84rKK8CMFn5SsIlzKKcKl7EJcyinGpexCXCkoRWFpOc6m5+Nsev51z6VRKeDnWhGA3PWm1qEmV10m83PVQaOy6SXYiIioEWHQsTBFPfvo2AoXnRot/FzQws/luscUlZYjNafIrEVIelyIS9lSMErLLUZpuTD1J7oehQLwddFd1SrkVO2+n5sOOrXquucgIiKqxKBjaeUlcldgcXqNCqHeBoR61zyCDABKy424nFt8VQgybxWq3F9aLpCWW4y03GL8fZ1h9QDg46JDoKnTtFPF/YrRZR56+LnqoeJlMiIih8egY2mljbtFp6FoVErTMPfrMRoFMvJLTJfKLl3VSpScVfW4pMyI9LxipOcVX3eOIZVSAX9XHQI8qsJP8FVzDAV5OHFoPRGRA3DIoGPN4eUKG+6jY2uUSgV8XXXwddWhbbB7jccIIZCZX1LRYbpqVFlKRYfq5Io+ROVGgeTsIiRnF+HAhSs1nstNr0aghzTjdOBVky0GeTihqZczJ1wkIrIDHF5u4eHlJd2eh3b4zAY7N91cuVHgcm6xKQBVzi109bxD2YWlNz2Pk0ZlWnesqZe0DEfVmmRsESIikguHl9uSMrboWJtKqTCN5kJIzcfkFZdJAeiqCRcrHyddKURqbhEKS8tvuBaZh7PGFHpCvAxo5u2MUG8Dmvk4w99Vz6H0REQ2gEHHwnjpyja56NSI8HdFhH/NQ+uLy8px8UohEq8UmtYhS8qUluFIvFKIzPySiiU4snHkYvV+Qjq1EqHezlUByEe6beZtQKCHEztKExFZCYOOpTnAqCt7pFOrEO7rgnDfmofV5xWXSaEnUwpCCZkFOJ8hrU6fmFmA4rLrr0yvUSnQ1NMZod7OaO7rguZ+LtKtrwFe7BdERNSgGHQsQKdzQkH7cXD+ewHUavbjsEcuOjVaN3FD6ybVrw2XlRuRnFVUEXzycT6jwHSbkFGAknKjaXLFzScvm73Ww1ljCj3SrRSEmno6Qc2JFImI6oxBxwIUSiWc/ZpL98tv3umV7ItapUSItzNCvJ0B+Jo9V24UuJRThAsVQefs5XzEX85D/OU8XMwqRFZBKQ5cuFJtpJhGpUAz74rw42dASz9XtGriiua+LtCqGYCIiK6HQcdS1DrptpHPjEwNS6VUSMPYPZzQs4WP2XOFJeU4l14VfOIv5yM+LQ9n0/NQVGrE6bQ8nE7LA45VvUatVCDc14BWTdzQuokrWvlLASjY04mXwIiIwKBjOSqtdMs+OlRLTloVogLdEBVofjnMaBRIzi40BZ8zl/NwOjUXJy7lIreozNQX6Le/q17jolOjVRMp9FQGoNZN3ODuzEupRORYHDLoWGXCQLVeuuXwcrpFSqUCwZ7OCPZ0Rr+IqkthQgikZBfh5CUp9Jy4lIOTl3IRfzkPecVlNV4Ca+KmR+sAV7QJdEPbIA+0DXZHoLuerT9EZLc4YaAFJgwEABxdDix7AgjtDTy+pmHPTXQDpeVGnL2cbwo+lUHoeoupehm0aBvkjrZB7ogOcmf4ISKbxwkDbUFliw7n0SEr06iUpstWV8stKsWp1FzEpeTiaJI0/8+p1Fxk5pdg66nL2HqqagSYt0ErhZ6K4NM2yB0BDD9E1Agx6FiKqrIzMoMO2QZXvQYxoV6ICfUy7SsqLceJS7k4cjEbR5KycORiDk6l5iLjBuGnXbA72gV7oFOIB7xddHL8KEREtcagYylqBh2yfXqNCh2aeqBDUw8AoQCk8HM8JQdHL0qtPv8kZeN0Wl6N4Sfcx4CYUE90buaJmFAvNPc1sNWHiGwKg46l8NIVNVJ6jQodQzzRMcTTtO/q8PNPUjYOJWbhTFqeaeLDnw8kAQA8nTWICZVCT+dmnmgb5A69RiXXj0JExKBjMWzRITtSU/i5kl+CgwlXsP/CFRw4fwV/J2XhSkEpNh5Pw8bjaQAArUqJ6CA3s/Djw8tdRGRFDDqWYhpezgkDyT55GrQYFOmPQZH+AICSMiOOJWfjwIUr2H9eCkDpecU4mJCFgwlZ+Gb7OQBAM29nU+jpGuaFcB9e7iIiy2HQsZTKFp1SBh1yDFq10tTqM76PNM9PQmaBKfQcuJCJU6l5OJ9RgPMZBVh+ULrcFeCuR68WPujdwgc9W3jDz1Uv809CRPaEQcdSrm7REQLg/1jJwSgUCoR6GxDqbcA9McEAgOyCUhxMkCYy3Hc+E4cSs5CSXYRlB5KwrKKfTyt/V/Rq4YM+LX3QNcwLBh1/TRFR/fE3iKVUtuhASMtAqNkvgcjdWYMBrf0woLUfAGl9r/0XMrHjTDp2nknHseQcnEzNxcnUXHy38xzUSgU6hXhKLT4tvdEu2AMaruJORHXgkEHHKktAaJyq7pcVMegQ1cBJq0Kflr7o01Ja2iIzvwS74qXQs/10OpKuFGLv+UzsPZ+JjzdKa3h1D/dG7xbe6N3SB819Xdi/h4huiEtAWGoJCCGAGZ4ABPDSKcDVv2HPT+QAEjIKsONMOnacuYxd8RnIKig1e97fTWfq39O7hQ/83Ni/h8gR1OX7m0HHUkEHAN5rApQVAi/+DXg2a/jzEzmQcqNAXHKO6TLX3vOZKCkzmh3TLtgdQyL9MbRNE0T4s7WHyF4x6NSSxYPOB82AwivAxL2Ab6uGPz+RAysqLceBC1ew/bQUfI4mZ+Pq32YhXs4YEuWPoVH+iAn1hJp9e4jsBhf1tBUaZynolBbIXQmR3dFrVOjVwge9WvgAANJyi7DpeBo2xKVix5l0JGQWYN6Oc5i34xw8nTUY2NofQ9v4o29LXzhpOVszkaNg0LGkyg7JpYXy1kHkAPxc9Xioawge6hqC/OIybDt1GRviUrHpRBquFJRi+cEkLD+YBJ1aiT4tfTE0yh8DI/04UzORnWPQsSQ1gw6RHAw6NYa3DcDwtgEoKzdi7/lMbIhLxYa4VCRdKcTG46nYeDwVCgXQOdQTQ6L8MSSqCcJ8DHKXTkQNjH10LNlH59shQNJe4IFFQOSIhj8/EdWJEALHU3Kl0HP8Eo5ezDF7vqWfi9Svp00TtAtyh1LJzsxEtoh9dGyFpmKoK1t0iGyCQqFAVKAbogLd8OLglriYVYiNFS09e85m4HRaHk6n5eHLLfHwc9VhcJQ/RrQLQPcwb4YeokaKQceSNM7SLTsjE9mkIA8njO3ZDGN7NkN2YSm2nEzDH3Gp2HryMtJyi7H4rwQs/isBTdz0GNUhEHd2CERUgBuHrRM1Igw6lsTOyESNhruTBnd2CMKdHYJQXFaO3fEZWHfkEtYeTcGlnCLM3XYWc7edRUs/F4zuGIRR7QPR1MtZ7rKJ6CYYdCxJU9GxsTRf3jqIqE50ahX6t/JD/1Z+eHd0G2w+cRm/Hr6ITSfScDotDx+tP4mP1p9ETKgnRncIxB3tAuFl0MpdNhHVwCGDjlXWugIAbeWlK7boEDVWOrUKt0U3wW3RTZBdWIr1Ry9h5eGL2H02AwcuSCuxz/gtDn0jfHFnh0AMifKHs9Yhf7US2SSOurLkqKsNbwE7PwW6TwRu+0/Dn5+IZHMpuwi//Z2MlYcv4lhy1egtZ60Kw9o0wagOgejTwoczMhNZAEdd2QrTpSt2RiayN03c9Xiqbzie6huOM2m5+PWwFHoSMwux4tBFrDh0Ed4GLUa2D8SDXZuidRML/GeKiG6KQceSKi9dlbCPDpE9a+HnipeGtsLUIRE4mJCFXw9fxOp/UpCRX4IFu85jwa7z6BjigYe7hmBEu0AuQUFkRQw6lqRliw6RI1EoFIgJ9URMqCemj4jCjtPp+N/+RGyIS8WhhCwcSsjCu6vjcFfHIDzcLYStPERWwKBjSVoX6bYkT946iMjqNColBrT2w4DWfkjLLcLP+5Pw074EJGYWYuHuC1i4+wI6hnjgoa4hGNEugB2YiSyE/7IsyRR0eOmKyJH5ueoxcUALPNuvOXbGp2PJ3gT8cayqleffv8Xhrk5BeKhrCCID2MpD1JAYdCyp8tJVMVt0iAhQKhXo09IXfVr64nJuMZYdSMKSvQlIyCwwtfJ0aFrRl6c9W3mIGgKHl1tyeHnSAeDbgYB7CDDlSMOfn4gaPaNRYFd8BpbsTcD6Y5dQZpR+Jbvq1BjdUWrliQpkKw/R1Ti83FZUtuiwjw4RXYdSqUDvlj7o3dIHl3OLsfyg1MpzIaMAP+y5gB/2SK08T/QOw/DoJtBwXh6iOmGLjiVbdLKTgI/bACotMP1yw5+fiOyS0Siw+2wGFu9NwB/HLqG0XPo1HeCux7iezfBg1xC4O2lkrpJIPnX5/mbQsWTQKcoG3g+R7r+ZBqh1Df8eRGTX0vOKsWhPAn7Ycx7peSUApNmX74sJxuO9wtDMxyBzhUTWx6BTSxYPOsZy4F0v6f4rZwGDd8O/BxE5hKLScqz6Oxnf7TiHE5dyAQAKBTA40h/je4eha5gXFAqFzFUSWQf76NgKpUpaBqI0HyjOYdAhonrTa1S4v3NT3BcTjJ1nMjBvx1lsPnkZG+JSsSEuFdFBbniydxjuaBsIrZr9eIgqsUXHki06APDfVkDeJeCZ7UBAO8u8BxE5pDNpefhu5zn8cjAJRaVGAIC/mw6P9WiGR7qFwMNZK3OFRJZRl+9vxn5L07lKtxx5RUQNrIWfC/5zV1vsfm0QXhnWCn6uOqTmFOOj9SfRY9afePe3OFzKLpK7TCJZMehYWmXQKcqWtw4islueBi0mDmiBHf8aiNn3t0dUgBsKS8vx3c5z6PvhZkxbcQSJmVxzjxyTQwad2NhYREVFoUuXLpZ/M31Fk1pRjuXfi4gcmlatxN2dgrHmhd5Y+ERXdGnmiZJyIxb9lYAB/92CV37+G+fSuSQNORb20bF0H53/PQbE/Qrc/l+g61OWeQ8iouv462wGPv/zDHacSQcAKBXAyPaBmDigBSL8XWWujqh+OOrKlugqW3SyZC2DiBxTt3BvdAv3xsGEK4j98ww2nUjDr4eT8evhZNzWpgmeH9gC0UHucpdJZDEOeenKqvQVv0B46YqIZNQpxBPzxnXB6km9MTy6CQDg92OXMOLzHXhiwT4cTLgic4VElsEWHUurDDrFDDpEJL/oIHfMeTQGp1Jz8eXmM1j1dzL+PJGGP0+koVcLb0wa2BLdwznnF9kPtuhYWmXQKcyStQwioqtF+Lvikwc7YtNL/XF/52ColQrsPJOBB+fuwZh5f+FIEkeKkn1g0LE006Ur/tIgItsT5mPAh/e2x5ZX+uPR7iHQqBTYfjodI7/YgYmLDiL+MucAo8aNQcfSnDyl20Je/yYi2xXs6Yz3RrfFpqn9cVfHICgUwJojKRj68Ta8tvwfpGQXyl0iUb0w6Fia3kO65agrImoEQryd8fEDHbDuxT4YHOmHcqPAT/sS0e+jLZi5Jg5X8kvkLpGoThh0LI0tOkTUCLVu4oZvx3bB8md7oGuYF0rKjPhmuzTT8mebTiO/uEzuEolqhUHH0pw8pNuiHMBYLmspRER1FRPqhaVPd8eCx7sgKsANucVlmL3hFPp+uBnzd55DcRl/r5FtY9CxtMpLVxDskExEjZJCoUD/Vn5YPak3Pn+oI8J8DMjIL8GM3+IwZPY2/H70Ehx4kn2ycQw6lqbWAtqKadZ5+YqIGjGlUoGR7QPxx5S++M9dbeHnqkNCZgEm/HgAD3/zF46ncL4wsj0MOtbgXNFPpyBD3jqIiBqARqXEw91CsPnl/pg0sAV0aiV2n83AHZ9tx+u/HEFGXrHcJRKZMOhYg3PFLKMFmfLWQUTUgAw6NV4a2gqbXuqHO9oFwCiAJXsT0P+jLfhm21mUlBnlLpGIQccqTEGHLTpEZH+CPZ0R+3An/O+ZHogOkjosz1x7HMM+2YaNcansv0OyYtCxBicv6ZZBh4jsWNcwL/w6sTc+vKcdfFx0OJeej/EL9+Ox7/biVGqu3OWRg2LQsQa26BCRg1ApFbi/S1NsfrkfJvRrDq1Kie2n0zH80+1469ejyC4olbtEcjAMOtbAoENEDsZVr8Frw1tjw9S+GNbGH+VGgYW7L2DQ7K349fBFXs4iq2HQsQaDj3Sbny5vHUREVhbqbcDXYzpj8fhuaO5rQHpeMV786TAe+24vLmTky10eOQAGHWsw+Eq3+Wny1kFEJJOeLXyw9sU+mDokAlq1dDlr6MfbELv5DEdnkUUx6FiDi590m39Z3jqIiGSkU6vwwqCW+P3FPujZ3BvFZUZ8tP4k7vhsO/ad5/QbZBkMOtZQ2aKTdxngdWkicnDhvi5YNL4bPn6gPbwNWpxOy8N9X+3Ga8v/QVYBV0enhuWQQSc2NhZRUVHo0qWLdd6wMuiUFQIlvCZNRKRQKHBXx2BseqkfHuzSFADw075EDPq/rVhxKImdlanBKIQD/23KycmBu7s7srOz4ebmZtk3mxkAlBYALxwCvMIt+15ERI3M3nOZmLbiCE6n5QEAerXwxszRbdHMxyBzZWSL6vL97ZAtOrKoHHmVx346RETX6hrmhTUv9MErw1pBp1Zi55kM3PbpNizYeQ5Go8P+f5waAIOOtbg0kW7zLslbBxGRjdKqlZg4oAX+mNIXPcK9UVRqxDu/xeHhb/cgMbNA7vKokWLQsRZXf+k2N1XeOoiIbFyotwGLxnfDu3e2gZNGhT1nMzHsk234cc8F9t2hOmPQsRbXAOmWLTpERDelVCrwWI9m+H1yH3Rt5oWCknK8ufIoxszbi4tZhXKXR40Ig461uFS26DDoEBHVVqi3AT893R3TR0RBp1Zix5l0DPt4G5buS2DrDtUKg461uFb00WHQISKqE6VSgSd7h2Hdi33QKcQDecVl+NfyI3h8wT5cyi6SuzyycQw61sKgQ0R0S8J9XfDzhJ544/bW0KqV2HLyMm77dBvWH+PvVbo+Bh1rcQ2UbnMuylsHEVEjplIq8HTf5lj7Qm+0DXJHVkEpnvnhAN5ceQRFpeVyl0c2iEHHWtyDpNuiLM6OTER0i1r4uWL5sz3xdF9pAtYf9yRg1Bc7cPJSrsyVka1h0LEWvTugdZXu5yTLWwsRkR3QqpV44/ZILHyiK3xcdDiVmodRX+zAD7vPs6MymTDoWFNlq052krx1EBHZkb4Rvvh9ch/0i/BFcZkR0389hqd/OIAr+VwglBh0rMutsp8OW3SIiBqSj4sO88d1wZt3REKjUmBDXCqGf7odu+Mz5C6NZMagY01uFS067JBMRNTglEoFxvcJx4rneiHcx4BLOUV4+Ns9+Gj9CZSVG+Uuj2TCoGNN7sHSbXaivHUQEdmx6CB3/DapN+7vHAwhgNjN8Xh03l+4nFssd2kkAwYda/IIkW6zGHSIiCzJoFPjw3vb47OHOsJZK62XNfLzHThw4YrcpZGVMehYkynoXJC3DiIiBzGqfSBWPd8LzX2lS1kPzt2NhRyV5VAYdKzp6hYdI68XExFZQws/V/z6fG/c3rYJSssF3vr1GKb+728UlnCCQUfAoGNNroGAUg0YS4HcFLmrISJyGC46NWIf7oRpt0dCpVRgxaGLuOvLnTifzglc7R2DjjWp1FUjr7IS5K2FiMjBKBQKPNU3HD8+2Q0+LlqcuJSLkV/swMa4VLlLIwti0LE2z1Dplv10iIhk0aO5N1ZP6oOYUE/kFpVh/ML9+HTjafbbsVMMOtbmURF0Ms/JWwcRkQNr4q7Hkqe6Y1zPZgCAjzeewos/HebCoHaIQcfavJtLt5nx8tZBROTgtGol3hnVBu/f3RZqpQKr/k7Gg3P3IC23SO7SqAEx6FibV0XQyWDQISKyBQ92DcHCJ7vC3UmDw4lZGP3FThxPyZG7LGogDDrWZmrROStvHUREZNKzuQ9WTpSWjkjOLsI9c3axk7KdYNCxNs8w6bYoCyjIlLUUIiKqEuZjwIrneqFnc28UlJTjqR/245ttZ9lJuZFj0LE2rbM0nw7Ay1dERDbG3VmD75/oioe7hUAIYOba43ht+RGUlHGS18aKQUcO7JBMRGSzNColZo6OxlsjoqBUAEv3J2Lc/L3ILSqVuzSqBwYdOXi3kG7TT8tbBxER1UihUOCJ3mGYN7YLDFoVdsVn4MG5e7gCeiPEoCMH31bS7eUT8tZBREQ3NKC1H356uge8DVocS87BvV/twoUMLhvRmDDoyMEnQrpNPyVvHUREdFNtg92x7NmeaOrlhAsZBbhnzm4cvZgtd1lUSww6cvBtLd1mxANlJfLWQkRENxXmY8DyCT0RGeCG9LxiPDh3D3bFp8tdFtVCow86iYmJ6N+/P6KiotCuXTv8/PPPcpd0c26BgNYVEOWcT4eIqJHwc9Nj6TPd0S3MC3nFZRj33T6sPZIid1l0E40+6KjVanzyySeIi4vDxo0bMWXKFOTn2/j1U4UC8Gkp3U8/KW8tRERUa256afj5bW2aoKTciImLD+KHPVyk2ZY1+qATEBCADh06AAD8/Pzg5eWFzMxGMBFf5eWrtOPy1kFERHWi16gQ+0gnPFIx1870lUfx7Xa2ztsq2YPOtm3bMHLkSAQGBkKhUGDlypXVjvnyyy8RFhYGvV6PmJgYbN++vcZz7d+/H0ajEU2bNrVw1Q3AP0q6TT0qbx1ERFRnKqUC742OxsQB0rxo7605jjlbODeaLZI96OTn56N9+/b44osvanx+6dKlmDx5MqZNm4ZDhw6hT58+GD58OBISEsyOy8jIwGOPPYa5c+dao+xb5x8t3V5i0CEiaowUCgVeHtoKkwdLXRE++P0EPt/E+dFsjULY0CIeCoUCK1aswOjRo037unXrhk6dOmHOnDmmfZGRkRg9ejRmzZoFACguLsaQIUPw1FNPYcyYMdc9f3FxMYqLqyZ7ysnJQdOmTZGdnQ03N7eG/4FuJD8d+KhihuTXLwI6F+u+PxERNZgv/jyN//4hTRny4qCWmDy4JRQKhcxV2a+cnBy4u7vX6vtb9hadGykpKcGBAwcwdOhQs/1Dhw7Frl27AABCCIwbNw4DBw68YcgBgFmzZsHd3d20yXqJy+ADuDSR7qfFyVcHERHdsucHtsRrw6W+l59uOo3//nGSi4HaCJsOOunp6SgvL4e/v7/Zfn9/f1y6dAkAsHPnTixduhQrV65Ehw4d0KFDBxw5cqTG873++uvIzs42bYmJiRb/GW6oSeXlq5rrJSKixmNCv+Z4845IAEDs5ni8v+4Ew44NUMtdQG1c2/wnhDDt6927N4zG2q0qq9PpoNPpGry+evOPBs5sZIdkIiI7Mb5PODQqJd5edQxfbzuL0nKB6SMieRlLRjbdouPj4wOVSmVqvamUlpZWrZWnUWrSVrpliw4Rkd0Y27MZZt4ltdh/t/McPlzP+dLkZNNBR6vVIiYmBhs2bDDbv2HDBvTs2VOmqhpQYEfp9tIRoLxU3lqIiKjBPNItFP+5S/rP7Jwt8fhyyxmZK3Jcsl+6ysvLw5kzVX8Bzp07h8OHD8PLywshISGYOnUqxowZg86dO6NHjx6YO3cuEhISMGHCBBmrbiBe4YDeHSjKljokB7SXuyIiImogD3cLQV5xKf6z9gQ+/P0kXHVqjOnRTO6yHI7sQWf//v0YMGCA6fHUqVMBAGPHjsWCBQvwwAMPICMjA++++y5SUlIQHR2NtWvXIjQ0tN7vGRsbi9jYWJSXl99y/bdEoZBadc5uAS4eZNAhIrIzT/dtjtyiMnz+5xlM//UYXPRq3NUxWO6yHIpNzaNjbXUZh28xG2cAO2YDnR4DRn0uTw1ERGQxQgjM+C0OC3adh0qpwJePdMKwNk3kLqtRs5t5dBxCUCfp9uJBeesgIiKLUCgUeGtEFO7pFIxyo8CkxYew80y63GU5DAYduQVWBJ2040BJgby1EBGRRSiVCnxwT1vTqudPLdyPw4lZcpflEBh05OYWKM2QLMqB5ENyV0NERBaiVinx6UMd0KelDwpKyvHkgn1IzOR/cC2NQUduCgUQ0l26n7Bb3lqIiMiidGoVvno0BlEBbsjIL8G4+XuRXcDpRSyJQccWhPSQbhP2yFsHERFZnEGnxnfjuiDAXY/4y/l45sf9KCmr3Qz/VHcOGXRiY2MRFRWFLl26yF2KpLJFJ3EvYJR5yDsREVlcE3c9vhvXBS46NfaczcRry//hulgW4pBBZ+LEiYiLi8O+ffvkLkXiHw1oXYDibKlTMhER2b3IADfEPtIJKqUCvxy6iE82npa7JLvkkEHH5qjUQHBF6xL76RAROYx+Eb74953SulifbjqNZQeSZK7I/jDo2Ar20yEickgPdwvBhH7NAQCv//IPdsVzjp2GxKBjKyr76VzYCfA6LRGRQ3l1WCvc0S4ApeUCExcd5LDzBsSgYyuadgPUeiA3BUg/JXc1RERkRUqlAv93X3u0DXLHlYJSPLvoAIpKOTilITDo2AqNvqpV5+wWWUshIiLr02tUmPNoJ3gZtDh6MQfTVhzlSKwGwKBjS8L7S7cMOkREDinY0xlfPNQRSgWw/GASfthzQe6SGj2HDDo2N49Opcqgc247UF4maylERCSPni188Nrw1gCAd3+Lw77zmTJX1Lg5ZNCxuXl0KjVpBzh5AiW5QDJXMyciclRP9QnHiHYBKDMKPLfoIFJziuQuqdFyyKBjs5QqIKyvdD9+s7y1EBGRbBQKBT68tx1a+bvicm4xnv3xAErLuUxEfTDo2JrwAdLtmY3y1kFERLJy1qrx9ZgYuOnVOJiQhY83cERufajr86KkpCSsWrUKCQkJKCkpMXtu9uzZDVKYw4oYJt0m7QPyLgMuvvLWQ0REsmnmY8AH97TDs4sOYs7WePRu4YOeLXzkLqtRqXPQ2bRpE0aNGoWwsDCcPHkS0dHROH/+PIQQ6NSpkyVqdCxugUBAeyDlb+D0eqDjo3JXREREMhreNgAPdW2KJXsTMeV/h7Huxb7wMmjlLqvRqPOlq9dffx0vvfQSjh49Cr1ej+XLlyMxMRH9+vXDfffdZ4kaHU/EcOn25Dp56yAiIpswfUQUmvsakJpTjFeXcaXzuqhz0Dl+/DjGjh0LAFCr1SgsLISLiwveffddfPDBBw1eoENqdZt0G78ZKGVPeyIiR+esVeOzhzpCq1Ji4/FU/Mj5dWqtzkHHYDCguLgYABAYGIj4+HjTc+npjWMhMpudR6dSQAfANQAozQfOb5e7GiIisgFtAt3xr4r5dd5bcxwnL+XKXFHjUOeg0717d+zcuRMAcMcdd+Cll17CzJkz8cQTT6B79+4NXqAl2Ow8OpUUCiCiolWHl6+IiKjCE72aoX8rXxSXGfHCkkNcD6sW6hx0Zs+ejW7dugEA3nnnHQwZMgRLly5FaGgo5s2b1+AFOqxWV/XTMXLuBCIikubX+e997eHjosPJ1Fx8tP6k3CXZvDqNusrNzUV8fDxKS0sRGBgIHx8ffPnll5aqzbGF9QN0bkBuMpD4FxDaQ+6KiIjIBvi46PDRfe3w+Px9+G7nOdzRLgCdQjzlLstm1bpF559//kHr1q1x2223YcSIEWjRogU2buSkdhaj0QOtR0j3jy6XtxYiIrIpA1r54Z5OwRACeHXZP7yEdQO1DjqvvfYaQkJCsH37duzfvx/9+vXD888/b8naKPpu6TZuJRf5JCIiM9NHRMLXVYczaXn4/M/Tcpdjs2oddPbv34/PP/8cPXv2RKdOnfDdd9/h9OnTyMvLs2R9ji28v7TIZ/5l4MIOuashIiIb4uGsxb/vjAYAfLX1LI5ezJa5IttU66CTnp6OkJAQ02Nvb284Ozvj8uXLFimMAKg0QNSd0n1eviIiomvcFt0Ed7QLQLlR4JVl/3DhzxrUOugoFArk5uYiJycHOTk5yM7OrrYvJyfHkrU6puh7pNvjvwFlJTc+loiIHM6MUW3g6azB8ZQcfLUl/uYvcDC1DjpCCERERMDT0xOenp7w8vJCXl4eOnbsCE9PT3h4eMDTk72+G1xoL8DFHyi8AsRvkrsaIiKyMT4uOrwzqg0A4LM/T+Ncer7MFdmWWg8v37x5syXroOtRqoB29wO7PgcO/lA1vw4REVGFUe0DseLQRWw5eRkzfjuG+eO6QKFQyF2WTah10OnXr99NjyksLLylYqwlNjYWsbGxKC9vJMPxOj4mBZ1TvwO5lwDXJnJXRERENkShUOCtEVHYeWYbtpy8jE3H0zA4yl/usmxCnWdGnjhxYo378/PzMXx442htsPklIK7lGwE07Q6IcuDwYrmrISIiGxTu64LxfcIBAO+ujuPcOhXqHHT++OMPvPnmm2b78vPzcdtttzWeFpLGqNNj0u3BhYAQ8tZCREQ26fkBLdDETY+EzAJ8s+2s3OXYhHoFnfnz5+Pjjz8GIC0LMWTIECgUCvz+++8NXiBVaDMa0LoCV85xRXMiIqqRQafG67dLK5zHbjmDi1mNo0uJJdU56ISFhWH9+vWYOXMmPv30UwwdOhRarRbr1q2DwWCwRI0EAFoD0PZe6f7BhfLWQkRENmtU+0B0DfNCUakR/1lzXO5yZFfnoAMA0dHRWL16NaZNmwZnZ2eGHGupvHwVtwooyJS3FiIiskkKhQIzRrWBUgGsOZKCnWfS5S5JVrUaddWxY8cah6npdDokJyejV69epn0HDx5suOrIXGBHwL8tkHoEOPQD0OtFuSsiIiIbFBnghjHdQ/H97guYueY4Vk/qDaXSMYeb1yrojB492sJlUK0oFEC3Z4BVzwN/zQW6PyctE0FERHSNyYMj8MvBi4hLycHqIykY1T5Q7pJkoRDCcYfw5OTkwN3dHdnZ2XBzc5O7nNopLQI+iZYW+rxnXlW/HSIiomt8vuk0/m/DKYR6O2Pj1H7QqOrVY8Xm1OX72z5+Ykei0QNdnpLu7/6CQ82JiOi6nugdBh8XLS5kFGDpvkS5y5EFg05j1OVJQK0Hkg8BCXvkroaIiGyUQafGpIEtAQCfbTqNwhLHm++OQacxMvgA7R+U7u/+Qt5aiIjIpj3UNQTBnk5Iyy3Ggl3n5S7H6hwy6MTGxiIqKgpdunSRu5T66/6cdHtiDZDJ2S+JiKhmWrUSU4dEAADmbDmD7IJSmSuyrloHnR49euCDDz7A8eONf/KhRrfWVU18WwEthwIQwM7P5K6GiIhs2J0dgtDK3xU5RWX4dodj/ee41kFnwoQJ2Lt3L7p27YqIiAi88sor2L59Oxx40Jb8ek+Rbg/9CFw5L2spRERku1RKBSYPlvrqfL/rPHKLHKdVp9ZBZ+zYsVi+fDnS09PxySefICcnBw888AD8/Pwwbtw4rFixAgUFBZasla4V2hMIHwAYS4GtH8ldDRER2bBhbZog3NeAnKIyLP4rQe5yrKbOfXR0Oh1uv/12fP3110hOTsbq1asRFBSEt956Cz4+PhgxYgR27txpiVqpJgMrVpL/ewmQES9vLUREZLOUSgUm9GsOAPh2xzkUlTrGCKxb7ozcrVs3zJw5E0eOHMGRI0cwaNAgpKSkNERtVBvBnYGWwwBRDmx5X+5qiIjIho3uEIQAdz0u5xbjl4MX5S7HKhp01FXz5s0xZcoU3HsvZ+u1qgFvSLdHfgbSTshbCxER2SytWonxfcIBAF9vi0dZuVHmiizPIYeX253ADkDrEQAEsGWW3NUQEZENe6hrU3g6a3AhowBrj16SuxyLY9CxFwPeAKAA4lYCKf/IXQ0REdkoZ60a43qGAQDmbIm3+9HTDDr2wr8NEH23dH/zf+SthYiIbNrYnqFw1qpwPCUHW05dlrsci7qloCOEsPsk2Kj0fx1QqIBT64D4zXJXQ0RENsrDWYuHuoYAAObvPC9vMRZWr6Azb948REdHQ6/XQ6/XIzo6Gt9++21D10Z15dMS6Fqxsvm6fwHljjMhFBER1c1jPUKhUADbTl3G+fR8ucuxmDoHnenTp+PFF1/EyJEj8fPPP+Pnn3/GyJEjMWXKFLz55puWqJHqov/rgLMPkH4S+OtruashIiIbFeptQL8IXwDAj3suyFyN5ShEHa89+fj44PPPP8dDDz1ktn/JkiWYNGkS0tPTG7RAS8rJyYG7uzuys7Ph5uYmdzkN5+BCYNUkQOsKTDoAuPrLXREREdmgP0+k4okF++GmV+OvNwbDSauSu6Raqcv3d51bdMrLy9G5c+dq+2NiYlBWVlbX05EldHgUCOwElOQCG9+RuxoiIrJR/SL80NTLCTlFZVj1t31OIFjnoPPoo49izpw51fbPnTsXjzzySIMURbdIqQRur1j76u/FQOJeeeshIiKbpFIq8Gi3UADAwt0X7HKA0S11Rh4/fjzGjx+P6OhofPPNN1AqlZg6dappIxkFd5ZadgBg7SuA0THWNCEiorq5v3NT6NRKHEvOwaHELLnLaXDqur7g6NGj6NSpEwAgPl5aRNLX1xe+vr44evSo6TiFQtFAJTa82NhYxMbGorzczr/8B78NHF8FpBwG9s4Fuj8rd0VERGRjPA1ajGwfiGUHkrD4rwR0CvGUu6QGVefOyPbEbjsjX23fPGDNVECtBybskIagExERXWX/+Uzc+9VuGLQq7HtzMJy1dW4HsSqLdkamRqbzE0DzgUBZEbDiGaCcHcaJiMhcTKgnQryckV9Sjj+OpcpdToNi0LF3CgUw6gtA5w5cPADs/FjuioiIyMYoFArc1TEIAPDLIfsafcWg4wjcg4DbP5Tub/mAi34SEVE1d3eSgs6O05eRmlMkczUNh0HHUbR7AGg9AjCWAismAGXFcldEREQ2JNTbgM6hnjAK4NfD9tOqw6DjKBQKYMQn0vIQaceALbPkroiIiGzMXRWtOr8cZNChxsjFFxj5iXR/56dAwl+ylkNERLZlRNtAaFVKnLiUi7jkHLnLaRAMOo4mciTQ7kFAGIGVE4AS+12xloiI6sbdWYNBkX4AgF8OJslcTcNg0HFEwz8AXAOBzLPA2lcBx51KiYiIrnF3p2AAwK9/J6Pc2Pi/Hxh0HJGTB3DXHEChBA7/CByYL3dFRERkI/pF+MJVr8bl3GIcTLgidzm3jEHHUYX3BwZOl+6vfZULfxIREQBAq1ZicKQ/AOD3o5dkrubWMeg4st5TgMhR0pDz/z0G5NrXbJhERFQ/w9o0AQCsP3ap0a9ozqDjyBQKYPSXgG9rIDcF+HksUFYid1VERCSzfhG+0GuUSLpSiGONfPQVg46j07kCDywCdG5Awm7gj2lyV0RERDJz0qrQL8IXAPDHscZ9+YpBhwCfFsBdX0v3984FDi+Rtx4iIpJd5eWr3xl0yC60vh3o9y/p/urJQPJhOashIiKZDWrtD7VSgVOpeTh7OU/ucuqNQYeq9HsNaDkMKCsCljwIXLkgd0VERCQTd2cNejT3BgD8Edd4B6sw6FAVpRK4ey7gGyl1Tv7hLiDvstxVERGRTAa2lmZJ3naq8X4XMOiQOScPYMwvgHsIkBkPLLoHKGrcPe6JiKh++lZ0SN5//goKSspkrqZ+GHSoOrdAYMwKaaXzlL+Bnx4GSovkroqIiKws3MeAIA8nlJQbsedshtzl1AuDDtXMpwXw6DJA6wqc3w4sfxIob5xpnoiI6kehUKBfK6lVZ+vJxnn5ikGHri+wI/DQYkClBU6slkZjNfIZMomIqG76tpSCzrbT6TJXUj8OGXRiY2MRFRWFLl26yF2K7QvrC9wzT1oA9NAPwKYZcldERERW1LOFN1RKBc6l5yMxs0DucurMIYPOxIkTERcXh3379sldSuMQNQoY8bF0f8fHwK4v5K2HiIisxk2vQacQDwDA1kY4+sohgw7VQ8w4YNBb0v0/pgGHFslaDhERWY/p8hWDDtm13lOB7hOl+79OBPbNk7ceIiKyisph5rvjM1BubFx9NRl0qPYUCmDoe0DnJwEIYM1UYPv/sYMyEZGdiw5yh4tOjdziMhxPaVxzqzHoUN0olcAd/wf0eVl6vOldYMN0hh0iIjumUirQKdQTALDvfKbM1dQNgw7VnUIBDJoute4AwK7PgVWTAGO5vHUREZHFdG3GoEOOpuckYNQXVUPPfx4HlBXLXRUREVlA1zBpgc+9565ANKJWfAYdujWdxgD3fS9NKnh8FbD4fqA4T+6qiIiogbULdodWpUR6XjHOZzSe+XQYdOjWRY0CHvkZ0BiAs1uAhXcCBY2raZOIiG5Mr1GhfVN3AMC+c43ndzyDDjWM8P7A2FWAkydwcT8w/3YgJ0XuqoiIqAF1aeYFANjbiPrpMOhQwwnuDDy+DnANAC4fB74bBmSelbsqIiJqIF3CKoIOW3TIYflFAk+sBzzDgKwLwDeDgPg/5a6KiIgaQEyoJ5QKICGzAGk5RXKXUysMOtTwPEOlsBPYESjMBH68B9g+m3PtEBE1cm56DVr6uQIA/k7Klrma2mHQIctw9Qce/x3o+CggjNKq5/8bAxQ1rhk1iYjIXNtgqUPykaQseQupJQYdshyNXppnZ8QngFIDHP8N+HYQcPmU3JUREVE9tasIOv9cZIsOkTSLcufHgSd+B1wDgfRTwDcDpdBDRESNTrtgDwDAkaTsRjFxIIMOWUdwZ+CZrUBob6AkF1j6KLDxHS4bQUTUyLRu4gq1UoGM/BIkZ9t+h2QGHbIeFz/gsV+BHs9Lj3d8LHVUzs+Qty4iIqo1vUaFVk2kDsmNoZ8Ogw5Zl0oNDJsJ3PsdoHEGzm4G5vYHkg/LXRkREdWSqZ9OIxh5xaBD8oi+Bxi/CfAKB7ITgHlDgUOL5K6KiIhqoW2QBwAGHaIb848CntoMRAwHyouBX58DVr3ARUGJiGxcVYtOls13SGbQIXk5eQAPLgYGTAOgAA5+D8zpIS0OSkRENinC3xValRI5RWVIulIodzk3xKBD8lMqgX6vSh2V3UOArARpBfTfJgPFuXJXR0RE19CqlQj3NQAATqXa9u9pBh2yHeH9gOd2AV3GS48PzAe+7AHEb5a3LiIiqqZy5NVJBh2iOtC5Anf8HzD2N8AjFMhOBH4YLfXd4fIRREQ2I8JfCjqnLjHoENVdWF/g2V1A16elxwe/l1p3zmySty4iIgIAtPKvbNGx7QEkDDpku3QuwO0fAePWAJ7NgJwk4Me7gV+fB4psf0gjEZE9q7x0FZ+Wh9Jyo8zVXB+DDtm+Zr2l1p1uE6THh36QWndOb5C3LiIiBxbk4QRnrQol5UZcyMiXu5zrYtChxkFrAIZ/ADy+TppkMOcisOheYOVEoDBL7uqIiByOUqlAy8rLV5ds9/IVgw41LqE9gQk7ge4TASiAwz8CX3YHTq2XuzIiIofT2t/2R14x6FDjo3UGbvsP8MTvgFdzIDcFWHw/sOxJICtR7uqIiBxGRBPbH3nFoEONV0h3YMKOitXQFcDRZcAXnYGNMzgUnYjIClr6uQAATqcx6FjUXXfdBU9PT9x7771yl0LWpnWWVkN/ejMQ2hsoKwJ2zAY+6wjs/QYoL5W7QiIiuxXmI82OnJhZiHKjba55ZRdB54UXXsDChQvlLoPkFNgRGLcaeHAJ4N0SKEgH1r4MzOkJnFwH2Piic0REjVGghxO0KiVKyo1IzrLNNa/sIugMGDAArq6ucpdBclMogNa3A8/tBm7/L+DsDaSfApY8CHw/Ekg+JHeFRER2RaVUIMTbGQBwLt02h5jLHnS2bduGkSNHIjAwEAqFAitXrqx2zJdffomwsDDo9XrExMRg+/bt1i+UGg+VBuj6FPDCIaDXZEClA85vB+b2B355BshOkrtCIiK70cxbunx13kbn0pE96OTn56N9+/b44osvanx+6dKlmDx5MqZNm4ZDhw6hT58+GD58OBISEqxcKTU6endgyAxg0n6g7f3Svn9+Aj6PATa9yw7LREQNIMyHLTo3NHz4cLz33nu4++67a3x+9uzZePLJJzF+/HhERkbik08+QdOmTTFnzpw6v1dxcTFycnLMNnIAHiHAPd8AT20GQntJHZa3/5/UYXnft0B5mdwVEhE1Ws0qOiSfZ9Cpu5KSEhw4cABDhw412z906FDs2rWrzuebNWsW3N3dTVvTpk0bqlRqDII6SetmPbgY8G4hdVhe8xIwpwdw8nd2WCYiqofKkVfnMwpkrqRmNh100tPTUV5eDn9/f7P9/v7+uHTpkunxsGHDcN9992Ht2rUIDg7Gvn37ajzf66+/juzsbNOWmMjJ5RyOQgG0vgN4bg8w/CPAyauiw/IDwMJRQMrfcldIRNSoVA0xL0CZDS7uqZa7gNpQKBRmj4UQZvvWr6/d9P86nQ46na5Ba6NGSqUBuj0NtLtfmndnzxzg3Dbg635AuweAPi8BvhFyV0lEZPP8XfXQa5QoKjUi6Uqh6VKWrbDpFh0fHx+oVCqz1hsASEtLq9bKQ1QvTh7AkHeB5/cD0fcCEFKH5diuwE+PAIk1tw4SEZFEqVTY9Mgrmw46Wq0WMTEx2LBhg9n+DRs2oGfPnjJVRXbJMxS4dx7w1J9Aq9sBCODEamDeYGD+7dKioezDQ0RUoyAPJwDARRucNFD2S1d5eXk4c+aM6fG5c+dw+PBheHl5ISQkBFOnTsWYMWPQuXNn9OjRA3PnzkVCQgImTJggY9Vkt4JigIeWAGkngF2fA/8sBS7slDa/KKDXi0D0PdKlLyIiAgAEeVYEnSsMOtXs378fAwYMMD2eOnUqAGDs2LFYsGABHnjgAWRkZODdd99FSkoKoqOjsXbtWoSGhtb7PWNjYxEbG4vy8vJbrp/slF9rYHQsMOANYM+XwIEFQFocsOIZYNO/gR4TgU6PAToXuSslIpKdLbfoKIRw3Pb4nJwcuLu7Izs7G25ubnKXQ7as8Aqw/ztgz1dAfpq0T+8BdH0a6PYMYPCRtTwiIjmt/icZzy8+hM6hnlj2rOW7ltTl+9um++gQ2QwnT2kk1uQjwIiPAa9woCgL2PYh8HEbaT6ezHNyV0lEJAtbbtFh0CGqC40e6PyENErrvu+lVdPLiqQZlj/vBCx7gnPxEJHDqeyjk5pThFIbm0uHQYeoPpQqoM1oaVmJx1YBzQcBwggcXQ583RdYOBo4u4UjtYjIIfgYdNCqlTAK4FJ2kdzlmGHQIboVCgUQ3g8Y8wvwzHZpLh6FEji7GVh4p7Ri+tFfACM7vhOR/VIqFabLV4lXbGspCAYdooYS0E6ai+eFQ1InZbUTkHIYWPa4tGL6vnlAqe1dvyYiagimfjo2NsTcIYNObGwsoqKi0KVLF7lLIXvk2Qy4/SNgylGg37+kjsxXzgFrpgL/1wpY8zL78RCR3bHVDskcXs7h5WRpJfnAwR+A3bFAdkLV/ibtpLl42t4rhSEiokbsk42n8MnG03ioa1PMurudRd+Lw8uJbInWAHSfALx4GHh0ORA1GlBqgEv/AGtfBv6vNbB8PHB2K2C0rdEKRES15e+mBwCk5hTLXIk52WdGJnIYShXQYrC05WdIy0sc+kGacfnIz9LmEQp0HAN0eBhwD5K7YiKiWvN30wEA0nI56oqIDN5Aj+eAZ3dJC4nGPA7o3ICsC8Dm94BPooEf7wWOrQTKSuSulojopvxc2aJDRNdSKKSFRINigGH/AeJ+lVp5LuwEzmyQNmdvoN2DQKcxgF+k3BUTEdWo8tJVel4xysqNUKtsoy3FNqogIkDrDHR4CHh8LTDpINB7CuDSBCjIAPbEAl92B74ZJC0wWpQjd7VERGa8DVqolAoIAaTn2U5LNIMOkS3ybg4MfgeYcgx4aCnQegSgVAMX9wO/vSgNU1/5HHBhN2dfJiKboFQq4Osi9dNJzbGdfjoOeekqNjYWsbGxKC/nbLVk41RqoNVt0paXBvy9RBqqnnEaOLxI2rxbAh0fBdo/BLj6y10xETkwfzcdLuUU2VTQ4Tw6nEeHGhshgMS/pMBzbAVQmi/tV6iAiGHSqK2WQ6WQRERkRU8t3I8Ncan49+hojOkearH3qcv3N38TEjU2CgUQ0l3ahr8vraV16EcgaS9wcq20ufgDUXcCkSOBkJ4MPURkFZVDzC/bUIsOf/sRNWY6VyBmrLSlnZBGbP39E5CXCuydK23O3kDrO4DIO4GwvoBaK3fVRGSnKoeYp+XazhBzBh0ie+HXGhg2Exj0trR6etwq4OQaadTWwYXSpnOX+vtEjgSaD5JGehERNRBvF+k/Uhn5tjPqikGHyN6otVJfnYhhQPkn0pw8cauAE6ullp5/lkqbxhloOQSIHCX16dGznxoR3RpvQ0XQyWOLDhFZg0oDhPeXttv/K/XjiVsFHP9NWmA07ldpU2mB5gOl0NNqOODsJXflRNQIeVcML2eLDhFZn1JZ1Yl52Ewg5XBF6FkFZJwBTv0ubQoVENZHCj2tR3DIOhHVWmWLTqYNTRjokEGH8+iQw1MogMCO0jboLeDyiaqWntQjwNkt0rbmJSkYRY4CIkcAHiFyV05ENszbILXo5BaXoai0HHqNSuaKOI8O59EhulZGvBR4jq8CLh4wfy6wY0XoGQX4tJCnPiKyWUIItJy2DmVGgV2vDUSgh5NF3ofz6BBR/Xk3B3pPlrbsJOD4ain0XNgFJB+Stk0zAL+oitAzEvBvI7USEZFDUygU8HbRIjWnGJn5JRYLOnXBoENE1+ceDHSfIG15acCJNVJrz7mtQFqctG19H/AKr2rpCerE0EPkwLwMOqTmFCPdRkZeMegQUe24+AGdH5e2wivAyd+l0HNmI5B5Ftj5ibS5BUsjt5oPBJr15rB1IgfjUzmXjo10SGbQIaK6c/IEOjwkbcV5wOk/pNBz+g8gJwnY9420KVRAcBcp9DQfAAR24nIURHbOq3LklY0MMedvHCK6NToXIPpuaSstkmZlPrMRiP9TaulJ3CNtW/4jzcwc1kcKPeEDpEtevMxFZFc8naWgk1XIoENE9kajly5btRouPb5yQQo+8X8CZ7cCRVnSDM0nVkvPe4RIrT3hA6R1uDhRIVGj5+6kAQBkFZTKXImEQYeILMczFIgZJ23GcmmSwvjN0pb4F5CVABxYIG2omNun8jJXcFcuQErUCHk4VwSdQgYdInIkShUQFCNtfV+W+vZc2FXR2rNZmrQw+aC0bf8voDEAzXpVtfj4tuJlLqJGoDLoZLNFh4gcms4FiBgqbQCQkyzNxhz/p3Sbf1nq3Hz6D+l518Cqvj3h/QEXX5kKJ6Ib8XBiHx3ZcQkIIhvkFgh0eFjajEYg7ZgUeuI3Awm7gdxk4PAiaQMA/7ZS8Gk+AAjpAWjkn5iMiAB3Z9vqo8MlILgEBJHtKy2Uwk78Zuky16Uj5s+r9dKaXJWXufyjpUVMicjqzl7Ow8D/2wpXnRpHZgyzyHtwCQgisi8ap4pOygOlx3mXKxYerejYnJtctRApADj7VF3maj5Aai0iIqvwqBhenltchtJyIzQqef/TwaBDRI2Piy/Q7j5pEwK4fLIq9JzfARSkA0d+ljYA8GkFhPYEmnYDQroBnmHs2ExkIZXDywEgp7AU3i46Gath0CGixk6hAPxaS1v3Z4GyEiBpb9VlrosHgfST0nZgvvQagx/QtGtF8OkOBLQH1PL+MiayFyqlAq56NXKLypDFoENE1MDUWmmNrWa9gUHTgYJM4MJOIGEPkLhXmssnP8184kKVTprDpzL8NO3GUV1Et8BNr0FuURlyi8rkLoVBh4jsnLMXEDlS2gBpmYqUw1XBJ/Ev6VJX5VIVlbzCgabdpfAT0l26/MUOzkS14qqX4kVukfwjrxh0iMixaCpGaIV0lx4LUbEm119V4efycWlf5lng78XScXp3abbmyn4+QTGA1iDfz0Fkw9z0Uj+dnEK26BARyUuhALybS1uHh6V9hVeApP0Vwecv4OIBoCgbOLNB2gBpZfYmbauCT9NugHuwfD8HkQ1hiw4RkS1z8gRaDpE2ACgvBVKPAgl/ScEn8S8g56J0CSzlMLD3a+k4t+CqS11Nu0qTGqr4a5Ycj1vFyKscBh0iokZApZE6Kwd2BLpPkPZlJ1Vc7qoIPpeOADlJwLEk4Ngv0jEaZ+kSV0h3qcUnuAvg5CHbj0FkLVUtOrx0RUTUOLkHS1v0PdLj4jxpQVJTq89eoDgbOL9d2ir5Rl7V6tNN6vTMOX3IzjDoyIxrXRFRg9O5AGF9pQ2Q1utKP3nV6K49Uufmy8el7eD30nHOPlX9fIK7SP1+dK7y/RxEDcDV1BlZ/ktXXOuKa10RkbXkXa7q45P4F5B8CCivYYVn7xZAk3bSRIaVm7OX9eslqqfFfyXgjRVHMDjSH9+O7dzg5+daV0REtsjFF4gcIW0AUFYMpPxdNbor+bDUzyfjjLRV9vUBAPcQIKAdENChIvy0A1ybyPFTEN2UQacCAOQVy9+iw6BDRCQXta5iNuauVfvy06Xwc/V25RyQnSBtlbM5A4CLv3mrT0B7wL0p+/yQ7Fx0UrzIL5a/iwiDDhGRLTH4AC0GSVulwixpVNelf6rCT/opIC8VOP2HtFVy8rzmslcHqcMzZ3UmK6oKOuyMTEREN+PkAYT1kbZKJflA6rGK4HMYSPkHSDsuTXZ4bqu0VdK6XBN+2klLWnCOH7IQQ0XQyWPQISKietEaql/2KiuWws7Vl71SjwIleUDCLmmrpNYD/m3ML3v5RXEVd2oQbNEhIqKGp9YBgR2krVJ5mXSZ6+rLXin/ACW50tIWFw9UHatUS/P8XB1+mkRzTS+qs8oWnfySchiNAkqlfP3GGHSIiOyZSg34R0lb+welfUaj1ME55bB560/hFSD1iLQd/rHiBArAJ6LqkldAe+kyGGd4phuobNEBgILScrPH1sagQ0TkaJTKqoVMK2d2FkJa1uLaEV95l6SJD9NPAkf+V3UOz2ZS4PFvA/hFSpe9PMPY74cAAHqNEkoFYBTS5SsGHSIikpdCAXg0lbbKeX4AIDe14rLX4arwk5UAXDkvbcdXVR2r0gG+EVLoqQw/vq2lIe8c9eVQFAoFDDo1covKkFdcBn8Za2HQISKi63P1B1yvWskdAAoyK8LPP8DlE0BaHHD5JFBaUDEM/oj5ObQuUuCpDD+Vty5+nPPHjhm0UtApLJF3Lh0GHSIiqhtnLyC8v7RVMhqBrAvSqK+0uIrb41JH6JI84OJ+abuak+dVweeqFiAud2EXnLXS7Mhyj7xi0CEiolunVAJeYdLW+vaq/eWlQEa8tJDp1SEo86zU+fnCTmm7mmuAefDxiwJ8W0kLp1Kj4VyxDERBKVt0iIjIXqk0gF9raWtzV9X+0kKptaey5adyy04AclOkLf5P83N5hJpf+vKLBHxacu4fG+WskSJGgczLQDDoEBGR9WmcqubquVpRjtTfx3T5K07qB5SXKl0ay7oAnFpXdbxCJY0e820ttfr4tJI6RHu34Pw/MjO16JTw0hUREZFE7wY07SJtV8vPqH75Ky0OKMqWWobSTwHHrzmXe4jU4uPbSrr1aSXNCWTwYSdoK6jso1PAzsjWFxsbi9jYWJSXy7+qKhER1YLBGzD0Bpr1rtonBJB7qarV5/JJIP20NOdPQUbViu/xm8zP5eQpBZ7KzbciAHmEAEqVdX8uO+ZUcemqkH10rG/ixImYOHEicnJy4O7uLnc5RERUHwoF4BYgbVev9g5ILUCVLT3ppypC0ClpDqDCK0DiX9J2NbVeuuTl0xLwbllx21y6r3ez3s9lJ5y00txJbNEhIiJqaAZvwNADCO1hvr+kAMiMN2/9uXwKyDgDlBVJi6CmHq1+Phd/KfB4NzcPQh6hnA36Opy1FS067KNDRERkJVpnoElbabuasVzq6Jx+uqIV6LQUfjLOSB2hK7cLO8xfp1RLS1/4tKxqDfJqLt138AkR9RrpMmBRqVHWOhh0iIiIlCrAK1zaIoaZP1eULQWe9Irgk3G66n5ZofQ443T1c2pdpPN5t5Bagrwq1hfzbuEQkyI6adgZmYiIyPbp3YGgGGm7mtEI5Cabt/6kn5YujWUlSDNCX/pH2qqd08M8/HiFS/e9wuwmBFWOuioqY9AhIiJqfJRKwD1Y2poPMH+urBi4ckEKPRlnpNmhM85IM0LnXASKsoCLB6TtWk6eVwWfcCkIeYZJ9529Gs3lML1G6ozMta6IiIjsjbpiJXffiOrPlRRIgSczXgpAmWerttwUaVTY9UKQzh3walZ1ma0yAHmFAS5NbGqV+Ko+Ogw6REREjkPrDDSJlrZrleQDmeeqglDmWSDjLHDlnNQSVJwNpPwtbddS6wHPZhXhJ+yq22bSHEFWXiqjMuhwHh0iIiKSaA3XD0GlhRWXwypaf66cq7qflSgNj798QtqqUUiX2DybXbNVBCELXBLjqCsiIiKqPY1T1QKp1yovBbITpdagK+cqbs9LW+Y5oDRfej47ETi/vfrrta4VwSe0ogUotCoMeYQAGn2dy9WrpctoxWzRISIiolui0lT127mWEEB+uhSAKoNP1oWqIJRzESjJBVKPSFtNXAMqwk9o9VvXwBonTWQfHSIiIrI8hQJw8ZW2pl2rP19aJA2Hvzr8XDkvXSa7cl4KQbkp0pa4p/rrlWrALagi/IQAHs0Az1B4lvvAH5koLvWz6I93Mww6REREjkyjv/4IMSGkUWBXB6CsC1IIyroAZCcB5SXS/awLZi8NAfCXHni77GkAw6qf20oYdIiIiKhmCoXUUdnZCwjqVP15o1Fq6TG1CF0w3S/LPA9kX8QlJVt0iIiIqDFSKgH3IGm7ZgFVNQCUl+FrWQq7pg4iIiKiBmcDK7vbzhSKRERERA2MQYeIiIjsFoMOERER2S0GHSIiIrJbDDpERERktxh0iIiIyG45ZNCJjY1FVFQUunTpIncpREREZEEKIYSQuwi55OTkwN3dHdnZ2XBzc5O7HCIiIqqFunx/O2SLDhERETkGBh0iIiKyWww6REREZLcYdIiIiMhuMegQERGR3ZJ/WVEZVQ44y8nJkbkSIiIiqq3K7+3aDBx36KCTm5sLAGjatKnMlRAREVFd5ebmwt3d/YbHOPQ8OkajEcnJyXB1dYVCoWjQc+fk5KBp06ZITEzkHD0WxM/ZOvg5Ww8/a+vg52wdlvqchRDIzc1FYGAglMob98Jx6BYdpVKJ4OBgi76Hm5sb/xFZAT9n6+DnbD38rK2Dn7N1WOJzvllLTiV2RiYiIiK7xaBDREREdotBx0J0Oh3efvtt6HQ6uUuxa/ycrYOfs/Xws7YOfs7WYQufs0N3RiYiIiL7xhYdIiIislsMOkRERGS3GHSIiIjIbjHoEBERkd1i0LkFX375JcLCwqDX6xETE4Pt27ff8PitW7ciJiYGer0e4eHh+Oqrr6xUaeNWl8/5l19+wZAhQ+Dr6ws3Nzf06NED69evt2K1jVdd/z5X2rlzJ9RqNTp06GDZAu1EXT/n4uJiTJs2DaGhodDpdGjevDm+++47K1XbuNX1s160aBHat28PZ2dnBAQE4PHHH0dGRoaVqm18tm3bhpEjRyIwMBAKhQIrV6686Wtk+R4UVC8//fST0Gg04ptvvhFxcXHixRdfFAaDQVy4cKHG48+ePSucnZ3Fiy++KOLi4sQ333wjNBqNWLZsmZUrb1zq+jm/+OKL4oMPPhB79+4Vp06dEq+//rrQaDTi4MGDVq68canr51wpKytLhIeHi6FDh4r27dtbp9hGrD6f86hRo0S3bt3Ehg0bxLlz58Rff/0ldu7cacWqG6e6ftbbt28XSqVSfPrpp+Ls2bNi+/btok2bNmL06NFWrrzxWLt2rZg2bZpYvny5ACBWrFhxw+Pl+h5k0Kmnrl27igkTJpjta926tXjttddqPP7VV18VrVu3Ntv3zDPPiO7du1usRntQ18+5JlFRUWLGjBkNXZpdqe/n/MADD4g333xTvP322ww6tVDXz3ndunXC3d1dZGRkWKM8u1LXz/qjjz4S4eHhZvs+++wzERwcbLEa7Ultgo5c34O8dFUPJSUlOHDgAIYOHWq2f+jQodi1a1eNr9m9e3e144cNG4b9+/ejtLTUYrU2ZvX5nK9lNBqRm5sLLy8vS5RoF+r7Oc+fPx/x8fF4++23LV2iXajP57xq1Sp07twZH374IYKCghAREYGXX34ZhYWF1ii50arPZ92zZ08kJSVh7dq1EEIgNTUVy5Ytwx133GGNkh2CXN+DDr2oZ32lp6ejvLwc/v7+Zvv9/f1x6dKlGl9z6dKlGo8vKytDeno6AgICLFZvY1Wfz/la//d//4f8/Hzcf//9lijRLtTncz59+jRee+01bN++HWo1f43URn0+57Nnz2LHjh3Q6/VYsWIF0tPT8dxzzyEzM5P9dG6gPp91z549sWjRIjzwwAMoKipCWVkZRo0ahc8//9waJTsEub4H2aJzCxQKhdljIUS1fTc7vqb9ZK6un3OlJUuW4J133sHSpUvh5+dnqfLsRm0/5/Lycjz88MOYMWMGIiIirFWe3ajL32ej0QiFQoFFixaha9euuP322zF79mwsWLCArTq1UJfPOi4uDi+88ALeeustHDhwAL///jvOnTuHCRMmWKNUhyHH9yD/K1YPPj4+UKlU1f5nkJaWVi2tVmrSpEmNx6vVanh7e1us1sasPp9zpaVLl+LJJ5/Ezz//jMGDB1uyzEavrp9zbm4u9u/fj0OHDuH5558HIH0hCyGgVqvxxx9/YODAgVapvTGpz9/ngIAABAUFwd3d3bQvMjISQggkJSWhZcuWFq25sarPZz1r1iz06tULr7zyCgCgXbt2MBgM6NOnD9577z22ujcAub4H2aJTD1qtFjExMdiwYYPZ/g0bNqBnz541vqZHjx7Vjv/jjz/QuXNnaDQai9XamNXncwaklpxx48Zh8eLFvL5eC3X9nN3c3HDkyBEcPnzYtE2YMAGtWrXC4cOH0a1bN2uV3qjU5+9zr169kJycjLy8PNO+U6dOQalUIjg42KL1Nmb1+awLCgqgVJp/JapUKgBVrQ50a2T7HrRoV2c7Vjl0cd68eSIuLk5MnjxZGAwGcf78eSGEEK+99poYM2aM6fjKYXVTpkwRcXFxYt68eRxeXgt1/ZwXL14s1Gq1iI2NFSkpKaYtKytLrh+hUajr53wtjrqqnbp+zrm5uSI4OFjce++94tixY2Lr1q2iZcuWYvz48XL9CI1GXT/r+fPnC7VaLb788ksRHx8vduzYITp37iy6du0q149g83Jzc8WhQ4fEoUOHBAAxe/ZscejQIdMQflv5HmTQuQWxsbEiNDRUaLVa0alTJ7F161bTc2PHjhX9+vUzO37Lli2iY8eOQqvVimbNmok5c+ZYueLGqS6fc79+/QSAatvYsWOtX3gjU9e/z1dj0Km9un7Ox48fF4MHDxZOTk4iODhYTJ06VRQUFFi56saprp/1Z599JqKiooSTk5MICAgQjzzyiEhKSrJy1Y3H5s2bb/j71la+BxVCsE2OiIiI7BP76BAREZHdYtAhIiIiu8WgQ0RERHaLQYeIiIjsFoMOERER2S0GHSIiIrJbDDpERERktxh0iIiIyG4x6BAREZHdYtAhIrpG//79MXny5Fs6hxACTz/9NLy8vKBQKHD48OEGqY2I6kYtdwFERDcybtw4NGnSBO+//77cpdTJ77//jgULFmDLli0IDw+Hj4+P3CUROSQGHSKyWUajEWvWrMGqVavkLqXO4uPjERAQgJ49e173mJKSEmi1WitWReR4eOmKyAEtWbIEer0eFy9eNO0bP3482rVrh+zs7HqfNzg4GF9++aXZvl27dsHZ2RkXLlyo8/l27twJpVKJbt261fh8//79MWnSJEyePBmenp7w9/fH3LlzkZ+fj8cffxyurq5o3rw51q1bZ/a64uJivPDCC/Dz84Ner0fv3r2xb9++69YhhMCHH36I8PBwODk5oX379li2bNl1jx83bhwmTZqEhIQEKBQKNGvWzFTv888/j6lTp8LHxwdDhgwBILX+9O7dGx4eHvD29saIESMQHx9vOp/RaMQHH3yAFi1aQKfTISQkBDNnzqztx0jk0Bh0iBzQgw8+iFatWmHWrFkAgBkzZmD9+vVYt24d3N3d633e7t27mwUGIQQmT56MyZMnIzQ0tM7nW7VqFUaOHAml8vq/qr7//nv4+Phg7969mDRpEp599lncd9996NmzJw4ePIhhw4ZhzJgxKCgoML3m1VdfxfLly/H999/j4MGDaNGiBYYNG4bMzMwa3+PNN9/E/PnzMWfOHBw7dgxTpkzBo48+iq1bt9Z4/Keffop3330XwcHBSElJMftMvv/+e6jVauzcuRNff/01ACA/Px9Tp07Fvn37sGnTJiiVStx1110wGo0AgNdffx0ffPABpk+fjri4OCxevBj+/v51/jyJHJIgIof022+/CZ1OJ2bOnCk8PT3F0aNHTc+NHj1aeHh4iHvuuadO5/zoo49EmzZtTI+///574e/vL3Jycup13oiICLFq1arrPt+vXz/Ru3dv0+OysjJhMBjEmDFjTPtSUlIEALF7924hhBB5eXlCo9GIRYsWmY4pKSkRgYGB4sMPPzSd98UXXzQdr9frxa5du8ze+8knnxQPPfTQdWv7+OOPRWhoaLV6O3TocOMfWgiRlpYmAIgjR46InJwcodPpxDfffHPT1xFRdWzRIXJQI0aMQFRUFGbMmIEVK1agTZs2pudeeOEFLFy4sM7n7N69O44fP468vDwUFBTgjTfewHvvvQdXV9c6n/f48eNISkrC4MGDb3hcu3btTPdVKhW8vb3Rtm1b077Klo+0tDQAUt+Z0tJS9OrVy3SMRqNB165dcfz48Wrnj4uLQ1FREYYMGQIXFxfTtnDhQrPLS7XVuXPnavvi4+Px8MMPIzw8HG5ubggLCwMAJCQk4Pjx4yguLsagQYPq/F5ExM7IRA5r/fr1OHHiBMrLy6tdBhkwYAC2bNlS53N27twZKpUKBw8exMaNG+Ht7Y0nnniiXuddtWoVhgwZAicnpxsep9FozB4rFAqzfQqFAgBMl4GEEGb7Kwkhqu27+nVr1qxBUFCQ2XM6na42P4oZg8FQbd/IkSPRtGlTfPPNNwgMDITRaER0dDRKSkpu+vMT0Y2xRYfIAR08eBD33Xcfvv76awwbNgzTp09vkPPq9Xq0b98ev/zyC/773/9i9uzZN+xfcyO//vorRo0a1SB1Xa1FixbQarXYsWOHaV9paSn279+PyMjIasdHRUVBp9MhISEBLVq0MNuaNm16y/VkZGTg+PHjePPNNzFo0CBERkbiypUrpudbtmwJJycnbNq06Zbfi8gRsUWHyMGcP38ed9xxB1577TWMGTMGUVFR6NKlCw4cOICYmJhbPn/37t3x2WefYcSIEfW+3JKWloZ9+/Zh5cqVt1zPtQwGA5599lm88sor8PLyQkhICD788EMUFBTgySefrHa8q6srXn75ZUyZMgVGoxG9e/dGTk4Odu3aBRcXF4wdO/aW6vH09IS3tzfmzp2LgIAAJCQk4LXXXjM9r9fr8a9//QuvvvoqtFotevXqhcuXL+PYsWOmer/44gusWLGCYYioBgw6RA4kMzMTw4cPx6hRo/DGG28AAGJiYjBy5EhMmzYNv//++03PsWDBAjz++OOmS0DX6tChA9RqNT766KN61/nbb7+hW7du8PPzq/c5buT999+H0WjEmDFjkJubi86dO2P9+vXw9PSs8fh///vf8PPzw6xZs3D27Fl4eHigU6dOps/wViiVSvz000944YUXEB0djVatWuGzzz5D//79TcdMnz4darUab731FpKTkxEQEIAJEyaYnk9PT69XfyEiR6AQ1/ttRUQObcuWLfjiiy+qzRfzzjvvYMuWLdftazNw4EC0a9cOn3zySZ3Oe7VRo0ahd+/eePXVV+tbPhERALboEFENhg0bhoMHDyI/Px/BwcFYsWIFunTpAkDqxPzpp5+aHW80GnH58mXMmzcPJ0+exIoVK+p83qv17t0bDz30UMP/YETkcNiiQ0S3bMuWLRg4cCBat26N+fPnX3cmYyIia2PQISIiIrvF4eVERERktxh0iIiIyG4x6BAREZHdYtAhIiIiu8WgQ0RERHaLQYeIiIjsFoMOERER2S0GHSIiIrJbDDpERERktxh0iIiIyG4x6BAREZHd+n8ERPkVdmDoSAAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "plt.plot(df['xL_0 / mole frac.'], df['pL / Pa']/1e3)\n", + "plt.plot(df['xV_0 / mole frac.'], df['pL / Pa']/1e3)\n", + "plt.yscale('log')\n", + "plt.gca().set(xlabel='$x_1,y_1$ / mole frac.', ylabel='p / kPa')\n", + "plt.title(f'Water(1) + Ammonia(2) @ {T} K')\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.10.6" + } + }, + "nbformat": 4, + "nbformat_minor": 5 +} diff --git a/doc/source/models/index.rst b/doc/source/models/index.rst index 39f5a3a..e723a4e 100644 --- a/doc/source/models/index.rst +++ b/doc/source/models/index.rst @@ -27,5 +27,6 @@ Models REFPROP_conversion GERG ECS + MultifluidPlusAssociation IdealGas \ No newline at end of file