{ "metadata": { "name": "Compute_seismic_velocities_beginner" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "heading", "level": 1, "metadata": {}, "source": [ "BurnMan - Compute seismic velocities - simple case" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Copyright (C) 2012, 2013, Heister, T., Unterborn, C., Rose, I. and Cottaar, S.\n", "Released under GPL v2 or later.\n", "\n", "This example script is intended for absolute beginners to BurnMan.\n", "We cover importing BurnMan modules, creating a composite material,\n", "and calculating its seismic properties at lower mantle pressures and\n", "temperatures. Afterwards, we plot it against a 1D seismic model\n", "for visual comparison. Run the seperate cells of the code with shift+enter\n", "\n", "requires:\n", "\n", "- geotherms\n", "\n", "- seismic models\n", "\n", "- compute seismic velocities\n", "\n", "- composite mineral helpers\n", "\n", "teaches:\n", "\n", "- changing composition\n", "- changing geotherm\n", "- computing self-consistent depth to compare to seismology" ] }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "1. Import modules" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Here we import standard python modules that are required for \n", "usage of BurnMan. In particular, numpy is used for handling\n", "numerical arrays and mathematical operations on them, and \n", "matplotlib is used for generating plots of results of calculations" ] }, { "cell_type": "code", "collapsed": false, "input": [ "import os, sys, numpy as np, matplotlib.pyplot as plt\n", "sys.path.insert(1,os.path.abspath('..'))" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Here we import the relevant modules from BurnMan. The burnman\n", "module imports several of the most important functionalities of\n", "the library, including the ability to make composites, and compute\n", "thermoelastic properties of them. The minerals module includes\n", "the mineral physical parameters for the predefined minerals in\n", "BurnMan" ] }, { "cell_type": "code", "collapsed": false, "input": [ "import burnman\n", "from burnman import minerals" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "2. Import seismic model" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Here we create and load the PREM seismic velocity model, which will be\n", "used for comparison with the seismic velocities of the \"rock\" composite" ] }, { "cell_type": "code", "collapsed": false, "input": [ "seismic_model = burnman.seismic.PREM()" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 3 }, { "cell_type": "markdown", "metadata": {}, "source": [ "We create an array of 20 depths at which we want to evaluate PREM, and then\n", "query the seismic model for the pressure, density, P wave speed, S wave\n", "speed, and bulk sound velocity at those depths" ] }, { "cell_type": "code", "collapsed": false, "input": [ "depths = np.linspace(750e3, 2800e3, 20)\n", "pressure, seis_rho, seis_vp, seis_vs, seis_vphi = seismic_model.evaluate_all_at(depths)\n", "\n", "plt.plot(pressure/1.e9,seis_vs/1.e3,'k',label='Vs')\n", "plt.plot(pressure/1.e9,seis_vp/1.e3,'b',label='Vp')\n", "plt.plot(pressure/1.e9,seis_rho/1.e3,'r',label='Vphi')\n", "plt.plot(pressure/1.e9,seis_vphi/1.e3,'g',label='rho')\n", "plt.xlabel('pressure (GPa)')\n", "plt.ylabel('velocity (km/s) density (kg/m^3)')\n", "plt.xlim(min(pressure)/1.e9,max(pressure)/1.e9)\n", "plt.title('PREM')\n", "plt.legend();" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAEWCAYAAACHVDePAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlYU1f+P/B3UFZZQ8KugqBsgiAIomJR6loQ67hUbXGt\n09axTjtT6886U6vWpdaOXew8/dZxwaVudUFFx6pFXKuiAqKCIAjKDmHf4f7+uJNLQtiTkEA+r+c5\nT25uQvLJxb45PTn3XB7DMAwIIYT0alqqLoAQQojyUdgTQogGoLAnhBANQGFPCCEagMKeEEI0AIU9\nIYRoAAp7onHs7e1hYGAAIyMjWFlZYdGiRaioqEBQUBD09fVhZGQEgUCAsLAwvHz5kvu5devWQVtb\nG0ZGRlzj8/nc41paWrC0tERDQwO3r66uDhYWFtDSov/UiGrRv0CicXg8Hs6ePYuysjLcv38f9+7d\nw8aNG8Hj8bBz506UlZUhNTUV1dXV+Pjjj6V+bu7cuSgrK+NaUVGR1Gvz+XycP3+eu3/+/Hnw+Xzw\neLxu+3yEtITCnmg0GxsbTJkyBY8ePZLab2JigrCwMCQmJnL7GIZBe+cgvvPOO4iIiODuR0REIDw8\nvN2fI0TZKOyJRhKHb2ZmJqKiojB8+HCp/YWFhThx4gT8/f079bphYWGIiYlBaWkpRCIRrl+/jrCw\nMMUWT0gXUNgTjcMwDKZPnw4zMzMEBgYiKCgI/+///T8wDIMPP/wQpqamEAqFKC8vx86dO6V+9ujR\nozAzM+NacHCw1ON6enoIDQ3F4cOHceTIEYSFhUFPT687Px4hLaKwJxqHx+Ph9OnTEIlESE9Pxw8/\n/AA9PT3weDx8//33KC4uRnx8PF68eIGoqCipn50zZw5EIhHXLl++LPPa4eHh2LdvH/bv309DOERt\nUNgTIkEczEOHDsWGDRuwevVqbh+Px+tQcAcGBiInJwd5eXkYPXq0UuslpKMo7AlpxYIFC1BZWYmj\nR48CQKd66GfOnEFkZKSySiOk0yjsCZEgOUVSW1sbK1euxFdffcU9duTIEal59sbGxigoKJD5WTc3\nN7i6urb4uoSoAo/WsyeEkN5PaT37xYsXw9LSEh4eHjKPbd++HVpaWjInpBBCCFEOpYX9okWLcOHC\nBZn9mZmZ+O233zBw4EBlvTUhhJBmlBb2gYGBMDMzk9n/8ccfc2OghBBCuke3fkF7+vRp2NnZwdPT\nszvflhBCNF7f7nqjyspKbNq0Cb/99hu3r7XvhmnmAiGEdE1rudptPfvU1FSkp6dj2LBhcHBwwMuX\nL+Hj44O8vLwWny9edEpd2+eff67yGlTd6BjQMaBjoF7HoC3d1rP38PBAbm4ud9/BwQGxsbFS64ET\nQghRDqX17OfOnYtRo0YhOTkZ/fv3x549e6Qep6EaQgjpPkrr2f/yyy9tPv78+XNlvXW3CAoKUnUJ\nKkfHgI4BQMcA6BnHQC3PoO3oglOEEEIAhgEqKgAjo9azs9vG7AkhpCv4fD5EIpGqy1Ar+vpmGDOm\nCPn5QEEBkJ8P9OnT9s9Qz54QotYoD2TxeDycP89AKASEQkAgAAwM2j5WFPaEELVGeSCrtWPS1rGi\nJY4JIUQDUNgTQogGoLAnhBANQGFPCCFdNHnyZHz++ecy+0+fPg1ra2s0NjaqoKqWUdgTQkgXLVy4\nEAcOHJDZv3//frz99tvQ0lKfiFWfSgghpIcJCwtDYWEhrl27xu0TiUQ4d+4cwsPDERUVBXd3dxgb\nG8POzg7bt29XWa10UhUhRK3U1QHp6cCzZ0BysqqraZu+vj5mz56NiIgIBAYGAgCOHj0KV1dXeHh4\nYOLEiTh+/DhGjx6NkpISlS4TQ2FPCOl2jY3Aq1dNgS7ZMjIAGxtg8GBgyJCOvZ4iFlbs6lz+BQsW\nICQkBDt37oSOjg4iIiKwYMECAICOjg4SExPh4eEBExMTeHt7y11nV9FJVYQQpWAYoLCwKcQlgz0l\nBTAxYcN8yJCmYB8yBBg0CNDVbXqdnpAHgwcPxsaNG+Hr6wtXV1e8evUKQqEQ9+7dw8aNGxETEwNP\nT09s2bIFI0eOlPv9unJSFYU9IUQu5eWyPXTxfYZpCnHJYB88GDAy6tjr94Q82LBhA27fvg0/Pz/E\nxsYiMjJS6vGGhgZ8//33+Oabb5CRkSH3+3Ul7GkYhxDSrvp6IC2tKcyTkppuRSLAyakpzIODgfff\nZ7fNzQFNuHRFeHg4NmzYgPj4eOzYsQMAUFdXh6NHjyIkJAQmJiYwMjJCn/ZWK1Mi6tkTQgCwvfDc\nXOkwF2+npwPW1oCzMxvikrd2doAyZxj2lDwYN24c4uPjkZOTA21tbdTV1WHatGn4448/0NDQABcX\nF/zrX//CqFGj5H4vGsYhhLSrsrIpxMVNHOw6Ok09dMlAd3QE9PRUUy/lgSwKe0IIAHa2y8uX0oGe\nlAQ8fcqufe7kJB3m4m11vCQ05YEsCntCNEx5eVMv/enTplB/9gwwNm4KcheXpu2BA9u/0IU6oTyQ\nRWFPSC/EMEBWFvDkCRvokq2oiO2lS4a5uBkbq7pyxaA8kEVhT0gPVlvLzj8XB7lkuPfrxwa6q2tT\nsLu4AAMGKPfLUXVAeSCLwp6QHqC4uOVAf/GCDW/JUBc3MzNVV606lAeyKOwJURPiaYxPngCPH7Pt\nyRO2lZXJBrqrKzvjRfLMUcKiPJBFYU9IN2MYIDOzKcwlg53HA9zcmpqrK9vs7DTjRCNFoTyQRWFP\niJI0NLBnkDYP9adP2S9CXV2bAl18KxRSqCsC5YEsCntC5NTYyIZ6YmJTe/SInd5oaSkb6K6ugKmp\nqqvu3SgPZCk87Ovq6nDx4kXExMQgPT0dPB4PAwcOxNixYzFp0iT07aucpXXol0uUrbGRXUq3eag/\nfQoIBIC7O9uGDmVvXV3ZGTGk+6lzHkyePBn+/v744osvpPafPn0a7733Hl69eqWUq1UpNOw3bNiA\nX3/9FQEBAfDz84ONjQ0aGxuRnZ2NO3fu4Pbt25g5cybWrl2r2E/RTsGEdAbDsGeSSoZ6YiI7BGNi\n0hTq4mB3c+v4aoyke6hzHhw+fBifffYZUlNTpfbPnDkTDg4O2LZtm1LeV6FhHxkZidDQ0FYvCtDY\n2IizZ89i2rRpcpTcMnX+5RL1VVwMJCRIt0eP2DVdxD10yUbDLz2DOudBVVUVrK2tcebMGe5KVSKR\nCDY2Nrhz5w62b98OPT09PH/+HLdv38bw4cMRERGBAQMGyPW+NGZPNEJNDTvc0jzYi4vZEPfwkG4C\ngaorJvJQ9zxYtmwZGIbBzz//DAD46aef8NNPP+H+/ftYuHAhTpw4gaioKPj5+WHVqlWIjY2VumZt\nVyg07EtKSvD9999DIBBg6dKl+PLLL3H37l14e3tjzZo10NfXl6vYtqj7L5d0j8ZG9kSj5qH+/Dng\n4CAb6vb2vf9sUk3UkTxQxKynrkbOjRs3EBISgtzcXOjo6GD06NGYPXs2Vq5ciYULF6K2thaHDh0C\nAFRUVMDExAQvXryAra1tl2tV6MVLFi5ciEGDBiErKwtBQUHw9PTEqlWrEBkZieXLl2P37t1dLpSQ\n5ior2SGXhw+BuDi2JSSw4+fiMH/jDWD1avYkJFUtt0vUkyr7hqNHj4ZAIMDJkyfh6+uLu3fv4tSp\nUwDY8LWzs+Oe269fP/D5fGRlZckV9l3Ratinpqbi5MmTYBgG1tbWiImJgZaWFgIDAzFs2LDurJH0\nIgzDXmhaHOjilpHBrvcybBjbZs4EPD3ZKx0Rou7Cw8MRERGBp0+fYvLkyRAKhQDYi5hnZmZyzysv\nL0dRURFsbGy6vcZWw148XYjH42HKlClS9wnpiNpadtZL82DX0moK9dBQYO1atreura3qignpmpYu\nSygWFRWFGzduYMSIEfjHP/6BgICAbu/VA22EvY+PD8rKymBkZIQ9e/Zw+1NSUmDcW9ZOJQpTVsYG\n+f37wIMH7G1yMju27uXFBvsnn7C3VlZ0ZinpXQYOHIjRo0cjPj5eaoYij8fDvHnz8MUXX+DWrVvw\n8fHBgQMHVFJjl2bjMAyj1B4+fUGr3goK2EAXh/qDB+xc9qFDAW9vYPhw9nboUECJ3+MTDdGT82DR\nokWws7PDhg0bFPq6Cv2CVqy+vl7mTNmOBv3ixYtx7tw5WFhYICEhAQDwySef4OzZs9DR0YGjoyP2\n7NkDExOTDr0e6V7ii2ZI9tbv3wdKStje+vDhwNSpTcMwSjqhmpAeS53+SLU5Ua2srAxhYWFdfvFF\nixbhwoULUvsmTpyIxMRExMXFYciQIdi8eXOXX58ojnj1xpMngc8+AyZPZteC8fYGdu4EqqqA+fOB\nK1cAkQi4ehX417+A8HC2B09BT4gsHo+nNt9ztvqfaHZ2NsLCwvDZZ591+cUDAwORnp4utW/ChAnc\ntr+/P3799dcuvz7putxc4O5d4N69psYwwIgRgK8v8MEHgI8PYGND4+uEdJXk952q1mrYjxkzBtu2\nbZOrZ9+e3bt3Y+7cuS0+tm7dOm47KCgIQUFBSqujtyssBGJjpcO9ooIN9REjgCVLgH//m9ZZJ6Sn\niY6ORnR0dIee2+oXtL6+vliwYAFWrFghVzHp6ekIDQ3lxuzFvvzyS9y/f7/Fnn1P/kJG1crKpHvr\nd++yYT98eFO4+/qys2Qo2ElPQHkgS6Ff0EZHR2P27NnQ0tLC8uXLFVclgL179yIqKgqXL19W6Otq\nmoYGdh777dvAH3+wLS2Nnd44YgQ7h339emDwYFpGgBBN12rYGxoaIjIyEu+9955C3/DChQvYtm0b\nrl69Cj06571TsrPZQBeHe2wsYG0N+PsDI0cC77/PnnVKJycRQppT6qqXc+fOxdWrV1FQUABLS0t8\n8cUX2Lx5M2pra8Hn8wEAAQEB+PHHH6WLov9tQ2UlO81R3GO/fZsdZ/f3b2p+fsD/DiMhvRblgSyl\nLnFcWlqK+vp67j5fiSmjib/cjAzgxg3g+nU22J88YZfrFffa/f0BJycaZyeapzfngZaWFlJSUjBo\n0CCZxw4ePIiIiAj897//lXlMKWH/008/4fPPP4eurq7U+jjPnz/v0Ifpit78ywXYsfaEBDbYb9xg\nW3U1MGYMMHo0EBDAzm+ns08JUe88kPeyhG2FfVuUEvZOTk64ffs2BN14BQh1/uV2RXk5OxQj7rn/\n8Qc7f330aLaNGUO9dkJao855IO9lCbsz7NudozFo0CClXqikN8rKAo4dA/76V3aao6Ul8M9/smPu\nf/kLkJrKDtPs2gUsWsTOlqGgJ6TnCQsLQ2FhodSVp0QiEc6dO4fw8HAsXLgQ7733HiZOnAhjY2ME\nBQUhIyND6jV+++03DBkyBGZmZvjLX/7C7d+7dy93qUNFaPck9y1btiAgIAABAQHQ0dEBwP71+O67\n7xRWRE/GMGx4x8Q0tZKSpl77t9+yZ6LSxCNCeh99fX3Mnj0bERERXDAfPXoUrq6u8PDwAAAcOnRI\n6rKE8+fPl/rjcO7cOdy7dw8lJSXw8fFBaGgoJk2apPBa2w37ZcuW4fXXX4eHhwe0tLSUvuKlumts\nBBITgWvXmsJdSwsYO5Ztn3wCuLrSvHZCupUKr0u4YMEChISEYOfOndDR0UFERAQWLFjAPR4SEoIx\nY8YAYE8mNTExwatXr7g17VevXg1jY2MYGxtj3LhxePjwoWrCvqGhAd98843C37inqK9nV3wUB/v1\n6+x0x7FjgSlTgM2b2WufavDfP0JUT4Vj+vJeltDKyop73MDAABUVFUqps92wnzJlCn766SdMmzYN\nurq63H5lTr1UpepqdokBcbjfugUMHMiG+7x57BoyKriiGCFEjfXoyxKKHTp0CDweD1u2bOH2KXvq\nZXeqr2fPRL18mW1//MEOw4wdy678eOgQXQeVENI2RV2WkGEYpc08ajXss7KyYGNjI7NEcU/HMOyY\n+5UrbLhfvQr07w8EB7OzZ8aOBehaKoSQzujqZQmbf/8puf69otfCb3We/ZQpU1BUVIRx48Zh8uTJ\nGDNmjMwVq5RF0fNq09Obeu5XrrAnKwUHs238eHZqJCFEPanzPPv29IjLEp4/fx5VVVWIjo7GiRMn\n8Pe//x39+/fHlClTMHnyZAwYMEBxlStYXh7w++9NAV9ezoZ6cDDw5Zfs8r6EEKJs6vRHqs2uur6+\nPqZMmYIpU6YAAJ4/f47z589j+fLlyMnJwd27d7ulyPYwDPDoEXtJvVOngOfP2eGY4GBgxQr2snk0\nW4YQ0t3U6bKEXV71sra2ljvJStE68r9tjY3sTJlTp9iQr68H3nwTmD6dPZmJrolKSO/Qk4dxlEWh\nwzhiRkZGMvtMTU3h6+uL7du3d3pNB3nU1rJj7idPAqdPAxYWbLgfOwZ4eVHvnRBCWtNu2K9cuRL9\n+/fnrhV7+PBhpKamwtvbG4sXL+7w9Q+7qqwMOH+eDfgLFwA3N7YHf/06u3gYIYSQ9rU7jOPp6Yn4\n+HipfV5eXnj48CGGDRuGuLg4xRfF42HXLgYnT7InNo0ezQb8tGmAxMlmhBANQMM4spQyjGNgYIAj\nR45g1qxZAIDjx49zlxNU5hcPFy8C8+cDBw/SvHdCCJFXuz371NRUrFy5Erdv3wYAjBw5Ejt27ICt\nrS1iY2O5BX4UWhT9JSeE/A/lgSylXLykqKhIZh2ctLQ0OChxsjr9cgkhYj0lD9LT0zFo0CDU19e3\neXUqRVDKxUtCQkJQUlLC3X/8+DFCQkLkKJMQQoiitfcHsd0x+88++wyhoaGIiopCUlISwsPDcfDg\nQYUVSAghvUF9fX23vt/cX+civyIfBZUFyK9kb9vSbti/8cYbqK2txYQJE1BeXo4TJ07A2dlZYQUT\nQkhPZW9vjw8++AAHDx5EUlISAODAgQP4xz/+gcrKSnz00UdYs2YNAKCmpgaffvopjh07BgCYPXs2\ntm7d2uWTU0OHhEJgIIDQQAiBgQACAwEM/mHQ6vNbHbNfsWKF1P0rV67A0dERAwcOVPplCXvKGB0h\nRPnUOQ/s7e3B5/Nx5swZlJWVwc3NDe+++y6+++47JCUlwc/PD3FxcXB2dsY///lPXLp0CZGRkQDY\n69cGBwdj/fr1nX5fhU699PHxkZpaKb6v6ZclJISoH94X8mcS83nn/6DweDx8+OGHsLW15ZaD//zz\nz6GrqwtPT0/uXCRnZ2ccOnQIP/zwAwQCAfe8P//5z10K+65oNewXLlzYLQUQQoi8uhLUitK/f3+p\n+80vM1heXg6AvUbIwIEDuccGDBiArKys7ikSbczGeeONN3Ds2DFUVlbKPFZZWYkjR45g6tSpSi2O\nEELUXUdHOppfDCojI6NbL0/Yas9+z549+OGHH/D555+jT58+sLa2BsMwyMnJQX19PebMmYN9+/Z1\nW6GEENKTzZ07Fxs3bsSIESMAAOvXr8c777zTbe/fathbWFhg/fr1WL9+PXJycvDixQsA7OW3rGiB\nGkIIkdFWL3/t2rUoLS2Fp6cnAHY2ztq1a7urtK6vZ69M6vztOyGk+zQyjeij1YfyoBmlLIRGCCHK\nVFVXhbTiNKQWpSJVlIrnoufcbXpxuqrL6zWoZ08IUSqGYVBYVciG+P8CPVWUitQiNtALKgsw0HQg\nBpkNgqOZI9v4jhhkNggOpg4w1DWkPGhGKQuhRUZGIiQkROkL+0iisCekZ2lobEBmaSYX4JKBnipK\nBQ88OPKbgtzRzJELdztjO/TR6tPqa1MeyFJK2M+fPx+3bt3CzJkzsXjxYri4uCim2jbQL5cQ9VNd\nX400URpSilKkwjylKAUZJRkQ9hM29cpNB0mFO1+f3/4btILP50MkEinwk/R8ZmZmKCoqktkvV9gD\nQElJCX755Rfs3bsXPB4PixYtwty5c1u8Pq0iUNgTohrF1cVNQy3/uxWHe35FPgaYDIAT34kLcie+\nExzNHOFg5gC9vnqqLl/jyR32AFBQUID9+/djx44dcHNzw7Nnz/Dhhx/iww8/VGixAIU9IcrCMAzy\nK/ORUpQi1cThXl1fDUd+U4hzgc53RH/j/m0OtxDVkyvsT58+jb179+LZs2cIDw/HwoULYWFhgcrK\nSri5uUmdESZp8eLFOHfuHCwsLJCQkACAvRDKnDlz8OLFC9jb2+Po0aMwNTXtVMGEkLYxDIPs8myZ\nQBc37T7acOI7cYEueWvRz4LWvurB5Ar7BQsWYMmSJRg7dqzMY5cuXcLrr7/e4s9du3YNhoaGCA8P\n58J+1apVEAgEWLVqFbZu3QqRSIQtW7Z0qmBCCDv//FXpKy7AnxU9k+qlG+oYcoHuZObUFO5yjp8T\n9SZX2K9atQpfffWV1L5PP/0UW7dubfeN09PTERoayoW9i4sLrl69CktLS+Tk5CAoKAhPnz7tVMGE\naArJQH9W9Ixtheztc9Fz8PX5rQa6sa6xqssnKiDXSVW//fabzL6oqKgOhX1zubm5sLS0BABYWloi\nNze3069BSG8iHnIRh7hkoKcWpcJEzwSD+YMx2HwwBvMHY6TnSAw2HwxHM0f00+mn6vJJD9Jq2P/7\n3//Gjz/+iNTUVHh4eHD7y8rKMHr0aLnfmMfjtTk2uG7dOm47KCgIQUFBcr8nIarAMAwKKguQXJiM\n5MJkqVBPKUqBgbYBF+aD+YPx1tC3MJg/GE58JxjpKmfGG+kdoqOjER0d3aHntjqMU1JSApFIhNWr\nV2Pr1q3c/xoYGRnB3Ny8Qy/e0jBOdHQ0rKyskJ2djXHjxtEwDuk1Kmor8KzoGRfqSYVJ3DYAOJs7\nS4W6eNtEz0TFlZPeokvDODweD/b29ti5c6dMD7yoqAh8fue/5Jk2bRr27duHTz/9FPv27cP06dM7\n/RqEqFJdQx3Si9O5EE8u+l+wFyShsKoQTnwnDDEfgiHmQzDefjze83kPQ8yHQGAgoFkuRKVa7dm/\n8cYbOHfuHOzt7Vv8R5qWltbmC8+dOxdXr15FQUEBLC0tsX79eoSFhWH27NnIyMigqZdErRVUFuBp\nwVMkFSThaSF7m1yYjPTidNgY2XCBLtloHjpRNYWcVNWdKOxJd6hrqMNz0XMkFSaxwS6+LUhCfWM9\nXAQucBG4wNncGc4CZzibO8OR70hnihK1JVfY37hxA8OGDYOhoSH279+PBw8eYOXKlVLXUlQ0Cnui\nSEVVRWwPXTLQC5OQJkqDrbEtF+iSt3RyEemJ5Ap7Dw8PxMXFISEhAQsXLsSSJUtw7NgxXL16VSnF\nAhT2pPMYhkFWWRaeFDzBk/wneFzwGE/yn+BJwRNU1VXBWcCGuIu5C7ftxHeiXjrpVeSaZ9+3b19o\naWnh1KlTWL58OZYuXYrdu3crvEhCOqKhsQFpxWlckD8peILH+Y/xtOAp9PrqwVXgClehK9wEbviT\n65/gKnCFjZEN9dKJxms37I2MjLBp0yYcOHAA165dQ0NDA+rq6rqjNqLBahtqkVyYzPbS8x9zwf6s\n8BmE/YRwE7rBVeCKUXajsMR7CVwFrjA36NiUYEI0UbvDONnZ2Th06BD8/PwQGBiIjIwM/P7771iw\nYIHyiqJhHI1R31iPlKIUJOYl4lHeIyTmJyIxPxHPRc8xwGQAF+riW2eBMwx1DFVdNiFqiWbjEJVr\naGzAc9FzNszz2EB/lPcIz4qewdbIFu4W7nAX/q9ZuMNF4ELj6YR0klxh/+uvv2L16tXIzc3lXoTH\n46G0tFTxlYqLorDvsRiGQWZpJuJz49neev4jJOYlIqkwCUIDIdwt3DHUYigX7K5CVxhoG6i6bEJ6\nBbnC3tHREWfPnoWrq6tSimsJhX3PUFJdgkd5j5CQl4D43Hgk5CUgITcBBtoG8LD0wFCLoRgqHAp3\nC3e4ClxpnRdClEyu2ThWVlbdGvRE/dQ31iO5MJkL9PjceCTkJqCgsgDuFu7wsPCAh4UHZrnNgoel\nBwQGAlWXTAhppt2e/cqVK5GTk4Pp06dDR0eH/SEeDzNmzFBeUdSzV5m8ijw8zHmI+Nx4LtyTCpJg\nZ2wHD0sPeFp4sreWnhhkNghaPC1Vl0wI+R+5hnEWLlzIvYikPXv2KKa6loqisFe6RqYRqUWpeJDz\nAA9zHnKtqr4KXlZe8LT0hKeFJzwtPeEmdKO10wnpAWg2joarqqvCo7xHbKDnPuR67gIDAbysvOBl\n6cXeWnlhgMkAOgGJkB5KrrBPSkrCBx98gJycHCQmJiI+Ph6RkZFYu3atUooFKOzlUVhZiPvZ96WC\nPU2UBmeBs1SwD7MaBlM92RVHCSE9l1xhP3bsWGzbtg3vvfceHjx4AIZhMHToUCQmJiqlWIDCvqPy\nK/IRmx2L2KxY9jY7FsXVxfC28oa3tTcX7K5CV+j00VF1uYQQJZNrNk5lZSX8/f2lXkxbW1tx1ZEO\nyS3Pxf3s+1yox2bForSmFMOth8PHxgez3Wdj6+tb4ch3pC9NCSEy2g17oVCIlJQU7v7x48dhbW2t\n1KI0XU55jlRvPTYrFpV1lVywzx06F19P+BqDzAbR+DohpEPaHcZJTU3FsmXLcPPmTZiZmcHBwQEH\nDx6Evb298orSoGGc8tpyxGbF4o9Xf7Dt5R+oqq+Cj7UPfGx8MNyKDXgHUwcKdkJImxQyG6eiogKN\njY0wMlL+WZC9NewbGhvwpOAJ/nj5BxfuKUUp8LT0hL+tP9vs/CnYCSFd0qWw3759u9QLNPfxxx8r\nqLwWiuolYZ9dli3VY7+XdQ9Whlbwt/Pnwn2Y1TD68pQQohBd+oK2rKwMPB4PSUlJuHv3LqZNmwaG\nYXD27Fn4+fkprdieqrq+GrFZsbj18hYX7hV1FVyorxq9Cn62fuDr81VdKiFEA7U7jBMYGIioqChu\n+KasrAxTp07FtWvXlFeUmvfsGYZBRkkGbr+8jVsvb+HWy1t4lPcIrgJXBPQP4ALeie9EwzGEkG4j\n19TLvLw8qamW2trayMvLU1x1PYBkr/3Wy1u4lXkLDUwDAuwCEGAXgG0TtsHXxpeW6iWEqK12wz48\nPBx+fn6YMWMGGIbBqVOnlHqVKlUTr8d+K/OWVK/dReCCALsAzHCZgW0TttGXqISQHqVDs3FiY2Nx\n7do18Hgmst24AAAcTElEQVQ8jB07Ft7e3sotqhuHcRqZRjzJf4KYFzGIyYjBtRfXUNtQi4D+AVzP\n3dfGlxYCI4SoPVoITUJ9Yz3icuKkwt1EzwRjB47F2AFjETgwEI5mjtRrJ4T0OBod9jX1NbibdRcx\nL2JwLeMabmbehJ2xnVS42xnbKeS9CCFElTQq7CtqK3Dr5S225/4iBvey7sFF4IKxA8cicEAgxgwY\nA2E/oYIrJoQQ1VNI2FdXV4PH40FXV1ehxbWkM2FfXF2M6xnXuXBPyEuAt5U323MfOBaj+o+Csa6x\nkismhBDV61LYNzY24tSpU/jll19w8+ZNNDY2gmEY9OnTBwEBAZg/fz6mT5+ulLHttgrOq8jDtRfX\nEJMRg6vpV5EqSoW/rT8X7v62/tDX1ld4TYQQou66FPZjx45FYGAgpk2bBi8vL65HX1NTgwcPHiAy\nMhLXr19HTEyMUgvOLMnkvkyNeRGD7LJsjBkwhgv34dbDabkBQghBF8O+pqam3SGbjjynK3g8HsJP\nhiPmRQwqaiu4YB87cCw8LDzQR6uPwt+TEEJ6OrnG7FNSUmBnZwc9PT38/vvvSEhIQHh4OExNlXdJ\nOx6Ph5/u/YSxA8fC2dyZpkESQkgHyBX2w4YNQ2xsLNLT0zF16lSEhYUhMTERUVFRSikWUP+1cQgh\nRB21lZ3tXr9OS0sLffv2xYkTJ7BixQps27YN2dnZCi+SEEKI8rQb9jo6Ojh06BAiIiIQEhICAKir\nq1N6YYQQQhSn3bDfvXs3bt26hc8++wwODg5IS0vDO++80x21EUIIUZBWx+yXLVuGKVOm4PXXX1f4\npQg3b96MAwcOQEtLCx4eHtizZ4/UrB4asyeEkM7r0pj94sWL8fDhQ0ydOhXjx4/H1q1bERcXJ3cx\n6enp+Pnnn3H//n0kJCSgoaEBhw8flvt1CSGEtK7V9exHjhyJkSNH4osvvkBBQQEuXryI7du3Iz4+\nHt7e3pgyZQpmz57d6Tc0NjaGtrY2Kisr0adPH1RWVsLW1lauD0EIIaRt7V68BAAEAgHmzZuHefPm\nAQDu3buH//73v116Qz6fj7/97W8YMGAA9PX1MWnSJLz++utdei1CCCEd027Yi0QiREREID09HfX1\n9QDYcaHvvvuuS2+YmpqKHTt2ID09HSYmJpg1axYOHjyI+fPnSz1v3bp13HZQUBCCgoK69H6EENJb\nRUdHIzo6ukPPbfekqoCAAAQEBMDDwwNaWlpgGAY8Hq/LlyY8cuQIfvvtN+zatQsAsH//fty+fRs7\nd+5sKoq+oCWEkE6T64LjNTU1+OabbxRWjIuLCzZs2ICqqiro6enh0qVL8PPzU9jrE0IIkdVuz/7r\nr7+GsbExQkNDpaZH8vn8Lr/pV199hX379kFLSwvDhw/Hrl27oK2t3VQU9ewJIaTT5Fob54cffsBn\nn30GU1NTaGlpcS/4/PlzxVcqLorCnhBCOk2usHdwcMDdu3chEAiUUlxLKOwJIaTz5FoIbfDgwdDX\npys/EUJIT9buF7QGBgbw8vLCuHHjuDF7eaZeEkII6X7thv306dOlrjUrnnpJCCGk52g37IcOHQpf\nX1+pfWfOnFFaQYQQQhSv3TH7ZcuWISEhgbv/yy+/YMOGDUotihBCiGK1Oxvn+fPnmDlzJg4dOoRr\n164hIiICZ8+ehYmJifKKotk4hBDSaXJNvQSApKQkTJ8+HQMHDsSJEydgYGCg8CKliqKwJ4SQTutS\n2Ht4eEjdz8vLg6mpKXR0dMDj8RAfH6/4SsVFUdgTQkindSns09PT23xRe3t7eetqFYU9IYR0XpcW\nQjM3N2/3coRlZWUKv2QhIYT0ZgzDoKamBqWlpaiurkZtbS1qampQW1sr0zq7vy2thv2bb74JZ2dn\nhIWFwdfXl1v4rLCwEPfu3cOpU6fw7NkzXLp0SbFHghBC1BDDMKiqqkJJSQlKS0vbbW09T0tLC0ZG\nRtDX14eOjg7XdHV1pe63tk9yv7GxMbevLW1+QXvlyhUcOnQIN27cQFZWFgDAxsYGY8aMwfz585V2\nQREaxiGEKArDMKisrGw3iJvva+k5Ojo6MDExgbGxsVxNcgVhRZJ7Nk53o7AnhIiHOyRDV7zd0r7W\nHi8tLYWurq5M4DYP7ZZCXHKfkZGR1FLs6ojCnhDSrRoaGlBWVoaSkhKpAG4exu09rqWlJRW6kuHb\n0nbzfUZGRjA2Nlb7kFYUCntCSIfV1dVxwdtekwxpyX2VlZUwNDSEiYmJVBB3ZFvyVlnDHb0VhT0h\nGkLcoy4uLkZJSQl3K7nd0mOSrba2ViqEO9KaP9/Q0JC72BHpPnKF/ccff4wlS5bA3d1dKcW1hMKe\naKra2loUFxdzTRzKzVtr4V1RUQFDQ0OYmprCxMSEu21ru3nr168frWzbQ8l1wXFXV1csW7YMdXV1\nWLx4MebOnavUdXEI6cnq6uqkQlkkEkndthXcxcXFqK+vh6mpqVRYN282NjatBriRkRH1qEmLOjyM\n8/TpU+zduxeHDh3CmDFj8O6772LcuHHKKYp69kRFxNP0RCJRi62tEBeJRKiuruZC2czMTObWxMSE\nu99SoOvr61OvmnSZ3GP2DQ0NOHPmDPbs2YOXL19i9uzZuH79OgwMDHDkyJFuLZiQjqiurkZRURFE\nIpHMbWtN/HifPn1gZmbWamspyMXbhoaGFNZEZeQK+48++ghnzpzB+PHjsXTpUvj5+XGPOTs7Iykp\nSbHVgsKesBiGQWlpKRfEkq29fQ0NDeDz+eDz+S0Gdmv7zczMoKenp+qPTkiXyBX2e/bswezZs9Gv\nXz+Zx4qLi2FqaqqYKiWLorDvdSorK1FYWIiioiIUFhZ2aFskEkFfX18qtMXbkq2l/TQcQjSRXGE/\nfvx4XLlyRWpfcHAwLl++rLgKmxdFYa+2GIZBWVkZCgoKUFhYiIKCgla3JcObYRiYm5vD3NwcfD6f\n225+v/m2ppwMQ4gidGk2TlVVFSorK1FQUICioiJuf2lpKV69eqX4KolKVFVVoaCgAPn5+cjPz+e2\nxcHdUpDr6urC3NwcAoGAuxVve3p6cmEt3sfn82FgYEA9bUJUqNWw/+mnn/Dtt98iKysLPj4+3H4j\nIyP85S9/6ZbiSOeIe915eXlceEsGeEv36+rqIBQKIRQKIRAIpLaHDRsmE+bm5uY0pk1ID9TuMM73\n33+PFStWdFc9AGgYR1J1dTUX3nl5ee02XV3dFsO7tftGRkbU4yakl+jSmP2VK1cwfvx4/Prrry2G\nwYwZMxRbpWRRvTjsGYZBeXk5cnJykJubK9XE+yTDu6amBhYWFhAKhbCwsGizCYVC6Ovrq/ojEkJU\npEtj9levXsX48eNx5syZbg/7nqiiogLZ2dnIycmRCfLm93k8HqysrGBpack1KysreHl5cffFAW5s\nbEw9b0KI3GghtDYwDIPCwkIuxLOzs1vdrqurg7W1NaysrFoMcsn7hoaGqv5ohJBeSK6pl2vWrMGq\nVau4+fQikQjbt2/Hxo0bFV+puCglhz3DMCgqKkJWVhZevXqFrKwsrr169YoL8tzcXBgaGnIhbm1t\n3eq2iYkJ9cAJISolV9h7eXnh4cOHUvu8vb3x4MEDxVXYvCg5wr6qqgqZmZkyQd481PX19WFjYwNb\nW1vY2Nhw29bW1rCxseF66DTzhBDSU8i16mVjYyOqq6u50Kuqqmr3Kubdob6+HqmpqUhISJBqL1++\nhK2trVSI9+/fH/7+/tx9GxsbGBgYqPojEEJIt2k37OfPn4/g4GAsXrwYDMNgz549CA8P747aALBD\nLllZWXj06JFUqD99+hTW1tYYOnQoPDw8MGfOHGzcuBGDBw+msy4JIaSZDn1Be/78eW55hAkTJmDS\npEnKLYrHw/Lly7lg79u3Lzw8PLhg9/DwgLu7O33RSQghEuQaxgHYMfr6+npuuzs4Ojpi+vTp8PDw\ngKWlZbe8JyGE9Fbt9uyPHj2KTz75BK+99hoAICYmBtu2bcOsWbO6/KbFxcVYunQpEhMTwePxsHv3\nbowcObKpKDWZekkIIT2JXLNxPD09cenSJVhYWAAA8vPzERwcjPj4+C4XtGDBArz22mtYvHgx6uvr\nUVFRIXWpQwp7QgjpPLmGcRiGgVAo5O6bm5vLFcQlJSW4du0a9u3bxxbQty9d05YQQpSs3bCfPHky\nJk2ahHnz5oFhGBw5cgRTpkzp8humpaVBKBRi0aJFiIuLg4+PD7799luZqZDr1q3jtoOCghAUFNTl\n9ySEkN4oOjoa0dHRHXpuu8M4DMPgxIkTuH79Ong8HgIDA/Hmm292ubh79+4hICAAN2/exIgRI/DX\nv/4VxsbGWL9+fVNRNIxDCCGdJvcFxxUpJycHAQEBSEtLAwBcv34dW7ZswdmzZ5uKorAnhJBO69KY\nvaGhYatrvfB4PJSWlnapGCsrK/Tv3x/JyckYMmQILl26BHd39y69FiGEkI5RyaqXcXFxWLp0KWpr\na+Ho6Ig9e/bQbBxCCJGT3MM4165dQ0pKChYtWoT8/HyUl5fDwcFB4YVyRVHYE0JIp8kV9uvWrcO9\ne/eQnJyM5ORkvHr1CrNmzcLNmzeVUixAYU8IIV3RVnZqtffDJ0+eRGRkJPr16wcAsLW1RXl5uWIr\nJIQQolTthr2uri60tJqeVlFRodSCCCGEKF67YT9r1iz8+c9/RnFxMf7v//4PwcHBWLp0aXfURggh\nREE69AXtxYsXcfHiRQDApEmTMGHCBOUWRWP2hBDCqqkBiouBkhLp1sI+3t69Xf+Cdvv27Xjrrbdg\na2urlM/REgp7QkivUF8PlJY2BXNrt22FeGMjYGICmJqyt+LW/L6JCXhLlnR9IbSysjJMnDgRZmZm\neOuttzBr1ixaX54Q0vsxDFBdzYavuLUW2K09VlkJGBtLh7NkSJuaAjY2gJubTHBzz9PTA1o5wVXG\nkiWtPtThk6ri4uJw9OhRHD9+HHZ2dtyVq5SBevaEELkxDFBR0XJYN28t7S8pYUPW1FQ6qJuHdlu3\nhoaAVrtfjSqM3FeqAgALCwtYWVnB3Nwc+fn5CiuOEEJaxDBAVRUgErUc0JKtpeeUlAA6OoCZmWxY\ni+8LhcDgwa0/rqen6qOgMO327H/88UccPXoUeXl5mDVrFubMmQM3NzflFkU9e0J6h7o62UBu7bal\n0O7blw1ecWC31Zo/x8QE0NZW9RHoVnL17DMzM7Fjxw54eXkpvDBCiJpjGHbcWRzEzZvk/pYCvKam\n5UCWvLW3lw1rMzM2rHV1VX0Eeg2VLITWHurZE6JADAOUl7cc1u0FuEjE9q4lw1myNd/XPMj79ev4\nl4tEbmq1nn1HUNgT0oy4hy0SAUVFTUEsud1SKypiw1tXVzaoW2vNA5x61z0GhT0h6kI8hi0O6aIi\n6e22bvv0Afh86SBufr954/PZ8NbRUfUnJ92Awp4QRauqagrq5oHdVquslA3qtm4lw7wXzQwhykFh\nT0hLxFP7ioqAwsKmQG5tW7IxTFMYd6YZGdEYNlEaCnvS+9XWssEsbuKglrzfUpADgLk5G8TNb9va\n1tdX7eclpAUU9qTnYBj2ZJjCQqCgoOXgbinUq6ubAlncmgd4S8FtYKDqT0yIwlDYE9Wor5cOanF4\nt3VfJGIDWDK0m7fmoW5uTsMjhIDCnihCQwMbxAUFLTdxWEu2sjJ2Jog4kAWC9rf5fJo5QkgXUdgT\naQzDLrsqDuX8/Ja3JVtJCXtGo0Ag3cRh3fy+uTk7g6RPH1V/WkI0BoV9b1df3xTUkq2tMNfVZReB\nEoe0eFtyn2Sj4CZE7VHY9zQ1NbLB3VYrLWWHP4RC6dZWkNOcbUJ6HQp7VRP3vPPy2Jaf3/Z2VVVT\nOHekUa+bEAIKe8UTTw/MywNyc6VvJZs4vMU9bwsLtgmFrW8LheyXmjSzhBDSSRT2HVFfz4Zzbq50\ncLe2rafXFNKWlk23QmHTfXGQ8/nderUaQohm0tywr6trCmnJlpMju6+4mA1lS8um1jzMJQOdzqAk\nhKiZ3hX2DQ3snO6cnLZbbi47fCIQSAe4pSVgZSW7TyCgcW9CSI/WM8P+229bDvH8fHZM28qq5SYZ\n5ubmNHxCCNEYCrngeLd79owN7dGjpcPcwkLjritJCCHyUt+evfqVRQghaq2t7KQxDkII0QAU9oQQ\nogEo7AkhRANQ2BNCiAZQWdg3NDTA29sboaGhqipBLtHR0aouQeXoGNAxAOgYAD3jGKgs7L/99lu4\nubmB10PXgOkJv1xlo2NAxwCgYwD0jGOgkrB/+fIloqKisHTpUppiSQgh3UAlYf/RRx9h27Zt0KKz\nWwkhpHsw3ezMmTPMBx98wDAMw/z+++9MSEiIzHMAUKNGjRq1LrTWdPsZtGvWrMH+/fvRt29fVFdX\no7S0FH/6058QERHRnWUQQohGUelyCVevXsXXX3+NM2fOqKoEQgjRCCofNO+ps3EIIaQnUWnYv/ba\na4iMjFRlCR3W/LyAoqIiTJgwAUOGDMHEiRNRXFys4gqVr7i4GDNnzoSrqyvc3Nzwxx9/aNxx2Lx5\nM9zd3eHh4YF58+ahpqam1x+DxYsXw9LSEh4eHty+tj7z5s2bMXjwYLi4uODixYuqKFnhWjoGn3zy\nCVxdXTFs2DDMmDEDJSUl3GPqeAxU3rPvKZqfF7BlyxZMmDABycnJCA4OxpYtW1RcofKtXLkSU6dO\nxZMnTxAfHw8XFxeNOg7p6en4+eefcf/+fSQkJKChoQGHDx/u9cdg0aJFuHDhgtS+1j7z48ePceTI\nETx+/BgXLlzABx98gMbGRlWUrVAtHYOJEyciMTERcXFxGDJkCDZv3gxAjY+BUqbc9DKZmZlMcHAw\nc+XKFW72kLOzM5OTk8MwDMNkZ2czzs7OqixR6YqLixkHBweZ/Zp0HAoLC5khQ4YwRUVFTF1dHRMS\nEsJcvHhRI45BWloaM3ToUO5+a59506ZNzJYtW7jnTZo0ibl161b3FqskzY+BpBMnTjDz589nGEZ9\njwH17DugpfMCcnNzYWlpCQCwtLREbm6uqsrrFmlpaRAKhVi0aBGGDx+Od999FxUVFRp1HPh8Pv72\nt79hwIABsLGxgampKSZMmKBRx0Cstc+clZUFOzs77nl2dnZ49eqVSmrsTrt378bUqVMBqO8xoLBv\nx9mzZ2FhYQFvb+/Wr+3I4/X6L5rr6+tx//59fPDBB7h//z769esnM1zR249DamoqduzYgfT0dGRl\nZaG8vBwHDhyQek5vPwYtae8z9/bj8eWXX0JHRwfz5s1r9TnqcAwo7Ntx8+ZNREZGwsHBAXPnzsWV\nK1fwzjvvwNLSEjk5OQCA7OxsWFhYqLhS5bKzs4OdnR1GjBgBAJg5cybu378PKysrjTkO9+7dw6hR\no2Bubo6+fftixowZuHXrlkYdA7HW/v3b2toiMzOTe97Lly9ha2urkhq7w969exEVFYWDBw9y+9T1\nGFDYt2PTpk3IzMxEWloaDh8+jPHjx2P//v2YNm0a9u3bBwDYt28fpk+fruJKlcvKygr9+/dHcnIy\nAODSpUtwd3dHaGioxhwHFxcX3L59G1VVVWAYBpcuXYKbm5tGHQOx1v79T5s2DYcPH0ZtbS3S0tLw\n7Nkz+Pn5qbJUpblw4QK2bduG06dPQ09Pj9uvtsdA1V8a9CTR0dFMaGgowzDsl3XBwcHM4MGDmQkT\nJjAikUjF1Snfw4cPGV9fX8bT05N58803meLiYo07Dlu3bmXc3NyYoUOHMuHh4UxtbW2vPwZvvfUW\nY21tzWhrazN2dnbM7t272/zMX375JePo6Mg4OzszFy5cUGHlitP8GPznP/9hnJycmAEDBjBeXl6M\nl5cX8/7773PPV8djoJYXHCeEEKJYNIxDCCEagMKeEEI0AIU9IYRoAAp7QgjRABT2hCjRnDlzkJqa\nCgAoLy/H+++/DycnJ/j4+MDX1xe7du0CwK67o6+vD29vb7i7u+P9999v85Kd3333Hfbv398tn4H0\nDhT2pFfozoWmOvpeKSkpqKiogKOjIwBg6dKlMDc3R0pKCmJjY3HhwgUUFRVxz3dycsKDBw8QHx+P\nx48f49SpU62+9qJFi/D999/L90GIRqGwJ2otPT0dLi4uePvtt+Hm5oZZs2ahqqoKAGBvb4/Vq1fD\nx8cHx44dw8WLFzFq1Cj4+Phg9uzZqKioAACsXr0a7u7uGDZsGFatWgUAOHbsGDw8PODl5YWgoCAA\n7NmQK1as4N47JCQEMTExAABDQ0P8/e9/h5eXF27duoUDBw7A398f3t7eeO+991r8A3D48GFMmzYN\nALvUwt27d7Fx40bucYFAwNUjqU+fPhg1ahRSUlKwa9cu+Pn5wcvLCzNnzuQ+u5GREczNzZGYmCjv\nISYagsKeqL3k5GQsX74cjx8/hrGxMX788UcA7HojAoEAsbGxCA4OxpdffonLly8jNjYWPj4++Oab\nb1BUVIRTp05xS9GuXbsWALBhwwZcvHgRDx8+5K6p0Hz9Esn7lZWVGDlyJB4+fAg+n4+jR4/i5s2b\nePDgAbS0tKROlxe7ceMGfH19AQCJiYkYNmxYhz5vZWUlLl++DE9PT8yYMQN37tzBw4cP4erqiv/8\n5z/c8/z8/Lg/RoS0h8KeqL3+/fsjICAAAPD222/j+vXr3GNz5swBANy+fRuPHz/GqFGj4O3tjYiI\nCGRkZMDExAR6enpYsmQJTp48CX19fQDA6NGjsWDBAuzatQv19fXt1tCnTx/86U9/AgDuD4qvry+8\nvb1x5coVpKWlyfzMixcvYG1tDUD2D8mmTZvg7e0ttWZKamoqvL29MWbMGISEhGDSpElISEhAYGAg\nPD09cfDgQamevI2NDdLT0ztyCAlBX1UXQEh7JIOSYRip+/369eO2J0yYgEOHDsn8/J07d3D58mUc\nP34cP/zwAy5fvox///vfuHPnDs6dOwcfHx/Exsaib9++UsMx1dXV3Laenp7U+y5YsACbNm1qt3bx\nl6yurq6Ii4vj6l+zZg3WrFkDIyMj7rmOjo548OCB1M8vXLgQkZGR8PDwwL59+xAdHd3qsSCkLdSz\nJ2ovIyMDt2/fBgAcOnQIgYGBMs/x9/fHjRs3uJkvFRUVePbsGSoqKlBcXIwpU6bgm2++QVxcHAC2\nF+3n54cvvvgCQqEQL1++hL29PR4+fAiGYZCZmYk7d+60WE9wcDCOHz+O/Px8AOwl+jIyMmSeN3Dg\nQGRnZwNgv3z19fXF2rVruT8o4gXV2lJeXg4rKyvU1dXJLKecnZ0Ne3v7Nn+eEDHq2RO15+zsjJ07\nd2Lx4sXctERAuscvFAqxd+9ezJ07FzU1NQDYdcaNjIwQFhaG6upqMAyDf/3rXwCAVatW4dmzZ2AY\nBq+//jo8PT0BAA4ODnBzc4Orqyt8fHy415d8L1dXV2zcuBETJ05EY2MjtLW18eOPP2LAgAFSdY8Z\nMwb37t3jXmfXrl345JNP4OTkBHNzc+jr62Pbtm0tvofYhg0b4O/vD6FQCH9/f5SXl3OP3blzB19/\n/XXXDirROLQQGlFr6enpCA0NRUJCgqpL6bTnz59jxYoVOHfunMJfu7S0FMHBwbh7967CX5v0TjSM\nQ9ReTx2XHjRoEIyMjLihJUXau3cvVq5cqfDXJb0X9ewJIUQDUM+eEEI0AIU9IYRoAAp7QgjRABT2\nhBCiASjsCSFEA1DYE0KIBvj/gPy9bkmdZbgAAAAASUVORK5CYII=\n" } ], "prompt_number": 4 }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "3. Input composition" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We define\n", "composite object and name it \"rock\". A composite is made by\n", "giving burnman.composite a list of minerals and their molar fractions.\n", "Here \"rock\" has two constituent minerals: it is 80% Mg perovskite\n", "and 20% periclase. More minerals may be added by simply extending\n", "the list given to burnman.composite" ] }, { "cell_type": "code", "collapsed": false, "input": [ "rock = burnman.Composite([0.8, 0.2], \n", " [minerals.SLB_2011.mg_perovskite(), minerals.SLB_2011.periclase()])" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 5 }, { "cell_type": "markdown", "metadata": {}, "source": [ "At this point we want to tell the rock which equation of state to use for\n", "its thermoelastic calculations. In general, we recommend the 'slb3'\n", "equation of state as the most self-consistent model. The parameters from\n", "the SLB_2011 mineral library are fit using this model." ] }, { "cell_type": "code", "collapsed": false, "input": [ "rock.set_method('slb3')" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 6 }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "4. Input temperature" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now we get an array of temperatures at which will be used for computing\n", "the seismic properties of the rock. Here we use the Brown+Shankland (1981)\n", "geotherm for mapping pressure to temperature" ] }, { "cell_type": "code", "collapsed": false, "input": [ "temperature = burnman.geotherm.brown_shankland(pressure)\n", "plt.plot(pressure/1.e9,temperature,'r')\n", "plt.xlim(min(pressure)/1.e9,max(pressure)/1.e9)\n", "plt.xlabel('pressure (GPa)')\n", "plt.ylabel('temperature (K)');" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEMCAYAAADeYiHoAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VOXdxvHvCAEpi5Y1OqOEJoHsi0CCKDQYFqEYkU2j\nFFHUFlsEFZe2WqFKQHBpUGNVRBDUuAHBBUxFomULGMIiUQiYlGxQBZEEAgnkef84Oi+RbEAmk+X+\nXJdX4eTMmd+ccs2d5zybzRhjEBER+YUL3F2AiIjUTwoIERGpkAJCREQqpIAQEZEKKSBERKRCCggR\nEamQywIiJyeHAQMGEBgYSFBQEPPmzQNg+vTpOBwOwsPDCQ8PZ+XKlc7XzJo1C19fX/z8/EhOTnYe\nT0tLIzg4GF9fX6ZMmeKqkkVE5DQ2V82D2L9/P/v37ycsLIyioiJ69uzJ8uXLeeedd2jbti333Xdf\nufMzMjK4+eab2bx5M3l5eQwcOJDMzExsNhsRERE8//zzREREMGzYMO655x6uvfZaV5QtIiI/cVkL\nwtPTk7CwMADatGmDv78/eXl5AFSUSUlJScTGxuLh4YGXlxc+Pj6kpqZSUFBAYWEhERERAIwfP57l\ny5e7qmwREflJnfRBZGdnk56eTp8+fQB47rnnCA0NZeLEiRw+fBiA/Px8HA6H8zUOh4O8vLwzjtvt\ndmfQiIiI6zR39RsUFRUxevRo4uPjadOmDZMmTeLvf/87AI8++ij3338/r7766nm/j81mO+9riIg0\nRZX1NLi0BVFaWsqoUaMYN24cI0aMAKBz587YbDZsNht33HEHmzZtAqyWQU5OjvO1ubm5OBwO7HY7\nubm55Y7b7fYK388YU6//e+yxx9xeg7v/0z3QPdA9qF/3oCouCwhjDBMnTiQgIICpU6c6jxcUFDj/\nvGzZMoKDgwGIiYkhMTGRkpISsrKyyMzMJCIiAk9PT9q1a0dqairGGBYvXuwMGxERcR2XPWJat24d\nS5YsISQkhPDwcADi4uJ466232Lp1KzabjW7duvHSSy8BEBAQwNixYwkICKB58+YkJCQ4HxslJCQw\nYcIEiouLGTZsmEYwiYjUAZcNc61rNput2uaSu6WkpBAVFeXuMtxK90D3AHQPoP7cg6q+OxUQIiJN\nWFXfnVpqQ0REKqSAEBGRCikgRESkQgoIERGpkAJCREQqpIAQEZEKKSBERKRCCggREamQAkJERCqk\ngBARkQopIEREpEIKCBERqZACQkREKqSAEBGRCikgRESkQgoIERGpkAJCREQqpIAQEZEKNXd3ASIi\n4gbFxfDuu1WeohaEiEhTsmcPTJsGl18Ob75Z5akKCBGRxu7kSVi+HIYMgb594YILYONGWLWqypfp\nEZOISGNVUACvvGL9d9llcPfdkJQEF15Yo5crIEREGhNjYM0aePFF+PRTGDsWPvgAwsLO+lI2Y4xx\nQYl1zmaz0Ug+iojI2Tt8GBYtgn/9y3qEdPfdMG4cXHRRlS+r6rtTLQgRkYYsLc1qLbz/Plx7Lbz0\nEvTrBzbbeV9aASEi0tAcOwZvv20Fw4ED8Ic/wDffQJcutfo2esQkItJQZGZaofD66xARYT1GGjoU\nmjU750tW9d3psmGuOTk5DBgwgMDAQIKCgpg3b165nz/99NNccMEFHDp0yHls1qxZ+Pr64ufnR3Jy\nsvN4WloawcHB+Pr6MmXKFFeVLCJS/5SVwccfW0Fw1VXQogVs2mQdGz78vMKhWsZFCgoKTHp6ujHG\nmMLCQtO9e3eTkZFhjDFm3759ZsiQIcbLy8scPHjQGGPMzp07TWhoqCkpKTFZWVnG29vblJWVGWOM\n6d27t0lNTTXGGDN06FCzcuXKM97PhR9FRKTuHT5szLPPGuPjY0x4uDELFhhz7Fitv01V350ua0F4\nenoS9tOwqjZt2uDv709+fj4A9913H3PmzCl3flJSErGxsXh4eODl5YWPjw+pqakUFBRQWFhIREQE\nAOPHj2f58uWuKltExL0yMqxHR15e1mS2hQutjujbboNWreq0lDrppM7OziY9PZ3IyEiSkpJwOByE\nhISUOyc/P58+ffo4/+5wOMjLy8PDwwOHw+E8brfbycvLq4uyRUTqxqlT8OGH8Nxz8NVXcNddsHMn\nXHqpW8tyeUAUFRUxevRo4uPjueCCC4iLi+Pf//638+emFjuWp0+f7vxzVFQUUVFRtXZtEZFad+gQ\nvPoqJCRA584weTKMGQMtW7rsLVNSUkhJSanRuS4NiNLSUkaNGsW4ceMYMWIEO3bsIDs7m9DQUABy\nc3Pp2bMnqamp2O12cnJynK/Nzc3F4XBgt9vJzc0td9xut1f4fqcHhIhIvbVjh9VaePddq6P57bet\nUUl14Je/PM+YMaPSc13WB2GMYeLEiQQEBDB16lQAgoODOXDgAFlZWWRlZeFwONiyZQtdunQhJiaG\nxMRESkpKyMrKIjMzk4iICDw9PWnXrh2pqakYY1i8eDEjRoxwVdkiIq5x8qQ1mS0qyprQdtll1tyF\nxYvrLBzOlstaEOvWrWPJkiWEhIQQHh4OQFxcHEOHDnWeYzttpl9AQABjx44lICCA5s2bk5CQ4Px5\nQkICEyZMoLi4mGHDhnHttde6qmwRkdr1/ffWYnkvvmgtsT15MowcCR4e7q6sWpooJyLiCtu2QXw8\nLFsGN9wAf/4zXHGFu6s6g9ZiEhGpC6dOwYoVVjDs2WMNV83MhI4d3V3ZOVFAiIicr8OHrdFIzz8P\nnp4wZQqMGtUgHiNVRQEhInKudu+GefOsrTuHDq3T0Uh1QQEhInI2jIHkZOsx0pdfWpPaduyASobf\nN2QKCBGRmjh61FpFdd48a8G8KVOsYat1vPxFXVJAiIhU5b//tfoWXnvN2ojnxRfht7+tlQ156juX\nTZQTEWmwjIH//AdGj7aGpp46ZS2xvWyZNdGtCYQDqAUhIvL/SkogMRH++U8oKoJ77rFaDm3bursy\nt9BEORGRwkJ4+WUrGHr0gPvus5bDuKDxP2TRRDkRkYrs3291Or/8MgwcCMuXQ8+e7q6q3mj88Sgi\n8ku7d1vDU/394cgRq38hMVHh8AtqQYhI07FpEzz5JHzxhbUMxu7d0KmTu6uqtxQQItK4GQOrVlnB\nkJ1t9S+8/jq0bu3uyuo9BYSINE6lpdZjo7lzrc7mBx+0dmtr4Osj1SUFhIg0LkVFMH8+PPMM+Ppa\nATF4cJOZu1CbFBAi0jj873/WNp7/+hcMGABLl0KvXu6uqkHTKCYRadi+/dbqcPbzg4MHYcMGeOcd\nhUMtUECISMO0Ywfccou1vHb79tb+zgkJ4OPj7soaDQWEiDQsGzZATAwMGQKhoVYL4oknoHNnd1fW\n6KgPQkTqv5/3YJg1C/btgwcesB4jXXihuytr1BQQIlJ/nTpldTbPnm0tpPfww3DjjdBcX111QXdZ\nROqfkhJYssSa3Na+PTz2GAwf3iQWz6tPFBAiUn8cPQqvvAJPP22tk/TSS01mc576SAEhIu73ww/W\nrm3PPQf9+1sb82iYqtupvSYi7lNQYHU4+/hYo5G++ALee0/hUE8oIESk7u3bB5MmQWCg1d+Qnm7t\n3Obn5+7K5DQKCBGpOz8HQ3g4XHwx7NoF8fFw+eXurkwqoIAQEderKBhmzdJeDPWcAkJEXKeyYOjY\n0d2VSQ24LCBycnIYMGAAgYGBBAUFMW/ePAAeffRRQkNDCQsLIzo6mpycHOdrZs2aha+vL35+fiQn\nJzuPp6WlERwcjK+vL1OmTHFVySJSW3JyrAX0wsPhoousdZIUDA2PqaHi4mJz/Pjxmp5uCgoKTHp6\nujHGmMLCQtO9e3eTkZFhjhw54jxn3rx5ZuLEicYYY3bu3GlCQ0NNSUmJycrKMt7e3qasrMwYY0zv\n3r1NamqqMcaYoUOHmpUrV57xfmfxUUTEVfbtM2bSJGPatzfmoYeM+d//3F2RVKOq785KWxBlZWUs\nXbqUMWPGYLfb6datG127dsVutzN69GiWLVuGde2KeXp6EhYWBkCbNm3w9/cnPz+ftm3bOs8pKiqi\n40+/USQlJREbG4uHhwdeXl74+PiQmppKQUEBhYWFREREADB+/HiWL19eC9EoIrXm5xZDWBi0a2e1\nGGbPVh9DA1fpRLmoqCj69evHtGnTCAsLo2XLlgCcOHGC9PR0VqxYwbPPPssXX3xR7ZtkZ2eTnp5O\nZGQkAH/7299YvHgxrVq1YtOmTQDk5+fTp08f52scDgd5eXl4eHjgcDicx+12O3l5eRW+z/Tp08vV\nHxUVVW1tInIecnKsIEhMhDvusIJBoVCvpaSkkJKSUqNzKw2ITz/9lBYtWpxxvGXLloSHh9OnTx9O\nnDhR7RsUFRUxevRo4uPjadOmDQAzZ85k5syZzJ49m6lTp/Laa6/VqNjqnB4QIuJCPwfDW2/BnXcq\nGBqQX/7yPGPGjErPrfQR0+zZsys8/uOPPzJ48GAAZ6uiMqWlpYwaNYpx48YxYsSIM35+8803s3nz\nZsBqGZzeYZ2bm4vD4cBut5Obm1vuuN1ur/J9RcRF8vLgT3+yHiW1aWONSnrySYVDI1VpQPznP//h\nr3/9a7lj+/fv57e//S3XXHNNtRc2xjBx4kQCAgKYOnWq83hmZqbzz0lJSYSHhwMQExNDYmIiJSUl\nZGVlkZmZSUREBJ6enrRr147U1FSMMSxevLjCsBERFzpwAO69F0JC4Fe/sloMCobGr7Le6+LiYjNs\n2DBz7733GmOM2b17t/H29jYvvvhijXrG//Of/xibzWZCQ0NNWFiYCQsLMx9//LEZNWqUCQoKMqGh\noWbkyJHmwIEDztfMnDnTeHt7mx49ephVq1Y5j3/55ZcmKCjIeHt7m8mTJ591T7yInKPvvzfmwQet\nUUn33GNMQYG7K5JaVtV3p+2nEypUUlLCTTfdRMuWLVm/fj3PPvssI0eOrLv0Ogs2m63KUVUichYO\nH4ZnnoEXXoAxY+Bvf4PLLnN3VeICVX13VhoQTz/9NDabjdLSUubMmcPVV19N//79nRe87777XFfx\nOVBAiNSCwkKYNw/++U9rg56//x26dXN3VeJCVX13VjqKqbCwENtPm3RMnjwZm81GUVGRayoUEfc6\ndgwSEmDuXIiOhrVroUcPd1clblblI6aGRC0IkXNw4gS8/LK1DEbfvjB9OgQFubsqqUNVfXdWOopp\n+vTpHDhwoNKLFhQU8Nhjj51/dSJS90pLrWDw9YXkZPjoI2ujHoWDnKbSR0y9evXipptuoqSkhCuu\nuIJLLrkEYwz79+9ny5YttGzZkmnTptVlrSJyvk6ehDfegBkzrF3c3nkHTlvBQOR01T5iysnJYd26\ndezbtw+Arl27ctVVV5Vb/qI+0CMmkSqUlVlh8Nhj0KULPPGEtfezNHnnNIqpoVFAiFTAGFi1Cv7y\nF2jZ0gqGgQPhpwEoIuc0iklEGrhNm+Chh6CgwOqEHjFCwSBnRTvKiTQ2u3bB6NEwciTccgt89RXc\ncIPCQc6aAkKksSgogD/+Ea6+Gnr1gt27rSW4m+tBgZybagNi165dREdHExgYCMD27dt54oknXF6Y\niNTQjz9aS2EEBUHbtlYL4uGHrUX1RM5DtQFx5513EhcX59wbIjg4mLfeesvlhYlINY4ft9ZL8vW1\nWg/p6dZM6Pbt3V2ZNBLVtj2PHTvm3AkOrB5vDw8PlxYlIlU4dQqWLLHWSQoLgzVr4KcWvkhtqjYg\nOnXqxJ49e5x/f++997jkkktcWpSIVMAY+Phj6/FRu3bWhLerr3Z3VdKIVTsPYu/evdx1111s2LCB\niy++mG7duvHGG2/g5eVVRyXWjOZBSKO2YYM1ZPXgQWvI6nXXaVSS1Ipzngdx6tQpXnzxRVavXk1R\nURFlZWW0a9fOJUWKSAV27bImuW3ebC2PMX68RiVJnamyk7pZs2asXbsWYwxt2rRROIjUlYMH4Z57\n4KqrrLWSdu+G229XOEidqvZfW1hYGNdffz1jxozhVz8Nm7PZbPV2ZzmRBq2kBJ5/3nqMNHYsfP21\n9n0Wt6k2II4fP0779u357LPPyh1XQIjUImNg+XJ44AFro57PP4eAAHdXJU2cFusTcbe0NLjvPjh0\nCJ5+GgYPdndF0oSc12J9t9122xkXA1iwYEEtlCbShOXmWjOgk5PhH/+w+hiaNXN3VSJO1QbE7373\nO2coFBcXs2zZMi699FKXFybSaBUVWTOen3/eWjtp925riQyReuasHzGVlZVx1VVXsWHDBlfVdE70\niEnqvVOn4PXX4ZFHICoK4uKga1d3VyVNXK3uB7F7926+++678y5KpElZs8bqZ/jVr2DpUjht+RqR\n+qragGjTpo3zEZPNZqNLly48+eSTLi9MpFHYvdsambRjBzz5pLVPg2ZASwNRbUAUFRXVRR0ijcsP\nP1gzn5cssZbIePttuPBCd1clclaqXe47Ojq6RsdEBKuf4ZVXwN/fWo7766+tFoTCQRqgSlsQxcXF\nHDt2jO+++45Dhw45jx85coS8vLw6KU6kQVm/3loeo1UrWLkSwsPdXZHIeam0BfHSSy/Rq1cvdu3a\nRc+ePZ3/xcTE8Oc//7naC+fk5DBgwAACAwMJCgpi3rx5ADzwwAP4+/sTGhrKyJEj+fHHH52vmTVr\nFr6+vvj5+ZGcnOw8npaWRnBwML6+vkyZMuV8Pq9I7SsosBbRGzvW6oj+4guFgzQOphrx8fHVnVKh\ngoICk56ebowxprCw0HTv3t1kZGSY5ORkc+rUKWOMMQ899JB56KGHjDHG7Ny504SGhpqSkhKTlZVl\nvL29TVlZmTHGmN69e5vU1FRjjDFDhw41K1euPOP9avBRRGrXiRPGzJljTIcOxvzlL8YUFrq7IpGz\nVtV3Z7Wd1Pfccw9fffUVGRkZHD9+3Hl8/PjxVb7O09MTT09PwBoJ5e/vT35+PoMGDXKeExkZyfvv\nvw9AUlISsbGxeHh44OXlhY+PD6mpqXTt2pXCwkIiIiKc77t8+XKuvfbas81CkdqzciVMnQrdu8PG\njeDj4+6KRGpdtQExffp0Pv/8c3bu3Mnvfvc7Vq5cydVXX11tQJwuOzub9PT0cluXgrVcR2xsLAD5\n+fn06dPH+TOHw0FeXh4eHh44HA7ncbvdXmkfyPTp051/joqKIioqqsY1itTInj1w773WPg3//CcM\nG+buikTOSkpKCikpKTU6t9qAeO+999i2bRtXXHEFr732GgcOHOCWW26pcTFFRUWMHj2a+Ph42rRp\n4zw+c+ZMWrRowc0331zja1Xn9IAQqVVFRdbM55dfhgcfhPffhxYt3F2VyFn75S/PM2bMqPTcagOi\nVatWNGvWjObNm/Pjjz/SuXNncnJyalRIaWkpo0aNYty4cYwYMcJ5fOHChXz88cesXr3aecxut5e7\nbm5uLg6HA7vdTm5ubrnjdru9Ru8vct6MgcREKxQGDIDt20FrkUkTUe08iN69e/PDDz9w55130qtX\nL8LDw+nbt2+1FzbGMHHiRAICApg6darz+KpVq5g7dy5JSUlceNrY8JiYGBITEykpKSErK4vMzEwi\nIiLw9PSkXbt2pKamYoxh8eLF5cJGxGW2boX+/eGpp6yQeP11hYM0KVUu1meMIScnh8svvxyArKws\njhw5QmhoaLUXXrt2Lf379yckJMS5VEdcXBz33HMPJSUltG/fHoArr7yShIQE588XLFhA8+bNiY+P\nZ8iQIYA1zHXChAkUFxczbNgw55DZch9Ei/VJbTl4EB591Foz6fHHtQy3NGpVfXdWGxDBwcF89dVX\nLiuutigg5LyVlcFrr1l7NIwZY+3R8Otfu7sqEZc659VcbTYbPXv2ZNOmTc5hpiKN0tatcPfdVp+D\nZkGLADXYD6JHjx7s2bOHrl270rp1a+tFNhvbt2+vkwJrSi0IOSc//gh//7vVxxAXB7fdBhdU2zUn\n0mic134Qn3zySa0XJOJ2xsBbb8G0aTB8OGRkQIcO7q5KpF6p9lclLy8vcnJyWLNmDV5eXrRu3Vq/\nqUvD9vXXEB1tjU5autSa26BwEDlDtQExffp05syZw6xZswAoKSlh3LhxLi9MpNYdPQoPP2wNXR05\nEjZvhtNm74tIedUGxLJly0hKSnL2P9jtdgoLC11emEitMcZqKQQEQF6etbvbn/+soasi1ai2D6Jl\ny5ZccFqn3dGjR11akEit2rsXJk+G//4XFi0Crc8lUmPVtiDGjBnDH/7wBw4fPszLL79MdHQ0d9xx\nR13UJnLujh+3tvyMjLRCIT1d4SBylqod5gqQnJzs3MBnyJAh5Zbsri80zFWcVq60Wg1hYfDss3DZ\nZe6uSKTeOq9hrgDBwcEUFxdjs9kIDg6u1eJEak1BgbXlZ3o6PP88aM8QkfNS7SOm+fPnExkZydKl\nS3n//feJjIzk1VdfrYvaRGqmrMwaqhoaCj16wFdfKRxEakG1j5i6d+/Ohg0b6PDTOPGDBw9y5ZVX\nsnv37jopsKb0iKmJ+uYbuOsuKCmBV14BtXBFzkpV353VtiA6duxYbqOfNm3a0LFjx9qrTuRclJRY\ni+n16wdjx8K6dQoHkVpWbR+Et7c3ffr04frrrwesvaNDQkJ4+umnsdls3HfffS4vUqSc9evhzjvB\n2xu2bFEntIiL1CggvL29nXs6XH/99dhsNoqKilxenEg5R47AX/4Cy5ZBfDyMHg0//bsUkdpXo2Gu\nDYH6IBq55cutoavXXgtz5mifBpFacl7DXDdv3kxcXBzZ2dmcPHnSecH6tty3NFL5+VYwfPUVLFkC\nv/2tuysSaTKqDYhbbrmFp556iqCgoHJLboi4VFmZNSrpkUfgj3+EN96A0/YwFxHXqzYgOnXqRExM\nTF3UImL5eehqaSmsWQNBQe6uSKRJqrYPIjk5mbfffpuBAwfSokUL60U2GyNHjqyTAmtKfRCNQEkJ\nzJ4N8+bB9OkwaZJWXBVxsfPqg1i0aBG7du3i5MmT5R4x1beAkAZu82Zru8/f/MZaKkNDV0XcrkZ7\nUn/zzTfOYa71lVoQDdTx4/DYY7BwoTV09cYbNXRVpA6d10zqvn37kpGRUetFibB+vbXi6rffwvbt\ncNNNCgeReqTaFoSfnx979+6lW7dutGzZ0npRPRzmqhZEA3LsmDU66a234LnnrAlvIuIW59UHsWrV\nqlovSJqwL76AiROhd29r60+t6yVSb1X7iMnLy4ucnBzWrFmDl5cXrVu31m/qcvaKiqwJb7Gx8NRT\n8OabCgeReq7agJg+fTpz5sxh1qxZAJSUlDBu3DiXFyaNyGefQUiItZbSjh3w08KPIlK/VfuIadmy\nZaSnp9OzZ08A7HY7hYWFLi9MGoHCQnjwQfjwQ/jXv+B3v3N3RSJyFqptQbRs2bLc/IejR4/W6MI5\nOTkMGDCAwMBAgoKCmDdvHgDvvvsugYGBNGvWjC1btpR7zaxZs/D19cXPz8+5BzZAWloawcHB+Pr6\nMmXKlBq9v7hZcrI1A7q01FpHSeEg0vCYasyZM8fcddddxsvLy7z00ksmMjLSxMfHV/cyU1BQYNLT\n040xxhQWFpru3bubjIwM8/XXX5tdu3aZqKgok5aW5jx/586dJjQ01JSUlJisrCzj7e1tysrKjDHG\n9O7d26SmphpjjBk6dKhZuXLlGe9Xg48ideGHH4y5/XZjLr/cmE8+cXc1IlKNqr47q33E9MADD5Cc\nnEzbtm3ZvXs3jz/+OIMGDao2eDw9PfH09ASsXej8/f3Jz88nOjq6wvOTkpKIjY3Fw8MDLy8vfHx8\nSE1NpWvXrhQWFhIREQHA+PHjWb58Oddqz+H656OPrIX1hg+3+hratXN3RSJyHqoNiIceeognn3yS\nwYMHn3GsprKzs0lPTycyMrLSc/Lz8+nTp4/z7w6Hg7y8PDw8PHA4HM7jdrudvLy8Cq8xffp055+j\noqKIioqqcY1yHg4fhqlTrSGsixbBNde4uyIRqURKSgopKSk1OrfagEhOTj4jDD7++OMaB0RRURGj\nR48mPj6+3N7WrnB6QEgd+eQTuOMOiImxZkO7+P9jETk/v/zlecaMGZWeW2lAvPjiiyQkJLB3716C\nT9sMvrCwkKuuuqpGhZSWljJq1CjGjRvHiBEjqjzXbreTk5Pj/Htubi4OhwO73U5ubm6543a7vUbv\nLy5UWAj3328FxGuvwcCB7q5IRGpZpaOYbr75Zj744ANiYmL48MMP+eCDD/jggw9IS0vjjTfeqPbC\nxhgmTpxIQEAAU6dOrfScn8XExJCYmEhJSQlZWVlkZmYSERGBp6cn7dq1IzU1FWMMixcvrjZsxMXW\nrLHmNZSVWX0NCgeRRslle1KvXbuW/v37ExIS4lwJNi4ujhMnTjB58mS+//57LrroIsLDw1m5cqXz\n5wsWLKB58+bEx8czZMgQwBrmOmHCBIqLixk2bJhzyGy5D6K1mFzv6FF4+GFYtgxefhmGDXN3RSJy\nnqr67nRZQNQ1BYSLrVsHEyZAnz7Whj6//rW7KxKRWnBei/VJE1dcDI8+aq2dlJAAerwn0mRUO5Na\nmrBNm+CKK2DfPmuEksJBpElRC0LOdOIE/OMfMH++9TjpxhvdXZGIuIECQspLT4dbb4Vu3WDbNvhp\nNryIND16xCSW0lKYMQOGDIFp02D5coWDSBOnFoTAN9/ALbdAp06wZQuctrSJiDRdakE0ZcZY/Qz9\n+sGdd8LKlQoHEXFSC6KpOnwY7rrLaj18/jkEBLi7IhGpZ9SCaIrWrYOwMOjSxRrKqnAQkQqoBdGU\nnDoFM2daE95eeQWuu87dFYlIPaaAaCpycmDcOGjeHNLSQCviikg19IipKVi6FHr1gqFDrb2iFQ4i\nUgNqQTRmx47BfffBv/8NK1ZAFTv6iYj8kloQjdX27dC7t7Wxz5YtCgcROWsKiMbGGHj+eYiOhoce\ngiVL4KKL3F2ViDRAesTUmHz/Pdx+O+Tnw/r14Ovr7opEpAFTC6Kx+Owza26Dn5/CQURqhVoQDd3J\nk/DYY/Daa7BwIQwe7O6KRKSRUEA0ZDk5EBsLrVtby3R36eLuikSkEdEjpobqgw+suQ3Dh1uL7Ckc\nRKSWqQVmDEevAAARY0lEQVTR0JSUwMMPw/vvWxPgrrrK3RWJSCOlgGhIvv3W2v7z0kutR0rt27u7\nIhFpxPSIqaF45x1rstu4cdZubwoHEXExtSDqu+Li/18uY+VKq99BRKQOqAVRn33zDfTpA4cOWSuw\nKhxEpA4pIOqr11+3tgL9058gMVHLZYhIndMjpvqmqAj+/GdITbVmRwcHu7siEWmi1IKoT35egdVm\ngy+/VDiIiFu5LCBycnIYMGAAgYGBBAUFMW/ePAAOHTrEoEGD6N69O4MHD+bw4cPO18yaNQtfX1/8\n/PxITk52Hk9LSyM4OBhfX1+mTJniqpLdxxh46SVrBda//MVaNqN1a3dXJSJNnXGRgoICk56ebowx\nprCw0HTv3t1kZGSYBx54wDz55JPGGGNmz55tHnroIWOMMTt37jShoaGmpKTEZGVlGW9vb1NWVmaM\nMaZ3794mNTXVGGPM0KFDzcqVK894Pxd+FNc6csSYG280JiTEmK+/dnc1ItLEVPXd6bIWhKenJ2Fh\nYQC0adMGf39/8vLyWLFiBbfeeisAt956K8uXLwcgKSmJ2NhYPDw88PLywsfHh9TUVAoKCigsLCQi\nIgKA8ePHO1/T4H39NUREQNu2sHGjtRKriEg9USd9ENnZ2aSnpxMZGcmBAwfo8tO6QV26dOHAgQMA\n5Ofn43A4nK9xOBzk5eWdcdxut5OXl1cXZbvWO+9A//7wwAPwyivQqpW7KxIRKcflo5iKiooYNWoU\n8fHxtG3bttzPbDYbNput1t5r+vTpzj9HRUURFRVVa9euNaWl8OCDkJQEn3wCV1zh7opEpAlJSUkh\nJSWlRue6NCBKS0sZNWoUv//97xkxYgRgtRr279+Pp6cnBQUFdO7cGbBaBjk5Oc7X5ubm4nA4sNvt\n5Obmljtut9srfL/TA6JeKiiAsWOhXTtrlJKWyxCROvbLX55nzJhR6bkue8RkjGHixIkEBAQwdepU\n5/GYmBgWLVoEwKJFi5zBERMTQ2JiIiUlJWRlZZGZmUlERASenp60a9eO1NRUjDEsXrzY+ZoG5Ysv\nrJnQgwdbS3UrHESknrP91Itd69auXUv//v0JCQlxPkaaNWsWERERjB07ln379uHl5cU777zDxRdf\nDEBcXBwLFiygefPmxMfHM2TIEMAa5jphwgSKi4sZNmyYc8hsuQ9is+Gij3J+jIGnn4annoJFi+Cn\nzyQiUh9U9d3psoCoa/UyII4cgdtvh//+F957D7p2dXdFIiLlVPXdqZnUrrJzpzUrumNHWLtW4SAi\nDY4CwhXeeguioqxZ0f/6F7Rs6e6KRETOmhbrq00lJTBtGnz0kbV/w08TBUVEGiIFRG3Jy4MxY6BD\nB2sI669/7e6KRETOix4x1YY1a6z+huHDrQlwCgcRaQTUgjgfxsBzz0FcHCxZAgMHursiEZFao4A4\nVydOwN13w+bNsGEDdOvm7opERGqVHjGdi4ICa5TS4cOwfr3CQUQaJQXE2dq82Vqie9gwePddaNPG\n3RWJiLiEHjGdjSVL4N57reW5G+J6UCIiZ0EBUROnTsHDD8PSpdaIpaAgd1ckIuJyCojq/PADxMbC\nyZOwaZM1z0FEpAlQH0RVvv4aIiPB3x9WrVI4iEiTooCozIcfwm9/a62n9Oyz0FyNLRFpWvSt90vG\nwOzZ8MILsGIF9Onj7opERNxCAXG6Y8es/Ru+/RZSU6GSrU1FRJoCPWL62b59cPXV0KKFtT2owkFE\nmjgFBFgb+vTpA+PGWduCXnihuysSEXE7PWJavBjuv9/6X+0XLSLi1HQDwhiYPt0KhpQUCAhwd0Ui\nIvVK0wyI48f/vzN6wwbo0sXdFYmI1DtNrw/iu++sfRtOnrSWzVA4iIhUqGkFxK5dcOWV0L8/JCZC\nq1burkhEpN5qOo+Y1qyBm26yJsHddpu7qxERqfeaRgti4UIrHN56S+EgIlJDjbsFUVYGjz5qPU76\n/HPw83N3RSIiDUbjDYjiYpgwAXJzYeNG6NTJ3RWJiDQojfMR0//+B9dcAxdcAKtXKxxERM6BSwPi\n9ttvp0uXLgQHBzuPbdu2jSuvvJKQkBBiYmIoLCx0/mzWrFn4+vri5+dHcnKy83haWhrBwcH4+voy\nZcqUqt80I8NaNmPQIHjjjXq1bEZKSoq7S3A73QPdA9A9gIZxD1waELfddhurVq0qd+yOO+5gzpw5\nbN++nRtuuIG5c+cCkJGRwdtvv01GRgarVq3i7rvvxhgDwKRJk3j11VfJzMwkMzPzjGs6ffopREVZ\nM6T/8Q+rBVGPNIR/EK6me6B7ALoH0DDugUu/Qfv168evf/3rcscyMzPp168fAAMHDuT9998HICkp\nidjYWDw8PPDy8sLHx4fU1FQKCgooLCwkIiICgPHjx7N8+fKK3/CWW+Ddd2H8eNd9KBGRJqLOf8UO\nDAwkKSkJgHfffZecnBwA8vPzcTgczvMcDgd5eXlnHLfb7eTl5VV88S++sHaBExGR81bno5gWLFjA\nPffcw+OPP05MTAwtWrSotWvbGsAw1hkzZri7BLfTPdA9AN0DqP/3oM4DokePHnzyyScA7N69m48+\n+giwWgY/tyYAcnNzcTgc2O12cnNzyx23V7CZz8/9FSIiUjvq/BHTd999B0BZWRlPPPEEkyZNAiAm\nJobExERKSkrIysoiMzOTiIgIPD09adeuHampqRhjWLx4MSNGjKjrskVEmhyXtiBiY2P5/PPP+f77\n77nsssuYMWMGRUVFvPDCCwCMGjWKCRMmABAQEMDYsWMJCAigefPmJCQkYLPZAEhISGDChAkUFxcz\nbNgwrr32WleWLSIiAEZc5uTJkyYsLMwMHz7cGGPMwYMHzcCBA42vr68ZNGiQ+eGHH9xcoev98MMP\nZtSoUcbPz8/4+/ubjRs3Nrn7EBcXZwICAkxQUJCJjY01x48fb/T34LbbbjOdO3c2QUFBzmNVfea4\nuDjj4+NjevToYT755BN3lFzrKroH06ZNM35+fiYkJMTccMMN5vDhw86f1cd7UL8mCjQy8fHxBAQE\nOFtCs2fPZtCgQezevZvo6Ghmz57t5gpdb8qUKQwbNoyvv/6a7du34+fn16TuQ3Z2Nq+88gpbtmxh\nx44dnDp1isTExEZ/DyqaA1XZZ65oDlRZWZk7yq5VFd2DwYMHs3PnTrZt20b37t2ZNWsWUI/vgbsT\nqrHKyckx0dHR5rPPPnO2IHr06GH2799vjDGmoKDA9OjRw50lutzhw4dNt27dzjjelO7DwYMHTffu\n3c2hQ4dMaWmpGT58uElOTm4S9yArK6vcb8+Vfea4uDgze/Zs53lDhgwxGzZsqNtiXeSX9+B0S5cu\nNbfccosxpv7eA7UgXOTee+9l7ty5XHDabO4DBw7Q5acd7Lp06cKBAwfcVV6dyMrKolOnTtx2221c\nccUV3HnnnRw9erRJ3Yf27dtz//33c/nll3PppZdy8cUXM2jQoCZ1D35W2WeubA5UY7dgwQKGDRsG\n1N97oIBwgQ8//JDOnTsTHh5e6fBbm83mfPTUWJ08eZItW7Zw9913s2XLFlq3bn3Go5TGfh/27t3L\nP//5T7Kzs8nPz6eoqIglS5aUO6ex34OKVPeZG/v9mDlzJi1atODmm2+u9Jz6cA8UEC6wfv16VqxY\nQbdu3YiNjeWzzz7j97//PV26dGH//v0AFBQU0LlzZzdX6loOhwOHw0Hv3r0BGD16NFu2bMHT07PJ\n3Icvv/ySvn370qFDB5o3b87IkSPZsGFDk7oHP6vs339Fc6AqmuvUWCxcuJCPP/6YN954w3msvt4D\nBYQLxMXFkZOTQ1ZWFomJiVxzzTUsXryYmJgYFi1aBMCiRYsa/XwOT09PLrvsMnbv3g3Ap59+SmBg\nINddd12TuQ9+fn5s3LiR4uJijDF8+umnBAQENKl78LPK/v1XNgeqMVq1ahVz584lKSmJC09babre\n3gN3d4I0dikpKea6664zxlgdltHR0Y12aGNFtm7danr16lVuWF9Tuw9PPvmkc5jr+PHjTUlJSaO/\nBzfddJO55JJLjIeHh3E4HGbBggVVfuaZM2cab29v06NHD7Nq1So3Vl57fnkPXn31VePj42Muv/xy\nExYWZsLCwsykSZOc59fHe2AzRmtUiIjImfSISUREKqSAEBGRCikgRESkQgoIERGpkAJCpJ658cYb\n2bt3LwBFRUVMmjQJHx8fevbsSa9evZg/fz5grfPUqlUrwsPDCQwMZNKkSVXuizJv3jwWL15cJ59B\nGgcFhDRZdbkYWk3fa8+ePRw9ehRvb28A7rjjDjp06MCePXtIS0tj1apVHDp0yHm+j48P6enpbN++\nnYyMjMr3a8daPO655547vw8iTYoCQhqd7Oxs/Pz8GDduHAEBAYwZM4bi4mIAvLy8ePjhh+nZsyfv\nvvsuycnJ9O3bl549ezJ27FiOHj0KwMMPP0xgYCChoaE8+OCDgLWHenBwMGFhYURFRQHWrNjJkyc7\n33v48OF88cUXALRp04Zp06YRFhbGhg0bWLJkCZGRkYSHh/PHP/6xwtBITEwkJiYGsJbp2Lx5M088\n8YTz5x07dnTWc7pmzZrRt29f9uzZw/z584mIiCAsLIzRo0c7P3vbtm3p0KEDO3fuPN9bLE2EAkIa\npd27d/OnP/2JjIwM2rVrR0JCAmCtb9OxY0fS0tKIjo5m5syZrF69mrS0NHr27MkzzzzDoUOHWL58\nuXNZ5kceeQSAxx9/nOTkZLZu3cqKFSuc1zvd6X8/duwYffr0YevWrbRv35533nmH9evXk56ezgUX\nXFBuqYWfrVu3jl69egGwc+dOQkNDa/R5jx07xurVqwkJCWHkyJFs2rSJrVu34u/vz6uvvuo8LyIi\nwhlgItVRQEijdNlll3HllVcCMG7cONauXev82Y033gjAxo0bycjIoG/fvoSHh/P666+zb98+Lrro\nIi688EImTpzIsmXLaNWqFQBXXXUVt956K/Pnz+fkyZPV1tCsWTNGjRoF4AyhXr16ER4ezmeffUZW\nVtYZr/nvf//LJZdcApwZPnFxcYSHh5dbo2fv3r2Eh4dz9dVXM3z4cIYMGcKOHTvo168fISEhvPHG\nG+VaDJdeeinZ2dk1uYUirt1yVMRdTv9yNcaU+3vr1q2dfx40aBBvvvnmGa/ftGkTq1ev5r333uP5\n559n9erVvPjii2zatImPPvqInj17kpaWRvPmzcs9Kjp+/LjzzxdeeGG597311luJi4urtvafO5r9\n/f3Ztm2bs/6//vWv/PWvf6Vt27bOc729vUlPTy/3+gkTJrBixQqCg4NZtGgRKSkpld4LkaqoBSGN\n0r59+9i4cSMAb775Jv369TvjnMjISNatW+ccMXT06FEyMzM5evQohw8fZujQoTzzzDNs27YNsH5b\nj4iIYMaMGXTq1Inc3Fy8vLzYunUrxhhycnLYtGlThfVER0fz3nvv8d133wFw6NAh9u3bd8Z5Xbt2\npaCgALA6oHv16sUjjzziDKGfF/2rSlFREZ6enpSWlp6xtHhBQQFeXl5Vvl7kZ2pBSKPUo0cPXnjh\nBW6//XbnEFAo37Lo1KkTCxcuJDY2lhMnTgDWOv1t27bl+uuv5/jx4xhjePbZZwF48MEHyczMxBjD\nwIEDCQkJAaBbt24EBATg7+9Pz549ndc//b38/f154oknGDx4MGVlZXh4eJCQkMDll19eru6rr76a\nL7/80nmd+fPn88ADD+Dj40OHDh1o1aoVc+fOrfA9fvb4448TGRlJp06diIyMpKioyPmzTZs28dRT\nT53bTZUmR4v1SaOTnZ3Nddddx44dO9xdyln79ttvmTx5Mh999FGtX/vIkSNER0ezefPmWr+2NE56\nxCSNUkN9zv6b3/yGtm3bOh971aaFCxcyZcqUWr+uNF5qQYiISIXUghARkQopIEREpEIKCBERqZAC\nQkREKqSAEBGRCikgRESkQv8HGrRzMiMqgSUAAAAASUVORK5CYII=\n" } ], "prompt_number": 7 }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "4. Calculate velocities" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " \n", "Here is the step which does the heavy lifting. burnman.velocities_from_rock sets the state of the rock at each of the pressures and temperatures defined,then calculates the elastic moduli and density of each individual phase. After that,it performs elastic averaging on the phases to get a single bulk and shear\n", "modulus for the rock. This averaging scheme defaults to Voigt-Reuss-Hill,\n", "but see example_averaging.py for other options. Finally, it calculates the seismic\n", "wave speeds for the whole rock. It returns a tuple of density, p-wave velocity\n", "s-wave velocity, bulk sound speed, bulk modulus, and shear modulus." ] }, { "cell_type": "code", "collapsed": false, "input": [ "density, vp, vs, vphi, K, G = burnman.velocities_from_rock(rock, pressure, temperature)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 8 }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "5. Plot results" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "All the work is done except the plotting! Here we want to plot the seismic wave\n", "speeds and the density against PREM using the matplotlib plotting tools. " ] }, { "cell_type": "code", "collapsed": false, "input": [ "# First, we plot the s-wave speed verses the PREM s-wave speed\n", "plt.figure(figsize=(15,5))\n", "plt.subplot(1,3,1)\n", "plt.plot(pressure/1.e9,vs/1.e3,color='b',linestyle='-',marker='o', markerfacecolor='b',markersize=4,label='computation')\n", "plt.plot(pressure/1.e9,seis_vs/1.e3,color='k',linestyle='-',marker='o', markerfacecolor='k',markersize=4,label='reference')\n", "plt.title(\"S wave speed (km/s)\")\n", "plt.xlim(min(pressure)/1.e9,max(pressure)/1.e9)\n", "plt.xlabel('pressure (GPa)')\n", "plt.legend(loc='lower right')\n", "plt.ylim(5,8.0)\n", " \n", "\n", "# Next, we plot the p-wave speed verses the PREM p-wave speed\n", "plt.subplot(1,3,2)\n", "plt.plot(pressure/1.e9,vp/1.e3,color='b',linestyle='-',marker='o',markerfacecolor='b',markersize=4)\n", "plt.plot(pressure/1.e9,seis_vp/1.e3,color='k',linestyle='-',marker='o',markerfacecolor='k',markersize=4)\n", "plt.title(\"P wave speed (km/s)\")\n", "plt.xlabel('pressure (GPa)')\n", "plt.xlim(min(pressure)/1.e9,max(pressure)/1.e9)\n", "plt.ylim(10,16)\n", " \n", "# Next, we plot the density verses the PREM density\n", "plt.subplot(1,3,3)\n", "plt.plot(pressure/1.e9,density/1.e3,color='b',linestyle='-',marker='o', markerfacecolor='b',markersize=4)\n", "plt.plot(pressure/1.e9,seis_rho/1.e3,color='k',linestyle='-',marker='o', markerfacecolor='k',markersize=4)\n", "plt.xlim(min(pressure)/1.e9,max(pressure)/1.e9)\n", "plt.xlabel(\"Pressure (GPa)\")\n", "\n", "plt.title(\"density (kg/m^3)\");" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAA2YAAAFPCAYAAAAm6d0pAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XlYVNX/wPH3oKggIyhuCChupSaoqaigiaWGu2IuaGpm\nJmZm2mJ9za1o+5UtWomSZi4ILrghorjgjmhaIG5p7vuGsqjIcH5/TEwMmysMDJ/X8/A8zNx7534u\nes7cc885n6NRSimEEEIIIYQQQpiMhakDEEIIIYQQQojiThpmQgghhBBCCGFi0jATQgghhBBCCBOT\nhpkQQgghhBBCmJg0zIQQQgghhBDCxKRhJoQQQgghhBAmJg0zUah5eXkxZ86cXLcfOnSI5s2bG167\nuLiwadOmggjN4JVXXiEiIqJAzymEMJ1Tp05hYWFBenp6rvt8/PHH/PjjjwBERUXh7OxcUOEBcPny\nZRo0aEBqamqBnlcIU3rttdeYOHFivp6jYcOGbNu27al93qxZsxg7dizwcHVLYSb1zpOThlk+27Fj\nBx4eHtjZ2WFvb0/r1q3Zt2+fqcMqMjQaDRqNJtftEydO5IMPPnjo/R/V4sWLGThwYJ77jB8/nk8+\n+eSpnVOIguDi4oK1tTVarZaqVasydOhQkpOTTR2WWbh69SoLFizAz88v387x8ssvs3Hjxly3V6lS\nhXbt2jF79ux8i0GIwuZp3wPk5ODBg7zwwgsATJkyhUGDBj32Z6WmpvL555/z4YcfPq3wHlg3PMit\nW7do1KgRNjY2REdHZ9v+4YcfUr16dcqVK4eTkxPjxo0jLS0NkHrnaZCGWT66ffs2Xbt2ZcyYMdy8\neZPz588zefJkSpcuberQzMLFixeJioqiZ8+e+XaOtWvX0qVLlzz3ad68Obdv3+aPP/7ItziEeNo0\nGg1hYWEkJiayf/9+9u3bh7+/v6nDMgvz5s2jS5cu+VbXJycn88cff9C2bds89xs4cCCzZs3KlxiE\nKKyUUqYO4aGtWrWK+vXr4+Dg8FQ+72HrhtzcvXuX7t2707hxY3755Rd69uzJ4cOHjfYZNmwYhw4d\n4vbt28TExLBhwwZ+/fVXw3apd56MNMzy0bFjx9BoNPTr1w+NRkOZMmXo0KEDrq6u2fa9e/cuVlZW\n3LhxA4DPP/8cS0tLkpKSAH3PUEZX99q1a2nSpAm2trZUr16dqVOnGj6nU6dO/Pzzz0af3ahRI1au\nXAnAkSNH6NChA/b29tSrV4+lS5fmGv+8efOoXbs25cqVo1atWgQFBRne9/T0ZPTo0djZ2VG/fn02\nb95sOO7WrVsMGzaMatWq4eTkxMSJE4265efOnUuDBg2oUKEC3t7enDlzxrAtMjKSevXqYWdnx+jR\no1FK5VrJRkZG0rRpU0qVKpXj9sOHD1OrVi1CQkIAfQ/Bt99+i5ubG1qtlmHDhnH58mU6deqEra0t\nHTp0ICEhwXB8eno6GzduxNvbm7t37/Lqq69SsWJFypcvj7u7O1euXDHs6+Xlxdq1a3P9WwpRmFWr\nVg1vb28OHjyYbdtvv/1G9+7dDa/r1q1L3759Da+dnZ2JjY0FYMyYMVSvXh1bW1uaNWvGjh07ALhw\n4QLW1tbcvHnTcNyBAweoVKkSOp0OyLteyCynsnj16lVAXw4//vhjWrRoga2tLT179jQ6Z3R0NB4e\nHpQvX57GjRuzdetWw7a86q309HTef/99KlWqRO3atR9Y1iMiIvK8MZo+fTrPPfccFy5cICoqCicn\nJ7755hsqV65MtWrVWLlyJeHh4TzzzDPY29vz1VdfGR2/adMmWrdujaWlJTExMTRr1gxbW1uqVq3K\ne++9Z9jP3d2df/75h7Nnz+YZrxBF1YEDB3j++ecpV64c/fv35+7du0bbw8LCaNy4MeXLl8fT05O4\nuDjDNhcXF6ZNm0ajRo2ws7Ojf//+3Lt3D4Br167RtWtXypcvj729vaGHLOO4TZs2ERERwZdffklI\nSAharZYmTZqwbNkymjVrZhTDd999l+sD5HXr1uVZVyxfvpyaNWty6NAhAObPn0+NGjWoWLEi/v7+\n2aZvZK4bpkyZQp8+fRg0aBDlypXDzc2Nv//+my+//JIqVapQo0YNIiMjDcfqdDp8fX2pX78+v//+\nO4MHD+bnn3+mW7dunDt3zrDfs88+i42NDaBvBFtYWBg1LKXeeUJK5Jvbt28re3t7NWTIELVu3Tp1\n48aNPPd/4YUX1PLly5VSSnXo0EHVqVNHrVu3TimlVJs2bdTKlSuVUkpFRUWpgwcPKqWUio2NVVWq\nVDFsmz9/vvL09DR8Znx8vLKzs1OpqakqKSlJOTk5qXnz5imdTqcOHDigKlasqA4dOpQtlqSkJFWu\nXDl17NgxpZRSly5dUvHx8UoppX777TdVsmRJ9cMPP6i0tDQVEhKibG1t1c2bN5VSSvXs2VP5+fmp\nlJQUdeXKFeXu7q5mzZqllFJq5cqVqk6dOurIkSNKp9Mpf39/5eHhoZRS6urVq0qr1arly5ertLQ0\n9f3336uSJUuqOXPm5Pj3ev/999Xbb79t9J6Li4vatGmT+uOPP1T16tXV2rVrjba1atVKXblyRZ0/\nf15VrlxZNWnSRP3555/q7t276sUXX1RTp0417L97925DbAEBAapbt27qzp07Kj09Xe3fv1/dvn3b\nsO93332nfHx8cvunFaLQcXFxURs3blRKKXXmzBn13HPPqUmTJmXb759//lF2dnZKKaXOnz+vatSo\noZydnZVSSp04cUKVL1/esO/ChQvVjRs3lE6nU9OmTVNVq1ZV9+7dU0op9eKLL6rAwEDDvu+//74a\nOXKkUirveiGrvMpi27ZtlaOjo4qPj1fJycmqd+/e6tVXX1VKKXXu3Dllb29vqFMjIyOVvb29unbt\nmlIq73pr5syZql69eurcuXPqxo0bysvLS1lYWCidTpdjjJUqVVL79u0zvN6yZYtycnJSSik1depU\n1bRpU8N5t2zZokqWLKk+++wzlZaWpgIDA5W9vb0aMGCASkpKUvHx8crKykqdPHnS8HkjRoxQs2fP\nVkop1bJlS7Vw4UKllFLJyckqOjraKBY3Nze1evXqHOMUoii7d++eql69uuFeZNmyZcrS0lJNnDhR\nKaXU/v37VeXKlVVMTIxKT09Xv//+u3JxcVGpqalKKX0d2KJFC3Xx4kV148YNVb9+fRUQEKCUUuqj\njz5Sfn5+Ki0tTaWlpakdO3YYzptxn6GUUlOmTFGDBg0yiqlChQrq8OHDhvcaN26sQkNDc7yG5s2b\nq2XLlhlenzx5Umk0GpWWlqbmzp2r6tSpo06cOKGU0t/P2djYqJ07d6rU1FT1/vvvK0tLS0MsShnX\nDZMnT1ZlypRRGzZsUGlpaWrw4MGqRo0a6osvvjDUNTVr1jQce+DAATVz5sxsMW7atEmFhIQYvffl\nl18qGxsbpdFo1Mcff5ztGKl3Hp80zPLZ4cOH1WuvvaacnJxUyZIlVffu3dXly5dz3HfixInqnXfe\nUWlpaapq1apq+vTp6qOPPlJ37txRVlZWuTbsxowZo8aOHauU0jcGy5Ytq86cOaOUUup///ufGjZs\nmFJKqeDgYNWmTRujY998802jxkiGpKQkZWdnp5YvX65SUlKMtv3222+qWrVqRu+5u7urBQsWqEuX\nLqnSpUurO3fuGLYFBQWpdu3aKaWU8vb2Nmpo6XQ6ZW1trU6fPq1+//131apVK6PPdXJyyrVhNnz4\ncPXRRx8Zvefi4qImTZqknJyc1NatW7NtCwoKMrzu3bu3euuttwyvZ8yYoXr27Gl4/cknnyh/f3+l\nlFJz585VHh4eKjY2NsdYZs+erV588cUctwlRGNWoUUPZ2NgoOzs7VaNGDTVq1Ch19+7dHPd1dnZW\n+/fvV4sXL1ZvvvmmatGihTpy5IiaO3eu6tGjR67nKF++vKHM/Prrr4Yykp6erpydndX27duVUrnX\nCxn1WGZ5lUUvLy+jm4RDhw6pUqVKKZ1Op7766iujGyillHr55ZfV77///sB6q127doZGmlJKbdiw\nQWk0mlwbZpaWluro0aOG11u2bFGOjo5q7Nixqk2bNkYPdbZs2aKsrKxUenq6Ukpfh2s0GhUTE2PY\np2nTpoaHb0opVb16dXXu3DmllP6B3uTJk9XVq1dzjMXT01MtWLAgx21CFGVbt27Ndi/i4eFhaJj5\n+fkZfs/w7LPPqm3btiml9PcEixYtMmz78MMPlZ+fn1JKqUmTJqkePXqo48ePZztv5obZ5MmTDQ9/\nMvj5+akJEyYopZQ6ePCgKl++vKExmFXdunXV+vXrDa8zGmbffPONatCggTp//rxh29SpU9WAAQMM\nr1NSUlSpUqWMGmaZ64bJkyerjh07GratXr1a2djYZKtrbt26lWNsD2P//v2qevXqhk6FDFLvPD4Z\nypjP6tWrx2+//cbZs2c5ePAgFy5c4N13381x37Zt2xIVFcX+/ftxdXWlffv2bN26lT179lCnTh3K\nly8PwJ49e2jXrh2VK1fGzs6OWbNmcf36dQC0Wi1dunRh8eLFAAQHBxuSV5w+fZo9e/ZQvnx5w09Q\nUBCXL1/OFkvZsmUJCQkhICCAatWq0bVrV44ePWrY7ujoaLR/jRo1uHDhAmfOnOH+/fs4ODgYzuHn\n52cYanT69GnGjBlj2GZvbw/A+fPnuXjxIk5OTkafm1cms/Lly5OYmGj0nlKKWbNm4enpaTT0IEOV\nKlUMv1tZWRm9LlOmjGHoKOiHGHTu3BmAQYMG8fLLL9O/f38cHR0ZP368YbIrQGJiInZ2drnGKkRh\no9FoWLVqFTdv3uTUqVP89NNPuc6Jyqibtm/fTtu2bWnbti1bt25l27ZtRsNwvv32Wxo0aICdnR3l\ny5fn1q1bXLt2DQAfHx92797NpUuX2LZtGxYWFrRu3RrIu17I6kFlMXOdUb16de7fv8+1a9c4ffo0\nS5cuNar/du7cyaVLlx5Yb128eDHb5+Ylp7opISGBX3/9lY8++gitVmu0zd7e3pCwwMrKCsheV2Uk\nZomLi8PW1tZQB8+ZM4djx45Rv3593N3dsw2zlLpJmKsLFy7keC+S4fTp00ybNs2ozJ87d44LFy4Y\n9qlatarhdysrK8M9wAcffECdOnXo2LEjtWvX5uuvv37ouIYMGWKY+rFgwQL69euHpaVljvuWL1+e\n27dvZ3t/2rRpjBo1imrVqhney3qPZGVlZagrIXvdAFC5cmWj/StWrJitrsl83/OomjRpwltvvcWC\nBQuM3pd65/FJw6wAPfvsswwZMiTHeRwArVq14ujRo6xYsQIvLy/q16/PmTNnCA8Px8vLy7DfgAED\n6NmzJ+fOnSMhIQE/Pz+jOVy+vr4sXryY3bt3c/fuXdq1awfobybatm3LzZs3DT+JiYnZ5qRl6Nix\nIxs2bODSpUvUq1eP4cOHG7ZlvWE6ffo0jo6OODs7U7p0aa5fv244x61btwzjuqtXr87s2bONYkhO\nTqZVq1Y4ODgYjUlWSuU5RtnNzY1jx44ZvafRaJg1axanT59m3LhxuR6b+Rw5uXTpEhcvXqRJkyYA\nlCxZkkmTJhEfH8+uXbsICwtj/vz5hv0PHz5M48aNH3g+IYqitm3bsmXLFrZv346Xl5ehobZ161ZD\nw2z79u188803LF26lISEBG7evImtra2hjJUvX56OHTsSEhJCUFAQvr6+hs/PrV5o2bJltlgeVBYz\nz007c+YMlpaWVKpUierVqzNo0KBs9d+HH36Ik5NTnvWWg4NDts/Ni5ubm9GDrIzrDwsLY+jQoeza\nteth//TZhIeHGyUkqlOnDkFBQVy9epXx48fzyiuvcOfOHQDS0tI4fvw4jRo1euzzCVFYOTg45Hgv\nkqF69epMmDDBqMwnJSXRr1+/HD8vczZHGxsbvv32W06cOMHq1av57rvv2LJlS57HZGjZsiWlSpVi\n27ZtLF68OM+sjTndxwBs2LABf39/QkNDja4381yvO3fuGB7KQ/a6Ib+zU2a4f/8+ZcuWNbyWeufJ\nSMMsHx09epTvvvvOUHGcPXuWxYsX06pVqxz3t7a2pmnTpvz888+Gmx0PDw8CAgKMnkonJSVRvnx5\nSpUqRUxMDEFBQUYFsHPnzpw+fZrJkyfTv39/w/tdu3bl2LFjLFy4kPv373P//n327t3LkSNHssVy\n5coVVq1aRXJyMpaWlpQtW5YSJUoYbZ8+fTr3799n6dKlHDlyhM6dO1O1alU6duzIuHHjSExMJD09\nnRMnThjW/PDz8+OLL74wTGS9deuWIQFJ586diY+PZ8WKFaSlpTF9+nQuXbqU69+3ffv27N+/P9t6\nGVqtloiICLZt28bHH3+c6/F5WbduHZ06dTK8joqKIi4uDp1Oh1arxdLS0ujvsW3bNqP9hTAnGQ2z\nu3fvUq1aNVq3bk1ERAQ3btwwPLxITEykZMmSVKxYkdTUVD799NNsT4IHDBjA77//zvLlyxkwYIDh\n/bzqhazyKotKKRYuXMjhw4dJSUlh0qRJ9OnTB41Gw6uvvsqaNWvYsGEDOp2Ou3fvEhUVxfnz53Fw\ncMiz3urbty/Tp0/n/Pnz3Lx5M1syjqw6d+5slFgkwwsvvMCiRYvw8fFh7969D/nXN7Zu3Tqjm6+F\nCxcaevZsbW3RaDRYWOi/2mNiYnBxcSnwNdSEKAgeHh6ULFnScC8SGhpqVK6GDx9OQEAAMTExKKVI\nTk5m7dq1ufYQZX5QGxYWxvHjx1FKUa5cOUqUKGEoV5lVrVqVU6dOZXvIO2jQIN5++21KlSqFh4dH\nrteQW13x3HPPERERwahRo1izZg2gXzN1zZo17N69m9TUVKZMmWJ03qx1Q24Pnp9ExqikhIQElFLE\nxMTwyy+/4OPjY9hH6p0nIw2zfKTVatmzZw8tWrTAxsaGVq1a4ebmxrRp03I9pm3btqSlpeHu7m54\nnZSUZDQs75dffmHSpEmUK1eOzz77LNvTn1KlSuHj48OmTZuMbn5sbGzYsGEDwcHBODo64uDgwMcf\nf5zjQoDp6el8//33ODo6Ym9vz/bt25k5c6Zhe4sWLfj777+pVKkSEydOZPny5YahlvPnzyc1NdWQ\nYa1Pnz6GBlbPnj0ZP348/fv3x9bWFldXV9avXw9AxYoVWbp0KR999BEVK1bk+PHjhqFOOalSpQov\nvviiIeNkZra2tkRGRrJu3TomT56c62dkbtBmXv9k7dq1hmGMoO9B69OnD7a2tjRo0AAvLy/DU7C9\ne/ei1WqzZWISwlzUrVsXrVZLmzZtAChXrhy1a9fG09PTUGa8vb3x9vbmmWeewcXFBSsrq2xD/rp3\n787x48dxcHAwyk6bV72QVV5lUaPRMGjQIF577TUcHBxITU1l+vTpADg5ObFq1Sq++OILKleuTPXq\n1Zk2bZphtEFe9dbw4cN5+eWXadSoEc2aNaN37955Po0ePHgw4eHhRhniMvZv3749c+fOpVu3bvz5\n559G27Lum1VCQgKHDh0yutFbv349DRs2RKvVMnbsWIKDgw1DUhctWsTIkSNzjVOIoszS0pLQ0FDm\nzZuHvb09S5YsoXfv3obtTZs2JTAwkLfffpsKFSpQt25d5s+fn2v5ynwPcPz4cTp06IBWq8XDw4NR\no0blmD2xT58+gH44cuZ7gEGDBhEfH8+rr76a5zV07dqVI0eOcPHiRaM4QN+bFhYWxvDhw1m/fj0N\nGjRgxowZ9O/fn2rVqqHVaqlcuTKlS5fOsW7IaU23h61r8rJy5Upq166Nra0tw4YNw9/f36hhJvXO\nk9Go/GhSC7M2b9485syZw/bt200dCocPH2bIkCHExMQ8tc9MS0vDwcGBkydPGlLC5uWVV17hjTfe\nwNvb+6nFIIR4dO3atWPQoEG8/vrrpg6FCRMmULlyZcaMGfPUPnPJkiWEhoYSHBz8wH2vXLmCl5cX\nf/75Z65Liggh8sedO3eoUqUKBw4coHbt2nnuGxgYyKFDh/j+++8f6RwZo6eOHz/Onj17HrpuyE9S\n7zy5B/aYffnllzz33HO4uroyYMAAwxoPmb3zzjvUrVuXRo0aceDAgXwJVIic1K9f/6k2ygBu3ryJ\nv7//QzXKAJYtWyaNMhN7/fXXqVKlSrY1AmfMmEH9+vVp2LAh48ePN1F0oiAVlmeNn3/++VNtlIF+\nnlrGepYPUrlyZQ4dOiQ3Rybi4uKCm5sbTZo0MYyAySoqKoomTZrQsGFDo3nkouibOXMm7u7uD2yU\ngb5H/mEbZWvWrCElJYXk5GTef/993NzcqFGjxiPVDflJ6p0nVzKvjadOnSIwMJDDhw9TunRp+vXr\nR3BwMEOGDDHsEx4ezvHjx/n777/Zs2cPI0eOJDo6Ot8DF6aTU/e4OalUqRIjRowwdRjiEQwdOpTR\no0czePBgw3tbtmxh9erVxMbGYmlpaZiHI8ybOddNHTp0MHUI4iFpNBqioqKoUKFCjtsTEhIYNWoU\n69evx8nJyZC9VBR9Li4uaDSaHKdZPKnVq1czePBglFI0b97c0EMmdYP5yLNhVq5cOSwtLUlJSaFE\niRKkpKRkS026evVqQ0OtRYsWJCQkcPnyZaNUv8K8DBkyxKhxLoSptWnThlOnThm9N3PmTD7++GND\nmuJKlSqZIDJRkHLKmiaEqeTVexsUFETv3r0N6c8rVqxYUGGJfJb1u+hpCgwMJDAwMN8+X5henkMZ\nK1SowHvvvUf16tWpVq0adnZ2tG/f3mif8+fPG2VecXJyMkrnKYQQpvD333+zbds2WrZsiZeXF/v2\n7TN1SEKIYkKj0dC+fXuaNWuW443033//zY0bN2jXrh3NmjXLtg6UEKJ4yrPH7MSJE/zwww+cOnUK\nW1tb+vTpw6JFiwwLFmfI+lQop6Ek5jy8RIjirLDM6ckqLS2NmzdvEh0dzd69e+nbty///PNPtv2k\nbhLCPJmybtq5cycODg5cvXqVDh06UK9ePUNWU9Cv/bR//342bdpESkoKrVq1omXLltStW9foc6R+\nEsI85VY/5dljtm/fPjw8PLC3t6dkyZL4+PhkWxjT0dHRaBHgc+fOZRvumDmIgv6ZPHmySc4r8Rf9\n+Ity7AUVf2Hm5ORkSOHbvHlzLCwsjBbjzMxc/30kfom/sP7kd/ym5uDgAOiHUPfq1StbkipnZ2c6\nduyIlZUV9vb2vPDCC/z11185fpY5/vtI/OYZu8T/cD95ybNhVq9ePaKjo7lz5w5KKTZu3EiDBg2M\n9unevTvz588HIDo6Gjs7O5lfJoQwuZ49e7J582YAjh07RmpqKvb29iaOSghh7lJSUkhMTAQgOTmZ\nDRs2ZMsY26NHD3bs2IFOpyMlJYU9e/Zku78SQhQ/eQ5lbNSoEYMHD6ZZs2ZYWFjw/PPPM3z4cGbN\nmgXAiBEj6Ny5M+Hh4dSpU4eyZcvy22+/FUjgQgiRwdfXl61bt3L9+nWcnZ359NNPef3113n99ddx\ndXWlVKlShgdIQgiRny5fvkyvXr0A/ZDqgQMH0rFjR6N7p3r16uHt7Y2bmxsWFhYMHz5cGmZCiIJb\nYFqj0Tyw+y4/REVFFen1QSR+0ynKsUPBxG+qcv00Sd30eCR+05L482YOdRNI/fS4inL8RTl2kPgf\nRl7l2uwbZkKI/GMO5docrkEIYcxcyrW5XIcQ4j95les855gJIYQQQgghhMh/0jATQgghhBBCCBOT\nhpkQQgghhBBCmJg0zIQQQgghhBDCxKRhJoQQQgghhBAmJg0zIYQQQgghhDAxaZgJIYQQQgghhImV\nNHUAQgghhBBCCGHOlFLExcXluY/0mAkhhBBCCCFEPlFK0bdvXzw9PfPcT6MKaEl5Wb1eCPNjDuXa\nHK5BCGHMXMq1uVyHEMVdbGwsHh4eJCcnA+RarqXHTAghhBBCCCHygVKKFStWGBpleZEeMyHEYzOH\ncm0O1yCEMGYu5dpcrkOI4iohIYERI0YQHx+Po6Mju3fvJjExUXrMhBBCCCGEEKIg7Nixg8aNG1O5\ncmX27t1LREQEO3bsyPMY6TETQjw2cyjX5nANQghj5lKuzeU6hChO0tLS8Pf3JyAggMDAQLp162a0\nPa9yLenyhRBCCCGEEOIJnT59moEDB2JlZcX+/fupVq3aIx0vQxmFEEIIIYQQ4gksWbKE5s2b06NH\nD9avX//IjTKQHjMhhBBCCCGEeCxJSUmMGTOGbdu2ER4eTrNmzR77s6THTAghhBBCCCEeQClFbGws\nsbGxKKXYv38/TZs2JT09nf379z9Rowykx0wIIYQQQggh8qSUom/fvkRERABQq1YtLly4wIwZM+jf\nv/9TOYdkZRRCPDZzKNfmcA1CCGPmUq7N5TqEMAexsbF4enqSlJQEgIWFBWFhYXTq1OmRPievci1D\nGYUQQgghhBDiAXQ6neH3smXL4ujo+FQ/X4YyCiGEEEIIIUQu7t+/z8KFC1FKYW1tTYkSJfD29sbV\n1fWpnkcaZkIIIYQQQgiRg5MnT+Lr64u9vT2nT5/m0qVLALi6uqLRaJ7quWQooxBCCCGEEEJksWzZ\nMlq0aEHfvn1Zs2YNlStXxs3NDTc3t6feKAPpMRNCCCGEEEIIgzt37jBu3Dg2bNjA2rVrad68eYGc\nV3rMhBBCCCGEEAI4dOgQ7u7uJCQkcODAgQJrlIE0zIQQQgghhBDFnFKKuXPn0rZtW959912CgoIo\nV65cgcYgQxmFEEIIIYQQxdbt27fx8/MjLi6OrVu30qBBA5PEIT1mQgghhBBCiGJBKUVsbCyxsbEo\npdi3bx/PP/885cqVIyYmxmSNMpAeMyGEEEIIIUQxoJSib9++REREAFC7dm3Onz/PL7/8Qp8+fUwc\nHWiUUqpATqTRUECnEkIUEHMo1+ZwDUIIY6Yu1y4uLpQrV44SJUpgaWlJTExMjvvt3buXVq1asWTJ\nEnx8fLJtN/V1CGFuYmNj8fT0JCkpCQALCwvCwsLo1KlTgcWQV7mWHjMhhBBCiKdIo9EQFRVFhQoV\nct1Hp9Mxfvx4vL29pfElRAFKT083/F62bFkcHR1NGI0xmWMmhBBCCPGUPaixNWPGDF555RUqVapU\nQBEJUbyFNLEGAAAgAElEQVSlp6ezbt06dDodVlZWaLVavL29cXV1NXVoBtJjJoQQQgjxFGk0Gtq3\nb0+JEiUYMWIEw4cPN9p+/vx5Vq1axebNm9m7dy8ajcZEkQpRPFy7do0hQ4Zw8+ZNjh49yq1btwBw\ndXUtVOVPGmZCCCGEEE/Rzp07cXBw4OrVq3To0IF69erRpk0bw/Z3332Xr776yjDXJK/etSlTphh+\n9/LywsvLKx8jF8L87Ny5E19fX3x9ffH398fS0rJAzx8VFUVUVNRD7SvJP4QQj80cyrU5XIMQwlhh\nKtdTp07FxsaG9957z/BerVq1DPFdu3YNa2trAgMD6d69u9Gxhek6hChq0tPT+eabb/j++++ZM2cO\nXbp0MXVIgCT/EEIIIYQoECkpKeh0OrRaLcnJyWzYsIHJkycb7fPPP/8Yfh86dCjdunXL1igTQjy+\na9euMXjwYG7dusXevXtxdnY2dUgP5YHJP44ePUqTJk0MP7a2tkyfPt1on6ioKGxtbQ37+Pv751vA\nQgiR1euvv06VKlWMJvBOmTIFJycnQ72UsWaJEELkp8uXL9OmTRsaN25MixYt6Nq1Kx07dmTWrFnM\nmjXL1OEJYfZ27NjB888/j6urK1FRUUWmUQaPOJQxPT0dR0dHYmJijC4yKiqK7777jtWrV+d+IumO\nF8LsFJZyvX37dmxsbBg8eDBxcXGAfviQVqtl3LhxeR5bWK5BCPH0mEu5NpfrEKIgpKen83//93/8\n8MMPzJ07l86dO5s6pBw9taGMGzdupHbt2jm2PKXiEEKYSps2bTh16lS296VeEkIIIcyTUsrwMNbB\nwYHBgwdz+/btIjV0MatHapgFBwczYMCAbO9rNBp27dpFo0aNcHR05Ntvv6VBgwbZ9pPMQkIUbY+S\nWagwmDFjBvPnz6dZs2ZMmzYNOzu7HPeTukmIoq2o1U1CiCejlKJv375ERESg0+kAGD16tEmyLj5N\nDz2UMTU1FUdHRw4dOpRtMcTExERKlCiBtbU169atY8yYMRw7dsz4RNIdL4TZKUzl+tSpU3Tr1s3w\n9OzKlSuGumrixIlcvHiROXPmZDuuMF2DEOLpMJdybS7XIcTTFhsbi6enJ0lJSQBYWVkRHR2Nm5ub\niSN7sLzK9QOTf2RYt24dTZs2zXGFeq1Wi7W1NQCdOnXi/v373Lhx4zHDFUKIJ1e5cmU0Gg0ajYY3\n3niDmJgYU4ckhBBCiKfg9u3b3Llzx/C6ZEnzSDT/0A2zxYsX4+vrm+O2y5cvG1p+MTExKKWoUKHC\n04lQCCEew8WLFw2/r1ixwihjoxBCCCGKpn379jFo0CBq1aqFVqtFq9Xi7e1tFt/zD9W8TE5OZuPG\njQQGBhrey0j5OmLECJYtW8bMmTMpWbIk1tbWBAcH50+0QgiRA19fX7Zu3cq1a9dwdnZm6tSpREVF\n8eeff6LRaKhZs6akqRZCCCGKMKUUAQEBTJ48mZkzZ+Lj42OYvuDq6opGozFxhE/ukdLlP9GJZJy0\nEGbHHMq1OVyDEMKYuZRrc7kOIZ5UUlISfn5+xMXFsXTpUp555hlTh/TYnsocMyGEEEIIIYQoSIcP\nH8bd3Z1SpUqxe/fuIt0oexBpmAkhHplSitjYWFOHIYQQQggzFhQUxAsvvMD777/P3LlzDckGzZV5\npDARQhQY/dohE4iIaGjqUIQQQghhhu7du8fYsWOJjIxk48aNNGrUyNQhPbHMC2LnRnrMhBCP5MCB\nOMLDG5KUlH2xeSGEEEKIR5ExCic2NhalFKdOnaJ169ZcuXKFffv2mU2jrG/fCXh6HsxzP+kxE0Lk\n6e5d2LsXtm+Hbdtgxw7ItHSIEEIIIcRj0TdY+hIREQFA48aNOXr0KP/73/8YM2aMWWRaBIiLiyMi\nIuOh9sBc95MeMyGKqaxPqDIkJsL69TBhArzwAlSsCO+9B9evg58fnDzpSu/eB9Fqg0wYvRBCCCGK\nOn2DJYKkpCSSkpLYuXMn3377Le+++67ZNMoA9u+H5OQH7yc9ZkIUQ1nnibm6huDu7s+OHRqOHIFm\nzaBNG5g4EVq1AhubzEdrWLLkc+Li4jCD0QVCCCGEKCTKli1L48aNTR3GU5OcrH/QHRLiSsuWwRw8\nGERiYu77S8NMiGLmyhWYPz+OVasacv++fp5YTEwQzZvH8eOPbjRrBqVL5/0ZGo0GNze3AohWCCGE\nEOYqKSkJnU5HqVKlKF26NN7e3ri6upo6rKdi61YYNkz/gPvgQQ0VKjz4obY0zIQo4jJn+XF1zb7y\n/dWr+sohKkr/c+4cNGoEFpkGMltb6ysPaWsJIYQQIr8ppfj555/59NNPCQ4OxsXFBcj5PqaoSUqC\njz6ClSth5kzo1i1jy4MfakvDTIgiLOuQxE6dQvjpJ3+2b9ewZct/DbHWrcHLC37/HRo3BgsLV/r2\nDWb9ev08MW/veFxdfU13IUIIIYQoFpKTkxkxYgQHDx5k9+7d1K5d29QhPTWbN8Mbb+jn6MfFQfny\nj3a8RmWe9Z+PNBoNBXQqIYqN2NhYPD0PGlLXazRBWFs3pG1bN9q10zfGGjeGkjk8gnlQT9vDMIdy\nbQ7XIIQwZi7l2lyuQ4gMx48fx8fHh8aNGxMQEGA2C0YnJsL48bBmDQQEQJcuue+bV7mWHjMhipiU\nFH3K+k2b9BVAUtJ/26yt9Sntn3/+wZ8j88SEEEIIUVBWr17NG2+8wdSpU/Hz8yvyQxYzbNyo7yV7\n6SV9L5md3eN/ljTMhCgE8uq9SkvTryO2aZO+8O/bB02a6CuAgABXfvwxmMjI/4YkNmkiQxKFEEII\nUTjodDomTZrEggULWL16NS1btjR1SI8t8/1ajRqufPCBhogImD0bvL2f/PNlKKMQJpZ9nlg8Eyf6\ns3mzho0b9Qs7u7joG2IvvaQft5w5ff3TGJL4ODHrMws1KvLlWuomIcyPuZRrc7kOUXxdu3YNX19f\ndDodwcHBVK5c2dQhPbbM92s6HUA8r77qz7ffaihX7uE/J69yLQ0zIUwsNjYWD4+DJCcP+PedIBwd\nG9K5sxsvvQQvvgiVKpk0RCP6iqmvYUHIol6upW4SwvyYS7k2l+sQxUPWB8V79+6lT58++Pr64u/v\nT8mcJrwXIVnn9VtZBREd3fCRp4XIHDMhCpmUFP1csA0bYPVq49XgbWwgPDz/U9c/bk/brl27CA8P\nJyUlJT/DE0IIIUQRkfmhLcAzzzzDmTNnmD17Nr169TJxdE/Hnj3G92v50c6UHjMhnqLcGjvp6RAb\nq2+IbdigL9zPPw8dO0KHDor/+78JbNigH8ro7R1PSIh/vg5JzFqBdurUiZCQEEM5vXjxIidOnMjx\nJyUlhbt37xrKc1Ev11I3CWF+zKVcm8t1CPOn703yJOnfjGQajYZVq1bR7b9FvIqs5GR9xsWVKxU1\na07gr7+e7H5NhjIKUQCyzhXz8ornlVf8iYzUEBmpz9LTsaP+x8sLtFrjYwtynljWCtTS0hIPDw+u\nXbvGyZMnsbGxoXbt2jn+VKpUiX79+rF+/XoSExOLfLmWukkI82Mu5dpcrkOYP/20DA+S/+1SsrGx\nYefOnUU++/OOHfDaa+DhAT/+CHZ2+bvUkAxlFOIp2b8/jrVrG3Lnjn7scVhYEHfvxtG3rxv+/voE\nHrkpqNT1aWlpbN++nV9//dVQeQJYWFjQvXt3XnrpJWrVqoU2c6sxB0uWLDEk/xBCCCFE8Xbt2jXS\n0tIoXbo0pUqVwtvbG1dXV1OH9dju3oVPPoGgIJg5E3r0yNiSv/dr0jAT4gmcOwfr1ul/IiPh3r3/\ntmm1MG1a/s8Ve5Dk5GTWr1/PypUrCQ8Pp2bNmvTo0YMOHTqwe/duALy9vRk7duxDP/mRNdCEEEII\noZTihx9+4Ouvv2bNmjVUqVIFKLgs0flh714YMgQaNtRPQ6lYseDOLUMZhchFTsML79+HXbv0yTnW\nrYMLF/RDEzt1go4dFW+/PYH16wturlhucV65coWwsDBWrlxJVFQULVu2pGfPnnTv3h0nJ6dcj3tU\n5lCuzeEahBDGzKVcm8t1CPOUkpLC8OHDOXToECtWrMAlr6FBRUBqKnz6KQQGwvTp0K9f/pxH5pgJ\n8YgyzxdTCurXj8fZWb+2WJ060LmzvjHm7g4lShgfV5BzxTIn8VBKUbt2bbRaLQcPHqRjx4707NmT\nzp07Y/cky9DnwRzKtTlcgxDCmLmUa3O5DmF+Tp48Sa9evXB1dWXWrFlYW1ubOqQn8tdf+l6y6tX1\ni0VXrZp/55KGmRCPQKeDxYtjGTbsIKmp+vliJUsGMWVKQ954w41/e+lNTilFSEgIQ4YMITU1FdAn\n8fjhhx8YNmwYpUuXzvcYzKFcm8M1CCGMmUu5NpfrEOYlMjKSQYMG8b///Y/Ro0cXuSGLmR+i16/v\nyv/9n4YffoBvv4XBgyG/L0eSfwjxAImJ+jT2YWH6YYpaLWQuM1ZW0K0bJm+UZSTvWLFiBStXrgSM\n09WXKVOG1q1bF0ijTAghhBDFh1KKb775hu+//57g4GC8vLxMHdIjyzoiqlSpEJo29Wf/fg3OzqaO\nDixMHYAQ+U0pRWxsLLGxsUaNmFOnYMYM/RyxatVg1iz92mLR0fD336706HEQrTYIrTYIb+/4fM0u\nlFuMoB/DvXLlSl577TWqVq3KBx98QJUqVVi3bh2nTp2iR48eaLVatFptkc+CJIQQ5sDFxQU3Nzea\nNGmCu7t7tu2LFi2iUaNGuLm54enpSWxsrAmiFOLhJSUl0b9/f5YuXUpMTEyRbJQBxMXFERHRkKSk\nASQnD+DOnef49tu4QtEoA+kxE2Yu69pizZuH4O7uT1iYhitXoEsX8POD5cuN1xUDDUuWfJ5pvphv\nvnXV57TY88yZM1m7di0rVqxg8+bNNGvWjJ49e/LZZ5/hnKX2yEhdr4+z6GZBEkIIc6HRaIiKiqJC\nhQo5bq9Vqxbbtm3D1taWiIgI3nzzTaKjows4SiFyl3m4n5WVFT4+PjRv3pzt27dTpkwZE0f3+C5e\nhDt3/nttaZn/QxcfhcwxE2Ztz55Y2rU7aFhbTKMJYtiwhgwb5oa7O1gUgj7jrIs9lyhRgtKlS9Ox\nY0d69epFly5dsLe3N3GUOTOHcm0O1yCEMGbqcl2zZk327dv3UHX3zZs3cXV15dy5c9m2mfo6RPGU\n+YGxTqcjPT2d7777jpEjRxbZh79KwcKFMHaswsFhAqdPF2wG7cxkjpkoVq5fhzVrYOVK/dpi/+bF\nAMDGBkaPNv3aYhmOHj3KnDlzSElJMbxXqlQpNm/eTIsWLUwYmRBCiMel0Who3749JUqUYMSIEQwf\nPjzXfefMmUPnzp0LMDoh8qYf7hdheGBsbW1N69ati2yj7No1/eioI0dg40YNjRoVzIioxyENM1Gk\n5JaO/vRpfUNs5UrYvx/at4fevWHOHFf8/IJZvz4I4N+5Yr4mjf/AgQOsWLGC0NBQEhIS6NGjB61b\nt+bAgQP/xuid45wEIYQQRcPOnTtxcHDg6tWrdOjQgXr16tGmTZts+23ZsoW5c+eyc+fOXD9rypQp\nht+9vLyK7NweUXTcuXOHO5nG+5XIvC5QERMeDsOHg6+vvsdMPwpTg1sBPqGPiooiKirqofaVoYyi\nyMg6X8zDIx4PD39WrdJw9qw+a2KvXvpGmZWV8XEFNQcrp3PpdDp27dpFaGgoK1aswNLSEh8fH3r1\n6oW7uzsWFhYFvv7Z02IO5docrkEIoZdRlzZq1KjQlOupU6diY2PDe++9Z/R+bGwsPj4+REREUKdO\nnRyPlfpJFLSMpGK3bt3ixo0bgP6BcUhISJG5NwFISoL334eICJg3DwrT8wwZyijMwl9/xREe3pCU\nFP18scjIICpXjuOHH9zw8ICSufxv1mgK5slI1iQeTZo04dlnn2X16tU4ODjg4+NDWFgYzz33XLbK\nraBiFEIIc5X14Z2ppKSkoNPp0Gq1JCcns2HDBiZPnmy0z5kzZ/Dx8WHhwoW5NsqEKGhbt26lf//+\njB8/nnfeeYeDBw8CReuBMcDu3fr1yDw99QtH29qaOqKHJw0zUajpdLBzpz5rYnCwcSYdGxv44IPC\nM18sJiaGtWvXGrr/d+3aRfPmzdm9eze1atUycXRCCGHeMqfBhoEmi+Py5cv06tUL0K89OXDgQDp2\n7MisWbMAGDFiBJ9++ik3b95k5MiRAFhaWhITE2OymEXxppQiICCAKVOmsHDhQjp06ABQ6B8YZx1t\ndP++hk8/hV9/hV9+AR8fEwf4GGQoozCJvIbupaXBtm2wbBmEhkLVqvDKK+Djo5g8eQLr15suk05W\nt2/fJjw8nOXLl7Nu3Tru3r2LTqcDQKvVsmPHjkJfsT0JcyjX5nANQhR3qanw00+xfPjhQXS6AYB5\nlGupn0R+S01NZfTo0ezYsYNVq1YVmR7crD3knp7xXL7sT7VqGubM0d87FlZ5lWtpmIkCl7UwdeoU\nz8KF/kRFaVi2TJ/Ao3p1fWOsd2+oW9f4WFPOFwN9auPVq1ezfPlyoqKiaNOmDb1796Zbt2689dZb\nrF+/HiiaY7IflTmUa3O4BiGKq8OHYc4cWLAAnnlGcffuBI4caUhS0kCzKNdSP4n8dPnyZXr37k3F\nihVZsGABWuMFXQs1/VJDB//tIQcIYuLEhkyd6lao1iXLicwxE4WK8XATWLkyiEqV4qhf341XXoE9\ne6BmzZyPNdV8sXbt2tGlSxdCQ0PZvXs3L730En379mX+/PnY2dkZjpPFnoUQIn8lJsKSJfoG2cmT\nMGQIbN8OzzyjQanP/03+YeoohSjc/vjjD3r16sXQoUOZPHkyFoVhYddHlLltU7as/oF+Ub/tkoaZ\nKFBpaRAdDffu/feehYX+S/bll00XV1ZxcXGsW7eO5ORkANasWUNqaipvvPEGy5cvx8bGJsfjJImH\nEEI8uawjFkDD7t36xlhoKLRtCx99BJ07Gyd+kjpYiAdbvHgx77zzDjNnzuSVV14xdTiP5dgxV9LS\ngildOohSpUy/HNLTIkMZRb5LT4ddu/TJO5YuBScnRVraBP75pyEaTeGYK5bh3LlzLF++nHnz5vHn\nn38a3i8O88UehzmUa3O4BiHMSdbh7rVqxXPvnj9KaRg2TJ9t7UHzR8ylXJvLdQjTyfyQo0GDBnzy\nyScsWbKElStXFsl7msREGDNGn4tg4UKFtXXRG6WUV7nOs9/y6NGjNGnSxPBja2vL9OnTs+33zjvv\nULduXRo1amRYJFcUH0opYmNjiY2NNfxHUwr27dOvIVGjBowcCQ4O+gyLf/yh4c8/P2fnzobs2NHQ\n5I2ys2fP8sMPP+Dp6YmbmxsHDhzgs88+w8fHB61Wi1arxdvb+9+ntqIwev3116lSpUqO/0bTpk3D\nwsLCsB6LEKJw+/PPONau1Q93T0oawOHDz/HRR3EcOQIffli4J/ULUZhkTMvw9PTE09MTZ2dn9uzZ\nQ0xMTJFslEVHQ5Mm+uGKBw5Ay5b6HnI3N7ci0yh7kDyHMj777LOGhlZ6ejqOjo6GFLAZwsPDOX78\nOH///Td79uxh5MiRREdH51/EolDJvuhzCM2a+RMSokGjgf79Yd06aJhlWZmCHG6SUxKPM2fOsGzZ\nMpYuXcqxY8fo3r07EyZMoH379pQqVQqALl26yHyxImLo0KGMHj2awYMHG71/9uxZIiMjqVGjhoki\nE0I8rPPnYe5cfZrrzMPdy5SB558v+nNHhCho+jn9ESQlJQFw7949pk2bRsWKFU0c2aNJS4MvvoCf\nf4aZM4tmGvyH9dBzzDZu3Ejt2rVxdnY2en/16tUMGTIEgBYtWpCQkMDly5epUqXK041UFEr6uVgN\nSU7+b9FnB4c4QkLcCsUXaeYkHkopateuTZkyZThx4gQ9evRg8uTJvPjii4bGWGYyV6HoaNOmDadO\nncr2/rhx4/i///s/evToUfBBCSEeSKfTP7ybPRt27IB+/SAszJWvvgpm/fogwHzmjghhChlL+ACU\nKVOGkiWLVnqJf/6BQYPA2hr27wdHR1NHlL8e+l8nODiYAQMGZHv//PnzRo01Jycnzp07Jw0zM5eQ\noF9nLCAA/s2PAegXfR43rvAs+rxx40bWrFnDvX8fvx4+fJgff/yRN954A0tLSxNHJ/LTqlWrcHJy\nksa1EIXQmTP63rE5c/Q3Wm++CUFB+u8Q0LBkyeeZRiz4yogFIR6RUoqtW7eSnp6OtbU1JUqUKFLT\nMpTSL4Px3nvwv//p55UVwcSRj+yhGmapqamsWbOGr7/+OsftWSew5VaBTpkyxfC7l5cXXl5eDxel\nKFC5rd+Vmgrh4bBwIURGQocO8MknrixYEExkZOF5snnp0iWWLVtGcHAwsbGx2Z4WeXp6SqPsMUVF\nRREVFWXqMB4oJSWFL774gsjISMN7eU2gl7pJiKcv63eJTqdh7Vp979ju3TBgAISFkWNq+0cdsVBU\n6iYhCsL9+/d555132LZtG/Hx8YYM04V1WkbWuiIhQYOfH8THw8aNOdcR5uqhsjKuWrWKmTNnGtZ0\nyszPzw8vLy/69+8PQL169di6dWu2HjPJLFQ05LT48zvv+LNokYalS/VzxV59Vb9WRMbyXYVh0edr\n166xfPlyQkJCOHDgAF27dqVfv3506NCBV199tVgt+lyQClO5PnXqFN26dSMuLo64uDjat2+PtbU1\noM+26ejoSExMDJUrVzY6rjBdgxDmIvN3iVJQvXo8CQn+1Kyp4c03oU8f/dCk/GIu5dpcrkMUnOvX\nr9OnTx+srKxYvHgx5cqVM3VIecp639m0aTwnTvjj46Phq6/AysrEAeaDJ15gevHixfj65twL0r17\nd3766Sf69+9PdHQ0dnZ2MoyxCMu6+POyZUHs2xfH8OFu/PGHPsNiVqZa9Pmll16iW7duLFmyhD17\n9uDt7c3o0aPx9vbGKlNJlkWfix9XV1cuX75seF2zZk3++OMPKlSoYMKohCg+MjIr3rmj/y45fjyI\n4OA4fHxkaLEQ+eXw4cN0796dHj168PXXX1OiRAlTh/RAWe87t20L4uef4xg5snjWFQ9smCUnJ7Nx\n40YCAwMN782aNQuAESNG0LlzZ8LDw6lTpw5ly5blt99+y79oRb66fRuWL4eUlP/es7aGFSsKRzdy\n1uxCq1atIikpiTfffJPQ0FDKli2b43GSxMP8+fr6snXrVq5fv46zszOffvopQ4cONWyXxrgQBePS\nJfj1V/jpp+yZFevUMV1cQpi7iIgIBg8ezNdff230/VcUZO48KlsWPD1NF4upyQLTxVx6OkRFwW+/\nwZo18OKLiitXJhAbq+9SLgyLP9+7d4+IiAhmzpxpGJIIYGNjw86dO6XRZULmUK7N4RqEMCWl9N8j\nM2fq5x/36QN+foovv5zA+vWm+S4xl3JtLtch8o9SiunTp/PVV1+xdOlSWrdubeqQHppSMG+ews9v\nAhpNQ0qVKhz3nfktr3ItDTMzl9ucrJMn4fff9T+2tjB0qH4idqVKBTtnLDc6nY5t27YRFBREaGgo\nDRs2xNfXl3Xr1rFlyxZA5osVBuZQrs3hGoTIT7l9J9y8qf8OCQiAEiVg5Eh9Wmtb27yPKwjmUq7N\n5TpE/khNTWXUqFFER0ezZs0aXFxcTB3SQ7t1C/z8IC4OgoIUUHymnEjDrJjKOqGyY8d4evTwZ948\nDXFx4Ourb5A1aWKa2LJ+YSul2L9/P0FBQQQHB1OlShUGDBhAv379DEsyFIZGo/iPOZRrc7gGIfJL\nTgmh3n/fn4AADaGh0KmTvkHWpo3p163MzFzKtblch3g6Mt8DOTg40KdPH8qVK8eiRYvQarUmju7h\nRUfrOwO8vWHaNPNM8JEXaZgVU7GxsXh6HjRMqIQgWrduyJgxbnTrBqVLmyaurEk8WrduTYsWLVi8\neDE6nY4BAwbg6+tL/fr1TROgeGjmUK7N4RqEyC9Zv0csLIKoUqUho0e78frrUFhzfZlLuTaX6xBP\nLvO9U8b/iVGjRvHFF18UiSQfoF9Q/uuv4ccfYdYs6NnT1BGZxhNnZRRFz82b+sU6sy7+/PPPpl/8\nOWsSj/Xr12Nvb8+CBQto3ry59IIJIUQhcfascRKP0qVh7VrTjLQQojjLeu9UpkwZBg4cWGQaZefP\n64c6p6fDH3+Ak5OpIyqcisEa2sWHUrBzJwwZAjVrwqlTrrRpcxCtNgitNohOneJNuuJ7SkoKixcv\n5q233jJULKBP4vHhhx/i7u4ujTIhhDCx9HRYtw66dIHBg12pWfMgNjb675GuXeNp3Nh03yNCFGf3\n7983/G5paWnCSB7N6tXQtCm8+CJs2iSNsrxIj1kRktv8qhs3YP58mD1b/4X65pv6MbsVK2pQ6vNM\nx/jmS8Mnr3lf6enpREVFsWDBAlauXIm7uzsjRoygcuXKbNy4EdAn8TBlg1EIIQQkJOgz9P7yi36E\nxejRsHSpBiur/P8eEULkLi0tjcDAQEqVKkXJkiWxsLAolPdOWe8H797V8MEH+l720FDw8DBxgEWA\nzDErInKagD1qlD+BgRrCwqBrV32DrKAnYGedL9apUydCQkI4fPgwCxYsYOHChdjb2zNo0CAGDBiA\ng4OD4ThJ4lH0mUO5NodrEOJh5FbvxsXph7mHhOiTebz9NrRqVbiSeTysjGts1KiRWZRrqZ9EYmIi\n/fv35/79+4SEhHD27Fmg8N07Zb1P9fSM5/x5f557TkNAANjZmTjAQkSSf5iBrBOwNZogXFwa8s47\nbgwaBPb2pozL0zA0sXTp0tSqVYtbt24xcOBABg0aVOie6IinxxzKtTlcgxAPkvWm6eWX4+nXz5+f\nf9bw998wYgQMHw7/PjsrkjI/KExKSjKLci31U/F25swZunbtioeHBzNmzCjUwxdzSjg3ZUpDJk1y\nK5IPefKTJP8wA4cPG0/AtrKCFSugUSPTxQT68c46nc7wOj09nXfffZdhw4YVmQmpQghh7vSJAxoa\nbg499IIAACAASURBVJpCQ4M4eTKO8ePd6NULCvH93kPLmhxBiKIsJiaGXr168d577zF27NhC1Tv2\nMMqWhV69imbPuylJw6wQu3cPli2Dn36C8+ddefbZYE6dCkKj0a+M7ubma7LY/vzzT3777TcWLVpE\n2bJlUUpRsmRJOnXqxPDhw4tcBSKEEObs8GG4e/e/19bW+vlkps7Sm5uHGe6enp7O0aNH2bVrFzt3\n7mTz5s3SKBNmYdmyZYwcOZI5c+bQvXt3U4fzUO7edSU9PRhLyyDKlNHfp7q6mu4+taiSoYwmlNsX\nz9mzEBAAc+aAqyuMGqWfQ1aiRMHMy8otrmvXrrFo0SLmzZvHjRs3GDJkCEOGDKFWrVoyX6yYMody\nbQ7XIEROdDpYtUq/ZtCJE4ry5Sdw6lRDw8O9kBD/Qllf5zZ3+c6dO+zdu9fQENu9eze2trZ4eHjg\n6elJq1at8Pf3Z8OGDSQmJppFuZb6qXhRSvHVV1/xyy+/sHr1apoUgXUplILvv4evvoKZMxV168r9\n4IPIHLNCKKdkHiNG6Mf7R0XBq6/CW29BvXqmiOu/L0Rvb28GDx7MvHnz2LRpE127dmXo0KG0a9cO\nCwtZbaG4M4dybQ7XIERmCQn6B3s//aSfM/buu/ohRSVLFo2kS1nnLltaWvLMM89w8uRJXF1djRpi\n1apVMzpWkn+Ioio1NZURI0bw119/sWbNGhwdHU0d0gNdvw6vvQZXr8LixfqlmsSDScOsEMopmUet\nWg157z03Xn0VtFpTxvXfF6JGo8HV1ZVRo0bRr18/bG1tTROYKJTMoVybwzWI4iW3UQ3HjsGMGbBo\nkT674pgx4O5uykgfnU6n4+eff2bcuHGG+culS5cmICCAfv36YWVl9VCfYy7l2lyuQ2SXuRw7OjrS\nu3dv7OzsDFNECrsdO2DAAOjfHz7/3DzmqRYUSf5RCF24kD2Zx/Llpk3mkZKSwqpVq0hJSTG8Z21t\nzYIFC3ArrBMRhBCiGMk+2iKEN97w58cfNezdq182JS4OisDDdiOXLl3i119/Zfbs2Tg4ONCkSROO\nHj0K6EduDBkypND28AnxqDKPTlJKodFoGDFiBF9//XWhT5ym0+mHLc6YAXPnQufOpo7IvEjDrIBF\nR+vH4kZGulKzZjAXLpg+mceBAwcIDAwkODiYVq1a0bJlS8NTnMK4gKEQQhRXWbMrLlsWxIEDcXz0\nkRvLlukf8hUVSimioqKYOXMmkZGR9OnTh1WrVtGkSRNZ61KYtawZRMuUKcPgwYMLfaPs0iUYNEjf\nsbBvHzg5mToi8yMNswKQlqZPbf/993Dxon54SWCgBq3280xfPL758sWT25fb7du3Wbx4MYGBgVy9\nepVhw4bx119/4ezsLF+IQghRSF2/Dqmp/722stJn7zX10ik5ye27JCEhgd9//52AgAAsLCwYOXIk\ngYGBRkPlNRqNjNQQZi0tLc3we2Fcnyxr+d20ScPgwfDGGzBpEpSUFkS+kD/rU5T1P3FiooY5c/QZ\nsZyc4L33oGdP+O+BSP5+8eSU2WrcuHEEBgYSGhpKu3bt+Oyzz+jYsaPRUxr5QhRCiMIlPh6++w5C\nQ11xcgrm8uUgLCxMv3RKbnL6/vnggw8ICAggNDQUb29vZs2aRZs2beThnyhWlFKEhYWh0WgoW7Ys\nFhYWhW50UtYh087OISQk+LNggYaXXjJxcGZOkn88JZn/EysF1arFc/WqP97eGsaONc0E7JwSeTg6\nOjJq1Chee+01qlatWvBBCbNiDuXaHK5BmCelYPNm+PZbOHBAv3TKyJFgb1/4RzVk/f6xsLCgSpUq\njB49mtdff50qVark6/nNpVyby3UIvbS0NN566y327t1LWFgY169fBwpfOc6aoK5EiSAiIxvSrp08\ntH8aJPlHAYiLiyM8vCEpKfr/xKdPB7F6dRwvv2y6/8SHDh3iXqYMI1ZW/8/efcdVVb8BHP8clkz3\nSEUl9wDUrEzJlSPE3GbuPVJ/ampaaaYm7lEuCDU1NAI1zYlbUCxXJqCWK3FlpLlABLnc8/vj5hUQ\ncd3F5Xm/Xr66cC/nPCc44/mO5+vEpk2bqFGjhtliEkIIofOkoX4PHkBYGMyZA6mpMHKkbji8o+PD\nn7T8UQ1xcXEZ7j958uRhy5YtOWJdJkPw8PAgb9682NraYm9vz+HDhx/7zLBhwwgPD8fZ2ZkVK1bk\nmv83uVVCQgIdO3YEYN++fbi5uVl0Sfz/iqICugXpCxUyXyy5iSRmL0lVYfdu+PxzuH//0ffz5NGt\nH2Nq9+/fZ/Xq1QQEBHDt2jUqVapEXFwciqLg6+tLdUuciCCEELlMVtUVv/nGnyVLFBYs0K1hOW0a\n+PqCBTWkZys1NZUNGzYQGBjIiRMnKFeuHFeuXNHff3JTo6CiKERERFCwYMEs39+6dSvnzp3j7Nmz\nHDp0iEGDBnHw4EETRylM5erVq7z33nu8+eabLFq0CDsLnqCVlgarV3uhqqE4O4dga6sbMu3lZXlD\npq2R5f5lWLi0NF0L5vTpkJQEo0d7UbJkKDt3hgCm/yM+d+4c33zzDd999x1vvPEG48ePp3nz5tjY\n2Fj8kBchhMhtMldX/OmnELZti6VNG282b4aclMNcuXKFxYsXs3TpUsqXL8+gQYNo164dDg4Oufr+\nk90QxI0bN9KzZ08Aateuze3bt4mPjzf6EE9herGxsbRo0YLBgwfzySefWPR58PffurXJFEXhwoUp\n/POPcQvUicdJYpaNrIaZpKRAcDDMmqXr1h0/Hlq2BBsbhV69jFtlMXM8Wq2WLVu2EBAQwK+//krv\n3r05ePAg5cqVy/Bzlj7kRQghciOt9tFrRYF166BJE/PF8yRZ3Qu1Wi27du0iMDCQyMhIunTpwo4d\nO/D09Mzws7n1/qMoCk2aNMHW1paBAwfSv3//DO9fvXqVUqVK6b92d3fnypUrWSZmEydO1L9u2LAh\nDRs2NFbYwsB27txJ165dmTdvHp07W3aPU0QEdO0KffvChAlga6vwyiu58/w1tIiICCIiIp7ps5KY\nPUHmYSZNmoTx1lu6RTxr1IClS6FevYxDTIxZzTBzhSsPDw/u3LlDiRIlGDx4MOvXr8cpJy1gI4QQ\nudDD4e/Tp+uGCuXJE4K9PTRvfpLGjS3vwS3zvadx48b4+PiwePFiXFxcGDRoECtXrsTV1dXMkVqW\nAwcOULx4ca5fv07Tpk2pXLky9erVy/CZzD1qT2rMTZ+YiZxj+fLlfPrpp6xdu5b69eubO5wn0mph\nxgxdBfHvvoN33zV3RNYnc4PKpEmTnvhZScyeIKthJikpsYSHe5tlvRhdcZGtJCUlAXD69GmCg4Pp\n1KmT6YMRQgjxXNLS4McfYeZM3fD3MWMUtmyZwunTlj1UKPNCuBs2bCAtLY3g4GDeeusti4zZEhT/\nb5J5kSJFaNu2LYcPH86QmJUsWZLLly/rv75y5YpFF4IQ2Uvfq+zp6cnEiRNZtWoVkZGRVK5c2czR\nPdm//0KPHnD7Nhw5Auk6cYWZ2Jg7AEt1/bpuZfOHXFx088lMnZRpNBrWrl1Lz549uZ+uuoijoyNV\nq1Y1bTBCCCGey/37EBgIFSvqWqS/+AJOnIBevSBPHt0oC29vb4tMcO7du8e6dev0DYIArq6uTJky\nhTp16lhkzJYgKSmJhIQEQPf/cMeOHY+tUdWqVSuCg4MBOHjwIPnz55f5ZTnUw15lHx8ffHx8KFu2\nLNu2beOXX36x6KTs0CGoVQuqVNENY5SkzDJIj1kmV6/qunRXrvTCwyOUa9dCUBTTF/O4efMmS5cu\nZdGiRZQuXZqxY8cSFhbGjh07ACxuMUIhhMjNMs/Dun1bISAAFizQrWMZHAw+PmYO8hn98ccffPPN\nN6xcuVL/sHn8+HFA7j3PIj4+nrZt2wK6xtWuXbvSrFkzgoKCABg4cCB+fn5s3bqV8uXL4+LiwvLl\ny80ZsngJmXuV79+/T1hYmMUm2qoK8+fDlCmweDG0aWPuiER6ssD0f65c0fWIhYRA794wejQUK2bc\nRTyzmlD9+++/M3/+fEJDQ2nZsiXDhw+nVq1aT/y8EOZk6ef1s7CGYxDmlX5OsqpCiRInuX7dn9at\nFUaPhmrVzB3h06WmprJx40YCAgI4efIkffv2ZcCAAZQpUyZH3nus5by2luOwZjExMdStW5d79+4B\n4ObmRlRUlMUUvkl//pYu7UW/fgoXLsCaNVC2rJmDy6Vkgen/ZHVzuXxZt1ZMaKiuEs3vv8OjRg7T\nFfOoWbMmTk5OREdHM3DgQE6dOqUfo66PxojFRYQQQryY2NhYwsM9uXdPNyc5Li6EzZtjadbMsq7X\nWd0D//rrL5YsWcLixYspV64cgwcP1pe6f0juPUI8mY2NDVqtljx58uDg4GBRvcrpG420WrCxCaNb\nN39WrVLSLVgvLEmuScwyV1ls0CAMd3d/Vq9W6N8f/vgDihY1XTyZu74PHDjAxIkT2bBhA45ytggh\nRI4QHQ1jxsB/jeUAODrCK6+YL6asZG4MrFWrFoULF2bPnj106tSJbdu2WczDpBA5xf79++nQoQPf\nfPONfgF1S+pVzlzIzskphEGDYnF0lIYWS5VrErPMf5xbtoTQp08sp097U6SIaWP5999/+eabb/RJ\nGYCLiwutW7eWpEwIIXKAgwd1czSOHoURI7xwdg5l9+4QwPRzkp9F5sbAffv2MXbsWJYvX46bm5uZ\noxMi5/nxxx8ZNGgQ33//PU2bNjV3OFlKSclYyM4u1zz151y55lf0778Z/zhdXWH4cEyalF24cIG5\nc+fy/fff06ZNG959911+/vlnQCZUCyGEpVNV2LtXl5CdP6/rKVu9GpycFFR1SrphgpZV+j42NpbJ\nkydnaAx0dXWlY8eOkpQJ8QIWLlzItGnT2L59OzVr1jR3OFm6eBE+/NCLokVDuXPHPIXsxPOz+sTs\n9m2YNQsCA70oUyaUv/82/R/nr7/+yqxZs9i1axf9+/fn5MmTFC9ePEdOqBZCCGuX+doMClu26BKy\nmzfhs8+ga1ewt3/0M5Y2D+vBgwesX7+eRYsWcf78eQYMGEBycjKRkZGANAYK8SJUVWXs2LGsW7eO\nqKgoXn31VXOHlKXt26FnT916iR99NIUTJyyz0Ug8zmqrMiYl6coUz5kDLVvq1o4pXdp0VRY9PT3Z\nsWMHs2bN4syZM4wYMYL+/ftL66SwKtZQMcwajkEYTub5yJ6eJ0lK8kdRFMaOhfbtwdbWzEFm4+rV\nqyxevJglS5ZQqVIlhgwZQuvWrbG3t89VjYHWcl5by3FYgwcPHtCvXz/Onj3Lpk2bKFy4sLlDeoxW\nq2tACgzUFbWrX9/cEYmsZHdeW11ilpoKS5eCvz/UrQuTJ4Mp1vdLP7Fao9Hg4OBA6dKlGTNmDJ06\ndcI+fdOqEFbCGh4arOEYhOHExMTg43NCPx/ZxiaEefM8GTLEG0vJYx7v0YPIyEgWLVrE7t276dy5\nM4MHD6ZaTqjTbyTWcl5by3HkdAkJCXTo0IE8efIQGhqKs7OzuUN6zK1b0L27bqTY6tVQooS5IxJP\nYrXl8tPfnKpV8yI0VOGLL6B8ediwAV5/3XSxHD16lE2bNpHy30Q2RVFYtWoV1atXN10QQgghXlhq\nKqxfn7HCoouLrtXZkpKy9NUVK1SoQEpKCqqqMmTIEL799lvy5s1r5iiFyNnSP18WKVKE9957j1q1\nahEQEICdBVbQOH5c15vfsqVu+o70BeRclvfX9YzSDzdJSwN7+zCqVPHn228VGjY0XRz37t0jKCiI\nadOmodFo9N+3s7Oz6mEiQghhLR48gOBgmDoVPDy8qF8/lGPHLLPCom7NtHD9YraxsbEEBgbSt29f\nuecIYQDpGz8e9mqMHj2aL774wiLPsRUrYPRo3fSdTp3MHY14WTk2MYuNjWXrVk+Skrr8950QgoJi\nqV7dNJOv79y5w8KFC5k/fz4NGjRgx44dTJ06le3btwMysVoIISzdgwe6h5qpU6FCBV1y9vbblllh\nMS0tja1btzJ16lR9Ugbg5OTEm2++aRExCmENMi8t4ejoSNu2bS3iHEvfk1exohcffaSwdy9EREAu\nHrlsVXJkYnbpEnz6Kdy//+h7dnaGH2qS1UTpf//9l3nz5hEQEICfnx8RERFUqVIFgNWrV+eaidVC\nCJETZHUdT0mB5cth2jTdHOSQEN2c5IcsqcLirVu3+PbbbwkICKBw4cL873//o3jx4uzatQuQRkAh\njCEtLU3/2lJqBKQfKaaq4OAQRsOG/hw5oiCjl63HUxOz27dv069fP06ePImiKCxbtoy33npL/35E\nRAStW7embNmyALRv357PP//cKMHevQvTp0NQEAwZ4kWePMZb0DPzOP5GjRpRqVIlli1bRvv27Tl0\n6BDlypXL8DOWdDMXIjfp06cPW7ZsoWjRovqH8PHjx7Nx40YURaFQoUKsWLGCUqVKmTlSYUqZKyy+\n+24YjRr5M2OGQrVqEBYG6W5nFiUmJoaFCxeyZs0aWrRowQ8//EDt2rUB6N69uzQCCmEksbGxaLVa\nnJ2dsbW1tZjGD11Pnqe+MJFGE8KECbHkzSvPndbkqYnZ8OHD8fPzY+3atWg0mgxDKB5q0KABGzdu\nNEqAABoNLFsGEyZAs2YQHQ3u7sYdbpK5K3vTpk24urpy/PhxebgTwsL07t2boUOH0qNHD/33xowZ\nw+TJkwFYsGABkyZNYunSpeYKUZhB5geZdetCiI+PZc0ab/7Lccwqc29eWloaP/30EwsWLODcuXN8\n+OGH/P7777zyyisZfk4aAYUwjvnz5zNr1iyOHDmin19mKY0fqqorUPSQg4PlFCUShpNtYnbnzh32\n79/Pd999p/uwnR358uV77HOGKuWa1ZCT7dth1CgoXBg2b4ZatR593pg3p+vXr+srLAK4urry6aef\nSlImhAWqV68ecXFxGb6Xfs3AxMREi1xzRhhXamrGBxlnZ1i0CCwhp8k8KsPDw4Nbt27h4eHB0KFD\nadeuncUMoRLC2qmqyhdffEFYWBj79+/Hw8PD3CFlkJQE06d74eQUir19CIpieYWJhGFkm5hduHCB\nIkWK0Lt3b6Kjo6lVqxbz5s3LsH6Doij8/PPPVK9enZIlSzJ79myqVq2a5fYmTpyof92wYUMapiuf\nmHnIiY9PGODP+fMKs2ZB69amaRm4fv06M2fOZOnSpXh4eHDt2jUURbGYrmwhzCkiIoKIiAhzh/HM\nxo0bx8qVK3F2dubgwYNP/Fx21yaR82g0sHIlTJrkRf78oSQmhmBjY1kPMroCVltJSkoC4PTp06xc\nuZIPPvjAzJHlTDnt2iQsR1paGkOGDOHo0aNERUVRtGhRc4eUwYUL0LYteHsrXL06hXPnLKswkTAw\nNRtHjhxR7ezs1MOHD6uqqqrDhw9Xx48fn+Ezd+/eVe/du6eqqqpu3bpVrVChQpbbesqu1OjoaNXV\n9XtV11mrqvC9OmZMtJqSku2PGczNmzfVcePGqQULFlQHDx6sXrlyRdVqtWp0dLQaHR2tarVa0wQi\nRA7ytPPalC5cuKB6enpm+d60adPUXr16ZfmeJR2DeDkajap+/72qVqigqg0aqOr+/arFXcc1Go26\nbt069bXXXlMVRVEBFVDd3NzU6Ohoc4dnNazlvLaW47BUycnJavv27dXGjRurd+/eNXc4j9mxQ1WL\nFVPVefNU1QIuX8JAsjuvbbJL2tzd3XF3d+eNN94AoEOHDhw7dizDZ9zc3PQ9aM2bNyc1NZWbN28+\nd4KYlpZxyImrK3TtqhtDa0x3795l8uTJVKhQgfj4eI4dO8aiRYsoWbKkfqikt7e3tEoIkYN16dKF\nI0eOmDsMYSSqCuvWQfXqurV8AgNh7154+20s5jp++/Zt5syZQ/ny5Zk5cyYff/wxbdu2xc3NDTc3\nNxmVIYSJJSQk4Ofnh6IobNmyJcPwd3NTVd1C0T16QGgoDBsm88lyi2yHMr7yyiuUKlWKM2fOULFi\nRXbt2kW1TAslxMfHU7RoURRF4fDhw6iqSsGCBZ8riIMHYfBgL9zcQrGzM96QEzXdHLayZcsSEBDA\n7Nmzeffddzl48CDly5c36P6EEOZz9uxZKlSoAMCGDRuoWbOmmSMSL0vNNA8ZFLZuhfHjde/PnAnN\nm5vnASZzbA+TwNOnTzN//nxCQkLw8/MjLCyMN998E4BOnTpJdUUhzOCff/7Bz8+P119/nUWLFmFr\na2vukPTu3YO+feHcOTh0CEqXNndEwpSeWpVxwYIFdO3alQcPHlCuXDmWLVtGUFAQAAMHDmTt2rUE\nBgZiZ2eHs7MzoaGhz7zz69d165GFh8OsWQqdO0/hxAnjjJ1V0020Tv2va65ly5ZEREQ8cU6cECJn\n6Ny5M5GRkdy4cYNSpUoxadIktm7dyunTp7G1taVcuXIEBgaaO0zxEtRM85Bfey2MlBR/EhMVvvxS\nNwfDXHmNmqmQR/PmzenTpw/z58/n119/ZcCAAZw8eZISJUpk+DmpriiE6cXFxdGsWTM6d+7MxIkT\nLapB5M8/ddeyGjVg/35wcjJ3RMLUlP/GOhp/R4qir96YlqZbi2zCBOjeHSZOxOiL4/3222/UrVuX\n5ORkAFxcXPj555/lpijES0h/XudU1nAMuUFMTAw+Pif0pe8VJYRp0zz5+GNvzN3YrYvNR7+8iqIo\nlC9fnk8//ZQuXbrg6Oho3gBzIWs5r63lOMwpfW+2oij4+fkxZswYhg4dajFxeXl5sXOnQvfuuhEA\nQ4bI0EVrlt15/dQeM0NSVZWDBxWGDAE3N9izB4w9pF5VVbZv387QoUN58OCB/vs2NtlOrxNCCGFB\nzp6F+/cffe3qqhu2aO6kDHTDotIvr+Lk5MSaNWuoXr26GaMSQqTvzdZqtaSlpbFs2TK6dOliAXE9\nGgHw6qthXL/uz5o1CvXrmzU0YWYmTczKlv2clBR/Zs1S6NLF+K0Bx44dY8yYMVy+fJmZM2cSEhLC\n9u3bAWSitRBC5ABxcbrRFeHhXnh6hvLnnyGAZZS+P378OHPnzmXTpk2UKVOGa9euYWNjg6+vr4zG\nEMIC6BaZ36bvzXZycsLT09PMUT2My1M/AuDkyRDCw2OpX1+uG7mdSROzq1erERERS926xv3Du3jx\nIuPGjWP37t1MmDCBvn37Ym9vT5s2bWSitRBC5ADXr8OUKbr1yIYMgXPnFNzcpqS7hptnDR+tVkt4\neDhz587l9OnTDB06lHnz5pE/f365vwhhgTQajf61nZ1JH3uzlX4km4sLvPKK+WIRlsOkf6GOjrrh\nJ4aSeXzu7du3mTp1KsuWLWPo0KEEBgZmKH8qE62FEMKyJSTA3Lkwfz506QKnTkGxYg/fNc01PKsK\ni0lJSaxcuZKvvvoKZ2dnRo0axfvvv49DujVd5P4ihGU5fPgwoKsr8LA32xJGS92+7YVGE0qePCE4\nOFjGCABhGUyamBnyDy9zFayyZcvy119/0b59e06cOEHx4sUNsh8hhBCGlzn5efBAISgIpk6FJk3g\nyBEoW9Y8caW/tzRq1Ijq1asTFBTEW2+9RVBQEPXr15ceMSEs3FdffcW8efOIjo7WF36zhN7sJUvg\n888VNmyYQvHi5h0BICyPSasyarVag/3hZa6CZWtry9q1a2nTpo1Bti+EeDprqBhmDceQ02Se+F61\n6kni4/2pVk1h6lTdQtHmkvneAtCxY0cmT55MxYoVzReYeC7Wcl5by3GYkqqqTJo0iR9++IGdO3dS\n2kIWAktNhZEjYedO2LgR5HKSe1lMVUZDtgb8/vvvJCUl6b92dnamrDmaV4UQQjyXzBPfjx4NYenS\nWHr3Nu9QQFVVOXLkCPfTlX90dXVl3LhxkpQJkQNotVpGjRrF3r172bdvH8UejYM2q3//hY4dIU8e\n3aLR+fKZOyJhqXJczfjr168zYMAAhg0bRo0aNXBzc8PNzc1ixg0LIYTI3unTGUvfu7hArVrmiyct\nLY21a9dSu3ZtZsyYkeHe0rx5c7m3CJEDpKWl0a9fPw4dOsTevXstJik7cQLefFN3jdu0SZIykT3L\nKU/zFKmpqSxatIgpU6bQrVs3Tp8+Tb58+aQKlhBC5BBXr+oWT92yxQtv71DOnTNv6fv79+/z3Xff\nMWfOHAoXLszYsWNp1aoViqLIvUWIHOTBgwd07dqV27dvs2PHDlwNWWnuJWzcCH37wldfQbdu5o5G\n5AQ5IjHbsWMHH330EaVKlWLfvn1UqVJF/55UwRJCCMuWkAAzZ0JAAAwcCGfOKOTNa5rS91lVWLx5\n8yYBAQEsXLiQN998k+XLl+Pj45MhBrm3CJEzJCUl0b59e5ycnNi8eTN58uQxd0ioKkybprvmbdmi\n6zET4llYVGKW+Qb6559/MmrUKGJjY5k7d66+JVMIIYTl02jg229h4kRo2hR++w0ezcM3fun7zBUW\n69evT/ny5Vm5ciVt2rRhz549VK1a1agxCCEMK/2zYunSpWnVqhUeHh4sW7bMbOuUpY+pXDkv+vZV\n+PNP3XyykiXNEpLIoSwmMct8A3V3d+eff/7h448/JjQ0FEdHRzNHKIQQ4knSP5h4enoRHq4werRu\n0dQtW+C110wfk67IyDZ9hcWtW7fSs2dPYmNjKSlPS8LI0tLSeP3113F3d2fTpk0Z3rtx4wbdunXj\n77//RqPR8PHHH9OrVy/zBJqDZH5WtLOzo0uXLixYsAAbG/OUTUhfZVZVwd4+DD8/fyIjFZyczBKS\nyMEsJjHLfAM9d+4c4eHhNGnSxMyRCSGEyE76BxOtFlxcwihUyJ/ZsxX8/MBcAx2io6Mfq7A4cuRI\nScqEScybN4+qVauSkJDw2HsLFy6kZs2aTJs2jRs3blCpUiW6detmth6fnCLzs6KDgwMDBgwwW1L2\nKKZHVWbz5AlhzJhYnJxkOLR4fhZVlVGr1epfOzk5UbRoUTNGI4QQ4lnExsYSHq57MElK6kJC16dK\nPgAAIABJREFUQjW+/z6WFi1Mn5SpqsquXbto1KgR48ePp3r16lJhUZjclStX2Lp1K/369ctyvaLi\nxYtz9+5dAO7evUuhQoUkKXtG6f9/5smTxyKmuGg0j147OJivMUrkfBZzFbh37x4ajQYnJyfs7Oyk\n/L0QQuQA9+/DkiVw796j79nbg6mfMbVaLRs3bmTq1KkkJiby2Wef0alTJ+zs7KTCojC5ESNGMGvW\nLH3ylVn//v155513KFGiBAkJCaxevfqJ25o4caL+dcOGDWnYsKGBo805HB0d0Wq15MmTBwcHB7M/\nK6oqrFvnhY1NKC4uIdjYmK/KrLBcERERREREPNNnFdVES8pnt8r1L7/8QuvWrVmxYgXu7u6A3ECF\nyAmyO69zCms4BnNQVfjxRxg9GmrWVElOHkdUlCegezAJC/M3+DU8qwqLGo2GsLAwpk2bhqOjI2PH\njqVNmzZmHdokzM+c5/XmzZsJDw9n0aJFREREMGfOnMfmmPn7+3Pjxg2+/vprzp8/T9OmTYmOjsbN\nzS3D5+T69MipU6do1qwZEyZMoHbt2oB5nxXv39eVwv/zT1i/XuX6dWkAEs8mu/Pa7D1mBw4coG3b\ntgQHB+Pr62vucIQQQjzF8eMwfDjcvg3LlkGjRgqqatzy95kn/b/77rs0bdqUmTNn4u7uzpw5c2jW\nrJk8EAmz+/nnn9m4cSNbt24lOTmZu3fv0qNHD4KDgzN8Zty4cQCUK1eOV199ldOnT/P666+bK2yL\nFh0dja+vLzNnzqR79+7mDof4eGjTBsqUgb17wclJoXhxmVMmXp5Ze8yioqJo164dq1atolmzZqYI\nQwhhQNbQmmsNx2Aq//wDn38OGzbAl19Cv35ga2uafcfExODj46Of9K8oCnXr1mX69Om8/fbbpglC\n5BiWcl5HRkYye/bsx3rMRo4cSb58+ZgwYQLx8fHUqlWLmJgYChYsmOFzlnIc5nT06FFatGjBggUL\n6Nixo7nDITYWWraEXr1gwgSZTyaen0X2mO3fv5/27dvz/fff07RpU3OFIYQQIpPMQwZTUxUWLIDp\n06FbN/jjDyhQwLQxJSUlkZqaqv/a2dmZgIAAWQhaWLyHvbhBQUEADBw4kLFjx9K7d2+qV6+OVqtl\n5syZjyVl4tFUlyVLltC6dWtzh8OWLbqEbN486NLF3NEIa2SWHrPIyEjef/99fvjhBxo3bmyK3Qsh\njMAaWnOt4RgMKX3pe4AaNU7y99/+VKyoMGcOVK5s2ngSExNZtGgRc+bMIU+ePNy6dQsbGxt8fX0J\nCwuToYsiS9ZyXlvLcbyIyMhIOnToQHBwMM2bNzdrLKqqS8ZmztTNra1Tx6zhiBzOonrM9u7dywcf\nfEBoaCjvvPOOqXcvhBAiG5nX5DlwIIQFC2IZMsS0PVMJCQksWrSIuXPn0rhxYyIiIqhSpYpUWBQi\nF9i5cyddunQhNDTULA346UcNVK7sxbBhCgcOwC+/6OaVCWEsJk3Mdu/eTefOnVm9enWuLvcqhBCW\n6v59SEl59LWrK9SrZ7r9JyQksHDhQr766iuaNGlCREQEVatW1b8vQxeFsG6bN2+mT58+rFu3jnqm\nvPj8J/OoAWfnMGrV8ufAAYW8eU0ejshlTDqU0cHBgR07dtCgQQNT7FIIYWTWMMzGGo7BEHTr8cBH\nH6koyjhu3fJEUUxX+v5hQvb111/TtGlTPv/8c6pUqWLQfYrcw1rOa2s5juykvxacPXuWwYMHs3Hj\nRn1JfFPTFRo6oR81YG8fwqFDntSsKY1CwjAsZiijnZ0dBUw9Y1wIIUS2zpyBoUPh6lVYtUqhfn3T\nlr738PDg77//plmzZkRGRkpCJkQukf5aoNFoUFWVAwcOUKtWLbPGpdU+eu3oaLrqs0KYNDGzlb9s\nIYSwGPfuwdSpEBQEY8fqkjN7ewDFqEMGdfPYtulL3//xxx+sXbvWIqquCSFMJ/O1wMXFBXvdRchs\nTp3yQqMJxckpBDs73agBL6/OZo1J5B4mTcx8fX3x8vIy5S6FECLXyzxsEBR++gk++gjefhtiYqBE\nCdPEkpKSQkhIiP5BDMDJyYlXX33VNAEIISxK+mUwbGxszBaHquqqLi5apHDo0BRsbIw3akCIJzFp\nYialjYUQwrQyT2SvVy8Mrdafy5cVvvsOTFWHSaPREBwczKRJk/D09KRJkyYcOnQIkEY7IXKrw4cP\nY2Njg4uLi34ZDHNcCzQaGDIEDh3SVV4sWVIBZE6ZMD2TJmaSlAkhhGllLn8fHh7CyJGxbNrkjSlG\nDGm1WtasWcMXX3xB8eLF+eGHH6hbt+5jvXhyfxAid1m6dCmTJk3i+PHjJCcnA+a5Fty9Cx07gqLA\n/v3g5mbS3QuRgcnXMRNCCGFaaWmPXru6Qs+eGDwpe3y4JGzZsoXPP/8cBwcHFi5cSJMmTfQPXYpi\n3HlsQgjLtWTJEr788kv27NlDhQoVzBbHlSvQooVuweiFC8FOnoqFmcmfoBBCWKn4eJg2zQtFMe5E\n9sxVFmvVqkVycjKJiYn4+/vTunVr6RETQgCwePFi/P392bt3L+XLlzdbHMePQ8uWMGwYfPyxrsdM\nCHMz6Tpm1r4WhxC5jTWc19ZwDJlptfDttzBuHPTpA+PHq5w/b7xhg7p1f3z0BT0URWHKlCmMGTNG\nqvEKs7CW89pajuOhoKAgpkyZwp49e8yalG3bBj16wKJF8P77ZgtD5FIWs46ZEEII4zp5EgYO1A1f\n3LULdKMFjTts8PLly/o5IgCurq60aNFCkjIhhN4333zD1KlT2bt3L+XKlTPZfjMPs16yRGHCBPjp\nJ6hb12RhCPFMJDETQggrcP8++PvD4sXw5Ze65MzYladv3LiBv78/K1eupFKlSly8eBFFUaTKohAi\ng8DAQKZPn26WpCx9VdqSJcNIS/Nn/34FM3bYCfFEkpgJIUQOk7kFePduhQ8/hNdeg+ho469Jdv/+\nfebNm8fs2bPp1KkTv//+O0WKFJEqi0KIxwQEBDBjxgz27t1L2bJlTbrvzFVpz50LYc+eWMqXl8JD\nwjJJYiaEEDlI5hbgggXDAH8CAhRatDDuvtPS0li1ahXjx4/njTfe4Oeff6ZixYr696XKohAifcPR\nvn37mD17tlmSsqw4O0P+/OaOQognk8RMCCFykMwtwCkpIezfH0vt2oZLirJaY2zHjh2MGTMGZ2dn\nQkNDqSuTM4QQmaSv0JqamoqtrS2xsbFmS8ry5/fCxiYUe/sQHB2NU5VWCEOSxEwIIXKQa9d088ke\ncnQEJyfDbT9z6fs6deqgqioXL15k+vTptG3bVoYpCiGypGs42qav0Orq6qp/bWonToCfn8L48VNo\n1uxhQ1NnuX4Ji/bUqeG3b9+mQ4cOVKlShapVq3Lw4MHHPjNs2DAqVKhA9erV+e2334wSqBBCZKdP\nnz4UK1YsQ9GJ0aNHU6VKFapXr067du24c+eOGSN8OVotBAZC165eVK58Aje3ENzcQv5rATZcoY30\nD1aJiYns3LmTN954g5MnT9KuXTt5qBFCZCs1NVX/2lzXi8hIaNwYZsyAjz/WVaX19vaW65eweE9N\nzIYPH46fnx+///47MTExVKlSJcP7W7du5dy5c5w9e5bFixczaNAgowUrhBBP0rt3b30vz0PNmjXj\n5MmTREdHU7FiRaZNm2am6F7OuXPwzjvw3Xewb59CbOwUoqI8iYryJCzM36APG8nJyTx48ED/taur\nK506dcLe3t5g+xBCWKcjR45gY2ODi4sLbm5uZqnQumaNbm2ykBDoLKMWRQ6T7VDGO3fusH//fr77\n7jvdh+3syJcvX4bPbNy4kZ49ewJQu3Ztbt++TXx8PMWKFTNSyEII8bh69eoRFxeX4XtNmzbVv65d\nuzY//vijiaN6OWlpMH8+TJkCY8fC8OGgWxrM8OuSqarKunXrGDVqFEWKFOH27dvY2NhI6XshxDNZ\ntWoVX3zxBb/99hspKSmA6Su0Llig6yXbsQNq1DDZboUwmGwTswsXLlCkSBF69+5NdHQ0tWrVYt68\neTg7O+s/c/XqVUqVKqX/2t3dnStXrkhiJoSwKMuWLaNzDmo+/f136NMH8uSBgwcx6po7MTExDB8+\nnH///Zdly5bRqFEjKX0vhHhmq1evZsyYMezatYtKlSqZfP9aLXz2GWzYAFFR4OFh8hCEMIhsEzON\nRsOxY8dYuHAhb7zxBh999BHTp0/nyy+/zPA5VVUzfP2km/jEiRP1rxs2bEjDhg1fLGohhFlEREQQ\nERFh7jCe25QpU3BwcKBLly5Zvm/ua1P6KoiVK3sxe7bC3LkwebJxF4q+ceMG48ePZ926dUycOJH+\n/ftjZ6e7LUjpe5GT5NRrkzVYv349w4YNY8eOHVStWtXk+3/wAPr2hfPn4cABKFTI5CEIYTCKmjmr\nSufvv/+mTp06XLhwAYCoqCimT5/O5s2b9Z/58MMPadiwIZ06dQKgcuXKREZGPtZjpijKYwmcECJn\ns7TzOi4ujpYtW+qTHIAVK1awZMkSdu/ejaOj42M/Y+5jSL8umVYLDg4neeMNf5YsUShTxnD7SN8D\nptFoCAgIwN/fny5dujBhwgQKFixomJ0JYQHMfV4biqUfx+bNm+nbty/h4eG89tprJtln+uuZh4cX\nHTooODnBDz/o1ikTwtJld15n22P2yiuvUKpUKc6cOUPFihXZtWsX1apVy/CZVq1asXDhQjp16sTB\ngwfJnz+/DGMUQliEbdu2MWvWLCIjI7NMyixB5nXJtNoQZs2KpUwZw/RYZS5/X7NmTa5fv06pUqWI\njIw0Swu3ECLn27FjB3369GHz5s0mTcoeNmSpKtjbh9Gxoz+LFinYyQJQwgo89c94wYIFdO3alQcP\nHlCuXDmWLVtGUFAQAAMHDsTPz4+tW7dSvnx5XFxcWL58udGDFkKIzDp37kxkZCQ3btygVKlSTJo0\niWnTpvHgwQN9EZA6deoQEBBg5kgzunABkpIefW1vD4ac0pV5XaGoqCjmzZvH//73P5k7JoR4IXv3\n7qVbt26sX7+eN99802T7zdyQlSdPCIMHx2JnJ0OvhXXIdiijQXdk4d3xQojnZw3ntbmOQavVVRD7\n8ksVD49xnD3rCYCv70mDlsA/cuQIb7/9tr4EvpubG1FRUTKHTFg1a7g2gWUeR1RUFO3atWPNmjU0\naNDApPuOiYmhTp0TJCXpEjM3txCiojzleiZylBceyiiEEMLwLl2CXr0gORkOHVIoV25KujlgnQ2W\nlO3cuZPBgwdTuHBh7t69i6IoUv5eCPHCDh06RLt27QgJCTF5UgZw44YXGk0oTk4h2NnpGrK8vHJO\ntV0hnkZ6zIQQL8wazmtTHoOqQnAwfPwxjBoFo0c/XJfMsP766y9GjhzJoUOHWLBgAS1atJDy9yJX\nsYZrE1jGcTwstnHq1CmGDRvGihUr8PPzM3kc69frqtSGhakUKiTXM5FzZXdeS2ImhHhh1nBem+oY\n/vlH91Bx/jysXAnVqxt+HxqNhkWLFjF58mQGDhzIuHHjMqw7KURuYQ3XJjD/cTwsHhQeHs69e/fw\n8fFh//79Jk+Gli2Dzz+HzZvBRHVGhDAaGcoohBAmlLk8/caNCh9+CD17QmiobtFoQ25fURQOHTrE\noEGDyJ8/P/v376dKlSovexhCiFwuNjZWn5SBbo5XbGysSed0zZwJgYEQEQEVK5pst0KYhSRmQghh\nQOnLOQMULhyGovizZo3C228bavuPyt83btyYYsWKsWnTJmbNmkWXLl1kaI8QFiAtLY3XX38dd3d3\nNm3a9Nj7ERERjBgxgtTUVAoXLmyRC2RfvXqVpPRlY01IVeGTT2DLFoiKgpIlzRKGECYliZkQQhhQ\n5nLOKSkh7NsXy1tvGaaFOXP5+w0bNvDBBx9w6tQp8ufPb5B9CCFe3rx586hatSoJCQmPvXf79m2G\nDBnC9u3bcXd358aNG2aIMHvXrl1j6NCh1KhRg3PnzgGYrHiQRqMb+n3yJOzbB4UKGX2XQlgEG3MH\nIIQQ1kSjgf8q0wPg6AiGnuaVfmy6i4sLY8eOlaRMCAty5coVtm7dSr9+/bKcSxISEkL79u1xd3cH\noHDhwqYOMVs3b96kWbNm9OnTh19//ZWoqCiioqIICwszeo98cjK8/z5cuQK7dklSJnIXScyEEMJA\nLl6EoUO9yJ//BK6uIbi5hfxXztkwLcxarZaIiAhSU1PJkycPbm5u+Pn5Sfl7ISzMiBEjmDVrFjY2\nWT9mnT17lps3b9KoUSNef/11Vq5caeIInywhIYHmzZvTvHlzPvvsMxRFwdvbG29vb6MnZXfvgp+f\nbh7upk3g6mrU3QlhcWQooxBCGMDq1fC//8GYMQojRkzh5EnDrkt25swZ+vTpA+gm4KekpPy3fSkX\nLYQl2bx5M0WLFqVmzZpPnDeWmprKsWPH2L17N0lJSdSpU4e33nqLChUqPPbZiRMn6l83bNiQhg0b\nGidw4P79+7Rq1YqaNWsyY8YMk1xbHhYzunkTRo70onZthYULjbOUiBDmEBER8cxzSKVcvhDihVnD\nef2yx3DvHgwbppsH8cMP8PrrBgwOXQn8r776ihkzZjBhwgSGDBnyxFZ4IYSOOa9NY8eOZeXKldjZ\n2ZGcnMzdu3dp3749wcHB+s/MmDGD+/fv65Oufv364evrS4cOHTJsy5THkZqaSrt27cibNy/BwcHY\nmiAzelgsKTzck6QkqFz5JCdO+GNjI41Nwnpld17L3V0IIV7QsWO6NXXS0nSvDZ2UnTx5krp167Jt\n2zYOHz7M0KFDJSkTwsJNnTqVy5cvc+HCBUJDQ3nnnXcyJGUArVu3JioqirS0NJKSkjh06BBVq1Y1\nU8S6CpI9evRAURRWrFhhkqQMHpbj9+TevS6oaheuXKnGiROxJtm3EJZIhjIKIcQzSL92WLVqXsyb\npzB9OsybB507G27bXl5eaDQaZsyYwbx585gyZQr9+/eX4YpC5FAPz92goCAABg4cSOXKlfH19cXb\n2xsbGxv69+9vtsRMVVUGDRpEfHw8W7duxd7e3mT7Pn0azFSNXwiLJEMZhRAvzBrO62c5hsxrk7m4\nnMTDw58fflB49dWX23/mdcnq1KnDP//8Q/HixQkKCqJ06dIvtwMhciFruDaB8Y9DVVXGjBnD/v37\n2blzJ25ubkbbV2aHDkHLlirlyo3j5EndtdXX9yRhYf7SECWsWnbntSRmQogXZg3n9bMcQ0xMDD4+\nJ/Rrkzk4hPDLL5689trLr02m27aPfl0ygMmTJzNu3Dh5OBHiBVnDtQmMcxzpe+g3bNjAmjVriIiI\noGDBggbdT3YiIqBjR1i+HPz8Mo4YkOuesHbZndcylFEIIZ4i89pkefKAnQGvnukv0K6urrRq1Uoe\nToQQBpe+hz41NRV7e3vOnDlj0qRs61bo2VNXybZRIwBdOX4hhBT/EEKIbP39t66Ec758hl+bTKPR\nsHnzZlJTU3F0dMTNzY3mzZvLumRCCKOIjY1l27ZtJCYm6pfcuH79usn2/+OP0Ls3bNz4MCkTQqQn\nPWZCCPEEERHQtSv066ewa9cUTp0y3NpkZ8+epUePHjg7O3PmzBnu3Lnz37ZlKI8QwnjS0tL0r015\nrQkOhk8+ge3boUYNk+1WiBxFEjMhhMhEq4UZM2D+fPjuO2jWDAw13Ear1RIYGMiECRNkXTIhhEnd\nuXOHtLQ0XFxcsLGxwdfX1yQ99IGBMHUq7N0LlSsbfXdC5FiSmAkhRDr//gs9esDt23DkCLi7G27b\nV65coU+fPty5c4cDBw5QqVIlw21cCCGyERsbS4cOHdi4cSPFixcHTNNDP2uWLjGLjISyZY26KyFy\nPGmmFUKI/xw+DLVqQZUqumGML5qUqapKTEwMMTExqKqKqqqsXLmS1157jfr160tSJoQwqYsXL+Ln\n58fXX3/Nu+++i7e3N97e3kZJyh5e/6KjYxg/XuXbb2HfPknKhHgW0mMmhMj1tFqVRYsUJk+GxYuh\nTZsX31bmdcneeecdbG1tOXPmDNu3b6dmzZoGiloIIZ7u+vXrNGvWjNGjR9O5c2ej7iv9mo8pKeDs\nHMbp0/4UKybzZoV4FpKYCSFyvTJlPqdwYX9++UWhXLmX21b6qmcAmzZtokePHhw9ehRHR0cDRCuE\nEM8mMTGRFi1a0KFDB4YNG2b0/emuf576NR/T0kKIj4+lWDEphy/Es5ChjEKIXC8+vhqLF8e+dFKW\nFScnJ0aOHClJmRDCpB48eEC7du2oXr06/v7+JtlnWhokJz/6WgrMCvF8pMdMCJHrOTrqFo02BHt7\ne2xtbbGzs8PR0VHWJRNCmJxWq6VXr164uLgQGBhokrL4Gg3Mnu1FgQKh3L8fgqLw35qPxh0+KYQ1\nkcRMCJHrGerhITg4mFGjRjFjxgzeeustFEWRdcmEECalqiojRozgypUrbN++HTs74z/qaTTQvTv8\n+6/ChQtTOH/ecGs+CpGbKKqqqibZkaJgol0JIUzEGs5rRVHQarUv9fCQmJjIkCFDOHz4MKtXr5Ye\nMiHMzBquTfBixzF16lTCwsKIjIwkf/78RorskdRU6NoVEhJg3TpwcjL6LoXI0bI7r6XHTAiR671M\nUhYTE0PHjh2pU6cOR48excXFxYCRCSFE9lRVJTZW10N16NAhli5dyoEDB0ySlD14AJ066f67fr1u\nWLgQ4sVJYiaEEC9AVVWCgoIYP348c+fOpXv37uYOSQiRy6RfniMtLQ1VVYmOjtYvIG1MKSnQsaPu\n9Y8/Gm6erhC5mSRmQgjxDNK3SpcqVYoBAwZw9uxZoqKiZLFoIYRZZF6ew8XFheT0ZRGNJDkZOnQA\nBwcIDdX9Vwjx8qRcvhBCPMXDVmkfHx/q1KlDqVKlKFKkCAcPHpSkTAhhVunnqtjYGP+xLjkZ2rbV\nzSULC5OkTAhDkh4zIYR4isyt0k5OTnz44YeyNpkQwqwKFy6Mqqo4Ojpib2+Pr6+vUYsP3b8PrVtD\noUKwciWYoOCjELmKnFJCCPEUWq2WlJQU/demKD8thBDZSUhI4L333uPTTz+ldevWAEZZnuPhMO77\n9+Gzz7woUUJhxQpJyoQwBimXL4R4YdZwXj/tGFJSUujevTt79+4lOTkZRVHw9fUlLCxM1ucRwkJZ\nw7UJnnwcqamptGrVCnd3dxYvXmy0a5FuGPc4tm3z5P59cHc/yblz/tjZybVPiBcl5fKFEOIF3L59\nm7Zt21KoUCEuXbrE2bNnAeO0SgshxLNQVZUPP/wQRVEIDAw06rVIN4zbk8TELgDcvBnCqVOxeHt7\nG22fQuRmkpgJIUQWrl69SvPmzWnQoAFff/01tra28jAihDC7yZMnc/z4cSIjI40+rPr+fd0/IYRp\nSFVGIYTI5NSpU9StW5euXbsyf/58bG1tzR2SEEKwYsUKli9fzpYtW3B1dTXqvh7OKStZ8gRubiG4\nuYXg63vSqMVFhMjtZI6ZEOKFWcN5nfkYoqKiaN++PbNnz5ZFo4XIoazh2gQZj2PHjh10796dyMhI\nKleubNT9Jic/qr4YHKxy6pRuDUcZxi3Ey8vu+iSJmRDihVnDeZ3+GNavX8+AAQNYtWoV7777rpkj\nE0K8KGu4NsGj4zh+/DjNmjVj3bp1vP3220bdZ0oKtGsHLi4QEiLVF4UwtOyuT880lNHDwwNvb29q\n1qzJm2+++dj7ERER5MuXj5o1a1KzZk38/f1fLmIhhHhOffr0oVixYhmG2axZs4Zq1apha2vLsWPH\nnvizqqoSEBDAkCFD2LZtmyRlQgiLcenSJd577z0CAgKMnpQ9eAAdO0KePPD995KUCWFqz3TKKYpC\nREQEBQsWfOJnGjRowMaNGw0WmBBCPI/evXszdOhQevToof+el5cX69evZ+DAgdn+bNWqVUlNTWX/\n/v2UK1fO2KEKIcQz8/X1ZdSoUXTo0MGo+0lNhS5dQFUhNBTs7Y26OyFEFp65LeRpQwKsYciAECLn\nqlevHnFxcRm+96zzMM6cOcPu3bslKRNCWJzU1FQ++ugjo+5Do4Hu3SEpCdavBwcHo+5OCPEEz9xj\n1qRJE2xtbRk4cCD9+/d/7P2ff/6Z6tWrU7JkSWbPnk3VqlUf287EiRP1rxs2bEjDhg1fKnghhGlF\nREQQERFh7jAMzs7OjrCwMCIiIuTaJEQOZK3XJoD4+HhiY423dlhaGvTqBTdvwsaNumGMQgjzeKbi\nH9euXaN48eJcv36dpk2bsmDBAurVq6d/PyEhAVtbW5ydnQkPD2f48OGcOXMm446sZCKuEOIRSzuv\n4+LiaNmyJbGxsRm+36hRI+bMmcNrr7322M8oisL7779PWFiYVBsTwkpY2rXpRSmKgpubG1FRUUZJ\nzLRa6NsXLl6EzZvB2dnguxBCZPLSxT+KFy8OQJEiRWjbti2HDx/O8L6bmxvO/53NzZs3JzU1lZs3\nb75MzCKHKliwIIqiyD8r+5fd/FJrIEmZEMJS+fr6GnTtMFVViYmJ4fjxGAYMUPnzT9i0SZIyISzB\nUxOzpKQkEhISALh37x47dux47AIRHx+vz/wOHz6MqqpW/yAnsnbr1i1UVZV/Vvbv1q1b5v7TemnZ\ntZ5LUiaEMLS0tDRq1qxJy5Ytn/iZI0eOYGdnx7p16574GUM2HKmqSseO4/DxOcGbb57gp58+Z9Mm\nFRcXg2xeCPGSnjrHLD4+nrZt2wKg0Wjo2rUrzZo1IygoCICBAweydu1aAgMDsbOzw9nZmdDQUONG\nLYQQmXTu3JnIyEhu3LhBqVKlmDRpEgULFmTo0KHcuHGDFi1aULNmTcLDw80dqhAiF5g3bx5Vq1bV\nN25nlpaWxieffIKvr6/JGo5iY2PZts2TxMQuAKSkhBAXZ7z5a0KI5yMLTAuDkt+zdXrS79Uaft/W\ncAxCiIzMfV5fuXKFXr16MW7cOObOncumTZse+8zXX3+Ng4MDR44c4b333qN9+/aPfcbDfwsBAAAg\nAElEQVTQxxETE8Mbb5zgwQNdYubmFkJUlKckZkKYUHbn9TPNMRNCCCGEEM9mxIgRzJo1CxubrB+z\nrl69yoYNGxg0aBBguuHU4eFeODicwNU1BDe3EHx9Txp0/poQ4uXImu5CWJFLly5RrVo17t69K/Om\nhBDCDDZv3kzRokWpWbPmE0v4f/TRR0yfPl3fcp5dr5ihlhr65hsIClL4/fcp3Lypq1zr5dVZ7hVC\nGNnzLOchQxmFQcnv+cXFxcVRtmxZNBrNE1tZM/Pw8GDZsmW88847Ro1NhjIKIXISc57XY8eOZeXK\nldjZ2ZGcnMzdu3dp3749wcHB+s+ULVtWH9+NGzdwdnZmyZIltGrVKsO2DHUcq1bBp5/Cvn1QtuxL\nb04I8RKyO68lMRMGJb/nF/cwMUtNTcXW1vaZfubVV19l6dKlNG7c2KixSWImhMhJLOW8joyMZPbs\n2VnOMXuod+/etGzZknbt2j32niGO46ef4MMPYfduqFbtpTYlhDAAmWMmLIKq6tZOiYmJeeEbzctu\n4/Lly7Rr146iRYtSuHBhhg4diqqq+Pv74+HhQbFixejZsyd3794FdMmSjY0NK1asoHTp0hQqVIhv\nvvmGI0eO4O3tTYECBRg6dKh++ytWrMDHx4ehQ4eSP39+qlSpwp49e/Tve3h4sHv3bv3XEydOpHv3\n7gDUr18fgPz58+Pm5sahQ4c4f/4877zzDoULF6ZIkSJ069aNO3fuANC9e3cuXbpEy5YtcXNzY/bs\n2fp4tVotAH/99RetWrWiUKFCVKhQgaVLl2bYd8eOHenZsyd58+bF09OTX3/99bn/nwohhHiyh0MF\ng4KC9BWtTWXnThgwALZskaRMiBxBNRET7kqY0ZN+z1qtVu3Q4TPV1fV71dX1e/X998eqWq32ubb9\nstvQaDSqt7e3OnLkSDUpKUlNSUlRo6Ki1G+//VYtX768euHCBTUxMVFt166d2r17d1VVVfXChQuq\noijqoEGD1JSUFHXHjh2qg4OD2qZNG/X69evq1atX1aJFi6qRkZGqqqrq8uXLVTs7O/Xrr79WNRqN\nGhYWpubLl0+9deuWqqqq6uHhoe7evVsf08SJE9Vu3bqpqqqqcXFxqqIoalpamv79c+fOqbt27VIf\nPHigXr9+Xa1fv7760Ucf6d/PvL2H8T7cRr169dQhQ4aoKSkp6vHjx9UiRYqoe/bsUVVVVSdMmKA6\nOjqq4eHhqlarVT/77DP1rbfeyvL/3ZN+r9ZwXlvDMQghMrKW8/pljuPAAVUtUkRV9+0zYEBCiJeW\n3XktPWbCJNKvnZKY2IVt26oRGxtr0m0cPnyYa9euMWvWLJycnHBwcMDHx4fvv/+eUaNG4eHhgYuL\nC9OmTSM0NFTf6wQwfvx4HBwcaNq0KW5ubnTp0oXChQtTokQJ6tWrx2+//ab/bNGiRRk+fDi2trZ0\n7NiRSpUqsWXLlixjUtP1+qlZ9ACWK1eOxo0bY29vT+HChRkxYgSRkZHPdLyXL1/m559/ZsaMGTg4\nOFC9enX69euXYZ5DvXr18PX1RVEUunXrRnR09DNtWwghhOX67Tdo2xZWroR69cwdjRDiWUlVRmEW\nCQlQvbpp93n58mXKlCnzWGGNa9euUaZMGf3XpUuXRqPREB8fr/9esWLF9K+dnJwe+/revXv6r0uW\nLJlh+2XKlOGvv/56oZjj4+MZPnw4UVFRJCQkoNVqKViw4DP97F9//UXBggVxcXHRf6906dIcPXpU\n/3X643B2diY5ORmtVvvMxUeEEEJYlt9/Bz8/CAiAd981dzRCiOchT1/CJLy8vPD1PYGbm27tlPff\nP4lW64Wq8sz/tFovOnR4tI3nXX+lVKlSXLp0ibS0tAzfL1GiBHFxcfqvL126hJ2dXYak5XlcvXo1\nw9cXL16kRIkSALi4uGRI4v7++2/966xKFo8dOxZbW1tOnDjBnTt3WLlyZYaevOzKHJcoUYKbN2+S\nmJio/96lS5dwd3d//oMSQghhkdR0c6///FOlWTOYPh2yWK9aCGHhJDETJqEoCqtXTyEqypOoKE/C\nwvyfe+2Ul91G7dq1KV68OJ9++ilJSUkkJydz4MABOnfuzFdffUVcXByJiYmMHTuWTp06PVevUfph\niP/88w/z588nNTWVNWvW8Mcff+Dn5wdAjRo1CA0NRaPRcPToUX788Uf9MRQpUgQbGxvOnz+v31Zi\nYiIuLi7kzZuXq1evMmvWrAz7LVasWIbPp1eqVCnq1q3LZ599RkpKCjExMSxbtoxu3bo983EJIYSw\nXKqq0rHjOHx8TlC37gm8vT9nzBiVnj3NHZkQ4kVIYiZMRlEUvL298fb2fuEFLV9mGzY2NmzatIlz\n585RunRpSpUqxZo1a+jTpw/du3enfv36lC1bFmdnZxYsWJBhn88S10O1a9fm7NmzFClShPHjx/Pj\njz9SoEABACZPnsz58+cpUKAAEydOpGvXrvqfc3Z2Zty4cfj4+FCwYEEOHz7MhAkTOHbsGPny5aNl\ny5a0b98+w74+++wz/P39KVCgAHPnzn0slh9++IG4uDhKlChBu3bt+PLLL/VrnimK8tixyUKjQgiR\nc6Sfe33vXhc0mmo0aPB887eFEJZD1jETBpXbf88rVqzg22+/Zf/+/eYOxaBkHTMhRE5iLef1044j\nJiYGH58TJCZ2AcDNLYSoKE+8vb1NFaIQ4jlld15L8Q8hhBBCiByocmUvXFxCSUkJwdGR/+ZedzZ3\nWEKIFySJmRAGlNXwQCGEEMLQtFro21fh9den8OWXsdjZgZdXZ7kHCZGDyVBGYVDye7ZOMpRRCJGT\nWMt5nd1xfPwx/PIL7NwJzs4mDkwI8cJkKKMQQgghhJWYNQu2bYN9+yQpE8KaSGImhBBCCJFDfPcd\nLFwIBw5AwYLmjkYIYUiSmAkhhBBC5ABbtsAnn8DeveDubu5ohBCGJomZEEIIIYSF++UX6NULNm+G\nKlXMHY0QwhhkgWkhhBBCCAt26hS0aQPBwVC7trmjEUIYiyRmQmQhPj6e+vXrkzdvXkaPHm3ucIQQ\nQuRSly9D8+Ywe7buv0II6yVDGYXIwuLFiylatCh37941dyhCCCFysWbNVIYNU+je3dyRCCGMTXrM\nRK6k0Wiyff/ixYtUecFB/E/bthBCCPGstNrPGTky56/JJoR4OknMhMmoqkpMTAwxMTEvvPDny2zD\nw8ODmTNn4u3tjZubGwcOHKBu3boUKFCAGjVqEBkZCUCvXr0IDg5m5syZuLm5sWfPHlRVZfr06ZQv\nX57ChQvzwQcfcOvWLQDi4uKwsbFh2bJllClThiZNmgCwbNkyqlatSsGCBfH19eXSpUv6WGxsbAgK\nCqJixYoUKFCA//3vfxliXbJkCVWrViVv3rxUq1aN3377DYC//vqL9u3bU7RoUcqWLcuCBQte6P+j\nEEKInOHatWrExsaaOwwhhAlIYiZMQlVVOnbsiI+PDz4+PnzwwQfPnVgZYhuhoaGEh4dz/vx5Wrdu\nzRdffMGtW7eYPXs27du3599//2XFihV07dqVTz75hISEBN555x3mz5/Pxo0b2bdvH9euXaNAgQIM\nGTIkw7b37ft/e3ceFOV9x3H8jWIGNXiEcKiAICiHcqwgqAkZE0SpQYzGeiQaxWLqMdHEaDyqrUYl\ntKZptKk1iremIiYqKoNMMMRqROQmUgcxEi4lKloFtAnw9A+HJ6yAtzy7+H3NPDPuPr/9PZ9nl/26\nv+c8ypkzZ4iPj2f//v189NFH7N27l8uXLxMQEMD48eP12h86dIjU1FSys7PZvXs3hw8fBiAmJoZl\ny5axfft2rl+/TmxsLBYWFtTW1jJ8+HB0Oh2lpaUkJiby6aefkpCQ8EDvgRBCCCGEMDwyMBPNIicn\nh/j4eCoqKqioqCA+Pv6BtwA+ah8mJibMmjWLbt26sX37doYNG0ZwcDAAgwcPxtfXl0OHDqnt6w/6\nPv/8c1asWEHXrl1p06YNf/rTn9izZw+1tbVqm6VLl9K2bVvMzMxYt24dCxcuxMXFhVatWrFw4UIy\nMzMpKipS2y9YsIAOHTpgZ2fHyy+/TFZWFgBRUVHMnz8fHx8fAJycnLC3t+fUqVNcvnyZxYsXY2pq\niqOjI+Hh4ezateuB3kchhBDGIzj4NB4eHlrHEEI0A7n4h9DEjRs38PLyavbl2tnZAbfPIYuJieHA\ngQPqvOrqal555ZVGX1dQUMDIkSNp1erXbRmmpqaUlZU16Luu/9mzZ/P+++/r9VNSUqK2s7GxUZ9v\n164dFRUVABQXF+Pk5NQgw48//khpaSmdO3dWn6upqeGll16694oLIYQwStHRKzAxMdE6hhCiGcjA\nTDQLDw8PgoOD1cP1goODiY6OfqD/bOoOZazfx4NuRaxbnr29PRMnTmT9+vX39Tp7e3s2b97MgAED\nGswrKCjQ67uu/ZIlSxocvng/7OzsyM/PbzSDo6MjeXl5D9ynEEII4ySDMiGeHnIoo2gWJiYm7N69\nm2PHjnHs2LEHHpQ9rj7qTJgwgQMHDpCQkEBNTQ23bt0iKSmJkpISgAbnrk2bNo1FixapF/C4dOkS\nsbGxTfY/bdo0IiIiyM3NBeC///0vMTExTbZXFEVdZnh4OB9//DHp6ekoikJ+fj6FhYX4+flhbm7O\nX/7yF27evElNTQ3ff/89qampD/UeCCGEEEIIwyEDM9FsTExM8PT0xNPT86EHVI+jDwBbW1v2799P\nREQEVlZW2Nvb89e//lUdHJmYmOj1P3v2bEJDQxkyZAgdOnRgwIABpKSk6OWq77XXXmP+/PmMGzeO\njh074uHhoe7pa6x9/eWNHj2aP/zhD7zxxht06NCBUaNGcfXqVVq1asXBgwfJzMykR48eWFpa8vbb\nb8u91oQQQgghWgAT5WGvW/6gCzIxeehLpAvjIZ9zy9TU59oSPu+WsA5CCH0t5XvdUtZDCPGru32v\nZY+ZEEIIIYQQQmhMBmZCCCGEEEIIoTEZmAkhhBBCCCGExmRgJoQQQgghhBAak/uYiceqc+fOcs+V\nFqj+Ta2FEEIIIcTjJ1dlFEI8tJbwvW4J6yCE0NdSvtctZT2EEL96pKsyOjg44OnpiU6nw8/Pr9E2\ns2bNomfPnnh5eZGRkfFoaR+zpKQkrSM8EsmvHWPODsaf/0FMmTIFa2trPDw81OfKy8sJCgqiV69e\nDBkyhGvXrmmYsCFj/3wkv7Ykv+GrqalBp9MxfPjwBvN27tyJl5cXnp6evPDCC2RnZ2uQsGnG/vkY\nc35jzg6S/1Hdc2BmYmJCUlISGRkZejfUrRMXF0d+fj5nz55l/fr1TJ8+/YkEfVhav8GPSvJrx5iz\ng/HnfxBhYWHEx8frPRcZGUlQUBB5eXkEBgYSGRmpUbrGGfvnI/m1JfkN3+rVq3F3d2/08P4ePXpw\n9OhRsrOzWbJkCW+//bYGCZtm7J+PMec35uwg+R/VfV3842670WNjY5k0aRIA/v7+XLt2jbKysseT\nTggh7kNAQECD8+Dq16ZJkyaxb98+LaIJIZ5CxcXFxMXFER4e3uhvqAEDBtCxY0fg9m+n4uLi5o4o\nhDBA97XHbPDgwfj6+rJhw4YG80tKSrCzs1Mf29raSoERQmiurKwMa2trAKytrWWDkRCi2bz33nus\nWrWKVq3uvf1748aNDBs2rBlSCSEMnnIPpaWliqIoyk8//aR4eXkpR48e1ZsfEhKiHDt2TH0cGBio\npKWlNegHkEkmmVrgZCjOnz+v9OnTR33cqVMnvfmdO3du9HVav38yySTTk5m0cuDAAWXGjBmKoijK\nN998o4SEhDTZ9siRI4qbm5tSXl7e6Hyt30OZZJLpyUxNuefl8rt06QKApaUlI0eOJCUlhYCAAHV+\nt27dKCoqUh8XFxfTrVu3Bv0oclUhIUQzsra25uLFi9jY2HDhwgWsrKwabSe1SQjxOH333XfExsYS\nFxfHrVu3uH79Om+99Rbbtm3Ta5ednc3UqVOJj49v8pYkUp+EeLrcdR97VVUVN27cAKCyspKEhAS9\nq54BhIaGqsUmOTmZTp06qYcPCSGEVkJDQ9m6dSsAW7du5bXXXtM4kRDiaRAREUFRURHnz59n165d\nvPLKKw0GZYWFhYwaNYodO3bg7OysUVIhhKG56x6zsrIyRo4cCUB1dTVvvvkmQ4YM4fPPPwfg97//\nPcOGDSMuLg5nZ2fat2/P5s2bn3xqIYSoZ/z48Xz77bdcvnwZOzs7PvzwQxYsWMCYMWPYuHEjDg4O\n7N69W+uYQoinUN1VGev/dvrwww+5evWqeiXrNm3aNHrlayHEU+YRDqM2SNXV1Yq3t7d6TPeVK1eU\nwYMHKz179lSCgoKUq1evapywaVevXlVef/11xdXVVXFzc1OSk5ONKn9ERITi7u6u9OnTRxk/frxy\n69Ytg84fFhamWFlZ6Z2XdLe8ERERirOzs+Li4qIcPnxYi8h6Gss/d+5cxdXVVfH09FRGjhypXLt2\nTZ1naPmfRlKftCG1qXlJbTI+Upu0I/Wp+RhDbbqvy+UbkzvvG2Lo9zKqb/bs2QwbNoz//Oc/ZGdn\n4+rqajT5CwoK2LBhA+np6eTk5FBTU8OuXbsMOv+D3PsqNzeX6OhocnNziY+PZ8aMGdTW1moRW9VY\n/iFDhnD69GmysrLo1asXH330EWCY+Z9GUp+an9Sm5ie1yfhIbdKG1KfmZRS1qVmGf82kqKhICQwM\nVI4cOaJu9XFxcVEuXryoKIqiXLhwQXFxcdEyYpOuXbumODo6NnjeWPJfuXJF6dWrl1JeXq788ssv\nSkhIiJKQkGDw+e+8kl9TeSMiIpTIyEi13dChQ5UTJ040b9hG3Jm/vq+++kp58803FUUx3PxPE6lP\n2pDapA2pTcZDapN2pD41P0OvTS1qj1lj9w0xlnsZnT9/HktLS8LCwujbty9Tp06lsrLSaPI/99xz\nvP/++9jb29O1a1c6depEUFCQ0eSv01Te0tJSbG1t1Xa2traUlJRokvF+bdq0Sb03jjHmb2mkPmlD\napPhkdpkWKQ2aUfqk2ExhNrUYgZmBw8exMrKCp1O1+TlZU1MTNTd9Iamurqa9PR0ZsyYQXp6Ou3b\nt2+w69qQ8587d45PP/2UgoICSktLqaioYMeOHXptDDl/Y+6V15DXZeXKlTzzzDO88cYbTbYx5Pwt\njdQn7UhtMixSmwyL1CZtSX0yHIZSm1rMwKzuviGOjo6MHz+eI0eOMHHiRPVeRsBd72WkNVtbW2xt\nbenXrx8Ao0ePJj09HRsbG6PIn5qaysCBA7GwsMDU1JRRo0Zx4sQJo8lfp6m/l/u9X58h2LJlC3Fx\ncezcuVN9zpjyt0RSn7QjtclwSG0yPFKbtCX1yTAYUm1qMQOzxu4bsn37dqO5l5GNjQ12dnbk5eUB\n8PXXX9O7d2+GDx9uFPldXV1JTk7m5s2bKIrC119/jbu7u9Hkr9PU30toaCi7du3i559/5vz585w9\nexY/Pz8tozYqPj6eVatWsX//fszMzNTnjSV/SyX1STtSmwyD1CbDJLVJW1KftGdwtemJn8WmgaSk\nJGX48OGKotw+sTIwMNAgLzl6p8zMTMXX11fvkp3GlP/Pf/6zesnXt956S/n5558NOv+4ceOULl26\nKG3atFFsbW2VTZs23TXvypUrFScnJ8XFxUWJj4/XMPltd+bfuHGj4uzsrNjb2yve3t6Kt7e3Mn36\ndLW9oeV/Wkl9an5Sm5qX1CbjJLVJG1Kfmo8x1CYTRWnioGIhhBBCCCGEEM2ixRzKKIQQQgghhBDG\nSgZmQgghhBBCCKExGZgJIYQQQgghhMZkYCaEEEIIIYQQGpOBmdDc2LFjOXfuHAAVFRVMnz4dZ2dn\nfHx88PX1JSoqCoCCggLatm2LTqejd+/eTJ8+vckbYgKsWbOG7du3N8s6CCFaHqlNQrRsrVu3RqfT\n4eHhwZgxY7h586bWke7LTz/9xKuvvqo+TklJYdCgQfTq1QsfHx9CQkL4/vvvAVi6dCm2trbqeh44\ncOCufQcGBnLjxo0nml80TQZmT5Ha2lqDW1Z+fj6VlZU4OTkBEB4ejoWFBfn5+aSlpREfH095ebna\n3tnZmYyMDLKzs8nNzWXfvn1N9h0WFsbf//73R1sRIcQTJ7VJCKGFdu3akZGRQU5ODs888wzr1q3T\nm19dXd1sWR5kWZ999hmTJ08GoKysjLFjxxIZGUleXh5paWksXLhQ3ahkYmLCnDlzyMjIICYmhilT\npty173HjxrFhw4aHXg/xaGRg1gIUFBTg6urKhAkTcHd357e//a261cfBwYEFCxbg4+NDTEwMCQkJ\nDBw4EB8fH8aMGUNlZSUACxYsoHfv3nh5efHBBx8AEBMTg4eHB97e3gwaNAi4fXf0d955R112SEgI\nR48eBeDZZ59l7ty5eHt7c+LECXbs2IG/vz86nY5p06Y1+oNo165dhIaGAnDu3DlOnTrFihUr1PnP\nP/+8mqe+1q1bM3DgQPLz84mKisLPzw9vb29Gjx6trru5uTkWFhacPn36Ud9iIcRDkNoktUkIYxEQ\nEEB+fj7ffvstAQEBjBgxgj59+lBbW8u8efPw8/PDy8uL9evXA3DhwgVeeukldU/U8ePHqa2tZfLk\nyXh4eODp6cnq1asBGDRoEGlpaQBcvnwZR0dH4HbdCg0NJTAwkKCgIKqqqpgyZQr+/v707duX2NjY\nRrPu2bNH3WNWN0jr37+/Ov+FF15gxIgR6uO6Pfiurq6Ymppy6dIlRo4cia+vL3369NEbiNXdWFlo\nQwZmLUReXh4zZ84kNzeXDh06sHbtWuD2lpLnn3+etLQ0AgMDWblyJYmJiaSlpeHj48Mnn3xCeXk5\n+/bt4/Tp02RlZbF48WIAli9fTkJCApmZmWpxMDEx0Vtu/cdVVVX079+fzMxMnnvuOXbv3s13331H\nRkYGrVq1YufOnQ1yHz9+HF9fXwBOnz6Nl5fXfa1vVVUViYmJeHp6MmrUKFJSUsjMzMTNzY2NGzeq\n7fz8/NQfZ0KI5ie1SWqTEIauurqauLg4PD09AcjIyGDNmjWcOXOGqKgoOnXqREpKCikpKWzYsIGC\nggL+9a9/ERwcrO4p9/LyIiMjg9LSUnJycsjOziYsLAy4XY/urFF1MjIy+PLLL/nmm29YsWIFgYGB\nnDx5kiNHjjBv3jyqqqr02l+8eJHWrVvTrl07AHJzc+nbt+99refJkydp3bo1lpaWbNq0idTUVE6d\nOsWaNWvUIwCsra25fPmyunFMNC8ZmLUQdnZ2DBgwAIAJEyZw7Ngxdd7YsWMBSE5OJjc3l4EDB6LT\n6di2bRuFhYV07NgRMzMzfve737F3717atm0L3N7iMmnSJKKiou5rF3vr1q15/fXXAdQfWL6+vuh0\nOo4cOcL58+cbvObHH3+kS5cuQMMfVhEREeh0Orp166Y+d+7cOXQ6HS+++CIhISEMHTqUnJwcAgIC\n8PT0ZOfOnXpbobt27UpBQcH9vIVCiCdAapPUJiEM1c2bN9HpdPTr1w8HBwemTJmCoij4+fnRvXt3\nABISEti2bRs6nY7+/ftTXl5Ofn4+/fr1Y/PmzSxbtozs7GyeffZZnJyc+OGHH5g1axaHDx/G3Nz8\nnhmCgoLo1KmTuqzIyEh0Oh0vv/wy//vf/ygqKtJrX7821al/Tqu/vz/u7u68++676ry//e1v6HQ6\n5s2bR3R0NACrV6/G29ubAQMGUFxczNmzZ9U+rK2tGyxXNA9TrQOIx6P+DwdFUfQet2/fXv13UFAQ\nX3zxRYPXp6SkkJiYyJ49e/jss89ITEzkn//8JykpKRw6dAgfHx/S0tIwNTXVO+zn1q1b6r/NzMz0\nljtp0iQiIiLumb2uoLi5uZGVlaXmX7RoEYsWLdIrbE5OTmRkZOi9fvLkycTGxuLh4cHWrVtJSkpq\n8r0QQjQvqU1Sm4QwVG3btm3wvQX92gS3DxcMCgpq0O7f//43Bw8eZPLkycyZM4eJEyeSlZXF4cOH\nWbduHbt372bjxo169al+bWpsWV999RU9e/a8a+76A7HevXuTnp6uHnp98uRJvvzySw4ePAj8eo7Z\nnDlz1NckJSWRmJhIcnIyZmZm6iCwfv9Sn7Qhe8xaiMLCQpKTkwH44osvCAgIaNDG39+f48ePqyeE\nVlZWcvbsWSorK7l27Rq/+c1v+OSTT8jKygJubwH28/Nj2bJlWFpaUlxcjIODA5mZmSiKQlFRESkp\nKY3mCQwMZM+ePVy6dAmA8vJyCgsLG7Tr3r07Fy5cAG6fPO/r68vixYvVAnbz5s27Xt0Mbl8tzcbG\nhl9++YUdO3bozbtw4QIODg53fb0Q4smR2iS1SQhjNnToUNauXavunc/Ly6OqqorCwkIsLS0JDw8n\nPDyc9PR0rly5Qk1NDaNGjWL58uXqoM/BwYHU1FTg9vlhd1vWmjVr1MeNDRq7d+/OxYsX1cczZ85k\ny5YtnDhxQn2usrJSHVgpitKgVl2/fp3OnTtjZmbGmTNn1Bpdp6ysDFtb2/t6f8TjJXvMWggXFxf+\n8Y9/MGXKFPVyzaC/tdrS0pItW7Ywfvx4dcvIypUrMTc3Z8SIEdy6dUvd5Q3wwQcfcPbsWRRFYfDg\nweqx146Ojri7u+Pm5oaPj4/af/1lubm5sWLFCoYMGUJtbS1t2rRh7dq12Nvb6+V+8cUXSU1NVfuJ\niopi3rx5ODs7Y2FhQdu2bVm1alWjy6izfPly/P39sbS0xN/fn4qKCnVeSkoKH3/88cO9qUKIRya1\nSWqTEIaqse/tneeDhYeHU1BQQN++fVEUBSsrK/bu3UtSUhKrVq2iTZs2mJubs23bNkpKSggLC1M3\n4ERGRgIwd+5cxowZw/r163n11VfV/u9c1pIlS3j33Xfx9PSktraWHj16NLgAiI2NDdXV1VRWVtK+\nfXusra2Jjo5m/vz5lJSUYGVlhaWlJX/84x8bXQZAcHAw69atw93dHRcXF/Vwc/hdbskAAAD5SURB\nVLh9DpuFhUWDPXmieZgo99rkJwxeQUEBw4cPJycnR+soD+yHH37gnXfe4dChQ4+97+vXrxMYGMip\nU6cee99CiHuT2tQ4qU1CiEexdOlS3Nzc1PN0H6f169dTWVnJe++999j7FvcmhzK2EMZ6LHCPHj0w\nNzdXD2F6nLZs2cLs2bMfe79CiPsntakhqU1CiEcxc+ZMtm7d+kT6jo6OZurUqU+kb3FvssdMCCGE\nEEIIITQme8yEEEIIIYQQQmMyMBNCCCGEEEIIjcnATAghhBBCCCE0JgMzIYQQQgghhNCYDMyEEEII\nIYQQQmMyMBNCCCGEEEIIjf0fZg+Xlwsa740AAAAASUVORK5CYII=\n" } ], "prompt_number": 9 } ], "metadata": {} } ] }