NACA profiledata.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "0b117020",
   "metadata": {},
   "source": [
    "### NACA 4-Series generator\n",
    "\n",
    "Philipp Schlatter, 2017-2024"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "id": "704d86d2",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pylab as plt\n",
    "from array import *\n",
    "import numpy as np\n",
    "%matplotlib notebook\n",
    "\n",
    "# NACA4412\n",
    "m=0./100\n",
    "p=0./10\n",
    "t=12./100.\n",
    "\n",
    "# wedge spreading orthogonal to camber line\n",
    "wd = 0.0005\n",
    "wd = 0.\n",
    "# wedge spreading in y direction\n",
    "#wd1 = 0.000\n",
    " \n",
    "\n",
    "# -----------------\n",
    "pp=max(p,0.1)     # make sure that you have some refinement at the leading edge even for p=0\n",
    "NpointsLE=200;      #Number of divisions before x/c=p\n",
    "dtheta=(np.pi/2)/NpointsLE;\n",
    "theta = np.arange(0.0, np.pi/2, dtheta)       #Half-cosine distribution up to x=p in order to refine LE\n",
    "x1 = (1-np.cos(theta))*pp\n",
    "dx = x1[NpointsLE-1]-x1[NpointsLE-2];     #Linear distribution from x=p to TE\n",
    "x2 = np.arange(pp,1,dx)\n",
    "x = np.concatenate((x1,x2))\n",
    "\n",
    "# original (open NACA)\n",
    "yt = 5*t*(0.2969*np.sqrt(x)-0.1260*x-0.3516*x**2+0.2843*x**3-0.1015*x**4)\n",
    "# modified to get zero at the end (closed NACA)\n",
    "yt = 5*t*(0.2969*np.sqrt(x)-0.1260*x-0.3516*x**2+0.2843*x**3-0.1036*x**4)\n",
    "\n",
    "# wedge orthogonal to camber line\n",
    "yt+=wd*x\n",
    "\n",
    "if p>0:\n",
    "  yc1 = m/p**2*(2*p*x1-x1**2)\n",
    "  yc2 = m/(1-p)**2*((1-2*p)+2*p*x2-x2**2)\n",
    "else:\n",
    "  yc1 = 0*x1\n",
    "  yc2 = 0*x2    \n",
    "  \n",
    "yc=np.concatenate((yc1,yc2))\n",
    "\n",
    "if p>0:\n",
    "  dd1 = 2*m/p**2*(p-x1)\n",
    "  dd2 = 2*m/(1-p)**2*(p-x2)\n",
    "  th1=np.arctan(dd1)\n",
    "  th2=np.arctan(dd2)\n",
    "else:\n",
    "  dd1 = 0*x1\n",
    "  dd2 = 0*x2\n",
    "  th1 = 0*x1\n",
    "  th2 = 0*x2     \n",
    "th = np.concatenate((th1,th2))\n",
    "\n",
    "# only move in y direction (as in xfoil)\n",
    "#th = 0*th\n",
    "\n",
    "xu = x - yt*np.sin(th)\n",
    "yu = yc + yt*np.cos(th)\n",
    "xl = x + yt*np.sin(th)\n",
    "yl = yc - yt*np.cos(th)\n",
    "\n",
    "#Multiply times chord before generating profile\n",
    "chord = 1; #chord in [mm]\n",
    "xu = xu*chord;\n",
    "yu = yu*chord;\n",
    "xl = xl*chord;\n",
    "yl = yl*chord;\n",
    "x1 = x1*chord;\n",
    "yc1 = yc1*chord;\n",
    "x2 = x2*chord;\n",
    "yc2 = yc2*chord;\n",
    "\n",
    "\n",
    "# add wedge in y\n",
    "#yu+=wd1*x\n",
    "#yl-=wd1*x\n",
    "\n",
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 51,
   "id": "8c8c7dc7",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "/* global mpl */\n",
       "window.mpl = {};\n",
       "\n",
       "mpl.get_websocket_type = function () {\n",
       "    if (typeof WebSocket !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof MozWebSocket !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert(\n",
       "            '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",
       "\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 = document.createElement('div');\n",
       "    this.root.setAttribute('style', 'display: inline-block');\n",
       "    this._root_extra_style(this.root);\n",
       "\n",
       "    parent_element.appendChild(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 (fig.ratio !== 1) {\n",
       "            fig.send_message('set_device_pixel_ratio', {\n",
       "                device_pixel_ratio: fig.ratio,\n",
       "            });\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 = document.createElement('div');\n",
       "    titlebar.classList =\n",
       "        'ui-dialog-titlebar ui-widget-header ui-corner-all ui-helper-clearfix';\n",
       "    var titletext = document.createElement('div');\n",
       "    titletext.classList = 'ui-dialog-title';\n",
       "    titletext.setAttribute(\n",
       "        'style',\n",
       "        'width: 100%; text-align: center; padding: 3px;'\n",
       "    );\n",
       "    titlebar.appendChild(titletext);\n",
       "    this.root.appendChild(titlebar);\n",
       "    this.header = titletext;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = (this.canvas_div = document.createElement('div'));\n",
       "    canvas_div.setAttribute('tabindex', '0');\n",
       "    canvas_div.setAttribute(\n",
       "        'style',\n",
       "        'border: 1px solid #ddd;' +\n",
       "            'box-sizing: content-box;' +\n",
       "            'clear: both;' +\n",
       "            'min-height: 1px;' +\n",
       "            'min-width: 1px;' +\n",
       "            'outline: 0;' +\n",
       "            'overflow: hidden;' +\n",
       "            'position: relative;' +\n",
       "            'resize: both;' +\n",
       "            'z-index: 2;'\n",
       "    );\n",
       "\n",
       "    function on_keyboard_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.key_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'keydown',\n",
       "        on_keyboard_event_closure('key_press')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'keyup',\n",
       "        on_keyboard_event_closure('key_release')\n",
       "    );\n",
       "\n",
       "    this._canvas_extra_style(canvas_div);\n",
       "    this.root.appendChild(canvas_div);\n",
       "\n",
       "    var canvas = (this.canvas = document.createElement('canvas'));\n",
       "    canvas.classList.add('mpl-canvas');\n",
       "    canvas.setAttribute(\n",
       "        'style',\n",
       "        'box-sizing: content-box;' +\n",
       "            'pointer-events: none;' +\n",
       "            'position: relative;' +\n",
       "            'z-index: 0;'\n",
       "    );\n",
       "\n",
       "    this.context = canvas.getContext('2d');\n",
       "\n",
       "    var backingStore =\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        this.context.webkitBackingStorePixelRatio ||\n",
       "        this.context.mozBackingStorePixelRatio ||\n",
       "        this.context.msBackingStorePixelRatio ||\n",
       "        this.context.oBackingStorePixelRatio ||\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        1;\n",
       "\n",
       "    this.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband_canvas = (this.rubberband_canvas = document.createElement(\n",
       "        'canvas'\n",
       "    ));\n",
       "    rubberband_canvas.setAttribute(\n",
       "        'style',\n",
       "        'box-sizing: content-box;' +\n",
       "            'left: 0;' +\n",
       "            'pointer-events: none;' +\n",
       "            'position: absolute;' +\n",
       "            'top: 0;' +\n",
       "            'z-index: 1;'\n",
       "    );\n",
       "\n",
       "    // Apply a ponyfill if ResizeObserver is not implemented by browser.\n",
       "    if (this.ResizeObserver === undefined) {\n",
       "        if (window.ResizeObserver !== undefined) {\n",
       "            this.ResizeObserver = window.ResizeObserver;\n",
       "        } else {\n",
       "            var obs = _JSXTOOLS_RESIZE_OBSERVER({});\n",
       "            this.ResizeObserver = obs.ResizeObserver;\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.resizeObserverInstance = new this.ResizeObserver(function (entries) {\n",
       "        var nentries = entries.length;\n",
       "        for (var i = 0; i < nentries; i++) {\n",
       "            var entry = entries[i];\n",
       "            var width, height;\n",
       "            if (entry.contentBoxSize) {\n",
       "                if (entry.contentBoxSize instanceof Array) {\n",
       "                    // Chrome 84 implements new version of spec.\n",
       "                    width = entry.contentBoxSize[0].inlineSize;\n",
       "                    height = entry.contentBoxSize[0].blockSize;\n",
       "                } else {\n",
       "                    // Firefox implements old version of spec.\n",
       "                    width = entry.contentBoxSize.inlineSize;\n",
       "                    height = entry.contentBoxSize.blockSize;\n",
       "                }\n",
       "            } else {\n",
       "                // Chrome <84 implements even older version of spec.\n",
       "                width = entry.contentRect.width;\n",
       "                height = entry.contentRect.height;\n",
       "            }\n",
       "\n",
       "            // Keep the size of the canvas and rubber band canvas in sync with\n",
       "            // the canvas container.\n",
       "            if (entry.devicePixelContentBoxSize) {\n",
       "                // Chrome 84 implements new version of spec.\n",
       "                canvas.setAttribute(\n",
       "                    'width',\n",
       "                    entry.devicePixelContentBoxSize[0].inlineSize\n",
       "                );\n",
       "                canvas.setAttribute(\n",
       "                    'height',\n",
       "                    entry.devicePixelContentBoxSize[0].blockSize\n",
       "                );\n",
       "            } else {\n",
       "                canvas.setAttribute('width', width * fig.ratio);\n",
       "                canvas.setAttribute('height', height * fig.ratio);\n",
       "            }\n",
       "            /* This rescales the canvas back to display pixels, so that it\n",
       "             * appears correct on HiDPI screens. */\n",
       "            canvas.style.width = width + 'px';\n",
       "            canvas.style.height = height + 'px';\n",
       "\n",
       "            rubberband_canvas.setAttribute('width', width);\n",
       "            rubberband_canvas.setAttribute('height', height);\n",
       "\n",
       "            // And update the size in Python. We ignore the initial 0/0 size\n",
       "            // that occurs as the element is placed into the DOM, which should\n",
       "            // otherwise not happen due to the minimum size styling.\n",
       "            if (fig.ws.readyState == 1 && width != 0 && height != 0) {\n",
       "                fig.request_resize(width, height);\n",
       "            }\n",
       "        }\n",
       "    });\n",
       "    this.resizeObserverInstance.observe(canvas_div);\n",
       "\n",
       "    function on_mouse_event_closure(name) {\n",
       "        /* User Agent sniffing is bad, but WebKit is busted:\n",
       "         * https://bugs.webkit.org/show_bug.cgi?id=144526\n",
       "         * https://bugs.webkit.org/show_bug.cgi?id=181818\n",
       "         * The worst that happens here is that they get an extra browser\n",
       "         * selection when dragging, if this check fails to catch them.\n",
       "         */\n",
       "        var UA = navigator.userAgent;\n",
       "        var isWebKit = /AppleWebKit/.test(UA) && !/Chrome/.test(UA);\n",
       "        if(isWebKit) {\n",
       "            return function (event) {\n",
       "                /* This prevents the web browser from automatically changing to\n",
       "                 * the text insertion cursor when the button is pressed. We\n",
       "                 * want to control all of the cursor setting manually through\n",
       "                 * the 'cursor' event from matplotlib */\n",
       "                event.preventDefault()\n",
       "                return fig.mouse_event(event, name);\n",
       "            };\n",
       "        } else {\n",
       "            return function (event) {\n",
       "                return fig.mouse_event(event, name);\n",
       "            };\n",
       "        }\n",
       "    }\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'mousedown',\n",
       "        on_mouse_event_closure('button_press')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'mouseup',\n",
       "        on_mouse_event_closure('button_release')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'dblclick',\n",
       "        on_mouse_event_closure('dblclick')\n",
       "    );\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    canvas_div.addEventListener(\n",
       "        'mousemove',\n",
       "        on_mouse_event_closure('motion_notify')\n",
       "    );\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'mouseenter',\n",
       "        on_mouse_event_closure('figure_enter')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'mouseleave',\n",
       "        on_mouse_event_closure('figure_leave')\n",
       "    );\n",
       "\n",
       "    canvas_div.addEventListener('wheel', function (event) {\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        on_mouse_event_closure('scroll')(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.appendChild(canvas);\n",
       "    canvas_div.appendChild(rubberband_canvas);\n",
       "\n",
       "    this.rubberband_context = rubberband_canvas.getContext('2d');\n",
       "    this.rubberband_context.strokeStyle = '#000000';\n",
       "\n",
       "    this._resize_canvas = function (width, height, forward) {\n",
       "        if (forward) {\n",
       "            canvas_div.style.width = width + 'px';\n",
       "            canvas_div.style.height = height + 'px';\n",
       "        }\n",
       "    };\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    canvas_div.addEventListener('contextmenu', function (_e) {\n",
       "        event.preventDefault();\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 toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'mpl-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'mpl-button-group';\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",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'mpl-button-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        var button = (fig.buttons[name] = document.createElement('button'));\n",
       "        button.classList = 'mpl-widget';\n",
       "        button.setAttribute('role', 'button');\n",
       "        button.setAttribute('aria-disabled', 'false');\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "\n",
       "        var icon_img = document.createElement('img');\n",
       "        icon_img.src = '_images/' + image + '.png';\n",
       "        icon_img.srcset = '_images/' + image + '_large.png 2x';\n",
       "        icon_img.alt = tooltip;\n",
       "        button.appendChild(icon_img);\n",
       "\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    var fmt_picker = document.createElement('select');\n",
       "    fmt_picker.classList = 'mpl-widget';\n",
       "    toolbar.appendChild(fmt_picker);\n",
       "    this.format_dropdown = fmt_picker;\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = document.createElement('option');\n",
       "        option.selected = fmt === mpl.default_extension;\n",
       "        option.innerHTML = fmt;\n",
       "        fmt_picker.appendChild(option);\n",
       "    }\n",
       "\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\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",
       "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",
       "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], msg['forward']);\n",
       "        fig.send_message('refresh', {});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function (fig, msg) {\n",
       "    var x0 = msg['x0'] / fig.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / fig.ratio;\n",
       "    var x1 = msg['x1'] / fig.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / fig.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,\n",
       "        0,\n",
       "        fig.canvas.width / fig.ratio,\n",
       "        fig.canvas.height / fig.ratio\n",
       "    );\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",
       "    fig.canvas_div.style.cursor = msg['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.handle_history_buttons = function (fig, msg) {\n",
       "    for (var key in msg) {\n",
       "        if (!(key in fig.buttons)) {\n",
       "            continue;\n",
       "        }\n",
       "        fig.buttons[key].disabled = !msg[key];\n",
       "        fig.buttons[key].setAttribute('aria-disabled', !msg[key]);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_navigate_mode = function (fig, msg) {\n",
       "    if (msg['mode'] === 'PAN') {\n",
       "        fig.buttons['Pan'].classList.add('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    } else if (msg['mode'] === 'ZOOM') {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.add('active');\n",
       "    } else {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    }\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",
       "            var img = evt.data;\n",
       "            if (img.type !== 'image/png') {\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",
       "                img.type = 'image/png';\n",
       "            }\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",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                img\n",
       "            );\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        } else if (\n",
       "            typeof evt.data === 'string' &&\n",
       "            evt.data.slice(0, 21) === 'data:image/png;base64'\n",
       "        ) {\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(\n",
       "                \"No handler for the '\" + msg_type + \"' message type: \",\n",
       "                msg\n",
       "            );\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(\n",
       "                    \"Exception inside the 'handler_\" + msg_type + \"' callback:\",\n",
       "                    e,\n",
       "                    e.stack,\n",
       "                    msg\n",
       "                );\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "};\n",
       "\n",
       "function getModifiers(event) {\n",
       "    var mods = [];\n",
       "    if (event.ctrlKey) {\n",
       "        mods.push('ctrl');\n",
       "    }\n",
       "    if (event.altKey) {\n",
       "        mods.push('alt');\n",
       "    }\n",
       "    if (event.shiftKey) {\n",
       "        mods.push('shift');\n",
       "    }\n",
       "    if (event.metaKey) {\n",
       "        mods.push('meta');\n",
       "    }\n",
       "    return mods;\n",
       "}\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * https://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",
       "        }\n",
       "        return obj;\n",
       "    }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function (event, name) {\n",
       "    if (name === 'button_press') {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    // from https://stackoverflow.com/q/1114465\n",
       "    var boundingRect = this.canvas.getBoundingClientRect();\n",
       "    var x = (event.clientX - boundingRect.left) * this.ratio;\n",
       "    var y = (event.clientY - boundingRect.top) * this.ratio;\n",
       "\n",
       "    this.send_message(name, {\n",
       "        x: x,\n",
       "        y: y,\n",
       "        button: event.button,\n",
       "        step: event.step,\n",
       "        modifiers: getModifiers(event),\n",
       "        guiEvent: simpleKeys(event),\n",
       "    });\n",
       "\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",
       "    // Prevent repeat events\n",
       "    if (name === 'key_press') {\n",
       "        if (event.key === this._key) {\n",
       "            return;\n",
       "        } else {\n",
       "            this._key = event.key;\n",
       "        }\n",
       "    }\n",
       "    if (name === 'key_release') {\n",
       "        this._key = null;\n",
       "    }\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.key !== 'Control') {\n",
       "        value += 'ctrl+';\n",
       "    }\n",
       "    else if (event.altKey && event.key !== 'Alt') {\n",
       "        value += 'alt+';\n",
       "    }\n",
       "    else if (event.shiftKey && event.key !== 'Shift') {\n",
       "        value += 'shift+';\n",
       "    }\n",
       "\n",
       "    value += 'k' + event.key;\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, { key: value, 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",
       "\n",
       "///////////////// REMAINING CONTENT GENERATED BY embed_js.py /////////////////\n",
       "// prettier-ignore\n",
       "var _JSXTOOLS_RESIZE_OBSERVER=function(A){var t,i=new WeakMap,n=new WeakMap,a=new WeakMap,r=new WeakMap,o=new Set;function s(e){if(!(this instanceof s))throw new TypeError(\"Constructor requires 'new' operator\");i.set(this,e)}function h(){throw new TypeError(\"Function is not a constructor\")}function c(e,t,i,n){e=0 in arguments?Number(arguments[0]):0,t=1 in arguments?Number(arguments[1]):0,i=2 in arguments?Number(arguments[2]):0,n=3 in arguments?Number(arguments[3]):0,this.right=(this.x=this.left=e)+(this.width=i),this.bottom=(this.y=this.top=t)+(this.height=n),Object.freeze(this)}function d(){t=requestAnimationFrame(d);var s=new WeakMap,p=new Set;o.forEach((function(t){r.get(t).forEach((function(i){var r=t instanceof window.SVGElement,o=a.get(t),d=r?0:parseFloat(o.paddingTop),f=r?0:parseFloat(o.paddingRight),l=r?0:parseFloat(o.paddingBottom),u=r?0:parseFloat(o.paddingLeft),g=r?0:parseFloat(o.borderTopWidth),m=r?0:parseFloat(o.borderRightWidth),w=r?0:parseFloat(o.borderBottomWidth),b=u+f,F=d+l,v=(r?0:parseFloat(o.borderLeftWidth))+m,W=g+w,y=r?0:t.offsetHeight-W-t.clientHeight,E=r?0:t.offsetWidth-v-t.clientWidth,R=b+v,z=F+W,M=r?t.width:parseFloat(o.width)-R-E,O=r?t.height:parseFloat(o.height)-z-y;if(n.has(t)){var k=n.get(t);if(k[0]===M&&k[1]===O)return}n.set(t,[M,O]);var S=Object.create(h.prototype);S.target=t,S.contentRect=new c(u,d,M,O),s.has(i)||(s.set(i,[]),p.add(i)),s.get(i).push(S)}))})),p.forEach((function(e){i.get(e).call(e,s.get(e),e)}))}return s.prototype.observe=function(i){if(i instanceof window.Element){r.has(i)||(r.set(i,new Set),o.add(i),a.set(i,window.getComputedStyle(i)));var n=r.get(i);n.has(this)||n.add(this),cancelAnimationFrame(t),t=requestAnimationFrame(d)}},s.prototype.unobserve=function(i){if(i instanceof window.Element&&r.has(i)){var n=r.get(i);n.has(this)&&(n.delete(this),n.size||(r.delete(i),o.delete(i))),n.size||r.delete(i),o.size||cancelAnimationFrame(t)}},A.DOMRectReadOnly=c,A.ResizeObserver=s,A.ResizeObserverEntry=h,A}; // eslint-disable-line\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis\", \"fa fa-square-o\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pgf\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\", \"webp\"];\n",
       "\n",
       "mpl.default_extension = \"png\";/* global mpl */\n",
       "\n",
       "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.binaryType = comm.kernel.ws.binaryType;\n",
       "    ws.readyState = comm.kernel.ws.readyState;\n",
       "    function updateReadyState(_event) {\n",
       "        if (comm.kernel.ws) {\n",
       "            ws.readyState = comm.kernel.ws.readyState;\n",
       "        } else {\n",
       "            ws.readyState = 3; // Closed state.\n",
       "        }\n",
       "    }\n",
       "    comm.kernel.ws.addEventListener('open', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('close', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('error', updateReadyState);\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",
       "        var data = msg['content']['data'];\n",
       "        if (data['blob'] !== undefined) {\n",
       "            data = {\n",
       "                data: new Blob(msg['buffers'], { type: data['blob'] }),\n",
       "            };\n",
       "        }\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(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 = document.getElementById(id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm);\n",
       "\n",
       "    function ondownload(figure, _format) {\n",
       "        window.open(figure.canvas.toDataURL());\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy, ondownload, element);\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;\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",
       "    fig.cell_info[0].output_area.element.on(\n",
       "        'cleared',\n",
       "        { fig: fig },\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function (fig, msg) {\n",
       "    var width = fig.canvas.width / fig.ratio;\n",
       "    fig.cell_info[0].output_area.element.off(\n",
       "        'cleared',\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "    fig.resizeObserverInstance.unobserve(fig.canvas_div);\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.innerHTML =\n",
       "        '<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 / this.ratio;\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] =\n",
       "        '<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 () {\n",
       "        fig.push_to_output();\n",
       "    }, 1000);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'btn-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'btn-group';\n",
       "    var button;\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",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'btn-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        button = fig.buttons[name] = document.createElement('button');\n",
       "        button.classList = 'btn btn-default';\n",
       "        button.href = '#';\n",
       "        button.title = name;\n",
       "        button.innerHTML = '<i class=\"fa ' + image + ' fa-lg\"></i>';\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message pull-right';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = document.createElement('div');\n",
       "    buttongrp.classList = 'btn-group inline pull-right';\n",
       "    button = document.createElement('button');\n",
       "    button.classList = 'btn btn-mini btn-primary';\n",
       "    button.href = '#';\n",
       "    button.title = 'Stop Interaction';\n",
       "    button.innerHTML = '<i class=\"fa fa-power-off icon-remove icon-large\"></i>';\n",
       "    button.addEventListener('click', function (_evt) {\n",
       "        fig.handle_close(fig, {});\n",
       "    });\n",
       "    button.addEventListener(\n",
       "        'mouseover',\n",
       "        on_mouseover_closure('Stop Interaction')\n",
       "    );\n",
       "    buttongrp.appendChild(button);\n",
       "    var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n",
       "    titlebar.insertBefore(buttongrp, titlebar.firstChild);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._remove_fig_handler = function (event) {\n",
       "    var fig = event.data.fig;\n",
       "    if (event.target !== this) {\n",
       "        // Ignore bubbled events from children.\n",
       "        return;\n",
       "    }\n",
       "    fig.close_ws(fig, {});\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (el) {\n",
       "    el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (el) {\n",
       "    // this is important to make the div 'focusable\n",
       "    el.setAttribute('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",
       "    } else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (event, _name) {\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which === 13) {\n",
       "        this.canvas_div.blur();\n",
       "        // select the cell after this one\n",
       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
       "        IPython.notebook.select(index + 1);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "};\n",
       "\n",
       "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(\n",
       "        'matplotlib',\n",
       "        mpl.mpl_figure_comm\n",
       "    );\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAAXNSR0IArs4c6QAAIABJREFUeF7tnQe4VNXV/hcqFlRAxYZSLHRsoAIqWD5B0U8Be0hUEEWTqBiNCqbBp4BRI1iiAkYpBjUYNbYYMAoWhGCNoVuwABZQ+KsoxPJ/1v5y7jf3MvfeM7NmvXdmz3ueZ57bzj7rrN967+x39jl7n3rff//998KNBEiABEiABEiABEigbAjUowEsm1ozURIgARIgARIgARIIBGgAKQQSIAESIAESIAESKDMCNIBlVnCmSwIkQAIkQAIkQAI0gNQACZAACZAACZAACZQZARrAMis40yUBEiABEiABEiABGkBqgARIgARIgARIgATKjAANYJkVnOmSAAmQAAmQAAmQAA0gNUACJEACJEACJEACZUaABrDMCs50SYAESIAESIAESIAGkBogARIgARIgARIggTIjQANYZgVnuiRAAiRAAiRAAiRAA0gNkAAJkAAJkAAJkECZEaABLLOCM10SIAESIAESIAESoAGkBkiABEiABEiABEigzAjQAJZZwZkuCZAACZAACZAACdAAUgMkQAIkQAIkQAIkUGYEaADLrOBMlwRIgARIgARIgARoAKkBEiABEiABEiABEigzAjSAZVZwpksCJEACJEACJEACNIDUAAmQAAmQAAmQAAmUGQEawDIrONMlARIgARIgARIgARpAaoAESIAESIAESIAEyowADWCZFZzpkgAJkAAJkAAJkAANIDVAAiRAAiRAAiRAAmVGgAawzArOdEmABEiABEiABEiABpAaIAESIAESIAESIIEyI0ADWGYFZ7okQAIkQAIkQAIkQANIDZAACZAACZAACZBAmRGgASyzgjNdEiABEiABEiABEqABpAZIgARIgARIgARIoMwI0ACWWcGZLgmQAAmQAAmQAAnQAFIDJEACJEACJEACJFBmBGgAy6zgTJcESIAESIAESIAEaACpARIgARIgARIgARIoMwI0gGVWcKZLAiRAAiRAAiRAAjSA1AAJkAAJkAAJkAAJlBkBGsAyKzjTJYG6JjBx4kQZOHCgbLHFFrJ48WJp0aJFpVM64ogjZNWqVfKvf/2r0u///e9/S7NmzeSjjz6SadOmySmnnJI1le+++07++Mc/yqRJk+S1116TtWvXynbbbScHH3ywnH/++XL88cfLJptsUqntG2+8Ifvuu69sttlm8t5778muu+5aI6bvv/9eDj/8cHnuuefkpz/9qdx6663V7r9gwQI54IADZMOGDTJv3jw58MADK/b94IMP5IYbbpBXX31VXn/99XCud999twwYMKDS8f7f//t/csstt8iMGTNk0aJF8sUXX8gee+whP/rRj2TIkCGy5ZZb1nVZGZ8ESKDECNAAlljBeLokUOoEEgOoeaiBmTJlSioD+NBDD8lJJ50U9j322GPlr3/960Yovv76a+nbt69Mnz5dzjjjDOnXr5/ssssu8sknn8iTTz4pkydPlvvvv1/69OlTqa2aqJtvvjn87tprr5Urr7yyRsxq+EaNGiUrV66s0QB+++23cuihh8r7778vK1as2MgAzpw5U0499VTZf//9Zccdd5R77703qwFUM3zkkUfKmWeeKWqQt9lmm2A+9Vz1+GoM69WrV+rS4PmTAAkACdAAAmEzFAmQgEhiANXEqVF75ZVXZL/99qtAU90I4H//938Ho6Mjb3//+9/l3Xffld13370S0p/85Cdy++23h9G/s846ayPcS5cula+++iqM9iXb+vXrZbfddgvH0pHHrbfeOoxMVrctW7ZM9tlnn2Am1ZDWNAKoo3tjx46VK664IozUVR0B1NHKZDTypZdekoMOOiirAfzyyy/D6ei5ZW56/MsvvzyYwcMOO4zyIgESIIHUBGgAU6PijiRAAoUgkBjAp59+Wk4//XTp1KlTGJ1LtmwGUEfPmjdvHkb39DJur1695Oqrr5Zf/vKXFe0+/PDDsM9RRx1V6Xi1nbOOCOpooV5i1RE9HdmryVBpbB2Be/DBB8OoW3UGUM2mGtv77rtPPv3003DZu6oBzDy3mgxgdTk8++yzwRBPnTpVfvCDH9SWKv9OAiRAAhUEaAApBhIgASiBxACqGZo9e3YYGdMRPTVuumUzgGrKfvGLX8jjjz8uvXv3Dve/6cjZW2+9VXHpUy+f9u/fP4wAXnDBBalzUkOnRkrNnxq1Vq1aydlnnx1G4qpud955p/z85z8Xva+vadOm1RpAvUdQ89h5553lT3/6U8WoZ6EN4PDhw2XEiBHh/sHMUc3UyXNHEiCBsiVAA1i2pWfiJFA3BDINoJqWdu3ayfbbby//+Mc/gqGqagDVTLVu3TpcutXLvptuuqkkxifTOP72t7+VoUOHhtG/Y445JlVyerw999xTTjvttHD/XWJAdTRODeG2225bcZzly5dLhw4d5LrrrpPBgweH31c3Aqj3CP7617+WhQsXBhOYmXPmJBDLCOA///lP6dq1a7gfUkcjuZEACZBALgRoAHOhxX1JgATMBKqaoWTkTi+V6iXhqgZQJ0roBIirrrpKRo4cGeKrcdNRQB3xu+eee8Lv8jGAiZHUexF79uwZjqP39ukI4IQJE+Tcc8+tyPeEE04QnY2r55NMuMhmAPXcOnbsGO79GzRoUGhfaAOo9yH26NFDttpqK3nxxReDgeZGAiRAArkQoAHMhRb3JQESMBOoaoZ0hE9HxXQJFB0xUyOWuQyMznxVk6eXT/fee++K+CeeeGL4nY7UNW7cOIzg5XIJWCdg6OifLi+jy8AkkzHWrVsXRhx1ooeaK90eeOCBcI/d888/L23atKk4B11e5rzzzgujgjpBo379+qKTVVavXl1plrLeo6f3Cj7zzDNhSZhGjRptxDHtPYBqMNUkq/nUS9dVJ8KYC8QDkAAJlAUBGsCyKDOTJIHiIZBtNOypp54Kxk8vneoaf4kBVFOoa/Lp5d/qtt///veis39znQSio361XSqeP3++tG/fvuKSc00UdZkanaTSsmXLMEJZ3abmb82aNXkZwMT8qWmeNWvWRmsoFk+VeSYkQALFToAGsNgrxPMjgcgIVHc5VCdj6GQGXexZ1/PTte90QoeaO53xm22ZE11DTxeS1tEz3WpbBkYnjeiSKnrvoV5u1pE9vX+u6oicLtCsI4+XXXZZWKhZL7nqq+qml6bV9OlEFr3s26RJE5kzZ044/8xN70vUS9R33HFHuI8wWy61jQDqAtU641fXFtTL0Dp6yY0ESIAE8iVAA5gvObYjARLIi0B1BlCfhtG5c2fR0S01SWoA9dLw22+/HRZRzva0CzVoN954Y3jihy65krkQtF6y1YWgdRKGjijqGoI6s1fvNVQDpmv/6czjJ554Imseei66gLNO/tBLu9m2mpaBydy/pnsA1YTqpnnqAtR6qVgv8eqWPO3k448/lm7duoVz+cMf/iB77bVXpdPRy8C8FJyXHNmIBMqWAA1g2ZaeiZNA3RCoyQz98Ic/DGvaqQHUr2rqLrnkEhkzZkzWk9UFm9u2bSsXXXRRxZM8dIQs81FwOnFD79VTM6mjejryp2v+6XEffvjhjZ4KkgQaN25cWE7mz3/+c8UTSKqeRCEMYE1P8FAzrFsyEaa6iv3mN78Jl6m5kQAJkEBaAjSAaUlxPxIgARIgARIgARKIhAANYCSFZBokQAIkQAIkQAIkkJYADWBaUtyPBEiABEiABEiABCIhQAMYSSGZBgmQAAmQAAmQAAmkJUADmJYU9yMBEiABEiABEiCBSAjQAEZSSKZBAiRAAiRAAiRAAmkJ0ACmJcX9SIAESIAESIAESCASAjSAkRSSaZAACZAACZAACZBAWgJRGcDbbrtNrr/++vBweF1IduzYsdK9e/daWbzwwgvhEUv6KCd9okDaTR8mr08o2HbbbcOD2bmRAAmQAAmQAAkUPwFdZP3zzz+Xpk2byiabbFL8J+xwhtEYwPvvvz+s8q8m8NBDDxVdxf/OO++UBQsWSPPmzatFpw+b79Spk+y9997y0Ucf5WQA9Xmh+txSbiRAAiRAAiRAAqVHQB/3WK6PUYzGAHbp0iUYOX14fLK1a9cuPKh99OjR1aryjDPOkFatWsmmm24aHguVywigmsfGjRuH54U2bNiwoMr/97//LdOnT5devXpV+xzSggYs04ORM6bw5EzOGAK4KNQ0hrUXZ31EpA7grFmzRho1aoRJpsiiRGEAN2zYIA0aNJBp06aFh78n25AhQ4KhmzVrVlbs+mB4HTF88cUX5ZprrsnZAKqAVDhqBD0MoD6k/rjjjqMBdPyn0TcXcnYE/J9Dk7M/Y41AzhjOZF36nD37bxwdW6QoDKDeh7fbbruJ3st3yCGHVBAZNWqUTJo0SfSB8VW3pUuXymGHHSbPPfectG7dOjxIvbYRwPXr14u+ki35BLFq1SoXAzhjxgzp2bMnDaBN4zW21g6TnB0BZxhAciZnfwK4CHzvwLD24qz9d5MmTVwGcDBk7FGiMoCzZ8+Wbt26VVAZOXKkTJkyRRYtWlSJ1Lfffitdu3aVQYMGyQUXXBD+lsYA6j4jRozYiPrUqVPDCCQ3EiABEiABEiCB4iewbt066d+/Pw1g8Zeq5jPM9RKwXvPfbrvtwn1/yaYzenVWkP5O77076qijNgrKEcBSV8rG5+/16TI+UraMyNnGL21rck5Lyr4fWdsZpjmCF2eOAIpEMQKoItJJIJ07dw739CVb+/btpU+fPhtNAlGzp7ODMzdt9/TTT8sDDzwge+yxh2y99da1atPzHgLey1Mr/oLsQM4FwVjrQci5VkQF2YGcC4Ix1UHIOhUm805enD37b3PSoANEYwCTZWDuuOOOcBl4/PjxMmHCBJk/f760aNFChg0bJsuXL5fJkydnRZvmEnDVhp4C8hI9SFclE4acMaUiZ3LGEMBFoaYxrL04e/bfGDL2KNEYQEWho3jXXXddWAhaF3UeM2aM9OjRI1AaMGCALFu2TGbOnEkDaNdNNEfwenOJBlCBEiHnAoGs5TDkjOGsUcgaw9qLMw1gRJeAMVKsHMVTQF6irwtOxRyTnDHVIWdyxhDARaGmMay9OHv23xgy9ihRjQDaceR2BE8BeYk+twzj35ucMTUmZ3LGEMBFoaYxrL04e/bfGDL2KDSABoaeAvISvSHdKJuSM6as5EzOGAK4KNQ0hrUXZ8/+G0PGHoUG0MDQU0BeojekG2VTcsaUlZzJGUMAF4WaxrD24uzZf2PI2KPQABoYegrIS/SGdKNsSs6YspIzOWMI4KJQ0xjWXpw9+28MGXsUGkADQ08BeYnekG6UTckZU1ZyJmcMAVwUahrD2ouzZ/+NIWOPQgNoYOgpIC/RG9KNsik5Y8pKzuSMIYCLQk1jWHtx9uy/MWTsUWgADQw9BeQlekO6UTYlZ0xZyZmcMQRwUahpDGsvzp79N4aMPQoNoIGhp4C8RG9IN8qm5IwpKzmTM4YALgo1jWHtxdmz/8aQsUehATQw9BSQl+gN6UbZlJwxZSVncsYQwEWhpjGsvTh79t8YMvYoNIAGhp4C8hK9Id0om5IzpqzkTM4YArgo1DSGtRdnz/4bQ8YehQbQwNBTQF6iN6QbZVNyxpSVnMkZQwAXhZrGsPbi7Nl/Y8jYo9AAGhh6CshL9IZ0o2xKzpiykjM5YwjgolDTGNZenD37bwwZexQaQANDTwF5id6QbpRNyRlTVnImZwwBXBRqGsPai7Nn/40hY49CA2hg6CkgL9Eb0o2yKTljykrO5IwhgItCTWNYe3H27L8xZOxRaAANDD0F5CV6Q7pRNiVnTFnJmZwxBHBRqGkMay/Onv03how9Cg2ggaGngLxEb0g3yqbkjCkrOZMzhgAuCjWNYe3F2bP/xpCxR6EBNDD0FJCX6A3pRtmUnDFlJWdyxhDARaGmMay9OHv23xgy9ig0gAaGngLyEr0h3SibkjOmrORMzhgCuCjUNIa1F2fP/htDxh6FBtDA0FNAXqI3pBtlU3LGlJWcyRlDABeFmsaw9uLs2X9jyNij0AAaGHoKyEv0hnSjbErOmLKSMzljCOCiUNMY1l6cPftvDBl7FBpAA0NPAXmJ3pBulE3JGVNWciZnDAFcFGoaw9qLs2f/jSFjj0IDaGDoKSAv0RvSjbIpOWPKSs7kjCGAi0JNY1h7cfbsvzFk7FFoAA0MPQXkJXpDulE2JWdMWcmZnDEEcFGoaQxrL86e/TeGjD0KDaCBoaeAvERvSDfKpuSMKSs5kzOGAC4KNY1h7cXZs//GkLFHoQE0MPQUkJfoDelG2ZScMWUlZ3LGEMBFoaYxrL04e/bfGDL2KDSABoaeAvISvSHdKJuSM6as5EzOGAK4KNQ0hrUXZ8/+G0PGHoUG0MDQU0BeojekG2VTcsaUlZzJGUMAF4WaxrD24uzZf2PI2KPQABoYegrIS/SGdKNsSs6YspIzOWMI4KJQ0xjWXpw9+28MGXsUGkADQ08BeYnekG6UTckZU1ZyJmcMAVwUahrD2ouzZ/+NIWOPQgNoYOgpIC/RG9KNsik5Y8pKzuSMIYCLQk1jWHtx9uy/MWTsUWgADQw9BeQlekO6UTYlZ0xZyZmcMQRwUahpDGsvzp79N4aMPQoNoIGhp4C8RG9IN8qm5IwpKzmTM4YALgo1jWHtxdmz/8aQsUehATQw9BSQl+gN6UbZlJwxZSVncsYQwEWhpjGsvTh79t8YMvYoNIAGhp4C8hK9Id0om5IzpqzkTM4YArgo1DSGtRdnz/4bQ8YehQbQwNBTQF6iN6QbZVNyxpSVnMkZQwAXhZrGsPbi7Nl/Y8jYo0RlAG+77Ta5/vrrZeXKldKhQwcZO3asdO/ePSul559/Xq688kpZtGiRrFu3Tlq0aCHnn3++/OxnP0tN1VNAXqJPnVyZ7EjOmEKTMzljCOCiUNMY1l6cPftvDBl7lGgM4P333y9nnnmmqAk89NBDZdy4cXLnnXfKggULpHnz5huRevXVV4P523fffWXrrbcWNYRqAMeMGSODBw9ORdZTQF6iT5VYGe1EzphikzM5YwjgolDTGNZenD37bwwZe5RoDGCXLl2kU6dOcvvtt1dQadeunfTt21dGjx6ditRJJ50UzOCUKVNS7e8pIC/Rp0qsjHYiZ0yxyZmcMQRwUahpDGsvzp79N4aMPUoUBnDDhg3SoEEDmTZtmvTr16+CypAhQ+S1116TWbNm1UpKRwR79+4t11xzjZx77rlZ91+/fr3oK9lUQM2aNZNVq1ZJw4YNa42Ryw4q+hkzZkjPnj2lfv36uTTlvjkQIOccYBl2JWcDvByaknMOsIy7krURYMrmXpy1/27SpImsXbu24P13ytTqfLcoDOCKFStkt912kxdeeEEOOeSQCqijRo2SSZMmyeLFi6sFvfvuu8snn3wi33zzjQwfPlx+9atfVbuv/n3EiBEb/X3q1KnBgHIjARIgARIgARIofgJ673///v1pAIu/VDWfYWIAZ8+eLd26davYeeTIkeFyrt7rV932zjvvyBdffCFz5syRoUOHyq233io/+MEPsu7OEcBSV8rG5+/16TI+UraMyNnGL21rck5Lyr4fWdsZpjmCF2eOAIpEMQJYiEvAKkS9/KuGsaYRw0zBet5D4HXfQ5p/uHLah5wx1SZncsYQwEWhpjGsvTh79t8YMvYoURhAxaCTQDp37hxmASdb+/btpU+fPqkngVx99dXyhz/8QZYtW5aKrKeAvESfKrEy2omcMcUmZ3LGEMBFoaYxrL04e/bfGDL2KNEYwGQZmDvuuCNcBh4/frxMmDBB5s+fH9b4GzZsmCxfvlwmT54cqP3+978Py8O0bds2/KzLwFxyySVy0UUXhZHANJungLxEnyavctqHnDHVJmdyxhDARaGmMay9OHv23xgy9ijRGEBFoaN/1113XVgIumPHjmFNvx49egRKAwYMCCN7M2fODD/fcsstYa1AvQdws802k7322kvOO++8sBbgJptskoqsp4C8RJ8qsTLaiZwxxSZncsYQwEWhpjGsvTh79t8YMvYoURlAO47cjuApIC/R55Zh/HuTM6bG5EzOGAK4KNQ0hrUXZ8/+G0PGHoUG0MDQU0BeojekG2VTcsaUlZzJGUMAF4WaxrD24uzZf2PI2KPQABoYegrIS/SGdKNsSs6YspIzOWMI4KJQ0xjWXpw9+28MGXsUGkADQ08BeYnekG6UTckZU1ZyJmcMAVwUahrD2ouzZ/+NIWOPQgNoYOgpIC/RG9KNsik5Y8pKzuSMIYCLQk1jWHtx9uy/MWTsUWgADQw9BeQlekO6UTYlZ0xZyZmcMQRwUahpDGsvzp79N4aMPQoNoIGhp4C8RG9IN8qm5IwpKzmTM4YALgo1jWHtxdmz/8aQsUehATQw9BSQl+gN6UbZlJwxZSVncsYQwEWhpjGsvTh79t8YMvYoNIAGhp4C8hK9Id0om5IzpqzkTM4YArgo1DSGtRdnz/4bQ8YehQbQwNBTQF6iN6QbZVNyxpSVnMkZQwAXhZrGsPbi7Nl/Y8jYo9AAGhh6CshL9IZ0o2xKzpiykjM5YwjgolDTGNZenD37bwwZexQaQANDTwF5id6QbpRNyRlTVnImZwwBXBRqGsPai7Nn/40hY49CA2hg6CkgL9Eb0o2yKTljykrO5IwhgItCTWNYe3H27L8xZOxRaAANDD0F5CV6Q7pRNiVnTFnJmZwxBHBRqGkMay/Onv03how9Cg2ggaGngLxEb0g3yqbkjCkrOZMzhgAuCjWNYe3F2bP/xpCxR6EBNDD0FJCX6A3pRtmUnDFlJWdyxhDARaGmMay9OHv23xgy9ig0gAaGngLyEr0h3SibkjOmrORMzhgCuCjUNIa1F2fP/htDxh6FBtDA0FNAXqI3pBtlU3LGlJWcyRlDABeFmsaw9uLs2X9jyNij0AAaGHoKyEv0hnSjbErOmLKSMzljCOCiUNMY1l6cPftvDBl7FBpAA0NPAXmJ3pBulE3JGVNWciZnDAFcFGoaw9qLs2f/jSFjj0IDaGDoKSAv0RvSjbIpOWPKSs7kjCGAi0JNY1h7cfbsvzFk7FFoAA0MPQXkJXpDulE2JWdMWcmZnDEEcFGoaQxrL86e/TeGjD0KDaCBoaeAvERvSDfKpuSMKSs5kzOGAC4KNY1h7cXZs//GkLFHoQE0MPQUkJfoDelG2ZScMWUlZ3LGEMBFoaYxrL04e/bfGDL2KDSABoaeAvISvSHdKJuSM6as5EzOGAK4KNQ0hrUXZ8/+G0PGHoUG0MDQU0BeojekG2VTcsaUlZzJGUMAF4WaxrD24uzZf2PI2KPQABoYegrIS/SGdKNsSs6YspIzOWMI4KJQ0xjWXpw9+28MGXsUGkADQ08BeYnekG6UTckZU1ZyJmcMAVwUahrD2ouzZ/+NIWOPQgNoYOgpIC/RG9KNsik5Y8pKzuSMIYCLQk1jWHtx9uy/MWTsUWgADQw9BeQlekO6UTYlZ0xZyZmcMQRwUahpDGsvzp79N4aMPQoNoIGhp4C8RG9IN8qm5IwpKzmTM4YALgo1jWHtxdmz/8aQsUehATQw9BSQl+gN6UbZlJwxZSVncsYQwEWhpjGsvTh79t8YMvYoNIAGhp4C8hK9Id0om5IzpqzkTM4YArgo1DSGtRdnz/4bQ8YehQbQwNBTQF6iN6QbZVNyxpSVnMkZQwAXhZrGsPbi7Nl/Y8jYo9AAGhh6CshL9IZ0o2xKzpiykjM5YwjgolDTGNZenD37bwwZe5SoDOBtt90m119/vaxcuVI6dOggY8eOle7du2el9OCDD8rtt98ur732mqxfvz7sP3z4cDnmmGNSU/UUkJfoUydXJjuSM6bQ5EzOGAK4KNQ0hrUXZ8/+G0PGHiUaA3j//ffLmWeeKWoCDz30UBk3bpzceeedsmDBAmnevPlGpC655BJp2rSpHHnkkdK4cWO5++675YYbbpC5c+fKAQcckIqsp4C8RJ8qsTLaiZwxxSZncsYQwEWhpjGsvTh79t8YMvYo0RjALl26SKdOncKoXrK1a9dO+vbtK6NHj05FSkcBTz/9dPn1r3+dan9PAXmJPlViZbQTOWOKTc7kjCGAi0JNY1h7cfbsvzFk7FGiMIAbNmyQBg0ayLRp06Rfv34VVIYMGRIu8c6aNatWUt999520bNlSrrjiCrnwwguz7q+XivWVbCqgZs2ayapVq6Rhw4a1xshlBxX9jBkzpGfPnlK/fv1cmnLfHAiQcw6wDLuSswFeDk3JOQdYxl3J2ggwZXMvztp/N2nSRNauXVvw/jtlanW+WxQGcMWKFbLbbrvJCy+8IIccckgF1FGjRsmkSZNk8eLFtYLWewevvfZaWbhwoey0005Z99d7BEeMGLHR36ZOnRoMKDcSIAESIAESIIHiJ7Bu3Trp378/DWDxl6rmM0wM4OzZs6Vbt24VO48cOVKmTJkiixYtqvEA9957r5x77rnyl7/8RY4++uhq9+UIYKkrZePz9/p0GR8pW0bkbOOXtjU5pyVl34+s7QzTHMGLM0cARaIYAbRcAtbJIwMHDgyXj48//vg0eqzYx/MeAq/7HnJKsAx2JmdMkcmZnDEEcFGoaQxrL86e/TeGjD1KFAZQMegkkM6dO4dZwMnWvn176dOnT7WTQHTk75xzzhH9qpNFct08BeQl+lxzjH1/csZUmJzJGUMAF4WaxrD24uzZf2PI2KNEYwCTZWDuuOOOcBl4/PjxMmHCBJk/f760aNFChg0bJsuXL5fJkycHamr6zjrrLLnpppvkpJNOqiC51VZbSaNGjVKR9RSQl+hTJVZGO5EzptjkTM4YArgo1DSGtRdnz/4bQ8YeJRoDqCh09O+6664LC0F37NhRxowZIz169AiUBgwYIMuWLZOZM2eGn4844oiss4PPPvtsmThxYiqyngLyEn2qxMpoJ3LGFJucyRlDABeFmsaw9uLs2X9jyNijRGUA7ThyO4KngLxEn1uG8e9Nzpig4qgRAAAgAElEQVQakzM5YwjgolDTGNZenD37bwwZexQaQANDTwF5id6QbpRNyRlTVnImZwwBXBRqGsPai7Nn/40hY49CA2hg6CkgL9Eb0o2yKTljykrO5IwhgItCTWNYe3H27L8xZOxRaAANDD0F5CV6Q7pRNiVnTFnJmZwxBHBRqGkMay/Onv03how9Cg2ggaGngLxEb0g3yqbkjCkrOZMzhgAuCjWNYe3F2bP/xpCxR6EBNDD0FJCX6A3pRtmUnDFlJWdyxhDARaGmMay9OHv23xgy9ig0gAaGngLyEr0h3SibkjOmrORMzhgCuCjUNIa1F2fP/htDxh6FBtDA0FNAXqI3pBtlU3LGlJWcyRlDABeFmsaw9uLs2X9jyNij0AAaGHoKyEv0hnSjbErOmLKSMzljCOCiUNMY1l6cPftvDBl7FBpAA0NPAXmJ3pBulE3JGVNWciZnDAFcFGoaw9qLs2f/jSFjj0IDaGDoKSAv0RvSjbIpOWPKSs7kjCGAi0JNY1h7cfbsvzFk7FFoAA0MPQXkJXpDulE2JWdMWcmZnDEEcFGoaQxrL86e/TeGjD0KDaCBoaeAvERvSDfKpuSMKSs5kzOGAC4KNY1h7cXZs//GkLFHoQE0MPQUkJfoDelG2ZScMWUlZ3LGEMBFoaYxrL04e/bfGDL2KDSABoaeAvISvSHdKJuSM6as5EzOGAK4KNQ0hrUXZ8/+G0PGHoUG0MDQU0BeojekG2VTcsaUlZzJGUMAF4WaxrD24uzZf2PI2KPQABoYegrIS/SGdKNsSs6YspIzOWMI4KJQ0xjWXpw9+28MGXsUGkADQ08BeYnekG6UTckZU1ZyJmcMAVwUahrD2ouzZ/+NIWOPQgNoYOgpIC/RG9KNsik5Y8pKzuSMIYCLQk1jWHtx9uy/MWTsUWgADQw9BeQlekO6UTYlZ0xZyZmcMQRwUahpDGsvzp79N4aMPQoNoIGhp4C8RG9IN8qm5IwpKzmTM4YALgo1jWHtxdmz/8aQsUehATQw9BSQl+gN6UbZlJwxZSVncsYQwEWhpjGsvTh79t8YMvYoNIAGhp4C8hK9Id0om5IzpqzkTM4YArgo1DSGtRdnz/4bQ8YehQbQwNBTQF6iN6QbZVNyxpSVnMkZQwAXhZrGsPbi7Nl/Y8jYo9AAGhh6CshL9IZ0o2xKzpiykjM5YwjgolDTGNZenD37bwwZexQaQANDTwF5id6QbpRNyRlTVnImZwwBXBRqGsPai7Nn/40hY49CA2hg6CkgL9Eb0o2yKTljykrO5IwhgItCTWNYe3H27L8xZOxRaAANDD0F5CV6Q7pRNiVnTFnJmZwxBHBRqGkMay/Onv03how9Cg2ggaGngLxEb0g3yqbkjCkrOZMzhgAuCjWNYe3F2bP/xpCxR6EBNDD0FJCX6A3pRtmUnDFlJWdyxhDARaGmMay9OHv23xgy9ig0gAaGngLyEr0h3SibkjOmrORMzhgCuCjUNIa1F2fP/htDxh6FBtDA0FNAXqI3pBtlU3LGlJWcyRlDABeFmsaw9uLs2X9jyNij0AAaGHoKyEv0hnSjbErOmLKSMzljCOCiUNMY1l6cPftvDBl7FBpAA0NPAXmJ3pBulE3JGVNWciZnDAFcFGoaw9qLs2f/jSFjj0IDaGDoKSAv0RvSjbIpOWPKSs7kjCGAi0JNY1h7cfbsvzFk7FGiMoC33XabXH/99bJy5Urp0KGDjB07Vrp3756Vku5z2WWXycsvvyxLly6Viy++OOyfy+YpIC/R55JfOexLzpgqkzM5YwjgolDTGNZenD37bwwZe5RoDOD9998vZ555pqgJPPTQQ2XcuHFy5513yoIFC6R58+YbkVq2bJmMGTNGOnfuHL4efvjhNIB2PZXcEbzeXEoOhPMJk7Mz4P8cnpwxnDUKWWNYe3GmARSJxgB26dJFOnXqJLfffnuFKtu1ayd9+/aV0aNH16jUI444Qvbff38aQMz/c1FF8XpzKaoki+BkyBlTBHLGcKYBLH3ONICRGMANGzZIgwYNZNq0adKvX78KZQ4ZMkRee+01mTVrVkEM4Pr160VfyaYCatasmaxatUoaNmxY0P8IfSOfMWOG9OzZU+rXr1/QY/Ng/0eAnDFqIGdyxhDARaGmMay9OGv/3aRJE1m7dm3B+28MGXuUKEYAV6xYIbvttpu88MILcsghh1RQGTVqlEyaNEkWL15cEAM4fPhwGTFixEbHmjp1ajCg3EiABEiABEiABIqfwLp166R///40gMVfqprPMDGAs2fPlm7dulXsPHLkSJkyZYosWrSoIAaQI4ClrpSNz9/r02V8pGwZkbONX9rW5JyWlH0/srYzTHMEL84cAeQl4KA/3gOY5t8wzn14zxSmruRMzhgCuCjUNIa1F2feAxiJAVQZ6iQQndGrs4CTrX379tKnTx9OAsH8n5ZkFK83l5KE4XjS5OwIN+PQ5IzhrFHIGsPaizMNYEQGMFkG5o477giXgcePHy8TJkyQ+fPnS4sWLWTYsGGyfPlymTx5coVqdYKIbueee660adNGLr/8ctl8881FjWOazVNAXqJPk1c57UPOmGqTMzljCOCiUNMY1l6cPftvDBl7lCgmgSQYdPTvuuuuCwtBd+zYMazv16NHj/DnAQMGiK79N3PmzApq9erV24igmkXdL83mKSAv0afJq5z2IWdMtcmZnDEEcFGoaQxrL86e/TeGjD1KVAbQjiO3I3gKyEv0uWUY/97kjKkxOZMzhgAuCjWNYe3F2bP/xpCxR6EBNDD0FJCX6A3pRtmUnDFlJWdyxhDARaGmMay9OHv23xgy9ig0gAaGngLyEr0h3SibkjOmrORMzhgCuCjUNIa1F2fP/htDxh6FBtDA0FNAXqI3pBtlU3LGlJWcyRlDABeFmsaw9uLs2X9jyNij0AAaGHoKyEv0hnSjbErOmLKSMzljCOCiUNMY1l6cPftvDBl7FBpAA0NPAXmJ3pBulE3JGVNWciZnDAFcFGoaw9qLs2f/jSFjj0IDaGDoKSAv0RvSjbIpOWPKSs7kjCGAi0JNY1h7cfbsvzFk7FFoAA0MPQXkJXpDulE2JWdMWcmZnDEEcFGoaQxrL86e/TeGjD0KDaCBoaeAvERvSDfKpuSMKSs5kzOGAC4KNY1h7cXZs//GkLFHoQE0MPQUkJfoDelG2ZScMWUlZ3LGEMBFoaYxrL04e/bfGDL2KDSABoaeAvISvSHdKJuSM6as5EzOGAK4KNQ0hrUXZ8/+G0PGHoUG0MDQU0BeojekG2XTtJy///57+frrr2XdunUVry+//LLSz/r3DRs2iB5Tv1b3fdW/f/fdd6IvjVHT91X/rs+y3mSTTSq9Nt10041+V3Uf/XmzzTaTzTffXOrXrx++Vvd9dX/X/bfaaivZcsstw9fMl7ap+pzttJyjFBkwKXLGwSZrDGsvzp79N4aMPQoNoIGhp4C8RG9ItySbqmn6/PPPZdWqVbJ69WpZu3ataN30q74+++wzee2112SHHXaQL774ouL3uk9Vg1eSAOrgpNVgVjWFahTVIO+6667SoEGDir9vvfXWss0226R6JfvqsasazDpIsyhD8n0DVxayxrD24uzZf2PI2KPQABoYegrIS/SGdIuiqRq6NWvWyMqVKyten3zySTB3avISo5f5vbIs5LbFFlsEE1P1pSanptG0bKNsOlqmo3aZo3nZvq/6u8QAJSOG+vXbb7+tGEHM/H3V73W/b775JueRysxRzfXr18tXX31V6VVIxjUdS3PPZhobNmwo+mrUqFGlV3W/09/rSGhMG983cNUkawxrL86e/TeGjD0KDaCBoaeAvERvSNe9qeb8/vvvy7vvvivvvfeeLF++vMLkrVixInz/4YcfhpGkXDc1a9tvv700bty4kklQI6Fm8YADDpDtttuuwjioOdC/VTV5OvoUm2nIlWW2/dWYZ5pCrVGmQdRR2Oeff17at28fjGfyN72kriOvVV86+prtd4U418xj6KhiYhirGkXVimpCdZPtq7bT0c5i2srxfaOu+JM1hrwXZ8/+G0PGHoUG0MDQU0Beojeka26qI1Fq8JYsWSLvvPOOLFu2LJi95KUmT/dJs2mHrJcT9bXTTjuFS7hNmjQJr8zvk5/VuGXbYuSchh96n0JwVm1UZxjVYOorubSf+TXzkn/y+3w+RFRlpiORahKrM4jJ76uaR9VndXq01qUQnK3nUC7tyRpTaS/Onv03how9Cg2ggaGngLxEb0g3ddOPP/5YFi9eHIyevpYuXRq+vvnmm2GUqKZNL6+2aNEivHbfffcKk5eYPf26yy67FKwDLWXOqQtSBDsWG2e9nF31ftCq5lFvNdB7RD/99NNKX/V3akQtm448Jh9Y0nxV06i3C9S2FRvn2s63lP9O1pjqeXH27L8xZOxRaAANDD0F5CV6Q7obNdVLdPPnz5c33nij4vWvf/1L9J686jbtxPbaa6/wUpPXsmXLCsOnP+toHvKyWilwLmTN6upYsXHWDzLZzGFiFmv6m95/mc+ml5xrM4s62rhw4UI59dRTw8gkJ8vkQzpdm9g0nS5r/F5enD37bzyl/CLSAObHLbTyFJCX6PNNV3N9+eWXZd68eeH1yiuvyNtvv531cNrpqLFr3bp1eLVq1aria/PmzYvqHrpi45xvfYq9HTn/b4X0Xkn9X0omKaX5qqYy7a0RmTrQSUn6gSrztfPOO1f8nPm9Gss0I4zFrjPk+VHTGNpenD37bwwZexQaQANDTwF5iT5NutrZ6Kie3rQ/d+7cYPj0kq52XlU3vRy7zz77VLw6duwYbvTXyROlsNUl51LgU6hzJOf8Ser/o16Ors0s6si73n6h99Lmc4+jXmauzSyqadT/eb2EXe4bNY1RgBdnz/4bQ8YehQbQwNBTQF6iz5auxtIRvWeffTa81Phph1N100u0Bx10UHgdeOCBsu+++4ZLUqW8ITmXMifruZOzlWC69gnnI444IvwPf/TRR8EU6qu679U45jrCuO222wYjmNybm3xf9auaSuQtHekoFWYvarowHGs7ihdnz/67tpyK5e80gIZKeArIS/RJurrMypNPPhleTz31VJhBmbnpEiiHHHJIeCWGT0cHYtu8OcfGK998yDlfcrm1y4ezrgupl5mzGcWqxlGXYtIlfNJuumSSjhomk7cyv2Z+r/voJetS2vJhXUr5Fcu5enH27L+LhV1t50EDWBuhGv7uKaBCi14v3/7jH/+QP//5z/L444/LggULKmWmN4/36NGj4rX//vsX1b16hjLV2LTQnL3Os9SPS86YCnpz1vcRnfyVLMSu63Im63NW/VrTZLBsNPQ9KHNUcbfddpPk1bRp0/C9msZiuVfRmzVGMcUfxYuzZ/9d/FT/9wxpAA2V8hRQIUSvb9azZ8+WBx54IBg/XYMv2fSyTNeuXeXYY48Nr06dOoUnUpTbVgjO5cYsn3zJOR9qubcpJs56LjqCWJtZ1L/rsjxpNp1gtuOOO1YYQzWFiTnM/B4xA7qYWKdhV6r7eHH27L9LhTUNoKFSngKyiF6N3qRJk2TixIny1ltvVWSol3VPOOEE6du3r/Ts2TM83aDcNwvncmeXS/7knAut/PctRc6Zj3dMRhR1Iou+9GlAyUt/1vzSbHo5WY1hNnOYaRQtl51LkXUadsW2jxdnz/672BhWdz40gIZKeQooV9Hrm+jf/vY3uemmm8LXZMaumr5+/frJKaecIr169Sq5+2wM5UnVNFfOqQ7KnTYiQM4YUcTMWSeq6DO/q5rCqj/rbOm0m44U6oLzzZo1y/pSs1idSYyZdVp+iP28OHv23wguhYhBA2ig6CmgtKLXSyeTJ0+WMWPGVLqv7/DDD5dzzjlHTj75ZC7ZUEON03I2yIRNRcLIzRNPPCHHHXdc0dzDFWNhyFnC04aS0cNso4iJYUy7VI5ecs5mEPV+RH3C0Y9+9KOSWfaqFDXvpWnP/rtUONMAGirlKaDaRK8z9+655x4ZPnx4eKaubro0w7nnnis/+clPZO+99zZkVj5Na+NcPiR8MyVnX77J0ck5HWe9QqJPa1EzqLfMfPDBB+Fr1Vcak6j3JerklepGEZNHWpbjPdbpqlHzXl6a9uy/C5E34hg0gAbKngKqSfTPPfec/PjHPw6PYdNNP4ledtllwfzp46K4pSfg9eaS/gzKY09yxtSZnAvHWU2iXnKuySDq39I82k/Nn15OTp5zXvWrPiFpq622KtzJR3QkL0179t+lgp8G0FApTwFlE70+rP7nP/+53HnnneGsdRLH0KFD5cILL+QliDzr6PXmkufpRNuMnDGlJWcMZ42irB977LGwTqpOXsk2gqgjizrKqFdsatt0ndVMY1j1Oenl+uHeS9Oe/XdttS6Wv9MAGirhKaCqotcndZx22mkVs3rPO+88ufbaa8MD37nlT8DrzSX/M4qzJTlj6krOGM6JAUxzX6uaPzWIuvj+u+++m/WlayvWtqkBrG4EUX+vBlIvR8e2eWnas/8ulRrQABoq5SmgTNHrm8zpp58ebm7Wf3S99++www4znDmbJgS83lxIuDIBcsYogpwxnHMxgLWdkV5q1iexVGcO9fd6Kbq2TWcr66VkHTncY489Kl577rln+B6xNmJt55jP37007dl/55NnXbQpqAG8++67RZcdOfXUUyvlMm3aNFm3bp2cffbZdZGjW0xPASWi10+GZ511VnhWp86gnDJlCkf9ClhRrzeXAp5iFIciZ0wZyRnDuZAGMM0Zaz9Q0wiiznZOlv6q7ng6STAxg4lBTH5W09igQYM0pwLfx0vTnv03HFKeAQtqANu0aSN33HGHHHnkkZVOZ9asWTJ48GBZvHhxnqdZnM08BaSi/93vfie/+c1vwir5AwYMkAkTJpTF49mQ1fZ6c0HmUAqxyBlTJXLGcEYbwNqy0j5C7zfU0UJdFeKdd96Rt99+O3zVlz5tpbZNZzJXNYbJzzqTWZ/rXBebl6Y9+++64JRPzIIaQB2CXrRoURiCztxUkO3atcvpIeL5JINu4ykgvSTQtm1b0edpnnTSSaKjqPr4Nm6FJeD15lLYsyz9o5EzpobkjOFcbAawtqy/+uqrCmNY1RyqUdS+rKZNzZ9eXlZDmIwa7rXXXmG5Mf3qOUHFS9Oe/Xdt9SiWvxfUAKpAbr31VjnxxBMr5feXv/xFfvrTn4ZPKJ7bbbfdJtdff334tNOhQwcZO3asdO/evdqQOjJ56aWXhuVU9JFBV1xxhVxwwQWpT9FTQFdddZWMHj06mOk33ngjXFrnVngCXm8uhT/T0j4iOWPqR84YzqVmAGuikqyJmM0Y6u90VLG2ZzU3adKkwgyqKUyMoX7Vv1kmp3hp2rP/xqnQFqmgBlAN1J/+9CfRewF79OgRzkxNlj6RQh9FdsMNN9jOtobW999/v5x55pmiJvDQQw+VcePGheVSFixYED65VN1U2B07dhSdTXv++efLCy+8EBZQvvfee8PTM9JsXgLSfzYdjteFSvV8zjjjjDSnw33yIOD15pLHqUTdhJwx5SVnDOeYDGBtxHQWs95jmFxOToyiPmf+zTfflI8//rjGQzRs2DCMElY1hvo7HXip7cqWl6a9+u/aeBbT3wtqANW4qAnTy5XJ/QI6eUEnMei9gZtvvrlb7l26dJFOnTrJ7bffXhFDLzv37ds3jKRV3a688kp55JFHZOHChRV/0tG/119/XV588cVU5+kloGeemSlHHXWkNGrUWN56a2VJPb9X7yMupZUIvN5cUgmojHYiZ0yxyRnDuZwMYG1EtR/Uy8hqBvWVGEP9qmsj1rTpbWNqBLMZRB24UR/hpWmv/rs2XsX094IawCQxfT6iGild2XyfffYJS5d4bmo8dQaTGs9+/fpVhBoyZIi89tprYRSy6qYjlAcccIDcdNNNFX966KGHwlp7OmO5fv36G7XRZVj0lWwqIH38jz58XD/lFGo755zz5Z577haRA0TklUId1vU4h8lzcr6Mk9NO/U7q6F7hvPLTDygrP/xQdt1ll1o/ieYVgI0CAXLGCIGcMZyp6XScv/rmG3nnyy/lrS++kLc+/1ze/M9X/XnZl1/Kt99/X+2BNq1XT7arX1+23HRTOXLffeWOp58u6HPEtf/Wy9P6gIVC9t/pyBTHXgUxgN26dQsjbXrvn466oTcdntbH7Ohl3EMOOaQi/KhRo2TSpElZZx+3bt06zKzVe+2Sbfbs2eHysR5PH69WddPn7o4YMWKj30+dOrWgU+iHDv2FLFqkj3lrLyL/+7i3Yt/OkkkySQYU+2ny/EiABEiABOqYwGciomuCzBWRV//z/Xsioqsd/t8Qy/+d5EE77ii/mDChoGetAz39+/enAbRSVZOll1OnT58ejFOfPn3CS82U5ebPtOeVGEA1cGpGk23kyJFh3TydmVx1UwM4cOBAGTZsWMWf1EDqAss6iUTvwau6oUYA1ZTq/ZJ77tlK5sx5raCfetIyzXW/TRYtkPpPzxAdOC2lS8B6f8viJUukTevWwoe151r19PuTc3pWlj3J2UIvt7ZkXT0vnVjysY72rV4tb61aJW+uXi1vr14tb65aFX736bp1NcJuvNVW0rxxY2my9dayxWabyb4dOsivJ08uaF/IEUCRgowAJpVUg/TUU08FM/joo4+Ga/fHH398MIPHHHNMQUfJMtWDugRcVbFe9xDoY986d+4c7v3TRwh5TrHP7S0vvr297i+Jj5QtI3K28UvbmpzTkrLvV+6s9XYDfc5x1Xv/kvsAa3u8nQ6yZC4lkzlJJPMRp16cvfpvu7JwRyioAax62nPnzg1mUF96Q+hRRx0VRtx0ZLDQm04CUdOks4CTrX379sF8VjcJRE2qzhJOth//+MfhnsG6ngSi/1j6z6CzrW688Ub52c9+VmhcPN5/CHi9uRBwZQLkjFEEOWM4a5RyYK0mqeryMMnPOvEj8574quT16p/eI5+5JEzyvRq/tEubeXGmASzwCKDeU6dLviRLwGQKQg2gGkEVhC4JU+gtWQZGZxvrZeDx48eHJ2foGn86CUWNp35amTx5cgidLAOjS8DoUjBq+nQWcDEsA6Pnp+spXnTRReHm1H/961+BG7fCE/B6cyn8mZb2EckZUz9yxnCOxQDq1TNd5y9Z4iXz6SH6vT6QoKZNZ+nq4tDZZvHqGrZ6Fcu6eWmaBrDABlDXz3v88ceDWdH76/TZvzo5A7Xp6N91110X7uHTNf7GjBlTYUbVnOoTSWbOnFlxOjo7WEfXkoWgdWmYYlkI+uuvv5Z9991Xli5dKvvtt58899xzos9y5FZYAl5vLoU9y9I/GjljakjOGM6lYgD1apLeRpRp7DJH9PThDLU9Q3iHHXao9AzhzCeBJEu1eFL30jQNYIENoIpg9erVcs8998jEiRPDyNXRRx8dRgV1lnC2pVU8heN9bE8Bqeh1Qe1f/epXYaHNrl27BnOdeW+Ed37lcHyvN5dyYJdLjuScC6389yXn/Nnl2rIYWKt5++ijj8IoXvLKfBawfl/TZVrNWZdry3zEW9Xv63rgwYuzZ/+dq5bqan/XewBfffVVueuuu8ITOfR6/w9/+MPwSLhWrVrVVb4FjespoET0O++8sxx33HHhqSDKTdc61BFBboUh4PXmUpizi+co5IypJTljOKNGAL/55pvwCNVMg5f5/XvvvVerwdPVDfSqXHUmb6eddoKs1pFvZbw07dl/55srup2bAdTLsHq/nRpAvfdOLw/r75555plwmTaGiQ2eAsoUvS6s3bt377Cqut5Tce2118qFF17IZUsK8N/i9eZSgFOL6hDkjCknOWM4F8oAfvXVV9WaOzV62nfqZdyaNn2Umj5STe91T17JZVo1fWr+Svnqm5emPftvnAptkQpqALVQOtFDL13qmoB6D9u5554bRv6SYeT77rtPdLatjmiV+uYpoKqi10vr+ki9J554ImDTGc864eXAAw8sdYx1ev5eby51mlQRBidnTFHIGcM5jQHUy7P6lCj94K4vHa2rOpJX23N0NY4+QlXvtcs0eJnf77777iVt8GqrmJemPfvv2nIqlr8X1ADqY1X008oPfvCDMLN2//333yhPNX76zF69EbXUN08BZRO9stXZzUOHDg2rl+t26qmnytVXXy1t2rQpdZx1cv5eby51kkwRByVnTHHIGcM5LHT88cdh1QgdZdOrW4nRS77qpdva7r/Ts9Xbo6ozdzqTVm8D0lG+ct28NO3Zf5dKrQpqAPWpG2pICjH1uxQAegqoJtHrrK7LL79c/vjHP4YZXPrmcPrpp8ull17KEcEcheP15pLjaUS/OzljSkzOheGs7+2Zhk7NXFWDp48SS7Ppgsc6SlfdKN52221X1PfgpcnRcx8vTXv23548CnnsghrAQp5YKRzLU0BpRP/Pf/4zzBLWy+7J1r1793B/oC6AvcUWW5QCxjo9xzSc6/QEIwlOzphCknPNnPUDs1490Xvrkpcau6oGT9/b02y6Tqveb6fmTk2e3m+X+dJl0PQSLrf8CXhp2rP/zj9bbEsaQANvTwHlInqdba1rHurlCJ01ppsuF6P3XuoSPNkuxRvSjqppLpyjShycDDljgJczZ13UWJ8Lr8Yu+Zrt+7Qjd40bN97I0Km5S4yeXprVSY26SkMpT7LAKDP/KF6a9uy/888W25IG0MDbU0D5iF7f7G6//fawBqN+n2zt2rULT1/R1z777MPLDRk1z4ezQTJl25ScMaWPkXMymSLT1GUzdp988klqyHrZVUfn9JVt5E5/V9ujymJknRogcEcvzp79NxCPKRQNoAGfp4Asov/2229lxowZYTb2ww8/LPrJONl0LUFdlFuXldFnMpf75QkLZ4N0yq4pOWNKXkqcdVKbrm6g9zTrJAp9Jd+r2csczct8D6uJpL6fqanTZVESg1f1Z/2bLn5s3UqJtTXXumzvxc7aF+MAACAASURBVNmz/65LXrnEpgHMhVaVfT0FVCjR6/0ujz32mDzwwAPy17/+tdKsNP2Ee9RRR8mxxx4rvXr1Cvey6AO8y2krFOdyYpZPruScD7Xc2xQDZ32MpT6doqqpS8xd8lX3SW5ZSZPpjjvuWGHqqjN4+tgy1HtYMbBOw63U9/Hi7Nl/lwpzGkBDpTwF5CH6zz//PKwjqK8nn3wyLGOQuembao8ePSpeeuk49uUHPDgbJBVtU3LGlNaLs47WrVmzJozQZRuxyzR3ua7xqsuH6UzZXXfdNbyS7zNH8PR3xTapzYs1RimlE8WLs2f/XSp0aQANlfIUkJfok3T1Df31118PRlBHBufMmSMaM3PTT9P6DOKDDjqo4qWfwmPavDnHxMqSCzlb6KVvmwtnXaNOPwQmLx2Rq+57vb8ul9E6vRSbaeYSU1f1qz6GrFRvQ8mFdfoKcs+qBLw4e/bfpVJFGkBDpTwF5CX66tLVRxLNnTtXnn322fCaPXu26O+qbrpgqRpCfQKJPulFJ5XoJ3XUZRdDubI2RXMu9PmXyvHI2b9SOllCjZre7qGj93p/XU3mLllMPpcz05mxVUfqspm7cljbjprORTn57+vF2bP/zj9bbEsaQANvTwF5iT5tunrT9SuvvCL/+Mc/ZN68eeG1ePHirM21U1AjqK+OHTuGr+3btw9L0RT7Vteci51Poc6PnHMjqWZOlyvRR4nV9FKTl/n3qqP4tUXV5Ut0FC556dIm2b5Pfleqo3W1ccjn79R0PtRyb+PF2bP/zj3LumlBA2jg7ikgL9Eb0g0LqL788svy0ksvha9vvPGGLFmyRHTWcbZNDWDr1q3DS2cfJ9/vvffetS6xYDnPXNoWI+dczr9U9i1nzmrmvvjiC/n000/DM9CrGrqqJi75u06myGdr0KBBmAWrZq46Q5f8Xj+8lerofT5sCtmmnDVdSI61HcuLs2f/XVtOxfJ3GkBDJTwF5CV6Q7pZm+p9RIsWLQpmMPOlK+vXtGkHVPX5l/rcy+R3usI+YisVzggWnjFi4KxaVwOXGLm0X7VNLvfPZdZBJz7oJIk0L71nVw3d008/zcWJPcX8n2PHoGkAJnMIL86e/bc5adABaAANoD0F5CV6Q7o5Nf3yyy/lzTffDCOES5cuDV+T73WEo7ZNOzJdjDW530i/6qhG5s/6vY52WLZS52zJHdm2rjnrpCcdhdP/WR3Jznxl+53+XWe9Zhq+tE+QqI6rXj7Ve+N0IlVVQ6fmLZvJ23rrrXMaoatrzkhN1XUsssZUwIuzZ/+NIWOPQgNoYOgpIC/RG9ItWFMdOVm2bFl4vfvuuxu99O9pt0aNGoV7ljI70KqdafKzftX9M5eTiJlzWoaI/fLhrG3UtKV9VTVymT/r93op1rrp5VI1cfrSWxxy+aqLD3tfbs2Hs5VJubYna0zlvTh79t8YMvYoNIAGhp4C8hK9IV1YU+3w1RjqkwB0Mdnka7KwbPJztlnKaU5SR2LUCOpr2223DZfn9thjj3D5LPm9XoLWhbJ1hDHNq1yeBar3e6o29XKo8q/upfevZf5Na6q3COiobmZbHVWrzuClffpDmpon+2y22WYVNU5qrV+13pk/J9+rycs0erpfMa+NWc7vG7nooBD7knUhKNZ+DC/Onv137VkVxx40gIY6eArIS/SGdIuqqY7mKH81hbr0ReaN9FVvqk9+znWB2lwSVmOx5ZZbhjXNkpeawmzf6++Sv+nXTTfdNIwMqbHQV7bvs/1OGeilzeSl5izz5+q+1/3U9KrG1GTpK9v32X6nx0RvykvNuF4O1a/ZXvq3TPOezcypeUOMwqH5ZMbj+waOPlljWHtx9uy/MWTsUWgADQw9BeQlekO6Jd9UjY8+DSXz0qBebp41a1Z4DJ7et5h5b5iOTtX00v0LcVmxVMGqKVXTq6aqtpdedteFhtu2bRuMXLK/jq7qKGxNxo5Lj6RXCN830rOy7knWVoLp2ntx9uy/02VW93vRABpq4CkgL9Eb0o2yqYWzmj+9nJmYRL3kqcfLdWQtGanLHNHL9n3VEb9kxDDzq44mZvt91d/pfmrKqhulzByhrG4UU3+f9p42C+coheeUFDk7gc1yWLLGsPbi7Nl/Y8jYo9AAGhh6CshL9IZ0o2xKzpiykjM5YwjgolDTGNZenD37bwwZexQaQANDTwF5id6QbpRNyRlTVnImZwwBXBRqGsPai7Nn/40hY49CA2hg6CkgL9Eb0o2yKTljykrO5IwhgItCTWNYe3H27L8xZOxRaAANDD0F5CV6Q7pRNiVnTFnJmZwxBHBRqGkMay/Onv03how9Cg2ggaGngLxEb0g3yqbkjCkrOZMzhgAuCjWNYe3F2bP/xpCxR6EBNDD0FJCX6A3pRtmUnDFlJWdyxhDARaGmMay9OHv23xgy9ig0gAaGngLyEr0h3SibkjOmrORMzhgCuCjUNIa1F2fP/htDxh6FBtDA0FNAXqI3pBtlU3LGlJWcyRlDABeFmsaw9uLs2X9jyNij0AAaGHoKyEv0hnSjbErOmLKSMzljCOCiUNMY1l6cPftvDBl7FBpAA0NPAXmJ3pBulE3JGVNWciZnDAFcFGoaw9qLs2f/jSFjj0IDaGDoKSAv0RvSjbIpOWPKSs7kjCGAi0JNY1h7cfbsvzFk7FFoAA0MPQXkJXpDulE2JWdMWcmZnDEEcFGoaQxrL86e/TeGjD1KFAbws88+k4svvlgeeeSRQOTEE0+UW265RRo3blwtoQcffFDGjRsnL7/8sqxevVpeffVV2X///XMi6ikgL9HnlGAZ7EzOmCKTMzljCOCiUNMY1l6cPftvDBl7lCgMYO/eveWDDz6Q8ePHByKDBw+Wli1byqOPPlotoSlTpsg777wjTZs2lfPOO48G0K6lkjyC15tLScJwPGlydoSbcWhyxnDWKGSNYe3FmQZQpOQN4MKFC6V9+/YyZ84c6dKlS1Ckft+tWzdZtGiRtGnTpkaVLlu2TPbYYw8aQMz/ctFF8XpzKbpE6/iEyBlTAHLGcKYBLH3ONIARGMC77rpLLr30UlmzZk0lRerl3zFjxsjAgQNpAHH/qyUXiR0mpmTkTM4YArgo1DSGtRdnGsAIDOCoUaNk4sSJsmTJkkpqbN26dTB/w4YNK5gBXL9+vegr2VRAzZo1k1WrVknDhg0L+t+gop8xY4b07NlT6tevX9Bj82D/R4CcMWogZ3LGEMBFoaYxrL04a//dpEkTWbt2bcH7bwwZe5SivQQ8fPhwGTFiRI0Zzps3T6ZPny6TJk2SxYsXV9q3VatWMmjQIBk6dGjBDGB15zR16lRp0KCBvRo8AgmQAAmQAAmQgDuBdevWSf/+/WkA3UnnEUBH1fRV06YTPdR8oS4BcwQwj0IWeROvT5dFnjb89MgZg5ycMZw1ClljWHtx5ghgBJeAk0kgc+fOlYMPPjgoUr/v2rUrJ4Fg/j9LOorX/SUlDcXh5MnZAWqWQ5IzhnNiAJ944gk57rjjeJuOI3YvTfMewAgMoOpOl4FZsWJFWNdPN10GpkWLFpWWgWnbtq2MHj1a+vXrF/b59NNP5b333gvtjj/+eLnvvvvCjOFddtklvNJsngLyEn2avMppH3LGVJucyRlDABeFmsaw9uLs2X9jyNijFO09gLmkpmau6kLQt956a6WFoOvVqyd33323DBgwIBxaJ45kmyH8m9/8RvRevzSbp4C8RJ8mr3Lah5wx1SZncsYQwEWhpjGsvTh79t8YMvYoURhAO4b8juApIC/R55dpvK3IGVNbciZnDAFcFGoaw9qLs2f/jSFjj0IDaGDoKSAv0RvSjbIpOWPKSs7kjCGAi0JNY1h7cfbsvzFk7FFoAA0MPQXkJXpDulE2JWdMWcmZnDEEcFGoaQxrL86e/TeGjD0KDaCBoaeAvERvSDfKpuSMKSs5kzOGAC4KNY1h7cXZs//GkLFHoQE0MPQUkJfoDelG2ZScMWUlZ3LGEMBFoaYxrL04e/bfGDL2KDSABoaeAvISvSHdKJuSM6as5EzOGAK4KNQ0hrUXZ8/+G0PGHoUG0MDQU0BeojekG2VTcsaUlZzJGUMAF4WaxrD24uzZf2PI2KPQABoYegrIS/SGdKNsSs6YspIzOWMI4KJQ0xjWXpw9+28MGXsUGkADQ08BeYnekG6UTckZU1ZyJmcMAVwUahrD2ouzZ/+NIWOPQgNoYOgpIC/RG9KNsik5Y8pKzuSMIYCLQk1jWHtx9uy/MWTsUWgADQw9BeQlekO6UTYlZ0xZyZmcMQRwUahpDGsvzp79N4aMPQoNoIGhp4C8RG9IN8qm5IwpKzmTM4YALgo1jWHtxdmz/8aQsUehATQw9BSQl+gN6UbZlJwxZSVncsYQwEWhpjGsvTh79t8YMvYoNIAGhp4C8hK9Id0om5IzpqzkTM4YArgo1DSGtRdnz/4bQ8YehQbQwNBTQF6iN6QbZVNyxpSVnMkZQwAXhZrGsPbi7Nl/Y8jYo9AAGhh6CshL9IZ0o2xKzpiykjM5YwjgolDTGNZenD37bwwZexQaQANDTwF5id6QbpRNyRlTVnImZwwBXBRqGsPai7Nn/40hY49CA2hg6CkgL9Eb0o2yKTljykrO5IwhgItCTWNYe3H27L8xZOxRaAANDD0F5CV6Q7pRNiVnTFnJmZwxBHBRqGkMay/Onv03how9Cg2ggaGngLxEb0g3yqbkjCkrOZMzhgAuCjWNYe3F2bP/xpCxR6EBNDD0FJCX6A3pRtmUnDFlJWdyxhDARaGmMay9OHv23xgy9ig0gAaGngLyEr0h3SibkjOmrORMzhgCuCjUNIa1F2fP/htDxh6FBtDA0FNAXqI3pBtlU3LGlJWcyRlDABeFmsaw9uLs2X9jyNij0AAaGHoKyEv0hnSjbErOmLKSMzljCOCiUNMY1l6cPftvDBl7FBpAA0NPAXmJ3pBulE3JGVNWciZnDAFcFGoaw9qLs2f/jSFjj0IDaGDoKSAv0RvSjbIpOWPKSs7kjCGAi0JNY1h7cfbsvzFk7FFoAA0MPQXkJXpDulE2JWdMWcmZnDEEcFGoaQxrL86e/TeGjD0KDaCBoaeAvERvSDfKpuSMKSs5kzOGAC4KNY1h7cXZs//GkLFHoQE0MPQUkJfoDelG2ZScMWUlZ3LGEMBFoaYxrL04e/bfGDL2KDSABoaeAvISvSHdKJuSM6as5EzOGAK4KNQ0hrUXZ8/+G0PGHoUG0MDQU0BeojekG2VTcsaUlZzJGUMAF4WaxrD24uzZf2PI2KPQABoYegrIS/SGdKNsSs6YspIzOWMI4KJQ0xjWXpw9+28MGXsUGkADQ08BeYnekG6UTckZU1ZyJmcMAVwUahrD2ouzZ/+NIWOPQgNoYOgpIC/RG9KNsik5Y8pKzuSMIYCLQk1jWHtx9uy/MWTsUWgADQw9BeQlekO6UTYlZ0xZyZmcMQRwUahpDGsvzp79N4aMPUoUBvCzzz6Tiy++WB555JFA5MQTT5RbbrlFGjdunJWQCuqXv/ylPPHEE/L2229Lo0aN5Oijj5Zrr71WmjZtmpqqp4C8RJ86uTLZkZwxhSZncsYQwEWhpjGsvTh79t8YMvYoURjA3r17ywcffCDjx48PRAYPHiwtW7aURx99NCuhtWvXyimnnCLnnXee7LfffqIG8pJLLpFvvvlGXnrppdRUPQXkJfrUyZXJjuSMKTQ5kzOGAC4KNY1h7cXZs//GkLFHKXkDuHDhQmnfvr3MmTNHunTpEojo9926dZNFixZJmzZtUlGaN2+eHHzwwfLuu+9K8+bNU7XxFJCX6FMlVkY7kTOm2ORMzhgCuCjUNIa1F2fP/htDxh6l5A3gXXfdJZdeeqmsWbOmEg29/DtmzBgZOHBgKkpPPfWU9OrVKxynYcOGqdp4CshL9KkSK6OdyBlTbHImZwwBXBRqGsPai7Nn/40hY49S8gZw1KhRMnHiRFmyZEklGq1btw7mb9iwYbVS+vrrr+Wwww6Ttm3byj333FPt/uvXrxd9JZsKqFmzZrJq1arUprHWk/nPDir6GTNmSM+ePaV+/fppm3G/HAmQc47A8tydnPMEl2Mzcs4RmGF3sjbAy6GpF2ftv5s0aSJ6S1jaQZ8cTrskdi1aAzh8+HAZMWJEjRD1su306dNl0qRJsnjx4kr7tmrVSgYNGiRDhw6t8RgqrlNPPVXee+89mTlzZo1CqO6cpk6dKg0aNCiJgvMkSYAESIAESKDcCaxbt0769+9PA1iMQtBRNX3VtOlEDzVf+V4CVvN32mmnhZnATz/9tOywww41xuMIYDEqxXZOXp8ubWcVX2tyxtSUnDGcNQpZY1h7ceYIoEjRjgCmlVYyCWTu3LlhEodu+n3Xrl1rnASSmL+lS5fKM888IzvuuGPakBX7ed5D4HXfQ85JRt6AnDEFJmdyxhDARaGmMay9OHv23xgy9iglbwAVgS4Ds2LFChk3blwgosvAtGjRotIyMHp/3+jRo6Vfv35huZeTTz5ZXnnlFXnsscdk5513riC5/fbby+abb56KrKeAvESfKrEy2omcMcUmZ3LGEMBFoaYxrL04e/bfGDL2KFEYwE8//XSjhaBvvfXWSgtB16tXT+6++24ZMGCALFu2TPbYY4+s9HQ08IgjjkhF1lNAXqJPlVgZ7UTOmGKTMzljCOCiUNMY1l6cPftvDBl7lCgMoB1DfkfwFJCX6PPLNN5W5IypLTmTM4YALgo1jWHtxdmz/8aQsUehATQw9BSQl+gN6UbZlJwxZSVncsYQwEWhpjGsvTh79t8YMvYoNIAGhp4C8hK9Id0om5IzpqzkTM4YArgo1DSGtRdnz/4bQ8YehQbQwNBTQF6iN6QbZVNyxpSVnMkZQwAXhZrGsPbi7Nl/Y8jYo9AAGhh6CshL9IZ0o2xKzpiykjM5YwjgolDTGNZenD37bwwZexQaQANDTwF5id6QbpRNyRlTVnImZwwBXBRqGsPai7Nn/40hY49CA2hg6CkgL9Eb0o2yKTljykrO5IwhgItCTWNYe3H27L8xZOxRaAANDD0F5CV6Q7pRNiVnTFnJmZwxBHBRqGkMay/Onv03how9Cg2ggaGngLxEb0g3yqbkjCkrOZMzhgAuCjWNYe3F2bP/xpCxR6EBNDD0FJCX6A3pRtmUnDFlJWdyxhDARaGmMay9OHv23xgy9ig0gAaGngLyEr0h3SibkjOmrORMzhgCuCjUNIa1F2fP/htDxh6FBtDA0FNAXqI3pBtlU3LGlJWcyRlDABeFmsaw9uLs2X9jyNij0AAaGHoKyEv0hnSjbErOmLKSMzljCOCiUNMY1l6cPftvDBl7FBpAA0NPAXmJ3pBulE3JGVNWciZnDAFcFGoaw9qLs2f/jSFjj0IDaGDoKSAv0RvSjbIpOWPKSs7kjCGAi0JNY1h7cfbsvzFk7FFoAA0MPQXkJXpDulE2JWdMWcmZnDEEcFGoaQxrL86e/TeGjD0KDaCBoaeAvERvSDfKpuSMKSs5kzOGAC4KNY1h7cXZs//GkLFHoQE0MPQUkJfoDelG2ZScMWUlZ3LGEMBFoaYxrL04e/bfGDL2KDSABoaeAvISvSHdKJuSM6as5EzOGAK4KNQ0hrUXZ8/+G0PGHoUG0MDQU0BeojekG2VTcsaUlZzJGUMAF4WaxrD24uzZf2PI2KPQABoYegrIS/SGdKNsSs6YspIzOWMI4KJQ0xjWXpw9+28MGXsUGkADQ08BeYnekG6UTckZU1ZyJmcMAVwUahrD2ouzZ/+NIWOPQgNoYOgpIC/RG9KNsik5Y8pKzuSMIYCLQk1jWHtx9uy/MWTsUWgADQw9BeQlekO6UTYlZ0xZyZmcMQRwUahpDGsvzp79N4aMPQoNoIGhp4C8RG9IN8qm5IwpKzmTM4YALgo1jWHtxdmz/8aQsUehATQw9BSQl+gN6UbZlJwxZSVncsYQwEWhpjGsvTh79t8YMvYoNIAGhp4C8hK9Id0om5IzpqzkTM4YArgo1DSGtRdnz/4bQ8YehQbQwNBTQF6iN6QbZVNyxpSVnMkZQwAXhZrGsPbi7Nl/Y8jYo9AAGhh6CshL9IZ0o2xKzpiykjM5YwjgolDTGNZenD37bwwZexQaQANDTwF5id6QbpRNyRlTVnImZwwBXBRqGsPai7Nn/40hY49CA2hg6CkgL9Eb0o2yKTljykrO5IwhgItCTWNYe3H27L8xZOxRaAANDD0F5CV6Q7pRNiVnTFnJmZwxBHBRqGkMay/Onv03how9Cg2ggaGngLxEb0g3yqbkjCkrOZMzhgAuCjWNYe3F2bP/xpCxR6EBNDD0FJCX6A3pRtmUnDFlJWdyxhDARaGmMay9OHv23xgy9ig0gAaGngLyEr0h3SibkjOmrORMzhgCuCjUNIa1F2fP/htDxh4lCgP42WefycUXXyyPPPJIIHLiiSfKLbfcIo0bN66W0PDhw+W+++6T999/XzbffHPp3LmzjBw5Urp06ZKaqqeAvESfOrky2ZGcMYUmZ3LGEMBFoaYxrL04e/bfGDL2KFEYwN69e8sHH3wg48ePD0QGDx4sLVu2lEcffbRaQlOnTpWddtpJ9txzT/nqq69kzJgxMm3aNHnzzTdlxx13TEXWU0Beok+VWBntRM6YYpMzOWMI4KJQ0xjWXpw9+28MGXuUkjeACxculPbt28ucOXMqRu/0+27dusmiRYukTZs2qSglYnjqqafkv/7rv3Jqs3btWmnYsGGqNml38hJ92vjlsh85YypNzuSMIYCLQk1jWHtxpgEUKXkDeNddd8mll14qa9asqaRGvfyro3oDBw6sVaUbNmyQm2++Wa655powAtikSZNa2+gOngLyEn2qxMpoJ3LGFJucyRlDABeFmsaw9uLs2X9jyNijlLwBHDVqlEycOFGWLFlSiUbr1q2D+Rs2bFi1lB577DE544wzZN26dbLrrrvKww8/LAcddFC1+69fv170lWwqoGbNmsmqVatcRgBnzJghPXv2lPr169srzSNkJaBvLuTsLw5y9mesEcgZw5msS5+z9t862ONxBQ9HxxapaA2gTtIYMWJEjdnNmzdPpk+fLpMmTZLFixdX2rdVq1YyaNAgGTp0aLXH+PLLL2XlypXBwE2YMEGefvppmTt3brg3MNtW3Tnp/YQNGjSwVYKtSYAESIAESIAEIAR04Kd///40gBDaOQZRU6avmjad6KHmy3oJOImhpvGcc86pdtSQI4A5FrEEdueICaZI5EzOGAK4KNQ0hrUXZ44ARnAPYDIJREfuDj744KBI/b5r1645TQLRdnvvvbf86Ec/Eh3pS7N53kPgdd9DmrzKaR9yxlSbnMkZQwAXhZrGsPbi7Nl/Y8jYoxTtJeBcUtNlYFasWCHjxo0LzXQZmBYtWlRaBqZt27YyevRo6devn+ilX13zT9cL1Hv/Vq9eLbfddpvcc8898vLLL0uHDh1ShfcUkJfoUyVWRjuRM6bY5EzOGAK4KNQ0hrUXZ8/+G0PGHiUKA/jpp59utBD0rbfeWmkh6Hr16sndd98tAwYMkK+//jpc+9eRQr3MvMMOO4TJH7/85S9rnARSFbengLxEb5dMXEcgZ0w9yZmcMQRwUahpDGsvzp79N4aMPUoUBtCOIb8jeArIS/T5ZRpvK3LG1JacyRlDABeFmsaw9uLs2X9jyNij0AAaGHoKyEv0hnSjbErOmLKSMzljCOCiUNMY1l6cPftvDBl7FBpAA0NPAXmJ3pBulE3JGVNWciZnDAFcFGoaw9qLs2f/jSFjj0IDaGDoKSAv0RvSjbIpOWPKSs7kjCGAi0JNY1h7cfbsvzFk7FFoAA0MPQXkJXpDulE2JWdMWcmZnDEEcFGoaQxrL86e/TeGjD0KDaCBoaeAvERvSDfKpuSMKSs5kzOGAC4KNY1h7cXZs//GkLFHoQE0MPQUkJfoDelG2ZScMWUlZ3LGEMBFoaYxrL04e/bfGDL2KDSABoaeAvISvSHdKJuSM6as5EzOGAK4KNQ0hrUXZ8/+G0PGHoUG0MDQU0BeojekG2VTcsaUlZzJGUMAF4WaxrD24uzZf2PI2KPQABoYegrIS/SGdKNsSs6YspIzOWMI4KJQ0xjWXpw9+28MGXsUGkADQ08BeYnekG6UTckZU1ZyJmcMAVwUahrD2ouzZ/+NIWOPQgNoYOgpIC/RG9KNsik5Y8pKzuSMIYCLQk1jWHtx9uy/MWTsUWgADQw9BeQlekO6UTYlZ0xZyZmcMQRwUahpDGsvzp79N4aMPQoNoIGhp4C8RG9IN8qm5IwpKzmTM4YALgo1jWHtxdmz/8aQsUehATQw9BSQl+gN6UbZlJwxZSVncsYQwEWhpjGsvTh79t8YMvYoNIAGhp4C8hK9Id0om5IzpqzkTM4YArgo1DSGtRdnz/4bQ8YehQbQwNBTQF6iN6QbZVNyxpSVnMkZQwAXhZrGsPbi7Nl/Y8jYo9AAGhh6CshL9IZ0o2xKzpiykjM5YwjgolDTGNZenD37bwwZexQaQANDTwF5id6QbpRNyRlTVnImZwwBXBRqGsPai7Nn/40hY49CA2hg6CkgL9Eb0o2yKTljykrO5IwhgItCTWNYe3H27L8xZOxRaAANDD0F5CV6Q7pRNiVnTFnJmZwxBHBRqGkMay/Onv03how9Cg2ggaGngLxEb0g3yqbkjCkrOZMzhgAuCjWNYe3F2bP/xpCxR6EBNDD0FJCX6A3pRtmUnDFlJWdyxhDARaGmMay9OHv23xgy9ig0gAaGngLyEr0h3SibkjOmrORMzhgCuCjUNIa1F2fP/htDxh6FBtDA0FNAXqI3pBtlU3LGlJWcyRlDABeFmsaw9uLs2X9jyNij0AAaGHoKyEv0hnSjbErOmLKSMzljCOCiUNMY1l6cPftvDBl7FBpAA0NPAXmJ3pBulE3JGVNWciZnDAFcFGoaw9qLs2f/IEdVgwAAEQ9JREFUjSFjj0IDaGDoKSAv0RvSjbIpOWPKSs7kjCGAi0JNY1h7cfbsvzFk7FFoAA0MPQXkJXpDulE2JWdMWcmZnDEEcFGoaQxrL86e/TeGjD0KDaCBoaeAvERvSDfKpuSMKSs5kzOGAC4KNY1h7cXZs//GkLFHoQE0MPQUkJfoDelG2ZScMWUlZ3LGEMBFoaYxrL04e/bfGDL2KDSABoaeAvISvSHdKJuSM6as5EzOGAK4KNQ0hrUXZ8/+G0PGHoUG0MDQU0BeojekG2VTcsaUlZzJGUMAF4WaxrD24uzZf2PI2KPQABoYegrIS/SGdKNsSs6YspIzOWMI4KJQ0xjWXpw9+28MGXsUGkADQ08BeYnekG6UTckZU1ZyJmcMAVwUahrD2ouzZ/+NIWOPEoUB/Oyzz+Tiiy+WRx55JBA58cQT5ZZbbpHGjRunInT++efL+PHjZcyYMXLJJZekaqM7eQrIS/SpkyuTHckZU2hyJmcMAVwUahrD2ouzZ/+NIWOPEoUB7N27t3zwwQfBxOk2ePBgadmypTz66KO1Enr44Ydl+PDh8sknn8jll19OA1grsbh28HpziYuSPRtytjNMcwRyTkOpMPuQdWE41nYUL840gCIlbwAXLlwo7du3lzlz5kiXLl2ClvT7bt26yaJFi6RNmzbV6mv58uWhzd/+9jc5/vjjg/njCGBt/45x/d3rzSUuSvZsyNnOMM0RyDkNpcLsQ9aF4VjbUbw40wBGYADvuusuufTSS2XNmjWVdKSXf/WS7sCBA7Pq67vvvpOjjz5a+vTpI0OGDAkjhjSAtf0rxvd3rzeX+EjZMiJnG7+0rck5LSn7fmRtZ5jmCF6caQAjMICjRo2SiRMnypIlSyppqXXr1sH8DRs2LKvGRo8eLc8880wY/atXr14qA7h+/XrRV7KpgJo1ayarVq2Shg0bptFy6n1U9DNmzJCePXtK/fr1U7fjjrkRIOfceOW7NznnSy63duScGy/L3mRtoZe+rRdn7b+bNGkia9euLXj/nT67ut2zaC8B6315I0aMqJHOvHnzZPr06TJp0iRZvHhxpX1btWolgwYNkqFDh250jJdffjlc8n3llVekadOm4e9pRgCrO6epU6dKgwYN6raSjE4CJEACJEACJJCKwLp166R///40gKlogXfSUTV91bSpaVPzlesl4LFjx4Y2m2yyScXhv/322/CzjugtW7Ysa1iOAIJFAAjn9ekScOolFYKcMeUiZwxnjULWGNZenDkCGMEl4GQSyNy5c+Xggw8OitTvu3btWu0kkNWrV8vKlSsrqfeYY46RM888M1w2rmniSGYjz3sIvO57wPzLlk4UcsbUipzJGUMAF4WaxrD24uzZf2PI2KMU7SXgXFLTZWBWrFgh48aNC810GZgWLVpUWgambdu2ovf99evXL+uh01wCrtpQ7x3QySbvv/9+we8hUNHr5e1evXrxHsBcxJDjvuScI7A8dyfnPMHl2IyccwRm2J2sDfByaOrFObmHXyeQNmrUKIczimfXKAzgp59+utFC0LfeemulhaB1osfdd98tAwYMKJgB1LUH9ZIxNxIgARIgARIggdIjoAM4u+++e+mdeAHOOAoDWAAOeR1Cl5LRkcdtt902zCQu5JZ8OvEYXSzkeZb6scgZU0FyJmcMAVwUahrD2ovz999/L59//nmYCJo5HwCTVXFEoQEsjjpsdBa8PwFTGHImZwwBTBTqGcNZo5A1hjU5+3GmAfRjazoyRW/Cl7oxOadGZdqRnE34Ujcm59SozDuStRlhqgOQcypMee1EA5gXNv9GFL0/Y36KxzAmZ3LGEcBF4ns0hjU5+3GmAfRjazqyrjmos5b1SSZbbLGF6VhsXD0Bcsaog5zJGUMAF4WaxrAmZz/ONIB+bHlkEiABEiABEiABEihKAjSARVkWnhQJkAAJkAAJkAAJ+BGgAfRjyyOTAAmQAAmQAAmQQFESoAEsyrLwpEiABEiABEiABEjAjwANoB9bHpkESIAESIAESIAEipIADWAdleW2226T66+/XlauXCkdOnSQsWPHSvfu3as9m1mzZsmll14q8+fPDyuXX3HFFXLBBRfU0dmXVthcWD/44INy++23y2uvvSY6+0xrM3z4cDnmmGNKK+k6ONtcOGee3gsvvCCHH364dOzYMXDnVjOBXDmrjv/nf/5H7rnnHvnwww/DY69+8YtfyDnnnEPUNRDIlfMf//hHue6662Tp0qXh2bLHHnus3HDDDbLDDjuQczUEnn322dAPvvzyy6EvfOihh6Rv37418mJfWDg50QAWjmXqI91///1y5plnir7BHHrooTJu3Di58847ZcGCBdK8efONjvPOO++EzvG8886T888/X7TD/MlPfiL33nuvnHzyyanjluOOubK+5JJLgsE+8sgjw7Ok9fnR+iY+d+5cOeCAA8oRYaqcc+WcHHTt2rXSqVMn2XvvveWjjz6iAayFdj6c+/TpE9hec801gfPHH38s33zzjRxyyCGpaluOO+XK+fnnnw8fYsaMGSMnnHCCLF++PHxAb9WqVTA13LIT+Otf/xr6M30P0L6sNgPIvrCwSqIBLCzPVEfr0qVLELyONCVbu3btwicfXfuv6nbllVfKI488IgsXLqz4k765vP766/Liiy+milmuO+XKOhsnHQU8/fTT5de//nW5Yqw173w5n3HGGaGT3HTTTeXhhx+mAayFdK6cn3zySVHGb7/9tmy//fa11pE7/C+BXDnrh0R9P3/rrbcqEN5yyy1hRFCf586tdgL16tWr1QCyL6ydYy570ADmQqsA+27YsEEaNGgg06ZNk379+lUccciQIaHz0+HtqluPHj3C6NNNN91U8Sf9pHTaaafJunXrpH79+gU4s/gOkQ/rqhS+++47admyZbjkfuGFF8YHqQAZ5ctZR1d1FFw/xOjoFA1gzcXIh7NeKViyZIkceOCBMmXKFNl6663lxBNPlKuvvlq22mqrAlQ/vkPkw3n27NnhqoG+L/fu3TuMsur7s36wv+OOO+KD5JBRGgPIvrCw4GkAC8uz1qOtWLFCdttttzDsnXkJZtSoUTJp0iRZvHjxRsdo3bq1DBgwQK666qqKv+kbjl4+1uPtuuuutcYtxx3yYV2Vk96fcu2114bR15122qkcMdaacz6c9T6pww47TJ577jlRfet9ljSANaPOh7PehzZz5kw5+uijwwj2qlWrwu0jRx11lNx111211rYcd8iHs3J64IEHZODAgfL111+HS+xqtPV3/ICeTkVpDCD7wnQs0+5FA5iWVIH2S95c1MB169at4qgjR44Mn9AXLVqU1QDqG4s+Fi7Z1EBqB6o3zu6yyy4FOru4DpMP60wCeo/lueeeK3/5y19CB8otO4FcOX/77bfStWtXGTRoUMVEJhrA2tWVK2c9Yq9evYLJ1skfOjFBN53odMopp8iXX37JUcAs2PPhrPdv63vEz372szBhTN+XL7/8cjnooIPkD3/4Q+3F5R6S1gCyLyycWGgAC8cy1ZHyubzAYe9UaDfaKR/WyUH0JnB9o9FL9ccff3x+J1AmrXLlvGbNGtluu+3CfX/Jppfav//++/C76dOnhxEqbpUJ5MpZW5999tnhasObb75ZcTAdzW7fvn24NKz3X3Kzc9ZJfTryp+8XyaYTQ3RlB16lSaewNAaQfWE6lmn3ogFMS6qA++kNxp07dw73PyWbviHrbL3qJoE8+uijYZZwsv34xz8O9wxyEkjNhcmVtR5NR/50iQz9WtuSBAWURUkfKhfOavYytayJ6//C008/HS6Z7bHHHuFeNW4bE8iFs7YeP3686Mx2vSdtm222CQfUEe2TTjpJvvjiC44AViOyXDnrDNbNNttM9INjsul7s97mozOCdWUBbjUTSGMAdRII+8LCKYkGsHAsUx8pWWJAbw7Wy8D6Jj1hwoSwxl+LFi3CpV5905g8eXI4ZjL1XZeA0aVg9I1FZwFzGZjakefKWpmeddZZYcKNdpLJpjfMJ5fQao9afnvkyrkqIV4CTqeZXDmrydOJCHrJfcSIEeEeQL2tQZcs0fccbtkJ5Mp54sSJ4b355ptvrrgErMZ7k002CUtIcctOQPWZjE7rRMcbb7wxTKbRGeu6JBr7Ql/l0AD68q326DrioUsE6L0iusafrh+lw9u66YSPZcuWhZu3k01nB+v9JclC0PpJiAtBpyteLqyPOOKIrDOx9VKavslzq55ALpxpAPNXUq6c9b7iiy66KFwK1kWJdXaqzrrmLOCaa5ArZ132RT/U6wd2XUNUb2P47W9/Gyb9cctOQPs4NXxVt+T9ln2hr3JoAH358ugkQAIkQAIkQAIkUHQEaACLriQ8IRIgARIgARIgARLwJUAD6MuXRycBEiABEiABEiCBoiNAA1h0JeEJkQAJkAAJkAAJkIAvARpAX748OgmQAAmQAAmQAAkUHQEawKIrCU+IBEiABEiABEiABHwJ0AD68uXRSYAESIAESIAESKDoCNAAFl1JeEIkQAIkQAIkQAIk4EuABtCXL49OAiRAAiRAAiRAAkVHgAaw6ErCEyIBEiABEiABEiABXwI0gL58eXQSIIESJqCP/9PHe82ZM6eEs+CpkwAJkMDGBGgAqQoSIAESqIbASSedJAcddFB4KD03EiABEoiJAA1gTNVkLiRAAqkIfPLJJ7LPPvvIxRdfLFdddVVoM3fuXOnevbs89thj0qtXL/n666+lSZMm4fcdOnSQ9evXy69+9Su599575eOPP5bmzZvL0KFDZdCgQalicicSIAESKCYCNIDFVA2eCwmQAIzAE088IX379pXZs2dL27Zt5YADDpDjjz9exo4dG87h8ccflyFDhsibb74Zfj799NPlxRdflJtuukn2228/eeedd2TVqlXh99xIgARIoNQI0ACWWsV4viRAAgUj8NOf/lSeeuqpcJn39ddfl3nz5smWW24Zjj948GDZdttt5Xe/+50sWbJE2rRpIzNmzJCjjz66YPF5IBIgARKoKwI0gHVFnnFJgATqnMBXX30lHTt2lPfff19eeukl2XfffcM5ff/999K0aVO577775PDDD5c//elP0r9/f9H969evX+fnzRMgARIgASsBGkArQbYnARIoWQLz58+XAw88UP7973/LQw89JCeccELIRe/7O+6448K9fptuuqk8+uij0q9fPxrAkq00T5wESKAqARpAaoIESKAsCWzYsEEOPvhg2X///cM9gDfeeKO88cYbsvPOO4eJIR988IFMnjw5sFm2bJnsueeeMn36dF4CLku1MGkSiI8ADWB8NWVGJEACKQhcfvnl8sADD4R7/7bZZhs58sgjwz1/OgtYLwuPGDFCTj755IojDRw4UP7+97/LzTffHCaBvPvuu2GE8LTTTksRjbuQAAmQQHERoAEsrnrwbEiABAAEZs6cKT179pRnnnlGDjvssBDxvffeC/cAjhw5Ui677LIww1eNYbLpsjA6Mqj3Ba5evTosA6M/qzHkRgIkQAKlRoAGsNQqxvMlARJwJaCXgnVmsC4Tw40ESIAEYiVAAxhrZZkXCZBAXgR0xu+uu+4aFoXmRgIkQAKxEqABjLWyzIsESIAESIAESIAEqiFAA0hpkAAJkAAJkAAJkECZEaABLLOCM10SIAESIAESIAESoAGkBkiABEiABEiABEigzAjQAJZZwZkuCZAACZAACZAACdAAUgMkQAIkQAIkQAIkUGYEaADLrOBMlwRIgARIgARIgARoAKkBEiABEiABEiABEigzAjSAZVZwpksCJEACJEACJEACNIDUAAmQAAmQAAmQAAmUGYH/D3AF9nldLh4DAAAAAElFTkSuQmCC\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# plotting----------------------\n",
    "\n",
    "plt.figure(1)\n",
    "plt.plot(x1, yc1, 'b')\n",
    "plt.plot(x2,yc2,'r')\n",
    "plt.plot(xu,yu,'k',xl,yl,'k')\n",
    "plt.axis('equal')\n",
    "plt.title('NACA4412')\n",
    "plt.ylabel('y/c')\n",
    "plt.xlabel('x/c')\n",
    "plt.grid(True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 52,
   "id": "ea55be68",
   "metadata": {},
   "outputs": [],
   "source": [
    "#write out current profile..........\n",
    "f = open('NACA4412_TEGap.txt','w')\n",
    "print(\"NACA4412 KTH Mechanics SimEx interpolated\",file=f)\n",
    "for i in range(xu.size-1,0,-1):\n",
    "    print(\" %18.9f %18.9f\"%(xu[i],yu[i]),file=f)\n",
    "for i in range(0,xl.size):\n",
    "    print(\" %18.9f %18.9f\"%(xl[i],yl[i]),file=f)\n",
    "f.close()\n",
    "\n",
    "f = open('Pressure_side_NACA4412.txt','w')\n",
    "#print(\" NACA4412 KTH Mechanics SimEx interpolated\",file=f)\n",
    "for i in range(0,xu.size):\n",
    "    print(\" %18.9f %18.9f\"%(xu[i],yu[i]),file=f)\n",
    "f.close()\n",
    "\n",
    "f = open('Suction_side_NACA4412.txt','w')\n",
    "for i in range(0,xl.size):\n",
    "    print(\" %18.9f %18.9f\"%(xl[i],yl[i]),file=f)\n",
    "f.close()\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}