{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import matplotlib.pyplot as plt\n", "import numpy\n", "import scienceplots\n", "from h5py import File\n", "\n", "import plt_utils\n", "\n", "\n", "%load_ext autoreload\n", "%autoreload 2" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [], "source": [ "with File(\"/mnt/extraspace/rstiskalek/csiborg_postprocessing/ACL/BORG2_0.25.hdf5\", 'r') as f:\n", " voxel_acl = f['voxel_acl'][...].flatten()\n", " voxel_dist = f['voxel_dist'][...].flatten()" ] }, { "cell_type": "code", "execution_count": 28, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnUAAAHWCAYAAAARl3+JAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAACNbUlEQVR4nOzdeVxU9foH8M+wzDCALILAgKgIqKEmpYGgRd4ovGJCdr1olksuZZgQLomJplmkZuJOWindNNBSKC26hqZ2RUwQf+4buDMsyiIIA8L5/UFMjswAg8Pq5/16zQvne55zzjNnLB/O+S4iQRAEEBEREVGbptfSCRARERHRo2NRR0RERNQOsKgjIiIiagdY1BERERG1AyzqiIiIiNoBFnVERERE7QCLOiIiIqJ2gEUdERERUTvAoo6IiIioHWBRR0RtRnZ2Nv71r3/BysoKIpEIUVFRGmNFIhE+/PBDnZ7/4WNu2bIFIpEIV65c0fpY3bp1w4QJE+qNe5RzNNSVK1cgEomwZcuWJjtHU2sPn4HoUbGoI3pMrV+/HiKRCJ6ennXGZWdnY9asWejVqxeMjY1hYmKC/v37Y8mSJSgoKFDGPf/88+jTp0+T5vzee+/h119/RXh4OP7zn/9g6NChTXo+IqK2xKClEyCilrF161Z069YNR48exaVLl+Di4lIr5s8//8SwYcNQXFyM119/Hf379wcAHDt2DJ9++ikOHjyI//73v82W8759+xAQEIBZs2bVG1taWgoDg6b9X9wbb7yB0aNHQyKRaL3v+fPnoafH36uJSHdY1BE9hjIzM3H48GHs3LkTb731FrZu3YqFCxeqxBQUFOCVV16Bvr4+jh8/jl69eqls//jjj7Fp06bmTBs5OTmwsLBoUKyRkVHTJgNAX18f+vr6jdq3MYUgEVFd+Gsi0WNo69atsLS0hL+/P/71r39h69attWK++OIL3Lx5E59//nmtgg4AbG1tMX/+fJ3kk5GRgVGjRqFjx44wNjbGwIEDsWfPHuX2mn5lgiBg3bp1EIlEEIlEdR7z4f5vH374IUQiES5duoQJEybAwsIC5ubmmDhxIu7du6eyr0KhwHvvvYdOnTqhQ4cOGDFiBG7cuFHrHA/3dxs+fDi6d++uNh8vLy8MGDBA+V5dn7rTp0/jH//4B6RSKTp37owlS5agqqqq3s+m6Zh37tzBrFmz0LdvX5iamsLMzAz//Oc/ceLECbU51uXYsWMQiUSIiYmpte3XX3+FSCTC7t27lW3Hjx/HP//5T5iZmcHU1BQvvPACjhw5oty+b98+6OnpYcGCBSrH2rZtG0QiETZs2KBsu3nzJt58803Y2tpCIpGgd+/e+Prrr+vNWS6XY+LEiejcuTMkEglkMhkCAgKatH8iUUvinTqix9DWrVsxcuRIiMVijBkzBhs2bMCff/6JZ555Rhnz448/QiqV4l//+leT5pKdnQ1vb2/cu3cPM2bMgJWVFWJiYjBixAh8//33eOWVV/Dcc8/hP//5D9544w28+OKLGDduXKPP9+9//xtOTk6IjIxEWloavvzyS9jY2GDp0qXKmMmTJ+Pbb7/Fa6+9Bm9vb+zbtw/+/v71HjsoKAjjxo2rdS2vXr2KI0eOYPny5Rr3lcvlGDJkCO7fv4+5c+fCxMQEGzduhFQqbfRnzcjIQHx8PEaNGgUnJydkZ2fjiy++gI+PD86cOQN7e/sGH2vAgAHo3r07tm/fjvHjx6tsi4uLg6WlJfz8/ABUF6fPPvsszMzMMGfOHBgaGuKLL77A888/jwMHDsDT0xP/+Mc/8M477yAyMhKBgYF4+umnkZWVhXfffRe+vr54++23AVT//Rg4cCBEIhGmT5+OTp064ZdffsGkSZNQVFSE0NBQjTm/+uqrOH36NN59911069YNOTk52Lt3L65du4Zu3bppfT2JWj2BiB4rx44dEwAIe/fuFQRBEKqqqoTOnTsLISEhKnGWlpZCv379GnxcHx8foXfv3lrnExoaKgAQDh06pGy7e/eu4OTkJHTr1k2orKxUtgMQgoODG3RcAMLChQuV7xcuXCgAEN58802VuFdeeUWwsrJSvk9PTxcACO+8845K3GuvvVbrmJs3bxYACJmZmYIgCEJhYaEgkUiEmTNnquy7bNkyQSQSCVevXlW2de3aVRg/fnyt65CSkqJsy8nJEczNzVXOoe6zaTpmWVmZyvUTBEHIzMwUJBKJsHjxYpU2AMLmzZtrHfNB4eHhgqGhoXDnzh1lm0KhECwsLFSua2BgoCAWi4XLly8r227duiV06NBBeO6555RtJSUlgouLi9C7d2+hrKxM8Pf3F8zMzFSu06RJkwSZTCbk5eWp5DJ69GjB3NxcuHfvntrPkJ+fLwAQli9fXudnImpP+PiV6DGzdetW2NraYsiQIQCqH+UFBQUhNjYWlZWVyriioiJ06NChyfP5+eef4eHhgcGDByvbTE1NMXXqVFy5cgVnzpzR6flq7gDVePbZZ3H79m0UFRUp8wGAGTNmqMTVdUeoRs3jze3bt0MQBGV7XFwcBg4ciC5dumjc9+eff8bAgQPh4eGhbOvUqRPGjh1b73k1kUgkysEYlZWVuH37NkxNTdGzZ0+kpaVpfbygoCBUVFRg586dyrb//ve/KCgoQFBQkPI8//3vfxEYGKjyKFomk+G1117DH3/8obzWxsbG2LJlC86ePYvnnnsOe/bswcqVK5XXSRAE/PDDD3j55ZchCALy8vKULz8/PxQWFmr8HFKpFGKxGL///jvy8/O1/qxEbRGLOqLHSGVlJWJjYzFkyBBkZmbi0qVLuHTpEjw9PZGdnY2kpCRlrJmZGe7evdvkOV29ehU9e/as1f7EE08ot+vSw4WVpaUlACj/4b969Sr09PTg7OysEqcuR3WCgoJw/fp1JCcnAwAuX76M1NRUZdGjydWrV+Hq6lqrvaHnVaeqqgorV66Eq6srJBIJrK2t0alTJ/zf//0fCgsLtT5ev3790KtXL8TFxSnb4uLiYG1tjX/84x8AgNzcXNy7d0/jd1pVVYXr168r2wYNGoRp06bh6NGj8PPzw5tvvqnclpubi4KCAmzcuBGdOnVSeU2cOBFA9eAZdSQSCZYuXYpffvkFtra2eO6557Bs2TLI5XKtPzdRW8Gijugxsm/fPmRlZSE2Nhaurq7K17///W8AUBkw0atXL1y4cAHl5eUtlW6T0DRa9cE7a4/i5ZdfhrGxMbZv3w4A2L59O/T09DBq1CidHL8uD95pBYBPPvkEYWFheO655/Dtt9/i119/xd69e9G7d2+1AzAaIigoCPv370deXh4UCgV+/PFHvPrqq42ePkahUOD3338HUF0APzhopSbH119/HXv37lX7GjRokMZjh4aG4sKFC4iMjISRkREiIiLwxBNP4Pjx443Klai140AJosfI1q1bYWNjg3Xr1tXatnPnTuzatQvR0dGQSqV4+eWXkZycjB9++AFjxoxpspy6du2K8+fP12o/d+6ccntz6tq1K6qqqnD58mWVu03qclTHxMQEw4cPx44dO/D5558jLi4Ozz77bL2DErp27YqLFy/Wald3XktLS5WJnwGgvLwcWVlZKm3ff/89hgwZgq+++kqlvaCgANbW1g36PA8LCgrCokWL8MMPP8DW1hZFRUUYPXq0cnunTp1gbGys8TvV09ODo6Ojsm3hwoU4e/YsPvvsM7z//vuYO3cuVq9erTxWhw4dUFlZCV9f30bl6+zsjJkzZ2LmzJm4ePEi3N3dsWLFCnz77beNOh5Ra8Y7dUSPidLSUuzcuRPDhw/Hv/71r1qv6dOn4+7du/jxxx8BVPc9k8lkmDlzJi5cuFDreDk5OViyZMkj5zVs2DAcPXpU+bgSAEpKSrBx40Z069YNbm5uj3wObfzzn/8EAGVhUaOuJckeFhQUhFu3buHLL7/EiRMn6n30ClRfhyNHjuDo0aPKttzcXLXTzTg7O+PgwYMqbRs3bqx1p05fX7/WHcgdO3bg5s2bDf4sD3viiSfQt29fxMXFIS4uDjKZDM8995zKOV966SUkJCSoTB2SnZ2Nbdu2YfDgwTAzMwMApKSk4LPPPkNoaChmzpyJ2bNnY+3atThw4IDyWK+++ip++OEHnDp1qlYuubm5GvO8d+8eysrKVNqcnZ3RoUMHKBSKRn9+otaMd+qIHhM//vgj7t69ixEjRqjdPnDgQHTq1Albt25FUFAQLC0tsWvXLgwbNgzu7u4qK0qkpaXhu+++g5eXl8oxcnNz1RZ6Tk5OGjv8z507F9999x3++c9/YsaMGejYsSNiYmKQmZmJH374odlXXXB3d8eYMWOwfv16FBYWwtvbG0lJSbh06VKDjzFs2DB06NABs2bNUhYm9ZkzZ45y6bOQkBDllCZdu3bF//3f/6nETp48GW+//TZeffVVvPjiizhx4gR+/fXXWnffhg8fjsWLF2PixInw9vbGyZMnsXXrVo1z6TVUUFAQFixYACMjI0yaNKnWd7RkyRLs3bsXgwcPxjvvvAMDAwN88cUXUCgUWLZsGQCgrKwM48ePh6urKz7++GMAwKJFi/DTTz9h4sSJOHnyJExMTPDpp59i//798PT0xJQpU+Dm5oY7d+4gLS0Nv/32G+7cuaM2xwsXLuCFF17Av//9b7i5ucHAwAC7du1Cdna2yp1FonalRcfeElGzefnllwUjIyOhpKREY8yECRMEQ0NDlekjbt26Jbz33ntCjx49BCMjI8HY2Fjo37+/8PHHHwuFhYXKOB8fHwGA2tcLL7xQZ26XL18W/vWvfwkWFhaCkZGR4OHhIezevbtWHHQwpUlubq5K3MPTkgiCIJSWlgozZswQrKysBBMTE+Hll18Wrl+/Xu+UJg8aO3asAEDw9fVVm9/D048IgiD83//9n+Dj4yMYGRkJDg4OwkcffSR89dVXtc5RWVkpvP/++4K1tbVgbGws+Pn5CZcuXVI7pcnMmTMFmUwmSKVSYdCgQUJycrLg4+Mj+Pj4KOMaOqVJjYsXLyq/2z/++ENtTFpamuDn5yeYmpoKxsbGwpAhQ4TDhw8rt7/33nuCvr6+yhQuglA95Y6BgYEwbdo0ZVt2drYQHBwsODo6CoaGhoKdnZ3wwgsvCBs3btT4GfLy8oTg4GChV69egomJiWBubi54enoK27dvb9BnJGqLRIKgo97BRERERNRi2KeOiIiIqB1gUUdERETUDrCoIyIiImoHWNQRERERtQMs6oiIiIjaARZ1RERERO0AJx9upKqqKty6dQsdOnSASCRq6XSIiIionRIEAXfv3oW9vX2dE7KzqGukW7duqaxfSERERNSUrl+/js6dO2vczqKukTp06ACg+gLXrGNIREREpGtFRUVwdHRU1h6asKhrpJpHrmZmZizqiIiIqMnV192LAyWIiIiI2gEWdURERETtAIs6IiIionaAfeqIiKjVq6qqQnl5eUunQdQkDA0Noa+v/8jHYVFHREStWnl5OTIzM1FVVdXSqRA1GQsLC9jZ2T3S3Lcs6oiIqNUSBAFZWVnQ19eHo6NjnROvErVFgiDg3r17yMnJAQDIZLJGH4tFHRERtVr379/HvXv3YG9vD2Nj45ZOh6hJSKVSAEBOTg5sbGwa/SiWv/IQEVGrVVlZCQAQi8UtnAlR06r5paWioqLRx2BRR0RErR7X2Kb2Thd/x1nUEREREbUDLOqIiIiI2gEWdURERETtAEe/EhFRm3OzoBT5Jc03GbGliRgOFtJmO19zuX37Np544gkcPXoU3bp1a+l0lGbNmoVLly4hPj5e7fbRo0fjmWeewcyZM5s3sVaORR0REbUpNwtK4bviAEorKpvtnFJDffw200frwm7dunVYvnw55HI5+vXrhzVr1sDDw0Onufn4+ODgwYMAqlcmcHZ2RkREBF577bVasRMnToSDgwOWLFkCAPj4448REBCgLOgmTJiAmJgYvPXWW4iOjlbZNzg4GOvXr8f48eOxZcsWnX6Gh/NKT0+Ht7e3xvj58+fjueeew+TJk2Fubq6zPCIjI7Fz506cO3cOUqkU3t7eWLp0KXr27KkS1xzfa2OwqCMiojYlv6QcpRWViApyh4uNaZOf71JOMULj0pFfUq5VURcXF4ewsDBER0fD09MTUVFR8PPzw/nz52FjY6NVDvn5+TA0NISpqernFQQBx48fx2effYaxY8eitLQUq1evxrhx4+Dl5QUnJydlbGVlJXbv3o09e/YAAO7du4evvvoKv/76q8oxHR0dERsbi5UrVyrnTysrK8O2bdvQpUsXrfJuiIfzAoATJ05g2rRpGvfp06cPnJ2d8e233yI4OFjrc2q6ngcOHEBwcDCeeeYZ3L9/H/PmzcNLL72EM2fOwMTEBIBuv1edE6hRCgsLBQBCYWFhS6dCj5lbd28Jp/NO1/u6dfdWS6dK9MhKS0uFM2fOCKWlpcq2kzcKhK7v7xZO3iholhwaez4PDw8hODhY+b6yslKwt7cXIiMjG7R/RUWFsHv3buFf//qXIJFIhPT09Fox58+fFwAIp06d+jvfkycFAMIvv/yiEnvw4EFBJpMJVVVVgiAIwo4dO4ROnTqpxIwfP14ICAgQ+vTpI3z77bfK9q1btwpPPvmkEBAQIIwfP17Z7uPjIwQHBwvBwcGCmZmZYGVlJcyfP195jprPvXTpUsHZ2VkQi8WCo6OjsGTJEo15Xb9+XQAgfP/994Kvr68glUqFHj16CEeOHFHJddGiRcLgwYMbdC0FoWHX82E5OTkCAOHAgQPKtkf9XjVR93e9RkNrDt6pI2pDsoqzEJAQgNL7pfXGSg2kSAhIgMy08UvOEFHjlJeXIzU1FeHh4co2PT09+Pr6Ijk5uc59T548iS1btmDr1q2oqKhAUFAQ9u/fj379+tWKTU1NhaWlJdzc3AAAN27cwAcffACJRIInn3xSJfbHH3/Eyy+/rJwP7dChQ+jfv7/aHN58801s3rwZY8eOBQB8/fXXmDhxIn7//fdasTExMZg0aRKOHj2KY8eOYerUqejSpQumTJkCAAgPD8emTZuwcuVKDB48GFlZWTh37pzGvNLT0wFUP+KMiIhA586d8c4772Du3LnYv3+/cj8PDw98/PHHUCgUkEgkOrmeDyssLAQAdOzYEcCjfa/NgUUdURuSr8hH6f1SRD4bie7m3TXGZRRmIPxQOPIV+SzqiFpAXl4eKisrYWtrq9Jua2urUtDUuH37Nr799lvExMTg9OnTGDZsGNavX4/hw4fXuZpGWloaCgsL0aFDB1RWVqKsrAxSqRTR0dGwt7dXiU1ISMDKlSuV769evVorpsbrr7+O8PBwXL16FQDwv//9D7GxsWqLOkdHR6xcuRIikQg9e/bEyZMnsXLlSkyZMgV3797FqlWrsHbtWowfPx4A4OzsjMGDB2vMKz09HR07dsT27dthbW0NABgxYgS++OILlfPa29ujvLwccrkcXbt21cn1fFBVVRVCQ0MxaNAg9OnTB4D232tzY1FH1AZ1N+8ONyu3lk6DiHRkzZo1WLRoEZ599llcunQJjo6ODdovLS0NwcHBmDFjBgoKCjBr1iwMGjQIEyZMUIk7e/Ysbt26hRdeeEHZVlpaCiMjI7XH7dSpE/z9/bFlyxYIggB/f39lgfWwgQMHqqyG4OXlhRUrVqCyshJnz56FQqFQOW99eaWnpyMgIEDlfJmZmXBxcVHZt6a/371792odt7HX80HBwcE4deoU/vjjD633bSmcp46IiEjHrK2toa+vj+zsbJX27Oxs2NnZ1YqfOnUqPvroI8jlcvTu3RsTJ07Evn37UFVVVed50tLS4O3tDRcXFwwYMADr16/H0qVLceXKFZW4H3/8ES+++KJKEWdtbY38/HyNx37zzTexZcsWxMTE4M0332zAp66tpvDSRF1e6enpGDhwoEpceno63N3dVdru3LkDoLoAfVhjr2eN6dOnY/fu3di/fz86d+6sbNf2e21uLV7UrVu3Dt26dYORkRE8PT1x9OjROuN37NiBXr16wcjICH379sXPP/+ssl0QBCxYsAAymQxSqRS+vr64ePGiSsyFCxeUvwWYmZlh8ODBKs/piYiIHoVYLEb//v2RlJSkbKuqqkJSUhK8vLxqxdvb22P+/Pm4cOECEhMTIRaLMXLkSHTt2hVz587F6dOna+2TkZGBgoIC5aNBAHBzc4OzszO2bdumEpuQkICAgACVtqeeegpnzpzR+BmGDh2K8vJyVFRUwM/PT2NcSkqKyvsjR47A1dUV+vr6cHV1hVQqVbkOdeV19+5dZGRk4KmnnlKJU1fUnTp1Cp07d1Z7B7Ex1xOoriGmT5+OXbt2Yd++fSqjhwHtv9fm1qJFXc2w4IULFyItLQ39+vWDn58fcnJy1MYfPnwYY8aMwaRJk3D8+HEEBgYiMDAQp06dUsYsW7YMq1evRnR0NFJSUmBiYgI/Pz+UlZUpY4YPH4779+9j3759SE1NRb9+/TB8+HDI5fIm/8xERPR4CAsLw6ZNmxATE4OzZ89i2rRpKCkpwcSJE+vcz9vbG1988QXkcjmWL1+O9PR09OvXDydPnlSJS01NhaGhIXr06KHS/sILL2DXrl3K9zk5OTh27BiGDx+uEufn54fTp09rvFunr6+Ps2fP4syZM9DX19eY77Vr1xAWFobz58/ju+++w5o1axASEgIAMDIywvvvv485c+bgm2++weXLl3HkyBF89dVXavM6ceIE9PX10bdvX2Xb1atXkZ+fX6uoO3ToEF566SWNedVo6PUEqh+5fvvtt9i2bRs6dOgAuVwOuVyO0tK/B6c19nttFo80/vYRaTss+N///rfg7++v0ubp6Sm89dZbgiAIQlVVlWBnZycsX75cub2goECQSCTCd999JwiCIOTm5goAhIMHDypjioqKBADC3r17G5w7pzShlnA677TQZ0sf4XTeaZ3EEbV2dU1psivthnDyRkGTv3al3Wj0FCpr1qwRunTpIojFYsHDw6PWtBwNdfPmzVr/3sydO1dwc3OrFfvDDz8IIpFIuH79uiAIgvDll18KgwYNUntcDw8PITo6Wvm+ZkoTTdRNafLOO+8Ib7/9tmBmZiZYWloK8+bNqzWlyZIlS4SuXbsKhoaGQpcuXYRPPvlEbV5r1qwRevfurdK2a9cuwcLCQqWttLRUMDc3F5KTkzXmWhd111MQBAGA2tfmzZtr5amL7/VBupjSRPTXh2h25eXlMDY2xvfff4/AwEBl+/jx41FQUICEhIRa+3Tp0gVhYWEIDQ1Vti1cuBDx8fE4ceIEMjIy4OzsjOPHj6tU9D4+PnB3d8eqVasgCAKeeOIJPPvss4iKioJEIkFUVBSWL1+Oc+fOwdLSskH5FxUVwdzcHIWFhTAzM2vsZSDSypnbZxC0Owhxw+PqHCjR0Dii1q6srAyZmZlwcnJS9rtqSytKtAYjRozA4MGDMWfOnFrb9uzZg9mzZ+PUqVPQ09P+4d3zzz8Pd3d3REVF6TSv+mzYsAG7du3Cf//7X633ba3U/V2v0dCao8VGvzZmWLBcLlcbX/PYtOZnXTEikQi//fYbAgMD0aFDB+jp6cHGxgaJiYl1FnQKhQIKhUL5vqioqIGflIiIdMnBQorfZvpw7dcGGjx4MMaMGaN2m7+/Py5evIibN282aoRoU+VVH0NDQ6xZs0bHGbV9j92UJoIgIDg4GDY2Njh06BCkUim+/PJLvPzyy/jzzz8hk6mf0ysyMhKLFi1q5myJiEgdBwtpmy2ymlt9d8IefPrVnBpzh67G5MmTdZhJ+9FiAyUaMyzYzs6uzvian3XF7Nu3D7t370ZsbCwGDRqEp59+GuvXr4dUKkVMTIzGfMPDw1FYWKh8Xb9+XbsPTERE1M78/vvvjXr0Sk2jxYq6xgwL9vLyqjUseu/evcp4Jycn2NnZqcQUFRUhJSVFGVMzSeHDfQf09PTqnL9GIpHAzMxM5UVERETUWrTo49ewsDCMHz8eAwYMgIeHB6KiolSGBY8bNw4ODg6IjIwEAISEhMDHxwcrVqyAv78/YmNjcezYMWzcuBFAdX+50NBQLFmyBK6urnByckJERATs7e2VgzG8vLxgaWmJ8ePHY8GCBZBKpdi0aRMyMzPh7+/fIteBHnMF14F7txsWW3Sl+mfuBUC/A2DRvH1giIio9WrRoi4oKAi5ublYsGAB5HI53N3dkZiYqBzocO3aNZU7at7e3ti2bRvmz5+PefPmwdXVFfHx8SoTL86ZMwclJSWYOnUqCgoKMHjwYCQmJipHklhbWyMxMREffPAB/vGPf6CiogK9e/dGQkJCgxb3JdKpguvAOg+govYyN2qJDQEHGbBzMiAYAsFHWdgREREAoMWmNGnrOKUJ6cStdGCjDzByE2Ddo97wM0VXEJQSgbiek+GWuACYegCwd68dxylNqJ2oa5oHovakTU9pQkQPsO6htjirRSKu/mnRtUnTISKitqfF134lIiIiokfHoo6IiIioHWBRR0RERNQOsKgjIiIiagdY1BERERG1Axz9SkREbY82k3brgrHVYzsn5O3bt/HEE0/g6NGj6NatW0unozRr1ixcunQJ8fHxarePHj0azzzzDGbOnNm8ibUgFnVERNS2aDtpty4YGms92XdkZCR27tyJc+fOQSqVwtvbG0uXLkXPnj2VMevWrcPy5cshl8vRr18/rFmzBh4eHjpP38fHBwcPHqz+KIaGcHZ2RkREBF577bVasRMnToSDgwOWLFkCAPj4448REBCgLOgmTJiAmJgYvPXWW4iOjlbZNzg4GOvXr8f48eOxZcsWnX6Gh/NKT0+Ht7e3xvj58+fjueeew+TJk2Fubq7TXGp8+umnCA8PR0hISK01cJvru30QizoiImpb7t2uLugaOGn3I8u7AOycUn1eLYq6AwcOIDg4GM888wzu37+PefPm4aWXXsKZM2dgYmKCuLg4hIWFITo6Gp6enoiKioKfnx/Onz8PGxsbrdPMz8+HoaEhTE1NVdoFQcDx48fx2WefYezYsSgtLcXq1asxbtw4eHl5wcnJSRlbWVmJ3bt3Y8+ePQCq10v/6quv8Ouvv6oc09HREbGxsVi5ciWkUimA6slzt23bhi5dumide30ezgsATpw4gWnTpmncp0+fPnB2dsa3336L4OBgrc+p6XrW+PPPP/HFF1/gySefrLVN199tQ7FPHRERtU01k3Y39auRhWNiYiImTJiA3r17o1+/ftiyZQuuXbuG1NRUAMDnn3+OKVOmYOLEiXBzc0N0dDSMjY3x9ddfN/gc9+/fx549ezBq1CjIZDJcvny5VszFixdx9+5dDB06FHZ2dnBycsKkSZNQWVmJ8+fPq8QePnwYhoaGeOaZZwAAP//8MyQSCQYOHKgS9/TTT8PR0RE7d+5Utu3cuRNdunTBU089pRL7/PPPY/r06Zg+fTrMzc1hbW2NiIgIPLigVVVVFZYtWwYXFxdIJBJ06dIFH3/8sca8bty4gby8PADAiy++CGNjY/Ts2RMpKSkq53755ZcRGxur0+sJAMXFxRg7diw2bdoES0vLWtt18d02Bu/UEVGdbhaUIr+kvEnPYWkihoOFtEnPQdTSCgsLAQAdO3ZEeXk5UlNTER4ertyup6cHX19fJCcn13uskydPYsuWLdi6dSsqKioQFBSE/fv3q13DPDU1FZaWlnBzq14y8MaNG/jggw8gkUhq3WX68ccf8fLLL0MkEgEADh06hP79+6vN4c0338TmzZsxduxYAMDXX3+NiRMn4vfff68VGxMTg0mTJuHo0aM4duwYpk6dii5dumDKlCkAgPDwcGzatAkrV67E4MGDkZWVhXPnzmnMKz09HUD1I86IiAh07twZ77zzDubOnYv9+/cr9/Pw8MDHH38MhUIBiUSik+sJVD9m9vf3h6+vr/JxcI1H/W4fBYs6ItLoZkEpfFccQGlFZZOeR2qoj99m+rCwo3arqqoKoaGhGDRoEPr06YNbt26hsrIStra2KnG2trYqxcyDbt++jW+//RYxMTE4ffo0hg0bhvXr12P48OEQi8Uaz52WlobCwkJ06NABlZWVKCsrg1QqRXR0NOzt7VViExISsHLlSuX7q1ev1oqp8frrryM8PBxXr14FAPzvf/9DbGys2qLO0dERK1euhEgkQs+ePXHy5EmsXLkSU6ZMwd27d7Fq1SqsXbsW48ePBwA4Oztj8ODBGvNKT09Hx44dsX37dlhbWwMARowYgS+++ELlvPb29igvL4dcLkfXrqrLKzb2esbGxiItLQ1//vmn2u15eXlaf7e6wqKOiDTKLylHaUUlooLc4WKjvl/Jo7qUU4zQuHTkl5SzqKN2Kzg4GKdOncIff/zR6GOsWbMGixYtwrPPPotLly7B0bFh/fvS0tIQHByMGTNmoKCgALNmzcKgQYMwYcIElbizZ8/i1q1beOGFF5RtpaWltRaXr9GpUyf4+/tjy5YtEAQB/v7+ygLrYQMHDlTeZQMALy8vrFixApWVlTh79iwUCoXKeevLKz09HQEBASrny8zMhIuLi8q+Nf397t2rPaimMdfz+vXrCAkJwd69ezVel5bEoo6I6uViY4o+Dk0zeoyovZs+fTp2796NgwcPonPnzgAAa2tr6OvrIzs7WyU2OzsbdnZ2ao8zdepUGBgY4JtvvkHv3r3x6quv4o033sDzzz8PPT3NXeTT0tIwZcoUZcGzfv16PPnkk5gyZYrKFCU//vgjXnzxRZVixdraGvn5+RqP/eabb2L69OkAqh+FNkZN4aWJurzS09MxZ84clbj09HQ899xzKm137twBUF2APqwx1zM1NRU5OTl4+umnlW2VlZU4ePAg1q5dC4VC0ajvVlc4UIKIiKgJCIKA6dOnY9euXdi3b5/KKFOxWIz+/fsjKSlJ2VZVVYWkpCR4eXmpPZ69vT3mz5+PCxcuIDExEWKxGCNHjkTXrl0xd+5cnD59utY+GRkZKCgoQJ8+fZRtbm5ucHZ2xrZt21RiExISEBAQoNL21FNP4cyZMxo/49ChQ1FeXo6Kigr4+flpjHt4AMORI0fg6uoKfX19uLq6QiqVqlyLuvK6e/cuMjIyag3ISE9Ph7u7u0rbqVOn0LlzZ7V3EBtzPV944QWcPHkS6enpyteAAQMwduxYpKenQ19fv1Hfra6wqCMiImoCwcHB+Pbbb7Ft2zZ06NABcrkccrkcpaWlAICwsDBs2rQJMTExOHv2LKZNm4aSkhJMnDix3mN7e3vjiy++gFwux/Lly5Geno5+/frh5MmTKnGpqakwNDREjx6qI3hfeOEF7Nq1S/k+JycHx44dw/Dhw1Xi/Pz8cPr0aY136/T19XH27FmcOXMG+vr6GvO9du0awsLCcP78eXz33XdYs2YNQkJCAABGRkZ4//33MWfOHHzzzTe4fPkyjhw5gq+++kptXidOnIC+vj769u2rbLt69Sry8/NrFXWHDh3CSy+9pDGvGg29nh06dECfPn1UXiYmJrCyslIpnB/lu30UfPxKRERtU96FVn2eDRs2AKie0uNBmzdvxoQJExAUFITc3FwsWLAAcrkc7u7uSExMrNXBvi5GRkYYPXo0Ro8ejVu3btWaUy0tLQ2urq61Ov77+voiOjoaN27cQOfOnfHTTz/Bw8Oj1h2tvn374umnn8b27dvx1ltvqc3BzMys3jzHjRuH0tJSeHh4QF9fHyEhIZg6dapye0REBAwMDLBgwQLcunULMpkMb7/9ttq80tPT0bNnT5XHscePH4eFhYXK4+SysjLEx8cjMTGx3vxq1Hc9G0oX321jsKgjIqK2xdiqeoWHnVOa75yGxtXn1cKD87BpUjN/my6oG6UaGRmJyMjIWu0jR45EVVWV8n1CQgJGjBih9rgLFizA7NmzMWXKFOjp6dW7UoS6ZbsMDQ0RFRWlLHQfpqenhw8++AAffPCBSvuIESNq5aXumgUGBiIwMFClbfPmzfDw8Kg1x15DaRr1+zB1o30B3X63DcWijqiNUzePXGZRMQDgck4xqsoKG33sSznFj5QbUZOwcKxesotrv+rM4MGDMWbMGLXb/P39cfHiRdy8ebPBI26bI6/6GBoaYs2aNTrOqHVjUUfUhuUUK+C7ofY8cnpGN2HiBITEpaOqLPeRziE11IelieY5m4hahIVjuy6ymtvDI0kfFhoa2jyJPKS+vOoyefJkHWbSNrCoI2rDikor1M4jl1l0HvOOAquC3OFk1rOOI9SPqz0Q0aPQ9HiSdI9FHVE78PA8cnpG1QWes40p3Kw4vxwR0eOARR0R0UOyirOQr9A84WoNS4klZKayZsiIiKh+LOqIiB6QVZyFgIQAlN4vrTdWaiBFQkACCzsiahVY1BERPSBfkY/S+6WIfDYS3c27a4zLKMxA+KFw5CvyWdQRUavAoo6ISI3u5t3hZuXW0mkQETUYlwkjIiIiagdY1BERERG1AyzqiIiIiNoBFnVERERE7QCLOiIiIqJ2gKNfiYiozWnoBNG68jhPNH379m088cQTOHr0KLp169bS6SjNmjULly5dQnx8vNrto0ePxjPPPIOZM2c2b2ItiEUdERG1KdpMEK0rjzrR9Kefforw8HCEhIQgKipK2b5u3TosX74ccrkc/fr1w5o1a+Dh4aGjrP/m4+ODgwcPAgAMDQ3h7OyMiIgIvPbaa7ViJ06cCAcHByxZsgQA8PHHHyMgIEBZ0E2YMAExMTF46623EB0drbJvcHAw1q9fj/Hjx2PLli06/QwP55Weng5vb2+N8fPnz8dzzz2HyZMnw9xct8sl3rx5E++//z5++eUX3Lt3Dy4uLti8eTMGDBigjGmu7/ZBLOqIiKhNaegE0bryqBNN//nnn/jiiy/w5JNPqrTHxcUhLCwM0dHR8PT0RFRUFPz8/HD+/HnY2NhofZ78/HwYGhrC1NRUpV0QBBw/fhyfffYZxo4di9LSUqxevRrjxo2Dl5cXnJyclLGVlZXYvXs39uzZAwC4d+8evvrqK/z6668qx3R0dERsbCxWrlwJqVQKACgrK8O2bdvQpUsXrXOvz8N5AcCJEycwbdo0jfv06dMHzs7O+PbbbxEcHKz1OTVdz/z8fAwaNAhDhgzBL7/8gk6dOuHixYuwtLRUxuj6u22oVtGnbt26dejWrRuMjIzg6emJo0eP1hm/Y8cO9OrVC0ZGRujbty9+/vlnle2CIGDBggWQyWSQSqXw9fXFxYsXldt///13iEQita8///yzST4jERHpVs0E0U39epTCsbi4GGPHjsWmTZtU/tEHgM8//xxTpkzBxIkT4ebmhujoaBgbG+Prr79u8PHv37+PPXv2YNSoUZDJZLh8+XKtmIsXL+Lu3bsYOnQo7Ozs4OTkhEmTJqGyshLnz59XiT18+DAMDQ3xzDPPAAB+/vlnSCQSDBw4UCXu6aefhqOjI3bu3Kls27lzJ7p06YKnnnpKJfb555/H9OnTMX36dJibm8Pa2hoREREQBEEZU1VVhWXLlsHFxQUSiQRdunTBxx9/rDGvGzduIC8vDwDw4osvwtjYGD179kRKSorKuV9++WXExsbq9HouXboUjo6O2Lx5Mzw8PODk5ISXXnoJzs7OyhhdfLeN0eJFXU01u3DhQqSlpaFfv37w8/NDTk6O2vjDhw9jzJgxmDRpEo4fP47AwEAEBgbi1KlTyphly5Zh9erViI6ORkpKCkxMTODn54eysjIAgLe3N7KyslRekydPhpOTk8qtUyIiokcRHBwMf39/+Pr6qrSXl5cjNTVVpV1PTw++vr5ITk6u97gnT57EzJkz0blzZ4wbNw6dOnXC/v370a9fv1qxqampsLS0hJtb9QopN27cwAcffACJRFLr7uGPP/6Il19+GSKRCABw6NAh9O/fX20Ob775JjZv3qx8//XXX2PixIlqY2NiYmBgYICjR49i1apV+Pzzz/Hll18qt4eHh+PTTz9FREQEzpw5g23btsHW1lZjXunp6QCqbwrNmzcPJ06cQJcuXTB37lyV83p4eODo0aNQKBRq86qhzfX88ccfMWDAAIwaNQo2NjZ46qmnsGnTJuX2R/1uH0WLF3XaVrOrVq3C0KFDMXv2bDzxxBP46KOP8PTTT2Pt2rUAqu/SRUVFYf78+QgICMCTTz6Jb775Brdu3VJ2phSLxbCzs1O+rKyskJCQgIkTJyr/whARET2K2NhYpKWlITIysta2vLw8VFZWqhQuAGBrawu5XK72eLdv38aqVavw9NNPY8CAAcjIyMD69euRlZWF9evXw8vLS+1+aWlpKCwsRIcOHSCVSuHo6Ii9e/ciOjoa9vb2KrEJCQkYMWKE8v3Vq1drxdR4/fXX8ccff+Dq1au4evUq/ve//+H1119XG+vo6IiVK1eiZ8+eGDt2LN59912sXLkSAHD37l2sWrUKy5Ytw/jx4+Hs7IzBgwdj8uTJGvNKT09Hx44dsX37dgwZMgSurq4YMWIEcnNzVc5rb2+P8vJytde0sdczIyMDGzZsgKurK3799VdMmzYNM2bMQExMDIDGfbe60qJ96mqq2fDwcGVbfdVscnIywsLCVNr8/PyUBVtmZibkcrlKhWxubg5PT08kJydj9OjRtY75448/4vbt2xp/wwAAhUKhUukXFRU16DMSEdHj5/r16wgJCcHevXthZGSkk2OuWbMGixYtwrPPPotLly7B0dGxQfulpaUhODgYM2bMQEFBAWbNmoVBgwZhwoQJKnFnz57FrVu38MILLyjbSktLNebfqVMn+Pv7Y8uWLRAEAf7+/rC2tlYbO3DgQJWbJl5eXlixYgUqKytx9uxZKBQKlfPWl1d6ejoCAgJUzpeZmQkXFxeVfWv6+927d6/WcRt7PauqqjBgwAB88sknAICnnnoKp06dQnR0NMaPH9+gYzSVFr1T15hqVi6X1xlf81ObY3711Vfw8/ND586dNeYaGRkJc3Nz5auhXz4RET1+UlNTkZOTg6effhoGBgYwMDDAgQMHsHr1ahgYGMDKygr6+vrIzs5W2S87Oxt2dnZqjzl16lR89NFHkMvl6N27NyZOnIh9+/ahqqqqzlzS0tLg7e0NFxcXDBgwAOvXr8fSpUtx5coVlbgff/wRL774okoRZ21tjfx8zVPHvPnmm9iyZQtiYmLw5ptv1nNV1KspvDRRl1d6enqtfn7p6elwd3dXabtz5w6A6gL0YY29njKZTPkou8YTTzyBa9euAai+Ztp+t7rS4o9fW9qNGzfw66+/YtKkSXXGhYeHo7CwUPm6fv16M2VIRERtzQsvvICTJ08iPT1d+RowYADGjh2L9PR0SCQS9O/fH0lJScp9qqqqkJSUpPGxn729PebPn48LFy4gMTERYrEYI0eORNeuXTF37lycPn261j4ZGRkoKChAnz59lG1ubm5wdnbGtm3bVGITEhIQEBCg0vbUU0/hzJkzGj/n0KFDUV5ejoqKCvj5+WmMe3gAw5EjR+Dq6gp9fX24urpCKpWqXIu68rp79y4yMjJqDchQV9SdOnUKnTt3VnsHsTHXEwAGDRpUa4DJhQsX0LVrVwDVXby0/W51pUWLusZUs3Z2dnXG1/xs6DE3b94MKysrlWf16kgkEpiZmam8iIiI1OnQoQP69Omj8jIxMYGVlZWywAoLC8OmTZsQExODs2fPYtq0aSgpKamzK1ANb29vfPHFF5DL5Vi+fDnS09PRr18/nDx5UiUuNTUVhoaG6NGjh0r7Cy+8gF27dinf5+Tk4NixYxg+fLhKnJ+fH06fPq3xbp2+vj7Onj2LM2fOQF9fX2O+165dQ1hYGM6fP4/vvvsOa9asQUhICADAyMgI77//PubMmYNvvvkGly9fxpEjR/DVV1+pzevEiRPQ19dH3759lW1Xr15Ffn5+raLu0KFDeOmllzTmVaOh1xMA3nvvPRw5cgSffPIJLl26hG3btmHjxo0q06Y8ynf7KFq0T92D1WxgYCCAv6vZ6dOnq93Hy8sLSUlJCA0NVbbt3btXWf06OTnBzs4OSUlJyi+3qKgIKSkpteazEQQBmzdvxrhx42BoaKjzz0dERE0nozCjTZ8nKCgIubm5WLBgAeRyOdzd3ZGYmFir+1BdjIyMMHr0aIwePRq3bt2qNadaWloaXF1dIRaLVdp9fX0RHR2NGzduoHPnzvjpp5/g4eFR645W37598fTTT2P79u1466231ObQkJsc48aNQ2lpKTw8PKCvr4+QkBBMnTpVuT0iIgIGBgZYsGABbt26BZlMhrffflttXunp6ejZs6fK49jjx4/DwsJCZcWLsrIyxMfHIzExsd78atR3PQHgmWeewa5duxAeHo7FixfDyckJUVFRGDt2rDJGF99tY7T45MNhYWEYP348BgwYAA8PD0RFRalUs+PGjYODg4Ny9FBISAh8fHywYsUK+Pv7IzY2FseOHcPGjRsBACKRCKGhoViyZAlcXV3h5OSEiIgI2NvbKwvHGvv27UNmZqbKCBuidq/gOnDvtnb7GFsBFuxHSq2DpcQSUgMpwg+F1x+sI1IDKSwllvUH1uH333+v1VYzf5suqBulGhkZqXb07ciRI1X6jj08uvRBCxYswOzZszFlyhTo6enVu1KEumW7DA0NERUVhQ0bNqjdR09PDx988AE++OADlfYRI0bUykvdNauZ3uxBNfPIPdz3rqE0jfoFgOHDh9e6q/kwXX63DdXiRV191ey1a9egp/f3U2Jvb29s27YN8+fPx7x58+Dq6or4+HiV/gJz5sxBSUkJpk6dioKCAgwePBiJiYm1RvB89dVX8Pb2Rq9evZrnwxK1tILrwDoPoKL2SLA6GRoDwUdZ2FGrIDOVISEggWu/6tDgwYMxZswYtdv8/f1x8eJF3Lx5s9kHCdaVV30MDQ2xZs0aHWfUurV4UQfUXc2q+81m1KhRGDVqlMbjiUQiLF68GIsXL67zvA93EiVq9+7dri7oRm4CrHvUHw8AeReAnVOq92VRR62EzFTWrous5jZnzpw6tz/Y5ak51ZdXXR7Hp3CtoqgjomZm3QOwd2/pLIjoMaDu5gw1jcd+ShMiIiKi9oB36ohIp7KKsxrU16m991EiImpuLOqISGeyirMQkBCA0vul9cZKDaRICEgAYNz0iRERPQZY1BGRzuQr8lF6vxSRz0aiu3l3jXEZhRkIPxT+1x09FnVERLrAoo6IdK67eXe4WbnVH0jUQIIgtHQKRE2qvjVnG4JFHRG1Cpdyipv0+JYmYjhY1L1wOLU+hoaGEIlEyM3NRadOnSASiVo6JSKdEgQB5eXlyM3NhZ6eXq3VP7TBoo6IWpSliRhSQ32ExqU36Xmkhvr4baYPC7s2Rl9fH507d8aNGzdw5cqVlk6HqMkYGxujS5cuKgsuaItFHRG1KAcLKX6b6YP8kvImO8elnGKExqUjv6ScRV0bZGpqCldXV1RUVLR0KkRNQl9fHwYGBo98J5pFHZEuabuuat6FpsulDXGwkLLYojrp6+tDX1+/pdMgatVY1BHpSiPXVa0ykOLiXTGyLp/H3fLCOmNvllwFAFzPvwc3AAVXT6G36D6M8swBkenfgUVXqn/mXgAUD9wBYxFJRNRusagj0pWH1lXNKVZg2n9SUXa/7hFN+WUdkPWf4zBxXgGRXv2Pl4QqQyzccxvP6kswIO197JEA2PVQkNgQcJABOycD5Q8d09AYMLbS7rMREVGrx6KOSNf+Wlc152YhUivyEBXkDhcb0zp3ySw6j3lHKxDcewEcTLrWGdtBbA7rl+xwvfgF6JfdgZnUEDamEtWgoitASgQw8kvArJvqNmMrwMJR+89FREStGos6oibmYmOKPg7mdcboGVUXfc859dZifrc6jin5a0h8px5APcerd1mvoiuA2BCWpXngol5ERK0Xizqix1iDl/VykEF6+H0k2O7meq1ERK0Uizqix1iDlvXKvYCMn6Yh3MYa+Yp8FnVERK0UizoiqntZL0U5wPnBiIhavcZPW0xERERErQaLOiIiIqJ2gEUdERERUTvAPnVE7VhGYcYjbW9T6liizSivGL1FmbVX3uCcfUTUjrCoI2qHLCWWkBpIEX4ovN5YqYEUlhLLZsiqCdWzRJsLoH7lDUNjIPgoCzsiahdY1BG1QzJTGRICEuqeVPgvlhLLhk9T8vBasg97eM3Z5roT9tASbQ+7lFuMkNh0rBrtDpdOf92py7sA7JxSvS+LOiJqB1jUETWxvDI5zty+WWdMUzwGlZnKdDOnnLEVYGBU/Wd1a8k+6OE1Z5v7TthfS7Q9rEwoxGmhEGXWfQH7ulf3ICJqq1jUETUhkUEBZh1+HYqqsnpjW+1jUAtH4N//AQ7MUL+W7IMeXHO2vJx3woiImhGLOqImJDIogaKqrO4VG/6i1WPQ5tbBtvpnfWvJPrjmbF2PaYmISOdY1NFjo96F6//SFMVVnSs2UMvKu6D6/uF+geoYWzVpSkREjcGijh4LDV64HtWPQRMCElrvXTPSDWOr6j5/O6eotj/cL1AdQ2Ng7Lamz5GISAss6uix0KCF61E9YCH8UDgXrm9Jdcw3p9HDd9sawsKxehDHw+d6sF+guv6DNaNmywq1PycRURNiUUePFT4GbQHqCi5Njzjv5QFxb2icb65OhsbaPxa1cKw9iOPBfoG6/LvSmGIV4ATJRNRgLOqIqGloerwJ1P2I09AYeP0HwNha+/O11uKnnsmR68QJkomogVjUEVGDabXsmKbHm0Ddjzhbc3HWWPVMjqwRJ0gmIi2wqCOiejV62TF1jzeBpnvE2dppmByZiEgXWNQRNVKtKVKKrlQ/Viy6AkjEyCwqhp44p6XS06kmW3aMiIh0psWLunXr1mH58uWQy+Xo168f1qxZAw8PD43xO3bsQEREBK5cuQJXV1csXboUw4YNU24XBAELFy7Epk2bUFBQgEGDBmHDhg1wdXVVOc6ePXuwePFi/N///R+MjIzg4+OD+Pj4pvqY1M5onCLFQVb9WPEvUgdAomfUOleK0JLOlh1rQS21ZBsRUXNo0aIuLi4OYWFhiI6OhqenJ6KiouDn54fz58/DxsamVvzhw4cxZswYREZGYvjw4di2bRsCAwORlpaGPn36AACWLVuG1atXIyYmBk5OToiIiICfnx/OnDkDI6Pq9St/+OEHTJkyBZ988gn+8Y9/4P79+zh16lSzfnZq29ROkZJ7obrj/8gvgU49cDmnGCFx6Vjzhk+bL4bag3axZBsRUR1atKj7/PPPMWXKFEycOBEAEB0djT179uDrr7/G3Llza8WvWrUKQ4cOxezZswEAH330Efbu3Yu1a9ciOjoagiAgKioK8+fPR0BAAADgm2++ga2tLeLj4zF69Gjcv38fISEhWL58OSZNmqQ8tpvbY9Svh3RGZYoURXn1SE6zboCVG6rKClFVlgtrI7sWzZGqtZsl24iINNDTdofExET88ccfyvfr1q2Du7s7XnvtNeTn19/fpkZ5eTlSU1Ph6+v7dzJ6evD19UVycrLafZKTk1XiAcDPz08Zn5mZCblcrhJjbm4OT09PZUxaWhpu3rwJPT09PPXUU5DJZPjnP/9Z7506hUKBoqIilRcRtT01hXhdLxZ0RNQWaX2nbvbs2Vi6dCkA4OTJk5g5cybCwsKwf/9+hIWFYfPmzQ06Tl5eHiorK2Fra6vSbmtri3PnzqndRy6Xq42Xy+XK7TVtmmIyMqr7y3z44Yf4/PPP0a1bN6xYsQLPP/88Lly4gI4dO6o9d2RkJBYtWtSgz0ZErdflnGJUPcJqEEZ5xXABcD3/Xq3jWZqI4WAh1UWaRERa07qoy8zMVD6q/OGHHzB8+HB88sknSEtLUxmw0FpVVVUBAD744AO8+uqrAIDNmzejc+fO2LFjB9566y21+4WHhyMsLEz5vqioCI6OnDeKqC2wNBFDYlD9YCIkLh1VZbmNPlZvUSb2SIC4//4BdANWx+2BWHEUAGBkoIcNb/SHjalEdafGLGNGRKQlrYs6sViMe/eqf0P97bffMG7cOABAx44dtXokaW1tDX19fWRnZ6u0Z2dnw85OfR8kOzu7OuNrfmZnZ0Mmk6nEuLu7A4Cy/cE+dBKJBN27d8e1a9c05iuRSCCRSDRuJ6LWy8FCig2v98e7B4BVQe5wMuvZ6GMZFjuhaocUs/S2IwgyrBKvgxseWBVjm6YdG7GMGRGRFrQu6gYPHoywsDAMGjQIR48eRVxcHADgwoUL6Ny5c4OPIxaL0b9/fyQlJSEwMBBA9V20pKQkTJ8+Xe0+Xl5eSEpKQmhoqLJt79698PLyAgA4OTnBzs4OSUlJyiKuqKgIKSkpmDZtGgCgf//+kEgkOH/+PAYPHgwAqKiowJUrV9C1a1dtLgURtSE2Hap/KXO2MYWblfkjHMkcmP4nIE9VWRXjUm4xQmLTsWq0O1w6mdberT2ulEFErYrWRd3atWvxzjvv4Pvvv8eGDRvg4OAAAPjll18wdOhQrY4VFhaG8ePHY8CAAfDw8EBUVBRKSkqUo2HHjRsHBwcHREZGAgBCQkLg4+ODFStWwN/fH7GxsTh27Bg2btwIABCJRAgNDcWSJUvg6uqqnNLE3t5eWTiamZnh7bffxsKFC+Ho6IiuXbti+fLlAIBRo0ZpezmI6HFk4QhU3q3+81+rYpQJhTgtFKLMui9g/yhFIxFR42hd1HXp0gW7d++u1b5y5UqtTx4UFITc3FwsWLAAcrkc7u7uSExMVA50uHbtGvT0/h6g6+3tjW3btmH+/PmYN28eXF1dER8fr5yjDgDmzJmDkpISTJ06FQUFBRg8eDASExOVc9QBwPLly2FgYIA33ngDpaWl8PT0xL59+2BpyXmpiIiIqG3SuqjT19dHVlZWrcmBb9++DRsbG1RWVmp1vOnTp2t83Pr777/Xahs1alSdd9REIhEWL16MxYsXa4wxNDTEZ599hs8++0yrXOkxcje7et45TYquVP/MvfB3HDvDa01XqzdwXjkiokYUdYIgqG1XKBQQi8WPnBBRq7D9DeDeXc3bxYbVS4LtnFw94XANdoZvEEuJJaQGUoQfCtfJ8aQGUiQEJLCwI6LHWoOLutWrVwOovhP25ZdfwtT0747AlZWVOHjwIHr16qX7DIlawv0yYOQmwLqH+u1FV1Q6ySuxM3yDyExlSAhIQL6i4ROWa5JRmIHwQ+HIV+SzqCOix1qDi7qaPnOCICA6Ohr6+vrKbWKxGN26dUN0dLTuMyRqKdY9AHt39dskf92V/quTPGlPZipjEUZEpEMNLuoyMzMBAEOGDMHOnTs5qICIiIioFdG6T93+/fubIg8iIiIiegRaF3WVlZXYsmULkpKSkJOTo1x2q8a+fft0lhwRERERNYzWRV1ISAi2bNkCf39/9OnTByKRqCnyImpRGYaG1YMhJOpHdOtqKg4iIiJd0bqoi42Nxfbt2zFs2LCmyIeoRVlKLCHVEyPcxrp6dGsdpAZSWErYt5SIiFoHrYs6sVgMFxeXpsiFqFndLChFfsnDEwwbY1WvhTDf+yauD1mNcgvNf9c7iM1xu9AYtwsL1W6/lFOsw2yJ/pZVnNWg6WA4KTPR40Xrom7mzJlYtWoV1q5dy0ev1GbdLCiF74oDKK2ovQJKb9E17JFUYPaPd3BayK3jKLkALtV5HqmhPixNOCk36U5WcRYCEgJQer+03lhOykz0eNG6qPvjjz+wf/9+/PLLL+jduzcMDQ1Vtu/cuVNnyRE1lfyScpRWVCIqyB0uNqYq24zyzIFdwKrR7tWLsz8CSxMxHCykj3QMogflK/JRer8Ukc9Gort5d41xnJSZ6PGjdVFnYWGBV155pSlyIWp2Ljam6ONgrtooqi7yXDqZAvbmavYianndzbvDjRNfE9EDtC7qNm/e3BR5EBG1STUjoTOLiqFndBOZReehZ6R695d924ioOWhd1BER0V8jpQ2kCD8UrmwzcQLmHa0dy75tRNQctC7qnJyc6hwgkZHB+buIqP2TmcqQEJCgHIV6OacYIXHpWBXkDucH+mmybxsRNReti7rQ0FCV9xUVFTh+/DgSExMxe/ZsXeVFRKSV+iaEbooJo2WmMmWhVlVWiKqyXDiZ9YSbFftiElHza9SKEuqsW7cOx44de+SEiIi0oe4xqCacMJqI2jOd9an75z//ifDwcA6kIKJm9fBj0LpwwAIRtWc6K+q+//57dOzYUVeHIyJqsAcfgxIRPa60LuqeeuoplYESgiBALpcjNzcX69ev12lyRETtRUP69PFOIhE9Cq2LusDAQJX3enp66NSpE55//nn06tVLV3kREbUL2vb549QnRNRYWhd1CxcubIo8iIjapYb2+ePUJ0T0qBrVp66yshLx8fE4e/YsAKB3794YMWIE9PX1dZocEVF7wD5/RNQctC7qLl26hGHDhuHmzZvo2bMnACAyMhKOjo7Ys2cPnJ2ddZ4kEREREdVNT9sdZsyYAWdnZ1y/fh1paWlIS0vDtWvX4OTkhBkzZjRFjkRERERUD63v1B04cABHjhxRmb7EysoKn376KQYNGqTT5IiIiIioYbS+UyeRSHD37t1a7cXFxRCLxTpJioiIiIi0o3VRN3z4cEydOhUpKSkQBAGCIODIkSN4++23MWLEiKbIkYiIiIjqofXj19WrV2P8+PHw8vKCoaEhAOD+/fsYMWIEVq1apfMEiYgee3kX/v5z0ZXqn7kXAEW55n3K6182jYjaF62LOgsLCyQkJODSpUvKKU2eeOIJuLi46Dw5IqLHmrEVYGgM7Jzyd5vYEHCQATsnA+UVdezbAbC1bPociajVaPTary4uLizkiIiakoUjEHwUuHf777aiK0BKBDDyS8Csm/r98i4Au99pjgyJqBXRuqh79dVX4eHhgffff1+lfdmyZfjzzz+xY8cOnSVHRPTYs3CsftWQ/DUgrVMPwMqtZXIiolZJ64ESBw8exLBhw2q1//Of/8TBgwd1khQRERERaUfrok7T1CWGhoYoKipqVBLr1q1Dt27dYGRkBE9PTxw9erTO+B07dqBXr14wMjJC37598fPPP6tsFwQBCxYsgEwmg1Qqha+vLy5evKgS061bN4hEIpXXp59+2qj8iYiIiFqa1kVd3759ERcXV6s9NjYWbm7aPwqIi4tDWFgYFi5ciLS0NPTr1w9+fn7IyclRG3/48GGMGTMGkyZNwvHjxxEYGIjAwECcOnVKGbNs2TKsXr0a0dHRSElJgYmJCfz8/FBWVqZyrMWLFyMrK0v5evfdd7XOn4iIiKg10LpPXUREBEaOHInLly/jH//4BwAgKSkJ3333XaP6033++eeYMmUKJk6cCACIjo7Gnj178PXXX2Pu3Lm14letWoWhQ4di9uzZAICPPvoIe/fuxdq1axEdHQ1BEBAVFYX58+cjICAAAPDNN9/A1tYW8fHxGD16tPJYHTp0gJ2dndY5ExEREbU2Whd1L7/8MuLj4/HJJ5/g+++/h1QqxZNPPonffvsNPj4+Wh2rvLwcqampCA8PV7bp6enB19cXycnJavdJTk5GWFiYSpufnx/i4+MBAJmZmZDL5fD19VVuNzc3h6enJ5KTk1WKuk8//RQfffQRunTpgtdeew3vvfceDAwaPSCYiAiXcoobtV9mUfV+BzNP43Idxyi4f7NRxyei9q9RFYy/vz/8/f0f+eR5eXmorKyEra2tSrutrS3OnTundh+5XK42Xi6XK7fXtGmKAYAZM2bg6aefRseOHXH48GGEh4cjKysLn3/+udrzKhQKKBQK5fvG9h8kovbJ0kQMqaE+QuPSG7W/yKAAJs6GWHd6cb2xEn0jWEo4Bx0RqXpsb0s9eLfvySefhFgsxltvvYXIyEhIJJJa8ZGRkVi0aFFzpkhEbYiDhRS/zfRBfkkdqzzUI69sAO6WF2rcfv3OPXy29wLWvOEDmams0echovapRYs6a2tr6OvrIzs7W6U9OztbY183Ozu7OuNrfmZnZ0Mmk6nEuLu7a8zF09MT9+/fx5UrV9CzZ89a28PDw1UKwaKiIjg6OtaKo1ao4Hr1kkqAcmklo7xi9BZlwijPHBCZqsY/uCQTkRYcLKRwsJA+whHM69x66mYhlpWVwNqIfYGJqLYWLerEYjH69++PpKQkBAYGAgCqqqqQlJSE6dOnq93Hy8sLSUlJCA0NVbbt3bsXXl5eAAAnJyfY2dkhKSlJWcQVFRUhJSUF06ZN05hLeno69PT0YGNjo3a7RCJRewePWrmC68A6D0BUobK0kguAPRIAuzTsZ2hcvUQTERFRG9Hij1/DwsIwfvx4DBgwAB4eHoiKikJJSYlyNOy4cePg4OCAyMhIAEBISAh8fHywYsUK+Pv7IzY2FseOHcPGjRsBACKRCKGhoViyZAlcXV3h5OSEiIgI2NvbKwvH5ORkpKSkYMiQIejQoQOSk5Px3nvv4fXXX4elJfuptCv3bgMV94Chi4HzXyqXVrqUW4yQ2HSsGu0Ol06mtfcztlKdxZ+IiKiVa/GiLigoCLm5uViwYAHkcjnc3d2RmJioHOhw7do16On9PZ2et7c3tm3bhvnz52PevHlwdXVFfHw8+vTpo4yZM2cOSkpKMHXqVBQUFGDw4MFITEyEkZERgOq7brGxsfjwww+hUCjg5OSE9957r9aoWmpHLLpW//xraaUyoRCnhUKUWfcF7Ot+5EVERNQWaF3UVVZWYsuWLUhKSkJOTg6qqqpUtu/bt0/rJKZPn67xcevvv/9eq23UqFEYNWqUxuOJRCIsXrwYixerH0X29NNP48iRI1rnSURERNRaaV3UhYSEYMuWLfD390efPn0gEomaIi8iIiIi0oLWRV1sbCy2b9+OYcOGNUU+RERERNQIWq/9KhaL4eLi0hS5EBEREVEjaV3UzZw5E6tWrYIgCE2RDxERERE1gtaPX//44w/s378fv/zyC3r37g1DQ0OV7Tt37tRZckRERETUMFoXdRYWFnjllVeaIhciIiIiaiSti7rNmzc3RR5ERERE9AhafPJhokeVVZyFfEW++o1FVwCxITJKbjZrTkRERM2tQUXd008/jaSkJFhaWuKpp56qc266tLQ0nSVHVJ+s4iwEJASg9H6p5iAHGXAqGlIDKSwlXAaOiIjapwYVdQEBAcrF7GvWTyVqDfIV+Si9X4rIZyPR3bx77YDcC8DOycDIL2HZ2RMyU1nzJ0lERNQMGlTULVy4UO2fiVqL7ubd4WblVnuDohworwDMugEs6IiIqB1rdJ+6Y8eO4ezZswAANzc39O/fX2dJEREREZF2tC7qbty4gTFjxuB///sfLCwsAAAFBQXw9vZGbGwsOnfurOsciYiIiKgeWq8oMXnyZFRUVODs2bO4c+cO7ty5g7Nnz6KqqgqTJ09uihyJiIiIqB5a36k7cOAADh8+jJ49eyrbevbsiTVr1uDZZ5/VaXJERERE1DBa36lzdHRERUVFrfbKykrY29vrJCkiIiIi0o7WRd3y5cvx7rvv4tixY8q2Y8eOISQkBJ999plOkyMiIiKihmnQ41dLS0uVCYdLSkrg6ekJA4Pq3e/fvw8DAwO8+eabnMeOiIiIqAU0qKiLiopq4jSIiIiI6FE0qKgbP358U+dBRERERI+g0ZMPExERNcbNglLkl5Q36TksTcRwsJA26TmIWhsWdURE1GxuFpTCd8UBlFZUQmRQAJFBSb37CPdNINy30Oo8UkN9/DbTh4UdPVZY1BERtTGXcorr3G6U9/f2yznFqCorbPCxm/oOV35JOUorKvHhK/ZYe+FDKKrK6t1HomeEz7y/hbWRXYPOcSmnGKFx6cgvKWdRR48VFnVERG2EpYkYUkN9hMal1xnXW5SJ5R2q/xwSl46qstwGn6O57nBZmVVAUVWGyGcj0d28u8a4jMIMhB8Kh41FJdyszJs0J6K2jkUdtR0F14F7t1Xbiq5U/8y9ACjU9NHJu9DkaRE1FwcLKX6b6VNvfzSjPHOU76n+86ogdziZ9awzvkZL3OHqbt4dblZuzXIuovauQUXdyJEjG3zAnTt3NjoZIo0KrgPrPICKe6rtYkPAQQbsnAyU117pBABgaAwYWzV9jkTNwMFCWn/BJTLFmb/+6GxjyjtcRI+JBhV15uZ//w9BEATs2rUL5ubmGDBgAAAgNTUVBQUFWhV/RFq5d7u6oBu5CbDu8Xd70RUgJQIY+SVg1k39vsZWgIVjc2RJ1PpououthmGxuImTIaKm1KCibvPmzco/v//++/j3v/+N6Oho6OvrA6he9/Wdd96BmZlZ02RJVMO6B2Dv/vd7yV//CHXqAfARDlE1YyvAwKj6z3XdxX6Iq4EU9ljahIkRUVPSuk/d119/jT/++ENZ0AGAvr4+wsLC4O3tjeXLl+s0QSIi0pKFI/Dv/wAHZtR9F/svWaV5yM89Dez7CCZGGcgsOg89I9NacZYSS8hMZU2UNBE9Kq2Luvv37+PcuXPo2VO14+25c+dQVVWls8SIiOgRdLCt/lnPXeys4iwE7J+K0vul1f1TsRPzjqrvGy01kCIhIIGFHVErpXVRN3HiREyaNAmXL1+Gh4cHACAlJQWffvopJk6cqPMEiYio6eQr8lF6vxSRfd5G918XIqQ8GDOC/OFso3qnrmZqkXxFPos6olZK66Lus88+g52dHVasWIGsrCwAgEwmw+zZszFz5kydJ0hERE2vu4kD3MorIFZ0gpNZT46YJWqDtC7q9PT0MGfOHMyZMwdFRUUAwAESRERERC1MrzE73b9/H7/99hu+++47iEQiAMCtW7dQXFz30jVERERE1DS0LuquXr2Kvn37IiAgAMHBwcjNrV5+ZunSpZg1a1ajkli3bh26desGIyMjeHp64ujRo3XG79ixA7169YKRkRH69u2Ln3/+WWW7IAhYsGABZDIZpFIpfH19cfHiRbXHUigUcHd3h0gkQnp6eqPyJyJqrTIKM3Dm9hmNr4zCDJV4F9FNGOWdBG6lq75y/1qdJfdC7W0F15vxExGRJlo/fg0JCcGAAQNw4sQJWFn9PUv/K6+8gilTpmidQFxcHMLCwhAdHQ1PT09ERUXBz88P58+fh42NTa34w4cPY8yYMYiMjMTw4cOxbds2BAYGIi0tDX369AEALFu2DKtXr0ZMTAycnJwQEREBPz8/nDlzBkZGRirHmzNnDuzt7XHixAmtcyciaq0sJZaQGkgRfii83lipgRSWZo6oMpBiFdYDu9bXDqpr9RZDYyD4KCf5JmphWhd1hw4dwuHDhyEWq8483q1bN9y8eVPrBD7//HNMmTJFOXI2Ojoae/bswddff425c+fWil+1ahWGDh2K2bNnAwA++ugj7N27F2vXrkV0dDQEQUBUVBTmz5+PgIAAAMA333wDW1tbxMfHY/To0cpj/fLLL/jvf/+LH374Ab/88ovWuRMRtVYyUxkSAhKQr8ivN7Zm/rnzo/YhbMs+rBrtDpdOD81Tp2n1lrwLwM4p1au+tLOiLqs4S6vrR9TStC7qqqqqUFlZWav9xo0b6NChg1bHKi8vR2pqKsLD//5NUk9PD76+vkhOTla7T3JyMsLCwlTa/Pz8EB8fDwDIzMyEXC6Hr6+vcru5uTk8PT2RnJysLOqys7MxZcoUxMfHw9jYuN5cFQoFFAqF8n3NIBEiotZKZirTqtioMHXAacEJZdZ9AfuHRr8+Zqu3ZBVnISAhoHr+vnpw/j5qLbQu6l566SVERUVh48aNAACRSITi4mIsXLgQw4YN0+pYeXl5qKyshK2trUq7ra0tzp07p3YfuVyuNl4ulyu317RpihEEARMmTMDbb7+NAQMG4MqVK/XmGhkZiUWLFjXocxERtVcP98FD0ZXqR7NFV5SFX3u4c6Wcv+/ZSHQ3764xjvP3UWuidVG3YsUK+Pn5wc3NDWVlZXjttddw8eJFWFtb47vvvmuKHHVuzZo1uHv3rsodwvqEh4er3CEsKiqCo2P7etRARKRJnX30HGTVj2b/ItUTI8F7KWRS61qhRnnF6C3KhEGp0JTp6kx38+5wewzuTFL7oHVR17lzZ5w4cQJxcXE4ceIEiouLMWnSJIwdOxZSqVSrY1lbW0NfXx/Z2dkq7dnZ2bCzs1O7j52dXZ3xNT+zs7Mhk8lUYtzd3QEA+/btQ3JyMiQSicpxBgwYgLFjxyImJqbWeSUSSa14IqLHhcY+enezge1vAPfLAAAZhoYIt7FGftxoyB4eUAHABcAeCXB6bwfAzrIZMid6fGhd1AGAgYEBxo4di7Fjxz7SycViMfr374+kpCQEBgYCqO6zl5SUhOnTp6vdx8vLC0lJSQgNDVW27d27F15eXgAAJycn2NnZISkpSVnEFRUVISUlBdOmTQMArF69GkuWLFHuf+vWLfj5+SEuLg6enp6P9JmIiNortX30rNyAqcnVAyUA5YCKDL9FgIlDrWNcz7+HuP/+gZF6u5o+YaLHjNZFnb6+Pp577jn88MMP6Nixo7I9Ozsb9vb2agdR1CUsLAzjx4/HgAED4OHhgaioKJSUlChHw44bNw4ODg6IjIwEUD2lio+PD1asWAF/f3/Exsbi2LFjKn38QkNDsWTJEri6uiqnNLG3t1cWjl26dFHJwdS0epSXs7MzOnfurO0lISJ6vFk4Kke+WprZQpoqRfipaM3x3YA/YQ2pnhiWEt6tI9IVrYs6QRCgUCgwYMAA/PTTT+jdu7fKNm0FBQUhNzcXCxYsgFwuh7u7OxITE5UDHa5duwY9vb/nSPb29sa2bdswf/58zJs3D66uroiPj1fOUQdUzz1XUlKCqVOnoqCgAIMHD0ZiYmKtOeqIiEi36ptK5XJOMVbH7cEq8TpYBsXqdHBBzRQkmUXF0DO6icyi89AzMq0V1x4GchCpo3VRJxKJ8MMPP+DTTz+Fl5cX/vOf/yjng6tZMkxb06dP1/i49ffff6/VNmrUKIwaNarOHBcvXozFixc36PzdunVrVEFKTedmQSku3r6Gu+WFAABxwSU4ig1x/eoZlBeW/x1XchVA9T8UVWWFDT7+pRwuaUfUVOqaSqWqrBBixVG4oQJQM5CisR6egsTECZinYXEiTkFC7VWj7tTp6+tj1apV6N27N4KCgjB//nxMnjy5KfKjx9DNglL4roqHfpflEOk90NHaQQZciKwVL1QZYsbWixDu52p1HqmhPixNxPUHElGr9+AUJILCBiFx6VgV5A5nG9U7dZyChNqzRg2UqDF16lS4urpi1KhROHjwoK5yosdcfkk5FMJdmOhVILj3AjiYdK2+U7d/Bq4PWY1yCxeV+A5ic1i/pH60dF0sTcRwsNBuxDYRtbCC638PynhQ0RUAQPfy+xCXl+MJRTmeKC+Hi6IcMLZqd6tdEKmjdVHXtWtX6OvrK98PGTIER44cwcsvv6zTxIgA4Dmn3tVzRN0SA+UVcOvqBti7t3RaRNQSCq4D6zyAinu1tz2wNq1LeQX2SADUDLDl2rT0mNC6qMvMzKzV5uLiguPHj9eaP46IiEhn7t2uLuhGbgKse6hue2Bt2ksKa4TEplevYSu61W7XpiV62CM9fn2QkZERunbtqqvDERERqWfdo/Yd+wfWpi0rc8BpobB6DVtR7dGvRO1Vg4q6jh074sKFC7C2toalpWWdo1zv3Lmjs+SIiKgdy7tQf8xffeWQewEoL68zlOhx16CibuXKlejQoQMAICoqqinzISKidi5f6IAqAyn0dk6pP/iBvnIor6juH2ds1fRJErVBDSrqxo8fr/bPRERE2roFa1wctQ89OzTgztsDfeVg1o0jWYnq0KCirqioqMEHNDMza3QyRET0eKgwdQDszesPfKCvHKzcmjYpojauQUWdhYVFvatFCIIAkUik9dqvRERERPToGlTU7d+/v6nzICIiIqJH0KCizsfHp6nzICIi0iijMOORthM9Dho9T929e/dw7do1lD80xPzJJ5985KSIiIgAwFJiCamBFOGHwuuNlRpIYSmxxO2yZkiMqBXSuqjLzc3FxIkT8csvv6jdzj51RESkKzJTGRICEpCvyK831lJiCZmpDLcLC2tvrJkT78F57xR1jL6tiSNqQ7Qu6kJDQ1FQUICUlBQ8//zz2LVrF7Kzs7FkyRKsWLGiKXIkIqLHmMxUBpmpTOv9LuUUw9BYDNcH58R7eN47Tf6Ku3o1A1VlDmpDLE3EWudE1JS0Lur27duHhIQEDBgwAHp6eujatStefPFFmJmZITIyEv7+/k2RJxERUYNYmoghNdRHaFw6AMAeS2EpugsAKEcugJ0IKQ+GWNFJ4zFsRWcAHMLnP/6Ji2XqZ3+QGupjwyTb2hsKrlevNasNzr9HOqB1UVdSUgIbGxsAgKWlJXJzc9GjRw/07dsXaWlpOk+QiIhIGw4WUvw20wf5JbUfr2YWnce8ozsxI8gfTmY9NR7j5rU9+PPCIczy6wmHLoNrbb+UU4zQuHTcLX3obl/BdWCdB1BxT7ukDY2B4KMs7OiRaF3U9ezZE+fPn0e3bt3Qr18/fPHFF+jWrRuio6Mhk2l/e5yIiEjXHCykcLCQ1mrXMzIFADjbmMLNSvPkx3p3jQEAjpbGcHNowCTJNe7dri7oRm4CrHs0bJ+8C8DOKdX7sqijR6B1URcSEoKsrCwAwMKFCzF06FBs3boVYrEYW7Zs0XV+RERErdb1O9V35C7nFKOqrBBGecVwAXBJsEeZ4NSgYxgJf+2TW4wyQXWQh6WJWG1xSqSO1kXd66+/rvxz//79cfXqVZw7dw5dunSBtbW1TpMjIiJqCvXOe1dys/oPBVeBW+m1ttuUKNDf8Cp27c0GnIDVcXsgVhyFi+gmVomBkNh0nBbUjMJVo7coE3sk6veRGurjt5k+LOyoQRo9T10NY2NjPP3007rIhYiIqElpNe9dlQDLvYsANVN12QD4QR84IzZEEGRYJV4HN1T3r6sykOLzMf+oXt+2AYzyzIFdwKrR7iiz7qtsr+m3l19SzqKOGkTrok4QBHz//ffYv38/cnJyUFVVpbJ9586dOkuO2qes4qw655zKLCqGnjinGTMioseFVvPeKUohg37dQUVXgJQIYOSXgFk3AICesRV6atM3TlTdz8+lkylgr0X/PaKHNGqeui+++AJDhgyBra0tRCL1Q72J1MkqzkJAQgBK75fWGSd1ACR6RrCUWDZTZkT0uGjsvHdqSf6aq65TD8DKTTfHJGokrYu6//znP9i5cyeGDRvWFPlQO5evyEfp/VJEPhuJ7ubd1cZczilGSFw61rzho7v/8RIREbVzWhd15ubm6N5d/T/GRA3V3bw73DT8VltVVoiqslxYG9k1c1ZERERtl9ZF3YcffohFixbh66+/hlTKjptERO3NpZziNnlsosed1kXdv//9b3z33XewsbFBt27dYGhoqLKdq0oQEbVNDy+v1VSkhvpcN5WoCWhd1I0fPx6pqal4/fXXOVCCGu9uNqCovYQPABjlFaO3KLN6mP9fo8KQd6EZkyN6PNW1vJYucUJdoqahdVG3Z88e/Prrrxg8uPZaeEQNtv0N4N5dtZtcAOyRANj10AZD4+pFr4moyWhaXouIWj+tizpHR0eYmZk1RS70OLlfpnFtxEu5xQiJTceq0e7V8zbVMLbiuohEREQaaF3UrVixAnPmzEF0dDS6devWBCnRY8O6B2DvXqu5TCjEaaGwemZ1TsRJRETUII1a+/XevXtwdnaGsbFxrYESd+7c0VlyRERERNQwWhd1UVFRTZAGERFR25VRmFFvjKXEkhOqU5PSqqirqKjAgQMHEBERAScnJ50lsW7dOixfvhxyuRz9+vXDmjVr4OHhoTF+x44diIiIwJUrV+Dq6oqlS5eqrHAhCAIWLlyITZs2oaCgAIMGDcKGDRvg6uqqjBkxYgTS09ORk5MDS0tL+Pr6YunSpbC3t9fZ5yIiovbNUmIJqYEU4YfC642VGkiREJCgdWHXlHP7cSRy+6JVUWdoaIgffvgBEREROksgLi4OYWFhiI6OhqenJ6KiouDn54fz58/DxsamVvzhw4cxZswYREZGYvjw4di2bRsCAwORlpaGPn36AACWLVuG1atXIyYmBk5OToiIiICfnx/OnDkDIyMjAMCQIUMwb948yGQy3Lx5E7NmzcK//vUvHD58WGefjYiI2jeZqQwJAQnIV+TXGZdRmIHwQ+FIzUlFd8VDqzIVXQHEhtU/JWLlHb3mmDdQaqiP32b6sLBrJ0SCIAja7DB+/Hi4u7vjvffe00kCnp6eeOaZZ7B27VoAQFVVFRwdHfHuu+9i7ty5teKDgoJQUlKC3bt3K9sGDhwId3d3REdHQxAE2NvbY+bMmZg1axYAoLCwELa2ttiyZQtGjx6tNo8ff/wRgYGBUCgUtfoJqlNUVARzc3MUFhZyNLAWztw+g6DdQYi7mQW3Cb+pHShx6mYhhq/5A7vfHYw+DhwoQURtX1ZxFgISAlB6v7Te2Afv6N0sKG2yeQMv5RQjNC6d/69tAxpac2jdp87V1RWLFy/G//73P/Tv3x8mJiYq22fMmNHgY5WXlyM1NRXh4X/fttbT04Ovry+Sk5PV7pOcnIywsDCVNj8/P8THxwMAMjMzIZfL4evrq9xubm4OT09PJCcnqy3q7ty5g61bt8Lb27tBBR0REZE26ryjl3sB2DkZGPklMsQGCD8UjnxFPmSmMs4bSFrRuqj76quvYGFhgdTUVKSmpqpsE4lEWhV1eXl5qKyshK2trUq7ra0tzp07p3YfuVyuNl4ulyu317Rpiqnx/vvvY+3atbh37x4GDhyocvfvYQqFAgqFQvm+qKionk9HRET0N5mpTH1/OkU5UF4BmHUDJFw+jRpPT9sdMjMzNb4yMuof/dOazJ49G8ePH8d///tf6OvrY9y4cdD0NDoyMhLm5ubKl6MjJ8ElIiKi1kPrO3UPqimAGrv+q7W1NfT19ZGdna3Snp2dDTs7O7X72NnZ1Rlf8zM7OxsymUwlxt3dvdb5ra2t0aNHDzzxxBNwdHTEkSNH4OXlVeu84eHhKo99i4qKWNgREVGT4BQp1BiNKuq++eYbLF++HBcvXgQA9OjRA7Nnz8Ybb7yh1XHEYjH69++PpKQkBAYGAqgeKJGUlITp06er3cfLywtJSUkIDQ1Vtu3du1dZiDk5OcHOzg5JSUnKIq6oqAgpKSmYNm2axlyqqqoAQOUR64MkEgkkEolWn4+IiEgbzTFFCrVfWhd1n3/+OSIiIjB9+nQMGjQIAPDHH3/g7bffRl5entajYsPCwjB+/HgMGDAAHh4eiIqKQklJCSZOnAgAGDduHBwcHBAZGQkACAkJgY+PD1asWAF/f3/Exsbi2LFj2LhxI4Dqu4ahoaFYsmQJXF1dlVOa2NvbKwvHlJQU/Pnnnxg8eDAsLS1x+fJlREREwNnZWe1dOiIioiaXdwEyAAkDP0F+Rd1z02WU3ET4qWjk551jUUdKWhd1a9aswYYNGzBu3Dhl24gRI9C7d298+OGHWhd1QUFByM3NxYIFCyCXy+Hu7o7ExETlQIdr165BT+/vrn/e3t7Ytm0b5s+fj3nz5sHV1RXx8fHKOeoAYM6cOSgpKcHUqVNRUFCAwYMHIzExUTlHnbGxMXbu3ImFCxeipKQEMpkMQ4cOxfz583k3joiImpexFWBoDOycAgCQ/fWqk9gQcJAB298ApiYDFuwORI0o6rKysuDt7V2r3dvbG1lZWY1KYvr06Roft/7++++12kaNGoVRo0ZpPJ5IJMLixYuxePFitdv79u2Lffv2NSpXIiIinbJwBIKPAvduN3yfoitASgRwv6x6PxZ1hEYUdS4uLti+fTvmzZun0h4XF6eyDBcRERE1kIWjdoUZpz4hNbQu6hYtWoSgoCAcPHhQ2afuf//7H5KSkrB9+3adJ0hERERE9dN6nrpXX30VKSkpsLa2Rnx8POLj42FtbY2jR4/ilVdeaYociYiIiKgejZrSpH///vj22291nQsRERERNZLWd+qIiIiIqPVp8J06PT29eleOEIlEuH///iMnRURERLqRVZyFfEV+rfbMomLoGd1EZtF56BmZcoWKdqDBRd2uXbs0bktOTsbq1auVqzIQERFRy8sqzkJAQgBK75eq3W7iBMw7Wv1nrlDR9jW4qAsICKjVdv78ecydOxc//fQTxo4dq3FeOCIiImoieRc0bsovuoLS+6WI7PM2ups4VDcamQMdbHE5pxghcelYFeQOkSQH4YfCka/IZ1HXhjVqoMStW7ewcOFCxMTEwM/PD+np6SorOhAREVEzMDBSrkSh1l8rT3T/dSHcyiuq2wyNgeCjqDJzQFVZLpzMekLPyLR58qUmpVVRV1hYiE8++QRr1qyBu7s7kpKS8OyzzzZVbkRERFSXf/8HEFtq3l6z8sTILwGzbtV39XZO+Wv1CrNmSpKaS4OLumXLlmHp0qWws7PDd999p/ZxLBERETWjDraAlZvm7TUrT3TqUXcctQsNLurmzp0LqVQKFxcXxMTEICYmRm3czp07dZYctS2aRlg9KKMwo5myISIierw0uKgbN25cvVOa0OOrvhFWD5LqiWFZyZHSRESPqr5flPmL9OOlwUXdli1bmjANauvyFfnVI6yejUR38+51xloW3ILs8qvNlBkRUftjKbGE1ECK8EPh9cZKDaSwlNTR747ajUaNfiXSpLt5d7jV129DUd48yRARtVMyUxkSAhLq7fICgJMKP0ZY1BEREbVBMlMZizVSwbVfiYiIiNoBFnVERERE7QCLOiIiIqJ2gEUdERERUTvAgRLUKDcLSpFf8vco1syiYgDA5ZxiVJUV1rmvUV4xXABcyi1GmVA79lJOsU5zJSIiehywqCOt3Swohe+KAyitqFS26RndhIkTEBKXjqqy3Dr37y3KxB4JEBKbjtNqijoAkBrqw9JErNO8iYiI2jMWdaS1/JJylFZUIirIHS42pgCAzKLzmHcUWBXkDieznnXub5RnDuwCVo12R5l1X7UxliZiOFhIdZ47ERFRe8WijhrNxcYUfRzMAQB6RtXFnbONKdyszOveUVQd69LJFLCvJ5aIiIgahAMliIiIiNoBFnVERERE7QCLOiIiIqJ2gEUdERERUTvAoo6IiIioHWBRR0RERNQOsKgjIiIiagc4Tx0REdHjJu8CjIRi9BZlwijPHOWSvOr23AuAolz9PsZWgIVj8+VIWmNRR0RE9LgwtgIMjYGdU+ACYI8EwC7gjNgQcJABOycD5RXq9zU0BoKPsrBrxVjUERERPS4sHKsLs3u3cSm3GCGx6Vg12h2Q5AEpEcDILwGzbrX3y7sA7JwC3LvNoq4VaxV96tatW4du3brByMgInp6eOHr0aJ3xO3bsQK9evWBkZIS+ffvi559/VtkuCAIWLFgAmUwGqVQKX19fXLx4Ubn9ypUrmDRpEpycnCCVSuHs7IyFCxeivFzDLWciIqL2wsIRsK9ee/u04FS9BnenHtXbOvUA7N1rv6x7tFi61HAtfqcuLi4OYWFhiI6OhqenJ6KiouDn54fz58/DxsamVvzhw4cxZswYREZGYvjw4di2bRsCAwORlpaGPn36AACWLVuG1atXIyYmBk5OToiIiICfnx/OnDkDIyMjnDt3DlVVVfjiiy/g4uKCU6dOYcqUKSgpKcFnn33W3JegbSu4Xv2bW9GV6vd19ceokXehydMiIiLtZRRmqN9QdAUQGwJFV2BpZguZqaxZ86KGEQmCILRkAp6ennjmmWewdu1aAEBVVRUcHR3x7rvvYu7cubXig4KCUFJSgt27dyvbBg4cCHd3d0RHR0MQBNjb22PmzJmYNWsWAKCwsBC2trbYsmULRo8erTaP5cuXY8OGDcjI0PAX+iFFRUUwNzdHYWEhzMzMtP3Ybdqpm4UYvuYP/DrRCT2/fwGouIczYkMEOcgQdzMLbpr6YzyIfTOIiFpUzf/Ld787GFbm9xCQEIDS+6X17ic1kCIhIIGFXTNqaM3RonfqysvLkZqaivDwcGWbnp4efH19kZycrHaf5ORkhIWFqbT5+fkhPj4eAJCZmQm5XA5fX1/ldnNzc3h6eiI5OVljUVdYWIiOHTtqzFWhUEChUCjfFxUV1fv52jv9sjtAxT1g5CZALK67P8bDOIqKiKjVkJnKkBCQgHxFvvqA3AvAzsnI8FuE8FPRyFfks6hrhVq0qMvLy0NlZSVsbW1V2m1tbXHu3Dm1+8jlcrXxcrlcub2mTVPMwy5duoQ1a9bU+eg1MjISixYtqvsDPa6sewAScfWfO/UArNxaNh8iItKazFSmuVBTlFePijVxaN6kSCst3qeupd28eRNDhw7FqFGjMGXKFI1x4eHhKncIi4qK4OjIO01ERNS2XcoprjfGKK8YLgCyr50HAFy/kA6x8Z1696s06ghTWyc4WEgfNU1qgBYt6qytraGvr4/s7GyV9uzsbNjZ2andx87Ors74mp/Z2dmQyWQqMe7u7ir73bp1C0OGDIG3tzc2btxYZ64SiQQSiaRBn4uIiKi1szQRQ2qoj9C49Hpj7ZGH3yQS2B5bDjjI4Lh/Blwa0H/6niDB8KrP8Z+Zr7KwawYtWtSJxWL0798fSUlJCAwMBFA9UCIpKQnTp09Xu4+XlxeSkpIQGhqqbNu7dy+8vLwAAE5OTrCzs0NSUpKyiCsqKkJKSgqmTZum3OfmzZsYMmQI+vfvj82bN0NPr1XM7kJERNQsHCyk+G2mD/JLGjad1/XiZ3D9TjpwIRLXh6yG2LhLnfGSgktw3B8C6f0C5JeUs6hrBi3++DUsLAzjx4/HgAED4OHhgaioKJSUlGDixIkAgHHjxsHBwQGRkZEAgJCQEPj4+GDFihXw9/dHbGwsjh07przTJhKJEBoaiiVLlsDV1VU5pYm9vb2ycLx58yaef/55dO3aFZ999hlyc3OV+Wi6Q0hERNTeOFhItSi2zFF5G8AFwLGHO1zq6z99yxTY/6gZkjZavKgLCgpCbm4uFixYALlcDnd3dyQmJioHOly7dk3lLpq3tze2bduG+fPnY968eXB1dUV8fLxyjjoAmDNnDkpKSjB16lQUFBRg8ODBSExMhJGREYDqO3uXLl3CpUuX0LlzZ5V8WniGFyIiIqJGafGiDgCmT5+u8XHr77//Xqtt1KhRGDVqlMbjiUQiLF68GIsXL1a7fcKECZgwYUJjUiUiIiJqldiRjIiIiKgdYFFHRERE1A6wqCMiIiJqB1jUEREREbUDLOqIiIiI2gEWdURERETtAIs6IiIionaARR0RERFRO8CijoiIiKgdYFFHRERE1A6wqCMiIiJqB1rF2q9ERETUdmQUZtQbY1maB1kz5EJ/Y1FHREREDWIpsYTUQIrwQ+H1xkr1xEjQ12+GrKgGizoiIiJqEJmpDAkBCchX5NcZl1GYgfBD4cjXZy+v5sSijuqVVZyl8h9wZlEx9Ixu4vK9QpSLDYGiK8gQ868SEdHjQGYqg8yUD1ZbI/5LTHXKKs5CQEIASu+XqrSbOAGzLgBwkAEpEQAAqYEUlhLLFsiSiIiIWNRRnfIV+Si9X4rIZyPR3bw7AOByTjFC4tKxwc8UjvtnACO/BDr1gKXEkr+9ERERtRAWddQg3c27w83KDQBQVVaIqrJcOBubw6W8AjDrBvy1jYiIiFoGezASERERtQMs6oiIiIjaARZ1RERERO0A+9QRERFRk8gwNES5JBeZReehZ2SqNoaD7HSHRR0RERHplKXEElI9McJtrAHsxLyjOzXGSg2kSAhIYGGnAyzqiIiISKdkpjIkeC9FftxohJQHY0aQP5xtat+pU648ochnUacDLOqIiIhI52RSa8jKKyBWdIKTWU+4WZm3dErtHgdKEBEREbUDLOqIiIiI2gE+fqW/FVwH7t1WbSu6Uv0z9wKgKAcAGOUVo7coE5IC9SOZiIiIqPmxqKNqBdeBdR5AxT3VdrEh4CADdk4GyisAAC4A9kgA7AdgaAwYWzV3tkRERPQQFnVU7d7t6oJu5CbAusff7UVXgJQIYOSX1Wu8AriUW4yQ2HSsGu0Ol65dAQvHFkmZiIiI/saijlRZ9wDs3f9+LxFX/+zUA7ByAwCUCYU4LRSizLovYMHRTERERK0BB0oQERERtQO8U0dEREQtKqMwo94YLidWPxZ1RERE1KQu5RSrbc8r04dEzwjhh8LrPYZEzwifeX8LayO7WtssTcRwsJA+cp5tXYsXdevWrcPy5cshl8vRr18/rFmzBh4eHhrjd+zYgYiICFy5cgWurq5YunQphg0bptwuCAIWLlyITZs2oaCgAIMGDcKGDRvg6uqqjPn444+xZ88epKenQywWo6CgoCk/IhER0WPLzTALm7bv0rjdymAkqvTLVNqKYIJcWCjf64lzAIc4vPmfA6gqc6h1DKmhPn6b6fPYF3YtWtTFxcUhLCwM0dHR8PT0RFRUFPz8/HD+/HnY2NjUij98+DDGjBmDyMhIDB8+HNu2bUNgYCDS0tLQp08fAMCyZcuwevVqxMTEwMnJCREREfDz88OZM2dgZGQEACgvL8eoUaPg5eWFr776qlk/MxER0WPB2AowNMbyirWARLtdqwykuDhqHypMqwu4zKLzmHc0DquC3OFk1lMl9lJOMULj0pFfUs6iriVP/vnnn2PKlCmYOHEiACA6Ohp79uzB119/jblz59aKX7VqFYYOHYrZs2cDAD766CPs3bsXa9euRXR0NARBQFRUFObPn4+AgAAAwDfffANbW1vEx8dj9OjRAIBFixYBALZs2dIMn5KIiOgxZOEIBB+tPal9ffIuQG/nFPTsUA7YV8+woGdUPdm9s40p15CtQ4sVdeXl5UhNTUV4+N/P0fX09ODr64vk5GS1+yQnJyMsLEylzc/PD/Hx8QCAzMxMyOVy+Pr6Krebm5vD09MTycnJyqKuMRQKBRQKhfJ9UVFRo49FRET0WLBw5FymzajFpjTJy8tDZWUlbG1tVdptbW0hl8vV7iOXy+uMr/mpzTEbKjIyEubm5sqXoyP/khIREVHr0eIDJdqK8PBwlbuERUVFbb6wyyrOQr4iv/pN0ZXqJcGKrvw94TAaNsyciIiIWl6LFXXW1tbQ19dHdna2Snt2djbs7GoPVwYAOzu7OuNrfmZnZ0Mmk6nEuLu7P1K+EokEEomWPT1bsaziLAQkBKD0funfjQ6y6iXBHiI1kMJSYtmM2REREZG2WqyoE4vF6N+/P5KSkhAYGAgAqKqqQlJSEqZPn652Hy8vLyQlJSE0NFTZtnfvXnh5eQEAnJycYGdnh6SkJGURV1RUhJSUFEybNq0pP06bk6/IR+n9UkQ+G4nu5t2B3AvAzsnVa7x26qESywkfiYiIWr8WffwaFhaG8ePHY8CAAfDw8EBUVBRKSkqUo2HHjRsHBwcHREZGAgBCQkLg4+ODFStWwN/fH7GxsTh27Bg2btwIABCJRAgNDcWSJUvg6uqqnNLE3t5eWTgCwLVr13Dnzh1cu3YNlZWVSE9PBwC4uLjA1NS0Wa9BS+tu3h1uVm6AohworwDMuinXeCUiIqK2o0WLuqCgIOTm5mLBggWQy+Vwd3dHYmKicqDDtWvXoKf391gOb29vbNu2DfPnz8e8efPg6uqK+Ph45Rx1ADBnzhyUlJRg6tSpKCgowODBg5GYmKicow4AFixYgJiYGOX7p556CgCwf/9+PP/88038qYmIiIh0r8UHSkyfPl3j49bff/+9VtuoUaMwatQojccTiURYvHgxFi9erDFmy5YtnKOOiIiI2pUWm9KEiIiIiHSHRR0RERFRO8CijoiIiKgdYFFHRERE1A6wqCMiIiJqB1jUEREREbUDLOqIiIiI2gEWdURERETtAIs6IiIionaARR0RERFRO8CijoiIiKgdaPG1X0n3bhaUIr+kvM6YzKJiAMDlnGJUlRXCKK8YLgAu5RajTCisc99LOcW6SpWIiEi9vAt//7noSvXP3AuAQvXfN6O8YvQWZcKw2AmAebOl1xqxqGtnbhaUwnfFAZRWVNYZp2d0EyZOQEhcOqrKctFblIk9EiAkNh2n6ynqAEBqqA9LE7Gu0iYiIqpmbAUYGgM7p/zdJjYEHGTAzslAeYVKuAuAPRKgaocUmP4nYOHYvPm2Iizq2pn8knKUVlQiKsgdLjamGuMyi85j3lFgVZA7nMx6wijPHNgFrBrtjjLrvvWex9JEDAcLqS5TJyIiqi7Kgo8C927/3VZ0BUiJAEZ+CZh1Uwm/lFuMNXF7sArrq/dhUUftjYuNKfo4aL4NrWdUXfA525jCzcocEFW/d+lkCtg/3reviYiohVk4qhZnkr+eDHXqAVi5qYSWCYW4JKQ3X26tGIs6IiIiahMyCjNqtWUWFaNckoszYkOg6AoszWwhM5W1QHYtj0UdERERtWqWEktIDaQIPxSuPsAJCIIMSImANFWKhICEx7KwY1FHRERErZrMVIaEgATkK/JrbbucU4zVcXuwSrwOGX6LEH4qGvmKfBZ1RERERK2RzFSmtlCrKiuEWHEUbqgATBxaILPWg5MPExEREbUDvFPXDokMCpBZdF45wlUddZ1NiYiIqO1iUdfO5JXJYer8GeYdvV9vrFRPDMuCW9Wzcz84czcRERG1OSzq2pnSwouA3n1E5uShe0VFnbGWlVWQXX717wZD4+qZvImIiKjNYVHXzugpigAAJv1C4dbHT7udja0e65m4iYiI2jIWde1UeQdHwN69pdMgIiKiZsLRr0RERETtAIs6IiIionaAj1+JiIiofSi4Wv0z90L1zA4N0Y76k7OoIyIiojYtX+iAKgMpsO8jwEEG7JwMlNc9A4SSoTEQfLRdFHYs6oiIiKhNuwVrXBy1DxAuACkRyPBbVOeSYZaGppBJravnaN05Bbh3m0UdERERUWtQYeoAK3NLSFOlCD8VXWes1ECKhIAEKFeSbcwE/K3wsS2LOiIiImoXZKYyJAQkIF+RrzEmozAD4YfCkZqTiu56JoBxB2D3O2pjLSurIKusVH+gVvjYlkVdG5JVnFXnX1QAuFEmb6ZsiIiIWh+ZqQwyU5nG7ZYSS0gNpAg/FF7dYGupMVaqJ0aC99LqR7UPaqWPbVtFUbdu3TosX74ccrkc/fr1w5o1a+Dh4aExfseOHYiIiMCVK1fg6uqKpUuXYtiwYcrtgiBg4cKF2LRpEwoKCjBo0CBs2LABrq6uypg7d+7g3XffxU8//QQ9PT28+uqrWLVqFUxNTZv0szZWVnEWAhICUHq/tN5YaVUVzAxa5+cgIiJqSQ25mwf8fUcv38IeMiu3Zsru0bR4URcXF4ewsDBER0fD09MTUVFR8PPzw/nz52FjY1Mr/vDhwxgzZgwiIyMxfPhwbNu2DYGBgUhLS0OfPn0AAMuWLcPq1asRExMDJycnREREwM/PD2fOnIGRkREAYOzYscjKysLevXtRUVGBiRMnYurUqdi2bVuzfv46FVyv/i0AQH7RFZTeL0Vkn7fRvY7On9nXzqNXylKUPN2xubIkIiJqU+q7m/egjMKM2o1FVwCxYfVPiRiWEssGH68piQRBEFoyAU9PTzzzzDNYu3YtAKCqqgqOjo549913MXfu3FrxQUFBKCkpwe7du5VtAwcOhLu7O6KjoyEIAuzt7TFz5kzMmjULAFBYWAhbW1ts2bIFo0ePxtmzZ+Hm5oY///wTAwYMAAAkJiZi2LBhuHHjBuzt7evNu6ioCObm5igsLISZmZkuLoWqguvAOg+g4h4A4IzYEEEOMsTdzIJbPcO07wkSXH/td/Ts2TZ+syAiImqsUzcLMXzNH9j97mD0cTDX2XG1ekJWM/CiiQq7htYcLXqnrry8HKmpqQgPD1e26enpwdfXF8nJyWr3SU5ORlhYmEqbn58f4uPjAQCZmZmQy+Xw9fVVbjc3N4enpyeSk5MxevRoJCcnw8LCQlnQAYCvry/09PSQkpKCV155RYefspHu3a4u6EZuAqx7VP82kBIBjPwSMOumcbdLucUY991lbDTVfDePiIiI6lbnY9rcC9Vz4Y38Ehlig+rHtIr8Fr9b16JFXV5eHiorK2Fra6vSbmtri3PnzqndRy6Xq42Xy+XK7TVtdcU8/GjXwMAAHTt2VMY8TKFQQKFQKN8XFhYCqK6em8TdYkAhABJ7wLQ7isvLUVlaiWIjexSZdte4W15RIW4oslB8twhFRaKmyY2IiKiVKL5bhCrFvSb5d88EJjAxNKm9QXQHuFsOiKxRrCeu/vf5bjGKDJumJqipNep7uNriferaisjISCxatKhWu6NjE496+fRZlbee8GzQbl5RTZALERFRK9Ui/+498G90Q/99fhR3796FubnmR8wtWtRZW1tDX18f2dnZKu3Z2dmws7NTu4+dnV2d8TU/s7OzIZPJVGLc3d2VMTk5OSrHuH//Pu7cuaPxvOHh4SqPfauqqnDnzh1YWVlBJOIdMV0qKiqCo6Mjrl+/3jT9FalOvP4th9e+5fDatyxe/7oJgoC7d+/W2+e/RYs6sViM/v37IykpCYGBgQCqi6WkpCRMnz5d7T5eXl5ISkpCaGiosm3v3r3w8vICADg5OcHOzg5JSUnKIq6oqAgpKSmYNm2a8hgFBQVITU1F//79AQD79u1DVVUVPD3VV9oSiQQSiUSlzcLCopGfnBrCzMyM/3G3IF7/lsNr33J47VsWr79mdd2hq9Hij1/DwsIwfvx4DBgwAB4eHoiKikJJSQkmTpwIABg3bhwcHBwQGRkJAAgJCYGPjw9WrFgBf39/xMbG4tixY9i4cSMAQCQSITQ0FEuWLIGrq6tyShN7e3tl4fjEE09g6NChmDJlCqKjo1FRUYHp06dj9OjRDRr5SkRERNTatHhRFxQUhNzcXCxYsAByuRzu7u5ITExUDnS4du0a9PT0lPHe3t7Ytm0b5s+fj3nz5sHV1RXx8fHKOeoAYM6cOSgpKcHUqVNRUFCAwYMHIzExUTlHHQBs3boV06dPxwsvvKCcfHj16tXN98GJiIiIdKjF56kjephCoUBkZCTCw8NrPfKmpsfr33J47VsOr33L4vXXDRZ1RERERO2AXv0hRERERNTasagjIiIiagdY1BERERG1AyzqqMV8+OGHEIlEKq9evXopt5eVlSE4OBhWVlYwNTXFq6++WmviaWqYgwcP4uWXX4a9vT1EIpFyreQagiBgwYIFkMlkkEql8PX1xcWLF1Vi7ty5g7Fjx8LMzAwWFhaYNGkSiouLm/FTtE31XfsJEybU+u9g6NChKjG89o0TGRmJZ555Bh06dICNjQ0CAwNx/vx5lZiG/H/m2rVr8Pf3h7GxMWxsbDB79mzcv3+/OT9Km9SQ6//888/X+vv/9ttvq8Tw+jccizpqUb1790ZWVpby9ccffyi3vffee/jpp5+wY8cOHDhwALdu3cLIkSNbMNu2q6SkBP369cO6devUbl+2bBlWr16N6OhopKSkwMTEBH5+figrK1PGjB07FqdPn8bevXuxe/duHDx4EFOnTm2uj9Bm1XftAWDo0KEq/x189913Ktt57RvnwIEDCA4OxpEjR7B3715UVFTgpZdeQklJiTKmvv/PVFZWwt/fH+Xl5Th8+DBiYmKwZcsWLFiwoCU+UpvSkOsPAFOmTFH5+79s2TLlNl5/LQlELWThwoVCv3791G4rKCgQDA0NhR07dijbzp49KwAQkpOTmynD9gmAsGvXLuX7qqoqwc7OTli+fLmyraCgQJBIJMJ3330nCIIgnDlzRgAg/Pnnn8qYX375RRCJRMLNmzebLfe27uFrLwiCMH78eCEgIEDjPrz2upOTkyMAEA4cOCAIQsP+P/Pzzz8Lenp6glwuV8Zs2LBBMDMzExQKRfN+gDbu4esvCILg4+MjhISEaNyH1187vFNHLerixYuwt7dH9+7dMXbsWFy7dg0AkJqaioqKCvj6+ipje/XqhS5duiA5Obml0m2XMjMzIZfLVa61ubk5PD09ldc6OTkZFhYWGDBggDLG19cXenp6SElJafac25vff/8dNjY26NmzJ6ZNm4bbt28rt/Ha605hYSEAoGPHjvj/9u41pKk+jgP4V7wMK9Rs5maZTDOldFBCY0hSGF4I6fbCDKGLBOYFCotQ6EVCEES9SnrTyoKIfFEIBoI1LYqhGM5L1GCiSbQlaN7SqOHvefHgoT36PJna9uz0/cBA9z/++Z3ftj9fjmfnAItbZ2w2G9LT05UL4gNAbm4uJiYm8ObNGx9WH/j+2f859+/fh1arRVpaGqqrqzE9Pa2Msf+/xu93lKA/l8lkQn19PVJSUuByuXDp0iXs2rULfX19cLvdCAsLm3d/3djYWLjdbv8UrFJz/fxx0Zz7fW7M7XZj/fr1XuMhISGIjo7m67FMeXl5OHToEAwGA/r7+1FTU4P8/HzYbDYEBwez9ytkdnYWZ86cQWZmpnIHosWsM263e8HPxtwYLc5C/QeAo0ePIiEhAXFxcejp6cGFCxfgcDjw6NEjAOz/r2KoI7/Jz89XfjYajTCZTEhISEBDQwPCw8P9WBmR7xw5ckT5OT09HUajEUlJSWhra0N2drYfK1OX8vJy9PX1eZ23S77zb/3/8dzQ9PR06PV6ZGdno7+/H0lJSb4uM+Dx36/0vxEVFYUtW7bA6XRCp9Ph27dvGBsb89rm06dP0Ol0/ilQpeb6+c9v/P3Ya51Oh+HhYa9xj8eD0dFRvh4rLDExEVqtFk6nEwB7vxIqKirQ1NSE1tZWbNy4UXl+MeuMTqdb8LMxN0Y/92/9X4jJZAIAr/c/+794DHX0vzE1NYX+/n7o9XpkZGQgNDQUz549U8YdDgeGhoZgNpv9WKX6GAwG6HQ6r15PTEygvb1d6bXZbMbY2Bhev36tbGO1WjE7O6sswrQyPnz4gJGREej1egDs/XKICCoqKvD48WNYrVYYDAav8cWsM2azGb29vV7BuqWlBREREdi6datvdiRA/az/C7Hb7QDg9f5n/3+Bv7+pQX+uqqoqaWtrk4GBAXn16pXs3btXtFqtDA8Pi4hIaWmpbNq0SaxWq3R2dorZbBaz2eznqgPT5OSkdHV1SVdXlwCQ69evS1dXl7x//15ERK5cuSJRUVHS2NgoPT09sn//fjEYDDIzM6PMkZeXJ9u3b5f29nZ5+fKlJCcnS1FRkb92KWD8V+8nJyfl3LlzYrPZZGBgQJ4+fSo7duyQ5ORk+fr1qzIHe780p0+flsjISGlraxOXy6U8pqenlW1+ts54PB5JS0uTnJwcsdvt0tzcLDExMVJdXe2PXQooP+u/0+mU2tpa6ezslIGBAWlsbJTExETJyspS5mD/fw1DHflNYWGh6PV6CQsLkw0bNkhhYaE4nU5lfGZmRsrKymTt2rWyatUqOXjwoLhcLj9WHLhaW1sFwLzHsWPHROTvy5pcvHhRYmNjRaPRSHZ2tjgcDq85RkZGpKioSNasWSMRERFy4sQJmZyc9MPeBJb/6v309LTk5ORITEyMhIaGSkJCgpw6dcrr8g0i7P1SLdR3AHLnzh1lm8WsM4ODg5Kfny/h4eGi1WqlqqpKvn//7uO9CTw/6//Q0JBkZWVJdHS0aDQa2bx5s5w/f17Gx8e95mH/Fy9IRMR3xwWJiIiI6HfgOXVEREREKsBQR0RERKQCDHVEREREKsBQR0RERKQCDHVEREREKsBQR0RERKQCDHVEREREKsBQR0RERKQCDHVEREREKsBQR0S0Amw2G4KDg7Fv374Fx91uNyorK5GYmAiNRoP4+HgUFBR43Uz++PHjOHDggI8qJiK1YagjIloBFosFlZWVePHiBT5+/Og1Njg4iIyMDFitVly9ehW9vb1obm7Gnj17UF5e7qeKiUhtQvxdABFRoJuamsLDhw/R2dkJt9uN+vp61NTUKONlZWUICgpCR0cHVq9erTy/bds2nDx50h8lE5EK8UgdEdEyNTQ0IDU1FSkpKSguLsbt27chIgCA0dFRNDc3o7y83CvQzYmKivJxtUSkVgx1RETLZLFYUFxcDADIy8vD+Pg4nj9/DgBwOp0QEaSmpvqzRCL6AzDUEREtg8PhQEdHB4qKigAAISEhKCwshMViAQDliB0R0e/Gc+qIiJbBYrHA4/EgLi5OeU5EoNFocOPGDSQnJyMoKAjv3r3zY5VE9CfgkToioiXyeDy4d+8erl27Brvdrjy6u7sRFxeHBw8eIDo6Grm5uairq8OXL1/mzTE2Nub7wolIlXikjohoiZqamvD582eUlJQgMjLSa+zw4cOwWCwoLS1FXV0dMjMzsXPnTtTW1sJoNMLj8aClpQU3b97E27dvlb8bHx+H3W73mmvdunWIj4/3xS4RUQALEp7wQUS0JAUFBZidncWTJ0/mjXV0dMBkMqG7uxtGoxEulwuXL19GU1MTXC4XYmJikJGRgbNnz2L37t0A/r748N27d+fNVVJSglu3bv3u3SGiAMdQR0RERKQCPKeOiIiISAUY6oiIiIhUgKGOiIiISAUY6oiIiIhUgKGOiIiISAUY6oiIiIhUgKGOiIiISAUY6oiIiIhUgKGOiIiISAUY6oiIiIhUgKGOiIiISAUY6oiIiIhU4C+BNbwiRV+gZQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "bins = numpy.linspace(0, 100, 10)\n", "\n", "\n", "plt.figure()\n", "\n", "mask = voxel_dist < 20\n", "plt.hist(voxel_acl[mask], bins=\"auto\", histtype='step', density=1, label=r\"$0 < R / (\\mathrm{Mpc} / h) < 20$\")\n", "\n", "mask = (voxel_dist > 20) & (voxel_dist < 40)\n", "plt.hist(voxel_acl[mask], bins=\"auto\", histtype='step', density=1, label=r\"$20 < R / (\\mathrm{Mpc} / h) < 40$\")\n", "\n", "mask = (voxel_dist > 40) & (voxel_dist < 60)\n", "plt.hist(voxel_acl[mask], bins=\"auto\", histtype='step', density=1, label=r\"$40 < R / (\\mathrm{Mpc} / h) < 60$\")\n", "\n", "# plt.scatter(voxel_dist.flatten(), voxel_acl.flatten(), s=0.1)\n", "plt.legend()\n", "plt.title(\"ACL of individual voxels\")\n", "plt.xlabel(r\"$\\mathrm{ACL}$\")\n", "plt.ylabel(r\"Normalized bin counts\")\n", "\n", "plt.tight_layout()\n", "plt.savefig(\"../plots/BORG_Stephen_ACL.png\", dpi=450)\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "venv_csiborg", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.4" } }, "nbformat": 4, "nbformat_minor": 2 }