temperaturePID.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import pyvisa\n",
    "import time\n",
    "from matplotlib.animation import FuncAnimation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Lens driver\r\n"
     ]
    }
   ],
   "source": [
    "rm = pyvisa.ResourceManager()\n",
    "rm.list_resources()\n",
    "inst = rm.open_resource(\"ASRL11::INSTR\")\n",
    "inst.read_termination = '\\n'\n",
    "inst.write_termination = '\\n'\n",
    "print(inst.query(\"*IDN?\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "24.0"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "float(inst.query(\":TEMP:PID:SET?\")[:-1])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "\n",
       "mpl.get_websocket_type = function() {\n",
       "    if (typeof(WebSocket) !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert('Your browser does not have WebSocket support. ' +\n",
       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "              'Firefox 4 and 5 are also supported but you ' +\n",
       "              'have to enable WebSockets in about:config.');\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent = (\n",
       "                \"This browser does not support binary websocket messages. \" +\n",
       "                    \"Performance may be slow.\");\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen =  function () {\n",
       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
       "            fig.send_message(\"send_image_mode\", {});\n",
       "            if (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
       "            }\n",
       "            fig.send_message(\"refresh\", {});\n",
       "        }\n",
       "\n",
       "    this.imageObj.onload = function() {\n",
       "            if (fig.image_mode == 'full') {\n",
       "                // Full images could contain transparency (where diff images\n",
       "                // almost always do), so we need to clear the canvas so that\n",
       "                // there is no ghosting.\n",
       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "            }\n",
       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "        };\n",
       "\n",
       "    this.imageObj.onunload = function() {\n",
       "        fig.ws.close();\n",
       "    }\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_header = function() {\n",
       "    var titlebar = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus () {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>');\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option);\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_message = function(type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function() {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
       "    }\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0, 0, fig.canvas.width / mpl.ratio, fig.canvas.height / mpl.ratio);\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch(cursor)\n",
       "    {\n",
       "    case 0:\n",
       "        cursor = 'pointer';\n",
       "        break;\n",
       "    case 1:\n",
       "        cursor = 'default';\n",
       "        break;\n",
       "    case 2:\n",
       "        cursor = 'crosshair';\n",
       "        break;\n",
       "    case 3:\n",
       "        cursor = 'move';\n",
       "        break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message(\"ack\", {});\n",
       "}\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = \"image/png\";\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src);\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig[\"handle_\" + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "}\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function(e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e)\n",
       "        e = window.event;\n",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\n",
       "\n",
       "    return {\"x\": x, \"y\": y};\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys (original) {\n",
       "  return Object.keys(original).reduce(function (obj, key) {\n",
       "    if (typeof original[key] !== 'object')\n",
       "        obj[key] = original[key]\n",
       "    return obj;\n",
       "  }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
       "    var canvas_pos = mpl.findpos(event)\n",
       "\n",
       "    if (name === 'button_press')\n",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.key_event = function(event, name) {\n",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, {key: value,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
       "    if (name == 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message(\"toolbar_button\", {name: name});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.close = function() {\n",
       "        comm.close()\n",
       "    };\n",
       "    ws.send = function(m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function(msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['data'])\n",
       "    });\n",
       "    return ws;\n",
       "}\n",
       "\n",
       "mpl.mpl_figure_comm = function(comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = $(\"#\" + id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm)\n",
       "\n",
       "    function ondownload(figure, format) {\n",
       "        window.open(figure.imageObj.src);\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy,\n",
       "                           ondownload,\n",
       "                           element.get(0));\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element.get(0);\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error(\"Failed to find cell for figure\", id, fig);\n",
       "        return;\n",
       "    }\n",
       "\n",
       "    var output_index = fig.cell_info[2]\n",
       "    var cell = fig.cell_info[0];\n",
       "\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
       "    var width = fig.canvas.width/mpl.ratio\n",
       "    fig.root.unbind('remove')\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable()\n",
       "    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n",
       "    fig.close_ws(fig, msg);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.close_ws = function(fig, msg){\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width/mpl.ratio\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message(\"ack\", {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () { fig.push_to_output() }, 1000);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>');\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items){\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) { continue; };\n",
       "\n",
       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
       "    buttongrp.append(button);\n",
       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
       "    titlebar.prepend(buttongrp);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(el){\n",
       "    var fig = this\n",
       "    el.on(\"remove\", function(){\n",
       "\tfig.close_ws(fig, {});\n",
       "    });\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
       "    // this is important to make the div 'focusable\n",
       "    el.attr('tabindex', 0)\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    }\n",
       "    else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager)\n",
       "        manager = IPython.keyboard_manager;\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which == 13) {\n",
       "        this.canvas_div.blur();\n",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.find_output_cell = function(html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i=0; i<ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code'){\n",
       "            for (var j=0; j<cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] == html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "}\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel != null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAyAAAAJYCAYAAACadoJwAAAgAElEQVR4XuzdC7hNdfrA8beQhGnSRcqdGapx718uNSRSEoUZokYuJbdRSUQGI2kijXH3EBplNGhCkkIajJFTucwouTPIbSiRa//nt85zzpxznHPsvd937bPXOt/1PJ75/7P2u9b7+b3Oed+91tr7kh9//PFHYUMAAQQQQAABBBBAAAEE4iBwCQNIHJQ5BAIIIIAAAggggAACCHgCDCAUAgIIIIAAAggggAACCMRNgAEkbtQcCAEEEEAAAQQQQAABBBhAqAEEEEAAAQQQQAABBBCImwADSNyoORACCCCAAAIIIIAAAggwgFADCCCAAAIIIIAAAgggEDcBBpC4UXMgBBBAAAEEEEAAAQQQYAChBhBAAAEEEEAAAQQQQCBuAgwgcaPmQAgggAACCCCAAAIIIMAAQg0ggAACCCCAAAIIIIBA3AQYQOJGzYEQQAABBBBAAAEEEECAAYQaQAABBBBAAAEEEEAAgbgJMIDEjZoDIYAAAggggAACCCCAAAMINYAAAggggAACCCCAAAJxE2AAiRs1B0IAAQQQQAABBBBAAAEGEGoAAQQQQAABBBBAAAEE4ibAABI3ag6EAAIIIIAAAggggAACDCDUAAIIIIAAAggggAACCMRNgAEkbtQcCAEEEEAAAQQQQAABBBhAqAEEEEAAAQQQQAABBBCImwADSNyoORACCCCAAAIIIIAAAggwgFADCCCAAAIIIIAAAgggEDcBBpC4UXMgBBBAAAEEEEAAAQQQYAChBhBAAAEEEEAAAQQQQCBuAgwgcaPmQAgggAACCCCAAAIIIMAAQg0ggAACCCCAAAIIIIBA3AQYQOJGzYEQQAABBBBAAAEEEECAAYQaQAABBBBAAAEEEEAAgbgJMIDEjZoDIYAAAggggAACCCCAAAMINYAAAggggAACCCCAAAJxE2AAiRs1B0IAAQQQQAABBBBAAAEGEGoAAQQQQAABBBBAAAEE4ibAABI3ag6EAAIIIIAAAggggAACDCDUAAIIIIAAAggggAACCMRNgAEkbtQcCAEEEEAAAQQQQAABBBhAqAEEEEAAAQQQQAABBBCImwADSNyoORACCCCAAAIIIIAAAggwgFADCCCAAAIIIIAAAgggEDcBBpC4UXMgBBBAAAEEEEAAAQQQYAChBhBAAAEEEEAAAQQQQCBuAgwgcaPmQAgggAACCCCAAAIIIMAAQg0ggAACCCCAAAIIIIBA3AQYQOJGzYEQQAABBBBAAAEEEECAAYQaQAABBBBAAAEEEEAAgbgJMIDEjZoDIYAAAggggAACCCCAAAMINYAAAggggAACCCCAAAJxE2AAiRs1B0IAAQQQQAABBBBAAAEGEGoAAQQQQAABBBBAAAEE4ibAABI3ag6EAAIIIIAAAggggAACDCDUAAIIIIAAAggggAACCMRNgAEkbtQcCAEEEEAAAQQQQAABBBhAqAEEEEAAAQQQQAABBBCImwADSNyoORACCCCAAAIIIIAAAggwgFADCCCAAAIIIIAAAgggEDcBBpC4UXMgBBBAAAEEEEAAAQQQYAChBhBAAAEEEEAAAQQQQCBuAgwgcaPmQAgggAACCCCAAAIIIMAAQg0ggAACCCCAAAIIIIBA3AQYQOJGzYEQQAABBBBAAAEEEECAAYQaQAABBBBAAAEEEEAAgbgJMIDEjZoDIYAAAggggAACCCCAAAMINYAAAggggAACCCCAAAJxE2AAiRs1B0IAAQQQQAABBBBAAAEGEGoAAQQQQAABBBBAAAEE4ibAABI3ag6EAAIIIIAAAggggAACDCAXqYGlS5fK3Xff7e319ddfS/ny5VNfcfLkSRk0aJC89dZbcvDgQe/vevbsKY8//jiVhQACCCCAAAIIIIAAApkIMIBkUxanT5+WKlWqyO7du+X777+/YAC5//77ZfHixdK9e3e5+eab5b333pN3331XXnnlFenduzcFhwACCCCAAAIIIIAAAhkEGECyKYlhw4bJH//4R2nTpo33v2mvgCxYsEAeeOABGTlypDz99NOpUZo1ayYffvih7Ny5U6699loKDgEEEEAAAQQQQAABBNIIMIBkUQ5ugHBXNcaMGeMNE4MHD043gLRt21bmzp0rR44ckQIFCqRGWbZsmdSvX18mTZrErVj8U0MAAQQQQAABBBBAgCsgkdWAu5Jx4MABWbVqlTd8ZBxAKlSoIFdddZWsXr06XUD3XMgVV1whnTt3lgkTJkR2MPZCAAEEEEAAAQQQQCCXCHAFJJOFdrdXuQFkzZo1UqNGDe9B84wDSOHCheWee+6ROXPmXBChSJEiUqdOHZk/f36WZbRv3z5xf9Ju27dvl3nz5kmtWrXExWdDAAEEEEAAAQQQSCyB7777TtavX+998JB7Q5otegEGkAxm7grGLbfcIg0bNpSJEyd6f5vZAJInTx5p1aqV9wlYGbcbbrjBu33ro48+ynJFUmJGv2S8AgEEEEAAAQQQQCCnBbp06SLjxo3L6dMI5PEZQDIs2wsvvOAV0+bNm+Waa67JcgDx4wrI7NmzxT347graXUFhQwABBBBAAAEEEEgsgZUrV8r48eO9ftH1bGzRCzCApDHbu3evlC1bVp566inp1KlT6t/86U9/ktGjR8uSJUukdOnS3j5+PAPy5ptvyiOPPCIzZswQ95A7GwIIIIAAAggggEBiCdCv6deDASSN4RdffCHVqlXLVrVgwYJy/Phx76N533nnHdNPwaKg9QVNBAQQQAABBBBAwE8B+jW9LgNIGsNjx4553+GRcXv77bflr3/9q3cVpHjx4vLggw96D5g3bdo0y+8B2bFjh1x33XVRrRAFHRUXOyOAAAIIIIAAAnEXoF/TkzOARGCY2UPo7mX33nuvd1tWjx49vIfO3adnuW9Cd89x9O3bN4LI6XehoKMm4wUIIIAAAggggEBcBejX9NwMIBEYZjWAnDhxQgYOHCgzZ86UgwcPSrly5byPZHPfARLLRkHHosZrEEAAAQQQQACB+AnQr+mtGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CMTar61du1ZmzJghS5cule3bt0vBggXllltukeeff14aNGjgy7kmalAGkARamVgLOoFS4FQQQAABBBBAAIFQC8Tar7Vs2VKWL18uLVq0kOrVq8vx48dl6tSpsnHjRhk3bpx06dIl1G5pk2MASaCljrWgEygFTgUBBBBAAAEEEAi1QKz92sqVK+XWW2+V/Pnzp/qcPHlSqlatKgcPHpQDBw5I3rx5Q22XkhwDSAItc6wFnUApcCoIIIAAAggggECoBaz7tV69esnIkSNl165dUqJEiVDbMYAk4PJaF3QCpsgpIYAAAggggAACgRaw7tcefvhhmT17thw9etR7LiQ3bFwBSaBVti7oBEqNU0EAAQQQQAABBEIhkNKvDRkyRBo3bpwup2LFion7E+m2adMm7xasJk2ayJw5cyJ9WeD3YwBJoCVkAEmgxeBUEEAAAQQQQACBTARS+rXMcAYOHCiDBg2KyO3YsWNSq1Yt2bdvn6xbt05KliwZ0evCsBMDSAKtIgNIAi0Gp4IAAggggAACCGQzgGiugLiHzxs1aiSffvqpLFq0SOrWrZurrBlAEmi5GUASaDE4FQQQQAABBBBAIJsBxH2nR9u2baM2On36tDRt2tT7PpC5c+d6t1/lto0BJIFWnAEkgRaDU0EAAQQQQAABBIwHkLNnz4r7PpD58+eL6/tat26dK40ZQBJo2RlAEmgxOBUEEEAAAQQQQMBwADl//ry0adNGZs2aJZMmTZLHH3881/oygCTQ0jOAJNBicCoIIIAAAggggIDhAPLMM8/Ia6+95j3v0alTpwsiN2zYUIoWLZorzBlAEmiZGUASaDE4FQQQQAABBBBAwHAAqVevnixfvjxL02XLlonbJzdsDCAJtMoMIAm0GJwKAggggAACCCBgOICA+T8BBpAEqgYGkARaDE4FAQQQQAABBBBgAPGlBhhAfGGNLSgDSGxuvAoBBBBAAAEEEIiXAP2aXpoBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBJI3hpk2bZPDgwZKUlCT79u2TSy+9VMqWLSuPPfaYdOnSRfLnz+/tvWPHDilTpkym+h07dpTJkyfHtDIUdExsvAgBBBBAAAEEEIibAP2anpoBJI3h4sWLZcSIEXL77bdL8eLF5dy5c7Jy5UqZOXOm3H///TJ//vx0A0izZs2kZcuW6VahfPnyUrNmzZhWhoKOiY0XIYAAAggggAACcROgX9NTM4BEYNi9e3cZO3asfPnll1KhQoXUKyD9+/eXF198MYIIke1CQUfmxF4IIIAAAggggEBOCdCv6eUZQCIwdFdFevfuLatXr/aujqTcguUGEPfHbQUKFIggUva7UNBqQgIggAACCCCAAAK+CtCv6XkZQDIxPHHihLg/33//vaxZs0a6desm+fLlky1btniDRsoAUqhQITl+/LgXwd169dRTT3n7xrpR0LHK8ToEEEAAAQQQQCA+AvRremcGkEwMBw0a5D2MnrK5qx6TJk2SypUre/9p165d0qFDB3nwwQelVKlSsnfvXu/B87Vr18qzzz4rw4cPv+jKuIfc3Z+028KFC2XAgAEyY8YMadu27UVjsAMCCCCAAAIIIIBAfAUYQPTeDCCZGG7btk3cn8OHD8vSpUtl/fr1MmzYMKlXr16W4u6B9fr168uKFStk8+bNUq5cuWxXJ+OQk3ZnBhB9YRMBAQQQQAABBBDwQ4ABRK/KABKB4WuvvSZ9+vSRdevWyU033ZTlK9ynZDVt2lQmTpwoTzzxRLaRuQISATy7IIAAAggggAACCSbAAKJfEAaQCAy/+eYbuf7666Vfv34ydOjQLF/hrpRUqVLF28ftG+1GQUcrxv4IIIAAAggggEB8BejX9N4MIBEYpjx07r6McNy4cVm+Yu7cudKiRQvveRD3hYTRbhR0tGLsjwACCCCAAAIIxFeAfk3vzQCSxvDAgQNy3XXXXaD63HPPeQ+WT5s2Tdq1aydHjhyRIkWKpNvv5MmTUqdOHdm4caNs3bpVSpQoEfXqUNBRk/ECBBBAAAEEEEAgrgL0a3puBpA0hg899JD34Ll72NwNEEePHpUPPvhAlixZInfccYcsW7ZM8ubNK82bN/c+ptd947n7xnT3KVjTp0/3Hlx3D6v37ds3ppWhoGNi40UIIIAAAggggEDcBOjX9NQMIGkMZ82a5V3lcM9yHDx4UPLnzy8VK1aUVq1aSY8ePbz/321TpkzxBo6vvvrKuxrivg+kevXq0rNnT+8h9Fg3CjpWOV6HAAIIIIAAAgjER4B+Te/MAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EWAfk3PygCiNzSLQEGbURIIAQQQQAABBBDwRYB+Tc/KAKI3NItAQZtREggBBBBAAAEEEPBFgH5Nz8oAojc0i0BBm1ESCAEEEEAAAQQQ8EVA068dP35cRowYIUlJSbJ27VrZv3+/tGvXTqZNm+bLuSZqUAaQBFoZTUEnUBqcCgIIIIAAAgggEFoBTb+2Y8cOKVOmjBQrVkxq1KghCxYsYAAJbaUEJDFNQQckRU4TAQQQQAABBBAItICmXzt16pQcOnRIbrzxRjl79qzky5ePASTQ1RCCk9cUdAjSJwUEEEAAAQQQQCDhBaz6NQaQhF/q3HGCVgWdO7TIEgEEEEAAAQQQiL+AVb/GABL/teOImQhYFTS4CCCAAAIIIIAAAv4IpPRrQ4YMkcaNG6c7iHu2w/2JZGMAiUSJfXwXYADxnZgDIIAAAggggAACKoGUfi2zIAMHDpRBgwZFFJ8BJCImdvJbgAHEb2HiI4AAAggggAACOgGugOj83Kv5GF69oVkEBhAzSgIhgAACCCCAAAK+CFj1a1wB8WV5CBqtgFVBR3tc9kcAAQQQQAABBBCITMCqX2MAicybvXwWsCpon0+T8AgggAACCCCAQK4VsOrXGEBybQklVuJWBZ1YWXE2CCCAAAIIIIBAeAS0/dqYMWPk6NGjcv78eXEPrVerVk2aN2/uATVt2lQqV64cHqwsMuEZkARaYm1BJ1AqnAoCCCCAAAIIIBBKAW2/Vrp0adm5c2emNlOnTpXHHnsslG5pk2IASaAl1hZ0AqXCqSCAAAIIIIAAAqEUoF/TLysDiN7QLAIFbUZJIAQQQAABBBBAwBcB+jU9KwNIGsNNmzbJ4MGDJSkpSfbt2yeXXnqplC1b1rsU1qVLF8mfP3/q3ufOnZPhw4fL5MmTZffu3VKiRAnp1KmT9O7dW/LkyRPTylDQMbHxIgQQQAABBBBAIG4C9Gt6agaQNIaLFy+WESNGyO233y7FixcXN2SsXLlSZs6cKffff7/Mnz8/de+uXbvK+PHjpX379lK7dm1ZtWqVuPv23H8fO3ZsTCtDQcfExosQQAABBBBAAIG4CdCv6akZQCIw7N69uzdUfPnll1KhQgXZsGGDVKlSRXr06CGjRo1KjdCzZ08ZPXq0rFu3TipVqhRB5PS7UNBRk/ECBBBAAAEEEEAgrgL0a3puBpAIDN1VEXdr1erVq72rI/3795eXXnpJtm3bJmXKlEmNsH37du+WrX79+snQoUMjiMwAEjUSL0AAAQQQQAABBHJQgAFEj88AkonhiRMnxP35/vvvZc2aNdKtWzfJly+fbNmyRQoUKCCNGjXyrnLs37//glcXLVrU+zznRYsWRb06FHTUZLwAAQQQQAABBBCIqwD9mp6bASQTw0GDBnkPo6ds7qrHpEmTUr8Yxt1eddlll3kPq2fcqlevLmfOnPFu08pucw+5uz9pt4ULF8qAAQNkxowZ0rZtW/3qEgEBBBBAAAEEEEDAVIABRM/JAJKJobu1yv05fPiwLF26VNavXy/Dhg2TevXqeXuXK1dO3JUO9+B5xs09kH7gwAHvakl2W8YhJ+2+DCD6wiYCAggggAACCCDghwADiF6VASQCw9dee0369Onj3XZ10003eQ+YcwUkAjh2QQABBBBAAAEEQibAAKJfUAaQCAy/+eYbuf7661MfLucZkAjQ2AUBBBBAAAEEEAihAAOIflEZQCIw3LFjh/dpV+7LCMeNG+cNIu6WLD4FKwI8dkEAAQQQQAABBEIkwACiX0wGkDSG7tmN66677gLV5557zvvW82nTpkm7du28W7HcJ11l9T0gX3zxReoD69EsEQUdjRb7IoAAAggggAAC8RegX9ObM4CkMXzooYe8B8/dw+YlSpSQo0ePygcffCBLliyRO+64Q5YtWyZ58+b1XvHkk0/KxIkTvW9Cr1OnjveN6e6b0Dt37iwTJkyIaWUo6JjYeBECCCCAAAIIIBA3Afo1PTUDSBrDWbNmeVc53KdeHTx4UPLnzy8VK1aUVq1aeVc73P+fsp09e1ZeeeUVmTx5suzZs0eKFy8unTp1Ene1JGVIiXZ5KOhoxdgfAQQQQAABBBCIrwD9mt6bAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNn5vveIAACAASURBVKMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBRG9oFoGCNqMkEAIIIIAAAggg4IsA/ZqelQFEb2gWgYI2oyQQAggggAACCCDgiwD9mp6VAURvaBaBgjajJBACCCCAAAIIIOCLAP2anpUBJI3h2rVrZcaMGbJ06VLZvn27FCxYUG655RZ5/vnnpUGDBql77tixQ8qUKZOpfseOHWXy5MkxrQwFHRMbL0IAAQQQQAABBOImoOnXzp07J8OHD/d6xd27d0uJEiWkU6dO0rt3b8mTJ0/ccsjpAzGApFmBli1byvLly6VFixZSvXp1OX78uEydOlU2btwo48aNky5dunh7pwwgzZo1E/eatFv58uWlZs2aMa2rpqBjOiAvQgABBBBAAAEEEIhKQNOvde3aVcaPHy/t27eX2rVry6pVq7xe0/33sWPHRnUeQd6ZASTN6q1cuVJuvfVWyZ8/f+p/PXnypFStWlUOHjwoBw4ckLx586YOIP3795cXX3zRbP01BW12EgRCAAEEEEAAAQQQyFIg1n5tw4YNUqVKFenRo4eMGjUqNX7Pnj1l9OjRsm7dOqlUqVKukGcAiWCZe/XqJSNHjpRdu3Z5l8pSroC4AcT9cVuBAgUiiJT9LrEWtPrABEAAAQQQQAABBBCISCDWfs31jC+99JJs27Yt3a387rb/smXLSr9+/WTo0KERnUPQd2IAiWAFH374YZk9e7YcPXrUey4kZQApVKiQd5uW29ytV0899ZR069YtgoiZ7xJrQcd8QF6IAAIIIIAAAgggEJVArP1ao0aNvKsc+/fvv+B4RYsWlWrVqsmiRYuiOpeg7swAcpGV27Rpk3cLVpMmTWTOnDne3u5KSIcOHeTBBx+UUqVKyd69e72HidxD7M8++6z3cNHFtn379on7k3ZbuHChDBgwQEqXXiFFitS5WAj+HgEEEEAAAQQQQCDOAkeOrJQdO+6QIUOGSOPGjdMdvVixYuL+ZLa526suu+wySUpKuuCv3bPHZ86cEXebVm7YGECyWeVjx45JrVq1vEHBTawlS5bMcm/3qQb169eXFStWyObNm6VcuXLZ1s+gQYNk8ODBWezjCrN6bqg/ckQAAQQQQAABBAIm8JmI1Mj0nAcOHCiux8tsc72hu9LhHjzPuLkH0t2zxlu2bAmYRWynywCShZt7+NxdKvv000+9y2F169a9qPD8+fOladOmMnHiRHniiSey3Z8rIBflZAcEEEAAAQQQQCDhBLgCol8SBpBMDE+fPu0NEu77QObOnevdfhXJtn79eu/TDdwDRO5Bomi3WO8pjPY47I8AAggggAACCCAQm0Cs/RrPgPzPmwEkQ+2dPXvW+24PdzXDFVjr1q0jrk43rLjvEHHPg7gvJIx2i7Wgoz0O+yOAAAIIIIAAAgjEJhBrv+benB42bBifgiUiDCBpau/8+fPSpk0bmTVrlkyaNEkef/zxTCvzyJEjUqRIkXR/527ZqlOnjvelhVu3bvU+rjfaLdaCjvY47I8AAggggAACCCAQm0Cs/Zp7nth90lVW3wPyxRdfSOXKlWM7qYC9igEkzYI988wz8tprr3nPe3Tq1OmCpWzYsKH38FDz5s3lxIkT3jeeFy9e3PsUrOnTp3sTrZts+/btG1MZxFrQMR2MFyGAAAIIIIAAAghELaDp15588knvWWH3TejujWv3Jdjum9A7d+4sEyZMiPpcgvoCBpA0K1evXj1Zvnx5lmu5bNkycftMmTLFGzi++uorcVdD3PeBuI9Pc99k6Z4diXXTFHSsx+R1CCCAAAIIIIAAApELaPo1d6v/K6+84t2uv2fPHu+NbPem93PPPSd58+aN/CQCvicDSAItoKagEygNTgUBBBBAAAEEEAitAP2afmkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhFSCvrKK6+UAgUKmMUlEAIIBFPAfT78n//8Z6lYsWIwE+CsEUAAgRAKMIDoF5UBRG9oFiGloM0CEggBBAIv0KdPH3n55ZcDnwcJIIAAAmERYADRryQDiN7QLEJKQVetWlVuvPFGs7gEQgCBYAmcP39e3n//fe+ku3XrJmPGjAlWApwtAgggEGIBBhD94jKA6A3NIlDQZpQEQiDQAufOnZO8efN6OTzxxBMyceLEQOfDySOAAAJhEqBf068mA4je0CwCBW1GSSAEAi+QJ08ecVdC2rdvL6+//nrg8yEBBBBAICwC9Gv6lWQA0RuaRaCgzSgJhEDgBS6//HI5deqUtG3bVmbMmBH4fEgAAQQQCIsA/Zp+JRlA9IZmEShoM0oCIRB4gcKFC8vx48fl17/+tcyaNSvw+ZAAAgggEBYB+jX9SjKA6A3NIlDQZpQEQiDwAldffbUcOXJEHnroIZk7d27g8yEBBBBAICwC9Gv6lWQA0RuaRaCgzSgJhEDgBa6//nr55ptv5P7775cFCxYEPh8SQAABBMIiQL+mX0kGEL2hWQQK2oySQAgEXqBEiRKyZ88eueeee+SDDz4IfD4kgAACCIRFgH5Nv5IMIHpDswgUtBklgRAIvEC5cuVk27Ztctddd8nSpUsDnw8JIIAAAmERoF/TryQDiN7QLAIFbUZJIAQCL1CxYkX56quvpE6dOrJixYrA50MCCCCAQFgE6Nf0K8kAojc0i0BBm1ESCIHAC1SqVEk2btwot912m/zzn/8MfD4kgAACCIRFgH5Nv5IMIHpDswgUtBklgRAIvECNGjXks88+k2rVqnn/y4YAAgggkBgC9Gv6dWAA0RuaRaCgzSgJhEDgBWrWrOld+bjlllu8KyFsCPgicOqUyEcfiRw8KNKggUjx4r4cJlcH3bdP5MMPRa6+WqRhQ5HLLouO48cfRVavFnE/B2rVEvnFLyJ7/ZkzIsuWiXzxhYiLkXa75JLkOG7Noz2fyI7uz15794q8957IkSORxXd5VqkiUr++SL58yQ7OcuVKkXPnIouRdq/ChUXuv1/eXLFCHnnkEe9LYt2XxbJFL8AAEr2Zb69gAPGNlsAXE3A/lNeuFfnHP0RuuknkrrtE8ua92KtETp8WWbxYJClJ5Pz55P1/8hORBx4Q+fnPs359yut27kz+BVihQvp9XSz3C+Lzz0Vq1Ej+pXvppRc/n6z2cM2V+yhb98uoSRORa65Jv+eJE8l/f/SoyH33iZQoceH5rFkj4h4Gdw1bHLY7p06VFbt3y8+vvlq+6tYt8yO6NbrjDpG6dZN99u8Xeeed5P/NuKXs+8tfiuTJk/5v3Tq8/77IVVd5v1ylUKH0f+9cnI9zyszH/SJ3Taz7xZ5SBy6Gi3XzzeljuVpbv17kk09EypVLbshcY5B227EjOQ93XD+3lObk3ntFChRIXluXpzu/eG7OoWlTkZ/+NP1Rf/gh+d+X+3eQ4pp2D2fr6rlgwfSv++675CZt06bss9i+XeTdd0W+/fZ/+9WuLVKvXmT//uNpFMRjpfwc+/jj/w0Abghp2VKkaNHIMnL/BtwauX+jKZtbd/dvy9VsVtuePSJ/+9vFG3X3b/7BBy/8mRfZ2cVvL/cz5u9/T/6TcZiK5CxcnnffLfLpp+ktI3ltJvu8+fOfyyObNzOAxOjnXsYAosCzfmnqAPLqq9LW/QP75hvrQ1wYr3z55B+GRYqk/7sDB0Tcl59t2+b/OVgewTW/7gdz1arJzeaXXyY3Mv/9r+VR4hurdOnkXxA33JD8jo37ZeYaYffulsXmmsqFC0VcM5KyXXttckOUsS7SHu/QoeRfjFm9E1WtWuaNjBsG3OvSrsltt4nceWdyE33ypMj8+el/Sbh3ZV2j5d59inZz7xq6Ji7l3S7XfLtfRJUrJ9eIe3fS/aI+fjw5svtvrkn/v/9L/r9dc+aaOfcLPY7b3SLiPvuqjIhc9F+h83FDnKuNi72rV6xY8tq6fyvu54wbOt2wl7K5psb9GyrjjiwiX3+dXB9uaEzZatYUcY2qs/z++2Q/985kZlulSsnebshwMdxHCrt/lymbqzF3Pq7m3Pm4BmH58jhKS3ID79bcDVA59bPCvQvtBqGUYXz37uRhKKUusxK54orkfxulSiXvsXmzyKJFcRuU47tQHA2BxBB4U0QeEWEAUSwHA4gCz/qlqQOIiMT1gp77xecaDtfgum3LluR3My/WyFgDWMb72c+S38V17xyGYXONsGv6du0S+c9/wpAROVxE4F4Rcd/+caOIxHf0YWlylYAb/n/9axH3Rod702nevIsPPbkKSJls/vzJA+KvfpX8Js/06ekH8EjCu+Gydevkod8N8H/9a/Itcxfb3Btxbm0bN77wKpl7o8dd9fzLX4LzezJl2HY5pbyBdDED9waJG8hnz06+yl+yZLJl8+bJt8RFu7mrim+9JW/OmSOPnDnDABKtX5r9GUAUeNYvzbEBxDoR4gVPwP1Qdj/UGzVKfhfY/VL6178unoe7t/bhh5Ov0KRcnfj3v0VmzRKZMyfrd5Pd69q0Sb4H2f0ydb8c0r7T697Rd/fVutuz3NWet98W2br14ueT2R7u1iP3zrI7nnuH3Z2b+8Wb9gqSa8IefVTkxhuTc3dNmLv9JWVzA627UvjQQ8n7xGF7oF07WfDRR3LdNdfIN+vWZX5EdxXK+b35ZvJg6tbP5eEaFTe0pt3cvm5NZs5MvqqRsrlf6s2aibRqlXw1yPm4KxAptzm4q1Lu/mm3zu5qhWtSM96245rXdu2SayjlNiJ3DBfLnZ87dsp23XXJ+7mrHu7heuft7lFP2VzD5s7HxXMNlJ+ba05cLbj6c1eB3O2H7rhunZ1LPDZXh67GncOSJenf+KlTJ7nuXAOZ8bY4d7uYa0bdeqa9LcVdlXLr5V7nbm9znllt7mqXuzUl7eauUqVdr3gYhPkYzjftrVLu35UbHs6ejSxr9+/P3a6V9t+ze3PQ3aWQ3a1Ibt0jbbAPHw7GFTP38+fyyyNzy2wv92/GveGa8WdjDBHfnDJFHunUiQEkBruUlzCAKPCsX5o6gNx1l7Tt0+d/t4BYHyglnvvF5375undk0t6j6m6VSGnY3K0TkTwL4Nc5RhvXvTvhGh7XJLl3eFwj4xpP1/Aa/NCJ9nTU+7tfUq65cA2SazZcc+caQZeTu5XGYnMuV155oY+79Si7q2CuLrK7Jcq9Nu295SnnmtnrXC2mvQUq4/m4X7QuVmb3wV/MwP3yz/hLy/0icreeuc3Ve8bmzjVhrjl1m2sA3O1Kca6fFi1ayNy5c+Wqq66SI5E+cHkxC/f3GS1d7hmfwXA2Kc+6OLuM95q7unTPGbgtq/pJORe3ZseO/e/M3NpmfJ7HxUppyFzjn13THEmOQd3H/cxKGXydQaRDUNr1yqzeg+rBeSOQoAI8s6tfGAYQvaFZhBwtaNcApjR3rkHM+JCqWZYEQgCBSARat24ts2bNkoIFC8rxiz0HEElA9kEAAQQQMBHI0X7NJIOcD8IAkvNrkHoGFHQCLQangkAOCzz66KPe5f3LLrtMTsXpk7dyOGUOjwACCARCgH5Nv0wMIHpDswgUtBklgRAIvEDHjh3l9ddfl0suuUTOnTvn/S8bAggggEDOC9Cv6deAAURvaBaBgjajJBACgRd48sknZeLEiV4eZ86ckbxBehYr8PokgAACCGQtQL+mrw4GEL2hWQQK2oySQAgEXqBHjx4yZswYL48TJ05Igey+dCzw2ZIAAgggEBwB+jX9WjGA6A3NIlDQZpQEQiDwAr169ZKRI0d6eRw7dkx+4j6Jiw0BBBBAIMcF6Nf0S8AAojc0i0BBm1ESCIHAC/Tt21f+8Ic/eHkcPHhQrrnmmsDnRAIIIIBAGATo1/SryACiNzSLQEGbURIIgcALDBgwQF588UUvj71790oxq+99CbwMCSCAAAI5K0C/pvdnANEbmkWgoM0oCYRA4AV+//vfy8CBA708du7cKSXdt9WzIYAAAgjkuAD9mn4JGED0hmYRKGgzSgIhEHiBl19+WZ5//nkvj6+//lrKly8f+JxIAAEEEAiDAP2afhUZQPSGZhEoaDNKAiEQeIFXX31Vnn32WS+Pf//733LTTTcFPicSQAABBMIgQL+mX0UGEL2hWQQK2oySQAgEXuBPf/qT9OzZ08tj3bp1Urly5cDnRAIIIIBAGATo1/SryACiNzSLQEGbURIIgcALTJgwQbp06eLl8emnn8qtt94a+JxIAAEEEAiDAP2afhUZQPSGZhEoaDNKAiEQeIEpU6ZIp06dvDxWrVoltWrVCnxOJIAAAgiEQYB+Tb+KDCB6Q7MIFLQZJYEQCLzAG2+8Ie3atfPyWL58ufzyl78MfE4kgAACCIRBgH5Nv4oMIHpDswgUtBklgRAIvMBf/vIXefjhh708PvzwQ2nQoEHgcyIBBBBAIAwC9Gv6VWQA0RuaRaCgzSgJhEDgBebMmSMtW7b08li4cKHcd999gc+JBBBAAIEwCNCv6VeRAURvaBaBgjajJBACgReYN2+eNGvWzMvj3XfflaZNmwY+JxJAAAEEwiBAv6ZfRQYQvaFZBArajJJACAReYNGiRalXPWbPni0tWrQIfE4kgAACCIRBgH5Nv4oMIHpDswgUtBklgRAIvMCSJUtSn/t46623Up8HCXxiJIAAAggEXIB+Tb+ADCB6Q7MIFLQZJYEQCLzAJ598InXr1vXymD59uvzmN78JfE4kgAACCIRBgH5Nv4oMIHpDswgUtBklgRAIvMDq1atTv/tj8uTJ0rFjx8DnRAIIIIBAGATo1/SryACiNzSLQEGbURIIgcALJCUlpX77+bhx41K/FT3wiZEAAgggEHAB+jX9AjKA6A3NIlDQZpQEQiDwAuvXr5cqVap4eYwaNUp++9vfBj4nEkAAAQTCIEC/pl9FBhC9oVkECtqMkkAIBF7gyy+/lJtuusnLY8SIEdKrV6/A50QCCCCAQBgE6Nf0q8gAojc0i0BBm1ESCIHAC2zdulXKly/v5fHSSy/J888/H/icSAABBBAIgwD9mn4VGUD0hmYRKGgzSgIhEHiBXbt2SalSpbw8Bg8eLL/73e8CnxMJIIAAAmEQoF/TryIDiN7QLAIFbUZJIAQCL7B//34pVqyYl8cLL7wgQ4YMCXxOJIAAAgiEQYB+Tb+KDCB6Q7MIFLQZJYEQCLzA4cOH5ZprrvHyeO655+QPf/hD4HMiAQQQQCAMAvRr+lVkANEbmkWgoM0oCYRA4AW+/fZbufLKK708nn76aRk5cmTgcyIBBBBAIAwC9Gv6VWQA0RuaRaCgzSgJhEDgBX744QcpUKCAl0f37t1l9OjRgc+JBBBAAIEwCNCv6VeRAURvaBaBgjajJBACgRc4d+6c5M2b18vjiSeekIkTJwY+JxJAAAEEwiBAv6ZfRQYQvaFZBArajJJACIRC4NJLL5Uff/xR2rdvL6+//noociIJBBBAIOgC9Gv6FWQA0RuaRaCgzSgJhEAoBC6//HI5deqUPPLII/LnP/85FDmRBAIIIBB0Afo1/QoygOgNzSJQ0GaUBEIgFAKFCxeW48ePy69//WuZNWtWKHIiCQQQQCDoAvRr+hVkANEbmkWgoM0oCYRAKASKFCki//3vf+Whhx6SuXPnhiInkkAAAQSCLkC/pl9BBhC9oVkECtqMkkAIhELg+uuvl2+++UaaNGki8+fPD0VOJIEAAggEXYB+Tb+CDCBpDNeuXSszZsyQpUuXyvbt26VgwYJyyy23yPPPPy8NGjRIp+0+oWb48OEyefJk2b17t5QoUUI6deokvXv3ljx58sS0MhR0TGy8CIHQCrifK3v27JF77rlHPvjgg9DmSWIIIIBAkATo1/SrxQCSxrBly5ayfPlyadGihVSvXt2793rq1KmyceNGGTdunHTp0iV1765du8r48eO9T6epXbu2rFq1ytvX/fexY8fGtDIUdExsvAiB0AqULVvWezPkrrvu8t4YYUMAAQQQyHkB+jX9GjCApDFcuXKl3HrrrZI/f/7U/3ry5EmpWrWqHDx4UA4cOOB9Lv+GDRukSpUq0qNHDxk1alTqvj179vS+LGzdunVSqVKlqFeHgo6ajBcgEGqBihUryldffSV33HGH/P3vfw91riSHAAIIBEWAfk2/UgwgERj26tVLRo4cKbt27fJuterfv7+89NJLsm3bNilTpkxqBPdOpXvHsl+/fjJ06NAIIqffhYKOmowXIBBqAfdGhrsCe9ttt8k///nPUOdKcggggEBQBOjX9CvFABKB4cMPPyyzZ8+Wo0ePes+FNGrUyLvKsX///gteXbRoUalWrZosWrQogsgMIFEj8QIEcpGAuxX0888/936mfPbZZ7koc1JFAAEEEleAAUS/NgwgFzHctGmTdwuW+xSaOXPmeHu7dyUvu+wySUpKuuDVrmE4c+aMd5tWdtu+ffvE/Um7LVy4UAYMGOA9CN+2bVv96hIBAQQCLXD77bfLmjVr5Be/+MVFf6YEOlFOHgEEEAiQAAOIfrEYQLIxPHbsmNSqVcsbFNwVj5IlS3p7lytXTtyVDvfgecbNPZDunhXZsmVLtqszaNAgGTx4cKb7MIDoC5sICIRB4M4775QVK1bIz3/+c+9ZEDYEEEAAgZwXYADRrwEDSBaG7uFzd6vVp59+6t1OVbdu3dQ9uQKiLzwiIIDAxQXq168vy5Yt8541c8+csSGAAAII5LwAA4h+DRhAMjE8ffq0NG3a1PvYS/ftw+72q7Qbz4DoC48ICCBwcQH3s2bx4sVSvHhx7/uG2BBAAAEEcl6AAUS/BgwgGQzPnj0r7vtA3LcOuwJr3br1BcruU66GDRvGp2Dp648ICCCQjcADDzwgCxYskOuuu877RnQ2BBBAAIGcF2AA0a8BA0gaw/Pnz0ubNm1k1qxZMmnSJHn88cczFXbPg7hPpcnqe0C++OILqVy5ctSrQ0FHTcYLEAi1QPPmzeWdd96Rq666So4cvJV1agAAIABJREFUORLqXEkOAQQQCIoA/Zp+pRhA0hg+88wz8tprr3nPe3Tq1OkC3YYNG3oPn7vtySeflIkTJ3rfhF6nTh1xX2Lovgm9c+fOMmHChJhWhoKOiY0XIRBagVatWsnbb78thQoVku+++y60eZIYAgggECQB+jX9ajGApDGsV6+eLF++PEtV9zCo28dt7latV155RSZPnix79uzx7tF2Q8tzzz3nfVt6LBsFHYsar0EgvAKPPvqo97Hc+fPnlx9++CG8iZIZAgggECAB+jX9YjGA6A3NIlDQZpQEQiAUAh06dPCurF5yySXibhFlQwABBBDIeQH6Nf0aMIDoDc0iUNBmlARCIBQC7pZO9zya29xV1zx58oQiL5JAAAEEgixAv6ZfPQYQvaFZBArajJJACIRCwH3QxZgxY7xc3HcTXX755aHIiyQQQACBIAvkRL/mPpJ99uzZkpSUJBs2bJAzZ87I9u3bpXTp0oGkZABJoGXLiYJOoPQ5FQQQyCCQ8sEY7j8fO3ZMfvKTn2CEAAIIIJDDAjnRrz322GMyc+ZMcV+GferUKdm4cSMDSA7XQWgOnxMFHRo8EkEghAJ9+vTxPuzCbYcOHZKrr746hFmSEgIIIBAsgZzo1/7zn//INddc430oyQsvvCBDhw5lAAlW2STu2eZEQSeuBmeGAAIDBgyQF1980YPYt2+fXH/99aAggAACCOSwQE73awwgOVwAYTt8Thd02DzJB4GgC/z+97+XgQMHemns3LlTSpYsGfSUOH8EEEAg8AI53a8xgAS+hBIrgZwu6MTS4GwQQGDYsGHSr18/D2LLli1Srlw5UBBAAAEEclggpV8bMmSING7cON3ZFCtWTNwfPzcGED91c2FsBpBcuOikjEA2Aq+++qo8++yz3h6bNm2SihUr4oUAAgggkMMCKf1aZqfhrloPGjQoyzN0n161devWiDK44oorMr3yzQASER87RSrAABKpFPshkDsE/vSnP0nPnj29ZNetWyeVK1fOHYmTJQIIIJDAAporIDt27JAyZcpElF3dunXl448/vmBfBpCI+NgpUgEGkEil2A+B3CEwfvx46dq1q5fs2rVrpUaNGrkjcbJEAAEEElhA06+573Ravnx5RNkVKVJEbrvtNgaQiLTYKWYBTUHHfFBeiAACCSswZcoU6dSpk3d+//jHP6RmzZoJe66cGAIIIJBbBHK6X+MKSG6ptDjlmdMFHac0OQwCCEQo8MYbb0i7du28vd07Zr/85S8jfCW7IYAAAgj4JZDT/RoDiF8rm0vj5nRB51J20kYgYQXct962adPGO7+PPvpI7r777oQ9V04MAQQQyC0COdGvrV+/XubNm5f6+8C9KdWrVy/56U9/6v3p3r17oPgv+fHHH38M1BmH+GRzoqBDzElqCAReYM6cOdKyZUsvj/fff1/uvffewOdEAggggEDQBXKiX5s2bZq0b98+U7pSpUqJe7g9SBsDSAKtVk4UdAKlz6kggEAGAfduV7Nmzbz/+u6770rTpk0xQgABBBDIYQH6Nf0CMIDoDc0iUNBmlARCIBQC7qpHypdczZ49W1q0aBGKvEgCAQQQCLIA/Zp+9RhA9IZmEShoM0oCIRAKgSVLlkiDBg28XNzzIK1btw5FXiSBAAIIBFmAfk2/egwgekOzCBS0GSWBEAiFwCeffCLui6jcNn36dPnNb34TirxIAgEEEAiyAP2afvUYQPSGZhEoaDNKAiEQCgH33R+1a9f2cpk8ebJ07NgxFHmRBAIIIBBkAfo1/eoxgOgNzSJQ0GaUBEIgFAJJSUly6623erm4b0V/8sknQ5EXSSCAAAJBFqBf068eA4je0CwCBW1GSSAEQiHgPve9SpUqXi6jRo2S3/72t6HIiyQQQACBIAvQr+lXjwFEb2gWgYI2oyQQAqEQ2LRpk9x8881eLiNGjPC+dIoNAQQQQCBnBejX9P4MIHpDswgUtBklgRAIhcCWLVvkZz/7mZfLsGHDpG/fvqHIiyQQQACBIAvQr+lXjwFEb2gWgYI2oyQQAqEQ2LVrl7hvuHXb4MGD5Xe/+10o8iIJBBBAIMgC9Gv61WMA0RuaRaCgzSgJhEAoBPbt2yc33HCDl8sLL7wgQ4YMCUVeJIEAAggEWYB+Tb96DCB6Q7MIFLQZJYEQCIXAoUOH5Nprr/Vy6dOnj7z88suhyIskEEAAgSAL0K/pV48BRG9oFoGCNqMkEAKhEPj222/lyiuv9HJ5+umnZeTIkaHIiyQQQACBIAvQr+lXjwFEb2gWgYI2oyQQAqEQOHnypFxxxRVeLt27d5fRo0eHIi+SQAABBIIsQL+mXz0GEL2hWQQK2oySQAiEQuDs2bOSL18+L5fOnTvLhAkTQpEXSSCAAAJBFqBf068eA4je0CwCBW1GSSAEQiHw448/Sp48ecT9b/v27eX1118PRV4kgQACCARZgH5Nv3oMIHpDswgUtBklgRAIjUD+/Pnl9OnT8sgjj8if//zn0ORFIggggEBQBejX9CvHAKI3NItAQZtREgiB0AgUKlRIvv/+e2nVqpX85S9/CU1eJIIAAggEVYB+Tb9yDCB6Q7MIFLQZJYEQCI1AkSJF5L///a80b95c5syZE5q8SAQBBBAIqgD9mn7lGED0hmYRKGgzSgIhEBqBokWLyoEDB6RJkyYyf/780ORFIggggEBQBejX9CvHAKI3NItAQZtREgiB0AgUL15c/vOf/0ijRo1k0aJFocmLRBBAAIGgCtCv6VeOAURvaBaBgjajJBACoREoW7asbN++XerXry9LliwJTV4kggACCARVgH5Nv3IMIHpDswgUtBklgRAIjUCFChVk8+bNcscdd8jf//730ORFIggggEBQBejX9CvHAKI3NItAQZtREgiB0Aj84he/kH/9619y++23y+rVq0OTF4kggAACQRWgX9OvHAOI3tAsAgVtRkkgBEIjUL16dfn888/F/W9SUlJo8iIRBBBAIKgC9Gv6lWMA0RuaRaCgzSgJhEBoBNyVjzVr1oi7ErJhw4bQ5EUiCCCAQFAF6Nf0K8cAojc0i0BBm1ESCIHQCLhnP1auXCnuWZAvv/wyNHmRCAIIIBBUAfo1/coxgOgNzSJQ0GaUBEIgNALu06+WLVsm7tOwtm7dGpq8SAQBBBAIqgD9mn7lGED0hmYRKGgzSgIhEBoB9/0fixcvFvd9ILt37w5NXiSCAAIIBFWAfk2/cgwgekOzCBS0GSWBEAiNgPsG9Pfee0/cN6Lv378/NHmRCAIIIBBUAfo1/coxgOgNzSJQ0GaUBEIgNALNmzeXd955R4oUKSKHDx8OTV4kggACCARVgH5Nv3IMIHpDswgUtBklgRAIjUCrVq3k7bfflkKFCsl3330XmrxIBAEEEAiqAP2afuUYQPSGZhEoaDNKAiEQGoFHHnlE3M+G/Pnzyw8//BCavEgEAQQQCKoA/Zp+5RhA9IZmEShoM0oCIRAagQ4dOsjUqVPl0ksvlXPnzoUmLxJBAAEEgipAv6ZfOQYQvaFZBArajJJACIRGoHPnzjJp0iQvn7Nnz0qePHlCkxuJIIAAAkEUoF/TrxoDiN7QLAIFbUZJIARCI9C9e3cZO3asl8/Jkyfl8ssvD01uJIIAAggEUYB+Tb9qDCB6Q7MIFLQZJYEQCI3AM888I6+99pqXz7fffiuFCxcOTW4kggACCARRgH5Nv2oMIHpDswgUtBklgRAIjUCfPn3klVde8fI5dOiQXH311aHJjUQQQACBIArQr+lXjQFEb2gWgYI2oyQQAqEReOGFF2To0KFePvv27ZPrr78+NLmRCAIIIBBEAfo1/aoxgOgNzSJQ0GaUBEIgNAK///3vZeDAgV4+u3btkhIlSoQmNxJBAAEEgihAv6ZfNQYQvaFZBArajJJACIRGYNiwYdKvXz8vny1btki5cuVCkxuJIIAAAkEUoF/TrxoDiN7QLAIFbUZJIARCIzBixAjp3bu3l8+mTZukYsWKocmNRBBAAIEgCtCv6VeNAURvaBaBgjajJBACoREYNWqUPPXUU14+69evl0qVKoUmNxJBAAEEgihAv6ZfNQYQvaFZBArajJJACIRGYPz48dK1a1cvn7Vr10qNGjVCkxuJIIAAAkEUoF/TrxoDiN7QLAIFbUZJIARCIzB58mR5/PHHvXz+8Y9/SM2aNUOTG4kggAACQRSgX9OvGgOI3tAsAgVtRkkgBEIjMH36dHnssce8fD755BO58847Q5MbiSCAAAJBFKBf068aA4je0CwCBW1GSSAEQiMwc+ZMadOmjZfPRx99JHfffXdociMRBBBAIIgC9Gv6VWMA0RuaRaCgzSgJhEBoBGbPni2/+tWvvHzef/99uffee0OTG4kggAACQRSgX9OvGgOI3tAsAgVtRkkgBEIj8O6778qDDz7o5TNv3jx54IEHQpMbiSCAAAJBFKBf068aA0gGw+PHj4v73P2kpCTvE2f2798v7dq1k2nTpqXbc8eOHVKmTJlMV6Bjx47iHhyNdqOgoxVjfwTCL+CuejRu3NhLdM6cOdK8efPwJ02GCCCAQAIL0K/pF4cBJINhymBRrFgx7+MuFyxYkO0A0qxZM2nZsmW6KOXLl4/pk2ooaH1BEwGBsAm45z4aNmzopeWeB2ndunXYUiQfBBBAIFAC9Gv65WIAyWB46tQpOXTokNx4441y9uxZyZcvX7YDSP/+/eXFF1/Ur4SIUNAmjARBIFQCy5cvl3r16nk5vfHGG/Loo4+GKj+SQQABBIImQL+mXzEGkGwMIx1A3BDitgIFCqhWhIJW8fFiBEIp4L77o3bt2l5uU6ZMkQ4dOoQyT5JCAAEEgiJAv6ZfKQYQ5QBSqFAhcc+NuM3devXUU09Jt27dYloZCjomNl6EQKgF3LNo//d//+fl2KVLF7nvvvtCnS/JIYAAAoku8PHHH8vIkSNlxowZ0rZt20Q/3YQ8PwaQGAeQXbt2ee9Euk+nKVWqlOzdu9d78Nw1C88++6wMHz482wXft2+fuD9pt4ULF8qAAQMo6IT8p8JJIZAzAuvWrZOqVavmzME5KgIIIIBAlgIMILEXBwNIjANIZi87d+6c1K9fX1asWCGbN2+WcuXKZRl90KBBMnjw4Ez/noKOvaB5JQJhEzhx4oSULFlSDh8+HLbUyAcBBBAItAD9WuzLxwBiOIC4UPPnz5emTZvKxIkT5YknnsgyOldAYi9aXolAbhNww4d7FuTHH3/MbamTLwIIIJBwAtyCpV8SBhDjAWT9+vVSpUoVGTp0qPTr1y+qFeIZkKi42BkBBBBAAAEEEIi7AP2anpwBxHgAmTt3rrRo0cJ7HsR9IWE0GwUdjRb7IoAAAggggAAC8RegX9ObM4DEOIAcOXJEihQpku7VJ0+elDp16sjGjRtl69atUqJEiahWiIKOioudEUAAAQQQQACBuAvQr+nJGUAyMRwzZowcPXpUzp8/LwMHDpRq1apJ8+bNvT3d8x2VK1f2/n/3cGjNmjWlePHi3qdgTZ8+XbZt2ybDhg2Tvn37Rr06FHTUZLwAAQQQQAABBBCIqwD9mp6bASQTw9KlS8vOnTsz1Z06dao89thj3heCuYHjq6++Enc1xH0fSPXq1aVnz57ekBLLRkHHosZrEEAAAQQQQACB+AnQr+mtGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CNCv6VkZQPSGZhEoaDNKAiGAAAIIIIAAAr4I0K/pWRlA9IZmEShoM0oCIYAAAggggAACvgjQr+lZGUD0hmYRKGgzSgIhgAACCCCAAAK+CMS7Xztx4oS88cYbMm/ePNmwYYMcPnxYSpcuLU2aNJF+/frJT3/6U1/y9DMoA4ifulHGjndBR3l67I4AAggggAACCOR6gXj3axs3bpTKlSvLnXfeKY0aNZLrrrtOkpKSZPLkyd4g4v7vn/zkJ4FaFwaQBFqueBd0AqXOqSCAAAIIIIAAAoEQiHe/dujQIdm7d683hKTdXn/9denYsaO8+uqr8swzzwTCLuUkGUASaLniXdAJlDqnggACCCCAAAIIBEIgUfq1b7/9Vq688krp0KGDTJkyJRB2DCAJuEyJUtAJSMMpIYAAAggggAACCSGQKP3aV199JRUrVpS+ffvKsGHDEsIm0pPgCkikUnHYL1EKOg6pcggEEEAAAQQQQCCQAin92pAhQ6Rx48bpcihWrJi4P/HYfvOb38iMGTPk888/lypVqsTjkGbHYAAxo9QHYgDRGxIBAQQQQAABBBDwUyClX8vsGAMHDpRBgwZlefgzZ87I1q1bIzq9K664QkqWLJnpvpMmTZLOnTt7z364Z0CCtjGAJNCKMYAk0GJwKggggAACCCCAQCYCmisgO3bskDJlykTkWrduXfn4448v2Pdvf/ubtGzZ0rv6MnfuXMmbN29E8RJpJwaQBFoNBpAEWgxOBQEEEEAAAQQQyGYAcbc/tW3bNiqjkydPyvLlyyN6TZEiReS2225Lt+/ixYuladOmUrt2bVm4cKFcfvnlEcVKtJ0YQBJoRRhAEmgxOBUEEEAAAQQQQMB4ANGAusHlvvvuk0qVKsmSJUukUKFCmnA5+loGkBzlT39wBpAEWgxOBQEEEEAAAQQQSJAB5J///Kc0aNDAu33LDSJXXXVVoNeGASSBlo8BJIEWg1NBAAEEEEAAAQQSYADZuXOnVKtWTb777jt5+eWXpWjRounOyv3/DRs2DNRaMYAk0HIxgCTQYnAqCCCAAAIIIIBAAgwg7kH0u+66K8u1yOph9URePAaQBFodBpAEWgxOBQEEEEAAAQQQSIABJIyLwACSQKvKAJJAi8GpIIAAAggggAACDCC+1AADiC+ssQVlAInNjVchgAACCCCAAALxEqBf00szgOgNzSJQ0GaUBEIAAQQQQAABBHwRoF/TszKA6A3NIlDQZpQEQgABBBBAAAEEfBGgX9OzMoDoDc0iUNBmlARCAAEEEEAAAQR8EaBf07MygOgNzSJQ0GaUBEIAAQQQQAABBHwRoF/TszKA6A3NIlDQZpQEQgABBBBAAAEEfBGgX9OzMoDoDc0iUNBmlARCAAEEEEAAAQR8EaBf07MygOgNzSJQ0GaUBEIAAQQQQAABBHwRoF/TszKA6A3NIlDQZpQEQgABBBBAAAEEfBGgX9OzMoDoDc0iUNBmlARCAAEEEEAAAQR8EaBf07MygOgNzSJQ0GaUBEIAAQQQQAABBHwRoF/TszKA6A3NIlDQZpQEQgABBBBAAAEEfBGgX9OzMoDoDc0iUNBmlARCAAEEEEAAAQR8EaBf07MygOgNzSJQ0GaUBEIAAQQQQAABBHwRoF/TszKA6A3NIlDQZpQEQgABBBBAAAEEfBGgX9OzMoDoDc0iUNBmlARCAAEEEEAAAQR8EaBf07MygOgNzSJQ0GaUBEIAAQQQQAABBHwRoF/TszKA6A3NIlDQZpQEQgABBBBAAAEEfBGgX9OzMoDoDc0iUNBmlARCAAEEEEAAAQR8EaBf07MygOgNzSJQ0GaUBEIAAQQQQAABBHwRoF/TszKA6A3NIlDQZpQEQgABBBBAAAEEfBGgX9OzMoDoDc0iUNBmlARCAAEEEEAAAQR8EaBf07MygOgNzSJQ0GaUBEIAAQQQQAABBHwRoF/TszKA6A3NIlDQZpQEQgABBBBAAAEEfBGgX9OzMoDoDc0iUNBmlARCAAEEEEAAAQR8EaBf07MygOgNzSJQ0GaUBEIAAQQQQAABBHwRoF/TszKA6A3NIlDQZpQEQgABBBBAAAEEfBGgX9OzMoDoDc0iUNBmlARCAAEEEEAAAQR8EaBf07MygOgNzSJQ0GaUBEIAAQQQQAABBHwRoF/TszKA6A3NIlDQZpQEQgABBBBAAAEEfBGgX9OzMoDoDc0iUNBmlARCAAEEEEAAAQR8EaBf07MygGQwPH78uIwYMUKSkpJk7dq1sn//fmnXrp1MmzbtAu1z587J8OHDZfLkybJ7924pUaKEdOrUSXr37i158uSJenUo6KjJeAECCCCAAAIIIBBXAfo1PTcDSAbDHTt2SJkyZaRYsWJSo0YNWbBgQZYDSNeuXWX8+PHSvn17qV27tqxatUqmTp0q7r+PHTs26tWhoKMm4wUIIIAAAggggEBcBejX9NwMIBkMT506JYcOHZIbb7xRzp49K/ny5ct0ANmwYYNUqVJFevToIaNGjUqN0rNnTxk9erSsW7dOKlWqFNUKUdBRcbEzAggggAACCCAQdwH6NT05A0g2htkNIP3795eXXnpJtm3b5l0xSdm2b98uZcuWlX79+snQoUOjWiEKOioudkYAAQQQQAABBOIuQL+mJ2cAiXEAadSokXeVwz0jknErWrSoVKtWTRYtWhTVClHQUXGxMwIIIIAAAgggEHcB+jU9OQNIjAOIu73qsssu8x5Wz7hVr15dzpw5I+42ray2ffv2ifuTdps9e7YMGzZMunTpInXq1NGvLhEQQAABBBBAAAEETAVWrlzpPQM8btw4r2dji16AASTGAaRcuXLirnS4B88zbu6B9AMHDsiWLVuyjD5o0CAZPHhw9CvGKxBAAAEEEEAAAQRyXMANH24IYYtegAEkxgHEjysg7vmRefPmSa1ataRw4cLRr2aIXuEsBgwYIEOGDEn3jE2IUow6FUwyJ8PlQhdMMInkBwx1Qp1EUiduH2olvdR3330n69evF/fBQxUqVIiUkf3SCDCAxDiA+PEMCJX5P4HPPvvM+xhkd4ubu6WNTQSTzKsAlwtdMMEkkp+Z1Al1EkmduH2olUil2C9SAQaQGAcQ9ylX7nkNy0/BinTRcsN+/LDjF2OkdU6tUCuR1Ap1Qp1QJ5EI8EZP7Eq8MhoBBpAYBxD3CVjuk66y+h6QL774QipXrhzNWrBvGgGaBZqFSP9BUCvUSiS1Qp1QJ9RJJAIMILEr8cpoBBhAMtEaM2aMHD16VM6fPy8DBw70Bo3mzZt7ezZt2jR1sHjyySdl4sSJ3jehu0+tcp+K4L4JvXPnzjJhwoRo1oF9MwjQLNAsRPqPglqhViKpFeqEOqFOIhFgAIldiVdGI8AAkolW6dKlZefOnZk6ugHjscce8/7OfVHhK6+8IpMnT5Y9e/ZI8eLFpVOnTvLcc89J3rx5o1kH9s0g4D6i2A13bpgrVqwYPiLexzZjcmEp4IJJJD8gqBPqhDqJRCDzffj3E7sdr8xcgAGEykAAAQQQQAABBBBAAIG4CTCAxI2aAyGAAAIIIIAAAggggAADCDWAAAIIIIAAAggggAACcRNgAIkbNQdCAAEEEEAAAQQQQAABBhBqAAEEEEAAAQQQQAABBOImwAASN2oOhAACCCCAAAIIIIAAAgwg1EDCCOzYsUPKlCmT6fl07NjR+7jjMG/Hjx+XESNGSFJSkqxdu1b2798v7dq1k2nTpl2Q9rlz52T48OGeye7du6VEiRLeR0D37t1b8uTJEyqmSF1yS/242pgxY4YsXbpUtm/fLgULFpRbbrlFnn/+eWnQoEG6tc8tdRKpSW6pkZQi2LRpkwwePNj7meI+RvXSSy+VsmXLeh8l36VLF8mfP39qveSWWonUJLfVSsZfGu7ny9133+3956+//lrKly+fusvJkydl0KBB8tZbb8nBgwe9v+vZs6c8/vjjofrdQzL+CjCA+OtL9CgEUn7gN2vWTFq2bJnule4HXM2aNaOIFrxdU/J333tSo0YNWbBgQZYDSNeuXWX8+PHel2DWrl1bVq1a5X0JpvvvY8eODV7y2ZxxpC65pX7cv43ly5dLixYtpHr16uIGNLf2GzdulHHjxnmNZcqWW+okUpPcUiMp67948WLvTY3bb7/d+54qN2S4L8ydOXOm3H///TJ//vxcVyuRmuS2Wkn7I/j06dNSpUoV782t77///oIBxNWOc+zevbvcfPPN8t5778m7777rfS+aexOMDYFIBBhAIlFin7gIpPzA79+/v7z44otxOWYiHeTUqVNy6NAhufHGG70vucyXL1+mA8iGDRu8Xw49evSQUaNGpabg3oEaPXq0rFu3TipVqpRIqanOJVKX3FI/roG89dZb07177d6RrFq1qvdu5IEDB7wvQs1NdRKpSW6pkYv9g3ONo3uj4ssvv5QKFSrkqlrJyiajSW6ulWHDhskf//hHadOmjfe/aa+AuDfGHnjgARk5cqQ8/fTTqZzujcMPP/zQ+xLna6+99mIlyN8jIAwgFEHCCKT9ge+GELcVKFAgYc4vnieS3QDibF566SXZtm1bulvW3O047vaKfv36ydChQ+N5unE7VnYuub1+evXq5TUFu3bt8m7Jy811klKQGU1ye42kuLirIu6d6tWrV3tXR6gV8a4UpTXJrbXiBgh3VWPMmDHeMOFu4Us7gLRt21bmzp0rR44cSff7edmyZVK/fn2ZNGkSt2LF7TdisA+k0IwMAAAIY0lEQVTEABLs9QvV2af8wC9UqJB3W4nb3K1XTz31lHTr1i1UuV4smewa7UaNGnlXOdwzIhm3okWLSrVq1WTRokUXO0Qg/z6SASS31s/DDz8ss2fPlqNHj3rPheTmOkkp7owmufVnzIkTJ8T9cbfTrFmzxvt56q6wbtmyxWsic2OtXMwkt9aKu5LhrqK623rd8JFxAHFXzK666ipveE27uauwV1xxhXTu3FkmTJgQyN8vnHR8BRhA4uvN0bIRcO/cdujQQR588EEpVaqU7N2713vI2j1g+uyzz3oPXeeWLbtG291eddlll3kPlmbc3DMBZ86c8W6pCOOWnUturh/3YK27BatJkyYyZ84cb+lzc524/DMzya014h4Ydo1kyuauerh3qitXrpxra+ViJrmxVtztVW4AcUOqew4xxSjtFZDChQvLPffck/pzJu3vmSJFikidOnXSPVsUxt9D5GQjwABi40gUnwTcQ5Pusu6KFStk8+bNUq5cOZ+OlFhhs2u0nYG70uHeocq4uQfS3btX7p3NMG7ZuWSWb26on2PHjkmtWrW8TzlyV8ZKlizpUeTmOsnKJLfWiLtd0/05fPiw9+lp69evF3eff7169XJtrVzMJLfViruC4T5Nr2HDhjJx4kQv/cwGEPcpi61atfI+ASvjdsMNN3i3b3300Udh/PVDTsYCDCDGoISzF3Cf1NK0aVPvh+ITTzxhf4AEjMgVkMwXJdoBxEUJc/24psHdPvPpp596t93VrVs3FS63XgHJzuT/27tDnTiXMAzA63oDhFQi4R56AdiikEg0Aiqa4Gua1JO0ghBAoBEI7gCLwGGa3kNz8k7ynxDY5XwnQxvYecY1nQy7Dx+7/zv/PzOL/tSXuUbmveevX7/ODg4OWmDd2NgY/m5ZjB6bjFYrnz9/brvoZaJvZWVlYQBxB+QVXjC80ZckgLzRX9xILzuzddn1KQurs8B6hGYNyMsFkGWtn2yVmWCeGe0sCs3jVw/biM/1/5fJos+OZa2RRe/358+fs/fv3/+7YcWItfLY5rHJSLWSx52zgUnWW+Y8qal9+/at7ax4dXU1W1tba32sARnhCuTvvEcB5O84+ykdArm4ypkHWQ+SAwlHaM8FkISwPD5hF6ynBzTOq41lrJ/UR86+yMz98fHxbHt7+8lbH61OKiaLPjuWsUae+5ycFljnzJjMeo9WK/NsHpuMVCs3Nzdt85LnWja2yOYw2Zr34uLCLlgjXIj84fcogPxhYMPXBbKtXxaxPWx5nCKL2nLI2t3dXdtedIT2XADJYxP5slh0Dki+TKbFpctm9ZzLKPXz+/fvdhFwenr67JaXI9VJ1WSUGpn+7rMebHV19cnHwP7+ftvU4/v37+2soZFqpWoyUq1kzVTO8Hjczs7OZufn5+0uSA6yzAYx0+OKi84BSZCbV3PL9l3k/fQLCCD9hkZ4IYGtra22VWROPM+HXW4L//jxo830Z8b/06dPL/STXu8w2Xs926jmgurw8LAFjbik5XGbKVjs7u62NTE5CT0BLQex5TTsZd0CseIySv3s7e2159Wz3uPh4xJTVWcRaTYpSBulTqomo9TIVAsfP35sC8+z2DyTN/lsuby8bI/UfPjwYZazG3Jo5Ui1UjUZrVbmfSvOW4Sefpubm62GMgmWRefZPSsnoY/yPf16ryDe1isTQN7W72upX+3R0VELHLe3t+32bs5zyLayOeE7F98jtDxnm8Of5rUEjJ2dnfZfuRPw5cuX9lja/f19C2y5GM3M5nRBsUxeFZdR6icXk9fX1wt/vbmonHY3GqVOqiaj1MhUHLlLlrscWePy69ev2bt372br6+ttF6NcPObfUxulVqomo9XK/wkgmSjMBNnJyUmrq+y4l+/pTIBpBKoCAkhVSj8CBAgQIECAAAECBLoFBJBuQgMQIECAAAECBAgQIFAVEECqUvoRIECAAAECBAgQINAtIIB0ExqAAAECBAgQIECAAIGqgABSldKPAAECBAgQIECAAIFuAQGkm9AABAgQIECAAAECBAhUBQSQqpR+BAgQIECAAAECBAh0Cwgg3YQGIECAAAECBAgQIECgKiCAVKX0I0CAAAECBAgQIECgW0AA6SY0AAECBAgQIECAAAECVQEBpCqlHwECBAgQIECAAAEC3QICSDehAQgQIECAAAECBAgQqAoIIFUp/QgQIECAAAECBAgQ6BYQQLoJDUCAAAECBAgQIECAQFVAAKlK6UeAAAECBAgQIECAQLeAANJNaAACBAgQIECAAAECBKoCAkhVSj8CBAgQIECAAAECBLoFBJBuQgMQIECAAAECBAgQIFAVEECqUvoRIECAAAECBAgQINAtIIB0ExqAAAECBAgQIECAAIGqgABSldKPAAECBAgQIECAAIFuAQGkm9AABAgQIECAAAECBAhUBQSQqpR+BAgQIECAAAECBAh0Cwgg3YQGIECAAAECBAgQIECgKiCAVKX0I0CAAAECBAgQIECgW0AA6SY0AAECBAgQIECAAAECVQEBpCqlHwECBAgQIECAAAEC3QICSDehAQgQIECAAAECBAgQqAoIIFUp/QgQIECAAAECBAgQ6BYQQLoJDUCAAAECBAgQIECAQFVAAKlK6UeAAAECBAgQIECAQLeAANJNaAACBAgQIECAAAECBKoCAkhVSj8CBAgQIECAAAECBLoFBJBuQgMQIECAAAECBAgQIFAVEECqUvoRIECAAAECBAgQINAtIIB0ExqAAAECBAgQIECAAIGqgABSldKPAAECBAgQIECAAIFuAQGkm9AABAgQIECAAAECBAhUBQSQqpR+BAgQIECAAAECBAh0Cwgg3YQGIECAAAECBAgQIECgKvAPnXYkI9oeBaAAAAAASUVORK5CYII=\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib notebook\n",
    "t0=time.time()\n",
    "\n",
    "fig, ax = plt.subplots()\n",
    "ax2=ax.twinx()\n",
    "timedata, tempdata,setdata,outdata = [], [], [],[]\n",
    "templn, = ax.plot([], [], 'r')\n",
    "setln, = ax.plot([], [], 'k')\n",
    "outln, = ax2.plot([],[],\"b\")\n",
    "\n",
    "def update(frame):\n",
    "    timedata.append(time.time()-t0)\n",
    "    tempdata.append(float(inst.query(\":TEMP:MEAS?\")[:-1]))\n",
    "    setdata.append(float(inst.query(\":TEMP:PID:SET?\")[:-1]))\n",
    "    outdata.append(float(inst.query(\":TEMP:PID:OUT?\")[:-1]))\n",
    "    templn.set_data(timedata, tempdata)\n",
    "    setln.set_data(timedata, setdata)\n",
    "    outln.set_data(timedata, outdata)\n",
    "    plt.xlim(timedata[0],timedata[-1])\n",
    "    return templn, setln,outln,\n",
    "\n",
    "ani = FuncAnimation(fig, update, blit=True)\n",
    "ax.set_ylim(10,40)\n",
    "ax2.set_ylim(-2.5,2.5)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(17, <StatusCode.success: 0>)"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "inst.write(\":TEMP:PID:SET 15\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1595511728.1071255"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "time.time()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}