{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/users/dae/demillan/anaconda3/lib/python3.6/site-packages/h5py/__init__.py:36: FutureWarning: Conversion of the second argument of issubdtype from `float` to `np.floating` is deprecated. In future, it will be treated as `np.float64 == np.dtype(float).type`.\n",
      "  from ._conv import register_converters as _register_converters\n"
     ]
    }
   ],
   "source": [
    "import h5py\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "from mpl_toolkits.mplot3d import Axes3D\n",
    "import math\n",
    "import gc\n",
    "import datetime\n",
    "import glob\n",
    "\n",
    "now = datetime.datetime.now()\n",
    "mfbins = np.logspace(np.log10(1e9),np.log10(1e16), 60, endpoint=True)\n",
    "vfbins = np.logspace(np.log10(1),np.log10(1000), 60, endpoint=True)\n",
    "vfbins = np.append(vfbins, 2000)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['out_35.list.p0.hdf5',\n",
       " 'out_35.list.p1.hdf5',\n",
       " 'out_35.list.p2.hdf5',\n",
       " 'out_35.list.p3.hdf5',\n",
       " 'out_35.list.p4.hdf5',\n",
       " 'out_35.list.p5.hdf5',\n",
       " 'out_35.list.p6.hdf5',\n",
       " 'out_35.list.p7.hdf5']"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "h5_files = glob.glob('out_35.list*')\n",
    "h5_files"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "out_35.list.p0.hdf5\n",
      "out_35.list.p1.hdf5\n",
      "out_35.list.p2.hdf5\n",
      "out_35.list.p3.hdf5\n",
      "out_35.list.p4.hdf5\n",
      "out_35.list.p5.hdf5\n",
      "out_35.list.p6.hdf5\n",
      "out_35.list.p7.hdf5\n"
     ]
    }
   ],
   "source": [
    "ymm_hdf5_distinct = np.zeros(len(mfbins) - 1) \n",
    "ymm_hdf5_subhalos = np.zeros(len(mfbins) - 1)\n",
    "\n",
    "yvm_hdf5_distinct = np.zeros(len(vfbins) - 1)\n",
    "yvm_hdf5_subhalos = np.zeros(len(vfbins) - 1)\n",
    "\n",
    "filei = 0\n",
    "\n",
    "for file in h5_files:\n",
    "    \n",
    "    print(file)\n",
    "    \n",
    "    f = h5py.File(file, 'r')\n",
    "    \n",
    "    PID = f['PID'].value    \n",
    "    Vmax = f['Vmax'].value    \n",
    "    M200c = f['M200c'].value\n",
    "    \n",
    "    pid_distinct_ix = np.where(PID == -1)\n",
    "    pid_subhalos_ix = np.where(PID > 0)\n",
    "    \n",
    "    M200c_distinct_hdf5 = M200c[pid_distinct_ix[0]]\n",
    "    M200c_subhalos_hdf5 = M200c[pid_subhalos_ix[0]]\n",
    "    \n",
    "    vmax_distinct_hdf5 = Vmax[pid_distinct_ix[0]]\n",
    "    vmax_subhalos_hdf5 = Vmax[pid_subhalos_ix[0]]\n",
    "    \n",
    "    ymm_hdf5_distinct_it = np.histogram(M200c_distinct_hdf5, bins = mfbins)[0]\n",
    "    ymm_hdf5_subhalos_it = np.histogram(M200c_subhalos_hdf5, bins = mfbins)[0]\n",
    "\n",
    "    yvm_hdf5_distinct_it = np.histogram(vmax_distinct_hdf5, bins = vfbins)[0]\n",
    "    yvm_hdf5_subhalos_it = np.histogram(vmax_subhalos_hdf5, bins = vfbins)[0]\n",
    "    \n",
    "    np.savetxt('mfunc_uchuu_results_hdf5_it_' + str(filei) + '_' + now.strftime(format = '%Y_%m_%d_%H_%M') + '.out', (mfbins[:-1], ymm_hdf5_distinct, ymm_hdf5_subhalos))\n",
    "    np.savetxt('vfunc_uchuu_results_hdf5_it_' + str(filei) + '_' + now.strftime(format = '%Y_%m_%d_%H_%M') + '.out', (vfbins[:-1], yvm_hdf5_distinct, yvm_hdf5_subhalos))\n",
    "    \n",
    "    ymm_hdf5_distinct += ymm_hdf5_distinct_it\n",
    "    ymm_hdf5_subhalos += ymm_hdf5_subhalos_it\n",
    "\n",
    "    yvm_hdf5_distinct += yvm_hdf5_distinct_it\n",
    "    yvm_hdf5_subhalos += yvm_hdf5_subhalos_it\n",
    "    \n",
    "    filei += 1\n",
    "    \n",
    "np.savetxt('mfunc_uchuu_results_hdf5_all.out', (mfbins[:-1], ymm_hdf5_distinct, ymm_hdf5_subhalos))\n",
    "np.savetxt('vfunc_uchuu_results_hdf5_all.out', (vfbins[:-1], yvm_hdf5_distinct, yvm_hdf5_subhalos))\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAzoAAANJCAYAAADEB0nbAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAMTQAADE0B0s6tTgAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3X1czff/P/DHqRy6zkXRUaQSmqtczbXlolOtYswUswzDGKLZJldl5GpyMRehrYy5nNHFcGomtrkcB/NlUq5Kp1wroajz+8Ov98epU4qOU6fH/XbrdnPe79f79X6ed++T9/O8rkQHDx5UgoiIiIiISIfoaTsAIiIiIiKiisZEh4iIiIiIdA4THSIiIiIi0jlMdIiIiIiISOcw0SEiIiIiIp3DRIeIiIiIiHQOEx0iIiIiItI5THSI6JUWLlyI+fPnazWGjIwMuLq64ubNm2/93AUFBRg1ahTOnTun8XP5+vrit99+A1D8PSckJODrr7/WeAylWbFiBTZt2iS8dnV1xalTp8pVx8SJExEVFSW8vnjxIkaNGoW+ffti4cKFrzy+oKAA/v7+uHDhQqnlFi5cCFdXV7i6uuLSpUvlirGqevTokfCeJ06cqO1wKkx57xFNCQgIwA8//KC18xNR+RhoOwAi0q6AgACcPXsWAKCnp4c6deqgZ8+eGDt2LMRiMQDgiy++eKNznDlzBlOmTMHvv/8OfX19YXtGRgb8/PywefNmNGzY8I3OoUmJiYkwNjZG69atNX6u8PBwGBoaAgAsLS2xa9cumJubAwB69+6NH3/8Ef/++y9atWpV7Ni8vDx8+OGHGDlyJAYMGFBs/7hx49CiRQtMnjwZqamp2LBhA86fP4+cnBzUq1cP77zzDqZNm4YaNWqojS0zMxMHDhzAli1bKvAdAxs2bICDgwMWLlwIQ0ND4b4oKjY2FiYmJtDT08OQIUPwww8/YOnSpaXW3bNnT0yePFm4hm/b0aNHsW7dOqSnp8POzg4BAQFwdnYusfzt27exatUqnDlzBvn5+WjevDnGjx8Pe3t7ocz9+/cRHh6OY8eOITc3F3Z2dpg/fz7q1q0LY2Nj7Nq1C9u3b39lIliVFL1HNC0/Px99+/bFsmXL0LZtW2H73LlzS/x8EFHlw0SHiPDhhx/Cz88P+fn5uHHjBhYtWgQjIyOMGjUKAGBiYqLlCLUrOjoabm5ub+VcFhYWwr/19fVRp04dlde9e/dGdHS02kRHLBbD1dUV8fHxxRKdGzdu4NKlSwgICMCzZ8/w1VdfoVmzZliwYAHMzMxw8+ZNHDp0CAUFBSXG9ttvv6FTp04Vfj+kp6ejX79+sLS0BPCiVQIAVq9ejQYNGgjljI2NhX/36tULK1euRGpqKmxtbUusWywWq1zDt+nGjRuYM2cOhg8fjh49eiAmJgbffPMNNm3aVGLitWDBAiiVSixduhRisRhRUVGYPn06tm3bBpFIhLy8PEydOhWNGjXCwoULYW5ujmvXrgkP3yKRCHXq1HkrycDbVPQe0RYzMzOtnp+Iyodd14gItWrVQp06dWBpaYn27dujZ8+euHz5srC/aNe1jIwMBAQEwM3NDaNHj8bBgwfh6uqKjIyMN47ljz/+wMiRI+Hm5obBgwerdJMCgKtXr+Lzzz+Hh4cHAgICkJmZKexT163k5a5gZ86cgaurK/Lz84X9+/fvx+DBg0uM5969ezh37hzeffddlW1z5szBwIED4enpicmTJyM5OVnYn5eXh++++w4ffPABpFIpPvnkE/z111/C/kuXLmHq1Klwd3eHj48PZs6cqTZedd31OnfujL///hvPnz9XG69UKsXFixeRlpamsl0mk6Fx48Zo3rw5rl+/joyMDEyZMgXNmjWDtbU1OnTogMDAQNSsWbPEa3Ho0CGV61Do9u3bCAwMhLu7O8aMGYOUlBRhn1KpREREBLy9vTFgwABs27ZN5VhXV1dkZmZi8eLFcHV1xf79+4V95ubmqFOnjvAjEomEfcbGxmjZsiUOHTpUYrzqBAQECF27Xv7RhNjYWDg5OWH48OGws7PDxIkTYWRkhISEhBKPuXjxIgYNGgRHR0c0atQIH3/8MW7duoX79+8DAOLi4vDs2TPMnj0bLVq0gEQiQdeuXd/4Abzws3Hy5EmMGDEC7u7uCA4ORl5eHnbv3o0PP/wQH3zwAbZu3Sock5eXh9DQUHz00Ufw8PDAmDFjcPr0aWG/UqnEhg0bMHjwYLi5uWHo0KGIiYkRji3tM/IydfdIVFRUsa55Rf9O+fr6YseOHQgODoaHhwf8/f1V4gOAf/75B59//jnc3NwwcOBALFu2DAAwbNgwAMCUKVPg6uoqdJcr+jcmLS0N06ZNg1QqxQcffIDw8HCVvy9liYGINIctOkSk4tatWzh16hR69epVYpn58+dDX18fa9aswcOHD7FixYoKOffJkyexYMECfPbZZ+jatSsePnyIa9euqZSJiorChAkTYGFhgcWLF2P16tWYO3duhZxfnfPnz6N27doq3yTn5eWhdevW8Pf3R40aNbBz507MmDEDmzZtglgsxq5du5CUlISFCxfCzMwMN27cEL5xf/DgAQIDA9GrVy+hS+Dx48fLHI+TkxPy8vKQnJyM5s2bF9vv7OwMW1tbxMfHY+TIkQBePHAeOHAA/fv3B/AigdDT08Off/6J999/XyWBKMnDhw9x48YNODk5Fdv3008/Yfz48QgICMDatWuxaNEirF+/HsCLBOvXX3/FV199hcaNGyMyMhIpKSlo3749AGDXrl0YM2YMfH190bt3bxgbGwsP9VOnTsXz58/RpEkTjBw5sliXr2bNmuH8+fNlvnbAi65HhUmiUqnE/PnzS0waAeDrr78udWxWv379MHXqVLX7/vvvP7Rr1054LRKJ4OLigosXL5ZYn7OzMw4ePIj27dvDwMAA8fHxcHBwQO3atQEAx44dQ4sWLbB48WKcOHEClpaWGDZsWKmfV+BFEpCRkYHly5eXWm7Lli2YPn068vLyMHPmTMyYMQOWlpYICwvDhQsXsGjRInTq1AkODg7Iz8+Hra0tPvroIxgaGuL333/HzJkz8fPPP6N27dpITEzEgQMHMHv2bFhaWkKhUCAnJwcASv2MFKXuHnk54SrNtm3bMHbsWIwePRrbtm1DaGgotm7diho1auDatWuYPn06hgwZgm+++Qa5ubnC/bRmzRoMGjQIISEhaNmypdCN92X5+fmYOXMmJBIJ1q5di9u3b2PRokUwMTHBxx9/XKYYiEizmOgQEbZt24ZffvkFBQUFyMvLQ4cOHVT+o37Z1atXcf78eZVxNYMHD0ZYWNgbx7F582Z4eXnho48+AgDY2NjgnXfeUSkzbNgwuLi4AHjR5W7lypVvfN7SZGZmFuv61KBBAwwaNEh4PXnyZHh5eeHixYto06YNbt++DUdHRzRr1gwAYG1tLZTdvXs3rKys8OWXXwoJxsvjL16lZs2aMDY2RmZmptpEBwDc3Nzw22+/4dNPP4VIJMKZM2dw+/Zt9O3bF8CLsT9jx47FypUrsX79ejg7O6N9+/Zwd3eHqamp2jpv3boFAKhbt26xfd7e3ujevTuAF7+fL774Ak+ePIGhoSGio6MxYMAAvPfeewCAadOmqbSgFbbUmJiYCNfZ0NAQX3zxBVq2bIn8/Hzs3bsXAQEBWL9+Pezs7FSOPXLkSJmvHaDa9Wjz5s1ITU3FunXrSiz/5ZdfIi8vr8T9RkZGJe578OCBSldE4EWSmZSUVOIxs2fPxuzZs4UEVCKR4LvvvhPulYyMDJw5cwZeXl5YvHgxzp49i7lz52L58uVquzMWqlu3bqndEguNGTNGuG979eqFQ4cOYf78+RCLxWjUqBG2bNmCs2fPwsHBAYaGhhg+fLhwrL+/Pw4cOIATJ05AKpXi9u3baNiwIVq2bAmRSKTSDbG0z0hR6u6RsurRowekUikA4NNPP8Vvv/2GtLQ0NGnSBFu3bkWHDh0wevRooXxhIl/YtdDMzKzEc546dQoKhQIrV66EmZkZ7O3tMWLECPz4448qfz9Li4GINIuJDhHh/fffx+DBg1FQUACFQoE1a9Zg1apVCAgIKFb25s2bMDY2Vpk8oPBh5U1du3YNAwcOLLXMy0lBnTp1kJWVhfz8fJVJDipSXl5esW9enz9/jsjISPz111+4e/cu8vPzkZubKyQD/fr1w7Rp05CcnIyOHTuiZ8+ewjW6evUq2rRpU6ZWlJLUrFkTubm5Je53c3NDZGQkzp07hzZt2iA+Ph7t2rVTaZX66KOP4O7ujlOnTuH//u//sGPHDuzYsQNr165FvXr11F4HAGq/hS76OwFeDJg3NDREamqqysQCpqamr5x4wtzcXCWRdHZ2RmpqKqKjozF58mSV61BaElKaEydOYPPmzVi2bFmpD89vMiZEqVSW+5gffvgB+vr6WLFiBcRiMXbs2IGZM2dizZo1qFGjBpRKJSwtLTFx4kSIRCI0bdoUZ8+exb59+0pNdD777LMynf/l32Xt2rUhkUhUWjNq166NBw8eCK+3b9+O+Ph43L59G8+ePUNeXh5u374N4MVEEDt27IC/vz/effdddOvWTRjYX9pnpCKVdG82adIEV69eRZ8+fV677hs3bsDGxkYleXZ2dsbDhw+RlZUlbC8tBiLSLI7RISLh4dPW1hadOnXCJ598gtjYWLUP06/z8Fb4rXdht5VChYPOS/tWvCgDg/99P1M0WdDT0ysW38v95QvLv1ymtG5LwItvdIvGvXXrVqFr2Pfff4+IiAiYmJgI52rRogW2bNmCDz/8EJmZmZg4cSK2b99e5vf4KtnZ2aXOImZlZYW2bdsiISEBubm5OHz4sNrJFMzMzODq6oovvvgCkZGReP78uTCGQl1Z4H+/s5ep+528fI3fJKkr5OTkBIVCobItOzv7tcam3Lx5E/PmzcOkSZPQokWLUst+/fXX8PDwKPGntJbMokkB8KILYGE3NHVxxcbG4uuvv0arVq3QrFkzTJ8+HWlpaTh58qRQZ8OGDVWuqa2trZBcvKmiv8uXXxduK/zdJiQk4KeffhJadCMiItC4cWPhM9WgQQNs2rQJo0aNwpMnTzBjxgyhBfZNPyMvx1FI3Wf5Vffm21AZYiCqrtiiQ0TF6Ovro6CgAM+ePSs2ON3GxgY5OTm4efOm8M18aV1xAEAikUAkEiE5OVllzMLly5dhbGwsPPjZ2dnhzJkzrxxvUBJzc3Pcu3dPeP3gwQOV14XdiO7duwcrKysAwJUrV0qt08HBAQqFAnl5ecI32xcuXECvXr2EOG/duoXs7GyV48zMzODm5gY3Nzds3boV+/btw5AhQ9CkSRP89ddfUCqVr5UAKBQK5ObmwsHBodRyUqkUK1euFMa19OjRo9Tyhd2Cnj59qna/RCKBoaEhbty4obb7WklsbW3x33//Ced/9OjRa62FlJKSgkaNGqlsu3bt2iuvQ1FPnjzB7Nmz0adPH3h6er6y/Jt0XWvevDnkcrnKNrlcjg8//FBt+cJrr6f3v+8gC++Rwm5nzs7Oxe6fmzdvCvfz23ThwgW4uLjA3d0dwItrW9iqWcjQ0FD4rLRv3x6LFi3CpEmTAJT8GSkLCwsLYSxXoatXr5arG2iTJk1w5swZtefU09ODnp6eyhclRTVq1AhpaWkqrTcXLlyAhYUFZ2cjqiR0pkUnISEBEyZMwIQJE8q9eB1Rdff06VPcu3cPd+/exb///ovNmzejVatWaqcRbtKkCVq2bImlS5ciJSUFp06dws6dO0ut38TERHjwLuzX/tdff+HHH39U6ar28ccfIy4uDjt37kRaWhouXryIffv2lfl9tG7dGn/++SdOnz6NK1euYMmSJSpdrRo2bIh69eohKioKN2/eREJCAhITE0ut08nJCUZGRioLTkokEpw4cQJJSUnCgOqXu/fs3LkTiYmJSEtLw5UrV/DPP/8IUyB/8MEHyMzMxHfffYcrV67g2rVr5fom+/z582jcuLHa7mUv69GjBwoKCrB27Vr06tULtWrVEvZdunQJs2fPxl9//YXU1FRcv34dERERuHr1Kjp37qy2Pn19fbi4uJR78L+Pjw/27NmDQ4cO4dq1a1iyZInKg7w6MpkMBw8eRGpqKq5evYrVq1fj33//hZeXl0q58+fPC5MalFVYWBhEIhH8/Pxw79494acklpaWaNiwYYk/JbXOAC/GLl26dAk///wzrl+/jlWrVuHx48fo16+fUOaTTz7Bn3/+CeDFg7O1tTWWLl2K5ORk3LhxA2FhYahRo4bQLc3b2xu3b99GREQE0tLS8Ntvv+HIkSN4//33S33fGzZsQGhoaHku1StJJBL83//9H86dO4erV69iwYIFKuOA9u/fj/379+P69etITU3FX3/9JXwOSvuMlEXr1q2Rnp6O6OhopKWlYcOGDeWe9dHPzw///PMPIiIicOPGDSQnJ2P37t0AXiSYlpaWkMvluH//Pp48eVLs+A4dOsDa2hqLFi3C1atXcfz4cURFRal0uyQi7dKJFp1Hjx5hy5YtCA8Px9OnTxEQEICIiAiN9dkn0jW//PILfvnlF4hEItSuXRvt2rXD2LFjSywfFBSEhQsX4vPPP4etrS38/PywaNEitTMTFZo8eTI2btyIZcuW4c6dO6hfvz4GDBig8m1qx44d8dVXX2Hz5s1Yv349ateuLcwUVhbvv/8+Ll++jNmzZ8PY2Bhjx45VmerYwMAAQUFBWLFiBUaNGoWOHTvCz88PO3bsKLFOAwMD9OnTB4mJicLD5vDhw5GWloZJkyahdu3a+Oyzz5CamiocU6tWLWzatAk3b95EzZo14eLiInyLbWFhgaVLl2LNmjUYN24cDA0N0aZNmzK/x8TEROEb9NIUfpO+f//+Yt3W6tevj3r16iEiIgKZmZkwMDBAo0aNMGfOHJUWt6KkUik2btyoMgD9Vdzd3ZGWloYlS5ZAX18fQ4YMKTWxAF5064mKikJmZibEYjHs7e3x3XffwdHRUShz+fJlPHr0CN26dStzLADw77//IjMzs9i3+AcPHixXPWXRqFEjhISEIDw8HBs3boSdnZ2w9k2h1NRUoWtkjRo1sGDBAoSHhyMwMBD5+flwdHRUOaZhw4YIDQ3FmjVr8Msvv6Bhw4aYNWtWqYuQAsDdu3eLtba8KR8fH/z333/45ptvYGRkBF9fX5VWFhMTE/z8889YuXIl9PT04OzsjFmzZgEo/TNSFk2aNMGECROwadMmbNiwAT4+PujevXuZJlwoVLjQakREBHbs2AETExNh0gwAGDt2LNatW4etW7eiX79++Oabb1SO19PTw7x587B8+XKMGzcORkZGkEqlahe7JSLtEB08eFDrHUUPHz6MPXv2ICkpCTk5OcVWTwdeTHn566+/4tGjR2jfvj0CAwOFQX0nT57E4cOHERgYCAD45ptv8Pnnn6Nx48Zv/b0QVUcJCQlYtmwZ4uLiXvltfVWUkZGB8ePHY9OmTSqLVr5tCoUCEyZMwE8//aSVRVzz8/MxatQoTJkypVzJmSYsXrwY1tbWpSZdCxcuRH5+PmbMmPEWI6scoqKicOrUKXz//ffaDoWISGsqxRNJbm4u2rVrV+K3IPv27cOmTZswadIkrFq1Cjk5OQgJCRH2P3z4UGVKVBMTEzx8+FDjcRNVV//88w+OHz8OhUKBo0ePIiIiAn369NHJJAd4Mah6woQJFbIg6pu4c+cOpk2bppUkB3jRfW3atGlqJyR4mwoKCtCwYcMSx7q87ODBg/Dw8FBZAFeXPXr0CB4eHvj555+1HQoRkdZViq5rhf2Vz5w5o3b/7t27MWjQIPTs2RMA8NVXX2HYsGFITk6Go6MjzMzMVAYCP3r0iAMBiTQoLy8P4eHhyMzMhIWFBbp3744xY8ZoOyyNepNpaCtKadMHvy1F1zXSBj09PWHl+tKMHTtWaPF5k2miqxIjIyNEREQAQKldSYmIqoNKkeiUJi8vDykpKSrjBSQSCRo0aIALFy7A0dERzs7OWLt2LfLy8vD06VNkZmaWa1AjEZVP165d0bVrV22HQVSq2rVrlzpZgC7S09N75TpFRETVRaVPdLKyslBQUFDsPysLCwthfQITExP4+voiICAAIpEIEyZMUDsRQUFBAe7evQtDQ8MKWdeBiIiI/icvL++Va1MREb1MqVTi2bNnsLa2rvAu8JU+0SnrolpSqRRSqbTUMnfv3sVHH31UEWERERFREWKxuNR1h4iI1BGLxfjxxx8rvEW60ic65ubm0NPTK7Yw2IMHD4TF/8rK0NAQwIvpPHV1DE9QUFCFr5VQmWKoqLrfpJ7yHlue8mUpW1FlqrLK8P54n79+ed7nZVMZ3l95YsjKyoKtra1O/x9LRBWv8G/Hy+veVZRKn+iIxWI4ODjgzJkzwsJwCoUCGRkZr1w3oKjC7mpmZmY6+0fY29tb6+9NkzFUVN1vUk95jy1P+bKULUsZsVis9ftAk3ifa74e3ufaV1Xvc13+P5aINEcTw0oqxTo6WVlZuHXrFi5duoTvvvsO4eHh0NfXR8OGDWFoaIi9e/di1apVmD59OqytrbFmzRrk5+djxYoV5TpPTk4OvLy88PDhQ/4RJp02depUhIWFaTsMIo3ifV65ZGVlwdzcnP/HElG5FP7t2LFjR4XPkFkpWnSOHDmCRYsWCa/HjRsHAFi2bBnatm0LT09P3L9/H8uXLxcWDP3yyy+1FS5Rpfeq8WpEuoD3ORERlaZStOi8LWzRISIi0gy26BDR69Bki45uLmNORERERETVWrVMdIKCgiCTybQdBhEREVViUVFREIlEwo++vj5sbGzg6+uLy5cvq5QNDg6GnZ1dhZ1bJBJh5syZaveNGDECNjY25aovODgYIpGoUq1zlJaWBmNjY5w6dUrY9t5776F79+6vXaednR2Cg4NVtv34449o2rQpxGKxMGPviBEjVH63hT8BAQEqxw4YMADjx48v07kTExNV6jI0NISNjQ08PT0RERGhdup1kUiEqKioMr+/xMREBAcHo6CgQGX7tWvXIBKJkJiYWOa6yiM4OBh//PGHRurWpEoxRudtCw0NZbM6ERERlcnOnTthY2OD/Px8pKSk4Ntvv0WfPn1w/vx54Xli9OjR6N+/v5YjrVpmzZqF3r17C7PqakJ6ejrGjBmDYcOGITIyErVq1RL2WVpaIiYmRqW8tbW1yuvg4GB07NgRkydPRrNmzcp0zpUrV6Jjx4549uwZ0tPTkZCQgAkTJmDVqlVISEhQ6Z519OhRODg4lPn9JCYmIiQkBDNnzlRZXNPa2hpHjx4t94zEZRUSEoIZM2agd+/eGqlfU6plokNERERUVm3btoWjoyMAoFu3bpBIJOjXrx+OHDkCd3d3AICNjU25W1mqs4yMDGzevBnR0dEaPU9SUhLy8/Ph7+9frKVILBajc+fOpR7ftm1btGvXDsuXL8fatWvLdM4WLVqo1DtkyBCMGjUKvXv3xsiRIxEbGyvse9X5y6pmzZoVVpcuqZZd14iIiIheV2ErzstdkdR1Xbty5Qo8PT1hZGQEKysrBAYGYv369RpZL+T27dsYP348bG1tUbNmTdja2mL48OHIzc1VKZeSkgKpVApjY2PY29sXW6qjpC54I0aMwHvvvSe8LuzWV1RZu/BFRUXB3Nwcbm5uavfLZDK0adMGRkZGcHFxUdtt6vvvv4ednR1q1aqFTp064ciRI8VidnV1BQD06dMHIpGoWLe2svDz88PPP/+MJ0+elPvYQl26dMHnn3+OuLg4pKSkCNuLdl07efIk+vXrh7p168LIyAj29vZC17ng4GCEhIQAAGrUqCF0kQPUd10r7AZYlmt56NAh9OvXD+bm5jA2NkabNm3www8/CDECwPz584Vzvs511AYmOkRERPRWyRVyhCSGQK6QV4m68/Pz8fz5c+Tm5uLixYsICgpCgwYNhIdodfLy8tCvXz+cPXsWa9asQVRUFK5evYp58+ZVWFyFHjx4gK5du2Lnzp0IDAzE3r17sXjxYjx79qzYuJCBAwfCzc0N0dHR8Pb2RkBAABISEio8plfZv38/OnfuDAOD4p2LUlJSMHXqVEybNg2//vor6tWrhwEDBuDevXtCmcjISEyaNAmurq7Ys2cPPvnkEwwZMgT3798XysyaNQsrV64EAKxevRpHjx7F6NGjhf23bt1CvXr1YGBgACcnJyxZsqTY2BcA6NGjB7Kzs4slUuXl6ekJAPj777/V7n/06BGkUin09fURFRWFvXv3Yvbs2cK4qtGjR2PUqFEAgL/++gtHjx7F0aNHSz1nWa5ldHQ0+vTpg+fPn2PdunWIjo7GyJEjcf36dQAQzjFixAjhnC9fx8qMXdeIiIjorZEr5PDa6oX07HSsP70ecX5xcLF2qdR1N2/eXOW1RCJBbGwsTE1NSzxm48aNuHLlCo4fP45OnToBADw8PNC2bVukpqa+cUwvW7ZsGa5cuYLTp0+jTZs2wnY/P79iZb/88kt8+umnAIC+ffviwIED2LFjB/r161ehMZVGqVTixIkTmDJlitr9d+7cweHDh9G0aVMAQLt27WBtbY19+/Zh2LBhKCgowJw5cyCVShEZGSkcZ2lpCV9fX+G1g4MDWrRoAQBwdnZW6drVtm1btG/fHu+88w6ePn2K3bt34+uvv0ZycjLWrVunEk/r1q2hp6eHY8eOoU+fPq/9vm1tbQEACoVC7f7//vsP9+/fx+LFi9G6dWth+4gRIwCodo9899131SaJRb3qWiqVSkyePBnt2rXDH3/8IbTe9O3bV6ij8Lo1bNiwynWPY4sOERERvTUxl2KQnp0OAEjPTkdsUuwrjtB+3bt378bJkydx4sQJ7NmzB87OzvDw8MDFixdLPObo0aNo1KiRkOQAL7oADRo0qEJiell8fDw6deqkkuSU5P3331d53apVqwpPvF7lwYMHePLkSYlrpjRt2lR4MAcAKysrWFlZCXGmpaUhNTUVH330kcpxgwYNKtPDPwAEBARg4sSJ6N27Nzw9PbFhwwZMnjwZGzZsQHJyskrZGjVqwNzcvMQEpayUyhdLV5bUdbFp06awsLDA2LFjsWXLFqSlpb3R+Qrsd3EFAAAgAElEQVTrLO1aJiUl4fr16xg1apRGulRqGxMdIiIiemt8mvlAYioBAEhMJfB28q70dbds2RIdOnRAx44d0b9/f8TExECpVJY6TkGhUMDKyqrY9vr165fpnPr6+sjPz1e7Lz8/X+WB/s6dO2jYsGGZ6q1Tp47K65o1a+Lp06dlOraiFJ6vZs2aavcXjbGwbOFxhQlH0WtpYGCAunXrvnZcQ4cOhVKpxMmTJ4vtMzQ0fKMxOgBw48YNAMVnditkbm6OgwcPokGDBhg7dixsbW3RqlUr/Prrr699zlddyzt37gBAme+fqqZaJjpcR4eIiEg7XKxdEOcXh5D3Qiq025qm636ZoaEh7O3tce7cuRLLWFtb49atW8W2Z2ZmlukcVlZWSE9PV7svPT1d5SG/Xr16uHnzZpnqfZVatWoVm8AA+N8D8cvlABQrW7ScOoXJyMvjacqjMFEoei2fP3+Ou3fvvladAITxOepaNu7du4d69eq9dt0AsHfvXgAvZu4rSdu2bbF79248ePAAR48eRePGjTF48GCcP3/+jc5dksL3VFH3T2VTLROd0NBQSKVSbYdBRERULblYu2B2r9kaSUQ0WXehx48fIyUlpcSuV8CLWbZu3LiBEydOCNuUSiV27dpVpnO4urpCJpMhJydHZXtGRgaOHDmiMhGCm5sbTpw4gbNnz5bznRTXuHFj3Lp1SyVJu3fvXrGB+I0bNwYA/Pvvv8K2Z8+eIT4+/pXnEIvFsLe3x5UrV14rRhsbG9ja2mLHjh0q23ft2vVGC6Ju2bIFIpFIpbsh8KIF6enTp8XGapXH0aNHsW7dOvj4+MDe3v6V5fX19dG5c2eEhoaioKBA6CZZ2Ar2pq1LhZycnGBnZ4cffvhB6FqnjlgsrrBzvk2cjICIiIioFGfOnMGdO3egVCqhUCiwatUq3Lt3DxMnTizxGH9/fyxcuBADBw7E/PnzYWlpiYiICNy/f79MYyFmzZqFmJgYdO3aFYGBgbC1tUVKSgpCQ0NhYWGBgIAAoeyUKVOwZcsW9O3bFzNnzkSrVq1w584dREdHIzw8vNRJE4ry8PCAqakp/Pz8MG3aNGRlZWHx4sXF6ujYsSMcHBzw2WefISQkBAUFBVi9enWxWd5K0rNnTxw/frzMcb1MT08Pc+bMwejRo/Hpp5/C19cXycnJWLBgQZkWhL9+/Tr8/f0xdOhQ2Nvb4/Hjx9i9ezc2btyI8ePHF0tECuPs2bNnmeK7ePEiTExM8Pz5cygUCsTHx2PTpk1wdnYWpmxWJy4uDuvXr8eAAQPQpEkT5OTkYOXKlTA1NUWXLl0AQFgQdOnSpfDw8IC+vj46dOhQprjUEYlEWL58OQYOHIjevXtj3LhxsLS0xMWLF3Hr1i1hOmtnZ2f89ttvcHd3R+3atSGRSCCRSF77vG8LEx0iIiKiUgwePFj4t6WlJVq2bIn9+/eX2jtELBYjPj4eEydOxLhx42BiYoKhQ4fi3XffxcKFC195zubNm+PYsWOYM2cOAgMD8eDBA9SrVw9ubm6YO3cuGjRoIJS1sLDA33//jZkzZ2LhwoW4e/cu6tevj969e0MsFpfrvVpYWCAuLg5TpkzBwIEDYWNjg9mzZ+P333/HtWvXhHIGBgaIjo7GhAkTMHToUNSpUwcBAQHo1q2byrowJRkyZIgw5XaTJk3KFSMAjBo1Co8ePUJYWBi2bt2Kli1bYtu2bfj4449feaypqSksLCwwf/58ZGZmQk9PDy1atMDq1asxbty4YuVjY2OFxK4sJk2aBOBF60vdunXRpk0brF69GsOHDy/199G0aVMYGhri22+/hUKhgKmpKTp27IiEhARhtjUvLy+MHz8ea9aswdy5c6FUKkttiSmL/v37IyEhAd9++60wfbWDg4NKMr1q1SpMmjQJ3t7eyM3NxZw5c6rEWjqigwcPvtnVqUJycnLg5eWFhw8flinjJyIiorLJysqCubk5/499BU9PT+Tm5uLAgQPaDkWrCgoK4OTkBH9/f8yaNUvb4ZToyZMnkEgkCAsLE6blpopV+Ldjx44dpXYHfR1s0SEiIiLSgLCwMJiYmMDR0RHZ2dnYvn079u3bh7i4OG2HpnV6enr49ttvMXnyZAQGBsLIyEjbIakVHh6O+vXrY/jw4doOhV4DEx0iIiIiDahZsybCwsJw48YNKJVKvPPOO/jll1+KrWVTXfn6+iI1NRXXrl0Txp5UNoaGhoiMjCzz+jxUubDrGhEREb0xdl0jotehya5r1XJ6aSIAkCvkCEkMgVwh13YoRERERFTBqmWiwwVDSa6Qw2urF4IPBcNrqxeTHSIiIiIdUy07HIaGhrJZvZqLuRSD9OwXK06nZ6cjNilWo4vLEREREdHbVS1bdIh8mvlAYvpioSuJqQTeTt5ajoiIiIiIKlK1bNEhcrF2QZxfHGKTYuHt5M3WHCIiIiIdw0SHqi0XaxcmOEREREQ6il3XiIiIiIhI5zDRISIiIlIjKioKIpFI+NHX14eNjQ18fX1x+fJllbLBwcGws7OrsHOLRCLMnDlT7b4RI0bAxsamXPUFBwdDJBLh+fPnFRFehUhLS4OxsTFOnTolbHvvvffQvXv3167Tzs4OwcHBKtt+/PFHNG3aFGKxGBYWFgBeXMOXf7eFPwEBASrHDhgwAOPHjy/z+Q8dOgR3d3dIJBLUqlULNjY2cHd3x88//1zu93Lt2jWIRCJERES8sqydnR0+/vjjcp+jtPqKXseqiF3XiCqYXCFHzKUY+DTzYdc4IiIdsHPnTtjY2CA/Px8pKSn49ttv0adPH5w/f16YxXX06NHo37+/liOtWmbNmoXevXujffv2GjtHeno6xowZg2HDhiEyMhK1atUS9llaWiImJkalvLW1tcrr4OBgdOzYEZMnT0azZs1KPdfu3bsxaNAg+Pj4YNWqVahTpw6uX7+OhIQE7N27F8OGDau4N0ZlwkSHqAIVrs+Tnp2O9afXI84vjskOEVEV17ZtWzg6OgIAunXrBolEgn79+uHIkSNwd3cHANjY2JS7laU6y8jIwObNmxEdHa3R8yQlJSE/Px/+/v7FWorEYjE6d+5c6vFt27ZFu3btsHz5cqxdu7bUsmFhYWjbti12794NkUgkbPf390dBQcHrvwl6bey6RlSB1K3PQ0RExcmSNbdwtybrBiC04uTl5Qnb1HVdu3LlCjw9PWFkZAQrKysEBgZi/fr1Kg/BFeX27dsYP348bG1tUbNmTdja2mL48OHIzc1VKZeSkgKpVApjY2PY29tjxYoVKvtL6oI3YsQIvPfee8Lrwm59RZW1C19UVBTMzc3h5uamdr9MJkObNm1gZGQEFxcX/PHHH8XKfP/997Czs0OtWrXQqVMnHDlypFjMrq6uAIA+ffpAJBK9VncsPz8//Pzzz3jy5Emp5e7fv4/69eurvS56ev975C7vtXv27BmmTZsGKysrGBkZwdvbG9evX1cbw+bNm+Hk5AQTExN069YNZ8+eVdkfHR0NqVSKBg0awNjYGK1atcKKFSvKlIjFx8ejS5cuMDQ0hIWFBQYOHFisC6dMJkPXrl1hbm4OExMTNGvWDHPnzn1l3ZpSLROdoKAgyGSa/SNI1RPX5yEiKhtZigYTnQquOz8/H8+fP0dubi4uXryIoKAgNGjQQHiIVicvLw/9+vXD2bNnsWbNGkRFReHq1auYN29ehcYGAA8ePEDXrl2xc+dOBAYGYu/evVi8eDGePXumkowBwMCBA+Hm5obo6Gh4e3sjICAACQkJFR7Tq+zfvx+dO3eGgUHxzkUpKSmYOnUqpk2bhl9//RX16tXDgAEDcO/ePaFMZGQkJk2aBFdXV+zZsweffPIJhgwZgvv37wtlZs2ahZUrVwIAVq9ejaNHj2L06NHC/lu3bqFevXowMDCAk5MTlixZovaBv0ePHsjOzi6WSBXVqVMnxMfHY9asWbhw4QKUSmW5r4s6oaGh+O+//xAZGYnVq1fjn3/+gVQqLTbe6vDhw1ixYgVCQ0Oxbds2PH78GP3791cpV5jo/vjjj/jtt9/w2WefYf78+Zg1a1apMcTHx8PT0xNmZmbYvn071qxZg3PnzqFbt25QKBQAXiT2Pj4+aNKkCbZv346YmBhMnToVOTk5FXIdXke17LoWGhoqfBtDVJG4Pg8Rke5p3ry5ymuJRILY2FiYmpqWeMzGjRtx5coVHD9+HJ06dQIAeHh4oG3btkhNTa3Q+JYtW4YrV67g9OnTaNOmjbDdz8+vWNkvv/wSn376KQCgb9++OHDgAHbs2IF+/fpVaEylUSqVOHHiBKZMmaJ2/507d3D48GE0bdoUANCuXTtYW1tj3759GDZsGAoKCjBnzhxIpVJERkYKx1laWsLX11d47eDggBYtWgAAnJ2dVbqptW3bFu3bt8c777yDp0+fYvfu3fj666+RnJyMdevWqcTTunVr6Onp4dixY+jTp0+J72vhwoVITk7GvHnzMG/ePJiZmaF37974+OOPMWjQoPJfqP/P1NQU0dHRQquQk5MTunfvjp9++gkjR44UymVnZ+Ps2bOoXbs2AKBBgwbo2LEjjh8/jm7dugEApk6dKpRXKpXo3r07ateujUmTJmHevHkltjbOmDED9vb22Lt3L/T19QEAXbp0gZOTE5YuXYrvvvsOp0+fRl5eHtauXSs8Z/fu3fu133dFqJaJDpEmcX0eIiL1ZMkyobVl2bFlwnapgxRSR2mlrXv37t2wsbGBUqlEeno6Vq1aBQ8PDxw+fFh4kC7q6NGjaNSokZDkAC9mUhs0aBDOnTv3RvEUFR8fj06dOqkkOSV5//33VV63atWqwhOvV3nw4AGePHkCS0tLtfubNm0qJDkAYGVlBSsrKyHOtLQ0pKamFuuGNmjQILUtROoUnV3N09MTJiYmWLFiBaZNmyaMyQKAGjVqwNzcXGi5KImVlRUOHz6MkydPYv/+/Th+/DgSEhKwZ88ejB07FuHh4WWKragPP/xQpetbt27dYGNjgyNHjqgkOl26dBGSHODF7xaAyu83IyMDc+fOxd69e3Hz5k2V1p7MzEw0aNCg2PlzcnJw6tQpzJgxQ0hyAKBJkybo1q0bEhMTAbxIHmvUqAFfX1989tln6N69e4m/47eFiQ4RERG9FVJH1aQjTBpWJepu2bKlyoOvm5sbbG1tERwcjO3bt6s9RqFQwMrKqtj2+vXrl+mc+vr6yM/PV7svPz9f5YH+zp07ZUpyAKBOnToqr2vWrImnT5+W6diKUni+mjVrqt1fNMbCsoXHFSYcRa+lgYEB6tat+9pxDR06FMuXL8fJkydVft8AYGho+MoxOoU6duyIjh07AniR1A0ePBjr1q3DF198gZYtW5Y7LnX3TP369XHz5k2Vbep+t8D/rrdSqYSPjw/u37+P2bNno2nTpjA0NMSJEycwYcKEEu+D+/fvQ6lUqk2CGjRogOPHjwMAHB0dIZPJsGDBAvj5+SEvLw+dOnXCokWL0KtXr3K/74pQLcfoUNUiV8gRkhgCuUKu7VCIiIhgaGgIe3v7UltmrK2tcevWrWLbMzMzy3QOKysrpKenq92Xnp6u8vBbr169Yg+9r6tWrVrFJjAAXiRTRcsBKFa2aDl1CpORl8fTlEfhFNBFr+Xz589x9+7d16oTgDA+R133rXv37qFevXrlrtPCwgITJ04EAFy4cAFA+a+dunsmMzMTDRs2LFcsKSkpOHnyJCIiIjBy5Ej06NEDHTp0KDHhLFS7dm2IRCJkZGQU25eRkaGSXLq6uiI+Ph5ZWVn4/fffoaenh/fff79M94UmMNGhSq1wuubgQ8Hw2urFZIeISEdIHd6sO5m26gaAx48fIyUlpdRuOV26dMGNGzdw4sQJYZtSqcSuXbvKdA5XV1fIZLJiA7kzMjJw5MgRlYkQ3NzccOLEiWIzbL2Oxo0b49atWypJ2r1794oNxG/cuDEA4N9//xW2PXv2DPHx8a88h1gshr29Pa5cufJaMdrY2MDW1hY7duxQ2b5r1643WhB1y5YtEIlEKt0NgRctSE+fPi02VqsodYkAAPz3338A/peglffa/fLLLyqTJPz9999IS0tDly5dXvGOVD1+/BiAakuaUqnEDz/8UOpxxsbG6NChA3bu3KnSynj9+nUcOXJEZTa+QmKxGL1798Y333yDnJwcXL16tVyxVhR2XaNKTd10zRz/QkRU9b3puJm3WfeZM2dw584dKJVKKBQKrFq1Cvfu3RO+qVfH398fCxcuxMCBAzF//nxYWloiIiIC9+/fL9P00rNmzUJMTAy6du2KwMBA2NraIiUlBaGhobCwsFAZYzJlyhRs2bIFffv2xcyZM9GqVSvcuXMH0dHRCA8PL3XShKI8PDxgamoKPz8/TJs2DVlZWVi8eHGxOjp27AgHBwd89tlnCAkJQUFBAVavXl1slreS9OzZU+jyVF56enqYM2cORo8ejU8//RS+vr5ITk7GggULyjTZ1PXr1+Hv74+hQ4fC3t4ejx8/xu7du7Fx40aMHz8e9vb2KuUL4+zZs2ep9UqlUtjZ2cHb2xtNmzbFkydPcPjwYSxduhRdunQRJgQo77XLzs7GgAEDMHbsWNy+fRvTp09H06ZN8cknn5TlcgmaN28OOzs7jBs3DiEhIRCJRAgPD1fb8ljUvHnz4OHhAS8vL4wfPx6PHj3CnDlzYGFhgcDAQABAeHg4Dh8+DE9PT9ja2uLOnTtYsGABJBLJa3XZqwhs0aFKjdM1ExGRtg0ePBhdunRB165dMW7cOAAvpkcePHhwiceIxWLEx8ejdevWGDduHPz9/WFra4sJEybA3Nz8leds3rw5jh07hqZNmyIwMBBubm6YNWsWevTogWPHjqmMl7CwsMDff/+NDz74AAsXLoS7uzsCAwNhYGAAsVhcrvdqYWGBuLg4ZGVlYeDAgZg5cyYCAgKKTaVtYGCA6OhomJubY+jQoZg0aRI8PDwwYsSIMp1nyJAhOH/+/Gt/0z9q1CgsX74cf/zxB/r374/IyEhs27ZNZTB+SUxNTWFhYYH58+fDy8sLvr6+OHfuHFavXo3vv/++WPnY2FghOSlNUFAQatasiQULFsDDwwMffPAB9uzZg8DAQMhkMmFCgfJeu+nTp8PR0REjRozA+PHj0a5dO8hkMtSoUePVF+olYrEYMTExwnnHjBkDJycnYQru0ri5uWHv3r3CmKNx48bhnXfewd9//y20VLVp0wY5OTmYPn063Nzc8MUXX6BJkyb4448/YGhoWK5YK4ro4MGDFTPJdxWQk5MDLy8vPHz4kNNLVyFyhZzTNRMRVXJZWVkwNzfn/7Gv4OnpidzcXBw4cEDboWhVQUEBnJyc4O/v/8o1XLTpyZMnkEgkCAsLE6blpopV+Ldjx44dFT5LG7uuUaXH6ZqJiKgqCgsLg4mJCRwdHZGdnY3t27dj3759iIuL03ZoWqenp4dvv/0WkydPRmBgIIyMjLQdklrh4eGoX78+hg8fru1Q6DUw0SEiIiLSgJo1ayIsLAw3btyAUqnEO++8g19++aXYWjbVla+vL1JTU3Ht2jU4OztrOxy1DA0NERkZWeb1eahyYdc1IiIiemPsukZEr0OTXdc4GQEREREREemcapnoBAUFQSaTaTsMIiIiIiLSkGrZ4TA0NJTN6kREREREOqxatugQEREREZFuq5YtOkRERKQZWVlZ2g6BiKoQTf7NYKJDREREb0wsFqNBgwawtbXVdihEVMWIxWKNTOHNRIeIiIjeWK1atXD16lXk5eVpOxQiqkKysrJga2sLsVhc4XUz0SEiIqIKUatWLdSqVUvbYRARAeBkBEREREREpIOY6BBVEXKFHCGJIZAr5NoOhYiIiKjSY6JDVAXIFXJ4bfVC8KFgeG31YrJDRERE9ApMdIiqgJhLMUjPTgcApGenIzYpVssREREREVVuTHSIqgCfZj6QmEoAABJTCbydvLUcEREREVHlxlnXiKoAF2sXxPnFITYpFt5O3nCxdtF2SERERESVGhMdoirCxdqFCQ4RERFRGbHrGhERERER6ZxqmegEBQVBJpNpOwwiIiIiompLJpMhKChIY/WLDh48qNRY7ZVMTk4OvLy88PDhQ5iZmWk7HCIiIiKiai0rKwvm5uaIi4uDsbFxhdZdLVt0iIiIiIhItzHRISIiIiIincNEh6iakyvkCEkMgVwh13YoRERERBWGiQ5RNSZXyOG11QvBh4LhtdWLyQ4RERHpDCY6RNVYzKUYpGenAwDSs9MRmxSr5YiIiIiIKgYTHaJqzKeZDySmEgCAxFQCbydvLUdEREREVDEMtB0AEWmPi7UL4vziEJsUC28nb7hYu2g7JCIiIqIKwUSHqJpzsXZhgkNEREQ6h13XiIiIiIhI5zDRISIiIiIincNEh4iIiIiIdA4THSIiIiIi0jlMdIiIiIiISOcw0SEiIiIiIp3DRIeINEaukCMkMQRyhVzboRAREVE1w0SHiDRCrpDDa6sXgg8Fw2urF5MdIiIiequY6BCRRsRcikF6djoAID07HbFJsVqOiIiIiKqTapnoBAUFQSaTaTsMIp3m08wHElMJAEBiKoG3k7eWIyIiIqLKRCaTISgoSGP1iw4ePKjUWO2VTE5ODry8vPDw4UOYmZlpOxwinSdXyBGbFAtvJ2+4WLtoOxwiIiKqZLKysmBubo64uDgYGxtXaN0GFVobEdFLXKxdmOAQERGRVlTLrmtERERERKTbmOgQEREREZHOYaJDREREREQ6h4kOERERERHpHCY6RERERESkc5joEBERERGRzmGiQ0REREREOoeJDhERERER6RwmOkREREREpHOY6BBRlSNXyBGSGAK5Qq7tUIiIiKiSYqJDRFWKXCGH11YvBB8KhtdWLyY7REREpBYTHSKqUmIuxSA9Ox0AkJ6djtikWC1HRERERJUREx0iqlJ8mvlAYioBAEhMJfB28tZyRERERFQZGWg7ACKi8nCxdkGcXxxik2Lh7eQNF2sXbYdERERElRATHSKqclysXZjgEBERUanYdY2IiIiIiHQOEx0iIiIiItI5THSIiIiIiEjnMNEhIiIiIiKdw0SHiIiIiIh0DhMdIiIiIiLSOUx0iIiIiIhI5zDRqYbkCjlCEkMgV8i1HQoRERERkUYw0alm5Ao5vLZ6IfhQMLy2elVossMEioiIiIgqi2qZ6AQFBUEmk2k7DK2IuRSD9Ox0AEB6djpik2IrpF5NJlBEREREpHtkMhmCgoI0Vn+1THRCQ0MhlUq1HYZW+DTzgcRUAgCQmErg7eRdIfVqKoEiIiIiIt0klUoRGhqqsfqrZaJTnblYuyDOLw4h74Ugzi8OLtYuFVKvphIoIiIiIqLXYaDtAOjtc7F2qbAE5+U64/ziEJsUC28n7wqvn4iIiIioPJjoUIXRRAJF9DbJFXLEXIqBTzMf3stERERVHLuuERGBE2oQERHpGiY6RETghBpERES6hokOERE4oQYREZGu4RgdIiJwQg0iIiJdw0SHiOj/44QaREREuoNd14iIiIiISOcw0SEiIiIiIp3DRIeIiIiIiHQOEx0iIiIiItI5THSIiIiIiEjnMNEhIiIiIiKdw0SHiIiIiIh0DhMdIiINkyvkCEkMgVwh13YoRERE1QYTHSIiDZIr5PDa6oXgQ8Hw2urFZIeIiOgtYaJDRKRBMZdikJ6dDgBIz05HbFKsliMiIiKqHpjoEBFpkE8zH0hMJQAAiakE3k7eWo6IiIioejDQdgBERLrMxdoFcX5xiE2KhbeTN1ysXbQdEhERUbXARIeISMNcrF2Y4BAREb1l7LpGREREREQ6h4kOERERERHpHCY6RERERESkc5joEBERERGRzmGiQ0REREREOoeJDhERERER6RwmOkREREREpHOY6BARERERkc5hokNERERERDqHiQ4REREREekcJjpERERERKRzmOgQEVVRcoUcIYkhkCvk2g6FiIio0mGiQ0RUBckVcnht9ULwoWB4bfViskNERFQEEx0ioioo5lIM0rPTAQDp2emITYrVckRERESVCxMdIqIqyKeZDySmEgCAxFQCbydvLUdERERUuRhoO4CKUFBQgICAAFy/fh0+Pj4YNWqUtkMiItIoF2sXxPnFITYpFt5O3nCxdtF2SERERJWKTiQ6enp6mDlzJk6fPo2bN29qOxwiorfCxdqFCQ4REVEJtJboHD58GHv27EFSUhJycnLw+++/Q19fX6XMli1b8Ouvv+LRo0do3749AgMDUadOHbX1WVlZvY2wiYiIiIioCtDaGJ3c3Fy0a9cOfn5+avfv27cPmzZtwqRJk7Bq1Srk5OQgJCTkLUdJRERERERVkdZadPr16wcAOHPmjNr9u3fvxqBBg9CzZ08AwFdffYVhw4YhOTkZjo6OiI6OxsGDB2FmZoa5c+e+tbhJt8iSZZA6SrUdBhERERFVsEo5RicvLw8pKSkYO3assE0ikaBBgwa4cOECHB0d0b9/f/Tv31+LUZIukKUw0SEiIiLSRZUy0cnKykJBQQFq166tst3CwgIPHjxQe8z8+fNx+fJl5Obm4sKFC1iyZAn09Dh7NhERERFRdVQpEx2lUlnuY2bMmKGBSEgXyZJlkKXIAADLji0TtksdpGzdISIiItIRlTLRMTc3h56eHu7fv6+y/cGDB7CwsHjj+oOCgiAWiwEAUqkUUikfbqsTqaNqQhMmDdNiNERERETVi0wmg0z24kvnvLw8jZ2nUiY6YrEYDg4OOHPmDNq3bw8AUCgUyMjIgLOz8xvXHxoaCjMzszeuh6gknOSAiIiISL2XGxqysrKwevVqjZxHa4NYsrKykJycLCzwmZycjOTkZDx58gQAMGDAAOzatQt//vknkpOTsWTJErRu3RqOjo7aCpm0SJYs00i9UgfNJCOFXeOIiIiISDu01qJz5MgRLFq0SHg9btw4AMCyZRVkgYkAACAASURBVMvQtm1beHp64v79+1i+fLmwYOiXX36prXBJyzQ1OxpbXYiIiIh0k9YSHXd3d7i7u5daZtiwYRg2bNhbiojozXCSA9IVcoUcMZdi4NPMBy7WLtoOh4iI6LVUyjE6REDVSxw4yQHpArlCDq+tXkjPTsf60+sR5xfHZIeIiKokJjpUaTFxIHr7Yi7FID07HQCQnp2O2KRYJjpERFQlVcsVNYOCgoQp7Yg0QVOTHBBpmk8zH0hMJQAAiakE3k7eWo6IiIh0lUwmQ1BQkMbqFx08eLD8q3NWUTk5OfDy8sLDhw85vXQVw+maid4euUKO2KRYeDt5szWHiIg0KisrC+bm5oiLi4OxsXGF1s2ua1QlMMkhentcrF2Y4BARUZVXLbuuERERERGRbmOiQ0REREREOoeJDhERERER6RwmOkREREREpHOY6BARERERkc6plokO19EhIiIiItIurqNTgbiODhERERFR5aHJdXSqZYsOERERERHpNiY6RERERESkc5joEBHRWyFXyBGSGAK5Qq7tUIiIqBpgokNERBonV8jhtdULwYeC4bXVi8kOERFpHBMdoipElszZAqlqirkUg/TsdABAenY6YpNitRwRERHpOiY6RFWILIWJDlVNPs18IDGVAAAkphJ4O3lrOSIiItJ1BtoOgIiIdJ+LtQvi/OIQmxQLbydvuFi7aDskIiLScdUy0QkKCoK3tzekUqm2QyF6JVmyTGjJWXZsmbBd6iCF1JH3MFUdLtYuTHCIiEggk8kQG6u5rsxcMJSoCpkqm4owaZi2wyAiIiKqEFwwlIiIiIiIqByY6BBVIVIHdlUjIiIiKgsmOkRViKbG5HDaaiIiItI1THSIiNNWExERkc5hokNERERERDqnWk4vTUSctpqIiIh0GxMdompK+v/Yu//4uuv6XuDvNLSCaZPaXYSEAo8tsXXFIacVtPde8YqU4+QkBcvkFthEneBF1A24wM6GUjbLusduAtriXR6bMuBBp6KU5Oz6+AoYZFeobjZwh2h704EKJxQmJUkztEh7/+BBLllrfrTn5Jx8z/P5ePh4NN9z+jmvPh7U9JXPr7bxhcax1QBAmli6BgAApE5NFp18Ph9JYvM1vMqx1QDATEuSJPL5fNnGr+vr69tfttGrzOjoaORyuRgaGorGxsZKxwGgBPoH+6Nne090LO2ITHOm0nEAmIbh4eFoamqKQqEQDQ0NJR27Jmd0AEiH/sH+yG3OxfXfvj5ym3PRP9hf6UgAVAlFBygrl5FSTj3be6I4UoyIiOJIMXp39FY4EQDVQtEBysplpJRTx9KOaFnQEhERLQtaon1Je4UTAVAtHC8NwKyVac5EYW0henf0RvuSdnt0ABij6AAl5zJSZlKmOaPgAHAARQcoOZeRAgCVZo9OFesf7I91D6xzihAAAEyTolOlHJlKWriMFACoBEWnSjkylbSwJwcAqARFp0o5MhUAAA5dTR5GkM/no729PbLZ0vykuX+wP3q290TH0o6SnfzjyFQAANIsSZLo7S3fqqW6vr6+/WUbvcqMjo5GLpeLoaGhaGxsLMmYr+6lKY4Uo2VBSxTWFpQSAACYguHh4WhqaopCoRANDQ0lHdvStcNkLw0AAFSfmiw6jz7zaMnGspcGAACqT00WnQ/c9YGSHdf86l6adf9lnWVrMIOSgaTSEQCAKlaTReeZkWdKusQs05yJT7/r00oOzKBkp6JDebm0GWB2q8mic+yCYy0xA+BXcmkzwOxXk8dLf+W8r5h9gVkoGUjGZnK6tnaNPc+2Zl1MSkkd7KAZ3zcAZpeaLDpvPfatlY4AHIJs2/hC05ntrGAa0qxjaUd0b+seuzrAKgCA2acmiw4ATMSlzQCzn6IDzErZVkvVKK9Mc0bBAZjFavIwAmD2sycHAJiIogMAAKSOolOjynXZokscAQCoBopOjSrXZYsucQQAoBrUZNHJ5/ORJP5BDgAAlZIkSeTz+bKNX5Onrq1fvz4aGxsrHWPGleuyRZc4AgAwXdlsNlauXBmbNm0qy/g1WXRqVbkuW3SJIwAA1aYml64BAADppujUqHJdtugSRwAAqoGiU6PKtXfGnhwAAKqBogMAAKSOogPwGi69pdz6B/tj3QPron+wv9JRAFJN0QF4DZfeUk79g/2R25yL6799feQ255QdgDJSdABghvRs74niSDEiIoojxejd0VvhRADp5R4doOa59JaZ0rG0I7q3dUdxpBgtC1qifUl7pSMBpJaiA9Q8l94yUzLNmSisLUTvjt5oX9IemeZMpSMBpJaiAwAzKNOcUXAAZoA9OgCv4dJbAEgHRQfgNezJAYB0UHQAAIDUUXQAAIDUUXQAAIDUUXQAAIDUqcmik8/nI0mSSscAAICalSRJ5PP5so1f19fXt79so1eZ0dHRyOVyMTQ0FI2NjZWOA9SQZCBxohsA/DvDw8PR1NQUhUIhGhoaSjp2Tc7oAMy0ZKdZZACYSYoOAACQOkdUOgBAWiUDydhMTtfWrrHn2dasZWyUXP9gf/Rs74mOpR2Rac5UOg5AxSk6AGWSbRtfaDqznRVMQ5r1D/ZHbnMuiiPF6N7WHYW1BWUHqHmWrgHALNezvSeKI8WIiCiOFKN3R2+FEwFUnqIDMAOyrZaqUT4dSzuiZUFLRES0LGiJ9iXtFU4EUHmWrgHMAHtyKKdMcyYKawvRu6M32pe0W7YGEIoOAKRCpjmj4AC8hqVrAABA6ig6AABA6ig6AABA6ig6AABA6ig6AABA6ig6AABA6ig6AABA6ig6AABA6ig6AABA6ig6ALNYMpBUOgIAVCVFB2AWS3YqOgBwMIoOAPAr9Q/2x7oH1kX/YH+lowBMyxGVDgDA9CQDydhMTtfWrrHn2dZsZNuylYpFCvUP9kducy6KI8Xo3tYdhbWFyDRnKh0LYEoUHYBZJts2vtB0ZjsrmIY069neE8WRYkREFEeK0bujV9EBZg1L1wCAg+pY2hEtC1oiIqJlQUu0L2mvcCKAqTOjAzCLZVstVaN8Ms2ZKKwtRO+O3mhf0m42B5hVFB2AWcyeHMot05xRcIBZqSaXruXz+UgSR7ICAEClJEkS+Xy+bOPX9fX17S/b6FVmdHQ0crlcDA0NRWNjY6XjAABATRseHo6mpqYoFArR0NBQ0rFrckYHAABIN0UHAABIHUUHAABIHUUHAABIHUUHAABIHUUHAABIHUUHAABIHUUHAABIHUUHgAMkA0mlIwDAYVF0ADhAslPRobz6B/tj3QPron+wv9JRgJRSdACAGdU/2B+5zbm4/tvXR25zTtkByuKISgcAoDokA8nYTE7X1q6x59nWbGTbspWKRQr1bO+J4kgxIiKKI8Xo3dEbmeZMhVMBaaPoABAREdm28YWmM9tZwTSkWcfSjuje1h3FkWK0LGiJ9iXtlY4EpJCiAwDMqExzJgprC9G7ozfal7SbzQHKQtEB4ADZVkvVKK9Mc0bBAcrKYQQAHMCeHABmO0UHAABIHUUHAABIHUUHAABIHUUHAABIHUUHAABIHUUHAABIHUUHAABIHUUHAABIHUUHgBmTDCSVjgBAjVB0AJgxyU5FB4CZoegAAKnRP9gf6x5YF/2D/ZWOAlTYEZUOAEC6JQPJ2ExO19ausefZ1mxk27KVikUK9Q/2R25zLoojxeje1h2FtYXINGcqHQuoEEUHgLLKto0vNJ3ZzgqmIc16tvdEcaQYERHFkWL07uhVdKCGWboGAKRCx9KOaFnQEhERLQtaon1Je4UTAZVkRgeAGZNttVSN8sk0Z6KwthC9O3qjfUm72RyocYoOADPGnhzKLdOcUXCAiLB0DQAASCFFBwAASB1FBwAASJ0J9+icf/750xqsrq4ubr755jjmmGMOKxQAAMDhmLDoPPfcc/Hxj388GhoapjTYTTfdFC+//HJJggEAAByqCYtOQ0NDnHnmmdHU1DSlwb74xS/GnDkzvxpu+/bt8fnPfz7q6+vj5ZdfjssuuyyWLVs24zkAAIDqMGHR6e3tndZgX/nKVw4rzKH6D//hP8SGDRuioaEhnnjiifiLv/iL+MIXvlCRLAAAQOVV7B6dBx98MLZs2RI7duyI0dHRuO+++6K+vn7ce+688874+te/Hnv27IkVK1bElVdeGYsWLTpgrF/7tV8b+/W8efMqMqsEAABUjyk1gsHBwXjggQfimWeeiYiIYrEYX/jCF6K7uzt+9KMfHdIH/+IXv4jly5fH2rVrD/r6N77xjbj99tvjk5/8ZGzcuDFGR0dj3bp1E4750ksvRWdnZ/ze7/3eIWUCAADSYdIZnX/8x3+MP/mTP4kjjzwy9u7dGzfccEPceOON0draGvv27Yu77rorNmzYEJnM9G4hXrVqVUREPPLIIwd9/e677441a9bE6aefHhERV199dVx44YUxMDAQbW1tcc8990RfX180NjbGDTfcEC+//HL86Z/+abz73e+Ot7/97dPKUs2SgcRN4gAAME2TFp0vfelLccEFF8QHP/jBeOCBB+KGG26I97///fGhD31o7PXbb7992kVnInv37o2dO3fGpZdeOvaspaUljj322Hj88cejra0tVq9eHatXr46IiH379sX69evjLW95S+RyuZLlqAbJTkUHAACma9Kla08++eTY7Mu73vWu+PnPfx7vfOc7x14/44wz4oknnihpqOHh4di3b1+84Q1vGPd84cKF8cILLxzw/gceeCC+853vxEMPPRR/8Ad/EH/8x39c0jwAVLdkIKl0BFKuf7A/1j2wLvoH+ysdBZiiSWd0Xj2yOeKVC0Hnzp0br3/968denzdvXrz44oslDbV///5pvf+MM86IM844o6QZKikZSCLZ+co37a6tXWPPs61ZszsAB2H2m3LqH+yP3OZcFEeK0b2tOwprC5FpLt1KFqA8Ji06ixcvjh//+Mdx/PHHR0TEV7/61XFF56mnnoo3vvGNJQ3V1NQUc+bMid27d497/sILL8TChQsPe/x8Ph/z5s2LiIhsNhvZ7OF/cyzlXpps2/hC05ntLMm4AMD09WzvieJIMSIiiiPF6N3Rq+jAYUiSJJLklR/q7927t2yfM2nRufjii2P+/PljXzc0NIx7/Sc/+Um8733vK2moefPmRWtrazzyyCOxYsWKiHjl5LdnnnmmJBeBrl+/PhobGw97nNfy00SAmWX2m5nSsbQjurd1R3GkGC0LWqJ9SXulI8Gs9tqJhuHh4di0aVNZPmfSojPZCWZr1qw5pA8eHh6OZ599Np5++umIiBgYGIj6+vo47rjj4qijjopzzjknNm7cGEuWLInm5ua45ZZb4uSTT462trZD+rzZKtvqmzXAwZj9ZqZkmjNRWFuI3h290b6k3WwOzBIVuzD0oYceig0bNox9/bGPfSwiIrq6uuKUU06J973vfbF79+646aabxi4MveqqqyoV96Bm4qeJfioJAJWXac4oODDL1PX19U1p5/9dd90V5513XrnzlNXo6Gjkcrm47GuXxesaXlfSQnJFcoWfJgJUiDvHAGan4eHhaGpqikKhcMAWmcM1pRmdW265Jf7pn/5p1hedV9145o0l36MDQOUoOQD8exMWnV/+8pexfv36eOKJJ6Kz02zFROylAQCA6jHhhaH//b//9/iXf/mX+B//438ccHnnbJbP58eOtCsVP00EAICpS5Ik8vl82cafsOg8+uijcfbZZ8eiRYvKFqAS1q9fX5K7cwAAgEOTzWZj/fr1ZRt/wqJzzTXXxN/8zd/EvffeW7YAAAAApTbhHp1sNhu/9mu/FuvWrYs5c+bEe97znpnKBQAAcMgmnNGJiHjb294WnZ2d0d3dPRN5AAAADtukRSci4k1velN87nOfK3cWAACAkphS0YmIOOaYY8qZAwAAoGSmXHQGBwfjqaeeOuD5U089Fc8880xJQ5VbOY6XBgAApq6ix0u/1o033hg//OEPD3j+ox/9KG688caShio3x0sDAEBlVfR46dcaGBiIZcuWHfB82bJlsXPnzpKGAoBqkAyY/af8+gf7Y90D66J/sL/SUSBVplx05s6dG0NDQwc8f/7556Ourq6koQCgGiQ7FR3Kq3+wP3Kbc3H9t6+P3OacsgMlNOWic9ppp8Utt9wSzz///Niz559/Pv7n//yf8fa3v70s4QAA0qxne08UR4oREVEcKUbvjt4KJ4L0mPDC0Nf6+Mc/Htddd12sXbs2Fi9eHBGvHESwdOnS+MQnPlG2gAAwk5KBZGwmp2tr19jzbGs2sm32d1JaHUs7ontbdxRHitGyoCXal7RXOhKkxpSLzsKFC+Pzn/989Pf3x5NPPhn79++P3/iN34hTTjmlnPkAYEZl28YXms5sZwXTkHaZ5kwU1haid0dvtC9pj0xzptKRIDWmXHRelclkIpPxlxAAoBQyzRkFB8pgWkUnSZLYsmVL/PSnP42IiOOPPz7OOeccRzUDkErZVt/fAGarKRedv/7rv44tW7bE+9///vjd3/3diHjlDp3Pf/7z8dOf/jR+//d/v2whSy2fz0d7e7uCBsCE7MkBKJ8kSaK3t3wHcNT19fXtn8obOzo64qqrrorTTz993PMHH3ww/vIv/zJ6enrKErCURkdHI5fLxdDQUDQ2NlY6DgAA1LTh4eFoamqKQqEQDQ0NJR17ysdLH3HEEXHiiSce8PzEE0+MOXOmPAwAAEDZTbmhnH/++fHFL34xRkdHx56Njo7Gl770pTj//PPLEg4AAOBQTHmPzne+853YuXNnnHfeeXHcccdFXV1dPPXUU1FfXx+7d++OrVu3jr335ptvLktYAACAqZhy0VmxYkWsWLGinFkAAABKYspF54Mf/GA5cwAAAJSMUwQAAIDUmXRGZ6oHDXz5y18+7DAAAAClMGnRee655+KNb3xjnHnmmbF48eKZyAQAAHBYJi06N998c3zzm9+Mnp6eOP744+Oss86KM844IxYsWDAT+QAAAKat/uKLL75+ojccc8wx8R//43+MNWvWxIIFC+KBBx6IjRs3xo4dO2Lu3LlxwgknzEzSEnjppZfizjvvjJGRkaivr4+2trZKRwIAKIv+wf7o/n53NL2uKZoXNFc6DhwgSZL4whe+EP/4j/8YF1xwQcybN6+k49f19fXtn+5v2rVrV6xfvz4ee+yxuPvuu6OxsbGkocpldHQ0crlcDA0NzZrMAADT1T/YH7nNuSiOFKNlQUsU1hYi05ypdCw4wPDwcDQ1NUWhUIiGhoaSjj2tU9cef/zxuOmmm+KSSy6Jn//853HZZZfF/PnzSxoIANIuGUgqHYGU69neE8WRYkREFEeK0bujt8KJYOZNukdncHAw7r333rj33nvjpZdeive85z3xuc99Lk488cSZyAcAqZPsTCLblq10DFKsY2lHdG/rHpvRaV/SXulIMOMmLToXXXRRHH300fGe97wnli9fHnV1dfGzn/0sfvazn4173/Lly8sWEgCAqcs0Z6KwthC9O3qjfUm7ZWvUpEmLzv79++PZZ5+NzZs3x+bNmw/6nrq6urj//vtLHg4A0iIZSCLZ+cqSta6tXWPPs61ZszuURaY5o+BQ0yYtOt/61rdmIgcApFq2bXyh6cx2VjANQPpNehhBV1dXPPzww7F3796ZyAMAAHDYJi06CxcujFtvvTVWr14d1157bWzZsiV27do1E9kAIJWyrZaqAZTblO/Ref755+Phhx+OrVu3xrZt2+KYY46JlStXxjve8Y54y1veEnV1deXOetjcowMAANWjnPfoTLpH51WLFi2Ks88+O84+++x46aWX4pFHHomtW7fGn//5n8eePXvi1FNPjfPOOy/e/OY3lzQgAADAdE256LzW3Llz49RTT41TTz01PvGJT8SPf/zj2Lp1a+zatUvRAQAAKu6Qik5ExL/927/FnDlz4sgjj4wTTzxxVl0gms/no729PbJZa6QBAKASkiSJ3t7eso0/6R6dF198MTZt2hSPPfZYnHzyyfGpT30q/uqv/iq+9rWvRUTEf/pP/ymuvfbaeP3rX1+2kKVijw4AAFSPcu7RmfTUte7u7njsscdi9erV8dRTT8VnPvOZ+P73vx8333xzdHV1xdNPPx133HFHSUMBAAAcjkmLzne+8534wz/8wzj33HPjT/7kT+Khhx6Kj370o/GWt7wlTj755Lj00kvjH/7hH2YiKwAAwJRMWnSGhoaiubk5Il45ee11r3tdHH/88WOvn3jiifHcc8+VLyEAAMA0TVp03vCGN8TPfvazsa/PPffcWLBgwdjXIyMjceSRR5YnHQAAwCGYtOgsW7Ystm/fPvb1JZdcMm4j/yOPPBJvetObypMOAICq0j/YH+seWBf9g/2VjgITmvR46U9/+tMTvn7qqafGf/7P/7lkgQAAqE79g/2R25yL4kgxurd1R2FtITLNmUrHgoOadEZnMieeeGIce+yxpcgCAEAV69neE8WRYkREFEeK0bujfHegwOGactEZGhoqZw4A4DAlA0mlI5ByHUs7omVBS0REtCxoifYl7RVOBL/alIrO008/HZ/4xCfKnQUAOAzJTkWH8so0Z6KwthDr/ss6y9aoepPu0dm+fXvk8/m48MILZyIPAABVLNOcUXCYFSYsOg8//HD82Z/9WXz4wx+O97///TOVCQCYomQgGZvJ6draNfY825qNbFu2UrEAKm7CovPpT386Lr744lizZs1M5QEApiHbNr7QdGY7K5gGoHpMuEfnxBNPjH/4h3+If/u3f5upPAAAAIdtwqLz+c9/PubPnx9XXnllvPjiizOVqezy+XwkiQ2bAKRLttVSNWD2SJIk8vl82cav6+vr2z/RG15++eXYsGFDPP3007Fp06ayBZkJo6OjkcvlYmhoKBobGysdBwAAatrw8HA0NTVFoVCIhoaGko496fHS9fX1kc/n45RTTinpBwMAAJTLlC8M/ehHP1rOHAAAACUz5aIDAAAwW0x6YehrPfHEE/Hoo4/GCy+8EPv27Rv32oc//OGSBgMAADhUUy46X/3qV+MLX/hCLF68OBYtWhR1dXVjr7321wAAAJU25aLzla98Ja644orI5XLlzAMAAHDYprxHZ+/evZHJZMqZBQAAoCSmXHRWr14d/+t//a9yZgEAoEb1D/bHugfWRf9gf6WjkBITLl1bv379uK8feuihePjhh+PXf/3Xo76+ftxr5bzVFACA9Oof7I/c5lwUR4rRva07CmsLkWm2kojDM2HR+fdl5p3vfGdZwwAAUHt6tvdEcaQYERHFkWL07uhVdDhsExada665ZqZyAABQozqWdkT3tu4ojhSjZUFLtC9pr3QkUmBa9+gAAECpZZozUVhbiN4dvdG+pN1sDiUx4WEEH/nIR2JkZGTKg1122WXx3HPPHXYoAABqS6Y5E59+16eVHEpmwhmdJ554Ivr7+2P+/PlTGuyJJ56Il156qSTBAIDqkAwkkW3LVjoGwLRMunTt+uuvn4EYAEC1SnYqOsDsM2HR+da3vjVTOQAAAErGYQQAwAGSgSSSnUlERHRt7Rp7nm3Nmt0BZoWaLDr5fD7a29sjm/V/1ABwMNm28YWmM9tZwTRAGiVJEr29vWUbvyaLzvr166OxsbHSMQAAoGZls9lYuXJlbNq0qSzjT3i8NABAttUKCGD2UXQAgAnZkwPMRpMuXdu2bduUBlq+fPlhhwEAACiFSYvOVVddFXV1dRERsX///oO+p66uLu6///7SJgMAADhEkxadlpaWeOmll2LVqlVx1llnxQknnDATuQAAAA7ZpEXnjjvuiB/84AeRJElcfvnlcdxxx0U2m40zzjjDyWUAAEBVmtJhBCeddFJcccUV8bWvfS0+8IEPxPe+9704//zz47rrrou9e/eWOyMAAMC0TOsenblz58Y73/nOmDNnTgwPD8d3v/vd2Lt3b8ybN69c+QAA4JD1D/ZHz/ae6FjaEZnmTKXjMIOmXHT++Z//Ob75zW/Gt7/97Vi8eHGsWrUqPvvZz8b8+fPLmQ8AAA5J/2B/5DbnojhSjO5t3VFYW1B2asikRefWW2+Ne++9N/bt2xdnnnlmbNy40YEEAABUvZ7tPVEcKUZERHGkGL07ehWdGjJp0bntttvijW98Y/zWb/1W7Nq1K+64446Dvi+fz5c8HAAAHKqOpR3Rva07iiPFaFnQEu1L2isdiRk0adE566yzxu7RAQCA2SLTnInC2kL07uiN9iXtZnNqzKRF59prr52JHAAAUHKZ5oyCU6OmdLw0AADAbDLlU9fWrl170CVsdXV1MXfu3GhpaYlVq1bFu9/97pIGBADSKRlIItuWrXQMIKWmPKPT0dERo6OjsWzZsjj33HPj3HPPjWXLlsWePXvijDPOiEWLFsWf//mfR6FQKGdeACAlkp1JpSMAKTblGZ3HHnssLrnkkjj77LPHPf/7v//7eOihh+Kzn/1sLFmyJO66667I5XIlDwoAADBVUy4627Zti//23/7bAc/f+ta3xsaNGyMi4m1ve1vccsstpUsHAKRKMpCMzeR0be0ae55tzVrGBpTUlIvO0UcfHffcc098/OMfH/f8nnvuiaOPPjoiIl544YVobGwsbUIAIDWybeMLTWe2s4JpgDSbctH51Kc+FZ/5zGfif//v/x1tbW1RV1cX//f//t8YGRmJdevWRUTET37yE8vWAACAipty0VmxYkX83d/9Xdx7771RLBZj//79sXz58jjzzDNj/vz5ERHx3ve+t2xBAYB0ybZaqgaUz5SLTkTE/Pnz49xzzy1XFgCghtiTA5TTtIrOD37wg7jrrrvipz/9aUREnHDCCXHeeefFsmXLyhIOAADgUEz5Hp1vfvOb8alPfSrmzJkTv/3bvx2//du/HXV1dfGpT30q7r333nJmBAAAmJYpz+jcfvvt8fGPf3zc0rU1a9bE3XffHbfddlusWrWqLAEBAACma8ozOrt27YpToIGgqwAAIABJREFUTz31gOennnpq7Nq1q6Shyi2fz0eSuI0ZAAAqJUmSyOfzZRt/ykXn+OOPj/vuu++A5/fee28cf/zxJQ1VbuvXr49s1gZIAAAOTf9gf6x7YF30D/ZXOsqslc1mY/369WUbf8pL1y655JK47rrr4nvf+97Y4QOPP/547Ny5M/70T/+0bAEBAKCa9A/2R25zLoojxeje1h2FtYXINGcqHYt/Z8ozOm9/+9vjtttui0wmE88++2zs2rUrMplM3HbbbXHaaaeVMyMAAFSNnu09URwpRkREcaQYvTt6K5yIg5nW8dLHHntsfPSjHy1XFgAAqHodSzuie1t3FEeK0bKgJdqXtFc6EgcxYdHZtm3blAdavnz5YYcBAIBql2nORGFtIXp39Eb7knbL1qrUhEXnqquumtIgdXV1cf/995ckEAAAVLtMc0bBqXITFp1vfetbM5UDAACgZKZ8GAEAAMBsMeGMzt/+7d9OeaAPfvCDhx0GAACgFCYsOt///vfHfb1z587Yt29ftLS0REREsViM+vr6aG1tVXQAgKqQDCSRbXMxONS6CYvO5z73ubFff/nLX46FCxfG1VdfHfPnz4+IiD179sRf/uVfxtKlS8ubEgBgipKdig4wjT06X/7yl+MjH/nIWMmJiJg/f3586EMfiq985StlCQcAAHAopnxh6C9/+cv48Y9/HCeeeOK4508++WS8/PLLJQ8GADBVyUASyc4kIiK6tnaNPc+2Zs3uQI2actHp6OiIv/iLv4gdO3bEm9/85oiI+NGPfhRbtmyJc845p2wBAQAmk20bX2g6s50VTANUgykXnd///d+PxYsXR09PT9xzzz0REbF48eK4/PLL473vfW/ZAgIAAEzXpEXn+eefj0WLFkVExHvf+16lBgCoatlWS9WAKRSdCy+8ME444YR4xzveEStXrhxbtgYAUI3syQEiplB0tmzZEtu2bYuHH344PvOZz8RLL70Ub3/72+Md73hHnHbaaXHUUUfNRE4AAIApm7TovO51r4uVK1fGypUrIyJiYGAgtm7dGl/96lfjxhtvjJNOOilWrlwZ73jHO2Lx4sVlDwwAADCZKR9G8Kq2trZoa2uLiy66KIaGhuK73/1ubN26NW677ba44IIL4r/+1/9ajpwAAABTNu2i81pNTU1x1llnxVlnnRUvv/xyjIyMlCoXAADUnP7B/ujZ3hMdSzsi05ypdJxZbcpF5/HHH4/HHnssnnnmmdi7d280NTVFa2trnHLKKbFo0aKor6+PhQsXljMrAACkVv9gf+Q256I4Uozubd1RWFtQdg7DhEXn5z//edx1111RKBRiaGgoWltbY9GiRfG6170udu7cGffff3/87Gc/i1NPPTUuvPDCOOmkk2YqNwAApErP9p4ojhQjIqI4UozeHb2KzmGYsOj87u/+bqxYsSKuuuqqyGQyUV9ff8B7nnvuuXjggQdiw4YNsWbNmli9enXZwgIAQFp1LO2I7m3dURwpRsuClmhf0l7pSLNaXV9f3/5f9eKuXbvimGOOmdJA+/fvj3/913+No48+umThSm10dDRyuVwMDQ1FY2NjpeMAAMA4/YP90bujN9qXtNfEbM7w8HA0NTVFoVCIhoaGko494YzOVEtORERdXV1VlxwAAKh2meZMTRScmTDpYQTbtm2b0kDLly8/7DAAAAClMGnRueqqq6Kuri4iXlmedjB1dXVx//33lzYZAADAIZq06LS0tMRLL70Uq1atirPOOitOOOGEmcgFAABwyCYtOnfccUf84Ac/iCRJ4vLLL4/jjjsustlsnHHGGTb0AwAAVWnOVN500kknxRVXXBFf+9rX4gMf+EB873vfi/PPPz+uu+662Lt3b7kzAgBUXDKQVDoCMA2Tzui81ty5c+Od73xnzJkzJ4aHh+O73/1u7N27N+bNm1eufAAAVSHZmUS2LVvpGMAUTbno/PM//3N885vfjG9/+9uxePHiWLVqVXz2s5+N+fPnlzMfAADAtE1adG699da49957Y9++fXHmmWfGxo0bHUgAANSEZCCJZOcrS9a6tnaNPc+2Zs3uQJWbtOjcdttt8cY3vjF+67d+K3bt2hV33HHHQd+Xz+dLHg4AoJKybeMLTWe2s4JpgOmYtOicddZZY/foAAAAzAaTFp1rr712JnIAAFS1bKulajCbTHi89P79+6c12HTfDwAwW9iTA7PLhDM6H/zgB+OCCy6Id73rXXHUUUf9yvcNDAzE3XffHccdd1xccMEFJQ85mdHR0bjmmmuivr4+fv7zn8d5550Xq1atmvEcAABAdZiw6PzxH/9xfPGLX4ybb7453vKWt8SSJUti0aJFMXfu3BgdHY2f/OQn8dhjj8Xo6Gj8zu/8Trz//e+fqdzjHHXUUXHzzTdHfX197NmzJy6++GJFBwAAatiERWfp0qWxYcOGePrpp+PBBx+Mxx57bOyS0MbGxmhtbY1LLrkkVq5cGUccMa27R+PBBx+MLVu2xI4dO2J0dDTuu+++qK+vH/eeO++8M77+9a/Hnj17YsWKFXHllVfGokWLDhhrzpz/vwLvxRdfjF//9V+fVhYAAEiz/sH+6NneEx1LOyLTnKl0nBkxpXZy3HHHxdq1a0v6wb/4xS9i+fLlsWLFivjrv/7rA17/xje+Ebfffnv80R/9UbS0tMTGjRtj3bp1cfPNNx90vOHh4bjuuuviySefjI985CMlzQoAALNV/2B/5DbnojhSjO5t3VFYW6iJsjO9aZgSenVp2SOPPHLQ1+++++5Ys2ZNnH766RERcfXVV8eFF14YAwMD0dbWFvfcc0/09fVFY2Nj3HDDDdHY2Bg333xz7N69Oy699NJ497vfHQsWLJixPw8AAFSjnu09URwpRkREcaQYvTt6FZ3X2rNnT9xxxx3x6KOPxu7duw84Ye3LX/5yyULt3bs3du7cGZdeeunYs5aWljj22GPj8ccfj7a2tli9enWsXr167P3z5s2LiFf268ybNy/mzp1bsjwAADBbdSztiO5t3VEcKUbLgpZoX9Je6UgzYspFZ8OGDfEv//IvcfbZZx90n0wpDQ8Px759++INb3jDuOcLFy6MF1544YD3P/HEE3HLLbdEXV1d/PKXv4wPf/jDceSRR5Y1IwAAzAaZ5kwU1haid0dvtC9pr4nZnIhpFJ1t27ZFV1dXLFmypJx5ImL69/EsXbr0V+7dOZh8Pj82A5TNZiObdS4+AADplWnOVE3BSZIkkiSJiFdWZpXLlIvO0UcfPWMXgjY1NcWcOXNi9+7d456/8MILsXDhwsMef/369dHY2HjY4wAAANPz2omG4eHh2LRpU1k+Z87kb3nFJz/5yfirv/qr2L59e+zduzf27ds37n+lNG/evGhtbR13UMHg4GA888wzsWzZspJ+FgAAkD5TntG56qqrIiLisssuO+jr999//7Q+eHh4OJ599tl4+umnIyJiYGAg6uvr47jjjoujjjoqzjnnnNi4cWMsWbIkmpub45ZbbomTTz452trapvU5AABA7Zly0ens7CzpBz/00EOxYcOGsa8/9rGPRUREV1dXnHLKKfG+970vdu/eHTfddNPYhaGvli0AAICJ1PX19c3MxpsqMDo6GrlcLoaGhuzRAQCAChseHo6mpqYoFArR0NBQ0rEnnNF59NFHpzzQW9/61sMOAwAAUAoTFp0//MM/HPd1XV1dRPz/459f/Tpi+nt0Kimfz0d7e7tjpQGAqpAMJJFt8+8SakuSJNHb21u28Sdcuvbyyy+P/fqf/umf4tZbb40PfehD8Zu/+ZsREfHDH/4wvvSlL8Xv/d7vxcqVK8sWslQsXQMAqtEVyRXRmS3tfmiYDSq2dK2+vn7s15s2bYprr7123PHOp512WsyfPz/Wr18/K4oOAABQG6Z86tqzzz47bqnaq+rq6uJf//VfSxoKACDtkoEkkp2v3A7ftbVr7Hm2NWsZG5TAlIvOaaedFhs2bIjLL7883vzmN0dExI9+9KPYuHFjnHbaaWULCACQRtm28YXG0jUorSkXnWuuuSY2btwYf/RHfxT79u2LiIg5c+bEe97znrj88svLFhAAAGC6plx0Ghoa4pprrolPfOITUSwWIyKipaUlXv/615ctHABALci2WqoGpTblovOq17/+9dHW1laOLAAANcmeHCi9OZUOUAn5fD6SJKl0DAAAqFlJkkQ+ny/b+BPeo5M27tEBAIDS6R/sj57tPdGxtCMyzZlp//5y3qNTkzM6AADA4ekf7I/c5lxc/+3rI7c5F/2D/ZWONI6iAwAATFvP9p4ojrxySFlxpBi9O3ornGg8RQcAAJi2jqUd0bKgJSIiWha0RPuS9gonGm/ap64BAABkmjNRWFuI3h290b6k/ZD26JSTogMAABySTHOm6grOqyxdAwAAUkfRAQAAUqcmi44LQwEAoLJcGFpCLgwFAIDq4cJQAACAaVB0AACA1FF0AACA1FF0AACA1FF0AABSKhlwyiy1S9EBAEipZKeiQ+1SdAAAgNQ5otIBAAAonWQgGZvJ6draNfY825qNbFu2UrFgxtVk0cnn89He3h7ZrL/sAEC6ZNvGF5rObGcF08CvliRJ9Pb2lm38miw669evj8bGxkrHAACAmpXNZmPlypWxadOmsoxvjw4AQEplW61eoXYpOgAAKWVPDrVM0QEAAFJH0QEAAFJH0QEAAFJH0QEAAFJH0QEAAFJH0QEAAFJH0QEAAFJH0QEAAFKnJotOPp+PJEkqHQMAAGpWkiSRz+fLNn5dX1/f/rKNXmVGR0cjl8vF0NBQNDY2VjoOAADUtOHh4WhqaopCoRANDQ0lHbsmZ3QAAIB0U3QAAIDUUXQAAIDUUXQAAIDUUXQAAIDUUXQAAJiWZMA1HVQ/RQcAgGlJdio6VD9FBwAASJ0jKh0AAIDqlwwkYzM5XVu7xp5nW7ORbctWKhb8SooOAACTyraNLzSd2c4KpoHJWboGAACkjqIDAMC0ZFstVaP6KToAAEyLPTnMBjVZdPL5fCSJYxEBAKBSkiSJfD5ftvHr+vr69pdt9CozOjoauVwuhoaGorGxsdJxAACgpg0PD0dTU1MUCoVoaGgo6dg1OaMDAACkm6IDAACkjqIDAACkjqIDAACkjqIDAACkjqIDAACkjqIDAACkjqIDAACkjqIDAACkjqIDAACkjqIDAACkjqIDAACkjqIDAACkjqIDAACkjqIDAEBVSAaSSkcgRRQdAACqQrJT0aF0FB0AACB1jqh0AAAAalcykIzN5HRt7Rp7nm3NRrYtW6lYpEBNFp18Ph/t7e2RzfrLAwBQSdm28YWmM9tZwTTMpCRJore3t2zj12TRWb9+fTQ2NlY6BgAA1KxsNhsrV66MTZs2lWV8e3QAAKgK2VarbSgdRQcAgKpgTw6lpOgAAACpo+gAAACpo+gAAACpo+gAAACpo+gAAACpo+gAAACpo+gAAACpo+gAAACpo+gAAACpo+gAAACpo+gAAACpo+gAAACpo+gAAACpo+gAAACpo+gAAACpo+gAAACpo+gAAACpo+gAAJBqyUBS6QhUgKIDAECqJTsVnVqk6AAAAKlzRKUDAABAqSUDydhMTtfWrrHn2dZsZNuylYrFDFJ0AABInWzb+ELTme2sYBoqwdI1AAAgdWqy6OTz+UgSm9IAAGpBttVStWqUJEnk8/myjV/X19e3v2yjV5nR0dHI5XIxNDQUjY2NlY4DAAA1bXh4OJqamqJQKERDQ0NJx67JGR0AACDdFB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AADgEyUBS6QhMQNEBAIBDkOxUdKqZogMAAKTOEZUOAAAAs0UykIzN5HRt7Rp7nm3NRrYtW6lYHISiAwAAU5RtG19oOrOdFUzDRCxdAwAAUkfRAQCAQ5BttVStmik6AABwCOzJqW6KDgAAkDqKDgAAkDqKDgAAkDqKDgAAkDqKDgAAkDqKDgAAkDqKDgAAkDqKDgAAkDqKDgAAkDqKDgAAkDqKDgAAkDqpKjpDQ0PR3t4ef//3f1/pKAAAQAWlquj87d/+bZx88smVjgEAAFTYEZX64AcffDC2bNkSO3bsiNHR0bjvvvuivr5+3HvuvPPO+PrXvx579uyJFStWxJVXXhmLFi066HhPPvlkvPjii/GmN71pJuIDAABVrGIzOr/4xS9i+fLlsXbt2oO+/o1vfCNuv/32+OQnPxkbN26M0dHRWLdu3a8c72/+5m/i4osvLlNaAABgNqnYjM6qVasiIuKRRx456Ot33313rFmzJk4//fSIiLj66qvjwgsvjIGBgWhra4t77rkn+vr6orGxMXK5XCxevDiOOeaYGcsPAABUr4oVnYns3bs3du7cGZdeeunYs5aWljj22GPj8ccfj7a2tli9enWsXr06IiLuuOOO+OEPfxhXX311PP300zFv3rxobm6O5cuXV+qPAAAAVFBVFp3h4eHYt29fvOENbxj3fOHChfHCCy8c8P6LLrooLrroooiIuPXWW+Poo49WcgAAoIZVZdHZv3//If9e+3QAAICqLDpNTU0xZ86c2L1797jnL7zwQixcuPCwx8/n8zFv3ryIiMhms5HNZg97TAAAKIVkIIlsW3r/fZokSSRJEhGvbFkpl6osOvPmzYvW1tZ45JFHYsWKFRERMTg4GM8880wsW7bssMdfv359NDY2HvY4AABQasnOdBed1040DA8Px6ZNm8ryORUrOsPDw/Hss8/G008/HRERAwMDUV9fH8cdd1wcddRRcc4558TGjRtjyZIl0dzcHLfcckucfPLJ0dbWVqnIAADALFGxovPQQw/Fhg0bxr7+2Mc+FhERXV1dccopp8T73ve+2L17d9x0001jF4ZeddVVlYoLAABlkwwkkex8ZTlX19ausefZ1myqZ3fKqa6vr+/Qd/7PMqOjo5HL5WJoaMjSNQAAqtIVyRXRme2sdIwZMTw8HE1NTVEoFKKhoaGkY88p6WgAAABVQNEBAIAqkm21VK0UarLo5PP5sSPtAACgmtTKnpwkSSKfz5dtfHt0AACAirBHBwAAYBoUHQAAIHUUHQAAIHUUHQAAIHUUHQAAIHUUHQAAIHVqsui4RwcAACrLPTol5B4dAACoHu7RAQAAmAZFBwAASB1FBwAASB1FBwAASB1FBwAASB1FBwAASB1FBwAASJ2aLDouDAUAgMpyYWgJuTAUAACqhwtDAQAApkHRAQCAGpAM1NbWDUUHAABqQLJT0QEAAJjVjqh0AAAAoDySgWRsJqdra9fY82xrNrJt2UrFmhGKDgAApFS2bXyh6cx2VjDNzLJ0DQAASB1FBwAAakC2Nd1L1f69miw6+Xw+kqS2Tp0AAKC2VduenCRJIp/Pl238ur6+vv1lG73KjI6ORi6Xi6GhoWhsbKx0HAAAqGnDw8PR1NQUhUIhGhoaSjp2Tc7oAAAA6aboAAAAqaPoAAAAqaPoAAAAqaPoAAAAqaPoAAAAqaPoAAAAqaPoAAAAqaPoAAAAqaPoAAAAqVOTRSefz0eSJJWOAQAANStJksjn82Ubv66vr29/2UavMqOjo5HL5WJoaCgaGxsrHQcAAGra8PBwNDU1RaFQiIaGhpKOXZMzOgAAQLopOgAAQOooOgAAQOooOgAAQOooOgAAQOooOgAAQOooOgAAQOooOgAAQOooOgAAwCFLBpJKRzgoRQcAADhkyU5FBwAAYEYcUekAAADA7JIMJGMzOV1bu8aeZ1uzkW3LVirWOIoOAAAwLdm28YWmM9tZwTQHV5NL1/L5fCRJda4lBACAWpAkSeTz+bKNX9fX17e/bKNXmdHR0cjlcjE0NBSNjY2VjgMAALNeMpAc8nK14eHhaGpqikKhEA0NDSXNVZMzOgAAQGlUy56cf0/RAQAAUkfRAQAAUkfRAQAAUkfRAQAAUkfRAQAAUkfRAQAAUkfRAQAAUkfRAQAAUkfRAQAAUkfRAQAAUkfRAQAAUkfRAQAAUkfRAQAAUkfRAQAAUkfRAQAAUkfRAQAAUqcmi04+n48kSSodAwAAalaSJJHP58s2fl1fX9/+so1eZUZHRyOXy8XQ0FA0NjZWOg4AANS04eHhaGpqikKhEA0NDSUduyZndAAAgHRTdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNRRdAAAgNSpyaKTz+cjSZJKxwAAgJqVJEnk8/myjV/X19e3v2yjV5nR0dHI5XIxNDQUjY2NlY4DAAA1bXh4OJqamqJQKERDQ0NJx67JGR0AACDdFB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1FB0AACB1jqh0gFI566yzYtmyZRER8ba3vS0uuuiiCicCAAAqJTVFp6mpKW666aZKxwAAAKpAxYrOgw8+GFu2bIkdO3bE6Oho3HfffVFfXz/uPXfeeWd8/etfjz179sSKFSviyiuvjEWLFh10vJGRkfiDP/iDOPLII+OSSy6J3/iN35iJPwYAAFCFKrZH5xe/+EUsX7481q5de9DXv/GNb8Ttt98en/zkJ2Pjxo0xOjoa69at+5Xj/d3f/V3cdNNN8eEPfzhuuOGGcsWGWSFJkkpHgLLz3zm1wH/ncOgqVnRWrVoVF110UZx00kkHff3uu++ONWvWxOmnnx5tbW1x9dVXx//5P/8nBv5fe/cf1VT9/wH8uU2HExygkghaGohKhoKZ4klrpE44hp4jp6PHk0cpozKVlJCOpoiHTulJ5chBSz3H8kcnTx6xH2JiR80s86RN4wzI2dKaTFTGxjjolPH5g6/3640fDraxsT0ff7H3+33f93XpFfLife/76nQAgMOHDyMzMxNr1qwBAISEhAAAYmJi0KtXL9TW1nbNhRB5If7DSP6AeU7+gHlO1HleueuazWbDlStXEB8fL7RFREQgPDwcWq0WADBz5kxs2bIFeXl5aGhoQGNjIwCguroaFosFSqXSI7F7mjf8QHRnDK6a25l5OnpsR8Y7MtYb/ht7mjd8D5jnnR/PPHeMN3wPmOedH888d4w3fA+Y550f7+157pWbEVgsFtjtdoSGhoraQ0JCWl2p+fvvv7Fp0yb07t0b9+/fx4oVKyCVtqzhmpqahPl91TfffIPExESfjcFVczszT0eP7ch4R8Y6MsZmszHPu3EMzHPmOcA874p5mOeexzx3/zzenucP8vvB7+muJDlx4oTrZ+0AjUaDd955R7QZwc2bN/Hyyy9j165dok0F3nzzTSQmJmL+/PmdOteDeYmIiIiIyHscOHAAYWFhLp3TK1d0goODIZVKYTKZRO21tbXCszid0a9fPxw4cAAKhQISicTZMImIiIiIyAlNTU1oaGhAv379XD63VxY6crkcUVFR0Gg0GDt2LACgqqoKRqNReCloZ0ilUpdXikRERERE1HlBQUFumddjhY7FYkF1dTUMBgMAQKfTQSaTITIyEgqFArNmzUJhYSFiYmIwcOBAFBUVIS4uDtHR0Z4KmYiIiIiIugmPPaNz9OhRfPTRRy3aN2/ejDFjxgAA9u3bJ3phaFZWVpsvDCUiIiIiInrA45sREBERERERuZpXPqPjKZ9//jl+/fVXyGQypKenCytLRN2Z3W5HZmYmrl69itTUVLz66qtCX2lpKYqLiwEA6enpwjNxRN1NW3neXv4TdTdt5XNlZSW2bt0KmUyGxsZGvPXWW04900zkSW3leX19PVauXAmZTIY7d+4gLS0NU6dObXcuruj8H51Oh6KiInz88cewWq1Yvnw5tm/fLmx5TdSdVVdX48KFCzAYDMIPDKvViiVLlmD79u24c+cOMjMzsXPnTuY8dVut5Xl77UTdUWv5fPv2bfTq1QuBgYHQ6/XYsGEDtm3b5uFIiTqvtTy32+1oamqCTCaD1WrFggUL8NVXX7U7j8+s6Pz4448oLi7Gn3/+ifr6etF7eR7Yv3+/6JmfFStWCM/8/PPPPxg2bBgkEgn69OmDnj17wmAw4PHHH/fE5RAJnM1tAHjsscdazFteXo5Ro0YhICAAAQEBGDBgAP7991888cQTbr8mov9yV563107U1dyV5w9vyyuXy1t9aTpRV3FXnj+c1w0NDRg6dOgjY/GZ/xPu3r2LhIQEzJ07t9X+kpIS7NmzB0uXLkVhYSHq6+uxbt06oX/o0KG4ePEibDYbjEYj9Ho9zGZzV4VP1CZnc7stZrMZffr0ET4HBQUx58lj3JXnRN7E3Xl+7949bNq0qdMvVidyBXfmucViwbJly/Daa69h0qRJjxzvMys6D+7R02g0rfYfOnQIs2fPxuTJkwEA2dnZmDdvHnQ6HaKjozFkyBCkpKTg3XffRf/+/TFs2DC3vLiIqKOcze22KJVK1NXVCZ+tViuUSqULIydynLvynMibuDPPGxsbsX79eqhUKowfP961gRN1gDvzXKlUoqCgACaTCRkZGVCpVKI/2v6Xz6zotMdms+HKlSuIj48NQDRBAAAJ/ElEQVQX2iIiIhAeHg6tViu0paamoqCgABkZGQgICEBERIQnwiVymKO53ZrY2FiUlZXBZrPBYrHgxo0bGDx4sLtDJuowZ/KcqLtwJs/tdjs++OADjBo1CjNmzHB3qESd5kye22w24WuFQgG5XI6ePXu2e4zPrOi0x2KxwG63IzQ0VNQeEhKC2tpa4fPKlStx7949KBQKLF++vKvDJOowR3M7Pz8fly9fxt27d6HVarFx40YEBQVhzpw5yMzMhEQiweLFi7kRAXklZ/JcKpW22U7kTZzJ85MnT+LMmTO4ffs2fv75ZwQGBiI/P7+rL4HokZzJc71ej6KiIkgkEty/fx/p6eno1atXu+fzi0KnqcmxjeVae4EpkTdzNLdXrVrVartarYZarXZlSEQu52yet9VO5E2cyfOkpCQkJSW5OiQil3Mmz4cPH46CgoIOnc8v/qQVHBwMqVQKk8kkaq+trUVISIiHoiJyHnOb/AHznPwB85z8QVfnuV8UOnK5HFFRUaKHoqqqqmA0GvlCLerWmNvkD5jn5A+Y5+QPujrPfebWNYvFgurqahgMBgDNLwCVyWSIjIyEQqHArFmzUFhYiJiYGAwcOBBFRUWIi4vjbj3k9Zjb5A+Y5+QPmOfkD7wpzyUnTpxw7GY5L3f06NFWn7HZvHkzxowZAwDYt2+f6OVEWVlZopcTEXkj5jb5A+Y5+QPmOfkDb8pznyl0iIiIiIiIHvCLZ3SIiIiIiMi/sNAhIiIiIiKfw0KHiIiIiIh8DgsdIiIiIiLyOSx0iIiIiIjI57DQISIiIiIin8NCh4iIiIiIfA4LHSIiIiIi8jksdIiIyOvt3r0bKpUKKpUKp06d8nQ4DtFoNFCpVGhsbGxzzOrVq4XrMhgMXRgdEZHvY6FDROQHMjMzoVKpcOTIEVH73bt3MWPGDNEv2jqdDrm5uUhLS0NycjIWLVrUanGxf/9+pKWlYfr06Vi1ahVqampE/VqtFhkZGZg2bRrmz5+PX375xalrGDFiBA4ePIjExEQAgNFohEqlQnJyMhoaGkRjS0tLoVKpsGTJEqfO2RnLli3DDz/84NDYnJwcbNu2zc0RERH5JxY6RER+IiwsDKWlpaK206dPo3fv3qK2y5cvIywsDGvWrMGuXbswffp05OXlQaPRCGNKSkqwZ88eLF26FIWFhaivr8e6deuEfrPZjJycHMTGxuLTTz+FWq3G2rVrce3atU7H36NHD/Tt2xdyuVzU3qdPnxaF2LFjxxAWFtbpc3VWXV0dysvLMW7cOIfGBwUFISQkxM1RERH5JxY6RER+YtKkSaioqMCNGzeEtu+//x5Tp04VjUtOTsbixYsRFxeHiIgIzJ49GwkJCThz5oww5tChQ5g9ezYmT56M6OhoZGdn49KlS9DpdACA48ePIyAgAEuXLsWQIUMwb948jBw5El9//bUwh8lkQn5+PlJTU5GcnIy33367U7dvTZkyBceOHRM+37x5E2VlZXj++edF43777TcsWrQIarUaM2fORE5OjtA3Z84cfPfdd6LxKpUK58+fB/D/t6GdP38eCxYsQEpKClavXo26ujrRMefOncOIESOgVCqFNo1G0+4xRETkHix0iIj8RO/evTFx4kRhVefWrVsoKyuDSqV65LFms1n45d1ms+HKlSuIj48X+iMiIhAeHg6tVgsAqKioQHx8PCQSiTAmISEB5eXlwuc1a9bg+vXryM/Px44dO/DSSy+1+zxLW1544QWUl5fj5s2bAJqLrAkTJiAwMFAY09jYiNzcXKjVanz22WfYtGkTnnnmmQ6fa+/evVi5ciU2b94MvV6PPXv2iPrPnj0r3Frn6DFEROQeLHSIiPyIWq0WCp3S0lKMHz8eQUFB7R5z6tQpXLt2DVOmTAEAWCwW2O12hIaGisaFhISgtrYWAFBbW9vilqzg4GCh//fff0dlZSVyc3Px9NNPY9CgQVCr1Xj88cc7fE2BgYFITEzE8ePHATSvUk2bNk00xmq1or6+HpMnT0Z4eDiioqKQlpbW4XO9/vrrGDlyJIYPH46UlBRcvHhR6GtsbMS5c+cwYcIEh48hIiL3YaFDRORHxo4dC6vVioqKCpSWlrYoCP6rrKwMGzZsQFZWFgYOHAgAaGpqeuR5HjVGr9cjMjLSZc/RTJs2DceOHUNlZSUsFgueffZZUX9wcDCSkpKQnp6OvLw8lJSUtNjAwBFPPvmk8HXfvn2Fwg1o3nxBoVBg6NChDh9DRETuw0KHiMiPyGQyvPjii9i2bRtMJlOLguBhFRUVyMnJQUZGhrCaAzQXDVKpFCaTSTT+4VWc0NDQFr/Qm81mod+RYqkjxo0bB7PZjKKiIiQlJUEmk7UY8/7772PDhg0YPHgwvvzySyxcuBBmsxkAIJVKRTHdv3+/1fP06NFD+FoikcButwufz54922I151HHEBGR+7DQISLyM2q1GpcuXYJKpRL9Ev6wy5cvIzs7G6+88gpSU1NFfXK5HFFRUaJd2KqqqmA0GhEbGwugeStojUYjKh4uXLiAkSNHAmhe5TAYDLh165ZLrulBAXfp0qV2V6liY2OxcOFC7NixA1arFRcuXADQXLw9vD32X3/91eEYWns+h4iIPIeFDhGRn4mKikJxcTHeeOONVvv1ej2ysrKQlJSEqVOnoqamBjU1NbBarcKYWbNm4eDBgzh9+jR0Oh02btyIuLg4REdHA2jeCe3OnTvYunUrrl69ii+++ALl5eVC0RQfH4/hw4dj7dq1+OOPP2AwGFBaWurU9tOLFi1CcXExYmJiWvRVVVVh586d0Gq1MBqNOHnyJBoaGjBo0CAAwOjRo1FSUoKKigpUVlbik08+6dC5b9y4gevXr4s2aCAiIs9q/U95RETk04KDg9vsO3XqFCwWCw4fPozDhw8L7Wq1WtiSOSUlBSaTCVu2bIHVasXYsWORlZUlmv/DDz9EQUEBvv32W4SHhyM3N1e02UBeXh4KCwvx3nvvwW63Izo6Gk899VSnr0kul7d4x84DAQEB0Ov1OHLkCKxWKyIiIpCdnY1hw4YBAObNmweDwYDly5ejf//+yMzMFFZ7HHH27FkkJCS0eX4iIup6khMnTrj2RmkiIiIX2717N86fP4+tW7d6OpRW5eTkYOLEiS1u83OE0WjE3LlzsXfvXkRGRrohOiIi/8Rb14iIqFvQarVITk7GTz/95OlQWhg9ejSee+65Dh+3fv16LFy40A0RERERV3SIiMjrWSwW1NXVAWjeolmhUHg4IteoqakRtrkeMGBAm5tDEBFRx/EnKhEReT2lUgmlUunpMFyub9++ng6BiMhn8dY1IiIiIiLyOSx0iIiIiIjI57DQISIiIiIin8NCh4iIiIiIfA4LHSIiIiIi8jksdIiIiIiIyOew0CEiIiIiIp/zP5SdZWFMeQjwAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 960x960 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "mfbins, ymm_hdf5_distinct, ymm_hdf5_subhalos = np.loadtxt('mfunc_uchuu_results_hdf5_all.out')\n",
    "\n",
    "norm = 1/(np.log10(mfbins[1:len(mfbins)]) - [np.log10(mfbins[0:(len(mfbins)-1)])])\n",
    "\n",
    "plt.style.use('classic')\n",
    "plt.figure(figsize=(12, 12))\n",
    "plt.plot(mfbins[0:34], norm[0,0]*ymm_hdf5_distinct[0:34]/(2000**3), 'g.', label = \"Big Uchuu (hdf5) Distinct\")\n",
    "plt.plot(mfbins[0:34], norm[0,0]*ymm_hdf5_subhalos[0:34]/(2000**3), 'g+', label = \"Big Uchuu (hdf5) Subhalos\")\n",
    "plt.title('Big Uchuu (ascii) VS (hdf5) [z = 0.86]: mass function', fontsize = 12)\n",
    "# plt.vlines(bins[:-1], 10, 40, linestyles = 'dashed', lw=.5)\n",
    "plt.legend(loc = 1, borderaxespad = 0.)\n",
    "plt.xlabel('M200c [Msun/h]')\n",
    "plt.ylabel('dn(M)/dlogM [h^3/Mpc^3]')\n",
    "plt.xscale('log')\n",
    "plt.yscale('log')\n",
    "# plt.show()\n",
    "plt.savefig(\"mfunc_ascii_vs_hdf5_biguchuu_z086.png\", bbox_inches='tight')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAysAAANJCAYAAAAMa6MNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAMTQAADE0B0s6tTgAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3XlcjXn/P/DXqZw6ooVKUiQKWUNuuymUqBjGKMswlmGyRWPmnsZS1jCyjGxjyTDGnhZLNSbM2Mdk+zIoW1otLWSJOr8/zq/rdjqnU1HOSa/n4+FR51rf5zrn5Hqd6/P5XKK4uDgpiIiIiIiINIyWugsgIiIiIiJShmGFiIiIiIg0EsMKERERERFpJIYVIiIiIiLSSAwrRERERESkkRhWiIiIiIhIIzGsEBERERGRRmJYIarigoKCsGDBArXWkJaWBicnJyQnJ3/wfRcUFGDMmDG4fPlyhe/Ly8sLBw8eBKD4nGNjY/Hdd99VeA2qrFy5Etu2bRMeOzk54cKFC2XaxuTJkxEaGio8vn79OsaMGYNevXohKCioxPULCgowcuRIXLt2TeVyQUFBcHJygpOTE27cuFGmGiurZ8+eCc958uTJH2y/5f35fPtzUF6ys7MxY8YM9OnTB15eXuW67bIIDQ39oK8NUVWgo+4CiKji+Pr64tKlSwAALS0t1KpVC927d8f48eMhFosBAJMmTXqvfVy8eBHTpk3D77//Dm1tbWF6WloavL29sX37dtSrV++99lGRjh07Bn19fbRq1arC97Vu3TpIJBIAgKmpKfbt2wdDQ0MAgLOzMzZv3owrV66gZcuWCuvm5eXhs88+w+jRozFgwACF+RMmTECzZs0wdepUJCUl4eeff8bVq1eRm5sLExMTNG/eHDNmzEC1atWU1paeno6jR49ix44d5fiMgZ9//hmNGjVCUFAQJBKJ8L4oKjIyEjVq1ICWlhaGDBmCTZs2YdmyZSq33b17d0ydOlU4hh/a6dOnsX79eqSkpMDa2hq+vr6wt7cvdvmHDx9i9erVuHjxIvLz89G0aVP4+PjAxsZGWCYzMxPr1q3DmTNn8OrVK1hbW2PBggWoXbs29PX1sW/fPuzatavEMKfJ3v4cJCcnY/jw4fjtt99gbm7+zts8cOAAHj16hE2bNkFfX7+8SlVp8uTJaNeuHUaNGiVMGzJkCAYOHPhB9k9UVTCsEH3kPvvsM3h7eyM/Px/379/H4sWLUb16dYwZMwYAUKNGDTVXqF7h4eFwcXH5IPsyMjISftfW1katWrXkHjs7OyM8PFxpWBGLxXByckJMTIxCWLl//z5u3LgBX19fvH79Gt9++y2aNGmCRYsWwcDAAMnJyTh+/DgKCgqKre3gwYPo0KFDub8fUlJS0Lt3b5iamgKQXR0AgJCQELmT07dPMHv06IFVq1YhKSkJVlZWxW5bLBbLHcMP6f79+5gzZw5GjBiBbt26ISIiAv/973+xbdu2YsPTokWLIJVKsWzZMojFYoSGhuL777/Hzp07IRKJkJeXh+nTp6N+/foICgqCoaEh7t69KwRMkUiEWrVqCSf6ldXbn4PykpqaCjs7O7V/MSKRSCr960OkadgMjOgjp6enh1q1asHU1BTt2rVD9+7dcevWLWF+0WZgaWlp8PX1hYuLC8aOHYu4uDg4OTkhLS3tvWv5448/MHr0aLi4uGDw4MFyTY4A4M6dO/j666/h5uYGX19fpKenC/N8fX2xadMmueXfbk5y8eJFODk5IT8/X5h/5MgRDB48uNh6njx5gsuXL+M///mP3LQ5c+Zg4MCB6Nu3L6ZOnYqEhARhfl5eHn788Ud8+umncHV1xRdffIG//vpLmH/jxg1Mnz4dffr0gaenJ2bOnKm0XmVNazp27IiTJ0/izZs3Sut1dXXF9evX8eDBA7np0dHRaNCgAZo2bYp79+4hLS0N06ZNQ5MmTVC3bl20b98efn5+0NXVLfZYHD9+XO44FHr48CH8/PzQp08ffPXVV0hMTBTmSaVSbNy4ER4eHhgwYAB27twpt66TkxPS09OxZMkSODk54ciRI8I8Q0ND1KpVS/gnEomEefr6+mjRogWOHz9ebL3K+Pr6Cs2k3v5XESIjI2FnZ4cRI0bA2toakydPRvXq1REbG1vsOtevX8egQYPQuHFj1K9fH8OHD0dGRgYyMzMBAFFRUXj9+jVmz56NZs2awcLCAp07d4aBgcE715mfn48BAwYgLi5ObvquXbswevRo4fGFCxcwfvx4uLq6YsSIEQgLC1O53ejoaAwbNgwuLi4YPXo0zp49Kze/tJ+D4cOHAwC8vb3h5OSE0NBQBAcHIyAgQG57ycnJcHZ2RmpqqkItvr6+iI6ORkxMDJycnBAUFFSqvweFf/s2bdoET09PfPbZZ9i9e7fCfn/44Qf069cP7u7u8PPzw9OnTxEUFISrV69i69atcHJyEpqeFW0G9uLFCyxduhQeHh5wc3PD7Nmz8eTJkzLVQFTV8coKURWSkZGBCxcuoEePHsUus2DBAmhra2PNmjXIzs7GypUry2Xf58+fx6JFizBu3Dh07twZ2dnZuHv3rtwyoaGhmDhxIoyMjLBkyRKEhIRg7ty55bJ/Za5evQpjY2PhW39AFkZatWqFkSNHolq1atizZw9++OEHbNu2DWKxGPv27cPNmzcRFBQEAwMD3L9/X/jmOysrC35+fujRo4fQvK7oSZwqdnZ2yMvLQ0JCApo2baow397eHlZWVoiJiRFONKVSKY4ePYr+/fsDkIUALS0t/Pnnn+jXr59cCChOdnY27t+/Dzs7O4V5v/zyC3x8fODr64u1a9di8eLF2LBhAwDZCev+/fvx7bffokGDBtiyZQsSExPRrl07AMC+ffvw1VdfwcvLC87OztDX1xdOzKdPn443b96gYcOGGD16tELzqSZNmuDq1aulPnYAMHfuXCHoSaVSLFiwoNjgBwDfffedyr5KvXv3xvTp05XO+/fff9G2bVvhsUgkgoODA65fv17s9uzt7REXF4d27dpBR0cHMTExaNSoEYyNjQEAZ86cQbNmzbBkyRKcO3cOpqamGDZsmMrPKyA74U1LS8OKFSsU5mlra6Nbt27Clw6Fjh07Jjy+f/8+Zs+eDR8fHzg4OODevXtYsmQJDA0N4ezsrLDNq1evYsmSJfDx8YGjoyN+//13zJo1C7/88gvMzc3L9DkICQnBxIkTsXbtWpiZmUEikeDOnTuYNm0anj17Jlzpi42NRcuWLVG3bl2FbcydOxfLli2DlpYWJk+eDLFYLPcFgyqnTp2Cp6cnQkJCcPHiRQQHB6Ndu3Zo1KgR8vLyhPd2cHAwJBIJ/vnnHxQUFGDSpEm4d+8eWrVqhSFDhkBLS/l3vyEhIbh06RLmz58PiUSCFStWYNGiRVi6dGmpaiAihhWij97OnTuxd+9eFBQUIC8vD+3btxe+zSzqzp07uHr1qlw/k8GDByM4OPi969i+fTvc3d3x+eefAwAsLS3RvHlzuWWGDRsGBwcHALLma6tWrXrv/aqSnp6u0IzI3NwcgwYNEh5PnToV7u7uuH79Olq3bo2HDx+icePGaNKkCQDInTyFhYXBzMwM33zzjRAS3u6PUBJdXV3o6+sjPT1daVgBABcXFxw8eBBffvklRCIRLl68iIcPH6JXr14AZH1hxo8fj1WrVmHDhg2wt7dHu3bt0KdPH9SsWVPpNjMyMgAAtWvXVpjn4eGBrl27ApC9PpMmTcKLFy8gkUgQHh6OAQMG4JNPPgEAzJgxQ+6b68IrJjVq1BCOs0QiwaRJk9CiRQvk5+fj0KFD8PX1xYYNG2BtbS237qlTp0p97ADIXYHYvn07kpKSsH79+mKX/+abb5CXl1fs/OrVqxc7LysrS6E5k6GhIW7evFnsOrNnz8bs2bOFEGlhYYEff/xReK+kpaXh4sWLcHd3x5IlS3Dp0iXMnTsXK1asUNo0sFDt2rVVNvFzcnLC999/j+fPn6N69epIS0vDv//+ix9++AEA8Ntvv8Hd3R39+vUDAFhYWGDQoEE4ePCg0rCyf/9+dOvWTficjB49Gn///TcOHDiACRMmlOlzUNhkzsjISHiP2Nvbw9zcHMeOHYO7uzsAWVgZOnSo0m0YGBhALBYrNK0sjcLPCwBYWVlhz549uHz5Mho1aoSjR4/i+fPnmD17NvT09AAA9evXF9bV0dGBRCIpdp/Pnz/HkSNHMH/+fLRu3RoA8N///hcjR47EnTt30LBhwxJrICKGFaKPXr9+/TB48GAUFBQgNTUVa9aswerVq+Hr66uwbHJyMvT19eXafReelL+vu3fvltjx9O0Tmlq1aiEnJwf5+flyHffLU15enkKH8zdv3mDLli3466+/8PjxY+Tn5+PVq1fCCX3v3r0xY8YMJCQkwNHREd27dxeO0Z07d9C6detSXc0ojq6uLl69elXsfBcXF2zZsgWXL19G69atERMTg7Zt28pdHfr888/Rp08fXLhwAf/3f/+H3bt3Y/fu3Vi7di1MTEyUHgcASjvfF31NAFkncIlEgqSkJLnO8jVr1iyxz4ChoaFcGLS3t0dSUhLCw8MxdepUueOgKkiocu7cOWzfvh3Lly9XefL69jErK6lUWuZ1Nm3aBG1tbaxcuRJisRi7d+/GzJkzsWbNGlSrVg1SqRSmpqaYPHkyRCIRbG1tcenSJRw+fFhlWBk3bpzK/bZp0wY1atTAyZMn0bt3b8TFxcHW1haWlpYAZO/b27dvIyIiQlgnPz9f6XsFAJKSktC7d2+5ac2bN0dSUpKwvff9HLi6uiI2Nhbu7u64cuUKHj16VOIVpndRGBgKGRsbC1f/7ty5g6ZNmwpBpaxSUlKQn58vd9Wwfv36qFGjBpKSkoR9q6qBiBhWiD56b59AWllZ4dmzZ1iwYAG+/vprhT4M73ICVvjtc25urtw324UdqVV9O12Ujs7//iQVPdHR0tJSqO/t9uiFy7+9jKomQIDsG9nc3Fy5ab/99htiYmIwadIk1K9fH2KxGF9//bWwr2bNmmHHjh04c+YMzp07h8mTJ2PMmDEYMmRIqZ+nKk+fPlU5upWZmRnatGmD2NhYNG3aFCdOnFAaPA0MDIQ+G6NGjcKIESMQEREh10/h7WUB2WtW9D2h7DV5+xi/zwlpITs7O+FEt9DTp0/fqa9GcnIy5s+fjylTpqBZs2Yql32fZmDGxsbIysqSm5adnS006VJWV2RkJHbu3Ik6deoAAL7//nt4eHjg/Pnz6Ny5M4yNjaGnpyd3TK2srOT6mL0LLS0t9OjRA8eOHRPCyttXTF68eIHBgwejb9++cusV9yXBu/ydKCsXFxds3rwZqampiImJQdeuXcs0+ENp/x68/f4uXK9wnQ/xPEuqgYgYVoiqHG1tbRQUFOD169cKJ6aWlpbIzc1FcnKyEHBUNWsBZE1GRCIREhIS5Nrw37p1C/r6+sLJm7W1NS5evPjO344aGhrKdUzNysqSe1zYJOfJkycwMzMDANy+fVvlNhs1aoTU1FTk5eUJQzlfu3YNPXr0EOrMyMjA06dP5dYzMDCAi4sLXFxc8Ntvv+Hw4cMYMmQIGjZsiL/++gtSqfSdTuJTU1Px6tWrEpt/uLq6YtWqVcI3tt26dVO5fGEzrJcvXyqdb2FhAYlEgvv37yttClYcKysr/Pvvv8L+nz179k734khMTJRrXgPIrsSVtRnMixcvMHv2bPTs2VPhxFuZ92kG1rRpU8THx8tNi4+Px2effaZ0+cJj/3bfhsL3SGETLnt7e4X3T3JysvB+fh9OTk7w8/PDjRs3kJCQINcXrFGjRkhKSir1SFr169dXGDr5//7v/4Thv8vyOSgMRG9/8QAAJiYmaNeuHQ4dOoTjx48LTdZK613+HhRlY2ODo0eP4uXLl0qvrujo6KhsfmdhYQFtbW1cu3YNHTt2BCDrH/Ts2TOF9zsRFe+jGA2soKAAU6ZMQf/+/RVGCyKq6l6+fIknT57g8ePHuHLlCrZv346WLVsq/ZayYcOGaNGiBZYtW4bExERcuHABe/bsUbn9GjVqCCfPFy5cQGpqKv766y9s3rxZrtnX8OHDERUVhT179uDBgwe4fv06Dh8+XOrn0apVK/z555/4559/cPv2bSxdulSu2VK9evVgYmKC0NBQJCcnIzY2FseOHVO5TTs7O1SvXl3upoIWFhY4d+4cbt68KXSkLwwyALBnzx4cO3YMDx48wO3bt/H3338Lw+t++umnSE9Px48//ojbt2/j7t272LVrV6mf49WrV9GgQYNim98U6tatGwoKCrB27Vr06NFD7kTqxo0bmD17Nv766y8kJSXh3r172LhxI+7cuSOcMBWlra0NBweHMndo9/T0xIEDB3D8+HHcvXsXS5cuLbajcaHo6GjExcUhKSkJd+7cQUhICK5cuSL0TSh09epVoaN+aQUHB0MkEsHb2xtPnjwR/hXH1NQU9erVK/ZfcVdJAFlfnhs3buDXX3/FvXv3sHr1ajx//lyuedQXX3yBP//8E4DsBL9u3bpYtmwZEhIScP/+fQQHB6NatWpCEy8PDw88fPgQGzduxIMHD3Dw4EGcOnVK6EtSnJ9//hkLFy5UuUyLFi1gbGyMxYsXC31CCg0ZMgRnzpzBpk2bcPfuXdy5cwdHjhzBgQMHlG5r4MCB+PPPPxEWFoakpCRs3rwZCQkJwpDaZfkc1KpVC9WqVcPff/+NrKwsuUDt6uqKnTt3olq1amjfvr3K51fUu/w9KKpnz56QSCSYO3cubty4gaSkJERERCA7OxsAUKdOHVy7dg0PHz5U+EIDkIVdNzc3hISE4PLly8Lfk/bt28v1zyIi1T6KKytaWlqYOXMm/vnnH7XcAZtIk+3duxd79+6FSCSCsbEx2rZtK3TmVMbf3x9BQUH4+uuvYWVlBW9vbyxevFjuhL2oqVOnYuvWrVi+fDkePXqEOnXqYMCAAXJNoxwdHfHtt99i+/bt2LBhA4yNjYURrEqjX79+uHXrFmbPng19fX2MHz9ebhhdHR0d+Pv7Y+XKlRgzZgwcHR3h7e2tchhQHR0d9OzZE8eOHRNOGEeMGIEHDx5gypQpMDY2xrhx4+SaKOnp6WHbtm1ITk6Grq4uHBwcMGXKFACyb3OXLVuGNWvWYMKECZBIJELH2tI4duwY+vTpU+JyEokEPXr0wJEjRxTuEVOnTh2YmJhg48aNSE9Ph46ODurXr485c+bIXfkqytXVFVu3bsWIESNKXW+fPn3w4MEDLF26FNra2hgyZIjKcADImtaEhoYiPT0dYrEYNjY2+PHHH9G4cWNhmVu3buHZs2fo0qVLqWsBgCtXriA9PV2hSV7RYXvLQ/369REYGIh169Zh69atsLa2Fu6NUigpKUloZlitWjUsWrQI69atg5+fH/Lz89G4cWO5derVq4eFCxdizZo12Lt3L+rVq4dZs2apvNEkADx+/FjoU1UckUiETz75BLt378bEiRPl5jVp0gRLly7Fpk2bsHv3buF1UXbzTkAWfL799lv88ssvWLNmDaysrDBv3jwhAJXlcyAWizF+/Hj88ssvWLlyJb744gvhJotdu3aFrq4uevbsWeZ+a+/y90BZbUuWLMFPP/0EX19faGtro1mzZkITuiFDhmDRokUYNmwYatWqpTB0NwD4+Pjgp59+gr+/P/Lz8+Ho6Ki02SYRFU8UFxen9oaRJ06cwIEDB3Dz5k3k5uYq3AkbAHbs2IH9+/fj2bNnaNeuHfz8/BQ6Th45cgTJycnCze6I6P3FxsZi+fLliIqKKvFb88ooLS0NPj4+2LZt2we787UyqampmDhxIn755Re13KgzPz8fY8aMwbRp08oUsCrCkiVLULduXZXBKSgoCPn5+WVuHvQxCA0NxYULF/DTTz+pu5QKlZmZicGDB2P9+vUcGYuoCtOIM49Xr16hbdu2xX6Lc/jwYWzbtg1TpkzB6tWrkZubi8DAwA9cJVHV8Pfff+Ps2bNITU3F6dOnsXHjRvTs2fOjDCqAbKjiiRMnlstNL9/Ho0ePMGPGDLUEFUDWFGzGjBnCwAjqUlBQgHr16hXb9+NtcXFxcHNze+8O6JXFs2fP4Obmhl9//VXdpVSogoICPHz4EJs2bULTpk0ZVIiqOI1oBlbYxvfixYtK54eFhWHQoEHo3r07AODbb7/FsGHDkJCQINd0gIjeX15eHtatW4f09HQYGRmha9eu+Oqrr9RdVoXq2bOnuktQOTTth1L0vjfqoKWlhWHDhpW43Pjx44UrL+8zBHFlUr16dWzcuBEAVDbLrOwyMjLg7e0NCwsLhTvZE1HVoxFhRZW8vDwkJibKtbG3sLCAubk5rl27xrBCVM46d+6Mzp07q7sMIpWMjY1VdoD/GGlpaZV6xK7KzNzcvEL6GRFR5aTxYSUnJwcFBQUK/ykZGRnJjXG/YMEC3Lp1C69evcK1a9eUjkpTUFCAx48fQyKRlMu9AYiIiOh/8vLySry/ERHR26RSKV6/fo26desqbXKu8WGltDdGKk0ny8ePH+Pzzz9/35KIiIhICbFYrPLeNUREyojFYmzevFnp1WONDyuGhobQ0tJCZmam3PSsrCzhpk+lJZFIAMiGk3yXOyPTh+fv71/i/QM+VpXtuWtSvR+6loreX0Vsv7y2WR7b0aT3Dqmm6rXKycmBlZUV/48lojIp/Nvx9r3T3qbxYUUsFqNRo0a4ePGicIOw1NRUpKWllTj2fFGFTb8MDAz4h7SS8PDwqLKvVWV77ppU74eupaL3VxHbL69tlsd2xGKxxrx3SLXSvN78P5aI3kVxXTQ0Iqzk5OQgIyNDuKFjQkICtLW1Ua9ePUgkEgwYMACrV6+GnZ0d6tatizVr1qBVq1bsXF8FuLq6qrsEtalsz12T6v3QtVT0/ipi++W1TU163ani8fUmog9NI8LKqVOnsHjxYuHxhAkTAADLly9HmzZt0LdvX2RmZmLFihXCTSG/+eYbdZVLRETliCfARERUHI24g/2HkpubC3d3d2RnZ/MSNRERUTnKycmBoaEh/48lojIp/Nuxe/dupffN+jhvSV0Cf39/REdHq7sMIiIiIiJSoUqGlYULF7LZAREREakUGhoKkUgk/NPW1oalpSW8vLxw69YtuWUDAgJgbW1dbvsWiUSYOXOm0nmjRo2CpaVlmbYXEBAAkUikUffBefDgAfT19XHhwgVh2ieffIKuXbu+8zatra0REBAgN23z5s2wtbWFWCwWRpIdNWqU3Gtb+M/X11du3QEDBsDHx6dU+z527JjctiQSCSwtLdG3b19s3LhR6bDeIpEIoaGhpX5+x44dQ0BAAAoKCuSm3717FyKRCMeOHSv1tsoiICAAf/zxR4VsuyQa0WeFiIiISFPt2bMHlpaWyM/PR2JiIubNm4eePXvi6tWrQpO3sWPHon///mqutHKZNWsWnJ2dhdFeK0JKSgq++uorDBs2DFu2bIGenp4wz9TUFBEREXLL161bV+5xQEAAHB0dMXXqVDRp0qRU+1y1ahUcHR3x+vVrpKSkIDY2FhMnTsTq1asRGxsr19Tp9OnTaNSoUamfz7FjxxAYGIiZM2fK3UCxbt26OH36dJlHyi2twMBA/PDDD3B2dq6Q7avCsEJERESkQps2bYQRSLt06QILCwv07t0bp06dQp8+fQAAlpaWZb7aUZWlpaVh+/btCA8Pr9D93Lx5E/n5+Rg5cqTCFRuxWIyOHTuqXL9NmzZo27YtVqxYgbVr15Zqn82aNZPb7pAhQzBmzBg4Oztj9OjRiIyMFOaVtP/S0tXVLbdtaZoq2QyMiIiI6F0VXk15u1mPsmZgt2/fRt++fVG9enWYmZnBz88PGzZsKPZ+Eu/j4cOH8PHxgZWVFXR1dWFlZYURI0bg1atXcsslJibC1dUV+vr6sLGxwcqVK+XmF9ecbdSoUfjkk0+Ex4VN5IoqbXO40NBQGBoawsXFRen86OhotG7dGtWrV4eDg4PSJkg//fQTrK2toaenhw4dOuDUqVMKNTs5OQEAevbsCZFIpNBErDS8vb3x66+/4sWLF2Vet1CnTp3w9ddfIyoqComJicL0os3Azp8/j969e6N27dqoXr06bGxshGZoAQEBCAwMBABUq1ZNaG4GKG8GVtikrjTH8vjx4+jduzcMDQ2hr6+P1q1bY9OmTUKNALBgwQJhn+9yHN8VwwoRERF9UPHxQGCg7Gdl2HZ+fj7evHmDV69e4fr16/D394e5ublwIqxMXl4eevfujUuXLmHNmjUIDQ3FnTt3MH/+/PIr7P/LyspC586dsWfPHvj5+eHQoUNYsmQJXr9+rdBPYuDAgXBxcUF4eDg8PDzg6+uL2NjYcq+pJEeOHEHHjh2ho6PYyCcxMRHTp0/HjBkzsH//fpiYmGDAgAF48uSJsMyWLVswZcoUODk54cCBA/jiiy8wZMgQZGZmCsvMmjULq1atAgCEhITg9OnTGDt2rDA/IyMDJiYm0NHRgZ2dHZYuXarQFwQAunXrhqdPnyqEobLq27cvAODkyZNK5z979gyurq7Q1tZGaGgoDh06hNmzZwv9jMaOHYsxY8YAAP766y+cPn0ap0+fVrnP0hzL8PBw9OzZE2/evMH69esRHh6O0aNH4969ewAg7GPUqFHCPt8+jhWNzcCIiIjog4mPB9zdgZQUYMMGICoKcHDQ7G03bdpU7rGFhQUiIyNRs2bNYtfZunUrbt++jbNnz6JDhw4AADc3N7Rp0wZJSUnvX9Rbli9fjtu3b+Off/5B69athene3t4Ky37zzTf48ssvAQC9evXC0aNHsXv3bvTu3btca1JFKpXi3LlzmDZtmtL5jx49wokTJ2BrawsAaNu2LerWrYvDhw9j2LBhKCgowJw5c+Dq6ootW7YI65mamsLLy0t43KhRIzRr1gwAYG9vL9dMqk2bNmjXrh2aN2+Oly9fIiwsDN999x0SEhKwfv16uXpatWoFLS0tnDlzBj179nzn521lZQUASE1NVTr/33//RWZmJpYsWYJWrVoJ00eNGgVAvqnhf/7zH6VBr6iSjqVUKsXUqVPRtm1b/PHHH8JVlF69egnbKDxu9erVU0tTsyp5ZYVDFxMREalHRISgrKOsAAAgAElEQVQsTACyn28139fYbYeFheH8+fM4d+4cDhw4AHt7e7i5ueH69evFrnP69GnUr19fCCqArDnNoEGDyqeot8TExKBDhw5yQaU4/fr1k3vcsmXLcg9PJcnKysKLFy+U3lMDAGxtbYWTawAwMzODmZmZUOeDBw+QlJSEzz//XG69QYMGleoEHgB8fX0xefJkODs7o2/fvvj5558xdepU/Pzzz0hISJBbtlq1ajA0NCw2ZJSWVCq7tWFxzQBtbW1hZGSE8ePHY8eOHXjw4MF77a9wm6qO5c2bN3Hv3j2MGTOmQponlocqGVY4dDEREZF6eHoCFhay3y0sAA8Pzd92ixYt0L59ezg6OqJ///6IiIiAVCpV2W4/NTUVZmZmCtPr1KlTqn1qa2sjPz9f6bz8/Hy5k/JHjx6hXr16pdpurVq15B7r6uri5cuXpVq3vBTuT1dXV+n8ojUWLlu4XmFoKHosdXR0ULt27Xeua+jQoZBKpTh//rzCPIlE8l59VgDg/v37ABRHHCtkaGiIuLg4mJubY/z48bCyskLLli2xf//+d95nScfy0aNHAFDq9486VMmwQkREROrh4CBrnhUYWL5NwCp622+TSCSwsbHB5cuXi12mbt26yMjIUJienp5eqn2YmZkhpfAyUREpKSlyJ+omJiZITk4u1XZLoqenp9ApH/jfSe3bywFQWLbocsoUBoq3+5eUReHJftFj+ebNGzx+/PidtglA6K+i7ArDkydPYGJi8s7bBoBDhw4BkI0oV5w2bdogLCwMWVlZOH36NBo0aIDBgwfj6tWr77Xv4hQ+p/J6/1QEhhUiIiL6oBwcgNmzKyZMVOS2Cz1//hyJiYnFNmMCZKM/3b9/H+fOnROmSaVS7Nu3r1T7cHJyQnR0NHJzc+Wmp6Wl4dSpU3Kd+11cXHDu3DlcunSpjM9EUYMGDZCRkSEXtJ48eaLQubxBgwYAgCtXrgjTXr9+jZiYmBL3IRaLYWNjg9u3b79TjZaWlrCyssLu3bvlpu/bt++9bnq5Y8cOiEQiuaZ7gOxKzsuXLxX6LpXF6dOnsX79enh6esLGxqbE5bW1tdGxY0csXLgQBQUFQpPDwqtR73uVp5CdnR2sra2xadMmoZmaMmKxuNz2WVbsYE9ERESkwsWLF/Ho0SNIpVKkpqZi9erVePLkCSZPnlzsOiNHjkRQUBAGDhyIBQsWwNTUFBs3bkRmZmap+gbMmjULERER6Ny5M/z8/GBlZYXExEQsXLgQRkZGcndanzZtGnbs2IFevXph5syZaNmyJR49eoTw8HCsW7dO5UAARbm5uaFmzZrw9vbGjBkzkJOTgyVLlihsw9HREY0aNcK4ceMQGBiIgoIChISEKL1LuzLdu3fH2bNnS13X27S0tDBnzhyMHTsWX375Jby8vJCQkIBFixYJw0qrcu/ePYwcORJDhw6FjY0Nnj9/jrCwMGzduhU+Pj4KYaKwzu7du5eqvuvXr6NGjRp48+YNUlNTERMTg23btsHe3l4YDliZqKgobNiwAQMGDEDDhg2Rm5uLVatWoWbNmujUqRMACDd9XLZsGdzc3KCtrY327duXqi5lRCIRVqxYgYEDB8LZ2RkTJkyAqakprl+/joyMDGGoZHt7exw8eBB9+vSBsbExLCwsYFHY5rKCMawQERERqTB48GDhd1NTU7Ro0QJHjhxR2f9VLBYjJiYGkydPxoQJE1CjRg0MHToU//nPfxAUFFTiPps2bYozZ85gzpw58PPzQ1ZWFkxMTODi4oK5c+fC3NxcWNbIyAgnT57EzJkzERQUhMePH6NOnTpwdnaGWCwu03M1MjJCVFQUpk2bhoEDB8LS0hKzZ8/G77//jrt37wrL6ejoIDw8HBMnTsTQoUNRq1Yt+Pr6okuXLnL3DSnOkCFDhOGcGzZsWKYaAWDMmDF49uwZgoOD8dtvv6FFixbYuXMnhg8fXuK6NWvWhJGRERYsWID09HRoaWmhWbNmCAkJwYQJExSWj4yMFMJZaUyZMgWA7CpI7dq10bp1a4SEhGDEiBEqXw9bW1tIJBLMmzcPqampqFmzJhwdHREbGyuMAubu7g4fHx+sWbMGc+fOhVQqVXlFpDT69++P2NhYzJs3TxgauVGjRnKBePXq1ZgyZQo8PDzw6tUrzJkz54Pda0UUFxf3fs+wEsnNzYW7uzuys7NLlbyJiIiodHJycmBoaMj/Y0vQt29fvHr1CkePHlV3KWpVUFAAOzs7jBw5ErNmzVJ3OcV68eIFLCwsEBwcLAz5TOWr8G/H7t27lTatrJJXVvz9/eHh4cERwYiIiKjCBAcHo0aNGmjcuDGePn2KXbt24fDhw4iKilJ3aWqnpaWFefPmYerUqfDz80P16tXVXZJS69atQ506dTBixAh1l1JlVcmwsnDhQn7rQ0RERBVKV1cXwcHBuH//PqRSKZo3b469e/cq3OukqvLy8kJSUhLu3r0r9MXQNBKJBFu2bCn1/Vuo/LEZGBEREb03NgMjondRUjMwDl1MREREREQaiWGFiIiIiIg0EsMKERERERFpJIYVIiIiIiLSSAwrRERERESkkRhWiIiIiIhII1XJsOLv74/o6Gh1l0FERERERCpUybCycOFC3r2eiIiIVAoNDYVIJBL+aWtrw9LSEl5eXrh165bcsgEBAbC2ti63fYtEIsycOVPpvFGjRsHS0rJM2wsICIBIJMKbN2/Ko7xy8eDBA+jr6+PChQvCtE8++QRdu3Z9521aW1sjICBAbtrmzZtha2sLsVgMIyMjALJj+PZrW/jP19dXbt0BAwbAx8en1Ps/fvw4+vTpAwsLC+jp6cHS0hJ9+vTBr7/+WubncvfuXYhEImzcuLHEZa2trTF8+PAy70PV9ooeR3Xh7TiJiIiIVNizZw8sLS2Rn5+PxMREzJs3Dz179sTVq1eFG2COHTsW/fv3V3OllcusWbPg7OyMdu3aVdg+UlJS8NVXX2HYsGHYsmUL9PT0hHmmpqaIiIiQW75u3bpyjwMCAuDo6IipU6eiSZMmKvcVFhaGQYMGwdPTE6tXr0atWrVw7949xMbG4tChQxg2bFj5PbEqhGGFiIiISIU2bdqgcePGAIAuXbrAwsICvXv3xqlTp9CnTx8AgKWlZZmvdlRlaWlp2L59O8LDwyt0Pzdv3kR+fj5GjhypcMVGLBajY8eOKtdv06YN2rZtixUrVmDt2rUqlw0ODkabNm0QFhYGkUgkTB85ciQKCgre/UlUcVWyGRgRERGpV0V2Ha3obqmFV1Py8vKEacqagd2+fRt9+/ZF9erVYWZmBj8/P2zYsEHuRLa8PHz4ED4+PrCysoKuri6srKwwYsQIvHr1Sm65xMREuLq6Ql9fHzY2Nli5cqXc/OKas40aNQqffPKJ8LiwiVxRpW0OFxoaCkNDQ7i4uCidHx0djdatW6N69epwcHDAH3/8obDMTz/9BGtra+jp6aFDhw44deqUQs1OTk4AgJ49e0IkEr1T0yZvb2/8+uuvePHihcrlMjMzUadOHaXHRUvrf6fcZT12r1+/xowZM2BmZobq1avDw8MD9+7dU1rD9u3bYWdnhxo1aqBLly64dOmS3Pzw8HC4urrC3Nwc+vr6aNmyJVauXFmqMBUTE4NOnTpBIpHAyMgIAwcOVGgOGR0djc6dO8PQ0BA1atRAkyZNMHfu3BK3rQrDChEREX1wlSms5Ofn482bN3j16hWuX78Of39/mJubCyfCyuTl5aF37964dOkS1qxZg9DQUNy5cwfz588v3+IAZGVloXPnztizZw/8/Pxw6NAhLFmyBK9fv5YLVAAwcOBAuLi4IDw8HB4eHvD19UVsbGy511SSI0eOoGPHjtDRUWzkk5iYiOnTp2PGjBnYv38/TExMMGDAADx58kRYZsuWLZgyZQqcnJxw4MABfPHFFxgyZAgyMzOFZWbNmoVVq1YBAEJCQnD69GmMHTtWmJ+RkQETExPo6OjAzs4OS5cuVXrS3q1bNzx9+lQhDBXVoUMHxMTEYNasWbh27RqkUmmZj4syCxcuxL///ostW7YgJCQEf//9N1xdXRX6H504cQIrV67EwoULsXPnTjx//hz9+/eXW64wrG7evBkHDx7EuHHjsGDBAsyaNUtlDTExMejbty8MDAywa9curFmzBpcvX0aXLl2QmpoKQBbOPT090bBhQ+zatQsRERGYPn06cnNz3+v5sxkYERERkQpNmzaVe2xhYYHIyEjUrFmz2HW2bt2K27dv4+zZs+jQoQMAwM3NDW3atEFSUlK51rd8+XLcvn0b//zzD1q3bi1M9/b2Vlj2m2++wZdffgkA6NWrF44ePYrdu3ejd+/e5VqTKlKpFOfOncO0adOUzn/06BFOnDgBW1tbAEDbtm1Rt25dHD58GMOGDUNBQQHmzJkDV1dXbNmyRVjP1NQUXl5ewuNGjRqhWbNmAAB7e3u5Jl9t2rRBu3bt0Lx5c7x8+RJhYWH47rvvkJCQgPXr18vV06pVK2hpaeHMmTPo2bNnsc8rKCgICQkJmD9/PubPnw8DAwM4Oztj+PDhGDRoUNkP1P9Xs2ZNhIeHC1dn7Ozs0LVrV/zyyy8YPXq0sNzTp09x6dIlGBsbAwDMzc3h6OiIs2fPokuXLgCA6dOnC8tLpVJ07doVxsbGmDJlCubPn1/sVb8ffvgBNjY2OHToELS1tQEAnTp1gp2dHZYtW4Yff/wR//zzD/Ly8rB27Vrh6qOzs/M7P+9CDCtERET0QURH/++qx/Ll/5vu6ir7p6nbDgsLg6WlJaRSKVJSUrB69Wq4ubnhxIkTwslwUadPn0b9+vWFoALIRvgaNGgQLl++/H4FFRETE4MOHTrIBZXi9OvXT+5xy5Ytyz08lSQrKwsvXryAqamp0vm2trZCUAEAMzMzmJmZCXU+ePAASUlJCk26Bg0apPRKjTJFR/3q27cvatSogZUrV2LGjBlCHyUAqFatGgwNDYUrCMUxMzPDiRMncP78eRw5cgRnz55FbGwsDhw4gPHjx2PdunWlqq2ozz77TK4ZWZcuXWBpaYlTp07JhZVOnToJQQWQvbYA5F7ftLQ0zJ07F4cOHUJycrLcVZf09HSYm5sr7D83NxcXLlzADz/8IAQVAGjYsCG6dOmCY8eOAZAFwGrVqsHLywvjxo1D165di32Ny4JhhYiIiD6IosEhOLhybLtFixZyJ68uLi6wsrJCQEAAdu3apXSd1NRUmJmZKUyvU6dOqfapra2N/Px8pfPy8/PlTsofPXpUqqACALVq1ZJ7rKuri5cvX5Zq3fJSuD9dXV2l84vWWLhs4XqFoaHosdTR0UHt2rXfua6hQ4dixYoVOH/+vNzrDQASiaTEPiuFHB0d4ejoCEAWzAYPHoz169dj0qRJaNGiRZnrUvaeqVOnDpKTk+WmKXttgf8db6lUCk9PT2RmZmL27NmwtbWFRCLBuXPnMHHixGLfB5mZmZBKpUqDjLm5Oc6ePQsAaNy4MaKjo7Fo0SJ4e3sjLy8PHTp0wOLFi9GjR48yP+9C7LNCREREVAYSiQQ2NjYqr5DUrVsXGRkZCtPT09NLtQ8zMzOkpKQonZeSkiJ3AmtiYqJw4vqu9PT0FDrlA7JAVHQ5AArLFl1OmcJA8Xb/krIoHF646LF88+YNHj9+/E7bBCD0V1HWFOrJkycwMTEp8zaNjIwwefJkAMC1a9cAlP3YKXvPpKeno169emWqJTExEefPn8fGjRsxevRodOvWDe3bty82NBYyNjaGSCRCWlqawry0tDS5gOjk5ISYmBjk5OTg999/h5aWFvr161eq90VxqmRY4R3siYiI1Ksi781c0fd9fv78ORITE1U2cenUqRPu37+Pc+fOCdOkUin27dtXqn04OTkhOjpaoXNyWloaTp06Jde538XFBefOnVMY+eldNGjQABkZGXJB68mTJwqdyxs0aAAAuHLlijDt9evXiImJKXEfYrEYNjY2uH379jvVaGlpCSsrK+zevVtu+r59+97rppc7duyASCSSa7oHyK7kvHz5UqHvUlHKTuYB4N9//wXwv5BV1mO3d+9euY7/J0+exIMHD9CpU6cSnpG858+fA5C/oiWVSrFp0yaV6+nr66N9+/bYs2eP3NW+e/fu4dSpU3KjxBUSi8VwdnbGf//7X+Tm5uLOnTtlqvVtVbIZ2MKFC4WOP0RERPThVaawcvHiRTx69AhSqRSpqalYvXo1njx5InxjrszIkSMRFBSEgQMHYsGCBTA1NcXGjRuRmZlZqqGLZ82ahYiICHTu3Bl+fn6wsrJCYmIiFi5cCCMjI7k+F9OmTcOOHTvQq1cvzJw5Ey1btsSjR48QHh6OdevWqRwIoCg3NzfUrFkT3t7emDFjBnJycrBkyRKFbTg6OqJRo0YYN24cAgMDUVBQgJCQEIXRx4rTvXt3oflQWWlpaWHOnDkYO3YsvvzyS3h5eSEhIQGLFi0q1fndvXv3MHLkSAwdOhQ2NjZ4/vw5wsLCsHXrVvj4+MDGxkZu+cI6u3fvrnK7rq6usLa2hoeHB2xtbfHixQucOHECy5YtQ6dOnYRO7mU9dk+fPsWAAQMwfvx4PHz4EN9//z1sbW3xxRdflOZwCZo2bQpra2tMmDABgYGBEIlEWLdundIrgEXNnz8fbm5ucHd3h4+PD549e4Y5c+bAyMgIfn5+AIB169bhxIkT6Nu3L6ysrPDo0SMsWrQIFhYW79T8rVCVvLJCREREVFqDBw9Gp06d0LlzZ0yYMAGAbOjdwYMHF7uOWCxGTEwMWrVqhQkTJmDkyJGwsrLCxIkTYWhoWOI+mzZtijNnzsDW1hZ+fn5wcXHBrFmz0K1bN5w5c0au/4CRkRFOnjyJTz/9FEFBQejTpw/8/Pygo6MDsVhcpudqZGSEqKgo5OTkYODAgZg5cyZ8fX0VhmnW0dFBeHg4DA0NMXToUEyZMgVubm4YNWpUqfYzZMgQXL169Z2/cR8zZgxWrFiBP/74A/3798eWLVuwc+dOuQ7mxalZsyaMjIywYMECuLu7w8vLC5cvX0ZISAh++uknheUjIyOFgKGKv78/dHV1sWjRIri5ueHTTz/FgQMH4Ofnh+joaKGTfFmP3ffff4/GjRtj1KhR8PHxQdu2bREdHY1q1aqVfKDeIhaLERERIez3q6++gp2dnTC8syouLi44dOiQ0AdnwoQJaN68OU6ePClcMWrdujVyc3Px/fffw8XFBZMmTULDhg3xxx9/QCKRlKnWt4ni4uLKZxDoSiA3Nxfu7u7Izs7mlRUiIqJylJOTA0NDQ/4fW4K+ffvi1atXOHr0qLpLUauCggLY2dlh5MiRJd7jQ51evHgBCwsLBAcHC0M+U/kq/Nuxe/dupU0rq2QzMCIiIqKKFhwcjBo1aqBx48Z4+vQpdu3ahcOHDyMqKkrdpamdlpYW5s2bh6lTp8LPzw/Vq1dXd0lKrVu3DnXq1MGIESPUXUqVxbBCREREVAF0dXURHByM+/fvQyqVonnz5ti7d6/CvU6qKi8vLyQlJeHu3buwt7dXdzlKSSQSbNmypdT3b6Hyx2ZgRERE9N7YDIyI3kVJzcDYwZ6IiIiIiDQSwwoREREREWkkhhUiIiIiItJIDCtERERERKSROLQBERERlZucnBx1l0BElUhJfzOqZFjx9/eHh4cHXF1d1V0KERHRR0EsFsPc3BxWVlbqLoWIKhmxWFzs8NBVMqwsXLiQwyoSERGVIz09Pdy5cwd5eXnqLoWIKpGcnBxYWVlBLBYrnV8lwwoRERGVPz09Pejp6am7DCL6iLCDPRERERERaSSGFSIiIiIi0kgMK0REREREpJEYVoiIiIiISCMxrBARERERkUZiWCEiIiIiIo3EsEJERERERBqJYYWIiIiIiDQSwwoREREREWkkhhUiIiIiItJIDCtERERERKSRGFaIiIiIiEgjMawQEREREZFGqpJhxd/fH9HR0eoug4iIiIioyoqOjoa/v7/KZURxcXHSD1SP2uXm5sLd3R3Z2dkwMDBQdzlERERERFVaTk4ODA0NERUVBX19fYX5VfLKChERERERaT6GFSIiIiIi0kgMK0REREREpJEYVoiIiIiISCMxrBARERERkUZiWCEiIiIiIo3EsEJERERERBqJYYWIiIiIiDQSwwoREREREWkkhhUiIiIiItJIDCtERERERKSRGFaIiIiIiEgjMawQEREREZFGYlghIiIiIiKNxLBCREREREQaiWGFiIiIiIg0EsMKERERERFpJIYVIiIiIiLSSAwrRERERESkkRhWiIiIiIhIIzGsEBERERGRRqqSYcXf3x/R0dHqLoOIiIiIqMqKjo6Gv7+/ymVEcXFx0g9Uj9rl5ubC3d0d2dnZMDAwUHc5RERERERVWk5ODgwNDREVFQV9fX2F+VXyygoREREREWk+hhUiIiIiItJIDCtERERERKSRGFaIiIiIiEgjMawQEREREZFGYlghIiIiIiKNxLBCRCrFxwOBgbKfRERERB8SwwoRFSs+HnB3BwICZD8ZWIiIiOhDYlghomJFRAApKbLfU1KAyEj11kNERERVC8MKERXL0xOwsJD9bmEBeHiotx4iIiKqWnTUXQARaS4HByAqSnZFxcND9piIiIjoQ2FYISKVHBwYUoiIiEg92AyM6CPFUbyIiIiosmNYIfoIaeIoXgxPREREVFYMK0QfIU0bxUsTwxMRERFpPoYVoo+Qpo3ipWnhiYiIiCoHhhWij1DhKF6BgbKf6u4gr2nhiYiIiCoHjgZG9JHSpFG8OAQyERERvQuGFSL6IDQpPBEREVHlwGZgRERERESkkRhWiIiIiIhIIzGsEBERERGRRmJYISIiIiIijcSwQkREREREGolhhYiIiIiINBLDChFVKvHxsptdxseruxIiIiKqaAwrRFRpxMcD7u5AQIDsJwMLERHRx41hhYgqjYgIICVF9ntKChAZqd56iIiIqGIxrBBRpeHpCVhYyH63sAA8PNRbDxEREVUsHXUXUF5iY2Nx4MABAMDo0aPRrl07NVdEROXNwQGIipJdUfHwkD0mIiKij9dHEVaePXuGHTt2YN26dXj58iV8fX2xceNGaGtrq7s0IipnDg4MKURERFWFRoSVEydO4MCBA7h58yZyc3Px+++/KwSNHTt2YP/+/Xj27BnatWsHPz8/1KpVCwBw/fp1tGjRArq6utDV1UWdOnXw4MEDNGjQQB1Ph4iIiIiIyoFG9Fl59eoV2rZtC29vb6XzDx8+jG3btmHKlClYvXo1cnNzERgYKMzPzs5GzZo1hcc1atRAdnZ2hddNREREREQVRyOurPTu3RsAcPHiRaXzw8LCMGjQIHTv3h0A8O2332LYsGFISEhA48aNYWBggKdPnwrLP3v2DAYGBhVfOBERERERVRiNuLKiSl5eHhITE+HwViN1CwsLmJub49q1awAAe3t7XL16FXl5ecjJyUF6ejqsrKzUVTIREREREZUDjbiyokpOTg4KCgpgbGwsN93IyAhZWVkAZM2+vLy84OvrC5FIhIkTJ7JzPRERERFRJafxYUUqlZZqOVdXV7i6upZqWX9/f4jF4jKvR0RERERE7yc6OhrR0dEAZK2oVNH4sGJoaAgtLS1kZmbKTc/KyoKRkdE7bXPhwoXs00JEREREpAZvXyzIyclBSEhIsctqfJ8VsViMRo0ayXW+T01NRVpaGuzt7dVYGVHFiI8HAgNlP4mIiIiqMo24spKTk4OMjAwkJycDABISEqCtrY169epBIpFgwIABWL16Nezs7FC3bl2sWbMGrVq1QuPGjdVcOVH5io8H3N2BlBRgwwbZ3dp5A0QiIiKqqjQirJw6dQqLFy8WHk+YMAEAsHz5crRp0wZ9+/ZFZmYmVqxYIdwU8ptvvlFXuUQVJiJCFlQA2c/ISIaVihQfLzvmnp48zkRERJpIFBcXV7oe7B+B3NxcuLu7Izs7m31WSCO9fWXFwoJXVioSjzUREZH65eTkwNDQEFFRUdDX11eYr/F9ViqCv7+/MAIBkSZxcJCdNAcG8uS5oim7ikVEREQfTnR0NPz9/VUuwysrRFQl8coKERGR+pV0ZUUj+qwQEX1ohVexIiMBDw8GFSIiIk3EsEJEVZaDA0MKERGRJquSfVaIiIiIiEjzMawQEREREZFGYlghIiIiIiKNVCXDCocuJiIiIiJSLw5dXASHLiYiIiIi0hy8KSQREREREVVKDCtERERERKSRGFaIiIiIiEgjMawQEREREZFGYlghIiIiIiKNVCXDCocuJiIiIiJSLw5dXASHLiaiihAfD0REAJ6egIODuqshIiKqPDh0MRFRBYqPB9zdgYAA2c/4eHVXRERE9PFgWCEieg8REUBKiuz3lBQgMlK99RAREX1MGFaIiN6DpydgYSH73cIC8PBQbz1EREQfEx11F0BEVJk5OABRUbIrKh4e7LNCRERUnhhWiIjek4MDQwoREVFFYDMwIiIiIiLSSAwrRERERESkkapkWOFNIYmIiIiI1Is3hSyCN4UkIiIiItIcvCkkERERERFVSgwrRERERESkkRhWiIiIiIhIIzGsEBERERGRRmJYISIiIiIijcSwQkREREREGolhhYhIQ8THA4GBsp9ERETEsEJEpBHi4wF3dyAgQPaTgYWIiKiKhhXewZ7KG78Rp/cVEQGkpMh+T0kBIiPVWw8REVFF4x3si+Ad7KkiFH4jnpICWFgAUVGAg4O6q6LKhu8jIiKqingHe6IKxm/EqTw4OMgCSmAggwoREVEhHXUXQFTZeXoCGzb87xtxDw91V0SVlYMDQwoREdHbGFaI3lPhN+KRkbKgwpNNIiIiovLBsEJUDviNOBEREVH5Y58VIiIiIiLSSAwrRERERESkkRhWiIiIiIhIIzGsEBERERGRRmJYISIiIiIijcSwQkREREREGqlKhhV/f39ER1lKwKsAACAASURBVEeruwwiIiIioiorOjoa/v7+KpcRxcXFST9QPWqXm5sLd3d3ZGdnw8DAQN3lEBERERFVaTk5OTA0NERUVBT09fUV5lfJKytERERERKT5GFaIiIiIiEgjMawQEX1k4uOBwEDZTyIiosqMYYWI6CMSHw+4uwMBAbKfDCxERFSZMawQEX1EIiKAlBTZ7ykpQGSkeushIiJ6HwwrREQfEU9PwMJC9ruFBeDhod56iIiI3oeOugsgIqLy4+AAREXJrqh4eMgeExERVVYMK0REHxkHB4YUIiL6OLAZGBERERERaSSGFSIiIiIi0kgMK0REREREpJEYVoiIiIiISCMxrBARERERkUZiWCEiIiIiIo1UJcOKv78/oqOj1V0GEREREVGVFR0dDX9/f5XLiOLi4qQfqB61y83Nhbu7O7Kzs2FgYKDucoiIiIiIqrScnBwYGhoiKioK+vr6CvOr5JUVIiIiIiLSfAwrRERERESkkRhWiIiIiIhIIzGsEBGRUvHxQGCg7CcREZE6MKxQlcaTMSLl4uMBd3cgIED2k58RIiJSB4YVqrJ4MkZUvIgIICVF9ntKChAZqd56iIioamJYoSqLJ2NExfP0BCwsZL9bWAAeHuqth4iIqiaGFaqyeDJGVDwHByAqStZMMipK9piIiOhD01F3AUTqUngyFhkpCyo8GSOS5+DAzwUREakXwwpVaTwZIyIiItJcbAZGREREREQaiWGFiIiIiIg0EsMKERERERFpJIYVIiIiIiLSSAwrRERERESkkRhWiIiIiIhIIzGsEBERERGRRmJYISIiIiIijcSwQkREREREGolhhYiIiIiINBLDChER/T/27j86r7LAE/g3TVsogaQUEZrKjEgpM4VBAovAGWAsWopsknIOzngQRgFHYBiXcQ4I+rrDtniogkiLU9BBcZlFgRXlR5KxvkM1hVVw3JEwDlsGaBd/TMOPEdrERqBD0/3jPWQtLU1Kk7w3fT+fczz3fe+93Hx7kAvfPM+9DwAUUk2WlVKplHK5XO0YADWjpydZvLiyBYAkKZfLKZVKOzynrru7e8s45am6gYGBtLa2pq+vL42NjdWOA1ATenqS1taktzdpbk66upKWlmqnAqAI+vv709TUlK6urjQ0NGxzvCZHVgAYPx0dlaKSVLadndXNA8DEoawAMKba2ysjKkll29ZW3TwATByTqx0AgN1bS0tl6ldnZ6WomAIGwEgpKwCMuZYWJQWAnWcaGAAAUEjKCgAAUEjKCgAAUEjKCgAAUEjKCgAAUEjKCgAAUEjKCgAAUEjKCgAAUEjKChNST0+yeHFlCwDA7klZYcLp6UlaW5NFiypbhQUAYPekrDDhdHQkvb2Vz729SWdndfMAADA2lBUmnPb2pLm58rm5OWlrq24eAADGxuRqB4Cd1dKSdHVVRlTa2irfgdrQ01MZXW1v988+QC1QVpiQWlr8hwrUmteeV+vtTW6+ufJLC/cBgN2baWAATAieVwOoPcoKABOC59UAao9pYABMCJ5XA6g9ygoAE4bn1QBqi2lgAABAISkrAABAISkrAABAISkrAABAISkrAABAISkrAABAISkrAABAIe0WZWVwcDCXXHJJFi5cmFtuuaXacQAAgFGwWywKOWnSpPzX//pf88gjj2TdunXVjgMAAIyCqpWVBx98MPfee2+efPLJDAwMZOXKlamvr9/qnNtvvz133313Nm7cmGOOOSaXXnppZsyYsd3rvfWtbx2P2AAAwDip2jSwV155JUcffXTOOuus7R5fsWJFbrvttlxyySVZvnx5BgYGsnjx4nFOCQAAVEvVRlbmz5+fJHn00Ue3e/yee+7JmWeemZNPPjlJcvnll+fss8/OmjVrMnv27Nx3333p7u5OY2NjrrrqqnHLDQAAjI9CPrOyadOmrF27NhdeeOHQvubm5hx44IFZvXp1Zs+enYULF2bhwoVVTAkAAIylQpaV/v7+DA4OZt99991q//Tp07Nhw4bt/jVXX311nnrqqbzyyitZvXp1Pv/5z2fSpN3iZWcAjLKenqSjI2lvT1paqp0GgDdSyLKyZcuWnf5rPv3pT4/43FKplKlTpyZJFixYkAULFuz0zwNgYurpSVpbk97e5Oabk64uhQVgPJXL5ZTL5SSVGVU7Usiy0tTUlEmTJmX9+vVb7d+wYUOmT5++y9dfsmRJGhsbd/k67Dy/zQSqraOjUlSSyraz0/0IYDz99mBBf39/brzxxjc8t5DzpKZOnZpDDjlkq4fvn3nmmTz77LOZO3duFZOxK177beaiRZVtT0+1EwG1qL09aW6ufG5uTtraqpsHgDdWtZGV/v7+PP/880OLOK5Zsyb19fWZNWtWpk2bljPOOCPLly/PnDlzMnPmzNx000058sgjM3v27GpFZhf5bSZQBC0tlalfnZ2VouI+BFBcVSsrDz30UK655pqh7xdddFGSZOnSpTnqqKNy+umnZ/369Vm2bNnQopCXXXZZteIyCtrbK/PDe3v9NhOorpYWJQVgIqjr7u7e+afZJ6iBgYG0tramr6/PMytV0tPjt5kAAFT09/enqakpXV1daWho2OZ4IR+wH2ulUiltbW3eAlYFfpsJAEBSeStYZ2fnDs8xsgIAAFTFcCMrhXwbGAAAgLICAAAUkrICAAAUkrICAAAUkrICAAAUUk2WlVKplHK5XO0YAABQs8rlckql0g7P8epiAACgKry6GAAAmJCUFQAAoJCUFQAAoJCUFQB4k3p6ksWLK1sARp+yAgBvQk9P0tqaLFpU2SosAKOvJsuKVxcDsKs6OpLe3srn3t6ks7O6eQAmGq8ufh2vLgZgtLw2stLbmzQ3J11dSUtLtVMBTCzDvbp4chUyAcCE19JSKSidnUlbm6ICMBaUFQB4k1palBSAsVSTz6wAAADFp6wAAACFpKwAAACFpKwAAACFpKwAAACFVJNlxaKQAABQXRaFfB2LQgIAQHEMtyhkTY6sAAAAxaesAAAAhaSsMCI9PcnixZUtAACMB2WFYfX0JK2tyaJFla3CAgDAeFBWGFZHR9LbW/nc25t0dlY3DwAAtUFZYVjt7Ulzc+Vzc3PS1lbdPAC7G1NtAbZvcrUDUHwtLUlXV2VEpa2t8h2A0fHaVNve3uTmmyv3W/dZgAplhRFpafEvT4CxsL2ptu63ABU1OQ3MCvYAFIWptkCtsoL961jBHoAi6ukx1RaoTcOtYG8aGABUmam2ANtXk9PAAACA4lNWAACAQlJWAACAQlJWAACAQlJWAACAQtrh28CWLFmy0xe8+OKLM3369DcdCAAAIBlmZGXlypV59dVXU19fP6L/rVq1Kr/5zW/GKzsAALAbG3adlf/yX/5L9t133xFd7IEHHtjlQAAAAMkwIys333xzmpqaRnyxG2+8MQcccMAuhwIAANhhWZk9e3YmTRr5M/gHH3xw6uvrdznUWCuVSimXy9WOAQAANatcLqdUKu3wnLru7u4tI7nYhg0bhh6c37x5c3784x+nvr4+RxxxRPbaa69dTzsOBgYG0tramr6+vjQ2NlY7DgAA1LT+/v40NTWlq6srDQ0N2xwf9pmVf/u3f8unPvWp9Pb25qCDDsrnPve5fOYzn8n//b//N0nS1NSUL3zhC5k1a9bopwcAAGrWsHO8vvSlL+Xggw/OV7/61bzrXe/KJz/5yey3337p7OxMR0dH5syZk//+3//7eGQFAABqyLBl5bHHHsu5556bgw8+OB/5yEfyy1/+Mn/yJ3+SyZMnZ8qUKTnrrLPyL//yL+ORFQAAqCHDlpVNmzYNPZOyxx57ZI899siMGTOGjs+YMSMbNmwYu4QAwIj19CSLF1e2ABPdsGXlrW99a3p7e4e+//Vf/3X222+/oe8vvPCCFesBoAB6epLW1mTRospWYQEmumHLymmnnZZf//rXQ99POOGE7LHHHkPff/CDH+SII44Ym3QAwIh1dCSv/X6xtzfp7KxuHoBdNezbwM4666wdHr/gggtGLQwA8Oa1tyc331wpKs3NSVtbtRMB7JphywoAMDG0tCRdXZURlba2yneAiWzEZWX16tWZO3fuWGYBAHZRS4uSAuw+hn1mJak8l/LXf/3XY50FAABgyLAjK/fee2++9rWv5bOf/ex45AEAAEgyTFn5yle+kvvuuy+f+9zncvjhh49XJgAAgB1PA7vjjjty4YUXejUxAAAw7nZYVk477bTccsstefrpp8crDwAAQJJhysrll1+ehQsX5tJLL92tCkupVEq5XK52DAAAqFnlcjmlUmmH59R1d3dvGe5Cf//3f59bb701d91116iFq4aBgYG0tramr68vjY2N1Y4DAAA1rb+/P01NTenq6kpDQ8M2x0e0zsp//s//OW95y1tGPRwAAMAbGdE6K0ly3HHHjWUOAACArYy4rAAAAIynEZeVK6+8Mnfeeec2+7/5zW/mv/23/zaqoQAAAEZcVv75n/95u1PB3vWud+Wf//mfRzUUAADAiMvKpk2btrt/y5Ytefnll0ctEAAAQLITZeX3f//3c/fdd2+z/9vf/nZ+7/d+b1RDAQDV1dOTLF5c2QJUy4heXZwkF154YS677LKsXr0673znO5MkP/3pT/P888/nuuuuG7OAAMD46ulJWluT3t7k5puTrq6kpaXaqYBaNOKRlcMOOyxf//rXc+KJJ+aFF17Ir371q5x44om57bbbMmfOnLHMCACMo46OSlFJKtvOzurmAWrXiEdWkqSpqSnnnXfeWGUBAAqgvb0yotLbmzQ3J21t1U4E1KqdKiu9vb3p6OjIL3/5yyTJQQcdlLa2tsyaNWtMwrFrenoqvx1rbzd8D8DItbRUpn51dlaKin+HANUy4mlgq1atyoc//OGsXr06s2bNyqxZs/L444/n3HPPzapVq8YwIm/Ga/ONFy2qbD0gCcDOaGlJrrxSUQGqa8QjK3/7t3+bD3/4wznnnHO22v+Nb3wjX/7yl/Pud797tLOxC7Y339i/cAAAmEhGPLKyYcOG7RaSP/qjP0pfX99oZmIUtLdX5hkn5hsDADAxjbis/NEf/dF2p3s98MADOemkk0YzE6PgtfnGixd75SQAABPTiKeBTZ8+PXfeeWd+9KMf5fd+7/dSV1eXxx9/PD/72c/S2tqar33ta0Pnnn/++WMSlp3T0qKkAAAwcY24rDzxxBM59NBDkyRr165NkkyZMiWHHnponnjiiaHz6urqRjkiAABQi0ZcVpYuXTqWOQAAALYy4mdWAAAAxtOwIytLliwZ0YVKpdIuhwEAAHjNsCMrK1euzE9/+tNMmjQp9fX1b/g/AACA0TTsyMp5552X+++/P4888kjmz5+fU089Nb/7u787HtnGTKlUSltbWxYsWFDtKAAAUJPK5XI6Ozt3eE5dd3f3lpFcbPXq1SmXy1m1alVmzpyZ+fPn573vfW+amppGJex4GBgYSGtra/r6+tLY2FjtOAAAUNP6+/vT1NSUrq6uNDQ0bHN8xA/Yz507N3/1V3+Vb33rWznrrLPS3d2dP/7jP87GjRtHNTAAAECyE68uTpJNmzbloYceSrlcztq1a3P88cdn6tSpY5UNAACoYSMqKz/96U9TLpfzwAMP5Hd/93dz6qmnplQqZZ999hnrfAAAQI0atqx88IMfzODgYN7znvfkpptuytve9rahY4ODg0OfJ02yZAsAsLWenqSjI2lvT1paqp0GmGiGLSvPPvtskuSOO+7InXfe+Ybnfe973xu9VADAhNfTk7S2Jr29yc03J11dCguwc4YtK9dff/145AAAdjMdHZWiklS2nZ3KCrBzhi0rL774Yt71rndl7733Ho88AMBuor29MqLS25s0NydtbdVOBEw0w5aV++67L9dee20OO+ywnHDCCTn++OPz9re/fRyiAQATWUtLZepXZ2elqBhVAXbWsGXlhhtuyMaNG/PjH/84P/rRj3LnnXdmr732yvHHH5/jjz8+LS0tmTJlynhkBQAmmJYWJQV480b06uK99947p5xySk455ZQMDg5m9erVefjhh/O3f/u3eeaZZ3L00Ufn+OOPz4knnpjp06ePdWYAAKAG7NSikEnlFcVHHHFEjjjiiHz0ox/Nc889l3/8x3/MD3/4wwwMDOQDH/jAWOQEAABqzE6Xldc74IAD0t7envb29tHIAwAAkGSEZWXlypV57LHHcuSRR+aUU07J//pf/ytf//rXU19fn/e85z0588wzxzonAABQY4YtK9/61rfyta99Lccee2xuuummPPPMM7nrrrvy/ve/P4ODg7n11lszZcoUIysAAMCoGrasdHR05BOf+ETmzZuXtWvX5oILLsgnPvGJnHbaaUkq08C+/e1vKysAAMComjTcCc8991zmzp2bJDnkkEMyefLk/P7v//7Q8SOOOCK9ry1PCwAAMEqGLSsNDQ0ZGBgY+n7ooYdmr732Gvq+efPmsUkGAADUtGHLytvf/vY88cQTQ9+XL1+e/ffff+j7U089lYMOOmhs0gEAADVr2GdWrrjiitTX17/h8T322CPnnnvuaGYCAAAYvqwccMABOzx+4oknjloYAACA1ww7DQwAAKAaRlRWNm3alCuvvHKsswAAAAwZtqz8+te/zmWXXZYpU6aMRx4AAIAkw5SV5557Lh/72Mfylre8JaVSabwyAQBspacnWby4sgVqxw7Lysc+9rG84x3vyKc//ekdvhEMAGCs9PQkra3JokWVrcICtWOHZeWll15KQ0ODogIAVE1HR9LbW/nc25t0dlY3DzB+dlhWvvjFL+bHP/5xli5dOl55AAC20t6eNDdXPjc3J21t1c0DjJ8dlpV3vOMdufHGG/PYY49l2bJl45UJAGBIS0vS1VV5ZqWrq/IdqA3Dvg1s//33zw033JBf/OIX45EHAGAbLS3JlVcqKlBrhl3BPkn23nvvXHvttWOd5U174okn8jd/8zepr6/P5s2bc/HFF2fu3LnVjgUAAOyCEZWVJJk8ecSnjru3vOUtueaaa9LQ0JCnn3461157bb70pS9VOxYAALALdqqBbN68Ob/85S+zYcOGDA4ObnXs6KOP3qkf/OCDD+bee+/Nk08+mYGBgaxcuXKbt47dfvvtufvuu7Nx48Ycc8wxufTSSzNjxoxtrrXffvsNfZ46dWomTRp2dhsAAFBwIy4rPT09WbJkSV544YVtjtXV1eV73/veTv3gV155JUcffXSOOeaYfPWrX93m+IoVK3LbbbflU5/6VJqbm7N8+fIsXrw4N9xwwxte8z/+4z9y/fXX50Mf+tBOZQEAAIpnxGXlhhtuyPHHH5/zzjtvu6MbO2v+/PlJkkcffXS7x++5556ceeaZOfnkk5Mkl19+ec4+++ysWbMms2fPzn333Zfu7u40NjbmqquuyubNm/OZz3wm8+bNy3HHHbfL+QAAgOoacVl57rnncvXVV49KURnOpk2bsnbt2lx44YVD+5qbm3PggQdm9erVmT17dhYuXJiFCxcmSQYHB7NkyZIcccQRaW1tHfN8AADA2BtxWTnuuOOyevXqzJo1ayzzJEn6+/szODiYfffdd6v906dPz4YNG7Y5f9WqVfnhD3+YF154IQ899FAaGhpy9dVXj3lOAABg7OywrHznO98Z+jx37tx8+ctfzuOPP563v/3t27wd7PTTTx+1UFu2bNmp80855ZSccsopIz6/VCpl6tSpSZIFCxZkwYIFO/XzAACAN6dcLqdcLiepzKjakR2Wldtuu22r71OnTs3DDz+chx9+eKv9dXV1o1pWmpqaMmnSpKxfv36r/Rs2bMj06dN3+fpLlixJY2PjLl8HAADYOb89WNDf358bb7zxDc/dYVm54447RjfZCE2dOjWHHHJIHn300RxzzDFJkmeeeSbPPvusxR4BAKBGVG2lx/7+/jz//PNZt25dkmTNmjWpr6/PrFmzMm3atJxxxhlZvnx55syZk5kzZ+amm27KkUcemdmzZ1crMgAAMI52WFa++MUv5qMf/WimTZs2oot9+ctfzgc/+MERTbF66KGHcs011wx9v+iii5IkS5cuzVFHHZXTTz8969evz7Jly4YWhbzssstGlAMAAJj46rq7u9/wafb3vOc9ueuuu0b0nMiWLVvS3t6er3zlK2lubh7VkKNlYGAgra2t6evr88wKAABUWX9/f5qamtLV1ZWGhoZtju9wZGXLli354z/+4zELVy2lUiltbW3eAgYAAFVSLpfT2dm5w3N2OLLyRqvL78jcuXOHXgtcNEZWAKC29fQkHR1Je3vS0lLtNMAujawcddRRYxYMAGA89fQkra1Jb29y881JV5fCAkU3qdoBAADGQ0dHpagkle0ws0+AAlBWAICa0N6evPYOoObmpK2tunmA4VVtnRUAgPHU0lKZ+tXZWSkqpoBB8SkrAEDNaGlRUmAi2amy0tfXl+eeey6vvPJKmpqa8ra3vS2TJk28mWReXQwAANW1y68uTpJnn3023/nOd7Jq1aqsW7cuW7b8/9P33HPPvPOd78z73ve+nHjiiYUvLl5dDAAAxbFLry5eunRpVq1alWOPPTYf+tCHcthhh2W//fbLHnvskf7+/vz85z/PY489lv/xP/5Hvva1r+UTn/hEDj/88DH7wwAAALVjh2Vl//33z+23377dlrPvvvtm3333zVFHHZVzzjkn//Iv/5Jf//rXYxYUAACoLTssK+ecc86IL/QHf/AHuxwGAADgNSN+yOT73//+dvcPDg7m61//+qgFAgAASHairFx//fVZvHhx+vr6hvb9/Oc/z8UXX5wVK1aMSTgAAKB2jbis3HLLLenv7895552XBx98MHfccUcuuOCCHHbYYfnqV786lhlHXalUSrlcrnYMAACoWeVyOaVSaYfnDPvq4tf70pe+lLvuuiv19fW56qqrcsIJJ+xSyPHk1cUAAFAcw726eKcWRrnnnnvS0dGR9773vZk1a1a+8pWv5Kmnnhq1sAAAAK8Z8Qr2l156aX7xi1/kyiuvzAknnJBNmzbllltuycc+9rF88IMfzIc//OGxzAkAANSYEY+sNDU15ZZbbhma9jV16tT8+Z//ea677rqsXLlyzAICAAC1acQjK1deeeV29//BH/zBhHvAHgAAKL4Rl5XXbNmyJS+++GJeffXVrfYfcMABoxYKAABgxGWlr68vN9xwQ37wgx9k8+bN2xz/3ve+N6rBAACA2jbiZ1a++MUv5oUXXsiyZcsyderUfPazn80VV1yRgw46KIsWLRrDiAAAQC0acVl55JFH8hd/8ReZO3duJk2alJkzZ+bUU0/NxRdfnNtvv30sM446i0ICALuqpydZvLiyBXbeSBaFHHFZefXVV7PPPvskSaZPn55f/epXSZJZs2bl6aef3oWY42/JkiVZsGBBtWMAABNUT0/S2posWlTZKiyw8xYsWJAlS5bs8JwRP7Ny6KGH5sknn8zMmTNz+OGH5+/+7u/ym9/8Jvfff38OOuigXQ4LADBRdHQkvb2Vz729SWdn0tJS3UywOxrxyMpHPvKRNDQ0JEk++tGPpr6+PldffXX+7d/+LZdddtmYBQQAKJr29qS5ufK5uTlpa6tuHthdjXhk5fDDDx/6vP/+++cLX/jCmAQCACi6lpakq6syotLWZlQFxspOr7Py0ksvpa+vL4ODg1vtb37t1wsAADWgpUVJgbE24rKydu3afP7zn89TTz2VpLI4ZF1d3dDWOisAAMBoGnFZ+exnP5vm5uZccsklmTFjxlhmAgAAGHlZ6e3tzeLFizNr1qyxzAMAAJBkJ94G9q53vSv/+q//OpZZAAAAhox4ZOWyyy7L1Vdfnf/zf/5P3v72t2fy5K3/0tNPP33Uw42VUqmUtrY2C0MCAECVlMvldHZ27vCcEZeVH/zgB/nJT36SyZMnp7GxMXV1dUPH6urqJlRZWbJkSRobG6sdAwAAataCBQtywgkn5MYbb3zDc0ZcVm6++eZ8+MMfzllnnZVJk0Y8ewwAAOBNGXHr2Lx5c9797ncrKgAAwLgYcfNobW3NihUrxjILAADAkBFPA3v++efz8MMP56GHHsrBBx+c+vr6rY6XSqVRDwcAANSuEZeVyZMn56STThrLLAAAAENGXFauuOKKscwBAACwlWGfWfnIRz6Sb37zm3nxxRfHIw8AAECSEZSVY489NnfddVc+8IEP5JOf/GRWrVqVTZs2jUc2AACghg07Deyiiy7KhRdemH/6p3/K/fffn2uvvTb19fV597vfndNOOy2HH374eOQEAABqzIieWamrq8uxxx6bY489Ni+99FIeeOCBrFy5Mpdcckmam5tz6qmn5k//9E/HOisAwG6npyfp6Eja25OWlmqngWLZ6RUep02bltNOOy3XXXddrrzyyvT19eXWW28dg2hjp1QqpVwuVzsGAFDjenqS1tZk0aLKtqen2olg/JTL5WGXPxnx28Be8+yzz+b+++/P/fffn3Xr1uWd73xnFixY8KZDVsOSJUvS2NhY7RgAQI3r6Eh6eyufe3uTzk6jK9SOBQsW5IQTTsiNN974hueMqKwMDAxk1apV+Yd/+Ic89thjmTVrVk499dTMnz8/BxxwwKgFBgCoJe3tyc03V4pKc3PS1lbtRFAsw5aVxYsX50c/+lGmTJmSefPm5YILLvBQPQDAKGhpSbq6KiMqbW1GVeD1hi0rL7/8cq644or84R/+YaZMmTIemQAAakZLi5ICb2TYsvLZz352PHIAAABsZaffBvbCCy/kmmuuSe9rT4MBAACMgZ0uKytWrEi5XE5XV9dY5AEAAEjyJsvKqaeemnK5nM2bN49FJgAAgJ0rKz09Pdm4cWM+/vGPp76+Pj/60Y/GKhcAAFDjdqqsfOc738kpp5ySPffcM/Pnz893vvOdscoFAADUuBGXlY0bN+YHP/hBTjvttCSVFSf/9//+33nxxRfHLBwAAFC7ZyQgEwAAIABJREFURlxW7r///jQ3N+ewww5LkvzO7/xODjvssJTL5TELBwAA1K4Rl5UVK1bkfe9731b7TjvttKxYsWLUQwEAAIyorPz7v/979t5778yfP3+r/fPmzcv++++fdevWjUk4AACgdg27gn2S7L///rn++uu32b/XXnvlC1/4wqiHAgAA2Ol1VnYHpVLJszYAAFBF5XI5pVJph+fUdXd3b3mjg9ddd13OPvvszJw5c4cX2bx5c77//e8nyTZTxYpkYGAgra2t6evrS2NjY7XjAABATevv709TU1O6urrS0NCwzfEdTgM76KCDcvHFF+eggw7Kcccdl0MPPTT77bdfpkyZko0bN+YXv/hFHnvssfzwhz/MnDlz8hd/8Rdj9gcBAABqyw7Lygc+8IEsXLgw3d3deeCBB3LHHXfkN7/5TZKkrq4uBx10UI455phce+21OfTQQ8clMAAAUBuGfcB+zz33zPve976h1xZv3LgxmzZtyj777JMpU6aMeUAAAKA2jehtYL9t7733HoscAAAAW6nJt4EBAOxuenqSxYsrW9hdKCsAABNcT0/S2posWlTZKizsLpQVAIAJrqMj6e2tfO7tTTo7q5sHRouyAgAwwbW3J83Nlc/NzUlbW3XzwGjZqQfs+/r68vjjj2fDhg0ZHBzc6tjpp58+qsEAABiZlpakq6syotLWVvkOu4MRl5Xvf//7ufbaa1NXV5empqbU1dUNHaurq1NWAACqqKVFSWH3M+Ky8pWvfCVnnXVWzjnnnNTX149lJgAAgJE/s9Lf35/58+crKgAAwLgYcVl5z3vek4cffngsswAAAAwZ8TSwhoaG3Hrrrfmnf/qnvOMd79hmhOX8888f9XAAAEDtGnFZ+dd//dfMnj07L7/8clavXr3Vsd9+2B4AAGA0jLisLF26dCxzAAAAbMWikAAAQCEpKwAAQCEpKwAAQCEpKwAAQCEpKwAAQCHVZFkplUopl8vVjgEAADWrXC6nVCrt8Jy67u7uLeOUp+oGBgbS2tqavr6+NDY2VjsOAADUtP7+/jQ1NaWrqysNDQ3bHK/JkRUAAKD4lBUAAKCQlBUAAKCQlBUAAKCQlBUAAKCQlBUAAIb09CSLF1e2UG3KCgAASSoFpbU1WbSoslVYqDZlBQCAJElHR9LbW/nc25t0dlY3DygrAAAkSdrbk+bmyufm5qStrbp5YHK1AwAAUAwtLUlXV2VEpa2t8h2qSVkBAGBIS4uSQnGYBgYAABSSsgIAABSSsgIAABSSsgIAABSSsgIAABSSsgIAABSSsgIAABSSsgIAABSSsgIAABSSsgIAABSSsgIAABSSsgIAABSSsgIAABSSsgIAABSSsgIAABSSsgIAABSSsgIAABSSsgIAABTS5GoHGA0DAwO54oorUl9fn5dffjnvf//7M3/+/GrHAgAAdsFuUVamTZuWG264IfX19dm4cWPOPfdcZQUAACa4qpWVBx98MPfee2+efPLJDAwMZOXKlamvr9/qnNtvvz133313Nm7cmGOOOSaXXnppZsyYsc21Jk36/7PZXnrppRx88MFjnh8AgDfW05N0dCTt7UlLS7XTMFFV7ZmVV155JUcffXTOOuus7R5fsWJFbrvttlxyySVZvnx5BgYGsnjx4je8Xn9/f/7yL/8yf/Znf5aTTjpprGIDADCMnp6ktTVZtKiy7empdiImqqqNrLw2TevRRx/d7vF77rknZ555Zk4++eQkyeWXX56zzz47a9asyezZs3Pfffelu7s7jY2Nueqqq9LY2Jgbbrgh69evz4UXXph58+Zln332Gbc/DwAAFR0dSW9v5XNvb9LZaXSFN6eQbwPbtGlT1q5dm5bf+n91c3NzDjzwwKxevTpJsnDhwixbtixXXXVVNm3aNHTetGnTMnXq1EyZMmXccwMAUJn61dxc+dzcnLS1VTcPE1chH7Dv7+/P4OBg9t133632T58+PRs2bNjm/Keffjo33XRT6urq8uqrr+b888/PnnvuOV5xAQD4LS0tSVdXZUSlrc2oCm9eIcvKli1bdur8ww47LDfccMOIzy+VSpk6dWqSZMGCBVmwYMFO/TwAAHaspUVJYfvK5XLK5XKSbDVDansKWVaampoyadKkrF+/fqv9GzZsyPTp03f5+kuWLEljY+MuXwcAANg5vz1Y0N/fnxtvvPENzy3kMytTp07NIYccstXD988880yeffbZzJ07t4rJAACA8VK1kZX+/v48//zzWbduXZJkzZo1qa+vz6xZszJt2rScccYZWb58eebMmZOZM2fmpptuypFHHpnZs2dXKzIAADCOqlZWHnrooVxzzTVD3y+66KIkydKlS3PUUUfl9NNPz/r167Ns2bKhRSEvu+yyasUFAADGWV13d/fOPc0+gQ0MDKS1tTV9fX2eWQEAgCrr7+9PU1NTurq60tDQsM3xQj6zMtZKpdLQGwgAAIDxVy6XUyqVdniOkRUAAKAqjKwAAAATkrICAAAUkrICAAAUkrICAAAUkrICAAAUUk2WFa8uBgCA6vLq4tfx6mIAACgOry4GAAAmJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAopJosK15dDAAA1eXVxa/j1cUAABNPT0/S0ZG0tyctLdVOw2jy6mIAACasnp6ktTVZtKiy7empdiLGk7ICAEBhdXQkvb2Vz729SWdndfMwvpQVAAAKq709aW6ufG5uTtraqpuH8TW52gEAAOCNtLQkXV2VEZW2Ns+s1BplBQCAQmtpUVJqlWlgAABAISkrAABAIdVkWbEoJAAAVJdFIV/HopAAAFAcFoUEAAAmJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAopJosK1awBwCA6rKC/etYwR4AAIrDCvYAAMCEpKwAAACFpKwAAACFpKwAAACFpKwAAACFpKwAAACFpKwAAACFpKwAAACFpKwAAACFpKwAAACFVJNlpVQqpVwuVzsGAADjqKcnWby4sqX6yuVySqXSDs+p6+7u3jJOeapuYGAgra2t6evrS2NjY7XjAAAwTnp6ktbWpLc3aW5OurqSlpZqp6K/vz9NTU3p6upKQ0PDNsdrcmQFAIDa0tFRKSpJZdvZWd08jIyyAgDAbq+9vTKiklS2bW3VzcPITK52AAAAGGstLZWpX52dlaJiCtjEoKwAAFATWlqUlInGNDAAAKCQlBUAAKCQlBUAAKCQlBUAAKCQlBUAAKCQlBUAAKCQlBUAAKCQlBUAAKCQlBUAAKCQlBUAAKCQarKslEqllMvlascAAICaVS6XUyqVdnhOXXd395ZxylN1AwMDaW1tTV9fXxobG6sdBwAAalp/f3+amprS1dWVhoaGbY7X5MgKAABQfMoKAABQSMoKAABQSMoKAABQSMoKAABQSMoKAABQSMoKAABQSMoKAABQSMoKAABQSMoKAABQSMoKAABQSMoKAABQSMoKAABQSMoKAABQSMoKAADshJ6eZPHiypaxpawAAMAI9fQkra3JokWVrcIytpQVAAAYoY6OpLe38rm3N+nsrG6e3Z2yAgAAI9TenjQ3Vz43NydtbdXNs7ubXO0AAAAwUbS0JF1dlRGVtrbKd8aOsgIAADuhpUVJGS+mgQEAAIVUk2WlVCqlXC5XOwYAANSscrmcUqm0w3Pquru7t4xTnqobGBhIa2tr+vr60tjYWO04AABQ0/r7+9PU1JSurq40NDRsc7wmR1YAAIDiU1YAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAIBCUlYAAKAKenqSxYsrW7ZPWQEAgHHW05O0tiaLFlW2Csv2KSsAADDOOjqS3t7K597epLOzunmKarcqK319fWlra8vf//3fVzsKAAC8ofb2pLm58rm5OWlrq26eoppc7QCj6e/+7u9y5JFHVjsGAADsUEtL0tVVGVFpa6t8Z1tVKysPPvhg7r333jz55JMZGBjIypUrU19fv9U5t99+e+6+++5s3LgxxxxzTC699NLMmDFju9f72c9+lpdeeimHHnroeMQHAIBd0tKipAynatPAXnnllRx99NE566yztnt8xYoVue2223LJJZdk+fLlGRgYyOLFi9/werfcckvOPffcMUoLAACMt6qNrMyfPz9J8uijj273+D333JMzzzwzJ598cpLk8ssvz9lnn501a9Zk9uzZue+++9Ld3Z3Gxsa0trbmbW97Ww444IBxyw8AAIytQj6zsmnTpqxduzYXXnjh0L7m5uYceOCBWb16dWbPnp2FCxdm4cKFSZKvf/3refzxx3P55Zdn3bp1mTp1ambOnJmjjz66Wn8EAABgFxWyrPT392dwcDD77rvvVvunT5+eDRs2bHP+Oeeck3POOSdJcuutt2b//fdXVAAAYIIrZFnZsmXLm/5rR/LcSqlUytSpU5MkCxYsyIIFC970zwMAAEauXC6nXC4nqcyo2pFClpWmpqZMmjQp69ev32r/hg0bMn369F2+/pIlS9LY2LjL1wEAAHbObw8W9Pf358Ybb3zDcwu5KOTUqVNzyCGHbPXw/TPPPJNnn302c+fOrWIyAABgvFRtZKW/vz/PP/981q1blyRZs2ZN6uvrM2vWrEybNi1nnHFGli9fnjlz5mTmzJm56aabcuSRR2b27NnVigwAAIyjqpWVhx56KNdcc83Q94suuihJsnTp0hx11FE5/fTTs379+ixbtmxoUcjLLrusWnEBAIBxVtfd3f3mn2afYAYGBtLa2pq+vj7PrAAAQJX19/enqakpXV1daWho2OZ4IZ9ZGWulUmnoDQQAAMD4K5fLKZVKOzzHyAoAAFAVRlYAAIAJSVkBAAAKSVkBAAAKSVkBAAAKSVkBAAAKqSbLilcXAwBAdXl18et4dTEAABSHVxcDAAATkrICAAAUkrICAAAUkrICAAAUkrICAAAUUk2WFa8uBgCA6vLq4tfx6mIAACgOry4GAAAmJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAopJosKxaFBABgd9HTkyxeXNlOJBaFfB2LQgIAsDvp6UlaW5Pe3qS5OenqSlpaqp1q5CwKCQAAu6mOjkpRSSrbzs7q5hltygoAAExQ7e2VEZWksm1rq26e0Ta52gEAAIA3p6WlMvWrs7NSVCbSFLCRUFYAAGACa2nZ/UrKa0wDAwAACklZAQAACklZAQAACklZAQAACqkmy4oV7AEAoLqsYP86VrAHAIDisII9AAAwISkrAABAISkrAABAISkrAABAISkrAABAISkrAABAISkrAABAISkrAABAISkrAABAISkrAABAIdVkWSmVSimXy9WOAQAANatcLqdUKu3wnLru7u4t45Sn6gYGBtLa2pq+vr40NjZWOw4AANS0/v7+NDU1paurKw0NDdscr8mRFQAAoPiUFQAAoJCUFQAAoJCUFQAAoJCUFQAAoJCUFQAAoJCUFQAAoJCUFQAAoJCUFQAAoJCUFQAAoJCUFQAAoJCUFQAAoJCUFQAAoJCUFQAAoJBqsqyUSqWUy+VqxwAAgJpVLpdTKpV2eE5dd3f3lnHKU3UDAwNpbW1NX19fGhsbqx0HAABqWn9/f5qamtLV1ZWGhoZtjtfkyAoAAFB8ygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBIygoAAFBINVlWSqVSyuVytWMAAEDNKpfLKZVKOzynrru7e8s45am6gYGBtLa2pq+vL42NjdWOAwAANa2/vz9NTU3p6upKQ0PDNsdrcmQFAAAoPmUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAoJGUFAAAopMnVDjBaTj311MydOzdJ8p/+03/KOeecU+VEAADArthtykpTU1OWLVtW7RgAAMAoqVpZefDBB3PvvffmySefzMDAQFauXJn6+vqtzrn99ttz9913Z+PGjTnmmGNy6aWXZsaMGdu93q9//et8/OMfz5577pkLLrgg73jHO8bjjwEAAIyRqj2z8sorr+Too4/OWWedtd3jK1asyG233ZZLLrkky5cvz8DAQBYvXvyG17vzzjuzbNmynH/++bnqqqvGKjYAo6xcLlc7AgC/pUj35aqVlfnz5+ecc87J4Ycfvt3j99xzT84888ycfPLJmT17di6//PL89Kc/zZo1a5Ik9913Xz7+8Y/nyiuvTJJMnz49STJnzpzsueee2bBhw/j8QQDYJUX6lyIAxbovF/JtYJs2bcratWvT0tIytK+5uTkHHnhgVq9enSRZuHBhli1blquuuiovvfRSNm/enCR5/vnn09/fn8bGxqpkZ3QV6R+W8TbR/uxFyjveWcb6543F9UfrmkX6+87Yq+W/3xPxz16kzO7L43fNIv19Hw2FfMC+v78/g4OD2XfffbfaP3369O2OmPzsZz/L9ddfn7322iuvvvpqLr300kyatG0P27Jly9D1mRg6OztzwgknVDtGVUy0P3uR8o53lrH+eWNx/dG65mhcZ9OmTe7LE0SR/jkfbxPxz16kzO7L43fNiXZffu3nvPbf6a9X193dvf0j4+TRRx/NX/3VX231gP2///u/50/+5E9yyy23bPWg/J//+Z/nhBNOyIc+9KE39bNeuy4AAFAc3/zmN7P//vtvs7+QIytNTU2ZNGlS1q9fv9X+DRs2DD2b8mbst99++eY3v5lp06alrq5uV2MCAAC7YMuWLXnppZey3377bfd4IcvK1KlTc8ghh+TRRx/NMccckyR55pln8uyzzw4t/PhmTJo0abuNDQAAqI699977DY9Vraz09/fn+eefz7p165Ika9asSX19fWbNmpVp06bljDPOyPLlyzNnzpzMnDkzN910U4488sjMnj27WpEBAIBxVLVnVr773e/mmmuu2Wb/0qVLc9RRRyVJvvGNb2y1KORll132hotCAgAAu5eqP2APAACwPYV8ZqUaBgcH8/GPfzw///nP097eno985CPVjgRQ05544on8zd/8Terr67N58+ZcfPHFu/TcIgBv3sDAQK644orU19fn5Zdfzvvf//7Mnz9/zH+ukZXf8vzzz+eRRx7JunXrlBWAKnvhhRey5557pqGhIU8//XSuvfbafOlLX6p2LICaNDg4mC1btqS+vj4bN27Mueeem29961tj/nN3m5GVBx98MPfee2+efPLJDAwMbLVuy2tuv/32rZ6BufTSS7d6Buatb33reMcG2G3t6n35t19jOXXq1O0u9gvAyOzqPfm378EvvfRSDj744HHJvdvc+V955ZUcffTROeuss7Z7fMWKFbnttttyySWXZPny5RkYGMjixYvHOSVA7Rit+/J//Md/5Prrr3/TCwIDMDr35P7+/vzlX/5l/uzP/iwnnXTSeMTefUZWXpsz9+ijj273+D333JMzzzwzJ598cpLk8ssvz9lnn501a9Z4HTLAGBiN+/LmzZvzmc98JvPmzctxxx03PsEBdkOjcU9ubGzMDTfckPXr1+fCCy/MvHnzss8++4xp7t1mZGVHNm3alLVr16alpWVoX3Nzcw488MCsXr26iskAatNI7suDg4NZsmRJjjjiiLS2tlYrKsBubyT35E2bNg0dmzZtWqZOnZopU6aMebbdZmRlR/r7+zM4OJh99913q/3Tp0/Phg0bhr5fffXVeeqpp/LKK69k9erV+fznP2+ONMAYGMl9edWqVfnhD3+YF154IQ899FAaGhpy9dVXVyMuwG5tJPfkp59+OjfddFPq6ury6quv5vzzz8+ee+455tlqoqxs2TKyF559+tOfHuMkACQjuy+fcsopOeWUU8YhDUBtG8k9+bDDDssNN9wwDmm2VhPDBk1NTZk0aVLWr1+/1f4NGzZk+vTpVUoFULvclwGKo8j35JooK1OnTs0hhxyy1QNFzzzzTJ599lkLjAFUgfsyQHEU+Z6820wD6+/vz/PPP59169YlSdasWZP6+vrMmjUr06ZNyxlnnJHly5dnzpw5mTlzZm666aYceeSR3gQGMEbclwGKY6Lek3ebFey/+93v5pprrtlm/9KlS3PUUUclSb7xjW9stdDNZZddttWikACMHvdlgOKYqPfk3aasAAAAu5eaeGYFAACYeJQVAACgkJQVAACgkJQVAACgkJQVAACgkJQVAACgkJQVAACgkJQVAACgkJQVACak7373u5k3b17mzZuX//k//+fQ/nnz5uUnP/lJ1XItX758KFc1cwDsDpQVAIZ88pOfzKJFi7Z77Mtf/nLOO++88Q00jP322y/f/va3s3DhwjH/WS+//HJOO+20/OpXv9rheeeff36+/e1vj3kegFqgrAAwZMGCBXn44YezcePGrfYPDg7me9/7XhYsWFClZNtXV1eXGTNmZM899xzzn/XII4/kd37nd/KWt7xlh+fttddemTFjxpjnAagFk6sdAIDi+MM//MNMnTo1q1atSmtr69D+n/zkJ3nxxRfz3ve+N0nyuc99Lq+++moOOOCAdHZ2ZvLkyfnoRz+ak046Kdddd13+8R//Mc3NzfnUpz6V2bNnJ0kee+yx3HLLLXnq/7V3dyFN/n0cx9/TUWDFFBJMU5oERiwQppTpQVYODbIHiDTJwPBgHUSiQkVPlkYOH8DQlGFE4JGEdBbbQUWCq5OkpKCJFpimUm2m9iBj90E07t3z5v63v3LbfX9eR+73++26vteZn31/13V5vRiNRqxWK6dPn8ZkMgFQV1eH3++nubkZg8HA5OQkJ0+epKqqij179kR1PYFAgIaGBiYmJnA4HMzNzVFaWkpdXR09PT28ffsWi8XCpUuXeP78OU6nk5mZGQoLCzl16hQGgyF0rIGBAXbs2AHAmzdvuHnzJsPDwxiNRtLT02loaGDt2rVR1SkiIotTZ0VEREJWrVpFfn4+LpcrbNzlcpGVlRXWVRgYGACgvb2dw4cP09LSQn19Pbm5uTidTlJSUnA4HKH1X79+pbi4mK6uLhobG5mamqK1tTU0f+bMGd69e0dfXx/BYBCHw4HVao06qCwsLHDlyhWmp6dpampi3bp1obm7d+9it9tpb29ncnKSuro63G439fX1nD9/nvv374eu75enT5+Gwsr169exWCx0d3fT1tYWCnEiIrK0FFZERCSMzWZjaGiIiYkJ4GfI6O/vj9gClpiYSGVlJampqZSWlmI0GtmwYQMFBQVs3LiRkpISvF4v8/PzAGRnZ5Ofn09KSgoZGRnY7Xb6+/sJBAIAmEwmqqurcTqd3Lp1i5GREaqqqqK6hh8/fnDx4kXm5uZwOBysWbMmbL6srIzMzEw2b97Mvn37GBwcpKamBrPZTE5ODpmZmQwODobWDw8Ps7CwwJYtWwCYmpoiJyeH5ORkzGYz+/fvV1dFRGQZaBuYiIiEsVgspKSk4HK5OHHiBI8fPyY2Npbc3NywdWazOfR3bGwsJpOJTZs2hcYSEhIA8Pl8xMXFMT09jdPpZGhoCJ/PRzAYJBAI8OnTJxITEwHYuXMneXl59Pb2cvnyZeLj46O6hpaWFhITE2ltbWX16tUR8+np6WF1xsfHh+r9Nebz+UKfPR4P27dvJybm5298hw4dora2luzsbKxWK7t37w5tZxMRkaWjzoqIiESw2Wy43W4A3G43u3btivin32iM/L3rn8d+3e8RDAaBn/e5TE5OUl1dTWdnZ+ipY786K/CzI+L1eomJiWF8fDzq+rOyshgdHeX169eLzv9rnbGxsWHzBoMhVDeE368CUFlZSWdnJ1u3bsXtdlNeXs7Y2FjU9YqIyOIUVkREJEJBQQHj4+M8evSIwcHBJXkK2KtXrzhy5AhWq5W0tDT8fn/Emu7ubgCuXr3KnTt3GB0djepce/fuxW63c+HChX8bWP4qv9+P1+slOzs7bNxsNnPs2DE6OjpISEjgyZMnf+s8IiISSWFFREQiJCUlkZmZSXNzM0lJSWzbtu1vHzM5ORmXy8XY2BjPnj2jp6cnbP7ly5f09fVx7tw5cnNzKSws5MaNG2Gdl99RXFzM8ePHOXv2LCMjI1HX7fF4sFgsoftevn//TltbGy9evODDhw94PB6mpqZITU2N+hwiIrI4hRUREVmUzWZjdnYWm822JMerra3l/fv3VFRUcPv2bSoqKkJz3759o7GxkZKSEjIyMgCw2+18+fIlItT8jqNHj3Lw4EFqamqi3qbl8XjCtoDFxMTg8/m4du0a5eXltLW1UV5eTl5eXtR1iojI4gwPHz4M/udlIiIiK8uDBw/o7u6mt7d32c4RCAQ4cOAAHR0dpKWl/dZ38/PzaWpqwmq1LlN1IiL/+9RZERGRP9bHjx8pKiri3r17y3L8mZkZysrKfiuodHV1UVRUtCz1iIj8v1FnRURE/kjz8/N8/vwZ+PmOlpXynhO/38/s7CwA69evX/TRySIi8tfoPSsiIvJHiouLIy4u7r9dRgSTyaR3roiILBFtAxMRERERkRVJYUVERERERFYkhRUREREREVmRFFZERERERGRFUlj7ClNzAAAAIElEQVQREREREZEVSWFFRERERERWJIUVERERERFZkf4Bmh/TK63PSYIAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 960x960 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "vfbins, yvm_hdf5_distinct, yvm_hdf5_subhalos = np.loadtxt('vfunc_uchuu_results_hdf5_all.out')\n",
    "\n",
    "norm = 1/(np.log10(vfbins[1:len(vfbins)]) - [np.log10(vfbins[0:(len(vfbins)-1)])])\n",
    "\n",
    "plt.style.use('classic')\n",
    "plt.figure(figsize=(12, 12))\n",
    "\n",
    "plt.plot(vfbins[0:59], norm[0,0]*yvm_hdf5_distinct[0:59]/(2000**3), 'b.', label = \"Big Uchuu (hdf5) Distinct\")\n",
    "plt.plot(vfbins[0:59], norm[0,0]*yvm_hdf5_subhalos[0:59]/(2000**3), 'b+', label = \"Big Uchuu (hdf5) Subhalos\")\n",
    "plt.title('Big Uchuu (ascii) VS (hdf5) [z = 0.86]: velocity function', fontsize = 12)\n",
    "# plt.vlines(bins[:-1], 10, 40, linestyles = 'dashed', lw=.5)\n",
    "plt.legend(loc = 1, borderaxespad = 0.)\n",
    "plt.xlabel('Vmax [km/s]')\n",
    "plt.ylabel('n (>Vmax) [h^3/Mpc^3]')\n",
    "plt.xscale('log')\n",
    "plt.yscale('log')\n",
    "plt.xlim([10, 1200])\n",
    "# plt.show()\n",
    "plt.savefig(\"vfunc_ascii_vs_hdf5_biguchuu_z086.png\", bbox_inches='tight')\n"
   ]
  }
 ],
 "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"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}