{ "cells": [ { "cell_type": "markdown", "metadata": { "toc": true }, "source": [ "

Table of Contents

\n", "
" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Simple contour plot with lines\n", "\n", "This produces a simple oldschool contour plot with labels" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "nbgrader": { "grade": false, "grade_id": "cell-e58a84ffce31f93a", "locked": true, "schema_version": 1, "solution": false }, "scrolled": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "setting context.data_dir to /Users/phil/repos/eosc213_students/notebooks/cookbook_examples\n" ] } ], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np\n", "\n", "import context" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "lines_to_next_cell": 0 }, "outputs": [], "source": [ "data = np.load(context.data_dir / \"contourline.npz\")\n", "x1, y1, psi = data[\"x1\"], data[\"y1\"], data[\"psi\"]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## make the plot with specified lines and labels\n", "\n", "They don't have to be the same -- you can label only some\n", "of the lines. You can control line thickness, dashed/solid/dotted etc.,\n", "see the matplotlib docs." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3QAAADKCAYAAADtuAIjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzs3XlYVdX6B/DvZp4UIecBBRQQAUEGZ8XEeTbTrExTK4fM7GqaZXkr7erNqcwhh2t6cypNvQ6kqIgDauKMihNOIIhaCIjA4Xx/fyj7h4pwwMM5B3s/z7Mf4Zy911ob3k7rZa29lkISQgghhBBCCCHKHjNjN0AIIYQQQgghRMlIQieEEEIIIYQQZZQkdEIIIYQQQghRRklCJ4QQQgghhBBllCR0QgghhBBCCFFGSUInhBBCCCGEEGWUJHRCCCGEEEIIUUZJQieEEEIIIYQQZZQkdEIIIYQQQghRRlkYuwEFqVChAuvWrWvsZgjxTBkZGbC3tzd2M4R4JolRYeokRoWpkxgVxhYTE3ObZKWizjPJhK5KlSo4cuSIsZshxDNFRkYiNDTU2M0Q4pkkRoWpkxgVpk5iVBiboihXdTlPplwKIYQQQgghRBklCZ0QQgghhBBClFGS0AkhhBBCCCFEGSUJnRBCCCGEEEKUUZLQCSGEEEIIIUQZJQmdEEIIIYQQQpRRktAJIYQQQgghRBklCZ0QQgghhBBClFGS0AkhhBBCCCFEGSUJnRBCCCGEEEKUUZLQCSGEEEIIIUQZJQmdEEIIIYQQQpRRktAJIYQQQgghRBklCZ0QQgghhBBClFEWRZ2gKMpSAF0B3CLp8+i1NQA8H51SAcBfJP0LuPYKgDQAuQA0JIP01G4hhBBCCCGE+NsrMqEDsAzAXADL814g2S/va0VRZgBILeT6NiRvl7SBQgghhBBCCCEKVmRCRzJKUZQ6Bb2nKIoCoC+Al/XbLCGEEEIIIYQQRXneZ+haAkgmeeEZ7xPAdkVRYhRFefc56xJCCCGEEEIIkY8uUy4L0x/AqkLeb04yUVGUygB2KIpyjmRUQSc+SvjeBYBKlSohMjLyOZsmROlJT0+XGBUmTWJUmDqJUWHqJEZFWaGQLPqkh1MuN+ctivLoNQsACQACSd7QoYzJANJJflvUuZ6enoyLiyuyXUIYS2RkJEJDQ43dDCGeSWJUmDqJUWHqJEaFsSmKEqPLopLPM+UyDMC5ZyVziqLYK4pSLu9rAO0BnH6O+oQQQgghhBBC5FNkQqcoyioA0QA8FUW5oSjKkEdvvYYnplsqilJdUZStj76tAmCfoignABwGsIVkuP6aLoQQQgghhBB/b7qsctn/Ga8PKuC1RACdH319GUDD52yfEEIIIYQQQohneN5VLoUQQgghhBBCGIkkdEIIIYQQQghRRklCJ4QQQgghhBBllCR0QgghhBBCCFFGSUInhBBCCCGEEGWUJHRCCCGEEEIIUUZJQieEEEIIIYQQZZQkdEIIIYQQQghRRklCJ4QQQgghhBBllCR0QgghhBBCCFFGSUInhBBCCCGEEGWUJHRCCCGEEEIIUUYVmdApirJUUZRbiqKczvfaZEVREhRFOf7o6PyMazsqihKnKMpFRVEm6LPhQgghhBBCCPF3p8sI3TIAHQt4fRZJ/0fH1iffVBTFHMAPADoB8AbQX1EU7+dprBBCCCGEEEKI/1dkQkcyCsDdEpQdAuAiycskswGsBtCjBOUIIYQQQgghhCjA8zxD976iKCcfTcl0KuD9GgCu5/v+xqPXhBBCCCGEEELogUUJr5sP4CsAfPTvDACDnzhHKeA6PqtARVHeBfAuAFSqVAmRkZElbJoQpS89PV1iVJg0iVFh6iRGhamTGBVlRYkSOpLJeV8rirIIwOYCTrsBoFa+72sCSCykzB8B/AgAnp6eDA0NLUnThDCIyMhISIwKUyYxKkydxKgwdRKjoqwo0ZRLRVGq5fu2F4DTBZz2B4B6iqK4KopiBeA1AJtKUp8QQgghhBBCiKcVOUKnKMoqAKEAKiqKcgPAFwBCFUXxx8MplFcAvPfo3OoAFpPsTFKjKMr7AH4HYA5gKcnYUrkLHWRlZSEpKQmJiYm4efMmypcvj8DAQDg5FfT4nxBCCCGEEEKYviITOpL9C3h5yTPOTQTQOd/3WwE8taWBIcTGxmLcuHG4ceMGEhMTcefOnQLPc3V1hb+/P9zc3FCnTh24urrC1dUVNWrUgIVF4T8eOzs7KEpBjwoKIYQQQgghROkr6aIoJm/16tX4/fff0a1bN7Ro0QLVqlVD9erVUa1aNVSrVg23b99GTEwMYmJicOrUKWzduhVZWVnFqsPZ2RmNGjVCYGAggoKCEBgYiDp16kiSJ4QQQgghhDCIFzahO3v2LOrWrYsNGzY885x27dqpX2u1WiQnJ+PKlSuIj49HYmIiyIIX5UxPT0dCQgJu3LiB+Ph4REZGQqPRAABsbW1Rp04d1K5dG25ubggICEBQUBAaNGgAS0tL/d6kEEIIIYQQ4m/thU3o7ty5gypVquh8vpmZmTp617Rp08fei46Oxo4dO3D06FEcPXoU169ff0YpwIMHD3Dx4kVcunQJ27dvh1arBQBYWlrC19cXISEhCAwMhL+/P7y9vWFnZ1eyGxRCCCGEEEL87b2wCV1WVhbs7e2fuxySePnll/HgwQN4eXmhZcuWaNSoEQICAmBvb4/U1FSkpqbir7/+Ur/OOxISEnDkyBGkpqYiJycHR48exbFjxx4b+XN1dUXTpk0RHByMwMBABAQEwMHB4bnbLYQQQgghhHjxvbAJXfny5fHnn38+dzmKomDs2LH4+uuvMW3aNHTv3r1Y15PElStXcPLkSZw6dQonT57EsWPHcPnyZWi1WsTHx+PKlStYuXKlek3t2rXRuHFjBAYGwtXVFWZmz95dQlEU1K1bF97e3kUu4iKEEEIIIYR4sbywGUClSpVw9uxZvZQ1adIkbNy4EWPGjEHHjh1hZWWl87WKoqgrZ/bo0UN9PSsrC+fOncPRo0cRExOD6OhonD59GtnZ2bh69SquXr2KtWvX6lyPpaUlvLy80Lx5c3VapyR5QgghhBBCvNhKtLF4WeDi4oKEhAR1sZLnYWVlhenTp+Py5ctYtGiRHloHWFtbo2HDhkhMTMTy5ctx6tQpZGdnl7i8nJwcnDp1CgsWLMDgwYPRsGFD2NjYoE6dOujWrRtmzpyJvXv3Ij09XS/tF0IIIYQQQhjfCzt84+rqitzcXNy4cQN16tR57vI6dOiApk2bYsaMGXjvvff0NvKVmZmJtLQ0vP3223BxcYGNjQ2sra3Vfy0sLArdBiEnJwdnz57FkSNHcObMGdy6dUt9Lzc3Vx3t27x5s/p61apVERQUhNDQUISEhMDPzw+Ojo56uR8hhBBCCCGE4bywCZ27uzsA4OLFi3pJ6BRFwbhx49C7d29s2LABffr0ee4yAeDtt9/G1KlTYW9vj8mTJz93eZmZmbhw4QLat2+P5OTkAs9JSkrC5s2bH0vyKlasiICAALRp0wZNmjRBo0aNJMkTQgghhBDCxL2wCV3dunUBAJcuXUJYWJheyuzevTvq1KmD7777Tm8Jnbu7O0aNGoXvvvsOr776Klq1avVc5dna2sLPzw9r165F//79kZiYiK5du8LPzw/79u1DbGws7ty589R1t2/fxo4dO7Bjxw71tYoVK8LPzw+tWrVCvXr1Ch0ptLS0RIMGDeDh4QFzc/PnugchhBBCCCGEbl7YhK5GjRqwsbHBhQsX9Famubk53n//fYwdOxbHjx+Hv7+/XsqdOnUqtmzZgjfffBPHjx+Hs7Pzc5fZqlUrnDhxAgMGDMDmzZvRrl077NmzBwCQmpqKkydP4ujRo4iKikJMTAyuX7+u7pmX5/bt29i1axd27dqlc73m5uZwc3NDkyZN0LZtWwQHB8PT01OSPCGEEEIIIUrBC7soipmZGerWrYvz58/rtdzBgwfDzs4O3333nd7KtLe3x+rVq5GUlIS33377sX3qnkfFihWxadMm9OrVC6NHj8bChQsBAI6OjmjZsiVGjx6NdevW4cqVK8jMzMSJEyewbNkyDB48GL6+vrC1tS12nbm5ubhw4QJWrFiBQYMGoUGDBrC2toaLiwt69uyJuXPn4vjx48jJydHLPQohhBBCCPF3VuQInaIoSwF0BXCLpM+j1/4NoBuAbACXALxN8q8Crr0CIA1ALgANySD9Nb1o9erV09vWBXmcnJwwcOBALFmyBN988w2qVKmil3KDgoIwffp0jBkzBrNmzcJHH32kl3ItLS2xevVq9O7dG8OHD0fFihXxyiuvPHWelZUV/Pz84Ofnh4EDBwJ4uIfezZs3ceLEiccWW3lSRkYGTp8+jZMnT+LSpUtISUlBbm6u+n5ubi6uX7+O69evY+PGjerr1atXR3BwMNq1a4fmzZvD29u7WFtCCCGEEEII8Xeny5TLZQDmAlie77UdAD4hqVEUZRqATwCMf8b1bUjefq5WlpCHhwc2b94MjUaj1/3YRo8ejfnz52P+/Pl6Wcgkf7l79uzB+PHj0bRpUzRt2lQv5VpZWWHt2rUICwvD66+/jg0bNqBTp05FXqcoCqpXr47q1asXel5YWBh27txZ7HYlJiZi48aNjyV5VapUQWBgINq1a4dWrVrBx8dHkjwhhBBCCCGeocgsh2SUoih1nnhte75vDwLQzwoheubh4YGcnBxcvXpVXfVSHzw9PdG1a1f88MMPGD9+fImmJhZEURQsXboUgYGB6NevH44dO4aXXnpJL2Xb2dlhy5YtCAsLQ69evbBhwwZ07NhRL2X/+9//RqtWrdQ97lq2bIlPP/0UTk5OAB6O9Gm1Wly6dAkRERE4ePAg4uPjC9x3Lzk5GVu3bsXWrVvV16pUqYJGjRohLCwMrq6uhbbF1tYWDRs2RLVq1fRyb0IIIYQQQpgyfQxbDQaw5hnvEcB2RVEIYCHJH/VQn848PT0BAHFxcXpN6ABg7NixCA0NxU8//YRhw4bprVwnJyesXbsWzZs3Vxc0MTPTz6OOTk5O2LFjB8LCwtCjRw+sXr0avXr1eu5yAwICkJaWhqNHj2LGjBlYs2YNBg4ciPnz5z9WftOmTfHmm28CeJjk3bhxAzExMYiIiMC8efOe+exgcnIytm3bhm3btuncJisrK3h4eKBly5bo0KEDQkJCJMkTQgghhBAvHEWXBTgejdBtznuGLt/rnwIIAtCbBRSkKEp1komKolTGw2mao0hGPaOOdwG8CwCVKlUKXLt2bTFv5Wmpqano2bMnhg8fjr59+z53efmRxIgRI5CWloaffvpJ76s4bty4EbNnz8bgwYMxYMAAvZadlpaGTz75BGfPnsW4ceP0NlKX58KFC5g+fTouXryIl19+GWPGjIGDg0Oh19y+fRvlypXDlClTcPnyZWRkZCAzMxPZ2dl6WyTGwsICVapUgY+PD0JCQuDj44NKlSoVuh3Ds6Snpxd5T0IYk8SoMHUSo8LUSYwKY2vTpk2MLmuQlDihUxRlIIBhANqSvK9DGZMBpJP8tqhzPT09GRcXV2S7dJG3CEjeCo/6tG7dOvTp0wdr167Fq6++qteySWLAgAFYuXIlwsPD0b59e72Wn56ejl69eiEiIgLTp0/HuHHj9Fp+Tk4Opk2bhsmTJ6NmzZpYuXIlmjVrVug1JNG5c2eEh4fDwsICXbp0gZWVFUgiPT0dKSkpyMrKKvC6Bw8eICMjA2lpacjMzHxqC4ZnMTc3h6urK5o3b45OnTqhWbNmqFmzZpFJXmRkJEJDQ3WqQwhjkBgVpk5iVJg6iVFhbIqi6JTQgWSRB4A6AE7n+74jgDMAKhVyjT2Acvm+PgCgoy71eXh4UF+aN2/Oli1b6q28/DQaDT08POjv70+tVqv38tPT0+nr60snJydeunRJ7+U/ePCA/fr1IwCOGTOGubm5eq8jOjqarq6uNDc355dffkmNRlPo+dnZ2Rw+fDgBsEWLFjx79myRdXz33XfEw+m9BEBHR0eamZk99lpxDnNzc7q6unLAgAH8+eefefXq1ad+v7t3736eH4sQpU5iVJg6iVFh6iRGhbEBOEIdciddti1YBSAUQEVFUW4A+AIPV7W0BrDj0UjGQZLDFEWpDmAxyc4AqgD47dH7FgBWkgwvMsPUs/r16z+2iqI+mZubY8KECRg8eDC2bduGzp0767V8e3t7/PbbbwgODkavXr1w4MAB2Nvb6618a2trrFy5ElWqVMGsWbOQkJCAn376CTY2Nnqro0mTJjh27BhGjBiBzz//HNu3b8eKFStQp06dAs+3tLTEvHnz0KRJE3z44Ydo2LAhJk2ahHHjxsHa2rrAa/r06YOlS5fi+PHjMDc3R/PmzTF48GB07NhRHbVLS0vDvXv3cO7cOWzatAlHjhzB7du3H9teIU9ubi7i4+MRHx+PFStWAHj4u3ZxcUGzZs3QsWNHJCUlFbqXXoUKFeDr66vXn6UQQgghhBBP0SXrM/ShzxG6GTNmEABTUlL0VmZ+2dnZdHFxYZMmTUpllI4kw8PDaWZmxj59+pRKHVqtlv/+97/VUbHbt2/rvQ6S/O9//8ty5crRwcGBixYtKvJekpKS2LdvXwKgh4cHw8PDCz3/ypUr/OSTT1izZk0CYP369bly5cpCRwVTU1P566+/8tVXX2XNmjVpYWFR4pG9gg5HR0e2aNGCkydP5sGDB5mZmVmin50QxSV/WRamTmJUmDqJUWFs0HGEzujJW0GHPhO6bdu2EQD37NmjtzKfNH/+fALg9u3bS62Ob7/9lgA4efLkUqtj9erVtLa2Zr169Xj+/PlSqSM+Pp5t2rQhAHbq1InXr18v8ppt27axbt26BMAePXrw4sWLhZ6v0Wi4Zs0aNmjQgADo5eXFX375RedkOCMjg1u3bmWfPn1YoUKF55q+WdBRrlw5BgcHc8KECYyKimJGRoZO7RKiOKQjIkydxKgwdRKjwtgkoXvk6tWrBMD58+frrcwnPXjwgDVr1mSzZs1KbZROq9Vy4MCBBMA1a9aUSh0kuW/fPlasWJHOzs6MjIwslTpyc3M5Z84c2trasnz58jqN1j148IBTp06lvb09raysOH78eN67d6/IetauXUtvb28CYEhISLHuKTc3l5GRkfT29ua5c+e4d+9evv/++/T19WXNmjVZr169Zx5OTk7FSvLs7e0ZFBTECRMmcN++fbx//77O7RSiINIREaZOYlSYOolRYWyS0D2i1WpZrlw5jhw5Um9lFmTevHkEUOS0wOfx4MEDtmjRgjY2NoyOji61ei5evEgvLy9aWlpy0aJFpVpP3mjdyy+/XOTIG0kmJCSoiW2VKlW4ZMmSIhda0Wg0XLp0qToVs3Pnzjx+/Pgzz8/JySFJZmVlkXw4bXfMmDGPnVPYh/wff/xBKyur5x7Js7e3Z2BgIMePHy8jeaLYpCMiTJ3EqDB1EqPC2HRN6PSzY7UJUxQFDRo0QGxsbKnWM3jwYNSuXRuTJk16mCmXAmtra6xfvx41atRAt27dcOnSpVKpx93dHdHR0WjTpg3eeecdjB49GhqNplTqiYiIwIIFC3DkyBH4+PjgX//6V6GLjVSvXh3Lli3DoUOH4ObmhiFDhiA4OBh79ux55jXm5uZ4++23cf78eUybNg0HDhxAQEAAZs+erZ5z/vx5/Prrr7hy5QqSk5MBPNyc/M8//0R0dDRq1Kih8315e3vjww8/RN26dR97vXLlyvjggw8we/ZsfPDBB/Dy8nrmQi8AkJGRgZiYGEybNg2tWrWCvb097O3tERQUhAkTJmDXrl04fPhwoce1a9dKLR6FEEIIIYQJ0CXrM/ShzxE6khw6dChfeumlUpsOmWfp0qUEwPXr15dqPXFxcXR2dma9evVKbbEX8uFI1YcffkgAbNu2banWdePGDfbq1YsA6OPjw3379hV5jVar5apVq1irVi0CYM+ePXV69u/u3bucMGECT58+rb42e/ZsKorCt956i23btmVQUBAHDhzI4OBgvvfee0xOTn6sDF3+aqfVahkbG8upU6eycePG6shbWFgYf/nlF2ZnZ5MkMzMzefDgQX7yySf09fWljY2NXp/ZUxSFNWrU4CuvvMJly5YVuA2DePHIX5aFqZMYFaZOYlQYG2TK5f/L26csISFBr+U+KScnh15eXvTy8lKn7ZWWffv20cbGho0bNy71qXhLly6ltbU1a9euzZiYmFKta+PGjXRxcSEADhkyRKck8v79+5wyZQodHBxoYWHB999//6kETBc9e/bk4sWLSZJHjhzhjh07+Mcff/DmzZsk+dg+fSX5kL9+/Tq/+uor9f6qVq3KSZMmFRiXWVlZjImJ4ZdffsmAgADa2dnpPcmrWrUqe/TowUWLFjE+Pl6SvBeMdESEqZMYFaZOYlQYmyR0+ezevbvUn2/Ls379egLgjz/+aJC6zMzM2LlzZ3W0p7QcPnyYtWrVorW1tZr0lJa0tDSOGzeOFhYWdHJy4oIFC4p8To4kb968yffee4/m5uZ0cHDg5MmTi1w4Jb+DBw/Sy8uLaWlphZ6XnZ39XB/yGo2GmzdvZpcuXagoCi0sLNi/f39GR0cXmlRpNBqePn2a//73v9m0aVM6OTnR0dGx0MPc3LxYiV7lypXZrVs3SfJeANIREaZOYlSYOolRYWyS0OVz584dAuC//vUvvZZbEK1Wy2bNmrFq1apFJgb68OOPPxIA+/fvr1PS8zxSUlLYrl07AuCgQYNKfWTw9OnTDA0NJQAGBgbywIEDOl137tw5vvLKKwTASpUqcfbs2Trv/zZ37lxGRUUV+J5Wq+WRI0f4/fffMyIiQuf7KMzFixc5ZswYli9fXr3PxYsX6/Vnm5qayoiICI4YMYL16tUr9oItlStXZteuXfnjjz9KkleGSEdEmDqJUWHqJEaFsUlC9wQXFxf2799f7+UWJDo6mgA4adIkg9T3zTffEADffffdUu9sazQaTpo0iYqisEGDBjxz5kyp1qfVarly5UpWr16dAPjmm2/yxo0bOl176NAhvvzyywTAWrVqcdGiRTqNZK5fv57nzp0r8L2UlBQOHTqUPXr0KNZ9FCUtLY0//PCDuneeo6MjP/jgg8ee89On7OxsHj9+nP/85z/p7+9PW1vbEiV5CxYs4OnTpxkXF1foIdswGJ50RISpkxgVpk5iVBibJHRP6N69O728vPRe7rO89tprtLGx4ZUrVwxS3yeffEIAHD16tEFGULZv385KlSrRzs6OixcvLvU609LS+Mknn9Da2pp2dnb84osvmJ6ertO1ERER6qIk7u7u/Omnnwp9xjE2NparV69Wv9dqtepBPkyG3NzcePjw4ee7qQJotVru3buXr7/+ujqS1rhxYy5cuJB//fWX3uvLL29K59SpUxkYGKj35/ZsbGzo5+fHsWPHMioqSpK8UiYdEWHqJEaFqZMYFcYmCd0TvvjiCyqKonMS8LyuXr1KW1tb9unTxyD1abVajh49mgA4btw4gyR1iYmJbNu2LQGwX79+/PPPP0u9zsuXL7Nfv34EwGrVqnHp0qU6TTXVarXcuHEjGzZsSACsW7culy1bVuiIXXx8PB88eFDgex06dHjm1Ex9uXXrFmfOnKmO2tna2nLAgAHcuXPnYwu0lKbc3FyeOXOG06ZNY3BwMO3t7fWe5Pn4+HD06NHcvXu37LWnR9IREaZOYlSYOolRYWyS0D1h48aNBKDTcvj68uWXXxIAd+zYYZD6tFotR4wYYdCkTqPRcOrUqTQ3N2etWrUM9uG3f/9+ddTNx8eHmzdv1ul+c3Nz+dtvv6mJnaurK3/88Ud1E/E8WVlZDA0Npb29PYcOHcpevXpx/PjxXLRoEd944w2WL1+e5cuX58cff8xr166V1m2SfPh7PXz4MIcNG6Y+a1erVi1OnDiRZ8+eLdW6C5Kbm8vz589zzpw5bNq0KStXrsxKlSo983jppZeKtTiLlZUV69evL0nec5KOiDB1EqPC1EmMCmOThO4JN27cIADOmTNH72U/S2ZmJt3d3enp6fnMkR59y5/UjRkzxmALWBw6dIj16tWjoigcO3aszouQPA+tVsu1a9eybt26BMDWrVvz4MGDOl+7adMmBgcHEwBr1qzJOXPmPJY8/Prrr6xfvz4TEhIYFRXF8ePHc9asWZw4cSJnzZrFPn360MzMjObm5uzbty/37t1b6j/v+/fvc9WqVezYsSPNzMzUhVRmzpxZ6ttyPK+0tDTu3LmT7777Ll1dXWlpaSlJXimSjogwdRKjwtRJjApj02tCB2ApgFsATud7zRnADgAXHv3r9IxrBz465wKAgbrUVxoJnVarZdWqVTlgwAC9l12Ybdu2EQC//PJLg9Wp1Wo5atQoAuCwYcMMNj0vPT2dw4YNIwB6e3vzyJEjBqk3Ozubc+fOZeXKlQmA3bt35/Hjx3W6VqvV8vfff2fLli0JgBUrVuRXX33FO3fukCQHDBjAqVOnPnXd7t27mZuby/j4eH700UesUKECAdDf359LliwxyPNhiYmJnDlzJoOCggg83Fuubdu2XLJkSak/b6cv2dnZPHr0KCdOnEhvb29aW1uXOMlLTk7mrVu3nnncuXPnb7VCp3REhKmTGBWmTmJUGJu+E7pWABo9kdBNBzDh0dcTAEwr4DpnAJcf/ev06OsCE7/8R2kkdOTDhVE8PT1LpezC9OvXj1ZWVgadHqfVajl+/Hh1ZcjS3qcuv/DwcFavXp3m5ub89NNPDTY6mZaWxq+//pqOjo4EwFdffbVYq0Tu27ePXbp0IQDa29vzgw8+4P79++nr68vbt2+TpJoQPPkhn56ezoULF9LHx4cA6OzszH/84x88f/683u6vMOfOnePnn39Od3d3Ndnp2bMn16xZY7DnRvUlJyeHx48f52effcYGDRoUK8nT5ahUqRK7dOnCxYsX8+rVqy9skicdEWHqJEaFqZMYFcam9ymXAOo8kdDFAaj26OtqAOIKuKY/gIX5vl8IoH9RdZVWQvf1118TgEEW78jv5s2bdHJyYosWLQw2WkY+TD7y7rlr164GnaZ29+5dDho0iABYv359RkdHG7Tuzz77jA55xnLXAAAgAElEQVQODlQUha+99hpjY2N1vv7kyZMcOHAgLS0taWZmxh49ejw1lfNZH/JarZa7d+/mq6++SgsLCwJgWFgYf/nll6ee0ysNWq2Whw4d4ujRo1m1alUCoJ2dHfv168d169aV2amKGo2GJ06c4GeffUYfH59SSfI6derE+fPnvzB77UlHRJg6iVFh6iRGhbHpmtApD88tmqIodQBsJunz6Pu/SFbI9/6fJJ2euGYsABuSXz/6fhKATJLfFlD+uwDeBYBKlSoFrl27Vqd2FUdMTAzGjh2L6dOnIzg4WO/lFyY8PBzTpk3DqFGj0Lt3b4PWvWnTJsyePRve3t6YMmUKHB0dDVb34cOHMWPGDKSkpKBXr14YMmQI7OzsDFJ3amoq1q5di/Xr1yMrKwutW7fGgAED4ObmptP1KSkp+O2337Bp0yZkZGTA29sbr7zyClq1aoUHDx7AwcGh0Ovv3LmDbdu2YfPmzUhOToaTkxM6dOiALl26oGbNmvq4xULl5ubi5MmTiIyMxJ49e5CamgobGxs0adIEoaGhCAkJga2tbam3o7Tk5ubi2rVr2Lt3L6KionD79u0iz8/IyICun3nlypWDl5cXmjdvjsaNG6NKlSpQFEUfTTeI9PT0ImNUCGOSGBWmTmJUGFubNm1iSAYVeaIuWR8LHqH764n3/yzgmnEAPsv3/SQA/yiqrtIaoUtNTaWiKJw8eXKplF8YrVbLTp060c7OjhcuXDB4/evWraO1tTU9PT156dIlg9Z97949jhw5koqisGbNmly3bp1BR0BSUlI4ceJElitXTn3Grjgjhvfu3eP333+vLr5Ss2ZNDhkyhMnJyTpdr9FouGXLFvbs2VNd7TE0NJTLly832IhZTk4OIyIiOGzYMPVZQzs7O77yyiv873//a/BRa2O6e/cuN27cyD59+rBq1arq4jK6HC+99BI7duxYJkby5C/LwtRJjApTJzEqjA0y5bJgfn5+bN++famVX5gbN26wQoUKbNq0aaEbW5eWqKgoOjs7s1KlSjxw4IDB6z948KC6XUCnTp148eJFg9Z/584dfvHFF3RyclJXxQwPD9e5U56bm8v//e9/6t57VlZWfOONNxgdHa1zGYmJiZwyZQrd3NwIgOXKlePQoUO5f/9+gyUHGo2Gu3bt4vDhw1mtWjUCoKWlJdu3b8/58+czMTHRIO0wJRkZGdy9ezeHDh3KWrVqFWubhfxJ3pUrV5idnc2cnJxnHoaadi0dEWHqJEaFqZMYFcZmiITu33h8UZTpBVzjDCAeDxdEcXr0tXNRdZVmQjd8+HA6ODgYJaEiyZUrVxIGXvUyv7i4OLq7u9Pa2po///yzwevPycnhrFmz6ODgQGtra3722WcGf64rLS2NM2fOZI0aNQiAfn5+XL58ebGecVu2bBnff/99ddQvICCACxYs4L1793S6Pjc3l5GRkRw4cCDt7OwIPNzs/Msvv+Tly5dLemvFlpuby+joaH788cesV6+emqSEhIRwypQpPHXqlEmPQpWm+/fvMzIyskRJXlGHpaUlvb29+dFHHzEqKqpUtvmQjogwdRKjwtRJjApj02tCB2AVgJsAcgDcADAEwEsAduLhdgQ78xI1AEEAFue7djCAi4+Ot3WprzQTuryEKiYmptTqKMobb7xBc3Nz7t+/3yj1p6SksFWrVgTAiRMnGnShljw3btzg66+/TuDhJtlr1qwxeOKQlZXF//znP2zQoAEBsEaNGpw2bZpO0w/zPuTv3bvHefPmqSOPDg4OHDZsGI8ePapzO+7du8elS5eyTZs2aoe/ZcuW/PHHH9XtEwxBq9Xy9OnT/PrrrxkSEqK2xdXVlaNGjeLvv/9usBVLTZUhkrzRo0czMjLyube+kI6IMHUSo8LUSYwKY9P7CJ0hj9JM6K5fv04AnDlzZqnVUZTU1FS6urrSxcWFd+/eNUobsrKyOHToUOLRCpjG2rds79699Pf3JwC2aNGChw8fNngbtFott23bpk6ltLe354gRIwrdZuLJD3mtVsvo6GgOGjSINjY26qjd3Llzi/U7vnr1KqdOnUovLy+1k9+tWzeuWrXK4NsPJCQkcMGCBezatat6T/b29uzevTsXLlzI69evG7Q9piozM5N79uzhO++8QxcXF1aoUKHQI29EVpfDwsKCXl5eJUrypCMiTJ3EqDB1EqPC2CShK4Sbmxu7d+9eqnUU5dChQ2pn3RgjZOTDJGTu3Lm0sLBgvXr1ePLkSaO0Q6PR8Mcff1QX6+jfv79Bpx3md+zYMQ4aNIhWVlYEwPbt23Pjxo3UaDSPnVfYh/zdu3c5d+5cBgQEEACtra3Zv39//v7770+V8yxarZZHjhzhP/7xD3VqqL29Pfv3788NGzYYfKQsIyODW7Zs4YgRI1i7dm014fDz8+OECRO4e/dug2zL8KLQarW8cuUKf/jhBzZv3lyduqtrkpe3oXphSZ50RISpkxgVpk5iVBibJHSFGDJkCCtUqKBz57q0zJkzhwA4depUo7Zj7969rFatGm1tbbls2TKjtSM1NZUTJ06kra0tLS0tOXr0aKakpBilLcnJyfz666/VZKp27dr85ptveOvWLZK6f8gfPXqUI0eOVBdiqVGjBsePH88zZ87o3Ja85+3effddOjs7EwDLly/Pt956i//73/8MntxptVrGxsZy+vTpDA0NVffbc3BwYLdu3Th37lyjrOT6IkhKSuLy5csZFhZGR0dHKopS4iRPOiLC1EmMClMnMSqMTRK6Qvz8888EwD/++KNU6ymKVqvla6+9RjMzM/7+++9GbcvNmzfVZ7gGDRrEtLQ0o7Xlxo0bHDp0KM3MzFiuXDlOnjyZqampRmlLdnY2f/31V/VnY2Vlxddff52zZs0q1jN/Dx484C+//MIuXbqoz10FBgZy5syZxVpVMjs7m+Hh4Xz77bdZoUIF4tFKmf379+cvv/xilN9bamoqN2zYwOHDh9PV1fWxZ+/ee+89/vrrr0abWvwiSE5O5vLly9m2bdtiJXm6HBUrVmTXrl25ePFiXr9+/W+7AI4wDuksC1MnMSqMTRK6QiQlJZnEyBhJpqen08/PjxUqVOD58+eN2pacnBxOmjSJiqLQ09OTx44dM2p7YmNj2bt3bwIPl4afPn26wZ8jy+/MmTMcNWoUHR0dCYAeHh6cPn06k5KSilXOzZs3OWPGDDZq1IgAaGZmxrCwMP7nP/8p1rOMWVlZ3Lp1K4cOHcqKFSsSAG1sbNijRw/+9NNPBl1QJY9Wq+X58+f5/fffs3v37upUQkVRGBQUxPHjx3P79u0GX9n0RZOcnMxly5axTZs2LF++vF6TvJdeeomdO3fmokWLeO3aNUnyRKmRzrIwdRKjwth0TeiUh+eaFk9PT8bFxZVqHf7+/qhQoQIiIyNLtR5dXLlyBUFBQahYsSKio6Ph5ORk1Pbs3r0bb775Jm7fvo0pU6bgo48+gpmZmdHaExMTg88++wzh4eGoXLkyxo4dixEjRsDe3t4o7bl//z6+/PJL7Nu3D/v374eFhQW6dOmCwYMHo1OnTrC0tNS5rLNnz2LlypX4+eefER8fD2tra3Tu3BmvvfYaunbtCjs7O53K0Wg02LdvH9avX4/169cjISEB5ubmaN26Nbp3745u3brBzc2tpLdcYjk5OTh06BB27tyJnTt34uDBg8jJyYGlpSWaNGmCNm3a4OWXX0aTJk1gbW1t8Pa9SG7duoXw8HAsX74cMTExePDgQaH/3WZnZ0Oj0ehUtrOzMxo3boxevXqhU6dOqFGjBhRF0VfTxd9UZGQkQkNDjd0MIZ5JYlQYm6IoMSSDijxRl6zP0Edpj9CR5Pjx42lhYWG0qXxPioqKoqWlJV9++WWTWFwiJSWFvXr1IgC2adOGV69eNXaTuG/fPoaFhalTxb755hujrc6Z91e7s2fP8uOPP2aVKlUIgFWqVOE//vGPYi8wo9VqeeDAAY4ePVrd7NvOzo6vvvoqV69eXayplLm5uTx06BAnTJigbskAgA0aNOCECRO4b98+oz0/mpaWxvDwcH788ccMCgqimZmZOrLYpk0bTp48mbt27ZIRPD3Q5S/LWVlZPHDgAEeOHEk3NzdaWlqWaCRPpmuKkpDRD2HqJEaFsUGmXBYuMjKSALh+/fpSr0tXP/30EwHwzTffNInOkVar5ZIlS+jg4MBy5cpx0aJFJtGu/fv3s0OHDgRAR0dHTpw4kcnJyQZtw5Mf8tnZ2dy0aRN79uypLhLSsGFDfvvtt0xISChW2RqNhrt37+bw4cPVlT/zplIuX75cp33y8rt48SJnzZrFNm3aqG1zdnbmG2+8wZUrVxplamaeP//8kxs3buSHH35If39/deqgpaUlmzVrxgkTJnDr1q0m84eXsqSkHZGcnBwePXqUY8eOpYeHhyR5otRIZ1mYOolRYWyS0BUhOzubjo6OfPvtt0u9ruL4+uuvCYDjxo0zdlNUly9fVhcFCQsL46VLl4zdJJLkkSNH2KdPHyqKQltbW44cOdJgbSvsQ/7WrVv87rvv1M25zczM2LZtWy5evLjYi4NoNBru2bOHH3zwgbripqWlJTt06MAFCxYUa0EV8mECtWbNGg4YMEB97s7MzIxNmjThF198wejoaKOu/vrnn39y8+bN/Pjjj9mkSRM1ATUzM6O/vz9HjhzJVatW8dq1a0ZrY1mhz45I/iSvbt26xUryFEUp9LCysqKfnx8//vhjRkdHm8QMBWEY0lkWpk5iVBibJHQ6eO2111ipUiWjb1+Qn1ar5YgRIwiA//rXv4zdHFVubi5/+OEHOjg40M7OjjNmzGBOTo6xm0Xy4bTHt99+m5aWljQzM2OfPn144MCBUh0d0PVD/ty5c5w0aRLr1q2rJmPdu3fnqlWrir0iZW5uLqOjozl27Fi6u7urHeaQkBB+/fXXPHnyZLHuWaPRMDo6mp9//jkbN26sjo45OTmxb9++XLJkidE3D09PT2dERAQ///xztm3blvb29up916pVi6+99hq///57Hj161GTi0VSUdkckJyeHMTEx/Oijj4qd5BV1SJL39yCdZWHqJEaFsUlCp4NVq1YRAPft22eQ+nSVm5vL1157jQA4b948YzfnMdeuXWOXLl0IgP7+/jx48KCxm6RKSEjghAkT1OX8g4ODuWLFilLpDBb3Q16r1fLw4cMcM2YMq1evTgC0tbVl7969uWrVKt67d6/Y5Z06dYpTpkxRRwIB0MXFhSNGjODWrVuZmZlZrDJv377N1atXc9CgQepzfADo5eXF999/nxs2bDDaM4t5cnJyeOTIEc6ZM4d9+/ZVRy2Bhxuvt27dmuPHj+dvv/1W7NHLF40xOiJ5v58PP/yQ7u7utLKyoqWl5TOPvBFYXZM8Hx8fjhs3TpK8F4R0loWpkxgVxiYJnQ5SU1NpZWXFMWPGGKS+4sjKymLXrl0JgIsWLTJ2cx6j1Wq5du1aVq9enYqi8L333jPqc1hPSktL4w8//EBPT08CYNWqVTl58mS9dvCf50Neo9EwMjKSI0eOZNWqVdVn5Hr27MkVK1YU+xk58mEyu2jRIvbo0YN2dnbqoirdunXjvHnzij0VVavV8sSJE/z222/ZsWNHtUwzMzM2btyYEydOZEREBO/fv1/stuqTVqvllStX+PPPP3PUqFEMDg5+LElwcXFh3759OWPGDO7bt+9vtdhKWeqIJCUlccmSJQwNDS3WNgxWVlb09fWVJK+MKksxKv6eJEaFsUlCp6OuXbuyVq1azM3NNVidusrMzGTHjh0JgAsWLDB2c56SmprKDz/8kObm5nR2dua8efNMatpbbm4ut23bpv4MLSws2LdvX+7evfu5p2Pq60Neo9EwKiqKo0aNUkfuLCws2L59e86fP7/YC6qQD+Nm27ZtHDlyJOvUqaN2fj08PDhq1Chu2bKl2InNgwcPGBkZyUmTJrFp06bq5uhWVlZs3bo1v/jiC+7atcvoCR5J3r9/n/v37+fMmTPZt29furi4qD8Dc3Nz+vn5cciQIVywYAFjYmKYnZ1t7CaXirLeEbl16xZ//vlntmvXjhUqVChWkifTNcuGsh6j4sUnMSqMrdQTOgCeAI7nO+4B+PCJc0IBpOY753NdyjZkQrd8+XIC4P79+w1WZ3FkZmaqUxxnz55t7OYU6OTJkwwNDSUA+vj4MDw83NhNesr58+f50UcfqdMxS7opeJ7S+JDPe0Zu3Lhx6jN3ABgUFMR//vOfPHr0aLETUa1Wy7i4OM6ZM4edOnWira0tAdDa2pphYWGcNm0aY2Jiiv0HjXv37nHLli0cO3YsGzVqpG4/YGVlxZYtW/LTTz/l77//XuyppKUlMTGRGzdu5GeffcYOHTrQ2dlZ/flaW1szJCSEw4cP55IlS3js2LEXIgl4ETsiKSkpXLFiBcPCwujo6KhzkmdhYUFra+tCjypVqqij5CX9XBDF8yLGqHixSIwKYzPoCB0AcwBJAGo/8XoogM3FLc+QCV1qaiptbGw4cuRIg9VZXFlZWezduzcB8IsvvjDJpcC1Wi1//fVXurm5EQDbt2/PEydOGLtZT8nIyOCyZcvYvHlztaP3yiuvcNu2bcVaHKe0P+S1Wi1Pnz7NKVOmsEmTJmrHtUaNGnzvvfe4efPmEk0fzMzM5Pbt2zlmzBj6+PioHd6XXnqJffr04YIFC3jhwoVix9hff/3FzZs3c+zYsQwJCVFH8MzNzRkcHMyPPvqI69atM5mOslar5aVLl7hmzRqOHTuWoaGhLFeu3GOjPIGBgXznnXc4b948Hjx40CRGH4vj79IRSUpK4tKlS9XpmvpamAXAY0meobdG+Tv4u8SoKLskRoWx6ZrQKQ/PfT6KorQH8AXJ5k+8HgpgLMmuxSnP09OTcXFxz90uXfXr1w+7du1CQkICrKysDFZvcWg0GrzzzjtYtmwZhg8fju+//x7m5ubGbtZTsrKyMH/+fHz55Zf466+/8MYbb+Cf//wn3NzcjN20p5w9exaLFy/G8uXLcfv2bdSoUQNvvfUWBg0aBA8Pj0KvjYyMRGhoqGEaCuDWrVvYunUr/ve//2H79u1IT0+HjY0NQkND0bFjR3Ts2BEeHh5QFKVY5d68eRM7d+7Ezp07ERERgRs3bgAAatWqhTZt2qhH7dq1i1Vueno6oqOjERUVhaioKBw+fBgPHjwAALi7u6NZs2bq0aBBA5OIZa1Wi0uXLuHo0aOPHXfv3gUAmJmZwcvLC/7+/mjYsCEaNmwIf39/VKlSxcgtL5ihY9SUJCUlYcuWLVixYgWOHj2K9PT0Iq/R9f+FlStXRrNmzdC7d2+0b9/eZH//ZcHfOUZF2SAxKoxNUZQYkkFFnqenhG4pgKMk5z7xeiiAdQBuAEjEw+Qu9hllvAvgXQCoVKlS4Nq1a5+7XbqKjo7GxIkT8dVXX6FFixYGq7e4tFotFi1ahNWrV6NZs2b47LPPYGtra+xmFSgtLQ2rVq3CunXrkJubi65du+LNN99ExYoVjd20p2RnZ+PAgQMIDw/HH3/8Aa1Wi/r166Ndu3YIDQ2Fk5PTU9ekp6fDwcHBCK192N4TJ07g0KFDOHz4MK5fvw4AqFq1KkJCQhAcHAx/f/9it48krl+/jmPHjuHYsWM4fvw4UlNT1bLzkhg/Pz9Ur169WMljdnY2Lly4gNOnTyM2NhanT5/Gn3/+CQCwt7dH/fr14e3tjfr166N+/fpwdHQsVttLC0kkJyfjwoULuHDhAi5evIhLly7h1q1b6jlOTk5wd3dH3bp14ebmBjc3N7i4uMDS0tKILTdujJZF2dnZOHfuHLZv344//vgDt2/fhlar1enaChUqwMfHB61bt0ZgYGCBnxniaRKjwtRJjApja9OmjWESOkVRrPAwWWtAMvmJ98oD0JJMVxSlM4A5JOsVVaahR+g0Gg1q1aqF4OBgbNq0yWD1ltTcuXPxwQcfICAgAJs2bUKNGjWM3aRnSkxMxFdffYXFixfD3Nwc7777LsaPH2+ybU5MTMTKlSuxYsUKnDx5Eubm5mjXrh3eeOMN9OjRA+XKlQNgWn+1i4+PR3h4OMLDw7Fz505kZGTA3NwcTZo0Qbt27RAWFoaQkJBiJxharRanT59GZGQk9uzZg6ioKNy+fRsAUK1aNbRq1QotW7ZEy5Ytiz3KRhLx8fHYv38/Dhw4gAMHDuD06dNqB9rd3R2NGzdWD39/f1hbWxer/aXp7t27OHnyJE6cOIHjx4/jxIkTiI2NRXZ2NgDAwsICnp6e8PX1hY+PD3x9feHr64vatWvDzMzMIG00pRgtq7KysnDkyBGsWLEC4eHhSEhIgEaj0enaqlWrqiN57dq1Q+XKlUu5tWWPxKgwdRKjwtgMNkKnKEoPACNJttfh3CsAgkjeLuw8Qyd0APDJJ59g+vTpuHbtmskmG/lt3rwZ/fv3R/ny5bFhwwYEBwcbu0mFio+Px9SpU7Fs2TKYm5tjyJAhGDduHOrUqWPspj3TqVOnsGrVKqxcuRJXr16Fra0tunbtin79+sHBwQEdOnQwdhOfkp2djejoaGzfvh3bt29HTEwMSMLBwQGtWrVSp1D6+/sXe5qjVqvFuXPnEBUVhb179yIqKkqdouno6IimTZuiWbNmaNq0KRo3bqwmv7pKT09HTEwMDh06pB4JCQkAAEtLSzRs2BCBgYEICgpCUFAQGjRoYPRRsPxycnIQFxeHU6dO4dSpUzh9+jROnTqFK1euqOc4ODjA29sbDRo0UP9t0KABatWqVezpskWRjkjpyEvyli9friZ5ubm5Ol1boUKFIv+7q1OnDrp27YoePXrA19cXFhYW+mi2SZIYFaZOYlQYmyETutUAfif5nwLeqwogmSQVRQkB8CseLpxSaKXGSOguX76MunXr4vPPP8fkyZMNWndJnTx5Et27d0dSUhLmzZuHwYMHG7tJRYqPj8c333yDZcuWQavVon///pgwYQIaNGhg7KY9E0lER0fj559/xq+//opbt27BxsYG3bt3R58+fdCpUyeTnZJx9+5dREZGYufOndi1axfOnTsH4GHHMi/BCw0Nha+vb7ETPJK4evUqoqKi1JG22NhYkISZmRl8fX3V5+SaNGkCd3f3YictCQkJOHz4MA4dOoQ//vgDMTEx6jRQa2trNGzYUE3wAgMDUb9+fZNK8gDg3r176jTTvETvzJkzSE7+/wkNDg4OqF+/vpro5U0/rV27domfL5SOiOE8ePAAhw4dwk8//YTt27cjKSlJ5ySvKHZ2dqhfv/4LmeRJjApTJzEqjM0gCZ2iKHYArgNwI5n66LVhAEBygaIo7wMYDkADIBPARyQPFFWuMRI6AOjcuTOOHTuGq1evmuziKE+6ffs2+vfvj4iICAwZMgTff/+9yT5Xl19CQgJmzpyJhQsXIiMjA126dMHYsWPRunVrvY9U6JNGo8GePXvw3Xff4eDBg2py165dO/Tq1QvdunUzyecE89y8eRO7d+/G7t27sWvXLly+fBnAwxG25s2bq1Mog4KCSjTF8a+//sKhQ4cQHR2NAwcO4ODBg0hLSwMAvPTSSwgJCVGP4OBgVKpUqVjla7VaXL58GUeOHFGPo0ePqnVYW1vDx8cHAQEB6uHn5wd7e/ti30tpu3PnDs6cOYMzZ84gNjZW/TcpKUk9x8bGBh4eHvDy8nrs8PDwKPKepCNiXJmZmTh48CCWLVuGiIgIJCUl6fxMXlFelCRPYlSYOolRYWwGXRRF34yV0IWHh6NTp05YsWIF3nzzTYPXX1IajQaff/45vvnmG/j6+uKXX36Bp6ensZulkzt37mDevHn4/vvvkZKSgoCAAIwePRr9+vWDjY2NsZv3TJGRkWjZsiX27t2L3377Db/99huuX78OMzMztGzZEj179kS3bt3g7u5u7KYW6vr169izZ486hTJvBM/GxgYhISFo2bIlWrRogSZNmqBChQrFLj83NxexsbE4ePAgDh8+jMOHDyM2Nlbt2Lq6uqJx48YICQlBYGAgAgICij1VU6vV4vz584iJicHx48fVRV3yVqdUFAUeHh4ICAhQF3Xx8/NDjRo1TPKPB3fv3sXZs2dx9uxZxMXF4dy5czh37hwuX778WEJQs2ZNeHp6wtPTEx4eHuq/eaN60hEpe0jiwoULWL9+PdasWYO4uDhkZmbqdK2dnR18fHzQpUsXdO/eHT4+Piaf5EmMClMnMSqMTRK6EiAJHx8fWFpa4tixYybZ2SvMtm3bMGDAAGRmZmL27NkYOnRombmHzMxM/Pe//8WcOXMQGxuLihUr4p133sGwYcPg4uJi7OY95ckPeZI4duyYmtzFxj5czNXLywtdunRB165d0bx5c5ObDviklJQU7Nu3D3v37sXevXtx7NgxdeqYt7e3+oxc06ZN4enpWaIFPvKelctL8A4fPoxr164B+P/kKzAwUD0CAgJQvnz5YtVBEjdu3Hhsxc680fc8Tk5OanKXdzRo0MAkR/OAh9P6Ll68qCZ4cXFxOH/+POLi4tRpqABgZWUFd3d3ODk5oUmTJqhXrx7q1q2LevXqoVatWgZblEXoB0lcvHhRTfLOnTv3wiR50lkWpk5iVBibJHQl9J///AeDBw/Gtm3b0LFjR6O04XkkJiZi4MCBiIiIQI8ePbBw4cIytU8SSezatQtz585VVxzt3r073nvvPbRv395kOqNFfchfvnwZW7ZswebNmxEZGYns7Gw4OjqiQ4cO6r5x1apVM1yDSyg9PR2HDx9+bApl3shXXsKQN8IWHBxc4ummycnJiImJeezIW3BFURTUq1cPAQEB8Pf3V4+qVasWu56//voLp06dwsmTJ3Hy5En164yMDLUuNzc3+Pj4qAuW+Pj4wNPT06RW2cyPJG7duoULFy7g/PnzapJ34sQJ3Lx5U93/D3g4JdXNzU1N8tzc3ODu7g53d3fUrl27zEw1/+Ze+m4AACAASURBVLsjiUuXLmHdunVYvXo1zp0799jvuTD29vaoUqVKoX/ss7CwQOPGjdG7d2+0atWq1LZhkM6yMHUSo8LYJKEroezsbLVzs3fv3jIzwpWfVqvFrFmz8Omnn6JcuXKYN28eXn31VWM3q9iuXr2K+fPnY+nSpUhJSUGdOnUwdOhQDBo0yOgrkRbnQz49PR0RERHYvHkztm7dips3bwIA/Pz80L59e7Rv3x4tWrQoE88+ksT58+dx4MABREdHIzo6Wl0IBQDc3NwQHBysJniNGjUq8YjXrVu3Hkvwnhxhq1KlymMJXsOGDVGvXr1ij0BotVpcuXIFJ06cUBcuiY2Nxfnz59Ul6s3NzVG3bt2nFi3x9PQ02d9bZGQkWrVqhcTERHUPvSf308s/0mNmZgYXF5fHkjw3Nze4urrC1dUVzs7OZfLz8O9Cq9UiLi4Oq1evxrp163Dx4kVkZWXppWxFUVC7dm20bt0avXr10luSJ51lYeokRoWxSUL3HH744Qe8//772LFjB8LCwozWjud15swZDBw4EEeOHEGPHj0wd+5c1KxZ09jNKrasrCz89ttvWLhwISIjI2FmZoZ27dph0KBB6NGjh1E61CX9kCeJkydPYtu2bdixYwf27duH7OxsWFtbo0WLFggLC0Pbtm3RqFGjEq9uaGhpaWmIiYnBH3/88dQUSjMzM3h7eyMwMBCNGjVCYGAgGjZsWOJVQf/888/H9n47fvw4YmNjkZOTA+DhCJS3t7e671veUa1atWInI9nZ2Th//jxiY2PVVSpjY2Nx6dIldRqqoiioU6eOuiG6l5eXmugZe3GcomKUJJKSknDp0iX1uHz5svp1SkrKY+eXK1dOTe7yH3Xq1EHt2rWLPS1WlL7c3FycOXMGa9aswfr163Hp0iV1r8TnpY8kTzrLwtRJjApjk4TuOWRlZcHDwwNVq1bFwYMHy/RfpTUaDWbOnInJkyfD3NwcU6ZMwYgRI0zqOYriuHjxIpYtW4bly5fj+vXrcHR0RL9+/fDWW2+hWbNmBvtd6etDPiMjA3v27EFERAQiIiJw6tQpAA+3FQgNDUWbNm3QunVr+Pr6msx0U10kJyerCV7eCFveMv2KosDT0xONGjVCo0aN1KmUzs7OJaorOzsbZ86cwYkTJ9Q94E6dOqWOhAKAs7Ozusl33lRKb2/vEiVdWVlZuHDhAs6ePYszZ87g7Nmz6jNt+ae9OTs7q4uW5D/c3d0NMn3zeWM0LS0N8fHxuHz5MuLj45867t+//9j5Tk5OqFOnjprg5f1bu3ZtuLi4yAifidBoNLh3716h52RmZmLn/7V33uFRlWn//55JJr33SgKkGUIKCaGFaugg6IUF17XsKtb11XVZG+7qq2u51vK6sj9X1tVVFxUFQRABIQYIRAgJBAg1CaTOpPc+5fn9kTxnz5w5M5lJm5nwfK7ruU6dk5NwOHO+577v752VhS+++AL5+flobW2FKc8KNMq7YMECrF27dlCRxx6WGdYOu0YZloYJumHyySef4Le//S22bduGO+64w6LnMhJcu3YNjz76KH766SckJSXhgw8+wNy5cy19WkNGq9UiOzsbn332GXbs2IGuri5ERkbi7rvvxvr165GQkDCqP3+0bvK1tbXIzs7GoUOH8PPPP+P69esA+gXe3LlzMX/+fMybNw8pKSk2J8qVSiUKCgpw+vRpflRWVvLbw8LCkJSUpDOioqKGHKlsbGzUEXg0lZK2OACAgIAAPoVyypQpfKRtsBojKTQaDcrLy3lnyitXrvDzwlYE9KE3Ojpab0RGRo5YHdtoPojQur2ysjKUl5ejrKxMb14s+FxdXTFhwgSdERERgfDwcISHhyMsLMxq6xQZ/femAwcO4IsvvsCpU6fQ1tY2bJHHHpYZ1g67RhmWhgm6YaLRaJCSkoKOjg5cvHjRqi30TYUQgh07duD3v/89Kisrcffdd+ONN96wShdJc2hvb8fOnTuxdetWHDp0CFqtFlOmTMEdd9yBdevWIT4+fsR/5ljd5Gnj7iNHjuDIkSMoKSkB0N+Ies6cOcjIyMDs2bORnp5utc3NjVFfX8+7UJ49exZnz57F5cuX+ZRG6tKXmJiIhIQEfgQEBAwp2kPdL2nfN2EPOGHUwsvLSyd9kvZ/mzRp0pCcSltbW1FcXMw7U9JatuLiYh2HSjs7O0RERCAqKoqvYxPWs5lTj2jJBxFCCBoaGlBeXo6Kigp+Skd5ebleSifQXxdJBR4VecIREhLCRJ8VUVNTg/379+Pzzz9HQUHBoJE/ikwmQ0REBFxcXODp6WlwP3d3dyxduhQrVqxAdHS0TWUpMMYHTNAxLA0TdCPAoUOHsHjxYrz66qvYtGmTpU9nxOjs7MSbb76Jt99+GwDw9NNP47nnnhsXNTC1tbXYvn07vv76axw/fhyEEMTFxWHdunVYt24dEhMTRyTty1I3eYVCgaNHj/KDtkews7NDUlISZs+ejdmzZ2POnDkIDw+3yRS3np4eXLx4EefOneNF3rlz59DY2Mjv4+fnpyPwaCrlUHrlAf0CpLq6mk+fFPaBE6Zu2tvbY9KkSYiJidEbISEhZv+9qfARGpYITUtaWlp09g8KCtIReNS0ZNKkSQgKCtJ54LX2B5Hu7m5UVFSgsrJSclRVVUkKBH9/f17ghYaG8lPh8PDwsMlrX4harUZ2djZiY2Nt5qUbIQSVlZXYvXs3vvrqK5w7dw4dHR0jdnwvLy8kJyfjlltuYSKPMSZY+32UMf5hgm6EWLduHX788UcUFRVh0qRJlj6dEaWiogIvvPACtm7dCl9fXzz77LN4/PHH4eLiYulTGxGUSiV27tyJ7du348iRI9BqtZg0aRLWrFmD1atXIyMjY8h94azlJt/c3IwTJ04gNzcXubm5OHnyJG/BHxISwrcUSE9PR1pams2KdkIIamtrUVRUxJuT0HnhA2NwcDCfNkmdKIeaQklpbW3F1atX+f5vNNJWXFys4xLp6urK93uLiorio2xRUVEICQkZ0oNnU1OTjmmJ0LikurpaJ+XN0dERkZGRvFGJVqvFokWL+Lq2oUY1LUlbWxsqKytRXV2NqqoqfiqcFwp9iqurK0JCQnRGcHAwP6Xz5jaxH0tKS0uxadMmXL58GTKZDO+8845V3HPMhRCC69ev8yKvqKhILx13OHh7e2PatGlYuXIlE3mMEcdavusZNy5M0I0QlZWViI+Px6xZs3DgwAGbeyAyhYKCAmzatAn79+9HYGAgXnjhBWzYsGFcpJlS6uvrsWvXLuzcuRNZWVno6+uDl5cXli9fjtWrV2PZsmVmObRZ601erVbj3LlzOH78OE6cOIG8vDw+TZPjOMTFxfECb/r06UhMTLTpFDZCCCoqKnD+/HneoIROhbVyXl5evAtlXFwcb1AyadKkIdesabVaVFdX873f6CgpKcH169d5500AcHJy4sWdOMIWERExpH+D3t5elJeX86YlYvOS5uZmnf2dnJx4kxKhYQk1LQkJCbG5ukygP9KnUCigUChQXV3ND6VSierqan69VJ82V1dXXuDRERQUpDMfFBQEPz+/MRcJ3d3dsLe3h1wux/bt27Ft2zZ8++23Y3oOo4VWq0VpaSmysrKQmppqcL+qqips3boVOTk5aGhogFarNen4VOStWrUKy5cvZyKPMWSs9buecePABN0IQtsY/POf/8SDDz5o6dMZNY4dO4ZNmzbhyJEjCAoKwjPPPINHHnnEJmuzjNHR0YGDBw9i9+7d2Lt3L+rr6yGTyTBr1iysWLECy5cvR1JSktEHAFu6yTc2NiI/P59vKZCXl4e6ujoAgFwuR0JCAu84OW3aNCQmJtp8lJYQAoVCwadOUpF35coVHYMSOzs7TJw4ETExMYiNjdUxKAkPDx/yQ6BarUZlZSVKS0tRUlKiM0pLS3XEBcdxCAsL00mhFLYECA4OHtJ57N27FxMmTOCNSsTmJeIaNplMhtDQUB3TElrHRmva/Pz8bPKlFiEEra2tUCqVUCgUOlPhfE1NjWSKoJ2dHQICAhAYGIigoCAEBgZKzgcEBMDX13fExINKpYJcLse2bduwa9cubNmyBe7u7iCEgOM41NbW4ttvv8Vnn32GWbNm4a677sLs2bNH5GePBebeR5uamnDkyBF8/vnnOH78OBoaGkwyZgH6XWeFIi8qKoqJPMag2NJ3PWN8wgTdCKLVapGZmYlTp07h7Nmz4y71UgghBIcPH8Zf/vIXZGVlwcfHB08++SQef/xxi/fVGg00Gg1OnjyJffv2Yd++fSgoKADQ7364ePFivi9ceHi4zuds+SZPo1qnTp1Cfn4+zpw5g4KCAj51TSaT4aabbkJKSgpSUlKQlJSExMRE+Pv7W/jMR4bW1lbenEQ8FaZQOjo6YvLkyToplDTCFh4ePuRollar5fu/0cgaHdevX4dCodDZ38HBARERETo934QjJCRE0gl0sGu0s7OTr2ETGpbQUVlZqdezzNHRka9fE5qWCGvZAgICbKaHohQdHR2oqamBUqlEbW0tampq+EGX6VQYhaXY2dnBz88PgYGBvAgMCAjQGf7+/oiIiEBQUJDB86BiTqPRYOnSpVizZg0eeeQRPk28ubkZL7zwApydnfHAAw9g+/btCAkJwcMPP8wLPmtnJO6jDQ0NOHz4MD7//HPk5uaiqanJZJHn6+uLlJQUBAQEGN0vMDAQq1atslnzKcbQseXvesb4YMwEHcdxZQDaAWgAqMU/lOv/VnkfwAoAXQDuJ4ScNnZMaxN0QH+9WWJiImJjY5GTkzNi1uLWzIkTJ/D6669jz549cHZ2xr333ounnnoKcXFxlj61UaOmpgY//fQTDhw4gIMHD/JRjJiYGNx88824+eabsXDhQpw7d25c3eSp+6OwpcDp06d1xEVwcDASExP5kZSUhNjY2HHzf0Gr1UKhUOiZkxQXF+tF1ezt7REREYFJkybppFDS6XBqs7q7u/lo2vXr1/WmDQ0NOvvb2dkhLCyMT52kQq+5uRmrVq1CeHj4kM5Hq9Wivr5ex6REbFyiUCigVqt1Pmdvb4/g4GAdsxJxPRutX7MF0WEIQghaWlp0RF9dXR3q6upQW1vLT2tra1FfX8/XtlIeeughbNmyxejPyMvLw+OPP44VK1bg+eefh5OTEy/W/va3v6Gurg4PP/yw3gsnW2G0HpYVCgV++uknbN261awWC6bg4OCAmJgYZGZmYvXq1UzkjXOYoGNYmrEWdGmEkAYD21cA+B36Bd0MAO8TQmYYO6Y1CjoA2LFjB9atW4ennnoK7733nqVPZ8y4cOEC/u///g9ffPEFent7sXTpUjz66KNYuXKlTdbcmIpWq0VRURGysrJw6NAhHD16FB0dHeA4DpMnT8aKFSswb948zJ07d9A3vLZKXV0dzp07x4+zZ8/i4sWLfORGLpcjLi6Ob9ZNx+TJk206SiOG1suVlJTwpiR0Wlpaqlev5uvrq5M2KZyPiIiAs7PzkM+ls7NTpx2AuDVAVVWVXq2Rp6enTjsAcVuA0NDQIT2UarVa1NXV8XVr1KxEOK9QKCTdKqlxibiGTTy8vb1tWvhROjs7UVdXh/r6etTX1yM4OBjTpk2T3Le+vh7btm3DgQMHcO+99+L2228HAF7M9fb24ne/+x26urqgVCrh5uaGt99+G9HR0WP5Kw2bsXpYpu6b+/btw5dffokzZ87o1NgOF0dHR0RHRyMzMxO33HIL0tPTzWoxwrBemKBjWBprEnQfAThMCPlqYPkKgAWEEKXU/oD1CjoAePLJJ/HBBx/gP//5D371q19Z+nTGlPr6evzjH//ARx99hOrqaoSHh2PDhg148MEHjaYOjRdUKhXy8vKQlZWFXbt24fLly3yKXlxcHObNm8cLPFuxGR8KKpUKV69e5QUedZ0sKyvj93F0dNQRetRtcvLkyUN2FrVmmpubeYFH0yhpvVpZWZle6mJgYKBONE08htp+Aeiv36uursbu3bvh6+ur1wqgsrJSsgecp6enwVYAISEhCA0Nhb+//5Dqjjo6OnRMSqhRibCOTalU6kWxgP6ICDUnoYPWrNFBUxvd3NwGFX89PT04dOgQ4uPjrTJ9vqenB7fccgsKCwvxz3/+E2vWrNHbp7i4GM8++yzCw8Px/vvv45133kFVVZXBF41ZWVm466674Ovryw8fHx9+KpwXrjPl7zkcLN0rsaSkBAcPHtTpBSlFUVERsrOzUVtba7Ixi5OTE6KiorBkyRKsXr0a06dPZyLPBmGCjmFpxlLQXQfQDIAA+IgQskW0/QcAbxJCjg0sZwF4lhCSb+iY1izoVCoVMjMzkZeXhyNHjiA9Pd3SpzTmqNVq7NmzBx9++CEOHjwIOzs7rF69Gg888ACWL18+Lh/YxRw+fBizZ89GQUEB3xPu2LFjfCQiNDQUs2bNwqxZszBz5kxMmzZtXLmGStHR0YFLly7hwoULOqOiooLfx97eHpMnT+bdJoXNu4cjYqwZrVYLpVKpkzpJI2t09Pb26nzG3d2dj6JRcxLx/GDXk7EHkZ6eHr1omrgtQE1Njd7Dq729PYKCgvTaAYiHv7//kCK0VPiJh7iWra6uTjKFztnZmTcpEdawCevY2tvb8cknn6CyspJvBzBv3jyzz3W00Gg02LNnDw4fPoyzZ8+ipqYGr732GqZOnQp/f394e3tDoVDgrbfewrp16zB37lx8++232L59O7744gvJFOhz587hww8/RGNjIxobG9HU1MRPpUQ0xd7eHt7e3vDx8eGndHh7exsdpkShbfFhuampCTk5Odi6dSuOHj2K+vp6k0Wes7MzYmJisHjxYqxevRppaWk2b0A13rHFa5QxvhhLQRdCCFFwHBcA4CCA3xFCjgq27wXwhkjQ/ZEQUiA6zgYAGwDA398/9ZtvvhnWeY0mLS0teOyxx9DT04MPPvgAoaGhlj4li1FZWYkff/wRBw4cQHNzM7y8vJCZmYnMzEzExMSMi1QpKTo6OvRS1DQaDa5du4bz58/zYqa2thZA/4NRdHQ03xstLi5uSI2obRFqviE24Kiurtapv/Ly8uJTAIWpgKGhocNKUbR2tFotWlpadGqwxEOc0gn0R9T8/f354efnpzN1dnaGv7//kK8xjUaDpqYmNDQ08IMKAjrf0NAgmbomk8l4ASCMCIlFgY+Pz5D+bTUaDVpbW9HU1MSP5uZmydHa2ir5wC2TyeDk5ASO4xATEwNPT094eXnBy8sLnp6ekMvlCAkJ4Zc9PT0tkkbc2NiIvr4+lJeXo6amBkuXLoWzszPeeOMNLF68GGlpaXjvvffg4+ODtWvXwtPT06zj9/X1ob29HW1tbWhra0N7eztaW1vR0dGB9vZ2fhudp8vGhCDQn47t7u4ONzc3uLm5wd3dHa6urvy8m5sb5HI5fH19+X3odrrNVmhpaUFhYSEOHTqEoqIis2r2HB0dER4ejqSkpEHFXWhoKFJSUob1/5phHlLf9QzGWLJw4cKxd7nkOO5lAB2EkLcF68ZVyiXlypUrmDNnDtzd3XHs2LEbWtQB/ZHLAwcO4N///jf27NmDvr4+REVF4e6778b69evHnZGKqW/tlEolTpw4gRMnTuCXX35Bfn4+n6bp7e2N1NRUpKWl8dOIiIgb5otarVbj+vXrOg27qeOksLUA0N8kXdjDTdjTzZz+gbZKT08PqqureUFM0yfpqKysNNhgW2hKIkyfFEbXhiOYe3p6eFdI4aDplLW1tVAqlairq5MUVi4uLjpRNXFkTRhdG0o9nUajQWNjI+rr63mBrFQq0dTUhKNHj6KsrAzBwcFoampCfX29pHimeHt76whn8fD19eWnvr6+8Pb2HjUR+OOPP+L555/nG8q/8cYbmDx58qj8LCmEotqQmG5paeGHcLm5uRkajcbo8Z2dnXmRTQW1qcPDwwOenp4WM2wihKC6uho//fQTvvzyS+Tn5w+a1mkOHMfxWSCrV6/GggULEBYWdsN8d4wlLELHsDRjEqHjOM4VgIwQ0j4wfxDA/xJC9gv2WQngCfzXFOVvhBCjeYq2IOgAID8/H4sWLUJoaCh+/vlnBAcHW/qUrILm5mZ89913+PLLL5GdnQ1CCJKTk3H77bfj1ltvxU033WTpUxw2Q73Jq1QqXLhwAXl5eSgoKEB+fj7OnTvHR6p8fX35fnC0ZUB0dPS4Mhgxhfb2dpSUlODq1au84yQ1IBGLPW9vb95dUtjDbeLEiYiIiBg3LpyDQRts07TJnJwcODk56dWsidM7gf7oqCFTEmHd2nAMSqiwEqZO0hYA4lFfXy8Z4bC3t4e/v79kGwAqtKZOnWpU2NB2AB0dHcjMzMQjjzyCe++9l68N/P777/Gf//wHmZmZCAwMxHvvvYeZM2fC2dlZJ1pJR319vV6NJIXjOHh5eenUrgkjloZq1zw8PEz+O5eVlfGpprYCIQSdnZ3Yt28fbrrpJrS2tuqIP+FobW3lh3BZ2GLEEE5OTrzAEw6pdR4eHnB3d9ebd3d3h7Oz87DFEiEEpaWl2L9/P7Zt24bCwkLJnodDhfaRnDNnDlauXIn58+czkTcCMEHHsDRjJegmAdg5sGgP4EtCyF84jnsEAAgh/xhoW7AZwDL0ty14wFj9HGA7gg4AcnJysHz5cl7U3eiROjFKpRLffPMNvv76a5w4cQIAEBsbi9tuuw233norUlNTbbK560je5Ht6enD+/Hle4J0+fRpFRUV8jytnZ2ckJCQgKSkJycnJSEpKwtSpU81OrRovdHZ24tq1a3yTbuo4ee3aNZSXl+v0BqNvsoUCT9jLLTw8HI6Ojhb8bUYPqWuUEIKmpiZe5BlqtK1UKiWFn4ODg44hibHI2nDFX0NDg47AE6ehCtcJH4xffvll/PnPfzZ6/OzsbDzzzDNYv349nnnmGchkMt5B8oEHHsCiRYuwZs0aeHh4ICMjA88//zxWrlwpeSwqTurr6/VSUsWjoaGBTxE15rJoZ2cnWbsmXGesds1WHuKHcx/t6+tDa2sr2traeJFnbJ6mktLUUrp9sEgh0P/vQVNFDQ3hdkPzdNBUX41Gg4sXL2L37t3Yv3+/ZF9DCiEEZWVlaGhoMLlmz87ODiEhIZgzZw5WrVqF+fPnIzQ01GauD2uACTqGpWGNxceQ48ePY/ny5fDx8cH+/fvHXXrhSFFdXY1du3Zh586dOHz4MDQaDYKDg7FixQqsWrUKmZmZNpOrPto3+d7eXly6dAlnz57lR2FhIZqamvh9JkyYoNcuID4+/oZ2UtNoNFAoFLh+/brkUCgUOg9DHMchKCiI7+EmZUBiq/Uqw7lGpXqsSUXWampqUF9fL/lQbG9vLxlJMxRdG46jYnd3Nx8t8/f3N9iXraamBl999RVycnLwm9/8BqtWrQLQX8cok8nQ2NiIe+65B6+//jpSUlIAAHPmzMFHH32EhIQEveN1dXWhvb0d3t7eZkeCVSoVL+6ERiXNzc16tYHC5ZaWFqPHdXBwgLe3N18PaGheWC8onB9L8yZLPywTQvh/Q2GtoLCeUFw/aGh0dHRIvgSRQigQpQSfIfHo6uoKFxcXtLS0ID8/H4cOHcLFixfNivTRvpVz587lW+/cKPXcQ8HS1yiDwQTdGFNQUIAVK1ZApVJh9+7dyMjIsPQpWTVNTU3Yu3cvfvjhB+zfvx9tbW1wcHDA3LlzsXTpUixevBiJiYlWG72zxE2e1mWcPXsW586d49sFXL58WedBIjIyUkfgxcXFjWsXSXNQqVSoqqrinSaFvdzKy8tRWVmp91Dm6OiIsLAwvf5tQlv/gIAAq7tWx+oa1Wq1aGxs1EmZFDfVFkbVDJlpODo66hm8COdpfZqwXs0c44yuri6sWrUKRUVF+Oijj3Drrbfq7XPs2DF8/PHHeOWVVxAREYFz587hpZdewvvvv4/IyEi9/WlvUgBwc3MzGFEzFlXz9PQ069rRaDRoa2vj69SkatiM1a6JG8GLcXR0NFi3Jl4vTF8UrqMRqMEYbw/LKpVKR+CJBZ+hbYb2M1UgCnn66adx+vRpFBYWmmzOYmdnhwkTJiAjIwPz588fNGshPj4eCQkJN0Q6+3i7Rhm2BxN0FuDatWtYvnw5rl+/js2bN2PDhg2WPiWbQKVS4dixY9i7dy/279+PCxcuAAD8/PywaNEiZGZm4uabb7aqnlHWdJNXq9W4du2aXruAy5cv66TwBAYG8uKOtgyIjY1FRETEDVejZwhCCBoaGnjjEeGg65RKpd5Dsb29PW82IjQhEQ9zaqOGizVdo0JoeqJQ5NXX1/MRNvG8sdRET09PHRMSoaumeHh5eeHUqVPIy8vD2bNnUVFRgddeew3x8fHw8fFBSEgITp8+jb///e/44x//iNjYWPz1r39FSUkJNm/eLCkeS0tLceDAAZ3omjCi1tzczLtUGoLjOHh6evLRM/FUOKgAFEbUzIls0ogUrUUT16oZql0T1q8N5m4J9DtcStWoiUdtbS2mTZtmsH7Nzc3thr43CQVie3s7urq6+NHd3a237Orqit/97nf8CwKapvnll19i165duHTpkkn/fqbi7e2N5ORk/iXseBR51nofZdw4MEFnIZqbm3H33Xdj//79eOihh/D++++Pa8v10aC6uhpZWVk4dOgQsrKyoFAoAAATJ07EggULMHfuXMydOxeTJ0+2WJqILdzkVSoVrl+/jitXruDy5cv89PLlyzqOiHK5HJMmTUJ0dDSio6MRFRWFqKgoREdHY8KECTf0A5UUGo0GdXV1Oj3chFNan0Z7EgpxcXFBSEgIgoKCeMMRqamfn9+w/+62cI2aQm9vL19/JmyfIDQlEfdXk/rbU+zt7eHj46PTioDjOKSmpiIgIADfffcdFi1ahDlz5uDll1/Gww8/jDvuuAPe3t5DstInhKC7u1svmkYjZoNF17q6uoweXyaT6UTPvLy8dCJmM4lDSQAAIABJREFUhhwghfNubm4mRwnVajVff2bK1NDo6ekx6ee5uLgMq25NmNro6upqU+0QRgNCCC5evIgvvvgCe/fuRUlJicn/Fqbg4+ODadOmYenSpbj55puRkJBg03/z8XIfZdguTNBZEI1Gg02bNuHNN9/ElClT8PXXX0vWXzAGhxCCy5cv8wIvJyeHryMLCgpCRkYGL/ASExPHTHzY+k2+oaEBV65cwZUrV/h2ASUlJSgpKdF5gJTL5YiMjMSkSZN4J0nhcHd3t+BvYd3QJtlU4AmHUqnkrf6lxIdMJoO/v7+O4YghExJ/f3/Y29vrHcPWr9HhQOvThIYkwjo0Ye0ajaI1NzcbjQYC/W0ghBE0U6JqdLi7uw8pLZeafwhFnjhqJo6yCSNrbW1tg5pocBzHR8WkHCCF0TPhOqETJB2Ojo4mvWjr6+vDvn37kJiYqFe/Zqj/naFhTtTJ0dGR73UnHLQHnrAXnnAdHYaWbTkypVarcfLkSXz++ecoKCgYdN+ysjKzeu1R9+YlS5YgMzMTU6ZMsRmRdyPfRxnWARN0VsCBAwdw3333oaWlBS+//DL+8Ic/SD54MUxHq9Xi8uXLyMnJ4UdFRQWA/hqW9PR0zJgxA+np6UhPT0dISMionMd4vckTQqBUKnlxR9sFUBdJcY8uX19fSfdIOn+jOnGaQ1dXFy/uhFOxlX9NTY3BmhpfX189o5HOzk5Mnz5dr/m4n58fuw8ZQKVS8aJJoVDoRdaEUTVhfVpzc7PRqCDQL5qE0TFxPZpwWSqK5unpOSRRSAhBR0eHnsgTR9HE80JzECqsTIE2FB8soubu7o6amhqkpKToNCAXRtTc3NxMcu3UarXo7Ozk68+kpoYGFYQdHR0oKyvjX2gZc5yUwt7enhd5Li4uOqJPvN7FxUVnXrhOvI0OJycnq6rTpemcX3/9NXbt2mW2OUtAQABSU1OxdOlSLFy4EPHx8VZ5Xxqv3/UM24EJOiuhrq4Ojz32GHbs2IFp06bhX//6F5KTky19WuOKiooK5OTk4JdffkFubi7Onz/P1ziFhoby4i4tLQ3Tpk2Dj4/PsH/mjXqTb25u5sXdtWvXePdIaioiTt3x8vLiRZ7YRTI8PBwhISFW+SVujRBC0NbWpmc8Irbwp0PoiCrG29vbYHNsoenIWDTIHi9QsxKpXmo0qmYoikaXTbHQF0fLDEXUDPVUo1NTI2kUKpqEgs+UyJohc5DBzFkoHMcZjKZJRdakpobElfCaJoTghx9+wLVr1/D//t//w+XLl9Hd3c0Lv87OTtTV1eHtt99GW1sbPD09sWDBAjg4OKCjowNdXV3o7Ow0OGitG102tfWAEGdnZzg7O/MiT2relKnwOMJlum6o4lGr1eLChQvYunUr9u3bh+LiYpP6BVICAgKQlpaGZcuWITk52eg52NnZITY2Ft7e3mafpzncqN/1DOuBCTorY8eOHXj88cdRX1+PDRs24NVXX4Wfn5+lT2tc0t3djcLCQuTl5SEvLw+nTp1CcXExvz0iIgLJyck6IyIiwqyHG3aT14cQgrq6OpSXl/MukmVlZSgrK+NNRVpbW3U+I5PJEBwcjPDwcISGhvLukcJpSEjImFqpjxeysrKQkJDAm4uIzUbosnCY0iBbbDxCXRylBq1RYwwONSwxpbeasdq09vZ2k1Lh7O3tjfZVE9efGZoKI2mmigBCCHp7e3HgwAEkJSXpRNHE8zR6Ziy6RvcxZj4jhZOTk170rK+vDwqFAgsWLNCJtNnZ2SE/v7+FblJSEkpKSqBWq/H73/9eL7pGhyHRTAhBX1+fnsgTG50Ih3i7lCkK3Ydu6+7uNjktUoyDg4Oe2KPDycnJ6HB0dOTn7ezscPHiRRw5cgTFxcUjWq8H9NdY3nTTTVi4cCGWLFmCtLS0ERV57LueYWmYoLNCmpub8corr2Dz5s3w8PDAn/70Jzz66KPjtrGxNdHU1ITTp0+joKAAhYWFKCwsxJUrV/gvO09PTyQnJyMxMRFTp07lLf8NWf2zm/zQaG1t1XOPrKys1DEWkUrb8fX11XGPDA4ORnBwMD9PjUZsuY5lpDH3GqWpeULTEXFjbHFdWmNjo9FUPKGD42ANscW90jw9PW2mzsaaoJE0Q1E0c/qqmRtJc3FxMak+jU6rq6uRkpIiGUkTrjP1pYBKpdKJigmja4NFz+i0uLgYjY2NCAoK0om8mRNpEv9Nhhs5k4qkOTk56a0TZzuIhSMVe0LBN9i6np4eyXU9PT3o7e3l5+kwVUS++uqrCA4Oxu7du5GXl4e6urohRS2lcHV1xZQpU3iX7NTU1CG37WHf9QxLwwSdFXPhwgU8/fTTOHjwIMLDw7Fp0ybcf//97GF0jOns7ERRUREv8M6cOYOioiKdAvuQkBDEx8fzAo9Oz549y27yo0RbW5uec2RVVRWUSiWqq6v5GjOpL38fHx8EBQUZHNRMZCRcJK2dsXoQoXVnxkxHDPVLG6xOyc3NTbL5tXjekKOjq6sra5g8THp7e/WiZ4PVpQmjZ+IoWmdnp9lRGgcHB8k6NHNr09LT0weN3rz44otQq9V46623oNFoYGdnxwuUbdu24bXXXoOTkxOCgoIwa9YszJs3j+8ZZyh6Jl5vKNJGBdNQsbe31xF7QtEnnqdDuGwo+iaMuEmtc3R0hFwuB8dxIIRApVKht7eXF3ziqUajwcyZM3Ve2BBCUFxcjJ07d+LHH3/ka+MNoVKpUFdXZ3Kto5ubGxISErBgwQKzRB4TdAxLwwSdDXDo0CG89NJLOHHiBCZMmICnnnoKDz74IHMOtCBarRbl5eW4ePEiLl68yPd0E/fv8fb2RkJCAmJiYhATE4PY2FjExMRg0qRJLOI6Bmg0GtTX1/OOkXRKzUPoUCqVkg9IMpkMfn5+vIkIdYykhiJ00GVPT0+bEwbW/iAi7Ikm1QBbqv5MXIM2WPRIJpPpGYxI1Z3RZXGtGZ03x9afMTgajQYdHR04ePAgkpKSDEbOxLVpUnVqwpREOjUU6cnNzcWsWbOMntu9996L9PR0PPHEE1Cr1ZDJZJDJZCgpKcGnn34KPz8/hIWF4bvvvsOTTz5p9HgFBQV6qYtUBBm6nmg6qrGombEIGl02dzrc6BjHcbzIk5oa22Zo38H2c3BwQEtLC37++Wf8/PPPKCwsRGNjo8lppp6enkhMTMTChQuxaNEipKSkwMPDQ2cfa7+PMsY/TNDZCIQQHDhwAG+++SaOHDkCT09PPPzww3jiiScQHh5u6dNjDKDValFRUcGLvMOHD6O9vR1XrlxBXV0dv59MJsPEiRP5fm6TJ0/mx8SJE1lPwjGGphFScSc2ExHO19bWGkwflMvl8PPz03OLNGQm4uPjAxcXF4uKwPH+ICJuki1lMkLTDo25OZoaLRI7NA6l/ozOu7q6MoGI0blGhWmG4tq0qVOnDvrCdN26ddi0aZOeedlLL70EFxcX/M///A9cXFywYcMGhIeH46WXXjJ4LAcHB4MRJCpMhhs9GyyiJjVPo2lCVCoVL+6kUimFKZVS0TepeWNROqmpuc6ig+Hs7IzFixcjLy8PtbW1Jgk9Ly8vJCYmYtGiRVi4cCGqqqqMinYHBweEhITY3As/hu3ABJ0NkpeXh3feeQfbt28HAKxatQqPPPIIlixZMu7Tw2wN4YNIS0sLiouLceXKFVy9epUfpaWlelbmwcHBfA83od1/ZGQkwsLCWNqthenp6dExDaHukcJ5YVNrY06Sjo6OOuYhQkMRoXEIrS+jwxSbdlMY74JupOjr69NzahTWmkm5OBrqkWZOKqG59WaDuTbSYUt1h9Z0jb733nv4+OOPcenSJSxYsACvvPIKAgIC4OzsjAkTJuCtt96CUqnEe++9B47jsHLlSqxduxYPPfSQwWPu3bvXaERNHF2jkTmxgBIKqeEYnQgxJZImHFLrhMPBwUFvHwcHB5OWHRwcIJPJoNVqBxWBQtGoUqnQ19fHT8Xz8+bNw8033wyg/95eWFiIrKwsHDx4EGfOnBm01YipyOVyxMTEYPbs2cjMzER6errZRmsMhiFGXdBxHBcO4HMAQQC0ALYQQt4X7bMAwPcArg+s+o4Q8r+DHftGFXSUsrIybNmyBf/6179QV1eH0NBQ3HPPPbjrrruQlJTEbhJWgCkPIoQQNDQ06PRxo6O0tBTV1dU6X8wcxyE0NJQXeUKL/7CwMISHh8PHx4f9+1sRarUazc3NqK+vR2NjIz+VMg8RrjOWKujg4KBjFkJFn7CZtaEG1h4eHnzkx5oelm8U1Gq1Tg2ZVE80c1wbOzo6zK6pksvlkr3QxDVnxnqeSQ2hRf5ItRqxpmtUpVKhrKwMVVVVOH/+PKZMmYLu7m7U19dj/fr10Gg0WL9+PTo7O+Hn54fu7m588MEHiIiIGNPzJIRArVbriD2pyJlUxMxQdM2QcBLvK143ktjb2+sIPDoVD0Pr5XK5wWVD8x0dHSgsLMTJkydx7tw5vgfhSCCXyxEbG4t58+YhMzMTaWlpCAsLY9/fDLMZC0EXDCCYEHKa4zh3AAUA1hJCLgr2WQDgD4SQVeYc+0YXdJTe3l788MMP+PTTT7F//35oNBrExsbizjvvxJ133on4+HhLn+INy0g8iPT19aGqqkrS5r+srAzV1dV6D/7Ozs46Ai80NJR3f6TzgYGBLKJrxRBC0NnZyRuHUMMQ4RCbh9DR2tpq9O28sHm1vb09QkNDDRqGCJtWi+vK2PVjPWg0Gj590BzXRqnaM6l1pvS+EyOXy81ybZTqdebs7IzS0lKkpaUZdW+kKYLWQFNTE4qLi1FWVoaZM2eOuZizJpqamnDx4kX87W9/Q1RUFJ544gl0dXVBrVbz4lGtVkOpVOKVV15BZ2cnnn32Wbi4uKC3txd9fX1GpzQFk0bb6DbhMt0unBd+ZqyIjo5GZWXloCLX0dERcXFxyMjIwKJFi5CWlobw8HAm8hhGGfOUS47jvgewmRByULBuAZigGxEaGhqwY8cObNu2DYcPHwYhBFOmTMGaNWuwcuVKzJgxgz2EjSFj8WZZo9GgtrZWx9pfbPdfU1Oj90Amk8kQFBTEW/oHBQXxU+F8YGAgq+mzMbRaLdra2vSMQqQaWZeWlkIul+vVlpnyAO/q6qrXqFqqSbWhITQTYY3jrRuVSqUj+MQ9z0x1bJSaCo09ent7h3yOdnZ2ZteU0e1SDo2m1JmJzTfYQ7c0GzZsQEpKCh599FFotVqd2lBCCN555x3s3LkTy5cvx/3334+wsLAxOS8ayRQLPeGUCkFDKZt9fX24cOECYmJioNVqdYZGo4FWq4W3tzfuu+8+AEB5eTlOnDiB/fv3Izc3F9euXRv0fuvs7IypU6di/vz5mDdvHlJTUxEcHDwWfyKGjTCmgo7juEgARwEkEELaBOsXANgBoAqAAv3i7sJgx2OCzjhKpRLbt2/Hzp07cfToUWg0Gnh7e2PJkiVYtmwZFi9ejNDQUEuf5rjGWlKFNBoN6urqoFAoUF1drTelbo91dXWSkR1PT0/e4VHo9ChcR81AvL29mZGDDSF1jQqbV4sbVIvNQsR9zMR1Zab2J3NycjJoDiI0CTFWSyaed3JyYg/YNoZWq9WrGTt69CimTp1q0KHRkCujuL5MyrhDuDwSmFNPNtgwVlsmTC0Upx+Kt1vDS9x169bh9ttvx5133skLOjp9++234e7ujtraWmi1Wjz66KMIDAy09CmbxXC+61UqFYqKinDixAkcOnQIeXl5eqUWUri5uSE+Ph4LFixARkYG0tPT4efnZ/QzMpmM3RPHKWMm6DiOcwNwBMBfCCHfibZ5ANASQjo4jlsB4H1CSLSB42wAsAEA/P39U7/55pthndeNQnt7O/Lz85GXl4e8vDzepCEsLAxJSUlISUlBUlLSoDcDhnl0dHTAzc3N0qdhMhqNRq9XGK3rohEdOgwVistkMj5VT9j7i9ZuSaXvWdrp8UZmNK9RqYbFUtEbU5sXd3d3m+VwJ5PJjEZqHB0ddSI0zs7Og0ZkxA/p7OXF6DMW91HaF00q+iKOyoi3iyM3UtEc4WfF82q1Wm95JJHJZLC3t4dcLodcLufnTZnSIbXezs4ODg4OWLhwIQICAoyewx//+EfccccdSEtL44UKx3G4evUqsrKy8Jvf/AbZ2dloaGjAbbfdBhcXlxH9G4w2I32Ndnd34+rVqzh37hzy8/NRUlIyIrV7MpkMoaGhiI+PR1JSEuLi4jBhwgSrEP2M4bFw4cLRF3Qcx8kB/ADgACHkXRP2LwOQRghpMLYfi9ANDUIIzp49i59//hnZ2dk4evQo/3A+ceJEzJkzBzNnzsT06dORlJTE+qUNA2uJ0I0GKpUKDQ0NvJ0/dXhsaGjQc4CkTo+G7iNyuVzH2dGQw6OU2QdL1RsetnaNqlQqg4YhUv3JDPUqk+pXNpQ6MScnJz0zEEN1YoPVjEkNYa2Yo6PjDfniw9au0eGi1Wp1asDokKofE9eOiadStWXiujPxoNvb29uhUCigVqv5Hot0m1B0bty4EcHBwVi1ahWio/Xfxff19WHRokV49913kZ6ezqchyuVyLFmyBBMnTsQ999yDzz77DN7e3ti4caNRgbh27VoUFxfzAlUsVAdbb2jZVIErta6wsBCzZs3SEcFSwpgK4aG8CFIoFMjPz0d2djays7Nx6dKlEakBtLOzQ3R0NDIyMjB//nxMnz4d0dHR7GWVjTEWpigcgM8ANBFCnjKwTxCAWkII4TguHcB2ABFkkB/KBN3IoNFocObMGRw9ehTHjx/H8ePHUVtbC6D/QTsxMRHTp0/H9OnTMW3aNMTHxzPbfBO50R5EjKHRaNDc3Kzj6NjQ0CDp8CiMDg7m3ufu7q7j7EijgsKpeJ0wQujk5DRGfwHrhF2j/dBoYnd3t2QjanHNmKFIo7F6MTo/nObMg9WIjWQfMqm0QUu8yWfXqGVQKBSoq6vD888/j3vuuQe/+tWvdNIle3p68OKLL6KmpgYymQw+Pj547rnnJEs5MjIycODAAbi6uuqs37x5M44dO4bOzk6cO3cO7e3tePXVV7FhwwaDJjcbN27E9evXdaKhVGSK1xlaT9dZqiUXx3F6Qk8o+AZbZ2dnB47jUFxcDIVCYdbPTkxMRGVlJZqbmw3u4+DggJtuugkZGRmYPXs2UlNTmcizcsZC0GUAyAFwHv1tCwDgBQATAIAQ8g+O454A8CgANYBuAL8nhOQOdmwm6EYHQgiqqqpw6tQp5OXl4dSpU8jPz+ejeHK5HHFxcUhISODHlClTMHHiRPafXQR7EBk+PT09fJqnOPVTal5oAtLW1jboF7aDg4NeGqjQ3ENs9iE08hCbfNhitJBdo2MLTe2TSieV6kNmqB+ZsVowYb8y4frhCEmKvb39oPVhUtuk6r2M1X8Jp+fPn8fs2bON2tKzlLHRY82aNXjwwQexevVqaLVacBwHjuPw/fffY/fu3XjllVcQFhaG5cuX49e//jXuvvtuEELAcRyee+45bNu2DeXl5YiNjcWf//xnuLu7g+M4LF68WEe0PfXUU5gxYwbWr18/Zr+bVqvVEXimTKXmCwsLERcXJ7mdLms0Gp11arWaXyfeLrWv1HrxstCIRWzMotVqERsbi4MHD8LOzg41NTU4efIkcnJykJWVhYsXLxqN+Dk6OvLN1KdPn47U1FTWR8+KMFXQDfkphRByDIDRf21CyGYAm4f6MxgjC8dxfF+z2267DUD/Ta+kpASFhYUoLCzE2bNnkZubi6+++or/nIuLC+Li4hAbG6szoqKi4O7ubqlfh2HjODk5ITg4eEiOXlqtFu3t7WhpaeGFntDNUWjsIVyuqKjQMfkwNa2FGnsITT3E5h7mNoYeqQbiDOuA4zhehHh6eo7pz6Z9yYT9yKT6kAlFoKE+Y1I9x+jo6OhAY2Oj5Daa0jcS4lIIrecyNsQ9yIbSl0yYqidcNjRvaH86bEGIdnZ28t/hHMfxYi03NxfJycn8Nh8fH96llN6zXn/9dWzcuBFtbW2oq6tDYGAgH+2Wy+W80LC3t4evry+cnZ35448FMpmMf7EwHDw8PGzuxVhQUBDWrFmDNWvWAOi/P1y4cAEFBQXIycnBsWPHcO3aNf7/am9vL06dOoVTp07xx3Bzc0NiYiLmzZuHtLQ0TJ8+nbVYsHJs77UzY0SRyWSIiYlBTEwM7rjjDn59W1sbLl68iKKiIhQVFeHSpUvIzc3F119/rRMZCQwMRFRUFCZPnoyJEyciMjKSb4wdFhZmNf2DGOMLoUHLcOjt7eWdHFtbW3UaQhsanZ2d/P7V1dU6TaHNrXsYbvNnqfouOt/d3Q2NRmMTD5aM4WFvb8+/ULA0tA+ZKb3GCgoKEBsba7A2zNSaMLrc1dWFlpYWk0xPRhuO4yTFnzm1X1LLptaDLV++HFFRUUb/nXp6eiTvoc3NzYiLi+Pb2lBzK2FbAplMBl9fX/j6+mLixIl6x5DJZPy+L7300kj8SRlDxN7eHklJSUhKSsJvfvMbAP1iPj8/H7m5ucjKykJeXh7a29v5z3R0dCA3Nxe5uf9NqvPy8tJLrRXj5+eHOXPmIC0tDampqYiPj7fJDBdbhP2VGZJ4eHhg5syZmDlzps767u5uFBcX48qVKygtLUVJSQlKSkqQlZUFhUKhI/ZkMhnCwsIQGRmJyMhITJgwgW9+HRYWhtDQUPj5+bF0TobFoG9wR8oFtq+vT8+4QzyVMu0QD6VSKdn7ayjQJtDmmHQIl6VquAzVcwnrtdiX+I0JFRSDPfgB/WnRloh+EEKg0WiM9igzxcXS2DpjtV/G9hH2RhNvN5QuKJV+vmPHDqOCjgpiKtoIIbxgs7Ozg1wuh1qthoODA+rq6gZtW3PrrbciKyvLaI2YOfPGpoNtE+9jbF5qCLefP3+erzGlpifGPjvYflToWjrS5erqivnz52P+/Pl4/vnnQQhBRUUFTp06hePHjyM7OxsXLlzQMcmhJQ/GqK6uxtmzZ/lluVyOqVOnYtasWUhPT8ecOXMwadIki//+4xH2jcswC2dnZyQmJiIxMVFvW19fHyorK1FWVoby8nKUlZXxIzs7G9XV1XrpOHK5HCEhIbzQow2vxdPAwEBm2MKwemgal7e394gfmxDCm3BQEw+xKYdQ/BUVFSEkJMRoy4DOzk40NDRI1nIN12Ldzs5O0oTDkCmHqf29xHVaptZryeVy9vKIAUDXuIIKGluG1osJ67UGi9bS1hxU9AkjajNmzMDx48dxzz33AACuXr2KSZMmGT3eihUrEBkZKVkDJp6n5ymsM+vp6TG4D91Paio8ti0hFnzCZTpv7jrhVDiE60ydT09PR3JyMl+KYw6rV69GeXk5Ll26BJVKhdOnT+PMmTP4+9//DqDf8Cw1NRXp6el8JG/ixIlM5A2TEWksPtIwU5TxiVqtRk1NDaqrqyWHQqFATU0NWltbJT/v7e2NwMBA+Pv7w9/fn294Tad0nlri0wLt0YAZTjCsneFeozR1Tqp5s1D8GavTotvFNV1SU6lhTm86U7C3t5c03hAbcgxWk2VKPZZU7ZVUrZWhuixhit14fdBh91HL8PTTT2PHjh2oqqrCxIkT8fLLL8PJyQmEENx2223o6urC+vXrMWPGDPT09ECpVOLDDz+0eudgahQiFoFSy+Jt4kG3nz59GgkJCSgtLcU//vEPKJVKLFy4EGvWrIGdnZ3OZ5RKJb7//nu0t7fD19cXK1asgKOjI79dfH7CITQ5kdpXaIRiaD+hYYrYNMWQoYqxfaS2E0L4aK7U/Nq1a/Htt9+C4zjeVCY3NxfHjx/HkSNHUFdXJ/lv5+HhgfT0dMydOxezZs1CamoqfHx8xvgKsk7GrLH4aMAE3Y1NT08PamtrUVNTozNVKpWoq6vj+6E1NDSgoaHBYBG+vb29Ts8z2gNN2OuMDuE66oZo7G0+exBhWDvj4RqlfbsMGXAY6+UlntJ0NkO1WFI1WcJ5YSqccPtQetyZi7l1WMZqrkztv2WoAbWxPlzG1kul4f3yyy9YuHAhv2wNqWg3AoQQvn9jS0sLvLy8+JcrMTExAIAjR45gz549UKlUePHFFwdtMD5eoffR9PR0vPTSS1i9ejUyMjLw2muv6dxfNRoNVq1ahdtvvx1RUVFobm7GkiVLxkUEeCSprKxEQUEBCgoKcPLkSeTm5qKzsxNA/31O+BJvwoQJvOsmjeTdiCKPCTrGDYFWq0VLS4tO42vaE432PBP2PqP2+LRVgzGopb24v5mHhwdaWloQHx+vZ28vHELXQ2ZMwRhrxoOgswXEtVim1mWZUoNlam2WoWVT6q/EduyWZrj1ViNdezXY+sHqsUytuTKnBstQ6p2dnR0TxCPM4cOHERERgUceeQSff/45AgMD8e6776KzsxN/+MMfeMG2ZcsWlJWV4fXXX7fwGdsWGo0G58+fx7Fjx5CTk4OcnBwolUoA/3UOpg6rAJCQkID09PRBo8VhYWFITU1FamoqfH19R/V3GG1GvW0Bg2EN0KanPj4+iI2NNflzarUabW1tfJGvsO8ZtbmnVvd02tzcjPLyct4mf/v27Sb/PCcnJ0lbe7GroXjekKMhXabD0dGRfZEzGBaAPkhbezqaKdD0KUN9uYQCUdwny9T+XML1V65c4euuBuvPZawmSzylNVmm1l4ZS7mzRYZThyVVizXYvCn7mrvOnO0jMVJSUhAZGSn596yqqkJISAj/f3zy5MnYv3+/zvWRm5uL6Oho/PrXv0ZbWxtefvllpKSkGPw3ampqQmdnp45JijlT8Tpb/f63s7NDcnIykpOT8cQTT4AQAoVCgYKCAvzyyy84evQoTp06BZVKBY7jUFFRgdLSUqN9aAkhOiIwMjKSF3cZGRmYMWPGuPRkYIKOcUMiTMccCocPH8a8efN4C3tD9vbU0p4Ouq69vR1dXV1oamrScTzs6uoasqUu5hCTAAAMEklEQVS22KVQ6FYoHOJ1QmfCwUwrBjOrsNUvFQaD0f9GnD7UjwW2EEUW1igNJv6M1UaZO6Q+Z6huytC64XzGUA0VraOiKcxSNVim1mIZ2ldq/WizZcsWPPTQQ5LbaOSa/r+g9WLCsgyVSoXs7Gy8+OKL6OrqwuOPP47vv/8e/v7+ksfcuHEjPvnkkxH9HaiwMyQEh7JOatnU/YYyxJ93dnbG9OnTkZeXx7+IN5W//vWv0Gg0fIrnjh07APSb+82aNYtP5xwvpixM0DEYQ0Qmk/HplSOJWq3mxZ0hF0Mpx0KpQY0p2tvbUV9fL2luIXyTNVyoyYMht0FTmgObYjIxWMNfcX3RYLVF49l4gsFgDB0aDWE9VS0LFYBSolBKDNLWFIQQnDhxAn/605/Q2dmJdevW8cKNbt++fTs++OADfPjhhwgICMDHH3+MsLAw/mcHBARAqVTy3xEKhQJ+fn46Lz5cXV2xfPlyLFy4EADwxhtvoLi42KCgu//++zFr1iyd8xQajAinUgYkwvXifQwdx9C2wY4xmBGKofmhDKnzIIQgLS1NZ73wupCav+uuu/DMM8/ofK83NjYiJycH2dnZyMnJwbvvvsvX7Hl7eyMtLQ0333wzli1bhsTERJt7JmCCjsGwMuzt7UekabapEEJ40wcpt0JDDoRiQwpDZhRSzYG7urrQ3NzMrxfuI56OFcM1g6DbhDU+UnU/I9mDyVgtz4ULF+Di4jLsOh+6zda+3BgMxvhhONHjxYsX45tvvkFiYiJSU1Nx3333ISkpid9eU1ODDRs2wNPTE5988gk2bdqEf/7zn7yInzJlCkpKSlBTU4PJkyfj3//+N9599104Ojryx1iyZAmOHj3KLyuVSqPta+bOnYu5c+ea/bswho6vry/Wrl2LtWvXAgB6e3tx/vx5PoJ34sQJPPfcc3juuecQFBSEOXPm8BE8W6jFY4KOwbjB4TiOT5v08PCw9OnoQN9eGjOVMGQcYcxMYjCzCGN1QVJ1PrRvm7A+SFwTRJfFPZjUarXRegBrQfhAZcyUwZyaHEPT4dTtjEX9jbCOxdz9pGpgRmN+KNPB1jEYtsbZs2cRERGBuLg4AMCdd96JvXv3YurUqXzKZGZmJn99L1u2DF999ZXePXnz5s1YvXo1VCoV1q5dizlz5uDNN9/E1KlTsXLlSqxZswY//vgjlixZAgC47bbbEB0dPYa/KcNcHB0dkZaWhrS0//qNKJVK/PTTTzhw4ADy8vL4NE0ASE9Px+23345169YZrLe0JEzQMRgMq0XYAHg8Q9OExCYNpswb2nbmzBlMmTJl2DU+xmp5hlKvY6jehgp3Y72RjB3HWB2OofQshnkYq6EZSh1OX18fnJychlSTI7V+sPodU+p7hlsDNBZD+G8xlP1Gc714nbHlkdgGAImJiUhISJC8ZquqqhAWFsZH9iIiInD06FEdwUaPAwD/+c9/MHfuXL3vnGXLlmH+/PlQqVSorKzEwYMHER0dDY1Gg4MHD/L71NbWQqPRICoqCkePHtU7Z6nfY6jTkTiG+G9gbNtQj2Xu/FA/Z+gYUvsZWpbL5Vi5ciVWrlwJjuPQ3NyM06dP4+TJk9i5cyc2btyIjRs3YvLkycjMzDSait3X14epU6fiscce06m3HC2G9ZTEcdwyAO8DsAPwMSHkTdF2RwCfA0gF0AjgTkJI2XB+JoPBYIw3RqNORy6XW73hhKURGjMYqscxVrMjXDZ0DKl9DNXAGKufMWcfYzU0xmpyTF03lHodqX0VCgWCgoJM/oxUvY0p60yp0THnGJYc9Lod6n7jkVdffdWgoNNqtTrihxCisyxk37592LNnD7KysiQfwKnR2JdffsnaE9zglJaWorS01KR9MzIykJycPMpnNAxBx3GcHYC/A1gMoArAKY7jdhNCLgp2+y2AZkJIFMdxdwF4C8CdwzlhBoPBYDBGAo4bW1dHhi624HI5XjFHAJq6XrzO2PJIbKPzxmqbQkJCUFVVxS9XV1cjKChIT9CdOnUKTz31FL7//vtB3a83bNiAlStX8j9f6pzMXWfudCSOQTFl21CPZe78UD9n6BhS+43EsqF1QrRaLVJSUsZEzAHDi9ClAyghhFwDAI7jvgawBoBQ0K0B8PLA/HYAmzmO48h4fU3EYDAYDAaDYeUYilKNN1JTU1FcXIyKigpERkZi69at+PTTT3V+95KSEjz44IPYtWsXX2tnjIiICERERIzmaTMYZjOcpM5QAJWC5aqBdZL7EELUAFoBWLdNDIPBYDAYDAZjXLB582YsW7YMUVFRWLNmDZKSkvDKK69gz549AIDHHnsM165dw3333Yfp06fjtttus/AZMxjmww01WMZx3O0AlhJCHhxY/jWAdELI7wT7XBjYp2pguXRgn0aJ420AsGFgMQFA0ZBOjMEYG/wANFj6JBgMI7BrlGHtsGuUYe2wa5RhaSIIIdINDQUMJ+WyCkC4YDkMgMLAPlUcx9kD8ATQJHUwQsgWAFsAgOO4fEJImtR+DIY1wK5RhrXDrlGGtcOuUYa1w65Rhq0wnJTLUwCiOY6byHGcA4C7AOwW7bMbwH0D8+sA/Mzq5xgMBoPBYDAYDAZjZBhyhI4QouY47gkAB9DftuATQsgFjuP+F0A+IWQ3gH8B+ILjuBL0R+buGomTZjAYDAaDwWAwGAzGMPvQEUJ+BPCjaN2fBPM9AG4fwqG3DOe8GIwxgF2jDGuHXaMMa4ddowxrh12jDJtgyKYoDAaDwWAwGAwGg8GwLMOpoWMwGAwGg8FgMBgMhgWxKkHHcdwyjuOucBxXwnHcc5Y+HwaD47hwjuOyOY67xHHcBY7j/mdgvQ/HcQc5jisemHpb+lwZNzYcx9lxHHeG47gfBpYnchx3cuAa3TZgXsVgWASO47w4jtvOcdzlgfvpLHYfZVgTHMc9PfA9X8Rx3Fccxzmx+yjDVrAaQcdxnB2AvwNYDiAewHqO4+Ite1YMBtQAniGE3ARgJoDHB67L5wBkEUKiAWQNLDMYluR/AFwSLL8F4L2Ba7QZwG8tclYMRj/vA9hPCIkDkIT+a5XdRxlWAcdxoQCeBJBGCElAv9nfXWD3UYaNYDWCDkA6gBJCyDVCSB+ArwGssfA5MW5wCCFKQsjpgfl29D+EhKL/2vxsYLfPAKy1zBkyGADHcWEAVgL4eGCZA7AIwPaBXdg1yrAYHMd5AJiHfudrEEL6CCEtYPdRhnVhD8B5oG+yCwAl2H2UYSNYk6ALBVApWK4aWMdgWAUcx0UCSAFwEkAgIUQJ9Is+AAGWOzMGA/8H4I8AtAPLvgBaCCHqgWV2P2VYkkkA6gF8OpAW/DHHca5g91GGlUAIqQbwNoAK9Au5VgAFYPdRho1gTYKOk1jHLDgZVgHHcW4AdgB4ihDSZunzYTAoHMetAlBHCCkQrpbYld1PGZbCHsA0AB8SQlIAdIKlVzKsiIH6zTUAJgIIAeCK/hIgMew+yrBKrEnQVQEIFyyHAVBY6FwYDB6O4+ToF3NbCSHfDayu5TgueGB7MIA6S50f44ZnDoBbOI4rQ3+q+iL0R+y8BlKHAHY/ZViWKgBVhJCTA8vb0S/w2H2UYS1kArhOCKknhKgAfAdgNth9lGEjWJOgOwUgesBRyAH9xai7LXxOjBucgVqkfwG4RAh5V7BpN4D7BubvA/D9WJ8bgwEAhJDnCSFhhJBI9N83fyaE/ApANoB1A7uxa5RhMQghNQAqOY6LHVh1M4CLYPdRhvVQAWAmx3EuA9/79Bpl91GGTWBVjcU5jluB/jfLdgA+IYT8xcKnxLjB4TguA0AOgPP4b33SC+ivo/sGwAT0fxHcTghpsshJMhgDcBy3AMAfCCGrOI6bhP6InQ+AMwDuIYT0WvL8GDcuHMclo9+0xwHANQAPoP+lMruPMqwCjuNeAXAn+t2tzwB4EP01c+w+yrB6rErQMRgMBoPBYDAYDAbDdKwp5ZLBYDAYDAaDwWAwGGbABB2DwWAwGAwGg8Fg2ChM0DEYDAaDwWAwGAyGjcIEHYPBYDAYDAaDwWDYKEzQMRgMBoPBYDAYDIaNwgQdg8FgMBgMBoPBYNgoTNAxGAwGg8FgMBgMho3CBB2DwWAwGAwGg8Fg2Cj/H+btAjBH1pFMAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1, 1, figsize=(15, 10))\n", "the_levels = [0.2, 0.6, 1.0, 1.4, 1.8, 2.0, 2.2, 2.6]\n", "cs = ax.contour(x1, y1, psi, levels=the_levels, colors=\"k\")\n", "ax.set_aspect(1)\n", "ax.clabel(cs, the_levels, fmt=\"%3.1f\")\n", "ax.grid(True)" ] } ], "metadata": { "jupytext": { "cell_metadata_filter": "all", "formats": "ipynb,python//py:percent", "notebook_metadata_filter": "all" }, "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.6.7" }, "toc": { "base_numbering": 1, "nav_menu": {}, "number_sections": true, "sideBar": true, "skip_h1_title": true, "title_cell": "Table of Contents", "title_sidebar": "Contents", "toc_cell": true, "toc_position": {}, "toc_section_display": true, "toc_window_display": false } }, "nbformat": 4, "nbformat_minor": 2 }