docs/notebooks/xspec_models.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Working with XSPEC models\n",
    "\n",
    "One of the most powerful aspects of **XSPEC** is a huge modeling community. While in 3ML, we are focused on building a powerful and modular data analysis tool, we cannot neglect the need for many of the models thahat already exist in **XSPEC** and thus provide support for them via **astromodels** directly in 3ML. \n",
    "\n",
    "For details on installing **astromodels** with **XSPEC** support, visit the 3ML or **astromodels** installation page. \n",
    "\n",
    "\n",
    "Let's explore how we can use **XSPEC** spectral models in 3ML. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-01-04T13:10:51.894290Z",
     "start_time": "2018-01-04T13:10:51.883115Z"
    }
   },
   "outputs": [],
   "source": [
    "%matplotlib notebook\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We do not load the models by default as this takes some time and 3ML should load quickly. However, if you need the **XSPEC** models, they are imported from astromodels like this:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-01-04T13:06:22.923013Z",
     "start_time": "2018-01-04T13:06:16.023150Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Loading xspec models...done\n"
     ]
    }
   ],
   "source": [
    "from astromodels.xspec.factory import *"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The models are indexed with *XS_* before the typical **XSPEC** model names."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-01-04T13:07:08.921540Z",
     "start_time": "2018-01-04T13:07:08.891013Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<ul>\n",
       "\n",
       "<li>description: The phabs model from XSpec (https://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/XspecModels.html)</li>\n",
       "\n",
       "<li>formula: $n.a.$</li>\n",
       "\n",
       "<li>parameters: \n",
       "<ul>\n",
       "\n",
       "<li>nh: \n",
       "<ul>\n",
       "\n",
       "<li>value: 1.0</li>\n",
       "\n",
       "<li>desc: (see https://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/XspecModels.html)</li>\n",
       "\n",
       "<li>min_value: 0.0</li>\n",
       "\n",
       "<li>max_value: 1000000.0</li>\n",
       "\n",
       "<li>unit: 10**22</li>\n",
       "\n",
       "<li>is_normalization: False</li>\n",
       "\n",
       "<li>delta: 0.001</li>\n",
       "\n",
       "<li>free: True</li>\n",
       "\n",
       "</ul>\n",
       "\n",
       "</li>\n",
       "\n",
       "</ul>\n",
       "\n",
       "</li>\n",
       "\n",
       "</ul>\n"
      ],
      "text/plain": [
       "  * description: The phabs model from XSpec (https://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/XspecModels.html)\n",
       "  * formula: $n.a.$\n",
       "  * parameters: \n",
       "    * nh: \n",
       "      * value: 1.0\n",
       "      * desc: (see https://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/XspecModels.html)\n",
       "      * min_value: 0.0\n",
       "      * max_value: 1000000.0\n",
       "      * unit: 10**22\n",
       "      * is_normalization: False\n",
       "      * delta: 0.001\n",
       "      * free: True\n"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "plaw = XS_powerlaw()\n",
    "phabs = XS_phabs()\n",
    "phabs\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The spectral models behave just as any other **astromodels** spectral model and can be used in combination with other **astromodels** spectral models."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-01-04T13:11:29.378937Z",
     "start_time": "2018-01-04T13:11:29.285614Z"
    }
   },
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "\n",
       "mpl.get_websocket_type = function() {\n",
       "    if (typeof(WebSocket) !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert('Your browser does not have WebSocket support.' +\n",
       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "              'Firefox 4 and 5 are also supported but you ' +\n",
       "              'have to enable WebSockets in about:config.');\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent = (\n",
       "                \"This browser does not support binary websocket messages. \" +\n",
       "                    \"Performance may be slow.\");\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen =  function () {\n",
       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
       "            fig.send_message(\"send_image_mode\", {});\n",
       "            if (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
       "            }\n",
       "            fig.send_message(\"refresh\", {});\n",
       "        }\n",
       "\n",
       "    this.imageObj.onload = function() {\n",
       "            if (fig.image_mode == 'full') {\n",
       "                // Full images could contain transparency (where diff images\n",
       "                // almost always do), so we need to clear the canvas so that\n",
       "                // there is no ghosting.\n",
       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "            }\n",
       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "        };\n",
       "\n",
       "    this.imageObj.onunload = function() {\n",
       "        fig.ws.close();\n",
       "    }\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_header = function() {\n",
       "    var titlebar = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus () {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option)\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_message = function(type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function() {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
       "    }\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0, 0, fig.canvas.width, fig.canvas.height);\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch(cursor)\n",
       "    {\n",
       "    case 0:\n",
       "        cursor = 'pointer';\n",
       "        break;\n",
       "    case 1:\n",
       "        cursor = 'default';\n",
       "        break;\n",
       "    case 2:\n",
       "        cursor = 'crosshair';\n",
       "        break;\n",
       "    case 3:\n",
       "        cursor = 'move';\n",
       "        break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message(\"ack\", {});\n",
       "}\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = \"image/png\";\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src);\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig[\"handle_\" + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "}\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function(e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e)\n",
       "        e = window.event;\n",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\n",
       "\n",
       "    return {\"x\": x, \"y\": y};\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys (original) {\n",
       "  return Object.keys(original).reduce(function (obj, key) {\n",
       "    if (typeof original[key] !== 'object')\n",
       "        obj[key] = original[key]\n",
       "    return obj;\n",
       "  }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
       "    var canvas_pos = mpl.findpos(event)\n",
       "\n",
       "    if (name === 'button_press')\n",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.key_event = function(event, name) {\n",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, {key: value,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
       "    if (name == 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message(\"toolbar_button\", {name: name});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to  previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"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",
       "    var width = fig.canvas.width/mpl.ratio\n",
       "    fig.root.unbind('remove')\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable()\n",
       "    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n",
       "    fig.close_ws(fig, msg);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.close_ws = function(fig, msg){\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width/mpl.ratio\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message(\"ack\", {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () { fig.push_to_output() }, 1000);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items){\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) { continue; };\n",
       "\n",
       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
       "    buttongrp.append(button);\n",
       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
       "    titlebar.prepend(buttongrp);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(el){\n",
       "    var fig = this\n",
       "    el.on(\"remove\", function(){\n",
       "\tfig.close_ws(fig, {});\n",
       "    });\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
       "    // this is important to make the div 'focusable\n",
       "    el.attr('tabindex', 0)\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    }\n",
       "    else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager)\n",
       "        manager = IPython.keyboard_manager;\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which == 13) {\n",
       "        this.canvas_div.blur();\n",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.find_output_cell = function(html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i=0; i<ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code'){\n",
       "            for (var j=0; j<cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] == html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "}\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel != null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAABQAAAAPACAYAAABq3NR5AAAAAXNSR0IArs4c6QAAQABJREFUeAHs3QeYXVd9KPo16r33YsuyLHfLVe6YXmIMhIROsA3G3HzJC3lweUDghhJKLiWEBHIfAYIxhICBBNMMDs02xtiWC7bkIsm2eu+S1Udz1zpzzszZZ0ZTz8yc8lvft7T3Wrut/duDxfy1SggSAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAgToSaKijd/WqlSUwPDbn7HyTtsZtY2U1T2sIECBAgAABAgQIECBAgEBNCAyObzE1/yaPxu2hmngrL9EtgSHdOtvJBMonkIJ/95fvdu5EgAABAgQIECBAgAABAgQIdCJwUTy+pJNzHK5BgUE1+E5eiQABAgQIECBAgAABAgQIECBAgACBvIAegH4UBkogDfvNpfvuuy/MnDmzULQlQIAAAQIECBAgQIAAAQIEyiSwcePGsHjx4sLdWn4XL1TY1oeAAGB9fOdKfMuWOf9S8G/OnDmV2EZtIkCAAAECBAgQIECAAAECtSTQ8rt4Lb2Ud+lcwBDgzo2cQYAAAQIECBAgQIAAAQIECBAgQKBqBQQAq/bTaTgBAgQIECBAgAABAgQIECBAgACBzgUEADs3cgYBAgQIECBAgAABAgQIECBAgACBqhUQAKzaT6fhBAgQIECAAAECBAgQIECAAAECBDoXEADs3MgZBAgQIECAAAECBAgQIECAAAECBKpWQACwaj+dhhMgQIAAAQIECBAgQIAAAQIECBDoXEAAsHMjZxAgQIAAAQIECBAgQIAAAQIECBCoWgEBwKr9dBpOgAABAgQIECBAgAABAgQIECBAoHMBAcDOjZxBgAABAgQIECBAgAABAgQIECBAoGoFBACr9tNpOAECBAgQIECAAAECBAgQIECAAIHOBQQAOzdyBgECBAgQIECAAAECBAgQIECAAIGqFRAArNpPp+EECBAgQIAAAQIECBAgQIAAAQIEOhcQAOzcyBkECBAgQIAAAQIECBAgQIAAAQIEqlZAALBqP52GEyBAgAABAgQIECBAgAABAgQIEOhcQACwcyNnECBAgAABAgQIECBAgAABAgQIEKhaAQHAqv10Gk6AAAECBAgQIECAAAECBAgQIECgcwEBwM6NnEGAAAECBAgQIECAAAECBAgQIECgagUEAKv202k4AQIECBAgQIAAAQIECBAgQIAAgc4FBAA7N3IGAQIECBAgQIAAAQIECBAgQIAAgaoVEACs2k+n4QQIECBAgAABAgQIECBAgAABAgQ6FxAA7NzIGQQIECBAgAABAgQIECBAgAABAgSqVkAAsGo/nYYTIECAAAECBAgQIECAAAECBAgQ6FxAALBzI2cQIECAAAECBAgQIECAAAECBAgQqFoBAcCq/XQaToAAAQIECBAgQIAAAQIECBAgQKBzAQHAzo2cQYAAAQIECBAgQIAAAQIECBAgQKBqBQQAq/bTaTgBAgQIECBAgAABAgQIECBAgACBzgUEADs3cgYBAgQIECBAgAABAgQIECBAgACBqhUQAKzaT6fhBAgQIECAAAECBAgQIECAAAECBDoXGNL5Kc4gQIAAAQIECBAgQKA7AseONYXDjcfCoSMxH20MjU1N4WhjUzgWt43xWNoejdvc/rGQO572U26Kx5ryD4u7udRSE8vHPRbPHNzQEAYNagiDYx4U95u3oWW/uL753HhNOj9/7sRRw3LX5x9vQ4AAAQIECNSIgABgjXxIr0GAAAECBAgQINAqkIJuu/cfCbsOHAk7nz2c26bynoNHwt6DR/P5SBgyeFBYfNLEcNKUMTFY1xiePXw0PHsobg/F7eHC9mjY31IX92P9oaMpsJcCfI3hcGE/PrM54HcsF/xrbU317I0fOTQsmjshjB42OGczdHBDGDpoUBg6pCEMidthQwbFbayLbrljcZsMh8Xz0ra4flgsDx8a85DBYcRxtsPj/VLAUiJAgAABAgT6VkAAsG993Z0AAQIECBAgQKAMAqlX3J4YuNu0+2DYtOdg2FzYxv2d+2OALwX7crk52JeCdF1N/3Hfmq6eWvPn7Y4B0zuXb+3X98wFCmMgcPjQwTFYOKglWJiChyNi8LCwTUHEkcOGhJHxvFExQDky5tw2lpv3h+TKI/LHC+c0nz8k19OxX1/MwwgQIECAQAUJCABW0MfQFAIECBAgQIBAPQocjUNlt+47lAvubY4BveYgXyofaA727Wk+diD2tpNqTyANlU55b+x12Zcp9V4sBAynjh0eZo0fGUYNHxzGDE+BwyFxmwKKaRvLcX903I6O5dFpP7dt3k8BxYY4ZFoiQIAAAQLVJCAAWE1fS1sJECBAgAABAlUmkObC2xADeWu2788F81LvvVyAL/bgywX7Ynnr3kNxTrwqe7EyNTc3/17RHHxpNGyqS6kQZCqEmlpjTg3xWHMDio+l+QILcww2zzMY5xZM8wrGA+k7pG06p15TGqqd8q5wJGyMP3+PrNvdI4pkXxoYTIHF5sBhayAxBRBTUHHM8KFhzIghYWzKsW7siOZyOj/lwvfuUWNcRIAAAQIEuiggANhFKKcRIECAAAECBAi0L5CG526P8+w9s+3Z8MzWZ8PTcbsq7aft9mdzc+W1f2X/1KbAS1rcIs1vl3Iqp8BLGjaa2njfMztygbLmHl/NPb9G5YIzRT3CCgGe2BMsBXsKc9o1z3HXPM9dGr5amPMut58fvprmxRsS58grLLSRFudI8+gNxNx36Vu1BgVbFx9JAcJtsRfmvdFiVxxSfSQuWHIk9spLC5WkbW4/1qWeemkxk+a6wjnx+NG4fyx7XmFuxIOx52aaLzGVayGlIOq+2Fsx5RAO9fqV0nyLhaBg4WdzXAoSpp/BfOAw7efqisq5gGL+vPQzWQgY97pBbkCAAAECNSkgAFiTn9VLESBAgAABAgT6RmB7DBI9uWlveCLm5Zubt09t3ZdbVKNvntj2rpNGDwvTx40IM8YNDzPGjwhTxgzPBfYmxCDfxFFDw4RcHhYm5AN+aXGKjlIKitVL8CS9ZwpGtvdLwMToesr0sR1R9epY8crIB1sWTGkMB/MrJRe2KVhYCBoWb9PxA3GRljS/4/4YVDyYtvn9VJ+GiB/I16X9FKirhpRbbCa2O+zpeWtTp9EUJEyBxBQYHBd/9lPAcFzLfvOxcSObA4npvNb95mNpiLREgAABArUr0N7f/bX7tt6MAAECBAgQIECgSwJpXr6nYm++R9fvDss27M4F+1Lgb9u+w126vicnpVVlmwN7I8L0GNibkQvyNe/PzJenxaBfWlW2nKlegn/lNOvJvVKPxxGDUu/JwWF8GNqTW3T5mhTUTYHEXIAwBgdTILF5v3ll59VxSHpa8GR/PLYvrvCctoXVn5vrmgONqZdfWhE6doSs6JTalxbJSbmnKc1t2Bo8bA4i5gKF+SBi7lgucJgPMsb98TGg2HzO0FyvWP9b6qm+6wgQIND3AgKAfW/sCQQIECBAgACBihZIQ0JXbtmXC/YtjQG/R9btCo9t3JPrmVWuhqfeSbMmxGBePqiXAnqFIF+uLpYnxR58AzEstlzv6D6VI5ACUSnQmHLqMdqbVAgmpkBgLkiYCxbG/dhrr7muOUhYKKdAYwocpkDi3hiQS+ekbaorbHvTnr66NteDMgZKt8Q5OXuSUgC/3YBhDBSm4GEafp96JhaG4hfvpzo9EHui7hoCBAh0XUAAsOtWziRAgAABAgQI1IRAmuPtwTU7w5JVzTn18ivHCrvD4lDbeVNGhXmTR4eTpo4O86fE7ZQxMY+Ow3SH1c0w25r4IfESLQLFwcTJY1qqe7yThkLvSz0PY1CwOSB4JNdzL5Wbg4RHcsdSb76WciF4GOua64+UNUDf45cpujDNG7kjzgWack9S6oGYhiUXAoSlAcNcOdfrMAYR4zD/4vNSoFciQIAAgY4FBAA79nGUAAECBAgQIFBRAqk30pod++Ow3D25efhS496weG6YOX5ku+1M56fhjktW7wwPrN6RC/qtiL39epPSAhcL41xxKZ82Y2ycN25MOHnqmNjDb6QVTXsD69q6EEi9XJvn5+vdMOi0EEshaLjnYHPQMAUU9x6KAcUDMVAYhzin+lSXtqlub9qmcv5YCtpVSir0QNy8p/s9EFPvwUJAMM17WNgvbIt7Gxbvp+MWUKmUnwDtIECgrwUEAPta2P0JECBAgAABAr0USD1q7l65Lfx2xbZw14qtYcPug5k7fuOeVeHbN14aTo3BuJTWxgDh757aFq/ZHu55envY2sMhfeles2NQ76zZ48LpM8flgn2nzhgXTpg0SqAv4UgEBlAgrS6dFm5JuSepMLS5ORhYCBJmA4bNAcR8QLEQTCwKLKbhzpWQ0grT6b9zPflvXVqRuzhQ2BogPH4gMS04lK5JKzib97ASfgK0gQCBrggIAHZFyTkECBAgQIAAgX4USD177l+1Iwb7moN+S+MiHB2taLpz/5Hw5q/eG154+rRc0C/1EOxJmhXn4Tt7zvhw9uyY50wIZ80aFybHFXYlAgRqT6B4aPO0cT17v/Tfqlyvw3wPw9KAYUtwMR80TAuvFHLqkViOqQd61vLWq47GIdnb4z+ypNzd1BI8TCuP5+c3LAQHU4AwrUhevB0/Mq5Onq9LAVyJAAEC/SkgANif2p5FgAABAgQIEDiOwO4YxPvN8i3hvx/bHO5YvjX3S/VxTm23OvV8+Y/71rZ7rL3KNIx30dwJ4cITJ4YLYk77UwT72qNSR4DAcQRSECststLThVYOHW3MDU1uDQo2BwhTIDH9N7FQX7xNAcdUTvMjDnTqTfAw9R7M9TyMvQlbg4f5uQ1zAcXmXoaFgGEhkJgWVNLrcKC/vOcTqE4BAcDq/G5aTYAAAQIECNSAwDPbng2/fHxzLuiX5uhLq/F2NaV5q7oz/C4twpECfReeOClcMG9i7N033qqbXcV2HgECfSIwfMjgMHVsyt3vaXw09j5M8xkWBwdTj8NCuXi/pS4fWNwbg4cd9aruk5ctuWlu1eg4hLp0SoeS09oUBxcNWS4EBXPbfA/E8ccJKKZzkrdEgED9CggA1u+39+YECBAgQIDAAAis3v5s+PEjG8OP/rChZRGPrjQjrZB5yfxJ4YpTpobnnDIlLJg2Jk72fzS88cu/D0vX72lzi4mxB8llJ08Jly2YnNvOmzxKr5E2SioIEKhWgSG96H2Y/rElLaBSCAymba7XYVEAsXAsN4y5pL4b/1ZTdt7U9p6utpz+HintUdgcRGzubZgNKObr4t8lY2Ovw7R4jUSAQHULCABW9/fTegIECBAgQKAKBNbt3B9+EoN+KfD36PrdXW5xmovvOQunhCsWTA3nnzihTe+NtJLoN956cfjwj5aFh9fuCidNGR0uzwf9To+LdfiFrcvUTiRAoI4Ecr3o0lx8MXc3HUvBw8MxeJgforyrsD1wOKT9FDAsrtsd5zrcvT8ei/Xd6bXd3XZ15fzcSsu7G8PGkoWkOrs2xf7S4ihpqPKMOFdsWgiqdK7DCXF+w5YAYj5oaKhyZ7KOE+hfAWH8/vX2tFaBOXE3N1HR2rVrw5w5qSgRIECAAIHaEdgZJ5S/9eH14Yexp9+Da3Z16cXSvHxXLJgSXnD69JinhenjRnTpOicRIECAQOULpNWKCz0Ld+cDhqlcCBg27x/OnZMChrlzY1Ax7XdniohKkEhB1nEjhrQECgs9D3PzHeZXUU77qb5wLC2SkoKIw+LfhVJ5BdatWxfmzp1buGnaWVco2NaPgB6A9fOtvSkBAgQIECDQxwLpF7S7VmwN312yLjev3+E4R1VnKS288YLTpoUXnjE9F/wbGef2kwgQIECg9gRSYCvNd9jdOQ+b4oSFac7AXbEnYXFQMBs8zAcOU8Aw3ysxHR+oxVLS34dphfqUu5tKF0hpCRCmYGFRT8PmYGJhpeVhIV2n12F3tZ1fTwICgPX0tb0rAQIECBAg0CWB9MtWWpTjP+5bE9Zs3x9efs7M8JZL5x13SO2quJjH9x5Yl8ub9hzs9BlpxcyXnTUjXLNoVrho3qSQekpIBAgQIECgPYEU1Eqr/6Y8Z2J7Zxy/7khaLCX1MizqUdgcOGwellwIKObq0nm5IGOaH/FwONLY9YWpjt+C7h/p6QIpQ+LfpSlYWBiunIYppyBhrpwLHjYP+84FEQvleDz1OkxzSkoEal1AALDWv7D3I0CAAAECBLoskIZnpcU5/u3uZ8KyDa0La6RgYFpt8q9ecErLvQ4dbQw/fXRj+PZ9a8O9z+xoqT/eThoK9bKzZoaXL5oZLp0/2S8bx4NST4AAAQJlExgaA1uTY0/zlLuT0j+EpTkDC70J0zYtPrJyy77Yqy/fEzEfMMwFF/O9Do/Gnn8DldKzt+07nMvdbUNa6CTNCVmYx7AQJMyVY4CwuRdi6zyHhV6JaWEVvQ67q+38gRIQABwoec8lQIAAAQIEKkYgrf74H/euCV+7e1U4Xg++z/1ieVg0d0I4Ja6+++/3rs4F/rbHX4Y6Smm414vj0N5Xnz87t5CHeY060nKMAAECBCpFIAW1Rg0bksuzJozsUrOON1S5JYgYexWmxVNay3Gew3wwMfX6G8i099DRkPK6nQe61YxhMcD62/c9L0wba87ebsE5eUAEBAAHhN1DCRAgQIBA9QukoNmmuJJgCohV679+b9h1IAb9nolDfdd2Ok9S7AwR/sc3HghpXr/OJmNPq/e+5sI54RVxiG8agiQRIECAAIFaF+jNUOXWBVKycxm2Dl1uHbLcEkDMBw8HsNNh7v8TjBvR/dWka/1nwftVpoAAYGV+F60iQIAAAQIVK5CGAH329idz8+Ol/9P9wrha7ZffcmFVBQHX7tgfvvCrleH7D64LHQ1XSnP1pfctpDQc6nhpYhw69KrzZofXXDA3nDFr3PFOU0+AAAECBAiUCPR0gZRj8f+I7Dsc5yzM9CyMwcL8kOSW+Q4z5TQn4uFw8EjnC3WVNLNNccTQQWFEHAYsEagGAQHAavhK2kiAAAECBCpAIE0k/o17Vod/jENh03x4hfSLx7eE+1ftDItPmlSoqthtWtDji7/uPPCX3uXtV84Pzzt1aviLbz0Yfr5s83HfKc3n9+ZLToyr+E4Lw4f4JeC4UA4QIECAAIEyCwyKC3+kHngpz+3m/w05GP9Rr3iBlBQ0bFlpOTe/YQoUNtcVzkvnpBEQaVRASmmOQIlAtQgIAFbLl+rbdqYJCz4R84UxL4g5/adzV8xPxfyVmL8Zc/fXb48XSQQIECBQGwL3PLU9/K9bl+Ym/27vjZ7ctKeiA4BdCfylhXhfGlfmTYG/805oXWbx069ZFJZvvjs8E1f6LaRRwwbn5vVLKwMvnD62UG1LgAABAgQIVIlA6rmX8rRx3Zu/L/U63Bv/IbRcvQirhEsza0BAALAGPmIZXmFMvMefx3xfzD+JeWvM6Tefl8X8bzG/Pr/f+z7S8UYSAQIECFSPwLZ9h8InfvJ4+M+H1nfY6PW7DnZ4fKAObt17KPzTL1fkhisfb6jv8LhQx2svnBtuuPKkcOLk0W2amnoVfOcdl4T/fduTIQ0lumLB5PDqC+bkehu0OVkFAQIECBAgUNMCqddhbsXgOPWHRKCaBAQAq+lr9V1bd8Rbj4+5dZKj5meln4//jvnFMadgYAoOSgQIECBQBwJpJb9v3782fPKnj2eG+x7v1dNiGpWU9sf5gL5y1zPhS3c8FY63smCat+fNF58Ybrxqfqer96XV/T772kWV9IraQoAAAQIECBAgQKDLAgKAXaaq6RNTz77S4F964TTB03/F/NyYF8QsESBAgEAdCKze/mx43/cfDfc8vb3dtz1r9rhw+oxx4bsPrGs5vr5CAoBH4zyFqV2f++/lYUvs/dde6k7gr73r1REgQIAAAQIECBCoNgEBwN5/sWnxFovz+aK4TXlyzCl9Pebr0k4X04nxvL+K+eqY58acfnNJ8/DdEvMXY94fc3+mQfFhL80/cGl/PtizCBAgQKD/BRrjnDZfu/uZ8Jm4wm97K+ONHT4kvOelp4Y3xV5zv3x8cyYAWAk9AO9asTV89EePhRVb9rWLl1YY/LO4WMc7utDjr90bqCRAgAABAgQIECBQpQICgL3/cMdfFrB7974mnv7NmMcVXTYq7qeFOVK+IeYUGFwZc1+lYfHGfxNznAY9F8R8QdyeFvPXYv5lzBIBAgQI1KhAWiTjXbc8HJas3tnuG75i0azwwatPb5koe9aEkZnzNu85GNIqwUMHp3876t+Ueh9+7MePhduWbmr3wQ3xb7U/OX9OeNeLFobSdrd7gUoCBAgQIECAAAECNSYgAFjeD7om3u6JmNOced1J58WTvxNz+m0qdVv4ZMy/zpfTAhxvj3lhzGkOvhQM3BtzX6QUAPxQ0Y2b4v5nYn5/UZ1dAgQIEKghgTTX33eXrAsf+dGydufKmx0DfR//47PCc09NHd5bU6ovTrHzYNi0+2CYOyn921X/pENHG3Pz/H3hVyvDgSON7T70qoVTw/tedlo4fWbxv6+1e6pKAgQIECBAgAABAjUrIADY+0/70XiL+/N5c9zOi/mZmLuTPh9PTr9JpTn3UvDwnpgL6VdxZ0XMn4o5BQHfHfOHYy5Nn40Vw0srOyinZ6b7FqcUfEy9/1L3jVkxp16Jn4j50pj/KOY9MUsECBAgUCMC2+MKv+//z0fD7Y+lv77apjRc9r0xeDYmDv0tTRPiynejhg0O+w+3Bt7SMOD+CgDeuXxr+PAPl4Wntz1b2rRc+cxZ48L7X3Z6uOKUKe0eV0mAAAECBAgQIECgngTa/j/6enr78rxrcY+5ntwxzR94Zf7Cr8ZtcfCvcL8U3Ls+5tNjfmfMH4/5SMzF6R2xMLq4opP978XjpQHAwiVpUZB1Mf+fmLfFfEvMH4j5vTFLBAgQIFADAvc8tT2889sPtbtQxtxJI8On/3RRuGR+YUrbti/cEMfVpuG0K4vm29uwu+9XAt4Wg5Yp8PfjRza2bVSsmRgDk+996WnhtRfODYMGpX/TkggQIECAAAECBAgQEAAc+J+BVxU14WtF+8W7KSB3c8yfjHlCzM+L+faYi9OY4kIZ9wvPeW4Z7+lWBAgQIDBAAmmhj3/65Yrwz79aEdKw3dL0+ovmhg++/Ix2e/2VnlsaAFy/s+8CgGmo8q0Pb8gNVd65v/TfwGL39Rjre+PiE8J7XnJqmDAqzWghESBAgAABAgQIECBQEBAALEgM3PaK/KPTGKYHOmjGHUXHLo/7hcBcUXWf7KahwCm1/W2rud6fBAgQIFAlAmmhjtTr7/dP72jT4smjh4W//5NzwovOmN7m2PEqZk8YkTm0ftfBTLlchTS0+IM/WBp+9cSWdm+5aO6E8HevPDOcMyf9G5lEgAABAgQIECBAgECpgABgqUj/l9Ow3pTS6r5pDsDjpbS4SCEVrimUe7s9I95gVcz7S240Kpb/IV/305JjigQIECBQRQL3Pr09/MW3Hgzb9h1u0+rnxIUyPvuaRWHq2OFtjnVUUboQSArUlTMdi10Uv3XfmvD3tz0R9h1q+1fk+JFDcwt8vM5w33KyuxcBAgQIECBAgEANCggADuxHTV0nCrOTr+ukKTvj8dRLMM3zN7eTc7t7+LXxgnfF/NuYV8WcFvuYHfPLYk4TQN0V8+di7k6a08nJMzo57jABAgQIlEEgDZ296Xerwsd/8ng4WjLmd3CcIy8Nmb3xyvk9mi8vDQEuTuvLGABMwcR33fJwu70V0zP/6OwZ4SOvOKvbQcvi9tonQIAAAQIECBAgUC8CAoAD+6XHFj0+rcDbWSoEAMs939+P44PTUN/LYk4r/qb77475kZi/HfO/xdy260Ws7CCt7eCYQwQIECDQDwIH4gq97//PR8IP4tx5pSn13vunN5wbLjhxUumhLpdLA4ApaJcCjmmBkN6kny3dGN77/UfD7gNtZ59IvRTTcN+XnjWzN49wLQECBAgQIECAAIG6EhAAHNjPXTx5UtsxWW3bdihfle1y0fa87tYsiRekLBEgQIBAjQhsjCvy3vD1JWHZhtSpO5uee+rU8I+vO7fXi2WUDgHeHwOOKWjX00U4UsDy737yWPjWvWuyDc6XXnPBnPDBq88I4+NKvxIBAgQIECBAgAABAl0XEADsulVfnFk8W3pXliwsTM5U3kmW+uLNOh+mnIYA3983j3ZXAgQI1LfAI+t25YJ/W/YW/t2o1eOvnr8gvPOFC0Ma/tvbNGP8iNzqu7HTX0tKw4B7EgB8YtOe8P9866GwYkvbDvEp0PjJV58d0lyFEgECBAgQIECAAAEC3RcQAOy+WTmv2Ft0s64M603z/6XU9rej5vpK+nNdJTVGWwgQIFAvAj99dGNu7ryDR45lXnnM8CHhH167KLz4zPJNwTp08KAwfeyIsCmuLlxI63ceCGfOGl8odrpNQ4a/8fvV4WNxjsLDR7NtThe//JyZ4eN/fHZIC35IBAgQIECAAAECBAj0TEAAsGdu5boq/ca0Pea00EZni2ZMjOcUAoDm14sYEgECBAi0CqRA2r/85qnw6Z8/2VqZ35s/ZXT48rUXhpOnduXfmtpc3mHFrAnZAGB3VgJOK/v+z1v+EH62bFObZ4wcOjgu8nFmeM2Fc3o9p2Cbm6sgQIAAAQIECBAgUGcCAoAD/8Efi024MuYFMafvcbzFNk6Lxwrp8cKOLQECBAgQaIyr+/7trUvDv7czd95lJ08O/+dNF/TZvHmzJ44KD67Z1fIRNuxu7Q3YUtnOztNb94Ubv/FAWNnOkN8zZo4L//zG8/okYNlOU1QRIECAAAECBAgQqHkBAcCB/8S/jU1IAcDUu++CmO+Nub10VVHl3UX7dgkQIECgjgUOHmkM7/z2Q+Hnyza3UXjD4rnho688K6Shun2VUg/A4pSGAHeWfvXE5tjmh8Peg23/zeutl58U3vuyU8PwIYM7u43jBAgQIECAAAECBAh0UUAAsItQfXjaD+K935+///Vx214AMP3m9pb8Oambxa/z+zYECBAgUMcCacXdt8eVfu9btSOj0BDX9/jAH50e3nbFSX0+fLZ0JeC0CMjx0rHYU/GLv14Z/uEXy0PxwiHp/DTH3+detyg8/7Tpx7tcPQECBAgQIECAAAECPRQQAOwhXBkvuy/e666YUy/At8X89Zjvibk4vTsWTs9XfD5ujxQftE+AAAEC9SewOS68ce2/3Ree2FS8nlQIw2Jvv8+//tzwsrNn9gtKaQDweHMAdjTf32kzxoZ//bMLwwmTR/VLmz2EAAECBAgQIECAQL0JCAD2/otfEW+xoOg2U4r2U/11ReW0e1NJORXfGXMa1jsy5ttj/kTMqZdfKr8+5htjTml5zJ/N7fmDAAECBOpWYN3O/eFNX7k3rN6+P2MwNq70mxb7uGR+Wluqf9KsCemvqta0Ze+hcOhoY2YI7zPbng033rwkrGhnvr+0yu+n/vScMGqY/0vSqmiPAAECBAgQIECAQHkF/L/t3nveEG9x7XFuc3msT7k43VRcyO8/FLevi/mbMY+LOQUAS1MK/l0dc7arR+lZygQIECBQ0wKrtz8b3vjle0PpUNupY4eHr1+/OJwxK/010n+pNACYnrwpLgRy4uTmhevve2ZHXOxjSdi1P9t5fVAcpvzel54WbnzO/D4fptx/Gp5EgAABAgQIECBAoDIFBAAr57v8KDblnJhTb8AU6JsT8+GYV8b83Zi/EHO2q0eskAgQIECgfgTSirlv+srvw+Y9hzIvfdKU0eHmty4Ocyf1/xDacSOGhDGx52Ea4ltIKTiZAoC3Prw+vOe7j4TDjccKh3LbNN/fF+Iqv1eeMjVTr0CAAAECBAgQIECAQN8ICAD23vW6eIuUy5FWx5u8K5/LcT/3IECAAIEaEXhi057w5jjsd9u+9G9DrWnh9DHhmzdcHKaNza7G23pG3+41xBVH0jyAT25u7aCeVgL+wq9WhM/cnjqvZ5P5/rIeSgQIECBAgAABAgT6Q0AAsD+UPYMAAQIECPRCYHkMrqVhvzuezQb/zozDfb/xtovDpNHDenH33l86a8KITADwszHwtykuUlKaXnDatPBPbzgvjI49BiUCBAgQIECAAAECBPpPwP8D7z9rTyJAgAABAt0WeHrrvnaDf+fOnZCb82/8qKHdvme5LyidB7C94N+1l54Y/vaaM8PgNPmfRIAAAQIECBAgQIBAvwoIAPYrt4cRIECAAIGuC6yJq/ymnn/b9mXn/Fs8b1L4t+svys291/W79d2ZsydmVwIuflIcIRw+ePUZ4a2Xz7PYRzGMfQIECBAgQIAAAQL9KCAA2I/YHkWAAAECBLoqkBbSeMOXf99mKO2FJ04MX4vBv0oaRpvmAGwvDRs8KPzj688Nf3T2zPYOqyNAgAABAgQIECBAoJ8EBAD7CdpjCBAgQIBAVwW2xPnz3hSDfykIWJwWxWG/lRb8S+0rHQKc6tLKwP/6lgvCZSdPSUWJAAECBAgQIECAAIEBFBg0gM/2aAIECBAgQKBEYM/BI+Har90fVsXhv8XpjJnjws3XLw5jRwz8nH/F7Ur7Z88en1mIZMqY4eHbN14i+FcKpUyAAAECBAgQIEBggAT0ABwgeI8lQIAAAQKlAgePNIYbb14SHt+4J3No4fQxcbXfxaESFvzINCxfGDF0cK59//Kbp8LY2PPvL5+/IMyZOKq9U9URIECAAAECBAgQIDAAAgKAA4DukQQIECBAoFSg8VhTeNctD4ffP70jc2je5FHhmzdcHCbHXnWVnM6cNT588Y3nV3ITtY0AAQIECBAgQIBA3QoYAly3n96LEyBAgEClCDQ1NYWP/GhZ+OmjmzJNSkNpb37rxWHa2BGZegUCBAgQIECAAAECBAh0R0AAsDtaziVAgAABAn0g8MVfrww337M6c+e0iMZNcbXfE2IPQIkAAQIECBAgQIAAAQK9ERAA7I2eawkQIECAQC8Fbn14ffjM7cszdxk6uCH8659dEM6Ki2tIBAgQIECAAAECBAgQ6K2AAGBvBV1PgAABAgR6KPDgmp3hPd97JHN1Q0MI//Dac8NlC6Zk6hUIECBAgAABAgQIECDQUwEBwJ7KuY4AAQIECPRCYP2uA3HF3wfC4aPHMnf54NVnhGsWzcrUKRAgQIAAAQIECBAgQKA3AgKAvdFzLQECBAgQ6IHAvkNHw9tuuj9s23coc/WbLzkhvPXyeZk6BQIECBAgQIAAAQIECPRWQACwt4KuJ0CAAAEC3RA42ngs/PW3HwpPbNqbuerKU6aED11zZmhIY4AlAgQIECBAgAABAgQIlFFgSBnv5VYECBAgQIBAicCzsbffA6t35nKa8++hNbtC6gFYnE6eOjp84Y3nh6GD/btcsYt9AgQIECBAgAABAgTKIyAAWB5HdyFAgAABAm0E7lqxNfz5Nx9sE/ArPnHCqKHhq9deFMaPHFpcbZ8AAQIECBAgQIAAAQJlE9DVoGyUbkSAAAECBFoF1u7YH/7i3zsO/g0d3BD+/zdfEOZNGd16oT0CBAgQIECAAAECBAiUWUAAsMygbkeAAAECBNLKvn/5Hw+FPQezQ32LZRZOHxO+fv3icMn8ycXV9gkQIECAAAECBAgQIFB2AUOAy07qhgQIECBQ7wKfvO3x8Ie1uzIMKeB31cKp4YITJ4bzT5gYpo0bkTmuQIAAAQIECBAgQIAAgb4SEADsK1n3JUCAAIG6FPjZ0k3ha3evyrz7SXGI7/f//LIwdoR5/jIwCgQIECBAgAABAgQI9IuAIcD9wuwhBAgQIFAPAmu27w/v+d4fMq86bMig8MW4wq/gX4ZFgQABAgQIECBAgACBfhQQAOxHbI8iQIAAgdoVOHS0MfzFtx4Me0vm/fvwNWeGM2aNq90X92YECBAgQIAAAQIECFS8gABgxX8iDSRAgACBahD41M+eDI+u351p6ivPnRXesHhupk6BAAECBAgQIECAAAEC/S0gANjf4p5HgAABAjUncPfKbeGrv30m817zp44On/jjs0NDQ0OmXoEAAQIECBAgQIAAAQL9LSAA2N/inkeAAAECNSWwe/+R8O5bsvP+DY/z/v3Lm84Po4dba6umPraXIUCAAAECBAgQIFClAgKAVfrhNJsAAQIEKkPgf926NGzaczDTmPe97LRw2gzz/mVQFAgQIECAAAECBAgQGDABAcABo/dgAgQIEKh2gVsfXh9++IcNmde48pQp4dpL52XqFAgQIECAAAECBAgQIDCQAgKAA6nv2QQIECBQtQIbdh0I/+sHSzPtHz9yaPj0ny4KgwaZ9y8Do0CAAAECBAgQIECAwIAKCAAOKL+HEyBAgEA1Chw71hT+53f/EPYcPJpp/sf/+KwwY/yITJ0CAQIECBAgQIAAAQIEBlpAAHCgv4DnEyBAgEDVCdz0u1Xhd09tz7T7VefOCi8/Z1amToEAAQIECBAgQIAAAQKVICAAWAlfQRsIECBAoGoE1mzfHz718ycy7Z0Ve/195JVnZeoUCBAgQIAAAQIECBAgUCkCAoCV8iW0gwABAgQqXqCpqSn8zX89Gg4eOZZp62deuyik+f8kAgQIECBAgAABAgQIVKKAAGAlfhVtIkCAAIGKFPjuA+vCb1duy7Ttzy45MVx28pRMnQIBAgQIECBAgAABAgQqSUAAsJK+hrYQIECAQMUKbNl7MHzsx49l2jczDv39/156aqZOgQABAgQIECBAgAABApUmIABYaV9EewgQIECgIgU+/MNl7a76O3aEob8V+cE0igABAgQIECBAgACBFgEBwBYKOwQIECBAoH2Bny/bFH766KbMwVcsmhWef9r0TJ0CAQIECBAgQIAAAQIEKlFAALASv4o2ESBAgEDFCOw5eCT87a1LM+2ZOGpo+NA1Z2TqFAgQIECAAAECBAgQIFCpAgKAlfpltIsAAQIEKkLg7297ImzecyjTlr+Nwb/JY4Zn6hQIECBAgAABAgQIECBQqQICgJX6ZbSLAAECBAZc4OG1u8K37l2TacdzT50aXnXu7EydAgECBAgQIECAAAECBCpZQACwkr+OthEgQIDAgAk0HmtqM/R31LDB4WOvOis0NDQMWLs8mAABAgQIECBAgAABAt0VEADsrpjzCRAgQKAuBG5ZsjY8sm535l3/+oWnhDkTR2XqFAgQIECAAAECBAgQIFDpAgKAlf6FtI8AAQIE+l1g1/7D4VM/eyLz3JOnjg7XXXZSpk6BAAECBAgQIECAAAEC1SAgAFgNX0kbCRAgQKBfBT57+/Kwc/+RzDM//Iozw7Ah/trMoCgQIECAAAECBAgQIFAVAn6TqYrPpJEECBAg0F8CS9fvDv9+7+rM41521oxw5SlTM3UKBAgQIECAAAECBAgQqBYBAcBq+VLaSYAAAQJ9LtDU1BQ+9MNlIa7/0ZJGDB0UPnD16S1lOwQIECBAgAABAgQIEKg2AQHAavti2kuAAAECfSbwnw+uDw+s3pm5/188d4GFPzIiCgQIECBAgAABAgQIVJuAAGC1fTHtJUCAAIE+Edhz8Ej45G3ZhT9OmDQqvP058/vkeW5KgAABAgQIECBAgACB/hIQAOwvac8hQIAAgYoW+OKvVoZt+w5l2viha84II4YOztQpECBAgAABAgQIECBAoNoEBACr7YtpLwECBAiUXWDN9v3ha3evytz3+adNCy84fXqmToEAAQIECBAgQIAAAQLVKCAAWI1fTZsJECBAoKwCn7zt8XC48VjLPYcMaggftPBHi4cdAgQIECBAgAABAgSqW0AAsLq/n9YTIECAQC8F7n16e7ht6abMXa69bF6YP3VMpk6BAAECBAgQIECAAAEC1SogAFitX067CRAgQKDXAseONYW/+8ljmftMHDU0/NXzT8nUKRAgQIAAAQIECBAgQKCaBQQAq/nraTsBAgQI9Erg+w+uC0vX78nc469fuDCMj0FAiQABAgQIECBAgAABArUiIABYK1/SexAgQIBAtwSePXQ0fPrnT2auWTBtTHjjxSdk6hQIECBAgAABAgQIECBQ7QICgNX+BbWfAAECBHok8KU7ngpb9h7KXPuBuPDH0MH+asygKBAgQIAAAQIECBAgUPUCfsup+k/oBQgQIECguwLrdx0IX7rz6cxlz1k4NTzv1GmZOgUCBAgQIECAAAECBAjUgoAAYC18Re9AgAABAt0S+Ewc+nvo6LGWawYPaggfjL3/JAIECBAgQIAAAQIECNSigABgLX5V70SAAAECxxVYsXlv+MHD6zPH37j4hLBw+thMnQIBAgQIECBAgAABAgRqRUAAsFa+pPcgQIAAgS4JfP6XK0JTU+upY4YPCf/vixa2VtgjQIAAAQIECBAgQIBAjQkIANbYB/U6BAgQIHB8geWx999PHt2YOeGtl88Lk0YPy9QpECBAgAABAgQIECBAoJYEBABr6Wt6FwIECBDoUKC099/Y2PvvbVfM7/AaBwkQIECAAAECBAgQIFDtAgKA1f4FtZ8AAQIEuiTw5Ka94aclvf+uv+KkMH7U0C5d7yQCBAgQIECAAAECBAhUq4AAYLV+Oe0mQIAAgW4JfP6XyzNz/40dkXr/ndSteziZAAECBAgQIECAAAEC1SggAFiNX02bCRAgQKBbAk9s2hN7/23KXJOCf+NH6v2XQVEgQIAAAQIECBAgQKAmBQQAa/KzeikCBAgQKBb4/C9WFBdD6v13/eV6/2VQFAgQIECAAAECBAgQqFkBAcCa/bRejAABAgSSwGMb9oTblmZ7/90QF/7Q+8/PBwECBAgQIECAAAEC9SIgAFgvX9p7EiBAoE4F/umX2d5/41Lvvyvm1amG1yZAgAABAgQIECBAoB4FBADr8at7ZwIECNSJQOr997Nl2d5/b79yfhg3wtx/dfIj4DUJECBAgAABAgQIEIgCAoB+DAgQIECgZgX+5TcrM++Whv1ed/m8TJ0CAQIECBAgQIAAAQIEal1AALDWv3DX3y/9LPxlzA/GvD/mPTHfGfMrYpYIECBQdQJrtu+PK/9uzLT7hrjy71i9/zImCgQIECBAgAABAgQI1L6AAGDtf+OuvGFDPOmWmP855nExfzXmb8d8asy3xpwCgxIBAgSqSuArv306HGtqbfLoYYPDWy6b11phjwABAgQIECBAgAABAnUiMKRO3tNrdizwJ/FwynfH/KKYD8Sc0t/EvCTmz8T845hXxSwRIECg4gW27zsUblmyNtPON158gpV/MyIKBAgQIECAAAECBAjUi4AegPXypTt+z1fmD38ibgvBv1S1LebPxTw85utjlggQIFAVAjffszocPHKspa1DBjWEt8bhvxIBAgQIECBAgAABAgTqUUAAsPdffVq8xctj/mjMt8WcgmZp0FnKN8XcnXRiPPmzMT8R87Mx74j5/pjfE/OomPsqzcjf+Jl2HlCoe347x1QRIECg4gT2Hz4abr5nVaZdrzx3dpg5fmSmToEAAQIECBAgQIAAAQL1ImAIcO+/9Obe3yJ3h2vin9+MOc3BV0gp6HdhPt8Qt1fHnF3SMlaUIaWgZUqpe8zjub3WPwpdZha2VtkjQIBA5Qp8d8m6sHP/kUwDb3zO/ExZgQABAgQIECBAgAABAvUkoAdgeb/2mni723twy/PiNd+JOQX/9sX8gZgvi/kFMX855pRSAO4nMY9NhTKn1HMxpffFPCK31/zH5Lj563x5QlG9XQIECFSkwNHGY+HLdz2dadvzT5sWTp3RF//pzDxGgQABAgQIECBAgAABAhUroAdg7z/NR+Mt0jDdlFNvwHkxF4bNxt0upc/Hs9LYtKMxvzjme2IupF/FnRUxfyrmFAR8d8wfjrk0paHDaa6+rqb0zHTflL4V83UxPy/mR2P+WcxDY35VzIUejq2TacVKiQABApUo8NOlm8K6ncVTmYbwDr3/KvFTaRMBAgQIECBAgAABAv0oIADYe+wP9fIWi+P1V+bv8dW4LQ7+FW6dgnvXx3x6zO+M+eMxH4m5OL0jFkYXV3Sy/714vBAATIHHl8WcegC+MeYbY94d83/F/JmYl8e8JWaJAAECFSvQ1NQUvnTHU5n2nTt3Qlh80qRMnQIBAgQIECBAgAABAgTqTcAQ4IH/4qmXXSF9rbBTsk29727O16WhuKmnXmkaEysaupF/U3KDQ7H8kZhPjTn1JEyLm6Sg4uyYU1rSvPEnAQIEKlPg7pXbw7INezKN+x9XzQ8NDek/jRIBAgQIECBAgAABAgTqV0AAcOC//RX5JqRVfx/ooDl3FB27vGi/r3fflH/At/v6Qe5PgACB3gh86c5s77+TpowOLzqjsMh5b+7sWgIECBAgQIAAAQIECFS3gADgwH+/NKw3pbS6bxqKe7z0RNGBwjVFVb3eLV59uHCzP407b405zW/4n4VKWwIECFSawLINu8NdK7ZlmvX2K+eHwYP0/sugKBAgQIAAAQIECBAgUJcC5gAc2M+eVtydkm/Cuk6asjMeT70E0zx/czs5tyeH740XrY358ZgPxpzmJnxuzGk5zdfE3Bhzd9KcTk7WLacTIIcJEOi6wE13r8qcPGXMsPDq8wszGGQOKRAgQIAAAQIECBAgQKDuBAQAB5alc4gAAEAASURBVPaTjy16/L6i/ePtFgKAab6/cqfvxBu+OuZLYk4rAKeVjD8W86djzk6qFSu6kFIwUSJAgECfC+x89nC49Q8bMs95y6XzwoihgzN1CgQIECBAgAABAgQIEKhXAQHAgf3yqQdgIR0u7HSwTQt1pDSyeVPWPz8c75ayRIAAgaoS+M6SteHw0bRWUnMaOrghvGHxCYWiLQECBAgQIECAAAECBOpeQABwYH8E0lDbQhpW2Olgm1bnTelA86ai/+xsmHIaApzmFpQIECDQY4HGY03hG/eszlx/9dkzw9Sxhf9cZg4pECBAgAABAgQIECBAoC4FBAAH9rPvLXp8V4b1pvn/UurKcOHmMwfuz3UD92hPJkCgXgR+/cSWsH5X9t9E3nLZvHp5fe9JgAABAgQIECBAgACBLglYBbhLTH12UuoBuD1/984WzZgYzysEAM2v12efxI0JEKgmga/fsyrT3LNnjw/nzZ2QqVMgQIAAAQIECBAgQIBAvQsIAA78T8Bj+SYsiNuOemSeVtTUtFKvRIAAgboWeGrrvnDXim0Zg7dcemJoaGjI1CkQIECAAAECBAgQIECg3gUEAAf+J+C3+Sak3n0XdNCcq4qO3V20b5cAAQJ1KVA699/EUUPDNYtm1aWFlyZAgAABAgQIECBAgEBHAgKAHen0z7EfFD3m+qL94t30nd6Sr9gVt78uPmifAAEC9Saw79DR8P0HslONvu6iE8KIoYPrjcL7EiBAgAABAgQIECBAoFMBAcBOifr8hPviE+7KP+VtcXtpO098d6w7PV//+bg90s45qggQIFA3Av/10PqwNwYBC2lQHPX7potPKBRtCRAgQIAAAQIECBAgQKBIoKM554pOs9uBwBXx2IKi41OK9lP9dUXltHtTSTkV3xlzGtY7MubbY/5EzKmXXyq/PuYbY05pecyfze35gwABAnUq0NTUFG7+3arM27/g9Olh7qRRmToFAgQIECBAgAABAgQIEGgWEADs/U/CDfEW1x7nNpfH+pSL003Fhfz+Q3H7upi/GfO4mFMAsDSl4N/VMe8tPaBMgACBehK45+ntYcWWfZlXvvbSeZmyAgECBAgQIECAAAECBAi0ChgC3Gox0Hs/ig04J+bPxZyCfftjTvP9LYn5vTGfF/PKmCUCBAjUtcDNv1udef/5U0eHyxdMztQpECBAgAABAgQIECBAgECrgB6ArRY93bsuXphyOVL6rfZd+VyO+7kHAQIEakpgw64D4fbHNmXeKfX+a2iIkwBKBAgQIECAAAECBAgQINCugB6A7bKoJECAAIFKFPjO/WvDsabWlo0eNji8+vzZrRX2CBAgQIAAAQIECBAgQKCNgABgGxIVBAgQIFCJAo0x8ve9B9ZlmvbHMfg3dsTQTJ0CAQIECBAgQIAAAQIECGQFBACzHkoECBAgUKECv3tqW1gfhwAXp9dfdEJx0T4BAgQIECBAgAABAgQItCMgANgOiioCBAgQqDyBNPy3OJ0xc1w4a/b44ir7BAgQIECAAAECBAgQINCOgABgOyiqCBAgQKCyBHbtPxxuX7Y506jXXjgnU1YgQIAAAQIECBAgQIAAgfYFBADbd1FLgAABAhUk8IOH1ofDjcdaWjRs8KDwqvMs/tECYocAAQIECBAgQIAAAQIdCAgAdoDjEAECBAhUhsAtS7KLf7z4zOlhwqhhldE4rSBAgAABAgQIECBAgECFCwgAVvgH0jwCBAjUu8DS9bvDYxv3ZBhed9HcTFmBAAECBAgQIECAAAECBI4vIAB4fBtHCBAgQKACBG5Zkl38Y/aEkeHyk6dUQMs0gQABAgQIECBAgAABAtUhIABYHd9JKwkQIFCXAgePNIY0/19x+tML5oRBgxqKq+wTIECAAAECBAgQIECAQAcCAoAd4DhEgAABAgMr8PNlm8Keg0dbGtEQ434pACgRIECAAAECBAgQIECAQNcFBAC7buVMAgQIEOhngdLhv2no79xJo/q5FR5HgAABAgQIECBAgACB6hYQAKzu76f1BAgQqFmBtTv2h7tXbs+832su1PsvA6JAgAABAgQIECBAgACBLggIAHYBySkECBAg0P8C331gXeah40YMCS85c0amToEAAQIECBAgQIAAAQIEOhcQAOzcyBkECBAg0M8CjceawvdKVv991Xmzw4ihg/u5JR5HgAABAgQIECBAgACB6hcQAKz+b+gNCBAgUHMCd6/cFjbsPph5r9deODdTViBAgAABAgQIECBAgACBrgkIAHbNyVkECBAg0I8C338wO/z3jJnjwlmzx/djCzyKAAECBAgQIECAAAECtSMgAFg739KbECBAoCYE9h8+Gm5ftjnzLhb/yHAoECBAgAABAgQIECBAoFsCAoDd4nIyAQIECPS1wC8e3xIOHGlseczgQQ3h5efMainbIUCAAAECBAgQIECAAIHuCQgAds/L2QQIECDQxwI/fHh95gmXL5gSpo4dnqlTIECAAAECBAgQIECAAIGuCwgAdt3KmQQIECDQxwK79h8OdyzfmnnKKxbp/ZcBUSBAgAABAgQIECBAgEA3BQQAuwnmdAIECBDoO4Hblm4KRxqbWh4wbMig8JIzp7eU7RAgQIAAAQIECBAgQIBA9wUEALtv5goCBAgQ6COBW0uG/77w9Glh7IihffQ0tyVAgAABAgQIECBAgEB9CAgA1sd39pYECBCoeIFNuw+Ge5/ZkWmn4b8ZDgUCBAgQIECAAAECBAj0SEAAsEdsLiJAgACBcgv8+JENoal19G8YO3xIeO6p08r9GPcjQIAAAQIECBAgQIBA3QkIANbdJ/fCBAgQqEyBWx/ekGnYS8+aEUYMHZypUyBAgAABAgQIECBAgACB7gsIAHbfzBUECBAgUGaBp7buC4+u35256yvOtfpvBkSBAAECBAgQIECAAAECPRQQAOwhnMsIECBAoHwCPyzp/TdlzPBw6fzJ5XuAOxEgQIAAAQIECBAgQKCOBQQA6/jje3UCBAhUgkBTnPjvR3/IDv99+Tkzw5DB/oqqhO+jDQQIECBAgAABAgQIVL+A366q/xt6AwIECFS1wNL1e8LT257NvIPhvxkOBQIECBAgQIAAAQIECPRKQACwV3wuJkCAAIHeCtz68PrMLeZOGhnOmzshU6dAgAABAgQIECBAgAABAj0XEADsuZ0rCRAgQKCXAo3H4vDfR7LDf1+5aHZoaGjo5Z1dToAAAQIECBAgQIAAAQIFAQHAgoQtAQIECPS7wH3P7Aib9xzKPNfw3wyHAgECBAgQIECAAAECBHotIADYa0I3IECAAIGeCvzwD9nhv6fNGBsWTh/b09u5jgABAgQIECBAgAABAgTaERAAbAdFFQECBAj0vcCRxmPhtqWbMg965bmzM2UFAgQIECBAgAABAgQIEOi9gABg7w3dgQABAgR6IHDv0zvCrv1HMldes2hmpqxAgAABAgQIECBAgAABAr0XEADsvaE7ECBAgEAPBH62bGPmqkVzxoc5E0dl6hQIECBAgAABAgQIECBAoPcCAoC9N3QHAgQIEOimwLG4+u/tyzZnrnrxmTMyZQUCBAgQIECAAAECBAgQKI+AAGB5HN2FAAECBLoh8NDaXWHL3uzqvy8RAOyGoFMJECBAgAABAgQIECDQdQEBwK5bOZMAAQIEyiRw+7Ls4h8Lpo0JKUsECBAgQIAAAQIECBAgUH4BAcDym7ojAQIECHQg0NTUFH5WEgB8yZnTO7jCIQIECBAgQIAAAQIECBDojYAAYG/0XEuAAAEC3RZ4cvPesHr7/sx1Lz3T6r8ZEAUCBAgQIECAAAECBAiUUUAAsIyYbkWAAAECnQv8fGl28Y/ZE0aGs2aP6/xCZxAgQIAAAQIECBAgQIBAjwQEAHvE5iICBAgQ6KlA6fDfF8fhvw0NDT29nesIECBAgAABAgQIECBAoBMBAcBOgBwmQIAAgfIJrIlDfx/fuCdzQ6v/ZjgUCBAgQIAAAQIECBAgUHYBAcCyk7ohAQIECBxP4Ocli39MGj0sXDRv0vFOV0+AAAECBAgQIECAAAECZRAQACwDolsQIECAQNcESgOALzp9ehg8yPDfruk5iwABAgQIECBAgAABAj0TEADsmZurCBAgQKCbAlv2HgwPrNmZueolZ03PlBUIECBAgAABAgQIECBAoPwCAoDlN3VHAgQIEGhH4L8f2xyamloPjBk+JFx28pTWCnsECBAgQIAAAQIECBAg0CcCAoB9wuqmBAgQIFAq8LOlmzJVzz11ahgxdHCmToEAAQIECBAgQIAAAQIEyi8gAFh+U3ckQIAAgRKB3QeOhHue2p6pfelZMzJlBQIECBAgQIAAAQIECBDoGwEBwL5xdVcCBAgQKBL49RNbwtFjreN/hw0ZFJ576rSiM+wSIECAAAECBAgQIECAQF8JDOmrG7svAQIECNSPwOGjx8Jff+ehkIb5zhw/Mlx5ypSYp4bLF0wOE0YNy9UXa1yxYEpIcwBKBAgQIECAAAECBAgQIND3An776ntjTyBAgEDNC3z3gbXhp482z/G3fteB8O371+byoIYQzpkzITyxaU/G4KVnGv6bAVEgQIAAAQIECBAgQIBAHwoIAPYhrlsTIECgXgR+vmxzu6+aRv0+vHZX5lgKCr7gdMN/MygKBAgQIECAAAECBAgQ6EMBcwD2Ia5bEyBAoB4EDh5pDPc+nV3go6P3XnzSpDB5zPCOTnGMAAECBAgQIECAAAECBMooIABYRky3IkCAQD0K3PvMjnAozgFYSKmH33MWTg0jhrb/V8wrFs0unGpLgAABAgQIECBAgAABAv0gYAhwPyB7BAECBGpZ4I4nt2Ze77wTJoab37o4pJ6BD6zeGe5cvjXcuWJb2Lj7QHjxGdPDay+ckzlfgQABAgQIECBAgAABAgT6VkAAsG993Z0AAQI1L3DnimwA8KrY+y+lEUMHx1WAp+Ty+2tewQsSIECAAAECBAgQIECgcgXaH59Vue3VMgIECBCoIIF1O/eHlVv2ZVqUhv9KBAgQIECAAAECBAgQIFA5AgKAlfMttIQAAQJVJ3Dn8m2ZNk8cNTScPXt8pk6BAAECBAgQIECAAAECBAZWQABwYP09nQABAlUtkOb3K05XnjI1DE6rgEgECBAgQIAAAQIECBAgUDECAoAV8yk0hAABAtUlcKTxWLh7ZbYHoOG/1fUNtZYAAQIECBAgQIAAgfoQEACsj+/sLQkQIFB2gYfW7Ap7Dx3N3Pc5p0zJlBUIECBAgAABAgQIECBAYOAFBAAH/htoAQECBKpS4I7lWzLtPmPmuDBt3IhMnQIBAgQIECBAgAABAgQIDLyAAODAfwMtIECAQFUKlC4AYvhvVX5GjSZAgAABAgQIECBAoA4EBADr4CN7RQIECJRbYNu+Q+HR9bszt71q4dRMWYEAAQIECBAgQIAAAQIEKkNAALAyvoNWECBAoKoE7lqRXf139LDB4YITJ1bVO2gsAQIECBAgQIAAAQIE6kVAALBevrT3JECAQBkFSof/XrZgShg2xF8pZSR2KwIECBAgQIAAAQIECJRNwG9rZaN0IwIECNSHwLFjTeHO5dkegOb/q49v7y0JECBAgAABAgQIEKhOAQHA6vxuXWn1m+NJX4p5ScyHYm6K+bqYO0qr4sF0Xnv5N7FeIkCAQFi2YU/Y/uzhjMRVp5j/LwOiQIAAAQIECBAgQIAAgQoSGFJBbdGU8gp8LN7uxJi3xbwxvx83naY0q/8/tnPWqnbqVBEgUIcCd5bM/zd/yuhwwuRRdSjhlQkQIECAAAECBAgQIFAdAgKA1fGdetLKG+JFK2JeHfP7Yv5kzF1Ju+JJH+7Kic4hQKA+Be540vDf+vzy3poAAQIECBAgQIAAgWoVEACs1i/Xebt/0fkpziBAgED3BPYcPBIeWLMzc9FVCw3/zYAoECBAgAABAgQIECBAoMIEBADbfpBpsWpxPl8UtylPjjmlr8d8XdrpYkpDcP8q5qtjnhtzmovvqZhvifmLMe+PudLS8Nig62KeFfOemO+P+d6YJQIECITfrdwWGuMiIIWUVv69eP6kQtGWAAECBAgQIECAAAECBCpQQACw7UfZ3LaqRzXXxKu+GfO4oqvTJFkX5vMNcZsCgytjrqQ0IzbmayUNSkHAN8ScgpcSAQJ1LHDH8jStaGtaPG9SGDXMXyWtIvYIECBAgAABAgQIECBQeQJWAe74m6yJh2/v+JR2j54Xa78Tcwr+7Yv5AzFfFvMLYv5yzCktjPknMY9NhQpJKfCX2jg95tExp/f4RsypF+QvY66ktsbmSAQI9KdAU1NTuHP51swjDf/NcCgQIECAAAECBAgQIECgIgV022j7WT4aq1KPt5RTb8B5MT8Tc3fS5+PJI2M+GvOLY74n5kL6VdxJi3N8KuYUBHx3zB+OuTR9Nlak4bhdTemZ6b69SR8pufjhWH5Lvu7P4vbtMf9DvmxDgECdCTy1dV9Yv+tA5q2vOnVqpqxAgAABAgQIECBAgAABApUnIADY9pt8qG1Vt2rS/IFX5q/4atwWB/8KN0rBvetjPj3md8b88ZiPxFyc3hELqRdeV9P34om9DQAe71lfigdSAPDymAUAj6eknkCNC5QO/505fkQ4ZdqYGn9rr0eAAAECBAgQIECAAIHqFzAEuPzf8FVFtyydS69w6FjcuTlfmBC3zyscKNqm36obupF/U3RtuXcLk351JyBZ7ja4HwECAyyQFgApTs85ZWpoaEj/mZIIECBAgAABAgQIECBAoJIFBADL/3WuyN/y2bh9oIPb31F0LPWsq+R0cb5xqyq5kdpGgEDfCaSVf+97ZkfmAZefMiVTViBAgAABAgQIECBAgACByhQQACz/d0nDelNKq/umOQCPl54oOlC4pqiq33dPi09MqxSXplT/v/OV3yo9qEyAQH0ILNuwO+w9lP1P2iUnTaqPl/eWBAgQIECAAAECBAgQqHIBcwCW9wOOiLcrdIlZ18mtd8bjqZdgGlY7t5Nze3L4hnhRoTfi2fkbpLrn5vd/G7dfye+nzetjflfMd8a8OubUtrRIyR/FPDTmT8acjnU1zenkxBmdHHeYAIEKEvj909szrZk/dXSYNi79J08iQIAAAQIECBAgQIAAgUoXGFLpDayy9o0tau++ov3j7RYCgH0xi34K/l1b8uA01Lh4uHFxAPDX8VjqiXhezGkRk9QbME349dOY/yXm22PuTlrbnZOdS4BAZQv8/uns8N9L5k+u7AZrHQECBAgQIECAAAECBAi0CAgAtlCUZae4O8zhLtzxUP6ckV04t7unXBcvSLmrKc1JWDwvYVevcx4BAjUucLTxWLi/ZP4/AcAa/+hejwABAgQIECBAgACBmhIQACzv5zxYdLthRfvH2x2eP3DgeCdUcX1nw5rTEOD7q/j9NJ1A3Qgs27Cn7fx/883/Vzc/AF6UAAECBAgQIECAAIGqFxAALO8n3Ft0u64M603z/6XUleHCzWdWz5/rqqepWkqAQEcCpfP/nZzm/xtb3OG5o6sdI0CAAAECBAgQIECAAIGBFrAKcHm/QOoBWJgpv7NFMCbGcwsBQPPllfc7uBsBAmUUKA0AGv5bRly3IkCAAAECBAgQIECAQD8ICACWH/mx/C0XxG1HPSxPK3r040X7dgkQIFAxArn5/1alRctbkwBgq4U9AgQIECBAgAABAgQIVIOAAGD5v9Jv87dMvfsu6OD2VxUdu7to3y4BAgQqRmBpnP9v36GjmfYIAGY4FAgQIECAAAECBAgQIFDxAgKA5f9EPyi65fVF+8W7yf0t+Ypdcfvr4oP2CRAgUCkCpcN/F0wbE6aOLaxfVCmt1A4CBAgQIECAAAECBAgQ6EhAALAjnZ4duy9edlf+0rfF7aXt3Obdse70fP3n4/ZIO+eoIkCAwIALlAYAL7H674B/Ew0gQIAAAQIECBAgQIBAdwU6mqOuu/eqlfOviC+S5u8rpCmFnbhN9dcVldPuTSXlVHxnzGlY78iYb4/5EzGnXn6p/PqYb4w5peUxfza35w8CBAhUmEBu/r9ndmRaZfhvhkOBAAECBAgQIECAAAECVSEgANj2M90Qq65tW52ruTz+mXJxuqm4kN9/KG5fF/M3Yx4XcwoAlqYU/Ls65r2lB5QJECBQCQKPrt8dnj3cmGmKAGCGQ4EAAQIECBAgQIAAAQJVIWAIcN99ph/FW58T8+diTsG+/TGn+f6WxPzemM+LeWXMEgECBCpS4PdPZ3v/nRLn/5syxvx/FfmxNIoAAQIECBAgQIAAAQIdCOgB2BbnuliVcjnS6niTd+VzOe7nHgQIEOg3gbbz/03ut2d7EAECBAgQIECAAAECBAiUT0APwPJZuhMBAgRqRuBI47GwZFW2B6DhvzXzeb0IAQIECBAgQIAAAQJ1JiAAWGcf3OsSIECgKwLtzf93sRWAu0LnHAIECBAgQIAAAQIECFScgABgxX0SDSJAgMDAC5QO/1043fx/A/9VtIAAAQIECBAgQIAAAQI9ExAA7JmbqwgQIFDTAqULgBj+W9Of28sRIECAAAECBAgQIFDjAgKANf6BvR4BAgS6K2D+v+6KOZ8AAQIECBAgQIAAAQKVLSAAWNnfR+sIECDQ7wKPrNsd9h9uzDz34pMmZcoKBAgQIECAAAECBAgQIFA9AgKA1fOttJQAAQL9IlA6/9+p08eGyWOG98uzPYQAAQIECBAgQIAAAQIEyi8gAFh+U3ckQIBAVQuUBgAvsfpvVX9PjSdAgAABAgQIECBAgIAAoJ8BAgQIEGgRaJ7/b2dLOe1YACTDoUCAAAECBAgQIECAAIGqExAArLpPpsEECBDoO4FH1u0KB46UzP83f3LfPdCdCRAgQIAAAQIECBAgQKDPBQQA+5zYAwgQIFA9Ar9/ekemsafNGBsmjR6WqVMgQIAAAQIECBAgQIAAgeoSEACsru+ltQQIEOhTgbbz/+n916fgbk6AAAECBAgQIECAAIF+EBAA7AdkjyBAgEA1CBxtPBaWrCqd/29SNTRdGwkQIECAAAECBAgQIECgAwEBwA5wHCJAgEA9CTyxaW+b+f8umicAWE8/A96VAAECBAgQIECAAIHaFBAArM3v6q0IECDQbYEHVmd7/500ZXSYPGZ4t+/jAgIECBAgQIAAAQIECBCoLAEBwMr6HlpDgACBARN4cE02AHjeCRMGrC0eTIAAAQIECBAgQIAAAQLlExAALJ+lOxEgQKCqBUoDgOefMLGq30fjCRAgQIAAAQIECBAgQKBZQADQTwIBAgQIhC17D4a1Ow5kJC44UQAwA6JAgAABAgQIECBAgACBKhUQAKzSD6fZBAgQKKfAg6t3ZW43ZviQsHD62EydAgECBAgQIECAAAECBAhUp4AAYHV+N60mQIBAWQUeKpn/b9Hc8WHwoIayPsPNCBAgQIAAAQIECBAgQGBgBAQAB8bdUwkQIFBRAqUrAF9g/r+K+j4aQ4AAAQIECBAgQIAAgd4ICAD2Rs+1BAgQqAGBw0ePhUfW7868yXnm/8t4KBAgQIAAAQIECBAgQKCaBQQAq/nraTsBAgTKIPDYxj0hBQGL0/lzLQBS7GGfAAECBAgQIECAAAEC1SwgAFjNX0/bCRAgUAaB0uG/C6aNCeNHDS3Dnd2CAAECBAgQIECAAAECBCpBQACwEr6CNhAgQGAABR4sWQDk/BMmDGBrPJoAAQIECBAgQIAAAQIEyi0gAFhuUfcjQIBAlQk8uHpnpsUXmP8v46FAgAABAgQIECBAgACBahcQAKz2L6j9BAgQ6IXAxt0HwsbdBzN3ON8KwBkPBQIECBAgQIAAAQIECFS7gABgtX9B7SdAgEAvBB5cvStz9bgRQ8LJU8dk6hQIECBAgAABAgQIECBAoLoFBACr+/tpPQECBHolULoAyHmx99+gQQ29uqeLCRAgQIAAAQIECBAgQKCyBAQAK+t7aA0BAgT6VaDtAiAT+/X5HkaAAAECBAgQIECAAAECfS8gANj3xp5AgACBihQ4eKQxLNuwO9M2C4BkOBQIECBAgAABAgQIECBQEwICgDXxGb0EAQIEui+wdP3ucKSxqeXChjjyd9Hc8S1lOwT+L3v3AmzZWReIfp3u08/0M93pbtJJOpBIgpoogszwEhAvIyIj5XiViw50RsRxHoUFl4sjU4qUMnUZkWLmjjqXV+twVdRR5qIyMgLyKubyUItYkIQQkqY76U46/X6/zv1/yd7ptdZ57bPP2mevx++r+vdez2993+87ObvOP2utjwABAgQIECBAgAABAu0QkABsxzjqBQECBBYsUH7895bt67P1q1csuB4nECBAgAABAgQIECBAgEC9BSQA6z0+WkeAAIGRCZRnAP6eXd7/NzJsFRMgQIAAAQIECBAgQGCMAhKAY8R3aQIECIxLYGpqKvvy3iOFy39PzACsECBAgAABAgQIECBAgED7BCQA2zemekSAAIF5BfYdOZM9cuJc4bjvuWFTYd0KAQIECBAgQIAAAQIECLRDQAKwHeOoFwQIEFiQQPn9f5vXrsievPWqBdXhYAIECBAgQIAAAQIECBBohoAEYDPGSSsJECBQqcDfPDD98d+JNA2wQoAAAQIECBAgQIAAAQKtE5AAbN2Q6hABAgTmF/ibvUcLB5kApMBhhQABAgQIECBAgAABAq0SkABs1XDqDAECBOYXOH3+YvbVh44XDjQBSIHDCgECBAgQIECAAAECBFolIAHYquHUGQIECMwv8JV9x7JLl6eeOHD5sonsu67f+MS6BQIECBAgQIAAAQIECBBol4AEYLvGU28IECAwr0B5ApCnPWl9tnbl5LznOYAAAQIECBAgQIAAAQIEmikgAdjMcdNqAgQIDC0w0wQgQ1fmRAIECBAgQIAAAQIECBCovYAEYO2HSAMJECBQncDU1FQ2bQKQGzZXdwE1ESBAgAABAgQIECBAgEDtBCQAazckGkSAAIHRCTzw6Ons8KnzhQs8Y5cEYAHECgECBAgQIECAAAECBFomIAHYsgHVHQIECMwl8OUHjhR2b123Krtu85rCNisECBAgQIAAAQIECBAg0C4BCcB2jafeECBAYE6Bv/1WMQH4jF2bsomJiTnPsZMAAQIECBAgQIAAAQIEmi0gAdjs8dN6AgQILEjgK/uOFY7/7us9/lsAsUKAAAECBAgQIECAAIEWCkgAtnBQdYkAAQIzCZy7eCn72kPHC7u+67qNhXUrBAgQIECAAAECBAgQINA+AQnA9o2pHhEgQGBGgXsOnMwuXJoq7PuOnRKABRArBAgQIECAAAECBAgQaKGABGALB1WXCBAgMJPAV/YfLWx+8tarso1rVhS2WSFAgAABAgQIECBAgACB9glIALZvTPWIAAECMwrcWXr/3+0e/53RyUYCBAgQIECAAAECBAi0TUACsG0jqj8ECBCYRaA8AchtHv+dRcpmAgQIECBAgAABAgQItEtAArBd46k3BAgQmFHg7IVL2T0HTxT23X7dpsK6FQIECBAgQIAAAQIECBBop4AEYDvHVa8IECBQEEiz/168fGUCkImJLPuOazcUjrFCgAABAgQIECBAgAABAu0UkABs57jqFQECBAoCd+4/Vli/+Zp12VWrJgvbrBAgQIAAAQIECBAgQIBAOwUkANs5rnpFgACBgsC09/+ZAKTgY4UAAQIECBAgQIAAAQJtFpAAbPPo6hsBAgR6Al/Zd7RgcbsJQAoeVggQIECAAAECBAgQINBmAQnANo+uvhEgQCAETp27mN378MmCxW0mACl4WCFAgAABAgQIECBAgECbBSQA2zy6+kaAAIEQ+GpMAJKb/yNbvmwi+/YnmQDEDwcBAgQIECBAgAABAgS6IiAB2JWR1k8CBDorUH7/37dtW5etWbm8sx46ToAAAQIECBAgQIAAga4JSAB2bcT1lwCBzgncWX7/nwlAOvczoMMECBAgQIAAAQIECHRbQAKw2+Ov9wQIdEDgK/uPFXp5u/f/FTysECBAgAABAgQIECBAoO0CEoBtH2H9I0Cg0wInzl7I7nvkVMHgdncAFjysECBAgAABAgQIECBAoO0CEoBtH2H9I0Cg0wJ/v/94of8rlk9kt+xYX9hmhQABAgQIECBAgAABAgTaLSAB2O7x1TsCBDoucOf+owWBW3dsyFZNmgCkgGKFAAECBAgQIECAAAECLReQAGz5AOseAQLdFijPAHybx3+7/QOh9wQIECBAgAABAgQIdFJAArCTw67TBAh0ReDO8gQgOzd2pev6SYAAAQIECBAgQIAAAQI9AQlAPwoECBBoqcDR0+ezBx49XeidOwALHFYIECBAgAABAgQIECDQCQEJwE4Ms04SINBFgfLdfysnl2VP3W4CkC7+LOgzAQIECBAgQIAAAQLdFpAAbOf474xu/XzExyL2RpyPOBDxXyP+QcRs5adix3+O+FLEuYipiN0RCgECDRQov//v25+0IVux3K/9Bg6lJhMgQIAAAQIECBAgQGBRApOLOtvJdRX419GwN0d8IyIlAR+J+LaIV/TiVfH5oYhy+dXYsCviUMRDveX4UAgQaKLAnfuOFZp9uwlACh5WCBAgQIAAAQIECBAg0BUBt4K0c6S/EN16YcTNEa+N+DcRPxbxoohLEb8VsSqiXNKxN0ZcE/HbEQoBAg0WKD8CfJsJQBo8mppOgAABAgQIECBAgACB4QXGmQBMSaZhyw8Oe2JHzvuT6OenZujrZ2LbJyM2R9w2w/6/im0PzLDdJgIEGiZw6OS5bP/RM4VW337dpsK6FQIECBAgQIAAAQIECBDohsA4E4B3BvEPLZA53bX2nyL+fIHnLeTwbXHwD0e8LeKjEelx2PQuvBR7IhZS0uO074y4K+JUxOGIL0a8KWJtxDjKhd5FL47j4q5JgMDSCJTv/luzYnl287Z1S3NxVyFAgAABAgQIECBAgACBWglMjrE1KdH2kYj0qOkbI85GzFWeHjv/n4hb5jqogn0HK6gjVfHyiA9GbEgrvZKSfs/sRXrc9mUR90YsVbkhLvQDEen9fikBqxAg0FKB8vv/vnPnhmz5somW9la3CBAgQIAAAQIECBAgQGAugXHeAZgeNU1/jf7ziL+JSAm+mUo6Jr3D7vMRKfmX1n8/YinK3rhImkRjoSX1JU2ykZJ/JyPeEvGciBdHvCciladGpDsZ16eVJSgr4hr/JWJVRJogJL0LUCFAoKUC5RmAb9vp8d+WDrVuESBAgAABAgQIECBAYF6Bcd4BeHu07jcjfjIiJfY+H/HWiP8zIj1um0p6hDYlrZ4bkRJ/xyL+ZcTvRYyqpEd/02O6KdLdgDdGfDNiIeXdcfCaiPSY7UsiUt/65ROx8PWId0SkJGC6+/GtEeWSHh1OybpBS7pmqnemkhK9eyK+L+I9EclUIUCgxQJ37j9a6J0ZgAscVggQIECAAAECBAgQINApgXEmAE+E9D+NSHfBpURguj3l1yLSBB+vjnhhxH+ISHfIpeTfpyLS9m9FjLL88iIrf1ac//xeHe+Lz3zyr191Su7dEfG0iNdHpH5fiMiXn42Vq/Ib5ln+49g/UwIwJf/eH/GqiA9GpDsuFQIEWixw8PjZ7ODxc4Ue3nbdxsK6FQIECBAgQIAAAQIECBDojkBKDo27/EE04LsiPh2REn0peXZPxAci0iO06S66X4j4/ohRJ//iEosur8jVkPowU7kcG3+3tyMlPl80w0Hpbf3JY9D46xnqSOOb2vCaiN+P2B2Rrq0QINBigfLjv+tWTWZP3rKQ/5/QYhxdI0CAAAECBAgQIECAQAcFUoKoDiUl9lIS7EMRKeGV3leXPtMjv/8wIj0uOxXRhPK8XiPTrL9fnqPBn8rte25uuarFfvIv3TWZXNPdlpeqqlw9BAjUV+DOfcXHf9MEIMtMAFLfAdMyAgQIECBAgAABAgQIjFigLgnA9Cjy2yN+LCIl+lLyL5V0B2B6PHZ7WmlISY/1ppJm9013L85W7srt6J+T27SoxTSu6bHflPz7o4ifipD8CwSFQBcEvrI//b+TK+X269KNxgoBAgQIECBAgAABAgQIdFVgnO8A7JuniTB+L+LpESnxdyAizfr7LyK+N+IfRdwZ8TMR/y2izmV1NG5rr4H75mnokdif7hJMz+VdP8+xC939S3FCeuw3zUCcHqf+txHl8uHY8Helja+N9ef1tt3W+0zbXthb/mx8vre3PN/HdfMcsGOe/XYTIDCEwNTUVHbnvmIC8Lad3v83BKVTCBAgQIAAAQIECBAg0BqByTH3JE1I8esRayJS8u//jUgJp0MRH4z45YiUDExJtT+JSHe1vT7idEQdS5qwpF9S8m2+0k8Apvf9VVlu7FWW6n3LLBXfH9vLCcCU/EuJw3xJjyfnH1EeNAHYhPc15vtpmUArBB48djZ79NT5Ql/MAFzgsEKAAAECBAgQIECAAIHOCYzzEeCPhPZ/ilgbcSYizXr7ioiU/EslPbKa7mT7vohvRqQE4T+LSEmrNNNuHUu6A7Bfin+B97cWP/vTdKYEaJVld1SWvOaKPbG/XHbHhrnOSfsVAgRqLFB+/9+G1ZPZDVenX7MKAQIECBAgQIAAAQIECHRVYJx3AL6sh/6l+PzJiK/PMgifj+23R/zHiN0RN0d8JmJVRN3K2VyDVuaWZ1vs9yElQNtW5nusOT0C/MW2dVp/CIxboDwDcHr/38REyusrBAgQIECAAAECBAgQINBVgXEmAC8H+r+L+JWIi/MMQHpUNt39l+4afE/E5og6lhO5Rg3yWG96/18qgzwu/PiRzfl3X3OaqqUE2iNw14H8r6Es+07v/2vP4OoJAQIECBAgQIAAAQIEhhQYZwLwBdHmzy2w3X8ax//PiPQuwDqWdAfgoxFbIuabBCMlMfsJQO/LCwyFAIHFC9z10PFCJU97Uv7VpIVdVggQIECAAAECBAgQIECgIwLjfAfgQpN//SF5KBZe2l+p4edXe21KjyrPlWC9Ndf2r+WWLRIgQGAogWNnLmRpEpB8uWWHBGDewzIBAgQIECBAgAABAgS6KDDOBGBbvT/b61i6u+8Zc3Qy3QHZL8MmQ/vn+yRAgEB2z8Hi47+Tyyayp2wd5G0E8AgQIECAAAECBAgQIECgzQISgNWP7odzVd6RW84vJvdX9zYcjc9P5ndaJkCAwDAC5ff/3XTNumzlpF/zw1g6hwABAgQIECBAgAABAm0SmOsR1VH3czHv8ZuKxv30qBs4ZP1fiPPSLMXPj0ht/J2Iz0fkyxtj5Wm9De+Ozwv5nZYJECAwjMDdB4rv//P47zCKziFAgAABAgQIECBAgED7BMaZANwdnCmRt9AyESeMMgH4vKj/5lyjtuaW0/bdufW0uKe0nlZfH5Ee610T8bGIt0eku/zS+isjXheRyj0R73xsyT8ECBBYpMBdDxUfAb7VBCCLFHU6AQIECBAgQIAAAQIE2iEwzgTg3iCcLwGY3qO3JaKf9DsUy6cjRlleG5W/ZpYLPDe2p8iXPfmV3vLfxudPRHwwYkNESgCWS0r+vSyi+Bd7+SjrBAgQGEBgamoqu7v0DsBbTQAygJxDCBAgQIAAAQIECBAg0H6BcSYAbxyQd3Mc979FvC0ivS/vH0fcHVH38pFo4O0R6W7AlOi7LuJ8xL0RfxTxf0WMOpkZl1AIEOiCQJr998TZi4Wu3rIj/f8HhQABAgQIECBAgAABAgS6LtCEt8MfiUH6zYh05922iI9GpKTgqMruqDjdcThozNWOB2LnGyJuiUh3M6Z2f2/EOyIk/wJBIUCgGoHy+//Wr57Mrt24uprK1UKAAAECBAgQIECAAAECjRZoQgKwD5zu+vsPETdGvDFCIUCAAIGeQHkG4Fu2r88mJtL/x1AIECBAgAABAgQIECBAoOsCTUoAprH6q96A/WjXB07/CRAgkBcwAUhewzIBAgQIECBAgAABAgQI5AWalgA82Wv8DflOWCZAgEDXBe4+UJxPyPv/uv4Tof8ECBAgQIAAAQIECBC4ItC0BODTe02/cKULlggQINBtgfMXL2ffeKT//0cetzADcLd/JvSeAAECBAgQIECAAAECeYEmJQCfHA1/a8RUxN9FKAQIECAQAvcdOpldvJx+NV4pT413ACoECBAgQIAAAQIECBAgQCAJTI6R4dUDXDslKNPMuc+M+JGItRHpr9zfjlAIECBAIATKj/+m2X83rlnBhgABAgQIECBAgAABAgQIPCYwzgTgnmhB8ZaVuQelP51lmgn4Q3Mfai8BAgS6I/C1h8rv/3P3X3dGX08JECBAgAABAgQIECAwv8A4E4Cpdf2k3nwtPRoHfDriNyM+Nt/B9hMgQKBLAncfOF7o7q1P2lBYt0KAAAECBAgQIECAAAEC3RYYZwIwvdNvvnI5Dki3tqQEoEKAAAECMwiUHwE2AcgMSDYRIECAAAECBAgQIECgwwLjTAA+0GF3XSdAgEAlAsfOXMgePHa2UNctOzwCXACxQoAAAQIECBAgQIAAgY4LNGkW4I4Ple4TIEBgusA9B4vv/5tcNpE9Zeu66QfaQoAAAQIECBAgQIAAAQKdFZAA7OzQ6zgBAm0QuOuh4vv/brpmXbZy0q/2NoytPhAgQIAAAQIECBAgQKAqAX8lViWpHgIECIxB4K4DxTsAb32Sx3/HMAwuSYAAAQIECBAgQIAAgVoLLMU7AO8bgcBU1HnTCOpVJQECBBolUJ4AxPv/GjV8GkuAAAECBAgQIECAAIElEViKBOCNI+hJSgAqBAgQ6LTA1NRUdnfpHYBmAO70j4TOEyBAgAABAgQIECBAYEaBpUgA/s6MV7aRAAECBBYlkGb/PXH2YqGOW3ZsKKxbIUCAAAECBAgQIECAAAECS5EA/NMe8yfi8yRyAgQIEKhGoDwByPrVk9m1G1dXU7laCBAgQIAAAQIECBAgQKA1AksxCciHQ+tPIm4oqb0/1t8X8aTSdqsECBAgMIDAtAlAdqzPJiYmBjjTIQQIECBAgAABAgQIECDQJYGlSAAmz5n+It0d21NsjlAIECBAYIECJgBZIJjDCRAgQIAAAQIECBAg0FGBpUgAnuvZruuosW4TIEBgJALTE4De/zcSaJUSIECAAAECBAgQIECg4QJLkQDc3zN6fsOtNJ8AAQK1ETh/8XL2jUeKr1U1A3BthkdDCBAgQIAAAQIECBAgUCuBpZgE5OPR45+JeHvEsyLuibgQ0S//IhYe7q8s4PNtCzjWoQQIEGiVQEr+Xbw8VejTU7evL6xbIUCAAAECBAgQIECAAAECSWApEoC/Gtf50YgtET8WkS/p3YA/l9+wgGUJwAVgOZQAgXYJlB//TbP/blyzol2d1BsCBAgQIECAAAECBAgQqERgKR4B/la09Hsi3htxf0S6+y8l/vq3rqTlYSJOUwgQINBNgWkzAD/J+/+6+ZOg1wQIECBAgAABAgQIEJhfYCnuAEytSEnA15WacznWUxLwtoivlvZZJUCAAIE5BO4+cLyw95YdHv8tgFghQIAAAQIECBAgQIAAgScEluIOwCcuZoEAAQIEqhGYdgegBGA1sGohQIAAAQIECBAgQIBACwWW6g7Ameju6G3cN9NO2wgQIEBgZoFjpy9kDx07W9jpDsAChxUCBAgQIECAAAECBAgQyAmMMwH4O7l2WCRAgACBAQXuPniicOTksonsKVvXFbZZIUCAAAECBAgQIECAAAECfQGPAPclfBIgQKAhAuX3/928bV22ctKv84YMn2YSIECAAAECBAgQIEBgyQX8xbjk5C5IgACBxQmU3//n8d/FeTqbAAECBAgQIECAAAECbReQAGz7COsfAQKtE5AAbN2Q6hABAgQIECBAgAABAgRGKiABOFJelRMgQKBagampqeyeA8V3AN5qBuBqkdVGgAABAgQIECBAgACBlglIALZsQHWHAIF2C+w/eiY7ce5ioZO37NhQWLdCgAABAgQIECBAgAABAgTyAhKAeQ3LBAgQqLnA3aW7/9avnsyu3bi65q3WPAIECBAgQIAAAQIECBAYp4AE4Dj1XZsAAQILFCi//y89/jsxMbHAWhxOgAABAgQIECBAgAABAl0SkADs0mjrKwECjRcoJwDNANz4IdUBAgQIECBAgAABAgQIjFxAAnDkxC5AgACB6gTKE4B4/191tmoiQIAAAQIECBAgQIBAWwUkANs6svpFgEDrBC5eupx989CpQr9u2b6+sG6FAAECBAgQIECAAAECBAiUBSQAyyLWCRAgUFOBfUfOZOcjCZgvN29bl1+1TIAAAQIECBAgQIAAAQIEpglIAE4jsYEAAQL1FPjGIycLDdu8dkV29VUrC9usECBAgAABAgQIECBAgACBsoAEYFnEOgECBGoqUE4APuUad//VdKg0iwABAgQIECBAgAABArUSkACs1XBoDAECBGYX+MbDxff/3XTNVbMfbA8BAgQIECBAgAABAgQIEOgJSAD6USBAgEBDBO47VHwE+CZ3ADZk5DSTAAECBAgQIECAAAEC4xWQAByvv6sTIEBgYIFvPFK+A9AjwAPjOZAAAQIECBAgQIAAAQIdFpAA7PDg6zoBAs0ROHzqfJYiX24yA3CewzIBAgQIECBAgAABAgQIzCIgATgLjM0ECBCok8B9pRmAVyyfyK7fvKZOTdQWAgQIECBAgAABAgQIEKipgARgTQdGswgQIJAXKM8AvGvLVdnkcr/C80aWCRAgQIAAAQIECBAgQGBmAX89zuxiKwECBGolMP39f2YArtUAaQwBAgQIECBAgAABAgRqLCABWOPB0TQCBAj0BcqPAJsBuC/jkwABAgQIECBAgAABAgTmE5AAnE/IfgIECNRAYPodgGYArsGwaAIBAgQIECBAgAABAgQaISAB2Ihh0kgCBLoscO7ipWzv4dMFAjMAFzisECBAgAABAgQIECBAgMAcAhKAc+DYRYAAgToI7H30dHbp8lShKU+5xjsACyBWCBAgQIAAAQIECBAgQGBWAQnAWWnsIECAQD0EyjMAX7N+VbZh9Yp6NE4rCBAgQIAAAQIECBAgQKD2AhKAtR8iDSRAoOsC09//5+6/rv9M6D8BAgQIECBAgAABAgQWIiABuBAtxxIgQGAMAuU7AM0APIZBcEkCBAgQIECAAAECBAg0WEACsMGDp+kECHRDYPodgGYA7sbI6yUBAgQIECBAgAABAgSqEZAArMZRLQQIEBiJwNTUVHbfwycLdZsBuMBhhQABAgQIECBAgAABAgTmEZAAnAfIbgIECIxT4JET57IT5y4WmvCUrd4BWACxQoAAAQIECBAgQIAAAQJzCkgAzsljJwECBMYrcO8jxbv/Vk0uy3ZuWjPeRrk6AQIECBAgQIAAAQIECDRKQAKwUcOlsQQIdE3gvkdOFbr8lGvWZcuWTRS2WSFAgAABAgQIECBAgAABAnMJSADOpWMfAQIExiwwfQZgj/+OeUhcngABAgQIECBAgAABAo0TkABs3JBpMAECXRIwA3CXRltfCRAgQIAAAQIECBAgMBoBCcDRuKqVAAEClQh8wwzAlTiqhAABAgQIECBAgAABAl0WkADs8ujrOwECtRY4c/5Stv/omUIbzQBc4LBCgAABAgQIECBAgAABAgMISAAOgOQQAgQIjEPgm4eKE4CkNjzlGu8AHMdYuCYBAgQIECBAgAABAgSaLCAB2OTR03YCBFotUJ4AZOemNdnalZOt7rPOESBAgAABAgQIECBAgED1AhKA1ZuqkQABApUIlBOA7v6rhFUlBAgQIECAAAECBAgQ6JyABGDnhlyHCRBoioAZgJsyUtpJgAABAgQIECBAgACBegtIANZ7fLSOAIEOC0ybAdj7/zr806DrBAgQIECAAAECBAgQGF5AAnB4O2cSIEBgZAKXL09l9x06Waj/pmvWFdatECBAgAABAgQIECBAgACBQQQkAAdRat4xO6PJPx/xsYi9EecjDkT814h/EDFT2R0bp+aJj890om0ECFQv8NDxs9nZC5cLFd+0TQKwAGKFAAECBAgQIECAAAECBAYSMJ3kQEyNO+hfR4vfHPGNiJQEfCTi2yJe0YtXxeeHIvLl72LlV/Ibcss/FsvfEfGXuW0WCRAYoUD58d91qyazbetXjfCKqiZAgAABAgQIECBAgACBtgpIALZzZL8Q3XphxKdK3Xt+rKe7+H4r4sMR5yL6JSUAU5TLytjwryIuRvxOead1AgRGI1CeAfimeP/fxMTEaC6mVgIECBAgQIAAAQIECBBotYBHgNs5vH8S3Son/1JPPxPxyYjNEbdFDFJeEQdtifiziIODnOAYAgQWLzA9Aejx38WrqoEAAQIECBAgQIAAAQLdFJAAnD7u22LTD0e8LeKjEYci+u/G2xPLCym74uB3RtwVcSricMQXI94UsTZiHOVC76Lpjr5Bymt7B713kIMdQ4BANQLfeDj9yrhSnmIG4CsYlggQIECAAAECBAgQIEBgQQIeAZ7OVdVdbi+Pqj8YsSF3iZT0e2YvUmLtZRH3RixVuSEu9AMRD0XcOcBFUwLzxRH7Iv77AMc7hACBigTMAFwRpGoIECBAgAABAgQIECBAIHMH4Nw/BHtjd5pEY6Hl6XFCmmQjJf9ORrwl4jkRKZn2nohUnhrx5xHr08oSlBVxjf8SkWYRSBOEXIqYr9wRB6SfkT0RgxwfhykECCxW4MTZC9nB4/lXdGaZGYAXq+p8AgQIECBAgAABAgQIdFdgsrtdn7Xn6dHf9JhuinQ34I0R34xYSHl3HLwmIj1m+5KIz0f0yydi4esR74hIScA3Rrw1olzSo8MLmfIzXTPVO1PpJ/G+L3amBGRKBM5X0jkpAZgef37/fAfbT4BAdQL3PVJ8/HfZRJbt2jKutwZU1y81ESBAgAABAgQIECBAgMB4BCQAp7v/8vRNC9ryrDg6zbabyvsi8sm/xzbGPym5l5JrT4t4fcSvRVyIyJefjZWr8hvmWf7j2D9TAjAl8lIC71URH4z45xGDlB+Ig26I+HjEQhOgg9TvGAIEZhEoTwByw9Vrs1WTy2c52mYCBAgQIECAAAECBAgQIDC3QEoOKdUKpFlz++UD/YXS5+VY/93etk3x+aLS/rSapvyM+34Gjr+OY8sljW9qw2sifj9id0S69iDltb2DTP4xiJZjCFQoUE4A3nRN+nWgECBAgAABAgQIECBAgACB4QQkAIdzm+us5/V2pmf4vjzHgZ/K7XtubrmqxX7y79VRYXof4T+NuDRg5VviuB+JSLMW/+mA5ziMAIGKBMwAXBGkaggQIECAAAECBAgQIEDgMQEJwOp/ENJjvamk2X3TOwBnK3fldvTPyW1a1GIa1/dHpOTfH0X8VMSgyb849LFk4cr4TI8MF2ciSHsVAgRGKmAG4JHyqpwAAQIECBAgQIAAAQKdE5jsXI9H2+HVUf3W3iX2zXOpI7E/3SWY3vN3/TzHLnT3L8UJ6bHfNAPxPRH/NqJcPhwb/q68sbf+073P986yf5DN181z0I559ttNoJMCFy9dzu4/dLrQdzMAFzisECBAgAABAgQIECBAgMACBSQAFwg2z+Hrc/tT8m2+0k8AVv2Crxt7F071vmWWRtwf22dKAKZJTL4z4gsRd0YMW7417InOI9BlgX1HzmTnIwmYL94BmNewTIAAAQIECBAgQIAAAQILFZAAXKjY3MenOwD75Xx/YY7P/uO1a+Y4Zphdu+OkFMOUlPhLk48oBAiMQaA8AcjmtSuyq69KT+QrBAgQIECAAAECBAgQIEBgOAEJwOHcZjvrbG7HIH+xr+odfyZ3XlsW53usOT0C/MW2dFY/CFQlUE4APsUMwFXRqocAAQIECBAgQIAAAQKdFZAArHboT+SqG+Sx3vT+v1QGeVz48SOb8+++5jRVSwmw4ReVAABAAElEQVTUR+C+R9KbAa6Um67p/5q4ss0SAQIECBAgQIAAAQIECBBYiIBZgBeiNf+x6Q7AR3uHzTcJxuY4rv+XvfflzW/rCAKdECjfAej9f50Ydp0kQIAAAQIECBAgQIDASAUkAKvn/Wqvypvjc647LG/NXfpruWWLBAh0WOAb0+4AHORm4g6D6ToBAgQIECBAgAABAgQIzCsgATgv0YIP+GzvjHR33zPmOPsFuX2fyy1bJECgowKHT53PUuTLTdskAPMelgkQIECAAAECBAgQIEBg4QISgAs3m++MD+cOuCO3nF9M7q/ubTgan5/M77RMgEA3Be57pPg60BXLJ7LrN1c9SXg3bfWaAAECBAgQIECAAAECXRaQAKx+9L8QVX6mV+1Px+ezZ7jEG2Pb03rb3x2fF2Y4xiYCBDomUH7/364tV2WTy/2a7tiPge4SIECAAAECBAgQIECgcoG53lFX+cUaUuHzop3p/X39srW/EJ9p++7celrcU1pPq6+PSI/1plt3Phbx9oh0l19af2XE6yJSuSfinY8t+YcAgc4LmAG48z8CAAgQIECAAAECBAgQIDASAQnA6ayvjU2vmb75sS3PjX9T5Mue/Epv+W/j8yciPhixISIlAMslJf9eFnGivMM6AQLdFLj/0VOFjj95q/f/FUCsECBAgAABAgQIECBAgMBQAp4tG4ptoJM+EkfdHvGuiJTsOx2R3vf3pYg3Rzw94t4IhQABAo8J7D18piCxa8vawroVAgQIECBAgAABAgQIECAwjIA7AKer7Y5NKaooD0Qlb+hFFfWpgwCBlgpMTU1le0t3AN5wtQRgS4dbtwgQIECAAAECBAgQILCkAu4AXFJuFyNAgMDMAo+eOp+dOn+psFMCsMBhhQABAgQIECBAgAABAgSGFJAAHBLOaQQIEKhS4IFH01sCrpQVyyeyazeleYMUAgQIECBAgAABAgQIECCwOAEJwMX5OZsAAQKVCOw9XJwA5LrNa7PlyyYqqVslBAgQIECAAAECBAgQINBtAQnAbo+/3hMgUBOB8h2AHv+tycBoBgECBAgQIECAAAECBFogIAHYgkHUBQIEmi+wt/QIsBmAmz+mekCAAAECBAgQIECAAIG6CEgA1mUktIMAgU4LPHC4+A5AdwB2+sdB5wkQIECAAAECBAgQIFCpgARgpZwqI0CAwHAC5UeAd225ariKnEWAAAECBAgQIECAAAECBEoCEoAlEKsECBBYaoFT5y5mh06eK1zWI8AFDisECBAgQIAAAQIECBAgsAgBCcBF4DmVAAECVQh860jx8d9Up0eAq5BVBwECBAgQIECAAAECBAgkAQlAPwcECBAYs0D58d/tG1Zlq1csH3OrXJ4AAQIECBAgQIAAAQIE2iIgAdiWkdQPAgQaK1CeAdjdf40dSg0nQIAAAQIECBAgQIBALQUkAGs5LBpFgECXBB44fKrQ3RuuNgFIAcQKAQIECBAgQIAAAQIECCxKQAJwUXxOJkCAwOIFyo8AmwBk8aZqIECAAAECBAgQIECAAIErAhKAVywsESBAYCwCew8XJwGRABzLMLgoAQIECBAgQIAAAQIEWisgAdjaodUxAgSaIHDx0uVs/5EzhaZ6B2CBwwoBAgQIECBAgAABAgQILFJAAnCRgE4nQIDAYgQePHo2u3h5qlDFri3eAVgAsUKAAAECBAgQIECAAAECixKQAFwUn5MJECCwOIHyBCDrV01mm9euWFylziZAgAABAgQIECBAgAABAjkBCcAchkUCBAgstUB5ApAbtqzNJiYmlroZrkeAAAECBAgQIECAAAECLRaQAGzx4OoaAQL1FzABSP3HSAsJECBAgAABAgQIECDQdAEJwKaPoPYTINBogb2PFmcAvuFq7/9r9IBqPAECBAgQIECAAAECBGooIAFYw0HRJAIEuiPwwOFiAnBXPAKsECBAgAABAgQIECBAgACBKgUkAKvUVBcBAgQWIDA1NZXtffRU4YxdV0sAFkCsECBAgAABAgQIECBAgMCiBSQAF02oAgIECAwn8Oip89mp85cKJ18vAVjwsEKAAAECBAgQIECAAAECixeQAFy8oRoIECAwlEB5BuAVyyeyazetGaouJxEgQIAAAQIECBAgQIAAgdkEJABnk7GdAAECIxbYe7j4+O91m9dmy5dNjPiqqidAgAABAgQIECBAgACBrglIAHZtxPWXAIHaCJTvALzB47+1GRsNIUCAAAECBAgQIECAQJsEJADbNJr6QoBAowT2PmoG4EYNmMYSIECAAAECBAgQIECgoQISgA0dOM0mQKD5Ag8cLiYA3QHY/DHVAwIECBAgQIAAAQIECNRRQAKwjqOiTQQIdEKg/Ajwri1XdaLfOkmAAAECBAgQIECAAAECSysgAbi03q5GgACBxwROn7+YHTp5rqCxa8vawroVAgQIECBAgAABAgQIECBQhYAEYBWK6iBAgMACBfaWHv9Np3sEeIGIDidAgAABAgQIECBAgACBgQQkAAdichABAgSqFSg//rt9w6ps9Yrl1V5EbQQIECBAgAABAgQIECBAIAQkAP0YECBAYAwC02YAvtr7/8YwDC5JgAABAgQIECBAgACBTghIAHZimHWSAIG6CTxw+FShSddf7f1/BRArBAgQIECAAAECBAgQIFCZgARgZZQqIkCAwOAC5UeATQAyuJ0jCRAgQIAAAQIECBAgQGBhAhKAC/NyNAECBCoRKE8CIgFYCatKCBAgQIAAAQIECBAgQGAGAQnAGVBsIkCAwCgFLl66nO0/cqZwCTMAFzisECBAgAABAgQIECBAgECFAhKAFWKqigABAoMIPHj0bHbx8lTh0F1bTAJSALFCgAABAgQIECBAgAABApUJSABWRqkiAgQIDCZQngBk/arJbPPaFYOd7CgCBAgQIECAAAECBAgQILBAAQnABYI5nAABAosVKL//74Yta7OJiYnFVut8AgQIECBAgAABAgQIECAwo4AE4IwsNhIgQGB0AnsfPV2o3AQgBQ4rBAgQIECAAAECBAgQIFCxgARgxaCqI0CAwHwCD5QSgDdc7f1/85nZT4AAAQIECBAgQIAAAQLDC0gADm/nTAIECAwl8MBhdwAOBeckAgQIECBAgAABAgQIEBhKQAJwKDYnESBAYDiBqampbO+jpwon77p6bWHdCgECBAgQIECAAAECBAgQqFJAArBKTXURIEBgHoFHT53PTp2/VDgqTQKiECBAgAABAgQIECBAgACBUQlIAI5KVr0ECBCYQaD8/r8VyyeyJ21cM8ORNhEgQIAAAQIECBAgQIAAgWoEJACrcVQLAQIEBhLYe7j4+O91m9dmy5dNDHSugwgQIECAAAECBAgQIECAwDACEoDDqDmHAAECQwqU7wC8wfv/hpR0GgECBAgQIECAAAECBAgMKiABOKiU4wgQIFCBwN5HzQBcAaMqCBAgQIAAAQIECBAgQGABAhKAC8ByKAECBBYr8MDhYgLQHYCLFXU+AQIECBAgQIAAAQIECMwnIAE4n5D9BAgQqFBgbykBuGvLVRXWrioCBAgQIECAAAECBAgQIDBdQAJwuoktBAgQGInA6fMXs0dOnCvUvWvL2sK6FQIECBAgQIAAAQIECBAgULWABGDVouojQIDALALlu//SYR4BngXLZgIECBAgQIAAAQIECBCoTEACsDJKFREgQGBugfIMwNs3rMpWr1g+90n2EiBAgAABAgQIECBAgACBRQpIAC4S0OkECBAYVGDaDMBXe//foHaOI0CAAAECBAgQIECAAIHhBSQAh7dzJgECBBYk8MDhU4Xjb/D+v4KHFQIECBAgQIAAAQIECBAYjYAE4Ghc1UqAAIFpAuVHgL3/bxqRDQQIECBAgAABAgQIECAwAgEJwBGgqpIAAQIzCZQnATED8ExKthEgQIAAAQIECBAgQIBA1QISgFWLqo8AAQIzCFy8dDnbf+RMYY87AAscVggQIECAAAECBAgQIEBgRAISgCOCVS0BAgTyAg8ePZtdvDyV35Tt2mISkAKIFQIECBAgQIAAAQIECBAYiYAE4EhYVUqAAIGiQPnx3/WrJrPNa1cUD7JGgAABAgQIECBAgAABAgRGICABOAJUVRIgQKAsMNMMwBMTE+XDrBMgQIAAAQIECBAgQIAAgcoFJAArJ1UhAQIEpgvsffR0YaMJQAocVggQIECAAAECBAgQIEBghAISgCPEVTUBAgT6At86UkwAXn/12v4unwQIECBAgAABAgQIECBAYKQCEoAj5VU5AQIEHhfYH5OA5Mt1myUA8x6WCRAgQIAAAQIECBAgQGB0AhKAo7NVMwECBJ4Q2H/kzBPLaWHnptWFdSsECBAgQIAAAQIECBAgQGBUAhKAo5JVLwECBHoCZy9cyg6dPFfw2LnJHYAFECsECBAgQIAAAQIECBAgMDIBCcCR0aqYAAECjws8dKz4+G/aeq07AP14ECBAgAABAgQIECBAgMASCUgALhG0yxAg0F2BB48WH//dsHoyW796RXdB9JwAAQIECBAgQIAAAQIEllRAAnBJuV2MAIEuCpTf/3ftpjVdZNBnAgQIECBAgAABAgQIEBiTgATgmOBHeNk0s8BvRHw64sGI9OzhgYjPRdwRMdNtR8OcE1UpBAgMIrC/dAfgdZslAAdxcwwBAgQIECBAgAABAgQIVCMgAViNY51qWReN+bmIqYg/j0jJwD+N2Bnx/og/iyiP+zDnRDUKAQKDCJQTgO4AHETNMQQIECBAgAABAgQIECBQlcBkVRWppzYCh6MlGyPOl1qUxvp/RLwk4qURKTnYL8Oc0z/XJwEC8wiU3wG40yPA84jZTYAAAQIECBAgQIAAAQJVCpTvBKuybnWNR+ByXLac/EstuRiR7gRM5ebHP574d5hznjjZAgECcwu4A3BuH3sJECBAgAABAgQIECBAYLQCEoBF322x+sMRb4v4aMShiPQobYo9EQspu+Lgd0bcFXEqIt1l98WIN0WsjVjqksb6B3sX/fsBLz7MOQNW7TAC3RC4fHkqe+hoehXnlbLTOwCvYFgiQIAAAQIECBAgQIAAgZELeAS4SHywuDr02svjzA9GbMjVkJJ+z+zFa+PzZRH3RoyqrIyKfzFiImJLxIsjbo34QMTHI2Yqw5wzUz22ESDQEzh08lx2/lK6yfZK8QjwFQtLBAgQIECAAAECBAgQIDB6AQnA2Y33xq509156Z95CytPj4A9FrIk4GfHvIj7ZW39lfP5MxFMj0jv4UkLwRMQoSkrm/XKu4qlY/vWIf5PbVl4c5pxyHdYJEMgJlB//XbF8Irtm3arcERYJECBAgAABAgQIECBAgMBoBSQAi77p0d/0mG6KdDfgjRHfjFhIeXccnJJ/6Z17KXn4+Yh++UQsfD3iHREpCfjGiLdGlEt6dHghGYJ0zVRvvqTkY7r7Lz3Ge21Euivx7RHPjvihiOMR5TLMOeU6rBMgkBN4sPT475M2rsmWLUv/aSoECBAgQIAAAQIECBAgQGBpBCQAi875O+aKewZbe1Yc9vzeoe+Lz3zyr19DSu7dEfG0iNdH/FrEhYh8+dlYuSq/YZ7lP4795QRg/5T07OG+iN+KSO80/MOIt0S8OWK2Msw5s9VlO4FOC+w/errQ/2s3rS6sWyFAgAABAgQIECBAgAABAqMWSHeHKdUJvCJX1Qdyy/nFlFz73d6GTfH5ovzO3vK6+Ey3CA0af907b76Pj/UOeOF8B+b2D3NO7nSLBLotUL4DcOem9DpQhQABAgQIECBAgAABAgQILJ2ABGC11s/rVZdm/f3yHFV/KrfvubnlUS+mR4FTKd9x+PjWmf8d5pyZa7KVQAcF9h05U+j1TncAFjysECBAgAABAgQIECBAgMDoBSQAqzVOj/Wmkmb3Te8AnK2kyUX6pX9Of32xn98eFcx0i1Ha9hu9yv+idJFhzilVYZUAgZkEHjxaSgBuTq8IVQgQIECAAAECBAgQIECAwNIJeAdgddbpxV5be9Wld+7NVY7EznSXYHrP3/VzHTjEvh+Pc94Q8dmI+yPSZB87I14asSXiMxHvisiXYc7Jnz/T8nUzbcxt25FbtkigtQLlWYCv3SQB2NrB1jECBAgQIECAAAECBAjUVEACsLqBWZ+rKs2mO1/pJwDT+/6qLH8WlaXHdp8TkWb8TfUfi/hKxB9EvD+ifHfiMOdENXOWb825104CHRA4ee5iduxM8Yn7nRKAHRh5XSRAgAABAgQIECBAgEC9BCQAqxuP/NSe5weo9lzvmKpvB/pS1JtiIWWYcxZSv2MJdFKg/PhvQnAHYCd/FHSaAAECBAgQIECAAAECYxWQAKyO/2yuqpW55dkWV/V2FF8QNtvRzds+36PN6RHgLzavW1pMYHCB8uO/W9etzFavWD54BY4kQIAAAQIECBAgQIAAAQIVCEgAVoDYq+JErqr02O18Jb3/L5VBHhd+/Mhm/buvWc3VWgLVC+wvzQDs7r/qjdVIgAABAgQIECBAgAABAvMLmAV4fqNBj0h3AD7aO3i+CTA2x3H9BKB35Q0q7DgCDRMoPwLs/X8NG0DNJUCAAAECBAgQIECAQEsEJACrHciv9qq7OT7nurvy1txlv5ZbtkiAQIsEyo8AuwOwRYOrKwQIECBAgAABAgQIEGiQgARgtYP12V516e6+Z8xR9Qty+z6XW7ZIgECLBNwB2KLB1BUCBAgQIECAAAECBAg0WEACsNrB+3Cuujtyy/nFZP7q3oaj8fnJ/E7LBAi0R8A7ANszlnpCgAABAgQIECBAgACBJgtIAFY7el+I6j7Tq/Kn4/PZM1T/xtj2tN72d8fnhRmOsYkAgYYLXLx0OTtwPD85eJZdt3lNw3ul+QQIECBAgAABAgQIECDQRIG53lPXxP4sts3PiwrS+/v6ZWt/IT7T9t259bS4p7SeVl8fkR7rTX/pfyzi7RHpLr+0/sqI10Wkck/EOx9b8g8BAq0TSMm/y1PFbnkHYNHDGgECBAgQIECAAAECBAgsjYAEYNH5tbH6muKmJ9aeG0sp8mVPfqW3/Lfx+RMRH4zYEJESgOWSkn8vizhR3mGdAIF2CDx4tHj33+oVy7LNa1e0o3N6QYAAAQIECBAgQIAAAQKNEvAI8GiG6yNR7e0R74pIyb7TEel9f1+KeHPE0yPujVAIEGipwEwTgExMTLS0t7pFgAABAgQIECBAgAABAnUWcAdgcXR2x2qKKsoDUckbelFFfeogQKBBAvuPnim01uO/BQ4rBAgQIECAAAECBAgQILCEAu4AXEJslyJAoDsC5QSgCUC6M/Z6SoAAAQIECBAgQIAAgboJSADWbUS0hwCBVgjsP1K6A3CjGYBbMbA6QYAAAQIECBAgQIAAgQYKSAA2cNA0mQCB+gtMewfgZgnA+o+aFhIgQIAAAQIECBAgQKCdAhKA7RxXvSJAYIwCU1NTWfkRYO8AHOOAuDQBAgQIECBAgAABAgQ6LiAB2PEfAN0nQKB6gWNnLmSnz18qVLxzkzsACyBWCBAgQIAAAQIECBAgQGDJBCQAl4zahQgQ6IrAvtL7/5ZNZNmOjau70n39JECAAAECBAgQIECAAIGaCUgA1mxANIcAgeYLlN//t33D6mzFcr9umz+yekCAAAECBAgQIECAAIFmCviLtJnjptUECNRYwPv/ajw4mkaAAAECBAgQIECAAIEOCkgAdnDQdZkAgdEKlO8A9P6/0XqrnQABAgQIECBAgAABAgTmFpAAnNvHXgIECCxYwB2ACyZzAgECBAgQIECAAAECBAiMUEACcIS4qiZAoJsC+4+eLXR852YzABdArBAgQIAAAQIECBAgQIDAkgpIAC4pt4sRINAFgf2lWYB3bjIDcBfGXR8JECBAgAABAgQIECBQVwEJwLqOjHYRINBIgbMXLmWHTp4rtH3nprWFdSsECBAgQIAAAQIECBAgQGApBSQAl1LbtQgQaL3AQ8eKj/+mDl/rDsDWj7sOEiBAgAABAgQIECBAoM4CEoB1Hh1tI0CgcQLlGYA3rJ7M1q9e0bh+aDABAgQIECBAgAABAgQItEdAArA9Y6knBAjUQMAMwDUYBE0gQIAAAQIECBAgQIAAgYKABGCBwwoBAgQWJzB9AhAzAC9O1NkECBAgQIAAAQIECBAgsFgBCcDFCjqfAAECOYHyI8A7N0sA5ngsEiBAgAABAgQIECBAgMAYBCQAx4DukgQItFfAI8DtHVs9I0CAAAECBAgQIECAQFMFJACbOnLaTYBALQWm3QG4yR2AtRwojSJAgAABAgQIECBAgECHBCQAOzTYukqAwGgFLl+eyh48erZwkWslAAseVggQIECAAAECBAgQIEBg6QUkAJfe3BUJEGipwKFT57Lzly4XeneddwAWPKwQIECAAAECBAgQIECAwNILSAAuvbkrEiDQUoHyDMArlk9k16xb1dLe6hYBAgQIECBAgAABAgQINEVAArApI6WdBAjUXqD8+O+TNq7Jli2bqH27NZAAAQIECBAgQIAAAQIE2i0gAdju8dU7AgSWUGD/0dOFq127aXVh3QoBAgQIECBAgAABAgQIEBiHgATgONRdkwCBVgqU7wDcuWltK/upUwQIECBAgAABAgQIECDQLAEJwGaNl9YSIFBjgX1HzhRat9MdgAUPKwQIECBAgAABAgQIECAwHgEJwPG4uyoBAi0UePBoKQFoBuAWjrIuESBAgAABAgQIECBAoHkCEoDNGzMtJkCgpgL7SwnAazetqWlLNYsAAQIECBAgQIAAAQIEuiQgAdil0dZXAgRGJnDy3MXs2JkLhfp3SgAWPKwQIECAAAECBAgQIECAwHgEJADH4+6qBAi0TKD8+G/qnjsAWzbIukOAAAECBAgQIECAAIGGCkgANnTgNJsAgXoJlB//3bpuZbZ6xfJ6NVJrCBAgQIAAAQIECBAgQKCTAhKAnRx2nSZAoGqB/aUZgN39V7Ww+ggQIECAAAECBAgQIEBgWAEJwGHlnEeAAIGcQPkR4Gs3mgAkx2ORAAECBAgQIECAAAECBMYoIAE4RnyXJkCgPQLlBODOzRKA7RldPSFAgAABAgQIECBAgECzBSQAmz1+Wk+AQE0Eyu8A9AhwTQZGMwgQIECAAAECBAgQIEAgkwD0Q0CAAIEKBB48erZQy85N7gAsgFghQIAAAQIECBAgQIAAgbEJSACOjd6FCRBoi8DFS5ezA8clANsynvpBgAABAgQIECBAgACBtglIALZtRPWHAIElFzh44lx26fJU4breAVjgsEKAAAECBAgQIECAAAECYxSQABwjvksTINAOgf1HzhQ6snrFsmzz2hWFbVYIECBAgAABAgQIECBAgMC4BCQAxyXvugQItEZg2gzA8f6/iYmJ1vRPRwgQIECAAAECBAgQIECg2QISgM0eP60nQKAGAmYArsEgaAIBAgQIECBAgAABAgQIzCogATgrjR0ECBAYTKCcALxusxmAB5NzFAECBAgQIECAAAECBAgshYAE4FIouwYBAq0WKL8D8NqNEoCtHnCdI0CAAAECBAgQIECAQMMEJAAbNmCaS4BA/QSmvQPQHYD1GyQtIkCAAAECBAgQIECAQIcFJAA7PPi6ToDA4gWmpqay8iPA18YkIAoBAgQIECBAgAABAgQIEKiLgARgXUZCOwgQaKTAsTMXstPnLxXavlMCsOBhhQABAgQIECBAgAABAgTGKyABOF5/VydAoOECDx49W+jBxESW7di4urDNCgECBAgQIECAAAECBAgQGKeABOA49V2bAIHGCxw8UUwAbrlqVbZiuV+tjR9YHSBAgAABAgQIECBAgECLBPyV2qLB1BUCBJZe4OCxYgJwx8ZVS98IVyRAgAABAgQIECBAgAABAnMISADOgWMXAQIE5hM4ePxc4ZDt6z3+WwCxQoAAAQIECBAgQIAAAQJjF5AAHPsQaAABAk0WOHC8eAfgdu//a/JwajsBAgQIECBAgAABAgRaKSAB2Mph1SkCBJZK4OFyAtAdgEtF7zoECBAgQIAAAQIECBAgMKCABOCAUA4jQIDATALlOwC9A3AmJdsIECBAgAABAgQIECBAYJwCEoDj1HdtAgQaL1B+B+C2Dd4B2PhB1QECBAgQIECAAAECBAi0TEACsGUDqjsECCydwIVLl7NHTxUnAdkhAbh0A+BKBAgQIECAAAECBAgQIDCQgATgQEwOIkCAwHSBR06cy6amitu3SwAWQawRIECAAAECBAgQIECAwNgFJADHPgQaQIBAUwXK7/9buXxZtnntiqZ2R7sJECBAgAABAgQIECBAoKUCEoAtHVjdIkBg9ALlGYC3bViVTUxMjP7CrkCAAAECBAgQIECAAAECBBYgIAG4ACyHEiBAIC9w4NjZ/Grm/X8FDisECBAgQIAAAQIECBAgUBMBCcCaDIRmECDQPIGD8Q7AfPH+v7yGZQIECBAgQIAAAQIECBCoi4AEYF1GQjsIEGicwMHSHYASgI0bQg0mQIAAAQIECBAgQIBAJwQkADsxzDpJgMAoBA6eKD4CvD3eAagQIECAAAECBAgQIECAAIG6CUgA1m1EtIcAgcYITHsH4MbVjWm7hhIgQIAAAQIECBAgQIBAdwQkALsz1npKgEDFAg8fL74DcNt6CcCKiVVHgAABAgQIECBAgAABAhUISABWgKgKAgS6J3Dq3MXsRES+7HAHYJ7DMgECBAgQIECAAAECBAjUREACsCYDoRkECDRL4ODx4vv/Uuu9A7BZY6i1BAgQIECAAAECBAgQ6IqABGBXRlo/CRCoVOBAKQG4fvVktnblZKXXUBkBAgQIECBAgAABAgQIEKhCQAKwCkV1ECDQOYHy+/+2b/D+v879EOgwAQIECBAgQIAAAQIEGiIgAdiQgdJMAgTqJVC+A3CHBGC9BkhrCBAgQIAAAQIECBAgQOAJAQnAJyhas5BuQ/qNiE9HPBiRXlR2IOJzEXdErIiYqaSfhX8V8TcRpyOOR6Q6/nGEQoBASaD8DsBtG1aVjrBKgAABAgQIECBAgAABAgTqISABWI9xqLIV66Kyn4uYivjziJQM/NOInRHvj/iziPK4T8S2P4z4jxEbIt4X8QcRt0T8t4iUGFQIEMgJlBOA7gDM4VgkQIAAAQIECBAgQIAAgVoJeGN9rYajksYcjlo2Rpwv1ZbG+n9EvCTipREpOdgv/yQWUqS7BP+XiDMRqfxixJcifj0iJQ7vj1AIEAiBg8fPFRy8A7DAYYUAAQIECBAgQIAAAQIEaiRQvhOsRk3TlCEFLsd55eRfqupiRLoTMJWbH/944t8f6S29PT77yb+06VDEuyLSs43p8WGFAIGewIFj6en6K0UC8IqFJQIECBAgQIAAAQIECBCol4AEYHE8tsXqD0e8LeKjESkBlh6lTbEnYiFlVxz8zoi7Ik5FpDvzvhjxpoi1EUtd0lj/YO+if1+6+I7e+jdL29Nqf9v3z7DPJgKdFJiamsoePlFOAHoHYCd/GHSaAAECBAgQIECAAAECDRDwCHBxkA4WV4dee3mc+cGI9D69fklJv2f24rXx+bKIeyNGVVZGxekR3vR+vy0RL464NeIDER+PyJeU6EzlyRFfe2zpyj9pWypPffzDvwQIHD51PrtwaaoAsWPj6sK6FQIECBAgQIAAAQIECBAgUBcBdwDOPhJ7Y9fHZt89656nx54PRaTk38mIt0Q8JyIl4N4TkUpKpqV38K1PKyMqKQH4yxG/FPEvI9KEHr8e8bqIckl3O6byCxH5LEZKHP582hFl0+Mf/iVAoPz+v4lIs29d5w5APxkECBAgQIAAAQIECBAgUE8BdwAWx+VtsZoe002R7ga8MaL/CGwsDlTeHUetiUjv3HtJxOcj+uUTsfD1iHdEpCTgGyPeGlEu6dHhhWQT0jVTvfmSko+Rlnhsxt9r4zPdlfj2iGdH/FDE8Yh++b1Y2B3xoog7I/57xIqIV0T074pM7xZUCBAIgfIMwCn5t2K5/5/ih4MAAQIECBAgQIAAAQIE6ikgAVgcl3TH3GLKs+Lk5/cqeF985pN//XpTcu+OiKdFvD7i1yIuROTLz8bKVfkN8yz/cewvJwD7p6TE3b6I34pIj/r+YcRbIt4c0S8pWfnSiHQH4KsiXhdxLCJNGpLuGrwn4uEIhQCBECgnALdvWEi+HiEBAgQIECBAgAABAgQIEFhaAbesVOud7pjrlw/0F0qfKSH3u71t6bHadNdduayLDenuvUHjr8sVzLL+sd72F86w/1xs+5WIWyJSNiNNiJISkTsjUvnS4x/+JUDgwPHiBCA7NuSfnOdDgAABAgQIECBAgAABAgTqJSABWO14PK9XXZr198tzVP2p3L7n5pZHvZgeBU7lwuMfA/37k72j/mCgox1EoAMC5XcAbpMA7MCo6yIBAgQIECBAgAABAgSaKyABWO3Ypcd6U0mz+6bHamcrd+V29M/JbVrU4rfH2WtnqCFt+43e9r+YYX9+xuL+7h+LhX8Wkd6J+Cf9jT4JdF2g/AiwOwC7/hOh/wQIECBAgAABAgQIEKi3gHcAVjc+6RnArb3q0jv35ipHYme6SzC95+/6uQ4cYt+PxzlviPhsxP0RabKP9BhvesdfmtX3MxHviiiX/y82fCviaxHp+cb0PsMXRtwX8b9GXIpYSLlunoN3zLPfbgK1FZAArO3QaBgBAgQIECBAgAABAgQIzCAgATgDypCb1ufOSzPwzlf6CcD0vr8qy59FZelR3+dEpBl/U/1pQo+vRKTHeN8fMdPdiR+K7T8a8Q8j0gzAafbjX4349xEpibjQkpKJCoFWCpQTgNtMAtLKcdYpAgQIECBAgAABAgQItEVAArC6kczPAnB+gGrTpBuprHn8o7J/02Qdw0zY8dY4L4VCgMAcAhcuXc4OnSz+J75jY/4//zlOtosAAQIECBAgQIAAAQIECIxBQAKwOvT02Gy/rOwvzPGZZtpN5czjH637d75Hm9MjwOndggqBRgk8cqKfu7/S7O3rJQCvaFgiQIAAAQIECBAgQIAAgboJSABWNyInclUN8lhvev9fKoM8Lvz4kc36d1+zmqu1BAYTOHA8n+vPspWTy7JNa9NT8woBAgQIECBAgAABAgQIEKingFmAqxuXlBV4tFfdfBNgbI7j+glA78qrbgzURGDkAg+XEoDb4/1/ExMTI7+uCxAgQIAAAQIECBAgQIAAgWEFJACHlZv5vK/2Nt8cn3PdXXlr7vQ0665CgEBDBA4cK94BuGODx38bMnSaSYAAAQIECBAgQIAAgc4KSABWO/Sf7VWX7u57xhxVvyC373O5ZYsECNRc4GDpHYDbJABrPmKaR4AAAQIECBAgQIAAAQISgNX+DHw4V90dueX8YjJ/dW/D0fj8ZH6nZQIE6i1w0B2A9R4grSNAgAABAgQIECBAgACBaQISgNNIFrXhC3H2Z3o1/HR8PnuG2t4Y257W2/7u+LwwwzE2ESBQU4GDJ4qPAKd3ACoECBAgQIAAAQIECBAgQKDOAnO9p67O7R5V254XFaf39/XL1v5CfKbtu3PraXFPaT2tvj4iPda7JuJjEW+PSHf5pfVXRrwuIpV7It752JJ/CBBojED5HYDbPQLcmLHTUAIECBAgQIAAAQIECHRVQAKwOPKvjdXXFDc9sfbcWEqRL3vyK73lv43Pn4j4YMSGiJQALJeU/HtZxInyDusECNRb4OHj5woNlAAscFghQIAAAQIECBAgQIAAgRoKeAR4NIPykaj29oh3RaRk3+mI9L6/L0W8OeLpEfdGKAQINEjg1LmL2YmIfDELcF7DMgECBAgQIECAAAECBAjUUcAdgMVR2R2rKaooD0Qlb+hFFfWpgwCBMQscPF58/19qjjsAxzwoLk+AAAECBAgQIECAAAEC8wq4A3BeIgcQIEDgcYEDpQTghtWT2ZqVy/EQIECAAAECBAgQIECAAIFaC0gA1np4NI4AgToJeP9fnUZDWwgQIECAAAECBAgQIEBgUAEJwEGlHEeAQOcFyncA7ti4uvMmAAgQIECAAAECBAgQIECg/gISgPUfIy0kQKAmAuV3AG5bLwFYk6HRDAIECBAgQIAAAQIECBCYQ0ACcA4cuwgQIJAXKCcAd2xcld9tmQABAgQIECBAgAABAgQI1FJAArCWw6JRBAjUUeDg8XOFZpkBuMBhhQABAgQIECBAgAABAgRqKiABWNOB0SwCBOoncODY2UKjJAALHFYIECBAgAABAgQIECBAoKYCEoA1HRjNIkCgXgJTU1PZwyckAOs1KlpDgAABAgQIECBAgAABAoMISAAOouQYAgQ6L3D41PnswqWpgsOODSYBKYBYIUCAAAECBAgQIECAAIFaCkgA1nJYNIoAgboJlN//t2wiy7auW1m3ZmoPAQIECBAgQIAAAQIECBCYJiABOI3EBgIECEwXKM8AvHXdqmxyuV+h06VsIUCAAAECBAgQIECAAIG6CfjrtW4joj0ECNRSoJwANAFILYdJowgQIECAAAECBAgQIEBgBgEJwBlQbCJAgEBZ4MBxE4CUTawTIECAAAECBAgQIECAQDMEJACbMU5aSYDAmAXK7wDcvmHVmFvk8gQIECBAgAABAgQIECBAYDABCcDBnBxFgEDHBcqPAJsBuOM/ELpPgAABAgQIECBAgACBBglIADZosDSVAIHxCZQTgN4BOL6xcGUCBAgQIECAAAECBAgQWJiABODCvBxNgEBHBaYlADeu7qiEbhMgQIAAAQIECBAgQIBA0wQkAJs2YtpLgMCSC1y4dDk7dPJ84breAVjgsEKAAAECBAgQIECAAAECNRaQAKzx4GgaAQL1EHj4xLlpDfEOwGkkNhAgQIAAAQIECBAgQIBATQUkAGs6MJpFgEB9BMqP/66cXJZtXLOiPg3UEgIECBAgQIAAAQIECBAgMIeABOAcOHYRIEAgCRw8drYAke7+m5iYKGyzQoAAAQIECBAgQIAAAQIE6iogAVjXkdEuAgRqI1C+A9D7/2ozNBpCgAABAgQIECBAgAABAgMISAAOgOQQAgS6LXDgePEdgNvjDkCFAAECBAgQIECAAAECBAg0RUACsCkjpZ0ECIxN4OHjxUeAJQDHNhQuTIAAAQIECBAgQIAAAQJDCEgADoHmFAIEuiVwoJQANANwt8ZfbwkQIECAAAECBAgQINB0AQnApo+g9hMgMHKB8jsAt21YNfJrugABAgQIECBAgAABAgQIEKhKQAKwKkn1ECDQWoGDpXcAugOwtUOtYwQIECBAgAABAgQIEGilgARgK4dVpwgQqErg5LmLWYp88Q7AvIZlAgQIECBAgAABAgQIEKi7gARg3UdI+wgQGKtA+fHf1BgJwLEOiYsTIECAAAECBAgQIECAwAIFJAAXCOZwAgS6JVBOAG5YPZmtWbm8Wwh6S4AAAQIECBAgQIAAAQKNFpAAbPTwaTwBAqMWKCcAd2xcPepLqp8AAQIECBAgQIAAAQIECFQqIAFYKafKCBBom0B5AhCP/7ZthPWHAAECBAgQIECAAAEC7ReQAGz/GOshAQKLEDhw7GzhbAnAAocVAgQIECBAgAABAgQIEGiAgARgAwZJEwkQGJ/AwyfKCcBV42uMKxMgQIAAAQIECBAgQIAAgSEEJACHQHMKAQLdESjfAbhjg3cAdmf09ZQAAQIECBAgQIAAAQLtEJAAbMc46gUBAiMSKL8DcJsE4IikVUuAAAECBAgQIECAAAECoxKQAByVrHoJEGi8wOXLU1n5EWB3ADZ+WHWAAAECBAgQIECAAAECnROQAOzckOswAQKDChw5fT67cGmqcLhJQAocVggQIECAAAECBAgQIECgAQISgA0YJE0kQGA8AgeOFycAWTaRZVvXrRxPY1yVAAECBAgQIECAAAECBAgMKSABOCSc0wgQaL/Aw8fPFTq5dd2qbHK5X5sFFCsECBAgQIAAAQIECBAgUHsBf8nWfog0kACBcQmU7wDcsdEMwOMaC9clQIAAAQIECBAgQIAAgeEFJACHt3MmAQItFzhYegR423oJwJYPue4RIECAAAECBAgQIECglQISgK0cVp0iQKAKgXICcMfGVVVUqw4CBAgQIECAAAECBAgQILCkAhKAS8rtYgQINEngYOkdgNvdAdik4dNWAgQIECBAgAABAgQIEOgJSAD6USBAgMAsAgeOFWcB3u4dgLNI2UyAAAECBAgQIECAAAECdRaQAKzz6GgbAQJjFXj4RCkBuME7AMc6IC5OgAABAgQIECBAgAABAkMJSAAOxeYkAgTaLnD+4uXs0MnzhW7ukAAseFghQIAAAQIECBAgQIAAgWYISAA2Y5y0kgCBJRY4dPLctCtu32ASkGkoNhAgQIAAAQIECBAgQIBA7QUkAGs/RBpIgMA4BNIdgOWyduVkeZN1AgQIECBAgAABAgQIECBQewEJwNoPkQYSIECAAAECBAgQIECAAAECBAgQGF5AAnB4O2cSIECAAAECBAgQIECAAAECBAgQqL2ABGDth0gDCRAgQIAAAQIECBAgQIAAAQIECAwvIAE4vJ0zCRAgQIAAAQIECBAgQIAAAQIECNReQAKw9kOkgQQIjENgahwXdU0CBAgQIECAAAECBAgQIDACAQnAEaCqkgABAgQIECBAgAABAgQIECBAgEBdBCQA6zIS2kGAQO0FJiZq30QNJECAAAECBAgQIECAAAEC0wQkAKeR2ECAAAECBAgQIECAAAECBAgQIECgPQISgO0ZSz0hQIAAAQIECBAgQIAAAQIECBAgME1AAnAaiQ0ECBDIsqkp04D4OSBAgAABAgQIECBAgACBdghIALZjHPWCAAECBAgQIECAAAECBAgQIECAwIwCEoAzsthIgAABAgQIECBAgAABAgQIECBAoB0CEoDtGEe9IECAAAECBAgQIECAAAECBAgQIDCjgATgjCw2EiBAYLrAxPRNthAgQIAAAQIECBAgQIAAgdoLSADWfog0kAABAgQIECBAgAABAgQIECBAgMDwAhKAw9s5kwCBFguYA7jFg6trBAgQIECAAAECBAgQ6JiABGDHBlx3CRAgQIAAAQIECBAgQIAAAQIEuiUgAdit8dZbAgQIECBAgAABAgQIECBAgACBjglIAHZswHWXAAECBAgQIECAAAECBAgQIECgWwISgN0ab70lQGARAhMT5gFeBJ9TCRAgQIAAAQIECBAgQGBMAhKAY4J3WQIE6i0wZRaQeg+Q1hEgQIAAAQIECBAgQIDAwAISgANTOZAAAQIECBAgQIAAAQIECBAgQIBA8wQkAJs3ZlpMgAABAgQIECBAgAABAgQIECBAYGABCcCBqRxIgAABAgQIECBAgAABAgQIECBAoHkCEoDNGzMtJkCAAAECBAgQIECAAAECBAgQIDCwgATgwFQOJECAAAECBAgQIECAAAECBAgQINA8AQnA5o2ZFhMgsCQC06cBnliS67oIAQIECBAgQIAAAQIECBCoVkACsFpPtREgQIAAAQIECBAgQIAAAQIECBColYAEYK2GQ2MIECBAgAABAgQIECBAgAABAgQIVCsgAVitp9oIECBAgAABAgQIECBAgAABAgQI1Epgslat0ZguCSzvd/ahhx7qL/okUBuBhw6dzC4eP1Roz759+7Jly7wJsIBihQABAgQIECBAgACBWguU/uZ+4m/xWjda4yoX8Jds5aQqHFDgmXHcFwc81mEECBAgQIAAAQIECBAgQIDA4gW+N6r40uKrUUPTBDwC3LQR014CBAgQIECAAAECBAgQIECAAAECCxBwB+ACsBxaqcCqqO22iA9FXIx4cUTV5eO9Cquqe0fU179rMf1fkwNVN1h9jReo+meu8SAL6EAX7JrUx7q1dVztWarrjvI6Vdftu3ABv9g6emjVP3NdYmy7XdP6V7f2jqs9S3XdUV6n6rqb+F2YHvu9pvcL9c74PNelX676+riAdwD6SRiXQPqFk247PttrwL4RNCQlFlMZRd0p+TeKeh9rsH8aKzDKn7nGogzY8C7YNamPdWvruNqzVNcd5XVGWbfvwgF/wXXssFH+zLWdsu12Tetf3do7rvYs1XVHeZ1R1t2k78IH2v5LVP/mFvAI8Nw+9hIgQIAAAQIECBAgQIAAAQIECBBotIAEYKOHT+MJECBAgAABAgQIECBAgAABAgQIzC0gATi3j70ECBAgQIAAAQIECBAgQIAAAQIECBAg0BGB66KfU71IywoBAgQIEOiagO/Cro24/hIgQIBAWcB3YVnEeiME3AHYiGHSSAIECBAgQIAAAQIECBAgQIAAAQLDCUgADufmLAIECBAgQIAAAQIECBAgQIAAAQKNEJAAbMQwaSQBAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgTqI/BT0ZT/HPGliHMRUxG7IxQCBAgQINAmge+NzvxFxNGIUxH/M+LHIxQCBAgQINAFAX/3dWGU9ZEAAQJzCNwf+1LS75GI/vLuWFYIECBAgEBbBF4UHTkfcTzi/454Z8T9Een7740RCgECBAgQaLvA/dFBf/e1fZT1jwABAnMI/EDs29Xb/wvxmb4UdvfWfRAgQIAAgaYLTEYH7o04G/Hduc5sjOW7I9Ld7/3vwdxuiwQIECBAoFUC/u5r1XA2pzPLmtNULSXQeoG/ih4+0Ppe6iABAgQIdFXg+6PjN0X8XsTf5RCOxfLbI1ZGvCa33SIBAgQIEGijgL/72jiqDeiTBGADBkkT5xXYFkf8cMTbIj4acSgi3T2XYk/EQkq68yA9jnRXRHov0eGIL0a8KWJthEKAAAECBJooUIfvyhf24D42A+Bf9ra9YIZ9NhEgQIAAgSoE6vBdWEU/1EGAAIHOCvSTfTN97lmAysvj2HQXwkz1pG3p8aSbI5aieAR4KZRdgwABAt0RmO27LW3fswCGxXxX/lFcJ13vGbNc70Rs3zvLPpsJECBAgMBiBerwXVjug7/7yiLWRybgDsCR0ap4TALpD4eZ7iyYrzlPjwM+FLEh4mTEWyKeE/HiiPdEpPLUiD+PWJ9WFAIECBAg0FCBcX1Xpnf9pZL+Z9tMJU0M0j9mpv22ESBAgACBqgTG9V1YVfvVQ2DBAullzAqBpgu8LTqQHtNNcTDixohvRiykvDsOXhNxMeIlEZ+P6JdPxMLXI94RkZKAaZbCt0aUS3p0eFV54xzr6ZqpXoUAAQIECIxaoC7flaPup/oJECBAgMBsAr4LZ5OxnQABAg0VuDHa3b+9e88AfXhW7vjfnuX4ZbH9q73jjsTnihmOS3cO9q87yOcLZ6ijv8mt4H0JnwQIECAwCoEbo9L+d9WeAS5QxXelR4AHgHYIAQIECCyZwI1xpaX+Lix3zt99ZRHrIxNISQ2FQNcFXpED+EBuOb94OVZ+t7dhU3y+KL+zt7wuPicWEH/dO88HAQIECBCou0AV35X9u96/bYbO7oht6Xu0f8wMh9hEgAABAgTGKlDFd+FYO+Di3RaQAOz2+Ov94wLP60GkWX+/PAfKp3L7nptbtkiAAAECBNouUMV3Zf97NL1qo1z+UW9D/5jyfusECBAgQGDcAlV8F467D67fYQEJwA4Pvq4/IfC03tK98Xnxia3TF+7Kbeqfk9tkkQABAgQItFag/723mO/Kj4fOfRGvivjunFSa+OMXI85H9O+2z+22SIAAAQIEaiFQxXdhLTqiEd0UMAlIN8ddr68IrI7Frb3VfVc2z7iU3v2X7hK8KuL6GY9Y3MbXxun9/6t0W6+qtO2FveXPxud7e8s+CBAgQIDAUglU9V2Z/idb+l77y4hPR/xBxImIfxKxK+J/j7g/QiFAgAABAnUTqOq7MPXL3311G92OtGeyI/3UTQKzCazP7UiTeMxX+gnA9J6iqktK/r2mVGl61Dj/uLEEYAnIKgECBAiMXKDK78pPRmvT992vRPxERJpU686IN0d8KEIhQIAAAQJ1FKjyu9DffXUc4Q60SQKwA4Osi3MKpP+T0y/p0aP5yrneAWvmO3CI/bvjnBQKAQIECBCok0DV35VfiM69tE4d1BYCBAgQIDCPQJXfhbvjWikUAksq4B2AS8rtYjUUOJtr08rc8myLq3o7zsx2gO0ECBAgQKBlAr4rWzagukOAAAECCxbwXbhgMifUTUACsG4joj1LLZDePdQvgzzWm97/l8ogjws/fqR/CRAgQIBAswV8VzZ7/LSeAAECBBYv4Ltw8YZqGLOABOCYB8Dlxy6Q/k/Oo71WXDdPazbH/n4C8FvzHGs3AQIECBBoi4DvyraMpH4QIECAwLACvguHlXNebQQkAGszFBoyRoGv9q59c3zO9V7MW3Nt/Fpu2SIBAgQIEGi7gO/Kto+w/hEgQIDAfAK+C+cTsr/WAhKAtR4ejVsigc/2rpPu7nvGHNd8QW7f53LLFgkQIECAQNsFfFe2fYT1jwABAgTmE/BdOJ+Q/bUWkACs9fBo3BIJfDh3nTtyy/nF9N/Kq3sbjsbnJ/M7LRMgQIAAgZYL+K5s+QDrHgECBAjMK+C7cF4iBxAgQGBpBW6My031Yk98DlI+HQelcy5EPHuGE97U25+OeesM+20iQIAAAQJNErgxGuu7skkjpq0ECBAgULXAjVGh78KqVdVXW4GJ2rZMwwgMLvC8OPTm3OFbY/nf99bTo7rvze1Li3tK62n16RHp2DURaYbft0d8srf+yvh8XUQq90Q8MyI/C1TarhAgQIAAgToL+K6s8+hoGwECBAgshYDvwqVQdg0CBAiMUGBP1N3/PzeDfM7WlJfHjmNz1HV37MsnGmerx3YCBAgQIFA3gT3RoEG+I/vHzNZ+35WzydhOgAABAnUX2BMN7H/PDfI5W398F84mYzsBAgRGLLAn6h/kF3j/mLmasyt2/kZESvadijgS8cWI/yNibYRCgAABAgSaKLAnGt3/Hhzkc64++q6cS8c+AgQIEKirwJ5o2CDfgf1j5uqH78K5dOwjQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQOD/b+/OQ3St6jiAY65J2r5oFpq2khlqmlo5QoVSWmALbaZFJmVJgRDRMllpoBGiBEXazagoS0gtMYyiMs1KS6+ZJjpm4cXUtExxafn+7Dlw7tP7zgx578yF93PgO+c855xneT/z3+FZCBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAJcAQlwAACj9JREFUAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQGDVBA7Kmb+S3JDck/wtuSo5OdkxmVTm0/nvITW+TXJ8cnny9yGXpT422SJZqtT+NfeHybrk/uTW5KLknclix1jIeF3LmqTKXsma5MbkvqTGxmXbdHw0uTL5R3J78rPkHclmyVzSfl+1W3lBGq3/Q61zkfp93fx9FplniAABAgQIECBAgAABAgQIECBAgMAGF6hFt28kbUFrUn13xg+dcOb5br8np31Ftz0+zrkZe0QyreyRgYVkvF+/XYuJdZ5JZSGdNXdNckzyQNLvW+2+7JSN65LxnLZ9XsZe0Y3Ppd2Xupaae03fOaVdC6I1d+2Ucd0ECBAgQIAAAQIECBAgQIAAAQIENopA3eV2ftIWvWqR7q3J/smLk/cnNyU1XnfR7Z30ZT4bbd+L0645pyYvT/ZM3pT8Lmlz3p32pLJbOu9Mat5dyYnJa5O6i++VyelJW9C7NO0tk3FZSEftf3XyYHJj8t5k3+SApL9Tr/b/bdKuqwxek9T5qv5eUmN1rjZnLu2+HJ2NNrZfPzBq18Jmm/fB0ZhNAgQIECBAgAABAgQIECBAgAABAhtV4F05ei1O1aO2B08502PTvzapefV4bF/ms9EWt+oYc8m4PC4d65KaV4tuk0otHtZ43Sn3hEkT0lfX98+k5tV1j8tCOtq1XJn2Y8YTuu3jurmf6/r75mndnDruXD+Y9nbJ3UmNfTGZVk7NQM0pnydOm6SfAAECBAgQIECAAAECBAgQIECAwIYWqLv/rk9qceqUJQ5+yDCv5j6zmzvf9X+26x83Txrm/Sv1o0eDLx3G6ti7j8bGm98c5taC4bgspKOOUaljLlbqsd2ad3Oy9ZSJj0z/n5N2zLkJ884cxu9KXfPHZat03JbUMc4ZD9omQIAAAQIECMyKwGLvgZkVA7+TAAECBAgQILAaAs/LSXcdTvztJS7gJ934tMddv9bNGTd/PXTUouMuo8HDhu1rU181Ghtvtut4UQa2GA8O27Wo99MpY9X91OQ5w/jZqe8b2uPq3nTU+GLlS8Pg9qkPnzDx0PQ9fuivxUKFAAECBAgQIDCTAhYAZ/Lf7kcTIECAAAECm4DA3t01XJJ2u9NtUl2PurbylNYY1b8fbfebd3Qb23XtarbreHbak87d950+7Ltl6nq0eFKpx38XK8/vBtvCZNe1XvNX623978bP01XvOKxy1H+r9f62vlvSe8F6IzYIECBAgAABAjMkYAFwhv7ZfioBAgQIECCwSQk86f+8mm2n7HfPlP7qrkd/W9m8NYZ6Q1/HX0fHH2/WOw1b+UtrTKmXGq/dzhj2PSj1zkO7qh2S9l7Fs9Ku9xcqBAgQIECAAIGZFJj26MZMYvjRBAgQIECAAIEVFOgX4upR1YVlnvvWZc5b7rR2HfWBkPoC8XJLvZ9vUlnphbZa3Dspqff9vT35RFLliKT9No//PkTiDwECBAgQIDCrAhYAZ/U/73cTIECAAAECqy1we3cBd6a9ttteyWa7jket0DX0dwgu9VXepcbL6bbku8nrk1oAPCGpx5aPTKpcnFz3UMsfAgQIECBAgMCMCngEeEb/8X42AQIECBAgsOoCV3RXcEDXXulmu45n5MTT3i+4Ia/p6u5ge3XtSc32fsJJY31f+xjILumcS/ZP2odG3P0XDIUAAQIECBCYbQELgLP9//frCRAgQIAAgdUTuDyn/tNw+qNTb7NKl3LucN7NUh+3AtdQv7ndkVd37W095ZzlUePLKRdl0k3DxKNSV6rcnXzroZY/BAgQIECAAIEZFrAAOMP/fD+dAAECBAgQWFWB+jDHicMV1N13ZyXTFsNq2vbJsdXYwOUHOd5lwzGPT/2GJY6/e8brnYUPp3xh2Hmn1J+ZcqCT07/jlLFxd1m2O/0OT/uNw4SzU9cioEKAAAECBAgQIECAAAECBAgQIEBgVQTqrrtzknpnXeX6pBbhDkxemLwsqbsDv57UQla9764v89lo+/b94/ZcN6/a47JrOupdgO1YdVfgW5J9knpM95Dkw8klSc05JRmXhXTU2JpkqbJVJlyVtPOdl/ZhyZ5Dff4w9ouhrnkHJouVWkysD5C0Y1b9ksV2MEaAAAECBAgQIECAAAECBAgQIEBgJQS2zEk+n9RdbP3i1aT2DaMLmu/2GQ2ttznXzav2pPKsdPaLcpPO3/o+NuEAC+mr8TXJcsrTM6kWPNsxx/WFGTu4G9837aXK9zOhHefapSYbJ0CAAAECBAjMioBHgGflP+13EiBAgAABApuqwAO5sPckeySnJbUId1dSd7NV/ZvkjOR1yXOTjVXqvXx11+Gbk+8kf0zuTe5Pbkl+nHwq2Ss5IXm4pY5fv/njydqkzlVfQ740KY+667DeA9hKWSxVvtpN+HLX1iRAgAABAgQIECBAgAABAgQIECBAYBMU+Eiuqe7oq0XSfjFw2qV+epj/YOodpk3ST4AAAQIECBAgQIAAAQIECBAgQIDA6gvU+xHbI8m/XMblbJ45Nye1YHjeMuabQoAAAQIECBAgQIAAAQIECBAgQIDARhTYOcfeYpHjfzJjtZhXWc7Xj4/o5r8qbYUAAQIECBAgQIAAAQIECBAgQIAAgVUUmM+5b0pOSl6d1BeA60MfRyY/Stri39Vpb51MKruls96L+LakfcW43plYdw8qBAgQIECAAAECBAgQIECAAAECBAisosB8zt0W+abV12TOzsm0Mt6vPliy37TJ+gkQIECAAAECBAgQIECAAAECBAgQWDmBp+VUH0guSP6Q1BeAawFvXXJhckyyVbJYaQuAd2TSRcn+i002RoAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQIAAAQIECBAgQGBFBP4Du/SfvucT5QcAAAAASUVORK5CYII=\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "Text(0,0.5,u'flux')"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from astromodels import Powerlaw\n",
    "\n",
    "am_plaw = Powerlaw()\n",
    "\n",
    "plaw_with_abs = am_plaw*phabs\n",
    "\n",
    "\n",
    "fig, ax =plt.subplots()\n",
    "\n",
    "energy_grid = np.linspace(.1,10.,1000)\n",
    "\n",
    "ax.loglog(energy_grid,plaw_with_abs(energy_grid))\n",
    "ax.set_xlabel('energy')\n",
    "ax.set_ylabel('flux')\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## XSPEC Settings\n",
    "\n",
    "Many **XSPEC** models depend on external abundances, cross-sections, and cosmological parameters. We provide an interface to control these directly.\n",
    "\n",
    "Simply import the **XSPEC** settings like so:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-01-04T13:13:49.868057Z",
     "start_time": "2018-01-04T13:13:49.861676Z"
    }
   },
   "outputs": [],
   "source": [
    "from astromodels.xspec.xspec_settings import *"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Calling the functions without arguments simply returns their current settings"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-01-04T13:13:57.429906Z",
     "start_time": "2018-01-04T13:13:57.420699Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'angr'"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "xspec_abund()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-01-04T13:14:08.255597Z",
     "start_time": "2018-01-04T13:14:08.246786Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'bcmc'"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "xspec_xsect()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-01-04T13:14:12.396361Z",
     "start_time": "2018-01-04T13:14:12.386707Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(70.0, 0.0, 0.7300000190734863)"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "xspec_cosmo()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To change the settings for abundance and cross-section, provide strings with the normal **XSPEC** naming conventions."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-01-04T13:16:15.789515Z",
     "start_time": "2018-01-04T13:16:15.778286Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'wilm'"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "xspec_abund('wilm')\n",
    "xspec_abund()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-01-04T13:16:41.944975Z",
     "start_time": "2018-01-04T13:16:41.935421Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'bcmc'"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "xspec_xsect('bcmc')\n",
    "xspec_xsect()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To alter the cosmological parameters, one passes either the parameters that should be changed, or all three:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-01-04T13:18:02.930968Z",
     "start_time": "2018-01-04T13:18:02.919884Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(68.0, 0.0, 0.7300000190734863)"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "xspec_cosmo(H0=68.)\n",
    "xspec_cosmo()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2018-01-04T13:18:20.371607Z",
     "start_time": "2018-01-04T13:18:20.361461Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(68.0, 0.10000000149011612, 70.0)"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "xspec_cosmo(H0=68.,q0=.1,lambda_0=70.)\n",
    "xspec_cosmo()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.14"
  },
  "toc": {
   "colors": {
    "hover_highlight": "#DAA520",
    "navigate_num": "#000000",
    "navigate_text": "#333333",
    "running_highlight": "#FF0000",
    "selected_highlight": "#FFD700",
    "sidebar_border": "#EEEEEE",
    "wrapper_background": "#FFFFFF"
   },
   "moveMenuLeft": true,
   "nav_menu": {
    "height": "38px",
    "width": "254px"
   },
   "navigate_menu": true,
   "number_sections": true,
   "sideBar": true,
   "threshold": 4,
   "toc_cell": false,
   "toc_section_display": "block",
   "toc_window_display": false,
   "widenNotebook": false
  },
  "varInspector": {
   "cols": {
    "lenName": 16,
    "lenType": 16,
    "lenVar": 40
   },
   "kernels_config": {
    "python": {
     "delete_cmd_postfix": "",
     "delete_cmd_prefix": "del ",
     "library": "var_list.py",
     "varRefreshCmd": "print(var_dic_list())"
    },
    "r": {
     "delete_cmd_postfix": ") ",
     "delete_cmd_prefix": "rm(",
     "library": "var_list.r",
     "varRefreshCmd": "cat(var_dic_list()) "
    }
   },
   "types_to_exclude": [
    "module",
    "function",
    "builtin_function_or_method",
    "instance",
    "_Feature"
   ],
   "window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}