Created
April 18, 2016 20:27
-
-
Save jrleeman/8bcafffd512feaa77ec2946dc92460ae to your computer and use it in GitHub Desktop.
Shear heating and diffusion modeling.
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| { | |
| "cells": [ | |
| { | |
| "cell_type": "markdown", | |
| "metadata": {}, | |
| "source": [ | |
| "# System Thermal Model\n", | |
| "\n", | |
| "We model the system as two steel sideblocks, an acrylic center block, and two layers of quartz representing the gouge. \n", | |
| "\n", | |
| "## Thermal Properties\n", | |
| "* __Steel__ $\\kappa$ = 1.172x10$^{-5}$ m$^2$/s\n", | |
| "* __Acrylic__ $\\kappa$ = 1.12x10$^{-7}$ m$^2$/s\n", | |
| "* __Quartz__ $\\kappa$ = 1.4x10$^{-6}$ m$^2$/s\n", | |
| "\n", | |
| "\n", | |
| "\n", | |
| "## Model\n", | |
| "* 1-D Model\n", | |
| "* Fixed room temperature at the right and left boundaries\n", | |
| "* Uniform temperature anomaly in the gouge layers at the beginning of simulation\n", | |
| "\n", | |
| "$\\frac{\\partial T}{\\partial t} = \\kappa \\frac{\\partial ^2T}{\\partial x^2} $\n" | |
| ] | |
| }, | |
| { | |
| "cell_type": "markdown", | |
| "metadata": {}, | |
| "source": [ | |
| "## Admin Code" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 1, | |
| "metadata": { | |
| "collapsed": false | |
| }, | |
| "outputs": [], | |
| "source": [ | |
| "from fipy import *\n", | |
| "import matplotlib.pyplot as plt\n", | |
| "import numpy as np" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 2, | |
| "metadata": { | |
| "collapsed": true | |
| }, | |
| "outputs": [], | |
| "source": [ | |
| "plt.rcParams['figure.figsize'] = 10, 7.5\n", | |
| "plt.rcParams['xtick.labelsize'] = 16\n", | |
| "plt.rcParams['ytick.labelsize'] = 16\n", | |
| "plt.rcParams['font.sans-serif'] = 'Helvetica'\n", | |
| "plt.rcParams['text.color'] = 'black'\n", | |
| "plt.rcParams['font.size'] = 18\n", | |
| "plt.rcParams['lines.linewidth'] = 2\n", | |
| "plt.rcParams['axes.facecolor'] = 'white'\n", | |
| "plt.rcParams['axes.edgecolor'] = 'black'\n", | |
| "plt.rcParams['axes.labelcolor'] = 'black'\n", | |
| "plt.rcParams['axes.axisbelow'] = True\n", | |
| "plt.rcParams['xtick.major.size'] = 10\n", | |
| "plt.rcParams['xtick.minor.size'] = 5\n", | |
| "plt.rcParams['xtick.major.width'] = 1\n", | |
| "plt.rcParams['xtick.minor.width'] = 1\n", | |
| "plt.rcParams['ytick.major.size'] = 10\n", | |
| "plt.rcParams['ytick.minor.size'] = 5\n", | |
| "plt.rcParams['ytick.major.width'] = 1\n", | |
| "plt.rcParams['ytick.minor.width'] = 1\n", | |
| "plt.rcParams['xtick.color'] = 'k'\n", | |
| "plt.rcParams['ytick.color'] = 'k'\n", | |
| "\n", | |
| "%matplotlib inline" | |
| ] | |
| }, | |
| { | |
| "cell_type": "markdown", | |
| "metadata": {}, | |
| "source": [ | |
| "## Setup Model Grid" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 3, | |
| "metadata": { | |
| "collapsed": false | |
| }, | |
| "outputs": [ | |
| { | |
| "name": "stdout", | |
| "output_type": "stream", | |
| "text": [ | |
| "Model Length : 0.146000 m\n", | |
| "Model Spatial Step : 0.000100 m\n", | |
| "Number of Cells : 1460\n" | |
| ] | |
| } | |
| ], | |
| "source": [ | |
| "# Setup Model Grid\n", | |
| "L = .146 # Model length in meters\n", | |
| "dx = 0.0001 # Spatial step in meters (0.1 mm)\n", | |
| "nx = L/dx + 1\n", | |
| "mesh = Grid1D(nx=nx, dx=dx)\n", | |
| "\n", | |
| "print \"Model Length : %f m\" %L\n", | |
| "print \"Model Spatial Step : %f m\" %dx\n", | |
| "print \"Number of Cells : %d\" %nx" | |
| ] | |
| }, | |
| { | |
| "cell_type": "markdown", | |
| "metadata": {}, | |
| "source": [ | |
| "## Set Spatially Varying Diffusion Coefficients" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 4, | |
| "metadata": { | |
| "collapsed": true | |
| }, | |
| "outputs": [], | |
| "source": [ | |
| "# Setup Diffusion Coefficients\n", | |
| "D = FaceVariable(mesh=mesh, value=1.172e-5) # Set all to Steel Values\n", | |
| "\n", | |
| "X = mesh.faceCenters[0]\n", | |
| "\n", | |
| "# Set 50> and <=53 mm as quartz\n", | |
| "D.setValue(1.4e-6, where=(X > 0.050) & (X<=0.053))\n", | |
| "\n", | |
| "# Set >53 and <=93 mm as acrylic\n", | |
| "D.setValue(1.12e-7, where=(X > 0.053) & (X<=0.093))\n", | |
| "\n", | |
| "# Set 93> and <=96 mm as quartz\n", | |
| "D.setValue(1.4e-6, where=(X > 0.093) & (X<=0.096))" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 5, | |
| "metadata": { | |
| "collapsed": false | |
| }, | |
| "outputs": [ | |
| { | |
| "data": { | |
| "text/plain": [ | |
| "(0, 146.00000000000003)" | |
| ] | |
| }, | |
| "execution_count": 5, | |
| "metadata": {}, | |
| "output_type": "execute_result" | |
| }, | |
| { | |
| "name": "stderr", | |
| "output_type": "stream", | |
| "text": [ | |
| "/Users/jleeman/anaconda/lib/python2.7/site-packages/matplotlib/font_manager.py:1282: UserWarning: findfont: Font family [u'sans-serif'] not found. Falling back to Bitstream Vera Sans\n", | |
| " (prop.get_family(), self.defaultFamily[fontext]))\n" | |
| ] | |
| }, | |
| { | |
| "data": { | |
| "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnEAAAHzCAYAAABPHUt2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xe4a1WZ+PHvy6UzIE0FVJoy/gZRUEdlLICoiChlRnQU\nBRS7Y28jYzv23h07KqLYG4oOolwUBRELFlQEBEGQInBBpHPf3x9rh5vk7lOSnZycJN/P8+TJyc5O\n1jrr7H32u1eNzESSJEnjZY1RZ0CSJEm9M4iTJEkaQwZxkiRJY8ggTpIkaQwZxEmSJI0hgzhJkqQx\nZBAnLVBEbBsRKyNi5ajzMk4i4lNVub1mhHk4r8rD7rO8v1dEfD8iVrT+xhFxaNv7W0TExyPigoi4\nqXp/+eL9Bh15nfN36fM7W7/z1oP6zmkUESd2HzsD+M6Z6js/Oajv1ORYc9QZkBZLg+DrB5n54LbX\nTq7Yn57KLSL2AE6oeesa4GrgXOBnwLcz8/gFpr9aHiLiQcB3gABuBi6p9ru2en/NKh//r9p+JXAj\ncHkvv8+A1f4uA/jODhGxC3AAcG5mHjng9OYVEecBrcDyN5m58zz7vxt4ftumbTPz/CFlby7D+B/h\n/x2txiBO06R1ce62GbAWcD2woub9UV6sVVwG3FL9vC6wBbAV8ADg+RFxDvD0zJytduxsSlB2bc17\nz6MEcF8AnpSZN3S9/3BKAHc5sGtmntPkFxmAuX6Xfp1JOTdu6tq+C/Bq4ERg0YO4LnePiJ0z81d1\nb1bB9kFtmwx6NPEM4jQ1MnPLuu0RcSKwG/D5zDxsUTOlhUjgPu01KhGxBrAz8CjgOcCdge9FxFMy\n81OrfUHmQ+f4/rtVz0fVBHDt7y9fAgHcfL9Lv9/5L4P+zgE7n1Ijdwjw4ln22Ru4LfBnYJtFypc0\nUvaJkzR2MnNlZv4yM18P7AScSqlN+0hE3KPHr1uvev5Hn+9r+D5LCeYfXwXwdQ6pnj+zOFmSRs8g\nTupTROwUEZ+PiIsj4vqI+H1EvDIi1prncw+sPveXiLghIi6PiOMj4nGz7L9H1bH53Or1IyLiOxFx\nabX9+dX2VqfqQyJio4h4W0ScExHXVc+vj4h12773IRFxXET8LSL+ERE/iIjdZsnDGlW6H4mIn0fE\nJRFxY0RcFBFfjYgH131uMWTmZcB/AFdRmsVf3b1P3WCAtkEqrVqb5W0d/Je3BmQArQEZh7a9f+sg\ngIV0Zp9rMEJE7B4RX66Ohxsj4qqIOCsivh4RT4+IWMDv8rFq25fmKquIOLza7xdd21cb2FD97p+o\nXu4Rnb/7yirfu1U/Xx8Rm86R7vZtn9thrjzO4nzgB5Rm9L1qvn9jYF9KoP3l+b4sIm4fEe+MiD9E\nxLVVmZ8aES+KiLXn+ezeEXFC9ZmrI+KUiHjiQn6JiFg7Ip4TESdFxBXV+f/niDgiIv7fQr5DamcQ\nJ/UhIvYCfgo8BlgbWAbcFXgd8MU5PvdW4IfAY4EtKf2abgM8BDg6Io7uvmi3yYh4MXAs8LAqzZtZ\nve/PplXeXkJpXgpgO+AVwNerfDwXOB7Ys/qedYEHAcdH6ejfbccq3adR+kmtTelDeHtKx/fvR8TL\nZ/u9hy0z/wq0Ru/tGxHr1+1GZ1ldAlwMtAa8XFG9vpjS/21F9XOrBu76tvcvZlUfvfbvnzOb3ftE\nxNOB5ZQgdEvgBsrfa3tgP+DDlMB0vu/6bPW8T0RsOEceHt+1/1z5v4QSGEPpK3dx1+OGzPwh8EfK\n8XAQs3ty9fyjzDxrjv3mclT1fEjNe48F1qEc39fM9SURcV/gd8ALgR0og1TWBO4DvAM4NSJuO8tn\nXwp8G9gD2IBSLvcBPh0R75gn3S0p5+X7KH05NwSuA+5IKZ9fRMS/z/UdUjeDOKl3AXwe+AawXWZu\nSgnEDqdcCPePiEes9qFSY/ZSygXwacDGmbkJ8E/A46rtjwP+e5Z0twDeAvwvsGVmbka5EHyla7/X\nUAKTB2bmRtU+T6MEfHtFxJuBdwJvAjar8rAdcAolYHhXTdo3AEdQakFuk5mbVN+9BfAqSkDzxuoC\nOSrfqZ7XBP6t5v2O4Dgzt8zMrYC/VJv+IzO3qh4HZuYLqvdbF+fPt72/VWZe2GP+umvU1qf8HaCU\n7daZuWFVrpsBjwCOpj447A70fwBcSGn6rQ0EIuJulKbnlcDn5sts1Yf0BdXLH3f97ltl5k+q9z5e\nPT959W+5tf9iq5byE3X7LEACX6IE0vvVBKqtwO7TrF427XnZhBLobQL8GrhvZm5MOUceQxl5vDM1\nQW5EPBB4a/XyKGCr6hzcDHgb8KLqs3XprkX5f3EP4HuU43PdKu07AO+h3EgdFRHbz1oKUheDOKk/\nP83Mx7c622fmtZn5VkptFcCB7TtXzT1voNx5Pzwzj8jMv1efvT4zv0ipjUngpVHfJLsu8IXMfG7V\nhEhm3lATTKwPPCozT672uSkzj6Bc4KAEiZ/OzFdm5tXVPuezqpbmXtE1X1hmnpWZT8vM72XmNW3b\nL8vMNwKvpVw8nzl/0Q3Nb6vnALYd4PfOGhQ0tBOlNucaysjaW/+OmbkiM4/LzCdmZveI0dVkZlJG\n18LsNWKtv+9JmXnRAvO4kN/9SMoNwi5R3x/xoZTapr8zRy31fKrj7muU4/vW86sKeu4PXEQJkOby\nHMqNx5XAXpn58+q7V2bmVyg3UQAPreki8Nrq+YTMPDQzL60+e1VmvpwSiN9mlnQPBf6VUgv/iMw8\nNTNvqT5/cWa+CPhI9bu9cJ7fQbqVQZzUu6TUiNX5evV8t67tj6ZcsL+Xmb+p/dJSs3EesDFw71nS\nffsC8velzPxTzfbWBS6BN9ekfz5l+oqoyf98vlU937/Hzw3SFW0/z9o/awlpNVWuBWw+gO87unre\nc5bmwMd17TcQ1Q3FMZTjpm50d6uG7kuZ2XRalNaNSHuTauvnz1bB7Fxawd/HW0FYu2q+wVOql49t\nba/6+z2Ycu68tftzlTfNkW6rJvK9reCtRuvvMvDRx5pcBnFSf06bZXurhmOTru2t4OYhUQZC1D4o\nNRYB3Knmu6+bbY6sLrVBImWuNYDr55gq45LqeePuNyJivYh4YdWR/9JYtXLBSqDVUX6rBeRvWIZV\nYzYsZ1WPdYBTIuIFEXHXfr8sM39Bme9tTdoCEICIuB+ln92NlGbJQftY9fyE9lrkqvnyAErwc8QA\n0jme0u1gt4honSMHV9//6Vk/VfKyNqX2Myn9EGfTmmD6nm3bWj+vBH5U96HMPJdVTfPt6a4JtLoZ\nfGyOc/+r1T6umqEFM4iT+pCZs003cX313N0c2pqjbn3KYIPZHmtSLjLrsbqFTjr811m2t2oALpnl\n/fZ9OvJfdco+ndKHazdKP6DrgEspF9VWgLjBAvM4DO2B8xWz7rVEZOZKStPnhZQA613A76OMVv5i\nROzbx9e2+ro9vmt76/VxmVk3oXVT36WMIN2MMkq05SBKkHpmZp5S98FeVGX2WUrAfnDVT2074PTM\nPGOej2/KqkB/rv6MrffaazNbP1+Vmdct4LPd6bbOp02Y/dxv1R6v2/0F0mwM4qTF0TrX3pOZyxbw\nqKtVmK0ZZjG8hzKS7xxK371NM3OjzNyi6vxfN5BgsbX3x6prTl5yqj5ZOwBPpNQknUOpBT0Q+EZE\nHBuzz4tWp9Ukt2tEbAO3Dix4bNf7A1U1Y7YGLbQPcGj9PMh1P1vnxsF0DmjoxWIGSq2/XwL3nOe8\nXyMzly1i3jTmDOKkxdGq/Rq7meSrZqj9KRehJ2Tm1zPzqq7dtlj8nK2mNSL4RuAnc+04YDdXz3MF\nBrdhlilIqoEtR2fmkzJzB8rqE2+u9n8EPQwWycyzKevJrsGqPnB7UP4+11BGSA7LJyh5fnhEbFEN\ncrgXpXx6DbJmVfUp/TVlSp9Dq+9fSHB6RZW/YO7z8I7V82Vt21r9524TEXW15C113QkuZ9U0NmN3\n/mtpM4iTFsfJ1fPu0Tbh7pjYnDIPGMAvZ9lnpJ2xI2Ir4EnVy2/M0+Q1aK3mybp+jETEXZh91OJq\nMvO8zHwFq0ab1k7APIdWQPP4rudvZOb1NfvPpRV8zNvfMDP/AhxH6RJwCKtq4b6TmXM14fejFRSu\nBXy3NVp7nvzdyKr+onNNTr1n9dw+IXLruF8DeGDdhyJiO2qOgWp08WmUMlxt6iGpCYM4aXF8iTJp\n7KbUrCjQruoMvpT8ve3n1aaQqPrLPXfxsrNa+rejdArfiDKf3esWOQu/rp73m+X91iTI3fPEzbmy\nB6v6V67TY34+T6lxuntE7EIZGQ39NaVeXT2vNtBlFq0BDoexaqqTQQxo6HYUpX/mO4A39vC51moO\nT4qI1WqPq0m8d6WU363ToWTmlcD3KX/Dl83y3XNNdv2ptnTnXBaumo5IWhCDOGkRZOYVlMmAAV4e\nER9tX34oItaPiAdHxEeBH48kk7Oo5rM7hXIB+0RE7Ay3LsX1EMpEs8PWHQAti4hdIuJVlNqV+1L6\nDD5tjg7uc00/Md/UFHP5MquCpvdExG2qPN4uIt5H6e9WN7XGI6slm54anctdrR8RTwOeUG06rpf8\nZubFlNGXQZmId2Pgb5TBB71qzb13twVO5PxNSteBf6Z01r+EVdPPDEw1P+FLM/NlPQ6Y+ABl4M96\nwP9FxL3h1uPp0ZQAGMpUQCd2fXaGUu4PibIk2+2qz94mIt5EmVC7u5tByxGUJv51gROqv/mtExZH\nxFYRcWhEnAQ8v4ffR1NuSQdxEXHHiHh/9Y/u2uha26/a56FRlir6U7XP2RHxwdmWTZFGJTM/QFnd\nIIGnAmdGxDURcQWlv9L3q+291rwshhdSRqPeHfhlRFxDyfPxlBF3Txli2gH8rG06hhWU5Y5+QZmA\n9baUqToekplzLX4+V5Ng39OTZObvKAM/AJ4HXBkRV1JG7T4TeDoliFrto8D9gI8C51X/v1rHwkco\nTYXHVu/3mt9Wrdu9qucvzTE/2azfWfWx+yGlifQn1cjZ8yLi3Graku79b6ZM/ttyVDWidBTqfp8V\nlClPrqTUKp8WEVdTyvxLlID3V6wKoNs/+2NWraZyCHBx9fe6nFIL907KCO7VVOWyP+UGbVPK33RF\nVZ7/oExN8knKVESjKi+NoSUdxAF3oSyFcjnlH0mdp1OGtb8BeDilQ/B+lH84o5zuQONjtTUtG37X\n7G+W1Q12pvwT/2O1eT3KP/H/oyzL1b126ULzNt/vsZDvqf2OzPwpZQTq1ykdxJdRgpQPU9ZSnWv+\nun7LN9ueN2XVVAxrUKZy+BHwXuBhmXnXLOt49pOHhZTbfH/XFwPPppTDdZRawe8Ae1Yjjeu+4wTK\nCMsjKU2y11CmaLmMUmt2cGbuWxMELaQ8v0JpWm7tu5Cm1Nm+8z+AD1JG/K5P6fe1NbPfbLQGT7SP\nWO1Xk3NztoEkp1HWAn43ZV69ZZTBMKdR1hu+X2bWBd1k5jso/dqWU5qa16Csh3pwZr50nnQvA3an\nBIjfptRSbkA5Vn5POQ4ey+qTCQ/y/5MmTMw/wfXoRES0ZuCOiKdSLnzbVjPLt/bZvPuEi7KA9w+A\np2TmIIe2S5LmEBGvAF4P/CQzR7mChzTxlnRN3AKWUGGWO6afVc+jnD1ekqZKRCyjdAmA+mZgSQO0\npIO4Bnavnn8/0lxI0pSoJhV+NWUutItZtXqEpCFZc9QZGLRqxM97gN+xajFySdIQRMSulFGdmwAb\nUvpv/U9m3jDSjElTYKKCuGqh4c9R1ql8QN2oqIhYup0AJWn8BfDJiLA/sjQgmVk7In1igriqKv9I\nymzbj8zM386271IezLEURYRl1gPLq3eWWW8sr95YXr2zzHozzPKKmH1GoYkJ4ihTHTwWeHRmLh91\nZiRJkoZpIoK4iHgnZbLRQzLzmFHnR5IkadiWfBAXEQdWP967et4nIv4GXJqZP4yI/6bMJv8J4Oyq\nk23LpZn5p0XMriRJ0qJY0pP9AkRE++CEZNVSKidm5p4RsRzYjfplaD6VmYd1fd9Cpp9TG/tG9Mby\n6p1l1hvLqzeWV+8ss94Mu0/cbAMblnwQN2gGcb3zZO6N5dU7y6w3lldvLK/eWWa9GVUQN6mT/UqS\nJE00gzhJkqQxtOQHNgzDJptsMuosjB3LrDdLsbw233xzvv3tb7PDDjuMOitST37zm9+w//77c+WV\nV/b9HUvxnFzqLLPejKK8pjKIW7FixaizMHYss94sxfJasWIFJ5100pIM4l7zmteMOgtjZdrKa/ny\n5Zx77rmNvmMpnpNLnWXWm1GU11QGcVdcccWosyAtquc+97l89rOfHXU2ZjUzMzPqLIyVaS2vpz3t\nabz1rW8ddTakRbXpppvO+t5UBnFWEWvarLPOOqPOgtTYuuuu6/9vqY0DG6Qp4pQBGkcet1I9gzhJ\nkqQxZBAnTYGI2nkipbHicSx1MoiTJEkaQwZx0hSxb5HGkcetVM8gTpIkaQwZxElTwL5EmgQex1In\ngzhJkqQxZBAnTRH7FmkcedxK9QziJEmSxpBBnDQF7EukSeBxLHUyiJOmiM1SGkcet1I9gzhJkqQx\nZBAnTQGboTQJPI6lTgZxkiRJY8ggTpoi9i3SOPK4leoZxEmSJI0hgzhpCtiXSJPA41jqZBAnSZI0\nhgzipCli3yKNI49bqZ5BnCRJ0hgyiJOmgH2JNAk8jqVOBnHSFLFZSuPI41aqZxAnSZI0hgzipClg\nM5Qmgcex1MkgTpIkaQwZxElTxL5FGkcet1I9gzhJkqQxZBAnTQH7EmkSeBxLnQziJEmSxpBBnDRF\n7FukceRxK9UziJMkSRpDBnHSFLAvkSaBx7HUySBOmiI2S0nS5DCIkyQtad58SPUM4iRJY8HmVKlT\nT0FcRKyMiFsaPl49rF9GUj0vfpI0edbs4zPnA3/uM73d+vycpAGwWUrjyONWqtdPEPfJzHxdP4lF\nxMp+PidJkqROo+gTZ7uOtMhsTtUk8DiWOvVaE3c74B8N0mv6eUmSJNFjEJeZf2uSWNPPS2rGvkUa\nRx63Ur1++sTNKyI2Bx4EXAt8LzNvGUY6kiRJ06pRn7iIeFZEnBoRm7ZtuzfwB+ArwHeAUyJig2bZ\nlNSEfYk0CTyOpU5NBzb8J0BmXtG27e3AxsAngGOBfwWe1TAdSQNgs5QkTY6mQdwOwK9aLyLitsDu\nwCcy86mZuS/wM+DxDdORJE0pbz6kek2DuM2AS9teP4AyhcjX2radBGzbMB1J0pSzOVXq1DSIuxLY\nvO31bsBK4OS2bQms2zAdSQ148ZOkydM0iPsdsG9EbB4RGwOPA07LzKva9tkGuLhhOpIGwGYpjSOP\nW6le0yDuvcCWwAXAX4AtgA927bMrbf3mJEmS1FyjIC4zjwGeSamROxN4cWYe1Xo/Ih4MbAgc18/3\nR8QdI+L9EXFKRFwbESsjYuua/TaJiI9HxGURcU1EHB8RO/X3W0mTx+ZUTQKPY6lT48l+M/OjwEdn\neW85ZbqRft0FeAxlhOsPgb26d4hyVn8T2Bp4DrACOBxYHhG7ZOaFDdKXJorNUpI0OXquiYuI90bE\n7rE4t0Q/yMwtMvNRwJdn2Wc/4P7AwZn5hcw8rtq2BvCyRcijJGmIvPmQ6vXTnPocYDlwSUQcERGP\nioi1B5wvAHJhZ+5+wIWZ+YO2z11NqZ3bfxj5kiRJGrV+gritKP3gfg48ETgG+FtEfCEiHhcRGw4y\ngwtwN+C3Ndt/B2wdEesvcn6kJce+RJoEHsdSp56DuMy8JDM/mpmPAG5LCeSOA/YBjgYujYhjI+Kp\n1QoOw7YpZb66bq2lwDZZhDxIY8FmKUmaHE1Hp16dmUdn5mMoAd3+wOeB+1EGO/w1In4YES+IiG2a\nZ7c+G0P6XknSEuDNh1Sv6Txxt8rM6zPzm5n5ZOD2wEMoc8ZtC7wLODcifjGo9NpcSamN67Zp2/sd\nImLex8zMzBCyKknql82pmiQzMzMLikfmMrAgrl1m3pKZyzPzeZQVG+4HvIXhLL91BqVfXLcdgT9n\n5rU1+Zv3YRCnSeLFT5KWlpmZmQXFI3MZShDXLovTMvN/MnPHISRxDHCHiNittSEiNgL2rd6TVLFZ\nSuPI41aq13iy32GLiAOrH+9dPe8TEX8DLs3MH1ICtVOAz0TES1k12W8Cb1vs/EqSJC2GJR/EAV9s\n+zlZtTbricCemZkR8SjgHdV76wInAw92tQapsDlVk8DjWOq05IO4zJy3yTczrwSeUj0kzcJmKUma\nHEPvEydJUhPefEj1DOIkSWPB5lSp04KDuIjYKiIOjIid27ZtExG7RsQ/DSd7kgbBi58kTZ4FBXHV\n9B1nUQYZ/DIi3l69dTGwJXDVcLInaZBslpKkybHQmrhXAocCGwM7AVtExFsy8wbK9B7e5kuShsKb\nD6neQoO4kzPzy9Vaqb/LzIOBsyLiybh2qSRpEdgtQOq00CDuaoCI2L61ITOPAC4FHjmEfEkaIC9+\nkjR5FhrE/Tgi3gycHRG7tjZm5rHAOcA1w8icpMGyWUqSJseCJvvNzFMj4jfA5zLz113v/SAidhlK\n7iRJU8+bD6negldsyMxrgV/P8t6fBpYjSZJq2C1A6uRkv9IUaF38rNGQpMnReO3UiLgT8EJgZ+CO\nwFp1+2Xm9nXbJUmaizcfUr1GQVxE7AF8B1gHuJkyWvXmml09AyVJjdicKnVqWhP3dkqT7CHA0Zm5\nsnmWJA2aFz9JmjxNg7idgM9n5mcGkRlJw2WzlCRNjqYDG1YAlw8iI5Ik1fHmQ6rXNIg7Fth9EBmR\nJGkudguQOjUN4g4HNomID0bEBoPIkKTB8+InSZOnUZ+4zLwsIvYBfgIcHBF/BK6aZd89m6QlqTmb\npSRpcjSdYmQn4ERgo2rTPZtmSJKkdt58SPWaNqe+C9gEeDWwDbB2Zq5R92icU0nSVLNbgNSp6RQj\nuwJfy8w3DCIzkobDZbckafI0rSG7CTh3EBmRJKmONx9SvaZB3HLgvoPIiCRJc7E5VerUNIj7b2DH\niDg8PLukJcvTU5ImT9M+ca8Efgu8EXhqRJzO7FOMHNYwLUkN2SwlSZOjaRB3aNvP21WP2RjESZJ6\n5s2HVK9pELf9QHIhSdI87BYgdWq6YsN5A8qHpCHy4idJk8dJeKUpYrOUJE0OgzhJ0pLmzYdUzyBO\nkjQW7BYgdTKIk6aAy25J0uQxiJMkLWnefEj1DOIkSWPB5lSpk0GcNAW8+EnS5DGIk6aIzVKSNDkM\n4iRJS5o3H1K9Ba3YEBGbZuYV1c93APYFlgFfy8yLhpg/SZIAuwVI3eatiYuIFwF/jYhtI+JhwG+B\nVwP/A/w+InYfch4lNeTFT5Imz0KaU68C/rlaJ/VhwB0yc6vMvANwF+CgIeZP0gDZLCVJk2MhQdwW\nmfnn6uefZea1rTcy8zLg90PJmSRJePMhzWYhQdwFEfGmiLgTcHFE3D8i1omI9SNiL+BuQ86jJEl2\nC5C6zDuwITM/HRGHAb8ANut6+0TgkCHkS9IAueyWJE2eBY1OzcxPRMRngPsDW1H6yZ2emRcOM3OS\nJEmqt6AgDiAzb6TUvEmStGhaNcg2p0qd+prsNyIeN+iMSJIkaeH6XbHh3gPNhaShsk+cJE0el92S\nJC1p3nxI9QziJEljwT5xUieDOGkKePGTpMljECdNEZulJGlyGMRJkpY0bz6kegZxkqSxYLcAqVO/\nQdyfBpoLSUPlFCOSNHn6CuIy80ODzogkSZIWbsHLbs0lItYCNgaWAX/PzH8M4nslSXLZLalev8tu\nrRURz4yI70fECuAG4BLgIuDqiLi8eu9ZETGQQHGe/DwoIo6PiEsj4uqI+HlEPHnY6UqSJI1KzwFW\nRNwO+D5wV+APwKnAP6rHWsAG1WMH4H+BZ0fE3pl54aAy3ZWfewLHAz8GngJcCzwGOCIi1snMDw8j\nXWmc2CdOkiZPP7VkbweOBj6cmVfOtWNE3Bl4LvBe4MA+0lqIx1bP+2bmtdXP34+IewCHAAZxkjTG\nvPmQ6vXTnPq3zHzzfAEcQGaek5kvAObdt4FlwE3AdV3brwbsQCFJE8I+cVKnfoK4fs6iv/fxmYU6\nArgFeF9EbBkRG0fE04A9gXcPMV1pbHjxk6TJ008Qt01EHBQRC/psRDwR+Oc+0lmQzDwTeDilH9yF\nwBXAB4BnZOYXh5WuNI5slpKkydFPn7j/oQxseFdEnAqcT2m6vBFoXSE2AbYG7kcZ5LBH45zOIiJ2\nAr4F/Ax4P6VZ9QDgIxFxQ2YePay0JUnD582HVK/nmriq5utulGbMfwH+CzgceA0wUz2eB2zT2icz\nfzWY7NZ6PbCCMrDh25m5PDOfD3yRMqBiNREx72NmZmaIWZYk9cpuAZokMzMzC4pH5tLXHG6ZeRXw\nCuAVEbEhJWDbiNI37TLgkkWc8HdH4NeZeXPX9tOAgyLidpl5afsb3tVp2jjFiCQtLTMzMwuqMJor\nkGs8EW9m/h34bdPvaeAvwM4RsVZm3tS2/X6UptUrRpMtSZKk4elrxYZeRcRBQ/z69wHbA9+MiP0i\nYq+I+ADwOOBDNTV0kqQx4rJbUr1FCeKApw7rizPzG5TRqWsDHwe+DNwfeDbwsmGlK0mSNEqLsa7p\nesBOw0wjM4+nLL0lqYZ94iRp8vSzduqvKaNTYf6Jf7PaxyuHJKkv3nxI9fqpidsV+AFwM3DmAvZf\nl1Xrm0qS1Bf7xEmdeg7iMvPaiHg28OzMfPJCPhMRd5t/L0nD4sVPkiZPXwMbMvM04LY9fOTCftKR\nNFg2S0nS5GgyOvUlPex7cIN0JElTzJsPqV7fQVxm/qGHfS/rNx1JksBuAVK3xZonTtIIOcWIJE0e\ngzhJkqQxZBAnSVrSXHZLqmcQJ0mSNIYM4qQpYJ84SZo8BnGSJEljyCBOkrSk2SdOqrcoQVxErLUY\n6UiSJE2LRkFcRHwsItadZ5/tgJOapCOpGfvESdLkaVoT9xTgtIj4l7o3I+JA4JfAfRumI0maUt58\nSPWaBnFvBHakBHKHtTZGxDoR8SHgi8DNwAEN05EkTTn7xEmdGgVxmfkqYC/g78DHI+IzEXEf4KfA\nM4AfA7svzVIVAAAgAElEQVRk5jGNcyqpbzanStLkaTywITO/D+wCfA84CPgJcDdKLd3umfmXpmlI\nkiSp06BGp14DXFr9HMBVwImZuXJA3y9JmlJOMSLVaxzERcQuwC8otXDfBZ4FrA0cFxFvigjnopMk\nSRqwplOMPBc4BdgOODwz987MjwD3An4NvBz4UURs3TinkvpmnzhJmjxNa8neC1xC6fv21tbGzDwL\n+Dfgf4FdgdMbpiNJkqQ2TYO4Y4B7ZuYp3W9k5g2Z+Vzg0YC3/5KkvtgnTqq3ZpMPZ+a8879l5tci\n4udN0pEkSVKnRRl0kJnnL0Y6kurZJ06SJk9PQVxE/D4int1vYk0/L0maPt58SPV6rYm7K7B5g/Sa\nfl6SNKXsEyd16qdP3B59nkiefdKI2JwqSZOnryCuevTLYE6SJKmhXoO4PQeQ5rkD+A5J0pRwihGp\nXk9BXGaeOKR8SJIkqQeuaypNAfvESdLkMYiTJEkaQwZxkqQlzT5xUj2DOEmSpDFkECdNAfvESdLk\nMYiTJC1p3nxI9QziJEljwT5xUqdGQVxE3H5QGZE0PDanStLkaVoTd0FEfDEiHjKQ3EiSJGlBmgZx\nZwIHAsdHxB8j4iURsdkA8iVJEuAUI9JsGgVxmXl34IHAp4E7Am8D/hIRR0fE7gPInyRJkmo0HtiQ\nmSdn5pOArYDnA2cDjwOWR8TvIuIFEbFJ03Qk9c8+cZI0eQY2OjUzV2Tm+7tq57YB3kWpnTsyIu4z\nqPQkSZKm2bCmGLkcuAK4vnq9LnAwcGpEfCMiNh1SupKkCWOfOKnewIK4iFg7Ih4fEScCvwdeAFwK\nvBDYDHgocBywL/DBQaUrSZI0jdZs+gURsQPwdOBJlGDtZuCrwAcz84S2XU8AToiILwMPb5qupIWz\nT5wkTZ5GQVxEnADsUb28CHgt8LHMvGiOj/0c+I8m6UqSJE27pjVxewDLKc2jX8/MWxbwmW9RAj5J\nkuZlnzipXtMg7lnATzLzV7PtEBF3B+6ZmZ8GyMzfAL9pmK6kPticKkmTo+nAhg8B+8+zz/7AJxum\nI6kBazAkafIMa4qRdssWIQ1J0oSyOVWqtxhB3A7AlYuQjiRJ0tTouU9cRHwSSKB1S3RARGxbs+sy\nyooNDwKO7TN/kgbAKUYkafL0M7Dh0K7Xu1SP2fyEMuGvJEmSBqSfIG57VtXE/Ql4L/AeVtXMtdwC\nXJmZ1zTKoSRpqtknTqrXcxCXmee1fo6I1wHLM/PPg8xUPyJiH+DlwD2BlcAfgZdl5vKRZkySJGkI\nGs0Tl5kzA8pHIxHxDOD91eO1lP54OwPrjTJf0lJhnzhJmjyN104dtWpQxXuAl2Tm+9re+u5IMiRJ\nkrQIegri2kamHp6Zl7S9nldmHtZH/hbiMOBm4MND+n5J0gjZJ06q12tNXGtk6luAS1h9pOpchhXE\nPRA4EzgoIl4FbA2cB7w7Mz84pDSlsWRzqiRNjl6DuO2r5790vR6lrYAtgbcBhwPnAI8FPhARa3Y1\nsUpTyRoMSZo8PQVx7SNT616PyBrAhsChmfn1atuJVV+5wwGDOEkaYzanSvUaLbsVERsPKiMNXE7p\nl3d81/bjgdtHxO27PxAR8z5mZmYWIeuSJGkazczMLCgemUvTtVMvjogvRsQjI2Ix1mGtcwarTzQ8\np8yc92EQp0niFCOStLTMzMwsKB6ZS9PA61zgQOCbwIUR8Y6IuHvD7+zVV6vnvbu27w1ckJmXLHJ+\nJEmShq5REJeZ/wLsCnwIWBt4EfCriPhFRDw/IjYfQB7ny8O3geXARyLiGRGxV0R8DHgY8Kphpy9J\nGi77xEn1GjeBZuZPM/O/KCNEHwN8C7g78G5K7dzXI+Lfm6YzjwOAz1NWa/gmcB/goMz89JDTlSRJ\nGomBrdiQmTcCXwG+EhG3Aw4CDgH2A/alLIU1FJn5d+A51UNSF/vESdLkGdZghMuA31WPm+hx4IEk\nSZLmNtC1UyPiXyirODyRMgkvwNnAkYNMR5I0PewTJ9VrHMRFxKbA4ynB279Wm68GPg58KjNPbpqG\npMGwOVWSJkejIC4ivgrsQxmZupIywe6ngK9l5vWNcydpIKzBkKTJ07Qm7gDK4vNHAkdl5oXNsyRJ\n0io2p0r1mgZx98/MnwwkJ5IkSVqwppP9GsBJY8ApRiRp8vRUExcRn6QsNn94Zl7S9npemXlYH/mT\nJElSjV6bUw+tnt8CXNL2eiEM4iRJPbNPnFSv1yBu++r5L12vJUmStIh6CuIy87y5XktamuwTJ0mT\nZ1jLbkmSJGmIBrrsVktEbA48CLgW+F5m3jKMdCRJk88+cVK9RjVxEfGsiDi1Wnqrte3ewB+ArwDf\nAU6JiA2aZVPSINicKkmTo2lz6n8CZOYVbdveDmwMfAI4lrKe6rMapiOpAWswJGnyNA3idgB+1XoR\nEbcFdgc+kZlPzcx9gZ8Bj2+YjiRpynkzInVqGsRtBlza9voBQABfa9t2ErBtw3QkSVPKbgBSvaZB\n3JXA5m2vdwNWAie3bUtg3YbpSBoAL4aSNDmaBnG/A/aNiM0jYmPgccBpmXlV2z7bABc3TEdSAzZD\nSdLkaRrEvRfYEriAsorDFsAHu/bZlbZ+c5Ik9cIpRqR6jYK4zDwGeCalRu5M4MWZeVTr/Yh4MLAh\ncFyTdCRJktSp8WS/mflR4KOzvLecMt2IpBFy2S1JmjwuuyVJkjSGGtXERcS5lNGnc+4GZGZu3yQt\nSdJ0sk+cVK9pc2pUj24bAxtVP18E3NQwHUkDYHOqJE2ORkFcZm4723sRcRfgfcAGwN5N0pHUjDUY\nkjR5htYnLjPPBh4N3AF4zbDSkSRNB29GpE5DHdiQmdcB36NMAixJUs/sBiDVW4zRqTdTJgSWNGJe\nDCVpcgw1iIuI2wIHUFZ0kDQiNkNJ0uRpOsXIa6ifYmRNYGtgf+A2wOFN0pEkTS+nGJHqNZ1iZL4B\nC1cDr8/MtzZMR5IkSW2aBnF7zrJ9JXAl8PvMvLlhGpIactktSZo8PQVxEfF84JTM/ClAZp44jExJ\nkiRpbr0ObHg3bRP3RsTKiHj1YLMkSdIq9omT6vUaxN0ArDOMjEgaPptTJWly9BrEnQs8PCK2GEZm\nJA2HNRiSNHl6DeI+DNwLuCgibqm2zUTELXM8VrbtK0lSX7wZkTr1NLAhM98XEZcCjwK2AvYA/lw9\n5vxoX7mTJE09uwFI9XqeYiQzPw98HsrABuBTmfnaQWdM0uB5MZSkydFTc2pEPD8i7tu26bXAiQPN\nkaSBsxlKkiZPoylGKCs27D647EiS1MkpRqR6TjEiSZI0hpxiRJoCLrslSZPHKUYkSZLGkFOMSJKW\nNPvESfWcYkSaIjanStLk6LU5tdvrcIoRacmzBkOSJk/PNXHtMnNmQPmQJGlO3oxInZrWxEmSNFR2\nA5Dq9VQTFxGfpAxSODwzL2l7Pa/MPKyP/EkaIC+GkjQ5em1OPbR6fgtwSdvrhTCIk0bEZihJmjy9\nBnHbV89/6XotSdJQeTMidep1nrjz5notSdKg2Q1AqtdodGpLRGwLbE7pH3dZZp4/iO+VNFheDCVp\ncvQ9OjUibhsR746IvwLnAKcCPwXOjYiLIuIdEbHpoDIqqX82Q0nS5OkriIuIHYCfAc8Hbg+sBC4D\n/lb9vAXwIuDnEWG/OUlS31x2S6rXcxAXEWsAnwXuRFmt4aHABpm5RWbeHtgQ2Av4AbBNte+iioj/\ni4iVEfH6xU5bWspsTpWkydFPTdxewL8CXwIekpknZOaNrTcz8/rM/B7wEODLwP0iYq+B5HYBIuLx\nwD1a2VmsdKWlzBoMSZo8/QRxjwZuBJ6bc9zWZ+ZK4DnATdVnhi4iNgHeBbxwMdKTJC0eb0akTv0E\ncfcCfpyZl863Y7XPj6rPLIa3Ar/JzC8sUnqSpCGzG4BUr58g7k7Ab3vY/wxK37ihiogHAgcD/zXs\ntKRx5cVQkiZHP0HcRsCKHvZfUX1maCJibeAjwNsz86xhpiWNI5uhJGny9BPErQ3c0sP+K6vPDNPL\ngHWANy5k54iY9zEzMzPUDEuSeuPNiCbJzMzMguKRuQxkxYZ5DLX9JiK2Bl4BPAVYLyLWa3t73Yi4\nDfD3aqBFyZBNSpI0NvyfrUk0MzOzoAqjuQK5fldseE1E3LKQB/AahhvIbU+phfsMcEXbA+AlwJXA\nTkNMXxobXgwlaXL0WxPXa532MOvAfwnsUZPecuAo4AjKsmDS1LIZSpImT89BXGb2vd7qMGTmVcAP\nu7dXF60/Z+Zq70mSxofLbkn1llRAJmm4bE6VpMmxGAMbRmKp1RhKo2QNhiRNHgMdSdJY8GZE6mQQ\nJ0la0uwGINUziJOmiBdDSZocBnHSFLAZSpImj0GcJGkseDMidTKIkyQtaXYDkOoZxElTxIuhJE2O\nnuaJi4hP0uc6qJl5WD+fk9SczVCSNHl6nez30AZpGcRJknrmsltSvV6DuO2HkgtJi8LmVEmaHD0F\ncZl53pDyIWmIrMGQpMnjwAZJ0ljwZkTqZBAnSVrS7AYg1eu1T1ytiNgKeAiwFbBO3T6Z+bpBpCWp\nf14MJWlyNA7iIuJ1wMsX8F0GcdKI2AwlSZOnUXNqRDwBeCXwQ+DAavORwBOAjwIrgS8AD26SjiRJ\n3oxInZrWxD0LuBB4RGbeVJ1g52bm54DPRcTXgG8Dn2uYjiRpStkNQKrXdGDD3YFvZ+ZNbduWtX7I\nzOOA44CXNExH0gB4MZSkydE0iFsL+Fvb6+uA23Tt81tgl4bpSGrAZihJmjxNg7iLgS3bXl8A3KNr\nny2BmxumI0mact6MSJ2aBnG/BHZqe/19YLeIOCQiNoiIR1EGPPyyYTqSBsDmVI0jj1upXtMg7pvA\nThGxXfX6rcAK4FPA1cAxQFBGsEqSJGlAGo1OzcxPUQK21uvzI+K+wIuAuwDnAh/MzN80SUdSMzZD\naRJ4HEudBrJiQ7vM/BPwnEF/ryRpOtmcKtVz7VRpingxlKTJMYhlt+4EvBDYGbgjZdqR1WTm9k3T\nktQfm6EkafI0CuIiYg/gO5RF728GLqV+OhFv/yVJjXgzInVqWhP3dkqT7CHA0Zm5snmWJElaxW4A\nUr2mQdxOwOcz8zODyIyk4fJiKEmTo+nAhhXA5YPIiKThsRlKkiZP0yDuWGD3QWREkqS5eDMidWoa\nxB0ObBIRH4yIDQaRIUnDY3OqxpHHrVSv6YoNl0XEPsBPgIMj4o/AVbPsu2eTtCRJkrRK0ylGdgJO\nBDaqNt2zaYYkDZ7NUJoEHsdSp6bNqe8CNgFeDWwDrJ2Za9Q9GudUkjSVbE6V6jWdYmRX4GuZ+YZB\nZEbScHkxlKTJ0bSG7Cbg3EFkRNLw2AwlSZOnaRC3HLjvIDIiSdJcvBmROjUN4v4b2DEiDg/PLknS\nENgNQKrXtE/cK4HfAm8EnhoRpzP7FCOHNUxLUkNeDCVpcjQN4g5t+3m76jEbgzhpRKwol6TJ0zSI\n234guZAkaR7ejEidmgZxWwNXZ+bpg8iMpOGyOVXjyONWqjeI0alPH0RGJEmStHBNg7jLgesGkRFJ\nw2MzlCaBx7HUaRA1cfcfREYkSapjc6pUr2kQ9yrgrhHxhohYaxAZkjQ8XgwlaXI0HdhwOGWeuP8B\nDouIXwEXA6tdKZwnThodm6EkafIMcp64LarHbAziJEl982ZE6uQ8cZKkJc1uAFK9RkFcZp43oHxI\nWgReDCVpcjQd2CBpDNgMJUmTp2lz6q0iYgPgrsAGmXnSoL5XkiTwZkTq1rgmLiLuFBFfBVYAPwNO\nbHvvQRHxu4jYo2k6kpqzOVXjyONWqtcoiIuILYGfAPsB3wJOAdpvlU4Fbg/8Z5N0JEmS1KlpTdxr\nKEHaXpn578Dx7W9m5o3AScADGqYjqQGboTQJPI6lTk2DuH2AYzLzhDn2OR/YqmE6c4qIAyPi6xFx\nfkRcGxF/iIg3RcQ/DTNdSZKkUWkaxN0e+OM8+9wEDDuYenGVzsuBvYEPAc8Cjg9v3aRb2bdI48jj\nVqrXdHTqlcCd5tlnB8pSXMP0qMy8vO31DyPiCuBIYA9g+ZDTlyRJWlRNa+J+BOxXDXBYTUTsQKkZ\nG2oQ1RXAtfyseh5qU640DqyQ1iTwOJY6NQ3i3g6sB/wgIh5R/UxE/FNE7EMZsZrAOxum04/dq+ff\njyBtSdKA2Jwq1Wu67NapEfF04MPAsW1vXUWZauQm4LDM/G2TdHoVEXcAXgccn5m/WMy0paXMi6Ek\nTY7GKzZk5ici4keUgQT/BmxGCeJOAT6QmWc2TaMX1YjUbwA3Ak9ezLSlpcpmKEmaPANZOzUz/5iZ\nL8zMXTNzh8z818x87ggCuPWAbwLbAg/PzItm2W/ex8zMzCLmXJI0H29GNElmZmYWFI/MZWBrp45a\nRKwFfBm4F/CwzDxjtn1tUtK08tjXOPK41SSamZlZUIXRXIHcINZO3SMijo2ISyPipoi4peuxMiJu\naZrOPHlYA/gsZTqRAzLzp8NMT5IkadQa1cRFxCMp/c/WAC6gTPx7c82uw76N+l/gQOCNwHURsWvb\nexdk5oVDTl9a0myG0iTwOJY6NW1OnaGMQN0/M7/bPDt925sSKL6ierSboYxUlSRJmhhNg7idgC+M\nOIAjM7cbZfrSuLBvkcaRx61Ur2mfuH8AdaslSJIkaYiaBnHfo8wNJ2kJsy+RJoHHsdSpaRD3cuDO\nEfGq8OySljybpTSOPG6lej31iYuIT7L6SNMzgNcCT46I04EVdZ/NzMP6yqEkSZJW0+vAhkPneG/b\n6jEbgzhpRKwo1yTwOJY69RrEbT+UXEiSJKknPQVxmXleRBwCnJ6Zvx5SniQNiX2LNI48bqV6/Qxs\n+BRwwIDzIUmSpB40XjtV0tJnXyJNAo9jqZNBnCRJ0hgyiJOmiH2LNI48bqV6/a6dunFEbN3LBzLz\n/D7TkiRJUpd+g7gXAM9f4L5BmSB4WZ9pSWrIvkSaBB7HUqd+g7irqsdCWRcuLQErV67kxhtvnHe/\nZcuWsWyZ910aroUci1COW0mr6zeIe09mvnagOZE0dCeffDLrrLPOvPutv/76HHvsseyxxx7Dz5Sm\n0pOe9CSOPPLIUWdDGmv9DmywZk0aI3e/+925y13uwlprrTXvIyK49tprOeWUU0adbU2w7373uwCs\nueaaCzoud9xxR+5617uOONfS0tJvTZykMbLxxhtz1llnLWjfV7ziFbzpTW+yCUtD1Tq+LrjgArbY\nYosR50YaT04xIqnDGmuUfwsGcRqm1vHVOt4k9a7fs8chQtKEMojTYjCIk5rruTk1Mz3jpAnWuqje\ncsstI86JJlnr+DKIk/rn2SOpgzVxWgzWxEnNefZI6mAQp8VgECc159kjqUNrkl+DOA1T6/hyUmmp\nfwZxkjpYE6fFYE2c1Jxnj6QOBnFaDAZxUnOePZI6ODpVi8HRqVJznj2SOlgTp8VgTZzUnGePpA4O\nbNCwZSaZZQlugzipf549kjpYE6dhawVwEUGECwBJ/TKIk9TBIE7DZlOqNBieQZI6OLBBw+agBmkw\nPIMkdbAmTsNmTZw0GJ5Bkjo4sEHD5moN0mAYxEnqYE2chs2aOGkwPIMkdTCI07AZxEmD4RkkqYNB\nnIbNIE4aDM8gSR0cnaphc3SqNBieQZI6WBOnYbMmThoMzyBJHRydqmFzdKo0GAZxkjpYE6dhsyZO\nGgzPIEkdDOI0bAZx0mB4Bknq4MAGDZsDG6TB8AyS1MGaOA2bNXHSYHgGSergwAYNmwMbpMEwiJPU\nwZo4DZs1cdJgeAZJ6mAQp2EziJMGwzNIUgcHNmjYHNggDYZnkKQO1sRp2KyJkwbDM0hSB4M4DZtB\nnDQYnkGSOjg6VcPm6FRpMAziJHWwJk7DZk2cNBieQZI6GMRp2AzipMHwDJLUwdGpGjZHp0qD4Rkk\nqYM1cRo2a+KkwfAM0rxmZmZGnYWxMu7lNYqBDeNeZott3MtrsQc2jHt5jYJl1ptRlVdk5kgSHpWI\nyGn7nZuKCCyzhRv38jr99NO55z3vyc4778zpp5++KGmOe5kttnEvr5NOOonddtuNBz7wgZx00klD\nT2/cy2sULLPeDLO8qu+OuvesiZPUweZUDZvNqdJgTMQZFBF3iogvR8SKiLgqIr4SEXcadb6kceTA\nBg2bAxukwVhz1BloKiLWB04ArgMOqTa/AVgeEffIzGtHljlpDLUurDfccAMXXnjhoqW7mGlNgnEu\nr0svvRQwiJOaGvsgDngasB3wz5n5J4CI+DVwFvAM4N0jzJs0dlqdzc855xzueMc7Llq6i5nWJJiE\n8nLFBqmZSQji9gNOaQVwAJl5XkT8GNgfgzipJ3e+85152MMexhlnnLFoaV500UVstdVWi5beuJuE\n8lq2bBkHHXTQqLMhjbWxH50aERcDX8vMZ3Vt/yBwYGbermu7o1N75Cil3lhevbPMemN59cby6p1l\n1htHp/ZvE+DKmu1XVO9JkiRNnEkI4iRJkqbOJPSJu5L6GrdNKbVxq4morZXUHCyz3lhevbPMemN5\n9cby6p1l1ptRlNckBHFnADvVbN8R+F33xtnalSVJksbJJDSnHgPsGhHbtTZExLbA/av3JEmSJs4k\njE5dH/gVZbLfV1abXw9sADjZryRJmkhjXxNXBWl7An8EjgI+A5wD7GkAJ0mSJtXYB3EAmXlBZh6Y\nmbfJzI0y8z8y8/zW+66tOruIODAivh4R50fEtRHxh4h4U0T8U9d+m0TExyPisoi4JiKOj4i6vohT\nJyL+LyJWRsTru7ZbZm0iYp+I+GFE/L06D0+LiAe3vW95VSLiQdXvf2lEXB0RP4+IJ3ftM5XlFRF3\njIj3R8Qp1f+slRGxdc1+CyqfiFg3It4eEX+tvu/kiHjQ4vw2w7eQ8oqIh0bE0RHxp2qfsyPigxFx\n25rvm+jygoUfY12f+XC131E17w2tzCYiiJtLrFpb9Z8pa6seDOxAWVt1/VHmbYl4MXAT8HJgb+BD\nwLOA46MaalM9fxPYC3gO8GhgLUoZ3mEUmV4qIuLxwD2ql9m23TJrExHPAL4OnAYcADwG+CKwfvW+\n5VWJiHsCx1P+Pz8F+HdKuR0REc+s9pnm8roL5fi5HPhh3Q49ls8RwFMp3XEeCfwVOC4idh5K7hff\nvOUFPB3YjLLu+MOBN1NWQ/pJRGzQte+klxcsrMxuFREPAJ4AXE3bdaDN8MosMyf6ATwfuBnYvm3b\ntpTA5YWjzt+oH8BmNdsOBlYCD65e71+93r1tn40oB/h7R/07jLDsNqlOxv+syud1be9ZZqt+720p\nfVafN8c+lteq3/vNwPXA+l3bTwZOnvbyourLXf381Koctu7aZ0HlA+xc7Xdo27ZlwB+Ab4z6d13E\n8tq85nMPqvZ98jSV10LLrO39tYDfAv8NnAt8uuv9oZbZxNfEMcvaqkBrbdWplpmX12z+WfXcWpxx\nP+DCzPxB2+euptzpTnMZvhX4TWZ+oeY9y2yVwyg3Uh+eYx/La5VllJvM67q2Xw20pkia2vLK6io4\nj4WWz36Usv5C2363AJ8HHh4Raw0k0yO0kPLKzL/VbO6+DsAUlBcs+BhreSnlvHwnq87PdkMts2kI\n4u5GiZK7/Y4yl5xWt3v1/Pvqea4y3Hoam6Uj4oGUGsv/mmUXy2yVBwJnAgdFxDkRcVNEnBURz27b\nx/Ja5QjgFuB9EbFlRGwcEU+jDOB6d7WP5TW3hZbP3YA/Zeb1NfutTWlWm1bd1wGwvDpExF2AVwDP\nzsybZ9ltqGU2DUGca6v2oOov8jrg+Mz8RbV5U2YvQ5iycoyItYGPAG/PzLNm2c0yW2UrSj/UtwFv\nAh5G6fP1gYh4XrWP5VXJzDMp/ZIeA1xIKYMPAM/IzC9Wu1lec1to+cy336YDztdYiIgNgfdQAo2v\nt71leXX6EPCVthrfuhq8oZbZJKzYoAGJMiL1G8CNQPtIuPGeTHDwXgasA7xxjn0ss1XWADak9Alp\nXRBOjDIp9+HA+0aUryWpGkH5LUpz1vspzaoHAB+JiBsy8+hR5m9MeP71KSLWBD4HbAk8IDNXjjhL\nS1JEPBG4N3DXUeZjGoK4ntdWnUYRsR6lv8i2lM7AF7W9fSX1dwubtr0/Faph5q+gjBpcryq3lnUj\n4jbANVhm7S4H7kypfWt3PLB3RGyB5dXu9cAKYN+2JprlEbEZ8N6I+ByW13wWWj5XAnVTR7T2m6pr\nRESsARxJabp/ZGZ2N0lbXtxa4fEuSuvCTRGxcfXWMmDt6jrwj+r8HWqZTUNzak9rq06jqmPll4F7\nAftk5hldu5xBadfvtiPw55yuSZW3p9TCfYZy8rUeAC+hnLA7YZm1O4P6Dr/d+1hexY7Ar2v62JxG\nmQbidlhe81lo+ZwBbBcR69bsdyNw9vCyuCR9GHgs8LjMXF7zvuVVbF493kTndeCOlPK7Etin2neo\nZTYNQZxrq86huvP6LLAHcEBm/rRmt2OAO0TEbm2f2wjYl+krw19Syqr90Zqw9qjq9dlYZu2+Wj3v\n3bV9b+CCzLwYy6vdX4Cda0at3Y/StHo5pduD5TW7hR5Px1CmiHhs235rUqYNOi4zb1qc7I5eRLyT\n0sLwpMyc7RiyvIq/Uv7v70HndeASSgvDHpQZMGDYZTbq+ViG/aBMJnoW8GvKUN/9KGutnk3XPEzT\n+KB0zFxJacLZtetxh2qfqA7I86sD7+HAicDfWvtM+4PV54mzzDrL5/vV7/4MygSsH6vK7BDLa7Wy\nas1x9n/V/6u9KAMbVgLvsLwS4MDq0fr/9czq9W69lg+l/9cVlADmIZRWiWuBXUb9ey5ief13tf3j\nlJuF9uvA9l3fNfHltZAym+Uz59E1T9ywy2zkBbVIf4w7VYV2FWWupa8yy8R90/agTE54S3WQdj9e\n3bbfJpSpDy4H/kG527j7qPO/VB50BXGW2Wrls2EViFwM3ACcTmmysbzqy+thlJVmLq3+Z/2iuois\nYbFyFMUAAAqxSURBVHnder61Hu3/v07otXyAdSlzfP2VUtN5ylwX6nF8zFdewPI5rgOfmLbyWugx\nVvOZ1Sb7HXaZRZWAJEmSxsg09ImTJEmaOAZxkiRJY8ggTpIkaQwZxEmSJI0hgzhJkqQxZBAnSZL+\nf3vnHmxVVcfxz1dEKRVRxEdQoPlAe2g+MpRByAfpKGYaY+oImvaczClMSwNMRw3zUTaW4gM0cyQt\nHMsXBtcHXp+Jo4GISoJoqAiCoiTcX3/81q7NZp97zuUeuJ7r7zOzZ13W+q21fmvvcznfu9b6rR00\nICHigiAIgiAIGpAQcUEQBEEQBA1IiLggCNqMpMGSWiSNaWc7I1M7I+rlWyt99Ut9Xb+u++psSJqQ\n7l12ndnRPtWKpFEF3+P5B52GEHFB0ADkvoBWSdqhFbtpOdt1LoyAer3ypaZ2cqIvfy2TNF/SFEnn\nStqxHn2V9F0X4drgXA6MBR7sYD/awnTc51+nf8drioJOw4Yd7UAQBDWzEv+d/SZwdrFQ0k7AATm7\nzvxlNQOYnH7+GLA1/rLunwNnS7oCGGVmq3J1XgH64+9Qbg+d+b5W43Izm9fRTrQFM2sGmiX1BX7Y\n0f4EQT0JERcEjcNC/AXKJ0kaXRAoAKek9A7gqPXq2fpnhpn9opgpaTAwAf+y7gZ8Nyszs5XA83Xo\nW3VoI1j/xHMLOh2xnBoEjYMB44FtgcPzBZK6AiPxpaOZlRqQtJOkGyQtkLQipRMrLUFK2kbStZIW\nSlou6SlJJ7bmpKQtJV0oaVaqs0TSfZIObuN424yZNQFDgf8A35K0R86v0j1xaYy/kjRb0juSFkt6\nTtL1krZPNhOAqanKmMJy7qBk013SGZKmSnol3d/XJd0u6Utl/qb60yT1lHS1pNckvS/pWUkjK41T\n0iGS7kjtvy9pnqTJkg4ssR0q6U5JbybbFySNk7R5W+5ta0hqSmPZUNJoSS9Kei/dx1Nzdt+T9Ez6\nXMyXNFaSCm397zlJ+rSkWyUtkrRU0r2SPpvsekm6Jt2z9yQ9nkR8EHxkiJm4IGgsbgYuxWfdbs/l\nDwN6AWcAO5dVlLQPcB+waao7E9gVOAE4UtJBZvZEzn4r4GFge3wP1EPAJ4DfA1Mq9NEXaAL6Ag8A\nd6b+DgfulvRtM7tmLcZdM2Y2W9IkfFzH4Uuvq5nk/P04Lnx3AO7F74uAfvg9/RMwF/hLqjcCH19T\nrr1/pXQ34Hzgfnw2dDF+H4YBh0o6wszuKXG5R/JhBTAJ2BgYDlwnqcXMbsgbSzoXXzZehi8pzwd6\nA/sBxwN/z9mOAcYAi5JPrwO7A6OAwyQNMLNlJT6tLbcAXwT+BnwAfB24StKq1O/xyY8pwJHAaGA5\nMK6krX7AI/jn9Dr8c3gU0CRpIP7ZWoz/TvQEjgXukrSzmc2v45iC4MOLmcUVV1wf8gtoAealn8fj\nX5C9c+V3419o3XAh0QKcmCsXMAtYBXyj0PbwZD8LUC7/6pR/ScF+L3ymqwUYXShrwvfkDS/kbw48\nhX9hb53LH1n0tcp9yOyvq2J3crKblsvrV6wLHFE2xlS2IbBp7t+Dy8acK+8ObFmS3xtYAMys8Fxb\n0r3O3/td0zP+Z8H+kGT/ArBdWV+5n4ck24eA7gW7Eans0hrv+4Rk/6kK5U2p/NF8X7jwWoHvQ3wx\n73P6TLyBC8suJc+pBfhpoZ9zUv4S4MpC2Qmtjans+ccVV6NfsZwaBI3HeKALLlSy2a+DgZvM7P0K\ndfYDdgGazezmfIGZTcK/6HcBBqY2u+KzJkvxyL68/ZPATcUOJO0ODAJuS23m67yd2ukGHF3zSNee\nV1Paq0b7Ne6bma00s3dq7dDMlprZWyX5C4DbgP6S+pRUfRf4kZlZrs4sfBa0f5otzPhBSn9sZq9V\n6CvjtJSeamZLC3YTgafxZ1xPzsr3ZWZz8VnGzYDz8j6nz8Rfga3wGd4ic4GLCnkTU9oFn3XO80f8\nD4jd2zOAIGgkYjk1CBoMM3tM0jPAyZLOx5dWhYu7SuyZ0qkVyqfhAm4PfOm0Px71+biVL7fdj8/m\n5BmQ0h6SxpbUyQTVrq34WS+yfVbVIkmb8FmysyTtCdyFC9oZZtbS5k6l/fGgigH4eDcqmPTGo2Tz\nzKkgFufj49gCn8EEj8BtwWdeqzEAn80bXtx3ltgI6CVpCzNbXEN71TDgiZL8TFA/WVKWic4++Hjz\nzMgL20QmAp83s3dX69ysRdLrqa0g+EgQIi4IGpPxwG+AQ4GTgCfM7OlW7LNN7GvM3hTyexTsF1aw\n/3dJXs+UHpyuMgzYpEJZPclmdt5ozcjMlqWgg3PxvWtDU9Gbkq4EzjePaq2KpKOAW3HBNQVfPnwX\nF11D8ONfNi6puqRCk1m/XXJ5PYDFZraiBpd6prqtnWtn+J7Feog4Kgj+bBxlR7tkZV1LytawN7OV\nSY9WOiZmZYW2gqBTEiIuCBqTG4FfAlfhgmVsFfvsS2/bCuXbFeyydJsK9mXtZHVOM7PfVvFnXTMk\npY9WM0xLkKcASNoN+DLwfXzT/QYprYXz8GXZvc1sdr5AUm9cxLWXJcAWkrq1snSe8TaAmW1Vh36D\nIPgQEnvigqABSfuJbsWX597BI/Ra4x8pHVKhfEjB7jngPWAPSd1L7AeX5DWndFAVX9YpkvrjUZEt\n+D6pmjGzmUmAZjOJR+aKs3P5ulDOjnjwQlHAbUDaa1gHmvH/t79So+2WSZgGQdAJCREXBI3LOcBX\ngaHF/UFFzGw6MBsYKGm1wAJJx+AiY7aZPZTsPwD+gEdcji3Y703JhvgU8PAg8DVJJ5X5IelzkmoN\nNmgzkg7A94t1BX5nZs9Usd9NUtlsYzbTuDyXtyilfSs0NxfYWVI2q0naizYW3wdYjzc9XJHSSySt\nEQxQyLsspePzPuVsN5G0bx18CoKgg4jl1CBoUMzPwmrLeVgj8L1at0i6HRd1u+BCcClQPMT3Z8CB\nwOlJuE3Hl12H4+eADSvp4zg8eOJaSacBj+FLgH2AzwOfwTfnt7pXrQa+kAue2Bhf9t0XF0urgEuA\nn9TQziHAxZIeBubgx130wWfgVgEX52yfwzfiHyvpA2AeLsxuMH8V1WX4GXpPSfozHlSwf/LpDvw4\nk3ZhZlNSMMs5wCxJk/FAiW1wId6M75HEzKZKOgu4EJgj6U78TLtNcSE6CBfdh7XXr0S8ESEI1jMh\n4oKg82GUzPqkqNZ9cAFwEC4q3sCPCznPzOYU7BelaMsLku3euJD5DvAyJSLOzBZI2gs/CuNoXNR1\nwQMnZuIvIX+2mq9VxgYuCLOjJJYDbyXfJgE3mtlLNbZ3N/BJXNAMw2ceXwXuwc8beyQ3tpYUvHAR\nvly7WSp6AD/D72pJK4DTcUG8HBdJI4BjKLxlo8axlj3H0ZKa8SNEDscDRRbikaETC7bjJE1PtgNx\ncboEF6NX0cbl5rb62s6ytfUjCD4yaM0I7iAIgiD4P/LXjp0IbG9mL3ewO2uFpH7AS8AEMzu5Y70J\ngvoQe+KCIAiCamR/7c9N7zU9s0O9aQOSRklqwQVcEHQqYjk1CIIgqMZkPHAj48GOcmQtmM7qwTnF\nd+kGQcMSy6lBEARBEAQNSCynBkEQBEEQNCAh4oIgCIIgCBqQEHFBEARBEAQNSIi4IAiCIAiCBiRE\nXBAEQRAEQQPyX8C8RCqlTBjUAAAAAElFTkSuQmCC\n", | |
| "text/plain": [ | |
| "<matplotlib.figure.Figure at 0x10b199b90>" | |
| ] | |
| }, | |
| "metadata": {}, | |
| "output_type": "display_data" | |
| } | |
| ], | |
| "source": [ | |
| "# Setup figure and axes\n", | |
| "# Generally plots is ~1.33x width to height (10,7.5 or 12,9)\n", | |
| "fig = plt.figure(figsize=(10,7.5))\n", | |
| "ax1 = plt.subplot(111)\n", | |
| "\n", | |
| "# Set labels and tick sizes\n", | |
| "ax1.set_xlabel(r'Model Distance [mm]', fontsize=20)\n", | |
| "ax1.set_ylabel(r'Thermal Diffusivity [x $10^{-6}$ m$^2$/s]', fontsize=20)\n", | |
| "\n", | |
| "# Plotting\n", | |
| "ax1.plot(X*1000., D.value*1e6, color='k')\n", | |
| "\n", | |
| "plt.title('Thermal Diffusivity Model', fontsize=24)\n", | |
| "\n", | |
| "# Set limits\n", | |
| "ax1.set_xlim(0,np.max(X.value*1000.))" | |
| ] | |
| }, | |
| { | |
| "cell_type": "markdown", | |
| "metadata": {}, | |
| "source": [ | |
| "## Boundary Conditions" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 6, | |
| "metadata": { | |
| "collapsed": false | |
| }, | |
| "outputs": [], | |
| "source": [ | |
| "# Setup model boundary conditions\n", | |
| "valueLeft = 20.\n", | |
| "valueRight = 20.\n", | |
| "\n", | |
| "phi = CellVariable(name=\"Temperature\", mesh=mesh, value=20.)\n", | |
| "phi.constrain(valueRight, mesh.facesRight)\n", | |
| "phi.constrain(valueLeft, mesh.facesLeft)\n", | |
| "\n", | |
| "\n", | |
| "# Introduce anomalies on the layers for testing\n", | |
| "\n", | |
| "X = mesh.cellCenters[0]\n", | |
| "# Set 50> and <=53 mm as quartz\n", | |
| "phi.setValue(21, where=(X > 0.050) & (X<=0.053))\n", | |
| "\n", | |
| "# Set 93> and <=96 mm as quartz\n", | |
| "phi.setValue(21, where=(X > 0.093) & (X<=0.096))" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 7, | |
| "metadata": { | |
| "collapsed": false | |
| }, | |
| "outputs": [ | |
| { | |
| "data": { | |
| "text/plain": [ | |
| "(19.8, 21.2)" | |
| ] | |
| }, | |
| "execution_count": 7, | |
| "metadata": {}, | |
| "output_type": "execute_result" | |
| }, | |
| { | |
| "data": { | |
| "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnkAAAHzCAYAAABcyguCAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XmcLFV58PHfc+GKLCJ7FOIFcSOCG+4bXNxAA6gBkxgV\nRXCLRt83ikB4lZEIxiUuiOIuKCoqURYTl5twERUQiYmRzaDsFxfkXkBk5z7vH1V9p6ep7unq7pmu\n6fl9P5/+9EzVqarTp09PP3OeqlORmUiSJGmyLBl3BSRJkjR6BnmSJEkTyCBPkiRpAhnkSZIkTSCD\nPEmSpAlkkCdJkjSBDPIkqYEi4sqIWBsRuw+4/Vnl9q8cYZ1eVe5z5aj2OQm6tXVELC+XXzHgfnco\nt187mppqsTHI04LV+uM3wMMvqHlSfslNRcQLx12XudLWr0YWTLXJ8tF+vO3LNn1LjX2M2lD7jIhH\nRMR7IuL8iPhtRNwZEWsi4qcRcWxEPHlUFZ1n3drlXssj4kXl+9hPEO+EthrI+uOugDSE31L9x29L\nYClwO3Bjxfob5rJSmmE58E7gBOC0sdZkbt0rGBuBXwK3lo92D6Zo0yuBj/TY/irgUuCmEddrYBGx\nFPhn4G+ZHmS4h6KOmwCPAR4LvCkiVgD7ZeYt46jriPwR+AVwbcW6FwEHAGuB73fZ/s5ye4M8DcQg\nTwtWZj6wanlEnAXsBpycma+e10ppsYpR7zAzn9NtVZ/bz8XI4sAiYn3gW8BzKV7DycBxwI8z856y\nzEOAFwN/DzwH2AJYsEFeZv4E+LMhtr9umO0l07WS5sPIg6BFbKG25T9SBHhrgVdl5t9k5jmtAA8g\nM3+VmR8AHkIx+rsYLNT3UwuAQZ4WnYhYEhGviIgVEXF9eT7QdRFxckQ8qcs2U+V5V58vf39jRPxX\nRNwSEb+JiC9ExIPayj88Ik6MiGsj4vaIuDAiXttl3zNOzo6IfSJiZXmO0i0RcU5EvLSP1/WM8jVc\nGxF3RMQN5Wv86z6P+/yI+HZE/K5c/pa2srtFxEci4sdlW91Zlvt2ROxXse8dypPF31kuemXFuZHL\nqurRpa5dT/hv29/2EfFnZbtfExF3RcQ3O8puXZ4L9vOybf9YvjfvjojNZ2vjOqLjpPmI2KV8f35T\n9olLIuL/lSnMqu3vdeFFRFwJnFn+ukPcu01f2Va264UXEbFrRPxTRPwwIq5u6y8rI+KgiBjpd0NE\nPJBidA7gY5n5xV7lM/O2zDwIuKZiXw+JiE9GxOVlO66JiO/3qndbWxwQERuWn+dfRMRtZT/+SkQ8\ndJbXsFdEnBkRN0XEzRFxbkS8fJZt7tW3W8soUrUAR3a+j21lZ73wIiIeFxEnlX3+joj4fUR8JyL+\nosc2rb61W0RsEREfjIgryu1XRcSnIuIBvV6bFojM9OFjoh7AWRSjBZ+rWHc/YEW5fi1wN7CG4ryg\n1u9vrNhuqlz/eYo001rgNuDmtn1dCfwJ8HSKcwHXAqvLfbbKHFqx7+XlusuB/9NWjxuAu9q2/WiP\n1/zetnL3lK+p/bhfBqLHcd/acdw7gTeX5Tbp2PeNHW22FvhEx77/FPg18Idy/a3AdR2P7Trr0eP1\nvaosc2bFula9Xk5xDtTaso5/BL7RVu4Z5Wtrlb+trfxainPYHj5Af2ttf0DH8h3ajvW8sg3uKftE\n+/v6zS77vbIsv1vbsvOB37e9V51t+pKKz8EBFfv+fVvd/lC2S/v7+S1gvTrvwyxt9A/ldncADxji\ns713+b616r6a4tzbVr2/B2zU42/C3wE/beuTt7Rt+3tgxy7HPaStXOsz0vp8faBbW1f1beCpFJ+N\nW8t1f+h8H6v6UJd6vbbtfbuH6c9uq65fAJb06FsvK39u1ePWtm0vBzYb9L3y0YzH2Cvgw8eoH/QO\n8r5ZrvsJxTk/9ymXb1Z+Ed1R/vF+Wsd2U+V2ayhOEn8psH657hnlH+dWEHgNxUUGO5Tr7wd8nOnA\ncMuOfbe+CG4pj/95YOu2er2/7Q/vSyte01vKddcBBwH3K5ffF/jLtrod1uW4t1IEHR9tO+4GTAdh\nGwJfBfZt/6MP3J/iBPpWoLt/Rd2O7PZeVNRjmCBvbVmPM4FHtq3bsXzenunA9DjavsyBnYHvlPu4\nkIovxVn622xBXivY/wqwrFy3EXAo01/Qz6/Y75Xlut06lu8+W3t1fA6qgrwvlX1jm7ZlG1F86bf6\ny9vqvA+z1KX1j9U5Q3yuH8J0UHYm8LBy+X2A1zAd/H26R1usBn5FkTaOts/v1eX6r1Zs+4y29/HE\nVpuV/f+fmP670FeQ17bu8+W6d/Z4za0+dK8gD3haW//5KrBtuXxjir9lrXVH9Ohbq4H/BJ5cLl8P\n2KdcvhZ476Dvl49mPMZeAR8+Rv2gS5BHEdStBS6mDIQqtj20LHNGx/Kptj/0r6jY7uVt6y+hI1Cg\nOO/mf2f5IlgLfKdLvVpfCP/bsXwziv/A/wg8qsu2T2H6v/ylXY570hDt3XrtVQFYq93mI8i7DNig\ny/YnlWWO7rJ+KfDfZZn9ar7+foK8bu/r6eX6z1asa30RdwZ5s7ZXx+fgXkHeLNs9o9v+e70Ps+zz\n2nK744foZ59tfQaA+1asfw3TI1oP6dIWt1AxWgf8BdP/8CztWPcf5bp/71KvT/foA13fK4pzDocJ\n8lr1OpuOUfpy/dFM//Nzv451rb51HbB5xbZ/X67/1aDvl49mPDwnT4tJ69ykT2fmH7qU+XL5vDwi\nqk6Iviarzyf697afP5CZa9tXZvGXc2X5685djp3Ae7qsO7p8fkhEPKZt+X4U/7n/e2b+vHKnmedR\n/FHfDHh8l+O+v8tx+/Gt8vnJXdpsvhyXmXd0LoyIjYCXUHz5f6hqw8y8C/iX8tduV7UOKilGfKqc\nWj536xPzLjN/SDFavX15Lt0obFk+rx5k47Jftc79/FBm3l5R7DPAKop/qPbvsqtTMvPyiuWnl88b\nAOvOzYuILYA9KN7D93bZ5zG9az96HfV6T/n3pdN7KTIDmwAv6LKrT2XmmorlrX65Q0RsOGx9NT5O\noaLF5Gnl8zsi4tBZym5M8cX0+47lF3cpf335nBQpvyq/K58367L+LuBHVSsy85cR8RvgAcCuwM/K\nVa3X9OxyfTebU3z5PQg4r2PdbZn5s3tvMi2K6S9eSREsPYZiaovOCwbuWx5noC/yISVwbpd1j6eo\nawIX9ohDW19my0ZbNaA4PaDKdeXzSC/66EdEvIQiPbsrsDVFgNPpgRTnj43bjsCmFO/hyqoCmZlR\nTJ/0MuBxXfZT+T5k5t0R8TuKdmj/fLb2sxb4YZdtr4iIaynOQ50vrXolXebYy8ybI+ICinOEH0eR\n0u00W7+Eoj1uG7CeGjODPC0mrVGJzZh9rrGkOEepU+UXXmbe0xY8dPtSbE0VUXk1JfD7zLy7R51W\nUQR5W7Uta72mjZgOUrrJLmV6Tg4dEZsA36U4Yby1n9uYPg+Jsl5QBMfjCPJgOtDu1GqjoPgS76Vb\nGw0lM//YZVVrRKpbnxi5MmD/GsVkvFC85jso2q/VR7ehmH1h4xEd9gZgW4p/DgbR/r6t6lGuta7b\n+9xtBB+K9yKY+V609nNTZvYKdFYxv0Fee706J8tuN1B7ZObtbX/P5q1vavQM8rSYtE5PeFFmnt6z\n5MLRek0fzsy/71myu3tmWf8OigDveoqrcL+TmetGOMtpK1rB6TjTtd1eR6uNbszMQYOMSfIaigDv\nj8BhFFf3to/cEBHXANsxuvfzYoog77Ej2Nd96R2sLSZVo6/SOp6Tp8Xkt+Xz9mOtRXdblaMs3Wxb\nPrePWM3Ha3pJ+fx3mXlSe4BXGnY+rVaAeN8eZe4/xP5baexNI2LTIfYzKVrv5z9m5scqArz1KEaL\nZxvtrqOVYt11wPP8ftf2c6++3hpN6zaqO+hx7z/LuWnb9lg3F1r12jAitupRbtTtoQXGIE+LyTnl\n8/PHWovuljJ9jt0M5UStD6T44v1p26rWa9o9InoFScP40/K4/9Vlfa8LFVrp3F4jQq0Tv7fpNjEw\n8MQe28/mAopRviXAXkPspyn6adNeWl/83d7PpzP6EaITKOZvW0oxelhLebHETRSveY+qMuWI8vLy\n159WlRlAq42WUFx1XHXcB1Oc61rXMO9jq1692uP+TF9oNar20AJjkKfF5ITyec+I2LNXwYjodnHE\nXArg8C7rWssvy8z/aVv+dYq02xZM312ieueD39Gh9eX66Ip9bgIc0WPbm8vnXu35vxQBwBKKObo6\nj/FQpq+srC2LG9yfUv56VFnnShGxfkSM6jy0uXJT+Tzo6GZr+6r3c33g3a1fB9z/vWTmr4EPl7++\nKSIO6FU+IjaOiBMion3UrnX181u6jKodTDGitpbiczG08srT/6Boi7d3KVY7aC21Phu1P5dlvc4s\nfz20y1Xth1IE638A/m2gGmrBM8jTopGZ3wW+QfEH+5sR8bb2VEdEbBUR+0fEvwIfHEMVb6W4Svaz\nEbF1WafNIuK9wIEUo2lT7Rtk5mqmA8DDytsRPay1PiI2iog9IuJTdLlytw/fK58/GBG7te37iRRf\ngL3Oc2tdafzMbreNKqcvaU3Z8KGIeHoUlkTE8ygm0u11cnk/DqO4IOThwDkRsWdr1LA81k4RcQjw\nC+AJQx5rlKpSppdRpLg363Xrqh5a7+c7ImLfcgSMiNgJOINi1LTbhSLDOILpgOnzEfGliHha+ykK\nEfHQiHgb8EuK2361v/5jynptC/xrRDy83GaDiHgNcGxZ7rOZecWAdaxq76ly+bPLwHOb8rj3j4hj\nKM5xvKliu9m0Pht7DXgLsXdQBLS7AidHxHZlvTaJiH9gOvj8p/IfnU6jTMerqcY9UR/FfEanUsw4\nfitwKcWHeZO2MpswfeuY1uz6u/e5/4dTzOR/MdO3jzkNePS4X7uPOetTZ9H9jhcbUQR6a9sea5i+\n/Vbr8dmO7aa67bOtTGvbZV3WV+6DmbcXezPTE7quZuYti47tcewjmHlbqluYnrW+9fhVt+PO0p4P\npjgHqLWf25i+88AtTE8yfa/XTnFx12Vt639HMWffFZR31Gg7xvVt5f7I9B0M/pPidlSzTYZc2e5t\n5Z7A9KS8aylGD39PcWVpa9k9wDNr9rfZJkOuvCXVbO8BXSZDLted0HbcG9vadL+2Mmd1qdfmHe/J\nnRRBSuvnA7odmwEnQ27bfinwMWbe1q3qdlz3UPyd3rhj+72ZeeutNR3bfQ/YsMffhK4TQ8/S3m9r\nO8ZaZt6u8P0U5xzWnQy5NUVT6/X+uvU+9tuHKG5rdneXet1DcVuzqomSu77Wup8rH81+NGEk760U\nH/jDKM6XOR54A7CibQh6K4qRjDuZ/i+03/9CnkdxzsLnKFJBf0txOfl5EbHrKF6AGifp0j8y89bM\n/AuKL4tvUHzp35diVPsyirmkXkURVPS1z4pytevVVr9jKW4f1pr76laK+d9enplv7rHd0RTz132K\nIv0JxVQg11LcsusQ4Jk16tq+7yuAJ1HcNeK3FCMxq8vfn5iZrYmg77W/LKaEeTbwRYrbvd2f4vyl\nZRS3UGo/xpMpbv31u/IYV1OkDp/OdGqrazX7eB0XADtRpLHOKfe5KUWg+hPgIxT/PP5gtn0NcvwB\ntuvVX15PMXH2JRSBU6tN21PNldtnkep7CsXf2muYDqq/SfH6v9Dj2EON/mTmXZn5RmAX4H0U7f57\nin/kb6E4d+wjwBMy84XZMfVMZn4LeBTFXSauoPjs3gL8gGJEbc+snuqkn89vr78bH6A4l3clRb9Z\nQnEf4Vdk5iFt21fts/pgmTdQfDd9g+JztSXT72NfMvNTFCOvX6aYLmUjisD3exT3MT4gM7vVa9i/\nZ1oAovr9n8cKRGxZdvb2Za+guEfgszNzZce651B04OWZefaA+9+U4j+ZMzLzlZUbSvMkIpZTnF9z\nZWbuOObqSJImxNhH8joDsNIF5fPQl6VX7T8zb6YYtZnvy94lSZLmxdiDvC52L58vmYudl/f922Wu\n9i9JkjRujQvyyiuEjgJWZOZcze3zUYpzDT48W0FJkqSFqFG3NSvnrzqN4gKLA+foGIcDLwVencUE\nm1VlPNlU4/Bg+54kaRCZea/5EhsT5JWTW55Bccn47tlxq50RHeP1wNHAEZl5Qq+y474gZSGJCNur\nJtusHturHturPtusHturvrlss+r5sBsS5JWTkp5CManjczPzojk4xiso5mf6QGa+Z9T7lyRJapKx\nB3nlbOtfopg0cu/MPH8OjvFiinnyPp2Z3W5NI0mSNDHGHuRRjK7tT5FGvS0intK27prMXAUQEc+n\nmOjzUeW65eXtZf6Ymd9ubRARdwMnZObB5e+7UUyw+jPgxI7935GZ3W7SLUmStGA1YTLkKyhm+K5K\nKE9l5lFt5Vo3q8628jMmkI2ItRRB3qvL348EjuzYhqpt2/bRZZJwVfHcjPpss3psr3psr/pss3ps\nr/rm+py8qgsvxh7kNZFBXj1+2OuzzeqxveqxveqzzeqxveobR5DXuHnyJEmSNDyDPEmSpAlkkCdJ\nkjSBDPI0tCOPPHLcVVhwbLN6bK96bK/6bLN6bK/6xtFmXnhRwQsvJEnSQuGFF5IkSYuIQZ4kSdIE\nMsiTJEmaQAZ5kiRJE8ggT5IkaQIZ5EmSJE0ggzxJkqQJZJAnSZI0gQzyJEmSJpBBniRJ0gQyyJMk\nSZpABnmSJEkTyCBPkiRpAhnkSZIkTSCDPEmSpAlkkCdJkjSBDPIkSZImkEGeJEnSBDLIkyRJmkAG\neZIkSRPIIE+SJGkCGeRJkiRNIIM8SZKkCWSQJ0mSNIEM8iRJkiaQQZ4kSdIEMsiTJEmaQAZ5kiRJ\nE8ggT5IkaQIZ5EmSJE0ggzxJkqQJZJAnSZI0gQzyJEmSJpBBniRJ0gQyyJMkSZpABnmSJEkTyCBP\nkiRpAhnkSZIkTSCDPEmSpAlkkCdJkjSBDPIkSZIm0FiDvIjYPyJOjYirI+LWiLg0Io6JiE3aymwS\nER+IiLMi4uaIWBsRu9c4RkTE4RFxZUTcFhH/HRF/MTevSJIkqRnGPZL3VuAu4DBgL+B44A3AioiI\nssxWwIHAncD3ymVZ4xjvBo4Eji2PcR7w9Yh4/tC1lyRJaqjIrBMvjfjgEVtm5g0dy14BnAg8OzNX\ndqx7DkWgtzwzz+5j/9sA1wDHZOa72pb/O7B1Zj6my3Y5znaRJEnqV0SQmdG5fKwjeZ0BXumC8nnb\nERxiT2ApcFLH8pOAR0XE9iM4hiRJUuOMO11bpXW+3SUj2NfOwB2Z+auO5ReXz48cwTEkSZIap1FB\nXkRsBxwFrMjMn45gl1sAayqWr25bL0mSNHHWH3cFWsorak+juMDiwFHueoT7kibWAQccwHe+852+\nyz/1qU/l1FNPZfoaKWl0zjjjDN74xjdy++2391V+00035eSTT+YJT3jCHNdMWjgaEeRFxIbAGcAO\nwO6Zed2Idr0G2KxieWsEb3XFuladZt35kUceydTU1EAVk5rk7rvv5otf/GKtbU4//XRuuOEGttpq\nqzmqlRazU089lWuuuabv8tdffz3f+973DPI0UaampnjXu941e8Euxh7kRcRS4BRgV+C5mXnRCHd/\nEbBBRDyk47y81rl4F1dsA4BX12oxWrJkCb/+9a9nLbfTTjuxZk3VmRDSaH3wgx/kZS97Wc8yRx99\nNMcee+w81UiaP1NTU30NJnUbmBprkBcRS4AvAcuBvTPz/BEf4tsU8/C9jOJcv5aXAz/PzKtGfDxp\nQYsIttlmm1nLLVnSqNN5NcHuf//7z9onN95443mqjbSwjHsk72PA/sDRwG0R8ZS2dddk5iqAcuLi\njYFHleuWl3Pg/TEzv93aICLuBk7IzIMBMvP6iPggcHhE/AH4L+CvgD2Afeb2pUkLx6Aj1454a64M\n0rfsj9JM4w7y9qK4e8UR5aPdFNOjbx8HWnPaZbkO4Epgx7ZtlnDvK4aPAG4B3gI8ALgUeElm/tuw\nlZckSWqqsQZ5mfngEZe7Vw4pM9dSjBQeXa920uLT75WyXlGr+dJPX7M/StU8sUaS6Vo1julaaXgG\neZIkSRPIIE/SOqZr1TSma6XBGeRJMl2rxjFdKw3PIE+SJGkCGeRJWsd0rZrGdK00OIM8SaZr1Tim\na6XhGeRJkiRNIIM8SeuYrlXTmK6VBmeQJ8l0rRrHdK00PIM8SZKkCWSQJ2kd07VqGtO10uAM8iSZ\nrlXjmK6VhmeQJ0mSNIEM8iStY7pWTWO6VhqcQZ4k07VqHNO10vAM8iRJkiaQQZ6kdUzXqmlM10qD\nM8iTZLpWjWO6VhqeQZ4kSdIEMsiTtI7pWjWN6VppcAZ5kkzXqnFM10rDM8iTJEmaQAZ5ktYxXaum\nMV0rDc4gT5LpWjWO6VppeAZ5kiRJE8ggT9I6pmvVNKZrpcEZ5EkyXavGMV0rDc8gT5IkaQIZ5Ela\nx3StmsZ0rTQ4gzxJpmvVOKZrpeEZ5EmSJE0ggzxJ65iuVdOYrpUGZ5AnSWocU6/S8AzyJHlOniaC\n/VGaySBPUm2mx9Qk9kepmkGepHX8slRTtEbl7JPS4AzyJJmu1USwP0ozGeRJqs3RFTWJ/VGqZpAn\naR2/LNUUpmul4RnkSTJdq4lgf5RmMsiTVJujK2oS+6NUzSBP0jp+WaopTNdKwzPIk2S6VhPB/ijN\nZJAnqTZHV9Qk9kepmkGepHX8slRTmK6VhmeQJ8l0rSaC/VGaaaxBXkTsHxGnRsTVEXFrRFwaEcdE\nxCYd5TaPiM9ExPURcUtErIiIXfo8xtYRcVxEXF4e4/KI+GhEbDU3r0qafI6uqEnsj1K19cd8/LcC\n1wKHlc+PA6aAPSLiaZmZUXx6zwCWAW8CbgQOB1ZGxGMzc1W3nbdtuyPwDuASYGfgKOAJwFPn6HVJ\nC5JflmoK07XS8MYd5O2dmTe0/X52RKwGTgSWAyuBfYGnAXtk5vcBIuJc4Arg7cBbeuz/4cCTgNdl\n5qfbjrEWOD4iHpaZl43yBUkLkelaTQL7ozTTWNO1HQFeywXl87bl877AqlaAV253M8UI3QtnOcR6\n5fNNHctbv3tOojQAR1fUJPZHqVoTg5zdy+dLyuedgQsryl0MLIuIjbrtKDMvBr4HvCMiHh8Rm0TE\nk4B3Av+Wmb8YYb2lBc8vSzWF6VppeI0K8iJiO4rz5VZk5k/LxVsAayqKry6fN59lty8GrgJ+AtwM\nnAf8Eth/6ApLE8J0rSaB/VGaqTFBXnlF7WnAncCBbasG/tRGxBLgFIoLOl4H7Aa8HngicEr4L6I0\nED86ahL7o1StEUFeRGxIcY7dDsCemXld2+o1FKN5nbZoW9/NPsDzgZdn5qcz84eZ+SngFcALyvXd\n6jTrY2pqqu/XKC0EflmqKUzXSjA1NdVXPNLN2IO8iFhKMdq2K/CCzLyoo8hFFOfldXokcFVm3tpj\n948sny/oWP6T8nmnbhtm5qwPgzxNCtO1mgT2R02aqampvuKRbsY9GfIS4EsU06W8KDPPryh2OrBd\nROzWtt2mFKNwp89yiGvL5yd2LH9y+dx1jj1J3Tm6oiaxP0rVxj2S9zGKCyA+CNwWEU9pe2xXljkd\nOBc4KSL+KiL2LJcl8L72nUXE3RHxmbZF3wSuAb4YEa+PiD0i4g3AF4Cry/WSSn5ZqilM10rDG3eQ\ntxdFsHYEcE7H4yCALD7pewMrgI8D3wDuopgcuXMkbgltrykzb6G4q8W/Aoe0PZ8GPHWWVK+0aJiu\n1SSwP0ozjfWOF5n54D7LraEI+g6apdy9gtbyIo7XDlRBSZUcXVGT2B+lauMeyZPUIH5ZqilM10rD\nM8iTZLpWE8H+KM1kkCepNkdX1CT2R6maQZ6kdfyyVNPYJ6XBdb3wIiKuYIhbipU+nJnHDrkPSXPM\ndK2aZpC+ZX+UZup1de32wE3lYxDLgM0G3FaSJElDmG0KlQ9l5lGD7Dgi1g6ynaTx6Tc1ZgpN86Wf\nvmZ/lKp5Tp4k07VqHNO10vB6jeQ9ieHu7Trs9pIkSRpQ1yAvMy8YZsfDbi9p/pmuVdOYrpUGZ7pW\nkulaNY7pWml4PYO8iFgSEV+PiJMj4j49yt0nIr4WEV8dfRUlSZJU12wjefuVj1Mz885uhcp1pwIv\niYj9R1g/SfPIdK2axnStNLjZgry/BK4F+hmhO5niQouXDlspSfPLdK2axnStNLzZgrwnAiuzj09O\nZq4FzgSeMIqKSZIkaXCzBXkPoBjJ69cq4E8Gr46kcTJdq6YxXSsNbrYg705ggxr72wC4a/DqSJJk\n6lUahdmCvF8Dj6mxv0cD1w1eHUnj4Dl5mgT2R2mm2YK8HwK7R8TDZttRRDwUWA78YAT1ktRgpsfU\nJPZHqdpsQd4nKe6K8fWI6HquXURsA3wNWA/41OiqJ2k++WWppmiNytknpcH1unctmXl+RHwSeB1w\nYUR8CvgPpi/G+FPg2cBrgS2BT2Tm+XNYX0lzwHStJoH9UZqpZ5BXejPFCN3BwOHlo/VJav8X69Nl\nWUkTztEVNYn9Uao2671rM/OuzHwt8Ezgi8AVwB3l4wrgC8AzMvN1mXn3XFZW0tzyy1JNYbpWGl4/\nI3kAZOaPgB/NYV0kjYnpWk0C+6M006wjeZLUydEVNYn9UapmkCdpHb8s1RSma6XhdQ3yIuK2iPiH\nQXc87PaS5o/pWk0C+6M0U6+RvA2occ7eHGwvqaEcXVGT2B+larMFYS+KiB0G2K+fOGkB8stSTWG6\nVhrebEHeY8vHoPx0SguA6VpNAvujNFOvIG/HEex/zQj2IalhHF1Rk9gfpWpdg7zMvHIe6yGpAfyy\nVFOYrpWG5xQqkkzXaiLYH6WZDPIk1eboiprE/ihVM8iTtI5flmoK07XS8AzyJJmu1USwP0ozGeRJ\nqs3RFTWJ/VGqZpAnaR2/LNUUpmul4RnkSTJdq4lgf5Rmqn1v2YjYBtgP+DNg48w8qFy+NfBg4MLM\nvHWktZTUKI6uqEnsj1K1WkFeRBwMHAvct1yUwEHlzw8AzgNeC3xmVBWUNH/8slRTmK6Vhtd3ujYi\nngt8EvjUqnuIAAAgAElEQVQF8GLgeNruTZuZPwcuAl444jpKmmOmazUJ7I/STHVG8g4FfgMsz8yb\nIuJxFWX+B3jKSGomqbEcXVGT2B+lanUuvHgC8K3MvKlHmWuBBw5XJUnj4pelmsJ0rTS8OkHefYBb\nZimzGXDP4NWRNA6mazUJ7I/STHWCvKuAx89S5kkU5+xJmmCOrqhJ7I9StTpB3qnAbhHxl1UrI+JA\n4DHAv4yiYpLmn1+WagrTtdLw6gR576cYzftyRHwVeCpARLwpIr4GfBq4DPhovzuMiP0j4tSIuDoi\nbo2ISyPimIjYpKPc5hHxmYi4PiJuiYgVEbFLjeNsFxGfi4hfR8TtEXF5RBzT7/bSpDNdq0lgf5Rm\n6vvq2sxcHRHLgROBl7StOrZ8/gHwN5k523l77d5KcbHGYeXz44ApYI+IeFpmZhT/xp0BLAPeBNwI\nHA6sjIjHZuaqXgeIiB2AHwG/Av4O+C3FpM0PqVFPSW0cXVGT2B+larUmQ87Mq4DlEfEYipG8LYGb\ngHMz8z8HOP7emXlD2+9nR8RqikByObAS2Bd4GrBHZn4fICLOBa4A3g68ZZZjfAK4pty+dVHIDwao\nqzTx/LJUU5iulYbXd5AXESuBH2bmOzLzZ8DPhj14R4DXckH5vG35vC+wqhXgldvdHBFnUEy83DXI\ni4iHAM8DXtEW4EnqYLpWk8D+KM1U55y8JwPrzVVF2uxePl9SPu8MXFhR7mJgWURs1GNfTy+fby/P\n47s9IlZHxIkRscWI6istOo6uqEnsj1K1OkHeL4EHzVVFoLhAAjgKWJGZPy0XbwGsqSi+unzevMcu\nW6OBnwMuBfaiuHPHnwPfDf8ySDP4kVBTmK6VhlfnnLxPA0dFxPbluXkjVV5RexpwJ3Bg26phxt9b\nQezKzPy78uezIuIm4GRgT+A7Q+xfmgimazUJ7I/STHVG8r5FccHCDyPi7yLiyRGxfUQs63zUrURE\nbEhxBe0OwJ6ZeV3b6jUUo3mdtmhb303rnL8VHctbvz+mR51mfUxNTfU4tDS5HF1Rk9gfNammpqb6\nike6qTOS96u2nz/So1xS49y9iFgKnALsCjw3My/qKHIRxcUTnR4JXJWZt/bYfdW5fH3xP0ItRn5Z\nqilM10pFkNfPgFK3z0mdIO8LfZbrOzqKiCXAlyimS9k7M8+vKHY6cGBE7JaZZ5fbbQrsA5w0yyHO\nA35DcS7ex9qW71U+/6TfukqTzHStJoH9UZqpzmTIr5qD438M2B84GrgtIp7Stu6acqLj04FzgZMi\n4hCmJ0NO4H3tO4uIu4ETMvPgss73RMRhwAkRcTzwTeChwLspztM7cw5ekzTxHF1Rk9gfpWp1zsmb\nC3tRBGtHAOd0PA4CyOJfs70pzqP7OPAN4C6KyY0773axhI7XlJlfAA4AnkERMB4BfJFiJFBSG78s\n1RSma6Xh1brjxahl5oP7LLeGIug7aJZylUFrZp7E7KldSZKkiVHnjhefp8/z7TLz1QPXSNK8q3su\nU2t0xXOg1AT2R6lanZG8V9Yoa5AnSRqYAZs0vDpB3o5dlm8GPAF4J8W5dIcOWylJ4+H5T2oa+6Q0\nuDpX117ZY/V/R8R3gf8B/h34zJD1kjSPTNdqIbM/StVGdnVtZl5DcVeMN49qn5KkxcmATRreqKdQ\n+S3w8BHvU9I8MTWmprFPSoMbWZAXEesBewA3jWqfkuaH6VotZPZHqVqdKVR267GPZcCBwOPwfDxJ\n0pAM2KTh1bm69qw+ypwNHDJYVSSNm6kxNY19UhpcnSDvqC7L1wJrgB9n5vnDV0nSfDNdq4XM/ihV\nqzOFytQc1kOSpHUM2KTh9X3hRUTsFhHLZimzrMe5e5IaztSYmsY+KQ2uztW1ZwGvmqXMAcDKQSsj\naTxM12ohsz9K1UY9T57/ckmShmbAJg1v1EHeMuAPI96npHliakxNY5+UBtfzwouIOBJIpkfolnf5\nwK0HbA/8NfDDUVZQ0twzXauFzP4oVZvt6tojO35fXj66WQUcNkR9JEkyYJNGYLYg71ltP58JnFg+\nOt0D3ABcmplrR1Q3SfPM1Jiaxj4pDa5nkJeZZ7V+jogvAKe2L5M0GUzXaiGzP0rV6kyG/Ko5rIck\nSZJGaNRX10pawEyNqSlao3L2SWlwtYK8iNg2Ij4eEb+KiNsi4p6Ox9qIuGeuKitpbgya5jI9piax\nP0oz9Z2ujYjtgJ8A2wAXAxsAVwF3AjtSTKPy38BNo6+mpCZxdEVNYn+UqtUZyXsn8CfA8zPz0eWy\nz2fmI4AHA98FNgT2G20VJc0XvyzVFKZrpeHVCfL2BL6bmSs6V2TmtcBLgI2Ad42obpLmielaTQL7\nozRTnSDvAcCFbb/fQzFyB0Bm3gKsAPYdTdUkNZWjK2oS+6NUrU6Q9wfgPm2/3whs11HmJopz9iQt\nQH5ZqilM10rDqxPkXQU8qO33nwHPioiNASJiCfBc4NrRVU/SfDBdq0lgf5RmqhPk/TtFULe0/P0E\nYFvgnIh4P3AOsAvw1ZHWUFLjOLqiJrE/StX6nkIF+BxFinZr4LrMPCkiHg+8GXhUWeZk4OjRVlHS\nfPHLUk1hulYaXp3bmv0v8E8dy/5vRLyHYp68KzLztyOun6R5YLpWk8D+KM1UZzLkVwK/yczvti/P\nzN8Bvxt1xSQ1l6MrahL7o1Stzjl5nwX2mquKSBo/vyzVFKZrpeHVCfJ+W7O8pAXCdK0mgf1RmqlO\n0PZtYI9yqhRJi5ijK2oS+6NUrU7AdgRwP+BzEbHVHNVH0hj5ZammMF0rDa/OFConAzcDBwB/FRFX\nAr8B7jU+npnPGkntJM0L07WaBPZHaaY6Qd7ubT9vADyifEhaZBxdUZPYH6VqdebJ81w8acL5Zamm\nMF0rDc/ATZLpWk0E+6M0k0GepNocXVGT2B+larWCvIhYLyLeHBE/joibI+KetnWPi4iPR8TDR19N\nSfPBL0s1helaaXh9B3kRcR9gBfBhinvV/gFo//RdCbwaePkI6ydJkqQB1BnJOwRYDrwLeADw6faV\nmbkG+AHwvFFVTtL8qHsuU2t0xXOg1AT2R6lanSDvZcA5mfmuzLynS5krgGXDV0uStJgZsEnDqxPk\nPRg4d5Yyq4EtB6+OpHHy/Cc1jX1SGlydIO8OYLNZyjwIuHHw6kgaB9O1Wsjsj1K1OkHefwHPi4gN\nqlZGxP2BPYHz+91hROwfEadGxNURcWtEXBoRx0TEJh3lNo+Iz0TE9RFxS0SsiIhdatS9tZ+/joi1\nEXFN3W0lSfPHgE0aXp0g71MUI3VfiohN21dExObACcAWwCdq7POtwF3AYcBewPHAG4AVUf5rVj6f\nQXFBx5uA/YClwMqI2K7fA0XEZhRXBlfeb1eSqTE1j31SGlyd25p9JSKeC7wK2IcyLRsRFwC7APcB\nPp6Z/1rj+Htn5g1tv58dEauBEymu5F0J7As8DdgjM79fHvNcios83g68pc9jvY9iNPI3wHNq1FGa\neKZrtZDZH6VqtSZDzsxXU8yFdzGwdbl4V+Ay4KDMfFPN/d1QsfiC8nnb8nlfYFUrwCu3u5lidO+F\n/RwnIp5OcXXwG5k5t58kqYEM2KTh9T2S15KZJwAnRMRGwObATZl5ywjrtHv5fEn5vDNwYUW5i4ED\nImKjzLy1284iYilFqvl9mXm5Q/9Sd34+1DT2SWlwtYO8ljKw6hpcDaI8x+4oYEVm/rRcvAVweUXx\n1eXz5rPU41CKc/jeM6p6SpPGdK0WMvujVK1WuhYgIu4XEQdExAcj4rPl8ys6r4gdYL+bAKcBdwIH\ntq0a+FMbEQ8F/gF4U2beWWefETHrY2pqatCqSZJ6MGCTYGpqqq94pJtaI3kR8ZcUV89WzZf3kYh4\nXWZ+vd5LgIjYkOIcux2A3TPzurbVayhG8zpt0ba+m2OBM4Efl1fXQnGByJJyypc7MvP2qg39A6PF\nyNSYmsY+qcVsamqqrwGlbp+TvoO88sraLwNrKa5+/T7FlaoPoLgS9mXAlyPixsxcUWO/S4FTKC7g\neG5mXtRR5CKq74f7SOCqXufjAX8GbE91ILiGYkqVv++3rtKkMl2rhcz+KFWrM5L3TopU6jMz8z87\n1p0QEccBPyjL9RXkRcQS4EsUQeLemVk1kfLpwIERsVtmnl1utynFNC4nzXKIvwbaJ28Oijn5Hg/s\nD6zqp56SpPllwCYNr06Q9zjgqxUBHgCZeUFEfJUieOrXx8ryRwO3RcRT2tZdk5mrKIK8c4GTIuIQ\nivn5Dqc4r+597TuLiLuBEzLz4LJOP+48YEQcSJGmPbtGPaVFwdSYmsY+KQ2uzoUXdwLXzVLm12W5\nfu1FEawdAZzT8TgIIIt/5/amGB38OPANirtk7FEGge2WMPtrSrzjhTSD6VotZPZHqVqdkbyzgafP\nUuZpZbm+ZOaD+yy3hiLoO2iWcrMGrZl54GxlJEnjZcAmDa/OSN5hwKMj4r0RsXH7iojYJCLeBzyK\nYl46SQuQqTE1jX1SGlydkbxDgf8BDgFeExE/BX4L/AnFlbGbUYziHdr5oSxvhyapoUzXaiGzP0rV\n6gR5r2z7eTPgWRVldisfnQzyJEl9M2CThlcnyNtxzmohqRFMjalp7JPS4PoO8jLzyjmsh6QxMl2r\nhcz+KFWrfe9aSZLmmgGbNLxa964FiIj1gG2BPwWWVpVxomFpYTI1pqaxT0qDq3Pv2qC4svZtwFY9\niiaw3pD1kjSPTNdqIbM/StXqjOQdSXFf2huAEynu+3p3RTk/ZZKkoRiwScOrE+QdBFwB7JqZN81R\nfSSNkakxNY19UhpcnQsvtgROM8CTJo/pWi1k9kepWp0g71fA5nNVEUmSWgzYpOHVCfKOB/aJiAfO\nVWUkjZepMTWNfVIaXJ3JkD8eEY8AfhgR/wj8J1CZus3Mq0dUP0nzwHStFjL7o1St7jx5PwNeBXyu\nRxmnUJEkDcWATRpenXnyXgN8ErgLOAu4DqdQkSaKqTE1jX1SGlydkby3Ar8DnpqZV8xRfSSNgela\nLWT2R6lanQsvtge+boAnSZprBmzS8OoEedfR5V61kiaDqTE1jX1SGlydIO9E4M8jYtO5qoyk8TBd\nq4XM/ihVqxPkHQOcD6yIiD0i4n5zVCdJ0iJnwCYNr86FF3e2/fwfQFYMoweQmekUKtICZGpMTWOf\nlAZXJ8g7u89y/vslLTCma7WQ2R+lanXueLF8DushSdI6BmzS8OqckydpwpkaU9PYJ6XB1b2tGQAR\nsTHwCGDjzPzBaKskqen84lWT2B+larVG8iLiQRHxDeBG4AKK25u11j0zIi6OiOUjraGkOTdoasyU\nmubKIH3L/ijN1HeQFxEPBM4D9gW+BZxLcTVty4+BPwH+apQVlCRJUn11RvKOpAjinpeZLwZWtK/M\nzDuBHwBPH131JM2nftNepsc0X/rpa/ZHqVqdIO8FwOmZeWaPMlcD2w5XJUnzzXStmsZ0rTS8OkHe\nnwD/O0uZu4BNBq+OJEmSRqFOkLcGeNAsZR4G/Gbw6kgaJ9O1ahrTtdLg6gR5PwT2LS/AuJeIeBiw\nF7ByFBWTNH9M16ppTNdKw+sZ5EXEKyPi0eWv7wc2BL4fEc8vfyYiNomIF1BccZvAP89hfSVJktSH\n2SZD/jwwBfxPZv44Il4LfAL417YyN1FMpXIX8OrMvHAuKipp7pmuVdOYrpUGV+uOF5n5uYj4IfAG\n4KnAlhRB3rnAcZn5i9FXUdJcM12rSWB/lGaqfVuzzPxf4P/OQV0kSQIM2KRRqHVbM0mTzbSXmsY+\nKQ2un5G8zSJiWZ2dZubVA9ZH0hjUHTVpffE62qImsD9K1foJ8v4P8JY+9xcUV9iuN3CNJEmLngGb\nNLx+grybyke//GRKC5SpMTWNfVIaXD9B3ocz811zXhNJY2O6VguZ/VGq1s+FF35qJEnzyoBNGp5X\n10pax9SYmsY+KQ3OIE+S6VotaPZHqVo/QZ7/RkmS5pUBmzS8nhdeZKYjfdIiYmpMTWOflAY39iAu\nIvaPiFMj4uqIuDUiLo2IYyJik45ym0fEZyLi+oi4JSJWRMQufez/4RHx0Yi4OCL+EBHXRcRpEfHo\nuXtV0sJiulYLmf1Rqjb2IA94K3AXcBiwF3A88AZgRZSf3PL5DOB5wJuA/YClwMqI2G6W/T8P2AP4\nHLAP8LfA1sB5EbHryF+NJGloBmzS8PqZJ2+u7Z2ZN7T9fnZErAZOBJYDK4F9gacBe2Tm9wEi4lzg\nCuDt9L4jx1cy87j2BRFxJnBlud0rR/MypIXP1Jiaxj4pDW7sI3kdAV7LBeXztuXzvsCqVoBXbncz\nxejeC+vuv9z2srb9S4ua6VotZPZHqdrYg7wudi+fLymfdwYurCh3MbAsIjaqs/OI2ALYpW3/kqQG\nMWCThte4IK88x+4oYEVm/rRcvAWwpqL46vJ585qH+SjFnTw+PFAlpQllakxNY5+UBteoIK+8ovY0\n4E7gwLZVI/uXLiIOB14KvCkzL+9RbtbH1NTUqKoljZXpWi1k9kdNqqmpqb7ikW6acOEFABGxIcU5\ndjsAu2fmdW2r11CM5nXaom19P8d4PXA0cERmntCrrH8sJGl8/BssFUFePwNK3QK9RozkRcRS4BRg\nV+AFmXlRR5GLKM7L6/RI4KrMvLWPY7wC+Bjwgcx8z5BVliaSqTE1jX1SGtzYg7yIWAJ8iWK6lBdl\n5vkVxU4HtouI3dq225Ri3rvT+zjGiynmyft0Zr59FPWWJonpWi1k9kepWhPStR8D9qdIo94WEU9p\nW3dNZq6iCOTOBU6KiEOAG4HDKc7Ve1/7ziLibuCEzDy4/H034CvAz4ATO/Z/R2b+19y8LEnSoAzY\npOE1IcjbiyJYO6J8tJsCjsrMjIi9gQ8AHwfuC5xDMTnyqo5tljBzhHIP4D7A44AfdZS9Ethx+Jcg\nTQZTY2oa+6Q0uLEHeZn54D7LrQEOKh+9yi3p+P1dwLsGrqC0CJiu1UJmf5Sqjf2cPEmSOhmwScMz\nyJO0jqkxNY19UhqcQZ4k07Va0OyPUjWDPElS4xiwScMzyJO0jqkxNY19UhqcQZ4k07Va0OyPUjWD\nPElS4xiwScMzyJO0jqkxNY19UhqcQZ6k2vziVZPYH6VqBnmSBk6NmVLTXBmkb9kfpZkM8iRJkiaQ\nQZ6kdfpNe5ke03zpp6/ZH6VqBnmSTNeqcUzXSsMzyJMkSZpABnmS1jFdq6YxXSsNziBPkulaNY7p\nWml4BnmSJEkTyCBP0jqma9U0pmulwRnkSTJdq8YxXSsNzyBPkiRpAhnkSVrHdK2axnStNDiDPEmm\na9U4pmul4RnkSZIkTSCDPEnrmK5V05iulQZnkCfJdK0ax3StNDyDPEmSpAlkkCdpHdO1ahrTtdLg\nDPIkma5V45iulYZnkCdJkjSBDPIkrWO6Vk1julYanEGeJNO1ahzTtdLwDPIkSZImkEGepHVM16pp\nTNdKgzPIk2S6Vo1julYankGeJEnSBDLIk7SO6Vo1jelaaXAGeZJM16pxTNdKwzPIkyRJmkAGeZLW\nMV2rpjFdKw3OIE+S6Vo1julaaXgGeZIkSRPIIE/SOqZr1TSma6XBGeRJMl2rxjFdKw3PIE+SJGkC\nGeRJWsd0rZrGdK00OIM8SaZr1Tima6XhjTXIi4j9I+LUiLg6Im6NiEsj4piI2KSj3OYR8ZmIuD4i\nbomIFRGxS5/HiIg4PCKujIjbIuK/I+Iv5uYVSZIkNcO4R/LeCtwFHAbsBRwPvAFYEeX4e/l8BvA8\n4E3AfsBSYGVEbNfHMd4NHAkcWx7jPODrEfH80b4UaeEzXaumMV0rDW79MR9/78y8oe33syNiNXAi\nsBxYCewLPA3YIzO/DxAR5wJXAG8H3tJt5xGxDfA24JjM/GC5+PsR8VDgn4Bvj/blSAuT6Vo1jela\naXhjHcnrCPBaLiifty2f9wVWtQK8crubKUb3XjjLIfakGPU7qWP5ScCjImL72pWWJElaAMadrq2y\ne/l8Sfm8M3BhRbmLgWURsVGPfe0M3JGZv6rYFuCRA9dSmkCma9U0pmulwTUqyCvPsTsKWJGZPy0X\nbwGsqSi+unzevMcuZ9t2i0HqKUmS1HTjPidvnfKK2tOAO4ED21YNe5LFQP/irb9+Y5pGmnNr164d\naLu3v/3tHHbYYSOujTTY+XWXXXaZf7ulNo34NETEhhTn2O0A7J6Z17WtXkP1iNsWbeu7WQNs1mPb\n1RXrALjnnnt67FaaPEuXLmX58uV9ld1tt9047rjjuP322/2saM485znP6SsVu2zZMnbaaScuvfRS\n+6PUZuxBXkQsBU4BdgWem5kXdRS5iGL6lE6PBK7KzFt77P4iYIOIeEjHeXmtc/EurtgGgLvuumvW\nukuTJCJYb731+iq7zz77cMstt3g1o+bUeuut11eQt8EGG3DxxRcb4GnRWrp0aeXysQZ5EbEE+BLF\ndCl7Z+b5FcVOBw6MiN0y8+xyu02Bfbj3VbOdvk0xD9/LKM71a3k58PPMvKrbhg75S731GxBK8yEi\n/LstdRj3J+JjwP7A0cBtEfGUtnXXZOYqiiDvXOCkiDgEuBE4nOJcvfe17ywi7gZOyMyDATLz+oj4\nIHB4RPwB+C/gr4A9KIJESZKkiRTjTLdExBXAMqovjpjKzKPKcpsDHwBeBNwXOAf4+8z8ecf+1lIE\nea9uW7aEIih8DfAA4FLgqMz8Ro96pWkoSZK0EEQEmXmvWGqsQV5TGeRJkqSFoluQ16h58iRJkjQa\nBnmSJEkTyCBPkiRpAhnkSZIkTSCDPEmSpAlkkCdJkjSBDPIkSZImkEGeJEnSBDLIkyRJmkAGeZIk\nSRPIIE+SJGkCGeRJkiRNIIM8SZKkCWSQJ0mSNIEM8iRJkiaQQZ4kSdIEMsiTJEmaQAZ5kiRJE8gg\nT5IkaQIZ5EmSJE0ggzxJkqQJZJAnSZI0gQzyJEmSJpBBniRJ0gQyyJMkSZpABnmSJEkTyCBPkiRp\nAhnkSZIkTSCDPEmSpAlkkCdJkjSBDPIkSZImkEGeJEnSBDLIkyRJmkAGeZIkSRPIIE+SJGkCGeRJ\nkiRNIIM8SZKkCWSQJ0mSNIEM8iRJkiaQQZ4kSdIEMsiTJEmaQAZ5kiRJE8ggT5IkaQIZ5EmSJE0g\ngzwNbWpqatxVWHBss3psr3psr/pss3psr/rG0WaRmfN+0KaLiLRd+hcR2F712Gb12F712F712Wb1\n2F71zWWblfuOzuWO5EmSJE2gsQd5EfGnEfHRiDg3Im6NiLURsayi3GMj4jsR8YeIuCkiTouIh/R5\njK0j4riIuLw8xuXlMbca/SuSJEkav7EHecBDgZcANwBnVxWIiIcBPwDuB/wN8GpgB+DsiNi6184j\nIoAzgL8E3gvsBbwf+OtyuSRJ0sRZf9wVAL6fmQ8AiIiDgedVlDkUuAt4fmbeXJY9D/gl8LZyfTcP\nB54EvC4zP10uOzsi1gLHR8TDMvOy0bwUSZKkZhj7SF6fVzg8BTi3FeCV260CLgJePMu265XPN3Us\nb/0+9jaQJEkatYUS4NwN3Fmx/A5gx4i4T7cNM/Ni4HvAOyLi8RGxSUQ8CXgn8G+Z+Ys5qbEkSdIY\nLZQg7xfAEyJiXXo5Iu4H7AwEsPks278YuAr4CXAz0Er17j8ntZUkSRqzhRLkHQtsB3wiIraNiO2B\nzwMbl+vXdtswIpYApwCPA14H7Aa8HngicEp5YYYkSdJEacKFF7PKzB9FxBuB91BcWQuwAjgReDmw\nusfm+wDPB56dmSvLZT+MiMsp0rj7AKd3bmTsV4/tVZ9tVo/tVY/tVZ9tVo/tVd98t9mCCPIAMvP4\niPgMxZQrN2fmqoj4NnBeZt7TY9NHls8XdCz/Sfm8Ex1BXtWs0ZIkSQvJQknXApCZd2XmJWWA9yjg\n2cDxs2x2bfn8xI7lTy6fV42yjpIkSU3QiHvXRkTrAohnU5w397fA74HfZebZEbFduewciitqnwAc\nBnw7M1/asa+7gRMy8+Dy900oplpZH/hHios4dgKOBG4HHpmZt87tK5QkSZpfTQny2i+cSIorZgHO\nysxnRcQ2wJeAx1Lc9eKXwGeBj2Tm2op9nZCZr25bti0wRRFEPhD4DcU5fVOZ+es5eVGSJElj1Ih0\nbWYuaXus1/bzs8r1v8vM52bm1pl538zcJTM/1Bngte3r1R3LrsvM12bmQzJzo8zcMTNf1x7gRcSD\nIuKUiLixvDfuv0TEg+b+1TdbROwfEadGxNXlfX8vjYhjyhHS9nKbR8RnIuL6iLglIlZExC7jqneT\nlPdcXhsR/9ix3DZrExEviIiz2+5P/ZOI2KNtve1Viohnlq//dxFxc0T8Z0Qc2FFmUbZXjfuh99U+\nEXHfiHh/RPy63N85EfHM+Xk186OfNouI50TEl9vuAf/LiPh41a1FJ73N+u1jHdt8oiz3xYp1c9Ze\njQjyxi0iNgLOpLgF2gHAK4CHASvLdYvZWyluKXcYxX1/jwfeAKyI8jKh8vkMilvSvQnYD1hK0X7b\njaPSTRERLwUeXf6abcttszYR8TrgVIoLol5EcT/rrwEblettr1JEPI4iE7EEOIhiHtCfAJ+NiNeX\nZRZze/VzP/Q67fNZ4GDg/wF/Dvwa+G5EPGZOaj8es7YZ8FpgS+DdwJ4Us13sC5wXERt3lJ30Nuun\nvdaJiKcDL6OYp7cqfTp37ZWZi/4BvIXirho7ti3bgSK4+b/jrt+Y22bLimWvoJibcI/y9xeWv+/e\nVmZTig/AR8b9GsbYdpuXH9a/KtvnqLZ1ttn0694BuA14c48yttf0634PxfnEG3UsPwc4Z7G3F+Vp\nSOXPB5ftsKyjTF/tAzymLPfKtmXrAZcCp437tc5zm21Vsd0zy7IHLqY266e92tYvBS4EDgWuAL7Q\nsX5O28uRvMK+FPfGvby1IDOvBH5E8cdg0crMGyoWt6aj2bZ83hdYlZnfb9vuZor/lBdz+70X+Hlm\nfgaUR4QAAA8nSURBVLVinW027dUU/2R9okcZ22vaehT/gN7Wsfxmps9nXrTtleW35Cz6bZ99Kdr6\nq23l7gFOBvaMiKUjqfSY9dNmmfn7isWd3wWwCNqszz7WcgjF5/Kfmf58tpvT9jLIK+xMEWl3upjp\nefY0bffy+ZLyuVf7LVuMKe+IeAbFiOcbuxSxzaY9g+Kq97+JiF9FxF0RcVlE/G1bGdtr2meBe4Bj\nI+KBEbFZRLwGeBbwobKM7dVbv+2zM3B5Zt5eUe4+FGm7xazzuwBss3Ui4qHAEcDfZubdXYrNaXsZ\n5BU2B9ZULF/N7PfFXVTK81WOAlZk5k/LxVvQvf1gkbVhRNwH+CTw/sy8rEsx22zathTnwL4POAZ4\nLsU5Z8dFxJvLMrZXKTN/QXFO1Eso5vlcDRwHvC4zv1YWs71667d9Ziu3xYjrtWBEcf/4D1MEI6e2\nrbLNph0P/EvbiHHVCOCctteCueOFxi+KK2pPA+4E2q/kG/88PM3ydmAD4OgeZWyzaUsopkZ6ZWa2\nvizOiogdgMMp7l2tUnkF6LcoUmUfpUjbvgj4ZETckZlfHmf9Fgg/f0OIiPWBr1BMSfb0rJjpYrGL\niJcDjwceMc56GOQV1lD9n+0W9L4v7qLx/9s792irqioOfz/xgS9E8RmUYKVIlmaWjxwE+cBMMdMY\nVg5By54jaxSlpSGGQ83UnqNSzFAzR6bmo0yj4PogzCdmgISFgi80BFFQAu7sj7l2Y7PvPvecC+dy\nYN/5jbHGvnetudeaa+59zp5nPeaWtCW+XmUgvlj5uVzxYsp/beyQK+8RpG30Z+O7HrdMdsvoLWk7\n4DXCZnkWAW/FR+/yTAaOkrQrYa88E4AlwLG5KaCpkvoBP5R0PWGvejRqn8VAWWiMTK7HPR8kbYK/\nN/6DwIfNrDjt3eNtlgZELsNnJ1ZK6puKegGbp+fAsvT57VZ7xXStMxMoix81BB+K7tGkhZ83AvsD\nR5vZzILITHxdQZEhwNPWs94osgc+ivcr/MOZJYCx+Ad6H8JmeWZSviC5KBP2coYAfy9Z4/MgHuJi\nZ8Je9WjUPjOBQZJ6l8j9Fw/M39P4OTAKOMnMppaUh81gx5QuYM3nwADcdouBo5Nst9ornDznNuAg\nSYOyjDRVdEgq67GkX23XAcOAj5jZAyVitwH9JQ3NndcHOJaeZ79HcVvlUxbQ99r0/5OEzfLcnI5H\nFfKPAhaY2QuEvfI8A+xbsuvuQHzqdhG+rCLsVZtG76fb8BAYo3Jym+Jhke4ys5XrR90NA0mX4rMU\nY8ys1n0UNvPQWcPp+BxYiM9QDMOjd0B326vV8WY2hIQHXJ0L/B3fzjwSeAx/GG/Vav1abJuf4TF8\nJgAHFVL/JKN0w85PN+YIoA1//3D/VvdhQ0h0jJMXNlvTPn9Jff8sHqB2YrLZKWGvDrbKYrzdmb6r\njsQ3XrQDl4S9DODElLLvr8+l/4d21T742rOXcefmMHxWYzmwX6v7uZ5tdmbKvxL/QZF/FuxRqKvy\nNqtnrxrnPEUhTl5326vlhtpQEvDmZNhX8HhTN1MjuGFPSnjwxtXpJi6mcTm57fHQDouAZfivlXe2\nWv8NJVFw8sJmHeyzbXJUXgBWADPw6aCwV7m9jsDf0vNi+r56JD1kNgl7/f/zlqX899eUrtoH6I3H\nOHseHymd3tmDfGNN9WwGTO3kWXBVT7NZI/dYyTkdgiF3t72UGgiCIAiCIAgqRKzJC4IgCIIgqCDh\n5AVBEARBEFSQcPKCIAiCIAgqSDh5QRAEQRAEFSScvCAIgiAIggoSTl4QBEEQBEEFCScvCIIgCIKg\ngoSTFwRBEARBUEHCyQuCoFuQNExSu6Rz17GeMame0c3SrZO2Bqa2ftndbVUNSZOS7bJ0Zqt1ahRJ\nYwu6x/UPKkE4eUFQEXIPqNWS9uhEbmpOttsdJ6BZr9VpqJ6cU5hPr0paIGmypPMkva0ZbZW03RTH\ndiPnB8B44N4W69EVpuE6/zD9H6+CCirBpq1WIAiCprIK/1x/Cji7WCjp7cAHcnJVfpjNAG5Jf28J\n7Iy/TP3bwNmSfgyMNbPVuXOeAQbj77BeF6ps13r8wMzmt1qJrmBm04HpknYHvtxqfYKgWYSTFwTV\nYiH+kutTJY0rODAAn07H24Hj16tm658ZZvadYqakYcAk/GHeG/h8VmZmq4B/NqFtNaGOYP0T1y2o\nFDFdGwTVwoCJwK7AMfkCSZsBY/CpqVm1KpD0dknXSHpW0op0vLrWFKekXST9QtJCScslPSrplM6U\nlLSDpAslzU7nLJH0Z0lHdLG/XcbM2oARwH+Bz0jaL6dX6Zq81MdLJM2R9JqkxZKekPRLSYOSzCRg\nSjrl3MJ08dAk00fS1yVNkfRMsu+Lkm6VdFCZvun8qZL6SbpC0vOS3pD0D0ljavVT0pGSbk/1vyFp\nvqRbJB1WIjtC0h2S/pNkn5R0saTtumLbzpDUlvqyqaRxkv4l6fVkx9Nzcl+Q9Hi6LxZIGi9Jhbr+\nf50kvVXSjZIWSVoq6U+S9klyO0m6MtnsdUkPJic/CHoEMZIXBNXjeuAyfNTu1lz+SGAn4OvAnmUn\nSnov8Gdgm3TuLGBv4GTgOEmHm9lDOfkdgb8Cg/A1WPcBbwJ+Dkyu0cbuQBuwO3APcEdq7xjgTkmf\nNbMr16LfDWNmcyTdgPfrE/jU7hoiOX23wh3jPYA/4XYRMBC36W+BecDv0nmj8f615ep7Kh2HAOcD\nd+OjqYtxO4wEPiTpWDO7q0TlvkmHFcANwBbAKOAqSe1mdk1eWNJ5+LT0q/iU9QKgP3AI8EngLznZ\nc4FzgUVJpxeBfYGxwNGSDjazV0t0Wlt+A7wP+AOwEvgYcLmk1andTyY9JgPHAeOA5cDFJXUNBO7H\n79Or8PvweKBN0qH4vbUY/0z0A04C/ihpTzNb0MQ+BcGGiZlFihSpAgloB+anvyfiD9D+ufI78Qde\nb9zRaAdOyZULmA2sBj5eqHtUkp8NKJd/Rcq/tCD/HnykrB0YVyhrw9cEjirkbwc8ij/Qd87ljynq\nWscOmfxVdeROS3JTc3kDi+cCx5b1MZVtCmyT+39YWZ9z5X2AHUry+wPPArNqXNf2ZOu87fdO13hm\nQf7IJP8ksFtZW7m/hyfZ+4A+BbnRqeyyBu0+Kcm/pUZ5Wyr/W74t3DFbga+D/Fde53RPvIQ7nr1K\nrlM78M1CO+ek/CXATwtlJ3fWp7LrHynSxpxiujYIqslEoBfuyGSjZ0cA15nZGzXOOQTYC5huZtfn\nC8zsBtwR2As4NNW5GT7qshTfmZiXfxi4rtiApH2BocBNqc78Oa+kenoDJzTc07XnuXTcqUH5DnYz\ns1Vm9lqjDZrZUjN7uST/WeAmYLCkASWnLgO+amaWO2c2Poo6OI02ZnwpHb9mZs/XaCvjjHQ83cyW\nFuSuBh7Dr3EzOSvflpnNw0cptwUm5HVO98TvgR3xEeIi84CLCnlXp2MvfNQ6z6/xHxj7rksHgmBj\nIaZrg6CCmNkDkh4HTpN0Pj51K9z5q8X+6TilRvlU3MHbD5+aHYzvWn3Qyqfz7sZHg/IcnI59JY0v\nOSdzuPbuRM9mka3zqrcTtg0fZTtL0v7AH3GHd4aZtXe5Uen9+KaPg/H+bl4Q6Y/v8s0zt4YzuQDv\nx/b4CCj4DuJ2fOS2Hgfjo4GjiuveEpsDO0na3swWN1BfPQx4qCQ/c7gfLinLnNIBeH/zzMg7vonM\nSfynmS1bo3GzdkkvprqCoPKEkxcE1WUi8CPgQ8CpwENm9lgn8tki+w6jP4X8vgX5hTXkXyjJ65eO\nR6RUhgFb1yhrJtnI0EudCZnZq2lTxHn42rkRqeg/kn4KnG++K7cuko4HbsQdssn49OQy3Ckbjoe3\n2aLk1CU1qsza7ZXL6wssNrMVDajUL53bWVw/w9dMNsPJo8YPgqwfZaFrsrLNSso6yJvZquSv1gqD\ns6pGXUFQOcLJC4Lqci3wXeBy3KEZX0c+eyjuWqN8t4JcdtylhnxZPdk5Z5jZT+ro090MT8e/1RNM\nU5yfBpA0BPgg8EV8U8Am6dgIE/Bp3wPMbE6+QFJ/3MlbV5YA20vq3cnUfMYrAGa2YxPaDYJgAyPW\n5AVBRUnrmW7Ep/9ew3cYdsYj6Ti8RvnwgtwTwOvAfpL6lMgPK8mbno5D6+jSrUgajO/qbMfXaTWM\nmc1KDmo2EnlcrjiLS9iLct6Gb64oOnibkNY6NoHp+Hf7UQ3K7pAc1yAIKkY4eUFQbc4BPgKMKK5P\nKmJm04A5wKGS1tj4IOlE3AmZY2b3JfmVwK/wHaPjC/IHULJgP23IuBf4qKRTy/SQ9E5JjW6G6DKS\nPoCvV9sM+JmZPV5HfoikstHKbKRyeS5vUTruXqO6ecCekrJRUdJauPH4OsRmvCnjx+l4qaQOmxUK\ned9Px4l5nXKyW0s6sAk6BUHQAmK6NggqjHkssK7EAxuNrxX7jaRbcadvL9xRXAoUgxx/CzgM+Epy\n7Kbh07qj8DhoI0va+AS+ueMXks4AHsCnGAcA7wLegW8e6HStXAO8O7e5Ywt8WvlA3JlaDVwKfKOB\neo4Evifpr8BcPJzHAHwEbzXwvZzsE/hGgZMkrQTm447bNeav+vo+HkPwUUk345se3p90uh0P17JO\nmNnktNnmHGC2pFvwjRy74I76dHyNJmY2RdJZwIXAXEl34DH9tsEd1aG4U370uuqViDdKBMF6JJy8\nIOiZGCWjRmlX7ntxB+Fw3Ol4CQ+HMsHM5hbkF6Xdohck2QNwR+dzwNOUOHlm9qyk9+ChPk7Anb5e\n+MaOWfhL4v9RT9c6fQN3GLNQGcuBl5NuNwDXmtm/G6zvTuDNuMMzEh+5fA64C4+3dn+ub+1pc8VF\n+HTwtqnoHjyG4RWSVgBfwR3m5bgTNRo4kcJbShrsa9l1HCdpOh4i5Rh8I8tCfGfr1QXZiyVNS7KH\n4s7rEtxZvZwuTmd3Vdd1LFtbPYKgR6COu8+DIAiCoGvIX+t2CjDIzJ5usTprhaSBwL+BSWZ2Wmu1\nCYJ1J9bkBUEQBM0gGzGYl94re2ZLtekCksZKascdvCCoDDFdGwRBEDSDW/CNJRn3tkqRtWAaa24e\nKr7LOAg2SmK6NgiCIAiCoILEdG0QBEEQBEEFCScvCIIgCIKggoSTFwRBEARBUEHCyQuCIAiCIKgg\n4eQFQRAEQRBUkP8BbGL38w9MrJsAAAAASUVORK5CYII=\n", | |
| "text/plain": [ | |
| "<matplotlib.figure.Figure at 0x10dc16a90>" | |
| ] | |
| }, | |
| "metadata": {}, | |
| "output_type": "display_data" | |
| } | |
| ], | |
| "source": [ | |
| "# Setup figure and axes\n", | |
| "# Generally plots is ~1.33x width to height (10,7.5 or 12,9)\n", | |
| "fig = plt.figure(figsize=(10,7.5))\n", | |
| "ax1 = plt.subplot(111)\n", | |
| "\n", | |
| "# Set labels and tick sizes\n", | |
| "ax1.set_xlabel(r'Model Distance [mm]', fontsize=20)\n", | |
| "ax1.set_ylabel(r'Temperature [C]', fontsize=20)\n", | |
| "\n", | |
| "# Plotting\n", | |
| "ax1.plot(X*1000., phi.value, color='k')\n", | |
| "\n", | |
| "plt.title('Temperature Initial Condition', fontsize=24)\n", | |
| "\n", | |
| "# Set limits\n", | |
| "ax1.set_xlim(0,np.max(X.value*1000.))\n", | |
| "ax1.set_ylim(19.8,21.2)" | |
| ] | |
| }, | |
| { | |
| "cell_type": "markdown", | |
| "metadata": {}, | |
| "source": [ | |
| "## State and Solve Equations" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 8, | |
| "metadata": { | |
| "collapsed": false | |
| }, | |
| "outputs": [], | |
| "source": [ | |
| "# Setup the Crank-Nicholson Scheme\n", | |
| "eqX = TransientTerm() == ExplicitDiffusionTerm(coeff=D)\n", | |
| "eqI = TransientTerm() == DiffusionTerm(coeff=D)\n", | |
| "eqCN = eqI + eqX" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 9, | |
| "metadata": { | |
| "collapsed": true | |
| }, | |
| "outputs": [], | |
| "source": [ | |
| "timeStepDuration = 0.01 # How long of a timestep we take\n", | |
| "steps = 100 # Number of timesteps" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 10, | |
| "metadata": { | |
| "collapsed": false | |
| }, | |
| "outputs": [], | |
| "source": [ | |
| "for step in range(steps):\n", | |
| " eqCN.solve(var=phi, dt=timeStepDuration)" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 11, | |
| "metadata": { | |
| "collapsed": false | |
| }, | |
| "outputs": [ | |
| { | |
| "data": { | |
| "text/plain": [ | |
| "(0, 145.94999999999999)" | |
| ] | |
| }, | |
| "execution_count": 11, | |
| "metadata": {}, | |
| "output_type": "execute_result" | |
| }, | |
| { | |
| "data": { | |
| "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnkAAAHzCAYAAABcyguCAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XmcXHWd7//Xu0kTQkIkYRGCCZu4sAgCAoJAorKMQkDF\n7TojIo6OI8q9w6B4+Y20zIW5LuO4jLgxCorjxiiLXGXiEEAgEDAMQgBFDFvYkxAIWTv9+f1xzkkq\nldpO96k61dXv5+NRj0rXOXXOt6qr0+/+fJejiMDMzMzMektf2Q0wMzMzs+I55JmZmZn1IIc8MzMz\nsx7kkGdmZmbWgxzyzMzMzHqQQ56ZmZlZD3LIMzPrMZKG0tuMstsCIGkgbc/3ym5LN5L0gfT9mVt2\nW6y3jCu7AWbdTNLQMJ96Q0TMKrQxVpOkmcBM4M6IuLLc1rRHjs/hthHxfPrvblwENVebJO0G/LnG\npvXAMuBu4HLg4ohYN9LGdYHN3h9JHwB2Ba6IiLs63iIb1RzyzBp7itq/mLYD+oHVwHM1ti9pZ6Ns\nEzOBzwCXAD0Z8iosB1Y12J6FwT+QfG57IfhklgJr039vTfIzOCu9nS7pTRGxvKzGtdEHgKOARYBD\nnuXikGfWQETsXOtxSdeT/Mf744j4YEcbZWPZmRHx/WY7RcSrO9GYDgrg7RFxY/aApJ2Ac4GPAQcC\nFwBnlNM8s+7kMXlm1itUdgOscyLiyYj4ODAnfei9ZbanA/z5ttwc8szaQFKfpL+SNEfSM5LWSnpc\n0o8lHVLnOZsMTpf0MUl3Sloh6UlJ35c0vWL/V0i6VNJjklZLukfSh+sce2Z67EXp1ydKmitpWXr8\nWyQ1/SUp6Q3pa3hM0hpJS9LX+J4Wz/sXkn4l6en08TMr9j1K0lck3Za+V2vT/X4l6R01jr1bOlbt\nM+lDp1ZMONhk4kF1O+q0te7g94rj7Srp1en7/qikdZJ+UbXvDpL+SdLd6Xv7Yvq9+T+SpjR7j4tQ\nb+JFjc/Yqen7/YKk5yVdJ+nNDY6b63vUIf+Z3m8rabtaO0jaV9J3JS1Kf1aek3STpI9IqtmjJWlH\nSV9Iv3cvps97NP1Z+WyN9/b69L09tV5DJT2U7nN0sxeVfR5JegwAvlf12a77WTbLuLvWrGCStgF+\nDrwpfWgIeAF4KfAu4BRJZ0bE1xsc48fpvmtIxlXtCPwlcJSkQ4GXA9cAk0nGBI4D9ga+KWlKRHyu\nzqFD0v8EvpS2azkwHjgMOEzS4Wl1pFabPgecnR0HeB54Sfo63yRpNvC+iKg1hjEknQV8oeK8g+lx\nkDQJuL7i2C8AL5KMuzoOOE7StyPibyqOOUgyZnISMJHa4yPXV7ej1mtrcZ8AjgS+BUwgef1rK/eX\n9AaScYFT0sfXpq937/T2V5KOiYg/ttCOkWr4WiVdDHyQ5H18keSzNJPkM/auiPh51f7D+R51QmWx\nYm31RklnAF8hqYQFsIJkTN/h6e3dkt4aEasqnrMrMA/YKX1okOT7vTOwC8nPy+Mkn4VqzT5j0cI+\nACtJPt9TScb/Vo/HfLqFY9gY50qeWfG+TxJ8fgccC2wdEVNIfhn+fyTB4yuSDq/z/JOBvwDeB2wT\nEZNJ/pp/EpgB/F/gx8ANwB4RMZUkVHwzff5AvYoGSVj8HHApsHNEbAfsAPxzuv1jtSp6acXt7LQN\nf00yi3MKScB6T/r4e4BP1TnvTmm7v15x3m2A/0i3rwd+lr727SIiO/5UknFWK4APSzolO2BEPJaO\nmfxi+tCPI2Ja1W1xnfYMh4CLgNuAfdM2TgT+HjYEg6uBbdP99oqICek++5FUnKYDP5c03P97i+qy\nOwn4H8DfAJPT93pP4EaS3wtfk7RF1XNyf4865Nj0/rGIeKFyg6STga+SBNKzgR0i4iUkfxQcDzxA\nEmz/peqY55F8Zh8AjoyILSNie5Jwvx/wj8ATw2xvS9/DiPhp+vmelz50ZtVn+9Bhnt/GkojwzTff\nct5IKhpDwHerHn9z+vi9JAGt1nM/le5zddXjA+njQ8Bf1XjeX1Zsvw/oq9ou4I/p9vdXbZtZ8dxf\n12nX99Ltf6x6fFs2Vm32q/Pcw0hCwBKgv855LxvB+5299utqbMvet+82eH7Wjj832OcDDc6RvYYH\ngPF1nn9Zus8Fdbb3A/+d7vOOnK8/O/9zJIG6+vYEcFzV/uuBGQ0+Y++tcZ6dSarH60nCTce+R3WO\nuVvFazm6Rlu/VrH9E1XbtwAeSrcdU+f4e5CE07XAThWP35se95052np9rZ+9qn0eSvc5Ksdnr+lx\nffOt3s2VPLNiZeNxvhNVVYUK/57ez5RU66/6RyPiBzUe/03Fv78YEZusnRYRAWTjyfapc+4A/qnO\ntgvS+z0l7V/x+DtIKh+/iYi7ax404laSX2DbAgfVOe8X6py3Fb9M7w+t8551yr9GxJrqByVtDbyT\nJFBUV4UAiGQdt6xyWXfcWxOTSSqv1bcdSbrdW/VwRPyoRhufAOaT/MFQ7zNUTzu/RwJ+oWRs6pOS\nngcWk8ysfREYiIivVj1nJknl+56ImEMNEfFnksrsuHT/TLbW4LTCXoFZCTwmz6xYWRfsP0iq13WZ\nmUjShfts1eP31tn/mfQ+gHvq7JON09m2zvZ1wM21NkTEnyQ9SdJNdSAb1+TKXtOb0u31TCH5ZTwd\nuLVq26pospBrOgD+VJKwtD8bxyJV2io9z9JGx2qTYGPXWbWDSNoawD0NMs6E9H64V6L4QLSwhEoL\n7miwLevi3mySSMnfo5fUefwKkiEI1bLP7SuafG4np/fTKx67BjgE+JykvUgWXL41IlbnaK9Z6Rzy\nzIqVrau3La0NwN66xuM1x/pExPqK8FBvPFA20aD6F2/m2YgYbNCmxSQhb/uKx7LXtDUbQ0o9UWef\nhotDp4P6rwVeX3GcVSRXNcgqltkg+ImUE/JgY9Culr1HIqmsNVLvPeqkelVmSCawQNVnqOTvUQCz\nIl0nT9KOJGPq/plk7Oo4Nl9CJfuejCf/9+RzJMF9NvC36W29pPnAL0gq9b248LL1GIc8s2JlQyBO\njoirSm1JcbLX9OWI+LthHqN6lmu1fyAJD88AZ5GMG9xQ4UwnKmThtMzu2nqvI3uPnotkIkwv6prv\nUUQ8DXxf0oMkk0XeLenyiPiPit2y78kVEfH2nMdfC5yczmR/G8nEpwNJXv/rgbPTWdK/H+lrMWsn\nj8kzK9ZT6f2upbaivu3rrQuWysYgVVasOvGa3pnefzwiLqsMD6mdqp+QUxY+tmqwT73uwFZk3YGT\nJU1uuOfo1e7vUW4RcTPJhBeAC6pmLWffk2F/biPitog4JyIOJ+mafi/wCEll8OKq3Vv9jHXjNYWt\nRznkmRXrlvT+L0ptRX39bByrtAlJLyfp4gpgQcWm7DUdLanRL7CReFl63jvrbG80USHrKmxUPVqW\n3u8oqV5X9usaPL+ZO0iqfH0k3Yi9aCTfo3a6gKRdrwDeXfF4Nn5yP0kjnkARESsj4idAtuD4gZIq\nu3izNRqnU0P68zWcPyRa+Xyb1eSQZ1asS9L74yQd12hHSfUmR7STgE/X2ZY9/kBVN9TPSGYwTmXj\n1SVqH3z4V3RYnrbtNTWOOYnkGqX1ZDMhG72ffyRZJqMPOLHGOV5OMot4WCJiBcngfIDz0zbXJGmc\npInDPVeJRvI9aptIFpbOhkZ8smLTfwGPkgxLajizu/pzK2nLBrtnYxYFVO6X/czMrvO8cyqel0f2\n+e7I1VKstzjkmRUoIq4ludpFtuTD30vaMIlB0vaSTpF0DclVJzptJcks2X+TtEPapm3Tq1mcRlIR\nGah8QkQsZWMAPEfSt9MZh6TP31rSLEnfps7M3RZkl6b6kqTsMk5Ieh3JL+tG49yymcZHpmFtM+ny\nJVekX/6LpCOU6JN0LMn1T1cOs+2Zc0gmG7wCuEXScVnVMD3XqySdDfwBOHiE5yrDSL5H7ZaFuP0l\nHQ+QTjA6g+Qz/V5Jv6hcGkjSlpIOk/TPwJ+rjnePpAskHZwFvvR7eAjJ2nwAt1dNvrg8Pdd+kr4s\n6SXp83aU9FWSdQSH8xnLPt9v7+GhANYuZS7SB5xC8h/vIyQf/vuBC4FJVftNIRn/8AzJwpVzSFac\nb+UcOwD/SvJDvDK9/xqwfdmLFPo2em/UWQw53bY1SdAbqrgtI5nRWPnYv1U9b6DeMSv2yZ47o872\nmsegYjFg4BNsXEB2aXqfff3VBuc+t2LfofRncWnVa3qw3nmbvJ+7kyz/kh1nVXr87Dxvps5rJ6nU\nPFCx/WmSNfsWAbtUneOZiv1eTM8zRHJ1ko/TfDHkmu97xX4HA49V7L+WZImcNRWPDWeh4ey5LS2I\n2+C9auUzdkm6z2cK/B41PW+dtuxWccyjmux7U7rf3KrHP0BSfcuOs5Jktvdg5fek6jnLKrYNpvuv\nrXjsKWr8DiKZ7Vv9c599Dt7P8BZDfmVF+9eRzIB/CPhtnvfSt7F5K7uSdxbJh/YcknEs3wA+CszJ\nFtNM768muXTNGSRdKv3AXEm7NDp4xXPfRTIl/niSv/jekz5uNlx1rz8ZydidtwMnkIS9x0gGY/eR\nhJGfkPynXn2N2Favadlon6bHiGTR2Nkkl0WD5JfePOAvI+ITDZ53AcnaaN8m6f6EZNmJx4Bfk1w2\n6sgcba089iKSdckuI/kFKpIAeRnwuojIFoLe7HiRVGzeBPyApHvuJSTjomaQXPWg8hyHAj8iCSsi\n+QPz/wBHsLFbrG4zW3gddwCvIrmqyS3pMSeThKDbSa6henRE/LbZsYZz/hb2b+UzVnOfkXyPWjxv\nszY18/n0/qi0upg8MeISkqD0ZZKq2DqSy/E9Q7J4+GfS7ZVOIlk0/CaSz/fWJEHrrvTxfSJis7Uq\nI+IskuVW7iIJweuBXwFvjGR9w3rvQ93XFxF/AI4h+RlbRrLw9XSSa+iaNaSI8ib6SNouIpZUPfZX\nJNfVfFNEzJV0Esm6RLMi4oZ0n8kkf6VfFhFnNjj+K0ku//SRiPhOxeMfIQmUr4yIB4p+XWbdRtJM\n4DrgoYjYo+TmmJlZB5RayasOeKlsJfZsNtRsYHEW8NLnPU9SiTupySmyv+KrF63Mvi67kmlmZmbW\nFt0Yco5O7+9L7/eh9iWc7gVmpNeMrCki7iUZLPwPkg6SNCkdOPsZ4P+lZXAzMzOzntNVIS8dY3c+\nMCcisnW6prJxjatK2SVzmk0rfxvwMMl4mOdJrqn5J5JJH2ZmZmY9qWsua5aus3QlySyk0yo2DXvQ\nYLr6+eXAAcBHSKqDewOfBS6XdGLUGJQoySuSW6/a3Z9vM7PeExGbrcHYFSEvXTX8apLp8kdHxOMV\nm5dRe/2lqRXb6zmR5MoDb4qIueljN0n6M0k37olsXERzE2VOSBltJPn9ysnvWT5+v/Lx+5Wf37N8\n/H7l1873LF2QZDOld9emi4VeTnLx57dExMKqXRaSjMurtjfwcEQ0Wlxy7/T+jqrHb0/vX5WzuWZm\nZmajQqkhL+1O/SHJgqknR8T8GrtdBexStcL6ZBpU4So8lt5XX5Py0PR+cd42m5mZmY0GZa+T9w2S\nsXIXANdUbX40IhanCxrfRLL449kkF4H+NLAvsH9ELK443iBwSUR8KP16EkklcBzwjySXE3oVcB7J\nwpZ716oESqo1VM/qcNk+P79n+fj9ysfvV35+z/Lx+5Vfu7tra43JKzvkLSJZlb5WZ/JARJyf7jcF\n+CJwMsmVA24B/i4i7q463hBJyPtgxWPTSC6p8yZgZ+BJksuiDUTEE3Xa5ZCXg3/Y8/N7lo/fr3z8\nfuXn9ywfv1/5jbmQ160c8vLxD3t+fs/y8fuVj9+v/Pye5eP3K78yQl7pEy/MzMzMrHgOeWZmZmY9\nyCHPzMzMrAc55NmInXfeeWU3YdTxe5aP3698/H7l5/csH79f+ZXxnnniRQ2eeGFmZmajhSdemJmZ\nmY0hDnlmZmZmPcghz8zMzKwHOeSZmZmZ9SCHPDMzM7Me5JBnZmZm1oMc8szMzMx6kEOemZmZWQ9y\nyDMzMzPrQQ55ZmZmZj3IIc/MzMysBznkmZmZmfUghzwzMzOzHuSQZ2ZmZtaDHPLMzMzMepBDnpmZ\nmVkPcsgzMzMz60EOeWZmZmY9yCHPzMzMrAc55JmZWde69tpr+fCHP8yjjz5adlPMRh1FRNlt6DqS\nwu+LmVm5nn32WaZPn87q1as5/vjj+dWvflV2k8y6kiQiQtWPu5JnZmZd6corr2T16tVAUtFbunRp\nyS0yG10c8szMrCvdfPPNG/4dEdx+++0ltsZs9HHIM7OWPf300yxYsKDsZtgYcccddwBw6KGHAviz\nZ5aTQ56ZtWTFihW87nWv46CDDuJb3/pW2c2xHjc0NMQDDzwAwNvf/nYA/vjHP5bZJLNRxyHPzFpy\n1VVX8cgjjwDwxS9+seTWWK97/PHHWb16NTvssAMHH3ww4JBnlpdDnpm15Lrrrtvw7z/96U9e0sLa\n6sEHHwTg5S9/ObvvvjvAhj8yzKw1Dnlm1pL77rsPgP7+fgDmz59fZnOsx/3pT38CYM8992TnnXcG\n4Mknn2RoaKjMZpmNKqWHPEmnSLpC0iOSVkq6X9KFkiZV7TdF0sWSnpG0QtIcSfu2cPwPSBpqcNux\nfa/OrHf84Q9/ADaOj8rGS5m1Q2Ulb6uttmLq1KkMDg7y7LPPltwys9Gj9JAHnAWsA84Bjge+AXwU\nmCNJAOn91cCxwBnAO4B+YK6kXZoc/5fAYVW3w4ElwPyIeLroF2TWa5599lmWLFnCpEmTOProowGP\nj7L2euqppwDYZZfkv/hp06YByVg9M2vNuLIbAJwQEUsqvr5R0lLgUmAmMBeYTRLMZkXEDQCS5gGL\ngE8CZ9Y7eEQ8C2zyp5+kI4HtgM8U9zLMeldWxXvVq17FXnvtBbiSZ+2VdcuOG5f8mpo2bRr33HMP\njz/+OAcccECZTTMbNUqv5FUFvMwd6f209H42sDgLeOnzniep7p00jNOeCqwBfjSM55qNOdn4qL32\n2sshzzoiC3l9fcmvKVfyzPIrPeTVcXR6f196vw9wT4397gVmSNq61QNLmgC8E/hlRDw3olaajRFP\nPvkkkHSdTZ8+nfHjx/PUU0+xYsWKkltmvapeyHviiSdKa5PZaNN1IS8dY3c+MCcisuXNpwLLauye\nXchwSo5TnAxsQ9IdbGYtyELeTjvtRF9fH9tssw0Aq1atKrNZ1sOqQ142w9aVPLPWdVXIS2fUXgms\nBU6r2BQFnuZU4Cng/xV4TLOeVhnyYOMvXi9nYe1SHfK23357AJYuXVr3OWa2qa4JeWk36tXAbsBx\nEVH559oykmpetakV21s5x87Am4B/j4iGv50kNb0NDAy0clqzUa9eyFu/fn1pbbLeln22ss/a1KnJ\nf/cOeTaWDAwMtJRH6umG2bVI6gcuBw4EjomIhVW7LCRZPqXa3sDDEbGyxVP9JbAFLXTVRhRZPDQb\n3apD3hZbbAG4kmftU13Jc8izsWhgYKClglK9oFd6JU9SH/BDkuVSTo6IWsvoXwXsIumoiudNBk5M\nt7Xq/cBdEfH74bfYbOxxd611mkOe2ciVHvKArwOnAF8CVkk6rOKWLXR8FTAPuEzSuyUdlz4WwOcr\nDyZpUNLF1SeRdCDJLF1PuDDLYc2aNTz33HOMGzeOKVOSOU4OedZuDnlmI9cNIe94krB2LnBL1e10\ngEj6Tk8A5gAXAT8nuUrGrIhYXHW8Pmq/rlPT5/yw+Jdg1ruyKw/suOOOG37hekyetVsW8rKhAZMn\nT6avr4/nn3+edevWldk0s1Gj9DF5EbF7i/stIwl9pzfZr2ZwjYgzaXBlDDOrLaucZLMbwZU8a7/q\nSl5fXx9TpkxhyZIlPPfcc+ywww5lNs9sVOiGSp6ZdbFly5LJ61lXLXjihbVfdcgDd9ma5eWQZ2YN\n1Qp5ruRZuznkmY2cQ56ZNZSFvG233XbDYw551m4OeWYj55BnZg01quR54oW1S62Ql30GHfLMWuOQ\nZ2YNeUyelaFWyHvJS14CwAsvvFBKm8xGG4c8M2vIY/KsDLVC3jbbbAPA888/X0qbzEYbhzwza8gh\nz8pQK+RNnjwZcMgza5VDnpk15DF5VgaHPLORc8gzs4ZcybMyOOSZjZxDnpk15IkXVoasSuyQZzZ8\nDnlm1tDy5csBr5NnneVKntnIOeSZWUPZchXZzEZwyLP2axTyvISKWWsc8sysrvXr17Nq1SoksfXW\nW2943BMvrN1cyTMbOYc8M6trxYoVAEyaNAlJGx73mDxrN4c8s5FzyDOzuipDXiV311q7eTFks5Fz\nyDOzumqNxwOHPGu/7LOVVY0BJk6ciCRWrlzJ4OBgWU0zGzUc8sysrizk1avkeUyetUutSp4kT74w\ny8Ehz8zqyrprXcmzTqsV8sDj8szycMgzs7rqVfI88cLarV7ImzhxIgAvvvhix9tkNto45JlZXa7k\nWVkc8sxGziHPzOpqNibPIc/apVnIW7lyZcfbZDbaOOSZWV3NKnmeeGHtUi/kZYtyu5Jn1pxDnpnV\n5SVUrCzurjUbOYc8M6ur3mLInnhh7ZZViR3yzIbPIc/M6nIlz8riSp7ZyDnkmVldzS5r5jF51i4O\neWYj55BnZnW5kmdlaTbxwrNrzZpzyDOzurwYspXFlTyzkXPIM7O6vBiylcUhz2zkHPLMrC4vhmxl\nccgzGzmHPDOry4shW1kc8sxGziHPzOpyJc/Kkn22svGfGV/xwqx1DnlmVlNEeDFkK42vXWs2cg55\nZlbTqlWrGBoaYvz48fT392+yzZU8azd315qNXKkhT9Ipkq6Q9IiklZLul3ShpElV+02RdLGkZySt\nkDRH0r45zrOLpO9KekLSakl/lnRh8a/IrHdklZLsl2olj8mzdnPIMxu5cSWf/yzgMeCc9P61wAAw\nS9LhERGSBFwNzADOAJ4DPg3MlXRARCxudAJJuwE3Aw8CHweeAnYH9mzD6zHrGatWrQJgwoQJm21z\nJc/azSHPbOTKDnknRMSSiq9vlLQUuBSYCcwFZgOHA7Mi4gYASfOARcAngTObnOObwKPp87Oyw28L\newVmPSqr5GUD3St5TJ61U0QQEQAkf+dv5IkXZq0rtbu2KuBl7kjvp6X3s4HFWcBLn/c8SXXvpEbH\nl7QncCzwtYqAZ2YtaBTyXMmzdsoCXnUVD1zJM8ujGydeHJ3e35fe7wPcU2O/e4EZkjb/DbTREen9\n6nQc32pJSyVdKmlqQe0160nurrWyZGM9a4W8CRMmIIk1a9Z4TKhZE10V8iTtApwPzImIBenDU4Fl\nNXZfmt5PaXDIrBr4XeB+4HjgU8BbgWtV3Q9gZhu0UsnzL1lrh3rj8SDpvs0+k15Gxayxrgl56Yza\nK4G1wGkVm2IEh81e39yI+HhEXB8R3wH+FjgIOK5Be5reBgYGRtA0s+7m7lorS6OQB+6ytbFjYGCg\npTxST1eEPEkTSMbY7QYcFxGPV2xeRlLNqza1Yns92Zi/OVWPZ1/vX++J2cDfRjeHPOtlWcir1V3r\niRfWTs1Cnidf2FgxMDDQUh6pp/SQJ6kfuBw4EHhLRCys2mUhybi8ansDD0dEo3p9rbF8ZtaCbEye\nK3nWaa7kmRWj7MWQ+4AfkiyXcnJEzK+x21XALpKOqnjeZODEdFsjtwJPkozFq5R9ffswmm02JnhM\nnpXFIc+sGGVX8r4OnAJ8CVgl6bCK2y7pPlcB84DLJL1b0nHpYwF8vvJgkgYlXZx9nS6bcg7wVknf\nkHSspL9Nzzs3Iq5r+ys0G6U8Js/K4pBnVoyyQ97xJGHtXOCWqtvpAJF0Np9AMo7uIuDnwDqSxY2r\nr3bRR9VriojvA+8H3kASDs8FfkBSCTSzOhotoeIxedZOzUJe9plcvXp1x9pkNhqVesWLiNi9xf2W\nkYS+05vsV/N/hIi4DLgsdwPNxjBX8qwsrYa87A8RM6ut7EqemXUphzwrS/a5yirG1RzyzFrjkGdm\nNbVyxQtPvLB2aFbJ22qrrQCHPLNmHPLMrCZX8qwsHpNnVgyHPDOrqVHI88QLayePyTMrhkOemdXU\nSnetQ561g0OeWTEc8sysJi+GbGXxmDyzYjjkmVlNHpNnZcn+ePCYPLORccgzs5o8Js/K4u5as2I4\n5JlZTR6TZ2Vxd61ZMRzyzKwmd9daWVzJMyuGQ56Z1eSJF1YWr5NnVgyHPDPbTES4u9ZK40qeWTEc\n8sxsM+vWrWP9+vWMGzeO/v7+zbZ74oW1k8fkmRXDIc/MNtOoqxZcybP2ciXPrBgOeWa2mUZdteCQ\nZ+3lMXlmxXDIM7PNtFrJ88QLawdX8syK4ZBnZpvJQl69Sp7H5Fk7ZZ+r7HNWzWPyzFrjkGdmm8m6\nwTwmz8rg7lqzYjjkmdlmsl+eWcWkmkOetZO7a82K4ZBnZpvJQt748eNrbveYPGunZiGvv7+fvr4+\nBgcHGRwc7GTTzEYVhzwz24wreVamZiFPksflmbXAIc/MNtMs5HnihbVTViGuF/LA4/LMWuGQZ2ab\ncSXPytSskgcel2fWCoc8M9uMQ56VySHPrBgOeWa2mTVr1gDNQ54nXlg7tBLyPCbPrDmHPDPbjMfk\nWZnyVPI8Js+sPoc8M9uMu2utTO6uNSuGQ56ZbcYhz8rkkGdWDIc8M9uMF0O2MnlMnlkxHPLMbDOu\n5FmZPCbPrBjj6m2QtAiIER7/yxHx1REew8w6zBMvrEzurjUrRt2QB+wKLE9vwzED2HaYzzWzErmS\nZ2Vyd61ZMRqFPIB/iYjzh3NgSf7f32yUcsizMmWfq6xiXEv22czWdDSzzZU+Jk/SKZKukPSIpJWS\n7pd0oaRJVftNkXSxpGckrZA0R9K+LZ7jIUlDNW6z2/OqzEY3L4ZsZWqlkpdNCnLIM6uvUSXvEGDx\nCI7d6vPPAh4DzknvXwsMALMkHR4RIUnA1SRdwGcAzwGfBuZKOiAimp0ngF+nx630x9ZeitnY4jF5\nVqY8Ic8TL8zqqxvyIuKOkRw4x/NPiIglFV/fKGkpcCkwE5gLzAYOB2ZFxA0AkuYBi4BPAme2cJ5n\nI2J+i20CtX32AAAgAElEQVQyG9PcXWtlyirErYzJcyXPrL7Su2urAl4mC4jT0vvZwOIs4KXPe56k\nundSC6dRejOzFrS6Tp5DnrWDu2vNitEw5Enqk/QzST+WtGWD/baU9FNJPymoXUen9/el9/sA99TY\n715ghqStmxwvgBMlvShptaR5kloJh2ZjUquVPI/Js3bIM7vW3bVm9TWr5L0jvV0REWvr7ZRuuwJ4\np6RTRtIgSbsA5wNzImJB+vBUYFmN3Zem91OaHPZqkrF8xwLvA1YDv5D0vpG01axXubvWyuRKnlkx\nmoW8d5FMhmilQvdjkokW7x1uY9IZtVcCa4HTKjaNaFHmiPhERFwWETdHxH8AbyLpEr5wJMc161We\neGFlcsgzK0azkPc6YG5ENA1ZETEEXAccPJyGSJpAUnHbDTguIh6v2LyMpJpXbWrF9palbb0cmC7p\npXXa0/Q2MDCQ57Rmo4YreVYmd9eaJQYGBlrKI/U0Wwx5J5JKXqsWAzVDUyOS+klC14HAMRGxsGqX\nhSRdrdX2Bh6OiJV5z9lMC7nWrGc55FmZXMkzSwwMDLRUUKoX9JpV8tYCtafX1TYeWJdjfyT1AT8k\nWS7l5DrLnFwF7CLpqIrnTQZOTLflImkc8G6SgPhU3ueb9TovhmxlcsgzK0azSt4TwP45jvca4PGm\ne23q68ApwAXAKkmHVWx7NF3o+CpgHnCZpLPZuBhyAJ+vPJikQeCSiPhQ+vV7gROAa9K27QR8DDiA\nEYwfNOtV69evZ926dUiiv7+/5j4ek2ft5O5as2I0q+TdBBwtaa9mB5L0cpJq3G9ztuF4krB2LnBL\n1e10gHRM4AnAHOAi4OckFcNZNa520cemr+vPJMHuS8B/At8AVgHHR8RPc7bVrOdVVvHqdgG4u9ba\nyJU8s2I0q+R9i2SW688kHVeva1PSjsBPgS2Ab+dpQETs3uJ+y0hC3+lN9uur+vo2ktm0ZtaCZgsh\ng0OetZcva2ZWjIYhLyLmS/oW8BHgHknfBv6LjZMxXkYSoD4MbAd805cOMxvdmk26AI/Js/bK013r\nSp5Zfc0qeQCfIKnQfYhkHFw2Fg42vVTYd9J9zWwUyxPyXMmzdnB3rVkxml67NiLWRcSHgSOBHwCL\ngDXpbRHwfeANEfGRiBhsZ2PNrP1aCXmeeGHtlH2uss9ZLZ54YdZcK5U8ACLiZuDmNrbFzLqAK3lW\nNlfyzIrRtJJnZmNLszXywCHP2isb6+mQZzYyDnlmtolWKnnZ0ioR4avDWOHyVvL8GTSrre5PkKRV\nkv73cA880uebWTlaDXmu5lm7tBLy+vr6NizWvXbt2o60y2y0aVTJG0+OMXtteL6ZlaCVkAfusrX2\naSXkgbtszZppFsJOlrTbMI5be5l8M+t6rSyGDA551j6thrytttqKFStWsHr1aiZPntyJppmNKs1C\n3gHpbbgc9sxGmbyVPC+IbEVzJc+sGI1C3h4FHH9ZAccwsw5yd62VLW/I81p5ZrXVDXkR8VAH22Fm\nXaLVkOcFka1d8nTXgit5ZvV4CRUz24QreVY2d9eaFcMhz8w20cpiyOCQZ+2Tt5Ln7lqz2hzyzGwT\nnnhhZXMlz6wYDnlmtgmPybOyOeSZFcMhz8w24TF5VjZ315oVwyHPzDbhxZCtbK7kmRXDIc/MNuEx\neVa2LORlQwLqccgzayz3tWUl7Qi8A3g1MDEiTk8f3wHYHbgnIlYW2koz6xh311rZsj8c3F1rNjK5\nQp6kDwFfBbL//QM4Pf33TsCtwIeBi4tqoJl1VqvdtZ54Ye3i7lqzYrTcXSvpGOBbwB+AtwHfoOLa\ntBFxN7AQOKngNppZB3mdPCubJ16YFSNPJe9TwJPAzIhYLum1Nfb5PXBYIS0zs1KsXbsWgC233LLh\nfg551i6u5JkVI8/Ei4OBX0bE8gb7PAbsPLImmVmZsl+Yrc6u9cQLK5pDnlkx8oS8LYEVTfbZFvD/\n+GajWKshz2PyrF3cXWtWjDwh72HgoCb7HEIyZs/MRqm8lTyHPCuaK3lmxcgT8q4AjpL0rlobJZ0G\n7A/8RxENM7NyeEyelc0hz6wYeSZefAF4D/Dvkt5B0jWLpDOAo4C3Aw8AXyu6kWbWOa7kWdncXWtW\njJZDXkQslTQTuBR4Z8Wmr6b3vwX+R0Q0G7dnZl3MEy+sbK7kmRUj12LIEfEwMFPS/sDrge2A5cC8\niPhdG9pnZh3miRdWtrwhz5U8s9paDnmS5gI3RcQ/RMRdwF3ta5aZlcVj8qxsebtrXckzqy3PxItD\ngcZXizazUW1oaIjBwUEA+vv7G+7rkGft4u5as2LkCXl/Aqa3qyFmVr7KrlpJDff1mDxrF0+8MCtG\nnpD3HeAESbu2qzFmVq5Wx+OBK3nWPtkfDq7kmY1MnpD3S5IZtDdJ+rikQyXtKmlG9a3VA0o6RdIV\nkh6RtFLS/ZIulDSpar8pki6W9IykFZLmSNo3R9uz47xH0pCkR/M+12wsaHU8HnjihbVP9pnKPmP1\nOOSZNZZndu2DFf/+SoP9gtbH7p1Fcr3bc9L71wIDwCxJh0dEKOkzuhqYAZwBPAd8Gpgr6YCIWNzK\niSRtC3wZeDJto5lVcSXPuoG7a82KkSfkfb/F/fIEqBMiYknF1zdKWkqyFt9MYC4wGzgcmBURNwBI\nmgcsAj4JnNniuT4P3EkS8t6co41mY4ZDnnUDT7wwK0aexZA/UPTJqwJe5o70flp6PxtYnAW89HnP\nS7oaOIkWQp6kI4D3AfsBnxlRo8162HBCnideWNGykNds8o9DnlljecbkdcrR6f196f0+wD019rsX\nmCFp60YHk9QPfBv4fET8ubBWmvWg7Jelx+RZmVodk+fuWrPGuirkSdoFOB+YExEL0oenAstq7L40\nvZ/S5LCfAvqBfyqkkWY9LJt44e5aK1Ors2uzP0ayz62ZbSrPFS++R4vj7SLig3kbks6ovRJYC5xW\nebi8x6o45suB/w2cHBGV/ws0PWazbgKA8847j4GBgeE2z6zreEyedYNWK3mVIS8iWvp/22w0GRgY\n4LOf/eywn59n4sWpOfbNFfIkTSCZQbsbcHREPF6xeRlJNa/a1Irt9XwVuA64LZ1dC7Al0CfpJcCa\niKhZ54/wBFwbezwmz7pBq5W8vr4+xo0bx+DgIOvWrWtpmIHZaDIwMNBSManeHzh5Qt4edR7fFjiY\nZELDLSTdoy1Lx8xdDhwIHBMRC6t2WQgcW+OpewMPR8TKBod/NbArtYPgMpIlVf4uT3vNepkredYN\nWq3kQfJZHRwcZM2aNQ55ZlXyzK59qMHm/5Z0LfB74DfAxa0cU1If8EOS5VJOiIj5NXa7CjhN0lER\ncWP6vMnAicBlTU7xHqDyt5VI1uQ7CDgFaGmNPbOxwoshWzdodQkVSELeiy++yJo1a9hmm23a3TSz\nUSVPJa+hiHhU0i+BT9BiyAO+ThK2LgBWSTqsYtuj6ULHVwHzgMsknc3GxZCDZO27DSQNApdExIfS\nNt1WfUJJp5F0096Y5/WZjQWu5Fk3aLW7Fjz5wqyRomfXPgW8Isf+x5OEtXNJunorb6cDRDI47gRg\nDnAR8HNgHcniyNWVuD6av6bAV7wwq8khz7pB3u5a8Fp5ZrUUVsmTtAUwC1je6nMiYvcW91tGEvpO\nb7Jf09AaEac128dsrPLEC+sGruSZFSPPEipHNTjGDJJlT15L6121ZtZlPCbPuoEreWbFyFPJu76F\nfW4Ezh5eU8ysbO6utW6Qp5KXfVZdyTPbXJ6Qd36dx4dIliO5rc7sWDMbJRzyrBvkmV2bVZ1dyTPb\nXJ4lVAba2A4z6wIek2fdwN21ZsVoeXatpKMkzWiyz4wGY/fMrMvlGZPnSp61iydemBUjzxIq1wMf\naLLP+4G5w22MmZUrTyXPEy+sXVzJMytG0evk+erQZqOYx+RZN/DEC7NiFB3yZgAvFHxMM+sQhzzr\nBnkqeZ54YVZfw4kXks4juTpEVqGbKdUs1m0B7EpyrdibimygmXXOcMbkeeKFFS3vtWvBIc+slmaz\na8+r+npmeqtnMXDOCNpjZiXymDzrBp54YVaMZiHvjRX/vg64NL1VWw8sAe6PCP+PbzZKubvWuoEn\nXpgVo2HIi4jrs39L+j5wReVjZtZbHPKsG7iSZ1aMPIshf6CN7TCzLuAxedYNXMkzK0bRs2vNbBRz\nJc+6wXAmXriSZ7a5XCFP0jRJF0l6UNIqSeurbkOS/Ge92SjliRfWDYbTXetKntnmWu6ulbQLcDuw\nI3AvMB54GFgL7EGyjMp/A8uLb6aZdYIredYN3F1rVow8lbzPAC8F/iIiXpM+9r2IeCWwO3AtMAF4\nR7FNNLNO8bVrrWyVn6c667JuwhMvzOrLE/KOA66NiDnVGyLiMeCdwNbAZwtqm5l12HAqeZ54YUXK\nU8UDV/LMGskT8nYC7qn4ej1J5Q6AiFgBzAFmF9M0M+s0j8mzsuUZjweeeGHWSJ6Q9wJQ2YfzHLBL\n1T7LScbsmdko5DF5VrY8M2vBEy/MGskT8h4Gpld8fRfwRkkTAST1AccAjxXXPDPrJI/Js7K5u9as\nOHlC3m9IQl1/+vUlwDTgFklfAG4B9gV+UmgLzaxjPCbPypa3u9YTL8zqa3kJFeC7JF20OwCPR8Rl\nkg4CPgHsl+7zY+CCYptoZp2wfv161q9fjyTGjWv+X4MredYOruSZFSfPZc3+CPzfqsf+l6R/Ilkn\nb1FEPFVw+8ysQyqreK0sXeGJF9YOnnhhVpw8iyGfCjwZEddWPh4RTwNPF90wM+usPOPxwJU8a4+8\nlTxPvDCrL8+YvH8Djm9XQ8ysXHnG44FDnrVH3tm1ruSZ1Zcn5D2Vc38zG0WGG/I88cKKNNyJF67k\nmW0uT2j7FTArXSrFzHpM3pDnMXnWDp54YVacPIHtXGAb4LuStm9Te8ysJFl3l7trrUxeQsWsOHmW\nUPkx8DzwfuDdkh4CngSieseIeGMhrTOzjskqIZ54YWVyJc+sOHlC3tEV/x4PvDK9mVkP8Jg86wae\neGFWnDzr5HksnlkP8+xa6waeeGFWHAc3MwPyj8nzxAtrh+Guk7du3Tp/Fs2qlB7yJJ0i6QpJj0ha\nKel+SRdKmlS13xRJF0t6RtIKSXMk7dvC8SdJ+qmkB9LnLZN0m6T3te9VmY0+HpNn3SBvJU8S/f3J\nJdXXrVvXtnaZjUa5Qp6kLSR9Ig1Jz0taX7HttZIukvSKnG04C1gHnEOy2PI3gI8Cc5ReWym9vxo4\nFjgDeAfQD8yVtEuT42+ZHv9C4ETgvcB9wA8knZmzrWY9y9211g3yVvLAky/M6slzWbMtgV8DM4El\nwAtAZbXtIeCDwLPAZ3K04YSIWFLx9Y2SlgKXpueaC8wGDgdmRcQNaXvmAYuATwJ1w1pELAWqq3a/\nTsPoB4Gv5GirWc/yxAvrBnkreZB8ZlesWOHJF2ZV8lTyziYJXZ8FdgK+U7kxIpYBvyWptrWsKuBl\n7kjvp6X3s4HFWcBLn/c8SXXvpDznq7AU8G8ns5TH5Fk3yDu7Fjz5wqyePCHvfcAtEfHZiKgXjhYB\nM0berA3LtdyX3u8D3FNjv3uBGZK2buWgksZJ2k7Sh0nC6FdH3FKzHuExedYNRtJd60qe2abyhLzd\ngXlN9lkKbDf85kA6xu58YE5ELEgfngosq3M+gCktHPcMYC3wDPB14KyIuGQkbTXrJR6TZ91gON21\nruSZ1ZZnMeQ1wLZN9pkOPDfcxqQzaq8kCWOnVWza7Koaw/Bj4BZge5Iu3n+RtCYivl3Asc1GPY/J\ns27giRdmxclTybsTOFZSzd8Akl4CHAfMH05DJE0gGWO3G3BcRDxesXkZSTWv2tSK7Q1FxLMRsSAi\n/jMiPgb8APiipJr/k0hqehsYGMjxCs26m69da91guBMvwN211nsGBgZayiP15Al53yap1P1Q0uTK\nDZKmAJeQhK5v5n0RkvqBy4EDgbdExMKqXRaSjMurtjfwcESszHtO4Hcks4NfWmtjRDS9OeRZL8k7\nJs8TL6wdhlPJc3et9aqBgYGW8kg9LYe8iPgRSZB7O8m4to8CSLoDeIKkC/SiiLgmzwuQ1Af8kGTm\n7skRUasSeBWwi6SjKp43mWTdu6vynK/C0STLwDw9zOeb9RSPybNuMJzZta7kmdWWZ0weEfFBSTeS\nrEu3f/rwgSSVti9FxPeG0YavA6cAFwCrJB1Wse3RiFhMEuTmAZdJOptk3N+nScbqfb7yYJIGgUsi\n4kPp1x8BDgV+AywmmRjyLpIFlT8VEYPDaLNZz3HIs27giRdmxckV8gDSGamXpMuWTAGWR8SKEbTh\neJKwdm56qzQAnB8RIekE4IvARcBWJJMoZqUhsFIfm1Yof0+yzt4XSbqTnyVZeuWtEfGrEbTbrKcM\nd0yeJ15YkTzxwqw4uUNeJh0HN5yxcNXH2b3F/ZYBp6e3Rvv1VX09D3jrsBtoNkZ4TJ51g5FU8txd\na7ap3CFP0jbA24ADgJcAy0lm3v5ihBU9MyuRu2utG7iSZ1acXCFP0rtIZs/WWi/vK5I+EhE/K6Rl\nZtZRDnnWDTzxwqw4LYc8SccA/w4MAZcCNwBPklzHdibJZc/+XdJzETGn+KaaWTt5nTzrBp54YVac\nPJW8z5BcieLIiPhd1bZLJP0r8Nt0P4c8s1FmuNeu9cQLK5KvXWtWnDyLIb8W+EmNgAdARNwB/CTd\nz8xGmbzdtZ54Ye3gSp5ZcfKEvLXA4032eSLdz8xGGY/Js27giRdmxckT8m4Ejmiyz+HpfmY2ynhM\nnnUDX7vWrDh5Qt45wGskfU7SxMoNkiZJ+jywH/CpIhtoZp3hMXnWDYYzu9bdtWa15Zl48SmSq0ec\nDfy1pAXAU8BLSS5tti1JFe9TkjZ5YkR8sJDWmlnbeEyedQNPvDArTp6Qd2rFv7cF3lhjn6PSWzWH\nPLMu5zF51g088cKsOHlC3h5ta4WZlc5j8qwbeOKFWXFaDnkR8VAb22FmJRvumDyHPCuSJ16YFSfP\nxAsz62HD7a71xAsr0nAqee6uNast17VrASRtAUwDXgb019onIryMitko44kX1g187Vqz4uS5dq1I\nZtb+PbB9g10DaP1PMDMr3YsvvsjQ0BDjxo1j3LjW/ltwd621gydemBUnTyXvPJLr0i4BLgUWA4M1\n9osC2mVmHbRgwQIA9ttvv5af45Bn7eAlVMyKkyfknQ4sAg6MiOVtao+ZlWD+/PkAvO51r2v5OR6T\nZ+3gSp5ZcfJMvNgOuNIBz6z3ZCHvkEMOafk5HpNn7eAlVMyKkyfkPQhMaVdDzKw8t99+O5Av5Lm7\n1trBEy/MipMn5H0DOFHSzu1qjJl13jPPPMOiRYuYOHEie++9d8vPc8izdhjJOnkrV65sS5vMRqs8\niyFfJOmVwE2S/hH4HVCz6zYiHimofWbWZlkV78ADD8zVReaQZ+0wnO7a6dOnM27cOB566CFWrFjB\npEmT2tU8s1El72LId5Esn/Ld9N8P1bgtKqhtZtYBwxmPB554Ye0xnErehAkTeM1rXsPQ0BC/+93v\n2tU0s1Enzzp5fw18C1gHXA88jpdQMRv1hjMeDzzxwtpjOJU8SD6/CxYsYP78+Rx99NHtaJrZqJNn\nCZWzgKeB10eEq3VmPSAiRlzJc8izIg2nkgdw6KGH8s1vfpPbbrutHc0yG5Xy/BTtCvzMAc+sdzz0\n0EM8++yzbL/99uy66665nuuQZ+0wnNm1sPGPlOyPFjPLF/Iep861as1sdKqs4iVXLmydx+RZOwy3\nkveqV72KbbbZhkcffZQnnniiHU0zG3Xy/BRdCrxV0uR2NcbMOmu4XbWwccyUQ54VKfs8tXoN5Uxf\nX9+GK7a4mmeWyBPyLgTmA3MkzZK0TZvaZGYdMtxJF5D8Us2qf+6ytaIMN+SBu2zNquX5KapcSvy/\ngKjRvSMgIiLftCgz67jBwcENy03kuWZtpS222ILBwUEGBwc3XD/UbCQGB5NFG/LOroVk8gXgyRdm\nqTwh78YW9/MSKmajwL333svKlSvZfffd2X777Yd1jHHjxjE4OOguWytMEZW822+/naGhodzj+sx6\nTZ4rXsxsYzvMrMNGMh4vk1VbsuqL2Uhln6XhhLxp06bxspe9jMcee4w//OEPvPrVry66eWajiv/M\nMRujRjIeL5P9InbIs6KMpLsWPC7PrNKwQp6kiZIOlHRk0Q0ys87IfgkOdzwebAx57q61ooykuxYc\n8swq5Qp5kqZL+jnwHHAHyeXNsm1HSrpX0sxCW2hmhVu5ciV33303fX19HHjggcM+jrtrrWgjreR5\n8oXZRi2HPEk7A7cCs4FfAvNIZtNmbgNeCrw7TwMknSLpCkmPSFop6X5JF0qaVLXfFEkXS3pG0gpJ\ncyTt28LxXyHpa2kAfUHS45KulPSaPO006yV33nkn69evZ99992XixInDPo4reVa0kVbyDjroICRx\n1113sXr16iKbZjbq5KnknUcS4o6NiLcBcyo3RsRa4LfAETnbcBawDjgHOB74BvBRkvX4BJDeXw0c\nC5wBvIPk6htzJe3S5PjHArOA7wInAn8L7ADcKmn4JQyzUayI8XjgSp4Vb6SVvG222YZ99tmHwcFB\n7rzzziKbZjbq5PlT6S3AVRFxXYN9HgHekLMNJ0TEkoqvb5S0lOQKGzOBuSTVw8OBWRFxA4CkecAi\n4JPAmQ2O/6OI+NfKByRdBzyUPu/UnO01G/WKGI8HnnhhxRvJ7NrMIYccwj333MP8+fN5/etfX1TT\nzEadPJW8lwJ/bLLPOmBSk302URXwMnek99PS+9nA4izgpc97nqS6d1Le46fPfaDi+GZjShHLp4C7\na614I+2uBU++MMvkCXnLgOlN9tkLeHL4zdng6PT+vvR+H+CeGvvdC8yQtHWeg0uaCuxbcXyzMWPJ\nkiU8+OCDTJgwgX322WdEx3J3rRVtpN214MkXZpk8Ie8mYHY6AWMzkvYiGVM3dyQNSsfYnQ/MiYgF\n6cNTSUJmtaXp/ZScp/kayZU5vjysRpqNYtl4vAMPPJD+/v4RHcuVPCtaEZW8ffbZhwkTJvDggw+y\nZEmtziKzsaFhyJN0asUs1C8AE4AbJP1F+m8kTZL0FpIZtwH883Abk86ovZLkOrmnVWwq7FJpkj4N\nvBc4IyL+3GC/preBgYGimmXWMUVNugBX8qx4RVTy+vv7NywNlH3ezUajgYGBlvJIPc0qed8DTgaI\niNuADwO7AdcAZ6f7LCcJeLsBH4yIWt2qTUmaQDLGbjfguIh4vGLzMpJqXrWpFdtbOcffABcA50bE\nJY32jYimN4c8G42KmnQBruRZ8Yqo5IG7bK03DAwMtJRH6sn1UxQR35V0E8kSJ68HtiMJefOAf42I\nPwznRUjqBy4HDgSOiYiFVbssJFkKpdrewMMRsbKFc/wV8HXgixHxT8Npp9loFxGFTboAz6614hVR\nyQNPvjCDnCEPICL+CPyvohogqQ/4IclyKSdERK2fyKuA0yQdFRE3ps+bTLLu3WUtnONtJOvkfSci\nPllU281Gm0ceeYSnn36aqVOnsscee4z4eO6utaIVsYQKbFrJi4iGXVpmvWpY164t2NeBU4AvAask\nHVZxyxY6voqkWniZpHdLOi59LIDPVx5M0qCkiyu+Pgr4EXAXcGnV8V/b/pdn1j0qx+MV8UvP3bVW\ntKK6a3fddVd22GEHlixZwqJFi4pomtmo08pP0baSZuQ5aEQ8kmP340nC2rnprdIAcH5EhKQTgC8C\nFwFbAbeQLI68uOo5fWwaXmcBWwKvBW6u2vchYOTlDLNRIhufVMR4PHAlz4pXVHetJA455BCuueYa\n5s+fX0jl2my0aSXk/U8aX1GikkgCW8s/nRGxe4v7LQNOT2+N9uur+vqzwGdbbY9ZL7v11lsBOOyw\nwwo5nit5VrSiKnmQdNlec8013HbbbbznPe8Z8fHMRptWfoqWp7dWFbbciZkVZ926ddxxR3IxmWy8\n0kh54oUVrahKHnjyhVkrIe/LaTXMzEaxu+66i9WrV/OKV7yC7bbbrpBjurvWilZkJS8blrBgwQLW\nrVs34sW/zUabViZeuDJn1gOK7qoFd9da8Yqs5E2dOpW99tqL1atXc/fdd4/4eGajTTfMrjWzDmhH\nyHMlz4pW1BIqGXfZ2ljmkGc2RmQh7/Wvf31hx3Qlz4pWZHct+MoXNra1EvK8gqTZKPf000/z4IMP\nsvXWW7PvvvsWdlxX8qxoRXbXgit5NrY1/FOpejkSMxudKtfHK6pCAp5da8UrupJ3wAEH0N/fz333\n3cfzzz/P5MmTCzmu2WjgEGc2BrRjPB64u9aKV3Qlb/z48RxwwAFExIYlhMzGCoc8szGgHePxwN21\nVryiK3ngLlsbuxzyzHrc+vXrN/xyK2oR5IwreVa0omfXgidf2NjlkGfW4xYuXMiKFSvYbbfd2Gmn\nnQo9tit5VrSiu2thYyXv1ltvJcJLv9rY4ZBn1uPa1VULnnhhxWtHd212lZcnn3ySRYsWFXZcs27n\nkGfW49o16QLcXWvFiogNn6UiK3mSOOKIIwC4+eabCzuuWbdzyDPrcfPmzQPaE/LcXWtFGhoaAqCv\nrw+p2CVaHfJsLHLIM+thS5Ys4f7772errbbigAMOKPz4ruRZkdoxHi+Thbybbrqp8GObdSuHPLMe\nllUtDjnkELbccsvCj+9KnhWpHePxMgcffDDjx49n4cKFLFu2rPDjm3UjhzyzHpZVLY488si2HN+V\nPCtSO5ZPyYwfP56DDz4YgFtuuaXw45t1I4c8sx7229/+FoA3vOENbTm+Z9dakdrZXQsbfw48Ls/G\nCoc8sx61cuVKfve739HX18fhhx/elnO4u9aK1M7uWvC4PBt7HPLMetTtt9/OunXreM1rXtO2i7K7\nu9aK1O5KXvbHzu23387atWvbcg6zbuKQZ9ajsq7ado3HA1fyrFjtruRtt912vPrVr2b16tUsWLCg\nLbAbv8YAACAASURBVOcw6yYOeWY9KuuSatd4PHAlz4rV7koeuMvWxhaHPLMetH79+g0zCNsZ8lzJ\nsyK1c3ZtxpMvbCxxyDPrQb///e954YUX2GOPPZg2bVrbzuPZtVakdnfXwqaVvIho23nMuoFDnlkP\navfSKRl311qROtFdu+eee7Lzzjvz7LPPsnDhwradx6wbOOSZ9aAbb7wRaO+kC3B3rRWrE5U8Scya\nNQuAuXPntu08Zt3AIc+sxwwNDXH99dcDMHPmzLaey5U8K1InKnmAQ56NGQ55Zj3mnnvuYcmSJUyf\nPp0999yzredyJc+K1IlKHmwMeTfccANDQ0NtPZdZmRzyzHrMddddByS/yCS19VyeeGFF6lQlb489\n9mD69OksXbqU3//+9209l1mZHPLMekzWBZVVK9rJ3bVWpE4soQIel2djh0OeWQ9Zv349N9xwA9CZ\nkOfuWitSp7prwePybGxwyDPrIXfeeSfLly9njz32YNddd237+VzJsyJ1qrsWNoa8G2+80Z9f61kO\neWY9JBuP98Y3vrEj53Mlz4rUyUrerrvuyu67787y5cu58847234+szKUHvIknSLpCkmPSFop6X5J\nF0qaVLXfFEkXS3pG0gpJcyTt2+I5/k7S1ZKekDQk6bz2vBqzcnVyPB544oUVq5OVPHCXrfW+0kMe\ncBawDjgHOB74BvBRYI7SqYHp/dXAscAZwDuAfmCupF1aOMeHgO2BX6Rf+1o21nPWrFmzYRHkToW8\nLbfcEoC1a9d25HzW29atWwds/Fy1W1bx/s1vftOR85l1Wvtr4s2dEBFLKr6+UdJS4FJgJjAXmA0c\nDsyKiBsAJM0DFgGfBM5sdIKI2Dt9zhbA3xT9Asy6wU033cTKlSvZb7/92HnnnTtyzvHjxwMOeVaM\n7HPU39/fkfMdc8wxQDIub9WqVUyYMKEj5zXrlNIreVUBL3NHep9dWX02sDgLeOnzniep7p2U43Tt\nXTTMrETXXnstAMcff3zHzplVXNasWdOxc1rv6nQlb8cdd+Sggw5i9erVG2alm/WS0kNeHUen9/el\n9/sA99TY715ghqStO9Iqsy7261//Gign5LmSZ0XodCUPNv68ZD8/Zr2k60JeOsbufGBORCxIH54K\nLKux+9L0fkon2mbWrRYvXszdd9/NxIkTOeKIIzp23qy71pU8K0KnK3kAxx13HLCxEm7WS7oq5KUz\naq8E1gKnVWzyRAmzBrJfULNmzdoQvDrBlTwrUhmVvMMOO4zJkydz//3389BDD3XsvGad0DUhT9IE\nkjF2uwHHRcTjFZuXkVTzqk2t2F50e5reBgYGij6t2bCUMR4PXMmzYpVRyevv7+fNb34z4GqedZ+B\ngYGW8kg9XRHyJPUDlwMHAm+JiIVVuywkGZdXbW/g4YhYWXSbIqLpzSHPusHg4CBz5swBOh/yXMmz\nIpVRyQN32Vr3GhgYaCmP1FN6yJPUB/yQZLmUkyNifo3drgJ2kXRUxfMmAyem28zGrNtuu41ly5ax\n5557sueee3b03F5CxYpURiUPNv5xNGfOHFelraeUHvKArwOnAF8CVkk6rOKWLXR8FTAPuEzSuyUd\nlz4WwOcrDyZpUNLFVY8dLOkU4O3pQ/ukV9o4Je0mNhu1rrzySgBmz57d8XNnFZd169YxNDTU8fNb\nb8lCXqcreTNmzGD//fdnxYoVvvqF9ZRuCHnHk4S1c4Fbqm6nA0RSizwBmANcBPyc5CoZsyJicdXx\n+tj8dX0M+Cnw/7d353FSVOf+xz/PzAADiBEEFBECLgkQo15jBI2oqChGRdzwJjEuSK7G5Gf8GRPj\nDT/BkOhNvHHL/UVR40bcIiouURMiIC5oxC2ijOAKIoIgiKzCzHP/OKeHsqdngx6qu+f7fr36VTNV\np6pPPdM9/fSpc07dHZ/r5Pj7PUC3vJ+RyFbi7jzwQLiRy4gRI7b685tZbatL5gNaZHNlWoS3dkse\nbHr/TJ48eas/t0hLST3Jc/e+7l7u7mU5Hr9KlFvu7me5+/bu3tHdh7r7azmOV+buo7LWnZk4ZnnW\nz/O3xnmKtIQ5c+bw1ltv0bVrVw444IBU6qAJkSVf0mrJAzjuuDCv/oMPPqhWaSkZqSd5IrL5Mq0O\nxx57LBUV6dylUIMvJF/SbMnbe++96d27Nx999BH//GeuruEixUdJnkgRyyR5aVyqzdA0KpIvabbk\nmVnt+yjTz1Wk2CnJEylSH3zwAS+88AIdOnSovdF6GtSSJ/mSZkseqF+elB4leSJFKtPacOSRR9K+\nfXqDxNWSJ/mSZksewODBg+ncuTNVVVXMmTOn8R1ECpySPJEidffddwNwwgknNFKyZaklT/Il7Za8\nioqK2ta8e+65J5U6iOSTkjyRIjR//nyefvppKisra0cFpkUteZIvad3xIuk73/kOAHfeeWeDdxIQ\nKQZK8kSKUKYVb/jw4XTq1CnVuqglT/Il80WhsrIytToMGTKE7t27M2/ePF566aXU6iGSD0ryRIrQ\nnXfeCWxqdUiTbm0m+bJu3Tog3SSvoqKCU045BYC77rortXqI5IOSPJEi88Ybb/Dqq6/ypS99iaOO\nOirt6mgyZMmbQkjyYNOXp7vvvlsTI0tRU5InUmQmTpwIwIknnljbipYmXa6VfCmEy7UAgwYNom/f\nvixcuJAnnngi1bqIbAkleSJFZMOGDdxyyy0AjBo1qpHSW4cGXki+ZFry0v7yYmaceeaZANx0002p\n1kVkSyjJEykif/3rX1m8eDH9+/dP7V612dSSJ/lSKJdrAc444wzKysp44IEHWLp0adrVEdksSvJE\nisiNN94IwOjRozGzlGsTqCVP8qWQkrxevXoxbNgwNmzYUNtFQqTYKMkTKRIffPABjz/+OG3atOG0\n005Luzq11JIn+VIoffIyRo8eDYQvV5ozT4qRkjyRInHddddRU1PD8ccfT9euXdOuTi215Em+FEqf\nvIxjjjmGHXfckTlz5jB16tS0qyPSbEryRIrA6tWruf766wE4//zzU67NF6klT/KlkC7XQrjzxo9+\n9CMArrzyypRrI9J8SvJEisDtt9/OJ598wsCBA9l///3Trs4XaDJkyYfq6mo2btxIWVkZFRUVaVen\n1jnnnENlZSWPPvooVVVVaVdHpFmU5IkUuJqaGq6++moALrjggpRrU5cmQ5Z8yLx+2rVrVzCDigC6\ndu1a2wf2mmuuSbk2Is2jJE+kwE2ePJm5c+fSu3dvTjjhhLSrU0f79u0BWLt2bco1kWK2evVqADp0\n6JByTerKdJG49dZbWbRoUcq1EWk6JXkiBaympoaxY8cC8POf/7ygLmNldOzYEYBVq1alXBMpZpkk\nL/N6KiT9+/fnxBNPZN26dVx++eVpV0ekyZTkiRSwe++9l9mzZ9OrV6/a6RwKzTbbbANs+pAW2RyZ\n10/m9VRoxo0bh5kxYcIEPvjgg7SrI9IkSvJEClR1dTXjxo0DYMyYMQUzrUS2TMuLkjzZEpmW4EJs\nyQPYY489GDlyJJ9//jm/+c1v0q6OSJMoyRMpUBMmTKCqqoq+ffvW3kezEOlyreRDobfkAYwdO5ay\nsjJuvPFGZs+enXZ1RBqlJE+kAC1dupQxY8YAcMUVV9CmTZuUa1Q/Xa6VfCj0ljwIffPOOeccqqur\nOe+883QXDCl4SvJECtCYMWNYvnw5hx12WEGOqE3S5VrJh0IeeJE0fvx4tt9+e6ZNm8akSZPSro5I\ng5TkiRSYGTNmcMMNN1BRUcG1115bUHOG5ZJpydPlWtkSmddPIV+uBejSpUttn7zzzjuPZcuWpVwj\nkfopyRMpICtXruT000/H3bnooosYMGBA2lVqlJI8yYdiuFybMXr0aA488EA++ugjfvjDH+qyrRQs\nJXkiBeT888/nvffeY5999uGSSy5JuzpNst122wGwYsWKlGsixSzz+uncuXPKNWlceXk5t912Gx07\nduTee+/ljjvuSLtKIjkpyRMpEBMmTOCWW26hXbt2TJw4sfZ2YYWuffv2tG3blnXr1umuF7LZPvnk\nE6A4kjyAXXbZhauuugqAs88+m1deeSXlGonUpSRPpADMmDGDH//4xwDccMMNRXGZNsPM6NKlCwDL\nly9PuTZSrDKvncxrqRiMHj2a0047jTVr1nDcccexZMmStKsk8gVK8kRS9vLLLzNixAg2btzIT3/6\n09qboReTTOtLpjVGpLmKrSUPqL0DxsCBA5k/fz5HHXWUui1IQVGSJ5Ki1157jaFDh7J8+XJGjBjB\nb3/727SrtFnUkidbqhhb8gAqKyt54IEH2HXXXXnppZcYNmwYn376adrVEgGU5ImkZsqUKQwePJhl\ny5Zx9NFHc/fdd1NeXp52tTZL165dAVi8eHHKNZFilbnUmXktFZMePXowdepU+vTpw/PPP8+BBx7I\n+++/n3a1RJTkiWxtNTU1XHXVVRx11FF8+umnHH/88UyaNKlg703bFDvvvDOAbtwum8Xda187PXv2\nTLk2m6d3795MmzaNfv36MXv2bPbbbz+eeOKJtKslrVyqSZ6ZnWRmk81svpmtMbMqM7vMzLbJKtfZ\nzG4ys4/NbJWZTTGzPZr4HGZmF5vZe2a21sxeMbPCvoWAlKx33nmHoUOHcsEFF1BdXc3FF1/MpEmT\nqKysTLtqW6RXr16AkjzZPEuXLuXzzz+nc+fORTFPXn369OnDzJkzOfzww1myZAmHH344P/nJT/js\ns8/Srpq0Umm35P0U2AD8AhgGXAf8EJhicZr/uHwYOAL4MXAi0AaYZmZN+cr3a2AscG18jueAe83s\nqPyeikj9Pv74Yy688EL69+/P1KlT6datG5MnT+ayyy6jrCztt+GWy7Tk6RKVbI7M6ybzOipm2223\nHY899hjjx4+vvWvN7rvvzoQJE1i/fn3a1ZPWxt1TewDb51j3faAGGBJ/Py7+fnCizLbAMuCaRo7f\nHVgPjM1a/w/g1Qb2c5EtVVNT488//7yPGjXK27Vr54ADfuqpp/rixYvTrl5evfjiiw54v3790q6K\nFKGbb77ZAT/llFPSrkpezZo1ywcNGlT73u/evbuPHTvW33333bSrJiUm5i118plUmxDcPddN/2bF\n5U5xORxY6O5PJvZbSWjdO66RpziS0Or356z1fwa+bmZfbnalRRqwZMkS7r//fs4//3z69u3LwIED\nufnmm1m/fj1HH300s2bNYuLEiXTv3j3tqubVgAEDKC8vZ+7cuaxcuTLt6kiRmTUr/Nvfc889U65J\nfn3jG9/g2Wef5Z577mGvvfZiyZIlXHrppfTt25dvfvObXHrppUydOpU1a9akXVUpURVpVyCHg+Ny\nTlx+DZido9wbwGlm1sHd63uHfA1Y7+5v59gXYACg60vSKHdnw4YNrFq1io8//pjFixezePFiPvro\nI+bNm0dVVRVVVVUsWLDgC/vttNNOjBw5knPPPZfdd989pdq3vMrKSg444ACeeuop7rvvPs4888y0\nqyRFYv369Tz00EMAHHrooSnXJv/MjJEjR3LyySfz5JNPcv311/PII48wa9as2uS2vLyc3XbbjQED\nBtC/f3969epFjx496NGjBzvssAPbbrstnTp1oqKiED+ypZCZF9CNlWMfu5eBl939yLhuLjDL3b+b\nVXY0cAPQy90X1nO8G4Bj3H2nrPW7AXOB77t7nZsOmpkfffTRtb/nilH2uqaU2dz9dOyte2yADRs2\nsGbNmi88qqur65TL1r59ewYNGsTgwYMZOnQoBxxwQEn0uWuKW265hVGjRtG2bVv2228/OnXqlOq5\nb+5rpBT2L8Q61bd+0aJFzJ8/n/79+zN79uxW8X5Zu3Ytf//735k+fTozZszglVdeoaamptH92rVr\nR6dOnejYsSNt2rShTZs2VFRU1FmWlZVhZsSu7XlbSuF65JFHcPc6f6iCSfLiiNrpwI7Afu7+YVz/\nJvDiFiR5x7p7j6z1jSZ5W35GUmrKy8vp2LEj3bp1Y4cddqh97LLLLvTr14/+/fvTp0+fVvttu7q6\nmrPOOovbbrst7apIkenZsyeTJ09m3333TbsqqVizZg1z585lzpw5VFVV8eGHH7Jo0SIWLVrE4sWL\n+eyzz1i1alWTEkFpvQo2yTOz9sCjwNcJAyxeT2x7Dljh7sOy9vk58F/ANvVdrjWz3wLnuXv7rPX7\nEUbZHu3uj+XYzx9++OHsdbmO3+wym7ufjr1161RRUUHHjh3p0KEDHTp0oH379rRp06bOflLXggUL\nePvtt1m9enW9LTr55O71tjRs7uuoFPYvxDrlWt+2bVv23HPPVvvlqKncnbVr17Jq1SpWrVrFxo0b\n2bhxIxs2bKizrKmpqX3v5WsphW348OE5k7zU31Vm1gaYBOwDDE0meNHrhOlTsg0A3m+gP15m33Zm\ntmtWv7zM3d/fyLEPAMccc0yjdReRunr16lU7b56I5IeZ1X7pLLWBW9Jy0p4MuQy4AzgEGOHu/8xR\n7CGgp5kdlNhvW+DYuK0hjxHm4fte1vpTgdfcXYMuREREpCSl3ZL3/4GTgN8Aa81sUGLbgtjX7iFg\nJvBnM/sZsAK4mDDv0O+SBzOzjcCt7j4awN0/NrMrgYvN7DPCoI5TgCGEJFFERESkJKXaJ8/M3gV6\nA7k6b4xz91/Fcp2B/wZGAJXAs8AF7v5a1vFqCEneqMS6MkJS+APCoI4q4Ffufn8D9XL1QxAREZFi\nYGaFO/Ci0CjJExERkWJRX5JX+hMSiYiIiLRCSvJERERESpCSPBEREZESpCRPREREpAQpyRMREREp\nQUryREREREqQkjwRERGREqQkT0RERKQEKckTERERKUFK8kRERERKkJI8ERERkRKkJE9ERESkBCnJ\nExERESlBSvJERERESpCSPBEREZESpCRPREREpAQpyRMREREpQUryREREREqQkjwRERGREqQkT0RE\nRKQEKckTERERKUFK8kRERERKkJI8ERERkRKkJE9ERESkBCnJExERESlBSvJERERESpCSPBEREZES\npCRPREREpAQpyRMREREpQUryREREREqQkjzZYuPGjUu7CkVHMWsexat5FK/mU8yaR/FqvjRiZu6+\n1Z+00JmZKy5NZ2YoXs2jmDWP4tU8ilfzKWbNo3g1X0vGLB7bsterJU9ERESkBKWe5JnZzmb2BzOb\naWZrzKzGzHrnKLe3mT1uZp+Z2adm9qCZ7drE5+hqZjeb2ZL4HM+Z2RH5PxsRERGRwpB6kgfsBpwM\nLANm5CpgZrsDTwGdgO8Co4A+wAwz69bQwc2sHTAVOAL4GXA8sAB4xMwOzs8piIiIiBSW1PvkWaID\nnJmNBm4A+rj7/ESZm4AT4vqVcV1P4C3gWne/qIHjnwrcDhzi7jMS618F1rn7wBz7qE9eM6hvRvMp\nZs2jeDWP4tV8ilnzKF7N1yr75DUxmxoEzMwkeHG/hcDrhJa5xvZdk0zwoinAN82sR3PqKyIiIlIM\nUk/ymmgj8HmO9euBXcysbQP7Vsf9c+0LsMcW1k1ERESk4BRLkvcmsK+ZVWRWmFkn4GuAAZ0b2LcK\n2NbM+mWt3z8uu+SzoiIiIiKFoFiSvGuBnsD1ZraTmX0ZuAXoGLfXNLDvncBS4DYz2yOOtP1PYHAT\n9hUREREpShWNF0mfuz9jZj8CLieMrIXQp+424FTgkwb2/dTMTohl/xVXvwWMA8YDi3LtZ1an/6I0\nQPFqPsWseRSv5lG8mk8xax7Fq/m2dsyKIskDcPfr4ijb3YCV7r7QzB4DnnP36kb2fRrYNc6rV+7u\nc83sImAN8GKO8nrlioiISFErmiQPwN03AHMAzOzrwGHA95ux/9tx322AHwAT3X1tC1RVREREJFUF\nkeSZ2Unxx2/E5bfNbCmwxN1nxDnxzgWeJYyK3Rf4BXCfu9+TdayNwK3uPjqx7nJgFmHC5d0IkyKv\nBy5uubMSERERSU9BJHnAXxI/O/DH+PN04FBgA7Af8B+Eu168BVwKXJPjWGXUHVDSHbg6LpcA9wNj\n3X1FfqovIiIiUlgKYnStu5clHuWJnw+N25e4+1B37+bule6+h7tf5e51RsbG/UZlrTvL3Xu5e7u4\n/El2gmdmvcxskpmtiPfGvc/MerXsmRc+MzvJzCab2fx4398qM7ssXvJOlutsZjeZ2cdmtsrMppiZ\n5iAE4j2Xa8xsfNZ6xSzBzL5tZjMS96d+wcyGJLYrXpGZDY7nv8TMVprZi2Z2ZlaZVhmvZtwPvUnx\nMbNKM7vCzBbF4z1rZoOzyxWzpsTMzA43szvN7J1Y5i0z+6PluLVoqcesqa+xrH2uj+Um5tjWYvEq\niCQvbWbWgXB/268ApxH6+e0OTIvbWrOfElpSfwEMA64DfghMsThMKC4fJtwf+MfAiUAbQvx6plHp\nQmFm3wH2jL96Yr1ilmBmZwOTgReAEYT7Wf8F6BC3K16Rmf0bYXaBMuAswl1/XgD+ZGbnxDKtOV5N\nuR96c+LzJ2A0MAY4mjAjw9/MbK8WqX06Go0Z4Ura9sCvgSMJs10MB54zs45ZZUs9Zk2JVy0z+xbw\nPWAlic+BhJaLl7u3+gfwE8JdMXZJrOtDSG7+b9r1Szk22+dY933C/IJD4u/Hxd8PTpTZlvAGuCbt\nc0gxdp3jm/WUGJ9fJbYpZpvOuw+wFjivgTKK16bzvhxYB3TIWv8s8Gxrjxfxnuzx59ExDr2zyjQp\nPsBesdzpiXXlhEn2H0z7XLdyzLrm2G9wLHtma4pZU+KV2N4GmA1cBLwL3J61vUXjpZa8YDjh3rjv\nZFa4+3vAM4R/Bq2Wuy/LsXpWXO4Ul8OBhe7+ZGK/lYRvyq05fr8FXvOswUGRYrbJKMKXrOsbKKN4\nbVJO+AKaPTPASsIdgKAVx8vjp2Qjmhqf4YRY35MoVw3cDRxpZm3yUumUNSVm7r40x+rszwJoBTFr\n4mss42eE9+Xv2fT+TGrReCnJC75GyLSzvQEM2Mp1KQYHx+WcuGwofr1b4yVvMzuQ0OL5o3qKKGab\nHEi4deF3zextM9tgZvPM7NxEGcVrkz8R7sl9rZn1MLPtzOwHhEFqV8UyilfDmhqfrwHvuPu6HOXa\nEi7btWbZnwWgmNUys92AXwLnuvvGeoq1aLyU5AWdgeU51n9Cw/fFbXVif5VfAVPc/aW4ugv1xw9a\nWQzNrC0wAbjC3efVU0wx22QnQh/Y3wGXAUMJfc7+x8zOi2UUr8jd3yT0iToZWEiIwf8AZ7t7ZqYC\nxathTY1PY+Va7b3PLdw//mpCMjI5sUkx2+Q6wlRvmRbjXC2ALRqvQplCRYqAhRG1DwKfA8mRfM1p\num4Nfg60A37TQBnFbJMywtRIp7t75sNiupn1IcxleW1K9SpIcQToI4RLZX8gXLYdAUwws/Xufmea\n9SsSev9tATOrAO4CegDf8hwzXbR2ZnYqYe7fr6ZZDyV5wXJyf7PtQgP3xW1NzKw9ob9KH0Jn5Q8T\nm5eT+9tGl8T2ViEOo/8lYdRj+xi3jEoz+xKwCsUsaRmwK6H1LmkKMMzMdkTxShoPrACOTVwCmmZm\n2wPXmNldKF6NaWp8lgO5psbIlGt1nw9mVka4F/yhwNHunn3Zu9XHLDaIXEm4OrHBzLaLm8qBtvFz\nYHV8/7ZovHS5NngdyDV/1ABCU3SrFjt+TgL2Ab7t7q9nFXmd0K8g2wDgfXdf08JVLCS7EFrx/kx4\nc2YeABcS3tB7oJglvU7uDsnZZRSvYADwrxx9fF4gTHHRHcWrMU2Nz+tAXzOrzFHuc8LE/K3N9cBI\n4N/dfVqO7YoZdI2Py/ji58DOhNgtB74dy7ZovJTkBQ8Bg8ysb2ZFvFR0QNzWasVvbXcAhwAj3P2f\nOYo9BPQ0s4MS+20LHEvri9/LhFglH5kJfSfG399CMUu6Py6HZa0fBixw949QvJI+APbKMepuIOHS\n7TJCtwrFq35NfT09RJgCY2SiXAVhWqS/ebifeqthZr8nXKU4w93rex0pZmHqrCHU/RxYTLhCcQhh\n9g5o6XilPd9MITwIE67OA/5FGM48HHiV8GHcIe36pRyb6whz+IwHBmU9esYyFl+w8+ML80jCLemW\nZsq09gd158lTzL4YnyfiuZ9NmKD2xhiz0xSvOrHKzPH2ePxfdQRh4EUN8N+KlwOcFB+Z/1/nxN8P\nam58CH3PPiEkN4cRrmqsAfZO+zy3cswuiutvInyhSH4W7JJ1rJKPWWPxqmef98iaJ6+l45V6oArl\nAfSKgf2UMN/U/dQzuWFrehAmb6yOL+LsxyWJcp0JUzssA1YTvq18Pe36F8qDrCRPMasTn04xUfkI\nWA+8QrgcpHjljtdQwl16lsT/Vy/FD5kyxav2/ZZ5JP9/TW1ufIBKwhxniwgtpTMb+iAv1kdjMQOm\nNfBZcHNri1lTXmM59qkzGXJLx8viE4iIiIhICVGfPBEREZESpCRPREREpAQpyRMREREpQUryRERE\nREqQkjwRERGREqQkT0RERKQEKckTERERKUFK8kRERERKkJI8EWkRZnaImdWY2dgtPM4Z8Tin56tu\nDTxXn/hct7T0c5UaM7s1xi7zuCjtOjWVmV2YVXf9/aUkKMkTKRGJD6hqM9ulgXLTEmVbPHEC8nVb\nnSYdJ5EUJh+fmdkCM5tiZpea2W75eK4cz52XxLbIXQ2MA55KuR7N8QyhztfE33UrKCkJFWlXQETy\naiPhfX0W8MvsjWa2O3Bwolwpf5i9AkyOP7cHuhNupv7/gF+a2R+AC929OrHPB0A/wj2st0Qpx7Ux\nV7v7/LQr0RzuPhOYaWZfBn6Sdn1E8kVJnkhpWUy4yfWZZnZJVgIDMDouHwaO36o12/pecfdfZa80\ns0OAWwkf5pXADzPb3H0jMDcPz215OIZsffq7SUnR5VqR0uLAjcCOwDHJDWbWBjiDcGnqjfoOYGa7\nm9ntZrbQzNbH5W31XeI0sx3M7E9mttjM1pjZy2Z2WkOVNLMuZna5mc2J+6wws3+Y2dBmnm+zNJ7x\nIQAACKtJREFUuft04Ejgc+A/zGzvRL1y9smL5/jfZvamma0ys+VmVmVmt5hZ31jmVmBq3GVs1uXi\ng2KZbc3sZ2Y21cw+iPFdYmYPmtmgXPWN+08zs+3N7AYzW2Rm68xstpmdUd95mtkRZvZwPP46M5tv\nZpPN7LAcZY80s0fNbGks+5aZ/c7MvtSc2DbEzKbHc6kws0vM7G0zWxvj+INEuXPN7LX4ulhgZuPM\nzLKOVft3MrNdzWySmS0zs5Vm9ncz2yOW62ZmN8WYrTWzF2KSL9IqqCVPpPTcBVxJaLV7MLF+ONAN\n+BnwlVw7mtk3gX8A28R93wD6A6cCx5nZ4e4+K1G+K/As0JfQB+tpYCfgemBKPc/xZWA68GVgBvBo\nfL5jgMfN7Gx3v2kzzrvJ3P1NM/sL4by+S7i0+4Uiifp2ICTGuwB/J8TFgD6EmN4LvAs8EPc7nXB+\n0xPHey8uBwC/Bp4ktKYuJ8RhOHCUmR3r7n/LUeXtYh3WA38B2gEjgZvNrMbdb08WNrNLCZelPyNc\nsl4A9AQOAL4HPJEoOxYYCyyLdVoC7AVcCHzbzPZ3989y1Glz3QPsB/wV2ACcDEwws+r4vN+L9ZgC\nHAdcAqwBfpfjWH2A5wiv05sJr8PjgelmdiDhtbWc8J7YHvh34DEz+4q7L8jjOYkUJnfXQw89SuAB\n1ADz4883Ej5Aeya2P074wKskJBo1wGmJ7QbMAaqB72Qde2QsPwewxPob4vrfZ5X/BqGlrAa4JGvb\ndEKfwJFZ678EvEz4QO+eWH9Gdl0biUOm/M2NlBsVy01LrOuTvS9wbK5zjNsqgG0Svx+S65wT27cF\nuuRY3xNYCLxRz9+1JsY6Gfv+8W/8elb5I2L5t4AeuZ4r8fOQWPZpYNuscqfHbVc2Me63xvK969k+\nPW5/PvlchMRsPaEf5NvJOsfXxMeExLM8x9+pBrg463nGxPUrgD9mbTu1oXPK9ffXQ49ifuhyrUhp\nuhEoJyQymdazocAd7r6unn0OAL4KzHT3u5Ib3P0vhETgq8CB8ZhtCK0uKwkjE5PlXwTuyH4CM9sL\nOAi4Lx4zuc+n8TiVwIlNPtPN92Fcdmti+Tpxc/eN7r6qqU/o7ivd/ZMc6xcC9wH9zGznHLuuBi5w\nd0/sM4fQitovtjZm/J+4/Km7L6rnuTLOi8sfuPvKrHK3Aa8S/sb59Ivkc7n7u4RWyk7A+GSd42vi\nEaAroYU427vAf2Wtuy0uywmt1kl3Er5g7LUlJyBSLHS5VqQEufs/zew1YJSZ/Zpw6dYIyV999onL\nqfVsn0ZI8PYmXJrtRxi1+oLnvpz3JKE1KGn/uNzOzMbl2CeTcPVvoJ75kunn1dhI2OmEVrZfmNk+\nwGOEhPcVd69p9pOafYsw6GN/wvm2zSrSkzDKN2lePcnkAsJ5dCa0gEIYQVxDaLltzP6E1sCR2f3e\norZANzPr7O7Lm3C8xjgwK8f6TML9Yo5tmaR0Z8L5Jr2STHyjTJI4191Xf+HJ3WvMbEk8lkjJU5In\nUrpuBK4FjgLOBGa5+6sNlM90sq/T+pO1frus8ovrKf9RjnXbx+XQ+MjFgY71bMunTMvQxw0VcvfP\n4qCISwl9546Mm5aa2R+BX3sYldsoMzsemERIyKYQLk+uJiRlQwjT27TLseuKeg6Zed7yxLrtgOXu\nvr4JVdo+7tvQvH5O6DOZjySPer4QZM4j19Q1mW1tcmyrU97dN8Z8tb5pcDbWcyyRkqMkT6R0TQR+\nC0wgJDTjGimf+VDcsZ7tPbLKZZY71FM+13Ey+5zn7v/TSH1a2pC4fL6xgvES52gAMxsAHAr8iDAo\noCwum2I84bLvvu7+ZnKDmfUkJHlbagXQ2cwqG7g0n/EpgLt3zcPzikiBUZ88kRIV+zNNIlz+W0UY\nYdiQl+JySD3bh2SVqwLWAnub2bY5yh+SY93MuDyokbq0KDPrRxjVWUPop9Vk7v5GTFAzLZHHJTZn\n5iUsJ7fdCIMrshO8MmJfxzyYSfjfPqyJZbvExFVESoySPJHSNgYYARyZ3T8pm7s/A7wJHGhmXxj4\nYGYnEZKQN9396Vh+A/BnwojRcVnl9yVHh/04IOMp4AQzOzNXPczs62bW1MEQzWZmBxP6q7UBrnP3\n1xopP8DMcrVWZloq1yTWLYvLL9dzuHeBr5hZplWU2BduHKEfYj7ulPGHuPy9mdUZrJC17qq4vDFZ\np0TZjmY2MA91EpEU6HKtSAnzMBdYc+YDO53QV+weM3uQkPR9lZAorgSyJzn+T+Aw4PyY2D1DuKw7\nkjAP2vAcz/FdwuCOP5nZecA/CZcYdwb2BL5GGDzQYF+5Jvi3xOCOdoTLygMJyVQ18Hvg5004zhHA\nFWb2LDCPMJ3HzoQWvGrgikTZKsJAgX83sw3AfELidruHW31dRZhD8GUzu58w6OFbsU4PE6Zr2SLu\nPiUOthkDzDGzyYSBHDsQEvWZhD6auPtUM/sFcDkwz8weJczptw0hUT2IkJR/e0vrFemOEiJbkZI8\nkdbJydFqFEflfpOQIBxOSDo+JkyHMt7d52WVXxZHi14Wy+5LSHTOAd4nR5Ln7gvN7BuEqT5OJCR9\n5YSBHW8QbhI/u7G6NnJuEBLGzFQZa4BPYt3+Akx093eaeLzHgV6EhGc4oeXyQ+BvhPnWnkucW00c\nXPFfhMvBneKmGYQ5DG8ws/XA+YSEeQ0hiTodOImsu5Q08Vxz/R0vMbOZhClSjiEMZFlMGNl6W1bZ\n35nZM7HsgYTkdQUhWZ1AMy9nN7euW7htc+sh0ipY3dHnIiIizWPhtm6nAX3d/f2Uq7NZzKwP8A5w\nq7uPSrc2IltOffJERCQfMi0G78b7yl6Uam2awcwuNLMaQoInUjJ0uVZERPJhMmFgScZTaVVkMzzD\nFwcPZd/LWKQo6XKtiIiISAnS5VoRERGREqQkT0RERKQEKckTERERKUFK8kRERERKkJI8ERERkRL0\nv7MhZlLVt+23AAAAAElFTkSuQmCC\n", | |
| "text/plain": [ | |
| "<matplotlib.figure.Figure at 0x10dbee110>" | |
| ] | |
| }, | |
| "metadata": {}, | |
| "output_type": "display_data" | |
| } | |
| ], | |
| "source": [ | |
| "# Setup figure and axes\n", | |
| "# Generally plots is ~1.33x width to height (10,7.5 or 12,9)\n", | |
| "fig = plt.figure(figsize=(10,7.5))\n", | |
| "ax1 = plt.subplot(111)\n", | |
| "\n", | |
| "# Set labels and tick sizes\n", | |
| "ax1.set_xlabel(r'Model Distance [mm]', fontsize=20)\n", | |
| "ax1.set_ylabel(r'Temperature [C]', fontsize=20)\n", | |
| "\n", | |
| "# Plotting\n", | |
| "ax1.plot(X*1000., phi.value, color='k')\n", | |
| "\n", | |
| "plt.title('Temperature Final Result', fontsize=24)\n", | |
| "\n", | |
| "# Set limits\n", | |
| "ax1.set_xlim(0,np.max(X.value*1000.))\n", | |
| "#ax1.set_ylim(19.8,21.2)" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": null, | |
| "metadata": { | |
| "collapsed": true | |
| }, | |
| "outputs": [], | |
| "source": [] | |
| } | |
| ], | |
| "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.11" | |
| } | |
| }, | |
| "nbformat": 4, | |
| "nbformat_minor": 0 | |
| } |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment