sql_load_netcdf_example.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Example showing how to get time-series data from a sql station\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "from soda.dataio.netcdfio import load_sql_ncstation\n",
    "\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "%matplotlib notebook"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "load_sql_ncstation?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "dbfile = '/home/suntans/Share/database/UWAOceanDynamics.sqlite'\n",
    "\n",
    "qstr = 'height_start < 30.'\n",
    "\n",
    "varname = 'waterlevel'\n",
    "station = 'KIM200'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "LOWER(Variable_Name) LIKE LOWER(\"waterlevel\") and StationName LIKE \"%KIM200%\" and height_start < 30.\n",
      "Querying database...\n",
      "LOWER(Variable_Name) LIKE LOWER(\"waterlevel\") and StationName LIKE \"%KIM200%\" and height_start < 30.\n",
      "/home/suntans/Share/ARCHub/DATA/FIELD/IMOS/IMOS_moorings_PILKIM.nc KIM200-1208_SEABIRD_SBE37SM_+_P_267\n",
      "/home/suntans/Share/ARCHub/DATA/FIELD/IMOS/IMOS_moorings_PILKIM.nc KIM200-1202_SEABIRD_SBE37SM_+_P_272\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "[<xarray.DataArray 'waterlevel' (TIME: 249355)>\n",
       " [249355 values with dtype=float32]\n",
       " Coordinates:\n",
       "   * TIME           (TIME) datetime64[ns] 2012-08-10T22:00:00.979200 ...\n",
       "     LATITUDE       float64 ...\n",
       "     LONGITUDE      float64 ...\n",
       "     NOMINAL_DEPTH  float32 ...\n",
       " Attributes:\n",
       "     standard_name:        depth\n",
       "     long_name:            actual depth\n",
       "     units:                m\n",
       "     reference_datum:      sea surface\n",
       "     positive:             down\n",
       "     valid_min:            -5.0\n",
       "     valid_max:            12000.0\n",
       "     comment:              depthPP: Depth computed using the Gibbs-SeaWater to...\n",
       "     ancillary_variables:  DEPTH_quality_control\n",
       "     _ChunkSizes:          249355\n",
       "     X:                    121.24265\n",
       "     Y:                    -15.534517\n",
       "     Z:                    25.0\n",
       "     time_start:           2012-08-10 22:00:00.979200\n",
       "     time_end:             2013-01-31 01:54:00.979200\n",
       "     StationName:          KIM200-1208_SEABIRD_SBE37SM_+_P_267\n",
       "     StationID:            KIM200,\n",
       " <xarray.DataArray 'waterlevel' (TIME: 272401)>\n",
       " [272401 values with dtype=float32]\n",
       " Coordinates:\n",
       "   * TIME           (TIME) datetime64[ns] 2012-02-02T00:00:01.036800 ...\n",
       "     LATITUDE       float64 ...\n",
       "     LONGITUDE      float64 ...\n",
       "     NOMINAL_DEPTH  float32 ...\n",
       " Attributes:\n",
       "     standard_name:        depth\n",
       "     long_name:            actual depth\n",
       "     units:                m\n",
       "     reference_datum:      sea surface\n",
       "     positive:             down\n",
       "     valid_min:            -5.0\n",
       "     valid_max:            12000.0\n",
       "     comment:              depthPP: Depth computed using the Gibbs-SeaWater to...\n",
       "     ancillary_variables:  DEPTH_quality_control\n",
       "     _ChunkSizes:          272401\n",
       "     X:                    121.243133\n",
       "     Y:                    -15.534617\n",
       "     Z:                    24.5\n",
       "     time_start:           2012-02-02 00:00:01.036800\n",
       "     time_end:             2012-08-09 04:00:00.979200\n",
       "     StationName:          KIM200-1202_SEABIRD_SBE37SM_+_P_272\n",
       "     StationID:            KIM200]"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data = load_sql_ncstation(dbfile, station, varname, otherquery = qstr)\n",
    "data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "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, fig.canvas.height);\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 overriden (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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4nOzdd5wkRcH/8S+goj8Vc0RcQRFJD2JEFAFFeZSgYs4o5seEaVHCgeScMwdHznm5yOWcc863l+/2bvfu9jZP/f6o6d3a2Z7Y3TM925/361Wv25vp6amp6fCd7upqCQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAoIqdLclI+lSB0x+dnr5D0vsynvtV+rl8ZUV6+svS/3+rpFMKfG2nz2v9LJH0qvP/j2TMp0NSg6Rpkm6QdLjPPPLV6SdZ3ruc9pFti1ckbZSt131ZpvXaLLPsKfC9JmR5fV3GdJ+UNFhSvaRW2XaeJOlHeea/r6SJ6XnelPHca7K8t5H0j4xpj5Z0p6TJkprT03whx/u+S9ItktZKapO0JV1/v2XrW5LGSdqdnvcCSedEPM98y+FtOT5bLhMkzSlw2tNk679d9vPUSxok6WM+03rL2WZJb/J5fr2kF3wef4uk/8iuk42y6+hmSUMk/VDS/s603vr8V2eehWw/fuLz2kznpZ//gPOYu+ynZL+r5ZKeknSW7Hro9zmz1eNVn+mjdLakJyUtk63/iizTnSLpAUlLJe1Vz3d1rM+0X5Q0UNJM2WUis81yeYukCySNlV039kiaJ+mf6v09+/lf9bRj5vr0G0kvyq53LbLf0e2S3ptlXsWsp0C/c7aKC4A3q2flq8147l2SjnPK59PTPZHx+MfT07sh7oCMaY6TtFV2A+E+9lmf1/rJFgBvTM/neNmd2oWSVsnubP6WMQ9vx/svn7odJ+mdWd67nF4ju/GcLOkO2c+RLwCeot6f49MFvtcE2R1IZjscljHdl2UD2I8lnSzpDNmdj5HduWbzV0kblDsAZi5Lx6nvD5Fz0vOpk/SycgfAD0haLWmxpF/K7tTOkg1V78qY9gJJXenn/le2Hf8k6Q8Rz9Nv3ThO0iPpz/blLJ8tn0ID4A3qCfrfkf08v5Zdx1okfSNjeveHxgCf+fkFwMNk22yXpOtkl5kTJH1P0v2yPyTceWWGuGPVu20ekP+y/k6f12bKFgDdZf/L6TYYnJ52jKQ3+3zOzO2XV/x+cEZptKT5kh6SDX/ZAuDzstvN30k6UdJ3JU2V1J7+v+tS2e/sOdnPX0wA/Ljs9v16SWdK+pKkS2S/52E5XvdmSevUs53I3P5vlvSw7A+GE9OfY4Psd5G57hWzngL90tkqPADuL3sEYI7sCrU0z/TeTjtzZ+7JF+KyHSko5LXZAqDfRv//SRqefv4rzuNeAPxmlveIi32dv1uVPwCW+uu2mCNGfmbIhm0/h8ge/fqWcgfAbMuSy22PHyh3AKyT3aHka5PPyAa1zB8J5Zpnpn0lrZG0Uv5HnwpRyPf5U9n2u9XnuTdJmiX7A+RDzuPecjZE9kjZuzNel7lev1Z2fW2Q9NEs9ThYNih48oW4XMt6qQEwW1t5Zz4ezXg81/ar3Nx1YqiyB8DM70qyoWtb+nXZ5unXZrm8SXa7m8mbz3FZXneX7BHiK+X//frV/zj5//gsdD0F+q2zVXgA/H562j9Kulz5T61VSwCUpINkj54Ndx6rlgDoinMAHCp7FCXTPpJGyZ5Oy7bMFBMAXbkC4IfTz51fwHwelg0z+U5PRTFPP18t4n2yKeT7XCL7o+/1WZ4/QT1H1j3ecvYJ2dNqmaeoM9frHyr3uuknTgFQsketujJeE6cA6MoVALMZJ2lhjueLDYDZfCk9n+/6PHey7JHIY1Tctmw/2dPetzuPFbOeAv3W2So8AA6XDRhvk12BUrKnWrKppgAoSdNl+714v2y9AHiW7GfJLPnsk+V1mWW/AuZVqEIC4CbZndVm2X5chW60J8ie8tsh2xdzhexpoGzhYF/Zz/du2dOanfLvL/e79Dzfo/wBcEe6Dm2yRxR/nqfOuQLgL9LP/VI2fDbLtt9o2aNzrrWyp8J+Lhtiu2SXz8tlj2BFOU8/T8q254F5psslX6g5SPazPJJnPttl+y163HXzFtmd9oed5zPX6/vT0x9SUK2tuAXA/0u/5gfOY+tl+6P5rfOFHLUtZNtR6LxcxQbAt8memn8qxzRhBUDvezsi4/H/J1vnKzKmKyQAettxt1tFMesp0G+drcICYI3sDupx57Exsqd/Mvu+eKotAD6Tnubt6f97G45sJVvHYs+vlPv1XunMNoMS5AqAP5f0b9m+ZifLbrR3yF48ktmPzs/lkn4r6STZ/pO3y9Z9pPx3Qvep5zO2ynbQznSQ7M7ll+n/Z1tm9pMNIj+SDXPfkd2RZetn5skVAC9IP9ck24/pVEnflu0rtVfSkc60Hel6Nkj6vWz7XSn7+R+MeJ6Z3iHbnpkX3xQrX6jx+vBemmc+M2Q/h8ddN9+Vfu4J5/nM9XpEevrMH1WZP6DcH0pxC4Cnp1/jns5fr+zrfK6+sFLui54yS7EXoxUbAB+XDfEfzzFNGAHwWNkfd35B8ybZH0nej81CA+BbZLsqrVbvU87FrKdAv3W2CguAF6tvH7mfpR/7VZbXVFsAfFb+AfDvsu2TWfIdBXxnltdllk/mmY9U+NHHXAHQz+dkj+ReX8RrXLWybXSGz3M1sp/v65Lulv0BcW7GNENkT/96ij3VO0T2aOA7sjyfKwBelH5unnr3Z/qA7I5okPNYV3ra72TM41bZ9js4wnlm+ovC6ZoQVgCcKbsj9WSumxfKfh5vG1NoAPyHegcdt65xC4BnyD8AjpH/Op/vx6OyvM6vvD3bDLIoJgB6fe1+l2e6oAHwENn2Wix7xNH1OdkfRSc5jxUSAN8gu23Zo777t2LWU6DfOlv5A6DX4XyD7MbmrenyfvVcgeonygB4cfq12Xb8K2TDgaeQADhD/qeAS93RhnUK2Ku7W7L1vSw2AEp2qISJRb7Gc2C6PpcXMO29skcSvO/sB+n/f0o9y9Q70vO7M/3/fCH7x+r7w8SVKwB6p+38wu802Z2DZ1t62syj3V9XTzeBqOaZaZ7safxCuiHkUo5TwJL0Rtn6jkz/P3O9Hij/U8DvUU/ImaPwAuCH1PPDzo93dOg9zmP52uqP8j8FHKQPYKVPAV8q+5n+VcC0QQLgwbLdIZbL7lNc+8iGwsfUs414q6Rr0+9XI/+hhl4v22Vpr3oHR08x6ynQb52t/AHQ63Ceq2T22ZCiDYC/T7/2f3ye21c2mA5yHsu3w/ig7K9M90q3oAEwrFPA+6vvr36/jZ5UWgBcIWl8ka/xfECFHSWS7LAZRj1HPL3vP1c5Pc88f5KeLttQKLkCoHcBg99OYLp67/BHyj+snabey0gU83R9Ov3cVT7PFStfqJGCXQTirtfe+nqq+q7X3neU68dZZl2DBMA3yh59vSHLa++U/WHi9sPM11bD0/N0A0yQAFjpU8Be+LugwHmWGgC98LdK9gdHpkLaYUbGa7zw1yK7DfdTzHoK9FtnK38AfFJ24/YN2V9TbvF2wNf5vC7KAPhR2dNKfkeevNMx7oYx1w7DGwYmpd5BImgADPMUcKGKDYCfl/3c15b4fv+RbaPTCpj2Mdmw652yOlh9l6cvp+f3TPr/2Y7weobKngLOdhosVwDcT7b/4wL1Pg10kGw73uU89of0fL6XMQ+vH6R3MUYU83TdlX5NtuFSilFIAPSGgbnF5zl3GJga53G/9fo1sn2x5sieSXgh47l8w8CEGQAle8R7pfpegf0G2UGuMwdqztVW3g+9zH6bQY8AVuoU8CWyn2dAEfMsJQB+SDb8rVHv5ce1j/puI06SvYLe+4Hobj9fL3tFdqtsX+dsillPgX7rbPUc5v+OTzlQdoUYnGMeM2UH9cy8cjHKACjZgY9Tsr/YvyF7dOFC2R3SlIz6eDsMdyDor8v+wvUGgv5Lxvy9AJhtIOigV7yF5ST1fF/tsjsv7/9egHqN7A7sH7Ibza/IDoHQKLvDc0937ScbQNwBWU+WPaX+G9l2OUN2I+lN556GGigbKL8nOxjrd2Q7dhv1XMWXTbZl5rz0fH+S/rzfV0/fscyhHN7ofP4b1XMk4zvqu1PwAuKLssvD92WHu9ih3n3wXifbfjtlr2g+RdI1sj+M3KNfUc1Tsj9UmmQHF87Ga79C7jQxQbZzvN96f4IznTsQ9LfVMxD0YuUeCDpzvf6Oeo7aZK7Xh8kGgSbZZed02dB+huyPvEbZoUg8QQPgF2R/OMyUDbkny14kNTf9mTKvBJ2g3gNBnywb/F5Jv89I9T0qv17ZB4L2u7NGlI5Uz3c7W/aUvPd/924uXp/eOvWt82fV27udeTySft1v1Hf58duevFd22dsrOwxQ5nvlu7o92/frDcx9ic88MwffLnQ9/WW6/u6djD6cfuzujHmukf0xA1SFs5X78LrX4TxzI+/6bXqazD5LUQfAfWSPonh999pkjzJcIRsCXJn96DplV/TpKv1WcBfnqFs5TVD2OnpHvvaVvRJzuWxAbpfdAPvdJskvRBwmu3FdL/uDoEV2Z3mebJBx/Ur2lPI22WC9U3Z4hXy3gnPfO3OZ+UZ6nlvT82yS7WCfefRM8u8z6RW/Ix/fkl2GWmWDxvPyv8XZOyTdI3u7qHbZDf3f5N8HK4p5ehdd/dTnOc9b09M8nGMaT67lJjNAnib7A6AhXU/vVnB+602u9XqK/AOgZK/YPF92ndwl+z1vkQ0Ov5M9OucJGgAle/TsednltFN22XpG/le7ZrbVHtllqdRbwa3JUa8o5Opu4Z7mzbVMZHZVybV9dJcfv+1Jvm1rvlPP2Y4y55qn34+iQtZT7wiv3xmlzLMt61X8GIsAAAR2puwRRL8+uQAAAOiHbpS93ysAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAqYh/Z0dMPoFAoFAqFUlXlQPkPKA7kdaByj4JOoVAoFAolviXfLfAAXwdIMvX19aapqYlCoVAoFEoVlPr6ei8AHlDhHIEqdYAk09TUZAAAQHVoamoiACIQAiAAAFWGAIigCIAAAFQZAiCCIgACAFBlCIAIigAIAECVIQAiKAIgAABVhgCIoAiAAABUGQIggiIAAgBQZQiACIoACABAlSEAIigCIAAAVYYAiKAIgAAAVBkCIIIiAAIAUGUIgAiKAAgAQJUhACKofhEAV2zdba54ZZHZtru10lUBACByBEAE1S8C4GEXDDY1tXXmFw9Mq3RVAACIHAEQQfWLAFhTW2dqauvMxy8ZVumqAAAQOQIggupXAfAYAiAAIAEIgAiqXwXAmtq6yN/risGLzL3jVkb+PgAAZEMARFAEwCIs3tRUtqAJAEA2BEAE1a8C4FEDhkb6PjPW7CAAAgAqjgCIoPpVADw64gA4cy0BEABQeQRABNWvAmDURwAJgACAOCAAIqj+FQAvIgACAPo/AiCCin0AbGppN6fdMs7cMXpF1mm8UHZkxAFwFgEQACquua2j0lWoOAIggop9ALxxxNK8oYsACADJ8PLcDaamts7cNSb7QYEkIAAiqNgHwKuGLC44AB5x4RBjjDGpVCqSusxet5MACAAV9NHzB7MdNgRABBf7AHh1EQGwprbO/OvpueaEq0dFcoqAAAgAleXd+z3p22ECIILqdwHQK09NXxd6XdwAGNVRRgBAdodfOIQAaAiACK7fBsBnZ9aHXhc3AHZ1VW8AXNfQbBr2tFW6GgBQtCMvGkoANARABBf7AHjN0NIC4HnPzgu9Lm4A7KzSALh1VysbTwBV6+gBBEBjCIAIrt8GwB/eMzn0uvSHADhu2VY2ngCq1jGXDGMbZgiACC72AfDaoUtKCoA/uW9K6HVxA2BHZ1fo849C4972Xv8nAAKoZkdxCtgYQwBEcP02AP5s4NTQ6+IGwPYqCIDnPTvP1NTWmS9dN7r7MQIggGrmbueTjADYv/1e0jxJu9JlsqSvOc/vL+lWSdslNUt6SdIHinyPWAfAKSu3F7Sye88f8u9Xuv8++/5oA2BbR/wDoF/bjV+2jY0ngKpFALQIgP3bGZK+Lumj6XK5pHZJR6afv1PSekmnSDpW0ihJcyTtV8R7xDoAZh7VyzedGwB/+cC00OszxwmArR2doc8/bARAAP0NAdAiACbPDknnSHqLbBj8vvPc+yV1STq1iPn12wD4qwenh16f2QRAAKgoAqBFAEyO/ST9QFKbpCMkfUn2i39bxnRzJV1SxHz7bQD8zUPRBsCW9uoMgBOWEwABVC8CoEUA7P+OlrRHUqekRtlTwpL0I9kwmGm4pLtzzG9/2YXFKweqnwbA3z08I/T6zCEAAkBFEQAtAmD/9zpJH5H0KUlXStomewQwWwAcIemuHPO7WHaB6VX6SwA8+Lyeaf/w6MzQ6+MGwL1tBEAAKJfNTS1md2sHATCNAJg8r8oe4Sv1FHC/PgLolj8+Niv0+rgBsLmtI/T5h+0I556Zr8zbaH770AwzZP4mNp4AqsqWXS19fuQnfRtGAEyekZIGqecikO85z71PCb8IxC1/fpwA6AZAr5x563g2ngCqivvD1W+fkEpV552ZgiAA9m9XSDpB0odk+wJeLhvwvpJ+/k5J9ZK+LDsMzEglfBgYt5z7xOzQ6zPZGZdw5dbdoc8/bH4BkF/PAKrNkPkbfbdhy7fsMpNWbDefvHSEGTJ/U6WrWVYEwP5toKQ1sn39tsqe/v2K8/zrZQeCbpC0V9LLkg4q8j1iGwDvGbsyb2B5dma9mVu/03fD8Lcn54Rep4krevrPLdm0K/T5hy1fANzdGv+jmADw0pwNObdlSfxRSwBEULENgPlW7udm1efcEPzjqfAD4KQVPUcAF22MX5tlOtK5Z6ZfaWppzz8TAKiw+yesIgBmIAAiqKoNgJ+8dETODUHtM3NDr5MbABduiF+bZSIAAugP7h3X94wQAZAAiGCqNgD++fFZOTcE5z07L7S6bNnVYoYu2NTrIpBtu1tDm39U8m0sCYAA4m5XS7v5xH+HEwAzEAARVFUFwO/eOcl0dtmrvfz6CLrlP8+FFwD/5+Jh6VA5t3v+DXvaQpt/VAiAAKrdX5+YXVD4IwACxamqAFhTW2emrmowxuQPgBc8Pz/0unzm8p7Tzv0hADbuLV8AfHnuBvP0jPqyvR+A6Lw8d0PZusEUGv4IgEBxqi4AeiEiXwAc8OKC0Ovy2ctf7f57ez84BVyuANjR2VVVp87X79xrPn/VSHP32BWVrgoSam9bpxnw4gIzeeX2SleljynOcFjlQAD0RwBEUJEGwCCDc2Zbwc971l7ckS8AXvLSwrA+hm8ArIYgc1Sei0Aam8sTAFs7Orvfc832PWV5zyDcU05AJVw/bElsl8FBE1cTAGOAAIigIguA/315ofn0ZSNKOlXa7hwxyize1b35AuA5g6aF9ln8AuCSTbvMHx6ZaSYs3xba+4Tt2Dwdp3c2t5nFm5pMW0dXpPWotgDoXmAEVMIfH4vvMviAMyRLORAA/REAEVRkAdBbIa8ftqT4BbulPXAADHNj4BcAz7xtQuw3OvmunDvv2XmmprbOHHr+4Ejr0dJeXQHwXI4AosLiHADvJwDGAgEQQcUyADY2Zw+A/3q6sAB45EVDQ/8sbgCsho3OJy8tbOiEqD/DlqaW7vfxu4XevPpG89+XF5b1opRc/vbknNh/t+jf/hTjADhwfPkCYOPe7PsCvxKXbUg5EAARVOQB8LoSAuCOPW2BA+BldeH3Aay2X535Bssu12dYsKGx+32emr6uz/Pec/98Ovy7t5TCbZdX5m00h/z7FbOiCu79XElB+vuirzgHwPvKGAA/d4X/j+5s5a4xyblwiwCIoGIZABsKCIC3vLos54bgqiGLQ/8sbnFPr8bVpy6LRwCcv74nAPptoL3nzrh1fKT1KJRf+xz6n2hPk1ezS19eaI6/cmSvi4pmrGkwJ1872oxZurWCNatecQ6A7l05olZM+CMAAsWJPABeO7T4ALh1V2veAHjbqOU5NwTXDI02AB7n/DKNq49fMix2AfDOHAHw9FviGwCPHhBel4L+xu9H10fPHxz79SPOfnzvlNi2HwEwHgiACCryAFhKENvs9BnLLN5pwrvHrsi5ISjlyGO+z+IWd1DouJm/vtF8+46JRW04o7Rk066CjgCedsu4SOtRKL/2OfGaURWpy40jlpqz7phoWto7K/L+hfDa6Cin3225lq3+Ks7t53a/iVqxATBJY3cSABFULAPg+p17s67gf3l8ljEmfx/AG4YvDf2zVEMA3Nmc/fR5pQKge0T3wUmr+zzvPXf8lSMjrUeh4hQAvfd/dMrairx/Ibw6Hn7hkD6PxW39qBZxbr84B8DbRy+PvE5xQQBEUJEHwKtL6Iu3rqE56wr+qwenG2N6n4bwKzeNWBb6Z3HLpy8rLAA2Nrebvz81x0xcUZ7xAhdvaopVANzV0m7+8OjM7vcZNHF1n2nceqxraDaLN5X3zjS3j15uLn6p584xfu1z0rWjy1qnzLo85BOcK827G45XxyOr6AhgKpUy/3p6bqjbiTDFuf3csy9RK2Vb9uzMZNxykgCIoCIPgKVcjLF6256sK/ev0wHQvRLNr9w6MtoA6F5hm8t5z84t64bcveI2DgHwgufn93qfEQs3m2ELNpkTrxll5qzbaYzxb9+o77QydVWD+cvjs3odnVyyaVfW+lQ6AD4yZU1F3j+bu8bYEOD+EDviwiHmnEHTzLVDl5Rl2QpiXn1jrOsX5/aLewD06rVi626zdVf879hUKgIggopdAEylUubWkdmv8L3lVRvsBuYJgLeNCu9UgN/8C70K+Pt3Tyrrhty94CIOAfB7d03q9T6ZRyjd0/1HOreum7l2R2R1MqbnO/3dwzO6/56dI5CeXOEA+NjUeJ0CjsOyFcTkMt/PNpeurr5D6MS5/bzwH+cAuLFxb2zbLywEQAQVeQC8cnBxAfC5WfV5V+4tu1p6jUbvV+4YHV5n4CAB8Dt3Tizrhmhu/c5YBUD389fU1plFG7OfonbvrvLqos3mgufnm/nrGyOpl/c+p90yrvvvXEckT75udCT1KLSeT0wjAIZp4optsajf8i27zTGXDOuzvYpj+6VSKTNu2VZzxSuLuuu2uanFpFIp0+kTYku1cEOTOfPW8Wbk4s0lB8ARCzfHrv3CRgBEUJEHwCsGLyrqdR/+9yt5V+5zBk3rdT9KvxLm1WB+8z+2wAD4rdvLe8u4WWt3xCoAftP5/DW1dWbqqoasdfiGEwCjvsjGm/fpt4zv/vuhyWt6BYM4BcAnfQbQrqQ4LFtBTFgejwD4s4FTfesRx/bz+3H+43unmJ8NnGqOuWSYaW7rCOV9ihm/lABIAETpYhUA3XvG5ipfuWGMGTRxdc5p7h23MvTP4pZjnDH2cin3PYPjFgDPuHV8r/fJDITZAmDUdfPmnVm/bOVLZQ6A3mlB7/2fmRGvju1xWLaCGLt0ayzqV00B0O+7dW+P+cCEVZG9DwGwLwIggoo+AL5SeAB078Gaq5x641jz4KTcAXDg+HA2Ru5nccv/XNwTAO8bv8os3bwr72vLYfrq7EfYyr2TzjWgt1/JFg6jEOcAOGzBJnPkRUN77cSem0UALMSqbXtMW0dX3ulGL9nSXb9K3sYurgFw/c69pqOzdzsW8l2H0fUmjAD46iICIJBP5AHw8iICYKEr99dvHmceyhMAw/o1Wky98r22HMYv8z+FWYmddLF1+PrN47LWLeydtDfvk64dXVDdvnz9GNPe2WVaO6IfkNnv/V+YvT7y9y1GpZctz6NT1prv3jnJNDa3dwfms+6Y2GuaQRNXm5fmbOj1mNu/zLv6uxJ+fn9PAJywfJv525NzTOPe9ooGwFHpcPzTgVN7PV6u75sAWBgCIIKqygB4xq3jzUOT1+Scxm/A4aCfpRoC4BTn6sZK76SLrUO2/p/XDVtiPnXZCLO5qSX0uh18XmF1+9J1o80n/jvcHHnR0D5HRsLm9/4vZgSYSqv0spVZj8tfWdQrTC3e1GSuGrK4z3AvnV0p8+N7p/Qax3NefTQXGuWzZvse3/Y6//l5Zd9uuPyOSnZ1pUL9vsct22qmr27wfa6U7VdmOWfQ9Iq1X7kQABFU5AHwsrqFRb8mX8l2pMgtYQ6cG2TDV84N+fItu0veYEYhjA25Wy56YX7F6vY5597P6xqaQ6tHoXWrm7sx0vcsVqWXrcx6nPfsXPOT+6b0ee9Trh/Tqx5+XSRmRTzkkDHGPDRpdZ8fw9nu1e0GsJpaG2ZvGL7U7GkN5yKLfNz39468t3d2hfZ9b9vdmnPasLcb/RUBEEFFHgAvfTn8AFhI+cOjM0M7UhNkQ1OuDdGMNaVd/BFV3Xa1tAeqj18Z8OKC/G9coGLf2w2AQ+aHH8ZuHLHUnHnreNPc1uH7/nG7FVxcdsDee/zz6fz9h1vaO82NI5b2eXzM0q3d81u7vdncPXZF6GHLey/3jkDFLoP/LWJbGoQbAE+5foxp6+gq+AK9Qr7vpZt77g3uXezU1tFlLqtb2OvinLBKf0UARFD9NgDW1NaZz18Vzr1lS9nQeL+cy70jjMNGsq2j8KMFldoBBqlHFP3xvHlnG9/yf28aZ/78+KzY3BEkLjtg7z3+/lT+APhl52igW9zuIodfOMTU1NaZfz83L7Q6uqH+nEHTTE1tXa8B4gst37x9Qmh1ysU9lV5TW2dGLdnSp19ikPXUPVPh/UjP16c7SPEbaLs/IAAiqMgDYDE77ShW/jA/S76yrqHZXDdsidnS1FKWHWFnV8oMX7i51ymVSraTZ1Oj/+cPWv7x1JzQ6hikHo+HeFeOzOFe7nTuslCO76rYuv7wnsnm3CdnV2zZyuS9R7F1cos7YoD32InXjAqtjrlubVnJ9TQb9zaXUdRzlbizueAAACAASURBVNMe3kVVVw9ZHMk2o6a2LvI+u5VCAERQkQfAS15aaLbvbi3oKs4oVv4wP0sc6uLyjhaFMXBqWDJv9RbH79OYYN/pPWP9x5icsHybuXbokoJ3ODPW7DBHDxhqnpy2rnvet49eXvblqFCz18XrLjPG9HyPv39kRkl1q6m1Y4Zu3WW3Ud5jx1850uxsbjO7WtqLqk8qlTIz1+4wTS3t5uW5G8xRA4aaYQs2xW75zyXKemb2JdzbZgPgbx8q/fvLV8px5X4lEAARVOQB0OvoXMjVwFGs/GF+lqAl7KsNf3jP5Fi1kzHRfIdxq+cp14/JOc9CjxB+2ie4xzkAvjJvY8nfWUt7Z0Hj8xXLe4+z7phYUt1qav2HAjrKuS91MYbMt230+atGxnr5zyWserZ3dpmrhyzu1e/RHX+xprauu69lmNuyzOKFzP6GAIigIg+AxWy84rrBDLM+v3loumlpD2eDRACsTD2PGjC0ez6pVMqsa2judfTomqGF3f+61CO3xtgd6eWvLCrr6a0vXjOqpPou2bSrV93D5M332wECYL5SzBiUXh+/uC//uYRRzyenr+s1WL/nhdnre03XlD7C+q0cdwgKWna3dpjnZ603CzeEv5+rJAIggur3AXDs0q0mlUqZ3zw03fz1idmhfZYg5a4x4dyn2L0yNWhJpVJmx562QPV5OM/YjGGUMLhXIZZSfvXgdFM3d6NZvKnJ3D3W9tm79OWF3c9fO3RJzvdfs32Pae3oNJ+8dHhJ73/H6J5+gmH2R8ymqytlWjs6zdEDhpbcZu5yFiZvvoUO6F1KKeYigmpY/sv1Gfxug3nD8N5XYe9sttucbLeBDKO8NGdDWduvXAiACKrfB8CaWntxhvd3Kf1Bwq7PlYMLO0JUznp95nJ7NOqJaaUHiqi+v2KWoUJ4V3pGVa4blj0ATlphB+r+2k35x7IspNQ+MzfyW5l9965J5pAsg3QXWzpDviKzHMtcMXWuhuU/l3y32CymuD9wPB/KGHh9++7WsrRbudqvnAiACCoRAXDt9p4AWMxpgLaOrkguarh6SPwCYJCNZCqV6jN0RFzqVs5280quABhFX6d/Ph3e1dFRt1d7SKesd7d2FHV3iiClsbnwC0GqYfnPZuKK0m4jma189vKeMxSpVMrU72juM83WXQTAUhEAEVQiAqB7BLDQAXXd0Bh2+duT4eyw47KRjPrKX7f8+fFZ5v6A93mOuo7XZwTAtdubu48ixeU7q1R7BbkQ5OW5G8zA8au6b6H243v73vkjijJ/ff6Lt1o7Oou6W0Ycv+ewf8S5F8L4DcBdU1tnNjbuNU0RDBpfifYrNwIggqp4AOzsSplxy7aaxuboNgKlBMAgVxWWY0PU2lH4yPxR123qqr6314rzhjzqul0/fGn3ez09o97U1NaZ/3t0ZqTvHaUw67loo93W7NjTZq4fvtSs2b6n6HqE2fe1kDK3fmfOenkDnx95UfA+kv3le66p7elWEqfSnxAAEVTFA+AD6bHsvnRddJ24V2ztGXm+0AB4wtWlXfFYrg1RlHXb2dxmGoq4ICSK2zdF2X5R1+364UvN7HU7zfz1jeZkZ7ketyy6dopSFHX99YPTTU1tnTnauaK63PUotOQ7ArhgQ2NVLPdxbd9ylv6EAIigKh4Aoz7SVlNrT7l6fxcaAI+/MppxvMLaEJVjY1nocDXuVXbVsCGPum6XvNRzRfDJEV6d6pbnZ603J107OrQhhqJsr8wfDHH53rKVBRtyB8Bpq8t3BPwrN4wxN41YVvR3uLFxb9aLhTq7UqHcTagaSn9CAERQZQ+A3kboiWlrzdAFm8wZt44v6wbg0Slr8141GWX/vzA2RMXcmD1Iybfjy/VdR11WbSv81KEr2236wix/faLntmQnljh2XqnlzFvHm47OLjN/fWOge6Bu391q/vn0HDNjTfThJtOijU3mFw9M67X8ueMslrusa2jO2VaV6AKxpaml4O/SO8ty+SuLzIINjWbRxibT2Nxuzn1ytlm7vTmUOwlVS+lPCIAIqqwB8PrhS81xV7xa1l/MmeWrN4w1H79kWHdfpELrHnb58vVjzMnXji66XVc6p7OjLgNeXGCenL7OLN+yu+jvOuoyZeX2gtvsyenrzOSV2805g6aXpW5/fnxWr++53G3jDUZ8y6vFHynylKutamrtTrm5rcP88bFZpm7uxu67cBxx4RAzt36n+fYdEytylNkrbj/FaasbzA3Dl/YagDvKU/vZSr510hXWED79ofQnBEAEVdYA6JWwxkALUr5x2wTT1NJu1u/c26veHWW4ks8t+Y4uuLzbTMVpw1nq7cGClqmrGgpqs1LvXxuk/OmxWWV/z2K/t2wenLTavFjmsPWFq0dmvUo0DmX5lt1m+MLNvY4ePzRptenqSpnGve3m0SlrK1KnQmWOvZfk0p8QABFURQLgV28YW/ENgVs2O6dTBo5fVdb3HrV4SyhtWqkNZ6Xqc8PwpeYXD0wzk/McCcy89VQ5ypll7taQrxR6t5DV2/ZUrI7lvrK3mHKuc0rfKwNeXGAO/c9gU1NbZz52QbQDi/uVAS8uMDW1daZ+h/8PyEUbm8z3755UllP41VTCHoi8kgiA/du/JU2XtFvSVkkvSDosY5oxsguAW54o4j0qEgCjvOK3lDJ84ebuekd9l4jMctOIZaZhT1vOfokNe9rMec/OrWgblfI9V7puxhhz/4TyBvqa2rpYHnEpxPQKds2glF42N7WY2evsUDUt7Z1l6yNcreUfT83JOrh3R2eXeW5WvZmycruZV19YH+hKIQD2b0MlnS3pSEnHSKqTtFbSG51pxki6R9J7nfKWIt6jIgHw0PMHV3wj4JaX5mwoqN5RlMMusG3xo3snZ23LPzw6s+JtZIwdu829sKCSHfMz61bKcpiksrFxrxm1eEvOHxqDJoZ3GzBK+cvE5eHeyaO/F1dHZ5fZ09rRZx1Yu73wLjrlRgBMlnfJftlfdB4bI+mmAPOsSACMW3E77VeybGlq8d1Jn3pj5U+Zu/3ajDHdd2OIQ8nm6iGLK163uBW/08GpVMq0dXSZXzwwreL1o1DKVVze1fqZg3kPW7AphL1hNAiAyfIR2S/7KOexMZK2SdouaaGk6yS9Occ89pddWLxyoAiAsSt/eXxWr7aM2xHTe8etrHgd3OJnVQX7s8W9ZIrjHRsolHKVzhz3lB4yf2Po+8awEACTYx9JL0kan/H4ryWdIhsKfyBptaQROeZzsfr2GSQAxrBMX91gxi/bFpvTrHEua7bvMT8bONVs2WUv5innUDnVWn738IzuoUwqXRcKpZIlV5/0KwcvDn3fGBYCYHLcLmmNpA/kme6TsgvEJ7I8zxHAKiuPTS3/EBPVXNz7PlPyl0krtle8DhRKXMtHzx8c+r4xLATAZLhVUr2kgwuYdh9J7ZK+X+C86QMY8xK3K6YpFAolSSWuCID92z6SbpO0QdKhBb7mKPW9UCQXAiCFQqFQKFlKXBEA+7c7JDVKOlG9h3l5Q/r5D0u6SNKnJH1I0tclLZY0S9J+Bb4HAZBCoVAoFJ/ytZvGhb5vDAsBsH/rc7FGupydfv4gSWMlNUhqk7RC0s2S3l7EexAAKRQKhULxKbePXh76vjEsBEAERQCkUCgUCsWnTFy+LfR9Y1gIgAiKAEihUCiUspYnplXH6AYTVxAA0X8RACkUCoVStjJ0wSbTlWPw5TgVAiD6MwIghUKhUMpSBry4oKr2EQRA9GcEQAqFQqEUXZ6eUV/0a64YvKiq9hEEQPRnkQVA7i9KoVAo/bcYY0wqlTKjFm8p+DUEwPAQABFUZAHQu4PF128eV/GVmEKhUCj5y6HnDy54Ws+23a0Fv+alORu6X1cNtyEkAKI/izwAjlu2tXtl+v7dkyq+QlMoFArFv3zqssLP3HiKuf92U0t7r/1EZ8wvBiEAoj+LPABOWbndpFIpY4z9pXjN0MVm7faeDcbpt4yv+EpOoVAolOKKZ8aaHXmnbWppN1t3tfruKyr9OXIVAiD6s7IEQD8t7Z1m2uoGk0qlzJPT15lzBk0zSzfv6l7xfnjP5Iqv/HEp5z07r+J1oBRXvn3HxIrXgUIJu3zlhjHdf3sKOQKYSyU/z+SV282e1o6szxMA0Z9VLABmM3VVg1mxdbfZuqvV3DhiqdnYuLfiG71KlIHjV3X//dKcDRWvT6XLnHU7e/1/+uoG87uHZ5ia2jpz6o1ju3dMcfnhcPFLC8yoJVuK6iBfrvKbh6abmto6c+I1o8yWXS3mjFs5Cl9s+fWD0yteh0qUsUu3mv88N89MdrbrW5pa8r4ul80FvD6KMmlFz2dwg61bCIDoz2IXAP1UeqNXrvKXx2eZSSu2m8emrjXGGHPes/PM+c/PM11dKfOT+6aYy+oWVryOlShfvGaUMcaYxr3t5rGpa83O5rbuZWNj416TSqVMS3unGbl4s9nb1lnx+tbU1pkbRyztrqMXuOJSdrd2JHo9K7Wcfst48+fHZ5kbRyw1TS3t5m9Pzql4ndxy86vLAs/jyIuG5nx+xpodfZab7c5FIL9/ZIbv6/Ipph9h0LKntcMs27yr1/uv3rbHd1oCIPqzRAfAPzwys+Ib7ZraOrN9d6t5Zd5G097ZVbG2iGNZvKnJDBy/ymxuailqefnIf14xNbV15jt3Vu407PbdPf2d9rR2mOuHLTELNzRVtD3bO7tMc5t/+CvXsnVUnoDxPxcPq/hyl63cOWZFr/ZauXV3xetkjL2QolzfoV8ANMaYm19dZm4fvdy0dXSZByas6tM2hYiivo9O6XvLuWxa2jtNKtX7ohQCIPqzRAfAbL/6KrERr3RbFFM+f9XIkl7X0t5pRi/pOSWa74hFscHP07Cnzcxet9O0d3aVvW06u1JmT5YjbMYYs2zzLnPTiGUFh4cRCzebAS8uKLk+XleCxZvyr+NBP/vO5ra809w5ZoVJpVJ9drReeWLaWnPpywvNHx6daXY2t5mz759a9u8ws1w1ZLH56g1juy9m8+TqOxZl+eSlPVfqZjrmktwB+rOXvxrovTc27i14PRy+cLM57opXzYqtuwuaPoq2ylwmfzZwalH1IACiP0tcADzh6lFm/LJt5uW5djyqdQ3NZvW2PebQ/xQ+/lUY5a4xK7JuxMvVFoUU75TOoecPNvdPWGU+fskws3hTk3l+1nrzqwenm1tH9oS48cu25ZyXF4zGLdtqRi/ZYjq7UubBSauzTt/Y3J6nNfJbtW2PuW3U8pJDayHlQ+fZfzOPEOWTLQS5xRh7ZOK8Z+cWXB/v1ORP7ptSVH1OvGZUoHYwxvTqu+pXxi7d2v1+7unG1dv2mGdm1Pc6muUp9zLv97lK/f6iKMs27zLjlm01Le2dfeo0ZeV2c8SFQ8wT03qOfH3vrp7ht8YutcNy/emxWUW95/z1jWbi8mjDUNjtdP+EVb3m+9OBU01bR3FnWQiA6M8SFQB/dO/krL9GF2xoLOtG/NVFm/scUSjE1UMWm5raOnNfnh1tWKWlvdM8P2u92ZY+pZlZZ++m7sdfObL7/5nf2a8fnG7++NisrJ9pyPxNZuqqhl6vuXtscWGqEIV0Vi+lNLW0m0krtnd/9mL4ze/UG8eanw6cam4YvjTvtH7lhdnrza6W9qLrs6e1w7y6aLM54epRBY8H514MYYwxgyb6B/r56xvNMzPqe72fd4HMoImri26jbOUn900xf38qvL55p944NmfdOjq7TGtHZ1n7A65raM5ZJ+97b9jTZh6butbsamk3SzbZ0GiMMa0dfU91euUZn9u7/bSAo2ZhCLONnpq+rntbdduo5eZrN40zjXsL+0E5Y03PtogAiP6sKgKg16cr3+mNXOX6jJ2pn/odzWZXS7s55Xr/K8LCLEHaxetfVo6dTancHUqh1m5vNj+/f2qvKwzD5J4OeijHkcdiynfvnBSoTpnzu+iF+VmnbdjTZlZu3W3++3Lfi4G+cduE7r+XbNqVdR7F2NlsA8R941eZVRndJW5+dZn5/SMzTGdXygyZv9Gs32lPDU5e6X93h2wKOSLz1RvG5vwObhi+1NTU1pmZa3v6pwX9XrftbjWPTV2b83S+a9HG6Pt3fvGaUd0/tMIwfXWDOfPW8Wbp5l2msytl9rbZI4oTlm8zm5taTGtHz1Bd5RCH7ZbHuyqYAIj+rCoC4K6WdrOuodn8+fHiTlt4JXP0+XzWNTQXNSJ+seXyVxaFslFdvKnJDF+42Zz75OxQ6/d4+irkIFo7Os2qbXsCzyds94xdaQaOt6eG8p2yzlfGLN1qWjv6noYrhju/Q88fXNBrUqmUueD5+WbSiu1m4YYmc96z88zmphZz/4RV5qohiwPVJxfvNmHfvmNizulemL3ezF/f+4h6EK0dnaZ+R7PZssv/CK4xps8FVPm+u8Mu6NvlY8j8TWbQxNVm1lr/Cx3yuWfsSvPinA1m7NKt5vRbxod+0U9/F7R9Trh6lLlt1PLu7j1BEACRBFURAD1/LLLfSpANZ1dXytQ+09Pv6sUQx+ILm3sxRWdXquTTUbtbO8ycdTtDr1+cuX0Yiyn1O3KfhiuUO8+GPW35X1BBLe2dZtyywkPvJS/ZI5XnPz8vtDos2thk1jU0mx172sypN47t1afQ9bkrei52mLl2h7lm6OJeFyC5dyC6+KUF5lu3TyjoKvxi+S07x/53eM5l6/ALh3T//eiUteaj5w/Oe9q3Pwi6Xf3aTeNCqwsBEElQVQGwlGFbiu0In82ulvbuX5kHn1f6Rsqv43ZQ23a3miMuHGLOfXK2MSb7abhyB9NqUcm2Gr5wszn9lvFm+ZbCrpSsJh2dXWbOup2mI4JglY97JxbPtNU9fbvWNTSb/3t0pplbH+0PHr9lJ99QN94p7aStk9cNWxIoAP7ygWmh1YUAiCSoqgDo3fmh0HLl4MXd/VrCsLet03R0dpnbRy83NbXF36Jt5OLNodUlU+ZO9uW5hR2xPHqAvRLzM5ePiKxucZfvylXCcvVZ19Bszr5/aq8duHsRULn4LTufuTx395KW9k5z77iV/fJHQS63jVpeUvDzuuuEua0nACIJqioAFntHhUI7cBcrlUqZ5Vt2dV8Bm60cfuEQ8/DkNea+8at8h7eIWiqVyhtSF29qMn94ZGbidjZ+6nf03I3gYxcM6dNW7hhqqE5z63ealQWOSxcGb3k5//l55ukZ9eaRKWvM/940znddvH7YEnPP2JVlq1vceD+s4/BjjACIJKiqADhxue20/5P7pphJK3qf5vQLh+WQ6/2vHBxdh/xCdXalzPz1jaYzS1hFDzcAnjNoWvffwxduNt+6fYJp2NNmtuxqCWV8QiTDjDU7zFVDFvfq+vE1nwB4x+jwhz2qNt7t4L51+wSzYEOj+e1DM8zJ6f0IAbAvAiCCqqoAaIwxmxpbeh1Nczes7p09ro7wakjXqm17zITl28w9Y1f22ij9+7l5ga8QDZtXN/fuCujR2tFzH+Ftu1vNP5+eY6auaqh0tdDPnHZLTwBcsKHR3DF6RSQXoFSjppb2PmdL/ELfP9JjPUa1nScAIgmqLgDGVXNbh/nRvZPNgBcXmGELNlW6Or68jeeLczaYc5+cnXcA3iTa2dxW8ICxQCkembLG1NTWmbPyDKcDK9tRv/U790Y2RiEBEElAAEyQC56fb75+87jYHZkEkiSVSpnZ63aGetFCf1aJrisEQCQBARAAEFsEQH8EQARFAAQAxFZm+Dv9lvGRvycBEElAAAQAxJYX/P78+Cyzt62zLENqEQCRBJEFwJMJgACAgLwA+OS0dWV7TwIgkiDyAMgwGgCAUhEA/REAERQBEAAQW14AfGLa2rK9JwEQSUAABADElhcAn5+1vmzvSQBEEhAAAQCxdevIZeanA6eato7y3S2FAIgkIAACAOAgACIJCIAAADgIgEgCAiAAAA4CIJKAAAgAgIMAiCQgAAIA4CAAIgkIgAAAOAiASAICIAAADgIgKu3fkqZL2i1pq6QXJB2WMc3+km6VtF1Ss6SXJH2giPcgAAIA4CAAotKGSjpb0pGSjpFUJ2mtpDc609wpab2kUyQdK2mUpDmS9ivwPQiAAAA4CICIm3fJftlfTP//LZLaJX3fmeb9kroknVrgPAmAAAA4CICIm4/IftlHpf//pfT/35Yx3VxJlxQ4TwIgAAAOAiDiZB/Z/n3jncd+JKnNZ9rhku7OMp/9ZRcWrxwoAiAAAN0IgIiT2yWtUe8LPLIFwBGS7soyn4tlF5hehQAIAIBFAERc3CqpXtLBGY+XcgqYI4AAAORAAESl7SPpNkkbJB3q87x3Ecj3nMfeJy4CAQCgZARAVNodkholnSjpvU55gzPNnbJHB78sOwzMSDEMDAAAJSMAotL69NVLl7OdaV4ve4q4QdJeSS9LOqiI9yAAAgDgIAAiCQiAAAA4CIBIAgIgAAAOAiCSgAAIAICDAIhi/LmIEifRBcBrCYAAgOpDAEQxVhdYVlWqgllEHgCnrSYAAgCqBwEQSUAABADAQQBEUK+TdJik11S6IjkQAAEAcBAAUar/J2mgpM50OST9+C2SzqtUpbIgAAIA4CAAolQ3S5oh6QuS9qgnAJ4paXalKpUFARAAAAcBEKVaK+m49N+71RMAPyJpV0VqlB0BEAAABwEQpdqrntDnBsBjJDVVpEbZEQABAHAQAFGqsZL+lP57t6SD03/fJmloRWqUHQEQAAAHARClOl72VO+dklok3SRphGx/wE9WsF5+CIAAADgIgAjiaEkPSlogaZGkR9KPxQ0BEAAABwEQSUAABADAQQBEqUZLOkfSWypdkQIQAAEAcBAAUapbJG2S7f/3rKRvyt4VJI4IgAAAOAiACGJfSV+VNEh26Jcdku6RdGIF6+SHAAgAgIMAiLC8XtJ3Jc2R1FXhumQiAAIA4CAAIgzvlfRX2VvDpSRNrWx1+iAAAgDgIACiVAdI+oXs2H8dkpZKGiB7K7i4IQACAOAgAKJULZI2yg4A/ekK1yUfAiAAAA4CIEr1VdmLQKoBARAAAAcBEEG8RtIpkn4r6c3px94v6U0Vq5E/AiAAAA4CIEpVI2mxpGZJnZIOST9+k6S7KlWpLAiAAAA4CIAo1QuSHpYd/Hm3egLgiZKWV6pSWRAAAQBwEABRqu2SDkv/7QbAD0naW4kK5UAABADAQQBEqXZIOiL9txsAvyBpS0VqlB0BEAAABwEQpXpS9rZvkg2AB8te/DFS0gOVqlQWkQXAkwiAAIAqRABEqd4vO/jzItmBoCfLnhZeIundFayXn8gD4HQCIACgihAAEcQbJP1S0m2S7pD0q/RjcUMABADAQQBEEhAAAQBwdAfA5QRA5HdmESVOCIAAADgIgChGqsDSVakKZkEABADAQQBEEhAAAQBwEABRitdKGi3po5WuSIEIgAAAOAiAKNU2SYdWuhIFIgACAOAgAKJU10u6qtKVKBABEAAABwEQpbpVUpOkmZLulnRDRokTAiAAAA4CIEo1OkcZVcF6+SEAAgDgIAAiCQiAAAA4CIAI6iOSTlXPLeD2qWBdsiEAAgDgIACiVO+QNFI9Az8fkn58oOwFIoX6oqSXJW2U/ZK/mfH8oPTjbplSZF0JgAAAOAiAKNVDkoZK+oCk3eoJgF+VtLCI+XxN0mWSzlL2ADhE0nud8vYi60oABADAQQBEqTZLOib9txsAD5a0p8R5ZguAL5Q4Pw8BEAAABwEQpdqtnoGg3QD4aUkNJc4zWwBslLRV0jJJ90p6d5757C+7sHjlQBEAAQDoRgBEqV6RdGn6792yR/72lfSUpGdKnKdfAPy+pNMkHSXpDElzJC2QDXnZXKy+/QYJgAAApBEAUaojZI/KDZHUJulpSYtkTw1/uMR5+gXATO+T1C7bZzAbjgACAJADARBBvFfSJZLqJA2WvZjjfQHmV0gAlKTlkmqLmC99AAEAcBAAUaoPKvuYfx8scZ6FBMB3SGqV9LMi5ksABADAQQBEqbrkfzHGO9LPFepNkj6eLkbSuem/P5h+7jpJn5P0IUknSZokab2kNxfxHgRAAAAcBECUKiXpXT6P10hqLmI+J8nngg3Zq3/fIGmYbF/Ddklr048fVGRdCYAAADgIgCjWDenSJeku5/83SLpZ9i4dEytWO3+RBcATrxlFAAQAVB0CIIo1Ol1SskFvtFOGSbpbPeMDxkXkAXDGGgIgAKB6EABRqgdUPV8IARAAAAcBEElAAAQAwEEARBCflnSNpCckPZdR4oQACACAgwCIUv1A9srcOtk7gbwsaYnsfXsfqGC9/BAAAQBwEABRqnmS/i/9925Jh8gODH2P7N1B4oQACACAgwCIUjXLDs4sSdslHZ3++3BJmypRoRwIgAAAOAiAKFW9ekLfXEk/TP/9OUlNFalRdgRAAAAcBECU6jFJf0v/fb7s3TrulbRGXAQCAECsEQBRqrdLen/6730l/UvSS7J3BHlbpSqVBQEQAAAHARClelTSryV9tNIVKQABEAAABwEQpbpbdtiXlKSNkh6X9DtJH6tkpbIgAAIA4CAAIqj3yo4JeJekxZK6xFXAAADEGgEQQb1R0qmSrpQ0WXZQ6NkVrVFfBEAAABwEQJTqaklTJLVImi7peklnSnprJSuVBQEQAAAHARClSknaIuk82cGf44wACACAgwCIUh0j6c+yY/5tk7RZ0pOSfq/4BUICIAAADgIgwnKMpAckdcheCBInBEAAABwEQARxrKRzJb0oaYekTkkzJF1byUr5IAACAOAgAKJUO2WP9s2QdJ2k0xXfL4gACACAgwCIUsU58GUiAAIA4CAAIgkIgAAAOL56w1gCIPo9AiAAAA4CIJKAAAgAgIMAiCSILAB+kQAIAKhCBEAkQRkC4I7Q5w0AQFQIgEgCAiAAAA4CIJKAAAgAgIMAiCQgAAIA4CAAIgkIgAAAOAiASAICIAAADgIgkoAARLuAMAAAEv9JREFUCACAgwCIJCAAAgDgIAAiCQiAAAA4CIBIAgIgAAAOAiCSgAAIAICDAIgkIAACAOAgACIJCIAAADgIgEgCAiAAAA4CIJKAAAgAgIMAiEr7oqSXJW2U/ZK/mfH8PpIuTj/fImmMpCOLfA8CIAAADgIgKu1rki6TdJb8A2CtpF3p54+S9IRsGHxzEe9BAAQAwEEARJxkBsB9JG2SDYGe/SU1SvptEfMlAAIA4CAAIk4yA+Ah6ceOzZjuRUkPFjFfAiAAAA4CIOIkMwAen37s/RnT3SNpWI757C+7sHjlQBEAAQDoRgBEnGQLgO/LmO5eSUNzzOfi9Ot6FQIgAAAWARBxEtYpYI4AAgCQAwEQcZLtIpB/OY+9TlwEAgBAIARAVNqbJH08XYykc9N/fzD9fK1s4PuW7DAwj4lhYAAACIQAiEo7ST799SQNSj/vDQS9SVKrpLGyQbAYkQXAE64mAAIAqg8BEEkQeQCcuZYACACoHgRAJAEBEAAABwEQSUAABADAQQBEEhAAAQBwEACRBARAAAAcBEAkAQEQAAAHARBJQAAEAMBBAEQSEAABAHAQAJEEBEAAABwEQCQBARAAAAcBEElAAAQAwEEARBIQAAEAcBAAkQQEQAAAHARAJAEBEAAAhxcAJxAA0Y8RAAEAcBAAkQQEQAAAHARAJAEBEAAABwEQSUAABADAQQBEEhAAAQBwEACRBARAAAAcBEAkAQEQAAAHARBJQAAEAMBBAEQSEAABAHAQAJEEkQXAL1w9kgAIAKg6BEAkQeQBcBYBEABQRQiASAICIAAADgIgkoAACACAgwCIJCAAAgDgIAAiCQiAAAA4CIBIAgIgAAAOAiCSgAAIAICDAIgkIAACAOAgACIJCIAAADgIgEgCAiAAAA4CIJKAAAgAgIMAiCQgAAIA4CAAIgkIgAAAOAiASAICIAAADgIgkoAACACAgwCIJCAAAgDgIAAiCQiAAAA4CIBIAgIgAAAOAiDi7mLZL98tm4ucBwEQAAAHARBxd7GkBZLe65R3FTkPAiAAAA4CIOLuYklzAs6DAAgAgIMAiLi7WFKzpI2SVkt6QtIheV6zv+zC4pUDRQAEAKAbARBx9zVJ35Z0tKRTJI2R7QP4jhyvuVh9+w0SAAEASCMAotq8UTYA/i3HNBwBBAAgBwIgqtEISXcWMX1kfQA/fxUBEABQfQiAqDb7S1ov6aIiXhN5AJy9bmfo8wYAICoEQMTddZJOlHSwpM9KelnSLkk1RcyDAAgAgIMAiLh7QvYK4HZJGyQ9K+mIIudBAAQAwEEARBIQAAEAcBAAkQQEQAAAHARAJAEBEAAABwEQSUAABADAQQBEEhAAAQBwEACRBARAAAAcBEAkAQEQAAAHARBJQAAEAMBBAEQSEAABAHAQAJEEBEAAABwEQCQBARAAAAcBEElAAAQAwEEARBIQAAEAcBAAkQQEQAAAHARAJAEBEAAABwEQSUAABADAQQBEEhAAAQBwEACRBARAAAAcBEAkAQEQAAAHARBJQAAEAMBBAEQSEAABAHAQAJEEkQXA468kAAIAqg8BEEkQeQCcQwAEAFQRAiCSgAAIAICDAIgkIAACAOAgACIJCIAAADgIgEgCAiAAAA4CIJKAAAgAgIMAiCQgAAIA4CAAIgkIgAAAOAiASAICIAAADgIgkoAACACAgwCIJCAAAgDgIAAiCQiAAAA4CIBIAgIgAAAOAiCSgAAIAICDAIgkIAACAOAgACIJCIAAADgIgEgCAiAAAI4wA2BnVyqEGvVFAERQBEAAABz5AuCkFdvNec/ONY1723POp6a2ztTU1pkXZq8PvY4EQARFAAQAwPGZy0eYmto6c9EL840xxvz1idndYc4t/35uXs75uNOGjQAISfqDpNWSWiXNlHRCEa8lAAIA4PALe9lKIfP52cCpodeRAIjvS2qX9CtJh0u6SdIeSR8s8PUEQAAAHGEHwH8+PSf0OhIAMVXSnRmPLZZ0ZYGvJwACAOAoJgAWUv719NzQ60gATLbXSeqU9K2Mx2+WNLbAeUQSAEcv2dK94BMAAQDVxC/EHfqfwTmfz1VqnyEAIlzvl/3yj894/D+SlmZ5zf6yC4tXDlQEAfDjlwyL9OonAACiku90b7EB8I+PzQq9jgTAZPMC4OcyHj9f0pIsr7k4/ZpeJewAuHBDU/eC39LeGeq8AQCIkrf/GvDiAt8AeOqNY3s91tbR1ev17vNcBYwolHIKuCxHAAEAqFZeaBu1eItviLvg+fk5g11m+Hts6trQ60gAxFRJd2Q8tkgxuAgEAIBqtH7nXjN6yRZjjDG7WtrNvPpG09rRczZrT2uHuWnEMrNs8y7f1//+kRnd4a+LO4EgIt4wML+UHQbmRtlhYGoKfD0BEACAEM2rb4zs1K+HAAjJDgS9RlKb7EDQXyzitQRAAABCNmrJFrN8y+7I5k8ARFAEQAAAqgwBEEERAAEAqDIEQARFAAQAoMoQABEUARAAgCpDAERQBEAAAKoMARBBEQABAKgyBEAERQAEAKDKEAARFAEQAIAqQwBEUARAAACqDAEQQREAAQCoMgRABEUABACgyhAAERQBEACAKkMARFAEQAAAqgwBEEERAAEAqDIEQAR1gCRTX19vmpqaKBQKhUKhVEGpr68nACKQA2UXIAqFQqFQKNVXDhRQgn1kF54DKAUXLzTTbrQd7VcdhXajDeNQomjDA2X34wDK4ADZlfiASlekCtF2wdB+paHdgqMNg6MNgSrHSlw62i4Y2q80tFtwtGFwtCFQ5ViJS0fbBUP7lYZ2C442DI42BKrc/pIuTv+L4tB2wdB+paHdgqMNg6MNAQAAAAAAAAAAAAAAAAAAAAAAACAoRlUHABSLfQdQxd4p6d2SXpP+Pyt0Yd4j6VD1jG9FuxWPoSFK83bZdfZ16f/vW8G6VCtve0fble4tkt6snm0f20CgitwqaYOkiZLGSTqqstWpGrdI2ilpvqTVkr4i6Q0VrVH1uVF2uXtPpStSRfaRdLOkJZJmSBov6aCK1qg63SxpWKUrUcX2kXSTpDmy6/DDskHQew5AjO0v6THZHcjxkk6VVCdpTfpvZFcru+E7QdJnJT0qu0P+YyUrVUUOkfS8pEWSUrLtifw+KWmapMmSvizpl5KmShqVfp4db36HS3pFdjuXkvTj9OMcBSzccZJmS5ok6euSLpK0WNIT6edpSyDmjpLdAZ+S8fhe2V/GHyt7jarHKEn3ZDx2r+wG8fjyV6fqnCjpDkmfl/R3SU2SPlLRGlWHAZJekvQO57HPSNotG6qR31mS7pN0suwR6E2SXlvRGlWX/SRdLelx9Rzxk6TvSVol6b2VqBSA3PbL+P8XJXVJer3z2Htlf8ktlnRBmepVbd4u++v3r+n/e792j5M0QtL9lahUlfDa6gBJH07/vY/s8jaoEhWqEl67fULSVzOe+4qk5ZLeJ44A+snc7r1D9iigJH1ItvvLlen/c+TK334Zf/+v7P7D9XNJCyW9USyHQKwMkD1NeYWkA2VX0HfL9l27R3allWx/wMdljwC+ItvBN8l+KOkI5//ehu1F2bDnPibZdh4ne1QGVmYb+jlDUqf67lSSLF+7eRcw/Fi2+8Hrc0ybVJnbvcyAt59st402SR9MP0Z46S2zDTN5y+G5kiaI9gNi4yDZjuLzJF0r2+9lluyvNcmeDmlPP79b9kjC2yV9SVKreq5uTZovyl7ckZJ0uXou7vB+CZ+Ufs4LLN5G8ChJ29PPJ122NsxmsGx/1KQHmULbzQsz90kamPFY0vlt92bK9pmUeoeUd0qaLtsnFT3yteG+Gf8+L+mq9N+EQCAGfiZ7scLb0v9/veyptrmSPp1+7FhJP1DvU0unSVop6eCy1DJePiC7Q71Rdge8Q/b0ruv/yW7wZqn3Dno/2f5sP4m+mrFWSBtmOlL2x8hPZftknS7bRzBJim23fWR30N91HjtG0lujqmCVyLXd+1T6sdc405+u3j/ovirpo5HXMt4KaUPvB/EbJK2VvTDJc7gAVNQA2V9t7sbuk5KGyB5xyeZ2Sc9GWK84O0DSt9UzHM5S2dPimTvVw2SHgblFPRu7s2Q3kB9UshXahplukLRVdsfTItu/LUmKbbdPSFon23f3cNmLk/aKC7hybfdecR7zjlS9QbadV0uaItuGSe/GUWgbSjYwr5ZtxyMkjZZtQy4IASroStkhIzL7bpwj22/oDOexD8uuvHfKXhn3zfTjST+cf4rs0YFvqO9gp2dKWiAbWl6S7Ut0lTgVl8mvDTN9WPZHR0rS3ep9hWFSZWs37++zZdfjm2X7UD4q2k0qfLvnracfkTRctq3vFW0oFbfv+Ht62mvFcghUnLdh+5jsRu1bGc8fJmmkeq5+k6QfyY4pNlnS/0RdwSrhteNg2SMDH/CZ5lDZjvgXSTq6TPWqJoW04ftld8BLZE8Fo7B2e1x2/R4tezQwKbL9iChlu3eY7NiKC5SsZS/MNhyZnnackrUcAhWTq7O8e9j+Kdm+au/MmGaE7K9dzwFKxl1ACm039/8flB0y50/q2UD67ZCTIqw29I4uvEHJGAcwrHbzuhicpJ4j9UnxVvUet8892l7Kdu/Nsv0mkyTMNtxXtg/5mSHXEYCP18negudZSQ+p99AZr82Y7iOyO4u9spfxe32J9pEd6uWaqCsbI4W222vUuwOz19H5etmBTn8o23YPK3lXrEbRhkm4fV6Y7TZcyWk31+tk+yZPll12atUTXNjuFYY2BKrYN2U7fY+WHb9qruzYS9/OmO7Pkpol/Sv9/19LWiHpVdlD+jdI2qzk3LWi2Hb7u3o2iN6pkoNkT3OkZG+Z9/Zoqxw7tGFpaLfgfiSpXtIY2e3XA7KDsf8yYzq2e9nRhkAV8zrJX+w89m7ZK/+8e6nuL3sRxxbZoUjcQ/uny169NVG2z0u+YTn6i1LaLbN/zI8kdci228cjrGtc0Yalod2Ce7ekRyT903ns7bLD33j38d1f0l1iu5cNbQhUKW+HcLjshQbe2HxeX43h6t0f41D1HsA584rUd0dQxzgK2m5ynjtL0m+iqWas0Yalod2C89rwtbIXpbn90D4me8HBaeo5GnqY2O5log2BKuWNQ5V5L0vXa2VPJ2Uexk8y2i042rA0tFtwXhtmG0rpStmLYmbInlqfKxtipNztniS0IVClvil7g/IG2ZuVS71XSvcU0RtlB4v9bFlqFm+0W3C0YWlot+D82tAvwNwoe8eJ18keOb1Ttm8baEOgqv1Ito/F47L3Rb0rz/SnStqontv1SNJ70v8maTBi2i042rA0tFtwhbSh1zaZfSTPkbRe3P2ENgSqlHe04LOyh+c/KNthd4nsWF/uNK67JD2W/vtY2SsNn1NydiS0W3C0YWlot+BKbUPX7eppzySiDYEqdbiyDwB7pKQX5X/vSsnuMF6Q9A9Jt8r26XhQvcdz6q9ot+Bow9LQbsEFaUNJeots0LlXdnzE07JM15/RhkCV+q6kNbJ9gebKjg+WOfaXJP1C0sL0v1LvowQfUM/YYBPVewDZ/op2C442LA3tFlwYbfgV2T5sm2SPnh4aXXVjiTYEqtiXZX9x/Z+kE2XHCGuVPWz/pvQ03i+5AyXdJ9u3w3vudel/j5T0hOwN45OAdguONiwN7RZcWG14kOwV1Cf///bu4EWnKAzA+BM1aSQ7ERYjWSg2/gZSNqwsNMN/wM5CKBspK4nYUDYKKwvKRlbT2LFD2bJTShaSxXtP96YRndf3zT3fPL/VnW/u1Omp0/Sd7j1n4iMeHxtKjSrfzi4SO7DPDX53lZiop1f5u2PAa2KyHyROAtg9sVGOj93ybFjHbnn/q+FTbGhDqXEPgUfddVm630r8k3hMf9B7eXB3nnhI9ydxIsBzYtf29cZueTasY7c8G+bZUGrEYeAGcJZ+Y06IsxW/0i/Vl4l8ijia5/jg3s3AOeAH8azGgQmOdyzslmfDOnbLs2GeDaVG7SCW2j8T5y++Ab7QT+R9xJ5LV7qfh8v5H4hl/mI/sAwsTnC8Y2G3PBvWsVueDfNsKDVsHrhPLNUvDD5fAe5111uAC8A3+ucxyhtaz1if+zHZLc+GdeyWZ8M8G0oz4A5wtLsuS/WXiG9j5WHeBeI80GX6yb4LeAucmM4wR8dueTasY7c8G+bZUGrccFPXMmkfAHd/u28n8I44jPsJsbT/Etg+4fGNld3ybFjHbnk2zLOhNINeAWe66w30y/Z7gZPAdWBp+sMaPbvl2bCO3fJsmGdDqWF7gE/AocFnc3+4Vz275dmwjt3ybJhnQ6lRZQl/iXhDq7gM3Aa2TX1EbbBbng3r2C3Phnk2lGbETeAasbfTR+IV/yNrOqI22C3PhnXslmfDPBtKDdsEvCd2Yv8OnF/b4TTDbnk2rGO3PBvm2VCaAS+AW8SE1r+zW54N69gtz4Z5NpQat/Hvt2gVdsuzYR275dkwz4aSJEmSJEmSJEmSJEmSJEkz5BdBUYg8HJQ9tgAAAABJRU5ErkJggg==\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# plot the first station\n",
    "plt.figure()\n",
    "data[0].plot()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "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, fig.canvas.height);\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 overriden (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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4nO3dd7gcVeHG8TeABAgdg0jovRdBQQQBqYI0RSkWQhURFRAJSAslgD86SBWkhC6dQOihV5ESaggkkkBI7+3m5p7fH2fHnZ2d2Tuze2ZnZvf7eZ7z3Ls7szNnz055d8pZCQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAoMD6SjKStow5/sal8edJ+nZg2BGlYd2V4aXxzy09XlrSzjFf2xny2jAfS3ra93itwHTmSZoo6Q1Jl0haP2Qa3dXpVxHzbqYesm3xqKSvZOt1Q8S4XpsFy4yY83op4vWDAuNtIekxSaMkzZFt51ckHdzN9BeQ9HJpmpcFhi0UMW8j6cTAuBtLukbSq5JmlsbZtsZ8e0u6QtJ/Jc2VNLZU/7Blaz9JL0iaXpr2+5IOT3ma3S2Hf6/x3mp5SdI7McfdU7b+E2TfzyhJN0taL2Rcbzn7WtLiIcNHS3ow5PmlJP1Vdp2cIruOfi1psKSDJPX0jeutz8f5phln+/GrkNcGnVwavpLvOf+y3yX7WX0q6R5JP5VdD8PeZ1Q9ng4ZP019Jd0taZhs/YdHjLezpJskfSJplsqf1eYh4/5Q0o2S3pJdJoJtVstSkk6T9LzsujFD0nuS/qLKzznM7iq3Y3B9OkrSQ7Lr3WzZz+gqSStETCvJegq0nL5KFgAvV3nl6xcY1lvS1r7yg9J4dwWe36w0vj/ELRkYZ2tJ42Q3EP7ntgp5bZioAHhpaTrbyO7UTpf0uezO5oTANLwd70khddta0jcj5t1MC8luPF+VdLXs++guAO6syvfx3Zjzekl2BxJsh3UD4+0kG8B+KWlHSXvJ7nyM7M41ynGSvlTtABhclrZW9ReRw0vTGSTpEdUOgCtJGiHpI0mHye7UfiobqnoHxj1N0vzSsN1l2/EPko5JeZph68bWkm4rvbedIt5bd+IGwEtUDvr7y76fI2XXsdmS9gmM7/+icWbI9MIC4LqybTZN0kWyy8x2kn4h6Z+yXyT80wqGuM1V2TY3KXxZ/2bIa4OiAqB/2d+p1AaPlcZ9TtISIe8zuP3yStgXzjQNkTRU0q2y4S8qAD4gu908WtL2kn4u6XVJHaXHfufIfmb3y77/JAFwM9nt+8WS9pb0I0lnyX7OT9R43RKSvlB5OxHc/n8taaDsF4btS+/jS9nPIrjuJVlPgZbUV/EDYE/ZIwDvyK5Qn3QzvrfTDu7MPd2FuKgjBXFeGxUAwzb6i0l6sjR8F9/zXgDcN2IeebGA7/856j4A1vvtNskRozD/lg3bYdaQPfq1n2oHwKhlyc/fHgeqdgAcJLtD6a5Nvicb1IJfEpo1zaAFJI2U9JnCjz7FEefz/LVs+10ZMmxxSf+R/QKymu95bzkbLHukbPnA64Lr9Tdk19eJktaJqMfqskHB012Iq7Ws1xsAo9rKO/Nxe+D5WtuvZvOvE48rOgAGPyvJhq7xpddFTTOszWpZXHa7G+RNZ+uI110re4T4fIV/vmH131rhXz7jrqdAy+qr+AHwgNK4x0oaoO5PrRUlAErSyrJHz570PVeUAOiX5wD4uOxRlKAekp6VPZ0WtcwkCYB+tQLgmqVhp8aYzkDZMNPd6ak0phlm1wTziRLn8/xY9kvfIhHDt1P5yLrHW86+I3taLXiKOrheH6Ta62aYPAVAyR61mh94TZ4CoF+tABjlBUkf1BieNABG+VFpOj8PGbaj7JHITZVsW7ag7Gnvq3zPJVlPgZbVV/ED4JOyAWMZ2RWoS/ZUS5QiBUBJelP2uhfvm60XAH8q+16CpTs9Il4XLAvGmFZccQLgGNmd1dey13HF3Wi/JHvKb5LstZjDZU8DRYWDBWTf3/KypzU7FX693NGlaX5L3QfASaU6zJU9onhIN3WuFQAPLQ07TDZ8zpRtvyGyR+f8/it7KuwQ2RA7X3b5HCB7BCvNaYa5W7Y9+3QzXi3dhZqVZd/Lbd1MZ4LsdYse/7p5hexOe03f8OB6/c/S+GvEqrWVtwD4+9JrDvQ9N1r2erSwdT7OUds424640/JLGgCXkT01f0+NcVwFQO9z2yDw/GKydT4vMF6cAOhtx/2XVSRZT4GW1VfxAuCqsjuoO33PPSd7+id47YunaAHw3tI4y5YeexuOqBJ1YbHnCNV+vVc6oyZQh1oB8BBJp8hea7aj7EZ7kuzNI8Hr6MIMkPRbSTvIXj95lWzdn1H4TugGld/jHNkLtINWlt25HFZ6HLXMLCgbRA6WDXP7y+7Ioq4z89QKgKeVhk2VvY5pN0k/k71WapakDX3jzivVc6Kk38m23/my7/+WlKcZtJxsewZvvkmqu1DjXcN7TjfT+bfs+/D4183epWF3+YYH1+unSuMHv1QFv0D5vyjlLQD+pPQa/+n80Ype52tdCyvVvukpWJLejJY0AN4pG+I3qzGOiwC4ueyXu7CgeZnslyTvy2bcALiU7KVKI1R5yjnJegq0rL6KFwD7q/oaud+Unjsi4jVFC4D3KTwA/lm2fYKlu6OA34x4XbBs0c10pPhHH2sFwDDflz2Se3GC1/j1k22jvUKGrSr7/vaQdJ3sF4jjA+MMlj3960l6qnew7NHA5SKG1wqAZ5SGvafK65lWkt0R3ex7bn5p3P0D07hStv1WT3GaQX+Sm0sTXAXAt2R3pJ7gunm67PvxtjFxA+CJqgw6/rrmLQDupfAA+JzC1/nuvjwq4nVhZdmoCURIEgC9a+2O7ma8RgPgGrLt9ZHsEUe/78t+KdrB91ycALio7LZlhqr3b0nWU6Bl9VX3AdC74PxL2Y3N0qWyosp3oIZJMwD2L702asc/XDYceOIEwH8r/BRwvTtaV6eAvbr7S9S1l0kDoGS7Sng54Ws8fUr1GRBj3H/IHknwPrMDS4+3VHmZWq40vWtKj7sL2b9U9RcTv1oB0DttFxZ+35DdOXjGl8YNHu3eQ+XLBNKaZtB7sqfx41yGUEszTgFLUi/Z+j5Tehxcr29U+Cngb6kcct6RuwC4mspf7MJ4R4e+5Xuuu7Y6VuGngBu5BjDrU8DnyL6nk2KM20gAXF32cohPZfcpfj1kQ+EdKm8jlpZ0YWl+qyq8q6FFZC9ZmqXK4OhJsp4CLauvug+A3gXntUrwmg0p3QD4u9JrNwkZtoBsML3Z91x3O4xVZL9l+u90azQAujoF3FPV3/rDNnpSfQFwuKQXE77Gs5LiHSWSbLcZRuUjnt7nX6v8pJtp/qo0XlRXKLUCoHcDQ9hO4E1V7vCfUXhY21OVy0ga0/T7bmnYBSHDkuou1EiN3QTiX6+99XU3Va/X3mdU68tZsK6NBMBeskdfL4l47TWyX0z812F211ZPlqbpDzCNBMCsTwF74e+0mNOsNwB64e9z2S8cQXHa4d+B13jhb7bsNjxMkvUUaFl91X0AvFt247aP7Lcpf/F2wBeFvC7NALiO7GmlsCNP3ukY/4ax1g7D6wamS5VBotEA6PIUcFxJA+APZN/3hXXO76+ybbRnjHHvkA273imr1VW9PO1Umt69pcdRR3g9j8ueAo46DVYrAC4oe/3j+6o8DbSybDte63vumNJ0fhGYhncdpHczRhrT9Lu29Jqo7lKSiBMAvW5grggZ5u8GZlXf82Hr9UKy12K9I3sm4cHAsO66gXEZACV7xPszVd+BvahsJ9fBjpprtZX3RS943WajRwCzOgV8luz7OTPBNOsJgKvJhr+Rqlx+/Hqoehuxg+wd9N4XRP/2cxHZO7LnyF7rHCXJegq0rL4qH+bfP6T0kV0hHqsxjbdkO/UM3rmYZgCUbMfHXbLf2PeRPbpwuuwO6bVAfbwdhr8j6D1kv+F6HUH/KTB9LwBGdQTd6B1vruyg8ufVIbvz8h57AWoh2R3YibIbzV1ku0CYIrvD85/uWlA2gPg7ZN1R9pT6UbLtspfsRtIbz38a6kbZQPkL2c5Y95e9sNuofBdflKhl5uTSdH9Ver8HqHztWLArh16+93+pykcy9lf1TsELiA/JLg8HyHZ3MUmV1+AtLNt+k2XvaN5Z0v/JfjHyH/1Ka5qS/aIyVbZz4She+8X5pYmXZC+OD1vvt/ON5+8I+mcqdwT9kWp3BB1cr/dX+ahNcL1eVzYITJVddn4iG9r3kv2SN0W2KxJPowFwW9kvDm/JhtwdZW+Serf0noJ3gr6kyo6gd5QNfo+W5vOMqo/Kj1Z0R9Bhv6yRpg1V/mzflj0l7z32/5qLd03vIFXXeStVWt43jdtKrztK1ctP2PZkBdllb5ZsN0DBeXV3d3vU5+t1zH1WyDSDnW/HXU8PK9Xf/0tGa5aeuy4wzZGyX2aAQuir2ofXvQvOgxt5v9+Wxgles5R2AOwhexTFu3ZvruxRhvNkQ4Bf8Dq6TtkV/U3V/1Nw/WvUrZleUnQdvSNfC8jeifmpbEDukN0Ah/1MUliIWFd24zpa9gvBbNmd5cmyQcbvCNlTyuNlg/Vk2e4VuvspOP+8g8vMPqVpjitNc6rsBfbBo2dS+DWTXgk78rGf7DI0RzZoPKDwnzhbTtL1sj8X1SG7oT9B4ddgpTFN76arX4cM8yxdGmdgjXE8tZabYIDcU/YLwMRSPb2fggtbb2qt168pPABK9o7NU2XXyWmyn/NY2eBwtOzROU+jAVCyR88ekF1OO2WXrXsVfrdrsK1myC5L9f4U3Mga9UpDrcst/Kd5ay0TwUtVam0f/ctP2Paku21rd6eeo44y15pm2JeiOOupd4Q37IxS8GzLaCXvYxEAgIbtLXsEMeyaXAAAALSgS2V/7xUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACZ6CHbe/qSFAqFQqFQClX6KLxDcaBbfVS7F3QKhUKhUCj5Ld39BB4QaklJZtSoUWbq1KkUCoVCoVAKUEaNGuUFwCUzzhEoqCUlmalTpxoAAFAMU6dOJQCiIQRAAAAKhgCIRhEAAQAoGAIgGkUABACgYAiAaBQBEACAgiEAolEEQAAACoYAiEYRAAEAKBgCIBpFAAQAoGAIgGgUARAAgIIhAKJRBEAAAAqGAIhGEQABACgYAiAaRQAEAKBgCIBoFAEQAICCIQCiUQRAAADqcfXVxgwcmMmsCYBoFAEQAICkRo0yRrIlAwRANIoACABAUu+/Xw6AXV1Nnz0BEI0iAAIAkNTQoQRAFBoBEACAJL76qhz+JGPmz296FQiAaBQBEACAJI49lgCIwiMAAgCQxDHHVAbAzs6mV4EAiEYRAAEASOJ3v6sMgPPmNb0KBEA0igAIAEASBEC0AAIgAABJBANgR0fTq0AARKMIgAAAJBEMgHPnNr0KBMDW9jtJ70maViqvSvqxb/hzsh++v9yVcB4EQAAAkjj6aAIgUrWXpD0krVMqAyR1SNqwNPw5SddLWsFXlko4DwIgAABx+TuA9sqcOU2vBgGw/UySdHjp/+ckXdbg9AiAAADE1bt3dQCcPbvp1SAAto8FJR0oaa6kDUrPPSdpvKQJkj6QdJGkJRJOlwAIAEBcwfBHAERKNpY0Q1KnpCmyp4Q9R0raWdJGsuFwhKSnupleT9mFxSt9RAAEACCesAA4a1bTq0EAbH0LS1pL0paSzpc94rdBxLhbyC4M36kxvf6qvnGEAAgAQBxhAXDmzKZXgwDYfp6WdF3EsB6yN4kcUOP1HAEEAKBeBEBk5BlJN0cM20h2YfhhgulxDSAAAHGFBcAZM5peDQJgaztP0naSVpO9FnCApPmSdpG0pqQzZE8NryZ7beBHkv4je8NIXARAAADiCguA06c3vRoEwNZ2o6SRsnf+jpM9/btLadjKkp6XNLE0fLikyyUtm3AeBEAAAOIiAKJFEAABAIgrLABOm9b0ahAA0SgCIAAAcYUFwAz2oQRANIoACABAXARAtAgCIAAAcYUFwClTml4NAiAaRQAEACCusAA4eXLTq0EARKMIgAAAxEUARIsgAAIAEFdYAJw0qenVIACiUQRAAADiIgCiRRAAAQCVXn7ZBpsrrsi6JvkTFgDHjWt6NQiAaBQBEABQyR9uUCksAN56a9OrQQBEowiAAIBKBMBoYQFQMubLL5taDQIgGkUABABUIgBGiwqA22zT1GoQANEoAiAAoJI/2Jx7rjHf+lYmXZ3kUlQA3GKLplaDAIhGEQABAJXCAs6ii2Zdq3yICoDf+U5Tq0EARKMIgACASlEhB9FtIxkzYULTqkEARKMIgACASgTAcE8+WTsAHnhg06pCAESjCIBAK5o505j338+6FigqAmC4WuFPMmbDDZtWFQIgGkUABFrRuuvaHdLjj2ddExQRATAcARAthAAItCJvh/Sb32RdExQRATBcdwFQMuaee5pSFQIgGkUABFqRtzM65JCsa4IiIgCGixMAJWM++ST1qhAA0SgCINCKCIBoBAEwXNwA2IS2IgCiUQRAoBV5O6G+fbOuCYqIABiOAIgWQgAEWhEBEI0gAIYjAKKFEACBVuTthA49NOuaoIgIgOEIgGghBECgFREA0QgCYDgCIFoIARBoRQRANIIAGI4AiBZCAARakbcTOuywrGuCIiIAhiMAooUQAIFWRABEva66igAYJScB8MOvpprrnx5KAERDCIBAK/J2QocfnnVNUDQZHtXKvZwEwFX7DTIrH3cPARANIQACSXV1GXP77ca8+27WNYlGAES9CIDRCIBoIQRAIKknnsj/DtGr3xFHZF0TFA0BMBoBEC2EAAgkdcEF+d8hEgBRLwJgNAIgWggBEEjq/PPzv0P06nfkkVnXBEVDAIxGAEQLIQACSZ13Xv53iF79dtst65qgaAiAxnR2GrPVVsYccEDl8wRAtBACIJDUgAH53yG2404bbhAAjXnllfJ7nj+//HySAJhiJ+wEQLhAAASSOvfc/O8Q23GnDTcIgMY8/3z5Pa+3njGzZ9vnkwRAyZiU9q0EwNb3O0nvSZpWKq9K+rFveE9JV0qaIGmmpIclrZRwHgRAIKlzzsn/DrEdd9pwgwBozLPPVr/3kSOTB8ApU5xX7dH3viIAtoG9JO0haZ1SGSCpQ9KGpeHXSBotaWdJm0t6VtI7khZMMA8CIJDUWWflf4fYjjttuEEANOapp6rf+5FHZn4EsHN+l1m13yACYJuaJOlwSUvJhsEDfMNWlDRf0m4JpkcABJLq3z//O8R23GnDDQKgMY8/njzshZVp05xWq6NzPgGwDS0o6UBJcyVtIOlHsh/8MoHx3pV0VoLpEgCBpM44I/87xHbcacONIgTATz4xZs01jbnhhnSm/+ijBEBkbmNJMyR1Spoie0pYkg6WDYNBT0q6rsb0esouLF7pIwIgkMzpp+dvhxiUx502iqEIAXDHHdOt0913uwmA06c7rRYBsL0sLGktSVtKOl/SeNkjgFEB8ClJ19aYXn/ZBaaiEABbzBVXGLPxxsaMGVN+bsgQG1zmzcusWi3jtNPyt0MMyuNOG8VQhAD4/e+nWyf/jV6NlBkznFZr7jwCYDt7WvYIX72ngDkC2GrmzjVm2LDK57yNz1FHVT937bXNrV8r+utf87dDDMrjThvFUIQAuPzy6dbJf5lHI2XmTKfVmjOvkwDYxp6RdLPKN4H8wjfs2+ImkPbjfRMeNKj8nLfx+c1vqp/785+bX8dWc8op5fYcO9Z2GdHVlXWtKuVxp41iKEIATLtOe+5JAESmzpO0naTVZK8FHCAb8HYpDb9G0ihJO8l2A/OM6AamtXV2GvPWW/avx9vQ7Lln9XO//nX1cwTAxp18crk9F1/c/r3//qxrVSmPO20UQ7sHwCeecBP+JGNmzXJWra6uLnPrKyMIgG3iRkkjZa/1Gyd7+ncX3/BFZDuCnihplqRHJK2ccB4EwCI54QS7UTnmmPJzYRtC7/GvflX93Akn2G+lm29ugwyS69evekN/+OFZ16qSv24PPWTMuHFZ1whF0e4B8LDD3AVA7xdEHPA6gCYAwhUCYJHUCntJAuA11+Rvg14kJ51UvaE/7LCsa1UpWL9VVsm6RiiKdg+AUad/998/0wB4xoNDCYBwqngBcPx4p4fVCyVpAPzlL6ufO+EEe5dw3jboRfKXvxQvAPJZI652D4BR732ffTINgBue8TgBEE4VKwCOHWtXqqWWyrom2XBxBPD44wmAjTruuOoN/aGHZl2rSnneeSPfagWak0825kc/MqajIz91THPajRaHAXC90wYTAOFUsQLgffe1987Mv2G59FJjHnuMAJiFFVckAKJ1xQk2Dz6YnzqmOW0CIFpYvgNgV5cx8+eXHxMAa5fgeFEB8PLL27sdG7XQQtVt37dv1rWqRABEveIEmwceyKZup59uzBJLtGUAPPaO/xAA4VR+A2BXlzE/+IExG21U7vaEAFi77L+/MVddVTsAHndcZQB87bXs3k9RLbhgddsfckjWtbKdgm+3Xfhdyu26ziC5OMHm0UfzU7dmzMNhAHxjxERz2VPDzLzO+SEzj8ZNIHAtvwGwo6O8Enm/dEEATFaiAuBvf0swaESPHtVt7e90Oyu33VZ7eQDiiLNtefLJ/NStGfNwGAC9AHfbayMTVcsf/giAcIEAWCQuAuCf/mTMj39MMGhEWFvnIQDecAMBEI2Ls2157rn81K0Z86i3hPRY4QW4sx/5IFG1CIBwLb8BcO5cAmCQiwD4xz+2VwA880xjbrnF7TTD2tr/qytZue46AiAaF2fbktVNIC0UAM96mACIbBEAiyKs77l6AuCxx7rdgD72mDGrr27Miy82Np16jRljzMUXGzNhQvWwN99MZ3nJawD0d/BNAES94mxbsroJJKwuH35ozDPPpDsPAiBaUDEC4Cef2OfaOQDWs/EJC4DLLus2GHjTWHjh8OG33mrMxhsbM3x4Y/OJsvHGdv677VY9zP+bni5MmWLM669339ZZOftsAiAaF2fb0q9f/ur2zjvpz8NhAOz/8PuJqkUAhGv5DYBz5hAA/erZ+IT9EkhY6epyU69aw3faqf55xJ3/yy/bgOZxHQBXWSVe2M5KWLgnACKJ996Lt23J6oh3d/X68svK8bu6jNl7b2N23z3+dq5JAfDMhwiAyFYxAiCngNMNgF43O43Wq9bwLbaofx5x5++VuXPtsCefdLu8xG3rrCyyCAEQjYm7bTnrrHzWb++9K8efMiU6HNY7DwIgWkR+A+DMmeWV6OOP7XMEwGRl//2NOeggYwYOrD3evHmN12uBBWoP33TT+ucRZ/7+4i0vzQyAu+ziZh6NWHhhAiAaE3fbctdd+azfd79bOf7kyeVho0e7mQcBEC0ivwFwxIjySvTqq/Y5fwD88ENjLrjABsV24HKjFCz/+lfj9VpoodrDN9ig/nnEmb+/vPCCHfbUU80LgHkIWGEdVOepfsi/uNuMp5/OZ/2CAXDSpNwGwDMeHJrorRMA4Vp+A+Bnn1V+2zzrLGOuv756BfvLX7KuaXOkGQDXXLPxekXdBOINX3/9+ucRZ/7+ctllxqyzjjEXXtheATCsg+qw+k2YYMzXX2dXT+RX3G3GJpsYM2pUPut3zDHl6/38ATBnp4BPTxAA/+/xjwiAcC6/AXD48Hgr2LbbZl3T5kgzAH7zm43Xa9FFaw9v5hHAsNLIjS5x55W1OPXr6io/njEj2/oif5JuOwYMyGf93nvPjj9xYvm5L75wO48GA+BpD8QPgMHwRwCEC/kNgJ9/Hm8F22abrGvaHGkGwOWXb7xevXrZx/fdZ4/KjhljA4Y3vJlHAMNKI9c5xp1X1uLUL6x7JbSvV16xN0543TTVs/1oprh18n7j3B8A77jD7TwaDICnPvBerOqMmzaHAIhU5DMA/ulPNlDEWcG23jrr2qbroYeMeeONdAPgN75Rf/28aSyxROVjyZiePSsfd3S4aZOw+XdXvDuD05xX1uLUzx/K33472/oie96ysMkmxgwdmk0AnDPHmCOOMObRR6uHBb+4JQ2AEyaUnxs4MF59chYAd7xwCAEQqchfADz33GQr2FZbZV3j9Hz0UbrBz8VG3Hv9UktVPg4rafQdFvf9hfwou/N5ZS1O/fx3Rd58c/m1AwYYc+ON2dS72To64p8ObHVZbjs8p5wSPq2TT7ZfTj/8MHl9vRsH/QEw7k9CNikA/vX+eAEwLPwRAOFCvgLgrFnJV7Dvfa/59Rw+vHynaZoefbQ4AXCZZSofNyskxX1/06enP6+sdVe/yy83Zq+9yo//8Af797LL8vMemmGbbex7bcY6nHdZbjs8664bPi3vuZ/+NHl9X3nFjj9+fPm5G24w5t13jZk/P/028UqNAHgKARAZSy8A3n+/7aYlrs5OY37xi+Qr2JZb2tePGWPMwQent1GfPNlO238R/QfJfssxsccfdx/0Gt2IB4+kea9fbrnKx1HlppucNU+s+XnlmGPSndfqqzc+/TTr16ydeRF479X/6y233GLMoYe6uVa0SFwsMyee2Fgd1lqrevkbPLj8nHd9cZL6PvusHX/s2PJz3/iG/Xvaaem3iVdqBMCT7yMAIlvpBUBvBXjppXjj3357fStYr172J8D23tv9Tmz8eGPWXtt2QbP66nbat95ans/tt7ubV5i8BcD337fjHnlk+Tnv9b17Vz5uVshoZriJM4977ml8PvU6+ODGl4N26FfTe6/LLlv9nP+0eDto5vYjymqrVU5n3rzq6Xd12ZvL4tbnwAPttL7+Onl9XW5XawbAd2M1DwEQaUk/AMbtZHiPPfKxMfI7+eTqaftPobVbAPQHDGMqv117dxI38/OJOz9X82328pfUn//c+HLw2WfZ1T9Ns2cb85OfGHPlleGflfe4la8pDpOHbe5KK1VOx/8zoF75v/9LVp/dd4WXEtMAACAASURBVLfTGjMmeX1dblcJgMix9ANgd0dEPv3Unp5zudK5ssEGtecTt1uBeuUtAP7855Xjh02jmZ9P3Pm5mm+zl7+kkhwhiSqt2jXMNdfU/qy8x2uvnV0ds5CHbe63vlU5ndmzG6/Pzjvbo4ZhX+K7q+8Pf+iuXWoEwH73vmvGTZvTbfMQAJGW9APg3XfHG89lccX/zTSs/OlP7uYV5okn0mmfettszTUrxw+bRjM/n//+t7nLhcv5PP20MY891nid/Pr1a3w5eC/edUmF8e67NgTsuGPtz8p7vNFG2dU1C3nY5ganU8/NgMGy/faVvwVeq75dXZXXNu+8s7t2qREAvXLfW7V/UYUAiLSkHwC7+9FwVyuayx29Z511up/Xrbe6m19Q3gJgcPywaTTz82n2cuFqPv5rnCZObLxenr/+tfHl4PXX7Q1ZrSLOe/7978v/b765my6DiiIP29zgdFwEwG22Cf/p0LD6HnaYfe6jj+zjXXZx1y4h19QGg9yW5z5Vs3kIgEhL+gHwzjvjjeeyuLLJJvHmN3asu3n6FS0Afvpp/HmdfHLjv+Pc7OXC1Xz8v8YxbFjj9fKcdlrjy8F++xmz2GI2CHomTrSXchQxGNXbDmeeWXu6X3/dnBtmRoww5uqr02v7PGxzg9OZObPx+my5ZfRp/x49wud/xBH28a67umuXGAFwi3MIgMhG+gGwu+vkXK1orjZGflttFW9+K6zgbp7G2CMwQ4fm6xrAf/yjevzgNJZfPt68pkwp/z9uXP3tVM/7bOTXL1y1ZWdneVzvqIMLq67qbnnYeOPydL/zHfvcH//orq5hurrssuFSGuvEl1/a4Usv7bauYRZayM6ru65L6tWs7UeYDz+svvnPGNtnZ5rbugUXDG+Do46yvyCy5Zbu5hUrAD5Zs5kIgEhL+gEw6k7Zrq7KPvVcFleirhtKc57GlDtG9V8cnXYxxvbbePLJlXXx+nNcYIHK8f0hJmnx/z7n6NHleQ0dmuwXGpq9fCSZ/rx59mjahRdWT2f+/PK4779ff308Dz1kzE47uV0eVlut+n1/85uN17UWbz4u+/JMYztyxx3prPejR9sjj199FV5/z8iR9ohsdx0ax+Fy+5GU17WWv5x6qjHTprldloNloYXC28B/jbOrEiMAbn42ARDZSD8A3nZb9bDOTmM228yYH/0onRXclSRd07iU5sYvqvjDuNcViP9nlFyWSZPK/48qXQA9alT5uY8+MuaEE+xptjTaKe3P5fDDjfnb38qPH3vM/pqBd7TT39Yufo83rWUiOH2vs2+X5syx11z528vl+uTi/Qe5DoDbbmun5X3x23hjY4YMMea55yrrc9dd5b44JTd9FrpeVlzM+7330lumpfJvn99/v9ujfWElRgBctd+g0OaZMWeemT+/iwCI1KQfAMN+gLveHx1Pc2MUNGxYsnm6vEYnzbaJKv4jeoMHV4Y018V/CvjZZ23Q9l+/tuCC9u8ee0S3USP9Rv7sZ/ZXH5r9ufzyl9XTOuWUxm66eO219Ncj77H3c39BxxxjzOmn11f/a6+tPe9GuXj/Qa4DYL318/+SSbPn7eLz6tEjvWW3Vll4YWMeeaQ586ozAH4xcaZZtd8g86sbXiMAIjXpB8Cwu2SHDEl3pWvEF18Yc/bZ9c3X/wsZjchiozhuXPn/dde1O/Qs6uEvffqk20ZJf/bL5bLpf27zze1zc+ca8/zz9m8elpWZMytPVweveZs/35jf/KY8vJ4bWvx34Ia10x132JDc1ZV82o22z5VX2nkH3XlndT3r1ehlMJ7XXzfmvPOM6ehINn9Xy0o9gpeVNKv07Nm8edUZAC996pPI4EcAhCvZBMCnnkp3pXvgAfsD943Uu97S3V3PzahDPSV4Pc4JJ2RTj2BJs41OPDFZsHD5fsKeP/xw+/9RR+VnWdlss/L/Sy5ZOe+BA6vHf/XV0L7P/mfcOPvrDHPm2GsXa302L7xQfvxU7Tslq0yZ4qY7Ee89GWOv0bv2WmNuvLH75TOu8ePdLk9//3uy+btaTiZMSP7evSP9zS7ejTXNKDED4Iw580zn/C7z5oiJ5qyHPzADHv2QAIjUpR8Ab7mlethddzVn5Yt7fVVHhz0S8eabbub7xhv1tdnPfmbM1ltns1EMFpe94TdSwmy4obvpDx4c//NxMb/NNrOdQIe9z+7ed1C9R6pdfR5nnFE9fIkl7N8//MGYgw4yZpFF7E0/3s0N9c63u/5E/bybCLy6NFrWWcdOt08f+3i55ZJ9TsG6DR1aDn777ttY3bxrrL3Hhx2WrD4ul42k3WE1M4hlVWIGwFX7DTJ7Xflit6GPAAiX0g+AYRcqX3ddc1a+J2vfXfU/Ln+H2CtJrwlM647oektaN+gkLd7RF89LL7md/k03xf+M0nyf/usipe7r4j9l38zi8d+M0IxyzjnxP6dXXnE//3vuqd0eWS5D/tP0Xp1GjrTtEGX+/MoOyZv9GaXVFnkrCQJg0kIAbG2nSHpT0nRJ4yQ9KGndwDjPyS4A/nJXgnmkHwDDdrAXX9yclS/u6dg05v3mm8nay8XvX7ZqMcb+Ru122xlzwAHupx/3SHGz33MtAwZk81kcdZS9Ns5ln4NxywMP2CNo3Z22P/ro/HxOzViGgl0y+ecT9tN+Tzxhj45GXX/ZSDn1VHvE9N57s2mLvBUCIOr0uKS+kjaUtKmkQZL+K6mXb5znJF0vaQVfWSrBPNIPgP/8Z/WwsFNHaZT1109WV5flzDONOessYyZPrj1v74aYc87JfmOV12JM5Wm3tOaRxXISVYYOzU9d8lK861T3379226yxRvPqVEtXl+3U2BhjLr88vaOmb79dXSfvf//298EHbYhOOv0ZM9y3TbsswwRAONJb9sP+oe+55yRd1sA00w+AN95YPWyttZq3Ag4c2P2vTTSrLlnOu8ilGafHJ08OP23f0WHD+ciRzX/ftW6myPozybpECZ4ObXY9urpsp9znnWdvmsmqbbz/vQBYb4jzpte7t7vPqJ2WYQIgHFlL9sPeyPfcc5LGS5og6QNJF0laosY0esouLF7pI6UcAG+4ofzctGnGXHRRdhvFoCFDjPn2t41ZdNHm1GHYMPsN3L9jz3oDVYQS9yfm0lhOsnzfO+8c3Udg1p9J1iXM+ec3vx5vvGHvah4zxl7uct552bfNQQeV/7/2WrvNifs73VFtXU+n8C+9VP5s7rnHdi3lnb7P2zXPaRUCIBzoIelhSS8Gnj9S0s6yofBASSMkPVVjOv2lqmsG0w2A/t9d9PcZ1uxSq47NKiusYP+GdQhMyUfJehkJK8GL67O6ASRPJQ/rczsUf0f+SV/bt2/1a5991j6u51R0EQsBEA5cJWmkpJW6GW8L2QXiOxHDm38EULKPs/7G5zdrVvNPFQXLAgvk505bSrl4d0caY8ymm2ZfH68MGmR/wSDreuSlTJtWXp+PP96Yk07Kvk6tVh59tHK7GbzWsLvy/e9Xb/svv9w+d+qp2b+/ZhQCIBp0paRRklaPMW4PSR2SDog57fSvAZTs41VWyXZFfOQRY3bbzfZVlfVGgZLf8oMfZF8HSvxyxRXZ16EVS9Sv8NQzre99r/LxRhvZO/qzfo/NKARA1KmHpL9L+lLS2jFfs5FUdaNILc0JgK77m6JQKBRKeiXqZ/2yrlfRCgEQdbpa0hRJ26uym5dFS8PXlHSGpC0lrSZpD0kfSfqPpAVjzqM5AZBCoVAoxSlRDj44+7o1s5x4YmOvJwCiTiai9C0NX1nS85ImSporabikyyUtm2AeBEAKhUKhVJYohx2Wfd2aWcaONWarrep//YwZBEDkFgGQQqFQKJUlyu9+l33dmlnGjbMdj9f7egIgcowASKFQKJTKEuWEE7KvWzPLjBn2buhGXk8ARE4RACkUCoVSWaKcckr2dWtW8fosbGR/RgBEjhEAKRQKhVJZopx5ZvZ1y6INVlqpvmkQAJFjBEAKhUKhVJYoZ5+dfd2yaIPrr69vGiEBcKsBTxMAkQsEQAqFQqFUlijtGgAHDqxvGiEBcPfLXjCr9htkvjfgKQIgMkUApFAoFEq5HHxw9Ha9XQPg7bfXN40aAfDhd740x97xHzPk47EEQGSCAEihUCiUcvnVr6K362edlX39mlX8Bg2qbxo1AuALw8b977nx0+eYT8dOM5+Nm/6/gLf2qY+ZVfsNMpv0f4IAiFQQACnNKddc0/04119vTEeHMS+/bMxCC9nnvv/97OtOobRT+fWvo7frrR4Ar7vO/j377Mr33dlpzDbbJJ9eSADc7dLnqwJgmMkz55qXPh1vpszsMPte9ZK56aXPzfF3v202PesJAiCcIABSmlNqLROff27M4MGVy09Hh+2F/8Ybs697nssHHxizxhrZ14PSOmX55aO36xdckH390ii9exvz2GO192n//W/y6dYIgC8OG59wp2pMV1eX6ZzfRQCEE+kFwKWXzn6lpuSnGGPMmDGVz+27rzH33Vd7OZo2rfl1/fprY4YMyb7NosrAgca8845tG2OMGT++cvhbb2VfR0qxS5QRI7KvW7Pea9Dcucmn7TgAegiAcCG9ALjMMuWV4M477d8ddsh+Zc9LOeSQ7OvQzOLp3ds+vuCC+MvSqFH2aOC8ec2ta9ZtFiyffWbXpfnzq9vIa5tdd7WPOSpIaaREmTkz+7rVU954I/z5BReMvx0yxpgHH0w235QC4HqnDSYAomHpB8ALL7SPZ860O65evbLfGGRZXn3Vtsf06cYccYQxTz6ZfZ2aUTwTJhjzyCM2sNSjmXX95BP7eN117Snp+fONGTrU/p92HXbd1f59443wwNed1Vevb749eoQ//8ADxpx0kjHf+Eb2yxIl3bLQQrWXrS++MObLL7Ovp1eMMWbYsPBhM2caM3myHcd/OcmJJ9q/V1yRbL16+eVkdUspAG557lMEQDQs/QD40UeVz3d2GjN8ePM3EvX24+SyBK9z82Rdr2YUV9Kq3777GrPoosY89FDl/Lq6quvQ1ZVePbbYwpjf/97Op57g51lttWTz3W238vy859ZYw5illjLmX/8qT9cLxZT8lqOPtn9fe62+wH7oofGWsTjT+uc/7biffmofb7CB/TKx1VbGvPmmfW6XXcJfe/DB8eZRqz5+XV3GbL+9Maeeav//4ovk65VX57iFAIgca34A9DRrYzh0aHTdmllGjYpuq2bXJYviSlr1e+SRZPU44oh8t9Oqqyabr/8CeO8573RyUJoBOGl57TVjTj/dmDlzjLnttuzq8dln2beFV776qvLzSvr6Qw6Jt4zVOh16wAHV40+aVH3kf9Ik+8Uj+PohQ+zwOXPi1zVsuGvvvpusLVMKgP8eOclsftqDRiIAon7ZBcD11mvOxjBMZ6e93f+LL5pTh1mzarfVbrtlu8PYfvtsPod6pFW/Z56prz71/kRUWFlvPXfttN9+3X/mv/ylXRdGj6587V/+YseZOzd6+p2ddrn2OgY+5pjmLrO9ehlzyinV9Wr2uiMZc+yxdt5TptgbcsaONWbbbY156aVs6hPcnk+caANV2JeWXr3sEUN/mDvzzGTL2rnnll975532CFvSSzy8159/frLP9Ljjqsfbf3/7d7XVktUhjqRBP6UAaIwxU6ZMMRIBEPXLLgDOmmWHp9WvVBwzZqS7IU5q8GD7ukUXbe4Owxi307v4YmMuuaR8nc3llydviyhbbpm8PmuvXa5LsBx5pDFbb13/NYnG2OuM6rnB6frrjbnnHnv90l/+Yu8+dmXcOGOOP94eAXexbEbp6rJ3ho4dm/5y+ve/278jR3Zfr2asN198YT/Dzs7oepx7rjE/+1lz6nPDDfYoaBTvaNo++9h1crHFKtvy6aeN+d3v7PKcxH/+0/hy1bOnff2wYdXDar3n004rjzd6tL1mtqPDmLvvtsuka7NnJ/tMUgyAU6dONRIBEPXLLgD6udj4bbyxMffea69fibODMCbdU1mPP15fu73zjr1BpBk7DP9G+7XX7P9HHpn89cE7TtM0enS8On3nO/bvTTeVX+sN+8c/bKeu997rrl577VVfuzdDM+addMdY73LayHtupIwfb8yll9p1c9ddjbnyymT18a5BvvRSt/XacUf79+ijk9XHtYED7RHPek2bZvsDDfOvf1Vf0jBggDGbbFK+waOZan0ed99tj/56jwmAyLHWCYBhF+vH0dXlLggedJD96w8d9Up7Z1prx5r09z6D9U3b7NndB66uLnt9kV9nZzpHBYwx5uabG2/3tBx9tDF77JH+vIcOjT7i2EjZfXdj/va3ZHVJex2ph3eTTWenLd6d2ksskbw+jz5qr9vs7DTm7bcbu2GoKLbbrvnrTphan8vIkZVd5UyfXvVyAiDyovgB8Oijo++uTcLFTmLiRNvNiQuud2CSMbfeWv3cySeHz7+ry/aRFXfn6HpnGcd776W7w07illuyCRVJePPt0SPd+cyaZQPKww+7WW7r8cEH5e1DVnXozvz59qaC+fONufpqG+hq1eHtt22XRC+8kE598s67BCCLdcfPq8PLL1f/xOXUqeXLmyQCIHItHwHQO42RtOyyi7v6jhtnLyofN852BRK1MfZ+4eS552wXGQ89ZMwrr5TvWnNlxoxy1wmuStiRzjFjouswdqy9w9Lf+/+oUZX1uugiO+5vf2sfb7KJ23boTlg/iuus09w6GBP96xu//rUx115bfjxggL3Q3fsVj2by6rDWWs2b5447GvOHP5TnvdRS8ZbVF16w3YB88knjdag1n8UWqz38+efr6y6kXtddZ8zmm1fW4ZNPjHn99ebVIa+eeCIfAXDyZPvl0/PRR3Y/4F32QwBEQeQjAHriBpnPP7c3j9S6+NqFBx6wNx2ccELlhqfW3ZGuuQyAYdOrFQD9hgyp7Cpl9uzq0+5ffVX/qfhG3HOPPZ3273/bYNXR0fw6GGPbx39Usndv+7w/eDczTAS98IIxO++cbJ10Zdgwe5Rr7lw7/1mzjPne92ovq64cfrgxyy4bfrf94MH2Ro4XX0y/HknkJezkSVeXvXnl2WezrkltBEAURPoB8OOP478mSZBppjlzbNcCN9zQ/HlPnGh3lEl7oHcdABGf9yse/v7YHnig9h2a7aizs/pO/OHDa/eZ2ci85s61R2qmTCnP78UX7fCwG6+WXtp9PZKYPLk9ru1rNQTAtvPHBCVP8hUAg6c+gqV//+pfamgnaQRAl12PoIwdd3z77GOXxc02a948l13WztN/Kv755+2drF9/bcw552R3JBnFRgBsOyNils+zqmCEfAXAMWNs32Uff2x/d/SPfzTm/vvttFzc6FF0TzxhzB132GsVCYBoFZMn299kbeayOGdOdYfJgAsEQBREegHQu1kiSQCMksV1ZXn32Wf2mqZbb7V9kvlD3Q9+YPv1e/pp+79kTyMbQwAEgDQRACFpYUnrSloo64rUUIwAiO795z/2TtwPP6wMzF1dtksM75RkMACm1S8eALQjAmBbW0zSjZI6S2WN0vNXSDo5q0pFIAC2m+efJwACQFoIgG3tckn/lrStpBkqB8C9Jb2dVaUiEADbkT8AcnodANwhALa1/0rauvT/dJUD4FqSpmVSo2gEwHYUvDEEAOAGAbCtzVI59PkD4KaSpmZSo2gEwHZEAASAdBAA29rzkv5Q+n+6pNVL//9d0uOZ1CgaAbAdEQABIB0EwLa2jeyp3mskzZZ0maSnZK8H3CLDeoUhALYjf8faAAB3CIBtb2NJt0h6X9KHkm4rPZc3BMB25G2czjor65oAQGuJGQBf+pQAiGwRANsRARAA0kEAbGtDJB0uaamsKxIDAbAdeRuns8/OuiYA0FoIgG3tCkljZK//u0/SvrK/CpLUKZLelL2RZJykB2V/WcSvp6QrJU2QNFPSw5JWSjAPAmA72mor+9l8+mnWNQGA1uIPgNOmVQ0mALa+BSTtKulm2a5fJkm6XtL2CabxuKS+kjaU7UJmkGwfg71841wjabSknSVtLulZSe9IWjDmPAiA7aiz05hJk7KuBQC0ntmzCYD4n0Uk/Vw2mM1vYDq9ZT/sH5YeLyWpQ9IBvnFWLM1jt5jTJAACAOAKARAlK0g6Tvan4bokvd7AtNaS/bA3Kj3+UenxMoHx3pV0VsQ0esouLF7pIwIgAABuEADb2pKSDpXt+2+epE8knSkb4OrVQ/b6vhd9zx0saW7IuE9Kui5iOv1lF5iKQgAEAMABAmBbmy3pK9kOoL/raJpXSRqpyhs8ogLgU5KujZgORwABAEgLAbCt7Sp7E4grV0oapfJPynnqOQUcxDWAAAC4QgBsewvJ3pn7W0lLlJ5bUdLiCabRQ/b3g7+UtHbIcO8mkF/4nvu2uAkEAIBsEADb2qqSPpLtl69T0hql5y9T9KnZMFdLmiLbdcwKvrKob5xrZI8O7iTbDcwzohsYAACyQQBsaw9KGijb+fN0lQPg9pI+TTCdqps1SqWvb5xFZE8RT5Q0S9IjklZOMI/0A+Ann7ifNgAAeUQAbGsTVP7FDn8AXE02pOUJARAAAFcIgG1tkqQNSv/7A+C2ksZmUqNoBEAAAFwhALa1u2V/9k2yAXB12Zs/npF0U1aVikAABADAFQJgW1tRtvPnD2U7gn5V9rTwx5KWz7BeYQiAAAC4QgBse4tKOky2G5erJR2hyrt38yK9ALjUUgRAAEB7IQCiIAiAAAC4QgBsO3snKHlCAAQAwBUCYNvpilnmZ1XBCARAAABcIQCiIAiAAAC4QgBsW9+QNETSOllXJCYCIAAArvgDYMi+lQDY2sZLWjvrSsREAAQAwJU5c2oGwF0vIQC2soslXZB1JWIiAAIA4ErMAPgyAbAlXSlpqqS3JF0n6ZJAyRMCIAAArhAA29qQGuXZDOsVhgAIAIArBEAUBAEQAABXCICQtJak3VT+CbgeGdYlCgEQAABXCIBtbTlJz6jc8fMapedvlL1BJE8IgAAAuEIAbGu3Snpc0kqSpqscAHeV9EFWlYpAAAQAwBUCYFv7WtKmpf/9AXB1STMyqVE0AiAAAK4QANvadJU7gvYHwO9KmphJjaIRAAEAcIUA2NYelXRO6f/pskf+FpB0j6R7s6pUBAIgAACuEADb2gaSxkkaLGmupH9J+lD21PCaGdYrTPoBcNgw99MGACCPCIBtbwVJZ0kaJOkxSedK+namNQpHAAQAwBUCYFtbRdF9/q3SzIrEQAAEAMAVAmBbmy9p+ZDnlysNyxMCIAAArhAA21qXpN4hz68qaWaT69IdAiAAAK4QANvSJaUyX9K1vseXSLpc0muSXs6sduEIgAAAuEIAbEtDSqVLNugN8ZUnJF2ncv+AeUEABADAFQJgW7tJxflACIAAALjiD4BTplQNJgAiL9ILgEsuSQAEALSXuXMJgG3uu5L+T9Jdku4PlDwhAAIA4AoBsK0dKKlDthPouZIekfSxpCmyp4fzhAAIAIArBMC29p6k35f+ny5pDdmOoa+X/XWQPCEAAgDgCgGwrc2UtFrp/wmSNi79v76kMVlUqAYCIAAArnQTAHe55DkCYAsbpXLoe1fSQaX/vy9paiY1ikYABADAlbgBcDgBsBXdIemE0v+nShon6R+SRirZTSA/lL1+8CvZD3nfwPCbS8/7y2sJ60oABADAFQJgW1tW0oql/xeQdJKkh2V/EWSZBNP5saRzJf1U0QFwsKQVfGXZhHUlAAIA4AoBsK3dLulISes4nGZUAHywwekSAAEAcIUA2Nauk+32pUv29O2dko6WtF4D04wKgFNkTzEPkz3NvHzC6RIAAQBwhQAI2VOyB0q6VtJHkuar/ruAwwLgAZL2lLSRpL0kvSPpfUk9a0ynp+zC4pU+IgACAOAGARCSeknaTdL5kl6V7RT67TqnFRYAg74t2wH1T2uM01/VN44QAAEAcIEA2Nb+Jns37mxJb0q6WNLekpZuYJpxAqAkfSqpX43hHAEEACAtBMC21iVprKSTZTt/diFOAFxO0hxJv0kwXa4BBADAFQJgW9tU0h9l+/wbL+lrSXdL+p2SBcLFJW1WKkbS8aX/VykNu0i2c+nVJO0g6RVJoyUtkWAeBEAAAFwhAMJnU0k3SZoneyNIXDso5Ho92bt/F5X0hOwdwB2S/lt6fuWEdUs/AH76qftpAwCQRwTAtre57BG7hyRNktQp6d+SLsyyUiEIgAAAuEIAbGuTZY/2/Vv2NO1PlN8PiAAIAIAr/gA4eXLVYAJga8tz4AsiAAIA4EpHBwEQhUAABADAFQIgCoIACACAKwRAFAQBEAAAVwiAKAgCIAAArhAAURAEQAAAXCEAoiAIgAAAuEIAREGkFwCXWIIACABoLwRAFAQBEAAAVwiAKAgCIAAArhAAURAEQAAAXOkmAO58MQEQ+UAABADAlZgB8JXhExqaDQEQjSIAAgDgCgEQBUEABADAFQIgCoIACACAKwRAFAQBEAAAVwiAKAgCIAAArhAAURAEQAAAXCEAoiAIgAAAuOIPgJMmVQ0mACIv0g+Aw4e7nzYAAHk0bx4BEIVAAAQAwBUCIAqCAAgAgCsEQBQEARAAAFcIgCgIAiAAAK4QAFEQBEAAAFwhAKIgCIAAALhCAERBEAABAHCFAIiCIAACAOAKARAFQQAEAMAVAiAKggAIAIArBEAUBAEQAABXCIAoCAIgAACuEABREOkFwMUXJwACANoLARAFQQAEAMAVAiAc+KGkRyR9Jfsh7xsY3kNS/9Lw2ZKek7RhwnkQAAEAcKWbALgTARAx/FjSuZJ+qvAA2E/StNLwjSTdJRsGl0gwDwIgAACuxAyAr35GAEQ8wQDYQ9IY2RDo6SlpiqTfJpguARAAAFcIgHAsGADXKD23eWC8hyTdkmC6BEAAAFzxB8CJE6sGEwCRVDAAblN6bsXAeNdLeqLGdHrKLixe6SMCIAAAbnR2EgDhVFQAsyX9CgAAEV1JREFU/HZgvH9IerzGdPqXXldRCIAAADhAAIRjrk4BcwQQAIC0EADhWNRNICf5nltY3AQCAEB2CIBwYHFJm5WKkXR86f9VSsP7yQa+/WS7gblDeewG5rPP3E8bAIA8IgDCgR0Ucr2epJtLw72OoMdImiPpedkgmAQBEAAAVwiAKAgCIAAArhAAURAEQAAAXCEAoiAIgAAAuEIAREEQAAEAcIUAiIIgAAIA4AoBEAVBAAQAwBUCIAqCAAgAgCsEQBQEARAAAFcIgCgIAiAAAK4QAFEQBEAAAFwhAKIgCIAAALhCAERBEAABAHCFAIiCSC8A9upFAAQAtBd/AJxQHfIIgMgLAiAAAK7Mn08ARCEQAAEAcKWbAPiji4YQAJELBEAAAFyJGQBfIwAiYwRAAABcIQCiIAiAAAC4QgBEQRAAAQBwhQCIgiAAAgDgCgEQBZF+APz8c/fTBgAgjwiAKAgCIAAArhAAURAEQAAAXCEAoiAIgAAAuEIAREEQAAEAcIUAiIIgAAIA4AoBEAVBAAQAwBUCIAqCAAgAgCsEQBQEARAAAFcIgCgIAiAAAK4QAFEQBEAAAFzxB8Dx46sGEwCRFwRAAABc6eoiAKIQCIAAALhCAERBEAABAHCFAIiCIAACAOAKARAFQQAEAMAVAiCaoL/sh+8vXyecBgEQAABXCIBogv6S3pe0gq/0TjiN9ALgYosRAAEA7YUAiCboL+mdBqdBAAQAwJVuAuCOBEA40F/STElfSRoh6S5Ja3Tzmp6yC4tX+ogACACAGzED4OufT2xoNgTA9vZjST+TtLGknSU9J3sN4HI1XtNf1dcNEgABAHCBAIgM9JINgCfUGIcjgAAApIUAiIw8JemaBONzDSAAAK4QAJGBnpJGSzojwWvSD4AjRrifNgAAeUQARBNcJGl7SatL2krSI5KmSVo1wTQIgAAAuEIARBPcJXsHcIekLyXdJ2mDhNMgAAIA4AoBEAVBAAQAwBUCIAqCAAgAgCsEQBQEARAAAFcIgCgIAiAAAK74A+C4cVWDCYDICwIgAAAuEQBRAARAAABcIgCiAAiAAAC4RABEARAAAQBwiQCIAiAAAgDgEgEQBUAABADAJQIgCoAACACASwRAFAABEAAAlwiAKAACIAAALhEAUQAEQAAAXCIAogAIgAAAuEQARAEQAAEAcIkAiAJILwAuuigBEADQfmoFwAsJgMgHAiAAAC7FCIBvjCAAIlsEQAAAXCIAogDSD4AjR7qfNgAAeUUARAEQAAEAcIkAiAIgAAIA4BIBEAVAAAQAwCUvAI4dWzWIAIi8IAACAOASARAFQAAEAMAlAiAKgAAIAIBLBEAUAAEQAACXCIAoAAIgAAAuEQBRAARAAABcIgCiAAiAAAC4RABEARAAAQBwiQCIAiAAAgDgEgEQBUAABADAJQIgCoAACACASwRAFAABEAAAlwiAKAACIAAALhEA0STHSBohaY6ktyRtl+C1BEAAAFzyAuBGG1UNIgDClQMkdUg6QtL6ki6TNEPSKjFfTwAEAMAlLwBKVYMIgHDldUnXBJ77SNL5MV9PAAQAwKUaAXAHAiAcWFhSp6T9As9fLun5iNf0lF1YvNJHaQTAjg5jFl6YAAgAaD/+ACgZ89VXxhhj7nj9v2bVfoPMqv0GmSfeH9PQLAiA7W1F2Q9/m8Dzf5X0ScRr+pdeU1GcB8C77iov+CNGuJ02AAB5FgyA559vjDFm0Ltf/S8A3vTS5w3NggDY3rwA+P3A86dK+jjiNc05Ati3b3nB/+ILt9MGACDPggHwmWeMMcZ0dXWZsx/5wPS7992GZ0EAbG/1nAIOSucawI4OY7bc0piVVnI7XQAAQACEXpd0deC5D5WHm0AAAEAqCIDwuoE5TLYbmEtlu4FZNebrCYAAABQMARCS7Qh6pKS5sh1B/zDBawmAAAAUDAEQjSIAAgBQMARANIoACABAwRAA0SgCIAAABUMARKMIgAAAFAwBEI0iAAIAUDAEQDSKAAgAQMEQANEoAiAAAAVDAESjCIAAABQMARCNIgACAFAwBEA0igAIAEDBEADRKAIgAAAFQwBEo5aUZEaNGmWmTp1KoVAoFAqlAGXUqFEEQDSkj+wCRKFQKBQKpXilj4A69JBdeJZ0VLxA6XKaRS+0Ce1E+9BOeSi0Teu1Ux/Z/TiQuSVlV5wls65IjtAm8dBOtdE+8dBO0WibeGgnoA6sONVok3hop9pon3hop2i0TTy0E1AHVpxqtEk8tFNttE88tFM02iYe2gmoQ09J/Ut/YdEm8dBOtdE+8dBO0WibeGgnAAAAAAAAAAAAAAAAAAAAAAAAAMVEj+MAkF+9xXYagGNLSVpC5Y0LGxm7sV1d5f6iaJNwK0n6paRlsq5IDq0h6Z+Sts26Ijm3mqRbJP0843rk1WqSHpN0X+nxAtlVJddWlLSTpHWzrghQBD0kXSbpHUkvSxooGwS9Ye2oh6QrJI2U9B9JwyStlWWFcmxlSeMldUnaJeO65EkPSdfKtsvNknplWpt8O0vSLEn3SNpH0kLZVidXeki6TtI8SaNl1zX6rwt3uaTJkt6UNEfS4ZIWybRGQI5tLeltSa9I2kPSGZI+knRXaXg7fsv8kaT3ZdtkB0m7S3pB0hMZ1inPviXpQdmw/Lik5TOtTT7sLWmi7JeHLQLD2vVLVZTNJL0k22ao9GdJUyW9LttOP5H0gaqXqXa3kKTrZbfZ35X0TUlXy7YVYRkIsaCkv0m6U+UjfpL0C0mfS1ohi0rlwEmSzpO0uO+5v0h6WLbNJHbifrtJekDSKrJHu/qqvNFt13a6XtIIlXfUm0k6WNL3ZC+1QNnfZU9tStI2smcjTpRtq3bWS/YL1aG+5zaRPRK4ZelxO35BD7OypHclHeF77vuyRwIXy6RGQA4tGPh/d0k/DIxziOw3p15qjx34goHHy8qGGU9v2W/gl6i9j1L426mHysvG1pL+Vfr/Nkkfyp528QfoVhdchtaW9Lykm2TD8QhJ/5Y9PfWKKpevduJvpwVKj/8p6VRJv5c0RtLtsmclvpINgu0kuBz5t78LyG6bPpLUr2k1yqdgO60jabrKAXBJ2cuZHpJdtlZrWs2AnDpTduN6nqQ+IcO9626Olz0l0w7hL9gmwff8C0nzZdvjXtnTerdIWrWJdcyDWu30R0mP+h7Pkm2v0ZJ+0KwKZijYNt7O6a+yIeYe2SOAa0haX9IXskG53W6YiWqnOyW9IekO2S9Y3rJ1meylF+1yXWnUOuYPO8tKekvSOc2tWq4E28nbb90g6RPZS3W6Sn9Plr1++1lJeza9pkAOrCx79OE9SRfKXqv1lqTDSsMXCPx9QNIFpf9bNQR21yaenVR5hPQ7kmZL2jn9KuZCnHY6XdKxpf/3kjRNUqda/+hNVNscWRq+uGwbrBl43Y8lzQ15vlV1105bSeqQNEOVd26uLXsksG+T6pmV7tYxbxvsbZ8fUvkLV6tun8NEtdNRpeGLyJ6tub803LOipCGywZFT5mg7v5G9w9c74rCI7F2J76p8LYn3LXNRSf+VDT6e9dOvYtPFaZMwi8ge4TqixjitJKqd3pHdcUv2SM2zsqc8J0n6k+wpz3+ota93q7UMeW2zZPXLtKls4Nk15frlRZx1baDs6fHgEePRav3TnUm2z5L9wvW27I0O7aRWO3239NwqsqFwm9JjLyB/IHtTCNB2zpRdKfxdK2whabAqT91Jdqc0QjYIbiD7zWmWWu+GkCRt4neo7OnNFdOrWq7EaaczZLum+IfKR7V+InsappVDTr3L0GmyYbld7k6MaqfHJQ0qPV5b0gTZU3ublZ7bXtJrkjZuTjUzE3c58sLMcZI+lr37vp2OANZqJ+8mol6y+6tf+8bZVvZmkL2aUEcgd86X9Kqqr/s7XHZD4l8x/lwa90LZ03i3q/IO4VaRpE3Wlu0D8O+Sxko6Qe2z4a3VTsNku8n5puyXhWCbHKvWDjlJlqF1VV6GviyNI7XHclSrnT5R+caqfWSP8EyQveN+pqSLVX3Bf6uJuxx5wWc92S9XmzSldvnRXTt5y9E5sndKPyR7hHC6bP+A9C2JtuJd7+BtMPYLDF9X0jOyK5bnmdK4L8he79ZqkrbJMpLOle0W5yXZ03ftIE47PavKZcfT6jvsepahC2VvCHlZ7bPjjttOF6gchFeTvfHqZEkbpV/FTNWzfZbs5QXXyZ6VaYcvEPW00/GSLpW9MaTVjyCjjdXq5dz/jece2U5pg9eNPCV76k6yK9qBKn5XJy7bRLJHt7ZR63HZTq22I3K9DG2i6u6WWgHLUG2ut8+tyvX6BrS0hWUvvr9P0q2q3Ll8IzDeWrIXx86SvXV+6dKwHrK3yv9f2pVtEtokHtopGm0TD+1UG+0TD+0EJLSvbF9iQ2Svs3pX9vTkzwLj/VH2OpqTSo+PlDRc0tOyh9EvkfS1WuPoFm0SD+0UjbaJh3aqjfaJh3YCElpT9ttSf99zy8tej+V1l9BT0jWyNyz8SpWnDn4ie2fZy7Kdr26dbnWbgjaJh3aKRtvEQzvVRvvEQzsBCXjXx6wv2+3G6qXH3vURT6ryWpG1VdkHWfD6keVTqGOz0Sbx0E7RaJt4aKfaaJ94aCcgAe8H0WvdWfkN2UPnwV+xaFW0STy0UzTaJh7aqTbaJx7aCUhgX9n+wiaq/EPW/pXHf8dcL9n+tLZSa6NN4qGdotE28dBOtdE+8dBOQEIHy17XcKekFyVd2834u8n2Meb/cflvlf62StcBtEk8tFM02iYe2qk22ice2glIwPtmtJVsp5arSPqLbC/nOwTG8btW0h2l/zeXvavqfrXGSkObxEM7RaNt4qGdaqN94qGdgATWV/XP1HiPN5T9SZuw34SU7MrxoKQTJV0pab6kW1TZh1IR0Sbx0E7RaJt4aKfaaJ94aCcggZ9LGil73cO7sn0heQu8f+U4VNIHpb9S5TeilWR/HqdL9pb49dOrblPQJvHQTtFom3hop9pon3hoJyChnWR/X/b3kraX7Q9pjuyh8sVL43jfnvrI/qbhG75hC5f+bijpLkk7p17j9NEm8dBO0WibeGin2mifeGgnIAHvG9Hpkt5WeQWQ7PUSb0g6JOR1e0p6U3YF20TSIEkrp1bL5qJN4qGdotE28dBOtdE+8dBOQAPukvSv0v/e4fKlZFeIe2UvnJXKF8suJukq2UPk8yQ9LttTeiuhTeKhnaLRNvHQTrXRPvHQTkANu0i6QtKfVO4MU7K/Zzhd5cPj3srzS0nvyfaf5Okl6ThJnbJ3Rm2cYn2bgTaJh3aKRtvEQzvVRvvEQzsBCXxb0iOyv2N4m+zKMEXllWcdSaMlnV167D+EPlz20LpnA0mvSfp1ivVtBtokHtopGm0TD+1UG+0TD+0EJLSYpJtlD4+v7nv+DUk3lf5fQtKpkmapfA2Ed1fUYJX7Q2oVtEk8tFM02iYe2qk22ice2gmo03WSdi/97x0eP0P2G5B3Ae3qsr99+JrKK9hKkoZK2q851Wwq2iQe2ikabRMP7VQb7RMP7QTUwd+BpbeiDJR0fWC8PpKGSfpC0n2yh9Ofk7RCyvXLAm0SD+0UjbaJh3aqjfaJh3YCHHlBUt/S/wuofKh8LUkHSLpI0m+aX61M0Sbx0E7RaJt4aKfaaJ94aCcgoTUkfS1pC99zC0eM2y5ok3hop2i0TTy0U220Tzy0E5CAd9j8N7J3RXnOlHSNpOWbXqPs0Sbx0E7RaJt4aKfaaJ94aCegAX+X9DfZ/pRGyN5Wv2umNcoebRIP7RSNtomHdqqN9omHdgISWkTSp7K9n8+R1C/b6uQCbRIP7RSNtomHdqqN9omHdgLq9JSkq2VXIli0STy0UzTaJh7aqTbaJx7aCajDgt2P0nZok3hop2i0TTy0U220Tzy0EwAAAAAAAAAAAAAAAAAAALr1/+G5FCGPiWlrAAAAAElFTkSuQmCC\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure()\n",
    "data[1].plot(color='r')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}