{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "e38e6eeb-ad51-4eb1-8964-cddbb934d7b2", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "import sys\n", "sys.path.append(\"..\") # This allows you to put the python file in the folder above. So you can have lots of sub folders\n", "import matplotlib.pyplot as plt\n", "import Thermobar as pt" ] }, { "cell_type": "code", "execution_count": 29, "id": "5d4729cf-0526-4e08-b845-80083967d1fc", "metadata": {}, "outputs": [], "source": [ "out=pt.import_excel('Waters2021.xlsx', sheet_name=\"Sheet1\")\n", "my_input=out['my_input']\n", "myLiquids1=out['Liqs']\n", "myCpxs1=out['Cpxs']" ] }, { "cell_type": "code", "execution_count": 36, "id": "2cf223e9-69e8-440f-ba78-74f544e22667", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'C:\\\\Users\\\\penny\\\\anaconda3'" ] }, "execution_count": 36, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import os\n", "import sys\n", "os.path.dirname(sys.executable)" ] }, { "cell_type": "code", "execution_count": 34, "id": "3717d5cf-20f6-419e-9bb7-d6a10d6571da", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 34, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAAAX3UlEQVR4nO3df2zc933f8eeL5GLgiqZIKWqoLZHnFnIg2a1dmRMaYFa6Fqs9IY0qFy1U3NY0HnytqBZSinWIwC7YVnCokgZFh0IKrlX6k7Hi1mnrLLMQ2UCU/hHVozorpaxoYmaepNhzmNjoCrCTRfK9P+5L6agcyeOPu+/3vvd6AMQdP/e505sfn1/3vc/n+0MRgZmZ5UtP2gWYmdnmc7ibmeWQw93MLIcc7mZmOeRwNzPLob60CwDYsmVLFIvFtMswM+soFy5c+FZEDDR6LBPhXiwWmZiYSLsMM7OOIqm63GOeljEzyyGHu5lZDjnczcxyyOFuZpZDDnczsxxyuJulaHx8nGKxSE9PD8VikfHx8bRLspzIxK6QZt1ofHyccrnM7OwsANVqlXK5DECpVEqzNMsBb7mbpWR0dPR2sC+anZ1ldHQ0pYosTxzuZim5du3amtrN1sLhbpaSwcHBNbWbrYXD3SwlY2NjFAqFJW2FQoGxsbGUKrI8cbibpaRUKlGpVBgaGkISQ0NDVCoVL6baplAWrqE6PDwcPnGYmdnaSLoQEcONHvOWu5lZDjnczcxyyOFuZpZDDnczsxxyuJuZ5ZDD3cwshxzuZmY55HA3M8shh7uZWQ41Fe6SjkialHRJ0tGk7WFJX5H0d5I+L+nddf2PSZqSdEXS4y2q3czMlrFquEt6CHga2AM8DHxA0g7g94GPRsQPAn8B/FrSfxdwEHgQeAI4Iam3NeWbmVkjzWy57wTOR8RsRMwB54ADwHuBLyd9zgI/ndzfD5yOiJsR8RowRe2DwczM2qSZcJ8E9krql1QA9gHbk/YPJn1+JmkDuA+4Xvf8G0nbEpLKkiYkTczMzKy3fjMza2DVcI+Iy8BxalvnZ4CLwBzwFHBY0gXgu4F3kqeo0cs0eN1KRAxHxPDAwMA6yzczs0aaWlCNiFMRsTsi9gJvAVcj4msR8RMR8SjwDPD1pPsN7mzFA2wDXt/Mos3MbGXN7i2zNbkdBJ4Enqlr6wF+HfhU0v154KCkeyTdD+wAXt7sws3MbHl9TfZ7TlI/cAs4HBFvJ7tHHk4e/xzwBwARcUnSs8Cr1KZvDkfE/GYXbmZmy/OVmMzMOpSvxGRm1mUc7mYpGRkZoa+vD0n09fUxMjKSdkmWI83OuZvZJhoZGeHkyZO3f5+fn7/9+4kTJ9Iqy3LEc+5mKejr62N+/jv3M+jt7WVubi6FiqwTec7dLGMaBftK7WZr5XA3S0Fvb+Nz6S3XbrZWDnezFJTL5TW1m62VF1TNUrC4aFqpVJifn6e3t5dyuezFVNs0XlA1M+tQXlA1M+syDnczsxxyuJuZ5ZDD3cwshxzuZmY55HA3M8shh7uZWQ453M3McsjhbmaWQw53M7MccribmeWQw93MLIeaCndJRyRNSrok6WjS9oik85JekTQhaU9d/2OSpiRdkfR4i2o3M7NlrHrKX0kPAU8De4B3gDOSvgB8HPhPEfGCpH3J7z8qaRdwEHgQuBd4UdIDEeFLzJiZtUkzW+47gfMRMRsRc8A54AAQwLuTPt8DvJ7c3w+cjoibEfEaMEXtg8HMzNqkmXCfBPZK6pdUAPYB24GjwCckXQd+CziW9L8PuF73/BtJ2xKSysl0zsTMzMwG/gTrFOPj4xSLRXp6eigWi4yPj6ddkllurRruEXEZOA6cBc4AF4E54BDwkYjYDnwEOJU8RY1epsHrViJiOCKGBwYG1lm+dYrx8XHK5TLVapWIoFqtUi6XHfBmLdLUgmpEnIqI3RGxF3gLuAp8CPhc0uXPuDP1coPalv2ibdyZsrEuNTo6yuzs7JK22dlZRkdHU6rILN+a3Vtma3I7CDwJPEMtsN+fdPkxaoEP8DxwUNI9ku4HdgAvb2bR1nmuXbu2pva73T2lMzIy4ikesxU0e4Hs5yT1A7eAwxHxtqSngd+R1Af8P6AMEBGXJD0LvEpt+uaw95SxwcFBqtVqw/bVLE7pLG75V6tVTp48efvxxSkegFKptEkVm3U2XyDb2uLugAYoFApUKpVVA7lYLDb8YLjb0NAQ09PTGy3VrGP4AtmWulKpRKVSYWhoCEkMDQ01FezQ/NRNs/3MuoHDPSO6YTfBUqnE9PQ0CwsLTE9PNz2F0szUzVr6mXUDh3sGeDfBlY2NjVEoFFbsUygUGBsba1NFZtnncM8A7ya4skZTOocOHVrXFI9Zt/CCagb09PTQ6L+DJBYWFlKoyMw6gRdUM265uWLPIZvZejncM6DRnLLnkM1sIxzuGbCR3QTNzBrxnLuZWYfynLuZWZdxuFvquuEALrN2a/bEYWYt0eikYD4JmNnGecvdUuUDuMxaw+Fuqdroed7NrDGHu6XKB3CZtYbD3VLlA7jMWsPhbqnyAVxmreGDmMzMOpQPYjIz6zIOdzOzHHK4m5mloNVHZjcV7pKOSJqUdEnS0aTts5JeSX6mJb1S1/+YpClJVyQ9vqkVm5l1uHZcWnPVBVVJDwGngT3AO8AZ4FBEXK3r80ng7yPiP0vaBTyT9L8XeBF4ICLml/s3vKBqZt2kWCxSrVa/o31oaIjp6emmX2ejC6o7gfMRMRsRc8A54EDdiwv4WWqBDrAfOB0RNyPiNWCKWtCbmRntOTK7mXCfBPZK6pdUAPYB2+sefwx4s25L/j7get3jN5K2JSSVJU1ImpiZmVlf9WZmHagdR2avGu4RcRk4DpylNiVzEZir6/Jz3NlqB1Cjl2nwupWIGI6I4YGBgTUVbWbWydpxZHZTC6oRcSoidkfEXuAt4CqApD7gSeCzdd1vsHTLfhvw+uaUa2bW+dpxZHZTR6hK2hoR35Q0CHwReF9EvC3pCeBYRLy/ru+DwGe4s6D6ErDDC6pmZptrpQXVZi/W8ZykfuAWcDgi3k7aD7J0SoaIuCTpWeBVatM3h1cKdjMz23zNTss8FhG7IuLhiHiprv0XIuJTDfqPRcQPRMR7I+KFzSzY1seXsjPrLr7MXhfwpezMuk9Hn35gZGSEvr4+JNHX18fIyEjaJWWSL2Vnlj2t/jbdsVvuIyMjnDx58vbv8/Pzt38/ceJEWmVlki9lZ5Yt7fg23bHnc+/r62N+/jvXaXt7e5mbm2vwjO61WYc6m9nmyMrpBzKpUbCv1N7NfCk7s2zJyukHMqm3t3dN7d3Ml7Izy5ZMnH4gqxbnp5pt73alUonp6WkWFhaYnp52sJulKDOnH8iiEydOcOjQodtb6r29vRw6dMiLqWaWeZk5/UCr+fQDZmZrl8sFVTMzW57D3cwsh7oq3H1+FTPrFh17hOpa+fwqZtZNumbL3edXMbNu0jXh7vOrmFk36Zpwb8cRYWZmWdE14e7zq5hZN+macPf5Vcysm/gIVTOzDuUjVM3MuozD3cwsh5oKd0lHJE1KuiTpaF37r0i6krR/vK79mKSp5LHHW1C3mZmtYNUjVCU9BDwN7AHeAc5I+gKwDdgP/FBE3JS0Nem/CzgIPAjcC7wo6YGI8CWSzMzapJkt953A+YiYjYg54BxwADgE/GZE3ASIiG8m/fcDpyPiZkS8BkxR+2AwM7M2aSbcJ4G9kvolFYB9wHbgAeAxSX8j6Zykf5b0vw+4Xvf8G0nbEpLKkiYkTczMzGzsrzAzsyVWDfeIuAwcB84CZ4CLwBy1KZ33AD8C/BrwrCQBavQyDV63EhHDETE8MDCw/r/AzDqaz9baGk0tqEbEqYjYHRF7gbeAq9S2yD8XNS8DC8CWpH173dO3Aa9vbtlmlgeLZ2utVqtExO2ztTrgN67ZvWUWF0sHgSeBZ4C/BH4saX8AeBfwLeB54KCkeyTdD+wAXt70ys2s4/lsra3T7H7uz0l6Ffg8cDgi3gY+DXy/pEngNPChZCv+EvAs8Cq1aZzD3lPGVuKv5d3LZ2ttnaYu1hERjzVoewf418v0HwN8Ri5blS+i0t0GBwepVqsN221jfISqpcpfy7ubz9baOg53S5W/lnc3n621dXxWSEtVsVhs+LV8aGiI6enp9hdk1kF8VkjLLH8tN2sNh7ulyl/LzVrD4W5ts9wuj6VSienpaRYWFpiennawdxnvCtsaTe0KabZR3uXRGvH7onW8oGpt4YVTa8Tvi43xgqqlzrs8WiN+X7SOw93aYrkjDn0kYnfz+6J1HO7WFt7l0Rrx+6J1HO7WFt7l0Rrx+6J1vKBqZtahvKBqZtZlHO5mZjnkcDczyyGHe5fwId5m3cWnH+gCPsTbrPt4y70L+GpHZt3H4d4FfIi3WfdxuHcBH+Jt1n2aCndJRyRNSrok6WjS9h8lfUPSK8nPvrr+xyRNSboi6fEW1W5N8iHeZt1n1XCX9BDwNLAHeBj4gKQdycO/HRGPJD//Pem/CzgIPAg8AZyQ1NuS6q0pPsTbrPs0s7fMTuB8RMwCSDoHHFih/37gdETcBF6TNEXtg+ErGy3W1q9UKjnMzbpIM9Myk8BeSf2SCsA+YHvy2C9L+qqkT0t6T9J2H3C97vk3krYlJJUlTUiamJmZ2cCfYGZmd1s13CPiMnAcOAucAS4Cc8BJ4AeAR4A3gE8mT1Gjl2nwupWIGI6I4YGBgXUVb2ZmjTW1oBoRpyJid0TsBd4CrkbEmxExHxELwO9Rm3qB2pb69rqnbwNe38yizfLCRw5bqzS7t8zW5HYQeBJ4RtL31XU5QG36BuB54KCkeyTdD+wAXt68ks3yYfHI4Wq1SkTcPnLYAW+boanzuUv6a6AfuAX8akS8JOlPqE3JBDAN/GJEvJH0HwWeojZ9czQiXljp9X0+d+tGvji0bdRK53P3xTrMUtLT00Oj//8ksbCwkEJF1ml8sQ6zDPKRw9ZKDnezlPjIYWslh7tZSnzksLWS59zNzDqU59zNzLqMw93MLIcc7mZmOeRwNzPLIYe7mVkOOdzNzHLI4W5mlkMOdzOzHHK4m5nlkMPdzCyHHO5mZjnkcDczyyGHu5lZDjnczcxyyOFuZpZDXRXu4+PjFItFenp6KBaLvsq8meVWX9oFtMv4+DjlcpnZ2VkAqtUq5XIZwFe+MbPcaWrLXdIRSZOSLkk6etdj/05SSNpS13ZM0pSkK5Ie3+Sa12V0dPR2sC+anZ1ldHQ0pYrMzFpn1S13SQ8BTwN7gHeAM5K+EBFXJW0H/iVwra7/LuAg8CBwL/CipAciYr4Vf0Czrl27tqZ2M7NO1syW+07gfETMRsQccA44kDz228C/B+ovxLofOB0RNyPiNWCK2gdDqgYHB9fUbmbWyZoJ90lgr6R+SQVgH7Bd0geBb0TExbv63wdcr/v9RtK2hKSypAlJEzMzM+ssv3ljY2MUCoUlbYVCgbGxsZb/22Zm7bZquEfEZeA4cBY4A1wE5oBR4GMNnqJGL9PgdSsRMRwRwwMDA2sqej1KpRKVSoWhoSEkMTQ0RKVS8WKqmeVSU3vLRMQp4BSApP8CvAmUgIuSALYBfytpD7Ut9e11T98GvL6JNa9bqVRymJtZV2h2b5mtye0g8CTwxxGxNSKKEVGkFui7I+L/AM8DByXdI+l+YAfwckuqNzOzhprdz/05Sf3ALeBwRLy9XMeIuCTpWeBVatM3h9PeU8bMrNs0Oy3z2CqPF+/6fQzwSqWZWUq66vQDZmbdwuFuZpZDDnczsxxyuJuZ5ZDD3cwshxzuZmY55HA3M8shh7uZWQ453M3McsjhbmaWQw53M7MccribmeWQw93MLIcc7mZmOeRwNzPLIYe7mVkOOdzNzHLI4Z4R4+PjFItFenp6KBaLjI+Pp12SmXWwZq+hai00Pj5OuVxmdnYWgGq1SrlcBqBUKqVZmpl1KG+5Z8Do6OjtYF80OzvL6OhoShWZWadzuGfAtWvX1tRuZraapsJd0hFJk5IuSTqatP2GpK9KekXSFyXdW9f/mKQpSVckPd6i2nNjcHBwTe3mNQqz1awa7pIeAp4G9gAPAx+QtAP4RET8UEQ8Avw34GNJ/13AQeBB4AnghKTe1pSfD2NjYxQKhSVthUKBsbGxlCrKnvow37JlCx/+8IepVqtExO01Cge82R3NbLnvBM5HxGxEzAHngAMR8X/r+nwXEMn9/cDpiLgZEa8BU9Q+GGwZpVKJSqXC0NAQkhgaGqJSqXgxNbG44LwY5t/+9re5devWkj5eozBbShGxcgdpJ/BXwPuAfwReAiYi4lckjQE/D/w98C8iYkbS71L7MPjT5PmngBci4s/vet0yUAYYHBx8tFqtbu5fZrlRLBZp5v0hiYWFhTZUZJYNki5ExHCjx1bdco+Iy8Bx4CxwBrgIzCWPjUbEdmAc+OXFf6/RyzR43UpEDEfE8MDAQFN/iHWnZheWvUZhdkdTC6oRcSoidkfEXuAt4OpdXT4D/HRy/wawve6xbcDrGy3Uulczoe01CrOlmt1bZmtyOwg8CTyTLKou+iDwteT+88BBSfdIuh/YAby8eSVbt2m04Pyud72L/v5+r1GYLaPZI1Sfk9QP3AIOR8Tbkn5f0nuBBaAK/BJARFyS9CzwKrXpm8MRMd+C2q1LLIb26Ogo165dY3BwkLGxMYe52QpWXVBth+Hh4ZiYmEi7DDOzjrKhBVUzM+s8DnczsxxyuJuZ5ZDD3cwshxzuZmY5lIm9ZSTNUNudcjlbgG+1qZyNcq2t4Vpbo5Nqhc6qtx21DkVEw0P8MxHuq5E0sdzuPlnjWlvDtbZGJ9UKnVVv2rV6WsbMLIcc7mZmOdQp4V5Ju4A1cK2t4Vpbo5Nqhc6qN9VaO2LO3czM1qZTttzNzGwNHO5mZjmUiXCXdETSpKRLko4mbZ+V9EryMy3plbr+xyRNSboi6fGs1iqpKOkf6x77VAZqfUTS+aSeCUl76vpnbVwb1prRcX1Y0lck/Z2kz0t6d13/1MZ1rfW2e2wlfVrSNyVN1rV9r6Szkq4mt++pe6zhWEp6NPlbpiT9V0mNrgiXlVq/lLQtjvHWza4VgIhI9Qd4CJgECtTOL/8isOOuPp8EPpbc30XtUn/3APcDXwd6M1prEZjM0rgCXwT+VdJnH/ClrI7rCrVmcVz/B/D+pM9TwG+kPa7rrLetYwvsBXbX/5vAx4GPJvc/ChxfbSypXRDofdQu8/nC4vsmo7V+CRhu9dhmYct9J7ULas9GxBxwDjiw+GDyCfyzwDNJ037gdETcjIjXgClgD+2x1lrTtFytASxuVX4Pdy6BmMVxXa7WNC1X63uBLyd9znLnspNpjut66m2riPgytUt31tsP/FFy/4+An6pr/46xlPR9wLsj4itRS88/rntOpmrd7JpWkoVwnwT2SuqXVKC2hVZ/DdbHgDcjYvG6rfcB1+sev5G0tcNaawW4X9L/lHRO0mNtqnOlWo8Cn5B0Hfgt4FjSP4vjulytkL1xnaR2uUmAn+HO+yLNcYW11wvpje2ifxoRbwAkt4vTFsuN5X3J/bvb22GttS76g2RK5j+0YgoJMhDuEXEZOE5t6+EMta8yc3Vdfo6lW8KNBqIt+3Ouo9Y3gMGI+GHgV4HP1M/FplTrIeAjEbEd+AhwKnlKFsd1uVqzOK5PAYclXQC+G3gneUpq4wrrqje1sW3CcmOZ6hgvY6WaShHxg9Q2Bh8D/k0rCkg93AEi4lRE7I6IvdS+9lwFkNRH7YLcn63rfoOlWxnbaOPX9bXUmnwl+3Zy/wK1ebcHUq71Q8Dnki5/xp2vilkc14a1ZnFcI+JrEfETEfEotQ/4ryfdUx3Xtdab9tgm3kymWkhuv5m0LzeWN5L7d7e3w1prJSK+kdz+A/AZWjVd0+pJ/WZ+gK3J7SDwNeA9ye9PAOfu6vsgSxcq/jftXaBaS60D3FlE+X7gG8D3plkrcBn40aT9x4ELWR3XFWrN4rgutvVQm/N9Kgvjuo562z623LWIC3yCpYuUH19tLKktEP8IdxZU92WxVmqL2luSPv8E+HPgl1pSazvfZCsM2F8DryaD8eN17X/Y6A8HRqltUVyhBavim1UrtUWqS0nfvwV+Mu1agX8OXEja/gZ4NKvjulytGR3XI8D/Sn5+k+To77THda31tntsqX1reAO4RW1r998C/cBL1L69vUTdh8tyYwkMU1tH+Drwu/Xjn6Vage9K3tNfTcb5d2jRh71PP2BmlkOZmHM3M7PN5XA3M8shh7uZWQ453M3McsjhbmaWQw53M7MccribmeXQ/wf4zyMo2+cT5wAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "T34=pt.calculate_liq_only_temp(liq_comps=myLiquids1, equationT=\"T_Put2008_eq34_cpx_sat\", P=10*my_input['P'])-273.15\n", "plt.plot(T34, my_input['temp_out'], 'ok')" ] }, { "cell_type": "code", "execution_count": 32, "id": "61c50f6b-3bd4-4dfb-a3b7-cd0fec7c40b8", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "C:\\Users\\penny\\anaconda3\\lib\\site-packages\\pandas\\core\\arraylike.py:358: RuntimeWarning: divide by zero encountered in log\n", " result = getattr(ufunc, method)(*inputs, **kwargs)\n" ] }, { "data": { "text/plain": [ "[]" ] }, "execution_count": 32, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAD4CAYAAAAEhuazAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAAAXEElEQVR4nO3db4xcV33G8e+zsWp1K4Kos65I7J1JkR3ZTknkLBa8iPlTFacWwnUQyGiipkqVgWhBdiqosLaiFLQVTkCIqgK04FAggx3TpG0ojYUTqaYvYtJ1FcOuHctLvbM2SZNFieiLLU7s/fXF3LVn7VnP7OyuZ+7c5yONdnzm3Nlz7PEz955zz72KCMzMLFu6Wt0AMzO79hz+ZmYZ5PA3M8sgh7+ZWQY5/M3MMmhZqxsAcMMNN0Q+n291M8zMUuXo0aO/ioieZrZti/DP5/MMDw+3uhlmZqkiqdzsth72MTPLIIe/mVkGOfzNzDLI4W9mlkEOfzOzDHL4m5k1oFQqkc/n6erqIp/PUyqVWt2kBWmLUz3NzNpZqVSiWCwyNTUFQLlcplgsAlAoFFrZtKZ5z9/MrI6BgYGLwT9jamqKgYGBFrVo4Rz+ZmZ1TExMzKs8DRz+ZmZ19Pb2zqs8DRz+ZmZ1DA4O0t3dPausu7ubwcHBFrVo4Rz+ZmZ1FAoFhoaGyOVySCKXyzE0NJTayV4AtcM9fPv6+sIXdjMzmx9JRyOir5ltvedvZpZBDn8zswxy+JuZZZDD38wsgxz+ZmYZ5PA3M8sgh7+ZWQY5/M3MMsjhb2aWQQ2Fv6SdkkYkjUralZTdJulZST+X9ENJ11fV3y1pTNJJSVuWqO1mZtakuuEv6VbgfmATcBvwAUlrgG8Bn4mIPwD+Cfh0Un89sAPYANwFfE3SdUvTfDMza0Yje/7rgCMRMRUR54HDwHbgFuAnSZ1DwIeS59uA/RFxLiJOA2NUvjjMzKxNNBL+I8BmSSskdQNbgdVJ+QeTOh9OygBuAs5UbX82KZtFUlHSsKThycnJZttvZmZNqBv+EXEC2ENl7/4gcAw4D9wH9Es6CrwJeD3ZRLXepsb7DkVEX0T09fT0NNl8MzNrRkMTvhGxNyI2RsRm4FXgVES8EBHvj4g7gH3AL5LqZ7l0FACwCnhxMRttZmYL0+jZPiuTn73A3cC+qrIu4K+AbyTVnwR2SFou6WZgDfDcYjfczMyat6zBeo9LWgG8AfRHxGvJ6Z/9yetPAN8GiIhRSQeA41SGh/oj4sJiN9zMzJrnO3mZmaWU7+RlZmbz4vA3M8sgh7+ZWQY5/M3MMsjhb2aWQQ5/M7MMcvibmWWQw9/MLIMc/mZmGeTwNzPLIIe/mVkGOfzNzDLI4W9mlkEOfzOzDHL4m5llkMPfzCyDHP5mZhnk8DczyyCHv5lZBjn8zcwyyOFvZpZBDYW/pJ2SRiSNStqVlN0u6Yik5yUNS9pUVX+3pDFJJyVtWaK2m5lZk5bVqyDpVuB+YBPwOnBQ0o+Ah4C/iYinJG1N/vweSeuBHcAG4EbgaUlrI+LCUnXCzMzmp5E9/3XAkYiYiojzwGFgOxDA9UmdNwMvJs+3Afsj4lxEnAbGqHxxmJlZm2gk/EeAzZJWSOoGtgKrgV3Aw5LOAF8Cdif1bwLOVG1/NimbRVIxGS4anpycXEAXLMtKpRL5fJ6uri7y+TylUqnVTTJLhbrhHxEngD3AIeAgcAw4DzwAPBgRq4EHgb3JJqr1NjXedygi+iKir6enp8nmW5aVSiWKxSLlcpmIoFwuUywW/QVg1oCGJnwjYm9EbIyIzcCrwCngXuCJpMoPuDS0c5bKkcGMVVwaEjJbNAMDA0xNTc0qm5qaYmBgoEUtMkuPRs/2WZn87AXuBvZRCfR3J1XeR+ULAeBJYIek5ZJuBtYAzy1moy07rjasMzExUXObcrnsvX+zOuqe7ZN4XNIK4A2gPyJek3Q/8FVJy4DfAEWAiBiVdAA4TmV4qN9n+lgzZoZ1ZvbuZ4Z1AAqFAr29vZTL5ZrbVtczsysp4orh+Guur68vhoeHW90MazP5fL5muOdyOcbHx6/4cpirnlmnknQ0Ivqa2bbRPX+za26uYZ2Z8pm9+nvuuWde25uZL+9gbay3t7dueaFQIJfLzWv7duHTVK2VHP7WtgYHB+nu7p5V1t3dzeDgYFP12olPU7WWi4iWP+64444wq+XRRx+NXC4XkiKXy8Wjjz66oHrtIpfLBZX1L7MeuVyu1U2zFAGGo8nc9YSvWQt0dXVR6/+eJKanp1vQIkujhUz4etjHrAUamc8wW0oOf7MWSOM8hXUWh79ZCxQKBYaGhsjlckgil8sxNDTkRWl2zXjM38wspTzmb2Zm8+LwN2sRL/KyVvLlHcxaoN5F68yWmvf8zVrA9yKwVnP4m7VAvYvWmS01h79ZC3iRl7Waw9+sBbzIy1rN4W/WAl7kZa3mRV5mZinlRV5mZjYvDn8zswxy+JuZLVAaV2s3FP6SdkoakTQqaVdS9pik55PHuKTnq+rvljQm6aSkLUvTdDOz1kvrLTnrTvhKuhXYD2wCXgcOAg9ExKmqOl8Gfh0Rn5e0HtiX1L8ReBpYGxEX5vodnvA1s7TK5/OUy+UrynO5HOPj40v6u5d6wncdcCQipiLiPHAY2F71ywV8hErgA2wD9kfEuYg4DYxR+SIwM+s4aV2t3Uj4jwCbJa2Q1A1sBVZXvX4n8HLVkcBNwJmq188mZbNIKkoaljQ8OTnZXOvNzFosrau164Z/RJwA9gCHqAz5HAPOV1X5KJf2+gFU621qvO9QRPRFRF9PT8+8Gm1m1i7Sulq7oQnfiNgbERsjYjPwKnAKQNIy4G7gsarqZ5l9ZLAKeHFxmmtm1l7Sulq7oRW+klZGxCuSeoEfA++KiNck3QXsjoh3V9XdAHyfSxO+zwBrPOFrZra4FjLh2+jNXB6XtAJ4A+iPiNeS8h3MHvIhIkYlHQCOUxke6r9a8JuZ2bXX6LDPnRGxPiJui4hnqsr/LCK+UaP+YES8LSJuiYinFrPBZmlcUGPWbnwbR0sV3/7QbHH48g6WKr79oXWKVh/Bes/fUiWtC2rMqrXDEaz3/C1V0rqgxqxaOxzBOvwtVdK6oMasWjscwTr8LVXSuqDGrFo7HME6/C11CoUC4+PjTE9PMz4+7uC31GmHI1iHv5nZNdYOR7C+gbuZWUr5Bu5mZjYvDn8zswxy+LexVq8ANLPO5RW+baodVgCaWefynn+baocVgGbWuRz+baodVgCaWedy+LepdlgBaGady+HfptphBaCZdS6Hf5tqhxWAZta5vMLXzCylvMLXzMzmxeFvZpZBDYW/pJ2SRiSNStpVVf5JSSeT8oeqyndLGkte27IE7TYzswWou8JX0q3A/cAm4HXgoKQfAauAbcDbI+KcpJVJ/fXADmADcCPwtKS1EXFhifpgZmbz1Mie/zrgSERMRcR54DCwHXgA+GJEnAOIiFeS+tuA/RFxLiJOA2NUvjjMzKxNNBL+I8BmSSskdQNbgdXAWuBOST+VdFjSO5L6NwFnqrY/m5TNIqkoaVjS8OTk5MJ6YWZm81I3/CPiBLAHOAQcBI4B56kMGb0FeCfwaeCAJAGq9TY13ncoIvoioq+np6f5HpjZkvLVZTtTQxO+EbE3IjZGxGbgVeAUlT36J6LiOWAauCEpX121+SrgxcVttpldCzNXly2Xy0TExavL+gsg/Ro922dmMrcXuBvYB/wz8L6kfC3wW8CvgCeBHZKWS7oZWAM8t+gtN7Ml56vLdq5Gz/N/XNJx4IdAf0S8BjwC/L6kEWA/cG9yFDAKHACOUxkm6veZPjZfHmpoD766bOdq6GYuEXFnjbLXgXvmqD8I+Apk1hTfyKZ99Pb2Ui6Xa5ZbunmFr7UdDzW0D19dtnM5/K3teKihffjqsp3LV/W0tpPP52sONeRyOcbHx699g8zalK/qaR3FQw1mS8/hb23HQw1mS8/h36HSfqpkoVBgfHyc6elpxsfHHfwpkvbPXlY0dKqnpYtPlbRW8WcvPTzh24E8YWqt4s/eteUJX5vFp0paq/izlx4O/w401+pLr8q0pebPXno4/DuQT5W0VvFnLz0c/h3Ip0paq/izlx6e8DUzSylP+JqZ2bw4/M3MMsjhb2aWQQ7/DuTl9WZWjy/v0GG8vN7MGuE9/w7ju2CZWSMc/h3Gy+vNrBEO/w7j5fVm1oiGwl/STkkjkkYl7UrKPifpl5KeTx5bq+rvljQm6aSkLUvUdqvBy+vNrBF1w1/SrcD9wCbgNuADktYkL38lIm5PHv+W1F8P7AA2AHcBX5N03ZK03q7g5fVm1ohGzvZZBxyJiCkASYeB7Vepvw3YHxHngNOSxqh8cTy70MZaYwqFgsPezK6qkWGfEWCzpBWSuoGtwOrktU9I+pmkRyS9JSm7CThTtf3ZpGwWSUVJw5KGJycnF9AFMzObr7rhHxEngD3AIeAgcAw4D3wdeBtwO/AS8OVkE9V6mxrvOxQRfRHR19PT01TjzcysOQ1N+EbE3ojYGBGbgVeBUxHxckRciIhp4JtUhnagsqe/umrzVcCLi9loSxevODZrP42e7bMy+dkL3A3sk/TWqirbqQwPATwJ7JC0XNLNwBrgucVrsqXJzIrjcrlMRFxccewvALPWauh6/pL+A1gBvAH8RUQ8I+l7VIZ8AhgHPhYRLyX1B4D7qAwP7YqIp672/r6ef+fyDb3Nls5Crufvm7nYkurq6qLWZ0wS09PTLWiRWefwzVysbXnFsVl7cvjbkvKKY7P25PC3JeUVx2btyWP+ZmYp5TF/MzObF4e/mVkGOfzNzDLI4W9mlkEOfzOzDHL4m5llkMPfzCyDHP5mZhnk8DczyyCHv5lZBjn8zcwyyOFvZpZBDn8zswxy+JuZZZDD38wsg1Id/qVSiXw+T1dXF/l8nlKp1OommZmlwrJWN6BZpVKJYrHI1NQUAOVymWKxCOC7RJmZ1dHQnr+knZJGJI1K2nXZa5+SFJJuqCrbLWlM0klJWxa5zQAMDAxcDP4ZU1NTDAwMLMWvMzPrKHX3/CXdCtwPbAJeBw5K+lFEnJK0GvgjYKKq/npgB7ABuBF4WtLaiLiwmA2fmJiYV7mZmV3SyJ7/OuBIRExFxHngMLA9ee0rwF8C1TcC3gbsj4hzEXEaGKPyxbGoent751VuZmaXNBL+I8BmSSskdQNbgdWSPgj8MiKOXVb/JuBM1Z/PJmWzSCpKGpY0PDk5Oe+GDw4O0t3dPausu7ubwcHBeb+XmVnW1A3/iDgB7AEOAQeBY8B5YAD4bI1NVOttarzvUET0RURfT0/PvBoNlUndoaEhcrkcksjlcgwNDXmy18ysAQ2d7RMRe4G9AJL+FngZKADHJAGsAv5L0iYqe/qrqzZfBby4iG2+qFAoOOzNzJrQ6Nk+K5OfvcDdwHcjYmVE5CMiTyXwN0bE/wBPAjskLZd0M7AGeG5JWm9mZk1p9Dz/xyWtAN4A+iPitbkqRsSopAPAcSrDQ/2LfaaPmZktTKPDPnfWeT1/2Z8HAc+8mpm1qVRf3sHMzJrj8DczyyCHv5lZBjn8zcwyyOFvZpZBDn8zswxy+JuZZZDD38wsgxz+ZmYZ5PA3M8sgh7+ZWQY5/M3MMsjhb2aWQQ5/M7MMcvibmWWQw9/MLIMc/mZmGdQx4V8qlcjn83R1dZHP5ymVSq1ukplZ22r0Hr5trVQqUSwWmZqaAqBcLlMsFgEoFAqtbJqZWVvqiD3/gYGBi8E/Y2pqioGBgRa1yMysvXVE+E9MTMyr3Mws6xoKf0k7JY1IGpW0Kyn7gqSfSXpe0o8l3VhVf7ekMUknJW1ZorZf1NvbO6/yuXjewMyyom74S7oVuB/YBNwGfEDSGuDhiHh7RNwO/Cvw2aT+emAHsAG4C/iapOuWpvkVg4ODdHd3zyrr7u5mcHCw4UCfmTcol8tExMV5A38BmFknamTPfx1wJCKmIuI8cBjYHhH/W1Xnd4BInm8D9kfEuYg4DYxR+eJYMoVCgaGhIXK5HJLI5XIMDQ0BNBzonjcwsyxRRFy9grQO+BfgXcD/Ac8AwxHxSUmDwJ8CvwbeGxGTkv6eypfFo8n2e4GnIuIfL3vfIlAE6O3tvaNcLi9uz4B8Pk+t983lcoyPj88q6+rqotbfhSSmp6cXvW1mZgsl6WhE9DWzbd09/4g4AewBDgEHgWPA+eS1gYhYDZSAT8y0p9bb1HjfoYjoi4i+np6eZtpe13wmghdr3sDMLA0amvCNiL0RsTEiNgOvAqcuq/J94EPJ87PA6qrXVgEvLrShzZhPoF9t3sDMrNM0erbPyuRnL3A3sC+Z9J3xQeCF5PmTwA5JyyXdDKwBnlu8JjduPoE+17yBF4mZWSdqdIXv45JWAG8A/RHxmqRvSboFmAbKwMcBImJU0gHgOJXhof6IuLAEba9rJrgHBgaYmJigt7eXwcHBOQO9UCg47M0sE+pO+F4LfX19MTw83OpmmJmlypJO+JqZWedx+JuZZZDD38wsgxz+ZmYZ5PA3M8ugtjjbR9IkldNF28ENwK9a3YhrxH3tTFnqK2Srv5f3NRcRTV0ioS3Cv51IGm721Km0cV87U5b6Ctnq72L21cM+ZmYZ5PA3M8sgh/+VhlrdgGvIfe1MWeorZKu/i9ZXj/mbmWWQ9/zNzDLI4W9mlkGZC39JOyWNSBqVtCsp+5ykX0p6Pnlsraq/W9KYpJOStrSs4U2o1dek/JNJf0YlPVRV3lF9lfRY1b/puKTnq+p3Wl9vl3Qk6euwpE1V9VPbV5izv7dJelbSzyX9UNL1VfVT019Jj0h6RdJIVdnvSjok6VTy8y1Vr9Xsm6Q7kr+LMUl/J6nWHRVni4jMPIBbgRGgm8q9DJ6mcrOZzwGfqlF/PZXbVi4HbgZ+AVzX6n4ssK/vTZ4vT+qt7NS+Xlbny8BnO7WvwI+BP07qbAX+Pe19rdPf/wTendS5D/hCGvsLbAY2AiNVZQ8Bn0mefwbYU69vVG6Y9S4qt9F9auazcLVH1vb811G5ufxURJwHDgPbr1J/G7A/Is5FxGlgDNh0lfrtZK6+PgB8MSLOAUTEK0n9TuwrAMle0EeAfUlRJ/Y1gJm93zdz6dapae4rzN3fW4CfJHUOcek2sqnqb0T8hMqtcattA76TPP8O8CdV5Vf0TdJbgesj4tmofBN8t2qbOWUt/EeAzZJWSOqmsoc0c7/hT0j6WXIYNnOYdRNwpmr7s0lZGszV17XAnZJ+KumwpHck9TuxrzPuBF6OiJl7T3diX3cBD0s6A3wJ2J3UT3NfYe7+jlC5fSzAh7n07532/gL8XkS8BJD8XJmUz9W3m5Lnl5dfVabCPyJOAHuo7CkcpHIIdR74OvA24HbgJSpDBFA5hLribZa8oYvgKn1dBrwFeCfwaeBAsmfciX2d8VEu7fVDZ/b1AeDBiFgNPAjsTTZJbV/hqv29D+iXdBR4E/B6skmq+1vHXH1rqs+ZCn+AiNgbERsjYjOVw61TEfFyRFyIiGngm1w6TDzL7D3IVVw6nG57tfpKpU9PRMVzVO7BfAOd2VckLQPuBh6rqt6Jfb0XeCKp8gM65DMMc/6ffSEi3h8Rd1D5Yv9FUj31/QVeToZySH7ODM3O1bezyfPLy6+u1RMe1/rBpQnOXuAFKnvBb616/UEq42oAG5g9wfLftPHkUYN9/Tjw+aR8LZXDSHViX5M/3wUcvqxux/UVOAG8Jyn/Q+BoJ/T1Kv2dKeuiMsZ9X1r7C+SZPeH7MLMnfB+q1zcqE+Dv5NKE79Z6v3dZrS+EDve4pBXAG0B/RLwm6XuSbqdyqDQOfAwgIkYlHQCOUznU7I+IC61pdlNq9fUR4JHk1LLXgXuj8unpuL4m5TuYPeTTqf+u9wNfTY50fgMUoSP6CrX7u1NSf/L6E8C3IX39lbQPeA9wg6SzwF8DX6QyHPvnwASVOY16fXsA+Afgt6mE/1N1f3fyrWFmZhmSuTF/MzNz+JuZZZLD38wsgxz+ZmYZ5PA3M8sgh7+ZWQY5/M3MMuj/AUIojLCq/ex+AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "T33=pt.calculate_cpx_liq_temp(liq_comps=myLiquids1, cpx_comps=myCpxs1, equationT=\"T_Put2008_eq33\", P=my_input['P'])-273.15\n", "plt.plot(T, my_input['temp_out'], 'ok')" ] }, { "cell_type": "code", "execution_count": 33, "id": "a0d09f26-3e04-4042-b6ee-9f7756d66e8c", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'equation 34 (liq-only sat temp)')" ] }, "execution_count": 33, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEGCAYAAACUzrmNAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAAAxPElEQVR4nO3deXxU1d3H8c8viGKKWGWxdSFRK1VA5QFcW2nVWi2P1rVqi48LaiyixaVaqT7V1uKGK+5BcCOCGwpaRQUVVFAMsm8VZZHHDRURTGXL7/nj3MgAmcxkmdyZyff9et1XZs7cmfs7mUl+c+655xxzd0RERGpSEHcAIiKS/ZQsREQkJSULERFJSclCRERSUrIQEZGUtog7gExp06aNFxcXxx2GiEhOmTJlyhfu3nbT8rxNFsXFxZSXl8cdhohITjGzxdWV6zSUiIikpGQhIiIpKVmIiEhKShYiIpKSkoWIiKSkZCEi0ojKysooLi6moKCA4uJiysrKGnT/TMnbS2dFRLJNWVkZJSUlVFRUALB48WJKSkoA6NWrV733zyTL1ynKu3fv7hpnISLZpLi4mMWLNx/GUFRUxKJFi+q9f0Mwsynu3n3Tcp2GEhFpJEuWLMloeSYpWYiINJL27dtntDyTlCxERBrJgAEDKCws3KissLCQAQMGNMj+maRkISLSSHr16kVpaSlFRUWYGUVFRZSWlibtrK7t/pmkDm4REfmeOrhFRKTOlCxERCQlJQsREUkpY8nCzIaa2edmNiuhbHsze8XM3o9+bpfwWH8zW2Bm883syITybmY2M3pskJlZpmIWEZHqZbJl8RBw1CZlVwDj3H0PYFx0HzPrCJwKdIqec4+ZNYuecy9QAuwRbZu+poiIZFjGkoW7TwC+2qT4WODh6PbDwHEJ5SPcfbW7LwQWAPub2Y+BVu4+ycNlW48kPEdERBJ9/TWMGJGRl27sPosd3P0TgOhnu6h8J+CjhP2WRmU7Rbc3La+WmZWYWbmZlS9btqxBAxcRyWqjRkHHjnD66fDRR6n3r6Vs6eCurh/CayivlruXunt3d+/etm3bBgtORCRrff45nHoqHHcctG0LEyfCLrs0+GEaO1l8Fp1aIvr5eVS+FEis3c7Ax1H5ztWUi4g0be5QVhZaE888A9deC+++C903G0/XIBo7WYwGzohunwGMSig/1cy2MrNdCR3Zk6NTVSvN7MDoKqjTE54jItI0LV0KxxwDp50GP/kJTJ0KV10FW26ZsUNm8tLZ4cAk4KdmttTMzgZuAI4ws/eBI6L7uPts4AlgDjAG6Ovu66OX6gM8QOj0/gB4MVMxi4hktcpKuP/+0Jp49VW49VZ4661wP8M0N5SISC5YsADOPRdefx0OOwwGD4bddmvww2huKBGRXLR+PdxyC+yzD7z3XkgSY8dmJFHURGtwi4hkq1mz4OyzYfLk0Edx772wU9LRAxmlloWISLZZswb+/nfo2hU+/BCGDw/jKGJKFKCWhYhIdnn3XejdO7Qqfv97uOOOMH4iZmpZiIhkg4oKuOwyOPBA+OorGD0aHnssKxIFqGUhIhK/8ePhnHM2XPE0cCBsu23cUW1ELQsRkbh88w306QO//GUYQzFuHJSWZl2igDSShZm1M7PjzayvmfU2s/3NTElGRKpVVlZGcXExBQUFFBcXU1ZW1qD7540XXoDOncMgu4svhhkzwviJLJX0NJSZHUpYb2J7YCphHqcWhCnCdzezp4Bb3P2bRohTRHJAWVkZJSUlVFRUALB48WJKSkoA6NWrV733zwtffgkXXQTDhoWR1xMnhn6KLJd0BLeZDQTudPcl1Ty2BXA00Mzdn85siHWjEdwija+4uJjFixdvVl5UVMSiRYvqvX9Oc4cnn4QLLoDly6F/f7jySthqq7gj20iyEdya7kNEGkxBQQHV/U8xMyorK+u9f8765BM4/3x49lno1g2GDg0jsrNQnaf7MLPWZnanmb1nZlPM7A4za52ZMEUkl7Vv3z6j5TnHHR58MJxuGjMGbrwR3n47axNFTdLpqB5B6K84ETgJWAY8nsmgRCQ3DRgwgMLCwo3KCgsLGTBgQIPsn1MWLYIjjwwD7PbeG6ZPh8svhy1ydMSCu9e4AVOqKStP9by4t27durmINL5hw4Z5UVGRm5kXFRX5sGHDGnT/rLd+vfsdd7j/4AfuLVu63313KMsRyf6/p+yzMLObgXLCehMQWhed3P3qzKWw+lOfhYg0unnzwsR/EyfCUUfBffdBUVHcUdVKfaYoPw94DFgTbSOAS8xspZnpslkRkbVr4frroUsXmDsXHn44jKPIsURRk5Qnz9x9m8YIREQkJ02dGloTU6fCSSfBXXfBDjvEHVWDS6unxcz2AYoT93f3kRmKSUQk+333HfzjH3DTTdCmDTz9NJxwQtxRZUzKZGFmQ4F9gNlA1YXPDihZiEjTNHFiaE3MmwdnnhlWstt++7ijyqh0WhYHunvmVwMXEcl2q1aFUdd33gm77BLGThx5ZNxRNYp0OrgnmZmShYg0ba+8EsZLDBoEffuGxYmaSKKA9FoWDxMSxqfAasAAd/fcG4IoIlJby5fDpZeGkdgdOsAbb8DPfx53VI0unWQxFPgfYCYb+ixERPLfs8+G9SaWLYMrroC//Q223jruqGKRTrJY4u6jMx6JiEi2+OwzuPDCMEvsvvvC88+HCQCbsHSSxTwzewx4jnAaCtClsyKSh9yhrAz69Qud2f/8Z5jPqXnzuCOLXTrJYmtCkvh1QpkunRWR/PLRR3DeefDii2ExoiFDwmyxAqQ3gvusxghERCQWlZVh3evLL4f16+H228MCRc2axR1ZVklnPYsOZjbOzGZF9/cxs6syH5qISIa9/35Y97pPH9h/f5g5M5yCUqLYTDrjLAYD/YG1AO4+Azg1k0GJiGTUunVw881hEaJp0+CBB8I4it12izuyrJVOsih098mblK2rz0HNrJ+ZzTKz2WZ2UVS2r5lNMrOZZvacmbVK2L+/mS0ws/lm1nRGwYhIw5s5Ew4+GC67DH79a5gzJ0zdYRZ3ZFktnWTxhZntTujUxsxOAj6p6wHNrDNwLrA/sC9wtJntATwAXOHuewPPAJdF+3cktGQ6AUcB95iZ2ogiUjtr1sDVV0PXrmEVuxEjwjiKHXeMO7KckE6y6AvcD+xpZv8HXAT8sR7H3At4290r3H0dMB44HvgpMCHa5xXCMq4AxwIj3H21uy8EFhASjYhIeiZPDkniH/+Ak08OrYlTTlFrohbSSRbu7r8C2gJ7uvvP03xeMrOAHmbW2swKgZ7ALlH5b6N9fheVAewEfJTw/KVR2WbMrMTMys2sfNmyZfUIUUTyQkUF/PnPcNBB8PXX8NxzYRxFmzZxR5Zz0vmn/zSAu3/r7iujsqfqekB3nwvcSGg9jAGmE/pAegN9zWwKsA1hVT4Ic1Ft9jJJXrvU3bu7e/e2bdvWNUQRyQevvx46sG+5Bc49F2bPhqOPjjuqnJV0nIWZ7UnoJ9jWzBJX9GgFtKjPQd19CDAkOs51wFJ3n0c08M/MOgD/He2+lA2tDICdgY/rc3wRyWMrVoQxE6WlsPvu8OqrcOihcUeV82oalPdT4Gjgh8AxCeUrCR3UdWZm7dz9czNrD5wAHJRQVgBcBdwX7T4aeMzMbgV2BPYANr06S0QE/vWvMAr7k0/gkkvg2muhsDDuqPJC0mTh7qOAUWZ2kLtPauDjPm1mrQljN/q6+/Loctq+0eMjgQejOGab2RPAHMLpqr7uvr6B4xGRXPbFF3DRRaE/olOnsMTpAQfEHVVeMfdqT//nvO7du3t5eXncYYhIJrnDE0+EGWKXLw+r2P31r7DllnFHlrPMbIq7d9+0PJ2JBEVEss/HH8P558OoUdC9O4wbF1ayk4xIZ26oXdMpExFpFO4bZoR96SUYOBAmTVKiyLC0L53dRJ0vnRURqbOFC+GII+Ccc8KiRDNmhHEUW+gkSabFcumsiEitrF8Pd90V+iOaNYN774WSEiioz/hgqY1YLp0VEUnb3Llhor9Jk+A3v4H774dddkn9PGlQSdOyu4+KFj462t3PStj+5O4TGzFGEcrKyiguLqagoIDi4mLKysriDkkybe1aGDAAunSB+fPh0UfDOIocSBT5+HlN50Tf1Gj8QycSTj+5e++MRSWSoKysjJKSEioqKgBYvHgxJSUlAPTq1SvO0CRT3nsPeveG6dPhd7+DO++EHXaIO6q05OvnNZ0Tfo8CPwKOJMwQuzPhVJRIo7jyyiu//8OrUlFRwZVXXhlTRJIx330H/fuHVes++wxGjgzjKHIkUUD+fl7TSRY/cff/Bb5194cJczbpGjVpNEuWLKlVueSoN98MVzjdcAOcfnqYRvz44+OOqtby9fOaTrJYG/38Olq4aFugOGMRiWyiffv2tSqXHLNyZRiB3aMHrF4NL78MQ4fCdtvFHVmd5OvnNZ1kUWpm2xEm9xtNmKPppoxGJZJgwIABFG4yGVxhYSEDBgyIKSJpMC+/DJ07w913wwUXwKxZYRxFDsvbz6u75+XWrVs3l/wxbNgwLyoqcjPzoqIiHzZsWNwhSX189ZX7mWe6g/tPf+r+5ptxR9SgcvnzCpR7Nf9TU04kaGb9CDPArgQGA10Ja2W/nPFMVg+aSFAkS40cCX37wrJlYd2Jv/0NWmicb7ZINpFgOqeherv7N4SFidoBZwE3NHB8IpLvPv00XAZ74onwox/Bu+/CddcpUeSIdJJF1bKmPYEH3X061S91KiKyOXd45JEw8d/o0WGg3eTJ8F//FXdkUgvpDMqbYmYvA7sC/c1sG6Ays2GJSF5YsiSsXDdmDBx8cJgtds89445K6iCdZHE20AX40N0rohXuzspoVCKS2yor4b774C9/CS2LQYPC2hPNmsUdmdRRymTh7pXAewn3vwS+zGRQIpLD3n8/TCE+YQL86lcweDAUF8cdldST5vcVkYaxbh3cdBPss0+Y02nIkDCOQokiL2jFEBGpvxkzwsR/U6bAcceFQXY77hh3VNKA0llW9WYz69QYwYhIjlm9OoyT6NYtdGY//ngYR6FEkXfSaVnMI0z5sQVhcN5wd1+R2bBEJOu9805oTcyZA6edBrffDq1bxx2VZEjKloW7P+DuPwNOJ0wgOMPMHjOzQzMdnIhkoW+/hUsugYMOgm++CQsSPfqoEkWeS6uD28yaAXtG2xfAdOASMxuRwdhEJNu8+mrowL7ttjB+YvZs6Nkz7qikEaTTZ3ErMJ8wgvs6d+/m7je6+zGAhmAmkY/LKuYzvV8prFgBJSVw+OFgBq+9BvfeC61axR2ZNJbqZhdM3IDeQGGSx7ZN9fy4tjhnnR02bJgXFhY68P1WWFiYUzNPNiV6v1J47jn3nXZyLyhw//Of3b/9Nu6IJINIMuts0paFmXU1s67ANGDPqvsJ5XgT7OhO5xtovi6rmIv0ftXDsmXwhz/AMceEhYjefhsGDoRN1mqQJqK6DBKSC6/VsL2a7HnZsmWiZZHuN1Az22ifqs3MGjwmSU7vVx1VVroPH+7epo178+bu11zjvnp13FFJIyFJyyKWf+RAP2AWMBu4KCrrArxNaMmUA/sn7N8fWEDoOzkynWNkIlkUFRVV+0+lqKioTvtJZun9qoOlS92POSb8a9hvP/eZM+OOSBpZvZIFcDDwB8Lls6cDp6fzvCSv1TlKFIWEcR5jgT2Al4HfRPv0BF6PbnckXH21FWHm2w+AZqmOk4lkke43UJ0Dzw56v2qhstK9tNS9VSv3Fi3cb77Zfd26uKOSGCRLFulcDfUocDPwc2C/aNtsFaVa2At4290r3H0dMB44Pvojrbq0Ylvg4+j2scAId1/t7gsJLYz963H8Okt3IfZevXpRWlpKUVERZkZRURGlpaX06tWrMcIUQl9FQUH1H+9cf78a/MqtDz8ME/6VlIQ1JmbOhEsv1QyxsrHqMohv3BKYC2H51YbYCMni30BrQutiEnBnVL4E+Aj4P6Ao2v8u4LSE5w8BTkry2iWEU1jl7du3b/CMq2+guaG69ylf3q8G/QyuW+d+223uhYXu22zjft997uvXN3jMkluo62ko4Engx6n2q81GWCPjPWACcB9wGzAIODF6/GRgbHT77mqSxYmpjpGpS2dzeSH2piJZH0SzZs1y/v1qsP6V2bPdDzww/Avo2dN9yZKMxCu5J1mysPBYcmb2GqHzeTKwOqFF8tsan5gmM7sOWApcD/zQ3d3MDFjh7q3MrH90vOuj/V8CrnH3STW9bvfu3b28vLwhQpQcU1BQQHWfazOjsjK3F3msd93WroUbboB//hO22QbuuCNcHmtaKVkCM5vi7pt1NaQzkeA1GQimnbt/bmbtgROAg4ALgV8ArwOHAe9Hu48GHotGku9I6Ayf3NAxSf5o3749ixcvrrY819WrblOmhIn/ZsyAU04Jq9e1a5eBKCUfpTOR4HjCzLPbRNvcqKw+njazOcBzQF93Xw6cC9xiZtOB6wj9D7j7bOAJYA4wJtp/fT2PL3lswIABFG4ycKywsJABAwbEFFHDqVPd/vMfuOIKOOCAMNDu2WdhxAglCqmd6s5N+cb9CycDi4GHgUeAhSTpYM6mLc7pPiR++dy3VKu6TZjg3qFD6Js4+2z35csbLU7JTdSjz2I6cIS7fx7db0vofN43gzms3tRnIU3aypWhNXHPPWFZ08GDw+WxIikk67NIZ4rygqpEEfkyzedJHsj32Vjzsn4vvQSdO4dZYfv1C+MmlCikntLp4B4TXYE0PLp/CvBC5kKSbFFWVkZJScn3k+wtXryYkpISgKwdsFYbeVe/r74KixI9/DDsuSe8+SYcfHDcUUmeSHkaCsDMTiCM4DZggrs/k+nA6kunoeqvuLi42itvioqKWLRoUeMH1MDyqn5PPw19+8IXX4TTT1ddBS1axB2V5KBkp6HSShYJL9LV3d9r0MgyRMmi/vJ5vALkSf0+/TQkiZEjw1QdQ4dCly5xRyU5rD59FokeaKB4JAekOxdWrsrp+rmH000dO4Y1sK+/Ht55R4lCMqa2yULDPJuQfB6vADlcv8WL4Te/gTPPDMli2rRw6ql587gjkzxW22Tx94xEIVkp12Zjra2cq19lJdx9d7jS6c034c47YcKE0JktkmHpjLPoWtPj2dqHoT4LySvz58M554QkccQRUFoaxk+INLD6zA11D9AVmEE4DbU3YW6mtYQZLw9rwDhFJNG6dXDLLXD11bD11vDgg3DGGZr4TxpdOqehFgHd3L27u3cjJI4F7n6ouytR5LG8HLBGDtVr+vQwn9MVV0DPnjBnTuinUKKQOFQ3B0jiBkxLpyzbNs0NVT/5utBTTtTru+/cr7rKfYst3Nu1c3/yybgjkiaEeswNNRz4FhgW/XGdBrR0999nKH81CPVZ1E9eDVhLkPX1mjQJzj4b5s6F00+HW2+F1q3jjkqakPqMszgLmA30Ay4iTBV+VoNGJ1lnyZIltSrPFVlbr2+/hYsugp/9DFatghdeCOMolCgkS6Ts4Hb37wjLnt6W+XAkW+TrAkJZWa9x4+Dcc2HhQujTJ6xk16pVfPGIVCNpy8LMnoh+zjSzGZtujReixCFnB6ylkFX1+vrrkCR+9SvYYgsYPz5MKa5EIdmouo6MqB/jx9HPouq2ZM/Llk0d3PWXrwsIZUW9Ro1y33FH94IC98svd6+oaPwYRKpBXTu4c5U6uCUrLVsGf/pTWNZ0773DxH/dN+tLFIlNrQflmdlKwtVPmz0EuLurrSySLncYPjwkim++gb//PYyf2HLLuCMTSUvSZOHu2zRmICJ5a+nS0HH9/PNhkN2QIdCpU9xRidRKTR3cLVM9OZ19RJos9zCHU6dO4YqnW2+Ft95SopCcVNM4i1FmdouZ9TCzH1QVmtluZnZ2tNTqUZkPUSQHffABHH44nHcedOsW1sG++GJo1izuyETqJGmycPfDgXHAecBsM1thZl8SRnL/CDjD3Z9qnDAbX87MHyTZZf360ILYe2+YMiW0LMaNg913b9DD6PMpja66S6TyYavPpbM5MX+QZJ9Zs9z3398d3I8+2v2jjzJyGH0+JZPQpbPpy/r5gyS7rFkTRl3/85+w7bYwaBCcemrGZofV51MyqT7rWTQ5WTt/kGSf8nLo3Tv0SZx6akgUbdtm9JD6fEocarusapOQbJ6gXJ8XSRrQf/4Dl18eLoX98ksYNSqMo8hwogB9PiUetUoWZvbbTAWSTbJq/iDJPhMmwD77wMCBoVUxezb8tvH+NPT5lFhU15ER9WOcsMl2IvBp1f1kz0tnI0x3Posw9flFUdnjwLRoW0TCAktAf2ABMB84Mp1j1HduqKyYP0iyy4oV7n36hA7sXXd1Hzs2tlD0+ZRMobYd3Ga2DhgDfE6Y4gPgJOCpkGO8d12Sk5l1BkYA+wNromP0cff3E/a5BVjh7v8ws47A8Gj/HYGxQAd3X1/TcTQ3lDSoF18MYyaWLoV+/UJn9g9+kPp5IjmmLosfHQRsDbwL9Hb3s4Av3P2suiaKyF7A2+5e4e7rgPHA8QmBGnAyIUEAHAuMcPfV7r6Q0MLYvx7HF0nfl1+GFet69oSWLWHiRLjtNiUKaXJqGpT3LnAEsCXwqpntT/UTC9bWLKCHmbU2s0KgJ7BLwuOHAJ8ltDR2Aj5KeHxpVLYZMysxs3IzK1+2bFkDhCpNljs8+SR07Bg6rq+6CqZOhQMPjDsykVjUeOmsu1cCd5jZk8DtDXFAd59rZjcCrwCrgOnAuoRdfs+GVgVsOAW20cskee1SoBTCaaiGiFeaoE8+gb594ZlnoGtXePll2HffuKMSiVVaV0O5+8fufrK779YQB3X3Ie7e1d17AF8B7wOY2RaEDvTHE3ZfysYtj52BjxsiDpGNuMODD4bWxAsvhIF277yjRCFCzbPOFprZ5WZ2mZm1MLMzzGy0md1U39lmzaxd9LM9ITlUtSR+Bcxz96UJu48GTjWzrcxsV2APYHJ9ji+ymUWL4KijwqWwnTvDjBnwl7+E5U5FpMaWxUPADsCuwL+A/YCbCaeF7q3ncZ82sznAc0Bfd18elZ/KxqegcPfZwBPAHMKVU31TXQklkrbKSrjzzpAgJk6Eu+4Ka2F36BB3ZCJZpaZLZ6e5e5fo6qRPCGtye3R/urvv05iB1pYunZWU5s+Hs88Oa0wceSTcfz8UFcUdlUis6nLpLBAGVAAvRD+r7qvzWHLX2rVw/fWhL2LOHHjooTCOQolCJKmaTsiWm1lLd1+VOK7CzHYHVmY+NJEMmDYt9EtMnQonnhhOO/3oR3FHJZL1ahpncY67r6qm/APCWAiR3PHdd3DlldC9O3z8MTz1VNiUKETSUqdLPTxZR4dINpo4MfRNzJsHZ5wRVrLbfvu4oxLJKZqiXPLXqlVhHqef/xwqKkK/xEMPKVGI1IGSRYJ8Xtc4n+tWrbFjwzrYgwbB+efDrFlhHIWI1Elap6HMbCegKHF/d5+QqaDiUFZWRklJCRUVFQAsXryYkpISAHr16hVnaPWWz3XbzNdfw6WXwtChsMceYe2JQ9TFJlJfKdfgjuZxOoUwKK5qMJy7e1YvhFTbcRb5vK5xPtdtI6NGQZ8+8Pnn8Oc/w9VXw9Zbxx2VSE6pzxrcxwE/dffVDR5VFsnndY3zuW5ASA4XXghPPBFWsHvuOejWLe6oRPJKOn0WHwLNMx1I3PJ5XeO8rZs7lJWFif+efRauvRbKy5UoRDIgnWRRAUwzs/vNbFDVlunAGls+r2ucl3X76CM45hg47bTQNzF1alhzonnef68RiUU6yWI0cC0wEZiSsOWVXr16UVpaSlFREWZGUVERpaWledEBnFd1q6wMczh16gSvvRZWrXvzzdC6EJGMSdnBDWBmWwJV03DOd/e1GY2qAWgiwTy0YAGcc06YFfaww2DwYNitQZZYEZFInScSNLNfEhYnuhu4B/i3mfVo6ABFklq/Hm6+OYybmDo1JImxY5UoRBpROldD3QL82t3nA5hZB8KaE+pFlMybNStM/Pfuu6GP4t57Yadql2AXkQxKp8+ieVWiAHD3f9MEro6SmK1ZA3//e1gDe+FCGD48jKNQohCJRToti3IzGwI8Gt3vRR52cEsWeffd0JqYNQv+8Ae4/XZo2zbuqESatHRaFn2A2cCfgH6Ekdx/zGRQ0kRVVISR1wceCMuXw+jRYRyFEoVI7FK2LKKR27dGm0hmjB8frnRasABKSuCmm2DbbeOOSkQiSZOFmT3h7ieb2UyqWUY129fglhzxzTfwl7/AffeFq5tefRUOPTTuqERkEzW1LPpFP49ujECkCXrhBTjvvLBy3SWXhOk6NhlpLiLZoaZlVT+Jbp7v7osTN+D8xglP8tIXX4RpOv77v6FVq7CS3S23KFGIZLF0OriPqKbsNw0diDQB7mFm2I4d4fHH4X//F957Dw44IO7IRCSFmvos+hBaELuZ2YyEh7YB3sp0YJJnPv4Y+vYNs8N26xZGYO+jbi+RXFFTn8VjwIvA9cAVCeUr3f2rjEYl+cMdHnww9EmsXh2ucrr4YtgirUUaRSRLJP2LdfcVwArg9wBm1g5oAbQ0s5bunicr50jGLFwYLoMdOzYsbfrAA9ChQ+rniUjWSWciwWPM7H1gITAeWERocYhUb/16GDQIOneGt9+Ge+6B119XohDJYel0cP8TOBD4t7vvChyO+iwkmblzoUcP6Ncv/Jw9O6yLXZDOR01EslU6f8Fr3f1LoMDMCtz9NaBLfQ5qZv3MbJaZzTazixLKLzSz+VH5TQnl/c1sQfTYkfU5tmTI2rVw3XXQpQvMmwePPBLGUeT60q0iAqQ3keDXZtYSmACUmdnnwLq6HtDMOgPnAvsDa4AxZvYvYGfgWGAfd18d9ZFgZh2BU4FOwI7AWDPr4O7r6xqDNLCpU8PEf9OmwUknwV13wQ47xB2ViDSgdFoWxwL/AS4GxgAfAMfU45h7AW+7e4W7ryP0gxxPmLDwhmguKtz984Tjj3D31e6+EFhASDQSt+++g7/+FfbbDz79FJ5+Gp58UolCJA+lTBbu/q27r3f3de7+sLsPik5L1dUsoIeZtTazQqAnsAth2dZDzOwdMxtvZvtF++8EfJTw/KVRmcTprbfCKafrr4f/+R+YMwdOOCHuqEQkQ1KehjKzlWyYSHBLwsJH37p7q7oc0N3nmtmNwCvAKmA64bTWFsB2hM70/YAnzGw3wKp7mSSxlgAlAO11rjwzVq0KrYm77gr9ES+9BL/+ddxRiUiGpdOy2MbdW0VbC+BE4K76HNTdh7h7V3fvAXxFWON7KTDSg8lAJdAmKt8l4ek7Ax8ned1Sd+/u7t3bag2Ehvfyy+Fy2LvuggsuCIsTKVGINAm1vp7R3Z8FDqvPQRM6r9sDJxDW9P7+daN1vrcEvgBGA6ea2VZmtiuwBzC5PseXWlq+HM46C448Elq0gAkTwjiKli3jjkxEGkk6p6EST0QXAN1JchqoFp42s9bAWqCvuy83s6HAUDObRbhK6gx3d2C2mT1BWKFvXbS/roRqLM88A+efD8uWwRVXwNVXh4QhIk1KOpfOJl75tI4wgvvY+hzU3Q+ppmwNcFqS/QcAA+pzTKmlzz6DCy8MVzd16QL/+hd07Rp3VCISk3SWVT2rMQKRLOEOw4bBRReFzuwBA+Cyy6B587gjE5EYpXMaalBNj7v7nxouHInVkiXwxz/Ciy/CQQfBkCGw115xRyUiWSCdDu4WQFfCFUvvE6b6WA9MiTbJdZWVcO+90KkTjB8Pt98Ob7yhRCEi30unz2IP4FB3XwtgZvcBL7v7xRmNTBrH++/DOeeEK5wOPxwGD4Zdd407KhHJMum0LHYkrI5XpWVUJrls3ToYODCsVjd9ejjl9MorShQiUq10WhY3AFPN7LXo/i+AazIWkWTejBlw9tlQXg7HHhvWm9hR+V9EkkvnaqgHzexF4ICo6Ap3/zSzYUlGrF4dphG/7jrYbjsYMQJOPhmsuhlVREQ2SGelPAN+Bezr7qOALc1Ms77mmnfegW7d4B//gFNOCRP/nXKKEoWIpCWdPot7gIOI1uIGVgJ3ZywiaVgVFXDppXDwwbBiBTz/fBhH0aZN3JGJSA5Jp8/iAHfvamZTAaKpObbMcFzSEF5/PVzp9MEHcN55cNNN0KpOkwWLSBOX1rKqZtaMaD4oM2tLmBFWstWKFSE5HHpouP/aa3DffUoUIlJn6SSLQcAzQDszGwC8CVyX0aik7p5/Pgyue+CBcPppxgz45S/jjkpEclw6V0OVmdkU4HDCQkTHufvcjEcmtfPFF9CvHzz2WEgWI0fC/roOQUQaRjp9Frj7PGBehmORunCHxx8PM8SuWBGmEP/rX2FLdSuJSMNJK1lIlvr4Y+jTB0aPhv32C6Ow99477qhEJA/VeqU8yQLuoU+iY8ew1OnAgTBxohKFiGSMWha55sMPoaQExo2DHj1C0thjj7ijEpE8p5ZFrli/PkwdvvfeMHlymFL8tdeUKESkUahlkQvmzg0T/02aBD17hjETu+wSd1Qi0oSoZZHN1q4Ny5p26QLz58Ojj4ZxFEoUItLI1LLIVu+9B717h7UmTj4Z7rwT2rWLOyoRaaLUssg2//kPXHFFGFD32WfwzDNhHIUShYjESC2LbPLmm6Fv4t//Dq2Km28O606IiMRMLYtssHIlXHABHHIIrFkTxk4MGaJEISJZQ8kibi+9BJ07h6VN//QnmDkTjjgi7qhERDaiZBGXr76CM8+Eo46CrbeGN96AO+6Ali3jjkxEZDNKFnEYOTJM1TFsWJj0b9o0+NnP4o5KRCQpdXA3pk8/DX0TTz8dxk68+CL813/FHZWISEpqWTQGd3jkkdCaeP55uO66MGWHEoWI5IhYkoWZ9TOzWWY228wuisquMbP/M7Np0dYzYf/+ZrbAzOab2ZFxxFxnS5aEKTrOOAP22iuccurfH5o3jzsyEZG0NfppKDPrDJwL7A+sAcaY2b+ih29z95s32b8jcCrQCdgRGGtmHdx9fSOGXXuVlWEOp7/8JbQsBg2C88+HZs3ijkxEpNbi6LPYC3jb3SsAzGw8cHwN+x8LjHD31cBCM1tASDSTMh5pXa1aFVoTb7wRLoMtLYXi4rijEhGpszhOQ80CephZazMrBHoCVTPjXWBmM8xsqJlVjUjbCfgo4flLo7LNmFmJmZWbWfmyZcsyFX9qP/hBmDp86NAwjkKJQkRyXKMnC3efC9wIvAKMAaYD64B7gd2BLsAnwC3RU6y6l0ny2qXu3t3du7dt27aBI68FszAC+6yzwm0RkRwXSwe3uw9x967u3gP4Cnjf3T9z9/XuXgkMJpxqgtCSSJyTe2fg48aNWESkaYvraqh20c/2wAnAcDP7ccIuxxNOVwGMBk41s63MbFdgD2ByY8YrItLUxTUo72kzaw2sBfq6+3Ize9TMuhBOMS0CzgNw99lm9gQwh3C6qm/WXwklIpJnYkkW7n5INWX/U8P+A4ABGQ1KRESS0ghuERFJSclCRERSUrIQEZGUlCxERCQlc692fFvOM7NlwOKYw2gDfBFzDJmmOuaPplBP1TG1InffbFRz3iaLbGBm5e7ePe44Mkl1zB9NoZ6qY93pNJSIiKSkZCEiIikpWWRWadwBNALVMX80hXqqjnWkPgsREUlJLQsREUlJyUJERFJSsqgHM+tnZrPMbLaZXRSVXWNm/2dm06KtZ8L+/c1sgZnNN7MjYwu8FqqrY1R+YVSP2WZ2U0J5XtTRzB5PeA8Xmdm0hP1zro6QtJ5dzOztqJ7lZrZ/wv45V88kddzXzCaZ2Uwze87MWiXsnxN1jFYP/dzMZiWUbW9mr5jZ+9HP7RIeq7ZeZtYt+j0sMLNBZrVYnc3dtdVhAzoT1twoJMzeO5aw1sY1wJ+r2b8jYVXArYBdgQ+AZnHXo451PDS6vVW0X7t8q+Mm+9wC/C1X65jivXwZ+E20T0/g9VytZw11fBf4RbRPb+DaXKsj0APoCsxKKLsJuCK6fQVwY6p6EdYCOoiwAumLVe99OptaFnW3F/C2u1e4+zpgPGHRpmSOBUa4+2p3XwgsYMNqgNkqWR37ADe4+2oAd/882j+f6ghA9M3rZGB4VJSLdYTk9XSg6pv2tmxYhTIX65msjj8FJkT7vAKcGN3OmTq6+wTCqqKJjgUejm4/DByXUL5ZvaIF5lq5+yQPmeORhOekpGRRd7OAHmbW2swKCd/KqpZ/vcDMZkRNx6qm4U7ARwnPXxqVZbNkdewAHGJm75jZeDPbL9o/n+pY5RDgM3d/P7qfi3WE5PW8CBhoZh8BNwP9o/1zsZ7J6jgL+G20z+/Y8P7mYh0T7eDunwBEP9tF5cnqtVN0e9PytChZ1JG7zwVuJHxTGUNo9q0D7gV2B7oAnxBOYUBo9m32MhkPtB5qqOMWwHbAgcBlwBPRN/B8qmOV37OhVQE5WEeosZ59gIvdfRfgYmBI9JScq2cNdewN9DWzKcA2wJroKTlXxzQlq1e96qtkUQ/uPsTdu7p7D0IT8X13/8zd17t7JTCYDc3apWz8jXVnNjT5s1Z1dSTUZaQHk4FKwuRl+VRHzGwLwhrxjyfsnpN1hKT1PAMYGe3yJHn4eXX3ee7+a3fvRkj8H0S752QdE3wWnVoi+ll1OjhZvZZGtzctT0/cHTe5vLGhY7c9MI/wbfvHCY9fTDh3CNCJjTudPiRLO9PSqOMfgX9E5R0ITV7LpzpG948Cxm+yb07WsYb3ci7wy6j8cGBKLtczSR2rygoI5+l752IdgWI27uAeyMYd3Delqhehs/9ANnRw90z3+LGswZ1Hnjaz1sBaoK+7LzezR82sC6F5twg4D8DdZ5vZE8AcQtO4r7uvjyfsWqmujkOBodFlfGuAMzx8EvOmjlH5qWx8CiqX30eo/r08F7gjakV9B5RATtezujr2M7O+0eMjgQcht+poZsOBXwJtzGwpcDVwA+EU8NnAEkJ/TKp69QEeArYmJIsX044hyjYiIiJJqc9CRERSUrIQEZGUlCxERCQlJQsREUlJyUJERFJSspC8YWY/NLPzE+7vaGZPNdBrDzGz6dE0Lk+ZWcuo/NiorGrW1p8neb6Z2auJM542NDN73cy6R7dfMLMf1uK5x5lZx0zFlnCcm83ssEwfRxqekoXkkx8C3ycLd//Y3U9qoNe+2N33dfd9CNe0XxCVjwP2dfcuhGklHkjy/J7AdHf/poHiqZG793T3r2vxlOMIs5Vm2p2EAWSSY5QspNGY2WlmNjn6Fn6/mTWLys8ys39HkxIONrO7ovKHzOykhOevin62NLNxZvZeNDf/sdEuNwC7R68/0MyKq+b/N7MWZvZgtP9UMzs0Kj/TzEaa2ZhoXYCbqEbVP/loDqytiebUcfdVvmGw0g9IPtdOL2BUQl1Oj1ok083s0YT63mdmb0S/j6Oj8kuigZCY2d4W1msoTPG7XmRmbaLbV1pY12CsmQ03sz9vsu/BhIn2Bka/u92jbYyZTYni2TMhxnvN7DUz+9DMfmFhwsy5ZvZQ4ntlZrdE79E4M2sb/b4WA63N7Ec1xS9ZKO4h7NqaxkaYPvo5oHl0/x7gdODHhG/qbYEtgbeAu6J9HgJOSniNVdHPLQhTLUOYk2oBYfqCYjaeDuH7+8ClwIPR7T2jY7YAziRMh7BtdH8xsEuSOjwIfAa8BhQmlB9PmFriK+CgJM9dDGwT3e4EzAfaRPe3T6jvGMKXuD0Ic/m0iO5PiI5TDvwsyTFeB7pHtxdFv5tuwEzCGg+tot9VdeutbPq7Hke0rgdwAPBqwn4jot/3scA3wN5RjFOALtF+DvSKbv+t6j2N7g8GToz7M6mtdptaFtJYDif843rXwqpzhwO7Ef4Rve7uy9x9DRtP2peMAdeZ2QzCAjc7ATukeM7PgUcB3H0e4Z93h+ixce6+wt2/I0yRUFTdC7j7WcCOhPmUTkkof8bd9yScyrk2yfG3d/eV0e3DgKfc/Yvo+YnrFDzh7pUepkT/ENjTw6SUZ0bxj3f3t1LUNdEhwDMe1nj4Bhid6glRf8zBwJPRe3U/IalXec7dnZCEPnP3mVGMswkJGsLkklXv5TDC77/K54Tfo+QQzQ0ljcWAh929/0aFZseR/NTNOqJTpdHpny2j8l6Elkg3d19rZosI38BTHT+Z1Qm311PD34W7rzezxwlTsz+4yWMTotM3baoSQWJdzKwg+qdqJK/zpuVV9/cAVpHwT9bMXiIkyXJ3PydZzDUcK5kC4GsP/TDVqfp9VbLx766S5L+7xBhaAP+pZUwSM7UspLGMA04ys3bw/frBRcA7wC8tLFjTnGgytMgiQmsEwimP5tHtbYHPo0RxKBtaAisJ6xVUZwIhyWBmHQizks5PJ/DoSqafVN0GjiGcdsLMfhKVYWZdCQnty2peZj6hJQXhd3GyhQnvMLPtE/b7nZkVmNnu0f7zzWxb4A7C0pqtq/px3P1Id++SIlFMAI43s63NbJso9up8/7uLWiALzex3CfXft4ZjVKcAqOpv+gPwZsJjHQgLEkkOUctCGoW7zzGzq4CXzayADbOCvm1m1wCTCItFvQc0i542GBhlZpMJ/2C/jcrLgOfMrByYRvSP292/NLO3ok7tF4G7E0K4B7jPzGYSWixnuvtqS2+9egMetnDZqxGmf+4TPXYicLqZrSV8Wz4lOkWzqX8RZg1d4GFW0AHAeDNbD0wlnGaCkFTGE1oMf3T378zsHuAed/+3hRlGXzOzCb5hOduk3P29qCU0jXDq7Y0ku44ABpvZnwj/5HsB90bvWfPo8empjpfgW6CThQWHVhCdtou+EPyE0PciOUSzzkpWMbMzCZ20F6TaN5dYWJzmEXc/ooZ9HgKed/cGGRuS5BjXEC4UuDlTx4iOs8rdW1ZTfjzQ1d3/N5PHl4an01AijcDDGsmDLYOD8nLEFmxYalhyiFoWIiKSkloWIiKSkpKFiIikpGQhIiIpKVmIiEhKShYiIpLS/wPhexyQn3xLxwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.plot(T33, T34, 'ok')\n", "plt.plot([950, 1000], [950, 1000], '-r')\n", "plt.xlabel('equation 33 (cpx-liq temp)')\n", "plt.ylabel('equation 34 (liq-only sat temp)')" ] }, { "cell_type": "code", "execution_count": null, "id": "2a4a964a-5f54-4c25-b244-18639fb3061b", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "6c8c3b21-3321-47e2-9b00-76d7b603cbbc", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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.8.5" } }, "nbformat": 4, "nbformat_minor": 5 }