{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "To obtain the data to run this tutorial, download the tar.gz file from https://zenodo.org/records/10668525 and unzip it. Then set simulation_directory to point to the resulting directory." ] }, { "cell_type": "markdown", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true }, "pycharm": { "name": "#%% md\n" } }, "source": [ "# Reading raw simulation files\n", "### Simulations do not have to be fully stitched into an h5 file before reading and visualzing their data" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [], "source": [ "from mayawaves.utils.postprocessingutils import get_stitched_data\n", "import matplotlib.pyplot as plt\n", "import numpy as np" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [], "source": [ "simulation_directory = \"D11_q5_a1_-0.362_-0.0548_-0.64_a2_-0.0013_0.001_-0.0838_m533.33\"" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Use get_stitched_data to stitch together any output file structured as columns of data with each row being a time/iteration step" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Stitching 17 files\n", "Stitching 17 files\n" ] } ], "source": [ "shifttracker0_data = get_stitched_data(simulation_directory, 'ShiftTracker0.asc')\n", "shifttracker1_data = get_stitched_data(simulation_directory, 'ShiftTracker1.asc')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### You can then plot the data based on the columns present in the data file. In the following, columns 2 and 3 (starting at 0) of the shift tracker files contain the x and y components of the positions." ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAASAAAAEKCAYAAACytIjQAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAAsTAAALEwEAmpwYAACFN0lEQVR4nO1ddXhT5xs9X720tAVa3J3i7u6MDdnGhG2wjY1tzF1+c3eHDQYbMBg23N21uHupU3dN8v3+OLmL9KZN2tTgnufJkza5ublJ7j3fK+d9XyGlhAYNGjSUBVzK+gA0aNBw+0IjIA0aNJQZNALSoEFDmUEjIA0aNJQZNALSoEFDmUEjIA0aNJQZypSAhBBzhBCxQogzZo9VFUJsEUJcNt5XsfHaEUKIi0KIK0KIN0vvqDVo0OAslLUF9BeAEVaPvQlgm5SyGYBtxv8tIIRwBfArgJEAggE8IIQILtlD1aBBg7NRpgQkpdwNINHq4TEA5hr/ngtgrMpLuwG4IqW8JqXMBbDI+DoNGjRUILiV9QGooIaUMhoApJTRQojqKtvUARBu9n8EgO6F7TgwMFA2bNjQKQepQYMGdRw9ejReShlkz7blkYDsgVB5TLWmRAjxJIAnAaB+/foICQkpyePSoOG2hxDihr3blnUMSA03hRC1AMB4H6uyTQSAemb/1wUQpbYzKeVMKWUXKWWXoCC7SFmDBg2lhPJIQKsBTDL+PQnAKpVtjgBoJoRoJITwAHC/8XUaNGioQCjrNPw/AA4AaCGEiBBCPA7gCwBDhRCXAQw1/g8hRG0hxHoAkFLqADwLYBOA8wCWSCnPlsVn0KBBQ9FRpjEgKeUDNp4arLJtFIBRZv+vB7C+hA5NgwYNpYDy6IJp0KDhNoFGQBo0aCgzaASkQYOGMkNF1QFpqOgw6IGMeCD9JpARB+RmGG/pQF4moMvldv8pvgTg5gV4+gIexpunL+BTHahcA/D0A4SaPExDeYZGQBpKDpmJQNwFIOEqkHjNdEuLJvmoa0eLBjdvwLc64F8XqNoIqNoEqNbEeN8UcPdy3ntpcBo0AtLgHKTHAuGHgOhTQMxpIOYUkBppet7FDQhoQHKo3RHwrUHC8K0O+AQBnpUBDx+jdeMDuHqYXislAAnosoGcdKOllAbkpAHpcUB6DJAWQ2sqJQK4vAVI/9vyvYNaAjXbAbXaAbU68Bg0UipzaASkoWhIDgdC9wA39gNhB4CEK3xcuACBzYEGvXjB1wimFeJfD3A1nm76PJJFWgytoZtngaxkICeVLlhOGklGGow3IwG5uAHu3oC7D+BRCXCvBFSqSjeseiujO1YL8AnkPhKv8bhizwHRJ4ErW4CTC3kMrh5AnS48zga9gPo9SHwaShXidhrL06VLF6nVghUR+jxaOJc20cKIO8/HvfyB+j2NF3FPoGZbkgRAqyj6FBB/iUSQcIXuWGqE+nu4eRljO0ZryMUVgDDGdgRg0AF5WYwRKTEjQ17+/bj7GN2wRkCVRkbrpy3vMxOAqGMkzhv7SUxST0Jq2AdoPgJoNoyv1VAkCCGOSim72LWtRkAabEKXC1zbCZz5F7i4AchJAVzcSTbNhgGNBwDVgwEXF7pG4YeAiCO8qKNOAGlm5Xme/kBgU6BaM6BKQ8CvFq2VyjUB35qAdxXAzUP9OGxBSlo6GXEku4xYIDUaSAql9ZN0nX/rjQFtF3eSUK12QL1uJMzKtYCIw8CV7cDlTSZLLrAF0Hos0OZuIKhFMb/I2wsaAdmARkB2QEogdC9weglwbjWQnUwrp+VooMVIoFF/wMsPyMsGbuwFru8BbuwDoo7TQoEAApsxzlKrPW/VWwGVqpVNlsqgJxnFnDbdok+QtADAuyqJqFE/oNlQPnZ5M3BhHb8HSKBGGxJR23uBgHq23kmDERoB2YBGQAUgPQ44sQA4NpcXrIcv0GIUL7wmg2idpMXQBbu0kZZRXibjMnU6Aw1604Wp140uVHmGlPyMSvzqxn5aSwBQtTGtu2bDaPmcX0sLMOIwAEGS6jwZaDbcFNPSYAGNgGxAIyAVhB0CDk7nim/IozXQaRIQPIaB3ox44OwK4PRSulgA4FcXaDGC8ZIGvW6N4G3iNeDyVgaqr+9mxs27KtDqTqD1OCCgPnByEXB8PgPnlWsBHR8Guj3BTJ6G/6ARkA1oBGSEQQ9cWAvs/4Uru1cA0PEhoNMjXPXzsoHza4BTi4Gr2xmkDWpFa6jFSKBG65Jxp/KygaxEIDsF0OUw8K3PNcZwJODqyWCxqzvvvfyZBVOC3k47jizg6g7g7HLGvnLTKRVoOwHoOBFIugEc/ZPBeFcPoP19QM9ntViRERoB2cBtT0D6PODEQmDvdwzOVmnIC6fDg7RiEq4CIXPoimUlMXXe9h7GPmq0LuZ76+jmxF8GksOAlHAg+QbT+RnxJJ68zKLt283bmI4PohAxoAFjNQH1KQmo0qjo7lJeFonm9FKSkSEPqNedVmKtdsbvayEtpuYjgP6v0yW9jaERkA3ctgSk1zGovOtLEk/tTkCfl4CWd1C3c3kLcOAX4PouxnRajga6PAY07MsMl6PIy2JQOiKEGp/Yc0DcRUCfY9rGzYsE4V+PokTvKkClKnR7vAP4vKuH6QYYrSGjVaTLpqWUlUjFdVYStUXJ4SQ3czJz9QSCmjNjV6M1ULcbULuD45ZTehxw8h/g2Dwg4TItx66P0zI8vwY49DuPp8UoYODbTP3fhtAIyAZuOwKSEji3Etj+CdPLtdoDA99hgNWgY3B1348kCL86QJdHgY6PsLbKEWSnMG4SupdxopjTxowYgMq1KUas3ooEENiCxOMTWHJZMSmp90m+AcRdAmLPArHneVPU2S7uJIh63YFGfUm2Xn727//GPuDgDMbOXN1pJXZ+lMH5/T9TshA8Bhj0HuUHtxE0ArKB24qAoo4DG99ilqd6MFfklqNpQRybC+z9gYLA6sFA7xe4iru627dvg4H7v7IVuLqNlo7UU5lcuxNQryutjLpdAd9y1oc7PY5apYjDQPgRIPIooMui5Ve3G9B0ENB0KMnaHoJMuEoiOrGAVlfrcUD3pxnMPjiDllq3qXTNvANK/OOVB2gEZAO3BQGlxQDbPmJcolI1YPC7zNYADCrv/JwxmHo9gL6vMK1sz4VmMNC6ObcKOL/aaEkI1lQ1GcRbvW72k1h5gS6Xn+vqNgbco0/y8SoNacEEjyGpFvYdZSYCB34FDv1GhXbbe+jGnlgIHP+bMaqB7zCF7+Ja0p+qTFHhCUgI0QLAYrOHGgN4T0r5g9k2A8CG9UYBB5ZLKT8qaL+3NAEZDMCxv4AtH3Al7vE00O9Vtqm4uB7Y+iEQf5Er++D3gCaD7SOe2AtMPZ9exqJPV0+g6WCg1V105XyqlfQnK12kxwGXNpBor+2kKxlQH2j/ANBhIlClQcGvz0gA9v8EHJ7JTF7XKUDLUcCur+i21eoA3PUzA9i3KCo8AZnDOIY5EkB3KeUNs8cHAHhVSjna3n3dsgQUdwlY8wIQtp+xjNE/MO4QdxHY+CZX9sDmwKD/kTgKI57sVMaHjv8NRIbQPWk+gu5F8+HlX2joLGQmMvN1einJCAAa96dF2epOwM3T9mvTY4GdXzBd7+kHDHiTAfbN73C/vZ8H+r/hfAlBOcCtRkDDALwvpext9fgA3O4EZNADe79ndsu9EjD8U67SOalccQ/9xsLMgW8zW1OYe5R4nZmc43+z3UVQK6DTw9S/lGQsx6CnuC89lin5jDjeclKpDdJlMbOmy6H7IlyN9y6UD3j58yL38mc2za8WxZI+QUXL4qkhOQw4voCxnpRwZu66PQF0fqxgK/DmOWDTWySwwObAkA9okR7/m6rrMdOBBj2dc4zlBLcaAc0BcExK+YvV4wMA/AsOKYwCyajA0Ty3FAElhQLLpwLhB4HgscCor6nIvbgBWPsSY0GdHqG75RNY8L5u7Gf84sI6Xtht7ga6PUk9izMzVVlJQMwZ4OYZ6oGSQqkNSg4zZc3MIVxIrG5exnsPtucw6E33eRm02NSam7m4MQtXrTGzb0HGW/VgxmSKAoOeIsWD0xk3cvOie9brOTZAU4OULF/Z8Do/a6dJQNMhwJZ3+X+fl4ABb1W8+JkN3DIEZBw6GAWgtZTyptVzfgAMUsp0IcQoAD9KKZup7MN8NHPnGzfsnhpbPiElg8nrXiU53PEtU8CZCcCGN4Azy4DqrRlnqFuIIO76HlpPoXvoHnR5jDELv1rFP05dDiviww4A4YfZoCwl3PS8V4CxXUZDCgWrNKBV4RNEwqwUSOvG3gB5bhrlAJkJrIhPjeQtJYIShLhLJCsFVRoBdbuQZOt0oS7IUQKIPU8iOrmYBNr+fsbdqjZW3z43g0mAA7/ycw75gHGh438z0D1+1i2Rsr+VCGgMgGlSymF2bBsKoIuUMt7WNhXeAsrNpHVzahFQvxcw/ncGSM+t4uPZqUC/17iiFtTaInQvsONzVrP71uD2nSax9quoMBhYZX55CwWNkUeZggbYErVWB6BmG2pvarR1XGtUXEhJQoq7QJ1S5FEg4qipZYhHZRbTNh7AW1AL+62/tJvUU4XMpsyh/QPAgDf426gh6gSw5nlm3No/wLjd5ndI2qN/YGlHBcatRECLAGySUv6p8lxNADellFII0Q3AMgANZAEfqEITUMJVYPHDFA0OfJspdF02g8zH5jEdPnYGBX+2EHeJZv+ljezB0+cloPOkogdCczPYuuLSJmqCMuIACGbaGvRmbKNej/KnBTJHahQttOu7GKdJvMbHqzRkoLnVXbSQ7Iklpd0E9v3A8gwpgZ7TgL4vqwft9XnA7m+A3V+xfGTIh8CRP2gRdXsSGPap4/2RygluCQISQlQCEA6gsZQyxfjYUwAgpfxNCPEsgKcB6ABkAXhZSrm/oH1WWAK6sA5Y8RTjM3fPZho8+hTw7+OMpfR5kRoTWy5ERgKw6wvgyGzGUvq9AnR/qmjEk5dNkd2Z5SSyvEy6b00HMy3fZFDhMafyjKQbzBpeXM9YjyGPZN16HGvm7Emfp0QC2z6kq+xTndnHjg+p63/CDwPLnzDFgvKy6NbV7QZMmAv41Xb+Zyxh3BIEVBKocAQkJbDnG5ZS1O4ITJhHsz7kTwY0K1UDxv3O1LAaDAaqnre+z46FnScz2FkUiyT6JHB0LlPSOal87+AxQOvxbMlxK4rrslPoUp5bSStPn0v3scODQLv7CtdARRxlBiz8EK2oO3+kG2qNnDRg/WusM2syiIr1ze8yBvbgogpX3KoRkA1UKALS5QJrX2Tat+29wF2/MCax/jWSStMhwLiZti+C2PPAmheZJWvYFxj1DVC9pWPHkJsBnFrC94s6zoxP8Fig3QR2RrydGnJlJlIbdWIBvwtXT34PPZ4uuFOAlCTtjW+xu2Sv51mWYW19Sgkc/YsLi28NWrQ7P6c04e5ZdAcrCDQCsoEKQ0BZSYz3hO4B+r9JEVtaNB+LDGH8Z+A76laHLpeZrX0/MPYw7FOu2I6k09NuUsl75A9eNNVbM1bUbgJ1Nrc7bp4DjswCTvxDjVKjfiSWpkNsf8+ZicDm/5HAqjYGxv4G1O+ef7vIo8CSSazsH/IBm8FFhADDPmbrlAowfFEjIBuoEASUGg3MH8ug85hfmNqNPgksmMDGWGNnAMF3qb829jzjCTGngXb3U5joSDwm4SqJ6+QiBklbjeZJX6972Z74SvN5XbaxUZmxSZmUVCO7utMicfOkINFZ4sPCkJlI6/DwLGbYanfiYtFsmO3v69ouYPWzlAf0e40369hdZiKw5BEuQL1fpFbq3CrG7YZ/Xnqfr4jQCMgGyj0BJYUC88ZQDXz/QsZ2Lm/lyVipKjBxqXqWy2Cg6nnrB7R67vqJvX7sft8bzMac+IcXQ4eJzODYEtY5EwYDkBxKFXZymLFJWRiJODOB/XWyktSFimoQLrTSvKvyO6tcixqjgAamwYhVGjn3ItbnMX6z+xsef+2OwMD/Ac2GqG+fnUpX6+Q/jO+Mn5X/u9blsrzm5EIq0b38aJF2eIi/bzmOuWkEZAPlmoDiLpF88jKBh5ZTRHhsHuM4NYKBB5eqCwQzEoDlU5i5aT6CAkR7exSnxbBk49g8XrhdHmMmpqQ0OrocWnORR6mGvnmOuhzz5mEubmxS5leHBFKpGu+9qzBu4uphbM1qtBrM27bqsklWmYkkr8wEptlTwk2jeQCWp1RvxdiN0hOoRuviX9T6PFqPu78mETUZDAz7hL+fGs4sp37LoAfGzcgf5zFPQjTsS3nDgV8Yhxs/q9ym6TUCsoFyS0A3zwFz7yQJPLKSF8Oeb9lWo8lgpmPVtCQRIYwXZMQBIz4ngdjjKuVlM9W751tetJ0eAfq+CvjXce7nys2g6DF0L9PNUcdNXRErBRoblbXmfbWmzPBVruX81d1gYAwt+QZlC7Hn2Kkx5jRjXACFiMqssKaDgFodi24l6XIZI9r1JV3HzpNpEaklDFIiaOFGHqW7Nejd/MH9U0sow6jdAWg8kKTUbBhw398FF8SWETQCsoFySUDxl4E/R/Gim7yOAcodn9ElajsBGDs9f4xASprjG9+iVTTBKEQsDFJSU7T5Hbp7Le5gcNNZrpaUjENd2QJc2cYyDH0urZZaHRh0rdedjcoq13TOexYHUtLdCz8EhB3kLfYcAEmCbDqE/ZKK2gEgM5EV8Uf+YDOy4Z8zkG+9SOhyKCgNmUNL5965+cnqwnpg6WT+VsFjgZ2f8febMLfc1ZBpBGQD5Y6AkkKBOSN5kT66gQP9trzLlp4dH6ZuxNoa0OcB616m29RsGHVA9hRWpkQAa1/m9M+glrSYmgwq/meQkhft2RXA2ZXslQyw4LPpYFpw9XsC7l7Ff6/SQEYC3dnLm6nuzkqk/KD5cKDNPfzOHf0sN8+x9CLiCC2Y0d+rj34+sZAut38dYOKy/AvDtZ3APw+QvIPHcrhA6/HA3X+Uq5iQRkA2UK4IKCUS+HMETfRJa+l2bXqbrlHXKcDIr/O7AFnJNNev76LLNPCdwt0Eg4E1Sls/ZNvUge8wm1JcDU9aDFPKJxdx9rtwYflF67FA85HOd+fKAgY9Xcezy0mwGXHMsrW9l+6umqjQ5r7MfgeDjhlKNZc5/DDwz/2s9r//n/ytOsIOAX+Pp0K6+Qg2P2v/IDDm13KTHdMIyAbKDQFlpwBzRnCCw6TVQJ1OwPZP6XZ1f5rWifWJmXgdWDiB93f9RG1PYUi8Dqx8mq5Q44HAnT+wxqmo0OtoGRyfT2Ww1LMotu3drJm6lQf06XVMi59cRDLS57BcouvjtELsDQinRDINf3U7iXrML/mlEonXgAX30j0c9zvQZrzl86H7gL/vpoXUsC9waAbjR0M/dMpHLS40ArKBckFAulxgwT0sOnzoX1Ze7/+ZIrWODzOLZU0+MWeoDdLnAfcvYNV2YTi5GFj3Ci2TkV+w6rqoWp7sVJLOwd+AlDDWN3V4kMd7C7SPcBiZiUyhh8xhqw+/OtRLdZ5k35RYc9mElz+7Gli7w5mJwKKJXDzu+pmN4cxxdTuw8H5jAL8Zxy7d8S2t5zKGRkA2UOYEJCUtkpP/UAnb4QHWV615nsWOd8/O78uHHwEW3M1Z7Q+v5HyrgpCdwl5Bp5cw9jJ+pu22EIUhNYou4dG5rP+q34ulBy1GlrvAZ5lASsaJ9n7PBcW7CtBjGr8jT9/CX3/zLPDvFAbuB78L9HnZcpHIzQQWTyTZjPwa6P6k5esvbuTzDXpxzNC1HdSPtRjp3M/pIDQCsoEyJ6BdXwE7PgUGvM1+MZc2A//cx0Dt/Qvzm/HXdgL/PEhdziOrCieSm2e5aiaHsei078tFC06mxQB7vmM/Y4OecZ2e0ypcUWSpIuwQiejSBjYb6/8G0++FEXVuBrD6OdaZtRxNpbv5fDJdDrD0UeDiOrbs6POi5euP/w2smkY3MPEa43GPbSrTpvcaAdlAmRLQxY0km3b30a+/eRaYM5zZkEc35l8xFRO7WhNaPoWJA88s54no6cfUbP0ejh9jehwvIqWxVseJDHYXNglCgwnhR9h94MY+SiqGfmycQFuA+yslLc3N71IPNXGp5XeuzwNWTCVJjfiCFpY5dnxGzVG3J4Hza5lgeHJX0dvOFhMaAdlAmRFQwlVg5kCeVI9vZkxl1iAGcadsy58x+i/I2JRB6oJOJIOesYT9P7H514S5jmtsdDmMSez6mm1L2z/AGiW1VLGGwiEl23hseQ+IO89s1cgvC08AXN8NLH6ISu+JSyy1XXodsHQScGEtOyOYx4TMXfuez7KQuEEvYOK/ZdKxQCMgGygTAspJB/4YzLYKT+4kOfw5kn7/oxuobjVHRAhLMvxqA5PXF9y7JzcDWPY4zf6uUyh0c0Se/58w8X8seGw2nKUDhcWZShu6XJZTpN+0LLPITjGWYigFqjpecEphqqsHg7xKj+lK1fi9Vq5VOilrfR6JfcfnTKv3fw3o9ULBpBB3Efj7HiAzHrj3L+qPFOhyqAO6uh24ZzaHB/z3XC4w7y6WunSeTIuq1/MUmpYyNAKygTIhoOVPsh/MQ8uBJgMpBgyZDUyYn7+qPeYM8NcoFlI+uqHg5vDpcUzLR58ARn7FETGOIDmMx3JlC4WJwz+l8rcskRHP8oibZyluTLzO8onUKKhPvXA3EY2bJ+vIFELSGWvDpD7/61w9aY1WachROTXbsSYssHnJWAwpkVQ6n1/NONq4mQVnD9Ni+NvGnKG73u5e03O5mcyihh8GHl7Bufb/ve4m8Hs/1szV6USX7YHFQIsRzv9MBUAjIBsodQI6tYTtMQa8xTYNJxcDK55UX5lSIoE/jATw+KaCA84JVylGS7vJldCRyneDnu0jthmHyA76H2MHpW2q63K5WocbSyAijtDCUeBTnfGvgAamana/Wsbi1Gok6cKa6EvJFiaZCVQ4Z8ZTEZ4UahoJFH/Z1Dzf1ZMWaYPeQMPedGntyWbZizPLqWLPy+bv33WK7dhQTjoFiaF7WY5jrvvKSgZmD+P3NWUrFfQKwg+ztKdBL37utGjg6f2lWvpySxCQccpFGgA9AJ31BxJCCAA/AhgFIBPAZCnlsYL2WaoElHgd+K0vFc6T1zE7MWsQV6ZHVlte8NmpdMuSbgCPbeBqbAs3z9FFk3rgwSUcLWMv4i+zqDEyBGg6FBj9XdFT9I5CSn4Hl7fQ6go7aLrwqzRiIWjNdvy+arQpvUb2eh3LR6JPcXRQ2EEWzUo9ByDW78m0dstRtsftOILUaGa9rmwxZr2m001UQ24msOhBZkPv/IGulYKkUGDWYBLklG2WYsbDs4D1r7KtytkV/G4fWlFqSulbiYBsjtkxzgJ7DiSg7uBcMJUWcyaUGgHpdSyziLsEPL2Xadnf+7Pyeuoey4yWPg9YeB9PsolLWT9lCzGnST4u7sCkNfbHaqRk46yNb9FVGfk10Paekm8yZjCw0POssYF9chgfD2xB4Z0yNaO0R/QUhpx0HnfoHiq+Y8/x8aBW/N7a389JFkWFlMDBGaz7C6hPd9xWWUdeNrDkYSrQx/zK5vYKwo8Ac0czWD1pjSnlL6UxYL2OnQ5C5qin8EsItwsB/Q5gp5TyH+P/FwEMkFJG29pnqRHQnu84FeHu2Txh178OHP6d6fQmAy233fgWA4Z3/cyTxRYijwHzx1GQOGm1/RXsmYlccS+sZR/ncb87Z/CgLUhJC+LMv1x9UyNZzNl4ICvLmw0tPavLWUi8zomz51dTmQzBNqwdH2JRaFH78tw4wAr37BRg7K+WQWVz6HLojl3bSbJqZTaN/PQyTkfpMQ0Y8Znp8cxE4Lc+jAdVacgM29Q9jvcFLwIcISBIKcvlDcB1AMcAHAXwpMrzawH0Mft/G0hYNvfZuXNnWeKIuyTlR0FSLnqI/1/eKuX7flKufyP/tieX2H7OHNGnpfy8npTft5EyMdT+Y4kIkfLbVlJ+WE3KfT9Jqdfb/1pHkZUs5aGZUk7vzc/0YTUpF9zHz5idWnLvW9pIuCbl9s+k/L4tP+fXzaTc+ZWU6XFF21/aTSlnD+e+dn8rpcGgvl12mpQzB/Hcurbb8rl1r/H1Z5ZbPn5tt5Tv+0u58AEpv2jI1+t1RTtOBwAgRNp7ndu7YWnfANQ23lcHcBJAP6vn16kQUGeV/TwJIARASP369Z31HatDr+fJ9Hk9KVNjpMxIkPKbFlL+0k3K3EzLbaNPS/lxDSlnj5BSl2t7nwlXeZJ/09Ix8jk6jyfrd22kjDxWtM9jD2LOSLniGX6W9/2k/K2vlIf/kDIzqeTeszxAr5fy8hYp54/n5/4oiESQGu34vnKzpFz6GPez6lnb50NGAs+lT+tIGXPW9HhejpSzhkj5aW0pYy9avmbze9zv8qd4v/8Xx4/PQThCQOXWBTOHEOIDAOlSym/MHit/LpgS/BsznSriFU+zJuuJ7WynqSA7lelSXTYVq7ZiIKnRVEvnpDEtb4/5rMsFNr5Bv7/xAOCeP52viJWSLUH2/QRc3cYWp+0msBjTnsZotxriLrKg+OQ/lAJ0eZytbR0JpBsMLNPZ8w0bjd37p3q3w5RIJjPcPIAndpiCzymRwO99WRg7ZZvJLczLpiumz2EmMSIEeGa/cwLqNuCIC1Y+GohYQQjhI4SorPwNYBiAM1abrQbwiCB6AEgpiHxKHOlxTG03HsCU6bVdbCje63lL8gFIUslhFJrZIp/sVKqhMxOAh5bZRz7ZKdSIhMxhe4aHljuXfKRkLGRmfwbDY06zhehLZ5iluR3JB+Ac+TG/AM+GMI5zaAbwcydg/y9cEOyBiwsLUkd+zbqvRRM5JdUa/nVYN5h2k72hlP3716FCOuYU54kpcPdi+5bkMJKVixuw8e3if2ZnwV5TqTRvABqDbtdJAGcBvGN8/CkATxn/FgB+BXAVwGkUEv+RJR0DWvWslB9WpQmcmyXljx2l/KF9ftdLifvs+Nz2vvQ6Kf++V8oPqkh5ZZt9758cIeWvPXkMxxcW+WPYxNUdUs4azGP/oZ2UIX/yc2rIj7hLUs6/m9/Vz13s/w0VHJnD2M3cu6TMyVDfRjmPVj9v+fiqZ6X8IEDK0P2Wj695kY8vfoSvu7TZsWNyALgVYkAlcSsxAoo8xhNm49v8f9sn/JGvbLfcLjFUys/qSvnHUCl1ebb3t/Ftvv7wH/a9f8xZxog+reP4yV7ovs9I+dedPJ5vW5F4CopZaTDh4kYpf+zA727lNAbq7cXxhTyn/r7X9vetxHdOLjY9lp3KBeL7NgxcK8hKkfLr5lLO6M3F8ceOjB2VABwhoHLpglUoSAlseJPmbf/XKRDb9wPbdpqn3KVktbqU7NFjS3l8bB5Hr3R/it32CkPUCZZvQFLE6Iw+zwDTuOtfY/wg5hSrsJ87Zl+LCQ1E8+HA0wcYDzqxAJjei3Vc9qDDA2wwdnkTzxuDIf82g96lUHLtS0D8FT7mWZm9ppLDLF0xLz9OWo05TVc58SrlH2UMjYCKi/NrWE4w+D0qWrd+SAXt0I8stzs+n8K24Z/YroqOOsEuhk0GcaRyYYg4ygJED18GqQtSUNsLKdmA7OfOnObQ5XEST4+nK05j+fIEdy9e+I9vYenI/HGc82VQqVGzRtfHWSpzajEnmVjD1Y1aM1cP6onyjMryBj2NBakzWO6ioN19rEUL3Qs06EO9WmaiEz5k0aERUHFg0PNkCmxB2bvSwLz386y6VpAWw4rzBn2AjjbEhtkpVK/6BAHj/yi8NivsEAPBXgHAo+ud0zoj8ToJbc3zHNw3dQ9wxzdl1lfmlkLdLsDU3Wxju/trtthNjy38dX1fpTV8cDqTC9bwrwOM+w24eRrY8Ynp8SEf8Hdb86KJ7FxcgBFfAukxtJRyUoF9PzrhwxUdGgEVB6eXAvEXgYFvs/fyprcB35rMfJljw+tcne78Ub0eR0pg1bNsUn/PHPUBduaIPsVsl28QLZ/iKosNBvZ7ntELiDzO45y8zrGpDxoKh7s3s2VjfuViNXMAK/8LghDA8M/YqWD9a7RerNF8OND5UeCAcb8A28OO+AKIOgYc/cu0bb2unDd3bSfV3Id+5wJZRtAIqKjQ5bITXc12nAhxcQMruge9Y1lBfWUbcG4Ve8HYasEQMocy/yEfFN7JUKmE96zMotbijr9Ji+FqvPENVoFPO0jzvaTrxG5ndHyILpk0ALOH8xwpCC6udLWqNAIWP2yqqTPH0I+oAVr5jCl93+Zu/qY7PqOsQ8GAN9k/yc0LMORxpn0ZQSOgouLUIvaqGfQuL9adn/MEaW/WNkGvo1VUpVF+q0hB4nW24mw8EOj1XMHvmRpNspAG1pUF1CveZ7i8BZjRm6vmnT+xGLY4RZYa7EetdhQMVmnAETwnFxW8vXcA8OBizhRb9jiLmM3h5Ue9T8JltmcFeF4O+4RtSPZ+Z9q2WhOjVm0H441H/zL2XCp9lH6/xlsBBgNVwDXbsbjy4gZmisbOsIzdHP0TiLsA3LdAXdVqMHDFcnGlaV6Q1ZGbASy8l0FDRyrh1aDXsVh2/0+czX7PnFIpUnQa9HnMNiZc4S01ir1x0mN5y02nFZCXRbW5iys7CLi6c9WvVJUdEn2qAZVrUxVcrTFQ1dh/qLQG/PnXoQu9eCLbpBh0ltXu1qjWhILPZY/RqhnyvuXzTQYBHR6iALLDRPYJqtOJwecD0zkIUXHX+79O0tPnsfXIwekkq1KGRkBFwaUNXGnuns3/Feun7QTTNllJlNY36me7YdihGUDYfhJXQZaHwcDOijfPsgdQnU5FP/asJE5ZuLaDJ+TwzxibKK/Iy2bqOOoYOwJEHSfpmHc69PDlUETfGlQle/rxM7l7Gd0MPV0NvQ7Iy+S45YwEtsW9vJV9sBV4+rMpWZ1OTHE36O3cpmTW8PLjb7roQWO6XWfZ98cabe4Gru7g8IDG/am8N8eQD+jOb3iDc+eEYIb23CrOqR9rTL0H1Gf3hWPz+DlD/mLA2zugRD6mLWgEVBTs+4k/YPBYzoWKOcX6L3PrZ//P7Fw3/DN1yyY5nBm05iPYBL4g7PiU7TSGf0aLq6iIu8i2Dsnhhbf/KCsY9JQjXNvBerOwQ6xjAkgwtTsBre5kw/7AZrQKvKsU/f2kpPWkjLSJOkGy2/8zL3IXd8blmg3l7LaSaCXi7s0xzIsfYtbK0y//NFRzjPySjdNWTmPMzrOy6TnfIHbg3PQWcHE9Fz//ugxSH55pOWyg17O00t28gNw0tgru+4rzP18BqBDFqM6CU4pRww8Ds4eyD3P3qcC8sVxJXzxtKgDMSAB+aMtevPeopE4B1vpc2QY8e7jgk/rsSqbnbU1NtRfXdvE93b2A+/4u2tiekoIuh1mZc6tZB5WVxMdrtOUqX78nLRJzaUNJIy+b+q6r24Er25nmBoB63SkybXef5fwup7xnFs+nqGOs4zPv92yNsEMsVO72BDDqa8vn9HnsxqnLBp49QtczNRr4sT2Lhsf8Ytp2ySP87oNaMh750tmi9zcyosIXo5ZrHJkNeFSmjx17nit1tycsf7R9PwC6LKD/m+r7uLyFFk3/1womn8RrbCZWpwtwx3dFJ5+zK5i2D6jHyRzlgXwMBpLi8qnAV03YhP38araKvXs28OoVdpMc/ikbcJUm+QAk6sYDmF16ei/wwkm6MjnpLCb+rhXT4ooC2Snv6Q088A/d+UUP0mK1hfrd2cv78CxaQ+Zwdeec+KTrHFwIsAld58ms2E+6Ydq21/PUoHn4ABmxPC9LEZoF5AgyE4FvWzJQOPo7ksOpJcBL50zanfRY4Id2nHgxfmb+fehygF+7syr56f22VxtdDhuPJ12nILCowwEPz+KFUr8HT+7iuCvOQEoEcHwBcOJvppM9/fldBY9hx8Zirr6lgsijwKGZFJ0adLSI+r9hf5fKwpAcDswaSJHpE9ttW1o56cD0HnTBpu6xDAFIyfMnNZJKdncvBut/bE8iMrea5oxgOw+A59nk4pGQZgGVFE7+w3hEl0dJRqeW0BQ3Fw4e+p2mb7/X1fcR8idJZcQXBV9sWz/gyJ2xM4pOPvt+5GrdYiRHuJQl+UQeY/r4h3bAzs+Yebp7NvDqRboEzYZWDPIBWM4w/ne6Kz2fpev4S1fWZDmjtCGgHvs4JV4DVj1DMlGDpy8txNhzwPF5ls8JwfYeqZEmBbVfbQaxjy9gfFJB1ylAShhQtSHLheIuFf8z2AmNgOyFlCSPul1Zc3V6GYmm6xTTNrkZDOS1vENddJidCuz+ipmxgprPX9/NtGi3qY6N3DHH/l84mbPN3ewjXBaZLimBS5uBOSO5ol/axJqyF05y1n3be8p3Bq4w+FbneJ0XTrJu6+hc4JcudHuK61k06kv37/waDje0hVZ3MVO3/VO6Uhb76Mfbvh9oUQMs68jLMLlmAKdzeFehW+zizsB0KUEjIHsReZSp906T+P/x+dQB1Wpn2ubEQgZQbQkK9//MBmNDPih4HtSqadSkDPmgaMd6YDqLF1uP4xC80p75JSVTxbOHUruUEs4M3svnuGIXNqK4oqFyDbo0U3cxO7dqGiedpMcVb789pzFLuvUD21aJEPxOMxNYXGqNPi8xy3d6Kf+v3QGo34tDEpQaMXcvCmjDDzHIfmoJJQulAI2A7MXppRxcF2wcfxtzyjKNbTDQaqnTmT+iNTISWKvTehy3sYWt7zMGMObXwgfvqeH430zBtroLGD+r9Mkn8ijw12gqtlOjgNHfMwbRc5rzs0blDTXbAo9uZIb02k5gRk/722+oQQgq1N0rcaCltfpZQe2OtCYPz+R0WXM0Hshs4v6fTS09ejzF+NuVrabtOk+iVkrqqZy+trPox+0ANAKyBwY9M0nNh7HlxvEFJCPzMSrXdtBn7/GMunVz6DeavrYyYwCzGUf+oJvSoKfjx3l5K7D6eZ50d88u3b496bFc+WcNZoHuyK9IPF0eqzixHWfAxYXyjCd3Um39990MWBcVlWsw4RF1vGBXrN9rTOMf+NXycSFokcdd4DBEAGgxisdm7oYFteDCmJXMc1yxmEoYGgHZg9A9NGPb3EMyOreSFcjmbSqOzeO44FZ35n99ThpN3pajbZc8GPQMGPvVYQ8YRxF1gpqOGsHAffNL76I36FlJ/3Nnjp7u9SyJp/vU27t/UI1gYMoWulAbXmMmUq2pmD0IHgs0G0Yls63K9aAWtK4Pz8wfCG8zHqhcixISgAtTuwksITLftvV4IO48LaYLazmZtYRRLglICFFPCLFDCHFeCHFWCPGCyjYDhBApQogTxtt7JXZAZ1dQ7t98OP3k9JtA67Gm5zMSOIWy/f3qNV8hfzJA2Odl2+9x9C+WHAz7hJoMR5AWw5hDpWrAxGWWytiSRNwljpTe+AbH/z5zkMd/q7ta9sKzMkWfPZ8lMax5vmgkJASzpvpcFi7bQr/XWAcXMtvycVd36taubDGl2ztMpMt1eplpu9bjeO/iyv1c3uz4sTqIcklAAHQAXpFStgLQA8A0IUSwynZ7pJQdjLePVJ4vPqQELm5kPxZ3b9bUuHlxRVJwahF/zI4P53+9Po+d6Rr1A+raiP1kJgLbPwYa9jWdBPZClwssmcTmUg8uBirXdOz1RYFBD+z9ge1a4y4y0D1xme12I7czXFxJyv3fYOJizXNFy5BVa0JX6vQSWrtqqBFM8WTIn/mDyB0fYheFEwv4f802TKKc/Me0jX8dqs7ToqlBurTR8eN0EOWSgKSU0VLKY8a/0wCcB1DMxjdFRPRJdpBrPoKr17lVJCNzK+PkP/Sfa6hw5MUNQFoU0P1p2++x93v63iO/dFztvPkdlgyM+UX9/Z0NpX/Q1vcZE5t2GGh/n9Y/qCAIwaZ1/V5n3MW8V7Mj6P0C4zMFvb7bVGp/Lq6zfLxqI5LTsfkmK6z1OJZ9JIebtms9jjVxgc1oAdnTOrYYKJcEZA4hREMAHQEcUnm6pxDipBBigxCidYkcwKVNAASFctEnuDqYx3kSrtJ1sjXXO2Q24FeX7psa0mKoVm53H1DDwY9wehlN+57P2n5/Z+Lqdlo94UeYpZsw3/ZcMw35MfBtWiK7vgRO/FP49tbw8qcVdGkj+4GroflwwL8+zylrdHiIgsMIY9fEVnfx/oIZWSnFzsKVqf2II44fpwMo1wQkhPAF8C+AF6WUqVZPHwPQQErZHsDPAFba2MeTQogQIURIXFwRdBmXNlJ86BPIKaAA0MRMRHh2Be+Dx+R/bcJVpjM7T6Yprobd39B9G1BAdkwNyeHA2peBut2AIR869lpHISWPc/54Zk+e3MELSbN6HIMQwOgf6GqvfQmIveD4Pro/RdGgeYMxc7i4Uh4SusfSsgFYHO3qyQJngC5zUCuKHRVUbUwtkzSwXOjiBseP0QGUWwISQriD5LNASrnc+nkpZaqUMt3493oA7kKIQJXtZkopu0gpuwQFOTAqF6CoMOo4XS6AVdE121mO3D23kgSl1s/n2DyuJLbaXqREMPjc8WHHmsobDMDKp6nZGP97yWp98rLZi2j7x7SyntjOhvUaigZXd+DuP5hoWDpZffppQfCszAXt4nr11qwANUEAcGaZ5eOelXkun1tlcsNa3cmeVBkJpu2aDWPooWZbtkQpQZRLAhJCCACzAZyXUqpSvRCipnE7CCG6gZ8lQW3bIiPsIAAJNOzNMoqIw5YlFEmhdL/UrB+DATiznNvbclMOzuBK07eA7Jjq66ZzhRv5ZYnO+EZGAjB3NAOfg/5nvHCKII7UYInKNblwxJ0vWj/mrlMACGrG1FC1ES3jUypanuAxjEkqrlWzYTwHzYmm6RBjDyZBIrIu8XAiyiUBAegN4GEAg8zS7KOEEE8JIZ4ybnMPgDNCiJMAfgJwv3R2aX/oXpqsdboAN/ax8tl88J+iJG0+Iv9rIw7T325zj/q+s1NYO+Rok6vkMDYoaz6SqdSSQmoUU+wxp4EJ85ji1Vwu56HpEKDd/SwYLqjthhr867JFybH5ttXR7SYAsWfzu3ktRtAqv7yJ/9fuyI4E5srn+j3pfgkXkpN1uw8nolwSkJRyr5RSSCnbmaXZ10spf5NS/mbc5hcpZWspZXspZQ8p5X6nH8iNfZzn5O7FH8HFje6WgivbGfCrppJ+Pr0McPMGWo5S3/exeexC1+tZ+49HSmC9scp+1NclRwiJ19jsKjWKbT3VLDwNxcfwT+mKbVIZOlgY2j/I1rK2Jmq0GMl7hWgUePnzHFZKRFzdWPh6bYdJHuBRCajVgQuuq4f6KCAnoVwSULlATjrNzwa9+H/EEaBWe1P1tj6PVetNB+UnAiVd33y4uijQYGD2qkEfrkD24sI69qMe8FbxJ2LYQuI1Vq/npAOTVgMN+5TM+2hgYqPPixQIhqkleQtAk0EMRtsqmfCvC9Row+Z3aq+NOmGK+zQeQMs66bppm/o92OajZlvghvPXdgUaAdnCzTM0P+t0IdlEHqNfrSDiCC0Y84yYgugT7C5nq5XG9V38wbs8av/x6HI5XbV6MGvFSgIpkZy2qs/hYMLiNL/XYB+6PclpuI5qg9w8WKJxcb3tQHazYUDYAcveP4AxjCBp9QCmRcacBBv0ovIagm64LVevmNAIyBaUmdq12vMH0GWx3EBB2AHeq1kIlzcDEOrkBFAR6xXA2jB7cWwuV6ihH5VMkWlGPAWGmUnsR1waokYNdMF6PE0ycDQW1PIOTvm4sU/9+aZD6EYp56qCOp3YVlh5PLA5/4806xaqLLZSzwUprgiSATugEZAtRJ/iylS5Ji0awNIiCD8CVGumPjf90iajdkhlxHJmInB+LYOE9hZr5qRRvNawr0kS4EzoctiwPjmM5RxlaPlIKRGemIkDVxOw6kQklh2NwIrjEdh5MRZXYtOQpy9iQWd5RsdH2AhMbfZ7QWjQm0kSW3GgOp24X+sgsosrUKejKROm/B9pJm70DWIBq4tR4hF13LFjsxPaWB5biD5J60cINp/38GXAGWCwLuKwevYrM5Hy9oE2AosX1nFF6fCg+vNqODwLyIgDHljs/MCzlBQ0hh/kBI+GvZ27fzuQkaPD1vM3seF0DLZfjEWuzjbJeLm7oGO9KhjYMgh3tq+NWv4VuKOiAt8g6nFOLQaGfWq/rsujEn+vK1sBqLhw7t5sQBauEl+q25UZuLwsblenM7to5mWbFsYabVjW4VGZMaMSGOOkWUBq0OtoctZow/9vnqP4TpmYmXiNMnXzjJgC5ce2Fbw9v5pp91od7DuWvGzqhZoMtl3MWhwcnMEG8f3fKJ1yDjPEp+fgiw0X0Pr9TXhh0QlsPBtTIPkAQHaeAQeuJeCz9RfQ64vtmDI3BCfCk0vngEsSrcdR+GrLnbKFhn1Yu2WrF3W97oxfKi1ZFdTpQvdMCTXU7kRFfuxZ0zY123LfQc0ZkC4BaASkhpQw/hiBzWghxJ5l8FeB8qOpuSphB5i6rK3yXHYKW5W2ust+S+bkQga0+7zo8McoFJFHgS3vMhZVUKM0J8NgkJh/8Aa6fLIVv+26WuT9SAlsPX8TY3/dh+f+OY7YtGwnHmUpo+lgdllwdCxOHePwiahjNp7vRIs7/nL+xwFTZX2QsU+VeevXmm1IUgY9p9GWADQCUkPCNd5XbULXJyvJkoBiz1OkFdgi/2vDDjK1rhbfubSZxKYUARYGg4GtNOt0ZvzHmchO5ZQK35qspC+leejJmbmY+MchvLvyjFP3u+5UFIZ8uwtbz9106n5LDR4+zDw5qrmp3RGAoJWjhiBj2UzsecvHfWtQExRvDHxXbcR4UbxZIDywOe+Fi/E6SHbs2OyAFgNSQ6JxVa7WxDTEzbxWK/YcycmaZHS5DNZ1fwqquLqdXRPVXDc1hO6muzd+lvNjP5veApJvAJPXl9q4nqjkLPT7agd0BtuCdU83F3SsHwBPN1foDAZk5eqRqzcgM1ePzBw9MnN1SM3O3zDdx9MN3h6umDIvBK8MbY5nBzWFqGjK7Xo9mI7PSrZ/RruXH8/TmNPqz1drykBynBUBCUGrR7F4XN25H/NMnDI8QCmkTrjq9DCARkBqSLjKwJtPEMWGgGW5ROx59dYZ8ZeonajVPv9zUlL/06iv/dbG0blM19trMdmL63vYl6b3i0XrPV0E3EzNRq8vbDdo71AvAPWrVkJ6jg5HQhORpkIyBSEtW4ecPAN6NK6Kb7dcQkauHm+MaFGxSKh+dwCS7pR5yU9hqNqEC5Ua3DxIQtYWEEALx7zaPdAq1uNZmd0PlExYwhWNgEoFqZFUGgthqjj2NyqPdbnU46gFbG8aA3hq5JRwlftt9Ip9x5CRwHhAl8ec21tZlwOsfREIaMDAcykgO0+PO35Sdy3qVvFGj8bVEBKaiBPhyQj09cDINjVRJ6ASUrPzkKc3IDUrDxm5eqRl5yElS4eUzFzEpuXks6Ry9QYcC0vG6Ha18Nuuq6jq444n+zlpWmlpQHGX4q84RkDVmrA4WUp1S7lKw/ytOQASzvH5JosroD41bOb7qdqIcSAASI1w4MPYB42A1JAWQx8ZIAFVqsYplAB/BGlQn1Yae5YBaLXasNA9vG/U375jOL+a1pSzC04PzuBK9tC/pVbZ/vHac4hPz8n3eO+m1ZCnk1h2NALNa/ji3dHBSMzIwZqT0QhL5MkuBFClkgdcBJCeo0N2HrNkLjYMm1ydAeGJmRjVtiY+33ABzWtUxoAW1UvsszkVvtUp93A04Fu1MQWJ6TfVW/L61VFPxfsbm4ymRpGA/Oty2GZmoknDFtCAAkVPf9sN8YsBjYDUkB5rIpH0mxRkKUgxrgL+KrVYsee5qqgplaOOMdZi7/zw82uAKo2YCnUWspLZ/rXp0JIRNKrgeFgSFhzK37emf/MgXIhJRVJGHl4Z2hxJmXn4auMF5OoN6N0kEE/0a4xuDauiUaAPPNxckJ2nh84gkZ5NF23T2RhsPBOjGk86GZGCO9vXxrW4DLy+7BS2vNQf/pVKcURRUSEErZUUFWulICiLZUacOgH512EiJTfTctHxMyOgGsFm/0eaCMi3BlXyfnU0AioVSMke0EoPn8wEWkAKFFNWrQFZUqgpnWmNqOPU/tgTk8hKZryo5zTnBp/3/wRkJwODS26AiDWenJ+/dWj7egE4GZEMD1cXvHtnMGbuvoqIpCyM71gXLwxuhtoBXjh8PRELD93A8fBkXLqZ9p/l4+HqgibVfdG9UVX88mAnrDsdjTUno/K9xyfrzmP+490w+c8j+HLTBXw2zolEXpLwruJ4tkk5P62HEiqoXJv3adGWC6Cf8fFU46QMcwJSJv76BHJChpdmAZUOslPo+vgYzfbMROohFPxnAVkRkJQkJ/NpGQrysmkd9XrevmO4spV+tyO1YoUhO5UD8lqPsxwnXYIICU1EXJql61XL3wsxKVnwdHPBo70b4eO151DDzxNLpvZEh3oBWHQkHDN3X0V4Yha83V3RsX4AHuzWANV8PeDqIpCUkYtz0an453AY/tofiv7NgzBtYBP8uuMqBAxoKcLRxuU66olY+Gz9FwtrSxw+BsTVGoag4AGAXy31gy0v8PIHEq8Xvp05FALKtNGPT8lyZidbPu5rtJYUYvE1nvPmROZj7P7p4moiKidCIyBr5BhbT3v5897aAsqM53PW878y4lmwqtZcLP4iCcXeCz90D+Dppy5mLCpOLGT1fk8bc+tLAL/uyB/L8HRzQVRKNp7u3wRfbbyAjvWr4I9HuiAiKQt3/rwXF2LS0LF+AN4c0QqDWlaHt4d6L+3U7DwsOhyGn7ddQXbYcSypvReNEnYhSLB7n14KpMb4wN/bHV1ckuG6cTmwUQBNBrJMpm6XEv3sRYabF88jR6C0fMlNL/j5nDSr9/JgzEnpeKjMczPvgKgQkC4n/+udAI2ArJGbwXsPH1o12cmWOpmsJHXdTIpVtswcSoq0WjP7juH6HorSnNXr2WDgWN963UumnEMFeXoDdly0HAIQUMkdoQmZuK9LPczacw3Btf0w77Fu2HLuJl7/9xSqVvLA7w93xrDgGoWmz/283PFkkxQ8cvEHeEXuR2aiJ056d8dnqa1xUjZBqKwJA1zw5wNdsf5EKMLOh2BB30S4HZ8L/DEY6PsqiaiUBJh2Q+rZsdARKGlyWyN0bBGQ8pyy6HpUBiBM/wNmeiTJ19vKtBURGgFZI8e4inhWpismDVyVFGQlU5tjDaW5k69KxkUhIHsaz6fFUAjpSK+gwhBxmNIBRydvFAPnoqyHmACZOXrUq+qNizfT4OYi8McjXbHpbAxeXnISPRpXxYyJnVHFx3Kk9PGwJKw4HolD1xJxPT4DeQYDGvm54DPfReiesApePoFI7vM+Ru1tiCqVgnA22fJ9lx+PxAPdGuPBE3HYXGMcRj3/ArDxTWCPcRrJ0JKZZ1lkGPS2J6jYgrK9LQJSJu0qi6s5PCubLCcXF1re5haQ0oDPxY1WvC7b9JgTUM7o3wQhxAghxEUhxBUhRL4rRxA/GZ8/JYRwjr+i/BgePqYCPnN3KztZXaWq+Ndq5JRwjf62PSOXldqcOk50Ec78SxK11SCtBBByIynfY7l6AxpW88GJ8GS8c0crRKdk4bVlp9CrSTXMfaybBfmExmfgkTmHMW76fiwJCUetAC9M7t0Qb3f3wHz5FnokrMRcw0gcvGMLAoa8jEcGtsfZqFS0q+tv8Z5rTkaha8Oq8HBzwfGwJF5wd/0CdHmc1eDXSnbqg8PISaVb5AiE8TKWNgp5C7JYPHwZH1TgWdnyf3dj1kyxsnJsuHlFRLm0gIQQrgB+BTAUQASAI0KI1VJK85LckQCaGW/dAcww3hcPOmNBo5uXsSMc2HNFQU66Ke1pjizjBafqnoXb30L1plFS7+iQQluQEji3msHx0poZD+DI9fzV2S6Ciuj6VSthTIc6GPXTHtT088JvD3eGp5tp1d95MRbPLjwOFwG8M6oVHuheH76ebiwC/vtpwFWPxHGLsHB7JSQuv4rdTevhns518cWGC8jT50/Lu7kItKrlh9ORxpVdCPZjvrwF2PkF0NhObVZpID1OXWNWEPIyeW/LMiloVoOrh0loCNDtN//fep95GQAcHG9VAAq1gIQQzwohSqdYyIRuAK5IKa9JKXMBLAJg3Rl9DIB5kjgIIEAIUfwUh5LlEsLMAjJzC5RG3dZQCMjLP/9zmYmUtNuDmDPUgigBweIi7iJlBWrZuRKEWouMWv7euBybjvGd6mDtqWhci8vAB3e1hp+XSaNzKiIZU+cfRYNqlbDhxX54ol9j6NLikRKyDPjrTi4Gj21C1fYj8e7oYMSn52D/lQQE+nqipp8X1AajpOXoEFzLD5dumq3e7t5A50mciZWaP41fZsiINQV+7YV53LJA2LKEzL4zZRKGAjcjASkLs5MHz9jjgtUELZAlRreoNIpr6gAwV2NFIP9seHu2cRzrXzX9rfwQwuxrMuSZzFFz6HJYTawWOM6MV++OqIaEK7a1REWBUsvWyMnV9IUgIzd/LVd6jg5SAn2bBeHfoxFoEuSDIa0sY2YfrD4Lf293zH+8O+oEeOPo4b0I+KUF/Nc+DukTCDy2kf1pADSrTosu1pjqr+ThCg+3/Kd0Tp4BNfw8kZiRa9lRsVE/3itub1kjO5ViQkctoFyzuKUalNiQWmxJCEtSEa4MhP8HaXrc/H8noVACklL+D3RzZgOYDOCyEOIzIURJFtmokZz1J7dnm+KNZlYsHfOG3Aa9utLZoFMnJimNqXw7LaCUCHWRY1ERGUIlt1LZXEpwU6mVSM3m99iwWiWE3EjEoJbVLbJdKVl5OBaWjId6NEBVHw9k5upYt2aEeGyjhSt7LY4XXk1/TxgMEjGp2fB2z3+RuboIDG9dEz/e38FyAVeEd+nlpIWHUomu1ITZCyUBYqurgZLVUiUoq9/JxdXSAlL+Vs5tJ1tAdsWApJRSCBEDIAaADkAVAMuEEFuklK879YiICADmQZO6AKztZHu2gZRyJoCZANClSxf7vz0pTURjTUBqRqBBr05ABh1jSZ52BBZz0hnMdiYBxV2w7GVUSgiu7Yd9VyyFccq5m5WnR55eol5Vy1o0b3dXVPZ0w8rjkQhLzMTuS3GonPckPuiQhr5jpuRzbxcdCYePhyu6NqyKc9GpyMzVq+qG/LzcUNXHD61q2XJrnXtRFRlKy4zqDlrASpGonw0HQMlqeap8fn0u4GZGTAadpQzgP+upZAjInhjQ80KIowC+ArAPQFsp5dMAOgMoqR6eRwA0E0I0EkJ4ALgfwGqrbVYDeMSYDesBIEVKGe20I9DnmQjIYEZAru7qI0oMOnVNiRLIdrGjFklRpFZ2klpXSnbCC1JpnFbC6N7ItsuZnceTOj3H0k3zcHPBTw92hIebC3ZfikO7uv74Yuo96HvfK/nIZ8/lOKw+GYVHejVEZS93/HssAq4uQtUCcnO1cZpblyCUNSKO8HMGNHTsdSmRJA21OjDApP9RiyvmplsujnlZlvViigWkBKadPJHFHgsoEMB4KeUN8wellAYhhBNrBSz2rRNCPAtgEwBXAHOklGeVsczG6ajrAYwCcAVAJgAnCmfAH0Zxwcz76bp7q89hcnVX12EoZKUWuM73nsYTRW2lKgpy05khsXViliAGtAjCd1suqT53JTYDDapVwv4rCXhmgGXngIEtqmNgIdXrZyJTMG3BMTSr7osXBjdDdEoWFh4Kw4g2NXHgqqXV5edVwCn+XzvS0idoVdzYD9Tv5bg4MiWCv7Et/ZBSoqEmEclJNwoQjcjNMKXeAeOMeJgWYXfndlAolICklDYrF6WUKl2OnAMp5XqQZMwf+83sbwlgmtPf2LuKsXI4g/ofNy9LZaiblykjYA53b17s1kpRZeVQc8+sYXc2w04ojcpLqeOhOYJtujvA6pORuLdzXXyz+RL2XI5D32b2Z33WnorCG8tOIaCSB/56rBs8XF3w6tKTcBECdQK8kZiRa7H9MwNVWqMouLqdsbFSjo+pIu0mExCdJjn+2rjzpvapakgtwELKTctvAZmTjGI9Keevk1u4lFshYplBMccVMlAISYEtC8jdm+aq3vICMMWRcvO/xhq5ip7DST+yQpROXrXsgZurC+pWUdelbDgTg26NqqFZdV88s+AYtl8oPAh8PjoVj/91BM8uPI7mNStj2dM9UdvfC++tPoN9VxIwbWATLD6Sv43F3Z1sxNPSYzkS2ZkFv8XBpY28bzzAsdcZ9AxeF6QbS4lk5bu1hWTQk2AU91afR4vHfAFURIkKAbk5TwUNlFMhYpnCvy7HMivMb01AXv7quhHlIlcsJ+vH7Skw/C/d76RA3381Qo61N3UWfrivA+757UC+x6UE3vz3FH56oCNeXXoSj/0Vgj5NA3FX+9poXccPVX08kKeTiEjKxPHwZGw7fxPHwpLh6+mGt0a2xKO9G0FC4vVlp7D0aASm9GmE/VcTkJJlGZtzEUBQZc987w8AOPoXv5fOk53/wYuC82toiTna/ynxGheaAgnIRmY1I46LpmIZZRizxObF14r1L1xIPk6undMIyBpVjeoC5cfwCrDsz+ITyKmp1lDcnKwky2mpru4kAlvzu83hiLVkD5QaNnveuwTQpaHK1FgjrsVn4NWlJ/HzAx2x/UIs5h24gdf/VfleAbSq5Yd3RrXCvV3qIqCSB05HpODN5adwNioV0wY2QWRSFvZbxX4AYPkzNoYsZiZy2kjzkRy9VNbITgGu7QR6POV4oacysVSZYaeGpOvqc+rSjDkbJemRHst783pGxQLS5fDcdzI0ArKGIhhMNsbcfQIt52L7BBlXDqtYj/Kjpcfm73ro7mNfDc1/uiMnEZBvdQDCdKKVAVY80wvjpu9Xfe5CTBrGTd+Ppwc0wYYX+yImJRuXb6YjNTsPrkKgThVvtKhZGYG+npBS4kR4Mv7YewbrT0cj0NcTX9/TDhvPxGDbhdh8+65SyR0d6gWoH9TOL2jhlmJjtgJxehmDvK3HOf7aG/uZtLBlAWUlMwZUXUVblKoQkJUFZF5qZK6RclShbQc0ArKGkilQpgj412ODMIVwfIJ4sli36VB+tHSVrnG+1e0Tuym+uLPmL7m6c3VLulH4tiWEjvWroHfTavk0QQDFih5uLvh600X8tO0y+jUPQsf6AWhQ1Qc+nq5Iy87DulPROBeVir1X4hGZnIXKnm6Y2q8JmlX3xbebLyIqRX0Y4eaXbNR33dgPHJ4JdHuCbUjLA47+BdRoW7T+Tzf2s82KrQyYsniqacEUGYLSMVGxgMyJJi2aWbK8LFMPaSdCIyBrKKtBjNEdCKjH7JbSqFshmrQYdQJKUyGayjXta2f5H4nlX9GLjJptgegTzttfEfDXo93Q7J0N+R7XGSQSM3LRKNAHVSq540JMKraoDBasUskdnepXwbSBTSEE8O/RiAInqv7+cGf12E9WErDyaZY6DPmgOB/JeYg8xnNt1DeOu18Z8Wx21/4+29sok1rULKCEK7TOlXM+JRyAsNShpUaxi2RGHFCno2PHZwc0ArKGtShNaTCWEkYCUlK2idctf1TvqgzSKWN8zFG5JhB+uPD39q7CeJGaFVVU1O0KXN5kLIi1HZMpSbi7uuD4u0PR8eMtqs9fj8/AdbBEY0TrmvD3doe3hyt8jcMGs/P0OBuVik/WnUNmro2eN0a8MrQ5hrdWSTcb9MC/U5gRmrzOeVKH4uLAL2yJ0fZex197eTPvGw+0vU30CVrWao3y4i8xBqYQX+I1BqvNx0ClRrE9ccIVUwtXJ0IjIGuYE1B2qqn2KDmMY3CrNub/SVZ9e11cOElDbaSKfz3g7ApLdbUaXFz4/mokVlQ0GQjs+IQna/v7nbdfB1HFx6NAEgKA0IRMhCZkFvk9/ndHK0zp2zj/E1ICm/9HV3r0D8YBgOUACVd5XvR6zv5JqOa4uJ7uU+0CLJOwQ0DdburWVfwVy+8i8Vr+pnmpkRzNIw32NdRzEJoOyBrmkf7Ea6asmDLC1ruKsXG4yiTKak2AhMv5Hw9szpSvPc3Gg1oCsRcK385e1OkM+NUFzix33j6LiCo+Hrjw8QjUCXCulgQA5j/eTZ18AAadD07nyGxndposLvb9wBKdHkXQ0+ZlA1e2Ay1G2nbdMhPpoqkRbm4GrXrzNsGJ10wLLMBAfVq0Kb5URSOgkof5j5l4lSrRgAamkbVC8EdKUIlBBDZjwFdnlcUyto9A/MX8r7FGUAuSmN5J2h0hgHYTKLpzdNpCCcDL3RX73hyEnx5wTjyhQ70AnHhvqLqaWkpg99fAri+ADg8Bwz93yns6BfGXOSig8yTTCChHcGUrm4O1HGV7G8Xtr68yfluRkihjxDMTWbJR1SyDG29cTBU9mWYBlRLqGVcMhWSqt7JMxVcPpljRujI4qCV7qShVzQoUmbw9lk3NtkzDm8/oLi66PUkp/oFfnLfPYuKu9rVx+dOR+N8dDraeMKJdXX/seHUAVk7rjYBKKnV2BgN7P2//BGh3H3DXT+WrAf2W9xgz7FfEZhIn/2FsptEA29tc38UGbmrZtUjjvLY6xueiT/LeXAiphBNcXCmoVesEWkyUo1+kHKGmcXzOzTO8r96Kq4FSWFqzHbMC1pkt5ceMPGb5uGdlxoeUH70gKKvVDXXtTJHgVwvo+BDTvXF2WGGlBHdXF0zp2xihX9yBvW8MxGO9G8HDVuU6gC4NquDPyV1x4eMRWP1sHzQKtBFIzkkHlk3mJJCezwJjf3O80XtJ4voexm/6vgz4FkFbk5EAXNpEy7agySmXN7MRnVr9VuRRxiYV/VqMlUUEMEgtXClCNA9WOxFaEFoNisVyeSvva7Sh9if2HH8gZb5XzCnLQXdVGjFGFHk0f6yhbjdmowobaxJQjydG2H4qY52FQf9jHGjty8Ck1eXrggRQt0olvHdnMN67s5janPjLwOKHePEM+4QB3vIEXQ6w7hXAvz7Q4+mi7eOMUbjY/gHb2yRcpQXT7Un15yOPWgavo0/xvDPPlMacYVwz7iJjTSUAzQJSgyLPz8ugTL5uV/4fcYT3iuxdMVsVCMEfVZHHm6NeV/rYasFrazTqR2m+Wt+hosInkI3Yb+xl8PNWg5TAqSXAzIG0Th9aXv7IBwB2f8NY4OjvizbexmCgkLJ2J8uJvda4bMw2Nh2S/7nkcCr9zWND0SdMlj/A7zPqGDVBmfGWzzkRGgGpwbw+KOYMp536VAfCjQTk5cd4j5q2p14Pir+UVhgKlB87dE/h79/yDhKfPds6go4PAa3HMy5ycaNz912WSI+l1bPcqG6eupvyg/KGm2eBvd8xJtVMhRjswZWttGx6PFPwdudW8hy1LgsCGBsCTNNA0uO4T/NpsWnRVO8rspGCyK4Y0AhIDeZaoKhjtGzqdTNZQAAnl4YdzJ+tajwAgDQ1g1cQ1JLp8Mu2dTD/ockgBv3Ory3qJ1CHEMBdP9ONXDrZuXGmsoBBz7jWr935vQ79GHh0g3Nb2joLeVnA8idZ6lOcbNzBX2mVtB5re5vkcCDsAND2HvXnr+1iuYVSnhFmPA/MC1YVK17p0OCsMVFW0AhIDUIAjYyrww1jO4m6XZiWV8okGvRmMycleKegTmfWzlzbmX+fzYbyces0vTXcvYHmIzhQME+91qnI8PQFHlzKup754yquJXRjPzBzALDmBZL7U3uA3s+Xu9jWf9j8PyY1xs6wf0KKNaJO8Pzp9kTBgtYz//K+jQoBSePi2Ki/KRYZuo8LXq0Opu0iQhiAzkknUamNm3ICNAKyBSVAd20n/W6FkJRJmg2MrR5u7LN8nasbMw9Xt+dP0zcbyjap1q9RQ+dJLHg9b90K2wnwDQIe3cgLd9GDwJ7v+BkrAiJCgL/vAf4cSTf3njnAo+vLT1tVNZxbDRz5gxm55sWYz7bjM1pQXafY3kZK4PRSTtZV0+1En2CpT9PBpsdu7Gec03z+XeheWsoxp9R1RE5CuSMgIcTXQogLxnHLK4QQATa2CxVCnBZCnBBChDj9QJSUel4GNUC12rPe6+p2Pu5XiyrSK9vyv7bZUAb5rLU8jQey7kdZoQpCw37MqoXMKd7nsAXfIGDyWiB4DLDtQ+Dv8UBSaMm8V3EhJReCv+8G/hjMDM7g94BnjwBt7i6R9LDTEHseWPkMg8aD3y/6fiJCmEXt/XzB1khECC2tDg+qP39+DS2b5iP4f9pNTuM1nxuXk8bvuFI1LpgNehX9uAtBuSMgAFsAtJFStgNwCcBbBWw7UErZQUrpxEHqRtTpbPr7xj6a9o37W1o2zYdzpTCfpQ2wzadwAc6tsnzcoxLQ6k6uiIW5Vi4uNLXDDjDWVBLwrEwLYvT3jG/92h3Y+aWpG2RZIzsFODyLxzVvDOMSQz4AXjwF9H3F6f2JnY6MBGDhfXSp75tvaWE4AimZOKhUDeg2teBtQ2ZzkWs3Qf3582sY61HS7UpBa7Phpm3CDlJQqyig6/co2nHbgXJHQFLKzVJKJbJ7EJz3Vfrwr2uqfbli1AM1GUTzVekV1GIk9RiKVaTAtzqnG1gTEMCq55wU0w9fEDpP5km3+5sif4xCIQTQ5TFg2mES6s7PgB/aAru+cm5bEHuRl83vbfHDwNfNOKnW3Zuxk5fOAX1eKtUZ90WGLhdY8gjFqg/8U7zA+KVNwLUdQN9XC54vl5lIrVe7+9S/o7iL1Ee1utNs3xuZdDFXQF/fxeZ4eZns/lCCQf1yR0BWeAxA/kYyhASwWQhxVAhhQ21VTChpyqvbmcVoOpT/X1jH+3o96JNfUgnkBo+h6xZrVZbRqD+zGEf/Kvz9PXyAntNYx2VPO4/iwL8OMGEeMGU7S1F2fAp81wpYNJGEYG3lOQtSUht1eBaw8H7gq8a8cMMOkICnbAee3EmXwrxNRHmGQU9JwI29wJhfLNPbjkKXw5KSwOa0iAvCsblsKt/1cfXnTy+lZa404tflAFd3cOExd2Mvb2EMNOII0KwYMSs7UCZKaCHEVnDmvDXekVKuMm7zDjiFdYGN3fSWUkYJIaoD2CKEuCCl3G29kZGcngSA+vXrO3agjfqRKPS5dLWaDeXFeX4V0P81Bpybj6CsXpdj2Yy+9Thg01vA8b8pAFTg6kaLY8enVO0W1pO421Tg0EyehI9vLfl6prqdgQcXs/r/+DzgxD/AhbWs2m7Qk6RbpzNjYpVrOhZ/0etYgZ1wlSLOyGMcHa10iwxowJYhLe8gURdUZlBeISXHSZ9bSSW2LVfIXhycztYvD/1bcOYrLws4MJ1WulrK3KDnb9l4oEm9f20XY5zNzVTOCVe5cDYbRgtIWXRLCGXyC0spC1RhCSEmARgNYLBx/pfaPqKM97FCiBUAugHIR0BFHs0MmDJfAK2cZkOBVncBm98xtS5oew9wahFXjVZmI158g+iinfyHwUdz/7/TJLo4h2cBo74q+Bg8fRn3WPkU38dWcNHZCGrOC2jwB0D4IX7+q9uBPd+YpmW6ebO7oF9t9iX28uNj0sAYgi4XyDJWWWfEUZ9iPmW2WlNeEHW78MKp2rh8B5QLg5TAlneBY/PoLhVXiZ14jedJi1HqimZzHP8byIhlbEwN13dzhPOwj0yPnV5KC958FNBFo8Ph4sahBmrN7J2IcrfECCFGAHgDQH8ppWp3KiGEDwAXKWWa8e9hAD5S27ZY8AlkOjMyBLiwHhj5NRBsJKBzq4E+L/ICqhQInF5iSUAAieb8GuDSBrpkCirXANqMB04sAAa8WXinwnb3MY27+V2uSEUpYCwqXN2Ahr15w8fsIxN9kqrepFDe0qJJLjmpjOG4uDDT4urBGFalqtSYBI+lMrdqE84/L4OBiSUGKYGNbwGHZgBdn2DtXXFgMACrn+f3OOrrgrfV5wH7fmK9oSIPscaJBcyetbiD/+dmMJTQ7l7LxfHiBup+Ys+TfEo40F/uCAjALwA8QbcKAA5KKZ8SQtQG8IeUchSAGgBWGJ93A7BQSlkyirpWo0lAaVFA+EGmJOt05urR+wVeoG3GA0fnMmtjniJtMogBviOzLQkIAHq/CJxaDBz4FRj8bsHH4OJCBfPM/sC6l4AJ88vOUvDw4XdQgqnZCgeDnm7XsXkskRj+WfF/n6NzWIpz50+FB4FPLKBrO+pr9ffNiOeC2elhUxzt4ga6X+atYNPjqIpuMphxxz4vFe8z2IFyF4SWUjaVUtYzptc7SCmVefBRRvKBlPKalLK98dZaSvlpwXstBlqaZQxOL+N9hweptVCKUdvdz+Df6aWWr3VxZeDw+q78has1gmkRHPo9f92YGmoEAwPfoUV18p8ifxwNToYuhwFnxe1yBvkkXAW2vE/rutMjBW+bmwHs+JyxyebD1bcJ+ZPnp3ll/KklbOda32whOfMv3WfhQhfMPFtWQih3BFTuENgUCDSqbM+tZCC1zT30j4//zcfrdGK18JE5+dXPnR+lLmP/z/n33f8NlnPs+9G+Y+n1HNCgD1tqxJwu8kfS4CRkJFCfdOZfYMiHtGSLSz66HGDZoySAu34ufH8Hp1MaMvQj9W11ucCRWYwhKWrx5DDKQDo8YJnUOLWY44HizhtDCyU/xEAjIHvQZjzvMxOox/AO4OpwegljHkJQHh97Nr9o0DuA6eQzy/M3m68RTOvp4HT72qW6uAL3/sl9Lppon+WkoWQQfwWYPYSZvHvmMB7oDGx+l9by2BmmgQi2kBEP7P2RcR1bYsGzK5hlNO89dHQu783HUsdfZuF1YFOep8o5X8LQCMgemE+TUPQ7HR9izOfcSv7f9h7A05+rjTW6P0Wzds93+Z8b8j5T3JvtDFr6VqdeJzWKJFRGY5dva1xYD/wxiL//pDUsB3EGzq0CDv/OOFJBvZ4VbHkP0GXZnnFmMLD3U2ALxnUAWkTH5tFdCzCTpZxcxHM0L5uZzJZ3FPfT2AWNgOxBlYam7MLFDVS3NurPH/bgdLpdHj4kpXOr8ls6AfXYIfH4/PzN7P1qszXnhbX5FdW2UK8bMP53ivWWPe68BvYaCoZex9jMogd4Tjyx3XkjfmJOAyueZoLDnqGJN/Yz+NzrOdPQA2ucX816xH6vmdyzC2uYru9iJlbU5/HcrN+LerfWY0us+t0aGgHZC6X9pdQz9iMEW6ZGnyQRAEDPZwAI9XhP31dp6ez8Iv9zPZ+lJmbNC/bNkAe46o78Eri4Dlj1jEZCJY2UCMZ79v1AIeljm01DKouL9FiqwL38gfsXWgpa1aDPYxzQv77tpvYGAzVE1ZqZ3CkpeW5WbWxZDX9+Dd20SlUZk+w0yTmfyw5oBGQvWo9lnx+AkneDnvEb7yq0ggCmS9vfTxPXuo6qcg0S1umlppEoCty9gDG/Ukuz9QP7j6n7VOpNTi1m4LKwPkMaHIeUdE+m92Ix7Ljfje1UnVQWkpfNbo6ZCawZU8YkF4R9PzBQPOor2zqdC2sYk+z/hqlH0vVd/Ay9X7Dsm3RkNlXoqZG06kuw+NQaGgHZC8/KJhVychhFXB6VmOW6sM7kWvV5iaUbB37Nv4/eL3CV2fB6/mxZ/R4MFB6ZZeo5ZA/6vcbU7/nVdA1KqmbrdkRGPLDkYWDFVE5GeXqvc6fL6nXAsseoNB83A6jdofDXRJ+kFd16nO1G8fo8YNvHrB8zDybv/Z7jlc2b2ceeZ81a3S5swdF5UqlqzDQCcgTmOgrFzer+FBW/SsV6tSbsu3x4JnutmMO7Cv37sAO0Wqwx6F2azMufcKwSvec04M4fWVg4ZziHI2ooOgwGZop+6cJK9KEfsemZ+dTQ4kJKYM3zdKFHfk1CKQx52cDyqVTe36GS0FAQMofDLYd+bLJ0Io+xp1LPZyxdvP2/MOiclUwLv+NDxflUDkMjIEcQ2NRUnBdxmCn3yjUY0Du12GQFDXybVtAulXhPh4dY3rH5XWZRzOFRCZgwl1bMv1Po5tmLzpNZsJgaCcwalL8lrAb7EHOaJL7meSCoFRvcW7ssxYUyq/7EAmDAW0B3O5s5bP+YrteYX2xrdLKSgJ2fM0liLkzc8SkXwM5m46JSIlhf2HQwz5dOj5Ra8FmBRkCOoqfZHO99P/G+9wusVN7zLf+v1oQ/9NG51IuYw8WFkvnMeJKQNWq05vPXd/FEcgRNBrJ9RaWqwLyxjCc5c7TPrYy0GGDNi8Dv/VkEOvY3Wj3Viza51SYU8jnwCzsd9H/Dvtdd3MjXdHmMRdG2sOsrWjPmiuzre9jTqu8rLBhWsN84KdfVHYBkTLGUoRGQo2g8wDQn7OI6rpiKFXRykWmedv/X2Uhrq0obzjqdmD49NtfU7MwcHR/ibffXlMw7gsCm7J/T6RH6/LOHsXBUgzqyUxkv+akjU9Fdp7DVa4cHnB8LkRLY8IaRfJ5kFtOe90gKBVY8SbV9QRM1ok5wGmznSaYxOlKy5W7l2pa9pDMSeP41G8a2wsFj2NmglKERkKMQwnLV2mE8Ifq8xMkCW97j/77VqY69sFa9b/SAt9kUftVzXLGs3+OO71l2sWqaaTKHvfDw4Sz0e/9ib+rf+wFbP9REi+bISibB/9ieLUaajyDxjPqqZEoQlILVw79TdjHyK/vIJy8bWDKJ7fcmzLOdfdPrKOOoFGipI7qwjo3FBrxpOQhx/488H1zc2MXAVhuPEoZGQEVB0yGmqRkX1zF74BsE9HuFzcmULFav59l6Yv2r+XtAu3tRbp9+k6N6rbNibh7sIxxQn5MrrDsr2oPW44BpR9jOY+93wK/dgFNLK84EjJJAeixd0+/bsM9y3S7AEztY4uLMILM5cjPZ5fHoX7zQh31iH/lICax7mZMsxv2mPuVCweGZ3G7kl6Y2J3nZdPcCmwMdJpq2TY1iEXSzoUxcBI+xbMlaitAIqCgQgpXpCrYbi/G7P01x2OZ3uOK5eQJ3fMOYwv6f8u+nTidg4Fuc9X30z/zPV6oKTFzKLNu8MflV1PbApxowdjpLBjz9geVT2Nbjytb8pHerQkpOtV0+Ffi+NbD3B04mnbqH368yAaUkkBEPzLuLlsjIrzjNw17Xbu93DFT3f6Pg0ozE6yTTZsMts2n7f2Y3xZFfWXaX3PUlz0/Pypx60f/Non02J0AjoKKi6RBTx8Sr27iSuHsBQz9gXOjYPD7XZBBPit3fsM2pNfq8wjqdDW/Qh7dG1cbAI6sAgw6Ye1fRU+yN+jGjM24m3Y+/7yYRnfn31lVRZ6cwETCzPwtHL6yjyvfZI3RPa5XMvPP/EHcJmD2U58N98x0L8p5dAWz7iP16BhQwGEav48RVFzfgjm9N5JYcxqRI8BjLMdXxl4Fj89nn6uJG6oRqBBft8zkBGgEVFUIAwz42/b/xLZ4MrcczdrP1fZMOaMSXTLGvfCr/xe7iAoyfRd996ST1CvfqLYGHV1Im/+dIdSKzBy4uQPv7gOdC2OgqN4NCuJ878WRNjS7afssTdDkcaa1M1VjzPBXid3wLvHKeFmlhfbidgfNrgFkDGeSetMax3jphB4EVT7H/9l2/FGwx7f2OkpDR31lWz280ktYwq1ZZW95jKxldLnsEDXjb/uMqAWgEVBzUam9SlcadpxslBHDnDwzwbTKeBJVrUDgWeZQyemv4VKP+JzWKF45aSUWtdsCktUyr/zmCwrKiws2TmZJpR4D7FgD+9bjaft+aNUmnl1UsRXVWEmNbyx4HvmkGLJ7IYs3Ok4Ep24BnDjADVBrjfAx6BvwXP8T+O1N3s3jYXkSdABbcy7Ke+xcUXPIRcZSq6LYTLOfAn13B5Ef/1yxJ6dJmxihbjGSb4G5TmTUtQwgbPd9vSXTp0kWGhDh5iGpKJPBLV7a39K4CPHeMsZudX1DHM/FfxhsAYOmjXBmf2K5u/p9aQhV0h4coNlNb+RKuAvPH0lK6d65p38VFwlWmoU/8wwZXrh50H1uM4ngiZxVeOgN52WyTG7qPzdbDDrBI2CeIaeXW4ymXKO2pGqlRLNu4vpvkN/KrwgtLzRF7gRauhy/w2IaCW7FmJtK1lBJ4ai97RAGMOf3ana+dss30HehygOk92HLDJ4gzwp4/ViJ9uYUQR+0dFlruCEgI8QGAJwDEGR96W0q5XmW7EQB+BOAK9opWkR1bokQICKCga7MxKN1pElPguhzgtz7MgDyznwrTzERgek9Ounhyp/qKvOMzBgkHvQv0e1X9/VKjuUrGnmVGpcczztOsGPScQXZ+DW8pxtYiAfU5Lrp+d1p+Qa2KPunT0eOJv8w4SsxJrvqRIVSaQ1Dv0mwYR8vU6VzyY4ts4dwqNpHX55J4Oj3s2OsTrgJ/3cGWqI9uoJjVFgx6YOEEEt2jGzlKSYGyyE3dbRnb2f0NldRdHuf01FHfFD5nrIi4FQgoXUppcxyoEMIVHNs8FEAEgCMAHpBSnrP1GqAECUivA2YO4IxtgK5So76c0z17GM3j8TP5XOheYO6dDEzfPTs/cUjJoOLpJQWfJDnpXG0vrKXFdMe3zh/cJyVnRF3fQ2V26F4gO5nPubhTJRzYjNZRlUa89wniNBHvKvaVL+hyqQrPiOPqnRLOrE7Sdd7HX2bTLQBw9SThNOjF/kz1e5T9ZI3sFGDT22zRUrsjMP4Px92am+eY5ZQGxosKCwori9To76mMVnB2BbB0MjDwf3S/FCRcBWb0pisYf5nf2dTdJWYhOkJA5XEqhj3oBuCKlPIaAAghFgEYA6BAAioxuLrxZJhtdIfWvAA8vY8ak36vsSas+QhmHBr2YQuNbR/xQjJXpwIkpLHTOZ99/au0ktQqsD19OR1j1xc8GaNPsDWo0vfXGRCCJFO9FeuVDAYSQ/RJ3mJOMa51diVdIMsX0+pz86I75+rOmz6X1qEuB9BlMw1sDRd3qnKrNOL3VbMdXdbA5gUP5yttXFhHDVdaDPU9A95y/PiiTgDzx9FVe2Rt4b/fhXX8vTtMtKzrSgoFVr/AOkPz9rAGA7D6Of4GnpU5Qum++eVm6GN5tYAmA0gFEALgFSllktU29wAYIaWcYvz/YQDdpZTPquzPfDJq5xs3SrBSfN2rppasfV6iIlWfB8wZASRcAZ7ezxHIBoPRhN5Fa0mtq15eNrDwXlod98wpuFr60iZg5dMMfI/8Euj4cOmO7dHnsbAx+QatmEzjMMKsRBKNPo/EY8gjubh58YJz8+Rq7BPILKBPIGMXfnWcW/zpbKTFAOtfYwuU6q3ZPN7cDbIXN/abGpFNWlW4EDLyGN20oBZ00xRls3KOxV8GntptGa87NBPY8BpjUkfnUgow8kvHj9UBlHsXrKDRzAAOAogHxecfA6glpXzMfCMhxL0AhlsRUDcpZYGjKEvMBVOQm0FTN8nYYP6xTXQTEq4Cv/Wl+zB5HVfJzETgj8HMNj2xXb0OJyedep2Iw2xYVtBU1NRoBrBD9zAecse3JDsNzoMuBzg4g5IFXQ4w4A2q3YtilZ1exkUjoAHwyMrCZ38lhwGzBpO8p2xlZlXBlvc4WeXevywXqqRQNlKr04kLQ04qMO1QiWcDHSGgMonYSSmHSCnbqNxWSSlvSin1UkoDgFmgu2WNCADmIwPqAogqjWMvEB4+lMwrWP4ECaZaE2DMz2w8pVTAV6oKPLiEVsHC+9TT3p6+wMPLKSJc+TRXM1vwq0XB4rBP2Fpheg/2hbmdyy6cBSkZZP61G/VdDXoxtd/3FcfJR0rWoP37OCeZPr65cPLJSgYWTCDpTVxqST7nVpN8Ok+2JB99HvDvE7QkK1WjTGT0D6UjRXAA5U4HJISoZfbvOABnVDY7AqCZEKKREMIDwP0AVpfG8RWK+j24KgJctdYbg4Ft7mapxqEZVB8DDODeOxeIv8SWqmqtMzx8gAcWc/TKhtdY/GrLanVxZZX9M/uZqVr7EnvbRJSg1XcrQ0rW9c0ZwVou90rAwyuABxcXnKWyhbwsYOUzLJtodx8Xl8IKX3PSmfFMuMLYTfWWpudunqVgsU4Xil3NsfNzWs5txnNyS5fHgebDHD/mkoaUslzdAMwHcBrAKZBUahkfrw1gvdl2o8BM2FUA79iz786dO8tSQW6WlDN6S/m+H28nl/BxXa6UfwyT8pNaUkafNm0f8he3W/qYlHq9+j51uVIuf4rbLZvC9ygIBoOUx+ZL+VVT02uSw53z+W4HXN8j5ZyR/O6+aSHl4T+k1OUVfX8J16Sc0Yf72/E5f5/CkJsp5Z93SPlBFSnPrrJ8Lj1eyu/bSvl1cylToiyfu7pDyvf9pVw0UcpvW0n5U2cpczKKfuwOAkCItPN6L3dB6JJEiceAzKGkPnVZzEBM3c1sUmo0OxYKQaGYn9Hg2/Md+7Z0ncL0u1oQWUrGH7Z/TJn+/QsYuC0IOWnc94Ffuc/OkzmX3q9Wwa+7HWHQc+zSgV85I923JkcmdZpUPInDpU10xwGW3dgaoWwOXS4V3Ze3sBF++/vMnsthbDD8MIPR5gHw9Fjqz7z8GV+6tgN4fEvJFtxaodzHgG4LVGsCjDU2ptfnUpqfncIL/8HF9Ov/uY+Ba4Aneu8XgCN/MEWvtjAIQXHiPX8y7T5zQOHulWdlDj989ghHSh+exR4461/jFA4NdHMOzQR+7syLPiUCGPEF8MIJZo2KSj66HLbDWDiBQs4nd9lHPnlZPF8ub6a8w5x8DAbGA0P35M++6XLpKman0i27sgUY/H6pko+j0Cygksbal6k8BYCWo6ndcXHhqvjP/cxY3Tef8Rsp2bTq6F8koyEf2k6nRx3nyZYazcBz96n2pd4Tr9OKOvkPhW8tR7OxfoNepZu6L2tISfI+Pp9js3PTGBTuOY3fSXF1MjfPUlB68wz1OiM+t2wIZgs5acA/D1B+Mfo7S6EhAGx6hx0Vh3yYfxz02peYeOg6BQj5ky08Jswv9d+13KfhywplQkC6HGo3Io7w//5vsgcQYNJodHwIuPNnEpPBQAFiyGwGrUd8bvsEykpiUPPieqDVXZyMYW83v+RwWlvH5nI/NdpQP9TmbjZXu1WRHEbh5IkFVHm7V2L2qPOjQL2uxd+/Qc9U/baP2H/5rl+AFiPse21mIgPOUceZTW03wfL5/T/Tour+FC008/Mi5E8uXu0fZHsYz8pstGbeA7qUoBGQDZQJAQFsyzFzAJBmVAqM+92kblZk9eYnlZSU9x+czgvjjm9tC/OkZLOzbR+TfO780fa8KDXkZnJYYshsqptd3NjrqN0E1liVs7RtkZB4naLBsyuBKGMXgbpdjYQ73nmfMfoUVfBRx2jZ3vWz/WSedINyjMSrdLFbjbZ8/vAsLkytxwF3z7Gsebu2izGh+j1oQcVfYnyxjPr8aARkA2VGQAAl97MGmUoWJq9jmYGUNKsP/sqyjUH/4/NSchXd+x1T8Hf/YXsKJsCTf8VTLFDt8BAw/FNThbS9uHmOY1pOLSVZunoADfuS0JoPZxyjIiA3g5XyV7cBV7fzggRYqxU8htZiUdLoBb3fzs+BA9O5CIz4gpakva5PxFG647ocuuON+1s+HzKH7lWLOyg2NC8CjjnDCnq/2pyqem0XkxMt73Dax3MUGgHZQJkSEGAqFgQATz8qWoNamIbUHZtHF23Am6aT99BMTlKt05nB64KyXrocWlN7v2dpw7CPqTdxNAZg0LMp1qUNwIX1XJUBZlUa9uGtfk9K/stD3CgthhmhiMOM60QeZeDfzYtFq00G0aJwdksRg4HtdLd+CKRGcBLJkA8da2p/bjVjRb7VKTK0rgU7+hetquYjGM8xJ5/kcHZcFC78TU4tZgOyXvkqkkoVGgHZQJkTEGDZuqNyLeCxjbwwDHq2czjxN8WEQz82Xdzn11I5W7kW8MAiSzGaGqKOs0gy8igvwFHfFN0cl5I1Rle3MTB6Yx9jRgBTvTXbUfRYsy0b8FdtROVtSYy0yUpimUvsBSp7Yy+wWX9qBLdx9eCx1O9B0qnfy/kdAhTcOEA3OeoY33PEFwzk2wuDnj2jdn9Fd/D+f/K7awd/Aza+QVf4vr8tewtlJlIgmRZDd/7w7wxY3/FdmS8KGgHZQLkgIIDlGEqT+oAGJCG/2lxRN7zOgtauUziyV/H1w49wOkZuBqvlW48t+D0MBgaYt33ItGzHiazW9qtdvGM3GHjxhx+i2xdzihkfndnUD4/KQNWG1NEo7Tl8Amn1uXsbC1G9mWky6KkAN+h4n53MW5bxPi2GafGUCMvKeVcPILAFybh2R2awarVzrAFYURB9Etj5JaehVK5NiUPbCY71IcpI4HCAq9tZ1X7Ht5YZMik5yXT312zlOv4PSyLNSmJ/8LiL/F1D/uSEi/sXlotuARoB2UC5ISApmb06uZD/B7bgFE6fQD635T0SVNsJ7IyoXFSpUUy9RxyhmHDwe4VXjWckcJU9Mpvbdp/KSn1n9tHR6+im/dfH5xoLIdNjWQSZEcf+w47AvRLgFUDXxL8u28b612UcqnortuoozZYSUSfo3l5cT8uv13NAj2kFx+XUEHGUvb/Tb3ICbqdJlhaLQU/r9eifdOlG/2D5G2encOrtzTMUlYbMAep1ByYuc/xYSggaAdlAuSEggKv90slsKAYA1YPZeL5yDZLQ3u+MPYP6APf/bSIMXQ4naBz9kwHicb/bV/WeFMo6slOL2fKzy6PspFgaimgpmZ3JTafITpdDhbheRxJxcWObDld3WkneASVvydgDKakkPjCdoj4vfw4V7D7V8Rnqeh1jc7u+oCs9YV5+gWBOGuNBF9dzkRj8viU55aQB88fT7ev+FGUUQS3ZxKwM0u22oBGQDZQrAgKoXF06meY8AFRrCjyy2kQop5YCq55hjGjiUssg6vEFVDO7urMFbPAY+97z5lmWZpxdzgu//QO8qIKaO/GDVXDkZbE/98EZdDd9qnOUcvcnHScegBbh8qnG4tB7OJnD2gJNDmNvoLjzLCzt/qTl8+YaoS6PUcfkX4+lGD7Viv5ZSwAaAdlAuSMgwEhCk7jqAXQxHlltmoIZuo+xHxc3Tu9s1M/02oSrwL9TuCJ2fBgY/pn9K2HidQrbjv9N96hhX1pFLe8snV7P5Q1Ssrzl+N/URWWnMLDeYxq1QkWxyAwG4OgcYMv7gHBlrKfdvfm3CzsILJpIq/jeP4Gmgy2fT41m18TEa+wJdWIhZQSPrKKLWs6gEZANlEsCAvJbQr41afEokzPirwCLHiDhDPmA8QfFNNfnUYOy5zua9qO/c0yImB7HcoSjf3IV9gli6r7tPUCtDmWeUSlxpEQYldELqaFy82Lgt/NkZhCL+vlvnmP6POIwJ3Tc9YvliByApHfod6qbA+qx7Yq1JZp4jf2iMxMZE1LcrodXljvLR4FGQDZQbgkIYIxgzfM0rQEGYe+bT1UyQP9/5TNU9AaPZXDaXMEbEcI0fuxZqmVHfGnZuKowGAzMyoTMYRGkIY9p9bb3ULhXo/WtQ0aKMvrcKkoVAOqsOkykgNBRAac5cjOBPd+wSZiXP61SNS1Wdgqw6lkeR/ORzGxa64fCj3DhMehphYX8yYzfQ8vKvhl/AdAIyAbKNQEBXBG3vs+TV8FdP3PlU57f/xOw9QNmge6exQtHgT6Pr931FVfy/q8zduGoS5WZyNEuZ5ZxIgYk4FeXqd7mw+kGevgU99OWHnLSqGG6uoMkm3CZj9fuSHINHlN8ZbTBwEkmWz+kirzDRGq51KyUqOO0eJPD81u0Ck4v44LjV4u/8Zl/gcYDuSiV8/IYjYBsoNwTkAJzsSLAIPGQD01p59C9DGqmx1A13edly1Rt/GVg45vAla20YoZ/SiVtUSyYtBhaRJc2sdVrbjrjGbXaU3hXvyfTwOWlgFVKulURR2gVRhzmBW/QUXvUsA9Fii3vUO/DXRSE7qMoMfoESW3450CDnvm30+exE8HurwHfGhw2UL9H/uPf9SXd6rrd6BJfXEdCu/PHcqHzKQwaAdlAhSEggBL9fx83DuADrY57/jSVYmQlsdXH2eUkgbEzTIFrBZe38MKIv8Qg86D/5T/hHYEuh5McQvcwcBoRYtL3+NZk0/2abVlZX60JRZbeVUrGdZOSkzcSr1ENHWdURceeo8YGoBVYuyM/c+OBvHdmej/iKNPqlzdzksfg94G296qLEmMvcI5b9Anqu0Z9ld+NsnDLRvD/sANA/zcoIq0gLrBGQDZQoQgIoNL4nwdMpQb+9Sz1I1JS17P+Na6uA96ktWQu0NPnMa6z+2sKApsMAga87ZzWE7ocWhcRIRTGxZymOtdg1tva05+Whl8dlmj4VGOdWqWqxvE8ys3DdLyKKlqXbVRFJ1EZnZVEMaaijFYGFgKMmQU2p0ixdifOZKvZtmQsBnPi8a5KF6r7U+pCQF0O3eZdX3PIwOjv1SUTMacpMk26wbjbjQMk0jt/KHgaSjlEhSYgIcRiAEpFXgCAZCllB5XtQgGkAdAD0NnzgSscAQFUEy+aSFdCwYgveMIrK2JKJEno4jpedHf9zJXfHLkZVEPv+5GTSJsMIlk1GeTclVWXyxhL4nWKH5VbegxV2ZnxJqvOESjKaL9aZsroeiS3oJa0tkpyLLPBAFzexHatoXtMxNPtCdsxmeu7aaUmXGbiYNQ3+V1VKZl4WPcKP1+r0dR4VarKeI95jK+CoEITkDmEEN8CSJFSfqTyXCiALlLKeHv3VyEJCOBFvfV99gdS0HQoMyfmOpBzq0lEGbFMIw98J3/1fE46a80OzuAKG9QK6PkM3YKSKtw0h5SMI2UmmqajKvcALRYXd1pxrsbBhWWpjM5JZ/fIgzNYbuJXh0roLo/ZJp7UKGp/Ti+heHTUt0CzIfm3y0wE1r3MLgl1uzIudGEti2gnzC2XGh97cEsQkBBCAAgDMEhKeVnl+VDcLgSk4PxaZkZyUvi/TxAwZrrluJWsZAYwD89ipqrfa7xgrC9gXQ5bkR74lTPtvasyXdzpYabcb2dIyWLb4/OBMyuAvAz2WO75DLNmtty67FRamAd+Zd+n3i+y17daK9Yr24BV0+gWB4+hCxZ/ia14B/6vQotBbxUC6gfgO1sfRAhxHUASOEH1dyml6tS+Uh3NXBpICgWWPWbSrwAspxj+maWOJO4iBW6XNzNlP/Adakmsi1elpEtxZDbnjhvyGEPpOJFuQ2FTN24lJFzlDK0T/9Bt8vClpqrTpIJjZrpcdh7Y+QVdzDZ3s1BYrf9QThpT9UdmMWZVrzuV155+bMNqrYKugCj3BFTQaGYp5SrjNjMAXJFSfmtjH7WllFFCiOoAtgB4Tkq5u6D3rfAWkAK9jsWqOz41PeYTxOrq4LGWMZ0rW4HN71GgGNiC44SDx6nHSzIS6DYcm8/tlUZXwWNYouGIsLGiIP4ySefcKlohAEcedXqY36Wnr+3X5mXTStr3I5ASzsLhYR/ZjtucX0sXOS2aavWMeMb2mgwm+VRQl8sa5Z6ACoMQwg1AJIDOUsoIO7b/AEC6lPKbgra7ZQhIQdQJtmGNO296rPkIWkPmwjqDgandnZ8zXV09mKZ+6/Hqpr6UzGqdW8UyhYTLAAQD200HM3Bdt2uF0KTkQ3YKxZVKu9akUD5er7uxXeudhbeezc1g98p9P5JM6naj6LPpEPWAfkoEsP51JgmqB/N2cT31VMM/yd+So4LjViCgEQDeklL2t/G8DwAXKWWa8e8tAD6SUm4saL+3HAEBpjase74DvVGwcLXXc0DfVy1XcIOeAc/dX5OIKteiUrrzZNttRKXktudW86KNOMJxPh6VKUSs15UXYJ3OBVsLZYWUSFoZ4cZb1HHGZzx8qa1SRIn2NGpLukHX6dg8ElnDvoyxNeqnTiC5GSz43fcjv8eWo5gdjDrGJMKdPxQ+F74C4lYgoL8AHJRS/mb2WG0Af0gpRwkhGgNYYXzKDcBCKeWn+fdkiVuSgBTEXWTK98Ze02O+NRmLaH+/ZexHSgZBD/zCfjfulSig6zSJGqOCVuOsZKaXr26jViX+Ih8XLsyo1WhNLU71YN771yvZ9LiC3ExmqW6eo/sYe54N25VJJG5ejG3V70Errm43+wK9Bj1wfRdjZBfXAxBA8F0cmVS/u43XGNjcf9tHtJAa9mUHx6vbGVMb9knRenVXEFR4Aiop3NIEBJBYTi5i8DnTLDkY2JxB6FZ35SeDmDNMMZ9dDuRlkjg6PcK0vD3V1llJFOZFHAYij1n2aAZ44fnVMXUzVASJ3lVodXlXIQG6etClc/WgBafPpcpal8u/c9OpfM5MBLISGT9JiWAFf0o4s0kKXNzZ3L16MAm1XjegRlvHMkvxV9ix8uQiIDWSWcIujwJdHrfdAE5K4NJGjlqKOcV+2YHNSdbZqbQ2B7xZvGLXCgCNgGzglicgBdkpwN4fmA42b4Vasx2JqPlwlersVBY8HptHF8HFDWjUn72nW452bNJDdoqxYfw5xlhSwll4mRLO2jI44Zxz8zKRmnKr0oikU61J0eJTidepwzm3iq6mcGFcp8ODrFi3pZNSiGfn5+wZHWAUR8ZdAJJv0AIa+VWZzekqbWgEZAO3DQEpSIkAtn9CIZ05qgcDvZ6n5F/tQo05w0r4sytIIAoZNRvGC7Jak6K7DwY9CSoz0VhikUjLS59nvOWyFMPVnUJEV3dqmDx8jJZTVd47o/+xQU/CuLyF1f83jVmwmm1pAbabwFlbtqDPI1nt/5k1XgENWKQbf4nkU6Mtm9bbCk7fotAIyAZuOwJSEHOaLTrOr7Z83K8u0ONpan7U+stIyQv07ApaBglX+HhAA15Ujfoxe1QafaWdASlZvHp9Fyv7r+82jhgSjA21utNYJd+w4P1kJVP3c+h3umdVGgJVG7NzYdx5WmKD/scsY2nEv8oZNAKygduWgBTcPMsM2NkVlo+7eVE81+UxZrNsrdaJ1xnPuLKNEzjzMvh4QH0SUd1uxmr44KL1TnYmpGQdXcxpIDLENLAwK5HP+9VhhXzjAZxEWpgGR0q+/tg89upR1NGVqpF0ksOos+rzIgP6FVGi4CRoBGQDtz0BKYi7yPjQqcWW87wAug0dHqRquiD3Q5fLizv8IFtzhB8ytcEAmP2qHgwENqOFEFCfllNAfeeNj5GSwWelOj4lglm5uIt0gZQBihDMyNXpzCr5Bn3sdyMz4tmg/vh8xrTcK7HdiKsHXbbsFO63z8tAi1G3pcVjDY2AbEAjICtkJADH/mLdWFq05XPChcHTtveyQruwFqBS0h25edZ0iz1Hl8ea5Dx8GdSuVI03rwBjSw5P001KU1sOQx4JLzuZF7zSmkNt3ph3FcoBgloY5QCt2NvakbE1SkfIsyvopkk9UL01s4I56cYgvTvT8V2nsB/TbRTjKQwaAdmARkA2oM9jHdiJBSzdkAbL54UrRYfNR7CEwJH2pQYDq/OTbjAjlBxmTKeb3bJTSDDmlfHChcFvV2VmmAfdOu8AEpZ3ADU1ysBC/7qMaVWq6jgZSMk6sCtbmc1SSCegAd/DoOfzuWl8jy6PUqpwi5ROOBsaAdmARkB2IDWK2pcTC0xBZ2tUa8oAdMM+dGcqYo1Yeixdx+u7mAVLNhYp+9cn0bl5UEWdHkPVd/AYZsUa9il8Gu1tDo2AbEAjIAdgXg92brVJ8awGpaq7dgegVkeqoUujt5C9yMsyKqNPsRwj7ABdQwBw92Fcys3DZAnlZfDxpoMYnG8+Qr2lhgZVaARkAxoBFQOxF1hMeXUHLQfztqvWcHGjEC+wOYPQ1ZryVrUR3aeS7BGdFGqaUR93gZqmhMsmt9K7CuNOLm6MMWWn0kUEWBvXYiTQ4g5aOuWJRCsQNAKyAY2AnIT/xtxsZz1Y7Nn8cSNbcPNmdq1yLd57V2FnQS8/9sTx8GHMSbiQqIRgHCkvk5ZMXiZvmQkMQmfE8z7tJmM05vCuauo3LSVFjtmpJvlApWokmoZ9KbQMbKYFk50ARwjIrfBNNGiwgmdlo6VgnMCanUqNTMQR3m6eZUZMDbosWidJ14t3DC5uJqW0QU9rxsufj7m4MpCtaH4APl6zLcWGtTuyMDWwuZY2L2NoBKSh+PDyA5oM5E1BVhLjLjfPsjQhOcyUCcvLLNr7uHoCkKZyDYMOMPcEddnsq+xX2+T+BTYHqjWjUrkCtzm9VaERkIaSgXcVpu4b9LJ8XEq6Tek3jWn4eOqRlJowXQ5dLV025QEubrRohCutFXcfEp6XP102L3+6cn61ST63sQK5IkIjIA2lCyE4mqa8TFLVUKbQHGANGjSUGTQC0qBBQ5mhTAhICHGvEOKsEMIghOhi9dxbQogrQoiLQojhNl5fVQixRQhx2XhfSKGSBg0ayiPKygI6A2A8AIsxOkKIYAD3A2gNYASA6UIINd37mwC2SSmbAdhm/F+DBg0VDGVCQFLK81JKNW3/GACLpJQ5UsrrAK4A6GZju7nGv+cCGFsiB6pBg4YSRXmLAdUBEG72f4TxMWvUkFJGA4DxXitL1qChAqLE0vD2TD9Ve5nKY8WqFbEazVycXWnQoMHJKDECklIOKcLLIgDUM/u/LoAole1uCiFqSSmjhRC1AMQWcBwzAcwEWAtWhGPSoEFDCaG8CRFXA1gohPgOQG0AzQActrHdJABfGO9tWVQWOHr0aLwQIgNAfKEbVzwEQvtcFQW34mcCTJ+rgb0vKJNqeCHEOAA/AwgCkAzghJRyuPG5dwA8BkAH4EUp5Qbj438A+E1KGSKEqAZgCYD6AMIA3CulTMz3RurvHWJvpW5Fgva5Kg5uxc8EFO1z3VbtOADtx69ouBU/1634mYCifa7ylgXToEHDbYTbkYBmlvUBlBC0z1VxcCt+JqAIn+u2c8E0aNBQfnA7WkAaNGgoJ7gtCKi4xa8VAUKID4QQkUKIE8bbqLI+puJACDHC+JtcEULcMrV+QohQIcRp429UYRuUCyHmCCFihRBnzB5zuEj8tiAgFL/4taLgeyllB+NtfVkfTFFh/A1+BTASQDCAB4y/1a2CgcbfqCJnwv4CrxlzOFwkflsQkBOKXzWULroBuCKlvCalzAWwCPytNJQTSCl3A7DW3jlcJH5bEFABsLf4taLgWSHEKaN5XJF7JN1qv4s5JIDNQoijxjrFWwkOF4mXt1KMIqO8FL+WJAr6jABmAPgYPP6PAXwLKsorIirU7+Igekspo4QQ1QFsEUJcMFoTtyVuGQIq4eLXcgF7P6MQYhaAtSV8OCWJCvW7OAIpZZTxPlYIsQJ0N28VArK7SFzB7e6CrQZwvxDCUwjRCLaLX8s9jD+4gnFg4L2i4giAZkKIRkIIDzBRsLqMj6nYEEL4CCEqK38DGIaK/TtZQykSB+wsEr9lLKCCYFX8uk4IcUJKOVxKeVYIsQTAObD4dZqUUl+Wx1oMfCWE6AC6KqEAppbp0RQDUkqdEOJZAJsAuAKYI6U8W8aH5QzUALBCcPyzG4CFUsqNZXtIRYMQ4h8AAwAECiEiALwPdqdYIoR4HMYi8UL3oymhNWjQUFa43V0wDRo0lCE0AtKgQUOZQSMgDRo0lBk0AtKgQUOZQSMgDRo0lBk0AtKgQUOZQSMgDRo0lBk0AtJQLiGE6GosrPUyKojPCiHalPVxaXAuNCGihnILIcQnALwAeAOIkFJ+XsaHpMHJ0AhIQ7mFsQ7sCIBsAL0qcJmMBhvQXDAN5RlVAfgCqAxaQhpuMWgWkIZyCyHEarAbYiMAtaSUz5bxIWlwMm6LangNFQ9CiEcA6KSUC409ovcLIQZJKbeX9bFpcB40C0iDBg1lBi0GpEGDhjKDRkAaNGgoM2gEpEGDhjKDRkAaNGgoM2gEpEGDhjKDRkAaNGgoM2gEpEGDhjKDRkAaNGgoM/wf+IZTTT3GCL0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.plot(shifttracker0_data[:,2], shifttracker0_data[:,3])\n", "plt.plot(shifttracker1_data[:,2], shifttracker1_data[:,3])\n", "plt.gca().set_aspect('equal', adjustable='box')\n", "plt.xlabel('x')\n", "plt.ylabel('y')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [], "source": [ "separation_vector = shifttracker1_data[:,2:5] - shifttracker0_data[:,2:5]\n", "separation_mag = np.linalg.norm(separation_vector, axis=1)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEGCAYAAABiq/5QAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAAsTAAALEwEAmpwYAAAkx0lEQVR4nO3de3xcdZ3/8dcnyeTW3JukDeklLW0ppUApKXcKgihUBBfkJqILunV3veD6c139uSq7v/1515+KuwKu3BTlIiCILCAgIBexLS3Q0jstbUrbpG3SpE3S3D6/P+akpKVNJpeZM5l5Px+P85iZM5Nz3pyUz5x8z/d8v+buiIhI+sgIO4CIiCSWCr+ISJpR4RcRSTMq/CIiaUaFX0QkzWSFHSAW5eXlXlNTE3YMEZFRZcmSJTvcveLg9aOi8NfU1LB48eKwY4iIjCpm9tah1qupR0Qkzajwi4ikGRV+EZE0o8IvIpJmVPhFRNKMCr+ISJpR4RcRSTOjoh//UG3e1crmxlaaWjvZtKuV46qLOfXIsZhZ2NFEREKT0oX/pmfXc9fLmw5Yd8a0cv7t4mM4sqIgpFQiIuGy0TARS21trQ/lzt2121vYsaeD4rwI44tzeXjZFn7wxzW0d3bzd2dO5VPzj6Q4PxKHxCIi4TOzJe5e+671qVz4D6WhZR/fenQlDyzdwpjsTD584gQuOLaKEyeXEsnUJQ8RSR0q/AdZubWZm59dz6PLt9HR1UMk05haXkBFYQ6FuVnkZWeSG8kkL5JJbiQjeIwuxXkRZo4vZGpFAZkZul4gIslJhf8w9u7r4rk1Dbxat5t19S3s2ttBc3sXbR3d7Ovqpq2jm/auHrp73n2cciMZzKoqYs7EUuZMKuGEiSVMKM3TxWMRSQoq/MPg7nR2O+1d3bR3drNzTwcrtzaz4u1mXqtr4vUtu2nv7AGgvCCbORNLgqWUaZUFlORHyI1khpZfRNLT4Qp/SvfqGSlmRnaWkZ2VQVFuhMrCXI6uKuKSudH3O7t7WL2thaWbm1i2qYmlmxt5cmX9AdvIycogq0+zkJlhBhlmZASPuZFMysZkU5IfoTQ/m9L8CEV5EYrzIhTlRijKy6IoN0JediZjcrLIi2SSHzzPycrQXxoiEhMV/hEQycxgdnUxs6uLueaUyQDsbu3k1bqm/fcRNLd1HtBc5ECPO+7Rx+4ep62jm12tHTQG9x007u2gZV8XsfxRZgb5kUzyc7LIz86kqjiXaZUFTK8s5LgJ0Wy6eC0ioMIfN8X5EebPeNfEN4PW0+Ps6eiiua2T5rYumts7aevoprWjm70dXfuft3Z0BY/d7N3XRV1jKw8te5uW9i4gej1izsQSTpoyljOmlTNnYgnZWfoiEElHKvxJLiPDos08uREoHdzPujvbmttZuqmJRRt3sWjjLn769Fp+8tRa8rMzOWXqWE6fVs4Z08qZMa5ATUUiaUKFP4WZGVXFeVQdm8eCY6uAaBPUS2/u5IV1O3hh3Q6eXhW9FlFRmMPpR47ljOkVnD5tLFXFeWFGF5E4Uq+eNLelqY0X1u3g+bXRL4KdezsAOLJiDGdMK+f0aeWccuTY6F8cIjKqqDunDKinx1m9vSX6RbBuBy+/uYu2zm6yMozTppVz/jHjOW/WOCoKc8KOKiIxUOGXQevo6mHppkaeXl3PY8u38dbOVsxgXk3Z/i+BiWX5YccUkcNQ4ZdhcXdWbWvhseXbeGz5NlZvbwHgqHGFvHdWJecePY45E0rI0BAWIkkj4YXfzG4FLgTq3X12sK4MuAeoATYCl7t740DbUuFPPht27OWpldt5cuV2Fm1spLvHKc6LUDu5lNqaMmprSjm6qoiCHPUfEAlLGIV/PrAHuLNP4f8usMvdv21mXwZK3f1fBtqWCn9y293ayTNr6nlx3U4WvbWLNxv27n9v8th8ZlUVcXSwzBxfSHVJnv4yEEmAUJp6zKwGeKRP4V8NnO3uW82sCnjG3Y8aaDsq/KPLzj37WLqpiZVbm3ljazMrtzazcWfr/vfHZGcyfVwhR40rZMb43scCKgpydC+ByAhKlsLf5O4lfd5vdPdD3pZkZguBhQCTJk068a233opbTom/vfu6WLWtmdXb9rBmewurt7Wwent0NNRepfkRZowr5KjxhcwYV8i0ygKmlI+hslBfCCJDMeoKf186409dO/bsY03wJbBme/RLYc22Flr2de3/TH52JjVjxzClfAw15fn7n08pH0PZmGx9KYgcRrKMzrndzKr6NPXUD/gTktLKC3Ion5bDadPK969zd97e3c76+j1s3LmXDTv2snHHXt7Y2szjK7bR1Wewu3FFOcyrKeOkKWXUTi5j5vhCXT8QGUCiC//DwMeBbwePDyV4/zIKmBnVJXlUl+QxnwMHuuvs7mFLYxsbduzlzR17eXVzdByiR17bCkS/SObPKOesGRWcOb2CsjHZYfwniCS1ePbq+Q1wNlAObAe+AfwOuBeYBGwCLnP3XQNtS0090h93p66xjZc37OK5NQ38eW0Dja2dmMHxE0o4Z2Yl588ez/RKDUQn6UU3cEna6O5xXqtr4tk1DTyzuoFX65pwh6nlY3j/7PGcf8x4jptQrC8BSXkq/JK26pvbeeKN7Ty+Yhsvrd9JV48zoTSPD584gctqJ1JdopFIJTWp8IsATa0dPLmynoeWbeH5dTsAOGNaOVfOm8R5s8ZpchpJKSr8IgfZvKuV+1+p477FdWxpaqO8IJvLaydy1UmTNPicpAQVfpHD6O5xnlvbwF1/2cTTq7bjwJnTK7j65EmcO7OSLM1VLKOUCr9IDN5uauOeRZu5Z9FmtjW3M64ohyvmTeLKeRM5QtcCZJRR4RcZhK7uHp5eVc+v/7qJZ9c0YMA5Myv5yMmTOGtGJZm6SUxGgWS5c1dkVMjKzOB9x4znfceMZ/OuVu5etIl7FtXx5MrFVJfkceW8iVw+byLjinLDjioyaDrjF4lRZ3cPf3xjO79+eRPPr9tBZobxnqMquGTuBM49upKcrMywI4ocQGf8IsMUycxgwbFVLDi2ig079nL3ok08+MoWnlxZT3FehAuPq+KSuROYO6lEN4dJUtMZv8gwdPc4z6/bwQOv1PH4im20d/YwpXwMl5xQzYdOqFa3UAmVLu6KxFlLeyf/s3wb9y+p4+UN0SGoTplaxiVzJ3DB7PEU5kZCTijpRoVfJIE272rld0u38MDSLWzYsZfcSAbnzRrPhcdVcdaMCnIjuh4g8Teswm9mpcARQBuw0d17Rj7i4anwy2jl7izd3MT9S+p49PWtNLZ2MiY7k/fOGscHjq1ivr4EJI4GXfjNrBj4NHAVkA00ALnAOOAvwH+5+5/ilrgPFX5JBV3dPbz05k4efX0rjy3fRmNrJwU5WZw3axwfPL6K+dMrdJewjKihFP4/AncCv3f3poPeOxG4Bnjd3X8x8nEPpMIvqaazu4eX1u/kD69t5fE3ttHU2sn4olwunzeRK+ZpxFAZGWrjF0lSncFdwr8J7hIGOGtGBVedpLGCZHiGcsY/t78NuvsrI5RtQCr8ki7qGlu5d3Ed9wZjBVWX5HH1KZO4ct4kTSMpgzaUwt8DrCDatg/Q944Ud/dzRjzlYajwS7rp6u7hqVX13PHiRl5cv5OcrAwuOv4IPn5aDbOri8OOJ6PEUAr/PwGXAruBu4EH3X1PXFMehgq/pLPV21q446WNPPjKFto6u6mdXMqnz5nG2TMqdIew9GvIbfxmNoVoz56LgbeAb7r7sniEPBwVfhHY3drJfUs2c9sLG9nS1MaciSV84bwZnDm9XF8AckiHK/wDXjVy9w3AQ8ATwEnAjJGPJyIDKc6P8Mkzp/KnL57NN//mWBpa9vGxW//K3962iPUNofwxLqNUf009U4EriZ7pbyba3POIu7cnLl6UzvhF3m1fVze/+ssmfvTkGto6urn29Bo+d+50DQ0h+w314u5rRM/2m4EDPujuP4xDzkNS4Rc5vB179vH9x1dzz+LNjCvM5T8+NJv3zhoXdixJAkNp6vl34EGgBygACg9aRCQJlBfk8O1Lj+OBfziN4rwIn7xzMZ/59Ss0tOwLO5okKd3AJZJCOrp6uPnZ9dz49DrysjP52oWzuHRutS7+pqlBn/Gb2b8Gg7Md7v1zzOzCkQooIsOXnZXBZ8+dzqPXn8H0ygK+eN+rXHv7IupbEn5pTpJYf238FwNfAtqBV3hnkLbpwBzgSaJdOxsOuYERpDN+kcHr6XHufGkj3/qfVYzJyeJblxzL+48ZH3YsSaDh9OOfDpwOVBEdlnkl8Jy7tw0jzD8BnyR6wfh14Nr+egup8IsM3br6Fq6/exkr3m7m8toJ3HDRMeRna9bVdJA0g7SZWTXwPDDL3dvM7F7gUXe//XA/o8IvMjwdXT386Mk1/OzZ9cyoLOTma06kpnxM2LEkzoZ8A1ecZAF5ZpYF5ANvh5RDJC1kZ2XwpfNncse1J7G9pZ0P/vR5nl61PexYEpKEF3533wJ8H9gEbAV2u/sTB3/OzBaa2WIzW9zQEPfLCCJpYf6MCn7/mTOYWJrPJ+5YzC3PrWc09OyTkZXwwh/0FLoYmEJ0OscxZvbRgz/n7re4e62711ZUVCQ6pkjKmliWz/3/cBoLZlfxzUdX8fWHVtDVndDZVCVkA17hMbMK4O+Amr6fd/frhrjP9wIbensDmdkDwGnAr4a4PREZpLzsTG686gQmlOZx83Nv8nZTGzd+5ARd9E0TsZzxPwQUE+2++Yc+y1BtAk4xs3yL3lVyLtGeQiKSQBkZxlcWHM3/+dBs/rS6nitu/ovu9k0TsXy957v7v4zUDt39ZTP7LdF7A7qApcAtI7V9ERmca06ZTHVJLp++aymX3fQiv/zEyUwsyw87lsRRLGf8j5jZgpHcqbt/w91nuvtsd7/G3XWaIRKic2aO41efPJnG1k4+fNOLrN3eEnYkiaNYCv/1RIt/u5m1BEtzvIOJSGKdOLmUez91Ku5w2c0vsXRTY9iRJE5imYil0N0z3D03eF7o7kWJCCciiXXU+EJ++/fRUT6v/u+XeX7tjrAjSRzE1J3TzC4ys+8HiwZmE0lhk8bmc9/fn8qksnyuu30Rz67RfTSpZsDCb2bfJtrc80awXB+sE5EUVVmYyz0LT2VaZQEL71zMS+t3hh1JRlAsZ/wLgPPc/VZ3vxU4P1gnIimsOD/CLz9xEpPK8vnEHYt4RW3+KSPWO3dL+jwvjkMOEUlCYwtyuOuTJ1NZmMN1ty9i4469YUeSERBL4f8WsNTMbjezO4AlwDfjG0tEkkVlUS53XHcSBlx3+yKaWjvCjiTDFEuvnt8ApwAPBMup7n53vIOJSPKYPHYMt3yslrrGNj71yyV0dGlsn9Gsv6kXZwaPc4lOwlIHbAaOCNaJSBqZV1PG9y47jpc37OI//vBG2HFkGPobsuELwELgB4d4z4Fz4pJIRJLWxXOqWb5lNz//8wZOnFzKxXOqw44kQ3DYwu/uC4OnFxw8LaKZ5cY1lYgkrS+dP5Nlm5v4ygOvM6uqiOnjCsOOJIMUy8XdF2NcJyJpIJKZwU8/Mpf87Ez+4a5XaOvoDjuSDFJ/bfzjzexEolMknmBmc4PlbKLTJYpImhpXlMuPrjiBdfV7+M5jq8KOI4PUXxv/+4G/BSYAP+yzvgX433HMJCKjwBnTy/nb02q4/cWNvPfocZwxvTzsSBIjG2i+TTO71N3vT1CeQ6qtrfXFixeHGUFEDqGto5sP3Phn2jq6eezz8ynOi4QdSfowsyXuXnvw+lj68d9vZh8wsy+Z2dd7l/jEFJHRJC87k/93+RzqW/bxXTX5jBqxDNJ2E3AF8FnAgMuAyXHOJSKjxPETS/jYqZP59V838VpdU9hxJAax9Oo5zd0/BjS6+78BpwIT4xtLREaTfzpvBuUFOXztd8vp7um/+VjCF0vh7+3D32pmRwCdwJT4RRKR0aYoN8JXFxzNq3W7eXDplrDjyABiKfy/N7MS4HtEJ0jfCPwmjplEZBS66PgjmF1dxI+eXKOxfJJcv4XfzDKAp9y9KejZMxmY6e66uCsiB8jIML74vqOoa2zjnkWbwo4j/ei38Lt7D33G6nH3fe6+O+6pRGRUOmtGBSfVlPGTp9fR3qk7epNVLE09T5jZpWZmcU8jIqOamfG5c6fT0LKPh5aprT9ZxVL4vwDcB+wzs2YzazGz5jjnEpFR6vRpY5lVVcTP/7yBHvXwSUqx3MBV6O4Z7p7t7kXB66JEhBOR0cfMWDh/Kuvq9/DMmvqw48ghxDTnrpmVmtlJZja/d4l3MBEZvT5wXBVVxbnc9sLGsKPIIcRy5+4ngeeAx4F/Cx5viG8sERnNIpkZXF47kefX7WDzrtaw48hBYjnjvx6YB7zl7u8BTgAahrNTMysxs9+a2SozW2lmpw5neyKSfC6fF73B/77Fm0NOIgeL6c7d3hm4zCzH3VcBRw1zvz8GHnP3mcDxwMphbk9Ekkx1SR7zp1dw7+I6DeOQZGIp/HXBnbu/A/5oZg8Bbw91h2ZWBMwHfgHg7h3u3jTU7YlI8rpy3kS2NbfzwrodYUeRPmLp1fM3wZ27NwBfI1qwPzSMfU4l2lR0m5ktNbP/NrMxB3/IzBaa2WIzW9zQMKyWJREJyXtmVlKQk8UfXtsadhTpI9ZePXPN7HPAcUCdu3cMY59ZwFzgZ+5+ArAX+PLBH3L3W9y91t1rKyoqhrE7EQlLbiST82aN47EV2+js1vg9ySKWXj1fB+4AxgLlRM/U/3UY+6wj+uXxcvD6t0S/CEQkBX3g2Cp2t3WquSeJxHLGfxUwz92/4e7fAE4Brh7qDt19G7DZzHovEJ8LvDHU7YlIcjtzRjmFuWruSSb9TbbeayOQyzvj8ucA64e5388Cd5lZNvAmcO0wtyciSSonK5OzZlTwzJoGenqcjAwN+xW2WAr/PmCFmf0RcOA84Hkz+wmAu39usDt192XAuyYAFpHUdPZRlTzy2lbe2NrM7OrisOOkvVgK/4PB0uuZ+EQRkVR11oxoB41nVter8CeBAQu/u99hZnnAJHdfnYBMIpJiKgpzmF1dxDOrG/jMOdPDjpP2YunV80FgGfBY8HqOmT0c51wikmLOmlHB0s1N7NnXFXaUtBdLr54bgJOAJtjfPq/J1kVkUE6aMpbuHmfppsawo6S9WAp/1yGmW9TAGyIyKHMnlZBhsGjDrrCjpL1YLu4uN7OPAJlmNh34HPBifGOJSKopzI0w64gi/rpRhT9ssZzxfxY4hmi3zt8AzcDn45hJRFLUvJoylm5qoqNLwzeEKZZB2lrd/avuPg84GfhO7zDNIiKDUTu5jH1dPazcqmm7wxRLr55fm1lRMILmCmC1mf1z/KOJSKo5bkK0D//ytw++bCiJFEtTzyx3byY6FPOjwCTgmniGEpHUNKE0j6LcLJZv0Rl/mGIp/BEzixAt/A+5eyfq1SMiQ2BmzK4uZoXO+EMVS+G/mehAbWOA58xsMtELvCIigza7uphVW1s0Pn+IYrm4+xN3r3b3Be7uwCbgPfGPJiKpaHZ1MR3dPazdvifsKGkrphm4+vIo3XMtIkNy1LhCANbWt4ScJH0NuvCLiAxHTXk+GQbr63XGHxYVfhFJqJysTCaV5bOuQYU/LLH04883s6+Z2c+D19PN7ML4RxORVDWtsoD19XvDjpG2Yjnjv43ocA2nBq/rgP+IWyIRSXlHVhawYcdeutSzJxSxFP4j3f27QCeAu7cBmjRTRIbsyIoCOrp72NzYFnaUtBRL4e8IZuByADM7kuhfACIiQzK1fAwAG3equScMsQzLfAPR2bcmmtldwOnAtfEMJSKpbWJZPgB1u1pDTpKeYplz9wkzWwKcQrSJ53p33xH3ZCKSsioKcsjOyqBOTT2hiKVXz1PuvtPd/+Duj7j7DjN7KhHhRCQ1ZWQYE0ry2NyoM/4wHPaM38xygXyg3MxKeeeCbhFwRAKyiUgKm1CWrzP+kPTX1PMpojNtHQEs4Z3C3wz8Z3xjiUiqm1Cax+t1TWHHSEuHLfzu/mPgx2b2WXe/MYGZRCQNTCzNp7G1kz37uijIiaWfiYyUWC7u3mhms4FZQG6f9XfGM5iIpLbq0jwA3m5qY0YwcJskRiwXd78B3Bgs7wG+C1w03B2bWaaZLTWzR4a7LREZfcYV5gBQ36zbghItlhu4PgycC2xz92uB44GcEdj39cDKEdiOiIxClUXRBoTtze0hJ0k/sRT+NnfvAbrMrAioB6YOZ6dmNgH4APDfw9mOiIxelb1n/C0640+0WK6oLDazEuDnRHv37AH+Osz9/gj4EnDYhj0zWwgsBJg0adIwdyciyWZMThYFOVnUt+iMP9FimXrxH929yd1vAs4DPh40+QxJMKRzvbsvGWC/t7h7rbvXVlRUDHV3IpLEKgtz1MYfgpj6UJlZNTC59/NmNt/dnxviPk8HLjKzBUR7CRWZ2a/c/aND3J6IjFKVRTk64w/BgIXfzL4DXAG8AXQHqx0YUuF3968AXwm2fTbwRRV9kfRUWZjLss1NYcdIO7Gc8X8IOMrd9feYiIyoykKd8Ychll49bwKReOzc3Z9xd03jKJKmygqyae/soa2je+APy4iJ5Yy/FVgWjMi5/6zf3T8Xt1QikhZK87MBaGztIC87L+Q06SOWwv9wsIiIjKjS/GhjQmNrB0eUqPAnSixj9dyRiCAikn72n/Hv7Qw5SXrpbzz+e939cjN7nWC+3b7c/bi4JhORlFc65p2mHkmc/s74rw8edfFVROKiJGjqaVLhT6j+xuPfGjy+lbg4IpJOept6dqmpJ6H6a+pp4RBNPL3cvSguiUQkbUQyMyjMyVJTT4L1d8ZfCGBm/w5sA35JdPrFq+lncDURkcEoGRNRU0+CxXID1/vd/b/cvcXdm939Z8Cl8Q4mIumhND+bxlY19SRSLIW/28yuDmbMyjCzq3lnzB4RkWEpzM1iz76usGOklVgK/0eAy4HtwXJZsE5EZNgKcrLY067Cn0j93sBlZpnAp9394gTlEZE0U5gboaVdTT2J1O8Zv7t3AycmKIuIpKGCnCxa1NSTULGM1bPUzB4G7gP29q509wfilkpE0kZR0Mbf0+NkZFjYcdJCLIW/DNgJnNNnnQMq/CIybAW5WbhDa2c3BTkxTQoowxTLIG1Dnl9XRGQgBTnRYRv2tHep8CfIgL16zGyCmT1oZvVmtt3M7jezCYkIJyKprzA3Wux1gTdxYunOeRvR8fiPAKqB3wfrRESGraC38OsCb8LEUvgr3P02d+8KltuBijjnEpE0UZjTe8avwp8osRT+HWb20eDO3Uwz+yjRi70iIsNWmPtOG78kRiyF/zqid+5uC5YPB+tERIZtTE4mAHv2qY0/UWLp1bMJuCgBWUQkDeVnR8tQa4eGAEsU9eoRkVDlZ0fP+Ns6VfgTRb16RCRUOVkZmEGbzvgTRr16RCRUZkZeJFNNPQmkXj0iErr87Ew19STQYHv1bEW9ekRkhOVlZ6qpJ4ES3qvHzCYCdwLjgR7gFnf/8UhtX0RGn/xIFq0d6sefKLH06rnDzEr6vC41s1uHsc8u4H+5+9HAKcCnzWzWMLYnIqNcbnYmbZ09YcdIG7E09Rzn7k29L9y9EThhqDt0963u/krwvAVYSbS3kIikqfxIJm0640+YWAp/hpmV9r4wszJiG8d/QGZWQ/RL5OVDvLfQzBab2eKGhoaR2J2IJCld3E2sWAr4D4AXzey3RCdguRz4v8PdsZkVAPcDn3f35oPfd/dbgFsAamtrfbj7E5HklZut7pyJFMvF3TvNbDHRGbgMuMTd3xjOTs0sQrTo36UpHEUk2tSjwp8oMTXZBIV+WMW+l5kZ8Atgpbv/cCS2KSKjm5p6EiuWNv6RdjpwDXCOmS0LlgUh5BCRJJGXnaWmngRK+ASX7v480SYjEREA8iKZdHT10N3jZGaoPMRbGGf8IiIH6B2hUzdxJYYKv4iELk9DMyeUCr+IhC4vEhR+tfMnhAq/iITunaYeFf5EUOEXkdCpqSexVPhFJHRq6kksFX4RCZ0mXE8sFX4RCV2eunMmlAq/iISut/C3q40/IVT4RSR0+RH16kkkFX4RCV2eunMmlAq/iIQuJyuDDFOvnkRR4ReR0JkZeRENzZwoKvwikhQ0NHPiqPCLSFLIz9aE64miwi8iSSEvonl3E0WFX0SSQp6mX0wYFX4RSQrRph4V/kRQ4ReRpKCmnsRR4ReRpKCmnsRR4ReRpKCmnsRR4ReRpBBt6lF3zkRQ4ReRpJCXnaWmngRR4ReRpJCfnUlnt9PZ3RN2lJSnwi8iSaEwNzoL1+62zpCTpD4VfhFJCpWFuQA0tOwLOUnqU+EXkaRQWZQDqPAnQiiF38zON7PVZrbOzL4cRgYRSS4VBdHCv725PeQkqS/hhd/MMoH/BC4AZgFXmdmsROcQkeRSXZpHXiST17fsDjtKyssKYZ8nAevc/U0AM7sbuBh4I4QsIpIkIpkZzJ9Rzi//8hZ/XruDDHvnPbN3XtghfjaVffOSY5lXUzai2wyj8FcDm/u8rgNOPvhDZrYQWAgwadKkxCQTkVB965LjqClfz9tN7fS4R1f6O+973xdpIi+YiH4khVH4D/WF/a7fprvfAtwCUFtbm36/bZE0VDYmm69ccHTYMVJeGBd364CJfV5PAN4OIYeISFoKo/AvAqab2RQzywauBB4OIYeISFpKeFOPu3eZ2WeAx4FM4FZ3X5HoHCIi6SqMNn7c/VHg0TD2LSKS7nTnrohImlHhFxFJMyr8IiJpRoVfRCTNmHvy3xtlZg3AW0P88XJgxwjGiRflHFnKObKUc+QkMuNkd684eOWoKPzDYWaL3b027BwDUc6RpZwjSzlHTjJkVFOPiEiaUeEXEUkz6VD4bwk7QIyUc2Qp58hSzpETesaUb+MXEZEDpcMZv4iI9KHCLyKSZlK68CfLpO5mNtHM/mRmK81shZldH6y/wcy2mNmyYFnQ52e+EuRebWbvT2DWjWb2epBncbCuzMz+aGZrg8fSMHOa2VF9jtkyM2s2s88nw/E0s1vNrN7MlvdZN+jjZ2YnBr+HdWb2E+s792D8cn7PzFaZ2Wtm9qCZlQTra8ysrc9xvSnknIP+PYeU854+GTea2bJgfWjHcz93T8mF6JDP64GpQDbwKjArpCxVwNzgeSGwhuhE8zcAXzzE52cFeXOAKcF/R2aCsm4Eyg9a913gy8HzLwPfCTvnQb/nbcDkZDiewHxgLrB8OMcP+CtwKtEZ6/4HuCABOd8HZAXPv9MnZ03fzx20nTByDvr3HEbOg97/AfD1sI9n75LKZ/z7J3V39w6gd1L3hHP3re7+SvC8BVhJdO7hw7kYuNvd97n7BmAd0f+esFwM3BE8vwP4UJ/1Yec8F1jv7v3d2Z2wnO7+HLDrEPuP+fiZWRVQ5O4vebQa3NnnZ+KW092fcPeu4OVfiM6Od1hh5exHUh3PXsFZ++XAb/rbRiJy9krlwn+oSd37K7YJYWY1wAnAy8GqzwR/Wt/apwkgzOwOPGFmSyw64T3AOHffCtEvMaAyCXL2upID/4dKtuMJgz9+1cHzg9cn0nVEzzh7TTGzpWb2rJmdGawLM+dgfs9hH88zge3uvrbPulCPZyoX/pgmdU8kMysA7gc+7+7NwM+AI4E5wFaifw5CuNlPd/e5wAXAp81sfj+fDfUYW3TqzouA+4JVyXg8+3O4XGEf168CXcBdwaqtwCR3PwH4AvBrMysivJyD/T2H/fu/igNPTkI/nqlc+JNqUnczixAt+ne5+wMA7r7d3bvdvQf4Oe80P4SW3d3fDh7rgQeDTNuDP0N7/xytDztn4ALgFXffDsl5PAODPX51HNjMkrC8ZvZx4ELg6qC5gaDpZGfwfAnRtvMZYeUcwu85zOOZBVwC3NO7LhmOZyoX/qSZ1D1o4/sFsNLdf9hnfVWfj/0N0Nsj4GHgSjPLMbMpwHSiF33inXOMmRX2Pid6sW95kOfjwcc+DjwUZs4+DjiTSrbj2cegjl/QHNRiZqcE/3Y+1udn4sbMzgf+BbjI3Vv7rK8ws8zg+dQg55sh5hzU7zmsnIH3AqvcfX8TTlIcz3hcMU6WBVhAtAfNeuCrIeY4g+ifbK8By4JlAfBL4PVg/cNAVZ+f+WqQezVxurJ/iJxTifaKeBVY0XvMgLHAU8Da4LEszJzBfvOBnUBxn3WhH0+iX0RbgU6iZ3CfGMrxA2qJFrT1wE8J7rKPc851RNvIe/+N3hR89tLg38OrwCvAB0POOejfcxg5g/W3A39/0GdDO569i4ZsEBFJM6nc1CMiIoegwi8ikmZU+EVE0owKv4hImlHhFxFJMyr8IgMwsxIz+8eD1j1mZtVm9oyZbeo7iqKZ/c7M9iQ+qUhsVPhFBlYC7C/8ZpZHtC/+lmBVE3B68F4J0dFYRZKWCr/IwL4NHBmMnf494GzgmT7v3030znCI3p7/QELTiQySCr/IwL5MdOjnOe7+z0THCHqsz/tPAfOD2/CvpM+4LCLJSIVfZPBOB57v87o7eH0FkOfuG8MIJRIrFX6RQQgG1drs0cl9+robuBG4N/GpRAZHhV9kYC1Ep8yEdzfz9Poz8C0GmGVJJBlkhR1AJNm5+04zeyGYSHsr8HeH+IwD3094OJEh0OicIjEysxzgBXevDTuLyHCo8IuIpBm18YuIpBkVfhGRNKPCLyKSZlT4RUTSjAq/iEiaUeEXEUkz/x9fe+VMNEldNQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.plot(shifttracker0_data[:,1], separation_mag)\n", "plt.xlabel('t/M')\n", "plt.ylabel('coordinate separation (M)')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Plot the real componenet of the $\\Psi_4$ data using columns 0 and 1 of the Ylm_WEYLSCAL4 files\n", "### Note that the initial ~(75 + extraction_radius) M will be junk radiation and should be cut off for most analyses. That time is marked in the following tutorial with a vertical dashed line." ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Stitching 17 files\n" ] } ], "source": [ "psi4_data = get_stitched_data(simulation_directory, 'Ylm_WEYLSCAL4::Psi4r_l2_m2_r75.00.asc')" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAakAAAEGCAYAAADfZmpgAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAAsTAAALEwEAmpwYAAA6hElEQVR4nO3de3yU1bnw/d81MzmShCQkQBLOgiJ4QEBArVSrtmBV3HXbam2hVqWtup+2++3exad1d+/92nfb47Nra7Xa2mJbi3b3IH2KWqV4FiUKgiBIgEgikAQIkAA5zMz1/jErySRMJjMhk5kw1/fzmc/MrHute66ZJHNlrXvd6xZVxRhjjElFnmQHYIwxxvTGkpQxxpiUZUnKGGNMyrIkZYwxJmVZkjLGGJOyfMkO4FRSUlKiEyZMSHYYvWpqagIgPz8/yZEYY0yXN998c7+qlkbaZklqAE2YMIHKyspkh9GrNWvWAHDppZcmORJjjOkiIu/3ts2G+4wxxqQs60mlkWnTpiU7BGOMiYslqTQyatSoZIdgjDFxseG+NNLY2EhjY2OywzDGmJhZkkojGzZsYMOGDckOwxhjYmZJyhhjTMqyJGWMMSZlWZIyxpjB0H4c1v8G7PJIcbEkZYwxg+H5/4In74D3nkl2JEOKTUFPI2effXayQzAmfTXVhe6P2wzbeFiSSiMlJSXJDsGY9CVu4EqDyY1jiLHhvjSyf/9+9u/fn+wwjElPlqT6xZJUGtm0aRObNm1KdhjGpCeR0L0lqbhYkjLGmMFgPal+SbkkJSILRGSbiFSJyLII20VE7nPbN4rIzL7aikixiDwrItvdfZErv0JE3hSRTe7+I2Ftnnf72uBuIxP93o0xpzBLUv2SUklKRLzA/cBCYBpwo4j0XLp7ITDF3ZYCD8TQdhmwWlWnAKvdc4D9wNWqejawBPh1j9e6SVVnuFv9wL1TY0za6Rjuw86TikdKJSlgDlClqjtVtQ1YASzqUWcR8KiGrAUKRaSsj7aLgOXu8XLgWgBVXa+qe1z5ZiBbRLIS9N6MMems8yReiVrNdJdqU9ArgJqw57XA3BjqVPTRdpSq7gVQ1b29DN1dB6xX1dawsl+KSAD4A3CP6omniovIUkI9OsaNGxf93SXZjBkzkh2CMUYsScUj1XpSkX56PRNDb3ViaRv5RUWmA98BvhBWfJMbBrzY3T4bqa2qPqSqs1V1dmlpaSwvlzRFRUUUFRUlOwxj0pQN8/VHqiWpWmBs2PMxwJ4Y60RrW+eGBHH3nceXRGQM8Cdgsaru6ChX1Q/cfRPwGKHhxCGtrq6Ourq6ZIdhTHqy4b5+SbUktQ6YIiITRSQTuAFY2aPOSmCxm+U3DzjshvKitV1JaGIE7v5JABEpBP4K3KWqr3S8gIj4RKTEPc4ArgLeGfB3O8i2bNnCli1bkh2GMWnKJSkb7otLSh2TUlW/iNwJPAN4gUdUdbOIfNFtfxBYBVwJVAHHgJujtXW7vhd4QkRuAXYD17vyO4HJwN0icrcr+yhwFHjGJSgv8BzwcOLeuTHmlGc9qX5JqSQFoKqrCCWi8LIHwx4rcEesbV35AeCyCOX3APf0Esqs2KM2xphY2bGpeKTacJ8xxpyiXHKy60nFxZKUMcYMho7cZCtOxCXlhvtM4syaZSOYxiSdJam4WJJKIwUFBckOwZg0ZsN9/WHDfWlkz5497NnT87QzY8yg6EhO1pOKi/Wk0si2bdsAKC8vT3IkxqQj7XFvYmE9KWOMGUzWk4qLJSljjBkMNtzXL5akjDFmUFiS6g9LUsYYM5gsScXFJk6kkblze16ayxgzaGy4r18sSaWR3NzcZIdgTBqzJNUfNtyXRnbv3s3u3buTHYYx6amzJ5XcMIYa60mlkR07Qtd0TPXL3BtzarKeVH9YT8oYYwaDHZPqF0tSxhgzGDqS0wv3wpr/Sm4sQ0jKJSkRWSAi20SkSkSWRdguInKf275RRGb21VZEikXkWRHZ7u6LXPkVIvKmiGxy9x8JazPLlVe517PLaRpj+i+8B/XCvcmLY4hJqSQlIl7gfmAhMA24UUSm9ai2EJjibkuBB2JouwxYrapTgNXuOcB+4GpVPRtYAvw67HUecPvveK0FA/dOjTFpx4b5+iXVJk7MAapUdSeAiKwAFgFbwuosAh51l5FfKyKFIlIGTIjSdhFwiWu/HHge+Lqqrg/b72YgW0SygGKgQFVfc/t6FLgWeGqA3++guvDCC5MdgjHpy5JUv6RUTwqoAGrCnte6sljqRGs7SlX3Arj7kRFe+zpgvaq2una1fcQBgIgsFZFKEalsaGiI8taSLysri6ysrGSHYUx6CgaSHcGQlGpJKtJxn55nFfRWJ5a2kV9UZDrwHeALccQRKlR9SFVnq+rs0tLSWF4uaXbt2sWuXbuSHYYx6UktSfVHqiWpWmBs2PMxQM+r9PVWJ1rbOjckiLuv76gkImOAPwGLVXVH2GuM6SOOIae6uprq6upkh2FMerLhvn5JtSS1DpgiIhNFJBO4AVjZo85KYLGb5TcPOOyG8KK1XUloYgTu/kkAESkE/grcpaqvdLyA21+TiMxzs/oWd7Qxxph+seG+fkmpJKWqfuBO4BngXeAJVd0sIl8UkS+6aquAnUAV8DBwe7S2rs29wBUish24wj3H1Z8M3C0iG9yt43jVl4Cfu9fZwRCfNGGMSTLrSfVLqs3uQ1VXEUpE4WUPhj1W4I5Y27ryA8BlEcrvAe7pZV+VwFnxxG6MMRH5W6G1KdlRDEkpl6SMMeaU84srYO/byY5iSLIklUYuvvjiZIdgTHqyBNVvlqTSiM9nP25jzNCSUhMnTGJVVVVRVVWV7DCMMSZmlqTSSE1NDTU1NX1XNMaYFGFJyhhjTMqyJGWMMSZlWZIyxhiTsixJGWOMSVk2JzmNXHrppckOwRhj4mI9KWOMMSnLklQa2bp1K1u3bk12GMakF43psnamF5ak0sjevXvZu3dvssMwJr0E/cmOYEizJGWMMYkUaEt2BEOaJSljjEkkS1InxZKUMcYkUsCG+05GyiUpEVkgIttEpEpElkXYLiJyn9u+UURm9tVWRIpF5FkR2e7ui1z5CBFZIyLNIvKTHq/zvNtXzyv2Dllerxev15vsMIxJL9aTOikplaRExAvcDywEpgE3isi0HtUWAlPcbSnwQAxtlwGrVXUKsNo9B2gB7ga+1ktIN6nqDHerH4C3mFTz589n/vz5yQ7DmPRiSeqkpFSSAuYAVaq6U1XbgBXAoh51FgGPashaoFBEyvpouwhY7h4vB64FUNWjqvoyoWRljDEDz2b3nZRUS1IVQPi1JGpdWSx1orUdpap7Adx9rEN3v3RDfXeLiESqICJLRaRSRCobGhpi3G1ybN68mc2bNyc7DGPSi/WkTkqqJalIiaDnmXC91YmlbTxuUtWzgYvd7bORKqnqQ6o6W1Vnl5aWnsTLJV59fT319UN+1NKYocWS1ElJtSRVC4wNez4G2BNjnWht69yQIO6+z29qVf3A3TcBjxEaTjTGmPgEA8mOYEhLtSS1DpgiIhNFJBO4AVjZo85KYLGb5TcPOOyG8KK1XQkscY+XAE9GC0JEfCJS4h5nAFcB75z82zPGpJ1Ae7IjGNJSahV0VfWLyJ3AM4AXeERVN4vIF932B4FVwJVAFXAMuDlaW7fre4EnROQWYDdwfcdrikg1UABkisi1wEeB94FnXILyAs8BDyfwrRtjTlVBS1InI6WSFICqriKUiMLLHgx7rMAdsbZ15QeAy3ppM6GXUGbFFvHQkZmZmewQjEk/NrvvpKRckjKJc9FFFyU7BGPSj604cVJS7ZiUMcacWtxw3zfaP5/kQIYmS1JpZOPGjWzcuDHZYRiTXtxwX2Xw9CQHMjTZcF8aOXDgQLJDMCb9uNl9fmzdzP6wnpQxxiSS60lZkuofS1LGGJNILkkFLEn1iyUpY4xJJDfc166WpPrDjkmlkZycnGSHYEz6seG+k2JJKo3Mmzcv2SEYk34sSZ0UG+4zxphEstl9J8WSVBpZv34969evT3YYxqSXoCWpk2FJKo0cOnSIQ4cOJTsMY9KLu1THCUlqr51YHwtLUsYYk0huuC/Q8+v2ZxcnIZihx5KUMcYkUrAd9WQQ+eLhpi+WpIwxJpGCfvDYROr+sk8ujeTn5yc7BGPST8AP3l6+alVBrIcVTcr1pERkgYhsE5EqEVkWYbuIyH1u+0YRmdlXWxEpFpFnRWS7uy9y5SNEZI2INIvIT3q8ziwR2eT2dZ/I0P9Nmj17NrNnz052GMakl6AfPBmRt9ml5fuUUklKRLzA/cBCYBpwo4hM61FtITDF3ZYCD8TQdhmwWlWnAKvdc4AW4G7gaxHCecDtv+O1FgzAWzTGpJtgO9rbcJ9dWr5PKZWkgDlAlaruVNU2YAWwqEedRcCjGrIWKBSRsj7aLgKWu8fLgWsBVPWoqr5MKFl1cvsrUNXX3OXqH+1oM5RVVlZSWVmZ7DCMSS/+VtSXHXlboG1wYxmCUi1JVQA1Yc9rXVksdaK1HaWqewHc/cgY4qjtIw4ARGSpiFSKSGVDQ0Mfu02upqYmmpqakh2GMenF3wLerMjbbLivT6mWpCId99EY68TSdiDjCBWqPqSqs1V1dmlpaT9f7iS89zfYumrwX9cYExt/K+qzJNVfqZakaoGxYc/HAHtirBOtbZ0bwusYyquPIY4xfcSRGh67HlbcmOwojDG98begrif1h8CHum/TQBICGlpSLUmtA6aIyEQRyQRuAFb2qLMSWOxm+c0DDrshvGhtVwJL3OMlwJPRgnD7axKReW5W3+K+2hhjTET+1s4k9f+0387Wc+/q2ha0JNWXlDpPSlX9InIn8AzgBR5R1c0i8kW3/UFgFXAlUAUcA26O1tbt+l7gCRG5BdgNXN/xmiJSDRQAmSJyLfBRVd0CfAn4FZADPOVuQ1phYWGyQzAm/fhbCGYUdD7V8KMJGkxCQENLXElKRK4APgncr6obRGSpqj40kAGp6ipCiSi87MGwxwrcEWtbV34AuKyXNhN6Ka8Ezoo17qHgvPPOS3YIxqQffyua3TW7Lxg+gGU9qT7F25O6nVDP5ZsiUgzMGPCIjDHmVNLWTNDblaQ0PEnZMak+xXtMqkFVD6nq14CPAucnICbTH43VfVZZu3Yta9euTXwsxpguxw8RzCrsfBoMX7zGelJ9ijdJ/bXjgaouI3SSq0kFP/twn1WOHz/O8ePHByEYYwwAwSC0HsGfNbyrKPxkFutJ9SneJHVMRB4WkRnueesAx2P6q+VQsiMwxvTUsBU0SHvBuK4ysWNS8bBjUsYYkyjVLwFwrOJCYCcAQZvdFxc7JmWMMYlycCdk5tOW17U2QFDDe1L+JAQ1tMTbk+p2TEpE/mmA4zEJNGLEiGSHYEx6adoH+aMIateBKLWJE3GJK0mp6pM9nv94YMMxiXTOOeckOwRj0ktzPeSNJhg2qtd9uM+SVF9SbVkkY4w5dTTvg7yR3XtSWE8qHieVpETk8yJyycCEYhLtlVde4ZVXXkl2GMakj+Z6yB9NWI5C1XpS8TjZntT/BXpZg96kmra2Ntra7CJrxgyK1mZoaz6hJ9V9WSSb3deXk1pgVlXrCS3oaowxJlxzXeg+b3S3JBWwY1JxiXeB2Ss7HgK3Ag+7RV2NMcaE60hSNrvvpMQ73PefwHSgBMh198YYY3pq2he6zxvVbSkktWWR4hJvkvowkAe0AJtV1dbuSyWN70fdPHLkSEaOHDlIwRiT5prdBcDzRhMMy1IavsrE3o09spbpKWqSEpHF7vYJEZmhqkdV9VvAfkIXHDSp5PGbom6ePn0606dPH6RgjElzzfvAkwE5Rd16UhKelF76Pmz586CHNpT01ZO61N0+BSwXkQ0iMltVV6vqNxMRkIgsEJFtIlIlIssibBcRuc9t3ygiM/tqKyLFIvKsiGx390Vh2+5y9beJyMfCyp93ZRvcLfW7IAd2JjsCY0yHpjrIGwkeD6q99KQA6rcOcmBDS9SJE6p6c/hzEfk48GcROUdVDw50MCLiBe4HrgBqgXUistJdzr3DQmCKu80FHgDm9tF2GbBaVe91yWsZ8HURmQbcQOg4WznwnIicrto5UHyTu0Lv0BCIPr38xRdfBGD+/PmDEY0x6a25DvJGAT0vz9EjSfnsLJ5ooiYpEfnnCMWHgO+KyBZV/eEAxzMHqFLVne71VwCLgPAktQh41F1Gfq2IFIpIGTAhSttFwCWu/XLgeeDrrnyFqrYCu0SkysXw2gC/r8ERbI+6ORCwg7TGDJbG+hoavCM5HbrN7jvh3ChbCT2qvob78iPcPIRm9uUlIJ4KoCbsea0ri6VOtLajVHUvgLvvGLrr6/V+6Yb67hYJnzfaRUSWikiliFQ2NDT09f6MMWkicGQflfszaWkPdD9PqmeSspXQo+pruO8/wp+LyBXAUuDLqpqIb+RIiaDn1Jfe6sTSNp7Xu0lVPxCRfOAPwGeJcCViVX0IeAhg9uzZNk3HGEOwrYUSOcI+LaahqbXbBL42f4+k1Mcwfbrra7jvEfcwB5jq6l+foAQFoZ7M2LDnY4A9MdbJjNK2TkTKVHWvGxqs7+v1VPUDd98kIo8RGgZMrSn3tqSKMSmpsb6GEcBeimlobu3Wk/L7ewy723BfVH0N973gbn8AbgHOUdVXExjPOmCKiEwUkUxCkxpW9qizEljsZvnNAw67IbxobVcCS9zjJcCTYeU3iEiWiEwkNBnjDRHxiUgJgIhkAFcB7yTiDZ+UOE8ELCsro6ysLEHBGGM6HD+wG4A6LaK5xU8gGGW4LxD9WHK662u4b/lgBeJezy8idxJaD9ALPKKqm0Xki277g8Aq4EqgitC5WjdHa+t2fS/whIjcAuwGrndtNovIE4QmV/iBO1Q1ICLDgGdcgvICzwEPJ/4TiFOc/4FNnTo1QYEYY8K1HwidWL9HR3C01U/4Ie1gzwlMNtwX1UktMJsIbi3AVT3KHgx7rMAdsbZ15QeAy3pp823g2z3KjgKz4o190Nm6X8akJE/DFtrUS7WOpqnVT35W11etntCTsiQVjV30cCiLc7hvzZo1rFmzJkHBGGMACPgZsftpNupp+PHR3OLvsXZfj4kTfktS0ViSGsrsgKsxqWfTE+QdreFB/9UAHGvzdz9PKmCz++JhSWooizTc19o0+HEYY0KCQXjph9QPO4PngqEV29r8wW5J6u2cOd3bWJKKypLUUBahJ3X8d0siVDTGDIpdz8OB7bxceiN5WRlk+Ty0+oPdzpPa461gYstvugosSUVlSWooi9CTyqlenYRAjDEAvPUo5BSzNvtDDM/pSlLhPan2QBAN/+p972l4/WdJCHZosCQ1lMV5TGrs2LGMHTu274rGmPi1HIFtT8FZn+BgKxTkZJCV4aXVH+g2ccIfiLAwzVP/OnhxDjEpNwXdxCHO2X2TJ09OUCDGGN75H/C3wNmf5PCqdobn+GhqaQ/1pMKylN9WiomL9aSGsPefvCeu+n6/H3/PdcOMMSfn2EF49cfw3H9AxWwYO4cjx/1Rhvtsic94WJIawsbv/F3kDa/dH7H4pZde4qWXXkpgRMakmcO18NCH4W/fhJLT4R9+BiIcaWmnIDuDLJ+X1vZg9+E+60nFxYb7hqJjB6H65d63P/O/Ye6XwGP/gxiTKEcaamj7+ZUUBg/iu/lpGH9B57amFj/52Rlk+jzumFTYcF9vPSlViHxFoLRmSWoIOvy7Wxle08csPn8LZOYOTkDGpIsDO+DdlbB/O2z+GzltR/iC79/42Zi5nV+mgaDS3OqnIMcXNgW9++y+iPwtkJGT+PcwxNi/2kPQodptfVdqP5b4QIxJB+0tULcFnrwDfnI+PPfvULWaHVrBp9u+weqjE9n4weHO6s0toeO++dkds/t6Dvf10pP6/ecS9x6GMOtJDUGBiNdq7MGSlDHxCwY5un83u9Y9xelHXiOz9lU46i6f582CObfBRV+hKbOE6/7jb3x67jjeXrubN3YdZOa4IgCOtIQuvZGf7SPT6+m24oRHogz3vfc0HNkDBeUJf5tDiSWpIcivfXeAG7a9Tunccd3KJkyYkKCIjBkC/K3QXA/H9ofOaWprhuY6OLgTDu6CgzvRg7sY5j/OWcB+KWL4WR8jo3QyFIyBSR/uTCBvvddAUGHB9DKe21LPe/u6liPrSFIF2RlkeAV/oKsn5fN4aOttuA9CcViS6saS1BAUy9yg0qduhbnXdyubOHFiYgIyJpH8rdDaDG1NtB8/Ql3DAUoy28kOHgslGrctdO+etzZB6xFoPULz4YMEjx2kgKOR9+/NguKJUDyJTVkzeWJnBpNmXMJ/vunja0VncOf8KSc0WbfrIF6PcN64Qk4fnc/WsCTV5Ib7CrJ9+Lwe/EHtPCbl9YSSVq82/BbGnA++zP5/XqeYlEtSIrIA+BGhiw3+XFXv7bFd3PYrCV308HOq+la0tiJSDDwOTACqgU+qaqPbdhehqw4HgP+lqs+48lnAr4AcQteo+rKGH/1MIo8GiTTi91zgPC73ru98ro9eS3D2rXinXQVAa2srAFlZWYMSpxlEqqEVSIIBCPpDJ3oH/e55WFmgHYJ+NNBO3aGjtLe3UZbvw0fQ1Y9wC4Q/b4dggEPNR9lZdxgNtDOxKIviHHH7DoTquNcJtW/vtr+jLa3sPdiEv72NwmwPI4d58KiLs6NuoL0r8QS7rlybAYzp7TPwZEBWHmTlQ1YBZObRQBEvN+WROexcatvz2c9wvnzNheQNHwGZeZA7AgoqwOMhEFTu/P7zjBqbxT3XX8jzTW/w6Gvv86VLJuP1dP+De6P6INPLCxiW5eP0kXm8vvMAqoqIcOR4x3BfqCcVPtzn9Ujvx6QglKQKKuAj3ziJX4ZTS0olKRHxAvcDVwC1wDoRWamqW8KqLSR0mfcpwFzgAWBuH22XAatV9V4RWeaef11EphG6zPx0oBx4TkROV9WA2+9SYC2hJLUAeCqxn0AMVBnB4Yib7vN/gheC5/L/ZvwKANm5Bu/ONbRc8FXqKGFzy2gCAT9nl2VRPrKETK+EZhR1fLmggIS+zDqWXNJgaGVn7biFtjUea2PfoeMcbQuQneGlaFgGZQXZeCIdLouY27Xb5vqmVvYdaaGlzU+G10Npfhblhdl4RbrV7WbfptAXUvbwUJWgH1D8AT/1ze3UHTpGa2sL+VleSvKyGDnMg6iG9tfxpR5+i1De5g9Qc6CZI8fbCAQC5GYII4ZlUpqXgYfwuj331/E6Hfvx09jcSnNrO8FAgCwv5GV5GZ7txSuEfgYaAG8meHyhHoF4XVIIdttPa7ufYMBPBn58BPFJEE/PaxT1QYDRcbXorhCYGfY8gBePLwPx+ELxezNCScPjA6+v8/Gh1iDvN7YREC9eXyZVBwPUHM3mnHEjycrMCtX3+EKfQ+awUNLJzOOgP4v7Xt5LuzeXy2ecxgvVx3itto2vfnwmC86bHKrn6/7PV3Orn6t/8AJFJZn8+Y4L2VF/lGt+8jLHd43lnmvPPuE9vfBePbsPHuNfF5wBwPWzxvBPv1tPZfVB5k4a0Vmv1R/g7ZpDfGbeeADGjcil1R+koamVkQXZXT2pHB8ZHg/+YNdwn9cTSlrhfuu/jArZzyXet0MFL343NOQ3++b+/4BOISmVpIA5QJWq7gQQkRXAIkKXd++wCHjU9WrWikihiJQR6iX11nYRcIlrvxx4Hvi6K1+hqq3ALhGpAuaISDVQoKqvuX09ClxLgpLU0UP7ObJvFxr0E3T//WrATzAYQIN+NODug36GffAao+UIVcFyvtp+O5t0Epd53uQXmT9gp5YxSfeesP/s1/4P44Gd478GwISN3z/pmIvcbSAIMMrd4m/sOWENQ48K5aKMUqEdH4ogKEclgwyflwyfF494Qm07b9LtcXtQONLip6k1CAqFXi8iHtqOKgcbhcNeL8V52RTnZePp1rb7/trVw+5DLew53E5QPeRk5ZGZ4aO+PcjRQ0EyvF4mluYzdmQeXq+vqyeROSz0vjwZ4PFyqCXAWzWH2dvcRobPw8jCPILio77Zz/7jQYbnZvORaeVUFOeBxxv6ohdv12OPl7ag8NfN+/n79oMUDMvl8unlZGZksm73EV7f3cRpIwv48kenUVowLNTOm9GVNDxemtqFbzy5hZd2HOLq88Zxy/zTycjM4uFX3ueXr+5mzoRifvbZWRQNizxU9evXqrn7yc3MnVjM/TfNpCQvizVb61n62FtU7M/h8S9cQHGEto1H2/jHB17lgJ7BH267kMkj87jIH2TxI6/z1WcamXx6JpOHnTg68N/PvkddUws//cxMsnxeppUX8Knzx/L4uhq+dMlkKgq7T/f+1avvM6ogi49ND6Xvj0wdSZbPw6pNe7slqXc+OEyrP8j5E4oBGFsUOtWjpvGYS1JdPSmfV/AHtLMn5fMIx8JO5q2XEXzDfwsA/zypkewdq1jq+yv8369Q//yDHB53BROnTGdvSyZ+TyZZOcNQbzbqCyX00O92b38b0SZYSd9V4pRXWMLwESfzr09kqZakKoCasOe1hHpLfdWp6KPtKNXQt7eq7hWRkWH7WhthX+3ucc/yhHj3xSeY/dZdMdf/vX8+X/cvJejOIFgdnMWElscAeFO7xs+/0nY7/535U/6auYD8s6/kSFsWATysPOtHNLYKG2sOs+twgHZ8jMjPYeb4YmaNG84540YgHg+NR9toPB5k894mntu6n+0NRwni4czy4Xz49FGcO66QwpwMjrYG2LW/mZerDvDi9v20+YOMH5HLR6eN5rIzS5kwYhgiHhRlT2MLr+3cz1837eOdPUfwCJw/sZgrzhzFtLICinIzOdoWYHt9M2u21fPS9v0cbwswfsQwrplRzmVnjqJ8eDYCtDU3cCz/NP76dg1/eOsD1n/QzLCsDK6ZUcH8ySXMnFBEaV4Wew63sG7XQX6z9n0q328kN9PLJ2ZW8Jl54zljVD7i/lIDQeWl7Q089vpunn23Dp9HWDSjglsvnsjU0QVA6ByXNVvr+fnLu3hj10FGHMvklosn8tl548nPzuj87I+2+ln+WjU/XbOD4+0Brp81htsvmcy4EV3nrr35fiM/+Ns2Xt1xgIqmHL5y+RQ+MXNMt6GlxqNt/ODZbTz29m6G52TwTwum8I/zxpHl8wKgqqx+t55vrdzM3WuPc+OcsSxbcCbDc7tiAXhtxwGW/XEj7x8oYfEF4/n6gqkMc5c0vwj4y9t7uOuPm/jr/7Rw341ncvGU0m7tdx84xi2/Wceu/co9n7iIG+Z0Tcr51jVnc974EXzt929z3QOv8qub53R7nwAr3tjN3U9u5vIzR3L/TTM747906kh+seR8PvfLN1jyyBs8dtvcbp9jS3uA2x6tpPbQcX5761wmj8wDINPn4Uc3nMfCH73EP/1uA3+6/UKyM7yd7d7de4RfvlrNDeeP65x1B3DHpZP5fWUt96+p4v/7h67e1M6GZl58r4F/vuJ0Mryhv6thWT4+fHopT2/ex7euno7H/VzWVTcCcP6E0H7HFoeSXc3B48waD0c6p6D7yPB6Qqueh/WkOpZF+ojnlzS2df2sf7itCLiJn/iv5XbfSi44splzt/wYtsBQWBZ67ZhbmHfrDwd8v6mWpCIOFsVYJ5a2sb5ezPsSkaWEhgUZN25cpCp9Gn3OZbzq+W88Yf/14vEh4gFvBuLx4vF6UY8Pf04JI3PLCD7yRsR91eiozoQFcNsNS1gw5XS8Xk/npeMvvfTSrvoHj/HCew28+F4DD713gOZNLcAHJ+z33DEV3HxVBQvOGk3Z8BNPOJx6Biy8CA4fb2fVpr38ef0H3PPSQe556RDDMr0U5GTQ1OKnuTX0Bzx1dCG3XTmdRTMqGFWQ3W1fJcD4CXD5nNCX/V837uWJyhr+bXUD/7a6gZwML8NzMth3pAWR91GFM8sKuPuqs7h+9phuX3IAFYU5VJxXwbXnVfDOB4f51avVPFFZy2/W7qZ4WCbTy0MJ6N29TexvbqV4WCZ3XDKZxReMZ2SP2DK8Hj46fTQfnT6addUH+cnfq/ju09t48PkdXDFtNKX5WdQ2HuP5bQ00t/q5/MxRLFs4tfPLNdys8UX89ta5vFy1n+89s41/+Z+N/PjvVXxk6kg8Iuw7cpy/b62nPaAsvmACX7l8CoW53XsaIsLl00Zx4eQR/Pdz2/nFy7t46p19fHrOOOZNGsHRVj+/f7OWv2+tZ/yIXFYsnce8sF5Bh6vPLeesiuF88ddvsviRN7jt4knc8qGJDM/J4M/rP+A7T28lqPDo5+dw4eSSE9pfc2455cOzufXRSv7hp6/w8yWzOW9cEcGg8uCLO/jeM9u45IzSbgmqwwWnjeCBz8xk6aNvcsvySh79/ByyM7y0+YPc8du3eHN3Iz+5cWZnz6XDqIJsvn/9OXz+V5V85+mtfOvq6QAEg8rdf36H4TkZfN0N3XUoL8zhU+ePZcW63dx+yWmMcb2g5a9Wk+EVbpjTPR0sOGs0f9tSx9u1hzjPJbt1uw4yqXQYI/JCvbeOfdQcDJ320dTSTk6GlwyvB587BtWxwKwv7B+QtqxCGo8dB2BS6TB2NhwlN9PLx2dM5d43hvHNj5/Jzkw/DXuqOaNIyJY22o4fwxtsxadteILtaG9fcVG++cRtHOgj7CWTZgzsDp1US1K1dP+nYQywJ8Y6mVHa1olImetFlQH1feyrlu7HZyPFAYCqPgQ8BDB79ux+/djHTDiDMRPO6Luic+hY6CJpHz69lBfea4hYRyQ0DDF96tSo+xpbnMtn5o3nM/PG0x4Isn73IdZVHyTT66EkP5MRw7KYMiovYmKKZHhOBjfOGceNc8ax9/BxXtjWwLa6JrdMjI9JJcO44LSSiF/akQzL8vHJ88fyyfPHsrOhmVeq9lN94FjomNjhFqaXF7BoRgVnVQyPaX9nVQzn+9efy10Lp/K3LXW8+X4j79U1ISJ8aPIILp82iiumjTrhizSS8ycUs/zzc9hUe5gHX9zBC+81cKSlndK8LBacNZpPz+3+X3wkIsLFU0r50OQSntm8j+Wvvs8TlTV4RMjL8nHdzDEsvmACZ4zOj7qf3Ewf//vKM1k0o5wfPbedB17YwU+f3wFAYW4G//KxM7j5ognkZvb+Jz+xZBh/uuNC/vMvW3j4pZ089OJOREJfZrPGF/H9689lYsmwXtvPnlDMH790IZ/75Tque+BVzp9QTH1TK7v2H+Wqc8r4/vXn9vq5fmTqKH74qRl8ecV6rr3/Fa6ZUc7T7+xjY+1h7rn2LD5+Tlmv7T534QR++Uo1s8YXcdU55fxo9XYq32/ke/94zglJHeD2S0/j8XU13L+miv/6xDkcPtbO79+s5ZpzKxiZ3/2fksumjsLnEZ7evI/zxhURCCrrqg92iyc7w0tpfhY1jaEkdeR4aLUJIDS7L6Bdx6S8XUkqJ6znd+6YQnY2HGX2hGLuufZsvnr56WH/IJ3e62eeDlItSa0DpojIREL/zt8AfLpHnZXAne6Y01zgsEs+DVHargSWAPe6+yfDyh8TkR8SmjgxBXhDVQMi0iQi84DXgcXAjxPyjvuhMDeTt//to1QfONprklpx2zzOHVvYrey0006Lut8Mr4c5E4uZM7E4ar1YlQ3P6TYsdLImleYxqTS25NaXEXlZncn0ZJ09Zjj3f3pm3xWjEBEWnFXGgrMifxnHanr5cB5aPJsDza1sr28mO8PLtLICMn2xLS6Tm+nj3uvO4daLJ/H3rXU0twaYNb6I+VNKOodFo5lUmsef77iIh1/ayStV+xlTlMNXrzidq88p67P9NeeWMyzTy3/8ZQvffXobZcOzuf/TM3tNUB2WLZzKxtpD/NPv1nPf6u28V9fMdTPH8I+zIs8DLBuew41zxvLb13dz09zxPFFZw/H2ALfNP/EUjeG5GVxw2gieeWcfyxZMZcueIxxp8Z/QGx1XnMtu15M6cLSNIpccM71CWyBIoGN2X9hnkJvZlaQunTqSP63/gEXnluP1yAk9+HSWUklKVf0icifwDKFp5I+o6mYR+aLb/iChmXZXAlWEpqDfHK2t2/W9wBMicguwG7jetdksIk8QmlzhB+5wM/sAvkTXFPSnSIWZfWGG52bgbez9j35kQXa3MXro/3CkGXpG5GV1Dkf1x+SReTH3dnsqHpbJ1xdE78H35rIzR/GRqSM50uInP8vXeRwomuwML7+9dR73/X07G3Yf4q6FY7j14klRk+KXLz+dZ7fUcc1PXiaocPNFEzqPO/a04KzRfONP7/BeXTMvV+0H4IIISeqNXQcBOHi0lRF5oSTlc8e3/IEgInR7Px29Wq9HuPqcMuZNLLbkFEFKJSkAVV1FKBGFlz0Y9liBO2Jt68oPAJf10ubbwLcjlFcCZ8UT+2Dree5GuILsE3+0x46F/tPLzbWFZ03qEhGG52T0XTFMTqY3rsRYPCyTx79wAT97cQcVhbncdnHvJ7pfMW0U3/zzO/zl7T08924d544ZfkIyGVecy583fECrP8DBo22cXVQIgM8N77X5g3ik+0y8jp7U8JwMRKz31JuUS1Imdp4o/ylGmsr7+uuvA90nThiTrsYW50Y8X6qnkfnZXHHmKH6ypgqA//rEiW3GFeeiCh80Hufg0TaK3ezKDHe5nFZ/EK9It7/ZbJekCuNMyOnGktQQ5o1ymCGW4wfGmNj8+zXTaW71U16Yw/URjnWNd1PudzQc5UiLn2J33lZ4T0qk+z+WuW44vsCSVFSWpIawSD2pu6+adsJ4uTHm5JQX5vDYbfN63d5xXthbu0PnUPU8JtUWcMN9YX+yHcN9hbmWpKKxJDWERTom9YnzKno9498YkxileVnkZHh5dccBAMYUhU7ZyHQ9qVZ/AE+PnlSOmzgR7/G3dGMXPRzCIvWkfF4b5jNmsIkI44pzebvmENCVpHzumFTHxAlP2Ddu+MQJ0zvrSQ1hkabnZkQ5UHXGGbGfMGyMic/Usny21TXhka5VKHydPanux6REIDsj9LcabQKUsZ7UkOaN8MsdLUmVl5dTXm4XVDMmETpO8C0bntN5jmLH32ObP4jHI50TmrwiVBSGEpn1pKKzntQQ5omQj6KdO3XkyBEACgoin7RojOm/q84pY83W+m4rXXSs1dc5ccKVezzC5dNG8s2Pnzmgq7KciixJDWHxDhO8+eabgJ0nZUwi5Gdn8NDi2d3KuvWkhM7rrfk8QpbPy60XTxrsMIccG+4bwiIN9xljUkdHkgodk+o6mdf+dmNnSWoIi2VdM2NM8nRfFqlr9MP+dmNnSWoIi3b8yRiTfBk91+5zf7L2txs7S1JDmP2eG5PaOs+TchMnOntSNtwXM5s4MYTF+4s+bdq0BEVijIkkfLgvN9PbOSM32rqbpjtLUkNYvEMGo0aNSlAkxphIus/us55Uf1g+H8LinSHU2NhIY2NjgqIxxvSU0W2B2a5yS1KxS5kkJSLFIvKsiGx390W91FsgIttEpEpElsXSXkTucvW3icjHwspnicgmt+0+caeDi8jnRKRBRDa4262JfO/9Fe8MoQ0bNrBhw4bEBGOMOYEv7G80vCdlEydilzJJClgGrFbVKcBq97wbEfEC9wMLgWnAjSIyLVp7t/0GYDqwAPip2w/AA8BSYIq7LQh7ucdVdYa7/XxA36kxJi2EL1MmYSfzWo6KXSolqUXAcvd4OXBthDpzgCpV3amqbcAK1y5a+0XAClVtVdVdQBUwR0TKgAJVfc1dkv7RXl7TGGP6JfyqBN2OSVmWilkqJalRqroXwN2PjFCnAqgJe17ryqK1761NhXscaV8A14nIRhH5HxEZ21vQIrJURCpFpLKhoaGv95hQ9ntvTGrJCFtg09tjgVkTm0FNUiLynIi8E+G2qO/WoV1EKNN+tom2r78AE1T1HOA5unpoJzZQfUhVZ6vq7NLS0j5CSayOlZeNMakhvCcVWhYp9NgmTsRuUKegq+rlvW0TkToRKVPVvW4orj5CtVogvFczBtjjHvfWvrc2te7xCftS1QNh5Q8D3+nzzaWAuROLo24/++yzBykSYwx0PyZlyyL1TyoN960ElrjHS4AnI9RZB0wRkYkikkloQsTKPtqvBG4QkSwRmUhogsQbbkiwSUTmuVl9izvauCTX4Rrg3YF4g4n0hy9dyP03zYxap6SkhJKSkkGKyBiT0fOYlJ3MG7dUOpn3XuAJEbkF2A1cDyAi5cDPVfVKVfWLyJ3AM4AXeERVN0drr6qbReQJYAvgB+5Q1YBr8yXgV0AO8JS7AfwvEbnG1T8IfC5h73qAzBofccZ+N/v37wewRGXMIBERvB4hEFQ8Ah1XlLLhvtilTJJyQ2yXRSjfA1wZ9nwVsCrW9m7bt4FvRyivBM6KUH4XcFcc4Q8JmzZtAux6UsYMJp9LUhK2wKwlqdhZp9MYYxKo47hU+DEpO5k3dpakjDEmgTpm+Hm6ze5LYkBDjCUpY4xJoK6elC0w2x+WpIwxJoEyXLdJhK6Tea0rFbOUmThhEm/GjBnJDsGYtOPr1pMKlVmSip0lqTRSVNT3NHVjzMDqPCbloXN2n9hwX8xsuC+N1NXVUVdXl+wwjEkrHev3dbtUh+WomFlPKo1s2bIFsCv0GjOYMnxdx6HsmFT8rCdljDEJ5HM9KZ+n65iUDffFzpKUMcYkUEa386TsUh3xsiRljDEJ1NmT8trsvv6wJGWMMQkUvuJExzCfdaRiZxMn0sisWbOSHYIxaadjxYnQMSmbOBEvS1JppKCgINkhGJN2Oo9JebpWQbdjUrGz4b40smfPHvbs2dN3RWPMgPF5bXbfybCeVBrZtm0bAOXl5UmOxJj00bF2n9fjCRvuS2ZEQ0vKfFQiUiwiz4rIdncfcQ0fEVkgIttEpEpElsXSXkTucvW3icjHwsq/LSI1ItLc4zWyRORx1+Z1EZmQgLc8IF78l0t5/muXJDsMY0wvOnpSXo8tMNsfKZOkgGXAalWdAqx2z7sRES9wP7AQmAbcKCLTorV3228ApgMLgJ+6/QD8BZgTIZZbgEZVnQz8H+A7A/IOE2DciFwmlAxLdhjGmF50HJPyeTw23NcPqZSkFgHL3ePlwLUR6swBqlR1p6q2AStcu2jtFwErVLVVVXcBVW4/qOpaVd3bRyz/A1wm9ltljOkHX8S1++zrJFaplKRGdSQMdz8yQp0KoCbsea0ri9Y+WpvedLZRVT9wGBgRqaKILBWRShGpbGho6GO3xph00zkF3U7m7ZdBnTghIs8BoyNs+kasu4hQpslso6oPAQ8BzJ49u6/9JtXcuXOTHYIxaSfDTuY9KYOapFT18t62iUidiJSp6l4RKQPqI1SrBcaGPR8DdMyp7q19tDa96WhTKyI+YDhwsI82KS83NzfZIRiTdnydx6S6zpOSiP8Hm0hSabhvJbDEPV4CPBmhzjpgiohMFJFMQhMiVvbRfiVwg5uxNxGYArwRRyz/CPxdVVO6lxSL3bt3s3v37mSHYUxa8XYek6LzmJT1pGKXSknqXuAKEdkOXOGeIyLlIrIKOo8P3Qk8A7wLPKGqm6O1d9ufALYATwN3qGrA7fu7IlIL5IpIrYj8u9vXL4ARIlIF/DMRZhoORTt27GDHjh3JDsOY9OL+v/WEncxrh6RilzIn86rqAeCyCOV7gCvDnq8CVsXa3m37NvDtCOX/CvxrhPIW4Po4wjfGmIiCbgzGG3ZMymNZKmap1JMyxphTTjCsJ9Vx1MCmoMfOkpQxxiRQR09KpHPkz8TBkpQxxiRQpi/0NdvSFugc5rNcFbuUOSZlEu/CCy9MdgjGpJ1Z40PLiE4ZlU9t43EAgkFLU7GynlQaycrKIisrK9lhGJNW5k8p4ZmvzOfqc8s7Z/UFbdwvZpak0siuXbvYtWtXssMwJq2ICGeMzge6lkOyjlTsLEmlkerqaqqrq5MdhjFpq2MKesCyVMwsSRljzCA5BRauGXSWpIwxZpDkZobmquVn25y1WNknZYwxg+Ta88rZ2dDMbfMnJTuUIcOSlDHGDJLcTB/fvGpa3xVNJ0tSaeTiiy9OdgjGGBMXS1JpxOezH7cxZmixiRNppKqqiqqqqmSHYYwxMbMklUZqamqoqalJdhjGGBMzS1LGGGNSliUpY4wxKcuSlDHGmJRlScoYY0zKEltLauCISAPwfj+blwD7BzCcRLE4B5bFObAszoE1WHGOV9XSSBssSaUIEalU1dnJjqMvFufAsjgHlsU5sFIhThvuM8YYk7IsSRljjElZlqRSx0PJDiBGFufAsjgHlsU5sJIepx2TMsYYk7KsJ2WMMSZlWZIyxhiTsixJpQARWSAi20SkSkSWJTGOsSKyRkTeFZHNIvJlV/7vIvKBiGxwtyvD2tzl4t4mIh8bxFirRWSTi6fSlRWLyLMist3dFyUzThE5I+wz2yAiR0TkK6nweYrIIyJSLyLvhJXF/fmJyCz3c6gSkftERAYhzu+JyFYR2SgifxKRQlc+QUSOh32uDyY5zrh/zkmK8/GwGKtFZIMrT9rn2Y2q2i2JN8AL7AAmAZnA28C0JMVSBsx0j/OB94BpwL8DX4tQf5qLNwuY6N6Hd5BirQZKepR9F1jmHi8DvpPsOHv8nPcB41Ph8wTmAzOBd07m8wPeAC4ABHgKWDgIcX4U8LnH3wmLc0J4vR77SUaccf+ckxFnj+0/AP4t2Z9n+M16Usk3B6hS1Z2q2gasABYlIxBV3auqb7nHTcC7QEWUJouAFaraqqq7gCpC7ydZFgHL3ePlwLVh5cmO8zJgh6pGW5Fk0OJU1ReBgxFeP+bPT0TKgAJVfU1D31yPhrVJWJyq+jdV9buna4Ex0faRrDijSKnPs4PrDX0S+F20fQxGnOEsSSVfBRB+kadaoieGQSEiE4DzgNdd0Z1ueOWRsGGgZMauwN9E5E0RWerKRqnqXgglXGBkCsTZ4Qa6//Gn2ucJ8X9+Fe5xz/LB9HlC/8l3mCgi60XkBRG52JUlM854fs7J/jwvBupUdXtYWdI/T0tSyRdpLDep5wWISB7wB+ArqnoEeAA4DZgB7CU0JADJjf0iVZ0JLATuEJH5Ueom9TMWkUzgGuD3rigVP89oeosr2Z/rNwA/8FtXtBcYp6rnAf8MPCYiBSQvznh/zsn++d9I93+kUuLztCSVfLXA2LDnY4A9SYoFEckglKB+q6p/BFDVOlUNqGoQeJiuIaikxa6qe9x9PfAnF1OdG4roGJKoT3aczkLgLVWtg9T8PJ14P79aug+1DVq8IrIEuAq4yQ054YbPDrjHbxI61nN6suLsx885mZ+nD/gE8HhHWap8npakkm8dMEVEJrr/uG8AViYjEDcm/QvgXVX9YVh5WVi1fwA6ZgatBG4QkSwRmQhMIXRANdFxDhOR/I7HhA6kv+PiWeKqLQGeTGacYbr9h5pqn2eYuD4/NyTYJCLz3O/O4rA2CSMiC4CvA9eo6rGw8lIR8brHk1ycO5MYZ1w/52TF6VwObFXVzmG8lPk8EzUjw25xzbi5ktBMuh3AN5IYx4cIdds3Ahvc7Urg18AmV74SKAtr8w0X9zYSOMOnR5yTCM2OehvY3PGZASOA1cB2d1+czDjd6+YCB4DhYWVJ/zwJJc29QDuh/4xv6c/nB8wm9OW7A/gJbhWbBMdZReiYTsfv6IOu7nXu9+Ft4C3g6iTHGffPORlxuvJfAV/sUTdpn2f4zZZFMsYYk7JsuM8YY0zKsiRljDEmZVmSMsYYk7IsSRljjElZlqSMMcakLEtSxpwiRKRQRG7vUfa0iFSIyPMisjt8tWoR+bOINA9+pMbEzpKUMaeOQqAzSYlIDqFznT5wRYeAi9y2QkKr3huT0ixJGXPquBc4zV3753vAJcDzYdtXEFrRBEJL4PxxUKMzph8sSRlz6lhG6HIgM1T1XwitGfh02PbVwHy31M0NhK3TZkyqsiRlzKnrIuDlsOcB9/xTQI6qVicjKGPiYUnKmFOQWxC0RkMX0gy3Avgx8MTgR2VM/CxJGXPqaALy3eOeQ30dXgL+iz6uvmpMqvAlOwBjzMBQ1QMi8oqIvENopevbItRR4PuDHpwx/WSroBtzihGRLOAVVZ2d7FiMOVmWpIwxxqQsOyZljDEmZVmSMsYYk7IsSRljjElZlqSMMcakLEtSxhhjUpYlKWOMMSnr/wc3YfeLlbYJXAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.plot(psi4_data[:,0], psi4_data[:,1])\n", "plt.plot(psi4_data[:,0], np.sqrt(psi4_data[:,1]*psi4_data[:,1] + psi4_data[:,2]*psi4_data[:,2]))\n", "plt.axvline(x=150, c='#a9a9a9', linestyle='--')\n", "plt.xlabel('t/M')\n", "plt.ylabel(r'$\\Psi_{4, 22}$')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "mayawaves_dev", "language": "python", "name": "mayawaves_dev" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.15" } }, "nbformat": 4, "nbformat_minor": 4 }