{ "cells": [ { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "\n", "[*NBBinder test on a collection of notebooks about some thermodynamic properperties of water*](https://github.com/rmsrosa/nbbinder)" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "\n", "\"Open\"Open\"View\"View " ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "\n", "[<- Introduction](01.00-Introduction.ipynb) | [Water Contents](00.00-Water_Contents.ipynb) | [References](BA.00-References.ipynb) | [Low-Dimensional Fittings ->](03.00-Low_Dim_Fittings.ipynb)\n", "\n", "---\n" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "# Reading the Data\n", "\n", "A table with the variation of density and viscosity in terms of the temperature, at a fixed pressure of $1$ atmosphere, is available in [Batchelor (2000)](BA.00-References.ipynb). The data has been digitized and saved into a local `csv` file. Here we load the table from the file and view and plot the data." ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Importing the libraries\n", "\n", "First we import the libraries used in this particular notebook." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "slideshow": { "slide_type": "fragment" } }, "outputs": [], "source": [ "import pandas as pd\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Using pandas\n", "\n", "The data has been digitized to the local file `water.csv`. An easy way to retrieve it is with the [pandas.read_csv()](https://pandas.pydata.org/pandas-docs/stable/reference/api/pandas.read_csv.html) function of the `pandas` library:" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "slideshow": { "slide_type": "fragment" } }, "outputs": [], "source": [ "water_pd = pd.read_csv('water.csv', header=[0,1])" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "### Viewing the data with pandas\n", "\n", "The data is diplayed nicely with pandas:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "slideshow": { "slide_type": "fragment" } }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
tempdensityviscosity
Temperature (C)Density (g/cm^3)Viscosity (cm^2/s)
000.99990.01787
151.00001.51400
2100.99971.30400
3150.99911.13800
4200.99821.00400
5250.99710.89400
6300.99570.80200
7350.99410.72500
8400.99230.65900
9500.98810.55400
10600.98320.47500
11700.97780.41400
12800.97180.36600
13900.96530.32700
141000.95840.29500
\n", "
" ], "text/plain": [ " temp density viscosity\n", " Temperature (C) Density (g/cm^3) Viscosity (cm^2/s)\n", "0 0 0.9999 0.01787\n", "1 5 1.0000 1.51400\n", "2 10 0.9997 1.30400\n", "3 15 0.9991 1.13800\n", "4 20 0.9982 1.00400\n", "5 25 0.9971 0.89400\n", "6 30 0.9957 0.80200\n", "7 35 0.9941 0.72500\n", "8 40 0.9923 0.65900\n", "9 50 0.9881 0.55400\n", "10 60 0.9832 0.47500\n", "11 70 0.9778 0.41400\n", "12 80 0.9718 0.36600\n", "13 90 0.9653 0.32700\n", "14 100 0.9584 0.29500" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "water_pd" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "### Plotting the data\n", "\n", "We may also visualize both variations of density and viscosity using `matplotlib.pyplot`:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "slideshow": { "slide_type": "fragment" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAApIAAAFRCAYAAADHHE3hAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdf5wbVb3/8dfpbmHXAi1QKtuCpkLpTAW0UhCBK9iqLcRCBa8C8kuvcrnfC9Iq1VSkeotItFxbVBQREfAHqFh7WQJU5IcgirZSEWFSihiluysFSqu0u7Rb5vvHmZTZbJKdZLP5sft+Ph55bHLm5MzJJJl89vwa4/s+IiIiIiKlGlXrCoiIiIhIY1IgKSIiIiJlUSApIiIiImVRICkiIiIiZVEgKSIiIiJlUSApIiIiImVRICkNzRhznDHGN8bsX+u6RGGMOcsY01PrelSCMeYwY8zvjTE9xph0Cc+73hhz91DWLc8+HzHGXFXNfQ6GMeaA4HN99CDKuMAY80Il61UuY0zSGLOmRvveyxjzc2PMPxvpXCHSKBRIRhScgIrdbqx1HYeKMaYleI3vq3VdpK5cCbwITAWOq3FdBnIS8D/ZB8aYfxhjLqxhfarhJmBarSsR+CIwO/vAGHOrMea2Ku37fOBI4B1AG7CxSvutK7UM5kutgzHmQGPMLcaYp4wxrxpjrq3m/qU0zbWuQANpC91/H/CdnLTu6lZn8IwxowDj+/7OKu5ztO/7O6q1PxlSBwM3+b7/t1pXZCC+72+qdR2qzff9burkvOT7/svAyzXa/cHAE77vP1HtHQ/H812VfjdasQH/FcBFQ7gfqQTf93Ur8QZ8wB66vNveCPwU2IxtrbkdmBzangTWAB8H/o49uX4bG9TPBzqAF4AvY7+s2ef9A7gUuBXYCnQCn8jZ9z7Ad4HngX8C9wFvDW2/ICj7FOBJoBd7kj0G+GVQ3y3Ag8CROfv2Q7d0+LXk1OEC4IU8r/d84K/Aq8FrHRW8nr9if+z+BHwwwrGfCzwVPOd+4JygTvuH8rwT+HWQ51ng68Aeoe2PAFcD3wxe74vYE1b4eLcA/xu8H1uB3wEzQ9vnBPs9IXh924I8h+XU92NBHbYBPw/e456cPKcCa4Ee4BngC8DonOP/GeAG4F9Bebnv/d7AdcBzQTlPBOWODfb9vjzHsQfYu8BxbgKWABuAV4A/AieFjo2fc0sUKKc5ONbZ78PS4HXcHcpT9LMAOME+Tgne823An4ETQnl2D97PrqC+fweW5LznV4Xuh+veM4jjNNB3J3usPhK8/9uAv5DzWce2lv0x2Nca4OTgeUcX2O/FwefA5KSvAH5c4Ls4GbgDeAn7mX4SODW0/Q3Y88uLQT3/ABwX2n4R9vO5HfsdPDdn3xcBTwfH/3ngztzzQOh+7ufnaOA32fco9Lx9g/JOKnJO+FDwWrYDfwM+nfO+h/dzd4EysufG04D1wftwD/CGfK9hqM93wHRgJ7BP8Hgs9nz985zj/Xjo8VeDuncH+7kC2C1Uz9xjfnqwrazfjQL1LqsOxW7Y79e1A+Urd/+89h39OJDCfvbT2F6WWLD/rdjvw2FR6jHSbjWvQCPeKBBIAnsGH97rgMMAF7gZe3LdPciTxAYDP8Z2O8WDD/3d2IDSAT4YnETiobL/gf2h+jRwCHAhsIPXftybgN8DK4EZwBTgK8AmYL8gzwXYk/KvsT9cU4ExwHuBM4N9u0E9ngfGBs87MPiinQXsD4wPvZYoJ9aXgy/odOBwwGCDtCeCfU8Gzg2Ow7uLHPeDsD8W/xvU/UxsQL0rkASOCL70n8AGye/Anth/ECrnkeA9CJfzMvD/Qnl+BjwUnEwOAj4ZHDs32J4NJH+DDVxdbJDzWKiMd2J/SMLv2SZCgSQ2YNiMDYjfBLwbG2h8Mee9fz44tgcDlwT7flvovV+N/XF6T1BOHJgbbL8JWJFzLH8O3FLkWC8K6vWh4Bh9Gft5y77+/bGf9S8G98cUKGdx8JpPDX22/knfQLLoZ4HXAsknsF3UhwC3BMelJciT/ZE+DvvP3LHAOTnveTaQ3Afb2rEoqPvrB3GcBvruZH+k/o790ToY+2PXA7QFecYGx+iHwJuD17ie4oHkhOD9OCGUNi4oN/u+534X78F+Dw8LPiMnhY7xXkAGeCA4dgdjz3P/Fmw/A/v5vyA4/p/EnqPeE2w/Fvvd/FBw/N8KfCrnPJANJPfAnqfag+O/PzAaG2w/R99/oi7G/jPTVOA4HIP9jl0a1Os8bCDw8dB7/UNscLQ/hf8hyJ4bf48Nao8Afgv8Pt9rGOrzHTbw3ATMCx7PxX6uNgGjQueor4ee83ns+S4W5O8ELg22tWL/of5j6Ji3MIjfjQL1LrkOhb5boTJLCSTLOQbh7+i/B5+j27D/qN2DPZdODeqxOko9Rtqt5hVoxBuFA8n/B/w5J200Nmg5OXicDSTHhPLcEXzgm0Npu374gsf/ANpzyv4B8Mvg/knBl3+3nDxpgtYrXvuP7M0DvL7sSewDwePsFy23tSbqibWH4D/rIG0c9kfnyJznXkvOD3nO9q8S+g88SPsifQPJnwDX5OQ5OsizV+jY5ivn6eD+NOyP5Otz8twNfDW4nw0kjw9tnxWkZQPtFQXes3Ag+XtgYU6e04FNOe/993LyPAtcEtyfi20lOKjAcTsuON7ZH4YJweNiQfuLhFp2Qsft+tDjpynQEplTTjigaMIGLHdH/SzwWiB5bmj7QUHajODxdYRawPLUI9/36cLBHqcSvjufD+XZHRsEZvN8AhsktITyfIwigWSQ5y7gOznPeYEgEKP/d/Ep4DMFyroI+4/DuALb/wB8MyftVl47/5wZ7Pt1BZ7f51wRPPe2nDyvw/6zPC+U9kdC/1TlKfdnue97sK+nQ4+vp0BLZChP9tx4RChtSpB2XL7XUOAYV/J893/A1cH9/wWWBZ/b6djg9HngtCLPn0/o96hA/Sv2u1FuHSKUETmQLPMY5PuOzgjSwo0L2XP+HuXUZTjfNNmmso4AHGPMy9kbthtpDPaHL+sZ3/e3hh4/B3i+7/fmpE3IKf+3eR5nB9MfQdCykbP/g3P23e3njBUyxrQFM2nXG2O2YFuMxmG7uirhr37fMWqHYQPs+3Pq+pFsXY0x94W2/SF4nkv+YxB2BPCxnHLvDbYdVOR5vwXeZIxpCcoYBfwlp5xZOWWAbQXM6gz+Zt+3ovU1xhjsj8KSnP3cAOxtjNm7wH6y+8ruZzrwN9/3/0Ievu//Ghv0nR0knY3tsr83X35jzARsS87DOZt+TQmTN4wxrw/K2fWafTuuanUo24CfhZBix/q7wDuMMeuMMV8zxswJjm9kpR6n4DVG/e7sqrvv+69gf7zDn5O1vu+HZ/Pnfm7y+QHwAWPM7sHjD2O7tQuNyVsOfNEY87AxZokx5q2hbdOBP/i+v7nAcx2Kfx7uxLbyZowx3w9WJxgT4TXs4vv+NuBHwEcBjDFHYFv0vlfkaW6BemW/y6XYDjwaqs96bHBc6oSlks93BTyAHTpD8Pd+4FfB/UOx3f6/ymY2xpxhjPlNMJHsZWzQNNA5vOzfjXzKrEPFDHL/4fPLc8Hfx/Ok5f4uj3iabFNZo7Dj5M7Nsy28DEfuid4vkFZKoD8K2wU0K8+2LaH72/Js/xG2NeAT2Ob9bDfGbgPs81Xsf8Zho/Pk25rzOPu65mD/ww7bHvw9B/ufYjgtSmAwCrgGO14u17MRnp8tYwf2x9XP2Zb7WsLvWzZv9vUNVF8T5P0ctvUh1z8L7Ce7r6j7ARtofRTbqvsRbAtn7msL1yu7j1yFnlOsnGKifBayCh5r3/d/Z4yJBeXMwn6mHzHGxIu8znxKOU4Q/bsz2Pcvn59jW7VOMsasxg6l+GyhzL7vf9MYcwdBlzaQMMYs9n0/WawOoYC84OfB9/3Nxpi3YIOcd2OHNFxhjDnS9/1SZklfj33f9se+Dw8U+gcpW70i9Srlfc8+p5ihPN/l8wDwv8aYg7DB6IPAJOBEbA/EE77vvwBgjDke+D72XPJL7Dn/37HvQzGD+d3oYxB1qIgK7D/f+aXY+V0CCiQr61HsjO7nfN//1xCUn7um3NGAF9r354BXfN+PGjBlfySOxXYb3hWkHUDf/7p6sV+ippynP48dZxL2Vgb2eFDmgUErUD++72/Ik/wk9kcqLPeYPIrtgnl6gDrkO5Z/9X2/xxjzKPYHYrzv+1Fahgp5ssB+APB9/1VjzB+BQyLUt5hHgcuMMQcV+dG9CfiSMeYibCvOjYUK833/OWPMi9iu3t+ENh1H3//Qi/J9/x/GmJd4bSIFxpgmbLfRuiDbgJ+FEva3BTv2+MfGmB9if4gPxAZ4ubbT//MMJRyniN+dKJ4ETjXG7B60VkL/z00/vu9vM8asxLZEHoz9/Bb9vPq+/3ds8HmtMebz2EkhSexnaJ4xZlxuq6Tv+76x64Qehw2cs44L6p7NtwM7puweY8wXsMMa5mDHiefaju3iz63fH4wxjwf1OgMboBfzJP2XnjoO2+vzSp78xeyO/efxUYAggBvPa+fYITvfFfAYdrjBZdihOJuNMQ9gl93aif18Zx0H/CX4p4Cg/rGc8vJ95sv63Sig3DpUSq33P2IpkKysm4AFwMrgRLoB26x+KnZs3WCXSTneGHMJtvXqPdixdO8Ptt2JPSncboz5DHY81ERs60O77/u/y1dg8COxHjjHGLMW281xFaFlQ3zf7zXGbADebYz5HXaM32bsAPavGmMWYscDvgc7eaQo3/c3GWOuBq42xozGdk3thR043+37/g0FnvpN4CJjzFJsy8V0gm6wkC8BvzHGfA3bRbwVGxDM9n3/v0P5JueUczFBa47v+48bY34G/NAY8ynsOK3xwEzgSd/32wd6jYGvAffmvGfxnDz/A/zMGNOBHe/1Krb14a2+7xdsXcpxJ/ZHZ0VQ379gB4yP9n3/juA1vRAEHVdhx7XlC67CrgI+a4z5a1D2R7HdYB+JWKesq4FLjTHPYMddXYztkiOoV7mfhT6MMZ/Gztj9I/YYno4dVpLbApSVwX6fbsN+nl8M6hP5OEX57kR0M/ZzcL0x5krsZJVPR3zuD7CfrUOxk0oKMsZ8A7uKxHrsLP/38FogeDOwEHvuuhQ7+/0t2PF/D2Fn298U/ONzP3Zc7r9jA0WMMadizze/xh7392J7FLJBWK4McJYxZkqQf3NoaM93sN3w3djvRDFXAQ8bYz6LXS3jGGzwefEAz8tnO/BNY8zF2Jaor2PH0z0UbK/q+S74R/Mh7CTH5UGaZ+wFDeZiJzZlPYU9p30QO541jp2BHpYBDjLGHI4dGvJPyvzdKKCsOvi+369VNvgn7S3Bwz2AfYOhGD2+7xe6+EG5x0AGq5oDMofLjeLL/0zENq8/z2vLuXyHYBA7+Qf79hsMjp1FF55pnF3+5zZscPQPYEHOc8Ziu3U7sSfFv2N/XN4YbO8zMDz0vCOwM5t7sD8yHyJnIgX2C/kX7Ak2HUq/CNtl/DL2R+0S8iyHkWefo7AzP9NBXTcCq4B3DXDs5wV168F29ZxL/+V/jsa2jGTXrvsTcFloe3b5n29juz82BfUcFcqzO3YCzl+D+nVhuxLfEmzvN/Ca1yaFHBpK+0/sPxTbsLNUF9B/+Z+TsC123UF9fg/8Z857nzsxJHfyyL7YwPmFoJw/E5q0EOR5b1C/KMssZZf/6Qhe/2OEVhEI8kSZbDMa+Ebwul7CdhnnW/6n4GehwHHtMwEM+G9sEPlysK/7gKOKHK9/C47RK3nej1KOU9HvTm49C72n2NaUx4L6rMV+zotOtgm9T9mVC6bkbMudCPJt7He4JzjGPyQ0oQw70/W24PhtDV7XsTnf9fDyP+eFtr0LO15vE/az/ifgrELnAewavPcG71ef14kNsrrJmdxT5BiEl//5O/0niUWdbPMCNjh+Jngf7gViOfmqfb6bn/v5wU5UepVgUl+Qlp0Z/gJ2MudPsAF1eGLfGOzvymb6Lv9T1u9GnrqWXYc8ZeVbYswn9NtTqf2T5zsKHJDnc/nWIC3v0kcj+WaCAyR1zhjzD+zsxW/Uui6NzhjzCPBr3/cvqXVdqskYcy72RDvJL73bb8TQcaqtoDvyGews5z8Uz12xfV6APb+Or8b+RIYTdW2LDHPGmNdh1w1MAN9ScJSfjlNtBd2+E7Frlj5crSBSRAZHs49Ehr/F2O7SDdiuN8lPx6m2ZmGHkhyO7U4VkQagrm0RERERKYtaJEVERESkLAokRURERKQsw3KyzahRo/zW1tZaV0NERERkQNu2bfN932/Ixr1hGUi2traydWvuVapERERE6o8xptQLGdSNhox+RURERKT2FEiKiIiISFkUSIqIiIhIWRRIioiIiEhZFEiKiIiISFkUSIqIiIhIWRRIioiIiEhZqrKOZCyRugF4H7Axk4wfmme7Aa4GTgK2AedlkvFHg23nAp8Lsn4xk4zfVI0619q9X7+Z3W68ln22vsSmMXuz/bwLmHXRObWuloiIiMgu1WqRvBGYU2T7icCU4HY+8C2AWCK1D/B54O3AUcDnY4nU3kNa0whWru3g2OR9TE6kODZ5HyvXdlS0/Hu/fjP7XHsV47e+xChg/NaX2Ofaq7j36zdXdD8iIiIig1GVQDKTjD8IbCqS5RTg5kwy7meS8UeAcbFEqg2YDdyTScY3ZZLxl4B7KB6QDrmVaztYtOJxOjZ34wMdm7tZtOLxigaTu914LS07d/RJa9m5g91uvLZi+8ga6qBYREREhq96GSM5CXg29HhDkFYovR9jzPnGmDXGmDW9vb1DVtGlq9bRvWNnn7TuHTtZumpdxfaxz9aXSkovVzWCYhERERm+6iWQNHnS/CLp/RN9/zrf92f4vj+juXnohn52bs5/OcxC6eXYNCZ/732h9HJVIygWERGR4ateAskNwIGhxwcAnUXSa2biuNaS0sux/bwL6Gka3Setp2k028+7oGL7gOoExSIiIjJ81UsgeTtwTiyRMrFE6mhgSyYZ7wJWAe+NJVJ7B5Ns3huk1czC2VNpHd3UJ611dBMLZ0+t2D5mXXQOmy64hBfG7M2rwAtj9mbTBZdUfNZ2NYJi0DhMERGR4cr4ft6e4oqKJVK3ACcA44HnsDOxRwNkkvFrg+V/voGdSLMN+EgmGV8TPPejwGeDoq7IJOPfG2h/Y8aM8bdu3Vrpl7HLyrUdLF21js7N3Uwc18rC2VOZNz3v0M26lh0jGe7ebh3dxJWnHlax11ONfYiIiDQyY8w23/fH1Loe5ahKIFltQx1IDidDHRQfm7yPjjxd5ZPGtfJwYmbF9iMiItKoGjmQrMqC5FK/5k2fNKQtgxqHKSIiMnzVyxhJGaaqNQ5TREREqk+BpAypakxOAk3oERERqQV1bcuQynabD+U4zNwJPdmF1cP7FxERkcrTZBtpeJrQIyIijayRJ9uoa1sanib0iIiI1IYCSWl4mtAjIiJSGwokpeFVa0KPiIiI9KXJNtLwqjGhR0RERPrTZBuREgyXy2OKiEj9iDLZxnPcG4D3ARvdtHdokXxHAo8AH3LT3m2VrWl/6toWiSi7zFDH5m58XltmSGtWiohIFdwIzCmWwXPcJuDLwKpqVAgUSIpEtnTVul1rVWZ179jJ0lXralQjEREZKdy09yCwaYBsFwE/AzYOfY0sBZIiEWmZIRERqVee404C3g9cW839KpAUiUjLDImIyBBpNsasCd3OL6OM5cBn3LS3c8CcFaRAUiQiLTMkIiJDpNf3/Rmh23VllDEDuNVz3AzwAeCbnuPOq2gt89DyPyIRaZkhERGpV27am5y97znujcAdbtpbOdT7VSApUoJ50ycpcBQRkarzHPcW4ARgvOe4G4DPA6MB3LRX1XGRYVpHUqTOaK1KEZGRJco6kvVKLZIidSS7VmV2maHsWpWAgkkREak7mmwjUke0VqWIiDQSBZIidURrVYqISCNRIClSR7RWpYiINBIFkiJ1RGtViohII9FkG5E6orUqRUSkkWj5HxEREZEaauTlf9S1LSIiIiJlUdf2CLelvZ2Ny5bT29VFc1sbExbMZ+zcubWulgwxLXouIiKVoEByBNvS3k7XZYvxe3oA6O3spOuyxQAKJocxLXouIiKVoq7tEWzjsuW7gsgsv6eHjcuW16hGUg1a9FxERCpFgeQI1tvVVVK6DA9a9FxERCpFgeQI1tzWVlK6DA9a9FxERCpFgeQINmHBfExLS58009LChAXza1QjqQYtei4iIpWiyTYjWHZCjWZtjyxa9FxERCpFC5KLiIiI1JAWJBcRERGREUeBpIiIiIiURYGkiIiIiJRFk21EZMjoUowiIsObAkkRGRK6FKOIyPCnrm0RGRK6FKOIyPCnQFJEhoQuxSgiMvwpkBSRIaFLMYqIDH8aIykiQ2Lh7Kl9xkiCLsUoIlIuz3FvAN4HbHTT3qF5tn8Y+Ezw8GXgv9y099hQ10stkiIyJOZNn8SVpx7GpHGtGGDSuFauPPUwTbQRESnPjcCcItv/Chzvpr3DgcuB66pRKV0iUURERKSGol4i0XPcGHBHvhbJnHx7A392096Q/+euFkkRERGR4eU/gLuqsSONkRQRERGprWZjzJrQ4+t83y+ra9pz3HdhA8njKlKzASiQFBEREamtXt/3Zwy2EM9xDweuB050096Lg6/WwNS1LSIiItLgPMd9A7ACONtNe09Va7+abCMiIiJSQ1Em23iOewtwAjAeeA74PDAawE1713qOez1wGvC34Cm9btobdCvnQKoWSMYSqTnA1UATcH0mGU/mbH8jcAOwH7AJOCuTjG8Itn0ZiAdZL88k4z8uti8FkiIiItIoos7arkdV6dqOJVJNwDXAicA04IxYIjUtJ9tVwM2ZZPxwYAlwZfDcOPA24K3A24GFsURqr2rUW0REREQKq9YYyaOApzPJ+DOZZHw7cCtwSk6eacC9wf37Q9unAb/KJOO9mWR8K/AYxRfkFBEREZEqqFYgOQl4NvR4Q5AW9hi2bx/g/cCesURq3yD9xFgi9bpYIjUeeBdw4BDXVypoS3s762fOwnOnsX7mLLa0t9e6SiIiIlIB1QokTZ603MGZlwDHxxKptcDxQAfQm0nGfwHcCfwGuAX4LdDbbwfGnG+MWWOMWdPb22+z1MiW9na6LltMb2cn+D69nZ10XbZYwaSIiMgwUK11JDfQtxXxAKAznCGTjHcCpwLEEqk9gNMyyfiWYNsVwBXBth8B63N3ECzceR3YyTaVfwlSjo3LluP39PRJ83t62LhsOWPnzq1RrURERKQSqhVIrgamxBKpydiWxtOBM8MZgm7rTZlk/FVgEXYGd3aizrhMMv5iLJE6HDgc+EWV6i2D1NvVVVK6iIiINI6qdG1nkvFe4EJgFeABP8kk40/EEqklsUTq5CDbCcC6WCL1FPB6ghZI7BpJD8USqSexLY5nBeVJA2huayspXURERBqHFiSXIZUdIxnu3jYtLbRdvkRd2yIiIjT2OpK61rYMqWywuHHZcnq7umhua2PCgvkKIkVERIYBtUiKiIiI1FAjt0hWa/kfERERERlm1LUtIjKAlWs7WLpqHZ2bu5k4rpWFs6cyb3ruNRVEREYeBZIiIkWsXNvBohWP071jJwAdm7tZtOJxAAWTItKQPMd9GxAH3gKMAzZjryR4l5v21pRSlsZIiogUcWzyPjo2d/dLnzSulYcTM2tQIxEZbqo1RtJz3PcCXwL2BH6FXZLxX8FjF3tlwZeBz7ppb1WUMtUiKSJSRGeeILJYuohIHftP4L/ctLe6UAbPcY8EPoNd+3tAAwaSsURqNHA0/Zs/H8kk4zui7EREpFFNHNeat0Vy4rjWGtRGRKR8bto7LUKe1cAHopZZsGs7uGRhAjgX2ASk6dv8uTdwE5DMJOMvRN1hNahrW0QqJXeMJEDr6CauPPUwjZEUkYqoxfI/nuPuB3S7ae9lz3GbgHOAncAP3LT3atRyirVIPgR8F3hrJhnvyN0YS6QmAh8GHgSmlVJ5EZFGkQ0WNWtbRIaZO4ALgLXYy1LPBXYA04EFUQspFki+JZOMby+0MZOMdwJLY4nU1VF3JiLSiOZNn6TAUUSGm0OAPwb3zwKOwU60eYISAknN2hYRERGpoRp1bb8ATMIGlLe6ae/NnuOOAra4aW/PqOUUnWwTS6RGAZ8ADga+DTwHfAt4E/BL4NJirZYiIiIiUpfuAn4C7AvcGqRNA/oNZyxmoFnbS4G3Aq9iA8drg52Oxk7E2Rn8FREREZHG8THshOodwPeDtPHAF0opZKBA8kPAoUAT8DxwcyYZ/wtALJF6FDtQU4GkiIiISAPwHPcHwJ3A3W7auy68zU17D5Ra3qgBtu+ZScY3Z5LxF4F/ZYNIgEwy/iSwX6k7FBEREZGaSQEnAk96jvuw57iXeo47vdzCBmqRfCGWSO2VScb/iZ0ivksskdoP0IwWERERkSHmOe4NwPuAjW7aOzTPdgNcDZwEbAPOc9Peo7n53LR3C3BLkP8o7DW3v+M57v7A3djWynvctPevKPUaqEXya9hBmGSS8Vtytp0c7Eyk5ra0t7N+5iw8dxrrZ85iS3t7raskIiJSSTcCc4psPxGYEtzOx06OLshNe76b9n7npr3FbtqbARwBPAycAfzFc9z/jFIpLf8jDW9Leztdly3G7+nZlWZaWmi7fAlj586tYc1EREQGFnX5H89xY8AdBVokvw08ELQ44jnuOuAEN+11lVqf4Eo3+7hp7/mB8g7UIilS9zYuW94niATwe3rYuGx5jWokIiJSdZOAZ0OPNwRp/XiOe7TnuAs8x31vnm0JN+3tjBJEQoRAMpZIGV29RupZb1f+f7YKpYuIiNSZZmPMmtDt/DLKMHnS+nU7e457NnZo4gnAjZ7jpjzH3SOU5bOl7LRoIBlLpJqBHwH7lFKoSDU1t7WVlC4iIlJnen3fnxG6XTfwU/rZABwYenwA0Jkn3yJgjpv2TgEOAl4A7vccd1ywPV9AWlDBQDKWSO2BXfW8GTivlEJFqmnCgvmYlpY+aaalhQkL5teoRiIiIlV3O3CO57jGc9yjsZc6zNc1N8lNe78HcNNet5v2zgUeAGV6t7AAACAASURBVB70HPf15GnFLKbY8j/zgdcBczLJ+M5SChWppuyEmo3LltPb1UVzWxsTFszXRBsRERk2PMe9BdsdPd5z3A3A57FXGsRNe9diu6tPAp7GLv/zkQJFPec57hQ37a3PJrhpb6HnuNuAh7JlRlVw1nYskZoF/Aw4PZOM311KobWmWdsiIiLSKKLO2q4Ez3GXAi+7ae9/8mxbBFzhpr3Ik7GLLv8TS6T+DfgxcGYmGX+g9OrWhgJJERERaRTVDCQrrWjEmUnGH8Iufvnt6lRHRERERBrFgE2XmWT8T0C/dYZEREREpDF5jvvJSpQTqQ88k4z/rRI7ExEREZHa8hz3K9hLKg5asVnbfcQSqdcBBwPhRSvJJOO/qURFRERERGToBJc+vAGYTPHrdkcWKZCMJVLnAN8AtgPdoU0+8IZKVEREREREhtQC4DDgnW7a21aJAqO2SH4FOC2TjN9TiZ2KiIiISNX9FYgBU4C1lSgw6jpB27GrnouISANbubaDY5P3MTmR4tjkfaxc21HrKolIlbhp72fAfwF3eI775kqUGTWQvAz4aiyRGl+JnYqISPWtXNvBohWP07G5Gx/o2NzNohWPK5gUGUHctPdj4FxgZSXKK7ogeVYskXoHcCv2AuC7ngv4mWS8qRIVqSQtSC4i0t+xyfvo2NzdL33SuFYeTsysQY1EBGqzILnnuEe4ae8Pgy0naovk94GbgbcAhwS3KcFfkRFjS3s762fOwnOnsX7mLLa0t9e6SiKRdeYJIouli8jwVYkgEqJPttkXWJxJxgduvhQZpra0t9N12WL8nh4Aejs76bpsMQBj586tZdVEIpk4rjVvi+TEca01qI2I1APPcfciZ2lHN+11Rn1+1BbJ7wFnl1AvkWFn47Llu4LILL+nh43LlteoRiKlWTh7Kq2j+45Gah3dxMLZU2tUIxGpFc9x3+057jPAS8CG0O3ZUsqJ2iJ5FHBhLJG6FHguvCGTjL+zlB2KNKrerq6S0kXqzbzpkwBYumodnZu7mTiulYWzp+5KF5ER5bvA5dg5MGWPb4kaSH4nuImMWM1tbfR29m/tb25rq0FtRMozb/okBY4iAtACfM9NezsHU0ikWduNRrO2ZSjkjpEEMC0ttF2+RGMkRUSkbDWatZ3ArsCTdNNe2cFg1OV/vgbcGr6udiyROgb4YCYZn1/uzoeKAkkZKlva29m4bDm9XV00t7UxYcF8BZEiIjIoNQokpwCrgPHAC+Ftbtp7U9RyonZtnwFckpP2B+xilnUXSIoMlbFz5ypwFBGR4eA24CHgp1RhjKRP/xneTXnSRERERKT+TQamu2nv1cEUEjUQfAj4YiyRGgUQ/P1CkC4iIiIijeX/gEFf0ipqi+TFwB1AVyyR+hvwBqALUB+fiIiISOPZHbjdc9yHyFna0U1750QtJFIgmUnGN8QSqbcBb8deb/tZ4PeZZHxQzaEiIiIiUhNPBLdBKTprO5ZIPQvcBdwJ3JNJxhtiKrRmbYuIiEijqMWs7UoZaIzkUcDvsJdHzMQSqXtiidSCWCJ1yNBXTURERESGgue4Cc9xj8xJO8pz3E+XUk7kBcljiVQz8E7gpOC2G7al8k7g/kwy/kopOx5KapEUERGRRhGlRdJz3DnA1dhVc653014yZ/sbgJuAcUGehJv27ixSXhdwsJv2tobS9gCectPexKh1jzrZhkwy3gvcF9wuiSVSk7EB5UXAocBVxZ4fS6T6HIBMMp7M2f5G4AZgP2ATcFYmGd8QbPsKEMe2oN4DXJxJxoffJXlE0KLnIiLSl+e4TcA1wHuADcBqz3Fvd9Pek6FsnwN+4qa9b3mOOw3b0BcrUuxuwI6ctO3YSydGFjmQzJVJxv+KfVHXDJQ3lkj1OwCxROr2TDIePgBXATdnkvGbYonUTOBK4OzgCjrHAocH+X4NHA88UG7dRepV7mUYezs76bpsMYCCSRGRkeso4Gk37T0D4DnurcApQDiO8oG9gvtjgc4ByvwD8P+A5aG0C4BHS6lYpEAymHSTrwXwFWxguAL4VtBqmc9RwNOZZPyZoLx8B2AasCC4fz/2qjkE+23BRs4GGE3ONHWR4WLjsuV9ruUN4Pf0sHHZcgWSIiIj1yTsijlZG7Ar6YR9AfiF57gXAWOAdw9Q5gLgHs9xzwb+AhwMvB7b6BdZ1AXJvwa8BPwP8DFgCfAi8D3gx8AngC8VeX6+AzApJ89jwGnB/fcDe8YSqX0zyfhvsYFlV3BblUnGvYj1FmkovV1dJaWLiMiw0GyMWRO6nZ+z3eR5Tm4D3xnAjW7aOwA79PD7nuMWjPPctPcEcAiwFFgNfAWYmtNdPnDFI+Y7D3hPJhnf1UwaS6TuAn6RScbfHEuk7gd+CRSa6RPlAFwCfCOWSJ0HPAh0AL2xROpgwMWuXwlwTyyRemcmGX+wzw7sQT8fYLfddov4skTqS3NbG72d/XsjmtvaalAbERGpkl7f92cU2b4BODD0+AD6d13/BzAHwE17v/UctwUYD2wsVKib9l4Gbi2rxoGoLZJtwMs5aVuB7Kyep7CzhAoZ8ABkkvHOTDJ+aiYZnw5cGqRtwbZOPpJJxl/OJOMvY9e1PDp3B77vX+f7/gzf92c0N5c99FOkpiYsmI9p6TvO2bS0MGHB/BrVSERE6sBqYIrnuJM9x90NOB24PSfP34FZAJ7juthhgc+HM3iOuyJ3yZ9cnuMe6TnuiqgVixpxtQP/F0ukrsAGhQcAi4J0gHcAmSLPXw1MCWZ6d2APwJnhDLFEajywKbhaziLsDG6wB+bjsUTqSmzL5vH0HRgqMmxkx0Fq1raIiGS5aa/Xc9wLgVXY1W9ucNPeE57jLgHWuGnvduBTwHc8x12A7fU9z017ub2/1wLf9Bx3L+BXwDrgX8Ce2G7uE4DN2BngkURaRzKWSLVgB3H+O7YVshP4KbAkk4xviyVS+wO7ZZLxvxcp4yRsANgE3JBJxq+IJVJLgDWZZPz2WCL1AexMbR/btf3fmWT8lWDG9zexa1j6wN2ZZPyTxeqrdSRFRESkUVT7yjae484ATgQOw/YovwT8CbjTTXtrSykr8oLkjUSBpIiIiDSKYXmJxFgi9ZYoBUTNJyIiIiLDS7ExktfEEql/At8HfpUzY7sNO1bxHGAPbLeziIiIiIwgRbu2Y4nU+7CrnM8EdvLagEyDXe7nukwyXvA6jrWirm0RERFpFI3ctR11ss1oYAqvDchcX+QqNjWnQFJEREQaxbAPJBuNAkmR4ra0t2uJIRGROlGLQNJz3LXATcAtbtor+9LTCiRFRpgt7e10Xba4zzW9TUsLbZcvUTApDWXl2g6WrlpH5+ZuJo5rZeHsqcybnnv1XZH6V6NA8jTgw8Bs7LKL3wd+7qa97lLKUSApMsKsnzkr/2UYJ05kyn331qBGIqVbubaDRSsep3vHzl1praObuPLUwxRMSsOpZde257j7AB8EzgIOBVYAP3DT3n1Rnh/1EokiMkz0dnWVlC5Sj5auWtcniATo3rGTpavW1ahGIo3JTXubgJuxV735O3AacJ3nuE95jvvugZ4fKZCMJVInxxIpXcBaZBhobmsrKV2kHnVuzt/7VihdRPryHHeU57izPcf9AfaKhR8GksD+bto7GHu56h8MVE7UFsnLga5YIvWNWCL19nIrLSK1N2HBfExLS58009LChAXza1QjkdJNHNdaUrqI9NMJ/C/20ojT3LR3opv2fpQdI+mmvZ8B3kCFRB4jGVzB5izgDGArdlDmDzLJeKas6g8hjZEUKU6ztqXRaYykDCc1mmwzw017awZbTsmTbWKJlAFmYaPYQ4GHgW8Dt2SS8VcHW6FKUCApIjL8ada2DBc1CiQ3uWlvnzzpG920NyFqOSWNe4wlUgdhWyXPAl4FFmMHZl6IHZx5ainliYiIlGve9EkKHEXKNzo3wXPc0UBTKYVECiRjidR/A2cDBwM/Ac7OJOOPhLb/DNhYyo5FREREpLo8x30I8IEWz3EfzNl8APCbUsqL2iJ5IrYr+/8yyfj23I2ZZHxbLJFSa6SIiIhIfbseMMCRwHdD6T7wHBBp/cisqNfaviSTjF+VJ/2TmWT8q6XssBo0RlKkPmhSj4jIwGo0RtJx0156sOVEbZFcDPQLJIHPAXUXSIpI7eVeirG3s5OuyxYDKJgUEakBz3HPdtPe94OHx3iOe0y+fG7auyFqmUUDyVgiNTObL5ZIvQvbFJr1JuBfUXckIiPLxmXL+1zPG8Dv6WHjsuUKJEVEauMM7PKNYOe+5OMDlQkkea3vfPecQn3gH8BFUXckIiOLLsUoIlJf3LR3Uuj+uypRZtFAMpOMTwaIJVI3Z5LxcyqxQxEZGZrb2ujt7MybLiIiteU57n5At5v2XvYctwk4B+gFfuimvcjrgke6RKKCSBEplS7FKCJS1+4ApgT3vwRcAnwKu0pPZAVbJGOJlJdJxt3g/rPY7ux+Msn4G0rZoYiMDNlxkJq1LSJSlw4B/hjc/zBwDPAy8ASwIGohxbq2Px66f1aptRMRGTt3rgJHEZH6tBPYzXPcQ4Atbtr7u+e4o4A9SimkYCCZScZ/Hbr/q7KrKSIiIiL15i7s1Qr3BW4N0qYBHaUUEvUSiZ8E7ssk43+MJVJHBzvuBT6cScZ/W8oORURERKQ0nuPOAa7GXgv7ejftJfPk+SDwBexwxMfctHdmkSI/BpwL7OC1JYHGB8+PLOqC5At4bSmgK7GLkP8LWA68vZQdioiIiEh0wazqa4D3ABuA1Z7j3u6mvSdDeaYAi4Bj3bT3kue4E4qV6aa9V4Drgu7s13uO+5yb9h4otW6RZm0DYzPJ+JZYIrUn8Bbg65lk/LvA1FJ3KCIiIiIlOQp42k17z7hpbzu2K/qUnDwfB65x095LAG7a21isQM9x9/Ic9yagG9ud3e057k2e444tpWJRA8lnY4nUMcDpwIOZZHxnLJHaCztQU0RERETK12yMWRO6nZ+zfRLwbOjxhiAt7BDgEM9xH/Yc95GgK7yYr2En1hwGtAZ/XxekR694xHwLgduA7cBpQdr7gN+XsjMRERER6afX9/0ZRbabPGm5yzI2Y9eFPAE4AHjIc9xD3bS3uUCZc4A3uWlvW/D4Kc9xPwL8JXq1IwaSmWT8TmBiTvJPg5uIiIiIDJ0NwIGhxwcAuZcO2wA84qa9HcBfPcddhw0sVxcoswfYD/hbKG088EopFYvaIkkskRqLHROZu77QfaXsUERERERKshqY4jnuZOx4xtOB3BnZK4EzgBs9xx2P7ep+pkiZ1wP3eI77VWww+Ubs5OrrSqlY1OV/zsPOFnoZ2Bba5ANvKmWHIiIiIhKdm/Z6Pce9EFiFXf7nBjftPeE57hJgjZv2bg+2vddz3Cexc1gWumnvxSLFXoFt1TwT2+vcCXwFuKGUuhnfz3vlwz5iiVQH8LFMMn5XKYXXypgxY/ytW7fWuhoiMkxsaW/XpR5FZMgYY7b5vj+m1vUoR9RZ283AL4ayIiIi9WhLeztdly2mt7MTfJ/ezk66LlvMlvb2WldNRKRsnuN+zXPcY3LSjvEcd3kp5UQNJL8MfC6WSEXNLyIyLGxcthy/p6dPmt/Tw8ZlJZ1rRUTqzRnAmpy0P9B/7GVRpVzZZn/g07FEqk9/eyYZf0MpOxQRaSS9XV0lpYuINAif/g2KTXnSiooaSJ5VSqEiIsNFc1ub7dbOky4i0sAeAr7oOe6n3bT3anCpxC8E6ZFFmmzTaDTZRkQqJTtGMty9bVpaaLt8iSbciEhF1GKyjee4BwB3AG3Y5X/eAHQBc920tyFqOVFnbe8OLMb2p++bScbHxhKp9wKHZJLxb5RR/yGlQFJEKkmztkVkKNVq1nbQCnkUdrHzZ4Hfu2nv1VLKiNq1vQx7TccPA9klgJ4I0usukBQRqaSxc+cqcBSRYScIGh8BHvEc913AccCDpZQRdUDl+4EzM8n4b4FXATLJeAf9LxguIiIiInXOc9xfeY57bHD/M8CtwC2e4362lHKitkhuz80bS6T2A4qtmC4iIiIRrFzbwdJV6+jc3M3Eca0snD2VedPVViND6lBsayTAx4ETsFcwfBj4UtRCorZI/hS4KZZITQaIJVJt2C7tW6PuSERERPpbubaDRSsep2NzNz7QsbmbRSseZ+XajlpXTYa3UYDvOe5BgHHTnuemvWeBvUstJIrPAhngcWAcsB57Tcb/KWVnIiIi0tfSVevo3rGzT1r3jp0sXbWuRjWSEeLX2EbBq4CfAwRB5QulFBKpazuTjG8H5gPzgy7tFzLJ+PBbN0hERKTKOjd3l5QuUiHnAZ8CngeWBmkOcHUphUQKJGOJ1DTg34B9gE3YxSqfLGVHIiJSe1rKqP5MHNdKR56gceK41hrURkYKN+29iO1xDqelSi2n6DqSsUTKAN8FzgU2YLuzJwETge8DH63HlkmtIyki0p8WV69P2TGS4e7t1tFNXHnqYZpwM0JUax1Jz3EvddPeFcH9JYXyuWlvcdQyB2qRPB87i+foTDK+OpsYS6SOBG4B/hO4NurORESkdjYuW94niATwe3rYuGy5AskaygaLmrUtVXBA6P6BlShwoBbJXwPJTDJ+R55t7wMWZZLxYytRkUpSi6SISH+eOw3ynfONwfU0WkmkVqp5ZRvPcUeVevWaYgaatT0N+FWBbb8KtouISANobmsrKV1EhqUOz3G/4jnuoZUobKBAsimTjP8r34YgPeryQSIiUmMTFszHtLT0STMtLUxYML9GNRKRGrgAmAys9hz3Uc9xL/Ycd79yCxuoa3sbEAdMgSztmWQ8UlNsLJGag51S3gRcn0nGkznb3wjcAOyHnRl+ViYZ3xBLpN6FvaZ3lgOcnknGVxbal7q2RUTy06xtkfpTza7tLM9xxwEfAs4GZgC/AG4CbnfT3o6o5QwUSGaAorOyM8n45IF2EkukmoCngPdgZ3+vBs7IJONPhvL8FLgjk4zfFEukZgIfySTjZ+eUsw/wNHBAJhnfVmh/CiRFRESkUdQikAzzHHcyNqD8GPA6N+2Nj/rcorO2M8l4bHBV2+Uo4OlMMv4MQCyRuhU4hb5rUU4DFgT37wfytTh+ALirWBApIiIiItF4jrs7cCTwduD1wG9KeX61xjhOAp4NPd4QpIU9BpwW3H8/sGcskdo3J8/p2GWHRERERKRMnuMe5znudcBzwBeBR4BD3LT3rlLKiXRlmwrIN8Yyt8v8EuAbsUTqPOBBoAPozW6MJVJtwGHAqrw7MOZ87LqX7LbbboOvsYiI1DWN9xQpnee4X8B2Y+8D/BSIu2nv4XLLq1YguYG+C18egL1Kzi6ZZLwTOBUglkjtAZyWSca3hLJ8EPh5JhnPOwDU9/3rgOvAjpGsXNVFRKTe5F6lp7ezk67L7MU4FEyKFHU0cCmw0k17PQNlHki1AsnVwJRYIjUZ29J4OnBmOEMskRoPbMok468Ci7AzuMPOCNJFRGSE01V6ZKTxHLfP6jdu2ksWyPcBbEvjkW7aW5O73U17cypZr6qMkcwk473AhdhuaQ/4SSYZfyKWSC2JJVInB9lOANbFEqmnsIM9r8g+P5ZIxbAtmoUWRxcRkRGkt6urpHSRRuY5bhNwDXAidnLyGZ7j9rsojOe4ewKfAH5XrboVXf6nUWn5HxGR4W39zFn0dnb2S2+eOJEp991bgxqJlG+g5X88x30H8AU37c0OHi8CcNPelTn5lgO/xM47uSRfi2Sl6co0IiLScHSVHhlmmo0xa0K383O2D7j6jee404ED3bR3xxDXtY9qjZEUERGpmOw4SM3almGi1/f9GUW2F139xnPcUdirAJ5X4XoNSIGkiIg0pLFz5ypwlJFioNVv9gQOBR7wHBdgf+B2z3FPHurubQWSIiIiNaY1MWUAq4EpwaUM+61+46a9LcCuyxp6jvsAGiMpIiIy/GXXxOzt7ATf37Um5pb29lpXTeqEm/b6rX7jpr0nPMdd4jnuycWfPbQ0a1tERKSGNANdBpq1Xc/UIikiIlJDWhNTGpkCSRERkRpqbmsrKV2kniiQFBERqSGtiSmNTIGkiIhIDY2dO5e2y5fQPHEiGEPzxIm0Xb6kIWdtb2lvZ/3MWXjuNNbPnKUJQyOAJtuIiIjIoGVnn/s9PbvSTEtLwwbF1aTJNiIiIjKibVy2vE8QCeD39LBx2fIa1UiqQYGkiIiIDJpmn49MCiRFRERk0DT7fGRSICkiIiKDptnnI5MCSRERERm04TT7HDQDPSrN2hYREREJqfYMdM3aFhERERkmNAM9OgWSIiIiIiGagR5dc60rICIiIsPDyrUdLF21js7N3Uwc18rC2VOZN31SratVsua2Nno7O/OmS19qkRQREZFBW7m2g0UrHqdjczc+0LG5m0UrHmfl2o5aV61kmoEenQJJERERGbSlq9bRvWNnn7TuHTtZumpdjWpUvuE2A30oqWtbREREBq1zc3dJ6fVu7Ny5ChwjUIukiIiIDNrEca0lpcvwoEBSREREBm3h7Km0jm7qk9Y6uomFs6fWqEZSDeraFhERkUHLzs4eDrO2JTpd2UZERESkhnRlGxEREREZcRRIioiIiEhZFEiKiIiISFk02UZERESkznmOOwe4GmgCrnfTXjJn+yeBjwG9wPPAR92097ehrpdaJEVERETqmOe4TcA1wInANOAMz3Gn5WRbC8xw097hwG3AV6pRN7VIioiIiNS3o4Cn3bT3DIDnuLcCpwBPZjO4ae/+UP5HgLOqUTG1SIqIiIjUVrMxZk3odn7O9knAs6HHG4K0Qv4DuKvSlcxHLZIiIiIitdXr+/6MIttNnrS8C4F7jnsWMAM4vhIVG4gCSREREZH6tgE4MPT4AKAzN5PnuO8GLgWOd9PeK9WomAJJERERkfq2GpjiOe5koAM4HTgznMFz3OnAt4E5btrbWK2KaYykiIiISB1z014vcCGwCvCAn7hp7wnPcZd4jntykG0psAfwU89x/+g57u3VqJuutS0iIiJSQ7rWtoiIiIiMOAokRURERKQsCiRFREREpCwKJEVERESkLAokRURERKQsCiRFREREpCwKJEVERESkLAokRURERKQsCiRFREREpCwKJEVERESkLAokRURERKQszdXaUSyRmgNcDTQB12eS8WTO9jcCNwD7AZuAszLJ+IZg2xuA64EDAR84KZOMZ6pVdxERERHpryotkrFEqgm4BjgRmAacEUukpuVkuwq4OZOMHw4sAa4MbbsZWJpJxl3gKGDj0NdaRERERIqpVovkUcDTmWT8GYBYInUrcArwZCjPNGBBcP9+YGWQdxrQnEnG7wHIJOMvV6nOIiIiIlJEtQLJScCzoccbgLfn5HkMOA3b/f1+YM9YIrUvcAiwOZZIrQAmA78EEplkfOeQ11pERERECqrWZBuTJ83PeXwJcHwskVoLHA90AL3YYPffgu1HAm8Czuu3A2PON8asMcas6e3trWDVRURERCSfagWSG7ATZbIOADrDGTLJeGcmGT81k4xPBy4N0rYEz12bScafySTjvdgu77fl7sD3/et835/h+/6M5uaqzSESERERGbGqFUiuBqbEEqnJsURqN+B04PZwhlgiNT6WSGXrswg7gzv73L1jidR+weOZ9B1bKSIiIiI1UJVAMmhJvBBYBXjATzLJ+BOxRGpJLJE6Och2ArAulkg9BbweuCJ47k5st/a9sUTqcWw3+XeqUW8RERERKcz4fu5QxcY3ZswYf+vWrbWuhoiIiMiAjDHbfN8fU+t6lENXthERERGRsmhWShm2tLezcdlyeru6aG5rY8KC+YydO7fW1RIREZEKWbm2g6Wr1tG5uZuJ41pZOHsq86ZPqnW16o4CyRJtaW+n67LF+D09APR2dtJ12WIABZMiIiLDwMq1HSxa8TjdO+yS1R2bu1m04nEABZM51LVdoo3Llu8KIrP8nh42LlteoxqJiIhIJS1dtW5XEJnVvWMnS1etq1GN6pdaJEvU29VVUrqIiIg0ls7N3SWlV4PnuHOwV/9rAq53014yZ/vuwM3AEcCLwIfctJcZ6nqpRbJEzW1tJaWLiIhIY5k4rrWk9KHmOW4TcA1wIjANOMNz3Gk52f4DeMlNewcDy4AvV6NuCiRLNGHBfExLS58009LChAXza1QjERERqaSFs6fSOrqpT1rr6CYWzp5aoxpxFPC0m/aecdPeduBW4JScPKcANwX3bwNmeY6b7xLVFaVAskRj586l7fIlNE+cCMbQPHEibZcv0UQbERGRYWLe9ElceephTBrXigEmjWvlylMPG8qJNs3GmDWh2/k52ycBz4YebwjS8uZx014vsAXYd6gqnKUxkmUYO3euAkcREZFhbN70SdWcod3r+/6MItvztSzmXlEmSp6KU4ukiIiISH3bABwYenwA0Fkoj+e4zcBYYNNQV0wtkiIiIiL1bTUwxXPcyUAHcDpwZk6e24Fzgd8CHwDuc9OeWiRFRERERrJgzOOFwCrAA37ipr0nPMdd4jnuyUG27wL7eo77NPBJIFGNuhnfH/JgterGjBnjb926tdbVEBERERmQMWab7/tjal2PcqhFUkRERETKokBSRERERMqiQFJEREREyqJAUkRERETKMiwn2xhjXgWqcWX1ZqC3CvuR6PSe1Ce9L/VH70l90vtSf6rxnrT6vt+QjXvDMpCsFmPMmgFWopcq03tSn/S+1B+9J/VJ70v90XtSXENGvyIiIiJSewokRURERKQsCiQH57paV0D60XtSn/S+1B+9J/VJ70v90XtShMZIioiIiEhZ1CIpIiIiImVRIFkGY8wcY8w6Y8zTxpiqXBRd+jPGHGiMud8Y4xljnjDGXByk72OMuccYsz74u3et6zrSGGOajDFrjTF3BI8nG2N+F7wnPzbG7FbrOo40xphxxpjbjDHp4DvzDn1XassYsyA4d/3ZGHOLMaZF35Xqh5LehgAACBZJREFUM8bcYIzZaIz5cygt73fDWF8Lfv//ZIx5W+1qXh8USJbIGNMEXAOcCEwDzjDGTKttrUasXuBTvu+7wNHAfwfvRQK41/f9KcC9wWOprosBL/T4y8Cy4D15CfiPmtRqZLsauNv3fQd4C/b90XelRowxk4BPADN83z8UaAJOR9+VWrgRmJOTVui7cSIwJbidD3yrSnWsWwokS3cU8LTv+8/4vr8duBU4pcZ1GpF83+/yff/R4P6/sD+Mk7Dvx01BtpuAebWp4chkjDkAiAPXB48NMBO4Lcii96TKjDF7Ae8Evgvg+/523/c3o+9KrTUDrcaYZuB1QBf6rlSd7/sPAptykgt9N04BbvatR4Bxxpi26tS0PimQLN0k4NnQ4w1BmtSQMSYGTAd+B7ze9/0usMEmMKF2NRuRlgOfBl4NHu8LbPZ9P3tlCH1nqu9NwPP/v727j5GrqsM4/n1aQIRuraVBbdqiIMIfYqQYxNpAA4gGwZcIAS1KqiQ00SAhNUGNkpI0aRCqIr40voVqNSggNpGQamrVhCBS61ugiViw3Za+020ptRb6+Mc9U67TfZmd7e4s7PNJJrv33nPuPXcmZ/LLOed3B/hhWXLwPUknkr7SMbY3AbcDG6gCyB5gDekro0VffSMxQJMEkoOnXvYl9b2DJE0A7gNutL2n0+0ZyyRdBmyzvaa+u5ei6TMj6xhgJvBt22cD+8g0dkeVNXcfBN4ETAVOpJo2bZa+Mrrk+6xJAsnB6wam17anAZs71JYxT9KxVEHkctv3l91bG1MN5e+2TrVvDHo38AFJT1Mt+7iQaoRyUpm+g/SZTugGum3/sWzfSxVYpq90zsXAU7a32z4I3A/MIn1ltOirbyQGaJJAcvD+BJxeMuuOo1ocvaLDbRqTytq77wNP2F5SO7QCuLb8fy3wy5Fu21hl+/O2p9l+I1XfWGV7LvBb4IpSLJ/JCLO9Bdgo6Yyy6yLgcdJXOmkDcJ6kE8p3WeMzSV8ZHfrqGyuAT5Ts7fOAnsYU+FiVB5K3QdKlVKMs44Ef2F7U4SaNSZJmA38A/s5L6/G+QLVO8mfADKov6yttNy+kjmEmaQ6wwPZlkk6lGqGcDKwFrrF9oJPtG2skvZ0qAeo4YD0wj2owIX2lQyQtBK6iegLFWuA6qvV26SsjSNJPgTnAFGArcAvwAL30jRL030WV5f08MM/2Y51o92iRQDIiIiIi2pKp7YiIiIhoSwLJiIiIiGhLAsmIiIiIaEsCyYiIiIhoSwLJiIiIiGhLAsmIiJeZ8uzBdZKmtFj+XEmrh7lZETEGJZCMiGEn6bna65Ck/bXtuZ1u31BI2lKeaTqSPg08ZHtHrR2zJK2U1CNpp6RHGu+t7UeBQ5LeM8LtjIhXuASSETHsbE9ovKge7nt5bd/yTrevL7Wfqhtt17ge+FHtHHOAlcBDVL/dPAW4AXh/rc7yUi8i4qhJIBkRHSdpvKQvSVovaYek5ZImlWNnSnpB0qckbSqjbZ+U9C5J/5C0W9KS2rnmS1olaamkPZIel3R+7fhkScvKSOJGSbdIGtdU95uSngVuLtdfLWmXpO2S7pbUVcr/HDgZWFlGV2+Q9D5JTzbd3+FRS0mLJf1E0j2S9gJX93f/vbxXbynX/HNt9+3AUttLbO9y5VHbH6uVWQ28V9L4Nj+miIgjJJCMiNHgc8AlwGxgGnAQ+Grt+HjgbcCpVD/t9w1gAXBB2T9P0jtr5c8H/gqcBCwGHpA0sRxbDvSUc50LfAj4eFPdv1CN6t1R9t0KvB44CzgD+CKA7SuBbcAlZXT1zhbv9yPA3cBrgPtauP+6s4B/uvwsWQk4zwHu7e+Ctv8FvAo4rcU2RkQMKIFkRIwG1wM3295s+z/AQuCq8ru2DbfaPmB7RdleZnun7Q3Aw8DZtbIbbX/L9kHby4BuqtG4U6gCxZtsP2/7GeBO4Opa3fW2v2v7Rdv7ba+zvcr2f21vAb5GFcAOxe9sP2j7kO39Ld5/wyRgb237pPL3mRauu7fUj4g4KoZ9/U9ERH9KsDQdeFCSa4fG8VKQ9KLtnbVj+4GtTdsTatvdTZf5NzAVOAU4Hthei9HGAfWp6I1N7ZsKfB2YBXSV8q0Ebf05fI0W7n9HU91nSzsaGu/LG4CnB7huF7C7jfZGRPQqI5IR0VFlinYTcKHtSbXX8fWs5EGa1rQ9A9hMFcA9B7y2dp2JtmfWm9RU9yvAPuCtticC1wHqp/w+4ITGhqRjgclNZQ7XaeP+/wa8uTFaaXs3sIZqurxPkk4DDgDr+ysXETEYCSQjYjT4DrBY0nQASSdLunwI55teEmeOkXQNVSC50vZTwCPAbZK6JI2TdPoAj+/pogo+90iaAdzUdHwr1XrLhieAyZIuKkHkQgb+rm35/m0/Wa5Zn8pfAMyXdGNJJpKkcyT9uFbmAuDXtl8YoC0RES1LIBkRo8FtwG+AVSWT+WFgZv9V+vV7qkBrF1VizIdt95RjH6VaJ7iuHL8HeF0/5/oyVRJMD/ALquSYukXAopI9/pkyivhZqqSebmALR05PNxvs/S+lliBkezVVss6lVNPbO4C7gF/V6sylClgjIo4alcS/iIhXBEnzgStsX9zptgwXSa8G1gKzW5n+l/QO4A7bQ00Sioj4P0m2iYh4mSmZ3mcOovxjDD3TPCLiCJnajoiIiIi2ZGo7IiIiItqSEcmIiIiIaEsCyYiIiIhoSwLJiIiIiGhLAsmIiIiIaEsCyYiIiIhoSwLJiIiIiGjL/wB95Vq/znkP4gAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax1 = plt.subplots(figsize=(10,5))\n", "\n", "color = 'tab:blue'\n", "ax1.set_xlabel(water_pd['temp'].columns[0], fontsize=12)\n", "ax1.set_ylabel(water_pd['density'].columns[0], color=color, fontsize=12)\n", "ax1.plot(water_pd['temp'], water_pd['density'], 'o', color=color)\n", "ax1.tick_params(axis='y', labelcolor=color)\n", "\n", "ax2 = ax1.twinx() \n", "\n", "color = 'tab:red'\n", "ax2.set_ylabel(water_pd['viscosity'].columns[0], color=color, fontsize=12)\n", "ax2.plot(water_pd['temp'], water_pd['viscosity'], 'o', color=color)\n", "ax2.tick_params(axis='y', labelcolor=color)\n", "\n", "plt.title('Temperature-dependency of density and viscosity of pure water at 1 atm', \n", " fontsize=14)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "\n", "\n", "---\n", "[<- Introduction](01.00-Introduction.ipynb) | [Water Contents](00.00-Water_Contents.ipynb) | [References](BA.00-References.ipynb) | [Low-Dimensional Fittings ->](03.00-Low_Dim_Fittings.ipynb)" ] } ], "metadata": { "file_extension": ".py", "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.7.3" }, "mimetype": "text/x-python", "name": "python", "npconvert_exporter": "python", "pygments_lexer": "ipython3", "version": 3 }, "nbformat": 4, "nbformat_minor": 4 }