{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "# add simpegseisc directory to path\n", "import os\n", "if '..' not in os.sys.path:\n", " os.sys.path.append('..')" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Populating the interactive namespace from numpy and matplotlib\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "WARNING: pylab import has clobbered these variables: ['linalg']\n", "`%matplotlib` prevents importing * from pylab and numpy\n" ] } ], "source": [ "from simpegseis import *\n", "%pylab inline\n", "from SimPEG.Utils import mkvc" ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "You are good to go:)\n", ">> Stability information\n", " dt: 8.00e-05 s\n", " Optimal dt: 1.09e-04 s\n", " Cell per wavelength (G): 1.60e+01\n", " Optimal G: 1.60e+01\n" ] } ], "source": [ "# Step1: set time\n", "dt = 0.8e-4\n", "ntime = 1700\n", "time = np.arange(ntime)*dt\n", "# Step4: set mesh\n", "cs = 0.5\n", "cs = 0.5\n", "hx = np.ones(500.)*cs\n", "hy = np.ones(200.)*cs\n", "mesh = Mesh.TensorMesh([hx, hy], 'CC')\n", "# Step2: set Tx and Rx\n", "tlag = 0.008\n", "options={'tlag':tlag, 'fmain':100.} # You need to set waveform to set Tx\n", "nrx = 99\n", "xtemp_rx = np.arange(nrx)*2.5+2.5+mesh.vectorCCx.min()\n", "xy_rx = np.c_[xtemp_rx, np.ones(nrx)*mesh.vectorCCy.max()]\n", "rx = AcousticRx(xy_rx)\n", "\n", "tx = AcousticSrc(np.r_[mesh.vectorCCx.min(), mesh.vectorCCy.max()], time, [rx], **options)\n", "# Step3: set survey (pass txlist)\n", "survey = SurveyAcoustic([tx])\n", "wave = tx.RickerWavelet()\n", "# Step5: set problem (pass mesh) and pair with survey\n", "prob = AcousticProblemPML(mesh)\n", "prob.pair(survey)\n", "# Step6: set boundary\n", "prob.setPMLBC(20, dt, bcflag='left', const=1.)\n", "prob.storefield = True\n", "# Step7: set velocity model and check stability\n", "v = np.ones(mesh.nC)*800.\n", "v[mesh.gridCC[:,1]<30.] = 2300.\n", "v[mesh.gridCC[:,1]<0.] = 1500.\n", "prob.stabilitycheck(v, time, 100.)" ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(,)" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEKCAYAAAD9xUlFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAADa1JREFUeJzt3X2MpXV5h/HrS7ekVqzBWhbLqruEF4H+oabZpjVp1lbe\nbAPUNBTTqFtq0ohakzZWF2wgpglqYm2jgT8KVmqshNq0oDGwEF6SmlhphIpdhE3jUtiwS9/SpGmT\nLnD3j/Nwz2Gd2R3WmfOcnXN9ksme8zvPztznmTl7zfOcszOpKiRJAjhh7AEkSfPDKEiSmlGQJDWj\nIElqRkGS1IyCJKltGnuAH0YSX08rScegqrLc+nEdhYnrxh4AuA9469hDzAn3xRL3xRL3xZJ52BfX\nrXiLp48kSc0oSJKaUVgTW8ceYI5sHXuAObJ17AHmyNaxB5gjW8ce4IiMwprYNvYAc8R9scR9scR9\nsWS+94VRkCQ1oyBJakZBktSMgiSpGQVJUjMKkqRmFCRJzShIkppRkCQ1oyBJakZBktSMgiSpzUUU\nkpyQ5NtJ7hiun5xkd5LHktyV5JVjzyhJi2AuogB8CNgzdf2jwD1VdTZwL7BrlKkkacGMHoUkW4C3\nAzdNLV8K3DJcvgW4bNZzSdIiGj0KwGeADwM1tba5qg4CVNUB4JQxBpOkRbNpzA+e5FeAg1X1cJId\nR9i0Vr7pvqnLW5n3X2AhSbP3fWDfqrYcNQrAW4BLkrwdeBnwiiRfBA4k2VxVB5OcCjyz8rt460wG\nlaTj1zZe/A3zAytuOerpo6q6uqpeV1WnA1cA91bVu4CvAjuHzd4D3D7SiJK0UObhOYXlfAI4P8lj\nwC8P1yVJ62zs00etqh5gOKapqv8A3jbuRJK0eOb1SEGSNAKjIElqRkGS1IyCJKkZBUlSMwqSpGYU\nJEnNKEiSmlGQJDWjIElqRkGS1IyCJKkZBUlSMwqSpGYUJEnNKEiSmlGQJDWjIElqqaqxZzhmSeq6\nsYeQpOPMdUBVZbnbPFKQJDWjIElqRkGS1IyCJKkZBUlSMwqSpGYUJEnNKEiSmlGQJDWjIElqRkGS\n1IyCJKkZBUlSMwqSpGYUJEnNKEiS2qhRSLIlyb1J/inJI0l+d1g/OcnuJI8luSvJK8ecU5IWxdhH\nCs8Cv1dV5wE/D7w/yRuAjwL3VNXZwL3ArhFnlKSFMWoUqupAVT08XP5v4FFgC3ApcMuw2S3AZeNM\nKEmLZewjhZZkK/BG4JvA5qo6CJNwAKeMN5kkLY65iEKSk4CvAB8ajhjqsE0Ovy5JWgebxh4gySYm\nQfhiVd0+LB9MsrmqDiY5FXhmpb9/39TlrcC2dZtUko5P3wf2rXLb0aMAfB7YU1V/OrV2B7AT+CTw\nHuD2Zf4eAG9d19Ek6fi3jRd/w/zAEbYdNQpJ3gL8JvBIkoeYnCa6mkkMbktyJfAEcPl4U0rS4hg1\nClX1DeBHVrj5bbOcRZI0J080S5Lmg1GQJDWjIElqRkGS1IyCJKkZBUlSMwqSpGYUJEnNKEiSmlGQ\nJDWjIElqRkGS1IyCJKkZBUlSMwqSpGYUJEnNKEiSmlGQJDWjIElqRkGS1IyCJKkZBUlSMwqSpGYU\nJEnNKEiSmlGQJDWjIElqRkGS1IyCJKkZBUlSMwqSpGYUJEnNKEiSmlGQJDWjIElqRkGS1OY6Ckku\nSvK9JI8n+cjY80jSRje3UUhyAvA54ELgPOCdSd4w7lSStLHNbRSA7cDeqnqiqg4BtwKXjjyTJG1o\nR41Ckg8mOXkWwxzmNODJqetPDWuSpHWyaRXbbAYeTPJt4PPAXVVV6zvW6t03dXkrsG2kOSRpXn0f\n2LfKbY8ahar6WJI/BC4Afgv4XJLbgJur6p+Pecqj2w+8bur6lmHtRbbVjS+6vm8dB5Kk41F48TfM\nD+R9K267qucUhiODA8Pbs8DJwFeSfOrYxzyqB4Ezkrw+yYnAFcAd6/jxJGnhHfVIIcmHgHcD/wbc\nBHy4qg4Nrw7aC/zBegxWVc8l+QCwm0m8bq6qR9fjY0mSJlbznMKrgHdU1RPTi1X1fJJfXZ+x+mPc\nCZy9nh9DkrRkNc8pXHuE2/zOXZI2kHn+fwqSpBkzCpKkZhQkSc0oSJKaUZAkNaMgSWpGQZLUjIIk\nqRkFSVIzCpKkZhQkSc0oSJKaUZAkNaMgSWpGQZLUjIIkqRkFSVIzCpKkZhQkSc0oSJKaUZAkNaMg\nSWpGQZLUjIIkqRkFSVIzCpKkZhQkSc0oSJKaUZAkNaMgSWpGQZLUjIIkqRkFSVIzCpKkNloUknwq\nyaNJHk7y10l+Yuq2XUn2DrdfMNaMkrRoxjxS2A2cV1VvBPYCuwCSnAtcDpwDXAzckCSjTSlJC2S0\nKFTVPVX1/HD1m8CW4fIlwK1V9WxV7WMSjO0jjChJC2denlO4Evj6cPk04Mmp2/YPa5KkdbZpPd95\nkruBzdNLQAHXVNVXh22uAQ5V1ZfXcxZJ0tGtaxSq6vwj3Z5kJ/B24JemlvcDr526vmVYW9ZD132t\nL5+64yxes+OsYxlVkjasp+9/nAP3P76qbVNV6zzOCh84uQj4NPCLVfXvU+vnAl8Cfo7JaaO7gTNr\nmUGT1M66cUYTS9LG8IW8j6pa9gU863qkcBSfBU4E7h5eXPTNqrqqqvYkuQ3YAxwCrlouCJKktTda\nFKrqzCPcdj1w/QzHkSQxP68+kiTNAaMgSWpGQZLUjIIkqRkFSVIzCpKkZhQkSc0oSJKaUZAkNaMg\nSWpGQZLUjIIkqRkFSVIzCpKkZhQkSc0oSJKaUZAkNaMgSWpGQZLUjIIkqRkFSVIzCpKkZhQkSc0o\nSJKaUZAkNaMgSWpGQZLUjIIkqRkFSVIzCpKkZhQkSc0oSJKaUZAkNaMgSWpGQZLUjIIkqY0ehSS/\nn+T5JK+aWtuVZG+SR5NcMOZ8krRINo35wZNsAc4HnphaOwe4HDgH2ALck+TMqqpxppSkxTH2kcJn\ngA8ftnYpcGtVPVtV+4C9wPZZDyZJi2i0KCS5BHiyqh457KbTgCenru8f1iRJ62xdTx8luRvYPL0E\nFPAx4Gomp45+KA9d97W+fOqOs3jNjrN+2HcpSRvK0/c/zoH7H1/VthnjVH2SnwHuAf6HSSi2MDki\n2A5cCVBVnxi2vRO4tqr+fpn3UzvrxlmNLUkbwhfyPqoqy902yumjqvpuVZ1aVadX1TbgKeBNVfUM\ncAfwG0lOTLINOAP41hhzStKiGfXVR1OKyREDVbUnyW3AHuAQcJWvPJKk2ZiLKFTV6Yddvx64fqRx\nJGlhjf2SVEnSHDEKkqRmFCRJzShIkppRkCQ1oyBJakZBktSMgiSpGQVJUjMKkqRmFCRJzShIkppR\nkCQ1oyBJakZhDTy9yl9ztwjcF0vcF0vcF0vmfV8YhTWw2t99ugjcF0vcF0vcF0vmfV8YBUlSMwqS\npJbj+dcfJzl+h5ekEVVVlls/rqMgSVpbnj6SJDWjIElqRuElSPLrSb6b5Lkkbz7stl1J9iZ5NMkF\nU+tvTvKdJI8n+ZPZT73+klyb5Kkk3x7eLpq6bdn9spEluSjJ94bP+UfGnmfWkuxL8o9JHkryrWHt\n5CS7kzyW5K4krxx7zvWQ5OYkB5N8Z2ptxfs+j48Po/DSPAL8GvDA9GKSc4DLgXOAi4EbkrzwJM6N\nwG9X1VnAWUkunOG8s/THVfXm4e1OOOp+2ZCSnAB8DrgQOA94Z5I3jDvVzD0P7KiqN1XV9mHto8A9\nVXU2cC+wa7Tp1tefM/ncT1v2vic5lzl8fBiFl6CqHquqvcDhn7hLgVur6tmq2gfsBbYnORV4RVU9\nOGz3F8BlMxt4tpb7Yl52v8x0qtnbDuytqieq6hBwK5P9sEjCD/7bcilwy3D5Fjbo46Cq/g74z8OW\nV7rvlzCHjw+jsDZOA56cur5/WDsNeGpq/alhbSP6QJKHk9w0dXi80n7ZyA6/zxv5c76SAu5O8mCS\n9w5rm6vqIEBVHQBOGW262Ttlhfs+l4+PTWMPMG+S3A1snl5i8kV+TVV9dZypxnek/QLcAHy8qirJ\nHwGfBt77g+9FC+ItVfV0kp8Cdid5jMnXyrRFfi38XN93o3CYqjr/GP7afuC1U9e3DGsrrR93XsJ+\n+TPghXhumPv/EuwHXjd1fRHu84tU1dPDn/+a5G+ZnBI5mGRzVR0cTqs+M+qQs7XSfZ/Lx4enj47d\n9Dn0O4ArkpyYZBtwBvCt4VDxv5JsH55Aejdw+wizrqvhC/0F7wC+O1xedr/Mer4ZexA4I8nrk5wI\nXMFkPyyEJD+e5KTh8suBC5i8QOMOYOew2XvYgI+DKeEH/33YOVyevu9z+fjwSOElSHIZ8Fng1cDX\nkjxcVRdX1Z4ktwF7gEPAVbX0X8XfD3wB+DHg6y+8MmeD+VSSNzJ51ck+4HcAjrJfNqSqei7JB4Dd\nTL7purmqHh15rFnaDPzN8CNoNgFfqqrdSf4BuC3JlcATTF51s+Ek+UtgB/CTSf4FuBb4BPBXh9/3\neX18+GMuJEnN00eSpGYUJEnNKEiSmlGQJDWjIElqRkGS1IyCJKkZBUlSMwrSGkrys8MvmDkxycuH\nX8p07thzSavl/2iW1liSjwMvG96erKpPjjyStGpGQVpjSX6UyQ/G+1/gF+bh59lIq+XpI2ntvRo4\nCXgFkx+EKB03PFKQ1liS24EvA9uAn66qD448krRq/uhsaQ0leRfwf1V1a5ITgG8k2VFV9488mrQq\nHilIkprPKUiSmlGQJDWjIElqRkGS1IyCJKkZBUlSMwqSpGYUJEnt/wHfU1/O9mZv7gAAAABJRU5E\nrkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "mesh.plotImage(v)" ] }, { "cell_type": "code", "execution_count": 16, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAD7CAYAAABqvuNzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFLtJREFUeJzt3X+MXWWdx/H3d2jL0NYWENpiCwXllxAsRYKsxjDrJiuw\nG9lE4yIx7mqySzaS3eiGiIkJNfuPxKybGE2AlSCsrLJrdOlGUTQ6S9wVUoEi1CJFpNBiCy39Ae1M\nOy3P/nFu2+nM/XHuPWdm7jnzfiUT7o/nnvNlpnc+832ec+6JlBKSJA3MdAGSpP5gIEiSAANBktRg\nIEiSAANBktRgIEiSAJgz0wVMFBEeBytJXUopRdFt9GWHkFKq5Nett9464zVY/8zXYf3V/Kpy/WXp\ny0CQJE0/A0GSBBgIpRoaGprpEgqx/pll/TOr6vWXIcqcfypDRKR+q0mS+llEkOq6qCxJmn4GgiQJ\nMBAkSQ0GgiQJMBAkSQ0GgiQJMBAkSQ0GgiQJMBAkSQ0GgiQJMBAkSQ0GgiQJMBAkSQ0GgiQJKCkQ\nIuKuiNgeEb9uM+arEbEpItZHxKVl7FeSVJ6yOoS7gQ+2ejIirgHekVI6D7gRuL2k/UqSSlJKIKSU\nfgHsajPkOuDexthHgcURsbSMfUuSyjFdawjLgZfG3d/aeGzKHDw4lVuXpPqZM9MFNLNmzZqjt4eG\nhrq+1umhQ3DiiXD//fDRj5ZbmyTNtOHhYYaHh0vfbmnXVI6IlcB/p5Te1eS524Gfp5Tub9x/Brgq\npbS9ydjC11R+6il417vgU5+Cu+4qtClJ6nv9eE3laHw1sxb4BEBEXAnsbhYGZXnySTjtNHjuuana\ngyTVTylTRhHx78AQ8NaIeBG4FZgHpJTSnSmlH0bEtRHxHLAP+GQZ+21l61YYGoJf/nIq9yJJ9VJK\nIKSUbsgx5qYy9pXHtm1w+eXw/e/D4cNwwgnTtWdJqq5anqm8bRusWAGnnAI7d850NZJUDbUMhO3b\nYelSOP10eOWVma5GkqqhloGwbRssWwZLlsCrr850NZJUDbUMhPEdgoEgSfnULhBSgt274eSTs6/d\nu2e6IkmqhtoFwhtvwOAgzJ0LixfDnj0zXZEkVUPtAmHPniwIwECQpG4YCJIkoIaBsHfv8YGwd+/M\n1iNJVVG7QNizBxYtym4vWmSHIEl51TIQnDKSpO4ZCJIkoIaB4BqCJPWmdoEwfg3BDkGS8qtdIOzb\nBwsXZrcXLMjuS5I6q2UgLFiQ3T7pJDhwAN58c2ZrkqQqqHUgRGShsH//zNYkSVVQ60AAp40kKa9a\nBsL8+cfuz59vIEhSHrUMBDsESepe7QJh//7JgeAagiR1VrtAsEOQpN4YCJIkwECQJDXUMhA8ykiS\nulerQBgbg5Rg3rxjj7moLEn51CoQjkwXRRx7zCkjScqnloEwnoEgSfkYCJIkwECQJDXULhDGH2EE\nHmUkSXnVKhAmfmwF2CFIUl61CoSRkez6B+PNnw+jozNTjyRVSe0DYXAwe1yS1F6tAmF0dHIgnHSS\ngSBJedQqEEZGso5gvMFBp4wkKY9aBYIdgiT1rlaB0KxDMBAkKZ9aBUKzDsEpI0nKp1aBYIcgSb2r\nVSC06hBGRrKPxZYktVarQGjWIcyZAwMDcOjQzNQkSVVRq0Bo1iGA00aSlEetAqFZhwCerSxJedQq\nENp1CB5pJEnt1SoQWnUIThlJUme1CoRWHYJTRpLUWa0CoV2H4JSRJLVXq0DwKCNJ6l2tAsGjjCSp\nd7UKBI8ykqTe1SoQ7BAkqXe1CgTXECSpd6UEQkRcHRHPRMSzEfG5Js9fFRG7I+LxxtcXytjveIcP\nw9gYzJs3+TmnjCSpszlFNxARA8DXgD8BXgbWRcQDKaVnJgx9OKX0oaL7a+XAgWxqKGLyc04ZSVJn\nZXQIVwCbUkqbU0pjwHeA65qMa/Krujyt1g/AKSNJyqOMQFgOvDTu/pbGYxP9UUSsj4gfRMRFJez3\nOK3WD8ApI0nKo/CUUU6PAWellPZHxDXAfwHntxq8Zs2ao7eHhoYYGhrquIN2HcLgIGzf3lW9ktS3\nhoeHGR4eLn27kQpeSiwirgTWpJSubty/BUgppdvavOb3wLtTSq81eS71UtPTT8P112f/nej22+GJ\nJ+COO7rerCT1vYggpVR4Wr6MKaN1wLkRsTIi5gHXA2vHD4iIpeNuX0EWRJPCoIhOawhOGUlSe4Wn\njFJKhyPiJuAhsoC5K6W0MSJuzJ5OdwIfiYi/A8aAEeAvi+53onZrCB5lJEmdlbKGkFL6EXDBhMfu\nGHf768DXy9hXKx5lJEnF1OZM5U4dglNGktRebQLBDkGSiqlNIHQ6D8FAkKT2ahMInc5DMBAkqb3a\nBIJnKktSMbUJhE4dgoEgSe3VKhBcQ5Ck3tUmEEZH7RAkqYjaBEK7DuFIIBT82CZJqrXaBEK7DmFg\nILuS2oED01uTJFVJbQKhXYcAriNIUie1CYR2HQK4jiBJndQmEOwQJKmY2gRCpw7Bk9Mkqb3aBEKn\nDsGPr5Ck9moTCHYIklRMrQLBDkGSelebQGj3WUZghyBJndQmEOwQJKmY2gSCHYIkFVObQMhzYpod\ngiS1VotASAkOHrRDkKQiahEIo6PZh9dFtB5jhyBJ7dUmENotKIMdgiR1UptAaDddBHYIktRJLQKh\n0xFGYIcgSZ3UIhDyTBnZIUhSe7UIBDsESSquFoFghyBJxdUmEPJ0CAaCJLVWi0DIM2XkJTQlqb1a\nBELe8xDsECSptdoEgh2CJBVTi0DodPlMsEOQpE5qEQh2CJJUXC0CIe95CHYIktRaLQIh73kIdgiS\n1FptAsEOQZKKqUUg5D0P4cCB7GI6kqTJahEIeaaMBgZg7twsFCRJk9UmEDp1COA6giS1U4tAyHMe\nAriOIEnt1CIQ7BAkqbhZFQh2CJLUWi0CIe+UkR2CJLVWi0CwQ5Ck4moRCHYIklRcLQLBDkGSiptV\ngWCHIEmt1SIQPA9BkoqrRSDYIUhScbUJBDsESSqm8oFw6FD2CaZz5nQea4cgSa2VEggRcXVEPBMR\nz0bE51qM+WpEbIqI9RFxaRn7hfzTRWCHIEntFA6EiBgAvgZ8ELgY+FhEXDhhzDXAO1JK5wE3ArcX\n3e8ReReUIRtnhyBJzZXRIVwBbEopbU4pjQHfAa6bMOY64F6AlNKjwOKIWFrCvrvqEAYH7RAkqZUc\nM+8dLQdeGnd/C1lItBuztfHY9qI7z3O1tCPK7BBSgrGxbHujo8euxDb+imytbktSXkuW5FsjLcM0\n7aY7a9asOXp7aGiIoaGhlmNHRmDBgnzbLdohbNgAt98Ow8Pwu99lgTA4mH0NjOu1IjrflqQ81q2D\n5cuPf2x4eJjh4eHS91VGIGwFzhp3f0XjsYljzuww5qjxgdDJ/v0wf36+sUU6hC9/Ofu66Sa4+244\n/3xYtKi3bUlSERP/UP7iF79YynbLCIR1wLkRsRL4A3A98LEJY9YCnwbuj4grgd0ppcLTRdBdIPTa\nIdx9N3zjG/D447BiRfevl6QqKBwIKaXDEXET8BDZIvVdKaWNEXFj9nS6M6X0w4i4NiKeA/YBnyy6\n3yP27ZvaDmH7drj5Znj4YcNAUr2VsoaQUvoRcMGEx+6YcP+mvNvr5q/+qe4QbrsNPv5xuOii7l4n\nSVXTl2cqf/vb+cdO5RrCyAjccw985jP5XyNJVdWXgXDfffnHTmWH8MADcPnlsHJl/tdIUlX1ZSA8\n9hjs3Jlv7P79+Q877bZDWLsWPvKR/OMlqcr6MhA+8AF48MF8Y6eqQzh8GH78Y7jmmnzjJanq+jYQ\nHn4439hu1xDyBsITT8AZZ3hkkaTZoy8D4f3vzx8I3R52un9/vrGPPALvfW++sZJUB30ZCJdckh3/\n/8orncd22yGMjWXTQZ088ghceWW+7UpSHfRlIJxwAqxenU3bdNJNIERkY/ft6zzWQJA02/RlIACs\nWgVPPtl5XDeBANkRSZ0CYccOePVVuPDC9uMkqU5qEQh5DzuFbOwbb7Qf8/TT2bTVQN9+dySpfH37\nK2/VKli/vvO4bjuEhQs7dwgbNvhRFZJmn74NhIsuguef73wiWTdHGUG+KaMNG+Dii/NvU5LqoG8D\n4cQT4eyzYdOm9uN6WUPoNGVkIEiajfo2ECC7CE3ZgdBpyiglA0HS7NT3gfDss+3HlH2U0c6dcOgQ\nLFuWf5uSVAd9HQjnnZevQ+j2KKN2gfDCC3DOOV7/WNLs09eB0KlDGBvLpnjmzs2/zU5rCC+8kK1d\nSNJs0/eB0K5DGBnpbroIOq8hGAiSZqu+DoQzzsj+mt+7t/nz+/Zln0/UjU5TRps3GwiSZqe+DoSI\n7Gplmzc3f/6NN+Atb+lum3mmjLxCmqTZqK8DAdoHwt693QeCU0aS1FzfB8LZZ2e/pJt5/XVYtKi7\n7bWbMkrJQJA0e/V9ILTrEF5/vdwpo127so/ePvnk7rYpSXVQ6UAoe8rI9QNJs1nfB8J0Thk5XSRp\nNuv7QJiKKSMDQZIm6/tAWLo0mxrav3/yc71OGbVaQ/AcBEmzWd8HwsAAnHkmvPji5Od6mTJatKj1\niW6uIUiazfo+ECD7Jd1sHaGXKaOFC7Nu49Chyc85ZSRpNqtEIJx9dvN1hF6mjAYGYPFi2LNn8nNO\nGUmazSoRCK0WlnuZMoLsPIPdu49/bPduePNNOOWU3mqUpKqrTCA0mzLqpUOArEOYGAhH1g+8DoKk\n2aoygdCqQ+glEE4+efKUkesHkma7SgfC3r3lTRm5fiBptqtEICxfDq+8AgcPHv/4a6/Bqad2v71m\ngWCHIGm2q0QgzJkDb3sbbNly7LGDB2F0tPcpo1ZrCJI0W1UiEGDywvKR7qCXReBWi8p2CJJms0oF\nwvh1hF6ni8AOQZKaqXQgvPWtvW1r4lFGu3ZlZy6fdlqxGiWpyiobCDt39t4hnHoq7Nhx7P7vfw9v\nf7vnIEia3SobCEWmjJYtg+3bj90/EgiSNJtVKhDGLyrv2NH7lNHSpccHwvPPwznnFCpPkiqvMoFw\n1lmwdSscPpzdf/nl7FDUXhwJhJSy+3YIklShQBgczDqCrVuz+1u3Zies9bqt+fOzxWSwQ5AkqFAg\nAFx4IWzcmN0uEghw/DqCgSBJFQuEiy+G3/wmu71lS7FAOHLm8+govPQSnHtuOTVKUlXNmekCunHR\nRfD449l1C7Zt630NAeC882DTJjj99Gz9YN688uqUpCqqXIewYUN2+OmSJdlaQK8uuAB++1t4+mm4\n5JLyapSkqqpUh7BqFTz1FPzqV1k4FHH++fDQQ9nt1auL1yZJVVepDmHx4uyv+S99CS6/vNi2Vq+G\ndevgwQfhqqvKqU+SqqxSgQBwww3ZOsKHP1xsO8uXZ+c2jIzAFVeUU5skVVmkI2dn9YmISO1qSik7\nKuiss4rva9eubIG61zOeJakfRAQppcKfxla5QJAkHa+sQCi0qBwRpwD3AyuBF4CPppT2NBn3ArAH\neBMYSyk5SSNJfaboGsItwE9TShcAPwM+32Lcm8BQSmm1YSBJ/aloIFwH3NO4fQ/wFy3GRQn7kiRN\noaK/pJeklLYDpJS2AUtajEvATyJiXUT8TcF9SpKmQMc1hIj4CbB0/ENkv+C/0GR4q9Xg96WU/hAR\np5MFw8aU0i9a7XPNmjVHbw8NDTE0NNSpTEmaNYaHhxkeHi59u4WOMoqIjWRrA9sjYhnw85TSOzu8\n5lbg9ZTSV1o871FGktSFso4yKjpltBb468btvwIemDggIuZHxMLG7QXAnwJPF9yvJKlkRTuEU4H/\nAM4ENpMddro7Is4A/jWl9OcRcQ7wfbLppDnAfSmlL7XZph2CJHXBE9MkSUD/TBlJkmrCQJAkAQaC\nJKnBQJAkAQZCqabiRJHpZP0zy/pnVtXrL4OBUKKq/4Oy/pll/TOr6vWXwUCQJAEGgiSpoS9PTJvp\nGiSpamp5prIkaWY4ZSRJAgwESVLDtAVCRFwdEc9ExLMR8bkWY74aEZsiYn1EXNrNa6dar/VHxIqI\n+FlEbIiIpyLi76e38qO19fz9bzw3EBGPR8Ta6an4uH0X+bezOCL+MyI2Nn4G75m+yo/WUKT+zzfq\n/nVE3BcR86av8qM1tK0/Ii6IiP+LiNGI+Gw3r50OvdZflfduu+9/4/n8792U0pR/kQXPc8BKYC6w\nHrhwwphrgB80br8HeCTva/u8/mXApY3bC4HfVqn+cc9/BvgWsLZKtQPfBD7ZuD0HWFSV+huveR6Y\n17h/P/CJPqz/NODdwD8Bn+3mtX1ef1Xeu03rH/d87vfudHUIVwCbUkqbU0pjwHeA6yaMuQ64FyCl\n9CiwOCKW5nztVOu5/pTStpTS+sbjbwAbgeXTVzpQ7PtPRKwArgW+MX0lH9Vz7RGxCHh/SunuxnOH\nUkp7p7F2KPa93wscBBZExBxgPvDytFWe6Vh/SmlHSukx4FC3r50GPddflfdum+9/1+/d6QqE5cBL\n4+5vYfI3ttWYPK+dar3Uv3XimIg4G7gUeLT0CtsrWv+/ADfT+prZU6lI7ecAOyLi7kbLfGdEnDSl\n1U7Wc/0ppV3APwMvNh7bnVL66RTW2kyR919V3rsd9fl7t52u3rv9vKhc+JjaftK4jOh3gX9o/LVR\nCRHxZ8D2xl9KQbV+LnOAy4Cvp5QuA/YDt8xsSflFxNvJ2v2VwNuAhRFxw8xWNfvMpvfudAXCVuCs\ncfdXNB6bOObMJmPyvHaqFamfRrv/XeDfUkqTrjs9DYrU/z7gQxHxPPBt4I8j4t4prHWiIrVvAV5K\nKf2q8fh3yQJiOhWp/3Lgf1NKr6WUDgPfA947hbU2U+T9V5X3bksVee+20v17d5oWRk7g2MLIPLKF\nkXdOGHMtxxbWruTYwlrH1/Zz/Y379wJfmc6ay6x/3JirmP5F5aLf+/8Bzm/cvhW4rSr1A6uAp4BB\nsr/uvgl8ut/qHzf2VuAfe3ltP9bfeKzv37vt6h/3XK737nT+j11Ntkq/Cbil8diNwN+OG/O1xv/8\nk8Bl7V47Az+Ybutf3XjsfcDhxg/yCeBx4OoK1H9Zk21MeyCU8G9nFbCu8f3/HrC4YvXfDGwAfg3c\nA8ztt/qBpWTz3LuB18jWPBa2em1V6q/Ke7fd93/cNnK9d/3oCkkS0N+LypKkaWQgSJIAA0GS1GAg\nSJIAA0GS1GAgSJIAA0GS1GAgSJIA+H9bwH19wsu7WQAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(time, wave)" ] }, { "cell_type": "code", "execution_count": 17, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">> Start Computing Acoustic Wave\n", ">> dt: 8.00e-05 s\n", ">> Optimal dt: 1.09e-04 s\n", ">> Main frequency, fmain: 1.00e+02 Hz\n", ">> Cell per wavelength (G): 1.60e+01\n", " Src at (-124.75, -124.75): 1/ 1\n", ">>Elapsed time: 4.35e+01 s\n" ] } ], "source": [ "# Step8: run forward\n", "U = prob.fields(v)" ] }, { "cell_type": "code", "execution_count": 18, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# Step9: project data\n", "data = survey.projectFields(U)" ] }, { "cell_type": "code", "execution_count": 19, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0.0333333333333\n" ] } ], "source": [ "print 50./1500." ] }, { "cell_type": "code", "execution_count": 20, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 20, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeUAAADTCAYAAAC/ZAtmAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnWnMLclZ3/9P99nf9733zpJ7HebeIfZ4yYADgYQJS0QI\nyMZEwiaJZcyHEEJQSBxIPkQIJlhCiiKNRURAIiEfwhKiEIyJEmKiiMyMEiEBtscWiz0eA6NRxp65\nmcxkzMy7nr0rH6qru7q6qvc+p885z+/q6u3TS1V1d1X963lqaRJCgGEYhmGY7eNtOwEMwzAMw0hY\nlBmGYRimI7AoMwzDMExHYFFmGIZhmI7AoswwDMMwHYFFmWEYhmE6QuuiTETvIqI/JKI/JqIfbjs+\nhmEYhtlVqM15ykTkAfhjAN8C4P8A+CSA9wsh/rC1SBmGYRhmR2nbUn4EwLNCiM8LIZYAPgzgPS3H\nyTAMwzA7Sdui/ACAF7TfL4b7GIZhGIYx4IFeDMMwDNMRei2HfxfAg9rv2+G+CCLixbcZhmGYg0MI\nQea+tkX5kwDeTERfCuAlAO8H8F3mSe997z+qGU2Tui6fEaUeVcxnP/txfPmXf236SsdFrv355+j7\nBISI/+vXeJ4X/vfh+x7k+DqB9XqN1WqFxWKB+XyOxWKB5XKJIAhARPB9H6PRCKPRCMPhEMPhEP1+\nH71eD0QEIg+eR/A8L9oWQiAIAiyXS0ynU1xdTXF5eYGrqytcXV1huVxivV4DAHy/h8Ggnwh7MBig\n1+tFaXbd96c//Vv4yq/8y6kj8bjE5LMIggBBEEAIAd/3o7Dn8zkuLi5wcXEJ3/dw33334fbt27h5\n8xZGoyEuL6/w0ksv4fnnn8fzz/9vvPzyy1gul7jnnnvw0EMP4W1vexvu3HkQw+EAp6eneOGFF/DS\nSy/h8vIS/X4fN27cg/vvvw8nJ9fg+z5WqyXOz89xfn6Oq6spVqtl9I5830e/34fv+/B9HwCBKHlP\n8X1upq369NMfw9vf/nUbiWtbuMqf+YzzfpcJX9/39NMfx9vf/rWZaQmv0mPX9qXzRZm05tU/Reqn\nrrOL+fhXfuWnrPtbFWUhxJqIfgDA45Cu8p8TQnwulYieXzeeUvuLHndfp1ekRcJL79fLARGlhFa/\nTi+IQsTXKuH0fS+q9GVYwHq9xnK5xGw2w3Q6xWw2w3y+gBBBJMYnJ9dwzz03cM899+DatesYjYYA\nCIvFHGdn53j99dfx+uuv4ezsDLPZDEEQoNfrR0I+Ho8wGg0xGo1w771SqFUDYLVaRUKpREjdp2oU\nSPGWlY7Z4FitVpjP5+E9UvRc4t9e+FeGZ4q8EAKr1Qrr9Rq+7+PoaIJ+X6ZdCGA6vcJ8Psd0OsVy\nuUCv52M0GuPo6Aiz2Ry9niwa8jkuQETRPenp8DyVluR7ZLZHnsjo78gUOdfr08uc+hvHU0zU9Dxu\nxmPWB3E80V4j/uQ96PneFXcyPrIe2weB3nXatpQhhPgNAG/LOsfzsru28yo6d0vYurdQPFkt0XR0\n5StiPXgZtgrUFlbcalaVvxIgJXaAFBAlRovFAovFArPZHMvlIhQnD77fx3g8xsnJCW7cuAf33nsP\nbty4gePjYwwGA6xWK5yfBxAiwOXlBV599VW8/PLLODs7w3K5RL/fx7Vr13Dffffh/vvvx2QywcnJ\nCQaDAYgotKD1hsAsSpcuyKpxoYRMHVMV1nodYLFYJO41FkIfvp9uzKiGSRAEkSALIdDr9SJvwGRy\nBCEELi8vwzikB6HX6+H4+BiLxRy93hUGgwGEEFgul5jP5wAQeRlUJRhb5K5Kzl0RMu2Q9YxdYiw3\ns8qw2SAkbZ89XpfI6/El0xOHI8uG3KHqRhk8pbb19AEiUSay6k2XELNAb5/WRbkIbb18Pdg4s7kF\nXE+H2dWtZ9abN+/A8/SWKVkaAGWF2t4gUAIg/yPhSvY8LxI05aqOBXmJ9XoViaASpfF4jOPj4/D/\nSWRByuukAJ2fn+H09BTn5xe4uLjA5eUlLi4uMJ9LCzIIAgwGAxwfH0MIgeFwiJOTaxgM+ggCgdls\nhsvLS5yfn+HiwsdsNsNisQjFUgqlElD17PVKTIgA9977BqxW60jwYkH24PsCgHT/KnyfNLcwEnEN\nBgMcHR1jMhmj1+thvV7j/PwqsugBoNfr4+TkGOv1Cp7nhxYwRd4GIQTm8wVWq5X2rsnpht8Vi/nm\nzdvbTkItilrFdkEG7OVUrwdMMZbHXWJs4+ZNc8JJtljqXRqq/Kq0yHPM+OO6La7HkgKth2+LU7Gr\nAr3r+VinE6JctgIrmkHKZSozDUp0zesF3vCGeOya27q2FYjsOG3PwRQkJcgAQtGRVuFyuYz+Kzer\n53no9/uhhThJCPFoNIpctJeXV1guF2Ef8RUuLi5wdnaGi4tzrFZr9Ho9jMfjSHw8z8N6rdzVsg/Z\n8wjD4Qi+72E4HKDX88PGQpy+2JKn6HhsbUqLwPME1mvg1q074fnquSjLQVnTsXdBNlb60TNTlvJq\ntUKv52M4HOHatZPQSg5wenoW9jNfQAgRueN7veuRyC6XSwBKlOcIAmlVSxFXjQovFPA4/bp7Um13\nuUK7efPOtpNQmjruaXNbCzUMOxmHaRmrc7IaYfGmfP/pZ5zlFbPfi9k4V8Kc7vayC3RW2K59pkBz\nPt4MOynK9v7XetjCUq1W40znb9tAcldlkLSw7WIsLbG4D1WxXq+j/0p8lstVwjJWg6pGoxGOjpRl\nfITxeILhcABAivpsNgsFSg5OUhaxGhTm+x6uX7+Oo6OjSFR9vxf2zQrMZlNcXFxiNBpFfbb9fh/j\n8QSr1TpMm3SlCyFdwFKYffT7a/T7fW1wmRyg5vs+1mvZqAiCIBJnQCAIYs+AQg0cU+lbLJahq3kN\noIfRaIijo2McHU0wm0l3+uXlJU5PTxMNF5WOIAhwdXUFz5P9yOqa+XwWNnpk3Low69ZN+IY1gRa5\nXTRMPmWt4nwhBnQxdgmxLdq0ezqrnJv3oOov1ZBLhOyMS22nGwpJ4VUWctKyLmc96/tt/dVdFuhd\npxOirAbPmJRpEbsplnlcUZWzsM3Ws+kSj89PVhyUKGyyAk9ep9y6SqSUe1rvOyUiDAYDDAaDyFU9\nHk8wmYwxGo0wGAzh+17k3p7Npjg/P8fZ2RnOzs5wfn6O6XSK1WoVCbs+clqJZNyv6mE2m+H09HUA\nsh/45OQYo9EInidd5sfHR1guF5FVvVotQyt0hSCI+5qJCL2eF4mX7/sQQrqalSt6vY5dx0EQW8/L\n5TK0bqUVO5vNw0aAQK/nR8+j3x9EfezSvX6Ofn+A69evYzgc4vj4BIPBMPRAIGrwyH5xOfBMpj3Z\np6y6EtbrIGEhJ/PCbrizu0IVb1gdq9jlnjbTYnN958Vr3ksyLpeVbW/k6QMjk2HZBVq/T/OcWLyz\n+6FtRhALdHt0QpTXa7so51E2L7gyjxpcYTsvnUldYSgRzk6fmpaWHl2ZdovpFUAQmNZxLFJKHPr9\nAUajYSjIk7AfeRBZoavVMhp1fHV1hcvLy7C/+BLT6VVkEQJyOlO/P8BkMsHRkbSw+305VUr1W89m\nM6xW62gQ2Hw+w3R6A9eunYTubjXK+wTrtbJ4pcAul4uogRE3ygj9fi+yfIHYKyCtZz+6Rgpg7FJT\nDZP1eo3ZbI7VahmF0+v14PteaKEvoz7v6XQaCajqIx+NhmEjYhEeR+R6l6IsG0EqbPnsZf+2mjIm\nRDxozZbXGDtV3NL6dp7FqcTWNYLadF2nLe84XL3hldfgMl2/evz63/T9mA13GX8cnRJTkej7jo0B\nu6VPhCg8071tWtJ5DRyb5c15vB6dEOUifSvWqwpdFhe2PDdN4qqMjGVr+aogLHPBLdfbw9IrCSVU\nUoQCBME67CtWI5gRWq8+er0e+v1e5K6W1u0oHPjkIQgElssp5vMFptOrSJCurq7CQVjS6lbWca/X\nw3AopzpNJpPI2lZ9xsrt7XleGMY8ml50dXWF6fRG6PI+DvujJ1rFJAvx1RUit7Y+6hoYJjwGvV4/\ntJpFJMbSkl+F7mmEoohoSpaykgeDQWThy2tn0XxqOTJcDlyTFnUPk8kY6/UgbLBcRX3zajDYfC7F\nWuUXOQrcD13YBOVRj/uUy+etQ6TsiGl9u64Qy33uhrjuns5yjzdD2oJPWrIiOkcfZyF/q3OTFnAc\nrt0IUNem+6STFHVxm+dyXi9PJ0S5Sy8uFldbZZo+ZrYUs6zxOI6keyzerywsQPafxpZk7KJWfag+\nej0/6jvu9wcYDPrR4hRyZHaA+XyJ+XyB2WwaiebV1TQU43lkGfd6vVDc+5ELfDAYYjiU4fp+L5qT\n2+97YfrixUmm0xlms1koylPM53PcuLHEyckJhsMBJpMjqEFqal71xYUU5tVqhel0FopZPJBLCWo8\nwMyH561ABHgeIQh60XNUC5ooq1beh3RZAwit+hXOzs5xdXUVuaJV40f1LUtxnmA8HkfTpZS1rtIa\nu9v90AqX1rIMS2gVZnJeuXrHLnfloVDXIpbb0ZYthpRrOo43LXxm+Da3cVUhtjXg9fiSx7MbcNn9\nwXHa1S5ToGPjIdk1lsyfcf+ziywXt+0ebfsZNwcpysULVtwyLRpWlotHd1kri1ovVLHrU4pTbEGK\nyHJULlNlBfb7A/T7vchqVm501W88n8/CBURm0bxhNUJbCDmoqtfrJYR4MEiuvCUt1SDsu42naKg+\nZyWGauqT6j9eLKSI3bhxPbS6x2H/qw+1+Mfl5SXm8znW6xVms3hlLnXvvV4/TAdBerSlVaymPknL\nOwgXRokHVckR52OMRkMIIXB5eYXp9Aqvv/46Li4uowFn6h2obbkKmeyDHgwGiS4DdT+yceFFrnHl\nHl+vPSi3YBymad3E+US5EeO8tr+UHR9i76dNnOGMo4oQJ8NPCrEtfVnY7tV1/0Vdv3bvXHmBVtfG\nLnrz2UQxJNzbtvhs6df32dzyLMz5dEKUmx6ZmlWA8jJG8lC5iiQvTv2vEEF0LO4rii1lddwUY2XN\nSquuj14vXlJS9jsHoYBIl7Jy0yrLULeMlfAqcR8OB5HIK8tYIUccrwEsI3efcvuORqPIoleLhqiR\n19IFvMaNGzcid7gaHCUXM/FxcXERrhi2xmKRnFI0GEjrU2ZVmZ54QQVZuagpWQCiZ3R8fIyTk5Mw\nbXIhFDn3Wq1OJjSPglpwZYnBAJEbv9/vY7FYAFCWuJzaJaeSUfTsYmt5VUh8ssSleZfodinqmjZ/\np8XSfC5poc0SYj0tbbinsyziMuRZlq548hovcVcKWQVat6BNgU6Og9HrsOyVxGz3wpZzPp0Q5aI0\nUWEVzQiulqGeluwBYfa41Hlq9LBpHeutTLWEploUQ1llsp9VHgMQjcaWIhj3qyprebVaRgtpqBXA\ndHGPBb4XuWR1C04K7jqyYAEljNJNKFe/igu4Wrry7Ow0vE7e340bN6I50vJaNb9XxiXnA69DKxbQ\nXdnyeXuh6zqeL61EWR0fDAYYjye4du06Tk5O4HkeLi8vcXV1hbOzs6ivWDZmpPgCCEejSxe6nI5F\nUX82gOjZLpfL8D6h9efL85bLVejiz89fpjWxT2JcRojNfdlibLpds/uJbWlJuqjr9xM3JcY2igq0\n7p5Wv+P7Md3jcYM37cGzC3Qszvp+Ke5lrGfXu3Dd36HSCVF2TYmqS96LrnLcVZkmW+l2YoEBVJ+x\n/lsJsRIeaR3H7lEpznE/quwTFdpc5YW2iIiaLrUKBRORyCSt7V6iHzqeD60aD0HUaFDWsJoLrX8I\nYzgchAXYiwaALRYLnJ+fR54Baf0LTCYTjEbDhCWunp0aAa4+5KDP5VReAbWt0qXuXzYQhqGVfBzO\nlZbH1ACvxWIO3/cwGo1BhHD6loflcoGrq8tolDuA6HkDCOcry+lUauCZ78fuazXoS92/cs8nhYQy\nKybb712prKq6p9PbQJYQq7iyhNhMT7qvuJ572gzf9rtpsgTM1k8dC65+/8m+Yts7iZ+rEu20OMtz\nAdmwLu7ezroHl8v7EOmEKG+aIi++aCEtWjiTmTPZCFFiplvHuss6/vpT7DZS/cJxf2dsKauVvvSv\nQfV6fvSVItMqVnHp6VPT1OLtZOFSrnZlSUphHkajkVXhXi6X4cpZyWuPjo4wGAxwcnJiPCdErmzV\n0FAVixD9KL0qHP3eiShcLGWCyWSCXq+H5XIRjQyfz+cIAhFO9ZIuejVlbLlc4upqivV6DSIvEl05\nlQrRlK/FYhGN5pb3nXxHqqGhLDk1h1l5KIr2Ldr2d63CqiPE5u8yLmqXe1pPk02I1f4uWsVlKCLQ\nNstUliMVBqA/EzN83Yo2xVn/a3N7F0l3Vr+zbf8hsTOiXPYlZVV+RVvtWWG6tm3hmhZy2qqOK2/z\nE4z6GtvKajVX8lL9tmpQkgpbWtlq2k7s/lbbygoE9ClY0sWlXML6fFvlolaLiKjrVMH1PC/6MIVs\nXMj+4sVigcvLC01cZZxHR0cYDgcQ4hjSegmi/1JA11GfrrpGF2Q14nqxmEMIYDDoa0tm9hAEa21E\n+CzqDx6PRzg6Oor65Pv9PtZr2R8eBHLktkq/aqyor1apqWCyf17v908uh6r+Ky+F7q5XjRZXfqva\nIGyTOi5p++9oyxlP2j3tXnM62S0Uh1vVPV2kP7cLuLrNbM8xKbYAQIl6SQ9PBatWKTQbQzDG26jw\nkseSFrTNmnZ5hlzvqYvvoGk6IcptLUFYpbJwYXNPZ7kiTUF2hZn8wIQfibOqgFShUX26+mhgfblN\nFY8+hUiJsRLh2JpTlnFckZmudCXuevrV/ar3FVeEseWvRoYr61xOfbrEfC7FMb53+XcyOQpd39dS\nz0qOyl5H845V+NKFT+EiJNKFr1YQG4/H6PcHWK8DzOezaKWy2WyK9TrAcOhHX7YaDkfh1CrpCZjP\nZxBCNjLU6HZp/Xvh2tdyQRHf9xOfp1TPzbRgkpayF1nSqotCCbZeWWaR13XSBmUs+2LlLdrSY0mJ\nbfKvXYzTlrYMt6oQR6kp0ODuMi5LWn+myXOAtEDr4pycA627rMMQNBHXLWb3+tu2vGu+L1f34S6+\nkzJ0QpTrUsT1l3WeDVdmNs6Cyy1md8fF4emfI9SFOY5DJMRytVonrGHdklVhKteuLsiqX1p3Ketu\n8FhQ1DSkON36OQpz5LNutao+JtUIiAd0SXGSbuQZiF5PPEtpMQ9xcnKiWToyfiV+QiC0UHvw/RXU\nhzhUX7Ic2X2E8XgcLiUqvwl9enoaju6Wn19UC6McHR1hNBoDkHOl1TxmIRCujOZjMOhjvV5H1r9q\nBKhpX6qrQHoL1mG/e5DwLihR0b0guos7dv/FzzxPoF3irO+rQx0Xe9Z52hkoIsRyO3mOXeSbdU3v\nU6XvErgsgVY/0wJtX0FMF2h9YJku0OYIbjN9tjTm3VOZa3aJTohy2QKUVTm4gsp7d/aKITpqiw1m\nqzxbjNPfQlbbegVlLh6SXEAkgD5/GYhHQSfd36oPU3eVUsIVnraIpTAn7yN+dnqaXH2jsh84tmil\nMCddu3JVsWmigeB5Po6OJhgOh7h+/VqUTnn/AovFPNXHTERRg2QwiJcDjT+UMQs/QXmGq6srBMFa\nWxxlGK5YNgwbFOuoASCf8zj6ulQQiMQXteKPbCzDpU7X0Qcq1uuk90K9JzVKXt6v+qu8DbqHggAk\nxxuUEWfz/KIVVt55TTZ6dTdn2j0tj9msYj18u7vV7ZFyp2U3XNRNkJUv3PcsG0/x9bqYulYQcwu0\n3Aeohrrp3q7i+dlHge6EKJdf+7po4csudLaWufztiNXRKjf7scw0EMXTeExx1jO9sjj1xUP0hTT0\nLyMpkVUVvHJ/633R6lhcIGKrWKY7yKjUkgVQxamegXKV689Rpjt2eythVi5ade5sJtfejp+HTON4\nPI6sa+UlUMKrXMVyhS3lHUC04MnR0REmkzF838d0OsXFxQVOT89weXkRDc5SYiznYUtPgnp3anEQ\nuazmEkReNKJ8PB5jOByi15OLs+hLn6r/q1U8ClwKt96/v4YQfuIZqHWy5aAyXZxVP3Tx0cF5/XKu\nyiqrEisqxHlpc8VZTYiBpt3Ttt/7Tp57W/9tf9/RVph3lefDJdDKjZ2sW9ODx7I/jFH0vnb9fXZC\nlIuLrImrssl21xQRYCDd/2W6fLPTnRZfPV1xv2y8iIgSYzXaWQm02lbpT1qZep+0h2Slp38aTjjT\nb7rBVfqTz8/2fJKLByjxVa51FZ7vS6tTv24+n+Py8lILX96P+pCGGvAlLXMP0+lVNLBNDf6Sfb4j\nrS+5jyAQmE5n0feS5dxjhPOXR4kvXqmKQ02rUq5ytdqX+pCHFOVRNEhOiWp87/EzlNay+q+eqRd9\nVCNufMWf1pN/PRCp96AquWQlVdQaLFI5ufrr8sIser4tLpcQm+mxCbHaX6fS3lcXdR1secXlQYjP\nVe9CPy9e710d1wVabgN6t0RsLKTr4zoCvevi3BFRtmUCx5k5LTszzKz3kuceyxIxd1h2MVbnqMym\nRFAPW7eKVcWtMPtykxZyWvRl2MkvFelx5VlhKixXBtcrSD1NNqteTZMaDkeRgAkhsFgscHFxGQm6\nGgSlRlCrbyerFcAuLy+xWFyFn0/UXeTxGtfz+Sz6NrQScrme9RFOTk4wmUyiaVzxsw6ir0B5nofZ\nbB5NjRoM/Gjp0Xg+t/qv5iMnn6+6P/WRj/Ua4TKcATxPTaUSCU+JrMgoes/pufvxYg1lXbV64ynr\nnWddX+VaXYT1xnAxMbZbw2wVt4OrnCcb3emGevxe7O9RF2gVnk2g42vdBlNVcbbdV5fphCj3en7G\n0eKC7QzBYRXEIgXoIgwkV9xS52eHq/cXJt3GZsVjimMswPZK1xTiMvHIbRmXfsw8x9XYKVKp6YVO\nTZfShVm5suWiHUPobq/ZTFq1sv9bjqqOl+Qch/cqp0HJucRXWK3WkbtezjPuhwPJFri4uMD5+Tku\nLi4wn89BFH/XWX6acQTf7yWeOSC7UNT3mBeLeWgtByDqRYuEyKVHvXCamRc1FuJGkXoXsQBJi1h6\nP4hWUSWlXOfJdxlflxbnpPVse4dZ5Off7GuqCzGQZRXrYZsNxqasYttvxk4R17YpsvIaW3dXeoyP\nbjkr97ZNoJP7KBG+mc4699VFOiHKVQuQrc/D3Q+StnZdlZurMrJl0mT/raxUbSKpsC2v6RJ/3UWa\nJcZmXDb3X9a9uZ63eb+2hoIa+KWfo+YSq0Fq+jW+74f9xtKSDwI5J/ni4lyPGTduECaTMY6OjkFE\nWK2kIKs1ueXUq37oiiYsFktMp1Ocnp7i9PQ0nHMcoN/vYzwe4fj4GJPJEQaDfphuOcBLR3UdxAuz\nBOj11DtWC7D0jEaSfOdSnON+fSXayv0OxN0Q8fOKV0YzPR0qbH36mfn8pSWS9nrkiZlu/RQh7zyX\ne1ruS55jhulqFFexiE1rLs99X0XsD4kyAp18nuqv/knbuK5R7zn5ztIriCXTEgu5Ho6ZzqL31WVh\n7oQox1/8ycZ06dmW59QHXZmCa2uJpyuHJLaKUt9OVs7VXeJxfElLyWaBq3S544nDzcqw6cLkxtXQ\nsTUG9OehD9YCEPXLSvd07KpfLJahxSzvTVmno5EcnDUey9HVcoCY/JrTeDxGr9cLB2md4eLiAq+9\n9houLy+xWq3g+z7G4zGOjo5wdCSnSxEhWn6UKP4Ws7oNNe1MnRMEsoioBpISW/Xc4jXL475mNfpd\nLkqyTgzyk3EIAEHiGaW9IemuAVN4VB6KrY1iFVbW+7YdstVf6TJRTISTcdg9Q2ozr94048ha78DM\n5yzI5XB5Wmxirec9ilzc6hgQ103quqRxoOcn01Ma5/XYgq7ybvMarNukE6Ks969mYRMc28vIcoOZ\nYSSJKxZbpSMrcU97cel+3Ky0FBFj87fNtamfmn5k+VZSliegqqvQVWh1q1ntUyOfx+NxlAYhLrBc\nrnB+fgHfj9fmVmIov8k8wfXr16P1qkejEYgIV1dTXF1dhlay/AoUAIxGQ0wmExwfy3Wwh8NBOK94\nGVnJqj0Yi2FsResNCplu9TUoOUI6HtgVf3RDnRevhW3OQzf7oNP9bboFraxhs/GnP/e48gOgzXm2\nveN4n+1Nut59Mk6z8nQJsR63y2tjF2s3rnKSdW5RFz9TjiKWtPqdtKKTAq1/xlaJeJzfEyGlGmwu\ny9pMX9l72qY4d0KUV6uV85irha2OuQQmW3wVzYiwzTIxxThOb3aabNa4ntbkvaXvMc8y0sPMK1RF\n3TxmnLqLVu9j1t3ZapDWZDKJpj8pC/fi4jz6RjRAGI2GIPIwmRzB93tRfhFCYLVaRhbya6+9Fgmy\nnF41xmRyFH4AY4Rez4/meut5Trncfd/DahVbr3JaUxAdV/OcVV+yELLfXH5n2o9c3oCIXOzxF62S\nq7Xp65+beUe9I73PWZ/brp9jvitTnHXhTr6nYhVWkcZhvlVs80wlz5FByAo7yzK3xWc7l4V4s9jE\nzPae3AIdb6vfujjLa/X4kvuyRnDvmjjvlCgn9xcJOX7R+rM1Kxj9mOlqVMKclT6b6Lpd5GbC05aG\nS4zT95ZMh+23q8I0ra4oNSVE2BR3/VrTAlTPRu9jVmtxTyZjqP7l6VSurHV+fg7flyKqvos8mYxx\ncnICImC5XOHq6hJf/OIXcXp6ildffRWnp2cQIoimR43H8ho5ncmHEkPlngZUo0FEwqfetfpS1WrV\nhxBB1IetVvfyPN8Q5SB0d6+t7mjd7S2fRyy6+qIt6tnJ5xc/Q7WUaGzR20bDJiumOEvYRm6nLWpX\nXnCVF/N9x9vRltMqzorPfU8sxl2nqYa+Hp6tzsoS4/SqYvWnV+n30DadEGX7pxuLPYCsVnX8+TwZ\nXpZFbBdEW1pcfdXxsXzL2G5pmH/zBLlo34mr8swjS3SLoq6NLWYlzASidWQxj8fjyCKdTqfRQC21\nnOV9992Hk5NjXLt2Hb2eXCBkvV6FAn6B119/HZeXV2FYk8Qc5n5ffR9aLoMZr5CmFvhQlrJcwpNI\nWcor7cOikRjrAAAgAElEQVQTFC7RKQepyTnMaunNFTxP/jXnaMu/cb7SB3V5nsovpPVNm3lLhSXP\nk2EG1sVhzPemvYXw/uW1ap8+Sja8UgsjbY2bZaUJIVZxZaWfXdS7iU3Q8uqQrHdoF+j0CmKxOAM2\no6xIXFlpa1ucOyHKpviUKaQud5rpagPM0czpYyZmpZOuNAG3VWx74W63nzvT2t3RroxRxH2dRZWM\nl/V+dBeULsxyFSx5jvqIBRCP6L64uMDZ2Tnm8zmEkJ969P0ejo+PMRwO4Hke/uRPXgu/czwN165e\nRp+kVMtoqvXAY0FeaW50Cucjy3ile7sXTlcKopXA9G81q5HjahS3coV7HmkfqCjyzBDdr7KK1dKF\nytK2Wc/q+SqBJgq0gWfuRmp8rQfXamFuqzTOt7Ic6O5Hs5xUE0cW4v3GVme1JdCmGOeJc15crnPb\nEudaokxEPw7g2wHMATwH4O8IIc7CY48C+F4AKwD/WAjxuCsctTawETbyLEr9PP352PqBXdebpEU3\nW4zNiiFLjIvdSyI1znDzfle1js0wbZ6DOha0DE9eu17LecHKCh2NZD5YLpeYTmfhnOQVhsMhptMp\npBvZ0wZbxVYvEYUjnnsYDAbRN6NVQ0AtfymtY9XHvQ5HV6s09CGEiFbeWiyWoSjLxUr6/V4kyvLz\njhS5r+X2SluRC4jzb7Jf1+aJkc8lOVVENQTVs067t02LNbnOeLr8JOPLw8w2ZoMz/ReA1SVu7yeW\nacrK//ZzWYh3kyzjoAmBjhu4cX4zXdlZi5PkxZN1bpMCXddSfhzAjwghAiL6EIBHATxKRF8G4H0A\nHgZwG8CTRPQW4bhjtdi/wuU+dlnQ5vQhM4w8UVHbpis6vepWWqDNcNLkNyxc5GUQ13HdMioaV5Vz\ny2AKuP7dZimwXtRvOxzKFbR6vV4kesvlEvO5/JpTEMhFR+QCH4i+SKWmTo3HY20+chBauhStTa0L\nhz5XWfY79yMhXyzmAES0rOdgMNBEeaAJvlwrOxZlaZF6nquvXeUn090Wi7M6V1nPNlFWv/V++iCI\n5zbbXObJMmN713rYSbG1paMKRStk234W5P3AVT/pdXPZ9+/y/pjubbnfLtCmh6ooTVrPtURZCPGk\n9vPjAP5muP1uAB8WQqwAPE9EzwJ4BMAnbOEUmaec51qT+xNps27H+xK/rNOy9Ja/3SqOj5dJcxnX\nswtbJnBl9KYFukwGTLpOkwM1lCtUffJQnavmIKsBWfJrTdKl/cUvfhG9Xg9XV1eYTqfwPMLx8TEA\nhF90ksKsrF35vWX5bvVPXuqNAylWQK/XDy3uHuTXqdQHKtaR+I/HsaWs3N9q2U9p/a+j+5NzoKUV\nbs+bZl+u6d1Rf/XFZuLnFH/IJZnX9YFmyY+GeGG6kp/xTKaLonnU8f4gkf9dZcueH9IzFmzn1bGU\nmN3G9m7zDKmshqGZP01DyHRxhzEk4tIbxlnpzIq7qkA32af8vQB+Odx+AMDHtGN3w31W+rR2HYpv\nUn8elmfjeR6CnK9NJR5YtC+83rLWLmnXWR8vxRVZEfec7XdlbK4XxM8qUdlCv+8G4i8bhH6+et4k\nIMI+IBILYLWCMvgG3hon4z5GvePoe8Z9WuPy9FWs5xdRP+98OsXQF7j32gQn49j13O/3MfAFsJoh\nEEuIyJUrQOs1+iQglBWLIHxeAfyegPClcPWwRLBYY7WS1rAv1hj3CQMvXkzEoxVoLRAEi+j2PAGQ\nCOB72V0KCYR5jNJ5nACQQIBQICHQ8+PvbdvC9MgDCQKF1rP67QkPBC8c+qVZEqoMkEAgguh3ALmt\n/4ujyctbFMXjegaFvEYFp3Axe0L4ulN5Q3mVohMKeCy1NjARRX+zjLu4R8ie1ws1ElPluhi5okxE\nTwC4pe+S0eFHhRC/Hp7zowCWQohftgSRy/3ilSqXJdH02F5JhS8O7vPKtooAWAWqjDuuUVzBtxxt\nJeIuT0nsLZXVfi9A4Kv+oRl8XMKbefAXHkCEvhAYrNc48QKsR2us+2v5jkPrkEDwglCUVAsZ2vs3\nLVdCLNSq4AYU5SsBKeRCK6SyQUHpwmdWHIUeh9F4s5wjIAVTIBZN9RuW/CqFV/7zEAoy6X9VzIa1\nLEQkw0C2KGflLWnsq2cfPZR0Q7VABmVRPkyceUMT42g76gbKyCtKkJEst+nSJ8zL7HFWTX8GuaIs\nhHhHZqRE3wPgrwH4Zm33XQB3tN+3w31WfufTn46279y6hTu3brlOzST1kHSRzTg3S5zLsjVB3hPi\nTE/wwr+662qtuW/V+yUgclfLS5JFIRoEkhO3+Y6y8oHt3CJjF5rCHDeQR5Q+IeXNlt6i4aWrrHq4\n+g8Zpihly0MmOeW+ajwvvPIKXnwl3wCtO/r6XQB+CMA3CiHm2qGPAvglIvpJSLf1mwE85Qrn67/i\nKyqnIcsqjn5mnN+kIDPNoUQ0UVlrgqK2w5NBoZUMJC20SJAbqvizwnDFUaQgV01f7HpGMt8X6KMV\nYeMlL+6UCCf8ewzTHll5s45AmvWC7IRMBN54Hr9z8ybu3LwZ/f7EZz9rPa9un/JPAxgAeCJ8cB8X\nQnxACPEMEX0EwDMAlgA+IBpSuyIiDGRbxq5w6iaRreT6mAPD9L8A7G5aVag0K3tb1BHmzHAt+1Jh\n1rjveOxG98SWG8pMk1RpADftHcqi7ujrt2QcewzAYxXDLXJSeleBsFxhsyB3B+eob7gLh+meLvrc\nN1nhN+piawrNw1CETt4Dc7CYjfhN5M2Ehd1CfJ1Y0Sv3xvJcfwXCzIqjjT7kvGNMNbK6IqqO7K3K\nLrzfKu72gEWX6RBlLduiYpkb7pa6aTohymVuPOvMslZwW9Zx3jGmWRorhA3Fs2mqzkkXykrW/1rC\nqXLHnPuZXcHZr7wlOiHKRQp92YnbVcPIouhiGUzz5L07UzCdLvCabqe677dpYS963wqz4ol+F2i4\nlnYVNuC56GIjiNk8RQd8Nd5w1hqrVfqVqxgC3RDlOgNgGhDrpmBB3i5d7Ecug16A81rt+vSmOnGZ\nlYYrvK56CRimCnVnZbTZr9wJUS5C2Rtv2hqpc5xhFE3ky9oVQcbgLhZeZhexeYwqzeO3B77RvuVO\niHITFUFbCzQ0cY7rfK4AmaaoNWcT3ehLY5hWaWDK5CamRnVClMuyqWHvdY4XOZ8FmukC5opnXcmL\nXUkHsweEXT22lf4Sp2H7DdTOiXIXCuImBDnrnC48A6YaTXZlFK0gmujfypsGZQvXOWdcJqpSOhjG\nRZ3ZE0L7SxndN12gE6LcFRFq2l3Nfc27R5mRzG3FX3falm27lfgK9rW1ng6GQYkuHNM6LpL3Ntiv\n7G0klh2gS4LMldN2aHvKXJ0+30OjKw11Zv/oet46eFEmMr+r6T6vTJh1YWHeHEKI2qu6bUWQKf4a\nVp38klgoxLGdjLaAVVEjHV2vNJluY10zvwBF850ZatN19UGLctGHuS03JtM+m1hQZtsUbXhWeRa7\ncP/M4bAPi9N0ok9507Rt9dbtl+xiRtlHuirI2+hfbSrP1U03532ms2yoX/mgRLmNUdN5ZLkAuQLa\nHkWW7axzvExcTWEuN5hHIfEvMpALSH/7mmG6TIdHYB+MKG9DkLNgQd4OTa2o1VW2kjatbzt9qLvP\nijkczOl7ReYsb4u971Mu2p9mXsPsH4UHcmS8/33osyrCtlbIY5hWKOPZajEZRdg7S5mnIjF12GQe\nKONqZhjmMNh5UW6qQuOKcb/pmtt6Y/3MG4mFYfaMLfY5d1qU2xZKFuLDYFuCXDTeuqtqUdYnHLeY\nx7l8MdtiEwNp24qjE6LM4ss0TdnCkvUB9abj33Q/s/Ueiny+0UjnrvaPM/tD2QbsLs5y6YQoNwWL\nLwM0J8ibir/r7Nv9MEyX2SlRZtFl8mhSQNrIb5sUuLIfpwDS6auaWjM+LrsMU4xOiDIXWKYuVcWu\nabf13sDWMcNshb2fp8wwLtoQ5K64eqvcQVfSzjBtsCvNbBZlZudp8kMKbQpya6JXcEWtTVn/B+9l\nYLrJjuRLFmVmZ6n6ycVd7yvOo8zdWT+XWGS96x2p4BimDF0oxyzKTGfYxEcemlxCs25aGo0ntJZJ\n+1033Cp3xWLN7AJdzqWNiDIR/RMiCojoXm3fo0T0LBF9joje2UQ8zH6h1iVX/4uIW1XrWMWXlY6y\nlEnLRtaSJor/a8dL3VuNdLIgM10nyqMtfgGwLrVHXxPRbQDvAPB5bd/DAN4H4GEAtwE8SURvEV3w\nDTBbZVuDqDa5OEhZ9kHM9uEemMOgsrdsQ/LVhKX8kwB+yNj3HgAfFkKshBDPA3gWwCMNxNUopqXG\nFUs7FHm+bQ6S6rIgdwrtvsw7tN1z0XLD5Yppko2Vvy2V81qiTETvBvCCEOIzxqEHALyg/b4b7usU\nyv2o/2eapUiFvEuCvC3q5k1ugDDMbpDrviaiJwDc0ndBNqQ/COCfQrqua/Gxp5+Otm/fvIk7N2/W\nDZLZMk2JXte+69t0eja97nvil77+NYsvs6d0Zf3rF155BS++8kruebmiLISwii4RvR3AnwHwByRL\n/m0Av0tEj0Baxg9qp98O91n5ure/PTehzO7QFUG2paNLgtwmdQZ3NXmXu+aRYBgrDZT9O4bB+YnP\nftZ6XmX3tRDiaSHEG4QQbxJCvBHAiwC+SgjxCoCPAvhOIhoQ0RsBvBnAU1XjYnaDKv3yzuk4NV3W\nuyDIVdejzgtjU6OtGWbXqTTrooV06DS59rVA6B0TQjxDRB8B8AyAJYAP8MhrxqQtQW6apr7HrIfT\nZANBD7vuCl56yFxkmUPBWVa2UAYaE2UhxJuM348BeKyp8Jlu09TCH10b1LXtqU9F429jUZS6sOua\nYcrTia9EMbtLEy7YIserpKErYrwr8bus5G0/C4YpSpm8uolGY5U4eJlNpnVYkDdP6b79LcfPMLvA\nJmoGtpSZVmh7+cm25h/vgyDrsDgyTLNjOtqGRZlpjDKC1jXreJ+oayVva6EShqlDGx+p2QYsykwt\nyhaEtkZWN1Gw6lj3XSvYRdkvvwDDSKqU0bzjrrLStHeN+5SZjVB3GdMuC3KXqD0liq1kZgcw8+mm\ny2Wb8bGlzLRKU3N8y+wvyz4KMsMcCptY12CT3xZnUWZaoSkBa1OQm+oD74IYNrEyEVvJzC7R9UZy\nVViUGQDZFWqRzL+JDzVsY03tvHM3LUR1+7PL9os1sXgJw2ySKnk2sW2GlxF+Gw0DFuUDZtufVbRR\n1TJWR5tKTVdb4a4lNW3nMQxjp6rbehP1AovyAdLUKlxNU0WQy0pPk1b/vggfW8nMIdO1b42zKB8I\nXf1Qg6IJQc5LzTbc8G1Qd2oHwzD5bMNKBjoiyus/fX3bSdhbqtTPeXlvcy5rx/kO+1hkyHKRJHdt\n3dwqlE2W65b34Vkwh0dZ746ZdVXdYqtLzKDbGm/SDVH+Ehblpint1s061nIL0cy8zrRX+LxakZSX\nvb8ui1AT750FmdlFSgty+oAeWLxZMa6dFmWmPlVefxPu3qpYR1dnX5DeV0OM92VJPkVTjTAWZGbX\nqJJnuyrIAIvyzlL1lXdNiIGSlnHeaPCcNOybGANsHTOHSyODEQsK8kH1KTPFaEOIq87ZLVMxl7KK\nK7iogXbc1DI53RYgFmTmUKk0Hzl5wAww3qwYlxlfFViUO0SdV9mUhVh38EIjbuk4Me5DmanYT4vY\npGhK6zbKdumZMPtNVXEs6q4GkuVlGw1WFuUtUff1NWUZNlEpN+KWjhPkPpSZisNZJrJMKlmQmUOk\nkCA3JMZmfHVhUd4QTbyyJl20TVXItUdOxwnKPpyTjk0sSt8FWJCZQ6VRd3WGIJel6XLCotwwTb6e\nTYsw0LAQy5OzA9vwdCadXRMdFmTmUGm8/1gPu2JcZnxNwaJcgbaqq6JZoemVqTbqnlZsyU29q2LD\ngswcKo3OP5YBxpsV49LjaxoWZQubqJKaFOAy5+k0OmirZSGWlx9Gn7FJEwO6gPZWIGKYNmhEIFsS\n5DY5OFHeVrXTBREGtiDEwFb7i4HdFZumrGMg+xnu6vNh9pfao6xz1jjoyqAuGzstyl2sSsq83rKZ\nYZNCDLQrxmwZZ9OkIGfGs8PPiNlPDlmQgY6IcterhbKVXhUxaVuASguvvKh4ArY4lUmxLwLTlLva\nGvaePCNm/6g9mEseNAONNyvGZ8bZNp0Q5W1SVw62IcBAtVHSiWPui8olZMvWMLBfQlPkTliMmX2j\nzdW5gO5bxzpe3QCI6AeJ6HNE9Bki+pC2/1EiejY89s668ZRBlPhfKDwhnP+rXFsVIor+lz2HjP/G\nRfH/IggR/zcPIf/ZNvEsAOQ+i12jLUFmmC5TaZBq4sf+CDJQ01Imom8C8O0A/pwQYkVE94f7Hwbw\nPgAPA7gN4EkieotwPJEuVDR1BaKNkXt1rOHcKxtyTQObs4gV+yTEirw76kIZYZimqf2Fp5bmH+vx\nbZq67ut/AOBDQogVAAghXg33vwfAh8P9zxPRswAeAfCJmvGVpkkx2PR3hcue11j/sILFeCNk3RWL\nMbOvtCLIYZi7KshAfff1WwF8IxF9nIj+FxH9hXD/AwBe0M67G+6zkuUervu/KE2GVQTdzZzlhs06\nr7BLuqxV7HBNA5t1Tyv2zU1dhDLdKwyzS5StG1KC7KrTLIJcpR7adl2TaykT0RMAbum7IO/7g+H1\n9wghvpaIvgbArwJ4U9lEfPKJ3462v+RNd/DAQw+WDSKTbU8Kr/uZw+hY/sWF40nQMYtYse3CsSlY\nfJlDoZXpTjLgWnHZ4myau899AXefeyH3vFxRFkK8w3WMiP4+gP8cnvdJIloT0X2QlrGurLfDfVa+\n5h3fkBV/XhI7RVMCDLQowkCuEAMsxgzDNEdtQXYHHG9WiMsWZxs88NCDCYPzU0/+jvW8uu7rXwPw\nzQBARG8FMBBCfBHARwF8JxENiOiNAN4M4ClXIJtyGTeFzaVcxMWada7piraGVNUlrSjhmt6ke1px\niG5qhmGSlJr2pLmsuyzIZag70OsXAPw8EX0GwBzAdwOAEOIZIvoIgGcALAF8wDXyels08RIKD8wq\nFlittCRowB0tg+nGwDaGYXaTKv25RQdzAfWmO6n4ukYtURZCLAH8LcexxwA8Vif8MrT1cMuGu3EB\nBhoTYRkUC3EZ1AALhmGStDndCdi9+cdF6cSKXtt8QFXiLnVFG/fWUH+wDKp9SelyAajCft0NwzRP\n7SUscwZ01bWQu0wnRLkpmqr8K4fSlvg0KMJxkGwRl8V2R/tVHTBMPepYrCR3ZAUu/9SIz4yzi3RC\nlDf5kGrHtIm07qAIA93P7FVx3RULMsPE1Bbk7MDlnxrxmXF2lU6Icllaf6ybenElM1WX+oYVu5DJ\n68CCzDAxrjEUVQW5yqCuqnXbrtRVnRDljT2qbb+UFkVYBs9CnIX+fOrMEWdBZg6RpgY1siBn0wlR\nLsQOPdSy4htdViqKzUrDLmVqk1KDTvLCqpcUhtlJsgS51mcXC6zQVbfM7Vrd1Q1R7vpDq/sFqdLR\nbb/q37WMrJP1/Kp+VatshbTLz49hFEqMm/jAQ6HpTgfUd+yiG6K8bRoSwTqhbFuIdzkTK/Keoe0e\ni9x1E31oDLNrNN1/rO3IPP+QBRnYd1FuoeKsG2KXKvNdz7w6ZQW56J1XsRD26bkyh4fKvY0P6NrA\ndKd9oBui3KGH33RKupix9k002hBkto6ZQ8QlyK68r/anyljFD0rkxZfHPtRt3RDlDbCJ6rSrlfY+\nZFQbZa3WNq1jhtl1sixkk6LjNqoM6DpkQQZ2SJS7Ui3uSgW9LxnURpXRnpaDKrBk2BXiYZhdJk+M\nzXJg/k6IcJk1rFmQrXRClLtW9e1iZbxPmTKLKn26RRe5b0KQD+U9MPtFEUG2lYlMQeYlMyvRCVHe\nJLsouCb7mBHzqGodZwpygy11W9yp43rUlWNhmOYoOsKaBXlz7Jwo74OolmVfM59O3cUJSlvHDS5O\nYIvfmg6G6RDOJWQrzjCwCTIRJcNjQc6lE6J8iEIL7HfGKoLel1VncFWdZfvMuMosPNL0QiQM0zZZ\n3pomrWMlxlEYhhhz95CbTojyIXAImakoZSqGzHDKjPK0xFlUjM24bL9T6ciJm2E2TZn5x1UFWe3P\nso4P1QgrCotyw7D4ZlN2HqQznDLfYZURWK3jplxnVec+M0zb5I1lKNI47YIgH0rdyqJcgUPJHE3S\n1CpBVeYdw1HpVFpIv6IgM8y22VVBPrT6lkU5g0PLDG3RtCDb+o7LVApVGwGpY4VDYSuZ2Q5Fy16e\nIGd1FVUVZNdqYLb4D4mDFeVDfNmbpikxBvIHlaTCtFQKdRfSr2SlgwWZ2Q5NCbIZXmbZKzDNsJHF\nf/aYvRHlQ32BXSTrTVTtw82qECyRyD8Z8Qkh5HlEdpe0I56yuYwFmdk0dUdYA+n8X6j85Ux3yhs8\nVmT/IdAJUT7kF7BPmG+xqcFcLndZKuwcd7V+TpF4XXEVhQWZ2TRtCLL2I7UvFWaB8K1hF9h/KHRC\nlJndJU+IgbQ4mgsKlBlIYgs3/GF1VdvcalpkVhd11c88JtJW4RqGqUpRMRZCpBf0UGEYZS2rDBbt\nP7bFx2KcDYsyU4kyVnFWC722IFtc1UUtY2vcNVzVUdoqXscwZalSDtsSZFc5Z0EuB4syU4qyLuoi\nK2SVEeNEmFnWcZ7rzPNS6Ugcz7zaDQsysw2qDuYCKgqyY3W8rPhYkIvh5Z/ihoi+hoieIqLfC//+\nRe3Yo0T0LBF9jojeWT+pzDYhZAuyObo5a7Sz7iJWg0hI/ohcyonj2ghPU5DVPv1/FUFW8djusygs\nyMwmqTu6GqgnyEKLW28QNzXl8FCpayn/OIAPCiEeJ6JvA/AvAPxVIvoyAO8D8DCA2wCeJKK3CF5f\nbeewFZmswRxFB3ckKgGb69ioDFzu6lKDuIz+YpvL2nVdVticqZlN05ogZ5TBcCMRb9m48vYxNS1l\nAC8BuB5u3wBwN9x+N4APCyFWQojnATwL4JGacTEbxmYZZ013yCuUif+IK4E8y9isEALTOg6CUoJs\nxqfulYzzo/8syExHyPLk2LxVznBcnioky4cRgfxTIHxbXHn7GEldS/lHAPw2Ef0E5Hv9+nD/AwA+\npp13N9zH7ABlrOMiBdOsAMIfqWOJ+CyjpytZx1pcrvhSYpxMiDNYFmRmk5QZYZ0Zjq3MVZxumBdf\n1gI8jJ1cUSaiJwDc0ndBvpsPAvhBAD8ohPg1InovgJ8H8I42Esq0i7P1rW+XdFUDcSH0VGHMcVWr\ncG2DuYIgiOMt2hNis4rN/jPtXCMhzmBZkJlN0oYg5w2sdM1BZkFul1xRFkI4RZaI/oM6LoT4T0T0\ns+GhuwDuaKfeRuzaTvHU478dbT/w0B088NCDecliGiKrmLjcVJWs44JiHP219B2XEuMwTpurGqgu\nyCzGzKZpTZALlMUm+5APnbvPfQF3n3sh97y67utnieivCCF+k4i+BbLvGAA+CuCXiOgnId3Wbwbw\nlCuQR975DTWTwRSliEUMlHdRA+6BI3kDqhoXY8QjrK2WQZguR2KcYbIgM5ukiBgXFWK1Xdg6dohx\nmfiy9h0iDzz0YMLg/NSTv2M9r64ofz+Af01EAwAzAH8PAIQQzxDRRwA8A2AJ4AM88nq7bFyMkRZG\na/wOV7UQInJXmxVEFK8tLRYxTp1bspLgjMtsmqKjqwuFVWC6UyJcFuStUkuUhRCfAvCXHMceA/BY\nnfCZeuQVhzqDt6I4dEs43tmodVw0NbmCnFdBWO6bBZnZNE1MdwLyy6Y1XBbkrcMreu0pbfUVR+Hn\n9BmXtY7LirEZOnlePeuYBZnpAFnTnbJ+p8IxBTlvdLXckH9QTvxtcTPVYVHeI6pYxrYCp4RRX/kq\niqOmGOvhV7WM9RgacVc7RpkyzKYwc2gdL9YmBZlHWDcPi/IeUMQqjn47ClxiMJVaNEOPowExBgAR\nBOkpTupvxn0AaaG1LgZinJMJW8dMB2hjdHUUbkFBLju62hYmC3IzsCjvKE30F1tHNTumLlUR40R6\nbNZxQTEGSlrHRdLE4w6ZDtDadCf5I3dMh02Qy0x5dP1mqsOivGM0NXhLCWIy8LTg6mKcVRCJ0t9o\n1d3SdQdxqfTlLgRSUZBZoplN08SUJ8AtyC5Ml3XV+Jh2YFFuCSFErcxb5Mqs0ZnmXyKyCzGQ+mqS\nTfhcgqzH47LCzbWqSUu7WQm4FvjIHMhV5jnraQzDYEFmNklT1jFgNJzljmi/+msVXK2hXCaurN9M\nM7AoN4BLWMoKc9EzXdaw/jtVAC2FLtMNTJQpxmYcqXiFSImxInBUAFZB1tKR129ss9a1BCZ/Gn8Z\nZhO0Lcgu71XWgC4W5G7BolySsgtq5J5XJm5HOpyF2ZVWU+QsYqzvN7ez4jb/uwaRpJNkF1zPbDho\nx8zrywoyw2yS1tzVji4dmyBX6T92xW89Bi5fdWFRdtB0Zo3OKZsOfbuOEANWMc4q1JkjNpEs2FXF\nOBGfJZ1FxNgJizHTEVyC3PR0J6sgG+M3ygqyq3GeOg9cxpqARRnVBBgo58KpahEDBcS4gBDLzXSB\ndgmw+buwGIfpEShRyRRMa9H0qTQwTBdoU5BtgmkT3SqCXMZlzYLcHAchylVF16QNl7ROYas4734y\nBmrluajVL5uVm2cZq+tKF3gjnVZXtu06SxrDHcmfualhmObJclfr1BnQpY6ZYbkayFW73/IEmWmO\nnRXlpoQ2i7YsYUUpi7jI/bpGThvHrMe1MFwjuFNibMwztqU/bxpVlIaM9NjCMeMLd6TPSe1hmPbJ\nE2Rz8GNmWA5BNgdiulzWmxBkLmfN0QlR3oTA5tG2ACvaEGL5J98q1s/zLJYqDDd0UctYPyc7qW5r\n3SXG+nU6RfLM9nMVc4gUGdBVOCy9bGSM+cjqQy4Tb1FBLuoFYMrTCVHeFFWH8dd1zxQWYsc8Ynui\n8mCIvB4AAAeKSURBVPtfbYXXKcaWdBV1UxcVY5W+zG8dFxhUwu5qZhco2wA315qvK8htWcgsyO3S\nCVHe9Jy3NmLLaxFXGi1tYva5OoQ4b0pTFKcQ0ZzhQtaxlv6sgu5yRetinNWIsMGuaqbruFy5LnHU\nRdRcwAcIG9A5YmwO9mxLjM09XNbaoxOi3CSbkHdXH1Hh3xWEWG66xczVV2skIkq/TWh1i9kU4yIF\n3TVozCnGxv25YEFmuozLcswaXZ0lyDbrOHO2QQ1BTtxHQeOIy1q7dEKUuz56b6MiHJLn4rUN6AKS\nYm3GnSXGgSHAZcTYmY7wPvLc1FnkPUeuIJht0qYguxrXbbireTBXd+iEKHcFV6azZfKmRdgltFlC\nnDg3PMdMg62w6i4vfSlMZSlbz7Um2W4Z29zsZfqMrfFyvzHTMYqIlbX8CJFd3h1lxFZ+9fJaxjpm\nQe4uBynKpQuRZX/h0dEmlv4bq7jl9AunxNhmARvbLhe1TYzNbTPN1jQY9+NqMBQewCV3pndZr2aY\nzZElVnkNWlv/cZ7Lui1BzjxPhV84ZKYJ9lKUi2aiogKc2F+yPzj5097vS5ZrnP3Cxj4hE+a03G0u\nrqxpTa54XPus/d2WfuOsMFmMmV0jz6tWRCSzBNkML9rOKbNF4nP9ThwDl7dt0UlRbjIzFMm0hTO2\n6zxH5i5k5RrXZ4lxnIzsxoRLiG3HipAlxNa/yZOt4TifueveCqaVYbZFWYHMm4Fgs46BtCCXaQC4\nfptwedsenRDlqhmgygjDKtdEFBBf277EtiO8ImJsUqTQ2sRYT4u+J7dvKUeMo/O0c23ptRxwxsuV\nA7MLNC3IiTAtjei2LGRm+3RDlOsI5Qbj8IzJ/SapzJ8+wXmuQgj5DWb9r+1YuEP+Ma7X/0bbrv6t\njPsB7P1f+nah4m3GHfaBp04Du82Y3SJv3IVrv6vx6prqlBVfnsiyCO8WnRDlTdBkxiwyWlE72Xpe\nViPB6pbS9plXNtaosQw6Sx7OaXQkE1X4mHBsM8w+kdeQtTWe65YHFuTd42BEuQpFMnTuGXohc4VX\npN/belmxIktE7rht52b8jvYXCs0OCy9zSGTWIy113bAY7y4HJ8pNZNbKIRQU0SZES91nFXdX0XOK\nwALMHCJ1GrNVywwL8X7QCVE+1My0CcEqPCex5DtgsWWYzTdeD7WuPCSyRy4xVl587gtyqlHN/3WJ\nluUrIbyu/13j7nNf2HYSDgJ+zu2jP+O26oJDZ5/yMYtyBe4+98K2k5AiS3C7KrxZdPEZ7yP8nNuH\nn3H77NMzZlFmGIZhmI7AoswwDMMwHYE2sXBHZgKIuEuFYRiGOTiEEOmlLbYtygzDMAzDSNh9zTAM\nwzAdgUWZYRiGYToCi3IGRPReInqaiNZE9NXGsUeJ6Fki+hwRvVPb/9VE9Gki+mMi+qnNp3q3IaIf\nI6IXieh3w//v0o5ZnzlTHiJ6FxH9YZhPf3jb6dkXiOh5IvoDIvo9Inoq3HcPET1ORH9ERP+DiK5v\nO527BhH9HBG9TESf1vY5n+su1xUsytl8BsBfB/Cb+k4iehjA+wA8DODbAPwMxROB/w2AvyuEeCuA\ntxLRt24wvfvCvxRCfHX4/zeA3GfOlICIPAD/CsC3AvhyAN9FRH92u6naGwIA3ySE+CohxCPhvh8B\n8KQQ4m0A/ieAR7eWut3lFyDzq471uRLRl2GH6woW5QyEEH8khHgW6SVr3wPgw0KIlRDieQDPAniE\niN4A4EQI8cnwvH8P4Ds2luD9wVaArM98o6naHx4B8KwQ4vNCiCWAD0M+X6Y+hHS9+h4Avxhu/yK4\nTiiNEOK3ALxm7HY913djh+sKFuVqPABAX0LmbrjvAQAvavtfDPcx5fgBIvp9IvpZzSXleuZMecxn\nyfm0OQSAJ4jok0T0feG+W0KIlwFACPF/AdzcWur2i5uO57rTdUUnPkixTYjoCQC39F2QBetHhRC/\nvp1U7TdZzxzAzwD4Z0IIQUT/HMBPAPi+dCgM00m+QQjxEhH9KQCPE9EfIb28Nc9DbYe9eK4HL8pC\niHdUuOwugDva79vhPtd+RqPEM/+3AFTDiJ9tc9wF8KD2m59lQwghXgr//j8i+jVIt+nLRHRLCPFy\n2MX1ylYTuT+4nutO1xXsvi6O3s/5UQDvJ6IBEb0RwJsBPBW6UE6J6JFwYMF3A/ivW0jrzhIWLsXf\nAPB0uG195ptO357wSQBvJqIvJaIBgPdDPl+mBkQ0IaLjcPsIwDshB4t+FMD3hKf9bXCdUBVCuh7+\nnnBbf647XVccvKWcBRF9B4CfBnA/gP9GRL8vhPg2IcQzRPQRAM8AWAL4gIiXRvuHAP4dgBGA/65G\nDzOF+XEi+vOQo1ifB/D9AJDzzJkSCCHWRPQDAB6HbJj/nBDic1tO1j5wC8B/CZcO7gH4JSHE40T0\nKQAfIaLvBfB5yJHBTAmI6D8C+CYA9xHRFwD8GIAPAfhV87nuel3By2wyDMMwTEdg9zXDMAzDdAQW\nZYZhGIbpCCzKDMMwDNMRWJQZhmEYpiOwKDMMwzBMR2BRZhiGYZiOwKLMMAzDMB2BRZlhGIZhOsL/\nB3jfi/P+7WJhAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "icount = 1699\n", "fig, ax = plt.subplots(1,1, figsize = (8, 5))\n", "extent = [mesh.vectorCCx.min(), mesh.vectorCCx.max(), mesh.vectorCCy.min()-mesh.vectorCCy.max(), mesh.vectorCCy.max()-mesh.vectorCCy.max()]\n", "ax.imshow(np.flipud(U[0][:,icount].reshape((mesh.nCx, mesh.nCy), order = 'F').T), cmap = 'binary', extent=extent) \n", "ax.imshow(np.flipud(v.reshape((mesh.nCx, mesh.nCy), order = 'F').T), cmap = 'jet', extent=extent, alpha = 0.2) " ] }, { "cell_type": "code", "execution_count": 21, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "(99, 1700)\n" ] } ], "source": [ "print data[0].shape" ] }, { "cell_type": "code", "execution_count": 30, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "
\n", " \n", "
\n", " \n", "
\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
\n", " Once \n", " Loop \n", " Reflect \n", "
\n", "
\n", "\n", "\n", "\n" ], "text/plain": [ "" ] }, "execution_count": 30, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from JSAnimation import IPython_display\n", "from matplotlib import animation\n", "fig, ax = plt.subplots(1,1, figsize = (8, 8))\n", "ax.set_xlabel('Easting (m)', fontsize = 16)\n", "ax.set_ylabel('Depth (m)', fontsize = 16)\n", "extent = [mesh.vectorCCx.min(), mesh.vectorCCx.max(), mesh.vectorCCy.min()-mesh.vectorCCy.max(), mesh.vectorCCy.max()-mesh.vectorCCy.max()]\n", "ax.set_xlim(extent[:2])\n", "ax.set_ylim(extent[2:])\n", "\n", "nskip = 79\n", "def animate(i_id):\n", " icount = i_id*nskip\n", " frame = ax.imshow(np.flipud(U[0][:,icount].reshape((mesh.nCx, mesh.nCy), order = 'F').T), cmap = 'binary', extent=extent) \n", " ax.imshow(np.flipud(v.reshape((mesh.nCx, mesh.nCy), order = 'F').T), cmap = 'jet', extent=extent, alpha = 0.2) \n", " return frame\n", "animation.FuncAnimation(fig, animate, frames=20, interval=40, blit=False)" ] }, { "cell_type": "code", "execution_count": 23, "metadata": { "collapsed": false }, "outputs": [], "source": [ "def traveltime(x, h, v):\n", " d = np.sqrt((0.5*x)**2+h**2)\n", " t = 2*(d/v)\n", " return t" ] }, { "cell_type": "code", "execution_count": 24, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0.0500975610684\n" ] } ], "source": [ "print traveltime(2.5, 20., 800.)" ] }, { "cell_type": "code", "execution_count": 25, "metadata": { "collapsed": false }, "outputs": [], "source": [ "def intercTime(x, h, v1, v2):\n", " out = 2*h*np.sqrt(v2**2-v1**2)/(v1*v2)" ] }, { "cell_type": "code", "execution_count": 26, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# print intercTime(0., 20., 800.)" ] }, { "cell_type": "code", "execution_count": 27, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0.05\n" ] } ], "source": [ "print (2*20)/800." ] }, { "cell_type": "code", "execution_count": 28, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 28, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUEAAAE4CAYAAADFI0E4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsvVuMZFl2nvetfTnnxC0z617dPT3dzW5SY4IkBF3IMUlZ\nBPREgbBeDBiwX6wXA4ZkGzCgB/uFLwb4IL9YNAyCgCSABgRQMGyLAiTAgG3JIgkPRYojkZwbZzjd\nPT0zVd11yazMuJzL3ssP+0TEyaysruru6qisqv0BWRkZceLEiQDir7X3WutfoqpkMpnMy4p51heQ\nyWQyz5IsgplM5qUmi2Amk3mpySKYyWRearIIZjKZl5osgplM5qXmsSIoIl8Qkf9bRP5ERP5IRP6r\nRxz390TkT0XkqyLy55/+pWYymczTxz3BMR3w36jqV0VkCvyBiPyfqvqN9QEi8ovA26r6oyLyM8Cv\nAV/+fC45k8lknh6PjQRV9ZaqfrW/fQJ8HXjtzGF/A/iN/pivAPsicuMpX2smk8k8dT7RnqCIvAn8\neeArZx56Dfje4O/v87BQZjKZzIXjiUWwXwr/r8B/3UeEmUwm89zzJHuCiIgjCeD/oqr/5JxDvg+8\nPvj7C/19Z8+TG5UzmczngqrKp3neE4kg8A+Ar6nq//iIx38L+FvAb4rIl4FDVb193oHfuHX0ya/y\nBedX/+6v8F/+nf/2WV/GhSN/Lg+TP5Pz+dLN/U/93MeKoIj8HPCfAn8kIn8IKPDfAW8Aqqq/rqr/\nTET+uoh8G5gDf/NTX1Emk8nskMeKoKr+DmCf4Li//VSuKJPJZHZI7hi5APz0z/78s76EC0n+XB4m\nfyZPH9mlqaqI6J/ePkIB1bSuzmQymc/Kl27uf+rEyM4jwcoJ3gjmU11uJpPJPF12LoIjZ/BWMJJV\nMJPJPHt2LoKlzZFgJpO5ODxpneBTw2iHUcgamMlkLgLPIDssZAnMZDIXhZ2LYBCLikFEMGQ5zGQy\nz5adL4ePm0hQJSpYA1HTTy6XyWQyz4KdR4LvH9XcX3Z0UfFWsEbIieJMJvOs2LkI/l/fvsvXPzrh\nuG7xRrCSl8SZTObZsXMR/Of/3/v8yXfucXRcU5gcCWYymWfLzkXw3a/+O+689z3qkzkikqPATCbz\nTNm5CD744Jss7v2QZnlCRFHycjiTyTw7di6CxfQyUk7oxFN3Soi6uZAshplMZtfsvETmrS//Bxy8\ndZ043uO4iVgBY8BaiBFCrpXJZDI7ZOdWWv/4qx/g+/7h0hkOKsdB5SidoQ1KF3PNYCaT+WQ8V1Za\nv7R3h9fjPe7eP+Tf/vCY7x+tWHURK8lZJmeKM5nMLtn5cvgbv/zL/OGVH+f3r/w4H914B/Mjl3hl\nWnBt7LOzTCaT2Tk7F8Hf+t+/yXs/NuG9H7tKG65xcnlEaEIqlxFNkWBeD2cymR2xcxH8/rLj7iow\nXwVMHWnb2O8Fpn5ieputrIOZTGYX7FwE//Llim9fGWGvjZlfGeHGnhZYdX3NoAgGzaYKmUxmJ+xc\nBP/af/5l3rz+47x5/W3uXX2FL+xXqLMc1R2lM5RWMAbI5TKZTGYH7LxEZvXD77CQirmUnETHvWXH\n/WWLAlfHnmsTjzdCF5U27uzSMpnMc8xnKZHZeST4Lw9H3OSYm/pDboYVdxYl788rVqaCq2MuVY7S\nGkRA+ra6TCaT+bzYuQj+g9/+M/7i4jv8xfl3uNHe5zujt/mT8du0l19lzxveOqgwYnMLXSaT2Qk7\nF8H33j3k2p1b3LjzLmF5l1uvTLhz4yrqLrNYdrRBWa/QRWCHq/VMJvMSsvOOEVUFEcQIxhmMEcQK\nGCECbVDaqERVhGyskMlkPl92HgnuX59w9dI1bn7xTV6XKxzOXmc+u0Z9acb+tKRDWHWKNWAFVPoZ\nJDkizGQynwM7F8G/8PZVvjSd8cXZj/GFkWJkj+vscSIV1gjLEGHVMi0sU2/7wumsgJlM5vNh5yL4\nU1+4xLVJQTfxHI08B1a4boR5G3j/qOb9oxXWCK9MS6aFwYgQI7lyOpPJfC7sXARf36s4aQJfv7NE\nu0Pe0Pu8offRruOwnvLuakoxmjJ2hpvTIu0ZAlkFM5nM58Hu9wRLxw+PVnz99jH3795jsXqXYvld\nXGz5YfEm7xVvMLlieGXiCVHxNmWJc6Y4k8l8HuxcBEsndCHyYNlx58GK+w8e8ODoLkVsOdm7wsms\nQWaBuot0qgRV1l0t2Vghk8k8bXYuglaSo/SktMzGBROZMCkOKOg4GE25Mq6oRh7XO007s51Bor0K\nZiHMZDJPi52LoIiwN/K8emnEzF/iNaNctXuMjPK2uYbYq8RywsG4YBUUkbix48/OMplM5mmzcxFU\nlFnleG2/op54rk6vUswcI294MxiudcJJGzlpAvcWLSNv2K8cVek2z897g5lM5mmxcxEEGHvL1XFB\nHSKdGD5YCqNGObA1N2xHReRBrfzgJLJXeQorXB45NKbC6Uwmk3la7D4SVCiswZbCsot8OK+5fdLg\n6hPe1juU8Q51sNzrLvFBe8DV/SkHhUPYZonzmjiTyTwtnkkkaAWsFdqgHC9avnd3Tnxwj2nzATea\n95hTcOi+yEfW43zBsq1yBJjJZD4Xdi6Ckd4UQcCIoECISgiRrgvELqASwUSMS2XSUZOpggBonykm\nB4SZTOaz80wiQSXV/FkjTErH5VkJZsqBXmOiESue6+Yab9gZs2lF4W2aTWwEO5xNnLPFmUzmM/JM\nRHCNFWGvcry6VyETy+VqRjV6g0Isb3Ylo1ASxeCMcLgKjJxh5A2VS+UyAbIKZjKZz8QzE0FNo+WY\nFJZrFECBKR1HlaW0hssCVwXurzo+nLfcOqm5VDms8Yyd66PJ7C+TyWQ+G88kO7zGoIy9pbSWOkSO\nm8DtectEGl5xNa+4mvlKWC4ctxYWQ8WstFhDLpfJZDJPhWe6HAbBChib5g7fmTd872jFXnOI4w5X\n9CMWccL9eIUfhiuMrOHa2A+fnpfDmUzmM/HMEiP9anhjoQ9Qt4GjRUNczFmEe7TdLVqzT20LVmaP\nuou0IRLWUaBmHcxkMp+NZxwJ9iIoYA0UzlAVlioUuFghYYw3I0amYmoLCmdQEZq4dZWRbKqQyWQ+\nA89cBNc4Y5iWjuvTkunkgH3vqYqr7MWSm92EVZiwX3mcMSzaramClVR7mM0VMpnMp+GZiqBqEjBI\n5TKz0hGBka0YTa6iE8ckwhdbuNzCsgssu8idRcOssMxKi3P9YjorYCaT+RRciEhQNbXSTb3FGwGE\nRRDePYGJF6ZeuDYWvn8ceLAKfLRogYLKWcZOeiHNKpjJZD45F0IEIe0LFlZw1nLSBG7PG+4tWq6Z\nFW/4BZf8krAqOF5VfNiUVM5wqXKI2M3eYLbYymQyn5RnLoLbTLGkBAkQo3K4aHn/cEkX7nHAHVq5\nw1IvcxyvcmQusxh52tB7C2bxy2Qyn5JnLoJDNtlekmlCFzSZKmiLUoO0IDGVxSgEVTpNQrgutcma\nmMlkPgkXRgRVtzNEnE3GClemBQdxwthcwokwZp/LOmUhJZPSgkAdFNtHkXlJnMlkPikXRgSHOJMy\nxdeCcslfoqr2oHqNWed5pfWUncMZg0bhpA5UzlA5k9rpyEKYyWSenAslgmvx8kbYK1KmeOQqYmm5\nX1pE4TpwJSr3lx33Vx3HTeSgcngLTgySBTCTyXwCLpYI9j/GCKPCUnpLEyJ3lx0/OGm54iPXisgl\nGzlpI4uVctSmrPKssJjeVya30mUymSflQongGtPv8QEc14EP5w23ThresieM/Zx9u6Spxxw3E+4z\nZlY6wjqM7PuR85I4k8k8CRdSBNdZYgPECPM6cPek4TLHrMxdAg9ouUKLpXUVbYh0UYl9pjjrXyaT\neVIupAgOMSZliwtv8DiMFIiUeApKcVTOYo0hKnR9D56wdZfJgpjJZD6Ox4qgiPx94JeA26r6U+c8\n/leBfwL8WX/X/6aq//1nuajYD1OKgDeGWeG4Ni7YL/YZlRXOX2HajrjajXBaMPaWqKlcxpmUXc5L\n4kwm8yQ8SST4D4FfBX7jY475f1X1P3w6l5RYC1hhhf3KYY2wPxrjxpZQWWYtvN4ql1ul7iJ1iDQx\nMvGWaWE3XoVZCDOZzMfxWBFU1d8WkTcec5g85vFPTWENs1IYeYsV4agV5kEZOcN+adkrlVsnDfdX\ngS4qVmDsDdZIXgtnMpnH8rT2BP99Efkq8H3g76jq1z7rCYc9xd4KzsCDuuVo1bFsOl4rO75QtoxM\nJC4s84WlEcfE2342saCiuXg6k8l8LE9DBP8A+KKqLkTkF4H/A/ixp3BeYJslFoG6i9xdtNydLxn5\nE666YyrT0nQzFt2M2o1pQqBTxedJdJlM5gn4zCKoqieD2/9cRP5nEbmsqvfOO/5X/+6vbG7/9M/+\nPD/zc3/lMS/QC6GAqtKGyKqNtLEmdEvU1KAlRgfGClGJfaZ4aKqQRTGTeTH4yu/8K37vd3/7qZxL\n9AnWiiLyJvBPVfUnz3nshqre7m//NPCPVfXNR5xHv3Hr6JNdIFsRvH3S8MPjmnuLmnf8nHf8nAPT\n8m6Y8W6YsrQjrow9l0c+9RL3RdfrJXEWwUzmxeRLN/dR1U+Vm3iSEpl/BPwCcEVE3gd+GSgAVdVf\nB/4jEfkvgBZYAv/xp7mQR6GA9AJWWNPPGRFmVYGr9hAHs85xs3Uso8EZoQkpDCytUDqDkMpt8t5g\nJpM5yxNFgk/txT5FJAjb4uc6xDR2MyrTwjIrLKUzBFVChFUXeVB3PKi73onGslcmnY+ah7VnMi8q\nn2skeBFY7+c5K1hjiao0Qbk1bzEol0vhcimMXORkFTipO4wxFNag5XaucSaTyZzluRDBNUYEkdRP\nfL9p+Wje0rQNWjbMqgaD0q4K5rXHupImWlQV+g6SXDydyWTO8lyJ4HpQO0ZognLcdCyXNYv2mFCf\nIAKhm9GFPdQoIaRlcrbWymQyj+K5EsEhKWMsWCNIP5hEBLwVSgTp7w+qiG4jQciCmMlktjxXIrge\n1q4KpTXMSkdlSsal4IoKa2DcVVzqSqJxOCu0AUBxRjDSx4S5XCaTyfQ8VyII9Oqlm7nDsbRMxjNk\n7DBWOGgU2yrLLm6WzKUVxt7ifbLcygKYyWTWPFciOOz6qLylcDYlOgTu1uBNqiW8MoJ5K3w0bzmu\nOzpv8FYQMTlBkslkTvFcieAaEcH1G3xdiBw3gQd1wBF4pVSmVaToIDTKcQ2KY1LYZ3vRmUzmQvJc\niiBsM8VBhHkT+PCkxoaaWVWjzQqipatLVm2Ft5YuaI7+MpnMQzyXIqi9qUK6rbQxsuoitu1oZUXU\nOeCgMwgVQG+/r5tM8tBYIZPJvLw8lyIIW0MERfDGMPIWZxxFWWJLxWGp2pJp56icQQTa3nTVCMjn\n5wObyWSeI55bEYQkhKbPFB9UDseYqqpgpFiEvdbxSueICFZg3kZKKxTW4M32HJlM5uXluRXBtXaJ\nCJPC4ozBCPjSsiwtRmBcKD4q8zZyUgcerDqmRbLpL60QNXeTZDIvO8+tCMLafh9KZygsBIVFpzxo\nms0Mkmsjh6HjeNVxVHeY/vhMJpOB51wEE4IhzSfuusjhsuXOomVmOkzVcakK0Art0rJoLWNvCGtP\nrTyRLpN56XkBRDAlOowBUFZd4KhuUWpqXSG6gq4gtiPaUNFFJapuOkey/mUyLzcvhAhu7PN1e4eI\nIv2dVlI3SemS87SSls7r7pH1vmAWxEzm5ePFEEHdjucsXTJWmJgSX1qoSlx0jF3B5c4z8hYhWXE5\nw2YOyfo8mUzm5eKFEEFIAmYFJt6iYxi5gnLsCGOHj8qlRnEttCHZ8x/XHZUzjLzp3WUymczLyAsh\nguulrDHCpHCMvEvF0QgfrpKNlneG64Vwf9lyfxU4aQJGHKVLFlu5XCaTeTl5IURwjRHB9j4JdRc5\nWnactIEDD9cq4UoFJyHSNJEHdcRbYbyums6Oq5nMS8kLJYIoiEnL4rkqh3XLreOGxndMmw7bdlA7\nusbTBEcblNBPocval8m8nLxYIjggqtIGZdVFGlqiqxFbI12FqAH85rjYPydnijOZl48XSgTXpTKx\n/8OI4G3aDzTWotZj1TPCMTUWbw2q0IZsrJDJvKy8UCIIvQhqKn2ZFJbL44KZ9/jRmDhSitay11m0\nSz3ECjQh4o2hsGmkZy6VyWReHnYugoOONT6PypThoPZZ4SisYewMUlpOCos4ZRaUMijzJjBvI8s2\nMikUZ2yaXkcWwkzmZeGZRYLrpSt8PoLojGC9UPm05J23ynHbMfHC1FtmBXQxcrgKLNqIM46xN0i/\nIM4RYSbzcnAhlsNrQVzr4NMQREGwJu3zPagD91Yt95cdr1aKHynTAlhB08AyGMZB6WJqr8val8m8\nPOxcBHUgMZ9nEiLNIElmql1UjleB2ycN09ByXTtMCNAUxFAQ1BNUUdWcGc5kXjKeaSR4VhA3/gdP\nISrU/vwRiOj2tVSRGJDYYdThRHH9IjgodLF3pSGXy2QyLwPPLDEyRADk/AhxuC/3SQRR+9dKFvxr\nYwVLVXpMYdAiUpiCmS1AHaVLe4chKtKrYE6QZDIvPrtfDj9C1Db3CwxjRPkMi+a14Dor7JUOI8Le\naIQfO2JlqVrlcqOUndJFpY1KiIHKGYyYzyV7nclkLhbPdjm8XqGSBFE2fzzd1ymtwVWGWenxBpZR\n6FZpz3DkDd5FjlaBeROSkYKAd4LJhdOZzAvPhcgOw9YTEE6bnSK60UXVJEqfNEIzIjhJz1u0gQer\nSBsjVyrL1cowtnBCoO4CncpmaYwM5xs/hTeZyWQuHDsXwY2BKY8XFj1zeyOMn+I1rQhWhDq03Fk0\nHK1azDiy10UKC7oydK2lFUeImkUvk3lJ2LkIOjMcnD74GYjOo6PCz/rqSheUZRc5aTpq2xFth3hF\nQoFoWpMr6zkkki34M5kXnGcggrIVwd7B5awonkeKBE9lT56YZKqwbU8xIngjWGsQkxIg3lpGhcNg\nsCIEVbp+Pon5FK+ZyWSeD3YugqWVTfnK6R/d3H5oqSynIzHZ6tkT7Q8Oy2WcSaUyVkrGY7AjAS9U\nrWG/Naw0OU13WysabH/NOQzMZF48noEImlS8rBA1RVxRIUQ2Bqdrj7+ze4Lop18WrwW0tIaDyjMr\nHOPSEkpLbQ3OKgc++Q8u2siiC1gRxt7ibY4CM5kXlZ2LYGEFVSHSi59CjBAkCWHX78WFmFyf1+K3\nEcS+jmYrS59MoAoreGs3S/CjOnIikZE3TLzBG6hDpO4iViRdLwYDaB5Cksm8cOxcBK2skw/SJyCE\nIEnwgoBVoYtJ5KTv4NCn6OhizXbE5uGq43DVEULglZEwGQkFQKs0nWKMYRztJjWdNTCTefHYuQiu\n9/DWdX8iYKxgVVM0qJp6dwXaqH1P77aDZOg4MyyZedL9wbWxgmpa+t5ftjRty54GIjFNnmsdXWcx\n1tPFSFD9fMwPM5nMM2fnIriJ8iQNS18LmSENSTJIb3OfWuZakl++9svnswmTT7U07s8TYt8uFyIh\nBAgdxhg8hsI6sIKsx3EOjBVyqUwm8+KwcxGctzGVnRjZzPVI4iIbkTGGgVD2QVhMohU534ThSUnZ\n5ySo1giVM3hxlIVBCocxhkIcU/EESRb8XdR1ncymxCeTybwY7FwEj+qIkyQm6ad3ge5HZZpe9KyR\nrQD2CGkoEpxeFhuefLW6fp6ilM5wMPKgjvHYIWOHWMO4jVxplDqkiXXLNlLY1E7n1q+lORrMZF4E\ndi6Ct+dNEhRrKK3B25SB9UbwJmVvTb9MdgbWS9yUT+lNT8NAzJTTIzPl4f3BcyM3hZGzeJNCPGOE\n405wUREM0wJsF3gQI8suENXgjCDOoDkUzGReGHYugt9/sKLso6rSGsbeMvaWyhsqK3SaylS8SaYH\ntg/1jKQbUSPRaCqrgVMJk0+yPygieIHCCFFT1PegbvFGuFxZLo0sVgwndWDVBgTo4kAAc6o4k3kh\n2LkIfjSv8dZQGKGwhlEvgpPCMvGGaWEZ+ySQpUtRYfohLUH7wuUW0EEtYVAwCvTL6iGPigpNvyxv\nA5w0LR/NWxwB18J+UCQKsRG6KHTREFWz7mUyLxi7T4zUAWsixqQlcNkGSmeYNJZZ6aiDUgfLxCtB\nLYXto0IDaoS4ifDWnSbbPcIIiAqRTR7joQ4TOXu7N0yoQzJVcNpS24i6CGqQ6EH9qVa/9SlzMJjJ\nPP/sXARHhU1tcVFpghJUqbvIqv+Zt4G9wnFQOYLC2BvUmrRPKEkQDYISN212IeopMbS9Aq5LbYbW\nC0OSsKVrgPV+oiR7fRGMGJw1lN7g+5N1UXtrru0YgCyEmczzy85F8PqkoI0p69rFSNfX6tVdJMZU\nwNx06f6gShssobAYMX1ZjWAtRDVEG/seZN0YMCRkkDZ+2Idws4/Yt+2pCoUxzEqHF0tZCVIJDsPE\nGi470zeNCHUX8TYJpeuLBrMIZjLPLzsXwTcORqy62Ed/YRMBNiH2CYqINknQ6hBZVa6Pvjylg9Km\nZIa32zRIVKE1QuyFMwYlmr7c2spwhtNGDYeehQiMC4vpy3bK0tIUFhAKGzkIyqqNG3OFSg3G2/Th\n5TVxJvNcs3MRfHVW9F0aKSN70nScNIF5G1i2kWUbUIVVF+hi7N2sBGcMU7W9FyC9ucG6BS4ti1ek\n6LINitGUJREREE3tecB5WWMDVM5Q9RFfp3B3GfBWqKxwUFoeoCw7Zd4GjEBhtxkYeYq9zZlMZrfs\nXAQvV5aud4vpgrJXWupeXB7UHUerlkUbaUOk7pRj6RBJ+35NcCgO13dyuL7rJMQ+g6u6ObeJghHF\nSO8JCDBo0xuauEo/g8QYaILyoO54UAcmFq6PDJeMYRkjGgJ1F6msIWhEN+mXTCbzvLJzEaxMRE2y\npwoORjFFhCMvVE6YeMNxE3oh6oiqzJtAE5I3tDNCYYTSQgU4gcJAdP05g7JS6GJMA911bZ8l28QH\nD3ei2L5zpe0iiyZwZ97Q2MAeJDub1qCdIajdeCDm4C+Tef7ZvYuMrotcIgbBGqE0hlG/HJ0UllEd\nNlngVRdoYmRRaxJAm5bGswKcCK7vODFGEEnL4ToITYAuQhsjQdeegKk9b+hks9HCfkndaXKvqUOk\n0UDoInQREx0Wn65LkidiiNuWvWyskMk8n+x++Lrx61uIKhawohgF64TS2pT4MMLIW45WHUd1SxM6\nll3k3rJNri9jjzee0lusKk4UY5XgDarCvI0s2sCyDcSYBNANfPmt2YaCIlvLf0gRYeUMpXU4B3jF\nGctIHHvGUliD6ctl7JlymUwm83yxcxE86dJOWnKPMZtBRkYEZ8H2j1mTWuusScXMTVBUlZMmCZsR\noXAGZ6DQjkJbnCqlOoJzhCisurQshnXyIpkmFE62y+L+uraGDMLIWS6NCiZeKMYOnXhcUCatEhsl\naCrhWXWBcl3D2G8P5gRJJvN8sXsRbOKm2DhlfdOydugiUxhBHHhjWRfiGRFO6sBx0zHvIiMfqFzA\nAlNpsdLiDYysxVsDotRBOOkNVOsuJVuCWsCesvGib4dTUnJkVjrGhcMZECPcb5K5qzfCpUo5bqAJ\nKZudjB7srj/GTCbzlNi5CC7buPUQ7EXQGkn7e6ZPYJAeKy3MipQJHnnHR6ahjcna6qTpAKXuDDdK\noagchYVCIkWsCUBtoakcy74ucdHGlAk2grfmjDlDwvTF2EIqpF60kcOVMvGG/dKyV1qaEDmKyrIL\nONNP0LNJCHO5TCbzfLFzEazD1ll6HRFak5Ie3qRSl6HXYGlT24czStNZVl0SoS4qdxct88bgTMls\nUjFyiqnn2HZOGQ0zGcFozGEd6UJL03sDWhEgUNq0pPZ9ILduMrH9a580gaNVx91ly5UCygCXVJBG\nCR20AUI0Z5xsMpnM88RjRVBE/j7wS8BtVf2pRxzz94BfBObAf6aqX/24c67336IqUVLE1UXoDHRG\n8VbwmswSBPB9+cpeaenUYwTuLlru1x0ndcfEG8aFgy4y6Wom3YrCOma+pCzTZt2qld4OK7Js+9f2\nFhm41Kwt/9d/x75o+2jVUcVAa1PyxQSLVYcViyCnjBUymczzxZNEgv8Q+FXgN857UER+EXhbVX9U\nRH4G+DXgy4862TAru2ZthUXcjuEMJgnLer/QG2Fa2H6ZbGiDcrTqWLWBu4sOkZq6VK6b1NtboZSx\noeiOqdVwYmFW2tSnHCN1E3uzBdlEfkLqLll3oETVzdJ2WFe4HsVZYrAmZZZD3A6Ikt6wIetiJnPx\neawIqupvi8gbH3PI36AXSFX9iojsi8gNVb193sFrr791NLi+DVurKu0TFRFNvoMiGAVvhbEk4Tlp\nHPPWE4F5Fzi+v+CwFNh37I8PKE2LtEvM8j4jKi7bGWY243CVulLmbcCJ9HV/qaC6sAbRfhRo7Ie1\nu2SsMCodrhRimeoUx9YQuzSISTV1mvg+kZIH02Uyzw9PY0/wNeB7g7+/39/3CBHc2k8p26zsMJkQ\nFbrQC5HdFiEbSd0irhAujzxdTFHXrQcrbh0tOXSGqZ9xeW+C0QXF8pDi+EOKap/9WUU1caDKsu14\nUCtNSLWECoyc2dT6rQc+WRGmhUt9y4XBVo6Vt4goYxuxNrnerELqVxZvcL3lV06OZDLPBztPjPxP\n/8OvbG7/9M/+PD/9s3/lIUEEtlPl1j6BmpbEaR4J7JcOa5LP33zVogoPVh3fP67x3nJsV1xvW64p\neO3wumKiJywFjq0y98klZtUlX0IAa0y/J9iX6xjDuEgOM0KK9u4suk0bXmmTNf+8Ta1/tp+RYmUb\nDWYxzGSePl/5nX/F7/3ubz+Vc8mTDA3ql8P/9LzEiIj8GvD/qOpv9n9/A/ir5y2HRUS/cevo3NdQ\nPR0Vrq8qZWtTVJZmk6RoMLXEJbODb91Z8M07c+4u2uQ3aITXR4GfmHX8xKxjJAFi+vnI7HFH9rgr\nUxZtYNEEgirTwjItHKPeQNVZs3G09kZYdpGTJjJvAlcqw5XKMPHC7UXg1qJj0cFeYZmVqaNk6FmY\nyWQ+X75HllpVAAAgAElEQVR0cx9V/VQbUU8aCZ71JR3yW8DfAn5TRL4MHD5qP/BxL7FxeBE9JYYh\nQkCTLVZIR1sDI2uwxvHqXkkHFIdL3r+74P17S473Sq5MD3jnygF28RFy53248z3Kgy9w5VLFaFrw\n4byhCYG6Sb3Ctksnj84gojiTulrSHBLlaNXy4bzBBphFZVaANAaCQdVsotese5nM88OTlMj8I+AX\ngCsi8j7wy0ABqKr+uqr+MxH56yLybVKJzN/8NBdyKpmgA5++gaSEqLQACqVLbXaVNVyfFFTOUopw\n9KChqzvuHgvvHTZcvtNwva45OF6xv5xTzVaMTMusirQtzBuh7lJyo+5Si50R8NZsMsRBkxFDMlUN\nNG0k+og47U0VPB5JJgqaRHu9rzi07cpkMhePJ8kO/ydPcMzfftIXPK+e7qy11fq+zfn730G1nzAn\niJhN6cyksJTO0HSBG3sFV/YrFm3gm3fn3Fo2/Dk/5ycLx6VL1zAWZHEHCUvGOuHATNCqoo2Rtp97\n4kzEdAFIZgyqphe2lEF21iBWwYJTxwiPRosflMukTpSz7ySTyVw0dm+l1f8eauHWvKA/5hwFXEdT\n6+FKbVRMRxIigdLB9YnjnetjghXeu7Pgg7sLvvH+Ic2NghtvXUfeegsefEQ8+ohw+3uM9l/nysHr\nlJOCw1XH4apL0WZQrKS+YNVtMXRhDfuVY1Q57MihlaXoYNqCDWnYUxcVFUUkqWDfmpzJZC4ouxfB\ngZff2TrBRzxjcPy2frDt9wZVlUoCViITIl+ceSaFZ6zCyf0Vf3q44t6k4HaY8oG/zigcU82XFEcf\nUY0vU/rIZOJBlboLLFrdzjohzS/RZE7NtEzzkStvaDE8CKkbpfKp9W/eBlad9m7Viu9rBnOmOJO5\nuOxcBIcMBfHjjlk/Htf3wXagkgrGGbwxeCdc9cL+GJbLlnf3CsZTz/0Q+eqHJxwa4UebBe80hjeK\nEaItZnUPscKk80xtgdCPBCW12CXzVKWwyfR1nfk9aZPoTQvDrDBUVqhDTNGg6mZpnMlkLjbPVAQ/\nKcOldNxkG5Q2GmpjKFXxBEo6rleRN65V/HvNZT46rvnmD475/W/f5a+9Juy9do13blxDQ4D5Pezi\nmPHoKleqa5Rl2ZuxRmKfllnXERpJfcyHdeCwb9l7dSxUY2HkBFqlC9Ahyc1aDSpZDDOZi8xzJYKn\n5mWyLZ+pg6Zh7BIpY02lDZcL5adujrlyMOWP3j3kD791h2+/d8T9S68wn71C88VLyO3vwoffRVd3\nKK8UyOQKrkxeg0oydViTluHpZ9kFDlctx6uWPQOtE8ZqkGgw2D5PLETADPc6yZniTOai8VyJ4Hpp\nvBaSdaa5Xc8bFsUYQ2E805HwxshzA0dz3PC9sccY4cNlyx/fWWK/57jxYMH1ec2BtlQERjZQekWi\nENVQB90kRdajPbs+g1yHSB0ibTQogpiUNS6dQTS5TccI0Whyns7ziTOZC8lzJYJnWQuLahqQJMZQ\n2IrOGzyBceyYhiXXx4HrVwtuvLHPXeBfvnuff3f7kJ/bW/KzszFX9iuwIMtDXOgopWJsR1iTXGe6\nXm2jQt2lnmaDUFiL9R6KAkpHYWBswIV0banGMCn3ZlA7OUGSyVwknlMRlMH+4La3OERl1UZUodSW\nsltQdnNuVvCX3pwxvbLHdz54wJ99cMTXvjvnR358j+VrX0ReK5F2Cc0c2zUU5SXGVYWVZNnVRjZi\n2IWIFZj2tYneeWo8c7VEUUqniGyjRavJcsvZj0v/ZDKZZ8VzKYKnymxI+3TS/66j0rYhtdmpUFjP\ntT1HtVfxtpb8i1XD3Q/u8P7hMYf1AR/ZAz4qplTNbar6EENN4SeITYYKK+nDTNblOcIIS+EMURVn\nDCdtpNVks19ZgxOli2lWslXpS2V695wcBWYyF4rnUgTPkmaS9G7VMbUXW2Mp/ZjgxlQGKoRrAt/a\nK7i0P6HYa7jdGf7N7QVzOt5kxVuiHFTJAMF5szGAVU3Fz7bvJ8b1oqvJwOHeSY0AN6cF40lB6ZIF\nfxsiISYT2KjbBux+MF1OkmQyF4AXQARPGy+sTQw6hBqLwVC0S4pmjm/n7NuG126WfNFc4UGE333v\nPt+6Db/wiuHqqzMuTwuis5h2iZGAxeJM+pisgF2LmaTscRMidxctXYhMHTASrDFo34ZnJS2lo+p2\noNN5PYGZTOaZ8NyL4NkWu3WkFTVNtmuCMg4BQoePHW9crviF6T5vfsHwR392jz/+7j1+8GHL/b2r\nLKtrhFmFxBYJDYaANRXWWMQKUWXzgUn/ejEqyybQdIG6tcTgQBRBN70uUZNgDsd8Zu3LZC4Gz70I\nnmU9xU574UnRmKV0I2LhuDYp2JOSVxvD/M6C7yAcrQJHjeGjruBS8IxCyyi0YCKCx5h0YrPORrOe\nUZyGQo18MnMovEdcgViDtS3ebmeUdL2pgpxxihj+mfcLM5nd84KJ4OmlcaDfJ7QFXZFKZ1yo8V2N\niSsuFTX7M8vdruBeiHz97oKGlteLli8UAe/tZj4y9OIn2709JzApLFfGBW2IlM7SxvTaxhimhaWL\n6RqWbaSwAs5g7eA6yVFhJvMseaFE8GxUtRaaNirzNu3RVe2KUXuCr+dcnShvvz4l7isLEX7/+0fc\nOjL85eueKzc83lnEGIxZjwjlVILDGmG/TDNIuqiUVlh1ignp8WlhkwdhE5m3gZE3GIGyH9S+uW6y\nEGYyz4oXSgSHrDPG0u+/1SHNASEoXoXSOX7kakm5X/LGifK1D4742gdHrA6Ft8YH1K9O0KLCGIcT\nQ4xKCKn+LwV0aT7x2FsqZ+hiivaOm2TBNS4MY29QVU5Qll1AJNlxhYHiDYU7L40zmd3zworg0IIr\nqm5s7ztbErxF7JhrUrAvnoMHHffuLPhaHZjXcNJa5jJiz1ZYA5X0g9hRmthbZCkYA4UBJyZFm03k\nQd0iQOUKxi6JoABNF3FG6GIk9PelNsBcRJ3JPEteWBEcWnDB1rg12oLWVzQObGiYhhUHumRmGkal\nUKthHuGH8w7ravZsZM9FCIIGQ4gWw3okZ5oz7PqMiapSdxFB0dDiVHEasZrMD9fW+11vv285vQzO\nUWEms3teWBE8y3p53MbISavUnTLuFoy7OUWz4vI48vrNCfc7wwOUP/romAcry49MhekUBI+oR8Qm\nU1dVYhCMMXgDimB6QTQoVgOmC9gAXpTSGYwRgirLNuCt9NPsthHrqevNAWImsxNeChEU+vIWUl1f\nHZWWiItQYhgVnrevjbB7nu+dRH5wtOLfvn/EvamherXkjVGJdQ4vqfOki5oiOk02/FZSD4gzhlnp\nsEQqn5bLzsLIw74KoV+SL7uYnGcwWKPbxIjKueKXo8JM5vNj5yJ43qClIZ93RLS24FKEYEuitRRV\n5JWJZx/P7N6Kw8MlH945wdeeB5crgptQFJ6q70A5aSLLJhmvBk2JDm/SvOP9yuGMUHlBvMFHZWYV\n45V5E5i3gUUbNqM8h9lmRAezVrZvPEeFmcznx85FUB9RDLJOEJx99Gzks04ofBKG+4NDkYmupPUj\nrIOxRqYamRfKRNLGXd0JS/WcmBEza/ECIwNtpxxryhSbNhm6ljYtbUcu9RwrMA8p3iusUFhAI6sO\n2qAUBqI9PWh+01fch4ZZ/DKZz58Lsxx+lDgOEWST4Hj4sScTjaExax2USKBuI6OwYBSW+OWSkW2Y\n7XncyHMc4f2jhmulcuA6DmyHaS0mWkQsXYwsWyVEYVpYvBWsQBOURat4o0wsTBzUJvJAthf/UKG0\nbh1xjOhGCHUwh3kztImzT85kMp+GC7Uc3tjQP0LMzhPKYQQ5FMfhKU5lXdn27gZVQgedRlwMVLFj\n4uHNKyP+kh+zCIZ5jPzrHxzx1jjyY9PI5UmkZMTUjojW0fQO0yEqpTW9wWsawnRcR0a0ONcytQ22\ns7joKKzD2u0+5en3uI0OjfT2/AMz1s17yQKYyTwVLkwkCNvv9aM2/88TyUdHkNuDhjONh/bOcfNa\nQmcKooOxV96pHJeueP7sfs3Xbx3zB7ePObls2H+14J1pkZa91qBiOK6VJqTe4DZGmpB8Bk/qwP1V\nR0fD1NfgaryWjIxjWqR9Q2sMj2IthOvI8GzJT06WZDJPhwslgo/jrEh+3BJ4KI7rZfQjEYMWY2Ix\noTDCNYVLCqu641ttx52PTrjnxyzjDB0fUBiPNY5SLRpTbWAavJ4yv11UFl1k2UUKo3QqqFictYzF\nEsSdmpFy9v09/Pdg45Dt0jiTyXx2di6C5y0B4eFo5kmCm/P2Bz9JtDg0MFgFRRsopaOMDSOtGTVz\nirBCUDocC1NxKFMK7fBdQ6WBUoVShK7P9K467f0Dk8GCdQ4pPDpSrHgqHGBoQuo+6eJ2ib5+T2tx\n3GSOdRsJDo/NZDKfnZ2L4CPbxETPHPcww0zqoxhGi4/bY1w/pqRERhuUlg6JK8q4YETDpanl1Rsz\nymnJ/U755t0lV23NVVuzZ1oKHTGyFdG4PkkSewNVmBSOUWGhdCxLixWhALympIm2kYBuPAaFNON4\naKjw8PvVzRJ++B6Gn1Emk3lyLkwkeHYP79zvsuhD4rj+0p93/McK4iCyWr9e0OQcHTCocezPLF+y\nhuKS5biJ3K9b/sV37/ITs8iPzwL7U8HbwMhCVMO8VboYUaByhspZrBHqoHw47xhbmHiYWGjZ2uwb\nAbseyxnTlUc9I4TrZAmDfULysjiT+azsXAQfNXRtWC7yqMhGVQbCoKeOXYvZMIwaCuOjki7r3HLU\nfj6JsUQ3JvoRs5Hw9lS4EQx//MNj/uD9Q/7NB0f410te8SVvHZQUrmBiHYIhRKUxkdgXT5dWaELk\nQRM4bgJXC+VGBbMSJFpQC1hEwJrB+47rXuReEwfXn2YgD+7JUWEm85nYvQia0xEfpO/xKbODR9S6\nqAz3AeWhc2yW1HL6qU+yjF4TEWqxnIihcMJ+IVwx8MM7gjQtR/cXLK5XNG6GzvZxYjFisTEVQTfB\n0Cl4mzpCFq1yVHfcPqmRkbBne+sZNQOVEgTB0u/9GcVoMls4e+3rPcPh/mAOBjOZT8/u9wSfxClF\nTgulDv5I7WWD5w0SCuui4oee95hltGraqTRmW0TdxcCIlomuKLTGLhc4rXHe0hrHA3XcbR0VHSU1\nTiOlWkbW0GL6MZvJxqsJkUUbqEubMsXGJhGUZMm1ubbN3qAkAUQxkVNL4/U1r++Ig/s/ruUwR4WZ\nzPlcmBKZj01enDnoVPTX/yFsI8WzzztvGf2o10l2WOvleaQgAC2TQrh5qeJHO081KbhTB75xZ84N\n33LTt4wteEoKUyVxIw2DV9J+X+Us3heYskKrEtqIaSOmn2kcFIipgNoaMGZrwR80PT4sqTktikOR\nz3FhJvNJ2H3v8CeMSB4njuuvf1pSbwcbPVYQZSCe68dYLzc11f0ZIdiCaA1XLik/UQr7V+Fo2XHr\npOaDwzk/uQ/jA5hOHM4WFDbNmUvF0+n6Rt5yeewZeUtQOGqS1T/93mGaRqdEgcIIBgEzuP7I5iqH\nEe46IoxnSmiyDGYyT87ORTAMwpmzk9fgyb7A59pNrf8dtJidF/1FTi+dzwoi9MkbEYz1aFkRSsve\nSHl7FnmlVb7y3n2++YMF796ZM/3iiLemI/AFzjgqa9CYHKSbkKz290rHrHD4fmzn4Srgei/B0qTa\nwjZoWv46g7Npf9CKIKYXtri+/nMyx+sbg4SJDP49+5nlpXEms2X3InhqHTcUxP73meOH5YMiT25G\nf8o5ZtAxYtYRFUMhkW2/bv9ksenvZacEDZQWphYue/g6AVYN8wc1q25KW+yhkz2cCAYDXdoHjL2D\ndGlNGsIUlEUbuLuMXKoclyrHyBnq0FGHSBsUI4qzYKxJWePhHmqfLAn953g2gl3vLz70I1n4MplH\n8QxE8PxvowySGkOhO21+oJsOiu3j25jnUUvn8wVxaGa6LZoeXl0kdZIkt5kax4qRrnDNAkeHLRyN\nWO63hru1UplIZVp8VCy9AYJJVlqVMyy7jpMmcG/ZUhi4VFnKfvxm1xsxFEEoezE0aVWcSmLWrg9C\nP/aOTe/z2eXx9gPcbp4Oo+6hIUMm87JzgVxktpv7p8RsKHh9FvesMKbfj44az2akz+4lKqc7jU1/\nvhA1JSWiUhCJEoHA3tjxhWsTFqXiK8+tRYu9M+dmFblZRpwxFOoZOQ+9eapI+g+g7iLzJtB0gRi6\ntMyNAe2XuV1MPoUi/UAn0xfQSJpJAvThqqbIkDO1lP0/ay3M/oSZzMdzoRIj+tCNs19cfWjjfy2M\nww6QtbydrsI7fb5TRdZnBXGQVNEIUQSMA1uh1vPKFfgLY7i5UO4tG947XPKDwxN+6rJh/5Jlf+Rx\nGErrCf0LtVH7bLGefk+Saie9MbQmSVcbQfp9AxHZ1AJur1kwpm+xi/2Qebb/wQy3B1UGd+TIL5N5\niAtTIvMozjNWkDP3bKLAU0I4XF4PxLM/33nL6LVO6Pp1Ja1HRcC6AvUVXWG4VCjjmXJt1fF77x3y\nx3cWHC9rrpgx7+yNObAOL45KLI1K6kkOqczFW8O0sJTeIdYRxeGdMCkixiQ3mlUX6KKganobfk4J\noRm+PwP0hg2Dt/fQ7U3C5Jxl8XmfcybzsvBMi6UfxeO+kI9ynBlGiaeXyXLKlVmGJ5FtsfTg0V5s\nkrIaSVHWolNaDZQEpgS8NFTtkna+ZL7oqLs9gp8iozEFBqOySZIsuzR2c+xtSpS4NKv4pEltdpUz\nGIF5G1g1StcPeff9XqARwZrh9Wnfe5yWzKHfHNT1z+CzWv/HcfYnk8lcKAOFLWc18KG/HyGSZ+d1\nbFaCa6v6wZ7iRhA1CebmeWdUemNxpclppg4KcYkPC2y9wHc1hVVs4amN5TBY9ltDaZXSKB0RjZG6\ni1gRxt7irUGAuoM2BEqXEifeGlYh0mlEFXzoZ5MI2L4G0qyvv/9tpI+CVZAIgaFZ7JnPTPobm0RR\njgozmWdgoLD9tj0yoht+ITmn5e3M70cxjIZkfa51kkC3MjAsNB6edSiIUbdJkqCKasSKcnnqeeum\nZ9JYfOm4NW8xKFfLyLUiYqLBRoMXi+1LZUoLJ21KkLRRuTZ2zLzBWsNxXyzd9s7XQVPrnEpyuNF+\nD3Gd+Blu9cmZeseHTFsHUSHn/M5kXkae+XJ4WMN2/pdRtqUf619nvuwPPX4OQ0Fc3zPcP3xYDDk/\nJBJAHGornLN84ZrF7Fk+rIXjuuO9oyWHc+XP7QmX9sA5T2lKxkWy1UpF0tDWkaO6Y9F0jE2g84HC\nGaxGjMTNfxZRT9dWan9dmyWubCPCdUZHFIjpkzvbd7x+Wyrbz2D7xvpbg+g3k3nReaaJkbOZ2vPY\neAEO9/MGX9jTbXLbL/UnMWDdZJLXS+ZTCZbB9SJJnPwIdWNwcC3CXoDL84Y/+N4R3751Qklg77WS\nt6YVE2MoxTDFEHtB76LSdCkBsmgDTWeIISKWNNbTWhqTouao66THQPL7yDDVEW6jQjsUcaObLE/U\nhwVt6ERzei81k3m5uPDZ4fME8pTT1jCzOzBKGDrHfJwgbjLBg3s2nRcyEAfZimMblZMO6gAlLVPt\nqLsVslpwfLxkYYRFNya4MaYoqMRhxdCE5FDThCS6RgRvDcY5xJeIt3iNjDVigxJJS2+i4oxgRIj0\nA9s1iZ4RTV6EsHGo3n4iW4GL/Z7io/9jGD6SpTDz8nCh6gTXPC6D/CSOM8Pe4PPmjTzOkXoTGbL+\n2YpqG1KSxGhkFuf4MMculpi2xZqIuJLWFSzMiIn1WFEmAkRltR7aLslYwRmDt5ZWhVVMUV3lDSKa\nRnmGZKygcvYat2UxxLX4yZnlcZ/46X/imahwYGd45rPb3p/lMPOicyEjwccJ5ce1xz1030PnPsdF\n5jGCuL6m9bmCpoSFxEjVJ0m8hWuzgrevlTSmwHvHnVVANLBnO/ZMgGCIwdIGizPCXuEAxRmz6VG2\nktrsUsSZZiIbBLveudTTbtPDH9DkQNO/L3Nmr3N71OnCaj0l9oPjh1mUTOYF5UKK4ON45FxiOf/2\n2eXzeukM5ydahs/bCszm3l440g5kdCXBGIoi8JozMBMOWyHEyLfuzrlfRN4cRabjgKHEkSz518tX\ngDZGTpqAEdgrLSNvMSIs2pjMVUXRyKBvbn0la8PVzbABtBfCTRmQpKclFxrZLI3hjD/hqahw+Hlk\nFcy82DyXIvgonkQcH3qM04J4qot4sBzcCGKvioIgRlA/IpRTCiPcGEX2g/KDBzVfv33C1z485mal\nTK45vlg6rFO8CJWRjTFCUFi0kfvLDkjO+5dKgzVpZOc2HbJuuUvP1TORXGqRU6wIavo9R7ZF4Kav\nDTqb7NFzznU2YbLdF33055jJPK+8UCL4KM5LpAw5O6zp4ehveGw6eG1x1QQ4qgOFEUYmctkrJ6Zl\ntVzx/kdzuqnjnf2C4MfYomAkHitmsz/X9BZa8zYQY+RyoWiX9gxFFbPuEhHZZHRNfx3DWSyRdUCr\naNymibdtdkkIRRh48uumXObRWwEPR4XZmivzIvFSiOCQTU3iI5fLMoj+9JQgnkm8opoGqNcBPAHD\nipIl5mSFLleEtqXVgtaOaKoDRt5SClQoq6Csur6OT5X4/7f3JjGybOl93+8750REDnWr7r1v7H7d\nr5u0mmzIGwIGWpZFgQa8MOmFuREgEwYEeUUIFuyVIcMQYO8MwxtblgxZAL2wAUMyvJAIyDK0sA2J\nJNwkJNEti91kk2z18F73G+5QQw4xnPN5ceJERGZlVtWd6lbdiv9F3sohMjIyMuKL/zf9P03Z25iC\nEWPIrDKxsWg6te5ZABNb+lQ4J7DqU0an/b+3oYNYYZrbSRxs1yVLttaVCsu7x/QXipEVjnhTcOeM\nYMJl7LB9Jf7fuprdWwZWsfZK1QQq9UyNRyWQGbg3dbxzNOVgllMi/HjR8CCrmZuGufEQHEEtTTAY\nE2eQALgsA5sjVsiywFxDp0DTBAWjiBqkTf+2+gkbm5WMmqassonfIYkwGFqJrjZrnFLhSXFmmAza\nzBzvuIKMGHHLcWeN4BAXxRLbypM+q6y6NfqTbiF1OeoMUzPhKxqoZkqtQgN8+/MFXywCX5p6DiYB\nKAgUBDXkxrSZ4lggvW4VZ4BYLtMEliFqEaoVjAx5WUqHDL5PurVMNpXQKL1LnNxqoDOSrVbr1qxn\nBhT4/JOjPRxx2zEawQuwbRxTaqFLNhBVXgyCNRaTT/CFoVD4cBJ4cN/z8UnJ9x4v+fZPTjk7EmZv\nOT4o+jSvEaKb7GxrpGKiJLNCYQ1zFxMhlQ9d90hifqmcRelnMW0bwyFLtGy69qYtBDyXAOky4Gz8\n3VhurCEc8YZgNILPiG3tQSRlXqEKykkZyIwwzwwPJpbFquJkUfIvPjrGNVO+cpCh2QzB4dSSD7pc\nmqCs6sBJEyis8HDiOMwdQWO3SBq7uYu4ysZGbS4ThnHC5BYbOnOekifd8u17dI+hiwbxXKR0xIhb\nidEIvgBE4lS46C5HDcCqUXIarJRkUuEWS8JiwXpZUTYzKjOhnh4hAfKgmNDOFA7gQ2DVBB6tagqj\nTGh4YA0mgAltb3Br6Lre30HSo4vtDbZxGCdMryqC3UpuDA3hRgnNDnaZ1rfJHkfXeMTtxGgEXxYk\nzQcJqNY0lMAS5ysmDu7Nc4yzLDx8svBMncRbblg3Som2hjSwrj0YpfGgITK2zAiFi4wQaSW2uvIX\n2diOYRZ3mO0dGsIkwpCKqrs5Jqn8pU2apJ7jwCZ0cK//9NEKjrh9GI3gFraZz1ULhZOic4QB41Ap\nODhwfPU95ayAg0lGCfzh4wXvTYT3JsIkF1BLCBbfaggCqBjUOshyDEKGMmNgyAZxwZT0UNks/N78\nfxcrbKfZsZsVJmPabdMOV3xYaD2ywhG3EaMR3MKwS+RZEbs1DM5kSJah2QGHE+VnpoEHD5RHq5rH\ny5pv/eSUrx9ZZkeG+2JRMhpyvBqMCIUz5NaAsVRqo9ttlImLs0rqENpyGenikbu/zH5WqJra57qq\nyIHxko3Yp9AXYvf/77pYjHHCEbcP5vJFQER+UUS+IyJ/ICJ/ZcfrvyAiT0Xkn7a3v/ryN/VmQwBr\nILeCMYYKx0nICDbnrXtTvv7OjKPc8ORszbd+cMwPHi05XftYuSyRd1kRps5wVGTMM4uqcFIHlk3k\nYrOsldunnV6n0YVOnSHDW6c12P3bNE1RtTq61fG2qTCTGKZNNxM3Ncl1bcQNt24B7WsV95QfjRhx\nU3ApExQRA/x14N8CPgZ+R0T+nqp+Z2vRf6Sq/+4r2MZbge1Ok7IJlA1MxHMgFQdS41ZLlmcLPjlZ\n8/TtCWtToMUB4gXnoRAwxjFxSqNKCMrTVcPUGe4Vlrk1ND76vnU7iClcZWgL7GeF2meCU/3j5hwT\nOecaJ9GGIRK7HJbQjBhxG3AVd/gbwHdV9fsAIvK3gV8Gto3gnTjsd8lubcfAlMjUaq+olkx0gbDE\nljUuNDhjCBjWwXIWHAhkDqy27q6BRRU4rRuO157DwlDYnHzSj99Mn6NsOqEyeG5TAGIgODswYcPn\nzsUKBwmT4Y8rEt3jZECH6xmuL22TDsQQx1jhiJuGqxjBD4AfDh7/iGgYt/GnReR3gY+A/0RVf2/X\nyvadBLfFbdoWW4DN7hGhnQsSYsmMJyY5EMvBVPjCg4yvAUfznHWATxY1h7nhMDdkTliLUHplJYGq\nUY7XNQZDWQjqQULqLgZUYmlNKm+RbcO370uwU4mmS2wMKqNlsE7ZeJ9E0YZznPA8K3yROOuIEa8a\nLysx8k+AD1V1KSK/BPxd4GeeZQW3ceTjRQZRRLBWMSYDZwluyn2jfG2iTO/HspOnZcPi0zM+PBDy\nuaEoLAaHtD+LV6UKgTrE+SQYh7FKZgMTFzbKZRBwrUu+s5h6yzjuSm6kMpo0tQ4d1A4OWGH3hcPm\n+4DdlwgAACAASURBVHZ97jCVIoP/R4y4KbiKEfwI+HDw+Evtcx1U9Wxw/x+IyH8vIg9V9fH2yv67\n//q/7O5/49/4ef7Un/mz5z7wNhrEIYyAQ7DWYq3DZ4aVE6ZW+cpUee9I+e7nC777+YJHizX27Yx3\nTMaRy1vqZBEB1worZDYmT+rWpbVGyF2U1gqhzfKa2FNsr+RvJuvd6wvuYrVbS2+6xPQJku0Smn3G\nUONHXrnsaMSIffjmb/5jfvu3fuOlrEv0EisjIhb4fWJi5MfAbwO/oqrfHizznqp+0t7/BvC/qupX\nd6xLf/+T4+fe2JtsEDdEBwbPd+6gwMzCrB3K/pvfP+E3/uUxf/xoxc9/OOfPfnnOTz+csJaMkoyT\nWnm8ani8qpk6w8NZxsOpA+2zumWjrJtYLlPYWFpjB4mSDZa3Y98NpcK2Y51DV7Yryk7iC23Dsg7X\noXE2c0qcbGPbNd4ecj9ixIvg6+8fobo9iedquJQJqqoXkb8M/EPiufBrqvptEfnV+LL+LeDPichf\nAmpgBfz5/eu7fKNu4/mxUWw8eN637XSNV4Q1E9ZY1uhZSb2uWDWBiow6P8BPZlHxJSgTDRxObGfU\nVOG4DO3w9mjwfPCoxrpBawxOdUNNep9rPNzSblv3jRigd3W7wmq0y4bHkpj+85L7vLOwevhXtS/d\nuYW/94g3B1eKCarq/wH87NZz/8Pg/t8A/sbL2qh9Wn9di+wNZoTbCG3Gd90ECjxeaoSSwiiHE8uD\nxiHOcFoHHq89uY0tclNnEIHMGFZ1YFl7VquG+xNHMXPMM8O6CW1McHOSXMfmduw7GCRxtmKaQ1Z4\njs1ufS9R7QqrDXRiCzGmeL7dbliSM5bRjLhJuPaOkb2u0hWyxreSMaggxFidMQ7sBIzlrQfwNQfz\ntXA0cTwpPXK85p2Z452pY2L7Uud1E1jVns+WNVaUw1xwGByxVjCxr9AWPvfFzM+ywzZZ4a53Ktrp\nDRrp56SkT2qTxoQ2oRLos8m7WaG2ccJ+a2/lbzziVuP65w7vcNCGJ0LCruD5bWGA22wsswZrlDyf\nIfkBITM8LDxfn3veXTd8tqj4+KTkJycrfva+41AzionDiSM3DiNRnGFZe9ZNdINRExMwIjgTDaa2\nJTOYgRDCFXGVOSu7kibd+wd/DXSyCqF7546MtI6McMTrx41hgikmFR/LOfcJBifqLTKMSZnFIvig\nnNVxqLoT4d2DjKOJ4fHZmo8eL1hVNQ/MjK/OBVPYrm1tQ+BABBULxuFsoMgE33YBq0Zj6UTOtbY9\n41b35TR7oophkFTpahSHMT7RLhPcdaWwrzxnVKIZ8fpw/UxwV5aSrbjTlkG86L030SBuFgf397wq\ndR1YqHIoNYemxtQlYbHk5HTJSQ2LGhqTRan+Ng3hDBTWcJBbMhMn1cVh7bGMpnBC1QorKGDUoCrP\n7VpelRWK9omRZAiH3zixwtDGDXexwuH60jrTe0fXeMR14DW4w/tf2C7PiA908Lyc4wk3xfBdBU07\nJ6T0gZwVc1li6zXOe3KjOOcIxrGWnLVkQNwPhbMcFkqjMHUGr3HMpzXRQFoxoIHaCyEoQRR9aW0a\nfbxumxXGMZ/aLWESo0v/STKE2scHh7cBO+zfNrLCEdeLGyOlpVsPLjKItzWQ3rmFqqgREItxGUcz\nxwcPlbl3ZLnjuFKKlWeWGeaZYeoMYeKwxsQuEq+UjeegiOwwd0LpQxRlDQFn7CWlMS/xO7V/o3ss\nHfsb2uB0P7JG2RR33bPOxAp3xYZHjHiZuDFGcIhLDeItYgvbSRJnBRGLc1PICowE3jXKn5wpT6r4\nvX50sua0qvnSQc7sXkZuhOBAxHC8bjirGk5Lz/tkHGRCYVwc+NQysL54eds1f3bsE2EYGrCh0Ups\nL2U9usuVpN9PB8ucL9QeS2hGXDeu3QgOj+ursJVtgzisAE4P03pvOlswIuRGUAOSOaroy3JolT8x\nU56sav748ZLvPVkxc1DolPcnEyaZozAGYwwnJawbz5N1zb0cmkaQDIwGrEinMqOp/KT79Fe/c1IJ\nDUS2nozZ0JAa2SyuHl7k+vUM17i5Prj5v/OI24VrN4JDGajh0X9V9y2NvO1PcO3t4g2XbEozfxGo\nA9RVwArMnfDWRLAB/qCq+dHnCwonfDC1NG8VkUG2EVFDVKgpfaAJEDCIWKwVchsIxDpBVTbKZV4G\ndrHCXbqCXvukiEnPSs/weibYMkvRnTL9I0ZcB67dCNo0EANQ6TOF20KfCcNTTLfuDLPKg1j8jcI5\nLT4AYlyv8koIHitrZqbELEv8WclqXRKKgkot3hSozRDi3GBnhMxG+X0Rg0eo24lJmZVWMVoJIRoo\n2Tc38xUjMdEotyWYzrptusYIGI3XrzBIQw9dY2AzQaZju92Il4drN4JukOnYiDFJ32a1UUd2CUXU\nDeqQ+AXcPHMYkb5Oo5HNee+ppUFNiQ0VU6fcnzpMnqFGOK0DzsVZxpmB3BrmmeV+4citoQ7KovYE\njfFGgLVXag2gYMQwHBj/stAlOwaDnc6xwh0Xs9h3HO8nxqpIlz0Oe9MlI0a8Gly7Eczs4KRJV/1B\nqURqy+rr084rnWwjLSs7bnA5Y9iloHKV9z0vIsORtpXOYG0GdkImjncfCj/roBZHnjs+X3kCcFRY\nisIyywwPpxnORMWY2sNpFSisoXCCE6iD4KMNxBolQ15LPeUwYaLSJk4kqVYPYxq62WZH/3sPvYMh\nKwRegWkfcRfxGphgbwRJTGD4l1hCMuwwCAOGeJExfFG86lNqaOytESaZRTBkeQ7ZPTKULxSB7MDz\ndO05qxq++2jB2zPHTx0VHGYFuYGD3GCMsKw9i8rTqPJw4pg6h7OGlQSk3afopjF5HS6kDrIfifWp\n9ko0Q3uYXg8DN2DXz6303UejazziRfDaEiMJKZfRB9p7w5cO9FR6EXTTZd42hpe5xp1h2MI2S9rF\nPl7qSdaxochm6iCcNIIVw6ywfGWSMTkt+b2zNd/99IzTecaRgw/vOZwYpk6wrRE8rRoWtWdihYdT\nS2biazL4Lp2xGGTWX9pXke4rXewap4/WgVo1aQfLoCRGezmu4ftGjHhFuP4SmS1r0mUbB1d9RfoM\np7DBBodzLXYJeF7mGu/Ctju8b5nh6y/CPjZ19JQqKGuvWAKH1nPPeo6rkuXZgo8eLcFPOXsYUHFk\n1mAAF6KxqH1gWXmqxhNCiAkUUZyRzm0M7RdMwqivE6FlhfGC1BZXDwypSF/fCPQMUjcNYjSa/TFz\nWwvoR7x+XLsR9AMqtn2sdgYrJZAlHtxKW1+mm7eUYbwqW1C92gkyXORlGT7YzBSn7xqAqgks64AJ\nNblZcc+skWWJVA0aPD4oVYCVV5DQCjIoZpBGim6koEQh1sKBGcwA0cGSr9rxj99td2E17XdOjNC0\nr/dMsJft71xnNmPDI0a8TFy7EaxD76YOGVFn+NgqKZGBy6W7b8nlu0i3rl/veQOwSzDgIrwMVphO\n6KBtSUtLd7XdiNwZjmYZX7yfcX+egxGerj1NJsyckhslE2XqDFUGxkjLKGOtYExAKXU79c5Im5R6\nhTTpWZIvKamRGL3Q/p4b26c717lr/4+scMTz4tqNYNUMjGAa9CP0Q77b+7ticUIvETW8JZc4XBBI\n38aLZEtf9ghJZwyFBWsznBNwcQD7BwaqGSAGsZaPziremhjenhiKIiZW7k8F55TcWUofmbEVKKxg\naUd/BsWIYE2//c/zvZ8XG0pA3W80jEFsGT4GF8Z0kWOgRLNH9HXEiOfBtRvB0uvGHNuhC5wMnGlP\ngI3hPoPlNd3fwQp36tZp/zkb2nz0z2/Hm3bZh2Ey5rnY33b8sn3OGcFkFmcsNitockOm8MVp4PAw\n8GhV89lZxe8fr/nwMCeXnIeTgmmmPMAycdHQretA7QMHmeXAWQzKqgldTDCoXKvqzmV5mLTPk1DC\nxvGQWGE6JtANCf+LGPswvjy22424DNfPBH0MVMWDXjaN3sAIbt/i86nebfA+Nm9wQXmF9CfIRo3Z\nDmZ3Vdf4WbG9vrQPkqe68sp65cmMMHGGo3uWqm74o2XJd358gvg5b08dGEvhQIzibODxsuHpuo7S\n93O4VxgyG0txUqihE1XgfJb+ZWNfiOGqrDAljdIzQ1ZIGyvceO+IEc+J648J+vYwbg/yzpCJxLkV\n5wzfrtvmGZwMaKIWEvrMcR+cp88J6GZ6YOOEHdy28aIxpr7LYrBN9OxUNQ5kKr1S0OBsTWZr9LTk\n7GTNJ8dr3p3lnJWeKsTvPLGCwfAIZVE1NEG5PzEIWRRUEPpZxEmKrNsVN4cabbPCtJ29p9CzwqES\nzb6h78M137TvOuJm4bVJaSVD4LsjWDGirUGLBtHuNIADIyn0yRWSUdNztG6jvKKLHw5l4XvsMnLb\nUk/wYgZxaHS7dUGnFbisPSHU1HaFmBW2brDBt2UvUAVlWQdyCRQSMN6jTU3tPVWQKKyg2oquGgqn\nnaiBDzHOsB0WuCnYZoW7tlG6/8Agbe3oblPYXwR1dI1H7MT1F0ub/v4w2wd9A31AEVH8wCAaAWt6\nA2g3DOHQoMi52FLqWOhOk5aF9mzw/BmxbRhfhWuc0GWKu60RxER1GIxjWhjeOYKfCvDgoABjeLKO\n9YTGeox6HHEmsUgc1dkEOsWazArio5H1bbW5vCxpmStgU33mvGu8r7A6Se13SjTd+lIxuHZdRun3\nHZ3jEc+KazeChTV94kL7uJ3q+Tqw1F4VNGY3g0YjaiUWUW8yxT6WZDiPrieZPlEydI2Hbnl3wtKP\nmLw0yv+c2HaNM2uYi1AYwbkMzebMM+UrecAexKxoHZTvP13x3lRwU+HAGSaZcFRYyhA7T9ZeCW3k\nrLAGIRCaWDKDiaouzlxvouR5kArmY6Jk1+WKls23r+h5o9qta/DKqEQzIuH6BRTMsIhW+iv4Rr9w\n2xGSsprtfQPY1p1LvcaJISq6ESscdiHEz9oM0nd9tQOjuT2hLWWgN4JUA8bxMlzjIaTdP651Y70z\nnBmDmSrv58rDe8pPzko+Oin5/qKCBwWHRc69LGMOvGNh7eMFYdUoTYilMoUVwFB7bWeQpIvQi2/z\nM3/HHaxwV7JkeD8JTpjBK9362IrlDlj7DbfvI24IbriKjPQzOdr3B1VCW/xrTTwxNLRzdgcucndq\nSd+VANEl7HqTZZPgxS4H3ThRzw0HGhjC58WuUhlgo1YSYtzvrApMLRxkwtuF8PjU8+R0xXc/XXJg\nhS8cFGBzChQxStYETivP2drjrPBw4jjIHEG1LS9pLzb6Eqz2NWDbMxjGfzcTZDt+kQt+p055eyys\nvvO4flHVjSI9HajIyEBFZtt9lY4ZhjA0iP2qVHtWaKErtu1ZR2t8u4/vWUZ/ckVe0hsi6RInXVxx\n97n2TKxwu1B7m30aicrTyzqwrD0HUpPbEmNLmpOas9Oaz88anixrziof+45FmGcxdnq8juUymRUO\nMoM1kGlkmFaE1NLmW8OYZpO8DgyZ+kXQDYu2O5bYs8L2eBmseWSFI/bhNQgobDw6173QxQtJBqpn\nbkYEFekHf2vL7FoWqNBaNOkyzelmTevahj6QnlgmxOlv2hpQht0rg+0OGllnev/L3x/9vghBqX2g\nbDyF1IRQIWGNDYHcKEVmoHV7z0rP1HhmJmBqj68bVrXHB0sdoilIFwhn+89sw4ODBNH1Yvjb9xGH\n3a4x9Gw8xQdhwAx3sMJhGc0+bLjfY/b4TuLGTJvbrNVLsyfoeoOjgEIa5LPN6jYZJMS5u5Z0ksQF\njcQge+o+iOtr51uE9IFtKcmAUZp2m+JH64a+4TbSdjzvyTR0va0RcmvJJMO4ABYO5soXvbDM4eFB\njkd4vK55YBsK5xEfkKBIF1elnUUiGDHktt/uMDQ+t+ik77P9+93hxMiTBxA4b1Q317kZ+h1xd3D9\nKjJ7Lsvb/cLdfQYMoKsJ1C2jF9F1EbRnv5hNgQYhsqGYbox/kmH1bawxmJZbtFbBDgzjdpdFx0i3\nv8tznkZDNiwiFNbgjGFmHS6f4nM4zAJfnXrmh54mKKvG89GJhynMppAZQ+HgoAAxFpE4kzi1HxZO\nqH1ss2taKnid5TK7sN1dAldxjwdKQrIZKzwnwkB7IRsTJiN24DVIae1+XgbXdRkss90/2ukM0guu\nDmNyG/V87VmSSmhS7EtEMWYQb9RYTOyDYoKiSWnApoJqoVdBFqQtaEyfcz7R8eKjIo0IxgpZ+76V\nCk0N4hzv3VPensP3n674/vGak3XNRHLenuZMMseBwjsmsj8nhmUdsG08cNIavCYEvCq23X4ZWPib\nUjZzkWu88VxyG9i9j/vnhP4t+9c3JkzuFl6znmB/oMEwS8vG4d8tPyyB6ZbTTmx1ePKmhEtKenTv\nM/Fz44SzmCIQUdQHGqIxrEPP8MQOP6uNoUnsVIDzCjYvC+mzko7iqg6cKhzYwKFTpjbww3rNk+Ml\nHy8a3p44lvdzjvKCqVEeuth+V/mYYS6sMM8NU2cIqpRIXy4jN8jwbbHCqzDCFCfs55TsS5iwkXwb\nWeEIeM0xweGVd3g0yoADDN2cDYNIcoNkcJAPag3p1UaiFUvBcunYnbZ+Ukq6CILXmJCoPfhgovIK\nhsy0ai8iYLRVrRGkS2XvP6mGjPEiZrEd7E9GMBmys8qjZs3UlsxlTXPqWSw9x0vheO05qTwPGsUa\nOCoMizqKtX6+qplnltxmZBODbeKGRB1D2RSpvaGs5zJW2D8rm6U0e8tono0VPm+Md8TNx41JjAxx\nzp1MMZ89J2rH0NrHvcJIMrDaxgGlG0aesqXpJFEEVaHy0ATFB0U18gpJBlRSBjHpGsZIpcrFrnG/\nnXJumX0n1dCIR03AEF1YPEqNSE1uYJ4Z5kWM/a2awFnVMDeBiQ3Y2lOVnpOVB4V6YqMqtRGciTWb\npqW5vQnp/3/d2Ndutwvda8kFkHOvxHWyyQpTQf7ICO8urr9EZs8JdmGrU/uSQTZOjDbEt+nupvcN\nTpphg71oKrLdFGiIPbeCSCA0Ps4EDoL6mEEOIXJII9qPjDQxPihtG5pvre5FrvEWF7kQnRGXOHB9\nmkFuC4wT1OQcHsGXjJCVwr1JxrJRHi1KyBqmmYdaCU3MDvsQOpZsJAq5ZqZX3/GDgvObYQKfHZ1r\nDBttdhexwudtt4ORFb4peK3T5jayoZe6O62wQnJDZbO0Q9p1dMKbbE6mC+1BLqF1n00vxGBFBokV\nofbR7fVBKZU4xChLbFDIjHQ9zNJSUAmJXeiFrvEQKTbVbf/WvknvNyJMnaVwllmWY7J7hEw4yhp+\neuJ5sPYsa89J2VBWHjdT7otiEHIrTJ0ha5M9tY8hg6Q+7dukUIwJ9Pv2pmG78L27wF0xYTJ4ol8n\ng1rCnUuMuAt4zcXS7XNbj4cH+S4k47H/8N7/Ht8aUgI402eNM5uyx6ab+VH6QNUmF+K2C9YYtC04\nNmYgiJqi8u1G7dK52ynouvG9B9u6kRDqRVe9wlmtrBrFiOG9A8v9SeCPHi/56LjCe8+RzajnOZm1\nHCC81TJJEeGsTsXTwiSDsgHvY6bY6C7piZuLy+r+YjxvByNkyAqHoY7NpMlFnxoTN23w4AZeNEZc\nHa81JrjNfIbPXyUWNDSGqRYwoc8UMuhFjs817UAjkejGJu1Ca2OnCUR2KJWnDp7KK9F8xvjcPLeA\nwWb2nIqNpv9C+syrcYvtpYaGUuhlxCqvLJsoo38/8zxwgSI06HrN8VnJ2gunhxPqbM504riXKSYP\n8TsDZ1Ugs1GlpjCGRjwaoPFgRLEqG2z9JiIdHxfNOU7YbrdLIZSkYCSxorqNE16uWK2Dv8mgjrjd\nuHYjuC8p8CIGsTM8G+hP5kAvidUPddcNhmAH2d8499jEUZc+sDaCD7CsfWxHU0VwWBEyG13OlGgB\naRPR0RcfZqv77R2wvCukIobqNrUPnNWeRdWQ1yWHdk0WKsI6UJfKSh1nledpGXAuCqs+mFjOKs9x\n6Tleew4LSz6xTJyh8oIxLcXU/oKxHWu9zRiyQoRzslxDVthdPK8Y0hhx+/F6S2T2HGGXGcRdXQXb\nrHDfepMhFCIjTO/NRXBEQ1ZgcCb1LCshwKrxrJvAugmd1JUzBpVWzaZ1q42JJ1Dcdu0M9FVOpssy\nxWmQVGiVcFQDqgFpBVUPJpagjiDCaeWZOOFeBlMHa/WsK8/nqwYjWTuDRDolnm4fDty9dquusOWv\nB5s1hVdjhd2MldbKb387IVUaSLffL0uYjIXVtxvXbgTTQQgXnPRXvPxuRnU237/NMjcC4K1R8qE1\nVDEngBVaAwcZguQGcBgxyDrW68WaPU/Q+PiwsFBYJOvLT9LAeG0Nl39G13jnPiElgoTcGsghdxOM\ncxjx3A/woROe1kLmLI9WDeprpAgc5AGtoKmgrKH2jhC0i2c6I7iWzSqRFKa+6ZuMdFF7FnTsTjcl\n1jbWO7i/ESt87i0dcZPxWphgF1fZOqquWnqwzRSFiw/QJLmv7dHv2+dT610IrWvsWyZgYwZY2kSJ\nNYIPgXUT53zUQTktG9aNRzWLElXGkLezfq1J2yNIa3U1XJwsSYxi1z4Yqm5bI8wyxyyDST6BzBCs\n8CBrkJnn6brh6brh0bKirjxzhWDjfA2DjQa93W++FY1wVsjVbPwuqbPihtvBDvs6TfbHk1Nh/XDi\n4b5Um3QPx3a7Nw+voU4wYm9sj96wwdUN4lWOt+ROG3px0U7QQQVtDYFIP6Vt4oTMaitLldEorOpY\nkrKoFWdi5tWIMM8MTizOxPfGFr30neLnXDYdbd8FIj20yeUWoQnKcekRiYrU784ypk44Xdc8XlSs\nTeCdPMObHJcZZkF4gDDJDIqwTllhgcxGo9iE1mWUm58g2c8CZXBh3J/gSHHCsd3ubuO1lchsn1/b\nSY/nMYhX/fxhQbVIkpuKMzm0TahYicXEmUBhILSucW4Nny0ryiaw8J6T0hOoKH3g4dRhRJhlpk+U\niPQnTUs994lIXLDV3T4YDplaNYFVHfAh8FYWuJd5fBUzxSeLkjp3rMhpigMyYzhyMVOcMtlnVYiF\n0yJYK6w0dOVBRg2qsnMw1k3HLlZ4ISMkHWODyYVju92dwY1pm7soC7xtEF/kgBpOo+uv7kqjoE3b\nJaLxsyZGKcSTEyjUcM8ZMuMoveesMixqoWx8zCDXAVHITTSAhTU4Fzc0G2QeVDZFJIbYlTXuH/d/\nU6F42QROy4a6aZjVa4IrMXWNVEBQmuBYNsrTCma5kDnhgTWUPhrPRR2YZYZ5HlWnax+3MQRunLDC\n86N3Ty+MyZ5j6JewQh3EmEfcatwYIzjEZWUxV+m9vcpndNqCBkKrseeD9ge7i/HB3BoyMVixzMRQ\nekfpA02ARdVwVsXSmdwaJpnFtn6kaw2iFRDbswxaDcPLXOOdLEJ6RezkuvdzVRQrMCssDw8MajOC\nGJ6UHhHhXm44LAyfrwJPGs/jVcPbs4xZJmTGnM8Ua5vI2pFFvS3YdSydj8sm5qZ9zellrLC9inYj\nXfcgxVVfxjE74tXgRhrBIS4yiM9zLHWGKJVItAEhL2kQUSxI9q14gjPRBXahwfgS62uOxMLUMrGW\nj07WnJWes6rhkalQVcom462Zw5qMiY1zgDMAk2ZfSCzLC/rMTEu1LenTtuvDGbAZhRMky3EhcF9i\npnjZikF8dLKmKsFM4Giq+DWsSzirhUNvUY0JF9u2A0YvXnckR2772XsxK9Rzd3Y/kY6hoUzb7hj3\nwDU+x+1H3BTceCM4xK46wWeNGXZXfeldSydtllBMp7q8akJXRwhQ4ClCTa4VU5NjJpYsdyzrONRo\n3QirxvPpQqlCQAQmziJ5LKbObB8bVG2zxgKN9oXcGwZnoEwzjC2lU0shMlQb2WaeTfAuGtyHeSCf\nBZ6saz5dVHxyVqF14AggU7SxhMbQeItPWevWqObWkmTHWu1YxLw57KVPmOzGMGFycbsdRJG1eOz0\nLPPWxw/uHG6VEUx4loww7GeQ6cC2EAeSiyCEthg5dossqsi8DowixpAbR2HirI5CleXUUDY5RoST\ndcPjRcWq8mRimDqHIBzkpusqcSbFJbXbtsQqrpQ1HuyDqGYTjVfZKOvGR2EEJ7w7t6CeHx97Pjle\nY2bwxSIjiMM6Q54ZZhon0XlV6hBQje60k0FRtlzcx32bsF1cfZHROtdu175XB9RdupXFvxex+lGJ\n5ubitRRL78OrPDAGx+rOV4fzipsQXcOgcezlqgloYckmE+ZTi9QrTLUgqyvukxPmOZl1/Eu/5POz\nknXtKWw8ccom5515RmYznEgUde2k/eNfnxRonjFh0rfqwbr2LOqAUc/bWeAw90yqGsqSVVmxKgoq\nOyFM5+QO7lvFZMrURUHBsmnVZVqXvQ4hlg+pRuEI2VSXuf0Jk8sxZIXQew5bIdq+BnWQMLkDu+eN\nwWuV0hpiW//vZRrEnSUTg0zz0O1xIhStW1k2ytoHyjr0qjHGMKkbiqYi8yWHmaPIHEVuWVY1T1YZ\nJ2XDaeVpTtZUIZahzDLTSlrFQmxHb0iEGB/kiifQcLshXlhKHwVVja859CU0FbZscCEq5IixlGo4\n8RY1wjRXMtePE6iDkpnUBy1Rg7DdWVYk5Y/eGDyrCMOudrtOkLc9mPqII+1vuX9926xwZISvDzfG\nHX4dGUjt/oswraHLBIxYwFOH2FFS1oFTie1yh4BgyayjECXXNSKG9yfK4ijnk4XhdF3z46cVde2Z\nGMO9whEKOMCQt0OPNM0ERUljplIR91WYRKc83bb+RcXrFIJXMms4nFre85ZpkVGr8OnCc5AbDnLD\nUWFjh0kVe6IPC8vEWTIjVD5mRVKGODH4zvi+IVTnqq133cVJhxeDTddYBssNy2j2rX5MldwM3Bgj\n+LqQWOFwCLkTwVlQNTRW8ZniAxyXDU/WDe9NLW52j2kOplpg1idMmoZ3zQw5mjFxGd9rAp+UyLV4\nHgAADW1JREFUJZ9rxTSzWCu8N8/x84zMxdPIisTstA7YQapTvMTKpBM3ldo4Y5jnsee5yCySFWSN\n8tAIX8mgDMKy9vzh4wVfmBrcXDicCKGEZamcNbE7RrBtG2C8+dQAvb3TbvkJ3Mf4NhNu27HXbYSW\nrovQCdCaDRq3FX3uknjn1xgvfX2fefuOkRVeM+6cEdynRKPaljxIywiB3AqaRwO2qDxnVWBReSYu\no6gN1hqKOjCp1mS+5uF0wnyS4XLhdF3x4xNDo8rjVY1v29yMEeZZZFuZjUXKw5Ovc870fEInvtJL\ns3YqNUDuDLmLXS42M1TOIBozxZN54LNFxY9PS358WpJ5eMsIYgRfG6rasg62m7KXTmwrycGTjgnG\n8pnbf5Z2RmeDCV7ebge9CIdhtwjDBiuUzfjvrt+z/+QRrwN3zgjuQxcn1LZ2UGKSYCKxkDjF3URg\n3SiP157Kex4EwZo8zgARz8SfUWL4YKacvDXl8aphWXkef3ZGaDyFEQ4LxywzzMQwsVHP79wPESLr\n2Ndit51ESWowAqxrZVl5nIGpM7w3c6zKmqpq+PjpivtiWU0czGObnzUSVXPoXd+kti1hUNqzsbfe\nnJN226ZvxI45fzFKi6to2zHy4u12w3e87Jj4iItxI+T1t/E6M4/KIP5Fm4G10ZjUmSWEWEf4aFnz\nGCXMMibTB0xsg1RnmLPPmDfCF90B8vYBPzjx/PGjBR8vKibWMM8dk8zycJbx1tThivgTGIkSXkMD\n06r6x1KVre3cLhNKytNBlUUd5b4K8byTB47yQFZW+HLNqmwoQ0bjpjCdkhu4ZxTTQOEsoU2UqEaj\nCvH7Nm1do5ibL6zwLLhMhCHe1Z0mP8UIvSaGfLkIQ6fXcW5dw/iivGGXmZuNG8kELzOUL8NIns8Y\nD3iOgqedRtf+nVgDhaOwhifrhtWqZtUEDvKcGQVWPHl9Rr44oxDhvYMZR/MpzjYcLyt+aIRlE/h0\nUaFGKH1ohx3FWj0n0k57iyM8SQo0gZ6hwTnXuWMf7VnjA6zbnuKGmvu+QpqSrAlMDMwLi80yaptz\nJhPEBe6JMvFRTaZuWwfTQCmhH/mpkiTJbuYwphfFru90FVYoLSOMzUcXs8KxuPrm4UYawcvwKurV\n0tUX+it8JD9Ra9BamEoseq58YFUbqhBFCD5Z1KyN52EjPHQFE6s44znQU86c8MFcePJwxqoJnKxr\njlcV+BCZYeZiQsJF+S7TuuFt116bLBl81z2sJChxVkhb6G0lahsm1lY4y8O55QMs80nOOsCni4aJ\nEwoXZcCWdWBZx6jj1MWSHhGofB8zfZPrAy9ihd3/e/b/7vefZ4WJFl60G0clmuvF9dcJ7nn+olKC\ni7DvwHiWk3WTFfZtaZ3r4vtEQdYmNuoivlj5wMenJY9NIEwnzO+9x0RKpFzC4yfc8zkfFveQLxzy\nw+OSHz1Z8fHxmlnmmBeOIrPcnzgeCBiJeyd1lhDiCNE4CjOeOaksZvj1umRJ+1puDQcFTE1BlmWQ\nTZlm8I4TmjzOIU7iqx/MHR8cOGa5YekDizLQAE4cB7lBbN9TnHZ1GLh/byp2xQlhPyuU7ZtsPt5e\nT6pGuKxzZcSrx41hgtsHyy48i6F8GWwxGUIFELquklluMSaKk362qDgpK5aqHM1mPJhMsGGJO1vg\nnj5i5mZ8+eiAtw7nZGI4Xdb80CtnZcOnZxXOGZqQ4ySPQqzSzyyxomgr8DAcyOR3uO/ppAShaGcU\n51YwmbByBpMpbxXKfK58fLLmRydrPltUTLzlbTFIsDSVsK4NDZYmj+UbaS6zHVwpNop939DI1fOy\nwn0zj4fH47aBHdvtXi9ujBG8CnYZSt36u/N9z2AQu/ga/QGo7RtjU31kRqnMZVkHJmVUmj5ee35w\nUnGmFQ8q5YFxTDLDJPM8zEqOp4HP71k+X02wVni8qqhCQINSWEPmDIUVik4eK7q0fUF1CwNRCXtH\nwgTaOr/4eFkrZ1VDbg0zJxzNLI/OYLGq+ejxkvdsRjnNYKKglm56ffuJRmmn6pmOGW+Wy1y8P98E\nXMQKhyGCnaxwixHK1jLJVl6lc2XEq8GVjKCI/CLw3xCP+19T1f9qxzJ/DfglYAH8RVX93V3reh4R\nyu3zbJd7scs47jqUhkWyez9PdgfEfcqaSpKeMhwWlspnuFJYNYHvPVnxxHh+anLI0TsFlhqaFXz2\nRzyoJ3xlNsd/cMSjs5JHZxUfL1dMjOEgj5JdB7mNQ5RaKxbb9UBD/w0Vje7xIL6U4k1CXN6aKAl2\nUjaclp5DG3CF8qAI2HJNU8Ue50oLvM0hK3BBmARoNI4L8Np2skhM3Hi0E1Z4E2oFr4rLM8gDVr7x\ndxi32J1M2pCK21zV7m0ZbFMysiNeDJe2g4qIAf468G8D/yrwKyLy9a1lfgn4V1T1a8CvAn9z3/qS\navNFt3Pv2bqFHbftZeBi97rrGb7wFlvQUr4vaOqpjZ9hBe7llvcPct6fF2RGOF03fF4ZTvIHrB7+\nNNXsXeqqofn0h9xbP+Gnpsq/9oVDPjyaMnWGVe35f3/7NzvJqyfrhlUdqL0SWi2rJJRgDRvuqUlb\n125vGgSfYoONV9Z14LTyLMuSZr3ArE7I/ZqZDRxMM7Jigs/mVNkBWTHlaDrh/jQjd4baQ9XSvsxG\n1xh68dlXnST55m/+41f7Ac+B7WMk7fNhDzr0oZTt264CeDO8Scuwd3y2onzzt/7RyBRfMq7SE/8N\n4Luq+n1VrYG/Dfzy1jK/DPxPAKr6TeBIRN7btbK9B8Yz3HZh2whuuybPi3SgQ5y94VVpNBkCyKzh\nXu54MHUc5JbCGTzweO353nHFj04rjtceDcrcCe/PLV97kPPlw4y35xnzwvHH3/ptnqxrfrIoebKq\nOatjL2/VfoaqtPG5tp1taBTTUCcZGmuoW0HVZKyNGEQMKsK0yHjncMKHD+dMJxkLr3y2bKiCMs2i\nArURKL1SttXahW2Lp6GbQxJesRX87d/6jVe6/ufB7uNQGF4uz72HaMDS8bORUBFpb1sXX9m9tt8Z\n7JPkoVx0Xoy4HFdxhz8Afjh4/COiYbxomY/a5z7ZXtmFir4t0tjCXbiM3V203mfFRjCbfnaxamtg\nJMbfcisIhofTjDrEpMdJ2fB4VfOe9fx0do/7730ZZy1UJ4TPv8uDOuNL84zGzPn/rGFZNjQ+UBjD\nLLNx0p0Kag2uLYROBj31qvZxJu0YsuqmdH9hDYeF48BZsrxAC+XAwhcycNPo2j5a1pyUDe/Pc94/\nyCis4Wkb5xTgILdMnYlFwV28sGc7d9Ut29lpMrwNjF3nvraPX1SEYVStfnm4kXqC5xr2B1B2G8ln\nMY773rczhtj+NUgXz1RiO1vQWFSsREb4YBrjep8vK358UvLJaYlOHW/de5fqrQJZPILjT+D4M44m\n7/LV6bvMDg75vzOLAGfrhuO8ZpZbRKLxkTxuQZpVQssaLFGYVTTurmT0hszAGsMsN8yAiRXIhDNn\nsEZ5O1MOpoEfHq/50VnJqmrIfMXbNna0NKWwKMFayzSL80dsO595uL/63tfLfoU3D7vZ16YKeP9/\n/0w/4U461/dF2+1GPD9kn5Bnt4DIvw78F6r6i+3j/xTQYXJERP4m8H+p6t9pH38H+AVV/WRrXeOv\nOGLEiFcC1QtcyAtwFSb4O8CfEJGvAD8G/j3gV7aW+XXgPwT+Tms0n24bwBfZyBEjRox4VbjUCKqq\nF5G/DPxD+hKZb4vIr8aX9W+p6v8uIv+OiPwhsUTmP3i1mz1ixIgRLweXusMjRowY8SbjlYyNEJFf\nFJHviMgfiMhf2bPMXxOR74rI74rIz72K7bhpuGy/iMgviMhTEfmn7e2vvo7tvE6IyK+JyCci8q0L\nlrmLx8qF++WOHitfEpH/U0T+hYj8cxH5j/Ys92zHi7aqIy/rRjSsfwh8hThz/HeBr28t80vA32/v\n/yng/3nZ23HTblfcL78A/Prr3tZr3i8/D/wc8K09r9+5Y+WK++UuHivvAz/X3j8Afv9l2JZXwQRf\nanH1G4Sr7Be4Y3UmqvobwJMLFrmLx8pV9gvcvWPlJ9q246rqGfBtYj3yEM98vLwKI7iruHp7Q/cV\nV7/JuMp+AfjTLY3/+yLyJ69n02407uKxclXc2WNFRL5KZMrf3HrpmY+XW6UicwfwT4APVXXZ9mP/\nXeBnXvM2jbiZuLPHiogcAP8b8B+3jPCF8CqY4EfAh4PHX2qf217my5cs86bh0v2iqmequmzv/wMg\nE5GH17eJNxJ38Vi5FHf1WBERRzSA/7Oq/r0dizzz8fIqjGBXXC0iObG4+te3lvl14C9A15Gys7j6\nDcOl+2UYuxCRbxBLmB5f72a+Flykc3EXj5WEvfvlDh8r/yPwe6r63+55/ZmPl5fuDutYXL0TV9kv\nwJ8Tkb8E1MAK+POvb4uvByLyvwD/JvCWiPwA+M+BnDt8rMDl+4W7eaz8GeDfB/65iPwzYlf1f0as\nuHju42Uslh4xYsSdxisplh4xYsSI24LRCI4YMeJOYzSCI0aMuNMYjeCIESPuNEYjOGLEiDuN0QiO\nGDHiTmM0giNGjLjTGI3giBEj7jT+f20vCJZMsyiZAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "extent = [0,2, 0, 2]\n", "fig, ax = plt.subplots(1,1, figsize = (8, 5))\n", "ax.imshow(data[0].T, cmap = 'RdBu', extent=extent) " ] }, { "cell_type": "code", "execution_count": 29, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtMAAAEACAYAAABxtvotAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcXFWZ//HPk6VJOmuHkIUQEiAk8kKgA4YlIBQENAyb\noEhYBsIoKMvAIPIjIoqKMoAjDjKOMAiIAwiKLAEMBsEiELYY0iwhGySdjSxk6exLJ31+f5yu6U6n\nuruq7umqW9Xf9+tVr65bdW7d00/drnr63Oeea845REREREQkex0K3QERERERkWKlZFpEREREJEdK\npkVEREREcqRkWkREREQkR0qmRURERERypGRaRERERCRHQZJpMxtrZrPNbK6Z3Zjm+TPN7D0zm2Fm\n/zCzk0JsV0RERESkkCzqPNNm1gGYC4wBPgWmAeOcc7MbtSl3zm2uv38I8LRzblikDYuIiIiIFFiI\nkekjgXnOuYXOuVrgceCsxg1SiXS97sCqANsVERERESmoEMn0IGBxo+Ul9Y/twsy+YmazgL8A1wTY\nroiIiIhIQeXtBETn3DPOuYOAM4D/zdd2RURERETaSqcAr7EU2LfR8j71j6XlnHvdzDqZ2Z7OudVN\nnzezaEXcIiIiIiIZcM5Z1NcIMTI9DRhmZkPMrAwYB0xs3MDMDmh0/3CAdIl0inNOtwC3W265peB9\nKKWb4ql4xvWmWCqecb4pnoplXG+hRB6Zds7tNLOrgcn45PwB59wsM/uWf9r9D/BVM7sY2A5sAs6L\nul1pXXV1daG7UFIUz7AUz3AUy7AUz7AUz3AUy3gKUeaBc+5FYESTx+5rdP9O4M4Q2xIRERERiQtd\nAbGEjR8/vtBdKCmKZ1iKZziKZViKZ1iKZziKZTxFvmhLaGbm4tYnERERESktZoaLyQmIElPJZLLQ\nXSgpimdYimc4imVYimdYimc4imU8KZkWEREREcmRyjyKwM6dMHo0DBoETz1V6N6IiIiIFD+VebQj\nr7wCmzfDW2/B3LmF7o2IiIiIpCiZLgKvvALnnAOnngp/+1vm66m2KizFMyzFMxzFMizFMyzFMxzF\nMp6UTBeBKVPg+ONh5Eh4//1C90ZEREREUlQzHXM7dkDPnrByJbz7LkyYAG+8UeheiYiIiBQ31Uy3\nE9XV0L8/dO8OhxwCH34I+l9DREREJB6UTMfc7Nnwuc/5+xUVUF4Oy5Zltq5qq8JSPMNSPMNRLMNS\nPMNSPMNRLONJyXTMNU6mAQYPhsWLC9cfEREREWmgmumYu+wyOPxwuOIKv3zOOXD++XDuuYXtl4iI\niEgxU810O7FwIey3X8PyvvtqZFpEREQkLpRMx9yiRb60I2XwYP9YJlRbFZbiGZbiGY5iGZbiGZbi\nGY5iGU9KpmPMOT8K3TSZ1si0iIiISDyoZjrG1qzxJR7r1jU8NnUq3HCD5poWERERiUI10+3A4sW+\nRrqxvfbyF3ARERERkcJTMh1jTeulAfr1g88+y2x91VaFpXiGpXiGo1iGpXiGpXiGo1jGk5LpGGta\nLw3Qqxds2QLbthWmTyIiIiLSQDXTMfa970G3bnDzzbs+PmgQvP027LNPYfolIiIiUuxUM90OLFmS\nPmHu10910yIiIiJxoGQ6xlauhAEDdn98r70yq5tWbVVYimdYimc4imVYimdYimc4imU8BUmmzWys\nmc02s7lmdmOa5y8ws/fqb6+b2SEhtlvqVq70iXNTGpkWERERiYfINdNm1gGYC4wBPgWmAeOcc7Mb\ntTkamOWcW2dmY4EfOeeObub1VDNdb5994M03dz8J8dprYf/9/U8RERERyV6caqaPBOY55xY652qB\nx4GzGjdwzr3lnEtdeuQtYFCA7ZY055ofme7dG9auzX+fRERERGRXIZLpQUDjC1wvoeVk+ZvApADb\nLWnr1kGXLv7WVEUF1NS0/hqqrQpL8QxL8QxHsQxL8QxL8QxHsYynTvncmJmdCFwKHNdSu/HjxzN0\n6FAAevfuTWVlJYlEAmjYkUp9ee+9E/Trl/755cth7dp49VfLWtZy4ZZT4tKfYl9OiUt/in05JS79\nKeblqqqqWPWn2JarqqqoqR+NrK6uJpQQNdNH42ugx9YvTwCcc+6OJu0OBf4MjHXOfdLC66lmGpg6\nFW64Ad54Y/fnnnkGHnoInn02//0SERERKQVxqpmeBgwzsyFmVgaMAyY2bmBm++IT6X9uKZGWBitX\n+lk70sm0zENERERE2lbkZNo5txO4GpgMzAQed87NMrNvmdnl9c1+APQB/tvMZpjZO1G3W+qaO/kQ\nMj8BsekhNolG8QxL8QxHsQxL8QxL8QxHsYynIDXTzrkXgRFNHruv0f3LgMtCbKu9+OwzjUyLiIiI\nxF3kmunQVDPtXXMNHHBA+rmk16+HQYNgw4b890tERESkFMSpZlraQEs10z16wJYtsGNHfvskIiIi\nIrtSMh1TLdVMm0GvXq2XeuRSW7VzJ1x8MVx+ub9wjDRQrVpYimc4imVYimdYimc4imU8KZmOqZZq\npqHt6qYnT4YZM/zUfH/9a/jXFxERESklqpmOqf794b33YMCA9M9/4Qtw773+Z0iXXgqHHQZlZfDq\nq/DEE2FfX0RERCQOVDNdwurqYPVq2HPP5ttkOj1etiZPhjPOgHPPhRdfVF22iIiISEuUTMfQmjW+\nJrpz5+bbZFLmkW1t1cqVsHkz7L+/r9fed19f8iGeatXCUjzDUSzDUjzDUjzDUSzjScl0DLU0k0dK\nW4xMV1VBZaU/wRHghBN8qYeIiIiIpKdkOoYySaZ79mx9nulEIpHVdlPJdMqRR8L06Vm9REnLNp7S\nMsUzHMUyLMUzLMUzHMUynpRMx1BL0+Kl9OgR/qIt8+bBiEbXsTzsMH8SpIiIiIikp2Q6hlqbFg8y\nS6azra1asAD2269h+aCD/GNbtmT1MiVLtWphKZ7hKJZhKZ5hKZ7hKJbxpGQ6hjIp8+jRw19WPKQF\nC/zJhyllZTB8OMycGXY7IiIiIqVC80zH0BVXwOc/D1dd1XybP/wBnn0WHn88zDZ37IBu3XyCvsce\nDY9fcgkcfzx84xthtiMiIiISB5pnuoSFKvPIxpIlfpuNE2mAQw6BDz4Itx0RERGRUqJkOoYyLfMI\nWTO9ZAkMHrz74yNGwNy5Gb9MSVOtWliKZziKZViKZ1iKZziKZTwpmY6hUMl0NpYtg7333v3xESNg\nzpxw2xEREREpJaqZjqE+ffxocN++zbeZNw9OPRU+/jjMNu++27/WPffs+nhtrU/c163bvQRERERE\npFipZrpE1db6Eec+fVpu1xYj0wMH7v54584wZAh88km4bYmIiIiUCiXTMbNqFey5J3Ro5Z0JXTPd\nXDINfno8lXqoVi00xTMcxTIsxTMsxTMcxTKelEzHTCb10gDl5bBtm5/SLoSWkmmdhCgiIiKSnmqm\nY+all+D22+Hll1tv26sXVFdDRUX07R5yCDzyiL+EeFP33w9vvgkPPhh9OyIiIiJxoJrpEpXJHNMp\nIeumW9ru8OEamRYRERFJR8l0zGRa5gGtJ9OZ1lY5B2vW+FrtdDQ9nqdatbAUz3AUy7AUz7AUz3AU\ny3gKkkyb2Vgzm21mc83sxjTPjzCzN8xsq5l9J8Q2S9XKlbDXXpm1DTUynbqEeFlZ+uf794etW2Ht\n2ujbEhERESklkWumzawDMBcYA3wKTAPGOedmN2rTFxgCfAVY65y7q4XXa9c105ddBqNGweWXt972\n5JPhxhvhlFOibXP+fBgzBhYsaL7NEUfAb34DRx4ZbVsiIiIicRCnmukjgXnOuYXOuVrgceCsxg2c\nc6ucc9OBQHNPlK6QZR6ZWr26+RKPFNVNi4iIiOwuRDI9CFjcaHlJ/WOSg0LUTKfmtm7JgQf6qy62\nZ6pVC0vxDEexDEvxDEvxDEexjKdOhe5AOuPHj2fo0KEA9O7dm8rKShKJBNCwI5Xq8sKFSebPh9Gj\nW2/fvTu8+26SIUOibf/116Fv35bbDx+e4IUXCh8fLWtZy7svp8SlP8W+nBKX/hT7ckpc+lPMy1VV\nVbHqT7EtV1VVUVNTA0B1dTWhhKiZPhr4kXNubP3yBMA55+5I0/YWYINqppvXsycsWgS9e7fe9oYb\noG9fXzcdxd13+8uF/+pXzbd55x244gqYPj3atkRERETiIE4109OAYWY2xMzKgHHAxBbaR+50qdqy\nBbZv9xdjyUT37rBpU/Ttrlrlk/KWpMo82vH/OSIiIiK7iZxMO+d2AlcDk4GZwOPOuVlm9i0zuxzA\nzPqb2WLgOuD7ZrbIzLpH3XapSdVLW4b/bnTr1nIy3fQQW3MyOQGxosJPn7diRWZ9K0WZxlMyU+h4\n1tb6f15LQaFjWWoUz7AUz3AUy3gKUjPtnHsRGNHksfsa3V8BDA6xrVK2YkXmJx9C68l0pjJJpqFh\nRo8BA6JvU6SQqqrgzDNh82Z47jk45phC90hERIpViDIPCWTlSn+BlEy1lkyniu5bk0mZB/hkuj3P\n6JFpPCUzhYrn9u0wbhzceis8+CCcf76/KFEx074ZluIZluIZjmIZT0qmYyTuI9MHHqi5pqX4Pf44\nDBoEl1ziR6cPOQQeeqjQvRIRkWKlZDpGVqwIOzKdaW3V2rW+Jro17f3CLapVC6tQ8fz1r+H66xuW\nr74aHnigIF0JRvtmWIpnWIpnOIplPCmZjpHQZR6ZWrcusxlE2nuZhxS/+fNhwQL40pcaHjv5ZFi8\n2D8nIiKSrcjzTIfWnueZPv98OP10uPDCzNq/+y78y7/4k6lyVVcHnTv7OtKOHVtuu3mzLwfZuLH1\ntiJxdPvtUF0N99676+Pjx8OoUXDVVYXolYiIFEKc5pmWQAoxMr1xo3+dTJLj8nJ/ouLixa23FYmj\n55+Hc87Z/fFTT4VJk/LfHxERKX5KpmMk9AmImdRW1dRkfpEYaN+lHqpVCyvf8dy61R/FGT169+e+\n9CV49dXinXda+2ZYimdYimc4imU8KZmOkdAnIGYi03rplPZ+EqIUr+nTYcQIf+XQpioq4IADYMaM\n/PdLRESKm2qmY2LHDujSxY+edcrwUjrbt/vSi9razK+a2NTrr8ONN8LUqZm1v+suWLgQ7r47t+2J\nFMqdd8LSpc3vu1de6ad/vO66/PZLREQKQzXTJWbVKujTJ/NEGqCsDDp0iHZoWiPT0l68/joce2zz\nz48eDW++mb/+iIhIaVAyHRPZnnyY0lKpRya1Vbkk03PmZN6+lKhWLax8xtM5eOON1pPpTI/QxI32\nzbAUz7AUz3AUy3hSMh0T2Z58mBK1bjrbZPqAA3xf16/PfZsi+TZnDvTo4a982Jz99vNHeZYty1+/\nRESk+CmZjonly2HAgOzXaymZTiQSra6fbTLdsSMcfDB88EHm65SKTOIpmctnPFsr8QB/3sHIkX7+\n9mKjfTMsxTMsxTMcxTKelEzHxNKlLY+aNSffI9MAlZXw3nu5b1Mk36ZObT2ZBp9Ma0YPERHJhpLp\nmPj0U9h77+zXy3fNNMBhh7XPZFq1amHlM55Tp8Jxx7XerliTae2bYSmeYSme4SiW8aRkOiY+/TT3\nkemNG3Pfbq7JdJRLmIvk04oV/gTfgw9uvW2xlnmIiEjhKJmOiaVLw49Mt0XNNMChh8LMmbBzZ3br\nFTvVqoWVr3i+8QYcc4yfRrI1Bx7op6lcu7bt+xWS9s2wFM+wFM9wFMt4UjIdE1FGpvNdM92zp5/G\nr71eVlyKS6YlHuAT7sMOK85SDxERKQwl0zFQV+dn8xg4MPt1C1EzDTBqFLz1VvbrFTPVqoWVr3hm\nMpNHY8V4gq32zbAUz7AUz3AUy3hSMh0Dq1b50d499sh+3UKMTINPTor1AhfSfmzZ4qdxPPLIzNfR\nOQEi+VFXB488AueeCxdcAM8/7y+wJFJszMVszzUzF7c+tbUZM2D8+NxGw370I//h8+Mf57bt3r1h\nwQKoqMhuvXffhYsugo8+ym27IvkwZQrccAO8/Xbm60ybBpddpoRapC1t2ADnnON/Xn01bN0Kv/wl\nfP7z8LvfQdeuhe6htAdmhnPOor6ORqZjINdp8SDayHRdnf8g69kz+3UPPRSWLIE1a3Lbtkg+ZFvi\nAf7LfM4cfzXEqOrqYPFiqK2N/loipaK2Fs44w191dOpUPzDzzW/CP/7hB4fGjWt/J7hLcVMyHQML\nF8KQIbmtG6VmeuNGv37Hjtlvt1MnGD0a/v737NctVqpVCysf8cz0Yi2Nde3qv+RnzYq27SlT/OuM\nGuX/Wb7nHp9ctwXtm2G1l3jW1cEnn8Czz8Kjj8KkSbBsWfjtNI3nhAlQXg733rvr90/Xrr7sY/Nm\nuP768P0oBe1l3yw2QZJpMxtrZrPNbK6Z3dhMm1+Z2TwzqzKzyhDbLRULFvgv3VxEGZnOtV465Z/+\nCV54Iff1pX3avBkWLfLTQUap929NXR28+Wb2yTREvzDRzJnw1a/Cfff5k4tffRUeewy+/nX/+0t+\nrFvn/+GfONG/nzpC4FVX+xLB/feHRAL+53/8Z/ldd/n52BMJePnlttl2Mgl/+pNPmtNNV1lW5p9/\n9ln/vokUg8g102bWAZgLjAE+BaYB45xzsxu1ORW42jl3mpkdBdztnDu6mddrdzXTX/uaPwHjvPOy\nX/fJJ+EPf4A//zn7dT/80G9z5szs1wX4+GP44hd9UpTJHL7S/jgH8+f7uZ5Tt3nzoE8f/6W5fLkv\nq7jkEn+YN5eTcJvz4Yfwla/4/TRbt9/uL/Ry113Zr+ucP2pz6aVw+eUNj2/b5n/HOXPguef89JLt\nQV2dP7di1iyfxG3c6D8v9tnHj9ofcghY5IrFXS1aBLfcAk8/7fevigo/AvvZZ3Dhhb6OPpepSIvZ\n5s3w1FPw0EP+H4sLLvD76MiRu7bbts23u+kmP6Xkvff6QZsQamv9bDm33urrpVvyxhu+zfTp7e+9\nagvO+bK33/zGl9GceWahexQPcaqZPhKY55xb6JyrBR4HzmrS5izg9wDOubeBXmbWTr5KWlddXZiR\n6VzrpVOGDfMj29Om5f4aUjqc88nK5Mnw05/6msgBA+CEE3zyOGIE3H8/1NT4evv58/3I4W23+VGx\nkSP9zBuh5FLikRJlery//MUnLt/85q6P77EH/P73/ojO0UeX7sm7O3b4z4Rf/ALOOgv22sv/U/P4\n476EoEMH32bKFDj7bD86etdd0a7k2tgDD8Dhh8O++/oE+vXX/f730Uf+RNTOnf05H9de2zYlDXFS\nV+dH5r/xDf/Py2OPwbe/7QdAfvWr3RNp8Pvp+ef7QZYOHfzf0MqVYfrzn//p35ezz2697ejRcNVV\nvp462/pp5+C11/zvPXKkf7/PPdePwNfU5Nb3YrV+Pfz61z4Gl10GRx3lB8EkrE4BXmMQsLjR8hJ8\ngt1Sm6X1j60IsP2gUl/yHTs23Dp12n158OBwo2htVeaRTCZbvFrShg3Qo0du20258EJ/5vVRR0V7\nnVw457+Ut2/3P3fs8B+6O3f6L5Gmt8aPO9fws+mBkNRjqTapdtOnJxk1KoGZH0lLNxqf7rUKvZx6\nrPHP5u6D/3137PCjSOl+bt/uk+C1a/0JqGvX+pPs5s716x96qN8fxo+H//5v/7eSTmr/PPlkGDPG\nJ5pjxvjRxFyT4MamTs39SyNV5uFc9qOmP/sZfP/76fcPs4bD68cf72tHr7pq15kLnPPnUcyY4W9L\nl/op/nr0gD339Mlpv37+1rev/0x6660khx2WYNs2PytC6qdzPnksK/O3bO537tzyESfn/GfPggU+\naZ0+3Y8mTpvmzwE5/ng/+vmb37R8gvU77/jE+667/M+vfz23keraWrjuOl+e8NprcNBBu7fZf3/4\n+c/hu9/1Rx8OPtiPzl5//a59TPfZuWOHf1+qq/1n56ZN/rEuXXa97bFHw/3OnRt+l+Z+du6863q5\nnMOSsnOnfy+mTYOXXvL/2Pbv7z+nf/KT7EZ4y8v9Z/stt8CJJ/q4DhiQW7+SySQHHJDgjjv89Qky\nfX9vusmXhdx6q/+7ycTatT5pfO89uPJK/89Dp07+n4Nnn/VHJU491b/vJ5/cfLyXLfMz+lRV+c+3\n9ev94716+Vvfvv5vMfWzZ0//d7ppk29bU+M/J1O3jRt9P1Lvd0WF/3tufOvb179OS/Fp7Xsd/N/m\nP/7h/7F84gn/e959t38fQx8FEi9EMh3c+PHjGTp0KAC9e/emsrLy/3aeVPF9Wy3ff3+Sp56CHj0S\n7NwJNTVJdu6E8nK/vGFDkh07oKYmwVFHwemnJzniiNy398ILSTZvhr59c1t/1qwky5cDZL/+hg2w\nZUuSZDL3/o8YkeTOO+EXv0hQXp5b/HfuhEGDEixYAK+8kmT9eujTJ8GaNTB7dpING/zvt2YNrF6d\nrE+aE9TWQocOSTp1gi5dEnTsCDt3JunYEbp2TdChA2zfnsQMunXzy9u2+e336OGT4k2b/PM9evj+\nbNzY8HyHDv55gF69EqxfD84lcQ66d09QV7dre/D7R+PXy2QZoGdPv7x+vX8+m+VU/wDWrfPPN13u\n3bthGaCiwv/+NTX++YoK/3xquV+/BJ06+eWOHf3706kTfPaZj/dBByXo0wdqa5PsvTdcd12CESPg\ngw/8+o3f308+yWx/uOQSWLEiyemnw9tvJxg+PNrf85QpMGZMbvv3CSf49//JJ5PstVfm6z/wQJJ5\n8+CrX225/cUXJxg9Gi6+OMmPfwyjRyfo0QPmzEmyYIF/f0eOhD59kgwYAJWV/u912rQkc+dCWVmC\nlSthwQK/P3bs6L+ct25N0rkzDByYoEsX/37t3On3h+3bYeVK//fTtatfrqnxyx07+r+nTZv8cl2d\nf77x31dZWcP+v2NHoj458P079NAElZVwyilJrr0Wzjyz4fedOxf23rvleDzxRIKpU308fv5z+OMf\nE+y/f+bv16GHJvja1/zn2c9/7vfP1t6vX/4Sjj02yWOPwcEH+3gPHpykXz8fz3fegTffTLJ0Kaxd\nm2DxYujd2/++Q4Yk6NYNVq1KUlvr368tW/z+u327f3+2bWv4e+vSxW9/yxa/vMceDcs7dgAk2LrV\nL3fo4N+fLl0AkpSV+b/fLl38+1tWBgMG+OX165P1//AlWLECPvkkSUWF35/GjIGTT06yzz65f76/\n+mqSk07yr3/qqXDbbUm6ds3+9cD/o3P66UmWLIFhwzJb/7XXklx1FVx9dYLjj/f7Y0vt77knyc9+\nBuedl+CRR/w/mZs2+edHjoR99kly8cWweHGC738fLrooSWUlHHecj+eMGUkWLYLqar//DxmSZNgw\nOP74BD17+u/bTZt8PrB4MUyalGTdOv/3sm6d//7p0gUGD07Qqxds3pykWzc45JAEffvC3Ln+/d57\n7wQffwwTJybrk/QEq1fDsmX+9xk6NMHgwdCxY5L+/f3n0b77wvLlST79tCrt7799O9x3X5Lp0/3n\nZ20tHH98kt/+tvXPo/a0XFVVRU394Ynq6mpCCVEzfTTwI+fc2PrlCYBzzt3RqM29wN+dc0/UL88G\nTnDO7TYyXSw101u2NPyXe/nl8IMf5PY677/vR28+/DC39WfP9odS58zJft2HH4ZXXvE/ozjrLDjp\nJH/YNFPz5/sRyL/9zY9m7bUXHHCAPxTZp0/DraJi1589evgRnLIMRs6kON1/vx+hnD4991rNhQv9\nhVqWL899JOZLX/L79GmnZb7Oddf5Pv/0p5mvs2yZ/123bPF/Bwcf7H8WmnN+pHP7dn+rrW2YLrB7\nd/97dgo8HFNb6+cavvNO+N73fPxb28acOXD66b504N//PbeR3a1b4a9/9SO6n37qf+++fX1JwgEH\n+JK2/fYLW9OfTupoW+OjC43vp3sM/FGKvfaCoUOjle611K/LLoMVK+CZZ7KP8Ysv+iMwH36Y2/zR\nL7/sR9cnT/ZHvpraudMfabjnHv/5ccYZmb3u7Nn+CNbChf7vb889/Xs9apR/7wsxirtunT9Cvnhx\nw23hQn8ewMKF/khVRYXvX9eufn9Ztsw/fvDB/rv47LN9mYxGoVsXqmYa51ykG9AR+BgYApQBVcBB\nTdr8E/BC/f2jgbdaeD1XTFascG7YMOfuvz+39Z9+2rnTT899+wsXOrfPPrmte889zl15Ze7bTnn/\nfef69XOupqbldvPnO/fTnzpXWenbX365c0895dxnn0Xvg5SWCy907uqrc1//4Yed+9rXovXhu991\n7mc/y7x9ba3fr+fNi7Zd8TE88UTnjjjCuenT07epq3PuwQed69vXuQceyG//2qPt250bM8a5a67J\nbr1Nm5zbbz/nJk2Ktv3HH3du4EDnXn1118fnzHHui1907oQTnFuyJNo2isGOHc4tXercG284l0w6\n99przs2d69y2bYXuWXGqzzkj58KRx/WcczuBq4HJwEzgcefcLDP7lpldXt/mL8ACM/sYuA+4Mup2\n46JfP39yy4QJ/j/HbM2aBZ/7XO7bjzrPdNSaafBn4597rp+Roek8unV1/oSs007z/+0vX+5PQvn0\nUz9t2Nln+xGgYtBaPCU7LcXznnv8CNiUKbm99quv+hMfo8j2suJTp/qa1GHDom03F6W2bw4b5kcj\nr7rKf3aMHQsPPujfjxkz/Ilkhx/uT6x65RX4l38Ju/1Si2cInTv72aNeesmfvJipH/8Yhg5NMnZs\ntO2fd57fB847z+8TN93kvz+OOcbPiPXyy+1j1o/XXvOldccc4z/jjjsODjzQH62VwglykM459yIw\noslj9zVZvjrEtuLoc5/zhyOvv97Pj5mNmTP9yQG5ijqbR/fuuW+7sbvugi9/2Zd8TJjgT7CYMqVh\nWqV//Vcfm/LyMNuT0lZR4Us9/vVffQlEtuUEU6ZkV3aUTmoKr0w9/XRmsxRIZsz8SWLnn++TuOef\n9ydRgZ/u7ic/8eUdOpSdP717+8GR0aP9SaZnNZ23q4k33/QnMf7mN2G2P3asP8HymWf8z3PO8a8f\n5XoJIiFErpkOrVhqppvatMmfKf7yy/6DPlMjR/pRllGjctuucz7R2LYt+4Tj2mt9n6MmHSnbtsF/\n/Iefo3RtvK4wAAAXt0lEQVTHDjjiCD+jwxe/qC88yZ5z/uzz886DK67IfL3Fi/3f1cqV0Wrqa2v9\nl/Rnn7Veu+2cr6l97jl/pEaklE2b5qd4fPppPzKazqpV/ujBf/2X5jSW+IrTPNOC/7L9znf8CTCZ\n2rnTnzyTbgqnTJn50d5crqoWYmq8xvbYw08JNn26n5bowQf99FhKpCUXZn4k8pZb/HRXmZo0yR8l\niXpyaufO/qhTJnNfV1X5f2az+UdapFiNGuUvFnb22emvgrt2rZ9+7p//WYm0tA9KpgO6/HJ/CCzT\nCe4XLPA111FLLZor9Wit7i9kmUd7oDrKsDKJ52GH+Qt+ZDM7xqRJ/os8hEwv3pIq8SjUP47aN8NS\nPFt38sl+Rqtvfxu+9S3/T+fq1b4k5/DD/SXJU3+3imc4imU8KZkOqKICvvpVP1F6Jj780E9lE1Wu\nddOhTkAUaUs/+Ymvi8zksuDbtvkT0r785TDbTl28pTXPPOOTfpH2ZPRof1SmosL/M3nAAb6s49e/\n9hfH0VFJaS9UMx3Yu+/6D5XUVRRbctNN/gzcTK/s1JxDDoFHH00//2ZLjjvOl6Xo0qISd7fd5v+2\nnnyy5XYTJ/q6/VxnAWnq9dfh3/7NX02sOZ984q/YuHRptKvXiYhIfqlmOqYOPxwGDkxfR9bUm2/C\n0UdH36ZGpqXUXXedP+nptddabvfYY/4iSKF84Qt++sqNG5tv88wzvi5UibSISPukZLoNXHll61MB\n1db6E/WOOir69lQznR+qVQsrm3h27eqPolx//e5zmaesW+frpb/2tTD9Az/F48iR8Pbbzbd55pnC\nT4mnfTMsxTMsxTMcxTKelEy3ga9/3R8W/uST5tu88QYMH+5rzaIqL89tZDr0bB4ibWncOP/z0UfT\nP3/fff4ywqEvAnTssb7cI50VK/yJVyedFHabIiJSPFQz3Ub+3//zc8/+/Ofpn58wwU+9lc1FIZoz\nbpyfPP/887Nbr7w8szl0ReJi2jR/oY5p02DffRse37LFXzXvL3/xJw2G9MIL/u843YDQb3/rrwj3\nxBNhtykiIm1PNdMx9+1v+xkI0tVa1tXBH//Y+tWjMpVLzfSOHX7mA12RUIrJqFF+Pvdx43bd52+5\nxY8gh06kwV845t130891/cQTfgYfERFpv5RMt5H994dTTvHTBDU1ZYpPYo84Isy2cqmZ3rTJ10tr\n6qLMqVYtrFzjecMN/mIqJ57oR42//31/qfp77gnbv5Tycj9n7qRJuz6+aJFPsuNwUQrtm2EpnmEp\nnuEolvGkZLoN/fCHcNddsHx5w2PO+Ynsr702XCKbS820Tj6UYtWhgy+vuPTShr+v11+H/v3bbptn\nnOEvzNLYww/7S5136dJ22xURkfhTzXQb++EP4Z134LnnfI30/ff7SyTPmOGXQ7j1Vl+ykc1V4mbP\n9heZmD07TB9EStnatf5o06xZMGAAbN7sL1AxebKf511ERIqPaqaLxA9+4C/Mcswxvs7zxz/2h6RD\nJdKQW820RqZFMldRARdd5K/GCP4iMscdp0RaRESUTLe5zp394eEf/hBOPtlfmvigg8JuI5eaaU2L\nlz3VqoVVbPH8yU/gxRdhzBj4/e/brkY7F8UWy7hTPMNSPMNRLOOpU6E70B507Ni2JymVl/vDztnY\nuFEj0yLZqKjwU/I9/zyceir061foHomISByoZroEPPUU/O//7n6CVEsefdTPhPDYY23XLxEREZG4\nUs20/J9caqY1Mi0iIiISnZLpEtDc1HiqmQ5LtWphKZ7hKJZhKZ5hKZ7hKJbxpGS6BHTrln3NtJJp\nERERkehUM10C5szxF5WYOzfzda6/HgYOhO9+t+36JSIiIhJXqpmW/5PrPNMamRYRERGJRsl0Ccil\nZlonIGZPtWphKZ7hKJZhKZ5hKZ7hKJbxFCmZNrMKM5tsZnPM7K9m1quZdg+Y2Qozez/K9iQ91UyL\niIiIFEakmmkzuwNY7Zy708xuBCqccxPStDsO2Aj83jl3aCuvqZrpLDnnr7S4ZUvmlyk/8UR/qfOT\nTmrbvomIiIjEUVxqps8CHq6//zDwlXSNnHOvA2sjbkuaYdZ8qUdzNDItIiIiEl3UZLqfc24FgHNu\nOaAL7BZIulKP1uaZVs10dlSrFpbiGY5iGZbiGZbiGY5iGU+dWmtgZi8B/Rs/BDjg5jTNg9RnjB8/\nnqFDhwLQu3dvKisrSSQSQMOOpOVdl7t1S7BpU+btN25M0KNHfPqvZS1rOffllLj0p9iXU+LSn2Jf\nTolLf4p5uaqqKlb9KbblqqoqampqAKiuriaUqDXTs4CEc26FmQ0A/u6cO6iZtkOA51Qz3TYOOwwe\nfhgqKzNr36sXLFwIvXu3bb9ERERE4iguNdMTgfH19y8Bnm2hrdXfpA1kUzPtnKbGExEREQkhajJ9\nB3CKmc0BxgC3A5jZQDN7PtXIzB4D3gCGm9kiM7s04naliWxqprdsgbIy6NRqkY801lw8JTeKZziK\nZViKZ1iKZziKZTxFSqecc2uAk9M8vgw4vdHyBVG2I63L5iqIGpUWERERCSNSzXRbUM10bi64AE47\nDS68sPW28+fDmDGwYEHb90tEREQkjuJSMy0xUV6e+VUQN27UHNMiIiIiISiZLhHpyjyaq63SHNO5\nUa1aWIpnOIplWIpnWIpnOIplPCmZLhGqmRYRERHJP9VMl4if/tSXedx2W+tt//xnePRReOqptu+X\niIiISBypZlp2kW5qvOaoZlpEREQkDCXTJUI1021PtWphKZ7hKJZhKZ5hKZ7hKJbxpGS6RKhmWkRE\nRCT/VDNdIp5+Gn73O3i2pQu617v5ZthjD/jBD9q8WyIiIiKxpJpp2UU2NdMbNqhmWkRERCQEJdMl\nIpuaaZV55Ea1amEpnuEolmEpnmEpnuEolvGkZLpElJerZlpEREQk31QzXSLmzoXTToN581pve9pp\ncMUVcPrpbd8vERERkThSzbTsIpvZPFQzLSIiIhKGkukSoZrptqdatbAUz3AUy7AUz7AUz3AUy3hS\nMl0iUjXTmVTIKJkWERERCUM10yWkc2efUJeVtdxu773hH//wP0VERETaI9VMy24yrZvW5cRFRERE\nwlAyXUKaJtPpaqvq6nybbt3y169SoVq1sBTPcBTLsBTPsBTPcBTLeFIyXULKy1u/CuKWLdClC3Ts\nmJ8+iYiIiJQy1UyXkMpKeOghGDmy+TYrVsChh/qfIiIiIu2VaqZlN5nUTKteWkRERCQcJdMlpOkl\nxdPVVmlavNypVi0sxTMcxTIsxTMsxTMcxTKeIiXTZlZhZpPNbI6Z/dXMeqVps4+ZvWJmM83sAzO7\nJso2pXndurVeM61kWkRERCScSDXTZnYHsNo5d6eZ3QhUOOcmNGkzABjgnKsys+7AdOAs59zsZl5T\nNdM5uvBCOPVUuOii5ttMmgR33w0vvpi/fomIiIjETVxqps8CHq6//zDwlaYNnHPLnXNV9fc3ArOA\nQRG3K2lkUjOtkWkRERGRcKIm0/2ccyvAJ81Av5Yam9lQoBJ4O+J2JQ3VTLct1aqFpXiGo1iGpXiG\npXiGo1jGU6fWGpjZS0D/xg8BDrg5TfNm6zPqSzyeBK6tH6Fu1vjx4xk6dCgAvXv3prKykkQiATTs\nSFrefblbN5g5M0ky2Xz7GTOSrFsHUPj+alnLWg6znBKX/hT7ckpc+lPsyylx6U8xL1dVVcWqP8W2\nXFVVRU1NDQDV1dWEErVmehaQcM6tqK+N/rtz7qA07ToBzwOTnHN3t/KaqpnO0W23+anv/v3fo7UR\nERERKXVxqZmeCIyvv38J8Gwz7R4EPmotkZZoVDMtIiIikl9Rk+k7gFPMbA4wBrgdwMwGmtnz9feP\nBS4ETjKzGWb2rpmNjbhdSaPp5cSbHmIDJdNRpIun5E7xDEexDEvxDEvxDEexjKdWa6Zb4pxbA5yc\n5vFlwOn196cCHaNsRzKT6ch0jx756Y+IiIhIqYtUM90WVDOdu2efhQcegIkTm29z7rn+9vWv569f\nIiIiInETl5ppiZGmU+OlozIPERERkXCUTJeQppcTV810WKpVC0vxDEexDEvxDEvxDEexjCcl0yVE\nNdMiIiIi+aWa6RLy8cfw5S/DJ58032bYMPjLX2D48Pz1S0RERCRuVDMtu8mkZnrDBujVKz/9ERER\nESl1SqZLSCY10+vWQc+e+etTKVGtWliKZziKZViKZ1iKZziKZTwpmS4hqZHp5qpktm2DnTuhS5f8\n9ktERESkVKlmusSUlcH69ekT5lWrYMQIWL06//0SERERiRPVTEtaTUs9Glu/XiUeIiIiIiEpmS4x\n3bv7kwxh99qq9et18mEUqlULS/EMR7EMS/EMS/EMR7GMJyXTJaZnz4ZkuimdfCgiIiISlmqmS8zo\n0fAf/+F/NvXcc3DfffD88/nvl4iIiEicqGZa0urZ05dzpKOaaREREZGwlEyXmMbJtGqmw1KtWliK\nZziKZViKZ1iKZziKZTwpmS4xGpkWERERyR/VTJeY73wHBg2C66/f/bmbbvJT533/+/nvl4iIiEic\nqGZa0tLItIiIiEj+KJkuMaqZbjuqVQtL8QxHsQxL8QxL8QxHsYwnJdMlRiPTIiIiIvmjmukS88c/\nwp/+5G9NnXQS3Hyz/ykiIiLSnqlmWtLSyLSIiIhI/iiZLjEt1UzrcuLRqFYtLMUzHMUyLMUzLMUz\nHMUyniIl02ZWYWaTzWyOmf3VzHY7vc3M9jCzt81shpnNNLPbomxTWtbSyHRNDVRU5Lc/IiIiIqUs\nUs20md0BrHbO3WlmNwIVzrkJadqVO+c2m1lHYCpwvXNuajOvqZrpCBYtgmOPhcWLd328rg7KymDr\nVujUqTB9ExEREYmLuNRMnwU8XH//YeAr6Ro55zbX392jfptrI25XmtHcyPSGDVBerkRaREREJKSo\nyXQ/59wKAOfccqBfukZm1sHMZgDLgaRz7qOI25Vm9OgBGzf6kejGtVVr1kCfPoXrVylQrVpYimc4\nimVYimdYimc4imU8tTpOaWYvAf0bPwQ44OY0zdPWZzjn6oCRZtYTmGxmJzjnXm1um+PHj2fo0KEA\n9O7dm8rKShKJBNCwI2k5/fJrryUpK4NNm3Z9vmfPBH36FL5/WtaylsMvp8SlP8W+nBKX/hT7ckpc\n+lPMy1VVVbHqT7EtV1VVUVNTA0B1dTWhRK2ZngUknHMrzGwA8Hfn3EGtrPMDYLNz7hfNPK+a6YgG\nDYK33oLBgxse+9vf4Pbb/U8RERGR9i4uNdMTgfH19y8Bnm3awMz6pmb5MLOuwClAVcTtSgv69IG1\nTarS16zRTB4iIiIioUVNpu8ATjGzOcAY4HYAMxtoZs/XtxkI/L2+ZvotYKJz7uWI25UW9Onjk+fG\nh9hUMx1d00OWEo3iGY5iGZbiGZbiGY5iGU+R5nZwzq0BTk7z+DLg9Pr7HwCHR9mOZCeVTDdOnpVM\ni4iIiIQXqWa6LahmOrpvfAOOOQa++c2Gx264Afr18z9FRERE2ru41ExLDKVGphtTzbSIiIhIeEqm\nS5BqptuGatXCUjzDUSzDUjzDUjzDUSzjScl0CWpuZFrJtIiIiEhYqpkuQX/6EzzxBDz5ZMNjhx4K\njzzif4qIiIi0d6qZlmb16QOrV+/62KpVGpkWERERCU3JdAnaay/47LOG2qqdO/1y//4tryctU61a\nWIpnOIplWIpnWIpnOIplPCmZLkEDBsDy5Q3Lq1dD797QuXPh+iQiIiJSilQzXYLq6qBLF9i4EcrK\n4P334cIL4YMPCt0zERERkXhQzbQ0q0MHX+qxcqVfXrbMj1aLiIiISFhKpkvUgAHwwgtJwJd8DBxY\n2P6UAtWqhaV4hqNYhqV4hqV4hqNYxpOS6RI1YEDDXNPLl2tkWkRERKQtqGa6RH3jG3D00XDZZXDN\nNbDffnDddYXulYiIiEg8qGZaWjRkCCxc6O9//DEMG1bY/oiIiIiUIiXTJWr4cHj99STgk+kDDyxs\nf0qBatXCUjzDUSzDUjzDUjzDUSzjScl0iRo+HJYsgR07YNEiX+YhIiIiImGpZrpEbdjgTzqsqoIv\nfQkWLCh0j0RERETiQzXT0qIePaCiAp5+Gj7/+UL3RkRERKQ0KZkuYSNGJLnxRjjhhEL3pDSoVi0s\nxTMcxTIsxTMsxTMcxTKeOhW6A9J2zj0Xtm6Fiy8udE9ERERESpNqpkVERESk3VHNtIiIiIhIgUVK\nps2swswmm9kcM/urmfVqoW0HM3vXzCZG2aZkTrVVYSmeYSme4SiWYSmeYSme4SiW8RR1ZHoC8Dfn\n3AjgFeB7LbS9Fvgo4vYkC1VVVYXuQklRPMNSPMNRLMNSPMNSPMNRLOMpajJ9FvBw/f2Hga+ka2Rm\n+wD/BPw24vYkCzU1NYXuQklRPMNSPMNRLMNSPMNSPMNRLOMpajLdzzm3AsA5txzo10y7XwI3ADqz\nUERERERKRqtT45nZS0D/xg/hk+Kb0zTfLVk2s9OAFc65KjNL1K8veVBdXV3oLpQUxTMsxTMcxTIs\nxTMsxTMcxTKeIk2NZ2azgIRzboWZDQD+7pw7qEmb24CLgB1AV6AH8JRzLu3sx2am0WsRERERaXMh\npsaLmkzfAaxxzt1hZjcCFc65CS20PwG43jl3Zs4bFRERERGJiag103cAp5jZHGAMcDuAmQ00s+ej\ndk5EREREJM5idwVEEREREZFikbcrIJrZWDObbWZz60tC0rX5lZnNM7MqM6vMZt32Jtd4mtk+ZvaK\nmc00sw/M7Jr89jx+ouyb9c/pgkSNRPxb72VmfzKzWfX76FH563k8RYzn9+rj+L6ZPWpmZfnreTy1\nFk8zG2Fmb5jZVjP7Tjbrtje5xlLfQ+lF2Tfrn9d3USMR/9az+y5yzrX5DZ+0fwwMAToDVcDnmrQ5\nFXih/v5RwFuZrtvebhHjOQCorL/fHZjTnuMZJZaNnr8OeASYWOjfp9C3qPEEfgdcWn+/E9Cz0L9T\nscazfp35QFn98hPAxYX+nYognn2BI4Bbge9ks257ukWMpb6HAsaz0fP6LgoUz2y/i/I1Mn0kMM85\nt9A5Vws8jr/gS2NnAb8HcM69DfQys/4Zrtve5BxP59xy51xV/eMbgVnAoPx1PXai7Ju6INHuco6n\nmfUEvuice6j+uR3OufV57HscRdk/1wPbgW5m1gkoBz7NW8/jqdV4OudWOeem42egymrddibnWOp7\nKK0o+6a+i3aXczxz+S7KVzI9CFjcaHkJu//hNNcmk3Xbm1ziubRpGzMbClQCbwfvYfGIGktdkGhX\nUeK5H7DKzB6qP1T5P2bWtU17G385x9M5txb4BbCo/rEa59zf2rCvxSDK94m+i3YVJB76Hvo/UeOp\n76JdRYln1t9FeauZzoEu7tKGzKw78CRwbf3IgGTJGl2QCL+/ap+NphNwOPBr59zhwGag2ak2pWVm\ntj/+sO8QYG+gu5ldUNheiTTQ91AY+i4KLuvvonwl00uBfRst71P/WNM2g9O0yWTd9iZKPKk/5Psk\n8L/OuWfbsJ/FIEosjwXONLP5wB+AE83s923Y12IQJZ5LgMXOuX/UP/4k/gOtPYsSzy8AU51za5xz\nO4GngNFt2NdiEOX7RN9Fu4oUD30P7SZKPPVdtLso8cz6uyhfyfQ0YJiZDak/m3wc0PRs04nAxQBm\ndjT+kOSKDNdtb6LEE+BB4CPn3N356nCM5RxL59xNzrl9nXP716/3imvmyp7tSJR4rgAWm9nw+nZj\ngI/y1O+4ivK3Pgc42sy6mJnh4zkrf12PpWy/TxqP8Om7aFdRYgn6Hmoq53jquyitKPHM/rsoj2dW\njsV/uM8DJtQ/9i3g8kZt/gt/9uV7wOEtrdvebznEc2T9Y8cCO/Fnts4A3gXGFvr3KbJYHp7mNU5A\nZ1BHjidwGP5DsAo/ktqr0L9PoW8R43kDMBN4H3gY6Fzo36fQt9biCfTH11rWAGvwNefdm1u3Pd9y\njaW+h8Lvm41eQ99FAeKZ7XeRLtoiIiIiIpKjOJ+AKCIiIiISa0qmRURERERypGRaRERERCRHSqZF\nRERERHKkZFpEREREJEdKpkVEREREcqRkWkREREQkR0qmRURERERy9P8BXmZaVwE7urQAAAAASUVO\nRK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1,1, figsize=(12, 4))\n", "ax.plot(time+tlag, data[0][0,:])\n", "ax.grid(True)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": false }, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "hide_input": false, "kernelspec": { "display_name": "pymain", "language": "python", "name": "pymain" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.9" }, "toc": { "toc_cell": false, "toc_number_sections": true, "toc_threshold": 6, "toc_window_display": false } }, "nbformat": 4, "nbformat_minor": 0 }