00-Godiva.ipynb

{
  "nbformat": 4,
  "nbformat_minor": 0,
  "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.6.5"
    },
    "colab": {
      "name": "Godiva",
      "version": "0.3.2",
      "provenance": [],
      "collapsed_sections": [
        "4x3hJX0SPOz2"
      ]
    }
  },
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "4x3hJX0SPOz2",
        "colab_type": "text"
      },
      "source": [
        "## Instalación de OpenMC"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "IatJFr3fPNkz",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "#\n",
        "# Executing this cell you will install OpenMC and the nuclear\n",
        "# data libraries in this instance of the Google Colaboratory virtual machine.\n",
        "# The process takes about 2 minutes.\n",
        "#\n",
        "\n",
        "def install_openmc():\n",
        "  #\n",
        "  # Clone source code from Github, make and install\n",
        "  #\n",
        "  \n",
        "  import os\n",
        "  \n",
        "  if not os.path.isdir('/content'):\n",
        "    print(\"Esta función instala OpenMC en una instancia de Google Colaboratory.\")\n",
        "    print(\"Para instalar localmente siga las instrucciones de la documentacion:\")\n",
        "    print(\"http://docs.openmc.org/en/stable/quickinstall.html\")\n",
        "    return\n",
        "  \n",
        "  %cd -q /content\n",
        "  print(\"Obtaining OpenMC source code from Github...\")\n",
        "  !git  --no-pager clone https://github.com/mit-crpg/openmc.git &> /dev/null\n",
        "  %cd -q openmc\n",
        "  !git --no-pager checkout master &> /dev/null\n",
        "  !mkdir build\n",
        "  %cd -q build\n",
        "  print(\"Running cmake...\")\n",
        "  !cmake ../ -DPYTHON_EXECUTABLE=/usr/bin/python3 -DCMAKE_INSTALL_PREFIX=/usr/local &> /dev/null\n",
        "  print(\"Running make...\")\n",
        "  !make -j &> /dev/null\n",
        "  print(\"Running make install...\")\n",
        "  !make install &> /dev/null\n",
        "  import sys\n",
        "  sys.path.append('/usr/local/lib/python3.6/dist-packages/openmc-0.10.0-py3.6-linux-x86_64.egg')\n",
        "  sys.path.append('/usr/local/lib/python3.6/site-packages/')\n",
        "  %cd -q /content\n",
        "\n",
        "def install_data_from_onedrive():\n",
        "  #\n",
        "  # Download preprocessed HDF5 files from OneDrive (faster).\n",
        "  #\n",
        "  import os\n",
        "  \n",
        "  if not os.path.isdir('/content'):\n",
        "    print(\"Esta función instala los datos nucleares de OpenMC en una instancia de Google Colaboratory.\")\n",
        "    print(\"Para instalar localmente siga las instrucciones de la documentacion:\")\n",
        "    print(\"http://docs.openmc.org/en/stable/quickinstall.html\")\n",
        "    return\n",
        "  %cd -q /content\n",
        "  print(\"Obtaining HDF5 files from OneDrive...\")\n",
        "  !wget -O nndc_hdf5.tar.gz \"https://onedrive.live.com/download?cid=22422A8EEA2A85B9&resid=22422A8EEA2A85B9%21105&authkey=AHL6xwYFXDwEzkk\" &> /dev/null\n",
        "  print(\"Uncompressing...\")\n",
        "  !tar xzvf nndc_hdf5.tar.gz &> /dev/null\n",
        "    \n",
        "from time import time\n",
        "t1 = time()\n",
        "install_openmc()\n",
        "install_data_from_onedrive()\n",
        "t2 = time()\n",
        "print(\"Installed OpenMC in {:.2f} minutes\".format((t2-t1)/60.0))\n"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "7ZTpRRpRPG1U",
        "colab_type": "text"
      },
      "source": [
        "# Ejemplo: benchmark Godiva."
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "sDM-e2nUPG1V",
        "colab_type": "text"
      },
      "source": [
        "Realizaremos la simulación del reactor [Godiva](https://en.wikipedia.org/wiki/Godiva_device). Los detalles de este reactor pueden encontrarse en el benchmark HEU-MET-FAST-001 de [ICSBEP](https://www.oecd-nea.org/science/wpncs/icsbep/) y también en el informe \"S. Frankle. A Suite of Criticality Benchmarks for Validating Nuclear Data. LA-13594 (1999)\", [disponible en la base de datos INIS de IAEA](https://inis.iaea.org/search/search.aspx?orig_q=RN:32027510). \n",
        "\n",
        " <p align=\"center\">\n",
        "<img src=\"https://upload.wikimedia.org/wikipedia/commons/2/2a/Godiva-before-scrammed.jpg\">\n",
        "</p>\n",
        "\n",
        "\n",
        "Godiva es una esfera de $8.7407$ cm de radio de uranio metálico de densidad $18.74$ gm/cm$^3$, $94.73$ % en peso de $^{235}$U, $5.27$ % en peso de $^{238}$U, y $1.02$ % en peso de $^{234}$U."
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "x0pWcw42PG1X",
        "colab_type": "code",
        "outputId": "87846f7f-5154-4d13-ab93-486bf492c118",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 1000
        }
      },
      "source": [
        "import openmc\n",
        "\n",
        "umet = openmc.Material()\n",
        "umet.add_nuclide(\"U235\", 0.9473, 'wo')\n",
        "umet.add_nuclide(\"U238\", 0.0527, 'wo')\n",
        "umet.add_nuclide(\"U234\", 0.0102, 'wo')\n",
        "umet.set_density(\"g/cm3\", 18.74)\n",
        "\n",
        "mats = openmc.Materials([umet])\n",
        "mats.cross_sections = \"/content/nndc_hdf5/cross_sections.xml\"\n",
        "mats.export_to_xml()\n",
        "\n",
        "esfera = openmc.Sphere(R=8.74, boundary_type=\"vacuum\")\n",
        "godiva = openmc.Cell(region = -esfera, fill=umet)\n",
        "simulacion = openmc.Universe(cells=[godiva])\n",
        "\n",
        "geom = openmc.Geometry(simulacion)\n",
        "geom.export_to_xml()\n",
        "\n",
        "run = openmc.Settings()\n",
        "run.run_mode = \"eigenvalue\"\n",
        "run.particles = 10000\n",
        "run.batches = 110\n",
        "run.inactive = 10\n",
        "run.export_to_xml()\n",
        "\n",
        "!rm statepoint*.h5\n",
        "!rm summary.h5\n",
        "openmc.run()"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "\n",
            "                               %%%%%%%%%%%%%%%\n",
            "                          %%%%%%%%%%%%%%%%%%%%%%%%\n",
            "                       %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
            "                     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
            "                   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
            "                  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
            "                                   %%%%%%%%%%%%%%%%%%%%%%%%\n",
            "                                    %%%%%%%%%%%%%%%%%%%%%%%%\n",
            "                ###############      %%%%%%%%%%%%%%%%%%%%%%%%\n",
            "               ##################     %%%%%%%%%%%%%%%%%%%%%%%\n",
            "               ###################     %%%%%%%%%%%%%%%%%%%%%%%\n",
            "               ####################     %%%%%%%%%%%%%%%%%%%%%%\n",
            "               #####################     %%%%%%%%%%%%%%%%%%%%%\n",
            "               ######################     %%%%%%%%%%%%%%%%%%%%\n",
            "               #######################     %%%%%%%%%%%%%%%%%%\n",
            "                #######################     %%%%%%%%%%%%%%%%%\n",
            "                ######################     %%%%%%%%%%%%%%%%%\n",
            "                 ####################     %%%%%%%%%%%%%%%%%\n",
            "                   #################     %%%%%%%%%%%%%%%%%\n",
            "                    ###############     %%%%%%%%%%%%%%%%\n",
            "                      ############     %%%%%%%%%%%%%%%\n",
            "                         ########     %%%%%%%%%%%%%%\n",
            "                                     %%%%%%%%%%%\n",
            "\n",
            "                   | The OpenMC Monte Carlo Code\n",
            "         Copyright | 2011-2018 Massachusetts Institute of Technology\n",
            "           License | http://openmc.readthedocs.io/en/latest/license.html\n",
            "           Version | 0.10.0\n",
            "          Git SHA1 | 455efffdbf6eacbc0b3b93ad61cfeb24490897d8\n",
            "         Date/Time | 2019-07-25 15:56:04\n",
            "    OpenMP Threads | 2\n",
            "\n",
            " Reading settings XML file...\n",
            " Reading cross sections XML file...\n",
            " Reading materials XML file...\n",
            " Reading geometry XML file...\n",
            " Building neighboring cells lists for each surface...\n",
            " Reading U235 from /content/nndc_hdf5/U235.h5\n",
            " Reading U238 from /content/nndc_hdf5/U238.h5\n",
            " Reading U234 from /content/nndc_hdf5/U234.h5\n",
            " Maximum neutron transport energy: 2.00000E+07 eV for U235\n",
            " Writing summary.h5 file...\n",
            " Initializing source particles...\n",
            "\n",
            " ====================>     K EIGENVALUE SIMULATION     <====================\n",
            "\n",
            "  Bat./Gen.      k            Average k         \n",
            "  =========   ========   ====================   \n",
            "        1/1    1.34507                       \n",
            "        2/1    1.15804                       \n",
            "        3/1    1.06264                       \n",
            "        4/1    1.03493                       \n",
            "        5/1    1.02419                       \n",
            "        6/1    1.01023                       \n",
            "        7/1    1.00572                       \n",
            "        8/1    0.99023                       \n",
            "        9/1    0.98896                       \n",
            "       10/1    1.00114                       \n",
            "       11/1    0.99771                       \n",
            "       12/1    1.00346    1.00058 +/- 0.00288\n",
            "       13/1    0.99742    0.99953 +/- 0.00197\n",
            "       14/1    0.99768    0.99907 +/- 0.00147\n",
            "       15/1    0.99863    0.99898 +/- 0.00114\n",
            "       16/1    0.98375    0.99644 +/- 0.00270\n",
            "       17/1    1.01690    0.99936 +/- 0.00371\n",
            "       18/1    1.02803    1.00295 +/- 0.00481\n",
            "       19/1    1.00030    1.00265 +/- 0.00425\n",
            "       20/1    1.01242    1.00363 +/- 0.00393\n",
            "       21/1    0.99448    1.00280 +/- 0.00365\n",
            "       22/1    0.99304    1.00199 +/- 0.00343\n",
            "       23/1    0.99102    1.00114 +/- 0.00327\n",
            "       24/1    0.99493    1.00070 +/- 0.00306\n",
            "       25/1    1.00905    1.00125 +/- 0.00290\n",
            "       26/1    1.00357    1.00140 +/- 0.00272\n",
            "       27/1    1.00703    1.00173 +/- 0.00257\n",
            "       28/1    1.00084    1.00168 +/- 0.00243\n",
            "       29/1    1.00458    1.00183 +/- 0.00230\n",
            "       30/1    0.99937    1.00171 +/- 0.00218\n",
            "       31/1    0.99828    1.00155 +/- 0.00208\n",
            "       32/1    1.00804    1.00184 +/- 0.00201\n",
            "       33/1    0.99480    1.00154 +/- 0.00194\n",
            "       34/1    1.00485    1.00167 +/- 0.00187\n",
            "       35/1    0.99129    1.00126 +/- 0.00184\n",
            "       36/1    0.99419    1.00099 +/- 0.00179\n",
            "       37/1    0.99293    1.00069 +/- 0.00174\n",
            "       38/1    1.00243    1.00075 +/- 0.00168\n",
            "       39/1    0.99508    1.00055 +/- 0.00164\n",
            "       40/1    0.98780    1.00013 +/- 0.00164\n",
            "       41/1    0.97851    0.99943 +/- 0.00173\n",
            "       42/1    1.00367    0.99956 +/- 0.00168\n",
            "       43/1    0.99973    0.99957 +/- 0.00163\n",
            "       44/1    1.02010    1.00017 +/- 0.00169\n",
            "       45/1    0.99747    1.00010 +/- 0.00164\n",
            "       46/1    0.98673    0.99972 +/- 0.00164\n",
            "       47/1    1.01014    1.00001 +/- 0.00162\n",
            "       48/1    1.00330    1.00009 +/- 0.00158\n",
            "       49/1    1.00637    1.00025 +/- 0.00155\n",
            "       50/1    0.99570    1.00014 +/- 0.00151\n",
            "       51/1    1.01667    1.00054 +/- 0.00153\n",
            "       52/1    1.00615    1.00068 +/- 0.00150\n",
            "       53/1    0.99965    1.00065 +/- 0.00146\n",
            "       54/1    0.99454    1.00051 +/- 0.00144\n",
            "       55/1    1.01277    1.00079 +/- 0.00143\n",
            "       56/1    1.00264    1.00083 +/- 0.00140\n",
            "       57/1    0.99655    1.00074 +/- 0.00137\n",
            "       58/1    0.99158    1.00054 +/- 0.00136\n",
            "       59/1    0.99355    1.00040 +/- 0.00134\n",
            "       60/1    1.00169    1.00043 +/- 0.00131\n",
            "       61/1    1.00165    1.00045 +/- 0.00128\n",
            "       62/1    0.99548    1.00036 +/- 0.00126\n",
            "       63/1    0.98849    1.00013 +/- 0.00126\n",
            "       64/1    0.99341    1.00001 +/- 0.00124\n",
            "       65/1    1.00728    1.00014 +/- 0.00123\n",
            "       66/1    0.98586    0.99989 +/- 0.00123\n",
            "       67/1    1.00130    0.99991 +/- 0.00121\n",
            "       68/1    0.99066    0.99975 +/- 0.00120\n",
            "       69/1    0.98908    0.99957 +/- 0.00119\n",
            "       70/1    1.00515    0.99966 +/- 0.00117\n",
            "       71/1    0.99426    0.99957 +/- 0.00116\n",
            "       72/1    0.99851    0.99956 +/- 0.00114\n",
            "       73/1    0.99706    0.99952 +/- 0.00112\n",
            "       74/1    0.99700    0.99948 +/- 0.00111\n",
            "       75/1    1.02953    0.99994 +/- 0.00118\n",
            "       76/1    1.00801    1.00006 +/- 0.00117\n",
            "       77/1    0.99137    0.99993 +/- 0.00116\n",
            "       78/1    1.00339    0.99998 +/- 0.00114\n",
            "       79/1    1.00177    1.00001 +/- 0.00113\n",
            "       80/1    1.00707    1.00011 +/- 0.00112\n",
            "       81/1    0.99500    1.00004 +/- 0.00110\n",
            "       82/1    1.00581    1.00012 +/- 0.00109\n",
            "       83/1    1.00248    1.00015 +/- 0.00108\n",
            "       84/1    0.99687    1.00011 +/- 0.00106\n",
            "       85/1    0.99779    1.00008 +/- 0.00105\n",
            "       86/1    0.99561    1.00002 +/- 0.00104\n",
            "       87/1    0.99045    0.99989 +/- 0.00103\n",
            "       88/1    1.01384    1.00007 +/- 0.00103\n",
            "       89/1    1.00163    1.00009 +/- 0.00102\n",
            "       90/1    0.99851    1.00007 +/- 0.00101\n",
            "       91/1    1.00222    1.00010 +/- 0.00099\n",
            "       92/1    0.99598    1.00005 +/- 0.00098\n",
            "       93/1    0.99877    1.00003 +/- 0.00097\n",
            "       94/1    0.99637    0.99999 +/- 0.00096\n",
            "       95/1    0.99264    0.99990 +/- 0.00095\n",
            "       96/1    0.98645    0.99975 +/- 0.00096\n",
            "       97/1    1.01379    0.99991 +/- 0.00096\n",
            "       98/1    1.01865    1.00012 +/- 0.00097\n",
            "       99/1    1.00764    1.00020 +/- 0.00096\n",
            "      100/1    1.02300    1.00046 +/- 0.00099\n",
            "      101/1    1.01066    1.00057 +/- 0.00098\n",
            "      102/1    0.99635    1.00052 +/- 0.00097\n",
            "      103/1    1.00169    1.00054 +/- 0.00096\n",
            "      104/1    1.00581    1.00059 +/- 0.00095\n",
            "      105/1    1.00616    1.00065 +/- 0.00094\n",
            "      106/1    1.00969    1.00075 +/- 0.00094\n",
            "      107/1    1.00227    1.00076 +/- 0.00093\n",
            "      108/1    1.01493    1.00091 +/- 0.00093\n",
            "      109/1    0.98667    1.00076 +/- 0.00093\n",
            "      110/1    0.99515    1.00071 +/- 0.00093\n",
            " Creating state point statepoint.110.h5...\n",
            "\n",
            " =======================>     TIMING STATISTICS     <=======================\n",
            "\n",
            " Total time for initialization     =  1.8327E-01 seconds\n",
            "   Reading cross sections          =  1.5266E-01 seconds\n",
            " Total time in simulation          =  4.3532E+00 seconds\n",
            "   Time in transport only          =  4.2522E+00 seconds\n",
            "   Time in inactive batches        =  3.9413E-01 seconds\n",
            "   Time in active batches          =  3.9591E+00 seconds\n",
            "   Time synchronizing fission bank =  4.5787E-02 seconds\n",
            "     Sampling source sites         =  3.1591E-02 seconds\n",
            "     SEND/RECV source sites        =  1.3946E-02 seconds\n",
            "   Time accumulating tallies       =  1.2344E-04 seconds\n",
            " Total time for finalization       =  1.0526E-05 seconds\n",
            " Total time elapsed                =  4.5843E+00 seconds\n",
            " Calculation Rate (inactive)       =  2.53722E+05 neutrons/second\n",
            " Calculation Rate (active)         =  2.52582E+05 neutrons/second\n",
            "\n",
            " ============================>     RESULTS     <============================\n",
            "\n",
            " k-effective (Collision)     =  1.00080 +/-  0.00098\n",
            " k-effective (Track-length)  =  1.00071 +/-  0.00093\n",
            " k-effective (Absorption)    =  1.00222 +/-  0.00113\n",
            " Combined k-effective        =  1.00137 +/-  0.00067\n",
            " Leakage Fraction            =  0.57186 +/-  0.00049\n",
            "\n"
          ],
          "name": "stdout"
        },
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "0"
            ]
          },
          "metadata": {
            "tags": []
          },
          "execution_count": 1
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "3d8jOVL3PG1u",
        "colab_type": "code",
        "outputId": "5425ab32-221e-4a96-e8f7-d7e8d635d136",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 301
        }
      },
      "source": [
        "% pylab inline\n",
        "estado = openmc.StatePoint(\"statepoint.110.h5\")\n",
        "keff = estado.k_generation\n",
        "a = plot(keff)\n",
        "a = xlabel(\"Nro. de generación\")\n",
        "a = ylabel(\"Estimador de keff\")"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "Populating the interactive namespace from numpy and matplotlib\n"
          ],
          "name": "stdout"
        },
        {
          "output_type": "display_data",
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEKCAYAAAAB0GKPAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3Xl4VOX1wPHvyWSDbCQk7ISwQ9gh\nAoIiilKsKO67IlWxra12sf3V1rV2UdtqW6tVqqhYxbqgUhcULYIiW9hX2cIStiQEEpJA1vP7Y27C\nhGQmA2QySTif58mTzL137pybSebc9z3vfa+oKsYYY0xdQoIdgDHGmKbBEoYxxhi/WMIwxhjjF0sY\nxhhj/GIJwxhjjF8sYRhjjPGLJQxjjDF+sYRhjDHGL5YwjDHG+CU02AHUp8TERE1JSQl2GMYY02Qs\nX748R1WT/Nm2WSWMlJQU0tPTgx2GMcY0GSKy099trUvKGGOMXyxhGGOM8UvAEoaITBeRLBFZ52X9\nJBFZIyKrRCRdRM7xWFfuLF8lIrMDFaMxxhj/BbKG8QrwD2CGl/VfALNVVUVkIPAW0MdZd1RVBwcw\nNmOMMScpYC0MVV0A5PpYX6DHb8YRBdiNOYwxphELag1DRK4QkU3AR8D3PFZFOt1Ui0Xk8jr2MdXZ\nNj07Ozug8RpjzJksqAlDVd9T1T7A5cBjHqu6qGoacCPwVxHp7mMf01Q1TVXTkpL8Gkpc2z4or7AG\njjHG+NIoRkk53VfdRCTRebzH+b4d+BIYEsjX7/vQHJ6csymQL2GMMU1e0BKGiPQQEXF+HgpEAAdF\nJF5EIpzlicBoYEMgYwlzhVBcVhHIlzDGmCYvYKOkRGQmMBZIFJFM4GEgDEBVnweuAm4VkVLgKHCd\nM2KqL/CCiFTgTmiPq2pAE0ZEaAgl5ZYwjDHGl4AlDFW9oY71TwBP1LL8G2BAoOKqTbgrhBJrYRhj\njE+NooYRbOGhljCMMaYuljCwhGGMMf6whIGTMKyGYYwxPlnCwGoYxhjjD0sYWJeUMcb4wxIGEB7q\noti6pIwxxidLGLi7pIpLy4MdhjHGNGqWMLAL94wxxh+WMLAahjHG+MMSBk4LwxKGMcb4ZAkDuw7D\nGGP8YQkDuw7DGGP8YQkDq2EYY4w/LGHgThhlFUqF3XXPGGO8soSBO2EAVscwxhgfLGHgrmEAdtc9\nY4zxwRIG7mG1gNUxjDHGB0sYHO+SKi6z6UGMMcabgCYMEZkuIlkiss7L+kkiskZEVolIuoic47Fu\nsohscb4mBzLOcGthGGNMnQLdwngFmOBj/RfAIFUdDHwPeBFARBKAh4ERwHDgYRGJD1SQ4S4XYEVv\nY4zxJaAJQ1UXALk+1heoauVY1iig8ufvAHNVNVdVDwFz8Z14TovVMIwxpm5Br2GIyBUisgn4CHcr\nA6AjsNtjs0xnWW3Pn+p0Z6VnZ2efUgzWJWWMMXULesJQ1fdUtQ9wOfDYKTx/mqqmqWpaUlLSKcVg\nCcMYY+oW9IRRyem+6iYiicAeoLPH6k7OsoCoGiVlNQxjjPEqqAlDRHqIiDg/DwUigIPAp8B4EYl3\nit3jnWUBUXnhnrUwjDHGu9BA7lxEZgJjgUQRycQ98ikMQFWfB64CbhWRUuAocJ1TBM8VkceAZc6u\nfquqXovnp8uK3sYYU7eAJgxVvaGO9U8AT3hZNx2YHoi4TmQ1DGOMqVujqWEEk00+aIwxdbOEgdUw\njDHGH5YwsLmkjDHGH5YwsBqGMcb4wxIG1iVljDH+sIQBiAjhoSF24Z4xxvhgCcMR4QqxFoYxxvhg\nCcMRHmoJwxhjfLGE4bCEYYwxvlnCcISHhtiFe8YY44MlDEe41TCMMcYnSxgO65IyxhjfLGE4rEvK\nGGN8s4ThCHeFUFxqCcMYY7yxhOGwC/eMMcY3SxiOCKthGGOMT5YwHO6it81Wa4wx3ljCcESEuqzo\nbYwxPgQsYYjIdBHJEpF1XtbfJCJrRGStiHwjIoM81u1wlq8SkfRAxejJrsMwxhjfAtnCeAWY4GN9\nBnCeqg4AHgOmnbD+fFUdrKppAYqvGrsOwxhjfAsN1I5VdYGIpPhY/43Hw8VAp0DF4g9LGMYY41tj\nqWHcDnzi8ViBz0RkuYhM9fVEEZkqIukikp6dnX3KAdiFe8YY41vAWhj+EpHzcSeMczwWn6Oqe0Sk\nDTBXRDap6oLanq+q03C6s9LS0vRU4wh3hVBarlRUKCEhcqq7McaYZiuoLQwRGQi8CExS1YOVy1V1\nj/M9C3gPGB7oWKru622tDGOMqVXQEoaIJAOzgFtUdbPH8igRian8GRgP1DrSqj5FWMIwxhifAtYl\nJSIzgbFAoohkAg8DYQCq+jzwENAaeE5EAMqcEVFtgfecZaHAG6o6J1BxVqpsYRSXVkBkoF/NGGOa\nHq8JQ0TuVdW/ichoVV14sjtW1RvqWH8HcEcty7cDg2o+I7DCXdbCMMYYX3x1SU1xvj/TEIEEW1UN\nw4bWGmNMrXx1SW0UkS1ABxFZ47FcAFXVgYENrWFFhLoASxjGGOON14ShqjeISDvgU+CyhgspOKyF\nYYwxvvmqYXyhquNE5FNV3dmQQQXD8WG1NmOtMcbUxleXVHsRGQVc6ox4qnY1m6quCGhkDayy6F1s\nLQxjjKmVr4TxEPAg7jmenjphnQIXBCqoYLAuKWOM8c1XDeMd4B0ReVBVH2vAmIIiwhKGMcb45M+V\n3r8XkZtF5CFwX6EtIgGfqqOh2dQgxhjjmz8J41ngbKDyQrwjzrJmperCPWthGGNMrfyZGmSEqg4V\nkZUAqnpIRMIDHFeDq5oaxBKGMcbUyp8WRqmIuHAXuhGRJKDZfapa0dsYY3zzJ2H8HfcU421E5PfA\n18AfAxpVEFjCMMYY3+rsklLV10VkOTAO97UYl6vqxoBH1sBs8kFjjPGtzoQhIrer6kvAJo9lj6vq\nrwIaWQOLsBqGMcb45E/R+yoROaaqrwOIyLM0wztGiAjhrhDrkjLGGC/8ShjAbBGpACYAh1X19sCG\nFRzhoZYwjDHGG1+TDyZ4PLwDeB9YCDwqIgmqmhvo4BpaeGiITT5ojDFe+GphLMc9lFY8vl/ifCnQ\nLeDRNTDrkjLGGO+8DqtV1a6q2u2E75VfdSYLEZkuIlkiss7L+ptEZI2IrBWRb0RkkMe6CSLyrYhs\nFZEGK65bl5Qxxnjnz3UYp+oV3DUPbzKA81R1APAYMA3AuUjwWeBiIBW4QURSAxhnFXeXlCUMY4yp\nTcAShqouALzWOVT1G1U95DxcjHsadYDhwFZV3a6qJcCbwKRAxenJuqSMMca7QLYwTsbtwCfOzx2B\n3R7rMp1ltRKRqSKSLiLp2dnZpxVEeGiIXYdhjDFe1JkwxC1g05uLyPm4E8b/ncrzVXWaqqapalpS\nUtJpxWIJwxhjvPOnhfEcAZreXEQGAi8Ck1T1oLN4D9DZY7NOzrKAi7CitzHGeOVPwhihqncDx8A9\nvTlw2tObi0gyMAu4RVU3e6xaBvQUka7ONOrXA7NP9/X8YTUMY4zxzp8rvU9penMRmQmMBRJFJBN4\nGAgDUNXncd8zvDXwnIgAlDldS2Ui8iPgU8AFTFfV9Sd7YKciIsxGSRljjDf+JIwTpze/Gnigriep\n6g11rL8D9xXkta37GPjYj9jqlbUwjDHGO5ve3INduGeMMd75O5dUFjDTc13znUvKEoYxxtTG37mk\nkoFDzs+tgF1A14BH18DCXS5rYRhjjBd1ziUFfA5cqqqJqtoamAh81lABNiTrkjLGGO/8GVY70ilC\nA6CqnwCjAhdS8FR2SalqsEMxxphGx5+EsVdEHhCRFOfrN8DeQAcWDHabVmOM8c6fhHEDkIR7aO0s\n52efQ2abqnCX+9dhhW9jjKnJn2G1ucC9DRBL0IU7LQyrYxhjTE2NZbbaRsEShjHGeGcJw0OEJQxj\njPHKZ8IQEZeI/LShggm2qhaG1TCMMaYGnwlDVctppgXu2lQVva2FYYwxNfgz+eBCEfkH8B+gsHKh\nqq4IWFRBEm7Dao0xxit/EsZg5/tvPZYpcEH9hxNcVvQ2xhjv/BlWe35DBNIYRFgNwxhjvPLnnt5x\nIvKUiKQ7X38RkbiGCK6hhbtcgLUwjDGmNv4Mq52O+z7e1zpf+cDLgQwqWCLCKmsY5UGOxBhjGh9/\nahjdVfUqj8ePisiqQAUUTC3D3S2MomJLGMYYcyJ/WhhHReScygciMho4WteTRGS6iGSJyDov6/uI\nyCIRKRaR+05Yt0NE1orIKhFJ9yPGehETEQbAkeKyhnpJY4xpMvxpYfwAeNWpWwiQC9zmx/NeAf4B\nzPCyPhe4B7jcy/rzVTXHj9epN1ER7hZGwTFLGMYYcyJ/RkmtAgaJSKzzON+fHavqAhFJ8bE+C8gS\nkUv8CzXwQl0htAhzUVBcGuxQjDGm0fF1T++feVkOgKo+FaCYwH2dx2ciosALqjrN24YiMhWYCpCc\nnHzaLxwdGUqBdUkZY0wNvloYMc733sBZwGzn8aXA0kAGBZyjqntEpA0wV0Q2qeqC2jZ0ksk0gLS0\ntNO+VV5MRChHrEvKGGNq8JowVPVRABFZAAxV1SPO40eAjwIZlKrucb5nich7wHCg1oRR36yFYYwx\ntfNnlFRboMTjcYmzLCBEJEpEYip/BsYDtY60CoToiFArehtjTC38GSU1A1jqnOmDe1TTq3U9SURm\nAmOBRBHJBB4GwgBU9XkRaQekA7FAhYj8BEgFEoH3nFpJKPCGqs45mYM6HdERoewqLGqolzPGmCbD\nn1FSvxeROUDltRhTVHWlH8/zOS26qu4HOtWyKh8YVNf+AyU60moYxhhTG39aGKjqchHZDUQCiEiy\nqu4KaGRBEh0RSmGJJQxjjDmRP5MPXiYiW4AMYL7z/ZNABxYslTUM1dMecGWMMc2KP0Xvx4CRwGZV\n7QpcCCwOaFRBFB0ZSlmF2k2UjDHmBP4kjFJVPQiEiEiIqs4D0gIcV9DERLh76ayOYYwx1flTwzgs\nItG4r4N4XUSy8LhVa3MTHen+lRQUl5EUExHkaIwxpvHwp4UxCffstD8F5gDbcF/t3SxFOzPW2rUY\nxhhTnT/DagsBnMkH/xvwiIIsurJLyiYgNMaYaupMGCJyF/AocAyowD3FuQLdAhtacMRUdklZC8MY\nY6rxp4ZxH9C/oe9NESyVLQybT8oYY6rzp4axDThj5sqoLHoXWsIwxphq/Glh3A98IyJLgOLKhap6\nT8CiCqLjNQxLGMYY48mfhPEC8D9gLe4aRrMWERpCaIhYDcMYY07gT8IIU9Va777XHImI3RPDGGNq\n4U8N4xMRmSoi7UUkofIr4JEFkd0TwxhjavKnhVE5Tfn9Hsua7bBacCcMq2EYY0x1/ly417UhAmlM\nYiKthWGMMSfymjBE5AJV/Z+IXFnbelWdFbiwgis6IpScgpK6NzTGmDOIrxbGebhHR9U2b5QCzTdh\nRIax4+AZc+mJMcb4xWvCUNWHnR9/q6oZnutEpM5uKhGZDkwEslS1fy3r+wAvA0OB36jqnz3WTQD+\nBriAF1X1cT+Opd5ER9goKWOMOZE/o6TerWXZO3487xVggo/1ucA9wJ89F4qIC3gWuBhIBW4QkVQ/\nXq/eREe4rIZhjDEn8FXD6AP0A+JOqGPE4tzb2xdVXSAiKT7WZwFZInLJCauGA1tVdbsTx5u4p1jf\nUNdr1pfoiDCOlpZTVl5BqMufnGqMMc2frxpGb9xdSq2oXsc4AtwZwJg6Ars9HmcCI7xtLCJTgakA\nycnJ9RLA8fmkyolraQnDGGPAdw3jA+ADETlbVRc1YEwnRVWnAdMA0tLStD72GeNxT4y4lmH1sUtj\njGny/Dl9vkJEYkUkTES+EJFsEbk5gDHtATp7PO7kLGswnrdpNcYY4+ZPwhivqvm4u6d2AD2AXwQw\npmVATxHpKiLhwPXA7AC+Xg1V98SwwrcxxlTxa/JB5/slwNuqmicidT5JRGYCY4FEEckEHq7cl6o+\nLyLtgHTcRfQKEfkJkKqq+SLyI+BT3MNqp6vq+pM7rNNT2cKw6UGMMeY4fxLGf0VkE3AU+IGIJOG+\nXatPqnpDHev34+5uqm3dx8DHfsQWEDHWwjDGmBrq7JJS1V8Bo4A0VS3Fffe9SYEOLJishmGMMTV5\nTRgi8kuPh+NUtRxAVQtxX3DXbFXWMOw2rcYYc5yvFsb1Hj/ff8I6X1dwN3lR4U4Nw7qkjDGmiq+E\nIV5+ru1xsxISIkSFu6xLyhhjPPhKGOrl59oeNzvRdk8MY4ypxtcoqUEiko+7NdHC+RnncZ1zSTV1\nNmOtMcZU52tqEFdDBtLYREeG2XUYxhjjwWbW8yImIpSCY6XBDsMYYxoNSxheWJeUMcZUZwnDCyt6\nG2NMdZYwvLAWhjHGVGcJw4vKhKHa7EcQG2OMXyxheBEdGUqFwtHS8mCHYowxjYIlDC8q55Oy6UGM\nMcbNEoYXbWPd1yZmHjoa5EiMMaZxsIThRWqHWAA27MuvY0tjjDkzWMLwokNcJHEtwtiw1xKGMcaA\nJQyvRITU9rHWwjDGGEfAEoaITBeRLBFZ52W9iMjfRWSriKwRkaEe68pFZJXzNTtQMdYltUMsm/bl\nU1ZeEawQjDGm0QhkC+MVfN9o6WKgp/M1Ffinx7qjqjrY+boscCH6lto+luKyCnYcLAxWCMYY02gE\nLGGo6gIg18cmk4AZ6rYYaCUi7QMVz6moLHyvtzqGMcYEtYbREdjt8TjTWQYQKSLpIrJYRC73tRMR\nmepsm56dnV2vAXZPiibMJWzcd6Re92uMMU1RYy16d1HVNOBG4K8i0t3bhqo6TVXTVDUtKSmpXoMI\nDw2hZ5sYK3wbYwzBTRh7gM4ejzs5y1DVyu/bgS+BIQ0dXKXUDrE2tNYYYwhuwpgN3OqMlhoJ5Knq\nPhGJF5EIABFJBEYDG4IVZGr7WHIKisk6cixYIRhjTKPg657ep0VEZgJjgUQRyQQeBsIAVPV54GPg\nu8BWoAiY4jy1L/CCiFTgTmiPq2rwEkblFd9782nTu9nfytwYY7wKWMJQ1RvqWK/A3bUs/wYYEKi4\nTlbf9senCBnbu02QozHGmOBprEXvRiOuRRgdW7WwOoYx5oxnCcMPqR1iWb83326mZIw5o1nC8MOY\nXklk5BSyeLuv6xCNMaZ5s4Thh2uGdSIxOoJn520NdijGGBM0ljD8EBnmYuqYrny9NYeVuw4FO5wa\nrKvMBENGTiG3vLSEvKLSYIdiGoglDD/dNKILrVqGNbpWxoxFOxjzp3mU2oy6poHN25TFV1tyWLCl\nfqfkMY2XJQw/RUWEMmVUVz7fmMXGRjJVyLHScv7+xVZ25x5le7bNqGsa1tbsAgCWZlht70xhCeMk\n3DYqheiIUJ6fvy3YoQDw7opMcgqKAdiwLy/I0ZgzzbYsd8JYtsMShqcjx0rZtL9xnFTWN0sYJyGu\nZRhXDe3IJ+v2k3c0uP225RXKvxZsp1+HWMJDQ2xGXdPgtmUXECKwaf8RDheVBDucRuNPn37Lpc98\nzYH85jedkCWMk3TVsE6UlFXw0Zp9QY1jzrr97DhYxN3n96BPuxi7sNCHp+Zu5tP1+4MdRrNyuKiE\nnIISLujjnv1g2Y7GNxgkGMrKK/h47T5Ky5WZS3cFO5x6ZwnjJA3oGEfPNtG8uyIzaDGoKs/P30bX\nxCi+068dfdu57z1uo6Vqyj9Wyj/+t4XpX2cEO5RmZZtTv7hqaCfCXSFnZLdU+o5cxv5pHv/bdKBq\n2dKMXHIKSmjVMow3luxqdoNRLGGcJBHhqmGdWL7zEBk5wSk0r9x9mLV78rjz3G64QoTUDrHkFpZw\nIL84KPE0ZssycqlQWJ15uF7+edftyWPiM1+xP6/5dTecjK1O/aJfhzgGdY5jyRlW+F6wOZtbXlrK\njoNFPDnn26qTtQ/X7qNluIs/XDGArCPFza5lawnjFFw+uCMhArOC1MqovGVsZXdA5Yy6jWX0VmOy\naNtBAI6VVtRLt92Ha/axbk9+oxn4ECxbswqICA2hY3wLhndNYP2ePAqLy4IdFkszcjlyLLD1xU/X\n7+eOV9NJSYzi/ov7sGn/Eb78Npuy8grmrNvPhX3b8p1+7eic0IIZ3+yst9fdnVvEmCfnsW5P8Aa4\nWMI4Be3iIhndI5FZK/ZQUdHw3UDbsgqICnfRNjYCgD7tYgCq3Rnw9SU7Wb7T+pUXbT9I96QogHr5\nfSzJcCegN5buIquOombe0VJuf2UZmw80vwEJW7MK6JoYhStEGN61NWUVyspdh4MaU25hCddPW8Qf\nPt4UsNdQVe6ftZbe7WJ4886RfO+crnSIi+S5L7eyaPtBcgtLuGRge1whwi0ju7B0R269ncj9e/FO\nduUW8dWWnHrZ36mwhHGKrh7WiT2Hj7Jo+8EGf+1t2QV0bxONiAAQExlGckLLqjPoHTmFPPD+Op6e\nu7nBY2tMDheVsGFfPpMGd6RjqxZeE8ahwhK/riUoLC5jbWYeEwe2p7xCeWHBdp/bf7HxAF9syuKP\nH288pfgbs23ZhfRoEw3A0ORWhAgsDXIdY+WuQ1QovL9yT8CuPt+bd4zcwhKuO6szcS3DCHOFcOeY\nbizbcYjHP9lEVLiL83q5bxV9bVpnIkJDmLHo5FsZxWXlHCstr/b4rfTdQHB7EixhnKLxqe1IjA7n\n3jdXsnp3w55ZbcsqoHtSdLVlqe1jq1oYMxbtRNV9Nnwy3QSqykMfrKv6wzxd5RVaVRwNhsXbc1GF\ns7u3ZliXeNJ35tY6MOAf87Zy7QuLqhUva7N85yHKKpRr0jpz+eCOvL5kJ9lHvNeNvvzWfQX0vG+z\nm1Vr71hpObsPFVUljJjIMPp1iGNpRsOfPHlasesQInC0tJy3l9fP3/CJ1jvdQZXdwADXn5VMQlQ4\n6/fmc1FqWyLDXAC0ahnOxIEd+HD13mof/nX5YuMBznliHte+sIiSMnfd7ZO1+zlUVEpidIQljKao\nRbiLN6eOJDLMxXXTFvFZAxW3CovL2Jt3rKqbpVJqh1h2HCwk68gx3k7fTUrrlpSWKwu3+t98Xb83\nnxmLdvLLd9bwzvLj9ZkFm7P5v3fWnPRY+2f+t4WLnprPliB1ySzefpAWYS4GdWpFWko8B/KL2XP4\naI3tKj/Mf/qf1WQeKqpavnLXIQ4VHj/mJRkHcYUIw7rEc/f53Skpq+DFr2pvZZRXKAu2ZHNx/3a0\njgpvkNZe5qEibpi2mH15NY+xPm3PLkSVaictw7smsGLXYYpK6j5BWb37ML95b22NbU93lN+KnYfp\n3yGOtC7xzFi0k/IAdBdv2JdPiEDfdscTRotwF7eNSgHgkoEdqm1/5dCOHCku4/ON3k9GjhwrZdfB\nItZkHuaX76zm9lfTaRHmYk1mHk9/7v67+ffinaS0bsn1Z3Vme07hSSWg+mQJ4zT0aBPDez8cTe92\nsdz17+WsaICJCStHZtXWwlCFP368iSPFZTxx1UCiI0KZ963/8/y8uyKTcFcII7om8Mt3VvPO8kx+\n9e4abp2+lP+k7+YH/15RbaTR+r15bM2qPRnkFBTzrwXbqVB4c1n1s715m7JYvjPw3ReLth0kLSWe\n8NAQhibHAzXrGMdKy1m/N4+L+7ejvEL50Rsr2ZpVwB2vLuOK577hp2+tqtp2yfZcBnSMIzoilG5J\n0Vw6qAOvfLOj1iLkmszDHC4q5eIB7fn+ed35emsOSwLcffnGkl0s2n6Qzzf4bimdKPtIMffPWut3\nN05lq7GyhQHuARglZRV19q+vyTzMzS8t4fUlu3h3xZ6q5arKrdOXMv7p+af0t1FeoazOPMzQ5FZM\nHpXCrtwi5m/OOun91GX93ny6JkbRItxVbfnUMd3463WDGden+l05R3ZrTbvYSN7zONZKZeUV/P6j\nDQx89DPG/Gkel/1jIe8sz+Tu87sz92djuP6szjw/fxuvLdpB+s5D3DSiC6kdYimvULYcCE7LPaAJ\nQ0Smi0iWiKzzsl5E5O8islVE1ojIUI91k0Vki/M1OZBxno6kmAhev2MELcJcvLv85EZNbc06wl2v\npZ/UGXjlP2v3NtUTRl+nifzeyj0M6BjH8K4JnNMjkS+/zfLrzK20vILZq/Yyrm8bXp5yFkOS47nv\n7dW8lb6bu87rxh+vHMCi7Qd56IP1lJVX8PTczVz6zNdM+sdC0mvpu3523laOlpYzqHMrZq3IpLjM\nfUaUlX+MH7y+nPveXuNXXMdKy/1u2SzYnM3V//yGtZl5HCwo5tsDRxjZrTXgHhgQFe6qkTDW782j\ntFy5fEhHnrhqIKt2H+bCp+azaNtB5/eXzdrMPI6WlLM68zAjuiVUPfeBS1JJiArnrteWk1tYPcYv\nv80mRODcHoncPLILSTER/GXu5oBdK1NeocxyPpROdojrW+m7mbl0Fy8t9O9ala1ZBYhA18Tjrdzh\nXROIjQxlrkeyUlUe/e96Hpm9nm+25rB692FueWkpcS3C6N02hlcWZlT9PpZk5PLVlhwyDx3l6ucX\n8fAH606qO/Xb/UcoKilnSHI8E/q3o01MBK/WMkKpvEJP6+x8w958UjvE1VgeGebi8iEdCQmRastd\nIcKkwR2YvzmbgwXHuy9zC0u4dfpS/vVVBtcM68SfrxnEv25N4/OfnccvvtOHiFAXD0xMpXN8Sx78\nYD3hoSFcPaxT1S2jg9UtFegWxivABB/rLwZ6Ol9TgX8CiEgC8DAwAhgOPCwi8QGN9DRER4RyQZ82\nfLp+P2VexvrP+zaLL789fsazL+8ot760lE/XH2DKK8t89oV72pblno6hS+uW1ZZ3iIskrkUYAJNH\npSAinN8niX15x/jWj4S0YHM2BwtLuHJoJ1qGhzL9trO4/ZyuvPODUdx/cV9uGJ7MD8Z2Z+bSXYx7\naj5/+2ILlw3qQNvYSCZPX1rtrDDzUBGvL97FtWmd+dlFvThUVFr1QfLcl9s4VlpBRk4h6XX066sq\nd85IZ9xf5tc5Igngg1V7Sd95iKv++Q0PzV4PuOsXAKGuEAYntyL9hCuSK0f2DEluxSUD23Pf+F5c\nl9aZL34+ln/ePJTYyFD+MW8LK3cdorRcGdm1ddVzk2IieP7mYWQXFPPjmSuqvfdfbs5mUOdWxEeF\n0yLcxY8v6MHSjFymL9xR53EIFCrkAAAZv0lEQVTszi3i2hcWMf7p+Uz46wKmvLy0RkI60Tfbctif\nf4zE6HCWZlSv1SzcmuPzZKbyWoFXFmZQ4MeH9NbsAjrHt6zqqwcIc4Uwrm9bvth4oOr3sHznIV5e\nuIMZi3Zw44tLmPTsQqLCXcy8cyR3ndeNbdmFfO10mT47byuJ0RF89cvzmXx2CjMW7+S2l5dytMS/\nD/eVu93v69DkeMJcIdw0ogvzN2fzwao9Vb+LjfvyufhvC7jwqfl+/795OlxUwp7DR+nnUb/wxxVD\nO1JWoXzozA6xI6eQS5/5mvSdh/jT1QN58upBXD2sExeltqWbR89BdEQoT183mBCBSwd2ID4qnC4J\nLWkZ7qo2IrIhBTRhqOoCwNfpziRghrotBlqJSHvgO8BcVc1V1UPAXHwnnqC7ZEB7cgpqH22Td7SU\nH/57Bbe9vIx731zJjpxCJk9fSv6xMp64agA5BcXcOSO9xpmPqrLn8NFqy7dlF5Kc0JKI0OpNYhGh\nf8dYEqLCmTiwPQBje7ubx/M21d0tNWvFHhKiwqtGeMS1COPBialVXTkAvxjfmwn92nGwoIS/XT+Y\nv14/hJlTR9ImNpLJ05cxbcE2lu3I5S+fbQaBey/syTk9EunYqgX/WbabPYeP8saSXVw2qANR4S7e\nrqO4/uXmbL7aksPBwhJ+8U7dLZI1mYc5KyWekd1b89GafUSFuxjQ8fjZ4LAuCWzan1/tQ3HFrkN0\nim9Bm5hIAH50QU+euHog7eIiiYkM47bRXfl0/QFmLNpJiEBaSvXzlkGdW/G7y/uzcOtBfvfRRlSV\ngwXFrMk8zNhex7snbh7RhQn92vG7jzb47DJSVX793lrW78mjW2I0yQkt+WbbQb7/2vKqVlpt3lme\nSVyLMH44tgdZR4rZefB4LeaR2ev5v3fXsDu3qMbzMg8VsSYzj0sGtCf/WBmvLz5+Vq5a+9n4tqyC\nat1RlS5KbcuhotKqE4GXv9lBXIswlv3mQp6/eRhTRqcwc+pIOie05JKB7UmMDueVhTtYvfswX23J\n4c5zu9I6OoJHLuvHMzcMIX3nIX70xgq/LrhcsfMwidHhdE5oAcCtZ3ehf8dY7n1zFbe8tJRn521l\n0rMLyS0sJaegmLteS6/2+9yeXVBn/aVyFGJq+5NLGH3axdK3fSzvrdzD7twibvzXYopKynjn+2dz\nTVpnn88d1iWej+89l99d3h+AkBChd7uYoLUwQoPyqsd1BDw/NTKdZd6W1yAiU3G3TkhOTg5MlH4Y\n27sNLcNdfLh2H6N6JFZb93b6bo6WlnPD8GTeTt/N7NV7CQ0RXp0ynFE9EolrEc4PXl/O5OlL6d4m\nmqLiMvYePsbGffkcKS7jyiEdeeq6wYAzpDap5j8rwO8uH0BRSVnVmV/b2EhS28cy79ssfjC2u9fY\n846WMnfjAW4cnkx4qPdziJAQ4bmbhlJcVlHVh9s2NpKZd47kjhnLqo1/v/PcrrSPc//zXpPWib9+\nvoXfvLcWgP+7uA8twlz8d81eHr60H1ERNf8MyyuUJz7ZRHJCSyaPSuGxDzcwY9FOJjvFxRMVFpex\nNbuAey7oyb3jejJ9YQbhoSGEuY4fT1qXeCoUVu06zDk93e/Ryl2HOSslodZ9AnxvdAovfbWdOev3\nM7BTHDGRYTW2uTatM5v2HWH6wgxcIcKAjnGownm9k6r97p6+bjDXTVvEPW+u5O3vn02/Wro23l+1\nh6+25PDoZf2qjvW/q/fy45kr+fWsdfz5moFVw6kr5R8rZc66/VyT1okxvdzHtTQjl5TEKDYfOMIW\n56rs5+dv4/dXDKj23E/Xu5PXL77Tm7yjpfzrqwwmj0phf94xpr6WTkZOIWldEhjTK4nkhJa4QmB7\nTiFjeiVxojG9kggPDWHuhgMkJ7Rkzrr93HGOOwlM6N+OCf3bVW0bEerixuHJPDNvKwcLS4hrEcZN\nI7tUrZ84sAOHi0p54P11/Oyt1QxPiWfj/iPkFZVy2eAOXNi3LS6P7p+Vuw4xuHN81e8mPiqcD+4+\nh9eX7ORPn37L11tzGNenDU9ePZAlGbn88PUV3D9rLd8b3ZUn5mziqy05xLcM49azU5g8KoWEqPAa\nx1d5Vp96ki0MgCuGdOAPH2/imucXUVRSxht3jqR/x5rvf236tIut8fijNXtR1Rp/C4HW5IveqjpN\nVdNUNS0pqeYfcUNpEe5yd0utq94tVVGhvLZ4J2ld4vnjlQOY/aNzGNsriWduGFKVWCb0b8cjl/Zj\n4758Plu/nxW7DlOuyqQhHTinRyIfr9tHYXEZ5RXK9pzCGvWLSl0To2p8CJ3fJ4nlOw/5nF3347X7\nKCmr4IohtebkakJCpEbBr11cJB/++FyW/eZCXpqcxgOX9OWecT2r1l+T1hkRd7/+DcM707FVC649\nqxNFJeV8tLb2SRxnrchk0/4j/OI7vfne6BTG9k7iDx9v9FrvWbcnD1UY1DmOkBDhjnO7cevZKdW2\nGdolnhZhrqor9PflHWVf3jGGJLfyerytWoZzi7OfEV29J5YHJ/bltlEpvPR1Bo/8dz0JUeEMPOED\noUW4ixdvTSOuRRhTXl5WNb1GpdzCEh77cCNDkltxs8eH56WDOnDvuJ68uyKT5+fXHJX18Zp9FJdV\ncPWwznRPiiYhKryqjvHhmn2IwPjUtrydnlljSpM56/bRp10MKYlR/PD87uQUFPPQB+uY9OxCso8U\nc9OILhwqKuGJOZu4+40VfP/fKygpq6i1WyY6IpTR3Vvz2Yb9vLZ4J6pa7ThOdNPILrhEWLX7cNWt\nAzzdPLIL943vxX9X7+XBD9bz0Zp9LMnI5a7XljPmyXlVr3GosITtOYUM7VL9fXSFCLeencL/fj6W\nV6acxYuT02gdHcF3B7Tnpxf2YtaKPUx85mvW7cnj5xf1YliXBP72xRZGP/4//u3s29P6vfm0jY0g\nMTrC6zF5M2lwR0SgoLiM124f4XeyqE1q+xjyj7lHS4K7FlpbHTEQgt3C2AN4tsk6Ocv2AGNPWP5l\ng0V1iiYObM+Hzh/1aCcZfLk5i50Hi7hvfG/AfXby8pThNZ47eVRKrWfPSzNyufaFRXy+8QBDOsdT\nUlZRY0itL+f3bsOz87Yx5eWljO/XjrG9k+jdNqbqzGTRtoP85bPN9GgTzcBOp/5HDO4+/XF929ZY\n3rFVC87rlcTi7Qe5+/wegLuvuVtSFG+n7+aaYZ34ZN1+Zi7dRdfEKIZ3TeCpuZsZ1CmOSwa0R0R4\n8uqBTPjrV9zz5irev3tUjS65NZnukUoDOnr/8I+OCOWmEcm8/M0OfnJhL9Y6o5s8u91qc8e5XVmS\ncZBLB3Xwuo2I8PClqYSIMH1hBpMGd6hRAAVoExvJq98bzo3/WsJ1Lyzi33eMoG/7WPKPlfLg++vI\nP1rKH68cUO3sGeAnF/ZkW3YBf/p0E0OSW1UV88HdHdU9KYpBneIQEYanJLB0x0FUlY/W7GVE1wQe\nnJjKF5u+5IUF23j40n4AZB05RvrOQ9zrJPezu7VmSHIr3krPpHfbGP51axrJTq0s+0gxh4pKKK9Q\nXCFCDy+t3PH92jFv1lqmf53BRalt6ZzQstbtwN06nTiwPZ9vzGLK6JRat7n7/B6M79eOmMhQ2sVG\nUl6hfL4xi+lfZ/Dg++vIzj/GEOf98/Y+JsVEVHXPVrpnXA+OlpYT7hLuGNONWKfluOXAEX774QYe\neH8dCzZn88RVA4l3Whsb9ubX2ir0R9vYSJ67cShdWkedUgvFU1Xhe28+HeIi+fWsdWzJOsLCX11A\ny/DAfqQHO2HMBn4kIm/iLnDnqeo+EfkU+INHoXs8cH+wgvRXZbfUR2v3VSWMV7/ZSdvYiGrN8ZOR\n1iWe9nGRzF61t+qP2luXVG2GJsfzi+/0ZvaqvTz+ySYe/2QTqe1juX54Z/KKSnn6882ktI7i2RuH\nBrR5+8RVA8nKL6ZNrLtWICJcM6wzT8zZxOXPLmR1Zl7V1diVV8Y+de3gqg/dNjGRPHHVQO6ckc5f\nPtvMr7/bt9r+1+zJo0NcJEkxvs/+7hzTjRmLd/LP+duICncRERpS9Q/oTWJ0BO/9cHSdxygiPDix\nL4M6x5Hmo5urV9sY3rprJDe9uITrpy3m3J6JzN1wgOKyCn56Ya8aXRCV+378qoFs2JvPvW+u5ON7\nziXBub4jfechfvPdvlXv3/CuCcxZv5/5m7PZll3IbaO70jmhJZcP7sjMpbv44dgeJMVE8Nn6A6jC\nxf3bV73GY5P68981e/nxBT2rnfEnxUTU+bsFGNe3DSJQXFbBlNFd69z+d1cM4OfjS2jVsmYXUGVM\nvdrGVD0OdQkT+rdjfGpbfjVrDX//31a6tG6JK0RO6oRHRPjVxX1qLO/ZNoZXpwxn+sIMnpiziYnP\nfM27PxhFq5ZhbM0u4KLUmidE/rp4QPtTfq6nPh4jpfKOlrJ0Ry6PXzkg4MkCApwwRGQm7pZCoohk\n4h75FAagqs8DHwPfBbYCRcAUZ12uiDwGLHN29VtVbfTTYUaGuRjXty0fr93HgI5xdGzVgvmbs/nZ\nRb2q9aWfjJAQ4dJBHXh5YUbVB9vJJIyQEOHu83tw9/k92J93jLkb9vPmst089IF7FNFlgzrwhysH\n1OgOqG9tYyNp6ySLSlcN7cjTczez5/AxHr9yANekdaZClXV78sg/VlY1wqnSRaltuXlkMtMWbGdM\nz6SqOgS4C94DO3lvXXjGcW1aJ/6zbDed4lsyoGOcz7rNyRIRJg2uu2uvW1I0b911Nje9uIQFm7O5\nNq0z16R1qlakP1F0RCjP3DiEK577hp+/vZruSdG89HUG16V15nvnHP9wHu50nT324QZCBCb0c5+s\n3H1+d2atzOT2V5dxxZCOfLRmH10To+jV9vjfU/+OcafVXdImJpLhKQkUlpT57MLzPKZT+dsLCREe\nv3IggvCf9N306xBbbx+YlV2aZ6UkcOO/FnPnjHQeuKQv5RV60iOkAiE6IpTkhJYszjjIpn1HGJrc\nimvrKJ7XF2lO91BIS0vT9PT0oMawevdh7pyRTpYzbC/MJSz81QVVo3BOhXtK7a+JiQwlzBXCigcv\nOu041+3JI+9oKaO6t27wwpmnHTmFJMVE1Fr4rs3RknImPvMVBcVlzLl3DPFR4eQVlTLot5/xywm9\n+eHYHnXuY3duEef/+UvKKpSpY7rVaK00pGOl5YhQo4vNl9cW7eBBJ+HfNiqFhyamVuv+Kq9QBj/6\nGUeKyxjVvTVv3Dmyat3rS3by8sIdVfWT75/XvdYz7dORV1SKol5bDfWpokL55/xtdE+KYkL/+jmD\n9/TFxgPcMSOdxOgIso8UM/8XY+nS2v8u4UC567V0Pl1/AFeI8N8fnXNa3VwislxV0/zZNthdUs3O\noM6tWPLrcWzLLmRJxkFaR4WfVrIA6Nchlm6JUWzPKeSslPq5HOV0ziLrU0riyf3ztQh38bfrh3DF\ncwt58tNN/PHKgazZ476WYqCP+oWnzgktuXxIR95ZnsmQzv49J1A8r2Xw180ju5B56CjxUeHcNaZb\njYTvChGGpcTz5bfZXDKw+ofoTSO6cNOILmTkFLJk+8FT7ir1Ja5lzZFkgVLZgg6UcX3b8uuL+/L7\njzcSExFK53jvNZmG1Ld9LJ+uP8Bto1JOuyZyMixhBICI0KNNdK1j1U91f5cO6sDfvthyUt1RzVX/\njnHcMjKFV77J4LZRXY8XvE+iD/unF/WiQpVzaxke2tiJCPfX0Soa0zOJJdtz+U6/2hNC18Soaldq\nG+/uOLcrWUeOUVqutQ5kCIZLBrQnI6eQn1zYs+6N65F1STUR27ILuPCp+Tw0MdWvYmJzd7iohDFP\nzmNwcjwtwkLYfKCAefeNDXZYjUZ5hZJbWOJXodqc2axLqhnqnhTNB3ePrjZi5EzWqmU494zrye8+\n2kh4aEhVYde4uULEkoWpd03+wr0zycBOrU6pz7u5uuXsLiQntKSkrOK0ryExxtTNEoZpsiJCXfz6\nu30RodqFbMaYwLAuKdOkTejfjhUPXFR1Na4xJnCshWGaPEsWxjQMSxjGGGP8YgnDGGOMXyxhGGOM\n8YslDGOMMX6xhGGMMcYvljCMMcb4xRKGMcYYvzSryQdFJBvYeYpPTwRy6jGcxqY5H19zPjaw42vq\nGvvxdVFVv6ZtblYJ43SISLq/MzY2Rc35+JrzsYEdX1PXnI7PuqSMMcb4xRKGMcYYv1jCOG5asAMI\nsOZ8fM352MCOr6lrNsdnNQxjjDF+sRaGMcYYv1jCMMYY45czPmGIyAQR+VZEtorIr4Idz+kSkc4i\nMk9ENojIehG511meICJzRWSL8z0+2LGeDhFxichKEfnQedxVRJY47+N/RKTJ3iRDRFqJyDsisklE\nNorI2c3p/RORnzp/m+tEZKaIRDbl909EpotIlois81hW6/slbn93jnONiAwNXuQn74xOGCLiAp4F\nLgZSgRtEJDW4UZ22MuDnqpoKjATudo7pV8AXqtoT+MJ53JTdC2z0ePwE8LSq9gAOAbcHJar68Tdg\njqr2AQbhPs5m8f6JSEfgHiBNVfsDLuB6mvb79wow4YRl3t6vi4GeztdU4J8NFGO9OKMTBjAc2Kqq\n21W1BHgTmBTkmE6Lqu5T1RXOz0dwf9h0xH1crzqbvQpcHpwIT5+IdAIuAV50HgtwAfCOs0mTPT4R\niQPGAC8BqGqJqh6mGb1/uG8N3UJEQoGWwD6a8PunqguA3BMWe3u/JgEz1G0x0EpE2jdMpKfvTE8Y\nHYHdHo8znWXNgoikAEOAJUBbVd3nrNoPtA1SWPXhr8AvgQrncWvgsKqWOY+b8vvYFcgGXna63F4U\nkSiayfunqnuAPwO7cCeKPGA5zef9q+Tt/WrSnzlnesJotkQkGngX+Imq5nuuU/dY6iY5nlpEJgJZ\nqro82LEESCgwFPinqg4BCjmh+6mJv3/xuM+yuwIdgChqduc0K035/TrRmZ4w9gCdPR53cpY1aSIS\nhjtZvK6qs5zFByqbvs73rGDFd5pGA5eJyA7cXYgX4O7zb+V0cUDTfh8zgUxVXeI8fgd3Amku79+F\nQIaqZqtqKTAL93vaXN6/St7eryb9mXOmJ4xlQE9nhEY47uLb7CDHdFqc/vyXgI2q+pTHqtnAZOfn\nycAHDR1bfVDV+1W1k6qm4H6//qeqNwHzgKudzZry8e0HdotIb2fROGADzeT9w90VNVJEWjp/q5XH\n1yzePw/e3q/ZwK3OaKmRQJ5H11Wjd8Zf6S0i38XdJ+4Cpqvq74Mc0mkRkXOAr4C1HO/j/zXuOsZb\nQDLuKeCvVdUTC3VNioiMBe5T1Yki0g13iyMBWAncrKrFwYzvVInIYNwF/XBgOzAF98lds3j/RORR\n4DrcI/pWAnfg7sdvku+fiMwExuKexvwA8DDwPrW8X06S/AfubrgiYIqqpgcj7lNxxicMY4wx/jnT\nu6SMMcb4yRKGMcYYv1jCMMYY4xdLGMYYY/xiCcMYY4xfLGEYcwYTkd4i0mTmbTLBZQnDNBkioiLy\nF4/H94nIIwF6rUdE5L5A7DtYROS3InKhx+MOwCPAl8GKyTQtljBMU1IMXCkiib428phiotk6lWNU\n1YdU9XOPx3tV9QZnNlxj6mQJwzQlZcA04KcnrhCRV0TkeRFZAjzp3MDmfecmNYtFZGBdOxeR34jI\nZhH5Gujtsby7iMwRkeUi8pWI9KnluUnOjXLWOzPM7qxMbCJys4gsFZFVIvKCcx8WRKRARH4vIqud\nGNt67OtdEVnmfI12lj8iIq+JyELgNRFJceJZ4XyN8ojn/0RkrbPvxz1+R1c7P49zZsNdK+4bAEU4\ny3eIyKPO/tbWdqzmzGUJwzQ1zwI3OfeNOFEnYJSq/gx4FFipqgNxT40yw9dORWQY7rmpBgPfBc7y\nWD0N+LGqDgPuA56rZRcP457Xqh/uCQOTnf32xT0NxmhVHQyUAzc5z4kCFqvqIGABcKez/G+4byZ0\nFnAVzn0/HKnAhap6A+4J7S5S1aHOa/zdec2Lcc8IO8LZ95MnHGsk7pv+XKeqA3DPkPsDj01ynH3+\n0zleYwD3H4oxTYaq5ovIDNx3bTt6wuq3VbXc+fkc3B+2qOr/RKS1iMSeONW7h3OB91S1CEBEZjvf\no4FRwNvuaYAAiKjl+ecAVzivN0dEDjnLxwHDgGXO81twfObSEuBD5+flwEXOzxcCqR6vF+vEATBb\nVSuPOwz4hzP3VDnQy+P5L1ceSy1zTvXGPWPsZufxq8DduOdUA/cMspUxXVnLsZozlCUM0xT9FVgB\nvHzC8sIAvFYI7pv7DD7F5wvwqqreX8u6Uj0+mVs5x/8fQ4CRqnqs2o7cCcTzGH+Ke7K7Qc5zqm1/\nGion/fOMyRjrkjJNj3PG/Ba+7/v8FU7XjzOrbY6P1gW4u4QuF5EWIhIDXOq8Vj6QISLXOPsSERlU\ny/MXAtc624wH4p3lXwBXi0gbZ12CiHSp4xA/A35c+cBpQdQmDtinqhXALbhnXAaYC0wRkZaVr3nC\n874FUkSkh/P4FmB+HTEZYwnDNFl/wT2dtDePAMNEZA3wOM69CUQkTURePHFj5z7o/wFWA5/gvldK\npZuA20VkNbCe2u/7/igwXkTWAdfgvi3nEVXdADwAfObEMheo6x7O9wBpTsF+A/B9L9s9B0x24uqD\n0/pQ1Tm477uQLiKrOKEO4bRcpuDuZqucBv/5OmIyxqY3N6Y+OKOMylW1TETOxn2L1VPtxjKmUbL+\nSWPqRzLwloiE4C5m31nH9sY0OdbCMMYY4xerYRhjjPGLJQxjjDF+sYRhjDHGL5YwjDHG+MUShjHG\nGL/8PyRm+VhdDTjgAAAAAElFTkSuQmCC\n",
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ]
          },
          "metadata": {
            "tags": []
          }
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "Fg3VcvTvyxes",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        ""
      ],
      "execution_count": 0,
      "outputs": []
    }
  ]
}