Comparing responses notebook

create_stationxml_from_nrl-L22.ipynb

create_stationxml_from_nrl-sts2.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 74,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "%matplotlib notebook\n",
    "#%matplotlib inline\n",
    "from obspy import UTCDateTime, read_inventory\n",
    "from obspy.clients.nrl import NRL\n",
    "from obspy.io.xseed import Parser"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### NRL "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 75,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "nrl = NRL()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 76,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "('STS-2, 120 s, 1500 V/m/s, generation 2 electronics',\n",
       " 'http://ds.iris.edu/NRL/sensors/streckeisen/RESP.XX.NS083..BHZ.STS2_gen2.120.1500')"
      ]
     },
     "execution_count": 76,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "nrl.sensors['Streckeisen']['STS-2']['1500']['2 - installed 09/94 to 04/97']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 77,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "('RT130, gain 1, 200 sps',\n",
       " 'http://ds.iris.edu/NRL/dataloggers/reftek/RESP.XX.NR010..HHZ.130.1.200')"
      ]
     },
     "execution_count": 77,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "nrl.dataloggers['REF TEK']['RT 130 & 130-SMA']['1']['200']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 78,
   "metadata": {
    "collapsed": true,
    "scrolled": false
   },
   "outputs": [],
   "source": [
    "response = nrl.get_response(\n",
    "    datalogger_keys=['REF TEK', 'RT 130 & 130-SMA', '1', '200'],\n",
    "    sensor_keys=['Streckeisen', 'STS-2', '1500', '2 - installed 09/94 to 04/97'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 79,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Channel Response\n",
      "\tFrom M/S (Velocity in Meters per Second) to COUNTS (Digital Counts)\n",
      "\tOverall Sensitivity: 9.38896e+08 defined at 0.050 Hz\n",
      "\t9 stages:\n",
      "\t\tStage 1: PolesZerosResponseStage from M/S to V, gain: 1500\n",
      "\t\tStage 2: CoefficientsTypeResponseStage from V to COUNTS, gain: 629129\n",
      "\t\tStage 3: CoefficientsTypeResponseStage from COUNTS to COUNTS, gain: 1\n",
      "\t\tStage 4: CoefficientsTypeResponseStage from COUNTS to COUNTS, gain: 1\n",
      "\t\tStage 5: CoefficientsTypeResponseStage from COUNTS to COUNTS, gain: 1\n",
      "\t\tStage 6: CoefficientsTypeResponseStage from COUNTS to COUNTS, gain: 1\n",
      "\t\tStage 7: CoefficientsTypeResponseStage from COUNTS to COUNTS, gain: 1\n",
      "\t\tStage 8: CoefficientsTypeResponseStage from COUNTS to COUNTS, gain: 1\n",
      "\t\tStage 9: CoefficientsTypeResponseStage from COUNTS to COUNTS, gain: 1\n"
     ]
    }
   ],
   "source": [
    "print(response)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "### Local resp file"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 80,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "filename = 'RESP.RT130-STS2G2-200SPS.WHIH.4A.HHZ.txt'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 81,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "with open(filename) as buf:\n",
    "    inv = read_inventory(buf, format='RESP')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 82,
   "metadata": {
    "collapsed": true,
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "chan = inv[0][0][0]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 83,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "response2 = chan.response"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 84,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Channel Response\n",
      "\tFrom M/S (Velocity in Meters per Second) to COUNTS (Digital Counts)\n",
      "\tOverall Sensitivity: 9.39316e+08 defined at 0.600 Hz\n",
      "\t10 stages:\n",
      "\t\tStage 1: PolesZerosResponseStage from M/S to V, gain: 1492.57\n",
      "\t\tStage 2: ResponseStage from V to V, gain: 1\n",
      "\t\tStage 3: CoefficientsTypeResponseStage from V to COUNTS, gain: 629330\n",
      "\t\tStage 4: CoefficientsTypeResponseStage from COUNTS to COUNTS, gain: 1\n",
      "\t\tStage 5: CoefficientsTypeResponseStage from COUNTS to COUNTS, gain: 1\n",
      "\t\tStage 6: CoefficientsTypeResponseStage from COUNTS to COUNTS, gain: 1\n",
      "\t\tStage 7: CoefficientsTypeResponseStage from COUNTS to COUNTS, gain: 1\n",
      "\t\tStage 8: CoefficientsTypeResponseStage from COUNTS to COUNTS, gain: 1\n",
      "\t\tStage 9: CoefficientsTypeResponseStage from COUNTS to COUNTS, gain: 1\n",
      "\t\tStage 10: CoefficientsTypeResponseStage from COUNTS to COUNTS, gain: 1\n"
     ]
    }
   ],
   "source": [
    "print(response2)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "### Overlay plot"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 85,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "MIN_FREQ = .001"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 86,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "UNIT = 'VEL'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 87,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "mpl.get_websocket_type = function() {\n",
       "    if (typeof(WebSocket) !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert('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",
       "            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",
       "        this.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 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);\n",
       "        canvas.attr('height', height);\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'];\n",
       "    var y0 = fig.canvas.height - msg['y0'];\n",
       "    var x1 = msg['x1'];\n",
       "    var y1 = fig.canvas.height - msg['y1'];\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;\n",
       "    var y = canvas_pos.y;\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\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\"];\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",
       "    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 + '\">');\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 dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\">';\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",
       "        // select the cell after this one\n",
       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
       "        IPython.notebook.select(index + 1);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\n",
       "\n",
       "\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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4XuydB3hUxRbHf0kooUhHQECl9yaoWBBBJQoiShHB9sTCQ1AUBSQgTYooRZQioqLYwAYIFtAHoiIgggZB6SqgIFJCDwnZvO/c3U02YZPc3bt3S3Lmffme7J36O7O7/z0zcyYKTUpACSgBJaAElIASUAL5ikBUvhqtDlYJKAEloASUgBJQAkoAFYA6CZSAElACSkAJKAElkM8IqADMZwbX4SoBJaAElIASUAJKQAWgzgEloASUgBJQAkpACeQzAioA85nBdbhKQAkoASWgBJSAElABqHNACSgBJaAElIASUAL5jIAKwHxmcB2uElACSkAJKAEloARUAOocUAJKQAkoASWgBJRAPiOgAjCfGVyHqwSUgBJQAkpACSgBFYA6B5SAElACSkAJKAElkM8IqADMZwbX4SoBJaAElIASUAJKQAWgzgEloASUgBJQAkpACeQzAioA85nBdbhKQAkoASWgBJSAElABqHNACSgBJaAElIASUAL5jIAKwHxmcB2uElACSkAJKAEloARUAOocUAJKQAkoASWgBJRAPiOgAjCfGVyHqwSUgBJQAkpACSgBFYA6B5SAElACSkAJKAElkM8IqADMZwbX4SoBJaAElIASUAJKQAWgzgEloASUgBJQAkpACeQzAioA85nBdbhKQAkoASWgBJSAElABqHNACSgBJaAElIASUAL5jIAKwHxmcB2uElACSkAJKAEloARUAOocUAJKQAkoASWgBJRAPiOgAjCfGVyHqwSUgBJQAkpACSgBFYA6B5SAPQRqA9cCJeypPtdazwA/A6sAR665836GcsCNwPlAtA3DPQ58A/xmQ91apRJQAkog4ARUAAYcqVaYzwkUK1as2GdJSUlXXXbZZUnnn3++iI2gv89OnTrlWLNmTYGUlJQjSUlJbYFf86ldogoXLvxsSkrKk/Xq1TtdvXr16Ojo6EDbI+3QoUPCOzY2NnbdiRMn4oBj+ZS3DlsJKIEIIRDoD8IIGbZ2UwnYQ6B48eILmzRpErdkyZLYUqVK2dOIyVrPnj3LyJEjHc8///zR5OTkysBpk0XzUrb/litX7oWvv/66cIMGDWwd17Fjx+jcufOZNWvWrDx58qSIQE1KQAkogbAloAIwbE2jHYtAAqWio6MPbtu2LaZGjRph0f20tDSqVq16+q+//roTWBAWnQpiJ4oXL/7LlClTGj7wwANBaXXv3r1cdNFFDofDUQk4EJRGtREloASUgB8EVAD6AU2LKIFsCLQ5//zzP/3nn3+KhBOhhx56KHX27NkTgafCqV9B6Issv5/dtWtXVLVq1YLQnLOJiy666OTu3bu7AZ8HrVFtSAkoASXgIwEVgD4C0+xKIAcC7atXrz5/586dxcOJ0sCBA5k8efI0h8PxSDj1Kwh9ESF+6t9//6VcOTkDEpzUsGHDE5s3b/4P8FFwWtRWlIASUAK+E1AB6DszLaEEsiPQvkaNGvN37NjhVQCOGDGCV199Fdkr1rx5c6ZPn052+9ISExPp27cvn332GVFRUXTo0IFp06ZRsmRJo+2VK1fSpk0bihd3NiVLvaVLl2b37t3n9G3QoEFMmjQp3wrAgwcPUrZs2XO4+GIPKfzGG2+IkGbXrl0G9zvuuIMXXnjhnHobNWp0YtOmTSoA9XNCCSiBsCagAjCszaOdizAC2QrA559/3hBwn3/+ObI/cNSoUcydO5dt27ZRtGjRc4Ypgi8lJYX58+cb4q579+4UK1aMhQsXpgvAtm3bIgc9RCDmlFQAnisAfbXHpEmTmDFjBm+99RYtW7bkzJkzbN26laZNm6oAjLA3qXZXCSgBJwEVgDoTlEDgCGQrAKtXr86AAQPo16+f0VpqaiqVKlViypQp3HmnnM/ISOLFu/jii9m4cSMNGzY0Hsh/i9iQZ1WqVDE8gCIAk5OTiYmJUQHonYCxBOzNA+iLPY4fP84FF1xgiPH27dvnOlvUA5grIs2gBJRAGBBQARgGRtAu5BkCXgWgLPlKSJjVq1dz+eWXpw82Li6ORo0aMXGinM/ISJ988onh8Tt9OnPUltjYWD788ENuvvnmdAEoYlC8UVLP008/zTXXXHMOTPUAZvYA+mqPpUuXGsJPPIDiCZTyzZo1Y8KECTRu3Pgc3ioA88z7WQeiBPI0ARWAedq8OrggE/AqACU0yIUXXshvv/1GnTp10rske8hKlCjBK6+8kqmbb7/9NnJwY9++fZler1ixorEHrWfPnvzzzz8cOHDA2EMoQvHll19m2LBhrF279hxRogIwswD01R7vvPMOd999N61atTK8gLLXUvYPyhK+LAOfd955meykAjDI7zptTgkoAb8IqAD0C5sWUgJeCQTMAyji8NSpU5ka8fQAemtdDoVcddVVjBkzJtNjFYDWPICLFy+mU6dOfPHFF7Rr185g63A4DOH38ccfI55cz6QCUD8dlIASiAQCKgAjwUrax0gh4NMeQNlXJh49b3sAJW5dQkJCpj2Asuz4559/GnsAvaXrrrvOOKAwduxYFYBOAj7tAczOHm6PoacAlD2c4r1VARgpb03tpxJQAlkJqADUOaEEAkcgWwEo+/zkFPCnn36KHEAQL517CdHbKeCOHTsap4Bl+VFOAYtIlHwLFjgv81i2bBm1a9eWoMMkJSUZy8hPPfUUq1at4pJLLlEBmIsA9NUeXbp04dChQ7z//vtGKB45xS0ngmVZ3x2Kxw1dPYCBe0NpTUpACdhHQAWgfWy15vxHIMc4gCNHjmTWrFnIqdIWLVqkxwHcs2cP9evXN5YYZQlXksQBlBPDS5YsMcK8iCAUASleJ0kiIGfPnm2IEhGGcghk+PDhtG7d+hzqugTsPQ6gL/Y4ceIEjz32mOHxk1PXl156qXF4R+yWNakAzH9vfB2xEohEAioAI9Fq2udwJZCjAAxVp1UAeheAdtlDBaBdZLVeJaAEAklABWAgaWpd+Z2ACsDwmgHZ7gG0s5sqAO2kq3UrASUQKAIqAANFUutRAqACMLxmgQrA8LKH9kYJKIEwIqACMIyMoV2JeALXV65ceeHevXuLhdNIHnnkEce0adOmAE+GU7+C0JcCQIqc4q1cuXIQmnM2UatWrRM7duzoCSwOWqPakBJQAkrARwIqAH0EptmVQA4EqkRHR/+5f//+6PLly4cNqEaNGp3atGnTI8DrYdOpIHWkWLFie+bOnVulc+fOQWlRDu+UK1cuNTU1VSJ+7wxKo9qIElACSsAPAioA/YCmRZRAdgSKFSu2vmfPno1nzZpVQE7vhjpJ2Jhu3bqdTU1NrQAcDnV/gt1+dHT0mIYNGz7x7bffxrpPUNvVBwnX89hjj6W+9tprW0+ePNnArna0XiWgBJRAIAiE/hsqEKPQOpRA+BC4ODY2dvUFF1xQolu3bkXLli1rhHEJZhIhItfDffnll6dWr15dMDU1tSvwSTD7EEZtxRYpUuTz2NjYlrfffnuBatWqFZAwLoFMwvvw4cN89NFHp3bv3n3qzJkzVwLbA9mG1qUElIASCDSB4H4zBbr3Wp8SCE8Ccjlsx0KFCl1fqFChsqHoYmpq6qnTp0+vAz4EdoeiD2HUpuwFvDYmJuaWIkWKyDUqgVWAQHJy8pHk5OT/uYT20TAau3ZFCSgBJeCVgApAnRhKQAkoASWgBJSAEshnBFQA5jOD63DDgUDsUCjcMTA9SfkHTt0GOAJTXz6opQDxRHNLQEbq4B/OIidMUgNSn1aiBJSAEggSARWAQQKtzSiBDAIlNsPd9aGFRShngP9KHaXl9jiLleWf4oXZRGMacIHFIYvkW2LUUQY4YrE2La4ElIASCCoBFYBBxa2NKQEhUOx9eLAzTLG4F20ZcEsinJGYM2eVrUkCBZnPJXThJot7ASXIy3sc5SzCP8Vk65pNCSgBJRAWBFQAhoUZtBP5jEBrOG8pHCgMsRaG3ikJPpsKZ5+yUEl+LHoNhfiSgRSioIXhv8cZtjMVB4Mt1KJFlYASUAIhIaACMCTYtdF8TiAKim6GTjXhrYL+HUp9E3jwDKTUA37P5zx9HX4UBdlEHWrSmUJE+1oc+Nk473sGB/WBXX7UoEWUgBJQAiEloAIwpPi18XxMoBIUWQ23XQBzfRSB6eLvJmBFPmZoZegViWE19anMbRT0SQQ6xV8yDoT/ciud0LJKQAkogVARUAEYKvLarhIAlwhsVRFeLQxVc2GSBIxxwHMpkNJexYflKVSJGL7nYipxC4UpmUt9ssvvGxys4iwOhL/E/dOkBJSAEohIAioAI9Js2uk8RKAcFHkR0rrC+ALwSJT3JWHRGvclwaHtcOoe1yJkHsIQsqGUpSAv4aArN1CAy4jy6g2URd5FJHGKHaRwt/IPmb20YSWgBAJEQAVggEBqNUrAIoG2UOQNqFEOphWB1q7qZHtffDJ85ICUQcAMjTlnkbT34m0oyJuUphw3UYSLAfl0lOAu/yOZX3G4DntMV/628NdKlYASCDIBFYBBBq7NKYEcCMRCgcEQPQiaRkGNGPggBgp+BKefAPYqPVsJyJHsgcTwFBWIoiwF2Ew0MXxMCgOUv63stXIloASCTEAFYJCBa3NKwASBkhD9CBS6EJKeB7abKKNZAkdAdgM+QgwXkcpzyj9wYLUmJaAEwoeACkBrthB+cp/AcWvVaGkloASUgBJQAhFB4DzgbyAtInqrncyWgApAa5Ojsi4LWQOopZWAElACSiDiCFQB/oq4XmuHMxFQAWhtQpQAju7Zs4cSJeQ/7U3x8fGMGzfOr0Z8LZtb/pyeZ/fM2+tmXuvTpw8zZ870a9y+Fspt3DnV52vZ3PIr43Np58bMs4SZvMpYGfv6GeHOb2Z+ZVe3L2XN5A3WPD527BhVqxrhqmSbxDF/2Wm58CCgAtCaHQwBePTo0aAIwKVLlxIXF+dXj30tm1v+nJ5n98zb62ZeGzBgAJMnT/Zr3L4Wym3cOdXna9nc8ivjc2nnxsyzhJm8ylgZ+/oZ4c5vZn5lV7cvZc3kDdY8FgFYsqQRMFMFoL8TJ4zKqQC0ZoygCkBrXY3c0sEUgJFLyVrPlbE1fmZKK2MzlKzlUcbW+OVWunDfuiTP2KoCMDdQEfJcBaA1QxkCsG/fvnTs2NFv75y1LuT90i+88AKPPfZY3h9oCEeojO2Hr4yVsf0E7GtBvIztB/fAkSDBMdUDaB/p4NWsAtAaa/UAWuNnqvSsWbPo3bu3qbyayT8Cytg/br6UUsa+0PIvrzL2j5vZUuoBNEsqMvKpALRmJxWA1vhpaSWgBJSAEogQAioAI8RQJrupAtAkqGyyqQC0xk9LKwEloASUQIQQUAEYIYYy2U0VgCZB5SQAdQ+gNYhaWgkoASWgBMKbgO4BDG/7+NM7FYD+UMsoox5Aa/xMlR4yZAjjx483lVcz+UdAGfvHzZdSytgXWv7lVcb+cTNbqsmAh9k4xYjJqmFgzEIL43wqAK0ZRwWgNX6mSu/fv5+KFSuayquZ/COgjP3j5kspZewLLf/yKmP/uJktpXEAzZKKjHwqAK3ZSQWgNX5aWgkoASWgBCKEgArACDGUyW6qADQJKptsKgCt8dPSSkAJKAElECEEVABGiKFMdlMFoElQOQlAPQRiDWJupc1chZRbHfo8ZwLK2P4ZooyVsf0E7GtB5u/ixYuZPn26NKJ7AO1DHbSaVQBaQ60eQGv8TJWeP38+3bt3N5VXM/lHQBn7x82XUsrYF1r+5VXG/nEzW0o9gGZJRUY+FYDW7KQC0Bo/La0ElIASUAIRQkAFYIQYymQ387MAlLGvBOoBcq59uAczcTf1B9KAqcD72fBUAWhyomk2JaAElIASiGwCGgYmsu2Xtff5WQAKi8rAdUBNDwEoom410AIQPuuAlsBxL6ZXARiE90NSUhKxsbFBaCn/NqGM7be9MlbG9hOwtwW9CcRevsGuPb8LQOF9L1DDQwC2A7oAvV3GmOXyAP5PBWCwp6ezvQEDBjB58uTQNJ5PWlXG9htaGStj+wnY24IKQHv5Brt2FYDnCsAeQCMg3mWMZ4GfgPkqAIM9PZ3tqefEfu7KWBnbT8D+FnQe28tYBaC9fINdeyQLQNmn1xdoAhQHCgKOLABHAQ8AslS73pV/c5Y83jyAnYH/engAPwC+UgEY7Omp7SkBJaAElEC4EFABGC6WCEw/IlkA3gCUAYoCr3oRgAOBfsBNwE5gBHAPUBs45YFPBKDsAXza9ZqIxe+Ay4BoYC1wpe4BDMyE01qUgBJQAkogMgmoAIxMu2XX60gWgO4xtQaWexGAuwDZODbNlTEG2Ac8Drzjeu1d13KviEgRiSIWUwE9BRxG81xje9lvDGWsjO0nYH8LOo/tZawC0F6+wa49rwpA8eIlAle4PHhurkuBX4AnAwTaOAXcrl07GjRoYFTZpk0bVqxYwbhx49JPrsqHUqlSpYiLizPyyIXlU6dOZfz48endmDVrFk2aNKFlSzlwDFu2bGHRokUMHjw4Pc+ECRPo1KkTdevWNV5bs2YNCQkJ9O7tPq8CQ4YMoX///lSsWNHII9HbExMT0wMpyx6Z+Pj4iOqf9Ld169bKD/vs++ijj9KzZ0+dfza+P+6//34GDhyY796/wfz8k4Dx8tmalz7/gsnP2/fHypUrjX3Ykk6dOoV8V+lNIAFSECGuJq8KwCrAbleMv60ejOcBx4CHAsTdEIB6FVyAaGo1SkAJKAElEJYE9Cq4sDSLpU7lVQEYVA/g0aNHKVFCmtSkBJSAElACSiBvEtCbQPKWXfOqABQredsD+LeElfPYA2jVmuoBtErQRHlZMncv6ZjIrln8IKCM/YDmYxFl7CMwP7IrYz+gmSyiHkCToCIoWyQLQDmhK6Ff5BDI58B5rgMcya4r3GSfn5wC7uASg8Ncp4DrZDkFbMVcehOIFXomy8q+FM/9kiaLaTYfCChjH2D5mVUZ+wnOh2LK2AdYfmRVD6Af0MK4SCQLQAnfMscl9gSxjEXu7m0DfONiPtJ1o4eIwx+ziQNoxTzqAbRCT8sqASWgBJRARBBQD2BEmMmnTkayAPRpoDZlVg+gTWC1WiWgBJSAEggvAk/P+ZgxveSmVEq6DlSGVwe1Nz4RUAHoE65zMqsH0Bo/La0ElIASUAIRQEA8gO0H98CRcEQFYATYy0wXVQCaoZR9HvUAWuNnqrTEnfKMdWiqkGbyiYAy9gmXX5mVsV/YfCqkjH3C5XNmDQTtM7KwLqAC0Jp5VABa42eqtAS8dgfINlVAM/lMQBn7jMznAsrYZ2Q+F1DGPiPzqYAKQJ9whX1mFYDWTKRLwNb4aWkloASUgBKIAAK6BBwBRvKxiyoAfQSWJbt6AK3x09JKQAkoASUQIQTUAxghhjLZTRWAJkFlk00FoDV+pkrLvcju+49NFdBMPhNQxj4j87mAMvYZmc8FlLHPyHwqoALQJ1xhn1kFoDUT6RKwNX6mSk+YMIHBgwebyquZ/COgjP3j5kspZewLLf/yKmP/uJkpJUvAIyZPY+2yJZJdw8CYgRbmeVQAWjOQegCt8dPSSkAJKAElECEE9CaQCDGUyW6qADQJKptsKgCt8dPSSkAJKAElECEEVABGiKFMdlMFoElQKgCtgdLSSkAJKAElENkEVABGtv2y9l4FoDV76h5Aa/xMldZ9PaYwWcqkjC3hM1VYGZvCZCmTMraEL8fCehewfWxDVbMKQGvkdQnYGj9TpfVknylMljIpY0v4TBVWxqYwWcqkjC3hy7WwegBzRRRRGVQAWjOXCkBr/LS0ElACSkAJRAgBFYARYiiT3VQBaBJUNtlUAFrjp6WVgBJQAkogQghoHMAIMZTJbqoANAlKBaA1UFZK6/2eVuiZK6uMzXGykksZW6Fnrqwy9s4pMSmROT/N4fMdnzO742wuKnURL78MmzdD375Qt645vjYJwPOBdsB55noR0bkcwO/Al0BaqEeiAtCaBfQQiDV+pkrPmjWL3r17m8qrmfwjoIz94+ZLKWXsCy3/8irjzNz+/hsKlzpEtanVOJ58nNgCsSzsvpC4mnH06wfTp0OBAjByVCrxQ6KJispeEli8C7gOUNQlfhJdvaxQtGjRdcnJyVWLFClCbGxsju37NyPCq5TD4eDo0aMULFjQkZSU9IrD4egTyh6qALRGX5eArfHT0kpACSgBJRBAAmlpsHQpTJkCy5fD77/D5M0DaFqxKbc3uN0QgZIk33ffwSOPOihx9310vLwhA68amGNPfPAAxgD9iaIHBaKakpJWwKg4JiqZ1LTCQFSRIkWONm/e/LyXX36ZBg0aBJBAeFeVlJTEwoULue+++0hKSpoMPBGqHqsAtEZeBaA1flpaCSgBJaAEAkDgdMppVu1ZRe0C11OtmrPCLl1g/HioUSP7Bo6dPknL1y9l++Ht/Prwr9QqWyvbzCYFYHNiopZTIK0EdaKgThpUB/4B3iGNFKKB7oULF5538OBBihcvHoDRR14V06dP56mnnjp94sQJ8YyGJKkAzB77UKADkAI8DXzjJasKwJBMW21UCSgBJaAEhMDBUweZ9P0kZm+YzZGkI/ze/3dWLLyQa6+Fiy4yx+jbP7/lmjeuoVfTXrzW6TUrAnASMQygXhTcnAZOZ6Mz/Qm8nS4AF9900003f/bZZ+Y6mAdz/fXXX1SpUkVGVg44FIohqgD0Tr0xIK7Z64FSwHKguZdNmyoAgzBrhwwZwnj5GavJNgLK2Da06RUrY2UcKAKyt2/bNgyR99exv6g6pSo1ytSg36X96NWsF+cV9v08xQtrXqDNxW1oUrFJtt1sMuBhNk6ZKc9LAsc8MkYRHbWeAmnN6AR4W9HNLAC/7dWr19WvvZa92AwUq3CtJzk5mcKFZTUc0Ru/hKKfKgC9U+8GXAoMcj1eA9wDbMuSXQVgEGbt/v37qVixYhBayr9NKGP7ba+MlbG/BFIdqfy8/2e++S6VtydcxoYNUKEC7N3rPMTxw18/0OKCFkRHyeqqfSnbOIBRLKYwN/MQUCab9rMIwAceeODq2bNnn5M5MTGRwYMH8+mnnyLtXXHFFbz44ovUqSPnSLJPqamptGzZkvXr17Njxw6qV5d158CnPXv28PDDD7Ny5Urj4Er37t2ZMmUKBcQQwKlTp3jiiSf45JNPOHnyJNWqVWP48OHcdtttmTpz9uxZOQyiAjDwJrJcoxyKfwO4BqgEbALigO9VAFpmqxUoASWgBJRANgQcDpDDuO4DucOWD+PFtS8ap3ibl4zj92e+IC4ORE/In0t3BIVnNgLwJQrQj17ABTl0w6QA7NSpEykpKbzzzjsUK1bMEIMfffQRW7duRU4LZ5eeeeYZvv/+e5YtW8b27dv9FoCjRo3izz//5PXXXz+nqbS0NJo0aULz5s2ZNm0ahw8f5uabb6Zt27aGCJQk4u/zzz83+iFLvB988AF33nknGzdupK5HvB0VgPZN2e5AX0B82bLDVGS2xN/xTKOABwDx4q135d/skeF+4C7gb6AycK9rF4NnHeoBtM+GWrMSUAJKIE8T2Hd8H6v3ribhjz9ZsOwAhb8bz2+/wS+/kH6QQ/b3fbLtE66sciVxNdrT6sJWxMgZWxvSokWwejU8+6z3yr0IwCeIZiI9gOzPjjgrMyEAT58+zXnnncfq1au59FJZhIMzZ84YB0Xmzp1Ljx7S0Llpw4YNhiduwYIFNG7c+BwPoJSdPHkyf/zxB1WrVmXYsGFGfm8pJwEoXr927doh3vzSpUsbxcXTJwJPxKB49G699VZq1KjBpEmT0qsvX748M2bMoFs3WVx0JhWANkxgV5U3uBzRcrrmVS8CUM669wNuAnYCI1xLvLXFg5ulW/KbRupo76W7KgDts2F6zRJ/Kk5+8mqyjYAytg2tzmP70YaM8ZkzaRxIPEXp8wpRvKixnJeeDpw8wJx177Hoy4McOvMvF+8Yy+SxZfGMeDJ/03zu+OiO9DIVXjtOw9rFmTqVTPmChfDuu+Htt+HgQShb9txWswjAOkTzAx2BZiZ6aEIAypJpiRIlDE/e5ZdfblTqFoWPP/44zz///DkNyV66Fi1aGAKvVq1axpKr5xLwG2+8wciRIw1x2KxZM6Pum266yfDSXXnllefUl5MAlKXomTNn8puodFfat28flStXNjx8DRs25Ouvv2bgwIHMnz+fiy++mHnz5vHYY4+xadMmzj9fYl47kwpAE3PGYpbWrgMcWT2Au1yHPKa56pffU/uAx42D6s4kx5Nkh+ZxI56R8/dL1qQC0KKBzBSXN1J2v9bMlNc8GQQk9pf8uZN7qSkrY8kjS1HuFB2dsSTlyVPynT2b8Yp4JiRv1iT5zpxxvupIcyBbX9xeDM99S5Lv5ElISZXD9xAbC4UKRhMTndnl4XCk8e8RV4WyI714QWILZ8mT5uBU8mmOuMLOxsbEUrJEDIUKZfTurOMsp1JOGUzkS694wfOMYLQlS4JzfzacOXvGyCPJGMdR5y9/SeIEkD7Kc3eemKgYShYujWzUdyf5Mp330Rza3+r8HVkophCli5Q2GO/e7cHvvEPEFHKOq0iBIkYeg5nDGc9N2j9waj9lyicTW8RBycIl0/O484n36I/j25CxVa6SRu1KlShTJPPGrOSUVOZ+uQFHWhrCMq5FHapVkvNuGUnGs3r3D/z4YxoO0qhd9DJaXV4cj+8w4wTq2r1rOZuaxver02hSLI4YCtGqFVzgWg78I/EPY4+aLJ+lOWI4vb6rMQ75a9cOqlaFhP0JbNi3gTTSKFG4BJ1qdWXmTGce+evUCf6M+ppNBzYZ9VQpUYXb6t1GcjJMnJiRb+eZx6kbV9GYY40qNOLm2jcbA5J8I0Y4822ImUGNpvsoWzaNay++luury1k/Zzp9Gu68E36ofDfJ0UepUTuZp9s9SvtamX////nPUS6eWgUKnYSoNMbVW8aQ28XvkJE2H9hMw5kN018oMDuBz+Y05gaPbHuP7WXlHyupUtdyHKkAACAASURBVKIq5WNqUa9KxZAGQx47KZFhr3zHuy80pMdNF2caj/zDQwCWJibqT+pTgi4mL7QwIQClDfeP/bfeestYAn7qqacM0XX//fcjgb6zJlkiln2D8kyWbmXvn+cSsCzZPvLIIzzwgCz4OdNDD8lmRXjllVfOqS8nAThmzBhjb6J4KN1J4voVLVqU7777zhCU4gkUwff2228b+wJl2Vo8kLK07ZlUAJ6DPuAveBOAItrkK+EKYK1Hi0tdJ3Ge9KEXhgAUl7A7kGWbNm1YsWIF48aNMzaISpIv11KlSqVPbHEfT506NdPJVpm8MlFlE6ukLVu2sGjRImP/gztNmDDBmETufQRy7VFCQkKmWzLkpGH//v3TD02IZ0feHG4BJZM1Pj4+z/bvkksvQa49+mfXYaZMWURc3GCOHoUbb4Snp95LTL0YClUsRM0yNbk6+mrWrUtg8eLexheECJSjZS6j4LVncBR30KFWB9qktjH4dezY3Yil5XAkcbB8HUrdnkKhYmn0v7w/1f6qlm5fES+lS+8npWltuCGVwrFpTOvwIqnrUjPZd/W6n7hy8GVwpXPj8Mj6HxO7Z2Mm+36w9AN6zujJ2cYuYTP7B7q2eo+XXsqw74Q3JzD0s6HENIohWXTDxF2QNImPPx7Hbbc551/f5/oya9MsomtFG6LiWPxWYCowniVLoEMHuPXJW1lybAlRVaOoVaYW71z5MZdcsghwzj/JN+bDK1hbfC1R5aNoW60tz9R6xuD36KMZt7RUvK4m+5vuNC51uqfJPfQs0jOdX7FiUlMSNK4KHQ8aGzOGthpKo38bpfM7dQqKFdsPLavBjUlG2zM7zCRqfdS5/OIvgaud744RNZdQZO+mTPzmfzGfO6bfAS1c76CZP9O11bxM/J5941mGfDoE3N/Tz/4FSRMz8Xt4wsPM3DQTajpF2elh8vsxM79OT3bik8RPoCrULlubeVctOoffMx+0NPhRHoPf2NpjvfCrwf6mu3LhVwU6HvLKT+Zf8eK581u1dgNXD2tumd/4N8YT/2l8jvz6TOjDy5tezpHfLU/cwuKjiw1+Mv/mXfUJzZtnnn+jP7icH4r/kInfDz8k0L+/x/xrW4P9zc7ld/PN3XGGm5P5l8Ev/up4Gh9snD7/RAAWLbofLq8GNzgoUqgQL3V4gbPrzmaafz8l/EKrfl0oc3ljisQU55kOj/H76qWZ5t/K71byyTeLuSGuF7Url+ei88sybOjQsP58vq5LD5ZXm8d/G77AzP/0R74/ZNlTvjckyQEHlwibRRF6G64Rz1AvOX15mhSAEhtQPGhfffUVcnPGgw8+aOyjk6XVsWPHZmpBvHmy/PrLL78Yy8SyxCvLr54CUMSZCDH3IQ35ASEHRlq3bs3ixYuRQx2ybOy+BUU8jpJHysn/y+tLliwxxJ0ZD6BoACkrewgrVKjAqlWrjAMgb775puF5dCcVgD4oLT+zehOAEnhHfnPXA+Sb0J3muY61O38amEt6FZw5TgHLlXI2ld9+jWHLFhHJzkCnsqQyc91MRq4ciSy7yBdwwgOyYTij2YULYVmBvsz4cYbxovzyX3HvCkP0iUdDvELy13DofWxL+drw0HSp14Vx140z8ku+665zere21r+H8+vuoGSJaEPkPNQ8Y8pIvq5dYX2NbjhiTlKzRhTxN/Q9x4Nw9GQSdYd3Te/gszcN597rL8vEScI79F7cxzjxJ6npvqk880Q1mktAIlcSz8nwr4fjSMXId9lfb1HQUQL53eDOt/z35cYmckkFogoR8/H76eUl3yWXwAebP+DdTe8ar1cqXokp18/gXtn16kqDBsEPjpdZtnOZ8Uqj8xsxqs0og8t992XkK9NpHPujNhgfmq0vak2/y2SnhZPfgw86860vO5Aq9XdTqjTcVvc27miYsQQm+f77X1hV9kHORp2kTl14ok2vTF4aqePYyTNcMirjF/3YDgPo3jrzOpTsrxr85RDWun7mtTg+mqf6XEijRhn9/e3f35j4/UTDw7ZmLbQ6MZWCacXp0wcaukTh93u+542f5UwYhner4JfuhQPxJDjzfbb9MxZvXWzkOb/Y+Qy9ahRPPZXRjjgfElLfMwL1SpI5+ujljxo/PIZKxFFXKtZ6FkdifjX+dWnlS7mrsWxDdnqwhg93ZvqhwEQuarCfcmWjuK76ddxY88b08ikpzv1b36SNJ5UkGjSI4t4rbjZOiHqmU0kpdJ86mSj5X1QUj9/YmWsbZ44WLD+k3vzpbX5Y68xzWbHbuS2unOGxc6e/j//N0h1LcTiiWL8+iitLdKdwTCxXXw2V5Pgc4PYASlvRFCD2j9uM95F4n2WOli8P2w9tNwIRS55ihYpxVZVrWLUq4zBE/fpwMG0bYlPxGJeKLWV4+MRuP//sfF9KfamltpFW8LjR37JFyhp33kqSfDt3OvPsObGLkmWSKVY0inJFy1G2aMZap3gIJY/7iz8TtHzwj39P/sv5E8/nP03/w5xOczKNWMTgXfGDOLhho7z+ANG8auyiz+ngh2cNJgVgVsz//vsvF110keF5E3HlmcRb99xzzxmCy2lnB0eOHKFMmTLIkvHQoUMNj+Do0aO56y7neym3lJMH8JtvvjH2AMqyb9Y9gIcOHaJQoULGErZ4Lz09fl26dDGcMhL82Z1UAOZmCevPg+IBlLv9xOia7CGQsDeBWQmzWPf3On7991eSRx/h7BnnOt4bb2AIlXmb5vHSDy8ZS0M1S9dkTNuxyCEuMYss5zVrBqlF9iPR8mWJqUjBIhQtGLIA7PaAslCr/MJ3e6wtVKNFcyCgjO2fHsrYOuMKEysYP/C+uuercyrLdBNIVNQiSnItD6eBx7aKbHtgUgBu27bNEFdycEL28knIFflsksMWWdOJEyeMZWl3Em+ehI0Rr1ujRo0Mr6B47UR4yaliOb0rewbFYxgdHc0l8us3S8rtFLDsI5RyUq8s94rQu/baa9NPAXfo0IGYmBgkxqGMQVbq5KSw7F+U699UAFqfo2Zr8GUPoOzYGeCxB9BMG+oBNEPJhzz79jnvsaxZE8OLIKnnQz15r/J7lC9a3vBkXPnPXKpXLIecqJc/148/H1rRrFkJDBgwwNhErck+AsrYPrbumpWxdcbHzxz3GkhaPIDtB/fAkXBEGpFA0GcoEHWYahSlRxrGBW85JZMCcM6cOUbcPPHklS1blp49eyKiTLxrkuSUsOzd83Yi2NseQCnz3nvvGadyd+3aZSwFy2EN8Qpe7f6S8eh3TgJQsonI7NOnT3ocQOnfxIkT3TH9OHDgAIMGDTLCwMihFvH8ifCTvYyeKb97ACUItUT3lcMXgU4yFeXghwjAz2XOAKmykuK6zUP2+cnalFz1Jht6hrlOAUukyayngHPqmx4CCYDljiYdZcGWBcS/9x77FjwG22/innvgzTedlR8/eZzDKYe5sOSFId0gHYChhm0V6jmx3zTKWBnbT8DeFrzcBVyTAvxGk6gCxtVvOV0tYVIA2juC8Kk9vwpAWXeTiImyw0hEmWwNl+MxsuMm8w5P/20ldcsGBvfxJJmW8t+ygcB9p+9IQHYPizj80UscQDOtqwfQDKUc8shemzrT6hj7f6IpSIN94xnc6gmuv94Z6V6TElACSkAJhJ6AFw+ge+31amL4htZRUVyTw4lgFYDqAQRkF6SEjJRAzJ+67tqVbcXiqcs4Mx/6+W6mB+oBNEPJI49s8neH13C//OoGCbMInet1PidchY/Va3YloASUgBKwiYAXD6C7pe5EM8+4L0vOsnnzBKoAVAEoS+iuGzoOA/LnDlAlGwsyAmzZNIEDXK16AE0CXfvrXzw29zV+eyWenzcU4OJzQ0xlW5PGATQJ2UI2ZWwBnsmiytgkKAvZlLEFeK6if/3lPHEuEUskooE75eABdGfpRwwvUi0qitvSnGt7nimzAPz6nnvuaS2hUfJrknAzrtPLEpFkSyg45LRib1d/9gMSikXCx7oFoATskBs5zB4ot6tvvtarHsBciJ1IPsHgz0cy48dpUOAMF61azPxnbsYV5N0Ub72lwhQmS5mUsSV8pgorY1OYLGVSxpbwGYW3/36a2vdMpf3ltfh0YpdMFWZzF7BnnlrERH9HQcf53CLxjiT2lOtxZgE4p2nTpv/56aefrHc4QmuQsV922WVyI4hEUsyIah/E8YRCAEqwrC8BCUzmFoAPAxKOvXMQxx6IplQA5kIxOTWZetPrcTrxPPrUGUN8t/bEeLsqIhDW0DqUgBJQAkrAEgGJEVns2cJUPtaFvZM/8FUASn7RFbOJoZfx3xdFQx0HnADWkEaKcV64fnR09Ga5UUNEUH5Lsvf93nvv5cMPP/zr9OnT4hALSQqFAKzrOoixTWKdAt+6bhKUmznktUhKugRswlq7j+424vN5XvllophmUQJKQAkogRAQKDKsEuULVmP3iO/TWxfvqtyc4QpmLGFgMgLwee+jCD1xId5HweiriCKGs2zD4TCC78XExHxSuHDhjn379qVt27bGtW/u2zhCMOSgNCk3kOzdu9e4Gu6bb75JS0pKEsfX8qA07qWRUAhA6YaEXr/HdRhEloTlxK7sDYy0pB5AD4s99/4K+na8mmJFMl+KbtWocnWexFLSZB8BZWwfW3fNylgZ208gMC10fK8jBaILsKD7gkwVmlgC9rUDTxcpUqR3TExMRYfDEZ3XBaBEI4mOjj574sSJNWlpaRKObp2vwAKZP1QCMJBjCGVdKgCBVZv20HX24+wv8xE3F5zE4niJpx24JPcbjx8/PnAVak3nEFDG9k8KZayM7Sdgbws2CEB7O6y150ggWAJQvH1m0lwzmcIojy4BSyDFx1tyotRaqh7vwse9p9CilsdloWFkLO2KElACSkAJ+EfAjyVg/xrSUkEjECwBmHV5V0L8yv4Ad+gXByBLwRcGbeSBaUg9gMBry1bz77GjPNU141L6wODVWpSAElACSiBcCDw952PG9DJOBpvZAxgu3dZ+ZEMgWALQs/lHgabAY65NpDKRJgIbgZcizFIqACPMYNpdJaAElIAS8I9ADoGg/atQS4WUQCgEoHgDJTrQaY+RS8hICYQYaWuH+UoAfpCwhMSz+3iw+YNBnbSzZs2id2+5tU+TXQSUsV1kM+pVxsrYfgL2tqAC0F6+wa49FALwINAI2OcxWAkA/YvrdHCwGVhpL1/sAdy+5wjtJj/GH6XmUrFYRXb130WRghK7OzhpzZo1tGzZMjiN5dNWlLH9hlfGyth+AoFpYWXC79zzxtPc3qgLz/e6zajUxE0ggWlcawkagVAIwNdcV8ENB34HqgEjXQLw/qCNPDAN5XkPYEICXHv39yR2vpoKx25ixROvUK9y5cDQ01qUgBJQAkog7Ah8vm4r7T+ry5Vn41n1zNj0/qkHMOxMZalDoRCAstz7AnAXUNh1BcrbwOM4Y4VHUsrzAjApCa6/Hm68fx3x97YgOjoUUyaSpoT2VQkoASUQ2QR2/HWIWq+Wo86JB9ny/CsqACPbnNn2PpTf5tJ2eeBfCY4YoXzzvAAUu6SlQVQIZ8qWLVuoW1cukNFkFwFlbBfZjHqVsTK2n0BgWjib6qDqoFtodUEc7z/xiArAwGANu1pC+LUediz86VCeEoDHzhxj0JeDeOKKJ6hVtpY/PGwpM2HCBAYPHmxL3Vqpk4Aytn8mKGNlbD8Be1vQMDD28g127aEQgHIKODuPX0TGAZS7DDt27EhcXFyw7Rew9l5c+B3P77ibvSf/YOCVA3nuhucCVrdWpASUgBJQApFNQANBR7b9vPU+FALw3iwdkRMFEldkFvBshCGOeA/gmTMwdChM+u0Ropu/xksdJ9KnRZ88fyl3hM0z7a4SUAJKIOQE9Cq4kJsgoB0IhQD0NgAJDC1HjToEdHT2VxbxAlCu2I2Ph0uvPMWYl/bQ7pI69lPTFpSAElACSiDiCKgAjDiT5djhcBGA0o+jgAiqcEkvApe5OjMfmOKlYxEvAE+fhldfhT59oECBcEGfuR+6d8p+uyhjZWw/Aftb0HlsL2MVgPbyDXbtoRCAcgewZ5KwMHLNg/yFy8mDJq5r6a4BYoDfgEtdItWz7xEnAP8+/jcz1s1gdJvRREdlNUWwp5+59vT0pDlOVnIpYyv0zJVVxuY4WcmljK3Qy1x21o+zeGXDKyzsvpCqJZ2XdKkADBzfcKgpFALQ4eUQyHFA9gYuCgcorhtJ5gHtAbnyYpVLACZl6V9ECcCPfv2Ih5Y8xOHTh/nizi+Iqxm5h1bCZJ5oN5SAElACeZLAnbPG8u7+YXzUbgOdr2imAjAPWjkUArB1Fo4i/raFYRBoCVYtd+AUcu1PnObF/hEhABMT4dEBSayoX4/TjuPM7jib2+o5r/fRpASUgBJQAkogK4E7Js1g/om+TGz0P57o3NZ4rDeB5K15EgoBeJXLo5aV5JXA9wHC2x3o67pyrjhQEBDPo2caBTzg2ne43pV/syvDDa6bSeRQitxW8h3QBfgzSx1hLwC//Rbuugt274b7h/7EM4MrUum8SgHCHJxq9A5V+zkrY2VsPwH7W9B5HDjGo19fy4h5HzL1P/fzaE9nIH4VgIHjGw41hUIAHsvmsMdhoEyAoIiAk7qKAq96EYADgX7ATcBOYARwD1AbOAW0cy1J3wkIo5XAw8CmSBOA/frBnDnwwgvwwAOhvdHDX9vOmjWL3r1li6gmuwgoY7vIZtSrjJWx/QQC18L69fDii84Dgi1bOutVARg4vuFQUygEoCz5npdl8OJl+9u19y6QXGS5ebkXAbgLmAy4l3XloMc+l9fvHZfokwsQ67sOgUgd8V46FpYewN1Hd/PLP7/QoXYH5JTvnj1QW6StJiWgBJSAElACfhJQAegnuDAtFkwB+K3r8McVwOosPOQGkK1AoE8leBOAItoSAenHWo9+LAV+AZ70wVaGAGzXrh0NGjQwirVp04YVK1Ywbtw4YmNjjdfmz59PqVKl0m8K2b9/P1OnTmW8BOFzJfEONGnShJaun1pymm3RokWZrkCTEAedOnVKvxdXljsSEhIyecdufeBWlldYjqOYgz8e+4P1364nMTGR7t1lVRySkpKIj48PWf+GDBlC//79qVixotEfiS6v/XNOAjP2VX5OV4S/7w/lp/z088+5nOvt+8Pb+2PlypXG94akqV/NwfGLfH1SUg4F+/BdqVnDkEAwBaAss0oST9o4DxayN2+/6CQbJpQ3AVgF2A3Uc4lOd1fk1K9M6Id8sFPYeADT0uCXAxtp8nITzi92PnM6zaF9LTnErEkJKAEloASUgHUC6gG0zjCcagimAHSPuyfwbpAgBMUDGOq7gFetcu7T+PhjWHF0NrfUuYUKxSsECbH9zcivUk9vqf0t5r8WlLH9NlfGyth+Ava1IKs1d8UP4uCGjdKIegDtQx20moMlAKWdNNeocoo+nPWkrlUQvuwBlD2IAwDZA2g2GR7AUAlA8fqNHu38K1wY3nkHbsuD0V1kydy9ZGzWMJrPNwLK2Dde/uRWxv5Q862MMvaNV265q0+tTosLWvB+t/eN7TqLFy9m+vTpKgBzAxchz4MlAD1P/noLBO3GJYcxApFEZEroFxGAn7sOnaQCyS4hKvv85BSwhHmRAyHDXKeA5SJcOQVsNoVsCXjP0T3ERMcw9qkLEA/ge+9BPVnU1qQElIASUAJKIAAEKkysQK0ytfiul0RC05tAAoA0rKoIlgC82hVLTwafNRC0JxAJtxKIJLeKzPHwOro9kG2Ab1wNjHRdPycnkn/MEgfQbB9C4gH88NcPeWjxQzSr1IzFXb8iJibK8ABqUgJKQAkoASUQKAIln65FjKMIh8duVA9goKCGUT3BEoBhNOSAdiXoHsDXf3qd+z+5n3JFy/H6La/TsU7HgA4oHCuTpYe4uEAfEA/HkYauT8rYfvbKWBnbTyCwLVx8y9ucOl6QAyucUST0LuDA8g11bcESgM57ZHJPEm8vklLQPIBynVupUnDszDEe/+JxxrQdE3E3evhrWAmj4w5j428dWi5nAsrY/hmijJWx/QQC20Lr1vDrr/Dvv86QXboHMLB8Q11bsASgmcMdckgkUHsAg8XVdg9gaiqMGweTJ8OPP0KNGsEamrajBJSAElAC+ZnApEnOa0SnTnVSUA9g3poNwRKAeYtaxmhsFYB/Jx6ix61l+GZlFI0awfvvQ11nDE9NSkAJKAEloASCSqDJgIfZOGWmtKlhYIJK3p7GVABa42rbEvCSbUu4b9F9NNg/jsYpD/Lcc+C6WMRajyOwtEShd9+qEoHdj4guK2P7zaSMlbH9BOxrQZaA2w/ugSPhiApA+zAHteZQCEBZ5h0E9AKqAnuB1+UmLEBCtURSssUDOPLrkYxaOYrSsaWZ3fE1utTPg8H9fLDygAEDmCxr4JpsI6CMbUObXrEyVsb2EwhsC2lpaZxMOUmRAkWMsGN6E0hg+Ya6tlAIwOeALsDzwO9AdeAJ4GOXMAw1E1/at0UAfrHjCyasmsDcW+dStaRo5Pyd1HNiv/2VsTK2n4D9Leg8DizjEStGMPqb0Wx+eDP1y9dXARhYvCGvLRQCUG7cuBbY5jF6CcD8NVAp5ER860BAloDlRo8tWzIHcpZfXlFRoTCPbwA0txJQAkpACeRNAs+vep5BXw1izf1rSNycqEvAeczMoVAY/wBVgBQPloVcS8HnRxhfyx7Affsd3N8rmm+/hZ9/1lO+EWZ/7a4SUAJKIM8SePnHl+nzaR++uvsrrqt+nXoA85ilQyEAx7muZBvluqlD+jAcEBE4NML4WhKAC35ezu2v9+fsm59xX5eqvPgiFC8eYQSC0F2Nn2Y/ZGWsjO0nYH8LOo8Dy3jeF38w/rXNjH7ocjrdUE4FYGDxhry2UAjAb4HLgUPAbuBCoAywNguNa0JOJ/cO+CUAHWkOhv5vqLHPr0BaMQZU/YBn778x99byaQ69QcF+wytjZWw/Aftb0HkcWMaLFsGtt8LcuXD33RoHMLB0Q19bKATgCJPDFg9huCe/9wDeu/Befv33V97r8h41y9QM93Fq/5SAElACSiCfEVi+HK67DqZPl+1JehNIXjN/KARgXmLolwdQAJxMPknBmIIUipGVb01KQAkoASWgBMKLwI4dGDFo77gD2rZVD2B4Wcd6b0IpAM8D5M8zyQnhSEqmBOCRI7BxI8i9ipp8J7B//34qVqzoe0EtYZqAMjaNyu+MythvdKYLKmPTqPzKqFfB+YUtbAuFQgC2BN4AanlQkX7kybuAX16yhoFz34EvXmTb1igqRVqgmzCYukOGDGH8+PFh0JO82wVlbL9tlbEytp9AYFtIdaSy88hOihYsSpUSVfQu4MDiDXltoRCAm4CvgLmyEpqFwNaQE/GtA9l6AOWN8+TCZ3khYQQ4CjCq8nqe7t0ADe3nG2DNrQSUgBJQAqEhkJiUSOkJpenVtBevdXpNBWBozGBbq6EQgMdcF0mLxy/SU7aHQA6dOkSjmY1wnCzD7Lh5dLy8YaSPVfuvBJSAElAC+YhASmoKhcYUokfDHtxb/F5GTJ7G2mVLhEBJQL7LNUUwgVAIwC+AAcCvEczN3fUc9wBuPrCZ6qWrU6RgkTwwVB2CElACSkAJ5DcCBZ8pSIdaHVh4x0KNA5jHjB8KASi3gLwJLAX2Z+Epy8KRlEwdAomkAYVjX2fNmkXv3r3DsWt5pk/K2H5TKmNlbD+BwLew8Z+NlI4tbdxLX7hvXZJnGDu11AMYeNRBrzEUArA/MBk4mmUPoCwJS1DocEiXAJM8DqZcAZR39dmzf4YAfGXRt+w+eoxn7m4fDn3Pc31Ys2YNLVvK2SFNdhFQxnaRzahXGStj+wkEvoUePaBUKZg5ExWAgccb0hpDIQAPAA8An4R05OYbvx54BOjkpYghABlYkKioUhyI/4NyJYuar1lzKgEloASUgBIIYwI1akCJEvDTTyoAw9hMfnUtFALwoMubFimHQN4GPgIWZCcAC/Srzptd5tPz2hZ+GUELKQEloASUgBIIRwKNGkFyMmzdqgIwHO1jpU+hEIAzgC+zEVRWxmJHWQlULWFragBnsxOA2//cS80LK9vRvtYJbNmyhbp16yoLGwkoYxvhuqpWxsrYfgKBb2HsWKcAHDVKBWDg6Ya2xlAIwHmu5dTvgH1Zhn9PgHB0B/oCTYDiQEHAkaVuuWtYlqJlGXe9K//mLHnkeSNA9i16S3oIJEAGy6maCRMmMHjw4CC0lH+bUMb2214ZK2P7CQS+hZ/2/cSBkweIqxnH03M+ZkyvLtKIHgIJPOqg1xgKATgnm1HKknCvABG4ASgDyIa8V70IwIFAP+AmYCcwAhDxWRs45dEHEamy/+8nFYABsoxWowSUgBJQAhFDoMO7HVj++3JODz2tgaAjxmrmOhoKAZi1Z+Kh+w/wXyDQ0ZLl9t3lXgTgLtdJ5GmuzsS4vJGPA++4XpNl349dXsTsaKoH0Nw801xKQAkoASUQgQS6fdCND3/9kNThqZw4foKSJcX5px7ACDTlOV0OpQCUpdWHgTtdvRJPnQSIDmTyJgBFtCUCEtplrUdjEpfwF+BJHzpgCMB27drRoEEDo1ibNm1YsWIF48aNIzY21nht/vz5lCpViri4OOPfcmH51KlTM91vKzHCmjRpkh7uRPYLLVq0KNPSpywhderUKX0/nISVSEhIyBQjT+4b7d+/PxUrVjTaWrp0KYmJiXTvLqvikJSURHx8vPYPUH4Z9yvr/NP3h36+6Oezt++PJ2c/yaa/NtHv8n6kJKUgnxUqAH1QCWGcNdgCUPbidXMJPxFgsvwqpydkr94OGzh5E4ASiHo3UA/wvHtY9ibK1TYP+dAP9QD6AMvfrLp3yl9y5sspY/Os/M2pjP0lZ76cMjbPymzOE8kniCKKogWLcvz4cfUAmgUXAfmCKQDHufb4lXKdAH4F+Br42yUAJT5gmuMPggAAIABJREFUoFNQPIB9+/alY8eO6R6+QA8iv9enpyftnwHKWBnbT8D+FnQe28dYVpMWL17M9OnT1QNoH+ag1hxMASincCUG4F3AMo9Ryklg8QAGSwBK0972AIoQlSVo9x5AM4YwPIAqAM2g0jxKQAkoASUQaQTefBNeew3uv38p69apAIw0++XU32AKwPuAPkBz4GdANhK8C2wDmgZYAEa7Dn6IB/BzQOL5pQLJruvdZJ+fnALu4BKDw1yngOtkOQWcm611CTg3QvpcCSgBJaAEIpbA+PEQHw/r18MVr+ldwBFrSC8dD6YAdDcv12VIjL7bXbH5CgEStuWbAIK9F5BwM+7bRmSc8t9tPNoZCfR2icMfs4kDmFuX1AOYG6EAPNc7VAMAMZcqlLEytp+A/S3oPA4846lT4bHHYOLEpQx6qweOhCPSiMYBDDzqoNcYCgHoHqTsBZS4fyLCarpuB7kx6ASsNageQGv8TJWWU2e9e8s00WQXAWVsF9mMepWxMrafQOBbeOfzHbz6v+UMvzOOG1+NI3mGcXZSBWDgUQe9xlAKQM/BtnPFAewcdALWGlQPoDV+WloJKAEloATCmMA7G9/hrgV3MbzycMZMe0k9gGFsK1+7Fi4C0Nd+h0t+9QCGiyW0H0pACSgBJRBwAgt+W0Dn9zvzTud3uO/l0eoBDDjh0FWoAtAaexWA1vhpaSWgBJSAEghjAl/u/JKuH3RlRvsZ9Jr1jArAMLaVr11TAegrscz5dQnYGj9TpSU6/Xg5iqbJNgLK2Da06RUrY2VsPwH7WpA4gHfFD+Lgho3SiO4BtA910GpWAWgNtXoArfEzVVquznNfbWeqgGbymYAy9hmZzwWUsc/IfC6gjH1G5lOBY8eO6U0gPhEL78wqAK3ZRwWgNX5aWgkoASWgBIJB4OhRKF4cYmIyWktLgyjzMkAFYDAMFbw2zFs+eH2KpJZ0CTiSrKV9VQJKQAnkNwJbt4IE85MrPdauhYYNMwiMGAHbtsGUKVCxolcyf/8NPXpAo0ZLAb0JJC9NHxWA1qypHkBr/EyVlr0ncXFxpvJqJv8IKGP/uPlSShn7Qsu/vMo4C7czZ+Dii+GGG5zRnC+5JHOGvXudr8s1H6tWwQUXnANeslSt6sw2apQuAfs3M8OzlApAa3ZRAWiNn6nS8+fPp3v37qbyaib/CChj/7j5UkoZ+0LLv7z5nnFqauYlXsF4+jQUKZI9UFkGHjIEevaExo3Pybdn/2ku7PEscU2b8v6o63QPoH9TMyxLqQC0ZhYVgNb4aWkloASUgBKwSmDPHpg5E956CzZsgPLlrdaYXv5A4nEqTC1B6xL/4ZP7p6oADBjZ0FekAtCaDXQPoDV+WloJKAEloASsEPjnH7joIpBtMo8+Cm3b+nSwI7emU1JTKDSmENeevZZNnx3WMDC5AYug5yoArRlLPYDW+JkqnZSURGxsrKm8msk/AsrYP26+lFLGvtDyL2+eZ/zvv7B7NzRvnhnQvn1QqZJ/0HIplZaWRuOXGxNXI46XFi3RQNC2UA5NpSoArXFXAWiNn6nSAwYMYPLkyabyaib/CChj/7j5UkoZ+0LLv7x5lvGHH8JzzzkPa3TpAu+/7x+g3Ep99pnz0Ej9+l5zFu5bVwVgbgwj6LkKQGvGUgFojZ+p0nn+V70pCvZmUsb28pXalbEyzpZASgps3AibNzv/nngCzj8/I/uaNbBjB1x/fbbhWgJC95574MILYcwYFYABARrelagAtGYfFYDW+GlpJaAElED+ILB8OUhMlf37oWVLuOaajHHL640aOWP0NWgAQ4c6Y68EO732mvMgyddfqwAMNvsQtKcC0Bp0FYDW+GlpJaAElEAGAYcDxBsmoUnkT/b+et5UkZTkDGsiz6KjoVSpzPQOHxZXK0g9RYtCmTKZn//+O5w96ywvHjbP8qdOOQMlJyc7/1q3hhLyEe9Kv/4KS5eC5JM6JIiyZ3r2WfjgAzhxAi67zCmkPFOnTiC3cVSoAPfeC+3bZzyV/kjy4VYOW6bNzz9D795ODl6SLgHbQj1klaoAtIZeTwFb42eqdMBie8kGatksffw4FCsGTZtmbv/jj0FO1MmHf4sWcNVVGc/ltWeecX5xyJeLLJXUrZvxXKLpT5/ufCZ/sl9H2nCnTz+FRYuc/5IvnqxLLC++CL/84nx+xRXQq1fmvkkUVvlikXTXXXDttRnPxaMwbFjGv+WLydN7IAFeX3/d+bxQIWe4CM80dy7zX3+d7tWrQ716MHBg5ufDh4OEmZDUsSN07pzxXALNPvRQxr8ffzwzV1nOGj8+4/nLLzuvo3KnBQucX5qS5CaCrHs95Uv1p5+cz8Vj0rdv5r7dd5/TnpKkH+3aZTyXKwz69cv4t9Qt+5vcSTwyckOCm4u7H+7nM2bA5587/yXx0caOzdy29EUEhaQ77nDOCXcSEXLrrc75Imn0aOb/8UdGPEsZ06BBGfk/+iiz2Hj77QybVakCc+dmbnvwYPj+e+dcu+mmzPaXnNddB4mJzvbj46Fr14zy0ucOHTJE1sKFmeey7Dd78skMASY3SXgmmR9yq4TUffXVIDb0TDJ/ZS7Lcxmjp1A6edI5TrfAW7LE2Rd3+uorZ9Bidzp0KLOIk/fN0087n9auDVn6Nr9ZM7qLiJEk9pB+eiaZY0eOOMXjSy/BAw9kPN21yzl/5D0if7LPTtpwJ5kv8h6X97UIR3nPeqZ16+Cvv5zzW9rxvHEjc86I+5fo1d27l/LNkR44Eo5I/0sCxyJuINrhTARUAFqbEOoBtMbPVGmv0f3dXyDyQe5O8qtVvtBlOUVEgVs4uJ8/8gjIF+t55znFhPy3Z3r4YafQkQ//bt2cX+ruJEJHykt78tenj3PJxp22b4dXXnE+k1/x8iXlKQBlSeV//3PmLlvWGVbfM0lfRERKEmHqKbLkNRFC4vmQdPPNcOmlGaXlS1K+zNzpv//NvE9IvhDd4rNAAefykmf69FOWzp9PXM2aToHkKWQk36xZcOCAs4R84bdpk1FavDUidt3p9tuhVq2Mf//xB7z3Xsa/ZdyeQWlFxHz3nfN56dLw4IOZ+ybCyC2yZGlMxI5nEpEmtpEkosczkK180Xt6YeQ+K8/4aL/9lmET4SLcPNPKlSBeH0mVK8Mtt2R+LoJRPE6SmjVzen3cSTxEb7yR8e927Vj6228ZN9qIOBVvkjtJ3zxPuickgFvIiGC67bbMba9Y4fwxI3OtWjXnkqJnEmElfZDnTZpkFr4iwoS5PJM/KSvvCXeSHxQyn+WZ3Bsrgs4ziU0PHnTOdSnnaW/Jt3On80eUlJe57slcBKucYnW3Lc8854PYUoSr+30kHjzP97g8d9ctfcsS4HjpZ58Rd+ONGfVn7rn+y08Cle4Yzem0o+yePULjAPrJMByLqQD0bpV6wPOAxB6RbxePn6iZCqgADOasluUTERDypbxlC8iXXKtWGT3YtMn5mnxZyxeLxMUK9ZJKMPloW0pACSgBGwgUf7wlSQX+4vDTm1UA2sA3VFWqADyXfAFAfprLmonh684hqQC0Y+bKL/zVq501yz4cd5JrjkaPdoYokOVX+Stc2I4eaJ1KQAkoASXgIlBzzLXsO7uZfQN2qgDMQ7NCBeC5xrwakE1QZwHZQSybSDzWcjIVUAEY6DfDyJEwcaJzH414+556iv3791NR9tRoso2AMrYNbXrFylgZ208gAC3IdpV333VuaXGl+P/Fs/PITma3m60CMACIw6UKFYDnWqI7IFGHZYNXEvAN0APY7sVoKgCtzGTZoyT7eDyXaeUey4IFnRuoXa8PGTKE8Z4HCay0qWW9ElDG9k8MZayM7ScQgBbk0Fj37s691FnSsWPHVAAGAHG4VJEXBaAIODkq2ASQ44YFAUcW4KMAOf4lAm69K/9mVx45gtYHcB91nAjIeuRHKgADNG1lc/6rrzoPLsgpPc9DBQFqQqtRAkpACSgBPwhIJARZcZHDQhJKxyOpAPSDZxgXyYsCUAScLN3KzH3ViwCU5V2JDSHHCXcCEsxJ4jfIef9TgByH+xaQI32yDLwCeATYqAIwADNZTgHKaVMJU9K/v/NkpZzA1KQElIASUAKhJyARFiQMkISWkpP5HunpOR8zplcXeUXDwITeUpZ7kBcFoBuKnB5Y7kUA7nIt8U5zZYwB9gGPA++4XpOYCwMAebYY8Ahklom5LgH7MwUlHp9naAh/6tAySkAJKAElEFQCGgg6qLhtbyy/CUARbIkSahfwDHUup34lCu+TPhI3BGC7du1oIDHKkNXMNqxYsYJx48YR64rrJYGMS5UqlR4DTDaDT506NdO+tlmzZtGkSRNauuJ5bdmyhUWLFjFYAr660oQJE+jUqRN1XQGI16xZQ0JCAr0lcrsryT6j/v37px+akBh6iYmJ6QFo5T7S+Pj44PRvwQISJk2i97x5IMFsAX/69+CDD9K1a9f8xy+I9r3zzjt55JFH8tb8CyI/M+/fzp07G++7iHn/hhk/M59/N954I2+88UZ4fP5FID9vn88rV6407rGWNPWrOTh+ka9Q9QD6qBXCMnt+E4CiQnYDEufPM7z9PFdUc48rDUzZSz2A3jDJEoJcZi43TnTpAhMmOGPz+ZlE6LqFsZ9VaLFcCChj+6eIMlbG9hOwp4V7przGwj/e4nTqXs5Ol51TKgDtIR3cWvObALTFA9i3b186duyYEeU/uDYMz9bkRo62bTPflhGePdVeKQEloASUgCcB+REvNw+5DoFcM2I43+56hqhfSpCWYNwAp3sA88CMyW8CUEzmbQ/g3649f+49gGZNq3cBu0nJB4beumF23mg+JaAElED4Epg2DZYtg08+Mfp4wzNj+WrbMBWA4Wsxv3qWFwWgXA4roV/kEIjc5C6nelOBZEBuZpd9fnIKWK53EzE4zHUKuI7rFLAvIHUJWO40HTDA6e17yNcVdHOoZT+ke9+UuRKay1cCythXYr7nV8a+M/O1hDL2lVg2+RcsgDFjYL1ESYPxHyxl3k+L+PXYUs5Ol69N9QAGiHRIq8mLAvBeYI5L7AlcGaMIP7nBXoI6SxoJyMkJEYc/ZokD6ItB8rcHcO5c6NcPOnZ03t5RqZIv7Eznlc3fnodhTBfUjKYJKGPTqPzOqIz9Rme6oDI2jSrnjGvXwq23wj4JkOFMcqBwxORprF22RAVggDCHupq8KACDyTR/ewDlyiCJ6yfeP01KQAkoASWQNwgcOwY//ADXX59pPBoIOm+Y1z0KFYDW7Jm/PYDW2GnpEBGQkA7JybIjQlMoCRQqVCg9VFQo+6FtKwEzBMQDuHjxYqZPn64eQDPAIiCPCkBrRso/HsDVq0GWfGfM0MMe1uZMSEuL+KtWrRoSi1JTaAlUrFiR33//XUVgaM2grftAQD2APsCKgKwqAK0ZKe8LQAkF8PTTzph+I0Y44/vFyAUpwUu6rydwrN0f4Hv27KFECZm+mkJBQOxQtWpVjh49qnYIoAH0syKAML1UpQLQXr7Brl0FoDXieX8J+OhR6NXLeSKsnsTPDn7Sk32BY+7+AFfhETim/tSkdvCHWu5l9LMid0amc8gewAsuMG5xWrN3DZ3HdabB4QZ8Nf8rqULjAJoGGb4ZVQBas03e9wBa46Olw4yACo/wMIjaITzsoL3IgUBcHHTtCg8+yJsrvuM/37TijuITmPekcT2pCsA8MHlUAFozYt4TgPKrT65ts3B1mzWkWtpOAuEqPC6++GKKFi3Kpk2biI6WUJ5w6aWXMmnSJK655hpGjRplbD6vUqWKcS9p7dq1mT17NuXLlzfyyh3cjz/+OLfcckuO+F5++WVmzpxptCH1tGjRgrfeeotmzZoRFRXFmTNn2Lp1K40bNyYtLc2IP/nee+8xevRoPvjgAwoUKGDkkZt/ZLkxa9q3bx/33Xcff/75J4ULF6ZWrVpIm2XLls2UNVztYOfc07ojjMBddzlXfYYO5c0vf+Q/318KR8rB1IMqACPMlNl1VwWgNUPmnSXgs2dh7Fjnvb1vvgndulkjE8DSeodq4GCGq/CQgykiyiTe40OugOJZBaAsW0+ePNmA0b17d0P8TZMbC0wKwPXr13P77bezYcMGSpYUBwb8/PPPNG3aNB2wCDcRg4cPH05/7aOPPjLaXbFiBXJy1+FwsHnzZho1anSOYQ4cOMCOHTu48sorjWeDBg3i33//Zc4cCU2akcLVDoGbaaGpST8rAsj98cdBbnh64QVWrNtP215PQczHkHBcBWAAMYeyKhWA1ujnHQ/g6NEwfz68/TY0a2aNSoBLz5o1i969JW63JqsEwlV4iAAUL1t8fDzbt283TsbmJADFq/bpp58aYSnMCsBFixYxcuRI1q1bZ3jyvCVvAnDq1Kl89dVX6W35YgMRj+K5XL58uQpAX8D5mVc/K/wE561YUhIULGgc+vvnHxg3DmYk1+Xsy1tVAAYQcyirUgFojX7eEYDHj4N8KRYpYo2Ilg5rArkJQFlB3bnT2hBq1Ei/QtR0RSIARaDJkq8suw4ZMiRbAShLsF26dKFHjx7ceeedpgXg6dOnadeuHbt27aJVq1bGn5QvVapUej+9CcB//vmHtm3bGkvGV199Nddee63RtojUnJJ4Cq+//npuvfVWHn30URWApmeDZgxXAoX71iV5hgrAcLWPr/1SAegrscz5I1cApqYGPZyLNdRaOhAEwl0AytLs5Zdfzq+//kpcXJzXPYCyxCp797744gtjSVaS2T2AkjchIYFvv/2WBQsWGMu18m+3CPQmAKWM7Adcu3Ytq1atYt68ecYysPw7O0+ilBGv9cGDBxEvYNaUmx0CYWutQwkEmoAKwEATDW19KgCt8Y/MPYCybDZwIHz3HZQrZ42Alo4oAuEqPNweQDl8MWDAAIPpd999x8SJE9MPgbj3ACYmJhqeNfl79tlnfRaAngarX78+48aNM7x0krITgJ5l5BaVChUqGHsC5XYEEYRygEQOhdxwww1GVvH4SZDnhQsXEuMlbma42iGiJrN2NqgEZK63H9wDR8IRaVdPAQeVvj2NqQC0xjWyPIDi9evTB95/H2TzvCyfRYX/FJDlwPHjx1uzlJY2CISr8PAUgIcOHaJevXqGd+/dd989RwDKOOTwhizH7ty50xBjZjyAcrpXxJv78Mbu3bsNT6Ls7xPh6RaAcijkyBHjS85IcnikdOnSVK9e3fi3eAyl7T/++OOc071u8Sf9EvFXUPZQeUnhaodIf5voZ0UALfj33/DqqzB8eHqlTQY8zMYpM1UABhBzKKsK/2//UNLJve3IEoAynuefh9tvh4suyn10YZJDri2Ta7M0WScQrsJDxJUIJrcQGzNmDCNGjDC8bO4wMJ6ngIXE/fffT/HixZFDGiIA5WSu7MuT5VrxyEnYFllOdqeffvrJCBUje/rc+/fEUydhW9zJmwdw2bJlRl+kfSkny77Dhw/n5ptvPscg33//vbG3UPYxupenZWxZl4HD1Q7WZ1hoa9DPigDy37YNmjSB06eN91Tx8cW5ofINLPrPIhWAAcQcyqpUAFqjH3kC0Np4tXSEE1DhER4GVDuEhx20FzkQOHAAKlQwBCCxsRQdW5RWFVqx7MFlKgDzyMRRAWjNkOEtAP/8E7Zvh+uvtzZKLZ1nCKjwCA9Tqh3Cww7aixwIJCfDFVfAl19CmTIUfLo05U81Yd/klSoA88jEUQFozZDhKwDnzYP//te4xsdY9o3gJJuP5USoJusEVHhYZxiIGtQOgaB4bh36WWEPV6m11NXzqVi8CFuXdlIBaB/moNasAtAa7vA8BSzXud10k3MD7223WRthGJSeP3++cfODJusEVHhYZxiIGtQOgaB4bh36WWEPV6m1YsWlpKYu5uDB6SoA7cMc1JpVAFrDHb4ewMRE8Ahwa22YWjqvEFDhER6WVDuEhx20F+YJjBkjQSOOMWyYcY2ihoExjy5sc6oAzN40ScD3rsdfAt7ikISHAHQ4IDo6bCeZdix8CKjwCA9bqB3Cww7aC98INBrQm01TXlEB6Bu2sM2tAjB70+wGLszFcqEXgHv2wN13wyOPQJcuYTvRrHRMruDK7dotK/Xnp7Ju4bFnzx5KlJDpqykUBMQOVatWNULLqB0CZwH9rAgcS6Om+Hjn3fDduvFn4p9cPKEaPJumAjDAmENVnQrA7MmfAH4ATsnbANjoJWtoBeDChSAxzLp2hRdegGLFQjWPbG1XboaYPHmyrW3kl8rlC1KCLku8NE2hJSCxLeW2EP1xEzg76GdF4FgaNd11F9SvbwhBiQVYuF89UvQu4ABDDl11KgCzZ18WOAQ0B94G6oWdAJQr3eSofh71/Ll5L168mI4dO4buXZLHWhYRKDdieKbPP/+cm+TgkCbbCGRlLIGiVfwFFrd+VgSWJw8/DCVLgusmJr0LOMB8Q1xdXhWAcmS0L9AEKA7IfUyOLKxHAQ8A4sVb78q/ORt7bADkkk8RhJ4ptB7AEE+eYDWvv+rtJ62MlbH9BOxvQedxgBkvWOC8LtR1V7YKwADzDXF1eVUAilgrAxQFXvUiAAcC/QBxeewERgD3ALVdS76ylnraJRrlzrT/ATVD6gFMTWXp6NHEjRzp1/29vsbHyi1/Ts+ze+btdTOvBfNDPbdx5/R+9bVsbvmV8bm0c2PmWcJMXmWsjP39DjYzv7Kr25eyZvIGax6rAPR3toRnubwqAN20WwPLvQjAXYBsKpvmyhgD7AMeB94BLgVmAceAAq49gN9kJwBt31Avl3I/9BDxGzcy7uefjajsvqb4+HjGjRtnulhu+XN6nt0zb6+bea1bt27Gva7BSLmNO6c++Fo2t/zK+FzauTHzLGEmrzJWxv5+rpiZX9nV7UtZM3mDNY/LPdGclFd3yLA0DIy/EyeMyuVHASjLtonAFcBaD1ssBX4BnvTBPpWBvT7k16xKQAkoASWgBCKdQBXgr0gfRH7vf34UgDJxJcSLHOrY6jEB5rk8fg/5MCmE3wXAcR/KaFYloASUgBJQApFK4Dzgb8CIB6MpcgnkRwEYSA9g5Fpee64ElIASUAJKQAnkWwL5UQCKsb3tAZRfNANcewDNTgj1AJolpfmUgBJQAkogLxBQD2BesCKQVwWg3IsmoV/kEMjngEzYVECCn4nbWvb5ySngDi4xOMx1CriO6xSwWfPqHkCzpDSfElACSkAJ5BUCugcwD1gyrwrAe4E5HnsUZJwi/NoA7tO8I4HeLnH4Yy5xALMztREHcOLbn7By52+M6NqNWlUkfnRG+vHvHxn65Ugcey+jdY3/t3cmcDKXfxz/sLvu+4oQOVeuKEWJUjlCmzsqSo6QZKu/Wyo5citFlFKSlBzJTSgslTNHzpzrtsuuvff/+v52Zs2u2fn9Zn7zm5nfzOd5vbww85zv55n5feb7PN/vUw/dH38Kd5cUbep80uINllWtzpZVy+8przMZT58+ffDZZ585D8yFEmrjdlSls2XV8pPxnbTVmNmW0JKXjMnYha8JpYiW9eWO72Mt7XhqHVuvMKQXsKurxrfK+asA9BRl1UDQC/cvRLefuyM+JRZIyAOMi8K5M8EoVcr5LmqJB5VVrc6WVcvvqbhTMh7GAbxzVhlrUf3zo7aGpQauY+1ryx5xMm6mvhDt5NDCzVpMS15PrWPrXeIUgC5Nu88VogDUNyWqAlCqj0tIwtJt/+C33SeQ4/hzmDYtY6OD1w3GllNbUDl3fZSNb4beTzVFGTGwMykEPCkAAxU5GRs/82RMxsYTMLYFBoI2lq+na6cA1EdckwBUa2LgqoGYu3suouKjgD0vAT/PQ926wF9yQZ0lJaUkITi7xKQOvDR16lS8+eabgTdwD46YjI2HTcZkbDwBY1ugADSWr6drpwDUR9wtAlC6kJKagvW7D2Pnzmz4769Q5M0LTJa7SuTwYmoqyk8rj+rFq6Nn3Z5oVaUVQoJcO0eob7jeKT1r1iz07i3HNZmMIkDGRpG9XS8Zk7HxBIxtgQLQWL6erp0CUB9xtwlAR92ISYhB+0XtsfbYWiSnJqPp1Z9RJ/dz6NEDqGTvhmJ9Y2JpEiABEiABEriDAAWgfy0KCkB98+kRAWjt4rkb5/D1rvmY2+tNHDmcZgGcPRuKEGQiARJwH4G4uDgkJEjUKCYSCCwCOXLkQK5cuewOmgLQv9YCBaC++fSoALR2NSUF2LgxTfyNGweULw9M2TYFJ6+fRI+6PVDzrpr6RsXSJBDABET83XvvvYiMjAxgChx6oBIoWbIkTpw4YVcEUgD616qgANQ3n14RgPa63Pzb5lh9bLXy1tsN3sawhyYgOBjIl0/fAH2h9JAhQzB27Fhf6Irf9oGMb0+tNdTF6dOnUaCAfMSZSCAwCFjj/EVFRdld+7XD+2LvFCUma0EA0YFBxX9HSQGob259RgCKo8i2M9sw5+856Fi9I7bOa44pU4DOnYGePYF69fQN1JulxRIjv0qZjCNAxncKwKwegsbNAmsmAe8SsP74yWrtMw6gd+fH3a1TAOoj6jMCMPMwvv8emDgxLZRMmzbA7G+vIOJsBJ6p/Iy+EbM0Cfg5AbWHoJ8Pn8MLYAJqa58C0L8WBwWgvvn0WQFoHdauXUBQEDDn7Bv4eMfHigCc0mwKqhStom/kLE0CfkpA7SHop8PmsEgAamufAtC/FgkFoL759HkBaB3e2eizGLRuEObvm4+21drip44/6Ru5B0truQrJg93xy6bI+Pa0qj0E/XIBcFAkIIf6oqNRsGBBcAs4MJYDBaC+eTaNALQO849Tf6BU/lKoULgCfv0VuHQJeOklIHt2fSCMLL1w4UJ06tTJyCYCvm4y9m8B+PXXX+OVV15BoUKFFA9PechbU3JyMkJCQjBq1CiMHDlSeXnTpk144okn0vMEBQXh7rvvRqtWrfDhhx8q9VjTf//9p3hNz5kzB927d3fqs7R7926MHj0aO3fuxIXSWHOdAAAgAElEQVQLF1CkSBGEhoaiTZs26N+/v9Jn6bujlC1bNmzcuBGNGjWCWn2O6tmwYQPmzp2Lbdu24dy5c8p4mzZtivfeew/Fixd3alxmzUwBaNaZc63fFICucbOWMp0AtHY8NRV4+GFg5860v6d/nIKH6vmwCtQ3TyxNApoJqD0ENVfkQxmtAlDE0qBBgzBmzBhVAdikSRN8/PHHePDBBxEbG4v169dj3LhxeOaZZ7B06VLdAlBEn4i2+vXro0+fPoqj15kzZ/D777/j8OHDSnsiVi/Jr1RLmj17Nr788kv88ccfyG7zq/W+++7DwYMHVetzNCUdO3ZETEwM5O8KFSrgyJEjiiCWmHh79+5Fnjx5fGhGjemK2trnFrAx3L1VKwWgPvKmFYAy7OvXgVGjgE8+ARqMewWVKgNjnxyLkvnocatvWbC0mQmoPQTNODarAGzWrJkisI4dO4YSJUooQ3FkAVy3bh1ECFpTr1698MUXX+D8+fPp5V21AHbt2hVS/6lTpxAsMas0JLHGvf/++0hMTMwgAKWoK/XZNnnlyhUULVo0Qy+2bNmCxo0bK6Lz5Zdf1tBDc2dRW/sMA2Pu+c3cewpAffNpagFoHfo/hxLx5s6WWHd8LfLnyI8F7RagZZWW+si4sbQE5s0qMr0bmwnoqsj49vSrPQTNuFBEAMr27ObNmyEi8NVXX8W0adOcFoAzZ85Ev379EBERoVgGJbkqAFu3bo3jx4/jn3/+0YzUkQB0pT61huVzIZa/Dz74AMOGDVPLbvr31dY+A0GbfoozDIACUN98+oUAFAQSR3Dp4aUYsXEEVr6wEmUKlNFHxo2lw8PDMXnyZDfWyKoyEyDjwBCAsq35+eefK+Lv33//RdmyZZ2yAMr28ZQpU5QzcsWKFdMlAEVUvfvuu+jdu7ciTuvWrQs5a+goORKArtSn9k2watUqZcv7p59+Us4l+nuiAPT3Gc44PgpAffPtNwLQikGEoJwTsqYzZ4AyXtaCtE7pW6RaSpNx4AjAwoULK2fc2rVrpzhuONoCXrNmjeIMcuvWLWW7tlu3bnjttdcwfvz4dGCuWgBlzb3wwgtYsmSJ8gM0d+7caNiwITp06KAIQtszftbGHAlAV+pz9Nm4efOmYuUUB5k9e/bY7Y+Wz5aZ8lAAmmm29PeVAlAfQ78TgLY4fv4ZeP55KDeK9OkD2OhCfdRYmgR8mIDaQ9C2688+Cxw7ZtxgKlYEli3TX791C1gsgCL+xPNWztLJ9qv8PysvYPkxKOLMmsQLeFmmDrkqAK11Hjp0CL/++qvifSvevFevXlW2qVeuXHnHwB0JQK31paSkZBiTPaujiOKwsDDIGcCtW7eievXq+ifBBDWorX1uAZtgEp3oIgWgE7DsZPVrAbhvX9otIseyrUKVtt9j1+hPkSfE/z3h9C0JljY7AbWHoD8IQPF2FeEnDh7ffvttlgLws88+Q7169ZS4cOKB+/333yvCcfjw4botgPbWSUJCAnr06IH58+dj+fLlyvarbdIiAG3z26tPQtaIaJUkAldCv4gDiTWJ4H3xxRfx888/K8L08ccfN/uS1tx/tbVPAagZpSkyUgDqmya/FoCCRjyFa47ohjPF5qH2XbWxuNNiJYagJxNj1BlPm4xvM1Z7CBo/G+5vIbMFUFqYOnUq3n77bcWhQ0SevTiAmb2ARTBu375dCZFSunRppaN6LYCZR7tv3z7Url0bH330kdI/PQJQymauT6ye8fHx6dWKIJRtcWsST+evvvpKOfcnjiWBlNTWPgWgf60GCkB98+n3AlDwJKekYMyWD/Hub+/i89afo0fdHvqoOVmat1Q4CcyF7GQceAJQrGOVK1eGxNCTs37ikGEbCFrE3tq1azOEgZF4eHXq1EHfvn2VGIF6BODZs2fTRaTtkv3hhx/w/PPP45tvvlHOCGoVgK7Ul/mj8tZbbykOMvPmzUOXLl1c+CSZu4iaAGQcQHPPb+beUwDqm8+AEIBWRHsv7EXNEjUzOInow8fSJOB7BNQegr7XY/Ue2bMASilxAhGLl2yFZhaA4vyR2QIoZSRQ8i+//KKEcJHgzVYLoDiHPPnkk3d0Rs7S2YvzJ9Y1uf3jpZdeQo0aNRRnlB07dmDChAlKjMFdu3bdEXxZLQyMs/XZdlYcW4YMGaKEyJFtaNskN4HIlrm/J7W1TwHoXyuAAlDffAaUANSHiqVJwBwE1B6C5hhFxl5mJQBFdIkF8OjRo5osgFKrOG3UrFkTr7/+uhISRgSgI3EkN3nIFW+Zk1gXv/vuO2VLWcLKiEWyTJkyyrk/OWNo7/o1EYAS7kXyZvYSdqU+2z6J4JU4ifaSeD9LMGh/T2prnwLQv1YABaC++QxoAbhkCbBmfQKmTw1GcJBx18hFRkYqlgYm4wiQ8W22ag9B42aBNZOAdwmorX0KQO/Oj7tbpwDURzRgBaBEhnimdQJWFWiPKqVLYN+YWcgR4jiIq6uoZVtm7NixrhZnOQ0EyJgCUMMyYRY/J0AB6OcTnGl4FID65jtgBaBguxp9C1XfbYvLhVbhqeIvYeVrXyI4u7Y7PfVhZ2kSMI6A2kPQuJZZMwl4l4Da2qcF0Lvz4+7WKQD1EQ1oASjoomPi8dTnHXAx215s77EdJfNxq1bfkmJpbxNQewh6u39snwSMIqC29kfMXYzR3dtJ8wXl69+ofrBezxCgANTHOeAFoOBLSE7AxZiLPnV/sL5pZelAJqD2EAxkNhy7fxNQW/uMA+hf808BqG8+KQD18dNUetasWcqF8UzGESDj22zVHoLGzQJrJgHvElBb+xSA3p0fd7dOAaiPKAWgPn6aSkuYiPr162vKy0yuESBjCkDXVg5L+RMBCkB/mk31sVAAqjNylIMCMAs6c9fswIKIdVgzYqg+wixNAh4moPYQ9HB32BwJeIyA2tqnBdBjU+GRhigA9WGmALTDLyUlFcXCn8a1wuvRLs90/PhOf32UWZoEPEhA7SHowa6wKRLwKAG1tU8B6NHpMLwxCkB9iCkAs+B3+PRl1Jz6KBLzH8GXj2zHK00fcpm03DwQGhrqcnkWVCdAxrcZqT0E1WkyBwmYk4Da2qcANOe8ZtVrCkB980kB6IDf5r0n8OGy77Fy6GBkz+76UpM7OgcNGqRvpljaIQEypgDkR4QE1AQgw8D41xpx/ansXxxcHQ0FoKvkWI4EfJSA2kPQR7vtsFtyF/Arr7ySnickJARly5ZFp06dMHLkSOTMmdOlYcn9uOPGjVPuA86bNy+uXr3qUj1mL7Rp0yaFw969exUGxYoVQ40aNdC1a1d07twZcs+w5HGUsmXLhhMnTuCee+5R8tqr76WXXkKXLl0c1vPTTz9hwYIF+PPPP3Hx4kWlvrZt22Lo0KHIly+fw7Jqa5+BoM2+UjP2nwJQ33xSAOrjx9Ik4HME1B6CPtdhDR0SAdi9e3f8+OOPKF26NG7cuIGff/4Zn376Kfr3749p06ZpqCVjlvPnzyviQkRJz549FRFZt25dp+sxe4ElS5agXbt2eO655xQWRYoUUQTx2rVrkZqaim+++QZyxELWlTW9//77ikBbvny5ksea6tSpgxUrVqjW54hZgwYNUK5cOYSFhaFMmTLYtWsX3n33XVSrVg1bt26lADT7gnNj/ykA7cMsDSAcQBHL2/KzVr4hT2XKTgHo5GKUL7vk1GReGeckN2b3HAF/FoBHjhxBhQoV0mE2bdpUEQU3b950GrBYqZo0aYL169fj8ccfd7p85gIpKSmKGAoKMuZOcd0dzKKCxo0bIzY2Fjt37tTchFhjhdupU5kfKYAr9dk2fOXKFRQtWjRDX0SEvvzyy6pzpbb2aQHUPMWmyEgBaH+aagM4CuARAMUALMhiNikAnVjmicmJ6LuiL6Ji4/B9x3mazwXyfJoTkF3MSsa3wak9BF1E7NViVgtgZgE4ePBgTJgwARcuXFC2La3p5MmTGDZsmGLFEh5iPRIrkli5JImAkTpl29JqwRKBIVvCkmbPno0ZM2bg8OHDyrajWKOkncKFC6e3kT17dmVbMn/+/JBA5KdPn1asYrVr14aImOHDhysWssuXL+Pee+9FeHi4Ymm0pq+++kqxam7btg0ff/yxklfaat++vdJWjhw50vOKQBOrm1hAz5w5o/SjYcOGigW0ePHiSj61MWc1gTVr1lQsoWK505ocCUBX6lNrVyyQ9913n2KNfOGFF7LMrrb2KQDVSJvrfQrArOfrVQCLAPQFMI4CUP/CTklNwVNznsPGc8vRPHgcVg7T5thBD1X97NVqIOPAFIByBlBEngguEXOSRCDJVm7JkiUxZMgQRRguXLgQc+fOxdKlS9GqVSvlrJoIngEDBigiSrYuRUiJUBNROXnyZLz55psQC+PZs2cVMSlnDsXaaG1HBODdd9+NihUrKnnlDKHUkytXLjz44IOIj49XzieWL18eq1evxqRJk5St6n79+in9tJ5rrFy5snLOrlGjRooYFKE6YsQI5W9JiYmJilVNzueJ4Hz44YcRFRWl1PnWW2+hSpUqmsac1Wfo1VdfVfoirDp27AgRcGrJkQB0pT619mbOnKlwEyulo216CkA1kv71PgVg1vM5G4D83FwIoBMFoHsW/rkrN1Dhw0cQH3Ie8x8+gi5tb1sE3NMCayEBfQTUHoL6avdOaasFUIS+bAHLGcDFixfjtddeU0RVnz590jsmAuSXX35RrHeFChVKf13EnFjj/v77b+U12cKU1zZu3KiIL0ly9k0E3XvvvaeIPmsSYfboo49Czss9++yzystWAXj8+PEM1roPPvgAY8eOxf79+zNsV/fq1UspHxkZqZS1CsBRo0YpQtGaWrduDbF0ylgliVVSLIfLli1Dy5Yt7U6A1jHbK3zp0iV06NABW7ZsUayhBQoUUJw+xFlDXreXHAlAV+pztKpEgIvoE3G9atUqhwtQbe3TAuidz69RrVIAZk32bgDnAJSX3QEKQPctwS37TiKsTSJaPVIZ8+a5r17WRALuIKD2ELRt49kFz+LYtWPuaNZuHRULV8Syzst015/ZC9haoViFZPvUNonjwNNPP405c+akvyzCZurUqUo4JrGeyVarPQEoZeTe7qNHjyrbotYk5cWS2KNHD0ycODFdAMoWrm078oZszYqX8rp16zL0S5xWxMImljzxsLWK2oiICMViaE1i5ZO+yravJLEObt68WbFEZpW0jtnRRPz111+KwJL+/Pbbb8q5ShGesr2dOTkSgNa8avUlJydnqNbe2cmYmBjF+ilb/NIvsbg6Smprn3EAdX8UfaoCCkB908EzgC7yO34cKF9erADqFfCeWnVGenOQ8W2Cag9BW9ZmEoAitsSCJl7AYmWSbVrZ/p03bx5efPHF9GHJ2TkRF7beqdY3xfJ27NgxxcvUngAcM2aMcnbPXpKtXwmLIlvJkqQusRKKxc82yZastGGvfaljw4YNiqjJ6lyjWB/lvJ9VIImV8vr169ixY0eWHxOtY9b6OZM1JJ7B0td9+/Yp5+9skxYBaJvfXn3Cz3oGU/62tcRKWdlCb9GihSKYRQBn7oO9saitfQpArSvAHPkoAO3PUwqA2775aXmElbxm66JGAeiBdS6/oMWqwGQcATJ2TQAaNyPurdmeWEpISECtWrUUi55sw+bOnVtptFSpUsqWrpzlsyfC5IybWOjsCUBZR3379lWEpe32sXU04p0q4tEqAEUsilizTRLGJDg4GNOnT7fbftWqVZXzgloFoGzFiseyIwug1jE7MyvilCJOM3J+UhxT9AhAKZu5PutWvLVeKxf5f1JSkuJ4I9vSMk/16tXT1HUKQE2Y/CYTBaC+qaQA1McvvfTl2Msomrto+gFxN1XLakjAaQJqD0GnK/SBAlmJJREVVg9dcYiQJNYpsQjv3r3bYYBoewJQhKQIkc8//zxD4Gl7CMSCZU8AigXvk08+wcGDBzN4JmeuQ6sAFG9hOePn6Ayg1jHbG4c4zcgWcub00UcfKY4hYn2T849aBaAr9dnWLaJdnHvkHOevv/7qVIgetbVPC6APfJjd2AUKQPswKwGoBaAkADlocechjrRyigCULYbq1asrL8jhXzHFy1aIeLNJkl+A8mu4WbNmyv/lELMcvJaDztYkv5wl/EH9+vWVl+QAs3jc2V6BJqE65Mvaei+ufEnv2bMng3VMvnDEM088+CSJp5tsf8gXgqS4uDjFE86X+rf75G5Mip6E4Y2Go1ftXj7XP1/nx/6lfYrc9fkQISSerWIZkwP9/pCyEksyNvGKlXh0EgZFgjlLOBZ5TUTN66+/rnjhXrt2TXHKEO9f65k9EYByVlDOu1mdQKQ+2daVM3hSVrZq5XtQ6pczfXImTl6TlJUAFBEiVkDZwh04cKAiKOUsm3wnikVLtrElaRWAYg2TNmUrVqyaMjZpY82aNUr9suWsdcz21oJYUUuUKKF8x0pfb926pVgcZYtdLG/S58zJ0RawK/XZ1i8OPfI8EXGd2elF5lSOAGSVMgtAeX7IWOS5IWnaurlI2Xdd/lkQwO3I1v7wIQnAMVAA2p90iU/ypzi6AZCDIw0sQjBzbloA3fChuR53HfXn1MeRq0cwofZKhIc1dUOtrIIEXCOgZgVxrVbvlnIkAGW7tnnz5opgkR+Pks6dOwfxrl25cqVyXlC2bsXxolu3bulXkdmzAFpHOX/+fCUOoIguOZ8mIWCefPJJRYBZHRHEaUFEilj8MicR37I1LGJPtm7lB7SIKzlX98Ybb6gKQDlXKMLPmsQhRNr54YcfIDeYyHjEKidC3xr/UMuY7c3iokWLINevidOG1C0WOAmH06ZNG2W8sl1tTwDK+UDxms6cXKnPtg5p216AackjoXFsPaYzt6229mkB9O7n2N2tUwA6JiqRRMVvvkkW2SgA3bQiD186gupTH0bQgc7YM3YGQkNvVyxWTVtrqZuaZDU2BMj4Ngy1hyAXDgn4KwG1tV87vC/2TvlMhk8LoB8sAgpAx5MoPzXXyI4sBaDxq/2bX07g5efKo8K92RARARSxXMQnW+bWLW3jexGYLZAxBWBgrnyO2paAmgBkHED/Wi8UgFnPZ3MA/wAQ17gjdryCpSQtgG7+PMyYAaxZA3z7LZA/v5srZ3UkoIGA2kNQQxXMQgKmJKC29ikATTmtWXaaAtA+mnYAhgK4YhF+ad4bdyYKQDd/HlJTAfmjJT6gm5tmdSSgEFB7CBITCfgrAbW1TwHoXzNPAahvPikA9fFzWFoOU++5sAcX9lxI96A2sLmArlq8/axe6gENggIw0Kc/oMdPARhY008BqG++KQD18XNYeuj6oZiwdQKGFBqC9/tnDBZrYLMBWbWEKrKGCgpIADaDVnsIBjofjt9/CaitfVoA/WvufV0AntKI+xaAqhrzujMbBaA7aWaq68ClA0p4mKDsQdjafSuqFa9mYGusmgTSCKg9BMmJBPyVgNrapwD0r5n3dQEowu41FeQyhmkWt3RPzw4FoMHE1xxbg5cWd0ON/Uvw7biHUaqUwQ2y+oAnoPYQDHhABOC3BNTWPsPA+NfU+7oAPA2grAbkJwGU15DP3VkoAN1N1E59c76+ip4vF0GdOsCmTfQONgK5RPq33lxjRP1mqlPtIWimsbCvJOAMAbW1z0DQztD0/by+LgB9nSAFoAdmKDw8HAULTsaoUYBcAjBN7L1MbiUgjOUmCCZuAXMNBC4BCsDAmnsKQH3zTQGoj5+m0mKdypkzF6ZPl4vqAT+5nlXT2D2ViRbA26TVHoKemhO2QwKeJqC29mkB9PSMGNuerwvAkRqH7y0XUQpAjRPkzmxHrhzBgv0LMKLRCOWeUSYScCcBtYegO9vyZF1yH/Ar8gvKkvLkyYPixYujTp066Ny5M9q3b5/+eZI7auVO2a+++gpdu3ZNLzNmzBjMmjVLuSu4Zs2a+Pvvv3HhwgX07NkTf/zxB65fv44pU6ak39fryfGxLf0E1NY+BaB+xr5Ug68/PbdkgvUwgGsAxDtYzgbKZWHbATTyElQKQC+A7/xTZ3y//3sMe2wYRjcZ7YUesEl/JqD2EDTr2EUAdu/eHT/++CNKly6N+Ph4nDp1CitWrMAPP/yAJ554Ar/88gty5syJhIQE7N69GxUrVkTRokWVIe/cuRMPP/wwBg0ahOeeew758uVD9erV0bdvXyxYsABSv1zZWL58eZQoUcKsmAK632prnwLQv5aHrwtAW9pjACQBGAUgBUB2AO8CCAYwzEvTQgHoAfCZY9TdiL+Bp795GhFnI/Bjhx/R7j65uIVJDwHGAbxNT+0hqIezN8taBeCRI0dQoUKFDF35+eefFQvg66+/jmlZHLK1lj927Jgi8qypSZMmSElJwW+//eaW4Yn4zJEjh1vqYiXOEVBb+xSAzvH09dxmEoAXAJQBkGgDVb4lxFP4Li+BpgD0AHh7t1RcvXUVk7dNxv8eehcrV4SgUycPdMSPm+BNIIEtAGX0bdu2xapVq3D16lVlW9d2C1isg5s2bVK2iOWGHvlbtoZFFFpfkzrk3ydOnMA999yDkydPYtiwYVi7dq0SW7FatWp49913FeuhNY0aNQrvv/8+9u3bh7feekvZRn7qqacgglTS4sWLMWHCBOzdu1cRhU8//TQmTZqEsmVvB4eQfjZs2BCtWrXCe++9p1g1pa2pU6fi0UcfzfCplTF8+OGH2LFjB5KSklCpUiUMGDAgw9b47NmzMWPGDBw+fFixcoaFhSl9KFy4sB9/A6QNTU0AMg6gfy0BMwnA8wAeB3DYZgok+PMmACW9NC0UgF4Cb212wAAoziGjRwPDvGUH9jIDNu9eAmoPQfe25rnaHFkApRdz5sxB7969FaEnAstWAB46dAjffPMNxo0bp4gz2eotVaoUzp8/j169eiE4OBifffaZIg7lTKEIyLp16yr5hgwZgmLFikGszHPnzsXSpUsVsSZJBJv8ka3mV199FfXr10f27NnRqFEjzJw5U9leltfbtWuHGzduKAJStq5FEObNm1epQ/op7d51110YPHiwsoU9fPhwRYDKnwIWrzFpV6ycjz32GF577TWlT//88w+uXLkCEaKSpLx4w7/55pto2rQpzp49q4hY4bF161a/P3OstvYpAD33efVES2YSgB8BaA9gIoAT8rkH8Jb8SATwjidg2WmDAtBL4K3NXr4MPPUUsGcP8OmnQJ8+Xu4Qmzc9AbWHoFkHqCYA16xZg+bNmytC7aGHHrrDCeSLL75QxJ7VwmflIIIqJCQEGzZsSEcjok3OE4oVrVChQumvi6i6fPmy4jxiFYBiAZRtZ9l+tqaYmBjlnGKHDh0gFjlrEute5cqVFYvcGxITyiIAZc6kX1ax99dff6FevXr47rvv8Pzzz6fnk7OJERERdqdQHF9EiIogFdFnTdu2bVMsiUuWLMGzzz5r1unX1G+1tU8BqAmjaTKZSQAGyQ80AN0sW8FnAHwNYLzlbKA3oFMAeoB6ZGSkYknIKl29CvToAXz8MbA7ZgVaVG6B7NnkiCiTVgJqjLXW4w/51B6CGcaYmgrcuGHssPPnl71V3W2oCUA5BtCiRQvFIUTEU2YvYGcEYJkyZZTtWrEqWpNY6WRbVpxIoqKilO1VEVsiAEV8SRlrWrduHZo1awb5W6yBtnU88MADiggUZxarABRnFBGc1iTnCCWwuVgs//e//ylCVLaFxYNZPJbtJasF9OjRo8oWtm2bYi3s0aMHJk4U+4P/JrW1TwHoX3Ov/1vFv3g4OxoKQGeJuZBftpDGjh2rWvKHf35Apx87oUvNLpgbNhc5gniQXBWaJYNWxlrrM3M+tYdghrFFRwMFCxo73KgotwS/VBOAVgG0efNmRYzpEYByXi85OVnZms2c5Jzg8ePHUa5cuXQBKIItKEh+46clsdy9+OKLdrlKeRGFGzduVN6XfooVct68eRnyy1aybO2OHDlS2b6Vc4IiEp955hm79UqIG9k6tpesZx5lC9ufk9rapwD0r9k3qwCUftv2XbyCvZEoAL1BPYs245LiICFilhxagnbV2uHHjmkWAiYScIaA2kMwQ11+ZAEU5wxx2JAzcZmdQGTMzlgA5XygiDQ5U2dPBEoMQdk2tloAExMTlbN/1mS1Roqou+++++6Yvvz58ytWQK0C8N9//0VoaKhDC6BYB+XMoTCw3ba2Ni7hcES0+nNSW/sUgP41+2YSgBKMajqAppb4f7Yzcfuno2fnhwLQs7xVW0tKScIbK99Aqyqt8Exl+7/0VSthhoAmoPYQNCscRxZA8baV83YDBw5UtjntBYJ2RgBKwOnt27crsQTFKSOrlJUAFIcPsUKKM4bkcZS0WAClvIS+kcDXWZ0BFKtk1apV8fnnn2fwCjbrfLvSb7W1P2LuYozuroTdErN3tCttsIzvEDCTABT7vvzkkzN/3wB4CcAQy78/8RJSCkAvgdfSbEoKIOfKJURM48ZaSjAPCaiHwjArI+tNIIsWLVLElWy7ilOFbIvKa+KgIY4Osn2rVwCePn1aCRot7Yhzh8QNvHbtGvbv3684a1jPBmYlAIWxCDEpKw4lcjaxYMGCileueClLWBpb5w61LWCpb9myZYo3sdULWMTgwYMHcfHixXQvYHH+kHOK0m7jxo2Vc4TCSM4iytlBec2fkyMBmJiciNxvl0XyVInIRgHoD+vATALwHAC5CUTi/l0HIK5lIgjlUEZDL00GBaAHwMvWjISncDYdPAg89BAQG5sWJmbQIMBml8nZ6vw6v6uM/RGKmhXErGO2WgCt/RdxI16xEq7lhRdeUOIAWpMIQLGYyZk361VwjiyAIhrXr1+fAY1cFydn8FauXIlLly4pN4rUqFED3bp1Q5cuXZS8IgA/+OADRYzabgFbK5K4hOLxK169sk0snsGytfz2228rW7qSpJ8i6mR8tknOFEr7I0aMSH9ZglVLe3KriSTx+hUro/TJmt+BdroAACAASURBVObPn6/EAZTYhHL2T0LAPPnkk8p29t13323W6dfUb0dr/8rNaBT7sCAwTqmKFkBNRH07k5kEYJRl0QlR26DQtq97mjYFoAeIy1aSxAdzJR0+DLRvD1yKP4PH3huMmWHTUDRP2tVWTLcJ6GHsbxz9VQD62zxxPO4n4Gjtnzh/DRWmF6EAdD92r9VoJgG4D0AbAEcB/AFgGoArli1gb/0sowD02tLV3rBYAEeunoBJe/+HsgXKYmH7hWhQtoH2CpgzoAhQAAbUdHOwNgQcrf1/z1xB1U+KA+MVz25aAP1g5ZhJAPYFILeByB1Bslfxg8UTWM4BSpBobyQKQG9Qd7HNubvmou+vfdHhvg6Y1yZjyAgXq2QxPyRAAeiHk8ohaSLgaO1fuACUHB4KzFEu46IA1ETUtzOZSQBmJlkaQL5MV8N5mjYFoAeIyzVU1vM+epvbf3E/yhcqj3w5ZOkA4ijCc4GAOxnrnSNvl6cA9PYMsH1vEXC09uW7Mlf/UCR+SgHorflxd7tmEYAhlu3eEgDi3A1BR30UgDrgaS06fvx45fYAdycJ4Sb30tesCcg5cQfRKtzdtM/VZxRjnxuohg5RAGqAxCx+SUBt7TMMjH9Nu1kEoFCX+3+rA4j1oSmgAPShyXC2K3KF3JNPArt3p6LAC73wYdcw9Hu6pd9f+O4sp0DLr/YQDDQeHG/gEFBb+wwE7V9rwUwCUM4A1gIwAEC8j0wDBaCPTISr3UhMBIZOPIpJMfcjNSQGzSs1x5zWc1C6gJwwYApEAmoPwUBkwjEHBgG1tU8B6F/rwEwCUOL/lQQgLkgX5fiWzVTcvrnbs/NDAehZ3oa1dvLqWYzYNBjrjq/D4dcPo0BOmVqmQCSg9hAMRCYcc2AQcLT2L8dexsQNEzG+tdzFQCcQf1gRZhKAtyN13kk+YwRQz80MBaAHWHvyfNq1W9dQOHfhDKNKSABy5PDAQL3YhCcZe3GYmpqmANSEiZn8kICjtX/g0gFUn1ydcQD9aN7NJAB9ETsFoAdmxZseqmfOAPfXTUbVV8dictdueLhaWQ+M2PNNeJOx50fruEXrQ1CuMytQgJZgX5sf9sc4ArL25eaTqKioO9b+9uP70GB2LQpA4/B7vGZfF4BiirmmgYpcCyfXw3k6UQB6mriH29uxA+g0bBVONmwBJAej431d8FnYFBTJXcTDPWFzniIQFxeHe++9F5GRkZ5qku2QgM8QKFmypHJfs1wVaJuW/n4Yz30dxjiAPjNT+jvi6wIwGoCWn+BXAXjjiUwBqH8N+nwNKSmpGL1wFeb++xFSCh7H0f5HERIkkYmY/JWAiEC5n5aJBAKNgNzrnFn8CYO//gIenMlA0P60HnxdAIq3r3LiVCW9qVEoqtXj7PsUgM4ScyG/L91TeyX2it27hM9duYG7i+Z3YXS+UcSXGPsGEff3gozdzzRzjWRsHGPZDXl4NgWgcYQ9X7OvC8DfLF6/amTkp3oztUwGvE8BaADUzFXOmjULvXv39kBLrjURG5eI/CPLIW9CBTxTtgs+6NQJlcsUda0yL5XydcZewuLWZsnYrTjtVkbGxjFOTgZyv8GbQIwj7PmafV0Aep6Icy1SADrHyy9z7z9xAa0+GYj/ci0FcsRiUq0NCG/zhF+OlYMiARIIXAI5+4UigVfB+c0CoADUN5UUgPr4+VXpyKs3MfanFZjwcnvkCAnKMLaT108iLjEBocWr+NWYORgSIIHAIUAB6F9zTQGobz4pAPXxC5jSb656E9MipiEkugrqBL2AYQ1H4tlnA2b4HCgJkIDJCRy8dBDVxzdA6qQoGUlBAOKkyWRiAhSA+iaPAlAfP02lhwwZgrFjx2rK66uZNp3chOELF2Db5RVIPtoE9SO/xrZtGXsblxSHXMEZQy94ajz+wNhTrFxth4xdJae9HBlrZ+VszrXH1qLpnKaMA+gsOB/OTwGob3IoAPXx01Ra4rFJbCp/SBJSZu+hWCTczIuHHro9otTUVJSeXBqFchVCkdj6aJS7D55vVA/VqwNBGXeTDcHgT4wNAeSGSsnYDRBVqiBj4xivOroKLb5oQQFoHGKP12xGASh9FjVw3uO07myQAtAHJsEfuhCbGIsey3pg25ltkPOCmP8LcKQlwsOBSZPSRih59l/cjxolaiBPSB5/GDbHQAIkYBICq/fvxOBfhmD3kPXSY24Bm2TeHHXTTAJQnnhTAMidwMkA8gIIA1ADwIdemgsKQC+B9+dm9x6PxKE9BfF3RG40bQo0aZI22s3/bUbjrxoje7bsqFykMh498hvKFCqJqlWBRo2AMmX8mQrHRgIk4E0Cy5YBYWFy7E+0HwWgN+fCXW2bSQDOAFAZwHsAVgCQ69/kYtaVFhHoLibO1EMB6AwtF/OuXr0azZp5I8yjix02qNjxa8fx3b7vsOfCHhy8eBgH39yFlKS0/eH584EuXYD5e+djzq45qFq0qmIp7Ffvdfz5J1CuHFC8OJAti088GRs0aTbVkjEZG0/AuBZ+/hlo25YC0DjCnq/ZTALwNIDaAOTaN9ur3+SuYLkz2J3pUQDPAdgPoAGAdwDcsNMABaA7qWdR18KFC9GpUycPtGSuJmJjgSNHgH//BR55BChdGpi0dRLe2/QebiTcQK27auG3TntQxHJJYu7cwIoVwIG8M5St5LIFy6JuqbpoXqk5FixYiKee6oSiRYHs2c3FwSy95To2fqbI2DjGP/4IdOhAAWgcYc/XbCYBKDezyyZXko0AzA3gGIC73YhO3DAPAKgJIAbAyxbhOZAC0I2UWZVhBMShJPJmJKLio1AyOBSffQb89x9w8iQwfTrwzt9hWHZ4mdJ+5xqd8V2775T3y5cHgoMB8bep90FPpOS6hJL5SqJ1ldZoWaWlkv/WLWDvXqBYMSAkXzTKFM+P7NnN9DViGHZWTAJ+TSApCYiMjEbZstwC9peJNtM393IAawFMtxGAfQE8BaCtGyekOYBhAB6z1FkJgJx6LUcB6EbKrMprBEQgXoy5iNPRpxVnkvuK34ezZ4EPPgDOnwfOnQNSuz+K3ZcikJyajPcffx8jGo9Q+rt/P1BTfhpJGlwQwbljUSRPEUxvPh2daty20ko9H02Ow978E1AoV0E0ebQA2tVurghK25SQmIyEpGTky53DazzYMAmQgDYCtcP7Yu+UzyQznUC0IfPpXGYSgKFyDh7AvwDqAdgCoI5li1Zec1fqBaAVAGuY3mIALgCQ7V6xCNombgG7i7qDeuLi4pArl3fi43lgeD7RhD3GKakpuBJ7BSFBIUp4GkmRkcC8ecDlK6lYktIDFWtexo2kKxjZeCSaVmyaPpaICKD+0+eAt0qnv7b55c14rJz1d1Xay/PW/Yluf9QDkoOBpDzY0mUPGtYon4FJxJkIfLBuCv79Jw9yBeXB00nTMXxYdhS2Ofix78I+7Di7AwlxOXHyUEE8Wqw1cuZMc46RrW9JF25eQHR8tBJrMUe2vAhOLIKQkDSrp+TJ6nykuyaI69hdJLOuh4yNY6zEAZz5HDAxlgLQOMwerdlMAlDAFAXQ1eIMIlvCcwHI2UB3piEWYdnRUqmIvOuWbWZpkwLQnbQ11BUeHo7JkydryMksrhJwN+OEBODk6QTsPnsQF6OiUbpCFB6v8CgK5854XHf1n/+i94JRiE+JRUJqDLa+tQBVy8pvrttp4f6FeP6n59NeSM0GvJeMM2eyKWcerWni1ol4Z60c1QUQVRaYckr55+nTt72jB64aiKkRU5XXaxRsgP0Dt6aXP3Mm7Qxlpx87YcmhJQjJHoL297VHv9JfoXVrpAvFrVuB8K2dceDSAQRnD0bXWl0xoP4AxTL65pu38xXq2hPxiFLE80u1XlLOWQrj/v0nK9vwwSGp2JpzOBo+Gox8uYPRoXoHhBaT37hp6cIFYPHyGOxJWIyQoGDUfygELUKboEhuy4FOS75zV25g3a7DyBkSjFwhIXi6bhXkyRWSgV9iciJiE+IRki0ncuUI9uste3evY1c/T/5Ybu6uuej+Q3fGAfSjyTWbAPQEetlWfgJAB0tjcr5QRKaYoBIzdYAWQA/MCH/VGw/Z1xknJCXhctQtXImKRf7sdykhb8RyZ00SO1Guqrp4NR4H9oegekhLxMcDL7wA5LGETFzx7wr8fup3xCfHI3dCWVxZMRCJiYCcbZo2DShUCBj3+zj8cfoPJKUk4bF7HkPLAkPRrx/S861eDfTZ2AF7IvcoeXrU7YGhjw2FCMNnnkmrS/4U+uBuXIhNC1Uq2+P9H+4PYbxjRy40bgwgeyIw8va29+KOi9GmWpv08WzZAjR69gTwZoX013b02IF6pWXz43b6ZPkW9P+7UfoLu148h/srlsqQZ8G+BeiyuEvaa8nByD0pAceOZkMpm2zf7PkGn/75KZLjc+HMgTJ46Mw3EKO7nB+1WlqXH16uWFnFglokewUEHeiMvHmh/GneHEr+U1GnkJCcgLwheVEwV0GPx6v09XVs/CfZuBbm/D0HPRf1pAA0DrHHazaTAOwH4A8AuwHUBbDEEg+wPYC/3EjuSQDvArB+q4rn8a8AbOwN6a0pArBp06aoLlc2iHJ84gls3LgRY8aMSd+2FM+0QoUKpYcykWj106ZNy3C92axZs1C7dm3Ur19fqefQoUNYunQpBg0alN7Y+PHjERYWhtDQNEvB9u3bsWfPHvTu3Ts9j1yFNGDAgPSbMyT0xPXr19O9aOULcujQoewf+XH9ATDy87FkyRKEvxOuxG0Ua6F8flu0CEORIqFISEjFL5t/xKWL/6BDt7YoV7CcIpisn98cOUpi09YY/LBpAmJu3kCHXg+hZejTyJMtT4bPb8TB03jtvcHIliMEBSveiwX93wISbmb4ftl+ZjteGx6OC/H5EFSyMGr/uxCjRh3Chg23v19m7JiBoe8PRUrl7IhNKIfsn+9GUtJ2TJy4B2+9lfb90ueXPpg5fibwMHB/mcexe+BGIahskFy82EkJM9Ty65b4deavwJNAl7pd0KvofDz11ELkylUIhQs3w99/A32WdUDETxG4/8X78XyN59GlZheMGTMLZ8/WRq1a9VGwIHCz4Nf4e3ME+g7si3sK3oMCOQso/Pj9J8cwPPv82LRpk/Lj5Uz0GSz6eymwK0GWA88AulF0eKsqMwnA4wBEHV20xAE8BOCmRaiJxc5dSewK4ll8PwAJMfMGAHEEkb8zJ1oA3UWd9ZAACfgUgeTktJBA1rOR4jgkZ0Llzuq4mFy4cbwaYmKg/OnYEciRA/h277dKiKGYhBg8ePeDqJXaDSNGpOWRPxs2AM8tfgq7I3crXuojGo1Qzo+uWgW0aHF7+AXGFER0goQcAWa1moVeD8jRbGDTJuCVV4AiJeLwX4M2eLJBMZQsUBR9HuyDqsWqplcgIZKuXE1G4aJ0MHLXopJIAZXGhSJp5mEKQHdB9XI9ZhKAUZZfHSLQLlvO5MlPEXHQkLOB7kwSdVi8gfcAkJiA4YwD6E68ztXF2F7O8XIlNxm7Qs25MmSckZd4o4ujUVD2IERHy64HEBVl+VNhLqLiryEqLgrPhT6HOqXE3w+QrfEBA4ALMZE4174ikENxSIDVwcjKWLm1ovcu4LW6yHarGKrfczdWvLBcsSbapjOXonH2chRqVyylnI9kyprAxx8Db+wNBeZQAPrLOjGTAJQDNfdZbv2YYLEGymnnKxYPXW/MCS2AHqDOGxSMh0zGZGw8Afe3IPdji1WyeN7iytlE6zretQuY9PVhbMo+EjdxHkXLn8POnjvvcELqP2sBPonsAqQEITimHGLG/oscIWm361iTWDNzBOVQHHoCOU2dCgwcyEDQ/rQGzCQAP4Vy+gT5ZFcAgLiFPgDgK0vQZm/MCwWgN6izTRIgARJwA4G5a3Zg2sb5uBh3Cgmpsbg8Vc40ZkzDNwzHmC1jEXSjHEreaoIu+eZg/Hg3NG6yKiZNAt5+mwLQZNPmsLtmEoDy86sbANn2/RZAisVbtwSAhV6aFApAL4FnsyRAAiTgCQJy//aU1T9i37kjSIyshKq7f8YBuSvKJr3323u4EHMBJVEHoXkaonX9aunxJz3RR0+0IaGdLl7kTSCeYO2pNswkAD3FxJl2KACdoeViXvF6Kyn3kzEZRoCMDUObXjEZm5+xOMZcuQKUELODTWrwRQOIt7WStoYjaP0kPPccIPfnWpOcecxmdLRxAxF//tfnmLRxEv59R7l3gV7ABrL2VNVmE4DVADQBIB8/276P9BSwTO1QAHoAvITGGDt2rAdaCtwmyNj4uSdj/2Us4k7CpHy3cRf+21MO53fXRpUqSN8qlmDc5aeVR80SNdHwnoZoVuol1K1QDkEZjxsaD0hHC6M3j8aIlSMYB1AHQ18raiYBKPH+vgMgxndxBpG/Jfje75atYG+wpQD0BnW2SQIkQAImInDuxjl0XNQRO8/tVAJlV9++FZE7G6BpU6BnT4kf6/uDkW3uUatHUQD6/lRp7qGZBKAEgJ5muf5N4vPJnVL9ARQHQAug5ilnRhIgARIgAW8QkBiKEWd2YO2Xj2DlimAlMPbMmYDE8pfYiLvO70Lbam2VoOC+ljad3IQRiz7Flv/9IF3jFrCvTZAL/TGTABT3IxF9yZa7eeV2erlLSQJEl3Fh7O4oQgugOyiyDhIgARIIQAJy57NcoSe3n7z2y2uY9dcs5AzKqQS+nt5ius8RydkvFAmfMg6gz02Mix0ykwCMBHAvgFsW0feY5aYOCQSd38Xx6y1GAaiXoIbyck2e7XV3Goowi5MEyNhJYC5kJ2MXoDlZxMyMr926hsUHF+PL3V/igVIPoHOh6QgPB+Q20LCw2zeyOInEbdklUHfNabwJxG1AfaAiMwnA5QC+sNwB/DkAufdHwsDLVe9yvbo3EgWgB6jLncfWO5I90FxANkHGxk87GZOxVgLiNPL13BD07w/ExQGNGwMbN3pXBL7zDjDxOm8C0TqHZshnJgEocUDEZ+osANn+lVCcIsCGW+7u9QZvCkBvUGebJEACJBAABCIjgY8+AgoXBpq+EoHRW0bj4xYfo3yh8h4fvVgjp9ygAPQ4eAMbNJMANBCDy1VTALqMjgVJgARIgAS0ElC8cDeNQu7g3IoIfLXuq1qLuiWf3ME8PZYC0C0wfaQSswnACpbr3zKf+fvSSzwpAD0A/tChQwgNDfVAS4HbBBkbP/dkTMZ6Caw7vg6vLnsVPev2xPBGsvnluTTg13BM3zIVGJcqjdIL2HPoDWvJTALwNQCfALgKIMaGiKxGEYbeSBSAHqA+fvx4DJKT0EyGESBjw9CmV0zGZOwOAtfjriN/jvwIyh6E779PiyF4113uqNlxHYqX8h+zGAfQeNQea8FMAvA/AAMBLPYYHfWGKADVGTEHCZAACZCAmwnIfcQ1agD33gusW5f2t5Gp57KemLNtDgWgkZA9XLeZBOB1i/OHhxE5bI4C0Jdmg30hARIggQAiMHt2WhDpCpWSsSniGkoXLmbY6P849Qd2ndyF/o3l/gVuARsG2oMVm0kAzgcg4V82eZCPWlMUgGqE+D4JkAAJkIBhBObMScUX0c/jRq5/sPXVrSiQUx5LxqTo6GgUlKjVFIDGAPZwrb4uAN+34SGruhuAnwCcy8SJV8F5eOF4sjmenTKeNhmTsfEEjG8hUNfx+N/HY/D6wWhXrR0WdViEbNmMebRTABq/hj3ZgjGrxH0j2KihKnECaaIhnxFZaAE0gmqmOuk9aTxkMiZj4wkY30KgruPU1FQ8/9Pz+P3U74joEYEyBdx/O+quXcCJE9Fo144WQONXsmda8HUB6BkKrrdCAeg6O5YkARIgARJwE4GbCTcRkxCDu/IZ4xLcsSOwaFG0ZfeXZwDdNG1ercYsArCRJf7fDgB/eJVYxsYpAH1oMtgVEiABEiABYwiEtY/BsiJ1gNlHpAHGATQGs0drNYMA7A5gDoArAAoDeBnAtx6llHVjFIAemAjeoWo8ZDImY+MJGN8C1zGQkgLMnw8UKACEhbmPeamBYYjMtYxhYNyH1Os1mUEA7rXc+ytewC8BeANAPa+TS+sABaAHJmLWrFnoLbEOmAwjQMaGoU2vmIzJ2HgCwLVrQIUKaQLw0CEgd273tNrim1ZYdWAFBaB7cPpELWYQgBL/Tyx/4uwRBOA8gBI+QY8C0Eemgd0gARIgARKwEpg4KRnvfDsPbTsk4aehPd0CpsX8Fli1fxUFoFto+kYlZhCAcurUNrCRXAVXxDfw0QLoI/PAbpAACZAACVgI3IhNQInRlZAzTzzODzqJ3CH6zYB/n/8bD497DknTT0srPAPoB6vNDAIw3rIFbMX9DoAJmdgzDqAfLEYOgQRIgARIwD0EpkdMx4BVAzDjmRnoW6+vWyrN2S8UCZ8epgB0C03vV2IGAfibZfs3K1qMA+j9dWRoD4YMGYKxY8ca2kagV07Gxq8AMiZj4wncbkFCwrRa0Arh9cPRumprtzRdO7wv9k75jALQLTS9X4kZBKD3KWXdAzqBeGB2IiMjUbJkSQ+0FLhNkLHxc0/GZGw8AWNb4E0gxvL1dO0UgPqIUwDq48fSJEACJEACJiCwbZt4GEejZUveBGKC6dLURQpATZiyzEQBqI8fS5MACZAACZiAQO0njuL8+WhcOvyA9JZOICaYM7UuUgCqEXL8PgWgPn6aSq9evRrNmjXTlJeZXCNAxq5xc6YUGTtDy7W8ZHwnNwkMLda7uDjgySdd4yql8g58EHFJF5HyCb2AXafoWyUpAPXNBwWgPn6aSi9cuBCdOnXSlJeZXCNAxq5xc6YUGTtDy7W8ZHwnt/h4oNg9l1HioY04tryDa2AB5BlYBwlJ15H8yUlaAF2m6FsFKQD1zQcFoD5+LE0CJEACJGAwgUpvv4xj+b/GumeP4Mk6lVxqrcaMWoi9GYMTg49TALpE0PcKUQDqmxMKQH38WJoESIAESMBgAu999ytGHWmJrqU/xNc9hrrU2olrJ9Bi9HAcnvwdBaBLBH2vEAWgvjmhANTHT1PpuLg45MqVS1NeZnKNABm7xs2ZUmTsDC3X8pKxfW6x8QkoO7UUWlRugW/bfusaXAAMBO0yOp8sSAGob1ooAPXx01Q6PDwckydP1pSXmVwjQMaucXOmFBk7Q8u1vGScNbfo+GgUyGl7q6rzjCkAnWfmyyUoAPXNDgWgPn6aSvNXvSZMujKRsS58mgqTsSZMujKRsS58qoUpAFURmSoDBaC+6aIA1MePpUmABEiABExCgALQJBOlsZsUgBpBZZGNAlAfP5YmARIgARIwAYGJP63H4N9eRfIn/0lvGQjaBHOm1kUKQDVCjt+nANTHT1NpxvbShElXJjLWhU9TYTLWhElXJjLWhc9h4aD/lUJK6nVgYhwFoHGYPVozBaA+3BSA+vhpKs3o/pow6cpExrrwaSpMxpow6cpExo7xRV69idofPI9yee7Djg8/cop19kElkD+6MqJnbqUAdIqc72amANQ3NxSA+vixNAmQAAmQgIcIpKSkIsegexCUnA/xkw861Wr2wUVRIKo6omZuoQB0ipzvZqYAtD83pQGEAyhiefsqgGkATmXKTgHou2ubPSMBEiABEshEoPLbr+Bo/q+ws/MZPFhFHnXaklgPY27cQKXyd1MAakPm87koAO1PUW0ARwE8AqAYgAVZzCQFoAeWeGRkJEqWLOmBlgK3CTI2fu7JmIyNJ6DewoSf1mPNge2Y+lJPVC9fQr2ATY7o6GgULCj+H3QCcQqcj2amAMx6Yl4FsAhAXwDjKAC9t4KHDBmCsWPHeq8DAdAyGRs/yWRMxsYTMLYFCkBj+Xq6dgrArInPBtATwEIAnSgAPb002R4JkAAJkIAvEaAA9KXZ0N8XCsCsGcpBh3MAygM4SQGof7GxBhIgARIgAfMSGDF3MUZ3bycDYBxA805jes8DUQAOBFDKztwJi1QAEwFc1Di3PAOoERSzkQAJkAAJmJNAamoqZv89G/3mjUbS9NMUgOacxjt6HYgCUMvUpVjEoG1eq0AMsnmRAlALTZ15Zs2ahd69e+ushcUdESBj49cHGZOx8QS0txAVF4XDVw7jodIPqRZKTEpGjg+DgRv5gUk3KABViZkjAwWgvnlSBGDTpk1RvXp1paYnnngCGzduxJgxY5ArVy7lNYlOX6hQITRr1kz5v3gDTps2LYNjgzwcateujfr16yt5Dh06hKVLl2LQoEHpPRw/fjzCwsIQGhqqvLZ9+3bs2bMngziSg+YDBgxI95qVwKjXr19Hp05pxxjlsvShQ4eaqn+jR49GvXr1yA+AUfMr66xNmzZcfwZ+Pvr164f+/fsH3OfXk99/Xbt2xUcffeRX339G8as1tBb+PXAEUYtuIGdIsPJ8sPf9smnTJkTfjMGMndOBuBBgdyIFoD7d4DOlKQDtT0UlALUASOyRZACzspgxWgB9ZimzIyRAAiRAAloJPPHBcPyW8iG+bBCBV5o6tgJeuxGHImNKAvEApkRRAGqF7OP5KADtT5CY3f4EsB7ADgANLEIwc24KQB9f4OweCZAACZDAnQQm/7wRb+1tglZBU7F8+ACHiG7eBPLnB7L1DkXqrMMUgH6yoCgAHU9kDgCrADShBdB7K162w63b3t7rhX+3TMbGzy8Zk7HxBLS3EHUzHiWqnEKjGpWwdo1jKZCaCogILDo4FImfUgBqp+zbOSkAHc/PGwDWyJE8CkDvLWQ5+2h7FtJ7PfHflsnY+LklYzI2noBzLWzfDsjxdbHuaUkMA6OFknnyUABmPVfNAfwDIDeAI3a8gqUkt4DNs9bZUxIgARIgAR0EGAhaBzwfLEoBaH9SJNLlUABXLMIvzX33zkQB6IOLml0iARIgARJwPwEKQPcz9WaNFID66FMA6uPH0iRAAiRAAj5O4HrcdUzaOgl1i9RF2zptpbe8CcTH50xL9ygAtVDKOg8FoD5+oTFNYwAAEKpJREFUmkrz7JQmTLoykbEufJoKk7EmTLoykbHz+Dae2Igvdn2BT575BIVyFbJbwfFrx1FxekUMqD0A09pMowB0HrNPlqAA1DctFID6+GkqTe9JTZh0ZSJjXfg0FSZjTZh0ZSJj5/HN2DEDr698HUs6LUFYaJjdCnafPoI6X1ZB93vfwpfdJlEAOo/ZJ0tQAOqbFgpAffxYmgRIgARIwIsE/rlwCDVmVkOzQv2xasB0uz1ZveM4mn8VhqcKvIR145XbqbgF7MU5c1fTFID6SFIA6uPH0iRAAiRAAl4kkJKSijy9mqPg5ea4sGSg3Z7s3QvUrg1k7x2KFAaC9uJsubdpCkB9PCkA9fHTVFruPLbekaypADM5TYCMnUbmdAEydhqZ0wXI2GlkSoGXXwa+/ho4eRIoV+7OOnbtAurWpQB0ja7vlqIA1Dc3FID6+GkqPWvWLPTu3VtTXmZyjQAZu8bNmVJk7Awt1/KSsWvcFi8GBg8GvvwSaNjwzjpSUoDYWKDIIN4E4hph3yxFAahvXigA9fFjaRIgARIgAS8TEIGXLVvaH0cpZ79QJPAqOC/PlvuapwDUx5ICUB8/liYBEiABEjAJAQpAk0yUxm5SAGoElUU2CkB9/FiaBEiABEjARwjEJ8VDgj7fle+uDD06cOkAJmydgG83rEfS9NPyHr2AfWTO9HSDAlAPPd4FrI+extJDhgzB2LFjNeZmNlcIkLEr1JwrQ8bO8XIlNxm7Qi2tTFRcFO6Zeg+eqfwMFrRbkKGiNcfWoNm3zVD6ZmOcnbiJAtB1zD5VkgJQ33TQAqiPn6bSkZGRKFmypKa8zOQaATJ2jZszpcjYGVqu5SVj17hZSzX4ogH2XtiLC29fQL4c+dIr+/XIr2j5XUvMaDID/Rr1owDUh9lnSlMA6psKCkB9/FiaBEiABEjARwhYbwWZH7YIXe5vn96rFbt34N3f3kW/+3uj+xNtKAB9ZL70doMCUB9BCkB9/FiaBEiABEjARwhcunkVT3U8glKpD2HVytvy4Pvvgc6dJVZgNLp1k+N/PAPoI1OmqxsUgLrwgQJQHz9NpVevXo1mzZppystMrhEgY9e4OVOKjJ2h5VpeMnaNm22p1q2BX38Fzp4FrCdvvvoKeOUVYOHCaHTqRAGon7Jv1EABqG8eKAD18dNUeuHChejUqZOmvMzkGgEydo2bM6XI2BlaruUlY9e42Zb68UegQwfgww+BoUPT3pk5E+jTB1i6NBphYRSA+in7Rg0UgPrmgQJQHz+WJgESIAES8CECSUnAF18AXboA+fOndSw5Gbh1C4iPj0axYhSAPjRdurpCAagLH7eA9eFjaRIgARIgAV8lkJqaimw214PUDu+LvVM+k+4yDqCvTpoT/aIAdAKWnay0AOrjp6l0XFwccuXKpSkvM7lGgIxd4+ZMKTJ2hpZrecnYNW72Sr295m1EnI3A5pc348tdX+L7f77Hph0nkPjJMQpA92H2ak0UgPrwUwDq46epdHh4OCZPnqwpLzO5RoCMXePmTCkydoaWa3nJ2DVu9koNXDUQUyOmYtnzy7DhxAbl3yGXKyDxk+MUgO7D7NWaKAD14acA1MdPU+nly5ejtbimMRlGgIwNQ5teMRmTsfEE3NfCuRvnUGFaBdxf8n5ULVYV8/bMQ8ilykiccYQC0H2YvVoTBaA+/BSA+vhpKs1f9Zow6cpExrrwaSpMxpow6cpExrrw3VH4q91f4cG7H8T20zuw/exWzFv9OxI/PUwB6F7MXquNAlAfeo8KQD0xrpwtq5bf0ftZvWfvdS2vefJLXW3cjpaLs2XV8pPxnbTVmNmW0JKXjMnY1UeAlvWVVd3OlNWS16h1nJoKhIWtxvLlzfDnn8AjX4YigQLQ1SXjc+UoAPVNiSIAT58+jQIF5J/GpqFDh2LMmDEuNeJsWbX8jt7P6j17r2t5rUOHDli0aJFL43a2kNq4HdXnbFm1/GR8J201ZrYltOQlYzJ29jvCml/L+sqqbmfKaslr5Dpu0GAoDhwYgxo1gMP1H0DinKO0ALq6aHysHAWgvgkpDeCMvipYmgRIgARIgARMRaAMgLOm6jE7ewcBCkB9i0L43Q3ghr5qWJoESIAESIAETEFAwkOfA5Bqit6yk1kSoADk4iABEiABEiABEiCBACNAARhgE87hkgAJkAAJkAAJkAAFINcACZAACZAACZAACQQYAQrAAJtwDpcESIAESIAESIAEKACNXwOFAKwAkAAgL4BpAOYb32xAtVAXwHQASQCCALwFYEdAETB+sPJdsQlANQByG/xI45sMmBY6ARhgOVQv3w8/BMzIPTNQrl3jOfM72HjGbm+BAtDtSO1WmB1ACgAJFngQgISPYXIfgbsAxAC4CaA6gC8BPOy+6lmThYCs2ycBVKIAdNuakO+EbQAeBCDfxzsB1GdkAbfxtVbEtet2pBkq5HewsXwNqZ0C0BCsWVYqX0JzALTwbLMB1VoVAHMBPBpQo/bcYLsBqEgB6DbgTQG0A9DbUuMsiwVwvdtaYEVWAly7nlkL/A72DGfdrVAA6kaoqYJiAH4CcB+AIRYRqKkgMzlFIMSy3T4BwFqnSjKzVgJ8iGolpS1fZwA1AQy1ZB8HYBeAhdqKM5cTBLh2nYDlYlZ+B7sIzhvFKADvpC7ncfoBqA0gHwBZ0LJ9a5veA9DDsqX7lyX/P5YMrwGQOq4DaJOpnAhByS91y/uBmoxgLNvsCwCstmwBBypb67iNYCx18yGacWXp5SwWwLYA5HtDklgA5d7DdYG+gG3Gr5extSqu3awXlTsY8zvYZB9aCsA7J+xpAEUA5LFY6jILwHcAvG7Zxj0G4F0AXQGI2TvWzvznsDiAyFtSp/y6rwUg3mRrxZ3ddTdjWcffWMT1FHd21MR1uZux7UNUzgCOMDEbd3ZdL2c5A/g7gIcAyAM0AsAjPAOYYYr0MubaVV/xehnzO1idsc/loADMekoaA9hgxwJ4HMBkAJ9YiorX6XkAA7Pw7q0HYKLFQzWnxQtYfuEzAe5iLL9ev7B4/sqatmd9DVTe7mIs/L6zbFfKDxn58SNnWZMDFWymcevhTC9gbYtID2OuXWMZ8ztYG1+fykUB6JwAlF/rIi4aWH6pW0vLtuM+AG/71Oz6fmfsfaGTsXvnjYzdyzOr2sjZeM5kTMbGEwigFigAnROAZQCcssRCO2xT9HsA0QB6BdDaccdQ7X2hk7E7yN6ug4zdy9MZAci17F72XMvu5WmvNjI2nrHPtEAB6JwApHXKvUuXv+jdy1PrFzrXsfu5cy27n2nmGsmYjI0nEEAtUAA6JwAlt70zgOcAhPOGD6c/Oc6c6SFjp/EqBcjYNW7OliJnZ4k5n5+MnWfmbAkydpaYifNTAN45eeKJJ56/8kFYCSC/5aC7XOWWajnnJ17ALS1icLjFC7hqFl7AJl4ehnWdjA1Dm14xGRvPWFogZ+M5kzEZG08gAFugALxz0iVWlNwkIWJPkjCSfz8BYLPltVGWyP0iDv/MFAcwAJeR00MmY6eROV2AjJ1G5lIBcnYJm1OFyNgpXC5lJmOXsJm7EAWgueePvScBEiABEiABEiABpwlQADqNjAVIgARIgARIgARIwNwEKADNPX/sPQmQAAmQAAmQAAk4TYAC0GlkLEACJEACJEACJEAC5iZAAWju+WPvSYAESIAESIAESMBpAhSATiNjARIgARIgARIgARIwNwEKQHPPH3tPAiRAAiRAAiRAAk4ToAB0GhkLkAAJkAAJkAAJkIC5CVAAmnv+2HsSIAESIAESIAEScJoABaDTyFiABEggAAhUBHAEQAyAZQBecOOYJ1puEsoL4HGbG4bc2ASrIgESIAHHBCgAuUJIILAI/AagAYB4m2H/bREigUXC8WhFAP4LoDyA0zZZtwBYC+D9TMW/AZAIoLtGiEGW/BSAGoExGwmQgHsJUAC6lydrIwFfJ7ARgIiYkRo7KkIlWWNef8pmFYD3AjhFAehPU8uxkAAJCAEKQK4DEggsAo4EoHXb81UAbwKoBOAxAGIhfMXyWjmLIBIL2I826OT9oQBKAFgN4BKAKgCetuQRK9owAPMs/7dnAXsEwBgANQBEW/J+YBGg1vyvA3gRQE0AxwD0AbAtUz/eACDCLRbAdwDetuRZCmCcTd7nAUwGUAZASqZloEcAPghAOKda6pTv2TwARljGJy/TAhhYnzuOlgR8jgAFoM9NCTtEAoYS0CIANwPoDCASQAiArhZx1wbAHgANAfxqEXcRABpZtkWftfzdGsD3FktjU40C8D4AUlc3AD8DKAtgOYBvAUywEUwiRttbROg0AC0sQlWaeQ3AuwBE2MkYRHTVAfA7gJcAjAIgws6aNlj6KGUyJz0CMHNdMiY59ydi+pADAWzoxLNyEiABErAlQAHI9UACgUVABGB9ALcsOwBipeoPYL5FHInjQxMAclbQmv6xiLCvbF77EkAcgL4A5N+5LaLRmmWJRYBpFYAzAIhTxMs2bYjwHAKgmo0AFGG60JKnFoBdAIoBuAbgIIBPAXxsZ0pzWs7yiTOHnOET6+R+y5htz/hZizoSgA9Yxm7NK9+jMn6xNmY+A9jKImKbWQSutQwtgIH1ueNoScDnCFAA+tyUsEMkYCgBNQugOD6I+Dlp0wtxGJE/SZbX5HtDBMx6AGIVXANgp2WL11psCoDqALQKQKlDLIsiKq0pu2X7t2gWW6aZRZqUDbNsQduD+BGAChYL4iSLCBRrpb3kDgugbGmvtFgk5W/bRAFo6DJn5SRAAmoEKADVCPF9EvAvAloEYGbHh/8ADLJs69qjIRbAXAC62LyZ2QJ4wGKZ+8yS5x6LyLR6wcrZQDmzJ9u49pI9wZRZpMn2qlgS7VkApU4Rf2LNDAXwp2W7WbayjRCAsqUt29ADAYiHcOZEAehfnyuOhgRMR4AC0HRTxg6TgC4CagJQtoAl9Imt52s4gJ4W5ws5g5cDgGy/StiT3QAaA1hlsb6tAyDbnnIGUM7eWS2AIvBEWLa0bD2LUJPt3CcsQulhSx2yhfqLxfIngk3KyJZtVgLQtr/iECKOFlKveDpbzwDKv61JLHHFLX9knFZHjcxQ9VgA5fziHwCmW87+aRW0uiaWhUmABEjAGQIUgM7QYl4SMD8BcXwQYWYvDExWokdGLZ63Ys0SQSZbwfsADLfxwBUvYPm/nMcTL+CLli1WqwAsbTkrKDEIzwEYDOAnGwEobYgIHA2gtkXwyTa0CEWxMIoATMiU315/ewAQL2ARdzct5+/+ZzNt4qgi1kkRih86mE6rR3RmMSxWPRG5juIAihf155Yg0tbvWBGaUkacQSQFW8bDOIDm/0xxBCRgSgIUgKacNnaaBHyewFgA4ixhFYC+0mEJMSNWTAlnc95Bp0ToynaxnH0Ui6R4EbsrjQfQy2JJFecQEeRMJEACJOBRAhSAHsXNxkggYAj4ogCUreuvLVvX4mHMRAIkQAIBS4ACMGCnngMnAUMJ+JoAtMYmFKuebANLjEMmEiABEghYAhSAATv1HDgJkAAJkAAJkECgEqAADNSZ57hJgARIgARIgAQClgAFYMBOPQdOAiRAAiRAAiQQqAQoAAN15jluEiABEiABEiCBgCVAARiwU8+BkwAJkAAJkAAJBCoBCsBAnXmOmwRIgARIgARIIGAJUAAG7NRz4CRAAiRAAiRAAoFK4P9I1Ry+xJrAcQAAAABJRU5ErkJggg==\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig = response.plot(MIN_FREQ, UNIT, label='NRL STS-2')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 88,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "fig = response2.plot(MIN_FREQ, UNIT, label='Reference STS-2', show=False, axes=fig.axes )"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Compare"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 89,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "SR = 200.0\n",
    "t_samp = 1 / SR\n",
    "nfft = SR / MIN_FREQ"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 90,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "cpx1, freq1 = response.get_evalresp_response(t_samp = t_samp, nfft= nfft)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 91,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "cpx2, freq2 = response2.get_evalresp_response(t_samp = t_samp, nfft= nfft)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 92,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.0"
      ]
     },
     "execution_count": 92,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sum(freq1 - freq2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 93,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "abs1, abs2 = abs(cpx1),  abs(cpx2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 94,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "percent_diff = (abs1 - abs2 + .000001) / (abs1 + .0000001)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 95,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "diff = abs1 - abs2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 96,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x117051da0>]"
      ]
     },
     "execution_count": 96,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "fig.axes[0].loglog(freq1, diff, label='Difference')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 97,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/lloyd/anaconda2/envs/krisher/lib/python3.5/site-packages/matplotlib/axes/_axes.py:531: UserWarning: No labelled objects found. Use label='...' kwarg on individual plots.\n",
      "  warnings.warn(\"No labelled objects found. \"\n"
     ]
    }
   ],
   "source": [
    "# Fix broken legend\n",
    "fig.axes[1].legend(*fig.axes[0].get_legend_handles_labels())\n",
    "fig.axes[1].legend(loc='lower center')\n",
    "# Change line style so very similar response can be seen\n",
    "for i, line in enumerate(fig.axes[0].lines[::2]):\n",
    "    line.set_dashes([3, 6-i])\n",
    "for i, line in enumerate(fig.axes[1].lines[::2]):\n",
    "    line.set_dashes([3, 6-i])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "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.5.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}

RESP.RT130-L22-250SPS.TOM.3C.EHZ

#
###################################################################################
#
B050F03     Station:     TOM
B050F16     Network:     3C
B052F03     Location:    ??
B052F04     Channel:     EHZ
B052F22     Start date:  2014,189,00:00:00
B052F23     End date:    2014,218,23:59:59
#
#                  +-----------------------------------+
#                  |    Response (Poles and Zeros)     |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B053F03     Transfer function type:                A
B053F04     Stage sequence number:                 1
B053F05     Response in units lookup:              M/S - velocity in meters per second
B053F06     Response out units lookup:             V - emf in volts
B053F07     A0 normalization factor:               +1.00195E+00
B053F08     Normalization frequency:               +8.00000E+00
B053F09     Number of zeroes:                      2
B053F14     Number of poles:                       2
#              Complex zeroes:
#              i  real          imag          real_error    imag_error
B053F10-13     0  +0.00000E+00  +0.00000E+00  +0.00000E+00  +0.00000E+00
B053F10-13     1  +0.00000E+00  +0.00000E+00  +0.00000E+00  +0.00000E+00
#              Complex poles:
#              i  real          imag          real_error    imag_error
B053F15-18     0  -8.88600E+00  +8.88600E+00  +0.00000E+00  +0.00000E+00
B053F15-18     1  -8.88600E+00  -8.88600E+00  +0.00000E+00  +0.00000E+00
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 1
B058F04     Sensitivity:                           +8.78286E+01
B058F05     Frequency of sensitivity:              +8.00000E+00
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 2
B058F04     Sensitivity:                           +3.20000E+01
B058F05     Frequency of sensitivity:              +8.00000E+00
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 3
B054F05     Response in units lookup:              V - emf in volts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  0
B054F10     Number of denominators:                0
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  3
B057F04     Input sample rate (HZ):                 2.5600E+05
B057F05     Decimation factor:                      00001
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +0.0000E+00
B057F08     Correction applied (seconds):          +0.0000E+00
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 3
B058F04     Sensitivity:                           +6.29330E+05
B058F05     Frequency of sensitivity:              +8.00000E+00
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 4
B054F05     Response in units lookup:              COUNTS - digital counts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  29
B054F10     Number of denominators:                0
#              Numerator coefficients:
#              i  coefficient   error
B054F08-09     0  +2.44141E-04  +0.00000E+00
B054F08-09     1  +9.76562E-04  +0.00000E+00
B054F08-09     2  +2.44141E-03  +0.00000E+00
B054F08-09     3  +4.88281E-03  +0.00000E+00
B054F08-09     4  +8.54492E-03  +0.00000E+00
B054F08-09     5  +1.36719E-02  +0.00000E+00
B054F08-09     6  +2.05078E-02  +0.00000E+00
B054F08-09     7  +2.92969E-02  +0.00000E+00
B054F08-09     8  +3.93066E-02  +0.00000E+00
B054F08-09     9  +4.98047E-02  +0.00000E+00
B054F08-09    10  +6.00586E-02  +0.00000E+00
B054F08-09    11  +6.93359E-02  +0.00000E+00
B054F08-09    12  +7.69043E-02  +0.00000E+00
B054F08-09    13  +8.20312E-02  +0.00000E+00
B054F08-09    14  +8.39844E-02  +0.00000E+00
B054F08-09    15  +8.20312E-02  +0.00000E+00
B054F08-09    16  +7.69043E-02  +0.00000E+00
B054F08-09    17  +6.93359E-02  +0.00000E+00
B054F08-09    18  +6.00586E-02  +0.00000E+00
B054F08-09    19  +4.98047E-02  +0.00000E+00
B054F08-09    20  +3.93066E-02  +0.00000E+00
B054F08-09    21  +2.92969E-02  +0.00000E+00
B054F08-09    22  +2.05078E-02  +0.00000E+00
B054F08-09    23  +1.36719E-02  +0.00000E+00
B054F08-09    24  +8.54492E-03  +0.00000E+00
B054F08-09    25  +4.88281E-03  +0.00000E+00
B054F08-09    26  +2.44141E-03  +0.00000E+00
B054F08-09    27  +9.76562E-04  +0.00000E+00
B054F08-09    28  +2.44141E-04  +0.00000E+00
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  4
B057F04     Input sample rate (HZ):                 2.5600E+05
B057F05     Decimation factor:                      00008
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +5.4688E-05
B057F08     Correction applied (seconds):          +5.4688E-05
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 4
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +8.00000E+00
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 5
B054F05     Response in units lookup:              COUNTS - digital counts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  13
B054F10     Number of denominators:                0
#              Numerator coefficients:
#              i  coefficient   error
B054F08-09     0  +2.44141E-04  +0.00000E+00
B054F08-09     1  +2.92969E-03  +0.00000E+00
B054F08-09     2  +1.61133E-02  +0.00000E+00
B054F08-09     3  +5.37109E-02  +0.00000E+00
B054F08-09     4  +1.20850E-01  +0.00000E+00
B054F08-09     5  +1.93359E-01  +0.00000E+00
B054F08-09     6  +2.25586E-01  +0.00000E+00
B054F08-09     7  +1.93359E-01  +0.00000E+00
B054F08-09     8  +1.20850E-01  +0.00000E+00
B054F08-09     9  +5.37109E-02  +0.00000E+00
B054F08-09    10  +1.61133E-02  +0.00000E+00
B054F08-09    11  +2.92969E-03  +0.00000E+00
B054F08-09    12  +2.44141E-04  +0.00000E+00
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  5
B057F04     Input sample rate (HZ):                 3.2000E+04
B057F05     Decimation factor:                      00002
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +1.8750E-04
B057F08     Correction applied (seconds):          +1.8750E-04
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 5
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +8.00000E+00
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 6
B054F05     Response in units lookup:              COUNTS - digital counts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  13
B054F10     Number of denominators:                0
#              Numerator coefficients:
#              i  coefficient   error
B054F08-09     0  +2.44141E-04  +0.00000E+00
B054F08-09     1  +2.92969E-03  +0.00000E+00
B054F08-09     2  +1.61133E-02  +0.00000E+00
B054F08-09     3  +5.37109E-02  +0.00000E+00
B054F08-09     4  +1.20850E-01  +0.00000E+00
B054F08-09     5  +1.93359E-01  +0.00000E+00
B054F08-09     6  +2.25586E-01  +0.00000E+00
B054F08-09     7  +1.93359E-01  +0.00000E+00
B054F08-09     8  +1.20850E-01  +0.00000E+00
B054F08-09     9  +5.37109E-02  +0.00000E+00
B054F08-09    10  +1.61133E-02  +0.00000E+00
B054F08-09    11  +2.92969E-03  +0.00000E+00
B054F08-09    12  +2.44141E-04  +0.00000E+00
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  6
B057F04     Input sample rate (HZ):                 1.6000E+04
B057F05     Decimation factor:                      00002
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +3.7500E-04
B057F08     Correction applied (seconds):          +3.7500E-04
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 6
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +8.00000E+00
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 7
B054F05     Response in units lookup:              COUNTS - digital counts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  13
B054F10     Number of denominators:                0
#              Numerator coefficients:
#              i  coefficient   error
B054F08-09     0  +2.44141E-04  +0.00000E+00
B054F08-09     1  +2.92969E-03  +0.00000E+00
B054F08-09     2  +1.61133E-02  +0.00000E+00
B054F08-09     3  +5.37109E-02  +0.00000E+00
B054F08-09     4  +1.20850E-01  +0.00000E+00
B054F08-09     5  +1.93359E-01  +0.00000E+00
B054F08-09     6  +2.25586E-01  +0.00000E+00
B054F08-09     7  +1.93359E-01  +0.00000E+00
B054F08-09     8  +1.20850E-01  +0.00000E+00
B054F08-09     9  +5.37109E-02  +0.00000E+00
B054F08-09    10  +1.61133E-02  +0.00000E+00
B054F08-09    11  +2.92969E-03  +0.00000E+00
B054F08-09    12  +2.44141E-04  +0.00000E+00
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  7
B057F04     Input sample rate (HZ):                 8.0000E+03
B057F05     Decimation factor:                      00002
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +7.5000E-04
B057F08     Correction applied (seconds):          +7.5000E-04
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 7
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +8.00000E+00
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 8
B054F05     Response in units lookup:              COUNTS - digital counts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  13
B054F10     Number of denominators:                0
#              Numerator coefficients:
#              i  coefficient   error
B054F08-09     0  +2.44141E-04  +0.00000E+00
B054F08-09     1  +2.92969E-03  +0.00000E+00
B054F08-09     2  +1.61133E-02  +0.00000E+00
B054F08-09     3  +5.37109E-02  +0.00000E+00
B054F08-09     4  +1.20850E-01  +0.00000E+00
B054F08-09     5  +1.93359E-01  +0.00000E+00
B054F08-09     6  +2.25586E-01  +0.00000E+00
B054F08-09     7  +1.93359E-01  +0.00000E+00
B054F08-09     8  +1.20850E-01  +0.00000E+00
B054F08-09     9  +5.37109E-02  +0.00000E+00
B054F08-09    10  +1.61133E-02  +0.00000E+00
B054F08-09    11  +2.92969E-03  +0.00000E+00
B054F08-09    12  +2.44141E-04  +0.00000E+00
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  8
B057F04     Input sample rate (HZ):                 4.0000E+03
B057F05     Decimation factor:                      00002
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +1.5000E-03
B057F08     Correction applied (seconds):          +1.5000E-03
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 8
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +8.00000E+00
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 9
B054F05     Response in units lookup:              COUNTS - digital counts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  13
B054F10     Number of denominators:                0
#              Numerator coefficients:
#              i  coefficient   error
B054F08-09     0  +2.44141E-04  +0.00000E+00
B054F08-09     1  +2.92969E-03  +0.00000E+00
B054F08-09     2  +1.61133E-02  +0.00000E+00
B054F08-09     3  +5.37109E-02  +0.00000E+00
B054F08-09     4  +1.20850E-01  +0.00000E+00
B054F08-09     5  +1.93359E-01  +0.00000E+00
B054F08-09     6  +2.25586E-01  +0.00000E+00
B054F08-09     7  +1.93359E-01  +0.00000E+00
B054F08-09     8  +1.20850E-01  +0.00000E+00
B054F08-09     9  +5.37109E-02  +0.00000E+00
B054F08-09    10  +1.61133E-02  +0.00000E+00
B054F08-09    11  +2.92969E-03  +0.00000E+00
B054F08-09    12  +2.44141E-04  +0.00000E+00
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  9
B057F04     Input sample rate (HZ):                 2.0000E+03
B057F05     Decimation factor:                      00002
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +3.0000E-03
B057F08     Correction applied (seconds):          +3.0000E-03
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 9
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +8.00000E+00
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 10
B054F05     Response in units lookup:              COUNTS - digital counts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  13
B054F10     Number of denominators:                0
#              Numerator coefficients:
#              i  coefficient   error
B054F08-09     0  +2.44141E-04  +0.00000E+00
B054F08-09     1  +2.92969E-03  +0.00000E+00
B054F08-09     2  +1.61133E-02  +0.00000E+00
B054F08-09     3  +5.37109E-02  +0.00000E+00
B054F08-09     4  +1.20850E-01  +0.00000E+00
B054F08-09     5  +1.93359E-01  +0.00000E+00
B054F08-09     6  +2.25586E-01  +0.00000E+00
B054F08-09     7  +1.93359E-01  +0.00000E+00
B054F08-09     8  +1.20850E-01  +0.00000E+00
B054F08-09     9  +5.37109E-02  +0.00000E+00
B054F08-09    10  +1.61133E-02  +0.00000E+00
B054F08-09    11  +2.92969E-03  +0.00000E+00
B054F08-09    12  +2.44141E-04  +0.00000E+00
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  10
B057F04     Input sample rate (HZ):                 1.0000E+03
B057F05     Decimation factor:                      00002
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +6.0000E-03
B057F08     Correction applied (seconds):          +6.0000E-03
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 10
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +8.00000E+00
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 11
B054F05     Response in units lookup:              COUNTS - digital counts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  101
B054F10     Number of denominators:                0
#              Numerator coefficients:
#              i  coefficient   error
B054F08-09     0  -7.15032E-07  +0.00000E+00
B054F08-09     1  -5.60109E-06  +0.00000E+00
B054F08-09     2  -2.62179E-06  +0.00000E+00
B054F08-09     3  -4.31403E-05  +0.00000E+00
B054F08-09     4  -4.64771E-06  +0.00000E+00
B054F08-09     5  +1.43006E-06  +0.00000E+00
B054F08-09     6  +2.34769E-05  +0.00000E+00
B054F08-09     7  +1.43006E-06  +0.00000E+00
B054F08-09     8  -5.27932E-05  +0.00000E+00
B054F08-09     9  -3.66693E-04  +0.00000E+00
B054F08-09    10  +3.76107E-04  +0.00000E+00
B054F08-09    11  +8.54225E-04  +0.00000E+00
B054F08-09    12  +3.05081E-05  +0.00000E+00
B054F08-09    13  -1.27621E-03  +0.00000E+00
B054F08-09    14  -9.10951E-04  +0.00000E+00
B054F08-09    15  +1.27669E-03  +0.00000E+00
B054F08-09    16  +2.15165E-03  +0.00000E+00
B054F08-09    17  -4.61554E-04  +0.00000E+00
B054F08-09    18  -3.33765E-03  +0.00000E+00
B054F08-09    19  -1.40933E-03  +0.00000E+00
B054F08-09    20  +3.77072E-03  +0.00000E+00
B054F08-09    21  +4.19414E-03  +0.00000E+00
B054F08-09    22  -2.64288E-03  +0.00000E+00
B054F08-09    23  -7.20121E-03  +0.00000E+00
B054F08-09    24  -6.44006E-04  +0.00000E+00
B054F08-09    25  +9.18400E-03  +0.00000E+00
B054F08-09    26  +6.08445E-03  +0.00000E+00
B054F08-09    27  -8.57824E-03  +0.00000E+00
B054F08-09    28  -1.27401E-02  +0.00000E+00
B054F08-09    29  +3.98225E-03  +0.00000E+00
B054F08-09    30  +1.86261E-02  +0.00000E+00
B054F08-09    31  +5.20520E-03  +0.00000E+00
B054F08-09    32  -2.09407E-02  +0.00000E+00
B054F08-09    33  -1.81629E-02  +0.00000E+00
B054F08-09    34  +1.66669E-02  +0.00000E+00
B054F08-09    35  +3.22448E-02  +0.00000E+00
B054F08-09    36  -3.46588E-03  +0.00000E+00
B054F08-09    37  -4.29528E-02  +0.00000E+00
B054F08-09    38  -1.93265E-02  +0.00000E+00
B054F08-09    39  +4.43090E-02  +0.00000E+00
B054F08-09    40  +4.97909E-02  +0.00000E+00
B054F08-09    41  -2.94164E-02  +0.00000E+00
B054F08-09    42  -8.26078E-02  +0.00000E+00
B054F08-09    43  -9.34166E-03  +0.00000E+00
B054F08-09    44  +1.07552E-01  +0.00000E+00
B054F08-09    45  +8.16604E-02  +0.00000E+00
B054F08-09    46  -1.03110E-01  +0.00000E+00
B054F08-09    47  -2.04208E-01  +0.00000E+00
B054F08-09    48  -3.12231E-05  +0.00000E+00
B054F08-09    49  +3.90433E-01  +0.00000E+00
B054F08-09    50  +5.89958E-01  +0.00000E+00
B054F08-09    51  +3.90433E-01  +0.00000E+00
B054F08-09    52  -3.12231E-05  +0.00000E+00
B054F08-09    53  -2.04208E-01  +0.00000E+00
B054F08-09    54  -1.03110E-01  +0.00000E+00
B054F08-09    55  +8.16604E-02  +0.00000E+00
B054F08-09    56  +1.07552E-01  +0.00000E+00
B054F08-09    57  -9.34166E-03  +0.00000E+00
B054F08-09    58  -8.26078E-02  +0.00000E+00
B054F08-09    59  -2.94164E-02  +0.00000E+00
B054F08-09    60  +4.97909E-02  +0.00000E+00
B054F08-09    61  +4.43090E-02  +0.00000E+00
B054F08-09    62  -1.93265E-02  +0.00000E+00
B054F08-09    63  -4.29528E-02  +0.00000E+00
B054F08-09    64  -3.46588E-03  +0.00000E+00
B054F08-09    65  +3.22448E-02  +0.00000E+00
B054F08-09    66  +1.66669E-02  +0.00000E+00
B054F08-09    67  -1.81629E-02  +0.00000E+00
B054F08-09    68  -2.09407E-02  +0.00000E+00
B054F08-09    69  +5.20520E-03  +0.00000E+00
B054F08-09    70  +1.86261E-02  +0.00000E+00
B054F08-09    71  +3.98225E-03  +0.00000E+00
B054F08-09    72  -1.27401E-02  +0.00000E+00
B054F08-09    73  -8.57824E-03  +0.00000E+00
B054F08-09    74  +6.08445E-03  +0.00000E+00
B054F08-09    75  +9.18400E-03  +0.00000E+00
B054F08-09    76  -6.44006E-04  +0.00000E+00
B054F08-09    77  -7.20121E-03  +0.00000E+00
B054F08-09    78  -2.64288E-03  +0.00000E+00
B054F08-09    79  +4.19414E-03  +0.00000E+00
B054F08-09    80  +3.77072E-03  +0.00000E+00
B054F08-09    81  -1.40933E-03  +0.00000E+00
B054F08-09    82  -3.33765E-03  +0.00000E+00
B054F08-09    83  -4.61554E-04  +0.00000E+00
B054F08-09    84  +2.15165E-03  +0.00000E+00
B054F08-09    85  +1.27669E-03  +0.00000E+00
B054F08-09    86  -9.10951E-04  +0.00000E+00
B054F08-09    87  -1.27621E-03  +0.00000E+00
B054F08-09    88  +3.05081E-05  +0.00000E+00
B054F08-09    89  +8.54225E-04  +0.00000E+00
B054F08-09    90  +3.76107E-04  +0.00000E+00
B054F08-09    91  -3.66693E-04  +0.00000E+00
B054F08-09    92  -4.10309E-04  +0.00000E+00
B054F08-09    93  +2.52645E-05  +0.00000E+00
B054F08-09    94  +2.61821E-04  +0.00000E+00
B054F08-09    95  +1.20602E-04  +0.00000E+00
B054F08-09    96  -9.99854E-05  +0.00000E+00
B054F08-09    97  -1.62312E-04  +0.00000E+00
B054F08-09    98  -9.79595E-05  +0.00000E+00
B054F08-09    99  -2.94355E-05  +0.00000E+00
B054F08-09   100  -3.09847E-06  +0.00000E+00
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  11
B057F04     Input sample rate (HZ):                 5.0000E+02
B057F05     Decimation factor:                      00002
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +1.0000E-01
B057F08     Correction applied (seconds):          +1.0000E-01
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 11
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +8.00000E+00
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        3C  TOM        EHZ         |
#                  |     07/08/2014 to 08/06/2014      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 0
B058F04     Sensitivity:                           +1.76873E+09
B058F05     Frequency of sensitivity:              +8.00000E+00
B058F06     Number of calibrations:                0

RESP.RT130-STS2G2-200SPS.WHIH.4A.HHZ.txt

#
###################################################################################
#
B050F03     Station:     WHIH
B050F16     Network:     4A
B052F03     Location:    ??
B052F04     Channel:     HHZ
B052F22     Start date:  2008,315,00:00:00
B052F23     End date:    2008,365,23:59:59
#
#                  +-----------------------------------+
#                  |    Response (Poles and Zeros)     |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B053F03     Transfer function type:                A
B053F04     Stage sequence number:                 1
B053F05     Response in units lookup:              M/S - velocity in meters per second
B053F06     Response out units lookup:             V - emf in volts
B053F07     A0 normalization factor:               +2.35878E+17
B053F08     Normalization frequency:               +6.00000E-01
B053F09     Number of zeroes:                      9
B053F14     Number of poles:                       14
#              Complex zeroes:
#              i  real          imag          real_error    imag_error
B053F10-13     0  +0.00000E+00  +0.00000E+00  +0.00000E+00  +0.00000E+00
B053F10-13     1  +0.00000E+00  +0.00000E+00  +0.00000E+00  +0.00000E+00
B053F10-13     2  -5.90700E+03  +3.41100E+03  +0.00000E+00  +0.00000E+00
B053F10-13     3  -5.90700E+03  -3.41100E+03  +0.00000E+00  +0.00000E+00
B053F10-13     4  -6.83900E+02  +1.75500E+02  +0.00000E+00  +0.00000E+00
B053F10-13     5  -6.83900E+02  -1.75500E+02  +0.00000E+00  +0.00000E+00
B053F10-13     6  -5.55100E+02  +0.00000E+00  +0.00000E+00  +0.00000E+00
B053F10-13     7  -2.94600E+02  +0.00000E+00  +0.00000E+00  +0.00000E+00
B053F10-13     8  -1.07500E+01  +0.00000E+00  +0.00000E+00  +0.00000E+00
#              Complex poles:
#              i  real          imag          real_error    imag_error
B053F15-18     0  -6.90900E+03  +9.20800E+03  +0.00000E+00  +0.00000E+00
B053F15-18     1  -6.90900E+03  -9.20800E+03  +0.00000E+00  +0.00000E+00
B053F15-18     2  -6.22700E+03  +0.00000E+00  +0.00000E+00  +0.00000E+00
B053F15-18     3  -4.93600E+03  +4.71300E+03  +0.00000E+00  +0.00000E+00
B053F15-18     4  -4.93600E+03  -4.71300E+03  +0.00000E+00  +0.00000E+00
B053F15-18     5  -1.39100E+03  +0.00000E+00  +0.00000E+00  +0.00000E+00
B053F15-18     6  -5.56800E+02  +6.00500E+01  +0.00000E+00  +0.00000E+00
B053F15-18     7  -5.56800E+02  -6.00500E+01  +0.00000E+00  +0.00000E+00
B053F15-18     8  -9.84400E+01  +4.42800E+02  +0.00000E+00  +0.00000E+00
B053F15-18     9  -9.84400E+01  -4.42800E+02  +0.00000E+00  +0.00000E+00
B053F15-18    10  -1.09500E+01  +0.00000E+00  +0.00000E+00  +0.00000E+00
B053F15-18    11  -3.70000E-02  +3.70000E-02  +0.00000E+00  +0.00000E+00
B053F15-18    12  -3.70000E-02  -3.70000E-02  +0.00000E+00  +0.00000E+00
B053F15-18    13  -2.55100E+02  +0.00000E+00  +0.00000E+00  +0.00000E+00
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 1
B058F04     Sensitivity:                           +1.49257E+03
B058F05     Frequency of sensitivity:              +6.00000E-01
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 2
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +6.00000E-01
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 3
B054F05     Response in units lookup:              V - emf in volts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  0
B054F10     Number of denominators:                0
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  3
B057F04     Input sample rate (HZ):                 1.0240E+05
B057F05     Decimation factor:                      00001
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +0.0000E+00
B057F08     Correction applied (seconds):          +0.0000E+00
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 3
B058F04     Sensitivity:                           +6.29330E+05
B058F05     Frequency of sensitivity:              +6.00000E-01
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 4
B054F05     Response in units lookup:              COUNTS - digital counts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  29
B054F10     Number of denominators:                0
#              Numerator coefficients:
#              i  coefficient   error
B054F08-09     0  +2.44141E-04  +0.00000E+00
B054F08-09     1  +9.76562E-04  +0.00000E+00
B054F08-09     2  +2.44141E-03  +0.00000E+00
B054F08-09     3  +4.88281E-03  +0.00000E+00
B054F08-09     4  +8.54492E-03  +0.00000E+00
B054F08-09     5  +1.36719E-02  +0.00000E+00
B054F08-09     6  +2.05078E-02  +0.00000E+00
B054F08-09     7  +2.92969E-02  +0.00000E+00
B054F08-09     8  +3.93066E-02  +0.00000E+00
B054F08-09     9  +4.98047E-02  +0.00000E+00
B054F08-09    10  +6.00586E-02  +0.00000E+00
B054F08-09    11  +6.93359E-02  +0.00000E+00
B054F08-09    12  +7.69043E-02  +0.00000E+00
B054F08-09    13  +8.20312E-02  +0.00000E+00
B054F08-09    14  +8.39844E-02  +0.00000E+00
B054F08-09    15  +8.20312E-02  +0.00000E+00
B054F08-09    16  +7.69043E-02  +0.00000E+00
B054F08-09    17  +6.93359E-02  +0.00000E+00
B054F08-09    18  +6.00586E-02  +0.00000E+00
B054F08-09    19  +4.98047E-02  +0.00000E+00
B054F08-09    20  +3.93066E-02  +0.00000E+00
B054F08-09    21  +2.92969E-02  +0.00000E+00
B054F08-09    22  +2.05078E-02  +0.00000E+00
B054F08-09    23  +1.36719E-02  +0.00000E+00
B054F08-09    24  +8.54492E-03  +0.00000E+00
B054F08-09    25  +4.88281E-03  +0.00000E+00
B054F08-09    26  +2.44141E-03  +0.00000E+00
B054F08-09    27  +9.76562E-04  +0.00000E+00
B054F08-09    28  +2.44141E-04  +0.00000E+00
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  4
B057F04     Input sample rate (HZ):                 1.0240E+05
B057F05     Decimation factor:                      00008
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +1.3672E-04
B057F08     Correction applied (seconds):          +1.3672E-04
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 4
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +6.00000E-01
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 5
B054F05     Response in units lookup:              COUNTS - digital counts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  13
B054F10     Number of denominators:                0
#              Numerator coefficients:
#              i  coefficient   error
B054F08-09     0  +2.44141E-04  +0.00000E+00
B054F08-09     1  +2.92969E-03  +0.00000E+00
B054F08-09     2  +1.61133E-02  +0.00000E+00
B054F08-09     3  +5.37109E-02  +0.00000E+00
B054F08-09     4  +1.20850E-01  +0.00000E+00
B054F08-09     5  +1.93359E-01  +0.00000E+00
B054F08-09     6  +2.25586E-01  +0.00000E+00
B054F08-09     7  +1.93359E-01  +0.00000E+00
B054F08-09     8  +1.20850E-01  +0.00000E+00
B054F08-09     9  +5.37109E-02  +0.00000E+00
B054F08-09    10  +1.61133E-02  +0.00000E+00
B054F08-09    11  +2.92969E-03  +0.00000E+00
B054F08-09    12  +2.44141E-04  +0.00000E+00
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  5
B057F04     Input sample rate (HZ):                 1.2800E+04
B057F05     Decimation factor:                      00002
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +4.6875E-04
B057F08     Correction applied (seconds):          +4.6875E-04
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 5
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +6.00000E-01
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 6
B054F05     Response in units lookup:              COUNTS - digital counts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  13
B054F10     Number of denominators:                0
#              Numerator coefficients:
#              i  coefficient   error
B054F08-09     0  +2.44141E-04  +0.00000E+00
B054F08-09     1  +2.92969E-03  +0.00000E+00
B054F08-09     2  +1.61133E-02  +0.00000E+00
B054F08-09     3  +5.37109E-02  +0.00000E+00
B054F08-09     4  +1.20850E-01  +0.00000E+00
B054F08-09     5  +1.93359E-01  +0.00000E+00
B054F08-09     6  +2.25586E-01  +0.00000E+00
B054F08-09     7  +1.93359E-01  +0.00000E+00
B054F08-09     8  +1.20850E-01  +0.00000E+00
B054F08-09     9  +5.37109E-02  +0.00000E+00
B054F08-09    10  +1.61133E-02  +0.00000E+00
B054F08-09    11  +2.92969E-03  +0.00000E+00
B054F08-09    12  +2.44141E-04  +0.00000E+00
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  6
B057F04     Input sample rate (HZ):                 6.4000E+03
B057F05     Decimation factor:                      00002
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +9.3750E-04
B057F08     Correction applied (seconds):          +9.3750E-04
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 6
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +6.00000E-01
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 7
B054F05     Response in units lookup:              COUNTS - digital counts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  13
B054F10     Number of denominators:                0
#              Numerator coefficients:
#              i  coefficient   error
B054F08-09     0  +2.44141E-04  +0.00000E+00
B054F08-09     1  +2.92969E-03  +0.00000E+00
B054F08-09     2  +1.61133E-02  +0.00000E+00
B054F08-09     3  +5.37109E-02  +0.00000E+00
B054F08-09     4  +1.20850E-01  +0.00000E+00
B054F08-09     5  +1.93359E-01  +0.00000E+00
B054F08-09     6  +2.25586E-01  +0.00000E+00
B054F08-09     7  +1.93359E-01  +0.00000E+00
B054F08-09     8  +1.20850E-01  +0.00000E+00
B054F08-09     9  +5.37109E-02  +0.00000E+00
B054F08-09    10  +1.61133E-02  +0.00000E+00
B054F08-09    11  +2.92969E-03  +0.00000E+00
B054F08-09    12  +2.44141E-04  +0.00000E+00
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  7
B057F04     Input sample rate (HZ):                 3.2000E+03
B057F05     Decimation factor:                      00002
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +1.8750E-03
B057F08     Correction applied (seconds):          +1.8750E-03
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 7
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +6.00000E-01
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 8
B054F05     Response in units lookup:              COUNTS - digital counts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  13
B054F10     Number of denominators:                0
#              Numerator coefficients:
#              i  coefficient   error
B054F08-09     0  +2.44141E-04  +0.00000E+00
B054F08-09     1  +2.92969E-03  +0.00000E+00
B054F08-09     2  +1.61133E-02  +0.00000E+00
B054F08-09     3  +5.37109E-02  +0.00000E+00
B054F08-09     4  +1.20850E-01  +0.00000E+00
B054F08-09     5  +1.93359E-01  +0.00000E+00
B054F08-09     6  +2.25586E-01  +0.00000E+00
B054F08-09     7  +1.93359E-01  +0.00000E+00
B054F08-09     8  +1.20850E-01  +0.00000E+00
B054F08-09     9  +5.37109E-02  +0.00000E+00
B054F08-09    10  +1.61133E-02  +0.00000E+00
B054F08-09    11  +2.92969E-03  +0.00000E+00
B054F08-09    12  +2.44141E-04  +0.00000E+00
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  8
B057F04     Input sample rate (HZ):                 1.6000E+03
B057F05     Decimation factor:                      00002
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +3.7500E-03
B057F08     Correction applied (seconds):          +3.7500E-03
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 8
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +6.00000E-01
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 9
B054F05     Response in units lookup:              COUNTS - digital counts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  13
B054F10     Number of denominators:                0
#              Numerator coefficients:
#              i  coefficient   error
B054F08-09     0  +2.44141E-04  +0.00000E+00
B054F08-09     1  +2.92969E-03  +0.00000E+00
B054F08-09     2  +1.61133E-02  +0.00000E+00
B054F08-09     3  +5.37109E-02  +0.00000E+00
B054F08-09     4  +1.20850E-01  +0.00000E+00
B054F08-09     5  +1.93359E-01  +0.00000E+00
B054F08-09     6  +2.25586E-01  +0.00000E+00
B054F08-09     7  +1.93359E-01  +0.00000E+00
B054F08-09     8  +1.20850E-01  +0.00000E+00
B054F08-09     9  +5.37109E-02  +0.00000E+00
B054F08-09    10  +1.61133E-02  +0.00000E+00
B054F08-09    11  +2.92969E-03  +0.00000E+00
B054F08-09    12  +2.44141E-04  +0.00000E+00
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  9
B057F04     Input sample rate (HZ):                 8.0000E+02
B057F05     Decimation factor:                      00002
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +7.5000E-03
B057F08     Correction applied (seconds):          +7.5000E-03
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 9
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +6.00000E-01
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Response (Coefficients)      |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B054F03     Transfer function type:                D
B054F04     Stage sequence number:                 10
B054F05     Response in units lookup:              COUNTS - digital counts
B054F06     Response out units lookup:             COUNTS - digital counts
B054F07     Number of numerators:                  101
B054F10     Number of denominators:                0
#              Numerator coefficients:
#              i  coefficient   error
B054F08-09     0  -7.15032E-07  +0.00000E+00
B054F08-09     1  -5.60109E-06  +0.00000E+00
B054F08-09     2  -2.62179E-06  +0.00000E+00
B054F08-09     3  -4.31403E-05  +0.00000E+00
B054F08-09     4  -4.64771E-06  +0.00000E+00
B054F08-09     5  +1.43006E-06  +0.00000E+00
B054F08-09     6  +2.34769E-05  +0.00000E+00
B054F08-09     7  +1.43006E-06  +0.00000E+00
B054F08-09     8  -5.27932E-05  +0.00000E+00
B054F08-09     9  -3.66693E-04  +0.00000E+00
B054F08-09    10  +3.76107E-04  +0.00000E+00
B054F08-09    11  +8.54225E-04  +0.00000E+00
B054F08-09    12  +3.05081E-05  +0.00000E+00
B054F08-09    13  -1.27621E-03  +0.00000E+00
B054F08-09    14  -9.10951E-04  +0.00000E+00
B054F08-09    15  +1.27669E-03  +0.00000E+00
B054F08-09    16  +2.15165E-03  +0.00000E+00
B054F08-09    17  -4.61554E-04  +0.00000E+00
B054F08-09    18  -3.33765E-03  +0.00000E+00
B054F08-09    19  -1.40933E-03  +0.00000E+00
B054F08-09    20  +3.77072E-03  +0.00000E+00
B054F08-09    21  +4.19414E-03  +0.00000E+00
B054F08-09    22  -2.64288E-03  +0.00000E+00
B054F08-09    23  -7.20121E-03  +0.00000E+00
B054F08-09    24  -6.44006E-04  +0.00000E+00
B054F08-09    25  +9.18400E-03  +0.00000E+00
B054F08-09    26  +6.08445E-03  +0.00000E+00
B054F08-09    27  -8.57824E-03  +0.00000E+00
B054F08-09    28  -1.27401E-02  +0.00000E+00
B054F08-09    29  +3.98225E-03  +0.00000E+00
B054F08-09    30  +1.86261E-02  +0.00000E+00
B054F08-09    31  +5.20520E-03  +0.00000E+00
B054F08-09    32  -2.09407E-02  +0.00000E+00
B054F08-09    33  -1.81629E-02  +0.00000E+00
B054F08-09    34  +1.66669E-02  +0.00000E+00
B054F08-09    35  +3.22448E-02  +0.00000E+00
B054F08-09    36  -3.46588E-03  +0.00000E+00
B054F08-09    37  -4.29528E-02  +0.00000E+00
B054F08-09    38  -1.93265E-02  +0.00000E+00
B054F08-09    39  +4.43090E-02  +0.00000E+00
B054F08-09    40  +4.97909E-02  +0.00000E+00
B054F08-09    41  -2.94164E-02  +0.00000E+00
B054F08-09    42  -8.26078E-02  +0.00000E+00
B054F08-09    43  -9.34166E-03  +0.00000E+00
B054F08-09    44  +1.07552E-01  +0.00000E+00
B054F08-09    45  +8.16604E-02  +0.00000E+00
B054F08-09    46  -1.03110E-01  +0.00000E+00
B054F08-09    47  -2.04208E-01  +0.00000E+00
B054F08-09    48  -3.12231E-05  +0.00000E+00
B054F08-09    49  +3.90433E-01  +0.00000E+00
B054F08-09    50  +5.89958E-01  +0.00000E+00
B054F08-09    51  +3.90433E-01  +0.00000E+00
B054F08-09    52  -3.12231E-05  +0.00000E+00
B054F08-09    53  -2.04208E-01  +0.00000E+00
B054F08-09    54  -1.03110E-01  +0.00000E+00
B054F08-09    55  +8.16604E-02  +0.00000E+00
B054F08-09    56  +1.07552E-01  +0.00000E+00
B054F08-09    57  -9.34166E-03  +0.00000E+00
B054F08-09    58  -8.26078E-02  +0.00000E+00
B054F08-09    59  -2.94164E-02  +0.00000E+00
B054F08-09    60  +4.97909E-02  +0.00000E+00
B054F08-09    61  +4.43090E-02  +0.00000E+00
B054F08-09    62  -1.93265E-02  +0.00000E+00
B054F08-09    63  -4.29528E-02  +0.00000E+00
B054F08-09    64  -3.46588E-03  +0.00000E+00
B054F08-09    65  +3.22448E-02  +0.00000E+00
B054F08-09    66  +1.66669E-02  +0.00000E+00
B054F08-09    67  -1.81629E-02  +0.00000E+00
B054F08-09    68  -2.09407E-02  +0.00000E+00
B054F08-09    69  +5.20520E-03  +0.00000E+00
B054F08-09    70  +1.86261E-02  +0.00000E+00
B054F08-09    71  +3.98225E-03  +0.00000E+00
B054F08-09    72  -1.27401E-02  +0.00000E+00
B054F08-09    73  -8.57824E-03  +0.00000E+00
B054F08-09    74  +6.08445E-03  +0.00000E+00
B054F08-09    75  +9.18400E-03  +0.00000E+00
B054F08-09    76  -6.44006E-04  +0.00000E+00
B054F08-09    77  -7.20121E-03  +0.00000E+00
B054F08-09    78  -2.64288E-03  +0.00000E+00
B054F08-09    79  +4.19414E-03  +0.00000E+00
B054F08-09    80  +3.77072E-03  +0.00000E+00
B054F08-09    81  -1.40933E-03  +0.00000E+00
B054F08-09    82  -3.33765E-03  +0.00000E+00
B054F08-09    83  -4.61554E-04  +0.00000E+00
B054F08-09    84  +2.15165E-03  +0.00000E+00
B054F08-09    85  +1.27669E-03  +0.00000E+00
B054F08-09    86  -9.10951E-04  +0.00000E+00
B054F08-09    87  -1.27621E-03  +0.00000E+00
B054F08-09    88  +3.05081E-05  +0.00000E+00
B054F08-09    89  +8.54225E-04  +0.00000E+00
B054F08-09    90  +3.76107E-04  +0.00000E+00
B054F08-09    91  -3.66693E-04  +0.00000E+00
B054F08-09    92  -4.10309E-04  +0.00000E+00
B054F08-09    93  +2.52645E-05  +0.00000E+00
B054F08-09    94  +2.61821E-04  +0.00000E+00
B054F08-09    95  +1.20602E-04  +0.00000E+00
B054F08-09    96  -9.99854E-05  +0.00000E+00
B054F08-09    97  -1.62312E-04  +0.00000E+00
B054F08-09    98  -9.79595E-05  +0.00000E+00
B054F08-09    99  -2.94355E-05  +0.00000E+00
B054F08-09   100  -3.09847E-06  +0.00000E+00
#
#                  +-----------------------------------+
#                  |            Decimation             |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B057F03     Stage sequence number:                  10
B057F04     Input sample rate (HZ):                 4.0000E+02
B057F05     Decimation factor:                      00002
B057F06     Decimation offset:                      00000
B057F07     Estimated delay (seconds):             +1.2500E-01
B057F08     Correction applied (seconds):          +1.2500E-01
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 10
B058F04     Sensitivity:                           +1.00000E+00
B058F05     Frequency of sensitivity:              +6.00000E-01
B058F06     Number of calibrations:                0
#
#                  +-----------------------------------+
#                  |      Channel Sensitivity/Gain     |
#                  |        4A  WHIH       HHZ         |
#                  |     11/10/2008 to 12/30/2008      |
#                  +-----------------------------------+
#
B058F03     Stage sequence number:                 0
B058F04     Sensitivity:                           +9.39316E+08
B058F05     Frequency of sensitivity:              +6.00000E-01
B058F06     Number of calibrations:                0