{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example #1: Creating a simple 1 tone simulation\n", "\n", " - This notebook contains a very simple SIS simulation:\n", "\n", " - The input consists of only one tone and one harmonic\n", " - The embedding circuit is ignored (i.e., no harmonic balance is required)\n", " - The DC current-voltage relationship is generated through a polynomial model\n", " \n", " - The simulation then calculates:\n", "\n", " - The DC tunneling current (i.e., the pumped I-V curve)\n", " - The AC tunneling current" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "\n", "import qmix\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "# https://github.com/garrettj403/SciencePlots\n", "plt.style.use(['science', 'notebook'])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Generate the response function\n", "\n", "- Here we will generate the response function using the polynomial model from Kennedy (1999).\n", "- The ``RespFnPolynomial`` class will automatically generate the Kramers-Kronig transform of the DC I-V curve, and then setup the interpolation." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Generating response function:\n", " - Interpolating:\n", "\t- DC I-V curve:\n", "\t\t- npts for DC I-V: 273\n", "\t\t- avg. error: 3.7010E-10\n", "\t\t- max. error: 0.0000 at v=-1.10\n", "\t- KK curve:\n", "\t\t- npts for KK I-V: 330\n", "\t\t- avg. error: 6.1923E-07\n", "\t\t- max. error: 0.0001 at v=1.11\n", "\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdwAAAFCCAYAAABbz2zGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdd3yUVfb48c9JQk3o0mtAZEVEQBSElRVRUSmKqMiCKKLSE9R10a8ooSm2XQIrgljwB4ro6gICroVmQcSVJiC9hxIInQTS7u+PO5M6mUzItCTn/Xo9r5l52pwHkTP3ee49V4wxKKWUUsq3QgIdgFJKKVUSaMJVSiml/EATrlJKKeUHQZlwReR+EflcRPaLSJKIbBeRV0SkggfHlhWR10XkiOPYn0Wkkz/iVkoppfIiwdhpSkTWAAeAhcAhoDUQA2wDOhhj0t0c+xHQDXgW2AMMB+4CbjLGbPBt5EoppZRrwZpwqxtjjudYNwD4EOhijFmex3HXARuAx4wxHzjWhQFbgO3GmJ6+jVwppZRyLShvKedMtg6/Ol7rujm0J5ACzM9yrlTgE6CriJTxWpBKKaVUAQRlws3DXxyvf7jZ5xpgrzEmMcf6LUBp4EpfBKaUUkrlJyzQAXhCROoC44HvjDH/c7NrVeCUi/Uns2x3df7gu6+ulFIqoIwx4s3zBX3CFZEIbOepVGBgfrsDrpJnvn9onjzLFpF89/Nkn0DtF8yxebqfxqaxBeN+GlvRjO3ChQvMmTOH2NhYtm3bBkDZsmXp378/7777br7nKaigvqUsImWBRUBjoKsx5lA+h5zEdSu2SpbtSimlSrjRo0dTr149hg4dyrZt26hbty4vv/wyBw8eZNasWT75zqBt4YpIKeBz4EbgNmPM7x4ctgXoJSLlczzHbQ4kA7sKE9PYsWO9so8v9gvEd+o1XJ5AxObN+AtyvkD8P+OpknIN+v+CZYzhp59+IjY2FhHhtddeA6B9+/ZER0fTu3dvSpUq5dH3XK5gHRYUgu1Z3BPoZoxZ5uFxrYD1wKPGmA8d68KA34FdxpgeeRxngvHPoSA8vZUSzIr6NRT1+EGvIVjoNXhPcnIy8+fPJzY2lt9++w2AsLAwHnzwQaKjo7nxxhtdHueIv0Q8w30LeACYBFwQkfZZth0yxhwSkYbAbmC8MWY8gDFmg4jMB6Y4Wsh7gaFAJNDPr1eglFIqYOLj45kxYwZvv/02R48eBaBatWoMGTKEoUOHUreuuxGmvhGsCfcux+sLjiWrcdiqUwKEkvs59EBsop4IVAY2AncaY9b5Kthg0KdPn0CHUGhF/RqKevyg1xAs9Bou3/r164mNjWXevHkkJycDcO211xIdHc1f//pXypUrF5C4IEgTrjGmkQf77MNF72NjTBLwtGMpMfr27RvoEAqtqF9DUY8f9BqChV5DwaSlpbFo0SKmTJnC999/D9hbwj179mTUqFHccsstiHj17vBlCcqEG4zOnj1LfHw8KSkpgQ7FpXr16vHHH+5qggQ/b11DqVKlqFGjBhUrVvRCVEqpYHX69Gnee+89/vWvf7Fv3z4AKlSowKBBgxgxYgRNmjQJbIA5aMJ1cP76GTt2LDExMdm2nT17lmPHjlG3bl3KlSsXFL+Ucjp9+jSVK1cOdBiF4o1rMMaQlJREXFwcgCZdpYqhHTt2MHXqVGbPns2FCxcAaNKkCVFRUTz66KOX/f99TEwM48aN82ao2WjCdXDXmy4+Pp66detSvnx5P0akLoeIUL58eerWrcvhw4c14SpVTBhj+Pbbb4mNjWXp0qUZ67t06UJ0dDR33303oaGhhfqOmJiYjAaXLxpWmnA9kJKSEtAH7argypUrF7S3/5VSnktMTMyoBuV85OSsBhUVFcW1114b4Ag9pwnXQ8F4G1nlTf97KVW0HThwgLfeeotZs2Zx6pQtkV+nTh2GDx/Ok08+yRVXXBHgCAtOE65SSqmgYIxh9erVxMbG8sUXX5CWlgb4txqUL2nCVUopFVDJycl8+umnxMbG8r//2QnhwsLC6Nu3L9HR0bRr1y7AEXqHJlyllFIBER8fz8yZM5k+fXq2alCDBw9m2LBhAakG5UuacFUuBw4cICoqilq1anHo0CEmT55MixYtAh2WUqqY2LBhQ0Y1qEuXLgHQokULRo0aFfBqUL4U1NPz+ZOIICK5xuCWREOHDqVv377MmDGDMWPG0K+flqFWShVOWloaCxYs4JZbbqF169bMnj2b5ORkevbsybJly9i0aRODBg0KaLKNiYnJyAW+oAnXwRiDMabEJ9yEhAS++uorunXrBtjOCocPH2bDhg0eHZ+amsr48eM5f/6812OLiYnh9OnTXj+vUsp3zpw5wz/+8Q+aNm1Kr169WLVqFRUqVCA6OpodO3awcOFCbr311qAYWRATE5ORC3xBE24JdeHCBVq3bk3lypUJCwujVatWzJw5k/3791O+fHkiIiIy9q1RowZ79+7N95xpaWncd999JCQkEBERwfHjx2nVqhUiQnh4OK1ateLAgQPZjpk8eTL169enevXq9O/f3+35GzZsSIcOHThx4sTlXbRSym8OHz7MyJEjqVu3Ls888wx79+6lSZMmTJkyhUOHDjFlyhSuvPLKQIfpV/oMt4QKDw9n/fr1tGvXjpSUFNats5MpOeeLvByTJ09m7969LFy4EIDq1auzYcMG6tWrB+Cylfzcc8+xYsUKxowZk+8A9oEDB/L9998zaNCgjO9QSgUPYwzfffcdsbGxLFmyJGP9rbfeyqhRo7xSDaoo0xZuCZaens7mzZtp1apVxrpGjRqRmJiY7ZZwfHw8kZGRbs+VkpLCm2++yVNPPZXr1lDLli2Ji4tz2TJdunQpERER3HzzzR7F/Oyzz7Jo0SI2b97s0f5KKd9LTExk5syZtGjRgjvuuIMlS5ZQunRpHn/8cTZt2sSyZcvo0aNHiU62oAm3RNu+fTuJiYnZEm61atW48847M36drlmzhtq1a2fbx5Xly5dz6tQpbrnlllzbrrvuOgA2btyYbX1aWhovvPACr776qscxN2/enFq1avHpp596fIxSyjcOHjzIc889R/369RkyZAhbt26lTp06TJo0iXfffZdZs2YVqdKLvqYJtwRz3uLNmUzffvttPv74Y4YMGcKECROYO3duvudasWIFYWFhNGrUKNe2li1bZvs+p1mzZtG5c+cCP8dp1qwZy5YtK9AxSinvcFaD6tOnD5GRkbz66qucPHmSG2+8kY8//ph9+/bxf//3fzpxiAv6DLcE27BhAyKS0QJ1atiwYYGfkR49epQqVaoQEpL7N5yrFu65c+eIjY1l9erVBY67evXqhXrWrJQquOTkZD777DNiY2P59ddfAVsN6qGHHiI6Opr27dsHOMLgpwn3MgVDF3ZwP61gfjZu3EhkZCSVKlUqdBzx8fGULVvW5bZmzZpRpkyZbC3cV155hcGDB1OlSpVs+3pSdKNs2bIcO3as0DErpfJ3/PjxjGpQR44cAYp3NShf0lvKDiWx8MWGDRvyfTbrqdDQ0DyTf2hoKNdccw3btm0jOTmZgwcPsnjxYoYPH55rX0+KbhhjCAvT34pK+dLGjRt57LHHqF+/Pi+++CJHjhyhRYsWzJo1i4MHDzJp0qRil2y18IWfFLTwhXP/QC+X69ixYxw7dizPhHvu3DkeeeQRl89kXalZsyZJSUl5bm/ZsiUpKSls2bKF559/nnHjxuWa9ePkyZMeFd1ITEykZs2aHsWllPKcsxpU586dadWqFR988AHJycn06NGD7777jk2bNvH4448X29KLWvhC+cT69euB3B2mnCpUqMC4ceM8Pl+9evU4depUnpO+O5/jzpo1i0OHDtGrV69c+xw8eNCjohtHjx6lfv36HsemlHLvzJkz/POf/8yoBrVy5UoiIiKIiopix44dLFq0iC5dugTNo7SiSu/LlVB59VCeMWMGc+bM4aqrrsrVe/jf//43n332GfXr12fr1q1ER0fTtWtXAO666y7GjRvHzp07ad68ea7vc/ZUnjFjBmvXrnUZkye/Ko0xbNu2jTFjxuR/kUopt3bu3Mm0adP44IMPMsbeN27cmJEjRzJw4ECv9O9QmTThljC7d++md+/exMXFAdCjRw/Gjh1Lr1692Lx5M2PGjGHPnj1UrFiRd955J+O47du3ExUVxe7duylXrhw//PADq1atyki4N954Iw0bNuSrr75ymXCdLdx+/frRtm1bl7E1aNAgo+iGs5Wbs+jGTz/9xJkzZ3jggQe88weiVAljjGHZsmVMmTKFpUuXZvzQ7dy5M6NGjaJbt24lvkCFr2jCLWGaNGmS50QEy5cvp23bthnj5zp16sTLL78MwLfffkvLli0znt3cfPPN2apDiQiTJk1i9OjRDBs2LNcznmrVquXbgq1atWpG0Y0+ffq4LLoxceJEhg8frreUlSqgxMRE5s6dy9SpU9myZQsAZcqUoX///kRFRWXchVK+owlX5Snr8xpjTL7Pb/r168eWLVu47777WLx48WX9Sn777beJiopixYoVHDx4MFvRjdGjRxMWFsYbb7xR4PMqVVIdOnSIt956i3feeYeTJ08CUKdOHYYNG8aTTz5J9erVAxxhyaEJV2W49dZbmTBhAufOnaNChQr8+OOPGdvuuOMOJk+ezMWLFylbtiwrV67kt99+45lnnsl2jpdffpnFixeTnJx8WT0Z3RXduOGGG3jllVdcFtdQSmUyxrBmzRpiY2P597//TVpaGmAf/YwaNYrevXtTunTpAEdZ8ui/XA4lcRxuTi1atGDChAl069aN4cOHs3fvXk6ePMnEiRNp1qwZsbGxPProozzzzDPMnj2bQYMGuTxP9+7dfTJs4P7779dkq5QbycnJfPzxx7Rr144OHTowf/58APr06cPPP//ML7/8Qt++fTXZ5sHX43C1hevgq3FXRc2QIUMYMmRIxueJEydmvL///vu5//77AxGWUsoNV9WgqlatmlENyjlFpnIvJiYmo9Hli6SrCVcppYqoTZs2ERsby0cffcSlS5cAuOaaa4iOjqZfv36UL18+wBGqrII24YpIPWA00Ba4DigHRBpj9nlwbF7N1dbGGNdddJVSqghIS0tj8eLFxMbGsmLFCsC2xrp37050dLQWqAhiQZtwgSuBB4HfgB+AOwp4/GxgZo51OwofllJK+d/Zs2d5//33mTZtGnv27AEgIiKCgQMHMnLkSJo2bRrgCFV+gjnhfm+MqQkgIo9T8IQbZ4xZ4/2wlFLKf1xVg4qMjCQqKkqrQRUxQZtwjTHpgY5BKaUCIWs1qCVLlmSrBhUdHU337t21GlQRFLQJ1wuGisizQBqwBhhrjPkhwDEppVSekpKSmDt3LhMnTuTAgQOArQbVr18/oqKiMkqkqqKpuCbcucBi4DDQEHgWWC4itxtjVgYyMKWUyunQoUNMnz6dd955h4SEBABq167NsGHDGDx4sFaDKiakKIw/dTzDnYWHvZRdHF8B2AwcNMb82cV2l38Iffr0oW/fvtSrV48mTZoU9GtVgO3evZtDhw4FOgyl8rR9+3YWL17M6tWrM6pBNW3alO7du9OhQ4dcc0Yr75s3b15GgZCcjDHe7e4d6AnUPZxk/XHAAI0KcY7pwKU8thl3tm7d6nZ7MDh16lSgQyg0b1+Dv/+7LViwwK/f5wt6Db536dIl89FHH5kbb7zROP5dM6GhoebBBx80q1evNunp6UF/DZ4o6tfgyAtezWXF9ZayK4L9y62UUn534sSJjGpQhw8fBqBKlSoZ1aB0Bqzir0QkXBGpCHQDfgl0LEqpkuX333/PqAZ18eJFAJo3b050dDT9+/fXalAlSFAnXBFxFu693vF6l4gcB44bY1aJSENgNzDeGDPecczfgGbACjI7Tf0NqAX082f8SqmSKS0tjSVLlhAbG8vy5csz1nfr1o3o6Ghuu+02rQZVAgV1wgU+y/F5uuN1FXAL9jZxKNlnPdoO9HIslYCzwE/AIGPMWl8GW1wcOHCAqKgoatWqxaFDh5g8eTItWrQIdFhKBb2zZ8/ywQcfMHXq1IxqUOHh4RnVoK666qoAR6gCKagTrsmnh5ixPZYlx7ovgS99GFaxN3ToUAYMGECfPn1Ys2YN/fr1Y+PGjYEOS6mgtWvXroxqUOfOnQOgUaNGREVF8dhjj2k1KAXofLgqh4SEBL766iu6desGQPv27Tl8+DAbNng250Nqairjx4/PKEHnTTExMZw+fdrr51XqchhHNaiePXty1VVXMXXqVM6dO8df/vIX/vOf/7Br1y6eeuopTbYqgyZch5I2Af2FCxdo3bo1lStXJiwsjFatWjFz5kz2799P+fLliYiIyNi3Ro0a7N27N99zpqWlcd9995GQkEBERATHjx+nVatWiAjh4eG0atUqo3qO0+TJk6lfvz7Vq1enf//+bs/fsGFDOnTowIkTJy7vopXygqSkJN59911atmzJbbfdxpdffknp0qUZOHAg69evZ+XKldx7771aerEI0gno/cQUgQIg3hQeHs769etp164dKSkprFu3DoDffvvtss85efJk9u7dy8KFCwGoXr06GzZsyJj82lUr+bnnnmPFihWMGTOGa6+91u35Bw4cyPfff8+gQYMyvkMpf4mLi2P69OnMnDkzoxpUrVq1MqpB1ahRI8ARqsLy9QT02sItwdLT09m8eTOtWrXKWNeoUSMSExOz3RKOj48nMjLS7blSUlJ48803eeqpp3L9RW3ZsiVxcXEuW6ZLly4lIiKCm2++2aOYn332WRYtWsTmzZs92l+pwvrll1/o27cvjRo14uWXXyYhIYG2bdsyZ84c9u/fz4svvqjJVnnEbcIVkdIi0l5E7hORfiLSVUQa+Sc05Wvbt28nMTExW8KtVq0ad955J0uWLAFgzZo11K5dO9s+rixfvpxTp05xyy235NrmLLies+NVWloaL7zwAq+++qrHMTdv3pxatWrx6aefenyMUgWVkpLCvHnzaN++Pe3bt+eTTz7BGMMDDzzATz/9xNq1a+nfvz+lS5cOdKiqCMl1S1lEQrFDah4H/gKUJntPYCMiccA8YJYxZpc/AlXe57zFmzOZvv3220RFRbFixQoOHjzI3Llz8z3XihUrCAsLo1GjRrm2tWzZMuP7unTpkrF+1qxZdO7cmSuvvLJAcTdr1oxly5Yxfvz4Ah2nVH5OnDjBO++8w/Tp04mLiwNsNagnn3ySYcOG0aBBgwBHqIqybAlXRHoDk4H6wNfAGGA9cBxIAqoCkUA74F7gaRGZDYwxxhzzX9hBwt09/pkz4ckn7ft33oHBg/PeN+vz4+uvB8fz1FyeeMKeC+C33+y+hbBhwwZEJNeUXw0bNizwM9KjR49SpUoVQkJy3zRx1cI9d+4csbGxrF69usBxV69evVDPmpXKafPmzcTGxjJ37tyMalBXX311RjWo8PDwAEeoioOcLdxpwGvAbGNMXuMv1gLzscm2HTAaeAKY6LMolU9s3LiRyMhIrwxbiI+Pp2zZsi63NWvWjDJlymTrNPXKK68wePBgqlSpkm1fT4pulC1blmPHSt7vO+Vd6enpLFmyhClTpmSrBnX33XczatQorQalvC5nwm1sjLno6cHGmF+A+0SkjHfDKiI87dn85JOZrd38eNpyK2TrFmwLt2PHjoU+D0BoaGiePb1DQ0O55ppr+P3330lOTubYsWMsXrzYZSvVk6IbxhjCwrSDvbo8Z8+eZfbs2UydOpXdu3cDttf+o48+ysiRI2nWrFmAI1TFVbb7fwVJtjmOu+SdcJS/HDt2jGPHjmV7fjt9+nTq1q3LypUrc+3vbhtAzZo1SUpKyvP7WrZsSUpKClu2bOH5559n3Lhxueb6PHnypEdFNxITE6lZs6aHV6qUtXv3bkaNGkW9evWIjo5m9+7dNGrUiDfffJNDhw7xr3/9S5Ot8ikdFuRQ0gpfrF+/HsjeYWrYsGE0bdrU5f7utgHUq1ePU6dOkZKS4nK78znurFmzOHToEL169cq1z8GDBz0qunH06FGdykx5xBjD8uXLueeee2jatCmxsbEZ1aC++OILdu3axdNPP03lypUDHaoKAr4ufJEt4YpIuoikebr4JKIAcU4QXFISbl49lJ2OHz/O9ddfT7du3fjqq6/y3XbXXXeRnp7Ozp07XZ7P2VN5xowZvPHGGy738aT4iDGGbdu2ZbSClXIlKSmJ9957j+uuu44uXbqwaNEiSpUqxaOPPppRDapXr15aDUplExMTk5ELfCHng7DxZE7SLsBjQDnsZADHsFPcdcf2WH7PJxEpn9q9eze9e/fOGPLQo0cPxo4dm6vF+dtvv3H33XczYcKEXOdwte3GG2+kYcOGfPXVVzRv3jzXMc4Wbr9+/Wjbtq3L2Bo0aJBRdMPZys1ZdOOnn37izJkzPPDAAwW8clUSxMXFMXfuXAYNGpStGtTQoUMZPHiwPopQAZUt4RpjYpzvRWQMsB/oaoxJzLI+HDtkKNVPMSovatKkSb4TESxYsID58+ezbds2j7eJCJMmTWL06NEMGzaMcuXKZdterVq1fH81Vq1aNaPohrPTVM6iGxMnTmT48OF6S1lls3btWqZMmcJnn31Gaqr9p+n6669n1KhRPPjgg1qgQgUFd89wBwOvZ022AMaYC8AbwBBfBqYCp1KlSvTp04eRI0cWaFu/fv0YMGAA9913H2lpl/fE4e233+bjjz9myJAhTJgwIVvRjdGjRxMWFpbnLWlVsqSkpPDJJ59w00030a5dO+bNm0d6ejodOnTgxx9/5Ndff9VqUCqouBtbcQW2ypQrpYFq3g9HBYPOnTvTvn172rZty+eff07v3r092gbw8ssvs3jxYpKTk3O1cj3hrujGDTfcwCuvvOKyuIYqORISEnjnnXd46623slWDeuKJJxg+fDjr16/32nA3pbzJXcL9HzBORH42xsQ5V4pIXSAG+NXHsSk/++CDD9i5cyexsbE0btyY2rVrM3jwYHbu3EnNmjXz3Pbcc89lO0/37t19Et/999/vk/OqomHz5s1MnTqVOXPmZFSD+tOf/kR0dDQPP/xwRjUoZw98pYKNu4QbBSwHdovIGmynqZpAeyAR+Kvvw1P+NHDgQAYOHJjx+dtvv821Pa9tSvlCeno6S5cuZcqUKSxbtixj/V133UV0dDS333673vFQRUaef1ONMeuBK4E3gTTgWsfrG0BTY4z7njdFTEkbh6tUMDt37hxTp06lWbNm9OjRg2XLllG+fHmGDRvGtm3bWLp0KV27dtVkq7wqoBPQG2MSgBd88s1BpqRNQK9UMNqzZw/Tpk3jvffe49y5c4B9rj9y5EgGDRqkBSqUT/l6Avp8C9KKyBXY28jVgC+NMSdFpCyQbIxJ93pESqkSxRjDypUriY2NZdGiRRk/fjt16kR0dDQ9e/bU2tmqWMjzb7HY9P4aMBLbK9kANwAngYXAj0DuqghKKeWBixcv8vHHHxMbG8umTZsAKF26NH379iU6OprWrVsHOEKlvMvdz8bngRHY6lPfAr9k2fYl8DCacJVSBXT48GGmT5/OzJkzOXHiBGAnvxg6dChDhgzRalCq2HKXcB8HxhtjXhGRnAVHdwFNfBdW8DHG6NyYRYg+kw8+a9euJTY2lk8//TSjGlSbNm0yqkGVKVMyZ/lUJYe7hFsXWJPHtmQg3PvhBKdSpUqRlJRE+fLlAx2K8lBSUlKu6f+U/6WmpvL5558TGxvLzz//DEBISAi9e/dm1KhRdOzYUX/IqhLDXcKNA1oAK1xsuw7Y62J9sVSjRg3i4uKoW7cu5cqV038ggpgxhqSkJOLi4vTWZACdPn2ad999l6lTp3Lw4EEAKleunFENqmHDhgGOUCn/c5dwPwNeEpF1ZLZ0jYhcBTwDvOPr4IJFxYoVAfvsKa/5XgMtMTGxyLfAvXUNpUqVombNmhn/3ZSHvv+exosXwx13wGWU5QQ7rCc2Npb333+f8+fPA9CsWTOio6MZMGBARjUopUoidwk3BugAfI+dNQhsEq4PrAYm+zQyP3O2WseOHeuy+EXFihWD+h/whQsXcs899wQ6jEIpDtdQZJ09Cz16cO3Zs/Drr7BmDZQt6/Hhq1ev5s0332TBggWkp9vRgrfeeitPP/00d911lxaoUEVCTEwM48aN89n53VWaSgJuAR7FJtjvsPWTnwRuN8Yk+yyqAChpE9ArlU2pUhAdbd9v3AjffOPRYT/++CNdunShY8eOfPHFF4SGhvLII4+wYcMGli1bRrdu3TTZqiLD3xPQZ2OMSQPmOBalVHFVrhyMH88fe/Zw9UcfwYIF0LNnnrvv3LmTESNG8I0jMVeqVInhw4czYsQIateu7a+olSpS8vzpKSLvi8jf8tjWWETe911YICL1RGSaiPwsIokiYkSkkYfHlhWR10XkiIgkOc7RyZfxKlUcHGnXzr5ZtAgcQ3eyMsbwxhtv0LJlS7755hsqVqzI2LFj2bdvH5MmTdJkq5Qb7u71PAq8KiLzRSTnALnqwCM+i8q6EngQOAX8UMBj3wOeAF4CugNHgK9FpJVXI1SqOPjtN3jxRfjpJ87Vrw+NGkFCAmzZkm231NRUHn74YZ599lkuXrzIgAED2LNnDzExMVrjWCkP5Pdw5QXgNmCliFT3QzxZfW+MqWmMuRvbWcsjInIddurAp4wxs4wxy7CJ+wC2apZSKquvvoKJE+HTT0EE2reH6tXhyJGMXYwxDB06lI8++oiIiAgWLlzIhx9+SLVq1QIYuFJFS34VwZcDnwOLgbUi0s0Ys9X3YUEhJkboCaQA87OcK1VEPgGeE5EyxphL3ohRqWJh7Vr76ryd/N579plulvHmn3zyCe+++y7lypXj66+/pkOHDgEIVKmiLd/ug8aYnUA7YA+wWkS6+jyqwrkG2GuMScyxfgt2EoYr/R+SUkHMMXEAbdrY1/LlsyXbs2fPMmLECABiY2M12Sp1mSSv7s8ikg60N8asdXwOBaYDA7G9lh81xuSsseybIEUeB2YBkcaYffns+w1Q0RjTPsf627CTMHQyxvyQY5vLP4Q+ffrQt2/fwoSuVFALuXSJ7g89hAkJYfH8+Zgs0+lAdrAAACAASURBVOCFXLqECQvj3wsWMGfOHK6++mpefvllrbSmipV58+Yxf/58l9uMMd79y+4cc5RzAdKBG12sfwpIBdLyOtbbC3YiBQM08mDfb4GfXay/3XGOm11sM0XdggULAh1CoRX1ayiS8W/caAwYc9VVxpgs13DvvcaEhJiU1atNjRo1DGC+/vrrAAbquSL53yEHvYbAc+QFr+Yyd89wOwN/uEjQ/xSRX4CmXsj3vnASaOBifZUs25VSANu329dmzbKvDw+H9HS2ffIJ8fHxXH311dx+++3+j0+pYsRdpalVxphzeWxbbYz50HdhFcoWIFJEchblbY6d5WiX/0NSKohdfTVce232da3sCLp4R2GLPn366K1kpQopWwtXRF4C3jXGHHa8d8cYY4JxAvpFwDjgAeBDABEJA/oA3xjtoaxUpgcesEtOjoQbvnMnAA8++KA/o1KqWMp5SzkG+C9w2PHeHQP4NOGKyP2Ot9c7Xu8SkePAcWPMKhFpCOwGxhtjxgMYYzaIyHxgioiUwk4jOBSIBPr5Ml6lio3rrgPgmrQ0GjVowNVXXx3ggJQq+rIlXGNMiKv3AZSz4MV0x+sq7MQKAoSS+9b4QGASMBGoDGwE7jTGrPNZpEoVNcbAiRNwxRXZhgEBUL065ytXJuL0aXq30gJtSnlDfoUvAsrk0yXb2CFCufYxdqajpx2LUsqVo0ehTh1o3Bh27861eUe5crQ5fZqutWoFIDilip9gaMUqpQLB2UO5uuuqra8lJ9MbaNC/v/9iUqoYy5ZwRSRdRNI8XHJPJVKEiQgiovPhqpIjryFBQEJCAvMTEvg6PJymHTv6OTClAiMmJiYjF/hCzlvK47GdoUoc46MJh5UKWm4S7hbHTEHNmzfXCeRViRETE5PR6PJF0s3ZaSrG69+glApOHiTcIaVKweOPw+TJtnOVUuqy6U9XpUoqDxLuHYcP29mD1q/3Z2RKFUtueymLSGngLqAZUDbH5mAtfKGUys+lS7B3L4SEwJW5J9ByJty0a66BfftgwwbQ0o5KFUqeCVdE6gA/Ao2wz3WdN7SzPuzUhKtUURQSAt99B4cOQdmcv6Vh61Y77XXEn/8MS5bAxo3+jlCpYsddC/d14DjQCTiAnRP3OPAYtkziHT6PTinlG6VKQefOLjdduHCB+Ph4SpcuTRXnPhs2+DE4pYond89wbwbexJZ5BEg3xuwzxrwE/BuY6uvglFL+t2/fPgAaNmxISMuWEBoKf/wBFy4ENjClijh3CbcacNgYkw5cIHN6O4Dl2NKKxYaOw1UlyltvwYsvgmNygqycCTcyMhLKlbMzCaWnw//+5+cglfIvf4/DzeoQ4BwHsBt7C/k7x+cbgYs+iShAdByuKlHmzIFffrEdoZpmn9p67969ADRq1MiuuOMOOyRIx+OqYs6v43BzWAH8BVgAzATeEpFWQArQ1bFOKVXUGAPbttn3LoYEOVu4GQn31Vf9E5dSxZy7hDsGqApgjHk7y5yy5YHXsFWplFJFzbFjcOYMVK4MNWrk2uxs4UZGRvo7MqWKtTwTrjHmBHAiy+dpwDR/BKWU8iFnwYs//Sn3tHy4aOECpKTAr79CkyZQs6bvY1SqGNKHMkqVNG5uJwMcOHAAsL2UMzz2GHTsCF984evolCq23CZcEblTRN4TkW9E5Pscyyp/BamU8qKsLdwcUlJSOHHiBCEhIdTIerv5z3+2rytX+j4+pYqpPBOuiPwdWAp0B8KBtBxLuj8CVEp5Wa1acM01dsnhzJkzANSoUYPQ0NDMDbfcYl9XrLBDhJRSBeau09QIbE/kEcaYND/Fo5Tytb//3S4unD59GoDatWtn33DVVdCgARw4AOvWQdu2vo5SqWLH3S3lisBnJSXZauELpeDkyZMA1KpVK/sGEbj7bvt+6VI/R6WUf/i68IW7hPs10N4n3xqEjDEYYzThquLt5Ek4fNiOxXXB2cLNlXAhM+EuWeKr6JQKqJiYmIxc4AvuEu4I4G4ReV5ErheRxjkXn0SklPKduXOhbl2Ijna5+dSpU4CLW8oAt95qZxaKj9e6ykpdBnfPcA1wDpgETMxjn9A81iulgtGmTfY1RzlHJ2fCddnCDQ+H33+3Y3F9dMtNqeLMXcKdDXQA/glsA5L9EZBSyofWrbOvrVq53JxnpyknF5PVK6U84y7h3oLtoTzbP6EopXzq4kXbQhWB1q1d7uK2hZvVgQN2Tt28ErNSKhd3z3BPAMf8FYhSysc2bYLUVLj6aoiIcLmLRwl32jSIjITXX/dFlEoVW+4S7lRgmIho+UeligPnfLZ5jKE1xrjvpex00022+MWHH9pWs1LKI+6SaRWgBbBVRKaJyPgcyzg/xegXOg5XFXvOhHv99S43nzlzhuTkZCIiIojIowUM2ITdpo0dYvTvf/sgUKUCI5AT0L+Q5f1VLrYbYKx3wwkcnYBeFXvTpsGAAbaXsQtHjx4F3HSYymrwYLtMnw79+mmvZVUs+HoC+jxbuMaYkHwWHRKkVFESHm5rItev73LzkSNHAA86TAH89a9QpQr8/DMsX+7FIJUqvlwmXBEpLSL/EZFO/g4oSwz1ReTfInJGRM6KyBci0sDDY00ei+uxEEqpjBauRwk3IgL+9jf7/qWX8qxcpZTK5DLhGmOSgdvy2u5rIlIeWA78CXgEeBhoCqwQkXAPTzMbuCnHssPrwSpVFEyeDD16wKq8Z9Us0C1lgJEjoUYNe4s6KckbUSpVrLl7hvsTtpbySv+Eks0TQGOgmTFmF4CIbAJ2AoOBf3hwjjhjzBrfhahUEbJokb39O3RonrsU6JYyQIUKsGuXfVVK5ctdC/YZYJCIjBCReiISKiIhWRcfxtUTWONMtgDGmL3YHwH3+PB7lSp+zp2DtWshNBRuvjnP3QrcwoXsyVZvKyvllruk+TvQBIgF9mNLO6ZkWXxZ6vEaYLOL9VuA5h6eY6iIXBKRRBFZLiJ5/0ujVHH2ww+QlgY33ui2NVqgZ7g5LV8O7drZmYiUUi65u6U8Hjv0JxCqAqdcrD+JHR+cn7nAYuAw0BB4FlguIrcbY1Z6K0ilioRly+zrrbe63a3At5SzmjIFfv0VRoyAzz/XYUJKuSDBOP5URJKBN40xz+dYPwkYbYxx90PB1fkqYFvMB40xf3ax3eUfQp8+fejbt29BvkqpoNNl2DAiDh/mx4kTSWjRIs/9BgwYwNmzZ5k9ezaVK1cu0HeUi4+nc3Q0pZKS2DB0KPu7di1s2Er5xbx585g/f77LbcYY7/5ydE62624BIoD6QLgn+xd2wdZwnuli/XTg+GWeczpwKY9tpqhbsGBBoEMotKJ+DUEZ/7ZtxoAxVaoYk5KS526XLl0ygAkJCTGpqamX910ff2y/q2xZYzZtusyACy8o/zsUkF5D4Dnygldzm9uOTyLSVUT+B5wG9gFnRGStiNzu1ayf2xbsc9ycmgNbL/OcQuBukSsVGFdcYStMPf88hOV9Yyg+Ph6ASpUqERp6mTVt+vaFxx6z9ZV79YITJy7vPEoVU3kmXBHpCizBtm4nAMOwE9FXAJb6OOkuAtqLSOMs8TQCOjq2FYiIVAS6Ab94KT6lioZq1exz1Wefdbubs8NUlSqedJFwY+pUO/Xf7t3w4IPac1mpLNw9C40BvgG6G2PSnStFZDy2Q9I44FsfxTULGAEsFJEx2JbpBOAgMDNLLA2B3cB4Y8x4x7q/Ac2AFWR2mvobUAvo56N4lSrSnB2mCvrsNpfwcPjyS7jjDhg9WjtPKZWFu4R7HfBA1mQLYIxJF5HpwKe+CsoYc0FEbgX+CczB3g5eBowyxpzPsqsAoWRvqW8HejmWSsBZ7PjdQcaYtb6KWamgs2SJrSzVvz+0bOl2V6+1cAHq1rVz717urWmliil3CfcSUDGPbRUc233GGHMA6J3PPvuwSTfrui+BL30XmVJFxKxZsHAh1KuXb8J1tnC9knAhe7L9739h/nyYMQPKlPHO+ZUqgtwl3JXABBFZY2yVJwAcEwjEYG/ZKqWC0YkTsHQphITYZ6n58GoLN6uzZ+Hhh208+/bBf/4Dhb1trVQR5a6X8mjsLdntIvK9iMwXkVXYesaVHduLDZ2AXhUr8+dDSop9lupBIQufJdyKFeHbb6F2bVi5Ev78Z5t4lQpCvp6A3t18uDuAlsBUoAzQBiiLLfXYyhiz0ycRBYhznJQmXFUszJljXx9+2KPdvdZpypVWrWDNGmjeHLZsgTZtbMcqpYJMTExM1voMXud2HK4x5ogx5m/GmHbGmKaO178bY474JBqlVOHt2AG//GLnrL33Xo8O8foz3JwaNIAff4Tu3eHUKejZE954wzffpVSQCsh8t0opH/rwQ/vauzeUL5/v7unp6Rx2TDpQtWpV38VVpYrtxPX663b4UD61nZUqbtzWJBaRR4C+QAPs7eSsjDGmia8CU0pdpgcfhCNHYNAgj3Y/ceIEKSkpVK1alTK+7kUcEgJ/+xs8+qitguU0Zw7cd59NxEoVU+4qTb0IfADUATYAq3Is3/sjQKVUAV13Hbz/PnTs6NHucXFxANSpU8eXUWWXNdl+8QUMGADNmsHcuZCenvdxShVh7lq4g4BYY8xT/gpGKeV/zoRbt27dwARQvz5cfz389pvt5DVtGkyaBF26aKUqVay4e4ZbDS0goVTRsXYt3HxzgXsAO5/fBizh3nCDjf2DD+wQprVr4fbb7RCiFTrcXxUf7hLuKmx5xxJBx+GqIu/1121P4O8L9rQnILeUcwoJsc91d+yAl1+2ky6sXm2vRyk/Cdg4XGAUMFBEBojIFSISknPxSUQBouNwVZH2xx/w+edQujQ8VbCnQAG/pZxVhQp2KsG9e+G11yAqKnPbjBkweTIkJAQuPlWsBXIc7g6gBbbj1DEgJceS7JOIlFIFN3mynQpv4EAoYEs1qBKuU4UKdkrBSpXs55QUGD/eJuO6dW1P7CVLIDU1sHEqVQDuOk2NRydsVyr47d0LH31kJwwYXfCKq4cOHQJswj148KC3o/OO0FD7jDc21k6G8NlndqlZE/r1g5EjoVGjQEeplFt5JlxjTIwf41BKXa5x4yAtzQ6tiYws0KHGGPbutXOTREZGBm/CDQmBrl3tcvCgHT704YewfTv84x+2cpUj4ZY+e9YOLQopVk+9VDGgfyOVKsrOnrWzApUqBWPHFvjw48ePc+HCBSpXruy7so7eVr++vbX8xx+2RvPTT9sezQ5tpkyxUxI++aR9rn3qVACDVSpTtoQrIk+JSM6KUm6JSBsRudO7YSmlPFKxIuzcaYtHNG5c4MP37NkDQOPLODbgRKBdO3jzzcz5d1NTCT961FbamjUL7r/fFtno0AFiYmySVipAcrZwBwD7RGSyiOQ5JEhEqojIwyLyDfAjeU9Ur5TytUqV7KQAl6FIJ1xXwsJY9tZbtojGpEnQqZO9tfzzz/bW+w8/ZO67dSssXw4XLgQuXlWi5HyG2wZ4GHgG+LuInAV+B44Dl4AqQGOgiePzfKC5MWafvwJWSmF7506fDo8/7tEEBXlxPr8tNgkXbMu3TRu7/N//wblztoDGN9/YZ8BOs2bBlCkQFmb3bd8e2ra1y1VXZbaalfKSbC1cY/0/Y8x1wE3AP4Fz2CTbGqgA/AA8BtQxxgwsLslWC1+oImXGDIiOhm7dCnUaZws3soCdrYqUChVsp6p//QsaNsxc36CBTbTp6ba61dSptuNZ8+bZZzJKTbXz+Kak+D925Ve+LnzhrpfyL8AvPvnWIOSrgc5Ked3x4/Dii/Z9dHShTrV7926gmLVwPfXUU3Y5e9bOH/y//2UuLVpk7rdzp/1cqhT86U/2fYsWcO219rVhQ+0RXUzExMRkNLp8kXTdTs+nlApCUVFw+jTccQfcc89ln8YYw5YtWwD405/+5K3oip6KFW3t5ttvz1yXtTV78iQ0aQJ79sDvv9slq40boWVL+37BAjhzxt6SvuoqW6JSKQdNuEoVJQsWwCef2Oe2b79dqNl0jhw5wokTJ6hcuTL169f3YpDFQKlSme87doRdu2znqq1bbcLdvNm+bttmW71OU6bAqlWZn6tUgaZN4cor4c477WxIYMdNp6RA2QINClFFnCZcpYqKkydhyBD7fvLkyxoGlNWmTZsAaNmypc+eWRUr4eF2ZqMbbsh7n3vugdq17SQMO3faMcBr19qlcuXMhLtxo52SsHZte0u6USO71K8PdeoQdvGiP65I+ZkmXKWKio8+gmPHbJGH4cMLfbqNGzcCcN11JWZSMN/LOnGEMfa/144d9nZ006aZ244ds72gjxyxy5o12U5TesaMzA8jR9rtderYOtJZXyMj7a1rVSRowlWqqBgxAmrUgNatvdJJJ2sLV/mAiJ3ft1YtOx44q7vugosXIS4O9u+HffvsEhcHhw9zMWvVr99/tx25XLn3XvjPf+z7I0egc2f7d6RGDahePfN9jRpwyy12HdgfA3pXw+/yTLgi0glYZ4w572JbBNDGGFOwiTeVUpdPBPr08drpfv75ZwBat27ttXOqAggLs7eTGzbMlZDTFy7M/PDRR7Z+tCMZZ7wePmyHNTkdPWprS2/f7vr7fvghM+E+9RS89x5UrZp7ueoqeOaZzON+/NHeDq9c2RZZCQ/XXtmXyV0LdwV2LO5aF9uaObYXm5HhzmdYY8eO1bG4KnikpNiawM88k32oSiHt2bOHvXv3UqVKFVq1auW18yofqFvXLvlp3tyOF46Pd71k7Rh36hScP2+XAweyn+fGGzMTbloa3Hxz9u0itmd3pUrw+ut2qkSwxUXmz7frK1Ykcv9+O+TKue9f/pLZqr54EcqUCbpWdkxMDOPGjfPZ+d0lXHd/EmWANC/HElA6DlcFpZdegtmz4aefbB1gL1U/+vbbbwHo0qULoVpRqXgoU8Ym3ebN89939myYNg0SEmxnvKxL5cqZ+yUl2T4DCQl2uNOZM7a3tvN9WpY0sG4dzJyZ8bEl2GpeYHtjJyVl7tumje3hHR5ul4iIzOWvf4Vhw+x+e/bYc0ZE5N4vPBxuuimz0trp0/a1XDkoXfqykrlfx+GKSCNsVSmnto7bx1mVw1aayvGzSCnlVV9+aXsjh4TY239eTIz//e9/Abjtttu8dk5VhDhbqRUrup/SMSIie/1psJW3zp2zCTfrOOPbb4e33rKt2jNn2LthA5FVq9rPOW9Bp6ba58jOVvaxY5nbsraod++G117LO75du+wYabB3gj77LPP6ypWzS/ny9vn1//t/dtv58/DQQ3Z91n2c7++/P/tQLy/K2cJ9BBiLnXjeANPI3tI1js+pQOG7SbohIvWxpSVvd3znd8AoY0y+id4x49EEoD9QGdgAjNZnzqrI2LrVTqwOMHFi7tt6hXDkyBEWL15MaGgo3S9z0gNVgoWF2fHFOadzbNkyswAIsGnhQiLzKsyyY4dtHV+4YBOg8/X8edv72qlJE/uj07kt637nz9tb1U5ly9rPSUmQnAyJiXZJSLCL0/nzsGRJ3tfXooXfEu5sYCU2wS3HJtWtOfa5BOwwxpz0SUSAiJR3fP8l7I8AA0wEVohIS2NMftN7vAd0A54F9mCv42sRuckYs8FXcSvlFadO2fGc587ZTlLPPefV07/zzjukpqZy3333UdeTZ4NK+UJoaGYrOy+NG8Po0Z6dz9mCBZvMk5LskphofyQ4VaoECxdm3571vQ+rrmVLuMaY/cB+ABHpjO2lfM5n3563J7C3tpsZY3Y54tkE7AQGA//I60DHtIJ/BR4zxnzgWLcK2AKMB3r6NnSlCsEY+wxr1y5o1creSvbis6Q9e/bwmuMW3YgRI7x2XqWCSmho5rPenMqVs5NZBECefbuNMasClGzBJsU1zmTriGcv8BOQX/HYnkAKdupA57GpwCdAVxEp4/1wlfISEXjkEahXz5ZxDA/32qkPHjxIz549SUxM5KGHHqJz585eO7dSKn/uxuGWBp4H+gINsD2TszLGGF8VzrgGWOhi/RbgAQ+O3WuMSXRxbGngSsf7bC7mUUotv97L7rYX5tiCnjsxMZGzZ88GXVwF2X7mzBmOHz8edHF5emxCQgJxcXGFP3fHjsjy5ZiQEEzOIRsFiCslJYX4+Hj27dvHqlWrmDt3LomJiVx99dVMmzbN7bFKKe9zlzBfxz77/Ar4Avs81V+qAqdcrD8JVHGx3tNjndtzKVeunMfBKeVtDwPbgF99/D333XcfM2fO5IorrvDxNymlcpK8fiWLSBww3Rgzyb8hgYgkA28aY57PsX4Strexu5b5t0CEMeamHOtvB74BOhljfsixzeUfQkhICKGhoYUaj5XfscXx3MXxmnx57nsvXWLm+fNcBG6qXJnDOYb/XO6ftYhQqVIlqlatSuPGjWnXrh2NGjVyG6dSJc28efOYP3++y23GGO8OxjXGuFyAM8CteW335QIcA2a6WD8dOJ7PsfOB7S7WP4jt7XyNi22mqFuwYEGgQyi0on4NlxX/F18YExpqDBgzbpz3gyqgov7fwBi9hmBR1K/BkRe8mtvcFcT8EujkZrsvbcE+i82pObmHKbk6NtIxtCjnscnArtyHKBUAX31lh/2kpcHzz8OLLwY6IqWUD7lLuNOAviLykoi0FZHGORcfxrUIaJ/1OxxVsDo6tuV3bCmydK4SkTCgD/CNMcafz6KVcm3hQjvTS0oKjBoFkyYFXV1ZpZR3ues09bPjNQZbfcoVXxVhnQWMABaKyBjsreAJwEEgo1iniDQEdgPjjTHjAYwxG0RkPjBFREoBe4GhQCTQz0fxKuW5+Hg71jY52c51+o9/aLJVqgRwl3AfwyY6vzPGXBCRW7GlHedgK18tw5Z2zDpdoGCTfs6W+kBgErY6VWVgI3CnMWadr2NXKl81atiqOOvXw4QJmmyVKiHyTLjGmNl+jMPV9x8Aeuezzz5czGpkjEkCnnYsSgWeMbYQ+5VX2s+9e9tFKVVi5DuLsIiEiEgLEfmLiHiv7I1SJUVqKowYYQu7r3U1vbRSqiRwm3BFZDhwFNiEnUygmWP9AhGJ8n14/iMiiIhOPq+86/x52zlq+nRIT4eDBwMdkVIqDzExMRm5wBfyTLgi8gQQCyzAjmHNGsEP5HO7t6hxjpPShKu8Ji7OTqu3ZImdN3TZMr2NrFQQi4mJyVqfwevctXCfxlZ7ehL4T45t23C0dpVSLqxeDTfcABs2QNOm8PPP0LFjoKNSSgWQu4QbCXydx7YL2N6/Sqmczp+3038dOQKdOtnk27RpoKNSSgWYu2FBJ4BGeWxrBuSeFkUpZefgfPddWLkSXn8dSpUKdERKqSCQX2nHl3JUlDIicgXwFPbZrlIKYNs26q5alfn53nthyhRNtkqpDO4S7hjslHybge+wRTCmAn8AacB4n0enVLAzBt57D66/njbTptliFkop5UKeCdcYkwC0BV7B1ibejb0F/S/gJmPMGb9EqFSwOn0aHnoIHn8cEhOJ69gRmjQJdFRKqSDl8hmuiJQGXgU+NsZMwNYxVko5/fADPPww7N9vn9lOn866ihWpX7FioCNTSgUply1cY0wyMBgo599wAkcLXyiPzZhhex/v3w/XXw/r1tnkq5Qq0gJW+AJYD1zrk28NQlr4Qnmsa1eoXBnGjNEhP0oVI74ufOFuWNAzwDwR2Q8sMb6KQKlgd/YsvP8+REVBSAhERsK+fVCpUqAjU0oVIe4S7mdAJWAhkCoi8WSfrs8YYxr6MjilAsoY+OQTePppOHoUypSBoUPtNk22SqkCcpdwlxGg+XCVCrg//oDhw2HFCvu5fXstzaiUKhR3CTcauGSMueivYJQKuHPnYNIkePNNO61etWrw6qswcKC9nayUUpfJ5b8gIhIGJAC3+zccpQLs449tgk1LgyefhO3bYdAgTbZKqUJz2cI1xqSKyDFsRSmlirf9+6GhozvCY4/ZMbYjR0K7doGNSylVrLj72T4XeNxfgQSajsMtgX79Fe64A669Fo4ds+tKlYK5czXZKlUC+XocrrtnuPuAv4rIr9ieykfI0YnKGPO+T6IKAB31VIJs2AAvvQRffmk/V6oEmzbB7foERamSLCYmJqPR5Yuk6y7hvuV4rQtc72K7AYpNwlUlwJYtMHYsfP65/Vy+PIwYAX//u+0cpZRSPuQu4Ub6LQql/CEqCpYvh7Jl7Xja0aOhZs1AR6WUKiHyTLjGmP3+DEQprzLGTgBfvTq0aGHXjR8PzZvD889DnToBDU8pVfLoWAdVvKSnwxdf2EIVt95qn9U6dewI06ZpslVKBUSeLVwR2Us+laaMMY29HpFSl+PMGZg9G/71L9i1y66rVg3atLGtXR/1OlRKKU+5e4a7itwJtxrQATgPLPdVUEoVyJdfQt++cOGC/dywITzzjB1TGx4e2NiUUsrB3TPcR12tF5HKwH+B73wUk1LupaTY2Xqc0+K1bg2XLkHnzrZgRY8eEObut6RSSvlfgZ/hGmNOA68DL+W3b1GihS+KgG3b4NlnoV49uOsu+7wW7Oe9e20P5F69NNkqpS5LIAtfuHMRqOfNQAJNC18EqfPn4dNP4b337GTvTldcYafMc3aAqles/joqpQIgkIUvcnFMatACiAG2eD0apbLassX2Nj5/3n6OiICHHrKTCbRrpx2hlFJFSp63lEUkXUTSsi7AJeA34ErgKV8FJSIhIvK8iOwTkYsislFEent47GwRMS6WKb6KV3mBMVTZtg1mzsxc16yZrQbVoYNt4R45ArNm2SSsyVYpVcS4a+GOJ3cv5YvAfuArY8wZn0UFE4C/AS9gE/xDwGci0t0Ys9SD448DPXOsO+LdEFWhGQO//WZLLc6bR6f9++3kAQ8+CFWq2Gexf/wBVasGOlKllCo0d72UY/wYRwYRqYFNtpONMW84Vq8QkSuByYAnCTfZGLPGVzGqQjp0CF57DRYsgIMHM1YnVatGuYEDJnKnsQAAGaNJREFUbS9kJ022SqliIlvCFZEQoBuw1xiz2dUBInIt0MgY86WPYuoKlMZOD5jVXOB9EYk0xuz10XcrXzh/3s45e8019nNIiK34BLbT0733Qp8+fJOQwD29egUuTqWU8qGcz3D7A/OAC26OOQfME5G+PorpGuyz4l051js7aTX34Bw1ROSEiKSKyA4RGS0ioV6NUuXNGNvh6Y03oEsX20rt3t2uB5tkp02DNWtsC/ett6BTJ5uIlVKqmJKsw2FE5BtguzFmpNuDRGKBZsaYO70ekMg7QE9jTK0c668EdgIDjDFz3Bw/CkjDJuiyQC9gEPC+MebxPI5xOSaoT58+9O3rq98VxU+FAwdo/OWX1Fi/nvInTmSsNyEhnGralDUvvEBKxYoBjFAppbKbN28e8+fPd7nNGOPd3pnGmIwFOAH0yLrO1QL0AE7kt59j39uwna/yW1Y69p8FHHFxnqaO/R725HtzHPtPx7FN89huiroFCxb49wtPnjRmwQJjli/PXLdihTG2HWtMzZrGPPKIMfPmGXPihEen9Ps1eFlRj98YvYZgodcQeI68UKBck9+Ss9NUBeCUB3n6lGNfT6wGrvZgv0TH60mgioiI46KdqmTZXlDzgFFAW2wrWRXU8eP2FvCqVbBiBaxfb1Pr3Xfbkopgh+tMmgR33gmtWuktYqWUyiJnwj0BNAR+zOe4Bo5982WMSQS2FSCmLUAZoAnZn+M6n91uLcC5nJy3BbSclCdSU20yLVXKfo6Kyuzk5FSqlC0+0alT5rqyZeH//s9/cSqlVBGSM+H+CDwCfJTPcY+Sf1K+XP8FkoF+wLgs6/sDm83l9VD+KzbZ/lr48IqZ9HQ7nd26dbbVunYt/PorzJ8P3brZfZo0sQUobrjBzinbubMtRlG+fGBjV0qpIiRnwp0C/Cgi/wRGG2OSs24UkVLAG8CtwJ99EZAxJt7x/c+LyDlgHdDH8Z335IhnGdDQGHOl43NDYA7wCbZ1XAbbaepRYKYxZrcvYi4y0tIg1NFZOzXV9iBety6zdGJWW7dmJtwnnoDhw3VSAKWUKoRs/4IaY34WkWeAN4F+jl7L+x2bGwK3Y+fEfcb4trDEC9g5d6OBWsB24EGTe+xvKNmv4Rz2Ge9ooCa2VfsHEAVM92G8wSU5GXbutENztm61r1u22MR64IDdJyzMlko8fx7q1rVT3LVpY5f27aFmzczzaUtWKaUKLVeTxRgzRUTWAc9hW4flHJuSgJXYClA/+DIoY0waMNGxuNvvlhyfTwL3+i6yIJKaapPnrl3QqFHm+g8/hMcft9tzEoHTp6FyZfv5s8+gdm2oUcMvISulVEnm8h6hMeZ74HtH5akrHKsTHIlQ+UNqavZbuLGxsGMH7NkDu3fb+V+dSTUmxvYKBltUIi2N/9/emcfbUVR5/PvLCiECQY0BgiQho5CAgEMYdZBkBCcOigkMChiWKKigoAEZB4GRF4ggKq44gKIGSExAWR5gQFkUEIzAuIBB4ggEDDtJCDhZJWf+qGpfp9P3vXvfXbrvy/l+PvXpd6tPdZ3q6nqn63QtjBkTVnYaPx7GjQvHXXfduLe6554tK47jOM7mTrcf5cxsA/B8i3QplGTvw7PPPrt5m9CbBWOYGNLFi8MUm+efD+Hpp8PKS0uXwurVoTeacMklYQP2NCNHhgFN6b1gJ04MbmJ3AzuO49RER0cHM2fO7Fmwl/gomMhGU37NwuhdqWsu6cqVsGxZMISrV8OaNV1/Dx4Mkyd3pT/5ZFi+PKRZuTIYzhUrglH95jfhhBOC3D33wIkn5ivUrx+88gq8Jk53njEj5DV6NIwdG3qwW27ZJd/ZGY6DBoXgOI7j1ESpNqDv0wweHHqfGzZ0rfl7661w4IHh73PPhQsvzE+7yy7hW2rClVcGQ5vH8tS6HXvsEUYADx8ewogRsNNOIYwYsbFL+eMf733ZHMdxnMJxg5uwbt3Gv6VgfBOGDw+Dk7bcMoQttuj6e/vtN0779a+H6TfbbBMGKCXH4cM37pVOmBCC4ziO0+dxg5uwZk1w4/bvH4xt1p3w2c+GUA3TpzdcPcdxHKe9cYObMHhw0Ro4juM4fRhfXd5xHMdxWoAbXMdxHMdpAW5wI5KQ1Lw5uI7jOE6p6ejo+LstaAb+DTey8da7juM4zuZGs+fheg/XcRzHcVqAG1zHcRzHaQFucB3HcRynBbjBdRzHcZwW4AbXcRzHcVqAG1zHcRzHaQFucB3HcRynBbjBjfjCF47jOJs3vvBFi/CFLxzHcTZvfOELx3Ecx+kDuMF1HMdxnBbgBtdxHMdxWoAbXMdxHMdpAW5wHcdxHKcFuMF1HMdxnBbgBjfi83Adx3E2b3webovwebiO4zibNz4P13Ecx3H6AKU0uJJOlXSjpGckmaSOGtNPlfRbSWskPSHpLEn9m6Su4ziO4/RIKQ0u8FFgOHB9rQklTQauAe4H/g34BnAWcF4jFXQcx3GcWijrN9zxZrZB0gDghBrTfhH4pZl9LP7+uaShwFmSvmZmzzZUU8dxHMepglL2cM1sQ2/SSdoJ2AuYkzl1JTCQ0ON1HMdxnJZTSoNbB+Pj8Q/pSDN7HFgFjGu5Ro7jOI5D3zO428XjipxzK1LnHcdxHKelqNnzTyUdCNxaheidZjYpk3YAsB6YaWYdVeQ1jeBO3tXMFmfOPQXcYmbH5aTLvQmHH344Rx55ZBWqO47jOO3IvHnzuOqqq3LPmVlDJ+O2YtDUvcBuVcitakBey+Mxrye7ber8JrT7whednZ1MmTKlaDXqot3L0O76g5ehLHgZWseUKVOYP3/+JvHNWPii6QbXzFYBjzQ7n8iieBwP/CqJlDQKGAI83CI9HMdxHGcj+tQ3XDN7Evg9MC1z6iiCa/rmlivlOI7jOJR0Hq6kfYBRdL0QjJN0WPx7Qew1I+l2YGczG5tKfgZwk6RLgXnA3oSFL77hc3Adx3GcoiilwQVOAo5N/f5ADACjgSXx7/5kymBmC6JxPhuYDjxHWGXqC81T13Ecx3G6p5QG18ymE4xlT3KTKsRfC1zbUKUcx3Ecpw761Ddcx3EcxykrbnAjvgG94zjO5o1vQN8i2n0eruM4jlMfvgG94ziO4/QB3OD2EebNm1e0CnXT7mVod/3By1AWvAx9k6avpdwOSLJ2vw+S2t4t3u5laHf9wctQFrwMxRP1b6hf2Xu4juM4jtMC3ODWQDUjmKsd5dxouSLy9DL0jiJ0a/To+yJ08zL07nreFnov12jcpUz1LuVqXCTVulGKkCuzbtXKuW6uWxnlXLc+q5u7lB3HcRyn3fAeLpU3oHccx3E2Xxrdw3WD6ziO4zgtwF3KjuM4jtMC3OA6juM4TgvoEwZX0k6SfixppaSXJV0r6Y1Vpt1C0pclPSNptaRfSdo/R66fpM9JWiJpjaTfS/r3ossgaR9J35H0iKRVkp6UNFfS6BzZJZIsJ0wtsgwxbZ5eJmmvjFxZ66GjmzKsycg2rR4kjZT0rfgcr4rXHVVl2rK0hV6VoWRtoZ56KEtb6G09lKUtHCbpGklPxOd5saTzJb2mirTNaQtm1tYBGAL8L/AHYCowBXgIeBTYqor0c4GXgI8CBxD20V0N7JWR+wKwFjgN+BfgUmADcFCRZQC+AtwDfAKYCHwI+COwDNgpI7sEuAV4WyYMK0E9GPCDHN2GtEk9jMzR/QBgPXB1C+thEvAcsAD4abyvo6pMW3hbqKcMZWkLDaiHwttCnfVQlrawELgamBafhxnx+V4I9CuiLdRdKUUH4NPAq8DYVNxo4G/AqT2k3TM+RB9OxQ0AFgM3pOKGx5s6M5P+duDBgsvw+py4nWOln5PzcM8pWz1EWQNm9SBT2nqocL2jY7ne28J66Jf6+/ga/kmWoi3UWYZStIV6yhDlC28L9ZYh51pFtIW85+GYqMe7uknXtLbQF1zK7wcWmtmfkwgze5zwpjulirTrgatSaf8GzAcmSxocoycDg4A5mfRzgD3yXFY10usymNkLOXFPAC8AO9apVy3UUw/VUtp6qMCxhB7CT+vUq2rMbEMvk5alLfS6DCVqC/XUQ7WUth4qUERb2OR5AO6Px+6eh6a1hb5gcMcTXIBZFgHjqkj7uJmtykk7CBibklsL/DlHjiry6Yl6yrAJknYjvH39Mef0wfF7zFpJCxv1zYrGlOHEqNcqSXdIemdOHm1RD5JGElxMc2NjzdKseugtZWkLDaWgttAIim4LDaNkbWFiPOY9DwlNawt9weBuB6zIiV8ODKsjbXI+Ob5k0V/QjVxvqacMGyFpAHAJ4a3+e5nTNwInE97MpgFrgOskHVWrwjnUW4Y5hG9vBwIfA14L3CFpUiaPtqgHggutH3B5zrlm1kNvKUtbaBgFtoV6KUNbaCSlaAuSdgTOAW4zswe6EW1aWxhQvbqlJm/1jmpWCFGVaauVq4dGXf8i4B2EbyUbPTRmdvJGF5euIwwgOJ9N3SK9oddlMLOjUz/vltRJ6G3OAvZLXatd6uEY4Ldm9uAmGTS/HnpDmdpCoyiyLfSaErWFRlF4W5A0FOgkjMf4cE/iNKkt9IUe7gry3+aGkf+WkmZ5N2mT88lxmKTsjczK9ZZ6yvB3JJ1PeCP+iJn9rCd5M3sV+BEwUtL21eZTgYaUIaXbK8BPgAmp6Haph32BXcl/o9+EBtdDbylLW2gIBbeFhlJQW2gIZWgLkrYAbgDGAJPNbGkPSZrWFvqCwV1E8KVnGQc8XEXa0ZKG5KRdR5dvfhEwGNglR44q8umJesoAgKQzgdOBT5vZlTXknTws9a7xWXcZcsi+QZa+HiLHEt6kf1hDmkbVQ28pS1uomxK0hWbQ6rbQKAptC5IGAtcA+xKm6jxURbLmtYVmDMduZSDMrfobMCYVN4owyuwzPaTdi1Cpx6biBhA+qN+YikuGf5+dSX8b8FCRZYiyn4rlOKPGfAcADwBPFF2GnOttDTwJ3Nku9RDlBxHmfXYWUQ+Z69YypaYUbaGeMkT5wttCvWXISd/yttCIMhTdFggdyqsJ34UPqCFd09pCQyuliABsRXjjeIgwdeP9wO+Bx4ChKbmd4z/Tz2fSzye4C48nTHD+caygt2bkvhjjTyVMCL+YML/v4CLLABwR9biZTSePj0vJHRnLegxhxOARwN3xwTqi4DKcBnyXsFDBJMJb8UOEt8l3tkM9pM4dGu/poRXyaGo9xDwOi+HieN0T4++JZW8L9ZShLG2hzjKUoi3U+yyVoS2kdJ6V8zyMLKItNKRSig7AGwlug5eBV4DrybyJEXoqBnRk4rcEvgo8G2/cr4FJOXn0B84CniC81TwIHFZ0GYDZMS4v/CIl9zbgDsJcuPXASsKb2OQSlOFgwlzXF6NuywjfXPZtl3pIneuM+g+qcP1W1EO3z0PZ20Jvy1CmtlBHGUrTFup5lsrQFgiLalTSv6OItuDb8zmO4zhOC+gLg6Ycx3Ecp/S4wXUcx3GcFuAG13Ecx3FagBtcx3Ecx2kBbnAdx3EcpwW4wXUcx3GcFuAG1ykESdMlWSq8KukpSVdLenNGtkNS0+evSRoo6QVJC7qROSDqO73Ka06K8pNScTMkHVq/xs1H0mckPZizXmyfQNKobH1Kmi1pScE6fEPST1qlg9Ma3OA6RfMB4O3A/sDngL2B2yVtk5K5LMo0FTNbD8wD/lXSGyqIHQOsIqw801tmEFbhKTWStgXOAM6xzWvC/rnAIQXr8EXgXZLeVbAeTgNxg+sUze/MbKGZ3WNmVxCWjtuRsK0aAGa21MwWtkifywmrx3woe0LSVgRDeY2Z/bVF+hTJcYQVgK5rZabR01BYj9rMHjWz3xaVf9ThGcJesacVqYfTWNzgOmXj5XgcmETkuZQlnSTpV5KWS3pJ0kJJ783IDJB0rqRHJa2R9KKkX0rajwqY2f8QdgE5Ouf0ocBQ4Ip4/a0lXSTpaUlrJS2WdEp3xiK6KncGpqXc6bPjubGSrpT0uKTVkh6TdLGkYTnX+bSkJbFc90l6R/w9OyM3WtLc6CpfK+l3kqrtvR0PXGVhy7Tkeon78+OSzpH0TLz/N0oamcl7oKRZUa918Tgr7uCSvd4nJH1J0tOEJfK2ja7dpZL2kXRvvCeLk3qWdGq85suSOiW9PpN/j89IHlmXcszD8kIm3aExj1Uxvx9JemNGZoik/5a0TNJfJd0AbHTfUswHJkvaqSednfagr2xA77Qv/SUNIPQqxwDnAc8Dv+gh3SiCq3kJ4Tk+GLhJ0kFmdnOU+U/gFOBM4HeEXVf2IX+vyzSXA1+SNN7MFqXijwaWAndI6kfYo/StwOcJC8y/l7D+6usJrtg8DgEWEDZF6IhxL8TjDvH6MwgLp4+J11lAyqUu6Xjg68D3CHuH7kLY/mzbdEbxH/WvCffzlJjP4cA1kqaa2Q2VbkA0FLsC/1VB5HPAvcBHCLumXAjMBSamZC4HPkio01/GMpwVy5X1IJwJ3E/Yw7Y/Yf1aCHV2BfAV4Okod42kbwNvAj4JvCHej2/H/BJG0fMzUg2HELZhSxhKuN/LkghJJxAWrv8BcA7wGkL93inpLRb2tAW4lFAHM2N5303lrevuInSK3g18vwZ9nbLSyIWuPXioNgDTyV9U/ClgQka2IzyqFa/Vj/AP9WektgIDbgKu7YVu2xN2ELkgFbcD8CpwXvz9vqjv9Ezaywg9tNfF35Oi3KSUzBJgThV6DAD2i+n3TpX1L8CCjGyyM8vsVNz3CEb2tRnZWwmu/O7yPjxe7x8y8aNi/J2Z+NNi/A7x9+7kLwp/Vox/S+Z6v4GwtntKdnY8t38q7i0xbjHQPxX/VYL7u3+F8lR6RpL8p2fyXdLNdTrjfR0T44YSFt7/fs69WgfMiL/fHJ+h0zNyF2d1SJ37C/CdVrZND80L7lJ2iuYQYAJhg+iphE2bF0jarbtEkv5R0k2SniMYx/WEnkB6hPP9wEGSviBpP0mDqlHIwvezWwlu36SNHEX4Z3tF/L0/YRuueZnkcwj7gNY8yEvSIElnSHpE0upYprvj6aRcI2P4USZ5J+E+pHkPoXe8MrrXB0Rvwk+BPSVt3Y06O8TjCxXOZ0fQJht7Jy7U/eNxTkYu+T0xE3+9RQuT4f/M7K7U70fi8TZLubpj/ADCyxJQ9TNSKxcAk4GpZvZYjHs7oSc+N3Ofl0a9knvxT3Tt0Zpmfjf5vUBXXThtjhtcp2j+YGYPmNn9ZtZJ2INWdLlbNyG6Sm8nuIZPJgywmgDcAmyREj0PODte825gmaQfSHpdFXpdThi8lYwSPRq4z8ySf/jbAcvNbG0m3bOp87VyPqHccwju6X3pGs2clCsxKM+nE0bj82LmesMJo6rXZ8KX4/nXdqNLkl+2fAnLM78TuSRdUv5nMnKV7k9WLuGl9A8zWxf/XJGRS+K3gJqekaqRdByhJ3+cmd2TOjU8Hm9j03u9B133Oam75zKXzv5Os5qwVZzTB/BvuE6pMLPVkh4juA4r8R5gG+CDZrY0iZQ0JHOt9YQeyQWSRhDcwF8FhhBcpt1xPcFNeLSkZQQX6SdT55cD20kalDICACPicRm1cwRwhZnNSiIkDc3IJIZpeDpSUn8g+yKxjPCicUGF/J7uRpdE/2GEf/q1khjkEcCjqfhK96fR046qekaqRdJEgut3ppnNzZxOyjKdMOAuS/L9Nqm7NwCPpc5XmoIG4YXhwZqUdUqL93CdUhH/Ie5CZVcmBIMJoQeRpHsT8M+VEpjZs2Z2GaEXsntPepjZGoLr71DgBEIPKu36u5PQfj6QSTotynY3jWkt+b2WIaTKFPlw5vfSGLL5TmXTF+hbCC8ui6IXIRsq9V6hy3U7phuZ7rgzHo/IxE+Lx7toLjU/I5WQNBa4FvixmXXkiNxLMKpjK9znxVHu14TPEB/MpM/eoyTf/sBOhO/VTh/Ae7hO0ewVXbwiuNxOIrzVf6ubNLcRvsldIenCmG4m8CSpl0hJnYTRwL8huCD3JvR8Lq1St8uBj8ZwnZml3ag3E0beXhKnoywCDiJMpTnfzLLu3TQPA++U9D6Ci/VFM1tCMJDHSnoI+DPB2L8jndDMNkiaCXxX0mWEb7ljgNMJPfINKfHPA/cBd0m6iDBYaxjhhWOMmX2kGx3vI7wY7BvLWRNmtkjSPKAjfs+8l/Ct87+AeWbW7F5bVc9IldxEWOzkUklvS5+wMIf8ZUn/AXw7Pgs3E+piR8K36l+Y2Q/NbLGkHwLnxLEBySjlgyrkuzuwFc1/OXFaRdGjtjxsnoH8UcrPA3cAkzOyHWRGKRN6CY8Qpo8sIvQSZpMaXQp8htDTXEZwiy6O1xpYg55/irpNyTm3NXARwVW4LsqeQmq0LfmjlHcluHpXkRpZTHAJzye8HKwgTLOZQP5o6BnAE7H8DxBGM68AvpaRG0kYOf1U1DEZEHZUFWW/Cvh5Jm5U1Of4THxeOQcCs6Ke6+NxVvr+V7pePDcbWJoTb8CsCs/T2BqfkST/6Zl8l2Tyyw0ZHQ4Cfk6YS76a8NL0fWBcSmYIwTW9HPgrcAOh151Xx2fG+hpQdHv10JigWLGO47QxkiYQeqXHmNmVDbrmJMIL0Cgze7IR13SqR9LDhFXNKs2FdtoMN7iO02ZIGk0YwHU3oTe1G2GBjHXA7ma2qoF5/Qz4k5md1KhrOj0jaQqhd7yLmb3Uk7zTHvg3XMdpP1YTvu8dQ/gmu4LwzfL0RhrbyKeAqZJk/nbeSrYkuP3d2PYhvIfrOI7jOC3ApwU5juM4Tgtwg+s4juM4LeD/AX8WeAnJwAoKAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "resp = qmix.respfn.RespFnPolynomial(p_order=50)\n", "\n", "# Plot\n", "fig, ax = plt.subplots(figsize=(7,5))\n", "resp.plot(ax=ax);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Define embedding circuit\n", "\n", "The ``EmbeddingCircuit`` class contains all of the information about the embedding circuit. This includes the Thevenin voltages, the Thevening impedances, and the frequencies of the applied signals. In this notebook, the Thevenin impedance is ignored. This means that the AC voltage across the junction is equal to the Thevenin voltage of the embedding circuit (no harmonic balance required)." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - Note:\n", "\n", " - All of the circuit properties are normalized: \n", " \n", " - voltages are normalized to the gap voltage (``vgap``), \n", " - resistances to the normal resistance (``rn``), \n", " - currents to the gap current (``igap = vgap / rn``), and \n", " - frequencies to the gap frequency (``fgap``)." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Running harmonic balance:\n", "Done.\n", "\n" ] } ], "source": [ "# Create an instance of the embedding circuit class\n", "cct = qmix.circuit.EmbeddingCircuit()\n", "\n", "# Set the normalized photon voltage\n", "# (equivalent to the normalized frequency)\n", "cct.freq[1] = 0.33\n", "\n", "# Set the Thevenin voltage\n", "cct.vt[1, 1] = 0.33\n", "\n", "# Calculate AC voltage across junction\n", "vj = qmix.harmonic_balance.harmonic_balance(cct, resp)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Calculate desired tunnelling currents\n", "\n", "``qmix.qtcurrent.qtcurrent`` will calculate the quasiparticle tunnelling current based on the voltage across the junction ``vj``, the simulation parameters contained within ``cct``, and the response function ``resp``.\n", "\n", " - Note: \n", "\n", " - The frequency that it solves for is the last argument passed to the function below. For example, \n", " \n", " - ``0`` corresponds to the DC quasiparticle tunnelling current\n", " - ``cct.freq[1]`` corresponds to the AC quasiparticle tunnelling current of the first tone and first harmonic" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Calculating tunneling current...\n", " - 1 tone(s)\n", " - 1 harmonic(s)\n", "Done.\n", "Time: 0.0041 s\n", "\n", "Calculating tunneling current...\n", " - 1 tone(s)\n", " - 1 harmonic(s)\n", "Done.\n", "Time: 0.0030 s\n", "\n" ] } ], "source": [ "# DC tunneling current\n", "idc = qmix.qtcurrent.qtcurrent(vj, cct, resp, 0)\n", "\n", "# AC tunneling current\n", "iac = qmix.qtcurrent.qtcurrent(vj, cct, resp, cct.freq[1])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Plot the pumped I-V curve" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdkAAAFCCAYAAAC95qeCAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOydd3gUVffHPzeEFHpASqihF+miiI2gSFMpthcFBRUU7PUFAWVFeMVe8OeLgCJFpQuhSJWmJLwCRkKvoUMILYH05Pz+mN2YstkUdnaymfk8zzxhZ+6dOd/DnT17Z+69R4kIFhYWFhYWFu7Hx2gDLCwsLCwsSipWkLWwsLCwsNAJK8haWFhYWFjohO5BVin1sFJqoVLqmFIqUSm1Xyn1gVKqfAHqBiilPlZKnbHXDVdK3eWknI9S6m2lVLRSKkkp9bdS6iF9FFlYWFhYWBQMT/Rk3wTSgVFAD+C/wHBgjVIqv+t/BwwF3gXuB84Aq5RSbXOUex+wAV8DPYEIYL5SqpebNFhYWFhYWBQapffoYqVUVRE5n2Pfk8AM4B4R+S2Pem2ASOBpEZlu3+cL7Ab2i0hv+75qwAlgooiMzVJ/HVBVRFrrIMvCwsLCwiJfdO/J5gywdv60/63lompvIBWYm+VcacAcoLtSyt++uzvgB8zOUX820EopVb8odltYWFhYWFwvRg186mz/u9dFmRuBoyKSkGP/brSg2ihLuWTgkJNyAC2uw04LCwsLC4si4+vpCyqlagHjgLUiss1F0crAJSf7L2Y57vh7WXI/985ZLqcd1iocFhYWFhbZEBHlzvN5tCerlCoHLAHSgKfyKw44C4Q5HVDQcrkQEdavX4+ION0cZfLbjCjnrnM59BdH2zxRLq///+Jgm6euGRoaWmxt84TfcrYBb9RwveVcfQ8abZtefkuMu8aGBwaRjkKAfcGOh6PuxWNBVikVAIQBDYDuInIynyoXcd4LDcpy3PE3SCmVM6jmLOeUm266KR8zSjaWfnPrB9i+fbvRJhiK1QbM54ODa/7geKOWdF46A1Gwse/TNDgUpcu1PBJklVKlgYXALUAvESmImt1AfaVUmRz7WwAp/PMOdjfgDzR0Ug5gj6uLuPqCGTt2bJ7HikM5d5zLod+dtrnT/sKcryjl8vr/L84a3O2Pgn7BGmGbJ+6FnG3AGzVcb7n8fmh5g4aClElPS2Pj0BHU6dmFJjFHOVGxBgfmLkNeeYLSZQIKdL1CU5Du+PVsaIF8HpCENmWnoPXaoj0GHpRlny/aYKmlWfZVQxv4NDZH/bVAlIvzi4hIRESEeCMO+68XI/W7S8P1cL36i4OG6+WWW24x2oTr5nr+H4rLd4CRbckdPiju98KJv/bKX/XbioAIyJa7+sjV2Esi8o9+uwa3xkBPDHz6P+ARYAJwTSl1a5ZjJ0XkpFKqHnAYGCci4+wRMFIpNRf4wt4TPoq2iEV9YIDjBCISo5T6HHhbKRUP7AD+BdwN9MnPuJCQEDdI9F4s/SFGm2A40dHRRptgKFYbKNk+kIwMtti+ouWHo6mdksCFwAoc/+BLOr0yOLOMnvo98bi4p/3vaCA8xzbEfkwBpZzY8xQwHRgPLAfqAD1EZEeOcqPtZV4BVgG3A4+KyNL8jNu719UsouKLux7fGKnf3Y9ki8L16i8OGq6X4cOHG23CdXM9/w/F5TvAyLbkDh8Ux3vhwrHTRHS4l9vff42KKQlsb3kb7IqiXZYAC/q2Ad1XfCquOKbwhISEMGjQIGw2m8EWGUN0dHSJ/hWbH2bXD5YPzK4fSqYPtk3+mVpvvkjwtYtc8/Vn1+tjueWDESif7H05m83GjBkzMp/oiJun8Jg6yIoIUVFRtGrVymmZuLg4YmJiSE1N9bB1niMlJQU/Pz+jzTAMPfT7+voSEBBA1apVCQjQaTCFG3F1D5gBs+uHkuWDa5fj2fbIEDqvnQfAnjotqLjgZ2rdkvcKuw79Sim3B1mPL0ZR3Lhw4YLT/XFxcZw7d45atWoRGBhI7hlCJYP4+HjKl883IVKJxd36RYS0tDSuXr3K8ePHqV69OhUrVnTb+fUgr3vALJhdP5QcH+xeso6Ap5+i88UTpKpS/G/gC9w67WNK5fNDWk/9pu/J5vUle+jQIWrWrEmZMjlnEJUs0tPTKVWqlNFmGIae+hMSEjh79iwNGjTQ5fzuwvqhZW794P0+SE1K4Y9n3uT2n7+htKRztEpt0mbMpPF9XQpU36Ffj56s6ZO25zU/LDU1lcDAQA9b43kSEnIuDW0u9NQfGBhIcnKybud3F2ZfjMLs+sG7fXA0PJJ9TdoR+tMkSks6W7r3J/jw3gIHWNBXv+mDbLly5fI8VlIfEWfFx8fcTUBP/d7SflzdA2bA7PrBO30gGRlsemMCVTt3otWJPZwrW5moqXO4beXPBFQsnB499Zv+nWzNmjWNNsFQSpcubbQJhmJ2/WDdA2bXD97ng5iDx4ju+zh37dkCwJ833U3TxT/SqnaNIp1PT/2m7sYopRgwYIBpp+8AXvE4U0/Mrh/gwIEDRptgKGbXD97lg4jPvke1acMte7Zwxa8Mf773BTf/uZYKRQywNpuNAQMG6PbkydRBVkSYOXOmKYLs4MGDqV27dq79/v7+bNiwAaUUa9euNcAyY/H39y90ncGDB5eoOYUNG+Zc9ttcmF0/eIcPrsRcZNOdfbj1jWeomniFqIZtSdq2g5vffQWuI0DabDZmzpyJXoOATR1koeQMXS8qaWlpRptgKGbXD9Y9YHb9UPx98PfPy7nS5Ebu+j2MpFKlCR/2Ni33b6N6q6ZuOb+e+k3/Tvby5ctGm2Ao6enpRptgKGbXD9Y9YHb9UHx9kJyQxJYBL3HXku8oJcLB6g0o/fOPdOpya/6VC4Ge+k3fkzVbHsWcOJsHHBISwsCBA5kzZw7NmzenbNmydOjQgd9//z1bOccj6C1btnDzzTcTEBBASEgIkyZNylbOZrM5fd+R87FrdHQ0SikmT57M22+/TY0aNShfvjwDBw4kISGBQ4cO0b17d8qVK0ejRo2YMWOG0+tERUXRpUsXypQpQ3BwMO+++y4ZGRnZysbGxjJ8+HCaNWuGv78/zZo1Y8qUKblsXLduHe3btycgIICGDRvy7bff5utTb8Ps94DZ9UPx9MGB37ZypGEruiyehhLY0ncwIYd3EeLmAAv66jd9kPXm+WHuIK95ops3b+bTTz/l/fffZ+7cuaSnp3P//ffn+sUXFxfHv/71LwYNGsTixYsJDQ3l5Zdf5ocffiiyTR988AGnT59mxowZjBs3jrlz5zJs2DD69evHfffdxy+//ELr1q156qmn2L17d676ffv2pWvXrixevJjHH3+c999/n3HjxmWz+fbbb2f58uWMHDmS5cuX88ADDzB8+PBsPxD27t1Lr169CAwMZM6cOfznP//hiy++YN26dUXWVhwx+z1gdv1QvHyQnpbG+uFjqNPtTpqfPcTpCtXY/+NibvtlOqXL6rN2gZ76Tf24WClF27Zt6dOnT4EGP5V/8if9jSoA8TMfd9u58lrtKC4ujsjISIKCggCoUaMGN998MytWrODxx/+5fnx8PFOmTKF///4A9OjRg1OnTjF27FgGDRpUpBF7DRs2zOyldu/enc2bNzNr1ixmzZrFwIEDAejQoQNhYWEsWLCAG2+8MVv9oUOHMnLkSAC6detGXFwcn376Ka+++iqVKlXiyy+/5NixY0RFRVGrVi3KlClD165duXz5Mu+99x7Dhw/H19eX8ePHU758eVavXk3ZsmUBuO2222jYsKHXTXlwRaVKlYw2wVDMrh+Kjw9ORR3k7IOP0eWQFvS23nEfLRfNpGbVyrpd02azsWTJEiIjI3U5v6l7siJCWFiYKUYX54Wvr/PfWZ06dcoMsEDm4uHHjx/PVq5UqVI89NBD2fb179+f48ePc+rUqSLZ1LNnz2yfmzVrBmgB10FQUBDVqlXjxIkTueo/+uijuey5evUqu3btAmDlypV07NiR+vXrA9rgp7S0NLp3786FCxfYs2cPAOHh4fTq1SszwALUqVOH22+/vUi6iitVqlQx2gRDMbt+MN4HIsLm97+mzM3tuOnQdi4GlOevT6bQcfMyyuoYYEELsmFhYbqNLjZ1Txbg8OHD1KlTp0Bl3dmD9DS+vr5OB/kkJydn7s8acCtXzt6wHVNdkpKSsu0PCgrKtaBD9erVATh16pTTaUP5kTW4A5lZcpztz2lP1us7swcgJiaGQ4cO5bkQhWOk4ZkzZ3Kdy3G+o0ePFkSKV1CYe6AkYnb9YKwPLpw4y94HB3HnttUARDbvSJ3FP9OuSX2P2aCnftMH2SZNmhhtgkeoVq0asbGxuVK7+fv7c/r0aSB3cCoIly5dIjU1NVvAOnfuHAC1atUCyEz3lvPaeg2bP3fuXLZF+XPaU6VKFapVq8aXX36Zy3aApk21aQHBwcGZdXOevyRhlnsgL8yuH4zzwbbvFlDz1eHccTWWxFJ+7Hz1HW75aFSunK96o6d+Uz8uBjIDTEmnS5cupKWlERYWlm1/amoqCxcuJDg4ODO4FIb09HQWLlyYbd+cOXOoW7duZlCrV68eQObjWtCGzG/ZsqXQ1ysI8+bNy2VPuXLlaNmyJaC9N963bx9169aldevWdOjQIdvmyEbSqVMnVqxYwbVr1zLPdeLECf744w9d7DYKs9wDeWF2/eB5H1y7cpX1PQfSfsij1Lway/5aTbnwewQdPxnj8QAL+uo3fU/26tWrRpvgEbp27cq9997L4MGD2bdvHx07diQ+Pp7Zs2ezZMkSpk+fXqTF8suXL8+///1vYmNjady4MT///DNr167lhx9+yBz01LNnTypWrMjQoUN57733SE5O5qOPPtJtUe6pU6eSkZHBzTffzKpVq5g2bRo2my1zcMdrr73G3LlzufPOOxk+fDht2rTh2rVr7Nu3j82bN7NkyRIAxowZw/z58+nWrRtvvfUWKSkpjB07tkg9/uKMWe6BvDC7fvCsD3Yv30jpwYPoEnuMNOXD1see55bvP6VUEVZfcxe66hcRU26adJG4uDhxxp49e5zu92YSExNl9OjR0rhxY/Hz85Ny5crJ7bffLosXL85Wrl69ejJgwIBc9QEZO3Zs5udBgwZJrVq15I8//pAOHTqIv7+/1K1bV7788stcdTdv3iwdOnSQwMBAady4scyaNUsGDRok9erVyyxz9OhRAWTq1KnZ6o4dO1YASU1NdWmno1xUVJSEhoZKQECAVK9eXcaMGSPp6enZ6l68eFFeffVVCQkJkdKlS0vVqlXljjvukM8//zxbuTVr1kjbtm3Fz89P6tevL5MnT85ltyu8oR3ldQ+YBbPrF/GMD1KTU2Tdk69Jko+vCMixoJpycMka3a9bEBz67XHBvbHG3Sf0lg0QQEJDQ7MFDgfe8OXoDq7n5nIE2eJCXsHYFXp/uXhDO1q/fr3RJhiK2fWL6O+DI1t3SmS9VlrIAQnv+rAkXryi6zULytixYyU0NFQcMUHcHGtM/bhYRIiKisqcnmJG8ponaxbMrh+Mn75hNGbXD/r5QDIy2DTqE9p/ZqN+aiLnywRx9rNJ3PrcAF2uVxRsNltmHNAjE4/ub5iVUrWVUpOUUuFKqQSllCilQgpQz2Yv62xLylE2Oo9yffO7jmOQi1kxe5Axu36w7gGz6wd9fHDu8Aki2nah84cjKJ+ayI62d+G3ZxetilGAdaBnG/DEMK5GwKPAJWBzIepNAzrl2LoCaUCYk/KrnJTfmN9FoqOjC2FSySMlJaXIdX/44QdOnjzpRmuuD5vNpj2eyWOBDWdcj/6SgtnvAbPrB/f7IPyrWajWrekUtYn40oFse+cT2u/YQMV6xXOlND3bgCceF28SkeoASqkhQLeCVBKRk0C2b3Cl1BNoNs9wUiVWRCIKa1zz5s0LW6VE4ZjDalbMrh+se8Ds+sF9PrgSe5m/Hn6G0I2LANhdvxU3LJxDh3Yt3HJ+vdCzDejekxWRjPxLFZhBwDm0XqtbMPuv2OTkZKNNMBSz6wfrHjC7fnCPDyIXrOZCk5aEblxEik8ptj7zBi0O7KB6MQ+woG8b8JrFKJRStYEuwI8i4izT9gP2d77JSqmIgryPBUhMTHSrnd6GiD7rdXoLZtcP1j1gdv1wfT5ITkxm7aMvcOOjvWhw6RRHqtbj9K/r6TjtE1QhXt0YiZ5twGuCLPAEmr3OHhUvBV4CugMDgCTgF6XUQFcnVErRu3dvlFKZm9mSBTjLJ2smzK4fimcuUU9idv1QdB8c2PwXBxq3pev8bygt6Wy9byB1juwhpNudbrbQvThyTzs2RxzQA28Ksk8Cf4nIzpwHROQlEZkpIptFZAFwD7AN+MDVCUUkM/uCYzNbkM0rn6xZMLt+KF65RI3A7Pqh8D7ISM9g3SvvU/PuTrQ6tY9z5aqwb8ZCOi6bRelyxf+Hq2OQpGMzfRYepdQtQDPg1YKUF5F0pdR84EOlVLCInMmrbLVq1dxkpXdSmJG4JRGz6wfrHjC7fiicD07tPcrxhwZwz95wALbfci9NF8+merD3+lHPNuAtPdlBaFN3CpM13dH3d/nzJGtWGDOi1yMSb8Hs+sG6B8yuHwrmAxFh04fT8LupHZ32hhPnV4bICZO4aetqynlxgAV920CxD7JKKT+gP7BCRM4XsI4v8AhwXETOuipbnOZ56oljwX7HVr58edq0acPXX39NWpqzcWTejUNvfqMGU1NTPWNQMcYs90BemF0/5O+DC2cusP7Ovtw1cihVE6+wu3E7Uv+KpO2oFz1kob7o2QY88qxMKfWw/Z+Ot+s9lVLngfMislEpVQ84DIwTkXE5qt8PVMb5gCeUUo8BfYAVwAmgOvCC/VqP5WebI/2ZWZg/fz61a9cmLi6O+fPn89ZbbxEXF8e4cTndbg6sebLmuwdyYnb94NoH/5u9lGovPsfdV86Q7ONL1LC3uOmr91ElaLU0PduAp15Izc/x+Rv7341AKNqj3VI471kPAi4Cy/I491GgGvAxWjBOAP4EeohIvvNpDxw4wA033JBfsRJD27ZtadSoEQDdunVj//79fPHFF6YNssnJybmStpsNs90DOTG7fnDug4SEZDY8+SrdFk3BVzI4XL0BAXN/okPnjgZZqR96tgGPPC4WEZXHFmo/Hm3/bHNSt4+IVBERp+vfiUiEiNwtItVFpLSIVBSRrgUJsEopRo4caboRxVlp37498fHxzJs3D6UUGzZsyHbc2WPXkJAQBg4cyKxZs2jatCmBgYHceeedHDx4kGvXrvHcc89RpUoVqlevzhtvvJHtcfSGDRtQSrFw4UIGDx5MUFAQFSpUYMCAAVy4cCHbtdPS0vjggw9o1qwZ/v7+1KxZkzfeeIOkpGxLV3PkyBHuu+8+ypQpQ9WqVXnllVcKvMiENU8W0tPTjTbBUMyuH3L74O+te9jZoiO9Fk7GVzLY1u8pQo7uolYJDLA2m42RI0fqNj7D1EMrRYTLly9nJvM2IydPnqRUqVKFTqC+adMmDh8+zIcffkhKSgqvvvoqDz30EA0aNKBRo0bMmTOHTZs2MX78eBo2bMjzzz+frf6rr75K165d+fnnnzl48CCjRo3i9OnTrF+/PrPMwIEDWbp0KSNGjOC2225j7969vPPOO0RHR7Nw4UJAW3v43nvvJTExkf/7v/+jWrVqfPvttyxatKhAOgIDAwuluyRi5ixUYOmHf3yQlp7B/PHfcffEt6iedIWLZSpy+b/T6PDkw/mcwXux2Wy8+uqrVKpUSZ9A6+7ced6yYU/anlceRad5QO25EA3fisD06dMFkH379klqaqpcvHhRJk+eLD4+PtKnTx9Zv369ALn84ah39OjRzH316tWToKAguXz5cua+L7/8UgB55plnstVv166dhIaGZn52XKd79+7Zys2ePVsAWbt2rYiIbNq0SQCZMWOG03J//fWXiIhMmTJFAAkPD88sk56eLi1atMhltzOsfLJWPlWz6xfRfHDwzBWZ1GOYpKNEQA41v0kSjx432jSP4GgD6JBPttiPLtab4OBgo03wKM2aNaN06dJUrlyZ559/nv79+/P9998X+jydOnWiYsWK2c4L0L1791zXO3HiRK76jz76aLbPjzzyCD4+PoSHa3PvVq5ciZ+fHw899BBpaWmZW7duWn6JTZs2ARAeHk6dOnW49dZbM8/l4+OT6/x5Yfb3sWC+eyAnZtcvImw55cParo/x4srJ+CAcG/YaDaO2EhBSx2jzPIKebcDUj4sLTQl4f/fLL79Qu3ZtypcvT7169YCijbANCgrK9tkxz8zZ/pzvUAGqV6+eq1xQUBCnTp0CICYmhpSUlDwfYzve3545cybXuZyd38LCIjfnLifyyuTNPDr5XR48tpW0Ur4kfzuVes8MNtq0EoPpg+yZM2do2rSp0WZ4jJYtW2aOLgaIj48nICAgM9DmzK+aczCSuzh37ly2zykpKVy6dIlatWoBUKVKFQICAti82XkK4po1tbyUwcHB7N69O9/z50Vqaqrpp/GY7R7IiVn1h207wZhv1jF5+URui9lPatnylA5bjO/ddxttmsfRsw2Y/nFx27ZtjTbBUBwDfxy92l27dmU7vmLFCl2uO2/evGyf58+fT0ZGBp06dQKgR48eJCUlceXKFTp06JBrcwTZTp06ceLECSIi/kklnJGRkev8eWENfLLuAbPpj0tMZfjUCEZNWMiChaO0ABtck9Lhf4AJAyzo2wZMHWQd2RfMPIXHkeIpODiYzp0788EHHzBz5kxWrlzJwIEDOXz4sC7X3b17N0899RSrVq1i0qRJDB8+nM6dO3PPPfcAEBoaymOPPcbDDz/M+++/z6pVq1izZg1Tp06lX79+HDhwAIBBgwbRoEEDHnzwQX744QdWrFhB3759iYuLK5R+MxMVFWW0CYZiJv1/7I/httErOLBoFb+tstEk7gzSujU7Jn0FJh1lbbPZdM3CY+rHxSLCli1buO2224w2xTCyNqzZs2czfPhwXn75ZQICAnj66acZM2YMQ4cOdft1v/zyS8LCwvjXv/5Feno6DzzwAF999VW2MrNnz2bSpEl8//33TJgwAX9/f0JCQujevXvmO1c/Pz/WrFnDiy++yPPPP0/ZsmV5/PHHue+++xg2bFi+dlhrF0OpErRyT1Ewg/7k1HQmLIriixV76H9oM19vnYZfehrccw9q4ULSnbxyMQs2m41u3bpx22236fJ9oKQEDOYpCkopERFiY2OdrvSxd+9emjdvboBlniU1NdWjI2w3bNhAly5dWLNmDV27dvXYdfNCb/3e0I7yugfMQknXv+fkZYZM3sLBIzF8tG0WTx38TTvwwgvw+edQunSJ90F+OPQrpRARt0ZaUz8uhtzvIM2Gs5G/ZsLs+sG6B0qq/owMYdKve7nz3ZWk/72T39e+pwVYf3+YMgW+/hrsPzBLqg8Kip76Tf24GKB27dpGm2AoZp8nanb9YN0DJVH/idhrDJsawe+7z/Dy3uWM3bkQ37RUaNQI5s2Ddu2ylS+JPigMeuo3fZDNOWXFbHj6dUFoaGixWi+4ONliFGa/B0qSfhFhXng0b8zcRtC5U6zdOoWbz+zVDj73HHzyCTiZe16SfFAU9NRv6sfFSileeOEFU48uLom5ZAuD2fWDtvCHmSkp+i9eTWbw//3BkP9uoe/fq/nfr6O0AFujBqxYAZMnOw2wUHJ8UBRsNhsvvPCCboMgTT/wKT4+nvLly+c67g0DVtxBenq6KUZX5oXe+r2hHeV1D5iFkqB/XdQZhk+LIP30GSZv/Y57T+zQDjzyCPz3v1Clisv6JcEH14NDvzXwSQe2b9+e5zEz/ABJSEgw2gRD0VO/t7QfV/eAGfBm/QnJabw5cxt9P17PzZEb2bHibS3AVqoEP/4Ic+fmG2DBu33gDvTUb/p3snmt+FO6dGkSExMpU6aMhy3yLGafJ6qn/sTERPz9/XU7v7sw+6pX3qp/x5ELDP02nHPHzjDtz5n868jv2oF774Xvv4dCDObxVh+4Cz31mz7IhoSEON1frVo1Tp06Ra1atQgMDCyxwcgbgoCeuFu/iJCWlkZ8fDyxsbFekaggr3vALHib/rT0DD5btocPFkdxx8koVmydSvX4WAgMhI8/huefh0J+X3mbD9yNnvpNH2T37t3r9IuwQoUKAJw+fZrU1FRPm+UxkpKSTL1Avh76fX19CQgIoG7dul7h27zuAbPgTfoPn4tn6OQt7Np/mv/smMPw/au1Ax07wsyZ0KRJkc7rTT7QAz31mz7IuvoFU6FChcxgW1KJjo429a9Ys+sHqxfjDfpFhOkbDvP2j9tpdvoAERHf0uDSKfD1hbFjYeRI7d9FxBt8oCdWT1YnlFK0atWKBx980LTTeOLj4402wVDMrh8sHxR3/TFXEnnhu62s3XGct6KWMGLXEkplpEOLFjBrFrRvf93XKO4+0BObzcaiRYt0SxRh6iArImzYsIHQ0FCjTTEMvfLFegtm1w+WD4qz/qXbT/DS9//jhpNH2LhlMq1jj2jvW19/HSZMADe9jijOPtAbm81GaGgooaGh3pkgQClVGxgBdADaAIFAfRGJLkDdvIxrJyKRWcr52K/xHFAD2A+ME5GFrs7tap6sWbD0m1s/WD4ojvrjE1MZ8eN2Zm88xLB9q3k/ci7+aSlQty7MmAFu7hgURx94Em+fJ9sIeBS4BGwuQv0fgE45tgM5yrwP2ICvgZ5ABDBfKdUrv5Nb88Ms/WbH7D4obvq37I/htjErWP/r/1i2biIfbZulBdjBgyEqyu0BFoqfDzyNt8+T3SQi1QGUUkOAboWsf0pEIvI6qJSqBrwJTBSRT+y71yulGgETgRWuTl4uj2XGzIKl39z6wfJBcdGfkqblfP182W76H/6dz7bPpFxyAlStqmXN6dtXt2sXFx8YhZ76dQ+yIpKh8yW6A37A7Bz7ZwPfK6Xqi8jRvCrXrFlTT9uKPZZ+c+sHywfFQf/ek5cZMjmc0weOMXvr9/Q+/qd2oE8fLcBWq6br9YuDD4xET/3esKzicKVUslIqQSn1m1LqzhzHbwSSgUM59u+2/23h6uQHDuR88mwuLP3m1g+WD6z3e2oAACAASURBVIzUn5EhfL1yH3eOXUmtP9axbflILcCWLw/Tp8Mvv+geYMFqA3rqL+5BdjbwPNAVeBaoAvymlArNUqYycFlyj+C6mOW4U5RSPPnkkyilMjezTeVp2LCh0SYYitn1g+UDo/SfvHCN3h/9xn9++J1PN33LvA2fUiXhivbONSpKewfroZXmzNYGbDZbtu99RxzQA5dBVilVSin1lFJqilJquf3vYKWUR9K2iMgTIjJXRDaLyGzgDuA0MD6rmYCzUcj5ekxECAsLQ0QyN7MFWTMP3QdLP1g+8LR+EWHelmhuHb2CtPUb2LpiNIMObQB/f/jsM1i3DurV86hNZmsDNpst2/e+Iw7oQZ7vZJVS9YBVQBPgJHAOaAUMAUYopXqIyDFdrMoDEYlXSi0Hnsmy+yIQpBxzcv4hKMvxPLl8+bKbrfQuLP3m1g+WDzyp/+LVZF6f8SfLthzincgFvLR3BT4i0K6dtrDEjTd6zJasWG1AP/2uBj59DVQA7hCRLY6dSqnbgfnAJKC3bpblTc6e627AH2hI9veyjnexe1yd7KabbnKrcd6Gpd/c+sHygaf0/7brDMOnRlDl8F42b5lM80snkFKlYNQoGDMG/Pw8YoczrDagn35Xj4vvBt7OGmABROQPYJT9uEdRSlUA7gO2Ztm9EkgBBuQoPhDY5WpkMVjzwyz95tYPlg/01p+Yksa/Z2/jwYlreXzTHDauHEvzSyegcWPUH3/AuHGGBliw2oBR82SvAjF5HIsBCpztWin1sP2fjp8LPZVS54HzIrLR/mj6MNoqTePsdd4EmgLr0d7D1kObD1uDLAFVRGKUUp8Dbyul4oEdwL/QfgT0yc+2SpUqFVRGicTSb279YPlAT/1/Hb3IkMlbSN9/gNVbJnPL+YPagRdegA8/hLJldbt2YbDagH76XQXZ2cAw4Fcnx54DZhbiOvNzfP7G/ncjEIr2CLgU2XvW+4F+9q0iEAf8ATwjIv/Lcb7RaD8KXuGfZRUfFZGl+RlWpUqVQsgoeVj6za0fLB/ooT8tPYPPl+/hP4t2MmjvWj746ycCU5OhVi0toXq3wq7Joy9WG9BPv6sgewh4RCkVBSxEG/hUHXgYKA/8qpR62lFYRL7P60T5rQVpX8dY5di3FMg3SNrLpqONOB6fX9msKKUIDQ2lc+fOphtV7ODw4cPUqVPHaDMMw+z6wfKBu/UfORfP0G/DObHzAPPCp3Lv6Z3agQEDYNIkCApyfQIDMHMbsNlsbNy4kQ0bNuhy/jwTBCilCrNSk4iIR6b1uAvHYOTTp0+berUTS7+59YPlA3fpFxF+2HCYt3/aQY99m/nizx+olHwVKleGyZPhkUfcYK0+WG1A069HggBXPdn67rxQccVqXJZ+M+sHywfu0O/I+bp16wG+3jqdh4/Zl1vv1QumTYPgYDdYqh9WG9BPf55B1tNzYI3i6tWrRptgKJZ+c+sHywfXq3/Z9pO89P1W2u77H39GTKV6wiVtQNPnn8OQIR5btel6sNqAfvp1zydbXLHyyWpY+s2tHywfFFV/fGIqI3/awcK1uxm/4yeGHFinHbjjDi3na4MGbrZUP6w2YFA+WaVUd6XUL0qpPUqpIzk3dxpiFNb8MEu/2TG7D4qiP/zAeW5/51f2z/+VLctHMeTAOsTPDz76CDZs8KoAC1YbMGSerD3h+VJgLdAMbdGHMsDtwDGKloC92GENXbf0mx2z+6Aw+lPS0vngl118veRvRkYu4NU9yyklGdCmDWrWLGjVSkdL9cNqA/rpd9WTfQf4P6CX/fMYEQlFSy1XCufzZ70KpRS9e/c27fQdwNSPiMDSD5YPCqp/36kr3P3eapbMXsO6Fe/yxu6l+Cjg7bdh61avDbBg7jZgs9no3bu3IVl4mqH1ZDPQ1gr2BRCRA4ANLQh7NSLC9OnTTR1ko6OjjTbBUMyuHywf5Kc/I0P4ZtU+7nhnBa3WLeb3FWNofekYNGyI2rwZ/vMfLYOOF2PmNmCz2Zg+fbrns/CgBdc0ERH7Eoh1AcdKS6fRFuT3epo3b260CYZi6Te3frB84Er/qYsJDJsSzvbIo3yz9XsejQ7XDgwcCN98oyVXLwFYbUA//a56svuBEPu/twGvKqWClVJVgTeAaN2s8iBm/gUHln6z6wfLB3npnx8eza2jlnM84m9+XzVWC7Bly2ojh2fNKjEBFqw2oKd+Vz3ZHwFHeB+LNgDqpP1zOvC4blZ5kMTERKNNMBRLv7n1g+WDnPovXUvh9Rl/siDiGJ1i9rNw8xeUT4jT3rnOnw9NmxpkqX5YbUA//QWeJ6uUqg30QBthvFZEXOZpLe5Y82Q1LP3m1g+WD7LqX7/rLMOmhnP6UiK3xR1j2arxlE5KgPvug59/LlG916xYbcCgebJZEZGTIjJNRL7y9gDrwBpdbM2PM7t+sHywfft2ElPSGDF7O70/+o3TlxJ5sFISKzZ/qgXYAQNgyZISG2DB3G1A79HFrh4Xl3hEhD179tCiRQujTTGMatWqGW2CoZhdP1g+uJhRjjvfXcn+03H4llKMeaAZr7/7BOriBejRA6ZPh1Jelf+k0Ji5Ddhstsw4oEegdbUYhWPqjjMEuIKWIP1jEVntdss8hJ+fn9EmGIql39z6wbw+SM/I4LNle/nPooOkZQhNgiswbdhttPvlB9i9W1u1af58KF3aaFN1x6xtwIGe+l09Ln4fOAGcB34APgRm2D+fBGYBVdHyyt6vm4U6c/LkyfwLlWAs/ebWD+b0wdGYq3SfsJZxC/4mLUMYdm8TNo/rQTv/JHC8PvrqKyhXzlA7PYUZ20BW9NTv6nFxEnAU6CkiSY6dSqlAtNWezgPtgeXAKGCZblbqSMuWLY02wVAs/ebWD+bygYgwc9MRRv64natJaQQHBTLx0eY8eHszrcCYMRAfD336aIOdTIKZ2oAz9NTvqic7DPg8a4AFEJFE4HNgmIhkANOA1rpZqDMHDhww2gRDsfSbWz+Yxwfn45Lo/8UmXvxuK1eT0njwlrpETOhFDXVRK3DlijaCGOCTT4wz1ADM0gbyQk/9rnqy1YC8Xkb4AY4VlWOB4p8wMQ/S09ONNsFQLP3m1g/m8MHyHSd58butxMYnU7FMaT578mYe6VQPpdQ/+ufNg6QkCA2FRo0MtdfTmKENuEJP/a56stsAm1IqOOtOpVRNtMUpttl31UNbZtHrsKbwQCsvXtTcHZhdP5RsH8QnpvLid1vp/8UmYuOT6dyiOhETevHobSGZI0kz9c+Yof0dPNgYYw2kJLeB/DAyQcArQG3gqFJqvVJqrlJqPXAEqAm8bC/XCPhJF+t0RkT45ZdfTB1kIyMjjTbBUMyuH0quDyIOajlfZ2w8jH9pHyY+3p6wf99N7Spls5WLjIyEQ4fgjz+0ZRMfesggi42jpLaBgmCz2fjll190SxCQZ5AVkR1oAfQztGQBrex/PwUai0ikvdy7IjI2r/MopWorpSYppcKVUglKKVFKheRnmFKqg1JqilJqn73ecaXUj0qp+k7KRtvPm3Prm991goOD8ytSorH0m1s/lDwfpKSlM27B33Qfv5ajMVdpXTeIze/14IUezfDxyd1bCQ4O1tYiBi3AmmREcVZKWhsoLHrqd7kYhYhcQBs5fD00Ah4FtqMleu9WwHr90XLXfgXsBmqhpdfbppRqKyIncpRfhZaCLyv7i2izhYWFF7Lv1BWGfruFyOhLKAWv39+C0Q+2ws83n8UkFi/W/g4cqL+RFqaiwMsqXgebRKS6iPQC5hei3ocicruIfCMiG0XkJ7S1k4OAoU7Kx4pIRI7tUn4XOXPmTCFMKnlY+s2tH0qGDzIyhP+u3s+d764kMvoS9W4oy8pRXXnv0bb5BtjY3bth504IDIS77vKQxcWLktAGrgc99WfrySqlZhairojIoAIUyii0VVq98072HbPntq1VlHM6o23btu46lVdi6Te3fvB+H5y6mMDwqRGs330WgCfuasDEATdRIbBgKzW1vXBB+8ddd3l98vWi4u1t4HrRU3/OnuxdwJ2F2DyKUqo52tSivU4OP2B/d5uslIooyPtYgKioKLfa6G1Y+s2tH7zbBwsjjnHrqOWs332WKuX9+fHlO/lmyK0FDrAA8Y5HxV276mRl8ceb24A70FN/tiArIiEiUr+AWwPdrHKCUsoXmIy20tR3OQ4vBV4CugMD0Far+kUp5fIFi1KKESNGoJTK3Mw20rhUCV/4PD/Mrh+80weXrqXwzH//YPA3f3A5IZXubWqydUIveneoU7gTiRC0zT4b8d573W+ol+CNbeB6sNls2b73HXFADwqcT9YtF1NqCDAVqC8i0YWsOxl4Brgvv4QESqlSQARQQ0Sc3nWOfLKxsbHccMMNhTGlRGHpN7d+8D4fbNh9lmFTIzh1MYEyfqX4YMBNPBXasGhfkvv3Q7NmULUqnD0LPp4YplL88LY24G4c+g3NJ2skSqkPgGeBpwuS8UdE0tEGWdXOuZhGTnbt2uUeI70US7+59YP3+CAxJY2RP27ngQ9/49TFBG5uWIUtE3rxdJdGRe+FrFmj/e3a1bQBFrynDeiFnvqLfT5ZpdRoYCTwsojMKkxV+1+XXfXatWsX1bQSgaXf3PrBO3zwd/RFhnwbzr5TV/AtpXi7bytev78FvqWuMzBmDbImxhvagJ7oqb9Y/3RTSr0MjAdGi8ikQtTzBR4BjovIWVdlU1JSrs9IL8fSX3T969evp02bNiilGDdunBut8izFuQ2kZ2Tw6dLddHlvNftOXaFxcAXWvdONf/dpef0BNi0N1q/X/m3yIFuc24An0FO/R4KsUuphpdTDwE32XT3t+zrbj9dTSqUppd7NUqc/8AWwEvhNKXVrlq1FlnKPKaXmKKWeVEp1sddbb7/WiPxsi4mJcZ9QL8TSX3T9Xbp0YenSpQC0aNEin9JFIyEhgcDAQDZv3qzL+aH4toGjMVfpMWEdtvl/k5qewbB7m/D7uB60b1Al/8oF4X//g/h4EurUgbp13XNOL6W4tgFPoat+EcncgE5ZP7trQ3tk62zbYD8eYv9sy1Lnh/zq2cvdCvwGnANSgSvAWqB7fjaJiMTFxYmZsfQXXf+bb74p9evXF0CaN28ubdq0kffee8+N1oksWbJEgoKCJC0tza3nzUpxawMZGRkyY8MhqTF0rpR74kdp/PIiWbvztPsv9N57IiDJQ4e6/9xeRnFrA57God8eF9wa/3L2ZH9XSp1RSn2rlOqhlCr4ZDPXgVzlsYXaj0fbP9uy1BmcXz17uQgRuVu0VaVKi0hFEekqIqvys8vKwgPbt2832gRDuR79H3/8MU8++SRBQUHs2bOHyMhI3n333fwrFoJly5bRo0cPXadYFKc2cD4uice+3MwL9pyv/ew5X+9ppcPasvb3sQdM3ouF4tUGPI2ns/DUAsahpa9bDMTas+/0V0pV0MUCAxERJk6caOogGxgYaLQJhnK9+nfu3Em7du2YOXMm3bt3p1q1anTv3p0LjlWE7ERGRtKxY0c+++yzXOeIiYmhevXqTJqUe9jBihUruP/++zM/p6en8+mnn9K+fXtat25NvXr16NOnD6dOnSqyhuLSBlb8dZKOo1awfMdJKpYpzbRhnZjxwu1ULqfDKkzx8RARAT4+JHXq5P7zexnFpQ0Ygc1mY+LEibpl4XH1OLU88BgwB7iMtsDDKmAYUNPdXWpPb9gfF589e9b58wOTYOm/Pv0NGjSQ+vXry9q1a0VE5NChQwLIxIkTM8vs2bNHKlasKLNnz87zPMuWLRMfHx/ZvXt35r4dO3ZIqVKl5MKFCyKiPUbt16+f9OzZM3PfpUuXpEaNGvLEE08UWYPRbSA+MUVe/C5Cyj3xo5R74kfp9Z+1cvz8VX0vunSpCIjceqvh+osDZveBQz8eeFycNfjGi8jPItIfqAr0A46iZcI5oZTaqpQaqU/o9xx79zpbodE8WPqLrv/q1ascPXqUfv36cc899wDQsGFD/Pz8iIuLyyz3zjvv0Lx5cwYMGJC5LyIighMn/kkkdeutt5KRkcE2x+pDaI+KO3XqROXKlQGYNGkSv/76KzNnzszcV7ZsWXx8fEhNTS2yDiPbwNaD57l9zK/8sOEwfr4+/OexdiwdcTd1biibf+Xr4ddftb/33mv6ewCs7wE99RdodLGIpIrIryIyTERqAXcAG4AndbPMQ4SEhBhtgqFY+kOKXDcqKgoR4aEsSb5jY2NJSUmhrv0937Vr1wgLC8sMwgBxcXF06tSJNY45msDZs9pMM58sCyIsW7Ys26Pir776il69enHDDTdw5coVVqxYQdeuXUlPT2f06NFF1mFEG0hNy+D9BX/TbfxajsRcpWWdSmwe14OXejZ3mvPVrSQlwc8/a//u18/09wBY3wN66i/SYhQiEg6EU4ApMsWd+Ph4o00wFEt/0fX//fff+Pr6ctNNN2Xuc6wc07JlSwDCw8NJTU3NluVj586dALRu3Tpz3549ewAyg3NMTAzbtm3j+++/zzzv4cOHSUpKok2bNvj7+9O4cWMGDBhA//79qVCh6EMmPN0GDpyJY+jkLew4ehGl4NX7mjPmwdb4l/bQ+rmLF8OlS9CuHbRrR7zJF8cH63tAT/3FfsUnvck5QMVsWPqLrn/nzp00adIE/yzp0Xbu3IlSKjOA7tu3D4B69eplK+Pj48ONN96YuW/jxo34+vpmBuMVK1ZQp06dzDInT54E4Msvv8zWc3YHnmoDIsLUdQcZM+cvElPSqXtDWb59thN3NKvmketn8p09v8iQIYB1D4DlAz31F+sVn/TGmsJDtl6YGbke/Tt37qRVq1bZ9kVFRRESEkL58uUBKFOmDKCNCnawfPly/Pz8Mkd0JiQkMG/evGw90pyPiqtV0wLRuXPnsl1v9erV/Op4v1hEPNEGzlxK4MFPNvDGzG0kpqTz+B312TK+p+cD7NGjsHYtBATA448D1j0A5vaBp6fwmAoRISwszNRB1szz46Do+kWEqKgop0G2TZs2mZ979+5NhQoVWGzPWbpo0SKOHz9OhQoV2LJlC2lpabz44otUqVKFL774AoDU1FTWrFmTLci2a9eOVq1aMW3aNK5evQrAwYMHeeWVV6hfv36RNDjQuw0s/t9xOo5awdqoMwSV9WP2S3fw7bOdqFjGT9frOuXLL7W/Dz0ElSoB1j0A5vaBzWYjLCxMtyk8pn9cXK5cOaNNMBRLf9H0R0dHExcXly3Iigi7d+/m8OHDPPXUU0yfPp0bbriB5cuXM3z4cBYsWECrVq1Yt24dO3bsYMCAAfj5+dGlSxfCw8OpZP/S37RpE2lpaYSGhmaeWynF4sWLee6552jYsCGNGjWicuXK/PTTTzRr1swQH+THlYQU3pq1nZ//OArAva2D+WbIrdSoZNCczC1b4KuvtGw7r72Wudvs9wBYPtBTf575ZO3rCE8TkdNOjgUDQ0XEa1dFd+STPX36NDVr1jTaHMOw9Bc//a+99hpHjhxhyZIlHrmeHj74fV8Mz00J53jsNQL9SjG+fzuG3tNYt0dy+ZKYqA102r8fRoyAiRMzDxXHNuBpzO4Dh35P55MdC+SV/6em/bjXc+DAAaNNMBRLf/HTv3z58myPivXGnT5ITk1nzJy/6PXBWo7HXqN9/cpsHteDZ7s2MS7AJifDwIFagG3eHHK8HiqObcDTmN0Heup39bjY1R0RBCS72RZDaNiwodEmGIqlv/jp9/QXnrt8sPvEZYZM3sKuE5fxUYp/97mREX1aUtrXwKEfcXHa+9e1a6FiRZg9Wxv0lIXi2AY8jdl9oKf+bEFWKRUK3J1l13NKqZw/qQOB+4DdulnlIZRStG3blj59+ph28NOFCxeoU6eO0WYYhtn1w/X7ICND+L9V+7DN/5uUtAwaVCvHlOc60bFxVTdaWQQ2bYJBgyA6GqpVg9WrIcugNAdWGzC3D2w2G0uWLCEyMlKX8+fsyXYGxtj/LcBTTuqkAHuAl3WxyIOICBs2bMg2wMRsXL582WgTDMXs+uH6fHDuciLPTgnnt13ailWDQxvywePtKRfglgReRePiRRgzBiZPBhFo3x7mzoVGjZwWt9qAuX1gs9kIDQ0lNDRUl1cargY+ZQC3isj/3H7VYoBj4FN8fHzmnEYzYuk3t34oug/W7jzNs1MiOB+XRJXy/vzfMx25r31ewzg8QGoqTJ0KY8dCbCz4+sLIkfDOO+CX93Qhqw1YPnDo9+jAJxHxKakBNitmnh8Gln6z64fC+yA1LYMxc/6i3ycbOB+XxF3NqxM+vqdxAVYEliyBli3hhRe0ABsaCpGR8P77LgMsWG0ALB/oqT/febJKqRpAXSAg5zER2aSHUZ7EMTfRrFj6za0fCueDC/HJPPH1ZjbvjaGUj2L0g614/f4WlPIxaHDTn3/Cm29q719BeyT80UfQty8U8NGf1QYsH+ipP88gq5SqBcwG7nJ2GO2drYdW9NaPKlWqGG2CoVj6za0fCu6DI+fi6fvxeo7GXKV6xQBmvXQnnZoYNLgpOhpGj4afftI+V6miPSZ+7rl8e645sdqA5QM99bv6+flfoCXwb6An2qhjx9aF7KOQvZbDhw8bbYKhWPrNrR8K5oPjsde4f+I6jsZcpW1IEBvf62FMgL18Gf79b2jaVAuw/v7a4hKHDsFLLxU6wILVBsDygZ76XT0uvhN4WURm6XZ1g1FKERwczLPPPmvaKTxNmjQx2gRDMbt+yN8H5y4ncv/EdZy4kEDHxjew+K0unh89nJKijRZ+7z1t9DDAgAEwYQJkyXBUFKw2YG4f2Gw2pkyZwpkzZ3Q5v6uebCIQo8tViwlWggBtOTEzY3b94NoH6RkZDJm8haMxV7mpQWUWvhHq2QArAgsXwo03wiuvaAG2c2ftXezs2dcdYMFqA2BuH+idIMBVkJ0KPHG9F1BK1VZKTVJKhSulEpRSopQKKWDdAKXUx0qpM0qpRPs5cr0jVkr5KKXeVkpFK6WSlFJ/K6UKlHTTkdHErFj6za0fXPtg+vrDbNhzjhvK+zPn1c6ezZwTEQF33AEPP6w9Dm7WDMLCYP166NDBbZex2oDlAz31u5on+ywwEogGVgAXc5YRke/zvYC2itRcYDvaQKluQH0RiS5A3R/RVpd6CzgCvID2friTiERmKTcBeBMYbb9Of2AocL+IrMjj3NY8Waz5cWbXD3n74PK1FNq8tZSLV5OZ+eId9LulrmcMOnwY3n4b5s/XPletCuPGaUnWfd2fOMxqA5YP9Jwni4g43YCMfLb0vOrmOI9Pln8PQRuVHFKAem3sZZ/Kss8X2A+EZdlXDW0d5fdy1F8H7HRxfhERWb9+vZgZS/96p/vPnz8v3333nWRkZHjWIB1Yu3at/Pnnn3kez8sHny/bLeWe+FG6jV/tGT9cuCDy2msipUuLgEhAgMjo0SJXruh6WbPfAyKWDxz67XEh37hWmM3V4+L6+WwNChjEMwpSzgm9gVS0XrDjXGnAHKC7Usrfvrs74Ic23Sgrs4FWSimXGa2toeuW/qzExsYyatQomjZtSmRkpG7vaTxJXFwcffv25YEHHmDbtm25jjtrA+kZGUxbdxCA1+5roW8GneRk+PRTaNgQPv8c0tK0NYcPHoTx46FCBf2ujXUPgOUDPfXn+exFRI7pdtWCcSNwVEQScuzfjRZUG9n/fSNaT/aQk3IALYCjeV3EzI9IwNKfVf+jjz7KsmXLaNGiBU8++SQVK1Zk3DivTZmcjcGDB7Njxw7uvvtugoODWblyJfXra78/nbWBDbvPcSz2GvWrlaNba53yjIrAvHnao+Gj9lv0nnvgk0+gbVt9rukEs98DYPlAT/35LtOilGqtlHpRKTXWvvoTSqlGSim9/1cqA5ec7L+Y5bjj72XJ3eXIWS4XSimeeuoplFKZm9lGGkdHRxttgqHkpb8k9GDzImev1JkPlu84CcAjt9bDx0eHXuzvv0OnTtC/vxZgb7wRVqyANWs8GmDBugfAfD6w2WzZvvcdcUAX8nqODPgD84F07O9ggfb2Y4uAiYV9Nk3h3smuAcKd7L/Xfo477Z+nAmeclGtsL/dEHucXEZGzZ88W5RF+icHSn11/bGysjBo1SipXriwvvviipKenG2SZ+1i0aJHUqlVLHnjgAdm2bVuu4zl9kJGRIc1e+UXKPfGj/HnovHuN2b9fpF8/7Z0riNSoITJ1qkhqqnuvUwjMfg+IWD5w6MfD72QnAF3RpvFUJ3sS91/R3oXqyUWc90KDshx3/A1SuX+G5CznFLP9gsuJpT862+cqVaowYcIEDhw4QLt27fR9F+khKlasyOLFiwkLC+Omm27KdTynD/acvMLJiwlUrxhA+/puelcVGwsvv6z1WH/5BcqU0ZZBPHhQt1HDBcXs9wBYPtBTv6uW/RgwRkR+UkrlXKP4KBCim1Uau4F+Sqkykv29bAu0nLaHspTzBxqS/b1sC/vfPa4ukpiY6B5rvRRLv3P9VapU4emnn/awNfpw992uV0DN6YOtB88D0LlF9et/VJyUBF99pa3MFBenLdr/zDPalJyaOr3rLSRmvwfA8oGe+l31ZKsAe13U88/jmLsIA0oDjzh2KKV8gX8Bq0Uk2b57JVrQHZCj/kBgl4jkOegJcPrL3kxY+s2tH3L7YNuRCwB0aHhD0U+akaGtLdy0qba2cFwcdO+upZ+bNq3YBFiw2gBYPtBTv6sgexTolMexW9DmqxYIpdTDSqmHAYeSnvZ9ne3H6yml0pRS7zrqiLbYxFzgC6XUEKXUPWjTd+oDY7OUiwE+B95WSr2ulApVSv0XLYHBqPxss/IoWvrNTk4fbLcH2ZsaFPFR8caN0LGjtrbw8ePQqhWsWgUrV0Lr1tdrrtux2oDlA6Pyyc4ERimlotEGOgGIUqoL8BpgK8R15uf4/I3970YgFO19bylyB/2n0N4NjwcqAX8DPURkF4JpXgAAIABJREFUR45yo4GrwCtADbQfAI+KyNL8DKtWrVqBRZRELP3m1g/ZfRCfmMreU1fwLaVoXTfIRS0n7N+v9VqXLNE+Bwdr81wHDYJSxTcrptUGLB/oqd9VkP0IbdWlWcA0+77f0ZK3zxGRSQW9iOSzTJVoSyzmKiMiicDr9s1V/XS0QDy+oDaBNpWhYcOGDBw40HRTdxz4FSE1WEnC7Pohuw+ijl9CBFrWqUSAXwEDY0yMlh3n228hPR3KltWC7euva/8u5lhtwNw+sNlszJ49W7d0d3k+LhaRdBHpD3QGPkULtF8Bd4tIzvefXomIMG3aNNMGWICTJ08abYKhmF0/ZPfBobPxADSrWTH/iklJMHEiNGoE33yjTcp59lltMf933vGKAAtWGwBz+8BmszFt2jTd5sY77ckqpfyACGCkiKwGNuty9WJAy5YtjTbBUCz95tYP2X3gCLINa+Sz1syvv2pJ0h2//nv1go8+0qboeBlWG7B8oKd+pz1ZEUlBG2CUptuViwkHDhww2gRDsfSbWz9k98Ghs3EANMoryB4/Dg8+qAXVw4e1oLpmDSxf7pUBFqw2AJYP9NTvanTxGrS0dCWa9PR0o00wFEu/ufVDdh8cPqf1ZBvVyLEov4g29aZ5c20xiXLltDWG//oLunb1pLlux2oDlg/01O9q4NMkYLZ9bupi4AzaMoWZiMgR3SzzEK1atTLaBEOx9JtbP/zjg4wMyQyyDatn6ckmJsLQofDjj9rnRx7RsuXUquVpU3XBagOWD/TU76onuxGogzaydyNwADiYY/N6IiMj8y9UgrH0m1s//OODkxcTSE7NoHrFAMoHltYOpqTAQw9pAbZsWZg1S8ucU0ICLFhtACwf6KnfVU/2Kd2uWkxQStGkSRMee+wx044wDg4ONtoEQzG7fvjHB8djrwIQUq2cdiA9XVtQ4tdf4YYb4LfftIUlShhWGzC3D2w2Gz///LNu72XzGl1cGrgM7MxvWUJvRkTYv38/TZs2NdoUCwvDOX1RW7+1VlAZbce338KCBVCxIqxeXSIDrIWFzWbjscceo2nTprokBMlrdHEqMA/9kwAYzpkzZ4w2wVAs/ebWD//44PQlLQ9HcFAgXLmiZckBbcBTu3ZGmac7VhuwfKCnflfvZI8AJX6trbYeThBd3LD0m1s//OODM5e0nmzNymXgww+19HS33669ky3BWG3A8oGe+l0F2Y+A0UqpqrpdvRgQFRVltAmGYuk3t374xweOnmyD9Hht9DDAp59q6elKMFYbsHygp35XA5/uRkuaflQpFUHuKTwiIoN0s8xDlCrGC5d7Aku/ufXDPz44be/JNt+2QVsy8f77tWw6JRyrDVg+0FO/q57sHUAqcB4tIfodwJ05Nq9GKUXv3r1NO7IYoEmTJkabYChm1w//+OCMvSdbI2KjdqBvX6NM8ihWGzC3D2w2G71799Zl0BO4ThBQP5+tgS4WeRARYcGCBaYOsrt27TLaBEMxu37QfJCRIZy5nIhfeiplfrcH2e7djTXMQ1htwNw+sNlsLFiwQLcEAa56sqagdu3aRptgKJZ+c+sHzQex8UmkpQvd4o6grl3TpuuYxDdWG7B8oKf+PN/JKqXq5ldZRI671xzPk5KSYrQJhmLpN7d+0Hxw9nISAPfH2AeA9OxpoEWexWoDlg/01O+qJxsNHM1n83piYmKMNsFQLP3m1g+aD2LjtSB7x7G/tJ0mCrJWG7B8oKd+lddzaKXUYHIkBACqAPcBDYD3ReR73SzTGaWUiAjx8fGUL59P7swSjKXf3PpB88HKXRcY88lS9i96Wcuwc+EC+PkZbZpHsNqA5QOHfqUUIuLWEVCuBj79ICIzcmyficg9wO9ogdbr2b59u9EmGIql39z6QfPBhfhkWl86pu3o2NE0ARasNgCWD/TUX9SBT7OBp91piBEopRgxYoSpRxcHBgYabYKhmF0/aD64cDWZRnFntR3NmhlrkIex2oC5fWCz2RgxYoTnp/DkQzUgoKCFlVJ1lFILlFJXlFJxSqlFBRlYpZSyKaUkjy0pR9noPMrlOdlPRAgLCzN1kA0JCTHaBEMxu37QfBAbl0zjOPv6rSabM2m1AXP7wGazERYWptsUHleji+9ystsPaAm8DWwuyAWUUmWA34BkYBDae97xwHqlVGsRueai+jRgZY59Ze37wpyUXwXYcuzb78q+vXv3Ur16dVdFSjSWfnPrB80HF676/hNkTZaVymoDlg/01O9qWcUN5B745OhPbwSGF/AaQ9He3zYVkUMASqmdaEnfnwM+y6uiiJwETmYzQKkn7HbPcFIlVkQiCmgXYO5fcGDpN7t+0HxwIfwwjUwaZK02YPlAT/2ugmwXJ/uSgGMicrYQ1+gNRDgCLICIHFVK/QH0wUWQzYNBwDm0Xut1Ex8f747TeC2WfnPrB80HiRcuEZx4mQw/f3zq1DHaJI9itQHLB3rqzzPIishGN13jRmCJk/27gUcKcyKlVG204P+FiKQ5KfKAUioBKAX8xf+3d+bxUVXn/38/SUjYIWwCsoQlbBUaRRDqQkTFBUQQUBABbbV+pVW/Ym2h+i2DVcHlR20VEKHFBVCEqhREcANEhVotkU02ISKIAdn3JXl+f5w7YTKZJJPM3EySez6v133NzLlneT7PPWeee855zjkwQVXfKSzPffv2FUeECgfL39v8weig1vdm2Xt2q1bEeWyzeFsHrA7c5F+g45OIdBORWwq4N0hEwj2eow5wIET4fiA5zDz8GIaROdRQ8QLgPuBaYCim1/22iNxeWIadO3cupggVC5a/t/kDXHTRRdT70WzeFucxz2KwdQCsDtzkX5h38XhMLzQU2jv3w0Uot62S+EsPB1ar6pp8Bajep6qvquoKVZ0HXAV8WZic/lN4RCT38pqnsV0f523+ACtWfUnLg7sAiG/vPSNr64D3dODz+fL878fkFB7g50BBTkRfAJ3CLOMApjcbjGRC93BDQkS6Au0I3YvNB1XNBuYCTUSkUQFxeOaZZ1DV3MtrRrZ69eqxFiGm8Dp/gLNxlUk95M3lO2DrAHhPBz6fL8//vt8OuIHCHJ8qU7ARjscspQkH6wndI+4AbAgzDzAOT2eB2cVI4381KVB7jRs3LkZ2FQ+Wv7f5A1SpWZcUj3oWg60DYHXgJv/CerLfYDyDQ6EvRaw/DcC/gG4ikrsNo4ikAJcSeq1rPohIIjAYWKSqe8NMk4BxrNpRmDf05s2bw8muwsLy9zZ/gG+2bD+325MHe7K2DlgduMm/MCP7InC3iDwjIm1EpKqIpIrIM8CvgMlhljENc6LPfBG5SUT6YryNvwem+iOJSHMROSsifwqRRx/MkHPIoWIRGSIib4jIcBG5UkQGA0uBzsAfChOuVatWYdKomLD8vc0f4DwVqmWf4nC1WlAn1MxOxYatA1YHbvIv7ICAaZg1rA9ierVHgI3O77+o6kvhFODs6NQT2Ay8BszCHJPXU1WPBkQVzDB0KJlGYLyRFxZQzHbMVo/PAO9jjPcp4DpVfaMw+azruuXvdZzcYpaw76vnzSFDWwesDtzkX9icLKr6OxGZAlyNOebuJ+BDVd1WnEKcw90HFBEnkwI8jlX1piLSrsIY8mLj4MGDJUlWYWD5e5s/gO4ym6odq9sgxpLEBrYOWB24yb9QIwugqt8C37omQQwhItSoUYNRo0Z5zqvYD7s+ztv8Ac47a/Z1OVnfm3vX2jrgbR34fD4mTpzo2q5PJT2Fp0LAnsLjvfVxwfA6f4Cz35ndnk6f1zDGksQGtg54Wwdun8LjaSMLULt27ViLEFNY/t7mD1D7kFmuntMw5HLyCg9bB6wO3OTveSNbt27dWIsQU1j+3uYPUOvIfgDk/PNjLElsYOuA1YGb/D1vZL/9tkJON4cNy9/b/AFqHjBLzxOaeNPI2jpgdeAmf88b2TYeXHwfCMvf2/wB6h01npWJzZrEWJLYwNYBqwM3+ecxsiLSUET+KSK9C0ogIjc4ceq5JlUp4ocffoi1CDGF5e9t/hw7Ro1TxzgZV4lqjb3pXez5OoDVgZv8g3uyDwAXAO8VkmYJ5hSe+90SqrQgIjz88MOe9i4+evRo0ZEqMLzOH+fPZXfV2tSqmhhjYWIDz9cBvK0Dn8/Hww8/7NopPBLotiwia4CXVXVioYlEHgRGqGqaK1KVAkREVZUjR45Qo0aNWIsTM1j+3uavy5YhV17JyvptuHDnBioneuvAdrB1AKwO/PxFBFWNqrUN7sm2AlaHkS4DSI2mILGCl9eHgeXvdf5nvje7PWVVS/akgQVbB8DqwE3+wUa2OKtx3Vm5W8qwruuWv5dxKvN7APbVrBAuFiWC1+sAWB2U5hKe7ZiTa4rCxZiTdco9vDxEApa/1/mf3bkLgMO1vWtkvV4HwOrATf7BRnYh8ICIFGjWHa/iBwjzLNiyjszMzFiLEFNY/pmxFiGm0F3GyB716OEAYOsAWB24yT/YyD6LOQlnpYj0F5HK/hsiUllE+gOfY4aKn3VNqlKCiHDrrbd62ru4ffv2sRYhpvA6f9m9G4ATdb25fAdsHQBv68Dn83Hrrbe65l2cx8iq6j6gF5AN/BM4LCK7RGQXcBiY59zrpar7XZGoFGEPCLBvsF7nX2nPjwCcbuDNwwHA1gHwtg7cPiAg31F3qrpBRC7AnP96FdDUufU98CHwlqpmuyJNDHDixIlYixBTWP4e5q9KZcfInm3oXSPr6TrgwOs6cJO/uGW9yzrsOlkDy9/D/A8cgDp1OFypMk/O+JQJQ715pqin64ADr+ugNNfJ5sKZg71MRAaJyEARuTRwjraiwK4Ps/w9C2e3px+rJHt2tyfweB1w4HUduMk/33CxiCQBTwN3A0kYRygwzk4nRWQK8EdVPe2aVKWIBg2861UJlr+n+fu3VKySTM0qlWIsTOzg6TrgwOs6cJN/HiMrxr1qIdATmA8sAnZgDG1ToA/wINABuME1qUoRiYnefYMHy9/T/J3lOz9UTaZmVe8aWU/XAQde14Gb/IOHiwcCVwIDVfVmVZ2uqu+r6hLnez9gENBLRG52TapSgohw9913e9q7eOfOnbEWIabwNH//cHFVbw8Xe7oOOPCyDnw+H3fffbdrS3iCh4uHAG+q6tsFJVDVt0RkLjAUeCucQkSkKfAX4BpMr/hD4H9VdUcYaQvyzLpQVTMC4sUBfwDuARoCm4DHVPWfhXDhp59+ol497+52c8EFF8RahJjC0/wDhosv8nBP1tN1wIGXdeDz+XLtgBuGNrgneyHwbhjpFgIXhVOAiFQFPgbaASOAYZjDBZaKSLUw5XwZ6B50bQ6K82fAB7wAXA+sAuaKSKHD2ps3B2fjLVj+Huafe8xdMjWreLcn6+k64MDrOnCTf3BPtj5mDrYo7ADCnSm+G2gJtFXVrZB7pN4WTK+z0GP1HOxS1VUF3RSRBsDvgAmq6t+JaqmItAYmYOaWQyI7u8Is+S0RLH8P8w90fPJwT9bTdcCB13XgJv/gnmxV4FQY6U4D4S7n6Qus8htYAFXdDnwG3BRmHkXhWiARmBkUPhPoKCItCkrYsWPHKIlQPmH5e5i/4/i0u2oytTxsZD1dBxx4XQdu8g+1TvZ8EWlZ2AU0KUYZPwPWhQhfj/FSDgf3isgpETkuIh+LyOUhyjgFbA0KX+98FlhORkZGQbc8AcvfHf47d+7kvvvuo3v37lStWhURKfbWdZmZmfh8PrZt2xZ9AXNyUGff4h+r1KaWh4eLvd4GwOrATf6hjOw8zFBuYdfcYpRRBzgQInw/kBxG+pnASOBq4NdAXeBjEUkPKuOg5t++an/A/ZBo1KhRGCJUXFj+7vDfunUrb775JsnJyVx+efA7YXjIzMxk3Lhx7hjZvXuR7Gz2JVVHkhKplFDgvjQVHl5vA2B14Cb/4DnZO10qJ5SHcFhuXKo6LODnChGZj+kZPw5cFpBXscsQEdq0aZNn0nvs2LGeXtJjER1cccUVZGVlATB9+nTef//9GEsUhID52OpJ8TEWxsKidOHz+Rg3blzu72A7EE0En8LzSnGuMMs4QOieZDKhe7iFQlWPYDyguwQE7weSJb//dXLA/VB5MXXqVFQ19/Kagd3tDBl6FW7xj4srumf4448/MmLECBo3bkxSUhKNGjWiT58+7Nmzh2XLlnHllVcCcM011yAiiAjLli2LjoABG1FUybfvm7fg9TYA3tOBz+fL87/vtwNuoDSa13rMnGkwOgAbSphncM91PWYLyFbknZf1z8UWWE5aWloJRagYsPxjx3/YsGF89913PPPMMzRt2pSsrCw++ugjjh8/zkUXXcSkSZP4zW9+w9/+9je6dDHvlB06hOvGUAQC9i2uX7t6dPIsp/B6GwCrAzf5l8ZEzL+Abo7DFAAikgJc6twrFkSkJtAb+HdA8GKMx/PQoOi3A+scb+aQWLt2bXFFqFCw/GPHf+XKlfz2t79l6NChXHHFFQwaNIgXX3yRlJQUatasmWtQ27dvT7du3ejWrRs1a9aMTuEBa2Tl7Mno5FlO4fU2AFYHbvIvjZ7sNOC3wHwReRTTA/0z5nzaqf5IItIc+BazS9NjTtjvgLbAUuAHoDlmPWxDAgyqqu4Rkb8AY0TkCPBf4FbMHsyFLhOKj/f2fJTlHzv+Xbp04ZlnnkFV6dmzJxdccIFrW7vlQ8CcbLUk7zo9gW0DYHXgJn/XW5eqHsMYu83Aa8AsYDvQU1WPBkQVID5Ipk2YId+/AR9gNq7YDlymqiuCinoE4wz1ALAE01O+RVUXFCZfmzZtSkasgsDyjx3/OXPm0LdvX55++mk6derE+eefz2OPPUZOTo77hQeskT2/QThO/hUXXm8DYHXgJv9SeYVV1R2qOkBVa6pqDVXtp6qZQXEyVVVU1RcQtkBVL1XVeqpaSVXrqmpfVf0iRBnZqvq4qjZX1SRV7aSq84qSbd26UEt4vQPLP3b8GzRowKRJk9i1axcbN27kjjvuYOzYsUydOrXoxJEiYLj41NGD7pdXhuH1NgBWB27y9/Q4kYhw1113ec6jOBBNmhRnX5GKh7LCv23btjz55JMkJyfnNvikpCQATpw4Ef0CA4aLG9WrHf38yxHKSh2IJbysA5/Px1133VVqp/B4CqrKhg0bouexWQ5x+vTpWIsQU7jJf948M5Dy1VdfAfDee+9Rv3596tevT1paGldffTVDhw6lXbt2VKpUifnz53PgwAF69eoFmCGshIQE/vGPf1CnTh2SkpJo27YtNWrUiEywM2dgzx5yJI69lWt6fgmP19sAeFsHPp+PW265hQ4dOrhiaD3evGDPnj2eNrKWv3v8Bw0alOf3yJEjAejRowdLlizhoosuYtq0aXz33XfExcXRtm1bZs2axU03GV+9unXr8sILL/DUU0/Ro0cPsrOzWbp0Kenp6ZEJ5uwgdSC5Ptlx8eScPhZZfuUcXm8DYHXgJn9xawFuWYeIqKpy5MiRyHsG5RiWvwf5v/MO9O9PRtsuXN71f5lxTxcGXpoaa6liBk/WgSB4XQd+/iKCqka1O+vpOVk4N5TnVVj+HuS/cSMA2+qYebjdO1zYG7kcwZN1IAhe14Gb/D1vZKtUqRJrEWIKy9+D/L/5BoBNtc8HoE4ND+ogAJ6sA0Hwug7c5O95I5uSkhJrEWIKyz8l1iKUPhwju6G6OXmkXatmsZQm5vBkHQiC13XgJn9PG1kRYfDgwZ5ewvON84frVXiOv2rucHFGlfMAO1zsuToQAl7Wgc/nY/DgwXYJjxtQVTIzMz39Fudl7uBB/j/8AEeOoHXr8h1miKx9aosYCxVbeK4OhICXdeDz+XLtgBuG1tM9WTBeZV6G5e8u/4yMDHw+H/v3hzxtsfTh9Fiy27ZDFapXTuD4saNFJKrY8HobAKsDN/l73sju27cv1iLEFJa/u/wzMjIYN25cmTOyJ1uZJTu1qibaOuBx/mB14CZ/zxvZzp07x1qEmMLy9xh/x8gebt4agNrVEr2ngyB4nT9YHbjJ3/NG1q4Ps/wjxebNm+nfvz8NGjSgcuXKNGvWjEGDBjF9+nTuvPNOAFJTUxERRITMzEwAzp49y/jx42nXrh1JSUk0btyYhx56iJMnz53vmpmZiYgwefJkRo0aRYMGDahatSp9+vTJzadYcJye9jUx87C1qybaOuBx/mB14CZ/Tzs+iQidO3emT58+nvUwrl69eqxFiCmiwb9Pnz7Url2bKVOmUK9ePXbt2sWiRYu48cYbefTRR3n88ceZO3du7ibsjRqZpTO33347CxYs4A9/+AO/+MUv+Oabb/i///s/MjMz+ec//5mnjPHjx5OWlsaMGTPYs2cPf/zjH+nVqxfr16+nUqVK4Qvr9GR/bJQCZFKraiWqV68WsQ7KM7zeBsDbOvD5fCxcuNA9Q6uqnrwMddVdu3apl2H5R8Z/7969Cuj8+fND3p8xY4YCumXLljzhn3zyiQL6yiuv5AmfOXOmArp69WpVVd2+fbsC2r59e83Ozs6N9+mnnyqg06dPD1/YzZtVQbVWLX3l481afdgs/fXUz20d8Dh/VasDP3/HLkTV1nh+uHjz5s2xFiGmsPwj41+3bl1atmzJ6NGjmTZtGlu2bAkr3eLFi0lMTGTAgAGcPXs29/KfwPPJJ5/kiT9w4EDi4s4110svvZQmTZqwcuXK8IWdO9d89u3LwZNnAUiulmjrgMf5g9WBm/w9b2RbtWoVaxFiCss/Mv4iwgcffMDFF1/MmDFjaNOmDS1btmTKlCmFptuzZw+nT5+mevXqVKpUKfdq0KABkN/b8bzzzsuXx3nnnceuXbvCF9ZvZAcO5NAxc7RZ7WqJtg54nD9YHbjJ39NzsmD+zJo2bRprMWIGyz9y/i1btuTVV19FVfn666954YUXGDlyZKEL/OvWrUvlypVZsWJFyPuNGzfO8zsrKytfnKysLNLS0sITcutWyMiAGjWgVy8OvrkWOLeEx9YB7/IHqwM3+Xu+J3vw4MFYixBTWP7R4y8ipKWlMXHiRADWrVtHUlISACdOnMgT97rrruPkyZMcOnSIiy++ON8VbGTnzZtHTk5O7u/PPvuMnTt30r179/CEcw6Qp29fqFyZgwE9WVsHvM0frA7c5O/5nqxdH2b5R4I1a9bwwAMPcOutt9K6dWuys7N5+eWXSUhIoGfPniQkmCY2adIkRowYQaVKlejUqRPp6ekMGTKEgQMHMmrUKLp27UpcXByZmZksWrSIp556ijZt2uSWc+TIEfr168c999zD3r17GTNmDKmpqQwfPrxoIXNy4PXXzfeBAwHOGdmqiXRubeuA1+F1HbjKP9qeVOXlAhTQ9PR0HTt2bPhuaBUMS5cujbUIMUWk/LOysnT48OGampqqVapU0eTkZL3iiit08eLFuXF8Pp82btxY4+LiFNDt27erqmp2drY+99xz2qlTJ01KStKaNWtqp06d9OGHH9aDBw+q6jnv4kmTJumDDz6o9erV0ypVqugNN9yg27ZtC0/IceOMV3H9+qrHj6uqas9xS7T6sFn62cYsWwc8zl/V2zoYO3aspqenq98maLRtTbQzDFkINAXmAYeAw8BbQLMw0l0MvARsBI4DO4BZQIsQcTP9Sgq6+hWQt6pq7lIJr8LyL9v8/UZ22rRpJcvgvfdURcy1ZElu8EW/X6DVh83SDd8fKPM6cBte569qdeDn74aRdX24WESqAh8Dp4ARjuF7HFgqIp1U9VghyQcDPwP+BqwHzgf+D/hSRNJU9fug+EsAX1DYpsLkq1u3bphMKiYs/wrKf/duePxxmDrVHG/32GPgLA8COHj83JzsWa2gOggTFbYOFANe14Gb/EvD8eluoCWmR/mOqs4H+gLNgXuKSPuUql6qqpNVdbmqzgauA5KdfIPxk6quCroOFFbAt99+W3xGFQjR4j9v3jwGDBhA8+bNqVKlCm3btmXMmDFRPd3i5ZdfzrMtYTRQIZ7/iROwbh3MmQN/+AN07QqNG8PkycbAPvAAPPJIbnRVzZ2TrVU1sWLoIAJ4nT9YHbjJvzQcn/oCq1R1qz9AVbeLyGfATcDEghKq6t4QYd+JyF5MrzZiBDqXeBHR4v/ss8/SrFkznnzySZo0acLq1avx+XwsXbqUzz//PM9GCmUJZf35p6SkmOmNM2dg82bYsiX/5/ffG2MaiKQkuP56+POf4YIL8tw6eSab02dzSEyIo0pifJnXgdvwOn+wOnCTf2kY2Z8B80OErwcGFTczEWkPNAC+CXH7RhE5DsQDq4EJqvpOYfn98MMP+ZZLRIIlS5bw1FNPsWHDBg4cOED9+vX5xS9+gc/no0OHDiXKMzMzkxYtWjBjxgzuuOOOqMkK0eO/YMEC6tevn/u7R48e1KlThxEjRrBs2TJ69uwZcRluINrPPyLk5BiDGcqYbt8O2dmh0yUkQMuW0LYtXHQRdOkC6elQLfSexPuOnAKgTvUkRKRs6SAG8Dp/sDpwk39pGNk6QKgh2/2YYd+wISIJwIvAXuDvQbcXAP8BtgPnAb8F3haRYao6s6A8jx6N7oHV+/fvp3PnzowcOZL69euzY8cOJkyYQLdu3Vi7di3NmzePanmRIlr8Aw2sH126dAHI3ZXo5Zdf5s4772TlypU8//zzLFiwgOrVqzNw4ECefvppKleunJt227Zt3HfffSxdupRq1apx2223lfglpTBE+/kXCVX48cfQhvTbb+HUqdDpRKB5c2jTBlJT8342bw7FOCTgJ8fI1qth1vCWug7KGLzOH6wO3ORfWutkNUSYlCCfF4BfAL2D51pV9b48mYu8DawCxgMhjayIUKNGjTzzhmPHjo3oRJ4hQ4YwZMiQPGFdu3alXbt2zJs3j4ceeihkulOnTuVuXFCacHN92PLlywFo3759nvBhw4YxZMgQ3nrrLVauXInP5yM5OZlx48YBcPr0aa655hpOnDjBpEmTaNDghFNNAAAgAElEQVSgAVOnTuWtt96Kuoyu8d+/3xjOYGO6ZQsU1qAbNsxrQP3fW7WCgJeQSPDTYXOUXv2aJj+7RtLb/MF7OvD5fLn/N0A+OxBVRNtdOfgCsoCpIcInA3uLkc94IAcYVow0v8cY+EYh7qlq6awP85/U8te//lVVz53Msnz5ch04cKDWqlVLf/7zn6uq6rFjx/Tee+/VOnXqaLVq1fTGG2/UFStWKKAzZsyIumxu8d+5c6fWr19fr7766twwP+8//elPeeL27t1bU1NTc3+/9NJLCujKlStzw7Kzs7VDhw551plGAxHxP3JE9b//VX3jDdU//1l12DDVbt1U69Qxq+MKupKTVS+5xMR/7DGT/quvVA8fjhqvwvD6p9u0+rBZeuekT1XV22skVS1/VasDP3/K4xIezNzrz0KEdwA2hJOBiDwCjAbuV9XXilG2v7ccqicNuOe6nZ2dTXZ2Nt999x2jR4+mYcOGDB48OE+coUOHMmTIEObNm8fZs+ZUlHvuuYc5c+YwduxYunTpwgcffMBtt93miozgDv+jR49y0003kZCQwIwZM/Ld7927d57fHTt25MMPP8z9vXLlSpo2bUq3bt1yw+Li4rjllluifu5vWPyzs2HTJli92uz/u3o1bNhglskUhGrV8g/t+r/HeLlE7nBxTTNyYpdveJs/WB24yb80jOy/gGdFpKWqbgMQkRTgUozhLBQicj9mXe0jqvp8uIU687eDgB2q+mNB8WrUqBFulsXCJZdcknsIcOvWrfn4449zT1jxwz8X6cemTZuYPXs2TzzxBKNHG9X06tWLo0eP8uKLL7oiZ7T5nzx5kr59+7Jt2zaWL1+ee1B5IOrUqZPnd1JSEqcC5iJ3795d4Kkz0Uah/L/4Al56CRYsgD178t9PSjLDuKHmSRs2NPOoZRA/HTHDxfVqmOFit9pAeYHX+YPVgZv8S2NdxTTMbkzzReQmEemL8Tb+HpjqjyQizUXkrIj8KSBsMPAcsBj4WES6BVwdAuINEZE3RGS4iFzppFsKdAb+UJhw0VxzGYjXXnuNVatWMXv2bGrWrMk111yTr6z+/fvn+f3vf/+bnJwcbrnlljzhwT3gaCKa/M+cOcOAAQP44osvWLRoER07dixRPo0aNSrw1JloIyT/H36AIUPgkkvg7383BrZJE+jXD3w+mD8ftm2DY8dg/Xp4+214+mm46y7o0QMaNSqzBhbgp8N5HZ/cagPlBV7nD1YHbvJ3vSerqsdEpCfwF+A1zBDuR8D/qmqgB4hglt4EGv7rnPDrnCsQy4F05/t2zLKeZzDezMcxnsbXqeqSwuQLdsqJFvz5XnLJJVx//fWkpKQwYcKEPD3SRo0a5Umz2xl+DO6xudGDC5YzUuTk5DB06FA++ugj3n333TxDvcVF9+7dmTFjBqtWrcrNJycnhzfffDMqsgYiH/89e+CKK4ynb1IS3HcfDB9u1pqWYcNZHPiHi+s6PVm32kB5gdf5g9WBm/xLxbtYVXcAA4qIk0mQx7Gq3gHcEUb+q4ASLcTMzMx01YgB1K5dm9atW7N169Y84RL0p+03ullZWbRs2TI33I0enB/R4v+b3/yGuXPn8sgjj1CtWjVWrVqVe69JkyYhh40LwogRI5gwYQI333wzTz75JA0aNODFF1/k8OHDEcsZjDz8jx2DPn2Mgb3wQtNDLWNLrqKBvY53sX9OtjTaQFmG1/mD1YGb/MvmNjylBBFh9OjRUXemCUZWVhYbN26kVatWhca75JJLiIuLy9dje+ONN1yTLfic05LivffeA+CJJ56ge/fuea7p06cXK6/ExEQ++OAD0tLSGDlyJCNGjKBFixY8+uijUZE1EHn4P/QQ/Oc/kJICixZVSAMLsC93Tjb0Wbdeg9f5g7d14PP5GD16dL5OT7QgGrwdm0cgIqqqHDlyJKqT3v379+eiiy6iU6dO1KxZk82bN/OXv/yFH3/8kS+++II2bdrkbsqwZcsWWrdunSf9sGHDmDNnDj6fL9e7eM6cOezYscOVHZ+izb+8IZf/Dz9AixZm+8I1a/JtRViR0OR/5nLo+Bm+mzyAOtWTbB3wOH+wOvDzFxFUNarW1tM9WSDXAzha6NatG++88w4jRoygd+/eTJw4kR49epCRkRHW/phTp07lV7/6Fc8++yz9+/dn48aNzJ49O6oyBiLa/Msbcvk/9xycPg0331yhDezps9kcOn6G+DihdtVEwNYBr/MHqwM3+Xu+J7thwwZXtusrL7D8N9ChcWNo1gyOHDHLdpztICsiMvcepeND/+L8OlXZ+Fw/wNYBr/MHqwM/f9uTdQGJiYmxFiGmsPwTYdo0Y2CvvLJCG1iAnfuOA3B+naq5YbYOeJs/WB24yd/zRnbnzp2xFiGmsPx3wuuvmx8PPBBbYUoBu/YfA6BJgJG1dcDb/MHqwE3+njeyF1Tg+bdw4Db/jIwMfD4f+/fvd7WckqJTrVpmm8Rq1eDaa2MtjuvI7cnWPWdkbRvwNn+wOnCTv6eNrIjQt29f15fwlGVs3rzZ1fwzMjIYN25cmTWyB1991Xy5/vqonXJTlrFrvzGygT1Zt+tAWYfX+YO3deDz+ejbt69rS3hK66i7MglVZcWKFVx++eWu5J+RkcE777zD/fffn2+/3rKC7IIOAvcI6jjH8dGvX2wFKSXs3O+fkz13oLvX64DX+YO3deDz+bjqqqu4/PLLXTG0nu7JAiXeXzcclPVeHESH/+bNm+nfvz8NGjSgcuXKNGvWjEGDBjF9+nTuvPNOAFJTUxERRCR3n9CzZ88yfvx42rVrR1JSEo0bN+ahhx7i5MmTuXlnZmYiIkyePJlRo0bRoEEDqlatSp8+fSLfb3TfPmqtWQMJCXDDDZHlVU6wyxkubhowXOxmGygP8Dp/sDpwk7/njWxGRkasRYgposG/T58+7Nq1iylTprBkyRImTJhAUlISN954Y+4uTXPnzmXlypWsXLkyd/vI22+/nccff5zbbruNd999lzFjxvD3v/+doUOH5itj/PjxbNmyhRkzZjBp0iS++uorevXqxZkzZ0ou+MKFSHa28SpOTi55PuUI3+9zHJ8CjKxtA97mD1YHrvKP9gG15eXCObR948aNoc7wDRubNm3Sfv36af369TUpKUmbNm2qAwcO1GnTpinmHNs8l//A8TNnzuiTTz6pbdu21cTERG3UqJGOGjVKT5w4kZv39u3bFdBJkybpgw8+qPXr19cqVapo7969o3ZweaT8/QfSz58/P+R9/0HtW7ZsyRP+ySefKKCvvPJKnvCZM2cqoKtXr1bVczpo3769Zmdn58b79NNPFdDp06eXXPh+/cwh6pMmlTyPcoSDx05p9WGztN4v39Ds7Jzc8EjrQHmH1/mrWh34+ePCoe2e78lGinLbi4sS6tatS8uWLRk9ejTTpk1jy5YtYaVbvHgxiYmJDBgwgLNnz+ZevXr1AuCTTz7JE3/gwIHExZ2rrpdeeilNmjRh5cqVJRP8+HFY4hzQ1LdvyfIoZ9i46xAAbRrVJC6uYpwoZGFR1uFpxycRIT09nR49epTIw/inn35iy5YtzJ8/n74Bf9S33XYbQO6BAGlpaXn2KF6xYgVz5szhlVdeYfjw4QBcffXV1KlTh9tvv52MjAzS0tJy49eoUYP58+fnGpk2bdpw2WWX8eqrr/KrX/2q2HIHYvfu3bRt27bE6UWEDz74AJ/Px5gxY9i3bx8tWrTg4Ycf5t577y0w3Z49ezh9+jTVq1cPeX/fvn15fhd0iPuuXbtKJvj778OJExxu146axTghqDzjG8fItm9SK094pHWgvMPr/MHbOvD5fCxfvpxly5a5kr+njayqcvDgQWrXrl2i9IG9uKysLNLT00lNTS0yXXAvzo/AXlygkS2sFxepkQ0sp6Ro2bIlr776KqrK119/zQsvvMDIkSNJSUkpME3dunWpXLkyK1asCHm/cePGeX4XdIh7ieV/5x0AKg0aVLL05RD+nmy78/Ma2WjUgfIMr/MHb+vA5/Pl2gHrXewC1q5dW+K0/l7cxRdfzJgxY2jTpg0tW7ZkypQphaYL7MVVqlQp92rQoAFQCr24AETCPxgiQlpaGhMnTgRg3bp1JCWFPk7tuuuu4+TJkxw6dIiLL7443xVsZOfNm0dOTk7u788++4ydO3fSvXv34gt69iwsWADA+jBeiioKcnuyQUY2mnWgPMLr/MHqwE3+nu7JAsTHx0eUvlz24gIQKf81a9bwwAMPcOutt9K6dWuys7N5+eWXSUhIoGfPniQkmCo2adIkRowYQaVKlejUqRPp6ekMGTKEgQMHMmrUKLp27UpcXByZmZksWrSIp556Ks+pRUeOHKFfv37cc8897N27lzFjxpCampo73F4sfPop7N8PbdpwumXLiPiXJ2wswMhGWgfKO7zOH6wOXOUfbU+q8nLheBfv3bs3PPezMHHo0CEF9Omnn9bZs2croGvWrMkTZ+nSpQrohx9+WGhernrWOoiUf1ZWlg4fPlxTU1O1SpUqmpycrFdccYUuXrw4N47P59PGjRtrXFxcHg/r7Oxsfe6557RTp06alJSkNWvW1E6dOunDDz+sBw8eVNX8Htb16tXTKlWq6A033KDbtm0rmdBDhxqv4t//PurPv6zi+5+OavVhs7Txr9/UswF1STX6baC8wev8Va0O/Pxxwbs45sYuVpffyC5dujTc55APX3/9taanp+uUKVP0gw8+0MWLF+vgwYM1ISFBv/zyS83IyFBA77nnHv3888/1P//5j546dUpVVYcMGaK1a9fWxx57TBcvXqzvv/++vvTSS9qvXz/dtGmTqp4zME2aNNEbb7xRFy5cqDNmzNCGDRtqamqqnj59usSy+xEJ/9KAXwfTpk2LToabN6vGxakmJKhmZpZ5/tHCzE++1erDZumgicvy3fOKDgqC1/mrWh34+bthZD0/XNwkAs/Shg0b0qxZMyZOnMjOnTupXLkyHTt2ZOHChXTu3Bkwk+ovvfQS06ZNIycnh+3bt5OSksLMmTN5/vnn+cc//sETTzxBUlISKSkpXHvttfnmYMeMGcPWrVu54447OHbsGFdeeSUvvPAClSpViog7RMa/XGL8eMjJgTvvhObNaVIGlkGVBpZvMFMOPdrnn9/3XB0Igtf5g9WBq/yjbbXLy4WzOUT79u117Nix4b/ylCKi3osLgfXr17uWdzQQVR18+61qfLy5tm5V1bLPPxrIzs7R1Pvf0urDZuna7/bnu+8FHRQGr/NX9bYOxo4dq+3bt8/dMEijbGtKxbtYRJqKyDwROSQih0XkLRFpFmbayiLyjIjsFpETIrJSRK4IES9ORMaISKaInBSRr0VkQGF5qyqTJ08ul6fwREvmPXv2RCWfkiAcDikpKagqd911V2SFff89XHcdZGfD0KHgrGGOlH95qDsfrt3N7gMnaF6vGh2a5F+u9uabb8ZAqugikucQyzYQiFjWpWjooDy0hVDw+XxMnjzZ3/mKPqJttYMvoCqwBVgH9ANuAtYC3wLVwkg/CzgI3A1cBbwFnADSguI9AZwCfgdcCUwFcoAbCshXVVUPHz4c4XuQeyisF+eXP1LEkn+0OBSKU6dU585VTUlRBdW0NNUAJ49I+ZcKhwgx8P8t1erDZumz/1oX8n6NGjVKWaLoI5LnUFb+A2JZl6Khg/LQFgqCnz/ldE72bqAl0FZVtwKIyBqM4b0HmFhQQhH5OXAb8EtVneGELQfWA48BfZ2wBhjjOkFVn3WSLxWR1sAEYFFBZXz11Vekp6dHws81+HtxbqJQ/qqm55edbdaW+r8HX6Hu5eSEjhsQfh3AwoVhxw8rPDsbjh6FAwdg+3ZYtw4OHjR8unaFxYvzHAZQlp9/NLA4YxeLM36gcqV4hvdoFTKO33/Aq6jodSAceF0HrvKPttUOvoCPgM9ChC8HlheR9v+A00DVoPBxmF5rkvN7GGY8PTUo3p1OeIsQeauq6qpVqwp8uwl3rjYW8QjzrbHQvE6d0rX/+Ifq88/rF507q151lWrHjqrNmqnWrm28cI2pLd9Xp06qf/2r6rFj+XRS0PMP91lF5TkUM164ef3q9xO06f/M1erDZulfF20oMF7Xrl1LXbZoxwvnORSUV3AdKI8cIo1X2P9guPnFoi2EG6+oOH7+uNCTFXW5pyQiPwLzVfWeoPDJwCBVrV9I2jeAC1W1bVD4LcAc4AJVXS8iE4D/BapoACER6Qr8G+ijqu8G5aG712xm//6fqJtcJ99fs6py+WWXseKTTxB/loH3cb4DV191FR9++GHub1TNQLX/u/N5/Q3X8967i5yffjHNZ2AZffv25V/z55s4gXkGfA4cMIB58+blv0+grDBkyGBen/0GcuokcQcPEp+1m0qZ20lct4bENRnEnTpVkPpNNnFxEB+PxsdDfAKHjx2jRnJtiPOHxaPxCRAfZ37HxbN527ektm1nwuKcOLmfcRAfB3HxfPL551yWng5x8RAXh8bH5Y0bHw8Sx7y332LALbc65Tl5Bsjl/z1pyhRG3n8/OdWqkVOjJmeaNuN0i9acbdQYArZL6/Czn7Fh/XoAftr3E3Xr1M3H+4ILLmDdunXn9FCAfjp27MjatWspqhl16tSJr7/+umA9O59paWkhj90KzP/CCy9k9erV+eKcPJPNnkMn2b7nCB+v+5EP1+4G4IYLz2f2A5cTHxfaBeO8884LueFJMESEov4vwokTq3gFxcnKysrj0V8eOUQaL1gHZUm2aMQrKo6fvxMvqnsrloaRPQ1MVNXRQeGPA6NVtcAhaxF5H6ipqt2Cwq8GPgCuUNUVIvIS0FdVGwbFa40Zlh6uqq8F3XOZefnA5pqN+LJea9YkN2dLrUbsrlKbw5WqcjixKkcqVeZsnOdXeZVbaPZpftcvjUdv7kRCfME+jldeeSVLly4tMr+K+ke7bNmyPEOF5ZFDpPGCdVCWZItGvKLi+Pm7YWRL6x80FLtwiEiYacONV2wBKjwO7zYXobd4tCjfGDfbzK0UhXA3Rg8nXjTzina8sixbuPGsbCWL58bm/+GgNIzsAaBOiPBk515h2A+EWuqTHHDf/5ksIqJ5X1eC4+Ui2m8rFhYWFhYWwSiNdbLrgZ+FCO8AbAgjbQsRqRoi7Wlga0C8JCDYfbKD81lUORYWFhYWFlFHaRjZfwHdRCT3uBMRSQEude4VlbYSkHvop4gkALcC76uq32tnMcboDg1KfzuwTlW3RyC/hYWFhYVFyRBtd+XgC6iG6XGuxWxE0Rf4GtgGVA+I1xw4C/wpKP0bmGHluzCbUcwDTgLXO98PAYeBTU74KCAdmILx8b0xhEyVgWeA3ZiNLVZinKiC48UBY4BMJ++vgQFR1E3TIA5vAc3CSHcx8BKwETgO7MBs2hFqqVImznZhQVe/WMnvpA0lk5J/k5Gy+gx8hXA4WVrPwMm/CfC8U4+PO3mnhJm2rLSFEnEoK20hCs8h5u0hgmdQltrCQOCfwHdOfd4EjAdqxKotRKWBhCF8M4f4YeAI8E7wwwNSHGX7gsKrYDas+NEh9G/gWkLvIrXPaWSngDXAwALkcWUXqWLqpMQ7YQHPAp8BI4EemA07vnH4Nw1RqRcD3YKu5FjJ76RXYEYIuYLXRJfVZ9AkhOxXAWeAN0vjGQTknw5kYTZdWULx/txj3hYi4VAW2kKUnkNZaA8lfQZlqS2sAt7EjGr2wCztPOiEx8WiLURMKhYX8ACQDbQOCGuB6QmPKiLtz53Kc2dAWALmjedfAWENHEWOC0r/EbAmxhzqhwhr7jzox0JU6pll6Rk4cRV4vIg4ZfYZFJCff1OU3qXxDALyjwv4flcx/hzLRFuIkEPM20KkHJz4ZaE9lFj+EHnFqi2Eqg/DHVl6FpLOtbZQKgcEuIC+wCp1tmkEUDPv+hmmN1JU2jOYzSz8ac9ihqWvFZEkJ/haIBGYGZR+JtBRRFpExCACDqq6N0TYd8Be4PwI5QoXkTyDcFFmn0EBGIHpCSyJUK5iQVVzio4VEmWlLZSYQxlpC/5yS/ocwoWrzyHK8seqLeSrD8B/nM/C6oNrbaG8GtmfYYb4grGecx7FhaXdrqrHQ6RNBFoHxDvFOQ/mwHiEUU5RiIRDPohIe8xb1jchbt8oIsdF5JSIrBKRfsXNPwSiIf+9jkzHReRjEbk8RBnl4hmISBPM0NEsp3EGw41nECnKSluIKmLQFqKFWLeHqKAMtoUezmeo+uCHa22hvBrZOoReY7ufc2tjS5LWf9//eVCdsYBC4pUUkXDIA8fj+kXM2/vfg24vAO7DvIENxcxrvy0itxdX4CBEKv9MzDza1cCvgbrAxyKSHlRGuXgGmOGxOOCVEPfcegaRoqy0haghRm0hGigL7SFaKDNtQUTOxxwm86GqfllIVNfaQnneMy+YJJSBXaSKiWjl/wLwC8z8R56Koqr35clc5G2ME8B48g95FBclll9VhwX8XCEi8zG9yseBywLyKi/PYDiwWlXX5CvA3WcQCcpSW4gWYtUWIkIZag/RQJloCyJSHZiP8bG4s6jouNQWymtPNtJdpApK67/v/0yW/HtxFbiLVDERCYdciMh4zJvvL1X1/aLiq2o2MBdoIiKNwi0nBKIif4BcR4B3gS4BweXlGXQF2hH6zT0fovgMIkVZaQtRQQzbQtQRo/YQMcpKWxCRyph9FloC16rqziKSuNYWyquRrQi7SEXCAQAReQQYDTygQQcgFJXU+Qz1RhYuIpY/BILfEsv8M3AwAvO2PLsYaaLxDCJFWWkLESPGbcEtlHZ7iAZi3hZEpBJmyWhXzLKatWEkc68tuOVK7eaFWft0FmgZEJaC8Q57qIi0aZiHOSLIVfsbYEEIV+2xQek/BNbGkoMT936Hxx+LWW4C8CXwXSzlD5FfTcwa5+UBYWX6GTjxEzFrMueX9jMoIO/iLH8pE20hEg5O/Ji2hWhwCJG+1NtDpPKXhbaA6Ti+iZnrvaoY6VxrC1GtWKV14d4uUhcFxZtAmLtIlSYHYLAjx3vkX9TdISDeEIfrcIy332DMcTsKDI6h/L8DpmE2DkjHvP2uxbwxXl4enkHAvZsdfd5cQBmuPYOgcgY61xQn73ud3z3KeluIhENZaAtR4FAm2kMk9aistIUAuR8PUR+axKItRKVixeIi+rtIpYcoIx54FLNFV6G7SJUmB+BlCt7GbFlAvG7Ax5j1amcwWwd+iJmjiKX8N2LWov7kyLUPM3/Stbw8g4B78x35EwvI39VnEFBOofWhrLeFknIoK20hQg5lqT2UqB6VlbZAwds25spc2m3B9UPbLSwsLCwsvIry6vhkYWFhYWFR5mGNrIWFhYWFhUuwRtbCwsLCwsIlWCNrYWFhYWHhEqyRtbCwsLCwcAnWyFpYWFhYWLgEa2QtKhRE5A4R0YArW0R2icibItI2KK5PRMrlGrZg2UWkthN2USzlChci0t15Pk1C3Ovv3LulkPSvichREanmrqQWFpHBGlmLiopBQHfgCmAMcCHwkYjUCogz3YlTEVAbGAuUCyML9AO+1NAbty/EbGowLMQ9/+kq/YF/quox90S0sIgc5fmoOwuLwpChqv5NvT8TkR+ADzDHoL0H4PzBF3U6h4U7uAkIuZG/qp4RkdeB/xGR+qq6NyjKAMyWmGGd9GJhEUvYnqyFV3DY+azkDwg1XCwivxWRlSKyX0QOisgqEekdFCdBRP4sIt+KyEkR+UlEPhWRyygAIjJZRLKcQ8UDw5NE5ICIPBcQ1lZE3nbKP+HIcF0heacA252f0wKGyu9w7vcSkUUisltEjovIOhF5SETig/KpKiJTRGSfiBxxZPhFYF4BcXuIyEdOvGMiskRELihIxqC07YC2mC0sC8IrmE7AkBD3hgPfA8sC8mzpcDwuIntE5P+JyK8d2VMC4g0WkY9FZK8z3LxaREaEkFFF5AkReUREdjrP4RMRSQuHo4WFH9bIWlRUxDvGMElE2gNPAnsI+GMuACmYYeRBwK2YE0IWisj1AXH+ADwI/A24FnMg9EeEPo/Sj1cxJ3j0CgrvgxnqfQ1ARBoDnwI/B34L3AIcBN4NkiEQuzGbs4M5/Lq7c73rhLV05Psl0BtjwHzAE0H5vOTEedbJbxMwK7gw56XjI+AocDtmY/samMPGmxaogXPoD2xV1fUFRVDVLzHHiuUZMnbmcNOB11Q1xwlLxIxS/BwYCdwBtAAeCZF1S8zG70MxQ9YLgOki8j8h4g4HbsA8hzuA8zBTDoU9ZwuLvIjm5tj2slesL8yfYajNwXcBXYLi+kwTKDCvOExv6n0Cju/CzBm+VQLZNgOvB4W9A2wI+P0s5oSQ1gFh8RiD99+CZOfcpud3FSGDOJwewZw4EueEt8WcJPL7oPh/c/K9IyBsK/BRULyamA3unwtDD/8Gngkj3u+dstsHhI12wtoGhP3aCesaxPNrCjmuLeD5TgO+DrqnDp9qQTo+A/w51vXcXuXnsj1Zi4qK/kAXzMHN/TCHKS9yerUFQkQ6i8hCEcnCGLszwDUYI+THf4AbnOHEy5yeVDiYCdwkIjWcsuoA12N6uX5cAazSc/PJqGo28DqQJiI1wywrkFMjEZkqIt9hjk87gzkKrDamdw1wCcYwzQ1KPi8or1TMgdWznJGCBGcI/Diw0pG/UFkwz6WwoWI/ZgLZ5O3NDgP+raqbAsK6ATtU9Qt/gKoq5nSl4PJTReR1EdmF0cMZzNFmbYPjAos0wLFKVTOBVVQcZzmLUoA1shYVFetU9UtV/Y+qzsecFSuYHmBIOEOd/mHf+zBOUl2AxUDlgKhPYjx5+2LOw9wnIjNEpF4RMr3m5DPQ+T0YM0ccOCRbBzP8G4wfHfmTiygjD0QkDnNsWh+MYe2J4eQfKvbzauR87gnKIivot98o/51zRsp/9QHqFiFSP2AvxiAXClX9AXMU2u1icDHQgfwOT41CyNHgVIIAAAO0SURBVJ1Pdscr2T+sPBq4HKOLfwBJRaUPCDu/KNktLPyw3sUWnoCqnhCRbUCnQqJdB9QCbtGApSUiUjUorzPAU8BTItIQY1wmAlUx87gFybBdRD7DzGPOcD6Xqer3AdH2Aw1DJG+IGcLcX4j8odAKuBgYpqozAzjdGBTPb9gbcM6JCsw8ZCD2OZ9jMAYwGKeLkKcf8C915lPDwCvAbMw8bD/M+Z1zguLsxhjfYATL3h1zYPflqvqpPzDYGa2Q9P6wXUVKbWHhwPZkLTwBx1C2wvSiCoLfmJ4JSNcGuLSgBKr6o6pOxxiccLxrXwPSRSQd86f/atD95UC3II/YeIzxXq2qRwrI95TzWSUoPBSnShjHn0D8G2PEBwWFB//ehDkY+2fOSEHwtaYA+XCGutMJb6jYj3cwh3v/EuNpvFBVg180VgHNRKRrQFmCWeoTiFC6SMYsJwqFGyRgswvnmXQjjF64hYUftidrUVGR5gzfCmY48beYodjnC0nzIWYe9lUR+X9OunHADgJeSEVkPsap5r8Y56ELMb3gqWHI9SbGmWgmcIL884Z/wThvfSAiYzFLj0YCbTCewQUhC9PLHCwia4BjmB7pN8B3wBMiko0xMA8GJ1bVTSIyG/izM8T8FWZo2d/jzXHiqYj8BpjvzEW/iXEQOg8zvL5DVScWIGNvTE83VA84JJwRiLnArzDPMtTa2JcxHt9vicgjmBepuzg3tO7vNX+O0eckR7fVgEcd+QM3KfHjBPC+iDyDGU4e56T/S7jyW1jE3PPKXvaK5kVo7+I9wMfAtUFxfQR5F2OWzGwETmKWkAzG/IlnBsR5CNN72of5I97k5FUpTBnnOnLNLuC+fw3pIUeOVcB1Ycjud/A6Q4BHMJCGWRZ0HLP5xmMYI5TH8xbT05uCGZI+ipnL7e3EuymorO4YL+sDjoyZwBtA90J4zwHmluCZXhbwHBMKiNMKWOQ8j73AXzGGV4FaAfF6AqudeN8C9xegS8XMW//R0dlJzPx7WqzruL3K1yWq5XLrVgsLi1KAiDyMmX9OUdUdEeSTiDF+I1U139pbNyAiCzHLf1qVIK0CT6jqo9GXzMJLsMPFFhYWAIhIH8y8cgZmiPVy4HfAm5EYWABVPU3oIdmoQERGYXrfWzAbYwzC9MLvdatMC4twYI2shYWFH0cwQ86jMfOVuzDzx2NjKVSYOIWZa27Guc077lLVv8dUKgvPww4XW1hYWFhYuAS7hMfCwsLCwsIlWCNrYWFhYWHhEv4/p4ik465DmL4AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(7,5))\n", "\n", "# vertical lines for photon steps\n", "for i in range(-10, 10):\n", " ax.axvline(1 - i * cct.freq[1], c='gray', ls='--', lw=0.5)\n", "\n", "# plot I-V data\n", "ax.plot(resp.voltage, resp.current, label='Unpumped')\n", "ax.plot(cct.vb, idc, 'r', label='Pumped')\n", "\n", "# label the steps\n", "for i, lbl in zip(range(3), ['1st', '2nd', '3rd']):\n", " vtmp = 1 - (i + 0.5) * cct.freq[1]\n", " itmp = np.interp(vtmp, cct.vb, idc)\n", " ax.annotate(\"{}\\nstep\".format(lbl), \n", " xy=(vtmp, itmp+0.15), \n", " xytext=(vtmp, itmp+0.15),\n", " va='center', ha='center', \n", " fontsize=16)\n", "# hw/e label\n", "ax.text(1-cct.freq[1]/2, 1.1, r'$\\hbar\\omega/e$', fontsize=16, ha='center', va='bottom')\n", "ax.annotate(\"\", xy=(1-cct.freq[1], 1), xytext=(1, 1), \n", " arrowprops=dict(arrowstyle=\"<->\", color='k'))\n", "\n", "# other labels\n", "ax.set(xlabel='Bias voltage / Vgap', xlim=[0,2])\n", "ax.set(ylabel='DC current / Igap', ylim=[0,2])\n", "ax.legend(frameon=True);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Plot the AC tunneling currents" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdwAAAE8CAYAAABq5wB3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdd3zN1//A8dfJTgRJiE3s2NSeLUUpSnWojm/R0n51fL/VraqJUqqDtlpditJ+f2ipUatWrdqzRewZRAgSsu99//44mTIEd2Sc5+Pxedzcz+dz733fhPu+53zOeR8lIhiGYRiGYV8uzg7AMAzDMIoCk3ANwzAMwwHcnB1AfqCUMv3qhmEYRiYiomz5fCbhpsjLtWyl1E3Py8s5zjovP8eW1/NMbCa2/Hieia1wxmZrpkvZMAzDMBzAJNxbEBISYpNz7HGeM17TvIfb44zYbBn/rTyfM/7P5FVReQ/m/8Ltn2drykwL0tdwC/rvIa9dKflZQX8PBT1+MO8hvzDvwflS4rdpv7Jp4RqGYRiGA5iEW0g4q4vElgr6eyjo8YN5D/mFeQ+Fk+lSpnB0KRuGYRi2Y7qUDcMwDKOAMgnXMAzDMBzAJFzDMAzDcACTcA3DMAzDAUzCTaGUQilFaGios0MxDMMwnCA0NDQtF9iDGaWMGaVsGIZhZGZGKRuGYRhGAWUSbhE3ePBglFK8+uqrWY6ldq8kJyc7JJY///wTpRR//vnnLT92/vz5TJgwwfZBGYZh2IhJuEVYXFwcv/zyCwA///yzwxKrPZiEaxhGfmcSbhH222+/ER0dTY8ePbhw4QLLli1zdkiGYRiFlkm4RdiPP/6Iv78/06dPx9vbmxkzZmR73oEDB+jUqRM+Pj6UL1+e9957D6vVmnb82rVrvPzyy1SpUgVPT0/Kli1Lly5dCAsLSzsnOjqal156iQoVKuDp6UlwcDATJ0686WoiVatWZeDAgVn2ZxxRPnDgQH788UfCw8PTRhhWrVo17dyLFy8ydOhQKlasiKenJ3Xq1OG7777L+y/KMAzDBtycHUBOlFKVgLeA5kBjwBuoJiIn8vBYL2A08BTgB+wG3hKRdXYLuIA5e/YsK1eu5LnnniMwMJAHH3yQefPmcfnyZfz9/TOd++CDD/LMM88wfPhwli9fzujRo3FxcUlLeMOGDWPhwoWMHTuWWrVqcenSJTZu3MiVK1cAsFqt9OzZk507d/L+++/TsGFDFi9ezKuvvkpkZCRjx469o/cycuRIIiMj2bZtGwsXLgTA09MT0Im+Xbt2xMXFERoaSrVq1Vi+fDlDhw4lISGBl19++Y5e2zAMI89EJF9uQEcgAlgCLAcEqJrHx/4MXAGGAJ2BeUAc0CSH8+V2pMSU7fbtt9+mnfftt9/mem5GTZs2zfG8IUOGpJ23ffv224o51YcffiiA/PXXXyIismzZMgHk66+/TjsnJCREABk3blymxw4ePFh8fX3l8uXLIiJSv359GTZsWI6vtWjRIgFk2rRpmfY/++yz4uHhIZGRkSIismbNGgFkzZo1aecEBQXJgAEDsjwnICEhIWn3BwwYIBUrVsxy3vvvvy+enp5y6NChLO+hVKlSkpSUlGPchmEUXSmfzTbNa/m5S3mdiJQVkR7AL3l9kFKqMfAEMExEvheRVUA/4BTwvn1CLXhmzJhBrVq1aNOmDQBdunShQoUK2XYr9+vXL9P9/v37c+3aNf755x8AWrRowfTp0xk7dizbt2/HYrFkOn/dunW4uLjw+OOPZ9r/1FNPkZiYyKZNm2z51jJZtmwZrVq1olq1aiQnJ6dt3bp149KlS+zfv99ur20YhpFRvu1SFhHrzc/KVm8gCZid4bmSlVKzgLeVUp4ikmCjGPN03nPPPcdzzz2Xp3N37NiRp/OaNWuWp/Oys23bNvbv389bb72V1u0L8NBDD/Hll19y6NAhateunba/bNmymR6fej88PByASZMmUa5cOaZOncqIESMICAjg6aef5oMPPsDHx4eoqCgCAgLSunlTlStXDoCoqKjbfi83c+HCBY4cOYK7u3u2xy9dumS31zYMw8go3ybcO1AfOC4isTfs3wd4ADVTfi6yfvzxRwDGjx/P+PHjsxyfMWMGY8aMSbsfERFB9erVM90HqFixIgC+vr6MGzeOcePGcfLkSX799VfefvttPDw8GD9+PAEBAURFRZGYmIiHh0fa85w/fx6AUqVK5Rirl5cXiYmJmfbdSoIuVaoUZcqU4fPPP8/2eHBwcJ6fyzAM407k5y7l2xUAXM5mf1SG40VWYmIis2bNolWrVqxZsybL1qRJE2bOnJmp9T5nzpxMzzFr1ix8fX1p0KBBlucPCgritddeo2HDhmldzvfccw9WqzVtzm+qn3/+GQ8PD1q3bp1jvEFBQWnPk+r333/Pcp6npydxcXFZ9nfv3p2wsDCqVKlC8+bNs2zFixfP8bWLmjNnzrBr1y5nh2EYhVZhbOEq9CCj7Pbn/KBsilWHhIQUusUMfv/9dy5dusSnn35Kx44dsxx//vnnGTp0aKZqT99//z1Wq5UWLVqwfPlypkyZQmhoKH5+fgC0adOG3r1707BhQ3x9fVm7di179uxhwIABANx///20b9+ef//730RGRlK/fn2WLFnClClTGD58OKVLl84x3v79+/PMM88wbNgwevXqxZ49e5g+fXqW8+rVq0dUVBRff/01zZs3x8vLi4YNGzJs2DBmz55Nhw4dGDZsGMHBwVy/fp2wsDDWr1/PggUL7uj3WViICF27diUsLIwVK1bQpUsXZ4dkGA4RGhrKqFGjHPNith6FZY8NGEweRymjr90ezGZ/v5TnqJ/NsbwOXCvwevfuLcWLF5fr169ne/zKlSvi7e0tAwYMSBul/Pfff0vHjh3Fy8tLypYtK++++65YLJa0x7z55pvSpEkTKVGihPj4+EiDBg3k888/z/S8V69elRdffFHKlSsn7u7uUqtWLZkwYYJYrda0c7IbpWyxWGTUqFFSpUoV8fb2lvvuu0+OHDmSZZTytWvXpH///uLn5yeABAUFpR2LioqSV155RapWrSru7u4SGBgo7du3l4kTJ97ZL7MQWbduXdpo+HPnzjk7HMNwOuwwSrlArBaklBoMfE8e5uEqpd4D3gX8JMN1XKVUKDAcKCE3DJoyqwUZRd3gwYP54YcfePvttxk3bpyzwzEMpzOrBeXNQsAdeDR1h1LKDXgM+OPGZGsYRV1sbGzadfrUywCQ91H4hmHkTb5OuEqpR5RSjwCpc2DuT9l3T8rxIKVUckqrFgAR2Y3uVv5MKTVYKdUZmAVUA0Ic/BYMI9+bP38+MTExtGzZkjp16jB79myaN2/ON9984+zQDKNQye+Dpm4seDE55XYtuhKVAlzJ+sVhEPABMAZd2nEP0F1EdtotUsMooH766ScgvXV77do1duzYQaVKlRg6dKgzQzOMQiVfJ9yb9Z+nXM/Nco6IxAGvpmyGYeSiffv2JCcn07dvXwC6du0KwOrVq0lKSsqxaIhhGLemQAyasjczaMowMqtTpw4HDx5kw4YNtGvXztnhGIbDmUFThmE4xH333QfAihUrnByJYRQeJuEaRhG2YcMGtmzZQkJC5sH7qd3Kf/zxhzPCMoxCyXQpY7qUjaKrQ4cObNiwgaVLl9K9e/e0/TExMQQEBCAiXLx4Ma2qmGEUFfboUs7Xg6YMw7Afi8XCzp164P6Nq08VL16cMWPGUK1atSyrPBmGcXtMl3IKpRRKqUJXO9kwchIWFkZsbCxBQUEEBgZmOf7WW2/Rr18/vL29nRCdYTheaGhoWi6wB5NwU6TWuiwqCXf69Olp/7CUUnh4eFCjRg3eeecd4uPj7fqaJ06csMvzG7dm+/btADRv3tzJkRhG/hAaGpqxxr7NmYRbxP3yyy9s2rSJxYsX061bN8aNG8cbb7zh7LAMB8hLwv3tt98YOHCg+ZJkGDZgEm4R16RJE1q3bk3Xrl2ZPHkyXbp04YcffsBqtTo7NMPOtm3bBuSecH/66Sd+/PFHM1rZMGzAJFwjk6ZNmxIXF8fFixcBXdj+rbfeolq1anh4eFCtWjU++OCDTAk5Pj6eYcOG0aBBA3x9fSlXrhwPPPAAYWFhznobxk1YLBZOnTqFUooWLVrkeF7qfFyTcA3jzplRykYmJ06coGTJkpQqVYrk5GS6devG/v37GTlyJA0bNmTz5s2MHj2aqKgoPv30UwASEhKIiYnh3XffpXz58kRFRTF58mRat25NWFgY5cqVc/K7Mm7k6upKeHg4J0+epGTJkjmelzofd9WqVSQmJuLh4eGoEA2j8LH1ArsFceM2FqCH/LHdrmnTpgkgYWFhkpSUJFFRUfLDDz+Iq6urTJo0SUREZsyYIYCsXbs202PHjBkj7u7uEhERke1zJycny/Xr18XX11cmTJiQ5TWPHz9++4EbDtegQQMBZP78+c4OxTAcBjssQG+6lIu4OnXq4O7uTkBAAM8++yzPP/88L730EgDLli0jKCiItm3bkpycnLbdd999JCUlsXnz5rTnmTNnDq1atcLPzw83NzeKFSvGtWvXOHjwoLPempELi8WS53MHDhwI6FHmhmHcPpNwb5Pz27Z6u1O//fYb27ZtY8mSJXTp0oXJkyczY8YMAC5cuMDJkydxd3fPtLVs2RKAS5cuAbBo0SIee+wx6taty//+9z+2bNnCtm3bCAwMtNsUI+PO1K5dm+bNmxMREXHTc5988klcXV35/fffiYyMdEB0hlE4mWu4RVyDBg2oWbMmAPfeey+NGjXijTfe4OGHH6ZUqVJUq1aNOXPmZPvYqlWrAjBr1ixq1qyZqQWUlJREVFSUvcM3bkN4eDjHjh2jRIkS2Ra8uFG5cuUYOnQoFSpUwM3NfGQYxu0y/3uMNJ6ennz88cf06dOHyZMn0717d+bOnYuvry916tTJ8XGxsbFZPohnzpx5S92WhuNs2bIFgJYtW+LikrdOrkmTJtkzJMMoEkzCNTLp3bs3LVq04JNPPuHIkSNMmzaNzp0789prr9G4cWMSExM5evQoCxcuZP78+fj4+NC9e3fmz5/PsGHD6NWrFzt27OCLL74wBe/zqdRr723atLmtx4uI3UrfGUZhZhKukcWYMWPo1q0bU6ZMYfny5Xz44Yd89913HD9+nGLFilGjRg169uyZNkVkyJAhnD59mqlTp/Ltt9/SokULFi1aRN++fZ38TozsrF27FoDWrVvf0uMSEhKYMGECCxcuZP369aZ72TBukVmeD708X+rPISEhRaaeslH0XLhwgXLlyuHh4cGlS5coVqxYnh+bmJhIvXr1OHr0KN999x1DhgyxY6SG4XihoaGMGjUq7b7YeHk+k3Ax6+EaRcfMmTN5+umn6datG8uWLbvlx8+ePZv+/fsTGBjIrl27qFixoh2iNAzns8d6uCbhYhKuUXRcuXKF5cuX4+/vn1a28VZYrVbuu+8+Vq1aRbt27VizZg3u7u52iNQwnMskXDsxCdcw8u7ChQvcddddnD17lhdffJFJkyaZQVRGoWOPhGsKXxiGcUvKlCnDnDlzcHNz46uvvkqrqW0YRu5MwjWMIuKTTz5h0KBB7Nq1646fq127dsybN48WLVrwzDPP2CA6wyj8TJcypkvZKPwsFgvVq1fn1KlTrFq1invvvdcmz2u1WtOKZ8TGxrJkyRIeeuihPBfUMIz8yh5dymYinWEUAStWrODUqVNUq1aNjh072ux5MybWt99+m0mTJlG/fn3eeOMNHnnkkVuadmRLFouFv/8+xuzZ59mxw5Vixdpy+TIkJUFc3AaqV1/BiBEPctdddzklPqNoMgnXMIqA77//HoBnn33Wbq3Pu+66iypVqrBv3z4GDhzIf/7zHx544AF69OhBp06dKF++vF1eF+D06dMsWLCQNWsusm1bAOHhDbBa2wK1sjm7Pbt2tWfhws9YsyaWdu3a2S0uw8jIdCljupSNwi0iIoJKlSohIpw6dYoKFSrY7bUSEhL46aefmDJlSqblG1999dW0wVX79+9nwYIFVK1aFT8/P4oXL46vry/FihVDRKhdu3ba4zZv3kxERARXr14lKiqKc+fOpW39+vXjoYeGsGIFTJlymlWrXICM84KtlCx5mDp1InjmmabUru1LcrLw5ZfnWbgwEBE3vL0fYe3at2jRooXdfidGwWSPLmWnLPie3zbuZCV3w8jnxo0bJ4D07t3boa976NAhmTBhgnTr1k0WLlyYtv/bb78VINvNxcUl03M0aNDghnNcBdoIhErp0kdFqcwLVpYocU26dz8r330XLZGROcc2fnxyymPOSrNm99nrV2AUYNhhAXrTwqVotnCnT5/OoEGDOHz4cNryfPmJUsqU2bSB5ORkqlatSnh4OEuXLqV79+7ODonNmzczd+5cTp48SXR0NDExMcTExBAbG4uLiwsHDx5Mm9f74osvcvBgLDExbbl0qRnh4XWJj/dOey4PD2jfHrp3h27doGFDyMuUYIsF7rnHwsaNrsA3hIV1Ijg42E7v2CiIzKApO0r9D24+5POHTZs2UalSJWeHUeC5ubmxYcMGZs6cSbdu3ZwdDqAXTcht4YS4OFi3DpYvhz///Ir9+zMfr1VLJ9fu3aFjR7idcVmurjB5siuNGwM8zI8/fsHYsaNv/YmMQuXGWsq2Zlq4mBZufmzhOkNCQgKenp7ODqPIEYGwMJ1gly2DtWshPj79uK8vdO6sk2y3blC9uu1et1KlBM6e9WTlyut07uycEdVG/mQqTRl207FjR9q3b8+yZcto0qQJ3t7e3HXXXWzZsoXk5GTeeecdypcvT0BAAAMHDuT69euZHh8SEkLTpk0pWbIkpUuX5t577800aCbVzp076dChA15eXlSuXJmxY8cSEhKSpTSgUipTT0NoaChKKQ4fPkzPnj3x9fUlKCiI999/H6vVmnZefHw8w4YNo0GDBvj6+lKuXDkeeOABwsLCMj3/9OnTUUqxbt06Hn30Ufz8/GjVqhWffPIJnp6eREZGZjpfRKhevTqPP/747f6KHW79+vVYLBZnh5GtK1dg7lx47jkICoJ69WDYMJ104+Phrrtg+HD480+4dAnmz4ehQ22XbEF3Pfftq79g/fmnSbaG/ZkuZSPNkSNHeOONNxgxYgS+vr68+eab9O7dm969e5OcnMz06dM5cOAAb7zxBmXKlOGjjz5Ke2x4eDjDhg2jUqVKXL9+nZ9++om7776b7du306hRIwAuXrxI586dqVChAjNmzMDDw4OJEydy4sSJPMfYt29fBg0axLBhw1i0aBEhISFUrlyZQYMGAbqVGhMTw7vvvkv58uWJiopi8uTJtG7dmrCwMMqVK5fp+Z588kkef/xxfv31V5KTk2nZsiUjR45k2rRpvPnmm2nn/fHHHxw/fpypU6fewW/YcbZu3UrHjh1p06ZNvlhgwGKBHTt0Ql2+HDZv1vtSBQbCfffpbuKuXaFsWcfE1asXfPUVLF4Mo02PsmFvth6FVRA3bmeUcsahkc7cbtO0adMEkMOHD4uIyD333CNubm5y9OjRtHMWLFgggHTu3DnTY/v27StVq1bN8bmTk5MlKSlJateuLf/5z3/S9g8fPlzc3d3l9OnTaftiY2OlTJkycuPfAJCQkJC0+yEhIQLI1KlTM53XoEED6dq1a66xXL9+XXx9fWXChAlZ3v8rr7yS5TEDBgyQGjVqiNVqzfSeg4ODc3yd/CQ+Pl7q168vgLz++utOi+PsWZHp00X69xcpVSrzP1s3N5G77xb54AOR7dtFLBbnxBgXJ+LjYxUQ6dLlaecEYeRL2GGUsulSNtLUrl2b6hn67OrUqQOQZbBNnTp1OHPmTOqXFQBWrlxJp06dKFWqFG5ubri7u3Po0CEOHjyYds7mzZtp06ZNpsFQ3t7e9OzZM88x3nhugwYNOHXqVKZ9c+bMoVWrVvj5+eHm5kaxYsW4du1aplhS9e3bN8u+F154gaNHj7Jq1SoAzp07x6JFi3j++efzHKczjRgxgn379lGrVi3ef/99h71uQgKsXg1vvQWNG0OFCjBwIMyapbuFq1aFf/8bfvtN31+7Ft55B5o1A2dVgvTygk6d9CWJlSvduXr1qnMCMYoE06V8uzIkm8LC398/030PD48c9ycnJ2OxWHBzc2Pnzp306NGDbt268cMPP1C+fHlcXV0ZPHgw8RlGv5w7d44GDRpked2yt9B/GBAQkOm+p6dnptdYtGgRjz32GAMGDCAkJITSpUvj4uJCjx49Mp2XKrvqRy1btqR58+Z88803dOnShSlTpuDm5saAAQPyHKezzJo1i08//RQ3NzemTZuGt7f3zR90B44c0QOdli+HNWsg46V9Hx89ijh1yk6tWnmbsuNovXq5sngxQBe2b99O586dnR2SUUiZhGvcsblz5+Lm5sa8efMyXSu8fPkyfn5+affLly/PhQsXsjw+IiLCZrHMmjWLmjVrMn369LR9SUlJREVFZXt+Tuu4Dh06lOeff57w8HCmTJnCo48+miXZ5ze7du1KW7ln4sSJdilZGBOjE2tqkj12LPPxhg3TRxO3b69bkPldq1apPzVm69b5JuEadpPnhKuUChSRyJufaRQ1sbGxuLq6Zkpeq1evTiuWn6p169Z88sknnDlzJq1bOS4ujsW6eWGzWNzcMv+znjlz5i2P1n388cd5/fXXeeKJJzh16hT//ve/bRajvUybNo24uDgGDBjAiy++aJPntFphz570KTt//aUXAEjl768HO3Xrpm8rVsz5ufKrunXBxcWK1VqbTZt2OzscoxDLNeEqpe4B3gdaAh5KqURgC/CeiKxzQHxGAdC9e3c+++wzBg4cyKBBgzh06BCjR4+m4g2fvq+++ipff/013bp1IyQkBE9PTyZMmICnp2eOLc3biWX+/PkMGzaMXr16sWPHDr744otMLe288Pb2ZuDAgUycOJGGDRvStm1bm8RnT5999hm1a9dmyJAhd/T7jIyEP/5IH1GcsVPCxQXatEnvJm7eXBeRKMi8vKBGjSQOH/Zk06ZoZ4djFGI5DlVQSj0KrAbKAB8D/wE+AcoCq5VSjzgkQiPf69atG1988QUbN26kV69eTJ06lRkzZmQpqFG6dGlWrVqFv78/Tz/9NC+88AJdunShb9++lCxZ0iaxDBkyhBEjRjB79mweeOABFi9ezKJFi27r+R999FGAfD1Y6p9//uHSpUuAXirvpZdeuuXiHUlJsH49jBihE2jZsvDUUzBzpk62lSrBs8/CnDlw8aJu5b73nu6KLejJNlWzZnq8wsWLFQgPD3dyNEZhlWOlKaXUAeAw8KCIWDPsdwEWAjVEpK5DorSzolhpKr+wWCw0bdo0LRnnJyNGjODzzz/n7NmzlChRwtnhZLFy5UoeeeQRgoODWbVqFb6+vnl+7IkT6d3Eq1dDdIaGnacn3HNPevnEunXz52AnWxo/Ht5+G1q0+IsFC6rZdSlBo2BwdC3lasCrGZMtgIhYlVKTgbm2DMTZTC1lxxg5ciQ1a9YkKCiIS5cuMWXKFPbu3cuSJUucHVqaXbt2cfDgQT7//HOee+65fJdsRYQJEybw5ptvYrVaqVSpUpbr1jeKjdVVm1KT7KFDmY/XqZPeTXz33XqEcVGSUpsFH5+2mFxbdNm7lnJuxSD+Bp7O4dgA4G9bTwp21oZZns9hRo4cKdWrVxcvLy/x9vaWVq1ayfz5850dViZBQUHi5eUlffr0kejoaGeHk8mFCxfkgQceSFuu7t133xVLNlUjrFaRv/8W+fhjkS5dRDw9MxeeKFFC5KGHRL79VuTECSe8kXwmPFz/Xvz99e/OMHDk8nxKqd7A50A/EdmWYX8rYBbwsoj8br+vAo5jupSNgmDevHm88MILRERE4Ofnx9SpUzMV7oiKgpUrdQv2jz8g46VIpfT12dRu4lat4CaN4iJFBAIDhUuXFOPH/x9vvllwamYb9mGPLuXcEu56oCZ60NRpIAI9YKpyys9HMpwuInKPLQNzJJNwjYLgyy+/5OWXX+aee+5hxowZVKxYha1b07uJt23T03hSlSuXPie2a1coXdp5sRcEnTsLq1croCfnzv2Qpe62UbQ4+hquBQhL2VIdT9kMw7CzqKgo/vnnH+6++25AF+NwdQ3Cza0nr7/uwsqVcPly+vnu7pkHOzVqVPgHO9lSo0aK1asBGrNnzx6TcA2byzHhikhHB8ZhGEaKixcvMmHCBCZNmoSIJ99/H8b27aVZvtyVffseyHRuzZqZF2O/hYHKxg30YvQAjdizZ0+WGuKGcafMVRzDyCeOHj3K5Mlf8/XXq4mL6wDMwcWlE088kV4fsVixzIux16jhvHgLm9SRytCY3bsXOjMUo5DK8Rpu2glK+QO1gCxVUcWO1aaUUpWBiUBXQAErgVdE5FSuD9SPzelN3SUiWWq3mWu4hjNdvJhI165j2b27HNAdqJrpeJMm6a3Ytm0hZU0Jw8bi48HXV7BYrNSp04IDB3Y6OyTDiRw9aMoLmAr0Qye8LETELnVmlFI+wB4gAXgXPQViDOADNBKR67k8PDXhTge+veHQXhGJze58k3ANR0lMTGLq1O1cvtyK5ctd2LgRkpPTj/v5JdGjh3tafWJzKdFx6tWzcuCAC0q15Pr1tXZfbcnIvxw9aGok0BE953Ym8CIQDwwEygP/tWUgNxgCVAeCReQIgFJqL7ry1fPAhDw8R7iIbLZfiIaRNyLC9u3H+e6746xa5caJE8GItEk77uICTZvG0b079O3rTdOm7k5bH7aoa9LEhQMHwMurFadOnSI4ONjZIRmFSG4J92H0wgWz0Al3i4jsBKYppX5B930ttVNcvYHNqckWQESOK6U2An3IW8I1DJu4du0abm5ueN1krTkRIT4+Hm9vbywW2LQpmeeem8eRIzVJSmqC/g6pubmdp3PnJAYPrkznzuDvb1pS+UGjRvB//weDB39BcLAZ4m3YVm7fo6sA+0TEAiQBxTIcmwo8Zse46gP/ZLN/H1Avj88xVCmVoJSKVUqtVkp1sF14RlGwf/9++vfvT4kSJZg/f37a/ilTplCuXDnq1q1Lo0aNaNiwIRUrVsTTsyrlyr1N/7JxWFoAACAASURBVP4QGAgdOrhx4EA/kpKaAkmUKbOXBx/cyMKFx0lMLMeyZZV55BG9xJ2RP6SOVN671yRbw/Zya+FeAlInGZwGGgPrU+6XBuz5lTwAuJzN/iggLx9PPwG/A2eBIOAN9ApHXUXkT1sFaRReixcvpnfv3lit1ixL3Z0/f56IiAgiIiLQ3w0fRHe8tCApCWbP1udVrw4NG56ma1cr//pXZUqUaISRv6WOVN67F65du46vb7HcH2AYtyKnmo/AL8A7KT+PRyfA4ejkdR5YbOs6kxleOxEYl83+D4Dk23i+4sBJYEMOxyW7LSQkJPsim0ahFhMTI5UqVRJAnn76aTl16lSm44cPx8urr16VypXjM9Un9va2Ss+eIpMmiRw+7KTgjTtitYoEBFgFRFxcgiQiIsLZIRl2FhISku3nP3aopZxbkmoOPCTpCWsuumvZCvwFVLF1MBleOwL4Npv9k4HI23zOyUBCDsdu/69lFDqvvfaaANK8eXNJTk5O2x8RITJ4sIiLS3qSLV1a5JlnRBYuFImNdWLQhs106pT69+0hv/zyi7PDMZzEHgk3t0pT24HtKT/HAA8rpTwBTxGJzulxNrIP3Vd3o3rA/tt8ToX+1mIYOTpw4ACfffYZLi4ufPPNN7imrLC+eDE88YReN9bNDR59FAYPhk6dCs8i7IbWqBGsWQPQhLVr1/LII484OySjkLilSlMikoCeG2tvC4FPlFLVReQYgFKqKtAOePtWn0wpVQLoCWyxYYxGIXTx4kXatGlD3bp1adasGaAXB3joIUhM1MUnPvsMzGyRwqtVq9Sf7mHt2tedGYpRyORW+OLpXB5nBa4Cu0TkjM2DUqoYuvBFHOmFL0aju7Ybici1lPOCgKPA+yLyfsq+14FgYA3pg6ZS93UWkfXcwBS+MG4kIiil2LED2rfXVYheegm++MIsCFDYRUZC2bKCbl8EcPHiaUqVKuXssAwHc3Thi+mkd8FmfNGM+6xKqdnAIBFJtFVQInJdKXUvurTjzJTXWoUu7Xgtw6kKcCXz9KaDQN+UrSQQDWwEnhWRrbaK0SjcUkcmDx+uk+3AgfD55ybZFgWBgdC0qWLHDi+gAxs2bKBPnz7ODssoBHJLuO2An4FFwK+kr4fbD+gFvAA0AEahRwC/Y8vARNdMfvgm55zghrKTIrIoJWbDuCUbNmzAx8eHu+66C6UUW7bAihVQvDh8+imm+lMR0q0b7NgB0I0lS5aYhGvYRG5dynOBgyKSJZEqpcYCdUWkr1JqNPCkiFTP8iQFhOlSNgDatm3Lpk2bmD9/Pn369KF3b1i0CN5+G8aNc3Z0hiOtW6fXFi5f/jJHjnji4+Pj7JAMB3P04gXRQF8RWZXNsS7APBEpoZTqCvwuIp62DMyRTMI1zp8/T4UKFfDw8ODixYscO+ZL48bg7Q0nTkCZMs6O0HCkpCQICIBr1+DUKahc2dkRGY5mj4SbWydZItAsh2PNUo6nPkeuq/cYRn73xx9/ICJ07twZX19ffvlF7x840CTbosjdHe69V/+8ZImup52YaLNhKkYRlVvC/QUYpZR6TSkVpJTyTrl9HQgFUgrY0QQ9UMkwCqy1a9cC0LlzZyB1Hibcf7+zIjKc7eGUESTDh0dToUINvvrqK+cGZBR4uXUpewPfA49nc/h/wBARiVdK9QRixI6L0dub6VI2atasydGjR9m2bRt16jTH3x+sVoiKgpIlnR2d4QwWCzRtqusqwwhcXD5kwYIF9OrVy9mhGQ7g0C5lEYkTkaeAuug1cIen3NYVkX+JSHzKeYsLcrJNpZRCKUVoaKizQzEcLDw8nKNHj1K8eHGaNGmStiB8s2Ym2RZlrq66yAmAu/t7WK2VeOyxx5g6dSoWi8W5wRl2ERoampYL7CHHFm5RYlq4Rdv69evp27cvLVq0YOnSpbz9NowfD2++qW+Noq1vX5g/Hzw9Y0hIeAL4neDgYJ588klGjhwJkFqT/ZY+qEWExMRE4uPjSUhIICEhAYtFKFu2CvHxev73unVbiY5O5Pr1ZK5dS+b6dQuxsVZiYy3UqFGf8uWrER8Phw+f5K+/dhIfr0hIUCQmupCYqEhKciUpyZV27TqTnOxKfDzs2bOXq1fP4OYWScWKEbRuncxrr/WgSZMm9vj1FVh2H6WslLqlmYYiYrVlMM5iEq5htVq5cuUKAQEBtGoFW7fC0qW6lKNRtEVG6jraK1fq+25up0lOXkW1ah68+OITAFy+fIWxYz/A3d0bd/cSuLgUQ8QTq9UTEU9atbqHYsUCiY+HgwdPEh5+CREv9CqnGW+9nPMmWUbjxr8zb96rVK9eYGd42pQjEq6VvBf4FxG5pVrM+ZVJuEaq6Gg9HUQpuHwZfH1v/hij8LNaYcIEPR87Ksq+r+XpqaejeXpCdHQESiXg5paMm5sFNzcL7u4WPDysVKkSSJUq5fDygvj4K5w7dwwvL4W3t8Lb2wUfHxeKFXPBx8eNevWq4+vrhpcXXL16lZgYxdGjFjZuTGDDBn+SkvSszjZt/mTduo64FYpP9jvjiIQbyi2sqCMio2wZjLOYhFt0xcfHk5ycjG9KZl26FHr0gDZt4K+/nBycke9YLLB7t+4BuXABrl7V+5UCESsWSxKurom4uyfj5QVeXoKXFwQG+lKypAdeXuDqmoinp1CihDs+Pi54e+sE6+Wlk6yjK5pdugRjx8YxcaIHIq506AD/939QsaJj48hvHFr4oigxCbfomjNnDk8++SRDhw7liy++YOxYGDECXnkFJk50dnSG4Tjr1kH//nDuHJQubaVy5RDGjGnL/fd3t9sgovzM0YUvDKPQW7FiBcnJyVRM+Tr/zz96f8OGTgzKMJzg7rt1671rV7h40YVdu0bTs6cvtWq9w/DhE1i7di3h4eEkJyc7O9QCy7RwMS3cokpEqFq1KqdOnWLHjh00bdqURo3g779hyxZo2dLZERqG41ks8PHHyYwenURsrHeGI2eBs5QsCW3bNk/p/k5i5crlaVNpdEs4/edGjRpRrlw5RODkyZMcPXo029d0d3ejQ4e70+5v3ryJuLj4lHuZP5uDgqpQs2ZNAC5fjmLXrt0old3nt9CmTdu0Otj//PM3Fy5EAKCUC0opXFz05ufnT+PGjfD0BA8PC5s3/8n27V1Ml7I9mIRbNB06dIjg4GBKly5NxObN8M679Jv7GHMtDxITYwZMGUVbTAxMmpTArFnR7N9fEovFw9khOZhj18M1jEJtxYoVADzVsiUuHTvCmTP8zDySKmzE17e5U2MzDGcrXhzeeceTd94JJClJX9s9dw6io4XEREVCAsTFWQkPP4vVakXEitVqxWqVlFsrFSpUoHjx4gBcvBhJREREtq/l4uJCvXr10u4fPBhGYmIS2bWDSpUqTfny5QGIiYnh2LHjWc5JfVxwcDCenp5pLeyrV68goge4JSfrGC0WCyVKFKdOnfokJEB0dAJr125lyZI7/AVmw7RwMS3courBBx9k84IFHCtRAp/oaBKK+eN5/TIR3kGUPb0DSpVydoiGYTiJ3QdNKaW+VUp1V0q52/JFDCO/sVgsrFmzhv8CPtHR0L49nw49yhZaUjbuJLz1lrNDNAyjkLlxlHJVYD5wUSk1WynVXylVwvFhOZ6ppVy0uLq6cnDnTl5LvVD78cdsP+rPAH7U92fN0ouhGoZRZDi8lrJSqjjQC+gDdEfXGlsL/AYsFJGzdonEiUyXchE1eTK8+CK0bg2bNlG7Nhw+DNeatKfY7o0wfToMGODsKA3DcAKHzMMVkRgR+T8R6Q8EAn2B48BI4LRSaotSarhSqq4tAzEMR0lKSsKSlJS+FMywYcTFwZEjeoUYj+cG6v3TpzsrRMMwCqFbGjSllGqDbvn2AYKBQyJSx06xOYxp4RYt//vf//jj5ZeZHhUFVarA0aPs3OtGs2ZQty7s33QVypXTy7UcPw5Vqzo7ZMMwHMzplaZEZJOIvC0idYF6wDRbBmMYjjBz5kyeTq1A/5//gJtbWoWpBg3Qi+A+9JDeMWOGU2I0DKPwue3SjiISJiJmtVCjQDl9+jQRy5dzL2AtVgwGDwbInHAh/drtjz+S7WRAwzCMW2RqKRtFyvvvv89/UxKoy7PP6tYs2STczp31cinHjsGGDU6I1DCMwsYkXKPIOHjwICumTuVxQJSC//437ViWhOvqCk8/rX82g6cMw7ABk3CNImPkyJF8YLXiAah+/aB6dUCvaXr6tF6LtEaNDA9I7VaeMweuX3d4vIZhFC4m4RpFQnJyMmWPH+dJQDw94cMP047t26dv69XTDds0wcF6ju61azBvnkPjNQyj8Mkx4SqljimlGudwrIFS6pj9wjIM23JTii9Sqseo117LNNUnS3dyRoMG6duvv7ZvgIZhFHq5tXCrAp45HPMCgmwejWHY2JUrV7h69SqMGoXatg3Kl4fhwzOdk2vCffJJ8PODTZtg61b7B2w41759ev61YdjBzbqUc5oP0Ry4YuNYDMOmLly4QKdOnfigXTsYPRpcXGDmzCwL3eaacIsVgyFD9M+ff27fgA3n+vJL/Y+gXj1YuNBMBzNs7sbVgoYppU4ppU6hk+2i1PsZtkjgK2CZMwK2F7N4QeGydetW2rZti/vu3bx74IDeOXq0nu5zg1wTLuh6yy4uevDU2UJXStwA/Xd95x398/Hj0KcPTJrk3JgMh3Po4gVKqT7Agyl3BwBLgMgbHpMA7AemiEisXaJyMFPasfCIjY3lo48+YsyYMbSxWFjq4oKv1QqPPqpXAHLJ3Klz4QKULasX2756FXL8f/bIIzB3LrzxBnz0kf3fiOFYTzwB//d/8MAD0L69Xp4xIEAn3xJFYsE04wb2KO2IiGS7ocs2VsvpeGHa9K/BKOjmzZsnFStWFEBeBElydRUBkccfF0lMzPYxy5bpU9q0ucmTb9umT/TyEgkPt33whvP89Vf63/b4cRGrVeTuu/W+9993dnSGk6TkBZvmmhyv4YrIIBE5btPsbhg2lpycnPbz5cuXcQ8PZ22JEnwJuFks8NJL+rqtu3u2j1+5Ut926HCTF2reXNdXjo/XXdNG4TEtpST8f/6jR68rBe+/r/d9+ilcvuy00IzCJdfVgpRS1YF+QBX0yOSMRESetWNsDmO6lAsOi8XCli1bWLJkCUuXLqVly5Z8/fXXcOECyR99hPriC1yTknQ34JQpuis5F3fdBbt3w4oV0KXLTV78wAF9odfFRf9cs6bt3pjhHElJemWoqCjYuxcaNkw/1rkzrF4NI0emJ2CjyLBHl3KOCTfleu4v6IFVF9DXbjMSEaluy2CcxSTcdNeuXSMxMRF/f3+7DRy4VTt37mTJkiVs3bqVDRs2cDlDi+P+smVZ3KcPasaM9OkcTz4J48frWsi5SL1+6+WlGzFeN36lzM4zz+gWUc+esGhRLhd9jQJh2TK4/369LuO+fZn/nhs36uu5xYvra7mlSjkvTsPh7JFw3XI5Ngb4E3hSRG4cOGUUIjvmzWPjmDHE/vMPJ5KSCAOu1K/P7tThu8A///xDjRo18Pb2tulrJyUlcf78ec6ePUt4eDinT5/m4MGD/Otf/6JNmzYALF26lJEjRwL621/vSpV4vmJFOkREUPzECfjuO/1kDzwAI0ZAq1Z5eu1Vq/Rthw55TLYAH3ygB08tXqyrTz38cN7frJH/zJmjbx97LOuXp3btoFs3WL4cPvkExo1zfHxGoZJbC/c60FdE/nBsSI5XJFu4Iuz/9FOSRo+mcXR0lsNX3N3xe+YZ6NeP2BYtKO7nh9VqJTAwkMqVK1O5cmXKlCmDr68vzz77LPXr1wdg/fr1bNq0CavVmrZdv36d6OhofH19GT9+fMrLC5UrV+bs2bNk97v/+OOPef311yE2ln2zZxM2YwYtr1+nwqFDuF69mn6iv79eZGDIEEiJIa+efRamTtWN4TffvIUHfv01vPACVKigu5bNKNaCKTFRd3FcuaL/jnXqZD1nyxZd3rNYMb1yVJkyjo/TcApHdynvAD4Rkf+z5QvmR0Uq4SYlwa+/6lrCe/cCEA2cqVWLyvffj29MDPz1F+rgwbSHWPz9+SM5meXXr/OP1coBIONs1CVLlnD//fcD8O677/LBBx9k+9KVKlXi9OnTme6fPXuWKmXK0CAwkEbFi1PX05Pa7u4EK4XfyZNw+DBYrZmfqFo1fX3t4Yfh3nvBw+OWfw0iEBSkFy3YuVNfy80zqxXattUfxs88Az/8cMuvb+QDv/+ue0UaNYI9e3I+r1cv3aMxZEh6b4pR6Dk64XYGPgP6iEihrptc2BOuJTGRfdOnc+m772i6fz8l4+IAkAoV2NymDfU//5wSGa93isDff+vuttmz4ciRLM+Z5O1NtL8/sV5elKpVC5/KlaFYMY6Fh3P8zBnE1RUlgrvVipeLC94uLhRzc6N6xYp6gEpkJJaICFyiolC5rcTj6qqvr911F3TsCJ066YR7h3buhGbNIDAQzp/PMj335vbt0yOX4+N1F/NDD91xTIaD/etf8NNP+jJBatGL7OzfD02a6C+rf/4J99zjsBAN53F0wl0P1ABKAYeBqBtOEREpFP/yClvCjT55kmMLFnBl9Wo8t22j1rlzlM7w/iy1a+P6+uu6K9Yzp3LZKUTg6FE9f2bHDggL091vly7ZLmB3d535qlbVS+ZVq6bXyWvYUJfZy/MF1rwRgfvu02/pv/+Fzz67zSeaNElPJQkI0C2kSpVsGqdhR/Hxuns4JkZ/ocy0LmM2QkL0SOXatfXf2sb/Jo38x9EJ909yrqUMgIh0smUwzlLQEq5YrVw+dozzW7ZwZfduEsPCqOHqSuWkJN0yzdBtm+q0mxvH69Uj8JVXqDNgAOqWm3Q3iIzU5fCionTyvXQJ4uL0dbHERN0acHXVCd3DI/NtQACULq23wEA9CtSBo31TexL9/XWP9W0PPhWBHj30SNdWrWDt2pt/gTHyh99+070SzZrB9u03Pz8hQbdyw8Jg4EB98d+MUC/UHDpKWUQ62vKF8rvUKTAhISEOradsTU4m5vx5rp89S+z588SeP09iZCQJZ89iPX+eDsHBev5KZCSHNmzAMzqaUsnJBAABOT2nlxcHXF25XLEitGlDpccfJ6hrVyrfaZLNKDBQbwXMxYvw2mv65/feu8OZHkrBjBm6a3nLFl1z+fvvzQdxQTB7tr597LG8ne/pqbufO3SA6dN1i/jdd+0WnuEcoaGhjBo1ym7Pn2vhi6JCKSXratZM/6BUClEKlKJR48b4+fuDUhw9fpwz4eFpxzJu3sWK0apVK93qSUhgw4oVkJCAa2IirklJuCUl4Z6cjJvFQgV/f0q6ukJ0NBIdjbqNv0EMcM7LiyslSxJXvjylmjalQa9eeqRujRo3rKRuWCy6ITpkCJw7p9eW37v3tsZbZbVzp55CEh8PkyfD0KE2eFLDbq5f193JsbFw4oQePZdXCxZA3776//n48bq2tvmCVSg5tEs55QUrAq8Bd6Ov5T4gIv8opV4BNonIFlsG4yxK3U7Ks51rwDUXF667uRHv7k68lxeJPj4kBgTQundvPCtVgsBATsXHI4GB+AcHU7xixTvvFi5grFbdax0bm7ft+nWIiNA1CzZu1L3foGsZ/PwzVKliw+B+/hmeegrc3GDNGv0iRv40c6Yev9C6tV7n+FZNnqx7M0DPLZs0CWw8P91wPkdfw60PrAcswCagJ9BCRHYqpSYCZUXkCVsG4yxKKVk3cKD+1pq6oIHVCiI0atAAv5IlQYSjR45w5vTptGNp51ut+Hh707J5c/1t18uLzbt3Y/HwwNXXF9dixXDz9cW9eHHcihenTJUqBFSqBCVL6uuXbrnVH8nfUpPgjVtsbO73czont+1O1wWvWVN/Pr7+up1+5a++ChMn6rmdW7faOKMbNtOmDWzerEt/Pnub1WnnztWjnOPioFYtfSmhMIxeTv1cK2Jf5rPj6IS7DCgOdAPigUSgeUrCfRQYb0o7Fm4xMXoWxM6deiDn+fO6RsDly/r22jU9lsSRvL3Bx0fXIfDxuflWqpQe/Ny4sf5ctKvkZF2ZaPVqPbp6wwY9MsvIP3bv1lPMSpaE8HD9D+l27dihW8r79+v7PXro67qtWxesbubISD0IbPly/Z4SEnQRkHbt9BeSpk2dHaFTODrhXgMeF5FFSilXIIn0hHs3sExEfGwZjLOYhJtZWJgeUDR/vh5sfDPe3pk3H5+s+/Ky/2bJ08urAHzxvnxZD6zZtw/uvlt/iJkpJPnH88/r4hV3NB8sg4QEfS33o4/0NQzQo5kHDdJD4W0wZ9wuRHQvzFdf6QFkiYk5n9u0qR788Pjj+otKEeHohBuNrqOcXcJ9CPheRApFNW+TcLWkJL3u9uef665iFxdo2VL3lNWqpaeZ+vuDn5++9fXVuaQgfZl3iNOndSvn7Fno108vbJ7vvykUAVev6nKcsbE5l3K8XRcv6qX8vvsufbAA6Nfo0UNXRmvaVK9M5CxJSbBrFyxdqqvNpdZKV0rHOGiQrqDm46OP/fKLHoWfuliIh4deUuv++3W3fKNGOS57WRg4OuGuBKJF5KFsEu4swEdEetsymBtevzIwEegKKGAl8IqInMrDY72A0cBTgB+wG3hLRNblcH6RT7jnz+uV7DZs0AOcBw/Wq5LdZMEdIyd79+qWbnQ0DBsGEyY4OyJj9GjdddO5c/pCyLYWH6+7hubOhT/+0H//jMqW1S3g+vX1tY7UrXx5/S32TmcXWK064Z86pUdgnzyprwft3auTbcaqbgEBust46NCcW+Lx8XqRjh9+0NeXMpZZ9faGFi3Sr9ekbpUqFYr56I5OuPegk9wa4H/AD8BwoD7QH7jbXqOUlVI+wB70koDvogtwjAF8gEYikkstQFBK/Ywe5PUGcAx4EbgfaCMiu7M5v0gn3E2bdFnic+d0A2DuXN1AM+7Q6tXQvbtuWUyYoBOv4RxXrujEdvWq48ozJiXp/1xLl+rb3bv16+dEKZ0ES5XS3Ueenumbl5e+tVjSi8ukbvHx6YMrLl/W3cU5qV1bl0jt2/fW65BfuAALF8L69fr9HD6c87mlSukPk/Ll028DAvSXitQt431f33zXVeaMaUE90fWUM9Y9OwG8KCJLbRnIDa/7X2ACECwiR1L2VUOXmHxTRHJsLiilGqNbtM+IyLSUfW7APuBgdq3yoppwo6P1ZawxY/Rnw9136/LJZcs6O7JC5H//0+vzgr5W1q+fc+MpqlJLM9qzdXszIrrVuWuXTlYnTujt+HE9f+3KFdu8TsmSeoR8UJDeqlXT3b+NG9t2taPISF3wJSxMv59Dh3Rr+tw5/cXgVri46KRbrFj6bXY/e3tn/iKSWr3uZvfd3HTvgYtL9rfZ7FOlSjkm4aZ0ITcAzopIpFKqJlAGuCQiB7M8wMaUUqsALxFpd8P+tQC51XBWSo0ERgJ+IhKbYf8o4G2ghIgk3PAYhyXcDLOJ8rTdMAPpjs5PStL/r48e1V3Hv/2WfrnplVf0uI9CfEnGeT76SF8c9/DQg6g6dnR2REVLZKQuBhMTo//ht2t388c4Q3JyeqnU1HlwCQnpW3y8ThweHumbu7u+TW0x+vk5f5qhxaKvaZ89q5Pv2bP6mlVUVHor/MYttwVMnESBwxKuC7o7t6cz1sNVSp0HFojI8zfsnww8KiI51hRMub58l4gE37C/HzAbaCAi+244Jg0aCFYred5Sk9vNzrlxy286dIBRo/QiPIadiOhRsZMm6dbH+vV6YQbDMZ55BqZN0937S+3WMWfciaQknXRTt2vXsv85NlZ/+UhMzPxl5Gb7kpP1h7LFkn6b8edsjqkrVxxTS1lErEqp08AdTFK7IwHA5Wz2RwE3m9iY22NTj2eROmDPUbKpDomLS/b7b/WcnM5zc9M9ShUr6lr7HTvqgZP57NJJ4aOULohx7pweHdq9u74GZgpj2N/GjTrZenjo4fdG/uTurlvnfn4Of2l710/OKLe+h2+BV5RSi0Ukl0ladpNdezAvqUHdzmP37tVJKi9bakK72Tk5JT+jCHJ11SUFL1yAdet00t2wQXcDGvaRlJRe1/qtt/SAIcO4QWhoaLYL1ig7fFjnlnCLowdLHUupOnWOzIlMRCTE5hFpl8m+JepP9q3XjKKA7JoO/hmOZ2F6+Ay78/LSxe87dNBdKg88oAfwmDq89jFqlF6uslo1GD7c2dEYRq7TgqzZHkgnImKXJWmUUqsBDxFpf8P+P9Ex5zZo6j30VKIbB02Foqc1OXXQlGFw5owuHHDmDPTpo7uZnT3QpbDZsEFP/RHR6xR36ODsiIwCxh7TgnIsfyMiLjfZ7Ln+20KgtVIqrVazUqoq0C7l2M0e6w48muGxbsBjwB83JlvDcLhKlfRagX5+usX70kv5c0RdQXX5sl5YwGqFt982ydbIN7JNuEopD6XUbyk1k53he/R83wVKqT5Kqd7AAuA0+toyAEqpIKVUckqrFoCUwhazgc+UUoOVUp2BWUA1wF5d4IZxa+rXh0WL9BzBb7/Vk6GNO2e16mUST5yAZs0gm2tzhuEs2SbclEFSXXI6bm8plaTuBQ4BM4GfgePAvSJyLcOpCnAla5yDgGno6lSLgcpAdxHZaefQDSPv2rdPr7P83nu6fJ5xZ0aPhiVL9GC0X3+9tUpKhmFnN1ue708R+dCxITmeuYZrONXXX8MLL+iRzMuX62pIxq379VddEFwpPd+2WzdnR2QUYA69hgu8BjyrlHpJKVVJKeWqlHLJuNkyEMMosoYOhTff1BPuH31Ul8gzbs3Wrfq6LcCHH5pka+RLeRmlnFPTT0SkUAytNC1cw+ksFnjoIV0cvnZt2LzZLF6fV4cP6+75Cxf0MlfffWcmvBt3zNGrBYWSc7IFQEQcU57DzkzCNfKFmBidOPbu1d3KS5ea4tY3c+aMro186hR07QqLF5vfvlfqQAAAIABJREFUmWETDl8tqKgwCdfIN06ehJYtdWvtxRfhyy+dHVH+dfGinvITFqbXk1yxQq8oYxg24OhruIZhOFpQkF7A3MMDvvoKfv7Z2RHlT9HRcP/9Otk2aKBbtibZGvlcbl3K72V7IJ2IyGjbh+R4Sqm0X0JISEi2dTUNw6G+/Rb+/W/w8dEDgurXd3ZE+ce1a9Crl64gVb26ripVvryzozIKgRsXMnDkNdzcSjtKSjD2rDblMKZL2ch3RODpp+GnnyA4GLZtg+LFnR2V88XEQM+eeonDChX0bfXqN3+cYdwip5d2BEoBA4F/gJq2DMQwjAyUgm++0S3bgwf16Nui/qUwOlqvsrR+vS6PmdrCNYwC4rYGTSmlhgFdRaSH7UNyPNPCNfKtgweheXPdjTppkq67XBRdvarn1m7ZotcRXrPGJFvDrvLToKk9gLPqLBtG0REcDFOn6p9ffx0OHHBuPM5w+bKe8rNlC1Stalq2RoF1uwm3FxBpy0AMw8jBo4/CoEGQkKCv6yYlOTsix4mKgi5d9DXs6tXhzz910jWMAii3QVNTs9ntATQAGgIhIlIoljgxXcpGvhcdDQ0b6gIPo0bpxQ4Ku4sXdct2926oWRNWr4bKlZ0dlVFEOLrS1AmyVpqKB06il7v7sbBkKZNwjQJhzRq49169WP3mzXr5ucIqMlJX2/r7b13qcvVqqFjR2VEZRYipNGUnJuEaBcZ//wtffAH16sGOHeDl5eyIbC8iQn+x2L8f6tTRydbMszUcLD8NmjIMwxnGjdMtvv37C+fi6ufOQceO+v3Vq6ev2ZpkaxQSuXUpvwVUEpGXszn2BXBaRD62c3wOYVq4RoGyZQu0aaPXz921S5c2LAzOnoVOnfTyhA0bwsqVUKaMs6MyiihHt3AHAXtzOLY75bhhGI7WqpUu+5icrNfSteZWFK6AOHNGt2wPHYJGjXQ3skm2RiGTW8KtAhzO4dgxIMj24TiPUgqllKmjbBQMY8fqhLRhA8yY4exo7sypU3DPPXpd2yZNdLItXdrZURlFUGhoaFousIfcupQvAv8VkSzLlSilngK+EJEAu0TlYKZL2SiQZs7U83LLldPJqiCulnPypO5GPn5cj7r+4w8IKBQfK0YB5+gu5fXAG0opzxuC8AReSzluGIazPPmkXjv3/Hn4uAAOpzh3Tk/9OX4cWrTQ69maZGsUYrm1cBsDfwEXgZ+AcKAi8BR6EYN2IrLHQXHalWnhGgXWxo3Qvj14e+vrn5UqOTuivLl0SXcj79sHTZvCqlXg5+fsqAwjjaNXC9oDdEIXungL+DLl9jjQsbAkW8Mo0Nq1g4cfhri4gjNNKHXVn3379NSf5ctNsjWKhDwVvlBKeQP+wGURibN7VA5mWrhGgXb4sC4QoZRu5ebnwv6xsXD//bBunY4zdV1bw8hnnFb4QkTiRORsYUy2hlHg1aoFTz0FFguMycflza1W+Ne/dLKtWFHPszXJ1ihCTGlHTAvXKASOHNGtXNBr6Nao4dx4svPmm3pwl58f/PUX1K3r7IgMI0emtKNhGNmrWVO3Hi0WGD/e2dFk9f33Otm6ucHcuSbZGkWSaeFiWrhGIXHwoG7lennpYhKBgc6OSFuxQl+3tVjghx/+v707j7Kjqtc+/n2SkGAYA8ogUxhDAF+CIDJKGO6NkAAJswIBkdervoqKXhmXRgXRJQ7c60W4igETZmSeNSjKECaZCUGQEOYpBDKQQbLfP351TKVyuvt0n7H7PJ+1zqruOruq9u46u39nV+3aG447rtk5MuuSW7hm1rFhw2D0aFiwAM47r9m5CU8/DYccEsH25JMdbK2tLdfClbQKsDiltKDsBtKKwAoppTkNyF9DuIVrfcYdd8RgEmuvHaM4DRrU9Tb1Mm9eDMzx1FMRdC+/HPr5O771DnVv4UoaCbwD7NbJNrsB70japZYZaTaPpWx9wp57wrbbxpyyl13WvHykBP/xHxFshw+HiRMdbK3lNXQsZUmXAiumlMZ1upF0FbAwpXRkXXLVYG7hWp/y29/C5z8fI1D9tUkjsE6cGJePBw+GBx6IAS7MepFG3MPdDbiigu2uAnavZUbMrEYOOwxWWilmEnrmmcYf/x//gBNOiJ/PPdfB1ixTDLhrAy9WsN1LWVozazUrrwyHHx4/X3hhY4/9wQcxg9HcuXHfdvz4xh7frIUVA+57QCXTdawBzK19dsysJj73uVhedFFMVN8ov/xlTKiw7rrRU7pO98LMeqNiwH0EGFPBdvtnac2sFe26awz5+MorMcdsI7z0Epx+evx83nmw5pqNOa5ZL1EMuBOBYyWN7WgDSQcBxwAX1DNjZlYFCY45Jn6+8srGHPOEE+JS8rhxcMABjTmmWS9S7jncK4GDgBuBG4jp+QA2Ilq2Y4CrUkqHNzCfdeVeytYnTZsWHZbWWCMeExowoH7Huv12GDUq7h9Pm9Z75uU160A9eimXC7j9gFOBE4HVgVICAbOBnwJnpZSW1DIjzeSAa33WllvGkI9TpsBee9XnGEuWxCTyjz4a4zh/+9v1OY5ZAzVkaMeU0pKU0hlEL+Rdgc9kr12BtVNKZ/alYGvWp43LHqm/5pr6HePiiyPYbrDB0seBzGw5PZq8QNIewDEppT4xMKpbuNZn3X8/fPKTMf/szJm1H+1pwYIYw3nmzHgEqXTf2KyXa+rkBZI2k/R9Sc8DfwIOq2VGzKwOdtghgu3LL8ODD9Z+/5MmRbDdZhs46qja79+sD+k04EpaTdIXJN0FTAdOI8Za/hLw0Qbkz8yq0a/f0h7Dt95a230vWQJnnx0/n3oq9O9f2/2b9THLBVxJ/STtJ+ky4FXgPGAo8D9Zkq+nlM5PKb3XuGzWnycvsD5r771j+ac/1Xa/118fQ0dutBEcemht923WBI2evOBs4EhgLWABcC1wEfBHYFVgFjAypfSXuuSmSXwP1/q0t96KyegHDYLZs2OC+lrYZRe491445xx3lrI+pxH3cE8kgu3NwIYppSNTSrdnvZIdkcx6ow9/OKbsW7gwAmQtPPRQ7GvIEE8qb1ahYsD9LTAHGA1Ml/RLSTs2PltmVlN77hnLWl1WnjgxluPHx2AXZtalZQJuSul4YB3gKOAh4IvAvZKmASfhVq5Z71TLgLtwIVxySfx87LHV78+sTXT6HK6kdYHxwNFAaVLLqcC5xPCOC+qewwbwPVzr82bPjskE+vePnwcP7vm+rroqOkmNGAEPP1y7PJq1kIY/h5tSejWl9OOU0jbAJ4lAuznwO6IHs5n1BquvDtttB4sXwz33VLev0uVkt27NuqXigS9SSg+klL5CPH97CHBnvTKVPZp0iqQZkhZIelTSwRVue6GkVOb1i3rl16xX2G23WN5/f8/38eab8TzvgAHw2c/WJl9mbaLb04eklBYDV2evevkB8C1ioI2HgCOAKyWNSSndXMH2bwLF+cHcIrf2tv32sXzooZ7v45ZbYsCLffaJR43MrGJ1nK+rZyStRQTbH6WUsmFs+JOkzYAfEY8sdWVRSmlqvfJo1it9/OOxrCbg3nhjLPffv/r8mLWZGo9kXhOjgIHA5ML6ycDHJG3c+CyZ9QFbbhmdpV54Ad5+u/vbL1oEt90WP48eXdu8mbWBVgy4WwMLgWcL65/MllvRtbUkvSXpn5KekXSSJA/0au2tf//oWQw9a+XedRe89x5svTVs7O+9Zt3VigF3DWB2med0ZuXe78wjwDeJ2YwOIDp3nQWc39lGpfEz8y+Pq2x9TjX3cW+6KZZu3Vofkh8/Of+qh7rfw5W0D/CHCpLemVIaCYjyA2xU9BdIKRV7I98saS7wdUk/Tin9vYPtKtm9We9WzX3c0v3bMWNqlx+zJpswYULZxlU9gm4jOk3dAwyvIN38bDkLGKLlR6MYknu/uy4Fvg7sAJQNuGZtoact3BdfjJmBVlsNdt659vkyawN1D7gppfnA093Y5ElgELApy97HLd27faoH2Sh9VXEz1trb8OHwoQ/BjBnRcWrNNSvb7u67Y7nLLvEMrpl1Wyvew70VWERME5h3FPBESun5Huzzs0SwfaDKvJn1bgMGxMxB0L1hGUsBd9dda58nszbRcl9VU0pvSPo5cIqkOcDfgMOBvYAD82klTQE2Siltlv2+ETAJuIxoHQ8CxgHHAuenlJ5rVDnMWtbHPgZTp8JTT8UAFpVwwDWrWssF3MxpwFzga8TsRdOBw1JKNxTS9WfZMswh7vGeBKxNtGqnAScQ40Cb2fCsS8W0aZWlnzMHHn00Wsc7erZOs55qyYCbUvoAOCN7dZZuZOH3WcDY+uXMrA/obsC9774YznH77aubZciszbXiPVwzq6fuBlxfTjarCQdcs3azwQbRUn3jDZhVwVN2DrhmNeGAa9Zu+vWLcZUBnu7iib0lS6KDFcQjQWbWYw64Zu2o0svKzz4bnabWWw8++tH658usD3PAzXj8ZGsrpRZuVwH3kUdiud129c2PWQvIj6tcDy3ZS7kZPJaytZVKW7ilwTFKswyZ9WH5cZXrEXTdwjVrR5UGXLdwzWpGbtnB8vMkmPVxixZFT+UlS2DevBhfuZx11oHXX4fnnoNNNmlsHs2aSBIppZo2c93CNWtHAwfCZptBSjB9evk0r70WwXbVVT3hvFkNOOCatauuLivn79/WqROJWTtxwDVrV10FXN+/NaspB1yzdtWdFq6ZVc0B16xduYVr1lDupYx7KVubmjsXVlklOlDNmxfT75XMmROdpVZYIdINHNi8fJo1gXspm1ntrLxyTGSwaBE8//yy7z32WCy32cbB1qxGHHDN2llHl5V9/9as5hxwzdpZRwHX92/Nas4BN+PJC6wtuYVr9i/1nrzAnaZwpylrY3/5C+yxB+y4I9x3X6xbvDju7y5aBO++G52nzNqMO02ZWW3lW7ilL53TpkWw3XRTB1uzGnLANWtnH/kIrLlmPAb0yiuxzvdvzerCAdes3ZUC6513xtL3b83qwgHXrN2NGRPLa6+NpVu4ZnXhTlO405S1uRkzYvq9lVeGJ56AYcNg4UJ49dWYD9esDbnTlJnV3tChsO22MYTjwQdHsB03zsHWrMYccM0Mxo6N5UMPxfK005qXF7M+ygHXzODAA5f+vO++sP32zcuLWR/le7j4Hq4ZKcHmm8Nzz8Fdd8GuuzY7R2ZNVY97uA64OOCaAfDMM/Dyy7Dnns3OiVnTudNUHXksZWt7W2zhYGttzWMpN4BbuGZmlucWrpmZWS/lgGtmZtYADrhmZmYN4IBrZmbWAA64ZmZmDeCAa2Zm1gAOuGZmZg3ggGtmZtYADrhmZmYN4IBrZmbWAA64ZmZmDeCAm/HkBWZm7c2TFzSAJy8wM7M8T15gZmbWSzngmpmZNYADrpmZWQM44JqZmTVASwZcSSdKukHSq5KSpAnd3H6spIclLZD0gqTTJfWvU3ZbQl/oXd3by9Db8w8uQ6twGfqmluylLGka8B7wN+CLwPdSShMq3HYUcDNwAXApsB3wQ+CclNJJHWzT63spZz3qmp2NqvT2MvT2/IPL0CpchuarRy/lVg24/VJKSyQNABbTvYD7MPBeSmmP3LrvAKcDG6aUXiuzjQNuC+jtZejt+QeXoVW4DM3XNo8FpZSW9GQ7SRsAI4DJhbcmASsA+1aTr0oukVR6GaXW6ZpxTJehZ5qRt1pf3mtG3lyGnu3PdaHn6WqtJVu4Jd1t4Ur6NHALsEtK6d7Ce/OAc1NK/1lmu4pauJV8Y6v0W10z0rVy3ipN57w5b62Yznnrs3nr+y3cKqyRLd8p8947uffNzMwaakC9DyBpH+APFSS9M6U0strDZctyX106/aZS6diZlaSr5b5qna6V81ZpOuetZ+mct/qmc956lq6V81ZrdQ+4wD3A8ArSza/BsWZly3It2dVz7y+j1pcNzMzMiuoecFNK84Gn632czJPZcmvgX/dwJQ0FBgNPNSgfZmZmy+hT93BTSjOBR4EjC28dRXS+uqXhmTIzM6NFA66kHSQdAhyUrdpK0iHZa3Au3RRJz0raQNJVkt4FtgBGSrpY0khJ3yCewT2ng2dwV5T0E8WoVu9LulfSp8qk6yfpFEkzFCNYPSrp4BqW+V9lkPSepKslbVjBdjtI+l9JT0uaL2lmVvaNy6SdoRi5q/ga28wyZNuWy1eSNKKQrlXPw4ROyrCgkLZu50HS+pL+O/scz8/2O7TCbVulLvSoDC1WF6o5D61SF3p6HlqlLhwi6feK0QbflzRd0lmSVqlg2/rUhZRSy72AC4mOT+VeQ3Pp/gy8APwdeAIYCxyYrVuYvWYC3wH6d3Csi4HZwP8F9gauBt4HRhTSnZnt71vAnsD5wBJgvxqUd3CZMjwOPAes1MW2ZwN3A18G9gA+C0wD3gY2KKSdAdwK7FR4DWlmGbLtEzCxTN4G95LzsH6ZvO9NXFm5ooHnYSTwOjHa2m3FOtPFtk2vC9WUoVXqQg3OQ9PrQpXnoVXqwlTgCuKK5x7A17PP91SgXzPqQtUnpdkv4GvAB8BmuXUbA/8ETuxi222zD9HncusGANOB63Pr1sr+qN8rbD8FeKzJZfhImXUbZSf9+2U+3JNb7TxkaRNwRhdpWvY8dLC/o7NyjW7geeiX+/n4bvyTbIm6UGUZWqIuVFOGLH3T60K1ZSizr2bUhXKfh/FZPvbqZLu61YWWvKTcTQcAU1NKz5ZWpJSeJ77pHljBtouBy3Pb/hO4DBglaVC2ehQwkOVHsJoMfKzcJatu6nEZUkpvlln3AvAmsF6V+eqOas5DpVr2PHTgGKKFcFuV+apY6uEobbROXehxGVqoLlRzHirVsuehA82oC8t9HoAHsmVnn4e61YW+EHC3Ji4BFj0JbFXBts+n6Eld3HYgsFku3ULg2TLpqOA4XammDMuRNJz49jWtzNv7Z/djFkqaWqt7VtSmDF/K8jVf0h2Sdi9zjF5xHiStT1xiujirrEX1Og891Sp1oaaaVBdqodl1oWZarC7skS3LfR5K6lYX+kLAXYPyI0vNAoZUsW3p/dJydsquF3SSrqeqKcMyFMNhnkd8q7+g8PYNwFeJb2ZHAguAayQd1d0Ml1FtGSYT9972Ab4ArAncIWlk4Ri94jwQl9D6AReVea+e56GnWqUu1EwT60K1WqEu1FJL1AVJ6wHfB/6YUnqwk6R1qwuNGPiiEbo9slQuTSXbVpquGrXa/y+BXYh7Jct8aFJKX11m59I1RAeCs1j+skhP9LgMKaWjc7/+VdJ1RGvzDGC33L56y3kYDzycUnpsuQPU/zz0RCvVhVppZl3osRaqC7XS9LogaWXgOqI/xue6Sk6d6kJfaOF2NEbyEMp/S8mb1cm2pfdLyyHScmOBFdP1VDVl+BdJZxHfiI9LKd3eVfqU0gfAlcD6ktat9DgdqEkZcnmbA9wEfCK3urechx2BLSn/jX45NT4PPdUqdaEmmlwXaqpJdaEmWqEuSFoRuB7YBBiVUnqpi03qVhf6QsB9kriWXrQVXY8s9SSwsXLP9ua2XcTSa/NPAoOATcuko4LjdKWaMgAg6TTgZOBrKaVJ3Th2Z+NPd0fVZSij+A2y5c9D5hjim/Ql3dimVuehp1qlLlStBepCPTS6LtRKU+uCpBWA3wM7Eo/qPF7BZvWrC/Xojt3IF/Fs1T+BTXLrhhK9zL7ZxbYjiJN6TG7dAOKG+g25daXu398tbP9H4PFmliFLe0JWjlO7edwBwIPAC80uQ5n9rUo8Q31nbzkPWfqBxHOf1zXjPBT2251HalqiLlRThix90+tCtWUos33D60ItytDsukA0KK8g7gvv3Y3t6lYXanpSmvECViK+cTxOPLpxADG84z+AlXPpNsr+mX6nsP1lxOXC44kHnK/KTtDHC+l+lK0/kXgg/FfE8337N7MMwBFZPm5h+YfHt8ql+0xW1vFEj8EjgL9mH6wjmlyGbwG/JgYqGEl8K36c+Da5e284D7n3Dsr+pgd1cIy6nofsGIdkr19l+/1S9vserV4XqilDq9SFKsvQEnWh2s9SK9SFXJ7PKPN5WL8ZdaEmJ6XZL2BD4rLBe8Ac4FoK38SIlkoCJhTWfwj4GfBa9oe7DxhZ5hj9iSEiS6NYPQYc0uwy0PmoXH/OpdsJuIN4Fm4x8C7xTWxUC5Rhf+JZ17eyvL1N3HPZsbech9x712X5H9jB/htxHjr9PLR6XehpGVqpLlRRhpapC9V8llqhLhCDanSU/wnNqAvKNjIzM7M66gudpszMzFqeA66ZmVkDOOCamZk1gAOumZlZAzjgmpmZNYADrpmZWQM44FqfJelYSSn3+kDSy5KukDSskHaCpF75jFwx75JWz9Z9vJn5qpSknbPzs36Z98Zl7x3WyfaTJM2VtFJ9c2pWHQdcaweHAjsDnwJOAbYDpkhaLZfmN1mavmB14LtArwi4wFjgwVR+UPkbicETji7zXmkWmHHA71NK8+qXRbPq9ZXp+cw680hKqTTg+N2SXgH+QEzddgtA9s++q1lErD4OBMpOMpBSWizpUuCLkj6SUnqzkORgYkjOimajMWsmt3CtHb2XLVcorSh3SVnSVyTdK2mWpNmSpkoaXUgzQNIPJD0naYGktyTdJWk3OiDpXEmvZxOk59cPkvSOpF/k1g2TdE12/PezPHy6k30PBZ7Pfv117nL6sdn7/y7pZkmvSpov6QlJ35TUv7CfwZJ+JeltSXOyPOyS31cu7R6SpmTp5km6TdI2HeWxsO2WwDBiCM2OXEQ0Dj5T5r3xwIvAn3P73CQr43xJb0j6qaQvZHkfmkt3hKQ7JL2ZXZJ+WNIxZfKYJJ0p6TRJL2Xn4S+SRlRSRrMSB1xrB/2zwDhI0nDgh8Ab5P5Jd2Aocan5UOBwYhaTGyXtm0tzEvAN4L+AUcTk1lMoP59mye+ImUb+vbB+DHE5eBKApI8CdwHbAl8BDgNmAzcV8pD3KjFoPMQk3jtnr5uydZtk+TsOGE0EswnAmYX9/G+W5uxsf9OBi4sHy76ATAHmAkcRg+6vQkycvkGHf4GlxgHPppSe7ChBSulBYiq0ZS4rZ/d8RwKTUkpLsnUDiasX2wJfBo4FNgZOK7PrTYhB6Y8kLmvfAPxG0hfLpB0P7Eech2OBtYnbEp2dZ7Nl1XKga7/8aqUX8Y+x3MDlLwOfKKSdENWhw331I1pZt5Obboy4x3h1D/L2DHBpYd21wFO5388mZjLZLLeuPxH8/tZR3lk6IPvxXeRBWZlOI2ZG6ZetH0bMePLtQvr/yvZ7bG7ds8CUQrpVicH3f1HB3+E+4CcVpPt2duzhuXUnZ+uG5dZ9IVu3Y6Gcj9LJ9HK58/tr4NHCeykrz0qFv/Fi4AfN/pz71XtebuFaOxgHfIKYhHosMTH0zVlrt0OStpd0o6TXicC3GPg3IiCVPADsl11y3C1rYVViMnCgpFWyY60B7Eu0fks+BUxNS+8/k1L6ALgUGCFp1QqPlS/TupLOl/QCMeXbYmL6stWJVjfAJ4kgdWVh86sK+9qcmHz74uwKwoDsMvl84N4s/53mhTgvnV1OLpkMfMCyrdyjgftSStNz63YCZqaU7i+tSCklYgao4vE3l3SppJeJv8NiYjq2YcW0wM0p1ykrpTQDmErf6WhnDeCAa+3giZTSgymlB1JK1xHz3IpoGZaVXQ4tXRr+KtHB6hPArcCKuaQ/JHoEH0DM5fm2pImSPtxFniZl+zkk+/0I4p5y/rLtGsQl4qLXsvwP6eIYy5DUj5jqbQwRZPciylS6nFwq17rZ8o3CLl4v/F4K0BewNGCVXmOANbvI0ljgTSI4dyql9AoxddtRCjsAW7F8Z6l1y+R7ubxnvZtLl55PBnYn/ha/BQZ1tX1u3Xpd5d2sxL2Ure2klN6X9A/g/3SS7NPAasBhKfe4iqTBhX0tBn4M/FjSOkSg+RkwmLjv21Eenpd0N3Hfc2K2/HNK6cVcslnAOmU2X4e4zDmrk/yXsymwA3B0Smlyrkz7F9KVgvxaLO2ABXHfMu/tbHkKEQyLFnWRn7HA9Sm7/1qBi4BLiPu2Y4n5Ry8vpHmVCMRFxbzvTEw+vntK6a7SymJHtk62L617uctcm2XcwrW2kwXNTYnWVUdKgXVxbrstgF072iCl9FpK6TdE8Kmkl+4kYKSkkUQA+F3h/TuBnQo9a/sTgfzhlNKcDva7MFt+qLC+XJlWIDoN5d1HBPRDC+uLv08nJvneOruCUHw91kH+yC6Hj6Syy8kl1xKTlB9H9Fi+MaVU/NIxFdhQ0o65Y4l4fCiv3N9iCPGIUjn7KTewRnZOdqKC1rlZiVu41g5GZJd4RVxy/Apxufa/O9nmj8R9299J+mm23feAmeS+qEq6juiQ8zei49F2ROv4/ArydQXREWky8D7L32f8OdHx6w+Svks8zvRlYAuih3FHXidan0dIegyYR7RUpwEvAGdK+oAINt8obpxSmi7pEuAH2WXoh4jLz6WW8JIsXZL0/4DrsnvXVxCdi9YmLsHPTCn9rIM8jiZawOVaxmVlVyauBD5PnMtyz95eSPQcv1rSacSXquNZevm91Jq+h/h7/k/2t10JOD3Lf35AlJL3gdsl/YS45Py9bPufV5p/s6b32vLLr3q9KN9L+Q3gDmBUIe0ECr2UicdwngYWEI+lHEH8Q5+RS/NNolX1NvFPeXq2rxUqzOOVWb4u6eD90jOq72b5mAp8uoK8lzqHLSbXsxgYQTxqNJ8Y6OP7REBapgcv0QL8FXHZei5x73d0lu7AwrF2Jnprv5PlcQZwGbBzJ+W+HLiyB+d0t9x5HNBBmk2Bm7Pz8SZwDhGEE7BaLt1ewMNZuueAEzr4WybiPvep2d9sAXG/fkSzP+N+9a6XUuqVw8eaWYNJ+k/ifvXQlNLMKvYzkAiEX04pLfdsbz1IupF4pGjTHmybgDNTSqfXPmfWTnxJ2cyWI2kMcR/6EeIy7O7At4CZdGKIAAAAZElEQVQrqgm2ACmlRZS/bFsTkk4kWuV/JwbhOJRonX+pXsc0q4QDrpmVM4e4LH0ycX/zZeJ+83ebmakKLSTuTW/I0oFCjk8pXdDUXFnb8yVlMzOzBvBjQWZmZg3ggGtmZtYA/x+y1HxtB5vn4AAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(7,5))\n", "\n", "ax.plot(cct.vb, np.abs(iac), 'k--', label=r'Absolute')\n", "ax.plot(cct.vb, np.real(iac), 'b', label=r\"Real\")\n", "ax.plot(cct.vb, np.imag(iac), 'r', label=r\"Imaginary\")\n", "ax.set(xlabel='Bias voltage / Vgap', xlim=[0,2])\n", "ax.set(ylabel='AC current / Igap', ylim=[-1,1.2])\n", "ax.legend(frameon=False);" ] } ], "metadata": { "kernelspec": { "display_name": "Python [conda env:qmix] *", "language": "python", "name": "conda-env-qmix-py" }, "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.7" } }, "nbformat": 4, "nbformat_minor": 2 }