diff --git a/MT Script - 3D.ipynb b/MT Script - 3D.ipynb new file mode 100644 index 00000000..65b6b8b3 --- /dev/null +++ b/MT Script - 3D.ipynb @@ -0,0 +1,210 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from SimPEG import *\n", + "from scipy.constants import mu_0\n", + "\n", + "\n", + "def omega(freq):\n", + " \"\"\"Change frequency to angular frequency, omega\"\"\"\n", + " return 2.*np.pi*freq" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 35 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "M = Mesh.TensorMesh([[(100.,32)],[(100.,34)],[(100.,18)]], x0='CCC')" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 18 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "print M" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + " ---- 3-D TensorMesh ---- \n", + " x0: -1600.00\n", + " y0: -1700.00\n", + " z0: -900.00\n", + " nCx: 32\n", + " nCy: 34\n", + " nCz: 18\n", + " hx: 32*100.00\n", + " hy: 34*100.00\n", + " hz: 18*100.00\n" + ] + } + ], + "prompt_number": 25 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "conds = [1e-2,1]\n", + "sig = Utils.ModelBuilder.defineBlock(M.gridCC,[0,0,-300],[500,500,-100],conds)\n", + "sig[M.gridCC[:,2]>0] = 1e-8\n", + "sigBG = np.zeros(M.nC) + conds[0]\n", + "sigBG[M.gridCC[:,2]>0] = 1e-8\n", + "colorbar(M.plotImage(log10(sig)))" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "metadata": {}, + "output_type": "pyout", + "prompt_number": 34, + "text": [ + "" + ] + }, + { + "metadata": {}, + "output_type": "display_data", + "png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAEKCAYAAAAPVd6lAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XtcVHX++PHXACNpiiYpCKNiiAmIiBe02lZbQw1X80Kk\ntuK90lyvpXYV259CtV00pXK/WGYmWrsJ3tg0b5WKF9RUTKcE5V6KIISCjOf3BzJJDDDAXNn38/E4\nj5hzPucz7zdDb898zjmfo1IURUEIIYTdcrB2AEIIIRpGCrkQQtg5KeRCCGHnpJALIYSdk0IuhBB2\nTgq5EELYOSnkQghh56SQC5v3+eef07t3b1q0aIGHhwehoaF8//33AJw/f54nnniCNm3a0KpVKwID\nA3n33Xe5deuWlaMWwnKkkAub9s477zB37lxeeeUVfvnlF9LT03nuuedISEjg559/pm/fvnTs2JHT\np0+Tn5/PF198wbFjxygsLLR26EJYjEru7BS2qqCgAI1GwyeffMLo0aOrbP/b3/5GQUEBW7ZssUJ0\nQtgOOSIXNuvgwYPcuHGDkSNHGtz+zTffEBYWZuGohLA9UsiFzbpy5Qr33nsvDg6G/0yvXLlCu3bt\nLByVELZHCrmwWa6urly+fLnaE5eurq5kZWVZOCohbI8UcmGzHnjgAZydnfnqq68Mbn/00Uf597//\nbeGohLA9UsiFzWrZsiWvv/46zz33HPHx8RQXF3Pz5k127NjBwoULWbJkCQcOHGDBggXk5uYC8NNP\nPzF+/HgKCgqsHL0QliOFXNi0efPm8c477/D//t//o23btnTo0IGYmBhGjhzJfffdx8GDB0lLS8Pf\n359WrVoRFhZGnz59aNGihbVDF8Ji5PJDIYQwscTERObMmYNOp2Pq1KksXLiwSptZs2axY8cOmjVr\nxieffEJQUFC930+OyIUQwoR0Oh0zZ84kMTGRlJQUNmzYwNmzZyu12b59Oz/99BNarZbVq1czffr0\nBr2nFHIhhDChw4cP07lzZ7y8vFCr1YwZM4b4+PhKbRISEpgwYQIAffv2JT8/X3+epz6kkAshhAll\nZmbSvn17/WuNRkNmZmatbTIyMur9nlLIhRDChFQqlVHt/nh60tj9DHGq954W1KNHD06ePGntMIQQ\ndiAwMJATJ040qI+mKhU3jGzbvHnzSpO0eXp6kp6ern+dnp6ORqOptM8f22RkZODp6VnveO2ikJcX\n8UhrhwHsAR4xSU+KshgAlWqJ1WJomKpx1D8n08VgDhV5VScyMpLIyEizx2HrMdhKHA05sq1wA/h/\nRrZ9paio0uvevXuj1WpJS0vDw8ODjRs3smHDhkpthg8fzsqVKxkzZgyHDh2iVatWuLm51Tteuyjk\nQghhaep67ufk5MTKlSsZPHgwOp2OKVOm4Ovry0cffQTAM888Q2hoKNu3b6dz587cfffdfPzxxw2K\nVQq5EEIY0JDi+Nhjj/HYY49VWvfMM89Uer1y5coGvENlUsjrxMvaAWAbMYBtxOFl7QAAGDBggLVD\nsIkYwHbiMIWm1g6gDqSQ10knaweAbcQAthGHLcRgG8XLFmIA24nDFOo7tGINUsiFEMIAeyqO9hSr\nEEJYjByRCyGEnbOn4mhPsQohhMXIEbkQQtg5KeRCCGHn5PJDIYSwc/ZUHO0pViGEsBgZWhFCCDtn\nT8XRnmIVQgiLkSNyIYSwc/ZUHO0pViGEsBh7OiKXR70JIYQBTY1c6iIvL4+QkBC6dOnCoEGDyM/P\nN9guKioKf39/AgICGDduHCUlJTX2K4VcCCEMUBu51EV0dDQhISGcP3+egQMHEh0dXaVNWloa//rX\nv0hOTubUqVPodDri4uJq7FcKuRBCGOBk5FIXCQkJTJgwAYAJEyawefPmKm1cXFxQq9UUFxdTVlZG\ncXFxrc/zrLGQT548GTc3NwICAvTrXnjhBXx9fQkMDGTUqFEUFBTot0VFReHj40PXrl35+uuv9euP\nHTtGQEAAPj4+zJ49W7++pKSEJ598Eh8fH/r168fFixdr+TUIIYRlqJ2MW+oiNzdX/2xONzc3cnNz\nq7Rp3bo18+fPp0OHDnh4eNCqVSseffTRGvutsZBPmjSJxMTESusGDRrEmTNnOHnyJF26dCEqKgqA\nlJQUNm7cSEpKComJicyYMQNFUQCYPn06sbGxaLVatFqtvs/Y2FhcXV3RarXMnTuXhQsXGvnrEEII\n83JyMrwcVMFbyu/LH4WEhBAQEFBlSUhIqNROpVIZfFD0zz//zHvvvUdaWhpZWVkUFRWxfv36mmOt\naePDDz9MWlpalSAr9O3bl3//+98AxMfHM3bsWNRqNV5eXnTu3JmkpCQ6duxIYWEhwcHBAERERLB5\n82aGDBlCQkICS5aUP3F99OjRzJw5s8ZghRDCUtSOhtf/xRH+csfr6N8qb9+5c2e1fbq5uZGTk4O7\nuzvZ2dm0bdu2SpujR4/y4IMP4urqCsCoUaM4cOAATz31VLX9NmiMfM2aNYSGhgKQlZWFRqPRb9No\nNGRmZlZZ7+npSWZmJgCZmZm0b98eKH/ydMuWLcnLy2tISEIIYRLVHZH/camL4cOHs3btWgDWrl3L\niBEjqrTp2rUrhw4d4vr16yiKwq5du/Dz86ux33oX8qVLl9KkSRPGjRtX3y6EEMJmqZ2NW+pi0aJF\n7Ny5ky5durB7924WLVoElB8IDx06FIDAwEAiIiLo3bs33bt3B+Dpp5+usd963RD0ySefsH37dr75\n5hv9Ok9PT9LT0/WvMzIy0Gg0eHp6kpGRUWV9xT6XLl3Cw8ODsrIyCgoKaN26dTXvuueOn72wlQfv\nCiGsa+/evezdu9f0HZvhdsnWrVuza9euKus9PDzYtm2b/vWCBQtYsGCB0f3WOdTExETeeust9u3b\nx1133aVfP3z4cMaNG8e8efPIzMxEq9USHByMSqXCxcWFpKQkgoODWbduHbNmzdLvs3btWvr168eX\nX37JwIEDa3jnR+oaqhDif8CAAQMYMGCA/nXFebcGs6P73msMdezYsezbt4/Lly/Tvn17lixZQlRU\nFKWlpfqTng888AAxMTH4+fkRHh6On58fTk5OxMTE6M/IxsTEMHHiRK5fv05oaChDhgwBYMqUKYwf\nPx4fHx9cXV1rvehdCCEsxo4KuUqpuEbQhpX/gxBp7TBMSlEWA6BSmejowQY0xpzg97yEfVCpVDS0\nrKlUKhQfI9tqafD7NZQd/ZsjhBAWZEfV0Y5CFUIIC6rjFSnWJIVcCCEMsaPqaEehCiGEBdlRdbSj\nUIUQwoKquUXfFkkhF0IIQ+yoOtpRqEIIYUF2VB3tKFQhhLAgO6qOdhSqEEJYkFx+KIQQds6OqqM8\ns1MIIQxxNHKpgy+++AJ/f38cHR1JTk6utl1+fj5hYWH4+vri5+fHoUOHauxXCrkBOTnz6dmzHQD7\n9k1kzJhu+m1+fm3YtCmMc+dmUlb2KqtXD6uyf//+HdHpXquyTJrUw2I5/FFDcwJwcFCxcOFD/Pjj\nc1y//jI5OfNZtSrUIvFXp6F5ffzx4wY/q7KyV3F1bWqxPIQNMsPTlwMCAvjqq6/485//XGO72bNn\nExoaytmzZ/nhhx/w9fWtNVRxB2/ve2jWTM3x49mo1Q707u3Bd99d0m9v2tSJtLQC4uPPMW/eAzVO\nlhMU9BHZ2YX619eulZg19uqYKqdPPnmcvn01LFiwkxMncmjRwhkvr1aWSqMKU+Q1a9YOFiz4/dFc\nKpWKzZufpKiolCtXrlskD2GjzFAdu3btWmubgoICvv32W/2ThCqenlYTKeR/8NBDHUhKykRRoE8f\nT65cKSYj45p++7Fj2Rw7lg3AlClBNfZ1+XIxv/5abNZ4jWGKnAYM8GLMmG507/4hP/54Wb/+9Olf\nzBt8DUyRV2FhKYWFpfrXPj6t6dtXwxNPfGHe4IXts9INQampqbRp04ZJkyZx8uRJevXqxfLly2nW\nrFm1+0ghv+3q1YUoioKzsxMODiry8hagVjvi7OxIXt4CFAVcXd+sU5/ffTeJZs3U/PRTHh99dIx1\n634wU/SGmTKn0aN9uXDhKoMGebNly1iaNHHk4MF0nn9+Z6XiaQnm+KwqPPNMb3Jyiti8+UcTRy3s\nTjXVcW8W7M2ufreQkBBycnKqrF+2bBnDhhketrxTWVkZycnJrFy5kj59+jBnzhyio6N5/fXX6xrq\n/57u3T9ApVJx6NAUnn12GydO5BAXN5rPPz9NfHzd/qfOyipkxoxtHD2axa1bCqGhPqxePYzOnVuz\nePFe8yRggClz8va+hw4dWvK3vwUwZUoCpaU6li79C7t3R9Ct2weUlurMlEVVpszrTk2aODJhQiAf\nfXSMW7dsfpp+YW53GV494L7ypcKS45W379y5k4bQaDRoNBr69OkDQFhYGNHR0TXuI4X8tvT0awQE\ntEWtdmTLlnM0b96EHj3cGT48jsuX6zY8otXmodXm6V8fP56Do6MDzz//AEuW7LNYkTBlTg4OKpyd\nnYiI2KwfWnnyyS/Jzp5PaKiPRY9gTZnXncLC/LjnnrtYvfqYCaMVdsvMQyvVnYtyd3enffv2nD9/\nni5durBr1y78/f1r7EsKOXD69HQ6dGiJk5MDarUjBQWL9IXrwoXy54v6+q4iM7Owlp6ql5SUwd13\nN6FNm2bk5v5mqtCrZeqcsrOLUBSl0vj45cvFXL5cTIcONZ+IMSVzflbPPtuL//73Zy5dKjB12MIe\nmaE6fvXVV8yaNYvLly8zdOhQgoKC2LFjB1lZWUybNk3/AOb333+fp556itLSUry9vfn4448tHar9\nGTJkPU2aOLJmzXB27PiJTZvOsHhxf0pKdERHfweUF7KG6NmzHcXFNxt0xFgXps5p//6LREQE0qWL\nK+fPXwGgdeum3HtvM9LS8s2SgyHm+qx8fe/loYc6MHLkRlOHLOyVGarjyJEjGTlyZJX1Hh4e+iIO\nEBgYyJEjR4zuVwo5kJFxDQcHFd27u/H001tJTc0nIMCNyMi9pKZWLlJOTg74+7cBoEULZ1xdmxIY\n6EZpqY6zZ8uPVufM6cfFi/mkpPyKosDgwd68/PLDrFx5BJ3OMsMqps5pw4bTvPzyw6xZM5xZsxK5\neVPHG288ilZ7hR07tBbJyRx5VXjmmd5kZRWyZcs5i+UibJxMY2t/goLcKSnRcf78FVxcnPH3b8P+\n/RertPP0bEFy8jNA+RhXz57tGDnSl7S0fLy9VwDg6Khi2bKBtG/vws2bt9BqrzBrViJr1hyv0p+9\n5HTjRhmPPrqO994bzN69EyguvsnevWk8+ug6bt68Zbd5Adx1lxPjx3fn/fcPY/uPIhcWY0fVUaVY\n+/HPRlCpVECktcMwqcb4xPnGmBP8npewDyqVqsFPtVepVCjVX+1Xue1r1Z+4tBQ7+jdHCCEsSGY/\nFEIIO2dH1bHGoZXJkyezbds22rZty6lTpwDIy8vjySef5OLFi3h5ebFp0yZatSqfbyMqKoo1a9bg\n6OjIihUrGDRoEADHjh1j4sSJ3Lhxg9DQUJYvXw5ASUkJERERJCcn4+rqysaNG+nYsWPVIE3wVUkI\n8b/BZEMry41sO9v6Qys1zn44adIkEhMTK62Ljo4mJCSE8+fPM3DgQP0dRykpKWzcuJGUlBQSExOZ\nMWOGPrnp06cTGxuLVqtFq9Xq+4yNjcXV1RWtVsvcuXNZuHChOXIUQoi6M8M0tuZS45eHhx9+mLS0\ntErrEhIS2LdvHwATJkxgwIABREdHEx8fz9ixY1Gr1Xh5edG5c2eSkpLo2LEjhYWFBAcHAxAREcHm\nzZsZMmQICQkJLFlSfmJs9OjRzJw5s9pYGusJtMaUV2PMCRpnXo0xJzDxiWk7Glqp83zkubm5uLm5\nAeDm5kZubi4AWVlZaDQafTuNRkNmZmaV9Z6enmRmZgKQmZlJ+/btgd+naszL+/3WdiGEsBozzEdu\nLg0KQ6VS3b40UAghGhkbGTYxRp0LuZubGzk5Obi7u5OdnU3btm2B8iPt9PR0fbuMjAw0Gg2enp5k\nZGRUWV+xz6VLl/Dw8KCsrIyCggJat25dzTvvueNnL6BTXUMXQjRKqUAaAJGRJjzpWM3sh7aozkMr\nw4cP1z+5Yu3atYwYMUK/Pi4ujtLSUlJTU9FqtQQHB+Pu7o6LiwtJSUkoisK6det4/PHHq/T15Zdf\nMnDgwBre+ZE7FiniQogKnaioDZGRkabrtrEMrYwdO5Z9+/Zx+fJl2rdvz+uvv86iRYsIDw8nNjZW\nf/khgJ+fH+Hh4fj5+eHk5ERMTIx+2CUmJoaJEydy/fp1QkNDGTJkCABTpkxh/Pjx+Pj44OrqSlxc\nnJnTFUIII9nR0Ircom8ljfGqgcaYEzTOvBpjTlCel8muI08wsu1w468jf+GFF9i6dStNmjTRT09b\n3fM4dTodvXv3RqPRsGXLlhr7rfPQihBC/E8ww9DKoEGDOHPmDCdPnqRLly5ERUVV23b58uX4+fkZ\ndUGJFHIhhDDEDDcEhYSE4OBQXnb79u1b6UKQO2VkZLB9+3amTp1q1NG+FHIhhDDkLiOXelqzZg2h\noaEGt82dO5e33npLX/RrYyPnXIUQwsbU82RnSEgIOTk5VdYvW7aMYcOGAbB06VKaNGnCuHHjqrTb\nunUrbdu2JSgoiL179xr1nlLIhRDCkGqq495jsDe5+t127txZY7effPIJ27dv55tvvjG4/cCBAyQk\nJLB9+3Zu3LjBtWvXiIiI4NNPP61rqEII8T+umuo4oG/5UmFJrPFdJiYm8tZbb7Fv3z7uusvwuMyy\nZctYtmwZAPv27eOf//xnjUUcZIxcCCEMM8NVK3//+98pKioiJCSEoKAgZsyYAZTPVTV06FCD+xhz\n1YockQshhCFmuCFIqzX8oHIPDw+2bdtWZX3//v3p379/rf1KIRdCCEPsqDraUahCCGFB8sxOIYSw\nc3ZUHe0oVCGEsCA7qo52FKoQQliQHVVHOwpVCCEsR7GjaWylkAshhAE6O6qOdhSqEEJYjhRyIYSw\ncyXOTYxsWWrWOIwhhVwIIQzQOdrPILkUciGEMEBnRw/tlEIuhBAGlEkhF0II+6azo/Io09gKIYQB\nOhyNWuri1VdfJTAwkB49ejBw4EDS09OrtElPT+eRRx7B39+fbt26sWLFilr7lUIuhBAGmKOQL1iw\ngJMnT3LixAlGjBjBkiVLqrRRq9W8++67nDlzhkOHDrFq1SrOnj1bY7/2891BCCEsqARjLz80XosW\nLfQ/FxUVce+991Zp4+7ujru7OwDNmzfH19eXrKwsfH19q+233kfkUVFR+Pv7ExAQwLhx4ygpKSEv\nL4+QkBC6dOnCoEGDyM/Pr9Tex8eHrl278vXXX+vXHzt2jICAAHx8fJg9e3Z9wxFCCJPS4WTUUlcv\nv/wyHTp0YO3atSxatKjGtmlpaRw/fpy+ffvW2K5ehTwtLY1//etfJCcnc+rUKXQ6HXFxcURHRxMS\nEsL58+cZOHAg0dHRAKSkpLBx40ZSUlJITExkxowZKIoCwPTp04mNjUWr1aLVaklMTKxPSEIIYVL1\nHVoJCQkhICCgyrJlyxYAli5dyqVLl5g4cSJz586t9v2LiooICwtj+fLlNG/evMZY6zW04uLiglqt\npri4GEdHR4qLi/Hw8CAqKop9+/YBMGHCBAYMGEB0dDTx8fGMHTsWtVqNl5cXnTt3JikpiY4dO1JY\nWEhwcDAAERERbN68mSFDhtQnLCGEMJnqxr+P7v2No3uLq91v586dRvU/btw4QkNDDW67efMmo0eP\n5m9/+xsjRoyota96FfLWrVszf/58OnToQNOmTRk8eDAhISHk5ubi5uYGgJubG7m5uUD5g0X79eun\n31+j0ZCZmYlarUaj0ejXe3p6kpmZWZ+QhBDCpKq7jrzHABd6DHDRv1695LLRfWq1Wnx8fACIj48n\nKCioShtFUZgyZQp+fn7MmTPHqH7rNbTy888/895775GWlkZWVhZFRUV89tlnldqoVCqjnv4shBC2\nyBxj5C+++CIBAQH06NGDvXv38vbbbwPlB7tDhw4F4Pvvv+ezzz5jz549BAUFERQUVOuQc72OyI8e\nPcqDDz6Iq6srAKNGjeLgwYO4u7uTk5ODu7s72dnZtG3bFig/0r7zesmMjAw0Gg2enp5kZGRUWu/p\n6VnNu+6542cvoFN9QhdCNDqpQBoAkZGKyXo1xy36X375pcH1Hh4ebNu2DYA//elP3Lp1q0791uuI\nvGvXrhw6dIjr16+jKAq7du3Cz8+PYcOGsXbtWgDWrl2rH9sZPnw4cXFxlJaWkpqailarJTg4GHd3\nd1xcXEhKSkJRFNatW1fDeNAjdyxSxIUQFTpRURsiIyNN1mspTYxabEG9jsgDAwOJiIigd+/eODg4\n0LNnT55++mkKCwsJDw8nNjYWLy8vNm3aBICfnx/h4eH4+fnh5ORETEyMftglJiaGiRMncv36dUJD\nQ+VEpxDCJtjTXCsqpeI6QBtWXvQjrR2GSSnKYgBUqqp3dtmrxpgTNM68GmNOUJ6XSqWioWVNpVKx\nQxlgVNvHVHsb/H4NJXd2CiGEAfY0ja3MtfIHOTnz6dmzHQD79k1kzJhu+m1+fm3YtCmMc+dmUlb2\nKqtXDzPYh49PaxITn6Ko6EV++eV5YmKG0rSpdf/NbGhebm5389lnIzl1ajqlpa/w9dd/s1jsNWlo\nXsOH38+2bePIyppHUdGLnDo1nb//Pdhi8RvS0Jx69HBnz54JZGfP5/r1l0lLm82KFY/h4uJssRz+\nyBT/X1Vwc7ub7Oz56HSv0a5dzTfKNIQ55loxFynkd/D2vodmzdQcP56NWu1A794efPfdJf32pk2d\nSEsr4PXX93HyZK7Br1N3363mm28iKC3V8cADsYSHf8mQId7Exg63ZCqVmCIvZ2cnrly5zttvH2TX\nrgvYwoCcKfLq378j33+fzogRG/H3j+Gttw4QFTWQF1540JKp6Jkipxs3yliz5jghIevo3HkFU6Yk\nMGjQfaxdW/uNJeZgipwqqFSwfv0okpIyqm1jKmU4GrXYAhlaucNDD3UgKSkTRYE+fTy5cqWYjIxr\n+u3HjmVz7Fg2AFOmVL2QH2DcuABcXZsxbtx/KCoqf5bfc89tZ+vWcbz44jdcvFhg/kT+wBR5XbpU\nwOzZ5dey9u/fEU/PFgbbWZIp8po//+tKrz/99CQ9e7YjPNyft946YL7gq2GKnH788TI//vj7TSqZ\nmYXExBxl8eL+5g2+GqbIqcKrr/bnxo0y3n33EMOG3W/WuEux3jeYupJCDly9uhBFUXB2dsLBQUVe\n3gLUakecnR3Jy1uAooCr65tG9fXQQ+05cCBdX8QBdu68wK1bCg8+2N6ihdyUedkSc+d1zz13Vfr8\nLMGcOWk0LoSF+bJjh9bEUdfM1DkNGODF1KlBBAV9RLdubc0YeTlbGTYxhhRyoHv3D1CpVBw6NIVn\nn93GiRM5xMWN5vPPTxMf/2Od+mrXrgU5OUWV1pWV3SIv7zrt2ln2KNaUedkSc+bVv39HxozpxsiR\nG00UrXHMkdP330+mRw937rrLif/+9yemTEkwcdQ1M2VObdvezbp1I4mI+IorV66bKeLKbGXYxBgy\nRg6kp1+jZUtn1GpHtmw5x9Wr1+nRw524uNOkp18jPf1a7Z3cZu3LkO5kyrxsibny6tvXk6++epLF\ni/eyfbtlj17NkVN4+BcEBX3E6NGb6NixFRs3hpkh8uqZMqf160fx6acn2bMnrdJ6c04DYq5pbM3B\nNqKwotOnp9OhQ0ucnBxQqx0pKFiEg4MKZ2cnLlyYBYCv7yoyMwuN6i87u4j27V0qrXNycqB166Zk\nZxvXhymYOi9bYa68+vfvSELCWJYt+5Y33vjeHKFXy1w5VbQ/f/4K2dmFHDgwha5d7600fm4ups7p\nL3/pRP/+HfUnoSsKeFrabP7v/44zY8Y2k+cgQyt2ZMiQ9TRp4siaNcPZseMnNm06w+LF/Skp0REd\n/R1QXpyN9f336SxfPoTmzZvox1lDQu7DwUHF999XfT6fuZg6rz+y1hcPc+QVGurDpk1hvPLKHt57\n75A5wq6RuT8rAEfH8i/fTk6W+RJu6py6dYup9Do42JM1ax5n0KDPOHv2V5PGXkEKuR3JyLiGg4OK\n7t3dePrpraSm5hMQ4EZk5F5SU/MrtXVycsDfvw0ALVo44+ralMBAN0pLdZw9W36U8/nnp3j11T/z\n+eejePnl3bi6NmPVqlDi4k5z6ZLlTnSaOi+AwMDyKYpbt25KixZN6N7dDZUKTp7Mtdu8wsL8WL9+\nFMuWfcvnn5/Cze1uAHQ6hcuXq59z2pZzmjIliKtXb5CS8is3bpTRrVtb3njjUY4dy+L06V/sMqc7\n/w6hfMwc4Ny5y+Tm/maWHKSQ25mgIHdKSnScP38FFxdn/P3bsH//xSrtPD1bkJz8DFA+Ft6zZztG\njvQlLS0fb+/yJ10XF9/k0Uc/5f33H+PgwSlcv17GF1+kMG/efy2aE5g2L0DfpqLd8ePPoCgKTk7/\nMH8ydzBlXjNm9MbRUcVrr/Xntdd+vzzvj7mbmylzKiu7xcsvP4y39z04OTmQnn6N//znrMUvpzT1\n398fmft8VIkdXX4oc61YSWOc66Ix5gSNM6/GmBOYdq6VN5W/G9V2ger9Or/f22+/zQsvvMDly5dp\n3bp1le2JiYnMmTMHnU7H1KlTWbhwYY39yVUrQghhgLlu0U9PT2fnzp107NjR8PvqdMycOZPExERS\nUlLYsGEDZ8+erbFPKeRCCGGAuW7RnzdvHm++Wf2NUIcPH6Zz5854eXmhVqsZM2YM8fHxNfZpN2Pk\nFV8FG5vGmFdjzAkaZ16NMSdTMcc14vHx8Wg0Grp3715tm8zMTNq3b69/rdFoSEpKqrFfuynkQghh\nSdUNm6TtvcjFvVVP2lYICQkhJyenyvqlS5cSFRXF11//Pr+PobH1+tzkZDeFfEkje5Dz4tsfYGPK\nqzHmBI0zr4qcGuPJTlOprpC3H3Af7Qfcp3+9f8l3lbbv3LnT4H6nT58mNTWVwMBAoPwZxb169eLw\n4cP65xu3d6DWAAAgAElEQVRD1Wccp6eno9FoaozVbgq5EEJYUomJn8fZrVs3cnN/v+eiU6dOHDt2\nrMpVK71790ar1ZKWloaHhwcbN25kw4YNNfYtJzuFEMIAc8+1cucQSlZWFkOHDgXAycmJlStXMnjw\nYPz8/HjyySfx9fWtsS85IhdCCAPMfWfnhQsX9D97eHiwbdvv88U89thjPPbYY0b3JYVcCCEMkFv0\nhRDCztnTfORSyIUQwgBbmWvcGPYTqRBCWJA9Da3U+6qV/Px8wsLC8PX1xc/Pj6SkJPLy8ggJCaFL\nly4MGjSI/Pzfp6uMiorCx8eHrl27Vrog/tixYwQEBODj48Ps2bMblo0QQphIKU2MWmxBvQv57Nmz\nCQ0N5ezZs/zwww907dqV6OhoQkJCOH/+PAMHDiQ6OhqAlJQUNm7cSEpKComJicyYMUN/R9P06dOJ\njY1Fq9Wi1WpJTEw0TWZCCNEA5pprxRzqVcgLCgr49ttvmTx5MlB+3WPLli1JSEhgwoQJAEyYMIHN\nmzcD5fMLjB07FrVajZeXF507dyYpKYns7GwKCwsJDg4GICIiQr+PEEJYkz09s7NehTw1NZU2bdow\nadIkevbsybRp0/jtt9/Izc3Fza38KTJubm76u5iysrIq3WKq0WjIzMysst7T05PMzMyG5COEECZh\nrmlszaFe/5yUlZWRnJzMypUr6dOnD3PmzNEPo1RQqVQmfcL1njt+9gI6maxnIYR9SwXSAIiMNN1z\ncmylSBujXkfkGo0GjUZDnz59AAgLCyM5ORl3d3f9rF/Z2dn6iWD+OAlMRkYGGo0GT09PMjIyKq33\n9PQ0+J6P3LFIERdC/K4TFdUhMjLSZL02+jFyd3d32rdvz/nz5wHYtWsX/v7+DBs2jLVr1wKwdu1a\nRowYAcDw4cOJi4ujtLSU1NRUtFotwcHBuLu74+LiQlJSEoqisG7dOv0+QghhTfY0Rl7vKN5//32e\neuopSktL8fb25uOPP0an0xEeHk5sbCxeXl5s2rQJAD8/P8LDw/Hz88PJyYmYmBj9sEtMTAwTJ07k\n+vXrhIaGMmTIENNkJoQQDWArlxYao96FPDAwkCNHjlRZv2vXLoPtX3rpJV566aUq63v16sWpU6fq\nG4YQQpiFrQybGEOmsRVCCAPMObTy9ttv4+DgQF5ensHtUVFR+Pv7ExAQwLhx4ygpKamxPynkQghh\ngLkuP0xPT2fnzp107NjR4Pa0tDT+9a9/kZyczKlTp9DpdMTFxdXYp22M1AuzWryy5u1LZlomDlNr\nrHkJ22Cuyw/nzZvHm2++yeOPP25wu4uLC2q1muLiYhwdHSkuLq72ar4KckQuhBAGmOOIPD4+Ho1G\nQ/fu3att07p1a+bPn0+HDh3w8PCgVatWPProozX2K0fkQghhQAnO9dovJCREfz/NnZYuXUpUVFSl\nSQMr5py6088//8x7771HWloaLVu25IknnmD9+vU89dRT1b6nFHIhhDCguqPt4r1HKN57tNr9du7c\naXD96dOnSU1NJTAwECi/AbJXr14cPnxYf/MkwNGjR3nwwQdxdXUFYNSoURw4cEAKuRBC1FV1hdx5\nQD+cB/TTv85b8qFR/XXr1k0//xRAp06dOHbsGK1bt67UrmvXrvzjH//g+vXr3HXXXezatUs/sWB1\nZIxcCCEMMPct+nfORZWVlcXQoUOB8nt0IiIi6N27t34s/emnn66xLzkiF0IIA8x9+/2FCxf0P3t4\neLBt2zb96wULFrBgwQKj+5JCLoQQBtjT7IdSyIUQwgAp5MKmNNYbYxprXsI2lJT+D0yaJYQQjZmu\nzH7Ko/1EKoQQFqQrk6EVIYSwa1LIhRDCzpXdlEIuhBB27ZbOfsqj/UQqhBCWJEMrQghh527YT3m0\nn0iFEMKSyqwdgPGkkAshhCFSyIUQws7ZUSGXaWyFEMKQm0YudRAZGYlGoyEoKIigoCASExMNtsvP\nzycsLAxfX1/8/Pw4dOhQjf1KIa+n+Tk5tOvZE4CJ+/bRbcwY/bbACRN4Taersng98oi1wjVKTTkB\nODVtysCoKGZduMDLN24wNz2dP7/yijVCrZOa8pqwZ4/Bz+rFwkJrhWuU2j6r4JkzmXHmDC8WFTEv\nM5PHP/6YZm3aWCPUOsnJmU/Pnu0A2LdvImPGdNNvc3RU8cILD3L27HMUF7/EuXMzmT69t/mC0Rm5\n1IFKpWLevHkcP36c48ePM2TIEIPtZs+eTWhoKGfPnuWHH37A19e3xn4bNLSi0+no3bs3Go2GLVu2\nkJeXx5NPPsnFixfx8vJi06ZNtGrVCoCoqCjWrFmDo6MjK1asYNCgQQAcO3aMiRMncuPGDUJDQ1m+\nfHlDQrKIe7y9UTdrRvbx4zio1Xj07s2l776r1OaWTsc7Hh5wx+TxN65etXSoRqstJ5WDA+O2baNJ\n8+ZsffppLp87RzNXV5rde68Vo65dbXltHDkSB7Va/1rl4MC0I0f4uZojJVtQW07dxoxh0Ntvs/XZ\nZ7mwaxct27dn6IcfMvLTT1n/2GNWjLxm3t730KyZmuPHs1GrHejd24Pvvruk375kySNMm9aTadO2\ncPJkDg8+2J7Vq4dRWqojNva46QMy09CKoed03qmgoIBvv/2WtWvXAuDk5ETLli1r3KdBR+TLly/H\nz89P/6SL6OhoQkJCOH/+PAMHDiQ6OhqAlJQUNm7cSEpKComJicyYMUOfzPTp04mNjUWr1aLVaqv9\nqmFLOjz0EJlJSaAoePbpQ/GVK1zLyKjSrvjyZYp//VW/3Cqz3UG32nIKjIigXc+erH/sMS7s2sW1\n9HRyTpzgwq5dVoy6drXldSM/v9Jn5BYQgIunJ0c/NO7xXdZQW06effuS+8MPnPj4Y66lp5N+4ADJ\nq1fjWcvjwqztoYc6kJSUiaJAnz6eXLlSTEbGNf32CRMC+ec/D5CQcI6LFwvYsOE0//d/ybz88sPm\nCeiGkUsdvf/++wQGBjJlyhTy8/OrbE9NTaVNmzZMmjSJnj17Mm3aNIqLi2vss95H5BkZGWzfvp2X\nX36Zd955B4CEhAT27dsHwIQJExgwYADR0dHEx8czduxY1Go1Xl5edO7cmaSkJDp27EhhYaH+eXQR\nERFs3ry52q8b1rbw6lUURcHJ2RmVgwML8vJwVKtxdHZmQV4eKApv3n5gqoOjI3//6SfUTZty+dw5\nDv7zn2i3b7dyBlUZm5Pv6NFkHj7MA3Pn0n38eHQ3b5L6zTfsWrTIJr9p1OWzulOvZ58lOzmZ7ORk\nK0RdM2Nz+mnHDoKmTKHjn//Mxf37udvNDb8nnuD81q3WTsGgq1cXoigKzs5OODioyMtbgFrtiLOz\nI3l5C1AUcHV9E2dnR0pKKo9l3LhRRseOrdBoXCoVfZOo53FXSEgIOTk5VdYvXbqU6dOn89prrwHw\n6quvMn/+fGJjYyu/bVkZycnJrFy5kj59+jBnzhyio6N5/fXXq33PehfyuXPn8tZbb3Ht2u+/vNzc\nXNzc3ABwc3PTP2g0KyuLfv1+f1ipRqMhMzMTtVqNRqPRr/f09CQzM7O+IZndB927o1KpmHLoENue\nfZacEycYHRfH6c8/58f4eH27yz/+SPykSeScPImTszP+4eGM3bKFhKlTOfHxx1bMoCpjc7rH25tW\nXl4oOh2bwsJo0rw5g999lzGbN/NJ//5WzMAwY/O6U3N3d+4fNoztzz1n4WiNY2xOP3/9Nf+dM4e/\n/fe/qBwccHBy4vzWrSRMnWrF6KvXvfsHqFQqDh2awrPPbuPEiRzi4kbz+eeniY//Ud9ux46fmDUr\nmG++ucCZM78SHOzJ5MlBKIqCh0cLyxXyU3vh9N5qd9u5c6dR3U+dOpVhw4ZVWa/RaNBoNPTp0weA\nsLAw/ehGdepVyLdu3Urbtm0JCgpi7969BtuoVKpKDxdtqD13/OwFdDJZz8a7lp5O24AAHNVqzm3Z\nQpPmzXHv0YO44cMpvnxZ3y4zKan8q2/F68OHuat1ax5auNDmCrmxOakcykfhvhwzhpKCAgASJk9m\n2pEjuAUGknvypFXir46xed0paPJkbl6/zqnPP7dwtMYxNqcuw4Yx+N13+e/cuVz89ltcNBpC3nqL\nx9es4avx462YgWHp6dcICGiLWu3Ili3naN68CT16uDN8eByXL/8+pDB7diIffjiUEyeeRVEUMjML\n+b//S2bRooe4desIcI7IyJrHn+ukukLuO6B8qRC3xOgus7Ozadeu/GTuV199RUBAQJU27u7utG/f\nnvPnz9OlSxd27dqFv79/jf3Wq5AfOHCAhIQEtm/fzo0bN7h27Rrjx4/Hzc2NnJwc3N3dyc7Opm3b\ntkD5kXZ6erp+/4yMDDQaDZ6enmTcMbaXkZGBp6enwfe09vUe00+fpmWHDjg4OeGoVrOooACVgwNO\nzs7Muv0Q1VW+vhRW840iMymJgHHjLBlyreqSU1F2No5qtb6IA/yakgJAq44dbaqQ1+uzUqnoOW0a\np9av52Yt45HWUJecHn7pJX747DP9OP+vZ85QWlTEpP372fPaa+SnplozlUpOn55Ohw4tcXJyQK12\npKBgEQ4OKpydnbhwYRYAvr6ryMwsJD//BmPG/BtHx//Qtu3dZGcX3b5qRcWFC96AB5GRi1myxPjC\nWqM6XlpojIULF3LixAlUKhWdOnXio48+AspHLaZNm6Z/APP777/PU089RWlpKd7e3nxcywFgvQr5\nsmXLWLZsGQD79u3jn//8J+vWrWPBggWsXbuWhQsXsnbtWkaMGAHA8OHDGTduHPPmzSMzMxOtVktw\ncDAqlQoXFxeSkpIIDg5m3bp1zJo1qz4hmd36IUNwbNKE4WvW8NOOHZzZtIn+ixejKynhu9tfe4qy\ns6vdv13PnhRculTtdmuoS04X9+/noQULaNKiBaW3L81zvf9+APLT0qwSf3Xq81l1HjKElh06cOz2\n/1i2pk45qVQouspjycqtW7c3me5bsikMGbKeJk0cWbNmODt2/MSmTWdYvLg/JSU6oqPLr8TJzi6q\ntI9Op+jXjR3bjX370sjLu2764Op4aaExPv30U4PrPTw89EUcIDAwkCNHjhjdr0nu7Kz441i0aBHh\n4eHExsbqLz8E8PPzIzw8HD8/P5ycnIiJidHvExMTw8SJE7l+/TqhoaE2e6LzWkYGKgcH3Lp3Z+vT\nT5OfmopbQAB7IyOrHOH0X7yYzKQkrmi1ODk74xcWRtDkyez4+9+tFL1hdcnpSEwMwTNnMvLTT9n9\n8suo776b0FWrSNu7l9wffrBSBobVJa8KvZ55hszDh20ulwp1yenH//yHP7/2GplHjnDp9tDK4Pfe\nI+fkSa7ePnq3FRkZ13BwUNG9uxtPP72V1NR8AgLciIzcS2pq5Ss6evVqR6dO95CcnE3btnczf/4D\ndO/uxp/+ZKbhStu9yKyKBhfy/v370//2ya7WrVuzq5rL0V566SVeeumlKut79erFqVOnGhqGRbgH\nBaErKeHK+fM4u7jQxt+fi/v3V2nn3KIFoatW0dzdnZvXr3P57Fm+eOIJfty82QpR18zYnH7LzWXt\nX/7C4HfeYdqRI1zPy0O7bRs7Fy60QtS1MzYvgBYeHviEhrL16actHGXdGJvT92++iaIo/OnFF2n5\nwQfcyM8ndfduvnnxRStEXbugIHdKSnScP38FFxdn/P3bsH//xSrtnJ2deO21P+Pt3ZrSUh379qXx\n4INrSEn51TyB1ePSQmtRKbVdnW4DVCoVkdYOwsQW3/61L7Gxr7oN0RhzgsaZV0VOKpWJxpNthKIs\nRqVS1XrTTW1UKhWsMrKP5xr+fg0lk2YJIYQh/0tDK0II0ShJIRdCCDtnhssPzcVuxsjtIEwhhA0w\n2Rj5UiP7eNn69UmOyIUQwhA7umrFbgp5Y7piABr3lRCNKSdonHk1xpzg97xMQsbIhRDCztnRGLkU\nciGEMMQMt+ibixRyIYQwRIZWhBDCzkkhF0IIO2dHY+QNemanEEI0WiVGLnX0/vvv4+vrS7du3VhY\nw6RzOp2OoKAgg08R+iM5IhdCCEPMMLSyZ88eEhIS+OGHH1Cr1fz6a/UzN1Y83L7w9vz/NZEjciGE\nMOSmkUsdfPDBB7z44ouo1WoA2rRpY7BdxcPtp06datRdo1LIhRDCEJ2RSx1otVr2799Pv379GDBg\nAEePHjXYruLh9g4OxpVoGVoRQghDqhtaubwXruytdreQkBBycnKqrF+6dCllZWVcvXqVQ4cOceTI\nEcLDw7nwh6c2GfNw+z+SQi6EEIZUV8hbDShfKpyv/HCOnTt3VtvlBx98wKhRowDo06cPDg4OXLly\nBVdXV30bQw+3j4iIqPZ5nyBDK0IIYZgZxshHjBjB7t27ATh//jylpaWVijiUP9w+PT2d1NRU4uLi\n+Mtf/lJjEQcp5EIIYZgZLj+cPHkyFy5cICAggLFjx+oLdFZWFkOHDjW4j8qIic1kaEUIIQwxw+WH\narWadevWVVnv4eHBtm3bqqy/8+H2NZFCLoQQhtjRnZ1SyIUQwhCZ/VAIIeycHU2aVa+Tnenp6Tzy\nyCP4+/vTrVs3VqxYAUBeXh4hISF06dKFQYMGkZ+fr98nKioKHx8funbtytdff61ff+zYMQICAvDx\n8WH27NkNTEcIIUykzMjFBtSrkKvVat59913OnDnDoUOHWLVqFWfPniU6OpqQkBDOnz/PwIEDiY6O\nBiAlJYWNGzeSkpJCYmIiM2bM0N92On36dGJjY9FqtWi1WhITE02XnRBC1JcZLj80l3oVcnd3d3r0\n6AFA8+bN8fX1JTMzk4SEBCZMmADAhAkT2Lx5MwDx8fGMHTsWtVqNl5cXnTt3JikpiezsbAoLCwkO\nDgYgIiJCv48QQliVmWY/NIcGj5GnpaVx/Phx+vbtS25uLm5ubgC4ubmRm5sLlF8j2a9fP/0+Go2G\nzMxM1Go1Go1Gv97T05PMzMyGhiSEEA1nI8MmxmhQIS8qKmL06NEsX76cFi1aVNqmUqmMupDdWHvu\n+NkL6GSynoUQ9iwVSLv9sxIZabqObWTYxBj1vrPz5s2bjB49mvHjxzNixAig/Ci8YrKY7Oxs2rZt\nC5Qfaaenp+v3zcjIQKPR4OnpSUZGRqX1np6eBt/vkTsWKeJCiAqd+L02RJqykJth9kNzqVchVxSF\nKVOm4Ofnx5w5c/Trhw8fztq1awFYu3atvsAPHz6cuLg4SktLSU1NRavVEhwcjLu7Oy4uLiQlJaEo\nCuvWrdPvI4QQVmVHV63Ua2jl+++/57PPPqN79+4EBQUB5ZcXLlq0iPDwcGJjY/Hy8mLTpk0A+Pn5\nER4ejp+fH05OTsTExOiHXWJiYpg4cSLXr18nNDSUIUOGmCg1IYRoABsp0sZQKcY8fsLKVCoVkdYO\nwsQW3/61LzHheQRra4w5QePMqzHmBOV5qVQqo56qUxOVSgVORvZR1vD3ayi5s1MIIQyxoyNymcZW\nCCEsZMyYMQQFBREUFESnTp30Q9N3qu7O+ZrIEbkQQlhIXFyc/ufnn3+eVq1aVWlTced8jx49KCoq\nolevXoSEhODr61ttv1LIhRDCwhRFYdOmTezZs6fKNnd3d9zd3YHf75zPysqSQi6EEHVnvjuCvv32\nW9zc3PD29q6x3Z13ztdECrkQQhhU3dnO/bcXw0JCQvQ3Rt5p2bJlDBs2DIANGzYwbty4Gt+9qKiI\nsLAwli9fTvPmzWtsK4VcCCEMqu6I/IHbS4Vllbbu3Lmzxl7Lysr46quvSE5Orv6db985/7e//c2o\nmySlkAshhEHXzdLrrl278PX1xcPDw+D26u6cr4lcfiiEEAaZZ0LyjRs3Mnbs2ErrsrKyGDp0KPD7\nnfN79uzRX6pY23Ma5IhcCCEMMs8dQR9//HGVdR4eHmzbtg2AP/3pT9y6datOfUohF0IIg+xnHlsp\n5EIIYZD93KMvhVwIIQyynyNyOdlZD/NzcmjXsycAE/fto9uYMZW2ewYHM/n773mpuJh5mZn8ZelS\nsINZ5mrKq42fH2GbNjHz3DleLStj2OrV1gqzzmrKq8ekSUTs3s3zv/zCooICph05Qrc/nIiyRTXl\n5D1oEJMPHOD5X37hpeJi/q7V8sjrr+PgZPvHbbX9v1XhXl9fXiwq4pXSUjNGc93Ixfps/5O1Mfd4\ne6Nu1ozs48dxUKvx6N2bS999p9/uotEwfudOUr74goQpU3Dt0oXha9agUqn45qWXrBh5zWrLy6lp\nUwrS0jgXH88D8+ZZfdpOY9WWl9cjj/DjV1+x8/nnuZ6XR9eRIxn56afcKisj5YsvrBh59WrL6UZB\nAYfefZdfTp+mtLCQdj178tfVq2nSogX/nTvXipHXrLa8Kjg1bcoTmzaR+s03dDbr8wtkaKXR6vDQ\nQ2QmJYGi4NmnD8VXrnDtjsfV9Z4+nRv5+SRMnQrA5R9/ZM+rrxLy5pvse/11ym7csFboNaotr+xj\nx8g+dgyAoClTrBVmndWW1+aIiErtD737Lh3798c/PNxmC3ltOWUmJZVvv+1aRgZeAwbQsX9/a4Rr\ntNryqhC6ahUX9+8nMymJzo89ZsaI7GdoRQq5kRZevYqiKDg5O6NycGBBXh6OajWOzs4syMsDReFN\nV1faP/QQP3/9daV9f/7vfwlduRL3oCAyDh60UgaGGZuXvWlIXk3vuYerFy5YOOLa1Tcn1/vvx3vI\nEH78z3+sEHXt6pJX9/Hj8ejVi3/16WOBITA5Im90PujeHZVKxZRDh9j27LPknDjB6Lg4Tn/+OT/G\nx+vbNXd359K331bat+j2vAst2rWzaMzGMDYve1PfvAKeegrPvn3ZMWuWBaM1Tl1zmpueTrN778Wx\nSRNOfPwxu195xQpR187YvO7t2pVB//wnnwwYgM6sY+MV7OeIXE52GulaejrOLVviqFZzbssWrl+9\ninuPHpyOi+NaejrX0tOtHWK9SF6/u3/4cIatXk3C5MnknjxphahrVtec1jz0EB8FBfHV+PF4Dx7M\nkOXLrRR5zYzJy7FJE5744gt2v/IKl8+etVBk9vP0ZTkiN8L006dp2aEDDk5OOKrVLCooQOXggJOz\nM7NufwVf5etLYWYmRdnZVY6873ZzA6AwO9visdekLnnZk/rk5f/kkzz+8cdsmTqVU59/bq3Qq1Wf\nnAouXQLKz9Pc0ukYtX4937z4IjeLi62SgyHG5uXg5EQbPz9CV60idNUqoPy5mioHB14pLWXPq6/y\n/RtvmDg6+zkil0JuhPVDhuDYpAnD16zhpx07OLNpE/0XL0ZXUsJ30dEAFN0u0unff0/38eMr7d95\nyBBKf/uNnOPHLR57TeqSlz2pa149p05lyIoVbI6IIOXLL60Vdo0a+lk5ODoCoLr9X1thdF4qFTHd\nulXat+uIEQxYsoQPAwP57ZdfzBCdbVxaaAwp5Ea4lpGBysEBt+7d2fr00+SnpuIWEMDeyEjyU1Mr\ntT3ywQf0mTmTYf/6F4fefZd7vL155PXXOfz++zZ3xUpd8nJwcqKNvz8Azi1a0NTVFbfAQHSlpRb8\nqmucuuTVb84cHn3zTbY/9xwXv/1W/+1JV1rKjatXrRG+QXXJ6YF58/j17FnytFoURcGjd28efeMN\nfty8mdLCQitlYFhd8vrj31lhcLDB9aYjR+SNjntQELqSEq6cP4+ziwtt/P25uL/q5PKFmZl8NmgQ\ng955h2lHj3IjP59jH31ksyeajM2rhacnz9yeP1lRFNr17InvyJHkp6WxopannFiDsXkFz5qFysGB\nv374IX/98EP9+rS9e/l04EBLhlwrY3NycHIi5M03aeXlhXLrFvlpaRxeuZJD771nhahrZ2xeBpn1\nfgbbGP82hkqxgTs7EhMTmTNnDjqdjqlTp7Jw4cJK21UqFZHWCc1sFt/+tS+xgzs+jdUYc4LGmVdj\nzAnK81KpVA2+YU2lUgExRraeYfT7HT58mJkzZ3Lz5k2cnJyIiYmhT58+VdrVVhP/yOpXreh0OmbO\nnEliYiIpKSls2LCBszb2Vb1Cau1NzM4WYgDbiMMWYgDbiMMWYgDbicM0TH/VyoIFC/jHP/7B8ePH\nef3111mwYEGVNvWpiVYv5IcPH6Zz5854eXmhVqsZM2YM8TZ6/XKatQPANmIA24gjzdoB3JZm7QCw\njRjAduIwDdM/WKJdu3YUFBQAkJ+fj6enZ5U29amJVh8jz8zMpH379vrXGo2GpDtuLxZCCOsw/Rh5\ndHQ0f/rTn3j++ee5desWBw3c6V2fmmj1Qq4ycoxusfWH8lEiI1kcGWnSPuualzliqI+a4rDUZ2Xp\n30V1ednCZ1LfGEz9WdnC78J06nf5YUhICDm37+a+09KlS1mxYgUrVqxg5MiRfPHFF0yePLnKw5qN\nrYmVKFZ28OBBZfDgwfrXy5YtU6Kjoyu18fb2VgBZZJFFlloXb2/vBtelurxf8+bNje63RYsW+p9v\n3bqluLi4VGljTE38I6sfkffu3RutVktaWhoeHh5s3LiRDRs2VGrz008/WSk6IcT/IsVM3yo7d+7M\nvn376N+/P7t376ZLly5V2hhTE//I6oXcycmJlStXMnjwYHQ6HVOmTMHX19faYQkhhMmtXr2a5557\njpKSEpo2bcrq2w9oycrKYtq0aWzbtq1eNdEmriMXQghRf1a//LA2iYmJdO3aFR8fH94w8aQ46enp\nPPLII/j7+9OtWzdWrFgBQF5eHiEhIXTp0oVBgwaRn5+v3ycqKgofHx+6du3K13fMO37s2DECAgLw\n8fFh9uzZdY5Fp9MRFBTEsGHDrBZDfn4+YWFh+Pr64ufnR1JSksXjiIqKwt/fn4CAAMaNG0dJSYlF\nYpg8eTJubm4EBATo15nyfUtKSnjyySfx8fGhX79+XLx40agYXnjhBXx9fQkMDGTUqFH6S9fMFUN1\ncVR4++23cXBwIC8vz+K/C4D3338fX19funXrVukmGXP9LuyG0aP0VlBWVqZ4e3srqampSmlpqRIY\nGK70m/sAAAYoSURBVKikpKSYrP/s7Gzl+PHjiqIoSmFhodKlSxclJSVFeeGFF5Q33nhDURRFiY6O\nVhYuXKgoiqKcOXNGCQwMVEpLS5XU1FTF29tbuXXrlqIoitKnTx8lKSlJURRFeeyxx5QdO3bUKZa3\n335bGTdunDJs2DBFURSrxBAREaHExsYqiqIoN2/eVPLz8y0aR2pqqtKpUyflxo0biqIoSnh4uPLJ\nJ59YJIb9+/crycnJSrdu3fTrTPm+q1atUqZPn64oiqLExcUpTz75pFExfP3114pOp1MURVEWLlxo\n9hiqi0NRFOXSpUvK4MGDFS8vL+XKlSsW/13s3r1befTRR5XS0lJFURTll19+Mfvvwl7YdCE/cOBA\npbO3UVFRSlRUlNne7/HHH1d27typ3H///UpOTo6iKOXF/v7771cUperZ48GDBysHDx5UsrKylK5d\nu+rXb9iwQXnmmWeMft/09HRl4MCByu7du5W//vWviqIoFo8hPz9f6dSpU5X1lozjypUrSpcuXZS8\nvDzl5s2byl//+lfl66+/tlgMqamplQqHKd938ODByqFDhxRFKf9H8t577zUqhjv95z//UZ566imz\nx1BdHGFhYcrJkycrFXJL/i6eeOIJ5ZtvvqnSzty/C3tg00Mrhi6MzzTT3NhpaWkcP36cvn37kpub\ni9vtWfDc3NzIzc0Fyk9IaDSaKvH8cb2np2ed4pw7dy5vvfUWDg6/fxyWjiE1NZU2bdowadIkevbs\nybRp0/jtt98sGkfr1q2ZP38+HTp0wMPDg1atWhESEmLx30UFU77vnX/LTk5OtGzZstLwhDHWrFlD\naGioVWKIj49Ho9HQvXv3SustGYdWq2X//v3069ePAQMGcPToUav8LmyRTRfyel0YXw9FRUWMHj2a\n5cuX06JFiyoxmDOOrVu30rZtW4KCgqq95MncMQCUlZWRnJzMjBkzSE5O5u677yb69nzQlorj559/\n5r333iMtLY2srCyKior47LPPLBpDdaz1vhWWLl1KkyZNGDdunMXfu7i4mGXLlrFkyRL9uur+Vs2p\nrKyMq1evcujQId566y3Cw8MtHoOtsulC7unpSfodj69KT0+v9C+sKdy8eZPRo0czfvx4RowYAZQf\nfVXcmZWdnU3btm0NxpORkYFGo8HT05OMO572nZGRYXAOBUMOHDhAQkICnTp1YuzYsezevZvx48db\nNAYoP4rRaDT6mdjCwsJITk7G3d3dYnEcPXqUBx98EFdXV5ycnBg1ahQHDx60aAx3MsVnUPH36unp\nyaXbT+wpKyujoKCA1q1bGxXHJ598wvbt21m/fr1+nSVj+Pnnn0lLSyMwMJBOnTqRkZFBr169yM3N\ntWgcGo2GUaNGAdCnTx8cHBy4fPmyxT8Pm2TtsZ2a3Lx5U7nvvvuU1NRUpaSkxOQnO2/duqWMHz9e\nmTNnTqX1L7zwgn7MLSoqqsoJppKSEuXChQvKfffdpz+pEhwcrBw6dEi5detWvU40Koqi7N27Vz9G\nbo0YHn74YeXcuXOKoijK4sWLlRdeeMGicZw4cULx9/dXiouLlVu3bikRERHKypUrLRbDH8dkTfm+\nq1atUp599llFUcrHaqs7ufbHGHbs2KH4+fkpv/76a6V25ozBUBx3MnSy0xK/iw8//FB57bXXFEVR\nlHPnzint27e3yO/CHth0IVcURdm+fbvSpUsXxdvbW1m2bJlJ+/72228VlUqlBAYGKj169FB69Oih\n7Nix4/+3b8cmFgJhEIAj7UEQQVYRZDEQwQLEAgwEA7EDyzCwCK3DTAswsAMrEA3FZC46MXtw3Ht3\nC/OFguzwIxMsv9i2DUmSwHVdpGmKfd/vd5qmgRACnudhGIb7+TzPkFJCCIG6rn+UZ5qme2vlLzIs\ny4IoihAEAbIsw3EcH8/Rti1834eUElVV4bquj2QoigKGYUDTNJimib7vf/Xc8zyR5zkcx0Ecx1jX\n9WWGruvgOA4sy7q/z+9Ni3dleObQdf2exZNt23eRv3sWzwzXdaEsS0gpEYYhxnF8+yxUwR+CiIgU\n96/vyImI6DUWORGR4ljkRESKY5ETESmORU5EpDgWORGR4ljkRESKY5ETESnuC6fXQ7k8AS+DAAAA\nAElFTkSuQmCC\n", + "text": [ + "" + ] + } + ], + "prompt_number": 34 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "Msig = M.getEdgeInnerProduct(sig)\n", + "MsigBG = M.getEdgeInnerProduct(sigBG)\n", + "\n", + "Mmu = M.getFaceInnerProduct(mu_0, invProp=True)" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 38 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "freq = 1e1\n", + "C = M.edgeCurl\n", + "A = C.T*Mmu*C - 1j*omega(freq)*Msig\n", + "ABG = C.T*Mmu*C - 1j*omega(freq)*MsigBG" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 39 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "def get1Dfields(sigHalf, freq):\n", + " hz = [(100.,18)]\n", + " M = Mesh.TensorMesh([hz],'C')\n", + " sig = np.zeros(M.nC) + 1e-8\n", + " sig[M.vectorCCx<=0] = sigHalf\n", + " G = M1D.nodalGrad\n", + " Mmu = Utils.sdiag(M.vol*(1.0/mu_0))\n", + " Msig = M.getFaceInnerProduct(sig)\n", + " A = G.T*Mmu*G - 1j*omega(freq)*Msig\n", + " Aii = A[1:-1,1:-1]\n", + " Aio = A[1:-1,[0,-1]]\n", + " bc = np.r_[0.0,1.0]\n", + " rhs = -Aio*bc\n", + " eii = Solver(Aii).solve(rhs)\n", + " e = np.r_[bc[0],eii,bc[1]]\n", + " return e" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 89 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "x0 = -h1.sum()/2" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 87 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "sigs = np.logspace(-3,1,4)\n", + "for s in sigs:\n", + " e = get1Dfields(s,1e1)\n", + " plot(e.imag,M.vectorNx)\n", + "ylim([-900,0])\n", + "ylabel('Depth, m')\n", + "legend(['$\\sigma = 1e%d$'%log10(s) for s in sigs],'best')" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "metadata": {}, + "output_type": "pyout", + "prompt_number": 101, + "text": [ + "" + ] + }, + { + "metadata": {}, + "output_type": "display_data", + "png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEACAYAAACUMoD1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XlYVOX7x/H3KLiL+1IzKSYogqikorZShrmkqZWGFVJZ\n5l6m2e5Siqa2mJFWal/J3EuslMCFckUFd1xQsWAAV1REZD2/P47wcwEZYGbODNyv6+JKD2fO3Egz\nn3mW8zw6RVEUhBBCiFKqoHUBQgghygYJFCGEEGYhgSKEEMIsJFCEEEKYhQSKEEIIs5BAEUIIYRZ2\nFyihoaG4ubnh6urKjBkztC5HCCHEDTp7ug8lJyeHli1bsmHDBvR6PR07dmTp0qW0atVK69KEEKLc\ns6sWyq5du3BxccHZ2RlHR0deeOEFQkJCtC5LCCEEdhYoRqOR++67L//vBoMBo9GoYUVCCCHy2FWg\n6HQ6rUsQQghRCAetCygOvV5PfHx8/t/j4+MxGAy3nOPi4sLJkyetXZoQQti15s2bc+LEiVJdw65a\nKB06dCA2NpbTp0+TmZnJ8uXL6dOnzy3nnDx5EkVR7PZr4sSJxTrfyUkhJUX7uu9We25uLuEXLvD0\ngQPU37qV906e5L/0dM3rLe2/va19Sf12UH92NtmvvsLe+6vzw6ZZmtd885c5PojbVQvFwcGBuXPn\n8tRTT5GTk8Nrr71W7md4VasG165B7dpaV3Knazk5BJ85w5yEBCrodIzW61nu7k61ihW1Lk0I68vO\nhldeoWJCAjU2beX9Fb50aNWVdo3baV2Z2dhVoAD06NGDHj16aF2GzcgLFFvy7/XrBBmNLExO5kEn\nJ75xdeXx2rVlDEyUX5mZ8OKLkJoKf/6JS7VqfPXUV7yw6gX2vLGHGpVqaF2hWdhVl1d54OPjU6zz\nbSVQFEWhbocOPHfoEA/s2UOWohD5wAOEeHryRJ06dhEmxf23tzVSv7YKrf/6dXjuOTVUQkLUFy3w\nYpsX6XJfF0auG2m9Ii3Mrm5sNIVOp6OM/Uh35e0N33wDnTpp8/zXc3JYdvYsc4xG0nJyGKXXM7hx\nY2o62F3jVwjzu3YN+vWDWrVgyRJwdLzl22mZabT/vj0fPvIhL7d9WaMiVeZ475RXvZ3TqoWSmJHB\nd4mJfJ+YiFeNGkxt1oyn6talgh20RISwitRU6N0bmjSBhQuhgA9Z1StVZ8XzK+i6uCudDJ1oUa+F\nBoWaj3R52TlrB8quK1d4MSYGj927uZiVxd/t2hHati096tWTMBEiz+XL8NRT0KIF/PRTgWGSp02j\nNkzxmcLAVQPJyM6wXo0WIIFi56wRKFm5uSw9c4bOUVEMjImhfc2axHXqxLctWuBWvbpln1wIe3Ph\nAnTtCh06wPz5UKHot9k3O7xJ8zrNGR8+3goFWo50edk5SwbKucxM5icm8l1iIi2qVeO9Jk3oXb8+\nFaUlIkTBzpwBX1/o0QOmTwcTXys6nY4fev+A13wvujbryjNuz1i4UMuQQLFz1apBerp5r7kvNZU5\nRiO/nT/Ps/Xrs75NG9rUKBvTGoWwGKMRnnwSBg6EiRNNDpM8darWYemzS+m7vC9e93jRpFYTCxVq\nORIodq5qVfO0ULJzc1l74QJfJyRwMj2dEXo9sd7e1K9UqfQXF6Ks+/dftZvr9ddhwoQSX6bLfV0Y\n23ksg1YPIiIgAocK9vUWLWModq60XV4pWVnM/O8/XCIjmR0fz3C9nrjOnXm/aVMJEyFMceIEPPYY\njBpVqjDJM/6h8VRzrMakiEmlr83K7Cv+xB2qVYO0tOI/LiYtjW+MRpadPcvT9eqx0sODjk5O5i9Q\niLLs6FF1zOTjj+GNN8xyyQq6CgT3C8ZrvhePOz9O1/u7muW61iAtFDtXnBZKrqLw54ULdNu/nyf2\n7aORoyMxHTsS3KqVhIkQxXXgADzxBEydarYwydOoRiMW91uM/xp/zqadNeu1LUlaKHbOlEH5K9nZ\n/JSczDdGI04VKzLGYGCgpyeVTZjOKIQoQFQU9OoFX3+tDsJbwJP3P8ngtoPx/82fdS+uo4LO9l+v\ntl+huKu7tVBOXLvGmNhYnHfuZOvly/zk5sae9u3xb9xYwkSIktqxQ50WPG+excIkz2SfyaRmpjJ7\n+2yLPo+5SAvFzt0+y0tRFDakpDDHaGTnlSsMuece9nfowH1VqmhXpBBlRUQEDBgAixdD9+4WfzrH\nio780v8XvH/05pGmj9DZ0Nniz1kaEih2Lq+FkpaTw8+y94gQlhMWpi5Bv3y5OnZiJU1rN2Ver3n4\nrfZj79C91K5ig5sf3SCrDdu5ZRHXGb/VSLpPEg/VqsUYg0H2HhHC3H7/HV57DX77DR56SJMSRq4b\nyZm0M6x4boVFXt/meO/UpCN95cqVeHh4ULFiRaKjo2/5XmBgIK6urri5uREWFpZ/PCoqCk9PT1xd\nXRkzZoy1S7YpiqLwz6VLPHvoEEPZQ1auwq727e1q7xEh7MbKlTBkCPz5p2ZhAjCr2yxiL8TyfdT3\nmtVQJEUDR44cUY4dO6b4+PgoUVFR+ccPHz6stG3bVsnMzFTi4uKU5s2bK7m5uYqiKErHjh2VyMhI\nRVEUpUePHsr69esLvLZGP5JVpGdnK4sSE5V2u3crLXbuVL6Jj1d27s9SPDy0rkyIMio4WFEaN1aU\nffu0rkRRFEU5eu6oUv/z+sqB5ANmv7Y53js1aaG4ubnRosWd6/6HhITg5+eHo6Mjzs7OuLi4EBkZ\nSVJSEqmpqXh7ewPg7+/PmjVrrF22ZhIzMvg4Lo6mO3ey7OxZpjVrxhFvb0YaDDSo4WATOzYKUeb8\n+KN65/vGjdC2rdbVANCyfktmd5vNwFUDScsswR3NFmZTc0cTExMxGAz5fzcYDBiNxjuO6/V6jEaj\nFiVa1a4rVxhUxN4j5lrLSwhxk7lzYcoUdVaXu7vW1dzCv60/He7twOj1o7Uu5Q4Wm+Xl6+tLcnLy\nHcenTZtG7969LfW0AEyaNCn/zz4+Pna5V/XRtDS6HzjAh02bEuTqSu3btg7NU6WK+VcbFqJcungR\nli6FBQvUT2l//w3NmmldVYGCegXhMseFExdP4FLXpUTXiIiIICIiwqx1WSxQwsPDi/0YvV5PfHx8\n/t8TEhIwGAzo9XoSEhJuOa7X6wu9zs2BYq+irl7Ft04d3rnvvrueV6EClKNJbUKYV26u2qW1cCGs\nX6/esDhjhjot2Ian3NeoVINmdZpxNu1siQPl9g/bkydPLnVdmnd5KTe9G/bp04dly5aRmZlJXFwc\nsbGxeHt707hxY5ycnIiMjERRFIKDg+nbt6+GVVvegatX8TRhN0QJFCFKIC5O3bOkWTN47z14+GE4\ndUptofj62nSY5KlTpQ4p6Slal3ELTQLlt99+47777mPnzp306tWLHj16AODu7s6AAQNwd3enR48e\nBAUF5U+BDQoKYsiQIbi6uuLi4kJ3K9ylqqWDaWl4mrCplU6nfsgSQhQhPR2WLFH3LenYEVJSICRE\nXZdrxAioW1frCoulTtU6pFy3rUCRGxtt1H07dvB3u3bcX7XqXc9LS4MGDWRgXogCKYoaGAsXqne4\nd+wIr74KffqoA5B2bOS6kbSs15JRnUaZ5XrmeO+UpVdsUEpWFpeys3E24X94nU66vIS4w7lzamtk\n4UK4elUNkX37oIgxSXtSp4rttVAkUGzQobQ0PKpVy58afDcSKELckJ2trre1cCFs2KC2QubMgUcf\nVQcby5g6VesQfzm+6BOtSALFBpk6fgISKEJw4oQaIv/7HxgMamtkwQKoVUvryiyqTpU6HDhzQOsy\nbiGBYoMOpqWZNMMLZJaXKKfS0mDVKjVIjh6Fl16Cv/6C1q21rsxqbHFQXgLFBh28epUBDRqYdK7M\n8hLlhqLAzp1qiKxapU71festdefESpW0rs7qbHHasASKjVEUhUPFaKFIl5co85KTIThYDZKcHLVL\n6/BhuPderSvTlLRQRJHiMzKoWrEi9U38xCWBIsqkrCz1zvWFC9UlUPr1gx9+UJePl+0ZAGmhCBMU\nZ/wEJFBEGXPkCCxapG6x6+KitkaCg6FmTa0rsznSQhFFOmjikit55MOasHtXrsCKFWprJC4OBg9W\nWyUtW2pdmU2r7lidrJwsMnMyqVTRNsaQJFBszMG0NHzr1DH5/LxAURQJF2FHFAW2blWn965ZA48/\nDu+/ry7O6CBvS6bQ6XTUrlKblPQUGtVopHU5gA0sDiluVZx7UPLITC9hN4xGCAyEFi1g6FDw9IRj\nx9S92nv3ljApJlvr9pLfng3Jys0lNj0d92rVivU4GUcRNi0zE37/Xe3S2rEDnn8efv4ZvL2lWV1K\ndarU4WL6Ra3LyCeBYkOOXbtGk8qVqVrMpbMlUIRNOnRIDZGff1Z3PXztNXWspBhjhOLu6lS1rZle\nEig2pLgzvPJIoAibcekSLFumBkliIgQEwPbt6owtYXa2tkCkBIoNKcn4CUigCI3l5qp7ry9cCH/8\nAd26qfux28lGVfbM1u5FkUCxIQfT0nilceNiP07W8xKa+O8/dUHGRYvU+0Reew2++grq19e6snLD\n1gblNZnlNX78eFq1akXbtm3p378/ly9fzv9eYGAgrq6uuLm5ERYWln88KioKT09PXF1dGTNmjBZl\nW1xx70HJI7O8hNVcv65uVPXUU+DlpS6LsnKlutfI6NESJlZmay0UTQKlW7duHD58mP3799OiRQsC\nAwMBiImJYfny5cTExBAaGsrw4cPzdxAbNmwYCxYsIDY2ltjYWEJDQ7Uo3WKuZGdzLiuryB0aCyJd\nXsLi9u6FUaPU5eF//FEdG0lIgG+/hfbtZbaWRqSFAvj6+lLhxoY3nTp1IiEhAYCQkBD8/PxwdHTE\n2dkZFxcXIiMjSUpKIjU1FW9vbwD8/f1Zs2aNFqVbzKG0NNyrV6diCV6YEijCIi5cgG++UVsi/fqp\nrY89eyA8HPz8oAQffoR5yaD8bRYuXIifnx8AiYmJdO7cOf97BoMBo9GIo6MjBoMh/7her8doNFq9\nVksq6QwvkEARZpSTo+52uGgRhIaqS8PPmqXeyV4Gdz20d+Vm2rCvry/Jycl3HJ82bRq9e/cGYOrU\nqVSqVIlBgwaZ9bknTZqU/2cfHx98fHzMen1LMGZk0KRy5RI9VpZdEaV2/Lg6wL54MTRqpC7K+N13\nUIxlgIT1HTp7iPrVSjZuFRERQUREhFnrsVighIeH3/X7P/30E+vWrWPjxo35x/R6PfHx/79HckJC\nAgaDAb1en98tlndcr9cXeu2bA8VeOOp0ZJRwZD0zE0qYRaI8u3JFHVBftAhiY9VdD9etU5dDETYv\nPSudwK2BrH1hbYkef/uH7cmTJ5e6Jk3asKGhocycOZOQkBCqVKmSf7xPnz4sW7aMzMxM4uLiiI2N\nxdvbm8aNG+Pk5ERkZCSKohAcHEzfvn21KN1iKlWoQGYJ+q1ycyE7W5ZAEibKzYXNm8HfH5o0Ue8b\nefdddYB99mwJEzsyP2o+He/tSPt722tdSj5N3oZGjRpFZmYmvr6+AHTp0oWgoCDc3d0ZMGAA7u7u\nODg4EBQUhO5GX05QUBABAQGkp6fTs2dPunfvrkXpFlNJpyOzBC2UzEx191Pp8hJ3FRendmn973/g\n5ASvvKIGiIlbTQvbkpaZxoxtM/jrpb+0LuUWOkUpW8O5Op0Oe/yRgoxGDqWlEdSiRbEed+WKOpPz\nyhULFSbsV1oarF6tdmkdOqTOzHrlFWjXTj6B2LnPt33OnsQ9rHh+hdmuaY73TukosRElbaFkZKgt\nFCEAdYbGtm1qiPz6q7pl7siR8PTTMtBWRqRmpDJr+yw2D96sdSl3kECxESUdQ8nr8hLlXHy8OkPr\np5/A0VFticTEwD33aF2ZMLM5kXPwbe6LR0MPrUu5gwSKjahUwlleMsOrHEtPV3c7XLQIoqJgwADZ\nZ6SMu3T9El9FfsW2V7dpXUqBJFBsRElbKNLlVc4oCuzapYbIypXQoYPaGgkJkTvXy4Gvdn5FL9de\ntKhXvLFWa5FAsRGlneUlyrikJAgOVru0srLUtbT27YP77tO6MmElF9MvMnfXXHa9vkvrUgolgWIj\nSjOGIl1eZVRGhrp17qJF6iZV/fvD99+rA+3SpVXuzNo+i/6t+nN/nfu1LqVQEig2QmZ5CUDt0tq7\nVw2RZcugdWu1S0u2zi3XzqWdY37UfKLfiNa6lLuSQLER0kIp586ehSVL1C6ty5fVLq1du6BZM60r\nEzbg822f84LHCzSt3VTrUu5KAsVGSAulHMrKUtfOWrRI3UK3Tx91x8PHHpOVfUW+5KvJLNi7gIPD\nDmpdSpEkUGyE3IdSjhw8qIbIkiXg6qp2aS1erC6JIsRtpm+djn9bf/ROhS+IayskUGxE5QoVSjzL\nS7q87MDFi/DLL2qXVnIyDB4MW7eqgSJEIRKuJBB8IJjDww9rXYpJJFBsRCWdTu5DKWuysyEsTA2R\nsDDo0QOmTYOuXaFiRa2rE3Zg2pZpvOb1Go1rNNa6FJNIoNiISqVooUig2JijR9UQCQ5WV+585RV1\num/t2lpXJuzIv5f+Zfnh5RwdcVTrUkwmgWIjStpCkS4vG3H5Mixfro6NnD4NL7+s7r3u7q51ZcJO\nffbPZ7zZ/k0aVLefLQYkUGxESVso0uWloZwc2LRJbY38+Sc8+SR89BE89ZTseCZK5eTFk/x29DeO\njzqudSnFIv/X2whpodiREyf+f7Oq+vXVLq2vv1b/LIQZTPlnCiO9R1K3al2tSykWTSa7f/zxx7Rt\n25Z27drRtWvXW/aRDwwMxNXVFTc3N8LCwvKPR0VF4enpiaurK2PGjNGibIuSFoqNS01Vu7MefRQe\nfBCuXlWXRYmOhlGjJEyE2Rw7f4x1set4u/PbWpdSbJoEyrvvvsv+/fvZt28fffv2ZfLkyQDExMSw\nfPlyYmJiCA0NZfjw4fk7iA0bNowFCxYQGxtLbGwsoaGhWpRuMRVvrM2UU8xWigzKW1BurnrDYUCA\nugjjmjXw9tvq/utffglt22pdoSiDJv89mbc7v02tKrW0LqXYNAmUmjVr5v/56tWr1L/x6S4kJAQ/\nPz8cHR1xdnbGxcWFyMhIkpKSSE1NxdvbGwB/f3/WrFmjRekWVZJWinR5WcDp0zBlCri4qLsdenrC\nsWPqEvH9+kmCC4s5dPYQG+M2Msp7lNallIhmYygffvghwcHBVK1alV271OWYExMT6dy5c/45BoMB\no9GIo6MjBoMh/7her8doNFq9ZkvL22SrajHuUZAuLzPJzYXffoOgIHVZ+BdeUBdkbN9eVvYVVqEo\nCh9t+ohxXcZRs3LNoh9ggywWKL6+viQnJ99xfNq0afTu3ZupU6cydepUpk+fzltvvcWiRYvM9tyT\nJk3K/7OPjw8+Pj5mu7YlOeh0ZBezy+vaNdlXqVRyctSNqj77DKpVg3fegWeegSpVtK5MlCOKovBW\n6FskXU1ihPcIqzxnREQEERERZr2mxQIlPDzcpPMGDRpEz549AbXlcfMAfUJCAgaDAb1eT0JCwi3H\n9frC17W5OVDsiQ4o7jyv1FSoaZ8fZrSVna0uD//ZZ1CnDsyapU73ldaIsLK8MNlp3MlfL/1FNcdq\nVnne2z9s541ll4YmYyixsbH5fw4JCcHLywuAPn36sGzZMjIzM4mLiyM2NhZvb28aN26Mk5MTkZGR\nKIpCcHAwffv21aJ0i9LpdMUOlKtXJVCKJStLvW+kVSuYPx/mzlU3r+reXcJEWJ2iKLz919vsSNjB\nXy/9Re0q9r2agiZjKO+//z7Hjh2jYsWKNG/enO+++w4Ad3d3BgwYgLu7Ow4ODgQFBaG78SIPCgoi\nICCA9PR0evbsSffu3bUo3aKkhWJBmZnqUijTpkHTpvDDD+oy8RIiQiN5YbI9fjthL4fZfZgA6BSl\nBHfT2TCdToe9/kiNtm1jf4cONC7GtK2OHeHbb+HGBDhxu4wMtUUSGAgtWsDHH8Mjj2hdlSjnbDFM\nzPHeKXfK2xBpoZjR9evw448wY4Y67XfpUujSReuqhEBRFMb+NZZt8dsIfzncJsLEXCRQbEhJxlAk\nUG5z7Zq6su/MmeqU319/VZtxQtiAvDDZGr+1zIUJSKDYFGmhlEJaGnz3HcyerbZEfv8dHnhA66qE\nyKcoCu+EvcPW+K2EvWQb3VzmJoFiQ3RQrD5MRVHfR2vUsFxNNi81VR1E+vJLdZD9r7+gTRutqxLi\nFnlh8s+//xD+cjh1qtbRuiSLkECxIcVtoaSlqffflcvN/y5fhm++UVf59fVVl5H38NC6KiHuoCgK\n48LGlfkwAQkUm1LcMZRy2d2VkqKGyNy50LMnbNkCbm5aVyVEgfLC5O9//y7zYQIa3dgoClbcFkq5\nCpQLF9Qpv66u8N9/sHMnLF4sYSJslqIojA8fT8S/EeUiTEACxaYUdwwlNbUcjJ+cOwfvv6/eQ3Lm\nDOzaBQsXqisBC2Gj8sJk8+nNbHh5Q7kIE5BAsSnFbaGU6WVXzpyB8eOhZUt1vCQ6Wp0OfP/9Wlcm\nxF3dHCblpWWSRwLFhsgYCpCYqG5i1aqVepf7gQPqkvJNm2pdmRBFUhSFd8PfZVPcJsJfDre7LXxL\nSwLFhpTrMZSEBHUr3dat1fW1Dh+GOXPgpn1whLBleWGyMW4jG/w3lLswARNmeWVnZ/Pnn39y+vRp\nsrOzAfWT9NixYy1eXHlTkjEUuw+Uf/+F6dNh+XIYMgSOHIFGjbSuSohiURSFCRsmlOswARMCpXfv\n3lStWhVPT08qVJAGjSWVqxbKqVPqgo2//gpDh6pb7DZooHVVQhRbXphsOLWhXIcJmBAoRqORAwcO\nWKOWcq9cjKHExqpLyP/+OwwfDsePQ716WlclRIkoisJ7G94j/FQ4G/03luswARPGULp168Zff/1l\njVrKvTLdQjl6FF5+WV1ny9lZDZYpUyRMhN3KC5OwU2FseLl8t0zyFBkoDz74IP369aNKlSrUrFmT\nmjVr4uTkZJYnnz17NhUqVODixYv5xwIDA3F1dcXNzY2wsLD841FRUXh6euLq6sqYMWPM8vy2pkyO\noRw+DH5+8Oij6k2IJ0/CxInqtrtC2ClFUXh/4/v5YVKvmnwwAhMCZezYsezcuZNr166RmppKamoq\nV65cKfUTx8fHEx4eTtObpoPGxMSwfPlyYmJiCA0NZfjw4flvsMOGDWPBggXExsYSGxtLaGhoqWuw\nNWWqhbJ/Pzz/PDzxBLRrpwbJhx9CrVpaVyZEqeSFyV8n/5IwuU2RgdKkSRM8PDzMPiA/duxYPv/8\n81uOhYSE4Ofnh6OjI87Ozri4uBAZGUlSUhKpqal439iW0N/fnzVr1pi1HltQJsZQoqOhXz91j/bO\nndXB9wkTbLBQIYovL0xCT4RKmBSgyEH5Zs2a8fjjj9OjRw8qVaoElH7acEhICAaDgTa3LTOemJhI\n586d8/9uMBgwGo04OjpiuOl+BL1ej9FoLPHz2yq7bqHs2gWffqoGyoQJ8MsvULWq1lUJYTaKovDB\nxg8IPRHKRv+NEiYFMClQmjVrRmZmJpmZmSZf2NfXl+Tk5DuOT506lcDAwFvGR+x1D3hzs8sxlORk\neP11tYvrvfdg5Up1TX0hypCM7AzGh4/nn3//YYO/tEwKU2SgTJo0qUQXDg8PL/D4oUOHiIuLo23b\ntgAkJCTQvn17IiMj0ev1xMfH55+bkJCAwWBAr9eTkJBwy3G9Xm9SzT4+Pvj4+JToZ7A2XTHPv3RJ\n47Htbdtg4EB45RVYtQoqV9awGCEs4+j5owxaPYgmtZqwafCmMjObKyIigoiICPNeVNGYs7OzcuHC\nBUVRFOXw4cNK27ZtlYyMDOXUqVPK/fffr+Tm5iqKoije3t7Kzp07ldzcXKVHjx7K+vXrC7yeDfxI\nJeYRGakcTE01+fxatRQlJcWCBRUmN1dR5sxRlAYNFOX33zUoQAjLy83NVebtnqfU/7y+Mm/3vPz3\norLKHO+dmm+wpdP9/+dyd3d3BgwYgLu7Ow4ODgQFBeV/PygoiICAANLT0+nZsyfdu3fXqmSbkJ2t\nrjZsphncpktLU+9sP3QIduyA5s2tXIAQlnf+2nle//11Tl86zT8B/9CqQSutS7ILuhvJVGbodDq7\nHZNpvWsXy9zdaW3CJifnzqkL8p4/b4XC8pw4Af37q9OA582DatWs+ORCWMeGUxsIWBPAC61fYOoT\nU6nsUD66cs3x3lnsFsq3335L/fr1efbZZ3Fw0LyBU25dvGjlm8zXrlUXb5w8Gd58U10RWIgyJDMn\nkw83fsgvh37hp2d+wre5r9Yl2Z1i31yiKApbtmyhX79+lqhHmOjCBahrjbHBnBz46CMYMUINlWHD\nJExEmXP0/FE6/9iZ4xePs//N/RImJVTsJsbIkSMtUYcoposXrRAo58/DoEHqgE1UFDRsaOEnFMK6\nFEXhh+gf+HDTh3z6+KcMbT/0lnFdUTxFBsr169dZvXr1HfuhfPLJJxYvThTO4l1ee/bAc8+p04Kn\nTgXp3hRlzIVrF3j999c5lXJKBt7NpMgur2eeeYa1a9fi6OhIjRo1qFGjBtWrV7dGbeIuLNrl9eOP\n0KMHzJ4NM2ZImIgyZ+OpjbSd15ZmtZsROSRSwsRMTNoPRZavtz0W6fK6fh1GjoTt22HLFnV1YCHK\nkMycTD7a9BFLDi5h0TOL6Na8m9YllSkmLV8vG2zZHrN3eZ0+DQ8/DFeuQGSkhIkoc46dP0aXBV04\nduEY+4bukzCxgEJbKJ6engDk5OSwaNEimjVrRuUbS2vodDoJGY2ZtYXy118weDC8+y68/bbM4hJl\niqIo/Bj9Ix9s+kAG3i2s0ED5/fffgYJvdpFfhvbMMoaSm6tuxxsUBMuXw2OPmaU2IWyFDLxbV6GB\n4uzsDMDLL79McHDwLd8r6JiwrlK3UC5dAn9/NZl274a7LLYphD3aeGojg9cMZoDHAJY+u7Tc3PGu\npSIH5Q8dOnTL37Ozs4mKirJYQcI0pRpDOXBAXUKlZ091leAb+9wIURbIwLt2Ch2UnzZtGjVr1uTg\nwYP5e8nLct+FAAAgAElEQVTXrFmThg0b0qdPH2vWKApQ4i6vn3+Grl3VJVTmzJEwEWXKsfPHeHDB\ngxw9f1QG3jVQ5OKQ7733HtOnT7dWPaVWHhaHzMpSN0PMzASTd2bOzIR33oHQUFi9Gm7bLVMIe3bz\nwPsUnym82eFNGestJqssDhkYGMjq1avZunUrFSpU4OGHH5Z1vDSWt7GWyWFiNMLzz0P9+up4Se3a\nFq1PCGu6eeD974C/cW/grnVJ5VaRb0nDhw9n/vz5tGnTBg8PD+bNm8fw4cOtUZsoRLG6uyIioGNH\nePppWLNGwkSUKZviNtFufjucazsTOSRSwkRjRbZQNm/eTExMDBVufBwOCAjA3V1+aVoyaYaXosAX\nX8DMmbB4MXSTvmRRdmTmZPLxpo/5+eDPLOyzkKdcntK6JIEJLRQXFxf++++//L//999/uLi4lOpJ\nJ02ahMFgwMvLCy8vL9avX5//vcDAQFxdXXFzcyMsLCz/eFRUFJ6enri6ujJmzJhSPb+9K3KGV2oq\nDBgAy5apd71LmIgy5PiF4zy44EFizsewb+g+CRMbUmSgXLlyhVatWvHYY4/h4+ODu7s7qamp9O7d\nu8SzvXQ6HWPHjmXv3r3s3buXHj16ABATE8Py5cuJiYkhNDSU4cOH5w8SDRs2jAULFhAbG0tsbCyh\noaEleu6y4K5dXkeOgLe3OsiyZQs0bWrV2oSwlLyB94cWPsSrXq+y9oW1NKjeQOuyxE2K7PKaMmVK\nod8rzSyKgmYThISE4Ofnh6OjI87Ozri4uBAZGUnTpk1JTU3F29sbAH9/f9asWVNu95UvtMtr1Sp1\nA6zp0+G116xelxCWcjH9Iq///jonLp4gYnAEHg09tC5JFKDIQPHx8eH06dOcOHGCJ598kmvXrpGd\nnY2Tk1Opnvibb75h8eLFdOjQgdmzZ1O7dm0SExPp3Llz/jkGgwGj0YijoyMGgyH/uF6vx2g0lur5\n7dkdgZKdDe+/rwZKaCi0b69ZbUKY2+a4zfiv8ee5Vs+xpP8SqjhU0bokUYgiu7y+//57nn/+eYYO\nHQpAQkKCSdOGfX198fT0vONr7dq1DBs2jLi4OPbt28c999zDO++8U/qfpAwwdQb4hQu3jaG89JJ6\n9/uePRImosy4cO0Cw/8czku/vcQPvX/gy+5fSpjYuCJbKN9++y27du3Kbzm0aNGCs2fPFnnh8PBw\nkwoYMmQIvXv3BtSWR3x8fP73EhISMBgM6PV6EhISbjmuv8vaU5MmTcr/s4+PDz4+PibVYgtM6UY8\nexYef/zGX379FfbtU7+qyItN2L/s3Gzm7ZnHlL+nMNBjIAeHHaRuVUvvd13+REREEBERYdZrFhko\nlStXzl+2HtS1vEp7B2pSUhL33HMPAL/99lv+Uvl9+vRh0KBBjB07FqPRSGxsLN7e3uh0OpycnIiM\njMTb25vg4GBGjx5d6PVvDpSy6MyZG9u7X7oEo0aps7kkTEQZsCluE2NCx9CwekM2Dd5E64attS6p\nzLr9w/bkyZNLfc0iA+Wxxx5j6tSpXLt2jfDwcIKCgvJbFCU1YcIE9u3bh06no1mzZsyfPx8Ad3d3\nBgwYgLu7Ow4ODgQFBeWHV1BQEAEBAaSnp9OzZ89yOyAPaqA0agSMHw/PPAOPPKJ1SUKUyulLpxkX\nNo6opChmd5tNP7d+snSKHSpyLa+cnBwWLFiQf0/IU089xZAhQ2z2l23Pa3l57NrFCg8PPKpXv+t5\ntWtD/OLN1BzhD4cOQa1aVqpQCPNKy0xjxrYZfLv7W97q9BbjHhxHVceqWpdVLlllLa+KFSvSt29f\n+vbtS8OGDUv1ZKL0rl+H3LR0arzzBnz7rYSJsEuKorD88HLeDX+Xh5o8xL6h+7iv1n1alyVKqdBA\nURSFyZMnM3fuXHJycgA1XEaNGsUnn3xisy2Usu7sWZhWZQo6Ly+QbQSEHdqbtJcxoWNIzUxlSf8l\nPNJUumzLikKnDX/55Zds27aN3bt3k5KSQkpKCrt27WLbtm18+eWX1qxR3CR1yz4GXV8A33yjdSlC\nFMv5a+d584836bGkBy+1eYk9r++RMCljCg2UxYsX88svv9CsWbP8Y/fffz9Llixh8eLFVilO3CY7\nG8PE1/jJbcaNUXkhbF9WThZzIufQ6ttWVHGowpERR3ij/RtUrFBR69KEmRXa5ZWdnU2DBneuk9Og\nQQOys7MtWpQoxFdfkVapDgceCNC6EiFMsuHUBsaEjuHemvfKkinlQKGB4ujoWOiD7vY9YSEnT8L0\n6ax9JZJGFWT8Sti2UymneCfsHfYn7+eLp77gmZbPyLhrOVBooBw4cICaNWsW+L309HSLFSQKoCgw\ndCi89x7HjM25aVkzIWzK1cyrTN86nXl75jG2y1iWPrtUlkspRwoNlLyZXcIG/PQTpKTAW29xxl+W\n6xK2R1EUlh5ayoQNE3i06aPse3MfBif55FPeFHkfitBYcjJMmABhYeDgwNmzMh4vbEt0UjSj148m\nPTudZc8u46EmD2ldktCIBIqtGzNG3dukXTvgpnW8hNDY2bSzfLjxQ34//jufPfEZr7R7pcQzt+rW\nrUtKSoqZKxQFqVOnDhcvXrTItSVQbNnatRAdrXZ53ZC/jpcQGsnKyWLurrlM2zqNlzxf4ujIo9Su\nUrtU10xJSbHbJZPsjSUnR0ig2KrLl2HECFi8GKqqaxvl5Kiba9Wvr3FtotwKOxnGmNAxNKnVhH8C\n/qFVg1ZalyRsiASKrXr/feje/aaNT+D8eXWreAf5rQkrO3nxJGPDxnL47GG+eOoLerfoLdOAxR3k\nrckWbd0KISFw+PAth6W7S1jb1cyrTP1nKj9E/8C4B8ex4rkVVHaoXPQDRbkkgWJrMjLg9dfVtbpq\n39ovffasDMgL61AUhSUHlzBhwwSeaPYEB4Yd4N6a92pdlrBxEii25ocfoFUr6N//jm9JC0VYw57E\nPYxeP5qs3CxWPb+KLvd10bokYScKXRzS0r755htatWpF69atmTBhQv7xwMBAXF1dcXNzy9/UCyAq\nKgpPT09cXV0ZM2aMFiVbXkYGLF8Oc+cW+G0JFGFJZ66e4bWQ1+i9tDdDHhhC5JBICRML2bdvH+PG\njTPb9UJCQliyZAlTpkwhKCjIbNctLk1aKJs3b2bt2rUcOHAAR0dHzp07B0BMTAzLly8nJiYGo9HI\nk08+SWxsLDqdjmHDhrFgwQK8vb3p2bMnoaGhZWsb4JwcSExU7zu5t+CuBenyEpaQmZPJN5HfMH3b\ndAa3HczREUepVUU2brOUL774gq1bt1LLTJvjXbp0iYEDB3Lp0iUqV65M/fr16dWrF02bNjXL9YtD\nkxbKd999x/vvv5+/yGTeqsYhISH4+fnh6OiIs7MzLi4uREZGkpSURGpqKt7e3gD4+/uzZs0aLUq3\nnLlzQacrsKsrj7RQhLmtj11Pm+/asCFuA1tf2cqsbrMkTCxs7NixPPPMM2a7Xu3atYmKiqJKlSro\ndDqys7M1u6dHkxZKbGws//zzDx988AFVqlRh1qxZdOjQgcTERDp37px/nsFgwGg04ujoiOGmFRH1\nej1Go1GL0i3j9Gn49FNYswYqFJ7xEijCXGIvxDI2bCxHzx/ly6e+pJdrL5kGXEKnTp3ihx9+KPT7\nnTt3viNAzP2G7+GhbguwdetWfHx8cHZ2Nuv1TWWxQPH19SU5OfmO41OnTiU7O5uUlBR27tzJ7t27\nGTBgAKdOnTLbc0+aNCn/zz4+Pvj4+Jjt2hYxezYMGwaVKt31NFl2RZRW/OV4PvvnM1YdWcWEhyaw\n6vlVdjEN2FxZV9L38ePHj/PRRx9x7tw59uzZg4+PD7169eLNN9/k/vvvJzAwsFjXKyi8Dx8+zOLF\ni3n00UeJiorik08+KdY1f/31V1auXMns2bNNOj8iIoKIiIhiPUdRLBYo4eHhhX7vu+++o/+Nrp2O\nHTtSoUIFzp8/j16vJz4+Pv+8hIQEDAYDer2ehISEW47r9fpCr39zoNiF3bth1qwiT0tKgrv82EIU\nKik1iWlbprHk4BLeaP8Gx0cep161elqXZTItV2W5ePEib775JuvWraNKlSr07duX//3vf6UaA7m9\nhXL27Fl69erF7t27adCgAdu2bQPg888/L3S7kMGDB9/SEunfvz/dunXDy8uL8PDwIlspt3/Ynjx5\ncol+lptp0uXVt29fNm3axGOPPcbx48fJzMykfv369OnTh0GDBjF27FiMRiOxsbF4e3uj0+lwcnIi\nMjISb29vgoODGT16tBalm192Nhw8CG3bwpEjdz3t/Hnp8hLFczbtLDO2zmDRvkUEtAvgyIgjNKoh\n/xMVx7fffsuIESOoUkXd1yUjI4Nq1arlf78kXV63t1BWrlxJ06ZN2bt3L+fOnWPUqFEAvPvuu0XW\n9+effzJt2jS2bdtGjRo1aNiwIatWrTLrLDJTaRIor776Kq+++iqenp5UqlQpf496d3d3BgwYgLu7\nOw4ODgQFBeX/wwcFBREQEEB6ejo9e/YsOzO8jh4FgwEK2cwsz9mzULeuLLsiTHMx/SKzts9i3p55\nDPIcxKHhh+TGxBJKTU3F3d0dULulPDw8btm1tiRdXre3UKpWrUqPHj3o1q0bAMnJyWRkZFC5ctHd\nkRUrVsxvaSiKQnx8PG3atClWPeaiU8rYEp86nc6+Vi0NDoZ162DpUjx27WKFhwce1avfcdqePfDG\nG+riw0IU5vL1y3y580vm7ppLP7d+fPToRzStbf3po8Vly6/buLg41q5di8FgICEhgREjRuBQik92\nc+fOZcWKFcTHxxMQEMDbb7+No6MjU6dOpXPnzuTm5pKTk0O/fv1MvmZQUBA5OTn8+++/uLq6MnTo\n0ELPLezf2hy/AwkUrb39NtxzD7z77l0DZe1a+P57+OMPDWoUNu9q5lXmRM7hy51f0tO1J588+gnN\n6zbXuiyT2d3r1o5ZMlA0u1Ne3BAdDQ88UORpiYmF3u8oyrFrWdeYtX0Wzec05+DZg2x5ZQv/6/s/\nuwoTUXZIj7yWcnNh3z7w8iryVAkUcbPr2df5Pup7pm+dTpf7urDRfyOtG7bWuixRzkmgaOnUKXWD\nk3pFT99MTIQbCwWIciwzJ5NFexcxdctU2jZuy5+D/sTrnqI/kAhhDRIoWoqONql1Auo9KNJCKb+y\nc7MJ3h/MlH+m4FrXlRXPr6CzoXPRDxTCiiRQtGTi+AlIl1d5lZObw7JDy5j892TurXkvi/su5pGm\nj2hdlhAFkkDR0t69YOINmhIo5UuuksuvR35lYsREalWuxXe9vuOJZk/IelvCpkmgaEVRTO7yysqC\nixfhxqLMogxTFIXfj//OJ5s/waGCA7N8Z9HdpbsEibALEihaSUiAihXVe1CKkJysLgpZsaIV6hKa\nUBSFv07+xSebPyEjJ4MpPlPo07KPBImwKxIoWtm7Vx0/MeENQ7q7yrZNcZv4ePPHpKSnMNlnMs+6\nP0sFndwiJuyP/F+rFRmQL/e2/reVx//3OEP/GMqwDsM4OOwgz3s8L2FSDph7C2BLXbO4pIWileho\nGDzYpFMlUMqWXcZdfLz5Y46dP8Ynj32Cf1t/HCrIS7G8MPcWwJa6ZknIRyGtSAul3NmXvI8+S/vw\n7Ipn6efWj+OjjvOq16sSJuWMubcAttQ1S0L+T9bC2bOQlgYmbtOZlAQPPWTZkoTlHD57mIkRE9kW\nv433HnqPFc+voIpDFa3LEmZiC1sAW+qaxSWBooW9e9XpwibO4JEWin06fuE4kyImsTFuI+O6jGNx\nv8VUc6xW9APFHXSTzTPbTZlYsjdde9gC2BZmBEqgaKEY3V2gBooJs4uFjTiVcoopf0/hz9g/eavT\nW8x/ej41K999AzVxdyUNAnOwly2Ay20L5YUXXuDYsWMAXLp0idq1a7N3714AAgMDWbhwIRUrVmTO\nnDn5O5hFRUUREBDA9evX6dmzJ19//bUWpZvH3r3Qt6/Jp0sLxT78d/k/PvvnM1YfWc3IjiOJHRVL\n7Sq1tS5LlJKtbwFc2DW1oEmgLFu2LP/P48aNo3Zt9UUXExPD8uXLiYmJwWg08uSTTxIbG4tOp2PY\nsGEsWLAAb29vevbsSWhoqP1uAxwdDVOmmHRqRgZcuQL161u4JlFiSalJTNsyjV8O/cIbD7zB8ZHH\nqVet6BWkhX2w9S2AC7umFjSd5aUoCitWrMDPzw+AkJAQ/Pz8cHR0xNnZGRcXFyIjI0lKSiI1NRXv\nG+u3+/v7s2bNGi1LL7lLl+DMGXB1Nen0pCRo3BgqyHw8m3M27Szv/PUOHkEeVKpYiSMjjhD4ZKCE\nSRkzbNgwwsLCWL16NRs2bGD69Omlut7cuXNZuHAhERERTJ48mStXruDn58fVq1f5448/WLt2LTt2\n7ChWmBR0TS1oOoayZcsWGjVqRPPm6u5yiYmJdO78/0tyGwwGjEYjjo6OGAyG/ON6vR6j0Wj1es1i\n3z5o08bkdVSku8v2XLh2gVnbZ/F99Pf4tfbj0PBD3FtTfkllVbNmzRgzZozZrjdy5EhGjhx5x/HP\nPvvM7Ne0NosFiq+vL8nJyXccnzZtGr179wZg6dKlDBo0yOzPPWnSpPw/+/j44OPjY/bnKLHERLjv\nvmKdLoFiGxKuJDB7+2z+t/9/POf+HHuH7qVJrSZalyVEiURERBAREWHWa1osUMLDw+/6/ezsbH77\n7Teio6Pzj+n1euLj4/P/npCQgMFgQK/Xk5CQcMtxvV5f6LVvDhSbU7eu2u1looQEuKlxJjQQeyGW\nGdtm8OuRXwloF8CBYQcwOMkvRdi32z9sT548udTX1KxnfsOGDbRq1Yp7b/r43adPH5YtW0ZmZiZx\ncXHExsbi7e1N48aNcXJyIjIyEkVRCA4Opm8xZknZlLp14cIFk09PSChWg0aY0b7kfQxcNZAHFz6I\nwclA7KhYvnjqCwkTIQqh2RjK8uXL8wfj87i7uzNgwADc3d1xcHAgKCgofypcUFAQAQEBpKen07Nn\nT/ud4VW3rrq5iYkSEqB9ewvWI+6w5d8tBG4NZP+Z/YztPJYfe/8o95EIYQKdYgtzzcxIp9PZxPS5\nQqWkwP33q/+9jceuXazw8MCjevX8Yw8/DIGB8Ijs+mpRiqKw/sR6ArcGkpSaxLsPvcvgtoOp7GD6\nTBtRcjb/ui1DCvu3NsfvQO6Ut7ZatSA1FbKzwaHof/74eBlDsaSc3BxWxawicGsguUou7z/8Ps97\nPC8LNgpRAvKqsbYKFaB2bXVgvoi7FXNy1PtQZJaX+WVkZxB8IJgZ22bQsHpDPnviM3q59rKJu42F\nsFcSKFrIG5gvIlDOnoU6daAY9zeJIlzNvMr3Ud/zxY4vaN2wNQv6LOCRJo9IkAhhBhIoWjBxYF6m\nDJvPxfSLfBP5Dd/u/hYfZx/W+q3lgXtMX6BTCFE0CRQt1KsngWIliamJfLHjCxbuXUg/t35seWUL\nLeu31LosUc7t27ePn3/+mVmzZt31vDVr1hATE0OFChXQ6/W8/PLLVqqwZCRQtGDivSgSKCV38uJJ\nPt/2OStjVuLf1p/9b+7nvlpyQ4/Qnqnb9V6+fJlPP/2UqKgoALp06UKPHj2ob8MrxcqSg1qQLi+L\nOXDmAINWD6LTj51oWL0hx0Ye46vuX0mYCJth6na9//zzT/4qxwBt27Zl8+bNliyt1KSFooV69Uxu\nobRubYV6yoDt8duZtmUaUUlRvN35beY9PQ+nyk5alyXKAUttAZyQkJC/tQdA7dq1iY2NLXmhViCB\nooW6deHIkSJPkxbK3SmKQtjJMKZtncZ/l//j3QffZdWAVbJfe1lkrll4Jbxxz9pbAO/Zs4dKlSrl\nb+oFUKlSJa5evVqi+q1FAkULxWihSKDcKSc3h9+O/kbg1kAysjN4/+H3Gdh6oNyMWJZpeBe9VlsA\n16tXjws3vU+kp6fTqFGjEj+nNcgrUAsmjKEoChiNcJdFlcudzJxMfj7wMzO2zaBOlTpMfGwiT7d4\nmgo6GQoUlqPFFsAjR45k//797NmzJ/+c8+fP88ADtj3VXQJFCyYEyvnzUL063PT/bbmVlpnGj9E/\nMnvHbNzquzGv1zx8nH3kZkRhFVptAfzYY4/dsqd8dHQ0M2bMKOmPYRWyOKQWTp2Crl0hLu6Wwzcv\nDrl3LwQEwP792pRoC1LSU/h297d8s+sbHm7yMO899B4d9R21LktYgC2/buPi4li7di0Gg4GEhARG\njBiBgwnr8BVm7ty5rFixgvj4eAICAnj77bdxdHRk6tSpdO7cmdzcXHJycujXrx/BwcH8+++/5Obm\n0rx5c1588cVS/zyWXBxSAkULly5B06Zw+fIth28OlN9/h/nz4Y8/NKpRQ8lXk/lyx5f8uPdHerfo\nzYSHJtCqQSutyxIWZBev2zJCVhsua2rVgrQ0yMqCm5rONyuPA/JxKXHM3D6TZYeW8aLni0S/EU3T\n2k21LksIYSIJFC3odOqqjykp0LBhgaeUp0A5fPYw07dNZ33set5o/wZHRx6lYfWC/12EELZLk+kx\nu3btwtvbGy8vLzp27Mju3bvzvxcYGIirqytubm6EhYXlH4+KisLT0xNXV1fGjBmjRdnmVcTAfHkI\nlJ0JO3lm2TN0XdwVjwYenBx9kmldp0mYCGGvFA089thjSmhoqKIoirJu3TrFx8dHURRFOXz4sNK2\nbVslMzNTiYuLU5o3b67k5uYqiqIoHTt2VCIjIxVFUZQePXoo69evL/DaGv1Ixdeli6Js3XrLIffI\nSOXQ1auKoijKE08oSni4FoVZVm5urhJ2Ikx5/KfHlSZfNlHmRs5VrmVe07osoTG7ed2WAYX9W5vj\nd6BJl9c999zD5RsD0pcuXUJ/42aLkJAQ/Pz8cHR0xNnZGRcXFyIjI2natCmpqal4e3sD4O/vz5o1\na+x3X3kody2UXCWXkKMhTNs6jbTMNN57+D38WvvhWLHgMSQhhP3RJFCmT5/Oww8/zLhx48jNzWXH\njh0AJCYm0rlz5/zzDAYDRqMRR0dHDDe9u+r1eoxGo9XrNqu7rDisKGqglIWbGjNzMvnl4C98vu1z\nalSqwQcPf8Azbs/IzYhClEEWCxRfX1+Sk5PvOD516lTmzJnDnDlz6NevHytXruTVV18lPDzcbM89\nadKk/D/7+Pjg4+NjtmubjbMznDhR4LdSUtTJXzVrWrckc7p8/TLzo+YzJ3IOrRq0Yk6POXRt1lVu\nRhTCRkRERBAREWHWa1osUO4WEC+99BIbNmwA4LnnnmPIkCGA2vKIj4/PPy8hIQGDwYBerychIeGW\n4/q7fHy/OVBslpcX/Phjgd+Kj4f77HS19YQrCXy982sW7ltId5fu/O73O173eGldlhDiNrd/2J48\neXKpr6lJv4OLiwt///03AJs2baJFixYA9OnTh2XLlpGZmUlcXByxsbF4e3vTuHFjnJyciIyMRFEU\ngoOD6du3rxalm4+XF+zdW+C37DFQDp45yOA1g2nzXRuyc7OJfiOaJf2XSJgIUY5oMoby/fffM2LE\nCDIyMqhatSrff/89AO7u7gwYMAB3d3ccHBwICgrK7yIJCgoiICCA9PR0evbsad8D8qDeKX/9Opw5\nA7etIPrff/YRKIqisPn0ZmZun8n+5P2M8h7FV6O/ok7VOlqXJoRNM3UL4OKeqzVNAqVDhw5ERkYW\n+L0PPviADz744I7j7du35+DBg5YuzXp0uv9vpdwWjvHx0KSJRnWZIDs3m1Uxq5i5fSbXsq4xrss4\n1gxcQ2WHylqXJoTNM3UL4OKeawvkTnkt3SVQbiw6alPSMtNYuHchX+z8AoOTgUmPTaJXi14yY0uI\nYhg7diz16tUzaUC8OOfaAgkULXl5QUjIHYdtbQzlzNUzzN01l3lR83i06aMsfXYpnQ2di36gEOWA\npbYALsm5WpNA0ZKXF0yceMdhWwmU4xeOM3v7bFbGrGSgx0C2v7od13quWpclyiGdmT6hKyW8hcDa\nWwBHRUXxySefFHqurZJA0VLLlpCUpC5jf6OPNDdX3alRy7vkt8dvZ+b2mWz7bxvDOgyTxRqF5koa\nBOag1RbAhZ1ryyRQtFSxInh6qrtoPfoooK7GUqsW3Nht1GpylVzWHlvLzO0zSb6azNjOY1nSfwnV\nHGXLSFG+abEF8KhRowo915ZJoGgtb2D+RqAkJVu3u+t69nUW71/M7B2zcarsxLsPvkv/Vv2pWKGi\n9YoQwoZptQVwRkYGlStXtqsWikzP0doDD9xyg2NyknWmDF9Mv8jUf6bi/JUzIcdCmP/0fHYN2cXz\nHs9LmAhxk2HDhhEWFsbq1avZsGED06dPL9X15s6dy8KFC4mIiGDy5MlcuXIFPz8/rl69yh9//MHa\ntWvZsWMHlStXLvBcWyZbAGttzx547TXYvx+PXbt4OtqD9JjqzJljmac7fek0X+74kuADwTzj9gzj\nuozDo6GHZZ5MCBPZ3evWjskWwGVZ69YQGwsZGYDaQmltgS6v6KRoZm6fSdjJMIZ4DeHgsIPoncrA\ncsZCCJshgaK1KlXAxQUOHQIgORl6PG6eSyuKwl8n/2Lm9pkcv3Cctzq9xfyn5+NU2ck8TyCEEDeR\nQLEFeQPzbdqQbIZB+aycLJYdWsbM7TMBGP/geAa2HkilipXMUKwQQhRMAsUWeHlBdHSpA+VKxhV+\niPqBryK/omW9lsz0nUm35t3satqhEMJ+ySwvW3CjhaIo6n0o995bvIcnpiYyIXwC9399P3uS9hDy\nQggb/DfwlMtTEiZCCKuRFootaNcODh4kOxvq1QMHE38rh88eZtaOWYQcDeHlNi+z5409ONd2tmip\nQghRGAkUW1CrFjRuTM61DBo3vvupiqLw979/M3P7TKKTohnZcSQnRp+gbtW61qlVCAuoU6eOtKat\npE4dy+1XpEmX1/79++nSpQtt2rShT58+pKam5n8vMDAQV1dX3NzcCAsLyz8eFRWFp6cnrq6ujBkz\nRlvpjjgAAAipSURBVIuyLcvLi9xr1wsNlOzcbFYcXkGnHzsx9I+h9G3Zl7gxcXz46IcSJsLuXbx4\nEUVR5MsKXxcvXrTY71GTQBkyZAiff/45Bw4coF+/fsycqc5GiomJYfny5cTExBAaGsrw4cPzb7QZ\nNmwYCxYsIDY2ltjYWEJDQ7Uo3XK8vNBlXIfr/9xy+FrWNb7d9S0tvmnBnMg5fPjIhxwZcYTX279O\nFQcrL/hVBHvZs6EwUr+2pH77p0mgxMbG8sgjjwDw5JNPsnr1agBCQkLw8/PD0dERZ2dnXFxciIyM\nJCkpidTUVLy9vQHw9/dnzZo1WpRuOV5eVMy8ztXzWwA4l3aOiZsn4vyVMxviNvBz/5/Z+upWnnF7\nxmY3tLL3F5TUry2p3/5p8s7k4eFByI2NpVauXEl8fDwAiYmJGG5at91gMGA0Gu84rtfrMRqN1i3a\n0ry8cMi+ToWqlxn2xzBazm1J8tVktr66ld8G/saD9z2odYVCCHFXFhuU9/X1JTk5+Y7j06ZNY+HC\nhYwePZpPP/2UPn36UKmS3HBH48bk6nI5dH4Vj1V7gyMjjtCoRiOtqxJCCNMpGjt27Jji7e2tKIqi\nBAYGKoGBgfnfe+qpp5SdO3cqSUlJipubW/7xX375RRk6dGiB12vevLkCyJd8yZd8yVcxvpo3b17q\n93NNpg2fO3eOBg0akJuby2effcawYcMA6NOnD4MGDWLs2LEYjUZiY2Px9vZGp9Ph5OREZGQk3t7e\nBAcHM3r06AKvfeLECWv+KEIIIW7QZAxl6dKltGzZklatWmEwGAgICADA3d2dAQMG4O7uTo8ePQgK\nCsqfmx4UFMSQIUNwdXXFxcWF7t27a1G6EEKIQpS5/VCEEEJowzbnnxbh4sWL+Pr60qJFC7p168al\nS5cKPC80NBQ3NzdcXV2ZMWNG/vHx48fTqlUr2rZtS//+/bl8+bLFay6slpuNHj0aV1dX2rZty96b\ndnE05bGWVtL64+Pjefzxx/Hw8KB169bMsdTOYUUozb8/QE5ODl5eXvTu3dsa5d6hNPVfunSJ5557\njlatWuHu7s7OnTutVTZQutoDAwPx8PDA09OTQYMGkXFj3yBrKqr+o0eP0qVLF6pUqcLs2bOL9Vhr\nKGn9JXrtlnoURgPjx49XZsyYoSiKokyfPl2ZMGHCHedkZ2crzZs3V+Li4pTMzEylbdu2SkxMjKIo\nihIWFqbk5OQoiqIoEyZMKPDx5nS3WvL8+eefSo8ePRRFUZSdO3cqnTp1Mvmxllaa+pOSkpS9e/cq\niqIoqampSosWLeyq/jyzZ89WBg0apPTu3dtqdecpbf3+/v7KggULFEVRlKysLOXSpUt2UXtcXJzS\nrFkz5fr164qiKMqAAQOUn376yWq1m1r/2bNnld27dysffvihMmvWrGI91pbrL8lr1y5bKGvXrmXw\n4MEADB48uMCbHHft2oWLiwvOzs44Ojrywgsv5N/74uvrS4UK6o/eqVMnEhISLFrv3Wop6Gfq1KkT\nly5dIjk52aTHWlpJ6z9z5gyNGzemXbt2ANSoUYNWrVqRmJhoN/UDJCQksG7dOoYMGaLJNrWlqf/y\n5cts2bKFV199FQAHBwdq1aplF7U7OTnh6OjItWvXyM7O5tq1a+j11t1l1JT6GzRoQIcOHXB0dCz2\nYy2tNPWX5LVrl4Fy5swZGjVS79Fo1KhR/gv/Zkajkftu2lgk7ybJ2y1cuJCePXtarlgTaynsnMTE\nRJN+Dksqaf23B/Xp06fZu3cvnTp1smzBtynNvz/A22+/zcyZM/M/hFhbaf794+LiaNCgAa+88goP\nPPAAr7/+OteuXbP52o1GI3Xr1uWdd96hSZMm3HvvvdSuXZsnn3zSarXfrTZLP9ZczFWDqa9dmw0U\nX19fPD097/hau3btLefpdLoCVyk1ZeXSqVOnUqlSJQYNGmS2ugti6iqqWnz6NUVJ67/5cVevXuW5\n557j66+/pkaNGmatryglrV9RFP744w8aNmyIl5eXZr+f0vz7Z2dnEx0dzfDhw4mOjqZ69epMnz7d\nEmUWqDT/7588eZKvvvqK06dPk5iYyNWrV1myZIm5S7yr0qyAbAurJ5ujhuK8dm12+frw8PBCv9eo\nUSOSk5Np3LgxSUlJNGzY8I5z9Hp9/pIuoA4w3bx8y08//cS6devYuHGjeQsvQFG1FHROQkICBoOB\nrKysIh9raSWtP697Iisri2effZaXXnqJvn37Wqfou9RWnPpXr17N2rVrWbduHdevX+fKlSv4+/uz\nePFiu6hfURQMBgMdO3YE4LnnnrNqoJSm9oiICB588EHq1asHQP/+/dm+fTsvvviidYovoLbivP5K\n81hzKW0NxX7tmm30x4rGjx+vTJ8+XVEU9e76ggbVs7KylPvvv1+Ji4tTMjIybhmMWr9+veLu7q6c\nO3fOKvXerZY8Nw9M7tixI39g0pTH2nL9ubm5yssvv6y89dZbVq35ZqWp/2YRERHK008/bZWab1ba\n+h955BHl2LFjiqIoysSJE5V3333XLmrfu3ev4uHhoVy7dk3Jzc1V/P39lblz51qtdlPrzzNx4sRb\nBrXt5bWb5/b6S/LatctAuXDhgtK1a1fF1dVV8fX1VVJSUpT/a9+OTSQEwiiOsw2IBajJyAYiarBV\nmNmA+WIfgh1YiaaCsWATsgWYmbzN5II77tid8074/0IZ4c3A8AI/JWlZFuV5vq/ruk7X61XGGNV1\nvT8Pw1BBECjLMmVZpvv9/uuZP8vStq3att3XVFUlY4ySJNE0Td/u40iv5h/HUZfLRWma7ufd9/1p\n8n80DMOfTHlJ7+Wf51m3201JkqgoikOnvN7N3jSNoihSHMcqy1Lbth2a/Sf5H4+HPM+T4zhyXVe+\n72td1y/fPUv+V+4uPzYCAKz4tx/lAQDnQqEAAKygUAAAVlAoAAArKBQAgBUUCgDACgoFAGAFhQIA\nsOIJ/Z33A0v7JjEAAAAASUVORK5CYII=\n", + "text": [ + "" + ] + } + ], + "prompt_number": 101 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [], + "language": "python", + "metadata": {}, + "outputs": [] + } + ], + "metadata": {} + } + ] +} \ No newline at end of file