0a_Minimum_working_example.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# OUTDATED, the examples moved to the gallery\n",
    "## See https://empymod.github.io/emg3d-gallery\n",
    "\n",
    "----\n",
    "\n",
    "# Minimum working example\n",
    "\n",
    "This is a simple minimum working example to get started, along the lines of the one given in https://emg3d.readthedocs.io/en/stable/usage.html#example.\n",
    "\n",
    "To see some more realistic models have a look at the other notebooks in this repo.\n",
    "\n",
    "-------------------------------------------------------------------------------\n",
    "This notebooks uses `discretize` to create meshes easily and plot the model as well as the resulting electric field, which also requires `matplotlib`. If you are interested in a basic example that only requires `emg3d` here it is:\n",
    "\n",
    "```py\n",
    "import emg3d\n",
    "import numpy as np\n",
    "\n",
    "# Create a simple grid, 8 cells of length 1 in each direction, starting at the origin.\n",
    "grid = emg3d.utils.TensorMesh([np.ones(8), np.ones(8), np.ones(8)], x0=np.array([0, 0, 0]))\n",
    "\n",
    "# The model is a fullspace with tri-axial anisotropy.\n",
    "model = emg3d.utils.Model(grid, res_x=1.5, res_y=1.8, res_z=3.3)\n",
    "\n",
    "# The source is a x-directed, horizontal dipole at (4, 4, 4), frequency is 10 Hz.\n",
    "sfield = emg3d.utils.get_source_field(grid, src=[4, 4, 4, 0, 0], freq=10.0)\n",
    "\n",
    "# Calculate the electric signal.\n",
    "efield = emg3d.solver.solver(grid, model, sfield, verb=3)\n",
    "\n",
    "# Get the corresponding magnetic signal.\n",
    "hfield = emg3d.utils.get_h_field(grid, model, efield) \n",
    "```\n",
    "-------------------------------------------------------------------------------\n",
    "\n",
    "**Requires**\n",
    "- **emg3d >= 0.9.1**\n",
    "- discretize, matplotlib\n",
    "\n",
    "First, we load `emg3d` and `discretize` (to create a mesh), along with `numpy` and `matplotlib`:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import emg3d\n",
    "import discretize\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "from matplotlib.colors import LogNorm"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib notebook\n",
    "plt.style.use('ggplot')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 1. Mesh\n",
    "\n",
    "First, we define the mesh (see `discretize.TensorMesh` for more info). In reality, this task requires some careful considerations. E.g., to avoid edge effects, the mesh should be large enough in order for the fields to dissipate, yet fine enough around source and receiver to accurately model them. This grid is too small, but serves as a minimal example."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "  <tr>\n",
       "    <td style='font-weight: bold; font-size: 1.2em; text-align: center;' colspan='3'>TensorMesh</td>\n",
       "    <td style='font-size: 1.2em; text-align: center;'colspan='4'>49,152 cells</td>\n",
       "  </tr>\n",
       "  <tr>\n",
       "    <th></th>\n",
       "    <th></th>\n",
       "    <th colspan='2' style='padding: 5px 20px 5px 20px;'>MESH EXTENT</th>\n",
       "    <th colspan='2' style='padding: 5px 20px 5px 20px;'>CELL WIDTH</th>\n",
       "    <th style='padding: 5px 20px 5px 20px;'>FACTOR</th>\n",
       "  </tr>\n",
       "  <tr>\n",
       "    <th style='padding: 5px 20px 5px 20px;'>dir</th>\n",
       "    <th style='padding: 5px 20px 5px 20px;'>nC</th>\n",
       "    <th style='padding: 5px 20px 5px 20px;'>min</th>\n",
       "    <th style='padding: 5px 20px 5px 20px;'>max</th>\n",
       "    <th style='padding: 5px 20px 5px 20px;'>min</th>\n",
       "    <th style='padding: 5px 20px 5px 20px;'>max</th>\n",
       "    <th style='padding: 5px 20px 5px 20px;'>max</th>\n",
       "  </tr>\n",
       "  <tr>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>x</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>48</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>-662.16</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>662.16</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>25.00</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>37.01</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>1.04</td>\n",
       "  </tr>\n",
       "  <tr>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>y</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>32</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>-857.96</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>857.96</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>50.00</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>63.34</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>1.03</td>\n",
       "  </tr>\n",
       "  <tr>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>z</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>32</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>-540.80</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>540.80</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>30.00</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>44.32</td>\n",
       "    <td style='padding: 5px 20px 5px 20px;'>1.05</td>\n",
       "  </tr>\n",
       "</table>\n"
      ],
      "text/plain": [
       "\n",
       "  TensorMesh: 49,152 cells\n",
       "\n",
       "                      MESH EXTENT             CELL WIDTH      FACTOR\n",
       "  dir    nC        min           max         min       max      max\n",
       "  ---   ---  ---------------------------  ------------------  ------\n",
       "   x     48       -662.16        662.16     25.00     37.01    1.04\n",
       "   y     32       -857.96        857.96     50.00     63.34    1.03\n",
       "   z     32       -540.80        540.80     30.00     44.32    1.05\n"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "grid = discretize.TensorMesh(\n",
    "        [[(25, 10, -1.04), (25, 28), (25, 10, 1.04)],\n",
    "         [(50, 8, -1.03), (50, 16), (50, 8, 1.03)],\n",
    "         [(30, 8, -1.05), (30, 16), (30, 8, 1.05)]],\n",
    "        x0='CCC')\n",
    "grid"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 2. Model\n",
    "\n",
    "Next we define a very simple fullspace model with $\\rho_x=1.5\\,\\Omega\\rm{m}$,  $\\rho_y=1.8\\,\\Omega\\rm{m}$, and  $\\rho_z=3.3\\,\\Omega\\rm{m}$."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "model = emg3d.utils.Model(grid, res_x=1.5, res_y=1.8, res_z=3.3)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can plot the model using `discretize`; in this case it is obviously rather a boring plot, as it shows a homogenous fullspace. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\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('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",
       "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 = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(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 (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.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 = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\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 nav_element = $('<div/>');\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\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",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option);\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\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",
       "\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",
       "\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]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.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, 0, fig.canvas.width / mpl.ratio, fig.canvas.height / mpl.ratio);\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",
       "    {\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.updated_canvas_event = function() {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message(\"ack\", {});\n",
       "}\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = \"image/png\";\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src);\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\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(\"No handler for the '\" + msg_type + \"' message type: \", msg);\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(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\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",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\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",
       "    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",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\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",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, {key: value,\n",
       "                             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",
       "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\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.close = function() {\n",
       "        comm.close()\n",
       "    };\n",
       "    ws.send = function(m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function(msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['data'])\n",
       "    });\n",
       "    return ws;\n",
       "}\n",
       "\n",
       "mpl.mpl_figure_comm = function(comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = $(\"#\" + id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm)\n",
       "\n",
       "    function ondownload(figure, format) {\n",
       "        window.open(figure.imageObj.src);\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy,\n",
       "                           ondownload,\n",
       "                           element.get(0));\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.get(0);\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",
       "\n",
       "    var output_index = fig.cell_info[2]\n",
       "    var cell = fig.cell_info[0];\n",
       "\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
       "    var width = fig.canvas.width/mpl.ratio\n",
       "    fig.root.unbind('remove')\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).html('<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/mpl.ratio\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<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 () { fig.push_to_output() }, 1000);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>');\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\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) { continue; };\n",
       "\n",
       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
       "    buttongrp.append(button);\n",
       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
       "    titlebar.prepend(buttongrp);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(el){\n",
       "    var fig = this\n",
       "    el.on(\"remove\", function(){\n",
       "\tfig.close_ws(fig, {});\n",
       "    });\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
       "    // this is important to make the div 'focusable\n",
       "    el.attr('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",
       "    }\n",
       "    else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\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",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which == 13) {\n",
       "        this.canvas_div.blur();\n",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\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('matplotlib', mpl.mpl_figure_comm);\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4nOy9fVRU56H/6x/3rnvXvet3173r3nXv77ymPT09p+f9tb/29Jz2nPa0iTVpkqbdwzAhoAIjDAgE6agYQQ2IRAO+TOILKiJQMscAptgTE0GjwgGaolGZ6IA6GREVKVJFoA7B7/0ji2kmjHFv3LPn2fv5ftb6/OHMnplnMn2e59N52cwDIYQQQgiRinnxHgAhhBBCCDEWBiAhhBBCiGQwAAkhhBBCJIMBSAghhBAiGQxAQgghhBDJYAASQgghhEgGA5AQQgghRDIYgIQQQgghksEAJIQQQgiRDAYgIYQQQohkMAAJIYQQQiSDAUgIIYQQIhkMQEIIIYQQyWAAEkIIIYRIBgOQEEIIIUQyGICEEEIIIZLBACSEEEIIkQwGICGEEEKIZDAACSGEEEIkgwFICCGEECIZDEBCCCGEEMlgABJCCCGESAYDkBBCCCFEMhiAhBBCCCGSwQAkhBBCCJEMBiAhhBBCiGQwAAkhhBBCJIMBSAghhBAiGQxAQgghhBDJYAASQgghhEgGA5AQQgghRDIYgIQQQgghksEAJIQQQgiRDAYgIYQQQohkMAAJIYQQQiSDAUgIIYQQIhkMQEIIIYQQyWAAEkIIIYRIBgOQEEIIIUQyGICEEEIIIZLBACSEEEIIkQwGICGEEEKIZDAACSGEEEIkgwFICCGEECIZDEBCCCGEEMlgABJCCCGESAYDkBBCCCFEMhiAhBBCCCGSwQAkhBBCCJEMBiAhhBBCiGQwAAkhhBBCJIMBSAghhBAiGQzAOBEKheD1ehEKheI9lEfCCs+Dz0EcrPA8+BzM//h6wuciLj6fD2VlZXA6nVAUBd3d3Z97fG9vLxRFmeXVq1fDx1y5cgUbN26Ey+WCoig4dOjQrPvxer2z7iMtLU335/cwGIBxYnx8HIqiYHx8PN5DeSSs8Dz4HMTBCs+Dz8H8j68nfC7icurUKTQ0NKCrq0tTAA4ODmJ0dDTs9PR0+Jj+/n7s378f7e3tSE9Pf2AA5ufnR9zH7du3dX9+D4MBGCesMpGs8Dz4HMTBCs+Dz8H8j68nfC7mQEsA3r17V9V9ulyuBwZgQUHBnMapJwzAOGGViWSF58HnIA5WeB58DuZ/fD3hczEHWgLQ5XIhPT0da9euxblz5x54/OcFYFJSEpxOJ1wuFyorK3Hjxo1Hfg5aYQDGiYmJCbjdbkxMTMR7KI+EFZ4Hn4M4WOF58DmY//H1hM/FWEKhEMbHxyNU851FNQE4ODiII0eO4NKlS/D7/aiqqoLNZoPP54t6/IMC8NSpU+js7EQwGMSZM2dQXFyMtLQ03LlzR92T1AkGICEWo6enB/PmzUNPT0+8h0KItHAePjr3p7W/0zg1NYWMjIyIH1h4vd6H3k5NAEajrKwMGzZsiHrdgwLws0xOTiItLQ0tLS2aH/9RYADqQO7bh/DUT/dTKoSt7e2YN28eWtvb4z4WKqe5bz9807P6WirDPJzr66yFj0dfxNTw06r8ePRFAJ+8UxmLdwCj0djYiLy8vKjXqQ1AAFi3bh127dql+fEfBQagDjz10/344pZNlArh4ZMnMG/ePBw+eSLuY6Fy+tRP90u/lsowD+f6Omvh3s2n8JtrX1TlvZtPzflx5hqAmzZtwpo1a6JepzYAQ6EQlixZggMHDmh+/EeBAagDVlq0qPmVYeOhYssAlGMeGhGAkzefxMS1L6hy8uaT2u57chKBQACBQACKoqClpQWBQADDw8MAgPr6emzbti18/KFDh9Dd3Y1r167hypUrqK+vh6Io6OrqCh8zNTUVvk+n04n9+/cjEAjg+vXr4WNqamrg8/kwNDSEvr4+lJWVITk5GTdv3nzE/1raYADqgJUWLWp+Zdh4qNgyAOWYh0YE4MTNBbh77Q9VOXFzgab7ftCJnT0eDwDA4/GguLg4fPzBgweRnZ0Nh8OBhQsXYvXq1bO+4zk0NBT1Pj99P5WVlXA6nbDb7XA6ndi4cSMGBgbm/N9orjAAdcBKixY1vzJsPFRsGYByzEMjAnDs5vdxe/APVDl28/sxH4+VYADqgJUWLWp+Zdh4qNgyAOWYh0YE4J2b8zE6+PuqvHNzfszHYyUYgDpgpUWLml8ZNh4qtgxAOeahEQE4evMJ/Grw91Q5evOJmI/HSjAAdcBKixY1vzJsPFRsGYByzEMjAvDWzcdxc/B3VXnr5uMxH4+VYADqgJUWLWp+Zdh4qNgyAOWYh0YE4PDQ47h+9XdUOTzEANSCKQPQ5XJF/ZVNVVUVAKC4uHjWdZWVlRH3MTY2hq1btyI5ORnJycnYunWr6j/w/FmstGhR8yvDxkPFlgEoxzw0IgBvDj2Owau/o8qbDEBNmDIAb9++jdHR0bBnzpyBoijo7e0F8EkA7tixI+KYz/7x6tLSUuTn58Pv98Pv9yM/Px9lZWVzGo+VFi1qfmXYeKjYMgDlmIdGBODgje8hMPDfVTl443sxH4+VMGUAfpbq6mpkZ2fj/v37AD4JwOrq6gcePzAwAEVR0NfXF77M7/dDURQMDg5qfnwrLVrU/Mqw8VCxZQDKMQ+NCMCBG9/DpYH/rsoBBqAmTB+AU1NTWLx4MRobG8OXFRcXY/HixVi8eDFefPFF1NTUYGJiInx9W1sbUlJSZt1XSkoKjh49qnkMVlq0qPmVYeOhYssAlGMeGhGAwRuPo2/gd1QZvMGPgLVg+gDs6OhAQkICRkZGwpcdOXIEZ86cQTAYRHt7O1wuF9atWxe+vrGxETk5ObPuKycnB01NTQ98rFAoFPEHpmei0kqLFjW/Mmw8VGwZgHLMw5nXeWJiImJvDIVCc3r9oxG4/jguXPkdVQauMwC1YPoALCkpeeh39y5dugRFUXDp0iUADw7ApUuXorm5+YH34/V6I35Y4na7AVhr0aLmV4aNh4otA1COeTjzOrvd7oi90ev1zun1j8al64/Dd+V3VXmJAagJUwfgzZs3YbPZ8Itf/OJzj7t//z7sdjs6OjoAzP0jYL4DSM2gDBsPFVsGoBzz0Ih3APuuP4GzV35PlX3XeSJoLZg6AL1eL9LT0/Hxxx9/7nHBYBCKosDn8wH47Y9A+vv7w8f09fXxRyDUEsqw8VCxZQDKMQ+N+A6g/9p8fBD8fVX6r/FPwWnBtAE4PT2NzMxM1NXVRVx+/fp1HDhwABcvXsTQ0BB6enqQl5cHt9uN6enp8HGlpaUoKCgInwZm2bJlPA0MtYQybDxUbBmAcsxDIwLww2vfR0/wD1T54bXvx3w8VsK0AfjBBx9EfcdueHgYRUVFWLRoERITE5GdnY29e/dibGws4rixsTFs2bIlfCLoLVu28ETQ1BLKsPFQsWUAyjEPjQjA3msL8IvgY6rsvbYg5uOxEqYNQJGw0qJFza8MGw8VWwagHPPQiAA8e+1JdH30BVWevfZkzMdjJRiAOmClRYuaXxk2Hiq2DEA55qERAfjBtafQ8dEXVfnBtadiPh4rwQDUASstWtT8yrDxULFlAMoxD40IwJ7Bp3Ei8CVV9gw+HfPxWAkGoA5YadGi5leGjYeKLQNQjnloRAD+cvAZvBf4sip/OfhMzMdjJRiAOmClRYuaXxk2Hiq2DEA55qERAdg9+CzaAn+iyu7BZ2M+HivBANQBKy1a1PzKsPFQsWUAyjEPjQjAzqvP4d3LX1Fl59XnYj4eK8EA1AErLVrU/Mqw8VCxZQDKMQ+NCMCOq8/h8OU/U2UHA1ATDEAdsNKiRc2vDBsPFVsGoBzz0IgAPHn1x/j55b9Q5cmrP475eKwEA1AHrLRoUfMrw8ZDxZYBKMc8NCIAjw/Y8LNLf6XK4wO2mI/HSjAAdcBKixY1vzJsPFRsGYByzEMjArDtih2NF/9WlW1X7DEfj5VgAOqAlRYtan5l2Hio2DIA5ZiHRgTgu1cc+I+Lf6/Kd684Yj4eK8EA1AErLVrU/Mqw8VCxZQDKMQ+NCMB3rjyPN/r/UZXvXHk+5uOxEgxAHbDSokXNrwwbDxVbBqAc89CIAPzPKy+gvv9/qPI/r7wQ8/FYCQagDlhp0aLmV4aNh4otA1COeWhEAB4KJmN/39dUeSiYHPPxWAkGoA5YadGi5leGjYeKLQNQjnloRAD+LJiCfX3/pMqfBVNiPh4rwQDUASstWtT8yrDxULFlAMoxD40IwIMfLcIe/z+r8uBHi2I+HivBANQBKy1a1PzKsPFQsWUAyjEPjQjApo8Wo8r/L6ps+mixpvv2+XwoKyuD0+mEoijo7u7+3ON7e3uhKMosr169Gj7mypUr2LhxI1wuFxRFwaFDh6Le1+HDh+FyueBwOOB2u/Hhhx9qGrseMAB1wEqLFjW/Mmw8VGwZgHLMQyMC8M1AKnZc+JYq3wykarrvU6dOoaGhAV1dXZoCcHBwEKOjo2Gnp6fDx/T392P//v1ob29Henp61ADs6OiA3W5Ha2srBgYGUF1djaSkJAwPD2sa/6PCANQBKy1a1PzKsPFQsWUAyjEPjQhAbyAdr134N1V6A+lzfhwtAXj37l1V9+lyuaIG4MqVK7Fr166Iy/Ly8lBfX69+wDrAANQBKy1a1PzKsPFQsWUAyjEPjQjANy47se38t1X5xmUnAGBiYgLj4+NhQ6HQQx9HSwC6XC6kp6dj7dq1OHfu3AOPjxaAU1NTSEhImPVYe/fuRVFR0UPHqScMQB2w0qJFza8MGw8VWwagHPPQiACsv7wEm8//uyrrLy8BALjd7ojv6Hm93oc+jpoAHBwcxJEjR3Dp0iX4/X5UVVXBZrPB5/NFPT5aAI6MjEBRFFy4cCHi8sbGRuTk5Dx0nHrCANQBKy1a1PzKsPFQsWUAyjEPjQjAusuZqPjwe6qsu5wJIHbvAEajrKwMGzZsiHrd5wWg3++PuLyxsRG5ubmaH/9RYADqgJUWLWp+Zdh4qNgyAOWYh0YEYM1lFzZ++Lgqay675vw4cw3AxsZG5OXlRb2OHwFLgJUWLWp+Zdh4qNgyAOWYh0YE4L5LWSj3zVflvktZc36cuQbgpk2bsGbNmqjXfd6PQKqqqiIu449ATIqVFi1qfmXYeKjYMgDlmIdGBOCeS9lY7/u+KvdcytZ035OTkwgEAggEAlAUBS0tLQgEAuHTsdTX12Pbtm3h4w8dOoTu7m5cu3YNV65cQX19PRRFQVdXV/iYqamp8H06nU7s378fgUAA169fDx8zcxqYtra2iNPA3Lx58xH/a2mDAagDVlq0qPmVYeOhYssAlGMeGhGAuy/moKT3SVXuvqjtRxQPOrGzx+MBAHg8HhQXF4ePP3jwILKzs+FwOLBw4UKsXr0aPT09Efc5NDQU9T4/fT/Ab08EnZiYCLfb/cAfksQSUwag1+ud9R83LS0tfP39+/fh9XrhdDrhcDhQXFyMK1euRNzH2NgYtm7diuTkZCQnJ2Pr1q2qz+3zWay0aFHzK8PGQ8WWASjHPDQiAHf252LtuR+ocme/sT+iMDumDcD8/PyIM3Hfvn07fH1zczOSk5PR1dWFYDCIyspKOJ1OTExMhI8pLS1Ffn4+/H4//H4/8vPzUVZWNqfxWGnRouZXho2Hii0DUI55aEQAvtaXj5fOPqvK1/ryYz4eK2HaACwoKIh63f3795Geno7m5ubwZaFQCCkpKXj33XcBAAMDA1AUBX19feFj/H5/+E+8aMVKixY1vzJsPFRsGYByzEMjAtDTtwyFZ36oSk/fspiPx0qYNgCTkpLgdDrhcrlQWVmJGzduAABu3LgBRVFw+fLliNuUl5eHv8zZ1taGlJSUWfebkpKCo0ePPvBxQ6FQxLmFZt5RtNKiRc2vDBsPFVsGoBzzcOZ1nst599Sypa8Ay8/8SJVb+qK/MUSiY8oAPHXqFDo7OxEMBnHmzBkUFxcjLS0Nd+7cwYULF6AoCkZGRiJus2PHDpSUlAB48Bm3c3Jy0NTU9MDH/ex3D91uNwBrLVrU/Mqw8VCxZQDKMQ9nXue5/OUNtWz2/wQ/+eDHqtzs/4lujysDpgzAzzI5OYm0tDS0tLSEA/DWrVsRx2zfvv2hAbh06dKIj44/C98BpGZQho2Hii0DUI55aMQ7gBX+5Vj2gU2VFf7luj2uDFgiAAFg3bp12LVrV0w/An4QVlq0qPmVYeOhYssAlGMeGvEdwE0XliPvdIIqN11gAGrBEgEYCoWwZMkSHDhwIPwjkIMHD4avn5qaivojkP7+/vAxfX19/BEItYQybDxUbBmAcsxDIwLwlfMrkHPKrspXzq+I+XishCkDsKamBj6fD0NDQ+jr60NZWRmSk5PDZ9Fubm5GSkoKuru7EQwGsXnz5qingSkoKAifBmbZsmU8DQy1hDJsPFRsGYByzEMjAnDDh4XI7nGocsOHhTEfj5UwZQDOnNfPbrfD6XRi48aNGBgYCF8/cyLo9PR0OBwOFBUVIRgMRtzH2NgYtmzZEj4R9JYtW3giaGoJZdh4qNgyAOWYh0YE4PoPVyGz53lVrv9wVczHYyVMGYCiYaVFi5pfGTYeKrYMQDnmoREBWOJbhSW/TFJliY8BqAUGoA5YadGi5leGjYeKLQNQjnloRAC+7FuN9PeTVfmyb3XMx2MlGIA6YKVFi5pfGTYeKrYMQDnmoREBuNZXhNT3U1S51lcU8/FYCQagDlhp0aLmV4aNh4otA1COeWhEABb3FmPhLxaqsri3OObjsRIMQB2w0qJFza8MGw8VWwagHPPQiABcfW4NkrsXq3L1uTUxH4+VYADqgJUWLWp+Zdh4qNgyAOWYh0YE4Kqza5HUlarKVWfXxnw8VoIBqANWWrSo+ZVh46FiywCUYx4aEYCFZ9fB0ZWmysKz62I+HivBANQBKy1a1PzKsPFQsWUAyjEPjQjAFWdehr0zXZUrzrwc8/FYCQagDlhp0aLmV4aNh4otA1COeWhEAC4/UwLbfy1R5fIzJTEfj5VgAOqAlRYtan5l2Hio2DIA5ZiHRgTgsg/W47mOTFUu+2B9zMdjJRiAOmClRYuaXxk2Hiq2DEA55qERAZh/ugzPtrtUmX+6LObjsRIMQB2w0qJFza8MGw8VWwagHPPQiADMO1WGp09mqTLvFANQCwxAHbDSokXNrwwbDxVbBqAc89CIAMzt2YCnTmSrMrdnQ8zHYyUYgDpgpUWLml8ZNh4qtgxAOeahEQG4tKccC44vVeXSnvKYj8dKMAB1wEqLFjW/Mmw8VGwZgHLMQyMCMKvnFcw/nqPKrJ5XYj4eK8EA1AErLVrU/Mqw8VCxZQDKMQ+NCEDXLzfi8fdyVen65caYj8dKMAB1wEqLFjW/Mmw8VGwZgHLMQyMCMOP9TfjusTxVZry/KebjsRIMQB2w0qJFza8MGw8VWwagHPPQiABc8otN+M7RF1W55BcMQC0wAHXASosWNb8ybDxUbBmAcsxDIwIwvftV/FtbvirTu1+N+XisBANQB6y0aFHzK8PGQ8WWASjHPDQiANO6K/CvrctUmdZdEfPxWAkGoA5YadGi5leGjYeKLQNQjnloRACmdlfim60FqkztrtR03z6fD2VlZXA6nVAUBd3d3Z97fG9vLxRFmeXVq1cjjuvs7EReXh4SExORl5c36349Hs+s+ygsLNQ0dj1gAOqAlRYtan5l2Hio2DIA5ZiHRgTgoq7N+OcjP1Hloq7Nmu771KlTaGhoQFdXl6YAHBwcxOjoaNjp6enwMX6/HwkJCWhqasLVq1fR1NQEu92Ovr6+8DEejwelpaUR9zE2NqbtP4wOMAB1wEqLFjW/Mmw8VGwZgHLMQyMCcGHnZnzjXbcqF3ZqC8BPoyUA7969+8BjKioqUFpaGnFZSUkJKit/++6kx+NBeXn8T1rNANQBKy1a1PzKsPFQsWUAyjEPjQjA5P/agq+/s1yVyf+1BQAwMTGB8fHxsKFQ6KGPoyUAXS4X0tPTsXbtWpw7dy7imIyMDLS0tERc1tLSgszMzPC/PR4PUlJSkJqaipycHGzfvh2//vWv1f4n0Q0GoA5YadGi5leGjYeKLQNQjnloSAB2bMXXDq9QZXLHVgCA2+2O+H6d1+t96OOoCcDBwUEcOXIEly5dgt/vR1VVFWw2G3w+X/gYu92OkydPRtzu5MmTSExMDP+7o6MDPT09CAaDeP/991FQUID8/HxVoaonpgzApqYmrFixAi+88AJSU1NRXl6OwcHBiGOKi4tnfcny02/BAsDY2Bi2bt2K5ORkJCcnY+vWrZ/71u6DsNKiRc2vDBsPFVsGoBzz0IgATOrYhq++vVKVSR3bAMTuHcBolJWVYcOGDeF/RwvAEydOwOFwPPA+bt26Bbvdjq6uLs2P/yiYMgBLSkpw7NgxXLlyBYFAAGVlZcjMzMTk5GT4mOLiYuzYsSPiS5bj4+MR91NaWor8/Hz4/X74/X7k5+ejrKxM83istGhR8yvDxkPFlgEoxzw0IgCfb9+Gf/zPlap8vn3bnB9nrgHY2NiIvLy88L/VfAQcjaVLl6K5uVnz4z8KpgzAz3L79m0oihLxNmxxcTGqq6sfeJuBgQEoihLxyxy/3x/+hY8WrLRoUfMrw8ZDxZYBKMc8NCIAE0968Hc/X6XKxJOeOT/OXANw06ZNWLNmTfjfFRUVWL9+fcQxpaWlsz6B/DR37tyBw+HAe++9p/nxHwVLBOD169ehKAqCwWD4suLiYixevBiLFy/Giy++iJqaGkxMTISvb2trQ0pKyqz7SklJwdGjRzU9vpUWLWp+Zdh4qNgyAOWYh0YEYMKJ1/A3h15SZcKJ1zTd9+TkJAKBAAKBABRFQUtLCwKBAIaHhwEA9fX12Lbtt+8qHjp0CN3d3bh27RquXLmC+vp6KIoS8dHthQsXkJCQgObmZly9ehXNzc0Rp4GZnJxETU0N/H4/hoaG0Nvbi1WrVmHJkiURjWIEpg/A+/fvY8OGDVi9enXE5UeOHMGZM2cQDAbR3t4Ol8uFdevWha9vbGxETk7OrPvLyclBU1NT1McKhUIR3yuYebGstGhR8yvDxkPFlgEoxzyceZ3n8p07tdiOv4a/bnlJlbbj2gLwQSd29ng+eSfR4/GguLg4fPzBgweRnZ0Nh8OBhQsXYvXq1ejp6Zl1v52dncjNzYXdbkdeXl5EIN67dw8lJSVITU2F3W5HZmYmPB5PODqNxPQBWFVVBZfLhV/96lefe9ylS5egKAouXboE4MEB+Hmfw3u93oj/kbjdbgDWWrSo+ZVh46FiywCUYx7OvM5z+dWtWmzHX8df/Wy1Km3HX9ftcWXA1AG4Z88eZGRkYGho6KHH3r9/H3a7HR0dHQDm9hEw3wGkZlCGjYeKLQNQjnloxDuAynvb8ZdvFalSeW+7bo8rA6YMwPv372P37t1wOp24du2aqtsEg8GIH4rM/Aikv78/fExfXx9/BEJNrwwbDxVbBqAc89CI7wD+6NgO/PnBYlX+6NiOmI/HSpgyAKuqqpCSkgKfzxdxmpd79+4B+ORHIQcOHMDFixcxNDSEnp4e5OXlwe12R/zNvtLSUhQUFIRPA7Ns2TKeBoaaXhk2Hiq2DEA55qERAfjc0R34s+ZiVT53lAGoBVMGYLQvbSqKgmPHjgEAhoeHUVRUhEWLFiExMRHZ2dnYu3fvrD+2PDY2hi1btoRPBL1lyxaeCJqaXhk2Hiq2DEA55qEhAdi2A19pWqPK59oYgFowZQCKhpUWLWp+Zdh4qNgyAOWYh0YE4A/bduJPG9eq8odtO2M+HivBANQBKy1a1PzKsPFQsWUAyjEPjQjAZ9t24k8a16ryWQagJhiAOmClRYuaXxk2Hiq2DEA55qERAfhM6y58+c11qnymdVfMx2MlGIA6YKVFi5pfGTYeKrYMQDnmoSEBeGQXvnxgnSqfOcIA1AIDUAestGhR8yvDxkPFlgEoxzw0IgCffrcKf/wfL6vy6XerYj4eK8EA1AErLVrU/Mqw8VCxZQDKMQ8NCcB3qvAl78uqfPodBqAWGIA6YKVFi5pfGTYeKrYMQDnmoREB+IN3duOP3ihR5Q/e2R3z8VgJBqAOWGnRouZXho2Hii0DUI55aEgAHt6NP2ooVeUPDjMAtcAA1AErLVrU/Mqw8VCxZQDKMQ+NCMCn3t6NL/60VJVPvc0A1AIDUAestGhR8yvDxkPFlgEoxzw0JAD/cze+WF+qyqf+kwGoBQagDlhp0aLmV4aNh4otA1COeWhEAD758z34Qu16VT758z0xH4+VYADqgJUWLWp+Zdh4qNgyAOWYh0YE4IJDe/HY/jJVLji0N+bjsRIMQB2w0qJFza8MGw8VWwagHPPQkABs2YvHaspUuaCFAagFBqAOWGnRouZXho2Hii0DUI55aEwAVuOxmg2qXNBSHfPxWAkGoA5YadGi5leGjYeKLQNQjnloSAD+rBqP7dugygU/q475eKwEA1AHrLRoUfMrw8ZDxZYBKMc8NCQA36rGY9UbVLngreqYj8dKMAB1wEqLFjW/Mmw8VGwZgHLMQ2MCcB8e21uuygVv7Yv5eKwEA1AHrLRoUfMrw8ZDxZYBKMc8NCQAD+7DH+4pV+WCgwxALTAAdcBKixY1vzJsPFRsGYByzENDArB5H/5wd7kqFzQzALXAANQBKy1a1PzKsPFQsWUAyjEPDQnApn34w6pXVLmgiQGoBQagDlhp0aLmV4aNh4otA1COeWhIADbuw2O7XlHlgkYGoBYYgDpgpUWLml8ZNh4qtgxAOeahIQH45j48tvMVVS54kwGoBQagDlhp0aLmV4aNh4otA1COeWhMANbgsR0bVbngzZqYj8dKMAB1wEqLFjW/Mmw8VGwZgHLMQ0MC8EANHtu+UZULDjAAtcAA1AErLVrU/Mqw8VCxZQDKMQ8NCcD/qMFjr29U5YL/YABqQfoAPHz4MFwuFxwOB9xuNz788ENNt+/p6UFrezsOnzxBqRDur63FvHnzsL+2Nu5joXLa2t6ueS222loqwzycy+uslQXeGjz22sqR2XcAACAASURBVEZVLvBqC0Cfz4eysjI4nU4oioLu7u7PPb63txeKoszy6tWrEcd1dnYiLy8PiYmJyMvLm3W/9+/fh9frhdPphMPhQHFxMa5cuaJp7How5wC8dOmSnuOICx0dHbDb7WhtbcXAwACqq6uRlJSE4eFh1fcxb948Simln1Er8R4vNeZ11sqCN/bjMc8mVS54Q9s7kqdOnUJDQwO6uro0BeDg4CBGR0fDTk9Ph4/x+/1ISEhAU1MTrl69iqamJtjtdvT19YWPaW5uRnJyMrq6uhAMBlFZWQmn04mJiQlt/3EekTm/ejabDYWFhThx4gSmpqb0HJNhrFy5Ert27Yq4LC8vD/X19arvw2r/r5WaXxneeaBiy3cA5ZiHRrwD+OQb+/GFbZtU+aTGAPw0WgLw7t27DzymoqICpaWlEZeVlJSgsrISwCfv/qWnp6O5uTl8fSgUQkpKCt599905j38uzDkA33vvPaxcuRI2mw1paWloaGjAr371Kz3HFlOmpqaQkJAw6wXfu3cvioqKNN2Xlb63Qs3v4ZPW/+4RFVt+B1COeWjEdwCfbNiPL2zdpMonG4wJQJfLhfT0dKxduxbnzp2LOCYjIwMtLS0Rl7W0tCAzMxMAcOPGDSiKgsuXL0ccU15ejm3bts15/HPhkd+/7e/vx7Zt2+BwOGC327Fx40b09vbqMbaYMjIyAkVRcOHChYjLGxsbkZOTE/U2oVAI4+PjYWferrXSokXNrwwbDxVbBqAc83DmdZ6YmIjYG0Oh0Jxe/0f938SjjEdNAA4ODuLIkSO4dOkS/H4/qqqqYLPZ4PP5wsfY7XacPHky4nYnT55EYmIiAODChQtQFAUjIyMRx+zYsQMlJSWq/pvohW4f4N++fRuNjY3IzMyEzWZDfn4+3n33Xdy7d0+vh9CVmQD0+/0Rlzc2NiI3Nzfqbbxeb8QXP91uNwBrLVrU/Mqw8VCxZQDKMQ9nXme32x2xN3q93jm9/lH/N1Ffiy9uflWVT9XXznk8agIwGmVlZdiwYUP439EC8MSJE3A4HAB+G4C3bt2KOGb79u3mDcDx8XEcOnQIS5cuhc1mQ1ZWFmw2GzIyMmZFlgjM5SNgvgNIzaAMGw8VWwagHPPQkHcA62rxxcpXVflUXe2cxzPXAGxsbEReXl7431J9BPzRRx9h586deOGFF5CYmIgtW7agv78/fF1hYSEKCgoeeaCxYOXKlaiqqoq4TOuPQABrLVrU/Mqw8VCxZQDKMQ+N+A7gU3W1+GLFq6qcCcC5MNcA3LRpE9asWRP+d0VFBdavXx9xTGlp6awfgRw8eDB8/dTUlLl+BNLR0YGioqLwj0C8Xi9GR0dnHXfu3DnY7fZHGmSsmDkNTFtbW8RpYG7evKnpfqy0aFHzK8PGQ8WWASjHPDQiAH9QW4s/evVVVf6gVlsATk5OIhAIIBAIQFEUtLS0IBAIhE8FV19fH/Gu3KFDh9Dd3Y1r167hypUrqK+vh6Io6OrqCh9z4cIFJCQkoLm5GVevXkVzc3PU08CkpKSgu7sbwWAQmzdvNt9pYNxuN44dO/a5p4G5ceMGXnvttbk+TMyZORF0YmIi3G53xJc51WKlRYuaXxk2Hiq2DEA55qEhAbi/Fn+06VVV/mC/tgB80ImdPR4PAMDj8aC4uDh8/MGDB5GdnQ2Hw4GFCxdi9erV6OnpmXW/nZ2dyM3Nhd1uR15eXkQgAr89EXR6ejocDgeKiooQDAa1/8d5ROYcgOfPn9dzHKbGSosWNb8ybDxUbBmAcsxDQwKwphZ/tPFVVf6gZu4fAcuI9H8KTg+stGhR8yvDxkPFlgEoxzw0IgCf3leLL73yqiqf3scA1AIDUAestGhR8yvDxkPFlgEoxzw0JgDr8KXyClU+va8u5uOxEgxAHbDSokXNrwwbDxVbBqAc89CQAKyuw5c2VKjy6WoGoBYYgDpgpUWLml8ZNh4qtgxAOeahIQG4tw5/XFahyqf3MgC1wADUASstWtT8yrDxULFlAMoxDw0JwD11+OP1Fap8eg8DUAsMQB2w0qJFza8MGw8VWwagHPPQiAB8ZncdvlxaocpndjMAtcAA1AErLVrU/Mqw8VCxZQDKMQ8NC8CSClUyALXBANQBKy1a1PzKsPFQsWUAyjEPDQnAqjp8+eUKVT5TxQDUAgNQB6y0aFHzK8PGQ8WWASjHPDQiAJ/dVYc/WVehymd3MQC1wADUASstWtT8yrDxULFlAMoxDxmA5oYBqANWWrSo+ZVh46FiywCUYx4aEYA/3FGHP11Tocof7mAAaoEBqANWWrSo+ZVh46FiywCUYx4aFoDFFapkAGqDAagDVlq0qPmVYeOhYssAlGMeGhKA2+vwp0UVqvzhdgagFhiAOmClRYuaXxk2Hiq2DEA55qERAfjc63X4yuoKVT73OgNQCwxAHbDSokXNrwwbDxVbBqAc89CQAHytDl95qUKVz73GANQCA1AHrLRoUfMrw8ZDxZYBKMc8NCQAPXX4s1UVqnzOwwDUAgNQB6y0aFHzK8PGQ8WWASjHPDQsAAsrVMkA1AYDUAestGhR8yvDxkPFlgEoxzw0JAC31eHPVlao8rltDEAtMAB1wEqLFjW/Mmw8VGwZgHLMQyMC8Edb6/DnKypU+aOtDEAtMAB1wEqLFjW/Mmw8VGwZgHLMQ0MCcEsd/nx5hSp/tIUBqAUGoA5YadGi5leGjYeKLQNQjnloSABursOfuytU+aPNDEAtMAB1wEqLFjW/Mmw8VGwZgHLMQyMC8Meb6/AXP6lQ5Y8ZgJpgAOqAlRYtan5l2Hio2DIA5ZiHhgRgZR3+oqBClT+uZABqgQGoA1ZatKj5lWHjoWLLAJRjHhoRgEpFHf5yWYUqlQoGoBYYgDpgpUWLml8ZNh4qtgxAOeahIQH4ah3+Mr9ClcqrDEAtmC4Ah4aG8Prrr8PlcsHhcCA7OxterxdTU1MRxyiKMsvTp09H3JfP54Pb7YbD4UBWVhbeeeedOY3JSosWNb8ybDxUbBmAcsxDQwJwUx3+8sUKVSqbGIBaMF0Anj59Gq+99ho++OAD3LhxA++//z7S0tJQU1MTPmYmAM+ePYvR0dGwn43EpKQkVFdXY2BgAK2trbDb7ejs7NQ8JistWtT8yrDxULFlAMoxD40IQNumOvzVixWqtDEANWG6AIzGW2+9haysrPC/ZwIwEAg88Da1tbXIy8uLuGznzp0oLCzU/PhWWrSo+ZVh46FiywCUYx4aFoB5FapkAGrDEgHY0NCA5cuXh/89E4AZGRlITU3FSy+9NOudvaKiIuzduzfisu7ubtjt9oh3Cj9NKBTC+Ph42ImJCQDWWrSo+ZVh46FiywCUYx7OvM4TExMRe2MoFJrT6x8N2yt1+OucClXaXmEAasH0AXj9+nUkJyejtbU1fNnt27fR0tKC/v5+XLx4EW+88QZsNhuOHz8ePiYnJweNjY0R93XhwgUoioJbt25FfSyv1xvxnUK32w3AWosWNb8ybDxUbBmAcszDmdfZ7XZH7I1er3dOr380bOV1+OulFaq0lTMAtSBMAH42rqJ58eLFiNuMjIxg6dKl2L59+0Pvf8+ePVi2bFn43zk5OWhqaoo45vz581AUBaOjo1Hvg+8AUjMow8ZDxZYBKMc8NOIdwIQNdfib7ApVJmzQFoA+nw9lZWVwOp1QFAXd3d2qb3v+/HkkJCSgoKAg4vKJiQlUV1cjMzMTDocDq1atQn9/f8QxHo9nVt/M5etnj4owAXj79m1cvXr1c7137174+JGREeTk5GDbtm2Ynp5+6P2fOHECDocj/O+5fAT8IKy0aFHzK8PGQ8WWASjHPDTiO4AJZXX4m6wKVSaUaQvAU6dOoaGhAV1dXZoCcHx8HNnZ2SgpKZkVgBUVFXjxxRfh8/lw/fp1eL1eJCcnY2RkJHyMx+NBaWlpxI9Ux8bGNI1dD4QJQC3MxF9lZaWq+AOAmpqaiB+KRPsRyK5du/gjEGp6Zdh4qNgyAOWYh0YEoL2sDn/rqlClXWMAfhotAVhZWYmGhgZ4vd6IALx37x4SEhLQ09MTcXxBQQEaGhrC//Z4PCgvL5/zWPXCdAE487Hv2rVrMTIyElHQMxw7dgwnT57EwMAABgcH8dZbb8Fut6OlpSV8zMxpYPbt24eBgQG0tbXxNDDUEsqw8VCxZQDKMQ8NCcD1dfjbzApV2td/EoBz+UhabQAePXoUK1aswMcffzwrACcmJsKnoPs0hYWFKC4uDv/b4/EgJSUFqampyMnJwfbt2/HrX/9a5X8R/TBdAB47duyB3xH89DF5eXlISkpCcnIyli9fHvEDkBlmTgSdmJgIl8vFE0FTSyjDxkPFlgEoxzw0IgATS+vwdxkVqkws/SQA5/KjFDUBeO3aNaSlpWFwcBAAZgUgAKxatQrFxcUYGRnB9PQ0jh8/DpvNhpycnPAxHR0d6OnpQTAYxPvvv4+CggLk5+fr+t1JNZguAEXESosWNb8ybDxUbBmAcsxDQwKwpA5/t6RClYklsXsHcHp6GitWrIh4oyhaAF6/fh1FRUVQFAUJCQlYsWIFtmzZMusrZ5/m1q1bsNvt6Orqeug49YQBqANWWrSo+ZVh46FiywCUYx4aEoAv1+HvnBWqTHw5dt8BvHv3bjjqZrTZbOHLzp07F3H85ORk+JRyFRUVWL9+/ec+/tKlS9Hc3Dzn8c8FBqAOWGnRouZXho2Hii0DUI55aEQAOl6uw9+nV6jSEcMAnJ6eRjAYjLCqqgq5ubkIBoOYnJyMeruxsTGkpKTgyJEjD7zvO3fuwOFw4L333pvz+OcCA1AHrLRoUfMrw8ZDxZYBKMc8NCQA19Xh79MqVOlYpy0AJycnEQgEEAgEoCgKWlpaEAgEMDw8DACor6/Htm3bHnj7aB8Bnz59GqdPn8bQ0BDOnDmDgoICrFy5Mnx6ucnJSdTU1MDv92NoaAi9vb1YtWoVlixZEj63sFEwAHXASosWNb8ybDxUbBmAcsxDQwJwbS3+IfVVVTrW1mq6797e3qg/KPV4PAA++bXup3+9+1miBWBHRweys7ORmJiI9PR07N69G+Pj4+Hr7927h5KSEqSmpsJutyMzMxMejyccnUbCANQBKy1a1PzKsPFQsWUAyjEPDQnANbX4h8WvqtKxRlsAyg4DUAestGhR8yvDxkPFlgEoxzw0IgCfL67FPy56VZXPFzMAtcAA1AErLVrU/Mqw8VCxZQDKMQ8NC8CFr6qSAagNBqAOWGnRouZXho2Hii0DUI55aEQAJhXV4qspr6oyqYgBqAUGoA5YadGi5leGjYeKLQNQjnloSACursVXk19VZdJqBqAWGIA6YKVFi5pfGTYeKrYMQDnmoREB+MJLtfgfL7yqyhdeYgBqgQGoA1ZatKj5lWHjoWLLAJRjHhoRgMmravG1519VZfIqBqAWGIA6YKVFi5pfGTYeKrYMQDnmoSEBWLgfX3NsUmVyYezHYyUYgDpgpUWLml8ZNh4qtgxAOeahIQG4cj++nrhJlckrGYBaYADqgJUWLWp+Zdh4qNgyAOWYh4YE4Ir9+Lp9kyqTVzAAtcAA1AErLVrU/Mqw8VCxZQDKMQ+NCMCU5fvxT7ZNqkxZzgDUAgNQB6y0aFHzK8PGQ8WWASjHPDQkAN378U/KJlWmuBmAWmAA6oCVFi1qfmXYeKjYMgDlmIdGBODCn+zHN368SZULf8IA1AIDUAestGhR8yvDxkPFlgEoxzw0JAALavCNH21U5cKCmpiPx0owAHXASosWNb8ybDxUbBmAcsxDIwJw0bIa/PNzG1W5aBkDUAsMQB2w0qJFza8MGw8VWwagHPPQsAD84UZVMgC1wQDUASstWtT8yrDxULFlAMoxD40IwMX5NfiXZzeqcnE+A1ALDEAdsNKiRc2vDBsPFVsGoBzz0JAAfHEf/uWZV1S5+MV9MR+PlWAA6oCVFi1qfmXYeKjYMgDlmIdGBGBq3j588+lXVJmaxwDUAgNQB6y0aFHzK8PGQ8WWASjHPDQkAHP34Zs/eEWVqbkMQC2YMgBdLhcURYmwrq4u4pjh4WGUlZUhKSkJixcvxp49ezA1NRVxjM/ng9vthsPhQFZWFt555505jcdKixY1vzJsPFRsGYByzEOjAvBbT72iSgagNkwbgAcOHMDo6GjYycnJ8PXT09PIz8/HmjVrcPnyZZw5cwZOpxN79uwJHzM0NISkpCRUV1djYGAAra2tsNvt6Ozs1DweKy1a1PzKsPFQsWUAyjEPDQnApfvwrQXlqkxdygDUgmkD8NChQw+8/tSpU7DZbBgZGQlf1t7eDofDgfHxcQBAbW0t8vLyIm63c+dOFBYWah6PlRYtan5l2Hio2DIA5ZiHhgRg9j586/vlqkzNZgBqwbQBmJ6ejkWLFqGgoACNjY0RH+++8cYbKCgoiLjN2NgYFEXBuXPnAABFRUXYu3dvxDHd3d2w2+2zPip+GFZatKj5lWHjoWLLAJRjHhoRgGlZ+/Cv88tVmZbFANSCKQOwpaUFPp8PH330EVpbW7F48WJs3749fP2OHTvw8ssvz7pdYmIiTp48CQDIyclBY2NjxPUXLlyAoii4detW1McNhUIYHx8POzExAcBaixY1vzJsPFRsGYByzMOZ13liYiJibwyFQnN6/aOR7qrGvz2+QZXprmrdHlcGhAlAr9c764cdn/XixYtRb9vZ2QlFUXDnzh0AnwRgSUnJrOPsdjva29sBfBKATU1NEdefP38eiqJgdHRU1RjdbjcAay1a1PzKsPFQsWUAyjEPZ15nt9sdsTd6vd45vf7RSM+sxr99b4Mq0zOrdXtcGRAmAG/fvo2rV69+rvfu3Yt625GRESiKgr6+PgCx+wiY7wBSMyjDxkPFlgEoxzw04h1AZ8ZefPu7Zap0Zux9+B2SMMIE4KPwy1/+EoqiYHh4GMBvfwTy6Y9yOzo6HvojkF27dvFHINT0yrDxULFlAMoxD434DqBzyV58+ztlqnQuYQBqwXQB6Pf70dLSgkAggKGhIXR0dMDpdKK8vDx8zMxpYNatW4fLly/j7NmzyMjIiHoamH379mFgYABtbW08DQy1hDJsPFRsGYByzEMjAnCJcw++8+31qlzi3PPwOyRhTBeAly5dQmFhIVJSUuBwOJCbmwuv14vf/OY3EcfNnAj6+eefx6JFi7Bnz55Zb0vPnAg6MTERLpeLJ4KmllCGjYeKLQNQjnloSACm78F3/m29KpekMwC1YLoAFBErLVrU/Mqw8VCxZQDKMQ8NCcC03fjOv5aqcknabk337fP5UFZWBqfTCUVR0N3drfq258+fR0JCwqzfG0xMTKC6uhqZmZlwOBxYtWoV+vv7I465f/8+vF4vnE4nHA4HiouLceXKFU1j1wMGoA5YadGi5leGjYeKLQNQjnloSACm7sa/f6tUlUtStQXgqVOn0NDQgK6uLk0BOD4+juzsbJSUlMwKwIqKCrz44ovw+Xy4fv06vF4vkpOTI/4wRXNzM5KTk9HV1YVgMIjKyko4nc7wD0uNggGoA1ZatKj5lWHjoWLLAJRjHhoWgN8sUaXWAPw0WgKwsrISDQ0N8Hq9EQF47949JCQkoKenJ+L4goICNDQ0APjk3b/09HQ0NzeHrw+FQkhJScG777475/HPBQagDlhp0aLmV4aNh4otA1COeWhEAGYs2o3v/nOJKjMWxT4Ajx49ihUrVuDjjz+eFYATExNQFAVnz56NuE1hYSGKi4sBADdu3ICiKLh8+XLEMeXl5di2bducxz8XGIA6YKVFi5pfGTYeKrYMQDnmoSEBuLAK3/3Gy6rMWFgFYG7nJVQTgNeuXUNaWhoGBwcBYFYAAsCqVatQXFyMkZERTE9P4/jx47DZbMjJyQHw27849umPhIEH/wGLWMIA1AErLVrU/Mqw8VCxZQDKMQ8NCcCUKnz3n15WZUbKJwE4l79M8rAAnJ6exooVKyLOFhItAK9fv46ioiIoioKEhASsWLECW7ZsCZ93+EF/cnb79u0MQDNipUWLml8ZNh4qtgxAOeahEQGYmbwL3/v6OlVmJu8CEJt3AO/evRuOuhltNlv4spm/MjbD5ORkOPIqKiqwfv16APwI2HJYadGi5leGjYeKLQNQjnloTADuxPe+tlaVmck75/w4at4BDAaDEVZVVSE3NxfBYBCTk5NRbzc2NoaUlBQcOXIEwG9/BHLw4MHwMVNTU/wRiFmx0qJFza8MGw8VWwagHPPQkABM2onvfXWtKjOTtAXg5OQkAoEAAoEAFEUJ/5WxmT8rW19f/7nvykX7CPj06dM4ffo0hoaGcObMGRQUFGDlypWYmpoKH9Pc3IyUlBR0d3cjGAxi8+bNPA2MWbHSokXNrwwbDxVbBqAc89CIAHQ9vwOP/+MaVbqe36Hpvnt7eyO+Kzijx+MBAHg8nvCvd6MRLQA7OjqQnZ2NxMREpKenY/fu3RgfH484ZuZE0Onp6XA4HCgqKkIwGNQ0dj1gAOqAlRYtan5l2Hio2DIA5ZiHhgSgYwce/4diVboc2gJQdhiAOmClRYuaXxk2Hiq2DEA55qExAbgdj/99kSpdju0xH4+VYADqgJUWLWp+Zdh4qNgyAOWYh0YEYFbC63jib1arMivh9ZiPx0owAHXASosWNb8ybDxUbBmAcsxDQwLQ9hqe+OuXVJlley3m47ESDEAdsNKiRc2vDBsPFVsGoBzz0JgA9OCJv1qlyiybJ+bjsRIMQB2w0qJFza8MGw8VWwagHPPQkABUPHjiLwtVmaUwALXAANQBKy1a1PzKsPFQsWUAyjEPDQnAH3nwxJ8XqjLrRwxALTAAdcBKixY1vzJsPFRsGYByzENDAvC5bZj/ZytVmfWcsX9KzewwAHXASosWNb8ybDxUbBmAcsxDYwJwK+Z/ZYUqs57bGvPxWAkGoA5YadGi5leGjYeKLQNQjnloSAA+uwXz/2S5KrOe3RLz8VgJBqAOWGnRouZXho2Hii0DUI55aEgAPrMZ87/sVmXWM5tjPh4rwQDUASstWtT8yrDxULFlAMoxDw0JwKc3Y/4f/0SVWU8zALXAANQBKy1a1PzKsPFQsWUAyjEPjQnASsz/UoEqs56ujPl4rAQDUAestGhR8yvDxkPFlgEoxzw0JACfqsD8Ly5TZdZTFTEfj5VgAOqAlRYtan5l2Hio2DIA5ZiHhgTgk5sw/wsvqjLryU0xH4+VMF0A9vb2QlGUqPb39wMAhoaGol5/+vTpiPvy+Xxwu91wOBzIysrCO++8M6cxWWnRouZXho2Hii0DUI55aEgALtiE+Y/lqTJrAQNQC6YLwKmpKYyOjka4fft2uFwu3L9/H8BvA/Ds2bMRx01NTYXvZ2hoCElJSaiursbAwABaW1tht9vR2dmpeUxWWrSo+ZVh46FiywCUYx4aEoDf34j5f5Cryqzvb4z5eKyE6QLws0xNTSEtLQ0HDhwIXzYTgIFA4IG3q62tRV5eXsRlO3fuRGFhoeYxWGnRouZXho2Hii0DUI55aEgAzn8F838/R5VZ81+J+XishOkDsLOzEzabDcPDw+HLZgIwIyMDqampeOmll2a9s1dUVIS9e/dGXNbd3Q273R7xTqEarLRoUfMrw8ZDxZYBKMc8NCQAnyjH/N9bqsqsJ8pjPh4rYfoAXL9+PdavXx9x2e3bt9HS0oL+/n5cvHgRb7zxBmw2G44fPx4+JicnB42NjRG3u3DhAhRFwa1bt6I+VigUwvj4eNiJiQkA1lq0qPmVYeOhYssAlGMezrzOExMTEXtjKBSa0+sfjazHN2D+72arMuvxDbo9rgwIE4Ber/eBP+6Y8eLFixG3+dWvfgWbzabqe3t79uzBsmXLwv/OyclBU1NTxDHnz5+HoigYHR1VNUa32w3AWosWNb8ybDxUbBmAcszDmdfZ7XZH7I1er3dOr380sr67HvP/v0xVZn13/cPvkIQRJgBv376Nq1evfq737t2LuM2BAweQmpqq6iPbEydOwOFwhP89l4+A+Q4gNYMybDxUbBmAcsxDQ94B/PdSzP9/M1SZ9e+luj2uDAgTgFq5f/8+srKyUFNTo+r4mpoaZGVlhf8d7Ucgu3bt4o9AqOmVYeOhYssAlGMeGvIdwO+UYP7/s0SVWd8pifl4rIRpA/Ds2bNQFAUDAwOzrjt27BhOnjyJgYEBDA4O4q233oLdbkdLS0v4mJnTwOzbtw8DAwNoa2vjaWCoJZRh46FiywCUYx4aEoDffhlP/N/pqsz69ssxH4+VMG0Abt68GS+99FLU644dO4a8vDwkJSUhOTkZy5cvj/gByAwzJ4JOTEyEy+XiiaCpJZRh46FiywCUYx4aEoD/ug5P/F9pqsz613UxH4+VMG0AioSVFi1qfmXYeKjYMgDlmIeGBOC31uKJ/zNVlVnfWhvz8VgJBqAOWGnRouZXho2Hii0DUI55aEgAfrMYT/wfi1SZ9c3imI/HSjAAdcBKixY1vzJsPFRsGYByzENDAvBfivHEf1uoyqx/KY75eKwEA1AHrLRoUfMrw8ZDxZYBKMc8NCIAXd9Yjcf/92RVur6xOubjsRIMQB2w0qJFza8MGw8VWwagHPPQmAB8CY//by+o0vWN6D8MJdFhAOqAlRYtan5l2Hio2DIA5ZiHhgTg11fh8f81SZWur6+K+XisBANQB6y0aFHzK8PGQ8WWASjHPDQkAL9WiMf/F4cqXV/T9occfD4fysrK4HQ6oSgKuru7Vd/2/PnzSEhIQEFBQcTlH3/8MRoaGuByueBwOJCVlYUDBw5geno6fIzH45n1p27n8kcoHhUGoA5YadGi5leGjYeKLQNQjnloRABmfnUlvvc/J6oy86srNd33qVOn0NDQgK6uLk0BOD4+juzsbJSUlMwKwMbGRixevBg9PT0YGhpCZ2cnXnjhBfz85z8PH+PxeFBaWorR0dGwY2NjmsauBwxA129frwAAHqBJREFUHbDSokXNrwwbDxVbBqAc89CYAFyO7/1PCarM/OryOT+OlgCsrKxEQ0MDvF7vrAAsKyvD66+/HnHZxo0bsXXr1vC/PR4PysvL5zxWvWAA6kDu24fw1E/3UyqEre3tmDdvHlrb2+M+FiqnuW8fkn4tlWEezvV11sL6pK3I+Iflqlyf9ElkTUxMYHx8PGwoFHro46gNwKNHj2LFihX4+OOPowZgc3MzXC4XBgcHAQCBQABpaWk4efJk+BiPx4OUlBSkpqYiJycH27dvx69//Wst/1l0gQFICCGEEEswNTWFjIyMiO/Xeb3eh95OTQBeu3YNaWlp4biLFoD3799HfX09bDYb7HY7bDYbmpqaIo7p6OhAT08PgsEg3n//fRQUFCA/P19VqOoJA5AQQohQTExMwO12Y2JiIt5DeWT4XIwlFApFvPun1zuA09PTWLFiBd55553wZdECsL29HRkZGWhvb0cwGMTx48exaNEiHDt27IH3fevWLdjtdnR1dT38CeoIA5AQQohQjI+PQ1EUjI+Px3sojwyfizl4WADevXsXiqIgISEhrM1mC1927tw5AEBGRgbefvvtiNu++eabyM3N/dzHX7p0KZqbmx/9iWiAAUgIIUQorBQafC7mQM07gMFgMMKqqirk5uYiGAxicnISALBo0aKIdwkBoKmpCTk5OQ+87zt37sDhcOC9997T58mohAFICCFEKKwUGnwu4jI5OYlAIIBAIABFUdDS0oJAIIDh4WEAQH19PbZt2/bA20f7CNjj8WDJkiXh08B0d3dj8eLFqK2tDT9mTU0N/H4/hoaG0Nvbi1WrVmHJkiWGf7TOACSEECIUoVAIXq/X8C/FxwI+F3Hp7e2ddUJmRVHg8XgAfBJzxcXFD7x9tACcmJhAdXU1MjMz4XA4kJ2djYaGBkxNTQEA7t27h5KSEqSmpsJutyMzMxMejyccnUbCACSEEEIIkQwGICGEEEKIZDAACSGEEEIkgwFICCGEECIZDEBCCCGEEMlgABJCCIk5Xq931q8t09LSwtffv38fXq8XTqcTDocDxcXFuHLlSsR9jI2NYevWrUhOTkZycjK2bt2Ku3fvGv1UVHP48GG4XC44HA643W58+OGH8R6SKlwu16zXqq6uLuKY4eFhlJWVISkpCYsXL8aePXvCv3Sdwefzwe12w+FwICsra9b58Uh8MWUAyriQEKKWnp4ezJs3Dz09PfEeCiFhvF4v8vPzMTo6Gvb27dvh65ubm5GcnIyuri4Eg0FUVlbC6XRGnButtLQU+fn58Pv98Pv9yM/PR1lZWTyezkPp6OiA3W5Ha2srBgYGUF1djaSkpLic7kMrLpcLBw4ciHitZk50DHxyUuT8/HysWbMGly9fxpkzZ+B0OrFnz57wMUNDQ0hKSkJ1dTUGBgbQ2toKu92Ozs7OeDwlEgXTBqBIC0nu24fw1E/3UyqEre3tmDdvHlrb2+M+FiqnuW8firpuf/acaTPcv38f6enpEX8KKxQKISUlBe+++y4AYGBgAIqioK+vL3yM3++HoigYHByc09odS1auXIldu3ZFXJaXl4f6+vo4jUg9LpcLhw7Nfg1nOHXqFGw2G0ZGRsKXtbe3w+FwhE8SXVtbi7y8vIjb7dy5E4WFhbEZNNGMaQNQpIXkqZ/uxxe3bKJUCA+fPIF58+bh8MkTcR8LldOnfro/6rqdlJQEp9MJl8uFyspK3LhxAwBw48YNKIqCy5cvR9ymvLw8/JcY2trakJKSMut+U1JScPToUc3rdiyZmppCQkLCrD8ttnfvXhQVFcVpVOpxuVxIT0/HokWLUFBQgMbGxoiPd994441Ze/DY2BgURQn/TdyioiLs3bs34pju7m7Y7fZZHxWT+GDaAIzHQhIKhTA+Ph525h1FBiAVSQYgjbczATgxMRFeLzs7O9He3o5gMIgzZ86guLgYaWlpuHPnDi5cuABFUSLeUQKAHTt2oKSkBADQ2NgY9e+p5uTkoKmpSeXuYQwjIyNQFAUXLlyIuPxBz0E0Wlpa4PP58NFHH6G1tRWLFy/G9u3bw9fv2LEDL7/88qzbJSYm4uTJkwA+eV0aGxsjrp95nW/duhXbJ0BUYcoAPHXqFDo7Ow1fSD773UO32w2AAUjFkgFI4+1MALrd7og10+v1htfTyclJpKWloaWl5YFhsH379oeu20uXLo34xEcEZgLQ7/dHXN7Y2Ijc3Ny4jCnad+c/68WLF6PetrOzE4qi4M6dOwAi99NPY7fb0d7eDiD6fnr+/HkoioLR0VGdnx2ZC6YMwM9i1ELCdwCpGWQA0ngb7R3A8fHxWX9Ddt26ddi1axc/AjaA27dv4+rVq5/rvXv3ot52JmhnvjbFj4CtgSUCEIjvQsIApCLJAKTxNtp3AD9LKBTCkiVLcODAgfB3tw8ePBi+fmpqKup3t/v7+8PH9PX1Cf0jkKqqqojLzPIjkM/yy1/+EoqihH/BPPMjkE+/0dLR0fHQH4Hs2rWLPwIRCEsEYLwXEgYgFUkGII230QKwpqYGPp8PQ0ND6OvrQ1lZGZKTk3Hz5k0An5y9ISUlBd3d3QgGg9i8eXPUszcUFBSEz96wbNky4U8D09bWFnEamJnnKyp+vx8tLS0IBAIYGhpCR0cHnE4nysvLw8fMnAZm3bp1uHz5Ms6ePYuMjIyop4HZt28fBgYG0NbWxtPACIYpA1C0hYQBSEWSAUjjbbQAnDkdl91uh9PpxMaNGzEwMBC+fub8renp6XA4HCgqKkIwGIy4j7GxMWzZsiV8/tYtW7YIff7WmRNBJyYmwu12w+fzxXtID+XSpUsoLCxESkoKHA4HcnNz4fV68Zvf/CbiuJkTQT///PNYtGgR9uzZM+sj/pkTQScmJsLlcvFE0IJhygAUbSFhAFKRZADSeKvmI2BCSHwxZQCKBgOQiiQDkMZbBiAh4sMA1AEGIBVJBiCNtwxAQsSHAagDDEAqkgxAGm8ZgISIDwNQBxiAVCQZgDTeMgAJER8GoA4wAKlIMgBpvGUAEiI+DEAdYABSkWQA0njLACREfBiAOsAApCLJAKTxlgFIiPgwAHWAAUhFkgFI4y0DkBDxYQDqAAOQiiQDkMZbBiAh4sMA1AEGIBVJBiCNtwxAQsSHAagDDEAqkgxAGm8ZgISIDwNQBxiAVCQZgDTeMgAJER8GoA4wAKlIMgBpvGUAEiI+DEAdYABSkWQA0njLACREfBiAOsAApCLJAKTxlgFIiPgwAHWAAUhFkgFI4y0DkBDxYQDqAAOQiiQDkMZbBiAh4sMA1AEGIBVJBiCNtwxAQsSHAagDDEAqkgxAGm8ZgISIDwNQBxiAVCQZgDTeMgAJER8GoA4wAKlIMgBpvGUAEiI+DEAdYABSkWQA0njLACREfBiAOsAApCLJAKTxlgFIiPgwAHWAAUhFkgFI4y0DkBDxYQDqAAOQiiQDkMZbBiAh4sMA1AEGIBVJBiCNtwxAQsSHAagDDEAqkgxAGm8ZgISIDwNQBxiAVCQZgDTeMgAJER8GoA4wAKlIMgBpvGUAEiI+DEAdYABSkWQA0njLACREfBiAOsAApCLJAKTxlgFIiPgwAHWAAUhFkgFI4y0DkBDxYQDqAAOQiiQDkMZbBiAh4sMA1AEGIBVJBiCNtwxAQsSHAagDDEAqkgxAGm8ZgISIDwNQBxiAVCQZgDTeMgAJER8GoA4wAKlIMgBpvGUAEiI+DEAdYABSkWQA0njLACREfBiAOsAApCLJAKTxlgFIiPgwAHWAAUhFkgFI4y0DkBDxYQDqAAOQiiQDkMZbBiAh4sMA1AEGIBVJBiCNtwxAQsSHAagDDEAqkgxAGm8ZgISIj/QBePjwYbhcLjgcDrjdbnz44Yea74MBSEWSAUjjLQOQEPGZcwAODw9jcnIy6nVTU1MYHh6e86CMoqOjA3a7Ha2trRgYGEB1dTWSkpI0j50BSEWSAUjjLQOQEPGZcwDabDYsWbIEgUBg1nV9fX2w2WyPMi5DWLlyJXbt2hVxWV5eHurr61XfR09PD1rb23H45AlKhXB/bS3mzZuH/bW1cR8LldPW9na9l2tCiM48UgBmZGQgOTkZZ86cibjODAE4NTWFhIQEdHd3R1y+d+9eFBUVRb1NKBTC+Ph42ImJCcybN49SSulnBICJiYmINTMUCsV8bSeEqOORArC3txfl5eVITEzEsWPHwteZIQBHRkagKAouXLgQcXljYyNycnKi3sbr9UJRlLBut5vvAFLh5DuANN7OvAPodrsj1kyv1xvztZ0Qoo5HCsD+/n5MT09j165dsNlsePPNNwGYKwD9fn/E5Y2NjcjNzY16m2jvAAL8DiAVy8Mn+R1AGl9nvgPIdwAJEZdHDsAZ3nzzTdhsNuzYsQPnz58XPgDn8hHwg2AAUpFkANJ4yx+BECI+ugUgABw7dgx2ux15eXnCByDwyY9AqqqqIi7T+iMQgAFIxZIBSOMtA5AQ8dE1AAHg9OnTeOGFF0wRgDOngWlra4s4DczNmzc13Q8DkIokA5DGWwYgIeITkxNBDw0NwefzxeKudWfmRNCJiYlwu91zGjcDkIokA5DGWwYgIeIj/V8C0QMGIBVJBiCNtwxAQsSHAagDDEAqkgxAGm8ZgISIDwNQBxiAVCQZgDTeMgAJER8GoA4wAKlIMgBpvGUAEiI+DEAdYABSkWQA0njLACREfBiAOsAApCLJAKTxlgFIiPgwAHWAAUhFkgFI4y0DkBDxYQDqAAOQiiQDkMZbBiAh4sMA1AEGIBVJBiCNtwxAQsSHAagDDEAqkgxAGm8ZgISIDwNQBxiAVCQZgDTeMgAJER8GoA4wAKlIMgBpvGUAEiI+DEAdYABSkWQA0njLACREfBiAOsAApCLJAKTxlgFIiPgwAHWAAUhFkgFI4y0DkBDxYQDqAAOQiiQDkMZbBiAh4sMA1AEGIBVJBiCNtwxAQsSHAagDDEAqkgxAGm8ZgISIDwNQBxiAVCQZgDTeMgAJER8GoA4wAKlIMgBpvGUAEiI+DEAdYABSkWQA0njLACREfBiAOsAApCLJAKTxlgFIiPgwAHWAAUhFkgFI4y0DkBDxYQDqAAOQiiQDkMZbBiAh4sMA1AEGIBVJBiCNtwxAQsSHAagDDEAqkgxAGm8ZgISIDwNQBxiAVCQZgDTeMgAJER8GoA4wAKlIMgBpvGUAEiI+DEAdYABSkWQA0njLACREfBiAOsAApCLJAKTxlgFIiPgwAHWAAUhFkgFI4y0DkBDxYQDqAAOQiiQDkMZbBiAh4sMA1AEGIBVJBiCNtwxAQsSHAagDDEAqkgxAGm8ZgISIDwNQBxiAVCQZgDTeMgAJER8GoA4wAKlIMgBpvGUAEiI+DEAdYABSkWQA0njLACREfBiAOsAApCLJAKTxlgFIiPgwAHWAAUhFkgFI4y0DkBDxYQDqAAOQiiQDkMZbBiAh4sMA1AEGIBVJBiCNtwxAQsSHAagDDEAqkgxAGm8ZgISIDwNQBxiAVCQZgDTeMgAJER9TBqDL5YKiKBHW1dVFHDM8PIyysjIkJSVh8eLF2LNnD6ampiKO8fl8cLvdcDgcyMrKwjvvvDOn8TAAqUgyAGm8ZQASIj6mDcADBw5gdHQ07OTkZPj66elp5OfnY82aNbh8+TLOnDkDp9OJPXv2hI8ZGhpCUlISqqurMTAwgNbWVtjtdnR2dmoeDwOQiiQDkMZbBiAh4mPaADx06NADrz916hRsNhtGRkbCl7W3t8PhcGB8fBwAUFtbi7y8vIjb7dy5E4WFhZrHwwCkIskApPGWAUiI+Jg2ANPT07Fo0SIUFBSgsbEx4uPdN954AwUFBRG3GRv7/9u726c0rj0O4P0D79ujy2ohEkQlCTIWTUxbjU8xTKb1qalGiyS3TVPaKm29mVtTjZ0Ig0zTxoexPmPxIcZadKxRK1a+90WHrRsw8aYRWPl+Zs6L7J4l58e48GWXc9iCEAJjY2MAgIqKCnR2dqr6BAIBSJIUd6v4ZRgA2dKpMQCypboxABKlP00GwJ6eHoyPj+OXX35BX18fjEYjWltblf1tbW2oqamJO06n02FwcBAAYLVa0dXVpdo/OTkJIQTW19cT/r+RSATb29tK29nZAcAAyJZejQGQLdUtFgB3dnZUr5mRSOS1vAcQ0T+XNgHQ7XbHTex4vs3OziY81u/3QwiB33//HcBfAbC2tjaunyRJ8Hq9AP4KgN3d3ar9ExMTEEJgY2PjWGMsKysDwADIll6NAZAt1S0WAMvKylSvmW63+9XeIIjotUubALi5uYmlpaUXtr29vYTHhsNhCCEwPT0N4ORuAfMKIJsWGgMgW6obrwASpb+0CYD/xKNHjyCEwNraGoC/J4EcvpXr8/leOgmkvb2dk0DYNN8YANlS3fgdQKL0p7kAODU1hZ6eHszPz2N1dRU+nw9msxkOh0PpE1sGprq6GsFgEKOjoygsLEy4DMytW7ewuLiI/v7+V14Gpvi//8G/Pv03G1tatD6vF2+88Qb6vN6Uj4UtM1vxf49epYGI0oPmAuDc3BwuX74Mg8EAWZZRXFwMt9uNP/74Q9UvthB0Tk4O8vLy4HQ6424/xBaC1ul0sFgsr7wQNBEREZGWaC4AnhY7OzsoKytTvkeoVaehDtaQPk5DHayBiLSAATBFtre3IYRQvpOoVaehDtaQPk5DHayBiLSAATBFTssL7GmogzWkj9NQB2sgIi1gAEyR0/ICexrqYA3p4zTUwRqISAsYAFMkEonA7XZrfl2s01AHa0gfp6EO1kBEWsAASERERJRhGACJiIiIMgwDIBEREVGGYQAkIiIiyjAMgEREREQZhgEwSX788UeUl5dDlmUYjUZcv35dtT/203W5ubkwGo1wOp3Y399X9Yn9dJ0syzh37lxKfrouEonAbrdDCIH5+XnVvlAohIqKCsiyDLPZjC+++ALRaFTVx+/3w2azQafTwWazIRAIJG3sq6ur+PDDD2GxWCDLMs6fPw+32x33PKd7HYl8++23Sl1lZWX4+eefUzqew7q7u3Hp0iW8+eabOHv2LBwOB5aXl1V9IpEInE4njEYjcnNzce3aNfz222+qPsc5R5Klu7sbQgi4XC5lmxZqCIfDaGpqQl5eHnJycmC32zE3N6fsj0ajcLvdMJvNkGUZlZWVWFhYUD3G1tYWmpubodfrodfr0dzcjGfPniWtBiJ6PRgAk8Dv9+PMmTPo7e3F8vIylpeX4ff7lf0HBwcoKSnBlStXEAwGMTIyArPZDKfTqfRZXV1Fbm4uXC4XFhcX0dfXB0mSVI+TDJ2dnbh69WpcANze3obJZEJDQwNCoRCGhoag1+vxzTffKH2mpqaQnZ2N7u5uLC0tobu7G5IkYXp6Oiljf/z4MW7cuIHh4WE8ffoUP/zwA0wmE27fvq2pOp7n8/kgSRL6+vqwuLgIl8uF3NxcrK2tpWQ8z6utrcXAwAAWFhYwPz+P+vp6FBUVYXd3V+nT3t6OgoICjIyMIBgM4sqVK7Db7Tg4OABwvHMkWWZmZmCxWGC321UBMN1r2NragsViwY0bNzAzM4PV1VWMjo5iZWVF6ePxeKDX6zE0NIRQKISGhgaYzWbVT8LV1dWhpKQEU1NTmJqaQklJCerr65NSAxG9PgyAJ+zPP/9EQUEB+vv7j+zz008/ISsrC+FwWNnm9Xohy7KyEOvHH38Mm82mOu7mzZu4fPnyyQz8iHHabDYsLi7GBcDe3l4YDAbVumEejwdms1m5evb++++jrq5O9Zi1tbVoaGhIyvgT+frrr3Hu3Dnl31qso7y8HO3t7aptNpsNd+7cScl4XmZzcxNCCIyPjwP4K3RLkgSfz6f0CYfDyMrKwuPHjwEc7xxJht3dXVitVoyMjKCyslIJgFqo4ZNPPsG777575P5oNIr8/Hx4PB5lWyQSgcFgwP379wFAOfcPf9iZmpqCECLuqi4RpTcGwBM2MzMDIQQePHiA0tJS5Ofno66uTnVb5bPPPoPdblcdt7W1BSEExsbGAAAVFRXo7OxU9QkEApAkKSm3kDY2NlBQUIDZ2Vmsrq7GBcCWlhY4HA7VMcFgEEIIrK6uAgAKCwvR09Oj6tPT04OioqITH/9R7t69i4sXLyr/1lod+/v7yM7OjrsF3dnZiYqKiqSP5zhWVlYghEAoFAIAjI2NQQiBra0tVT+73Q632w3geOdIMrS0tCih73AA1EINNpsNLpcL7733Hs6ePYvS0lJ89913yv6nT59CCIFgMKg6zuFwoKWlBQDQ398Pg8EQ99gGgwEPHjw40fET0evFAHjCvF4vhBAoKiqC3+/H3NwcGhsbYTQalTeLtrY21NTUxB2r0+kwODgIALBarejq6lLtn5ychBAC6+vrJ1pDNBpFXV0dvvzySwBIGABramrQ1tamOi4cDkMIgampKQCAJElKPTGDg4PQ6XQnOv6jrKysQK/Xo6+vT9mmtTpiY5ucnFRt7+rqgtVqTfp4XiYajeLatWuqK1FHPXc1NTW4efMmgOOdIyfN6/WipKQEe3t7ANQBUAs1yLIMWZZx584dBINB3L9/H7Is4/vvvwfw9+vJ4SuUsXHX1tYCOPrvymq1oru7+8RrIKLXhwHwFbndbgghXthmZ2cxODgIIYTqk3YkEoHRaFRuqxx+gT1MkiR4vV4AiV9gJyYmIITAxsbGidZw7949vPPOO8p3mY4KgLE3upjjBKeHDx9CluVXGv//W8fzY7tw4QJaW1tV21NZx6t4fmwxXV1dKC4uTvp4XqajowMWi0U1OeKo8FRdXa0KTy87R07S2toaTCaT6m/+OAEwnWqQJAlvv/22apvT6VS+RnLUB8rW1taXBsALFy6obh0TUfpjAHxFm5ubWFpaemHb29tTbg1NTEyoji8vL8enn34KIHW3gI9bg8PhQFZWFrKzs5UmhEB2drZyayiVt06PW0dMOByG1WpFS0uLEmpjeAv45DidThQWFirPY4wWbp8GAgHlb/7wORA7L0ZHR9O+hqKiorgPPL29vTCbzQB4C5go0zAAnrDt7W3IsqyaBLK/vw+TyaRcFYx9OfzwJ2+fz/fSSSDt7e1JmQSytraGUCiktOHhYQgh4Pf7lSs5vb29OHPmjCqMJpo8cfXqVdVj19XVJXXyRCz8NTQ0xIU/QDt1HFZeXo6Ojg7VtnSaBBKNRvHRRx/BbDbjyZMncfsTTaBYX19POIHiRefISdrZ2VGdA6FQCJcuXUJzczNCoZAmamhsbIybBOJyuZSrgrFJIF999ZWyf39/P+EkkJmZGaXP9PQ0J4EQaRADYBK4XC4UFBRgeHgYy8vLaG1thclkUq4WxJaHqK6uRjAYxOjoKAoLCxMuA3Pr1i0sLi6iv78/JcvAxMZy1DIwjY2NCIVCCAQCccunTE5OIjs7Gx6PB0tLS/B4PEldPiV227eqqgrhcBgbGxtK01Idz4stA9Pf369aBubXX39NyXie19HRAYPBgPHxcdVzfviqbHt7OwoLCzE6OopgMIiqqqqES6i86BxJtsO3gIH0r2FmZgaSJKGrqwsrKysYHBxEbm4uHj58qPTxeDwwGAwIBAIIhUJobGxMuAyM3W5XloF56623uAwMkQYxACbB/v4+bt++DZPJBL1ej+rq6rjFVWMLxObk5CAvLw9Op1O1FAnw90LQOp0OFoslJQtBA4kDIKBeQDk/Px+ff/55wgWUi4uLIUkSbDYbhoaGkjbugYGBI78jqKU6EoktBK3T6VBWVqYssZIOjnrOBwYGlD57e3twOp3KAsX19fVx6xge5xxJpucDoBZqePToEUpKSiDLMmw2m+q7ycDfC0Hn5+dDlmVUVFQos7Vjtra20NTUpCwE3dTUxIWgiTSIAZCIiIgowzAAEhEREWUYBkAiIiKiDMMASERERJRhGACJiIiIMgwDIBEREVGGYQAkIiIiyjAMgEREREQZhgGQiIiIKMMwABIRERFlGAZAIiIiogzDAEikUXt7eygtLcX58+exvb2tbN/Y2IDJZEJlZSUODg5SOEIiIkpXDIBEGvbkyRPo9Xpcv34dAHBwcICqqiqYTCasr6+neHRERJSuGACJNM7n80EIgXv37sHtdiMrKwsjIyOpHhYREaUxBkCiU6CjowOSJCErKwt3795N9XCIiCjNMQASnQKzs7MQQkCn02FzczPVwyEiojTHAEikcbu7u7DZbLBarTAYDHA4HKkeEhERpTkGQCKNa25uRm5uLhYWFuD3+yGEQE9PT6qHRUREaYwBkEjD+vr6IITAwMCAss3pdEKSJMzMzKRuYERElNYYAIk0KhQKIScnBx988IFqeyQSwcWLF2GxWPDs2bMUjY6IiNIZAyARERFRhmEAJCIiIsowDIBEREREGYYBkIiIiCjDMAASERERZRgGQCIiIqIMwwBIRERElGEYAImIiIgyDAMgERERUYZhACQiIiLKMAyARERERBmGAZCIiIgow/wPKmnrNZWxweAAAAAASUVORK5CYII=\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "grid.plot_3d_slicer(np.ones(grid.vnC)*model.res_x)  # x-resistivity"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 3. Source field\n",
    "\n",
    "The source is an x-directed dipole at the origin, with a 10 Hz signal of 1 A (`src` is defined either as `[x, y, z, dip, azimuth]` or `[x0, x1, y0, y1, z0, z1]`; the strength can be set via the `strength` parameter)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "sfield = emg3d.utils.get_source_field(grid, src=[0, 0, 0, 0, 0], freq=10)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 4. Calculate the electric field\n",
    "\n",
    "Now we can calculate the electric field with `emg3d`:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "* WARNING :: ``emg3d.solver.solver()`` is renamed to ``emg3d.solve()``.\n",
      "             Use the new ``emg3d.solve()``, as ``solver()`` will be\n",
      "             removed in the future.\n",
      "\n",
      ":: emg3d START :: 07:31:17 :: v0.9.3.dev17\n",
      "\n",
      "   MG-cycle       : 'F'                 sslsolver : False\n",
      "   semicoarsening : False [0]           tol       : 1e-06\n",
      "   linerelaxation : False [0]           maxit     : 50\n",
      "   nu_{i,1,c,2}   : 0, 2, 1, 2          verb      : 3\n",
      "   Original grid  :  48 x  32 x  32     => 49,152 cells\n",
      "   Coarsest grid  :   3 x   2 x   2     => 12 cells\n",
      "   Coarsest level :   4 ;   4 ;   4   \n",
      "\n",
      "   [hh:mm:ss]  rel. error                  [abs. error, last/prev]   l s\n",
      "\n",
      "       h_\n",
      "      2h_ \\                  /\n",
      "      4h_  \\          /\\    / \n",
      "      8h_   \\    /\\  /  \\  /  \n",
      "     16h_    \\/\\/  \\/    \\/   \n",
      "\n",
      "   [07:31:18]   2.623e-02  after   1 F-cycles   [1.464e-06, 0.026]   0 0\n",
      "   [07:31:18]   2.253e-03  after   2 F-cycles   [1.258e-07, 0.086]   0 0\n",
      "   [07:31:18]   3.051e-04  after   3 F-cycles   [1.704e-08, 0.135]   0 0\n",
      "   [07:31:19]   5.500e-05  after   4 F-cycles   [3.071e-09, 0.180]   0 0\n",
      "   [07:31:19]   1.170e-05  after   5 F-cycles   [6.531e-10, 0.213]   0 0\n",
      "   [07:31:19]   2.745e-06  after   6 F-cycles   [1.532e-10, 0.235]   0 0\n",
      "   [07:31:20]   6.873e-07  after   7 F-cycles   [3.837e-11, 0.250]   0 0\n",
      "\n",
      "   > CONVERGED\n",
      "   > MG cycles        : 7\n",
      "   > Final rel. error : 6.873e-07\n",
      "\n",
      ":: emg3d END   :: 07:31:20 :: runtime = 0:00:02\n",
      "\n"
     ]
    }
   ],
   "source": [
    "efield = emg3d.solver.solver(grid, model, sfield, verb=3)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The calculation requires in this case seven multigrid F-cycles and taken just a few seconds. It was able to coarsen in each dimension four times, where the input grid had 49,152 cells, and the coarsest grid had 12 cells."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 5. Plot the result\n",
    "\n",
    "We can again utilize the in-built functions of a `discretize`-grid to plot, e.g., the x-directed electric field."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\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('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",
       "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 = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(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 (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.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 = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\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 nav_element = $('<div/>');\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\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",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option);\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\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",
       "\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",
       "\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]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.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, 0, fig.canvas.width / mpl.ratio, fig.canvas.height / mpl.ratio);\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",
       "    {\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.updated_canvas_event = function() {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message(\"ack\", {});\n",
       "}\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = \"image/png\";\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src);\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\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(\"No handler for the '\" + msg_type + \"' message type: \", msg);\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(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\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",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\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",
       "    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",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\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",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, {key: value,\n",
       "                             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",
       "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\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.close = function() {\n",
       "        comm.close()\n",
       "    };\n",
       "    ws.send = function(m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function(msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['data'])\n",
       "    });\n",
       "    return ws;\n",
       "}\n",
       "\n",
       "mpl.mpl_figure_comm = function(comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = $(\"#\" + id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm)\n",
       "\n",
       "    function ondownload(figure, format) {\n",
       "        window.open(figure.imageObj.src);\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy,\n",
       "                           ondownload,\n",
       "                           element.get(0));\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.get(0);\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",
       "\n",
       "    var output_index = fig.cell_info[2]\n",
       "    var cell = fig.cell_info[0];\n",
       "\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
       "    var width = fig.canvas.width/mpl.ratio\n",
       "    fig.root.unbind('remove')\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).html('<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/mpl.ratio\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<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 () { fig.push_to_output() }, 1000);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>');\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\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) { continue; };\n",
       "\n",
       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
       "    buttongrp.append(button);\n",
       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
       "    titlebar.prepend(buttongrp);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(el){\n",
       "    var fig = this\n",
       "    el.on(\"remove\", function(){\n",
       "\tfig.close_ws(fig, {});\n",
       "    });\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
       "    // this is important to make the div 'focusable\n",
       "    el.attr('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",
       "    }\n",
       "    else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\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",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which == 13) {\n",
       "        this.canvas_div.blur();\n",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\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('matplotlib', mpl.mpl_figure_comm);\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4nOy9eVRUZ/7/yR8zZ+bMnN+cmTNzZn7f5dfZu5N0Ouk2SWfpztLpRI1JTGJyEUoEZdUqhRJIKS7gArIIGJe0Csqi2IY2aBvScQm4QoQkmpiIETeCxCRKE9so8BWin/nD1M29t+7nLsVTCFXv1znvP4DLrec+93me+6rlfiqMAAAAAABASBF2sxsAAAAAAAAGFwggAAAAAECIAQEEAAAAAAgxIIAAAAAAACEGBBAAAAAAIMSAAAIAAAAAhBgQQAAAAACAEAMCCAAAAAAQYkAAAQAAAABCDAggAAAAAECIAQEEAAAAAAgxIIAAAAAAACEGBBAAAAAAIMSAAAIAAAAAhBgQQAAAAACAEAMCCAAAAAAQYkAAAQAAAABCDAggAAAAAECIAQEEAAAAAAgxIIAAAAAAACEGBBAAAAAAIMSAAAIAAAAAhBgQQAAAAACAEAMCCAAAAAAQYkAAAQAAAABCDAggAAAAAECIAQEEAAAAAAgxIIAAAAAAACEGBBAAAAAAIMSAAAIAAAAAhBgQQAAAAACAEAMCCAAAAAAQYkAAAQAAAABCDAggAAAAAECIAQEEAAAAAAgxIIAAAAAAACEGBBAAAAAAIMSAAAIAAAAAhBgQQAAAAACAEAMCCAAAAAAQYkAAAQAAAABCDAggAAAAAECIAQEEAAAAAAgxIIAAAAAAACEGBBAAAAAAIMSAAAIAAAAAhBgQQAAAAACAEAMCCAAAAAAQYkAAAQAAAABCDAggAAAAAECIAQEEAAAAAAgxIIAAAAAAACEGBPAm0dfXR9XV1dTX13ezmzIgguE4cAxDh2A4DhzD8H98keBYwFAFAniT6O7uJkmSqLu7+2Y3ZUAEw3HgGIYOwXAcOIbh//giwbGAoQoE8CYRLBMpGI4DxzB0CIbjwDEM/8cXCY4FDFUggDeJYJlIwXAcOIahQzAcB45h+D++SHAsYKgCAbxJ9PT0kMfjoZ6enpvdlAERDMeBYxg6BMNx4BiG/+OLBMcChioQQACCjEOHDlFYWBgdOnToZjcFgJAF83DgXL+GVxoDCQRQAIvnvUNTo9f4nSkxpeaZZJykuLXGiTdI4jo2iVPK2CQ4y9nEuyrYxE1nksLErZ/YGfqZnFZpO5PSmbyx3lZiPExmDl52722gsLAw2r23Qfz+meOz209cf/tz7rhxwI0bdpwx49JoLHPj32jeGM03w3kabzLH49aarhNm68xA1jFlFs97x6+1NNdVQa5R+XKcjtWqaI9X27fxrnJVVOdSuXakVqrCrgsDmNu256GFuWU4DxTHoxrvPuP55/7xGZua/tX2v9OxWnV+cl0Vgq+mvvx4cQb1d461lB8vzgh4e4IJCKAApkavoeceWeB3nn1skXkeN86fn8g2zDNP5fD502I2f3o2l83TI/PYPDU6n82TY5i8WKCbJ8bq548v6+cPry7Rzzg+j7/G5PVCW3lMYhI+ePmgfj+FhYXRB/X7xe+fOT67/cT1N3uOuHP66hJ2HHDjhhtn3Lg0Gsvc+DeaN0bzzXCePpVjOs/N1gmzdWYg65gyU6PX+LWWukbl0+j/mC5n5INZqmiPV9u3T43KU0V5HlVrxytLVOHWiYHMbdvz0MLcMpoXyuNRjXfteFb0j7b/tP2r7f+RD2apzo9rVL7gq6kvVy+8SP/1zW2WcvXCiwFvTzABARQABBACCAGEAEIAIYAQQPH0XniBer651VJ6L7wQ8PYEExBAAUAAIYAQQAggBBACCAEUT8+FMXTlm19YSs+FMQFvTzABARQABBACCAGEAEIAIYAQQPFcvvA8XTr3Pyzl8oXniYiosLCQ3G43bd++PeDtG85AAAUAAYQAQgAhgBBACCAEUDw/XBhNF8/9p6X8cGF0wNsTTEAABQABhABCACGAEEAIIARQPBcvjKJ/nvsPS7l4YVTA2xNMQAAFAAGEAEIAIYAQQAggBFA8318YSRfO/bulfH9hZMDbE0xAAAUQl1LBX1QshLtQGS1YPguYgdxoFzIfaTFYlB6N5PPIhCI2v5/I5+EY/Tw0ST8PxjKJ08+IhGLd/C7RIElMpujnt1Ntxmk/D7j8y649Ny48u/bsN9zOnzbZPW6u/9j+Zs4Pd05HJBSz44AbN9w448al0Vjmxr/RvDGab2bibDbPzdYJs3VmIOuYMnEp/tWHczpWq2TjyefzVdEer7ZvteeHW0dGxBero1wnmPlvZT7bnYdW5pnRvFHNBcXxGI13Zf9o+0/bv9r+f/L5fNX5cTpWC76a+tJ5fiR9+/W/WUrneQigHYalADqdTpIkySelpaVERJSVleXzt6VLl6r2cfnyZVq+fDlFR0dTdHQ0LV++nK5cueJXeyCAEEAIIAQQAggBhACK58L5kXTu63+zlAtDUADPnz9P8+fPJ7fbTampqdTb23uzmyQzLAXw0qVLdPHiRTlHjhwhSZLo6NGjRHRDAFevXq3aRvvl1Tk5OZSamkqtra3U2tpKqamplJub61d7IIAQQAggBBACCAGEAIrn3HfPUVvHf7eUc989F/D22CUzM5OOHTtGRDdeePrxxx9vcot+ZlgKoJby8nKaNm0aXb9+nYhuCGB5eTm7fUdHB0mSRCdOnJB/19raSpIk0blz52w/PgQQAggBhABCACGAEEDxdHz3HJ3u+O+W0jHEBPDs2bO0cOHCm90MlmEvgP39/RQbG0s1NTXy77Kysig2NpZiY2NpxowZVFlZST09PfLf6+vrKSYmxmdfMTExtHv3btttgABCACGAEEAIIAQQAiie9u9G0omOf7OU9u/svQXc0tJCubm5lJiYSJIkUXNzs882O3bsIKfTSQ6Hgzwej/xqnhWam5spPz+f8vLyyOPxqDxlKDDsBbCxsZHGjx9PXV1d8u8++OADOnLkCLW3t1NDQwM5nU6VhdfU1FBycrLPvpKTk2nLli3sY/X19VF3d7ccr1RCACGAEEAIIAQQAhiqAtjT06O6Nvb19fl1/vVo+3YkHT/7b5bS9q09ATx8+DBt2rSJmpqadAWwsbGRIiIiqK6ujjo6Oqi8vJyioqKos7NT3sbj8VBqaqpPurq66MMPP6TJkydTZ2cn9fX10fz58+nIkSNC+kUEw14As7OzTT+7d/r0aZIkiU6fPk1EvABOnz6dtm7dyu6nurpadWOJx+MhIgggBBACCAGEAEIAQ1cAPR6P6tpYXV3t1/nX4/S3I6nl7L9byumfBNAfIdUTwIyMDCopKVH9zu1208aNGy21vbW1lbKzs+Wft23bRtu2bbP0v4PBsBbACxcuUHh4OH300UeG212/fp0iIiKosbGRiPx/CxivAEIAIYAQQAggBBACOHivAJ74dhR9fvY/LOXEtzcKQfsjpFoB7O/vp/Hjx/tIYVlZGWVmZlpq+48//khvvPEGXb58ma5du0a5ubn0ySef2Dj6wDKsBbC6upoSEhJM76ppb28nSZKopaWFiH6+CeTkyZPyNidOnMBNIBBACCAEEAIIAYQADqHPALZ+M5o+a/9PS2n95sZXwYl4BbCrq4skSaLjx4+rtuPeQeQ4fPiw/LZwRYV/8yJQDFsBvHbtGk2dOpWqqqpUv//2229p8+bNdOrUKTp//jwdOnSI3G43eTweunbtmrxdTk4Opaeny2Vg0tLS/C4DMym9kv8mCQvfJsFWmFfk0QjjPOIwjj9C9nBMET00mY/PIqqIoWzZlCpWYqbp5/5kJil8fuNmMkM/93FJZZImLr9ON87OvTcuPDv37qdfp4t9bPb4mP7g+o/rb/Yccec0mR8H7EXXrqwajGVu/BvNG6P5ZjRPfz+xyHSem60TZuuM2Vr1+OuFxmvdT5mUXunXWpoUv1b1zSfa/Wr7Qyv32vPDrSP3T9eEWQsGMre189Cbgcwtw/miOB6j8a7sH23/aftX79wqz09S/Fq/zrMdjn3zPB1q/x+Wcuyb54mIqLCwkNxuN23fvt3y43AC2NraqtqupqaGUlJSxBzcTWbYCuBnn32m+4pdZ2cnZWZm0uTJkykyMpKmTZtGZWVldPnyZdV2ly9fpmXLlsmFoJctW+Z3IWgIIAQQAggBhABCACGA4jn6zRj6qP0WSzn6zRi/HycQbwEPdYatAA4lIIAQQAggBBACCAGEAIrn829eoKavbrWUz795gYjEvAJIdOMmEO83jHmxcxPIUAcCKAAIIAQQAggBhABCACGA4vnsmxep8avbLOWzb160te/e3l5qa2ujtrY2kiSJamtrqa2tTS7z4i0DU19fryoDc+HChUAc6qADARQABBACCAGEAEIAIYAQQPEcOjeW9rfdYSmHzo0lIuuvAB49elR1t7A3K1eulLfxFoKOjIwkj8cj30waDEAABQABhABCACGAEEAIIARQPJ+ce5n2tt1lKZ+cezng7QkmIIACgABCACGAEEAIIAQQAiie5nOvUH3bLy2l+dwrROTfZwBDEQigACCAEEAIIAQQAggBhACK5+DX42jXmbst5eDX4wLenmACAiiAmJnrjYspmyzKRkVj5QKnUUWGeTjaOH6LHFewN4m/kBpKm8u+oLESxiyYrBy9wedeD5NZ+rknw2ZmM5kjPsoLj1/74No62/5xc/3H9Td7jgyElxsH3LixK5lGY5mVSYN5YzTfjObpQ5PN57nZOmG61pitVRH8OqdMzMz1fq2lCc5yenpknhxtgW1tsW9tv2vFjpU57ZMq5VowUxFuzluYz5bmoYV5ppo3M9VRtlt1TEbSqOgfbf9p+1evyLny/CQ4y8VeTHVo/Hoc7Thzj6U0QgBtAQEUAAQQAggBhABCACGAEEDxHPj6dfrHmV9byoGvXycivAVsFQigACCAEEAIIAQQAggBhACKZ19HOL17+jeWsq8jPODtCSYggAKAAEIAIYAQQAggBBACKJ76sxFUc+q3llJ/NiLg7QkmIIACgABCACGAEEAIIAQQAiieXWcd9LdTIyxl11lHwNsTTEAABQABhABCACGAEEAIIARQPDvPTqC3Tz5kKTvPTiAifAbQKhBAAUAAIYAQQAggBBACCAEUz/tnJ9LGk7+3lPfPTgx4e4IJCKAAIIAQQAggBBACCAGEAIrnvfZoWn/iEUt5rz064O0JJiCAAojxrA9oEWezQs6/n1hED00yjt+S5+TjU0yVWXi1YYsoc4WPOWmbqR9WSAyk5+65+vlVJpMsJvP188sFfmShf9m576cLz779xtv60Sbu+Nj+YPqP62/2HBlJJjMOuHHDFuJlxqXRWObGv9G88VcOR8QXm87zQBeSfjSCX+eUifH4J4BxKRX0xEsFcrTH5yN8GlnXnjtuvTBaG+6e93PYOc/NbbvzkJtbzBxStu3uefwcMZoHyv7R9p+2f/XGmPL8xKVUCL6a+vJuewxVnHjMUt5tjwl4e4IJCKAAIIAQQAggBBACCAGEAIrn719NpnWtf7CUv381OeDtCSYggAKAAEIAIYAQQAggBBACKJ4tX8VSaesfLWXLV7FEhJtArAIBFAAEEAIIAYQAQgAhgBBA8bzTFkerjz9pKe+0xQW8PcEEBFAAEEAIIAQQAggBhABCAMVT3ZZAbx1/2lKq2xIC3p5gAgIoAAggBBACCAGEAEIAIYDieftMIq348k+W8vaZxIC3J5iAAAoAAggBhABCACGAEEAIoHg2nkmiN7/8s6VsPJMU8PYEExBAAUAAIYAQQAggBBACCAEUT9WZqVR87DlLqTozNeDtCSYggAKYnFZJf3h1CZvHXzOOlUXXrMCrWQHZEQnFbIwKOvsteUxx3vtS/RA6pliqbWnj5GU+L0l3LWKSzSRHP3cuZpLL5448/7Jz/08Xnv37DbdjH5tr62L++Nj+YPqP62/2HHHn1EgmmXFjWxgNxrJdMdQrvquM0TwdkWBeKNq0ELTJOmO2Vj3+Gr/OKTM5rdL/tXTcEjna4/cRPo3EG60dKrHTPFHhJE41brn5zMwjbh5amXequWU0b5QCqTgelSQajH9t/2n7V28MKs+Pv+fZDpVnnLTk2EhLqTzjJCLcBWwVCKAAIIAQQAggBBACCAGEAIqn4rSL8ltGW0rFaVfA2xNMQAAFAAGEAEIAIYAQQAggBFA8605Po8Utz1vKutPTAt6eYAICKAAIIAQQAggBhABCACGA4ll7Kpmyj75gKWtPJQe8PcHEsBTA6upqkiRJlfj4ePnv169fp+rqakpMTCSHw0FZWVl09uxZ1T4uX75My5cvp+joaIqOjqbly5fTlStX/GoPBBACCAGEAEIAIYAQQPGsOZlCC754yVLWnEwJeHuCiWErgKmpqXTx4kU5ly5dkv++detWio6OpqamJmpvb6elS5dSYmIi9fT0yNvk5ORQamoqtba2UmtrK6WmplJubq5f7YEAQgAhgBBACCAEEAIonrdOpNLcz1+xlLdOpAa8PXY4d+4cpaeny3E4HNTc3HyzmyUzbAUwPT1d92/Xr1+nhIQE2rp1q/y7vr4+iomJoV27dhERUUdHB0mSRCdOnJC3aW1tJUmS6Ny5c7bbAwGEAEIAIYAQQAggBFA8K0+k0ewjr1rKyhNpAW+Pv/T29lJsbCz19vbe7KbIDFsBjIqKosTERHI6nbR06VL67rvviIjou+++I0mS6MyZM6r/yc/PpxUrVhARUX19PcXExPjsNyYmhnbv3s0+bl9fH3V3d8vxvqIIAYQAQgAhgBBACGCoCmBPT4/q2tjX1+fX+ddj2Yl0mnnkNUtZdkL/haGhwIEDB6i4uPhmN0PFsBTAw4cP08GDB6m9vZ2OHDlCWVlZFB8fTz/88AMdP36cJEmirq4u1f+sXr2asrOziYiopqaGkpN9PyyanJxMW7ZsYR9X+9lDj8dDRERJCevomacXs3lqVJ5h/vhygWnMFu6Ho4sMY1Rc1l8BZAs6zyimX6fzuXeWfrgiwNyF324BZVZSsg0kLJ9JQZFubi9kUqSf24r9yFLj7PjpwrNj//4bv7O5f66ttxfxx8f1B9d/XH9z54c7p9qLoOqCyIwbbpxx49JoLHPj328BNCkEbTbPzdYJs3XGbK16alSe4VrnTVLCOr/W9piZ6+mx8EI52iLaWgHXCp/Pk0Gl5BmsBUr5Uj1RUs7tJYpYmMM+89AbK/NO8Viq+aR9EqeURma+aGVX2T/a/tP2r14hc+X5iZl5o+C3x+NRXRurq6v9Ov96vNn6Br3x2euW8mbrG7b23dLSQrm5uZSYmEiSJOm+Pbtjxw5yOp3kcDjI4/HQsWPH/DqO/Px8ampq8ut/A8WwFEAtvb29FB8fT7W1tbIAfv/996ptVq1aZSqA06dPV711rIV7BRACCAGEAEIAIYAQwFAVwEC+AljcOpPSPgu3lOLWmbb2ffjwYdq0aRM1NTXpCmBjYyNFRERQXV0ddXR0UHl5OUVFRVFnZ6e8jcfjodTUVJ8oX4Tq7u6m2NhYunr16sA6QzBBIYBERAsXLqSSkpKAvgXMAQGEAEIAIYAQQAhgqApgICk8PpPcn463lMLj9gRQiZ4AZmRkUElJiep3brebNm7caGvf+/bto2XLlvndtkARFALY19dHSUlJtHnzZvkmkL///e/y3/v7+3VvAjl58qS8zYkTJ/y+CQQCCAGEAEIAIYAQQAigeAq+nEXJhyMspeDLWUTk3yuSWgHs7++n8ePH+0hhWVkZZWZm2jqG3Nxc+vjjj239z2AwLAWwsrKSWlpa6Pz583TixAnKzc2l6OhounDhAhHdKAMTExNDzc3N1N7eTm+++aZuGZj09HS5DExaWprfZWAggBBACCAEEAIIAYQAiifv2GyadshhKXnHZhORf59J1ApgV1cXSZJEx48fV23HfYSMo7u7m+Lj46m/v9/y/wwWw1IAvXX9IiIiKDExkZYsWUIdHR3y372FoBMSEsjhcFBmZia1t7er9nH58mVatmyZXAh62bJlfheChgBCACGAEEAIIAQQAiiexcfm0NRDEyxl8bE5RES0ZMkSSklJoW3btvn9CqBXAFtbW1Xb1dTUUEpKcBScHpYCONSAAEIAIYAQQAggBBACKJ7sljmU9EmUpWS3zPH7cQL5FvBQBQIoAAggBBACCAGEAIoVQD35gADqy18wC+CilnmU8HG0pSxqmUdERIWFheR2u2n79u2WH4e7CaS0tFT1O39uAhmqQAAFEDujwnBRfWKsSV4yj1nh1UcjCw3z0KQiNr9LKmbzwDQ+hgLIFNX99RvFdE+GfpTFS5XhLvBckWGuiDErc/nFrNiwovSmfm5dXqifFfq5ZaVB3lriV3Yc+OnCc2C/8bbM43JtvXUFc2zLC9n+YC92nEhywmhUnJorNs2MG26ccePSaCxz499o3hjNN6N5+tCkItN5brZOmK41ZmvVWPMnq398uYBiZ1T4tZZOnLNBVbhaW5j4Xo86PsKnKSDOrQs+818hXNw858Y+N4/YecjNOytzSzuPlKLIzBftvFD2j7b/tP2r7f/7k4tV52finA2Cr6a+LGjJpLiPYyxlQYu9V+Z6e3upra2N2traSJIkqq2tpba2NrnMi7cMTH19vaoMjPd+g+EOBFAAEEAIIAQQAggBhABCAMWTdTSLJn00yVKyjmYRkfVXAI8ePaq6WcSblStXytt4C0FHRkaSx+OhlpaWQB7uoAIBFAAEEAIIAYQAQgAhgBBA8cz7Yj5FN8dayrwv5ge8PcEEBFAAEEAIIAQQAggBhABCAMUz5/MFFNUUZylzPl8Q8PYEExBAAUAAIYAQQAggBBACCAEUz+zPF5KjKd5SZn++kIj8uwkkFIEACgACCAGEAEIAIYAQQAigeGYdWUQRBxMsZdaRRQFvTzABARQABBACCAGEAEIAIYAQQPHMPJJN4R8mWcrMI9kBb08wAQEUwISsDfTQ5CL/Y7LoPzSpiB6OMYlJfTCj/30wjo9hjcAUPka107gLLVvvj6vrl6MfVvIM6tvZFTpWqlYxWVOgm1+UGmRtvl/Z0fDThadhv/G2zONybb1lTQF/fEx/2BVG7vwYyTs3Drhxw9YH5ATQYCxz499o3hjNt0DO84djLKw1A1nHFJmQ5Z8YTMjaoGqPVqzvma2Oj/Bp6vsp6+axkldcRLctK5TDznNmPnNz1tI85OYdM7d85pGi3awYaopHK/tH23/a/tV7cqM8P/6eZzukfbaYxjVOtZS0zxYTEd4CtgoEUAAQQAggBBACCAGEAEIAxZP6aS690uC0lNRPcwPenmACAigACCAEEAIIAYQAQgAhgOJxH86lsQdcluI+DAG0AwRQABBACCAEEAIIAYQAQgDFk3Ioj17cP81SUg7lBbw9wQQEUAAQQAggBBACCAGEAEIAxTP9UD6N2TfdUqYfyicifAbQKhBAAUAAIYAQQAggBBACCAEUj+tQAY3el2wprkMFAW9PMAEBFAAEEAIIAYQAQgAhgBBA8Tg/WUIj96ZYivOTJQFvTzABARTA+Lwqw7pfQyIu8bl/Op/fzOCjrTXlzT1z9MPVc+PqwikXP9Wiz9Wqe7OIr9PHCA8rT2X5+qnI009lLptbNyz2Kzsa9v104dlnuB372FxbK/LY42P7g+k/rr9ZMWTO6e1LithxwI0bbpxx49JoLHPjPxDz7QHXEFhHLGZ8XpVfa6ljwQaVEGvrLmqlXVvjTis7qpqeBk/sVKK3+uewc5uZw3bnITvvmLmlbNstq3k5VM0dba1NpQxqawRq+lev9qXy/DgWBF4Ap3xcSM/ucVvKlI8LA96eYAICKAAIIAQQAggBhABCACGA4kn6qJCe2T3DUpI+uiGA+AygNSCAAoAAQgAhgBBACCAEEAIonoTmInq6PtVSEpqLAt6eYAICKAAIIAQQAggBhABCACGA4olvLqan6tIsJb65OODtCSYggAKAAEIAIYAQQAggBBACKJ645qX0RF26pcQ1Lw14e4IJCKAAIIAQQAggBBACCAGEAIpnctOb9IcP3rCUyU1vBrw9wQQEUAAQQAggBBACCAGEAEIAxTPp4Jv0+C6PpUw6CAG0AwRQAFJxFd2XVux3jOqMydGZiMpwFy85Mw0yy79wddP8DnNh5uq2aRcv+SLA1IszEkC2rl9JgX7K83TDLfS3/TVHN3dUL2Jz1+aFuvllzQLD7Gy8ceHZ2biPflmzgN0P97hcW2/7aw5/IWP6g+0/pr+582NYB5AZB2y9P2aciRzL/s4pw3k603yem60TZuvMQNYxZaRi/wQwMruKfpdULEfbr9qajnfmqqOVHe5Jns+TEuXYVIpeZZ6cW6sWy7l9U44cbs5q56E33BxU7lP5WMo2+Dz5UrabeYLl82RK0T/a/tP2r97YVp6fyGz/zrMdoj9cRo/unGkp0R8uC3h7ggkIoAAggBBACCAEEAIIAYQAiie6cTk9smOWpUQ3LicilIGxyrAUwC1bttCsWbNo4sSJFBcXR/n5+XTu3DnVNllZWSRJkipLl6o/IHr58mVavnw5RUdHU3R0NC1fvpyuXLliuz0QQAggBBACCAGEAEIAxRPVuIIe3p5hKVGNKwLenmBiWApgdnY27dmzh86ePUttbW2Um5tLU6dOpd7eXnmbrKwsWr16NV28eFFOd3e3aj85OTmUmppKra2t1NraSqmpqZSbm2u7PRBACCAEEAIIAYQAQgDFM6FhBT30foalTGiAANphWAqglkuXLpEkSdTS0iL/Lisri8rLy9n/6ejoIEmS6MSJE/LvWltbSZIkn1cTzYAAQgAhgBBACCAEEAIonsgDK+l3/5hjKZEHVga8PcFEUAjgt99+S5IkUXt7u/y7rKwsio2NpdjYWJoxYwZVVlZST0+P/Pf6+nqKiYnx2VdMTAzt3r3b1uNDACGAEEAIIAQQAggBFM/4/W/RA+/NtZTx+98KeHuCiWEvgNevX6e8vDyaN2+e6vcffPABHTlyhNrb26mhoYGcTictXLhQ/ntNTQ0lJyf77Ohl7/8AACAASURBVC85OZm2bNmi+1h9fX3U3d0txyuUEEAIIAQQAggBhACGqgD29PSoro19fX1+nX89wve9RffXzrWU8H1DTwBra2tpxowZ5Ha7ad26dXT9+vWb3SSZYS+ApaWl5HQ66Z///KfhdqdPnyZJkuj06dNExAvg9OnTaevWrbr7qK6uVt1U4vF4iIho3Ioq4wvCbJNwFyRF7p5nHO5iJ2c+H+2kV+aubIPk8NEuLMpoa3XJMbjIi8jtRUVstHXB5PpgnMAwNfRufztbN5ys3bctkw33LHfE+7MNs+vDvRQWFka7PtxLI96fze6He1wjueSOj60pyPQf19/c+Qn02Lgjnx+XRmOZnQMG88ZovhnN01/NN5/nZuuE6VpjtlbNtibA41b4J4ARuVX0W2exnLvnqqPtyzsKilS5bak6nPT51KtUPmlZnytH+eTnzr8tkvOrmgVyuDmrnYfecHNQuU/lYynboGzbLetzVe3m6m/61NpU9I+2/7T9q+3/u+cWq85PRO6N8+zxeFTXxurqar/Ovx7h+/5Cv3l3nqWE7/uLsMcVwaVLl2jatGl09epVunbtGs2dO5daW1tvdrNkhrUArlu3jqZMmULnz5833fb69esUERFBjY2NROTfW8DcK4AQQAggBBACCAGEAIaqAAbyFUBp7yrDdVIZae8qYY8rgkuXLtHUqVPpypUrdPXqVZo1axZ9++23N7tZMsNSAK9fv05r166lxMRE+uabbyz9T3t7u+pGEe9NICdPnpS3OXHihF83gUAAIYAQQAggBBACGKoCGEhe27Oa7v17lqW8tme1rX23tLRQbm4uJSYmkiRJ1Nzc7LPNjh07yOl0ksPhII/HQ8eOHbP1GO+//z5FR0fTpEmTaOPGjbb+N9AMSwEsLS2lmJgYamlpUZV5uXr1KhHduClk8+bNdOrUKTp//jwdOnSI3G43eTweunbtmryfnJwcSk9Pl8vApKWl+VUGBgIIAYQAQgAhgBBACKB4xu1eTfdszbKUcbvtCeDhw4dp06ZN1NTUpCuAjY2NFBERQXV1ddTR0UHl5eUUFRVFnZ2d8jYej4dSU1N90tXVRZcvX6bs7Gy6fPkyXb16lbKyslTVSm42w1IAtQWevdmzZw8REXV2dlJmZiZNnjyZIiMjadq0aVRWVkaXL19W7efy5cu0bNkyuRD0smXL/CoEDQGEAEIAIYAQQAggBFA84+pX091b5lvKuPobAujPW9J6ApiRkUElJSWq37ndbsuv5H344YdUWloq/7xt2zb6+9//bul/B4NhKYBDDQggBBACCAGEAEIAIYDiebV+japvjPJq/Roi8u+mFK0A9vf30/jx432ksKysjDIzMy21vbW1ld544w35JpDFixfTRx99ZOPoAwsEUAAQQAggBBACCAGEAEIAxfNK/RrT0lfevPKTAIp4BbCrq4skSaLjx4+rtuMqiHD89a9/JbfbTTNmzEAZmGDk1VVVxgtzlknMFn0TSTMVtWwTISvwXQi8MZImo7p63AX+1uWFdMtKJn9Zop/VTBi50C7uVnLLGpvR1uP6KZwgcbJl9MXmT9en6mb0vmTD1B28ceGpO7iXRu9LZvfDPa6RlLICyPSH3X7169xxks6NG26cMePSaCxz499o3hjNN6N5emeu+Tw3WydM1xqztSrLwpPNzGJ6dZV/YjA+r4oemFYsR7vfOxerY/ZETnl+VWOmIk8VZU1L5bi+652FcpQ3GyiF7vFdHjnP7J4hRzsPvVFuo/xf5T6Vj6Vsg3beqWpxKo5HdayacW70hEvbv3rnVnl+xucFXgBfritR9YFRXq678XZtYWEhud1u2r59u+XH4QRQW7alpqaGUlJSxBzcTQYCKAAIIAQQAggBhABCACGA4nn5gxK2mL02L39QYr5DhkC8BTzUgQAKAAIIAYQAQgAhgBBACKB4xu4qVb0lbpSxu27ccCHiFUCiGzeBKG/iILJ3E8hQBwIoAAggBBACCAGEAEIAIYDiGbuz1PArM5UZu7PUfIcKent7qa2tjdra2kiSJKqtraW2tja5zIu3DEx9fb2qDMyFCxcCcaiDDgRQABBACCAEEAIIAYQAQgDF89LOtey6o81LO9cSkfVXAI8ePapbUm7lypXyNt5C0JGRkeTxeIZUHb+BAgEUAAQQAggBhABCACGAEEDxvLRjLd2+KcdSXtqxNuDtCSYggAKAAEIAIYAQQAggBBACKJ4Xt69VlcIxyovbIYB2gAAK4JWSKtOF11DeFplHOzG1GUgNPEOR09SQ4upq+YS7KJfw8nRLZZ5ubq1arBvuWSD3AWGj+lFcZXnuK4e47dk71JiyBU/UpbOJOJigm2mHHIbZ07SbwsLCaE/Tbpp2yMHuh3tcw1ILzPGJ6j/u/Bh98JsbB9y44cYZOy6NxjIz/o3mjdF8M5qntxeZ10M0WyfM1pmBrGPKvFLinxiEF1TR/cnFcrSCqj1ebd0/H4HnnrRp6ukpx41yvCuf/Dz0foacp+rS5Iw94JIT3RwrRzsPvVFuo/xf5T6Vj6Vsg3Y+KtutOiajJ2GK/tH2n7Z/9Z4kKM9PeMEgCOD7a+m2jTmW8uL79t4CDnUggAKAAEIAIYAQQAggBBACKJ4X/rFO9UqnUV74x7qAtyeYgAAKAAIIAYQAQgAhgBBACKB4xrxX5vsNKEzGvFcW8PYEExBAAUAAIYAQQAggBBACCAEUz5jaMv6rJjUZU3tDAPEWsDUggAKAAEIAIYAQQAggBBACKJ4xteX8nNVkTG15wNsTTEAABQABhABCACGAEEAIIARQPGPeLfe5a5vLmHfLA96eYAICKAAIIAQQAggBhABCACGA4hmzrZxuKc+zlDHbyokIbwFbBQIogJdLq4wXVrMafSaL9kAFz0zy2Lp8Kwv5Wmqrl9Av1uazMXqWxgodU2OOkxFlrSxlHnhvrm6Mau4pF19luFp7rzQ4dcPJlqMpXjfKC4A2yguFMn87NcIwDc27KCwsjBqad9HfTo1g98M9LtdWR1M8e3xcf3D9x/U3d364c/rAe3PZccCNG26cscJoMJbZ8W8wb4zmm9E8HRRBNFmr7sq29oT15VI/BbCwin7jLpajFUtt3cTblhWqoq15x9X+05YPUT6hUD5RUdbmUz5JUo7xpE+i5CxueV6Odh56o9xG+b/KfSofS9kG7RMpZbuVx8PVBNTWBdT2n7Z/9eReeX7CCwdDACv4J2eajNlWEfD2BBMQQAFAACGAEEAIIAQQAggBFM+Yv1fQL9blW8qYv0MA7QABFAAEEAIIAYQAQgAhgBBA8YzZWmF4rVFmzFYIoB0ggAKAAEIAIYAQQAggBBACKJ4xWyosfy3kmC0QQDtAAAUAAYQAQgAhgBBACCAEUDxjaiqMb8RSZEzNDQHETSDWgAAKAAIIAYQAQgAhgBBACKB4xrxT4XsnM5Mx7+AVQDtAAAUAAYQAQgAhgBBACCAEUDxj3qk0nk+KjHmnMuDtCSYggAJ4eW2V8YKZY5w78sxz+5Iiw9xWbJxblxey0S6aqhi95G5Uk8ng+xq5um12Re93/5ijm8d3eXTz7B43m3GNU3Uz6aNJukn7LFw3ysVdmYoTj+lm15m72Rw/+2+6ufbtXYb55KOdFBYWRp98tJOufXsXux/ucbm2Vpx4jD0+rj+4/uP6mzs/3Dl9fJeHHQe2xZAZl4bfP8qNf6N5YzDfjObprcsLTee52Tphts6YrVV35ViTxJfXBkYAtcfjs5Zp6jGqpFzxlWG3/TVHFa72n/JJyMi9KXKimuLkLPjiJTk7ztwjRzsPvVFuo/xf5T6Vj6Vsg7Jt2rqAyuNRHqvPkxNlvUpN/2n7d0gI4OZK45qziozZDAG0AwRQABBACCAEEAIIAYQAQgDFM+ZvlcbXKEXG/A0CaIeQF8AdO3aQ0+kkh8NBHo+Hjh07Zuv/Dx06RHUNDbRz337/s988OwaaAwFIg1H22c5OLo362dW4VzcffKifuoN86g/u0c3upt262dtUr5t9zXW6aWj+QDcfNu9k09y8QzeffLTTMFVV6yksLIyqqtbTJx/x++Eel2trQ/MH7PFx/cH1H9ff3PnhzukHH+qPgV2Ne9lxw40zf8YsO/4DMd8ODHwdMF1rBrKOKVLX0GB7LT506BDV72+gXXv3y9Hu1/a6ZnFNUo4Dbi1RjkflWFbOA+U80s5Db5TbKP9XuU9u7PuMcUvj18bY1I4XnXOrPD/1++2fZ7uMqa6kW95aYiljqiGAdvBbAE+fPi2yHTeFxsZGioiIoLq6Ouro6KDy8nKKioqizs5Oy/sICwtDEARBNLHLzW4vMjjn2S5j3l5v/Dl1Rca8vT7g7Qkm/D574eHhNHv2bNq/fz/19/eLbNOgkZGRQSUlJarfud1u2rhxo+V94BVAvAKIVwDxCiBeAcQrgHgFMDC88PZ6w+/WVuaFnwQQZWCs4bcA7t27lzIyMig8PJzi4+Np06ZN9M9//lNk2wJKf38/jR8/npqbm1W/Lysro8zMTFv7wmcA8RlAfAYQnwHEZwDxGUB8BlA8L2xabzovvHlhE14BtMOAX789efIkrVixghwOB0VERNCSJUvo6NGjItoWULq6ukiSJDp+/Ljq9zU1NZScnKz7P319fdTd3S2np6eHiCCAEEAIIAQQAggBDF0B7OnpUV0b+/r6/Dr/erz41/U+5Wq4vPhXCKAdhL2Bf+nSJaqpqaGpU6dSeHg4paam0q5du+jq1auiHkIoXgFsbW1V/b6mpoZSUlJ0/6e6upokSZLj8XiICAIIAYQAQgAhgBDA0BVAj8ejujZWV1f7df71eHHjBtP6mN68uHGDsMcNBYQJYHd3N7333ns0ffp0Cg8PJ5fLReHh4TRlyhQfyRoK+PMWMF4BhABCACGAEEAIIARwEF8BrNpAty0tspQXqyCAdhiwAH711Ve0Zs0amjhxIkVGRtKyZcvo5MmT8t9mz55N6enpA25oIMjIyKDS0lLV7+zeBEIEAYQAQgAhgBBACGDoCmAgebFqg+m492YoCuC2bdtoxowZlJqaSvv27bvZzVHhtwA2NjZSZmamfBNIdXU1Xbx40We7L774giIiIgbUyEDhLQNTX1+vKgNz4cIFW/uBAEIAIYAQQAggBBACKJ6XNmww/QpEb17aMLQEsL29nTweD129epWuXr1Kc+fOpStXrtzsZskMqAyMx+OhPXv2GJaB+e677+itt97y92ECjrcQdGRkJHk8HmppabG9DwggBBACCAGEAEIAIYDieWn9Brq9sMhSXlo/tASwsbGR1q5dK/9cWlpKDX6USAoUfgvgl19+KbIdwxoIIAQQAggBhABCACGA4nmpcoPpuPbmpUp7AtjS0kK5ubmUmJhIkiT53BNANLBvC+vo6KAZM2bQlStX6MqVK+R2u+ndd9+11cZAEvJfBScCCCAEEAIIAYQAQgAhgOIZW7GB7igospSxFfYE8PDhw7Rp0yZqamrSFUAr3xbm8XgoNTXVJ11dXUREtGvXLvJ4PDR//nxasWIF/eMf/xh4pwgCAiiAl0ur6K5FxXxMFss7F5vnjnzjmH02wujWecOv11m9hI12YVEtMhV5bG6tWqyb29/O1o1dMXzgvbm6US6k2jxVl6ab0fuSdfNKg1M3EQcTdONoitfN2AMuNtHNsbr526kRhmlo3kVhYWHU0LyL/nZqBLsf7nG5tjqa4tnj4/qD6z+uv7nzw53TB96ba1/0mHHGjUujscyOf4N5YzTfzMpcmM1zs3XCbJ2xIneGa91Pebk0MAKoveBr68BphfoXpQVylOftto05qtz5t0Vy7tmaJWfE+7PlPFGXLkc5xpM+iZKjfFKknYfeKLdR/q9yn8rHUrZB2bZ7tmap2q08HtUYVfTBL0rVT0C0/aft36EhgFWm49qbsRX+35WsJ4Aivi1MyapVq+jQoUN+/W8ggAAKAAIIAYQAQgAhgBBACKB4xpZXWXqX7I68Yhpb7n9dQq0Aivq2sH/9619ERHTu3DlKS0ujH3/80fL/BhoIoAAggBBACCAEEAIIAYQAimdsWRXdmVtsKWPLxL0C6M+3hekxZ84ccrvdNGvWLDp9+rTl/xsMIIACgABCACGAEEAIIAQQAiieseuqLF0j71xcTGPX3WhPYWEhud1u2r59u+XH4QTQzreFDTcggAKAAEIAIYAQQAggBBACKJ6X11ZZuzkpx/+bj4gC9xbwUAYCKAAIIAQQAggBhABCACGA4jGtsqFz97mIVwCJxH1b2FAFAigACCAEEAIIAYQAQgAhgOIxvb4OYNz19vZSW1sbtbW1kSRJVFtbS21tbXKZF1HfFjZUgQAKAAIIAYQAQgAhgGIFUE8+7Aigj/gwAmhF/oaiAFqRP60AagXZrgAq5W+wBPCVkirdsaCXV0rsvQJ49OhR1d3C3qxcuVLeRsS3hQ1VIIACsDNA9WJlQQ2kIBpdcG5dUchGW3RVFaNCuGX5+qnM0w0rjEzhXuXCqMwvaxawuXvLfN1on3F7w21/1+aF+mFkRLnQa8PJ1rRDDsPsadpNYWFhtKdpN0075GD3wz0u19a73mGObfNCYf3HnR/unN75t0XsOGCFjhln7Lg0GsvM+DeaNzdV8EzWmYGsY3oXYruEF1TR/cnFcn41Xx3t8d62tEgVH6leU/BzlOfUoDg9VxT6ofcz5CiftHBP2rTz0BvuSZhyn8rH4go/a4uXq45JeazKPlhToH7Coek/bf9q+/9X84tV5ye8YGgKILAGBFAAEEAIIAQQAggBhABCAMXz6uoq3Xbo5dXV/n8GMBSBAAoAAggBhABCACGAEEAIoHheXV1Fv8oqthSvAAJrQAAFAAGEAEIAIYAQQAggBFA8r66qol9lFlvKq6sggHaAAAoAAggBhABCACGAEEAIoHjG/aWK7p5XbCnj/oK3gO0AARQABBACCAGEAEIAIYAQQPGMe6uK7p5bbCnj3sIrgHaAAAoAAggBhABCACGAEEAIoHjGrayie+YUW8q4lRBAO0AABWD6GQWzD69auLvJVCLNag0afIG2tvaTMn6L4/JCNmwdNE2NKjlcPTXmoqyt/2UlPoukWRhZ4GrMKRdxZYxqEz5dn6obrraeN3UH91JYWBjVHdxLo/cls/vhHpdr633bMtnjY+XJZr/6de44OePGDTfOmHFpNJb9ETmj+Wb2Zfdm89xsnTBda6x82D6An8Uan1dFD0wrlqPdr1ZotX3rs9ZYKAp9S0Ue3bphsRyuBqmyxqSyNt/juzxyntk9Q452Hnqj3Eb5v8p9cvUsfWpXKtpttfafsn+0/aftX71zqzw/4/MGSQBnF1sKBNAeEEABQAAhgBBACCAEEAIIARTPuBVVdE9GsaWMW4HPANoBAigACCAEEAIIAYQAQgAhgOJ5bXkV3Tur2FJeW45XAO0AARQABBACCAGEAEIAIYAQQPG8tqyK7p1ZbCmvLYMA2gECKAAIIAQQAggBhABCACGA4nntzSq611NsKa+9CQG0AwRQABBACCAEEAIIAYQAQgDF8/qbVfTrN4ot5XUIoC0ggAKAAEIAIYAQQAggBBACKJ7Xl1bRr9OLLeX1pbgJxA4QQAGY3qVkdvu6hfpGZhXQTRdlP+XS8IKTw8dQOPOYmNQ4G2iMLsqsrHJyUZmrG06QuPp2RrL1wHtzdaO8UOhl14c3Ljy7PtxLI96fze6He1yjeomsADL9wfUf19+sNAV4bNyRz49LQyHj5oCfkmYqaCbz3GydMF1rrJTbsHE3pl0icqvot85iOdpCv9q+1Aq0tq4dV7fU50lEed7PUdTTu+2vOXKUtSd/VbNADjdntfPQG24OKvepfCxlG7T1C5XtVh2PQU1KVd0/Tf9p+1ev2LLy/ETkBl4ApeIqui+t2FKkYrwCaAcIoAAggBBACCAEEAIIAYQAikcqqqL7UostRSqCANph2Ang+fPn6S9/+Qs5nU5yOBw0bdo0qq6upv7+ftU2kiT55NNPP1Xtq6WlhTweDzkcDnK5XLRz506/2gQBhABCACGAEEAIIARQPFJhFd03o9hSpEIIoB2GnQB++umn9NZbb9Fnn31G3333HX388ccUHx9PlZWV8jZeAfz888/p4sWLcrSSGBUVReXl5dTR0UF1dXUUERFBBw8etN0mCCAEEAIIAYQAQgAhgOIJL6yi38wotpRwCKAthp0A6rFt2zZyuVzyz14BbGtrY/9nw4YN5Ha7Vb9bs2YNzZ492/bjQwAhgBBACCAEEAIIARRPeGEV/cZdbCkQQHsEhQBu2rSJZs6cKf/sFcApU6ZQXFwczZ071+eVvczMTCorK1P9rrm5mSIiIlSvFCrp6+uj7u5uOT09PUQEAYQAQgAhgBBACGDoCmBPT4/q2tjX1+fX+dcjvKCK7k8utpTwAtwFbIdhL4DffvstRUdHU11dnfy7S5cuUW1tLZ08eZJOnTpFb7/9NoWHh9O+ffvkbZKTk6mmpka1r+PHj5MkSfT999/rPlZ1dbXqM4Uej4eIIIAQQAggBBACCAEMXQH0eDyqa2N1dbVf51+P8Pwqun96saWE5+MVQDsMGQHUypVeTp06pfqfrq4umj59Oq1atcp0/+vWraO0tDT55+TkZNqyZYtqmy+//JIkSaKLFy/q7oN7BdDObep6sVTjyKQApmmVdKOv0LH4PYvaWP2CbsthLkrcxe6uRfrhLvBGdQtveWuJfjgBVF4wFFHW5VJGuYgrc0f1IjZ3bV6oGyNB+2XNAtrZuI/CwsJoZ+M++mXNAnY/3ONybb3trzns8XH9wfYf09/c+TESN24ccOOGlR+BY9nfOWX6dVcm89y0WK7JOjOQdUxEOY7I7Cr6XVKxHG2/aoVZK+S3FxapoqpNysjgLas081xZx7IyT86tVYvl3L4pRw43Z7Xz0BtuDir3qXwsZRt86mwq281In0+NSkX/aPtP2796Y1t5fiKzA/8KoLY2pFEGoy5hMDFkBPDSpUv09ddfG+bq1avy9l1dXZScnEwrVqyga9eume5///795HA45J/9eQuYAwIIAYQAQgAhgBDAUBXAQDI+t4oecBVbyvhB+ExiMDFkBNAOXvlbunSpJfkjIqqsrFTdKKJ3E0hJSYlfN4FAACGAEEAIIAQQAggBFI/2YwFGGYybUoKJYSeA3rd9FyxYQF1dXaoyL1727NlDBw4coI6ODjp37hxt27aNIiIiqLa2Vt7GWwamoqKCOjo6qL6+3u8yMBBACCAEEAIIAYQAQgDFE7G4in47tdhSIhYPPQEsKCigSZMmUWFhoc/fPvnkE0pOTqbp06er7mMYLIadAO7Zs4f9jKByG7fbTVFRURQdHU0zZ85U3QDixVsIOjIykpxOp9+FoCGAEEAIIAQQAggBhACKJzKnin43pdhSInOGngB+8cUX9PHHH/sI4I8//kjJycnU1dVFPT09NH36dLp8+fKgtm3YCeBQBAIIAYQAQgAhgBBACKB4tGPCKIPRHn84evSojwAeP36cCgoK5J/LysrowIEDg9ouCKAAIIAQQAggBBACCAGEAIonclEV/S6x2FIiF9lrT0tLC+Xm5lJiYiJJkkTNzc0+2+zYsUP+6lmPx0PHjh2zfQx6Anjw4EFau3at/PO2bdto27Zttvc9ECCAArBzm/pNi8W7qOzEqB6T0df1cBcv7sL8ywX6uXOxfm5fUqQbIwHU1spiLxRMHTG5/pZ2gfamIk8/XP28ylxWtsyyo+HGhWdHwz7D7djH5tpaoXMB+ilsfzD9x/U3d364c3r7kiJ2HHDjhhtn3Lg0Gsvc+A/EfHvANQTWkQCX43As2EAPxhXJ0Qqs2ZNAbQ1HVga19T6VT0RW/xx2bjNz2O48ZOcdM7eUbbtltfrJk/J4OOG7vbBI1T9mT5r0nkQoz49jwQbBV1OdMbGoikYkFFuKw6YAHj58mDZt2kRNTU26AtjY2EgRERFUV1dHHR0dVF5eTlFRUdTZ2Slv4/F4KDU11SddXV3yNnoC+OGHH/oI4Lvvvmur/QMFAigACCAEEAIIAYQAQgAhgOJxLKyiEfHFluJY6P8rknoCmJGRQSUlJarfud1u2rhxo6194y3gIAYCCAGEAEIAIYAQQAigeLRjwije9vhTmForgP39/TR+/HgfKSwrK6PMzExbx6AngD/++CNNnz5ddRPIDz/8YGu/AwUCKAAIIAQQAggBhABCACGA4nHM30APxhZZimP+jfb489V0WgHs6uoiSZLo+PHjqu1qamooOTnZcvuzs7MpNjaWJkyYQElJSXTy5En5bx9//DElJyfTtGnT6IMPPrC8T1FAAAUAAYQAQgAhgBBACCAEUDwTsjbQQ5OLLGVC1o32LFmyhFJSUmjbtm1+vwLoFcDW1lbVdjU1NZSSkiL2IG8SEEABQAAhgBBACCAEEAIIARTPhKwN9NCkIkvxCqA/BPIt4KEKBFAAEEAIIAQQAggBhABCAMUTlbmBHo4pspSozBvtKSwsJLfbTdu3b7f8ONxNIKWlparf+XMTyFAFAigAOy9R68bCMxvTwR9tEoP/NfpQrdHX7tyfwseo1hhXO42r26atTeXNXTn64erF3V5UxIYVQ82FQveCoQwjPLesKdANJ06/KC2gX6zN9ys7Gvb/dOHZb7wtJ21MW29Zwwsd1x9c/7Gix5wfwzqAzDjgxg1bH5AZl0ZjmRv/RvPGaL4Fcp4/HGNhrRnIOqbzVpxfa6miPffNKFblntnq/Gq+Ondlq6OSQeWTwWJNlhXKYec5M5+5OWtpHnLzjplbPvNI0W7l8SiPVSvFyv7R9p+2f7X9f9+MYmGvuFklat4G83H/U6Lm2WtPb28vtbW1UVtbG0mSRLW1tdTW1iaXefGWgamvr1eVgblw4UIgDnXQgQAKAAIIAYQAQgAhgBBACKB4Js7dQL+fWGQpE+faewXw6NGjul8ru3LlSnkbUG+llAAAIABJREFUbyHoyMhI8ng81NLSEtDjHUwggAKAAEIAIYAQQAggBBACKJ7oORvokQlFlhI9J/DtCSYggAKAAEIAIYAQQAggBBACKJ7o2evpEUehpUTPXh/w9gQTEEABQAAhgBBACCAEEAIIARRPdMZ6ejSy0FKiM24IoD83gYQiEEABQAAhgBBACCAEEAIIARRP9Kz19GhEoaVEz8IrgHaAAAoAAggBhABCACGAEEAIoHhiZq6nx8ILLSVmJgTQDhBAAcTOqKA/vlzA5omxJnnJPH94dYlhzF4aN1r0f5dUzMao1pfe4uCNXv0ob7gL7d3z9MPVc9PWsPKGqwtnJBHaWlly7TDthcIbm8LI1hlcaRBOMk2y48BPF54D+423ZR6Xa+utK5hjMxA6rv+4/ubOD3dO71zMjwNu3HDjjBVAg7HMjX+jeWM038wEzWyem60TpmuN2Vo1ll/nlImdUeHXWjpxzgb6fVSRnPuTi1XR1mi8e646PkLIrAs+85+TQytPBpl5xM5Dbt5ZmVvaeaQUPWa++NT6U/SPtv+0/avt//uTi1XnZ+Ig3HQR41lPj0mFlhLjwVvAdoAACgACCAGEAEIAIYAQQAigeCa9sZ4ef73QUia9gVcA7QABFAAEEAIIAYQAQgAhgBBA8UxKr6THX1tiKZPSKwPenmACAigACCAEEAIIAYQAQgAhgOKZnFZJfxi3xFImp0EA7QABFAAEEAIIAYQAQgAhgBBA8UxOqzQd1954BRCfAbQGBFAAEEAIIAQQAggBhABCAMUTm1pJf3xliaXEpuIVQDtAAAUAAYQAQgAhgBBACCAEUDxm11cR4y5UgQAKIClhHT3z9GI2T43KM4yVgW1WANOsPtiI+GI2hrX+pvMxFECD2mn3ztLPPXP0w9YHZC7wbN3AbD535uqHlcaCIt1wYsPWH+QE0yhLjbNj/08Xnv37b/zO5v6N6iWy4sb0Byt0TH9z54c7p79cyI8Dtt4fM864cWk0lrnxbzRvjOab0TwdEV9sOs/N1gmzdcZsrXpqVJ7hWudNUsI6v9ZSbc233zqLVbkvVZ17Z6rjI4RZP8doLVDJobJunnJuK2TLyhz2mYfeWJl3SrFTzidNTT+V6DHzRdkHv8rSCJ+m/7T9q+3/3zqLB73uXpy7wtITkyfGFlCcGwJoBwigACCAEEAIIAQQAggBhACKJy6lwtK7ZE+8VEBxKTcEEJ8BtMawFECn00mSJKlSVVWl2qazs5Nyc3MpKiqKYmNjad26ddTf36/apqWlhTweDzkcDnK5XLRz506/2gMBhABCACGAEEAIIARQPHEpFfTkiwWW4hVAYI1hK4CbN2+mixcvyunt7ZX/fu3aNUpNTaX58+fTmTNn6MiRI5SYmEjr1v28KJ0/f56ioqKovLycOjo6qK6ujiIiIujgwYO22wMBhABCACGAEEAIIARQPHHTK+jJMfmWEjcdAmiHYSuA7733Hvv3w4cPU3h4OHV1dcm/a2hoIIfDQd3d3UREtGHDBnK73ar/W7NmDc2ePdt2eyCAEEAIIAQQAggBhACKJ25aBT35fL6lxE2DANph2ApgQkICTZ48mdLT06mmpkb19u7bb79N6enpqv+5fPkySZJEX3zxBRERZWZmUllZmWqb5uZmioiI8Hmr2AwIIAQQAggBhABCACGA4ol3VdBTo/MtJd4FAbTDsBTA2tpaamlpoa+++orq6uooNjaWVq1aJf999erVtGjRIp//i4yMpAMHDhARUXJyMtXU1Kj+fvz4cZIkib7//nvdx+3r66Pu7m45PT09RAQBhABCACGAEEAIYOgKYE9Pj+ra2NfX59f51yPBWU5Pj8yzlARnubDHDQWGjABWV1f73NihzalTp3T/9+DBgyRJEv3www9EdEMAs7OzfbaLiIighoYGIrohgFu2bFH9/csvvyRJkujixYuW2ujxeIgIAggBhABCACGAEMDQFUCPx6O6NlZXV/t1/vVImFpOTz+XZykJU8uFPW4oMGQE8NKlS/T1118b5urVq7r/29XVRZIk0YkTJ4gocG8Bc68Amn1VjdkXWJst2o9GFKqKb+rlocnGGZFQzMZfAfyNm492IVGGK6qrXYy8uWe2frSLvbzoMxd+bZFYZVhpZIoMszKZox+2kDEjQnfmFvss9laz86cLz879+w23Yx/bqOgyc3xsf3BFmpn+Zs8Rd04z+XHAjRtunLHFng3GMjf+/RZAg3k6IqHYdJ6brRNm64zZWvX4a/a+kssu2u991R6/tijxfWnqGK0dyqLfWininkCqxi03n5l5xM1DK/NONbeM5o3yCY9S8pQFzg3Gv7b/tP2rNwb1vns3kK8AJk4poz89m2spiVNuXNNRBsYaQ0YAB8Inn3xCkiRRZ2cnEf18E4jyrdzGxkbTm0BKSkr8ugkEAggBhABCACGAEMBQFcBAkphURn96JtdSEpPKzHcIZIadALa2tlJtbS21tbXR+fPnqbGxkRITEyk/P1/exlsGZuHChXTmzBn6/PPPacqUKbplYCoqKqijo4Pq6+v9LgMDAYQAQgAhgBBACCAEUDxJievomT8ttpSkRP8+ehCqDDsBPH36NM2ePZtiYmLI4XBQSkoKVVdX03/913+ptvMWgp4wYQJNnjyZ1q1b5/OytLcQdGRkJDmdTr8LQUMAIYAQQAggBBACCAEUj9ln7EV89jRUGXYCOBSBAEIAIYAQQAggBBACKJ6k+LX0zFM5lpIUvzbg7QkmIIACgABCACGAEEAIIAQQAiiepLi19OcncywlKW7oCWBBQQFNmjSJCgsLbf1tMIAACgACCAGEAEIAIYAQQAigeJLi1tKfn8i2lKEogF988QV9/PHHupJn9LfBAAIogBjPenpMKmRjtug+4jDP7ycWGeahScYxqi32uyQ+enWgvPFbDrn6gWn6sS2MGUyY+m/3zPFDJrOYcILJ1Sw0CiNJZtm576cLz779xtv60SZW0Lj+sCtt3DnizmmGH0LHjDNuXBqNZVbyDOaN0XwzqwNoNs/N1gmzdcbKk1Gjtc6bGI9/3xARl1JBT7xUIEd7fD5PULVCqK1HyqwXRmuDShS5Oc/NbbvzkJtbzBxStu3uefwcMZoHqnqVmv7T9q/eGFOen7iUwH/zxpTJa+nZP2RbypTJQ08AiYiOHj3KSp7R3wINBFAAEEAIIAQQAggBhABCAMUzZVIpPfv4IkuZMqnU1r5bWlooNzeXEhMTSZIkam5u9tlmx44d5HQ6yeFwkMfjoWPHjtk+BghgEAMBhABCACGAEEAIIARQPFNiSunZxxZZypQYewJ4+PBh2rRpEzU1NekKYGNjI0VERFBdXR11dHRQeXk5RUVFyTWHiW58C0pqaqpPurq65G0ggEEMBBACCAGEAEIAIYAQQPFMjS6h5x5daClTo0uIyL9vJtETwIyMDCopKVH9zu1208aNG20dAwQwiIEAQgAhgBBACCAEEAIonqnRa+i5RxZYytToNUTk33cTawWwv7+fxo8f7yOFZWVllJmZaesYIIBBDAQQAggBhABCACGAEEDxTI1aQ889vMBSpkbdEEARrwB2dXWRJEl0/Phx1XY1NTWUnJxsuf3Z2dkUGxtLEyZMoKSkJDp58qSlvw0GEEABQAAhgBBACCAEEAIIARSPc8JqGvnQfEtxTlhNRESFhYXkdrtp+/btlh+HE8DW1lbVdjU1NZSSkiLm4G4yEEABQAAhgBBACCAEEAIIARSP07GaRj6YZSlOx2q/HyeQbwEPVSCAAoiZuZ4eCy9kY7qoRprHrMDrw9HGMSwS7a8cGhS0fcDFR7voyEnRz29m6IctNJ3OhBGCX79RTPd6mMzSj5GQ6IYpSmwkpf5m596fLjx79/u3D66ts+0fN9d/XH+z54g7p+n8OODGDTfOuHFpNJa58e+v5JkVejab56aFoM3WGisCaLDWeRMz0z8BTHCW09Mj8+Q8MqFIlQfj1DErXM8+8TR4kql6IsHNeQvz2dI8tDDPVPPG4EkO92TG6EmLtv+0/avt/0cmFKnOT4KzXOzFVAenYxWNHJFpKU7HKiIS8wog0Y2bQEpL1XcW+3MTyFAFAigACCAEEAIIAYQAQgAhgOJxjf8LjXpgnqW4xv/F1r57e3upra2N2traSJIkqq2tpba2NrnMi7cMTH19vaoMzIULFwJxqIMOBFAAEEAIIAQQAggBhABCAMXjCn+LRt0/11Jc4W8RkfVXAI8ePaq6W9iblStXytt4C0FHRkaSx+OhlpaWgB7vYAIBFAAEEAIIAYQAQgAhgBBA8bjCV9Ko38yxFFf4SvMdAhkIoAAggBBACCAEEAIIAYQAisclraRR9822FJcEAbQDBFAAEEAIIAQQAggBhABCAMXjem0ljbp3tqW4XrshgP7cBBKKQAAFAAGEAEIAIYAQQAggBFA8rnEraPQ9GZbiGrci4O0JJiCAAoAAQgAhgBBACCAEEAIoHte45TT67lmW4hq3nIjwCqBVIIACgABCACGAEEAIoFgB1JMPCKC+/AW1AL6yjEb/cqaluF5ZFvD2BBMQQAFMSq+kx19bwuf1QsNYqaxvtigP5JtEHo7h47c4Jhpkin5sy+Q0/dgVzPtTfBdJOZx8cmFkhP32CT9iJEO/TldfeH6dLvax2eNj+oPrP/abNWzK2f3J/DiwLW3MuDQay/6InNF8C/Q3eZitM2Zr1eOvFxqvdT9lUnqlX2tpUvxaeuapHDna/Wr748FYdbTnh1tHfL69hVkLBjK3tfNQfsIygLllOF8Ux2M03pX9o+0/bf/qnVvl+UmKXyv4auqL6+U3afRdHktxvfxmwNsTTEAABQABhABCACGAEEAIIARQPK6xb9LoO9+wFNdYCKAdIIACgABCACGAEEAIIAQQAige19ilNPqOdEtxjV1KRPgMoFUggAKAAEIAIYAQQAggBBACKB7Xi8U0+rY0S3G9WBzw9gQTEEABQAAhgBBACCAEEAIIARSP64VCGn3rDEtxvVAY8PYEE8NOALnv7pMkiU6ePElEROfPn9f9+6effqraV0tLC3k8HnI4HORyuWjnzp1+tQkCCAGEAEIAIYAQQAigeFxjCmn0LW5LcY2BANph2Algf38/Xbx4UZVVq1aR0+mk69evE9HPAvj555+rtuvv75f3c/78eYqKiqLy8nLq6Oiguro6ioiIoIMHD9puEwQQAggBhABCACGAEEDxuJ5fQqP/R4qluJ5fEvD2BBPDTgC19Pf3U3x8PG3evFn+nVcA29ra2P/bsGEDud1u1e/WrFlDs2fPtt0GCCAEEAIIAYQAQgAhgOJxjS6g0f+ZbCmu0QVEhJtArDLsBfDgwYMUHh5OnZ2d8u+8AjhlyhSKi4ujuXPn+ryyl5mZSWVlZarfNTc3U0REhOqVQivEpVTQky8W+J0nxprnj68sMcwfxhnHUECNilgbFI3VKxLqjT/C+dAk/WgXKTlx+hmRUKwbQynlCvfalFU2Tvthxdcku/bcuPDs2rPfcDt/2mT3uFmp4vqbEy3mnI5I8C1eK4cZN9w480fKuPFvNG+M5puZfJnNc7N1wmydGcg6pkxcSoVf67nTsZpGPpgl58nn81XRHq9PwXzN+eHWER9pV64TzPy3Mp/tzkMr88xo3qjmguJ4jMa7sn+0/aftX23/P/l8vur8OB2r/TrPdnCNyqfR/zHdUlyj8gPenmBi2Avg4sWLafHixarfXbp0iWpra+nkyZN06tQpevvttyk8PJz27dsnb5OcnEw1NTWq/zt+/DhJkkTff/+97mP19fVRd3e3nJ6eHiKCAEIAIYAQQAggBDB0BbCnp0d1bezr6/Pr/OvhGplHo/99mqW4RuYJe9xQYMgIYHV1NXtzhzenTp1S/c8///lPCg8Pt/S5vXXr1lFaWpr8c3JyMm3ZskW1zZdffkmSJNHFixcttdHj8RARBBACCAGEAEIAIYChK4Aej0d1bayurvbr/OvhenYxjf7/plqK69nF5jsEMkNGAC9dukRff/21Ya5evar6n82bN1NcXJylt2z3799PDodD/tmft4DxCiAEEAIIAYQAQgAhgIP4CuCfc2j0/zvFUlx/zhH2uKHAkBFAu1y/fp1cLhdVVlr7sHFlZSW5XC75Z72bQEpKSvy6CQQCCAGEAEIAIYAQwFAVwEDieiabRv8/SZbieiY74O0JJoatAH7++eckSRJ1dHT4/G3Pnj104MAB6ujooHPnztG2bdsoIiKCamtr5W28ZWAqKiqoo6OD6uvr/S4DAwGEAEIAIYAQQAggBFA8rj8tolH/d4KluP60iIhwF7BVhq0AvvnmmzR37lzdv+3Zs4fcbjdFRUVRdHQ0zZw5U3UDiBdvIejIyEhyOp1+F4KGAEIAIYAQQAggBBACKB7XUwtp1P8VbymupxYGvD3BxLAVwKEEBBACCAGEAEIAIYAQQPG4nlxAo/7POEtxPbkg4O0JJiCAAoAAQgAhgBBACCAEEAIoHtcTWTTq/5hsKa4nsgLenmACAiiAqdFr6LlHFvidZx9bZJ7HjfPnJ7INo6ze7pM/LWbzp2dz2Tw9Mo/NU6Pz2Tw5holNQf7jy/r5w6tL9GMkyH5+i4uPTHPfsmBw0RedD+pvXHg+qN8vfv8D+QYJZZj+Zs8Rd05fXcKOA9uSw4xLo7HMjX+jeWM03wzn6VM5pvPcbJ0wW2cGso4pMzV6jV9rqbbor1I2Rj6Y5XO82r59alSeKsrzaPSEmlsnBjK3bc9DC3PLaF5wou8znhX9o+0/bf9q+3/kg1mDXnjZ9ccsGvXfJlmK649ZAW9PMAEBFAAEEAIIAYQAQgAhgBBA8Tgfn0cj//doS3E+Pi/g7QkmIIACgABCACGAEEAIIAQQAige5+NzaeT/NtFSnI/r3xh6MykoKKBJkyZRYWGh6vednZ2UlZVFbreb0tLS6MMPPxz0tkEABQABhABCACGAEEAIIARQPM5H59DI/zXKUpyPzgl4e+zyxRdf0Mcff+wjgN9//z21tbUREdG//vUvSkpKot7e3kFtGwRQABBACCAEEAIIAYQAQgDF43xkNo38XxyW4nzE/hc5DAZHjx71EUAtaWlp1NnZOUgtugEEUAAQQAggBBACCAGEAEIAxTP14Qx67n+OtJSpD2fY2ndLSwvl5uZSYmIiSZJEzc3NPtvs2LGDnE4nORwO8ng8dOzYMdvHYCaAp06dotTUVNv7HSgQQAFAACGAEEAIIAQQAggBFM/Uh2fSc//TeEuZ+vBMW/s+fPgwbdq0iZqamnQFsLGxkSIiIqiuro46OjqovLycoqKiVK/UeTweSk1N9UlXV5e8jZEA/vDDD+R2u+n48eO22i4CCKAAFs97h6ZGr/E7U2JKzTPJOElxa40Tb5DEdWwSp5SxSXCWs4l3VbCJm84khYlbP7Ez9DM5rdJ2JqUzeWO9rcR4mMwcvOze20BhYWG0e2+D+P0zx2e3n7j+9ufcceOAGzfsOGPGpdFY5sa/0bwxmm+G8zTeZI7HrTVdJ8zWmYGsY8osnveOX2tprquCXKPy5Tgdq1XRHq+2b+Nd5aqozqVy7UitVIVdFwYwt23PQwtzy3AeKI5HNd59xvPP/eMzNjX9q+1/p2O16vzkuvyr92iHxVHLacqDMy1lcdRyIiLq6emh7u5uOX19faaPoyeAGRkZVFJSovqd2+2mjRs32joGTgD7+vooMzNT95vKBgMIIAAAAACCgv7+fpoyZQpJkiSnurra9P+0Atjf30/jx4/3kcKysjLKzMy01SY9Abx+/TotXbrUUtsCBQQQAADAkKKnp4c8Hg/19PTc7KYMGBzL4NLX16d69c/fVwC7urpIkiSft2ZramooOTnZcnuys7MpNjaWJkyYQElJSXTy5EkiIvryyy8pPDyc0tPT5bS3t1verwgggAAAAIYU3d3dJEkSdXd33+ymDBgcy/CAE8DW1lbVdjU1NZSSkjLYzQsIEEAAAABDimASDRzL8CCQbwEPVSCAAAAAhhTBJBo4luEBdxNIaWmp6nf+3AQyVIEAAgAAGFL09fVRdXW1pc9uDXVwLEOX3t5eamtro7a2NpIkiWpra6mtrU0u8+ItA1NfX68qA3PhwoWb3HIxQAABAAAAEHIcPXpUdbewNytXrpS38RaCjoyMJI/HQy0tLTexxWKBAAIAAAAAhBgQQAAAAACAEAMCCAAAAAAQYkAAAQAAAABCDAggAACAgFNdXe3zYfv4+Hj579evX6fq6mpKTEwkh8NBWVlZdPbsWdU+Ll++TMuXL6fo6GiKjo6m5cuX05UrVwb7UCzjvYHA4XCQx+OhY8eO3ewmWcLpdPqcq6qqKtU2nZ2dlJubS1FRURQbG0vr1q2j/v5+1TYtLS3k8XjI4XCQy+WinTt3DuZhABOGpQCG4kICgFUOHTpEYWFhdOjQoZvdFABkqqurKTU1lS5evCjn0qVL8t+3bt1K0dHR1NTURO3t7bR06VJKTExUfe1YTk4OpaamUmtrK7W2tlJqairl5ubejMMxxVtCpK6uTlVCxFtiZCjjdDpp8+bNqnPV29sr//3atWuUmppK8+fPpzNnztCRI0coMTGR1q1bJ29z/vx5ioqKovLycuro6KC6ujqKiIiggwcP3oxDAjoMWwEcSgtJ5pvvUYxnPZ+ZxomeZSEZJpmzwTAT5/KJmsdnQhYfx3yDLDDIoirdRHLJ0U/EYia5+hlvkPB8mylgsoRJ4eClfn8DhYWFUf3+BvH7546P6w+b/cqdH+6cRuTy44AbN9w448al4Vhmxr/RvDGab0bzdOJc4zkePWeD+Tphss6YrVUxMw3WOUUy33xPd91OT0/XXUOvX79OCQkJtHXrVvl3fX19FBMTQ7t27SIioo6ODpIkiU6cOCFv09raSpIk0blz5/xauwNJRkYGlZSUqH43XIoIO51Oeu8933Po5fDhwxQeHk5dXV3y7xoaGsjhcMhFojds2EBut1v1f2vWrKHZs2cHptHANsNWAIfSQhLjWU+PSYV8wo3zaISFRBrnkQlFhvn9RD4PR/N5aBKfB2MNEsdnREKxbn73/7d3rl9RXHm/5984b55EZ3KdJHM5cyYzc94861lnzVrnxZwnE1dh0yJEBESMQAi214DcRIJykRgQgmiCicQAyTCZBQiIXIIk5oIOEUQxBI0SBlmGgBFGv+eFdg8t9dvV1RTQBd/PWt8Xdld3167atfnYvfevooXE6OdXG4XE6ueXm+Q8t9lk4oTE6+fZhIVLQ0srgoKC0NDSavl7S+0Tj4fJ4yqdH+mc/ipW7gdSv5H6mdQvVX1Z6v+q60Z1vamu09+tVV/jz6/ZbzhOGI0zRmPVH4IV49yMhLve1h23Q0NDER0djdjYWOTl5eHGjRsAgBs3bkDTNAwMDHi9Jjs7G4WFhQCApqYmhIeHzx5/w8PR3NxsetyeT+x+G7HY2FhERUVh3bp1SEpKQlVVldfPu8ePH5/1N3h8fByapuH8+fMAgOTkZBw+fNhrm66uLjgcjlk/FZPFwbYCuBgDydTUFCYmJjxxf6NIAaQAUgApgBTA2QI4OTnpGS87OzvR3t6OwcFBdHd3IyUlBZGRkfjhhx/Q29sLTdO8vlECgOLiYmRkZAAAqqqqEBcXN2tcjouLQ3V1tY9/PRaG0dFRaJqG3t5er8elNgQatbW16OnpwTfffIPGxkZERESgqKjI83xxcTHS09NnvS4kJARtbW0A7p+Xqqoqr+fd5/nmzZvz2wDiE7YUwC+++AKdnZ0LPpA8PPfQ5XIBoABSACmAFEAKoJ4AulwurzGzsrLSM57evn0bkZGRqK2tFcWgqKjIcNzevHmz1y8+gYBbAPv6+rwer6qqQnx8/KLsk97c+Ydz6dIl3dd2dnZC0zT88MMPALz/ns7E4XCgvb0dgP7f0wsXLkDTNIyNjVncOuIPthTAh1mogUT6BtCxp0KWEV+i+uPmlpeX1TH84yr94Y7PxbOvyHkmUZFX5fxii5yntwrZpp+ndgjZqZ8ndwl5Tc4TyUJShOzWz+OpQtIWLvWn7wtg/elW699faJ90PMTjJxxv8RxJ53SX3A+kfiP1M6lfqvqyeA0orhvV9aa6Tp+L90GgDcYJw7FmLuPYjDj23F8xOvMbwImJiVn3kE1LS0NJSQl/Al4Abt26hatXrypz584d3de6hdY9bYo/AS8NloQAAos7kFAAKYAUQAogBXC2AKqYmprChg0bcOLECc/c7Q8//NDz/PT0tO7c7f7+fs82Fy9eDOhFIKWlpV6P2WURyMOcPXsWmqZ5VjC7F4HM/KKlo6PDcBFISUkJF4EEEEtCABd7IKEAUgApgBRACqBaAI8ePYqenh4MDw/j4sWLyMrKQlhYGL7//nsA96s3hIeHo6urC4ODg8jPz9et3pCUlOSp3vDqq68GfBmYpqYmrzIw7vYGKn19faitrcWVK1cwPDyMjo4OREdHIzs727ONuwxMWloaBgYGcO7cOcTExOiWgTly5AiGhobQ1NTEMjABhi0FMNAGEgogBZACSAGkAKoF0F2Oy+FwIDo6Gjk5ORgaGvI8767fGhUVBafTieTkZAwODnq9x/j4OAoKCjz1WwsKCgK6fqu7EHRISAhcLhd6enoWe5cMuXz5Mnbs2IHw8HA4nU7Ex8ejsrISP/30k9d27kLQa9aswbp161BWVjbrJ353IeiQkBDExsayEHSAYUsBDLSBJHhfhXryvGLAf/aVXDzjS1QilpiLXySp87RLEekPokq+FAL21E4D2RKkwKw8PZYuJENIppyfZ+nnZ3uFZJvM6/tNZ2WOf6lrvS+Ada2tyu382SfT7RaOn3S8xXMkndMMuR+YllihX6r6sln5VAroNoPr1GV8nRuNE0bjjNFY9ewrvi0WCt5n/BMwIWRxsaUABhoUQAogBZACSAGkABJiJyiAFkABpABSACmAFEBXR70YAAAgAElEQVQKICF2ggJoARRACiAFkAJIAaQAEmInKIAWQAGkAFIAKYAUQAogIXaCAmgBWm6FekWswaCtWmVouHLWYAWtLyKnWmEpruJUCZvBqldR3KySM3+Ear9+VuQKyROSL6Rgn+k8esC/1LU9EMC2VuV2/uyT2D7peAjHTzrefsmqVZJpViT9kMknUgxWNKv+w2Ukj4qVzL6saP7FFh/GqiR5nJsZLZcCSEigQwG0AAogBZACSAGkAFIACbETFEALoABSACmAFEAKIAWQEDtBAbQACiAFkAJIAaQAUgAJsRMUQAtYlV9hOHnb30UYHoHbbhDFggx/F2U8keKnyClkTil0wkIC09JmVs7y95sXrUL9PPKGkIM5+nnTjxSpM1MAHyny4/2lfT2YI7ZPOh5mj6tpwczzQyaFfmZaGP2URqUcKq5T5aITdwzGCcOxZi7j2IysyqcAEhLoUAAtgAJIAaQAUgApgBRAQuwEBdACKIAUQAogBZACSAEkxE5QAC2AAkgBpABSACmAFEBC7AQF0AJWFVTMbYGGkdzNUfAMJU+xmMNfyRMXaKgWaezTj7gQw6zMSZJSqBA0SbSKhRx6XTf/UepH3sr2K3XtDwSwvVW9rR/7JLVPPB7S8ROOtz8iKUqjtABF6GdmF5n8bK95MXwsXX29qa7TgBBEHxaaPL01F6sKKICEBDoUQAugAFIAKYAUQAogBZAQO0EBtAAKIAWQAkgBpABSAAmxExRAC7CFAO6Wo5znp7gH68/3yFHdH9a06JmskyfOw1PNnTMrbmXZunnksJAje/Vz1J9kKVPXfvqBAJ5+8JjJ95f29chesX3S8TAtktL5kc7pGwoxlGoZmhVDRV+W+r/qulHOD1Rcp0/spgASQqyDAmgBFEAKIAWQAkgBpAASYicogBZAAaQAUgApgBRACiAhdoICaAEUQAogBZACSAGkABJiJyiAFrDqQIVvA+dSE0ChbtrPs+SFHj97XVG/z+SiDtOiV/K6GGmhxCPle/UjydPbWbp5tGKPblYcy5Tzrn+ZKYDKbYXPlfb10Yo9YvvE4yEcP/F4S+dIJe9mxVDoZ1K/VPVlsXbgUhZAH7LqAAWQkECHAmgBFEAKIAWQAkgBpAASYicogBZAAaQAUgApgBRACiAhdoICaAFGAmgodzuM8+QudZ5INshSFECpdp80n09RE08WPf15dqLQCbK18niGbn7+frqYx06k6ecDdeofCGB9++n7jwnvI32utK8rj2eI7ROFUZqraFYMpTmDh16XawouQQE0us6NxgnDscaH/4xSAAlZGlAALYACSAGkAFIAKYAUQELsBAXQAiiAFEAKIAWQAkgBJMROUAAtgAJIAaQAUgApgBRAQuwEBdACXiysmJvgGU3c9mWRh8E9RBf6fr8rc/aLMX1vX+n+vW/m6EcSQKlWXVm2uIhBEpuV72XqRpKqx6tSdfNUTYqYpz/UzzMfJSvT0NGCoKAgNHS04JmPksX3kT5X2tfHq1JlaRSOhyyG+sdbPD8qART6gdRvTN87WNGXTd8jeA73CV6QewX78J9RXyTxxUIKICGBDgXQAiiAFEAKIAWQAkgBJMROUAAtgAJIAaQAUgApgBRAQuwEBdACXnyjQi1vBvNyjAbtuQreYkieNJ9qUev9SbdpO5wtzlWTauVJIvREVapuJFn79cc7xfz279t183zdNmVOfnJfAE9+0oLn67aJ7yN9rkoupfZJx0M6fuLcQOn8KGo4LlZ9wJX7zYvhXORwIQTRcKza5dt/WF98gwJISKBDAbQACiAFkAJIAaQAUgAJsRMUQAugAFIAKYAUQAogBZAQO0EBtIC/FFWoyzMYDMqGpR8MSrUYlWuZj5ItKpFTyZxfQifN6SvVjygRR/bKEW5xZlYApXl1kmz9scEl5j+bEnXzp5Z4ZRo77wtgY2cL/tQSL76P9LkqKZXaZ1oApVvKSedHIe9iP5DmDFokjP5Koz9lZXwpL2NUYsanMjNGY1WKYpybkb8UUQAJCXQogBZAAaQAUgApgBRACiAhdoICaAEUQAogBZACSAGkABJiJyiAFvDCoQr13ByDQdmXKOcNGczV87ssS66fIvfGPjxSnCPGtLhJJVre2aMbs3fjUN2Rw+zdOCRBkubh/VfzK2Je7Niom9Az65VpPtOMoKAgNJ9pRuiZ9eL7SJ8r7etv/75dbJ94ZxKL7kIindMV72aK/UDqN6ZFUtGX/ZJGxfWmnFdrMLfQcH5hpjVjkdE8xMdTc/HCIQogIYEOBdACKIAUQAogBZACSAEkxE5QAC2AAkgBpABSACmAFEBC7AQF0AJeKKlQy5tBaQejQXuugrcyZz9W5MlRSt7BHDmKuzNYetcNQQrMSthztbvE/ObvO3QjlVr5ff1WU5GE6v+cShAT9VmYblLP/1mZ012NCAoKwumuRqSe/7P4PtLnqgTQbLul4ycdb+n8qO6YIvUDqd9YdXcS5R1KFNeNUg4V1+mKvPkXRMOxKt03SXyhhAJISKBDAbQACiAFkAJIAaQAUgAJsRMUQAugAFIAKYAUQAogBZAQO0EBtAAKIAWQAkgBpABSAAmxExRAC/jvtyrU9bkMBl2j2l9GtfqM6vWt3L8fKwr2iVEu5rBY8h45ulesDfezynTdSLcfe/avr+lGkpf/fXKLGEmG/l/bJt04z0TqZv1n4brZ+Pka3Ujbr/8sHGFdEbqpG3hKmU+66hEUFIRPuupRN/CU+D5m91W1v9LxkI6fdLyl86OSUqkfSP1G6mfybevkvmxaDA+9rrzeVNfpioJ9hte50ThhWGfQh/+MGtUifCwjF//9FgWQkEBn2QtgXV0dYmNj4XQ64XK58PXXX5t+DwogBZACSAGkAFIACbETfgvgyMgIbt++rfvc9PQ0RkZG/N6phaKjowMOhwONjY0YGhpCeXk5QkNDTe87BZACSAGkAFIAKYCE2Am/BTA4OBgbNmzAlStXZj138eJFBAcHz2W/FoTt27ejpKTE67GEhAQcO3bM5/f4/PPP0djejvrTrf6n1Th1c03bPKRdldOmUy+lQz8NHS26OfmJfho75TR1ntJN85lm3bScadLN6a5G3bR3ndTNJ131Yrq66nRz9tN6ZSoq3kZQUBAqKt7G2U/l95E+V9rX9q6TYvuk4yEdP+l4S+dHOqcnP9HvAw0dLWK/kfqZP31W7P/zcb21zX0cMBxr5jKOzUhje7vVwzUhxGLmJIAxMTEICwtDd3e313N2EMDp6WmsXr0aXV1dXo8fPnwYycnJuq+ZmprCxMSEJ5OTkwgKCmIYhmEeCgBMTk56jZlTU1PzPrYTQnxjTgL4j3/8A9nZ2QgJCcGpU6c8z9lBAEdHR6FpGnp7e70er6qqQlxcnO5rKisroWmaJy6Xi98A8htAfgPIbwD5DaDwDaDL5fIaMysrK+d9bCeE+MacBLC/vx93795FSUkJgoOD8cEHHwCwlwD29fV5PV5VVYX4+Hjd1+h9AwjMfQ6gUekGn8q8KG4vZXhLtzdz5JS8Lqd8r5y3s8SsfC9TN1I5j6c/TNHNrz/eqZs/Nrh0oyq5It0u7aVPX9LNq18F62ZPz//VzZGLf9BNw8CTYnq//R+6uXv9MWXOfnp/DuDZT+tx9/pj4vtInyvt65GLfxDbJx0P6fhJx1s6P9I5/WODS+wHUr8Ry8YI/VLVl8X+r7puFNeb6jo1uo3cQtxKzuwcQH4DSEjgMmcBdPPBBx8gODgYxcXFuHDhQsALoD8/AUtQACmAFEAKIAWQi0AIsROWCSAAnDp1Cg6HAwkJCQEvgMD9RSClpaVej5ldBAJQACmAFEAKIAWQAkiIvbBUAAHgyy+/xNq1a20hgO4yME1NTV5lYL7//ntT78MyMCwDwzIwLAPDMjAUQELsxLwUgh4eHkZPT898vLXluAtBh4SEwOVy+bXfFEAKIAWQAkgBpAASYieW/Z1ArIACSAGkAFIAKYAUQELsBAXQAl4oqVDeGN3wBus+DLpzniOouMG8cn6g4qb2VsvhoxV7dLPyeIZupLlcT9Wk6Oa52l1ifvP3Hbp5vm6bbn5fv9VUJHlRzUuM+ixMN6nn/6zM6a5GBAUF4XRXI1LP/1l8H+lzVbJltt3S8ZOOt3R+pHP6VI1iTp/Qb6R+ZlbylKKnuG6Uc/wU1+mKvAWY42c0VqXL49zMvFBCASQk0KEAWgAFkAJIAaQAUgApgITYCQqgBVAAKYAUQAogBZACSIidoABaAAWQAkgBpABSACmAhNgJCqAFvHCoAo+n5srxYcA0ynwKot/1A1WLR4pzxPxH6eu6eeRwtn4kYXxnj25WvJupG0kIVh7PwM/fT9fNYyfS9GNSPiWh+q/mV8RItfJCz6xXpvlMM4KCgtB8phmhZ9aL7yN9rkoAzUqYdPyk4y2dH+mcrng3U+wH4uINoZ+J/VLRl6X+728tv/kWPCvGIuVY9yAvHKIAEhLoUAAtgAJIAaQAUgApgBRAQuwEBdACKIAUQAogBZACSAEkxE5QAC3gL0UVeCI5V06KQXYbx1AQDepy+VtjUFlbMF+OX+JYlKMfaZ6hWZE8sleOcKsvqTacJDCSIJm9bd0fG1z4z6ZE3fypJV6Zxs4WBAUFobGzBX9qiRffx+yt1X798U6xfdLxEGvrSbdWk86PdE5V4ibNz5P6mT8yJ/R/f2v1GdbpM7jOjcYJw7HGaKxKUYxzM/KXIgogIYEOBdACKIAUQAogBZACSAEkxE5QAC2AAkgBpABSACmAFEBC7AQF0AJefKMCT+3MFfPkLoO8ZhyjQdloTo6/ZWj8Lj3jhzhaJowlr+tHIRGPHM3SjVkBlG4/9sxHybpRyZY0D08qreLOyU/uC+DJT1rwfN028X2kz5X29ZmPksX2mRZA4XiL50c6pyXWCZ0/MufPXL25lGExus6NxgnDscZorNolj3Mz8+IbFEBCAh0KoAVQACmAFEAKIAWQAkiInaAAWgAFkAJIAaQAUgApgITYCQqgBbxYWIGntufK2WEQHwbU+RTEhZZDsQyGJIaF+3TzyJs5+lmIu5O8l6kbSYQer0rVjeoOF09/qB+VoD3zUTIaOu4LYENHC575KFl8H+lzpX19XCF60vGw6q4bqrvOSP1A6jfi3FU/SrNYLXnzLXiGY43RWLVDMc7NyIuFFEBCAh0KoAVQACmAFEAKIAWQAkiInaAAWgAFkAJIAaQAUgApgITYCQqgBaw6UIGnt+WKMRwwfRh0jeblGJZm8LPEzEKWljE9N/Bgjn4kAXwrW8wj5Xv1I8xVk8TG7F1IJKHy5y4k7tS3n0ZQUBDq209bejcO1R05ZNET5voJx1s8PyoBFPqBVXP9/Cnb4m+pFsMyLQbXueH8PQvkTjXWubPqAAWQkECHAmgBFEAKIAWQAkgBpAASYicogBZAAaQAUgApgBRACiAhdoICaAEUQAogBZACSAGkABJiJyiAFmAkgEbxZdCd8yKRpSiAVtUHLJHlUBZDIcItzkRhFOrkrTimL1q+pO6BANa1n1ZvK3yutK+PVuyRb+EmHQ+zome21t8C1PsLJAGc8yIPC+SOAkjI0oACaAEUQAogBZACSAGkABJiJyiAFkABpABSACmAFEAKICF2ggJoAasKKvD01lw5S1UA98j5Wbaclfv260asD1iwTzeWiWGRYt5gqRCpXp10K7Mje/UjiZMywry6B5kpgPcfM/n+0r4e2Su2T6yxKBw/cT6fSclTip7Qb8R6f0K/VPVlqf8vaQFUjXUPsqqAAkhIoEMBtAAKIAWQAkgBpABSAAmxExRAC6AAUgApgBRACiAFkBA7QQG0AFsIYADdJk6cG2hWDE3OGZRuDfZooaKkjCQkxULMiqQqirI1qtS1tz4QwFb1tn7skyhu0vGQjp9UukU6R9I5VczpMy160lw/RV+2+nZv836rNwogIeQBFEALoABSACmAFEAKIAWQEDtBAbQACiAFkAJIAaQAUgAJsRMUQAugAFIAKYAUQAogBZAQO0EBtIBV+RV42pXrf3wYdOdTEJVyqFog4qccirUDpUUjOfv1I9UTzBMiiYJKGk3KpLhYQRLMN/2IajFLUQ7q2h4IYFvr/cfMvr+0rwdz5IUYfoibKZmTzmme3A+kfiMu6JBq+in6sih5qoUegSx4cxnHZmRVPgWQkECHAmgBFEAKIAWQAkgBpAASYicogBZAAaQAUgApgBRACiAhdoICaAEUQAogBZACSAGkABJiJyiAFqDlVuCZV3PF/CLJIFuMYzjx2kggd8h5cpccfxeP+CWOZoXR7CITSSRVMiktQDErmVJRYkVMS+mDzBRApWz5sU+mBU1aiGFW2hT345X6gWmh80Pm/FnMobreVNfpUzt8EDiDccJwrDEaq5LkcW5mtFwKICGBDgXQAiiAFEAKIAWQAkgBJMROUAAtgAJIAaQAUgApgBRAQuwEBdACKIAUQAogBZACSAEkxE5QAC0geF8Fnk3IlfOKOs/4kkR1jAZtfxehKP8gWbzoxKxMiiKZIcSfhSmSZAoLCcSoBMYfYVWkrvWBALa2KrfzZ59Mt9sqOZPOaYYf4mZS2vxalOGvxBksrjCUM4NxwmicMRqrnn1FMc7NSPA+CiAhgQ4F0AIogBRACiAFkAJIASTETlAALYACSAGkAFIAKYAUQELsBAXQAiiAFEAKIAWQAkgBJMRO2FIAY2NjoWmaVyoqvAeckZERZGVlITQ0FBERESgrK8P09LTXNj09PXC5XHA6ndi0aRPq6+v92h/Hngr8amOu/4k1zi9fVue5zQaJl6OUU9UfE9XCF38WtJiVUEk+pUn2KilNFiJN5t+tH38WxVid+tP3BbD+dKv17y/Jk3A8xOMnHG/xHKkWTpiUMLMLKFR9WbwGFNeN6npTXafPxRtf50bjhOFYM5dxbEYceyiAhAQ6thXAEydOYGxszJPbt297nr979y4SExOxe/duDAwMoLu7G9HR0SgrK/NsMzw8jNDQUJSXl2NoaAiNjY1wOBzo7Ow0vT8UQAogBZACSAGkABJiJ2wrgH/729/E57/44gsEBwdjdHTU81h7ezucTicmJiYAAO+88w4SEhK8Xnfo0CHs2LHD9P5QACmAFEAKIAWQAkiInbCtAEZFRWHdunVISkpCVVWV18+7x48fR1JSktdrxsfHoWkazp8/DwBITk7G4cOHvbbp6uqCw+GY9VOxERRACiAFkAJIAaQAEmInbCmAtbW16OnpwTfffIPGxkZERESgqKjI83xxcTHS09NnvS4kJARtbW0AgLi4OFRVVXk939vbC03TcPPmTd3PnZqawsTEhCeTk5MAKIAUQAogBZACqCeAk5OTXmPm1NSUJX8DCCFzJ2AEsLKyctbCjodz6dIl3dd2dnZC0zT88MMPAO4LYEZGxqztHA4H2tvbAdwXwOrqaq/nL1y4AE3TMDY25tM+ulwuAEC46238QdsnJ1id3zt8SIg6z6/Zr8zv1sr5X2FyfvuSnP8Zoch6Ob+JytXNr6OFxOjHrFD/cpMcQ4F+OHFCJMn2YeWkVWlouS+ADS2tlr+3KCbS8TB5XKXz44+0SP1G6mdSv1T1Zan/q64b1fWmuk5/t1Z9jT+/Zr/hOGE0zhiNVX8IVoxzMxLuehsA4HK5vMbMyspK//5AEEIsJ2AE8NatW7h69aoyd+7c0X3t6OgoNE3DxYsXAczfT8DSN4AUQAogBZACSAGcLYD8BpCQwCVgBHAunD17FpqmYWRkBMC/F4HM/Cm3o6PDcBFISUmJX4tAKIAUQAogBZACOFsACSGBi+0EsK+vD7W1tbhy5QqGh4fR0dGB6OhoZGdne7Zxl4FJS0vDwMAAzp07h5iYGN0yMEeOHMHQ0BCampr8LgOTnP83hLvelrNVnbBtPmS7QXa+o8zaXXJCX5OzJkWOc7ciqYqkV+gmREqmfhx7hGTpZ7Uiwdkm87qQHCH7Fi5Nre0ICgpCU2u79e8vtU86HiaPq3R+pHPqyJL7gdRvpH4m9UtlXxb6v+q6UV1vqut07S71NR628x3jccJgnDEaq8K3Ksa5GUnOl6s0EEICA9sJ4OXLl7Fjxw6Eh4fD6XQiPj4elZWV+Omnn7y2cxeCXrNmDdatW4eysrJZPz+4C0GHhIQgNjbW70LQhBBCCCF2wnYCuFSYnJyEy+XyzCO0K0uhHWxD4LAU2sE2EELsAAVwkZiYmICmaZ45iXZlKbSDbQgclkI72AZCiB2gAC4SS2WAXQrtYBsCh6XQDraBEGIHKICLxFIZYJdCO9iGwGEptINtIITYAQrgIjE1NYXKykrb18VaCu1gGwKHpdAOtoEQYgcogIQQQgghywwKICGEEELIMoMCSAghhBCyzKAAEkIIIYQsMyiAhBBCCCHLDArgAvH5559j+/btcDqdiIiIQE5Ojtfz7lvXhYaGIiIiAmVlZZienvbaxn3rOqfTiU2bNi3KreumpqaQlJQETdNw5coVr+cGBweRnJwMp9OJ6OhonDhxAvfu3fPaprOzEwkJCQgJCUFCQgK6uroWbN+Hh4fx5ptvIjY2Fk6nEy+//DIqKytnHedAb4cedXV1nna5XC58/fXXi7o/M6mursa2bduwdu1arF+/HtnZ2bh27ZrXNlNTUygrK0NERARCQ0Oxd+9e/POf//TaxpdrZKGorq6GpmkoLy/3PGaHNoyOjqKgoADr1q3DmjVrkJSUhMuXL3uev3fvHiorKxEdHQ2n04mUlBR8++23Xu8xPj6OAwcOICwsDGFhYThw4AB+/PHHBWsDIcQaKIALQGdnJ1566SXU19fj2rVruHbtGjo7Oz3P3717F4mJidi9ezcGBgbQ3d2N6OholJWVebYZHh5GaGgoysvLMTQ0hMbGRjgcDq/3WQgOHz6MPXv2zBLAiYkJREZGIi8vD4ODgzhz5gzCwsLw17/+1bNNX18fVq9ejerqaly9ehXV1dVwOBy4ePHiguz7l19+iYMHD+Krr77CjRs38NlnnyEyMhJHjx61VTsepqOjAw6HA42NjRgaGkJ5eTlCQ0MxMjKyKPvzMBkZGTh16hS+/fZbXLlyBVlZWdi4cSNu377t2aakpAQbNmxAd3c3BgYGsHv3biQlJeHu3bsAfLtGFor+/n7ExsYiKSnJSwADvQ3j4+OIjY3FwYMH0d/fj+HhYZw7dw7Xr1/3bFNTU4OwsDCcOXMGg4ODyMvLQ3R0tNct4TIzM5GYmIi+vj709fUhMTERWVlZC9IGQoh1UADnmX/961/YsGEDmpqaxG2++OILBAcHY3R01PNYe3s7nE6npxDrO++8g4SEBK/XHTp0CDt27JifHRf2MyEhAUNDQ7MEsL6+HuHh4V51w2pqahAdHe359iw3NxeZmZle75mRkYG8vLwF2X89PvroI2zatMnzbzu2Y/v27SgpKfF6LCEhAceOHVuU/THi1q1b0DQNPT09AO5Lt8PhQEdHh2eb0dFRBAcH48svvwTg2zWyENy+fRtxcXHo7u5GSkqKRwDt0IaKigq89tpr4vP37t1DVFQUampqPI9NTU0hPDwcDQ0NAOC59mf+Z6evrw+aps36VpcQEthQAOeZ/v5+aJqG5uZmbNmyBVFRUcjMzPT6WeX48eNISkryet34+Dg0TcP58+cBAMnJyTh8+LDXNl1dXXA4HAvyE9LY2Bg2bNiAS5cuYXh4eJYAFhYWIjs72+s1AwMD0DQNw8PDAICYmBjU1tZ6bVNbW4uNGzfO+/5LvPfee9i6davn33Zrx/T0NFavXj3rJ+jDhw8jOTl5wffHF65fvw5N0zA4OAgAOH/+PDRNw/j4uNd2SUlJqKysBODbNbIQFBYWeqRvpgDaoQ0JCQkoLy/H/v37sX79emzZsgUnT570PH/jxg1omoaBgQGv12VnZ6OwsBAA0NTUhPDw8FnvHR4ejubm5nndf0KItVAA55n29nZomoaNGzeis7MTly9fRn5+PiIiIjx/LIqLi5Genj7rtSEhIWhrawMAxMXFoaqqyuv53t5eaJqGmzdvzmsb7t27h8zMTHzwwQcAoCuA6enpKC4u9nrd6OgoNE1DX18fAMDhcHja46atrQ0hISHzuv8S169fR1hYGBobGz2P2a0d7n3r7e31eryqqgpxcXELvj9G3Lt3D3v37vX6Jko6dunp6Th06BAA366R+aa9vR2JiYm4c+cOAG8BtEMbnE4nnE4njh07hoGBATQ0NMDpdKKlpQXAv8eTmd9Quvc7IyMDgNyv4uLiUF1dPe9tIIRYBwXQTyorK6FpmjKXLl1CW1sbNE3z+p/21NQUIiIiPD+rzBxgZ+JwONDe3g5Af4C9cOECNE3D2NjYvLbh448/xq5duzxzmSQBdP+hc+OLOLW2tsLpdPq1/2bb8fC+bd68GUVFRV6PL2Y7/OHhfXNTVVWF+Pj4Bd8fI0pLSxEbG+u1OEKSp7S0NC95MrpG5pORkRFERkZ69XlfBDCQ2uBwOLBz506vx8rKyjzTSKT/UBYVFRkK4ObNm71+OiaEBD4UQD+5desWrl69qsydO3c8Pw1duHDB6/Xbt2/Hu+++C2DxfgL2tQ3Z2dkIDg7G6tWrPdE0DatXr/b8NLSYP5362g43o6OjiIuLQ2FhoUdq3fAn4PmjrKwMMTExnuPoxg4/n3Z1dXn6/MxrwH1dnDt3LuDbsHHjxln/4amvr0d0dDQA/gRMyHKDAjjPTExMwOl0ei0CmZ6eRmRkpOdbQffk8Jn/8+7o6DBcBFJSUrIgi0BGRkYwODjoyVdffQVN09DZ2en5Jqe+vh4vvfSSl4zqLZ7Ys2eP13tnZmYu6OIJt/zl5eXNkj/APu2Yyfbt21FaWur1WCAtArl37x7eeustREdH47vvvpv1vN4Cips3b+ouoFBdI/PJ5OSk1zUwODiIbdu24cCBAxgcHLRFG/Lz82ctAikvL/d8K+heBJIgpjsAAAPSSURBVPLhhx96np+entZdBNLf3+/Z5uLFi1wEQogNoQAuAOXl5diwYQO++uorXLt2DUVFRYiMjPR8W+AuD5GWloaBgQGcO3cOMTExumVgjhw5gqGhITQ1NS1KGRj3vkhlYPLz8zE4OIiurq5Z5VN6e3uxevVq1NTU4OrVq6ipqVnQ8inun31TU1MxOjqKsbExT+zUjodxl4FpamryKgPz/fffL8r+PExpaSnCw8PR09PjdcxnfitbUlKCmJgYnDt3DgMDA0hNTdUtoaK6RhaamT8BA4Hfhv7+fjgcDlRVVeH69etoa2tDaGgoWltbPdvU1NQgPDwcXV1dGBwcRH5+vm4ZmKSkJE8ZmFdffZVlYAixIRTABWB6ehpHjx5FZGQkwsLCkJaWNqu4qrtA7Jo1a7Bu3TqUlZV5lSIB/l0IOiQkBLGxsYtSCBrQF0DAu4ByVFQU3n//fd0CyvHx8XA4HEhISMCZM2cWbL9PnTolzhG0Uzv0cBeCDgkJgcvl8pRYCQSkY37q1CnPNnfu3EFZWZmnQHFWVtasOoa+XCMLycMCaIc2nD17FomJiXA6nUhISPCamwz8uxB0VFQUnE4nkpOTPau13YyPj6OgoMBTCLqgoICFoAmxIRRAQgghhJBlBgWQEEIIIWSZQQEkhBBCCFlmUAAJIYQQQpYZFEBCCCGEkGUGBZAQQgghZJlBASSEEEIIWWZQAAkhhBBClhkUQEIIIYSQZQYFkBBCCCFkmUEBJIQQQghZZlAACbEpd+7cwZYtW/Dyyy9jYmLC8/jY2BgiIyORkpKCu3fvLuIeEkIICVQogITYmO+++w5hYWHIyckBANy9exepqamIjIzEzZs3F3nvCCGEBCoUQEJsTkdHBzRNw8cff4zKykoEBweju7t7sXeLEEJIAEMBJGQJUFpaCofDgeDgYLz33nuLvTuEEEICHAogIUuAS5cuQdM0hISE4NatW4u9O4QQQgIcCiAhNuf27dtISEhAXFwcwsPDkZ2dvdi7RAghJMChABJicw4cOIDQ0FB8++236OzshKZpqK2tXezdIoQQEsBQAAmxMY2NjdA0DadOnfI8VlZWBofDgf7+/sXbMUIIIQENBZAQmzI4OIg1a9bgjTfe8Hp8amoKW7duRWxsLH788cdF2jtCCCGBDAWQEEIIIWSZQQEkhBBCCFlmUAAJIYQQQpYZFEBCCCGEkGUGBZAQQgghZJlBASSEEEIIWWZQAAkhhBBClhkUQEIIIYSQZQYFkBBCCCFkmUEBJIQQQghZZlAACSGEEEKWGRRAQgghhJBlxv8HtNLlTsuDqzQAAAAASUVORK5CYII=\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "grid.plot_3d_slicer(efield.fx.ravel('F'), view='abs', vType='Ex', pcolorOpts={'norm': LogNorm()})"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table style='border: 3px solid #ddd;'>\n",
       "  <tr>\n",
       "     <td style='text-align: center; font-weight: bold; font-size: 1.2em; border: 2px solid #fff;' colspan='6'>Sat Jan 11 07:31:20 2020 CET</td>\n",
       "  </tr>\n",
       "  <tr>\n",
       "    <td style='text-align: right; background-color: #ccc; border: 2px solid #fff;'>Linux</td>\n",
       "    <td style='text-align: left; border: 2px solid #fff;'>OS</td>\n",
       "    <td style='text-align: right; background-color: #ccc; border: 2px solid #fff;'>4</td>\n",
       "    <td style='text-align: left; border: 2px solid #fff;'>CPU(s)</td>\n",
       "    <td style='text-align: right; background-color: #ccc; border: 2px solid #fff;'>x86_64</td>\n",
       "    <td style='text-align: left; border: 2px solid #fff;'>Machine</td>\n",
       "  </tr>\n",
       "  <tr>\n",
       "    <td style='text-align: right; background-color: #ccc; border: 2px solid #fff;'>64bit</td>\n",
       "    <td style='text-align: left; border: 2px solid #fff;'>Architecture</td>\n",
       "    <td style='text-align: right; background-color: #ccc; border: 2px solid #fff;'>7.7 GB</td>\n",
       "    <td style='text-align: left; border: 2px solid #fff;'>RAM</td>\n",
       "    <td style='text-align: right; background-color: #ccc; border: 2px solid #fff;'>Jupyter</td>\n",
       "    <td style='text-align: left; border: 2px solid #fff;'>Environment</td>\n",
       "  </tr>\n",
       "  <tr>\n",
       "     <td style='text-align: center; border: 2px solid #fff;' colspan='6'>Python 3.7.5 (default, Oct 25 2019, 15:51:11) \n",
       "[GCC 7.3.0]</td>\n",
       "  </tr>\n",
       "  <tr>\n",
       "    <td style='text-align: right; background-color: #ccc; border: 2px solid #fff;'>0.4.11</td>\n",
       "    <td style='text-align: left; border: 2px solid #fff;'>discretize</td>\n",
       "    <td style='text-align: right; background-color: #ccc; border: 2px solid #fff;'>1.17.4</td>\n",
       "    <td style='text-align: left; border: 2px solid #fff;'>numpy</td>\n",
       "    <td style='text-align: right; background-color: #ccc; border: 2px solid #fff;'>1.3.2</td>\n",
       "    <td style='text-align: left; border: 2px solid #fff;'>scipy</td>\n",
       "  </tr>\n",
       "  <tr>\n",
       "    <td style='text-align: right; background-color: #ccc; border: 2px solid #fff;'>0.46.0</td>\n",
       "    <td style='text-align: left; border: 2px solid #fff;'>numba</td>\n",
       "    <td style='text-align: right; background-color: #ccc; border: 2px solid #fff;'>0.9.3.dev17+g698fe17</td>\n",
       "    <td style='text-align: left; border: 2px solid #fff;'>emg3d</td>\n",
       "    <td style='text-align: right; background-color: #ccc; border: 2px solid #fff;'>7.10.2</td>\n",
       "    <td style='text-align: left; border: 2px solid #fff;'>IPython</td>\n",
       "  </tr>\n",
       "  <tr>\n",
       "    <td style='text-align: right; background-color: #ccc; border: 2px solid #fff;'>3.1.1</td>\n",
       "    <td style='text-align: left; border: 2px solid #fff;'>matplotlib</td>\n",
       "    <td style= border: 2px solid #fff;'></td>\n",
       "    <td style= border: 2px solid #fff;'></td>\n",
       "    <td style= border: 2px solid #fff;'></td>\n",
       "    <td style= border: 2px solid #fff;'></td>\n",
       "  </tr>\n",
       "  <tr>\n",
       "     <td style='text-align: center; background-color: #ddd;border: 2px solid #fff;' colspan='6'>Intel(R) Math Kernel Library Version 2019.0.5 Product Build 20190808 for Intel(R) 64 architecture applications</td>\n",
       "  </tr>\n",
       "</table>"
      ],
      "text/plain": [
       "\n",
       "--------------------------------------------------------------------------------\n",
       "  Date: Sat Jan 11 07:31:20 2020 CET\n",
       "\n",
       "             Linux : OS\n",
       "                 4 : CPU(s)\n",
       "            x86_64 : Machine\n",
       "             64bit : Architecture\n",
       "            7.7 GB : RAM\n",
       "           Jupyter : Environment\n",
       "\n",
       "  Python 3.7.5 (default, Oct 25 2019, 15:51:11)  [GCC 7.3.0]\n",
       "\n",
       "            0.4.11 : discretize\n",
       "            1.17.4 : numpy\n",
       "             1.3.2 : scipy\n",
       "            0.46.0 : numba\n",
       "0.9.3.dev17+g698fe17 : emg3d\n",
       "            7.10.2 : IPython\n",
       "             3.1.1 : matplotlib\n",
       "\n",
       "  Intel(R) Math Kernel Library Version 2019.0.5 Product Build 20190808 for\n",
       "  Intel(R) 64 architecture applications\n",
       "--------------------------------------------------------------------------------"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "emg3d.Report(discretize)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}