animation_example.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "f40a1ce5",
   "metadata": {},
   "source": [
    "### Imports"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "fdfd49a3",
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "%matplotlib notebook\n",
    "import time\n",
    "#import multiprocessing\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "from cycler import cycler\n",
    "import geo.caustic_frame as caustic_frame\n",
    "import monte_carlo_simulation as mcs\n",
    "from bandstructure.delafossite_bandstructure import delafossite\n",
    "from matplotlib.animation import FuncAnimation"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "49eafcb1",
   "metadata": {},
   "source": [
    "### Simulation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "8eda4c61",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "10.0 T: [3811]\n"
     ]
    }
   ],
   "source": [
    "k = delafossite()\n",
    "bar_frame = caustic_frame.Frame('geo/bar.dxf')\n",
    "field = 10.0\n",
    "phi = 0\n",
    "n_inject = 1\n",
    "np.random.seed(42)\n",
    "bar_sim = mcs.Simulation(bar_frame, k, phi, field, 0)\n",
    "edge_to_collisions, trajectories = bar_sim.run_simulation(n_inject)\n",
    "print(f'{field:.1f} T: {list(map(len, trajectories))}')"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "43ab8527",
   "metadata": {
    "scrolled": true
   },
   "source": [
    "### Animation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "ed6171f2",
   "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",
       "\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 (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",
       "        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",
       "    rubberband_canvas.addEventListener(\n",
       "        'dblclick',\n",
       "        on_mouse_event_closure('dblclick')\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",
       "            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",
       "// 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.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 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\", \"pgf\", \"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.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",
       "    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,iVBORw0KGgoAAAANSUhEUgAAAyUAAAG/CAYAAABG9niTAAAgAElEQVR4nOzdeXgUVaI28A4RdVgERC864q0QISIogoiIDhJcceTiuKAzKjoOfgo4M4KKFUgCAUVERAEVhMgioEjYobKHLaxZCGETkrAkYQlJyE62Tne/3x+dNIQEUt2p7qrqfn/P088MbXXX4XDqVL1ddc4xgIiIiIiISEUGtQtARERERESejaGEiIiIiIhUxVBCRERERESqYii5wtq1a3H//fejVatW6Ny5M3777Te1i0RERERE5PYYSmpt2bIFd911F3bs2AGTyYT8/HxkZGSoXSwiIiIiIrfHUFLr8ccfx4IFC9QuBhERERGRx2EoAWAymdCyZUvMmDEDfn5+uPPOO/HWW2+hoKBA7aIREREREbk9hhIA586dg8FgQO/evXHmzBkUFxfjxRdfxGuvvdbo9iEhIfDy8rK9DAZDvT974stgYFMiIiIiIsfwShJAUVERDAYDfv75Z9t7ycnJaN26NSwWS5Of9/LycmbxdIF1QERERESOYiipdffdd2PRokW2PycnJ6NVq1YMJTKxDoiIiIjIUQwltaZMmYI+ffogJycHZWVlePnll6/5+NbVeEHOOiAiIiIixzGU1KqpqcFHH32EDh064LbbbsObb74pe6A7L8hZB0RERETkOIYSBfCCnHVARERERI5jKFEAL8hZB0RERETkOIYSBfCCnHVARERERI5jKFEAL8hZB0RERETkOIYSBfCCnHVARERERI5jKFEAL8hZB0RERETkOIYSBfCCnHVARERERI5jKFEAL8hZB0RERETkOIYSBfCCnHVARERERI5jKFEAL8hZB0RERETkOIYSBfCCnHVARERERI5jKFEAL8hZB0RERETkOIYSBfCCnHVARERERI5jKFEAL8hZB0RERETkOIYSBfCCnHVARERERI5jKFEAL8hZB0RERETkOIYSBfCCnHVARERERI5jKFEAL8hZB0RERETkOIYSBfCCnHVARERERI5jKFEAL8hZB0RERETkOIYSBfCCnHVARERERI5jKFEAL8hZB0RERETkOIYSBfCCnHVARERERI5jKFEAL8hZB0RERETkOIYSBfCCnHVARERERI5jKFEAL8hZB0RERETkOIYSBfCCnHVARERERI5jKFFAsy7IQzoAk9sB0+9W7zW5nbUcatUBEREREXk0hhIFNOuCfHI7YPItKoeSW6zlUKsOiIiIiMijMZQooFkX5HXBQE0KlIGhhIiIiIgcxVCigOZckBcHtENxQPPuUjSX0N4bQnvvZn0HQwkREREROYqhRAHNuSDPGncLssbdomBp7Nf+Ty3Q/k8tmvUdDCVERERE5CiGEgU054K8ffv2aN++vYKlUacMDCVERERE5CiGEgUwlDCUEBEREZHjGEoUwFDCUEJEREREjmMoUQBDCUMJERERETmOoUQBDCUMJURERETkOIYSBTCUMJQQERERkeMYShTAUMJQQkRERESOYyhRAEMJQwkREREROY6hRAEMJQwlREREROQ4hhIFMJQwlBARERGR4xhKFMBQwlBCRERERI5jKFEAQwlDCRERERE5jqFEAQwlDCVERERE5DiGEgUwlDCUEBEREZHjGEoUwFDCUEJEREREjmMoUQBDCUMJERERETmOoUQBDCUMJURERETkOIYSBTCUMJQQERERkeMYShTAUMJQQkRERESOYyhRAEMJQwkREREROY6hRAEMJQwlREREROQ4hhIFMJQwlBARERGR4xhKFMBQwlBCRERERI5jKFEAQwlDCRERERE5jqHkKnl5eejYsSP69u0r+zMMJQwlREREROQ4hpKrvPnmm3jiiScYSuzEUEJEREREjmIouUJ0dDQGDhyIxYsXM5TYiaGEiIiIiBzFUFKrvLwc9957L44cOYIlS5YwlNiJoYSIiIiIHMVQUmv8+PEICAgAgCZDSUhICLy8vGwvg8HxamQoISIiIiJPx1ACIDU1Fd26dUNFRQWApkPJ1XinhKGEiIiIiBzHUALgu+++Q6tWrdCpUyd06tQJt9xyC2644QZ06tQJBQUFTX6eoYShhIiIiIgcx1AC63iSnJwc22v27Nl48MEHkZOTA4vF0uTnGUoYSoiIiIjIcQwljeDjW/ZjKCEiIiIiRzGUKIChhKGEiIiIiBzHUKIAhhKGEiIiIiJyHEOJAhhKGEqIiIiIyHEMJQpgKGEoISIiIiLHMZQogKGEoYSIiIiIHMdQogCGEoYSIiJyHwaDAV5eXnzxxZeTXwbD5SjCUKIALy+GkubUARERkZbwnEbkGlceawwlCmAoYQdORETug+c0ItdgKFEYQwk7cCIich88pxG5BkOJwhhK2IETEZH74DmNyDUYShTGUMIOnIiI3AfPaUSuwVCiMIYSduBEROQ+eE4jcg2GEoUxlLADJyIi98FzGpFrMJQojKGEHTgREbkPntOIXIOhRGEMJezAiYjIffCcRuQaDCUKYyhhB05ERO6D5zRtOl9cgZijFzAnLh0frUzB3xfsxZPfbMOAL+PwyLRY9PsiFk9+sw2vzt+NMSv2Y1ZMGjalnsPp/EuwWCxqF58awVCiMIYSduBEROQ+eE7TBqPJjC3HLiBw/SEM+norBFFq9OU7IRzdJkag28SIa27Tf1ocPlqZgs0Hz+FSVY3afzWqxVCiMIYSduBEROQ+eE5Tj8Viwf6sQgSuP4TeU6LrBYshs+MRsPYQlu3NxN6TF5F58RIqqk31Pl9RbcK5ogokZxbg98QsTNl0FMN+2AXfCeG27/ELjMCYFfuxKyOfd1BUxlCiMIYSduBEROQ+eE5zPaPJjA0HzmLY9ztt4aHrxHC890sSNhw4i/yyqmZ9f1lVDWKOXsC4VQdw/+Qo2z4Gz9yGRTtPobyad0/UwFCiMIYSduBEROQ+eE5zHaPJjJUJWXhs+pZ6d0R+3ZeFovJqp+yzqsaEDQfOYvj8PbZ9PjQ1Bj9uy0BppdEp+6TGMZQojKGEHTgREbkPntOcz2y2YMOBs/XGioxcmojdJ1z7SNXxnFKM+/2A7fGuB6dEY9HOU6iuMbusDJ6MoURhDCXswImIyH3wnOZcyZkFeGFuvC2MvLM4AYfPFqtapsyLlzB+dSq6BFjLNOjrrYg8fJ5jTpyMoURhDCXswImIyLmMRiNGjx6N9u3b49Zbb8X48eOve8E4f/58dO3aFa1atYKvry+2bdsme188pzlHbkklxv1+wBZGXp63G4mnC9QuVj3pF0rx7pJEWxnfCN2L0/mX1C6W22IoURhDCTtwT/D999+jb9++uPHGG/HKK69cd9uSkhK8/vrraNOmDe644w588803LiolEbmrSZMmoV+/fsjNzUVWVhb8/Pwwd+7cRrddsmQJ7rvvPqSmpsJiseDs2bPIzs6WvS+e05RVYzIjNP4kegRHQhAl9PsiFutTzmr6LsSujHw8++0O22xdP2zN4CNdTsBQojCGEnbgnmDt2rVYv349PvzwwyZDydtvv41hw4ahpKQEhw4dwu23345Nmza5qKRE5I46d+6MjRs32v68cOFC9O7du8F2ZrMZd999N6Kjox3eF89pyjlyrhhD5+60zab1ZcQfKNPJOiFGkxk/bM2AX6B1/ZNnv92B5MxCtYvlVhhKFMZQwg7ck0yePPm6oaS8vBw33ngjDh48aHtv4sSJ+Nvf/uaK4hGRGyosLITBYEBmZqbtvcTERLRs2bLBr+3Z2dkwGAyYO3cuBEFA586d8Z///AcVFRWy98dzWvNVVJvwZcQftgHkr8zbjYzcUrWL5ZDT+ZfwZug+CKIEnwAJgesPoYSzdCmCoURhDCXswD1JU6EkJSUF3t7eMJsv3+YOCwtD165dZe/DYDDAy8uLL774csHLYND+pUBd0CgqKrK9l56eDoPBgMrKynrb7t69GwaDAc8++ywKCgqQk5ODRx99FJ999tk1vz8kJER3daJlO9PzMXCGdVat+ydFYfneTJjN2n1USw6LxYJ1KWfQZ2qM7RG0iEMcCN9cXl4MJYq6skLtxVBCetNUKImPj0e7du3qvRcTE4NOnTpd8zO8ICBSjx7677o7JVlZWbb3kpKSGr1TcuDAARgMBsTFxdneW7NmDe677z7Z+9NDnWhR4aVqfLwq1TZI/P1lScgprmz6gzpScKkan4Sl1pvG+FyR/LtwVB9DicIYStiBexK5d0quvFBYvXq1XXdK2J6IXEcvx1vnzp3rjU0LDQ1tdExJeXk5brrpJmzZssX23po1a9C9e3fZ+9JLnWiFxWLB2v317yJEHj6vdrGcandGvm2NlfuCI/HzzlMw6fxukBoYShTGUMIO3JPIHVNy6NAh23uBgYF2jSlheyJyHb0cb8HBwejfvz/y8vKQnZ2N7t27X3P2rXfffRdDhgxBcXEx8vPz8dhjj1338a2r6aVOtOBU/iW8EbrXNt5i4rpDKK7wjPEWlUYTZkYdxz2142aGzt2JQ2fUXW9FbxhKFMZQwg7cE9TU1KCyshKBgYF46aWXUFlZierq6ka3HTFiBF588UWUlpbi8OHD6NSpk12zb7E9EbmOXo43o9GIUaNGoV27dujQoQM+/fRT2x3ZIUOGYNq0abZty8rKMGLECNxyyy244447ONDdCaprzPh+Szq6cWYqpF0oxSvzdkMQJXQJkDB54xGUciC8LAwlCmMoYQfuCSZPngyDwVDvNWjQIAANLwhKSkrw2muvoU2bNujUqZPd65SwPRG5Do+3hlgn15d0ugDPfLvdtobHj9syYDR59hoeZrMFKxOy0CskGoIoof+0OK4ILwNDicIYStiBk7LYnohch8dbQ6yTxhWVVyNg7SHbIO+3ft6HzItc7fxK+WVVGHvFqvX/WpKIM4XlahdLsxhKFMZQwg6clMX2ROQ6PN4aYp3UZzJb8Ou+LPSeYr0L8NDUGGw4oO0V2dW2KyMf/jO3QRAl3BsUgblx6ag0mtQuluYwlCiMoYQdOCmL7alxxRVG7M8qxIYDZ7Fo5ynMjDqOgLWHMPb3A/j3bykYtTwZ7/2ShP/3SxL+uzIF4pqDmLzxCL6KPIY5celYuOMklu05jbCkbGw+eA5xf1zArox8JGcW4si5YpzMK8O5ogoUXqpGRbUJVTUmGE1m1JjMMJstilyAWCwWmMwWVNeYcamqBvllVcguKEdGbikOnSlGwqkCbE/LQ+ThHGw4cBYrE7KwaOcp/LA1AzOjjmPq5qOYsM76dx69Ihmjlifjg2WXX+N+P4DgDYcxI/IYFuw4gc0Hz+FAdhEKLzU+/ol4vDWGdXJZcmahbUX2uoHsReU8nuSoNJrwbUyabdzN419tQeThHIa5KzCUKIyhhB04KYvtyfrL5KEzxQiNP4mRSxPR74tY2+MAar+6BEjwnRCOrhPD0W1iBPwCI3BvUATuC45Ej+BI3D8pCj0nRaF7UCT8AiPQdWI4fCeEwydA3XI/Nn0LPliWjPnbT+DouRJeGNTi8dYQ6wTILa2st+bIK/N24/BZzizliOyCcnywLNlWl2+E7kXaBX2ubq80hhKFMZSwAydleXJ7OnimCCGbjjQIIV0nhuPpWdsxankypkccQ2j8SaxPOYsdaXlIOl2AlKxCHD5bjD/Ol+Dw2WIkZxZgV0Y+Yo9ewKbUc1iVlI1f9pzGgh0nMDs2HV9G/IHJG49AXHMQ/12ZgveXJWHEogQM/2kP/u/7nXjm2+148pttGDxzG/xnbsMTX2/FwBlb8fhXW/DY9C0Y8GUcHv0yDo9Mi0W/L2LR9/NY9P08Bn2mxuDBKdHoMzUGfT+PwcNfxKL/tDgM+DIOj03fgoEztmLQ11sxeOY2PPvtDgz7YRdeX7AH/1ycgNErkjFu1QFMXHcIUzcfxcyo4/h+Szp+3nkKvyVkYX3KWUQezsGOtDwkni7AwTNFOHy2uN4r8XQBth7LxcbUc1iy6xS+kI5i9IpkPPvtDnS5KhT1nxaHzzcfxbGcErX/2VXlycfbtXhynRhNZoTGn0TPSVEQRAkPfxGLdSlnGOIVsCsj3zZBgO+EcEzeeATF5Z49SxdDicIYSjy7AyfluXN7MpvNCAsLg7+/PwRBgL+/P37/fRWk1LN4uXZKybq7EX9fsBc/bM1A0ukCVNd49sw2SiivrkFyZgF+2JqB4fP3wLd2bQFBlPDyvN3YcuyCR154ufPx5ihPrZMdaXl4apb1ovmeCeGYFv4Hp7ZVWI3JjCW7TuGBydbQ12dqDH7dl+WxCy8ylCiMocRzO3ByDndtT2azGW+99Ra8vb3h5eUFg8EALy8veLVogdY9B+N/P9uEp2dtx9Ldp5FfVqV2cd1eflkVQuNP4unaizBBlPB/3+/0uLUW3PV4aw5Pq5OM3DK8uySx3qxaGbllahfLrV0sq8KEdYdsj7X+dU48Ek8XqF0sl2MoURhDied14ORc7tqewsLC4O3t3WC9F4PBAK8WLRAy+2eP/KVebWazBbFHL2DYD7tsF2WfhqWixEN+IXbX4605PKVOisqrEbLpiG1F8kFfb0XMUc+8Y6iWw2eL8er8y3fJ//NbCs4Xy1/oU+8YShTGUOI5HTi5hru2J39/f9sdkgahxMsL/v7+ahfRo1ksFqxLOYOHa8fzPP7VFo+4a+Kux1tzuHudGE1mLN19Gg/WTvF7/+QohMaf5GOiKrFYLNhw4Cz6T4uDIEroHhSJ77d4xhTCDCUKYyhx/w6cXMtd25MgCI0GkrqXj4+P2kUkWKdeHvPrfttz9WFJ2WoXyanc9XhrDneuk23Hc23jRroESAhafxgX+bioJpRX12Bm1HHbFMJ/mbEF0kH3XhWeoURhDCXu3YGT67lje7JYLPDt9QgMBt4p0QOLxYLfErJsj7XMiklz2wsDdzzemssd6yQjtwzvLE6wPSb0Zug+HM/htLRalHWxHP/vlyTbv9Wr83fj4JkitYvlFAwlCmMocc8OnNTjbu3JYrFg0obDuO3FABi8WjQaSry9vREWFqZ2Uekqu0/k22bJmRV9XO3iOIW7HW9KcKc6KSqvxuSNR2yzzQ2euQ1xf3DciB7sPpGP52fH28LJuN8PuN14E4YShTGUuFcHTupzt/Y0KyYNgijhoSnRGPbq3xvMvuXt7Y0RI0bAbObz3Fp0+GyxLZjMiUtXuziKc7fjTQnuUCfG2qlne4VYx408MDkKP+88xXEjOmMyW7AqMds21u3eoAjMiklDeXWN2kVTBEOJwhhK3KMDJ+1wp/b0674s20DS4zml9dYp8fHxgb+/P8LCwhhINC41uwj31y4mt3b/GbWLoyh3Ot6Uouc6sVgs2HosF09+s822SF/whsMouFStdtGoGS5V1eCb6OPwqx1v8si0WIQlZet+fROGEoUxlOi7AyftcZf2tO/kRfhOCEe3iRHYe/Ki2sWhZtqdkW/99wyMQEqW+8zK5S7Hm5L0WifJmYUY/tOeeuuNpF3guBF3craoAv9dmWL7N3722x26nsaZoURhDCX67cBJm9yhPeWWVtput69Jdq9f1j3Zsr2Ztl8pC93kl2d3ON6Uprc6SbtQiveuGBj93Hc7sPVYrm4vVKlpKVmFeO2KAPrSj7uwT4c/fjGUKIyhRH8dOGmb3tuT2WzBG6F7IYgSJqw7pHZxSGHjV6dCECV8sCzZLS769H68OYNe6uRMYTk+XpVqWxX88a+2YF3KGd0/0kPyWCwWbE/LqzcY/p3FCThyrljtosnGUKIwhhL9dOCkD3pvT8trf01/7rsdHrH4lae5VFWDJ77eCkGU3GINE70fb86g9Tq5WFaFKZuOottE6/iCh6bGYMmuU6iqYX/jicxmCzamnrP1S3Urw5/Ov6R20ZrEUKIwhhLtd+CkL3puT2cKy9EjOBL3TAjH4bP6+bWK7LM/qxC+E8LRc1IUcoor1S5Os+j5eHMWrdZJSaURs6KPo0dwJARRQo/gSMyOTUdZlXvMxETNYzSZsXxvpu3RYd8J4QhYexDnirQ7jTBDicIYSrTbgZM+6bU9WSwWvL3IujjZN266pgVdNj3iGARRwoe/7le7KM2i1+PNmbRWJ5VGExbsOIEHp1in9+0WGIGpm49yJXZqVEW1CT9uy7BNB90tMAIhm44gr1R77YWhRGEMJdrrwEnf9NqeYo5egCBKeGrWdj5G4QHKq2vw2PQtEEQJO9Pz1S6Ow/R6vDmTVurEaDJjxb5MPDLt8i/f4hpt//JN2lFcUf/OWvegSMyIPIbicqPaRbNhKFEYQ4l2OnByD3psT1U1JtvzvPHpeWoXh1wk8nCObZVsvQZRPR5vzqZ2nZjNFmw4cLbeGIEPf92Pk3llqpaL9OliWRW+kI7a1ji5f3IU5sal45IGHvtjKFEYQ4n6HTi5Fz22p3nbTkAQJYxcmqR2UciFLBYL3llsfWRv8a5TahfHIXo83pxNrTqxWCyI++MCnvtuhy2M/HNxAsenkSIulFQiaP1hdJ0YDkGU0O+LWKxSeQFGhhKFMZTwpEbK0lt7ulhWhR7Bkeg2MUIXs52QstIulMInwDoDkhZ+ebSX3o43V1CjTtIulOL1BZfXnXh1/m4knCpweTnI/WUXlGPs7wdsbe352fHYnaHOI6gMJQpjKOFJjZSlt/Y0ZdNRCKKEqZuPql0UUsm4VdYT/Ny4dLWLYje9HG9GoxGjR49G+/btceutt2L8+PFNrhNz7Ngx3HTTTXjllVfs2pcr66S00ojPNx+F74TwywsfHufCh+R8qdlFGD7/chAetTwZF0pcO5sgQ8lVqqqq8N5778HHxwdt2rRB9+7dsWTJEtmfZyjRz0mN9EFP7elsUQW6TYxAj+BIzoTjwbILytF1YjjunxSlu5Xe9XK8TZo0Cf369UNubi6ysrLg5+eHuXPnXnN7i8WCgQMH4oknntBsKAk/dB79aqdvvX9SFJbsOoUak9kl+yYCrMdJ5OHzePyrLbZ2uGxvJswueqSLoeQqly5dQnBwME6cOAGLxYK9e/eiffv22LJli6zPM5To56RG+qCn9vTZ6oMQRAnfxqSpXRRS2aQNhyGIEr4M/0PtothFL8db586dsXHjRtufFy5ciN69e19z+wULFuDNN9/E5MmTNRdKisqr8Z/fUmy/UH+8KlWT07WS5yivrsGX4X/Y7ti9On83sgvKnb5fhhIZXnrpJUyZMkXWtgwl+jmpkT7opT2dyCuD74Rw9J4SjdJK7UyxSOrIK61C96BI+AVG6GpBRT0cb4WFhTAYDMjMzLS9l5iYiJYtWzb6mFNOTg66dOmC3NxczYWSPScu2qb4fWz6Fuw+od/ppMn9HDlXjKFzd9pm6dqUes6p+2MoaUJlZSXuuusurFmzRtb2DCX6OKmRfuilPY35dT8EUcLCHSfVLgppxNdR1gUVA9YeUrsosunheMvOzobBYEBRUZHtvfT0dBgMBlRWNgyAw4cPx08//QQAskJJSEgIvLy8bC+DQfnLI4vFgoU7Ttp+if4kLJU/ZpAmGU1mfB11DD4B1jt5E9YdgtFJjxUylFyHxWLBm2++CX9/f5jNjf8DKNl5MZQQNaSH9nT4bDEEUUL/aXGoNOpzfQpSXnGFEb1CouE7IRyndDITmx6Ot7o7JVlZWbb3kpKSGr1TIkkSHn30Udv7WrhTUmk04cPaHzHuDYrA+pSzin4/kTPsOXHRNubp7wv2oqhc+fFyDCXXYLFY8MEHH+Dhhx9GcbH8OcF5p0QfJzXSDz20p7q1KX7dl9X0xuRR5m+3rlnz799S1C6KLHo43gDrmJJNmzbZ/hwaGtromJKPPvoIbdu2RadOndCpUye0bt0aN910E3x9fWXvS8k6Kak04rWfrDMc/WXGFhw9V6LYdxM5W05xJV6YGw9BlOA/c5vi40wYShphsVgwevRo9OnTB4WFhXZ9lqFEPyc10gett6eEUwUQRAmDvt7qtFvapF8V1SbbmAE9LHqn9eOtTnBwMPr374+8vDxkZ2eje/fujc6+VVpaipycHNvrk08+wdChQ5Gbmyt7X0rVSVF5NZ6fbb2gG/b9ThTobGY2IsA6CP6DZcm2cVCZF5W7C8xQ0ogxY8agV69euHjxot2fZSjRz0mN9EHL7cliseDV+bshiBI2HOAjGNS4FfsyIYgS3lmcoHZRmqTl4+1KRqMRo0aNQrt27dChQwd8+umntke0hgwZgmnTpjX6ObUe3yqrqsGwH3ZBECW8GbpPlwtrEtUxmS34eFWq7bHlM4XK3DFhKLlKZmYmDAYDbrrpJrRu3dr2+uCDD2R9nqFEPyc10gctt6etx3IhiBKGzI532TzupD9GkxmDvt4KQZSw76T9P3a5kpaPN7U0t06MJjP+sXAvBFHC8J/2oKKa485I/8xmCz4NswaTRyZvwsCnnoMgCPD390dYWNg1x2JfD0OJwhhKeFIjZWm1PZnNFgypfRRjy7ELaheHNG5j6jkIooSX5+3W9OrcWj3e1NTcOpm88QgEUcILc+M5wxa5lWpjDfp8tACCKKHT61/A4NUCXl5e8Pb2xogRI+wOJgwlCmMo4UmNlKXV9rRJJxeZpA1ms8U2niD2qHZDrFaPNzU1p07WpZyBIEro+3kMzhVVKFgqIvWFhYXhhptb485/zoEgSmj32N9hMBhgMBjg7e2NsLAwu76PoURhDCU8qZGytNiejCYz/Gdu08XjOKQd245bH/d79tsdMGn0cT8tHm9qc7ROzhdX4P5JUegSwH6C3JO/vz+8vLxwQ4c/4+6xYfjf8Rtx0109YDAY4OXlBX9/f7u+j6FEYQwlPKmRsrTYnn5LyIIgSnh7kfYHLpN2WCwWDK+dDnbt/jNqF6dRWjze1OZInVgsFry9yDpV+FeRx5xQKiL1CYJguzPSuudgCKKEO0fOg6HFDTAYDPDx8bHr+xhKFMZQwpMaKUtr7amsqgYPf6GfKV5JW5IzC2xTaWpxoU2tHW9a4EidRB/JgSBKeGrWdlTVaO/fmUgJdXdK6oLJ/7w2FYIo4Zb+r8m3hAoAACAASURBVPBOiRYwlPCkRsrSWnv6Jvo4BFHCuFUH1C4K6dToFdY5/r/fkq52URrQ2vGmBfbWSY3JjMHfWB/v3JGW56RSEakvLCwM3t7etlByQ4c/438/WY+7x4bhhtbtOKZEbQwlPKl5CqPRiNGjR6N9+/a49dZbMX78+GsO+B40aBBuvPHGetNsV1VVydqPltrT2aIK+AVG4N6gCJwv5qBVckx2QTm6TYzAfcGRuFBSqXZx6tHS8aYV9tbJ6mTr4Pa3ft7npBIRaYPZbMaIESPg7e1tu2Ny69MfQBAlDBzzFWffUhtDCU9qnmLSpEno168fcnNzkZWVBT8/v0ZXVAasoeT77793aD9aak8frUyBIEr4NiZN7aKQzn0VeQyCKOGTsFS1i1KPlo43rbCnTiwWC577bgcEUUJyZoETS0WkDWazGWFhYfD394ePjw/+8uxQ+AZsxv2To1Bm5yKhDCUKYyjhSc1TdO7cGRs3brT9eeHChejdu3ej27pDKDmQXWRdJGpaLMqruRozNU9ppRF9P7eOTTp0Rjtjk7RyvGmJPXWy+0Q+BFHCsB92capw8lif1C6q+Htill2fYyhRGEMJT2qeoLCwEAaDAZmZmbb3EhMT0bJly0ZPxIMGDULHjh1x6623om/fvli7dq3sfWmhPVksFrw8bzcEUcLqZG3OmkT6s7J2FrdX52tnrRstHG9aY0+d1F2MaXV2NSJXqJvQ428/7rLrcwwlCmMo4UnNE2RnZ8NgMKCoqMj2Xnp6OgwGAyorGz4jv3fvXpSUlMBoNGLTpk1o3bo1tm/f3uh3h4SEwMvLy/YyGNTvmjYfPGdbkdms0fUlSH9MZguG1C6oGH7ovNrFAcD+uzFy66SqxoT7J0fBLzDC7sdWiNyJxWKxTfaQXVAu+3MMJQpjKOFJzRPU3SnJyrp8azYpKemad0qu9v777+PDDz+UtS+121Ol0YTHv9oCQZSwlwugkcLqHvd5/CttTBGs9vGmRXLrZOsx6+KYo5YnO7lERNo3o3bc3NLdp2V/hqFEYQwlPKl5is6dO2PTpk22P4eGhl5zTMnVRo0ahTFjxsjaVu329P2WdAiihPeXJalaDnJf7y9LgiBK+HFbhtpFUf140yK5dfL55qMQRAlhSdlOLhGR9tU9wjXCjkWGGUoUxlDCk5qnCA4ORv/+/ZGXl4fs7Gx079690dm3ioqKEBERgYqKCphMJkRERKB169aIjY2VtR8129O5ogp0D4pEt8AIu25BE9kj8+IldJ0Yjh7BkcgtVXeKYPbfDcmtk+drH8U7U8i+gshktqBHcCR6ToqCSeZjzwwlCmMo4UnNUxiNRowaNQrt2rVDhw4d8Omnn9oe3RoyZAimTZsGAMjLy0O/fv3Qtm1b3HLLLejduzdWrlwpez9qtqd//2adAvib6OOqlYE8w5fhf0AQJXy2+qCq5WD/3ZCcOqmoNsEnQMJj07e4oERE+vCPhXshiBKO55TK2l6zoSQhIQFBQUEYOnQoHn/8cQwdOhRBQUHYt0/bixExlPCkRspSqz0lnLLeen70yzhOAUxOV1JpxENTY+ATIOHwWfWmCGb/3ZCcOkmtnTJ85FI+5klU5+so67iSVYnyHmnUXCiJiYlBnz594Ovri5EjR2LWrFkIDQ3FrFmzMHLkSPj6+qJPnz6IiYlRu6iNYijhSY2UpUZ7unJWpI2p51y+f/JMK/ZlQhAlvPbTHtWmCGb/3ZCcOlmVmA1BlDAzindVieqsSzkDQZQwPeKYrO01F0qGDBmC+Pj4626zc+dOPP/88y4qkX0YSnhSI2Wp0Z7qLg6Hz1fv4pA8T43JbFsNfPNBdcIw+++G5NTJtNrH7/gjBtFlKVmFEEQJHyyTNyOd5kKJ3jGU8KRGynJ1eyoqr0bvKdHoEiDhyDntrLRNnmHvyYuqPjbI/rshOXVSN/4sObPQBSUi0oei8moIooTnvtsha3tdhBKLxYLy8vJ6L61iKOFJjZTl6vY0acNhCKKEiesOuXS/RHU+/HW/ahMssP9uSE6dDP9pD2feImrEvUER6Pu5vCEXmg4lCQkJeOihh9CyZUu0aNECLVq0gJeXF1q0aKF20a5J7VBiNgNhYYC/PyAI1v8NC7O+76oy8KRGSnJlezqWUwLfCeHoFRKNgkvVLtsv0ZXOF1+eijrromsvctl/NySnTgZ9vRWCKKGqRv0FMIm0ZMCXcfCdEA6zjGmBNR1K7r33XgQFBeGPP/5AZmZmvZdWqRlKzGbgrbcAb2/AywswGKz/6+0NjBghL5gwlJDWuKo9WSwW/H2BdfrCX/acdsk+ia7lh60ZqszmxP67ITl18tDUGHQPinRBaYj05a9zrJPGFJcbm9xW06Gkffv2uhtkqmYoCQuzBhCDoeHL29v6351dBoAnNVKWq9pT+KHztmdfa0wyby0SOUml0YSBM6y/vm87nuuy/bL/bkhOndwXHIk+U7U5KyiRml5fYH20Uc4CxJoOJaNGjcLGjRvVLoZd1Awl/v6X75Bc/fLysv53Z5cB4EmNlOWK9lRRbcJj07dAECXsPpHv9P0RyRF79AIEUcLgmdtQXeOaoKyX/ttoNGL06NFo3749br31VowfP77RHzFzc3Pxxhtv4K677kLbtm3Rp08fSJJk176aqhOLxcKFE4mu4a2f90EQJZzOv9TktpoOJYWFhejatSsGDhyI4cOH13tplZqhRBAaDyR1Lx8f55cB0M9JjfTBFe3pu9g0CKKEMSv2O31fRHJZLBa8szgBgihh3rYTLtmnXvrvSZMmoV+/fsjNzUVWVhb8/Pwwd+7cBtudPHkSM2fOxJkzZ2A2m7Fp0ya0bt0ax4/Ln0SgqTqpNJogiBKe/GabvX8NIrdX14dl5JY1ua2mQ8nQoUPRs2dPfPbZZwgJCan30ireKdHPSY30wdnt6WReGboFRuDeoAjOnEOacyr/ErpNtLZPOY8/NJde+u/OnTvXe5Ji4cKF6N27t6zP9unTB7/88ovsfTVVJxXV1lDy9Kztsr+TyFOMXJoIQZRwLKekyW01HUpat26NkpKm/xJawjEl+jmpkT44sz1ZLBb8Y+Fel/4STWSvOXHpEEQJIxYlOH2cpR7678LCQhgMhnqT3iQmJqJly5ZN1k9ubi5uvvlmJCXJn0CgqTqpqrGGksG8U0LUwAfLkiGIEg6fbXrdL02Hkr59++LcOX2tjqr27FsjRnD2LXIvzmxPa/efsQ1uN3JwO2lUdY0ZT83aDkGUsOHAWafuSw/9d3Z2NgwGA4qKimzvpaenw2AwoLKy8pqfq6qqwuDBg/H2229f9/tDQkLg5eVlexkM1788qjGZIYgSBn291b6/CJEHqLtTcjyntMltNR1KZsyYgYceegiLFi1CeHh4vZdWaWmdEh8frlNC+ues9lRwqRp9psbAJ4CrMJP2JZ4ugCBKeGhqDIrKnbeGjh7677o7JVlZWbb3kpKSrnunpLq6GkOHDsULL7yA6mr76k/OQHdBlPCXGRzoTnS1N0KtTyPIWXNJ06HEx8en0VeXLl3ULto1qR1KmouhhLTGWe3p07BUCKKEwPVcuZ30IWDtIQiihM9WH3TaPvTSf3fu3BmbNm2y/Tk0NPSaY0qqq6sxbNgwPPfcc6iqqrJ7X3LqpEuAhAFfxtn93UTu7qUfd0EQJeSWXvsuZh1Nh5LG1NTUYPbs2WoX45oYSvRzUiN9cEZ72nPiIgRRwsNfxKKksukFnYi0oLjCiIe/iIUgSth78qJT9qGX/js4OBj9+/dHXl4esrOz0b1790Zn3zIajXjxxRfx1FNPXffRruuRUyc9giPRe0q0Q99P5M6GzLYunlgq41yr6VBy9TTAw4cPx3PPPYeOHTuqXbRrYijRz0mN9EHp9lRVY8Lgb7ZBECVIB88r+t1Ezrb54DnboOpKo0nx79dL/200GjFq1Ci0a9cOHTp0wKeffmp7dGvIkCGYNm0aAGD79u0wGAy4+eab0bp1a9ur7r/LIadOHv4iFn6BEY79ZYjcmP9M6/lWzrhNTYeSq6cBHjt2LDp37oyFCxeqXbRrYijRz0mN9EHp9vRtjHVNkn8udv5MRkRKs1gseHeJdeDorJg0xb+f/XdDcurkia+3QhAlmMzsU4iu1H9aHLpOlDcWXNOhpDEZGRno2bOn2sW4JoYSntRIWUq2p4zcMnSbGIHuQZEuWfOByBnOFJaje1Akuk4MR0Zu0zPa2IP9d0Ny6sSeR1SIPMkDk6PwwOQoWdvqLpTU1NTgjjvuULsY18RQwpMaKUup9mQ2WzD8pz0QRAkLd5xU5DuJ1BIafxKCKOHV+bthVvDXefbfDcmpk5fn7YYgSrhQ4ti4FSJ31W1iBB6ZFitrW02HkqNHj9Z7JScnY8yYMRgzZky997WEoYQnNVKWUu3p98QsCKKEv86JRw3XJCGdqzGZ8cJc66/zKxOymv6ATOy/G5JTJ2/9vA+CKOFU/iUXlIhIH+xdw0fTocTLywstWrSot4jR1a8WLVqoXcx6GEp4UiNlKdGe8suq0CskGl0CJBw8U9T0B4h04NCZYnQJkPDA5Cjkldo/1W1j2H83JKdO3l+WJHvVaiJPUVZVY1ugWA5NhxI9YijhSc0VEhISEBQUhKFDh+Lxxx/H0KFDERQUhH379qldNMUp0Z7+uzIFgighZNMRBUpEpB1TNh2FIEr4aGWKIt/H/rshOXUy7vcDEEQJiacLXFAiIn3IK62CIEr424+7ZG3PUKIwhhKe1JwpJiYGffr0ga+vL0aOHIlZs2YhNDQUs2bNwsiRI+Hr64s+ffogJiZG7aIqprntaXtaHgRRwqNfxqGsqkahUhFpQ1lVDR79Mg6CKCE+Pa/Z38f+uyE5dTJhnXVhy+1pzf83IHIX2QXlEEQJf1+wV9b2mgsl77zzDjIyMq67TXp6Ot555x3XFMhODCU8qTnTkCFDEB8ff91tdu7cieeff95FJXK+5rSnimoT/jJjCwRRQszRCwqWikg7Yo5egCBKeOLrraiqad7aJXKON0+6UwvIq5O6O1ZRR3JcUCIifUi/UApBlPDukkRZ22sulCxfvhw+Pj4YMGAAAgMDsXz5cmzcuBHLly9HYGAgBgwYAB8fH6xYsULtojaKoYShhJTVnPY0q3ZNkg+WJStYIiLtGbnUOqZhwY4Tzfqe6x1vnninFpDXB82IPAZBlLDhwFkXlIhIHw6dKYYgShi9Qt45WHOhBADMZjM2bNiA9957D4888gi6deuGRx55BO+99x7Wr18Pk0n5VWyVwlDCUELKcrQ9XSipRPegSPgFRuBcUYXCpSLSlpN5ZbhnQjjunxSFi2WOD3q/3vHmiXdqAXl90Jy4dAiihFWJ2S4oEZE+JJ4ugCBKGPf7AVnbazKU6BlDCUOJK4miqHYRnM7R9iSuOQhBlDAj8pjCJSLSprpHiALXH3L4O9h/NySnThbsOAFBlPDLntPOLxCRTsSnW8d0Bqw9KGt7hhKFMZTwpOYMa9eubfBas2YNHnzwQbWL5nSOtKe0C6XoEiChz9QYlHCFZfIQxeVGPDjFOvX1aQfXy7D3eLNYLCgvL6/3cjdy6uSXPacVeXyOyJ3E1o53m7xR3syXDCUKYyhhKHGGDh06YMqUKQgJCan3uueee9QumtM50p7G/Lofgihh8a5TTigRkXbN22b9xf6z1fJ+mbya3OMtISEBDz30EFq2bIkWLVrY1hTT2tphSpBTJ6sSsyGIEubEpbugRET6sPngOQiihOkR8p5YYChRGEMJQ4kzPP7447hwoeHsUa+99poKpXEte9tT5sVL6BIgoe/nsag0anf8GZEzlFYa0SskGvdMCMdZB8ZSyT3e7r33XgQFBeGPP/5AZmZmvZe7kVMnGw6c5eOiRFdZnXwGgijhu9g0WdszlCiMoYShxBlqajx3fQ1721Pgeut6AT9svf7U4kTu6rvYNLsembiS3OOtffv2sFgsdn+/Hsmpk8jDORBECVM2HXVBiYj0YcW+TAiihHnb5D3WyFCiMIYShhJXiIiIULsILmNPeyouN8IvMAI9giNRXMGxJOSZisqr4RcYgfsnRaGi2r67hXKPt1GjRmHjxo2OFE935NTJlmPWZ+cnbTjsghIR6cPS3daxVqHxJ2Vtr+lQkp+fj08++QSDBw9Gv3796r20iqGEocQVBgwYoHYRXMae9rRsb2azZx8icgfjfj9gnaI2yb4pauUeb4WFhejatSsGDhyI4cOH13u5Gzl1svV4LgRRQtB6hhKiOkt2nXKfUPLss8/C398f8+bNw9KlS+u9tIqhhKHEFR599FG1i+Ay9rSnF3/YBUGUcCC7yIklItK+uvUBXp2/267PyT3ehg4dip49e+Kzzz5rMAGHu5FTJ9vTrFOfTlzHH0SI6iyuDSU/75Q36YymQ0nbtm1RWVmpdjHswlDCUOIKvFPSUHZBOQRRwpPfbPOYZ92JrsViseCx6VvgEyDZtZii3OOtdevWKCkpcbR4uiKnTnam59u1HgORJ/h55ym7ZsLUdCjp168fzpw5o3Yx7MJQwlDiCgwlDS2vfXRL7tSDRO5u0obDdj/CJfd469u3L86dO+do0XRFTp3szshv1lTMRO4oNP4kBFHCEncIJTNnzkTv3r2xePFihIeH13tpFUMJQ4kr8PGtht77JQmCKGH3iXwnl4hIH+p+vR+9Iln2Z+QebzNmzMBDDz2ERYsW6eb87Cg5dbLnxEUIooRPw1JdUCIifagLJUt3n5a1vaZDiY+PT6OvLl26qF20a2IoYShxhcLCQrWL4DJy2pPFYsEDk6Nwb1AEqmvMLigVkfZVVJvgOyEcj34ZJ/szcvtvPZ6fHSWnThJOWcfwfLyKoYSozoId1sVcf9lzWtb2mg4lesRQwlCiJovFgqioKLWLoSg57aluPMlLP+5yQYmI9OOFufEQRAk5xfLGZ7L/bkhOnSTVTiww9vcDLigRkT78tN0aSpbtlbeoKkOJwhhKeFJTQ2pqKj755BPceeedaNu2rdrFUZSc9hR9JIcz3xA14pOwVAiihF0Z8h5rbE7/bTKZMHv2bIc/r1Vy6mR/ViEEUcJ/V6a4oERE+jC/NpQsd5dQEhUVhTfeeAODBg3CG2+84ZJfgY1GI0aPHo327dvj1ltvxfjx42XP5sNQwlDiKufPn8fMmTPRq1cveHl5oUWLFpg2bRqKi4tdsn97jpOSkhK8/vrraNOmDe644w588803svcjpz0ttnMudCJPMScuHYIo4beELFnby+2/X3311Qbrkzz33HPo2LFjc4orm6v6H0BenaTUhpJ//8ZQQlRn3jZrKFmxzw1CycKFC3HbbbchICAAP/30EwICAnD77bdjwYIFTt3vpEmT0K9fP+Tm5iIrKwt+fn6YO3eurM8ylDCUOMOSJUvQq1cvvP766wgNDcUzzzwDb29vtG/fHu+//z527tyJO+64A7m5uS4rkz3Hydtvv41hw4ahpKQEhw4dwu23345NmzbJ2o+c9vRV5DEIooQNB87a9Xcgcnerk89AECXMij4ua3u5/ffVa5OMHTsWnTt3xsKFC5tTXNlc1f8AvFNC5KgftmZAECX8us/+H0U0F0r8/Pywf//+eu+lpKSgW7duTt1v586dsXHjRtufFy5ciN69e8v6LEMJQ4kz3HnnnSgtLUVeXh68vLxw++23Y/ny5aiqurz+gKtDidzjpLy8HDfeeCMOHrw8VebEiRPxt7/9TdZ+5LSngLUHIYgS4tPzZH0nkadYu98aSj6SebHcnP47IyMDPXv2dPjz9nBV/wPIq5PkzAK76pnIE9SFkpUO3KnVXCjp0KEDjEZjvfeMRqNTL9wLCwthMBiQmXn5VlNiYiJatmwp6xEuhhKGEmd48cUXMW7cOHz88cd44okn4O3tjY4dO2LMmDFISEgA4NpQYs9xkpKSAm9vb5jNl2fFCgsLQ9euXWXtS057+rT2ufk9Jy7K/BsQeYYZtXcRBVGStX1z+u+amhrccccdDn9eLlf2PwAHuhM56vst1sdHf090g1DyzDPPIDg4GCaTCQBgNpsxadIkPP30007bZ3Z2NgwGA4qKimzvpaenw2AwNLq6fEhICLy8vGwvg8HxamQooWsxGo3YsGEDjh49CgA4e/YsvvrqK9x///1o0aIF7r33XrRq1QrHj8t7RKO57DlO4uPj0a5du3rvxcTEoFOnTo1+tyPH1MerrKFk70mGEqIr1Q00VTqUHD16tN4rOTkZY8aMwZgxY+q97wzO7H8Ax/qgxNpQMo6hhMimbkzbqkR5C7hqOpRkZGTAz88PHTt2RK9evdCxY0f4+fkhLS3Nafus+wUmK+tyqktKSuKdEjswlLhWcnIyPvroI3Tq1Ak33XQT3nrrLafv057jpO6XyivfX716taJ3SiZvPAJBlBB79ILMvwGRZwg/dB6CKOHzzfICgtz+u25yjSsv3q9+tWjRojlFvyZX9j+AfeuUjFvFUEJUZ3ZsbShJcoNQAlhvB+/atQurVq3Crl27UFNT4/R9du7cud4guNDQUI4psQNDiTpMJhM2b96M119/3SX7k3uc1D3TfejQ5el6AwMDFR1TYuv4ZP4aQ+QpltTOTLdgxwlZ2+ul/3ZV/wPIq5O9Jy9y8USiq3wXmwZBlLA6+Yys7TUXSqKjo23/Pzw8/JovZwoODkb//v2Rl5eH7OxsdO/enbNv2UEvJzVqHnuOkxEjRuDFF19EaWkpDh8+jE6dOik6+9bviVnWGYZinHcXlUiPvoz4A4IoYWPqOVnb66X/dlX/A8irkz0nrKHk0zCGEqI638ZYQ8kavYaSK2fu8PHxafTVpUsXp5bBaDRi1KhRaNeuHTp06IBPP/2U65TYQS8nNT165513kJGRcd1t0tPT8c477zi9LNc7ToYMGYJp06bZti0pKcFrr72GNm3aoFOnToqvU5KcaZ2Oc/SKZPv/IkRu7F9LEiGIElKzi5reGNc/3jyx/wHk9UG7T+RDECWMX81QQlSn7k5JmLs8vqU3DCUMJc60fPly+Pj4YMCAAQgMDMTy5cuxceNGLF++HIGBgRgwYAB8fHywYsUKtYuqGDntqaTSCEGUMOjrrS4oEZF+9PsiFvdMCEel0SRr++sdb57Y/wAyQ0mGNZR8tvpgk9sSeYq57jTQfdy4cY2+/8knn7i4JPIxlDCUOJvZbMaGDRvw3nvv4ZFHHkG3bt3wyCOP4L333sP69etts9W5C7nt6Zlvt0MQJeSWNJwlj8gTZReUQxAl/HVOvOzPNHW8eVr/A8jrg3amW0OJuIahhKhO3Tolv7nDOiVt27Zt9P0OHTq4uCTyMZQwlJCy5Lan4A2HIYgS1qdwVXciAFi2NxOCKGF6xDHZn2H/3ZCcOolPz4MgSghYy1BCVGfeNuuU5Cv2ZTa9MTQaSuoGs7dq1QoRERH1BrjPnTsXPj4+ahfxmhhKeFIjZcltT3F/XIAgShi1nONKiIDL40n22bF+D/vvhuTUyY40ayiZsO5Qk9sSeYoFO6yhZNme07K212QoqRvQ3qJFi3oD3H19fTFgwABs3rxZ7SJeE0MJT2qkLLntqarGhPsnR8EvMAKXqpw/dTiRlhWXG9EtMAIPTI6C0WRu+gO12H83JKdOtteGkokMJUQ2ofEnIYgSlu4+LWt7TYaSOv/4xz/ULoLdGEp4UiNl2dOe6lZ2X7tf3vSDRO5q2Z7TEEQJgevtu0hm/92QnDrZejzXofomcmeLdlrXSVq085Ss7TUXSrSwTklzMJTwpKaWoiJ5U37qjT3tqW6tgL/9uMuJJSLSNovFghfmxkMQJRw8Y1+/4Gj/7a79DyAzlByzhpKg9YddUCIifahbvDU0/qSs7TUXSrSwTklzMJQwlKhl8ODBmDNnjtvNfmNPe7JYLHj22x0OXYwRuYt9tauLD5kdL3uNrTqO9t/u2v8A8uok9qh1TNvkjUdcUCIifai7Y7tgxwlZ22sulOgdQwlDiVosFgtCQ0MxaNAgbNy4Ue3iKMbe9rRin3XGoTG/7ndSiYi07e1FCRBECRsO2D8TnaP9t7v2P4C8OpEOnocgSvgy/A8XlIhIH1YmZEEQJfyw9fqLrtZhKFEYQwlDidqKi4vx0ksv4cknn0Rqqv5XF7a3PVUaTej3RSx8AiQcyylxUqmItOnw2WIIooSBM7aixo4B7nWa23+7W/8DyKuTtfvPQBAlzIpJc0GJiPRhfcpZCKKEb6KPy9pe06Fk1apVOHbMOr/68ePHMXDgQAwePBhpado96BlKGErU8tFHH2HQoEHo378//vWvf2HBggX45z//ia+++krtojWLI+1p6e7TnB6YPNLIpdZpgJfvlbcuwNUc7b/dtf8B5NXJr/vs+0WYyBNEHs6BIEr4fPNRWdtrOpTcc889uHDhAgBg6NCh+PDDD/HZZ5/hqaeeUrlk18ZQwlDiKpmZ9S864uLiUFLS8M7Aiy++6KoiOYUj7anSaMKjX8ZxbAl5lN0nrKuKP/H1VlTX2H+XBJB3vJlMJvzzn/+E0Wi0veeu/Q8gr04W1w7o/VnmLENEnmCbnbPSaTqU1K3oXllZiVtuuQUVFRWoqanhiu5OxFCiHw8//DCqqqoavJ+Tk1Pvz+fPn3dVkZzC0fb0e6L1l8uX5+22e7Avkd6YzZdn3Ao/5PgxL/d469Spk6xB7XrvfwB5dVK3crWjd6iI3NHe2kk3Pl4l71FOTYeSu+++GxkZGdiwYQMGDRoEAKiqqsItt9yibsGug6GEocRVAgMD8e6779Z7LzU1Ff/7v/+rUomcw9H2ZLriIs2RAb9EerIm+YwiIVzu8fbJJ59g+fLlDu9HT+TUyXexaRBECWFJ2S4oEZE+pGYX2TXxjKZDycyZM9GqVSu0atUKYWFhAKy3iPv3769yya6NoYShxFUsFguefvppLFiwAACwefNm3H777Vi5cqXKJVNWK9aa5AAAIABJREFUc9pT4ukCCKKE/tPiUF7NVd7JPRWVV+OhqTEQRAn7swqb9V1yj7eHH34YN998M/7zn/9g586dqKioaNZ+tUxOnUyPOAZBlLAp9ZwLSkSkD8dzSiGIEkYuTZS1vaZDCQCkpaXh5MmT9f586JB2V0xlKGEocaaoqKh6j0Pk5+fD19cXY8eOxZ///Gfs3btXxdI5R3Pb079/S7FO1RnBqTrJPX22+iAEUULA2uafG+Ueb2vWrEFgYCD++te/4s9//jNuuOEG9OjRA2+99Vazy6A1cupk8sYjEEQJ0UdymtyWyFNkXSyHIEp4I1TetYnmQ0lUVBTeeOMNDBo0CG+88QYiIyPVLtJ1MZQwlDjTAw88gJYtW+L222/Hk08+ibFjx2LkyJG48847kZ6ernbxnKK57elCSSXunxQF3wnhOHKuWKFSEWlD3UKJfT+PRXG5sekPNMHR4y0vLw/R0dFuMdvW1eTUScBaazDcmZ7vghIR6UNuSaXtsVI5NB1KFi5ciNtuuw0BAQH46aefEBAQgNtvv932uIoWMZQwlDhbdXU1Dhw4gKVLl2LcuHF46qmncNttt+HGG29Er1698Pbbb6tdREUp0Z6W77UuqPh/3++EycxB7+QeqmpMGPzNNkUfG2L/3ZCcOvlopfWObHJmgQtKRKQPJZVGCKKE52fHy9pe06HEz88P+/fXHxyTkpKCbt26qVSipjGU8KSmlnPnziEiIgLTp09XuyiKUqI9mc0WvDxvNwRRwiJO2Uluom5w9TuLExSbYY79d0Ny6uT9ZUkQRIl3Y4muYDSZIYgS/Gduk7W9pkNJhw4d6s2DDgBGo1H1C/frYSjhSY2UpVR7Sr9Qiq4Tw3FfcCTOFrnvoFzyDBm5Zeg2MQLdgyKRXVCu2Pey/25ITp2MWJQAQZRwMq/MBSUi0o9uEyPQ74tYWdtqOpQ888wzCA4Ots2FbjabMWnSJDz99NMql+zaGEp4UiNlKdmeZsVYf1l+d0ki1y4h3TKbLRj+0x4IooSFO042/QE7sP9uSE6dDJ9v/fc4X8wfPIiu1CskGvdPipK1raZDSUZGBvz8/NCxY0f06tULHTt2hJ+fH9LS0tQu2jUxlPCkRspSsj1VGi8/gy8d1P+ibuSZ6hYG/euceNSYHFu5/VrYfzckp06Gzt0JQZRQeKnaBSUi0o9Hv4yD74RwWT8EajqUAEBNTQ127tyJVatWYdeuXaip0fZaAwwlPKmRspRuT3UrzD78RSyKK5o/WxGRK+WXVaFXSDS6BEg4eKZI8e9n/92QnDp5atZ2CKKESmPTq9wTeZK6HwKrapo+NjQfSgCgsrISp06dQmVlpdpFaRJDCU9qpCxntKfxq1MhiBImrNPumkdEjflv7SxPIZuOOOX72X83JKdOHpu+BYIo8bFQoqu8MDcegiihqLzpu4iaDiUFBQV4+eWX4e3tDS8vL3h7e+OVV15Bfr525wFnKOFJjZTljPZ05QrYSac5hSfpw/a0PAiihEe/jENZlXOeGtBD/71u3Trcc889+NOf/gR/f3+cOnXtGfVmzpyJBx54AG3btsVdd92F8ePHN5hApyly6uShqTHoHqTtddSI1FA33krOBDOaDiXDhg3DK6+8grS0NFRVVSEtLQ3Dhw/HsGHD1C7aNTGU6OOkRvrhrPa0PuUsBFHC07O2o7pG2efyiZRWUW3CX2Zscfqq4Vrvv48fP442bdogOjoaFRUVGDt2LPr27XvN7b/66iskJibCaDTizJkz6Nu3L4KDg+3ap5w6uS84En2mxtj1vUSeoG5muhMyZqbTdChp3749ysvrT3V46dIltGvXTqUSNY2hRPsnNdIXZ7Uni8WCt37eB0GUMDcu3Sn7IFLKV5HHIIgS3l+W5NT9aL3/DgoKwksvvWT7c2lpKW666SakpqbK+vycOXMwePBgu/bZVJ1YLBb4BEh4bPoWu76XyBP8c7E1lGTklja5raZDSe/evRvclj158iR69+6tUomaxlCi/ZMa6Ysz21PmxUvwC4xAt8AIri9AmnUspwT3TAhHz0lRyCl27thKrfffw4YNQ0hISL33evTogRUrVsj6/Msvv4x///vfdu2zqTqpqjFBECU8+c02u76XyBP8a0kiBFFC2gWdh5Lp06fDz88Pc+bMwbp16zB79mx0794d06dPR3h4uO2lJQwl2j+pkb44uz3N23YCgijh7wv2cpAqaY7ZbMFLP+6CIEpYsuvaYyeUovX++8knn8R3331X773HHnsM8+fPb/KzCxYsQKdOnXD+/PWnAw8JCYGXl5ftZTBc//KovLoGgijhmW+3N/0XIPIwI5cmQRAlHMspaXJbTYcSHx+fJl9dunRRu5j1MJRo/6RG+uLs9mQ0mfHcdzsgiBI2HDjr1H0R2Wt18hkIooT/+34nTGbnh2Y1++9BgwbBYDBc8wVY75RMmTKl3ud69uzZ5J2SFStW4H/+539w8OBBu8vVVJ1U15ghiBIG804JUQP/7xdrKDl6TuehRI8YShhKSFmuaE8JpwogiBL6T4vDJSfNakRkr9JKI/p+HgtBlJCSVeiSfWq9/w4KCsLLL79s+3NZWRluvvnm644pWbFiBW677TakpKQ4tM+m6sRstkAQJTzx9VaHvp/IndU9vnU8R+ePb+kRQ4n2T2qkL65qT3XrP8yIPOaS/RE1ZVr4HxBECZ+EyRvErQSt9991s2/FxsaisrISH3/88XVn3/rtt99w6623IinJ8QkC5NRJFw50J2rU3xfshSBKyC4ob3JbhhKFMZRo/6RG+uKq9pRTXIn7giPRbWIETudfcsk+ia7lRF6ZbXB7bqnrFg7WQ/+9du1a+Pr64uabb8agQYPqTYizYsUK9OjRw/ZnHx8f3HDDDWjdurXtdeV/l0NOnfQIjkTvKdF2fS+RJxj2/U4IooSLZVVNbstQojCGEn2c1Eg/XNme6ga9v7sk0WX7JLqaxWLB27Vz+y/YccKl+2b/3ZDcFd27BEgwu2DcD5GePD1rOwRRQkW1qcltGUoUxlDCkxopy5XtqarGBP+Z2yCIErYcu+Cy/RJdKe6PC9aB0zO3uXxhT/bfDcmpk+dnx0MQJRRX2LdaPJG76/t5LLoFRsia3VLToWTgwIFYvHgxysr0s34AQwlPaqQsV7enrcdzIYgSBn29FVU1Tf+yQ6SkqhoTnvh6KwRRwrbjuS7fP/vvhuTUyZuh+2SvWk3kKepmppM7CYSmQ8n8+fMxYMAAtGnTBiNGjEBcXJzaRWoSQwlPaqQsNdrTyKXW2UJ+2Jrh8n2TZ/txWwYEUcLIpeo8Qsj+uyE5dRKw9iAEUcJWFYIkkVZlF5RDECW89tMeWdtrOpTUOX78OCZMmIC7774bgiAgODgYGRnavFhgKOFJjZSlRnvKvHgJ3SZGoEdwJHJLXDfImDybFiZbYP/dkJw6mb/9hMsWuCTSi90Z+RBECWN/PyBre12EkjoJCQno3bs3vLy80KpVKzzzzDM4cuSI2sWqh6GEJzVSllrtaXrEMQiihPGrXTcdK3m2jzQwLTX774bk1EnEofMQRAmTN2rrmoRITaHxJyGIEn7aLm/CDs2HkvPnz2PGjBno2bMnOnfujAkTJiA9PR1lZWWYOnUq7rnnHrWLWA9DCU9qpCy12pN14boY+ARIOHy2WJUykOdIOm1dwPORabGqLuDJ/rshOXWSefESBFHCy/N2u6BERPrw8apUCKKEHWl5srbXdCh55pln0KpVKwwfPhzh4eEwm+vPQmKxWNCmTRuVStc4hhKe1DzBunXrcM899+BPf/oT/P39660TcLUlS5agRYsW9dYJ+Pnnn2XvS8329FtCFgRRwvCf9siaOYTIESazBS/Mtc7etOHAWVXLwv67ITl1YrFY8OCUaPgFRsBocu2MaURa9dx3OyCIEvJKm16jBNB4KJkzZw4uXrx43W2uDipqYyjhSc3d1a2oHB0djYqKCowdO/a6KyovWbLkuv+9KWq2J5PZYutUIw+fV60c5N7qwu8r83arHn7Zfzckt07eWWxdWyY5s8DJJSLSvrzSKgiiBP+Z22R/RtOhRI8YSnhSc3dBQUF46aWXbH8uLS3FTTfdhNTUxsde6DmUAMDuE9aBen+ZsYVTBJPitPaYoNrHmxbJrZPFu05BECV8E33cySUi0r4NB85CECUErT8s+zMMJQpjKOFJzd0NGzYMISEh9d7r0aMHVqxY0ej2S5YswZ/+9Cfcdttt8PX1xccff4xLl+TPLKSF9vTeL0kQRAnzZQ7WI5JLaxMqaOF40xq5dXIq3zqu5PnZ8U4uEZH2jVmxH4IoIfpIjuzPMJQojKGEJzV39+STT+K7776r995jjz2G+fPnN7r9yZMnceLECZjNZmRkZODxxx/Hv/71r2t+f0hICLy8vGwvg0H9rulU/iV0nRiOnpOiZD8bS9SU7IJydJsYgfs0NPU0+++G7KmTZ7+1Pu55PKfUiSUi0rbSSiP8AiNw/+QoVBrlP2HAUKIwhhKe1PRs0KBBMBgM13wB1jslU6ZMqfe5nj17XvNOydX27NmDtm3byn52Xivt6Qvp6P9v797joioT/4GPmgaColLRljVy0a+iKcWaphleupDbUuqWtav1W/1umX6/XvjqooKilWZl2mZhXlq8sGm6QMJwE/FKomKoIDoCiuAVBLkzMMzM5/cHcpZR0JnhzMyZ4fN+veaPOZw555mHc/vMOc/zQB6kwMLIM9YuCtmJmf9q/CXx23051i6KQCr7m5QYUycbDzV2gfqZItuMJSKSth132sn9fbdx50uGEpExlPCkZu9CQkIwceJE4X1VVRUcHBxabVNyt2PHjsHZ2dnmQkl5rRrPfrIX7gsVyL5WYe3ikI1r6gL4hZX7jPol0dyksr9JiTF1UlxZB89FcRi8LAlVVuzamchadDodXv76IORBCvxWcNuozzKUiIyhhCc1e9fU+1ZycjJUKhUCAwPv25A9Pj4e16839lyVn5+PUaNGYcqUKQavT0rb0/a0y5AHKfDuhjSr95JEtkur1SHgu1TIgxSIzrBuF8B3k9L+JhXG1sn8XY1jM2w8dNFMJSKSrv3ni4TeBI3FUCIyhhKe1NqDyMhIeHh4wMHBAX5+fnrjlERERMDb21t4P3/+fLi5ucHR0RFPPfUUZs+ejcpKw5+3ltL21KDR4pU1B41uvEfUXHRGY680Ad+lQquVVriV0v4mFcbWSV5xFfosVMD302TeLaF2RaPV4fVvDpt8jmQoERlDCU9qJC6pbU+Hc4ohD1LgpS/3s4tgMppKrcHwlfskO56F1PY3KTClTppGsl4Zd84MJSKSpp/TCyEPUuBP600bc4mhRGQMJTypkbikuD1NCz8BeZAC4amtj2RP1JINh/IgD1JgZsRv1i5Ki6S4v1mbKXVSVKnCwKWJ8FwUhws32RMX2b9bVXV47pO9kAcpcLqwzKRlMJSIjKGEJzUSlxS3J+WNyjuPZ+xFNR/PIANVqNQYsjwJHovicOmW4WP1WJIU9zdrM7VOfjzSOJii/zeHJdWZAZHYdDqdMJ6XMYMl3o2hRGQMJTypkbikuj3N+/mU5LpzJWlbnaSUfLfSUt3frMnUOtFqdfjzpjTIgxQI3XNW5FIRScfWo/mQByng9+V+1NSb/kMdQ4nIGEp4UiNxSXV7KiytgdfiOAxamojb1fXWLg5JXHFlHQYsSUDf4HhcL6+1dnFaJdX9zZraUic3K1TwWZ4EeZACPx0vELFURNLwa+4teCyKQ9/geJMf22rCUCIyhhKe1EhcUt6elv6SBXmQAivYmJUeIHTPWciDFPg0VtqD6kl5f7OWttbJsYsl8FocB49FcTh4oVikUhFZ37nrFRi8rDF0R2VcafPyGEpExlDCkxqJS8rbU9Ov3/0k/us3WdeV2zXouzgeA5cmolTid9WkvL9Zixh1EpVxBfIgBfoFxyM195YIpSKyrtyiKqFh+1eJSlGWyVDSzJYtWzB06FB0794djz/+OP77v/8b5eXlRi2DoYQnNRKX1LenrxKl306ArOv/7gymtzb5grWL8kBS39+sQaw6aXruvl9wPPYri0RZJpE1ZF0th++nyZAHKbAs5qxogwkzlDQTFhaGAwcOQKVSobS0FK+//jqmTp1q1DIYSnhSI3FJfXtq3qPSxeIqaxeHJCbnZiXcFyrw7Cd7bWIgPanvbwAQFRUFT09PODo6YvTo0XqDt7ZGp9PBz88PMpkMVVXG7adi1klTMHFfqMC2o/miLZfIUg5dKIb3kgShAwexAgnAUHJfe/bsgbu7u1GfYSixjZMa2Q5b2J5+ONg49sSsf0lz7AmynhnbT0IepMCmwxetXRSDSH1/UyqVcHZ2RlJSEmprazF37lz4+vo+8HObNm3CSy+9ZPVQAgAxp6+hb3C8cFFX36AVdflE5qDT6RB2IA/uCxWQBynw/YFcUQMJwFByX4GBgXjjjTeM+gxDifRPamRbbGF7Uqk1eH5F463srKvGPfJJ9uvMlTLIgxQYvnKfzYxTIfX9LSQkBBMmTBDeV1ZW4uGHH8bp06db/czNmzfh6emJEydOSCKUAMDJy7eF5/EDvktFQUmN6OsgEkt5rVoYh6R/SAL2nL5mlvUwlLQiPj4e3bt3R2Zm5n3nW7ZsGTp06CC8ZDLTq5GhhOhetrI9bU+7DHmQAh/887i1i0ISMfXH45AHKbDDhrqClfr+FhAQgGXLlulN8/b2RkRERKufmTx5MjZs2ID8/HzJhBIAuFpWi0lhv0IepMCgpYmIyrgi+i/PRG118EIxhq3YB3mQAmO+OoALNyvNtq52E0qaniVt7dVcSkoKevXqhX379hm9Ht4pkf5JjWyLrWxPao0WL325H/IgBY5fKrV2ccjK0i6WQB6kwOivDqBBYzuP50h9fxs7dizWrl2rN23EiBFYv359i/PHx8djxIgR0Ol0BocSMX9sfJAGjRZfJynR584jMe//eByFpbxrQtZXXqvGwsgzkAc1bpvzfj6FSpXarOtsN6HEUCkpKejZsycSExNN+jxDifRPamRbbGl7+uXUVciDFJgU9it/8WzHdDqd8At4jJkeczAXa+5vhvx4GBAQgOXLl+t9buDAgS3eKampqYGXl5fwxIPU7pQ0l55finFfHxQej9l46CLUNhRmyX5otTrsSi8UHi/0/TQZe7NvWmTdDCXNHDhwAD169IBCoTB5GQwltnURSdJnS9uTVqvDa2sPQR6kQMp5yxzESXr2ny+CPEgB/28OQ6u1rXAq9f0tJCQEEydOFN5XVVXBwcGhxTYlp06dwkMPPQQ3Nze4ubnhkUcegUwmw2OPPYbY2FiD12mpOqlr0GBt8gV4LY4THpVJzr7JHzjIYjKvlGPinR9U+ixUIOjfZ3DbgmMrMZQ0M3r0aHTs2BFOTk56L2MwlEj/pEa2xda2p33nbtrsBSm1XYNGa9PBVOr7W1PvW8nJyVCpVAgMDGy1962GhgbcuHFDeDU1dL948SLq6uoMXqel6yS3qFJojyQPUuDPm9Jw7nqFRctA7cvlkmr8z08Zwjb3x3VHcKqwzOLlYCgRGUOJ9E9qZFtsbXvS6XTCL03RGVetXRyysKYODyZvOGqTv3Dbwv4WGRkJDw8PODg4wM/PT2+ckoiICHh7e7f4OSk/vtWSA8oi4ZGuPgsVmLvzFC6XVFulLGSfblXVIXTPWeHu3PMrkrHjeAE0VvpBjaFEZAwltnFSI9thi9tTen4p5EEKDP0s2SYGzCNxlNeq8ewne+G+UIHsa7b5y7Yt7m/mZs06adBose1oPnw/bXy+33NRHBZFZeJ6ea3VykS2r1KlxtrkC8IgiIOWJuK7/bmorbdu1+UMJSJjKOFJjcRlq9vTvJ2nIA9S4DNFtrWLQhbySWw25EEKLIy8f1fyUmar+5s5SaFOauob8N3+XDwTmgh5kAJ9g+PxSWw2SqoMfwyNSKXWYMOhPPgsT2rcjhbH49PYbIu2G7kfhhKRMZRI4wBO9sNWt6eiShUGLU2Ex6I4s/brTtKQeaUcHoviMCg0Ebds+ELRVvc3c5JSnZTXqrE6SYkBd37hHrAkAauTlKgwc1etZNvqGxrvuA39rHGQX49FcViw+zSu3JZW99MMJSJjKJHWAZxsny1vT/9MvQR5kAJvfZ9qU2NVkHHqG/7TuN2WBkpsiS3vb+YixTopqarDp7HZ6BscD3mQAoOXJeGHg3lQqa37+A1JS4NGi5/TCzHi8xShEfv//JSBi8XGtauyFIYSkTGUSPMATrbLlrenBo0Wb32fCnmQAutScqxdHDKTf+zLgTxIgfc2ptlk4/bmbHl/Mxcp18n18losjMyEx6L/NFT+17ECjnHSzmm1OsScvoYxqw8IYWT6lnTJ9+LGUCIyhhJpH8DJ9tj69nTpVjX6hyTAc1Ecsq6WW7s4JLL0/FJ4LIpD/5AEFJRI61EIU9j6/mYOtlAnF4urMOtfvwkXoKO/OoCY09fYLXk7o9PpcEBZBP9vDgvbwl82HUNGwW1rF80gDCUiYyixjQM42Q572J6auol96cv9KKuRRoNCaruSqjoMW7EP8iAFdqUXWrs4orCH/U1stlQnWVfL8X6zMU7G/+MwDiiLbP4OHj3Y2Wvl+MumY8L/fmLYrziaV2LtYhmFoURkDCW2dQAn6bOH7Umn02FmROOvmFM2H2P7EjtQ36DFuxvSIA9S4P923TuauK2yh/1NbLZYJ2kXSzDhzqOjTY8Wnr3GO7X26Hp5Lf5v12n0Wdj4vx67+gCSs2/aZBBlKBEZQ4ltHsBJuuxle6qpbxAaQwdHZ9rkCYMa6XQ6zNmRIfwSXVNvP2PR2Mv+JiZbrROdToe92TfxcrMBGP+++wyKKlTWLhqJoKquAV8lKvFfIY2dHfh+uhfb0y7b9I9eDCUiYyix3QM4SZM9bU+FpTXw/bSxS8YVcecYTGyQTqfDp3fGIxm5KgVFlfZ1gWdP+5tYbL1OGjRabE+7jGc/2St0I7wuJYc9ddmopv9n04Ca/xUSj9VJSrsYqJehRGQMJbZ/ACdpsbftSXmjUhi46vP48wwmNkSr1WHpL1mQByngszwJuUXS7FazLextfxODvdRJhUqNlfHn0Hdx4y/rL6zch19OXeUxyIYczinWu/M1f9dp3Ci3nx9GGEpExlBiPwdwkgZ73J6yrpZj8LLGYDJnRwbqGviLpdTV1msw807vRr6fJkN5wz4HxLTH/a2t7K1OCkpq8HHESaG9yZvfpeI3G+mdqb3Kv1WN6VvS9XrUyr4m7e59TcFQIjKGEvs7gJN12ev2lFtUiZGrGge0Clh3xC66k7VX+beqMf4fjV1sjvpiP/IkOvCYGOx1f2sLe62TE/ml+OO6I8KF7rydp3CT7U0kpfLO3S2vxXFCD457bbQRuyEYSkTGUGK/B3CyDnvenooqVfjT+l8hD1Jg0NJE7DhewHEFJESr1WFb2mX0D0kQejC6XW3fXTrb8/5mKnuuE61Wh8jfrmDoZ41t3byXJCDsQB7v3lqZVqvDz+mFQhtE7yUJWH/Q/v8vDCUiYyix7wM4WZ69b08NGi2+TlIK3Tm+9X0qTl7moxTWduxiCd74tvFXZM9FcViXkmPTvdoYyt73N1O0hzqpqmvA5/HnhfYmfl/ux75z9vuLvJSdvFwqHHvkQY3tRuytQ43WMJSIjKGkfRzAyXLay/Z05koZApo9SjFl8zEczinmnRML0mp12H++SBh/pCkkZl1tP+M7tJf9zRjtqU4u3arGtPATwvb/wT+P2/XjilJyvbwWs+90NS4PUmDC96k4XVhm7WJZFEOJyBhK2tcBnMyvPW1PGq0Ou09eweivDggnppGrUvD13gvIvFLOXy3NQKvVIetqOVYlnMeIz1OEeh/91QHsOX2t3dV5e9rfDNUe62S/sghj7hyHPBfFYUXcOVSq1NYull1SqTX4x74c4THRYSvab69oDCUiYyhpnwdwMp/2uD01aLSIOX0N7238zy/28iAFfv9ZMmZsP4n1B/Pwa+4tXC2rhUakOylarQ619RqUVtfjyu0a5BZVIvNKOY5fKsXBC8VIOX8Tydk3kXj2BhKyriMu8zpiTl/DL6euIirjCv59svEVlXEF0RlX8cupq9hz+hpiTl9D7JlriMu8jvjM60jIuoHEszeQdPYGkrNvYt+5m9h/vgiHc4pxIr8UmVfKkXOzEoWlNSiurEOlSg21SI9N6XQ6XC2rxX5lETYcysOH29KFXtCauth8d0Maks7eaLd3qNrj/vYg7bVO6hu02HjoIgYuTRR6nduVXthu9w2x6XQ6xGVeF34M6Rccj6+TlHY1GKuxGEpExlDSfg/gZB7tfXsqKKnBhkN5ePuHo/BYFKcXUuRBCngtjsOLX6TgD98exnsb0/DhtnTM3pGBOXdec3eewv/+lIG/bU3HlM3H8Kf1v2L8Pw5jzOoDeGHlPvgsTxJGBJbyy3NRHAYtTYTvp8l48YsUvLLmIP647gje/uEopv54HB9uS8ecHRkI/Pk0Zu/IwMyI3/C3remY+uNxjP/HYTy/IhmeLdRfv+B4/HlTGjYcysO1slpr/7utrr3vby1p73VSVKnC/F2nhX0m4LtUZLAL4TYpLK3B+z8eF+p0ZsRvKCxlD4wMJSJjKOEBnMTF7ek/auobcCK/FJsOX8TcnacwKexXodccU17uCxt7dWl+oR+w7gje+eEoPvjncczYfhLzdp7CwshMLPklC6F7zmJ5TDY+U2RjZdw5rEo4jy8Tz+PrJCXW7L2ANXsv4Ou9F7A6SYmvEpX4MvE8ViWcx+fx57Ey/hxWxJ3DZ4psfBKbjeUx2VgWcxahe85iyS9ZWBiZiXk7T+HjiJP4f/88jnc3pOGt71Px2tpDGP3VAQxfuQ9D2hCgPBfF4fkVyXjr+1T8ffcZbDp8EUfzSjiq9V24v92LddLoVGEZ3vxHwpE+AAAgAElEQVQuVdinAn8+jSJ2IWwUtUaLsAN5wnFs7OoDOJpXYu1iSQZDicgYSngAJ3Fxe3owtUaL29X1KCipQdbVcpy8XHrX6zayr1Ug/1Y1blaoUF6rRn2D1mafWW561Ox2dT2uldUir7gKZ681fu+0iyU4ebkUZ66UIftaBXKLqnC7up6PnBiI+9u9WCf/ob3T7u33d34MGbg0EeGpl0R7jNSenbx8G6+uOQR5kAJ9g+Px7b4cu+/i11gMJSJjKOEBnMTF7YnIcri/3Yt1cq9KlRor484Jj0T+cd2RdtVLnTHKa9RYFJUp3GH686Y0XLpVbe1iSRJDicgYSngAJ3FxeyKyHO5v92KdtE55oxITw34VHgf9JDYb1XXtt6F2czqdDntOXxMGQHzuk72IzmifvWoZiqFEZAwlPICTuLg9EVkO97d7sU7uT6vV4V/HCvBMaGMvXS+s3IekszesXSyrulxSjSmbjwl3R4L+fQZlNfXWLpbkMZSIjKGEB3B7l5WVhVdffRWurq6QyWSoqnrwwFphYWF48skn4eTkhDfffBMlJYY37OP2RGQ5trC/RUVFwdPTE46Ojhg9ejQuXbp03/kzMzMxbtw4ODs7w9XVFYGBgUatzxbqRAqKK+v0Bv/729Z0XC9vXz3a1Tdo8d3+XPQLbmzIPu7rgzh+qdTaxbIZDCUiYyjhAdzeKZVKbN68GbGxsQaFkpSUFPTq1Qu//fYbKisr8fbbb2PixIkGr4/bE5HlSH1/UyqVcHZ2RlJSEmprazF37lz4+vq2On9paSkeffRR/Pjjj1CpVKitrcXp06eNWqfU60RqDl0oxktf7hcawm/5Nb9dNIRPzy/FK2sOCg3Z16XkoL5BnDGW2guGEpExlPAA3l7k5+cbFEqmTJmCefPmCe9zcnLQqVMnlJYa9usRtyciy5H6/hYSEoIJEyYI7ysrK/Hwww+3GjSCg4Px3nvvtWmdUq8TKVKpNfgi4bzQEP7N71Jx7nqFtYtlFuU1aiyMPCPcIZqy+Rjy2ZDdJAwlImMo4QG8vTA0lAwePBhbtmzRm9a1a1ekpqYatB5uT0SWI/X9LSAgAMuWLdOb5u3tjYiIiBbnHzt2LGbPno1Ro0bB1dUVY8aMQWZmplHrlHqdSNn5GxXC2Caei+KwKuG83YwNpNPp8Mupq/D9dO+dEe/34pdTbMjeFgwlImMo4QG8vTA0lHh4eCA6Olpv2hNPPIGEhIQW51+2bBk6dOggvGQymd77+72MmVeqL3v4Dvwe0noZuw9J2dixY7F27Vq9aSNGjMD69etbnL9v377o1q0b0tLSUF9fjxUrVkAul0Olan3Qv5aOQWQ6jVaHrUfzMXBpY0P4l77cjyM5t6xdrDbJv6XfkH1hZCbKa9TWLpbN69CBoURUzSvUWAwlZG1+fn6QyWStvpoz5k7J1q1b9aY5OTkZfKfEGPaw7dnDdwD4PaTEVr6DIcefgIAALF++XO9zAwcObPVOyZAhQzBlyhThvU6nQ7du3XD8+HGDy2Ur9Sd1N8pV+HBbunAhP2/nKZRU1Vm7WEZRqTX4OkmJvosbG7K/suYg0vPZkF0sDCUiYyjhAby9MKZNSfPebnJzc41qU2IMe9j27OE7APweUmIP36FJSEiIXkcZVVVVcHBwaLVNydSpUzF16lThvU6nQ/fu3XHs2DGD12lP9ScFiWdvYNiKfZAHKeCzPAm70guhtYGG8MnZNzFyVQrkQQr0D0nA9wdy2ZBdZAwlImMo4QHc3ul0OqhUKiiVSshkMpSUlEClUrX6HG1KSgpcXV2RkZGBqqoqTJ482ajet4xhD9uePXwHgN9DSuzhOzRp6n0rOTkZKpUKgYGB9+196+DBg+jevTvS09PR0NCAL7744oGPb93NnupPKipVaiz9JQt9FjbeNXnj2yNIu2h4V/GWlFtUiWnhJ4Q7PDO2n8TVsvbV1bGlMJSIjKGEB3B713SH5O5Xfn4+ACAiIgLe3t56nwkLC8MTTzyBrl27Gj1OiTHubgBri+zhOwD8HlJiD9+hucjISHh4eMDBwQF+fn5645S0dPzZvHkznn76aXTv3h1+fn5s6C4hGQW38db3qcIF//Qt6cgtevDYV5ZwvbwWf999Bu53gpPfl/txQFlk7WLZNYYSkTGU8ABORET2g+c089LpdIg9cw0vftH4aJTHojjM2ZGBCzcrrVKea2W1WB6TLQyA6LM8CZuPXEJdg330GiZlDCUiYyjhAZyIiOwHz2mWUdegwabDF/HsJ3v1RoU/knPLIm1Osq6WY97Pp4SxVfqHJOCrRCUqVOxVy1IYSkTGUMIDOBER2Q+e0yyrtl6DH49cEhrDN3UjHHYgD4WlNaKuq7S6HuGpl/D6N4eFdQ1ZnoSv916wuZ7B7AFDicgYSngAJyIi+8FzmnXUNWjwy6mreOeHo0JgkAcp8Po3h/FVohKHLhSjqq7BqGWq1BqcyC/Fd/tzMSnsV6G9iDxIAf9vDmPLr/moNnKZJB6GEpExlPAATkRE9oPnNOvLK67Cl4nnMe7rg3oBxX1hYwP0v4afQOies1ibfAGhe84i4LtUzN15CutScvBJbDb+e2s6Xl1zSHg0q+n1wsp9CN1zFllXy639FQkMJaJjKOEBnCwjKysLr776KlxdXQ0aLwVo7AXsySefhJOTk1l7ATNWVFQUPD094ejoiNGjR+v1JnS38PBwdOzYEU5OTsJr8+bNFixtI7VajY8//hg9evRAr169sGDBgla7ha6oqMDkyZPh7OyMxx9/HKtXr7ZwaVtnzPfw8/NDly5d9Oq+rs76j3isW7cOvr6+6NKlCyZNmnTfeaX8v5AqntOkJa+4CtvSLuN/f8oQxg0x9DX0s2R8uC0d6w/mIftaRav7OlkHQ4nIGEp4ACfLUCqV2Lx5M2JjYw0KJSkpKejVqxd+++03VFZW4u233zbbeCnGaBp3ISkpCbW1tZg7d+59x10IDw+/798tZenSpRg6dCiKiopQUFCAfv364dtvv21x3vfffx8BAQGoqKhAZmYmHn30UcTExFi4xC0z5nv4+flh3bp1Fi7hg0VGRiI6OhqzZs16YCiR8v9CqnhOk7aa+gacvVaOpLM3sCu9EH/ffQYDliTgo20nseN4AeIyryPzSjnKa9hgXeoYSkTGUMIDOFmWMSPLz5s3T3ifk5NjtpHljRESEoIJEyYI7ysrK/Hwww+3OkK1VEJJ7969sWfPHuH9xo0b4ePjc898NTU16NKlC86cOSNMW7x4Md566y2LlPNBDP0egHRDSZPQ0ND7hhKp/y+kiuc0IstgKBEZQwkP4GRZhoaSwYMHY8uWLXrTunbtitTUVHMW74ECAgLuGdzO29sbERERLc4fHh4OR0dHPPLII/Dw8EBgYCCqq6stUVTB7du3IZPJcPnyZWHaiRMn0Llz53seh8jIyECnTp2g1WqFabt27YKXl5fFytsaY74H0BhKXF1d0atXL/j6+iIyMtKSxX2gB4USKf8vpIznNCLLYCgRGUMJD+BkWYaGEg8PD0RHR+tNe+KJJ5CQkGDO4j3Q2LFjsXbtWr1pI0aMwPr161uc/+LFi8jLy4NWq0Vubi5GjhyJadOmWaKogsLCQshkMpSVlQnTcnJyIJPJoFKp9OY9fPgwXFxc9Kbt3bsXbm5uFinr/RjzPQAgLS0NFRUVUKvViImJgZOTEw4ePGjJIt/Xg0KJlP8XUsZzGpFlMJSIjKGEB3BqOz8/P8hkslZfzRlzp2Tr1q1605ycnMx6p8SQ7xEQEIDly5frfW7gwIGt3im529GjR9GtWzeLNthsusNQUFAgTEtPT7/vnZLm03fv3i2JX+eN+R4t+fDDDzFr1ixzFtEoht4pkeL/Qsp4TiOyDIYSkTGU8ABOlmVMm5LAwEDhfW5urmTalDRvcF9VVQUHB4dW25Tc7dixY3B2drZ4LzK9e/fWayC9adOm+7YpyczMFKYFBwdLph2Dod+jJTNmzMDMmTPNVTSjGdqmRKr/C6niOY3IMhhKRMZQwgM4WYZOp4NKpYJSqYRMJkNJSQlUKlWrF+cpKSlwdXVFRkYGqqqqMHnyZEn1vpWcnAyVSoXAwMD7NmSPj4/H9evXATQGslGjRmHKlCmWKq5gyZIlGDZsGIqLi1FYWIj+/fu32mvV1KlT8eabb6KyshJZWVlwc3OTTI9Phn6PsrIyxMfHo7a2FhqNBvHx8XByckJycrIVSq2voaEBKpUKwcHBmDBhAlQqFerr61ucV8r/C6niOY3IMhhKRMZQwgM4WUbTHZK7X/n5+QCAiIgIeHt7630mLCwMTzzxBLp27SqpcUoiIyPh4eEBBwcH+Pn56Y1Tcvf3mD9/Ptzc3ODo6IinnnoKs2fPRmVlpcXLrFarMWPGDLi4uKBnz56YP3++EAj9/f2xYsUKYd6Kigq88847cHZ2hpubm6TGxjD0exQXF2Po0KHo1q0bunfvDh8fH+zYscOaRReEhobesx/4+fkBsK3/hVTxnEZkGQwlImMo4QGciIjsB89pRJbBUCIyhhIewImIyH7wnEZkGQwlImMo4QGciIjsB89pRJbBUCIyhhIewImIyH7IZDJ06NDBIi9LrsveX6xL26vH5l3+M5SIoEMHhpK21AEREVF7xfOneFiX4rBWPTKUiIChhAcCIiIiU/D8KR7WpTgYSmwYQwkPBERERKbg+VM8rEtxMJTYMIYSHgiIiIhMsWzZMmsXwW6wLsVhrXpkKGnm/fffh0wmQ1ZWllGfYyhhKCEiIiIi0zGU3JGcnAw/Pz+GEhMxlBC1b/X19QgICICbmxt+/PFHgz+nVqsxYsQIuLi4YPfu3WYsIRERSRlDCYDa2loMGDAA2dnZDCUmYighat+io6MxZ84clJeX45lnnjH4czqdDtevX0doaChDCRFRO8ZQAiAoKAiLFi0CAIYSEzGUENmOQYMGIT4+3uTPb9++HVOnTtWbdvbsWcydOxdFRUV4++23jV5mS6HktddeQ2JiosnlJCIi29HuQ8mZM2fQr18/1NbWAjAslCxbtqzVgV+MxVBCRJZUV1eHhx56CFevXjXp8xqNBnK5HEql8p6/jRo1CjKZDBs2bDB6uS2FkkOHDsHHx8ekchJJVVRUFDw9PeHo6IjRo0fj0qVLrc4bHh6Ojh07wsnJSXht3rzZgqWVFrVajY8//hg9evRAr169sGDBAuh0uhbnraiowOTJk+Hs7IzHH38cq1evtnBppcuYevTz80OXLl30tsG6ujqzlMuuQ0lTG5HWXlqtFs8//7zeL4a8U2IahhIi25CRkQFXV1eTPx8dHY1hw4bdM/3KlSvo2LEjXF1dMW7cOL2/lZaWYtiwYfe8FAqFME9LoUSn06FPnz5IS0szubxEUqJUKuHs7IykpCTU1tZi7ty58PX1bXX+8PDw+/69vVm6dCmGDh2KoqIiFBQUoF+/fvj2229bnPf9999HQEAAKioqkJmZiUcffRQxMTEWLrE0GVOPfn5+WLdunUXKZdeh5EHKysrQoUMHuLm5CS+ZTAZXV1eEhYUZvByGEoYSIlsRHh6O0aNHAwC0Wi2Cg4Px3HPP3ffX2uamT58uPO7a3KpVq9CzZ09s3boVHTt2NPpOTGttSqZNm4bFixcbtSwiqQoJCcGECROE95WVlXj44Ydx+vTpFudnKNHXu3dv7NmzR3i/cePGFu+m1tTUoEuXLjhz5owwbfHixXjrrbcsUk6pM7QeAYYSi9HpdLhx44beSyaT4eDBg6ipqTF4OQwlDCVEtmLu3LmYPXs2iouLMW7cOHzwwQfC46uGGDp0KLZu3XrP9IEDB+Kjjz6CSqWCi4sLvvjiC4OX+fbbb8Pd3R3PPPMMFixYoPe31atXY/z48QYvi0jKAgIC7hkDwtvbGxERES3OHx4eDkdHRzzyyCPw8PBAYGAgqqurLVFUybl9+zZkMhkuX74sTDtx4gQ6d+58z6NHGRkZ6NSpE7RarTBt165d8PLyslh5pcqYegQaQ4mrqyt69eoFX19fREZGmq1s7TqUtISPb5mGoYTINowdOxbTpk2Du7u7UXeEm3h5eSE6Olpv2smTJyGTyXD06FEAwN/+9jcMGjRIlPJu3LgRL7zwgijLIrK2sWPHYu3atXrTRowYgfXr17c4/8WLF5GXlwetVovc3FyMHDkS06ZNs0RRJaewsBAymQxlZWXCtJycHMhkMqhUKr15Dx8+DBcXF71pe/fuhZubm0XKKmXG1CMApKWloaKiAmq1GjExMXBycsLBgwfNUjaGEhEwlDCUENkKV1dXuLm5Yfbs2SZ9vqU7JXPmzEHfvn2F90eOHIFMJkNGRkabygrwTgnZjge1YwUa75QsX75c73MDBw5s9U7J3Y4ePYpu3bq12ijZnjX9wl9QUCBMS09Pv++dkubTd+/ezTslMK4eW/Lhhx9i1qxZZikbQ4kIGEoYSohswdWrV/HQQw9BqVTid7/7nV6jT41GA19fX9y6dQsZGRnw9/dv8bGuu9uUNDQ04LHHHoODg8M97fPmzZvX5jKzTQnZk5CQEEycOFF4X1VVBQcHh1bblNzt2LFjcHZ2bpehBGhsC9H8uLVp06b7tinJzMwUpgUHB7NNyR2G1mNLZsyYgZkzZ5qlXAwlImAoYSghsgVxcXHCY1VHjx6Fq6srTp06Jfx9165dmDJlCsaOHYvS0tIWlxEVFYXhw4cL72NjY9G5c2dkZ2frtc/77LPP8Pjjj0Oj0bSpzO7u7sJjYUS2rqn3reTkZKhUKgQGBt63IXt8fDyuX78OAMjPz8eoUaMwZcoUSxVXcpYsWYJhw4ahuLgYhYWF6N+/f6u9Rk2dOhVvvvkmKisrkZWVBTc3N/a+dYeh9VhWVob4+HjU1tZCo9EgPj4eTk5OSE5ONku5GEpEwFDCUEJkC1auXKl3QfPTTz+hd+/euHbtGoDGZ407d+6M1NTUVpfR0NCAp59+Whin5J133sF77713z3xlZWVwcnJq0yCNR44cwZAhQ0z+PJEURUZGwsPDAw4ODvDz89Pr+S4iIgLe3t7C+/nz58PNzQ2Ojo546qmnMHv2bFRWVlqj2JKgVqsxY8YMuLi4oGfPnpg/f75w18jf3x8rVqwQ5q2oqMA777wDZ2dnuLm5cZySZgytx+LiYgwdOhTdunVD9+7d4ePjgx07dpitXAwlImAoYSghsnVlZWUYPXo0wsLCMH369PvOu23btntGdDcHf39/JCQkmH09RERkfQwlImAoYSghsmV1dXXw9/fH8ePHAQBjxoxpccR2IiIic2EoEQFDCUMJEREREZmOoUQEDCUMJURERERkOoYSETCUMJQQERERkekYSkTAUMJQQkRERESmYygRAUMJQwkRERERmY6hRAQMJQwlRERERGQ6hhIRMJQwlBARERGR6RhKRMBQwlBCRERERKZjKBEBQwlDCREREVnG9u3bMXXqVLMt/7XXXkNiYqLZlk8tYygRAUMJQwkRERGZn0ajgVwuh1KpNNs6Dh06BB8fH7Mtn1rGUCIChhKGEiIiIjJdQ0ODQfNFR0dj2LBhZi2LTqdDnz59kJaWZtb1kD6GEhEwlDCUEBERkeF2794NX19fBAYG4rHHHsP7779v0OemT5+ORYsW6U2TyWTIysoS3oeGhmLSpEnCe7lcjlWrVuG5555D165dMX78eJSWluKjjz6Ci4sL+vbte08AmTZtGhYvXtyGb0jGYigRAUMJQwkREREZbvHixejSpQu2bNkCtVqNuro6gz43dOhQbN26VW+aIaHEx8cHhYWFKCsrg7e3Nzw9PREVFQWNRoPg4GA8++yzestcvXo1xo8f34ZvSMZiKBEBQwlDCRERERnu9ddfx6xZs4z+nJeXF6Kjo/WmGRJKNm3aJLxfsGABRo4cKbzPzs5Gp06doNFohGkbN27ECy+8YHT5yHQMJSJgKGEoISIiIsP97ne/w6FDh4z+nKl3SmJjY1v9e35+PmQyGaqqqoRpvFNieQwlImAoYSghIiIiwxQVFaFjx46orq4Wpmk0Gvj6+uLWrVvIyMiAv78/amtr7/lsa21KTpw4IbyfM2dOm0MJ25RYHkOJCBhKGEqIiIjIMImJiejfv/8903ft2oUpU6Zg7NixKC0tbfGzUVFRGD58uN40mUyGGTNmQK1WIzc3F08++SReeeUVqNVqAKaFEnd3dxw9erRN35OMw1AiAoYShhIiIiIyzOeff44///nP90wvLCxE586dkZqa2upnGxoa8PTTT+uNU9IUStzc3DBgwAD88MMP6NmzJ9asWQPA+FBy5MgRDBkypM3fk4zDUCKCtlyQd3JwRicHZxFLY7ynZobjqZnhbVoGQwkRERGZqqysDKNHj0ZYWBimT59+33m3bdumN6L73W1K2srf3x8JCQmiLY8Mw1AigrZckHd2eQydXR4TsTTGk8/bBfm8XW1aBkMJERERmaKurg7+/v44fvw4AGDMmDFGjdgudigh62AoEUFbLsifCU3EM6GJIpbGOmVgKCEiIiJrYCixDwwlImjLBXmfhQrIgxRCMLDGSx6kQJ8ghdXqgIiIiIjaN4YSEbTlgtxjYRz6WDmU9FmogMfCOKvVARERERG1bwwlIuAFOeuAiIiIiEzHUCICXpCzDoiIiIjIdAwlIuAFOeuAiIiIiEzHUCICXpCzDoiIiIjIdAwlIuAFOeuAiIiIiEzHUCICXpCzDoiIiIjIdAwlIuAFOeuAiIiIiEzHUCICXpCzDoiIiIjIdAwlIuAFOeuAiIiIiEzHUCICXpCzDoiIiIjIdAwlIuAFOeuAiIiIiEzHUCICXpCzDoiIiIjIdAwlIuAFOeuAiIiIiEzHUCICXpCzDoiIiIjIdAwlIuAFOeuAiIiIiEzHUCICXpCzDoiIiIjIdAwlIuAFOeuAiIiIiEzHUCICXpCzDoiIiIjIdAwlIuAFOeuAiIiIiEzHUCICXpCzDoiIiIjIdAwlIuAFOeuAiIiIiEzHUCICXpCzDoiIiIjIdAwlIuAFOeuAiIiIiEzHUCICXpCzDoiIiIjIdAwlIuAFOeuAiIiIiEzHUCICXpCzDoiIiIjIdAwlIuAFOeuAiIiIiEzHUCICXpCzDoiIiIjIdAwlAAoKCjBhwgR069YNPXr0wOTJk436PC/IWQdEREREZLp2H0rUajX69u2LFStWoKqqCmq1GidPnjRqGbwgZx0QERERkenafSjZtGkThg8f3qZl8IKcdUBEREREpmv3oWTatGn4y1/+gj/+8Y/o1asXnn/+eRw6dMioZfCCnHVARERERKZr96HklVdeQadOnRATEwO1Wo1t27ahe/fuuHnzZqufWbZsGTp06CC8ZDKZ3ntLv6y9/qYyEBERERGZwq6vJP38/CCTyVp9AcCbb76JF198Ue9zgwYNws8//2zweqx9l8Da65dKGYiIiIjINtl1KDHEkiVLMGrUKL1pgwYNws6dOw1ehrUvyK29fqmUgYiIiIhsU7sPJbm5uXB0dER8fDy0Wi1++uknuLi44MaNGwYvw9oX5NZev1TKQERERES2qd2HEgCIjY1Fv3794OTkBF9fX6Mbui9btsxMJbON9UulDERERERkmxhKiIiIiIjIqhhKiIiIiIjIqhhKiIiIiIjIqhhKiIiIiIjIqhhKRJKVlYVXX30Vrq6ukMlkqKqqMvs61Wo1Pv74Y/To0QO9evXCggULoNPpzL7eJuvWrYOvry+6dOmCSZMmWWy9RERERGRfGEpEolQqsXnzZsTGxloslCxduhRDhw5FUVERCgoK0K9fP3z77bdmX2+TyMhIREdHY9asWQwlRERERGQyhhKR5efnWyyU9O7dG3v27BHeb9y4ET4+PmZf791CQ0MZSoiIiIjIZAwlIrNUKLl9+zZkMhkuX74sTDtx4gQ6d+5s0Ue4AIYSIiIiImobhhKRWSqUFBYWQiaToaysTJiWk5MDmUwGlUpl1nXfjaGEiIiIiNqCocQAfn5+kMlkrb6as/SdkoKCAmFaeno675QQERERkc1hKBGZpduUxMTECO83bdrENiVEREREZHMYSkSi0+mgUqmgVCohk8lQUlIClUpl1rsWS5YswbBhw1BcXIzCwkL079/for1vNTQ0QKVSITg4GBMmTIBKpUJ9fb3F1k9ERERE9oGhRCRNd0jufuXn55ttnWq1GjNmzICLiwt69uyJ+fPnW/TRrdDQ0Hu+r5+fn8XWT0RERET2gaGEiIiIiIisiqGEiIiIiIisiqGEiIiIiIisiqGEiIiIiIisiqGEiIiIiIisiqGEiIiIiIisiqGEiIiIiIisiqGEiIiIiIisiqGEiIiIiIisiqGEiIiIiIisiqFEgry9vZGQkGDy39tq+/btmDp1qlmW/dprryExMdEsyyYiIiIi28RQInFyuRyxsbEWW59Go4FcLodSqTTL8g8dOgQfHx+zLJuIiIiIbBNDicRZOpRER0dj2LBhZlu+TqdDnz59kJaWZrZ1EBEREZFtYShp5vLly+jZsydSUlIAALW1tRgwYAC++eabFueXy+VYuXIlhgwZgu7du+ONN97ArVu3hL/n5OTg5ZdfRo8ePeDl5YUNGzYIf/vqq6/Qu3dvODs7w93dHTt27NBbbmxsLN5991106NABDg4OcHJywpw5c/T+/qB1NM375Zdf4ve//z2cnZ3x8ssvo6ioqNU6mD59OhYtWqQ3TSaTISsrS3gfGhqKSZMmCctftWoVnnvuOXTt2hXjx49HaWkpPvroI7i4uKBv3773BJBp06Zh8eLFrZaBiIiIiNoXhpK77NixA08++SRKSkrw0Ucf4fXXX4dOp2txXrlcDi8vL+Tm5qKqqgoTJ07EhAkTAABqtRp9+/bF4sWLUVdXh5MnT8LV1RXx8fFQKpVwdHQUHpG6fv06srOz9ZbbFDpaulPSNO1+62g+7+DBg3H58mXU1H04lEIAAAN4SURBVNTAz88P8+bNa/X7Dx06FFu3btWb9qBQ4uPjg8LCQpSVlcHb2xuenp6IioqCRqNBcHAwnn32Wb3lrV69GuPHj2+1DERERETUvjCUtOCDDz7AM888g8cff/y+dxXkcjnWrFkjvFcqlejQoQNqampw5MgR9OzZE2q1Wvh7UFAQ3n33XeTl5cHBwQFRUVGora1tcbmGhJL7raP5vJs2bRLef//99/Dz82v1O3l5eSE6Olpv2oNCSfPlL1iwACNHjhTeZ2dno1OnTtBoNMK0jRs34oUXXmi1DERERETUvjCUtOD48eOQyWRYuHDhfeeTy+X497//LbxXqVSQyWTIy8vDzp07MWjQIL35169fj1GjRgEAdu7ciZdeegndunXD+PHjcf78eb3lGhJKHrSOlj4fHh4OX1/fVr+TKXdKmi+/+d8AID8/HzKZDFVVVcI03ikhIiIiouYYSu5SV1eHwYMH469//StcXFxw7ty5Vue9+07JhQsX7rlT0tDQIPz97rsYAFBdXY2ZM2fixRdf1Ftu04W+u7v7A++U3G8dxoaS1tqUnDhxQng/Z86cNoUStikhIiIiouYYSu4ye/ZsjBs3DjqdDqtWrcKQIUNQV1fX4rxyuRz9+vVDXl4eqqur8ac//QlvvfUWgMY2JV5eXggJCUF9fT0yMjLwyCOPIC4uDkqlEsnJyVCpVGhoaMCCBQv0HqlqfqE/fPjwexraN29T0to6WloW8OBQEhUVheHDh+tNk8lkmDFjBtRqNXJzc/Hkk0/ilVdegVqtNimUuLu74+jRo62WgYiIiIjaF4aSZuLj4/HII4/g2rVrAACtVouxY8cKvV7dTS6XY8WKFRgyZAi6deuGP/zhD3ptUM6fP49x48bBxcUFnp6eCAsLAwCcOXMGzz//PJydneHi4oIxY8a0+vhWTEwM5HI5XFxchAbqzf/e2jpaWhbw4FDS0NCAp59+Wm+ckqZQ4ubmhgEDBuCHH35Az549sWbNGqNDyZEjRzBkyJBW109ERERE7Q9DSRtYegwRS9m2bZveiO53tylpC39/f7OORk9EREREtoehpA3sNZTcTcxQQkRERER0N4aSNmAoISIiIiJqO4YSIiIiIiKyKoYSIiIiIiKyKoYSIiIiIiKyKoYSIiIiIiKyKoYSIiIiIiKyKoYSIiIiIiKyqv8PdPba9/TVM0cAAAAASUVORK5CYII=\" width=\"898.9166666666666\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "FRAME_SKIP = 3\n",
    "trajectory = trajectories[0]\n",
    "ns, xs, ys, state, edge = zip(*trajectory)\n",
    "\n",
    "plt.rcParams['figure.figsize'] = [9, 5]\n",
    "fig, (ax1, ax2, ax3) = plt.subplots(1, 3, constrained_layout=True)\n",
    "\n",
    "plt.subplot(1, 3, 1)\n",
    "ax1.set_xlim(-1, 1)\n",
    "ax1.set_ylim(-6, 6)\n",
    "for edge in bar_frame.edges:\n",
    "    x, y = [edge.xs, edge.ys]\n",
    "    if edge.layer == 0:\n",
    "        plt.plot(x, y, color='k')\n",
    "    else:\n",
    "        plt.plot(x, y, color='C'+str(edge.layer-1))\n",
    "plt.xlabel(r'x position ($\\mathrm{\\mu m}$)')\n",
    "plt.ylabel(r'y position ($\\mathrm{\\mu m}$)')\n",
    "line1, = ax1.plot([], [], 'bo')\n",
    "ax1.set_aspect('equal')\n",
    "\n",
    "plt.subplot(1, 3, 2)\n",
    "plt.plot(k[0], k[1])\n",
    "plt.xlabel(r'$k_x\\ (\\mathrm{\\AA^{-1}})$')\n",
    "plt.ylabel(r'$k_y\\ (\\mathrm{\\AA^{-1}})$')\n",
    "line2, = ax2.plot([], [], color='black', linestyle='dashdot', marker='o')\n",
    "ax2.set_aspect('equal')\n",
    "\n",
    "plt.subplot(1, 3, 3)\n",
    "line3, = ax3.plot([], [], color='black', linestyle='dashdot', marker='o')\n",
    "plt.plot(bar_sim._bandstructure.r[0], bar_sim._bandstructure.r[1])\n",
    "plt.xlabel(r'$r_x\\ (\\mathrm{\\mu m})$')\n",
    "plt.ylabel(r'$r_y\\ (\\mathrm{\\mu m})$')\n",
    "ax3.set_aspect('equal')\n",
    "\n",
    "def animate(frame_num):\n",
    "    frame_num_adapted = frame_num*FRAME_SKIP\n",
    "    line1.set_data((xs[frame_num_adapted], ys[frame_num_adapted]))\n",
    "    line2.set_data(k[0][ns[frame_num_adapted][0]], k[1][ns[frame_num_adapted][0]])\n",
    "    line3.set_data(bar_sim._bandstructure.r[0][ns[frame_num_adapted][0]], bar_sim._bandstructure.r[1][ns[frame_num_adapted][0]])\n",
    "\n",
    "anim = FuncAnimation(fig, animate, frames=int(len(xs)/FRAME_SKIP), interval=15)\n",
    "plt.show()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d705d8eb",
   "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.7.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}