examples/weather/vismet.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 54,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Requirements\n",
    "\n",
    "#!pip3 install numpy\n",
    "#!pip3 install pandas\n",
    "#!pip3 install geopandas\n",
    "#!pip3 install matplotlib\n",
    "#!pip3 install scipy\n",
    "#!pip3 install networkx\n",
    "#!pip3 install descartes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "# Libraries used\n",
    "\n",
    "import pandas as pd\n",
    "import geopandas as gpd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import networkx as nx\n",
    "from scipy.spatial import distance_matrix"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Constants and Configurations\n",
    "\n",
    "FILEPATH_DATA = 'data/mediciones.csv'\n",
    "FILEPATH_DDBB = 'data/estaciones.csv'\n",
    "\n",
    "# Plot in Dark Mode like a Hacker\n",
    "plt.style.use('dark_background')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Utils Functions\n",
    "\n",
    "def delete_nan(row):\n",
    "    try:\n",
    "        row[[\"Latitud\", \"Longitud\"]] = pd.to_numeric(row[[\"Latitud\", \"Longitud\"]])\n",
    "        return row\n",
    "    except:\n",
    "        return np.nan"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Script\n",
    "\n",
    "coords = pd.read_csv(FILEPATH_DDBB).apply(delete_nan, axis=1).dropna()\n",
    "D = distance_matrix(coords[[\"Latitud\", \"Longitud\"]], coords[[\"Latitud\", \"Longitud\"]])\n",
    "\n",
    "D[D == 0] = 0.0001 # Deal with 1/0\n",
    "D = 1/(D**2)\n",
    "\n",
    "graph = nx.convert_matrix.from_numpy_matrix(D)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "metadata": {},
   "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_dpi_ratio', { dpi_ratio: fig.ratio });\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(\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",
       "    );\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('style', 'box-sizing: content-box;');\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",
       "    if (this.ratio !== 1) {\n",
       "        fig.send_message('set_dpi_ratio', { dpi_ratio: this.ratio });\n",
       "    }\n",
       "\n",
       "    var rubberband_canvas = (this.rubberband_canvas = document.createElement(\n",
       "        'canvas'\n",
       "    ));\n",
       "    rubberband_canvas.setAttribute(\n",
       "        'style',\n",
       "        'box-sizing: content-box; position: absolute; left: 0; top: 0; 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",
       "            canvas.setAttribute(\n",
       "                'style',\n",
       "                'width: ' + width + 'px; height: ' + height + 'px;'\n",
       "            );\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 (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",
       "        return function (event) {\n",
       "            return fig.mouse_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousedown',\n",
       "        on_mouse_event_closure('button_press')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseup',\n",
       "        on_mouse_event_closure('button_release')\n",
       "    );\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousemove',\n",
       "        on_mouse_event_closure('motion_notify')\n",
       "    );\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseenter',\n",
       "        on_mouse_event_closure('figure_enter')\n",
       "    );\n",
       "    rubberband_canvas.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",
       "    this.rubberband_canvas.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",
       "    var cursor = msg['cursor'];\n",
       "    switch (cursor) {\n",
       "        case 0:\n",
       "            cursor = 'pointer';\n",
       "            break;\n",
       "        case 1:\n",
       "            cursor = 'default';\n",
       "            break;\n",
       "        case 2:\n",
       "            cursor = 'crosshair';\n",
       "            break;\n",
       "        case 3:\n",
       "            cursor = 'move';\n",
       "            break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = 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",
       "            /* 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",
       "            evt.data.type = 'image/png';\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",
       "                evt.data\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",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function (e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e) {\n",
       "        e = window.event;\n",
       "    }\n",
       "    if (e.target) {\n",
       "        targ = e.target;\n",
       "    } else if (e.srcElement) {\n",
       "        targ = e.srcElement;\n",
       "    }\n",
       "    if (targ.nodeType === 3) {\n",
       "        // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "    }\n",
       "\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    var boundingRect = targ.getBoundingClientRect();\n",
       "    var x = e.pageX - (boundingRect.left + document.body.scrollLeft);\n",
       "    var y = e.pageY - (boundingRect.top + document.body.scrollTop);\n",
       "\n",
       "    return { x: x, y: y };\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://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",
       "    var canvas_pos = mpl.findpos(event);\n",
       "\n",
       "    if (name === 'button_press') {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * this.ratio;\n",
       "    var y = canvas_pos.y * 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",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\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.which === this._key) {\n",
       "            return;\n",
       "        } else {\n",
       "            this._key = event.which;\n",
       "        }\n",
       "    }\n",
       "    if (name === 'key_release') {\n",
       "        this._key = null;\n",
       "    }\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which !== 17) {\n",
       "        value += 'ctrl+';\n",
       "    }\n",
       "    if (event.altKey && event.which !== 18) {\n",
       "        value += 'alt+';\n",
       "    }\n",
       "    if (event.shiftKey && event.which !== 16) {\n",
       "        value += 'shift+';\n",
       "    }\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\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 icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\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.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",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['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",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager) {\n",
       "        manager = IPython.keyboard_manager;\n",
       "    }\n",
       "\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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4Xu2dCbyP1fb/lyFS/TWgDNHsIs3cWzdCgwxRhoSraB400CgpStI84VZKEUXS3C1KhkqhUaWkMoQIhQyR8H+tXd/zO45zzveZv8/z7Pd+vbwq9rTe62mdj7332ruEiGwTCgQgAAEIQAACEICANQRKIACt8TWGQgACEIAABCAAAUMAAciHAAEIQAACEIAABCwjgAC0zOGYCwEIQAACEIAABBCAfAMQgAAEIAABCEDAMgIIQMscjrkQgAAEIAABCEAAAcg3AAEIQAACEIAABCwjgAC0zOGYCwEIQAACEIAABBCAfAMQgAAEIAABCEDAMgIIQMscjrkQgAAEIAABCEAAAcg3AAEIQAACEIAABCwjgAC0zOGYCwEIQAACEIAABBCAfAMQgAAEIAABCEDAMgIIQMscjrkQgAAEIAABCEAAAcg3AAEIQAACEIAABCwjgAC0zOGYCwEIQAACEIAABBCAfAMQgAAEIAABCEDAMgIIQMscjrkQgAAEIAABCEAAAcg3AAEIQAACEIAABCwjgAC0zOGYCwEIQAACEIAABBCAfAMQgAAEIAABCEDAMgIIQMscjrkQgAAEIAABCEAAAcg3AAEIQAACEIAABCwjgAC0zOGYCwEIQAACEIAABBCAfAMQgAAEIAABCEDAMgIIQMscjrkQgAAEIAABCEAAAcg3AAEIQAACEIAABCwjgAC0zOGYCwEIQAACEIAABBCAfAMQgAAEIAABCEDAMgIIQMscjrkQgAAEIAABCEAAAcg3AAEIQAACEIAABCwjgAC0zOGYCwEIQAACEIAABBCAfAMQgAAEIAABCEDAMgIIQMscjrkQgAAEIAABCEAAAcg3AAEIQAACEIAABCwjgAC0zOGYCwEIQAACEIAABBCAfAMQgAAEIAABCEDAMgIIQMscjrkQgAAEIAABCEAAAcg3AAEIQAACEIAABCwjgAC0zOGYCwEIQAACEIAABBCAfAMQgAAEIAABCEDAMgIIQMscjrkQgAAEIAABCEAAAcg3AAEIQAACEIAABCwjgAC0zOGYCwEIQAACEIAABBCAfAMQgAAEIAABCEDAMgIIQMscjrkQgAAEIAABCEAAAcg3AAEIQAACEIAABCwjgAC0zOGYCwEIQAACEIAABBCAfAMQgAAEIAABCEDAMgIIQMscjrkQgAAEIAABCEAAAcg3AAEIQAACEIAABCwjgAC0zOGYCwEIQAACEIAABBCAfAMQgAAEIAABCEDAMgIIQMscjrkQgAAEIAABCEAAAcg3AAEIQAACEIAABCwjgAC0zOGYCwEIQAACEIAABBCAfAMQgAAEIAABCEDAMgIIQMscngtzq1atKmvXrs3F0IwJAQhAwEoC/+///T/56aefrLQdo50RQAA640QtjwRU/C1ZssRja5pBAAIQgIBXAtWqVUMEeoVnQTsEoAVOzqWJ+rfQ3377TTQQsQqYS08wNgQgYAsBjbv6F+/y5csTd21xugc7EYAeoNHEOYGMACQQOWdGTQhAAAJ+CBB3/dCzpy0C0B5f58RSAlFOsDMoBCBgMQHirsXOd2E6AtAFLKq6J0Agcs+MFhCAAAT8ECDu+qFnT1sEoD2+zomlBKKcYGdQCEDAYgLEXYud78J0BKALWFR1T4BA5J4ZLSAAAQj4IUDc9UPPnrYIQHt8nRNLCUQ5wc6gEICAxQSIuxY734XpCEAXsKjqngCByD0zWkAAAhDwQ4C464eePW0RgPb4OieWEohygp1BIQABiwkQdy12vgvTEYAuYFHVPQECkXtmtIAABCDghwBx1w89e9oiAO3xdU4sJRDlBHsiBi0pIg1FpIqILBWR90RkayJmziQhEG8CxN14+ycus0MAxsUTKZ0HgSiljvVpVhsReUhEqufrZ5GIXCUiL/nsm+YQsJ0Acdf2L8CZ/QhAZ5yo5ZEAgcgjuBQ3U/E37m/7dBUwUzKrf+0RgSn2PqZFQYC4GwXl5I+BAEy+D2NtAYEo1u6JfHIq+BaISDURyS/+8ovAxSJyANvBkfuGAdNDgLibHl+GaQkCMEy69C0EIj6C/AQaicgUB0gai8hUB/WoAgEI7EiAuMtX4YQAAtAJJep4JkAg8owulQ07ishoB5Z1EpExDupRBQIQQADyDXgjgAD0xo1WDgkgAB2CsqQaK4CWOBozc0qAuJtT/IkZHAGYGFclc6IEomT6LaxZcwYwLLL0C4H/I0Dc5WtwQgAB6IQSdTwTIBB5RpfahmQBp9a1GBYTAsTdmDgi5tNAAMbcQUmfHoEo6R4MZ/6F3QP4o4j04AqYcIDTq1UEiLtWuduzsQhAz+ho6IQAgcgJJTvr8BKInX7H6vAJEHfDZ5yGERCAafBijG0gEMXYOUwNAhBIJQHibirdGrhRCMDAkdJhfgIEIr4HCEAAAtESIO5GyzupoyEAk+q5hMybQJQQRzFNCEAgNQSIu6lxZaiGIABDxUvnBCK+AQhAAALREiDuRss7qaMhAJPquYTMm0CUEEcxTQhAIDUEiLupcWWohiAAQ8VL5wQivgEIQAAC0RIg7kbLO6mjIQCT6rmEzJtAlBBHMU0IQCA1BIi7qXFlqIYgAEPFS+cEIr4BCEAAAtESIO5GyzupoyEAk+q5hMybQJQQRzFNCEAgNQSIu6lxZaiGIABDxUvnBCK+AQhAAALREiDuRss7qaMhAJPquYTMm0CUEEcxTQhAIDUEiLupcWWohiAAQ8VL5wQivgEIQAAC0RIg7kbLO6mjIQCT6rmEzJtAlBBHMU0IQCA1BIi7qXFlqIYgAEPFS+cEIr4BCEAAAtESIO5GyzupoyEAk+q5hMybQJQQRzFNCEAgNQSIu6lxZaiGIABDxUvnBCK+AQhAAALREiDuRss7qaMhAJPquYTMm0CUEEcxTQhAIDUEiLupcWWohiAAQ8VL5wQivgEIQAAC0RIg7kbLO6mjIQCT6rmEzJtAlBBHMU0IQCA1BIi7qXFlqIYgAEPFS+cEIr4BCEAAAtESIO5GyzupoyEAk+q5hMybQJQQRzFNCEAgNQSIu6lxZaiGIABDxUvnBCK+AQhAAALREiDuRss7qaMhAJPquYTMm0CUEEcxTQhAIDUEiLupcWWohiAAQ8VL5wQivgEIQAAC0RIg7kbLO6mjIQCT6rmEzJtAlBBHMU0IQCA1BIi7qXFlqIYgAEPFS+cEIr4BCEAAAtESIO5GyzupoyEAk+q5hMybQJQQRzFNCEAgNQSIu6lxZaiGIABDxUvnBCK+AQhAAALREiDuRss7qaMhAJPquYTMm0CUEEcxTQhAIDUEiLupcWWohiAAQ8VL5wQivgEIQAAC0RIg7kbLO6mjIQCT6rmEzJtAlBBHMU0IQCA1BIi7qXFlqIYgAEPFS+cEIr4BCEAAAtESIO5GyzupoyEAk+q5hMybQJQQRzFNCEAgNQSIu6lxZaiGIABDxUvnBCK+AQhAAALREiDuRss7qaMhAJPquYTMm0CUEEcxTQhAIDUEiLupcWWohiAAQ8Ub78579eolbdu2lVq1asnvv/8uH3zwgdxwww0yd+7cvImXLVtW7rvvPunYsaPov0+YMEEuu+wyWb58uSPjCESOMFEJAhCAQGAEiLuBoUx1RwjAVLu3eOPefPNNGTNmjHz00UdSunRpueOOO6Ru3bpSp04d2bBhg2n83//+V1q2bCndunWTNWvWyODBg2Xr1q3SoEEDR+QIRI4wUQkCEIBAYASIu4GhTHVHCMBUu9edcRUrVpQVK1bICSecIO+9956UL1/e/Hfnzp3lhRdeMJ394x//kDlz5sixxx4rM2bMyDoAgSgrIipAAAIQCJQAcTdQnKntDAGYWte6N+yggw6S77//3qwCzp49W5o0aSKTJk2SPfbYw6z+ZcqCBQvkwQcfNL+yFQJRNkL8OQQgAIFgCRB3g+WZ1t4QgGn1rEu7SpQoIa+++qoRew0bNjStO3XqJE899ZTsvPPO2/WmK3+TJ08WPUNYsJQpU8acFcwUDURLliwxq4lr1651OSuqQwACEICAWwIIQLfE7KyPALTT7ztYrWf9mjdvbs72qWDzKgD79u0r/fr126F/BCAfGgQgAIFoCCAAo+Gc9FEQgEn3YADzHzRokJx++unm7J9u72aKly1gVgADcAhdQAACEPBBAAHoA55FTRGAFjm7MFNV/LVp00YaN25szv/lL5kkEN0KfvHFF80f1axZU7799luSQCz/bjAfAhCILwEEYHx9E6eZIQDj5I2I5zJkyBCT4aurfyrqMkUTPjZu3Gj+U7eGW7RoYa6B+e2330QFo5bjjz/e0WwJRI4wUQkCEIBAYASIu4GhTHVHCMBUu7d447Zt21ZoBRV7I0aMMH+WuQhaVwHzXwT9888/OyJHIHKEiUoQgAAEAiNA3A0MZao7QgCm2r25N45AlHsfMAMIQMAuAsRdu/zt1VoEoFdytHNEgEDkCBOVIAABCARGgLgbGMpUd4QATLV7c28cgSj3PmAGEICAXQSIu3b526u1CECv5GjniACByBEmKkEAAhAIjABxNzCUqe4IAZhq9+beOAJR7n3ADCAAAbsIEHft8rdXaxGAXsnRzhEBApEjTFSCAAQgEBgB4m5gKFPdEQIw1e7NvXEEotz7gBlAAAJ2ESDu2uVvr9YiAL2So50jAgQiR5ioBAEIQCAwAsTdwFCmuiMEYKrdm3vjCES59wEzgAAE7CJA3LXL316tRQB6JUc7RwQIRI4wUQkCEIBAYASIu4GhTHVHCMBUuzf3xhGIcu8DZgABCNhFgLhrl7+9WosA9EqOdo4IEIgcYaISBCAAgcAIEHcDQ5nqjhCAqXZv7o0jEOXeB8wAAhCwiwBx1y5/e7UWAeiVHO0cESAQOcJEJQhAAAKBESDuBoYy1R0hAFPt3twbRyDKvQ+YAQQgYBcB4q5d/vZqLQLQKznaOSJAIHKEiUoQgAAEAiNA3A0MZao7QgCm2r25N45AlHsfMAMIQMAuAsRdu/zt1VoEoFdytHNEgEDkCBOVIAABCARGgLgbGMpUd4QATLV7c28cgSj3PojjDA455BD57rvvtptaiRIlZNu2bXGcLnOCQKIIEHcT5a6cTRYBmDP0dgxMILLDz06sLFWqlGzZskVOPvlkUbH39ttv5zXL/FnJkiVl69atTrqjDgQgUAQB4i6fhhMCCEAnlKjjmQCByDO6VDVUwXfggQfKTjvtJHPmzNnBNv3zY489Vj788MNU2Y0xEMgFAeJuLqgnb0wEYPJ8lqgZE4gS5a5QJlu6dGnp3bu3tG7dWurVq7fDGJnVv759+8rq1avl888/l19++UW++uqrUOZDpxBIOwHibto9HIx9CMBgONILWxHWfwO6fas/ePbYYw/RM36XX365LF26VE444QSZNWuWdO3aVcqXLy8tW7aUp59+2vCqXbu2dOzYUdq1ayd16tQxW8Na9Czgrbfean5RIAABdwQQgO542VobAWir5yOym0AUEegYDNOhQwd56qmnZJdddtluNgsXLpSaNWsasTd+/Hj54IMPjOBTobhq1apiZ96iRQuZMGGC6Cphpvz5559GIGZWDmNgOlOAQKwIEHdj5Y7YTgYBGFvXpGNiBKJ0+DGbFbvttpusWbNGdBWwYLniiitk8ODBMm/ePKlRo4Y0adJE3nvvPSlXrpysXLlyB8GYv/306dPlxRdfNIkhmkCy++67S5UqVWSfffaRJ554Qv73v//tMN6RRx4pw4YNM6uRu+66q9x8883y5JNPZjOBP4dAaggQd1PjylANQQCGipfOCUTp/wZ02/akk07aLqs3v9XHH3+8/PDDD7Js2TK577775LrrrjMreCoWO3fuLCNHjvQESUXhww8/LK+++qo5L/jHH39Inz59RAVn2bJl8/rUVca6deua8TOrhmQae0JOo4QQIO4mxFE5niYCMMcOSPvwBKK0e1hEkzyuuuoquffee3cwVoWWrto1btxYXnvtNdFEjwEDBpjVPBWO9evXlxkzZgQCSUVl5gxhwQ5VgL7wwgvy888/my1oXVmkQCCtBIi7afVssHYhAIPlSW8FCBCI0v1JZC5v1i3X8847bwdj7777brnhhhukZ8+ectttt8mGDRukevXqZrVOVwBVDPbq1StSSCoUdb4qWnU+FAikjQBxN20eDcceBGA4XOn1bwIEonR/Cueff77cfvvtUrly5UIN1fN3AwcONCt+ehWMCr6rr75aHnjgAZME8sknn5j7AXNRLr74YnM2UJNKKBBIEwHibpq8GZ4tCMDw2NKziDmI/9tvv5nrP9auXQuThBPIZOOqoNOVPRV3RW275jf1gAMOMPX0jj9NzPj111+lYsWKjtqGhWzo0KFyySWXmHOBiMCwKNNvLggQd3NBPXljIgCT57NEzZhAlCh3OZrslVdeKQ0bNpT27ds7qq+Vxo0bZ+77063ge+65x3G7MCtqUohuP+u5QL1qRgUqySFhEqfvqAgQd6MinexxEIDJ9l/sZ08gir2LXE1QV8s0kUO3dt2WatWqyUEHHSTvvvuu26ah1derZLp37y4777yznHHGGfLMM88YIciKYGjI6TgCAsTdCCCnYAgEYAqcGGcTCERx9o67uWUSPh555BGzdeq26AXRugoYlzv5NBFFt6P13kDNENY7Cr/55hs55ZRTzKogItCth6kfFwLE3bh4It7zQADG2z+Jnx2BKPEuzDOgUqVK5nm3m266abuXOZxaqGLr999/N5dA692AWjZu3CibN282Z0WjLp9++qkcc8wxcuedd5rzjJny/fffy3HHHWfeJUYERu0VxguCAHE3CIrp7wMBmH4f59RCAlFO8Qc6uAo4fdNXt3G9FN36bdWqlckI1u9CBZYKQC377befnHPOOebtXydJJV7GL9hG7wY8+OCDzW9/9913ef+u/63X12jWss6VAoGkESDuJs1juZkvAjA33K0ZlUCUHlfrhc9HHXWUzJw507NRH330kVx66aXmzJ2uKOqTblqmTp0qc+bMMQki1157ref+nTY87LDDzOshWvTs33PPPSdlypTJa65JK2eeeabT7qgHgVgRIO7Gyh2xnQwCMLauScfECETp8GN+K/TNX73WJ+ii2676KsioUaPMO8H6nNuee+4pP/74o8ydO9esHup5vSCK3kGo2796HlFXHwsrS5cuNc/U5d8eDmJs+oBA2ASIu2ETTkf/CMB0+DG2VhCIYusazxPTM3wNGjTw3L6ohpqUoYkYmoBRsGTe8J04caJ5d9hPuf/++80qo55l7N+/f9auVAR269atyCtiMnPL2hEVIBARAeJuRKATPgwCMOEOjPv0CURx95D7+en2aLt27dw3dNBi+fLl8tBDD0mbNm1M7YULF0qdOnVk+PDh5lyebhP/4x//cNBT0VXOPfdc05+bbOZTTz1V3n77bbMqmTm3qCOoGNVXUPQ84ccff2ySRhCEvtxD4wAIEHcDgGhBFwhAC5ycSxMJRLmkH/zYem5vyZIlstNOOwXfeRE9quDae++9RV8T0SQUv0Uvez7iiCPk66+/liFDhkiFChWkUaNGZoyiSs2aNc25RX3JREXp7NmzzWrl4Ycfntdk/vz5RpxqVjMFArkkQNzNJf3kjI0ATI6vEjlTAlEi3VbkpDVD95dffjFn86Iq9913n9mynTJlihFqQRTdbt53331lxYoVprvq1asbQbjbbrtt172KT734Wt8u7ty5s1k5LKpoFvOzzz5L5nAQDqIPXwSIu77wWdMYAWiNq3NjKIEoN9zDHPWVV16R1q1bhznEdn2fd9558tRTT8mmTZu2y9T1M4Ft27aZa170DkBdDdSVRV3B23///bfrtnbt2mbb+dhjjxV9NeTQQw8tdFhd9cufRexnbrSFgF8CxF2/BO1ojwC0w885s5JAlDP0oQ2sV6Z06NAhtP4Ldty8eXNzTcz69et93RGook/P6Y0dO1ZGjx5ttrL11Q891/f666/L6aefLi+99NJ2Y+jzdbr1q2f8shX91tetW5etGn8OgdAJEHdDR5yKARCAqXBjfI0gEMXXN15nps+l1apVy2tz1+00+UOvYvnss8+KvAZGt3QbNmxo7ihs3LixHH300dKyZUuTraxXvbz22msyffp0+fXXX835xUMOOUQuuugiueqqq0TbXnjhhfL000+beqeddpqZo15UXbFiRenXr5/06dMn67x1JfHLL78UFZpFlZIlSxaZTZx1ACpAwCEB4q5DUJZXQwBa/gGEbT6BKGzC0favAkbPzmkiRFRFV9X03r5OnToZkVbYSyEq4jQBY8GCBWZauh2rv5f5pyZ66HauXvrcvXt3s6qXv7zxxhtG+N14440yYMAA80cbNmww47Zo0UJefvnlrOZqnUz2cmGV9SJtFaT66oiuPlIgEBYB4m5YZNPVLwIwXf6MnTUEoti5xNeEVHz17NlTNDEjyqKrdY8//rjZuu3YsWOhQ2t2r2br6iXSeqZPEzL07OCJJ55o7vwrrmTO/ukKXv53iXWF8OGHHzbCsGvXrlmvoNGsYH1hpLBVQJ2XnifULXRd0SxupTBKtoyVPgLE3fT5NAyLEIBhUKXPPAIEovR8DLqCpUkQus0adcJD5u4+3c7Vfy+q6OqkZvR++OGHrsDvsssuZhtY7yDMX1RU6tvHuuKpZxH/97//FdnvsmXLpGrVqkUKOxWSV1xxhemrYLKJq8lSGQJZCBB3+UScEEAAOqFEHc8ECESe0cWuoQpAvQevqEzYsCas7wNff/31JknjxRdfFN2GLqrUq1fPXFOjGb1Oi14+rW8SDxs2TDTjuLDSpEkT+eSTT2TlypVFit9XX33VzLGwUvD+RH3S7osvvmAV0KmTqOeKAHHXFS5rKyMArXV9NIYTiKLhHPYomjihd+Hp+bmoy/vvvy/NmjUzL3C8++678u9//7vIKfTo0UMGDx4se+21l7mr8Oabb5YuXboUO2VN9lCBdumll5rt3sKKXn2j5/tUiOrVMZoYopnDmiWsySn6rrCeAVSRrK+B5C/6MogmkWibTBk4cKD5PV1hpEAgaALE3aCJprM/BGA6/RobqwhEsXGF54loFu7VV19tnjjLVdFzdZrZq4Jr3rx5RV4H8/3334uKK90q1rLffvvlJYYUN/eTTz7ZiEvNEi54GbS20ytodHVv0qRJ5kLqtm3bmroqNHVOek2MJp2o+Mt/tk+ZHXjggWZLWhNRMuXbb7+NNJM6V35j3NwQIO7mhnvSRkUAJs1jCZsvgShhDitkurrtWlx2a1AWaoKEiisVTfl/qahSAbpq1SozlG5D65UrRZWjjjrK1FFx1r9/f7nsssuyTlGzizXJQ9/7VTGYKbpCp28f67lDzQrW7Wdd5VOxl/m2M3VV5OkKol5Ho6uWH3zwgeiWtD43p1vMBYtmLatw3LJlS9b5UQECbggQd93QsrcuAtBe30diOYEoEsyhDlLc/XtBDnzmmWcasVVY0VU5FUoqvPRpNt3aLVh+++03GTlypFx++eVy5ZVX7pDQUdxcta0mcOgYjz76qKiIvPfee2Xy5MmyePFis8rXtGlT0ZU73Q7XhA+9DPu2227zjEBXNXVbPbNa6bkjGkKgAAHiLp+EEwIIQCeUqGNWUa677jqpXLmyuWJDsxk/+uijrGQIRFkRxbKCXveiW5n6zzp16sgLL7yQ9QoUv4aoANQsW11Fy1z1ooJPf+mFzJqlO2LECJPlq5nImVU1PYs3aNAgsz2rq4Uq5PQ6F10BdFM0GeS4444zW8xaMrbrt3/JJZcUm3ziZpz8ddU2fYmEAoEgCRB3g6SZ3r4QgOn1bWCW6UqHbpHpD8EZM2aIHrTXH9a6hbVixYpixyEQBeaGSDrKXLKsCRR69m/NmjVmxe3333+XGjVqyDPPPGNEUhjljjvuyLuv76STTpKJEyduN8yPP/5oVuZ0tU7P9ulderoSp6tweseertDpX1L+9a9/eZ6ebvnqXYP6lxx9UUQvgQ6z6IrlOeecE+YQ9G0hAeKuhU73YDIC0AM025roE1q62qerfpmVkUWLFplVl7vuugsBmKIPQu/3061PPQ9Xvnx5Y5mKfBVZml2rq256wbK+aPH8888bIab/7rfotS36dJs+M6dFxdyECRMK7fadd94xr3ZoVrBux6pY1dc1MvP1O5eo2utqZ+vWrckEjgq4ReMgAC1ytg9TEYA+4NnQVH/A6uH39u3bi16FkSnDhw83z2Tp01rFFQJRsr4SXUnTjNfCnlsbO3as2ZrVreHMm7b6z1atWsmYMWN2eF7NjeWZZ9S6detmtm51tVnf6y2q6NaprtYVfNLNzZi5rKvZynXr1jXb27wIkktPpHNs4m46/Rq0VQjAoImmrL8qVarITz/9ZLb9dCUwU3Tlr1GjRnLsscduZ7GuIOU/06SBSN891dWZtWvXpoxOusxRMXfttdcWu6qr38Ls2bPNL311Q7Nd9fWMXXfd1awGagKGirfCBGRxtF566SVztYomRkR90XQuvKhb17qdraurKgC5DzAXXkjvmAjA9Po2SMsQgEHSTGFfbgWgnhfLf+FtBgkCMP4fh95Tp1ux+d/CdTJrfSFDrzvJFN2eHTVqlOy+++5Omps6KihVEOmTbHrcoLjVP8edxrziggULjHhWUb106VLzTxWEesaRAgE/BBCAfujZ0xYBaI+vPVnqdguYFUBPmHPeSO+2K+4ljGwTVNGnmbmaqavJI5o4pAkjboreo6d38KkQ1JVFW4smV82dO9fchajX0mT+aSsP7HZPAAHonpmNLRCANnrdpc269as/3PVuNS26vacZmZoUQBKIS5gxrJ55vmzq1KlywgkneJqhng3U61jWrVtnhKTepXfxxReb5JG9997bcZ+dOnUyV8HYvAqm7wRrFrKuyCpHvSBb3xn++OOPtzt/6RgqFa0jgAC0zuWeDEYAesJmVyNdzdH71zaeRm4AACAASURBVPQHugpBvQZGf69WrVqid6cVVwhEyfhWdCsyI/Ddzlj/MqDPnTVv3lxee+01c55Nk4b03/X39VyfriQ7Kfpyx+233y6bNm1yUj2VdTQbWpNENAt6//33z7Pxl19+MYlYTz31lFkhzSTipBICRvkiQNz1hc+axghAa1ztz9Du3bvnXQStz2ypWFAxmK0QiLIRyv2f6wrgLbfcUujrGtlmp/cEnnrqqaLnAPVqIL0oPFOee+45kzWsT6E5eY5N2/Xs2VOGDh1qMpFtLZp1rz7R4xRFFT0/qOcs9RURPbdJgUB+AsRdvgcnBBCATihRxzMBApFndJE11JUk3XpVQbF69WpzvY+ToitSemZPV/hUtOl7uQVLzZo1TWKDipRMv7q9qVnDBx98sLkGRVcHM1nDKia/+OILkxRByU5AM4g7d+4sekUPmcTZedlSg7hri6f92YkA9MeP1lkIEIiS8Ynoxc+1a9eWPn36mOt+jjnmmGIn/vPPP0uTJk3MBcwqHM8666xC63/66adSv35984qIXhmkma7vvvuuEXx6NlCPEOi4KiCPP/5484ybXiejCSEUZwR0tVTFtT6DR0k+gZIi0lBEqoiI/jXoPRHZ6tIs4q5LYJZWRwBa6viozCYQRUU6uHF0RamwDF7NSH3vvfdk3Lhx5rk0vddR/11fsyiu3H///TJs2DAj/rToSqGuHqp41DNumjyifZ133nlmS1PPmmryA8UZgYLX8Ohbyeor5a1CW1d49b8p8SfQRkQeEpHq+aa6SESuEpGXXEyfuOsClsVVEYAWOz8K0wlEUVAOdowBAwZI7969t+tUk4D0fN6qVavMRd/6isWDDz4YyDNwmjmsz7llVrD0fdwuXboEa1SKe9NVVM221uxhfUP5gAMOMNnCmkjy1ltvyfjx42XSpEnmRZ/CCtfMxOPjUPE37u+p6CpgpmRW/9q7EIHE3Xj4NO6zQADG3UMJnx+BKDkO1MSDo48+WqZNm2aSEDJFz+TplrC+A6zi8OyzzzarSkGVNm3amOzWzMqfbi+7uTomqHmkuZ/NmzfLjBkzzKs8KuD1lyaZ6MrrrbfeKprYpauFPEuXm69A/29aICLVRKSw/7NUBC4WkQMcbgcTd3Pjx6SNigBMmscSNl8CUTIcllkFUiGWf0t348aN5mJmvepFxYPTBBGnVqso0bOBhx9+uDRr1kwefvhh+f333502p14ABDR55NlnnzVZ/iq+EYEBQHXZRSMRmeKgTWMRmeqgHnHXASSqCAKQjyBUAgSiUPEG0nn+++R0W/eqq/TE0V9FX/W48cYbTZZpu3btAhkvfyd6NvCCCy4wv1WtWjWTFZztbsnAJ0GHhoAKbz2vqc858j5xtB9FRxEZ7WDITiIyxkE94q4DSFRBAPINhEuAQBQuX7+968qfXjysmbl6+bJuw+qFzqeccorZEvzXv/4lixcvNqt/QRZNStDrY+677z7RF0h0hVHLEUccYbYjKbkjMHHiRDn99NPN90DySDR+YAUwGs6Msj0BVgD5IkIlgAAMFa/vzlXk6buzeg7v66+/Nmf/dBWuX79+JtHjzDPPNPfM6VUvQZRly5aZBJPnn3/ePBunq4+atKDj6vbj8OHDjfig5JaAnvts2rSpydbmepnwfcEZwPAZM8KOBBCAfBWhEkAAhorXd+eabKGre5ltYP2nrv5p1qgW3ZbVFbmKFSv6GkuTEPRFEL1nUFeWGjZsaETnf/7zH9lrr7189U3jcAjoyq+KQH2LWAtnA8PhnOmVLOBw+dI7ApBvIGICCMCIgXsYTt+W7datW15LTQpo3LixWZnTP/Ob8atXx+glz998840ceuihJuP3oIMO8jBTmkRNQLOE9Rt44okn5Msvv8wbnqtjwvFEYfcA6uGIHi6ugNGZEXfD8U/aemUFMG0ejZk9BKKYOaTAdHQLeN9995UffvjBPMkWRtHLonUr+bHHHpOLLroojCHoMwIC+o3os34qCvXspt4xyKpg8OB5CSR4pvRYOAEEIF9GqAQQgKHiDaRzXeHTH+xhbMWqQDjjjDPkzTffFL1Sxu9qYiAG04lvApqprfdCUuJJgLgbT7/EbVYIwLh5JGXzIRDF36H16tWTCRMmhCIAZ82aZV6o4Hm3+H8Hbmeo2/maOESJHwHibvx8EscZIQDj6JUUzYlAFH9n6nmu3XbbzbzXG/Rdf/PmzTPn/a699lq555574g+DGTomcM0115iLu8kSdowssorE3chQJ3ogBGCi3Rf/yROI4u8jnWHmGbDbb7/dXNOi/+23DB061LwuoWfG+vfvLzfddJPfLmkfIwK6rd+iRYsYzYipZAgQd/kWnBBAADqhRB3PBAhEntHlrGH37t1l8ODBvsevVauWLF26VMaPHy/HHXec7/7oIF4E1q9fb54GZAUwXn7R2RB34+eTOM4IARhHr6RoTgSi5DlTV/8041OvgfFa9G3Zc845R5o0aSJvv/22125oF3MCnAOMp4OIu/H0S9xmhQCMm0dSNh8CUfIcqmcC9T1gverDa9Hzftdff71Mnz7dPCdHSR8BfTu4fPnyrADG0LXE3Rg6JYZTQgDG0ClpmhKBKJne1OtaXnvtNc9nvFavXi1Vq1Y1GcAffPBBMiEw62IJTJkyxazwUuJHgLgbP5/EcUYIwDh6JUVzIhAl05m6DbzzzjvL999/b4Scl6KZvw899JC5YubEE0/00gVtYkzgtttuE/21ZcuWGM/SzqkRd+30u1urEYBuiVHfFQECkStcsaqsW8E9e/b0fH3LypUr5Z///KdJBPnoo4+kbt26sbKPyfgjoKt/U6dO5TUQfxhDaU3cDQVr6jpFAKbOpfEyiEAUL3+4nY2uAi5ZssTzJdELFiyQI444wpwVW7RokdvhqR9TAnoBtG7vb968OaYztHtaxF27/e/UegSgU1LU80SAQOQJW2wa6VnAvn37yi233OJpTosXL5b99tvPvDPcuHFjGTNmjLk6JFP0CpELL7xQjj/+eLngggs8jUGjaAn89ttvcvTRR4uKe7Z/o2XvdDTirlNSdtdDANrt/9CtJxCFjjj0AXbffXdRIaevhbgtug2sz8DpiyBffvmlHH744fLpp5+abvSfTZs2lV9++UUOOeQQmTt3rqPu33//fXnxxRdlzpw5snDhQvOOsYrHO+64w1F7KnknsGHDBmnVqpXZ+kX8eecYdkvibtiE09E/AjAdfoytFQSi2LrG8cQ0IUSvhNHzgH7KAw88IFdffbXssssuoiuLGzdulD333NNsI1avXl2++OKLYrvX1cIzzzxTXn31Vdm6davsuuuuZnVRhaOuSM2YMcPP9GibhYCKP33547333jP8KfElQNyNr2/iNDMEYJy8kcK5EIjS4dS9997brLbpmUCvRUXDkCFDZP78+bJp0yYj4PSuQBVvmmmsdwYWVlQo6ruzI0aMEH19QoWoPiu31157yR9//GEEZZcuXcyfU8IhoNybN28u06ZNQ/yFgzjQXom7geJMbWcIwNS6Nh6GEYji4Qe/s9BVwHvvvdes4AVdKlSoYBIK3nnnHdFVpm3bthlxqEV/r2XLlkYwnnzyyeYsYsOGDfOmoOfRdBWxU6dOMmrUqKCnRn8ism7dOjn11FPNCivbvsn4JIi7yfBTrmeJAMy1B1I+PoEoPQ7Ws4C6Cqj/DLLo2UJdHdTkkOXLlxuRUbZsWSldurQRftWqVZOXX37ZiMTCimYZ6wqV3llICZaA+kWve9GVP8RfsGzD7I24Gybd9PSNAEyPL2NpCYEolm7xNCk9t3fjjTfK7bff7ql9UY2GDx8uo0ePFl3N23fffeWFF14wF0fXqVPHPGp/7rnnysEHH1zkmDVq1JDKlSvLzJkzA50XnYmMHDnSvOlMSRYB4m6y/JWr2SIAc0XeknEJROlytJ4B1IzeKlWqhGKYnulTsaeJIXqBtJNSqVIl8x7trFmzRMUgJRgCuqp60EEHmVVZ3ZanJIcAcTc5vsrlTBGAuaRvwdgEonQ5WV8Hueiii+S///1vaIbpamDnzp1Nwshll12WdZy3335b2rZta7YoX3nlFTnllFOytqFCdgL6zNutt95K0kd2VLGrQdyNnUtiOSEEYCzdkp5JEYjS48uMJSoCP/vsMznssMNCMU7Pnel5vx9//FF+/fVXc2VMtvLtt9+aK0q0zdlnn23+XS+erlixYram/HkhBPTMn57946WPZH4exN1k+i3qWSMAoyZu2XgEovQ5XJMzNPFCs0JVDIZRnn76aenatatZ0WvdurWjIVatWiUdO3aUt956K6++JpPoSyZ6dpHijMAPP/wg9evXlzVr1rD65wxZ7GoRd2PnklhOCAEYS7ekZ1IEovT4sqAl999/v+/LoYui06xZM3n33XfN+TO3L5DoyyJ6OfR3330nKiQnT54sxxxzjElm0O3rMmXKpNcpPi179NFH5e677zbvNuu5SkoyCRB3k+m3qGeNAIyauGXjEYjS63BdXfv666/lwAMPDNxITT7Qa2E++eQTX31rUoleNq1CUFcIVUzquTa9ykZXMevVq+er/zQ1/uabb0zmtW6589JHsj1L3E22/6KaPQIwKtKWjkMgSq/jVSh0795dHn744cCNzH85dBCdq6DRFcXzzz/fZDFr0cuta9asKePGjZO6desGMUyi+7jhhhvMk3/c95doN5rJE3eT78MoLEAARkHZ4jEIROl2/uDBg40IDLro825NmzY1F0AHWfRZuTlz5phVrokTJ5rzgXr3oK5+2V5q165t2FCST4C4m3wfRmEBAjAKyhaPQSBKr/M1AeTf//63WVkLuug5vTPOOEPGjh0bdNfb9demTRvzyoWeNbS96Cro7NmzbceQCvuJu6lwY+hGIABDR2z3AASi9Pv/9ddfN+/1BlX0RY9//etfcumll4Z636DOV1cA9WUTfXNYr42xuShzXlNJxxdA3E2HH8O2AgEYNmHL+ycQpfsD0FVAPUunW8EXX3yxb2N///13k1SyYcMGWbJkiesMYLcT0POAmnCiz809+eSTbpunqv5JJ50kkyZNSpVNthpD3LXV8+7sRgC640VtlwQIRC6BJbi63tv3wAMPyJ577unZitWrV8vee+8tDRs2NKtyYRe9cFqfkOvZs6fotTY2l9NPP11effVVmxGkxnbibmpcGaohCMBQ8dI5gcieb0BXA8uXL2+ElIpBXRl0W1asWGHeAt5vv/3kiy++cNvcdX3NCtaVP818dfLiiOsBEtRAL9F+/vnnuQImQT4raqrE3RQ4MQITEIARQLZ5CAKRXd5X0bdt2zY57rjj5LHHHnP9XJwmZDRo0EBuvvlm0bdowy6aCaxvB48cOVK6dOkS9nCx7l+TbnQFUP1HSTYB4m6y/RfV7BGAUZG2dBwCkZ2O1+fi9O493Vq95ZZbzMpg1rJunWzr0kXmjh8vc//8U1rqO8BO2mXtuOgK+tqFvhesokfPHu60007ml2Yh6y+94sbpU3Q+phGLprrqqm8pU5JPgLibfB9GYQECMArKFo9BILLY+SJmW1W/Ab1k+Iorrig6qaN+fZGPP94Olq5DldDfnzkzVIh6bnHYsGGir4boLxWF+k/djtbXQj7//PNQx49D52vXrnUm0uMwWeaQlQBxNysiKmh8FRHW+/kUQiNAIAoNbWw6LikiDUWkiogsFZH3RGRrgdmpENRVwN69e8tll10mu+666//VOPhgkR9+2MEeIwD1dyMQgQUH//nnn6Vy5cpy+eWXy6BBg2LDOqyJfPjhh+ZOR0o6CBB30+HHsK1AAIZN2PL+CUTp/gDaiMhDIlI9n5mLROQqEXmpENP1jKBmCd90003mnr9y48aJnHNOdkhr14rstlv2egHVWLZsmVSpUkXOOeccGTFiREC9xrcbPa+p/uD8X3x95GZmxF03tOytiwC01/eRWE4gigRzTgZR8Tfu71W6/Pm+unKnv9oXIQJ1sioEhzVvLue+8YazuetZwDVrnNUNqNaJJ54okydPNk/FffLJJ+Z6mrQW3Z5XEbh58+a0mmiVXcRdq9zt2VgEoGd0NHRCgEDkhFLy6ui2788iUiGzTVvABBWAK0WkciHbwVpVxeMLRbQtkkbEq4B6FnDUqFGiV8WceeaZMmbMmOQ5yuGMGzVqFMqTfg6Hp1rABIi7AQNNaXcIwJQ6Ni5mEYji4olg59FERJy8GXGiiEwuMLSKR801repWAB5zzA6JIsFaVXhvp512mhFHv/32WxTDRT7Gc889J3oHICU9BIi76fFlmJYgAMOkS98mA1R/cGoCgGYaUtJB4FYRucWBKXqTX98C9RqJyBQHbQut8sILIm3bem3tqZ1eZfPoo4+KPlOXxlKrVi359ttv02iatTYRd611vSvDEYCucKWnst75pZft6jknzXb86aefzHbXgAEDtjsHdNhhh8mQIUOkfv365loMzYi85557HIMgEDlGlaiKfgSgrjWN9mpt9eoi8+eLlCrltQfX7erVqye//PKLzNdxU1Y06WPnnXc2195Q0kOAuJseX4ZpCQIwTLox7vvUU0+Vs846S0aPHi3ff/+91K1bVx5//HHzIsJ1111nZq5BZO7cuaKvJQwcONC86qDPZvXo0cPUdVIIRE4oJa+Ony1gXyuAimryZJHGjQOFpuf99Kqags/B6fdfp04dcwZQ/19JW1m6dKlUraqb8ZQ0ESDupsmb4dmCAAyPbeJ6vvbaa81VEAcddJCZ+yWXXGJWBHWFMJMdqEJQn4yqXbu2I/sIRI4wJa6SnyQQbauHAXbxavWzz4p06uSptR5H0Pdu9ck5fWtYX77Q39u0aZPpT7OT9QWQzp07m23f9u3by1tvvSWLFy82L4akrWSe3kubXbbbQ9y1/QtwZj8C0BknK2r1799fmjVrZrZ7tej9Z3p2r00bzdn8qzRu3NhcjaF3ua1evTorFwJRVkSJrZDJ5DXCKZ8VmZvl2xVzDYyeC+zn1XKPK4CZZ9/WrFljhJ5e76IvfRx66KHmGTj9Pa0zadIkGTdunLmsWrdG99lnH1m0SG83TF+58sorZfDgwdz/lzLXEndT5tCQzEEAhgQ2ad3qqp/edaargE888YSZ/oQJE8y5J10JzBRd+fv666/NCuCcOXN2MFNXT8qWLZv3+xqIlixZQhJI0j4Ih/Mt7CJozfDtUYz4065PFpG3HY6xXTUfZwD1rjsVOyrumjdvLrvsUvQa5IwZM8z1L7oFrN9zGhOYxo4da46BUNJHAAGYPp+GYRECMAyqOexTt2h79epV7AwKZv3pGaCpU6fKlClT5MILL8xr60UA9u3bV/r123FthyzgHH4UIQ/t5Cm4glPQtzUcvP+x48xbtxZ55RXXFuk2r17kfMopp8hrr73mqL1uFXfo0MGcedX3gtNU9C9xmtyyceNGVv/S5Ni/bUEAptCpIZiEAAwBai671HNKFSro9bxFl3nz5uWd6dPnrlT4TZ8+Xbp167bdDwMvW8CsAObS+8kZWy+B9nyZiyYp3X23K2M1e123O2fOnJl3xCFbB+vWrZP999/f/P+UpmtSdDXz6KOPlgULFpgtb0r6CCAA0+fTMCxCAIZBNSF96sqfnufTrd8uXbrI1q1bt5t5JglEz0BlflBoUkjbtm1JAkmIj+M6Td0i9rymVqKEyMaNImXKODZPrzy6/fbbRbNeNanJadG/FL388suOzrs67TPX9ZSDrtJv2bIl11Nh/JAIIABDApuybhGAKXOoU3NU/OnK38KFC6Vr167b/TD4+Wd95EvMuT1d+dAsyLvuustcFaPXwOjFuFwD45Q09QojUFpENO9WA1D+BBLHtJo0EZnk5C2Sv3r897//bc6urly5UkqX1tGdFb3+Rf8/0Tsw01I0s/mll17a4S98abEPO/66wosL+PkSshFAAGYjlNI/V9E3fPjwQq3TbMhMyX8RtP7w1K20u11svxGIUvoB+TRLJdhGEdHzg14EoGYaf/rBB1LnyCOlXLlyxc5GV7Y1o/c///lPXoKT0+nrfZlfffWVSWRKSznkkEPM3Z+U9BIg7qbXt0FahgAMkiZ97UCAQMRHURiBq0TkQZ9ortE+SpY091bq1UX6nm2rVq3kmWeeMX9R0SxfFX6ara6C58033zTXHLkpxx57rGn/8MMPi74J/Pnnn8vQoUNNZrxeiZSkoncdPvbYY2YFv+BxjyTZwVyzEyDuZmdEjb/+8p25tgseEAicAIEocKSp6PB/ItLCpyVrdt9dup92msyaNctcV6Rv9b7yyivmSINmt2ryht7jlznSoMcZatas6WrUiy66yKwa6pNpWnR1PPPvlSpVMmdo9R7BuBdN8tLtbF3JzMw/7nNmft4JEHe9s7OpJQLQJm/nwFYCUQ6gx3xI3fZdKSJ7+p1nyZIi+oZtqVJG8OnFzr/++qvpVS9zzqzQ6Xm3Tp06mXNRep71qKOOcjXy+vXrZfbs2eblEBWA+hLO66+/bq6HUZGZhO3U448/3mT6s/LnyvWJrUzcTazrIp04AjBS3PYNRiCyz+fZLPb9FnD+AfK9CvLDDz+YRCV91UMFW/6i/33ccceZbHbN6m3atGm2aWb98xtuuMFsDevKYxyKruzpKue7775rfumqqBb9/Q8//DAOU2QOEREg7kYEOuHDIAAT7sC4T59AFHcPRT+/i0Xk0aCGLfAusK4EapZvYZm+y5cvF01q+uWXX2TYsGFmq9hP0VdF9G5BFZUldTUyR2XixInm3WLdjtYVUF2lLFWqFHf85cgfcRiWuBsHL8R/DgjA+Pso0TMkECXafaFMfoOIFJ+362LYiRNFTjrJcQMViEceeaRZKdM7LW+88UYjmLyUF198Udq1a2eyhHNxDlDP811//fXy7LPPGsHHvX5evJjONsTddPo1aKsQgEETpb/tCBCI+CAKEtgsIs5v4svCz6UA1N70HNxJJ51k7vfTpJB//etf5td5552X9UqZ/LNZtGiR1KhRwzwxp3euab+aHfzQQw+F4nTtX59nzPzS7GSEXyioE98pcTfxLozEAARgJJjtHYRAZK/vi7J8rYjsFhSWAlvAbrrV10HGjRsnP/30k+izb3o5+n//+19zlYzToitwKsj0fW3NNtazdnrdStAlc2G7vtmt29s84RY04XT1R9xNlz/DsgYBGBZZ+jUECER8CAUJnCEiLwWFJV8SiJ8uVcTpc4h62XnLli3NPYIHHHCAqy71ibX+/fsHKgA1gUPPK1511VXmShuEnyuXWFuZuGut610ZjgB0hYvKbgkQiNwSS3/9jiIy2oGZmQtKizyhV726iGa6lirloLfsVXSLNZPZu3nzZtFXQM4//3yzIli2bNmsHZx99tnmHkLdDg6iqNjr0KGDebaNAgE3BIi7bmjZWxcBaK/vI7GcQBQJ5kQNcoqIvOVgxouOPVb2nT698KfiNHFj3DiRtm0d9OSuyrJly+Taa6+V1157zYi5KlWqmHeE99hjjyI70pU6vYdQzwN+9tln7gYspLb21717d5Pdy8XNvnFa1wFx1zqXezIYAegJG42cEiAQOSVlTz29ke5YB+b+evDBctH338voPfeUnVat+r8WuvL34IOhiL+C09Jn5XRreMSIEXLOOecUOetPPvlE6tWrZ7aNdcu2SZMmDiwsusqDDz5onmyjQMALAeKuF2r2tUEA2ufzSC32E4g0U7S7iBwkIj+IyBAR+TPS2TNY0AT0tjy9NrmMg451C7idiFz33nty3J9/iixdKlKlikjDhoFt+2abhp4F1Lv+Ro4caYRgUeWWW26RgQMHmkzgcuXKmWQQvXPQS9FXTE4++WRW/rzAo40h4CfugtAeAghAe3ydE0u9BqI7ReSaAteFqPi7T0R65cQSBg2CgJtXQFQALhKR7996S048RTeOoy16qXK1atXMvYHTpk0r9rJnTRzRp9a+/PJLc7XMPvvsI/oyiZeiV8noKiIJH17o0QYByDfglAAC0Ckp6nki4EUAqvi7/u/R8icAZJIC7kYEevJFHBo5TQDJP9f1r78uu7ZsGfn09em4jz/+2Ig6vealqKJn9CpXrmyeoNOVvzZt2sj48eM9PxF39NFHB3KOMHJgDBgbAl7ibmwmz0QiI4AAjAy1nQO5DUS67atbhJrXWVj2p4rALX+/JMF2cPK+KX2zY6LLaX/Vu7fUHTDAZSt/1fXsn2b16tUuvXv3LrIzvfNPM4W1vm4B9+rVS0455RSZMWOG52xgvVxaL5mmQMArAbdx1+s4tEs2AQRgsv0X+9m7DURXiciDDqzqISLhvLfgYHCqeCagfnvAZesXr7hC2j78sMtW/qrvtddesmbNGrn//vulRYsWcvDBB+/wZNyKFSukdevWMnPmTJOxe+edd5orY/RN3kaNGpl/eikXX3yxPPnkk2wBe4FHG0PAbdwFm50EEIB2+j0yq90GIv0xf4WD2Q0SkSsd1KNKvAi8LiJuNnM3lyghO23eHFnSR4bWhRdeaF4JURGoW7z6HWuG7z/+8Q+pUKGCfPrppzJr1iyT9KGvh2jp0aOH2fbVVUB9ZWTnnXf2BP+mm26Su+++GwHoiR6NEIB8A04JIACdkqKeJwJuBSArgJ4wJ6KRZgCv+3v73tWEN2wQKVfOVZOgKuszcUOHDjUJHgsWLBBd9dNt34oVK5osX135U+H37bffyhFHHCFPPfWUHHXUUb6G14SSN998kyxgXxTtbuw27tpNy17rEYD2+j4Sy90GIs4ARuKWnAxyooi842Xk7t1FBg/20jLUNroKeMIJJ5hVwPvuu08uuugiKRXAqySaUKLvClMg4JWA27jrdRzaJZsAAjDZ/ov97L0EIrKAY+9WTxP8QkQ83YzXtKnIhAmexgyr0dixY02SiJ4V1IxfXf0Lorz//vvSUO85pEDABwEvcdfHcDRNKAEEYEIdl5Rpew1E3AOYFA87yTN5dAAAIABJREFUm6eu7G4SEd0Gdl1itgKo2b56xq9OnTpG/FWtWtW1SYU10C1mTTiZO3eubNmiue4UCHgj4DXuehuNVkklgABMqucSMm8/gYiXQBLiZAfT9JL9m9dtDs8AFjRtwIABoq9+NG3aVJ5//nnZbbfdHFifvcrEiROlffv2sn79epI/suOiRhYCfuIucO0hgAC0x9c5sZRAlBPssRv0BRFp62VWJUuKxGQ17J577jEZvs2aNZOXXnpJypRx8qBddqN15a927dryxx9/mPOEFAj4JUDc9UvQjvYIQDv8nDMrCUQ5Qx+rgUeJyH+8zmjtWpGAVtq8TuGhhx6Sq6++2rzR++qrr0rZsmW9drVdO71ipnnz5vLOO++w8hcIUTpRAsRdvgMnBBCATihRxzMBApFndKlq+L6IHO/VotNPF3n5Za+tfbfTbd8+ffqYflT86WXPQZXnnntOOnbUB/IoEAiOAHE3OJZp7gkBmGbvxsA2AlEMnBCDKawVEc+n5Q44QGTevJxYoVuyeqHz5s2bpWTJkmaLVsVgv379fF/5smrVKvPCiP5TVwIpEAiKAHE3KJLp7gcBmG7/5tw6AlHOXRCLCei7zfq+s6dSr57IRx95auq30ejRo6Vz584ycuRIs1J32mmnyVtvvSWNGzeWW2+91WQAH3TQQZ6G0XsD9ck3Mn494aNRMQSIu3weTgggAJ1Qoo5nAgQiz+hS1XCjiHg+NffGGyLNm+eEh77pqxc+//rrr7LTTjuZOeilzzfeeGPeqqC+BKJi0E1GsGYQd+jQISc2MWj6CRB30+/jICxEAAZBkT6KJEAg4uNQAr+KyJ5eUJQuLbJxY+RvAWemuvvuu8upp54qevFz/rJw4UL56KOPZNiwYWZFcJ999pFHH31UWrdundXKRx55xDwhp4Wt36y4qOCBAHHXAzQLmyAALXR6lCYTiKKkHd+xPAvA558Xad8+J4YtXrxYatSoYYSdbtcWVd59910588wzZfny5dK1a1cZNGiQycIsWPT8YN++feX222/PiT0Mag8B4q49vvZjKQLQDz3aZiVAIMqKKHEV9DUPfaysiogsFZH3RCTb7XWetoD1nr1N+n5Ibope+Ny/f3/57rvvTLJGcUXFnYrE4cOHS7Vq1Uy2cP7n4fSuvy5dusi0adNyYwyjWkWAuGuVuz0biwD0jI6GTggQiJxQil+dokReGxF5SESq55vyIhG5SkReKsYMz0kgOcyObdeunbz++uuyceNGKVFCQ2X2MmXKFJMosssuu5izg/vuu6/Mnz9fDjvsMNm0aRN3/WVHSI0ACBB3A4BoQRcIQAucnEsTCUS5pO9t7KJE3rMict3fXeZ/0zez+qcbtYWJQH0vQ1cAnUmoAnNevlykUiVvhvhspdu6r7zyinmhw0354osvpH79+nLIIYfIzJkzzQXSTzzxBNm+biBS1xcB4q4vfNY0RgBa4+rcGEogyg13r6Oq+BtXhMjLCLjChJyKwMUickAh28GDRORyrxOqUEFk5UqvrX2102tfXnjhBZPt66bodvABBxwgeoZQX/g45ZRTWPlzA5C6vgkQd30jtKIDBKAVbs6dkQSi3LF3O7Ku6i0QkWoikn+Fz00/jUVkaoEG40XkVDedFKz755+RZwGvW7dOKlWqJNWrV5e5c+e6mv0ll1wijz32mAwZMkT+97//yYQJE1j9c0WQyn4JEHf9ErSjPQLQDj/nzEoCUc7Qux64kYhMcd1q+wadRGRMgT6Wi4ivTdzJk0Uaq7SMrrzxxhvSsmVLc8WLruC5KZnLo7VN5vUQN+2pCwG/BIi7fgna0R4BaIefc2YlgShn6F0PrC/SjnbdavsGBVcAfZ3/y3T97LMinVRaRlcGDhwovXv3Nnf91dOXSFyUP//8UypUqCC//fabi1ZUhUBwBIi7wbFMc08IwDR7Nwa2EYhi4ASHU/CzAljUGUBf5/8y887BCqC+9tGrVy+pWLGivPbaa2Y7uHTp0uaKl2xFr47R+/645DkbKf48LALE3bDIpqtfBGC6/Bk7awhEsXNJkRPKdgZQRV4mYDjNAv5cRI7wg0DvAtywIfIzgPfff79cc801UqpUqe3O7+mrIPoM3AknnCArVqyQcuXKmUufVexNnTrVCD+9GJoCgVwSIO7mkn5yxkYAJsdXiZwpgShZbisuC1gtuUdEOhe4B/BHEelRyBUwKhL1ApVSfhD06SPSv7+fHjy11SSQnj17mutb9A7AzGqergLqFm/58uXNFq8KxH/+85/mz6dPn76DYPQ0OI0g4JMAcdcnQEuaIwAtcXSuzCQQ5Yq893ELuwcwv8hz+hLILSJyq/dp/NUyBxnA+ac8fvx46datm6xcubLITF4ViJrssWXLFr/W0h4CgRAg7gaCMfWdIABT7+LcGkggyi1/r6M7FXlF9a/t9RG30l4noO2uv17krrv89BBI29WrV0uDBg3km2++Eb3jjwKBuBMg7sbdQ/GYHwIwHn5I7SwIRKl1bbGG6ashvvN2W7cWeeWVWACsUaOGLFqkj95RIBB/AsTd+PsoDjNEAMbBCymeA4Eoxc4twjRd9dPVP6+XSed1e9hhIl98EQuAmg28Zs0aXvSIhTeYRDYCxN1shPhzJYAA5DsIlQCBKFS8sey8l4gMDGJmZ5wh8lJhrwsH0bm7PhYuXChdunSR999/311DakMgBwSIuzmAnsAhEYAJdFqSpkwgSpK3gpmrrv7pBdC+y9q1Irvt5ruboDr4+uuv5dBDDw2qO/qBQGgEiLuhoU1VxwjAVLkzfsYQiOLnkzBnFMjLHzrB+vVFZs4Mc6qu++7atauMGjWKRBDX5GgQNQHibtTEkzkeAjCZfkvMrAlEiXFVIBP9r4hc6ren8uVF1qzx20ug7YcPHy7nnntuoH3SGQTCIkDcDYtsuvpFAKbLn7GzhkAUO5eEOqGlIlLZ7wjdu4sMHuy3l8Daz549W4455hjZtEk3tykQiD8B4m78fRSHGSIA4+CFFM+BQJRi5xZi2ioR2cOvyTE5+6eve3z22WfSqVMn+eGHH7jo2a9faR8ZAeJuZKgTPRACMNHui//kCUTx91GQM/xYRI7x2+HEiSInneS3F8/tf/zxR3nmmWdEt33nzp3L826eSdIwVwSIu7kin6xxEYDJ8lcosy1TpozMmDFDjjzySPNr1qxZeeMcdthhMmTIEKlfv76sWLFCBg0aJPfcoy/COisEImec0lLrZRE53a8xIbz/+9NPP5kEjj333FP22WcfqVSpkugLH0uWLBH9M/2n/po/f75otq8+7carH34dSftcESDu5op8ssZFACbLX6HM9sEHH5RDDjlEWrRosZ0A1CCiKyATJ06UgQMHiorBJ598Unr06CGPP/64o7kQiBxhSk2lN0WkmV9rQhCAo0ePls6dOxc6s9Kl/3qw7k99d5gCgRQQIO6mwIkRmIAAjABynIdo1qyZ3H///dKuXTuz8pF/BfCSSy6RAQMGSOXKlWXz5s3GDBWCZ5xxhtSuXduRWQQiR5hSUymQLOAQtoBXrVole++9NyIvNV8ahhRHgLjL9+GEAALQCaWU1tEfiJ988okRdCtXrpQFCxZsJwBHjBgh5cuXlzZt2uQRaNy4sUyePNlspekWWrZCIMpGKF1/PkRELvNj0s47i6xbJ1KqlJ9edmirf4HRb3b9+vWB9ktnEIgjAeJuHL0SvzkhAOPnk8hm9MYbb8i0adPMKt9+++23gwCcMGGCOROlK4GZoit/ulKo/5wzZ84Oc9XzhGXLls37fQ1EerZKheRaze6kpJrApyJylB8Lu3UTeeopPz0U2la/8wYNGgTeLx1CII4EEIBx9Er85oQAjJ9PfM1It2h79dLXWIsutWrVkqZNm0qHDh2kUaNG5rB7UAKwb9++0q9fvx0GRwD6cmsiGpcUEZX4u/iZrd61VyaQh+S2mwUXOftxCm2TRgABmDSP5Wa+CMDccA9t1IoVK0qFChWK7X/evHkyduxYadWqlehdZ5mih+H1ILxegdGtWzfxsgXMCmBoro19x41EZIqfWV5zjci993ruYfr06bLHHnuI/gWnYNEs92OPPdZz3zSEQJIIIACT5K3czRUBmDv2OR25evXqZls2U6pWrSpvvfWWSQbRH5a6bZtJAtFrMzIZkrpd3LZtW5JAcuq9eA7eUURG+5na5MkijRt76uGCCy6QYcOGmbZ6ZdH5559vvm/9y863334rY8aMMX/ZIdPXE14aJYwAAjBhDsvRdBGAOQIft2EL2wLWH6D6w1OF4V133SV169Y118D07NmTa2Di5sAYzOcWEbnVzzyefVakUyfXPWzYsEF23XXXvHb57/DjihfXOGmQAgIIwBQ4MQITEIARQE7CEIUJQJ13/ougNVNYL4K+++67HZtEIHKMKtEV9fzfryKyux8rfKwA7rXXXqJXvVAgAAER4i5fgRMCCEAnlKjjmQCByDO6RDVsIiKT/MxYjyP8+qvn61/++c9/ykcffeRnBrSFQGoIEHdT48pQDUEAhoqXzglEdnwDuvWrW8CeS9++IoVkjzvt74477pCbbrpJdMtXz/pt2bKFN3ydwqNe6ggQd1Pn0lAMQgCGgpVOMwQIRHZ8C74FoM8MYKX82WefyQMPPCDvvvuu2QL76quv7ICPlRAoQIC4yyfhhAAC0Akl6ngmQCDyjC5RDX1vAevLHxs2+LoDUF+yOfPMM+Xjjz9m9S9RXw+TDZoAcTdoounsDwGYTr/GxioCUWxcEepENAlkk4iU9jPKAw+I9OjhqYfx48dLx44dzVNvXPXiCSGNUkSAuJsiZ4ZoCgIwRLh0TTaaTd/AOBFp58fgyy8XGTTIdQ+64qdJICVKlDCv2lAgYDsBBKDtX4Az+xGAzjhRyyMBApFHcAlspmt3D/iZt8cVwMGDB8uVV1653as2fqZBWwgknQBxN+kejGb+CMBoOFs7CoHIHtd3FpFn/JirZwDLlXPdg4q/Rx99VDZv3uy6LQ0gkEYCxN00ejV4mxCAwTOlx3wECET2fA6+3wL2eBF0kyZNZMoUX68Q2+MkLLWCAHHXCjf7NhIB6BshHRRHgEBkz/ehiSALRGRfEdHA4ro4fArujz/+kFmzZplkj19//VVatWrF9q9r2DRIMwHibpq9G5xtCMDgWNJTIQQIRHZ9Fn1EpL9XkydOFDnppGJb62sf55xzjsyZMyevXqlSpczFzxQIQOAvAsRdvgQnBBCATihRxzMBApFndIlseJ+IXO115hMmiDRtWmTrO++8U3r37i0lS5ZE8HllTDsrCBB3rXCzbyMRgL4R0kFxBAhE9nwfbUREr4LRrWBPpU8fkf6Frx8uW7ZMqlWrxjUvnsDSyDYCxF3bPO7NXgSgN260ckiAQOQQVMKrZc7/VfMjAHv3FhkwoFAS9957r9xwww0IwIR/J0w/GgLE3Wg4J30UBGDSPRjz+ROIYu6ggKZ3YsmS8o7fS5jvu0/k6sI3kDt37izPPfccAjAgf9FNugkQd9Pt36CsQwAGRZJ+CiVAILLjw+goIqP9mjpqlMh//lNoL6tXr5Zzzz1XXn75Zb+j0B4CqSdA3E29iwMxEAEYCEY6KYoAgciOb8P3HYCKKcs9gNu2bZMhQ4ZIz549zbUv+oun3+z4vrDSHQHirjtettZGANrq+YjsJhBFBDrHw+gZwKUisrfXeey7r8iCBSKlSmXtYfbs2fLiiy/K9OnT5Y033shanwoQsI0Acdc2j3uzFwHojRutHBIgEDkElYJqugWsW8GeygsviLRt66rptGnTpEGDBq7aUBkCNhAg7trgZf82IgD9M6SHYggQiOz4PNqVKCFjt21zfwXMzjuLPPOMa/G3atUqcyH0+PHjzYsgFAhA4P8IEHf5GpwQQAA6oUQdzwQIRJ7RJaahryfgTjtN5LXXHNv6888/y/333y+DBg2STZs2cQbQMTkq2kSAuGuTt73bigD0zo6WDggQiBxASngVXwkgBxwgMm+eIwJ67q9Tp07mFRCefnOEjEqWEiDuWup4l2YjAF0Co7o7AgQid7ySWNvXFTBHHy3yySdZzd6wYYMceOCBsnz5cpP9S4EABIomQNzl63BCAAHohBJ1PBMgEHlGl5iGvlYAdftXt4GzlAEDBsjNN9+M+MsGij+HgIgQd/kMnBBAADqhRB3PBAhEntElpqGeAfyxRAnRZ+DEzepc6dIiGzdmvfpF3wHW1b/ff/89MUyYKARySYC4m0v6yRkbAZgcXyVypgSiRLrN9aT1ApdxJUqIBhTHIrBePZGPPipyrIULF5qLnx977DFZv3495/5ce4UGthIg7trqeXd2IwDd8aK2SwIEIpfAElx91bBhssf114v88otzK665RuTee3eoP2rUKOnatauUKFEC4eecJjUhYAgQd/kQnBBAADqhRB3PBAhEntElqmHJkiXlnP/8R56aNElkyRJ3c9+wQaRcubw2esFzo0aNEH7uKFIbAnkEiLt8DE4IIACdUKKOZwIEIs/oEtfQczJI2bIil16alwzy8tCh8tDYsfKuiGxNHAUmDIHcEyDu5t4HSZgBAjAJXkrwHAlECXaey6l3LlFCnnGTBJKl/+UicpmIvOByHl6qly5d2qw4Zq6Y0a1n/ffMP730SRsI5IoAcTdX5JM1LgIwWf5K3GwJRIlzmecJn1iypLyzNdg1O73x724R6eV5VsU3LFWqlOy0007SoUMHmTVrlvmlYvD222+XSpUqybXXXitr1qzhxZGQ+NNtOASIu+FwTVuvCMC0eTRm9hCIYuaQEKdTSkR+32cf2ennnwMbJXPl81ki8nxgvYoRefqG8BVXXCFVqlSRhx9+WPbZZx/zzJw+Nzd//nzRs4gTJ0409bYGLGwDNIWuILADAeIuH4UTAghAJ5So45kAgcgzusQ1VFF17YEHyh1z5/51HUyAZYuIqAgMejt41113NVfM6Favrgaq2NOi/65bwAi/AJ1IV5ERIO5GhjrRAyEAE+2++E+eQBR/HwU5QxWBrf7806zW6YpgkEVXA9uJyEtBdkpfEEghAeJuCp0agkkIwBCg0uX/ESAQ2fk17FSypCzeulX2DtB8PV24WEQOIDs4QKp0lUYCxN00ejV4mxCAwTOlx3wECET2fg76RNwGESkbMILGIjI14D7pDgJpIkDcTZM3w7MFARgeW3rmRnprvwE9Q6fXqty1yy5yvV70HGA5W0RGBdgfXUEgbQQQgGnzaDj2IADD4UqvfxMgENn7KWhiRYlt22Rd6dKy859/BpYYomcB3xcRXQkM9tIZe32F5ekiQNxNlz/DsgYBGBZZ+jUECER8CG1FZNzfGILMDtZ1xS4khfCBQWAHAsRdPgonBBCATihRxzMBApFndKlq2EZEhopIxQCt0pVA/dUeERggVbpKAwHibhq8GL4NCMDwGVs9AoHIavdvZ7wmhfQRkX4BrgaqANTM4P3ZDuZDg0AeAeIuH4MTAghAJ5So45kAgcgzutQ21NXAh0SkeoAWnlSypEzitY4AidJVkgkQd5PsvejmjgCMjrWVIxGIrHR7VqN1NbChiFwtIi0DuDS6k4iMyToqFSBgBwHirh1+9mslAtAvQdoXS4BAxAeSjUBpEekuIgeLSDcR2VXEdcYwdwNmo8yf20SAuGuTt73bigD0zo6WDggQiBxAokoeAd0e1oxhDUxOMoZ5HYSPBwI7EiDu8lU4IYAAdEKJOp4JEIg8o7O2odOM4cwdgGQBW/upYHgRBIi7fBpOCCAAnVCijmcCBCLP6KxuqGcEG/192bNuEZ8oIkeISLl8VH4UkR5cAWP1d4LxhRMg7vJlOCGAAHRCiTqeCRCIPKOjYQECmcSRKiKyVETe4+oXvhEIFEqAuMuH4YQAAtAJJep4JkAg8oyOhhCAAAQ8ESDuesJmXSMEoHUuj9ZgAlG0vBkNAhCAAHGXb8AJAQSgE0rU8UyAQOQZHQ0hAAEIeCJA3PWEzbpGCEDrXB6twQSiaHkzGgQgAAHiLt+AEwIIQCeUUlynRYsWcsstt8jhhx8uGzdulKlTp0qbNnoRx1+levXq8sgjj0iTJk1k3bp1MmLECLnxxhtly5YtjqgQiBxhohIEIACBwAgQdwNDmeqOEICpdm/xxrVt21Yef/xx6d27t0yaNElKly4tdevWleeff940LFmypHz++eeybNkyue6666RKlSry9NNPmzY33XSTI3IEIkeYqAQBCEAgMALE3cBQprojBGCq3Vu0caVKlZIFCxZI37595cknnyy0YrNmzeT111+XqlWryvLly02diy++WO666y6pVKmSbN68OSs9AlFWRFSAAAQgECgB4m6gOFPbGQIwta4t3rD69evLzJkz5dxzz5Urr7xSKleubFb7dKVv9uzZpvGtt94qrVu3lqOOOiqvs/3331/mz59vfk/rZysEomyE+HMIQAACwRIg7gbLM629IQDT6tksdp111lkyZswYWbhwoVx99dVmNfCaa66Rpk2bSs2aNWXVqlXy2GOPyX777Se6Epgp5cqVkw0bNkjz5s1l/PjxO4xSpkwZKVu2bN7vayBasmSJlC9fXtauXWspbcyGAAQgEB0BBGB0rJM8EgIwyd4rZO4DBw6UXr16FWtVrVq15Oijj5Znn31WLrroInOmT4uKt8WLF0ufPn1k6NChngSgbin369dvh/ERgCn70DAHAhCILQEEYGxdE6uJIQBj5Q7/k6lYsaJUqFCh2I7mzZsnxx9/vEyePFkaNGgg06ZNy6s/ffp0mThxohGBXraAWQH070N6gAAEIOCHAALQDz172iIA7fH1dpZqgNDEju7du+clgWgWsK4A3nzzzWZVMJMEotm/K1asMO0vvPBCueeee2TvvfeWP/74Iys9AlFWRFSIEQHeG46RM5iKZwLEXc/orGqIALTK3dsb+8ADD0j79u3lvPPOM2cBNQGkVatWolvEq1evzrsG5qeffpLrr7/eJIqMHDlSnnjiCa6Bsfi7SavpevvlQ3r3ZT4DF4nIVSLyUlqNxq5UEkAAptKtgRuFAAwcaXI61BU/PTN49tlniyZ3zJgxQ3r06CFff/11nhE1atQwF0E3btxY1q9fby6C1jOGXASdHD8z0+wEVPyN+7uargJmyta//6U9IjA7RGrEhgACMDauiPVEEICxdk/yJ0cgSr4P026BCr4FIlJNLz8vxFgVgYtF5AARyQjCtDPBvmQTIO4m239RzR4BGBVpS8chEFnq+ASZ3UhEpjiYb2MRmeqgHlUgkGsCxN1ceyAZ4yMAk+GnxM6SQJRY11kz8Y4iMtqBtZ1EZIyDelSBQK4JEHdz7YFkjI8ATIafEjtLAlFiXWfNxFkBtMbV1hhK3LXG1b4MRQD6wkfjbAQIRNkI8ee5JsAZwFx7gPGDJkDcDZpoOvtDAKbTr7GxikAUG1cwkWIIkAXM55EmAsTdNHkzPFsQgOGxpWcRIRDxGSSFQGH3AP4oIj24AiYpLmSefxMg7vIpOCGAAHRCiTqeCRCIPKOjYQ4I8BJIDqAzZOAEiLuBI01lhwjAVLo1PkYRiOLjC2YCAQjYQYC4a4ef/VqJAPRLkPbFEiAQ8YFAAAIQiJYAcTda3kkdDQGYVM8lZN4EooQ4imlCAAKpIUDcTY0rQzUEARgqXjonEPENQAACEIiWAHE3Wt5JHQ0BmFTPJWTeBKKEOIppQgACqSFA3E2NK0M1BAEYKl46JxDxDUAAAhCIlgBxN1reSR0NAZhUzyVk3gSihDiKaUIAAqkhQNxNjStDNQQBGCpeOicQ8Q1AAAIQiJYAcTda3kkdDQGYVM8lZN4EooQ4imlCAAKpIUDcTY0rQzUEARgqXjonEPENQAACEIiWAHE3Wt5JHQ0BmFTPJWTeBKKEOIppQgACqSFA3E2NK0M1BAEYKl46JxDxDUAAAhCIlgBxN1reSR0NAZhUzyVk3gSihDiKaUIAAqkhQNxNjStDNQQBGCpeOicQ8Q1AAAIQiJYAcTda3kkdDQGYVM8lZN4EooQ4imlCAAKpIUDcTY0rQzUEARgqXjonEPENQAACEIiWAHE3Wt5JHQ0BmFTPJWTeBKKEOIppQgACqSFA3E2NK0M1BAEYKl46JxDxDUAAAhCIlgBxN1reSR0NAZhUzyVk3gSihDiKaUIAAqkhQNxNjStDNQQBGCpeOicQ8Q1AAAIQiJYAcTda3kkdDQGYVM8lZN4EooQ4imlCAAKpIUDcTY0rQzUEARgqXjonEPENQAACEIiWAHE3Wt5JHQ0BmFTPJWTeBKKEOIppQgACqSFA3E2NK0M1BAEYKl46JxDxDUAAAhCIlgBxN1reSR0NAZhUzyVk3gSihDiKacaOQEkRaSgiVURkqYi8JyJbYzdLJhRHAsTdOHolfnNCAMbPJ6maEYEoVe7EmIgItBGRh0Sker7xFonIVSLyUkRzYJjkEiDuJtd3Uc4cARglbQvHIhBZ6HRM9kVAxd+4v3vQVcBMyaz+tUcE+uJrQ2Pirg1e9m8jAtA/Q3oohgCBiM8DAs4JqOBbICLVRCS/+MsvAheLyAFsBzuHamFN4q6FTvdgMgLQAzSaOCdAIHLOipoQaCQiUxxgaCwiUx3Uo4qdBIi7dvrdrdUIQLfEqO+KAIHIFS4qW06go4iMdsCgk4iMcVCPKnYSIO7a6Xe3ViMA3RKjvisCBCJXuKhsOQFWAC3/AAIyn7gbEMiUd4MATLmDc20egSjXHmD8JBHgDGCSvBXfuRJ34+ubOM0MARgnb6RwLgSiFDoVk0IlQBZwqHit6Jy4a4WbfRuJAPSNkA6KI0Ag4vuAgHsChd0D+KOI9OAKGPcwLWxB3LXQ6R5MRgB6gEYT5wQIRM5ZURMC+QnwEgjfg1cCxF2v5OxqhwC0y9+RW0sgihw5A0IAApYTIO5a/gE4NB8B6BAU1bwRIBB540YrCEAAAl4JEHfv88IHAAALHElEQVS9krOrHQLQLn9Hbi2BKHLkDAgBCFhOgLhr+Qfg0HwEoENQVPNGgEDkjRutIAABCHglQNz1Ss6udghAu/wdubUEosiRMyAEIGA5AeKu5R+AQ/MRgA5BUc0bAQKRN260ggAEIOCVAHHXKzm72iEA7fJ35NYSiCJHzoAQgIDlBIi7ln8ADs1HADoERTVvBAhE3rjRCgIQgIBXAsRdr+TsaocAtMvfkVtLIIocOQNCAAKWEyDuWv4BODQfAegQFNW8ESAQeeNGKwhAAAJeCRB3vZKzqx0C0C5/R25tJhBVq1ZN1q5dG/n4DAgBCEDANgIad5csWSLly5cn7trmfBf2IgBdwKKqewJVq1Y1gYgCAQhAAALREtC/eP/000/RDspoiSGAAEyMq5I7URWBSVr9y/zt2bZVS+y2b5Uan6fX5+pbxF9yf25GMXMEYBSUGSNRBGw9P4Pd9m2X4XP7fJ6oYMxkQyWAAAwVL50nkQA/FO36oWirv/X/TVttt9XuJMZj5hweAQRgeGzpOaEEbP3hgN12CV8E4G8kSSQ0RjPtYAggAIPhSC8pIlCmTBm58cYbZeDAgfLHH3+kyLLiTcFuu/ytXwM+t8/n1gQ0DM1KAAGYFREVIAABCEAAAhCAQLoIIADT5U+sgQAEIAABCEAAAlkJIACzIqICBCAAAQhAAAIQSBcBBGC6/Ik1EIAABCAAAQhAICsBBGBWRFSAAAQgAAEIQAAC6SKAAEyXP7HGJYHevXtLy5Yt5cgjjzQZv3vuuecOPVSvXl0eeeQRadKkiaxbt05GjBhhsoS3bNmSV7dRo0Zy//33y6GHHiqLFi2S22+/3dRLSjnkkEPknnvukeOPP95khn7xxRdy8803y5QpU/JMcMIhKfbmn2eLFi3klltukcMPP1w2btwoU6dOlTZt2qTe7oyB6u8ZM2aY/wf016xZs/JsP+yww2TIkCFSv359WbFihQwaNMh8J0kt++23n/muTzzxRKlcubJ5KWPUqFEyYMAA2bx5c2rtTqq/mHe4BBCA4fKl95gT6Nevn6xevVr23XdfOf/883cQgCVLlpTPP/9cli1bJtddd51UqVJFnn76aXn88cflpptuMtbtv//+8tVXX8mjjz4qTzzxhJx00kny4IMPGmH51ltvxZzAX9P79ttv5bvvvjPC9vfff5cePXpIt27d5KCDDpKff/5ZnHBIhKEFJtm2bVvjS/2LwKRJk6R06dJSt25def75503NtNqdH4N+q/oXABXC+QWg3gs5d+5cmThxorkSScXgk08+ab4NZZbEcuqpp8pZZ50lo0ePlu+//974Wm0ZOXKk+f9bSxrtTqKvmHP4BBCA4TNmhAQQ6Nq1qxFtBVcAmzVrJq+//rroe8bLly83llx88cVy1113SaVKlcyqwZ133mnEnv6AzBT9AbPHHntI8+bNY299hQoVZOXKldKwYUN5//33zXx32203837zySefLO+884444RB7QwtMsFSpUrJgwQLp27evETaFlTTand9OtU9Xrtu1aydff/31dgLwkksuMStjulKWWR1TIXjGGWdI7dq1k+buIud77bXXyqWXXmr+sqPFFrtT40AM8UwAAegZHQ3TRKAoAXjrrbdK69at5aijjsozV1f85s+fb35PVwd1y/DTTz+Vnj175tXR1TMVlCoCk1C++eYbee+998zqzqZNm8w/dUWkVq1aZoXUCYck2Jl/jrqtOXPmTDn33HPlyiuvNEJH/al2z54921RNo90ZBnvvvbd88sknRtDpXwBUDOdfAdQjDOXLl99uO7xx48YyefJk8xcl/S7SUPr372/+gqPfgxZb7E6D77DBHwEEoD9+tE4JgaIE4GOPPSZ6bkh/QGRKuXLlZMOGDWZ1b/z48Wb79KmnnjIrgZmif/bGG2+I1tVzZXEv1apVk5dfflmOPvpo2bp1q1nt1FVNFURanHCIu40F56dbgWPGjJGFCxfK1VdfbQTQNddcI02bNpWaNWvKqlWrUml3hoN+n9OmTTOrfPqNFxSAEyZMMH/R0RWxTNGVP10p1H/OmTMnaS7fYb666qciWFcB9fiGFhvsTrzjMCAQAgjAQDDSSZwI6DZVr169ip2SrmypcMuUNApANxxU/O20005GDOgZwAsuuMCsfOqqiJ5/TJIAdGq3it1nn31WLrroorwzbZoQsXjxYunTp48MHTo0UXbrt+zUdhW5HTp0EE1eUsGfdAHo1O78/8/rsQ5dvddEpwsvvDAvFiAA4xTNmUuYBBCAYdKl75wQqFixoui5tuLKvHnztsv6S+MWsFMOevZPk1V0W0/P/WWKJgAMGzbMnHdM0laoU7s141m3Mxs0aGBWwjJl+vTpJvFBRWCS7Nb5O7V97Nix0qpVK9m2bVue3ZoA8+eff8ozzzxjEoCStBXq1O7MWUZN5lLhp75WW/NzSJLdOQmwDJoaAgjA1LgSQ/wQyJYEoj8w9BoMLbpaoFdh6BkqvTpGt341g1KvEckU/SG61157JSIJ5LTTTjPbv7vvvrusX78+zwbd4tMfhrq6kkmGKI6DH/65aKvZnrrV3b1797wkEBVBugKoV4Vodmga7VbWeqWPnu/LFF0N078EaDKIXgmzZMmSvGSIffbZxwhDLbpCrJnTSU4CUVtV+OvWb5cuXcwKaP6SSQJJm925+H+MMeNNAAEYb/8wu5AJ6A9CFWq63amH/3U1TIteEaFiKHMNiN4Xdv3115tEAb0yQs8LFbwGRu9L02xSvWPs4YcfTsw1MLpaqmJPt8Nuu+02swWsIveqq64yW8B6J6ATDiG7KpTuH3jgAWnfvr2cd9555iygfgO6MpZJfkmr3QVhFrYFrAJRt0xVGOoqsF6Zot+3Jjsl9RoYFX+68qe+1r/05b/LU6870pJGu0P5n4dOE08AAZh4F2KAHwKavKFbQAWLZjuqINJSo0YNcxG0/p6KQl0V0zOGBS+CVjFRp04ds4KkmYVJugj6mGOOMas79erVM2cBNQtWxaAmuWSKEw5+fJGLtrripyucZ599tknY0dUvzYDWRIc02+1EAGqd/BdBa6awXgR9991358JVgYypom/48OGF9lWihP44/Kukze5A4NFJ6gggAFPnUgyCAAQgAAEIQAACxRNAAPKFQAACEIAABCAAAcsIIAAtczjmQgACEIAABCAAAQQg3wAEIAABCEAAAhCwjAAC0DKHYy4EIAABCEAAAhBAAPINQAACEIAABCAAAcsIIAAtczjmQgACEIAABCAAAQQg3wAEIAABCEAAAhCwjAAC0DKHYy4EIAABCEAAAhBAAPINQAACEIAABCAAAcsIIAAtczjmQgACEIAABCAAAQQg3wAEIAABCEAAAhCwjAAC0DKHYy4EIAABCEAAAhBAAPINQAACEIAABCAAAcsIIAAtczjmQgACEIAABCAAAQQg3wAEIAABCEAAAhCwjAAC0DKHYy4EIAABCEAAAhBAAPINQAACEIAABCAAAcsIIAAtczjmQgACEIAABCAAAQQg3wAEIAABCEAAAhCwjAAC0DKHYy4EIAABCEAAAhBAAPINQAACEIAABCAAAcsIIAAtczjmQgACEIAABCAAAQQg3wAEIAABCEAAAhCwjAAC0DKHYy4EIAABCEAAAhBAAPINQAACEIAABCAAAcsIIAAtczjmQgACEIAABCAAAQQg3wAEIAABCEAAAhCwjAAC0DKHYy4EIAABCEAAAhBAAPINQAACEIAABCAAAcsIIAAtczjmQgACEIAABCAAAQQg3wAEIAABCEAAAhCwjMD/B4oXRsLjKOBLAAAAAElFTkSuQmCC\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:>"
      ]
     },
     "execution_count": 59,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# To see if data has sense\n",
    "gdf = gpd.GeoDataFrame(coords, geometry=gpd.points_from_xy(coords.Longitud, coords.Latitud))\n",
    "world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))\n",
    "ax = world[(world.continent == 'South America')|(world.continent == 'Antartctica')].plot(color='white', edgecolor='black')\n",
    "gdf.plot(ax=ax, color='red')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}