show_bug.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "5f58fad0-b414-4dbb-a7b3-5413a883c1db",
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib widget\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "import numpy.fft as npf"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "dd8f7141-ede6-4d31-abfa-eb6396abae14",
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_data(fn, Nt, Ny, Nx):\n",
    "    data = np.fromfile(fn)\n",
    "    return (data[::2] + 1j*data[1::2]).reshape(Nt, Ny, Nx)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "4cdc6ab8-5bda-453d-a991-8b74ac79033b",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "1b8caa9f68104c178178cace2f76f6af",
       "version_major": 2,
       "version_minor": 0
      },
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAABYXElEQVR4nO3dbWwU57k//u/MPtkmywbHXT8EcH2qOCEyyTmYBEweSlG7wacmEF4UidQi+qXoUDCSa6oKxIvQqAd0KiWNlJSoD2kaBVJ4k/SkCseNUzAJxU6Ig0WpCUkq+NskXgyOWdvY7NNc/xezO97xAyTA7nh3vh9ptN6Ze2bve5jLXB5f91gREQERERER2YZqdQeIiIiIKLOYABIRERHZDBNAIiIiIpthAkhERERkM0wAiYiIiGyGCSARERGRzTABJCIiIrIZJoBERERENsMEkIiIiMhmmAASERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGSaAN2D37t2oqKhAXl4eqqur8d5771ndJSIiIqJrYgJ4nfbv34/GxkZs374dx48fx0MPPYTa2lp0d3db3TUiIiKiq2ICeJ2effZZPPnkk/jRj36EefPm4bnnnsOcOXPw4osvWt01oozj3XCiMYwHygZOqzuQjSKRCDo6OrB161bT+kAggKNHj066TzgcRjgcNt5rmoYvv/wSt912GxRFSWt/s5GIYGhoCGVlZVBV/pwynSXvhu/evRsPPPAAfvOb36C2thZdXV2YO3fuNffXNA1ffPEFvF4vY2ESjIXswnhIL8bDTST0tX3++ecCQP7+97+b1v/3f/+3VFZWTrrPU089JQC4fM2lp6cnE/+kdAPuv/9+2bBhg2ndXXfdJVu3bv1K+/f09Fh+nWXDwljIDowHxkO24B3AGzD+pzMRmfIntm3btqGpqcl4HwqFMHfuXPx/H30TM2+Z+qeYxyrn35zOZpkYojiCA/B6vVZ3ha7ieu6Gj5f8N2YsTI6xkD0YD+nHeLh5mABeh6KiIjgcDgSDQdP6vr4+FBcXT7qPx+OBx+OZsH7mLSpmeqcOcqfiurHOZivRX/grkOnt4sWLiMfjE6774uLiCfGRNL4cYmhoCABjYUqMhazBeMgAxsNNw1+gXwe3243q6mq0tLSY1re0tGDJkiUW9YrIOl/nbviuXbvg8/mMZc6cOZnoIlHGMB4oGzABvE5NTU34/e9/jz/84Q84deoUfvKTn6C7uxsbNmywumtEGXM9d8O3bduGUChkLD09PZnoKlHaMR4om/BXwNdpzZo16O/vx9NPP43e3l5UVVXhwIEDKC8vt7prRBmTejf8scceM9a3tLRg5cqVk+4zVTkEUbZjPFA2YQJ4AzZu3IiNGzda3Q0iSzU1NaG+vh4LFy5ETU0Nfvvb3/JuONkW44GyBRNAIrohvBtONIbxQNmCCSAR3TDeDScaw3igbMAE0GKPVc6373R+ohSMBaIxjAdKN84CJiIiIrIZJoBERERENsMEkIiIiMhmmAASERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGSaARERERDbDBJCIiIjIZpgAEhEREdkME0AiIiIim2ECSERERGQzTACJiIiIbIYJIBEREZHNMAEkIiIishkmgEREREQ2wwSQiIiIyGaYABIRERHZDBNAIiIiIpthAkhERERkM0wAiYiIiGyGCSARERGRzTABJCIiIrIZJoBERERENsMEkIiIiMhmmAASERER2QwTQCIiIiKbYQJIREREZDM5lQDu2rUL9913H7xeL/x+P1atWoXTp0+b2jzxxBNQFMW0LF682NQmHA5j8+bNKCoqwowZM/Doo4/i3LlzpjYDAwOor6+Hz+eDz+dDfX09Ll26lO4hEhEREd2wnEoADx8+jE2bNqG9vR0tLS2IxWIIBAK4fPmyqd3y5cvR29trLAcOHDBtb2xsxBtvvIF9+/bhyJEjGB4eRl1dHeLxuNFm7dq16OzsRHNzM5qbm9HZ2Yn6+vqMjJOIiIjoRjit7sDN1NzcbHr/8ssvw+/3o6OjAw8//LCx3uPxoKSkZNJjhEIhvPTSS3j11Vfx3e9+FwCwZ88ezJkzB++88w4eeeQRnDp1Cs3NzWhvb8eiRYsAAL/73e9QU1OD06dP484770zTCImIiIhuXE7dARwvFAoBAAoLC03rW1tb4ff7UVlZifXr16Ovr8/Y1tHRgWg0ikAgYKwrKytDVVUVjh49CgBoa2uDz+czkj8AWLx4MXw+n9FmvHA4jMHBQdNCNN29++67WLFiBcrKyqAoCv785z+btosIduzYgbKyMuTn52Pp0qX45z//aU1nidKM8UC5JGcTQBFBU1MTHnzwQVRVVRnra2trsXfvXhw8eBDPPPMMjh07hmXLliEcDgMAgsEg3G43Zs2aZTpecXExgsGg0cbv90/4TL/fb7QZb9euXUa9oM/nw5w5c27WUInS5vLly7j33nvxwgsvTLr9l7/8JZ599lm88MILOHbsGEpKSvC9730PQ0NDGe4pUfoxHiiX5NSvgFM1NDTgxIkTOHLkiGn9mjVrjK+rqqqwcOFClJeX46233sLq1aunPJ6IQFEU433q11O1SbVt2zY0NTUZ7wcHB5kE0rRXW1uL2traSbeJCJ577jls377diJ1XXnkFxcXFeO211/Bf//VfmewqUdoxHiiX5OQdwM2bN+PNN9/EoUOHMHv27Ku2LS0tRXl5OT799FMAQElJCSKRCAYGBkzt+vr6UFxcbLQ5f/78hGNduHDBaDOex+PBzJkzTQtRNjtz5gyCwaCpXMLj8eDb3/72lKUQRLmK8UDZJqcSQBFBQ0MDXn/9dRw8eBAVFRXX3Ke/vx89PT0oLS0FAFRXV8PlcqGlpcVo09vbi5MnT2LJkiUAgJqaGoRCIXzwwQdGm/fffx+hUMhoQ5TrkuUO43/oSS2XmAzrYSkXMR4o2+RUArhp0ybs2bMHr732GrxeL4LBIILBIEZHRwEAw8PD+OlPf4q2tjacPXsWra2tWLFiBYqKivDYY48BAHw+H5588kls2bIFf/vb33D8+HH88Ic/xPz5841ZwfPmzcPy5cuxfv16tLe3o729HevXr0ddXR1nAJPtjC97uFopBMB6WMptjAfKFjmVAL744osIhUJYunQpSktLjWX//v0AAIfDgX/84x9YuXIlKisrsW7dOlRWVqKtrQ1er9c4zq9+9SusWrUKP/jBD/DAAw+goKAAf/nLX+BwOIw2e/fuxfz58xEIBBAIBHDPPffg1VdfzfiYiaySfJTS+LsbqeUSk9m2bRtCoZCx9PT0pLWfRJnAeKBsk1OTQETkqtvz8/Px17/+9ZrHycvLw/PPP4/nn39+yjaFhYXYs2fP1+4jUa6oqKhASUkJWlpa8B//8R8AgEgkgsOHD+N//ud/ptzP4/HA4/FkqptEGcF4oGyTUwkgEd1cw8PD+Oyzz4z3Z86cQWdnJwoLCzF37lw0NjZi586duOOOO3DHHXdg586dKCgowNq1ay3sNVF6MB4olzABJKIpffjhh/jOd75jvE8+ymjdunX44x//iJ/97GcYHR3Fxo0bMTAwgEWLFuHtt982lVQQ5QrGA+USRa71e1NKi8HBQfh8PizFSjgVl9XdmXZiEkUr/hehUIiPzMlxjIWrYyzYC+Ph6hgPN09OTQIhIiIiomtjAkhERERkM0wAiYiIiGyGCSARERGRzTABJCIiIrIZJoBERERENsMEkIiIiMhmmAASERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGSaARERERDbDBJCIiIjIZpgAEhEREdkME0AiIiIim2ECSERERGQzTqs7YFciAgCIIQqIxZ2ZhmKIAhg7T5S7GAtXx1iwF8bD1TEebh4mgBbp7+8HABzBAYt7Mr0NDQ3B5/NZ3Q1KI8bCV8NYsAfGw1fDeLhxTAAtUlhYCADo7u7mRZxicHAQc+bMQXd3NxRFQVlZmdVdojRjLEwtGQ9dXV2MBZvItnhIXqM9PT2YOXNm2j9PRDA0NMR4uAmYAFpEVfXyS5/Pl5GgyTY8L/bBWLi222+/3ThPlNuyNR5mzpyZsf5mQ2KcDfgdhYiIiMhmmAASERER2QwTQIt4PB489dRT8Hg8VndlWuF5sR/+m0+N58Z+su3fPNv6S2MU4VxqIiIiIlvhHUAiIiIim2ECSERERGQzTACJiIiIbIYJoAV2796NiooK5OXlobq6Gu+9957VXUqrXbt24b777oPX64Xf78eqVatw+vRpU5snnngCiqKYlsWLF5vahMNhbN68GUVFRZgxYwYeffRRnDt3LpNDoTSwUzwwFuhqrIiFTF6TAwMDqK+vh8/ng8/nQ319PS5dupTuIdIUmABm2P79+9HY2Ijt27fj+PHjeOihh1BbW4vu7m6ru5Y2hw8fxqZNm9De3o6WlhbEYjEEAgFcvnzZ1G758uXo7e01lgMHzH8KqbGxEW+88Qb27duHI0eOYHh4GHV1dYjH45kcDt1EdosHxgJNxapYyOQ1uXbtWnR2dqK5uRnNzc3o7OxEfX19WsdHVyGUUffff79s2LDBtO6uu+6SrVu3WtSjzOvr6xMAcvjwYWPdunXrZOXKlVPuc+nSJXG5XLJv3z5j3eeffy6qqkpzc3M6u0tpZPd4YCxQ0nSJhXRdk11dXQJA2tvbjTZtbW0CQD7++OObPxC6Jt4BzKBIJIKOjg4EAgHT+kAggKNHj1rUq8wLhUIAxv7mZVJrayv8fj8qKyuxfv169PX1Gds6OjoQjUZN566srAxVVVW2One5hPHAWCDddIqFdF2TbW1t8Pl8WLRokdFm8eLF8Pl8vG4twgQwgy5evIh4PI7i4mLT+uLiYgSDQYt6lVkigqamJjz44IOoqqoy1tfW1mLv3r04ePAgnnnmGRw7dgzLli1DOBwGAASDQbjdbsyaNct0PDudu1xj93hgLFDSdImFdF6TwWAQfr9/wmf6/X5etxZxWt0BO1IUxfReRCasy1UNDQ04ceIEjhw5Ylq/Zs0a4+uqqiosXLgQ5eXleOutt7B69eopj2enc5er7BoPjAUaz+pYSPc1OdlYeN1ah3cAM6ioqAgOh2PCTzt9fX0TfvLLRZs3b8abb76JQ4cOYfbs2VdtW1paivLycnz66acAgJKSEkQiEQwMDJja2eXc5SI7xwNjgVJNh1hI9zVZUlKC8+fPTzjWhQsXeN1ahAlgBrndblRXV6OlpcW0vqWlBUuWLLGoV+knImhoaMDrr7+OgwcPoqKi4pr79Pf3o6enB6WlpQCA6upquFwu07nr7e3FyZMnc/rc5TI7xgNjgSZjZSxk6pqsqalBKBTCBx98YLR5//33EQqFeN1axZKpJza2b98+cblc8tJLL0lXV5c0NjbKjBkz5OzZs1Z3LW1+/OMfi8/nk9bWVunt7TWWkZEREREZGhqSLVu2yNGjR+XMmTNy6NAhqampkdtvv10GBweN42zYsEFmz54t77zzjnz00UeybNkyuffeeyUWi1k1NLpBdosHxgJNxapYyOQ1uXz5crnnnnukra1N2traZP78+VJXV5fW8dHUmABa4Ne//rWUl5eL2+2WBQsWmKbb5yIAky4vv/yyiIiMjIxIIBCQb3zjG+JyuWTu3Lmybt066e7uNh1ndHRUGhoapLCwUPLz86Wurm5CG8o+dooHxgJdjRWxkMlrsr+/Xx5//HHxer3i9Xrl8ccfl4GBgbSPkSaniIhk+q4jEREREVmHNYBERERENsMEkIiIiMhmmAASERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGSaARERERDbDBJCIiIjIZpgAEhEREdkME0AiIiIim2ECSERERGQzTACJiIiIbIYJIBEREZHNMAEkIiIishkmgEREREQ2wwSQiIiIyGaYAN6A3bt3o6KiAnl5eaiursZ7771ndZeIiIiIrokJ4HXav38/GhsbsX37dhw/fhwPPfQQamtr0d3dbXXXiIiIiK6KCeB1evbZZ/Hkk0/iRz/6EebNm4fnnnsOc+bMwYsvvmh114gyjnfDicYwHigbMAG8DpFIBB0dHQgEAqb1gUAAR48etahXRNbg3XCiMYwHyhaKiIjVncg2X3zxBW6//Xb8/e9/x5IlS4z1O3fuxCuvvILTp09P2CccDiMcDhvvNU3Dl19+idtuuw2KomSk39lERDA0NISysjKoKn9Omc4WLVqEBQsWmO5+z5s3D6tWrcKuXbuuub+mafjiiy/g9XoZC5NgLGQXxkN6MR5uHqfVHchm44NTRKYM2F27duHnP/95JrqVU3p6ejB79myru0FTSN4N37p1q2n917kb/sUXX2DOnDnp6F5OYSxMf4yHzGE83DgmgNehqKgIDocDwWDQtL6vrw/FxcWT7rNt2zY0NTUZ70OhEObOnYsH8Z9wwpXW/majGKI4ggPwer1Wd4Wu4uLFi4jH4xOu++Li4gnxkTT+bnjylxCMhckxFrIH4yH9GA83DxPA6+B2u1FdXY2WlhY89thjxvqWlhasXLly0n08Hg88Hs+E9U644FQY5BMkChP4K5DscDPuhjMWpsBYyDqMhzRiPNw0/AX6dWpqasLvf/97/OEPf8CpU6fwk5/8BN3d3diwYYPVXSPKmOu9Gx4KhYylp6cnE10lSjvGA2UT3gG8TmvWrEF/fz+efvpp9Pb2oqqqCgcOHEB5ebnVXSPKmJt5NxzJn+gVBRAZe5+UXDf+dfw21QFocf1VNEBJ/JwrGhSHAxKLQXG5IbEoFKcLUBUgHodaUID44CAcd1ci4r8FzuEIlJOfAS4XEI1Ci0ShuJyQcHhsf4cDoqXMo0t+XvLztfjkffy6Y+RUvaxwU+NhGokvXQBH60dWd4NuMiaAN2Djxo3YuHGj1d0gslRTUxPq6+uxcOFC1NTU4Le//S3vhpNtMR4oWzABJKIbwrvhRGMYD5QtmAAS0Q3j3XCiMYwHygacBEJE00OiVk9xOPS6N0U1luQ6xeGYos3EbYrDMba4nPpx3G69jdulr3e7oLjdUJxOKL6ZAIDhO25F34I8DMzzQplRAHVGAZT8fCguJ9T8PH3/PI++v9M58TNczq/Vx9RxprZJvjdqGImsoCjoW5BndS8oDfidhYiIiMhmmAASERER2QwTQCIiIiKbYQJIREREZDNMAIloWjBNfki8V1QFipoyESL1VVHN240JFsmJI4l1qqJP1ki8AtAnaiQmbigeN+ByQW7JBwAMzXZi6O4IBv9NgZKXBxTk65M+nE4g8cBexe3SP9PlguJQE0uivxP64Rib7JEcY8pYk2OYbPz6W/7JK7JIYhLS0N0Rq3tCacAEkIiIiMhmmAASERER2QwTQCIiIiKbYQJIRNNDstYtWdOnpjxEOfG14kjWx6XUAqY+XBnQa/tURX9NeSC0XrOX+ONHiQc4w+mEkqjj0wr0h0SPlAjmfesLXJkdBTxuSL4H8LgTD3t2GfsnawvhcIwtic/Xx6HX9yG5AFCUlNcpxjjhPPBB0GSF5LWqKpj3rS8s7gylA7+zEBEREdkME0AiIiIim2ECSERERGQzTACJiIiIbIYJIBFNC4mac/NECZjXQVWNVyV1coWqjB1ATXmQssOhr088EBpqYsKI0wkoiUkciYkcWp4+wSNyWxwBfxdmFQ9CXE5oeU5I4sHRyUkkitOp9ynxAGioqvkB1MDY5BRF0RdVMfXfGFdiHOPHnTp+IisF/F1Wd4HSgAkgERERkc0wASQiIiKyGSaARERERDbDBJCIppfUOjnjQcqquc4vIVlfl7pfshYvuV/yOPpDoVPq89TEQ6QdKqA6oLn1berMKJYWnMY3b/0ScDogLgfg1GsIxZHaN9WoLVSSdYaA8eDmsYc9p3ybHSt0NNcEjh+36XywFpAyzLhO9VhaWnDa2v5QWjABJCIiIrIZJoBERERENsMEkIiIiMhmmAASERER2QwTQBrDB8+SlVImPyjXmviQeq2OnzQxoalimnwx4TUxSURz6Os8eVFUuhSU5A9BHImJH6o61hbQJ44kJ6AkJ5wAY5NWrjI20/vUSSzJcadMarnmeSDKgEoXr8NcxASQzJgEEhER5TwmgKRL3skgIiKinMf/8Umn6H/TVEk+y4yIiIhyltPqDtidWpAPjMagzpgBCYehuN2QeByA/sfsJRKB4vFAu3x5rI3TCagqJBqDmp+H+OAgHLf6oF0ehXrLDEgsptcVuV2Q4ctQZ92KWG8QzjmzIaFBKL6ZgMsJuTyCK/fOhevtD9H7eiVO3P8nbOldgK7/Nw8xrwfOS1egBC8Ahbci/sm/4PzmXGgX+qF6b4GEI4BogNMJbWgY6i0zEO//0uiH4tIvLYlEoOTnQxsagur1QkZHoXg8QGKMUFXzGAsK9PeqCoSt+lchS11PGYKimmrpkusAbcrPMNUGAkbtntMZR4HqRr4aARQFogCiKLhmr8Z//lepBTR2ZekFTV8FqtvqLlAaMAG0mDYyClVxQbt8GQAgsZixTZKviXWTtYlHI/rrpZD+OhCZ+Bm9QQBArOecvmJw0NjmersPAFC66hQewb9D/w/zFEz/dV7s1/c/223qR6p4OGzqh0TH+iFDQ/p+idfU/httkmMcGdHfS3RCGyIiIro5+CtgIiIiIpthAkhERERkMzmVAO7atQv33XcfvF4v/H4/Vq1ahdOnzX/E+oknnjD+gHxyWbx4salNOBzG5s2bUVRUhBkzZuDRRx/FuXPnTG0GBgZQX18Pn88Hn8+H+vp6XLp0Kd1DJCIiIrphOZUAHj58GJs2bUJ7eztaWloQi8UQCARweVzN2vLly9Hb22ssBw4cMG1vbGzEG2+8gX379uHIkSMYHh5GXV0d4smJCwDWrl2Lzs5ONDc3o7m5GZ2dnaivr8/IOIlymsi120zYR4OM30+mmACS+AwRMX+Wpn8dizkwokUwqrkBESgCKF+lT+PbaJPso03epwl9J5pGRrSJteWU/XJqEkhzc7Pp/csvvwy/34+Ojg48/PDDxnqPx4OSkpJJjxEKhfDSSy/h1VdfxXe/+10AwJ49ezBnzhy88847eOSRR3Dq1Ck0Nzejvb0dixYtAgD87ne/Q01NDU6fPo0777wzTSMkIiIiunE5dQdwvFBIn5FaWFhoWt/a2gq/34/KykqsX78efX19xraOjg5Eo1EEAgFjXVlZGaqqqnD06FEAQFtbG3w+n5H8AcDixYvh8/mMNkRERETTVc4mgCKCpqYmPPjgg6iqqjLW19bWYu/evTh48CCeeeYZHDt2DMuWLUM48RiTYDAIt9uNWbNmmY5XXFyMYDBotPH7/RM+0+/3G23GC4fDGBwcNC1E0927776LFStWoKysDIqi4M9//rNpu4hgx44dKCsrQ35+PpYuXYp//vOf1nSWKM0YD5RLcjYBbGhowIkTJ/CnP/3JtH7NmjX4/ve/j6qqKqxYsQL/93//h08++QRvvfXWVY8nIuY/2j7Jg1vHt0m1a9cuY8KIz+fDnDlzrmNURJl1+fJl3HvvvXjhhRcm3f7LX/4Szz77LF544QUcO3YMJSUl+N73voehxDMfiXIJ44FySU4mgJs3b8abb76JQ4cOYfbs2VdtW1paivLycnz66acAgJKSEkQiEQwMDJja9fX1obi42Ghz/vz5Cce6cOGC0Wa8bdu2IRQKGUtPT8/1DI0oo2pra/GLX/wCq1evnrBNRPDcc89h+/btWL16NaqqqvDKK69gZGQEr7322tf/sJQJEjLZBArzh0+63+RNUyZ7TPaqaYAmUOP6uvAVFz6JCoKjXihxgRLX9DbJtgAQ11ImkmhjEz40+WqTP5Lvk8dIHXfyvaZd+zxQRmU0HqaRT6K8DnNRTiWAIoKGhga8/vrrOHjwICoqKq65T39/P3p6elBaWgoAqK6uhsvlQktLi9Gmt7cXJ0+exJIlSwAANTU1CIVC+OCDD4w277//PkKhkNFmPI/Hg5kzZ5oWomx25swZBINBU72sx+PBt7/97avWwrIcgnIR44GyTU4lgJs2bcKePXvw2muvwev1IhgMIhgMYnR0FAAwPDyMn/70p2hra8PZs2fR2tqKFStWoKioCI899hgAwOfz4cknn8SWLVvwt7/9DcePH8cPf/hDzJ8/35gVPG/ePCxfvhzr169He3s72tvbsX79etTV1XEGMNlGst51/F3v1HrZybAcgnIR44GyTU4lgC+++CJCoRCWLl2K0tJSY9m/fz8AwOFw4B//+AdWrlyJyspKrFu3DpWVlWhra4PX6zWO86tf/QqrVq3CD37wAzzwwAMoKCjAX/7yFzgcDqPN3r17MX/+fAQCAQQCAdxzzz149dVXMz5mIquNr3u9Wi0swHIIym2MB8oWOfUcwGs9TDU/Px9//etfr3mcvLw8PP/883j++eenbFNYWIg9e/Z87T4S5YrkszSDwaBRQgGY62Un4/F44PF4pj5wsj5OS6mtU/R6uPEPZJZJagElro3V5RnH0mvtJJ6s4Yvr2+JxQFUBLQ41om/TBl1oHbkTZy8VoiQWhxKNA7E4oCXqAY2+JfZPHFtJPig+8QDqsfrAlPq/lNpD0QRK6rbUcZvOB+uvskHa4sEKIoCi6NevKGgd4W+2clFO3QEkosypqKhASUmJqV42Eong8OHDU9bCEuUqxgNlm5y6A0hEN9fw8DA+++wz4/2ZM2fQ2dmJwsJCzJ07F42Njdi5cyfuuOMO3HHHHdi5cycKCgqwdu1aC3tNlB6MB8olTACJaEoffvghvvOd7xjvm5qaAADr1q3DH//4R/zsZz/D6OgoNm7ciIGBASxatAhvv/22qaaWKFcwHiiXMAEkoiktXbr0qrW1iqJgx44d2LFjR+Y6RWQRxgPlEtYAEtG0MDY/Qkyvpq9TJkpIykOXxz9AGQAkHtcnaYgYkzWg6RM1JBbTJ2LEYkAsBsTjUK9EAQDufgfe7rsbA+dnQonGoF6JQYnqbRCNGfuLCBDX9M/RNGPyiWmiSaLvkpj0MX6ih4gYD48eP+5rTWojypS3++62uguUBkwAiYiIiGyGCSARERGRzTABJCIiIrIZJoBEND0kH3isydiDnEUzfZ2srzPq7EQbq7VLPog5Htfr7eLxxMOf4/o20YwaPsRi+rpYDBKNQuIa1JEIAKAgqODUv8qQd84FhCNQRsNAOKIfJxo19oeWqCFM1hqmfL4+jpQ6RW2S+r4pxjjhPMi4B0MTZYKM1dee+leZxZ2hdGACSERERGQzTACJiIiIbIYJIBEREZHNMAEkIiIishn+JRAimhb0SRyqMenBmNQBQFFTJn0kX0WDaOrYe2MySMoDmUUDFFV/cHNy0gYASTzYWX/VgGgUyvAoAMB7Loa4x438CwK5cgVKPA4JRyCxGJRwWN8/Ek1MKommTEgRY/KJuR+Y0EfjVRITRfSVpjEmxz+2nSjDRABo8Ha5re4JpQHvABIRERHZDBNAIiIiIpthAkhERERkM6wBJKLpQTQAauKhzXrtkbEpDr3GLvWhz6ltJtlmPPxZEb28TjRIRH/Ys0Si+vZIFFAVvR4wNAgAuOXTS3CHboFzOAK5PAJxufRav2gMGq7o+18JT6zRS35G4mHTX7mPpnOQsk60CeeBKONE4P/oitW9oDTgHUAiIiIim2ECSERERGQzTACJiIiIbIY1gBZJ/lH4GKIAH/M1QQxRAGPniXKXEQsShR4MSqL2TRnf0LwttU3qOtH0osHEswKNn3NFgyIaRGJQRIFIFIoAEAWQOFSJIC5RSDyMWMwJxCNQJJI4ZhSaxKCIJPZL7q+Zr9Hk5xmfH5+8j19jjDHETOeJctt0/L8hHrsCkajV3QDA/xtuJiaAFunv7wcAHMEBi3syvQ0NDcHn81ndDUqjsVh4S1+R/L4+2ff38dtkkm3auNfUNrHEa3TcKwAMJl4/TixJqfXv4XH7xTDR+H5M1sevM8YExoI9TMv/G977X6t7MAHj4cYxAbRIYWEhAKC7u5sXcYrBwUHMmTMH3d3dUBQFZWVlVneJ0oyxMLVkPHR1dTEWbCLb4iF5jfb09GDmzJlp/zwRwdDQEOPhJmACaBFV1X8t5fP5MhI02YbnxT4YC9d2++23G+eJclu2xsPMmTMz1t9sSIyzAb+jEBEREdkME0AiIiIim2ECaBGPx4OnnnoKHo/H6q5MKzwv9sN/86nx3NhPtv2bZ1t/aYwinEtNREREZCu8A0hERERkM0wAiYiIiGyGCSARERGRzTABJCIiIrIZJoAW2L17NyoqKpCXl4fq6mq89957VncprXbt2oX77rsPXq8Xfr8fq1atwunTp01tnnjiCSiKYloWL15sahMOh7F582YUFRVhxowZePTRR3Hu3LlMDoXSwE7xwFigq7EiFjJ5TQ4MDKC+vh4+nw8+nw/19fW4dOlSuodIU2ACmGH79+9HY2Mjtm/fjuPHj+Ohhx5CbW0turu7re5a2hw+fBibNm1Ce3s7WlpaEIvFEAgEcPnyZVO75cuXo7e311gOHDD/LczGxka88cYb2LdvH44cOYLh4WHU1dUhHo9ncjh0E9ktHhgLNBWrYiGT1+TatWvR2dmJ5uZmNDc3o7OzE/X19WkdH12FUEbdf//9smHDBtO6u+66S7Zu3WpRjzKvr69PAMjhw4eNdevWrZOVK1dOuc+lS5fE5XLJvn37jHWff/65qKoqzc3N6ewupZHd44GxQEnTJRbSdU12dXUJAGlvbzfatLW1CQD5+OOPb/5A6Jp4BzCDIpEIOjo6EAgETOsDgQCOHj1qUa8yLxQKARj7o+dJra2t8Pv9qKysxPr169HX12ds6+joQDQaNZ27srIyVFVV2erc5RLGA2OBdNMpFtJ1Tba1tcHn82HRokVGm8WLF8Pn8/G6tQgTwAy6ePEi4vE4iouLTeuLi4sRDAYt6lVmiQiamprw4IMPoqqqylhfW1uLvXv34uDBg3jmmWdw7NgxLFu2DOFwGAAQDAbhdrsxa9Ys0/HsdO5yjd3jgbFASdMlFtJ5TQaDQfj9/gmf6ff7ed1axGl1B+xIURTTexGZsC5XNTQ04MSJEzhy5Ihp/Zo1a4yvq6qqsHDhQpSXl+Ott97C6tWrpzyenc5drrJrPDAWaDyrYyHd1+RkY+F1ax3eAcygoqIiOByOCT/t9PX1TfjJLxdt3rwZb775Jg4dOoTZs2dftW1paSnKy8vx6aefAgBKSkoQiUQwMDBgameXc5eL7BwPjAVKNR1iId3XZElJCc6fPz/hWBcuXOB1axEmgBnkdrtRXV2NlpYW0/qWlhYsWbLEol6ln4igoaEBr7/+Og4ePIiKiopr7tPf34+enh6UlpYCAKqrq+FyuUznrre3FydPnszpc5fL7BgPjAWajJWxkKlrsqamBqFQCB988IHR5v3330coFOJ1axVLpp7Y2L59+8TlcslLL70kXV1d0tjYKDNmzJCzZ89a3bW0+fGPfyw+n09aW1ult7fXWEZGRkREZGhoSLZs2SJHjx6VM2fOyKFDh6SmpkZuv/12GRwcNI6zYcMGmT17trzzzjvy0UcfybJly+Tee++VWCxm1dDoBtktHhgLNBWrYiGT1+Ty5cvlnnvukba2Nmlra5P58+dLXV1dWsdHU2MCaIFf//rXUl5eLm63WxYsWGCabp+LAEy6vPzyyyIiMjIyIoFAQL7xjW+Iy+WSuXPnyrp166S7u9t0nNHRUWloaJDCwkLJz8+Xurq6CW0o+9gpHhgLdDVWxEImr8n+/n55/PHHxev1itfrlccff1wGBgbSPkaanCIikum7jkRERERkHdYAEhEREdkME0AiIiIim2ECSERERGQzTACJiIiIbIYJIBEREZHNMAEkIiIishkmgEREREQ2wwSQiIiIyGaYABIRERHZDBNAIiIiIpthAkhERERkM0wAiYiIiGyGCSARERGRzTABJCIiIrIZJoBERERENsMEkIiIiMhmmAASERER2QwTwBuwe/duVFRUIC8vD9XV1Xjvvfes7hIRERHRNTEBvE779+9HY2Mjtm/fjuPHj+Ohhx5CbW0turu7re4aERER0VUpIiJWdyIbLVq0CAsWLMCLL75orJs3bx5WrVqFXbt2WdgzIiIioqvjHcDrEIlE0NHRgUAgYFofCARw9OhRi3pFZB2WQxCNYTxQNnBa3YFsdPHiRcTjcRQXF5vWFxcXIxgMTrpPOBxGOBw23muahi+//BK33XYbFEVJa3+zkYhgaGgIZWVlUFX+nDKdJcshdu/ejQceeAC/+c1vUFtbi66uLsydO/ea+2uahi+++AJer5exMAnGQnZhPKQX4+EmEvraPv/8cwEgR48eNa3/xS9+IXfeeeek+zz11FMCgMvXXHp6ejLxT0o34P7775cNGzaY1t11112ydevWr7R/T0+P5ddZNiyMhezAeGA8ZAveAbwORUVFcDgcE+729fX1TbgrmLRt2zY0NTUZ70OhEObOnYsH8Z9wwpXW/majGKI4ggPwer1Wd4WuIlkOsXXrVtP6q5VDjL8bLokyZMbC5BgL2YPxkH6Mh5uHCeB1cLvdqK6uRktLCx577DFjfUtLC1auXDnpPh6PBx6PZ8J6J1xwKgzyCfTvgfwVyDR3PeUQu3btws9//vMJ6xkLU2AsZA3GQwYwHm4a/gL9OjU1NeH3v/89/vCHP+DUqVP4yU9+gu7ubmzYsMHqrhFl3PhvxiIy5Tfobdu2IRQKGUtPT08mukiUMYwHyga8A3id1qxZg/7+fjz99NPo7e1FVVUVDhw4gPLycqu7RpQx11MOMdXdcCT/g1QUQGTsfVJy3fjX8dtUB6DF9VfRACXxc65oUBwOSCwGxeWGxKJQnC5AVYB4HGpBAeKDg3DcXYmI/xY4hyNQTn4GuFxANAotEoXickLC4bH9HQ6IJil9THxe8vO1+OR9/LpjFFAWuKnxMI3Ely6Ao/Ujq7tBNxnvAN6AjRs34uzZswiHw+jo6MDDDz9sdZeIMiq1HCJVS0sLlixZYlGviKzBeKBswjuARHRDmpqaUF9fj4ULF6Kmpga//e1vWQ5BtsV4oGzBBJCIbgjLIYjGMB4oWzABJKIbtnHjRmzcuPHGDqKogECvq4vHx2r3ACiqAonHx2r4JrTRzNsS9X4QRW+jKpBoDIrbrbdxu/SaQLcLUFUgGoXimwkMDmL4jltx6VtO5F/w4LYzBVBcLkg4AkUEan4e4uEwlDwPMKpBcTqBuKZ3IfkZLickHDf6MWUfk+NJrSFMaaPvpwFQWQOYZW5KPEwXioK+BXkobbW6I3SzsQaQiIiIyGaYABIRERHZDBNAIiIiIpthAkhERERkM0wAiWhaUByOxBeq8V5RFSiqMjaRIvVVUc3bE/srjsQ2R2KdqkBxOvW2Tn3em+JyAg4HFJcTiscNuFyQW/IBAEOznRi6O4LBf1Og5OUBBflQ8jz6vokH9ipul/6ZLhcUh5pYEv2d0A+Hvi6lj0gZa3IMk41ff8s/eUUWUfTYGro7YnVPKA2YABIRERHZDBNAIiIiIpthAkhERERkM0wAiWh6SNa6qWMPb07W+iW/VhzJ+riUWkCj5m6svk5RFf1VUcw1eK7Es++dTn2d0wklUcenFbgBACMlgnnf+gJXZkcBjxuS7wE8bv04Lpexf7K2EA7H2JL4fH0cen0fkgsARUl5nWKME86Dwm/TZIHktaoqmPetLyzuDKUDv7MQERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgEU0LiZpz80QJmNdBVY1XJXVyhaqMHUBNeZCyw6GvTzwQGmpiwojTCSiJSRyJiRxanj7BI3JbHAF/F2YVD0JcTmh5TkjiwdHJSSSK06n3KfEAaKiq+QHUwNjkFEXRF1Ux9d8YV2Ic48edOn4iKwX8XVZ3gdKACSARERGRzTABJCIiIrIZJoBERERENsMEkIiml9Q6OeNByqq5zi8hWV+Xul+yFi+5X/I4+kOhU+rz1MRDpB0qoDqgufVt6swolhacxjdv/RJwOiAuB+DUawjFkdo31agtVJJ1hoDx4Oaxhz2nfJsdK3Q01wSOH7fpfLAWkDLMuE71WFpacNra/lBaMAEkIiIishkmgEREREQ2wwSQiIiIyGaYABIRERHZDBNAGsMHz5KVUiY/KNea+JB6rY6fNDGhqWKafDHhNTFJRHPo6zx5UVS6FJTkD0EciYkfqjrWFtAnjiQnoCQnnABjk1auMjbT+9RJLMlxp0xqueZ5IMqAShevw1zEBJDMmAQSERHlPCaApEveySAiIqKcx//xSafof9NUST7LjIiIiHKW0+oO2J1akA+MxqDOmAEJh6G43ZB4HID+x+wlEoHi8UC7fHmsjdMJqCokGoOan4f44CAct/qgXR6FessMSCym1xW5XZDhy1Bn3YpYbxDOObMhoUEovpmAywm5PIIr986F6+0P0ft6JU7c/yds6V2Arv83DzGvB85LV6AELwCFtyL+yb/g/OZcaBf6oXpvgYQjgGiA0wltaBjqLTMQ7//S6Ifi0i8tiUSg5OdDGxqC6vVCRkeheDxAYoxQVfMYCwr096oKhK36VyFLXU8ZgqKaaumS6wBtys8w1QYCRu2e0xlHgepGvhoBFAWiAKIouGavxn/+V6kFNHZl6QVNXwWq2+ouUBowAbSYNjIKVXFBu3wZACCxmLFNkq+JdZO1iUcj+uulkP46EJn4Gb1BAECs55y+YnDQ2OZ6uw8AULrqFB7Bv0P/D/MUTP91XuzX9z/bbepHqng4bOqHRMf6IUND+n6J19T+G22SYxwZ0d9LdEIbIiIiujn4K2AiIiIim2ECSERERGQzOZUA7tq1C/fddx+8Xi/8fj9WrVqF06fNf8T6iSeeMP6AfHJZvHixqU04HMbmzZtRVFSEGTNm4NFHH8W5c+dMbQYGBlBfXw+fzwefz4f6+npcunQp3UMkIiIiumE5lQAePnwYmzZtQnt7O1paWhCLxRAIBHB5XM3a8uXL0dvbaywHDhwwbW9sbMQbb7yBffv24ciRIxgeHkZdXR3iyYkLANauXYvOzk40NzejubkZnZ2dqK+vz8g4iXKayLXbTNhHg4zfT6aYAJL4DBExf5amfx2LOTCiRTCquQERKAIoX6VP49tok+yjTd6nCX0nmkZGtIm15ZT9cmoSSHNzs+n9yy+/DL/fj46ODjz88MPGeo/Hg5KSkkmPEQqF8NJLL+HVV1/Fd7/7XQDAnj17MGfOHLzzzjt45JFHcOrUKTQ3N6O9vR2LFi0CAPzud79DTU0NTp8+jTvvvDNNIyQiIiK6cTl1B3C8UEifkVpYWGha39raCr/fj8rKSqxfvx59fX3Gto6ODkSjUQQCAWNdWVkZqqqqcPToUQBAW1sbfD6fkfwBwOLFi+Hz+Yw2RERERNNVziaAIoKmpiY8+OCDqKqqMtbX1tZi7969OHjwIJ555hkcO3YMy5YtQzjxGJNgMAi3241Zs2aZjldcXIxgMGi08fv9Ez7T7/cbbcYLh8MYHBw0LUTT3bvvvosVK1agrKwMiqLgz3/+s2m7iGDHjh0oKytDfn4+li5din/+85/WdJYozRgPlEtyNgFsaGjAiRMn8Kc//cm0fs2aNfj+97+PqqoqrFixAv/3f/+HTz75BG+99dZVjyci5j/aPsmDW8e3SbVr1y5jwojP58OcOXOuY1REmXX58mXce++9eOGFFybd/stf/hLPPvssXnjhBRw7dgwlJSX43ve+h6HEMx+JcgnjgXJJTiaAmzdvxptvvolDhw5h9uzZV21bWlqK8vJyfPrppwCAkpISRCIRDAwMmNr19fWhuLjYaHP+/PkJx7pw4YLRZrxt27YhFAoZS09Pz/UMjSijamtr8Ytf/AKrV6+esE1E8Nxzz2H79u1YvXo1qqqq8Morr2BkZASvvfba1/+wlAkSMtkECvOHT7rf5E1TJntM9qppgCZQ4/q68BUXPokKgqNeKHGBEtf0Nsm2ABDXUiaSaGMTPjT5apM/ku+Tx0gdd/K9pl37PFBGZTQeppFPorwOc1FOJYAigoaGBrz++us4ePAgKioqrrlPf38/enp6UFpaCgCorq6Gy+VCS0uL0aa3txcnT57EkiVLAAA1NTUIhUL44IMPjDbvv/8+QqGQ0WY8j8eDmTNnmhaibHbmzBkEg0FTvazH48G3v/3tq9bCshyCchHjgbJNTiWAmzZtwp49e/Daa6/B6/UiGAwiGAxidHQUADA8PIyf/vSnaGtrw9mzZ9Ha2ooVK1agqKgIjz32GADA5/PhySefxJYtW/C3v/0Nx48fxw9/+EPMnz/fmBU8b948LF++HOvXr0d7ezva29uxfv161NXVcQYw2Uay3nX8Xe/UetnJsByCchHjgbJNTiWAL774IkKhEJYuXYrS0lJj2b9/PwDA4XDgH//4B1auXInKykqsW7cOlZWVaGtrg9frNY7zq1/9CqtWrcIPfvADPPDAAygoKMBf/vIXOBwOo83evXsxf/58BAIBBAIB3HPPPXj11VczPmYiq42ve71aLSzAcgjKbYwHyhY59RzAaz1MNT8/H3/961+veZy8vDw8//zzeP7556dsU1hYiD179nztPhLliuSzNIPBoFFCAZjrZSfj8Xjg8XimPnCyPk5Lqa1T9Hq48Q9klklqASWujdXlGcfSa+0knqzhi+vb4nFAVQEtDjWib9MGXWgduRNnLxWiJBaHEo0DsTigJeoBjb4l9k8cW0k+KD7xAOqx+sCU+r+U2kPRBErqttRxm84H66+yQdriwQoigKLo168oaB3hb7ZyUU7dASSizKmoqEBJSYmpXjYSieDw4cNT1sIS5SrGA2WbnLoDSEQ31/DwMD777DPj/ZkzZ9DZ2YnCwkLMnTsXjY2N2LlzJ+644w7ccccd2LlzJwoKCrB27VoLe02UHowHyiVMAIloSh9++CG+853vGO+bmpoAAOvWrcMf//hH/OxnP8Po6Cg2btyIgYEBLFq0CG+//bapppYoVzAeKJcwASSiKS1duvSqtbWKomDHjh3YsWNH5jpFZBHGA+US1gAS0bQwNj9CTK+mr1MmSkjKQ5fHP0AZACQe1ydpiBiTNaDpEzUkFtMnYsRiQCwGxONQr0QBAO5+B97uuxsD52dCicagXolBieptEI0Z+4sIENf0z9E0Y/KJaaJJou+SmPQxfqKHiBgPjx4/7mtNaiPKlLf77ra6C5QGTACJiIiIbIYJIBEREZHNMAEkIiIishkmgEQ0PSQfeKzJ2IOcRTN9nayvM+rsRBurtUs+iDke1+vt4vHEw5/j+jbRjBo+xGL6ulgMEo1C4hrUkQgAoCCo4NS/ypB3zgWEI1BGw0A4oh8nGjX2h5aoIUzWGqZ8vj6OlDpFbZL6vinGOOE8yLgHQxNlgozV1576V5nFnaF0YAJIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGf4lECKaFvRJHKox6cGY1AFAUVMmfSRfRYNo6th7YzJIygOZRQMUVX9wc3LSBgBJPNhZf9WAaBTK8CgAwHsuhrjHjfwLArlyBUo8DglHILEYlHBY3z8STUwqiaZMSBFj8om5H5jQR+NVEhNF9JWmMSbHP7adKMNEAGjwdrmt7gmlAe8AEhEREdkME0AiIiIim2ECSERERGQzrAEkoulBNABq4qHNeu2RsSkOvcYu9aHPqW0m2WY8/FkRvbxONEhEf9izRKL69kgUUBW9HjA0CAC45dNLcIdugXM4Ark8AnG59Fq/aAwaruj7XwlPrNFLfkbiYdNfuY+mc5CyTrQJ54Eo40Tg/+iK1b2gNOAdQCIiIiKbYQJIREREZDNMAImIiIhshjWAFkn+UfgYogAf8zVBDFEAY+eJcpcRCxKFHgxKovZNGd/QvC21Teo60fSiwcSzAo2fc0WDIhpEYlBEgUgUigAQBZA4VIkgLlFIPIxYzAnEI1AkkjhmFJrEoIgk9kvur5mv0eTnGZ8fn7yPX2OMMcRM54ly23T8vyEeuwKRqNXdAMD/G24mJoAW6e/vBwAcwQGLezK9DQ0NwefzWd0NSqOxWHhLX5H8vj7Z9/fx22SSbdq419Q2scRrdNwrAAwmXj9OLEmp9e/hcfvFMNH4fkzWx68zxgTGgj1My/8b3vtfq3swAePhxjEBtEhhYSEAoLu7mxdxisHBQcyZMwfd3d1QFAVlZWVWd4nSjLEwtWQ8dHV1MRZsItviIXmN9vT0YObMmWn/PBHB0NAQ4+EmYAJoEVXVfy3l8/kyEjTZhufFPhgL13b77bcb54lyW7bGw8yZMzPW32xIjLMBv6MQERER2QwTQCIiIiKbYQJoEY/Hg6eeegoej8fqrkwrPC/2w3/zqfHc2E+2/ZtnW39pjCKcS01ERERkK7wDSERERGQzTACJiIiIbIYJIBEREZHNMAEkIiIishkmgBbYvXs3KioqkJeXh+rqarz33ntWdymtdu3ahfvuuw9erxd+vx+rVq3C6dOnTW2eeOIJKIpiWhYvXmxqEw6HsXnzZhQVFWHGjBl49NFHce7cuUwOhdLATvHAWKCrsSIWMnlNDgwMoL6+Hj6fDz6fD/X19bh06VK6h0hTYAKYYfv370djYyO2b9+O48eP46GHHkJtbS26u7ut7lraHD58GJs2bUJ7eztaWloQi8UQCARw+fJlU7vly5ejt7fXWA4cMP8tzMbGRrzxxhvYt28fjhw5guHhYdTV1SEej2dyOHQT2S0eGAs0FatiIZPX5Nq1a9HZ2Ynm5mY0Nzejs7MT9fX1aR0fXYVQRt1///2yYcMG07q77rpLtm7dalGPMq+vr08AyOHDh41169atk5UrV065z6VLl8Tlcsm+ffuMdZ9//rmoqirNzc3p7C6lkd3jgbFASdMlFtJ1TXZ1dQkAaW9vN9q0tbUJAPn4449v/kDomngHMIMikQg6OjoQCARM6wOBAI4ePWpRrzIvFAoBGPuj50mtra3w+/2orKzE+vXr0dfXZ2zr6OhANBo1nbuysjJUVVXZ6tzlEsYDY4F00ykW0nVNtrW1wefzYdGiRUabxYsXw+fz8bq1CBPADLp48SLi8TiKi4tN64uLixEMBi3qVWaJCJqamvDggw+iqqrKWF9bW4u9e/fi4MGDeOaZZ3Ds2DEsW7YM4XAYABAMBuF2uzFr1izT8ex07nKN3eOBsUBJ0yUW0nlNBoNB+P3+CZ/p9/t53VrEaXUH7EhRFNN7EZmwLlc1NDTgxIkTOHLkiGn9mjVrjK+rqqqwcOFClJeX46233sLq1aunPJ6dzl2usms8MBZoPKtjId3X5GRj4XVrHd4BzKCioiI4HI4JP+309fVN+MkvF23evBlvvvkmDh06hNmzZ1+1bWlpKcrLy/Hpp58CAEpKShCJRDAwMGBqZ5dzl4vsHA+MBUo1HWIh3ddkSUkJzp8/P+FYFy5c4HVrESaAGeR2u1FdXY2WlhbT+paWFixZssSiXqWfiKChoQGvv/46Dh48iIqKimvu09/fj56eHpSWlgIAqqur4XK5TOeut7cXJ0+ezOlzl8vsGA+MBZqMlbGQqWuypqYGoVAIH3zwgdHm/fffRygU4nVrFUumntjYvn37xOVyyUsvvSRdXV3S2NgoM2bMkLNnz1rdtbT58Y9/LD6fT1pbW6W3t9dYRkZGRERkaGhItmzZIkePHpUzZ87IoUOHpKamRm6//XYZHBw0jrNhwwaZPXu2vPPOO/LRRx/JsmXL5N5775VYLGbV0OgG2S0eGAs0FatiIZPX5PLly+Wee+6RtrY2aWtrk/nz50tdXV1ax0dTYwJogV//+tdSXl4ubrdbFixYYJpun4sATLq8/PLLIiIyMjIigUBAvvGNb4jL5ZK5c+fKunXrpLu723Sc0dFRaWhokMLCQsnPz5e6uroJbSj72CkeGAt0NVbEQiavyf7+fnn88cfF6/WK1+uVxx9/XAYGBtI+RpqcIiKS6buORERERGQd1gASERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGSaARERERDbDBJCIiIjIZpgAEhEREdkME0AiIiIim2ECSERERGQzTACJiIiIbIYJIBEREZHNMAEkIiIishkmgEREREQ2wwSQiIiIyGaYABIRERHZDBNAIiIiIpthAngDdu/ejYqKCuTl5aG6uhrvvfee1V0iIiIiuiYmgNdp//79aGxsxPbt23H8+HE89NBDqK2tRXd3t9VdIyIiIroqJoDX6dlnn8WTTz6JH/3oR5g3bx6ee+45zJkzBy+++KLVXSPKON4NJxrDeKBs4LS6A9koEomgo6MDW7duNa0PBAI4evTopPuEw2GEw2HjvaZp+PLLL3HbbbdBUZS09jcbiQiGhoZQVlYGVeXPKdNZ8m747t278cADD+A3v/kNamtr0dXVhblz515zf03T8MUXX8Dr9TIWJsFYyC6Mh/RiPNxEQl/b559/LgDk73//u2n9f//3f0tlZeWk+zz11FMCgMvXXHp6ejLxT0o34P7775cNGzaY1t11112ydevWr7R/T0+P5ddZNiyMhezAeGA8ZAveAbwB4386E5Epf2Lbtm0bmpqajPehUAhz587Fg/hPOOFKaz+zUQxRHMEBeL1eq7tCV3Ez7oaLCAAwFqbAWMgejIf0YzzcPEwAr0NRUREcDgeCwaBpfV9fH4qLiyfdx+PxwOPxTFjvhAtOhUE+gf49kL8CmeYuXryIeDw+4bovLi6eEB9Ju3btws9//vMJ6xkLU2AsZA3GQwYwHm4a/gL9OrjdblRXV6OlpcW0vqWlBUuWLLGoV0TW+bp3w0OhkLH09PRkootEGcN4oGzAO4DXqampCfX19Vi4cCFqamrw29/+Ft3d3diwYYPVXSPKmJt5NxzJ/yAVBRAZe5+UXDf+dfw21QFocf1VNEBJ/JwrGhSHAxKLQXG5IbEoFKcLUBUgHodaUID44CAcd1ci4r8FzuEIlJOfAS4XEI1Ci0ShuJyQcHhsf4cDoklKHxOfl/x8LT55H7/uGAWUBW5qPEwj8aUL4Gj9yOpu0E3GO4DXac2aNXjuuefw9NNP49///d/x7rvv4sCBAygvL7e6a0QZw7vhRGMYD5RNeAfwBmzcuBEbN260uhtEluLdcKIxjAfKFkwAieiGrFmzBv39/Xj66afR29uLqqoq3g0n22I8ULZgAkhEN+ym3A1XVECg19XF42O1ewAUVYHE42M1fBPaaOZtiXo/iKK3URVINAbF7dbbuF16TaDbBagqEI1C8c0EBgcxfMetuPQtJ/IveHDbmQIoLhckHIEiAjU/D/FwGEqeBxjVoDidQFzTu5D8DJcTEo4b/Ziyj8nxpNYQprTR99MAqKwBzDI59dshRUHfgjyUtlrdEbrZWANIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGSaARDQtKA5H4gvVeK+oChRVGZtIkfqqqObtif0VR2KbI7FOVaA4nXpbpz7vTXE5AYcDissJxeMGXC7ILfkAgKHZTgzdHcHgvylQ8vKAgnwoeR5938QDexW3S/9MlwuKQ00sif5O6IdDX5fSR6SMNTmGycavv+WfvCKLKHpsDd0dsbonlAZMAImIiIhshgkgERERkc0wASQiIiKyGSaARDQ9JGvd1LGHNydr/ZJfK45kfVxKLaBRczdWX6eoiv6qKOYaPFfi2fdOp77O6YSSqOPTCtwAgJESwbxvfYErs6OAxw3J9wAet34cl8vYP1lbCIdjbEl8vj4Ovb4PyQWAoqS8TjHGCedB4bdpskDyWlUVzPvWFxZ3htKB31mIiIiIbIYJIBEREZHNMAEkIiIishkmgEREREQ2wwSQiKaFRM25eaIEzOugqsarkjq5QlXGDqCmPEjZ4dDXJx4IDTUxYcTpBJTEJI7ERA4tT5/gEbktjoC/C7OKByEuJ7Q8JyTx4OjkJBLF6dT7lHgANFTV/ABqYGxyiqLoi6qY+m+MKzGO8eNOHT+RlQL+Lqu7QGnABJCIiIjIZpgAEhEREdkME0AiIiIim2ECSETTS2qdnPEgZdVc55eQrK9L3S9Zi5fcL3kc/aHQKfV5auIh0g4VUB3Q3Po2dWYUSwtO45u3fgk4HRCXA3DqNYTiSO2batQWKsk6Q8B4cPPYw55Tvs2OFTqaawLHj9t0PlgLSBlmXKd6LC0tOG1tfygtmAASERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgjeGDZ8lKKZMflGtNfEi9VsdPmpjQVDFNvpjwmpgkojn0dZ68KCpdCkryhyCOxMQPVR1rC+gTR5ITUJITToCxSStXGZvpfeokluS4Uya1XPM8EGVApYvXYS5iAkhmTAKJiIhyHhNA0iXvZBAREVHO4//4pFP0v2mqJJ9lRkRERDnLaXUH7E4tyAdGY1BnzICEw1Dcbkg8DkD/Y/YSiUDxeKBdvjzWxukEVBUSjUHNz0N8cBCOW33QLo9CvWUGJBbT64rcLsjwZaizbkWsNwjnnNmQ0CAU30zA5YRcHsGVe+fC9faH6H29Eifu/xO29C5A1/+bh5jXA+elK1CCF4DCWxH/5F9wfnMutAv9UL23QMIRQDTA6YQ2NAz1lhmI939p9ENx6ZeWRCJQ8vOhDQ1B9Xoho6NQPB4gMUaoqnmMBQX6e1UFwlb9q5ClrqcMQVFNtXTJdYA25WeYagMBo3bP6YyjQHUjX40AigJRAFEUXLNX4z//q9QCGruy9IKmrwLVbXUXKA2YAFpMGxmFqrigXb4MAJBYzNgmydfEusnaxKMR/fVSSH8diEz8jN4gACDWc05fMThobHO93QcAKF11Co/g36H/h3kKpv86L/br+5/tNvUjVTwcNvVDomP9kKEhfb/Ea2r/jTbJMY6M6O8lOqENERER3Rz8FTARERGRzTABJCIiIrIZJoBERERENpNTCeCuXbtw3333wev1wu/3Y9WqVTh9+rSpzRNPPAElUQCeXBYvXmxqEw6HsXnzZhQVFWHGjBl49NFHce7cOVObgYEB1NfXw+fzwefzob6+HpcuXUr3EIlyn8i120zYR4OM30+mmACS+AwRMX+Wpn8dizkwokUwqrkBESgCKF+lT+PbaJPso03epwl9J5pGRrSJteWU/XIqATx8+DA2bdqE9vZ2tLS0IBaLIRAI4PK4SQvLly9Hb2+vsRw4cMC0vbGxEW+88Qb27duHI0eOYHh4GHV1dYgnZ64CWLt2LTo7O9Hc3Izm5mZ0dnaivr4+I+MkIiIiuhE5NQu4ubnZ9P7ll1+G3+9HR0cHHn74YWO9x+NBSUnJpMcIhUJ46aWX8Oqrr+K73/0uAGDPnj2YM2cO3nnnHTzyyCM4deoUmpub0d7ejkWLFgEAfve736GmpganT5/GnXfemaYREhEREd24nLoDOF4opD+SpLCw0LS+tbUVfr8flZWVWL9+Pfr6+oxtHR0diEajCAQCxrqysjJUVVXh6NGjAIC2tjb4fD4j+QOAxYsXw+fzGW3GC4fDGBwcNC1E0927776LFStWoKysDIqi4M9//rNpu4hgx44dKCsrQ35+PpYuXYp//vOf1nSWKM0YD5RLcjYBFBE0NTXhwQcfRFVVlbG+trYWe/fuxcGDB/HMM8/g2LFjWLZsGcKJ59gFg0G43W7MmjXLdLzi4mIEg0Gjjd/vn/CZfr/faDPerl27jHpBn8+HOXPm3KyhEqXN5cuXce+99+KFF16YdPsvf/lLPPvss3jhhRdw7NgxlJSU4Hvf+x6GEs98JMoljAfKJTn1K+BUDQ0NOHHiBI4cOWJav2bNGuPrqqoqLFy4EOXl5XjrrbewevXqKY8nIqan9U/25P7xbVJt27YNTU1NxvvBwUEmgTTt1dbWora2dtJtIoLnnnsO27dvN2LnlVdeQXFxMV577TX813/919f7sJQJEqLJ1f/yhmnyxlUmeyT6aUzimOxV0wBNoMb1deErLnwSFQRHvVDiAiWu6W2SbQEgrqVMJNH3hwP661eZ/JF8nzxG6riT7zUNokkO/5iefTIaD9PIJ1FOUspFOfmtZfPmzXjzzTdx6NAhzJ49+6ptS0tLUV5ejk8//RQAUFJSgkgkgoGBAVO7vr4+FBcXG23Onz8/4VgXLlww2ozn8Xgwc+ZM00KUzc6cOYNgMGgql/B4PPj2t789ZSkEUa5iPFC2yakEUETQ0NCA119/HQcPHkRFRcU19+nv70dPTw9KS0sBANXV1XC5XGhpaTHa9Pb24uTJk1iyZAkAoKamBqFQCB988IHR5v3330coFDLaEOW6ZLnD+B96UsslJsN6WMpFjAfKNjmVAG7atAl79uzBa6+9Bq/Xi2AwiGAwiNHRUQDA8PAwfvrTn6KtrQ1nz55Fa2srVqxYgaKiIjz22GMAAJ/PhyeffBJbtmzB3/72Nxw/fhw//OEPMX/+fGNW8Lx587B8+XKsX78e7e3taG9vx/r161FXV8cZwGQ748serlYKAbAelnIb44GyRU4lgC+++CJCoRCWLl2K0tJSY9m/fz8AwOFw4B//+AdWrlyJyspKrFu3DpWVlWhra4PX6zWO86tf/QqrVq3CD37wAzzwwAMoKCjAX/7yFzgcDqPN3r17MX/+fAQCAQQCAdxzzz149dVXMz5mIqskH6U0/u5GarnEZLZt24ZQKGQsPT095gbJ+rhEbZ6+JOrhxj0wWVJr6BL7SVwbq8vTNOM4IqJvA4B4XF8fj+v1fFocakTfpg260DpyJ85eKgRicSjROBCLA1qiHtDoW2L/xLGRfE5o4gHUY/WBKfV/KbWHool5W+q4TeeD9VfZIG3xYAXjOtVjqXWENzZyUU5NArnW0/Tz8/Px17/+9ZrHycvLw/PPP4/nn39+yjaFhYXYs2fP1+4jUa6oqKhASUkJWlpa8B//8R8AgEgkgsOHD+N//ud/ptzP4/HA4/FkqptEGcF4oGyTUwkgEd1cw8PD+Oyzz4z3Z86cQWdnJwoLCzF37lw0NjZi586duOOOO3DHHXdg586dKCgowNq1ay3sNVF6MB4olzABJKIpffjhh/jOd75jvE8+ymjdunX44x//iJ/97GcYHR3Fxo0bMTAwgEWLFuHtt982lVQQ5QrGA+USJoBENKWlS5detbRCURTs2LEDO3bsyFyniCzCeKBcklOTQIgoe43NjxDTq+nrlIkSkvLQZdMEEWMySFyfmJGcoKEJoOkTNSQW0ydixGJALAbE41CvRAEA7n4H3u67GwPnZ0KJxqBeiUGJ6m0QjRn76xM/NP1zNM2YfGKaaIKxiSqmSR9aykSRxDjGj/taNc1EmfJ2391Wd4HSgAkgERERkc0wASQiIiKyGSaARERERDbDBJCIpofkA481GXuQs2imr5P1dUadnWhjtXbJBzHH43q9XTyeePhzXN8mmlHDh1hMXxeLQaJRSFyDOhIBABQEFZz6VxnyzrmAcATKaBgIR/TjRKPG/tASNYTJWsOUz9fHkVKnqE1S3zfFGCecBxn3YGiiTJCx+tpT/yqzuDOUDkwAiYiIiGyGCSARERGRzTABJCIiIrIZJoBERERENsO/BEJE04I+iUM1Jj0YkzoAKGrKpI/kq2gQTR17b0wGSXkgs2iAouoPbk5O2gAgiQc7668aEI1CGR4FAHjPxRD3uJF/QSBXrkCJxyHhCCQWgxIO6/tHoolJJdGUCSliTD4x9wMT+mi8SmKiiL7SNMbk+Me2E2WYCAAN3i631T2hNOAdQCIiIiKbYQJIREREZDNMAImIiIhshjWARDQ9iAZATTy0Wa89MjbFodfYpT70ObXNJNuMhz8ropfXiQaJ6A97lkhU3x6JAqqi1wOGBgEAt3x6Ce7QLXAORyCXRyAul17rF41BwxV9/yvhiTV6yc9IPGz6K/fRdA5S1ok24TwQZZwI/B9dsboXlAa8A0hERERkM0wAiYiIiGyGCSARERGRzbAG0CLJPwofQxTgY74miCEKYOw8Ue4yYkGi0INBSdS+KeMbmreltkldJ5peNJh4VqDxc65oUESDSAyKKBCJQhEAogAShyoRxCUKiYcRizmBeASKRBLHjEKTGBSRxH7J/TXzNZr8POPz45P38WuMMYaY6TxRbpuO/zfEY1cgErW6GwD4f8PNxATQIv39/QCAIzhgcU+mt6GhIfh8Pqu7QWk0Fgtv6SuS39cn+/4+fptMsk0b95raJpZ4jY57BYDBxOvHiSUptf49PG6/GCYa34/J+vh1xpjAWLCHafl/w3v/a3UPJmA83DgmgBYpLCwEAHR3d/MiTjE4OIg5c+agu7sbiqKgrKzM6i5RmjEWppaMh66uLsaCTWRbPCSv0Z6eHsycOTPtnyciGBoaYjzcBEwALaKq+q+lfD5fRoIm2/C82Adj4dpuv/124zxRbsvWeJg5c2bG+psNiXE24HcUIiIiIpthAkhERERkM0wALeLxePDUU0/B4/FY3ZVphefFfvhvPjWeG/vJtn/zbOsvjVGEc6mJiIiIbIV3AImIiIhshgkgERERkc0wASQiIiKyGSaARERERDbDBNACu3fvRkVFBfLy8lBdXY333nvP6i6l1a5du3DffffB6/XC7/dj1apVOH36tKnNE088AUVRTMvixYtNbcLhMDZv3oyioiLMmDEDjz76KM6dO5fJoVAa2CkeGAt0NVbEQiavyYGBAdTX18Pn88Hn86G+vh6XLl1K9xBpCkwAM2z//v1obGzE9u3bcfz4cTz00EOora1Fd3e31V1Lm8OHD2PTpk1ob29HS0sLYrEYAoEALl++bGq3fPly9Pb2GsuBA+a/hdnY2Ig33ngD+/btw5EjRzA8PIy6ujrE4/FMDoduIrvFA2OBpmJVLGTymly7di06OzvR3NyM5uZmdHZ2or6+Pq3jo6sQyqj7779fNmzYYFp31113ydatWy3qUeb19fUJADl8+LCxbt26dbJy5cop97l06ZK4XC7Zt2+fse7zzz8XVVWlubk5nd2lNLJ7PDAWKGm6xEK6rsmuri4BIO3t7UabtrY2ASAff/zxzR8IXRPvAGZQJBJBR0cHAoGAaX0gEMDRo0ct6lXmhUIhAGN/9DyptbUVfr8flZWVWL9+Pfr6+oxtHR0diEajpnNXVlaGqqoqW527XMJ4YCyQbjrFQrquyba2Nvh8PixatMhos3jxYvh8Pl63FmECmEEXL15EPB5HcXGxaX1xcTGCwaBFvcosEUFTUxMefPBBVFVVGetra2uxd+9eHDx4EM888wyOHTuGZcuWIRwOAwCCwSDcbjdmzZplOp6dzl2usXs8MBYoabrEQjqvyWAwCL/fP+Ez/X4/r1uLOK3ugB0pimJ6LyIT1uWqhoYGnDhxAkeOHDGtX7NmjfF1VVUVFi5ciPLycrz11ltYvXr1lMez07nLVXaNB8YCjWd1LKT7mpxsLLxurcM7gBlUVFQEh8Mx4aedvr6+CT/55aLNmzfjzTffxKFDhzB79uyrti0tLUV5eTk+/fRTAEBJSQkikQgGBgZM7exy7nKRneOBsUCppkMspPuaLCkpwfnz5ycc68KFC7xuLcIEMIPcbjeqq6vR0tJiWt/S0oIlS5ZY1Kv0ExE0NDTg9ddfx8GDB1FRUXHNffr7+9HT04PS0lIAQHV1NVwul+nc9fb24uTJkzl97nKZHeOBsUCTsTIWMnVN1tTUIBQK4YMPPjDavP/++wiFQrxurWLJ1BMb27dvn7hcLnnppZekq6tLGhsbZcaMGXL27Fmru5Y2P/7xj8Xn80lra6v09vYay8jIiIiIDA0NyZYtW+To0aNy5swZOXTokNTU1Mjtt98ug4ODxnE2bNggs2fPlnfeeUc++ugjWbZsmdx7770Si8WsGhrdILvFA2OBpmJVLGTymly+fLncc8890tbWJm1tbTJ//nypq6tL6/hoakwALfDrX/9aysvLxe12y4IFC0zT7XMRgEmXl19+WURERkZGJBAIyDe+8Q1xuVwyd+5cWbdunXR3d5uOMzo6Kg0NDVJYWCj5+flSV1c3oQ1lHzvFA2OBrsaKWMjkNdnf3y+PP/64eL1e8Xq98vjjj8vAwEDax0iTU0REMn3XkYiIiIiswxpAIiIiIpthAkhERERkM0wAiYiIiGyGCSARERGRzTABJCIiIrIZJoBERERENsMEkIiIiMhmmAASERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGSaARERERDbDBJCIiIjIZpgAEhEREdkME0AiIiIim/n/AXi4qKym9bQDAAAAAElFTkSuQmCC",
      "text/html": [
       "\n",
       "            <div style=\"display: inline-block;\">\n",
       "                <div class=\"jupyter-widgets widget-label\" style=\"text-align: center;\">\n",
       "                    Figure\n",
       "                </div>\n",
       "                <img src='data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAABYXElEQVR4nO3dbWwU57k//u/MPtkmywbHXT8EcH2qOCEyyTmYBEweSlG7wacmEF4UidQi+qXoUDCSa6oKxIvQqAd0KiWNlJSoD2kaBVJ4k/SkCseNUzAJxU6Ig0WpCUkq+NskXgyOWdvY7NNc/xezO97xAyTA7nh3vh9ptN6Ze2bve5jLXB5f91gREQERERER2YZqdQeIiIiIKLOYABIRERHZDBNAIiIiIpthAkhERERkM0wAiYiIiGyGCSARERGRzTABJCIiIrIZJoBERERENsMEkIiIiMhmmAASERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGSaAN2D37t2oqKhAXl4eqqur8d5771ndJSIiIqJrYgJ4nfbv34/GxkZs374dx48fx0MPPYTa2lp0d3db3TUiIiKiq2ICeJ2effZZPPnkk/jRj36EefPm4bnnnsOcOXPw4osvWt01oozj3XCiMYwHygZOqzuQjSKRCDo6OrB161bT+kAggKNHj066TzgcRjgcNt5rmoYvv/wSt912GxRFSWt/s5GIYGhoCGVlZVBV/pwynSXvhu/evRsPPPAAfvOb36C2thZdXV2YO3fuNffXNA1ffPEFvF4vY2ESjIXswnhIL8bDTST0tX3++ecCQP7+97+b1v/3f/+3VFZWTrrPU089JQC4fM2lp6cnE/+kdAPuv/9+2bBhg2ndXXfdJVu3bv1K+/f09Fh+nWXDwljIDowHxkO24B3AGzD+pzMRmfIntm3btqGpqcl4HwqFMHfuXPx/H30TM2+Z+qeYxyrn35zOZpkYojiCA/B6vVZ3ha7ieu6Gj5f8N2YsTI6xkD0YD+nHeLh5mABeh6KiIjgcDgSDQdP6vr4+FBcXT7qPx+OBx+OZsH7mLSpmeqcOcqfiurHOZivRX/grkOnt4sWLiMfjE6774uLiCfGRNL4cYmhoCABjYUqMhazBeMgAxsNNw1+gXwe3243q6mq0tLSY1re0tGDJkiUW9YrIOl/nbviuXbvg8/mMZc6cOZnoIlHGMB4oGzABvE5NTU34/e9/jz/84Q84deoUfvKTn6C7uxsbNmywumtEGXM9d8O3bduGUChkLD09PZnoKlHaMR4om/BXwNdpzZo16O/vx9NPP43e3l5UVVXhwIEDKC8vt7prRBmTejf8scceM9a3tLRg5cqVk+4zVTkEUbZjPFA2YQJ4AzZu3IiNGzda3Q0iSzU1NaG+vh4LFy5ETU0Nfvvb3/JuONkW44GyBRNAIrohvBtONIbxQNmCCSAR3TDeDScaw3igbMAE0GKPVc6373R+ohSMBaIxjAdKN84CJiIiIrIZJoBERERENsMEkIiIiMhmmAASERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGSaARERERDbDBJCIiIjIZpgAEhEREdkME0AiIiIim2ECSERERGQzTACJiIiIbIYJIBEREZHNMAEkIiIishkmgEREREQ2wwSQiIiIyGaYABIRERHZDBNAIiIiIpthAkhERERkM0wAiYiIiGyGCSARERGRzTABJCIiIrIZJoBERERENsMEkIiIiMhmmAASERER2QwTQCIiIiKbYQJIREREZDM5lQDu2rUL9913H7xeL/x+P1atWoXTp0+b2jzxxBNQFMW0LF682NQmHA5j8+bNKCoqwowZM/Doo4/i3LlzpjYDAwOor6+Hz+eDz+dDfX09Ll26lO4hEhEREd2wnEoADx8+jE2bNqG9vR0tLS2IxWIIBAK4fPmyqd3y5cvR29trLAcOHDBtb2xsxBtvvIF9+/bhyJEjGB4eRl1dHeLxuNFm7dq16OzsRHNzM5qbm9HZ2Yn6+vqMjJOIiIjoRjit7sDN1NzcbHr/8ssvw+/3o6OjAw8//LCx3uPxoKSkZNJjhEIhvPTSS3j11Vfx3e9+FwCwZ88ezJkzB++88w4eeeQRnDp1Cs3NzWhvb8eiRYsAAL/73e9QU1OD06dP484770zTCImIiIhuXE7dARwvFAoBAAoLC03rW1tb4ff7UVlZifXr16Ovr8/Y1tHRgWg0ikAgYKwrKytDVVUVjh49CgBoa2uDz+czkj8AWLx4MXw+n9FmvHA4jMHBQdNCNN29++67WLFiBcrKyqAoCv785z+btosIduzYgbKyMuTn52Pp0qX45z//aU1nidKM8UC5JGcTQBFBU1MTHnzwQVRVVRnra2trsXfvXhw8eBDPPPMMjh07hmXLliEcDgMAgsEg3G43Zs2aZTpecXExgsGg0cbv90/4TL/fb7QZb9euXUa9oM/nw5w5c27WUInS5vLly7j33nvxwgsvTLr9l7/8JZ599lm88MILOHbsGEpKSvC9730PQ0NDGe4pUfoxHiiX5NSvgFM1NDTgxIkTOHLkiGn9mjVrjK+rqqqwcOFClJeX46233sLq1aunPJ6IQFEU433q11O1SbVt2zY0NTUZ7wcHB5kE0rRXW1uL2traSbeJCJ577jls377diJ1XXnkFxcXFeO211/Bf//VfmewqUdoxHiiX5OQdwM2bN+PNN9/EoUOHMHv27Ku2LS0tRXl5OT799FMAQElJCSKRCAYGBkzt+vr6UFxcbLQ5f/78hGNduHDBaDOex+PBzJkzTQtRNjtz5gyCwaCpXMLj8eDb3/72lKUQRLmK8UDZJqcSQBFBQ0MDXn/9dRw8eBAVFRXX3Ke/vx89PT0oLS0FAFRXV8PlcqGlpcVo09vbi5MnT2LJkiUAgJqaGoRCIXzwwQdGm/fffx+hUMhoQ5TrkuUO43/oSS2XmAzrYSkXMR4o2+RUArhp0ybs2bMHr732GrxeL4LBIILBIEZHRwEAw8PD+OlPf4q2tjacPXsWra2tWLFiBYqKivDYY48BAHw+H5588kls2bIFf/vb33D8+HH88Ic/xPz5841ZwfPmzcPy5cuxfv16tLe3o729HevXr0ddXR1nAJPtjC97uFopBMB6WMptjAfKFjmVAL744osIhUJYunQpSktLjWX//v0AAIfDgX/84x9YuXIlKisrsW7dOlRWVqKtrQ1er9c4zq9+9SusWrUKP/jBD/DAAw+goKAAf/nLX+BwOIw2e/fuxfz58xEIBBAIBHDPPffg1VdfzfiYiaySfJTS+LsbqeUSk9m2bRtCoZCx9PT0pLWfRJnAeKBsk1OTQETkqtvz8/Px17/+9ZrHycvLw/PPP4/nn39+yjaFhYXYs2fP1+4jUa6oqKhASUkJWlpa8B//8R8AgEgkgsOHD+N//ud/ptzP4/HA4/FkqptEGcF4oGyTUwkgEd1cw8PD+Oyzz4z3Z86cQWdnJwoLCzF37lw0NjZi586duOOOO3DHHXdg586dKCgowNq1ay3sNVF6MB4olzABJKIpffjhh/jOd75jvE8+ymjdunX44x//iJ/97GcYHR3Fxo0bMTAwgEWLFuHtt982lVQQ5QrGA+USRa71e1NKi8HBQfh8PizFSjgVl9XdmXZiEkUr/hehUIiPzMlxjIWrYyzYC+Ph6hgPN09OTQIhIiIiomtjAkhERERkM0wAiYiIiGyGCSARERGRzTABJCIiIrIZJoBERERENsMEkIiIiMhmmAASERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGSaARERERDbDBJCIiIjIZpgAEhEREdkME0AiIiIim2ECSERERGQzTqs7YFciAgCIIQqIxZ2ZhmKIAhg7T5S7GAtXx1iwF8bD1TEebh4mgBbp7+8HABzBAYt7Mr0NDQ3B5/NZ3Q1KI8bCV8NYsAfGw1fDeLhxTAAtUlhYCADo7u7mRZxicHAQc+bMQXd3NxRFQVlZmdVdojRjLEwtGQ9dXV2MBZvItnhIXqM9PT2YOXNm2j9PRDA0NMR4uAmYAFpEVfXyS5/Pl5GgyTY8L/bBWLi222+/3ThPlNuyNR5mzpyZsf5mQ2KcDfgdhYiIiMhmmAASERER2QwTQIt4PB489dRT8Hg8VndlWuF5sR/+m0+N58Z+su3fPNv6S2MU4VxqIiIiIlvhHUAiIiIim2ECSERERGQzTACJiIiIbIYJoAV2796NiooK5OXlobq6Gu+9957VXUqrXbt24b777oPX64Xf78eqVatw+vRpU5snnngCiqKYlsWLF5vahMNhbN68GUVFRZgxYwYeffRRnDt3LpNDoTSwUzwwFuhqrIiFTF6TAwMDqK+vh8/ng8/nQ319PS5dupTuIdIUmABm2P79+9HY2Ijt27fj+PHjeOihh1BbW4vu7m6ru5Y2hw8fxqZNm9De3o6WlhbEYjEEAgFcvnzZ1G758uXo7e01lgMHzH8KqbGxEW+88Qb27duHI0eOYHh4GHV1dYjH45kcDt1EdosHxgJNxapYyOQ1uXbtWnR2dqK5uRnNzc3o7OxEfX19WsdHVyGUUffff79s2LDBtO6uu+6SrVu3WtSjzOvr6xMAcvjwYWPdunXrZOXKlVPuc+nSJXG5XLJv3z5j3eeffy6qqkpzc3M6u0tpZPd4YCxQ0nSJhXRdk11dXQJA2tvbjTZtbW0CQD7++OObPxC6Jt4BzKBIJIKOjg4EAgHT+kAggKNHj1rUq8wLhUIAxv7mZVJrayv8fj8qKyuxfv169PX1Gds6OjoQjUZN566srAxVVVW2One5hPHAWCDddIqFdF2TbW1t8Pl8WLRokdFm8eLF8Pl8vG4twgQwgy5evIh4PI7i4mLT+uLiYgSDQYt6lVkigqamJjz44IOoqqoy1tfW1mLv3r04ePAgnnnmGRw7dgzLli1DOBwGAASDQbjdbsyaNct0PDudu1xj93hgLFDSdImFdF6TwWAQfr9/wmf6/X5etxZxWt0BO1IUxfReRCasy1UNDQ04ceIEjhw5Ylq/Zs0a4+uqqiosXLgQ5eXleOutt7B69eopj2enc5er7BoPjAUaz+pYSPc1OdlYeN1ah3cAM6ioqAgOh2PCTzt9fX0TfvLLRZs3b8abb76JQ4cOYfbs2VdtW1paivLycnz66acAgJKSEkQiEQwMDJja2eXc5SI7xwNjgVJNh1hI9zVZUlKC8+fPTzjWhQsXeN1ahAlgBrndblRXV6OlpcW0vqWlBUuWLLGoV+knImhoaMDrr7+OgwcPoqKi4pr79Pf3o6enB6WlpQCA6upquFwu07nr7e3FyZMnc/rc5TI7xgNjgSZjZSxk6pqsqalBKBTCBx98YLR5//33EQqFeN1axZKpJza2b98+cblc8tJLL0lXV5c0NjbKjBkz5OzZs1Z3LW1+/OMfi8/nk9bWVunt7TWWkZEREREZGhqSLVu2yNGjR+XMmTNy6NAhqampkdtvv10GBweN42zYsEFmz54t77zzjnz00UeybNkyuffeeyUWi1k1NLpBdosHxgJNxapYyOQ1uXz5crnnnnukra1N2traZP78+VJXV5fW8dHUmABa4Ne//rWUl5eL2+2WBQsWmKbb5yIAky4vv/yyiIiMjIxIIBCQb3zjG+JyuWTu3Lmybt066e7uNh1ndHRUGhoapLCwUPLz86Wurm5CG8o+dooHxgJdjRWxkMlrsr+/Xx5//HHxer3i9Xrl8ccfl4GBgbSPkSaniIhk+q4jEREREVmHNYBERERENsMEkIiIiMhmmAASERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGSaARERERDbDBJCIiIjIZpgAEhEREdkME0AiIiIim2ECSERERGQzTACJiIiIbIYJIBEREZHNMAEkIiIishkmgEREREQ2wwSQiIiIyGaYAN6A3bt3o6KiAnl5eaiursZ7771ndZeIiIiIrokJ4HXav38/GhsbsX37dhw/fhwPPfQQamtr0d3dbXXXiIiIiK6KCeB1evbZZ/Hkk0/iRz/6EebNm4fnnnsOc+bMwYsvvmh114gyjnfDicYwHigbMAG8DpFIBB0dHQgEAqb1gUAAR48etahXRNbg3XCiMYwHyhaKiIjVncg2X3zxBW6//Xb8/e9/x5IlS4z1O3fuxCuvvILTp09P2CccDiMcDhvvNU3Dl19+idtuuw2KomSk39lERDA0NISysjKoKn9Omc4WLVqEBQsWmO5+z5s3D6tWrcKuXbuuub+mafjiiy/g9XoZC5NgLGQXxkN6MR5uHqfVHchm44NTRKYM2F27duHnP/95JrqVU3p6ejB79myru0FTSN4N37p1q2n917kb/sUXX2DOnDnp6F5OYSxMf4yHzGE83DgmgNehqKgIDocDwWDQtL6vrw/FxcWT7rNt2zY0NTUZ70OhEObOnYsH8Z9wwpXW/majGKI4ggPwer1Wd4Wu4uLFi4jH4xOu++Li4gnxkTT+bnjylxCMhckxFrIH4yH9GA83DxPA6+B2u1FdXY2WlhY89thjxvqWlhasXLly0n08Hg88Hs+E9U644FQY5BMkChP4K5DscDPuhjMWpsBYyDqMhzRiPNw0/AX6dWpqasLvf/97/OEPf8CpU6fwk5/8BN3d3diwYYPVXSPKmOu9Gx4KhYylp6cnE10lSjvGA2UT3gG8TmvWrEF/fz+efvpp9Pb2oqqqCgcOHEB5ebnVXSPKmJt5NxzJn+gVBRAZe5+UXDf+dfw21QFocf1VNEBJ/JwrGhSHAxKLQXG5IbEoFKcLUBUgHodaUID44CAcd1ci4r8FzuEIlJOfAS4XEI1Ci0ShuJyQcHhsf4cDoqXMo0t+XvLztfjkffy6Y+RUvaxwU+NhGokvXQBH60dWd4NuMiaAN2Djxo3YuHGj1d0gslRTUxPq6+uxcOFC1NTU4Le//S3vhpNtMR4oWzABJKIbwrvhRGMYD5QtmAAS0Q3j3XCiMYwHygacBEJE00OiVk9xOPS6N0U1luQ6xeGYos3EbYrDMba4nPpx3G69jdulr3e7oLjdUJxOKL6ZAIDhO25F34I8DMzzQplRAHVGAZT8fCguJ9T8PH3/PI++v9M58TNczq/Vx9RxprZJvjdqGImsoCjoW5BndS8oDfidhYiIiMhmmAASERER2QwTQCIiIiKbYQJIREREZDNMAIloWjBNfki8V1QFipoyESL1VVHN240JFsmJI4l1qqJP1ki8AtAnaiQmbigeN+ByQW7JBwAMzXZi6O4IBv9NgZKXBxTk65M+nE4g8cBexe3SP9PlguJQE0uivxP64Rib7JEcY8pYk2OYbPz6W/7JK7JIYhLS0N0Rq3tCacAEkIiIiMhmmAASERER2QwTQCIiIiKbYQJIRNNDstYtWdOnpjxEOfG14kjWx6XUAqY+XBnQa/tURX9NeSC0XrOX+ONHiQc4w+mEkqjj0wr0h0SPlAjmfesLXJkdBTxuSL4H8LgTD3t2GfsnawvhcIwtic/Xx6HX9yG5AFCUlNcpxjjhPPBB0GSF5LWqKpj3rS8s7gylA7+zEBEREdkME0AiIiIim2ECSERERGQzTACJiIiIbIYJIBFNC4mac/NECZjXQVWNVyV1coWqjB1ATXmQssOhr088EBpqYsKI0wkoiUkciYkcWp4+wSNyWxwBfxdmFQ9CXE5oeU5I4sHRyUkkitOp9ynxAGioqvkB1MDY5BRF0RdVMfXfGFdiHOPHnTp+IisF/F1Wd4HSgAkgERERkc0wASQiIiKyGSaARERERDbDBJCIppfUOjnjQcqquc4vIVlfl7pfshYvuV/yOPpDoVPq89TEQ6QdKqA6oLn1berMKJYWnMY3b/0ScDogLgfg1GsIxZHaN9WoLVSSdYaA8eDmsYc9p3ybHSt0NNcEjh+36XywFpAyzLhO9VhaWnDa2v5QWjABJCIiIrIZJoBERERENsMEkIiIiMhmmAASERER2QwTQBrDB8+SlVImPyjXmviQeq2OnzQxoalimnwx4TUxSURz6Os8eVFUuhSU5A9BHImJH6o61hbQJ44kJ6AkJ5wAY5NWrjI20/vUSSzJcadMarnmeSDKgEoXr8NcxASQzJgEEhER5TwmgKRL3skgIiKinMf/8Umn6H/TVEk+y4yIiIhyltPqDtidWpAPjMagzpgBCYehuN2QeByA/sfsJRKB4vFAu3x5rI3TCagqJBqDmp+H+OAgHLf6oF0ehXrLDEgsptcVuV2Q4ctQZ92KWG8QzjmzIaFBKL6ZgMsJuTyCK/fOhevtD9H7eiVO3P8nbOldgK7/Nw8xrwfOS1egBC8Ahbci/sm/4PzmXGgX+qF6b4GEI4BogNMJbWgY6i0zEO//0uiH4tIvLYlEoOTnQxsagur1QkZHoXg8QGKMUFXzGAsK9PeqCoSt+lchS11PGYKimmrpkusAbcrPMNUGAkbtntMZR4HqRr4aARQFogCiKLhmr8Z//lepBTR2ZekFTV8FqtvqLlAaMAG0mDYyClVxQbt8GQAgsZixTZKviXWTtYlHI/rrpZD+OhCZ+Bm9QQBArOecvmJw0NjmersPAFC66hQewb9D/w/zFEz/dV7s1/c/223qR6p4OGzqh0TH+iFDQ/p+idfU/httkmMcGdHfS3RCGyIiIro5+CtgIiIiIpthAkhERERkMzmVAO7atQv33XcfvF4v/H4/Vq1ahdOnzX/E+oknnjD+gHxyWbx4salNOBzG5s2bUVRUhBkzZuDRRx/FuXPnTG0GBgZQX18Pn88Hn8+H+vp6XLp0Kd1DJCIiIrphOZUAHj58GJs2bUJ7eztaWloQi8UQCARweVzN2vLly9Hb22ssBw4cMG1vbGzEG2+8gX379uHIkSMYHh5GXV0d4smJCwDWrl2Lzs5ONDc3o7m5GZ2dnaivr8/IOIlymsi120zYR4OM30+mmACS+AwRMX+Wpn8dizkwokUwqrkBESgCKF+lT+PbaJPso03epwl9J5pGRrSJteWU/XJqEkhzc7Pp/csvvwy/34+Ojg48/PDDxnqPx4OSkpJJjxEKhfDSSy/h1VdfxXe/+10AwJ49ezBnzhy88847eOSRR3Dq1Ck0Nzejvb0dixYtAgD87ne/Q01NDU6fPo0777wzTSMkIiIiunE5dQdwvFBIn5FaWFhoWt/a2gq/34/KykqsX78efX19xraOjg5Eo1EEAgFjXVlZGaqqqnD06FEAQFtbG3w+n5H8AcDixYvh8/mMNkRERETTVc4mgCKCpqYmPPjgg6iqqjLW19bWYu/evTh48CCeeeYZHDt2DMuWLUM48RiTYDAIt9uNWbNmmY5XXFyMYDBotPH7/RM+0+/3G23GC4fDGBwcNC1E0927776LFStWoKysDIqi4M9//rNpu4hgx44dKCsrQ35+PpYuXYp//vOf1nSWKM0YD5RLcjYBbGhowIkTJ/CnP/3JtH7NmjX4/ve/j6qqKqxYsQL/93//h08++QRvvfXWVY8nIuY/2j7Jg1vHt0m1a9cuY8KIz+fDnDlzrmNURJl1+fJl3HvvvXjhhRcm3f7LX/4Szz77LF544QUcO3YMJSUl+N73voehxDMfiXIJ44FySU4mgJs3b8abb76JQ4cOYfbs2VdtW1paivLycnz66acAgJKSEkQiEQwMDJja9fX1obi42Ghz/vz5Cce6cOGC0Wa8bdu2IRQKGUtPT8/1DI0oo2pra/GLX/wCq1evnrBNRPDcc89h+/btWL16NaqqqvDKK69gZGQEr7322tf/sJQJEjLZBArzh0+63+RNUyZ7TPaqaYAmUOP6uvAVFz6JCoKjXihxgRLX9DbJtgAQ11ImkmhjEz40+WqTP5Lvk8dIHXfyvaZd+zxQRmU0HqaRT6K8DnNRTiWAIoKGhga8/vrrOHjwICoqKq65T39/P3p6elBaWgoAqK6uhsvlQktLi9Gmt7cXJ0+exJIlSwAANTU1CIVC+OCDD4w277//PkKhkNFmPI/Hg5kzZ5oWomx25swZBINBU72sx+PBt7/97avWwrIcgnIR44GyTU4lgJs2bcKePXvw2muvwev1IhgMIhgMYnR0FAAwPDyMn/70p2hra8PZs2fR2tqKFStWoKioCI899hgAwOfz4cknn8SWLVvwt7/9DcePH8cPf/hDzJ8/35gVPG/ePCxfvhzr169He3s72tvbsX79etTV1XEGMNlGst51/F3v1HrZybAcgnIR44GyTU4lgC+++CJCoRCWLl2K0tJSY9m/fz8AwOFw4B//+AdWrlyJyspKrFu3DpWVlWhra4PX6zWO86tf/QqrVq3CD37wAzzwwAMoKCjAX/7yFzgcDqPN3r17MX/+fAQCAQQCAdxzzz149dVXMz5mIquNr3u9Wi0swHIIym2MB8oWOfUcwGs9TDU/Px9//etfr3mcvLw8PP/883j++eenbFNYWIg9e/Z87T4S5YrkszSDwaBRQgGY62Un4/F44PF4pj5wsj5OS6mtU/R6uPEPZJZJagElro3V5RnH0mvtJJ6s4Yvr2+JxQFUBLQ41om/TBl1oHbkTZy8VoiQWhxKNA7E4oCXqAY2+JfZPHFtJPig+8QDqsfrAlPq/lNpD0QRK6rbUcZvOB+uvskHa4sEKIoCi6NevKGgd4W+2clFO3QEkosypqKhASUmJqV42Eong8OHDU9bCEuUqxgNlm5y6A0hEN9fw8DA+++wz4/2ZM2fQ2dmJwsJCzJ07F42Njdi5cyfuuOMO3HHHHdi5cycKCgqwdu1aC3tNlB6MB8olTACJaEoffvghvvOd7xjvm5qaAADr1q3DH//4R/zsZz/D6OgoNm7ciIGBASxatAhvv/22qaaWKFcwHiiXMAEkoiktXbr0qrW1iqJgx44d2LFjR+Y6RWQRxgPlEtYAEtG0MDY/Qkyvpq9TJkpIykOXxz9AGQAkHtcnaYgYkzWg6RM1JBbTJ2LEYkAsBsTjUK9EAQDufgfe7rsbA+dnQonGoF6JQYnqbRCNGfuLCBDX9M/RNGPyiWmiSaLvkpj0MX6ih4gYD48eP+5rTWojypS3++62uguUBkwAiYiIiGyGCSARERGRzTABJCIiIrIZJoBEND0kH3isydiDnEUzfZ2srzPq7EQbq7VLPog5Htfr7eLxxMOf4/o20YwaPsRi+rpYDBKNQuIa1JEIAKAgqODUv8qQd84FhCNQRsNAOKIfJxo19oeWqCFM1hqmfL4+jpQ6RW2S+r4pxjjhPMi4B0MTZYKM1dee+leZxZ2hdGACSERERGQzTACJiIiIbIYJIBEREZHNMAEkIiIishn+JRAimhb0SRyqMenBmNQBQFFTJn0kX0WDaOrYe2MySMoDmUUDFFV/cHNy0gYASTzYWX/VgGgUyvAoAMB7Loa4x438CwK5cgVKPA4JRyCxGJRwWN8/Ek1MKommTEgRY/KJuR+Y0EfjVRITRfSVpjEmxz+2nSjDRABo8Ha5re4JpQHvABIRERHZDBNAIiIiIpthAkhERERkM6wBJKLpQTQAauKhzXrtkbEpDr3GLvWhz6ltJtlmPPxZEb28TjRIRH/Ys0Si+vZIFFAVvR4wNAgAuOXTS3CHboFzOAK5PAJxufRav2gMGq7o+18JT6zRS35G4mHTX7mPpnOQsk60CeeBKONE4P/oitW9oDTgHUAiIiIim2ECSERERGQzTACJiIiIbIY1gBZJ/lH4GKIAH/M1QQxRAGPniXKXEQsShR4MSqL2TRnf0LwttU3qOtH0osHEswKNn3NFgyIaRGJQRIFIFIoAEAWQOFSJIC5RSDyMWMwJxCNQJJI4ZhSaxKCIJPZL7q+Zr9Hk5xmfH5+8j19jjDHETOeJctt0/L8hHrsCkajV3QDA/xtuJiaAFunv7wcAHMEBi3syvQ0NDcHn81ndDUqjsVh4S1+R/L4+2ff38dtkkm3auNfUNrHEa3TcKwAMJl4/TixJqfXv4XH7xTDR+H5M1sevM8YExoI9TMv/G977X6t7MAHj4cYxAbRIYWEhAKC7u5sXcYrBwUHMmTMH3d3dUBQFZWVlVneJ0oyxMLVkPHR1dTEWbCLb4iF5jfb09GDmzJlp/zwRwdDQEOPhJmACaBFV1X8t5fP5MhI02YbnxT4YC9d2++23G+eJclu2xsPMmTMz1t9sSIyzAb+jEBEREdkME0AiIiIim2ECaBGPx4OnnnoKHo/H6q5MKzwv9sN/86nx3NhPtv2bZ1t/aYwinEtNREREZCu8A0hERERkM0wAiYiIiGyGCSARERGRzTABJCIiIrIZJoAW2L17NyoqKpCXl4fq6mq89957VncprXbt2oX77rsPXq8Xfr8fq1atwunTp01tnnjiCSiKYloWL15sahMOh7F582YUFRVhxowZePTRR3Hu3LlMDoXSwE7xwFigq7EiFjJ5TQ4MDKC+vh4+nw8+nw/19fW4dOlSuodIU2ACmGH79+9HY2Mjtm/fjuPHj+Ohhx5CbW0turu7re5a2hw+fBibNm1Ce3s7WlpaEIvFEAgEcPnyZVO75cuXo7e311gOHDD/LczGxka88cYb2LdvH44cOYLh4WHU1dUhHo9ncjh0E9ktHhgLNBWrYiGT1+TatWvR2dmJ5uZmNDc3o7OzE/X19WkdH12FUEbdf//9smHDBtO6u+66S7Zu3WpRjzKvr69PAMjhw4eNdevWrZOVK1dOuc+lS5fE5XLJvn37jHWff/65qKoqzc3N6ewupZHd44GxQEnTJRbSdU12dXUJAGlvbzfatLW1CQD5+OOPb/5A6Jp4BzCDIpEIOjo6EAgETOsDgQCOHj1qUa8yLxQKARj7o+dJra2t8Pv9qKysxPr169HX12ds6+joQDQaNZ27srIyVFVV2erc5RLGA2OBdNMpFtJ1Tba1tcHn82HRokVGm8WLF8Pn8/G6tQgTwAy6ePEi4vE4iouLTeuLi4sRDAYt6lVmiQiamprw4IMPoqqqylhfW1uLvXv34uDBg3jmmWdw7NgxLFu2DOFwGAAQDAbhdrsxa9Ys0/HsdO5yjd3jgbFASdMlFtJ5TQaDQfj9/gmf6ff7ed1axGl1B+xIURTTexGZsC5XNTQ04MSJEzhy5Ihp/Zo1a4yvq6qqsHDhQpSXl+Ott97C6tWrpzyenc5drrJrPDAWaDyrYyHd1+RkY+F1ax3eAcygoqIiOByOCT/t9PX1TfjJLxdt3rwZb775Jg4dOoTZs2dftW1paSnKy8vx6aefAgBKSkoQiUQwMDBgameXc5eL7BwPjAVKNR1iId3XZElJCc6fPz/hWBcuXOB1axEmgBnkdrtRXV2NlpYW0/qWlhYsWbLEol6ln4igoaEBr7/+Og4ePIiKiopr7tPf34+enh6UlpYCAKqrq+FyuUznrre3FydPnszpc5fL7BgPjAWajJWxkKlrsqamBqFQCB988IHR5v3330coFOJ1axVLpp7Y2L59+8TlcslLL70kXV1d0tjYKDNmzJCzZ89a3bW0+fGPfyw+n09aW1ult7fXWEZGRkREZGhoSLZs2SJHjx6VM2fOyKFDh6SmpkZuv/12GRwcNI6zYcMGmT17trzzzjvy0UcfybJly+Tee++VWCxm1dDoBtktHhgLNBWrYiGT1+Ty5cvlnnvukba2Nmlra5P58+dLXV1dWsdHU2MCaIFf//rXUl5eLm63WxYsWGCabp+LAEy6vPzyyyIiMjIyIoFAQL7xjW+Iy+WSuXPnyrp166S7u9t0nNHRUWloaJDCwkLJz8+Xurq6CW0o+9gpHhgLdDVWxEImr8n+/n55/PHHxev1itfrlccff1wGBgbSPkaanCIikum7jkRERERkHdYAEhEREdkME0AiIiIim2ECSERERGQzTACJiIiIbIYJIBEREZHNMAEkIiIishkmgEREREQ2wwSQiIiIyGaYABIRERHZDBNAIiIiIpthAkhERERkM0wAiYiIiGyGCSARERGRzTABJCIiIrIZJoBERERENsMEkIiIiMhmmAASERER2QwTwBuwe/duVFRUIC8vD9XV1Xjvvfes7hIRERHRNTEBvE779+9HY2Mjtm/fjuPHj+Ohhx5CbW0turu7re4aERER0VUpIiJWdyIbLVq0CAsWLMCLL75orJs3bx5WrVqFXbt2WdgzIiIioqvjHcDrEIlE0NHRgUAgYFofCARw9OhRi3pFZB2WQxCNYTxQNnBa3YFsdPHiRcTjcRQXF5vWFxcXIxgMTrpPOBxGOBw23muahi+//BK33XYbFEVJa3+zkYhgaGgIZWVlUFX+nDKdJcshdu/ejQceeAC/+c1vUFtbi66uLsydO/ea+2uahi+++AJer5exMAnGQnZhPKQX4+EmEvraPv/8cwEgR48eNa3/xS9+IXfeeeek+zz11FMCgMvXXHp6ejLxT0o34P7775cNGzaY1t11112ydevWr7R/T0+P5ddZNiyMhezAeGA8ZAveAbwORUVFcDgcE+729fX1TbgrmLRt2zY0NTUZ70OhEObOnYsH8Z9wwpXW/majGKI4ggPwer1Wd4WuIlkOsXXrVtP6q5VDjL8bLokyZMbC5BgL2YPxkH6Mh5uHCeB1cLvdqK6uRktLCx577DFjfUtLC1auXDnpPh6PBx6PZ8J6J1xwKgzyCfTvgfwVyDR3PeUQu3btws9//vMJ6xkLU2AsZA3GQwYwHm4a/gL9OjU1NeH3v/89/vCHP+DUqVP4yU9+gu7ubmzYsMHqrhFl3PhvxiIy5Tfobdu2IRQKGUtPT08mukiUMYwHyga8A3id1qxZg/7+fjz99NPo7e1FVVUVDhw4gPLycqu7RpQx11MOMdXdcCT/g1QUQGTsfVJy3fjX8dtUB6DF9VfRACXxc65oUBwOSCwGxeWGxKJQnC5AVYB4HGpBAeKDg3DcXYmI/xY4hyNQTn4GuFxANAotEoXickLC4bH9HQ6IJil9THxe8vO1+OR9/LpjFFAWuKnxMI3Ely6Ao/Ujq7tBNxnvAN6AjRs34uzZswiHw+jo6MDDDz9sdZeIMiq1HCJVS0sLlixZYlGviKzBeKBswjuARHRDmpqaUF9fj4ULF6Kmpga//e1vWQ5BtsV4oGzBBJCIbgjLIYjGMB4oWzABJKIbtnHjRmzcuPHGDqKogECvq4vHx2r3ACiqAonHx2r4JrTRzNsS9X4QRW+jKpBoDIrbrbdxu/SaQLcLUFUgGoXimwkMDmL4jltx6VtO5F/w4LYzBVBcLkg4AkUEan4e4uEwlDwPMKpBcTqBuKZ3IfkZLickHDf6MWUfk+NJrSFMaaPvpwFQWQOYZW5KPEwXioK+BXkobbW6I3SzsQaQiIiIyGaYABIRERHZDBNAIiIiIpthAkhERERkM0wAiWhaUByOxBeq8V5RFSiqMjaRIvVVUc3bE/srjsQ2R2KdqkBxOvW2Tn3em+JyAg4HFJcTiscNuFyQW/IBAEOznRi6O4LBf1Og5OUBBflQ8jz6vokH9ipul/6ZLhcUh5pYEv2d0A+Hvi6lj0gZa3IMk41ff8s/eUUWUfTYGro7YnVPKA2YABIRERHZDBNAIiIiIpthAkhERERkM0wAiWh6SNa6qWMPb07W+iW/VhzJ+riUWkCj5m6svk5RFf1VUcw1eK7Es++dTn2d0wklUcenFbgBACMlgnnf+gJXZkcBjxuS7wE8bv04Lpexf7K2EA7H2JL4fH0cen0fkgsARUl5nWKME86Dwm/TZIHktaoqmPetLyzuDKUDv7MQERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgEU0LiZpz80QJmNdBVY1XJXVyhaqMHUBNeZCyw6GvTzwQGmpiwojTCSiJSRyJiRxanj7BI3JbHAF/F2YVD0JcTmh5TkjiwdHJSSSK06n3KfEAaKiq+QHUwNjkFEXRF1Ux9d8YV2Ic48edOn4iKwX8XVZ3gdKACSARERGRzTABJCIiIrIZJoBERERENsMEkIiml9Q6OeNByqq5zi8hWV+Xul+yFi+5X/I4+kOhU+rz1MRDpB0qoDqgufVt6swolhacxjdv/RJwOiAuB+DUawjFkdo31agtVJJ1hoDx4Oaxhz2nfJsdK3Q01wSOH7fpfLAWkDLMuE71WFpacNra/lBaMAEkIiIishkmgEREREQ2wwSQiIiIyGaYABIRERHZDBNAGsMHz5KVUiY/KNea+JB6rY6fNDGhqWKafDHhNTFJRHPo6zx5UVS6FJTkD0EciYkfqjrWFtAnjiQnoCQnnABjk1auMjbT+9RJLMlxp0xqueZ5IMqAShevw1zEBJDMmAQSERHlPCaApEveySAiIqKcx//xSafof9NUST7LjIiIiHKW0+oO2J1akA+MxqDOmAEJh6G43ZB4HID+x+wlEoHi8UC7fHmsjdMJqCokGoOan4f44CAct/qgXR6FessMSCym1xW5XZDhy1Bn3YpYbxDOObMhoUEovpmAywm5PIIr986F6+0P0ft6JU7c/yds6V2Arv83DzGvB85LV6AELwCFtyL+yb/g/OZcaBf6oXpvgYQjgGiA0wltaBjqLTMQ7//S6Ifi0i8tiUSg5OdDGxqC6vVCRkeheDxAYoxQVfMYCwr096oKhK36VyFLXU8ZgqKaaumS6wBtys8w1QYCRu2e0xlHgepGvhoBFAWiAKIouGavxn/+V6kFNHZl6QVNXwWq2+ouUBowAbSYNjIKVXFBu3wZACCxmLFNkq+JdZO1iUcj+uulkP46EJn4Gb1BAECs55y+YnDQ2OZ6uw8AULrqFB7Bv0P/D/MUTP91XuzX9z/bbepHqng4bOqHRMf6IUND+n6J19T+G22SYxwZ0d9LdEIbIiIiujn4K2AiIiIim2ECSERERGQzOZUA7tq1C/fddx+8Xi/8fj9WrVqF06fNf8T6iSeeMP6AfHJZvHixqU04HMbmzZtRVFSEGTNm4NFHH8W5c+dMbQYGBlBfXw+fzwefz4f6+npcunQp3UMkIiIiumE5lQAePnwYmzZtQnt7O1paWhCLxRAIBHB5XM3a8uXL0dvbaywHDhwwbW9sbMQbb7yBffv24ciRIxgeHkZdXR3iyYkLANauXYvOzk40NzejubkZnZ2dqK+vz8g4iXKayLXbTNhHg4zfT6aYAJL4DBExf5amfx2LOTCiRTCquQERKAIoX6VP49tok+yjTd6nCX0nmkZGtIm15ZT9cmoSSHNzs+n9yy+/DL/fj46ODjz88MPGeo/Hg5KSkkmPEQqF8NJLL+HVV1/Fd7/7XQDAnj17MGfOHLzzzjt45JFHcOrUKTQ3N6O9vR2LFi0CAPzud79DTU0NTp8+jTvvvDNNIyQiIiK6cTl1B3C8UEifkVpYWGha39raCr/fj8rKSqxfvx59fX3Gto6ODkSjUQQCAWNdWVkZqqqqcPToUQBAW1sbfD6fkfwBwOLFi+Hz+Yw2RERERNNVziaAIoKmpiY8+OCDqKqqMtbX1tZi7969OHjwIJ555hkcO3YMy5YtQzjxGJNgMAi3241Zs2aZjldcXIxgMGi08fv9Ez7T7/cbbcYLh8MYHBw0LUTT3bvvvosVK1agrKwMiqLgz3/+s2m7iGDHjh0oKytDfn4+li5din/+85/WdJYozRgPlEtyNgFsaGjAiRMn8Kc//cm0fs2aNfj+97+PqqoqrFixAv/3f/+HTz75BG+99dZVjyci5j/aPsmDW8e3SbVr1y5jwojP58OcOXOuY1REmXX58mXce++9eOGFFybd/stf/hLPPvssXnjhBRw7dgwlJSX43ve+h6HEMx+JcgnjgXJJTiaAmzdvxptvvolDhw5h9uzZV21bWlqK8vJyfPrppwCAkpISRCIRDAwMmNr19fWhuLjYaHP+/PkJx7pw4YLRZrxt27YhFAoZS09Pz/UMjSijamtr8Ytf/AKrV6+esE1E8Nxzz2H79u1YvXo1qqqq8Morr2BkZASvvfba1/+wlAkSMtkECvOHT7rf5E1TJntM9qppgCZQ4/q68BUXPokKgqNeKHGBEtf0Nsm2ABDXUiaSaGMTPjT5apM/ku+Tx0gdd/K9pl37PFBGZTQeppFPorwOc1FOJYAigoaGBrz++us4ePAgKioqrrlPf38/enp6UFpaCgCorq6Gy+VCS0uL0aa3txcnT57EkiVLAAA1NTUIhUL44IMPjDbvv/8+QqGQ0WY8j8eDmTNnmhaibHbmzBkEg0FTvazH48G3v/3tq9bCshyCchHjgbJNTiWAmzZtwp49e/Daa6/B6/UiGAwiGAxidHQUADA8PIyf/vSnaGtrw9mzZ9Ha2ooVK1agqKgIjz32GADA5/PhySefxJYtW/C3v/0Nx48fxw9/+EPMnz/fmBU8b948LF++HOvXr0d7ezva29uxfv161NXVcQYw2Uay3nX8Xe/UetnJsByCchHjgbJNTiWAL774IkKhEJYuXYrS0lJj2b9/PwDA4XDgH//4B1auXInKykqsW7cOlZWVaGtrg9frNY7zq1/9CqtWrcIPfvADPPDAAygoKMBf/vIXOBwOo83evXsxf/58BAIBBAIB3HPPPXj11VczPmYiq42ve71aLSzAcgjKbYwHyhY59RzAaz1MNT8/H3/961+veZy8vDw8//zzeP7556dsU1hYiD179nztPhLliuSzNIPBoFFCAZjrZSfj8Xjg8XimPnCyPk5Lqa1T9Hq48Q9klklqASWujdXlGcfSa+0knqzhi+vb4nFAVQEtDjWib9MGXWgduRNnLxWiJBaHEo0DsTigJeoBjb4l9k8cW0k+KD7xAOqx+sCU+r+U2kPRBErqttRxm84H66+yQdriwQoigKLo168oaB3hb7ZyUU7dASSizKmoqEBJSYmpXjYSieDw4cNT1sIS5SrGA2WbnLoDSEQ31/DwMD777DPj/ZkzZ9DZ2YnCwkLMnTsXjY2N2LlzJ+644w7ccccd2LlzJwoKCrB27VoLe02UHowHyiVMAIloSh9++CG+853vGO+bmpoAAOvWrcMf//hH/OxnP8Po6Cg2btyIgYEBLFq0CG+//bapppYoVzAeKJcwASSiKS1duvSqtbWKomDHjh3YsWNH5jpFZBHGA+US1gAS0bQwNj9CTK+mr1MmSkjKQ5fHP0AZACQe1ydpiBiTNaDpEzUkFtMnYsRiQCwGxONQr0QBAO5+B97uuxsD52dCicagXolBieptEI0Z+4sIENf0z9E0Y/KJaaJJou+SmPQxfqKHiBgPjx4/7mtNaiPKlLf77ra6C5QGTACJiIiIbIYJIBEREZHNMAEkIiIishkmgEQ0PSQfeKzJ2IOcRTN9nayvM+rsRBurtUs+iDke1+vt4vHEw5/j+jbRjBo+xGL6ulgMEo1C4hrUkQgAoCCo4NS/ypB3zgWEI1BGw0A4oh8nGjX2h5aoIUzWGqZ8vj6OlDpFbZL6vinGOOE8yLgHQxNlgozV1576V5nFnaF0YAJIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGf4lECKaFvRJHKox6cGY1AFAUVMmfSRfRYNo6th7YzJIygOZRQMUVX9wc3LSBgBJPNhZf9WAaBTK8CgAwHsuhrjHjfwLArlyBUo8DglHILEYlHBY3z8STUwqiaZMSBFj8om5H5jQR+NVEhNF9JWmMSbHP7adKMNEAGjwdrmt7gmlAe8AEhEREdkME0AiIiIim2ECSERERGQzrAEkoulBNABq4qHNeu2RsSkOvcYu9aHPqW0m2WY8/FkRvbxONEhEf9izRKL69kgUUBW9HjA0CAC45dNLcIdugXM4Ark8AnG59Fq/aAwaruj7XwlPrNFLfkbiYdNfuY+mc5CyTrQJ54Eo40Tg/+iK1b2gNOAdQCIiIiKbYQJIREREZDNMAImIiIhshjWAFkn+UfgYogAf8zVBDFEAY+eJcpcRCxKFHgxKovZNGd/QvC21Teo60fSiwcSzAo2fc0WDIhpEYlBEgUgUigAQBZA4VIkgLlFIPIxYzAnEI1AkkjhmFJrEoIgk9kvur5mv0eTnGZ8fn7yPX2OMMcRM54ly23T8vyEeuwKRqNXdAMD/G24mJoAW6e/vBwAcwQGLezK9DQ0NwefzWd0NSqOxWHhLX5H8vj7Z9/fx22SSbdq419Q2scRrdNwrAAwmXj9OLEmp9e/hcfvFMNH4fkzWx68zxgTGgj1My/8b3vtfq3swAePhxjEBtEhhYSEAoLu7mxdxisHBQcyZMwfd3d1QFAVlZWVWd4nSjLEwtWQ8dHV1MRZsItviIXmN9vT0YObMmWn/PBHB0NAQ4+EmYAJoEVXVfy3l8/kyEjTZhufFPhgL13b77bcb54lyW7bGw8yZMzPW32xIjLMBv6MQERER2QwTQCIiIiKbYQJoEY/Hg6eeegoej8fqrkwrPC/2w3/zqfHc2E+2/ZtnW39pjCKcS01ERERkK7wDSERERGQzTACJiIiIbIYJIBEREZHNMAEkIiIishkmgBbYvXs3KioqkJeXh+rqarz33ntWdymtdu3ahfvuuw9erxd+vx+rVq3C6dOnTW2eeOIJKIpiWhYvXmxqEw6HsXnzZhQVFWHGjBl49NFHce7cuUwOhdLATvHAWKCrsSIWMnlNDgwMoL6+Hj6fDz6fD/X19bh06VK6h0hTYAKYYfv370djYyO2b9+O48eP46GHHkJtbS26u7ut7lraHD58GJs2bUJ7eztaWloQi8UQCARw+fJlU7vly5ejt7fXWA4cMP8tzMbGRrzxxhvYt28fjhw5guHhYdTV1SEej2dyOHQT2S0eGAs0FatiIZPX5Nq1a9HZ2Ynm5mY0Nzejs7MT9fX1aR0fXYVQRt1///2yYcMG07q77rpLtm7dalGPMq+vr08AyOHDh41169atk5UrV065z6VLl8Tlcsm+ffuMdZ9//rmoqirNzc3p7C6lkd3jgbFASdMlFtJ1TXZ1dQkAaW9vN9q0tbUJAPn4449v/kDomngHMIMikQg6OjoQCARM6wOBAI4ePWpRrzIvFAoBGPuj50mtra3w+/2orKzE+vXr0dfXZ2zr6OhANBo1nbuysjJUVVXZ6tzlEsYDY4F00ykW0nVNtrW1wefzYdGiRUabxYsXw+fz8bq1CBPADLp48SLi8TiKi4tN64uLixEMBi3qVWaJCJqamvDggw+iqqrKWF9bW4u9e/fi4MGDeOaZZ3Ds2DEsW7YM4XAYABAMBuF2uzFr1izT8ex07nKN3eOBsUBJ0yUW0nlNBoNB+P3+CZ/p9/t53VrEaXUH7EhRFNN7EZmwLlc1NDTgxIkTOHLkiGn9mjVrjK+rqqqwcOFClJeX46233sLq1aunPJ6dzl2usms8MBZoPKtjId3X5GRj4XVrHd4BzKCioiI4HI4JP+309fVN+MkvF23evBlvvvkmDh06hNmzZ1+1bWlpKcrLy/Hpp58CAEpKShCJRDAwMGBqZ5dzl4vsHA+MBUo1HWIh3ddkSUkJzp8/P+FYFy5c4HVrESaAGeR2u1FdXY2WlhbT+paWFixZssSiXqWfiKChoQGvv/46Dh48iIqKimvu09/fj56eHpSWlgIAqqur4XK5TOeut7cXJ0+ezOlzl8vsGA+MBZqMlbGQqWuypqYGoVAIH3zwgdHm/fffRygU4nVrFUumntjYvn37xOVyyUsvvSRdXV3S2NgoM2bMkLNnz1rdtbT58Y9/LD6fT1pbW6W3t9dYRkZGRERkaGhItmzZIkePHpUzZ87IoUOHpKamRm6//XYZHBw0jrNhwwaZPXu2vPPOO/LRRx/JsmXL5N5775VYLGbV0OgG2S0eGAs0FatiIZPX5PLly+Wee+6RtrY2aWtrk/nz50tdXV1ax0dTYwJogV//+tdSXl4ubrdbFixYYJpun4sATLq8/PLLIiIyMjIigUBAvvGNb4jL5ZK5c+fKunXrpLu723Sc0dFRaWhokMLCQsnPz5e6uroJbSj72CkeGAt0NVbEQiavyf7+fnn88cfF6/WK1+uVxx9/XAYGBtI+RpqcIiKS6buORERERGQd1gASERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGSaARERERDbDBJCIiIjIZpgAEhEREdkME0AiIiIim2ECSERERGQzTACJiIiIbIYJIBEREZHNMAEkIiIishkmgEREREQ2wwSQiIiIyGaYABIRERHZDBNAIiIiIpthAngDdu/ejYqKCuTl5aG6uhrvvfee1V0iIiIiuiYmgNdp//79aGxsxPbt23H8+HE89NBDqK2tRXd3t9VdIyIiIroqJoDX6dlnn8WTTz6JH/3oR5g3bx6ee+45zJkzBy+++KLVXSPKON4NJxrDeKBs4LS6A9koEomgo6MDW7duNa0PBAI4evTopPuEw2GEw2HjvaZp+PLLL3HbbbdBUZS09jcbiQiGhoZQVlYGVeXPKdNZ8m747t278cADD+A3v/kNamtr0dXVhblz515zf03T8MUXX8Dr9TIWJsFYyC6Mh/RiPNxEQl/b559/LgDk73//u2n9f//3f0tlZeWk+zz11FMCgMvXXHp6ejLxT0o34P7775cNGzaY1t11112ydevWr7R/T0+P5ddZNiyMhezAeGA8ZAveAbwB4386E5Epf2Lbtm0bmpqajPehUAhz587Fg/hPOOFKaz+zUQxRHMEBeL1eq7tCV3Ez7oaLCAAwFqbAWMgejIf0YzzcPEwAr0NRUREcDgeCwaBpfV9fH4qLiyfdx+PxwOPxTFjvhAtOhUE+gf49kL8CmeYuXryIeDw+4bovLi6eEB9Ju3btws9//vMJ6xkLU2AsZA3GQwYwHm4a/gL9OrjdblRXV6OlpcW0vqWlBUuWLLGoV0TW+bp3w0OhkLH09PRkootEGcN4oGzAO4DXqampCfX19Vi4cCFqamrw29/+Ft3d3diwYYPVXSPKmJt5NxzJ/yAVBRAZe5+UXDf+dfw21QFocf1VNEBJ/JwrGhSHAxKLQXG5IbEoFKcLUBUgHodaUID44CAcd1ci4r8FzuEIlJOfAS4XEI1Ci0ShuJyQcHhsf4cDoklKHxOfl/x8LT55H7/uGAWUBW5qPEwj8aUL4Gj9yOpu0E3GO4DXac2aNXjuuefw9NNP49///d/x7rvv4sCBAygvL7e6a0QZw7vhRGMYD5RNeAfwBmzcuBEbN260uhtEluLdcKIxjAfKFkwAieiGrFmzBv39/Xj66afR29uLqqoq3g0n22I8ULZgAkhEN+ym3A1XVECg19XF42O1ewAUVYHE42M1fBPaaOZtiXo/iKK3URVINAbF7dbbuF16TaDbBagqEI1C8c0EBgcxfMetuPQtJ/IveHDbmQIoLhckHIEiAjU/D/FwGEqeBxjVoDidQFzTu5D8DJcTEo4b/Ziyj8nxpNYQprTR99MAqKwBzDI59dshRUHfgjyUtlrdEbrZWANIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGSaARDQtKA5H4gvVeK+oChRVGZtIkfqqqObtif0VR2KbI7FOVaA4nXpbpz7vTXE5AYcDissJxeMGXC7ILfkAgKHZTgzdHcHgvylQ8vKAgnwoeR5938QDexW3S/9MlwuKQ00sif5O6IdDX5fSR6SMNTmGycavv+WfvCKLKHpsDd0dsbonlAZMAImIiIhshgkgERERkc0wASQiIiKyGSaARDQ9JGvd1LGHNydr/ZJfK45kfVxKLaBRczdWX6eoiv6qKOYaPFfi2fdOp77O6YSSqOPTCtwAgJESwbxvfYErs6OAxw3J9wAet34cl8vYP1lbCIdjbEl8vj4Ovb4PyQWAoqS8TjHGCedB4bdpskDyWlUVzPvWFxZ3htKB31mIiIiIbIYJIBEREZHNMAEkIiIishkmgEREREQ2wwSQiKaFRM25eaIEzOugqsarkjq5QlXGDqCmPEjZ4dDXJx4IDTUxYcTpBJTEJI7ERA4tT5/gEbktjoC/C7OKByEuJ7Q8JyTx4OjkJBLF6dT7lHgANFTV/ABqYGxyiqLoi6qY+m+MKzGO8eNOHT+RlQL+Lqu7QGnABJCIiIjIZpgAEhEREdkME0AiIiIim2ECSETTS2qdnPEgZdVc55eQrK9L3S9Zi5fcL3kc/aHQKfV5auIh0g4VUB3Q3Po2dWYUSwtO45u3fgk4HRCXA3DqNYTiSO2batQWKsk6Q8B4cPPYw55Tvs2OFTqaawLHj9t0PlgLSBlmXKd6LC0tOG1tfygtmAASERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgjeGDZ8lKKZMflGtNfEi9VsdPmpjQVDFNvpjwmpgkojn0dZ68KCpdCkryhyCOxMQPVR1rC+gTR5ITUJITToCxSStXGZvpfeokluS4Uya1XPM8EGVApYvXYS5iAkhmTAKJiIhyHhNA0iXvZBAREVHO4//4pFP0v2mqJJ9lRkRERDnLaXUH7E4tyAdGY1BnzICEw1Dcbkg8DkD/Y/YSiUDxeKBdvjzWxukEVBUSjUHNz0N8cBCOW33QLo9CvWUGJBbT64rcLsjwZaizbkWsNwjnnNmQ0CAU30zA5YRcHsGVe+fC9faH6H29Eifu/xO29C5A1/+bh5jXA+elK1CCF4DCWxH/5F9wfnMutAv9UL23QMIRQDTA6YQ2NAz1lhmI939p9ENx6ZeWRCJQ8vOhDQ1B9Xoho6NQPB4gMUaoqnmMBQX6e1UFwlb9q5ClrqcMQVFNtXTJdYA25WeYagMBo3bP6YyjQHUjX40AigJRAFEUXLNX4z//q9QCGruy9IKmrwLVbXUXKA2YAFpMGxmFqrigXb4MAJBYzNgmydfEusnaxKMR/fVSSH8diEz8jN4gACDWc05fMThobHO93QcAKF11Co/g36H/h3kKpv86L/br+5/tNvUjVTwcNvVDomP9kKEhfb/Ea2r/jTbJMY6M6O8lOqENERER3Rz8FTARERGRzTABJCIiIrIZJoBERERENpNTCeCuXbtw3333wev1wu/3Y9WqVTh9+rSpzRNPPAElUQCeXBYvXmxqEw6HsXnzZhQVFWHGjBl49NFHce7cOVObgYEB1NfXw+fzwefzob6+HpcuXUr3EIlyn8i120zYR4OM30+mmACS+AwRMX+Wpn8dizkwokUwqrkBESgCKF+lT+PbaJPso03epwl9J5pGRrSJteWU/XIqATx8+DA2bdqE9vZ2tLS0IBaLIRAI4PK4SQvLly9Hb2+vsRw4cMC0vbGxEW+88Qb27duHI0eOYHh4GHV1dYgnZ64CWLt2LTo7O9Hc3Izm5mZ0dnaivr4+I+MkIiIiuhE5NQu4ubnZ9P7ll1+G3+9HR0cHHn74YWO9x+NBSUnJpMcIhUJ46aWX8Oqrr+K73/0uAGDPnj2YM2cO3nnnHTzyyCM4deoUmpub0d7ejkWLFgEAfve736GmpganT5/GnXfemaYREhEREd24nLoDOF4opD+SpLCw0LS+tbUVfr8flZWVWL9+Pfr6+oxtHR0diEajCAQCxrqysjJUVVXh6NGjAIC2tjb4fD4j+QOAxYsXw+fzGW3GC4fDGBwcNC1E0927776LFStWoKysDIqi4M9//rNpu4hgx44dKCsrQ35+PpYuXYp//vOf1nSWKM0YD5RLcjYBFBE0NTXhwQcfRFVVlbG+trYWe/fuxcGDB/HMM8/g2LFjWLZsGcKJ59gFg0G43W7MmjXLdLzi4mIEg0Gjjd/vn/CZfr/faDPerl27jHpBn8+HOXPm3KyhEqXN5cuXce+99+KFF16YdPsvf/lLPPvss3jhhRdw7NgxlJSU4Hvf+x6GEs98JMoljAfKJTn1K+BUDQ0NOHHiBI4cOWJav2bNGuPrqqoqLFy4EOXl5XjrrbewevXqKY8nIqan9U/25P7xbVJt27YNTU1NxvvBwUEmgTTt1dbWora2dtJtIoLnnnsO27dvN2LnlVdeQXFxMV577TX813/919f7sJQJEqLJ1f/yhmnyxlUmeyT6aUzimOxV0wBNoMb1deErLnwSFQRHvVDiAiWu6W2SbQEgrqVMJNH3hwP661eZ/JF8nzxG6riT7zUNokkO/5iefTIaD9PIJ1FOUspFOfmtZfPmzXjzzTdx6NAhzJ49+6ptS0tLUV5ejk8//RQAUFJSgkgkgoGBAVO7vr4+FBcXG23Onz8/4VgXLlww2ozn8Xgwc+ZM00KUzc6cOYNgMGgql/B4PPj2t789ZSkEUa5iPFC2yakEUETQ0NCA119/HQcPHkRFRcU19+nv70dPTw9KS0sBANXV1XC5XGhpaTHa9Pb24uTJk1iyZAkAoKamBqFQCB988IHR5v3330coFDLaEOW6ZLnD+B96UsslJsN6WMpFjAfKNjmVAG7atAl79uzBa6+9Bq/Xi2AwiGAwiNHRUQDA8PAwfvrTn6KtrQ1nz55Fa2srVqxYgaKiIjz22GMAAJ/PhyeffBJbtmzB3/72Nxw/fhw//OEPMX/+fGNW8Lx587B8+XKsX78e7e3taG9vx/r161FXV8cZwGQ748serlYKAbAelnIb44GyRU4lgC+++CJCoRCWLl2K0tJSY9m/fz8AwOFw4B//+AdWrlyJyspKrFu3DpWVlWhra4PX6zWO86tf/QqrVq3CD37wAzzwwAMoKCjAX/7yFzgcDqPN3r17MX/+fAQCAQQCAdxzzz149dVXMz5mIqskH6U0/u5GarnEZLZt24ZQKGQsPT095gbJ+rhEbZ6+JOrhxj0wWVJr6BL7SVwbq8vTNOM4IqJvA4B4XF8fj+v1fFocakTfpg260DpyJ85eKgRicSjROBCLA1qiHtDoW2L/xLGRfE5o4gHUY/WBKfV/KbWHool5W+q4TeeD9VfZIG3xYAXjOtVjqXWENzZyUU5NArnW0/Tz8/Px17/+9ZrHycvLw/PPP4/nn39+yjaFhYXYs2fP1+4jUa6oqKhASUkJWlpa8B//8R8AgEgkgsOHD+N//ud/ptzP4/HA4/FkqptEGcF4oGyTUwkgEd1cw8PD+Oyzz4z3Z86cQWdnJwoLCzF37lw0NjZi586duOOOO3DHHXdg586dKCgowNq1ay3sNVF6MB4olzABJKIpffjhh/jOd75jvE8+ymjdunX44x//iJ/97GcYHR3Fxo0bMTAwgEWLFuHtt982lVQQ5QrGA+USJoBENKWlS5detbRCURTs2LEDO3bsyFyniCzCeKBcklOTQIgoe43NjxDTq+nrlIkSkvLQZdMEEWMySFyfmJGcoKEJoOkTNSQW0ydixGJALAbE41CvRAEA7n4H3u67GwPnZ0KJxqBeiUGJ6m0QjRn76xM/NP1zNM2YfGKaaIKxiSqmSR9aykSRxDjGj/taNc1EmfJ2391Wd4HSgAkgERERkc0wASQiIiKyGSaARERERDbDBJCIpofkA481GXuQs2imr5P1dUadnWhjtXbJBzHH43q9XTyeePhzXN8mmlHDh1hMXxeLQaJRSFyDOhIBABQEFZz6VxnyzrmAcATKaBgIR/TjRKPG/tASNYTJWsOUz9fHkVKnqE1S3zfFGCecBxn3YGiiTJCx+tpT/yqzuDOUDkwAiYiIiGyGCSARERGRzTABJCIiIrIZJoBERERENsO/BEJE04I+iUM1Jj0YkzoAKGrKpI/kq2gQTR17b0wGSXkgs2iAouoPbk5O2gAgiQc7668aEI1CGR4FAHjPxRD3uJF/QSBXrkCJxyHhCCQWgxIO6/tHoolJJdGUCSliTD4x9wMT+mi8SmKiiL7SNMbk+Me2E2WYCAAN3i631T2hNOAdQCIiIiKbYQJIREREZDNMAImIiIhshjWARDQ9iAZATTy0Wa89MjbFodfYpT70ObXNJNuMhz8ropfXiQaJ6A97lkhU3x6JAqqi1wOGBgEAt3x6Ce7QLXAORyCXRyAul17rF41BwxV9/yvhiTV6yc9IPGz6K/fRdA5S1ok24TwQZZwI/B9dsboXlAa8A0hERERkM0wAiYiIiGyGCSARERGRzbAG0CLJPwofQxTgY74miCEKYOw8Ue4yYkGi0INBSdS+KeMbmreltkldJ5peNJh4VqDxc65oUESDSAyKKBCJQhEAogAShyoRxCUKiYcRizmBeASKRBLHjEKTGBSRxH7J/TXzNZr8POPz45P38WuMMYaY6TxRbpuO/zfEY1cgErW6GwD4f8PNxATQIv39/QCAIzhgcU+mt6GhIfh8Pqu7QWk0Fgtv6SuS39cn+/4+fptMsk0b95raJpZ4jY57BYDBxOvHiSUptf49PG6/GCYa34/J+vh1xpjAWLCHafl/w3v/a3UPJmA83DgmgBYpLCwEAHR3d/MiTjE4OIg5c+agu7sbiqKgrKzM6i5RmjEWppaMh66uLsaCTWRbPCSv0Z6eHsycOTPtnyciGBoaYjzcBEwALaKq+q+lfD5fRoIm2/C82Adj4dpuv/124zxRbsvWeJg5c2bG+psNiXE24HcUIiIiIpthAkhERERkM0wALeLxePDUU0/B4/FY3ZVphefFfvhvPjWeG/vJtn/zbOsvjVGEc6mJiIiIbIV3AImIiIhshgkgERERkc0wASQiIiKyGSaARERERDbDBNACu3fvRkVFBfLy8lBdXY333nvP6i6l1a5du3DffffB6/XC7/dj1apVOH36tKnNE088AUVRTMvixYtNbcLhMDZv3oyioiLMmDEDjz76KM6dO5fJoVAa2CkeGAt0NVbEQiavyYGBAdTX18Pn88Hn86G+vh6XLl1K9xBpCkwAM2z//v1obGzE9u3bcfz4cTz00EOora1Fd3e31V1Lm8OHD2PTpk1ob29HS0sLYrEYAoEALl++bGq3fPly9Pb2GsuBA+a/hdnY2Ig33ngD+/btw5EjRzA8PIy6ujrE4/FMDoduIrvFA2OBpmJVLGTymly7di06OzvR3NyM5uZmdHZ2or6+Pq3jo6sQyqj7779fNmzYYFp31113ydatWy3qUeb19fUJADl8+LCxbt26dbJy5cop97l06ZK4XC7Zt2+fse7zzz8XVVWlubk5nd2lNLJ7PDAWKGm6xEK6rsmuri4BIO3t7UabtrY2ASAff/zxzR8IXRPvAGZQJBJBR0cHAoGAaX0gEMDRo0ct6lXmhUIhAGN/9DyptbUVfr8flZWVWL9+Pfr6+oxtHR0diEajpnNXVlaGqqoqW527XMJ4YCyQbjrFQrquyba2Nvh8PixatMhos3jxYvh8Pl63FmECmEEXL15EPB5HcXGxaX1xcTGCwaBFvcosEUFTUxMefPBBVFVVGetra2uxd+9eHDx4EM888wyOHTuGZcuWIRwOAwCCwSDcbjdmzZplOp6dzl2usXs8MBYoabrEQjqvyWAwCL/fP+Ez/X4/r1uLOK3ugB0pimJ6LyIT1uWqhoYGnDhxAkeOHDGtX7NmjfF1VVUVFi5ciPLycrz11ltYvXr1lMez07nLVXaNB8YCjWd1LKT7mpxsLLxurcM7gBlUVFQEh8Mx4aedvr6+CT/55aLNmzfjzTffxKFDhzB79uyrti0tLUV5eTk+/fRTAEBJSQkikQgGBgZM7exy7nKRneOBsUCppkMspPuaLCkpwfnz5ycc68KFC7xuLcIEMIPcbjeqq6vR0tJiWt/S0oIlS5ZY1Kv0ExE0NDTg9ddfx8GDB1FRUXHNffr7+9HT04PS0lIAQHV1NVwul+nc9fb24uTJkzl97nKZHeOBsUCTsTIWMnVN1tTUIBQK4YMPPjDavP/++wiFQrxurWLJ1BMb27dvn7hcLnnppZekq6tLGhsbZcaMGXL27Fmru5Y2P/7xj8Xn80lra6v09vYay8jIiIiIDA0NyZYtW+To0aNy5swZOXTokNTU1Mjtt98ug4ODxnE2bNggs2fPlnfeeUc++ugjWbZsmdx7770Si8WsGhrdILvFA2OBpmJVLGTymly+fLncc8890tbWJm1tbTJ//nypq6tL6/hoakwALfDrX/9aysvLxe12y4IFC0zT7XMRgEmXl19+WURERkZGJBAIyDe+8Q1xuVwyd+5cWbdunXR3d5uOMzo6Kg0NDVJYWCj5+flSV1c3oQ1lHzvFA2OBrsaKWMjkNdnf3y+PP/64eL1e8Xq98vjjj8vAwEDax0iTU0REMn3XkYiIiIiswxpAIiIiIpthAkhERERkM0wAiYiIiGyGCSARERGRzTABJCIiIrIZJoBERERENsMEkIiIiMhmmAASERER2QwTQCIiIiKbYQJIREREZDNMAImIiIhshgkgERERkc0wASQiIiKyGSaARERERDbDBJCIiIjIZpgAEhEREdkME0AiIiIim/n/AXi4qKym9bQDAAAAAElFTkSuQmCC' width=640.0/>\n",
       "            </div>\n",
       "        "
      ],
      "text/plain": [
       "Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# not working for 3d FFT\n",
    "Nt, Ny, Nx = 2**4, 2**9, 2**13\n",
    "Nt, Ny, Nx = 2**4, 2**13, 2**9\n",
    "Nt, Ny, Nx = 2**8, 2**8, 2**12\n",
    "\n",
    "# not working for both 3d and 2D+1D\n",
    "Nt, Ny, Nx = 2**8, 2**10, 2**10\n",
    "Nt, Ny, Nx = 2**8, 2**12, 2**8\n",
    "\n",
    "# working\n",
    "Nt, Ny, Nx = 2**4, 2**9, 2**9\n",
    "Nt, Ny, Nx = 2**4, 2**9, 2**10\n",
    "Nt, Ny, Nx = 2**4, 2**9, 2**11\n",
    "Nt, Ny, Nx = 2**4, 2**9, 2**12\n",
    "Nt, Ny, Nx = 2**4, 2**12, 2**9\n",
    "Nt, Ny, Nx = 2**4, 2**10, 2**10\n",
    "Nt, Ny, Nx = 2**8, 2**9, 2**9\n",
    "Nt, Ny, Nx = 2**8, 2**9, 2**10\n",
    "Nt, Ny, Nx = 2**8, 2**8, 2**11\n",
    "Nt, Ny, Nx = 2**4, 2**12, 2**8\n",
    "\n",
    "# eval\n",
    "Nt, Ny, Nx = 2**4, 2**12, 2**8\n",
    "\n",
    "z = get_data('/tmp/Etyx.bin', Nt, Ny, Nx)\n",
    "Z3d = get_data('/tmp/Efkykx3d.bin', Nt, Ny, Nx)\n",
    "Z21d = get_data('/tmp/Efkykx21d.bin', Nt, Ny, Nx)\n",
    "Zp = npf.fftn(z)\n",
    "\n",
    "fig, ((ax01, ax02, ax03),\n",
    "      (ax31, ax32, ax33),\n",
    "      (ax211, ax212, ax213),\n",
    "      (axp1, axp2, axp3)) = plt.subplots(ncols=3, nrows=4)\n",
    "ax01.imshow(abs(z[0]), aspect=Nx/Ny, interpolation='none')\n",
    "ax02.imshow(abs(z[:, 0]), aspect=Nx/Nt, interpolation='none')\n",
    "ax03.imshow(abs(z[:, :, 0]), aspect=Ny/Nt, interpolation='none')\n",
    "ax31.imshow(abs(npf.fftshift(Z3d[0])), aspect=Nx/Ny, interpolation='none')\n",
    "ax32.imshow(abs(npf.fftshift(Z3d[:, 0])), aspect=Nx/Nt, interpolation='none')\n",
    "ax33.imshow(abs(npf.fftshift(Z3d[:, :, 0])), aspect=Ny/Nt, interpolation='none')\n",
    "ax211.imshow(abs(npf.fftshift(Z21d[0])), aspect=Nx/Ny, interpolation='none')\n",
    "ax212.imshow(abs(npf.fftshift(Z21d[:, 0])), aspect=Nx/Nt, interpolation='none')\n",
    "ax213.imshow(abs(npf.fftshift(Z21d[:, :, 0])), aspect=Ny/Nt, interpolation='none')\n",
    "axp1.imshow(abs(npf.fftshift(Zp[0])), aspect=Nx/Ny, interpolation='none')\n",
    "axp2.imshow(abs(npf.fftshift(Zp[:, 0])), aspect=Nx/Nt, interpolation='none')\n",
    "axp3.imshow(abs(npf.fftshift(Zp[:, :, 0])), aspect=Ny/Nt, interpolation='none')\n",
    "fig.tight_layout()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "c3001a05-c308-478a-8f70-2995d743e72d",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.11.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}