From 6df9bfcda6061160b054418831a7caf71823ed3c Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Georg=20=C2=B4Brantegger?= Date: Wed, 15 Jun 2022 14:59:48 +0200 Subject: [PATCH] commit for sync 15.06. --- Main_Program.ipynb | 30 +++++++++++++ e-RK4-Test.ipynb | 109 +++++++++++++++++++++++++++++---------------- 2 files changed, 100 insertions(+), 39 deletions(-) diff --git a/Main_Program.ipynb b/Main_Program.ipynb index e69de29..3557ac5 100644 --- a/Main_Program.ipynb +++ b/Main_Program.ipynb @@ -0,0 +1,30 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3.8.13 ('DT_Slot_3')", + "language": "python", + "name": "python3" + }, + "language_info": { + "name": "python", + "version": "3.8.13" + }, + "orig_nbformat": 4, + "vscode": { + "interpreter": { + "hash": "4a28055eb8a3160fa4c7e4fca69770c4e0a1add985300856aa3fcf4ce32a2c48" + } + } + }, + "nbformat": 4, + "nbformat_minor": 2 +} diff --git a/e-RK4-Test.ipynb b/e-RK4-Test.ipynb index 147f48f..9acfd87 100644 --- a/e-RK4-Test.ipynb +++ b/e-RK4-Test.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "code", - "execution_count": 6, + "execution_count": 8, "metadata": {}, "outputs": [], "source": [ @@ -13,63 +13,92 @@ }, { "cell_type": "code", - "execution_count": 7, + "execution_count": 9, + "metadata": {}, + "outputs": [], + "source": [ + "# define constants\n", + "initial_level = 5.\n", + "initial_influx = 3. \n", + "initial_outflux = 0. \n", + "initial_pipeline_pressure = 0. \n", + "\n", + "total_min_level = 0.01\n", + "total_max_time = 100\n", + "\n" + ] + }, + { + "cell_type": "code", + "execution_count": 10, "metadata": {}, "outputs": [ { - "ename": "KeyboardInterrupt", - "evalue": "", - "output_type": "error", - "traceback": [ - "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m", - "\u001b[1;31mKeyboardInterrupt\u001b[0m Traceback (most recent call last)", - "\u001b[1;32mv:\\georg\\Documents\\Persönliche Dokumente\\Arbeit\\Kelag\\Coding\\Python\\DT_Slot_3\\Kelag_DT_Slot_3\\e-RK4-Test.ipynb Cell 2'\u001b[0m in \u001b[0;36m\u001b[1;34m()\u001b[0m\n\u001b[0;32m 12\u001b[0m \u001b[39mwhile\u001b[39;00m V\u001b[39m.\u001b[39mlevel \u001b[39m>\u001b[39m \u001b[39m0.005\u001b[39m:\n\u001b[0;32m 13\u001b[0m t \u001b[39m=\u001b[39m t\u001b[39m+\u001b[39mV\u001b[39m.\u001b[39mtimestep\n\u001b[1;32m---> 14\u001b[0m V\u001b[39m.\u001b[39;49me_RK_4()\n\u001b[0;32m 15\u001b[0m V\u001b[39m.\u001b[39mlevel \u001b[39m=\u001b[39m V\u001b[39m.\u001b[39mupdate_level(V\u001b[39m.\u001b[39mtimestep)\n\u001b[0;32m 16\u001b[0m V\u001b[39m.\u001b[39mset_volume()\n", - "File \u001b[1;32mv:\\georg\\Documents\\Persönliche Dokumente\\Arbeit\\Kelag\\Coding\\Python\\DT_Slot_3\\Kelag_DT_Slot_3\\Ausgleichsbecken_class_file.py:81\u001b[0m, in \u001b[0;36mAusgleichsbecken_class.e_RK_4\u001b[1;34m(self)\u001b[0m\n\u001b[0;32m 79\u001b[0m Y3 \u001b[39m=\u001b[39m yn \u001b[39m+\u001b[39m dt\u001b[39m/\u001b[39m\u001b[39m2\u001b[39m\u001b[39m*\u001b[39mFODE_function(Y2, h_hs, alpha, p_hs)\n\u001b[0;32m 80\u001b[0m Y4 \u001b[39m=\u001b[39m yn \u001b[39m+\u001b[39m dt\u001b[39m*\u001b[39mFODE_function(Y3, h_hs, alpha, p_hs)\n\u001b[1;32m---> 81\u001b[0m ynp1 \u001b[39m=\u001b[39m yn \u001b[39m+\u001b[39m dt\u001b[39m/\u001b[39m\u001b[39m6\u001b[39m\u001b[39m*\u001b[39m(FODE_function(Y1, h, alpha, p)\u001b[39m+\u001b[39m\u001b[39m2\u001b[39;49m\u001b[39m*\u001b[39;49mFODE_function(Y2, h_hs, alpha, p_hs)\u001b[39m+\u001b[39m \\\n\u001b[0;32m 82\u001b[0m \u001b[39m2\u001b[39m\u001b[39m*\u001b[39mFODE_function(Y3, h_hs, alpha, p_hs)\u001b[39m+\u001b[39m FODE_function(Y4, h, alpha, p))\n\u001b[0;32m 84\u001b[0m \u001b[39mself\u001b[39m\u001b[39m.\u001b[39moutflux \u001b[39m=\u001b[39m ynp1\u001b[39m*\u001b[39m\u001b[39mself\u001b[39m\u001b[39m.\u001b[39marea_outflux\n", - "\u001b[1;31mKeyboardInterrupt\u001b[0m: " - ] + "data": { + "text/plain": [ + "[]" + ] + }, + "execution_count": 10, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAD4CAYAAADM6gxlAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAS0UlEQVR4nO3da4zldX3H8ff33PbM7IXddQddubiQIpXaVHBqoRhq0VKqxj4xjSZatG1I00vUNrEQ+8SHbRpjjYmFeKmpeCtqa0kUDWKUBwWHggouCALKCrgD7MLuzO7O5fz64PzPzJnZYefssmfO/3fO+5Vs9pz//Gfm+9tdPnzn+79FSglJUj4qgy5AknRyDG5JyozBLUmZMbglKTMGtyRlptaPL7pr1660Z8+efnxpSRpKd99999MppYle9u1LcO/Zs4epqal+fGlJGkoR8fNe93VUIkmZMbglKTMGtyRlxuCWpMz0FNwRsT0ibo6IByJib0Rc1u/CJElr6/Wskn8FvplSentENIDxPtYkSTqBdYM7IrYBVwDvAUgpzQFz/S1LkvRCehmVnA9MA5+JiHsi4pMRsXn1ThFxbURMRcTU9PT0iyrq9gf288PHD76oryFJw6qX4K4BlwCfSCldDMwA163eKaV0Y0ppMqU0OTHR08U/a5qdW+C9//4DrvnMXaf8NSRpmPUS3PuAfSmlO4v3N9MO8r546FeHATg4O8/R+cV+fRtJyta6wZ1Segp4PCIuLDa9EfhJvwp67JmZpddPHDzSr28jSdnq9aySvwVuKs4oeQR4b78KenZm+bjnvgNHOH9iS7++lSRlqafgTindC0z2t5S2A7PzS6+7Q1yS1Fa6KycPzi6HtcEtSccrXXAfmJ3nnJ1jVAIOzBrckrRaX+7H/WIcnJ1j5+ZNzB5btOOWpDWULrifPzLPGeMNZjY37LglaQ2lG5UcmV9kvF7ljLE6B7sOVEqS2koX3LNzi4w1qmzZVGPm2MKgy5Gk0ildcB+dL4K7WeOQwS1JxyndjPvI3CJj9SqtVuLwUYNbklYrVXCnlJidX2S8USWAw3bcknScUgX3sYUWKUGzXqVaCWbnFllsJaqVGHRpklQapZpxH5lr3w1wvDg4CXbdkrRauYK7uI3rWL3K1qbBLUlrKdWoZLbouMcaVWqV9v9TPEApSSuVKriPdnXcm+pVAA4f8yIcSepWquCeXZpx1xhrtIP7eTtuSVqhnDPuRoXxIrg7BywlSW3lCu4ipDfVqkvBPWtwS9IKpQruucUWAM16ZWlUcmTOUYkkdStVcM8vtIO7Ua0y3miP34/4pHdJWqFUwd3puOu1YKzuqESS1lKu4F7quCtUK0GjVvHgpCStUqrgnl/quNtljTeqdtyStEqpgvtYV8cNMF43uCVptVIFd6fj7gT3WKO6dDWlJKmtVME9t9CiVgkqxW1cxxs1Zj0dUJJWKFVwzy+2aNSWSxpzVCJJxylVcM8ttKhXu4K7UfU8bklapVzBvarj9qwSSTpeuYJ7IS0dmISi4za4JWmFcgX3Gh23oxJJWqlUwT2/0KJeXX4wsGeVSNLxShXcqzvuZr3K0fkWrVYaYFWSVC7lCu6F1ooZ99LDFByXSNKScgX34srTAQ1uSTpeuYJ74fgLcMDHl0lSt1IF9/zi6lGJD1OQpNVKFdzHddyN9msvwpGkZaUK7tX3Kmk6KpGk45QquI+7V0kR3N7aVZKW1XrZKSIeAw4Bi8BCSmmyH8WsPo97zLNKJOk4PQV34fdTSk/3rRLWOI+73i7PGbckLSvXqGT1jLs4OGnHLUnLeu24E/CtiEjADSmlG1fvEBHXAtcCnHvuuadUzA3vnuSs7c2l90szbjtuSVrSa3BfnlJ6IiLOBL4dEQ+klL7XvUMR5jcCTE5OntLNRX7vlRMr3i9dgGPHLUlLehqVpJSeKH7fD3wNeF0/i+qoVSs0qhWDW5K6rBvcEbE5IrZ2XgNXAff1u7COZr3iedyS1KWXUclLga9FRGf/z6eUvtnXqrr4FBxJWmnd4E4pPQL81gbUsqaxuk/BkaRupTodcC1Ng1uSVih9cI83ql7yLkldSh/cY42qV05KUpfyB3fdg5OS1K30wd1+YLDBLUkdpQ9uzyqRpJVKH9zjzrglaYXSB3ezYcctSd1KH9xj9SpzCy0WW6d03ypJGjpZBDf4+DJJ6ih9cI/7+DJJWqH0we2T3iVppdIHtw8MlqSVyh/cdtyStEI+wW3HLUlADsHdsOOWpG75BLcdtyQBOQS3M25JWiGf4LbjliQgh+BueOWkJHUrfXB3LsDxDoGS1Fb64K5XK9Sr4ahEkgqlD24onvRuxy1JQCbBPebjyyRpSRbB7VNwJGlZFsHd9LmTkrQki+AeazgqkaSOPILbg5OStCSf4LbjliQgl+Bu2HFLUkcewW3HLUlL8gjuhsEtSR15BLcHJyVpSR7B3ahybKFFq5UGXYokDVwewe09uSVpSR7B7ePLJGlJFsHd9PFlkrQki+DujEq87F2STiK4I6IaEfdExC39LGgt4w2fgiNJHSfTcb8P2NuvQk7Eg5OStKyn4I6Is4G3AJ/sbzlra3pwUpKW9NpxfxT4INB6oR0i4tqImIqIqenp6dNR25KlGbejEklaP7gj4q3A/pTS3SfaL6V0Y0ppMqU0OTExcdoKhOUZtx23JPXWcV8OvC0iHgO+CFwZEZ/ra1WrdDpuD05KUg/BnVK6PqV0dkppD/AO4DsppXf1vbIunRm3pwNKUmbncXsBjiRB7WR2Til9F/huXyo5gXq1Qq0SzrgliUw6bmjfr8QZtyTlFNx1n/QuSZBTcPsUHEkCcgpun4IjSUBGwd30gcGSBGQU3OMNO25JgoyCe8yOW5KAjIK76cFJSQIyCu6xetW7A0oSGQX3eKPKrB23JOUT3J4OKElt2QR3s17l2EKLVisNuhRJGqhsgnusc2vXBbtuSaMtn+D21q6SBOQU3A2fgiNJkFNw130KjiRBhsHtRTiSRl0+wd1wxi1JkGNw23FLGnH5BLdnlUgSkGNw23FLGnH5BLejEkkCMgrupqMSSQIyCu5xzyqRJCCj4K5XK9Qq4ahE0sjLJrjBx5dJEmQW3M1G1UveJY28rILbhylIUmbBPd6oendASSMvq+BuOuOWpLyCe6zujFuS8gruhh23JGUX3M64JY26vIK7XuWowS1pxGUX3I5KJI26rIJ7vFFlxo5b0ojLKri3NmvMLbQ4tmB4SxpdWQX3lk01AA4fXRhwJZI0OFkF99ZmHYBDBrekEbZucEdEMyLuiogfRsT9EfHhjShsLVub7Y7b4JY0ymo97HMMuDKldDgi6sAdEfGNlNL/9rm24yx33PMb/a0lqTTWDe6UUgIOF2/rxa/Uz6JeyFLHfcyOW9Lo6mnGHRHViLgX2A98O6V05xr7XBsRUxExNT09fZrLbHNUIkk9BndKaTGl9BrgbOB1EfHqNfa5MaU0mVKanJiYOM1ltjkqkaSTPKskpXQQ+C5wdT+KWU/ndEA7bkmjrJezSiYiYnvxegx4E/BAn+taU6NWYVOtwmFn3JJGWC9nlewGPhsRVdpB/+WU0i39LeuFbW3WHZVIGmm9nFXyI+DiDailJ9uaNZ53VCJphGV15SS0zyxxxi1plGUX3FuaNUclkkZadsG9dVPdm0xJGmn5BbejEkkjLsPg9qwSSaMtu+DeNlZjZm6R+cXWoEuRpIHILrh3bm4AcHDWrlvSaMouuLePd4J7bsCVSNJgZBfcO8bbN5o6YMctaURlGNztjvuAHbekEZVdcG8vOm5HJZJGVXbBvdxxOyqRNJqyC+7xRpVGteKoRNLIyi64I4Lt43UOzthxSxpN2QU3tMcldtySRlWewb257gU4kkZWnsFtxy1phOUZ3JsNbkmjK8vg3rVlE8/MzLHgjaYkjaAsg/vMrZtICZ6ZseuWNHqyDW6A/c8fG3AlkrTx8gzubU0A9h86OuBKJGnj5RncnY77kB23pNGTZXDv2uKoRNLoyjK4G7UKOzc3HJVIGklZBje0xyWOSiSNonyDe1uT/c/bcUsaPdkG91nbx9h34Migy5CkDZdtcJ+zc4xnZuaYObYw6FIkaUPlG9w7xgF4/MDsgCuRpI2VbXCfu7MI7mcdl0gaLdkG9zlLwW3HLWm0ZBvcO8brbG5UHZVIGjnZBndE8IqXbOaR6ZlBlyJJGyrb4Aa48GVb+emvDg26DEnaUNkH95PPHeU5nz8paYRkH9wAD9p1SxoheQf3S4vgfur5AVciSRsn6+DefUaTXVsa3POLg4MuRZI2zLrBHRHnRMTtEbE3Iu6PiPdtRGG9iAgmX7GTux57dtClSNKG6aXjXgD+PqX0KuBS4K8j4qL+ltW73z5vJ/sOHOHJ57yCUtJoWDe4U0pPppT+r3h9CNgLnNXvwnr1O+ftBOD7Dz094EokaWOc1Iw7IvYAFwN3rvGxayNiKiKmpqenT1N56/uNl2/jrO1j3HrfUxv2PSVpkHoO7ojYAnwFeH9K6bjTOFJKN6aUJlNKkxMTE6ezxvXq4upXv4zvP/S053NLGgk9BXdE1GmH9k0ppa/2t6ST9/bXns3cYoub7vr5oEuRpL7r5aySAD4F7E0pfaT/JZ28V+3exut/bRefvuNRnjti1y1puPXScV8OvBu4MiLuLX69uc91nbR/uPrXeXZmjg//z/2klAZdjiT1TW29HVJKdwCxAbW8KL959hn8zZUX8LHbHmLLphr/+JaLaNSyvr5Ikta0bnDn5ANvuoCZYwt86o5HuW3vfv5k8hwuPX8n501sZtfmTVQqpf//jyStK/oxVpicnExTU1On/ev26ns/nebjtz/MDx57lu7lbdlUY7xRpVYJKpWgEkEloBJBRPsMFZ1+/qlqVOwYb/Dlv7zslD43Iu5OKU32su9QddwdV7xygiteOcEzh4/x418+x2NPz/Ds7DyHjy4wc2yBxZRopURK0EqJVoJWy7l4PyT8c9Xo2Nasb8j3Gcrg7njJlk284cIz4cJBVyJJp49H7yQpMwa3JGXG4JakzBjckpQZg1uSMmNwS1JmDG5JyozBLUmZ6csl7xExDZzqzbF3AaP2HDLXPPxGbb3gmk/WK1JKPT2Fpi/B/WJExFSv1+sPC9c8/EZtveCa+8lRiSRlxuCWpMyUMbhvHHQBA+Cah9+orRdcc9+UbsYtSTqxMnbckqQTMLglKTOlCe6IuDoiHoyIhyPiukHXczIi4pyIuD0i9kbE/RHxvmL7zoj4dkQ8VPy+o+tzri/W+mBE/GHX9tdGxI+Lj30siuepRcSmiPhSsf3OiNiz4QtdQ0RUI+KeiLileD/Ua46I7RFxc0Q8UPx9XzYCa/5A8e/6voj4QkQ0h23NEfHpiNgfEfd1bduQNUbENcX3eCgirump4JTSwH8BVeBnwPlAA/ghcNGg6zqJ+ncDlxSvtwI/BS4C/hm4rth+HfBPxeuLijVuAs4r1l4tPnYXcBntRzV+A/ijYvtfAf9WvH4H8KVBr7uo5e+AzwO3FO+Hes3AZ4G/KF43gO3DvGbgLOBRYKx4/2XgPcO2ZuAK4BLgvq5tfV8jsBN4pPh9R/F6x7r1Dvo/hKL4y4Bbu95fD1w/6LpexHr+G/gD4EFgd7FtN/DgWusDbi3+DHYDD3RtfydwQ/c+xesa7auzYsDrPBu4DbiS5eAe2jUD22iHWKzaPsxrPgt4vAiWGnALcNUwrhnYw8rg7vsau/cpPnYD8M71ai3LqKTzj6NjX7EtO8WPQBcDdwIvTSk9CVD8fmax2wut96zi9ertKz4npbQAPAe8pC+L6N1HgQ8Cra5tw7zm84Fp4DPFeOiTEbGZIV5zSumXwL8AvwCeBJ5LKX2LIV5zl41Y4yllX1mCO9bYlt15ihGxBfgK8P6U0vMn2nWNbekE20/0OQMREW8F9qeU7u71U9bYltWaaXdKlwCfSCldDMzQ/hH6hWS/5mKu+8e0RwIvBzZHxLtO9ClrbMtqzT04nWs8pbWXJbj3Aed0vT8beGJAtZySiKjTDu2bUkpfLTb/KiJ2Fx/fDewvtr/QevcVr1dvX/E5EVEDzgCePf0r6dnlwNsi4jHgi8CVEfE5hnvN+4B9KaU7i/c30w7yYV7zm4BHU0rTKaV54KvA7zLca+7YiDWeUvaVJbh/AFwQEedFRIP28P7rA66pZ8WR408Be1NKH+n60NeBzlHia2jPvjvb31EcaT4PuAC4q/hx7FBEXFp8zT9d9Tmdr/V24DupGIoNQkrp+pTS2SmlPbT/vr6TUnoXw73mp4DHI+LCYtMbgZ8wxGumPSK5NCLGi1rfCOxluNfcsRFrvBW4KiJ2FD/dXFVsO7GNPgBwggMDb6Z9NsbPgA8Nup6TrP31tH+8+RFwb/HrzbRnWLcBDxW/7+z6nA8Va32Q4shzsX0SuK/42MdZvrq1Cfwn8DDtI9fnD3rdXTW/geWDk0O9ZuA1wFTxd/1ftM8EGPY1fxh4oKj3P2ifTTFUawa+QHuGP0+7C/7zjVoj8GfF9oeB9/ZSr5e8S1JmyjIqkST1yOCWpMwY3JKUGYNbkjJjcEtSZgxuScqMwS1Jmfl/L1mIquoTiucAAAAASUVORK5CYII=", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" } ], "source": [ "V = Ausgleichsbecken_class(1.,0.5,0.,10.,timestep=0.001)\n", - "V.set_initial_level(5.) \n", - "V.set_influx(1.)\n", - "V.set_outflux(0.)\n", + "V.set_initial_level(initial_level) \n", + "V.set_influx(initial_influx)\n", + "V.set_outflux(initial_outflux)\n", "\n", - "V.p0 = 0.\n", + "V.p0 = initial_pipeline_pressure\n", "\n", "outflux_vec = []\n", "level_vec = []\n", "\n", "t = 0\n", - "while V.level > 0.005:\n", + "while V.level > total_min_level:\n", " t = t+V.timestep\n", " V.e_RK_4()\n", " V.level = V.update_level(V.timestep)\n", " V.set_volume()\n", " outflux_vec.append(V.outflux)\n", " level_vec.append(V.level)\n", - " \n", + " if t > total_max_time:\n", + " break\n", + "\n", "plt.plot(level_vec)\n" ] }, { "cell_type": "code", - "execution_count": 31, + "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/plain": [ - "[]" + "[]" ] }, - "execution_count": 31, + "execution_count": 11, "metadata": {}, "output_type": "execute_result" }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAD8CAYAAABw1c+bAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAk9ElEQVR4nO3deZyU1Z3v8c+vqvem6YamWWRfVFxR7LCoUYkjUZMZzJhFjWbTMEw0jpPJTfSVyczkZu6d69y8vJlMNIRJnHGJGpOIIXE3Go2DKI0BQQVFlgCNdrMKNNDdVb/7Rz0FRdNNPU1Xd3U//X2/rFc9dZ5zqs4B/NWp8zznHHN3RESk/4jluwIiItKzFPhFRPoZBX4RkX5GgV9EpJ9R4BcR6WcU+EVE+plQgd/MNpjZSjNbbmZ1x8j3ITNLmNknM9IuNbM1ZrbWzG7NRaVFROT4FXQi7yx339bRSTOLA7cDT7VJuxO4BNgMLDWzRe7+5nHWV0REuiiXQz1fBX4FNGSkTQPWuvs6d28GHgLm5PAzRUSkk8L2+B142swc+LG7L8g8aWYjgU8AHwE+lHFqJLAp4/VmYHq2DxsyZIiPGzcuZNVERGTZsmXb3L0mTN6wgf88d683s6HAM2a22t1fzDj/feCb7p4ws8xyR7wItLtGhJnNBeYCjBkzhrq6Di8liIhIG2a2MWzeUIHf3euD5wYzW0hqCCcz8NcCDwVBfwhwuZm1kurhj87INwqo7+AzFgALAGpra7WAkIhIN8ka+M2sHIi5+57geDbwPzPzuPv4jPz/BfzW3R81swLgRDMbD2wBrgKuyWH9RUSkk8L0+IcBC4PefAHwgLs/aWbzANx9fkcF3b3VzG4idadPHLjb3d/oerVFROR4WW9clrm2ttY1xi8iEp6ZLXP32jB5NXNXRKSfUeAXEelnFPhFRPqZyAf+1zfv4t6XN7C/OZHvqoiI9AqdWaunz2ncc5CrFiyhqTnB79c08tPP19JmgpmISL8T6R7/ohX1NDUnuHraaJ5b3cAzb76f7yqJiORdpAP/knXbGVtdxnfnnM6EmnLueOZteuPtqyIiPSnSgf+NLbs5e3QVBfEY8y6cyOr39vDyu9vzXS0RkbyKbOA/0JKgfvcBxg8ZAMBfTDmB6vIi/nPxhvxWTEQkzyIb+DftaAJg3JAyAEoK41wzfQzPvvX+oXMiIv1RZAP/xu2p4D5mcNmhtGumjyFmxv2vhF69VEQkciIb+Bv3HgRg2MCSQ2kjKku55JRhPLx0EwdadF+/iPRPkQ38O/Y1AzC4vOiI9M/NHMvOphYee31rPqolIpJ3kQ382/c2U14Up6QwfkT6zInVTKgp574lGu4Rkf4psoF/x76DDB5QdFS6mXHdjLEs37SLlZt356FmIiL5FdnAv31fM4PLi9s9d+U5oygtjHPfkg09WykRkV4gsoF/x75mqsuP7vEDDCwp5IqzR/Lr5fXsbmrp4ZqJiORXqMBvZhvMbKWZLTezo7bGMrM5ZvZ6+ryZnR+2bHfZvb+FqtLCDs9fN2MsB1uT/GLZpp6qkohIr9CZ1Tlnufu2Ds79Dljk7m5mZwIPA5NDlu0W+w62Ul7ccfNOPWEgtWMHcf+SjXzpvPHEYlq1U0T6h5wM9bj7Xj+8+lk5kPeV0PYdTBwz8ANcN3MsG7Y38Ye1PfqdJCKSV2EDvwNPm9kyM5vbXgYz+4SZrQYeA77UmbK51tyapDmRZEBx/Jj5Lj19OEMGFHHfy7q1U0T6j7CB/zx3nwpcBtxoZhe0zeDuC919MnAF8N3OlAUws7nB9YG6xsbGTjWirX0HWwGy9viLC+J85kOjeW71+2zeqfV7RKR/CBX43b0+eG4AFgLTjpH3RWCimQ3pTFl3X+Dute5eW1NT06lGtLU3ZOAHuGb6WAAeeOVPXfpMEZG+ImvgN7NyM6tIHwOzgVVt8kyyYE9DM5sKFAHbw5TtDvuaU4F/QIjAP7KqlItPGcbPl27iYKvW7xGR6AvT4x8GvGRmK4BXgcfc/Ukzm2dm84I8VwKrzGw5cCfwmeBib7tlc96KNsIO9aRdN2Ms2/c188TK97qzWiIivULWyOju64Ap7aTPzzi+Hbg9bNnutvdgquee7eJu2vmThjCuuoz7lmzkirNHdmfVRETyLpIzdzvb44/FjGtnjGXZxp28Ua/1e0Qk2qId+IvCz0/71DmjKSmM8V//vaGbaiUi0jtEMvDvDzZZKS0KN9QDUFlWyKdrR/Po8i28t/tAd1VNRCTvIhn4D7YkAY5aiz+bL394AkmHn760rjuqJSLSK0Qz8Ae3ZRYXdK55oweX8bEzRvDAK3/Sqp0iElkRDfxJYgYFx7Hw2l9dOIF9zQltyC4ikRXZwF9SGCeYU9Ypp51QyYUn1XD3S+u1IbuIRFIkA/+BlkSnh3kyzbtwItv3NfOLZZtzWCsRkd4hkoH/YEuS4oLOXdjNNGPCYM4aXcWPX3iX5tZkDmsmIpJ/0Qz8rQmKC4+/aWbGzRdPYvPO/fzqNfX6RSRaIhr4k10a6gGYdfJQpoyu4ofPrVWvX0QiJcKB//iHeiDV67/lz05ky6792pdXRCIlooG/axd30y46qYazRldx53NrtWSziERGJAP/gZZkl8b408yMv73kJOp3H+DhOo31i0g0RDLwH2xNUNLFoZ60C04cwtQxVdz1/Frd1y8ikRDNwJ+jHj+kev1f/+jJbN19gHsWb8jJe4qI5FM0A38OLu5mOnfiEC46uYY7n1/LrqbmnL2viEg+hAr8ZrbBzFaa2XIzq2vn/Bwzez193szOzzh3qZmtMbO1ZnZrLivfkVxd3M1062WT2XOwlTufX5vT9xUR6WmdiY6z3P0sd69t59zvgCnufhbwJeAnAGYWJ7UH72XAqcDVZnZq16qcXS7u429r8vCBXDl1FPcs3simHU05fW8RkZ6Uk+jo7nuDzdUByoH08TRgrbuvc/dm4CFgTi4+81gOtCQo7uRa/GF87ZKTMIM7nnk75+8tItJTwgZ+B542s2VmNre9DGb2CTNbDTxGqtcPMBLInP20OUjrVs2tSYriub98cUJVKV88bzyPLt/Cqi3am1dE+qaw0fE8d59KasjmRjO7oG0Gd1/o7pOBK4DvBsntrYvs7aRhZnOD6wN1jY2NIat1tGTSSToUxDu/JHMYf33RRAaXFfGPi97g8I8cEZG+I1Tgd/f64LkBWEhqCKejvC8CE81sCKke/uiM06OA+g7KLXD3WnevrampCVn9o7UkU+vqFHZDjx+gsrSQb146mWUbd7Lwj1u65TNERLpT1uhoZuVmVpE+BmYDq9rkmWTBridmNhUoArYDS4ETzWy8mRUBVwGLctuEI7UkUr3wwm7q8QN88pxRTBldxb88sZo9B7RFo4j0LWG6xcOAl8xsBfAq8Ji7P2lm88xsXpDnSmCVmS0ndRfPZzylFbgJeAp4C3jY3d/IeSsytCZSPf6CWPdNUYjFjO/8xWk07jnIvz+n2ztFpG8pyJbB3dcBU9pJn59xfDtwewflHwce70IdO6UnevwAZ42u4tO1o7j7pfV8unY0k4YO6NbPExHJlcjN3G0NxvgLummMP9M3Lp1MaVGcbz+6Shd6RaTPiFzgb2lN9/i7v2lDBhRz22Wn8PK67fxCq3eKSB8RvcB/6K6e7h3qSbvqQ6OZNn4w//zYmzTsOdAjnyki0hWRC/ytwRh/d17czRSLGf/yl2dwoDXJd37zZo98pohIV0Qu8Lek7+rpoR4/wMSaAdz8kUk89vpWnn3z/R77XBGR4xHZwN8dSzYcy9wLJjJ5eAXfenQlu5t0b7+I9F6RC/ytyWCopwd7/ABFBTG+96kpbN/bzLd/vSp7ARGRPIlc4G/pgQlcHTl9ZCV/c/GJLFpRz6IV7a5MISKSd5EL/K09NIGrI3990UTOHlPF3y9cyXu7dZePiPQ+kQv86R5/T9zH356CeIw7Pn0WLQnnf/xyBcmkJnaJSO8SwcCfnzH+TOOHlPOtj53CH97Zxk9fWp+3eoiItCdygb+1m5dlDuuz08dw6WnDuf3J1SzbuDOvdRERyRS9wH9oAlf+evwAZsbtnzyTEVUlfPWB19i5rzmv9RERSYtc4G/O8xh/psrSQu665hy27W3maw8v13i/iPQK+Y+OOXb4rp7e0bQzRlXy9x8/hefXNPKjF97Nd3VERCIY+JM9v2RDNtfNGMvHzxzB955ew3OrtaSDiORX5AL/oY1Y8jCBqyNmxr9+8kxOHTGQv3lwOWsb9uS7SiLSj4WKjma2wcxWmtlyM6tr5/xnzez14LHYzKaELZtrh+7jL+g9PX6AsqICFnyuluLCGF++d5nW8xGRvOlMt3iWu5/l7rXtnFsPXOjuZwLfBRZ0omxO9cSeu8drZFUp8689h807m7jpwdcO1VVEpCflJDq6+2J3T9+svgQYlYv3PR49tefu8aodN5h/vuJ0/vDONv5eWzaKSB6EDfwOPG1my8xsbpa81wNPHGfZLmtNJonHDLPeGfgBPvOhMXz1I5N4aOkm/u137+S7OiLSzxSEzHeeu9eb2VDgGTNb7e4vts1kZrNIBf7zj6PsXGAuwJgxYzrdkLSWhPfa3n6mr11yElt3H+D7z77DsIElXD3t+NssItIZoXr87l4fPDcAC4FpbfOY2ZnAT4A57r69M2WD8wvcvdbda2tqajrbjkNaEsledUdPR8xSWzZeeFIN31q4kme0c5eI9JCsEdLMys2sIn0MzAZWtckzBngEuM7d3+5M2VxrTXivuof/WArjMe767FTOGFnJjT97jRffbsx3lUSkHwjTNR4GvGRmK4BXgcfc/Ukzm2dm84I8/wBUA3e1uW2z3bI5bsMRWpNOvA/0+NPKiwu450vTmDh0AF++t46X392evZCISBdkHeN393XAlHbS52cc3wDcELZsd0omPe8LtHVWVVkR918/jav/YwnX37OUe780jdpxg/NdLRGJqL7TNQ4p1ePvW4EfoHpAMfffMJ3hA0v4wn8uZemGHfmukohEVOQCf9KdPjTSc4ShFSU88OUZDB1YzHU/fYUXNOYvIt2gj4bIjiWS3itn7YY1vLKEh/9qJhOGDOCGe5by5Kqt+a6SiERM342QHUgknT440nOEIQOKeXDuDM4YWclXfvYav1y2Od9VEpEIiWTg78s9/rTK0kLuu346504cwtd/sYJ/e/YdLe8gIjnR9yNkG61JJ9bXu/yB8uIC7v7Ch7hy6ij+37Nv83e/WEFzqxZ2E5GuCbtkQ5+R9L53O+exFBXE+N6nzmRsdRl3PPM29bv2M//ac6gqK8p31USkj1KPvw8wM26++ES+/5mzeG3jLubc+d+8tfWDfFdLRPqoyAX+vjiBK6wrzh7Jg3NncKAlwSfu+m8e/eOWfFdJRPqgyAX+1mSSeC9ekrmrzhk7iN989XzOHFXFLT9fzj8tekPj/iLSKZEL/MkkfXLmbmcMrSjhZzdM5/rzx/NfizfwqfmLWb9tX76rJSJ9ROQCf3ojlqgrjMf49sdPZf61U9mwvYmP/eAPPLx0k275FJGsIhf4Ex79Hn+mS08fwZO3fJgpo6r4xq9e5ys/e42d+5rzXS0R6cUiF/iTfXSRtq4YUVnK/TdM59bLJvPsW+/zZ3e8wG9W1Kv3LyLtilzg76urc3ZVPGbMu3Aii246n5GDSvnqg3/ky/cu473dB/JdNRHpZSIX+JNJj/RdPdmcMmIgj/z1uXzr8lN4aW0jl9zxAvcs3kBrQnf+iEhK5AJ/azJJvI9svdhdCuIxvnzBBJ665QKmjK7iHxe9wcd+8BKL392W76qJSC8QKvCb2QYzW9lmW8XM8581s9eDx2Izm5Jx7lIzW2Nma83s1lxWvj1Jp1/3+DONrS7nvuunMf/aqexrbuWa/3iFr/xsGZt2NOW7aiKSR51Zq2eWu3fUZVwPXOjuO83sMmABMN3M4sCdwCXAZmCpmS1y9ze7VOtjaE0mIztz93iYGZeePoKLTh7KghfXcdfv1/Lsmw1cM30MN86aRE1Fcb6rKCI9LCdDPe6+2N13Bi+XAKOC42nAWndf5+7NwEPAnFx8ZkeSSSK3Vk8ulBTGufniE3n+6xdx5TkjuW/JRi78v8/zvafWsHt/S76rJyI9KGzgd+BpM1tmZnOz5L0eeCI4Hglsyji3OUjrNlFfsqGrRlSW8i9/eSbP/O0FXHzKMH74/Fo+fPtz3PH0GrbvPZjv6olIDwgb+M9z96nAZcCNZnZBe5nMbBapwP/NdFI72dq9udzM5ppZnZnVNTYe/16ziST9/uJuGBNqBvDvV5/NYzefz8yJ1fzgubWcd/tz/NOiN9iya3++qyci3ShU4Hf3+uC5AVhIagjnCGZ2JvATYI67bw+SNwOjM7KNAuo7+IwF7l7r7rU1NTXhW9BGQj3+TjnthEp+fF0tz37tAj5+5gncv2QjF/7r89zy0B957U87NQlMJIKyBn4zKzezivQxMBtY1SbPGOAR4Dp3fzvj1FLgRDMbb2ZFwFXAolxVvj2JfjqBq6smDa3ge5+awovfmMV1M8fy7FsN/OVdi/nzH77Ew3WbONCSyHcVRSRHwtzVMwxYaKledAHwgLs/aWbzANx9PvAPQDVwV5CvNei9t5rZTcBTQBy4293f6IZ2HJLsZ2v15NoJVaX845+fxt/NPpmFf9zCvYs38I1fvs7/fvwtPnH2SK6cOorTThiI6VeVSJ9lvfGnfG1trdfVHTVdIJTJ336Cz88cx22Xn5LjWvVP7s6SdTu4b8kGnn2zgeZEksnDK7hy6ijmnH0CQytK8l1FEQHMbJm714bJG709d3U7Z06ZGTMnVjNzYjW7mpr5zYp6fvnaFv7X42/xf55czbkTq7n8jBHMPnUY1QM0J0CkL4hc4NcEru5TVVbEdTPHcd3Mcaxt2MOvXtvCY69v5bZHVvKthSuZMaGay84YwUdPHcbQgfolINJbRSrwuztJh5jGn7vdpKEVfPPSyXzjoyfz5tYPeGLlezy+aivffnQV3350FaePHMhFJw3lopNrOGt0FQXxyC0LJdJnRSrwJ5Kp6xXq8fccM+O0Eyo57YRK/m72SbzTsJdn3nyf369p4EcvvMsPn1/LwJICPnxSDR+eNIQZE6oZW12mi8MieRStwB9cqNYYf36YGScNq+CkYRXcOGsSu/e38NI72/j9mgZeeLuRx17fCsCwgcXMmFDN9PHVzJgwmPFDyvVFINKDohX41ePvVSpLC/nYmSP42JkjcHfebdzHknXbeWX9Dha/u51fL0/N5asuL2LK6CqmjKpiyuhKzhpdRVVZUZ5rLxJdkQz8uo+/9zEzJg0dwKShA7h2xljcnfXb9vHK+h28tnEnyzft4vk1DaTvLh5XXcaU0VWcOmIgk0cM5JThFdRUFOuXgUgOKPBLXpgZE2oGMKFmAFdPGwPAngMtrNyym+WbdrFi0y5eXb/j0K8CgEFlhUwePpDJIyo4ZfhAJg4tZ/yQAQwu168Dkc5Q4Jdeo6KkkHMnDuHciUMOpe1qamb1e3tYvfUD1ry/h7e27uGhVzexP2MJiaqyQsYPKWf8kHImDClnQs0AxlWXM3JQKZWlhfloikivFq3A7wr8UVNVVsSMCdXMmFB9KC2ZdDbtbGJd4z7WbdvHusa9rN+2j8Vrt/PIa1uOKF9RUsDIqlJGDSpj1KBSRg0qZWRVKSMHlTK8soTq8mL9e5F+J1qBP93j1zhwpMVixtjqcsZWlzOrzbl9B1tZv20fG7c3sWVXE5t37mfLzv1s2tHEknXb2Xuw9Yj88ZhRXV7E0IHFDKsoYejAYmoqShhaUczQimJqKooZXF7EoPIiKooLdI1BIiGagV89uH6rvLiA00dWcvrIyqPOuTu797ekvgx27ef9Dw7Q8MFBGvYcoGHPQep3H2DF5l1s39dMe0tYFcSMQeVFDC4roqqs8NAXQvp1ZWkhFSWFDCwpoKKkkAElBVQEj+KCeA+0XiQcBX7pN8yMqrIiqsqK2v1iSGtJJNm+t5mGPQdo3HOQnU0t7NzXzI6mZnY1NbNjXzM797XwTsNedjU1s7Op5dC/vY4UFcQOfSGkvwzKiwooK4pTWhSntPDwcVnwKCmMU1Z0ZHppYeq4OB6nqCBGUUFM/96l0xT4RdoojMcYXlnC8Mpw6w0lk86eA618cKCFDw60sOdAa/BoOeL5gzZp2/c209ScoKk5wf7mVppaEu3+0simIGYUFcQoDr4IUsdxiuIxigtjFMUPpx3KE49REDcK46kvjoK4URAzCmIxCuNG/NCzURCPURhLHafzF8ZTeeNxozB2OC0eM2KWejbj0OvUI+N1LHhthlm6HEG6BemZ5dEwWw4p8It0USxmVJYVUlnWtTuI3J2DrUn2Nydoagm+DJoTGa+DL4mWBM2tSZpbkxxszTxOPTcnDqcfDB57DrSyvbU5lZ5I5WtNOK1JpzWRpCXpJIJHb2XBF0Xqi4NDXxCk/sOCL4tDx0EZOJwea5Mn/b6ptKPLZ753LCMPbT6jbXkyzx3xuRmxyY54wiw16fHH14VaWblLohX4XRd3pe8yM0oKU0M8g/JUh2TSSbjTmnBakkkS6edkkJZIHbckUl8SLcn0F0jqOeFOMplaLDGRdNyDNCdI9yCdIP3I/Ml0Wmb5JBnpqdcevE/C/dCvJHfHAXdwPHgmOB+89tR7ZebjUL7M9MOvyXi/pHfwvh2UJ+NzksnDf87BOx+uO+m6deNfboZIBf7WhHr8Il0RixkxjMI4lKIL0lEVaq1cM9tgZivNbLmZHbU1lplNNrOXzeygmX29M2VzKan7+EVEsupMj3+Wu2/r4NwO4GbgiuMomzOtGuMXEckqJ7tjuHuDuy8FWnLxfscrqcAvIpJV2MDvwNNmtszM5nbyM7pStlM0c1dEJLuwQz3nuXu9mQ0FnjGz1e7+Yi7LBl8KcwHGjBkT8q2PlL4irvt9RUQ6FqrH7+71wXMDsBCYFvYDwpZ19wXuXuvutTU1NWHfvu17ABrqERE5lqyB38zKzawifQzMBlaFefOulD0eh7ZeVNwXEelQmKGeYcDCYPikAHjA3Z80s3kA7j7fzIYDdcBAIGlmtwCnAkPaK5vzVgTSQz3ac1dEpGNZA7+7rwOmtJM+P+P4PWBUO8U/aK9sd0ke6vEr8IuIdCQnt3P2FunbOdXhFxHpWLQCf3qoRz1+EZEORSrwJ5Ia6hERySZSgT99O2csUq0SEcmtSIXI9FCPZu6KiHQsUoE/fR+/Zu6KiHQsUoFfM3dFRLKLVOBPauauiEhWkQr8iWBrM93VIyLSsUgF/kM9fnX5RUQ6FK3Ar5m7IiJZRSvwa+auiEhWEQv8mrkrIpJNRAN/nisiItKLRSvwa60eEZGsIhX4E9qIRUQkq0gFftdQj4hIVqECv5ltMLOVZrbczOraOT/ZzF42s4Nm9vU25y41szVmttbMbs1VxduT1JINIiJZhdlzN22Wu2/r4NwO4GbgisxEM4sDdwKXAJuBpWa2yN3fPI66ZqWZuyIi2eVkqMfdG9x9KdDS5tQ0YK27r3P3ZuAhYE4uPrM9yUOrc3bXJ4iI9H1hA78DT5vZMjOb24n3Hwlsyni9OUg7ipnNNbM6M6trbGzsxEdkVDI91KPILyLSobCB/zx3nwpcBtxoZheELNdeBPb2Mrr7AnevdffampqakG9/JA31iIhkFyrwu3t98NwALCQ1hBPGZmB0xutRQH1nKtgZGuoREckua+A3s3Izq0gfA7OBVSHffylwopmNN7Mi4Cpg0fFWNpukOzHTDlwiIscS5q6eYcDCIJgWAA+4+5NmNg/A3eeb2XCgDhgIJM3sFuBUd//AzG4CngLiwN3u/kY3tANIB34FfRGRY8ka+N19HTClnfT5GcfvkRrGaa/848DjXahjaEnXrF0RkWwiNXM3mXTN2hURySJagV9DPSIiWUUq8CeSuodfRCSbSAX+pLtu5RQRySJSgd/ddXFXRCSLSAX+hLuGekREsohU4E+6Jm+JiGQTqcDv7sQj1SIRkdyLVJhMJHU7p4hINpEK/EnXypwiItlEK/AnnVikWiQiknuRCpOauSsikl3EAr+GekREsolU4E+4FmkTEckmUoHfNdQjIpJVpAJ/MqmhHhGRbCIV+BNaq0dEJKswWy9iZhuAPUACaHX32jbnDfg34HKgCfiCu78Wpmwuucb4RUSyChX4A7PcfVsH5y4DTgwe04EfBc9hyuZMIunEFflFRI4pV0M9c4B7PWUJUGVmI3L03qFpkTYRkezCBn4HnjazZWY2t53zI4FNGa83B2lhygJgZnPNrM7M6hobG0NW60hJDfWIiGQVdqjnPHevN7OhwDNmttrdX8w431649ZBlU5ndFwALAGpra73t+TCSWo9fRCSrUD1+d68PnhuAhcC0Nlk2A6MzXo8C6kOWzRndzikikl3WwG9m5WZWkT4GZgOr2mRbBHzOUmYAu919a8iyOZPQnrsiIlmFGeoZBiwMLpoWAA+4+5NmNg/A3ecDj5O6lXMtqds5v3issjltQQZ3p0A7sYiIHFPWwO/u64Ap7aTPzzh24MawZbuLFmkTEckuUt3jRFIzd0VEsolU4NfMXRGR7CIV+DXUIyKSXaQCvzZbFxHJLlKBXzN3RUSyi2DgV+QXETmWiAV+tDqniEgWEQv8mrkrIpJNtAK/Lu6KiGQVrcCvoR4RkawiFvg11CMikk20An9S6/GLiGQTrcCvmbsiIllFKvAn3IlFqkUiIrkXqTDpmsAlIpJVpAK/hnpERLKLVOBPLdKW71qIiPRuoQK/mW0ws5VmttzM6to5b2b2AzNba2avm9nUjHOXmtma4Nytuax8W0nXRiwiItmE2XM3bZa7b+vg3GXAicFjOvAjYLqZxYE7gUuAzcBSM1vk7m92oc4d0sxdEZHscjXUMwe411OWAFVmNgKYBqx193Xu3gw8FOTtFpq5KyKSXdjA78DTZrbMzOa2c34ksCnj9eYgraP0o5jZXDOrM7O6xsbGkNU60kdPG8YpIyqOq6yISH8RdqjnPHevN7OhwDNmttrdX8w4314324+RfnSi+wJgAUBtbW27ebL5/lVnH08xEZF+JVSP393rg+cGYCGpIZxMm4HRGa9HAfXHSBcRkTzJGvjNrNzMKtLHwGxgVZtsi4DPBXf3zAB2u/tWYClwopmNN7Mi4Kogr4iI5EmYoZ5hwEJL3S1TADzg7k+a2TwAd58PPA5cDqwFmoAvBudazewm4CkgDtzt7m/kvBUiIhKauR/XcHq3qq2t9bq6o6YLiIhIB8xsmbvXhskbqZm7IiKSnQK/iEg/o8AvItLPKPCLiPQzvfLirpk1AhuPs/gQoKM1haJKbY6+/tZeUJs7a6y714TJ2CsDf1eYWV3YK9tRoTZHX39rL6jN3UlDPSIi/YwCv4hIPxPFwL8g3xXIA7U5+vpbe0Ft7jaRG+MXEZFji2KPX0REjiEygb8n9/bNNTMbbWbPm9lbZvaGmf1NkD7YzJ4xs3eC50EZZW4L2rrGzD6akX5OsD/y2mAfZAvSi83s50H6K2Y2rscb2g4zi5vZH83st8HrSLfZzKrM7Jdmtjr4+57ZD9r8t8G/61Vm9qCZlUStzWZ2t5k1mNmqjLQeaaOZfT74jHfM7POhKuzuff5BauXPd4EJQBGwAjg13/XqRP1HAFOD4wrgbeBU4F+BW4P0W4Hbg+NTgzYWA+ODtseDc68CM0ltgvMEcFmQ/hVgfnB8FfDzfLc7qMvXgAeA3wavI91m4B7ghuC4CKiKcptJ7bi3HigNXj8MfCFqbQYuAKYCqzLSur2NwGBgXfA8KDgelLW++f4fIUd/6DOBpzJe3wbclu96daE9vya1Qf0aYESQNgJY0177SC17PTPIszoj/Wrgx5l5guMCUpNELM/tHAX8DvgIhwN/ZNsMDCQVBK1NepTbnN5+dXBQn9+S2tMjcm0GxnFk4O/2NmbmCc79GLg6W12jMtQTem/f3i74CXc28AowzFMb2hA8Dw2yHWuP483tpB9Rxt1bgd1Adbc0IrzvA98AkhlpUW7zBKAR+M9geOsnltrcKLJtdvctwPeAPwFbSW3S9DQRbnOGnmjjccW+qAT+0Hv79mZmNgD4FXCLu39wrKztpGXb47hX/RmZ2ceBBndfFrZIO2l9qs2kempTgR+5+9nAPlJDAB3p820OxrXnkBrSOAEoN7Nrj1WknbQ+1eYQctnG42p7VAJ/n9/b18wKSQX9n7n7I0Hy+2Y2Ijg/AmgI0o+1x/GodtKPKGNmBUAlsCP3LQntPOAvzGwD8BDwETO7n2i3eTOw2d1fCV7/ktQXQZTb/GfAendvdPcW4BHgXKLd5rSeaONxxb6oBP4+vbdvcOX+p8Bb7n5HxqlFQPoq/edJjf2n068KrvSPB04EXg1+Tu4xsxnBe36uTZn0e30SeM6DQcF8cPfb3H2Uu48j9ff1nLtfS7Tb/B6wycxODpIuBt4kwm0mNcQzw8zKgrpeDLxFtNuc1hNtfAqYbWaDgl9Xs4O0Y+vpCyDdeGHlclJ3w7wLfCvf9elk3c8n9fPsdWB58Lic1Bje74B3gufBGWW+FbR1DcGV/yC9FlgVnPshhyfplQC/ILUv8qvAhHy3O6POF3H44m6k2wycBdQFf9ePkroTI+pt/g6wOqjvfaTuZolUm4EHSV3DaCHVC7++p9oIfClIXwt8MUx9NXNXRKSficpQj4iIhKTALyLSzyjwi4j0Mwr8IiL9jAK/iEg/o8AvItLPKPCLiPQzCvwiIv3M/wcgaaxUFozYfAAAAABJRU5ErkJggg==", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD4CAYAAAD8Zh1EAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAh2UlEQVR4nO3deXhV1b3/8feXhATCkAGSEBJCCDMymyKIRRRnbR2x0quiYtEO3mq9rWhbO9hr6a/Wtuptb60TzuII4kVFcKgTyAwhDGEKIXMgA0NChvX7I9uWIpAEkuyccz6v58lz9l7nnJzvCuGTddZeex9zziEiIoGvg98FiIhIy1Cgi4gECQW6iEiQUKCLiAQJBbqISJAIb8sX69mzp0tLS2vLlxQRCXgrVqwocc7FN/a4Ng30tLQ0li9f3pYvKSIS8MxsZ1MepykXEZEgoUAXEQkSCnQRkSChQBcRCRIKdBGRIKFAFxEJEgp0EZEg0abr0EVEQkVNXT3bS/azqaCSTQWVXDOuDymxUa36mgp0EZGT4Jwjr7yKTQUVbCyoZHNBJRsLKtlWvJ9DdfUAhHUwTu0bq0AXEWkvqmrq2FK4j8y8cjbkV7Ahr4JNhZVUVtX+8zG9ozsxuFc3Jg9OYHCvrgxO7E7/hC5Ehoe1en0KdBGRoyg7cOifob0hr4LMvAqyi/dRV9/wKW9dIsIYmtSdS0f3ZnCv7gzp1Y1Bid2I7tzRt5oV6CIS0pxz5JdXsW53eUN4eyG+u+zgPx+T2D2SYUndOWdYAqf0jmZYUndS46Lo0MF8rPyrFOgiElJK91WzNrecNbllrM0tZ21uOSX7qgEwg/SeXRjbN5brJvRlWFJ3hvXuTs+ukT5X3TQKdBEJWhVVNazPLWdNbjlrvQD/cuRtBgPiuzJpUE9GpcQwPDmaoUndiIoI3FgM3MpFRA5TU1fPxvxKVubsZVXOXtbmlrOtZP8/70+Ni2JMagw3nJ7GiJRohidH0zUyuCKw0d6Y2WDgpcOa0oF7gae99jRgB3C1c25vy5coIvJVpfuqWZlTxsqcvazcuZc1uWVU1TQsE0zsHsmolBiuGJvMyJQYRqZEExMV4XPFra/RQHfObQJGA5hZGLAbeB2YBSx2zs02s1ne/l2tV6qIhKq6esemgkpW5Oxl1c69rMzZy47SAwCEdzBOSY5m2rhUxqbGcmrfWHrHdPa5Yn809/3GFGCrc26nmV0KTPba5wAfoEAXkRZQVVPHml1lLNu+h2U79rBy5172H6oDoGfXSMamxjQEeN9YRiRH06lj66/xDgTNDfRrgBe87UTnXD6Acy7fzBKO9gQzmwnMBEhNTT3ROkUkiO2rrmXlzr0NAb59D6t3lXGorh4zGJzYjSvGpnBq31jvbMvOmLWv5YLthTnnmvZAswggDzjFOVdoZmXOuZjD7t/rnIs93vfIyMhw+kxRESk7cIgvduxl2fZSlm3fw/q8CurqHWEdjOHJ0ZzWL45xaXFkpMWGxNx3Y8xshXMuo7HHNWeEfiGw0jlX6O0XmlmSNzpPAopOpFARCX77q2tZtmMPn2aX8El2KVkFFTgHEeEdGN0nhu9N7s+4fnGMTY2lS5CtPGlLzfnJTeNf0y0A84HpwGzvdl4L1iUiAexQbT1rcsv4eEsJn24tYfWuMmrqHBFhHRjbN4Y7zhnE+PQejEzR/HdLalKgm1kUcC5wy2HNs4G5ZjYDyAGmtnx5IhII6usdWQUVfJpdyidbS1i2fQ8HDtVhBiOSo5lxRjoTB/Qgo28cnSMU4K2lSYHunDsA9DiirZSGVS8iEoJK91Xzjy0lfLi5mI82F1O6/xAA6fFduHJsChMH9GB8eg/NgbchTVaJSJPU1TtW7yrjw01FfLi5mLW7y3EO4rpEMGlgT84YGM/EAT1Iig7NNeDtgQJdRI6pqKKKDzcX88HmYj7eUkL5wRo6GIzuE8PtUwYxeXA8w5OjCWtnVx0MVQp0Efkn5xyZeRUs2lDI4o2FrN9dAUB8t0jOHZbImYPi+frAnppGaacU6CIhrrq2js+2lvJeViGLs4rIL6/CDMamxvLj8wczeXA8w5K662SeAKBAFwlBe/YfYsnGIhZnFfLR5mL2H6qjc8cwJg3qyY/OHcRZQxIC5hrg8i8KdJEQsbvsIG+vL+Dt9fms2LmXetdwVcLLxiRzzrBEJqT30JrwAKdAFwliO0v3s3B9AQvXF7BmVxkAQ3p14wdnD+TcoYkMT9ZUSjBRoIsEma3F+1i4Lp+F6wvIzGs4qDkyJZq7LhjChcN7kdazi88VSmtRoIsEgeyiSt5ck8/C9flsLtwHwKl9Y/nZxUM5/5Re9ImL8rlCaQsKdJEAlVd2kDfX5DFvdR4b8iswg3Fpcfzqm6dw/im96BXdye8SpY0p0EUCyJ79h3hrXT5vrs5j2Y49QMNJPr/4xjAuHplEQjeFeChToIu0c/uqa1m0oYB5q/P4eEsJtfWOgQld+a/zBvGNUb3p20Nz4tJAgS7SDtXXOz7dWsqrK3NZuD6fqpp6kmM6c/PX07l0dG+G9Oqm1SnyFQp0kXZke8l+Xl2Ry2src8krr6Jbp3CuGJvCFWOSGZsaSwddM0WOQ4Eu4rOKqhreWpvPKytyWbFzLx0MJg2K556Lh3LO0ESd7CNNpkAX8UF9veOzbaXMXb6Lt9cXUF1bz4CErsy6cAiXj0kmsbsObkrzKdBF2lDJvmpeWZHLi8ty2FF6gOjOHbk6ow9XnZrCyJRozYvLSVGgi7SyL0fjzy/N4d0NBdTUOcb1i+P2cwZxwfBemlKRFqNAF2klJfuqeXl5Li9+kcPO0gPERHXk+glpTBvXhwEJ3fwuT4KQAl2kBTnnWLFzL099uoN3Mv81Gr9Do3FpA00KdDOLAR4DhgMOuAnYBLwEpAE7gKudc3tbo0iR9q6qpo75a/KY8+kOMvMq6NYpnOvGp/Ht01IZkNDV7/IkRDR1hP5n4G3n3FVmFgFEAfcAi51zs81sFjALuKuV6hRpl3aXHeTZz3fy4rIc9h6oYVBiV/778uFcPiaZqAi9AZa21ehvnJl1ByYBNwA45w4Bh8zsUmCy97A5wAco0CUEOOf4fNse5ny6g3c3FABw7rBEpp+exoT0HlqpIr5pyhAiHSgGnjSzUcAK4IdAonMuH8A5l29mCUd7spnNBGYCpKamtkjRIn6oqavnrbX5PPrRNjbkVxAT1ZGZk/pz7fhUUmJ1eVrxX1MCPRwYC9zmnFtqZn+mYXqlSZxzjwKPAmRkZLgTqlLER5VVNby4bBdPfLKd/PIqBiR05XdXjuDS0ck6yCntSlMCPRfIdc4t9fZfoSHQC80syRudJwFFrVWkiB/yyg7y1Kc7eGFpDpXVtUxI78H9l4/gzEHxuqaKtEuNBrpzrsDMdpnZYOfcJmAKsMH7mg7M9m7ntWqlIm1kQ14Ff//HNt5ck4cDLhqRxMyvpzMiJdrv0kSOq6mH4W8DnvNWuGwDbgQ6AHPNbAaQA0xtnRJF2saKnXt5ZMkW3t9UTFREGNdPSOPGiWn6+DYJGE0KdOfcaiDjKHdNadFqRNqYc47Ptpby8JJsPttWSmxUR+48dxDXT0gjOqqj3+WJNIsWykpIcs7x/qYiHlmSzcqcMhK6RfKzi4cybVwqXSL130ICk35zJaTU1zveySzgkfezycyrIDmmM/ddNpypp6ZoxYoEPAW6hATnHO9kFvLHRZvZVFhJv55d+P1VI7lsTDIdwzr4XZ5Ii1CgS1D7cmrlwUWbWb+7gvSeXfjzNaO5ZGRvwrT0UIKMAl2CknOOT7JL+cOiTazKKaNPXGcemDqKy0b3JlwjcglSCnQJOku3lfKHRZtZtn0PSdGduP/yEUzNSNHUigQ9BboEjaz8CmYv3MiHm4uJ7xbJL78xjGvGpepgp4QMBboEvN1lB3nw3c28tiqXbpHh3H3hEK6fkEbnCAW5hBYFugSs8gM1/OXDbJ78ZAcA3/l6Ot+b3J+YqAh/CxPxiQJdAk5VTR3PfLaTR97PpqKqhsvHJPOjcwfpErYS8hToEjCccyxYm8/shRvZXXaQSYPimXXBEIb17u53aSLtggJdAsK63HJ+vSCTL3bsZVhSd3535UjOGNjT77JE2hUFurRrxZXVPPDOJuau2EVcVAS/vWIEV2f00UlBIkehQJd2qbq2jqc+2cHDS7Kpqqnj5jP6cduUgXTvpCsgihyLAl3ancVZhdy3YAM7Sg9w9pAEfnrxUPrHd/W7LJF2T4Eu7cauPQf41ZsbeC+rkP7xXXjqxq8xefBRP3tcRI5CgS6+O1Rbz2Mfb+OhxVswjFkXDmHGGf10qr5IMynQxVefbS3l5/PWk120j/NPSeTeb5xCckxnv8sSCUgKdPFFcWU19/9fFq+v2k1KbGcen57BlKGJfpclEtAU6NKmnHO89MUu7v+/LA7W1PGDswbw/bMG6LorIi2gSYFuZjuASqAOqHXOZZhZHPASkAbsAK52zu1tnTIlGOws3c+sV9fx2bZSTusXx39fPoIBCVq9ItJSmjNCP8s5V3LY/ixgsXNutpnN8vbvatHqJCjU1tXz5Cc7+MOiTXTs0IHfXjGCb2X0oYNODhJpUScz5XIpMNnbngN8gAJdjrCxoIK7XlnLmtxyzhmayG8uG06v6E5+lyUSlJoa6A5418wc8Dfn3KNAonMuH8A5l29mR10wbGYzgZkAqampLVCyBILq2jr+5/2t/OX9bKI7d+SRb4/h4hFJmGlULtJamhroE51zeV5oLzKzjU19AS/8HwXIyMhwJ1CjBJjMvHLunLuGjQWVXDEmmZ9fMozYLrpGuUhra1KgO+fyvNsiM3sdGAcUmlmSNzpPAopasU4JALV19fzto2386b3NxERF8MQNGZw9REsRRdpKo4FuZl2ADs65Sm/7PODXwHxgOjDbu53XmoVK+7ateB93vryGVTllXDIyifsuHa5RuUgba8oIPRF43Zv7DAeed869bWZfAHPNbAaQA0xtvTKlvaqvdzzz+U5+uzCLyPAwHpo2hm+O6u13WSIhqdFAd85tA0Ydpb0UmNIaRUlgyC8/yI9fXsvH2SVMHhzP764cSWJ3rWAR8YvOFJUT8vb6Au56dS01dfXcf/kIpo3roxUsIj5ToEuzHDxUx31vbeD5pTmMSI7moWlj6Nezi99liQgKdGmGjQUV3Pb8KrYU7eOWSenced5gIsJ1iVuR9kKBLo1yruHA52/eyqJ7p448fdM4Jg2K97ssETmCAl2Oa+/+Q/z4lbW8l1XIWYPj+f3UUfTsGul3WSJyFAp0OabVu8r4/nMrKaqs4t5LhnHjxDQd+BRpxxTo8hXOOZ7+bCe/eWsDCd068cqtpzOqT4zfZYlIIxTo8m/2Vddy92vreHNNHmcPSeDBq0cRE6UzPkUCgQJd/mlzYSW3PruCHSX7+ckFg7l1Un9ds1wkgCjQBYDXV+Vyz2vr6RIZznM3j2dC/x5+lyQizaRAD3G1dfXc/38beeKT7YzrF8cj08aQoNP3RQKSAj2E7d1/iB+8sJJPsku5cWIa91w0lI5hOlFIJFAp0EPUxoIKvvP0cgrLq/n9VSOZmtHH75JE5CQp0EPQwnX53PnyGrp1CuelW8YzJjXW75JEpAUo0ENIfb3jj+9t5uEl2YxJjeFv156q+XKRIKJADxEHD9Vx+0ureCezkG9l9OHXl51CZHiY32WJSAtSoIeAoooqbn56Oet2l/PzS4Zxk07hFwlKCvQgl5VfwYynvqDsYA1/vy6Dc4bpQ5tFgpUCPYi9v7GIHzy/km6dOjL3lgkMT472uyQRaUVNXnRsZmFmtsrMFnj7cWa2yMy2eLdaKtGOPP3ZDmbM+YK0nl144/sTFeYiIaA5Z5H8EMg6bH8WsNg5NxBY7O2Lz+rrHfct2MC98zI5e0gCc2+ZQK9orWQRCQVNCnQzSwEuBh47rPlSYI63PQe4rEUrk2arrq3jhy+t5vGPt3PD6Wn87boMukRqVk0kVDT1f/ufgJ8A3Q5rS3TO5QM45/LNLKGFa5NmqKyq4dZnV/BJdimzLhzCLZPStZJFJMQ0OkI3s0uAIufcihN5ATObaWbLzWx5cXHxiXwLaURxZTXXPPo5n2/bwx+mjuLWM/srzEVCUFNG6BOBb5rZRUAnoLuZPQsUmlmSNzpPAoqO9mTn3KPAowAZGRmuheoWz46S/Vz/xDKKK6t5bHoGZw3WGyWRUNXoCN05d7dzLsU5lwZcAyxxzl0LzAemew+bDsxrtSrlqNbllnPlXz9lX3UtL8wcrzAXCXEnc8RsNjDXzGYAOcDUlilJmmLptlJmzFlOTFRHnr5pHOnxXf0uSUR81qxAd859AHzgbZcCU1q+JGnMh5uLueWZ5STHdOa5m8drWaKIADpTNOC8k1nAbc+vYkBCV56eMY6eXSP9LklE2gkFegB5Y9Vu7nx5DSNTonnqhnFER3X0uyQRaUcU6AHi+aU5/PSNdZzWL47Hpn+NrjphSESOoFQIAI9/vJ37FmzgrMHx/PXaU+nUUdcxF5GvUqC3c4/9Yxu/eSuLC07pxUPTxhARrg9xFpGjU6C3Y49/vJ3fvJXFhcMbwrxjmMJcRI5NCdFOPflJwzTLlyNzhbmINEYp0Q499cl2fvXmBs4/JZGHv60wF5GmUVK0M3M+3cEv39zAucMSeXjaWIW5iDSZ0qIdefbznfxifibnDkvkf749VgdARaRZlBjtxOurcvnZG+s5e0iCwlxETohSox14J7OA/3p5LRPSe/CX/1CYi8iJUXL47OMtJdz2/CpGJEfz9+kZOmlIRE6YAt1HK3bu4TtPLyc9vgtP3ajT+UXk5CjQfZKZV84NT35BYvdInp4xjpioCL9LEpEAp0D3wfaS/Vz/+DK6RYbz7M2nkdBN1zMXkZOnQG9jxZXVTH9iGQ545ubTSImN8rskEQkSCvQ2tL+6lpue+oLiymoen55Bf31snIi0IB2FayM1dfV897mVbMiv4O/Xn8qY1Fi/SxKRIKMRehtwznHXq2v5aHMx/33ZcM4ekuh3SSIShBoNdDPrZGbLzGyNmWWa2a+89jgzW2RmW7xbDTmP4YF3N/Hayt3ccc4grhmX6nc5IhKkmjJCrwbOds6NAkYDF5jZeGAWsNg5NxBY7O3LEZ75fCf/8/5Wpo3rw39OGeB3OSISxBoNdNdgn7fb0ftywKXAHK99DnBZaxQYyD7cXMwv52cyZUgC9106HDPzuyQRCWJNmkM3szAzWw0UAYucc0uBROdcPoB3m9BqVQagzYWV/OC5lQxK7MZD08YQrsvgikgra1LKOOfqnHOjgRRgnJkNb+oLmNlMM1tuZsuLi4tPsMzAUrKvmpue+oJOEWE8Pj2DLjqlX0TaQLOGjc65MuAD4AKg0MySALzbomM851HnXIZzLiM+Pv7kqg0AVTV13PLMCkr2VfPY9Rn0junsd0kiEiKassol3sxivO3OwDnARmA+MN172HRgXivVGDC+XJ64YudeHrx6NKP6xPhdkoiEkKbMBSQBc8wsjIY/AHOdcwvM7DNgrpnNAHKAqa1YZ0B4eEk281bn8ePzB3PRiCS/yxGRENNooDvn1gJjjtJeCkxpjaIC0TuZBTy4aDNXjE3me5P7+12OiIQgLb1oAdlFlfzopdWM6hPD/ZeP0PJEEfGFAv0kVVTVMPPpFXSOCON/rx2rTxwSEd9oPd1JqK933PHianL2HOD574wnKVorWkTEPxqhn4Q/vbeZxRuLuPcbwxjXL87vckQkxCnQT9Db6wt4aEk2U09N4brxff0uR0REgX4ithbv4865DQdB77tM12gRkfZBgd5MBw/V8f3nVhIR3oG//ocOgopI+6GDos30i/nr2VRYyZM3fE2n9YtIu6IRejO8vHwXc5fn8oOzBjB5sC4uKSLtiwK9iTYVVPLzeeuZkN6D288Z5Hc5IiJfoUBvgv3VtXz3uRV0jezIn6eNJqyDDoKKSPujQG+Ec457Xl/HjpL9PDRtNAndOvldkojIUSnQG/Hyilzmrc7jjnMGcXr/nn6XIyJyTAr049hesp9fzs/ktH5xfO8sfcCziLRvCvRjqKmr5/YXV9ExrAN//JbmzUWk/dM69GP446LNrMkt56//MVbrzUUkIGiEfhSfbS3lrx9u5VsZfbhQnzwkIgFCgX6EsgOH+NHc1aT16MK93xjmdzkiIk2mKZfDOOf46evrKa6s5rXvnU6XSP14RCRwaIR+mPlr8nhrXT4/Om8QI1Ni/C5HRKRZGg10M+tjZu+bWZaZZZrZD732ODNbZGZbvNvY1i+39RRWVHHvvEzGpsZwyyR9yLOIBJ6mjNBrgTudc0OB8cD3zWwYMAtY7JwbCCz29gOSc467X1tHdW0dD0wdpSWKIhKQGg1051y+c26lt10JZAHJwKXAHO9hc4DLWqnGVvfy8lyWbCziJ+cPIT2+q9/liIickGbNoZtZGjAGWAokOufyoSH0gYC8nuzusoP8esEGTusXxw2np/ldjojICWtyoJtZV+BV4HbnXEUznjfTzJab2fLi4uITqbHVOOe465W11DvH768aRQdNtYhIAGtSoJtZRxrC/Dnn3Gtec6GZJXn3JwFFR3uuc+5R51yGcy4jPj6+JWpuMc8tzeHj7BJ+evFQUntE+V2OiMhJacoqFwMeB7Kccw8edtd8YLq3PR2Y1/LltZ788oPMXriRiQN68O1xqX6XIyJy0ppy5sxE4DpgnZmt9truAWYDc81sBpADTG2VCluBc46fv5FJbX09v718JA1/s0REAlujge6c+xg4VuJNadly2sbC9QW8l1XI3RcO0VSLiASNkDtTtPxADb+Yn8kpvbsz44x+fpcjItJiQu5iJbPfzmLP/kM8ecPXCA8Lub9nIhLEQirRPt9WygvLdjHjjH4MT472uxwRkRYVMoFeXVvHPa+to09cZ+44Z5Df5YiItLiQmXL5+0fb2Faynzk3jaNzRJjf5YiItLiQGKHv2nOAR97P5sLhvThzUPs6uUlEpKWERKDft2ADhvGzS/QJRCISvII+0N/fVMS7Gwq5bcoAkvVhzyISxII60Ktq6vjl/EzS47tw8xnpfpcjItKqgvqg6N8/2sbO0gM8M2McEeFB/bdLRCR4R+i5exsOhF40ohdfH6gDoSIS/II20Gcv3IgZ/OxiHQgVkdAQlIG+YuceFqzNZ+ak/vTWgVARCRFBF+j19Y5fL8gioVskt56pA6EiEjqCLtDfXJvHml1l/Pj8wURFBPUxXxGRfxNUgV5VU8fvFm5keHJ3rhyb4nc5IiJtKqgC/bF/bCOvvIqfXTxMH/gsIiEnaAK9qLKKv3ywlfNPSWR8eg+/yxERaXNBE+gPL87mUG09sy4c6ncpIiK+CIpAzyk9wAvLcvjW1/rQr2cXv8sREfFFo4FuZk+YWZGZrT+sLc7MFpnZFu82tnXLPL4/vreZ8DDjP6cM9LMMERFfNWWE/hRwwRFts4DFzrmBwGJv3xcbCyp4Y/Vubji9H4ndO/lVhoiI7xoNdOfcR8CeI5ovBeZ423OAy1q2rKZ74J1NdI0M57tn9verBBGRduFE59ATnXP5AN5twrEeaGYzzWy5mS0vLi4+wZc7uhU79/BeVhG3ntmf6KiOLfq9RUQCTasfFHXOPeqcy3DOZcTHt+xVDx94ZzM9u0Zy48S0Fv2+IiKB6EQDvdDMkgC826KWK6lpvtixh8+2lfLdyf11ir+ICCce6POB6d72dGBey5TTdA8t3kLPrhF8e1xqW7+0iEi71JRliy8AnwGDzSzXzGYAs4FzzWwLcK6332ZW5ezlH1tKuPnr6XSOCGvLlxYRabcanatwzk07xl1TWriWJnt4STaxUR25bnxfv0oQEWl3Au5M0cy8cpZsLGLGGf3oEqm5cxGRLwVcoD/+8XaiIsK4bkKa36WIiLQrARXoRRVVvLkmj6sz+hDdWevORUQOF1CB/sznO6mtd9xweprfpYiItDsBE+hVNXU8tzSHKUMSSdMVFUVEviJgAv29rEL27D/E9RO0skVE5GgCJtDnLs8lOaYzEwf09LsUEZF2KSACPa/sIP/YUsyVY5MJ02eFiogcVUAE+msrc3EOrjq1j9+liIi0WwER6AndO/GtjD6k9ojyuxQRkXYrIE61vDqjD1dnaHQuInI8ATFCFxGRxinQRUSChAJdRCRIKNBFRIKEAl1EJEgo0EVEgoQCXUQkSCjQRUSChDnn2u7FzIqBnSf49J5ASQuWEwjU59CgPoeGk+lzX+dcfGMPatNAPxlmttw5l+F3HW1JfQ4N6nNoaIs+a8pFRCRIKNBFRIJEIAX6o34X4AP1OTSoz6Gh1fscMHPoIiJyfIE0QhcRkeNQoIuIBImACHQzu8DMNplZtpnN8rue5jCzPmb2vpllmVmmmf3Qa48zs0VmtsW7jT3sOXd7fd1kZucf1n6qma3z7nvIzMxrjzSzl7z2pWaW1uYdPYKZhZnZKjNb4O0HdX8BzCzGzF4xs43ev/eEYO63md3h/U6vN7MXzKxTMPbXzJ4wsyIzW39YW5v008yme6+xxcymN1qsc65dfwFhwFYgHYgA1gDD/K6rGfUnAWO97W7AZmAY8P+AWV77LOB33vYwr4+RQD+v72HefcuACYABC4ELvfbvAf/rbV8DvNQO+v0j4Hlggbcf1P31apkD3OxtRwAxwdpvIBnYDnT29ucCNwRjf4FJwFhg/WFtrd5PIA7Y5t3Getuxx63V7/8ETfhhTgDeOWz/buBuv+s6if7MA84FNgFJXlsSsOlo/QPe8X4GScDGw9qnAX87/DHedjgNZ6OZj31MARYDZ/OvQA/a/np1dKch4OyI9qDsNw2BvssLm3BgAXBeEPc3jX8P9Fbv5+GP8e77GzDteHUGwpTLl784X8r12gKO91ZqDLAUSHTO5QN4twnew47V32Rv+8j2f3uOc64WKAd6tEonmuZPwE+A+sPagrm/0PAOshh40ptqeszMuhCk/XbO7QYeAHKAfKDcOfcuQdrfo2iLfjY7+wIh0O0obQG31tLMugKvArc75yqO99CjtLnjtB/vOW3OzC4BipxzK5r6lKO0BUx/DxNOw9vyvzrnxgD7aXgrfiwB3W9vzvhSGqYVegNdzOza4z3lKG0B099maMl+Nrv/gRDouUCfw/ZTgDyfajkhZtaRhjB/zjn3mtdcaGZJ3v1JQJHXfqz+5nrbR7b/23PMLByIBva0fE+aZCLwTTPbAbwInG1mzxK8/f1SLpDrnFvq7b9CQ8AHa7/PAbY754qdczXAa8DpBG9/j9QW/Wx29gVCoH8BDDSzfmYWQcNBg/k+19Rk3pHsx4Es59yDh901H/jyqPV0GubWv2y/xjvy3Q8YCCzz3tZVmtl473tef8RzvvxeVwFLnDfp1tacc3c751Kcc2k0/Fstcc5dS5D290vOuQJgl5kN9pqmABsI3n7nAOPNLMqrcwqQRfD290ht0c93gPPMLNZ7R3Se13ZsfhxgOIEDEhfRsDpkK/BTv+tpZu1n0PA2aS2w2vu6iIY5ssXAFu827rDn/NTr6ya8I+Feewaw3rvvEf51pm8n4GUgm4Yj6el+99urazL/OigaCv0dDSz3/q3foGFlQtD2G/gVsNGr9RkaVnYEXX+BF2g4TlBDw6h5Rlv1E7jJa88GbmysVp36LyISJAJhykVERJpAgS4iEiQU6CIiQUKBLiISJBToIiJBQoEuIhIkFOgiIkHi/wO4VYwMW4kEFgAAAABJRU5ErkJggg==", "text/plain": [ "
" ] @@ -83,47 +112,47 @@ "source": [ "V = Ausgleichsbecken_class(1.,0.1,0.,10.,timestep=0.001)\n", "\n", - "V.set_initial_level(5.) \n", - "V.set_influx(0.75)\n", - "V.set_outflux(0.)\n", + "V.set_initial_level(initial_level) \n", + "V.set_influx(initial_influx)\n", + "V.set_outflux(initial_outflux)\n", "\n", - "V.p0 = 0.\n", + "V.p0 = initial_pipeline_pressure\n", "\n", "outflux_vec = []\n", "level_vec = []\n", "\n", "t = 0\n", - "while V.level > 0.01:\n", + "while V.level > total_min_level:\n", " t = t+V.timestep\n", " V.e_RK_4()\n", " V.level = V.update_level(V.timestep)\n", " V.set_volume()\n", " outflux_vec.append(V.outflux)\n", " level_vec.append(V.level)\n", - " if t > 100:\n", + " if t > total_max_time:\n", " break\n", - " \n", + "\n", "plt.plot(level_vec)\n" ] }, { "cell_type": "code", - "execution_count": 18, + "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "text/plain": [ - "[]" + "[]" ] }, - "execution_count": 18, + "execution_count": 12, "metadata": {}, "output_type": "execute_result" }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAD4CAYAAAD1jb0+AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAgSklEQVR4nO3deXiU9b3+8fdnZrKQFUJCCEvYV9kJiKJYRetet6qoIG5FLVrrsadHbc9pz/m1tdXWrcUqRUTEtSqute77giTsEPZ9CwkJCSGQkOT7+yPTliprZpJnJnO/rmuuDE8m89z5XuPtk2f7mnMOERGJPj6vA4iISOOowEVEopQKXEQkSqnARUSilApcRCRKBZpzZZmZma5r167NuUoRkahXUFBQ4pzL+ubyZi3wrl27kp+f35yrFBGJema24WDLtQtFRCRKqcBFRKKUClxEJEqpwEVEopQKXEQkSqnARUSilApcRCRKNet54LGguraOnZU1lO6pYU91LXtqaqmsrqOqupY9NXXU1tVT76DeOerrHfUOzCAh4Gt4xPlJjPOREPCTlhhH66SGR5ukeJLi/ZiZ17+iiEQIFfgxqKt3bCvfy8bSKjaVVgW/7mV7xT5KKqsp2V1Nxb7aJlt/vN9HRnI8Oa0T6ZDeig6tE8lJb0WH1q3olplMl7ZJJMb5m2z9IhJZVOAHUV/v2Fy2l+XbK1hZtJsVRZWs3L6btSWV7K/71wQYfp81lGhaK/q2TyWzZyZZKQlkpibQJime1MQAyQkBkuP9wa8B4gKGzwwz8FvD83rnqKmrp3p/PdW19VTX1rF3fx0Ve2spq6phV1UNZVX7KauqoWR3DdvK97JsWwXvFRZRXVv/zzxm0KlNK3pkpdA9M4Xe2SkM6JhOr+wUEgIqdpGWRgUO7Ni9jwUbd7FgU8Nj0eZyKqv/tSXdOaMVfbJTObVvO7q2TSI3I4nOGUnkpCcS8Id+GMGHEfD7SIo/tp9zzlFWtZ+tu/aytmQPa3ZUsrZkD2uLK5mztpS9++sAiPMbvbNTOa5DGgM6pjO8Sxv6tk/D79PuGJFoFpMFvrmsii/W7OSrNTuZs66ULbv2AhDwGX1zUrlwaAcGdEinT/tUemWnkpIQmcNkZmQkx5ORHM+Ajun/9r36esemsiqWbKlgydZylm6t4L3CHbyQvxmAlIQAw7q0YUSXNuR1zWBobmvtfhGJMtacc2Lm5eU5L25mVV61n49XFfP5qhK+WFvCptKGwm6bHM/x3TMYltuGIZ1bM6BjeosuMeccW3btpWBDGXPXl5K/vowVRbtxDuIDPkZ2zWBM70zG9M6iT3aqDpiKRAgzK3DO5X1reUsscOcca0v28EHhDt5fXsTc9WXU1TvSW8VxfLcMTuzRlhN6ZNI7OyXmS6q8aj8FG0v5YvVOPllVzMqiSgDapSZwcq8sxvZrxym9s0iO0L9CRGJBiy9w5xxLt1bwxqJt/H3JNtbvrAKgb/tUxvZrx2l9sxnSubX2+x7BtvK9fLqqhE9WFvPZ6hJ2Ve0nIeDj5F5ZnHlcNqf3y6ZN8jHurBeRkLTYAl9ZtJvXF27ljUXbWFeyh4DPOLFnJmf0a8dp/bLp2LpVWNcXS2rr6pm7voy3l27nnaXb2Vq+D7/POL5bBhcO6chZA9uTlhjndUyRFq/RBW5m04HzgB3OuQHBZRnA80BXYD1wmXOu7EghwlXgJZXVvDJ/C3/N38yKot34DEZ1b8v5gztw1nHttYXYBJxzLN5SzttLt/Pmooa/cOIDPk7v144Lh3TkO33aER/Qhb0iTSGUAh8DVAIzDyjwe4FS59xvzexOoI1z7r+OFCKUAq+tq+fjlcW8kL+J9wt3UFvvGJrbmouGduTsATlkpSY06n3l2DnnWLi5nFfmb+H1hVvZuaeG9FZxnD84hytG5nJch/Qjv4mIHLWQdqGYWVfgjQMKfAXwHefcNjPLAT5yzvU50vs0tsBnfL6OKR+toXh3NZkp8Vw8rBOXDu9Er+zUY34vCa/9dfV8trqEV+Zv4e9LtlNdW8/gzq25cmRnzh/cgaR4HfwUCVW4C3yXc671Ad8vc861OcTPTgImAeTm5g7fsOGgU7sd1uOfrePLNTu5LK8Tp/ZtR1wYLp6R8NtVVcPs+Vt4Zs5GVu2oJCUhwIVDOzBhVFf6tNf/bEUay7MCP5BX54FL83LOUbChjGfmbOSNxduoqa3npJ6ZXH9SN07pnYVPZwKJHJNDFXhj/74tMrOcA3ah7AgtnrQkZkZe1wzyumbw3+f155mvNzLzy/VcO2Mu3bOSuW50Ny4Z1olW8S33oimR5tDYfRGvARODzycCr4YnjrQ0bZLjmXxqTz796Wk8ePkQkuMD/PyVJZzw2/d54N2V7Kqq8TqiSNQ6mrNQngW+A2QCRcAvgFeAF4BcYCNwqXOu9Egr0y4Ucc6Rv6GMqZ+s5d1lRSTH+xl/QhduOKm7ziQSOYQWeyGPRK/l2yuY8uEa3ly0lTi/jytG5jJpTHc66OIrkX+jApeItba4kj9/tIbZ87dgBleOzGXyaT1pl5rodTSRiKACl4i3uayKKR+u4YX8TcT7fVw7uis3julBepIu15fYpgKXqLG+ZA8PvLeS1xZuJSUhwE2n9OCaE7vqjogSs1TgEnUKt1Xwh3dW8l5hEZkp8fzHGX24fERn3VFSYs6hClyXNErE6peTxrSJebz8wxPplpnM3bMXc+7Dn/LZqhKvo4lEBBW4RLxhuW144cYTeOSqYeypqWX843O4fsZc1hRXeh1NxFMqcIkKZsY5A3N49/ZTuPPsvsxZV8qZD3zCL19bSvne/V7HE/GEClyiSmKcn5tO6cFH//kdLhvRmZlfrmfsHz5i9vzNNOfxHJFIoAKXqJSZksBvLhrIa7ecRMc2Sdz+/ELGTf2KVUW7vY4m0mxU4BLVBnRMZ/bNJ3LPxQNZvn03Zz/0Kfe8Vcie6lqvo4k0ORW4RD2fz7hiZC4f3HEKFw/ryGMfr+WM+z/mg+VFXkcTaVIqcGkx2qYkcO/3B/PiTSeQkhjguhn5/Pi5+ZTu0R0PpWVSgUuLk9c1gzduPZnbxvbizcXbOP3+j3lt4VYd5JQWRwUuLVJ8wMftZ/Tm9VtPonObVvzo2fn8YGYBRRX7vI4mEjYqcGnR+rZP46WbT+Tuc/ry6apiTr//Y/6av0lb49IiqMClxQv4fUwa04O3fzyGfu3T+M8XF3HTrAJ2VlZ7HU0kJCpwiRldM5N5dtIo7jq7Lx8uL+bMBz/l/UKdqSLRSwUuMcXvM248pQev3jKazJR4rn8yn7teXqzzxiUqqcAlJvXLSePVW0Zz4yndeW7uRs55+FMKNpR5HUvkmKjAJWYlBPzcdXY/np90AnX1jsse+5IpH66mvl4HOCU6qMAl5o3slsFbt53MOQNzuO/tFVw9/Wt27NbphhL5VOAiQGpiHA+PG8LvLhlI/oZSznnoUz5ZWex1LJHDUoGLBJkZl4/I5fVbTqJtcgJXT/+a3761nP119V5HEzkoFbjIN/TKTuXVW0Zz5fG5PPrxGi577Eu27trrdSyRb1GBixxEYpyf31w0kClXDmNVUSXn/fEzzcUpEUcFLnIY5w7K4bXgOeNXT5+js1QkoqjARY6ge1YKr0wezXmDOnDf2yuY9FSB5uGUiKACFzkKSfEBHho3hF+e35+PVuzge3/6jMJtFV7HkhgXUoGb2e1mttTMlpjZs2aWGK5gIpHGzLhmdDeev3EU+/bXcdEjn/PK/C1ex5IY1ugCN7OOwI+APOfcAMAPjAtXMJFINbxLw4QRgzu15sfPL+Cetwqp035x8UCou1ACQCszCwBJwNbQI4lEvqzUBGbdcDxXHZ/LYx+v5Qcz89m9T/vFpXk1usCdc1uA3wMbgW1AuXPunW++zswmmVm+meUXF+vKNmk54vw+fn3RQP7fhQP4ZGUxFz3yBetL9ngdS2JIKLtQ2gAXAN2ADkCymY3/5uucc1Odc3nOubysrKzGJxWJUBNGdWHm9SMpqazmgimf63xxaTah7EI5HVjnnCt2zu0HXgZODE8skehyYo9MXpt8Eu3TEpn4xNc88fk6TdsmTS6UAt8IjDKzJDMzYCxQGJ5YItEnt20SL/3wRE7t047/fX0Z//3qEmp1HxVpQqHsA58DvAjMAxYH32tqmHKJRKWUhABTJwxn0pjuzPpqIz+YmU+lZvuRJmLN+WdeXl6ey8/Pb7b1iXhp1lcb+MVrS+mTncr0a0bQPl2XSUjjmFmBcy7vm8t1JaZIExk/qgvTJuaxYeceLpzyOcu26spNCS8VuEgTOrVPO/56U8Ox/Usf/YKPVuzwOJG0JCpwkSbWv0MasyefSG7bZK5/Mp+n52zwOpK0ECpwkWaQk96Kv950Aif3yuRns5dw/zsrdJqhhEwFLtJMUhICTLs6j8vyOvHwB6u5e/YS3UNFQhLwOoBILAn4ffzukkFkpiTwyEdrKN1TzUPjhpIY5/c6mkQhbYGLNDMz46dn9eV/zuvP20uLuHr615ogQhpFBS7iketO6sbDVwxl/sYyLn/sS4oq9nkdSaKMClzEQ98b3IHp14xgU2kVl/z5C9YWV3odSaKIClzEYyf3yuLZSaPYW1PHpY9+ydKt5V5HkiihAheJAIM6tebFm08kIeDjiqlfUbChzOtIEgVU4CIRoltmMn+9+UQykuOZ8Pgcvlit+4rL4anARSJIx9ateOHGE+jcJolrZszlg+VFXkeSCKYCF4kw7dISeW7SKPq2T2XSzALeWKSpZuXgVOAiEahNcjxP33A8w3Lb8KNn5/NC/iavI0kEUoGLRKjUxDievG4ko3tm8tMXFzHj83VeR5IIowIXiWCt4v1Mm5jHmcdl88vXlzHt07VeR5IIogIXiXAJAT9/unIY5w7M4VdvFjL1kzVeR5IIoZtZiUSBOL+Ph8YNwQx+87fl1Du46ZQeXscSj6nARaJEwO/jwcuHYGb89q3l1DvHD7/T0+tY4iEVuEgUCfh9PHDZYHwG9/59Bc7B5FNV4rFKBS4SZQJ+H/dfNgSfGfe9vYL6esetY3t5HUs8oAIXiUJ+n/H7SwdjwB/eXUm9g9tOV4nHGhW4SJTy+4z7Lh0MBg+8txKfoS3xGKMCF4lifp9x3/cH41zDlnhinJ8fjOnudSxpJipwkSjXUOKDqKmt59d/KyQhzsfVJ3T1OpY0AxW4SAsQ8Pt4cNwQqmvr+Z9XlxLv9zFuZK7XsaSJ6UpMkRYizu9jylVDOaV3FnfNXszs+Zu9jiRNTAUu0oIkBPw8NmE4o7q15Y4XFvLmom1eR5ImFFKBm1lrM3vRzJabWaGZnRCuYCLSOIlxfh6/Jo9huW247bn5vLtMk0K0VKFugT8E/N051xcYDBSGHklEQpUUH+CJa0dwXMd0Jj89j49W7PA6kjSBRhe4maUBY4DHAZxzNc65XWHKJSIhSk2MY+a1I+nZLoWbZhUwd32p15EkzELZAu8OFANPmNl8M5tmZsnffJGZTTKzfDPLLy4uDmF1InKs0pPimHn9SDqkt+K6GXNZtrXC60gSRqEUeAAYBvzZOTcU2APc+c0XOeemOufynHN5WVlZIaxORBojMyWBp244npSEAFdP/5p1JXu8jiRhEkqBbwY2O+fmBP/9Ig2FLiIRpmPrVjx1/fHUO8f4aXPYXr7P60gSBo0ucOfcdmCTmfUJLhoLLAtLKhEJu57tUnjy2pGU793P+MfnULqnxutIEqJQz0K5FXjazBYBQ4DfhJxIRJrMwE7pTJuYx8bSKq594msqq2u9jiQhCKnAnXMLgvu3BznnLnTOlYUrmIg0jVHd2zLlymEs2VrBpJn57Ntf53UkaSRdiSkSg87on8193x/EF2t28qNn51NbV+91JGkEFbhIjLp4WCd+cX5/3llWxN2zF+Oc8zqSHCPdjVAkhl07uhtle2p4+IPVZKclcsd3+xz5hyRiqMBFYtztZ/SmqKKaPwZLfPyoLl5HkqOkAheJcWbGry8aQHFlNf/z6hKyUhM487j2XseSo6B94CJCwO/jT1cOZVCn1vzo2fm6b0qUUIGLCNBwB8Pp14ygY+tWXD9jLquKdnsdSY5ABS4i/5SRHM+T140kIc7PxOlfs618r9eR5DBU4CLybzpnJPHENSOo2FfLNdPnUr53v9eR5BBU4CLyLQM6pvPo+OGsLanU1ZoRTAUuIgd1Uq9Mfn/pYOasK+WOvy6kvl4X+kQanUYoIod0wZCObC/fxz1vLadzmyTuPLuv15HkACpwETmsSWO6s6G0ikc/XkOXtklcMTLX60gSpAIXkcMyM/7ve8expWwvP39lCR1at+KU3ppdKxJoH7iIHFHA72PKVcPonZ3K5KfnUbhNc2tGAhW4iByVlIQA06/JIznBz3Uz5lJUoWnZvKYCF5GjlpPeiunXjKBi736umzGXPZrRx1MqcBE5Jsd1SOdPVw6jcFsFt2oyCE+pwEXkmJ3atx3/e8EAPli+g/97Y5kmg/CIzkIRkUaZMKoLm0qrmPrJWnIzkrjh5O5eR4o5KnARabQ7z+rLptIqfv23QnIzkviu7iPerLQLRUQazeczHrh8CIM6pvPj5xewbKtOL2xOKnARCUlinJ+/XJ1HWmIcNzw5lx27dXphc1GBi0jI2qUlMm1iHmVV+7nxqQLdvbCZqMBFJCwGdEzngcsHM3/jLv7rpUU6M6UZqMBFJGzOGpDDf57Zh1cXbGXKh6u9jtPi6SwUEQmrH36nB6uKdvP7d1bSIyuFswfmeB2pxdIWuIiElZnx20sGMTS3Nbe/sIAlW8q9jtRihVzgZuY3s/lm9kY4AolI9EuM8zN1Qh4ZSfHc8GS+bnzVRMKxBX4bUBiG9xGRFiQrNYFpE0dQsW8/P5iZz94anZkSbiEVuJl1As4FpoUnjoi0JP07pPHQuKEs3lLOT15cqDNTwizULfAHgZ8Ch7wdmZlNMrN8M8svLi4OcXUiEm3O6J/Nf53VlzcXbePh93VmSjg1usDN7Dxgh3Ou4HCvc85Ndc7lOefysrI0DZNILLpxTHcuHtaRB95bydtLt3sdp8UIZQt8NPA9M1sPPAecZmazwpJKRFoUM+M3Fw1kcKd0/uP5BazYvtvrSC1CowvcOXeXc66Tc64rMA74wDk3PmzJRKRFSYzz89iEPJISAvxgZj67qmq8jhT1dB64iDSb9umJPDp+ONvL93HLM5rNJ1RhKXDn3EfOufPC8V4i0rIN79KGX100gM9Wl3DPW8u9jhPVdCm9iDS7y/I6s2xrBY9/to5+OWl8f3gnryNFJe1CERFP/OzcfpzYoy13z17M/I1lXseJSipwEfFEnN/HlCuHkZ2WwI1PFehy+0ZQgYuIZ9okx/OXq/OorK7VRBCNoAIXEU/1bZ/G/ZcNZsGmXfz8lSW63P4YqMBFxHNnDcjhtrG9eLFgMzO+WO91nKihAheRiHDb2F58t382v3qzkM9Xl3gdJyqowEUkIvh8xv2XD6FHVjK3PDOPTaVVXkeKeCpwEYkYKQkBpk7Io7becdMsHdQ8EhW4iESUrpnJPDRuCMu2VXD3y4t1UPMwVOAiEnFO65vN7af35uX5W3hSBzUPSQUuIhHpllN7cnq/hoOac9bu9DpORFKBi0hEajioOZjcjCQmPzOPbeV7vY4UcVTgIhKx0hLjmHr1cPbW1HHzrHlU1+qg5oFU4CIS0Xq2S+UPwSs1f/naUq/jRBQVuIhEvLMG5DD51B48+/Umnpmz0es4EUMFLiJR4T/O6MOY3ln84rUlzNPtZwEVuIhECb/PeHjcEHLSW3HzrAJ27NbtZ1XgIhI1WifF89iE4VTsrWXy0/OoqY3tOTVV4CISVfrlpPG77w9i7voyfv3mMq/jeEpzYopI1Pne4A4s2rSLaZ+tY1Cn1lwSo3NqagtcRKLSnWf3/eecmos3l3sdxxMqcBGJSgG/jz9eMZTMlARumlXAzspqryM1OxW4iESttikJPDp+OMWV1dz67Hxq62LroKYKXESi2sBO6fz6wgF8sWYn9769wus4zUoHMUUk6l2a15nFW8qZ+slaBnZM5/zBHbyO1Cy0BS4iLcLPz+1PXpc2/PTFRSzfXuF1nGahAheRFiE+4OOR8cNITQwwaWYB5VX7vY7U5Bpd4GbW2cw+NLNCM1tqZreFM5iIyLFql5rIn8cPZ1v5Xm57fj519S17OrZQtsBrgTucc/2AUcBkM+sfnlgiIo0zvEsbfnH+cXy0opgH31vpdZwm1egCd85tc87NCz7fDRQCHcMVTESksa46PpfL8zrzxw9W8/bS7V7HaTJh2QduZl2BocCccLyfiEgozIz/veA4BndK544XFrJ6R6XXkZpEyAVuZinAS8CPnXPfOvRrZpPMLN/M8ouLi0NdnYjIUUmM8/Pn8cNJCPiY9FQ+u/e1vIOaIRW4mcXRUN5PO+dePthrnHNTnXN5zrm8rKysUFYnInJMOrRuxZSrhrFhZxV3vLCQ+hZ2UDOUs1AMeBwodM7dH75IIiLhM6p7W+4+px/vLCvikY9Wex0nrELZAh8NTABOM7MFwcc5YcolIhI2143uyoVDOvCHd1fy4YodXscJm1DOQvnMOWfOuUHOuSHBx9/CGU5EJBzMjHsuHkS/9mnc9ux81pfs8TpSWOhKTBGJCa3i/Tw2YTg+n3HTrAKqamq9jhQyFbiIxIzOGUk8PG4oK4t289MXF+FcdB/UVIGLSEwZ0zuLn5zZhzcWbWPap+u8jhMSFbiIxJybT+nBOQPbc89bhXyxusTrOI2mAheRmGNm3Pv9wfTISmHyM/PYXFbldaRGUYGLSExKSQjw2ITh1NY5bppVwL79dV5HOmYqcBGJWd2zUnhw3BCWbKngZ7OXRN1BTRW4iMS0sf2yuW1sL16at5mnvtrgdZxjogIXkZh329hejO3bjv97fRlz15d6HeeoqcBFJOb5fMYD44bQOSOJm2fNY3v5Pq8jHRUVuIgIkJYYx2MThlNVUxs1BzVV4CIiQb2zU7n/siEs2LSLO1+K/Cs1VeAiIgc4a0B7fvLd3ryyYCuPfLTG6ziHFfA6gIhIpJl8ak9W7ajkvrdX0LNdCmce197rSAelLXARkW8wM353ySAGd27N7c8vYNnWb80WGRFU4CIiB5EY5+cvE4aT3iqOG56cS/Huaq8jfYsKXETkENqlJfKXq/MorarhxqfyI+7MFBW4iMhhDOiYzv2XDWHexl3c/fLiiDozRQUuInIE5wzM4fbTe/Py/C08/H7kTIyss1BERI7Cj8b2ZEPpHh54byUdWidyaV5nryOpwEVEjoaZ8duLB1G8u5q7Xl5Mu7RETumd5Wkm7UIRETlK8QEfj1w1jN7ZqfxwVgFLtpR7mkcFLiJyDFIT43ji2hG0Torn2hlz2VTq3Ww+KnARkWOUnZbIjGtHUL2/jolPfE1JpTfniKvARUQaoVd2Ko9fM4Ktu/Yy4fGv2VVV0+wZVOAiIo00omsGf7k6jzU7Kpk4/Wt279vfrOtXgYuIhODkXlk8ctUwlm6t4LoZc6mqqW22davARURCdHr/bB4aN5SCDWWMnzan2XanqMBFRMLg3EE5PHLVcJZsqeDyx76iqKLpp2ULqcDN7CwzW2Fmq83sznCFEhGJRmcNaM8T145gc1kVFz/yBQs37WrS9TW6wM3MD0wBzgb6A1eYWf9wBRMRiUaje2by3KQTALj00S+Z9ula9tfVN8m6QtkCHwmsds6tdc7VAM8BF4QnlohI9BrYKZ03bj2Jk3tl8qs3Cxlz74d8sbok7OsJpcA7ApsO+Pfm4LJ/Y2aTzCzfzPKLi4tDWJ2ISPRokxzPtIl5PHHNCI7rkEbnjKSwryOUm1nZQZZ960a5zrmpwFSAvLy8yLmRrohIEzMzTu3bjlP7tmuS9w9lC3wzcOD9FDsBW0OLIyIiRyuUAp8L9DKzbmYWD4wDXgtPLBEROZJG70JxztWa2S3A24AfmO6cWxq2ZCIiclghTejgnPsb8LcwZRERkWOgKzFFRKKUClxEJEqpwEVEopQKXEQkSplzzXdtjZkVAxsa+eOZQPivRW15NE5HR+N0ZBqjo9Mc49TFOZf1zYXNWuChMLN851ye1zkincbp6GicjkxjdHS8HCftQhERiVIqcBGRKBVNBT7V6wBRQuN0dDROR6YxOjqejVPU7AMXEZF/F01b4CIicgAVuIhIlIqKAtfkyf9iZuvNbLGZLTCz/OCyDDN718xWBb+2OeD1dwXHbYWZneld8qZlZtPNbIeZLTlg2TGPi5kND47vajN72MwONnFJ1DrEOP3SzLYEP1MLzOycA74Xc+NkZp3N7EMzKzSzpWZ2W3B55H2enHMR/aDhVrVrgO5APLAQ6O91Lg/HYz2Q+Y1l9wJ3Bp/fCfwu+Lx/cLwSgG7BcfR7/Ts00biMAYYBS0IZF+Br4AQaZpx6Czjb69+tGcbpl8BPDvLamBwnIAcYFnyeCqwMjkXEfZ6iYQtckycf2QXAk8HnTwIXHrD8OedctXNuHbCahvFscZxznwCl31h8TONiZjlAmnPuS9fwX9/MA36mRTjEOB1KTI6Tc26bc25e8PluoJCG+X4j7vMUDQV+VJMnxxAHvGNmBWY2Kbgs2zm3DRo+fMA/JuCL9bE71nHpGHz+zeWx4BYzWxTcxfKPXQMxP05m1hUYCswhAj9P0VDgRzV5cgwZ7ZwbBpwNTDazMYd5rcbu4A41LrE6Xn8GegBDgG3AH4LLY3qczCwFeAn4sXOu4nAvPciyZhmnaChwTZ58AOfc1uDXHcBsGnaJFAX/XCP4dUfw5bE+dsc6LpuDz7+5vEVzzhU55+qcc/XAX/jXbraYHSczi6OhvJ92zr0cXBxxn6doKHBNnhxkZslmlvqP58B3gSU0jMfE4MsmAq8Gn78GjDOzBDPrBvSi4aBKrDimcQn+WbzbzEYFzxa4+oCfabH+UUpBF9HwmYIYHafg7/Q4UOicu/+Ab0Xe58nrI75HeVT4HBqOBK8BfuZ1Hg/HoTsNR7sXAkv/MRZAW+B9YFXwa8YBP/Oz4LitoAWdKXCQsXmWhj//99Ow5XN9Y8YFyKOhwNYAfyJ4tXJLeRxinJ4CFgOLaCijnFgeJ+AkGnZ1LAIWBB/nROLnSZfSi4hEqWjYhSIiIgehAhcRiVIqcBGRKKUCFxGJUipwEZEopQIXEYlSKnARkSj1/wESJRTNgSVxTwAAAABJRU5ErkJggg==", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAhc0lEQVR4nO3deXxU1cHG8d+ZyUbIRhYCJIGwQ9ghbKLWDUVLcV8Q0doqKqhttVqX19a3y6uV1lq1ogi4Iu5L645WUVCBhH2VPewkJAFCQtbz/jFDRQqSwEzuncnz/Xzmk8nMJPMwJE/unHvuPcZai4iIuJfH6QAiIvLDVNQiIi6nohYRcTkVtYiIy6moRURcLiIY3zQ1NdVmZ2cH41uLiISl/Pz8Imtt2pHuC0pRZ2dnk5eXF4xvLSISlowxm452n4Y+RERcTkUtIuJyKmoREZdTUYuIuJyKWkTE5VTUIiIup6IWEXG5oMyjlv9WVVPH7v2V7C6roqS8iv2VtVRU11BRVUd5VQ2VNXUcespZYwwAMZFeYqMOXiJoHuUlMTaStLhokptHEeHV31qRcKeiDpA9FdVsKNrP5uJyCorL2VLi+7i99ABFZZXsPVAT8Oc0BlrERpHSPIpWiTG0TY79zyUrOZZ2KbHEx0QG/HlFpHGpqBuots6yZtc+Vmzby+qd+1i9w3fZvufA9x6XGhdFZotYurdOIDUuipS4aFLiokiNi6ZFbNT3tpKbRXmJifTg8W9FH9ywrrOWyuo6yqtrfFvgVbWUVdawp6KKwrIqdpdVUlRWSdG+KrbvqeD9pdspKa/+Xo6MpGZ0bRVP11bxdGsVT/fWCXRMi8PrMY3yeonIiVNRH0Px/ioWFpSwsKCUBQUlLN5cyv6qWgCivB46toxjSIcUuraKp2NaHG2TY8ls0Yzm0YF5aWMivSRS/63ifQeq2VxcQUFxOesKy/7zh+SLbwupqfP9BYiLjqBPViL9slrQr20SfbOSSImLDkheEQk8FfVh9lfWMG9jMXPWFDFn3W5Wbt8LgNdj6N46nov6Z9KvbRK9MhLJTm1OpMvGiONjIslpE0lOm4Tv3V5VU8f6ojKWb93Los2lLNxcwqRZ66j1l3fX9HhO6pTCsI6pDO6QrCETERcxwVgzMTc314bSSZnWF5Yxc8VOPl21i4UFJVTXWqIiPOS2a8GwTqkMzE6mV0YizaK8TkcNqIqqWpZu3UPepmK+Wrub+RuLqaypw+sx9M5M5PSuLTmrezrdW8f/Z+emiASHMSbfWpt7xPuaYlHX1VkWbi5l5oqdzFyxg3WF+wHIaZ3AqV3SOLlTKrnZLYiJDK9iPpYD1bUsKCjhq7W7+XJtEUu2lGKtb5x7eE46Z3VPZ1D7ZKIi3PUuQiQcqKj9Vm7fyzuLtvGvxdvYWlpBhMcwuEMyZ+e04qycdDKSmjkd0VV27TvAZ6t2MXPFTr5cU0RlTR1JsZGc27MVo/pkMKh9snZKigRIky7qbaUVvLVwK/9ctI3VO/fh9RhO7ZzKqL5tOKNbOonNNBZbHxVVtXy5ppD3l27n4xU7Ka+qJT0hmpG923B+3zb0ykjU8IjICWhyRV1TW8fnqwt5aV4Bn6/eRZ2F3HYtOL9vG87r1VozHE5QeVUNn67cxT8Xb2PW6kKqauvo1iqeKwZmcWG/TBJj9cdPpKGaTFFvK63g5fmbeXX+ZnbsPUBafDSX52ZxWW4WbVNiGz1PU7CnvJp3l27j5XmbWbp1D1ERHs7r2YorBrVlcPtkbWWL1FPYF/WizaVM+XI9HyzbQZ21nNo5jdGD2nJm95aumz4XzpZt3cMr8zfz9sKt7KusoUt6HNcOa8+F/TKa3I5ZkYYKy6KurbPMXLGDKV9uIG9TCfHREYwe3JaxQ9qRlaytZydVVNXyryXbeHbORlZs30uL2EiuHNyWsUOyaZUY43Q8EVcKq6KurKnljfytPDlrHQXF5WQlN+Pak9pz2cAs4gJ0NKAEhrWWuRuKmTZ7AzNX7sRrDCN7t+am0zrRtVW80/FEXOWHijpkmu1AdS2v5W1m0ufr2LbnAH2ykrj73G6c3aOVpoi5lDGGIR1SGNIhhYLd5Tz71UZemV/A24u2cXZOOjef0YnemUlOxxRxPddvUR+ormXGvAKenLWOnXsryW3XglvP7MwpnVO1oyoElZZX8cycjTwzZwN7D9Rwapc0bj69E4PaJzsdTcRRJzz0YYzZCOwDaoGao32zgwJR1DW1dbyxYAt/m7mGHXsPMKRDMree2ZmhHVJU0GFg34FqXvymgKmz11NUVsXQDincMaIr/du2cDqaiCMCVdS51tqi+jzhiRS1tZaPV+xk4kerWburjH5tk/jNiG4M6ZByXN9P3K2iqpaX5hXwxGdr2b2/iuE56fz67K4aw5YmJ2SKet6GYh78YCULCkrpmNacO0d04+ycdG1BNwH7K2uYNnsDk79YT1lVDRf0zeCXZ3WmXUpzp6OJNIpAFPUGoASwwFPW2sk/9PjjKeo9FdUMfeBT4mMi+NVZXbhkQKaWmWqCSsurmDRrHc/O2UidtVwzNJtbzuysQ/0l7AWiqNtYa7cZY1oCM4FbrLVfHPaYccA4gLZt2w7YtGlTg4PO31hMzzbhdzpRabidew/w8Mff8mr+ZlrERnHb8C5cMTBLf7wlbAV0HrUx5n6gzFr7l6M9xulzfUj4WLZ1D79/dwXzNhTTNT2e/xnZnVM6pzkdSyTgfqioj7l5YoxpboyJP3gdOBtYFtiIIkfWMyORV8YNYdKY/pRX1zB26jyuey6PLSXlTkcTaTT1eR+ZDsw2xiwG5gHvWWs/DG4ske8YYzi3V2tm/upH3DmiK3PWFjH84S+Y9Pk6qmvrnI4nEnSuP+BF5HBbSyu4/5/LmbliJ13S4/jjBb10wIyEvBMa+hBxm4ykZjx9dS5PX53L/spaLnvqa+54bTHF+6ucjiYSFCpqCVnDc9KZedup3Pijjry1cCvDH57F+0u3Ox1LJOBU1BLSYqMiuOvcbvzrlpNpnRTD+OkLGD89n6KySqejiQSMilrCQvfWCbw9fhh3nNOVT1bsYvjDs3hn0VaCsQ9GpLGpqCVsRHg9TDi9E+/dejJtU5rzi5cXMe6FfHbtO+B0NJEToqKWsNM5PZ43bhzK3ed2Y9a3hYx45Es+WbHT6Vgix01FLWEpwuvhhh915L1bTqZVQgzXPZ/HPW8tpbyqxuloIg2mopaw1jk9nrcmnMQNP+rAjHkFjHx0Nku2lDodS6RBVNQS9qIjvNx9bnemXzeYiupaLnriK/7x2Vpq67SjUUKDilqajJM6pvLhL05lRM9WTPxoNVdNmasdjRISVNTSpCTGRvLY6H5MvKQ3CzeXcN7fZzNnbb3WwxBxjIpamhxjDJfmZvHPm08mKTaSq6bO5ZFPvtVQiLiWilqarC7p8bwzYRgX9s3gkU/WMHaqhkLEnVTU0qQ1j47gr5f14aGLe5O/yTcU8vW63U7HEvkeFbU0ecYYLhuYxTs3DyOhWQRXTZ3LtNkbdPi5uIaKWsSvW6sE3pkwjNO7tuT3767gtlcXU1FV63QsERW1yKHiYyKZPHYAtw3vwtuLtnLJk1+xuVjLfomzVNQih/F4DLee2ZkpV+dSsLucUY9rCp84S0UtchRndk/nnZuHkRoXzdipc5ny5XqNW4sjVNQiP6BDWhxvTRjG8Jx0/vjeSu55a6kW1JVGp6IWOYa46AgmjRnA+NM6MmPeZq6ZNo/Scq3PKI1HRS1SDx6P4c4R3fjrpX2Yv7GYC5/4ivWFZU7HkiZCRS3SABcPyOSl64ewp6KaC5/4iq/WaSejBJ+KWqSBBmYn8/b4YaTFR3P11Hm8PK/A6UgS5lTUIsehbUosb44/iZM6pXLXm0t5+OPVmhEiQaOiFjlOCTGRTLsml8tzs3j032u58/UlmhEiQRHhdACRUBbh9fDgxb1olRjD3z9dw659lTwxpj/No/WrJYFT7y1qY4zXGLPQGPNuMAOJhBpjDL8a3oUHL+rF7LVFXDH5Gwr3VTodS8JIQ4Y+fgGsDFYQkVB3xaC2PH31ANbuKuOiSXM0fU8Cpl5FbYzJBH4MTAluHJHQdka3dGaMG8L+ylounvQVCwtKnI4kYaC+W9SPAHcC2lMicgx9s5J486aTiI+JZMyUuTqhk5ywYxa1MWYksMtam3+Mx40zxuQZY/IKCwsDFlAkFGWnNuf1G4eS1SKWa5+Zz4fLdjgdSUJYfbaohwGjjDEbgZeBM4wxLx7+IGvtZGttrrU2Ny0tLcAxRUJPy4QYXrlhCDltEhg/PZ/X87c4HUlC1DGL2lp7t7U201qbDVwB/Ntae1XQk4mEgaTYKKZfN5iTOqby69cW88ycDU5HkhCkA15Egqx5dARTf5rLiB6t+N9/reCRT77VUYzSIA0qamvt59bakcEKIxKuoiO8PH5lPy4ZkMkjn6zhD++uVFlLvenwKZFGEuH18NDFvX2Hns/ZQGVNLX84vycej3E6mricilqkEXk8hvtGdic60sOkz9dRU2t54KJeKmv5QSpqkUZmjOHOc7oS6fXw6KdrqK6rY+IlffCqrOUoVNQiDjDGcNvwLkR4DA/P/JbaOstfL+1DhFf79+W/qahFHHTrmZ2J9Hr484erqKm1PHJFXyJV1nIYFbWIw246rSORXsMf31tJTV0dj43uT1SEylq+o58GERe47pQO3P+THD5avpPx0xdQVaPT6sh3VNQiLvHTYe35w/k9+GTlTm6dsVCrxch/qKhFXGTs0GzuG5nDh8t3cNuri6lRWQsaoxZxnZ+f3J6a2joe+GAVkR7DxEs1da+pU1GLuNANP+pITZ1l4kerifAaHryotw6KacJU1CIuNeH0TlTV1PH3T9cQ4fXwpwt6YozKuilSUYu42C/P6kx1bR1PfL6OKK+H3/0kR2XdBKmoRVzMGMMd53SluraOp7/cQHSEh7vO7aaybmJU1CIuZ4zhnvO6U1lTx1NfrCehWSQTTu/kdCxpRCpqkRBgjOH+n/Rg34EaJn60moSYCMYOzXY6ljQSFbVIiPB4DA9d0pt9B2q4753lxMVEcGG/TKdjSSPQAS8iISTS6+HxK/sxtEMKv35tCTNX7HQ6kjQCFbVIiImJ9PL0Nbn0zEhkwksL+GpdkdORJMhU1CIhKC46gueuHUh2SizXP5fHos2lTkeSIFJRi4SopNgoXvj5YFLiovnpM/NYvWOf05EkSFTUIiEsPSGG6dcNJjrCwzXT5rG1tMLpSBIEKmqREJeVHMtzPxvE/qoarpk2j9LyKqcjSYCpqEXCQLdWCTx9dS4FxeX8/Lk8KqpqnY4kAaSiFgkTQzqk8PfL+7KgoIRbZizUuazDiIpaJIyc26s1vx/lWyXmvneWYa11OpIEgI5MFAkzY4dms3NvJY9/tpaW8TH8angXpyPJCTpmURtjYoAvgGj/41+31v4u2MFE5PjdfnYXdu07wN8/XUNafDRXDWnndCQ5AfXZoq4EzrDWlhljIoHZxpgPrLXfBDmbiBwnYwz/d2Evisqq+O07y0iNi2ZEz1ZOx5LjdMwxautT5v800n/RwJeIy0X4zwvSOzOJW19eSP6mEqcjyXGq185EY4zXGLMI2AXMtNbOPcJjxhlj8owxeYWFhQGOKSLHIzYqgmk/HUibxBiufz6PTbv3Ox1JjkO9itpaW2ut7QtkAoOMMT2P8JjJ1tpca21uWlpagGOKyPFKbh7FM9cOos5arn1mPiX7dUBMqGnQ9DxrbSnwOTAiGGFEJDjapzbn6atz2VJSwQ0v5FNZowNiQskxi9oYk2aMSfJfbwacBawKci4RCbCB2cn85bI+zNtYzJ2vL9Ec6xBSn1kfrYHnjDFefMX+qrX23eDGEpFgGNWnDZuLy5n40WraJsdy+9ldnY4k9XDMorbWLgH6NUIWEWkE40/rSMHuch7791qyWsRy2cAspyPJMejIRJEmxhjDHy/sybY9Fdzz1lLaJDXj5M6pTseSH6BzfYg0QZFeD/8Y05+OaXHc9GK+Fh1wORW1SBOVEBPJtGsHEhPl5WfPzqdwX6XTkeQoVNQiTVhGUjOmXTOQ3fsrufHFfA5Ua9qeG6moRZq4XpmJ/PXSvuRvKuGeN5dq2p4LqahFhB/3bs1tw7vw5sKtTJq1zuk4chjN+hARAG45oxNrdpXx0Ier6ZgWxzk9dLY9t9AWtYgAvml7Ey/pTZ/MRH71yiKWb9vjdCTxU1GLyH/ERHp5+upcEmIiuf65PHbtO+B0JEFFLSKHaZkQw5Rrcikpr2bc85oJ4gYqahH5Lz0zEvnb5X1YtLmUu97QCZycpqIWkSMa0bM1vz67C28v2sYTn2smiJM060NEjmrC6b6ZIBM/Wk3nlnGcrZkgjtAWtYgclTGGP1/cm97+mSBrduqcIE5QUYvID4qJ9PLU2AE0i/Jy/fN57KmodjpSk6OiFpFjap3YjElXDWBraQW3zlhIbZ12LjYmFbWI1MvA7GTuH9WDWd8WMvGj1U7HaVK0M1FE6m3M4HYs37aXJ2etI6dNAqP6tHE6UpOgLWoRaZD7f9KD3HYtuPP1xTrMvJGoqEWkQaIiPDxxVX+SmkUx7vl8dpdpwYFgU1GLSIO1jI9h8tUDKCyrZMJLC6iurXM6UlhTUYvIcemdmcQDF/bim/XF/Om9lU7HCWvamSgix+3iAZks37aXaXM2kNMmgctys5yOFJa0RS0iJ+Se87oxrFMK//P2MpZsKXU6TlhSUYvICYnwenhsdH/S4qK56cUFFO+vcjpS2FFRi8gJS24exaSr+lNYVqkjF4NARS0iAdE7M4k/nN+D2WuL+OvHOnIxkI5Z1MaYLGPMZ8aYlcaY5caYXzRGMBEJPZcPbMvoQVk88fk6Plq+w+k4YaM+W9Q1wO3W2u7AEGCCMSYnuLFEJFT97ic96J2ZyK9fXcz6wjKn44SFYxa1tXa7tXaB//o+YCWQEexgIhKaYiK9TLpqAJERHm54IZ/9lTVORwp5DRqjNsZkA/2AuUe4b5wxJs8Yk1dYWBigeCISijKSmvHY6H6sKyzjN1pz8YTVu6iNMXHAG8AvrbV7D7/fWjvZWptrrc1NS0sLZEYRCUHDOqVyxzndeHfJdqbO3uB0nJBWr6I2xkTiK+np1to3gxtJRMLFjT/qwDk90nngg1V8s36303FCVn1mfRhgKrDSWvtw8COJSLgwxvCXS/vQLiWWm19awI49B5yOFJLqs0U9DBgLnGGMWeS/nBfkXCISJuJjInnqqgGUV9Uyfno+VTU6015D1WfWx2xrrbHW9rbW9vVf3m+McCISHjqnxzPxkj4sKCjlj++tcDpOyNGRiSLSKH7cuzXXn9Ke57/exDuLtjodJ6SoqEWk0dw5ohsDs1tw95tLWbtrn9NxQoaKWkQaTaTXw+NX9ic2ysuNLy7QwTD1pKIWkUaVnhDD36/ox/rCMu59a6kOhqkHFbWINLphnVK5bXgX3l60jelzC5yO43oqahFxxPjTOnFa1zR+/68VWhnmGFTUIuIIj8fwt8v6khbvWxmmtFwrwxyNilpEHNOieRT/GNOfXfsOcPuri6nTyjBHpKIWEUf1zUrivpE5fLpqF09+sc7pOK6kohYRx40d0o6f9GnDXz5azdfrdPKmw6moRcRxxhgeuKgX7VObc8uMhezaq5M3HUpFLSKuEBcdwaSrBrC/soabX1pITa1O3nSQilpEXKNLejwPXNSLeRuLmaiVzP9DRS0irnJBvwzGDG7LU7PW87FWMgdU1CLiQveNzKFXRiK3v7aYgt3lTsdxnIpaRFwnJtLLE2P6Y4CbpudzoLrW6UiOUlGLiCtlJcfy8GV9Wb5tL//7r+VOx3GUilpEXOusnHRuOq0jM+Zt5o38LU7HcYyKWkRc7fbhXRjcPpl7317Kqh17nY7jCBW1iLhahNfDY1f2Iz4mkvEvLqCsCS42oKIWEddrGR/DY6P7sXH3fn7zxpImt9iAilpEQsKQDinccU433luynee+2uh0nEalohaRkHHDqR04q3tL/vT+ShYUlDgdp9GoqEUkZHg8hr9e2pf0hBhunr6A4v1NY7EBFbWIhJTE2EgmjRlAUVkVv3xlUZNYbEBFLSIhp1dmIr8blcMX3xby+GdrnY4TdCpqEQlJVw5qy4X9MvjbJ98ye02R03GC6phFbYyZZozZZYxZ1hiBRETqwxjDny7sSeeWcdz68kK276lwOlLQ1GeL+llgRJBziIg0WGxUBE+MGcCB6lpufmkh1WG62MAxi9pa+wVQ3AhZREQarFPLOB68uDf5m0r48wernI4TFAEbozbGjDPG5Blj8goLCwP1bUVEjmlUnzZcM7QdU2Zv4MNl252OE3ABK2pr7WRrba61NjctLS1Q31ZEpF7u+XF3+mQlccdrS9hQtN/pOAGlWR8iEhaiI7z848p+eL2Gm14Mr8UGVNQiEjYyW8Tyt8v7smrHPn77TvhMVKvP9LwZwNdAV2PMFmPMz4MfS0Tk+JzetSW3nNGJV/O28Or8zU7HCYiIYz3AWju6MYKIiATKL8/qQv6mEu57Zxk9MxLJaZPgdKQToqEPEQk7Xo/h0dH9SIqNZPz0fPYeqHY60glRUYtIWEqNi+bxK/uzuaSCO18L7cUGVNQiErYGZidz14hufLh8B1Nnb3A6znFTUYtIWLvulPacnZPOgx+sIm9jaB5kraIWkbBmjGHipX3IaNGMm19aSFFZpdORGkxFLSJhL7FZJE+M6U9JeRU3v7Qg5E7epKIWkSahR5tEHrioF9+sL+b/3l/pdJwGOeY8ahGRcHFR/0yWbd3LtDkb6NkmkYsHZDodqV60RS0iTco953VjaIcU7n5rKUu2lDodp15U1CLSpER4PTx+ZT/S4qK54YX8kNi5qKIWkSYnJS6ap8YOoHh/FeOnu3/noopaRJqknhmJ/Pni3szbUMwf313hdJwfpJ2JItJkXdAvg2Vb9zBl9gZ6ZiRyaW6W05GOSFvUItKk3XVuN4Z1SuHet5eRv8mdRy6qqEWkSYvwenh8dH/aJMYw7vl8NheXOx3pv6ioRaTJa9E8iqk/HUh1bR0/e3a+606LqqIWEQE6psXx5NgBbCjaz4TpC6hx0UwQFbWIiN9JHVP5vwt78eWaIn73z+WuOYe1Zn2IiBzisoFZrCsq46lZ62mf2pzrTungdCQVtYjI4X5zTjc2FZXzp/dX0jIhhlF92jiaR0MfIiKH8XgMj1zRl4HZydz+6iK++LbQ2TyOPruIiEvFRHqZck0unVrGc+OL+SwsKHEsi4paROQoEmIiee5nA0mNi+baZ+ezZuc+R3KoqEVEfkDL+Bhe+PkgIr0eRj89l28dKGsVtYjIMbRLac6M64fgMTB68jes2rG3UZ9fRS0iUg+dWsbx8rghRHgNoyd/06iLDqioRUTqqUNaHK+MG0psVASXP/UNHy/f0SjPW6+iNsaMMMasNsasNcbcFexQIiJulZ3anLcnDKNLq3hueDGff3y2ltq64B7BeMyiNsZ4gX8A5wI5wGhjTE5QU4mIuFhafDQvXz+Ekb3bMPGj1Vz+1Ncs37YnaM9Xny3qQcBaa+16a20V8DJwftASiYiEgGZRXh69oi+PXN6XNbvK+PGjs7nsqa+DsqxXfQ4hzwA2H/L5FmDw4Q8yxowDxgG0bds2IOFERNzMGMMF/TI4vVtLps/dRMHuciK9gd/1V5+iNke47b8GZKy1k4HJALm5ue445ZSISCNIbBbJ+NM6Be3716f6twCHLiSWCWwLThwRETlcfYp6PtDZGNPeGBMFXAH8M7ixRETkoGMOfVhra4wxNwMfAV5gmrV2edCTiYgIUM/zUVtr3wfeD3IWERE5Ah2ZKCLicipqERGXU1GLiLicilpExOVMMJZDN8YUApuO88tTgaIAxgkU5Wo4t2ZTroZRroY7nmztrLVpR7ojKEV9IowxedbaXKdzHE65Gs6t2ZSrYZSr4QKdTUMfIiIup6IWEXE5Nxb1ZKcDHIVyNZxbsylXwyhXwwU0m+vGqEVE5PvcuEUtIiKHUFGLiLica4rayQV0jTFZxpjPjDErjTHLjTG/8N+ebIyZaYxZ4//Y4pCvudufdbUx5pwg5/MaYxYaY951Wa4kY8zrxphV/tduqBuyGWN+5f9/XGaMmWGMiXEilzFmmjFmlzFm2SG3NTiHMWaAMWap/75HjTFHWswjENkm+v8vlxhj3jLGJDV2tiPlOuS+XxtjrDEm1S25jDG3+J97uTHmoaDlstY6fsF3+tR1QAcgClgM5DTi87cG+vuvxwPf4lvI9yHgLv/tdwF/9l/P8WeMBtr7s3uDmO824CXgXf/nbsn1HHCd/3oUkOR0NnxLx20Amvk/fxX4qRO5gFOB/sCyQ25rcA5gHjAU32pLHwDnBinb2UCE//qfnch2pFz+27PwnWp5E5DqhlzA6cAnQLT/85bByuWWLWpHF9C11m631i7wX98HrMT3C38+vjLC//EC//XzgZettZXW2g3AWv+/IeCMMZnAj4Eph9zshlwJ+H54pwJYa6ustaVuyIbv9L3NjDERQCy+FYkaPZe19gug+LCbG5TDGNMaSLDWfm19v+nPH/I1Ac1mrf3YWlvj//QbfKs5NWq2o7xmAH8D7uT7ywA6nesm4EFrbaX/MbuClcstRX2kBXQznAhijMkG+gFzgXRr7XbwlTnQ0v+wxsz7CL4f0EOXNnZDrg5AIfCMf1hmijGmudPZrLVbgb8ABcB2YI+19mOncx2ioTky/NcbK99BP8O3xed4NmPMKGCrtXbxYXc5/Zp1AU4xxsw1xswyxgwMVi63FHW9FtANeghj4oA3gF9aa/f+0EOPcFvA8xpjRgK7rLX59f2SI9wWrNcxAt9bwUnW2n7Afnxv5Y+msV6zFvi2aNoDbYDmxpirnM5VD0fL0ej5jDH3AjXA9IM3HSVD0LMZY2KBe4HfHulup3L5RQAtgCHAHcCr/jHngOdyS1E7voCuMSYSX0lPt9a+6b95p//tCv6PB9/aNFbeYcAoY8xGfMNBZxhjXnRBroPPtcVaO9f/+ev4itvpbGcBG6y1hdbaauBN4CQX5DqooTm28N0QRNDzGWOuAUYCY/xvz53O1hHfH93F/t+DTGCBMaaVw7nwP8+b1mcevne9qUHJdSID7IG64PvLtB7ff8jBnYk9GvH5Db7xokcOu30i39/x85D/eg++v7NgPUHcaed/ztP4bmeiK3IBXwJd/dfv9+dyNBswGFiOb2za4BsHvsWpXEA2398B1eAc+BaYHsJ3O6DOC1K2EcAKIO2wxzVqtsNzHXbfRr7bmehoLuBG4Pf+613wDXeYYOQKyi/wcb4I5+GbbbEOuLeRn/tkfG9BlgCL/JfzgBTgU2CN/2PyIV9zrz/ragKwF74eGU/ju6J2RS6gL5Dnf93exvc20PFswP8Cq4BlwAv+X5hGzwXMwDdOXo1va+rnx5MDyPX/W9YBj+M/ojgI2db6y+bg78CTjZ3tSLkOu38j/qJ2Ohe+jcoX/c+zADgjWLl0CLmIiMu5ZYxaRESOQkUtIuJyKmoREZdTUYuIuJyKWkTE5VTUIiIup6IWEXG5/weyiGKGItOrkQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] @@ -137,24 +166,26 @@ "source": [ "V = Ausgleichsbecken_class(1.,0.9,0.,10.,timestep=0.001)\n", "\n", - "V.set_initial_level(10.) \n", - "V.set_influx(3.2)\n", - "V.set_outflux(0.)\n", + "V.set_initial_level(initial_level) \n", + "V.set_influx(initial_influx)\n", + "V.set_outflux(initial_outflux)\n", "\n", - "V.p0 = 0.\n", + "V.p0 = initial_pipeline_pressure\n", "\n", "outflux_vec = []\n", "level_vec = []\n", "\n", "t = 0\n", - "while V.level > 0.005:\n", + "while V.level > total_min_level:\n", " t = t+V.timestep\n", " V.e_RK_4()\n", " V.level = V.update_level(V.timestep)\n", " V.set_volume()\n", " outflux_vec.append(V.outflux)\n", " level_vec.append(V.level)\n", - " \n", + " if t > total_max_time:\n", + " break\n", + "\n", "plt.plot(level_vec)\n" ] }