{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Two color MOT\n", "\n", "This example covers calculating the forces in a two-color type-II\n", "three-dimensional MOT. This example is based on Fig. 1 of M.R. Tarbutt and T.C. Steimle, “Modeling magneto-optical trapping of CaF molecules” *Physical Review A* **92**, 053401 (2015). http://dx.doi.org/10.1103/PhysRevA.92.053401" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import scipy.constants as cts\n", "import pylcp\n", "import pylcp.tools" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Specify the problem\n", "Here, we use $F=2 \\rightarrow F'=1$ with $g_l=1$ and $g_u=0$. For details about the units, chosen here see `01_F0_to_F1_1D_MOT_capture.ipynb`. This notebook uses the *hybrid* unit system that requires us to neglect the magnetic forces in the rate equation." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "Hg, Bgq = pylcp.hamiltonians.singleF(F=2, gF=0.5, muB=1)\n", "He, Beq = pylcp.hamiltonians.singleF(F=1, gF=0, muB=1)\n", "\n", "dijq = pylcp.hamiltonians.dqij_two_bare_hyperfine(2, 1)\n", "\n", "# Define the full Hamiltonian:\n", "hamiltonian = pylcp.hamiltonian(Hg, He, Bgq, Beq, dijq)\n", "\n", "# Define the magnetic field:\n", "magField = pylcp.quadrupoleMagneticField(1)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Run the detuning of the blue detuned beam\n", "\n", "At every detuning point, calculate the trapping frequency and damping coefficient for the MOT. " ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Define the detunings:\n", "dets = np.linspace(-5, 5, 101)\n", "s = 3.6\n", "\n", "# The red detuned beams are constant in this calculation, so let's make that\n", "# collections once:\n", "r_beams = pylcp.conventional3DMOTBeams(\n", " delta=-1, s=s, pol=+1,\n", " beam_type=pylcp.infinitePlaneWaveBeam\n", ")\n", "\n", "it = np.nditer([dets, None, None])\n", "for (det_i, omega_i, beta_i) in it:\n", " # Make the blue-detued beams:\n", " b_beams = pylcp.conventional3DMOTBeams(\n", " delta=det_i, s=s, pol=-1,\n", " beam_type=pylcp.infinitePlaneWaveBeam\n", " )\n", " \n", " all_beams = pylcp.laserBeams(b_beams.beam_vector + r_beams.beam_vector)\n", "\n", " trap = pylcp.rateeq(all_beams, magField, hamiltonian, include_mag_forces=False)\n", " omega_i[...] = trap.trapping_frequencies(axes=[2], eps=0.0001)\n", " beta_i[...] = trap.damping_coeff(axes=[2], eps=0.0001)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Plot up the result:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA74AAAGqCAYAAAAsppilAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOzdd5xcdb3/8ddne3bTK6QXICSUUAKEzhUQMPQiClhQQdQr+sN+vV7x2rDrRbEg0hT1amiCcgUkhBBIgEASkpCQ3ssm2U12N9s/vz/O7MyZZbPZMjNnduf9fDzmMXPOnDnnM0p25jPfz/fzNXdHREREREREpLfKizoAERERERERkXRS4isiIiIiIiK9mhJfERERERER6dWU+IqIiIiIiEivpsRXREREREREejUlviIiIiIiItKrKfEVERERERGRXk2Jr4iIiIiIiPRqSnxFRERERESkV1PiKyIiIiIiIr2aEl8RERERERHp1ZT4ioiIiIiISK+mxFdERERERER6tYKoA5CAmW0DSoGNUcciIiIpMwaocfdDog4kl+gzVUSkV+rWZ6q5e4rjka4ws73FxcX9Jk2aFHUoIiKSIqtXr6aurm6fu/ePOpZcos9UEZHep7ufqTkz4mtmlwBfBSYDrwKfdfelbRx3EXAbcDRQArwUO3Zl6JgzgTltXKbE3eu6GOLGSZMmTV269B0hiYhID3XUUUexbNkyjTpmnj5TRUR6me5+pubEHF8zmwk8DDwIzADKgTlmNrKNw8cCjwNnAucDo4B/mFlR6JgJwDJgYPjWjaRXRERERERE0iRXRnxvB2a5+y8AzOxGYCtwc+y5OHf/dXjbzL4J/AWYCrwR2z0ReNvdK9MatYiIiIiIiHRbrx/xNbMJwHTgkZZ97l4LzAZmduAULTXke0L7JgCrUhSiiIiIiIiIpFEujPiOi92vabV/A0HZ8wGZ2WnAN4B73X196KmJwJVm9gFgPXB/y2jywZjZgSYcqQOHiIiIiIhIGvT6EV+gpd11Rav9VcCgtl5gZl8ys/3Ai8BS4NZWh3wCOBG4iqDJ1c/M7M6URSwiIiIiIiIpkwsjvuWx+9ZtrwcA2w7wml8RlEZPJegE/ZqZne7u5QDu/mbsuJXAXDNz4HNmdoe7b24vGHc/qq39sZHgqQd7MyIiIiIiItI5uTDiuyl2P77V/nGh55K4e6W7r3T3Rwk6O08A/r2da8wBDBjRvVBFREREREQk1XIh8V0Ru13WssPMyoBzgFkdeP1+oA4oaueY6bFjlnU5ShEREREREUmLXp/4ursD3wGuN7ObzGwycC9QCdxrZu82s11mdiGAmX3RzM4ys5Fmdjzwx9ipHow9f6SZ3WJmU8xsrJl9Gvgy8B+xbtEiIiIiIiKSRXJhji/u/oCZFQKfBb4LzANOc/cKM2smGK1tMLNi4DDgJmA0wfzgBcAZ7r48drqBwHXAN4FiguZX17v7XzP5nkRERERERKRjciLxBXD3e4B72tj/DDAytOvmg5znZeCs1EYnIiIiIiIi6dLrS51FJLfU1Dcye8UOVu+sijoUERGRA9pX28Bzb+1gxz7NlBPJhJwZ8RWR3qu+sZlnlm/nycVbefat7dQ2NJNncNv5R/DJcw4jL8+iDlFERASAjbtruPfFdfzvqxupqmtkaN9inrntLAaWttdHVUS6S4mviPRotQ1NXP/b+by2fk/S/maHH/5zJfPX7ubH7z2OYf2KI4pQREQEmpudb/xtKQ++vJ5mT+wvr6pj1sLNfPSMCdEFJ5IDVOosIj3ab+aseUfSG/bC2+W8539eYNHGigxGJSIikmzWwk3c/1Jy0tviz69sIFiIRETSRYmviPRYWyr2c9fsVfHtMw8fyj0fms6r/3keFx19SHz/zn11fPIPC6mua4wiTBERyXHuzj1z18a3xw4u5bbzj4hvr9xexev6gVYkrZT4ikiP9b2n3qK2oRmAoX2L+MX1J3DulBEM7VvMXdefwDcvO4qC2PzezRX7+cH/rYgyXBERyVHzVu/irW374tu/+/B0bj33cGZMHBzf9+cFG6MITSRnKPEVkR7p1XW7eeyNLfHtL1wwmf4lhfFtM+MDp47nE+dMiu+7/6V17ZZFi4iIpEN4tPfsI4Zx2PB+ALzvpLHx/X9bvIUqVSaJpI0SXxHpcZqbndv/tjS+ffSo/lx94pg2j/33dx3GpGFlALjDl2Ytpq6xKSNxioiIrN5Zxb/e2hHfDjexuvDoQ+hXEvSaralv4olFW97xehFJDSW+ItLjPLZoM29u3hvfvv2So8g/wJJFxQX5fO+qY7HY06t2VPGLf61q81gREZFUu/fFxGjvESP6cubhQ+PbJYX5XHH8qPj2n15RubNIuijxFZEe59HXE7+IXzJtJNPHD27naJg+fjAfnDEuvn3X7NVs3F2TtvhEREQAKmrqmfXa5vj2R06fgFnyD7XXnpSoWHpjYwUrQnOBRSR1lPiKSI+yv76Jl9bsim+//+S2S5xb+8KFRzJyQAkAjc3Ob+asSUt8IiIiLR5asIH9DcH0msFlRVweGt1tcdTIARwzakB8+5HXN7/jGBHpPiW+ItKjzFtdTn1j0Mm5X3EBJx1ktLdF3+ICPvlvh8W3//zqRnbsq01LjCIiIgCPhSqUbjhlLCWF+W0eN/PYQ+OP39q2t81jRKR7lPiKSI8SbhBy5hFDKczv+J+xq08czbB+xQDUNzbzu7nrUh2eiIgIABt317Bie6Js+coTRh/w2IlDy+KP1+/SVByRdFDiKyI9hrvzXCjx/bfJwzv1+pLCfD4W6qb5+5fXU7m/IWXxiYiItHhm+fb440nDyhgfSm5bCz+3aU8NjU3NaY1NJBcp8RWRHmPF9n1sqUyUJ589eVinz3H9jHH0jy0dUVXXyO9fXp+y+ERERFo8uzzxQ+15U0a0e+zYwaXxxw1NztZKTcURSTUlviLSY4TLnI8dPYDh/Uo6fY6+xQV86LTx8e3fzV3L/nqt6ysiIqmzr7aB+WsTjRjPPUjiW1KYzyH9E59pKncWST0lviLSY8x+a2f8cWfLnMNuPH0CfWINRnZV1/OX17RuooiIpM6cleU0NDkAg0oLOWHswIO+ZtyQxKjvul3VaYtNJFcp8RWRHqGypoHXNuyJb7/ryK4nvoPLinhfaBmkB19aj7t3Kz4REZEWz4bm9/7b5OEUdKARYzjx3aC15kVSTomviPQIz7+9k6bmIDkd2rcoac3DrvjQqePjj9/eUcWr6/cc+GAREZEOamxq5l8rElNzDlbm3GLckESDq3XlGvEVSTUlviLSI4S7OZ99xHDy8qxb5xs/tIzTDxsS335o/oZunU9ERARg4YYKKmqCFQMK842zjhjaodeFR3w1x1ck9ZT4ikjWc3fmrEzM7+1OmXPYdSePiz9+cslW9lTXp+S8IiKSu8JlzqdMGEK/ksIOvW58aMR3/e5qTcERSTElviKS9bZU1rIrlJSGR2q74/ypIxjatwiA+sZmZi3clJLziohI7gqv33vulI7/UDs2NOJb29DMjn11KY1LJNcp8RWRrLdsy97441ED+zCwtCgl5y0qyOO90xNNrh5asEG/sIuISJdt2FXD6p2J+bkHW783rH9JIYPLEp9vKncWSS0lviKS9cKJ75RD+6f03O8/eSwWmy68Zmc1L6/ZndLzi4hI7pi3ujz+eOKwMsYMLm3n6HfSkkYi6aPEV0Sy3vKticR36sjUJr5jBpdy5uHD4tsPLVCTKxER6Zp5q3fFH58+qWNNrcLGDQ43uFLiK5JKSnxFJOstCye+KR7xBbju5LHxx//35jYqatTkSnKXmV1iZi+b2R4ze9rMjmrn2LFmdp+ZNZnZ19t4/kwz8zZuxaFjRpjZX81sp5ktNrNPpuu9iaSTuyclvqdN6nw/ivCSRip1FkktJb4iktX21TawYXfiwz8die+5U4YztG/wPby+qZm/Ld6a8muI9ARmNhN4GHgQmAGUA3PMbGQbxw4FVgB9OPD3iQnAMmBg+ObudbFzlALzgGLgTODHwP+Y2S0pfFsiGfH2jirKq4KGVGYwY2JXEl8taSSSLkp8RSSrvbVtX/xxv+ICRg/qk/JrFObncdlxie/1s15Td2fJWbcDs9z9F+6+AriR4LvCza0PdPdyYIa7X9vO+SYCb7t7ZfgWev4jwCjgA+7+lrvfB/wW+KKZ6TuK9CjzViXm9049tD+DyjrfiDE84rtul5Y0EkklfaiISFZr3dgqL8/Scp2rThgdf/zGxgpW76xKy3VEspWZTQCmA4+07HP3WmA2MLOt17j7ooOcdgKwqp3nrwGec/eK0L5/xF435eBRi2SP7pY5A4wPjfjuq22koqah23GJSECJr4hktXDim+rGVmFTR/ZP6hj9sNb0ldwzLna/ptX+DcBoumYi8HEz225mC8zsU21cs63r0ZFrmtnStm7ApC7GK9IlTc3Oy2vCiW/nG1sBDC4rom9xQXx7/W6VO4ukihJfEclqy7eFR3z7pfVaV50wKv74kYWbaW5WiZnklENi9xWt9lcBg7p4zk8AJwJXAXOAn5nZna2u2db16MY1RTJu2Za97K1tBKAgzzhpwuAuncfMWs3zVWdnkVRR4isiWauxqTlpju/UQwek9XqXHTeK/Fgp9ZbK2qRf70VyQMsExdalFQOAbV05obu/6e4r3X2uu38e+AnwKTNr+ZWp/ADXoyPXdPej2roBq7sSr0hXhdfvnTZmYNKobWclreVbrhFfkVRR4isiWWtNeTX1jc0A5OcZh4/om9brDetXzNlHJNb0/avKnSW3tPwHP77V/nGh57prDmDAiNA127peOB6RrJeK+b0tkpY02q0RX5FUyZnEt6PrEprZRbHnt8aO/buZHdHGcV8xs+Vmts3MHjSz9A5FieSg5aH1eycNK6OkMD/t17wyVO781JvbqK5rTPs1RbLEitjtspYdZlYGnAPMStE1pgN1BEscATwGvMvMwr9qXQwscff2mmKJZI36xmZeWbc7vn1qNxPf8VrSSCQtciLx7cy6hMBY4HGC9QTPJ1hm4R9mFu9Jb2Z3AJ8HPk3wAT0NeCqd70EkFyU1tkrD+r1tOW/KCPqVBCVqNfVNPPVmlyo8RXocD9ZN+Q5wvZndZGaTgXuBSuBeM3u3me0ysws7cj4zO9LMbjGzKWY21sw+DXwZ+I9Yt2iAuwnm9D5gZpPN7CPADcC3U/3+RNJl0aYKauqbACgqyOOEsd2bnj52cGjEV4mvSMrkROJL59Yl/LW73+nuq9z9VeCbBF0ppwKY2TDgVuCr7v5M7JibgRlmdk5G3o1Ijli2NTMdncNKCvO5+NjEb2KPLdqSkeuKZAN3f4DgM+1W4EWgBDgtttxQM8FobUfXVxkIXEdQ3vxm7PH17v7j0PXKCX5o7gfMJ2iG9UF3/3NK3pBIBsxblShznj5uULerk8JzfMur6qhtaOrW+UQk0PWZ9z1EaF3CH7bsc/daM5tNsC7h7Qc5Rcu37T2x+4uBPoTWOST4sN4dO9/s7sYsIuDurUZ8Mzeb4PLjRvLHBcGKKnPf3snOfXUM61ecseuLRMnd7wHuaWP/M0BblVK4+zsW2Hb3l4GzOnC9lQQVViI90vy1icT31IndK3MGGNK3KGl7T009hw7o0+3ziuS6XBjx7fK6hGZ2GvAN4F53Xx86X427b285LlYetulg54udU2sOinTAzn117Kquj2+neymjsJPGD2bkgBIAmh2eWKxRXxEReaf6xmYWbtgT3z4lBYlvcUE+pUWJUeM91R0tshCR9uRC4tvpdQnN7Etmtp+gzGspQclX+Hytz9Xu+USk85aGypxH9C9mSN/Mjbjm5RmXHJcY2Hr0DSW+IiLyTks2V1DbEKw+UFSQx7QxqalOGlSaGPWtqKlv50gR6ahcSHy7si7hrwgaVl0BDAFeM7OhofO1NdmwQ+scas1BkY5ZvaMq/njyIZmZ3xt2+XGJ7s6LNlawtlxLSoiISLL5axPdnI8fM5DigtSsPjCwtDD+eE+NRnxFUiEXEt9Or0vo7pXuvtLdHyWYdzQB+PfQ+fqaWbyWxczyCcqcteagSIps2J3oZDkh1OgjU6Yc2p/JIxLl1Y+9sTnjMYiISHabvyaR+KaizLlFeMR3j0Z8RVIiFxLf7q5LuJ+gi2XLX6Angabw+YAzCEZ8H+5+uCICyYnvmMGZT3wBLjs+Ue78+BtbCKbzi4iIQGNTM6+tD83vnTA4ZeceEBrxVamzSGr0+sS3s+sSmtkXzewsMxtpZscDf4yd6sHY+TYC9wN3mNk5ZjYd+CXwd3dfmOG3J9JrhRPfsRElvpdOSyS+a8qrWbK5MpI4REQk+yzbupequkYACvKs2+v3hg1SqbNIyvX6xBc6vi6hmRUDhxEs47AaeBxw4Ax3Xx465ceB3wMPEIwazwauysibEckBTc3Opt3749tjIyh1Bhg9qJSTxie+yDz6uppciYhIYEFofu+xowfQpyg183uhdXMrJb4iqdDr1/Ft0Yl1CW/uwLkagdtiNxFJse17a6lvao5vjxkUTeILcNlxo3hlXVDK9rfFW/jqzCnk571jyVIREckxL6dpfi/AQHV1Fkm5nBjxFZGeJVzmPLRvMWXF0f1GN/OYQymIJbo799Xx8ppdkcUiIiLZobnZeWVdKPFN4fxeaF3qrMRXJBWU+IpI1kme39snwkhgUFkRZx0xLL6t7s4iIrJi+z4q9wclyHkGJ45L3fxeUKmzSDoo8RWRrLNhV/SNrcIuOy4xG+Ifb26jtqEpwmhERCRq80PVP0ePGkC/ksJ2ju68gRrxFUk5Jb4iknWyoaNz2HlTRtCnMGhasq+2kdkrdkYckYiIRGlBqMz55PGpLXOG5BHfyv0NNDdrOT2R7lLiKyJZJxvW8A0rKy7g/Kkj4tt/W6TuziIiucrdkzo6p7qxFSQnvs0Oe2tV7izSXUp8RSTrbAwlvuOGlEUYSUK43PmZ5dvZpy8hIiI5afXOasqrgvJjM5KWvUuVfiUFhBcQ0Fq+It2nxFdEskpVXSO7qhPzmbKh1BngzMOHxedc1TU288+l2yOOSEREohAe7Z08ol/S0kOpkpdnDOijeb4iqaTEV0SySrixVVFBHsP7FUcYTUJRQR7vOebQ+PZjKncWEclJC9YmGlulehmjsEFay1ckpZT4ikhWSZrfO6gPeeFar4hdOi1R7vziqnJ27quLMBoREck0d2d+aMT35Ampn9/bIqmzc7VKnUW6S4mviGSVbJzf2+Lk8YM5dEAJAE3NzpOLNeorIpJLNu3Zz9bK2vj2yZka8d2vxFeku5T4ikhWybaljMLy8oxLQqO+KncWEckt4dHeicPKGJbG6TgDVeosklJKfEUkq6zPsqWMWgt3d359QwXrd1VHGI2IiGRSpub3AgwqVXMrkVRS4isiWWVjFo/4Akw9tD+HD+8b337sDY36iojkigVJ83vTnPiWJUZ8tZyRSPcp8RWRrNHU7Gzak92Jr5lx+fGj4tuPvrEZd48wIhERyYTte2tZF1p54JQ0NraC5OZWKnUW6T4lviKSNbbtraWhKZFEZmPiC8ndndfsrGbplr0RRiMiIpkQnt87elAfRg7sk9brhZtbqauzSPcp8RWRrBGeLzusXzF9ivIjjObAxgwu5cRxg+Lbj76+OcJoREQkE8Lze9Nd5gwwsI9GfEVSSYmviGSNbJ/fG3Z5qMnV44u20NSscmcRkd4sPL833Y2tILmrs+b4inSfEl8RyRrZvJRRazOPHUlBngGwY18d89fsOsgrRESkp9pdXc/K7VXx7XTP7wUYVJYY8d3f0ERtQ1ParynSmynxFZGssWH3/vjjbE98B5cVcdYRw+Lbj6jcWUSk1wqP9g7vV8y4Ien/jArP8QWo0KivSLco8RWRrLEhy9fwbS28pu9Tb27Tr/EiIr3Uy2uS5/eaWdqvWVKYT0lh4qt6xX7N8xXpDiW+IpI1tlUmRnxHDiyJMJKOeffUQyiLNeDaV9fI08u2RxyRiIikQzjxPXVS+sucW6izs0jqKPEVkazQ1Ozs3FcX3z6kf/Ynvn2K8rnw6EPj2yp3FhHpffZU1/PWtn3x7RkTM5f4hhtcqbOzSPco8RWRrFBeVUe4MfKIHpD4Alx5wqj44+dX7qS8qq6do0VEpKcJr987rF8xE4eWZezag0oTDa7U2Vmke5T4ikhW2FZZG3/cr7iAsuKCCKPpuBkTh8RHp5uanScWbYk4IhERSaVwmfOMiUMyMr+3RVKps0Z8RbpFia+IZIXtexOJ74gBPWO0FyA/z5KaXKncWUSkd0lOfNO/fm/YwNCIr0qdRbpHia+IZIVw4tsT5veGXREqd160qZLVO6vaOVpERHqK3RHO74XkxFelziLdo8RXRLLCtlDiO7x/cYSRdN6Rh/RnyqH949uPLNSor4hIb7BgbWK0N9PzeyG51FkjviLdo8RXRLLC9r09q6Nza1cenxj1feT1zTSHO3WJ9CBmdomZvWxme8zsaTM7qp1jx5rZfWbWZGZfb/VcoZndYmYvmVmFma01s5+aWWmr4541M291+2663p9IZ7y8JtHYKtPzeyG5q7NGfEW6R4mviGSFpFLnHjTHt8Wlx40kL/Z9aHPF/qQuoCI9hZnNBB4GHgRmAOXAHDMb2caxQ4EVQB/a/j5RBFwHfB84Fvgs8AHgrlbHTQA+CAwM3b6Rgrcj0m1Rzu+F1l2dNeIr0h1KfEUkK4QT3+H9el7iO6J/CWcePiy+/dfXNkUYjUiX3Q7McvdfuPsK4EaC7wo3tz7Q3cuBGe5+bVsncvdqdz/L3R9x9w3u/hjwC+B6MysCMLMCYAywxN0rQ7fats4pkklRz++F5BHfSo34inSLEl8RyQrh5Yx64ogvwNUnjo4//sebW6mua4wwGpHOMbMJwHTgkZZ9sQR0NjCzrde4+6JOXqYcKABayp3HxLZXd/I8ImkXnt87PIL5vZA84luxvwF3TaMR6SolviISuf31TeytTSSJPXGOL8D5U0fQryRYf7imvom/L9kacUQinTIudr+m1f4NwGhS43RgrbtXxLYnxu5Xm9kaM/uDmR3R0ZOZ2dK2bsCkFMUrOSzq+b2Q3NyqqdmTPitFpHOU+IpI5MJlznkGQ/sWtXN09iopzOfSaYmpkCp3lh7mkNh9Rav9VcCg7p7czE4ErgZ+GNr9BnAM8C7gPwjmAs81syO7ez2R7pq3ujz++JQI5vcC9O9TSDjfVmdnka7LmcS3o10qzewUM3vYzDaEjj2m1TEfaKMDpYZ2RLoovJTR0L7FFOT33D9N4XLn+Wt3s2FXTYTRiHRKy7f8/q32DwC2defEZjaGoIT6/4DftOx3913u/mbs9ifgPIIy6P/XkfO6+1Ft3VDptHTTjr21rNyeWJP99ElDI4kjP88Y0Edr+YqkQs/9dtkJnelSCbwfWABcAJwFODDbzMI/9U0AniK5A+XktL0BkV6up3d0DjtuzEAmDUvMA5u1UKO+0mO0/Mc6vtX+caHnOs3MhhAkvFuAa9z9gLWa7r6doFP0oV29nkgqzFudmN87amAfxg0pbefo9BrYR52dRVIhJxJfOtel8rPufoe7L3f3JcAtwGDg8tBhE4EVrTpQ7k3/2xDpnXp6R+cwM+PqE8fEt2ct3KQ1faWnWBG7Xdayw8zKgHOAWV05Yez1TxJ85l7s7tUHOX4AcDjweleuJ5Iqc1clypxPPyya+b0twp2dVeos0nW9PvHtSpfKVlr+8vUL7ZsArEpRiCI5b1tlXfzxIQOKI4wkNa44flR8Td9Ne/bzcqgzqEi28qBd7HcIlhu6ycwmA/cClcC9ZvZuM9tlZhd25HxmVgo8DowCrgEazGxA7FYYO+aTZnaOmY00sxmx4/cAP039OxTpGHdnXlLiG02Zc4ukzs4qdRbpsl6f+NL9LpWnx+7Dvz5PBL5lZlvN7AUze183YxTJaUmlzj20o3PYIQOS1/T9y6sqd5aewd0fIKiGuhV4ESgBTot1YW4G6oCOfvP+DEHTqtHAYoKmWS2362LHHAH8lmBO7kPASuBkd9+Tivcj0hVry6vZElpi77SI5ve26FeSSHyr1NVZpMsKog4gA7rcpdLM8oFvAa+4+5zQU5cC1QS/Yn8Q+KOZjXL3Hx0smNgyC23R0guSs5JKnXtB4gtw7UljeH7lTgD+vmQrt19yFANCv9qLZCt3vwe4p439zwBt9cbA3d9RB+ru3wW+e5BrfRb4bNciFUmPF0PzeyeP6MewftFWIpUV58cfV9c3RRiJSM+WCyO+3elSeTdBWfON4Z3u/rq7r3T359z9RuCvwDdaSrdEpHO29bIRX4DzpoxgcFkwL6uusZnHFm2OOCIREemIF9/OnjJngLKixDhVdZ1GfEW6KhcS3y51qTSz7xN0eL7U3Q80SttiDlBG8jzgNmnpBZFk7s6OveE5vr0j8S0qyOOqE0bFt/+4YCPBFEoREclWTc3OS2sSI76nHzYkwmgCpcWhxLdeia9IV+VC4tvpLpVm9nngNuD97j63A9eYDmxw993djlYkx+ypaaC+qTm+PaKHd3UOu/akRHfn5Vv38uZmNX8XEclmS7dUUrk/mMaen2ecMjH6xLdvuNRZI74iXdbrE9/Odqk0s08A3yNoyvFcqANlWez5083sBjM7zMwmmtntwAcIEmUR6aRtoQYiJYV59O/Te1oPHDa8H9PHJVoJ/OmVDRFGIyIiB/PiqsRo73FjBtK3OPrPpNJQqXON5viKdFmvT3yh410qzWwk8AuC/11+TnIHyidipxsUO88rsdtZwLnu3qU1DkVy3fZ9yfN7o1wrMR3Co76Pv7GFGpWpiYhkrRezaBmjFuHku0ojviJdFv3PWBnSiS6V7f4Y4O5PkEiCRaSbtlf2vo7OYTOPPZT//tsy9tU1sq+ukScXb+Wa6WMO/kIREcmo2oYmXlmXmLV2+qToy5wBSosSpc41dRrxFemqnBjxFZHs1Rs7OoeVFhVw6XGJ39b+9MrGCKMREZEDeXXdHuoag54TfQrzOX5su6teZoxGfEVSQ4mviERqey/s6Nza+04aG3/82vo9vLVNTa5ERLLN8yt3xB+fOmkIRQXZ8X0JX7YAACAASURBVDW5rDg8x1eJr0hXZce/aBHJWdtDI77D+xVHGEn6HDN6AMeMGhDf/v3L6yOMRkRE2vL8yp3xx2cfMSzCSJKVJXV1VqmzSFcp8RWRSIUT39464gtww4zEqO8jCzerXE1EJItsqdjPyu1V8e1zJmdT4psY8a1vaqa+sbmdo0XkQJT4ikiktvfyOb4tLpk2kn4lwZeX6vomHntjc8QRiYhIi/Bo7/ghpYwbUhZhNMnCyxmByp1FukqJr4hEpr6xmfKq+vj2iF6c+JYWFXDVCaPj279/eQPBMuMiIhK151dkZ5kzQFmoqzOowZVIVynxFZHI7KyqS9oe3r93zvFtcf0piXLn5Vv3snBDRYTRiIgIQENTc9L6vedMHh5hNO9UkJ9HcajRVk295vmKdEW3E18z078+EemScJnzoNJCigvy2zm65zt8RD9OmTA4vv2H+WpyJSIStYXr97AvNopaVJDHKRMHH+QVmacljUS6LxUjvpaCc4hIDtoVKnMe2rd3j/a2uGHGuPjjJxZvZU91fTtHi4hIuoXn954yYfA75tRmg6QljdTZWaRLUpH4apKaiHTJrlCp85C+RRFGkjkXHHVIPMmvb2zmT69sjDgiEZHclq3LGIWVhub5asRXpGs0x1dEIrMrNNo5JEdGfIsK8rju5DHx7QdfWkdjk5amEBGJwo59tSzdsje+nU3LGIWFS53V1Vmka5T4ikhkykMjvkPLcmPEF+D6GeMoyAtmiWyprOXpZdsjjkhEJDfNWZloajVqYB8mDesbYTQHVhpKfKs14ivSJUp8RSQy4Tm+uTLiC8GyTe855tD49r3z1kUXjIhIDntuxY7447OOGIpZdrau6VscLnXWHF+Rrsh44mtmJWb2dTP7m5n93sy+aWZXmNm4g79aRHqTXdW5N8e3xYdPHx9/vGDtbpZuqYwuGBGRHFTf2Jy0fu+5R46IMJr2hRtuqdRZpGuiaFt3L3Btq30OYGYVwEJgobt/KdOBiUhmJY34luXOiC/A8WMGMm30ABZtChLe+15cxw+umRZxVCIiuWP+2l3xRlElhXmcftjQiCM6MC1nJNJ9UZQ6zwSeBCYD/YFpwIeB/wGWACcBn48gLhHJsPKk5Yxya8TXzJJGfR9btCWpy7WIiKRXuL/CGYcNo09R9q4lH+7qrOWMRLomisS3Apjl7m+7e5W7L3H3B939/7n7Oe4+EDg8grhEJIOam53dSaXOuTXiCzDzmJEM65dY2uiPCzZEHJGISG5wd54JJb7nTx0eYTQHF17Ht0qlziJdEkXiOwu4sL0D3H1NhmIRkYhU7G+gObQKeK7N8YVgaaPrTxkb377/pfXUNeqXfBGRdFu2dS9bKmsBMIN3ZfH8Xmi1nJFKnUW6JIrE92vAFDP7agTXFpEsES7rLcrPo19xFC0Honf9KeMoKgj+FO/cV8djb2yJOCIRkd7vmWWJbs7HjRkYr77JVuFS52qVOot0SRSJbwPBqO/tZva6mX3OzE4ws9wb7hHJYeVJSxkVZe0SEuk2rF8xV50wOr5995w1uHs7rxARke56Znm4zDm7R3shecS3WqXOIl0SReL7APB1IJ+gsdUPgFeAKjNbbGb3m9lnI4hLRDIol5cyau1jZ06IP357RxWzV+5s52gREemOrZX7WbI5sYTc+VOyP/EtDSe+KnUW6ZIoEt8LgGeB44GxwNnA/wMeij1/HfCjCOISkQzK5aWMWps0rC/nhb543T1HbQ5ERNLl2eWJMudxQ0o5bHjfCKPpmL7FoVLnepU6i3RFFJPqVgF/dPdFse1NwAstT5pZCXBsBHGJSAaF5/jm+ogvwM1nTYyX3s1bvYs3N1dy9KgBEUclItL7hMucz5syokdMtSkt0oivSHdFMeL7FMGob5vcvdbdF2QwHhGJQHl1eA3f3B7xBThp/CCmjRkY3/6NRn1FRFJub20D81btim+f1wPKnKFVV+f6Jpqb1QtCpLOiSHy/AxxmZl+J4NoikiWSRnzLNOJrZtx85sT49pNLtrJxd02EEYmI9D7PLt9OfVMzAINKCzlp/KCII+qYcFdngJoGlTuLdFYUie9CoB74lpm9aGY3mtmYCOIQkQglzfHViC8AFxw1grGDSwFoanZ+9fzqiCMSEeldnly8Lf74gqMOoSA/iq/CnVfWask/lTuLdF4U/9rfAgYCjcCpwD3AOjPbZWb/MrOfmNmHIohLRDJoV3XyckYCBfl53HL2pPj2X17dxI69tRFGJCLSe+yrbWDO24mu+TOPPTTCaDqnuCCPgrzEXGQlviKdl/HE190vd/epQClwNHADQRfn12LbnwF+l+m4RCSzykOlzkNzvKtz2FUnjmJE/+B/j/qmZu5+QXN9RURS4dnlO6hvTJQ5nzpxSMQRdZyZJZU7V9ep1Fmks9Ke+JpZaVv73b3J3Ze5+0Pu/kV3f7e7DwdGA5ekOy4RiU5dYxP7ahO/VmvEN6G4IJ+bQnN9/zB/A3tCo+Mi6WZml5jZy2a2x8yeNrOj2jl2rJndZ2ZNZvb1AxzzFTNbbmbbzOxBMxvQ6vkRZvZXM9tpZovN7JOpfk8iEPROaNGTypxbhBtcVddrxFekszLxL36bmd1lZtM6crC7b3H3v6c7KBGJzu5WidxgNbdKct0pYxlUWggE3Tvvnbcu2oAkZ5jZTOBh4EFgBlAOzDGzkW0cOxRYAfThAN8nzOwO4PPAp4GLgWkEqzu0PF8KzAOKgTOBHwP/Y2a3pO5diQRlzs+vTJQ5v+eYnlPm3KK0WEsaiXRHJhLfrwFnAQvNbH6smVWfDFxXRLJUuLFVv+ICSgrz2zk695QWFfCR0yfEt+97cS1V+pIjmXE7MMvdf+HuK4AbCb4r3Nz6QHcvB2a4+7VtncjMhgG3Al9192fc/dXYeWaY2Tmxwz4CjAI+4O5vuft9wG+BL5pZzxqOk6wWLnMeWFrIqZN6Tplzi7KkEV+VOot0Vto/VNz9Z+5+NEHy+xbwc2Crmd1pZsek+/oikn3C83tV5ty2D546Pl7Wtre2kQdeWhdpPNL7mdkEYDrwSMs+d68FZgMz23qNuy9q55QXE4wGPxLaNx/YHTrfNcBz7l4ROuYfwARgSufegciBJZU5Tz2Ewh5W5gxQljTHVz+GinRWxv7Vu/uL7v4hYCTBKPDZwBtmNs/MPmhmJZmKRUSipaWMDm5AaSEfOHVcfPs3c9awr7YhwogkB7T8B9e6o9oGgv4bXTlfjbtvb9nh7g5sCp1v3AGuR0euaWZL27oBkw72Wskd7yhz7kHdnMPKVOos0i1RdHWudPc73f1Y4AyC+UF3AVvM7GdmNjXTMYlIZu2qDo34an7vAd105sT4L/wVNQ3cr7m+kl6HxO4rWu2vAgZ18Xytz9X6fG0dUxW778o1Rd6hdZnzaT2wzBlaj/iq1FmksyKt83D3l9z9RpJHgZeY2VwzOy2V1+pol0ozO8XMHjazDaFjj2l1TIGZ/dDM1saOu9PM9O1dpIM04tsxg8uK+PDp4+Pbv5mzhr0a9ZX0KY/d92+1fwCwrYvna32u1udr65iWrs8Hvaa7H9XWDVjdhXill3r49c3xxxce1TPLnKH1HF+N+Ip0Vlb8y3f3vbFGGscBpxHMBT41VefvTJdK4P3AAuACgnnJDsw2s8GhY34PXAlcC3wg9viBVMUr0tuVhxNfjfi266YzJybN9b137rpoA5LebFPsfnyr/eNCz3X2fH3NLD68Zmb5BCXMm0LHtHW9cDwiXbZjby1z306UOV9+/KgIo+mevip1FumWtCW+ZnaGmT0bm7/b4SEdd5/v7h9z9x+lMJzb6XiXys+6+x3uvtzdlwC3AIOBywFio7/XAp909wXu/jzwOeC9Zjax9flE5J2SSp3V3KpdA0uL+MgZiQ7Pv527hsoajfpKWqyI3S5r2WFmZcA5wKwunO9JoCl8PoIpTgMIfowGeAx4l5n1DR1zMbDE3Vd14ZoiSR5ftIVmDx6PGtiHk8cPbv8FWay0SImvSHekc8Q3n6Aj433AJjP7rpmNa/8lqdeVLpWttJR+9YvdXw3sA54JHfMPwICLuhmuSE5QqXPnfPSMCfQrCb7w7Ktt5LdzW/cCEum+WOOp7wDXm9lNZjYZuBeoBO41s3eb2S4zu7CD59sI3A/cYWbnmNl04JfA3919Yeywuwnm9D5gZpPN7CPADcC3U/vuJFfNWpgoc77i+FHk5VmE0XRPWXFojq+WMxLptLQlvu7+vLuPBI4l+CCdBiw2s8fM7Lx0XbcN3e1SeXrs/vXQ+Ta6e/ynNnevBPZ25HzqQCkCu0LLGQ1VqfNBDehTyMfOSBSU/G7u2qQloURSxd0fIKiGuhV4ESgBTostN9QM1AGdKTn4OMH0oAcIRo1nA1eFrlcOnEnw4/J84BPAB939z919LyLLt+5l+da98e0rTui5Zc6grs4i3VVw8EO6x93fBN4EfmJmhQSJ5HSSR0zTqctdKmNzkb4FvOLuc0LnO1iXShE5AHenvFojvp114xnjuXfeWipqGqiub+Ln/1rF7Ze22aNPpFvc/R7gnjb2P0PQjLKt17Q5jBb7kfi22O1A11sJnN+lYEXa8UioqdW0MQOZNKxvO0dnv+TmVhrxFemsjDa3cvcGd5/t7j/M4GW706XybmACwZzg8PkO1qXygNSBUnJdVV1jfFkJ0BzfjupfUsinzjksvv2H+evZuLsmwohERLJXU7Pz2BuJxPfKHtzUqkXyckYa8RXprKzo6pxmXepSaWbfJ+jwfKm7L211vtGx0eCWY4cAZe2dT0QC4fm9ZjCoVIlvR33g1HGMHFACQEOT8+OnV0YckYhIdpq3upzte4MpIQV5xiXT2ixW6FHCI741SnxFOi0jia+ZDTGzj5vZ58xsSiauGdLpLpVm9nmCsqz3u/vcVk8/BgwkWHO4xcUEc56eSFnUIr1UuKPz4NIi8ntwo5FMKynM57PnHxHffvSNzUnz10REJPBwqKnVOZOHM7gX9JMIL2dUpcRXpNPSnvia2QjgVeAu4AfAG2b24dhzw83sJ2Y2x8z+2tFOkZ3R2S6VZvYJ4HvAZ4DnzGxA7FYWO99LwL+AX5rZdDM7O3b83e5+0FJnkVyXtIavypw77aoTRnPEiGCemjt8/6m3Io5IRCS7VNY08PclW+PbV/bwplYtSouSuzoHX3FFpKMyMeL7dYKGUF8CriFY1++XZnY0MIcgwTwDuBJ4MjbamlId7VJpZiOBXxD87/JzgiZWLbfwaO5lwEux9/Ib4LfAp1Mdt0hvlLSUUZkaW3VWfp7xhQuOjG8/t2In81aXt/MK6QnM7BAzu8HMrjSz4W08f2cUcYn0RA+/vom6WC+JwWVFnDvlHf+keqTwiG9Ts8ffo4h0TCYS33OBu9z9h+4+y92vBJ4GHidY/ud6gtLhGQSJ8B2xUdmUcvd73P0Ydx/q7pe6+4bY/mfcfaS7P+vuW9w9z92tjdu/hc5V5e4fdvcR7j7Z3f/T3fXXR6QDwksZacS3a86bMpzp4xJN5L/1xHKamvXLf09lZucSNDj8GcHSP+vM7FtmFp4HcF0kwYn0MO7OHxdsiG9fc+Joigvy23lFz1FanLwYixpciXROJhLfMcDCVvt+StBs6hfu/kd33+vuC4CLgC0EI7Mi0gvtCi1lNFRLGXWJmfEfMxPtEpZt3cus19Rbrwe7A/iSuw8B+gI3ECS6j5hZSewYTYYX6YDX1u9h5faq+Pb7Th4bYTSpVVqYnMDXaEkjkU7JROJbzTvXC14eu38hvNPda4EHgVMzEJeIRGB3KPHtDc1GonLC2EFcdlyiS+kP/rlCzU56rkkEfTBw92Z3f5igCmokMCdW+qwhfZEOeGh+YrT3tElDmDC0LMJoUisvz5Lm+epvvkjnZCLx3QZ83sx+aGbXxcqYW2odK9s4fiVwaAbiEpEI7KlJJL6DSgsjjKTn++KFR1JcEPwZ37mvjl/OXhVxRNJFmwmm/sS5+w6C1QO2Ac8DKo8QOYiKmnqeCDW1uu6U3jPa2yJpSaN6Jb4inZGJxPclYALB8kC/B5YRrHfrwKfM7EYzO97MWoZ+6oHBGYhLRCJQUdMQfzxAa/h2y6iBfbj5rInx7btfWMumPTURRiRd9AfgE613uvt+4HJgLtAn00GJ9DQPL9xMfazh05CyIt499ZCII0q9sqQRX5U6i3RG2hNfd78Z6AdMIZiz9EOCEuddwHuBewiWO9pnZouAzwH70h2XiESjYr9GfFPplrMnMbxfMBhY39jMd/+h5Y16oB+4+1faeiJW+nwTweeliBxA66ZWV08fTVFBJsZ3MitpxFelziKdkpG/CB5Y4e5/cvcvufsF7j6coPHVpcDtBMsF9QNOAPZkIi4RybyK6sSI7yCN+HZbWXEBX7gg0Qj/ycVbeXGVljfqSdz9oMM27j4rE7GI9FTz1+7m7R2JplbvP6n3lTlDcuKrOb4inRPpT2Huvtndn3D3b7r7Ve4+ERiEftkW6ZUamprZF/qgHtBHI76pcNUJozluzMD49tcfXxov95OezczaHAkWkWT3zF0bf3zGYUMZ34uaWoWFS53V1VmkcyKvATGzk8zsQjPLB3D3Snd/Peq4RCT1Kvc3JG0PUlfnlMjLM7552dG0rPq6akcV9764tv0XSU/xbTO75UBPWuCGTAYkkm3WlVfzzPLt8e2PnjEhwmjSSyO+Il0XaeJrZh8H5gF/Bv4e2/dlM3vdzB42s7OjjE9EUqsi1NG5IM+SfrmW7jlm9ACuC61X+bNn32ZbZW2EEUmK3Af83Mze3/oJM7sWWArcn+mgRLLJvS+uxWMLfk0cVsbZRwyLNqA0KitKJL7VSnxFOiXqEd/bgPcQlDevMbP7gZsIEuGtwBNmpjV9RXqJcEfngaVFWMsQpaTEFy6YHG8YVlPfxLeeXBZxRJICHyXo+nyfmc0EMLMrzWwx8BDBWr/fjjA+kUhV7m/gL69tim9/5PQJ5OX13s+W5OWMVOos0hlRJ74jgGfdvRn4b+ADwEfc/Q53/xTwaUDzm0R6iT014cZWmt+bagNLi/jyRUfGt59YvJXnV+6MMCLpLDNLWq/X3R34MMEPwn+JrX7wV4J1f/8bGOfu/5XpOEWyxZ8WbIgngANLC7nqhNEHeUXPVlYcXs5II74inRF14lsLnG5mZe6+FagmWK+wxUPAtEgiE5GU2xMqdR6oxDctrjlxDMePTTS6+uojS6ip15ejHmSvmb1mZr82s4+Z2TSCz+oPAbOAY4AfAOPd/RvuXhllsCJRamxq5v556+Lb158ylj69fApN8oiv/raLdEbUie9zBMsYVZjZm0Ah8Ekze5eZDXP3eoIljkSkF6hsVeosqZeXZ3znimMoiJX6bdqzn588vTLiqKQTngVGEUz7+Q2wENgLvECw1N8GgmX/qqMKUCRb/OPNbWyJ9TIozDc+eOr4aAPKgHBvjKo6lTqLdEbBwQ9JH3d/P4CZTQSOi93OB74AjDKz3YB+zhLpJcIjvip1Tp8ph/bn5rMmctfs1UCwzMdlx43i6FEDIo5MDsbd3wNgZqOA6cCJofvTYoeNBXbHyp5fB1539/syH61IdNydX8b+xgFcfOxIRvQviTCizEga8VWps0inRJr4tnD3NcAa4OGWfWY2mCAR7p0rkIvkoIr9GvHNlFvPPZy/L9nKul01NDt8adZiHvvU6RTkR13oIx3h7puBzcBjLfvMbAxBAhxOhs8AnKD7s0jOeHb5DpZt3RvfvunMiRFGkzlazkik6zL+DcjMro6N8LbL3Xe7+7/0K7ZI71GhOb4ZU1KYz3euOCa+vXTLXn47V2v79mTuvtHdH3X3r7n7Re4+HBgPXB1xaCIZ5e7c+dyq+PZ5U0YwdWT/CCPKnPByRurqLNI5Ufz0/7/A22a2x8z+ZWY/NLPrzWyqmfU1s5mx9X1FpJfZUx0a8e2jEd90O+2woVxzYqLD6Y+fXsnb2/dFGJGkmrtvcPdHoo5DJJNeeLucRRsr4tu3nntYhNFkVp+ixFf3/Q1KfEU6I4rE9+PAr4DlwMkEa/k+ACwBKoHHgTsiiEtE0ixc6qw5vpnx1ZlTGN4vWCGnvrGZz/9lEY1NzRFHJR1lZkeY2fNmVmNmO2Mdn+82s1vM7GQz6/2TGkVC3J07//V2fPvsI4Zx7OiB7byidykpTDS3qtWIr0inZDzxdfe73f1T7n4aQcfmqcDNwEuAAfcDl2Y6LhFJv+RSZ434ZsLA0iK+e2Wi5HnRpkp+PWdNhBFJR5jZZ8ysFPgdcCawE1gPjAE+CtxF8Lm5N9bkSiQnvLxmN6+s2xPfzqXRXoDScKlzQxPBUt8i0hGRdjnxwFvufo+7nwF8nWAUWBPRRHohreMbjXOnjOCqExIlzz99ZiUrtqnkOcv9FzAEOB74ubuPc/fpsXm9Y4HLgW8B/wSGRhemSGaFR3tPnTiEE8cNjjCazOsTGvFtanYampT4inRUVrX3dPdvAgsIfskWkV6ktqGJ2oZEie0gjfhm1H9dMpUR/YOS54Ym57b/fYP6RpU89wAbgcXhHe6+yd0fd/evu/vF7j4qothEMmru2+XMW70rvv3pHBvtheTEFzTPV6QzoujqbAc55CXg7EzEIiKZU1HTkLStEd/MGtCnkDuuOja+vXTLXn789MoII5K2mFlxq12/Ba6NIhaRbNLc7Hzvqbfi26dMGMypE4dEGFE0+hS1Snw1z1ekw6IY8a0ys3lmdpeZfczMTjSzIoDY/TWAvo2J9DIV+xNlziWFeUkNOiQz/m3ycK47JbE0+q/nrOal0OiJZIW9ZrYQ6Au8n6CUebuZ/cTM9I9GctaTS7ayZHNlfPvLFx3JwcdSep/CfCM/L/G+NeIr0nFRJL5zgUnALcBvCEqbq8xsI7CbYLT3O/qAF+ldwksZqcw5Ov85cwoTh5YB4A63/e8bVLYajZdIPQ0MAwoJVjh4HbgCuBVYaGa3mtmRHaieEuk16hub+eE/V8S3Lzr6EI4fOyjCiKJjZknlzhrxFem4KLo6X+DuI4DRwCUEDa3+BjQCpQQf9n8Fas1sg5m9aGZ/zHScIpJa4Y7OA/qozDkqpUUF/Ox9x1MQGzHYWlnLfzy6RJ1Bs0Rszu4YoAL4CEGTq38CW4BjgJ8CS4F94eqpyAIWyYA/vbKB9btqAMjPMz5/weSII4pWuNx5f0NjhJGI9CwFBz8ktczsK+7+XXffQvBB/mTouYEEHSyPB6YBRwCHAUdmOk4RSa3kNXw14hulY0YP4LZ3H8H3nwpGUJ5cvJWzDh/KtSeNPcgrJYO2A8+6+8aWHWY2ApgOnBS7nw7MAJxgLrBIr1NV18j/PJvo5HztSWOYNKxvhBFFL3nEV00KRToq44kv8G0z2+Puv2r9hLtXmNlsYJS7fyjzoYlIuoSXMhpUphHfqH38rEk8v2In89fuBuDrjy/luDGDmHxIv4gjy21mVuruNe4+pfVz7r6d4Mfi8A/GowkSYJFe6a7nVlFeFXx+9CnM57PnHh5xRNFLSnw1x1ekw6KY43sf8HMze3/rJ8zsWoISrvszHZSIpFe4q/OAPhrxjVp+nvHT9x3H4LLg/4vahmY++YfXqK5T2VzEtsXKl6d15ODY0kaPpjsokSis2VnF3S+siW9/7MwJDO9fEmFE2SFc6lxTr7/ZIh0VReL7UeAPwH1mNhPAzK40s8XAQ8BI4NsRxCUiaRSe4ztISxllhUMH9OHH703kV6t3VvO1R9/UfN9ofQ04i6CR1Xwzu9HM+kQdlEimuTu3/20ZDU3B36ORA0r4xDmTIo4qO4RHfGs14ivSYWlPfFuvSejBN6oPA38G/mJmiwiaWY0G/hsY5+7/le64RCSz9oRGfLWGb/Y4Z/LwpC+TD7++mf99dWM7r5B0cvefufvRBMnvW8DPga1mdqeZHZOJGMzsEjN72cz2mNnTZnZUO8fOMLPnzazCzF4ws7NDz51tZn6A2+9Cx93TxvNqapnj/rlsO3NW7oxv/+fFUyktimKGXvZJam6lrs4iHZaJEd+9Zvaamf06tm7vtNh1PwTMIuhS+QNgvLt/w90r2zuZiPRMlUmJr0qds8nnzj+C6eMSS4N87bGlLN5UEWFE4u4vxnpdjCQYBT4beCPWyfmDZpaWes9YJdbDwIMEjbPKgTlmNrKNY08A/gXMIWi49SrwlJkdHztkLjCw1e0EoBl4PHSqCcB/tjruo6l+b9Jz1DY08d9/WxbfPuOwoVx09CERRpRdkkqdNeIr0mGZSHyfBUYBNxGs27sQ2Au8AOwBNhB8EFZnIBYRiUhScyslvlmlID+PO687niGx+b71jc3c8uBr7KqqizgycfdKd7/T3Y8FzgBWAHcBW8zsZ2Y2NcWXvB2Y5e6/cPcVwI0E3xVubuPY/wAWu/vX3P1t4DZgLfCZWOxNsfjjt9h5VgKPhc4zEVja6tiaFL8v6UHumr2azRX7ASjIM26/dCpaujohqdRZI74iHZb2xNfd3+PuhwBjgCsI5u8+DxwO/DswFjgX2G1mc2If5B9Od1wiklkqdc5uhw7ow53XHU9seV+2VNby6T++TmOTlsrIFu7+krvfSPIo8BIzm2tmp3X3/GY2gaBD9COha9YCs4GZrY4tBC5vdawD/9f62NBrDiVIpL8fO7blPKOBVd2NX3qHt7bt5ZezE/85fPSMCRw2XN3mw9TVWaRrMtbcyt03u/tj7v5fsWR4BDAOuJIgGZ5HsF7vp4F7Un39Ts5ZGmtm95lZk5l9vY3nzzzAnKXits4nkuvcncr9am6V7U6bNJSvXJRYRWfe6l18//9WRBiRtMXd98ZGZI8DTiOYC3xqCk49Lna/ptX+DQTJadgoIP8Axw49QCn2bQSl038I7RsbO89sM1tvZo+a2UkdDdjMlrZ1A9QFqQdqbGrmi39dnNTQ6tNavugdSpO6OivxFemoKLo6DWk0cQAAIABJREFUx7n7Rnd/NFYmdZG7DwfGA1en8jqdnLM0lKCUrA8H/t9nArCMVnOX3F11gSJtqK5vin+RAc3xzWYfO3MCFx97aHz7N3PW8MjrmyKMKPeY2XAze5+ZXWxmpa2eO8rMrmn58dbd57v7x9z9Rym4dMskytYTvKuAQZ04FoLPxTgzGwzcAvzI3etDT20HpgFnArcCQwmS4DM7Hb30eL+du5bFmxKtXr5z5TH0LVZDq9ZKNOIr0iVZ99fE3TcQ/GKcSrcTm7MEYGY3AlsJ5hrd3ur65WY2w90Xmdl7D3C+icDbasQl0jF7quuTtgf00YhvtjIzvn/1sby9vYoV2/cB8KW/LmHs4DJOHNc695FUM7MZwN+BAYABO2Pr3j8H/A74YOjYtcAV7r4kRZcvj933b7V/ALCtE8c2h55v8WmgAbg7vNPdq4DFsc3lZvY0sJyg2dUFBwvY3dus3oqN+qZ6/rOk0aodVfz46ZXx7atPHM05k4dHGFH2Cje30nJGIh0X6YhvJnRmzlLo+UUHOe0ENB9JpMMq9yfm9/YrLqAwv9f/6enRSosKuPuD0xnc0uyqqZmPP/hqvNmMpNUPgDrgy8BVwFMEn19fJFgN4QngswTLHA0H5rZVvdRFLUP741vtHxd6rsVmwA9w7FZ3b2zZYWZ9CUZzfx5LdA8o1tTqVeDQ9o6T3qWp2fnSrMXUNwY9BYb1K+ZrM/W7xYGo1Fmka1Lx7dPM7EexphXZqDNzljpqIvBxM9tuZgvM7FMdfaHmI0kuCnd0HqD5vT3C2CGl/OqGEynMD7pdlVfV87H7X6W6rvEgr5RumgZ8091/4O6PxJY0+jnwHeAP7n5ZrMvzZ4CTCUaFv5aia6+I3S5r2WFmZcA5BMsPxrn7fuDpVsfmARe1Phb4OMH0oTsPFkCs2dU04PWuvAHpmX45exWvrd8T3/725Ufrs6IdSaXOSnxFOiwViW8JwQflM2b2q9gIazbpzJyljvoEcCLBr/FzgJ+Z2UE/0EVyVUWoo7OWMuo5Tp4wmG9fcUx8e/nWvfz7QwvV6Tm9thMkiWE/JUhwHwjvdPe3gF8SrIzQbbFOy98Brjezm8xsMnAvUAnca2bvNrNdZnZh7CXfA04zs6+Z2WHAzwjm6MY/D2NNHz8H/Nbdd7a+ppndYGYXmtloM5sG/BkYRqtpSNJ7vbZ+z/9n777j4yqv/I9/jnqxJbn3DrZxxcbGxoAhoQRCX0LvKYRsNpCQ7GZTSEjyS5bNbhppG0roIYQWQggQerNNMQY3bIx775bV6/n9MWPpjpBsydbMHc1836+XXtK9c+fOGVme5577PM95+MXzK5q2z548kFPHa83e/YlZzkhDnUXa7ZATX3evdffbgIlE1uZ93Mzui8PaggerI3OW2sXdF7v7h+7+urt/A/gF8GUzG9SO545v7QtYeTCxiHQFewI9vlrKqGu5cNoQrp09smn7peXbuemJJURXo5HO9zTw1WB7Ek0Yv02kqGJLi4kMee4U7n4vkfoX1wNvELm5Pcvd9xCZu1tDZK4u7v4ikWUKzwbeBo4ATnD34FSgq4kksm0V3xpMpA1dTmQYdy0wxd1Xd9Z7kuRVWlXH9Q8uoKEx8nkyqCSfH507IeSokp+GOoscnE6baOfuje7+QHR5hceAP5rZ4x1ZliBOOjJn6WC9SuRufL9OOp9ISoldw1c9vl3NN08by2mBHpgH31rH717Wvbo4+T5QASwxsx9Ee0Fx91vcfWMrxx9O6wnxQXP3O919orv3dvezo0Uncffn3X2gu78QOPZJd5/u7j3c/eSWNTLc/Q/unu3ua9t4rVvc/Qh3L3T3Ie5+sbu3nJokKcjd+c7ji5pqB2RmGLdecqSKH7aD1vEVOThxqTATnZc0k8gQrJ+a2T/N7MR4vFY7tHvO0iGYRuQueKdefIikitihzrqo6WoyM4xfXnwkU4c2r1DzP88u57F3tcxRZ3P33cB04Fkivbzvmtl2M3s6mgifaWb9AMxsIpGpN38ML2KRg/PQ2+v5+8LNTdtfO/lwjhrWM8SIuo48VXUWOShxLa3q7v90908QmavzDTN7I7qmbsIcxJyl/TKzsWZ2nZkdYWZDzewrRKpvfjtaLVpEWogd6qwe364oLzuTO66azojehU37/uORhby0bFuIUaUmdy9z94uIVDb+EjCfyBr0NwFPAJvMbB2R0UY1wBYzGx5OtCIdt2Ddbr73xJKm7Zkje/KlEw8LMaKuRUOdRQ5OQtYUcfc57n4m8GXgSjObb2YXmZkl6PXbPWepHUqAS4lccCyO/nyZu/+8s+MWSRXBqs4lGsbWZfUszOHua6bTK7rMUX2jc93983l7za6QI0tN7r7D3W9z99PcvQeRYc0XE1nyaDmR9msgkWR4pZlVmNlSM3s9vKhF9m9bWTVfuv9daqNF8noV5vCLi44kMyMhl4QpoeVQZ9VcEGmfrES+mLu/B1xkZocD3wJuMrOfA/cG1/yL02vfCdzZyv7niVw4tPacj30Ku/s8YHanByiSwvYE1vHtUajEtysb1quQu685mktun0d5TT019Y189u63eejaYxg3sGUNQelM7r6SSCHEh/fti/b0HgVMBcYDo2mjTRMJW219I19+4F227I0MkMvKMH532VQGFLcsZC77E0x83aGmvjFmiSMRaV1CenxbcvcV7v5ZIuv9TSFSxON6M8sLIx4Ria/gHN+SfA117uomDi7m9iunkZMVaULKquu58o9vsWp7eciRpR93X+Puj7r7d9z9XHcf5+4jD/xMkcRyd3749yW8vaZ5vd7vnnEEM0b2CjGqrik/JzbJ1Vq+Iu0TSuK7j7uvd/evAMcTuUO9yMy+ZWbdw4xLRDqXljNKPceM6sVvLpnSNDxxR3kNl97+Jmt3VoQcmYgkoztfX83989Y1bZ8/dTBXzRoeXkBdWMveXVV2FmmfUBPffdx9m7v/J5FKllnAAjP7sZn1Djk0ETlEjY1OaXCos4pbpYxTx/fnp+dPatresreaS29/kw27K0OMSkSSzVMLN/P/nvqgaXvS4GJ+fN4EElTqJeVkZ2aQndn8u1PiK9I+SZH47uPue9z9R8AkYCcwx8x+YWaDQg5NRA5SWU09jYG6G1qjMbWcf9RgfnLexKbtjXuquOT2eWwurQoxKhFJFm+v2cXX/vJe0/bgHvnccdU0zUk9RDEFrjTUWaRdkirxDegB9Il+vwF4P9xwRORgldfE1q3rlpfQmnqSAJfOGMoPzh7ftL1+VxUX3zaPjXuU/IqksxVby/j8Pe9QWx+p4Fycn83d1xxN3+4q6XKogvN81eMr0j5Jlfia2XFm9hdgNfBNIksM/YBIpUoR6YLKqpuHOedlZ5CdmVQfO9JJrpo1nO+ecUTT9tqdlVz0h7ms36VhzyLpaNX2ci69482mqS45WRnccdU0DuvbLeTIUoN6fEU6LvQrUDPLMbOrzWw+8ArwGWA+cDkwzN1/4O5bQw1SRA5aWXVzj2/3PA1zTmWfP34k3/702KbtDburuOgPc1mzQwWvRNLJup2VXHr7m2wvqwEgw+CXFx3J9OE9Q44sdeTnNI+eqlTiK9IuoSW+ZjbQzP4fsJ7I+rrjgQeAGe5+jLv/yd3r9nsSEUl6wR7f7hrmnPKunT2K7581rml7U2k1F/5hLh9uLQsxKhFJlE17qrj0jnlNa/UC/PQzk/n0xAEhRpV68rObL+GrNdRZpF1CSXzN7DdEhjN/G6gHbgaGuvuV7v52GDGJSHyoxzf9XHPsCH507oSm7W1lNVz4h7ksWLd7P88Ska5u3c5KLrptLht2N8/v//F5E/jMUYNDjCo1aY6vSMeF1eP7V2ALcDGR4cw/cvdtIcUiInEUk/jmqsc3XVwxcxg/PX8S0WV+2VNZx2V3vMnrK3aEG5iIxMWHW8v4zP/NYf2u5qT3e2eO47IZw0KMKnXlZ2uos0hHhZL4uvvzwK3AJnevP9DxItJ1xfb4KvFNJxdOH8JvLp3atN5kZW0Dn737bZ5auDnkyESkMy3aUMpFf5jLtuicXoDvnnEEnz1uRIhRpbZgj6+GOou0T2hzfN39Z0BRWK8vIomhOb7p7dMTB3DnVdObKpDWNjTybw++y52vrw45MhHpDK98uJ1Lbp/H7srIZ70Z/Pf5E/n88SNDjiy1Bef4qqqzSPuEWtXZ3f8R5uuLSPxpjq/MHt2H+z8/g+L8yL+/O/zo70v5f39fSmOjhxydiBysP725js/e/XbTeu1ZGcavL5nCRdOHhhxZ6itQVWeRDot74mtmBfF+DRFJXvsuiEA9vunsqGE9ePRLsxhUkt+0747XV/OVBxdomJ5IF9PY6PzX0x/w7ccX0RC9eVWQk8ntV03jzEkDQ44uPeRlq7iVSEclosd3i5n9zswmJ+C1RCTJBIc6d1Nxq7R2WN9uPP6vsxg/sHmWy1OLNnPRbfPYVla9n2eKSLIorarj2vve4Q+vrGra168ol7988Rg+MaZviJGll/xszfEV6ahEJL43AbOBd83sTTO7xszyD/QkEUkNewNDnYs01Dnt9S3K46EvHsPs0X2a9r2/fg/n/uYNlm7aG2JkInIgy7bs5ZzfvM7zHzQvxDG2f3f++uVjmTCoOMTI0k9BoLhVZa3qxIq0R9wTX3f/lbtPIJL8LgN+A2w2s1+b2cR4v76IhEtVnaWlbrlZ/PGqaVw+s3ke4KbSaj7zf3N4ZvGWECMTkbY88d5GzvvtHNbsrGzad9LYvjx83TEMKFZ/RqLlxazj2xhiJCJdR8KKW7n7G+5+FTCQSC/wCcB7ZjbHzK40s7xExSIiiRNb1Vk9vhKRlZnBj86ZwM1njWta67eytoHr7p/Pz5/7UEWvRJJEeU0933j4fW7483tNc0nN4MZTRnP7ldP0uR6SmKHOKm4l0i4Jr+rs7qXu/mt3nwQcBywHfgdsMrNfmdm4RMckIvGj4lbSFjPj6mNHcOfV0+kemP996wsruPa++TE3TUQk8d5fv4czb32NR+ZvaNpXlJfFH6+ezvUnHU7GvrtWknAxQ53rNNRZpD3CXs5orrtfQ2wv8CIze93MZoUZm4gcOnePGercTYmvtOITY/ry+JePZWSfwqZ9z3+wlXN+8wbLt5SFGJlIeqqtb+QXz33I+b+PHdo8ZWgJT11/vIpYJYFgj6/W8RVpn1AT333cfa+7/9bdjwRmEZkLfEzIYYnIIaqqa2ha6gLU4yttO6xvN/765WM5aWzzBfWqHRWc+9s3eHzBhv08U0Q60+KNpZz9m9f51QsrqI9+fpvBVz55GH/54jEM6alVKpNBXkxVZ83xFWmPUK5CzawIGAcMAtYCK919N4C7vwm8GUZcItK5gr29oKrOsn9FedncfuU0fvnCCm59YQUQuXnytYfe5+01u/nemeNiLvZEpPNU1tZz6wsfcftrq2JuWA4qyednF05m5sheIUYnLamqs0jHJTzxNbNfAtcB2S32vwrc4u7PJjomEYmP4BzN7EwjNyspBplIEsvIMG48ZTRThpbwtYfeY09l5G/oT2+u4921u/nNpVM5rG+3kKMUSR3uzrNLtvDDJ5eyqTR2Pe3LZgzlW58+QmuwJ6H8mKrOGuos0h5hXIV+HrgXuBG4CvgO8BfgaOAfZnZrCDGJSBzELmWUjZkKoUj7fGJMX578t+OYNLh5bdBlW8o469ev8/A763FX1WeRQ/Xh1jKuuuttrrv/3Zikd0jPfP70+Rn8+LyJSnqTVH6Loc6qhC9yYGF8mo1y960td0aHP/8U+LKZrXT3XyU+NJGOa2h0Vm4vZ+GGUtburGDD7io27K5kR3ktdQ2NNDQ69Y1OSX42g3rkM6gknxG9C5k9ug+H9+2W0slgTGErXTxJBw3pWcDD1x3DLU8v46431gCRno1/f2Qhr3y4nR+fO5HiAg2fF+mo7WU1/OL5D/nzW+sI5kvZmca1s0fyb584PKZHUZJPy3+f6voGCnLUzorsT8L/h7SW9Eb37wWuM7PewOWAEl9JSnUNjSxYt4fXV2xn3qpdLN5USmU7KipuL6thxbby5h1PfcCQnvmcNLYfnzlqMBMGFbf95C4qtsdXDbJ0XG5WJt8/azyzRvXm3x95v2no898Xbmb+2t387MLJzBrVO+QoRbqGvdV13Pnaau58fXXMUnMAxx3Wmx+cM55RfTSVoCvIb1HvoKpWia/IgYT2P8TMBgINrSTCS4FTQghJpE17Kmt54YNtPLtkC3NW7vzYBcPBWr+rirvnrOHuOWs4fUJ/vn7qaA7r271Tzp0MgnN8lfjKoThlXD+evuF4bvjze7y1ehcAm0urueyON/ncsSP4xqfGqPCVSBsqauq5Z+4a/vDKKkqrYtfHHtG7kG+eNpZPje+X0iOQUk3LzzvN8xU5sDCvRC8G/sfMNgLvAB8C/YCLgEdCjEsEgNKqOp5ZvJm/L9zM3JU7m5Z1aE2Gweh+3RnbvzuDexQwuEc+/YrzyM3KICsjgwyDHeW1bNhdyYbdVbyzdheLN+6NOcfTi7fw7JItnD91MN89cxzF+V1/CGfwBkF3VXSWQzSgOJ8HvzCT215dxc+fW05dg+MOd7y+mpeWb+N/L5jMlKE9wg5TJGmUVtZxz9w13PXGanZXxia8JQXZ3HDS4Vw2Yxg5KjzY5WRmRApG1tRHljLSWr4iBxZm4vsMUAxMBqYA5wJOpPDVtSHGJWmsrqGRl5Zt47F3N/Li8m3U1re+Nl5OVgbTh/fguMP6MH14D8YPLO7wfKgtpdW8uGwb989by9LNkSS40eHh+RuYv3Y3t181rcsPOduroc7SyTIzjC+dOIrjD+/NDX9ewMrtFQCs3F7B+b+fw3UnjOL6kw5X76+ktU17qrhnzhrun7eWihYJUbfcLD533Ag+d/wILTHXxeXnZDYnvurxFTmgMJYzynD3RndfCnw/sP8w4MvANcCDwD8THZukrxVby3h4/gYee3cDO8prWz2md7dcTh3fj1PH9WPmyF6HfGHdvziPS2cM5eLpQ3hq0WZ+/tyHrN4RuYhftaOCc3/7Br+5dConjO5zSK8TppihzipuJZ1owqBinrr+eP732eXc+cZq3CM3jn738kqeXbKF/z5/EtOG9ww7TJGEWrBuN398Yw3/WLQ5Zi1eiKz7esUxw7hu9ih6FOaEFKF0pvzsTPYQaWfV4ytyYGFciT5uZpe4e2Vwp7t/BHzNIhNMfgmMCyE2SSPVdQ38Y9Fm/vTmOt5Zu7vVY3oW5nDmpAGcNXkgU4f2IDOj8+c/ZWQYZ00eyOkT+nPrix9x6wsrgEhhqGvueosfnjOBy2cO6/TXTYSWyxmJdKa87Ey+e+Y4Th3fn39/5H3W7ow0Kyu3V3DBH+Zyxcxh/MdpY1VRvIPM7CwiSw2OITIV6avuvqSNY2cC/01k9NYi4Lvu/krg8eOBV1t5ap6710SP6Qf8FjgB2Az8n7v/rvPeUWqrqm3gyfc38cCba3l/Q+nHHi/Oz+aaY4dz1THDlfCmmOBIs0r1+IocUBhXA0cBy83sP4GH3b1l91oWMDjxYUm6WL2jggfmreXh+Rs+VuQDIsOYTx3Xj/OnDua4w3uTnZmYuU9ZmRnceMpoxvbvztf/8j5VdQ00Otz0xGIGluTxybH9EhJHZ1JxK0mEo0f05Okbjud/nl3O3XPW4A7ucO/ctfxzyVZuPnscnxrfX4V72sHMzgAeA74KPA/cDLxqZhPdfVOLY6cCLwI/Az4P/CvwjJnNcvcF0cNGEClaOSv43EDSWwDMiR5zPDATuMPMGt39/+LyJlOAu7Nk014emb+BR9/dEHOTcZ9BJflcc+xwLj56qG7+pKiYtXzV4ytyQGF8Eh4HPADcB9xqZs8BHwHlwNHA2cCzIcQlKayh0Xlx2TbunbuG11bsaPWYcQOKuGj6EM45ciAlBeHdFf/0xAEM61XA5+95h82l1bjD9Q++x2P/OovR/bpWxWcVt5JEKcjJ4vtnjeesyQP55iMLm5YO27K3muvuf5eTxvbl5rPHM6RnQciRJr2bgUfd/bcAZnYNkV7Ya6OPBX0bWOjuN0WPvRH4FHADcHX0mJHACnf/eFdkxGeBQcBR7r4HWBbtRf4PM7vN3VsvtJCmtpVV8+T7m3n4nfUs21LW6jHThvXgc8eN4JRx/chK0I1bCUcw8dUcX5EDC2Md3zXAsWZ2KZGG9AIgeBt+LnEobtXBoVtDgR8CVwA/dPcftHLMt4ArgR7Ac8C/7adhl5DsrqjlL++s5755a9mwu+pjj+dnZ3L25IFcOmMokwYXJ02P0PiBxdx1zXTO/90cKmobKK+p53P3vM0TXz6Onl1oqJrW8ZVEmzq0B3+//jh+//JKfvfSSmobInnTC8u28cbKHfzriYdx7eyRKn7VCjMbAUwD/nffPnevNrOXgTMIJL5mlk2kKOV3Ase6mT0LXB447QgiN7fbcgHwUjTp3edp4IvAEUCr7XQ6Ka2q49klW3jy/U288dEOWltgoDAnk/OmDuLSo4cxbmBR4oOUUMQMdVaPr8gBhXYl6u5/Av5kZsVE7gj3BlZH5/p2qg4O3eoNLAf+BrR6q9TMbgG+QGTppT3AH4lUqT6ms2OXg/PB5r3cM2cNjy/Y2FTxMOiwvt248phhnDtlUNJWtRzbv4hfXTyFL9z3Du6RNX+vu38+D3x+RsKGXx+qYOLbTYmvJEhuViZfPXk0Z00eyE1/XcyclTsBqK5r5OfPfcgj8zfwvTPHcfK4rjd9IM72FRNY1WL/OiJDkIMGAZltHNvbzPLcvZpI+/4vZnYFsBa4Z19vcuA1n2rlHBCZ9pSWie/uilqeW7qVfyzezBsf7aCuofXl9KYOLeGCaUM4a/JADWdOQzFDndXjK3JAoX9KRntJFxzwwENzM+0cuuXuO8xspru/b2YXtjyRmfUBrgdudPfno/uuBeaa2Ynu/nI834i0rb6hkec/2Mbdc1Yzb9Wujz2emWGcOq4fVx4znJkjeyZN7+7+nDyuH988bSy3PL0MgLdW7+KuN1Zz7exRIUfWPsE5vkVKfCXBRvXpxgOfn8Ff39vIj59axo7yGgDW7ark8/e+w+zRfbjpjCM4vItNIYij/tHve1rsLycyuqm9xwKUAFuALwG1QF8iPcS/MrOx7v6VwHnaOscBF2U2s7YS467xIRmwans5L3ywjec+2Mo7a3a12rMLkbm7Z00eyGeOGsxhfbv2kndyaII9vqrqLHJgKX8l2pGhW4HH39/PKc8E8oHHA/veBHZFz/fyocYsHbOnspaH3l7PvXPXsnHPx4cz9yrM4ZKjh3LZzKEMKM4PIcJD88XZI/lg816eeC8yOOHWFz7i3CmD6Ns9L+TIDmyvqjpLyMyM86YM5qQj+vGr51dw95w1Tcu8vPrhdk77aAeXzRjKV08e3aWmEcTJvgIILcfKFhNJYtt7bOO+x919cXT/h8DrZubA183sFnffGD2utXPQymumlKraBuat2snLy7fxyofbWbOzss1jexbm8OmJ/TnnyEEcNbQHGXFYYUC6ngJVdRbpkJRPfOnY0K32nq/S3bfu2xGd17SBdlSjTqW702FbvqWMu+es4fEFG6iu+/hw5omDirl61nDOnDyA3KyuO5/PzLj5rPG8vHw7pVV1lNfU8z/PLOd/Lpgcdmj7VVPfQG1gmLnm+EqYivKyuenMcVw0fQg3/21J0/Dnhkbn3rlreXzBRr78icO4etbwdJ7/uyH6fTgwP7B/WOCxfTYCHj2WFsdudvePlxmOeBX4BtAveo4NbZwjGE+b3H18a/ujbW1SLYtY39DIoo2lzFm5k9dX7GD+2t1Nc9Bb07d7LqdN6M9pE/pz9PCeKlQlHxP8rFKPr8iBpcOVaEeGbrX3fC3PdSjnkw5oaHSe/2Ar98xZ03ThGpSVYZw2oT/XHDuCqUNLusRw5vboUZjDjaeM5vt/i9w3eXj+Bi6fOYzJQ0pCjqxt5S2W11CPrySD0f2688DnZ/Dc0q38+B8fNK39W1Zdzy1PL+PeOWv4+qljOHfKoLis253klke/zgEeBTCzQuBE4HvBA929KroqwzlEljPCzDKA0/c9tw3TgBoiyxcBPAF818y6ufu+Ic5nAoviUfMjkWrqG1i0oZS31uxi3qpdzF+zi4oDJCdj+3fnlHH9OOmIfkwaVKyeXdkvzfEV6Zh0SHw7MnSrvedrrWRiMbDiQE/uSnenk8nuiloeemc997UxnLlnYQ6XzRjKZTOG0b84+YcAH4zLZgzlT2+uY/nWyBIWNz+5hEevm5W0F0bBwlZmUJC+vWiSZMyMU8f354Qxfbhnzhp+/eJHTX+vm0qr+frD73Pbq6v4+qmjOWVcv5S5gXYg0dFLPwHuMrPXiPTO/ggoje47FXgQuMzdnwH+G/inmd0U3X8DkUKVvwYws7FEkuZXgAoiSfJ/At+KFr4CuJ1I4cl7o6slHEukKnSwMnSXsHVvNQvW7WHB+t3MX7ObhRtK99ujC9AtN4vjDuvNCWP6cMLoPgws6XrTcSQ8BarqLNIh6ZD4dmToVnvP183Mern7TgAzyyQyzPmvhxCntGLxxlLum7uWv77XenXm8QOLuObYEZw5aUDKD0/Myszg+2eN49I73gRgwbo9/PW9jfzL1AOOsA9FTEXn3KykTdAlfeVmZXLt7FFccNQQfvPSR9w7d01T9dzlW8u49r75TBlawr9/agyzRvUON9gEcfd7o0sVfRX4L2AOMMvd95hZI5He2rrosS+a2XlEeoNvJNLGnhDoqS0BLiWSPOcSqdB8mbs/Eni9HWZ2PPBbIvUyVgBXuvtD8X+3B29PZS2LNpaycEMpizaUsnDDHjaVVh/wedmZxpQhPTju8N4ce1hvJg0u7jJV+iX55GkdX5EOSYfEt91Dt9rpKaAher4/RvcdR6TH97FDjFWA2vpGnl68mfvmruWdtbs/9nhmdDjz1bOxzTZbAAAgAElEQVSGM21Yj7TpjQGYdVhvTp/Qn6cXRwYr/PrFjzj3yEFJmVTGVnTWMGdJXj0Kc7jpzHFcdcxw/vefy/nb+82r3C1Yt4dLb3+TO66cljbLH7n7ncCdrex/HhjYYt+TwJNtnGceMLsdr/chcMpBBRtnjY3Oht1VfLBlLx9s3suSTXtZumlvqyOPWpOXncHkwSXMGNmLmSN6MmVoj5hKvCKHIvi3pKHOIgeW8onvQQzdOtD51pvZPcAtZraKyNze3wP/cPd34/ZG0sDGPVU8+OY6/vz2OnaU137s8V6FOVya4sOZ2+M/Tx/LM0u24A6rd1QwZ+VOjjs8+XqjymqCFZ1T/qNGUsDQXgXceskUrjthFP/7z+W8uGwbAMN7FXDCmD4hRyeJ9PeFm/jj66tZvqXsgPNygwaV5HPkkBKmDC1h2vCejB9YpB5diRsNdRbpmLS4Gu3I0K12+iKRxPleIlUtnyIyzEs6qLHReXXFdu6ft44Xl21tdd3CI4eUcNWsYXx6YteuztxZhvUq5BNj+jZdlN83b01yJr7VSnylaxo3sIg/Xj2dd9bs4qfPLufymcOUvKSZypoG3l3XWh3LZj0Kspk0uIRJg4uZMKiYI4eU0K8ofW/KSuLlq6qzSIekzdVoR4ZuBR5rdfxodJmGG1Gye9C2l9Xw8Pz1/Pmt9azb9fG1C3OzMjh78kCuOGYYkwYnb+XisFwxc1hT4vvc0q1sLq1KujWKg0Odu+WmzUeNpJBpw3vy0LUHs+qddHVjB3SP2R7Wq4Ax/bozbmAR4wcWM35gEQOK89Jqqo0knzxVdRbpEF2NSsI0NjpzVu7kwbfW8c+lW5qKyAQN7VnA5TOHcsFRQ+hRmBNClF3D7NF9GNIzn/W7qmh0ePCt9dx4yuiww4oR2+OrOb7SNSmxSU+j+3Xnx+dNYGz/Isb0766bd5KUCnKa/y411FnkwPRJLnG3bW81D8/fwENvt967m2Fw0hH9uHzmMI4/rHdSFmpKNpkZxmUzhnHL08sAePCtdXzlk4cl1XDMcs3xFZEuKi87k8tmDAs7DJH9yldV51BV1zVQWlXH3qo6KmsbqKxtoLqugdqGRuoaGqlvcBoanWA3T4ZFruHMjOwMIzszg+ysDHIyM8jJyiA3K4O87AzysjPJz84kPyeTvKxMXRt3El2NSlzUNTTy0rJt/OWd9by0fDsNrUze7V+Ux0XTh3DR9CFau/AgXDhtCD9/7kNq6xvZXlbDP5ds5YxJA8IOq0lwqLN6fEVERDpXfk7zzW4lvp2rpr6B9bsqWb2jkjU7Kli/u5ItpdVsLath295qdlXUtrrMZrzkZWdQkJNFQU4mhTlZFORGv+dk0i03i8LcLLrlZdEtN6t5OzeL7nnN20V5kWPyszPTdjSTEl/pVCu2lvHI/A08+u5GdpTXfOzxDIMTRvfh0hnD+MSYPmQlUQ9lV9OzMIczJg7g8QUbgUiRq2RKfPequJWIiEjc5AeGOtfWN9LQ6GSqZ7DD3J3lW8t4e/UuFm0sZfHGvXy4tYz61iquhqS6rpHqulp2VRz6uTKMaFKcHUmU85qT5O5NyXM23aLb3Vsk1ft+LszJ6nI90boalUNWWlnHkws38fD8Dby/vvUqmINK8rlw2hAumDZYvbud6PKZw5oS33mrdrFmRwXDexeGHFWEqjqLiIjET3CoM0R6fTUfvX1KK+v459ItvLpiB3NX7mh1Gc32ysnKoDAnk4KcLPKyM8jJyiQ708jKMLIyoh08Bjg0utPoToNDfXRIdF2DU1vfSE19IzX1DdTUNVLbEL/e5EaPdE4EOygOVmFOZpu9zYW50cdysijIzaJbdLt/UR7ThvfshHfScfrfIQelrqGRV5Zv57EFG3h+6bZW/4PmZGZw6vh+XDR9CLNG9dZdyDiYOrSEw/t2Y8W2cgBeXLaNzx43IuSoImKHOuujRkREpDN9LPGtVeK7P1W1DTy9eDN/X7iZ11Zsb7XIalBmhjGkRz7DexcyvFch/Yvz6F+UR9+iXHp3y6UkP5ui/OyY6tqdpaHRqa5roKqugarayPfK2gYqa+qprG2goraeipoGKmrqKa+pb/oe+3MD5TV1lFXXU15dH5ce7IraBipqG9hW9vFRnm2ZMaInD33xmE6PpT30v0Pazd15b/0e/rpgI08u3Myuitbvjk0YVMRnpg7mnCMHqTJznJkZJx3RrynxfWl58iS+5cEe31zN8RUREelMedmx08W0lm/rNpdWce/ctTz41jr2VNa1eowZjB9YxNShPZgwqJgJA4s5rG83crLCmZKXmWEURntPO4O7U1PfSHlNPWXVkeS4rLqesuq6poS5rLqevdV1kcS5OrJdHkio923XHuLc5jBvzijxlQNatb2cv763ib+9t5E1Oz9elRmgV2EO5xw5iAumDeaIAUUJjjC9fWJMH/7vlZUAvLlqF5W19TFLHIRFQ51FRETix8zIz85sKmylAlex1u+q5BfPfcgT729qtchqr8IcPjWhP7MP783Mkb0oKUjdzhozIy87k7zsTHp3yz2kc9XUNzT1Nu9Lhlvrca6I9jpX1jbvq6xtYGSf8Kbk6WpUWrVxTxVPLdzEk+9vZtHG0laPycnK4OQj+nL+1MHMHt0nqZbSSSdTh/Wge14WZdX11DY0MuejnZw8rl/YYamqs4iISJzl5yjxbam0qo7fvfwRd72x5mO9kwU5mZw9eSBnTR7IjBE9VWT1IORmZZKblUnPLjiqU4mvNNm6t5qnF0XmPryzdnebx80Y0ZPzpgzi9IkDKM5XQhO27MwMZo/uw1MLNwPw4vJtoSe+DY1ORWDIlXp8RUREOl9wnm9l7aEXK+rqnnhvIzf/bQm7WwxpHtwjn6tnDeeCaUN07ZrGdDWa5rburebZJVv4+8LNvL1mF97GvPex/btz7pRBnDV5IINUlTnpfGJM36bE9+Vl23D3UNdoK29RKVCJr4iISOcryGlOfNN5jm9pZR3ffWIxT76/KWZ//6I8vn7qaP5l6mAVWRUlvulow+5Knlm8hWcWb2H+ut1tJrvDexVw1uSBnDlpIGP6d09skNIhJ4zu0/TzptJqPtxaHuq/WVlN7J1WVZkUERHpfAWB9rUiTRPfuSt3cuNf3mNzaXXTvsKcTL504ig+d9xI8nM6v+qydE26Gk0D7s5H28p5dskWnl2ytc05uxAZCnLGxAGcMWkAEwcVh9prKO3Xp3sukwYXs3BD5N/2peXbwk18Az2+BTmZmkMjIiISB4WBpK6yJv2GOv/lnfV8+7FFMUv1HD28Jz+7cDJDehaEGJkkIyW+Kaqh0Xl33W6eW7qV55ZuZfWOijaPHdIzn9MnDOCMiQOYNFjJbld14pi+zYnvsm1cd8Ko0GJRRWcREZH4C67iUJlGPb7uzi+fX8GvXljRtC8707jxlDFcO3ukhjVLq3RFmkLKa+p57cPtPP/BNl5avq3NdXYBRvYu5FMT+nPGxAGMH1ikZDcFfGJMH26NNgDvrN1NaVVdaAUcghWdNcxZREQkPoJzfNOluFV9QyPfemwRD8/f0LSvd7dc7rxqGpOHlIQYmSQ7XZGmiJeWb+OL986ntqHtRaUnDCriU+P6c9qE/hzWt5uS3RQzeXAJvQpz2FlRS0Oj8/qKHZwxaUAoscT2+Kp6ooiISDwU5jYnvukwx9fd+eaji3j03eakd1SfQu6+5mgNbZYDUuKbIsYPKPpY0puVYcwY2ZNTx/Xn5HH9VI05xWVkGCeM7sNjCzYC8PpHISa+NRrqLCIiEm8xQ53TYI7vLU8vi0l6jx7ek9uuPIqSgq63pqwknq5IU0TfojwmDylh9fZyPjm2Lycd0Y8TxvShSL1taWX6iJ5Nie/STW0XMYu3cs3xFRERibtgcatU7/G97dWV/OHVVU3bM0b05J7PHk1etqo2S/voijSF/P6yqfTpnku2KuimrfEDi5p+XraljPqGxlAqKlcE7joX5uhjRkREJB6Cyxml8jq+j727gZ/8Y1nT9hEDirj9qmlKeqVDlCGlkIEl+Up609zoft2bKhnW1Deycnvb1bzjqTyY+Kq4lYiISFwUxPT4puZQ52Vb9vKtxxY1bQ/tWcA9n52uUY3SYcqSRFJIXnYmh/ft1rS9dHM4w52DPb6q6iwiIhIfsXN8U6/Ht7K2nn/70wJq6iN1bHoV5nDf546mb/e8kCOTrkiJr0iKGTegebjzko17Q4lBPb4iIiLxV5jiPb4/fHIpH20rb9r++UVHMqxXYYgRSVemxFckxYwLzPNdsin8xLebiluJiIjERXCOb2WKzfH92/ub+PPb65u2v3jCSE4Y3SfEiKSrU+IrkmLGDyxu+nnJplLcPeExxA51VuEJERGReAjO8a1MoR7fTXuq+E5gXu/kISV849QxIUYkqUCJr0iKCfb47q2uZ8PuqoTHUBGYZ6SqziIiIvERk/im0BzfHz/1AWXRm+jdc7P49cVTVMBVDpn+gkRSTHF+NkN65jdthzHcuVzFrUREROIueHO5orY+lFFene2Nj3bw1KLNTds3nTmOob0KQoxIUoUSX5EUFCxwtXRz4hPfYIENFbcSERGJj4LAdKJGp6n6cVdV19DI9/+2pGl7ytASPnPU4BAjklSixFckBQXn+S7dlPgljSpU1VlERCTuWk4n6uoFru6Zs6apirMZ/PDsCWRkWMhRSapQ4iuSgsaHWNm5pr6BuobmoVYa6iwiIhIf+dmxBSSDN567mm17q/nl8yuati85eigTBxfv5xkiHaPEVyQFBXt8N5dWs6uiNmGvXdGiuEahqjqLiIjERUaGxSS/XbnH93//ubypRkhxfjb/rirO0smU+IqkoH5FufQqzGnaXpLA4c4t7zarqrOIiEj8BG8wV3TRJY3W76rksXc3Nm1/49TR9Ahcx4h0BiW+IinIzGKWNVqawOHOwYrOBTmZmpsjIiISRwWBG8xddUmjP7y6kvrGyDSpwT3yufjooSFHJKlIia9IihoX0jzfCi1lJNKlmdlZZjbPzHab2XNmNn4/x840s1fMbI+ZvWZmJwQeyzaz68xsbvTx1Wb2SzMraHGOF8zMW3z9Vzzfo0gqiVnLtwv2+G7bW81f3tnQtH3dCaO0Zq/ERdr8VXVWQx59/PhWGmk3s9z4vxOR9gnO803kUOcyJb4iXZaZnQE8BtwHzAR2AK+a2cBWjp0KvAi8CkwH3gGeMbMp0UNygEuBnwKTgK8CVwC/a3GqEcCVQEng6wed+sZEUlhs4tv1enxvf20VtdFlmPoV5Wr5IombtEh8O7khh0gjvZTYRrrE3Wvi+T5EOiK4lu+qHRXUNyRmbT8tZSTSpd0MPOruv3X35cA1RK4Vrm3l2G8DC939JndfAdwIrAZuAHD3Cnef7e6Pu/s6d38C+C1wmZnlAJhZFjAEWOTupYGv6ni/UZFUEWxru9oc310Vtdw/b13T9heOH0letopiSnykReJLJzbkUSOBFS0a6cQvliqyH/2KmgcguMPe6sQ0hrGJrxovka7CzEYA04DH9+2LJqAvA2e0ODYbOLfFsQ482/LYFnYAWcC+4c5DotsrD/kNiKSpmB7fLjbH9643VlNVF4m5Z2EOl87Q3F6Jn5RPfOPUkI8APopb0CKdoFtuFpmBwlJ7KhOzpFF5oNHVUGeRLmVY9PuqFvvXAS3HHg4CMts4treZ5bXxGscCq919T3R7ZPT7SjNbZWYPmNno9gZsZkta+wJGtfccIl1dcPWErtTjW15Tz91z1jRtf+64ETGFukQ6Wzr8de2vIZ/ZYt8BG/Jo0jwS+BczuwJYC9zj7r9tTzDRBrk1aqSlU5kZJfnZ7Iyu4VtaVZeQ19VQZ5Euq3/0+54W+8uBHh04FiJTgLYEHzCzo4DPAF8J7H4PmBj9eQLwHeB1M5vt7ss6FL1ImsoP9PhWdaE5vk++v4my6Gi07rlZXHHMsAM8Q+TQpMNVaTwa8i8BtUBfIj3EvzKzse7+FUSSSHEg8d2jxFdE9m9H9HtRi/3FtEhiD3BsY+BxAMxsCJHRVM8Ct+3b7+47gZ3RzcVm9hKRYc9fA754oIDdvdVCldGbzOMO9HyRVNBV5/j++a3mub3/MnUQRXnZIUYj6SAdrko7vSF398XR/R8SuTPtwNfN7BZ338h+qJGWRCouaG5ESisTk/iWq6qzSFe1bz2R4cD8wP5hgcf22Qh49FhaHLvZ3Zs+CMysF5GEdxNwQfCxltx9q5ktBwYcRPwiaakrzvFdumkv729oLo+jdXslEVJ+ji+xDXnQITXkLbwKGNDvoKMUiYOS/ObEN1FzfGN6fDVXR6QrWR79OmffDjMrBE4EHg0e6O5VwHMtjs0ATg8eG33+U0SuN85094r9BWBmxcDhwIJDeysi6aMrzvH989vNvb2Th5RwxICWfU4inS8dEt9Ob8hbMQ2oIbLEkUjSKA4kvqVViWkMg8WtVNVZpOuIFnP8CZHlhr5gZmOAu4BS4C4zO9XMdprZadGn/Dcwy8xuMrPDgF8BvYFfA5hZAfA3IvUzLgDqzKw4+pUdPeZfzexEMxtoZjOjx+8GfpmwNy7SxeV3sXV8q2obeHxB8wDJS6YPCTEaSScpn/jGoSEfa2bXmdkRZjbUzL4C/Cfwba07KMmmpCCn6ec9VYnv8dVQZ5Guxd3vJbLU3/XAG0AeMCtahbmRyE3euuixLwLnAWcDbwNHACe4+75VD24APkmkIvRCIvUz9n1dGj1mNHAHkXm9fyIyhehod98d1zcqkkKCN5m7QuL7j0Wbm4paFeZkctbkgSFHJOkiLa5K3f3e6N3lrwL/Bcwh2pCb2ccacjM7D/gekTV85xPbkJcQabB/BOQCS4DL3P2RRL4nkfaI6fFN0Bzf4DCrbnlp8REjklLc/U7gzlb2Pw8MbLHvSeDJNs7zX0Ta3P291leJtM0icpCCSwAFbz4nq+Aw57OPHKRCmJIwafOX1okN+TxgdjxiFOlsJcHiVgmq6lyuqs4iIiIJE5zjm+w9viu2lvH2muYBHZccrWHOkjgpP9RZJJ0Fe3wTtZxRebWGOouIiCRKQcxQ5+Tu8X3k3ea6suMGFDFxUHGI0Ui6UeIrksKCPb6q6iwiIpJ6CrpIcSt35+lFzSuJXjBtMGYWYkSSbpT4iqSw4vzm4laJGOrc2OhUBBpd9fiKiIjEV8uhzo2NHmI0bVu6eS/rdlU2bZ8+Qct1S2Ip8RVJYS3n+EaKnMdPZV3snWYtZyQiIhJfwR5fgKq65Oz1fWZxc2/v1KEl9C/OCzEaSUdKfEVSWHCOb12Dx30IVMtqkipuJSIiEl8t29qKJJ3n+49Fm5t+Vm+vhEGJr0gKCya+EP8CV8GKzlkZRm6WPmJERETiKTcrg+BU2cqa5OvxXbG1jJXbK5q2T5vQP8RoJF3pqlQkhWVnZsTMs413gauKFksZqWiFiIhIfJlZ0i9p9HRgmPPEQcUM6VkQYjSSrpT4iqS4YK9vvAtcBXt8VdhKREQkMWIrOyffUOfgMGf19kpYlPiKpLiYxLcyvolvRWB4lQpbiYiIJEZwnm9FkvX4rtlRwbItZU3bpyvxlZAo8RVJcTFr+ca5x7dCPb4iIiIJF9PjW5NcPb7BYc5j+3dnZJ9uIUYj6UyJr0iKi0l849zjW95ijq+IiIjEXzDxTbYe36cXa5izJAclviIprjg/p+nneM/xVY+viIhI4hUEiltVJdEc3+1lNSzcUNq0rcRXwqTEVyTFxRa3im9VZ/X4ioiIJF6wrkYy9fjOWbmj6ecBxXmM6dc9xGgk3SnxFUlxYQ11Vo+viIhIYgR7fJNpju+rHzYnvscd1lvLHEqolPiKpLiS/MQlvrHr+Kqqs4iISCIk4xxfd+f1j7Y3bR8/uk+I0Ygo8RVJecEe3/jP8Q0uZ6QeXxERkUSI6fFNksT3o23lbN1b07R97KheIUYjosRXJOUV5Scu8dVQZxERkcQrDC5nlCTFrV5d0TzMefzAInp1yw0xGhElviIpryRQ1XlPZXyLW8UMdc5R4isiIpIIBYGbzcHRV2F6fUVgmPPhGuYs4VPiK5LigkOdK2obqGtojNtrqaqziIhI4iVbj29NfQPzVu1q2j7+8N4hRiMSocRXJMUFE1+I73DniloNdRYREUm0/CQrbvXu2j1U1UXiyMvO4KhhPUKOSESJr0jKy8/OJDuzefmAeFZ2Dg6v6panxFdERCQRgtOLqpKgxzdYzfnoEb3Iy9ZKDxI+Jb4iKc7MKA7M8y2tit8839jiVmrkREREEqEg0OYmwxzf1wKFrY4/TMOcJTko8RVJA8HhzvHq8a1raKS2vnn+sOb4ioiIJEZhzHJG4fb47q6oZdHG0qbt40cr8ZXkoMRXJA2UJGBJo2BFZ1DiKyIikigFSTTHd+6qnbhHfu7TPZcx/bqHGo/IPkp8RdJAcX78e3zLqlskvlrOSEREJCGCyxnV1jdSH8cVHA7k7TXN1ZyPGdkLM9vP0SKJo8RXJA0UB4c6x6vHNzC0Kj87k8wMNXQiIiKJEFzOCKCyLrxe3/lrdzf9PG24qjlL8lDiK5IGSoLFrSrjU9yqQmv4ioiIhKKgxSirypAKXFXW1rNk096mbS1jJMlEia9IGggWt4rXHN/y4FJGqugsIiKSMDlZGTFLF1aEVODq/fWlNDRGJvgW5mQytn9RKHGItEaJr0gaiJnjm4DiVurxFRERSaz8wFq5YfX4zl/bPL93ytAemvYkSUWJr0gaSMRyRuVKfEVEREITbHvDWtIoOL9Xw5wl2SjxFUkDxQlezqi7El8REZGECi5pVBnCkkaNja7CVpLUlPiKpIGSgkBxKw11FhERSTnBtjeMOb4fbS9nb3RpwwyDI4eUJDwGkf1R4iuSBmLX8a2lMVp4ojMFi1sp8RUREUmssOf4Bnt7x/Qvonte9n6OFkm8tEl8zewsM5tnZrvN7DkzG7+fY2ea2StmtsfMXjOzE1o55ltm9oGZbTGz+8ysOL7vQOTglQQS30aH8jjcCQ72+Kqqs0jXlej20sz6mdkjZrbdzBaa2b/G432JpLqwe3zfWRMY5qz5vZKE0iLxNbMzgMeA+4CZwA7gVTMb2MqxU4EXgVeB6cA7wDNmNiVwzC3AN4CvAGcCk4Fn4vw2RA5aUX7sXdfSOBS40lBnka4v0e2lmRUAc4Bc4Hjg58CtZnZdPN6fSCoLe45vsKKz5vdKMkqLxBe4GXjU3X/r7suBa4i892tbOfbbwEJ3v8ndVwA3AquBGwDMrA9wPfAdd3/e3d+JnmemmZ0Y93cichAyM4yivOZkNB7zfMtjenyV+Ip0UTeT2Pbys8Ag4Ap3X+budwN3AP9hZulyjSLSKQpzwqvqvKO8hjU7K5u2pw5V4ivJJ+UbFTMbAUwDHt+3z92rgZeBM1ocmw2c2+JYB54NHHsmkB88BngT2NXyfCLJpDjOSxppOSORri2k9vIC4CV33xM45mlgBHDEob4nkXRSEJhmVJHgOb7B+b39inIZ3CM/oa8v0h7pcHU6LPp9VYv964gM4woaBGS2cWxvM8uLnq/S3bfue9Dd3cw2AIM7LWqRTlaSn8N6qgD42l/eozCnc+fhbiqtbvpZia9IlxRGezkMeKqVcxA9ZklH34RIugoOdX50/gZeXr4tYa8dHEk2bVhPzCxhry3SXulwddo/+n1Pi/3lQMtxGPs7FqAkekzLx9s638eYWVuN+KgDPVfkUJQEeny3l9WwPY6v1T0vHT5aRFJOGO1la8fsO4faVJEO6Jbb3M6X1dRTVpP4AlcAU1XYSpJUyg91JlKYA6Coxf5iYEsHjm2MPr6jlcfbOp9I0jhvyqCEvE7/ojyOGdkrIa8lIp0qjPaytWP2VX1WmyrSAadP6E9uVriX9kV5WZw1aUCoMYi0JR26ZTZEvw8H5gf2Dws8ts9GwKPH0uLYze5eHx2i1c3Mern7TgAzyyQyJOuvBwrG3VtdFiJ613rcgZ4vcrD+Zepgpg3rydpdFXF7jcwMY+rQHuRlazkjkS4ojPZyQxvnCMbTJrWpIs2G9y7krW+fzMKNrQ20iL8MMyYNLtb6vZK00iHxXR79Ogd4FMDMCoETge8FD3T3KjN7Lnrsz6LHZgCn73sukblIDdFj/hjddxyRO9SPxfF9iByyob0KGNqrIOwwRCQ5hdFePgF818y6ufu+Ic5nAovc/aNOfn8iKa+4IJvjD+8TdhgiSSnlhzpHq0z+BLjMzL5gZmOAu4BS4C4zO9XMdprZadGn/Dcwy8xuMrPDgF8BvYFfR8+3HrgHuMXMTjSzacDvgX+4+7uJfXciIiKdI6T28nYic3rvNbMxZvZZ4HLgxwl50yIikjZSPvEFcPd7iawdeD3wBpAHzIoun9AI1AB10WNfBM4DzgbeJrKcwgkt7jx/EbgfuJfIne2XgfMT8V5ERETiJdHtpbvvAI4HuhNZ6uhLwJXu/lDc3qSIiKQli9zglbCZ2ZJx48aNW7JEKzeIiKSK8ePHs3Tp0qVtzUWV+FCbKiKSeg61TU2LHl8RERERERFJX0p8RUREREREJKUp8RUREREREZGUpsRXREREREREUpoSXxEREREREUlpSnxFREREREQkpSnxFRERERERkZSmxFdERERERERSmrl72DEIYGZ7c3Nzu48aNSrsUEREpJOsXLmSmpqaMncvCjuWdKI2VUQk9Rxqm6rEN0mY2RagAFgfdiydZN/VxspQo0ge+n3E0u8jln4fsVLp9zEEqHT3/mEHkk5SrE1Npf8PnUW/k1j6fcTS7yNWKv0+DqlNVeIrcWFmSwDcfXzYsSQD/T5i6fcRS7+PWPp9iDTT/4eP0+8kln4fsfT7iKXfRzPN8RUREREREZGUpsRXREREREREUpoSXxEREREREUlpSnxFREREREQkpSnxFRERERERkZSmqs4iIiIiIiKS0tTjKyIiIiIiIilNiYZ2kCkAAAg0SURBVK+IiIiIiIikNCW+IiIiIiIiktKU+IqIiIiIiEhKU+IrIiIiIiIiKU2Jr4iIiIiIiKQ0Jb4iIiIiIiKS0pT4SsKY2ZfMzM3s+2HHEgaLuNjMXjazHWa20czuMbNeYceWSGZ2lpnNM7PdZvacmY0PO6YwmdlpZvaMmW02s+1m9riZDQs7rrCZ2enRz4u7wo5FJBmpTVWbCmpTg9Setk1taoQSX0mI6AfPj4HysGMJ2WeBO4CjgMuB2cCjoUaUQGZ2BvAYcB8wE9gBvGpmA0MNLFxXAU8AxwFnAMOA58wsO9SoQmRmxcBtwO6wYxFJRmpTm6hNVZsapPa0FWpTm5m7hx2DpDgzM+A54EPgdOBud/9BuFElBzP7LHAnMMrdV4UdT7yZ2dvASne/OLqdB2wGfuXuN4cZW7Iws08CLwAnufuLYccTBjO7ExgA5ANr3P2akEMSSRpqU9umNlVtapDa0wi1qc3U4yuJcC0wGbg55DiS0Y7o9+6hRpEAZjYCmAY8vm+fu1cDLxO5MysRafM30RozOw24FPh62LGIJCm1qW1Lm89PtantkjZ/D21RmxpLia/ElZkNBf4H+KK7bws7niR0LFAFLAs7kATYN8+m5V34dcDgBMeSzI4FHHgv7EASzcyKiAzH+ra7fxB2PCLJRm3qAalNVZsalLbtKahNbY0SX4m3O4DH3f2xsANJNmY2GPgy8Bt3rwk7ngToH/2+p8X+cqBHgmNJSmbWDfgO8Ii7rw07nhD8jMhF3C/DDkQkSalNbYPa1CZqU1F7GqU2tQUlvnJQzOyiaHW4/X19HxgLXB92vPHWzt/HCYHji4G/A8uBH4YWeGLtG3JU1GJ/MbAlwbEknWjxjUeBeuCrIYeTcGZ2KnAhcLWr+ISkGbWpsdSmtova1Dake3sKalPbouJWEjdmtgIYDlQHdhcCdUCtu6flnIto8YlngUHAse6+NeSQEsLMxgIfAJ9x90cD+58Eerj7caEFFzIzywAeAE4h8jexPOSQEs7MngM+QWSY4j75RIapVQPj3X1dGLGJJAO1qa1Tm6o2NUjtaYTa1NZlhR2ApLST+fjf2CtE7sLdmvhwwmdmmcCfidy1n5UuDXTU8ujXOUSXmzCzQuBE4HvhhZUUbgXOBj6Zro00kWUo8lvsewDYCPwHsCnhEYkkF7WpLahNVZvaCrWnEWpTW6HEV+KmtTkVZlYH7HL3lSGEFKpoA30PkUbpbGBHdHgWQHWqz0lydzeznwB3mdlrwKvAj4BSIG0XVDez/wK+AFwCLAv8TdS6e1Xbz0wt7v6xRtjMqoC96fh5IdKS2tRYalPVprak9rSZ2tTWaY6vSOJcCFxGZP7NK0QKUuz7+s8Q40oYd7+XyFIc1wNvAHlE7tK3LM6RFszsGCL/9jlE7tgH/yZ+F2JoIiLJTm2q2tQmak+lPTTHV0RERERERFKaenxFREREREQkpSnxFRERERERkZSmxFdERERERERSmhJfERERERERSWlKfEVERERERCSlKfEVERERERGRlKbEV0RERERERFKaEl8RERERERFJaUp8RUREREREJKUp8RUREREREZGUpsRXREREREREUpoSXxEREREREUlpSnxF5JCY2Xtm5mZ2Wyee8wYzKzOzbmb2/ej52/rq1VmvKyIiEia1qSLxkxV2ACLSdZnZkcBkYC9wkZnd4O5VnXDqa4E/uXu5me3b93vglVaOLeuE1xMREQmV2lSR+FLiKyKH4mqgGvgO8GvgXODBQzmhmR0LjAOuaPHQm+7+0KGcW0REJIldjdpUkbjRUGcROShmlg1cBjwD3A/UAld1wqmvBd5x93c74VwiIiJJT22qSPwp8RWRg3UG0Bt4xN33AM8Bp5jZwIM9oZmVABcAf+icEEVERLoEtakicabEV0QO1tVABfDX6Pb9RD5TLj+Ec14B1HOIQ7tERES6mKtRmyoSV0p8RaTDzKwP8GngUXeviO5+gkhBjqsP4dTXAvcHzikiIpLS1KaKJIYSXxE5GJcB2cB9+3ZEK08+DBxhZtP37TezLDP7lpktNrPdZjbfzE5ueUIzOwaYAHTaEg4iIiJdgNpUkQRQ4isiB+NqYCPwYov990a/BwtyNAAlwFlAH+Bu4G+tzFu6lkiVyfc6O1gREZEkdjVqU0XiTomviHSImU0hss7gA+7e2OLh14DVwCVmlgPgEd9099XuXu/uvyayXMPJgXMWAxeiAhwiIpJG1KaKJI4SXxHpqKuj3+9r+YC7e3R/T+DM1p5sZplAPlAZ2H05UAdoTUEREUknV0e/q00ViTMlviLSbtF1Bi8FFrj74jYO29d4X93G4+cABrwe2LevAEdl608RERFJLWpTRRIrK+wARKRLOZPIOoMfmNl39nPcLuB0M+vr7tv27TSzvsDPgF+4+5bovhnAJA5tyQYREZGuRm2qSAIp8RWRjpgU/X589OtAjgC2AZhZAfB3YDHw3cAx1wJz3X1RJ8YpIiKS7NSmiiSQRaYPiIjET3QO0uPAIGB2cE1BM/t3Io306209X0RERCLUpoocHCW+IhJ3ZnY7MBs4Yd9wLJH/394d20AIxUAU3N8jXVxfJBSILvqkpCAhxGomdeLMepEBuM5NhXuEL/CoMcaSZEuyJ/mfRuuc8/fOVgDwPW4q3Cd8AQAAqOadEQAAANWELwAAANWELwAAANWELwAAANWELwAAANWELwAAANWELwAAANWELwAAANWELwAAANWELwAAANWELwAAANWELwAAANWELwAAANWELwAAANWELwAAANWELwAAANUOWYKa0nDvTHAAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1, 2, figsize=(6.25, 2.75))\n", "ax[0].plot(dets, it.operands[1])\n", "ax[1].plot(dets, it.operands[2])\n", "\n", "[ax_i.set_xlabel('$\\Delta_2/\\Gamma$') for ax_i in ax];\n", "ax[0].set_ylabel('$\\\\omega/\\sqrt{k \\mu_B B\\'/m}$')\n", "ax[1].set_ylabel('$\\\\beta/\\hbar k^2$')\n", "\n", "fig.subplots_adjust(left=0.08, wspace=0.25)" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.6" } }, "nbformat": 4, "nbformat_minor": 4 }