OT-2_BalanceAutomation.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "2f39c2fa",
   "metadata": {
    "pycharm": {
     "name": "#%% md\n"
    }
   },
   "source": [
    "# Opentrons Mass Balance Automation \n",
    "\n",
    "## Aniket Chitre"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9da6671f",
   "metadata": {
    "pycharm": {
     "name": "#%% md\n"
    }
   },
   "source": [
    "### 0. Import Packages"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "f44fb10b",
   "metadata": {
    "pycharm": {
     "name": "#%%\n"
    }
   },
   "outputs": [],
   "source": [
    "%matplotlib notebook\n",
    "import matplotlib.pyplot as plt\n",
    "import seaborn as sns\n",
    "import pandas as pd\n",
    "import time\n",
    "import serial\n",
    "import csv"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "910e04af",
   "metadata": {
    "pycharm": {
     "name": "#%% md\n"
    }
   },
   "source": [
    "### 1. Port for Mass Balance\n",
    "\n",
    "Find the port number on Windows through Device Manager or on Mac through running this command (`ls /dev/tty.*`) on the Terminal"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "1b2454b1",
   "metadata": {
    "pycharm": {
     "name": "#%%\n"
    }
   },
   "outputs": [],
   "source": [
    "#port = '/dev/tty.usbserial-14130'\n",
    "port = 'COM5'"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "dda217b1",
   "metadata": {
    "pycharm": {
     "name": "#%% md\n"
    }
   },
   "source": [
    "### 2. Function to Read Mass from Balance"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "6692b754",
   "metadata": {
    "pycharm": {
     "name": "#%%\n"
    }
   },
   "outputs": [],
   "source": [
    "def get_mass(port):\n",
    "    try:\n",
    "        balance = serial.Serial(port, 9600, xonxoff=True) # check COM port and xonxoff (handshake) = True default for MT balance\n",
    "        balance.write('S\\r\\n'.encode('utf-8'))\n",
    "        time.sleep(1)\n",
    "        result  = balance.read_until('\\r\\n'.encode('utf-8'))\n",
    "        value = str(result[7:14].decode(\"utf-8\"))\n",
    "        #balance.close()\n",
    "        if (value.strip() != '-------') and (value.strip() != ''):\n",
    "            return (True,float(value))\n",
    "        return (False,-1)\n",
    "    except:\n",
    "        return (False,-1)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "4a9656db",
   "metadata": {},
   "source": [
    "### 3. Real-time Data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "06a5b1b9",
   "metadata": {
    "pycharm": {
     "name": "#%%\n"
    },
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "0.0\n",
      "-0.001\n",
      "-0.001\n",
      "-0.001\n",
      "-0.001\n",
      "0.0\n",
      "-0.002\n",
      "-0.002\n",
      "0.98\n",
      "0.983\n",
      "0.985\n",
      "0.982\n",
      "0.983\n",
      "1.966\n",
      "1.966\n",
      "1.966\n",
      "1.966\n",
      "2.946\n",
      "2.949\n",
      "2.95\n",
      "2.95\n",
      "2.949\n",
      "3.93\n",
      "3.932\n",
      "3.933\n",
      "3.933\n",
      "3.932\n",
      "4.913\n",
      "4.916\n",
      "4.917\n",
      "4.917\n",
      "4.917\n",
      "5.902\n",
      "5.901\n",
      "5.901\n",
      "5.902\n",
      "6.588\n",
      "6.591\n",
      "6.592\n",
      "6.591\n",
      "6.591\n",
      "6.591\n",
      "7.275\n",
      "7.273\n",
      "7.275\n",
      "8.256\n",
      "8.261\n",
      "8.261\n",
      "8.26\n",
      "8.26\n",
      "9.244\n",
      "9.246\n",
      "9.248\n",
      "9.247\n",
      "9.246\n",
      "10.23\n",
      "10.232\n",
      "10.234\n",
      "10.231\n",
      "10.234\n",
      "11.222\n",
      "11.226\n",
      "11.226\n",
      "11.225\n",
      "11.225\n",
      "12.209\n",
      "12.211\n",
      "12.213\n",
      "12.213\n",
      "12.213\n",
      "13.195\n",
      "13.197\n",
      "13.198\n",
      "13.203\n",
      "13.202\n",
      "13.758\n",
      "13.761\n",
      "13.761\n",
      "13.761\n",
      "14.319\n",
      "14.322\n",
      "14.323\n",
      "14.323\n",
      "14.323\n",
      "15.311\n",
      "15.313\n",
      "15.314\n",
      "15.314\n",
      "15.313\n",
      "16.297\n",
      "16.299\n",
      "16.301\n",
      "16.299\n",
      "16.299\n",
      "17.286\n",
      "17.286\n",
      "17.286\n",
      "17.286\n",
      "18.273\n",
      "18.273\n",
      "18.272\n",
      "18.271\n",
      "19.258\n",
      "19.258\n",
      "19.258\n",
      "19.257\n",
      "20.248\n",
      "20.248\n",
      "20.247\n",
      "20.247\n",
      "20.936\n",
      "20.936\n",
      "20.936\n",
      "20.936\n",
      "21.65\n",
      "21.649\n",
      "21.646\n",
      "21.648\n",
      "22.637\n",
      "22.637\n",
      "22.634\n",
      "22.635\n",
      "23.624\n",
      "23.624\n",
      "23.623\n",
      "23.623\n",
      "24.618\n",
      "24.618\n",
      "24.617\n",
      "24.617\n",
      "25.605\n",
      "25.605\n",
      "25.603\n",
      "25.602\n",
      "26.591\n",
      "26.59\n",
      "26.59\n",
      "26.59\n",
      "27.577\n",
      "27.576\n",
      "27.575\n",
      "27.575\n",
      "28.469\n",
      "28.469\n",
      "28.468\n",
      "28.468\n",
      "29.34\n",
      "29.34\n",
      "29.339\n",
      "29.339\n",
      "30.33\n",
      "30.33\n",
      "30.329\n",
      "30.329\n",
      "31.319\n",
      "31.319\n",
      "31.317\n",
      "31.317\n",
      "32.302\n",
      "32.302\n",
      "32.302\n",
      "32.302\n",
      "33.287\n",
      "33.287\n",
      "33.286\n",
      "33.285\n",
      "34.276\n",
      "35.26\n",
      "35.258\n",
      "35.258\n",
      "35.26\n",
      "36.041\n",
      "36.041\n",
      "36.04\n",
      "36.817\n",
      "36.817\n",
      "36.816\n",
      "36.816\n",
      "37.808\n",
      "37.808\n",
      "38.793\n",
      "38.791\n",
      "38.79\n",
      "38.79\n",
      "39.776\n",
      "39.776\n",
      "39.775\n",
      "39.774\n",
      "40.759\n",
      "40.759\n",
      "40.758\n",
      "40.758\n",
      "41.754\n",
      "41.753\n",
      "41.752\n",
      "41.752\n",
      "42.755\n",
      "42.755\n",
      "42.754\n",
      "43.266\n",
      "43.266\n",
      "43.264\n",
      "43.264\n",
      "43.777\n",
      "43.774\n",
      "43.773\n",
      "43.773\n",
      "43.774\n",
      "43.774\n",
      "43.772\n",
      "43.772\n",
      "43.772\n",
      "43.774\n",
      "43.775\n",
      "43.775\n",
      "43.774\n",
      "43.772\n",
      "43.772\n",
      "43.772\n",
      "43.772\n",
      "43.772\n",
      "43.772\n",
      "43.771\n",
      "43.768\n",
      "43.769\n",
      "43.769\n",
      "43.769\n",
      "43.769\n",
      "43.769\n",
      "43.769\n",
      "43.769\n",
      "43.769\n",
      "43.768\n",
      "43.768\n",
      "43.768\n",
      "43.768\n",
      "43.768\n",
      "43.768\n",
      "43.768\n",
      "43.768\n",
      "43.768\n",
      "43.768\n",
      "43.768\n",
      "44.39\n",
      "44.39\n",
      "44.39\n",
      "44.41\n",
      "44.41\n",
      "44.41\n",
      "44.41\n",
      "44.409\n",
      "44.439\n",
      "44.462\n",
      "44.462\n",
      "44.506\n",
      "44.507\n",
      "44.525\n",
      "44.524\n",
      "44.53\n",
      "44.53\n",
      "44.53\n",
      "44.53\n",
      "44.53\n",
      "44.528\n",
      "44.528\n",
      "44.528\n",
      "44.528\n",
      "44.528\n",
      "44.527\n",
      "44.527\n",
      "44.527\n",
      "44.527\n",
      "44.527\n",
      "44.527\n",
      "44.528\n",
      "44.528\n",
      "44.528\n",
      "44.527\n",
      "44.527\n",
      "44.526\n",
      "45.045\n",
      "45.045\n",
      "45.045\n",
      "45.045\n",
      "45.045\n",
      "45.045\n",
      "45.07\n",
      "45.07\n",
      "45.094\n",
      "45.094\n",
      "45.118\n",
      "45.118\n",
      "45.142\n",
      "45.142\n",
      "45.165\n",
      "45.165\n",
      "45.165\n",
      "45.164\n",
      "45.164\n",
      "45.164\n",
      "45.163\n",
      "45.163\n",
      "45.163\n",
      "45.161\n",
      "45.16\n",
      "45.161\n",
      "45.161\n",
      "45.161\n",
      "45.161\n",
      "45.161\n",
      "45.161\n",
      "45.161\n",
      "45.161\n",
      "45.161\n",
      "45.161\n",
      "45.162\n",
      "45.16\n",
      "45.64\n",
      "45.64\n",
      "45.64\n",
      "45.64\n",
      "45.64\n",
      "45.64\n",
      "45.663\n",
      "45.663\n",
      "45.688\n",
      "45.688\n",
      "45.711\n",
      "45.711\n",
      "45.734\n",
      "45.734\n",
      "45.753\n",
      "45.758\n",
      "45.758\n",
      "45.758\n",
      "45.757\n",
      "45.757\n",
      "45.756\n",
      "45.755\n",
      "45.754\n",
      "45.754\n",
      "45.755\n",
      "45.756\n",
      "45.756\n",
      "45.756\n",
      "45.754\n",
      "45.754\n",
      "45.754\n",
      "45.754\n",
      "45.753\n",
      "45.753\n",
      "45.753\n",
      "45.753\n",
      "45.753\n",
      "45.753\n",
      "45.751\n",
      "45.751\n",
      "45.751\n",
      "45.751\n",
      "45.751\n",
      "45.751\n",
      "45.75\n",
      "45.75\n",
      "45.75\n",
      "45.752\n",
      "45.752\n",
      "45.748\n",
      "45.749\n",
      "45.749\n",
      "45.749\n",
      "45.749\n",
      "45.748\n",
      "46.344\n",
      "46.345\n",
      "46.346\n",
      "46.346\n",
      "46.346\n",
      "46.347\n",
      "46.344\n",
      "46.378\n",
      "46.401\n",
      "46.394\n",
      "46.417\n",
      "46.418\n",
      "46.443\n",
      "46.443\n",
      "46.437\n",
      "46.441\n",
      "46.443\n",
      "46.443\n",
      "46.442\n",
      "46.442\n",
      "46.44\n",
      "46.439\n",
      "46.44\n",
      "46.44\n",
      "46.44\n",
      "46.439\n",
      "46.439\n",
      "46.438\n",
      "46.439\n",
      "46.438\n",
      "46.438\n",
      "46.438\n",
      "46.439\n",
      "46.44\n",
      "46.44\n",
      "46.439\n",
      "46.439\n",
      "46.439\n",
      "46.437\n",
      "46.959\n",
      "46.959\n",
      "46.958\n",
      "46.958\n",
      "46.958\n",
      "46.958\n",
      "46.987\n",
      "46.987\n",
      "47.008\n",
      "47.008\n",
      "47.029\n",
      "47.029\n",
      "47.028\n",
      "47.031\n",
      "47.031\n",
      "47.03\n",
      "47.03\n",
      "47.029\n",
      "47.029\n",
      "47.029\n",
      "47.028\n",
      "47.028\n",
      "47.028\n",
      "47.028\n",
      "47.027\n",
      "47.027\n",
      "47.027\n",
      "47.027\n",
      "47.028\n",
      "47.029\n",
      "47.027\n",
      "47.028\n",
      "47.027\n",
      "47.027\n",
      "47.027\n",
      "47.026\n",
      "47.516\n",
      "47.516\n",
      "47.516\n",
      "47.514\n",
      "47.543\n",
      "47.543\n",
      "47.543\n",
      "47.565\n",
      "47.565\n",
      "47.606\n",
      "47.606\n",
      "47.609\n",
      "47.609\n",
      "47.609\n",
      "47.609\n",
      "47.607\n",
      "47.607\n",
      "47.606\n",
      "47.606\n",
      "47.605\n",
      "47.606\n",
      "47.606\n",
      "47.608\n",
      "47.605\n",
      "47.602\n",
      "47.602\n",
      "47.604\n",
      "47.604\n",
      "47.604\n",
      "47.604\n",
      "47.604\n",
      "47.604\n",
      "47.604\n",
      "47.604\n",
      "47.603\n",
      "47.603\n",
      "47.603\n",
      "47.603\n",
      "47.603\n",
      "47.603\n",
      "47.602\n",
      "47.602\n",
      "47.601\n",
      "47.601\n",
      "47.602\n",
      "47.602\n",
      "47.602\n",
      "47.602\n",
      "47.602\n",
      "47.603\n",
      "47.602\n",
      "48.398\n",
      "48.415\n",
      "48.415\n",
      "48.415\n",
      "48.415\n",
      "48.436\n",
      "48.457\n",
      "48.457\n",
      "48.482\n",
      "48.481\n",
      "48.481\n",
      "48.48\n",
      "48.497\n",
      "48.496\n",
      "48.495\n",
      "48.494\n",
      "48.494\n",
      "48.494\n",
      "48.493\n",
      "48.493\n",
      "48.493\n",
      "48.493\n",
      "48.492\n",
      "48.492\n",
      "48.492\n",
      "48.492\n",
      "48.491\n",
      "48.491\n",
      "48.491\n",
      "48.491\n",
      "48.491\n",
      "48.491\n",
      "48.491\n",
      "48.491\n",
      "48.491\n",
      "48.49\n",
      "48.491\n",
      "48.491\n",
      "48.49\n",
      "48.49\n",
      "48.49\n",
      "49.29\n",
      "49.29\n",
      "49.289\n",
      "49.334\n",
      "49.334\n",
      "49.355\n",
      "49.355\n",
      "49.374\n",
      "49.374\n",
      "49.394\n",
      "49.394\n",
      "49.393\n",
      "49.416\n",
      "49.415\n",
      "49.415\n",
      "49.414\n",
      "49.414\n",
      "49.414\n",
      "49.413\n",
      "49.412\n",
      "49.412\n",
      "49.412\n",
      "49.412\n",
      "49.411\n",
      "49.411\n",
      "49.411\n",
      "49.411\n",
      "49.411\n",
      "49.411\n",
      "49.411\n",
      "49.411\n",
      "49.41\n",
      "49.41\n",
      "49.41\n",
      "49.411\n",
      "49.411\n",
      "49.411\n",
      "49.411\n",
      "49.41\n",
      "50.192\n",
      "50.192\n",
      "50.192\n",
      "50.192\n",
      "50.217\n",
      "50.217\n",
      "50.239\n",
      "50.238\n",
      "50.259\n",
      "50.279\n",
      "50.279\n",
      "50.3\n",
      "50.3\n",
      "50.321\n",
      "50.321\n",
      "50.32\n",
      "50.32\n",
      "50.32\n",
      "50.32\n",
      "50.32\n",
      "50.319\n",
      "50.317\n",
      "50.317\n",
      "50.317\n",
      "50.318\n",
      "50.32\n",
      "50.32\n",
      "50.32\n",
      "50.318\n",
      "50.318\n",
      "50.317\n",
      "50.317\n",
      "50.316\n",
      "50.316\n",
      "50.316\n",
      "50.316\n",
      "50.316\n",
      "50.316\n",
      "50.315\n",
      "50.315\n",
      "50.315\n",
      "50.315\n",
      "50.315\n",
      "50.315\n",
      "50.315\n",
      "50.315\n",
      "50.315\n",
      "50.315\n",
      "50.882\n",
      "50.882\n",
      "50.882\n",
      "50.881\n",
      "50.881\n",
      "50.897\n",
      "50.897\n",
      "50.897\n",
      "50.896\n",
      "50.895\n",
      "50.895\n",
      "50.895\n",
      "50.895\n",
      "50.895\n",
      "50.895\n",
      "50.896\n",
      "50.896\n",
      "50.895\n",
      "50.895\n",
      "50.895\n",
      "50.894\n",
      "50.894\n",
      "50.894\n",
      "50.894\n",
      "50.894\n",
      "50.892\n",
      "50.892\n",
      "50.892\n",
      "50.892\n",
      "50.892\n",
      "50.892\n",
      "50.892\n",
      "50.892\n",
      "50.892\n",
      "50.892\n",
      "50.892\n",
      "50.892\n",
      "50.891\n",
      "50.891\n",
      "51.514\n",
      "51.514\n",
      "51.514\n",
      "51.514\n",
      "51.533\n",
      "51.533\n",
      "51.533\n",
      "51.533\n",
      "51.533\n",
      "51.533\n",
      "51.559\n",
      "51.559\n",
      "51.559\n",
      "51.558\n",
      "51.558\n",
      "51.558\n",
      "51.558\n",
      "51.558\n",
      "51.558\n",
      "51.557\n",
      "51.557\n",
      "51.555\n",
      "51.555\n",
      "51.554\n",
      "51.556\n",
      "51.557\n",
      "51.557\n",
      "51.556\n",
      "51.556\n",
      "51.555\n",
      "51.555\n",
      "51.554\n",
      "51.553\n",
      "51.553\n",
      "51.553\n",
      "51.553\n",
      "51.553\n",
      "51.553\n",
      "51.553\n",
      "51.553\n",
      "51.553\n",
      "51.553\n",
      "51.553\n",
      "51.553\n",
      "51.553\n",
      "51.553\n",
      "51.553\n",
      "51.552\n",
      "52.093\n",
      "52.093\n",
      "52.093\n",
      "52.112\n",
      "52.112\n",
      "52.112\n",
      "52.112\n",
      "52.112\n",
      "52.112\n",
      "52.111\n",
      "52.111\n",
      "52.111\n",
      "52.111\n",
      "52.11\n",
      "52.11\n",
      "52.11\n",
      "52.112\n",
      "52.112\n",
      "52.111\n",
      "52.11\n",
      "52.11\n",
      "52.11\n",
      "52.109\n",
      "52.109\n",
      "52.109\n",
      "52.11\n",
      "52.103\n",
      "52.103\n",
      "52.104\n",
      "52.106\n",
      "52.107\n",
      "52.107\n",
      "52.107\n",
      "52.11\n",
      "52.105\n",
      "52.105\n",
      "52.109\n",
      "52.697\n",
      "52.697\n",
      "52.697\n",
      "52.697\n",
      "52.697\n",
      "52.696\n",
      "52.696\n",
      "52.696\n",
      "52.696\n",
      "52.696\n",
      "52.696\n",
      "52.696\n",
      "52.711\n",
      "52.711\n",
      "52.711\n",
      "52.712\n",
      "52.712\n",
      "52.712\n",
      "52.712\n",
      "52.712\n",
      "52.71\n",
      "52.709\n",
      "52.708\n",
      "52.707\n",
      "52.71\n",
      "52.709\n",
      "52.712\n",
      "52.711\n",
      "52.709\n",
      "52.708\n",
      "52.708\n",
      "52.707\n",
      "52.707\n",
      "52.707\n",
      "52.707\n",
      "52.707\n",
      "52.707\n",
      "52.707\n",
      "52.707\n",
      "52.707\n",
      "52.707\n",
      "52.707\n",
      "52.707\n",
      "52.707\n",
      "52.707\n",
      "52.707\n",
      "52.709\n",
      "52.706\n",
      "53.023\n",
      "53.023\n",
      "53.023\n",
      "53.023\n",
      "53.023\n",
      "53.022\n",
      "53.052\n",
      "53.052\n",
      "53.079\n",
      "53.079\n",
      "53.106\n",
      "53.106\n",
      "53.106\n",
      "53.106\n",
      "53.106\n",
      "53.107\n",
      "53.107\n",
      "53.106\n",
      "53.106\n",
      "53.106\n",
      "53.106\n",
      "53.105\n",
      "53.104\n",
      "53.104\n",
      "53.104\n",
      "53.104\n",
      "53.104\n",
      "53.104\n",
      "53.104\n",
      "53.104\n",
      "53.103\n",
      "53.103\n",
      "53.104\n",
      "53.104\n",
      "53.104\n",
      "53.103\n",
      "53.103\n",
      "53.103\n",
      "53.515\n",
      "53.515\n",
      "53.515\n",
      "53.515\n",
      "53.515\n",
      "53.515\n",
      "53.542\n",
      "53.542\n",
      "53.572\n",
      "53.572\n",
      "53.599\n",
      "53.599\n",
      "53.627\n",
      "53.627\n",
      "53.627\n",
      "53.627\n",
      "53.627\n",
      "53.627\n",
      "53.627\n",
      "53.626\n",
      "53.626\n",
      "53.623\n",
      "53.623\n",
      "53.623\n",
      "53.624\n",
      "53.624\n",
      "53.626\n",
      "53.626\n",
      "53.624\n",
      "53.623\n",
      "53.623\n",
      "53.623\n",
      "53.623\n",
      "53.623\n",
      "53.622\n",
      "53.622\n",
      "53.622\n",
      "53.622\n",
      "53.622\n",
      "53.621\n",
      "53.621\n",
      "53.621\n",
      "53.621\n",
      "53.621\n",
      "53.62\n",
      "53.62\n",
      "53.62\n",
      "53.62\n",
      "53.621\n",
      "53.621\n",
      "53.621\n",
      "53.62\n",
      "53.62\n",
      "53.62\n",
      "54.521\n",
      "54.521\n",
      "54.52\n",
      "54.53\n",
      "54.53\n",
      "54.53\n",
      "54.53\n",
      "54.53\n",
      "54.53\n",
      "54.53\n",
      "54.529\n",
      "54.529\n",
      "54.529\n",
      "54.529\n",
      "54.529\n",
      "54.529\n",
      "54.528\n",
      "54.528\n",
      "54.528\n",
      "54.528\n",
      "54.528\n",
      "54.528\n",
      "54.525\n",
      "54.525\n",
      "54.525\n",
      "54.526\n",
      "54.528\n",
      "54.528\n",
      "54.526\n",
      "54.526\n",
      "54.525\n",
      "54.525\n",
      "54.524\n",
      "54.524\n",
      "54.524\n",
      "54.524\n",
      "54.524\n",
      "54.524\n",
      "54.524\n",
      "54.524\n",
      "54.524\n",
      "54.524\n",
      "54.524\n",
      "54.523\n",
      "54.523\n",
      "54.523\n",
      "54.523\n",
      "54.523\n",
      "54.523\n",
      "55.118\n",
      "55.118\n",
      "55.118\n",
      "55.118\n",
      "55.117\n",
      "55.117\n",
      "55.117\n",
      "55.117\n",
      "55.117\n",
      "55.117\n",
      "55.117\n",
      "55.116\n",
      "55.116\n",
      "55.116\n",
      "55.116\n",
      "55.115\n",
      "55.117\n",
      "55.117\n",
      "55.117\n",
      "55.116\n",
      "55.115\n",
      "55.114\n",
      "55.114\n",
      "55.114\n",
      "55.114\n",
      "55.114\n",
      "55.114\n",
      "55.114\n",
      "55.113\n",
      "55.113\n",
      "55.113\n",
      "55.113\n",
      "55.113\n",
      "55.113\n",
      "55.113\n",
      "55.114\n",
      "55.114\n",
      "55.114\n",
      "55.113\n",
      "55.115\n",
      "55.731\n",
      "55.731\n",
      "55.731\n",
      "55.731\n",
      "55.731\n",
      "55.73\n",
      "55.73\n",
      "55.729\n",
      "55.729\n",
      "55.729\n",
      "55.729\n",
      "55.729\n",
      "55.729\n",
      "55.728\n",
      "55.728\n",
      "55.74\n",
      "55.739\n",
      "55.741\n",
      "55.741\n",
      "55.741\n",
      "55.74\n",
      "55.74\n",
      "55.739\n",
      "55.738\n",
      "55.737\n",
      "55.738\n",
      "55.739\n",
      "55.74\n",
      "55.739\n",
      "55.739\n",
      "55.737\n",
      "55.737\n",
      "55.737\n",
      "55.736\n",
      "55.736\n",
      "55.736\n",
      "55.736\n",
      "55.736\n",
      "55.736\n",
      "55.736\n",
      "55.736\n",
      "55.736\n",
      "55.736\n",
      "55.735\n",
      "55.735\n",
      "55.734\n",
      "55.734\n",
      "55.734\n",
      "55.734\n",
      "55.734\n",
      "55.734\n",
      "55.734\n",
      "55.734\n",
      "55.735\n",
      "55.735\n",
      "55.735\n",
      "55.734\n",
      "55.734\n",
      "55.734\n",
      "55.734\n",
      "55.732\n",
      "55.73\n",
      "55.778\n",
      "55.838\n",
      "55.838\n",
      "55.838\n",
      "55.837\n",
      "55.837\n",
      "55.837\n",
      "55.837\n",
      "55.837\n",
      "55.855\n",
      "55.855\n",
      "55.855\n",
      "55.855\n",
      "55.855\n",
      "55.855\n",
      "55.855\n",
      "55.855\n",
      "55.855\n",
      "55.856\n",
      "55.854\n",
      "55.854\n",
      "55.854\n",
      "55.853\n",
      "55.853\n",
      "55.853\n",
      "55.852\n",
      "55.852\n",
      "55.852\n",
      "55.852\n",
      "55.852\n",
      "55.851\n",
      "55.851\n",
      "55.851\n",
      "55.851\n",
      "55.851\n",
      "55.851\n",
      "55.851\n",
      "55.851\n",
      "55.851\n",
      "55.85\n",
      "55.85\n",
      "55.85\n",
      "55.85\n",
      "55.85\n",
      "55.85\n",
      "55.85\n",
      "55.85\n",
      "55.85\n",
      "55.85\n",
      "55.85\n",
      "55.849\n",
      "55.848\n",
      "55.848\n",
      "55.848\n",
      "55.847\n",
      "55.847\n",
      "55.892\n",
      "55.892\n",
      "55.939\n",
      "55.985\n",
      "55.985\n",
      "56.032\n",
      "56.032\n",
      "56.032\n",
      "56.032\n",
      "56.032\n",
      "56.032\n",
      "56.061\n",
      "56.061\n",
      "56.062\n",
      "56.062\n",
      "56.062\n",
      "56.062\n",
      "56.062\n",
      "56.062\n",
      "56.061\n",
      "56.061\n",
      "56.06\n",
      "56.062\n",
      "56.062\n",
      "56.061\n",
      "56.061\n",
      "56.06\n",
      "56.06\n",
      "56.06\n",
      "56.059\n",
      "56.059\n",
      "56.059\n",
      "56.059\n",
      "56.058\n",
      "56.058\n",
      "56.058\n",
      "56.058\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.057\n",
      "56.056\n",
      "56.098\n",
      "56.098\n",
      "56.161\n",
      "56.161\n",
      "56.16\n",
      "56.16\n",
      "56.16\n",
      "56.179\n",
      "56.179\n",
      "56.179\n",
      "56.179\n",
      "56.178\n",
      "56.178\n",
      "56.178\n",
      "56.178\n",
      "56.178\n",
      "56.177\n",
      "56.177\n",
      "56.178\n",
      "56.177\n",
      "56.177\n",
      "56.176\n",
      "56.176\n",
      "56.176\n",
      "56.176\n",
      "56.176\n",
      "56.175\n",
      "56.175\n",
      "56.175\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "56.175\n",
      "56.175\n",
      "56.174\n",
      "56.174\n",
      "56.174\n",
      "56.173\n",
      "56.173\n",
      "56.173\n",
      "56.173\n",
      "56.173\n",
      "56.173\n",
      "56.172\n",
      "56.172\n",
      "56.172\n",
      "56.172\n",
      "56.172\n",
      "56.172\n",
      "56.172\n",
      "56.172\n",
      "56.172\n",
      "56.172\n",
      "56.172\n",
      "56.172\n",
      "56.172\n",
      "56.172\n",
      "56.171\n",
      "56.171\n",
      "56.216\n",
      "56.216\n",
      "56.261\n",
      "56.308\n",
      "56.308\n",
      "56.352\n",
      "56.351\n",
      "56.351\n",
      "56.351\n",
      "56.351\n",
      "56.351\n",
      "56.378\n",
      "56.378\n",
      "56.377\n",
      "56.377\n",
      "56.377\n",
      "56.377\n",
      "56.377\n",
      "56.377\n",
      "56.377\n",
      "56.377\n",
      "56.378\n",
      "56.378\n",
      "56.377\n",
      "56.377\n",
      "56.376\n",
      "56.376\n",
      "56.375\n",
      "56.375\n",
      "56.375\n",
      "56.375\n",
      "56.375\n",
      "56.375\n",
      "56.375\n",
      "56.374\n",
      "56.374\n",
      "56.374\n",
      "56.374\n",
      "56.373\n",
      "56.373\n",
      "56.373\n",
      "56.373\n",
      "56.373\n",
      "56.373\n",
      "56.372\n",
      "56.372\n",
      "56.372\n",
      "56.372\n",
      "56.372\n",
      "56.372\n",
      "56.372\n",
      "56.372\n",
      "56.372\n",
      "56.372\n",
      "56.372\n",
      "56.372\n",
      "56.372\n",
      "56.371\n",
      "56.529\n",
      "56.529\n",
      "56.529\n",
      "56.529\n",
      "56.529\n",
      "56.529\n",
      "56.528\n",
      "56.552\n",
      "56.552\n",
      "56.552\n",
      "56.552\n",
      "56.552\n",
      "56.552\n",
      "56.552\n",
      "56.552\n",
      "56.552\n",
      "56.552\n",
      "56.551\n",
      "56.551\n",
      "56.55\n",
      "56.55\n",
      "56.55\n",
      "56.55\n",
      "56.55\n",
      "56.55\n",
      "56.55\n",
      "56.55\n",
      "56.549\n",
      "56.549\n",
      "56.549\n",
      "56.549\n",
      "56.548\n",
      "56.548\n",
      "56.548\n",
      "56.548\n",
      "56.547\n",
      "56.547\n",
      "56.547\n",
      "56.547\n",
      "56.546\n",
      "56.546\n",
      "56.546\n",
      "56.546\n",
      "56.546\n",
      "56.546\n",
      "56.546\n",
      "56.546\n",
      "56.547\n",
      "56.547\n",
      "56.547\n",
      "56.547\n",
      "56.546\n",
      "56.546\n",
      "56.545\n",
      "56.544\n",
      "56.589\n",
      "56.589\n",
      "56.635\n",
      "56.681\n",
      "56.681\n",
      "56.727\n",
      "56.773\n",
      "56.795\n",
      "56.795\n",
      "56.795\n",
      "56.795\n",
      "56.795\n",
      "56.795\n",
      "56.813\n",
      "56.813\n",
      "56.813\n",
      "56.813\n",
      "56.813\n",
      "56.813\n",
      "56.813\n",
      "56.813\n",
      "56.812\n",
      "56.808\n",
      "56.812\n",
      "56.812\n",
      "56.811\n",
      "56.811\n",
      "56.811\n",
      "56.809\n",
      "56.809\n",
      "56.808\n",
      "56.81\n",
      "56.809\n",
      "56.811\n",
      "56.81\n",
      "56.809\n",
      "56.809\n",
      "56.809\n",
      "57.207\n",
      "57.207\n",
      "57.237\n",
      "57.237\n",
      "57.237\n",
      "57.236\n",
      "57.237\n",
      "57.237\n",
      "57.236\n",
      "57.235\n",
      "57.237\n",
      "57.597\n",
      "57.597\n",
      "57.621\n",
      "57.63\n",
      "57.63\n",
      "57.631\n",
      "57.631\n",
      "57.63\n",
      "57.63\n",
      "57.63\n",
      "57.63\n",
      "57.954\n",
      "57.954\n",
      "57.982\n",
      "57.983\n",
      "57.983\n",
      "57.984\n",
      "57.983\n",
      "57.983\n",
      "57.983\n",
      "57.983\n",
      "58.214\n",
      "58.216\n",
      "58.241\n",
      "58.241\n",
      "58.241\n",
      "58.242\n",
      "58.242\n",
      "58.242\n",
      "58.24\n",
      "58.24\n",
      "58.569\n",
      "58.569\n",
      "58.597\n",
      "58.597\n",
      "58.597\n",
      "58.598\n",
      "58.598\n",
      "58.598\n",
      "58.596\n",
      "58.596\n",
      "58.597\n",
      "58.607\n",
      "58.607\n",
      "58.622\n",
      "58.622\n",
      "58.622\n",
      "58.618\n",
      "58.622\n",
      "58.622\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.621\n",
      "58.62\n",
      "58.62\n",
      "58.62\n",
      "58.62\n",
      "58.62\n",
      "58.62\n",
      "58.62\n",
      "58.62\n",
      "58.62\n",
      "58.62\n",
      "58.619\n",
      "58.619\n",
      "58.619\n",
      "58.619\n",
      "58.619\n",
      "58.619\n",
      "58.619\n",
      "58.618\n",
      "58.618\n",
      "58.618\n",
      "58.618\n",
      "58.618\n",
      "58.618\n",
      "58.618\n",
      "58.617\n",
      "58.617\n",
      "58.617\n",
      "58.617\n",
      "58.616\n",
      "58.616\n",
      "58.616\n",
      "58.616\n",
      "58.616\n",
      "58.616\n",
      "58.616\n",
      "58.616\n",
      "58.615\n",
      "58.615\n",
      "58.615\n",
      "58.615\n",
      "58.615\n",
      "58.615\n",
      "58.615\n"
     ]
    },
    {
     "ename": "KeyboardInterrupt",
     "evalue": "",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mKeyboardInterrupt\u001b[0m                         Traceback (most recent call last)",
      "Cell \u001b[1;32mIn[4], line 17\u001b[0m\n\u001b[0;32m     14\u001b[0m \u001b[38;5;28mprint\u001b[39m(curr_val)\n\u001b[0;32m     15\u001b[0m mass\u001b[38;5;241m.\u001b[39mappend(curr_val)\n\u001b[1;32m---> 17\u001b[0m \u001b[43mtime\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43msleep\u001b[49m\u001b[43m(\u001b[49m\u001b[38;5;241;43m1\u001b[39;49m\u001b[43m)\u001b[49m\n\u001b[0;32m     18\u001b[0m i \u001b[38;5;241m+\u001b[39m\u001b[38;5;241m=\u001b[39m \u001b[38;5;241m1\u001b[39m\n",
      "\u001b[1;31mKeyboardInterrupt\u001b[0m: "
     ]
    }
   ],
   "source": [
    "x = [] #store time here \n",
    "mass = [] #store mass here\n",
    "i = 0\n",
    "\n",
    "while True:\n",
    "    x.append(i*2)\n",
    "    \n",
    "    status, val = get_mass(port)\n",
    "    if status == True:\n",
    "        curr_val = val\n",
    "    else:\n",
    "        curr_val = mass[-1]\n",
    "        \n",
    "    print(curr_val)\n",
    "    mass.append(curr_val)\n",
    "    \n",
    "    time.sleep(1)\n",
    "    i += 1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "d91fbc1e",
   "metadata": {
    "pycharm": {
     "name": "#%%\n"
    },
    "scrolled": false
   },
   "outputs": [],
   "source": [
    "results = {'Time': x, 'Mass': mass}"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e51054ed",
   "metadata": {},
   "source": [
    "### 4. Plotting and Exporting Data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "b7f2fb0f",
   "metadata": {
    "pycharm": {
     "name": "#%%\n"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Time</th>\n",
       "      <th>Mass</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>0.000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>0.000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>4</td>\n",
       "      <td>0.000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>6</td>\n",
       "      <td>0.000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>8</td>\n",
       "      <td>0.000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1505</th>\n",
       "      <td>3010</td>\n",
       "      <td>58.615</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1506</th>\n",
       "      <td>3012</td>\n",
       "      <td>58.615</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1507</th>\n",
       "      <td>3014</td>\n",
       "      <td>58.615</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1508</th>\n",
       "      <td>3016</td>\n",
       "      <td>58.615</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1509</th>\n",
       "      <td>3018</td>\n",
       "      <td>58.615</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>1510 rows × 2 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      Time    Mass\n",
       "0        0   0.000\n",
       "1        2   0.000\n",
       "2        4   0.000\n",
       "3        6   0.000\n",
       "4        8   0.000\n",
       "...    ...     ...\n",
       "1505  3010  58.615\n",
       "1506  3012  58.615\n",
       "1507  3014  58.615\n",
       "1508  3016  58.615\n",
       "1509  3018  58.615\n",
       "\n",
       "[1510 rows x 2 columns]"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mass_df = pd.DataFrame(results)\n",
    "mass_df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "78fb75a2",
   "metadata": {
    "pycharm": {
     "name": "#%%\n"
    },
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "/* global mpl */\n",
       "window.mpl = {};\n",
       "\n",
       "mpl.get_websocket_type = function () {\n",
       "    if (typeof WebSocket !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof MozWebSocket !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert(\n",
       "            'Your browser does not have WebSocket support. ' +\n",
       "                'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "                'Firefox 4 and 5 are also supported but you ' +\n",
       "                'have to enable WebSockets in about:config.'\n",
       "        );\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure = function (figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = this.ws.binaryType !== undefined;\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById('mpl-warnings');\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent =\n",
       "                'This browser does not support binary websocket messages. ' +\n",
       "                'Performance may be slow.';\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = document.createElement('div');\n",
       "    this.root.setAttribute('style', 'display: inline-block');\n",
       "    this._root_extra_style(this.root);\n",
       "\n",
       "    parent_element.appendChild(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen = function () {\n",
       "        fig.send_message('supports_binary', { value: fig.supports_binary });\n",
       "        fig.send_message('send_image_mode', {});\n",
       "        if (fig.ratio !== 1) {\n",
       "            fig.send_message('set_device_pixel_ratio', {\n",
       "                device_pixel_ratio: fig.ratio,\n",
       "            });\n",
       "        }\n",
       "        fig.send_message('refresh', {});\n",
       "    };\n",
       "\n",
       "    this.imageObj.onload = function () {\n",
       "        if (fig.image_mode === 'full') {\n",
       "            // Full images could contain transparency (where diff images\n",
       "            // almost always do), so we need to clear the canvas so that\n",
       "            // there is no ghosting.\n",
       "            fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "        }\n",
       "        fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "    };\n",
       "\n",
       "    this.imageObj.onunload = function () {\n",
       "        fig.ws.close();\n",
       "    };\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_header = function () {\n",
       "    var titlebar = document.createElement('div');\n",
       "    titlebar.classList =\n",
       "        'ui-dialog-titlebar ui-widget-header ui-corner-all ui-helper-clearfix';\n",
       "    var titletext = document.createElement('div');\n",
       "    titletext.classList = 'ui-dialog-title';\n",
       "    titletext.setAttribute(\n",
       "        'style',\n",
       "        'width: 100%; text-align: center; padding: 3px;'\n",
       "    );\n",
       "    titlebar.appendChild(titletext);\n",
       "    this.root.appendChild(titlebar);\n",
       "    this.header = titletext;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = (this.canvas_div = document.createElement('div'));\n",
       "    canvas_div.setAttribute('tabindex', '0');\n",
       "    canvas_div.setAttribute(\n",
       "        'style',\n",
       "        'border: 1px solid #ddd;' +\n",
       "            'box-sizing: content-box;' +\n",
       "            'clear: both;' +\n",
       "            'min-height: 1px;' +\n",
       "            'min-width: 1px;' +\n",
       "            'outline: 0;' +\n",
       "            'overflow: hidden;' +\n",
       "            'position: relative;' +\n",
       "            'resize: both;' +\n",
       "            'z-index: 2;'\n",
       "    );\n",
       "\n",
       "    function on_keyboard_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.key_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'keydown',\n",
       "        on_keyboard_event_closure('key_press')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'keyup',\n",
       "        on_keyboard_event_closure('key_release')\n",
       "    );\n",
       "\n",
       "    this._canvas_extra_style(canvas_div);\n",
       "    this.root.appendChild(canvas_div);\n",
       "\n",
       "    var canvas = (this.canvas = document.createElement('canvas'));\n",
       "    canvas.classList.add('mpl-canvas');\n",
       "    canvas.setAttribute(\n",
       "        'style',\n",
       "        'box-sizing: content-box;' +\n",
       "            'pointer-events: none;' +\n",
       "            'position: relative;' +\n",
       "            'z-index: 0;'\n",
       "    );\n",
       "\n",
       "    this.context = canvas.getContext('2d');\n",
       "\n",
       "    var backingStore =\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        this.context.webkitBackingStorePixelRatio ||\n",
       "        this.context.mozBackingStorePixelRatio ||\n",
       "        this.context.msBackingStorePixelRatio ||\n",
       "        this.context.oBackingStorePixelRatio ||\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        1;\n",
       "\n",
       "    this.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband_canvas = (this.rubberband_canvas = document.createElement(\n",
       "        'canvas'\n",
       "    ));\n",
       "    rubberband_canvas.setAttribute(\n",
       "        'style',\n",
       "        'box-sizing: content-box;' +\n",
       "            'left: 0;' +\n",
       "            'pointer-events: none;' +\n",
       "            'position: absolute;' +\n",
       "            'top: 0;' +\n",
       "            'z-index: 1;'\n",
       "    );\n",
       "\n",
       "    // Apply a ponyfill if ResizeObserver is not implemented by browser.\n",
       "    if (this.ResizeObserver === undefined) {\n",
       "        if (window.ResizeObserver !== undefined) {\n",
       "            this.ResizeObserver = window.ResizeObserver;\n",
       "        } else {\n",
       "            var obs = _JSXTOOLS_RESIZE_OBSERVER({});\n",
       "            this.ResizeObserver = obs.ResizeObserver;\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.resizeObserverInstance = new this.ResizeObserver(function (entries) {\n",
       "        var nentries = entries.length;\n",
       "        for (var i = 0; i < nentries; i++) {\n",
       "            var entry = entries[i];\n",
       "            var width, height;\n",
       "            if (entry.contentBoxSize) {\n",
       "                if (entry.contentBoxSize instanceof Array) {\n",
       "                    // Chrome 84 implements new version of spec.\n",
       "                    width = entry.contentBoxSize[0].inlineSize;\n",
       "                    height = entry.contentBoxSize[0].blockSize;\n",
       "                } else {\n",
       "                    // Firefox implements old version of spec.\n",
       "                    width = entry.contentBoxSize.inlineSize;\n",
       "                    height = entry.contentBoxSize.blockSize;\n",
       "                }\n",
       "            } else {\n",
       "                // Chrome <84 implements even older version of spec.\n",
       "                width = entry.contentRect.width;\n",
       "                height = entry.contentRect.height;\n",
       "            }\n",
       "\n",
       "            // Keep the size of the canvas and rubber band canvas in sync with\n",
       "            // the canvas container.\n",
       "            if (entry.devicePixelContentBoxSize) {\n",
       "                // Chrome 84 implements new version of spec.\n",
       "                canvas.setAttribute(\n",
       "                    'width',\n",
       "                    entry.devicePixelContentBoxSize[0].inlineSize\n",
       "                );\n",
       "                canvas.setAttribute(\n",
       "                    'height',\n",
       "                    entry.devicePixelContentBoxSize[0].blockSize\n",
       "                );\n",
       "            } else {\n",
       "                canvas.setAttribute('width', width * fig.ratio);\n",
       "                canvas.setAttribute('height', height * fig.ratio);\n",
       "            }\n",
       "            /* This rescales the canvas back to display pixels, so that it\n",
       "             * appears correct on HiDPI screens. */\n",
       "            canvas.style.width = width + 'px';\n",
       "            canvas.style.height = height + 'px';\n",
       "\n",
       "            rubberband_canvas.setAttribute('width', width);\n",
       "            rubberband_canvas.setAttribute('height', height);\n",
       "\n",
       "            // And update the size in Python. We ignore the initial 0/0 size\n",
       "            // that occurs as the element is placed into the DOM, which should\n",
       "            // otherwise not happen due to the minimum size styling.\n",
       "            if (fig.ws.readyState == 1 && width != 0 && height != 0) {\n",
       "                fig.request_resize(width, height);\n",
       "            }\n",
       "        }\n",
       "    });\n",
       "    this.resizeObserverInstance.observe(canvas_div);\n",
       "\n",
       "    function on_mouse_event_closure(name) {\n",
       "        /* User Agent sniffing is bad, but WebKit is busted:\n",
       "         * https://bugs.webkit.org/show_bug.cgi?id=144526\n",
       "         * https://bugs.webkit.org/show_bug.cgi?id=181818\n",
       "         * The worst that happens here is that they get an extra browser\n",
       "         * selection when dragging, if this check fails to catch them.\n",
       "         */\n",
       "        var UA = navigator.userAgent;\n",
       "        var isWebKit = /AppleWebKit/.test(UA) && !/Chrome/.test(UA);\n",
       "        if(isWebKit) {\n",
       "            return function (event) {\n",
       "                /* This prevents the web browser from automatically changing to\n",
       "                 * the text insertion cursor when the button is pressed. We\n",
       "                 * want to control all of the cursor setting manually through\n",
       "                 * the 'cursor' event from matplotlib */\n",
       "                event.preventDefault()\n",
       "                return fig.mouse_event(event, name);\n",
       "            };\n",
       "        } else {\n",
       "            return function (event) {\n",
       "                return fig.mouse_event(event, name);\n",
       "            };\n",
       "        }\n",
       "    }\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'mousedown',\n",
       "        on_mouse_event_closure('button_press')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'mouseup',\n",
       "        on_mouse_event_closure('button_release')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'dblclick',\n",
       "        on_mouse_event_closure('dblclick')\n",
       "    );\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    canvas_div.addEventListener(\n",
       "        'mousemove',\n",
       "        on_mouse_event_closure('motion_notify')\n",
       "    );\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'mouseenter',\n",
       "        on_mouse_event_closure('figure_enter')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'mouseleave',\n",
       "        on_mouse_event_closure('figure_leave')\n",
       "    );\n",
       "\n",
       "    canvas_div.addEventListener('wheel', function (event) {\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        on_mouse_event_closure('scroll')(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.appendChild(canvas);\n",
       "    canvas_div.appendChild(rubberband_canvas);\n",
       "\n",
       "    this.rubberband_context = rubberband_canvas.getContext('2d');\n",
       "    this.rubberband_context.strokeStyle = '#000000';\n",
       "\n",
       "    this._resize_canvas = function (width, height, forward) {\n",
       "        if (forward) {\n",
       "            canvas_div.style.width = width + 'px';\n",
       "            canvas_div.style.height = height + 'px';\n",
       "        }\n",
       "    };\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    canvas_div.addEventListener('contextmenu', function (_e) {\n",
       "        event.preventDefault();\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus() {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'mpl-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'mpl-button-group';\n",
       "    for (var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'mpl-button-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        var button = (fig.buttons[name] = document.createElement('button'));\n",
       "        button.classList = 'mpl-widget';\n",
       "        button.setAttribute('role', 'button');\n",
       "        button.setAttribute('aria-disabled', 'false');\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "\n",
       "        var icon_img = document.createElement('img');\n",
       "        icon_img.src = '_images/' + image + '.png';\n",
       "        icon_img.srcset = '_images/' + image + '_large.png 2x';\n",
       "        icon_img.alt = tooltip;\n",
       "        button.appendChild(icon_img);\n",
       "\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    var fmt_picker = document.createElement('select');\n",
       "    fmt_picker.classList = 'mpl-widget';\n",
       "    toolbar.appendChild(fmt_picker);\n",
       "    this.format_dropdown = fmt_picker;\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = document.createElement('option');\n",
       "        option.selected = fmt === mpl.default_extension;\n",
       "        option.innerHTML = fmt;\n",
       "        fmt_picker.appendChild(option);\n",
       "    }\n",
       "\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.request_resize = function (x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', { width: x_pixels, height: y_pixels });\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.send_message = function (type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function () {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({ type: 'draw', figure_id: this.id }));\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function (fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] !== fig.canvas.width || size[1] !== fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1], msg['forward']);\n",
       "        fig.send_message('refresh', {});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function (fig, msg) {\n",
       "    var x0 = msg['x0'] / fig.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / fig.ratio;\n",
       "    var x1 = msg['x1'] / fig.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / fig.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0,\n",
       "        0,\n",
       "        fig.canvas.width / fig.ratio,\n",
       "        fig.canvas.height / fig.ratio\n",
       "    );\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function (fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function (fig, msg) {\n",
       "    fig.canvas_div.style.cursor = msg['cursor'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_message = function (fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function (fig, _msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function (fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_history_buttons = function (fig, msg) {\n",
       "    for (var key in msg) {\n",
       "        if (!(key in fig.buttons)) {\n",
       "            continue;\n",
       "        }\n",
       "        fig.buttons[key].disabled = !msg[key];\n",
       "        fig.buttons[key].setAttribute('aria-disabled', !msg[key]);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_navigate_mode = function (fig, msg) {\n",
       "    if (msg['mode'] === 'PAN') {\n",
       "        fig.buttons['Pan'].classList.add('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    } else if (msg['mode'] === 'ZOOM') {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.add('active');\n",
       "    } else {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function () {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message('ack', {});\n",
       "};\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function (fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            var img = evt.data;\n",
       "            if (img.type !== 'image/png') {\n",
       "                /* FIXME: We get \"Resource interpreted as Image but\n",
       "                 * transferred with MIME type text/plain:\" errors on\n",
       "                 * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "                 * to be part of the websocket stream */\n",
       "                img.type = 'image/png';\n",
       "            }\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src\n",
       "                );\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                img\n",
       "            );\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        } else if (\n",
       "            typeof evt.data === 'string' &&\n",
       "            evt.data.slice(0, 21) === 'data:image/png;base64'\n",
       "        ) {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig['handle_' + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\n",
       "                \"No handler for the '\" + msg_type + \"' message type: \",\n",
       "                msg\n",
       "            );\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\n",
       "                    \"Exception inside the 'handler_\" + msg_type + \"' callback:\",\n",
       "                    e,\n",
       "                    e.stack,\n",
       "                    msg\n",
       "                );\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "};\n",
       "\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * https://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys(original) {\n",
       "    return Object.keys(original).reduce(function (obj, key) {\n",
       "        if (typeof original[key] !== 'object') {\n",
       "            obj[key] = original[key];\n",
       "        }\n",
       "        return obj;\n",
       "    }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function (event, name) {\n",
       "    if (name === 'button_press') {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    // from https://stackoverflow.com/q/1114465\n",
       "    var boundingRect = this.canvas.getBoundingClientRect();\n",
       "    var x = (event.clientX - boundingRect.left) * this.ratio;\n",
       "    var y = (event.clientY - boundingRect.top) * this.ratio;\n",
       "\n",
       "    this.send_message(name, {\n",
       "        x: x,\n",
       "        y: y,\n",
       "        button: event.button,\n",
       "        step: event.step,\n",
       "        guiEvent: simpleKeys(event),\n",
       "    });\n",
       "\n",
       "    return false;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (_event, _name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.key_event = function (event, name) {\n",
       "    // Prevent repeat events\n",
       "    if (name === 'key_press') {\n",
       "        if (event.key === this._key) {\n",
       "            return;\n",
       "        } else {\n",
       "            this._key = event.key;\n",
       "        }\n",
       "    }\n",
       "    if (name === 'key_release') {\n",
       "        this._key = null;\n",
       "    }\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.key !== 'Control') {\n",
       "        value += 'ctrl+';\n",
       "    }\n",
       "    else if (event.altKey && event.key !== 'Alt') {\n",
       "        value += 'alt+';\n",
       "    }\n",
       "    else if (event.shiftKey && event.key !== 'Shift') {\n",
       "        value += 'shift+';\n",
       "    }\n",
       "\n",
       "    value += 'k' + event.key;\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, { key: value, guiEvent: simpleKeys(event) });\n",
       "    return false;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function (name) {\n",
       "    if (name === 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message('toolbar_button', { name: name });\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function (tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "\n",
       "///////////////// REMAINING CONTENT GENERATED BY embed_js.py /////////////////\n",
       "// prettier-ignore\n",
       "var _JSXTOOLS_RESIZE_OBSERVER=function(A){var t,i=new WeakMap,n=new WeakMap,a=new WeakMap,r=new WeakMap,o=new Set;function s(e){if(!(this instanceof s))throw new TypeError(\"Constructor requires 'new' operator\");i.set(this,e)}function h(){throw new TypeError(\"Function is not a constructor\")}function c(e,t,i,n){e=0 in arguments?Number(arguments[0]):0,t=1 in arguments?Number(arguments[1]):0,i=2 in arguments?Number(arguments[2]):0,n=3 in arguments?Number(arguments[3]):0,this.right=(this.x=this.left=e)+(this.width=i),this.bottom=(this.y=this.top=t)+(this.height=n),Object.freeze(this)}function d(){t=requestAnimationFrame(d);var s=new WeakMap,p=new Set;o.forEach((function(t){r.get(t).forEach((function(i){var r=t instanceof window.SVGElement,o=a.get(t),d=r?0:parseFloat(o.paddingTop),f=r?0:parseFloat(o.paddingRight),l=r?0:parseFloat(o.paddingBottom),u=r?0:parseFloat(o.paddingLeft),g=r?0:parseFloat(o.borderTopWidth),m=r?0:parseFloat(o.borderRightWidth),w=r?0:parseFloat(o.borderBottomWidth),b=u+f,F=d+l,v=(r?0:parseFloat(o.borderLeftWidth))+m,W=g+w,y=r?0:t.offsetHeight-W-t.clientHeight,E=r?0:t.offsetWidth-v-t.clientWidth,R=b+v,z=F+W,M=r?t.width:parseFloat(o.width)-R-E,O=r?t.height:parseFloat(o.height)-z-y;if(n.has(t)){var k=n.get(t);if(k[0]===M&&k[1]===O)return}n.set(t,[M,O]);var S=Object.create(h.prototype);S.target=t,S.contentRect=new c(u,d,M,O),s.has(i)||(s.set(i,[]),p.add(i)),s.get(i).push(S)}))})),p.forEach((function(e){i.get(e).call(e,s.get(e),e)}))}return s.prototype.observe=function(i){if(i instanceof window.Element){r.has(i)||(r.set(i,new Set),o.add(i),a.set(i,window.getComputedStyle(i)));var n=r.get(i);n.has(this)||n.add(this),cancelAnimationFrame(t),t=requestAnimationFrame(d)}},s.prototype.unobserve=function(i){if(i instanceof window.Element&&r.has(i)){var n=r.get(i);n.has(this)&&(n.delete(this),n.size||(r.delete(i),o.delete(i))),n.size||r.delete(i),o.size||cancelAnimationFrame(t)}},A.DOMRectReadOnly=c,A.ResizeObserver=s,A.ResizeObserverEntry=h,A}; // eslint-disable-line\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis\", \"fa fa-square-o\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pgf\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\", \"webp\"];\n",
       "\n",
       "mpl.default_extension = \"png\";/* global mpl */\n",
       "\n",
       "var comm_websocket_adapter = function (comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.binaryType = comm.kernel.ws.binaryType;\n",
       "    ws.readyState = comm.kernel.ws.readyState;\n",
       "    function updateReadyState(_event) {\n",
       "        if (comm.kernel.ws) {\n",
       "            ws.readyState = comm.kernel.ws.readyState;\n",
       "        } else {\n",
       "            ws.readyState = 3; // Closed state.\n",
       "        }\n",
       "    }\n",
       "    comm.kernel.ws.addEventListener('open', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('close', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('error', updateReadyState);\n",
       "\n",
       "    ws.close = function () {\n",
       "        comm.close();\n",
       "    };\n",
       "    ws.send = function (m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function (msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        var data = msg['content']['data'];\n",
       "        if (data['blob'] !== undefined) {\n",
       "            data = {\n",
       "                data: new Blob(msg['buffers'], { type: data['blob'] }),\n",
       "            };\n",
       "        }\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(data);\n",
       "    });\n",
       "    return ws;\n",
       "};\n",
       "\n",
       "mpl.mpl_figure_comm = function (comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = document.getElementById(id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm);\n",
       "\n",
       "    function ondownload(figure, _format) {\n",
       "        window.open(figure.canvas.toDataURL());\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy, ondownload, element);\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element;\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error('Failed to find cell for figure', id, fig);\n",
       "        return;\n",
       "    }\n",
       "    fig.cell_info[0].output_area.element.on(\n",
       "        'cleared',\n",
       "        { fig: fig },\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function (fig, msg) {\n",
       "    var width = fig.canvas.width / fig.ratio;\n",
       "    fig.cell_info[0].output_area.element.off(\n",
       "        'cleared',\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "    fig.resizeObserverInstance.unobserve(fig.canvas_div);\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable();\n",
       "    fig.parent_element.innerHTML =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "    fig.close_ws(fig, msg);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.close_ws = function (fig, msg) {\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function (_remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width / this.ratio;\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function () {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message('ack', {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () {\n",
       "        fig.push_to_output();\n",
       "    }, 1000);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'btn-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'btn-group';\n",
       "    var button;\n",
       "    for (var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'btn-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        button = fig.buttons[name] = document.createElement('button');\n",
       "        button.classList = 'btn btn-default';\n",
       "        button.href = '#';\n",
       "        button.title = name;\n",
       "        button.innerHTML = '<i class=\"fa ' + image + ' fa-lg\"></i>';\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message pull-right';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = document.createElement('div');\n",
       "    buttongrp.classList = 'btn-group inline pull-right';\n",
       "    button = document.createElement('button');\n",
       "    button.classList = 'btn btn-mini btn-primary';\n",
       "    button.href = '#';\n",
       "    button.title = 'Stop Interaction';\n",
       "    button.innerHTML = '<i class=\"fa fa-power-off icon-remove icon-large\"></i>';\n",
       "    button.addEventListener('click', function (_evt) {\n",
       "        fig.handle_close(fig, {});\n",
       "    });\n",
       "    button.addEventListener(\n",
       "        'mouseover',\n",
       "        on_mouseover_closure('Stop Interaction')\n",
       "    );\n",
       "    buttongrp.appendChild(button);\n",
       "    var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n",
       "    titlebar.insertBefore(buttongrp, titlebar.firstChild);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._remove_fig_handler = function (event) {\n",
       "    var fig = event.data.fig;\n",
       "    if (event.target !== this) {\n",
       "        // Ignore bubbled events from children.\n",
       "        return;\n",
       "    }\n",
       "    fig.close_ws(fig, {});\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (el) {\n",
       "    el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (el) {\n",
       "    // this is important to make the div 'focusable\n",
       "    el.setAttribute('tabindex', 0);\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    } else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (event, _name) {\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which === 13) {\n",
       "        this.canvas_div.blur();\n",
       "        // select the cell after this one\n",
       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
       "        IPython.notebook.select(index + 1);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "};\n",
       "\n",
       "mpl.find_output_cell = function (html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i = 0; i < ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code') {\n",
       "            for (var j = 0; j < cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] === html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "};\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel !== null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target(\n",
       "        'matplotlib',\n",
       "        mpl.mpl_figure_comm\n",
       "    );\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAyAAAAJYCAYAAACadoJwAAAAAXNSR0IArs4c6QAAIABJREFUeF7s3XeYFFXahvFnZsg5IzmroGRESQICghEVE4prwhUDi7BmUEEJLvipiAoqrjnnRQXJQUBAkAySc84ZJn1XnWGGIU53n9M10zN37bIM01XP6f6d94/zbnVVRSUmJiaKDQEEEEAAAQQQQAABBBDwQSCKBsQHZYZAAAEEEEAAAQQQQAABI0ADQiEggAACCCCAAAIIIICAbwI0IL5RMxACCCCAAAIIIIAAAgjQgFADCCCAAAIIIIAAAggg4JsADYhv1AyEAAIIIIAAAggggAACNCDUAAIIIIAAAggggAACCPgmQAPiGzUDIYAAAggggAACCCCAAA0INYAAAggggAACCCCAAAK+CdCA+EbNQAgggAACCCCAAAIIIEADQg0ggAACCCCAAAIIIICAbwI0IL5RMxACCCCAAAIIIIAAAgjQgFADCCCAAAIIIIAAAggg4JsADYhv1AyEAAIIIIAAAggggAACNCDUAAIIIIAAAggggAACCPgmQAPiGzUDIYAAAggggAACCCCAAA0INYAAAggggAACCCCAAAK+CdCA+EbNQAgggAACCCCAAAIIIEADQg0ggAACCCCAAAIIIICAbwI0IL5RMxACCCCAAAIIIIAAAgjQgFADCCCAAAIIIIAAAggg4JsADYhv1AyEAAIIIIAAAggggAACNCDUAAIIIIAAAggggAACCPgmQAPiGzUDIYAAAggggAACCCCAAA0INYAAAggggAACCCCAAAK+CdCA+EbNQAgggAACCCCAAAIIIEADQg0ggAACCCCAAAIIIICAbwI0IL5RMxACCCCAAAIIIIAAAgjQgFADCCCAAAIIIIAAAggg4JsADYhv1AyEAAIIIIAAAggggAACNCDUAAIIIIAAAggggAACCPgmQAPiGzUDIYAAAggggAACCCCAAA0INYAAAggggAACCCCAAAK+CdCA+EbNQAgggAACCCCAAAIIIEADQg0ggAACCCCAAAIIIICAbwI0IL5RMxACCCCAAAIIIIAAAgjQgFADCCCAAAIIIIAAAggg4JsADYhv1AyEAAIIIIAAAggggAACNCDUAAIIIIAAAggggAACCPgmQAPiGzUDIYAAAggggAACCCCAAA0INYAAAggggAACCCCAAAK+CdCA+EbNQAgggAACCCCAAAIIIEADQg0ggAACCCCAAAIIIICAbwI0IL5RMxACCCCAAAIIIIAAAgjQgFADCCCAAAIIIIAAAggg4JsADYhv1AyEAAIIIIAAAggggAACNCDUAAIIIIAAAggggAACCPgmQAPiGzUDIYAAAggggAACCCCAAA0INYAAAggggAACCCCAAAK+CdCA+EbNQAgggAACCCCAAAIIIEADQg0ggAACCCCAAAIIIICAbwI0IL5RMxACCCCAAAIIIIAAAgjQgFADCCCAAAIIIIAAAggg4JsADYhv1AyEAAIIIIAAAggggAACNCDUAAIIIIAAAggggAACCPgmQAPiGzUDIYAAAggggAACCCCAAA0INYAAAggggAACCCCAAAK+CdCA+EbNQAgggAACCCCAAAIIIEADQg0ggAACCCCAAAIIIICAbwI0IL5RMxACCCCAAAIIIIAAAgjQgFADCCCAAAIIIIAAAggg4JsADYhv1AyEAAIIIIAAAggggAACNCDUAAIIIIAAAggggAACCPgmQAPiGzUDIYAAAggggAACCCCAAA0INYAAAggggAACCCCAAAK+CdCA+EbNQAgggAACCCCAAAIIIEADQg0ggAACCCCAAAIIIICAbwI0IL5RMxACCCCAAAIIIIAAAgjQgFADCCCAAAIIIIAAAggg4JsADYhv1AyEAAIIIIAAAggggAACNCDUAAIIIIAAAggggAACCPgmQAPiG3XaA+3YsUO//fabKlasqNy5c6d9AHsggAACCCCAAAII+Cpw+PBhrVmzRm3btlWxYsV8HTuzDEYDkoFm8rPPPlOnTp0y0DvirSCAAAIIIIAAAgicSeDTTz/VnXfeCU4IAjQgIaCF65CpU6eqadOm8gq6evXq4RqGXAQQQAABBBBAAIEQBZYsWWL+D+Pff/9dTZo0CTElax9GA5KB5n/OnDmqX7++Zs+erXr16mWgd8ZbQQABBBBAAAEEEPAEWK/Z1wENiL2hswQK2hklQQgggAACCCCAQFgEWK/Zs9KA2Bs6S6CgnVEShAACCCCAAAIIhEWA9Zo9Kw3IGQy3bNmil156ST///LO8n4sUKaIGDRrorbfeUvny5VOOeO+99zR48GCtWLFCxYsX1x133KHevXuHfAcrCtq+oElAAAEEEEAAAQTCKcB6zV6XBuQUw+XLl+vyyy9Xzpw5de+996pcuXLauXOnZsyYoeeee061a9c2RwwcOFBPPfWU2rdvr2uuuUaLFy/WkCFD1K5dO9O4hLJR0KGocQwCCCCAAAIIIOCfAOs1e2sakFSGiYmJuvTSSxUbG6vJkycrf/78ZxTevn27KlSooDZt2uinn35K2ee1115Tjx49TAPiNSXBbhR0sGLsjwACCCCAAAII+CvAes3emwYkleH48ePVqlUr/e9//9N1112nI0eOKDo6Wjly5DhJ2vvq1T//+U95+7ds2TLltUOHDqlo0aLq0KGDuZVusBsFHawY+yOAAAIIIIAAAv4KsF6z96YBSWX45JNPatCgQebsR8+ePTVlyhRFRUWZsyKvvvqqGjVqZPZ+8MEH9e6778prOE59Yrm37969e7V06dKgZ4eCDpqMAxBAAAEEEEAAAV8FWK/Zc9OApDK84YYbzFeqvAvKvWbDe7qld/1H3759tWvXLs2cOVM1a9Y0Z0emTZtmXjt1u+mmmzRmzBjt37//nLOzefNmeX9Sb8kPtuE5IPaFTQICCCCAAAIIIBAOARoQe1UakFSGrVu31rhx49S8eXNNnDgx5RXvbMbFF18sr7n4+uuvzde0vGZh06ZNp82Adycsb5+4uLhzzo53t6w+ffqccR8aEPvCJgEBBBBAAAEEEAiHAA2IvSoNSCpD78yGdwH58OHDdf/995+k26RJE3O73a1bt3IGxL7uSEAAAQQQQAABBCJSgAbEftpoQFIZdunSRe+8845GjBiha6+99iTdm2++2VycfuzYMa4Bsa87EhBAAAEEEEAAgYgUoAGxnzYakFSG77//vjp37qyhQ4fKa0ZSb97F5Rs3btSGDRvMBejehehnuwuW91Wtzz77LOjZoaCDJuMABBBAAAEEEEDAVwHWa/bcNCCpDHfs2GGe71G9enX98ccfypYtm3l11qxZ5k5Y9913n/l61rZt28x+bdu21Y8//piSkPwckOTb+AY7PRR0sGLsjwACCCCAAAII+CvAes3emwbkFMPBgwfrscceU+PGjXX77bfLa0q833nPAvEuDveejO5tL7/8sp555hl5d87yHjq4aNEivfnmm/IuZB85cmRIM0NBh8TGQQgggAACCCCAgG8CrNfsqWlAzmDoPUTQe+7H4sWLlSdPHvPE8wEDBqhy5con7e1dL+I1JytXrjS37u3YsaO5s5V3TCgbBR2KGscggAACCCCAAAL+CbBes7emAbE3dJZAQTujJAgBBBBAAAEEEAiLAOs1e1YaEHtDZwkUtDNKghBAAAEEEEDAZ4Ht+4/qwNE4xcUnKDY+UfuPxGr/kRPPRYuKSnpDKX8r+RfHf3/8/UYd3+H4q6ftf+L4VB8wSqpXvrByZY8J+6dmvWZPTANib+gsgYJ2RkkQAggggAACCPgksH7XId3/0Swt23rApxHPPMzUp69QmUK5w/4eWK/ZE9OA2Bs6S6CgnVEShAACCCCAAAI+Cbzy2996c8IKn0Y7+zA0IOk+BQG/ARqQgKnCvyMNSPiNGQEBBBBAAAEE7AX+3rJfv8zfpL2HYzVlxQ6t2n7QhD7R9gLlzBatbNFRypMjmwrnzWG+aJXo/Un0/jfpZ/N38g/Hf5P871NfTzzr6ykBJq9NjZJmzHBvrNfshWlA7A2dJVDQzigJQgABBBBAAIEwCWzac1gtXpmoY3EJJ41wcZkC+rlrszCNmnFiWa/ZzwUNiL2hswQK2hklQQgggAACCCDgUGD51v36ePpardl5UDsPHNPizftOSs+VPVp9b6ipm+uXdThqxoxivWY/LzQg9obOEihoZ5QEIYAAAggggIAjAe+rU80HTdS6XYdOS/zjmVYqlCe7YqKjlD0m2tGIGTuG9Zr9/NCA2Bs6S6CgnVEShAACCCCAAAKWAkdi47Vy+wHtPRSrO4bPSEkrXTCXYmKi1LFheT3coqrlKJF3OOs1+zmjAbE3dJZAQTujJAgBBBBAAAEELASe/2mhvpi5zjzPI/X28k01dXvD8hbJkX8o6zX7OaQBsTd0lkBBO6MkCAEEEEAAAQRCFNi674gu7T/utKOjo6SxPZqrcvF8ISZnjsNYr9nPIw2IvaGzBAraGSVBCCCAAAIIIBCEgHdHq8nLtmvXoWNat/NQynM9KhXLq/uaVjJJtcoUVO1yhYJIzZy7sl6zn1caEHtDZwkUtDNKghBAAAEEEEAgQAHvWo+2r0/W2p2nX2T+1T8v06WViwaYlDV2Y71mP880IPaGzhIoaGeUBCGAAAIIIIBAgAIzV+/Sre9MP23v8wrk0oTHWyh3jpgAk7LGbqzX7OeZBsTe0FkCBe2MkiAEEEAAAQSytEB8QqK8p5UfjYvX7kPHdOBovHkiubdFRUlR3n/M39Kfa3fr/d9Xm9cealFFrauXMD9XL1XAlyeLR9pEsV6znzEaEHtDZwkUtDNKghBAAAEEEMiyAn+u2aUnv52vVTsOBm0w+YmWKl80T9DHZaUDWK/ZzzYNiL2hswQK2hklQQgggAACCGR4gQNH4zRrzS7FxiVoz6FYbT9w9IxnKLwzFd6WfNbC/Bzl/SvpbEbSa1JcQqKmLN+hScu2h/TZKxfLqzE9mpuHCrKdXYD1mn110IDYGzpLoKCdURKEAAIIIIBAWARi4xMUF5+ozXsPm682pbQGx7/OdObmIOnrTslb8leg/vHfGdpxIDnD/dtteUFxXXnRecoRE60SBXKaBibR+0+ilPx0D+8p597P0VFRqle+kPLnyu7+jWSyRNZr9hNKA2Jv6CyBgnZGSRACCCCAAAJOBTbsPqT+vy7RuCXbdDQuwWm26zDv4vFnr6mu62uXdh1NniTWa/ZlQANib+gsgYJ2RkkQAggggEAWFvCeZ/HdnA06dCxe2/YdUXxi4slfX0o6TZHyFaakrzIlnaJI/lqTefX4WY39R+K0ePO+sIl6F343rFhEZQrnVs5s0SlnKJLPTngDe2ctvPMWyWcvkv5OOo9hfj7+75zZYlSleN6UzxO2N52Fg1mv2U8+DYi9obMECtoZJUEIIIAAAhEi4D2D4oOpazR91U4dPBqnNTsO6nBsfEozkNQrnGgOkl84cUen468d/8XR2ISU48NBUCxfDl1xYQkVzpND5YrkMV9vSv5aU1KLcKIZSGkOTnQQJ15P6ih0/nn51bhKsXC8VTLDJMB6zR6WBsTe0FkCBe2MkiAEEEAAgTAIrNi2Xws3Jp0J8C6g3ns4NmWU1Ld2NU3DKbd7Td4x9cXT2/cf1bezN2jb/qSLr8OxlS2cWwVyZT/pmocTjcGZr4fwdk5qJI6fYZCUP1c2XVurtO5vWknZY6LD8VbJjBAB1mv2E0UDYm/oLIGCdkZJEAIIIJAlBZZv3a8vZ62Xd1Yh0C31xdHnOmbltoPmLEU4tzKFcqtQnuyqWCyvuSg6dQNgxj1+5fRJZxzO8Lu8ObOpadViuqleGb6KFM4Jy6LZrNfsJ54GxN7QWQIF7YySIAQQQCBLCtz09lTNWbcn4j57u4vOU+/rL9J5BXNF3HvnDWc9AdZr9nNOA2Jv6CyBgnZGSRACCCCQJQVq9xl90teiXCN4d1fqdFl51a9QRN6jIkoWyGWeGZF8AfSJC6STL48++8XTx7/dJC+zcN4crt8qeQiETYD1mj0tDYi9obMECtoZJUEIIIBAlhM4Fpeg83uNNJ/79kvKqe8NFzs38L4WFc1D6py7EhhZAqzX7OeLBsTe0FkCBe2MkiAEEEAgYgW86x5W7TioQ0dPXMdxpus0TnqwnaK059Ax3TF8hvncj7asqsfbXhCxBrxxBDKyAOs1+9mhAbE3dJZAQTujJAgBBBCIWIGXRy7VsEkrrd7/89fW0H1NK1llcDACCJxZgPWafWXQgNgbOkugoJ1REoQAAghErECr/5uoldsPWr3/b7s0UoOKRawyOBgBBGhAwlUDNCDhkg0hlwYkBDQOQQABBDKZQPKF5DVKFdAdl5Y3n+74nWaTPmmqZ1Mkf/TkC7q9f1crkU+Nq/Jgu0xWFnycDCTAes1+MmhA7A2dJVDQzigJQgABBCJSIPWF5B0bltOAm2pF5OfgTSOQmQVYr9nPLg2IvaGzBAraGSVBCCCAgO8Cvf+3SP+bt0kJqU9HnDhhcdr7SX7IXuoXvEP3H40zv+JCct+nkAERCEiA9VpATOfciQbE3tBZAgXtjJIgBBBAwFeBbfuPqGG/cU7HHHRzLd3SoJzTTMIQQMBegPWavSENiL2hswQK2hklQQgggICvAsu27teVr002Y154Xn5VKJrnpPGjFHXyv0/+p069zW7Fonn1r1bVlCt7jK+fg8EQQCBtAdZraRultQcNSFpCPr5OQfuIzVAIIICAQ4HZa3epw9DpJnFgh1q69RLOXDjkJQqBDCXAes1+OmhA7A2dJVDQzigJQgABBHwVmPD3Nt37wSwz5tA76+mqmqV8HZ/BEEDAPwHWa/bWNCD2hs4SKGhnlAQhgAACvgp4F5//64u/zJif3n+pmlbjNri+TgCDIeCjAOs1e2waEHtDZwkUtDNKghBAAIGgBI7Exmvayh1auHGfvFvhetdkpFymEZV0BUfS76JOes38LipKCzbs1ahFW8yYPz7SRHXKFQpqfHZGAIHIEWC9Zj9XNCD2hs4SKGhnlAQhgAAC5xTwboH7xrgVmrNutw4cjdPc9XsUn3DS4/5CFhz37+aqUjxfyMdzIAIIZGwB1mv280MDYm/oLIGCdkZJEAIIIHBOgVlrdumWYUkXjbvcLi5TQD890lQx0afc5srlIGQhgEC6CrBes+enAbE3dJZAQTujJAgBBBA4p8DnM9bp2R8WmH2K5cuhSsXyqnGVYubaDe8WuIk6fjYkUeYn7wGB3u+S/k7avLMoyc8cTP67bOHciqb5oPoQyNQCrNfsp5cGxN7QWQIF7YySIAQQyEICuw4e0+a9h7Vt31FzHcf8DXt1JC7BCCSfh0h9TYd3zcaWvUe0cc9hs8+MZ1upZIFcWUiMj4oAAjYCrNds9JKOpQGxN3SWQEE7oyQIAQSyiMDstbt16zvTQ75+I0+OGC3s3ZazFlmkXviYCLgQYL1mr0gDYm/oLIGCdkZJEAIIZBGBfr8s1ntTVp/0aXNki1bFonnMHauSv0qV8lWp41+d8g7IHhOtexpX1O0Ny2cRLT4mAgi4EGC9Zq9IA2Jv6CyBgnZGSRACCGRwgf1HYtX9q7lavu3AyddXJF964V1fcdZrL7wPl3T9xf6jcea2uflzZtOTV12oC0rmV/0KhbkIPIPPP28PgUgWYL1mP3s0IPaGzhIoaGeUBCGAQAYX+HrWej353Xxn77JJ1aL6rPNlzvIIQgABBM4mwHrNvjZoQOwNnSVQ0M4oCUIAgQwu8PrYZXp97HLzLr0zFjliopMe8JfqYX/ea94F4yceApj0oVL/zrvMPH+ubHqweWVdeF6BDP6peXsIIJAZBFiv2c8iDYi9obMECtoZJUEIIJDBBZ77caE++WOteZdLX2qnXNljMvg75u0hgAACSQKs1+wrgQbE3tBZAgXtjJIgBBDwUWDi39s0edkObd135MSox+9/e+I2uCcezOf95N29yrsNbr6c2bSwT1sf3y1DIYAAAnYCrNfs/LyjaUDsDZ0lUNDOKAlCAAGfBNbvOqTmgyYoIfnpfEGOW7lYXo1/vEWQR7E7AgggkH4CrNfs7WlA7A2dJVDQzigJQgABnwTGL92q+z78M2U0725UqZ8U7r1w4t9JuyXfGjd39hg9fdWFuu0SboPr03QxDAIIOBBgvWaPSANib+gsgYJ2RkkQAgj4JPDVrHV66rsFZrSvH2ykhpWK+DQywyCAAALpI8B6zd6dBsTe0FkCBe2MkiAEEAhAwPv61PRVO7V17xEdiYvXlr1HtefQMXOkdzcq7w5TSXelOtPdqcyrWrXjoJZs3meOmfREC1UomjeAkdkFAQQQiFwB1mv2c0cDYm/oLIGCdkZJEAIIpCFw8GicGg0Yp31H4pxYeY3K4j7tlDsHd7NyAkoIAghkWAHWa/ZTQwNib+gsgYJ2RkkQAgikIbBw415dO+T30/bynqmRM5vXRCQ9aTzpaeSpnkp+/GfvheTXoqOj1OmyCnqq3YW4I4AAAplegPWa/RTTgNgbOkugoJ1REoQAAmkIeLfOveeDWWavR1pW0d2NKypnTIwK5smOHQIIIIDAOQRYr9mXBw2IvaGzBAraGSVBCCCQhsD3czaox9fzzF5fPHCZGlUpihkCCCCAQAACrNcCQEpjFxoQe0NnCRS0M0qCEMgyArPX7tKLPy8xF48nfWXq+FenvJ9Tf3Xq+O+TnteRqMPH4nXwWLxxGt39cp1fMn+WMeODIoAAAjYCrNds9JKOpQGxN3SWQEE7oyQIgSwjcP+HszRu6baQP2/2mCj92auNCubmq1chI3IgAghkKQHWa/bTTQOSynDNmjWqVKnSGVXvv/9+DR8+POW1+Ph4DRo0yPxu/fr1KleunDp37qwnnnhCMTGh3QWGgrYvaBIQyGoCVw2eYm6Dmy9nNl1cpoC5Na53R6roqKS/zf/TFBWl6JTb6Xp7JP0uJlq6umYpta9TJqux8XkRQACBkAVYr4VMl3IgDcgZGpD27dvr5ptvPkm3atWquuyyy1J+9/DDD2vo0KG699571bhxY02dOlUffvihvN+/9dZbIc0MBR0SGwchkKUFGvQdqx0HjqrFBcX14b0Ns7QFHx4BBBDwQ4D1mr0yDcgZGpCePXuqb9++Z9VdsGCBateura5du2rw4MEp+3Xr1k1DhgzRvHnzVLNmzaBnh4IOmowDEIh4gYSERC3fdkAz1+zSn2t2afnWA4qNTzj+AMATZzGSz2SkPBTQnNGI0sJNe821H7fUL6tBt9SOeA8+AAIIIJDRBViv2c8QDchZGhCvCfG23Llzn6bsvda/f3+tWrXqpK9srV69WpUrV1ZaDczZpo2Cti9oEhCIFIG56/fo6e/ma/WOgzoal2D9tru1qqbubc63ziEAAQQQQODcAqzX7CuEBuQMDUi+fPl04MAB84r31avu3bubr1Ylb23btjVnObZs2XLaDJQsWVJ169bVqFGjgp4dCjpoMg5AIGIFenw9V9/P2Xja+y9bOLeK58+Z8hBAs8MZ7mblnfUwL0kqXTCX+t1YU+cVzBWxHrxxBBBAIFIEWK/ZzxQNSCrDdevW6b777tONN96o8uXLa9OmTXr33XflFdrjjz9uLjr3Nu/rVTly5NDs2bNPm4F69eopNjZW3te0zrVt3rxZ3p/U25IlS9SpUyeT6+WwIYBA5hW4+78zNWnZdvMBOzetpLrlC6tJ1aIqlCdH5v3QfDIEEEAgEwjQgNhPIg1IGobe3a6aN2+u6dOna9myZapSpYr5453pmDZt2mlHexekb9u2TStWrDhncu/evdWnT58z7kMDYl/YJCCQ0QVueGuqvK9hnV8yn0Z3b57R3y7vDwEEEEDguAANiH0p0IAEYDhixAhdf/31euedd/TPf/6TMyABmLELAgicW6DlKxPN9R+XVCysb7o0hgsBBBBAIEIEaEDsJ4oGJADD+fPnm7te9evXT88++6y4BiQANHZBIAsKHImN16rtB7VuV9KF5cnP4jB/m2dvJD1/w/v58W/mad+ROLWuXkLD774kC2rxkRFAAIHIFKABsZ83GpAADL///nt16NDBPHTQeyCh14QMGDDgrHfB8l73mpVgNwo6WDH2RyD9BeZv2KMvZq7T7LW7tWLbASUcvzg80Hd2U70yevXWOoHuzn4IIIAAAukswHrNfgJoQFIZetdulChR4iTVw4cPq0mTJlq4cKFWrlxpnnju3QHLu9PV2Z4DMnfuXNWqVSvo2aGggybjAATSVWDb/iNq+vIEHYsP/Ta6L99UU7c3LJ+un4PBEUAAAQQCF2C9FrjV2fakAUkl4939aufOnbriiitUtmxZcxesjz76yJzp8M54PP300yl7d+nSxVwT4j0J3WtQvCehf/DBB3rwwQc1bNiwkGaGgg6JjYMQSDeBaSt36I73ZqSMf2WNkqpRuoDKFc5jbonr3So34fgtdL2/vXvmmn8fP0tSokBO1SxT0Hwtiw0BBBBAIDIEWK/ZzxMNSCrD999/3zQcf//9t3bt2iXveSDe7XC9J5x7F6Gn3uLi4jRw4EDztayNGzeqTJky6ty5s5588klly5YtpJmhoENi4yAE0k1g1MLN6vLpHDP+O3fVV9uLzku398LACCCAAAL+CLBes3emAbE3dJZAQTujJAgBXwS+nrVeT34334z1bZdGalCxiC/jMggCCCCAQPoJsF6zt6cBsTd0lkBBO6MkCAFfBN6bvEr9fl1ixhrT/XJVK5nfl3EZBAEEEEAg/QRYr9nb04DYGzpLoKCdURKEQFAC+4/EasPuw1q6ZZ9mrt6lPYditfPAsaSMKCnau32uohQdnfR38iUb63cd0pqdh8xuM55tpZIFcgU1LjsjgAACCESeAOs1+zmjAbE3dJZAQTujJAiBgAXmrNutu4bP0MFj8QEfc+qOMdFRWtSnrXJljwk5gwMRQAABBCJDgPWa/TzRgNgbOkugoJ1REoRAwAIDRi7RO5NWnbR/jpgyzEjfAAAgAElEQVRoFc+fUzmzRXs3rlJiYqJ5vkeiEpVw/I673u+817yHDN5xaXk90rJqwGOyIwIIIIBA5AqwXrOfOxoQe0NnCRS0M0qCEAhY4Onv5uvLWevN/gM71FKTasVUplDugI9nRwQQQACBrCXAes1+vmlA7A2dJVDQzigJQiBggYc+na2RC7cof85sWtCnbcDHsSMCCCCAQNYUYL1mP+80IPaGzhIoaGeUBCEQsMAd7/2haSt3mrMeU5++IuDj2BEBBBBAIGsKsF6zn3caEHtDZwkUtDNKgrKQgHctxvwNe/X31v0pnzr5ueKpnzB+4ndJuyXfyerVMcu0ftdh1ShVQL92a5aF5PioCCCAAAKhCLBeC0Xt5GNoQOwNnSVQ0M4oCcrkAlv3HdF/Ri7Vsm37tWr7QR2yuINVMtVllYvoy382yuRyfDwEEEAAAVsB1mu2ghINiL2hswQK2hklQZlc4OWRSzVs0kqnn/Lfbc5X11bVnGYShgACCCCQ+QRYr9nPKQ2IvaGzBAraGSVBmVygyyezNWrRFvMpG1Ysouql8qt+xSLH716VqETv/rje7XOPO6T8+/gPp/6+YO7sJiP1V7YyOSEfDwEEEEAgRAHWayHCpTqMBsTe0FkCBe2MkqAIFUhISNTqnQe1ff9R8wlSX8uRfM2G97uePyw013xUKZ5X4/7dIkI/LW8bAQQQQCASBViv2c8aDYi9obMECtoZJUERJPDT3I365s8NOngsTiu2HtD+o3EBv/umVYvp086XBrw/OyKAAAIIIGArwHrNVpBrQOwFHSZQ0A4xiYoIgWNxCarV5zcdiT3+ePEg3/W/WlVTjzbnB3kUuyOAAAIIIBC6AOu10O2Sj+QMiL2hswQK2hklQREisHHPYTV5ebx5t4XyZFe5wnnUsFIRnV8yn7LHRJvfn349R9IVHIXz5NDl5xdXjmxJ+7EhgAACCCDghwDrNXtlGhB7Q2cJFLQzSoIygID3fI7DsfGmgfCu34g+fhGH91eU958oaf6GPeowdLp5t/1uvFh3XlohA7xz3gICCCCAAAJnF2C9Zl8dNCD2hs4SKGhnlAT5LOBdPD5rzS5t239Uew/HauLf2zV3/W7tOHAs4Hcy/B8N1LpGyYD3Z0cEEEAAAQTSQ4D1mr06DYi9obMECtoZJUE+C3w2Y625M5XNNuHxFqpULK9NBMcigAACCCAQdgHWa/bENCD2hs4SKGhnlAT5LND9q7n64a+NJ41atUQ+1S1XSPlzZVdC8vM3EhPNszm8f3o/Jf0tXVa5qK6vXdrnd81wCCCAAAIIBC/Aei14s1OPoAGxN3SWQEE7oyTIZ4G73p+hKct3qGjeHBp8e11VKp73+EMBfX4jDIcAAggggECYBViv2QPTgNgbOkugoJ1REhSiwLb9R3ToaPzpF41HRSk61cXj3sMAvaeGJ11QLt05fIaWbtmvi0oX0C//ahbi6ByGAAIIIIBAxhdgvWY/RzQg9obOEihoZ5QEpSGw59AxzVm3W0djE7Rky37NXrtLCzbs1b4jgT8E8ExDeLfF/fi+hvgjgAACCCCQaQVYr9lPLQ2IvaGzBAraGSVB5xA4Ehuv5oMmaOu+o86dHmhWST2vqeE8l0AEEEAAAQQyigDrNfuZoAGxN3SWQEE7oyToHAKLNu3VNW/8ftIe2aKjdHGZgqpbvpCK58954uF/iYlK8C4UT33R+CkXkie/XiRvdt3WoLwK5smOPwIIIIAAAplWgPWa/dTSgNgbOkugoJ1REnQOgXFLtur+j/40ezzcoopaVS+hmmUK8URxqgYBBBBAAIEABFivBYCUxi40IPaGzhIoaGeUWSroWFyC4hMS9ffW/fpzzS79uWa3lmzZZ67vSL5I3Ltg3Nu8vw4fi9fOg0kPCPzh4caqW75wlvLiwyKAAAIIIGAjwHrNRi/pWBoQe0NnCRS0M8osEzRg5BK9N3mV+ZpUsJt3V6uZPVurWL6cwR7K/ggggAACCGRZAdZr9lNPA2Jv6CyBgnZGmSWCEhMTdcFzo+SdATl1q1I8r8oWznP8oX9J3UnydRwJCUlnQq6qWUp3XVYhS1jxIRFAAAEEEHAlwHrNXpIGxN7QWQIF7YwySwTtOnhM9V4aYz5rrbIF1axaMdWvUFj1KxRRwdxcCJ4lioAPiQACCCDguwDrNXtyGhB7Q2cJFLQzykwRtGH3If29Zb/iEhKP35Xq5O9Zbd57RH1GLDaftdc11dW5WeVM8bn5EAgggAACCGRkAdZr9rNDA2Jv6CyBgnZGGXFB3teopq7YoRXbDmjfkViNWbzVPFk80G1Ix7q6rnbpQHdnPwQQQAABBBAIUYD1WohwqQ6jAbE3dJZAQTujjLigXj8u0Kd/rAvpfeeIidbEJ1qodKHcIR3PQQgggAACCCAQuADrtcCtzrYnDYi9obMECtoZZcQFtXt98klnPGKio3R5tWJqVb2kShfKpShFyfuvtx3/y1xg7v1P9VIFdF7BXBH3mXnDCCCAAAIIRKIA6zX7WaMBsTd0lkBBO6OMuKBL+o3V9v1HlT9XNv36r2YqlCe78ufiQvKIm0jeMAIIIIBAphdgvWY/xTQg9obOEihoZ5QRFZSQkKhqvUaahwneVK+MXr21TkS9f94sAggggAACWUmA9Zr9bNOA2Bs6S6CgnVGmS9DEv7ep5w8LzUXk2aKjFJ3q6eNJbyjKPH8j6aekZ3Ekb1v3HTU/PtCsknpeUyNd3j+DIoAAAggggEDaAqzX0jZKaw8akLSEfHydgvYROwxD/eO/MzV52Xar5BfbX6R/NKpolcHBCCCAAAIIIBA+AdZr9rY0IPaGzhIoaGeU6RLU5tVJWr7tgArnya4GFYuc9OwO7ynk3ub95T3BPPln8/fx18oVya2n2l3ItR/pMnsMigACCCCAQGACrNcCczrXXjQg9obOEihoZ5S+B+08cFT1+441415Zo6Te/UcD398DAyKAAAIIIIBA+AVYr9kb04DYGzpLoKCdUfoatPdwrC7tP1ZHYhPMuJ0uK6++N9T09T0wGAIIIIAAAgj4I8B6zd6ZBsTe0FkCBe2M0tegaSt26I7hM1LGfKNjXV3PU8l9nQMGQwABBBBAwC8B1mv20jQg9obOEihoZ5S+Bv3w1wZ1/2qeGfOJthfo4RZVFJX6Fle+vhsGQwABBBBAAIFwCrBes9elAbE3dJZAQTuj9C3Iu6C88cvjtXnvETPm2B7NVbVEPt/GZyAEEEAAAQQQ8FeA9Zq9Nw2IvaGzBAraGaVvQTNW7dRt7/6RMt6C3ldyFyvf9BkIAQQQQAAB/wVYr9mb04DYGzpLoKCdUfoW9O3sDXr8m6SvX9UoVUC/dmvm29gMhAACCCCAAAL+C7BeszenAbE3dJZAQTuj9C3oo2lr9ML/Fpnxfu7aVBeXKejb2AyEAAIIIIAAAv4LsF6zN6cBsTd0lkBBO6P0LeitCSs06Le/zXiTnmihCkXz+jY2AyGAAAIIIICA/wKs1+zNaUDsDZ0lUNDOKH0LGjhqqd6euNKM92ev1iqWL6dvYzMQAggggAACCPgvwHrN3pwGxN7QWQIF7YzSt6Dnf1qoj6evNeMtfamdcmWP8W1sBkIAAQQQQAAB/wVYr9mb04DYGzpLoKCdUfoStOvgMdV7aYwZKyY6Siv6XcXzP3yRZxAEEEAAAQTST4D1mr09DYi9obMECtoZpS9BXT6ZrVGLtpixCuTKpvm92/oyLoMggAACCCCAQPoJsF6zt6cBsTd0lkBBO6P0JajlKxO1esdBM9ZtDcrpPzfX8mVcBkEAAQQQQACB9BNgvWZvTwNib+gsgYJ2RulLUIO+Y7XjwFFVK5FPo7tfztevfFFnEAQQQAABBNJXgPWavT8NiL2hswQK2hmlL0HVnxulw7Hxuq52aQ3pWNeXMRkEAQQQQAABBNJXgPWavT8NiL2hswQK2hll2IPi4hNUtedIM07HhuU04Ca+fhV2dAZAAAEEEEAgAwiwXrOfBBoQe0NnCRS0M8qwB+09HKvafUabcR5oVkk9r6kR9jEZAAEEEEAAAQTSX4D1mv0c0IDYGzpLoKCdUYY9aNrKHbrjvRlmnMdaV9Njrc8P+5gMgAACCCCAAALpL8B6zX4OaEDSMBw/frxatWpl9lq+fLmqVq2acsThw4fVu3dvff7559q+fbt5rVu3bnrggQdCmhkKOiQ23w9at/OQLh80IWXcXtdUV+dmlX1/HwyIAAIIIIAAAv4LsF6zN6cBOYdhbGysatWqpfXr1+vgwYOnNSDXXHONRo8erUcffVQ1atTQL7/8op9++kkDBw7UE088EfTsUNBBk6XLASPmbVLXL/5KGfunR5qodrlC6fJeGBQBBBBAAAEE/BVgvWbvTQNyDsMBAwbo9ddf1x133GH+Tn0G5Oeff9Z1112nV199Vd27d09Jad++vcaMGaO1a9eqePHiQc0QBR0UV7rt/OXMdXr6+wVm/MG311H7OmXS7b0wMAIIIIAAAgj4K8B6zd6bBuQshuvWrVP16tX15ptvmmaiT58+JzUgd955p77//nvt2rVLuXPnTkmZMGGCrrjiCr333nvq3LlzUDNEQQfFlW47D5+ySn1/WWLGH9P9clUrmT/d3gsDI4AAAggggIC/AqzX7L1pQM5ieOONN2rz5s2aPn26aT5ObUAuuOACFS5cWH/88cdJCYcOHVLevHn14IMPatiwYUHNEAUdFFe67fzamGUaPG65GX/6M1eoVMETDWi6vSkGRgABBBBAAAFfBFiv2TPTgJzB0LuW4/rrr9eMGTPUoEEDc6H5qQ1I/vz5deWVV+q77747LaFIkSJq0qSJRowYcdYZ8pob70/qbcmSJerUqZNmz56tevXq2c8uCWER6PvzYg3/fbXJnt/7ShXIlT0s4xCKAAIIIIAAAhlPgAbEfk5oQE4xPHLkiC666KKUr1F5L5+pAYmJidFtt91m7oB16la6dGlzUfrYsWPPOkPJmWfagQbEvrDDlRCfkKgrX5ukldsPmiFW9r9aMdFR4RqOXAQQQAABBBDIYAI0IPYTQgNyiuELL7ygIUOGaNmyZSpWrJh5lTMg9oWWWRJ6/rBAn81YZz5OnhwxWvxiu8zy0fgcCCCAAAIIIBCAAA1IAEhp7EIDkgrI+0pUpUqVzF2tvGs4kjfvDliDBw/WxIkTzevly5cX14DYF18kJrR5dZKWbztg3nqTqkX1WefLIvFj8J4RQAABBBBAIEQBGpAQ4VIdRgOSCmPu3LmqW7fuOVW9C8wPHDhgbs37ww8/nPUuWO+++27QDySkoO0LOtwJjQaM0+a9R8wwi19sqzw5soV7SPIRQAABBBBAIAMJsF6znwwakFSGe/fulXcb3VO3L7/8Ul999ZWGDh2qsmXL6tprrzUXmHsXqp/tOSBr1qxRiRIlgpohCjoornTZuWbv37T/SJyuuvg8De1UP13eA4MigAACCCCAQPoJsF6zt6cBCcDwTNeAeIe1a9dO48aNU9euXc1F597DCb0noXsPMHz66acDSD55Fwo6aDJfD0hMTFTlZ39VYqJ0S/2yGnRLbV/HZzAEEEAAAQQQSH8B1mv2c0ADEoDh2RoQ75kf3kXrX3zxhbZv364qVaqoW7duJ10/EkB8yi4UdDBa/u978GicLnrhNzPwPY0rqvf1F/n/JhgRAQQQQAABBNJVgPWaPT8NiL2hswQK2hllWIKWbd2vK1+bbLK7XlFV/77ygrCMQygCCCCAAAIIZFwB1mv2c0MDYm/oLIGCdkbpPGjbviNq2H9cSu7TV12oLs2rOB+HQAQQQAABBBDI2AKs1+znhwbE3tBZAgXtjNJ50MgFm/XQZ3NScj/vfKkaV016TgwbAggggAACCGQdAdZr9nNNA2Jv6CyBgnZG6Tzomz/X64lv55tc7+zHg5dXVlQUT0B3Dk0gAggggAACGVyA9Zr9BNGA2Bs6S6CgnVE6D/pg6mr1GbHY5I7s1kzVSxVwPgaBCCCAAAIIIJDxBViv2c8RDYi9obMECtoZpfOgN8cv1yujl5ncKU+2VLkieZyPQSACCCCAAAIIZHwB1mv2c0QDYm/oLIGCdkbpPGjAyCV6Z9Iqk/vXc21UOG8O52MQiAACCCCAAAIZX4D1mv0c0YDYGzpLoKCdUToPeub7Bfpi5jqTu6zvVcqRLdr5GAQigAACCCCAQMYXYL1mP0c0IPaGzhIoaGeUToO8J6BXeuZXk+k1Hl4DwoYAAggggAACWVOA9Zr9vNOA2Bs6S6CgnVE6Ddpz6JjqvDgmJXPNy9c4zScMAQQQQAABBCJHgPWa/VzRgNgbOkugoJ1ROg3asPuQmv5ngsn8R6MKerH9xU7zCUMAAQQQQACByBFgvWY/VzQg9obOEihoZ5ROg5Zu2ad2r08xmc9dW0P3N63kNJ8wBBBAAAEEEIgcAdZr9nNFA2Jv6CyBgnZG6TRo9tpd6jB0uskc2KGWbr2knNN8whBAAAEEEEAgcgRYr9nPFQ2IvaGzBAraGaXToEnLtuvu/840mW/fWU9X1yzlNJ8wBBBAAAEEEIgcAdZr9nNFA2Jv6CyBgnZG6TToh782qPtX80zmR/c1VPPzizvNJwwBBBBAAAEEIkeA9Zr9XNGA2Bs6S6CgnVE6C5q9drc6DJ2WkvfdQ41Vv0JhZ/kEIYAAAggggEBkCbBes58vGhB7Q2cJFLQzSmdBL45YrP9OXZ2SN7NnK5XIn8tZPkEIIIAAAgggEFkCrNfs54sGxN7QWQIF7YzSWVCPr+fq+zkbTd5nnS9Vk6rFnGUThAACCCCAAAKRJ8B6zX7OaEDsDZ0lUNDOKJ0Fdf5olsYu2aYieXNoznNtnOUShAACCCCAAAKRKcB6zX7eaEDsDZ0lUNDOKJ0F3TJsmmat2a3KxfJq/OMtnOUShAACCCCAAAKRKcB6zX7eaEDsDZ0lUNDOKJ0E7T0Uq9ovjjZZdcoV0o+PNHGSSwgCCCCAAAIIRK4A6zX7uaMBsTd0lkBBO6N0EtTvl8V6b0rSBejerXe9W/CyIYAAAggggEDWFmC9Zj//NCD2hs4SKGhnlE6CvIcPeg8h9LbBt9dR+zplnOQSggACCCCAAAKRK8B6zX7uaEDsDZ0lUNDOKJ0E3fDWVM1dv0dlC+fW709d4SSTEAQQQAABBBCIbAHWa/bzRwNib+gsgYJ2Rukk6IpXJmrVjoNqUKGwvn2osZNMQhBAAAEEEEAgsgVYr9nPHw2IvaGzBAraGaWToPovjdHOg8fU6sISev+eS5xkEoIAAggggAACkS3Aes1+/mhA7A2dJVDQziitghITE/Xkt/P1zewNJuemumX06m11rDI5GAEEEEAAAQQyhwDrNft5pAGxN3SWQEE7o7QKWr3joFq+MjEl419XVFWPKy+wyuRgBBBAAAEEEMgcAqzX7OeRBsTe0FkCBe2M0iror3W7dePb00xGruzRmvrUFSqaL6dVJgcjgAACCCCAQOYQYL1mP480IPaGzhIoaGeUVkFTlm/XXe/PNBlDOtbVdbVLW+VxMAIIIIAAAghkHgHWa/ZzSQNib+gsgYJ2RmkVNGrhZnX5dI7J+ODeS9TyghJWeRyMAAIIIIAAAplHgPWa/VzSgNgbOkugoJ1RWgV98+d6PfHtfJPxTZdGuqRiEas8DkYAAQQQQACBzCPAes1+LmlA7A2dJVDQziitgj6culq9Ryw2GSO7NVP1UgWs8jgYAQQQQAABBDKPAOs1+7mkAbE3dJZAQTujtAp6c/xyvTJ6mcmY8mRLlSuSxyqPgxFAAAEEEEAg8wiwXrOfSxoQe0NnCRS0M0qroJdHLtWwSStNxl/PtVHhvDms8jgYAQQQQAABBDKPAOs1+7mkAbE3dJZAQTujtArq9eMCffrHOpOxrO9VypEt2iqPgxFAAAEEEEAg8wiwXrOfSxoQe0NnCRS0M0qroO5fzdUPf200jYfXgLAhgAACCCCAAALJAqzX7GuBBsTe0FkCBe2M0iqo80d/auySrSqSN4fmPNfGKouDEUAAAQQQQCBzCbBes59PGhB7Q2cJFLQzSqug29+drj9W7VL5Ink0+cmWVlkcjAACCCCAAAKZS4D1mv180oDYGzpLoKCdUVoFXTfkdy3YuNfcfte7DS8bAggggAACCCCQLMB6zb4WaEDsDZ0lUNDOKEMO+nLmOj39/QJzfMOKRfR1l0YhZ3EgAggggAACCGQ+AdZr9nNKA2Jv6CyBgnZGGVLQsbgE1ez9m47GJZjj2110nobdVT+kLA5CAAEEEEAAgcwpwHrNfl5pQOwNnSVQ0M4oQwradfCY6r00JuXYnx5potrlCoWUxUEIIIAAAgggkDkFWK/ZzysNiL2hswQK2hllSEHrdh7S5YMmmGMfv/J8PXpFtZByOAgBBBBAAAEEMq8A6zX7uaUBsTd0lkBBO6MMKWjxpn26+o0p5tje19XQPU0qhZTDQQgggAACCCCQeQVYr9nPLQ2IvaGzBAraGWVIQbPW7NItw6abYwfdXEu3NCgXUg4HIYAAAggggEDmFWC9Zj+3NCD2hs4SKGhnlCEFTVi6Tfd+OMscO6xTPbW7uFRIORyEAAIIIIAAAplXgPWa/dzSgNgbOkugoJ1RBh3kXYDe5OXxOhwbb4795P6GalateNA5HIAAAggggAACmVuA9Zr9/NKA2Bs6S6CgnVEGHTR04kr9Z9TSlON+7tpUF5cpGHQOByCAAAIIIIBA5hZgvWY/vzQg9obOEihoZ5RBB/X+3yJ9OG2NOe6GOqX12m11FBUVFXQOByCAAAIIIIBA5hZgvWY/vzQg9obOEihoZ5RBBz3+zTx9O3uDoqOklf2vpvkIWpADEEAAAQQQyBoCrNfs5znTNSCxsbE6cOCAChcubK/jcwIF7TN4quG6fDJboxZtUf5c2bSgd9v0eyOMjAACCCCAAAIZWoD1mv30RGwD8u2332r69On6v//7vxSFl19+WS+88ILi4uJ01VVX6ZtvvlHu3LntlXxKoKB9gj7DMHe9P0NTlu9Q6YK5NO2ZVun3RhgZAQQQQAABBDK0AOs1++mJ2AakcePGqlatmj766COjMHfuXNWvX1/e7y+44AJ9+OGH6t27t3r16mWv5FMCBe0T9BmGueGtqZq7fo+qlcinMT2ap98bYWQEEEAAAQQQyNACrNfspydiG5BixYqpZ8+e6t69u1F49tln9e6772rjxo3KmTOnHnroIU2ePFmLFi2yV/IpgYL2CfqUYcYs3qoHPv7T/LZu+UL64eEm6fNGGBUBBBBAAAEEMrwA6zX7KYrYBiRXrlwaOnSo7r33XqPQsGFDVa1aVZ9//rn59/vvv6/HHntM+/fvt1fyKYGC9gk61TDxCYmq8+Jo7T8SZ37b8oLi+uDehv6/EUZEAAEEEEAAgYgQYL1mP00R24BUqlRJd9xxh/r166ft27erVKlSeuutt/Tggw8alddee00vvviidu/eba/kUwIF7RN0qmH2Ho5V7T6jU37DAwj9nwNGRAABBBBAIJIEWK/Zz1bENiCPPPKIPvvsMz333HMaP368xo0bp1WrVql06dJGpXPnzpo1a5bmzZtnr+RTAgXtE3SqYTbtOazGL483v3msdTU91vp8/98EIyKAAAIIIIBAxAiwXrOfqohtQHbu3KkOHTqY6zy8az68Mx5dunQxIocPHzaNyD333GN+HykbBe3/TC3ful9tXptsBu51TXV1blbZ/zfBiAgggAACCCAQMQKs1+ynKmIbkOSPvm/fPnOr3ezZs6doeA3IsmXLVK5cORUpUiRgpSVLlqhPnz6aPXu2Nm/erOjoaFWpUsVcZ+I1Nzly5EjJio+P16BBgzR8+HCtX7/ejOWddXniiScUExMT8Jipd6SgQ2KzOmjOut266e1pJmPATTXVsWF5qzwORgABBBBAAIHMLcB6zX5+I74BsSc4kTB69Gi98soruvTSS1W2bFl5TcbUqVP1xRdf6Prrr9ePP/6YsvPDDz+cchG8d+tfbz/v1r/e771rUULZKOhQ1OyOmbJ8u+56f6YJGdKxrq6rnfQVPjYEEEAAAQQQQOBMAqzX7OsiYhuQiRMnyiuAHj16pCh8+umn5rkf3oXnd911l9544w1zFsN2e/TRR01TsXTpUvOMkQULFqh27drq2rWrBg8enBLfrVs3DRkyxFx3UrNmzaCHpaCDJrM64NCxON0ybLoWbdpncj645xK1vLCEVSYHI4AAAggggEDmFmC9Zj+/EduAXHHFFSpatKh52rm3rVixQhdddJHKly+vihUrmgvTvebAax5sN++siPfVqj/++MOcHfGeP9K/f39z0bt3N67kbfXq1apcubJ5vW/fvkEPS0EHTWZ1wMfT1+j5n048J+a7hxqpfoXAv7JnNTgHI4AAAggggEBECrBes5+2iG1ASpYsqccff9w0Bt7mLfgHDhxorscoWLCg7rzzTnnXdHhFEux26NAheX8OHjyomTNnyrvjVrZs2bRy5UpzvUnbtm3NWY4tW7acFu29r7p162rUqFHBDmveq/c0d+8alHr16gV9PAcEJ9D358Ua/vtqc9DFZQrox4ebKFuM/Rmz4N4FeyOAAAIIIIBAJAmwXrOfrYhtQLwHEQ4bNszc6crbmjVrpuLFi+v77783/37vvfdMc7Jnz56glXr37m0uRk/eLrnkEvOU9Tp16phfeV+v8i5I9xqFUzevcYiNjTVf0zrX5l3k7v1JvXkNU6dOnWhAgp6x0A545vv5+mLmenPwyv5XKyY6KrQgjkIAAQQQQACBLCNAA2I/1RHbgHh3nXrggQf0/PPPm6edFytWzJwB8a7D8DbvWgzvq1DeXbKC3byvVnl/vFv9el/lWrhwoXngYYsWLUyUd2cs70zHtGlJd09KvXkXpG/bts18Jexc26lNTup9OQMS7IyFtv+jn8/Rz/M3K3f2GC15qa5j+gUAACAASURBVF1oIRyFAAIIIIAAAllKgAbEfrojtgHxzhR4zcGbb76pX375RR999JG59a53DYa3edd+eK8vXrzYWsl7lshTTz1lvnZVvXp1zoBYi2aMgHs/mKkJf29X8fw5Natn64zxpngXCCCAAAIIIJChBWhA7KcnYhuQtWvXqk2bNilnGryzHS+99JIRiYuLM8/laN++vfmalu22detWnXfeeXr22WfNmRCuAbEVzRjH3zJsmmat2a1KxfJqwuNJZ7fYEEAAAQQQQACBcwnQgNjXR8Q2IMmNhneGo1ChQubuV8mb97WrCRMmmFvlenfEst3WrFlj7nb10EMP6e233zaNyIABA856F6zkRiXYcSnoYMVC3//vLfvV9vWkJ6DXLFNQI7o2DT2MIxFAAAEEEEAgywiwXrOf6ohuQOw//skJ3rUbJUqc/hyIJ5980jz13HvQ4N13322+iuXd6epszwGZO3euatWqFfTbo6CDJgvpgMTERDX9zwRt3HPYHN+oclF98c/LQsriIAQQQAABBBDIWgKs1+znO1M0IAcOHDB3u0pISDhNJPWZkbS4brzxRnPhuXexufcVLi/zt99+07hx49S0aVNzVsW7Ha+3denSRe+8847uvfdeNWnSxDwJ/YMPPtCDDz4Y8te+KOi0ZsjN694DCGs8/1tK2MAOtXTrJeXchJOCAAIIIIAAAplagPWa/fRGdAPinZH4z3/+Yy4+P9sWHx8fsNJXX31lznLMnz9f27dvV86cOXXhhRfqtttuM2c7vH8nb951Jt5dt4YPH66NGzeqTJky6ty5s7yzJclNSsADH9+Rgg5WLLT9t+47okv7jzMH39ekkp6/rkZoQRyFAAIIIIAAAllOgPWa/ZRHbAPy6aef6h//+Ie8J6K3atXK3HK3e/fupknwzkSULVvW3AnL+8pUpGwUtD8ztXzrfrV5Len6j17XVFfnZkl3TmNDAAEEEEAAAQTSEmC9lpZQ2q9HbAPiPfCvQIECmjhxovnalPcQwrFjx5qGxPu3d43GM888Yy4cj5SNgvZnpmav3a0OQ5Oe4fKfDjV12yUnbmDgzztgFAQQQAABBBCIVAHWa/YzF7ENSJ48efTyyy/rX//6l3bv3q2iRYtq5MiR5ha53ubdkvfLL7/UokWL7JV8SqCg/YH+YOpq9RmR9HyYt++sp6trlvJnYEZBAAEEEEAAgYgXYL1mP4UR24AULlzY3ArXuxg8NjZWuXLlMl+98r6W5W3vvvuuHnvsMR06dMheyacECjr80JOWbdfd/52ZMtCn91+qptWKhX9gRkAAAQQQQACBTCHAes1+GiO2AWnQoIG5W9Urr7xiFLwnlJ9//vn66aef5N1m9ZprrtHSpUvNszoiZaOgwz9TA0ct1dsTV6YM9Gev1iqW78TNBcL/DhgBAQQQQAABBCJZgPWa/exFbAPiPezPuxDde0hgdHS0Xn/9dfXo0UOVKyddULx69Wo9//zzeuGFF+yVfEqgoMMP/dyPC/XJH2vNQOP/3VyVi+cL/6CMgAACCCCAAAKZRoD1mv1URmwDsn//fm3atElVqlRJue2td1tc77qPmJgY3XDDDXr66afNz5GyUdDhn6nHvvxLP87dpBwx0VrW76rwD8gICCCAAAIIIJCpBFiv2U9nxDYg9h894yVQ0OGfk/s/nKVxS7epaN4cmv1cm/APyAgIIIAAAgggkKkEWK/ZT2dENSAff/xx0J84+aL0oA9MhwMo6PCj3zpsumau2aWKRfNo4hMtwz8gIyCAAAIIIIBAphJgvWY/nRHVgHjXekRFRZlP7V1ontbm7RvMk9DTygv36xR0eIXf/321Xvo56fa7NcsU1IiuTcM7IOkIIIAAAgggkOkEWK/ZT2nENSC5c+dW+/btdeuttyp//vxpCnhPSY+UjYIO30x5DeuFz43S0bgEM0jr6iU1/O4G4RuQZAQQQAABBBDIlAKs1+ynNaIakK+++kre17BGjx6tnDlzmkbknnvuUevWrVPOjNiTpF8CBR0++0PH4lTj+d9SBhjZrZmqlyoQvgFJRgABBBBAAIFMKcB6zX5aI6oBSf64W7du1SeffGKakYULF6p06dLq1KmTeQhhjRo17FXSKYGCDh/8tn1H1LD/ODNA1yuq6t9XXhC+wUhGAAEEEEAAgUwrwHrNfmojsgFJ/bG9Ivjoo4/0xRdfaOfOnapXr5682/G2bBl5FxhT0PYFfbaEldsPqNX/TTIvP3PVhXqweZXwDUYyAggggAACCGRaAdZr9lMb8Q1IMsHmzZvNGZDx48ebhw96DyGMtI2CDt+Mzd+wR9e/OdUM0PeGi9XpsgrhG4xkBBBAAAEEEMi0AqzX7Kc2ohuQuLg4jRgxwpwBGTVqlNG4+uqr1atXL3MmJNI2Cjp8MzZtxQ7dMXyGGWDw7XXUvk6Z8A1GMgIIIIAAAghkWgHWa/ZTG5ENyMyZM831H95Tz3ft2qUGDRro7rvvVseOHVWkSBF7lXRKoKDDA7/jwFE16Ds2JXz4PxqodY2S4RmMVAQQQAABBBDI1AKs1+ynN6IakAEDBpjGY9myZSkXnnuNx4UXXmgvkQESKOjwTMJ7k1ep369LUsJ/eqSJapcrFJ7BSEUAAQQQQACBTC3Aes1+eiOqAfEeROg9B+T666+X93wP799pbffdd19au2SY1yno8EzFgJFL9M6kVSb8lvplNfDmWpnits3h0SIVAQQQQAABBM4lwHrNvj4irgFJ/sjeU87Teho6T0K3L5DMkNDzhwX6bMY681GW9b1KObKl3bhmhs/NZ0AAAQQQQAAB9wI0IPamEdWATJqUdBvVYLbmzZsHs3u67ktBh4e/25d/6ae5m0zj4TUgbAgggAACCCCAQKgCrNdClTtxXEQ1IPYfN2MnUNDhmZ/OH83S2CXbVDRvDs1+rk14BiEVAQQQQAABBLKEAOs1+2mmAbE3dJZAQTujPCnotnema8bqXSpfJI8mPxl5D6gMjwqpCCCAAAIIIBCKAOu1UNROPoYGxN7QWQIF7YwyJeiX+Zv1yOdzzL9rlCqgX7s1cz8IiQgggAACCCCQZQRYr9lPNQ2IvaGzBAraGWVKUJ0XR2vPoVjz76ZVi+nTzpe6H4REBBBAAAEEEMgyAqzX7KeaBsTe0FkCBe2M0gQdjYvXBb1GpYR+1vlSNalazO0gpCGAAAIIIIBAlhJgvWY/3TQg9obOEihoZ5QmaNfBY6r30hjz80MtquipdpnjgZVulUhDAAEEEEAAgWAEWK8Fo3XmfWlA7A2dJVDQzihN0Lqdh3T5oAnm5yfaXqBHWlZ1OwBpCCCAAAIIIJDlBFiv2U85DYi9obMECtoZpQlavGmfrn5jivm593U1dE+TSm4HIA0BBBBAAAEEspwA6zX7KacBsTd0lkBBO6M0QTNX79Kt70w3P79yS23dXL+s2wFIQwABBBBAAIEsJ8B6zX7KaUDsDZ0lUNDOKHXwaJyuf/N3rdx+0IQO61RP7S4u5W4AkhBAAAEEEEAgSwqwXrOfdhoQe0NnCRS0M0p9+sda9fpxYUrgFw9cpkZVirobgCQEEEAAAQQQyJICrNfsp50GxN7QWQIF7YxSg35bqrcmrDSBdcoV0rddGilbTLS7AUhCAAEEEEAAgSwpwHrNftppQOwNnSVQ0M4o1ft/i/ThtDUmcMmL7ZQ7R4y7cJIQQAABBBBAIMsKsF6zn3oaEHtDZwkUtDNK/fvrefpuzgZFR0kr+1+tqKgod+EkIYAAAggggECWFWC9Zj/1NCD2hs4SKGhnlOryyWyNWrRFBXJl0/zebd0Fk4QAAggggAACWVqA9Zr99NOA2Bs6S6CgnVHqrvdnaMryHSpdMJemPdPKXTBJCCCAAAIIIJClBViv2U8/DYi9obMECtoN5aJNe3XNG7+bsPNL5tPo7s3dBJOCAAIIIIAAAllegPWafQnQgNgbOkugoN1QNh80QWt3HjJh9SsU1ncPNXYTTAoCCCCAAAIIZHkB1mv2JUADYm/oLIGCtqeMT0hUlWd/TQka2KGWbr2knH0wCQgggAACCCCAgCTWa/ZlQANib+gsgYK2p9x3JFa1eo82QXc3qqA+7S+2DyUBAQQQQAABBBA4LsB6zb4UaEDsDZ0lUND2lJv3HlajAeNNULdW1dS9zfn2oSQggAACCCCAAAI0IM5qgAbEGaV9EA2IveHyrfvV5rXJJqjXNdXVuVll+1ASEEAAAQQQQAABGhBnNUAD4ozSPogGxN7wr3W7dePb00zQgJtqqmPD8vahJCCAAAIIIIAAAjQgzmqABsQZpX0QDYidoXcB+lPfzde3szeYoCEd6+q62qXtQjkaAQQQQAABBBBIJcB6zb4caEDsDZ0lUNB2lL/M36xHPp+TEvLBPZeo5YUl7EI5GgEEEEAAAQQQoAFxWgM0IE457cJoQOz8Xh+7TK+PXZ4SsrBPW+XLmc0ulKMRQAABBBBAAAEaEKc1QAPilNMujAbEzq/fL4v13pTVJmTOc21UJG8Ou0CORgABBBBAAAEEThFgvWZfEjQg9obOEihoO8pnvp+vL2auNyEr+1+tmOgou0CORgABBBBAAAEEaECc1wANiHPS0ANpQEK3847s+sVfGjFvk3Jnj9GSl9rZhXE0AggggAACCCBwBgHWa/ZlQQNib+gsgYK2o7z3g5ma8Pd2Fc+fU7N6trYL42gEEEAAAQQQQIAGJCw1QAMSFtbQQmlAQnPzjtq674gu7T/OBFQqllcTHm8RehhHIoAAAggggAACZxFgvWZfGjQg9obOEijo0CmvHTJFCzfuMwEXlymgn7s2Cz2MIxFAAAEEEEAAARqQsNUADUjYaIMPpgEJ3sw7wnsAYdWevyoxMen4bq2qqXub80ML4ygEEEAAAQQQQOAcAqzX7MuDBsTe0FkCBR0a5e6Dx1T3pTHm4HYXnadhd9UPLYijEEAAAQQQQACBNARYr9mXCA2IvaGzBAo6NMpV2w/oiv+bZA5+qt2FeqhFldCCOAoBBBBAAAEEEKABCXsN0ICEnTjwAWhAArdKvedXs9bpqe8WmF8NuKmmOjYsH1oQRyGAAAIIIIAAAjQgYa8BGpCwEwc+AA1I4FbJe/6xaqduf/ePlAOHdaqndheXCj6IIxBAAAEEEEAAgQAEWK8FgJTGLjQg9obOEijo4CnfGLdcr45ZlnLgH8+00nkFcwUfxBEIIIAAAggggEAAAqzXAkCiAQkc6c8//9Snn36q8ePHa/Xq1cqbN68uuugiPfPMM2rd+uQH28XHx2vQoEEaPny41q9fr3Llyqlz58564oknFBMTE/igqfakoINn6zNikT6YusYcOLZHc1UtkS/4EI5AAAEEEEAAAQQCFGC9FiDUOXbjDEgqnJtvvlmTJk1Shw4dVK9ePR04cEAffPCBFi5cqLffflsPPfRQyt4PP/ywhg4dqnvvvVeNGzfW1KlT9eGHH8r7/VtvvRXSzFDQwbN1/2qufvhro7LHRGlZ36sUFRUVfAhHIIAAAggggAACAQqwXgsQigYkMCiviWjQoIFy5syZcsDhw4dVp04dbd++Xdu2bVO2bNm0YMEC1a5dW127dtXgwYNT9u3WrZuGDBmiefPmqWbNmoENmmovCjo4sj2HjqnOi0m33y2RP6dm9jz5LFVwaeyNAAIIIIAAAgikLcB6LW2jtPbgDEhaQpL+/e9/69VXXzVftSpbtqx69uyp/v37a9WqVapUqVJKgve1rcqVK5vX+/btG0DyybtQ0MGR3Tx0mv5cu9scdEHJ/Pqt++XBBbA3AggggAACCCAQpADrtSDBzrA7DUgAhh07dtS3336r3bt3K1++fGrbtq05y7Fly5bTji5ZsqTq1q2rUaNGBZBMAxI0UqoDzu81UsfiEsxv/nl5ZT17dXWbOI5FAAEEEEAAAQTSFKABSZMozR1oQNIgWrJkifkK1rXXXqvvvvvO7O19vSpHjhyaPXv2aUd7147Exsaar2mda9u8ebO8P6k3b6xOnTqZXC+H7ewCR+PidUGvpCbvkoqF9fWDjbj+g4JBAAEEEEAAgbAL0IDYE9OAnMNw7969atSokTZt2mTOeFSoUMHsXaVKFXlnOqZNm3ba0d4F6d61IitWrDjn7PTu3Vt9+vQ54z40IGkX9q6Dx1TvpaTrP7o0r6Knr7ow7YPYAwEEEEAAAQQQsBSgAbEElEQDchZD7+Jz76tWM2fO1MiRI9WyZcuUPTkDYl94tgnrdx1Ss4ETTMzjV56vR6+oZhvJ8QgggAACCCCAQJoCNCBpEqW5Aw3IGYiOHTum66+/XuPGjTPXfrRv3/6kvbgGJM26CvsOizft09VvTDHj9L6uhu5pcuJmAGEfnAEQQAABBBBAIMsK0IDYTz0NyCmGcXFx8p4HMmLECH388ce68847T1N+9tlnNWDAgLPeBct7vV+/fkHPDgUdONk3f67XE9/ONwcMurmWbmlQLvCD2RMBBBBAAAEEEAhRgPVaiHCpDqMBSYWRkJCgO+64Q1999ZWGDRumBx988IzC3vUg3p2uzvYckLlz56pWrVpBzw4FHRjZ7LW71GHo9JSdh95ZT1fVLBXYweyFAAIIIIAAAghYCLBes8A7figNSCrDHj166LXXXtPll1+uBx544DTdNm3amIvPva1Lly565513zJPQmzRpYp6E7j013WtavOYllI2CDkxt2KSVennk0pSdf3+qpcoWzhPYweyFAAIIIIAAAghYCLBes8CjATkdr0WLFpo0adJZVSdMmCBvH2/zvqo1cOBADR8+XBs3blSZMmXUuXNnPfnkk+Zp6aFsFHRgaq/89rfenJB0l7HvHmqk+hWKBHYgeyGAAAIIIIAAApYCrNcsAbkLlj2gywQKOjDN3v9bpA+nrTE7L36xrfLkCK3hC2w09kIAAQQQQAABBE4IsF6zrwa+gmVv6CyBgg6M8vFv5unb2RsUHSWt7H81DyAMjI29EEAAAQQQQMCBAOs1e0QaEHtDZwkUdGCUD306WyMXblH+XNm0oHfbwA5iLwQQQAABBBBAwIEA6zV7RBoQe0NnCRR02pTxCYlq0HeMdh+KVemCuTTtmVZpH8QeCCCAAAIIIICAIwHWa/aQNCD2hs4SKOi0Kbt8MlujFm0xO1YrkU9jejRP+yD2QAABBBBAAAEEHAmwXrOHpAGxN3SWQEGnTXnR86N08Fi82fGmumX06m110j6IPRBAAAEEEEAAAUcCrNfsIWlA7A2dJVDQ56b0vn5V5dlfzU75c2bTH8+2Ut6c3AHLWQEShAACCCCAAAJpCrBeS5MozR1oQNIk8m8HCvrc1vuOxKpW79Fmp85NK6nXtTX8mxxGQgABBBBAAAEEJLFesy8DGhB7Q2cJFPS5KTfvPaxGA8abnbq1qqbubc53Zk8QAggggAACCCAQiADrtUCUzr0PDYi9obMECvrclMu37leb1yabnXpdU12dm1V2Zk8QAggggAACCCAQiADrtUCUaEDslXxKoKDPDT1l+Xbd9f5Ms9OAm2qqY8PyPs0MwyCAAAIIIIAAAkkCrNfsK4EzIPaGzhIo6LNTrti2X61fTTr74W1DOtbVdbVLO7MnCAEEEEAAAQQQCESA9VogSufehwbE3tBZAgV9dsrhU1ap7y9LUnYY3f1ynV8yvzN7ghBAAAEEEEAAgUAEWK8FokQDYq/kUwIFfXbo18cu0+tjl5sdXruttm6sW9anWWEYBBBAAAEEEEDghADrNftq4AyIvaGzBAr67JT9flms96asNjvMe+FKFcyd3Zk7QQgggAACCCCAQKACrNcClTr7fjQg9obOEijos1M+8/18fTFzvdlhZf+rFRMd5cydIAQQQAABBBBAIFAB1muBStGA2Ev5kEBBnx256xd/acS8TcqdPUZLXmrnw2wwBAIIIIAAAgggcLoA6zX7quAMiL2hswQK+syUCQmJunHoNM1bv0fF8+fUrJ6tnZkThAACCCCAAAIIBCPAei0YrTPvSwNib+gsgYI+nfLA0ThdN+R3rd5x0LxYuVhejX+8hTNzghBAAAEEEEAAgWAEWK8Fo0UDYq8V5gQK+nTgsYu3qvPHf6a80Lp6CQ2/+5IwzwTxCCCAAAIIIIDAmQVYr9lXBmdA7A2dJVDQp1P+8NcGdf9qnnmhWbViGnRzbZ1XMJczc4IQQAABBBBAAIFgBFivBaPFGRB7rTAnUNCnA3/yx1o99+NC88KIR5uqZtmCYZ4F4hFAAAEEEEAAgbMLsF6zrw7OgNgbOkugoE+nHDpxpf4zaql5YcLjLVSpWF5n3gQhgAACCCCAAALBCrBeC1bs9P1pQOwNnSVQ0KdTvvLb33pzwgrzgnf3K+8uWGwIIIAAAggggEB6CbBes5enAbE3dJZAQZ9O6X39yvsalrctebGdcueIceZNEAIIIIAAAgggEKwA67VgxU7fnwbE3tBZAgV9gjIxMVGdP/pT45ZuM7/0nny+ot9VioriCejOCo4gBBBAAAEEEAhagPVa0GSnHUADYm/oLIGCPkG5Ytt+tX51csoveP6HszIjCAEEEEAAAQQsBFivWeAdP5QGxN7QWQIFfYJy9trd6jB0mvlFvpzZ9OU/L9PFZbgDlrNiIwgBBBBAAAEEQhJgvRYS20kH0YDYGzpLoKBPUE5atl13/3em+cWbd9TVtbVKO3MmCAEEEEAAAQQQCFWA9VqocieOowGxN3SWQEGfoPx1wWY9/Nkc84sP771ELS4o4cyZIAQQQAABBBBAIFQB1muhytGA2MuFIYGCPoH69az1evK7+eYX33ZppAYVi4RBnEgEEEAAAQQQQCA4AdZrwXmdaW/OgNgbOkugoE9Q/vf31Xrx58XmF6Mea6YLzyvgzJkgBBBAAAEEEEAgVAHWa6HKnTiOBsTe0FkCBZ1E+ca45Xp1zLIU19+faqmyhfM4cyYIAQQQQAABBBAIVYD1WqhyNCD2cmFIoKClw8fiddELo5SQmAScM1u05jzXRnlzZguDOJEIIIAAAggggEBwAqzXgvM6096cAbE3dJZAQUvrdh7S5YMmpJi+fWc9XV2zlDNjghBAAAEEEEAAARsB1ms2eknH0oDYGzpLoKCl2Wt3qcPQ6cb05Ztq6vaG5Z35EoQAAggggAACCNgKsF6zFaQBsRd0mJDVC3r9rkNqNvDE2Y/3726gVtVLOhQmCgEEEEAAAQQQsBPI6us1Oz3OgLjwc5qR1Qv6ld/+1psTVqSYju1xuaqWyO/UmDAEEEAAAQQQQMBGIKuv12zsko/lK1guFB1lZPWC/vfX8/TdnA1Gs9c11dW5WWVHssQggAACCCCAAAJuBLL6es2FIg2IC0VHGVm9oP/x35mavGy7SuTPqZk9WztSJQYBBBBAAAEEEHAnkNXXay4kaUBcKDrKyKoFfeBonLp/NVdjFm81kheVLqBf/tXMkSoxCCCAAAIIIICAO4Gsul5zJ8hF6C4trbOyakF/P2eDenw9L8Xv6prn6e0761t7EoAAAggggAACCLgWyKrrNZeOnAFxqWmZlVULeviUVer7yxKjV6NUAb3RsQ4Xn1vWEocjgAACCCCAQHgEsup6zaUmDYhLTcusrFrQr49dptfHLjd6056+QqUL5baU5HAEEEAAAQQQQCA8All1veZSkwbEpaZlVlYt6L4/L9bw31cbvfm9r1SBXNktJTkcAQQQQAABBBAIj0BWXa+51KQBcalpmZVVC/qZ7+fri5nrjd7K/lcrJjrKUpLDEUAAAQQQQACB8Ahk1fWaS00aEJealllZtaAf/XyOfp6/WXlyxGjxi+0sFTkcAQQQQAABBBAIn0BWXa+5FKUBcalpmZXVCjoxMVH9flmS8vUrnv9hWUAcjgACCPx/e3cD51OZNnD8GgZZpO3FZIfCVJLQjIgZhdBU29amty0KmYy0Vm9mQz3MZtiwtRJmarYZm1a2Um09eXlStA0lRJMGTQwa5LUXLzGG53OfdmYZM+N/zrnP+Z9z/r/z+fiUmftc576/9+U/9zXnDQEEEHBcINLWa06AUoA4oWoxZqQl9Geb98rN05aUa6knYL07jPd/WEwfdkMAAQQQQAABFwQibb3mBCkFiBOqFmNGWkIvWLNdBr20olzr7/d2lKsuOseiHrshgAACCCCAAALOC0Taes0JUQoQJ1Qtxoy0hH7js2/kodk/v4Bw1n2dpHPcWRbl2A0BBBBAAAEEEHBHINLWa06oUoA4oWoxZqQl9EtLi+SJt9YYWu8M7SKXxja0KMduCCCAAAIIIICAOwKRtl5zQpUCxAlVizEjLaGnL/panpq31tBa9Gg3aXZ2PYty7IYAAggggAACCLgjEGnrNSdUKUCcULUYM9ISeuL8tTL1g68NrU9H9ZRzGtSxKMduCCCAAAIIIICAOwKRtl5zQpUCxAlVizEjKaFfWbZZHpuTXy619slr5bRaNS3KsRsCCCCAAAIIIOCOQCSt15wSpQBxStZC3EhJ6N37DkmHjPfk6LGfkerWUi8gTJaoKN6AbiFt2AUBBBBAAAEEXBSIlPWak6QUIE7qmowdKQm9bvuPkvzXD8t10m9sLf0Sm5nUojkCCCCAAAIIIOC+QKSs15yUpQBxUtdk7EhJ6BWb9sot039+AeFTt7SROzqcZ1KK5ggggAACCCCAQHgEImW95qQuBUgF3X379smkSZNkxYoVsnz5ctm+fbv069dPcnNzT5qH0tJSmThxomRnZ8uWLVukadOmkpKSIsOHD5eaNc3fzxApCb14/U7p9+Iyw3PqXQny67aNncxxYiOAAAIIIIAAAtoEImW9pg2skkAUIBVQioqKpHnz5tK4cWNp3769vPPOO1UWIEOGDJHp06fLgAEDJDExUfLy8oxCRX196tSppuctUhL63fxtMuTllYZP7oAO0q1lI9NW7IAAAggggAACCIRDIFLWa07aUoBU0D106JDs2rVLYmNjT6haeAAAIABJREFU5ciRI1KrVq1KC5D8/Hxp166dDB06VCZPnlweZdiwYTJlyhRZvXq1tGnTxtTcRUpC//PTLZL2+ueGzev3d5b2559pyonGCCCAAAIIIIBAuAQiZb3mpC8FSDW61RUgo0aNknHjxsmGDRuMMyZl28aNG6VFixaivj927FhTcxcJCb28aI/cmrm03GX+g1dJy3MbmHKiMQIIIIAAAgggEC6BSFivOW1LAWKxAElOTjbOcqh7RCpuMTExEh8fL/PmzTM1f0FP6GPHjknHcQtl54+Hyl2WjrhaGjesa8qJxggggAACCCCAQLgEgr5ec8OVAsRiAaIur6pdu7Zxs3rFLSEhQUpKSkRdplXVtm3bNlF/jt8KCgqkb9++RkwVI2jbTyWlcvET/y3Kbm3fRCbd1i5ow2Q8CCCAAAIIIBBgAQoQ+5NLAWKxAImLixN1pmPJkp8fJ3v8pm5I37FjhxQWFlYZfcyYMZKenl7p94NagOzad0guH/ueMeYh3eIk7dqL7WcwERBAAAEEEEAAARcFKEDsY1OAWCxAOANiPvmKdu2XbpMWGTumXdtShnS7wHwQ9kAAAQQQQAABBMIoQAFiH58CxGIBwj0g5pPvi+Lv5YYpHxk7/umm1nJPZ95+bl6RPRBAAAEEEEAgnAIUIPb1KUAsFiAjR46U8ePHV/kULPX9jIwMUzMU9IT+eMNu+d3zHxsmT9/eTnonNDHlQ2MEEEAAAQQQQCDcAkFfr7nhSwFisQBRT8BST7qq6j0gq1atkrZt25qawyAn9ObdB2TwzBXy5bYfDJOsu9tLcutzTfnQGAEEEEAAAQQQCLdAkNdrbtlSgFQi/dxzz8l3330nR48eldGjRxuFRu/evY2WN954Y3lhMXjwYMnKyjLehJ6UlGS8CT0nJ0dSU1MlMzPT9BwGOaHvyFoqn2zcU27yj/uukMS4s00bsQMCCCCAAAIIIBBOgSCv19xypQCpRLpZs2ayadOmSudAFRj9+/c3vqdeVDhhwgTJzs6W4uJi4+3pKSkpkpaWJtHR0abnMMgJ3THjPdnxn/d/XBp7urx+f6LUia5p2ogdEEAAAQQQQACBcAoEeb3mlisFiFvSIRwnyAnd+n/myf7DpdKt5TmS07+DREVFhSBCEwQQQAABBBBAwFsCQV6vuSVNAeKWdAjHCWpClx49JnEj3zUEftehqfz5FnP3xoRARxMEEEAAAQQQQMAVgaCu11zB+89BKEDc1D7FsYKa0D/8VCJtxywwRj+wS3N54oZLPKROVxBAAAEEEEAAgdAFgrpeC13AfksKEPuG2iIENaG3fX9QOo9/33Aa1uNCeajXRdrMCIQAAggggAACCLgpENT1mpuGFCBuakfgGRB19uPpBesld0mRMfrHf91KUq5s4SF1uoIAAggggAACCIQuQAESulVVLSlA7BtqixDEhB71Rr68/MnmcqM/924jv+t4njYzAiGAAAIIIIAAAm4KBHG95qafOhYFiNvi1RwviAl9W+YS+bRorzHqM+vVlneGdpFfnVHXQ+p0BQEEEEAAAQQQCF0giOu10EevpyUFiB5HLVGCmNDXTf63FGz7Qc6qV1uWjughtaNraLEiCAIIIIAAAgggEA6BIK7X3HakAHFbPMLOgFw54X3ZsuegdG5xlswa1MlD2nQFAQQQQAABBBAwL0ABYt6s4h4UIPYNtUUIYkLH/2mB7D1QIj1bxUh2v8u1WREIAQQQQAABBBAIh0AQ12tuO1KAuC0eYWdALho1Vw6XHpWb42PlmTsu85A2XUEAAQQQQAABBMwLUICYN+MMiH0zxyIELaG3f/+TdBq/0PC6u9P58uRvL3XMjsAIIIAAAggggIAbAkFbr7lhRgESDuUQjxm0hL49c6ksK9pjjH5w1zh57LqLQ5SgGQIIIIAAAggg4E2BoK3XwqHMJVjhUK/imEFL6I4Z78mOHw8Zo/1namfp2PxMD2nTFQQQQAABBBBAwLxA0NZr5gXs70EBYt9QW4QgJfSxY8ek5ePzjPs/urU8R3IHdNTmRCAEEEAAAQQQQCBcAkFar4XLkAIkXPKVHDdICb3v0BG5dPR8Y5T9E5vJmBtbe0iariCAAAIIIIAAAtYEgrResyZgfy8KEPuG2iIEJaFXbflO3lpVLDl5RYbNw70ukj/0uFCbE4EQQAABBBBAAIFwCQRlvRYuP3VcCpBw6lc4dhAS+uDhUuk47j358acj5aNTT79ST8FiQwABBBBAAAEE/C4QhPVauOeAAiTcM3Dc8YOQ0Bt37ZfukxaVjyrm9Dry5gNJ0rhhXQ9J0xUEEEAAAQQQQMCaQBDWa9ZGrm8vChB9lrYjBSGhV27eK72nLTEsnrjhEunX+XyJrlnDtg0BEEAAAQQQQAABLwgEYb0WbkcKkHDPQMDOgCws+FYGzlhujCqzb3u59tJzPSRMVxBAAAEEEEAAAXsCFCD2/NTeFCD2DbVF8HtCf3+wRH797L/lm70HDZPZgzrJFS3O0uZDIAQQQAABBBBAINwCfl+vhduPAsQLMxCgMyDPLvxKnv6/9eUjeu/hq+SCRg08pkx3EEAAAQQQQAAB6wIUINbtyvbkDIh9Q20R/J7Qj766Wl5b8Y3h0euSGHn+7vYSFRWlzYdACCCAAAIIIIBAuAX8vl4Ltx9nQLwwAwE6AzL4pRUyb812aXBatOSPSfaYLt1BAAEEEEAAAQTsC1CA2DfkDIh9Q20R/J7QfbM/kY8Kd8mvGp4mS0b00OZCIAQQQAABBBBAwCsCfl+vecGRAsQLs/CfPvg9oX87NU/UW9AviqkvCx7q6iFZuoIAAggggAACCOgR8Pt6TY+CvSgUIPb8tO7t94Tu+fRiKdyxTxLOO0PmDEnSakMwBBBAAAEEEEDACwJ+X695wZACxAuz4PMzIEdKj8obnxXL8Nc+N0Zy5YVny0sDr/CQLF1BAAEEEEAAAQT0CFCA2HekALFvqC2CXxP63fxtMuTlleUOv27TWKb2SdDmQiAEEEAAAQQQQMArAn5dr3nFT/WDAsRDs+HXhJ6y8Cv5y3Hv/3jhnsuNx/CyIYAAAggggAACQRPw63rNS/NAAeKh2fBrQo9/t0CyPtxgSC4b2UManX6ah1TpCgIIIIAAAgggoE/Ar+s1fQL2I1GA2DfUFsGvCT3yjXz5xyebDYevMq6TWjVraDMhEAIIIIAAAggg4CUBv67XvGRIAeKh2fBrQg975TN5a9VWOa1WDVn75HUeEqUrCCCAAAIIIICAXgG/rtf0KtiLRgFiz0/r3n5N6IG5n8rCtTvk7Pq1ZfnjvbSaEAwBBBBAAAEEEPCSgF/Xa14ypADx0Gz4MaE/XL9T7nlxmaHY7KxfyKLh3T0kSlcQQAABBBBAAAG9An5cr+kVsB+NAsS+obYIfkvoL7f+INc/++/y8bf+1enyv3+4UpsHgRBAAAEEEEAAAa8J+G295jU/1R8KEA/Nit8S+q1VxTLslVXlgiOvv1gGXRXnIVG6ggACCCCAAAII6BXw23pN7+j1RKMA0eOoJYrfEvrlTzbJqDe+MMaeO6CDdGvZSIsDQRBAAAEEEEAAAa8K+G295kVHChAPzYrfEjpr8dcyfu5aQ3DhI10l7pz6HtKkKwgggAACCCCAgH4Bv63X9AvYj0gBYt9QWwS/JfRfFqyTKe8XGuPnBYTa0oBACCCAAAIIIOBhAb+t17xISQHioVnxW0Knv71GcvKKDME16clSr060hzTpCgIIIIAAAgggoF/Ab+s1/QL2I1KA2DfUFsFPCb1u+4+S/NcPjbFHRYlsGHe9RKn/YUMAAQQQQAABBAIs4Kf1mlengQLEQzPjl4Teu/+wdP7zQvmp5Kih16BOtOSnJ3tIkq4ggAACCCCAAALOCPhlvebM6PVEpQDR46glil8SeuXmvdJ72pLyMd/WvolMvK2dFgOCIIAAAggggAACXhbwy3rNy4YUIB6aHb8k9PFvPx+e3FKGdIvj8isP5RFdQQABBBBAAAHnBPyyXnNOwH5kChD7htoi+CWh5+Zvk/tfXmmMO2dAB+nO+z+05QCBEEAAAQQQQMDbAn5Zr3lZkQLEQ7Pjl4T+5/Itkvba54bcq4M7S4dmZ3pIka4ggAACCCCAAALOCfhlveacgP3IFCD2DbVF8EtC5+RtlPS3vzTGPXfYldKq8enaDAiEAAIIIIAAAgh4WcAv6zUvG1KAeGh2/JDQh46Uys1Tl8iX234w5P6d1l2anvkLDynSFQQQQAABBBBAwDkBP6zXnBu9nsgUIHoctUTxQ0JPmr9Onvvg57efq23lE73kzHq1tYyfIAgggAACCCCAgNcF/LBe87ohBYiHZsgPCd0/Z5ksWrfTUIs/7wyZc38iT8DyUA7RFQQQQAABBBBwVsAP6zVnBexHpwCxb6gtgh8S+rbMJfJp0V5jzIUZ10l0zRraxk8gBBBAAAEEEEDA6wJ+WK953ZACxEMz5IeEvvavH8ra7T9K2yYN5V+/7+IhPbqCAAIIIIAAAgg4L+CH9ZrzCvaOQAFiz0/r3n5I6C5PvS/f7D0onVucJbMGddI6foIhgAACCCCAAAJeF/DDes3rhhQgHpohPyR0/J8WyN4DJdLrkhh54Z7LPaRHVxBAAAEEEEAAAecF/LBec17B3hEoQOz5ad3bywl97NgxWbl5r9wyfakx5t7xsfL0HZdpHT/BEEAAAQQQQAABrwt4eb3mdbuy/lGAeGimvJzQb35WLA/OXlWudXen8+XJ317qIT26ggACCCCAAAIIOC/g5fWa86PXcwQKED2OWqJ4OaFHv/WFzFi6qXyc0/skyHVtGmsZN0EQQAABBBBAAAG/CHh5veYXQwoQGzNVWloqEydOlOzsbNmyZYs0bdpUUlJSZPjw4VKzZk3Tkb2c0A//c5XMWVlsjOmDR7tJ87PrmR4fOyCAAAIIIIAAAn4X8PJ6zS+2FCA2ZmrIkCEyffp0GTBggCQmJkpeXp7k5uaK+vrUqVNNR/ZyQqe+tFzmr/lWGtatJatHX2N6bOyAAAIIIIAAAggEQcDL6zW/+FKAWJyp/Px8adeunQwdOlQmT55cHmXYsGEyZcoUWb16tbRp08ZUdC8ndJ/sjyWvcLfEnlFX8h672tS4aIwAAggggAACCARFwMvrNb8YU4BYnKlRo0bJuHHjZMOGDdK8efPyKBs3bpQWLVqI+v7YsWNNRfdyQt80NU9Wb/lOWsY0kPkPXWVqXDRGAAEEEEAAAQSCIuDl9ZpfjClALM5UcnKycZZj+/btJ0WIiYmR+Ph4mTdvnqnobib0wcOlsnj9zhP6FxV1YneP/+uYf62Rrd//JO3P/6W8fn+iqXHRGAEEEEAAAQQQCIqAm+u1oJhVHAcFiMWZVZdX1a5dW1asWHFShISEBCkpKRF1mVZV27Zt20T9OX4rKCiQvn37GjFVDCe3b/YekC5PfWD6EF0vOkdm3NvR9H7sgAACCCCAAAIIBEGAAsT+LFKAWDSMi4sTdaZjyZIlJ0VQN6Tv2LFDCgsLq4w+ZswYSU9Pr/T7Xi5AHul1kQztcaFFNXZDAAEEEEAAAQT8LUABYn/+KEAsGvr9DMhPJaWydMPu/47+2PH/+9+/HDvu62f8opYknPdLiap4rZZFQ3ZDAAEEEEAAAQT8JkABYn/GKEAsGvr9HhCLw2Y3BBBAAAEEEEAgogUoQOxPPwWIRcORI0fK+PHjq3wKlvp+RkaGqegktCkuGiOAAAIIIIAAAq4LsF6zT04BYtFQPQFLPemqqveArFq1Stq2bWsqOgltiovGCCCAAAIIIICA6wKs1+yTU4DYMBw8eLBkZWUZb0JPSkoy3oSek5MjqampkpmZaToyCW2ajB0QQAABBBBAAAFXBViv2eemALFheOTIEZkwYYJkZ2dLcXGxxMbGSkpKiqSlpUl0dLTpyCS0aTJ2QAABBBBAAAEEXBVgvWafmwLEvqG2CCS0NkoCIYAAAggggAACjgiwXrPPSgFi31BbBBJaGyWBEEAAAQQQQAABRwRYr9lnpQCxb6gtAgmtjZJACCCAAAIIIICAIwKs1+yzUoDYN9QWgYTWRkkgBBBAAAEEEEDAEQHWa/ZZKUDsG2qLQEJroyQQAggggAACCCDgiADrNfusFCD2DbVFIKG1URIIAQQQQAABBBBwRID1mn1WChD7htoikNDaKAmEAAIIIIAAAgg4IsB6zT4rBYh9Q20RSGhtlARCAAEEEEAAAQQcEWC9Zp+VAsS+obYIJLQ2SgIhgAACCCCAAAKOCLBes89KAWLfUFsEElobJYEQQAABBBBAAAFHBFiv2WelALFvqC0CCa2NkkAIIIAAAggggIAjAqzX7LNSgNg31BaBhNZGSSAEEEAAAQQQQMARAdZr9lkpQOwbaotAQmujJBACCCCAAAIIIOCIAOs1+6wUIPYNtUXIy8uTLl26yMyZM6VVq1ba4hIIAQQQQAABBBBAQI9AQUGB9O3bVz766CNJSkrSEzTColCAeGjCX375ZSOh2RBAAAEEEEAAAQS8LaB+YdynTx9vd9KjvaMA8dDE7Nq1S+bPny/NmjWTunXrOt6zsgqeMy6OU59wANzd9VZHwxxz9wXCc0Ry3X13zCPP/ODBg1JUVCTJycly9tlnuw8QgCNSgARgEq0OgWsYrcrZ2w93e35W9sbcipq9fTC352d1b9ytylnfD3Prdlb3xNyqnHf2owDxzly43hP+AbtObhwQd/fdMcfcfYHwHJFcd98dc8zdF/D/ESlA/D+HlkfAh6ZlOls74m6Lz9LOmFtis7UT5rb4LO+Mu2U6yztibpnO8o6YW6bzzI4UIJ6ZCvc7wj9g9805A4J5eATcPyqfL+6b8/mCeXgE3D8qny/um+s+IgWIblEfxdu2bZtkZWVJamqqNG7c2Ec993dXcXd//jDH3H2B8ByRXHffHXPM3Rfw/xEpQPw/h4wAAQQQQAABBBBAAAHfCFCA+Gaq6CgCCCCAAAIIIIAAAv4XoADx/xwyAgQQQAABBBBAAAEEfCNAAeKbqaKjCCCAAAIIIIAAAgj4X4ACxP9zyAgQQAABBBBAAAEEEPCNAAWIb6aKjiKAAAIIIIAAAggg4H8BChD/z6HpEZSWlsrEiRMlOztbtmzZIk2bNpWUlBQZPny41KxZ03S8SNyhqKhImjdvXunQBw4caNiWbWa8zbQNuvu+fftk0qRJsmLFClm+fLls375d+vXrJ7m5uScN3YybU22DMB+hmpvJf+WCedXZoXJ75syZ8v7778vGjRulXr160rp1axkxYoT07NnzhB2dcjQTNwh5Hqo5ea53tgsKCiQ9Pd34TFePLq5Ro4bExcXJgAEDZPDgwVK7dm3Hf25GWq7rnUG90ShA9Hr6ItqQIUNk+vTpxj/6xMREycvLMxZ16utTp071xRjC3cmyH0w33XST3HrrrSd054ILLpBOnTqVf82Mt5m24TZw+vhlxuodNe3bt5d33nmnygLEjJtTbZ32cCN+qOZm8l/1G/OqZ099fixevFhuueUWSUhIEFUE5uTkyBdffCHTpk2T+++/3/HPEjPz40YeOn2MUM3Jc70zsWDBAuOXSldccYU0adLE+MWEWn/MmjVLbrzxRnnzzTfJdb3kno5GAeLp6dHfufz8fGnXrp0MHTpUJk+eXH6AYcOGyZQpU2T16tXSpk0b/QcOWMSyH0yjRo2SsWPHVjk6M95m2gaMs9LhHDp0SHbt2iWxsbFy5MgRqVWrVqUFiBk3p9oGZT5CNQ81/5UL5tVnh1qAXX755VKnTp3yhgcPHpTLLrtMdu7cKTt27JDo6GjHHM3MT1DyPFRz8tydGf/9739v/PJz7dq10rJlS3LdHfawH4UCJOxT4G4H1IJ53LhxsmHDhhMuIVKn/lu0aCGnWlC721vvHu34H0zKTG1169Y9qcNmvM209a6MMz2rrgAx4+ZUW2dGHd6o1ZmHmv9qBJhbm8dHHnlEnn76aeMyWfXbYqcczcS1NhL/7FXRnDx3Z+7UWRF1CfjHH39snB0xk5NOtXVn5JF9FAqQCJv/5ORk4yyHup6+4hYTEyPx8fEyb968CFMxP9yyH0z169c3LplQm7r06qGHHjIuNynbzHibaWu+x/7eo7rFsBk3p9r6W7fy3odSgJwq/1VkzK1lx5133imvvfaa7N27V5SzU45m4lobiX/2qmge6uc8eW5ujg8cOCDqz/79+2XZsmXywAMPGGf5vv76a+MXeWZy0qm25kZEaysCFCBW1Hy8j7q8St3opW4Cq7ip649LSkqM059s1Qts3rxZ7r33Xrn55pvlvPPOk61bt8rzzz8vK1eulEcffdS4yV9tZrzNtI20+aluMWzGzam2QZyP6sxDzX/+DVjLDHWzrroE64YbbpDXX3/d0c8SM/8mrI3GH3tVZk6eOzN3Y8aMMW5GL9s6dOhg/PxUOe/kZwa57sx8Wo1KAWJVzqf7qSdOqDMdS5YsOWkE6oZ0db1xYWGhT0cX3m6rG+q6du0qS5culfXr1xtP9zDjbaZteEfq/tGrWwybcXOqrfsizh+xOvPKjl5Z/qt2mJubq++//146d+5s/FJDna0+//zzjQBOOZqJa24k/mldlTl57swcqkvA1Z/du3cbT39TD1zIyMiQbt26kevOkHsyKgWIJ6fFuU7xGwDnbFXkt99+23iaR1ZWlgwaNIgzIJq4OQOiCdJEGLMFSGX57+RvM00MxTdN1c3n6pISdVnK3LlzpXv37uV9N/PZ7VRb30Ca6Gh15lWFqfg5T56bAK+k6TPPPCN//OMfjYK7VatWjv3cNPPvwt6I2DsUAQqQUJQC1MbM9ZIBGrZrQ/n888+Np4yp3+aMHDnSsWtZXRuQRw7EPSDuT4SVAqRi/qtem/nMMdPWfRFnj3j48GHjlxcLFy407v1Qj/g+fjNj41RbZwXcj34q86p6RJ7rnatvv/1Wzj33XONnpvrZ6VT+momrd4REq0yAAiTC8kL9Ax8/fnyVT8Eq+wCIMBZtw50zZ47xPH/1IkL1QkIz3mbaauuwTwJVtxg24+ZUW58wmuqmlQKkYv6rA2J+anZlrd5NoX6z/ve//1369Olz0k5OOZqJe+qR+KdFKOZVjYY81zvPZTf7q3feqHffmMlJp9rqHSHRKEDIAeMUp3rSVVXvAVm1apW0bdsWqVMIqHtlGjVqdEIrdSo/KSnJuJ5VPc1DvWHejLeZtpE2QdUths24OdU2iPNRnXmo+a9cMK8+O44ePSp33XWXzJ49WzIzMyU1NbXSHZxyNBM3KHkeqjl5rnfGK/NUR0hLSzMe3KJeiNyvXz/HPjMiMdf1zqDeaJwB0evpi2iDBw827lFQb0JXC2b1Uib15l31g0/9AGQ7tYB6+pW6ge7qq682ntGvbhidMWOGcWZJnWF67LHHyoOY8TbT9tS99H+L5557Tr777jtRC4bRo0cbxXPv3r2NganLVcqKZTNuTrX1v/bPIwjF3Ez+q5iYV50dDz/8sKhr4K+66iq57777TmrYq1cv48EhTjqamZ8g5Hmo5uS53tku81Q3m6tf0KnP9vnz5xuXHXbp0kU++OAD43G85Lped69GowDx6sw42C/1m80JEyYYlwkVFxcbb5pOSUkxfgtR9o/fwcMHIvTf/vY3o+BYt26d7Nmzx3hOv3qMsXqjvFoYH7+Z8TbTNhCQpxhEs2bNZNOmTZW2UkVz//79je+ZcXOqbVDmIxRzM/nv5PwEwVwtxhYvXlzlUNSirOzpQE7lrpm4kWROnuudbXWWT53lUPfQ7Ny5U+rUqSMXX3yx3HHHHcZVGervZZuZnHSqrd7RE62iAAUIOYEAAggggAACCCCAAAKuCVCAuEbNgRBAAAEEEEAAAQQQQIAChBxAAAEEEEAAAQQQQAAB1wQoQFyj5kAIIIAAAggggAACCCBAAUIOIIAAAggggAACCCCAgGsCFCCuUXMgBBBAAAEEEEAAAQQQoAAhBxBAAAEEEEAAAQQQQMA1AQoQ16g5EAIIIIAAAggggAACCFCAkAMIIIAAAggggAACCCDgmgAFiGvUHAgBBBBAAAEEEEAAAQQoQMgBBBBAAAEEEEAAAQQQcE2AAsQ1ag6EAAIIBEegqKhImjdvLqNHj5YxY8YEZ2CMBAEEEEDAcQEKEMeJOQACCCDgD4GoqKiQOtq1a1fJzc2lAAlJi0YIIIAAAhUFKEDICQQQQAABQ2DmzJknSMyZM0feeOMNmTRpksTExJR/T/1/z5495dChQxIdHW38YUMAAQQQQCBUAQqQUKVohwACCESYgLq0Kj09Xb766iu54IILImz0DBcBBBBAwCkBChCnZImLAAII+FygugKksntAFi1aJN27d5cXXnhB9u/fL88++6xs3bpVEhISZPr06dK2bVt55ZVXZOzYsVJYWChxcXEyZcoUufrqq0+QKi0tNfbNycmR9evXS926daVbt24ybtw4adWqlc9V6T4CCCCAAAUIOYAAAgggUKmA1QIkPj5eDh8+LAMGDJADBw7IU089JQ0bNpSMjAzjjEpqaqpx2Zb6eklJiWzatMn4ftl22223yZtvvil33323dOzYUXbv3i3Tpk0zippPP/1ULrzwQmYMAQQQQMDHAhQgPp48uo4AAgg4KWC1AGnSpIkUFBRI/fr1je5NnjxZHnzwQaPIWLduXfn9JK+++qrcfvvtkpWVJYMGDTLaln1N/ffWW28tH15xcbFccsklcv3118usWbOcHDaxEUAAAQQcFqAAcRiY8AgggIBfBawWICNGjDAulyrbli9fLh06dJB77rlHZsyYUf71b775Rpo2bSppaWnG2RC19e7dWz755BNZvXr1SWx9+vSRlStXys6dO/1KSr8RQAABBESEAoQ0QAABBBCoVMBqAZKZmWlcZlW2qfs91GVTjz/+uDz55JPlX9+3b580aNBAUlJSjPtG1KbOcqhbrRJhAAACW0lEQVSzJ9Vt6h6RGjVqMGsIIIAAAj4VoADx6cTRbQQQQMBpAasFiComVFFRsQCp+NLCsgJk4MCBkp2dbTRXN5mrx/s+//zzVQ6vR48eEuo7S5w2Ij4CCCCAgHkBChDzZuyBAAIIRIRAOAqQ3/zmN6KeprVnzx6pVatWRDgzSAQQQCDSBChAIm3GGS8CCCAQokA4ChB1g/ldd90lFc+WlHV5x44d0qhRoxBHQDMEEEAAAS8KUIB4cVboEwIIIOABgXAUIMeOHZM777xTZs+eLddcc4306tXLuE9EPap37ty50rp165Pe2O4BKrqAAAIIIGBCgALEBBZNEUAAgUgSCEcBonxVEaLuAXnxxRdlzZo1xt9jY2MlKSlJ7rvvPklMTIykaWCsCCCAQOAEKEACN6UMCAEEEEAAAQQQQAAB7wpQgHh3bugZAggggAACCCCAAAKBE6AACdyUMiAEEEAAAQQQQAABBLwrQAHi3bmhZwgggAACCCCAAAIIBE6AAiRwU8qAEEAAAQQQQAABBBDwrgAFiHfnhp4hgAACCCCAAAIIIBA4AQqQwE0pA0IAAQQQQAABBBBAwLsCFCDenRt6hgACCCCAAAIIIIBA4AQoQAI3pQwIAQQQQAABBBBAAAHvClCAeHdu6BkCCCCAAAIIIIAAAoEToAAJ3JQyIAQQQAABBBBAAAEEvCtAAeLduaFnCCCAAAIIIIAAAggEToACJHBTyoAQQAABBBBAAAEEEPCuAAWId+eGniGAAAIIIIAAAgggEDgBCpDATSkDQgABBBBAAAEEEEDAuwL/DwWiufXKuaHRAAAAAElFTkSuQmCC\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.lineplot(x='Time', y='Mass', data=mass_df);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "4023c188",
   "metadata": {
    "pycharm": {
     "name": "#%%\n"
    }
   },
   "outputs": [],
   "source": [
    "mass_df.to_csv('mass_data/MassProfile_050923_S817-822.csv')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "90ac9fc0",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}