{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "scrolled": true }, "outputs": [], "source": [ "from devito import *\n", "from examples.seismic.source import DGaussSource, TimeAxis\n", "from examples.seismic import plot_image\n", "import numpy as np\n", "\n", "from sympy import init_printing, latex\n", "init_printing(use_latex='mathjax')" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "scrolled": true }, "outputs": [], "source": [ "# Initial grid: 1km x 1km, with spacing 100m\n", "extent = (2000., 2000.)\n", "shape = (81, 81)\n", "x = SpaceDimension(name='x', spacing=Constant(name='h_x', value=extent[0]/(shape[0]-1)))\n", "z = SpaceDimension(name='z', spacing=Constant(name='h_z', value=extent[1]/(shape[1]-1)))\n", "grid = Grid(extent=extent, shape=shape, dimensions=(x, z))" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Timestep size from Eq. 7 with V_p=6000. and dx=100\n", "t0, tn = 0., 200.\n", "dt = 1e2*(1. / np.sqrt(2.)) / 60.\n", "time_range = TimeAxis(start=t0, stop=tn, step=dt)\n", "\n", "src = DGaussSource(name='src', grid=grid, f0=0.01, time_range=time_range, a=0.004)\n", "src.coordinates.data[:] = [1000., 1000.]" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "scrolled": true }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3Xt4XVWd//H3N0lzT06Tpm3Spm16g1quxVCKgKIoIgpVRh28IOPgoL/f4Og48zjMzG8UHX0e787o+OgPgZGfl8G7FEUBuSgg0JYKtIXek9KmSZtb07TNPd/fH2eHnrZJmjTnnH0un9fznOecvc/K3t/snOSbtdZea5m7IyIiMlU5YQcgIiKZQQlFRETiQglFRETiQglFRETiQglFRETiQglFRETiQglFRETiQglFRETiQglFRETiIi/sAJKpqqrK6+rqwg5DRCStPPvss23uPvNU5bIqodTV1bF+/fqwwxARSStmtnsi5dTkJSIicaGEIiIicaGEIiIicaGEIiIicaGEIiIicaGEIiIicaGEIiIicZFV41BE4mnHgcM8tvUAC6tKWLmwkrLCaWGHJBIqJRSR07CpqYv33fEMXT0DAJw5u4x7b7mEwmm5IUcmEh41eYlM0q7Ww7zvjmcoLcjjdx+/jK//5Xls3d/N1x/aFnZoIqFSDUVkkr760DaGhp17bl7FvMpillWXs66xk9sf38Wbls+mvq4y7BBFQqEaisgkbN/fzf0bm7nxNQuYV1n8yv5/ufpVzCgp4PY/7goxOpFwKaGITMI3HtlB0bRcbrp00XH7SwvyWH3+HB7b2krX0YGQohMJlxKKyATt6TjKr1/Yxw0XL6CyJP+k91efP4f+oWF+t7k5hOhEwqeEIjJBa57fhzvcsGrBqO+fMzfCwqoS7n1uX5IjE0kNSigiE+Du/OrPTdQvqKC2onjUMmbGtefN4ald7ew/1JvkCEXCp4QiMgFbWrrZfuAwq8+fM265t51bgzs8suVAkiITSR1KKCITsOb5feTmGFefUzNuuSWzSqkqzWddY0eSIhNJHUooIqfg7tz3/D4uXVLFjNKCccuaGfULKlnf2Jmk6ERShxKKyClsP3CYvZ09XHV29YTK19dV8HLHUfWjSNZRQhE5hce2RvtDLj9z5oTKr1wYHSm/tkHNXpJdlFBETuGxra0sqy6jJlI0ofLLa8opzs9lvfpRJMsooYiM43DfIOsaO3jdBGsnAHm5OVwwv4K16keRLKOEIjKOJ3e0MTDkXH7GrEl9XX1dBVtaDnGoV9OwSPZQQhEZx2NbWyktyKO+rmJSX3de7XTcYUtzd4IiE0k9Sigi43hyRxsXL57BtNzJ/aosqykDYGvLoUSEJZKSlFBExtB0sIeXO45y8aIZk/7a6vJCIkXTeKlFNRTJHkooImN4amc7ABcvnnxCMTOWVZexpVk1FMkeSigiY3hqZzsVxdM4c3bZaX39q2rK2drSzfCwxzkykdSkhCIyCnfn6V3trFo0g5wcO61jnFldxpH+IfZ29sQ5OpHUFGpCMbOrzGyrme0ws1tHeb/AzH4cvP+MmdWd8P58MztsZv+YrJglO+zp6KHpYM9pNXeNWFYdrdlsUce8ZInQEoqZ5QLfAt4CLAfeY2bLTyh2E9Dp7kuArwNfPOH9rwG/TXSskn2e2tUGcFod8iPOmF2GWXTqe5FsEGYNZSWww913uXs/cA+w+oQyq4G7g9c/A64wMwMws7cDDcDmJMUrWeSZXR1UleazZFbpaR+jpCCPBZXFqqFI1ggzocwF9sRs7w32jVrG3QeBLmCGmZUC/wR8JglxShZat7uD+gWVBP+/nLYzq8s0uFGyRrp2yt8GfN3dD5+qoJndbGbrzWx9a2tr4iOTtNfS1cuejp5Jj44fzdJZZezuOMrA0HAcIhNJbWEmlCZgXsx2bbBv1DJmlgdEgHbgIuBLZtYIfBz4FzO7ZbSTuPvt7l7v7vUzZ058gj/JXiOrLY5MQz8VdVUlDA07ezqOTvlYIqkuL8RzrwOWmtlCoonjeuC9J5RZA9wIPAW8E3jE3R24bKSAmd0GHHb3/0pG0JL51jd2UJyfy/Ka8ikfa2FVCQCN7UdYNPP0+2NE0kFoNZSgT+QW4AHgJeAn7r7ZzD5rZtcGxe4k2meyA/gEcNKtxSLxtraxkwvmV5A3yfm7RjOSUHa1HpnysURSXZg1FNz9fuD+E/Z9KuZ1L/CuUxzjtoQEJ1npUO8AW1oO8bErlsbleBXF04gUTaOxXQlFMl+6dsqLJMSzuztxhwvrpt5/AtE5vRZWldDQpoQimU8JRSTG+sYOcnOMFfOnx+2YC6tKaGxTp7xkPiUUkRjrGjo5e045xfnxaw1eWFVC08EeegeG4nZMkVSkhCIS6Bsc4rm9B+PW3DWiLuiY392uWopkNiUUkcCmpi76B4epj3NCWRQklIa2U47DFUlrSigigbUNnQBcGIcR8rHqXkkoqqFIZlNCEQmsb+xg0cwSZpQWxPW4pQV5zCwrYFeraiiS2ZRQRIDhYWf97k4uXBDf5q4RdTOKeVnTr0iGU0IRAbYd6KarZyAuE0KOpraiWCs3SsZTQhEB1jVEJ4RcNYUFtcYzr6KI5q4ezTosGU0JRQR4pqGDmkghtRVFCTl+bUUxwx6dGl8kUymhSNZzd9Y2dLBy4dQX1BrLSKLSNPaSyZRQJOvtbj/Kge6+uKx/MpZ5lcUA6keRjKaEIllvbdB/clECE0p1pJAcgz2dqqFI5lJCkaz3TEMHlSX5LE7gAljTcnOoiRSphiIZTQlFst7axnYurKtIWP/JiNqKIvaqhiIZTAlFslpzVw97OnpYuTAxtwvHmldZzJ4O1VAkcymhSFZLRv/JiNqKIvZ399I3qGnsJTMpoUhWW9vQQWlBHq+qKU/4ueZVFOMO+w5qLIpkJiUUyWprGzqor6sgNyex/SdwbCyK+lEkUymhSNZqP9zH9gOHEzr+JNbIWBT1o0imUkKRrLWuMbr+STL6TwBmlxeSl2M0HVQNRTKTEopkrbUNHRTk5XDO3OlJOV9ujjG7vJBm9aFIhlJCkaz11K52Xr2ggvy85P0azJleSNNBNXlJZlJCkazUdriPl5oPccmSqqSetyZSRLNmHJYMpYQiWempne0AyU8o0wtp6epleNiTel6RZFBCkaz05I42ygrzOGduJKnnnTu9iP6hYdqO9CX1vCLJoIQiWenJnW1cvGhGUsafxKqJRMeiqGNeMlGoCcXMrjKzrWa2w8xuHeX9AjP7cfD+M2ZWF+x/k5k9a2Ybg+c3JDt2SV8vtx9lT0dP0pu7AGoihQDsU8e8ZKDQEoqZ5QLfAt4CLAfeY2bLTyh2E9Dp7kuArwNfDPa3Ade4+znAjcD3kxO1ZIIndrQBye8/gWiTF8A+dcxLBgqzhrIS2OHuu9y9H7gHWH1CmdXA3cHrnwFXmJm5+5/dfV+wfzNQZGYFSYla0t4fth1gTqSQxTNLkn7u6cXTKJyWQ7NqKJKBwkwoc4E9Mdt7g32jlnH3QaALOHGe8b8ANri7ejnllPoHh3liexuXL5uV8PVPRmNmzIkUsa9LCUUyT17YAUyFmZ1FtBnsynHK3AzcDDB//vwkRSapan1jB0f6h3j9mbNCi2HO9CLNOCwZKcwaShMwL2a7Ntg3ahkzywMiQHuwXQv8EviAu+8c6yTufru717t7/cyZM+MYvqSjR7ceID83h9csTvyCWmOpiRTSrBqKZKAwE8o6YKmZLTSzfOB6YM0JZdYQ7XQHeCfwiLu7mU0HfgPc6u5PJi1iSXuPbm3lokWVlBSEVzmvmV7Ege4++geHQ4tBJBFCSyhBn8gtwAPAS8BP3H2zmX3WzK4Nit0JzDCzHcAngJFbi28BlgCfMrPngkd4bRiSFvZ0HGXHgcNcHmJzF8Dc6YW4w/5DavaSzBJqH4q73w/cf8K+T8W87gXeNcrXfQ74XMIDlIzy4Iv7AXjDsnATysjgxn0He15ZI0UkE2ikvGSNBza1sKy6jIVVyb9dONbI4MYW1VAkwyihSFZo7e5j3e4Orjq7OuxQqB5JKBrcKBlGCUWywoMvtuBOSiSUssJplBbkaRp7yThKKJIVfrephYVVJZw5uyzsUIBoLUU1FMk0SiiS8TqO9PPUznbefFZ1KKPjR1MTKVQfimQcJRTJePc9v4/BYeftK+aEHcorZperhiKZRwlFMt4vNuxleU05y6rLww7lFTWRQg509zI4pMGNkjmUUCSj7TjQzfN7u7jughPnHQ1XdaSQYYfWw5rTVDLHKROKmRWb2b+Z2XeD7aVm9rbEhyYydb/Y0ERujnHt+anT3AUxY1HU7CUZZCI1lP8G+oCLg+0mNEpd0sDA0DA/37CX1y6tYlZZYdjhHGd2uRKKZJ6JJJTF7v4lYADA3Y8CqXGrjMg4Htjcwv5Dfbx/1YKwQznJK2vLK6FIBplIQuk3syLAAcxsMdEai0hKu/tPjcyvLA59MsjRVBRPIz8vR7cOS0aZSEL5NPA7YJ6Z/RB4GPhkQqMSmaLN+7pY19jJDasWkJuTehVqM4uORVENRTLIKWcbdveHzGwDsIpoU9fH3L0t4ZGJTMGdTzRQOC2Hd9fPO3XhkFRrLIpkmDETipldcMKu5uB5vpnNd/cNiQtL5PQ1th3h3uf2cePFdUSKp4UdzpiqI4VseLkz7DBE4ma8GspXg+dCoB54nmgN5VxgPcfu+hJJKd98ZAd5OcZHLl8Udijjqo4Usr+rj+FhJycFm+VEJmvMPhR3f727v55ozeSCYF32VwMrOHntd5GU0NB2hF8918T7Vy1IuVuFT1RTXkj/0DAdR/vDDkUkLibSKX+mu28c2XD3TcCrEheSyOn7/G9eoiAvhw+/LrVrJwDVwa3D6keRTDGRhPKCmd1hZpcHj+8CLyQ6MJHJevil/fz+pf187IqlKV87AS20JZlnImvKfxD4X8DHgu0/At9OWEQip6F3YIjP3PciS2aV8sFLFoYdzoSMTL/SrLEokiEmcttwL/D14CGSkj73mxd5ueMoP/qbi8jPS485T6tKC8jNMVq6esIORSQuTplQzKyBYJR8LHdP/UZqyQq/3djMD55+mQ+/dhGvWVwVdjgTlptjzC4roKVLE09IZphIk1d9zOtC4F1AZWLCEZmcLS2H+OTPX+C8edP5hyvPDDucSauOFNJySDUUyQynbBtw9/aYR5O7/wfw1iTEJjKuPR1H+cCdaynOz+Vb712RNk1dsaojhZogUjLGRJq8YkfM5xCtsUykZiOSMDtbD3PjXWvpHRjipx95DbUVxWGHdFqqy4t4bGsr7p4y692LnK6JJIavxrweBBqAdycmHJFT+9PONv73DzeQl2N8/6aLOLO6LOyQTltNpJCj/UN09w1SXpi608SITMREEspN7r4rdoeZpcd9mZJRevqH+OqDW7nzyQYWVpXwvb9ayfwZ6VkzGRE7FkUJRdLdRBqdfzbBfZNmZleZ2VYz22Fmt47yfoGZ/Th4/xkzq4t575+D/VvN7M3xiEdSU0//EN9/qpHXfvlR7niigfdftIBff/TStE8mcCyhqB9FMsF4sw0vA84CImZ2Xcxb5UTv9poSM8sFvgW8CdgLrDOzNe7+Ykyxm4BOd19iZtcDXwT+0syWA9cH8c0Bfm9mZ7j70FTjktTQNzjE2oYOHtjcwr3P7aO7d5AL6yr49vsuoL4uc24yrH5lKWDd6SXpb7wmrzOBtwHTgWti9ncDfxOHc68Edow0p5nZPcBqIDahrAZuC17/DPgvi/Zcrgbucfc+oMHMdgTHeyoOcUmSDA87HUf7aenq5UB3L3s7e9jcdIhN+7rYtr+bgSGncFoOVy6v5n0XzWflwsqM67geWVteNRTJBGMmFHe/F7jXzC5290T8oZ4L7InZ3gtcNFYZdx80sy5gRrD/6RO+dm4CYgTgtjWbaQ7+g/SYIZ6xoz39uKGfPur+48v7GPtPXf64Mx1XfozzTqTMcQedXPkJfS/u9AwMcbR/5DFI78DwSd/P9OJpnDM3wk2XLqJ+QQWXLKmiKD/3pHKZIj8vh6rSAvZn6fQrfYNDtHb3vfLoPNrP0f4hegaG6A2eewaGGBqOfoaG3Rl2GHbHg+dj287wyR+p0PjJ48FD9c33XJDwW+vHa/L6pLt/CXivmb3nxPfd/e8SGlmcmNnNwM0A8+fPP61jNB3sYU/H0bGOf+z1cfvHeB1T6vj9o3/xmMecbAyM/sUW82wxGxbzzmS/l1gjsRlQMy2X4vxcivJHnvOYUZLP7PICZpcXUhMpYnZ5QcbVQk6lJkvGovQPDrPh5U6e3d3JS82H2NLSTUPbEYaGx/7DWzQtl8JpOeTl5pBjkGNGjkU/k9HXHLdtdsJnPWSp9FFORoIbr8nrpeB5fYLO3QTErs9ay8nrrIyU2WtmeUAEaJ/g1wLg7rcDtwPU19ef1hX97gfqT11I5DTNLi9kb+fo/7CkO3fnmYYOfrp+Lw9sbuFw3yAAtRVFLKsu56qzqplXWURVaQEzywqoLMmnJD+PovxcCvJysu6fi3Q3XpPXfcHz3Qk69zpgaXALchPRTvb3nlBmDXAj0b6RdwKPuLub2RrgR2b2NaKd8kuBtQmKUyShaiKFrGvsCDuMuPvjtla+9tA2nttzkNKCPK4+p5o3LJvNxYtnECnSLdKZaLwmr/sYZVLIEe5+7VROHPSJ3AI8AOQCd7n7ZjP7LLDe3dcAdwLfDzrdO4gmHYJyPyHagT8I/K3u8JJ0VR0ppKtngJ7+oYzoL2rt7uMz923m1y80M3d6EZ9/x9lct6I2I743Gd94TV5fSfTJ3f1+4P4T9n0q5nUv0ckoR/vazwOfT2iAIkkwsi5Ky6FeFlaVhBzN1Gx4uZMPf/9Zuo4O8A9vOoMPv25xWs6xJqdnvCavP4y8NrN8YBnRGstWd9ci2CJxUv3KrcM9aZ1QfvNCM3//4+eYHSng+x+9hGXV5WGHJEk2kckh3wp8B9hJ9GadhWb2YXf/baKDE8kGmbAU8O82tfB39/yZC+ZP5/Yb6qkoyQ87JAnBRCeHfL277wAws8XAbwAlFJE4SPfpV57a2c5H/2cD59ZG+O8PrqS0QJORZ6uJNG52jySTwC6io+VFJA6K8/OIFE1Ly8GNzV093PKjDcyvLOZ7SiZZbyI//fVmdj/wE6J9KO8iOu/WdQDu/osExieSFdJxcOPA0DB/+8MN9A4M8X9vWKVbgWVCCaUQ2A+8LthuBYqIzu/lgBKKyBTNLi9Muz6UOx5vYMPLB/nme1awZFb6rkkj8XPKhOLuH0xGICLZrCZSyOZ9h8IOY8Ia247wH7/fxlVnVXPNeXPCDkdSxETu8loIfBSoiy0/1YGNInJMdaSQ9iN99A8Op/y4DXfn3+7dRH5uDp9ZfVbY4UgKmUiT16+Ijli/D0ihuTxFMkdNpBB3ONDdS21Fai8c9sftbTy+vY1PX7P8len3RWBiCaXX3b+R8EhEstjs8mNjUVI5oQwPO19+YAu1FUW876IFYYcjKWYiCeU/zezTwINA38hOd9+QsKhEskxNpAhI/bEov93UwqamQ3zt3eelfNOcJN9EEso5wA3AGzjW5OXBtojEwcjgxlQei+LufPOR7SydVcrq8xO2np2ksYkklHcBizR/l0jilBfmUZyfm9I1lCd2tLGlpZsvv/NccnO0TomcbCJ11k1E15UXkQQxM6ojqT0W5Y7HG5hZVsC15+s2YRndRGoo04EtZraOY30o7u6rExeWSPapLi+kuasn7DBGtW1/N3/Y1so/XnkGBXla10RGN5GE8umY1wZcRrDQlYjET3WkkKd3tocdxqh+8PRu8vNyeK/u7JJxnLLJK1gX5RDwNuB7RDvjv5PYsESyT02kkAPdfQwNj7lQaih6B4b41Z+beMvZ1VRqWnoZx3hLAJ8BvCd4tAE/BszdX5+k2ESySnWkiMFhp/1wH7NSaMDggy/u51DvIO+unxd2KJLixquhbCFaG3mbu1/q7t8EtG67SILUlKfmuig/Xb+HudOLuHjRjLBDkRQ3XkK5DmgGHjWz75rZFUT7UEQkAVJxoa2mgz08saONd9XXkqNbheUUxkwo7v4rd7+e6FryjwIfB2aZ2bfN7MpkBSiSLVJxcON9z+/DHa5bURt2KJIGJtIpf8Tdf+Tu1wC1wJ+Bf0p4ZCJZprI4n/zcnJSqody/sZnzaiPMn5G684tJ6pjUZDzu3unut7v7FYkKSCRb5eQYsyMFtKTIWJQ9HUd5YW8XV59TE3YokiY0u5tICokObkyNGspvNjYDKKHIhCmhiKSQ6khRyvShjDR3zatUc5dMjBKKSAqpiURrKO7hDm5sOtij5i6ZNCUUkRRSXV5I3+AwB48OhBrHIy/tB+CNy2eHGoekFyUUkRRSkyJjUR7ecoCFVSUsnlkaahySXkJJKGZWaWYPmdn24LlijHI3BmW2m9mNwb5iM/uNmW0xs81m9oXkRi+SOLODhNJyKLw7vY72D/Knne28Ydms0GKQ9BRWDeVW4GF3Xwo8HGwfx8wqic50fBGwEvh0TOL5irsvA1YAl5jZW5ITtkhijdRQWrr6TlEycZ7c0U7/4DBXKKHIJIWVUFYDdwev7wbePkqZNwMPuXuHu3cCDwFXuftRd38UIFhFcgPRAZciaW9maQE5RqhjUR5+aT9lBXnU11WGFoOkp7ASymx3bw5etwCj9fzNBfbEbO8N9r3CzKYD1xCt5YikvbzcHGaVhTcWxd15bGsrl51RRX6eulhlciaywNZpMbPfA9WjvPWvsRvu7mY26XskzSwP+B/gG+6+a5xyNwM3A8yfP3+ypxFJutmRQlpCGouy48BhWg718tqlM0M5v6S3hCUUd3/jWO+Z2X4zq3H3ZjOrAQ6MUqwJuDxmuxZ4LGb7dmC7u//HKeK4PShLfX19aq1cJDKKmvJCdrYeDuXcf9zeBsClS6tCOb+kt7DqtGuAG4PXNwL3jlLmAeBKM6sIOuOvDPZhZp8DIkRnQBbJKNWRQlpCavJ6Ynsri6pKqK3Q6HiZvLASyheAN5nZduCNwTZmVm9mdwC4ewfw78C64PFZd+8ws1qizWbLgQ1m9pyZfSiMb0IkEWoihXT3DdLdm9zBjX2DQzy9q0O1EzltCWvyGo+7twMnzVjs7uuBD8Vs3wXcdUKZvWihL8lgseuilBVOS9p5N+w+SM/AEJcuUUKR06PbOERSTHV5OGNRntjRSm6OsWqxlvqV06OEIpJiaiJFADQneSzKn3a2c25thPIk1ooksyihiKSYWeUFAEntmD/SN8gLe7u4eJFqJ3L6lFBEUkzhtFxmlOTTnMSxKOt3dzI07Fys5i6ZAiUUkRQ0uzy5tw4/vaudvBzj1QtGnadVZEKUUERSUE2Sx6I8tbOd8+ZNpzg/lBs/JUMooYikoOokTr9yuG+QjU3qP5GpU0IRSUE1kUI6jvTTOzCU8HOtb+xgaNhZpYQiU6SEIpKCqoNbh/cnoZbyTEMHeTnGBQumJ/xcktmUUERS0MhCW00HEz8WZX1jB2fPjaj/RKZMCUUkBc2vjE7OuLcjsQmld2CI5/d0sXKhFtOSqVNCEUlBNZFCcnOMlzuOJvQ8G5u66B8apl63C0scKKGIpKC83BzmTC9kT2diE8rahg4ALfcrcaGEIpKi5lcWJ7yGsr6xgyWzSqksyU/oeSQ7KKGIpKh5FcXsSWBCGRp21u/u5ELVTiROlFBEUtS8ymLaDvdzpG8wIcfftr+b7t5BLqxT/4nEhxKKSIp65U6vzsTc6bWuMdp/ohqKxIsSikiKGkkoiepHWdfYSXV5IbUVRQk5vmQfJRSRFDUvgQnF3VnX0EF9XQVmWlFb4kMJRSRFVRRPo7QgLyEd83s7e2g51KsBjRJXSigiKcrMmFeZmDu9RvpP6hcooUj8KKGIpLB5FUUJafJa19hJWWEeZ1aXxf3Ykr2UUERS2MjgxuFhj+tx1zV28OoFFeTmqP9E4kcJRSSF1VWV0Dc4HNfFtjqO9LPjwGHdLixxp4QiksIWzSwBYFfrkbgdc2T+rovUIS9xpoQiksKWzCwFYGfr4bgdc11jBwV5OZxTG4nbMUVACUUkpc0sK6C0II9dcUwoaxs6WDF/OgV5uXE7pgiElFDMrNLMHjKz7cHzqJMJmdmNQZntZnbjKO+vMbNNiY9YJBxmxuKZJeyMU5NXd+8Am/d1sXKh1o+X+AurhnIr8LC7LwUeDraPY2aVwKeBi4CVwKdjE4+ZXQfE7982kRS1aGZp3Gooz+7uZNjVfyKJEVZCWQ3cHby+G3j7KGXeDDzk7h3u3gk8BFwFYGalwCeAzyUhVpFQLZ5Zwr6uXo72T33W4bUNHeTlGCvmT49DZCLHCyuhzHb35uB1CzB7lDJzgT0x23uDfQD/DnwVSOzqQyIpYFHQMR+PO73WNnRwTm2E4vy8KR9L5EQJSyhm9nsz2zTKY3VsOXd3YMKjtszsfGCxu/9yguVvNrP1Zra+tbV1ct+ESApYHKc7vXoHhnh+70HN3yUJk7B/U9z9jWO9Z2b7zazG3ZvNrAY4MEqxJuDymO1a4DHgYqDezBqJxj/LzB5z98sZhbvfDtwOUF9fH9/hxiJJsGBGMTk29RrKn18+yMCQq/9EEiasJq81wMhdWzcC945S5gHgSjOrCDrjrwQecPdvu/scd68DLgW2jZVMRDJB4bRcaiuKp1xDWdvQgRm8WhNCSoKElVC+ALzJzLYDbwy2MbN6M7sDwN07iPaVrAsenw32iWSdxTNL2HFgigmlsZ1XVZcTKZoWp6hEjhdKz5y7twNXjLJ/PfChmO27gLvGOU4jcHYCQhRJKcvnlPP49jb6BodOa0Bi/+Awz+7u5PoL5ycgOpEojZQXSQNnzYkwOOxsazm9WsqmfV30Dgyr/0QSSglFJA2cNacciCaG0zEyIeSFSiiSQEooImlgfmUxZYV5bD7NhPLkjjaWzCqlqrQgzpGJHKOEIpIGzIzlNeVsajo06a/tHRhibUMHr106MwGRiRyjhCKSJs6aE2FLyyGGJrl647rGDvoGh7nsjKoERSYSpYQikibOmlNO78DwpCeK/OO2VvJzc9QhLwnOLpMkAAAJP0lEQVSnhCKSJs6eG10Qa7Id849vb6O+rkLzd0nCKaGIpInFM0soyMth496J96McONTLlpZuLlP/iSSBEopImsjLzeH8edN5pqF9wl/z2LbohKiXLVX/iSSeEopIGrlkSRUvNh+i80j/hMo/sKmFudOLXhnHIpJISigiaeSSJTNwh6d2nbqW0t07wOPb27jq7GrMLAnRSbZTQhFJI+fWTqckP5c/7Ww7ZdlHthygf2iYq8+pTkJkIkooImllWm4OFy2awZ92nLqGcv/GZmaXF7BiXkUSIhNRQhFJO69ZPINdbUdo7uoZs8yRvkH+sK2Vq86qJidHzV2SHEooImnmkiXRO7Ye3TL2kta/eq6J3oFhrjlvTrLCElFCEUk3y6rLOHN2GT9auxv3k6dhGR527nqigXPmRnj1AjV3SfIooYikGTPj/avms6npEM/vPXnU/B+3t7Kz9Qh/fWmd7u6SpFJCEUlDb18xl+L8XH7w9O6T3rvziQZmlhXw1nPU3CXJpYQikobKCqfxjhVzue/5fTQdPNY5/7Nn9/L49jZuunQh+Xn69Zbk0idOJE196LJF5OfmcMOdz9B+uI9NTV38n19tZNWiSj506cKww5MspOlHRdLUwqoS7vyrC7nhzme4+AuP0D84TFVpAd+4fgV5ufpfUZJPCUUkja1cWMn3PriS+17Yx5KZpbxp+WxmlReGHZZkKSUUkTR38eIZXLx4RthhiKgPRURE4kMJRURE4kIJRURE4kIJRURE4kIJRURE4kIJRURE4kIJRURE4kIJRURE4sJGW08hU5lZK3Dy9KwTUwWceiHv5FNck6O4JkdxTU6mxrXA3WeeqlBWJZSpMLP17l4fdhwnUlyTo7gmR3FNTrbHpSYvERGJCyUUERGJCyWUibs97ADGoLgmR3FNjuKanKyOS30oIiISF6qhiIhIXCihnIKZXWVmW81sh5ndGmIc88zsUTN70cw2m9nHgv23mVmTmT0XPK4OIbZGM9sYnH99sK/SzB4ys+3Bc0WSYzoz5po8Z2aHzOzjYV0vM7vLzA6Y2aaYfaNeI4v6RvCZe8HMLkhiTF82sy3BeX9pZtOD/XVm1hNz3b6TiJjGiWvMn5uZ/XNwrbaa2ZuTHNePY2JqNLPngv3JvF5j/W1I/ufL3fUY4wHkAjuBRUA+8DywPKRYaoALgtdlwDZgOXAb8I8hX6dGoOqEfV8Cbg1e3wp8MeSfYwuwIKzrBbwWuADYdKprBFwN/BYwYBXwTBJjuhLIC15/MSamuthyIVyrUX9uwe/A80ABsDD4fc1NVlwnvP9V4FMhXK+x/jYk/fOlGsr4VgI73H2Xu/cD9wCrwwjE3ZvdfUPwuht4CZgbRiwTtBq4O3h9N/D2EGO5Atjp7qc7qHXK3P2PQMcJu8e6RquB/+dRTwPTzawmGTG5+4PuPhhsPg3Uxvu8pxPXOFYD97h7n7s3ADuI/t4mNS4zM+DdwP8k4tzjGedvQ9I/X0oo45sL7InZ3ksK/BE3szpgBfBMsOuWoOp6V7KblgIOPGhmz5rZzcG+2e7eHLxuAWaHENeI6zn+Fz3s6zVirGuUKp+7vyb6n+yIhWb2ZzP7g5ldFkI8o/3cUuVaXQbsd/ftMfuSfr1O+NuQ9M+XEkqaMbNS4OfAx939EPBtYDFwPtBMtNqdbJe6+wXAW4C/NbPXxr7p0Xp2KLcTmlk+cC3w02BXKlyvk4R5jUZjZv8KDAI/DHY1A/PdfQXwCeBHZlaexJBS8ucW4z0c/09L0q/XKH8bXpGsz5cSyviagHkx27XBvlCY2TSiH5gfuvsvANx9v7sPufsw8F0SVN0fj7s3Bc8HgF8GMewfqUYHzweSHVfgLcAGd98fxBj69Yox1jUK9XNnZn8FvA14X/CHiKBJqT14/SzRvoozkhXTOD+30H9HzSwPuA748ci+ZF+v0f42EMLnSwllfOuApWa2MPhP93pgTRiBBG20dwIvufvXYvbHtn2+A9h04tcmOK4SMysbeU20U3cT0et0Y1DsRuDeZMYV47j/HMO+XicY6xqtAT4Q3I2zCuiKabpIKDO7CvgkcK27H43ZP9PMcoPXi4ClwK5kxBScc6yf2xrgejMrMLOFQVxrkxVX4I3AFnffO7IjmddrrL8NhPH5SsZdCOn8IHpHxDai/2H8a4hxXEq0yvoC8FzwuBr4PrAx2L8GqElyXIuI3mXzPLB55BoBM4CHge3A74HKEK5ZCdAORGL2hXK9iCa1ZmCAaJv1TWNdI6J333wr+MxtBOqTGNMOou3rI5+x7wRl/yL4+T4HbACuSfK1GvPnBvxrcK22Am9JZlzB/u8BHzmhbDKv11h/G5L++dJIeRERiQs1eYmISFwooYiISFwooYiISFwooYiISFwooYiISFwooYiISFzkhR2ASKozs5H7+QGqgSGgNdg+6u6vScA5VwC3uPtNUzzOLURjvCs+kYmMTeNQRCbBzG4DDrv7VxJ8np8Cn3P356d4nGLgSY/OKSWSUGryEpkCMzscPF8ezCp7r5ntMrMvmNn7zGytRRcfWxyUm2lmPzezdcHjklGOWQacO5JMLLq41N1m9riZ7Taz68zsS8FxfxfM40RwzheDGXm/AuDR6VMazSzMOcskSyihiMTPecBHgFcBNwBnuPtK4A7go0GZ/wS+7u4XEp2e445RjlPPyXOMLQbeQHTm5B8Aj7r7OUAP8NagWe4dwFnufi7wuZivXU90enWRhFIfikj8rPNgkj0z2wk8GOzfCLw+eP1GYHl0Pj8Ays2s1N0PxxynhmN9NCN+6+4DZraR6AqUv4s5dh3wa6AXuNPMfh1sjzgALJvi9yZySkooIvHTF/N6OGZ7mGO/aznAKnfvHec4PUDhaMd292EzG/BjnZ/DRJfsHQyata4A3gncQrRGQ3CsntP4fkQmRU1eIsn1IMeavzCz80cp8xKwZDIHDRZXirj7/cDfE21+G3EG4U7TL1lCCUUkuf4OqA86zl8k2udyHHffAkRG1pmZoDLg12b2AvAE0VUCR1wCPDSFmEUmRLcNi6QgM/t7oNvdR+u0n8xxVgCfcPcb4hOZyNhUQxFJTd/m+D6Z01UF/FscjiNySqqhiIhIXKiGIiIicaGEIiIicaGEIiIicaGEIiIicaGEIiIicfH/AYCRK7Ane6haAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "#NBVAL_SKIP\n", "\n", "src.show()" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "scrolled": true }, "outputs": [], "source": [ "# Now we create the velocity and pressure fields\n", "p = TimeFunction(name='p', grid=grid, staggered=NODE, space_order=2, time_order=1)\n", "v = VectorTimeFunction(name='v', grid=grid, space_order=2, time_order=1)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "scrolled": true }, "outputs": [], "source": [ "from devito.finite_differences.operators import div, grad\n", "t = grid.stepping_dim\n", "time = grid.time_dim\n", "\n", "# We need some initial conditions\n", "V_p = 4.0\n", "density = 1.\n", "\n", "ro = 1/density * dt\n", "l2m = V_p*V_p*density * dt\n", "\n", "# The source injection term\n", "src_p = src.inject(field=p.forward, expr=src)\n", "\n", "# 2nd order acoustic according to fdelmoc\n", "u_v_2 = Eq(v.forward, v + ro * grad(p))\n", "u_p_2 = Eq(p.forward, p + l2m * div(v.forward))" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle \\left[\\begin{matrix}\\operatorname{v_{x}}{\\left(t + dt,x + \\frac{h_{x}}{2},z \\right)}\\\\\\operatorname{v_{z}}{\\left(t + dt,x,z + \\frac{h_{z}}{2} \\right)}\\end{matrix}\\right] = \\left[\\begin{matrix}\\operatorname{v_{x}}{\\left(t,x + \\frac{h_{x}}{2},z \\right)} + 1.17851130197758 \\frac{\\partial}{\\partial x} p{\\left(t,x,z \\right)}\\\\\\operatorname{v_{z}}{\\left(t,x,z + \\frac{h_{z}}{2} \\right)} + 1.17851130197758 \\frac{\\partial}{\\partial z} p{\\left(t,x,z \\right)}\\end{matrix}\\right]$" ], "text/plain": [ "Eq(Vector(v_x(t + dt, x + h_x/2, z), v_z(t + dt, x, z + h_z/2)), Vector(v_x(t, x + h_x/2, z) + 1.17851130197758*Derivative(p(t, x, z), x), v_z(t, x, z + h_z/2) + 1.17851130197758*Derivative(p(t, x, z), z)))" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "u_v_2" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle p{\\left(t + dt,x,z \\right)} = p{\\left(t,x,z \\right)} + 18.8561808316413 \\frac{\\partial}{\\partial x} \\operatorname{v_{x}}{\\left(t + dt,x + \\frac{h_{x}}{2},z \\right)} + 18.8561808316413 \\frac{\\partial}{\\partial z} \\operatorname{v_{z}}{\\left(t + dt,x,z + \\frac{h_{z}}{2} \\right)}$" ], "text/plain": [ "Eq(p(t + dt, x, z), p(t, x, z) + 18.8561808316413*Derivative(v_x(t + dt, x + h_x/2, z), x) + 18.8561808316413*Derivative(v_z(t + dt, x, z + h_z/2), z))" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "u_p_2" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "op_2 = Operator([u_v_2, u_p_2] + src_p)" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "Operator `Kernel` run in 0.04 s\n" ] } ], "source": [ "#NBVAL_IGNORE_OUTPUT\n", "\n", "# Propagate the source\n", "op_2(time=src.time_range.num-1)" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAATUAAAD8CAYAAAAWjzPMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztnXusXdV95z/f2LwD5mGCwfbEDnZAhgIpFqVKVaUhAaeN4sxMHmaalrakqCqZptOOktBIoSWDRNRJ03SSZmoRN+TROAwJipUhdR0I6lRqwNflaYPjyyvYNQZjMITwiOE3f+y1zv3d47V99vE99/refX4f6cprr/1ae5/j3/n91u+xZGYEQRC0hdcd6gEEQRAMkhBqQRC0ihBqQRC0ihBqQRC0ihBqQRC0ihBqQRC0ihBqQRC0igkJNUkrJG2VNCrpE4MaVBAEwcGigw2+lTQL+DHwTmA7sBG41My2DG54QRAE/TF7AudeAIya2cMAktYCK4FaoXbUUUfZnDlzJnDLIAgOxN69e3nxxReVt1esWGG7d+9udO6mTZvWm9mKSRvcFDERoTYfeNxtbwd+6UAnzJkzhw996EMTuGUQBAfi61//+rjt3bt3MzIy0uhcSXMnY0xTzUSEWiMkXQFcAXDsscdO9u2CIOhi2PK7JyLUdgAL3faC1DcOM1sNrAaYN2/ecL3dIJgGvPbaa4d6CFPKRLyfG4GlkhZLOhxYBawbzLCCIBgEZtb4rxe9oh0kHSHpW2n/HZIWuX1Xpf6tki7pOm+WpLskfW8Aj3zwQs3M9gEfAdYDDwA3mtnmQQwqCILBMQihlqIdvgi8C1gGXCppWddhlwPPmNkS4HPAZ9K5y6iUnrOAFcDfputlPkolQwbChOLUzOwWM3uzmZ1uZtcOalBBEAyOAWlqnWgHM3sFyNEOnpXADal9E3CRJKX+tWb2spk9Aoym6yFpAfAbwPUDeVgioyAIWk8fQm2upBH3d4W7TCnaYX7XrTrHJEtuL3BSj3P/GvgYMLCJv0n3fgZBcGjpw/u528yWT+ZYPJLeDTxpZpskvW1Q1w2hFgQtxswG5f1sEu2Qj9kuaTYwB3j6AOe+B3iPpF8HjgSOk/R1M5tQMGuYn0HQcgY0p9Yk2mEdcFlqvw+4zaoLrwNWJe/oYmApcKeZXWVmC8xsUbrebRMVaBCaWhC0nkEE35rZPkk52mEWsMbMNku6Bhgxs3XAl4GvSRoF9lAJKtJxN1KlUO4DrjSzVyc8qBpCqAVByxlURoGZ3QLc0tX3Kdd+CXh/zbnXArUREmZ2O3D7IMYZQi0IWkzTwNo2EUItCFrOsKVJhVALgpYTmloQBK0hzM8gCFpHCLUgCFpFCLUgCFpFCLUgCFrDANOkZgwh1IKg5YSmFgRBqwihFgRBqwihFgRBqxg2odaz9JCkNZKelHS/6ztR0gZJ29K/J0zuMIMgOBiyo6DJX1toUk/tK1SLJXg+AdxqZkuBW9N2EATTkEGtJjVT6CnUzOyfqWojefwCCzcA7x3wuIIgGBDDJtQOdk7tFDPbmdpPAKcMaDxBEAyYNgmsJkzYUWBmJqn2raUVaa4AOPbYYyd6uyAI+qBtWlgTDnaNgl2STgVI/z5Zd6CZrTaz5Wa2/Oijjz7I2wVBcLAMm/l5sELNL7BwGfDdwQwnCIJBE97PLiR9E/hX4AxJ2yVdDlwHvFPSNuAdaTsIgmnIoDQ1SSskbZU0Kmm/iIe0WtS30v47JC1y+65K/VslXZL6jpR0p6R7JG2W9BeDeN6ec2pmdmnNrosGMYAgCCaPQZmWkmYBXwTeSbXC+kZJ68xsizvscuAZM1siaRXwGeCDkpZRrSx1FnAa8ANJbwZeBt5uZj+VdBjwL5K+b2Y/mshYY93PIGg5A9LULgBGzexhM3sFWEsV2uXxoV43ARdJUupfa2Yvm9kjwChwgVX8NB1/WPqbsAQOoRYELWdAQm0+8Ljb3p76iseY2T5gL3DSgc6VNEvS3VTOxg1mdsdBPmaHEGpB0HL6EGpzJY24vyumYGyvmtl5wALgAklnT/SakdAeBC2mzyKRu81sec2+HcBCt70g9ZWO2S5pNjAHeLrJuWb2rKQfUqVk3s8ECE0tCFrOgMzPjcBSSYslHU418b+u6xgf6vU+4DarLrwOWJW8o4uBpcCdkk6WdDyApKOonBAPTvR5Q1MLgpYzCO+nme2T9BFgPTALWGNmmyVdA4yY2Trgy8DXJI1S5YuvSudulnQjsAXYB1xpZq+mwP0bkmf1dcCNZva9iY41hFoQtJxBZQuY2S3ALV19n3Ltl4D315x7LXBtV9+9wFsGMjhHCLUgaDltSoFqQgi1IGgxfToKWkEItSBoOaGpBUHQKkKoBUHQKkKoBUFDes3VvO51EQZ5qGlbrbQmhFALgpYTQi0IujhY71loctOD8H4GQdAqQlMLgqA1xJxaMGOoau8NHv8fYLLNltL1wyQdPMMm1JqsUbBQ0g8lbUl1xD+a+k+UtEHStvTvCZM/3CAI+iVWk9qffcCfmtky4ELgylRz/BPArWa2FLg1bQdBMM0YNqHWZOGVncDO1H5e0gNUpXhXAm9Lh90A3A58fFJGOeSUTE3fV2eKlvpLX17f1+vL3c+xBzO+OpM3zNKDI3I/e5CWvHoLcAdwShJ4AE8Ap9ScEyu0B8EhpE1aWBMa//xJej3wbeCPzew5vy9Vtyy+OYsV2oPgkBLmZ4G0Jt+3gW+Y2XdS9y5Jp5rZzlTB8snJGuQw0cvUzGZYnflZMtPqzI/8Rfb7e5ms/tjc7vUfoslYc3+dqZrvFWZo/7RJYDWhifdTVGV6HzCzv3K7fD3yy4DvDn54QRBMlNDU9uetwG8B96X1+QD+DLgOuFHS5cBjwAcmZ4jtp6TJeI3Et2fNmnXA/f5aJU3Mt1999dX9zilRp6nl831fL0dCr+fK/3aPK7fDkdAf4SgoYGb/AtR96y8a7HCCIBg0g9LCJK0APk+18Mr1ZnZd1/4jgK8C51MtjfdBM3s07bsKuBx4FfgjM1svaWE6/hSqOfnVZvb5iY4zft6CoOUMwvxMKz59EXgXsAy4NMWrei4HnjGzJcDngM+kc5dRrSx1FtW6nn+brlcXAzshIk3qEFE3YZ7bs2ePfTS+nc0z3+fxX859+/bt11dq9zJP6saazU+P7yuZpx5vaubn6fXc4UjonwFpahcAo2b2MICktVSxqlvcMSuBP0/tm4AvpDn5lcBaM3sZeCQtoXeBmf0r5RhYf82+iW9AELSYplpaA8E3H3jcbW9PfcVjzGwfsBc4qcm5XTGwEyI0tSBoOX1oanMljbjt1Wa2ehKGNI4DxcAeDCHUppiSd9ObYYcffvi4f7vb/tjMz3/+8077lVde2a/t+/yxvchjrDPpDjvsMGDMzK27v9/v2/4/W77HEUcc0ekrvQP//L1i8sIUrejD+7nbzJbX7NsBLHTbC1Jf6ZjtkmYDc6gcBrXn1sTAToj41IOg5QzI/NwILJW0WNLhVBP/67qO8bGr7wNuS9lG64BVko6QtBhYCtx5gBjYCRGa2hTgJ7dLE/1eIznqqKOA8RqLPzZPvr/44oudvp/+9KfFdj6mTjsrOSX8ffNYs0bW3c77S84Jf38/pueff36//f65/FiPPPLITjtrG3XvJbS2MoMKrDWzfZI+AqynCulYY2abJV0DjJjZOioB9bXkCNhDJfhIx91I5QDYB1xpZq9K+hUKMbBmdstExhpCLQhazqDi1JKwuaWr71Ou/RLw/ppzrwWu7eo7UAzsQRNCLQhaTptSoJoQQm2SKJmcMGYyedPKVy/J5qc38/zkezbfnnrqqU7f008/3Wl7Uy+bgiVHhB9DKbXK93szz487t/01vZmXTckXXnihOFb/DM8++yww3iT1pmg2T/34So6EYTUzD0QItSAIWkPkfgZB0DpCUwsmRCkOrWS+HXPMMZ2+17/+9Z12NqNeeumlTp830x577DEAduwYCxF67rlyvGI2ZX3FYW+KZhPXm7reFM2/8HVxaPm5vEfSP1e+v7/mvHnzOm3/DI8++igAu3bt6vSVvKNe6/DvreTJ9Z/BMKdRhVALgqBVhFAL+qaU8F03OZ+dAl578pPvWUP793//907fli1j+b2jo6PAeIeAP//kk0/er3388cd3+rL2BGMaWp2jIGtoXlP72c9+1mnnMfhr+nvNnTu3dkzd7fw+vCb1k5/8pNPOzoaSI8O366rsluqxDYvWFkItCILWEI6CIAhaR2hqQd+UzM9S6hOMTW772DQfj5XNznvuuafT59vZ5MumHcCiRYs67YULx/KGsynozVNPKTWplHDu+0qxY57du3d32jkmzZusp59++gHH6p0W/vrZQeKvVYoFrCsHXqo/Nyym6LAJtSYLrxwp6U5J90jaLOkvUv9iSXdIGpX0rZTkGgTBNGPYFl5p8vP0MvB2MzsXOA9YIelCqlK9n0ule5+hKuUbBME0YoBFImcMTRZeMSC72g5Lfwa8Hfgvqf8GqjK+Xxr8EKcndSWue8VulRZ09nFo999/PwCbNm0q7s8m27nnntvpW7p0aad93HHH7TfGUo0zGDMlvZlX8i56k6705X/55Zc7bR9btnfvXmB8apRPmTr//PP3e6669LB83hNPPFG8Vj6vLk6tV224NtMmgdWEposZzwI2AUuoFl94CHg2leyFcmnffO4VwBUwPowhCIKpIbyfBczsVeA8SccDNwNnNr1BKge8GmDevHkz+iejyUrjWWOoS1jPmkTWYmAs9gzgrrvuAmD79u2dPh+Fv3x5VZjUa2reEVHSmnx2gtfKeq2wXtJqvAOklJFQur93Hnjtyt83O1AWLFjQ6TvvvPM67T179gDj4/P8O8wOhNL4oLxwSyl7AtqlzbXNtGxCX95PM3tW0g+BXwaOlzQ7aWul0r5BEEwDhk2oNfF+npw0NCQdBbwTeAD4IVXJXqhK+H53sgYZBMHBE46C/TkVuCHNq70OuNHMvidpC7BW0v8A7qIq5Ts09Fo4xZuE3iTKMV87d+7s9G3evLnTfvjhh4HxE97nnHNOp51NMu8QKJlhUHYElKhLk2qKH2uv2LFcNw3g7rvv7rRPPPFEYHz83Zve9KZOOzsI/Hvzpmg2e72p7Z012RSti2NrM20SWE1o4v28l2o9vu7+h6kWOA2CYJoyjGlS7ZkRDYKgyKDMT0krJG1NAfefKOw/IgXij6bA/EVu31Wpf6ukS1z/GklPSrp/QI8baVL90Cs2DcZMHm/6+GOzqZjTfmC89zP/qi5evLjT572AubLFM8880+nzHsUS/v7+yztRsySfX3ed/L6899fHnnmvaDZFfRrVL/zCL3TaORbvoYce6vT52mv5fXjz008B5M+jzlQueULb4gUdhPmZpp++SDWnvh3YKGmdmW1xh10OPGNmSyStogrQ/6CkZVQrS50FnAb8QNKbU1TFV4AvAF+d8CAT7fjUgiCoZUCa2gXAqJk9bGavAGuBlV3HrKQKxAe4Cbgore25ElhrZi+b2SPAaLoeZvbPVMvpDYwQakHQcgYk1OYDj7vtUsB955gU6rUXOKnhuQMjzM8GlAoQ1nk/SwGpvspFThnK5athfOpP5swzx+KbvSma8Z6/usKJ2etZF2Ra+iL3E4Sa99fdPz+3v47PKvGe0Pw+tm7d2ulbsmRJp50DkN/4xjd2+rZt29Zp55Lm3rz1gcDZFPXP5z+3ttJnuMZcSSNue3UKnp9RhFALgpbTh/dzt5ktr9m3A1jotksB9/mY7ZJmA3OApxueOzBCqPVBnaZWWn/SawF+8jprar5ct+ekk04CxmspvkR2nlz3ieN+8ttrhfnL7Pf7X+38DHXaWdP/DHWaWr6vv05d+lhO2vfal3ckZA1t/vwxq8XHtOX4Nf+u/bvIbR+z5z+jNmttA4pT2wgslbSYSiCtYqygRWYdVSD+v1IF5t9mZiZpHfAPkv6KylGwFLhzEIMqEUItCFrOIISame2T9BFgPTALWGNmmyVdA4yY2TqqAPyvSRqlmvxflc7dLOlGYAuwD7gyeT6R9E3gbVSm73bgajObUCB/CLUgaDGDTIEys1uAW7r6PuXaLwHvrzn3WuDaQv+lAxmcI4RaA/pxFJTMGG8qZvPTm1aePCHuzSxvPmYHgTetSmlY3ePNlEzNXpPndSZp7q9Lwyqt4OTH6tftzM4Svxaod6BkZ0k2z7vb2az177pkfvZylLSRYXnOTAi1IGg5IdSCIGgVw5b7GUKthpLp1o/56X8dvUcux2Y9//zzxfu+4Q1vAMaqVsB4885X4SiNq2QKTvSXuq6KR8ks9+R34WPHPD6NKV8jF4OE8WXMs/no49xKizT72Df/LnK7zvwsPVcbCke2raxQE0KoBUHLCaEWNKK0qIdv+195H9meJ/r9+pmerH3UrQuatZ467agUh9YP/dQYK/1nqYtZK1Eqve01Ua+15ef2cW5+QZucsF6nVfYqXd5mhu2ZQ6gFQcsJoRYEQWsYxiKRjYVaqqc0Auwws3endIm1VFn4m4DfSiVJWkcvM6pUl8v/OpbMx7rzs9npE+L9+aUkdU9pcnuyfqnzGOrqtZXMX7/fO1hyLJ6vDedjzvJ7884Fb4rm83vVjos4tfbTz6TLR6kWXMnECu1BMAMYtoVXGgk1SQuA3wCuT9uiWqH9pnTIDcB7J2OAQRBMjGETak3Nz78GPgbkIKGTaLhCe5uo8yb2Mk+9yVOKIyvFvPUyo5qsAFUyhQdJqZx3P95P/4yl2LDSu6qLFezl6R22eSVPmwRWE5qs+/lu4Ekz23QwN5B0haQRSSOlwNEgCCaPplpamwRfE03trcB7JP06cCRwHPB5Gq7QnipnrgaYN2/ejH5zdb/2vb4QddkHpeuWIt/rFgg5UF+TcU2UXhkFve5fmrTv9a5KsWfd7RIzNSNgEAybltrzkzazq8xsgZktoqqPdJuZ/SaxQnsQzAiGTVObyM/Xx4E/SQXhTmLIVmgPgpnCsAm1voJvzex24PbUHpoV2nt94KWJcm/u+JgznxpUOj/PO/rYtJJJVlfDrGRm9ZP61A+l2K+6RV5K+/0z5IR1/658TFp+b75GWql0t79naY3PknOizbRNYDUhMgqCoOWEUAuCoFWEUAsaUed5K63glCtIwFgJ67o0qFwPzIe/+LLX2QyrC48pxa9NtPJGHf14Ykv4lLH8Dny9NF9TLj/33r17O30+pSpXQqmLk8tm5zCYnN0MyvspaQVV5MMs4Hozu65r/xHAV4HzqZbG+6CZPZr2XUWVdfQq8Edmtr7JNQ+G4fVzB8EQMKg4tZT7/UXgXcAy4FJJy7oOuxx4JqVOfo4qlZJ03CrgLGAF8LeSZjW8Zt+EplZD6Re/LkaqVGHVawQ+8TrXS/Maia8b9uSTT+7XN2fOnE47J7x7Ta0uSbxX8ntT6rSfUkaBp9f9fcJ6vobXzk4++eROO2u+vmKwr3Kbr1X3LkqZGr3i/9oS2zYg8/MCYDQ5CJG0FlhJtexdZiXw56l9E/CFlFK5ElhrZi8Dj6SIiexk7HXNvmnHpxYEQS0DCumYDzzutkupkZ1jUlD+Xqpwr7pzm1yzb0JTC4KW04emNlfSiNtenTKCZhQh1BpQMrN6mZ8eH2+V16qcO3dup8+bmqX1LxcuXNhpZ6eBN73q4sRKX+aSSVWXPF+6fqm/V5qUv75/P7m0OYyZh36907wGqr9XXje1u51j1vxYvLOmVG9tGJwG1l+RyN1mtrxm3w5godsupUbmY7ZLmg3MoXIYHOjcXtfsmzA/g6DlDMj83AgslbRY0uFUE//ruo5ZR5UyCVUK5W1WXXgdsErSEam47FLgzobX7JvQ1IKg5QzCUWBm+yR9BFhPFX6xxsw2S7oGGDGzdVSpkl9LjoA9VEKKdNyNVA6AfcCVZvYqQOmaEx1rCLU+qDM/fbxVbnszy8ekZfPztNNO6/T9+Mc/7rSzSTU6OtrpO/PMMzvt7P30Jq2Pc/MmVx5DnfeyZJaUyoH3Ml/qzLicuuRNPp/a5D24J5xwAgBLly7t9HkTPcekebN89+7dnXaOU+tlfnqPaFu8m70YVPCtmd0C3NLV9ynXfgl4f8251wLXNrnmRAmhFgQtJzIKgv3ox1GQtSavPZUcBYsWLer0+Qnx7Ch48MEHO31nnXVWp33GGWcA47MMvNPAU1ot3msnJQ2sH+2ltJZmL0dD3cr0+X3k54Px7+2RRx4B4LHHHuv0eQdL/gz8OT6To+QoGAYioT0IgtYxbEUiQ6gFQcsJTS2opc789DW+8oR1/hfGm0E5PeqNb3xjp2/JkiWddk6T8mbW3Xff3WnnNKLjjjuu0+dN3VKi+2T9UpfKeXtTPL8vnw7lx+cdAeeddx4w/r34OLZt27YB8Oijj3b6nnvuuU47m5U+Jc3XrusnTq1tJmoItSAIWkPMqQVB0DpCqBWQ9CjwPFUtpH1mtlzSicC3gEXAo8AHzOyZyRnm9KPO+5nj1LzJ5U2i3D711FM7fd67mb2fPk7t3nvv7bSz+Xnuued2+rzHrxQT5mPDSubhRL/0JfMbxt6B7/MVR/wznHPOOcD4mL5scgJs2VIVbti5c2enzz9Lfgd15mf2BA9DalQ3wybU+pk8+DUzO8/lhn0CuNXMlgK3pu0gCKYZr732WqO/tjCRGdGVwA2pfQPw3okPJwiCQTKoIpEziaZzagb8kyQD/i6VIznFzLIt8ARwymQMcDpR98H7X7nsiaxLB8peOF8k0ns/c0CpP/+pp57qtEdGRvYbi08t8l7RnFLl04V6pXT1WljZm5q57c1L74nNpp73cvqxnn/++Z12TpPavn17p897fXMqWamyh39Wb356U7aX+dk2j6enTQKrCU2F2q+Y2Q5JbwA2SHrQ7zQzSwJvPyRdAVwB4/8jB0EwNYRQK2BmO9K/T0q6maoU7y5Jp5rZTkmnAk/WnLsaWA0wb9681rzdXjFrXnvxC4RkrcmnOfmy1WeffTYwXqO65557Ou2cEuX7/PV97bVcOtxrL3XrkZaeK7dL63P6fq8JHnPMMZ121r78mE4//fT99sOYA8A7Re67775O2yevl+5V0tT8uNqsifVi2IRaz09a0jGSjs1t4GLgfsbXTroM+O5kDTIIgoMjF4kcJkdBE03tFODmNBcxG/gHM/tHSRuBGyVdDjwGfGDyhhkEwcEybJpaT6GWVno5t9D/NHDRZAxqplEyRb356E2f0qpG3mTKddbq6rHl+DU/Ye5Th7wpms3abIbC+Ji2fN26EuC539+/1PbX9PfKDgJvXvvYsccfH1tz46GHHgJg8+axGoG+dlo2e7OZCeNN+DyGknOgjmExSUOoBUHQKkKoBUHQKkKoBX1T5wnNeFM0x6zVmXzZPPPlvr2ZlMNivGnmq1X49Kwc3+Zj3np5DEvmpTcZfcWRfL6/Zskk9CZxTgPrfoZsQu/atavT503wfA9/L9/O4/ImZ5sXKG5K2wJrmzBcn3AQDCFT4f2UdKKkDZK2pX9PqDnusnTMNkmXuf7zJd0naVTS36SV3ZH0fkmbJb0mqW75vnGEpjZg8q9iXRS+15q6z4HyhLefaM/9PjHcr3/pHQilhU9Kpce9dlOK7fJ9pRplfvxeK8sxdb7Pj9VnSuRjS0nqMKYVekdBr3ppnmHT0DxTpKnlXPDrJH0ibX/cH5CKYFwNLKfKUtokaV0qhPEl4PeBO6gWYlkBfJ8qfOw/AX/XdCDD+0kHwZAwRbmfTXLBLwE2mNmeJMg2ACtS8P5xZvajtE7oV/P5ZvaAmW3tZyChqQVBi5nCObUmueDzgcfd9vbUNz+1u/sPihBqk4T/InmTKuPN0FISeWklJBhzFJT6YLz5mZ0GPsnck02yusn1/Ax15nO+rn9Wf2y+vx+TX03KOzUy3uQs1aGrey/DbF72og+hNlfSiNtendIcAZD0A2De/qfxya771eaCTwUh1IKg5fThBNjt6iXuh5m9o26fpCa54DuAt7ntBcDtqX9BV/8ODpL4eQuCFjOF9dSa5IKvBy6WdELyjl4MrE9m63OSLkxez9+uOb8RoalNASVTtC62Lbe9yVgqS+29oz42zB+bK4X4OLk6U7REaaz+Wnms3uQs1WvzFUv8sSWvqjcv/bPktjeVe5mcYZJWTNGc2nUUcsFTGMYfmNmHzWyPpE8DG9M515hZXpH6D4GvAEdReT2/n87/j8D/Ak4G/q+ku83skgMNJIRaELScqRBqdbngZjYCfNhtrwHW1Bx3dqH/ZuDmfsYSQm2KKcWxebJWU4ongzFNxWs5Hq/JlLQbr+Hle5QWY6nbX6ryWxeTl8/zGlNJ6/TP45/Lt5sunBLa2f4MW0ZBCLUgaDkh1IIgaA25SOQwEULtEFHnKCitxenNv2yGeZPUm1ylOLNSjTTfrkuj6jXWEr0m8n1fqc5cXcxcmJ0HT2hqQRC0ihBqQRC0ihBqBSQdD1xP5XI14PeArcC3gEXAo8AHUpJq0CelL12dxzG360yzkhlWZzKWzNMSva7fy/RrMtbcH2bmYIl6avV8HvhHMzuTar2CBxgrNbIUuDVtB0EwzYgV2ruQNAf4VeB3AMzsFeAVSSsZy+O6gSqH6+P7XyHoh15frqzB1Wk/XsMraT2l65fW+vT00q56UXdc0/O77xv0x7B5P5t8UxYDTwF/L+kuSden9T+blBoJguAQM2yaWhOhNhv4ReBLZvYW4AW6TM1U2K34ViRdIWlE0kiuzx8EwdQwhQnt04YmjoLtwHYzuyNt30Ql1JqUGiHVY1oNMG/evPa8uSmk1xduKr+Q/cSOHQxhZg6eNgmsJvT8BpnZE8Djks5IXRcBW2hWaiQIgkNMaGpl/ivwDUmHAw8Dv0slEPcrNRIEwfRi2BwFjYSamd1NtQJMN/uVGgmmhqn4Zc2m4GT9pwhTc/JpmxbWhMgoCIKWE0ItCIJWEUItCLqoMxN7maVhXk4PQqgFQdAqQqgFQUNCE5v+DGORyPhWBkHLmYo4NUknStogaVv694Sa4y5Lx2yTdJnrP1/SfZJGJf1NWioPSX8p6UFJ90q6OVUMOiAh1IKg5UxR8G3Pqj2STgSuBn4JuAC42gm82g+oAAAIIklEQVS/LwG/DyxNfytS/wbgbDM7B/gxcFWvgYRQC4KWM0VCbSVVtR7Sv+8tHHMJsMHM9qTaixuAFSnN8jgz+1HKI/9qPt/M/snM8hJlP2L8Su5FYk4tCFpMnwJrrqQRt7065W43oUnVnvnA4257e+qbn9rd/d38HlVh2gMSQi0IWk4fQm23mZUyhwCQ9ANgXmHXJ7vuZ5IG6nKV9ElgH/CNXseGUAuCljMo76eZvaNun6QmVXt2MFZYFipT8vbUv6Crf4e79u8A7wYusgYSOubUgqDlTNGcWpOqPeuBiyWdkBwEFwPrk9n6nKQLk9fzt/P5klYAHwPeY2aNCjKGUAuCFjOFRSKvA94paRvwjrSNpOWSrk9j2QN8GtiY/q5JfQB/SLW40yjwEPD91P8F4Fhgg6S7Jf3vXgMJ8zMIWs5UZBSY2dMUqvaY2QjwYbe9BlhTc9zZhf4l/Y4lhFoQtJxIkwqCoFUMW5pUCLUgaDFRJDIIgtYxbEKtp/dT0hnJ65D/npP0x00TWIMgOLQM28IrTVaT2mpm55nZecD5wM+Am2mQwBoEwaEnhNqBuQh4yMweo1kCaxAEh5hhE2r9zqmtAr6Z2k0SWIMgOIREkcgDkNb8fA/wf7r3pXysoqiXdIWkEUkjP/tZoyyHIAgGyLBpav2Yn+8C/s3MdqXtXSlxlQMksGJmq81suZktP/rooyc22iAI+iaEWj2XMmZ6QrME1iAIDjEh1ApIOgZ4J/Ad111MYA2CYPowhQnt04ZGjgIzewE4qauvmMAaBMH0ok0CqwmRURAELWfYvJ8h1IKg5YSmFgRBa2jbfFkTQqgFQcsJoRYEQasIoRYEQasYNkdBLLwSBC1mquLUmpYik3RZOmabpMtc//mS7pM0Kulv0qpSSPq0pHtT2bN/knRar7GEUAuCljNFwbc9S5FJOhG4Gvgl4ALgaif8vgT8PrA0/a1I/X9pZuek0mffAz7VayAh1IKg5UyRUGtSiuwSYIOZ7TGzZ4ANwIqUO36cmf0oFcf4aj7fzJ5z5x9DTeEMT8ypBUHL6UNgzZU04rZXm9nqhuc2KUU2H3jcbW9PffNTu7sfAEnXUi1wvBf4tV4DCaEWBC2nD6G228yW1+2U9ANgXmHXJ7vuZ5IG5nI1s08Cn5R0FfARKhO2lhBqQdBiBlkk0szeUbdP0i5Jp5rZzgOUItsBvM1tLwBuT/0Luvp3FM7/BnALPYRazKkFQcuZojm1JqXI1gMXSzohOQguBtYns/U5SRcmr+dv5/MlLXXnrwQe7DWQ0NSCoOVMUfDtdcCNki4HHgM+ACBpOfAHZvZhM9sj6dPAxnTONWa2J7X/EPgKcBTw/fQHcJ2kM4DX0nX/oNdAQqgFQcuZCqFWV4rMzEaAD7vtNcCamuPOLvT/537HEkItCFpMJLQHQdA6QqgFQdAqIvezgKT/JmmzpPslfVPSkZIWS7oj5Wp9Ky2hFwTBNGPY1ijoKdQkzQf+CFhuZmcDs6gWNf4M8DkzWwI8A1w+mQMNgqB/hnHhlaZxarOBoyTNBo4GdgJvB25K++tyvYIgOMSEUOvCzHYA/xP4CZUw2wtsAp41s33psHG5Wp5YoT0IDi0h1LpIkb8rgcXAaVSZ8isOeJLDYoX2IDikvPbaa43+2kIT7+c7gEfM7CkASd8B3gocL2l20tbqcrXGsWvXrt2f/exnXwB2T2DM05W5xHPNFNr4TFA91zFdfetTfxNa8U6aCLWfABdKOhp4kSpqeAT4IfA+YC31uV7jMLOTJY0cqBLATCWea+bQxmeCznMt8n1m1tiqagtN5tTuoHII/BtwXzpnNfBx4E8kjVKt3v7lSRxnEARBIxoF35rZ1exf7uNhqpK8QRAE04ZDUXqoaSXNmUY818yhjc8E7X2uvlCbXLlBEARRJDIIglYxpUJN0gpJW1O+6H5LaM0EJC2U9ENJW1I+7EdTf6N1D6c7kmZJukvS99L2jM/xlXS8pJskPSjpAUm/PNM/r8jHrmfKhJqkWcAXgXcBy4BLJS2bqvsPkH3An5rZMuBC4Mr0HD3XPZwhfBR4wG23Icf388A/mtmZwLlUzzdjP6/Ixz4wU6mpXQCMmtnDZvYKVXzbyim8/0Aws51m9m+p/TzVf5D5NFv3cFojaQHwG8D1aVvM8BxfSXOAXyWFHJnZK2b2LDP/84p87BqmUqjVrfk3Y5G0CHgLcAfN1j2c7vw18DGqevBQxR82yvGdxiwGngL+PpnV10s6hhn8eU00H7vthKPgIJH0euDbwB93rSJNWmV6RrmVJb0beNLMNh3qsQyY2cAvAl8ys7cAL9Blas60z2ui+dhtZyqF2g5godtulC86HZF0GJVA+4aZfSd170rrHXKAdQ+nM28F3iPpUaqpgbdTzUUdn0wcmJmf2XZge8qMgco8+0Vm9ufVycc2s58D4/Kx0zEz8bMaCFMp1DYCS5OH5nCqic11U3j/gZDmmb4MPGBmf+V2NVn3cNpiZleZ2YKUO7gKuM3MfpOxHF+Ymc/1BPB4WmYNqtzlLczsz6uTj52+j/mZZvRnNSimNPhW0q9TzdvMAtaY2bVTdvMBIelXgP9HlQeb557+jGpe7UbgP5DWPXRrGs4oJL0N+O9m9m5Jb6LS3E4E7gI+ZGYvH8rx9Yuk86icH4dTpff9LtUP+oz9vCT9BfBBKm/8XVTL0M1nhn9WgyAyCoIgaBXhKAiCoFWEUAuCoFWEUAuCoFWEUAuCoFWEUAuCoFWEUAuCoFWEUAuCoFWEUAuCoFX8f3zh7JWG9NgdAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAATUAAAD8CAYAAAAWjzPMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztnXuQXdV1p78fEu+HhCSMQNJYspFxAQE5qAgpp1KOsUFOXJZnxg+YOCEJDpUKnjiTTNkmrjIJHqpwZRzHGTueqLBi/IhlBpuyyoOjyGAqk6oYaAUMSCCQeRgpQiAkxMM8LFjzx9n79uruffqe23271X3u+qq6ep99Xvuce3v1Wns9tsyMIAiCtnDYoR5AEARBPwmhFgRBqwihFgRBqwihFgRBqwihFgRBqwihFgRBqwihFgRBq5iUUJO0RtJ2STskfaJfgwqCIJgommjwraQ5wIPAO4GdwJ3AJWa2rX/DC4Ig6I25kzj3PGCHmT0MIGkDsBaoFWpHH320zZs3bxK3DIJgPA4cOMCLL76ovL1mzRrbu3dvo3O3bNmyyczWTNngponJCLUlwONueyfwS+OdMG/ePD70oQ9N4pZBEIzH17/+9RHbe/fuZWhoqNG5khZNxZimm8kItUZIuhy4HOD444+f6tsFQTCKQcvvnoxQ2wUsc9tLU98IzGwdsA5g8eLFg/V2g2AG8Nprrx3qIUwrk/F+3gmslLRC0hHAxcDG/gwrCIJ+YGaNf7rRLdpB0pGSvpX23y5pudt3ZerfLumiUefNkXSXpO/14ZEnLtTM7CDwEWATcD9wg5lt7ceggiDoH/0Qaina4YvAu4AzgEsknTHqsMuA/WZ2GvA54DPp3DOolJ4zgTXA36brZT5KJUP6wqTi1MzsZjN7k5m90cyu6deggiDoH33S1DrRDmb2CpCjHTxrgetT+0bgAklK/RvM7GUzewTYka6HpKXAbwDX9eVhiYyCIGg9PQi1RZKG3M/l7jKlaIclo27VOSZZcgeAhV3O/WvgY0DfJv6m3PsZBMGhpQfv514zWz2VY/FIejfwpJltkfS2fl03hFoQtBgz65f3s0m0Qz5mp6S5wDzg6XHOfQ/wHkm/DhwFnCDp62Y2qWDWMD+DoOX0aU6tSbTDRuDS1H4fcKtVF94IXJy8oyuAlcAdZnalmS01s+XperdOVqBBaGpB0Hr6EXxrZgcl5WiHOcB6M9sq6WpgyMw2Al8GviZpB7CPSlCRjruBKoXyIHCFmb066UHVEEItCFpOvzIKzOxm4OZRfZ9y7ZeA99ecew1QGyFhZrcBt/VjnCHUgqDFNA2sbRMh1IKg5QxamlQItSBoOaGpBUHQGsL8DIKgdYRQC4KgVYRQC4KgVYRQC4KgNfQxTWrWEEItCFpOaGpBELSKEGpBELSKEGpBELSKQRNqXUsPSVov6UlJ97m+BZI2S3oo/T5xaocZBMFEyI6CJj9toUk9ta9QLZbg+QRwi5mtBG5J20EQzED6tZrUbKGrUDOzf6aqjeTxCyxcD7y3z+MKgqBPDJpQm+ic2slmtju1nwBO7tN4giDoM20SWE2YtKPAzExS7VtLK9JcDnD88cdP9nZBEPRA27SwJkx0jYI9kk4BSL+frDvQzNaZ2WozW33MMcdM8HZBEEyUQTM/JyrU/AILlwLf7c9wgiDoN+H9HIWkbwL/Cpwuaaeky4BrgXdKegh4R9oOgmAG0i9NTdIaSdsl7ZA0JuIhrRb1rbT/dknL3b4rU/92SRelvqMk3SHpx5K2SvqLfjxv1zk1M7ukZtcF/RhAEARTR79MS0lzgC8C76RaYf1OSRvNbJs77DJgv5mdJuli4DPAByWdQbWy1JnAqcAPJL0JeBl4u5k9L+lw4F8kfd/MfjSZsca6n0HQcvqkqZ0H7DCzh83sFWADVWiXx4d63QhcIEmpf4OZvWxmjwA7gPOs4vl0/OHpZ9ISOIRaELScPgm1JcDjbntn6iseY2YHgQPAwvHOlTRH0t1UzsbNZnb7BB+zQwi1IGg5PQi1RZKG3M/l0zC2V81sFbAUOE/SWZO9ZiS0B0GL6bFI5F4zW12zbxewzG0vTX2lY3ZKmgvMA55ucq6ZPSPph1QpmfcxCUJTC4KW0yfz805gpaQVko6gmvjfOOoYH+r1PuBWqy68Ebg4eUdXACuBOySdJGk+gKSjqZwQD0z2eUNTC4KW0w/vp5kdlPQRYBMwB1hvZlslXQ0MmdlG4MvA1yTtoMoXvzidu1XSDcA24CBwhZm9mgL3r0+e1cOAG8zse5Mdawi1IGg5/coWMLObgZtH9X3KtV8C3l9z7jXANaP67gHe0pfBOUKoBUHLaVMKVBNCqAVBi+nRUdAKQqgNMFVc5OQZNE1gtjFon08ItSBoOSHUgiBoFSHUgtZQMi+7mZzd9pf+QPw5pf2D9kc1k2hbrbQmhFALgpYTQi2YFdRpVL4/t0t9410jc9hhVcJJN++Z/6Mptev2110j6C/h/QyCoFUM2j+MEGpB0GJiTi2Y8ZRMymwm1rV72d8L2azx5k2p3W2/Z9D+AKeDQXunTdYoWCbph5K2pTriH039CyRtlvRQ+n3i1A83CIJeidWkxnIQ+FMzOwM4H7gi1Rz/BHCLma0EbknbQRDMMAZNqDVZeGU3sDu1n5N0P1Up3rXA29Jh1wO3AR+fklEOIHVeypLJOGfOnE577ty5Y/qPOOKI4rG57a9Vum+d9zKbj6+++mqnz7dfeeWVMX0HDx4sHlsyVZt4TYPxidzPLqQlr94C3A6cnAQewBPAyTXnxArtQXAIGbR/CI1niCUdB3wb+GMze9bvS9Uti2/OYoX2IDikhPlZIK3J923gG2b2ndS9R9IpZrY7VbB8cqoGOUhkk6/OY5nNy8MPP7zT583LI488stPOx0zU/CwxEfPz5z//eafv5ZdfHrPfH+PN05KntE1/fNPFoL2zJt5PUZXpvd/M/srt8vXILwW+2//hBUEwWUJTG8tbgd8C7k3r8wH8GXAtcIOky4DHgA9MzRDbT8kR4DWqkiZ21FFHdfqOPvroMfthWFPz1/LkL3Ld5HyJktZYp91lra1OU3vxxRc77ZdeemnMfq/JZXoZaxCOgiJm9i9AnU1yQX+HEwRBv+mX8Je0Bvg81cIr15nZtaP2Hwl8FTiXamm8D5rZo2nflcBlwKvAH5nZJknL0vEnU83JrzOzz092nLFEXhC0nH6Yn2nFpy8C7wLOAC5J8aqey4D9ZnYa8DngM+ncM6hWljqTal3Pv03Xq4uBnRSRJnWI8CZbKc7Mm5zevDz22GMB8J5kf6w3D/OkezbtoLvJ503FbLb4a3oHRTdTOB9bZx778372s58B8MILLxSfJZuidXFuYYrW06d3cx6ww8weBpC0gSpWdZs7Zi3w56l9I/CFNCe/FthgZi8Dj6Ql9M4zs3+lHAPrr9kzoakFQYtpqqU1EHxLgMfd9s7UVzzGzA4CB4CFTc4dFQM7KUJTC4KW04OmtkjSkNteZ2brpmBIIxgvBnYihFCbZrLZWZfalE0yb14ed9xxnXbu9+d4k+zAgQOd9tNPPz3i9+j2M888A8Dzzz/f6fMex2ze1Xli87jmz5/f6Vu4cOGYtu/zWSXe/MzPU0rzgmHz1JvSnjzWMEPH0oP3c6+Zra7ZtwtY5raXpr7SMTslzQXmUTkMas+tiYGdFGF+BkHL6ZP5eSewUtIKSUdQTfxvHHWMj119H3BryjbaCFws6UhJK4CVwB3jxMBOitDUpoFSHFpJO4NhR4DXaErpZV4j2717d6f92GOPddqPPvooAP/+7//e6du7d2+n/dxzzwEjnQO9kB0BfqyLFi3qtE899VQAli9f3ul7/etf32mfcsopnXa+htdKu2U6eK1tIjF3g0C/AmvN7KCkjwCbqEI61pvZVklXA0NmtpFKQH0tOQL2UQk+0nE3UDkADgJXmNmrkn6FQgysmd08mbGGUAuCltMv4Z6Ezc2j+j7l2i8B76859xrgmlF948XATpgQakHQcgZNYw2hNkXUldsu1TgrOQW8SerNw6eeegqAHTt2dPq2bt3aafv+J554oi/jrjPp8rj27dvX6fPtBx98EIDFixd3+k477bRO+8wzzxzTf9JJJ3X6/DvI9y0l1I9ul/oG7Q/bM2jPHkItCFpM5H4GQdA6QlMLJkW31Z6y2elTh3w77/eePe+9vO+++wC46667On0PP/xwp136r+zjxLwp+LrXvQ4YGWfmTeFsKvt0pBwvBsNxbk8+OVxKz5u8OSbO99Udm83Ws846q9OXvadQfm+lcuA+5av0xzxof+AweM8cQi0IWk4ItaAvlOqOwbDGUap7BsMamtfOfvzjH3faW7ZsAWDnzp3F669YsaLTfvOb3wyMnJxfsmQ45W7BggXASO3Mj6XkKPBOi6y1eefArl3DQebZafHAAw90+nwcnXdq5OcuZTTAsNZWSqj346pLcvf9g0YItSAIWkM4CoIgaB2hqQU9U1ors5v56fu8yZXj0LZtGy4p5c3PvN9P+J999tmd9qpVqzrtbIp6R4C/b6nctm/nP4Y6p0eOqZs3b16nb9my4bzlbP76eLS77767077nnnvGPJd/1pIp7OPYSu/Vj7/kSPAMyh/7oDxnpsnCK0dJukPSjyVtlfQXqX+FpNsl7ZD0rZTkGgTBDGPQFl5pUqXjZeDtZnYOsApYI+l8qlK9n0ule/dTlfINgmAG0ccikbOGJguvGJALbh2efgx4O/BfUv/1VGV8v9T/Ic4uSqtBldreNMrVMmDYO+g9g77eWTbvVq8eLnvlTU5vnmV8ZQ5/rey9rPM4lszPUm21brXfTj/99E5f9riObg8NVbUJc+wbjHwHuYqHN0l9dZA8rrr3XlqjdFBok8BqQtPFjOcAW4DTqBZf+AnwTCrZC+XSvvncy4HLYeSXMAiC6SG8nwXM7FVglaT5wE3Am5veIJUDXgewePHi1vzLKDkHfLvOUZDxC6DkSXIYjvPyGpWvUXbOOeeM+A1wwgkndNr79+/vtPM1/L18ldymz+IpZReUsgxgOPrfa29+rP4ZsjbhHQX+HeT34p0S/r1mDc73eadBKdPD01Ztpm2mZRN68n6a2TOSfgj8MjBf0tykrZVK+wZBMAMYNKHWxPt5UtLQkHQ08E7gfuCHVCV7oSrh+92pGmQQBBMnHAVjOQW4Ps2rHQbcYGbfk7QN2CDpfwB3UZXyHUi6mZ+ebBJ508ovhvLss9ViOr6WmC+HvXLlSmBkYrcv7e3Xzcxf1FI81+hxlyg5CsY7DsopVd4k9Sahf4b8XH78uRw5DL8X/678HK03cTOlcuD+Wdr0hzweg/KcmSbez3uo1uMb3f8w1QKnQRDMUAYxTSpWkwqCltMv81PSGknbU8D9Jwr7j0yB+DtSYP5yt+/K1L9d0kWuf72kJyXd16fHjTSpfpNNnjozJ8eEefPTtzM+3synHmXvoa8b5r2PnjwG7/H07dJamd3+q5dMOh8P5p+75PX1Y/XXys/ln9WbotmDW/feslldN5Z8r4hTmxhp+umLVHPqO4E7JW00s23usMuA/WZ2mqSLqQL0PyjpDKqVpc4ETgV+IOlNKariK8AXgK9OepCJ0NSCoOX0SVM7D9hhZg+b2SvABmDtqGPWUgXiA9wIXJDW9lwLbDCzl83sEWBHuh5m9s9Uy+n1jRBqQdBy+iTUlgCPu+1SwH3nmBTqdQBY2PDcvhHmZx8oeQf9l8SbfNn8rAuIzR5Bb376Khv5Xv78usWIs6nlTcrSWL0ZWDIv677wub9bYUZvEnp8yfKcUuWf1b+DHKDs35V/B9kc997d0ri7eXLbRo/hGoskDbntdSl4flYRQi0IWk4P3s+9Zra6Zt8uYJnbLgXc52N2SpoLzAOebnhu3wih1mfyF6huzcmsqXntymsyOfbKayw+Zi2f77UcrymVNK3ShP3ocZX6umk6TVOP6hwR/hnyGP2z+neQj/XP599hfi/+XfbiAGkzfYpTuxNYKWkFlUC6mOGCFpmNVIH4/0oVmH+rmZmkjcA/SPorKkfBSuCOfgyqRAi1IGg5/RBqZnZQ0keATcAcYL2ZbZV0NTBkZhupAvC/JmkH1eT/xencrZJuALYBB4ErkucTSd8E3kZl+u4ErjKzSQXyh1ALghbTzxQoM7sZuHlU36dc+yXg/TXnXgNcU+i/pC+Dc4RQ6wOlL02d+Vma1PeT29l88ilEnmxm1cVblUyu6TTDSuZvXcqYf4b8XP5Y/w6OPfbYMed48nv19da6OXAGhUF75hBqQdByQqgFQdAqBs1JEkKtz3SL3cp4M8t7/LL55M0of61sZtV9UUvey7o4tKmm5An2lLzC/rl9O8exefO99Cz+XQ1iRY7RtK2sUBNCqAVBywmhFvRMaSLeawklTcXHjvljs3bi+0prWfp7douS9/efTkdBN0rZBx7/XFmbrVtYpXQdP5ZSzNygMGjPHEItCFpOCLUgCFrDIBaJbCzUUj2lIWCXmb07pUtsoMrC3wL8VipJMhDUpRM1TZg+8sgjxz2/LiG+W2qTJ3+ZvRlWMtkmmuTdbSyl+9dRiqnz7dJqUd3i0Eom7aBpLTB4z9yLK+yjVAuuZGKF9iCYBfSr8u1soZFQk7QU+A3gurQtqhXab0yHXA+8dyoGGATB5Bg0odbU/Pxr4GNAXr5nIQ1XaG8rE52nyKZYXdnpfN06k7Nb6lHpXodqTqUXT2g3r242P0seTX/+RMt1T2f83nTTJoHVhCbrfr4beNLMtkzkBpIulzQkaaiuln4QBFNDUy2tTYKviab2VuA9kn4dOAo4Afg8DVdoT5Uz1wEsXrx4Vr+5JhpPKR6qtACIj5b3k99ZQ6tbLKWuimxTSgnnk6UuYb4XDTEf6zUt/45y0n/pXfkx1GVydKv95sfaNq1t0LyfXT89M7vSzJaa2XKq+ki3mtlvEiu0B8GsYNA0tcn8S/o48CepINxCBniF9iCYyQyaUOsp+NbMbgNuS+1YoT1RMr+6mTPetCqV6/bzjz5NKt/L12Drhbo6b5OhH2Zsyez27yjH9fnn9uXA83urM3/zZzAV5vdMpm0CqwmRURAELSeEWhAErSKEWjCGbt6jkvnpvXAlj573YnrzM3v3nn/++eL5o68z3lhKZm+3dT1LNcjq6pKVzi+9q27rjsKw2Vnn/cz11Lz305vl+R51a6CW7t/NE9oWL2i/vJ+S1lBFPswBrjOza0ftPxL4KnAu1dJ4HzSzR9O+K6myjl4F/sjMNjW55kRox6cWBEGRfsWppdzvLwLvAs4ALpF0xqjDLgP2p9TJz1GlUpKOuxg4E1gD/K2kOQ2v2TOhqfWZ/OWo+++YNRL/JfKT33lCvG7hlbwSeV2cWykhvC7Ju9tq7SVNrfRcdZpeN03QOwXyc/lnKS28Upfwnq9V997bonVNhD6Zn+cBO5KDEEkbgLVUy95l1gJ/nto3Al9IKZVrgQ1m9jLwSIqYyE7GbtfsmcH9pINgQOhTSMcS4HG3XUqN7ByTgvIPUIV71Z3b5Jo9E5paELScHjS1RZKG3Pa6lBE0qwihNkXUmXw5nsqbXt40ymbW/PnzO3179+7ttF988UVg2FyD+gnvbqlBpbHWtUv3ahrnVfcu/DPkCf78/DDyHRx33HEAvPDCC50+/w5L73WyKWVtwHorErnXzFbX7NsFLHPbpdTIfMxOSXOBeVQOg/HO7XbNngnzMwhaTp/MzzuBlZJWSDqCauJ/46hjNlKlTEKVQnmrVRfeCFws6chUXHYlcEfDa/ZMaGpB0HL64Sgws4OSPgJsogq/WG9mWyVdDQyZ2UaqVMmvJUfAPiohRTruBioHwEHgCjN7FaB0zcmONYRaAyZSl6yb+ZnNSBgZW5U9fosWLer0Pf300532gQMHxpxfV5stm2J1JbDzeaUqIr6/WxxaXWUMbwqW9vtnyPc68cQTO33+HeT38swzzxTPz++1bjWpXmibp7RfwbdmdjNw86i+T7n2S8D7a869BrimyTUnSwi1IGg5kVEQ9Ey3yXOv3WTtxWcM+MnvfP5JJ53U6fPJ7Vlr884Dr6nlyHsoa5i+nbWbyVKn6ZXu77Ur7yjIWtmyZcPzxv4d5Hfk35V/h/ladVpWU6dJ24iE9iAIWsegFYkMoRYELSc0tWBSdFsYJZufzz33XKfPOwIWL14MjDS93vjGN3baJTPMT557k/L446t1cnzCfLc6bA1yAMfd7yfqc70z/6zelJ43b16nvXLlSmDks/qxPvHEE8DId+Wvm9+rP2fQaqfVEUItCILWEHNqQRC0jhBqBSQ9CjxHVQvpoJmtlrQA+BawHHgU+ICZ7Z+aYc4MvDlTVyOsVDbaeyfzF8x7AZ966qlOe9euKkvEm5/eI3juueeOuA7A3Xff3Wl7r2g2S71HNKcbwXDslzfZ/FhL3stuMXfeI1laEtHHnp1zzjljnsvHqT3++HCuc34v/l35++ZYvNL4fbvODG1bbJpn0IRaL5/kr5nZKpcb9gngFjNbCdyStoMgmGG89tprjX7awmT+Pa0Frk/t64H3Tn44QRD0k34ViZxNNJ1TM+CfJBnwd6kcyclmtjvtfwI4eSoGONvI5o03g3yaUjZzvBnnvZePPvooMOy5hJHVKrIp6s3IBQsWdNreFM3X8iZb9iL6sdYVnCylSfnUp5ze5ff7586m5PLlyzt9q1at6rTPPvvsMcfu3r270/eTn/xkzLP4d+UpLXZclz42aLRJYDWhqVD7FTPbJel1wGZJD/idZmZJ4I1B0uXA5TDyDzUIgukhhFoBM9uVfj8p6SaqUrx7JJ1iZrslnQI8WXPuOmAdwOLFi1vzdrs5Deo0tVyu2yex+wnvPXv2jLm+16Rye+nSpZ0+P/nu47y2b98OwEMPPdTpyxPuAPv27QNGTuj7+LfSAiR+LPmflNcUlywZLlyaY89OP/30Tt/rX//64rV27twJwD333NPp27p1uGBDfi9ew/XlvvN77UVTa7NzwDNoQq3rpyrpWEnH5zZwIXAfI2snXQp8d6oGGQTBxMhFIgfJUdBEUzsZuCn9p5sL/IOZ/aOkO4EbJF0GPAZ8YOqGGQTBRBk0Ta2rUEsrvZxT6H8auGAqBtUGvLnjTaI8oe1Tl0qm6E9/+tNOXyk2zE+4v+ENb+i0f+EXfqHTPu2004CRsWveUZAdCNkM9ff39/VmnDf5stnpY+pymhcMm8X+HB/H5s3i7OC49957O33eVM5j8dcqpX/V1Y4bZEKoBUHQKkKoBUHQKkKoBT1T8oSWymbDsJlUl3qU294L+dhjj3Xaud+bjN6kzB5HGDYFvcdxxYoVnXaOOfOVPboVjvQpVSWTzz9XHusjjzzS6fMm57Ztw2vWPvjgg8BIU9nHxOVVpnzKV8n8DI/nSNoWWNuEwfqEg2AAmQ7vp6QFkjZLeij9PrHmuEvTMQ9JutT1nyvpXkk7JP1NWtkdSe+XtFXSa5Lqlu8bQWhqfaaUBF7S2nJcFZQXNvF9Po4sa2V+wr0uCj9raD52bOHChZ12qd6an4jvllGQF4Gpqw2XJ/q9ppkzA0aP2z9PxmtlJU3Nv8PSIjKeQdPQPNOkqeVc8GslfSJtf9wfkIpgXAWspspS2iJpYyqE8SXg94HbqRZiWQN8nyp87D8Bf9d0IIP7SQfBgDBNuZ9NcsEvAjab2b4kyDYDa1Lw/glm9qO0TuhX8/lmdr+Zbe9lIKGpBUGLmcY5tSa54EuAx932ztS3JLVH90+IEGrTTDaDSrFrMJyoXrdCVZ58z6YfjDTdcjoRDE/K+zQqb37mRPls2sFIU7RkSucS3X4sPsncm5950t87NZ599tlO2ztIsvnox+KT9nO/f1elQgHBWHoQaoskDbntdSnNEQBJPwAWjz2NT466X20u+HQQQi0IWk4PToC9rl7iGMzsHXX7JDXJBd8FvM1tLwVuS/1LR/XvYoLEv7cgaDHTWE+tSS74JuBCSScm7+iFwKZktj4r6fzk9fztmvMbEZraFFFXxaO0v7QCUl2cW65s4T2iflHg/fuHK6pnU897GUueTu9FLJl0pcWY/X19apU3T/N+b2b65/ae1uzV9N7NblU4upmcYZJWTNOc2rUUcsFTGMYfmNmHzWyfpE8Dd6ZzrjazPDfxh8BXgKOpvJ7fT+f/R+B/AScB/1fS3WZ20XgDCaEWBC1nOoRaXS64mQ0BH3bb64H1NcedVei/Cbipl7GEUJsGetHasiZSp6nl/V6789qRb+fsAN/ntao8wV9aOKYOP/78x+L/aEoOjrok9FLb9/l6a/kdhHbWO4OWURBCLQhaTgi1IAhaQy4SOUiEUJtmSpPvpf11cWp5vzfN/ES/N/XyRL2f3PftPIHv/5N3+69eqldWV7o8t/34Sgnx/nkmunBKmJ31hKYWBEGrCKEWBEGrCKFWQNJ84Doql6sBvwdsB74FLAceBT6QklSDBtSZS6V6bKXYrFLsGpRN0W7mZ8mj6fv9WEumcDfzs9Q3+rzwbk4NUU+tns8D/2hmb6Zar+B+hkuNrARuSdtBEMwwYoX2UUiaB/wq8DsAZvYK8IqktQzncV1PlcP18bFXCHphIo6EOk0pX8NH9Pvrlmq3dfOUlbQ23+fbJe2rTusLR8DUMWjezybflBXAU8DfS7pL0nVp/c8mpUaCIDjEDJqm1kSozQV+EfiSmb0FeIFRpmYq7FZ8K5IulzQkacjnKwZBMPVMY0L7jKGJo2AnsNPMbk/bN1IJtSalRkj1mNYBLF68uD1vborpllpVZ7qVzFJvnnaLSev25S6Zid1i6rqdX0eYnP2hTQKrCV2/NWb2BPC4pNNT1wXANpqVGgmC4BATmlqZ/wp8Q9IRwMPA71IJxDGlRoIgmFkMmqOgkVAzs7upVoAZzZhSI0H/KZlhdV/UknnnzU9Pv/4792JS1hGm5tTQNi2sCZFREAQtJ4RaEAStIoRaMCvolmbVhH6Yjb0SZub0E0ItCIJWEUItmNX0ogn1yysW2tfMZRCLRMa3MQhaznTEqUlaIGmzpIfS7xNrjrs0HfOQpEtd/7mS7pW0Q9LfpKXykPSXkh6QdI+km1LFoHEJoRYELWeagm+7Vu2RtACTYt3dAAAINElEQVS4Cvgl4DzgKif8vgT8PrAy/axJ/ZuBs8zsbOBB4MpuAwmhNsAcdthhffkJZjbTJNTWUlXrIf1+b+GYi4DNZrYv1V7cDKxJaZYnmNmPUh75V/P5ZvZPZpaLAP6IkSu5F4k5tSBoMT0KrEWShtz2upS73YQmVXuWAI+77Z2pb0lqj+4fze9RFaYdlxBqQdByehBqe82slDkEgKQfAIsLuz456n4mqa8uV0mfBA4C3+h2bAi1IGg5/fJ+mtk76vZJalK1ZxfDhWWhMiVvS/1LR/Xvctf+HeDdwAXWQELHhEgQtJxpmlNrUrVnE3ChpBOTg+BCYFMyW5+VdH7yev52Pl/SGuBjwHvMrFFBxhBqQdBiprFI5LXAOyU9BLwjbSNptaTr0lj2AZ8G7kw/V6c+gD+kWtxpB/AT4Pup/wvA8cBmSXdL+t/dBhLmZxC0nOnIKDCzpylU7TGzIeDDbns9sL7muLMK/af1OpYQakHQciJNKgiCVjFoaVIh1IKgxUSRyCAIWsegCbWu3k9JpyevQ/55VtIfN01gDYLg0DJoC680WU1qu5mtMrNVwLnAz4CbaJDAGgTBoSeE2vhcAPzEzB6jWQJrEASHmEETar3OqV0MfDO1mySwBkFwCIkikeOQ1vx8D/B/Ru9L+VhFUS/pcklDkoZ+9rNGWQ5BEPSRQdPUejE/3wX8m5ntSdt7UuIq4ySwYmbrzGy1ma0+5phjJjfaIAh6JoRaPZcwbHpCswTWIAgOMSHUCkg6Fngn8B3XXUxgDYJg5jCNCe0zhkaOAjN7AVg4qq+YwBoEwcyiTQKrCZFREAQtZ9C8nyHUgqDlhKYWBEFraNt8WRNCqAVBywmhFgRBqwihFgRBqxg0R0EsvBIELWa64tSaliKTdGk65iFJl7r+cyXdK2mHpL9Jq0oh6dOS7kllz/5J0qndxhJCLQhazjQF33YtRSZpAXAV8EvAecBVTvh9Cfh9YGX6WZP6/9LMzk6lz74HfKrbQEKoBUHLmSah1qQU2UXAZjPbZ2b7gc3AmpQ7foKZ/SgVx/hqPt/MnnXnH0tN4QxPzKkFQcvpQWAtkjTktteZ2bqG5zYpRbYEeNxt70x9S1J7dD8Akq6hWuD4APBr3QYSQi0IWk4PQm2vma2u2ynpB8Diwq5PjrqfSeqby9XMPgl8UtKVwEeoTNhaQqgFQYvpZ5FIM3tH3T5JeySdYma7xylFtgt4m9teCtyW+peO6t9VOP8bwM10EWoxpxYELWea5tSalCLbBFwo6cTkILgQ2JTM1mclnZ+8nr+dz5e00p2/Fnig20BCUwuCljNNwbfXAjdIugx4DPgAgKTVwB+Y2YfNbJ+kTwN3pnOuNrN9qf2HwFeAo4Hvpx+AayWdDryWrvsH3QYSQi0IWs50CLW6UmRmNgR82G2vB9bXHHdWof8/9zqWEGpB0GIioT0IgtYRQi0IglYRuZ8FJP03SVsl3Sfpm5KOkrRC0u0pV+tbaQm9IAhmGIO2RkFXoSZpCfBHwGozOwuYQ7Wo8WeAz5nZacB+4LKpHGgQBL0ziAuvNI1TmwscLWkucAywG3g7cGPaX5frFQTBISaE2ijMbBfwP4GfUgmzA8AW4BkzO5gOG5Gr5YkV2oPg0BJCbRQp8nctsAI4lSpTfs24JzksVmgPgkPKa6+91uinLTTxfr4DeMTMngKQ9B3grcB8SXOTtlaXqzWCPXv27P3sZz/7ArB3EmOeqSwinmu20MZnguq5jh3Vtyn1N6EV76SJUPspcL6kY4AXqaKGh4AfAu8DNlCf6zUCMztJ0tB4lQBmK/Fcs4c2PhN0nmu57zOzxlZVW2gyp3Y7lUPg34B70znrgI8DfyJpB9Xq7V+ewnEGQRA0olHwrZldxdhyHw9TleQNgiCYMRyK0kNNK2nONuK5Zg9tfCZo73P1hNrkyg2CIIgikUEQtIppFWqS1kjanvJFxyyhNRuQtEzSDyVtS/mwH039jdY9nOlImiPpLknfS9uzPsdX0nxJN0p6QNL9kn55tn9ekY9dz7QJNUlzgC8C7wLOAC6RdMZ03b+PHAT+1MzOAM4HrkjP0XXdw1nCR4H73XYbcnw/D/yjmb0ZOIfq+Wbt5xX52OMznZraecAOM3vYzF6him9bO4337wtmttvM/i21n6P6A1lCs3UPZzSSlgK/AVyXtsUsz/GVNA/4VVLIkZm9YmbPMPs/r8jHrmE6hVrdmn+zFknLgbcAt9Ns3cOZzl8DH6OqBw9V/GGjHN8ZzArgKeDvk1l9naRjmcWf12TzsdtOOAomiKTjgG8DfzxqFWnSKtOzyq0s6d3Ak2a25VCPpc/MBX4R+JKZvQV4gVGm5mz7vCabj912plOo7QKWue1G+aIzEUmHUwm0b5jZd1L3nrTeIeOseziTeSvwHkmPUk0NvJ1qLmp+MnFgdn5mO4GdKTMGKvPsF5ndn1cnH9vMfg6MyMdOx8zGz6ovTKdQuxNYmTw0R1BNbG6cxvv3hTTP9GXgfjP7K7erybqHMxYzu9LMlqbcwYuBW83sNxnO8YXZ+VxPAI+nZdagyl3exuz+vDr52On7mJ9pVn9W/WJag28l/TrVvM0cYL2ZXTNtN+8Tkn4F+H9UebB57unPqObVbgD+A2ndQ7em4axC0tuA/25m75b0BirNbQFwF/AhM3v5UI6vVyStonJ+HEGV3ve7VP/QZ+3nJekvgA9SeePvolqGbgmz/LPqB5FREARBqwhHQRAErSKEWhAErSKEWhAErSKEWhAErSKEWhAErSKEWhAErSKEWhAErSKEWhAEreL/A6Uf7PshhFStAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAATUAAAD8CAYAAAAWjzPMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztnXuw3VWV5z/LhGd45UUIIVNJIIU8qggQQFC6bHkYaBEtFYk9DDow6a7RUbt7SqG7Shx6qEKrR9upsa3OYLppHzwGTIkSxYik2ikkJMhDQswkBuwkJCQh4amAgTV/nN/+3XVufr+c38k99yZ3n++n6tbZZ/9ee+9zzrpr7bX3WubuCCFELrxtXzdACCF6iYSaECIrJNSEEFkhoSaEyAoJNSFEVkioCSGyQkJNCJEVQxJqZjbPzNaY2Tozu65XjRJCiL3F9nbxrZmNAf4fcBGwEVgBzHf3p3rXPCGE6I6xQ7j2bGCdu68HMLPbgcuBWqF2+OGH++TJk4fwSCHEnti2bRsvv/yypffz5s3z7du3N7r2kUceuc/d5w1b40aIoQi1acCG8H4jcM6eLpg8eTI33njjEB4phNgTX/jCF9reb9++nZUrVza61swmDUebRpqhCLVGmNkCYAHAxIkTh/txQohB9Nv+7qEItU3A9PD+uKKuDXdfCCwEmDVrVn+NrhD7AW+99da+bsKIMhTv5wpgtpnNNLMDgSuBe3rTLCFEL3D3xn+5sNeamrvvMrNPAfcBY4BF7r6qZy0TQvSEnARWE4Y0p+buS4AlPWqLEGIYkFATQmSFhJoQIisk1IQQ2eDufef9lFATInOkqQkhskJCTQiRFRJqQohsyG1hbRMk1ITIHDkKhBBZIU1NCJENMj+FENkhoSaEyAoJNSFEVkioCSGyQdukhBDZIU1NCJEVEmpCiKyQUBNCZEW/CbWOiVfMbJGZbTWzJ0PdBDNbamZri9fxw9tMIcTekBwFTf5yoUk2qX8GBmdtvg64391nA/cX74UQ+yH9lk2qo1Bz938Fdgyqvhy4tSjfCnygx+0SQvSIfhNqezunNsXdNxflLcCUHrVHCNFjchJYTRiyo8Dd3cxqR83MFgALACZOnDjUxwkhuiA3LawJe5uh/TkzmwpQvG6tO9HdF7r7XHefe8QRR+zl44QQe0uvzE8zm2dma8xsnZntNo9uZgeZ2R3F8eVmNiMcu76oX2Nm7w31z5jZr8zsMTNb2Yv+7q1Quwe4uihfDXy/F40RQvSeXng/zWwM8HXgEuBkYL6ZnTzotGuAne5+AvBV4EvFtScDVwKn0HI6/kNxv8Qfu/scd5/bi/42WdJxG/AL4EQz22hm1wA3AxeZ2VrgwuK9EGI/pEea2tnAOndf7+5vALfTchhGogPxLuACM7Oi/nZ3f93dnwbWFfcbFjrOqbn7/JpDF/S4LUKIHtPDObVpwIbwfiNwTt057r7LzF4EJhb1Dw26dlpqIvCTYl7+H9194VAbqh0FQmROF0Jt0qB5rYW9EDIdeJe7bzKzo4GlZvbrYhnZXiOhJkTmdCHUtu9hXmsTMD28P66oqzpno5mNBY4Ent/Tte6eXrea2WJaZumQhNreOgqEEKOEHs2prQBmm9lMMzuQ1sT/PYPOiQ7EDwM/89aN7wGuLLyjM4HZwMNmNs7MDgcws3HAxcCTDBFpakJkTK+CRBZzZJ8C7gPGAIvcfZWZ3QisdPd7gG8C3zKzdbR2IV1ZXLvKzO4EngJ2AZ909zfNbAqwuOVLYCzwXXf/8VDbKqEmROb0avGtuy8Blgyq+0IovwZ8pObam4CbBtWtB07rSeMCEmpCZE6/7SiQUBMicyTUhBBZIaEmhMgGZZMSIvC2t7VW/PTbjyI3pKkJIbJCQk0IkRUSaiIbkvlYR7HosSNjxozZ4/FOPxqZr/uOfgwSKaEmROZIqIlRR9TIqrSvWFdVbqqxDSb9WOKPpqocNb14XBrcyNBv4yyhJkTmSFMTQmSD5tTEfk8yNetMymiKpnJVXZN7VVFlXkbzpqpcdzw+q+peojf0m1BrkqNgupk9YGZPmdkqM/tMUT/BzJaa2dridfzwN1cI0S39lsy4SZDIXcBfufvJwDuATxbZYa4D7nf32cD9xXshxH5Gvwm1JolXNgObi/LLZraaVtKEy4F3F6fdCiwDPj8srexD6jyaqRw9irE8duzY3cqxruq6Ts+q824mU/HNN98s62J5165dba+Dy1XXVZmk8VmiO7T3swNFctLTgeXAlELgAWwBptRcowztQuxDctLCmtA4R4GZHQbcDXzW3V+Kx4o45JUjpwztQuxbZH5WYGYH0BJo33H37xXVz5nZVHffbGZTga3D1ch+opPHMpmSBxxwQFl34IEHluWDDjpot/p4vMoUjSZfpy1R0WRMP4QqkxPgjTfeaHsFeP3113c7DvCHP/xht+urzKZ+M6V6QU4CqwlNvJ9GK6HCanf/SjgUM8dcDXy/980TQgwVaWq7807gKuBXZvZYUffXwM3AnWZ2DfBb4IrhaWL+RE0saUpRY4qa1sEHHwzAIYccUtbFcjoOA1pb1M46bXLv9OXutI4talJJ64ra2WuvvVaWf//73+9WjsejJhe1wapniWrkKKjA3f8vUPdNvqC3zRFC9JqctLAmaEeBEJkjoSaGlaqtSdEkTA6AaEYeeuihZfmwww5re4V28zM6EKrWkUVTMJWjmVc1UV/lqIABszg6J2I5tWvcuHFlXXIIQLv5+corr7S9Avzud78ry9EsTWibVTMk1IQQ2ZCbE6AJEmpCZI6Emug5VduQonczmozJ1Dz88MPLurho+cgjjwTazbz4pX311VfL8s6dOwHYvn17Wbdt27ay/MILLwDtJl8085LZGtsazeJkAh911FFl3eTJk8vypEmTABg/fiDWQTRFo9mc6uP9X3ppYI33yy+/DLSbpNGUTW2NYy1TtEW/jYOEmhCZI01N9JyqFftx7Vl0BCStZ8KECWVd1NTSRH3UyJ599tmy/Jvf/KYsr1u3DoBnnnmm8tznn38eaNfUqrSfOq0yaWpxT++xxx5blmfMmAHACSecUNYdf/zxlecmTS1qoFFrS/2OYxm1tuTsqFrP1s9oTk0IkR0SakKIrJBQEz2hbkN6Mt+iyRkn0tPkepx8jzz33HMAPPnkk2XdihUryvKjjz5alteuXQvA1q0DsQaiyZboFA68Lhx31Y8l9uvoo48GYPbs2WXd6aefXpbPOuussnzqqacCMGXKQASraNZWbb6PdMps1W+T5REJNSFENmjvpxAiO6SpiSGRTLZO69CieZlMThgwReN/1+i9/PnPfw7A0qVLy7qHH364LMd1aInoRUweSRgw9aKnNa6PS32IHsW0Xgxgx44dwIBJDLBly5bd2h3b/9hjj5XlVatWleWLLroIgPPPP7+yrdFET1SZmp1+wP2mtYCEmhAiMyTURNdUTa7XrcJPmlDUjqLWljSJtMYM4Ec/+lFZXrx4MQCPP/54ZVvivebMmQO0T86fdNJJZXn69OlA+4R8XOWf1obFTe5xE3pa57Zhw4aybvXq1WU5OS2idhY1yXvvvbcsb9y4ERjY5QBwySWXlOW01i32ryribqzr5NToF62tV0LNzOYBXwPGALe4+82Djh8E/AtwJvA88FF3f6Y4dj1wDfAm8Gl3v6/JPfeGxjkKhBCjj+QoaPK3J8xsDPB14BLgZGB+kSozcg2w091PAL4KfKm49mTgSuAUYB7wD2Y2puE9u0ZCTYjM6VE477OBde6+3t3fAG6nlSYzcjmtdJkAdwEXFOkALgdud/fX3f1pYF1xvyb37BqZnz2gap1X3TaotOWpLrNWmlSPJue3v/3tshzN0sQ555xTli+99NKyfN555wHt68TihHtsY6Jq8r0qF2gkxmNLm+hhYJ3cgw8+WNYtWbKkLC9fvrwsJ3M6bv+KXHbZZQDMnDmzrItjmDbid4oN12/zS9BVnyeZ2crwfqG7LyzK04AN4dhG4BzaKc9x911m9iIwsah/aNC104pyp3t2TUehZmYHA/8KHFScf5e732BmM2lJ1onAI8BVhbQVQuxHdCHUtrv73OFsy0jQxPx8HXiPu58GzAHmmdk7aNnLXy3s55207GkhxH5EU9OzgeDbBEwP748r6irPMbOxwJG0HAZ11za5Z9c0SbziQArjcEDx58B7gI8V9bcCXwS+MdQGjRaqYqTBgMewKu4YDMRDi2Gx4zqvtA4teTmh2uT80Ic+VJbnz59fls8444yynDyFMfJGfFZacxZNvqoMTnWZrVJkjbi2LfUP4MwzzwTaI3OccsopZfm2224ry3fffTfQ3tc4BqkvVduw4nPrslElUzROiPdL7LUemdwrgNmFhbaJ1sT/xwadk9Jm/gL4MPAzd3czuwf4rpl9BTgWmA08TCuhU6d7dk3TZMZjaJmYJ9DyVvwGeMHd06RFtJEHX7sAWADtSweEECNDLwR2MUf2KeA+WssvFrn7KjO7EVjp7vfQyg/8LTNbB+ygJaQozrsTeArYBXzS3d8EqLrnUNvaSKgVDZhjZkcBi4G3N31AMdG4EGDWrFnZzNLWZTVPuwfieq+oqaV4YVE7ipvT006BunVoSUNbsGBBWZc2g0O7VrZ+/XqgPfJt3BHQKQZZlaMgUhUbLmptaadE3DFx7rnnVp6bSBobtI9BGpdp0wb+d0YHSRrXONZxTV1KMlOVYT5nehlPzd2XAEsG1X0hlF8DPlJz7U3ATU3uOVS68n66+wtm9gBwLnCUmY0ttLWe2MJCiN7TD8I70tFRYGaTCw0NMzsEuAhYDTxAy26Glh39/eFqpBBi7+mRo2DU0ERTmwrcWsyrvQ24091/aGZPAbeb2X8HHqVlT2dNnXOgaiI9mp+xnL48Max2jIcWN6cnopmVnALR5IwJSjZv3lyWk4kbTdKqvJ2xX92Q5mrinE2M15baEts3derUshz7kPqVtktB+zq2NC7R0ZC2ecHA5ve6cU/tqgpXDnk7DXISWE1o4v18Aji9on49rRXBQoj9FMVTE0JkhzQ1UUtd2Oto0iUvXFynFuOpJZMwZn2KIbhTFIsYjSJufUrr0KIZFU3O5OUb3K4q0j2qthPVUWWCx/5F0vNjm2Jbo3mY+hX7umbNmrKcxiWO1WmnnVaWU77ROO6xnD6X2JbY75x/+Dn3rQoJNSEyR0JNCJEVEmpiN6q8g3XmZ/IoxqS80aRLUSzidqAUzSKSAjzCQLQNGDBL08JaaF/IG9uSFp/G7UKRZH7Ga6JJVhUksurcOvOzyhMc2xoXBc+aNQto7+sDDzxQlpctWwa0j1UcwxNPPBFoz0YVP4PUltj++BlWeUJzmGDPbblGEyTUhMicHIRzN0iodUETR0GVRhC1gKSdxGQkMS9nmtyOIbhjPLSkXUUtJzoNIklDi22pWl8X19nFc1N9VV2sj2MRtYL0/KjJxbbGPqQ1Z3U5Qh96qBWOK45VHMN0r7glq9PnUrfuMDfNJrf+dEJCTYjMkVATQmSD5tTEHqkzPzuZZHFtVFpvFbdJxa1FabtPzPoUQ3CneGgx2kan9Wix3dFpkL7sdVuHqvJ+xr5W3b8qRHjdubEPL774ItA+0R/H4JhjjgHaTc44hmlcY7jv6ChIYxTbX2d+5oaEmhAiKyTUhBBZIe+n2I1kmjQxP6vMs2h+pmS9KRHwYJL5FSNQRJMumWzRjIzHoymZyrFN8b92MsliXTw39beqT/F43Tar9PzYpniv2IfUr9jvWE7jEs3POIZpXONYR/OzytNbZ37mZIpqTk0IkR0SaqKWuv/mnSaco0byyiuvtL0OvmbChAlAfT6HtCK/LhZY1er/qJ1ErSl92evCeXeiSoOrcprEH1XdKv6qfJ9xDNK4xLGKY5jKdbsn0nX9oJ0NRkJNCJEVEmpCiGxQkMg9UITzXglscvf39XuG9jrTJZlf8b9jNAlTXspoBkaTKGVYqgoBDvXmVSJ+gdNz67ZsdTIf073qTMaquliuen4dqV+xr3EM0rjEtsQxTONaFyMt9atfTM5Iv2lq3QSn/wythCsJZWgXYhTQb4lXGgk1MzsO+BPgluK90crQfldxyq3AB4ajgUKIodFvQq2p+fn3wOeAlH12Ig0ztI9Wusmw1MmMiSZhMs/qvJdVkTHiFy5d1+RLWOX9rFqHVtfWTnMxVevY6iJ+VNGpX1XRQepM4VRu2uYm5JJhKieB1YQmeT/fB2x190f25gFmtsDMVprZypgqTQgx/DTV0nISfE00tXcC7zezS4GDgSOAr9EwQ7u7LwQWAsyaNWvUjFz8z1y3oj7R6QtRpYlVTcjDgMYRJ7yrHBFNNI69cRR0ivIbr69a59aNo6BTv+IYVGliVRpiJw27mx/vaNbOIrn0oykdNTV3v97dj3P3GcCVwM/c/U9RhnYhRgX9pqntXWruFp8H/tLM1tGaY8s+Q7sQo5F+E2pdLb5192XAsqLc1xna45egasI7mlHR/ErhumOI62gepI3dKWnK4Ht1ildWFWa804b7SKd1aFV0chQ0cbqkfsW+xjFI4xLHKo5hGte60OVVjoicfsh15CawmqAdBUJkjoSaECIrJNRELXWmSzSJqr5A0WQ87LDD2l4HX7Njxw6gPt7auHHjgHrvaZX3sJP3sVO8tLofRZUHtsr7WZcDNT439SsSxyCNS2xLHMNUrjPP03V1n1XOP/x+835KqAmRMZpTE5WkL0Wddla3TisRNZWUYb0uXlpKrLJhw4ayLm5iTxu7o0ZSN3meynUxxOI6sKr2703ilfis9PzYpvjM2IfUr9jXOAZpXCJxDNO4xrGOVO3k6BetLbf+dGIoSzqEEKOAkVjSYWYTzGypma0tXsfXnHd1cc5aM7s61J9pZr8ys3Vm9j+L/eWY2RfNbJOZPVb8XdqpLRJqQmTOCK1Tuw64391nA/cX79swswnADcA5tJaD3RCE3zeA/wTMLv7mhUu/6u5zir8lnRoi87MLmpifybyKx6NJNHnyZACOPfbYsu7QQw8ty1u2bAFg9eqBKE87d+4sy0ceeSQwYK5Be97QTu2uMlujSdmN+Znqo8nZ6ccRj8c+pH7FvsYxSOMSxyqOYRrXOqdE+lyamJ854SMXJPJy4N1F+VZa61k/P+ic9wJL3X0HgJktBeaZ2TLgCHd/qKj/F1pRf360Nw2RpiZE5oyQpjbF3TcX5S3AlIpzpgEbwvsU3WdaUR5cn/iUmT1hZovqzNqIhJoQmdOFUJuUIuoUfwvifczsp2b2ZMXf5YOe50CvVN9vAMcDc4DNwP/odIHMzy6oMz+jRy9572JdDEs9adIkAGbMmFHWHX300WU55bV89NFHy7q1a9eW5TPPPLPtPgCbN28uy3FNWjI160KAp3DY8ZqqbFRVdbE+ejcj6fl16+RiH9I9nnjiibIujkEK1x3HLZbTveryilZ9Lv1gfkJXfdvu7nP3cJ8L646Z2XNmNtXdN5vZVGBrxWmbGDBRoRXdZ1lRf9yg+k3FM0u3t5n9b+CHnTohTU2IzBkh8/MeWtF6oD5qz33AxWY2vjAjLwbuK8zWl8zsHYXX8z+k6wsBmfgg8GSnhkhTa0DVhHonTS1mCo+r5cePb00JnHDCCWXd7Nmzy3LS1B577LGy7sEHHyzLxx9/PNCu5cTgm/G5SUOs06SSpha/0FVaVayLToE0KV93/6p7xbGIfUgZ1mNf4xgk4ljFMUzjGtsXx6KTphbJaQX+CC6+vRm408yuAX4LXAFgZnOBP3f3a919h5n9LbCiuObG5DQA/jPwz8AhtBwEyUnwZTObQ8ucfQb4s04NkVATInNGQki7+/PABRX1K4Frw/tFwKKa806tqL+q27ZIqAmROTnPF1YhodYFTRwFyeRJE9vQnp8ymV/JjAQ4/fTTy3IyubZt21bWLVkysN7wlFNOAeDcc88t66ZOHZh2iE6DV199dbfnV4XQ7ibJTCSNQV2uzUQ0OWNbo9n6i1/8AmjvazJJYWAdWhyrOIbpGXHcYzl9LnIU5I+EmhAZow3tQojskFCrwMyeAV4G3gR2ufvcYh/XHcAMWl6JK9x9Z909cqAublncepO8bDEUdSwnj2Tc4nPWWWeV5VWrVgFw7733lnXLly8vy7fddhvQvsXo1FMH5lfjmrjt27e3vcJAWOzY1rpw3enHUJe5KpmvVdE2YMC7WbUeDeDJJwe886lfsa+Rs89uRY6PYxXHMLWxbtyr+lo3nZAb/SbUuplM+eNiQ2lanNdxA6sQYt/z1ltvNfrLhaEsvr2c1sZVitcPDL05Qohe0nThbU7aXNM5NQd+YmYO/KO3EhQ32cCaLfFLEE2a5GmMps8rr7xSlpOXLnoEo/l40UUXAbBx48D+3scff7ws33333bu1Zf78+WX5jDPOKMuzZs0CYPr06WXdiy++WJaTKZq8pNC+tSj1Ky46jqZm6kNVtA0YMDWjFzN5OWHA5Kzr12mnnVaW07jEsYpjmPoQxzp+BulzqTM/c6Zf+ploKtTe5e6bzOxoYKmZ/ToedHcvBN5uFJtiF0B9tFchxPAhoVaBu6fNpVvNbDGtAG9NNrBSaHULAWbNmpXN6HZKdhLXSEXtIeWnjHG/pkwZUHLPP/98oF27iZrUunXrgHbNJmp1l146EBj0vPPOA9q3FsVnRQ0uURXWOvavymkQtbsYDy1tTo9bn+I6tCqnQNz69MEPfrAsp3GJ7Y/r75IGGsc6fgbpc+kX50Ck34Raxzk1MxtnZoenMq1NqE/SbAOrEGIfkoJE9pOjoImmNgVYXPyHHgt8191/bGYrqNjAKoTYv+g3Ta2jUHP39cBpFfWVG1j7kSqTJppkMdx2iqiRzFBon2tMMcIuueSSymctXrwYaHceRDNuzZo1ZfmBBx4A2rcWnXTSSWU5mZ/x+XGdW1U8tTj5nvJyxqxPMQR3iocWo21EszqSnALR5IxjEGOnJWJ0klSOYx0/g/S59NsPHPqvz9pRIETmSKgJIbJCQk10TdUka1wPFb1waW1YDKAY14FVBZG87LLLynJK2rt06dKy7uGHHy7LMbrHsmXLAHjooYfKumOOOaYsJ0/ihAkTyrq45qwqm1TcZrVjRyu+X0w0nLI+QXu/EynaBgxsfYKBdWjJywntJmeKJBK9q+n5sV3xmVWZo3KaEG9CbgtrmyChJkTm9Jsgl1DrMZ2+QGkiu9Par6SRAcycObMsp7yX06YNZBBLMdagOmHL1q0DSwhTuPBYjs+PsdVSOfapUwyymJczaVpxnVx0WsTN6WmnQFyHFkkaWtycH50OaVzj2rW6HJ/9hjQ1IURWSKgJIbJBc2qiZ9Rto0rmUVxPFUlfwGg6HXHEEWU55Qg955xzyrq43SluAk9bqqLJ+eyzz5bltM4sbi2qMt+iIyPGQzvssMOA9nVuMcZZMj+j0yOG4I7nps3p8flxHVpyClSZnPG6OlO5n5FQE0JkRb8Jdwk1ITJG5qcYFqpir8UtPFWRMeLWpLj2KsUri1E+4nquuA7sxBNPBNo9hnEdWzLl6iJbVJmfcXtXMj+jpzY+P4XxTmvvoD0GWjTLUySSGO8tmp9pHVqVyRnb2m8/4Cb025hIqAmRORJqoufEOY209qsuAmuVoyBqdUmTSloStG9Cj5pUWvMVE5/ENW8pF2a8f1VezLh2Le6ESFFwo9YYy1V5RaMmWBUdOGqNUStL10XtrGrNXL/NHzVBQk0IkRUSakKIbEhBIvsJCbURpuoLVvWfNJ4XTcKqvKJ15mcyBaPJWGcq7qktkbocoFXtTm1NZi7Um5+pHI9XJYHR1qfukaYmhMgKCTUhRFZIqFVgZkcBtwCn0soB+h+BNcAdwAzgGeAKd99ZcwuxBzqZpFWmaDTpopcwmpTJOxlzdVbFcYsmZVyTVkWV1zbWVZnK0YyM7Y71yatZ5X2NZZmc3dGPi2+bZmj/GvBjd387rXwFq4HrgPvdfTZwf/FeCLGfoQztgzCzI4E/Aj4O4O5vAG+Y2eXAu4vTbgWWAZ8fjkb2E021tqqs8NCuCSWtrC7KbtU6sqo4b1Xr6OraUqW1Re0rlquu6/Qs0T39NnZNNLWZwDbgn8zsUTO7pcj/OcXdNxfnbKGVSk8IsZ/Rb5paE6E2FjgD+Ia7nw68yiBT01sjUjkqZrbAzFaa2cq4l08IMfw0FWhDFWpmNsHMlprZ2uJ1fM15VxfnrDWzq0P9TWa2wcxeGXT+QWZ2h5mtM7PlZjajU1uaOAo2AhvdPSWXvIuWUHvOzKa6+2YzmwpsrbrY3RcCCwFmzZqVz7+DYabOZEimYtUmeWg379K5VSG6Y7lTaPFIlXlYF8OsanK/Lhy4tjkNHyOkhaU59pvN7Lrifdt0lJlNAG4A5tJSgh4xs3sKB+MPgP8FrB1032uAne5+gpldCXwJ+OieGtJRU3P3LcAGMzuxqLoAeAq4B0iS9mrg+53uJYQYeUbI/Lyc1tw6xesHKs55L7DU3XcUgmwpMK9o40NhOqvuvncBF1iH/7xN16n9F+A7ZnYgsB74BC2BeKeZXQP8Frii4b2EECPICGm/TebYpwEbwvuNRd2eKK9x911m9iIwEdhed0Ejoebuj9FSGQdzQZPrRe+o+oJGk7LKLK0zL1O5k8lZR3pWncey6r+/PJojS5da2CQzWxneLyymjwAws58Cx+x+GX8z6JluZvtsqkk7CoTInC6E2nZ3r1Je0n0urDtmZk3m2DcxsAwM4DhaS8H2xCZgOrDRzMYCRwLP7+mCpotvhRCjlBGaU2syx34fcLGZjS+8oxcXdU3v+2HgZ96hsdLUMqAbT2lkb83OwXT6QcjM3LeMkPfzZirm2M1sLvDn7n6tu+8ws78FVhTX3OjuO4rzvgx8DDjUzDYCt7j7F4FvAt8ys3XADuDKTg2RUBMic0ZCqLn781TMsbv7SuDa8H4RsKjivM8Bn6uofw34SDdtkVDLmKFqSEnTk6Y1elGQSCFEduS0BaoJEmpCZI6EmhAF/Wa25IqEmhAiG3KLwNEECTUhMkdCTQiRFf02jSChJkTmSFMTQmSD5tSEENkhoSaEyAoJNSFEVshRIITIBs2pCSGyo9+EWscgkWZ2opk9Fv5oTVZpAAAGuElEQVReMrPPNk2JJYTYtyjv5yDcfY27z3H3OcCZwO+AxQykxJoN3M+gXKBCiP0DCbU9cwHwG3f/Lc1SYgkh9jH9JtS6nVO7EritKDdJiSWE2If0Y5DIxppakfPz/cD/GXysSIRQKerNbIGZrTSzlS+99NJeN1QIsXf0m6bWjfl5CfBLd3+ueP9ckQqLPaTEwt0Xuvtcd597xBFHDK21QoiukVCrZz4Dpic0S4klhNjHSKhVYGbjgIuA74Xqm4GLzGwtcGHxXgixH9FUoOUk1Bo5Ctz9VWDioLrKlFhCiP2LnARWE7SjQIjM6Tfvp4SaEJkjTU0IkQ25zZc1QUJNiMyRUBNCZIWEmhAiK+QoEEJkg+bUhBDZIaEmhMgKCTUhRFZIqAkhskJCTQiRDf0YJFJCTYjM6TdNrdscBUKIUcZIhB5qml3OzK4uzllrZleH+pvMbIOZvTLo/I+b2baQze7aTm2RUBMic0YonlrH7HJmNgG4ATgHOBu4IQi/HxR1VdyRMtq5+y2dGiKhJkTGjGCQyCbZ5d4LLHX3He6+E1gKzCva+VBI5DQkJNSEyJwuhNqklCSp+FvQxWOaZJebBmwI7zcWdZ34kJk9YWZ3mdn0TifLUSBE5nTh/dzu7nPrDprZT4FjKg79TXzj7m5mvfJO/AC4zd1fN7M/o6UFvmdPFzQSamb2F8C1tNLg/Qr4BDAVuJ1WmO9HgKvc/Y29b7sQYjjolffT3S+sO2Zmz5nZVHffvIfscpuAd4f3xwHLOjzz+fD2FuDLndrZ0fw0s2nAp4G57n4qMIZWUuMvAV919xOAncA1ne4lhBhZRnBOrUl2ufuAi81sfOEguLioqyWl4Sx4P7C6U0OazqmNBQ4xs7HAocBmWirgXcXxuolBIcQ+ZoSEWmV2OTOba2a3FO3YAfwtsKL4u7Gow8y+bGYbgUPNbKOZfbG476fNbJWZPU5Lufp4p4Z0ND/dfZOZ/R3wb8DvgZ/QMjdfcPddxWm1E37FZOMCgIkTJ1adIoQYRkZi8W1ddjl3X0lr6iq9XwQsqjjvc8DnKuqvB67vpi1NzM/xtNy1M4FjgXEUbtgmuDK0C7FPeeuttxr95UITR8GFwNPuvg3AzL4HvBM4yszGFtracbQmAffI008/vf2qq656Fdg+hDbvr0xC/Rot5NgnaPVr3KC6+4r6JmQxJk2E2r8B7zCzQ2mZnxcAK4EHgA/T8oDWTQy24e6TzWzlntzGoxX1a/SQY5+g7NeMWOfuja2qXOhofrr7cloOgV/SWs7xNmAh8HngL81sHa1lHd8cxnYKIUQjGq1Tc/cbaO3Ziqynfq+WEELsE/bFNqmF++CZI4H6NXrIsU+Qb7+6wvot1pIQIm+0oV0IkRUjKtTMbJ6ZrTGzdWa2W7yl0YCZTTezB8zsqWKl82eK+kZB8vZ3zGyMmT1qZj8s3s80s+XFZ3aHmR24r9vYLWZ2VBHh4ddmttrMzh3tn5eZ/UXx/XvSzG4zs4Nz+Kx6wYgJNTMbA3wduAQ4GZhvZieP1PN7yC7gr9z9ZOAdwCeLfnQMkjdK+Azt++ty2OP7NeDH7v524DRa/Ru1n5f2Y++ZkdTUzgbWufv6IprH7bR2Kowq3H2zu/+yKL9M6wcyjWZB8vZrzOw44E9oRUPAzIxRvsfXzI4E/ohiyZG7v+HuLzD6Py/tx65hJIXa3gaI228xsxnA6cBymgXJ29/5e1r779KemYk03OO7HzMT2Ab8U2FW32Jm4xjFn5e7bwLSfuzNwIt0sR87d+Qo2EvM7DDgbuCz7v5SPOYtl/Kociub2fuAre7+yL5uS48ZC5wBfMPdTwdeZZCpOdo+r6Hux86dkRRqm4AYirfRftH9ETM7gJZA+467f6+ofi7FftpDkLz9mXcC7zezZ2hNDbyH1lzUUYWJA6PzM9sIbCx2xkDLPDuD0f15lfux3f0PQNt+7OKc0fhZ9YSRFGorgNmFh+ZAWhOb94zg83tCMc/0TWC1u38lHGoSJG+/xd2vd/fjir2DVwI/c/c/ZWCPL4zOfm0BNpjZiUXVBcBTjO7Pq9yPXXwfU59G9WfVK0Z08a2ZXUpr3mYMsMjdbxqxh/cIM3sX8HNa+2DT3NNf05pXuxP4d8BvgStSALzRhpm9G/iv7v4+M5tFS3ObADwK/Ht3f31ftq9bzGwOLefHgbS2932C1j/0Uft5mdl/Az5Kyxv/KK2YZdMY5Z9VL9COAiFEVshRIITICgk1IURWSKgJIbJCQk0IkRUSakKIrJBQE0JkhYSaECIrJNSEEFnx/wG9+beSYxkc6wAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "#NBVAL_SKIP\n", "\n", "# Let's see what we got....\n", "plot_image(v[0].data[0])\n", "plot_image(v[1].data[0])\n", "plot_image(p.data[0])" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [], "source": [ "norm_p = norm(p)\n", "assert np.isclose(norm_p, .35098, atol=1e-4, rtol=0)" ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "scrolled": true }, "outputs": [], "source": [ "# # 4th order acoustic according to fdelmoc\n", "p4 = TimeFunction(name='p', grid=grid, staggered=NODE, space_order=4, time_order=1)\n", "v4 = VectorTimeFunction(name='v', grid=grid, space_order=4, time_order=1)\n", "src_p = src.inject(field=p4.forward, expr=src)\n", "u_v_4 = Eq(v4.forward, v4 + ro * grad(p4))\n", "u_p_4 = Eq(p4.forward, p4 + l2m * div(v4.forward))" ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "scrolled": true }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "Operator `Kernel` run in 0.01 s\n" ] } ], "source": [ "#NBVAL_IGNORE_OUTPUT\n", "op_4 = Operator([u_v_4, u_p_4] + src_p)\n", "# Propagate the source\n", "op_4(time=src.time_range.num-1)" ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "scrolled": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAATUAAAD8CAYAAAAWjzPMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztnX+wXVWV5z9fEn4LCSSYQBJMaOKPQCm0KaRLq4sWhWhbpmfGH8F2pLvR1FTjtD3dUyptlYzYVGn1KGqr9KQw3fijDQwtZcpBYxQoZ6oUeBhaCSEQgpBECIT84HcwZM0fZ+/71rvvnNxz373v5b1z16fq1dtnn1/7nHvfemvttddaMjOCIAiawhGHewBBEAT9JIRaEASNIoRaEASNIoRaEASNIoRaEASNIoRaEASNIoRaEAS1kLRM0mZJWyR9smT/0ZJuSPvvkLTQ7bsi9W+WdHHbedMkbZD0g36Msyeh1ukhgyBoBpKmAV8D3gEsAS6RtKTtsMuAPWZ2JnAN8Pl07hJgBXAWsAz4erpe5mPApn6NdcxCreZDBkHQDM4DtpjZVjN7CVgDLG87ZjlwfWrfBFwoSal/jZntN7OHgS3pekiaD/wxcF2/Bjq9h3NbDwkgKT/kfVUnHHvssTZjxowebhkEwaHYt28fL7zwgvL2smXLbNeuXbXOvfvuu9eZ2bKK3fOAbW57O/CmqmPM7ICkfcCs1P+LtnPnpfaXgI8DJ9QaZA16EWp1HnIEM2bM4IMf/GAPtwyC4FB8+9vfHrG9a9cuhoaGap0r6bWS/MGrzGxVH4fXfr93AU+Y2d2SLujXdXsRarWQtBJYCXDCCX0TxkEQ1KSL+O5dZra0Yt8OYIHbnp/6yo7ZLmk6MAN46hDnvht4t6R3AscAJ0r6tpn1pPn04iio85CY2SozW2pmS4877rgebhcEwVg4ePBgrZ8O3AUslrRI0lEUE/9r245ZC1ya2u8BbrVCoq4FViTv6CJgMXCnmV1hZvPNbGG63q29CjToTVNrPSSFMFsBfKDXAQVB0D/MrBtN7VDXOSDpo8A6YBqw2sw2SroKGDKztcA3gG9J2gLsppAJpONupJhvPwBcbmYv9zyoCsYs1Koesm8jC4KgL/QrvZiZ3QLc0tb3add+EXhvxblXA1cf4tq3A7f3Y5w9zamVPWQQBJOLQcuZOO6OgiAIDi8h1IIgaBQh1IIgaAxmVsez2ShCqAVBwwlNLQiCRhFCLQiCRhFCLQiCxtCvxbdTiRBqQdBwwlEQBEGjCE0tCILGEOZnEASNI4RaEASNIoRaEASNIoRaEASNIcKkgiBoHKGpBUHQKEKoBUHQKEKoBUHQKAZNqHWsJiVptaQnJN3r+k6WtF7Sg+n3SeM7zCAIxkJ2FPShmtSUoU6JvH8B2qs2fxL4qZktBn6atoMgmITkqIJOP02ho1Azs59RlLvyLAeuT+3rgT/p87iCIOgT/RJqkpZJ2ixpi6RRikyq63lD2n+HpIVu3xWpf7Oki1PfMZLulPTvkjZK+kw/nnesxYznmNljqf04MKcfgwmCoP/0Q6hJmgZ8DXgHsAS4RNKStsMuA/aY2ZnANcDn07lLKGqAnkVh9X09XW8/8FYzewNwDrBM0vm9Pm8vFdoBSBWYK9+IpJWShiQNPf/8873eLgiCLqgr0GpoaucBW8xsq5m9BKyhsNg83oK7CbhQklL/GjPbb2YPA1uA86zg2XT8kemnZzt4rEJtp6RTAdLvJ6oONLNVZrbUzJYed9xxY7xdEARjpQuhNjsrIOlnpbvMPGCb296e+ig7xswOAPuAWYc6V9I0SfdQyJD1ZnZHr8871iUda4FLgc+l39/vdSBBEIwPXXg2d5nZ0vEcSztm9jJwjqSZwM2SzjazezuddyjqLOn4LvBz4DWStku6jEKYvV3Sg8Db0nYQBJOQPpmfO4AFbnt+6is9RtJ0YAbwVJ1zzWwvcBujV1p0TUdNzcwuqdh1Ya83D4JgfOnjco27gMWSFlEIpBXAB9qOyRbcz4H3ALeamUlaC/yrpC8CpwGLgTslnQL8zsz2SjoWeDvJudALEVEQBA2nH0LNzA5I+iiwDpgGrDazjZKuAobMbC3wDeBbkrZQLANbkc7dKOlG4D7gAHC5mb2c5uOvT57QI4AbzewHvY41hFoQNJx+Law1s1uAW9r6Pu3aLwLvrTj3auDqtr5fAef2ZXCOEGpB0HCaFC1QhxBqQdBgIklkEASNIzS1IAgaRQi1IAgaRQi1IAgaRQi1IAgaQzgKgoGnmz+AI47oOclLMAGEphYEQaMIoRYEQaMIoRZMCYrce72Rv+xjnXOpe16YqYePptUfqEMItSBoOCHUgkmB18TK2lX7y/Bf6k7aVdmxNfLXl7azhlZ1z9DgJobwfgZB0ChCUwuCoDHEnFpwWCgz37xp5tvTp0+vPKedbHZ48+Pll19utQ8cOFB5jj/Wn+P/QHK7zOT07WnTppXuLzOLwiTtP4Mm1OrUKFgg6TZJ96WCox9L/SdLWi/pwfT7pPEfbhAE3RIV2kdzAPhbM1sCnA9cnoqTfhL4qZktBn6atoMgmGQMmlCrU3jlMeCx1H5G0iaKmn3LgQvSYdcDtwOfGJdRNpBO5uWRRx7Z6vPmW5n56b+Q3qTMZmM3X9gy72eVyVpm3vpx5XFXPUtZv79WmKK9E7GfHZC0kCKn+B3AnCTwAB4H5lScsxJYCXDCCSeMdZxBEIyRJmlhdaj9r1DSK4B/A/7azJ72+6x4a6VvLiq0B8HhJczPEiQdSSHQvmNm30vdOyWdamaPpVJXT4zXIKc6Zd5Bb4YdffTRo9pHHXVUqy+bnDD8X9ebgfv372+1ff9LL70EjDQfvSmSr1Vm/voxVnkv8339Pcs8qv58/1zHHHPMqH5///CO9od+CSxJy4AvU5TIu87MPte2/2jgm8AbKYoYv9/MfpP2XQFcBrwM/JWZrZO0IB0/h0IpWmVmX+51nHW8n6Ko57fJzL7oduXCpaTf3+91MEEQ9J9+aGqpNufXgHcAS4BLksPQcxmwx8zOBK4hFSZOx60AzqKowP71dL0qJ2RP1NHU3gz8Z+DXku5JfX8HfA64UdJlwCPA+3odTJMomzCH4cnxY489ttXn21l78RrJ7373u1b7+eefB+CZZ55p9T333HOt9osvvthqZ02n01j8hL3XpHx/xmt9WVPz9y8blx/TCy+8MOp8gDw1Uaa9QXnIVWhtnemjo+A8YIuZbQWQtIbCWXifO2Y58D9S+ybgq0kpWg6sMbP9wMOp2PF5ZvZzyp2Q/ppdU8f7+f+AquDCC3u5eRAE408X5udsSUNue5WZrUrtecA2t2878Ka281vHpIru+4BZqf8XbefO8ye2OSF7IiIKgqDhdCHUdpnZ0vEcSxmHckKOhRBqfSabelWOgGxm+eUt3vzMX0Bv0u3atavVfvLJJwF4+unhz77M5IRh883f35t3eYzejCs71veVOS38/fft2zdq3E899VSr79lnn221vfmZnRqveMUrWn3HH3/8qGfx7zVM0Xr0yVGwA1jgtuenvrJjtkuaDsygcBhUnlvhhOyJ+CYEQYOp6ySoIfjuAhZLWiTpKIqJ/7Vtx3jn4XuAW9Nyr7XACklHS1oELAbuPIQTsidCUwuChtMPTS3NkX0UWEexpGO1mW2UdBUwZGZrKQTUt5IjYDeF4CMddyOFA+AAcLmZvSzpLZQ4Ic3sll7GGkKtD5R5F73J5k2qGTNmAMNmKIz0KGZTbdu24TnZ7du3t9p79uwBRnpEvZfSXzePwZtsZWvOqsKs8j28F9K3s9ns73/qqae22rt37wZgx45hK8U/lzer9+7dO+q5vHmZ36F/r2WmaJiho+lXmFQSNre09X3atV8E3ltx7tXA1W19h3JCjpkQakHQcJoULVCHEGpjpCqHWNZkvHZ28sknt9q532skO3fubLXvv/9+AB5++OFWX9bOYFg78Y4G386aIAxrbWXrvTxeU/ST/nmMeRIfRk7un3jiiQCcdNJw1qmZM2e22rNmzQLglFNOKR3rAw880GpnbdQ7QPy4Mv69e60tNLRymhYCVYcQakHQcEKoBUHQKEKoBbXwZpCfKM9rq7xJls00GJ60/e1vf9vqu+eee1rtTZs2ASPNMD/5n025OXOGMz35e/k1b34ivf3+vl0V8J7NTh/a5NeZ5TH6NXV+XHPnzgVgwYLhJUrePC5bn/foo4+Ouj4Mm5fezCwzRWPt2mhCqAVB0BgiSWQQBI0jNLWgEm/u+HAhH3qUzStvZnkzMJud3uTcsGFDq53Xa3mPqTffctt7Gb0X0FNmXnZDWe427wnNntKcOaS9nb2nCxcubPX5dWz+vWW8d9Wvacvvxb/3srbvCwpCqAVB0ChCqAWV+Iln7xwoWyfmtZAchA6wceNGYKR29thjj7Xas2fPBmDRokWtPt/OGpq/f1Vdz9xfVquzvZ3plKXXXz+3vXZWtqatLMgeRmptGe+I8LnZcvRB1thgpIaar1vlHBhkp0EItSAIGkM4CoIgaByhqQWjKMuR5tdYefMzr1PzE95bt25ttbPZ6ddj+fPPOOMMYKTJmcONYHgivGy9GXRnfpbhzc9MlemWr1UVZuXzqGW8yeif+1WvehUw0vz0ZnvZmji/ji2b+94sH2ST0zNoQq1O4ZVjJN0p6d8lbZT0mdS/SNIdkrZIuiHlWAqCYJIxaCXy6vwr2w+81czeAJwDLJN0PkWlmGtS5Zg9FJVkgiCYRPQxSeSUoU7hFQOyTXBk+jHgrcAHUv/1FFVkru3/EA8/2YzxnjsfuuTNqGwe+hxi2eMJ8OCDDwIjTTa/Di2bn9kL2n7ffJ4/368jK6vrOdYvbDZFq8y4TmFY2RPq86Z5U9yvtXv9618PwKtf/epWn3+HOYvHI4880urz3tGc/cRPC5TVEB1Ek7RJAqsOdYsZTwPuBs6kqP33ELDXzPJf06jqMO7clcBKGPnHHwTBxBDezxLM7GXgHEkzgZuB19a9QSqxtQpg7ty5U+ZfRtl6La8x+aIgXjvIDgKvUeQgdRguTHLaaae1+n7v936v1c4B4X5C3X8ps1bmNbVO69D6+Z+6Ko9c2Vhz2ztNfO44r7XNm1f8T5w/f36r77WvHf6a5Txz/nwfaJ8dCD6PXVWV+0GiaaZlHbr6pM1sr6TbgD8AZkqanrS1ssoyQRBMAgZNqNXxfp6SNDQkHQu8HdgE3EZRMQaKCjLfH69BBkEwdsJRMJpTgevTvNoRwI1m9gNJ9wFrJP09sIGikkxjKAte96FP3vz0pk2eFPemlQ/Mzngzy5ui2QHhv2RloUllDoH2dhll+70ZWTb/UrZ2zZ/nr1lWDMWP368zKyvIks1QgNNPP73Vzg6UzZs3t/p8eFkO1fKmbpXTYNBoksCqQx3v568oysG3928FzhuPQQVB0B/6GSYlaRnwZYoSedeZ2efa9h8NfBN4I0UR4/eb2W/Svisoln29DPyVma1L/auBdwFPmNnZ/Rjn4Pm3g2DA6If5mSy1rwHvAJYAl0ha0nbYZcCetHb1Goq1rKTjVgBnAcuAr6frAfxL6usbg6uTd6Cslqf3pnlT1H8hsvnpvZ8+i0X2zvm1aT53WjbpfLWpMvOzjsk53mZH2To2f8+yFNz+uXyVrLJqUn6tXg6j8u/Ke0JzeJYP0yrzEA9iuu8+fQ/OA7YkCw1Ja4DlFAWKM8sp1qsC3AR8NVVhXw6sMbP9wMOp2PF5wM/N7GeSFvZjgJnB+FSDYIDpQlObLWnI/ax0l5kH+MnhsrWprWPSqoh9wKya5/aN0NSCoOF0oantMrOl4zmWiSCEWgWdvJ/em+ZNqpxZwnvmPNl8euUrX9nq89fN5lFVFo5MNx7PfuLfSydPajbbvUfUh3R5s/yJJ54AYPfu3a0+7xXOCSX9e8shZzDs9fRJKv3nkt9hWWhXk+njco0dwAK3XbY2NR+zXdJ0YAaFw6DOuX0jzM8gaDgHDx6s9dOBu4DFKTvPURQT/2vbjllLsWYVijWst6bY8bXACklHS1oELAbu7NsDthGamsNrIWWams/V5TU1H1idNQ0/Ce7JmlpVYZasyVRpauOhlfXq8q9ax1a2v8ppkNN0e03NP2uubeodBT6ULGtoXlMb7/c2VejHs5vZAUkfBdZRLOlYbWYbJV0FDJnZWoq1qt9KjoDdFIKPdNyNFE6FA8DlKfQSSd8FLqCYz9sOXGlmPa15DaEWBA2nXwLdzG4Bbmnr+7Rrvwi8t+Lcq4GrS/ov6cvgHCHUgqDBNC0Eqg4h1CroVFXJm1E+NKcs7bTnxBNPBEbmYyubfPdfxH7mSOsnedxVYylbp1a1Nixn3PDr1LxTIb+v/P5g5LrB/I78OVWZTAaNQXv2EGpB0HBCqAVB0CgiSWQwijLz0+NNnmxGVf13zJkjvPd0vBM7Hi7KvKJVntLsCa0Kc8qeTp95w3ujM1WVtZrwPsdCzKkFQdA4QqgFwMgJ7axdVGkZ/kvjtbZDXbdqwnzQTIVMfodldUuhs7Om/TrBMIP2TkKoBUHDCaEWBEFj6GeSyKlCCLUKulkbVhZS1em6VV+0Qcnx1U5+h1XmZSfztP06wTCDpqnV/guSNE3SBkk/SNuLJN0haYukG1KQaxAEk4xBK7zSjVrwMYoqUpnPA9ek1L17KFL5BkEwyRg0oVa3Qvt84I8pAlL/JqXofSvwgXTI9RRpfK8dhzEedsqqInm8yZnXUVXlHcvr2LyXtCw7SBPMqLI/lKo/nrzmzOeW86ZoDjvzBYx9lo9MWT43aMb7HCtNElh1qDun9iXg48AJaXsWsDel7IVDpOdNKYFXApxwwgllhwRBME40TQurQ0ehJimXr7pb0gXd3sDMVgGrAObOnTtl3m5ZQHnVJLXP65UDrn1d0GeffbbVzgHbZcVYoFxTqypscqi+iaDTfetm8YVhDdcHrHsNOL8vH/Duc6fld+TPCU2tILyfo3kz8G5J7wSOAU6kqP03U9L0pK2Na3reIAjGzqBpah0dBWZ2hZnNN7OFFJksbzWzPwVuo0jZC0UK3++P2yiDIBgz4SiozyeANZL+HthAkcp3SlNVzCRP6vuJaT/R783PnG46p5+GkeZnTle9b9++Vt+sWbNa7U45yLIp0c8vYdn1u6FTEZiq3HD+vjNnzgRGpuv2JmNOj+7Tffs8dvlaPsdaWajboNE0gVWHroSamd0O3J7aWykKkgZBMIkJoRYEQaMIoRYA5eanz/VVZX6ecsopwHCdSoBt24aLU2fzKde5hJH1LfO1OqXArqp8Nd5f4LoeTxj2FlflNfMpzXM9T29+eu9mrqPq31vZZ+DNT59vbVDDz2DwvJ+D+0kHwQBQ10lQ55+hpGWSNqfQyE+W7D86hUxuSSGUC92+K1L/ZkkX173mWAhNrQL/IWdNw2sOXmvzq+Bnz54NwKte9apW37333ttqZ6eB195OP/30VnvOnDlA9cr43PbaT5Wm1qkwSq/k61atQytbp+afyztT5s+fD4xcp+Zrpz7yyCPASEeBv27+DKoiEvK7GESNrR+fv6RpwNeAt1Mstr9L0lozu88ddhmwx8zOlLSCIpTy/ZKWUKycOAs4DfiJpFenczpds2sG7xMOggGjT5raecAWM9tqZi8Ba4DlbccspwiZBLgJuDCFVC4H1pjZfjN7GNiSrlfnml0TQi0IGk4XQm22pCH3s9JdZh6wzW2XhUa2jkmL8vdRhFRWnVvnml0T5mcFnRwFvq6nD3OaMWMGAGeccUarb8GCBa325s2bAdi+fXur77e//W2rneNjfYGRMvOzam1YJ6dB2f5OE8lV9ypbM+fN4rIEAD58bN684e9vfkd+fI8++mirvXXrVmDk+j4fEpWdDt5p0ym33SDQZZLIXWa2dDzHMxHEpx4EDadPc6o7gAVuuyw0Mh+zXdJ0YAbwVIdzO12za8L8DIKG06c5tbuAxSk57FEUE/9r245ZSxEyCUUI5a1WXHgtsCJ5RxcBi4E7a16za0JTq6DMzPLeT29++hxf2bzy3s/Xve51rfbjjz8OwFNPPdXqe+ihh1rtHC7k17n5tVfZpKoKPeoUDtTrf+1O2UvK1qR5kzB7d2GkiZ5DxbzH8/7772+1s/fTTwF4Uza3/b2q6rQOGv3Q1MzsgKSPAuuAacBqM9so6SpgyMzWUoRKfkvSFmA3hZAiHXcjcB9wALjczF4GKLtmr2MNoRYEDadfS3rM7Bbglra+T7v2i8B7K869miLJbMdr9koItRqUaWo+H9ozzzzTamdNIUcWAJx11lmtdnYQ3Hff8FIcv2Ytr93yq+F9wHvW1LwWUubU8OP2dPMF77QOrSxiwO/PGqYfv9fO/Pq8fI0HHnig1bdx4/A/7SeffBIYOfnvk45mR4F/b1U56QaJCGgPgqBxDFqYVAi1IGg4oakFo8hfCm9meeeANz/z+jK/ds2bXOeeey4wMsdaDtaG4fVY3jngyY6EqmDtssnxKkdCp9xnZeZnp7qbPkwpj3XRokWtPt/2k/rZEbBhw4ZWn3egZNPfOwd8SFV+7+EcGE0ItSAIGkPMqQVB0DhCqJUg6TfAM8DLwAEzWyrpZOAGYCHwG+B9Zran6hpNwJtbPrV3mffTe+myGQbDnlBvfnr27t0LwMMPP1x63xxO5K/pzbgyj19VGFXZmrZO69D8/mzqeVO5LPPGwoULW33eLPdmdzY7vcczvwsYfkb/3N4UzWPolIduEBk0odbNp/5HZnaOiw37JPBTM1sM/DRtB0EwyTh48GCtn6bQy78yn2bkeuBPeh9OEAT9pJ9JIqcKdefUDPixJAP+lxUFiueYWbYfHgfmVJ7dEKoWufrQnZxFwpufOXMHwNy5cwE455xzSu+xadMmYGTRXr84N5u9PtzIm3xV2T3GQjZPvafVXzN7Ov39/bjys3rz2Juc+VkB7rnnHgB27txZeq9sdnqPZ1VCyGAkTRJYdagr1N5iZjskvRJYL+l+v9PMLAm8UaScTCth5ArwIAgmhhBqJZjZjvT7CUk3U2Ss3CnpVDN7TNKpwBMV564CVgHMnTu3MW/Xf1G80yAHulfVnMwT5b7YiidrH95R4IO8czrrsnvCSK0whw5V1cIs6yvL3ebP98VSstbkNTX/jys7GLym6dee+ZCoHTtGZ5zxWlnW1MqcAzDY6bo7MWhCreM3QNLxkk7IbeAi4F5GphmJCu1BMAnJSSIHyVFQR1ObA9yc/hNOB/7VzH4k6S7gRkmXAY8A7xu/YQZBMFYGTVPrKNSsqMT+hpL+p4ALx2NQU4GqfGY5nKdqHVo2ybwZ5yfXs4PBm5E+9Xc2RavWyfk1Zfv37wdGOg/85HqewPcmpzfp8v6yvGX+ut689mvLsqnsTUtviu7atYt2vMnp30E22zuZ0sFoQqgFQdAoQqgFQdAoQqgFXVNWTSmbfu1kU9Wvc/PmYVkWDm+S5WSJfh2bXyfnTeGy7CJlWTa8+ejX1+Ux+P0+UWZ+xrJ1ejBsXvrU5d4s92ZvNi99GFWZp7NqPVqYouU0bWFtHUKoBUHDaZJnsw4h1PpMmXbktbayIHG/P0/k+wlxnxo8ay/eOeDXqZVpbV7T8ppO1m78l75srD51edm4/f3LxuW1Uh9d4J0W2XHi+7pZXxdUM2iaWnwrgqDhTETsp6STJa2X9GD6fVLFcZemYx6UdKnrf6OkX0vaIukrSv+JJb1X0kZJByXVKrQcQi0IGswEBrR3zNqT0pVdCbyJIirpSif8rgU+QlETdDGwLPXfC/xH4Gd1BxLm5zhR5jzw/b7PrznLJl1ZrU8oX1vmHQ1+Ij+3q1Jwl6Xr9uZrTlnuzceysXqT1R+bx+1NTv9cZaamf9ayfG9hcnbPBJmfy4ELUvt64HbgE23HXAysN7PdAJLWA8sk3Q6caGa/SP3fpMj680Mz25T6ag8khFoQNJwJchTUydozD9jmtrenvnmp3d4/JkKoBUGD6dK0nC1pyG2vSgkpAJD0E2BuyXmfartnZdaeiSCE2gRQZopWhVll864qh1k2z6oK9ZaZb94kLGt787Rsf9X5+RnK1rb5cVc9S1V/2XMFY6cLobbLZbYuu87bqvZJqpO1ZwfDJirAfAozdUdq+/7RaVtqEt+aIGg4E+QoqJO1Zx1wkaSTkoPgImBdMluflnR+8np+qOL8WoSmdpgom7D3/V4jKtPEOhVQ8dfq5gtbpgGWrW3z163SGjudH46AiWGCHAWfoyRrT1qG8V/M7MNmtlvSZ4G70jlXZacB8JfAvwDHAj9MP0j6D8A/AqcA/0fSPWZ28aEGEkItCBrORAi1qqw9ZjYEfNhtrwZWVxx3dkn/zcDN3YwlhFoQNJicJHKQCKE2CfD/SXO7ypGQJ/XrmJ9l1+/0BR9LOFLVWDqdF6bmxDBoYVIh1IKg4YRQC4KgUYRQK0HSTOA6iok8A/4C2AzcACwEfgO8z8z2VFwi6JIyk9TTTdhInXuUXbdftTTDzDx8DGI+tbrfti8DPzKz11LUK9hEjQDWIAgOP1GhvQ1JM4A/BP4MwMxeAl6SVCeANRgn+vkl9JpUN56y0MCmBoPm/azzrVwEPAn8s6QNkq5L9T/rBLAGQXCYGTRNrY5Qmw78PnCtmZ0LPEebqWnFGyl9K5JWShqSNOQzqAZBMP5MYD61SUMdobYd2G5md6TtmyiE3M4UuMohAlgxs1VmttTMlvpal8Hk5Igjjqj9E0wNQqi1YWaPA9skvSZ1XQjcR70A1iAIDjODJtTqrlP7r8B3JB0FbAX+nEIgjgpgDYJgcjFojoJaQs3M7gHK8iyNCmANgmDy0DQtrA4RURAEDSeEWhAEjSKEWhAEjSKEWhAEjSKEWhAEjSGSRAZB0DhCUwuCoFGEUAuCoFEMmlCLAL4gaDATFdAu6WRJ6yU9mH6fVHHcpemYByVd6vrfKOnXkrZI+kqq/4mkf5B0v6RfSbo5Jaw9JCHUgqDhTFDsZ8eksZJOBq4E3gScB1zphN+1wEeAxelnWepfD5xtZq8HHgCu6DSQEGpB0HAOHjxY66dHllMkiyX9/pOSYy4G1pvZ7pT6fz2wLGX5OdHMfpHSmH0zn29mPzazXNn7F8By7GkJAAAHy0lEQVT8TgOJObUgaDgTNKdWJ2nsPGCb296e+ualdnt/O39BURflkIRQC4IG06VpOVvSkNteZWar8oaknwBzS877VNs9TVJfJamkTwEHgO90OjaEWhA0nC6E2i4zK8vGk6/ztqp9knZKOtXMHjtE0tgdDNc1gcKUvD31z2/r3+Gu/WfAu4ALrcbDxJxaEDScCXIU1Ekauw64SNJJyUFwEbAuma1PSzo/eT0/lM+XtAz4OPBuM6tVDyCEWhA0nAlyFHwOeLukB4G3pW0kLZV0HYCZ7QY+C9yVfq5KfQB/SVFbeAvwEPDD1P9V4ARgvaR7JP1Tp4GE+RkEDWaikkSa2VOUJI01syHgw257NbC64rizS/rP7HYsIdSCoOFEREEbkl6T1L7887Skv667gjgIgsPLoBVeqVNNarOZnWNm5wBvBJ4HbqbGCuIgCA4/IdQOzYXAQ2b2CPVWEAdBcJgZNKHW7ZzaCuC7qV1nBXEQBIeRQUwSWVtTSzU/3w387/Z9aUFcqaiXtFLSkKSh55+vtcwkCII+MmiaWjfm5zuAX5rZzrS9M60c5hAriDGzVWa21MyWHnfccb2NNgiCrgmhVs0lDJueUG8FcRAEh5kQaiVIOh54O/A91126gjgIgslDXYHWJKFWy1FgZs8Bs9r6SlcQB0EwuWiSwKpDRBQEQcMZNO9nCLUgaDihqQVB0BiaNl9WhxBqQdBwQqgFQdAoQqgFQdAowlEQBEFjiDm1IAgaRwi1IAgaRQi1IAgaRQi1IAgaxaAJtSiRFwQNJieJHO8SeXVrlki6NB3zoKRLXf8bJf1a0hZJX0n1P5H0WUm/SvVRfizptE5jCaEWBA1ngrJ0dKxZIulk4ErgTcB5wJVO+F0LfARYnH6Wpf5/MLPXpxopPwA+3WkgIdSCoOFMkFCrU7PkYmC9me02sz3AemBZSjJ7opn9ImXR/mY+38yeducfT0WGbU/MqQVBw+lCYM2WNOS2V5nZqprn1qlZMg/Y5ra3p755qd3eD4Ckq4EPAfuAP+o0kBBqQdBgutTCdpnZ0qqdkn4CzC3Z9am2e5qkvnknzOxTwKckXQF8lMKErSSEWhA0nH55P83sbVX7JO2UdKqZPXaImiU7gAvc9nzg9tQ/v61/R8n53wFuoYNQizm1IGg4E+H9pF7NknXARZJOSg6Ci4B1yWx9WtL5yev5oXy+pMXu/OXA/Z0GUrdGwX+TtFHSvZK+K+kYSYsk3ZFcsDekEnpBEEwyJshRUFqzRNJSSdelcewGPgvclX6uSn0AfwlcB2wBHgJ+mK+b5M6vKITgxzoNpKP5KWke8FfAEjN7QdKNFEWN3wlcY2ZrJP0TcBmFWzYIgknCRAW0V9UsMbMh4MNuezWwuuK4s0v6/1O3Y6lrfk4HjpU0HTgOeAx4K3BT2l/lwg2C4DAzaNWkOgo1M9sB/E/gUQphtg+4G9hrZgfSYSNcsJ6o0B4Eh5cQam2kCb3lwCLgNIoFcMsOeZLDokJ7EBxWJshRMGmos6TjbcDDZvYkgKTvAW8GZkqanrS1KhfsCHbu3LnrC1/4wnPArh7GPFmZTTzXVKGJzwTFcx3f1rcu9dehEe+kjlB7FDhf0nHACxSTgUPAbcB7gDVUu3BHYGanSBo61AK/qUo819Shic8Ereda6PvMrLZV1RTqzKndQeEQ+CXw63TOKuATwN9I2kJRvf0b4zjOIAiCWtSKKDCzKxm9incrRaR9EATBpOFwRBTUDZCdasRzTR2a+EzQ3OfqCjXJlRsEQRCxn0EQNIoJFWqSlknanOJFR2XGnApIWiDpNkn3pXjYj6X+WumMJzuSpknaIOkHaXvKx/hKminpJkn3S9ok6Q+m+ucV8djVTJhQkzQN+BrwDmAJcImkJRN1/z5yAPhbM1sCnA9cnp6jYzrjKcLHgE1u+/MUMb5nAnsoYnynGl8GfmRmrwXeQPF8U/bzcvHYS83sbGAaRTx2Ez6rnplITe08YIuZbTWzlyjWty2fwPv3BTN7zMx+mdrPUPyBzKNeOuNJjaT5wB9TZEsgpYGZ0jG+kmYAf0hacmRmL5nZXqb+5xXx2BVMpFCrSuU7ZZG0EDgXuIN66YwnO18CPg7kmJlZ1IzxncQsAp4E/jmZ1ddJOp4p/Hn1Go/ddMJRMEYkvQL4N+Cv24pDkIpHTCm3sqR3AU+Y2d2Heyx9Zjrw+8C1ZnYu8BxtpuZU+7x6jcduOhMp1HYAC9x2rXjRyYikIykE2nfM7Hupe2dKY8wh0hlPZt4MvFvSbyimBt5KMRc1M5k4MDU/s+3A9hQZA4V59vtM7c+rFY9tZr8DRsRjp2Om4mfVFyZSqN0FLE4emqMoJjbXTuD9+0KaZ/oGsMnMvuh21UlnPGkxsyvMbH6KHVwB3Gpmf8pwjC9Mzed6HNgm6TWp60LgPqb259WKx07fx/xMU/qz6hcTuvhW0jsp5m2mAavN7OoJu3mfkPQW4P9SxMHmuae/o5hXuxE4HXgEeJ9LVTylkHQB8N/N7F2SzqDQ3E4GNgAfNLP9h3N83SLpHArnx1EU4X1/TvEPfcp+XpI+A7yfwhu/gSK77Dym+GfVDyKiIAiCRhGOgiAIGkUItSAIGkUItSAIGkUItSAIGkUItSAIGkUItSAIGkUItSAIGkUItSAIGsX/B6PSeL09jZFHAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAATUAAAD8CAYAAAAWjzPMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztnXuQXNV95z9fS7wfEpJAAklYIsgmgrIhVmFSTqWIsUF2XFZ21w/heE0SbNWW8cbZZMs2cZXZ4FBlV9YhTmKTVWEl+BELlpiyyoutyAYqu1U2METERghZQhgkRUgIPXgjC/32j3tOz2967lXfVveMZm7/PlVTc/rcR597eubX53d+L5kZQRAETeF1x3oAQRAE/SSEWhAEjSKEWhAEjSKEWhAEjSKEWhAEjSKEWhAEjSKEWhAEtZC0VNImSVskfabk+AmSbk/H75e0wB27PvVvknRV23VTJK2X9L1+jLMnodbpIYMgaAaSpgBfAd4FLAaulrS47bRrgX1mdj5wM/DFdO1iYDlwIbAU+Gq6X+aTwMZ+jfWohVrNhwyCoBlcCmwxs61mdhBYDSxrO2cZcFtq3wlcIUmpf7WZvWpmTwBb0v2QNA/4beDWfg10ag/Xth4SQFJ+yEerLjjppJNs2rRpPbxlEARH4sCBA7z88svKr5cuXWp79uypde1DDz201syWVhyeC2xzr7cDb606x8wOSToAzEz9P2m7dm5q/xXwKeC0WoOsQS9Crc5DjmDatGl8+MMf7uEtgyA4Et/85jdHvN6zZw9DQ0O1rpV0gSR/8kozW9nH4bW/33uA3Wb2kKTL+3XfXoRaLSStAFYAnHZa34RxEAQ16SK+e4+ZLak4tgOY717PS31l52yXNBWYBjx7hGvfC7xX0ruBE4HTJX3TzHpa+fRiKKjzkJjZSjNbYmZLTj755B7eLgiCo+Hw4cO1fjrwILBI0kJJx1Ns/K9pO2cNcE1qvw+4xwqJugZYnqyjC4FFwANmdr2ZzTOzBel+9/Qq0KC3lVrrISmE2XLgQ70OKAiC/mFm3azUjnSfQ5I+AawFpgCrzGyDpBuBITNbA3wN+IakLcBeCplAOu8Oiv32Q8B1ZvZaz4Oq4KiFWtVD9m1kQRD0hX6lFzOzu4G72/o+59qvAO+vuPYm4KYj3Ps+4L5+jLOnPbWyhwyCYGIxaDkTx9xQEATBsSWEWhAEjSKEWhAEjcHM6lg2G0UItSBoOLFSC4KgUYRQC4KgUYRQC4KgMfTL+XYyEUItCBpOGAqCIGgUsVILgqAxhPoZBEHjCKEWBEGjCKEWBEGjCKEWBEFjiDCpIAgaR6zUgiBoFCHUgiBoFCHUgiBoFIMm1DpWk5K0StJuSY+4vhmS1knanH6fMbbDDILgaMiGgj5Uk5o01CmR9w9Ae9XmzwA/MrNFwI/S6yAIJiA5qqDTT1PoKNTM7F8oyl15lgG3pfZtwO/0eVxBEPSJfgk1SUslbZK0RdKohUyq63l7On6/pAXu2PWpf5Okq1LfiZIekPRvkjZI+rN+PO/RFjOebWY7U/tpYHY/BhMEQf/ph1CTNAX4CvAuYDFwtaTFbaddC+wzs/OBm4EvpmsXU9QAvZBC6/tqut+rwNvN7M3AxcBSSZf1+ry9VGgHIFVgrpwRSSskDUkaeumll3p9uyAIuqCuQKuxUrsU2GJmW83sILCaQmPzeA3uTuAKSUr9q83sVTN7AtgCXGoFL6Tzj0s/PevBRyvUdkk6GyD93l11opmtNLMlZrbk5JNPPsq3C4LgaOlCqM3KC5D0s8LdZi6wzb3envooO8fMDgEHgJlHulbSFEkPU8iQdWZ2f6/Pe7QuHWuAa4AvpN/f7XUgQRCMDV1YNveY2ZKxHEs7ZvYacLGk6cBdki4ys0c6XXck6rh0fBv4MfBGSdslXUshzN4paTPwjvQ6CIIJSJ/Uzx3AfPd6XuorPUfSVGAa8Gyda81sP3Avoz0tuqbjSs3Mrq44dEWvbx4EwdjSR3eNB4FFkhZSCKTlwIfazska3I+B9wH3mJlJWgP8o6S/BM4BFgEPSDoT+KWZ7Zd0EvBOknGhFyKiIAgaTj+EmpkdkvQJYC0wBVhlZhsk3QgMmdka4GvANyRtoXADW56u3SDpDuBR4BBwnZm9lvbjb0uW0NcBd5jZ93odawi1IGg4/XKsNbO7gbvb+j7n2q8A76+49ibgpra+nwKX9GVwjhBqQdBwmhQtUIcQakHQYCJJZBAEjSNWakEQNIoQakHjKCJV+sOg/YM0gUH7zEKoBUHDCaEWBEFjCENBMGnwKmWndtXxzOteNxwtV/YP4L/py9qdjgfHlkH7HEKoBUHDCaEWBEGjCKEWTDi8yphVRa8yTp06/DFOmTJlVL8/7q/LbX9//w+QVVGvkh46dGhU+7XXXis9Xnb9oP2DHWuaVn+gDiHUgqDhhFALjil51eRXVH71dfzxx4/4DXDiiSeOOu77jzvuuNJ75XbVSi2vwPxK7Je//GWr/corrwBw8ODBUX2+3x/394oV3PgQ1s8gCBrFoH1hhFALggYTe2rBuFHlO5Y39b3KeNJJJ7XauXjNKaec0urzba+KdjIU1KWTocCrnC+++OKotq8i9vLLL7faXpXN9wo/t/4zaPNYp0bBfEn3Sno0FRz9ZOqfIWmdpM3p9xljP9wgCLqlX8WMJwt1vrYPAX9iZouBy4DrUnHSzwA/MrNFwI/S6yAIJhiDJtTqFF7ZCexM7eclbaSo2bcMuDyddhtwH/DpMRllgyizbnpVM6uPp512WquvrO1VUq9e+j/OrN49//zzrb5XX3211S5T+cpU4RNOOKHVV9b2Ku+pp57aamdV079/VTursF4lDeto70TsZwckLaDIKX4/MDsJPICngdkV16wAVsDIf84gCMaHQftCqL1rLOlU4J+APzKz5/wxK2atdOaiQnsQHFtC/SxB0nEUAu1bZvad1L1L0tlmtjOVuto9VoOc7HiVLju8eidZb72cNm3aiN/tx7NK6NXIPXv2lLafeeYZAPbu3dvqe+654e+jrB56i6ZXZbOKe/rpp7f6ZsyY0WqfeeaZAMyaNWvU+P24vcrq1Wbff+DAAWCk9bTMabdJ/3zjRb/mTNJS4MsUJfJuNbMvtB0/Afg68BaKIsYfNLNfpGPXA9cCrwF/aGZrJc1P58+mWBStNLMv9zrOOtZPUdTz22hmf+kO5cKlpN/f7XUwQRD0n36s1FJtzq8A7wIWA1cng6HnWmCfmZ0P3EwqTJzOWw5cSFGB/avpflVGyJ6os1J7G/CfgZ9Jejj1/SnwBeAOSdcCTwIf6HUwTaJsdQbDqxO/oX7GGcPeMHml4zff/apsx44dADz55JOtvq1bt7bavn/nzmLLc9++fa0+vxKq8YcMjFwp+rGeffbZALz+9a9v9Z133nmtdu7PKzoY+dx+VZjb3oDywgsvtNp5DnyYVazaOtNHQ8GlwBYz2wogaTWFsfBRd84y4H+k9p3A36ZF0TJgtZm9CjyRih1famY/ptwI6e/ZNXWsn/8PqEpyf0Uvbx4EwdjThfCfJWnIvV5pZitTey6wzR3bDry17frWOami+wFgZur/Sdu1c/2FbUbInoiIgiBoOF0ItT1mtmQsx1LGkYyQR0MItT6TVbYylROG1S+/4e434vN1eZMfRqqXGzZsAGDjxo2tvm3bhr9AfUhSr+R/Bq8G+nZ+30ceeaTVN3/+/Fb7V3/1VwG48MILW31ePZ0+fXqrndXuTpWvvCoeqmg9+jQ3O4D57vW81Fd2znZJU4FpFAaDymsrjJA90X0gYBAEk4a6RoIagu9BYJGkhZKOp9j4X9N2jjcevg+4J7l7rQGWSzpB0kJgEfDAEYyQPRErtSBoOP1YqaU9sk8AaylcOlaZ2QZJNwJDZraGQkB9IxkC9lIIPtJ5d1AYAA4B15nZa5J+gxIjpJnd3ctYQ6j1mWy986FP3uk4q1neCuitU//+7/8ODKuZAOvXr2+1N2/eDAz7dbXj75tV3CpVN/uMVVWTyn5s3rfN+7zltldJN23a1Go//fTTAGzfvr3Vd8kll7TaXi2dM2fOqPGXJZSsCp3y5wYj6VeYVBI2d7f1fc61XwHeX3HtTcBNbX1HMkIeNSHUgqDhDNp+Ywi1PlCWett7zvuY17xq84HbeXUG8PDDxSrcr86eeuqpVjuvSM4555xW37x581ptv1F/1llnASO9/P2qscw3rCx3mjc++BXi7t1FEIk3VPhV2bPPPgvAo48Oux1VGR0uvvjiUc/lx1pW5KXMUDBowdudaFoIVB1CqAVBwwmhFgRBowihFtSirBYnDPukefXTt7PKtGvXrlZfVjlhWO3MIU4wUn3N6uWv/MqvtPq8yuaNAjnUyvvMVVWOKiOf6zfvZ86cOep9zz333FafV6Uff/xxYKR66p+rE7NnD2ezynPo/dTKcq9FOvDRDNo8hFALggYTSSKDIGgcsVILalGW9hrKU1z7P6psEXzsscdafT7kaf/+/cDIHGU+tCi3vWrmrYRllkyfL61MPfPf5P76/IydVG0/Fq8q5zAon9nDh3zlZ4XhOSirhuXvVZW9JKuikQJ8NIM2DyHUgqDhhFALKqlasfjogZzR1h/3OczypvkTTzzR6vMe+3mjf+HCha0+384rOB8kX1ZsBTp74dc1FFQ9d277Pm8UyfnWyuYHyufA93n/unwPf71v5+OdVqWDyKA9ewi1IGgwYSgIgqBxxEotqKRKDfN+YHlz26uBvhhKDiPyKbb9Rn/2Q/Mqp/cNyyqX//atCh0qUz+P5lu7TOX0bf/8vp3H6sfv8WPJarmfFx9ylQPxfWpwb0jI71sV8jVo/9ieQXv2OoVXTpT0gKR/k7RB0p+l/oWS7pe0RdLtKcdSEAQTjEErkVcnSeSrwNvN7M3AxcBSSZdRVIq5OVWO2UdRSSYIgglEH5NEThrqFF4xIKdTOC79GPB24EOp/zaKKjK39H+IE4cq3zTfzn8cPrOFT82d1SuvpnmVKqufPtW1v39WL72Vr0r9bB9Te7su/poy9bVK/czjLvM3g5EZRbK67vO1eVU0z6GvbOUtwGXvFTnWCpoksOpQt5jxFOAh4HyK2n+PA/vNLP9njaoO465dAayAkY6ZQRCMD2H9LMHMXgMuljQduAu4oO4bpBJbKwHmzJkzKb8yyvy1qtp5BfX888+3+rwfWl6ReAHvPfLzSsb7dpUZBcoMAtB5Vdbrt7a/Pj93VQGU3ParN/9cftWW58AbWHw7z6Gf17JA/U6f0aCtWpqmWtahK+unme2XdC/w68B0SVPTaq2sskwQBBOAQRNqdayfZ6YVGpJOAt4JbATupagYA0UFme+O1SCDIDh6wlAwmrOB29K+2uuAO8zse5IeBVZL+nNgPUUlmUZTpdr4P4gcZO1Do1555ZVWO6tfXv30Ad9587vKtyyrt1Vq5liEBlXdp0zlKzMq+L6ygHgYngM/b17VzHPoj/uQrHzfqs9okGmSwKpDHevnTynKwbf3bwUuHYtBBUHQH/oZJiVpKfBlihJ5t5rZF9qOnwB8HXgLRRHjD5rZL9Kx6yncvl4D/tDM1qb+VcB7gN1mdlE/xhnFjIOg4fRD/Uya2leAdwGLgaslLW477VpgX/JdvZnCl5V03nLgQmAp8NV0P4B/SH19I8Kk+oD3Gcttr3L6b8ocEuUzUHg1KtMp9GkiZKDI71tmEfVtPz9VqcXzHPh5KXtuP68HDx5stb1/WjCSPv19XApsSRoaklYDyygKFGeWUfirAtwJ/G2qwr4MWG1mrwJPpGLHlwI/NrN/kbSgHwPMxEotCBpOFyu1WZKG3M8Kd5u5wDb3usw3tXVO8oo4AMyseW3fiK+3IGg4XazU9pjZkrEcy3gQQu0oqcqSkVUif9wnM8wWP5+Zw6tknTJrlKl8nmOhinaqUFW1UV2W6cTPi0/XXZZ63KufeQyD5j3fiT66a+wA5rvXZb6p+ZztkqYC0ygMBnWu7RuhfgZBwzl8+HCtnw48CCxK2XmOp9j4X9N2zhoKn1UofFjvSbHja4Dlkk6QtBBYBDzQtwdsI1ZqXVAV2O3786rNr168P1ZekfjVmydf38n3bCL5HnUzlk45zvy8lNVL9fPqV8id/PcGmX7Mg5kdkvQJYC2FS8cqM9sg6UZgyMzWUPiqfiMZAvZSCD7SeXdQGBUOAdel0EskfRu4nGI/bztwg5n15PMaQi0IGk4fnbDvBu5u6/uca78CvL/i2puAm0r6r+7L4Bwh1IKgwTQtBKoOIdS6oMo4UHaO3/z3tSo7pZ3OdKNGTaQ/2jKjQSc/No+fFz9veb7LjCow0heu7PggM5H+PsaDEGpB0HBCqAVB0CgGbcUaQq0GZd90ZeqOP9cnQ/Qqk1evMp3Uz8lEN6pypywa3mqc58jPn79X1edRd1xNJfbUgiBoHCHUglF0s3zPKwm/UitbXXQqltINZas/Tz/Vj7F+r6p6qmURA3UMN0e6/6AQQi0IgkYRQi0IgsbQzySRk4UQahV02ryvCvfJeb3Kgtihc761TDdpqcfzj7bX9+rkf+dVTj+HeV59kLufwzyvfq68qlmmvg6KKjpoK7Xan6qkKZLWS/peer1Q0v2Stki6PQW5BkEwwRi0wivdfFV9kqKKVOaLwM0pde8+ilS+QRBMMAZNqNWt0D4P+G2KgNQ/Til63w58KJ1yG0Ua31vGYIwThiprm+/P6pO3fvpwn1yg1xfq9XnB8h9XlWo0GSskVamcft7yc3n106vteT79NS+//HKrndXPqnnz9x00miSw6lB3T+2vgE8Bua7bTGB/StkLR0jPm1ICr4CRZeGCIBh7mrYKq0NHoSYpl696SNLl3b6Bma0EVgLMmTNnUs9uWd609nbZisCvOPLmt1+dla04qqIQylYiE2n11inzbdVqN6/EvHHglFNOabXzM7700kutvrKCN52iNwaRsH6O5m3AeyW9GzgROJ2i9t90SVPTam1M0/MGQXD0DNpKreNXmZldb2bzzGwBRSbLe8zsd4F7KVL2QpHC97tjNsogCI6aMBTU59PAakl/DqynSOU7MFQFU2eVx6tWviZlVq+8v9ULL7zQame1tFPYTx2Vc6zV0rJ/hDKjQJXK6c8tK7zi1c+yefHGlk6GgkGlaQKrDl0JNTO7D7gvtbdSFCQNgmACE0ItCIJGEUItALpLt+2PZ7WyrGYlDFdIOv3001t9zz33XKudQ386qWlVFtFONTiPRiXtJkdamaWzk08fDPvy+Xnx1aTK5rVsjjuNdRDV00Gzfg7eJxwEA0RdI0Gd1ZykpZI2pdDIz5QcPyGFTG5JIZQL3LHrU/8mSVfVvefRECu1LqgKMi8rAPLiiy+2+nzgdfbHOuOMM1p9/tzsh+X9sToZDfwfZCc/tl5VkU5+aGW+fFXj9z5peT78vPhx5zn0c+VXatkY000igEGhH+qnpCnAV4B3UjjbPyhpjZk96k67FthnZudLWk4RSvlBSYspPCcuBM4BfijpDemaTvfsmlipBUHD6dNK7VJgi5ltNbODwGpgWds5yyhCJgHuBK5IIZXLgNVm9qqZPQFsSferc8+uCaEWBA2nC6E2S9KQ+1nhbjMX2OZel4VGts5JTvkHKEIqq66tc8+uCfWzCzrl6oJh9fP5559v9R04cKDVPvvsswGYPn16q2/27NmtdlY7fRhVWe61qg33Tqpor5QZBar80HLbH/fB/V7VzHPg44P379/fauc59PPq/QOzWl/1uQwqXSaJ3GNmS8ZyPONBCLUgaDh9cunYAcx3r8tCI/M52yVNBaYBz3a4ttM9uybUzyBoOH3aU3sQWJSSwx5PsfG/pu2cNRQhk1CEUN5jxY3XAMuTdXQhsAh4oOY9uyZWal1Q5RtWlpnDW+n27NnTau/duxeAmTNntvrmzJnTaufQH2/Ze/bZZ1vtrJZWZb7oNK5u1LOyMCdPmXpZ5l/nrZxe7Z43b16rnefAX5/nCobn0M+rDz/r5L83yPRjpWZmhyR9AlgLTAFWmdkGSTcCQ2a2hiJU8huStgB7KYQU6bw7gEeBQ8B1ZvYaQNk9ex1rCLUgaDj9iigws7uBu9v6PufarwDvr7j2Jooksx3v2Ssh1GqQv/GrNufLqrH7zX2/0tqxo9gyOPPMM1t98+cPbyssWLBg1Ht58orFr+T86qZspVZl4KhLJz+0qrHmPHJ+Vbpw4cJWOz+rP3fbtmFjWJ4rGJ5Dbxwoy4xbleF2UFdtEdAeBEHjGLQwqRBqQdBwYqUW1KJqQz5vivsU3T5fWlavvD/WtGnTWu3sx1ZVN3Tr1q0A7Nq1q9XnN8875Rjr1VBQpor6e/ocaNn37Lzzzmv1efXz1FNPbbV37twJwOOPP97q86ponsOy1OhQbigICkKoBUHQGGJPLQiCxhFCrQRJvwCeB14DDpnZEkkzgNuBBcAvgA+Y2b6xGebEoCrHWlmNz6q8X9l6+fOf/7zV5/OG5euzGgojVdXs55XVUBhpJdy3b/gjyCFXXiXtFT8HOfW2D3eaO3c4dC+rneeee26rz6uPWeUE2LixqJPt58X792X10odZ+Xb+DHq19DaRQRNq3Xzqv2VmF7vYsM8APzKzRcCP0usgCCYYhw8frvXTFHr5KvNpRm4Dfqf34QRB0E/6mSRyslB3T82Af5ZkwP+yokDxbDPL+sPTwOzKqxtOmfXTV0XyGTdy5ont27e3+jr9Qb3+9a9vtd/0pjcBI9U8byX09929e/eI94SRVtmslvr39+pbVum8euzDnM466yxgZLiTdyTOTrfeOfjJJ59stdevX99qP/zww8BIVdqTLaV+Xsusn8FomiSw6lBXqP2Gme2QdBawTtJj/qCZWRJ4o0g5mVbAyL2hIAjGhxBqJZjZjvR7t6S7KDJW7pJ0tpntlHQ2sLvi2pXASoA5c+Y0ZnarjAY5yNpvYnt/rLw68sVWnnrqqVY7GxW8b5tvv+ENRRZkvzryqzZ/3xwQ7gPDOxV56VQMZcaMGaPa/rhf6WWjhd/896uzDRuGY5e9313G3zfPoZ9XH9CeCePAaAZNqHX8C5B0iqTTchu4EniEkWlGokJ7EExAcpLIQTIU1FmpzQbuSt/AU4F/NLMfSHoQuEPStcCTwAfGbphBEBwtg7ZS6yjUrKjE/uaS/meBK8ZiUJONMlXUb2L70KGyWpheJcyb/j5t9TPPPNNq5430Cy64oNXn/cBmzZrVap9zzjnAyD9qn+WiLMtGWfiXV/O8epkNIN43zqvSjz1WbL16NdOHQXkDRn4vr3L68LE8h35ewyetHiHUgiBoFCHUgiBoFCHUgr7g1SGvMnlLaNm5WSXzIUJePc1+aFm1g5FZMLxPWw618mFM3s8rhyz59/f/ADn7hy+s7FXNHObkfc98+Fbu9+qz99nzIVPZ/837wfm5ynMYamZ3NM2xtg4h1IKg4TTJslmHEGp9plPq77w6qdrkzpvyfhPd50vLqx/v17Vp06ZW2/uRZY9/3+c34nOkgB+fN2Dk6IMy3zcYjljwfb7GafaD84YGbzTxq7I8Ln+8U8RArNrqMWgrtfirCIKGMx6xn5JmSFonaXP6fUbFedekczZLusb1v0XSzyRtkfTXSt/6kt4vaYOkw5JqFVoOoRYEDWYcA9o7Zu1J6cpuAN5KEZV0gxN+twAfo6gJughYmvofAf4j8C91BxLq5xhRFUaV1Si/Se5V0ayq+eNe/cv+az4w3ecl82rp5s2bR93Lq3Q5YL1qrDmky2/u+9xw2eetLEwMhlVJH/Pr1V+vaubwJz++MvUyVM7uGSf1cxlweWrfBtwHfLrtnKuAdWa2F0DSOmCppPuA083sJ6n/6xRZf75vZhtTX+2BhFALgoYzToaCOll75gLb3OvtqW9uarf3HxUh1IKgwXSpWs6SNORer0wJKQCQ9ENgTsl1n217z8qsPeNBCLVxoEy9831ePczqm1fDyjJ+eIuo9yPzRZSzqujVx6P51q6qRpXH7cfn/eCyeun7fG62MutmlZoRaufR04VQ2+MyW5fd5x1VxyTVydqzg2EVFWAehZq6I7V9f3livRrEX0oQNJxxMhTUydqzFrhS0hnJQHAlsDaprc9JuixZPT9ScX0tYqU2zpStOMo22v15vrBLXul4b3u/ee9XanmF5ldqPqA9v29V5ts8Br/571dXue1Xan7Vmdt+/N7fLAwB48M4GQq+QEnWnuSG8V/M7KNmtlfS54EH0zU3ZqMB8HHgH4CTgO+nHyT9B+BvgDOB/yPpYTO76kgDCaEWBA1nPIRaVdYeMxsCPuperwJWVZx3UUn/XcBd3YwlhFoQNJicJHKQCKE2AeikcpWFUXk10G+++zCn7Gfm/6j9cd8uI6uKVSpjViv9ca+qxub/xGDQwqRCqAVBwwmhFgRBowihVoKk6cCtFBt5BvwBsAm4HVgA/AL4gJntq7hF0CVVqllWJb2a59v+D7jMulnVLiOrjV59LLOOdgphCTXz2DGI+dTq/rV9GfiBmV1AUa9gIzUCWIMgOPZEhfY2JE0DfhP4PQAzOwgclFQngDXoM51WPVV53MaSWIlNbAbN+lnnr3Eh8Azw95LWS7o11f+sE8AaBMExZtBWanWE2lTg14BbzOwS4EXaVE0rZqR0ViStkDQkacjHKAZBMPaMYz61CUMdQ8F2YLuZ3Z9e30kh1OoEsJKi/FcCzJkzpzkzN0EJVTBop0kCqw4d/wPM7Glgm6Q3pq4rgEepF8AaBMExJlZq5fxX4FuSjge2Ar9PIRBHBbAGQTCxGDRDQS2hZmYPA2V5lkYFsAZBMHFo2iqsDhFREAQNJ4RaEASNIoRaEASNIoRaEASNIoRaEASNIZJEBkHQOGKlFgRBowihFgRBoxg0oRaBgkHQYMYroF3SDEnrJG1Ov8+oOO+adM5mSde4/rdI+pmkLZL+OtX/RNJfSHpM0k8l3ZUS1h6REGpB0HDGKfazY9JYSTOAG4C3ApcCNzjhdwvwMWBR+lma+tcBF5nZm4CfA9d3GkgItSBoOIcPH6710yPLKJLFkn7/Tsk5VwHrzGxvSv2/DliasvycbmY/SWnMvp6vN7N/NrNcgfsnwLxOA4k9tSA5djXyAAAHvklEQVRoOOO0p1YnaexcYJt7vT31zU3t9v52/oCiLsoRCaEWBA2mS9VylqQh93plyocIgKQfAnNKrvts23uapL5KUkmfBQ4B3+p0bgi1IGg4XQi1PWZWlo0n3+cdVcck1Ukau4PhuiZQqJL3pf55bf073L1/D3gPcIXVeJjYUwuChjNOhoI6SWPXAldKOiMZCK4E1ia19TlJlyWr50fy9ZKWAp8C3mtmteoBhFALgoYzToaCLwDvlLQZeEd6jaQlkm4FMLO9wOeBB9PPjakP4OMUtYW3AI8D30/9fwucBqyT9LCkv+s0kFA/g6DBjFeSSDN7lpKksWY2BHzUvV4FrKo476KS/vO7HUsItSBoOBFR0IakN6ZlX/55TtIf1fUgDoLg2DJohVfqVJPaZGYXm9nFwFuAl4C7qOFBHATBsSeE2pG5AnjczJ6kngdxEATHmEETat3uqS0Hvp3adTyIgyA4hgxiksjaK7VU8/O9wP9uP5Yc4kpFvaQVkoYkDb30Ui03kyAI+sigrdS6UT/fBfyrme1Kr3clz2GO4EGMma00syVmtuTkk0/ubbRBEHRNCLVqrmZY9YR6HsRBEBxjQqiVIOkU4J3Ad1x3qQdxEAQTh7oCrUlCrZahwMxeBGa29ZV6EAdBMLFoksCqQ0QUBEHDGTTrZwi1IGg4sVILgqAxNG2/rA4h1IKg4YRQC4KgUYRQC4KgUYShIAiCxhB7akEQNI4QakEQNIoQakEQNIoQakEQNIpBE2pRIi8IGkxOEjnWJfLq1iyRdE06Z7Oka1z/WyT9TNIWSX+d6n8i6fOSfprqo/yzpHM6jSWEWhA0nHHK0tGxZomkGcANwFuBS4EbnPC7BfgYsCj9LE39f2Fmb0o1Ur4HfK7TQEKoBUHDGSehVqdmyVXAOjPba2b7gHXA0pRk9nQz+0nKov31fL2ZPeeuP4WKDNue2FMLgobThcCaJWnIvV5pZitrXlunZslcYJt7vT31zU3t9n4AJN0EfAQ4APxWp4GEUAuCBtPlKmyPmS2pOijph8CckkOfbXtPk9Q364SZfRb4rKTrgU9QqLCVhFALgobTL+unmb2j6pikXZLONrOdR6hZsgO43L2eB9yX+ue19e8ouf5bwN10EGqxpxYEDWc8rJ/Uq1myFrhS0hnJQHAlsDaprc9JuixZPT+Sr5e0yF2/DHis00Dq1ij4b5I2SHpE0rclnShpoaT7kwn29lRCLwiCCcY4GQpKa5ZIWiLp1jSOvcDngQfTz42pD+DjwK3AFuBx4Pv5vknu/JRCCH6y00A6qp+S5gJ/CCw2s5cl3UFR1PjdwM1mtlrS3wHXUphlgyCYIIxXQHtVzRIzGwI+6l6vAlZVnHdRSf9/6nYsddXPqcBJkqYCJwM7gbcDd6bjVSbcIAiOMYNWTaqjUDOzHcD/BJ6iEGYHgIeA/WZ2KJ02wgTriQrtQXBsCaHWRtrQWwYsBM6hcIBbesSLHBYV2oPgmDJOhoIJQx2XjncAT5jZMwCSvgO8DZguaWparVWZYEewa9euPV/60pdeBPb0MOaJyiziuSYLTXwmKJ7rlLa+tam/Do2YkzpC7SngMkknAy9TbAYOAfcC7wNWU23CHYGZnSlp6EgOfpOVeK7JQxOfCVrPtcD3mVltraop1NlTu5/CIPCvwM/SNSuBTwN/LGkLRfX2r43hOIMgCGpRK6LAzG5gtBfvVopI+yAIggnDsYgoqBsgO9mI55o8NPGZoLnP1RVqkik3CIIgYj+DIGgU4yrUJC2VtCnFi47KjDkZkDRf0r2SHk3xsJ9M/bXSGU90JE2RtF7S99LrSR/jK2m6pDslPSZpo6Rfn+yfV8RjVzNuQk3SFOArwLuAxcDVkhaP1/v3kUPAn5jZYuAy4Lr0HB3TGU8SPglsdK+/SBHjez6wjyLGd7LxZeAHZnYB8GaK55u0n5eLx15iZhcBUyjisZvwWfXMeK7ULgW2mNlWMztI4d+2bBzfvy+Y2U4z+9fUfp7iH2Qu9dIZT2gkzQN+myJbAikNzKSO8ZU0DfhNksuRmR00s/1M/s8r4rErGE+hVpXKd9IiaQFwCXA/9dIZT3T+CvgUkGNmZlIzxncCsxB4Bvj7pFbfKukUJvHn1Ws8dtMJQ8FRIulU4J+AP2orDkEqHjGpzMqS3gPsNrOHjvVY+sxU4NeAW8zsEuBF2lTNyfZ59RqP3XTGU6jtAOa717XiRSciko6jEGjfMrPvpO5dKY0xR0hnPJF5G/BeSb+g2Bp4O8Ve1PSk4sDk/My2A9tTZAwU6tmvMbk/r1Y8tpn9EhgRj53OmYyfVV8YT6H2ILAoWWiOp9jYXDOO798X0j7T14CNZvaX7lCddMYTFjO73szmpdjB5cA9Zva7DMf4wuR8rqeBbZLemLquAB5lcn9erXjs9PeYn2lSf1b9YlydbyW9m2LfZgqwysxuGrc37xOSfgP4vxRxsHnv6U8p9tXuAM4FngQ+4FIVTyokXQ78dzN7j6TzKFZuM4D1wIfN7NVjOb5ukXQxhfHjeIrwvt+n+EKftJ+XpD8DPkhhjV9PkV12LpP8s+oHEVEQBEGjCENBEASNIoRaEASNIoRaEASNIoRaEASNIoRaEASNIoRaEASNIoRaEASNIoRaEASN4v8D4V9dKQOkU2UAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAATUAAAD8CAYAAAAWjzPMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztnX+wXVWV5z/LhB8SQn4RICSBFyBA0RQYOiKghT8QjLZFtEoRpsfJdOPkj9FpbXtKsa1qbB2msKtH2plyrMkgPXQX8mNoU2QsRoyoMBQCCUQwQAKPkHReCIQkENIoYMiaP+7Z56373rm557133699v5+qV2+ffc89d5993lt3rb32WsvcHSGEyIV3jPcAhBCik0ioCSGyQkJNCJEVEmpCiKyQUBNCZIWEmhAiKyTUhBBZMSKhZmbLzGyzmfWa2TWdGpQQQgwXG+7mWzObAjwDXAr0AeuAq9z9qc4NTwghhsbUEbz3fKDX3bcAmNltwHKgpVCbPn26z5kzZwQfKYQ4FHv27GH//v2WjpctW+a7d++u9d5HH330HndfNmqDGyNGItTmA9vDcR/wnkO9Yc6cOXz9618fwUcKIQ7Fdddd13S8e/du1q9fX+u9ZnbsaIxprBmJUKuFma0EVgLMnj17tD9OCDGAbovvHolQ2wEsDMcLir4m3H0VsAqgp6enu2ZXiAnAwYMHx3sIY8pIvJ/rgMVmtsjMDgeuBNZ0ZlhCiE7g7rV/cmHYmpq7HzCzLwD3AFOAm9z9yY6NTAjREXISWHUY0Zqau98N3N2hsQghRgEJNSFEVkioCSGyQkJNCJEN7t513k8JNSEyR5qaECIrJNSEEFkhoSaEyIbcNtbWQUJNiMyRo0AIkRXS1IQQ2SDzUwiRHRJqQoiskFATQmSFhJoQIhsUJiWEyA5pakKIrJBQE0JkhYSaECIruk2otS28YmY3mdkuM9sY+mab2Voze7b4PWt0hymEGA7JUVDnJxfqVJP6X8DAqs3XAPe6+2Lg3uJYCDEB6bZqUm2FmrvfD+wd0L0cuLlo3wx8osPjEkJ0iG4TasNdUzve3XcW7ReB4zs0HiFEh8lJYNVhxI4Cd3czazlrZrYSWAkwe/bskX6cEGII5KaF1WG4FdpfMrN5AMXvXa1OdPdV7r7U3ZdOnz59mB8nhBgu3WZ+DleorQFWFO0VwF2dGY4QotN0m/ezrflpZrcCHwCONbM+4FrgeuAOM7sa2AZcMZqDFEIMn5y0sDq0FWruflWLly7p8FiEEB0mN9OyDoooECJzuk2oDXdNTQgxSeiUo8DMlpnZZjPrNbNBG+7N7Agzu714/WEz6wmvfa3o32xmHwn9W83sN2b2azNb34n7laYmROZ0QlMzsynA94BLgT5gnZmtcfenwmlXA6+4+2lmdiXwbeAzZnYWcCXwB8CJwM/M7HR3f7t43wfdffeIB1kgTU2IjOlg7Of5QK+7b3H3t4DbaEQWRWKk0Z3AJWZmRf9t7v6muz8P9BbXGxUk1ITInA6Zn/OB7eG4r+irPMfdDwD7gDlt3uvAT83s0WKj/oiR+SlE5gzB/Dx2wLrWKndfNQpDirzP3XeY2XHAWjPbVMSbDxsJNUHDQmjQbZ6ybmAIz3S3uy9t8doOYGE4XlD0VZ3TZ2ZTgRnAnkO9193T711mtpqGWToioSbzU4jM6ZD5uQ5YbGaLzOxwGgv/awacEyONPgX83BsXXgNcWXhHFwGLgUfMbJqZTQcws2nAZcBGRog0NSEyplPVpNz9gJl9AbgHmALc5O5Pmtk3gfXuvgb4AfCPZtZLI13ZlcV7nzSzO4CngAPA5939bTM7HlhdWApTgR+6+09GOlYJtQyI5mNVu9XrdYnf4lXtVq+LiUGnnom73w3cPaDvr0L7DeDTLd57HXDdgL4twLkdGVxAQk2IzOm2LxoJNSEyR0JNTDiqzMd3vKPfxzNlypSyHfunTp066PXYrjJPq0zJt99+u+yL7QMHDgA0rdnE11O/zNPxQwHtQojskFAT40qV9hS1q6R9HXHEEWVfbB955JFl+/DDD2/6DXDYYYeV7aTVtdLUkqb1+9//vux76623BrXfeOONsu/NN98c1E4aHTRrctLgxoacEkDWQUJNiMzpti8MCTUhMkZramLMaLV3LJma0WSMJuW0adOafgPEgjZHHXXUoPdFk7OVU6GKZLZEkzGaosns/O1vf1v27d+/v2y//vrrTb/je6DZlE2fIZO083TbPLYNkzKzhWb2CzN7ysyeNLMvFv2zzWytmT1b/J41+sMVQgwVVZMazAHgL9z9LOAC4PNF0rdrgHvdfTFwb3EshJhgdJtQq1N4ZSews2jvN7OnaeRCWk6jyhQ0EsP9EvjqqIwyI6r2mUXz8J3vfCfQbFLOnDmzbM+aNWvQ69E8Td5R6Dcfo/cxtqOnsoo0rqrxxXHFa0bzMpmir7zyStn36quvDnod4He/+x3QbN5Gr11O/3RjSadiPycTQ1pTK3KOLwEeBo4vBB7Ai8DxLd6jCu1CjCPd9oVQO/WQmR0N/BPwJXd/Lb5WpBepnDlXhXYhxhWZnxWY2WE0BNot7v6jovslM5vn7jvNbB6wa7QGOdmp8m7GDbNR2Cfzcu7cuYP6oN+7GU2K6F3cs2fPoPbu3f01Lfbt21e2k9cyejejdzR91owZM8q+Y489tmzPmTOn6TdUe2Wj+RxN0ZdffnlQfzRJo3lc5R0V9ei2Oavj/TQaeZKedvfvhJdiQrgVwF2dH54QYqRIUxvMe4HPAr8xs18XfX8JXA/cYWZXA9uAK0ZniJOTqJ3Fxfu00B61n+OOO25QO2o3kaR9bd26tezbtGlT2d68eXPZTufs3Lmz7IsL9XU1tTiWefPmle2enh4AzjjjjLLvzDPPHPR61OTiPrqo1aX+Xbv6Ff6oVSZHQnRK5PSPOFrIUVCBuz8AtMoseElnhyOE6DTdJvwVUSBE5kioiRGRzM4qkxP6t7VEM+6EE04o28kkiwvm0aR85JFHAHjooYfKvieeeKJsb9u2rWy324c2UpKz4+STTy77zjnnnLJ9wQUXAHD++f11a6OpGs3StNcuhodFU3jv3r1AvxkKMkXr0m1zI6EmRMbk5gSog4SaEJkjoSaGTNU+tCqTE2D+/PkAnHjiiWVf3LP2wgsvAPCrX/2q7PvpT39atu+/v1HnNZqZrUhjiHvLqva8RVM5mnTJOxr3lsU9b8kUfOaZZ8q+2F63bh0AF198cdl32WWXle0LL7ywbKf5WLBgQdlXlUUkmaFxfKB9bIdC3k8hRFZ0m6CXUBsmUTuLGkXSuuI+tOgUSBpJXBB/7rnnynbSylavXl32Pfjgg2W76lv3lFNOKdtxIf60005r+kxo1tSSJhfHH6+fNLGoqSVNEqC3txdodmRs2bKlbCdt8pZbbin7nn/++bId96QlDS7eSxx3VW632E6B9AqCb0ZrakKI7JBQE0JkhYSaqEU0P2O+sRTEHUOf4j60ZJ5Gk3PNmjVl+9ZbbwVg48aNlZ+bFtLf/e53l31Lliwp24sWLSrbyUEQnRZV6bzbVZOKZl7cJ5acBtGk3LBhQ9lOjoK+vr6y74EHHijbMWQrmY+XX3552XfqqaeW7TSHMQV4zL2Wxhhf77Z/5lZ02zxIqAmRMYr9FEJkhzQ10ZJWBYajeZe8i9H8jNkokvcw7j1LJidUm53nnXde2X7/+98PNIcjRe9qHEsabzQfY+hU+gaPf/RVXt14r1WZNWLut4ULF5btlKXjvvvuK/see+yxsh3vNc1BTE2+fPnysp08oXFe4z61dF8qljyYbrt3CTUhMkdCTbQkajFxn1nMXJu0lpiDLAanp0iBuA+tSju76KKLyvall15ats8++2ygORg87jOLi+cpOmC4xUzaFYlJkQhRu4rB7UcffTQAxxxzTNkXz43779IcxHmJWlnaxxbnNWqIKftvVS1R6L5/7Ei33buEmhAZI0eBECI7pKmJQSQzLC6YRzOqqi5nJIYRJQdBNL0iySkQTc5zzz33kNePi//R/ErmZzTDqszPVo6CKvMzBrwn8zOa4rGd9snF8UdijdDkQIjzEvfcJbM27smLc5H2vMUiNFV71rrtHxy6757rFF450sweMbPHzexJM/vron+RmT1sZr1mdruZHd7uWkKIsafbCq/Uqfv5JvAhdz8XeBewzMwuAL4N3ODupwGvAFeP3jCFEMOhrkDLSajVKbziwL8Uh4cVPw58CPhXRf/NwDeA73d+iONPVYruuF+rKkdZrL+ZUnBDfz60aAbGHGJpH1rycg68fiKabtH8rPJ+xs8azqJxND9jO81HK/M2mehx/PG+XnutvyZ2ytgRQ6rSXEG/CXvSSSeVfdEDnD6jqlpWHGNO/7x16bZ7rlWh3cymFOXxdgFrgeeAV909LbD0AfNbvHelma03s/Vxa4MQYmw4ePBgrZ9cqOUocPe3gXeZ2UxgNXBmm7fE964CVgH09PRMyq+MpJ3EDLVVlcihX1OJdTljkZSqjLUxOD1FCkQtJJI0tKiptdqbFduJ4Xxrt9L00rXa7X2LTpV4XzEqIs1X1NTiXKU5XLp0adkXNcD0DOJziZpg0mBz+uetQydNSzNbBnwXmALc6O7XD3j9COAfgD8E9gCfcfetxWtfo7FE9TbwZ+5+T51rDodamlrC3V8FfgFcCMw0syQUFwA7RjoYIUTn6cSamplNAb4HfBQ4C7jKzM4acNrVwCvFOvsNNNbdKc67EvgDYBnw3wvrr841h0wd7+fcQkPDzN4JXAo8TUO4fao4bQVw10gHI4ToPB1yFJwP9Lr7Fnd/C7gNWD7gnOU01tcB7gQuscaC9HLgNnd/092fB3qL69W55pCpY37OA24upOo7gDvc/cdm9hRwm5n9J2AD8IORDmYiURW8Hs3PaHJG8yrtk9q0aVPZF+tyJmLa6rj3KgWnxwX56AhI7TomZ5V5GKkbJhWJ46oyb+O4qva5xTmMgfhpDh5//PGyL6YGT3MY5/X0008v2+kZxOcSi7SkPHDdGDrVofucD2wPx33Ae1qd4+4HzGwfMKfof2jAe9MafLtrDpk63s8ngCUV/VtoSFohxARliGFSx5rZ+nC8qlgTn1QookCIzBmCprbb3Ze2eG0HsDAcV62jp3P6ivX2GTQcBod6b7trDhkJtRZU5RWLZmbajwbN+9fS/rQYGlXl8YxVn2I4UMqHFvebVaWwbmVytvNEDsUUSefGuYjXrzJPI2msMbysVR66NAdxXqoqU8V5jXsBUxhVfC7xebVLXZ4zHbrPdcBiM1tEQ/BcSf8+1cQaGuvrv6Kx3v5zd3czWwP80My+A5wILAYeAazGNYeMhJoQmdMJoVaskX0BuIfG9oub3P1JM/smsN7d19BYV/9HM+sF9tIQUhTn3QE8BRwAPl9sE6PqmiMdq4SaEJnTKY3U3e8G7h7Q91eh/Qbw6RbvvQ64rs41R4qEWguimVKVjSKaNtEkSyZR3HwbvZfJ5EqFhqE/m0X83JgNo6o9lLTVrf6o63o/W2XxqMp8EcfV7l6qMnrEeYnmafJexnmN5mdKIx6fS7x+eobRlO8GcovrrIOEmhCZ021RFBJqLajapxa/+WNa66h9pFqYO3furLxu0khSIRFo1kiSplNV0xKq/0Db7UMb7jd1XU0ufn50BFTVDW11X2kO4rxEDXb79sZ2pjivaa6h/xnE5xKfVxpXO+dGjkhTE0JkhYSaECIbtKYmSqrMz2jaRDMrmp/79u0DmvN6RVJmiZhhIl4rORVaZcaoqtUZGcs/4Kp9bPHz01hb3Us0P1P4VJyX2E7mZ5zXNNfQ/wziXFY9L5mf+SOhJkTmSKgJIbJC3k8xiKpsE62yaKQU0jGVdCSF8USPZ7xWlXnZ7pt2uKbocLJ0VO1ZG8rnVJmn0D8HcV5iyFMizmtsJ69qq3ntRrMTtKYmhMgQCbUuptW3eeqv822fFr+rco1B/872VlrESIPQJyJ1NMkqbTgmCki0S1ceqaph2ur1yT7HhyLne6tCQk2IzJFQE0JkwxCTRGZBbaFWpPNeD+xw948XOZBuo5Gu91Hgs0We8UlLOzOpzjde2g8V90tFqmpxtgoYP1TfZKKdWQ/Vqcfj/r9Eq9xsVQw30D83uuU+E0OpJvVFGgVXEqrQLsQkoEOFVyYNdYsZLwD+CLixODYaFdrvLE65GfjEaAxQCDEyuk2o1TU//w74CpBK9cyhZoX2HKgyjWI7huOkvVVVe6ygf29Vyg828FpVaafbmZ+tvHjt9pENx6wdjnnc6l6q9ufFeana6xfnNbbTM2hXeLkb6bZ7r1P38+PALnd/dDgfYGYrzWy9ma3fv3//cC4hhBgmdbW0nARfHU3tvcDlZvYx4EjgGBpl4mea2dRCW2tZBaYosbUKoKenZ9LMXFU211a5wKKmNmPGDABmzpxZed1XXnml6ffAa6XF71bRC6ndyqM1lnuv2mltaayt7iUu9Kc5iPMS24k4r2muoX9PW3xGVc8rp3/eunSb97OtpubuX3P3Be7eQ6OQws/d/Y9RhXYhJgXdpqkNxfs5kK8CXy4qx8whswrtQuRCtwm1IW2+dfdfAr8s2llXaK8yP6vqb0JzEHVKQT1v3rzK66YU1C+88ELZFxfHp02bBrROF15lflYtuEP7BfzhBLRX9cXPr+pvl+MM4PXXXwea5yWm607EeY3pvpP5GZ0L8Xl1q/mZm8CqgyIKhMgcCTUhRFZIqAmg+Q8hmX/RnHnjjTfKdvTIzZkzB4Cenp6yL6Wqhn5Ts7e3t+yLZlbaexUzVFS129XajOe028fWjqFkL4kmZRprq3uJY0lzEOclmuVpDuO8prmGflM3Ppf4vNIz7LZ/cOg+76eEmhAZozU1UVKVoTVqAXFBOi7kJ+3hjDPOKPtOPvnksv3MM88AsHnz5rLv+eefL9tz584FWlcaT1pZnW/fqkiIoexja+cUSK+3CjJPToGqSunQrImlOYjzEklzGOc1amrpGcTnEp9Xu4I1OdNt9yyhJkTmSKgJIbJCQk0A1QvxscBKjGONZk7aZ3bmmWeWfeecc07ZTubnli1byr4NGzaU7YULFwLNJms036rMqLggHknjjibhcP7AWzkKqnLHxbGmxf2q8QPs3LmzbKc5iPMSSXMY5zXNNfTvc4vPJT6vbt6nJkeBECIruk2QS6gJkTkSamIQSX2P5kwyd6DZ5Jk+vZFyLu6nuuCCC8r2unXrANi2bdugvvi+o48+uuyL4UDRK5qI5mEM30r9rXKM1aVVlo3kyYyhT3FPXtVY9+zZU7afeOKJsh3nIBFN8DSHcV7jWNIziM8lPq9uM8EiEmpCiKyQUBODSH8UcT9a1Ahi3q8UXRD3UJ1/fn/c/8UXXwzALbfcUvb19fWV7fvuuw+AY445puw799xzy/asWbOAZi2oVY6yNN4YcVCVDbZVlEFVLc6qiIHoCIjtRJyfjRs3lu10r9A/B/Gz0lxB/xzGeY170tJnxOcSn1e3/WMntPlWCJEd3WZ6S6gJkTnS1MQg0h9FNOPi3rRXX321bCczKBYFiaE9l112GdAcGvXAAw+U7cceewyoXmQHOPvss4FmMywuzleZn9F5MJRiJFXmZ3QKVAWsx+snp0A0OdeuXVu2071GLrroorKd5gqa5zARzdr0DOJzic+r2/6xI9127xJqQmRMN66pjSSdtxBiEjAW6bzNbLaZrTWzZ4vfs1qct6I451kzWxH6/9DMfmNmvWb2X4vawpjZN8xsh5n9uvj5WLux1NLUzGwrsB94Gzjg7kvNbDZwO9ADbAWucPfB5X8yolVoUtyn9vLLLwPNITzRVLzwwgsB2LVrV9kXzddkqj344INlXzSpXnvtNaA59CqmuI6pxZMnMpphVZ7QVt7PqnTcVSFXMdtGDH1K+9Cil7PK5IR+s/qTn/xk2ZfmCvr3/8V9bmmuof8ZxOfSbRpKK8ZoHq4B7nX3683smuL4q/GEQmZcCywFHHjUzNYUcuP7wL8DHgbuBpYB/7d46w3u/rd1BzIUTe2D7v4ud1864CYWA/cWx0KICcbBgwdr/YyQ5cDNRftm4BMV53wEWOvuewtBthZYZmbzgGPc/SFvSOB/aPH+WozE/KxzE0KIcWQMixkf7+5JTX8ROL7inPnA9nDcV/TNL9oD+xNfMLMnzOymVmZtpK6jwIGfmpkD/6MoUFznJrKiVQrtaH5VeT+jJ/PEE08Emj170by89dZbgWaPYTTZktm6devWsm/JkiVle9GiRWU7hVdFkzR6SpN52SpxZPr2jvcaN7emFNzRkxszjqTQp7i5OJJMToCrrroKaJ6XNFfxc6PZHr2f6RnI4zmYIczDsWa2PhyvKv7XATCznwEnVLzv6wM+zwtZ0Qm+D3yLhgz6FvBfgD891BvqCrX3ufsOMzsOWGtmm+KLh7oJM1sJrASYPXt2zY8TQnSKIQi13WF5qeo6H271mpm9ZGbz3H1nYU7uqjhtB/CBcLyARsnNHUU79u8oPvOl8Bn/E/hxu5uoJdTcPX3ALjNbTaPeZ52boJD0qwB6enqy+eqMfyhxH1hasI4aRQwdWrCg8exOOeWUsu/yyy8v20mrW716ddkXnQZJ64naz+OPP162436u0047DWjWeFKYFfRrcK3qhibtJ2pEsS5nKpISU3BX5UOL14/70KJTIGlocV5iQPqLL74INM9rdNCkZyDtbDBjNCdrgBXA9cXvuyrOuQf4z8GEvAz4mrvvNbPXzOwCGo6CfwP8N4AkY4rzPwlsHHjRgbRdUzOzaWY2PbWLgWwMN8EhbkIIMY6kJJFj4Ci4HrjUzJ4FPlwcY2ZLzezGYix7aZiQ64qfbxZ9AP8euBHoBZ6j3/P5N8VWjyeADwJ/3m4gdTS144HVxbrLVOCH7v4TM1sH3GFmVwPbgCtqXEsIMcaMhabm7nuASyr61wOfC8c3ATe1OO/siv7PDnUsbYWau28Bzq3or7yJbqFqQR36TaZ9+/aVfXFvVzLFokl46qmnlu3ly5cDcNxxx5V9cfH//vvvB5rzsUWTL7aTeRnzsUXzs6rGaMxskbJgRPMz1iiNDpIqUj60mG0jOgLiPrQ0H9HkjKZu2v8W57UqX5rMz8F025woTEqIzJFQE0JkhYSaGDJV+9eiabZ3795B74km6wkn9G/9qdrHFtNap4SRDz30UNkX02JHszSNYfv2/v2OsT1S0p63OL4YvpVScMckmdE7m0KfoH8fWvJyQnPIVZrDOK/ak9aebgxol1ATInOUJFKMiKrU31VaW9QyYhB2chCktODQHDFw0kknAbB0af8eyU2b+vdCxz1jKeogajwxeD45AuJYolMjORLiWGLwfCqCErWvWJczvR4D+iMxOD3tP4v70KJTIM2hUnQPnW6bJwk1ITJHQk0IkQ1aUxMdI/4hVe39iiZfDLNKr8+dO7fsi3vLkikX+04//fSyHU261I57y6JJV9f8nDFjRtkX97ylsUTzMuaRS3vyqqo+QXM+tNQfQ5/iPrQ0xm77B+0E3TZnEmpCZI4cBUKIbJD5KUaFqn1sraoeJZMr5i2LHstkdsY9XjFfW9wztnDhQqDZ/I3taPZWkSpHxTCq2E7mZasiz8mUrKr6FF+Hfu/mcCtfidZ029xJqAmRORJqYlRJf2BRC4n71JLWFvuqtLa4IB+1tqqMu7FWZ3QExIy4VaQxVjkyoF/bjH1R+0rjjuOPGmrVfcd/wG77Zxwtum0eJdSEyBwJNSFENqQkkd2EhNo40crMSu1WYVTJ1Ev1P6E5YD46DVIa8ZhOPJqidQuvRPMzjiW1o0kZ95aldnQktApC7zZtYizptrmVUBMicyTUhBBZIaFWgZnNpFEU4Wwa9ff+FNgM3A70AFuBK4qqy2IEpD/AVqZZMt+iSRizgMTKTWlPWfR4xnYyO1uZn1WmcGwnszKu2cTXq1Jsd9s/2HjTjZtv61Zo/y7wE3c/k0a9gqeBa4B73X0xcG9xLISYYIxRhfYJQ1tNzcxmABcD/xbA3d8C3jKz5fQXJr2ZRlHSr47GILudKk2nlXYUta60kB/7Ynsknx/b0sQmNt3m/ayjqS0CXgb+3sw2mNmNRf3P40OR0RdplNITQkwwuk1TqyPUpgLnAd939yXA6wwwNb0xI5WzYmYrzWy9ma2Pu82FEKNPXYGWk1Cr4yjoA/rc/eHi+E4aQu2lVBLezOYBu6re7O6rgFUAPT09+czcBGKk5l8rR4HIg257pm01NXd/EdhuZikR/SXAU8AaYEXRtwK4a1RGKIQYEdLUqvkPwC1mdjiwBfgTGgLxDjO7GtgGXDE6QxRCjIRucxTUEmru/mtgacVLl3R2OGI8yOlbWjSTmxZWB0UUCJE5EmpCiKyQUBNCZIWEmhAiKyTUhBDZoCSRQojskKYmhMgKCTUhRFZIqAkhskGbb4UQ2SGhJoTIim7zftZN5y2EmKSMRZYOM5ttZmvN7Nni96wW560oznnWzFaE/uvMbLuZ/cuA848ws9vNrNfMHjaznnZjkVATImPGMElk25olZjYbuBZ4D3A+cG0Qfv+n6BvI1cAr7n4acAPw7XYDkVATInPGSKgtp1GrhOL3JyrO+Qiw1t33FpXn1gLLijE+FMoDtLruncAl1qbQhtbUhMicIQisY81sfTheVWSurkOdmiXzge3huK/oOxTle9z9gJntA+YAu1u9QUJNiMwZgqNgt7tX5U0EwMx+BpxQ8dLX44G7u5mNm8tVQk2IjOnkPjV3/3Cr18ysTs2SHfSX1QRYQKO05qHYASwE+sxsKjAD2HOoN2hNTYjMGaM1tTo1S+4BLjOzWYWD4LKir+51PwX83NsMtq1QM7MzzOzX4ec1M/tSXReuEGJ8GSOhdj1wqZk9C3y4OMbMlprZjcU49gLfAtYVP98s+jCzvzGzPuAoM+szs28U1/0BMMfMeoEvU+FVHUhb89PdNwPvKj54Cg11cDX9Ltzrzeya4lgV2oWYYIxFRIG776GiZom7rwc+F45vAm6qOO8rwFcq+t8APj2UsQzV/LwEeM7dt1HPhSuEGGdUIu/QXAncWrTruHCFEONINyaJrK2pFTU/Lwf+98DXioVu/9O5AAAFF0lEQVS7SlFvZivNbL2Zrd+/f/+wByqEGB7dpqkNxfz8KPCYu79UHL9UuG45hAsXd1/l7kvdfen06dNHNlohxJCRUGvNVfSbnlDPhSuEGGck1Cows2nApcCPQnelC1cIMXEYw4D2CUMtR4G7v04j3ir2VbpwhRATi5wEVh0UJiVE5nSb91NCTYjMkaYmhMiG3NbL6iChJkTmSKgJIbJCQk0IkRVyFAghskFrakKI7JBQE0JkhYSaECIrJNSEEFkhoSaEyIZuTBIpoSZE5khTE0JkhYSaECIrJNSEENmgzbdCiOyQUBNCZEW3eT/r1ij4czN70sw2mtmtZnakmS0ys4fNrNfMbi9K6AkhJhjdVqOgrVAzs/nAnwFL3f1sYAqNosbfBm5w99OAV4CrR3OgQoih042FV+qWyJsKvNPMpgJHATuBDwF3Fq/fDHyi88MTQowUCbUBuPsO4G+Bf6YhzPYBjwKvuvuB4rQ+YH7V+1WhXYjxRUJtAGY2C1gOLAJOBKYBy+p+gCq0CzG+HDx4sNZPLtTxfn4YeN7dXwYwsx8B7wVmmtnUQltbAOxod6Ft27btXrly5evA7hGMeaJyLLqvyUKO9wSN+5o2oO+eor8OWcxJHaH2z8AFZnYU8DsaBYzXA78APgXcBqwA7mp3IXefa2br3X3p8Ic8MdF9TR5yvCco76sn9rl7basqF+qsqT1MwyHwGPCb4j2rgK8CXzazXhrV238wiuMUQoha1Np86+7XAtcO6N4CnN/xEQkhxAiou6Wjk6wah88cC3Rfk4cc7wnyva8hYTm5coUQYjw0NSGEGDXGVKiZ2TIz21zEi14zlp/dKcxsoZn9wsyeKuJhv1j0zzaztWb2bPF71niPdTiY2RQz22BmPy6OJ32Mr5nNNLM7zWyTmT1tZhdO9ueleOzWjJlQM7MpwPeAjwJnAVeZ2Vlj9fkd5ADwF+5+FnAB8PniPq4B7nX3xcC9xfFk5IvA0+E4hxjf7wI/cfczgXNp3N+kfV6Kxz40Y6mpnQ/0uvsWd3+Lxv625WP4+R3B3Xe6+2NFez+Nf5D5NO7l5uK0SRkLa2YLgD8CbiyOjUke42tmM4CLKbYcuftb7v4qk/95KR67BWMp1OYD28Nxy3jRyYKZ9QBLgIeB4919Z/HSi8Dx4zSskfB3wFeAFDMzh5oxvhOYRcDLwN8XZvWNZjaNSfy8RhqPnTtyFAwTMzsa+CfgS+7+WnzNGy7lSeVWNrOPA7vc/dHxHkuHmQqcB3zf3ZcArzPA1Jxsz2uk8di5M5ZCbQewMBzXihediJjZYTQE2i3u/qOi+yUzm1e8Pg/YNV7jGybvBS43s600lgY+RGMtamZh4sDkfGZ9QF8RGQMN8+w8JvfzKuOx3f33QFM8dnHOZHxWHWEshdo6YHHhoTmcxsLmmjH8/I5QrDP9AHja3b8TXlpDIwYWasbCTiTc/WvuvqCIHbwS+Lm7/zH9Mb4wOe/rRWC7mZ1RdF0CPMXkfl5lPHbx95juaVI/q04xpptvzexjNNZtpgA3uft1Y/bhHcLM3gf8PxpxsGnt6S9prKvdAZwEbAOucPe94zLIEWJmHwD+o7t/3MxOoaG5zQY2AP/a3d8cz/ENFTN7Fw3nx+E0wvv+hMYX+qR9Xmb218BnaHjjNwCfo7GGNqmfVSdQRIEQIivkKBBCZIWEmhAiKyTUhBBZIaEmhMgKCTUhRFZIqAkhskJCTQiRFRJqQois+P/Ok9yYxZOn6QAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "#NBVAL_SKIP\n", "\n", "# Let's see what we got....\n", "plot_image(v4[0].data[-1])\n", "plot_image(v4[1].data[-1])\n", "plot_image(p4.data[-1])" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [], "source": [ "norm_p = norm(p)\n", "assert np.isclose(norm_p, .35098, atol=1e-4, rtol=0)" ] } ], "metadata": { "hide_input": false, "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "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.7.6" }, "latex_envs": { "LaTeX_envs_menu_present": true, "autoclose": false, "autocomplete": true, "bibliofile": "biblio.bib", "cite_by": "apalike", "current_citInitial": 1, "eqLabelWithNumbers": true, "eqNumInitial": 1, "hotkeys": { "equation": "Ctrl-E", "itemize": "Ctrl-I" }, "labels_anchors": false, "latex_user_defs": false, "report_style_numbering": false, "user_envs_cfg": false }, "widgets": { "state": {}, "version": "1.1.2" } }, "nbformat": 4, "nbformat_minor": 2 }