From 57b7e949c0b0bb1cd65f469d0b745958923091c6 Mon Sep 17 00:00:00 2001 From: jmframe Date: Sun, 24 Mar 2024 10:22:25 -0700 Subject: [PATCH] number of years --- README.md | 4 +++- figure1.ipynb | 13 ++++++---- figure3.ipynb | 66 +++++++++++++++++++++++++++++++++++++++++++++++---- table1.ipynb | 11 ++++++--- 4 files changed, 80 insertions(+), 14 deletions(-) diff --git a/README.md b/README.md index ef64f77..56c5a84 100644 --- a/README.md +++ b/README.md @@ -25,6 +25,8 @@ To use the analysis notebooks: 4. Activate the environment: * ``conda activate nwm-cnn`` 5. Open the Jupyter notebooks in Jupyter Lab or Notebook +6. Monify the paths to your locally downloaded data, replacing the existing path: +* ``LOC_DATA_DIR = "../2023_nwm_cnn_california/"`` # Floodbase -For more information please visit: (https://www.floodbase.com/about) \ No newline at end of file +For more information please visit: (https://www.floodbase.com/about) diff --git a/figure1.ipynb b/figure1.ipynb index cb46e5c..8c17a47 100755 --- a/figure1.ipynb +++ b/figure1.ipynb @@ -22,7 +22,11 @@ "outputs": [], "source": [ "# Set up variables, etc.\n", - "LOC_DATA_DIR = \"/mnt/d/2023_nwm_cnn_california/\"\n", + "\n", + "#### User needs to modify this line to correspond to their local data, which can be downloaded from Hydroshare:\n", + "#### http://www.hydroshare.org/resource/dbf8e4c2a39a4c228db867b04f9c21ed\n", + "LOC_DATA_DIR = \"../2023_nwm_cnn_california/\"\n", + "\n", "loc_max_daily_nwmcnn = f\"{LOC_DATA_DIR}NWM-CNN_predictions/v01_daily_max/\"\n", "loc_huc_shp = f\"{LOC_DATA_DIR}shapefile/WBD_18_HU2_Shape/Shape/\"\n", "loc_max_daily_huc8 = f\"{LOC_DATA_DIR}csv_files/NWM-CNN_mean_huc8\"\n", @@ -293,7 +297,7 @@ ], "metadata": { "kernelspec": { - "display_name": "Python 3", + "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, @@ -307,9 +311,8 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.8.18" - }, - "orig_nbformat": 4 + "version": "3.9.7" + } }, "nbformat": 4, "nbformat_minor": 2 diff --git a/figure3.ipynb b/figure3.ipynb index c2269a6..88edc53 100644 --- a/figure3.ipynb +++ b/figure3.ipynb @@ -19,7 +19,10 @@ "metadata": {}, "outputs": [], "source": [ - "LOC_DATA_DIR = \"/mnt/d/2023_nwm_cnn_california/\"\n", + "#### User needs to modify this line to correspond to their local data, which can be downloaded from Hydroshare:\n", + "#### http://www.hydroshare.org/resource/dbf8e4c2a39a4c228db867b04f9c21ed\n", + "LOC_DATA_DIR = \"../2023_nwm_cnn_california/\"\n", + "\n", "with open(f\"{LOC_DATA_DIR}csv_files/sac_county_annual_water_year_max_ensemble_water_area_km2.csv\") as f:\n", " df_annualmax = pd.read_csv(f)\n", "with open(f\"{LOC_DATA_DIR}csv_files/ncdc_sac_damage.csv\") as f:\n", @@ -33,9 +36,9 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhwAAAFrCAYAAACaMVCYAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAA9hAAAPYQGoP6dpAACBH0lEQVR4nO3dd1hT59vA8W9YghMUBVkq7j1Q60YcKG7QirV1W+vPveqssypV61bcim0dOGrde7fW2br3RBBFloCCMvL+YclrBCSBhOX9uS6vNuc8ec6dw0ly5znPUCiVSiVCCCGEEHpkkNkBCCGEECLnk4RDCCGEEHonCYcQQggh9E4SDiGEEELonSQcQgghhNA7STiEEEIIoXeScAghhBBC7yThEEIIIYTeScIhhBBCCL0zyuwAPkfh4eGEhYUhk7wKIUT2oVAosLCwwNzcPLNDyZYUMrV5xnr+/DkxMTHkzp0bhUKR2eEIIYTQkFKp5M2bN5iammJtbZ3Z4WQ7ckslg71584Y8efJIsiGEENmMQqEgT548vHnzJrNDyZYk4RBCCCGE3knCIYQQQgi9k4RDCCGEEHono1SyiC+/PIS//+sMO56dXR62bnXNsOOJrGPXrl1ERkZmyLHy5ctHu3btMuRYQncGDx7Mq1ev+OWXXwDo0KEDlSpVYvr06ZkcWeqKFCmCj48PrVq1yuxQxEck4cgi/P1f8+BBRGaHkaLBgwfj6+sLgJGREebm5lSoUAEPDw+6dOmCgYE0ln3KmTNnmDRpEnfu3MHKyopBgwbRs2fPTIklMjKSV69eZcqxNZF4rf3www8MGTJEtX3fvn307NmToKCgTIzu87Ru3TqMjY0zOwyRzcm3hNBYkyZNuHbtGpcuXWLz5s00aNCACRMm8PXXXxMXF5fZ4enNq1ev0tUi8OTJE7p27UqdOnU4evQow4YNY8KECezevVuHUeYspqamLF68mPDw8MwO5ZPi4+NJSEhIsv3du3eZEI3+WFhYkDdv3swOQ2RzknAIjZmYmGBlZUXRokWpUqUKw4YN45dffuHo0aNs3rxZVW7ZsmU4OztTvHhxqlWrxujRo4mKilLt37x5M6VKleLQoUPUrVuXYsWK0bt3b16/fs3mzZtxcnKidOnSjBs3jvj4eNXztm7dSvPmzSlRogQVK1akf//+vHz5Ui3GAwcO8MUXX+Dg4IC7uzubN2+mSJEiar/oz58/T7t27XBwcKBatWqMHz+e169Tvp1148YNKlWqxP/+9z9OnDiR7BfMp6xfvx5bW1umT59OmTJl+Oabb/jqq6/w9vbWqp7PSaNGjShSpAgLFy78ZLndu3fTsGFD7OzscHJy0uicHjhwgObNm2Nvb0+5cuXUWprCw8MZOHAgpUuXplixYnTp0oWHDx+q9n947TZo0AA7OzuePn2Kk5MT8+bNY/DgwZQsWZIRI0YAqV9ra9eu5YsvvsDe3p4KFSrQu3dv1b4OHTowduxYxo4dS6lSpShbtixeXl5qEwZqGu+xY8eoX78+xYsXx9PTkxcvXqjKxMfHM3HiRNUxpk6dmmRSwg4dOvDDDz+oHjs5ObFgwQKGDh1KiRIlqF69uur2S6Lz58/j4uKCvb09zZs3Z9++fRQpUoRr164l+3eZPn06bm5uSbY7Ozsza9YsAP799186depEuXLlKFmyJO3bt+fq1avJ1gfw119/JXn/X7t2jSJFiuDn56cWqzafCSJtJOEQ6dKwYUMqVqzI3r17VdsMDAyYMWMGJ0+eZPHixfz5559MmzZN7XnR0dGsWrWKlStXsnnzZs6cOUOvXr04evQoGzduZOnSpfz6669qrQCxsbGMHTuW48ePs379ep48eaLW5O7n50efPn1wc3Pj2LFjdO/eHS8vL7Xj3rx5E09PT1q3bs2JEydYuXIl586dY9y4cSm+xrp167Jp0yZy5cpFnz59qFGjBjNmzOD+/fsanaOLFy/SuHFjtW0uLi5cuXKF2NhYjer43BgYGDB+/HjWrFnDs2fPki1z5coVvv32Wzp06MDJkyf5/vvvmTVrllry+7HDhw/Tq1cvmjVrxtGjR9m2bRvVqlVT7R8yZAhXrlzh119/Ze/evSiVSr766iu1v1N0dDQLFy5k3rx5nD59GktLSwCWLl1KuXLlOHz4MCNGjEj1Wrt8+TITJkxgzJgxnDlzBl9fX+rWrasWr6+vL0ZGRhw4cIAZM2awYsUKfvvtN63j9fb2ZunSpezatYuAgAAmT56s2u/t7c2mTZuYP38+u3fvJjw8nH379qX6N1q2bBlVq1bl6NGj9OrVi9GjR3Pv3j0AoqKi6NatG+XLl+fIkSOMHTuWH3/88ZP1dezYkUuXLvHo0SPVttu3b3Pr1i06deqkqtfT05Ndu3axf/9+HB0d+eqrr9R+0GgrLZ8JIm2kD4dIt9KlS3Pz5k3V4++++071/8WKFWPs2LGMHj2a2bNnq7bHxsYye/ZsSpQoAUCbNm3Ytm0b169fJ2/evJQtW5b69evz559/0qFDBwC6du2qen7x4sWZOXMmLVq0ICoqirx587J+/XpKlSrFlClTAChVqhS3b99m/vz5quctXboUDw8PVYyOjo7MmDGDDh06MHv2bExNTZO8PoVCQb169ahXrx5eXl7s37+fLVu2sHTpUqpWrUrnzp3p2LEj+fPnT/b8BAUFUbhwYbVthQsXJi4ujtDQUKysrDQ5zZ+d1q1bU7FiRWbPns2CBQuS7F+2bBkNGzZk5MiRAJQsWZI7d+6wdOlSunTpkmyd8+fPp0OHDowZM0a1rVKlSgA8fPiQAwcOsGfPHmrXrq06RvXq1dm/f7+q82tsbCyzZs1SPS9RgwYNGDhwoOrxwIEDP3mt+fv7kzt3blxdXcmbNy/29vZUrlxZrU5bW1t+/PFHFAoFpUqV4tatW6xYsYJu3bppFe+cOXNU77XevXszd+5c1TFWrlzJkCFDaNu2LQBz5szh+PHjKf5dEjVt2lTVIjN48GCWL1/OX3/9RenSpdm2bRsKhYJ58+ZhampK2bJlCQwMVLX8JKd8+fJUqFCB33//XfU33b59O9WrV6dkyZLA+x84H/r5558pXbo0Z86cwdU1bZ3g0/KZINJGEg6RbkqlUm3m1D///JMFCxZw9+5dIiMjiY+PJyYmhtevX5MnTx4AcufOrfoAhPdfwPb29mr3iQsXLkxwcLDq8bVr15gzZw7Xr19XW4smICCAsmXLcv/+fbVfqwDVq1dXe3z16lUePXrE9u3b1bYnJCTg5+dHmTJlPvlazczM8PDwwMPDgwcPHvDdd98xZswY3r17p5ZofezjmWVlRQHNTJo0CQ8PD/73v/8l2Xfv3j1atmyptq127dqsXLmS+Ph4DA0Nkzznxo0bdOvWLdlj3b17FyMjI5ycnFTbChYsSMmSJbl7965qm4mJCRUrVkzy/I+vvdSutcaNG2NnZ0etWrVwcXGhSZMmtGrVity5c6vKOjk5qV07NWvWZNmyZcTHx2sc78fvNSsrK9X7KiIighcvXlCrVi3VfiMjI6pWrZrqNVqhQgXV/ysUCooUKaKq98GDB1SoUEHty/rj92JyOnbsyKZNmxg5ciRKpZIdO3bQr18/1f6XL18ya9Ys/vzzT16+fEl8fDzR0dH4+/unWndK0vuZIDQnCYdIt7t37+Lg4ADA06dP6dq1Kz169GDs2LFYWFhw7tw5hg0bptax1MhI/dJTKBTJbkvsL/H69Ws6d+6Ms7Mz3t7eFCpUCH9/fzw9PVUd9D5OfBK3fSghIYHu3bvTt2/fJK/Dzs4u1dcaFxfHyZMn2bJlCwcOHKBYsWJMnDiRjh07pvicIkWKJBlZERwcjJGREQULFkz1mJ+zunXr4uLiwowZM5K0Wmjy9/7Yp36tpvTcj49jamqa7NIEHyYKkPq1ZmJiwtGjR/nrr784ceIEs2bNYs6cORw6dIgCBQp88nVoE29y7ytdJLwfj1r58P2alr8NvE84pk+fztWrV4mOjiYgIEDVwgnvbyGFhIQwffp01Tls3bp1ircmE0fPfXjsjzu4p/czQWhOEg6RLqdPn+bWrVv0798feH9fOi4ujqlTp6re7Dt37kz3ce7fv09ISAgTJ07E1tZWdawPlS5dmiNHjqhtu3LlitrjypUrc+fOHRwdHbU6/rVr19iyZQu///478fHxuLu7s3PnziS/apNTs2ZNDh06pLbtxIkTVK1aVYYaauCHH36gSZMmqmb1RGXKlOHcuXNq2y5cuEDJkiWTbd2A97/KT506xVdffZVkX9myZYmLi+PSpUuqWxShoaE8fPgwTb9yNbnWjIyMcHZ2xtnZmVGjRlG6dGlOnz5NmzZtALh06ZJa+UuXLuHo6IihoaFO4s2fPz9WVlZcvHhR1X8kLi6Oq1evJrm9o41SpUqxfft23r59S65cuYCk78Xk2NjYULduXbZt20ZMTIyq83Cis2fPMmvWLJo1awa8b90MCQlJsb5ChQoB8OLFC9UKr9evX1crk9bPBKE96TQqNPbu3TtevHhBYGAgV69eZcGCBXTv3h1XV1c6d+4MvO9bERcXx+rVq3n8+DFbtmxh/fr16T62ra0tJiYmqnoPHDjAvHnz1Mp0796d+/fvM23aNB48eMDOnTtVHQgTf20NHjyYixcvMmbMGK5du6a6D/6pDmJnz56lZcuWPHnyhFmzZnH16lW8vLw0SjYAevTogb+/PxMnTuTu3bts3LiRjRs3MmDAgLSdjM9MhQoV6NixI2vWrFHbPmDAAE6fPs3cuXN58OABmzdvZu3atZ88r6NGjWLHjh3MmjWLu3fvcvPmTRYvXgy8v3ffsmVLRo4cydmzZ7l+/ToDBgzA2to6ya0bTaR2rR06dIhVq1Zx7do1nj59ypYtW0hISKBUqVKqOgICApg4cSL379/n999/Z/Xq1Xz77bc6jffbb79l8eLF7N27l3v37jFmzJh0z9PSsWNHEhISGDlyJHfv3uXYsWOqEUSpLVzZsWNH/vjjD3bv3q3qLJqoRIkSbN26lbt373Lp0iX+97//YWZmlmJdJUqUwNbWljlz5vDgwQMOHz7MsmXL1Mqk5TNBpI20cGQRdnZ5svzxjh07RuXKldUm/po5cyaenp6q1ozKlSszbdo0Fi9ezIwZM6hTpw4TJkxg0KBB6YrX0tKSRYsWMXPmTFavXk3lypWZMmWK2v34YsWKsWbNGiZPnsyqVauoWbMmw4YNY/To0ZiYmABQsWJF/vjjD2bOnEm7du1QKpUUL15crdn2Y2XKlOHKlSuq0QjaKlasGBs3bmTixImsW7cOa2trZsyYoeqkl9Hy5cuX7Y41duxYdu3apbatSpUqrFq1itmzZzNv3jysrKwYPXp0ih1GAerXr8/q1auZN28eixcvJl++fNSpU0e1f9GiRUyYMIFvvvmG2NhY6tSpw6ZNm9LUEpXatZY/f3727t3LnDlzePv2LSVKlGDFihWUK1dOVUfnzp2JiYmhRYsWGBoa0rdvX7p3767TeAcMGEBQUBBDhgzBwMCAr776ilatWhERkfaJCPPly8evv/7K6NGjadKkCeXLl2fkyJH0798/1U6Y7dq1Y/z48RgaGiYZJrtw4UJGjhxJ06ZNsbW1ZcKECapO4skxNjZm+fLljB49GhcXF6pVq8a4cePo06ePqkxaPhNE2iiU0nstQz18+DBDP/A/d/Pnz2f9+vVJbr8IkdVlp+nENbFt2zaGDh3K/fv3P9kqkR1ERkbKLZg0kBYOkaOsXbuW6tWrY2Fhwfnz51m6dKnarxkhRMbw9fWlWLFiFC1alBs3bvDjjz/Srl27bJ9siLSThEPkKI8ePWL+/PmEh4dja2vL//73P4YOHZrZYQnx2QkKCmL27NkEBQVhZWVFu3btpF/EZ05uqWQwuaUihBDZm9xSSRsZpSKEEEIIvZOEQwghhBB6JwmHEEIIIfROEg4hhBBC6J0kHEIIIYTQO0k4hBBCCKF3Mg9HFnEtIo63CRl3vFwGUDm/5n/+wYMH4+vrC6A2tbmHhwddunRRTW0uknfz5k3GjRvHv//+i7m5Od27d2fkyJGprivxOZJrTYicSRKOLOJtAkRnYMKRFk2aNGHhwoUkJCTw8uVLjh07xoQJE9i9eze//vprkmWwxXuRkZF8+eWXNGjQgIMHD/LgwQOGDBlC7ty5s8UCboYPH2K6cSMGT5+SYG9PTNeuxOt5DgK51oTIeeSngtCYiYkJVlZWFC1alCpVqjBs2DB++eUXjh49qlqVFWDZsmU4OztTvHhxqlWrxujRo4mKilLt37x5M6VKleLQoUPUrVuXYsWK0bt3b16/fs3mzZtxcnKidOnSjBs3jvj4eNXztm7dSvPmzSlRogQVK1akf//+vHz5Ui3GAwcO8MUXX+Dg4IC7uzubN2+mSJEiaqtfnj9/nnbt2uHg4EC1atUYP348r1+/1tt527ZtG2/fvmXRokWUL1+eNm3aMGzYMJYvX05Wn3cv16ZNWNSrh9nSpeTauROzpUuxqFePXB/8vfVBk2tNrjMhshdJOES6NGzYkIoVK7J3717VNgMDA2bMmMHJkydZvHgxf/75J9OmTVN7XnR0NKtWrWLlypVs3ryZM2fO0KtXL44ePcrGjRtZunQpv/76K7t371Y9JzY2lrFjx3L8+HHWr1/PkydPGDJkiGq/n58fffr0wc3NjWPHjtG9e3e8vLzUjnvz5k08PT1p3bo1J06cYOXKlZw7dy7V5emLFy/+yX8LFixI8fkXL16kXr165MqVS7XNxcWF58+f4+fnl+o5ziyGDx+Sb/hwFAkJKOLjVf8lIYF8w4Zh8PBhhsbz8bWW064zIXI6aZcU6Va6dGlu3rypevzdd9+p/r9YsWKMHTuW0aNHM3v2bNX22NhYZs+eTYkSJQBo06YN27Zt4/r16+TNm5eyZctSv359/vzzT9Uy0V27dlU9v3jx4sycOZMWLVoQFRVF3rx5Wb9+PaVKlVItV12qVClu377N/PnzVc9bunQpHh4eqhgdHR2ZMWMGHTp0YPbs2ckunV21alWOHTv2yXNgYWGR4r6goCAcHBzUthUuXFi1r1ixYp+sO7OYbtwIyfQxUQBKhQKzjRt5/cMPGRrTh9daTrvOhMjpJOEQ6aZUKtU6P/75558sWLCAu3fvEhkZSXx8PDExMbx+/Zo8efIAkDt3btWXALz/Ara3tydv3rxq24KDg1WPr127xpw5c7h+/TphYWGq2xEBAQGULVuW+/fvU61aNbXYqlevrvb46tWrPHr0iO3bt6ttT0hIwM/PjzJlyiR5fWZmZjpfNyEx9qzcadTg6VNI6ZaPUvl+fwb78FrLadeZEDmdJBwi3e7evav6Bf/06VO6du1Kjx49GDt2LBYWFpw7d45hw4YRFxenes7Hnf4UCkWy2xIS3vekff36NZ07d8bZ2Rlvb28KFSqEv78/np6evHv3Dkia+CRu+1BCQgLdu3enb9++SV6HnZ1dsq/v7NmzdOnS5ZPnYNiwYQwbNizZfUWKFEnSByDxCy6xpSMrSrC3T7aFAwCF4v3+DJZ4reXE60yInE4SDpEup0+f5tatW/Tv3x+Ay5cvExcXx9SpU1XDF3fu3Jnu49y/f5+QkBAmTpyIra2t6lgfKl26NEeOHFHbduXKFbXHlStX5s6dO1q1WKT3lkrNmjWZOXMm7969w8TEBIDjx49jbW2d5FZLVhLTtStmS5ag5P1tlERKAKWS6A9uPWSED6+1nHidCZHTSadRobF3797x4sULAgMDuXr1KgsWLKB79+64urrSuXNn4P0977i4OFavXs3jx4/ZsmUL69evT/exbW1tMTExUdV74MAB5s2bp1ame/fu3L9/n2nTpvHgwQN27typGtGQ+It08ODBXLx4kTFjxnDt2jUePnzIgQMHPtmZL/GWyqf+fSrh6NixIyYmJgwZMoRbt26xd+9eFi5cSP/+/bP0LZV4R0ciFywAAwOUhoYo//svBgZELlhAgh6/TFO71nLidSZETicJh9DYsWPHqFy5Mk5OTnTp0oU///yTmTNn8ssvv2BoaAi8/2U3bdo0Fi9ejLOzM9u3b2fChAnpPralpSWLFi1i9+7dNGzYkEWLFqk67SUqVqwYa9asYe/evTRu3BgfHx/VbY7EloWKFSvyxx9/8PDhQ9q1a0eTJk346aefsLKySneMKcmfPz9bt27l2bNnuLq6MmbMGPr378///vc/vR1TV9526ULomTNEDxzI2/btiR44kNAzZ3ibyi2m9ErtWpPrTIjsR6HM6hMB5DAPHz4kX758SbZn9ZlGs6v58+ezfv36JM3iQuiSXGefl8jISLldlgY5/xsnm/gcvvwzwtq1a6levToWFhacP3+epUuX0qdPn8wOS+Qwcp0JoT35lhM5yqNHj5g/fz7h4eHY2tryv//9j6FDh2Z2WCKHketMCO2l65bKy5cviY6OTrI9K/e8z2wp3VIRQgiRPcgtlbRJUwvH9OnTWbRoESEhIcnu/3BdAiGEEEIIrUeprF27lp9++okhQ4agVCoZP34848aNw87OjtKlS7N69Wp9xCmEEEKIbEzrhGPp0qWqJAPA3d2d6dOnc/v2bfLly6c2RbBIKivPuyCEECJ18jmeNlonHPfv36dOnTqq2f0Sp/s1MzNj5MiRrFy5UrcR5jAfTqMshBAie0lISJCEI420TjgS1yFQKBTkz58ff39/1T5LS0sCAgJ0F10OZGVlRWRkpCQdQgiRzSQkJBAZGSkTuKWR1p1GS5cuzdP/VomsVasWq1aton379hgYGLBy5UqKFy+u6xhzFDMzM2xtbXnx4gVKpTLJok9CCCGyHoVCgUKhwNbWFjMzs8wOJ1vSOuFo1aoVp06dokePHowbN44WLVpgbm6OkZERUVFRrF27Vh9x5ihmZmaSmAkhhPispHtq8/Pnz+Pr64tCoaB169a4uLjoKjYhhBBC5BCylooQQggh9C7NU5vfvn2bkydPEhwcTJ8+fbC2tubZs2dYWFjI/S0hhBBCqNG6hSM+Pp5+/frh4+ODUqlEoVBw4cIFatSoQdu2balevTrTpk3TV7xCCCGEyIa0HhY7Y8YMNm7cyJw5c7h+/braKAs3NzcOHDig0wCFEEIIkf1pfUvFx8eHiRMnMmLEiCRrppQoUYJHjx7pLDghhBBC5Axat3AEBARQt27dZPeZmpoSGRmZ7qCEEEIIkbNonXAUKVKEhw8fJrvvzp072NnZpTsoIYQQQuQsWiccrVq1YsaMGWpTmCsUCl69esWiRYto27atTgMUQgghRPan9SiVFy9eUKtWLV69eoWLiwu7d+/G1dWV69evY2xszMWLFylYsKC+4hVCCCFENqR1C4eVlRUXLlzgq6++4tKlSxgaGnLlyhXc3Nw4c+aMJBtCCCHER3799VeGDx/+WY/k1KqFIyYmhmnTptGxY0ecnJz0GVe2EhUVRWxsbGaHkSoLCwvCwsIyOwyNZJdYs0ucILHqi8SqHzklVmNjY/LmzZvBEWVNWg2LNTU1Zf78+bRs2VJf8WRLsbGxREdHZ3YYn6RQKID3SWNWn80+u8SaXeIEiVVfJFb9kFhzJq1vqZQvX17m2hBCCCGEVrROOCZOnMj06dN58OCBPuIRQgghcpy5c+d+9n04tJ5pdN26dbx584by5ctTpUoVihYtqmpSgvfNSzt37tRpkEIIIUR2NnLkyMwOIdNpnXBcvXoVExMTbG1tCQkJISQkRB9xCSGEEFma4cOH5PH1heBg8lla8trTk3hHx8wOK8vSOuF4/PixHsIQQgghsg8zX1/MR40ChQKUSvIoFOTx9ib855+J9vTM7PCyJK37cHxKfHw8f/zxhy6rFEIIIbIUw4cPMR81CkVCAor4ePjgv+ajRmEoAyuSpZOE4/bt24wePRpbW1s6duyoiyqFEEKILCm3r+/7lo2PKAAUCnJv3pxkn3QaTcMtlUSvX7/G19eXNWvWcPbsWZRKJTVq1GDatGm6jE8IIYTIUgyfPoWU5txQKt/v/4h0Gk1DwnHmzBnWrFnD1q1bef36Nblz5wbgt99+o2vXrjoPUKTdqlWrUtz37bffZmAkQojPna3txhT3BQRkr++OeHv7ZFs4AFAo3u8XSWh0S+XFixfMmTOH8uXL07BhQ3x8fHBycsLHx4c7d+6gVCplWXohhBCfhTeenqBU8nEbhxJAqeRNly6ZEFXWp1ELh4ODA3Fxcdja2jJu3Dh69+6N439Df169eqXXAIUQyctJvxiFyE7iHR0J//lnzEeNQqlQoFAqUf43WiX855+JL1Eis0PMkjRKOBIXJrO0tMTGxoZChQrpNSghhBAiK4v29ORd7drk2byZvMHBvLa05HWXLpJsfIJGCcfVq1dZtWoVGzduZNCgQYwaNQoPDw969+5N9erV9R1jjpRS/wrpWyGEENlDfIkSRI4fT96iRYkMDPzk4m1z587F39+fFi1afLYLoGqUcFSqVImFCxcyZ84cduzYwZo1a9i8eTObNm3CxsYGhUJBRESEvmMVQmhJbrsIkTXIKBUt5+EwMTHB09OTQ4cO8fDhQyZOnIiRkRFKpRJ3d3fatm3Lnj179BWrEEIIIbKpNE/85eDgwJQpU3j48CEHDx7Ew8ODI0eO0L59e13GJ4QQQogcIM0TfyVSKBQ0b96c5s2bExoaym+//aaLuIQQQgiRg+h0LZWCBQsyZMgQXVYphBBCCC3Nnz+fihUrUqFCBYYMGfLJDq0ZJd0tHELIjKZCCJF1vHz5kiVLlnDjxg2MjY1p1KgRZ8+epW7dupkalyQcQgghRA4TFxdHTEwM8H4urSJFimRyRDq+pSKEEEKI9Dl16hRt27ZVTTvxxx9/JCnj7e1NiRIlMDU1xcnJidOnT6v2FS5cmFGjRuHg4ICNjQ3NmjWjZMmSGfgKkqd1C8e7d+8wNjZGkdLCNUIIIUQOpu/5bV6/fk3VqlXp1asXHTt2TLLf19eXYcOG4e3tTf369VmxYgVubm7cvHkTBwcHwsLC2LNnD48fP8bMzAw3NzdOnTpFo0aN0h1bemjVwhETE4OZmVmy2ZYQQgghkhcdHc2bN29U/xKXDEmOm5sb06dPx8PDI9n98+bNo0+fPvTt25fy5cuzYMEC7O3tWbZsGQBHjhyhVKlSFCxYEDMzM1q3bs3Zs2f18rq0oVULh6mpKYUKFSJPnjz6ikcIIYTIcaZMmcKjR49Ujzt16kTnzp21rufdu3dcunSJsWPHqm13dXXlzJkzANjb23PmzBliYmIwNjbmxIkT9OvXL30vQAe0vqXStm1bduzYgaurqz7iEUIIIXKcKVOmqA1NNTY2TlM9wcHBxMfHY2VlpbbdysqK58+fA1CnTh1atWpF9erVMTAwoGnTprRr1y7tweuI1glHly5d6NOnD71798bDw4OiRYsm6c9Ro0YNnQUohBBCZHfe3t46Xbzt4+9dpVKptm3GjBnMmDEj3cfRJa0TjhYtWgDg4+PD+vXr1fYlvuD4+HjdRCeEEELkALpavM3S0hJDQ0NVa0aioKCgJK0eWY3WCce6dev0EYcQQgiRJZyNTHkUZp18aZuxU1fL05uYmODk5MThw4dxd3dXbT98+HCWX8tM64SjR48e+ohDCCGEyLG0aeGIiori/v37qsePHj3i8uXLFCxYEAcHB0aMGEG3bt2oWbMmdevWZeXKlfj5+dG/f399hK4z6Zpp9M6dOwQHB1OtWjUZuZIJZEpxIYTIHrRp4bh48SIuLi6qxyNGjADe/+D38fHB09OTkJAQpk2bRmBgIJUqVWLfvn0UK1ZMr68hvdKUcPzyyy+MHz+ewMBAAC5cuECNGjXo3LkzzZs3ly87IYQQ4gPatHA0btw41cXWBgwYwIABA9IbVobSemrzrVu30rNnT2rUqMGSJUvUTkqNGjXYsmWLTgMUQgghsru5c+cyfPhwDhw4kNmhZBqtWzi8vLzo1asXa9asIT4+noEDB6r2lS9fnsWLF+s0QE3s2LGD8+fPExAQgImJCWXKlOGbb77Bxsbmk8+7efMm69evx9/fHwsLC9q1ayfziwghhNA5XY1Syc60Tjhu3brFrFmzkt1XsGBBQkJC0h2Utm7evEmLFi0oWbIk8fHxbN68menTpzNv3jxMTU2TfU5QUBBeXl40bdqUwYMHc+fOHVavXk3+/PmpU6dOBr8CkZ3oou+M9L8RQnxutE44cufOzatXr5LdFxAQgIWFRbqD0taECRPUHg8YMIC+ffvy8OFDKlSokOxzDh06hKWlJT179gTAzs6OBw8esHv3bkk4hBBC6JSuhsVmZ1onHPXr12fJkiXJrmDn4+ND48aNdRFXurx58waAvHnzpljm3r17VKlSRW1btWrVOH78OHFxcRgZpWsAjxBCCKEit1TSkHBMmjSJBg0aULt2bbp27YpCoeD3339n8uTJnDp1ivPnz+sjTo0plUrWr19PuXLlcHBwSLFceHg4BQoUUNtWoEAB4uPjiYyMTLalJjY2Vm2FPwMDA9Utm4+nmU2rjKpHV8dJbxyfek5GxahLuog5o+rIqudX339/G5sNyW5/9uxrrevKTtdqVo7145iyU6za7v+caZ1w1KxZk/379zNgwABVxjZz5kxKly7Nvn37qFSpks6D1MaaNWvw8/Nj2rRpqZZNbi765LYn2rFjB9u2bVM9rl+/PkOHDtXpbaSiRYtmSD26Ok564/gUa2trHUaSMXRxXjOqjoy6BtIqo//+n8u1mhVjTencZ1qsEYEp7soqn63ZUZruG7i4uHDr1i0ePHjAixcvsLS0pEyZMrqOTWtr167l0qVLTJ06lUKFCn2yrLm5OeHh4WrbIiIiMDQ0TPFWjLu7O23atFE9NjB4P6o4LCyMmJiY9AX/n8S5TfRdj66Ok944kqNQKLC2tub58+epjkXPanRxXjOqjoy6BrSVWX//nH6tZuVYPz732SnW1PabmppiYWEhfThIQ8Ixbdo0+vbti42NDSVLlqRkyZKqfYGBgaxatYpJkybpNMjUKJVK1q5dy/nz55kyZQpFihRJ9TmlS5fm0qVLatuuXLmCo6Njiv03jI2NU1xSWFdvioyqJ6PexOk5jlKpzHIfNqnRRbyJdaRnPQdN4sjq5zaj//6fy7WaFWNNKZ7MizXl915aP1ulD0caJv6aOnUq/v7+ye579uwZU6dOTXdQ2lqzZg2nT59m6NChmJmZER4eTnh4OO/evVOV2bhxI0uWLFE9dnV1JTg4WDUPx7Fjxzh27Bht27bN8PiFEEKInE7rFo5PZXdRUVEptgDo06FDhwCYMmWK2vYBAwaoRs2EhYURHBys2lekSBHGjRvH+vXrOXjwIBYWFvTq1UuGxArxEVvbjSnuCwjomoGRiMwi14DQBY0SjqtXr3L58mXV43379nH79m21MtHR0WzYsEHtFktG0WQ69Q9nRE1UoUKFFCcxE0IIkXFyelIjfTg0TDh27NihulWiUChSHAFiZmbGunXrdBedEEIIkQNIHw4NE45+/frRpk0blEoltWvXZt26dUmGv+bKlYuSJUtiZmaml0CFEEIIkX1plHAULVpUNbb4+PHjODk5fXIWTyGEEOmTnlFKQmRFWncadXZ21kccQgghhMjB0jTx171791ixYgW3bt0iOjpabZ9CoeDo0aM6CU4IIYRIlNM7luZ0Wicc169fp06dOtja2nL//n2qVKlCcHAwAQEB2NvbZ8oolaxGlh4XQgiRHT169Ih9+/bx119/ERAQQHR0NJaWllSoUIEmTZrQvHnzNE9/ofXEX+PHj6dFixbcuHEDpVLJmjVrePr0Kbt37yYmJobp06enKRAhhBAip5o7dy7Dhw/nwIEDmR1Ksk6cOEHLli0pXbo0gwcP5vTp06q5tR49esTy5ctp06YNdnZ2TJo0iYiICK2PoXULxz///IO3t7dqHZGEhAQAWrduzahRoxg3bhwnT57UOhAhsgJpnRI5jVrnU9WiZArpeJrBsvKwWHd3d/bu3UvLli3ZtGkTjRs3pnDhwmpl4uPjuXr1Kjt27OC3335jxYoVbNiwgWbNmml8HK0TjrCwMAoWLIiBgQHGxsaEhYWp9tWsWVOjVVqFEEIIkTXky5eP27dv4+jomGIZQ0NDqlevTvXq1ZkyZQq//vorAQEBWh1H64TD1tZWNUV4qVKlOHXqFM2bNwfez0gqw2WFEEKI7OOXX37RqryBgQE9evTQ+jhaJxwNGjTgzJkzdOjQga+//prJkycTGBiIiYkJPj4+fPPNN1oHIYTQn1K8pDcXKE4Yj7FgLbW4T+HUnyiEEDqkdcIxYcIEnj17BsCYMWN4/vw5GzZsQKFQ0LlzZ37++WedB/k5kr4EQhd6coHVbEWJAgVKlCgYzQn68CWQccMIZRIrITLOnTt38PT0VHu8adMmOnTo8MnnxcbG4uvry6lTp1Tf8zY2Nri4uPDll19iZJSmmTRUNHr2kiVLcHFxoWLFipQsWVI19NXQ0JBFixaxaNGidAUhhNA9w4cPWc1WDFECiV/qSpTAGrYS/GgI8SVKZGKE2VdKHTFBEiiR+cqWLatacDUqKorixYuruj6k5OrVq7Rv3x4/P78kq8KvWbOGiRMnsnv3bsqXL5/muDQaFjtkyBCqVKmClZUVnp6erFixgjt37qT5oEII/cvt64uSpC0LCkCJgtybN2d8UEKIDLVr1y6aNm1Knjx5UiwTFRVF27ZtCQ8PZ9asWdy+fZvo6GjevHnDrVu38PLyIiQkhNatW/P69es0x6JRwnHq1CmmTJlC5cqV2bNnD//73/+oUKECNjY2fP3116xevZoHDx6kOQghhO4ZPn2KgpR+bSsxfPo0Q+MRQmjm1KlTtG3bFhsbGxQKBX/88UeSMt7e3pQoUQJTU1OcnJw4ffp0snVt2bJF7fZKcnx8fHj58iWnTp1i1KhRlClThly5cmFqakrZsmUZPXo0J06c4Pnz5/j4+KT5dWmUcDRo0ICJEydy5MgRwsPDOXnyJJMnT6ZcuXLs2LGDfv36UaZMGRwcHNLUc1UIoXvx9vbJtnC8pyDe3j5D4xFCaOb169dUrVqVJUuWJLvf19eXYcOGMWHCBP79918aNmyIm5sbfn5+auUiIiL466+/aNWq1SePt2fPHrp160blypVTLFO1alW6devGnj17tH9B/9F6plFjY2MaNmzIpEmTOHbsGOHh4Rw/fhx3d3cCAgL47bff0hyMEEJ33nh6/tdRVJ0SUKDkTZcumRGWEJ+lxFsUif9iY2NTLOvm5sb06dPx8PBIdv+8efPo06cPffv2pXz58ixYsAB7e3uWLVumVm7nzp20aNECU1PTT8Z248YNmjZtmupraNq0KdeuXUu1XErS3OX01q1bHD9+nBMnTnDixAmCg4MpVKgQDRs2THMwQgjdiXd0pA9fsoatJPw3SoX//tuHL5mZzTqMajLSRUbDiKxqypQpPHr0SPW4U6dOdO7cWet63r17x6VLlxg7dqzadldXV86cOaO2bcuWLfTr1y/VOkNDQ7GxsUm1nI2NDaGhodoF/AGNE467d+9y/Phxjh8/zsmTJ3nx4gVFihShUaNGTJ48GWdnZypVqpTmQIQQureeWvxJCfpwXjUPxxpq8wBLZurwOPJFL8SnTZkyRW30R1oXQAsODiY+Ph4rKyu17VZWVjx//lz1+NWrV5w/f57t27enWmd0dDQmJiapljM2Nubt27faB/0fjRIOW1tbnj9/TtGiRWnUqBFTpkzB2dmZcuXKpfnAImuQ+T5yvgdYMp5P38MVQuiXt7c3/v7+tGjRgpYtW6a7PoVCPclXKpVq2woUKMCLFy80ri8iIiLV1otXr15pF+RHNEo4AgMDyZMnDy4uLjRu3BhnZ2dKlSqVrgMLIYQQQjuWlpYYGhqqtWYABAUFJWn10EaLFi3SG1qqNEo4/vnnH1VfjdGjRxMeHk7RokVxdnZW/Stbtqy+YxVC5GAymZbIyXS1WqyJiQlOTk4cPnwYd3d31fbDhw/Tvn37NNU5efJkncSWGo0SjmrVqlGtWjWGDRuGUqnk33//5eTJk5w4cYKxY8cSHh6OlZUVzs7ONG7cmP79++s7biGEyNZyUr8XW9uNKe4LCMi4KfRziqioKO7fv696/OjRIy5fvkzBggVxcHBgxIgRdOvWjZo1a1K3bl1WrlyJn59fmr97s1TC8SGFQkGNGjWoUaMGw4cPR6lUcv78eebMmcPWrVvZtm2bJBxCCCHEB+bOnatxH46LFy/i4uKiejxixAgAevTogY+PD56enoSEhDBt2jQCAwOpVKkS+/bto1ixYnp9DemVpmGxCQkJXLhwgRMnTnD8+HHOnDlDVFQUQLruIQkhhBA5kTa3VBo3bpxkPZOPDRgwgAEDBqQ3LAACAgJ4+PBhkmktTp8+zffff8/169extbVlwoQJdO/ePc3H0SjhSEhI4OLFi6p+HH/++SevX79GqVRiaWlJixYtcHFxwcXFJV0LuwghhBA5kTYtHBlt6tSpXLx4kX/++Ue17cmTJ7i5uRETE0OVKlV4+vQpvXr1wtraGldX1zQdR6OEw9zcXJVgWFhY0LRpU1WC8ampUIUQ+leKl/TmgmqejbXUyuyQsjw5ZyKj6arTqD6cPXuWrl3V+9osWrSI6OhofH196dSpE1FRUbi4uLBw4UL9JhyNGjWiSZMmuLi4UK1atSTjf4XILnLavCM9ucBqtqL8bwZRJQpGc4IIX2uiU1mw6XMl50wIdQEBAVSoUEFt24EDByhVqhSdOnUCIG/evAwcODDJDKfa0CjhSM9iLUII/SjFS1azFUOUoFox5f3aKeajRvGudu1U68hJIyU0YfjwYarnLD6bTfkusoesfEslOjqaAgUKqB5HRkZy+/ZtevfurVbO0dExY6Y2F0JkLb258N9qsOqJgQJQKhTk3rwZyNq91jNabl/fVM9Z5LhxmRJbWnxuCWN2lpVvqdjb23Pnzh2cnZ0B+Pvvv1EqldSsWVOt3MeJiba0Xi1WH6Kjo3n06BE3b94kKCgos8MRIlsoTth/C7IlQ6nE8OnTjA0oGzB8+lTOmRAfadq0KXPnzsXPz4/o6GjmzZuHoaFhkmXtL1++jL29fZqPk2ktHAEBAaxatYq9e/dy+fJlEhISVPsKFSqEs7Mz33zzDW3btsXAIEvkRUJkKY+xSPbXOgAKBfHp+GDIqeLt7eWciUyRlW+pjB8/nq1bt1KiRAkMDAyIj4+nf//+SZILX19fGjRokObjZHjCERgYyPjx49mwYQN58uShXr16jB07liJFimBqakpoaCgPHz7k7NmzuLu7U6xYMby8vOjSpUtGhypElraWWozmxH+Lzv8/JYBSyZsuXWDJ35kTXBb1xtOT3EuWfvqcCaEHWfmWip2dHZcvX2blypWEhoZSt27dJKNWnj9/Tv369enWrVuaj5PhCUeZMmWoXbs2mzdvpm3btp9covfhw4esW7eOgQMHEhAQkKX/YEJktPsUpg9fsoatJPw34oL//hvx88//dX6UhOND8Y6OGpwzIT4/tra2TJ06NcX91tbWLF68OF3HSHPCcfv2bU6ePElwcDB9+vTB2tqaZ8+eYWFhgZmZWYrP27lzJ02aNNHoGI6Ojvz444+MGjWKR48efbLszZs32bVrF48ePSIsLIxRo0ZR+xO99G/cuJHsyZ0/fz62trYaxSdEZltPLf6kBH04r5pTYg21OSXDO1Mk50yIzKF1whEfH0+/fv3w8fFBqVSiUChwc3PD2tqa7777jurVqzNt2rQUn69psvGhAgUKUK1atU+Wefv2LcWLF8fFxYW5c+dqXPeCBQvInTu36nH+/Pm1jk+IzPQAS8bTKvWCQkXOmRD/71Pfy8bGxlhZWdG8eXO++uorjIzSfmNE696YM2bMYOPGjcyZM4fr16+rzffu5ubGgQMHNK6rSZMm3L59O9l9d+/e1So5qV69Ol26dOGLL77Q+DnwPpkxNzdX/ZMOqkKI9DobqUjxn/g8zZ07l+HDh2v1HZlRgoKCePnyZbL//Pz82LVrFz169KBBgwa8efMmzcfROlXx8fFh4sSJjBgxgvj4eLV9JUqUSPXWx4dOnDhBREREsvsiIyM5efKktuFpbfTo0cTGxmJnZ4eHhweVKlXS+zGFEEJ8XrJyH8Tr16+nWmb//v14enoyY8YMZsyYkabjaJ1wBAQEULdu3WT3mZqaEhkZmaZAPhYYGKh2q0PXLCws6NevH46OjsTFxXHq1Cl+/PFHJk+enGSK10SxsbHExsaqHhsYGGBqagqg8XTvqZXTpB5dlNHVcfRRR+JzMnoK/ex0zjIi1ux0LWZUHZrIKrHK69WPjDqvWY2bmxtjxozh119/zbiEo0iRIjx8+BAXF5ck++7cuYOdnd0nn79z50527typevzjjz9SuHBhtTLR0dGcOHGC6tWraxuexmxsbLCxsVE9LlOmDMHBwezevTvFhGPHjh1s27ZN9bh+/foMHToUCwsLjY9btGjRdO3XVRldHUfTOlLq/Tx58uQUn2ttbZ3u42sjq50zIgLTfByNY03HMdTKpFaPJsfJiDo+Qau/XXZ6vSmU0cW1qkk92enzSiWjrqNspnbt2kyfPj3Nz9c64WjVqhUzZsygZcuWqi8EhULBq1evWLRoEW3btv3k82/evMnWrVtVzzt27FiSfhO5cuWicuXKLFy4UNvw0qVMmTKcPn06xf3u7u60adNG9Tgx7rCwMGJiYjQ6RmBgyheyJvt1VUbTOlauXJni/n79+uklDoVCgbW1Nc+fP1frI6RvGfW30UUdGRFrVrsWs0Idmsgqscrr1Q9tYzU1NdXqR2lWFh0dna5Oo1o/c9q0aezfv58KFSrg4uKCQqFg/PjxXL9+HWNjYyZOnPjJ548bN45x/61VYGBgwPHjxz85fDUjPXr0CHNz8xT3GxsbpzhviKZfjKmV06QeXZTJDnUolcoMTTgy6m+jeR0pN83qLta0H0O9TGr1aHKcjKgjZdr97bLT602+jK7eW5n9OfFxmZQ65mq3tkxGXUfZy759+1K8A6AJrRMOKysrLly4wOTJk9m7dy+GhoZcuXKFNm3aMG3aNAoWLKhxXR9OZ55eMTExPH/+XPU4KCiIx48fkzdvXiwtLdm4cSOhoaEMGjQIgL1791K4cGHs7e2Ji4vj9OnTnDt3Lkt37BFCCCF07VMrwMbHxxMQEICvry+rV6/G29s7zcdJU9uIlZUVy5cvT/NBk/Py5Uuio6OTbHdwcNDo+Q8ePFDrJ/DLL78A4OzszMCBAwkLCyM4OFi1Py4ujl9//ZXQ0FBMTEywt7dn7Nix1KhRI52vRAhYtWpVstu//fbbDI5ECJEVZOW1VCwtLTXq7Dpo0CCNbqWnROuEIzY2lnfv3pEnT54k+16/fo2Jicknpyv/UEREBCNGjGDTpk0p9oH4eOhtSipWrMiWLVtS3D9w4EC1x+3bt6d9+/Ya1S2EEEKkR1ZuPZ80aVKKCYeRkRFWVlY0adIER0fHdB1H64Tj22+/5e3bt2zatCnJvn79+mFmZsbq1as1qmv48OFs3LiRPn36UKVKFXLlyqVtOEIIIYTGPjX5mnb9PHKOKVOmZMhxtE44jh8/zk8//ZTsvrZt26o6hGpi7969/PTTTwwdOlTbMIRIIqXbGCC3MoQQIrNpPY/3ixcvUhxnnDiUUVMxMTFUrlxZ2xCEEEIIoSOzZ89Otg/lp1y6dIm9e/dq9RytWzjMzc25f/8+jRs3TrLv/v375MuXT+O6WrVqxenTp9O0oJv4/EgLhhBCaObRo0f07t2bFy9eYGhoyNmzZ5Ptewmwbt065s+fT58+ffjmm28oV65csuViYmLYu3cva9eu5ciRI6xbt06rmLROOFxcXPDy8sLDw0NtCGxoaCg//fSTVsnDDz/8QKdOnciXLx9t27alUKFCScpoM8xWCCGEENCzZ0+mT59Ow4YNCQ0N/WQfyWvXrrF06VJ+/vlnvLy8KFKkCDVq1KBIkSKYmpoSGhrKgwcPuHbtGnFxcbRu3Zp//vmHihUrahWT1gnHlClTqFWrFqVLl8bT0xNbW1v8/f3ZunUrsbGxKU5hnZzEhdK+//57vv/++2TLaDpKRQghhBBw48YNjI2NadiwIZD6D3cjIyOGDh3KoEGD2LlzJ/v27ePvv//mzJkzREdHY2lpSbly5Zg4cSJdu3ZN82gVrROOsmXLcvr0aUaMGMGqVauIj4/H0NAQZ2dn5s2bR9myZTWu61NDcYQQQuifjNrIek6dOsWcOXO4dOkSgYGB7Nixgw4dOqiV8fb2Zs6cOQQGBlKxYkUWLFigSjDu3btH3rx5adeuHf7+/nTq1Inx48enelxDQ0M8PDzw8PDQx8tK28RfVatW5ejRo0RHRxMWFkbBggVVq6ZqI6OG4gghhBDZxevXr6latSq9evWiY8eOSfb7+voybNgwvL29qV+/PitWrMDNzY2bN2/i4OBAbGwsp0+f5vLlyxQpUoSWLVtSq1Ytmjdvngmv5v9pPUrlQ2ZmZtjY2KQp2fhYdHQ0AQEBxMXFpbsuIYQQIiuJjo7mzZs3qn+xsbEplnVzc2P69OkptjTMmzePPn360LdvX8qXL8+CBQuwt7dn2bJlANjZ2VGrVi3s7e3JlSsXrVq14vLly/p4WVpJUwtHfHw8+/fv59atW0mG0igUilQXcPvQ8ePHGT9+PBcuXADg/Pnz1KhRg4EDB9K0aVO9Ne2IrEVGoAghcrIpU6bw6NEj1eNOnTrRuXNnret59+4dly5dYuzYsWrbXV1dOXPmDAC1atXixYsXhIWFUaBAAU6dOsV3332XvhegA1onHCEhITRs2JDbt2+jUChUK+N92BdD04Tj2LFjtGjRgkqVKjFq1Chmz56t2mdpaYmPj48kHEIIIbK9woUL8/btW5o0aUKzZs00XgLkY8HBwcTHx2NlZaW23crKSjUPlpGRETNnzqRRo0YolUpcXV1p06ZNul9DemmdcEyYMAFTU1OePHlCsWLFOHfuHAULFmT58uXs2bOHI0eOaFzXpEmTaNWqFTt37iQuLk4t4ahatarWY3yFEEKIrGjUqFE6re/jARdKpVJtm5ubG25ubjo9Znpp3Yfj6NGjjBgxAhsbm/cVGBhQsmRJ5syZQ7NmzbQ6qf/++6+qmefjk1e4cGGCgoK0DU8IIYTIcubOncvw4cM5cOBAuuqxtLTE0NAwyazeQUFBSVo9shqtWzj8/f0pXrw4hoaGGBgY8Pr1a9W+tm3b0rVrV80PbmSUYseZoKAgrWYtFUIIIbIqXa0Wa2JigpOTE4cPH8bd3V21/fDhwzpdAf327ducPHmS4OBg+vTpg7W1Nc+ePcPCwgIzM7M01al1wmFpacmrV68AsLGx4fr16zRq1Ah4P9uoNqNMatWqxa+//prsSdq2bRt169bVNjwhhBAiy5k7dy7+/v60aNGCli1bfrJsVFQU9+/fVz1+9OgRly9fpmDBgjg4ODBixAi6detGzZo1qVu3LitXrsTPz4/+/funO874+Hj69euHj4+P6jaNm5sb1tbWfPfdd1SvXp1p06alqW6tEw4nJydu3LhB69atadWqFdOmTSN//vyYmJgwfvx46tSpo3FdY8eOpUWLFri7u9O9e3cUCgXnzp1j7dq1bNu2jePHj2sbnhBCCJHlaNPCcfHiRVxcXFSPR4wYAUCPHj3w8fHB09OTkJAQpk2bRmBgIJUqVWLfvn0UK1Ys3XHOmDGDjRs3MmfOHFq2bKmaERze9wvx8fHJuIRj0KBBPHjwAIAff/yRs2fP0r17dwBKlizJwoULNa6rWbNmrF+/nmHDhrFz504ABg4ciLm5OT4+PjRo0EDb8IQQQohsrXHjxqoRoCkZMGAAAwYM0PmxfXx8mDhxIiNGjEiytEiJEiXUhvZqS6OEIyIignz58qFQKGjWrBnNmjUD3nfs/Pfff7l+/ToKhYJy5cphZKRdDvPNN9/QsWNHzpw5w4sXL7C0tKR+/foprmonhBBCZDfa3FLJTAEBASl2ZzA1NSUyMjLNdWuUHVhYWPD3339Tu3ZtevfuzcSJEylRogTwfnRJ5cqV0xwAvJ+xtGnTpumqQwghdEHWFhH6oKtOo/pWpEgRHj58qHZLJ9GdO3ews7NLc90aDYs1MjJSNa34+Pjw8uXLNB/wQ+vWrUtxPZUpU6bwyy+/6OQ4QgghhEhdq1atmDFjBgEBAaptCoWCV69esWjRItq2bZvmujVq4XBwcGD9+vWqmdHu3LnzyVsnNWrU0OjgixYtomfPnsnus7S0ZNGiRar+IUIIIUR2lV1uqUybNo39+/dToUIFXFxcUCgUjB8/nuvXr2NsbKzV0iUf0yjhGDJkCEOHDmXVqlUoFIoUk4TEITQfdzRJyf3799V6wH6oQoUK3Lt3T6N6hBBCiKwsu9xSsbKy4sKFC0yePJm9e/diaGjIlStXaNOmDdOmTaNgwYJprlujhGPw4ME0atSI69ev061bN3744QdKliyZ5oN+KHFOj+S2y8qxQgghRMaysrJi+fLlOq9X4yElVatWpWrVqqxevZquXbtSrly5dB+8cuXKbN68OdkF2jZt2pTuzqhCCCFEVpBdbqnok1ZjWKOjo7lz5w737t3TScIxaNAgvvnmG3r06MGAAQOws7PD39+fZcuWsX37duk0KoQQIkfILrdUevfuneI+AwMDzM3NqVWrFu7u7piYmGhVt1YJh5mZGTExMTqbI6Nr167cvn0bLy8vfvvtN9V2AwMDfvjhB77++mudHEcIIYQQqTt+/DivXr0iPDwcIyMjChUqREhICHFxcZibm6NUKpk3bx5ly5blxIkTWi0Yp/VqsU2bNtVqCfpPeffuHVOmTOHevXssX76cH3/8keXLl3P//n2mTp2qk2MIIYQQQjPbt28nX758bNq0iejoaAIDA4mOjmbjxo3ky5ePgwcP8ueffxIWFsb48eO1qlvrqc3Hjx9Px44dMTU1xcPDg6JFiyZZWl6TXqyJLSXbtm3D3d2db7/9VttQhBBCCKFDI0aMYNSoUXh6eqq2GRoa0qVLF168eMGIESP4888/GTNmDD///LNWdWvdwuHk5MTjx4+ZMmUKVatWpUiRIhQuXFjtnyZMTU0pVKiQTGEuhBAix5s7dy7Dhw/nwIEDmR3KJ124cIEKFSoku69SpUr8+++/AFSrVo3g4GCt6ta6hWPSpElJWjTSqm3btuzYsQNXV1ed1CeEEEJkRdml02j+/Pk5fvx4ssuNHDt2jPz58wPvB5Hky5dPq7q1TjhSmoo8Lbp06UKfPn3o3bt3irdnNJ21VAghhBDp07VrV2bNmoVSqeTLL7/EysqKFy9e4Ovry9y5cxk6dCgAly5donz58lrVrXXCoUstWrQA3q/Psn79erV92s5aKoQQQoj08fLyIjAwEC8vL3766SfVdqVSyVdffcXMmTMBqFu3ruo7XFNaJxzTpk375H6FQqHxXOvr1q3T9vBCCCGE0BMTExM2btzIxIkTOXnyJCEhIRQqVIhGjRqp9e1o1qyZ1nXr/JaKNglHjx49tD28EEIIIfSsfPnyWt8ySY3WCUdCQkKSbaGhofzxxx8sWLCAvXv3pimQO3fuEBwcTLVq1dI0cuXmzZvs2rWLR48eERYWxqhRo6hdu3aqz1m/fj3+/v5YWFjQrl076cAqhBBC57Lj1OYvX74kOjo6yXYHB4c01aeTPhwFCxakd+/eBAUFMWTIEHbs2KHxc3/55RfGjx9PYGAg8H5ITo0aNejcuTPNmzfXeH6Ot2/fUrx4cVxcXJg7d26q5YOCgvDy8qJp06YMHjyYO3fusHr1avLnz0+dOnU0jl8IIYRITXYZpQIwffp0Fi1aREhISLL709q3Uut5OD6ldu3aHD16VOPyW7dupWfPntSoUYMlS5agVCpV+2rUqMGWLVs0rqt69ep06dKFL774QqPyhw4dwtLSkp49e2JnZ0fTpk1xcXFh9+7dGh9TCCGEyEnWrl3LTz/9xJAhQ1AqlYwfP55x48ZhZ2dH6dKlWb16dZrr1mnCceXKFfLmzatxeS8vL3r16sWuXbvo16+f2r7y5ctz8+ZNXYan5t69e1SpUkVtW7Vq1Xj48CFxcXF6O64QQgiRVS1dulSVZAC4u7szffp0bt++Tb58+bSe7OtDWt9SSW4F17dv33L16lXWrl3LN998o3Fdt27dYtasWcnuK1iwYIrNOboQHh5OgQIF1LYVKFCA+Ph4IiMjsbCwSPKc2NhYYmNjVY8NDAwwNTUF0HgytNTKaVKPLsrkpDoy6jhZpY6MOs7nFqu83qx7nKxSh66Ok5Xdv3+fOnXqYGDwvj3i3bt3wPvFW0eOHMnEiRP5/vvv01S31glHz549k91uamrKN998o9Xc6rlz5+bVq1fJ7gsICEj2S1+XPr4wEm/ppHTB7Nixg23btqke169fn6FDh2oVZ9GiRdO1X1dlclIdGXWcDK8jIlD/x0nHMdTKpFaPJsfJiDo+Qas6stPrTaGMLs6ZJmWyyvtGp39fTeLIZEZGRlSqVAmAmjVranwrxMjofVqgUCjInz8//v7+qn2WlpYEBASkPSZtn/Do0aMk20xNTbVaojZR/fr1WbJkCR07dkyyz8fHh8aNG2tdp6bMzc0JDw9X2xYREYGhoWGKt4Xc3d1p06aN6nFiBhgWFkZMTIxGx03sHJvW/boqk5PqyKjjZJU6Muo4n1us8nqz7nGySh1pOY6pqanefzwnx9zcnMuXL2v9vNKlS/P06VMAatWqxapVq2jfvj0GBgasXLmS4sWLpzkmrROOYsWKpflgH5s0aRINGjSgdu3adO3aFYVCwe+//87kyZM5deoU58+f19mxPla6dGkuXbqktu3KlSs4OjqqMryPGRsbY2xsnOy+Dzu8fkpq5TSpRxdlclIdGXWcjK8j5aZZ3R0n7cdQL5NaPZocJyPqSJl2dWSn15t8GV2cM03KZJX3jW7/vprEkT21atWKU6dO0aNHD8aNG0eLFi0wNzfHyMiIqKgo1q5dm+a609Vp9N27dyxfvpxBgwYxc+ZMXrx4odXza9asyf79+4mKimLkyJEolUpmzpzJ3bt32bdvn6o5SBMxMTE8fvyYx48fA++HvT5+/FjVwWXjxo0sWbJEVd7V1ZXg4GDVPBzHjh3j2LFjtG3bVqvXIIQQQujSqVOnaNu2LTY2NigUCv74448kZby9vSlRogSmpqY4OTlx+vRptf0RERE4OTnRoEEDTp48qfGxJ02axJo1awBo0qQJZ86coX///nz77bccOXIkXRN2atTCMWvWLLZv367W4hAXF0eDBg24dOmSKqNbunQpFy5cwMbGRuMAXFxcuHXrFg8ePODFixdYWlpSpkwZLV8GPHjwgKlTp6oeJ3ZudXZ2ZuDAgYSFhan1ri1SpAjjxo1j/fr1HDx4EAsLC3r16iVzcAghhMhUr1+/pmrVqvTq1SvZLge+vr4MGzYMb29v6tevz4oVK3Bzc+PmzZuqSbkeP36MjY0N169fp3Xr1ly7dk210qs2atWqRa1atdL9mkDDhOPAgQNJDrhy5UouXryIp6cn48aN4+bNm/Tv35+ZM2eqtSQkZ//+/SxYsIAnT55QuHBhvvnmG7777jtKliyZ5hdSsWLFT87bMXDgwCTbKlSokOIoGSGEEEJXoqOj1W63fOoWvZubG25ubinWNW/ePPr06UPfvn0BWLBgAQcPHmTZsmV4eXkBqH74V6pUiQoVKnD37l1q1qypq5eTJholHHfv3lW9sEQ7duygQIECrF27FjMzM6pUqcLdu3f57bffPlnXvn37aNu2LUqlksKFC3P//n3OnDlDSEgI48ePT/srEUIIIbKoKVOmqA266NSpE507d9a6nnfv3nHp0iXGjh2rtt3V1ZUzZ84A7wcy5M6dm1y5cuHv78/NmzdxdHTUqP74+HgWL17Mhg0bePLkSZIBEQqFIsXRpanRKOEICQnB3t5e9TghIYEzZ87QvHlzzMzMVNvr1Kmjyq5SMmfOHMqUKcPevXtxdHQkLCyMr7/+mp9//plx48Zl+zHMQgghxMemTJmSpIUjLYKDg4mPj08yMtTKyornz58D7+e4+u677zAwMEChULBw4UIKFiyoUf1jx45l7ty5VKtWjebNm2NiYpKmOJOjUcJRqFAhQkNDVY+vX79OdHR0kuYZU1PTFEd4JPr333/x9vZWZVsWFhbMmzePihUr8uDBA0qVKqXtaxBCCCGyNG9vb50u3pbcPFKJ2+rVq8e1a9fSVO+GDRsYO3YsM2fOTHeMH9Mo4ahUqRKbNm2iQ4cOwPvbKQqFIsk8GQ8ePMDa2vqTdUVERCRp2ilVqhRKpZLIyEjNIxdCCCGyCV0t3mZpaYmhoaGqNSNRUFBQmubD+lh0dDTNmjVLdz3J0SjhGDFiBG5ubrx8+RJra2u2bdtGpUqVaNCggVq5/fv3U61atVTrS5wwK1FiVpbdxy8LIYQQydHV8vQmJiY4OTlx+PBh3N3dVdsPHz5M+/bt0x2nq6sr586do0mTJumu62MaJRwtWrTA29ubn376iQsXLtC4cWOWL1+uVub58+ccPXpUo1Efc+fOVcvEEpuC5syZQ+HChVXbE+89CSGEENmZNi0cUVFR3L9/X/X40aNHXL58mYIFC+Lg4MCIESPo1q0bNWvWpG7duqxcuRI/Pz/69++f7jgXLVpE69atyZMnD61atUq274em/UE+pvFMo/379//ki7G2ttZ4sbWtW7cmu93X11ftsSQcQgghcgJtWjguXryIi4uL6vGIESMA6NGjBz4+Pnh6ehISEsK0adMIDAykUqVK7Nu3TyczgefPn5+yZcsyfPhwhg8fnmyZ+Pj4NNWt9dTm6ZWQkJDRhxRCCCEylTYtHI0bN061i8GAAQMYMGBAesNKon///mzZsoX27dtTvnz5jB+lIoQQQoi001UfDn37448/8PLyYtSoUTqvO8MTjkqVKvHjjz+qdXb5lMDAQLy8vLCxsUky0YkQQgiRHehqlIq+GRsbU716db3Una7F29Kic+fOdO/eHQcHB8aNG8fBgwd5+fKlqvkoOjqa69evs3r1atq2bUuxYsW4dOkS7dq1y+hQhRBCiM+Kh4cHBw8e1EvdGd7CMWnSJL799lsWLFjA6tWrmTVrFgqFAoVCgbGxMe/evQPej1xp2LAhmzdvxsPDI6PDFEIIIXQmu9xS+eqrr/j222+JjY2ldevWyY5IqVGjRprqzpQ+HEWLFmXWrFlMnz6dc+fO8ffff/Ps2TOio6OxtLSkXLlyNG7cGDs7u8wITwghhNCp7HJLpWnTpgAsXLiQRYsWqe1LnMIiw0apXL16lfDwcBo1agS8Hy88evRo/vnnH1xdXZk6darG66EYGxvToEGDJBOICSGEECLjrVu3Tm91a51wjBgxgho1aqgSjgkTJrBq1SoqV66Ml5cXhQsXZvDgwToPVAghhBD61aNHD73VrXWn0evXr1OvXj3gffPKhg0bmDp1Kv/88w9jxoxh7dq1Og9SCCGEyM7mzp3L8OHDOXDgQGaHkmm0buEIDw/H0tISgCtXrhAWFkbnzp2B9/d+Fi9erNsIhRBCiGwuu/ThAAgNDWXjxo3cunWL6OhotX0KhYI1a9akqV6tE45ChQrx9OlTAI4fP46VlZVqSfl3797JAmxCCCFENuXn50etWrV48+YNb968wdLSktDQUOLj47GwsKBAgQJprlvrhKNhw4ZMmTKF4OBg5s+fT+vWrVX77t27h729fZqDEUJkTaV4SW8uUJwwHmPBWmpxn8KpP1EIHZNrUb/Gjh1LxYoV2bNnD3nz5mX//v1UqlSJVatWMXPmTPbu3ZvmurVOOLy8vHBzc2Po0KGULFmSSZMmqfZt3bqVOnXqpDkYIUTW05MLrGYrShQoUKJEwWhO0Icvga6ZHZ74jMi1qH9///03s2fPxtTUFHjfV9PExISBAwfy4sULvv/+e/bs2ZOmurVOOEqUKMHt27cJDQ1NMiHIkiVLsLa2TlMgQoisx/DhQ1azFUOUQOLtUiVKYA1bCX40hPgSJTIxQvG5yO7XYnaZ+OvFixcULVoUAwMDDA0NiYiIUO1zdnZOMjeHNtI8tXlys49VrlyZwoWlaUuInCK3ry9Kks6rowCUKMi9eXPGByU+S9n9Whw5ciTz58/P0skGgJWVFaGhoQAUL16cixcvqvY9fvwYI6O0zxeq9TN/+eWXVMt07949TcEIIbIWw6dPUZBSR3Alhv91IBdC3+RazBh16tTh33//pV27dnh4eDBt2jTevn2LiYkJc+bMoUmTJmmuW+uEo2fPnslu/3B2UUk4hMgZ4u3t//tVmdwHvYJ46SQuMohcixlj1KhRPH78GHi/9tmtW7eYPHkySqWSRo0asXDhwjTXrXXC8ejRoyTbgoOD2blzJ76+vmzO4s1aQgjNvfH0JPeSpShBrTH7/WMlb7p0yaTIxOdGrsWM4eTkhJOTEwB58uRh165dREREoFAoyJcvX7rq1jrhKFasWLLbnJyciI2NZeHChfj4+KQrKCFE1hDv6EgfvmQNW0n4b2QA//23D18yMwt30hM5i1yLmSd//vw6qUenq8U2bdpUNeuoECJnWE8t/qQEfTivmvtgDbV5gCUzMzs48VnJztdidhil8vLlS1asWMGpU6d49uwZADY2Nri4uNCvXz8KFSqUrvp1mnA8efIEQ0NDXVYphMgCHmDJeFpldhgim9PFpF3Z9VrM6lObHz16lI4dOxIREYGhoSGWlpYolUru3LnDkSNH+Pnnn9mxY4dq4da00DrhOHXqVJJtb9++5erVq3h5edG0adM0ByOEECJnkkm7sq6XL1/i6elJgQIFWL16Na1atSJ37twAvHnzhj179jBq1Cg6derErVu30tzSoXXC0bhxY7URKYBq/ZRmzZrJ4m1CCCHUaDJpl8g8a9asIT4+nr/++gs7Ozu1fblz56Zz587UqVOHqlWrsmbNGkaPHp2m42idcBw/fjzJNlNTU4oXL46VlVWaghBCCJFz/f+kXepDWtUn7Uo6IEFkjEOHDtG7d+8kycaHHBwc6NWrFwcOHMi4hMPZ2TlNBxJCCPF50mzSLkk4dO3NmzeUL1+eL7/8kp9//jnFcrdu3WLw4MGp1tewYUM2bdqU5njS3Gn0/v37HDt2jJCQECwtLXFxcVEtUy+EEEIkkkm7MseMGTP44osvUi0XHh5OkSJFUi1XpEgRwsPD0xyP1gmHUqlk8ODBLF++nISEBNV2AwMDBgwYkK6FXYQQQuQ8Gk3ateRvIPmRLEJ79+7d4/bt27Rt25br169/suzbt28xNjZOtU4jIyPevXuX5pi0Xrxt/vz5eHt7891333Hu3DmePn3KuXPn6N+/P97e3syfPz/NwQghhMh5EiftSkBBLAbEoSAOAxJQ0IcvVau89uQCt5nD95ykM1f4npPcZg5mvr6Z/Aoy1qlTp2jbti02NjYoFAr++OOPJGW8vb0pUaIEpqamODk5cfr0abX9o0aNwsvLS+Nj3rlzh3/++eeT/27fvp2u16V1C8fq1asZPHiw2nzqtra21KpVC0NDQ1atWsXw4cPTFVRaHTx4kF27dhEeHo6dnR09e/akfPnyyZa9ceMGU6dOTbJ9/vz52Nra6jtUIYT4rKQ2aVcpXqY4ksV81Cje1a6decFnsNevX1O1alV69epFx44dk+z39fVl2LBheHt7U79+fVasWIGbmxs3b97EwcGBnTt3UqZMGcqUKcOZM2c0OmZK66R9SKlUJhmlqg2tE46HDx/Spk2bZPe1adOGFStWpDmY9Dhz5gw+Pj707duXsmXLcuTIEWbOnMn8+fOxtLRM8XkLFixQjTcG3U3hKoQQQt2nJu3qzYWUR7Iosv9IlujoaNUUEgDGxsYp3sZwc3PDzc0txbrmzZtHnz596Nu3L/D+e+zgwYMsW7YMLy8vzp49y+bNm9m6dStRUVHExsaSP39+Jk2alGx969atS8cr05zWCUeBAgV48uRJsvuePHmSaV/Ye/bsoUmTJqqJx3r27MmVK1c4dOgQXbumPKlMgQIFyJMnT0aFKYQQIhnFCUt5JIsy+49kmTJlitrip506dUrTUiDv3r3j0qVLjB07Vm27q6urqjXDy8tLdTvFx8eH69evp5hsAPTo0UPrONJC64SjefPm/PDDD1SvXl21ohzA5cuXmTx5Mi1atNBpgJqIi4vj4cOHdOjQQW17lSpVuHPnziefO3r0aGJjY7Gzs8PDw4NKlSqlWDY2NpbY2FjVYwMDA0xNTQE0bmZKrZwm9eiiTE6qI6OOk1XqyKjjfG6xyuvN3OM8xiLlkSwKBQmpjGTJ6q93ypQpSVo40iI4OJj4+Pgk815ZWVnx/PnzNNWZUbROOLy8vDhx4gS1a9emQoUKFC1alMDAQG7evImNjY1WnVR0JSIigoSEBAoUKKC2vUCBAikO4bGwsKBfv344OjoSFxfHqVOn+PHHH5k8eTIVKlRI9jk7duxg27Ztqsf169dn6NChWFhYaBxr0aJF07VfV2VyUh0ZdZwMryMiUP/HSccx1MqkVo8mx8mIOj5Bqzqy0+tNoYwuzpkmZTStYy21GM2J5EeyKJXkHToUlvyR+nF08Xr18N7z9vbW6eJtyc34nVyyo0nfjIyidcJhb2/P5cuXmTt3LsePH+fRo0cUKlSIsWPHMnz48E/2l9C35E52StmmjY0NNjY2qsdlypQhODiY3bt3p5hwuLu7q/VfMTB4P8gnLCyMmJgYjWIMDEz5QtZkv67K5KQ6Muo4WaWOjDrO5xarvN7MPc59Cqe4/Hzk3LlEp3LrO6u+XlNTUywsLHS2eJulpSWGhoZJWjOCgoKy/GzfaZr4y9LSMlNaMlKSP39+DAwMkrRmvHr1Kkmrx6eUKVMmydCiD32qk8+HTWWfklo5TerRRZmcVEdGHSfj60i56VZ3x0n7MdTLpFaPJsfJiDpSpl0d2en1Jl9GF+dMkzLa1JHSSJZTnTtDOv9+uv37ahKHOl0tT29iYoKTkxOHDx/G3d1dtf3w4cO0b98+zfVmBJ0uT59ZjIyMcHR05OrVq9T+YOjU1atXqVVL80ljHj16hLm5uR4iFEIIoYnsuvx8arRp4YiKiuL+/fuqx48ePeLy5csULFgQBwcHRowYQbdu3ahZsyZ169Zl5cqV+Pn50b9/f32ErjMaJRy9e/fWuEKFQsGaNWvSHFBatWnThsWLF+Po6EiZMmU4cuQIwcHBNG/eHICNGzcSGhrKoEGDANi7dy+FCxfG3t6euLg4Tp8+zblz53TW7CWEEEIk0qaF4+LFi7i4uKgejxgxAng/msTHxwdPT09CQkKYNm0agYGBVKpUiX379lGsWNYexaNRwnHs2DGdjcLQl3r16hEZGcn27dsJCwvD3t6ecePGUbhwYeB9P4vg4GBV+bi4OH799VdCQ0MxMTHB3t6esWPHUqNGjUyJXwghRM6lzY/Zxo0bp3rrZsCAAQwYMCC9YWUojRKOx48f6zkM3WjRokWKw3IHDhyo9rh9+/ZZ/n6XEEKInEFXfTiysxzRh0MIIYTIyuR2fRoWbxNCCCGE0FaaEo7ffvuNmjVrkidPHgwNDZP8E0IIIcT/mzt3LsOHD+fAgQOZHUqm0fqWyq5du+jVqxc9e/bkn3/+oXfv3sTExLBr1y5sbGz46quv9BGnEEIIkW3JLZU0tHD89NNPjBgxguXLlwPve8r+9ttv3L17l/j4eOxTme9eCCGEEJ8frROOO3fu0KxZM9Xw17i4OACsra354YcfmDdvnm4jFEIIIUS2p/Utlfj4eExMTDAwMCBPnjxq87k7ODjw8OFDnQYohBBCZHcyLDYNCUeJEiV49uwZAFWrVmXTpk20a9cOgG3btmm0MqAQQgjxOZE+HGm4pdK0aVOOHDkCwNChQ/H19aVUqVJUqFCB5cuXZ/m53IUQQgiR8bRu4ZgxYwZv374F4Msvv8TQ0JANGzagUCgYPXo0PXv21HWMQgghhMjmtE44cuXKRa5cuVSPPTw88PDw0GlQQgghRE4ifThkanMhhBBC76QPRxoTjj/++IMNGzbw5MkTYmJi1PYpFAquXLmik+CEEEIIkTNonXDMmTOHMWPGULhwYUqVKkWePHn0EZcQQgghchCtEw5vb2969+7NihUrZN0UIYQQQmhE62GxISEhdO3aVZINIYQQQmhM64Sjfv363Lp1Sx+xCCGEECKH0vqWyoIFC3B3d8fe3p6WLVtiYmKij7iEEEKIHEOGxaYh4ShVqhTNmjXD3d0dhUJB7ty51fYrFApevXqlswCFEEKI7E6GxaYh4Rg9ejRLliyhWrVqlC9fXlo4hBBCCJEqrRMOHx8fxowZg5eXlz7iEUIIIUQOpHWn0fj4eJo3b66PWIQQQgiRQ2mdcLi6unL27Fl9xCKEEEKIdIqMjKRWrVpUq1aNypUrs2rVqswOCUjDLZWJEyfi6elJnjx5aN26NQULFkxSJrltQgghhNC/3Llzc/LkSXLnzs2bN2+oVKkSHh4eFCpUKFPj0jrhqFq1KgAjRoxgxIgRyZaJj49PX1RCCCGylFK8pDcXKE4Yj7FgLbW4T+FUy4iMZ2hoqBpBGhMTQ3x8PEqlMpOjSkPCMWnSJBQKhT5iEUIIkQX15AKr2YoSBQqUKFEwmhP04Uug6yfLRPhaE+3pmbkvIJs5deoUc+bM4dKlSwQGBrJjxw46dOigVsbb25s5c+YQGBhIxYoVWbBgAQ0bNlTtDw8Px9nZmXv37jFnzhwsLS0z+FUkpXXCMWXKFD2EIYQQIisyfPiQ1WzFECWQ+CtZiRJYw1aCHw2hFC9TLGM+ahTvatfOjNCzrdevX1O1alV69epFx44dk+z39fVl2LBheHt7U79+fVasWIGbmxs3b97EwcEBAHNzc65cucKLFy/w8PCgU6dOWFlZZfRLUaN1p9EPxcTEEBgYmGSJeiGEEDlDbl9flCRt1VYAShTk3ryZ3lxIsQyK92U+d9HR0bx580b1LzY2NsWybm5uTJ8+HQ8Pj2T3z5s3jz59+tC3b1/Kly/PggULsLe3Z9myZUnKWllZUaVKFU6dOqWz15JWaUo4zpw5Q8OGDcmXLx92dnbky5cPZ2dn/v77b13HJ4QQIhMZPn2KgpTu/ysxfPqU4oSlXEb5vsznbsqUKfTs2VP1b8eOHWmq5927d1y6dAlXV1e17a6urpw5cwaAFy9eEBERAUBERASnTp2ibNmy6XsBOqD1LZWzZ8/SpEkTzM3N6devHzY2NgQEBPD777/TpEkTTpw4wRdffKGPWIUQQmSweHv7/1ovkksoFMTb2/OYiJTLKN6X+dwVLlyYt2/f0qRJE5o1a4axsXGa6gkODiY+Pj7J7RErKyueP38OgL+/P3369EGpVKJUKhk0aBBVqlRJ92tIrzR1Gq1SpQrHjx8nT548qu1z5szBxcWFSZMmcfDgQZ0GKYQQInO88fQk95KlKEHtpsn7x0redOnC2iW7GM2JZMugfF+GJZ93C/ioUaN0Wt/HgzeUSqVqm5OTE5cvX9bp8XRB61sqZ8+eZfTo0WrJBkCePHn4/vvv5baKEELkIPGOjvThSxJQEIsBcSiIw4AEFPThS+JLlOA+hVMsE/7zz8SXKJHZLyPTzZ07l+HDh3PgwIF01WNpaYmhoaGqNSNRUFBQpncKTY3WLRzx8fHkypUr2X2mpqYyB4cQQuQw66nFn5SgD+dVc2ysoTYPsGRmKmVOyZBYQHerxZqYmODk5MThw4dxd3dXbT98+DDt27fXyTH0JU0Tfy1btoy2bdsm2bdixQrVxGBCCCFyjgdYMp5W6S7zuZo7dy7+/v60aNGCli1bfrJsVFQU9+/fVz1+9OgRly9fpmDBgjg4ODBixAi6detGzZo1qVu3LitXrsTPz4/+/fvr+2Wki9YJx9ixY+nQoQPVq1fnm2++oWjRogQGBrJx40YuX77MH3/8oYcwNXPw4EF27dpFeHg4dnZ29OzZk/Lly6dY/ubNm6xfvx5/f38sLCxo165dkp6/QgghRHpp08Jx8eJFXFxcVI8TZ/Xu0aMHPj4+eHp6EhISwrRp0wgMDKRSpUrs27ePYsWK6TxuXdI64WjXrh2//fYbo0eP5vvvv1dtt7W15bfffku25SMjnDlzBh8fH/r27UvZsmU5cuQIM2fOZP78+cnOsBYUFISXlxdNmzZl8ODB3Llzh9WrV5M/f37q1KmTCa9ACCGEgMaNG6c6FfmAAQMYMGBABkWkG1onHABdu3blq6++4s6dO4SEhFCoUCHKli2bqVOe79mzhyZNmtC0aVMAevbsyZUrVzh06BBdu3ZNUv7QoUNYWlrSs2dPAOzs7Hjw4AG7d++WhEMIIYROaXNLJadKU8IB74fklCtXTpexpFlcXBwPHz5MMtd8lSpVuHPnTrLPuXfvXpJxydWqVeP48ePExcVhZJT01MTGxqrNDmdgYICpqWmS8dSVKlVKMVYzM7NPlkltvzZlskodnyojrzflMnYplgAzM/jqq5Tff4l1pFYmtWNoUgekHqsmx8mIOiDl16Ppa9GkTFap41NldHHO3pcx08m1qIs6QDevVxfvvUSJ3w+66jSanSmUaVhCLjIykv379/PkyROio6PVK1QomDhxos4C1ERoaCj9+/fnxx9/VJtN7ffff+fkyZMsXLgwyXOGDh2Ks7Oz2tSxd+7cYeLEiaxYsQILC4skz9myZQvbtm1TPe7Vqxdubm46fjVCCCFEzqP1PBznzp2jRIkSdOnShTFjxjBlypQk/zJLcrd0PnWbJ7mJUz71HHd3d3x8fFT/ateuzf79+7PFWjIxMTEsXLhQYtWh7BInSKz6IrHqh8SaM2mdcAwfPhxbW1vOnz9PTEwMCQkJav8yYx6O/PnzY2BgQHh4uNr2V69eUaBAgWSfY25unqR8REQEhoaG5M2bN9nnGBsbkzt3btU/MzMz1q1bR0JCgi5ehl4lJCTw119/Saw6lF3iBIlVXyRW/ZBYcyatE45r164xffp0atasiYmJiT5i0pqRkRGOjo5cvXpVbfvVq1dTXLCmdOnSScpfuXIFR0fHZPtvCCGEECLttE44ChcurI840q1NmzYcPXqUY8eO4e/vj4+PD8HBwTRv3hyAjRs3smTJElV5V1dXgoODVfNwHDt2jGPHjmXasF4hhBAiJ9P6p/zgwYNZvnw5bdq0ydRhsB+rV68ekZGRbN++nbCwMOzt7Rk3bpwqQQoLCyM4OFhVvkiRIowbN47169dz8OBBLCws6NWrl1ZDYo2NjenUqVOaV/3LSBKr7mWXOEFi1ReJVT8k1pxJ61Eqc+fOxdvbm3z58tG6dWsKFSqkXqFCwfDhw3UapBBCCCGyN60TDgODT9+FUSgUsoCbEEIIIdRofUvl0aNH+ohDCCGEEDlYmib++pSXL19m2Y6lQgghhMgcWo9SSY5SqWTfvn107NgRO7tPTQorhBBCiM9RuiacePDgAWvXrmX9+vUEBgZiYmJCx44ddRWbEEIIIXIIrROOmJgYtm7dypo1azh9+jRKpRKFQsGIESMYO3ZsklErQgghhBAaJxwXLlxgzZo1bN68mcjISPLkyUPPnj3p2LEjbdq0oW3btpJsZKKEhIRURxAJ7cl51Y3EdSbevn1LgQIFsvR5jYiIICwsDIBixYplcjSfJudVP7LTec1ONEo4qlSpwo0bNwCoW7cuvXv3xtPTkzx58vDq1Su9BpiRgoKCuH79OqGhoTRo0AArK6ssNbnZh8LDw3nx4gWRkZHUrFkTAwODLPumkPOqH9nlvPr7+7N+/XqioqIIDw+nf//+VK1aNbPDSpafnx8LFixAoVDg7+9Px44d6dy5c2aHlSw5r/qRnc5rdqPRKBUDAwMUCgWtW7fmp59+okKFCqp9r169wsLCghMnTtCoUSO9BqtPfn5+zJgxA2tra549e4aRkREzZsygYMGCmR1aEn5+fixcuBCFQkFISAh2dnb8+OOPmR1WsuS86kd2Oa9+fn5MnjwZZ2dnihUrxt27d7l8+TI///wzefLkUd2SzQqeP3/OxIkTadq0KQ0bNuTevXssW7YMb2/vLNd6K+dVP7LTec2ONGrhWLBgAevWrWPPnj3s3buX2rVr06dPHzw9PfUdX4YIDQ3l559/xtXVlTZt2pArVy4GDx7MzZs3adCgQWaHp8bf35+pU6fSrFkzGjVqRGRkJHPnzuX27duUK1cus8NTI+dVP7LLeQ0ODmbhwoU0b96crl27AmBpaUlERARKpZLQ0NAslSCdOHGCkiVL0qVLFwAKFSrE33//zatXrwgODsba2jrF1aczkpxX/chu5zU70ijhGDJkCEOGDOHixYuqfhz9+vVj2LBhtG7dGoVCka2zvqdPn2JmZoaLiwu5cuUCwM7OjqCgINauXUv16tVxdHTM9DdFVFQUy5Ytw9nZma+++gqA+Ph47O3tCQ0N5cSJE1SrVg1zc/NMjTORnFf9yC7nNTw8HDs7O5o1a6baduPGDW7evMnkyZMJCwujRYsWtG/fHlNT00yM9L3Q0FAMDAyIi4vDyMiIffv2ceXKFSIiInj27BmVKlXC3d2dUqVKZWqccl71I7ud1+xIq1EqNWvWpGbNmsyfP181UmXbtm0olUr69OnDd999R8+ePbNcM1lqwsPDef78OVFRURQsWJDdu3dz+fJlcufOTVBQEJcvX6ZZs2a0adMmU+/l582blzp16qi9Mf/44w9u375NdHQ0b968YcOGDYwaNYqyZctmevNfdjqvdevWpWTJkqptcl7Tr1SpUnTv3l31eXD06FF2797Nt99+i4ODA8+ePWPx4sWUKFGC2rVrZ1qcicqVK8eKFStYtmwZSqWSs2fPMnLkSCpWrEhQUBDz58/nn3/+yfQvxlKlStGjRw/Vr205r7qR3c5rdpTumUYfPHjAmjVr+OWXX3j27Bmmpqa8efNGV/HpTWRkJAkJCapfgePGjSMoKAhHR0du3LjB2LFjqVKlCgDr16/nwoUL/PTTT+TNmzczw1bzzz//sGbNGnr16kW5cuXImzcvU6dOBWDy5MmZHN17Y8eO5eXLl1n2vMbHx2NoaKi2Laue1w8Tnex2vcbHx3PixAns7OwoW7asavuYMWMoX748PXv2zLzgPnDkyBHCwsJ48uQJ5ubm9O3bV9Vp2Nvbm5cvX/LDDz8kuWYyQuJH9YfJblY9rx/HeuzYMUJCQrLkef1QYkxZ9bxmd+ma+AugZMmSzJw5k+nTp7Nv3z7Wrl2ri7j06unTp8yePZtu3bqpMlUvLy8uX75MZGQkcXFxlClThnfv3mFiYkKVKlW4fPkyb9++zfAP8IiICF6+fImRkRFWVlZqTXnW1taMGzcOOzs71Ru8dOnS3Lt3L0NjTC7WwoULkzt3bn766acseV7fvn2LsbFxsglHVjuvibHGxcVhYmICvL9e//nnH16/fp1lzuuzZ884cuQIUVFRFClShEaNGlGkSBEADA0Nadq0qVr5qKgo8uTJQ4kSJTIsxpRiTRzlk9ic7u3tTf78+YH///KMjY3F3t4+w1u3wsPDMTc3R6FQ8PHvw6x2Xj+MNSEhQXWumjRpAmSt8xoUFMTFixeJiorC2tqaRo0aqUalZbXzmlOkO+FIZGBgQJs2bWjTpo2uqtSLx48fM3nyZOLi4ti5cyeVKlUid+7cAFSrVo1Tp04RFRWl9sV+5coV8uXLh5mZWYbG6ufnx/z58zE0NOTp06d4eHiohpIpFApsbGxUZRPfrOHh4djb26ve7Bn1Jk4u1o4dO2JkZJQlz+v69etJSEggNDSU1q1bU61aNdWXY1Y7rx/HWqVKFaytralRowYnTpzIEufV39+fCRMmUL58eczMzNi/fz9Xr16lUaNGNGnSRPWr8cPkbs+ePYSEhKiNesvMWBs2bEjTpk0xMDDA2tqa33//nWrVqmFkZMTFixe5evUqU6dOzdDbVP7+/owcOZIaNWowZswYVdLx4fX38ePMPK8fxprckHIbGxu2bduW6efVz8+PmTNnYmtry5s3b/Dz8yMoKIhOnTqp4sgq5zUn0fnibVnZ48eP+eGHH2jVqhVly5bFx8eHgQMHUq5cOdUbIyQkhDFjxlCyZEmqV6+On58fZ8+eZfLkyRk6WU3iULLGjRvj6urKP//8w9q1a1m2bFmyPaXj4+PZunUrR48eZerUqWpfmlkh1pCQEL7//ntKly6dqec1MDCQiRMn0qBBA8qWLcvTp0/Zvn07tWvXpm3btpQpU0atfGaeV01iDQkJYfTo0ZQqVSrTzmtcXBxLlizB1NSU/v37A+9bu1avXk1ISAj169enZcuWqg/y27dv89dff/Hnn38yadKkDP3FmFqsdevWpVWrVkRGRrJu3Tr+/vtvbGxsMDQ0ZNCgQRQvXjzDYg0PD2fu3LkYGhry7NkzSpcuzffffw8k/TKEzD2vqcUK75P34OBgNmzYwJkzZzLtvL58+ZJp06bxxRdf8PXXXxMTE8Nff/3Fvn37GDt2rOqHR6LMPK85jc5aOLK6hw8fMmHCBNq3b0+XLl1QKpVs2LCBffv2Ua5cOdWHoYWFBUOGDGHjxo0cP34cCwsLpkyZgoODQ4bGmziU7OuvvwagRYsW/PPPPwQHBxMcHIy5ubnqjXHjxg2OHj3KjRs3GD9+fIZ+KWoT67Bhw9iwYUOmnteDBw9SqVIl1T3YunXr8vLlS/7++28UCgUdOnTA0dERgJs3b3LkyJFMO6+axjp06NBMvV6NjIx48+aN6vZNQkIC+fPnp1+/fvj4+PD3339jZWWFk5MT4eHh+Pn58ezZM6ZOnZrhf//UYj179iy2trZUr16doUOH0rJlS0xNTbGwsMjwUT/37t2jUKFCtGjRgvj4eBYuXMicOXP4/vvvVbcsEj+3IiIiMvW8ahKrQqHA0tKSoUOH4urqipmZWYaf14SEBM6cOUPRokXx8PBAoVBgZmaGo6MjERERvHv3Tq18Zp/XnOazSTjOnj1Ly5Yt6dKli+qN6u7ujq+vL7du3aJ8+fIolUoMDAyoUqUKlStXJiYmBkNDQ9V984wUHR2NUqkkKiqKvHnzsn37di5fvkxERATBwcEUL14cd3d3KlSoQK5cubC0tGTSpEnY2tpmyVg7deqUJc5raGio6gMuOjoaMzMzihYtSvny5fHz8+P8+fM4Ojqq+kNk5nlNLdZz587h6OiYqec1ISGBhIQETExMCA0NBVANgcybNy89evRg9uzZHDp0CCcnJ8zNzWnQoAENGjRQ3crMarHu37+f6tWro1AoMnUOlgoVKmBsbEz58uUBGDp0KAsXLmT27NmMHj0aAwMDVUtH/vz5ady4caacV01j/fB2ZGK5jGZgYEDJkiVRKpWq86RUKnFwcMDMzIzXr1+rlc+fP3+mXa85Udabr1lPunTpQo8ePQBUvwpKlSpFbGwst27dUiub+CY2MzPLlC9FAFtbW27evMnKlStZtGgR27ZtY+TIkUyePJnvv/+ed+/ecfnyZQAcHR3x9PTMlC9FTWJ9+/Ytly5dUn3gZOZ5LVSoEOfPnycmJgYzMzPCw8PZvXs3bdu2pVWrVuzdu5fg4GBMTEwoVapUpp5XTWMFMu28GhgYYGRkRLt27bh06RJ79uwB3rckxMXFkS9fPvr06cP169d5+PAhALlz586UD29NY71x44Yq1syUJ08eqlWrpnpcsWJFhg0bxr1795g9ezbw/u9+6NAh7t69i4mJSaZ9KWoSq4GBAUeOHOHu3buZEmOicuXK0aFDB+D/P+uNjN7/9v6whePq1askJCRk2vWaE302LRyJScaHzZBFixalZcuW7N27ly+++EL1xZIVJjFzdXUlISGBd+/ecf/+fZo0aaIaUVOmTBmsrKy4fft2lljnI7VYra2tuX37dqbGmKh169bcv3+fPn36ULFiRW7dukWDBg2oUqUKJUqU4Pfffyc4OBhLS0uATB2mp22sGSU4OBg/Pz/CwsKoXr06ZmZmlClTBk9PTzZs2ICRkREtW7ZUfYgnJCSoRi1ltOwaa40aNcidOze5cuVSvccTWwaGDRvGggUL+Pnnn7GwsODQoUMsWrRIYk0l1tDQUJycnDAzM8PIyEhtCGxsbKyq9Qtg8+bN7NixI8U+cyJtcmzC8ezZM06cOEFoaCgODg5UqVKF4sWLJ+k5XblyZU6dOsXt27extbXNlC/wj2OtVKkSjo6OtGzZEgAfHx/VG+HDzmL29vZJhslJrCnHWq1aNRwcHJgwYQIHDhxAqVTSsGFDGjZsCLz/YMqVK1emfNlkl1ifPHnC9OnTKViwIEFBQWzbto169erRqlUrOnTowLt37/Dx8SE0NJQmTZpgZmbG+fPnSUhIyPDRSDkh1hYtWlCkSBHV55KBgQEVK1Zk8ODBzJgxgzx58uDl5YWVlZXEqkGs27dvV4v14/lCcuXKxY4dO9i3bx8zZ86UZEPHcmTC4e/vz8SJE6lcuTL58uVj7969/PnnnzRr1gxXV1e1pKNUqVI4Ojry+++/q4bEZYVYmzZtSosWLYD3M2Hu3LmT6tWrY2hoyOXLl7l48SLTpk3L0F/g2T3W06dP4+rqSvPmzenQoUOS5PKvv/4iV65cGT6FeXaJ9fXr16op4Dt06EDevHnZtm0bV69eZc2aNfTu3ZsuXbpgbW2Nj48Pp0+fxsTEhLdv3zJ69OgM7RyYE2K9du0agYGBdO/eHWtra1UCn5CQwN9//42JiQnTpk3Dzs5OYk1HrIaGhhgaGpI3b15WrFjBkydP+PHHH9VmHha6keMSjpiYGNavX0/Tpk355ptvAHB3d+f7779n27ZtREVF4eHhoTYvQJMmTXj8+DFhYWFYWFhkiVi3b99OVFQUHTt2xMPDA39/f6ZMmYK1tTW5cuVi0qRJGfrmzSmxbtmyhcjISNU1AHDr1i3+/vtvTp48ydSpU1UTE0ms6qKjo4mMjKRy5cqqUR6dOnWicOHCHDt2jI0bN9K9e3caN25MlSpV8PPzQ6FQYGdnl+HLHeSkWLds2UK3bt1Un023b9/m/v37TJkyJUPfVzkx1u7du2Nubk5UVBQvXrzg7du3zJo1S0aj6EmOSzgUCgVRUVGqcd1v377F0tKSypUr8/r1a/79919KlCih+gUO7ztdTpo0KcOHvaUW6+XLl3F0dKR69eoMGzaMFi1akCdPHgoUKCCxpiPWD6+BxPIJCQlMnz4de3t7iTUFBgYGmJiYEBYWBvz/tPDOzs7ExsZy4MABrly5grOzMwULFszU5uicFOvBgwe5evUqzs7OwPvPq4kTJ2ZosplTY028BvLmzUufPn1wdHTMtE7in4McNUpFqVQSExNDaGioathbrly5CAkJwd/fn0aNGhETE8O5c+fUnmNmZpbhX4ppibV8+fI4ODhIrDqOtVy5cnTv3j3Dv8CzU6wABQsWxNramn379vH69WsMDQ2Jj48HoFmzZtjY2HD48OEMjys5OSnWokWLqmJVKpWYmppmyhd4To4VoH79+pJs6FmOSDgSEhKA97/+ChQogLu7Oxs2bGDZsmVs3ryZESNGULZsWZydnenYsSPXrl1TLd6W0SNS0hNrRvscYk38AMrouSuyQ6wxMTGq1XIT/e9//+PNmzfMnz+fuLg4tX45VatWRalUEhcXlyHxfY6xZvTn1ecQa2xsLECmj/b7HGT7WyrPnj3j0qVLNGjQQHXf0NXVFVNTUw4dOkRYWBgeHh60b98eeD8Fb968ecmbN2+GvyEkVok1u8Tq7+/P+vXriYiIIDw8nG+++Yb69euTP39+hgwZwvz585k+fTr9+vXD0tISExMT7t+/j5mZWYaPRpJYJdb0xCoyTrZeS+X58+eMHz+e169f0759e9q0aaPWfPfu3TsUCgXGxsaqbWvXriU8PJxBgwZhbGycYR/iEqvEml1i9ff3Z/LkyTRq1IiSJUvy8OFD9u/fz8yZM1XrSPj5+bFo0SLVirTm5ubcvHmTqVOnZui6GBKrxJqdYv3cZduEIyYmhnXr1qFUKnF0dGTdunW0bduWdu3aqS1/nPgBHRAQwOHDhzl+/Dg//vhjhvZCllgl1uwSa1RUFAsXLsTGxoZevXqptieuI9GrVy+1OA8cOEBoaCgmJibUq1cvQ9ebkVgl1uwUq8jGt1QMDAxwdHQkX7581KtXj/z587Nw4UIA1Yd44kUWHR3N1atXefz4caYswCOxSqzZJda4uDhev35NnTp1gP+fmdfKyorIyEjg/0fKGBgYqCZ8ywwSq8SanWIVgDIbi46OVnv8119/KTt37qz85ZdflBEREUqlUqmMj49XhoeHK+Pi4pSRkZGZEaZSqZRY9UVi1b1nz56p/j82NlapVCqVvr6+ysWLF6uVe/Pmjer/ExISMia4j0is+iGxCn3Iti0cAKampgCq0Sb16tVDqVSyaNEiFAoFrVq1Yvfu3QQFBTF06FDV5C8Sq8QqsaasaNGiqjg/XGPk1atXqjI7duzA2NgYNzc3DA0NM239IYlVYs1OsX7usnXCkShxmeaEhATq16+PQqFg8eLFXLx4kRcvXjBz5sxMW530YxKrfkis+osz8cM5cdigr68vv//+O7NmzcrUxe0+JLHqh8QqdCnHDDxWKBQoFAqUSiX16tWjfPnyREREMGvWLFVP5axCYtUPiVX3lP/1KTcwMKBQoULs2rWLXbt24eXlleV690us+iGxCl3JES0ciRI7B/3yyy/cuHGD2bNnZ9k58SVW/ZBYdSvxV6KRkRFHjx7FzMyMadOm4ejomMmRJSWx6ofEKnQlx7RwfMje3p5Zs2ZRrFixzA4lVRKrfkisulW1alUApk+fnuVX0ZRY9UNiFemVbefh+JQP7+NldRKrfkisuhcTE6Pq+JrVSaz6IbGK9MiRCYcQQgghspYceUtFCCGEEFmLJBxCCCGE0DtJOIQQQgihd5JwCCGEEELvJOEQQgghhN5JwiGEUPnpp5/o2bMnwcHBSfZFRUXRr18/Jk6cSEJCQiZEJ4TIziThEEKo9O/fH0NDQ1asWJFk35o1a4iOjmbgwIGqGR2FEEJT8qkhhFAxNzenT58+XLlyhcOHD6u2nz9/nr/++otu3bphbW2t1xgSEhKIjY3V6zGEEBkvR62lIoRIv3r16nH+/Hl+/fVXqlatipmZGatWraJKlSq4urry4MEDtm3bxu3bt3n37h22trZ06NCBevXqqeqIiIhgy5Yt3Lhxg+DgYHLlyoW9vT2dO3emfPnyqnJBQUEMGjSIr7/+mri4OI4dO0ZISAhjx46lWrVqmfDqhRD6IjONCiGSiIqKYuTIkdjY2JA/f36uXr3Kzz//TGBgIDNnzqR06dK0aNGC3Llzc+bMGU6cOMGAAQNo3LgxAM+ePWPfvn1UqFABc3NzYmJiOH/+PMePH2fSpElUrFgR+P+Eo2DBghQtWlRVp7W1NUWKFMnEMyCE0DVp4RBCJJE3b1769++Pl5cXAIMGDaJQoUJMnz4de3t7Jk2ahKGhIQDVqlUjIiKCTZs20ahRIwwMDLCxsaFv376q+hISEqhatSovX75k//79qoQjkbGxMRMmTMDISD6ShMip5N0thEhW9erVKV26NFFRUTRq1Ijnz58TEBBAt27dAIiPj1eVrVGjBv/88w/Pnj3Dzs4OgEOHDnH06FH8/f3V+mTY2tomOVbNmjUl2RAih5N3uBAiRcbGxqpEIDw8HIBff/2VX3/9NdnykZGRAOzZs4dffvmF5s2b4+npSb58+TAwMMDX15eAgIAkz7OwsNDPCxBCZBmScAghNJI/f34AOnTowBdffJFsGRsbGwBOnz5NxYoV+fbbb9X2R0dH6zdIIUSWJQmHEEIjNjY2FC1alCdPntC1a9dUy398i+TJkyfcvXsXS0tLfYUohMjCJOEQQmjs22+/xcvLixkzZuDs7EzBggWJiooiICCAR48eMWLECACcnJzYvn07W7ZsoUKFCjx79oxt27ZRpEgRmaVUiM+UJBxCCI1VqlSJmTNn8vvvv7N+/XqioqLIly8fdnZ21K1bV1XOw8ODt2/fcuzYMXbu3ImdnR3ffvst58+f5+bNm5n4CoQQmUXm4RBCCCGE3snU5kIIIYTQO0k4hBBCCKF3knAIIYQQQu8k4RBCCCGE3knCIYQQQgi9k4RDCCGEEHonCYcQQggh9E4SDiGEEELonSQcQgghhNA7STiEEEIIoXeScAghhBBC7yThEEIIIYTe/R9evRI2C49p/gAAAABJRU5ErkJggg==", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEGCAYAAABlxeIAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABRo0lEQVR4nO3deXxM1/vA8c9kJYglocgklsS+74029lpib1FqKUGlaGlrCaW2atGqtVTtpfalaieWRO3EUlIkBJFQElKSkGVmfn/km/kZ2W5iJpPE8369vF6ZM3fu89xMzDP3nnPPUel0Oh1CCCHeeBbmTkAIIUTOIAVBCCEEIAVBCCHE/0hBEEIIAUhBEEII8T9SEIQQQgBgZe4EchqdTsejR49ISEgwdypCiBzG2tqa4sWLo1KpzJ2KSajkPgRDDx8+JDExEWtra3OnIoTIYRISErCysqJEiRLmTsUk5JLRKxISEqQYCCFSZW1tnaevHkhBEEIIAUhByJFKlixJ8+bN8fDwoFmzZixevBitVmvutIxGp9Mxfvx4GjZsSNOmTbl8+bK5UzKaEiVK8M033+gf//zzz8yaNcuMGZnfrFmz+PnnnwGYMWMGfn5+ZsmjXr16REZGmiV2biGdyhn48ssT3Lz51Gj7c3W156efGqe7Tb58+Thy5AgAjx49wtvbm6dPnzJ27Fij5WFsUVFRFClSRNG2hw4d4tatW5w+fZrz588zZswY9u3bZ/Scjhw5QlRUlNH2V6RIEZo3b57uNra2tuzZs4cRI0bg4OBgtNhKaDQaLC0t03ycE/j4+Jg7BZEOKQgZuHnzKSdP/mu2+MWLF2f27Nm0adOGMWPGEBoayrBhw4iNjQXg+++/p2HDhhw/fpxZs2ZRvHhxrly5Qvv27alSpQpLly7lxYsXrFq1inLlyrF//37mzJlDfHw8RYsWZfHixZQoUYKIiAi8vb158uQJtWvX5siRIxw8eBAHBwc2b97MsmXLiI+Pp27dusyaNSvFB824ceN48OABvXv3pkOHDuTLly/NY9q7dy89evRApVJRv359/vvvP/7991/eeusto/7uoqKiuH//vlH3mRFLS0v69u3LkiVLGD9+vMFzoaGhjBw5koiICBwdHZk3bx5qtdpgm+joaMaPH8+lS5cAGDVqFB07dmTbtm3MmzcPnU5Hq1at9GchZcuWxdvbm6NHjzJ58mR69uxp8Dg0NDTFewcwcuRILl26hEqlolevXnh7e9OlSxeqV69OQEAA0dHRzJ07l7p16/LkyRNGjBjBnTt3sLOz48cff6RatWrMmjWLsLAw7ty5w7179xgyZAiDBw8GYM6cOWzatInSpUvj6OhIzZo1Afjss89o3bo1HTt2pF69evTo0YMDBw6QmJjIsmXLqFChQrp/i8lWrlzJ3bt3mTRpEgAbNmzg0qVLfP/99/Tr14/w8HDi4uIYPHgw/fr1M/gd3717lz59+uDv7w8kncXFxMQwZswYQkJC8PHxITIykvz58/PTTz9RoUIFo/xt5AZyySgXKFu2LFqtlkePHuHo6MjmzZs5dOgQS5cu5euvv9Zvd/XqVaZPn46fnx+bN2/m5s2b7N+/n969e7N8+XIAGjVqxN69ezl8+DBdu3Zl4cKFAPz44494eHhw6NAh2rdvz7179wC4ceMGO3bsYNeuXRw5cgRLS0u2bNmSIsfFixczefJkzp49i4eHB+PGjePKlSupHs+DBw8oXbq0/nHp0qWz/YPblLy8vNi6dStPnxqeWY4bN47u3bvj5+fHBx98kKJgAPz000/Y29vj5+eHn58fHh4ePHjwgGnTprF161YOHz7MxYsX2bNnDwCxsbFUqVKFffv28fbbbxs8LlasWKrv3ZUrV7h//z7+/v74+fnRq1cvffzY2Fj27NnDzJkzGTlyJJB0yadGjRr4+fkxfvx4hg8frt8+KCiIjRs3sn//fn788UcSEhK4dOkSf/zxB4cOHWLVqlVcuHAhzd+Vg4MDhw4d4uOPP2bRokVA2n+LL+vYsSO7d+/WP/7jjz/o3LkzAPPmzcPX15cDBw6wbNkyHj9+nNFbpjdq1Ci+//57fH19mTx5co4+KzcFOUPIJZJHBycmJuLj48PVq1exsLDg1q1b+m1q166t/5ZdtmxZmjVrBkCVKlX466+/AAgPD2fw4ME8fPiQ+Ph4XFxcADh9+jSrVq0CoEWLFvrLP8eOHePSpUu0bt0agBcvXuDo6JhqjrVq1aJWrVq8ePGC3377jbZt2/L111/z6aefpnosL8tL47oLFSpE9+7dWbp0qcGZ0rlz51i5ciUA3bt3Z+rUqSle6+/vz6+//qp/XKRIEfbu3Uvjxo31v/cPPviAkydP4unpiaWlJR06dNBv//LjtN67Nm3acOfOHcaNG8d7772n/zsB6Nq1KwDu7u48e/aM//77j9OnT7NixQoAPDw8ePLkib7Yvffee9ja2mJra4ujoyOPHj3i1KlTtGvXDjs7OwDatGmT5u+qffv2QNLfTvIHfFp/iy9zdHSkTJkynDt3jvLlyxMcHEyjRo0AWLp0qb5ghoWFcevWLYoVK5ZmDsmio6M5e/YsAwcO1LfFx8dn+Lq8RApCLnD79m0sLS0pXrw4P/zwA8WLF+fIkSNotVqcnZ3129na2up/trCw0D+2sLBAo9EAMH78eLy9vWnbti3Hjx/nhx9+AFL/kE5u//DDD5kwYUKGeSYmJuLr68v69esJCQlhzJgxdOvWLcV2pUqVIjw8XP84PDyckiVLKvhN5B5DhgyhVatW9OzZM81tUiuCqb0P6d0qZGtra3D57uXH6b13R44c4ciRI6xYsYIdO3Ywb968VHNSqVTpFnAbGxt9m6WlJYmJiWkeW2qSX29paan/G1V6a1SXLl3YsWMHFSpUwNPTE5VKxfHjx/H392fPnj3Y2dnRpUsX4uLiDF5nZWVlMEgj+XmdToe9vb2+/+5NJJeMcriIiAhGjx6Nl5cXKpWKZ8+e8dZbb2FhYcHmzZv1/4mUevr0KaVKlQJg48aN+vZGjRqxY8cOwLAz1sPDg507d/Lo0SMAnjx5QmhoaIr9Ll68mLfffptdu3bh7e2Nv78/n3/+OcWLF0+xbdu2bdm0aRM6nY5z585hb29v9P4DcytatCidOnVi3bp1+rYGDRqwfft2ALZu3UrDhg1TvK5Zs2b6b+OQ1A9Sr149Tp48SWRkJBqNhm3bttG4cfoDEyDt9y4yMhKdTkfHjh0ZO3aswSivP/74A4BTp05hb2+Pvb097u7ubN26FYDjx49TrFgxChUqlGZcd3d39uzZw/Pnz4mOjubAgQMZ5vqytP4WX9W+fXv27t3L9u3b6dKlC5D0912kSBHs7OwICgri/PnzKV5XvHhxIiIiePz4MXFxcRw8eBBIOrNzcXHhzz//BJIKRFqXPfMqOUPIgV68eEHz5s1JTEzE0tKS7t276y+7DBgwgAEDBrBz507eeecd/Wm5UqNHj2bgwIGUKlWKevXqcffuXSDp2qm3tzc7duygcePGvPXWWxQsWBAHBwfGjRtHjx490Gq1WFtbM2PGDIMzE4CqVaty5MiRdD8okrVq1QpfX18aNmyInZ2d/ttpXvPpp58afLhPnz6dkSNH8vPPP+s7lV/1xRdf4OPjQ5MmTbCwsGDUqFF06NCBr7/+mvfffx+dTkfLli1p165dhvErVaqU6nuXL18+RowYof+W/PIZRJEiRfD09NR3KkPS38yIESNo2rQpdnZ2LFiwIN24NWvWpEuXLrRo0QK1Ws3bb7+t5Nell9bf4quKFClCpUqVuH79OnXr1gWSLjGtXr2apk2b4ubmRr169VK8ztramq+++oq2bdvi4uKCm5ub/rnFixczZswYfvrpJxITE/Ud7W8KmbriFWFhYQanweYYdmoOcXFxWFpaYmVlxdmzZxkzZkyuP3U2x7DT3KxLly5MnjyZ2rVrmzWPnP63GB8fj5OTk7nTMAk5Q8hATvzwNoWwsDAGDRqEVqvFxsaGn376ydwpvba8/OGdl+XFv8XcQs4QXvHqGYIQQrwsL58hSKeyEEIIQAqCEEKI/5GCIIQQApCCIIQQ4n9klFEOVLJkSapUqaK/D+HDDz9kyJAhWFjknfo9b948fv/9dywtLZk+fTotWrQwd0pG8Sa8dyLvkoKQgRvRGp5rjDcQK7+liooF05+SODdOf50Z169fZ/v27Rw7dowHDx7QrVs3Tp06lf1TNWs02Bw6hNXff5NYowbxLVvCa+aQ1987kbdJQcjAc42Op5mbHSIDmSsuuWX668zYt28fXbt2xdbWljJlylCuXDkCAgJo0KBBlveZaRoNhXv0wCogAFVsLDo7OxLr1uW/TZteuygky4vvncjb5Dw2F8gN019PnDiR5s2bp/g3f/78FNvev38/xfTXDx48MOrvLCM2hw5hFRCARUwMKp0Oi5gYrAICsDl0yKhxcsN7J0QyOUPIJXL69NfTpk3L9LGYk9Xff6P63zf1ZKrYWKyuXCH+f8dqLDn9vRMimRSEXCA3TH89ceJE/QfXy7p27crnn39u0Fa6dGmzT3+dWKMGOjs7VDEx+jadnR2JRp7ILDe8d0Ikk0tGOVxumf562rRp+jn2X/73ajGApAVTtm/fTlxcHHfu3OHWrVv62SqzS3zLliTWrYu2QAF0KhXaAgVIrFs3qWPZSHLLeydEMjlDyIFy4/TXmVG5cmU6d+7Mu+++i5WVFTNnzsz+jk5LS/7btCmpL+HKFRKrVzfKKKO8/t6JvE0mt3vFq5PbmWPYqTnk9CmHRdrkvcteeXlyOzlDyEBO/PA2BZlyOPeS904Yi6IzhKdPnzJ58mT8/PyIiIgw6MRKPm3NK2T6ayFEevLyGYKiTuWhQ4cSEBDAN998w+PHj1mwYAEuLi588cUXps5PCCFENlF0hlCiRAn++ecfHBwcKFKkCFFRUYSFhdGxY0cCAgKyI89sI2cIQoj0vPFnCFqtlsKFCwNQsGBBoqKiKFWqFMHBwSZNzhysra1JSEgwdxpCiBwoISEBa2trc6dhMoo6lWvVqoWfnx8tW7bEw8ODYcOGUbBgQSpWrGjq/LJd8eLFefToEfHx8eZORQiRw1hbW1O8eHFzp2Eyii4Z3bp1C51Oh6urK48ePWLcuHE8e/aMSZMmUbVq1ezIUwghhInJfQhCCCEAhX0IOp2OpUuX0qJFC2rWrAmAv78/mzZtMmlyQgghso+igvDNN9+wfPlyPvnkE/19B2q1mpkzZ5o0OSGEENlH0SUjZ2dnLly4gKOjI0WLFuXJkyfodDqKFSvGkydPsiNPIYQQJqboDEGj0VCwYEEAVCoVANHR0fo2IYQQuZ+igtCuXTu+/PJL4uLigKQ+hYkTJ9KxY0eTJieEECL7KCoIc+bM4f79+xQuXJj//vuPggULcufOHelDEEKIPCTDG9M0Gg1btmxh/fr1PH36lDt37uDs7JztK1wJIUROt2TJEurVq0f9+vXNnUqWKOpUTp6/KK95eRnH7OLo6EhERITEzaNxzRlb4po3bunSpbM9F2NTdMmoY8eO7Ny509S5CCFErrZkyRLOnTtn7jSyTNFcRi9evKBbt264u7vj7OysH2kE8Ntvv5ksOSGEyE2GDBli7hRei6KCUL16dapXr27qXIQQIufSaLA9fBiLkBBsy5UjrkWL116DO6dRVBAmTZqUante7FcQQogUNBocPvoI6wsXUMXGUtTOjoQ6dYhct86gKOT2TmVFfQgv02g07Nq1i+7du1OqVClT5CSEEDmK7eHDWF+4gEVMDCqdDouYGKwDArA9fNhguyFDhuTaYgAKzxAALl68yOrVq1m3bh0RERH06tULf39/U+aWpyxdujRF2+DBg82QiRA5n5PTuhRtYWEfmSGTJNZXrqCKjTVoUz1/jvXVq8S9956ZsjK+dM8QHjx4wOzZs6lRowYNGjTg2rVr/PDDDxQrVoyffvqJBg0aZFeeQghhNgnVq6OzszNo0+XPT0K1ambKyDTSPUNwdnamSJEifPPNN3z44YeUKFECgLFjx2ZLckKkJad9gxR5W1yLFiTUqYN1QACq58+TikHdukkdy3lIugWhd+/ebNu2jR9//JHw8HA++ugjatSokV255QhyqUcIgaUlkevWYXv4MIVv3+a/smVTHWWU2zuV0y0Iq1atYtGiRWzZsoXffvuNWbNmUbVqVZ4+fUpERIT+jEGInELOHITJWFoS9957aB0diUvjDuncfh9ChqOM7Ozs6NevH76+voSEhNCzZ0+cnZ2pXbs2PXr0yI4chRBCZINMDTt1cXHh66+/5tq1a/j5+eHg4GCqvIQQQmQzxcNOX+Xu7o67u7sxcxEKSJ+GEMJUslwQhBBCmM7du3cZPnw4jo6OVKxYER8fH5PHzPSdykIIIbLGy8uLEiVKpJgbbt++fVSqVAk3NzdmzJgBwI0bN2jfvj0rVqwgMDAwW/LLsCBotVoOHz5MfHx8duQjhBA5jpPTOv0/W9v5qY5mU6J///7s27fPoE2j0TBs2DD27t1LYGAg69evJzAwkDp16rBhwwZatGhB8+bNjXEYGcrwkpGFhQWdO3fm2bNn2ZGPEELkWi9f1mnVqhWtWrUyeL5Jkybcvn3boO3MmTO4ublRvnx5AHr27MmOHTuwtrZmypQpNGnShG7dujFgwACT56+oD6FJkyacOnWKt99+29T5CCFErlWuXLlM35gWFhaGs7Oz/rFareb06dN4e3szefJk1q1bR9myZU2QbUqKCkKZMmVo164dnTt3TrFAztSpU02WnBBC5CZZuTEttVWMVSoV1atXZ8uWLcZISzFFBeH58+d06dIFgHv37pkyHyGEMJtTz1Qp2t4ulOGy83pZmbpCrVYTGhqqf3zv3j2zrc+sqCCsXLnS1HnkGXKfgBAiMxo0aEBQUBAhISE4OTmxYcMG1q3LWqf161I87PSff/5h2rRpDB8+HIDr169z+fJlkyUmhBC5TUYL5PTq1Qt3d3euX7+OWq1m+fLlWFlZsXDhQtq0aUOVKlXo0aMH1cw0rbaiM4TNmzczdOhQPvjgA9atW8fChQt59uwZPj4++Pr6mjTBiIgIfv75Z6KiolCpVLRq1QpPT0+DbXQ6HStXruTChQvY2toydOhQfY+9EEJkl4wuGa1fvz7Vdk9PzxSfa+agqCB88803HDx4kNq1a7Nx40YAatWqxaVLl0yaHIClpSV9+/alfPnyPH/+HB8fH2rWrIlardZvc+HCBR48eMD8+fMJCgpi2bJlfPfddybPTaQts5fO5FKbyAvy/GynAA8fPqRWrVoA+hFGKpXKYLSRqRQtWlT/bT9//vw4OTnx+PFjg23OnTtHkyZNUKlUVKxYkZiYGJ48eWLy3IQQ4mVLlizh3Llz5k4jyxSdIdSrV481a9bQr18/fduGDRto2LChyRJLzcOHDwkJCcHNzc2g/fHjxzg6OuofOzg48PjxY4oWLZqt+QnxumQ9h9wtt58hKCoI8+fPp3Xr1ixfvpyYmBjatGnDjRs3OHDggKnz03vx4gWzZ8+mf//+2L26tmka43hf5evrq+/zmDFjhkERyYzMvi6t7bMaPyv7sbKyMlq8rDLW7y2z22f3cRv7d610X+Z6j7Mr7qsxTBL3WWSGcTNqz80UFYTKlStz7do1du3aRYcOHXB2dqZDhw4ULFjQ1PkBkJiYyOzZs/Hw8KBRo0YpnndwcCDipRWMIiMjUz07ePVW8og0Vj3KSGZfl9b2WY2flf04OjoaLV5WZeX3lplx4ab+PStl7N+10n2Z6z3OrrivxjBN3JR/b0r/rkqXLp3rl9BU1Ifw+eefY2dnR48ePRg9ejQ9e/akYMGCjBw50sTpJX37/+WXX3BycqJDhw6pblO/fn38/f3R6XTcuHEDOzs7uVwkhMh2GQ07zekUnSGsWrWK+fPnp2hfs2YNc+fONXZOBq5fv46/vz8uLi6MHj0aSBrLm1ydW7duTZ06dQgICODzzz/HxsaGoUOHmjQnkXfINfvsIb/n3CHdgrBixQog6ZJN8s/Jbt26lS3X0CpXrsymTZvS3UalUjFo0CCT5yKEMK+cXlhy+yWjdAvCmjVrAIiPj9f/DEkfwG+99RarV682bXZCCJGL5OlRRkeOHAFgwoQJfPvtt9mSkBAi+73upG4ib1DUh/ByMdDpdAbDPC0sZBVOIYR55fRLSbmFooIQHh7OsGHD8Pf3JyoqyuA5jUZjiryynUydIIR4XebqQ0hISODUqVNcunSJqKgoihQpQq1atXj77bextrZWvB9FBWHIkCHY2dlx6NAhmjZtir+/P5MnT84RkzEJ85JC+mYyuMT0LBJQySUmsr8PISIighkzZrB69WqKFStG5cqVKVSoEM+ePWP+/Pk8efKEjz/+GB8fH2U3sCoJeuLECe7evUuBAgVQqVTUqlWL5cuX07hxY/nPL4QQZuLh4cHAgQO5ePEiTk5OKZ4PDw/n999/p0mTJgQGBma4P0UFwdLSEiurpE2LFCnCo0ePsLe3JywsLJPpC2FaFmhpxzXqEMYFnNhLZXOnJITJXLp0CRsbmzSfL126NKNHj2bEiBGK9qeoIDRq1Ig9e/bQtWtX2rRpw4cffkj+/Plz7VhbY5BLJTmQRsN+ltKIu9gRTyw2nMYFNL3A0tIoIWQ0jshJXi0Gz58/1/ch5M+fP83t0qJoiNCaNWto2rQpAHPnzqVFixZUr17dbMu8CZEa28OHacRdChGPJVCIeBpxF9vDh82dmsmceqbS/9sVEplqwRK507Fjx/D29mbQoEE0btw43W2PHDlCw4YNKVSoEGq1mkKFCtGwYUMOHTqUqZjpFoQqVarg7e3Nnj17iI2NBZLWJJgwYQIzZ86kVKlSmQomhClZX7mCHfEGbXbEY331qpkyEsKQl5cXJUqUoHr16gbt+/bto1KlSri5uTFjxgwgqX/gl19+oUOHDnz88cdp7vPcuXN4enrSqFEjDh48SGBgIAcOHKBhw4Z07NiRs2fPKs4v3YIwYcIEAKZPn06ZMmVwdXVlwIABrFy5kps3byoOIkR2SKhenVgMT41jsSHBTOvTCvGq/v37s2/fPoM2jUbDsGHD2Lt3L4GBgaxfv96gA3jdunX06tUrzX3+8MMPjBkzhgULFtC8eXMqVapEixYtWLhwIWPHjuWHH35QnF+6fQi9e/emd+/eQNLwpr/++gt/f38WLVrEkCFDeOuttwgNDVUcTAhTimvRggu4pOhDqNqihdlySqvPQfoi8iYfHx/9z69Otw/QpEkTbt++bdB25swZ3Nzc9CtD9uzZkx07dlC1alXu3r1L4cKFsbe3TzPmyZMnmTNnTqrPDR48OFMLmSnqVIakucddXV25d+8e9+7d486dO+kmKUS2s7SkDYNpxzVqE85FSrOXyoRmoUNZPrBFViRf7smMsLAwnJ2d9Y/VajWnT58GYPny5QwYMCDd10dFRVG6dOlUnytdujT//fef4lzSLQjnzp3Dz88PPz8/zpw5Q/ny5XnnnXfo06cPv/zyC8WKFVMcSKRORisZlxYLdlOV3VQ1dyriDZSVO5XTW/FxypQpr51TaqtHpiXdgtCwYUOqVKnC2LFj2bRpE/ny5Xvt5ITIi+TOXZFVarXa4NL7vXv30vzGn5qYmBhcXFxSfU6n0+kHBCmR4fTXx44dY+bMmUyaNIl3330XDw8PPDw8qFKliuIgQoicIaddCpNJ6aBBgwYEBQUREhKCk5MTGzZsyNSQ/sNGHFadpU7lZcuWERoayttvv8327duNlowQQuRmGc1l1KtXL44ePUpERARqtZopU6YwcOBAFi5cSJs2bdBoNHh5eVEtEyPjku8RM4YsdyqHhISwZ88eoyUihEjy6vQbaDRGu9NamFZGfQjr169Ptd3T0zPLk4Xu27cPe3t7/c1rN2/epF+/fly5cgV3d3dWrlyp+J6xdO9DOHfuHLNnz6Zz5844OjrSqFEjtmzZQqVKlVi3bh2PHz/O0gEI8bos0NKeQCZwkPYEYoHW3CkZhQVa9rOU9fzOZA6wnt9x+OijpKIgcrwhQ4Zk+5Q+EydONOg4HjhwIIULF2bdunUUKFCAUaNGKd5XumcIzZs3p3Hjxnh4ePDVV1/RqFEjbG1ts565yLVy1Gio9OYsyuXacU0//QYkTb+hDQjA9vBh4t57z8zZiYyYYz2Emzdv0qBBAwAePnzIX3/9xZ07d3BycqJRo0bUrFlT8b7SLQhRUVFYyqmqyGFenrMI/n/Oovh0OtdyWmdqWuoQlmL6DdXz51hfvWqWgpBbfm85hbnXVD558iTlypXTT4Xt4OBAdHS04tenWxCMVQxCQ0NTrOTz8o0YQmRGWnMW6a5eBUqYJykjuYATsdjoix2ALn9+mX5DpKlBgwbMnz+fQYMGsWzZMtq1a6d/7tatW4oWxkmmuFM5sxISEliyZAlLlizh1q1buLm56VfyCQ4Oply5cnh7e/PJJ58onppVCPj/OYte/tD8/zmLHpkvMSPYS2VOvzL9hk3dusSZcfoNoZw5LhnNmTOHjh07MmrUKCpUqMCSJUv0z61Zs4YmTZoo3pfJCkKtWrVo0aIFS5YsoVGjRgZnGxqNhjNnzvD7779Tp04drspslCIT0p+zaKO503stWixSTL/x67opMsoolzDHJaOqVaty8+ZNIiMjcXBwMHhu5MiRmfrCbbKCcPToUUqUSP303dLSEnd3d9zd3Xn0KP1vdIsWLSIgIIDChQsze/bsFM9fvXqVWbNm6WM1atSIbt26vf4BiJzLiHMW5UQppt/II8clTOPWrVv6n5PnLbK2tsbJyYkiRYpkal+KFsiJi4vj66+/pnz58hQuXBiAAwcOsHDhwjRf83Ix+PHHH1Pd5qeffqJ48eLpxm7WrBnjx49Pd5sqVarwww8/8MMPP0gxeEMkf2hOpxW7qYpW2Z/yG+PlhXOS/wnTW7JkCefOncvWmG5ublSoUAE3Nzf9vzJlypAvXz569+6dqcntFP0v+uKLL7hy5Qq///67frxrtWrVWLx4saIgU6dOTbX922+/zfC1VatWpWDBgoriCCGEOZnjPgStVotGo0Gr1er/xcfHc+XKFeLj4xk9erTifSm6ZLR9+3aCg4MpUKAAFhZJNcTJyYmwsLB0X5c8x4ZGo+HIkSMGs/rdunWLQoUKKU40PTdu3GD06NEULVqUvn37pjmCydfXF19fXyBpmtqMet/Tet5c7WlxdHTk+++/T9E+btw4/c9WVlaZ3q+SuMZoT3f/zyKNEzez+0lre1O3p5VPWsyVZybyT0tO+38EGO99MTMrKysqVqzIkiVLjHcfQjIbGxsSExMN2h49epSiA+NVAwcOBODFixd4eXnp21UqFSVLlmTBggWKE01LuXLlWLRoEfny5SMgIIAffviB+fPnp7rtqwtWREREpLvvtJ7PjvbM3AimZP+Ojo4ZHm9mGet4099/yksdWYub2f2ktb2p29PKJy3mylN5/mkx1/+v9GX9fcnMDKXZxd7e3niznSbr3r07H3/8sX5Vnvv37zNy5Eh69uyZ7utCQkIA6NevH7/99pvipDLDzs5O/3PdunVZvnw5T58+lcV7hBBv/I11GzduzNREeYoKwnfffceYMWOoUaMGsbGxVKhQgcGDB/PNN98oCvJyMdBqDeecSb4ElVVRUVEULlwYlUpFcHAwWq3WaJeihPHlqCkwhDAyc9yH0Ldv3xSL4CQkJHD79m2uX7+eqUlIFV8ymjt3LnPnzuXRo0c4OjpmahWegIAAhg0bxuXLl3nx4gWQtHCDSqVCk8GkXXPnziUwMJBnz57h7e1Njx499JevWrduzalTpzhw4ACWlpbY2NgwcuTITOUmhHjzmOrMwRz3Ibi5uaVos7KywtPTk7Zt22Y4ktPgdUo2+u2336hduzY1a9bU7/zSpUtcvnyZvn37Zvj6jz/+mI4dO7JixQqDSzxKjBw5Mt3n27ZtS9u2bTO1T5E2+QYvRO4yadIko+1L0fWaiRMnphi54+zszIQJExQFuXPnDtOnT6dKlSqUKVPG4J8QQois+fPPP426naIzhNQ6aQsXLkxUVJSiIF27duXAgQO0adNG0fbC9ORMQIicTavVMnHiRJ4+fUr9+vX5+OOPU2yzYcMGxo8fT+/evWnatCmVKlXSzxl348YN/Pz8WLt2LbVr16ZTp04ZxlR0hlC1alW2bt1q0LZ9+3bF6yq/ePGCrl270rp1a/r162fwTwgh3hReXl6UKFGC6tWrG7Tv27ePSpUq4ebmxowZMwDYsWMHYWFhWFtbo1arU93funXrWL9+PWFhYfTt25fixYuTP39+SpQowccff8yDBw/YuHEja9euVZSfojOEmTNn4unpycaNG3F1dSU4OJhDhw4p7r2uWrUqVatWVbStECL3etOHeWakf//+DB8+3ODLsEajYdiwYRw8eBC1Wk2DBg3o1KkT169fx93dnSFDhtCtWzdatmyZ6j5r1Kihn0YoNjZWv8xAZvtrQWFBePfdd7ly5Qrr1q0jNDSUhg0bMm/ePMVrGhiz00MIIXKrJk2acPv2bYO2M2fO4ObmRvny5QHo2bMnO3bswNnZWT9TqdK1aezs7LJUCJIpnu3UxcUFHx+fLAc6ePAgGzZs4OHDh+zcuZNz587x9OlTWsg87yYlfQVCZJ+XPyNfnRkhLWFhYQZfrtVqNadPn2bEiBF89tlnHDt2LFNrGrwOxQXhzz//xM/Pj4iICIM5iZTcgbxgwQLmzZvHoEGD2LJlCwD58+fn888/58SJE1lIWwghcp5y5cpl+sa0lz9Pk6lUKuzs7Fi+fLkx08uQok7lKVOmMGTIELRaLZs3b8bBwYH9+/crnmt77ty5+Pr64uPjo78zuXLlyly/fj3LiQshRF6gVqsJDQ3VP753757Z5kVSVBBWrFjBwYMHmTNnDjY2NsyZM4edO3emuBaWlmfPnulPiZLvIk5ISJClM4UQeUpWpr9u0KABQUFBhISEEB8fz4YNGxQNEU1NQkICx44dY+PGpJUDY2JiiImJUfx6RQUhKipKP0zKxsaGhIQEGjZsiJ+fn6IgTZo00Q+lSjZ//nyaN2+uOFEhhMjpMlogp1evXri7u3P9+nXUajXLly/HysqKhQsX0qZNG6pUqUKPHj0yNSFdsr///puKFSsyePBg/UzTfn5+BjNNZ0RRH4KrqytXr16lWrVqVK9encWLF1O0aFGKFi2qKMiCBQvo2LEjS5cu5dmzZ1SqVAl7e3t27typOFEhhMjt1q9fn2q7p6cnnp6er7XvTz/9lKlTp9K3b1/9Z3PTpk0zNYhEUUH49ttviYxMWjhixowZfPTRR0RHR7No0SJFQUqVKsXZs2c5e/Ysd+7cwdnZmYYNG772TKdCvClkfH/uYI7J7ZJdvXqVPn36AP9/ab5AgQI8f/5c8T7SLAgLFy5k+PDhAFSsWFE/o17Dhg0JDg7OVKIXL17EwcGBhg0b0rBhQwBCQ0N5/PgxtWrVytS+hBBCpFS2bFnOnz9v0IeRfI+DUml+Rf/666/1P9etWzeLKSbp06cPCQkJBm3x8fGKZkoVQojcIqM+BFOaNm0a7du3Z9KkScTHx/P999/TvXt3RWvXJ0vzDMHV1ZWvvvqKatWqkZCQwIoVK1LdTkmHxd27d/V34b28f6WjlIQQIjcw5yWjDh06sHfvXpYtW0bTpk25c+cO27Zto169eor3kWZB2LBhA7NmzWL9+vUkJCSwZs2aFNuoVCpFBUGtVhMQEGBwphEQEJAj1yAVQoisMseKaS+rW7eu4r7d1KRZECpWrMiyZcuApFuwfX19sxzkiy++oHPnzowZMwZXV1du3rzJjz/+aHBZSgghcjtzniGktaSxra0tarWatm3b8tZbb6W7jwxHGWk0Go4fP05cXBy2trZZSnTw4MEUKVKE5cuXExoairOzM7Nnz6Zbt25Z2p8QQghDN27cYPv27TRs2BBnZ2dCQ0M5c+YMHTt2ZOfOnQwdOpStW7emu8JkhgXB0tKSihUrEhkZmaVLPBqNhpYtW7J//366d++e6dcLIURuYc5LRlqtlg0bNtC1a1d9244dO1i3bh2nTp1i9erV+Pj4vF5BAOjduzcdOnRgxIgRqNVqg0XsM5qt1NLSkpCQkFQncBJCiLzEnJeM9u/fn+LGtw4dOuhHc/bp00d/K0FaFBWExYsXAzB58mSDdpVKxa1btzJ8/aRJk/D29mbKlCkpCorcnCaEEK/P1dWVxYsXG3zo//LLL7i6ugIQERFBgQIF0t2HooIQEhLyGmnCoEGDAAxGKul0OlQqFRqN5rX2LYQQApYtW8b777/PzJkzcXJyIiwsDEtLS7Zt2wbA9evXmTZtWrr7ULwewut43YIihBAifXXr1iUoKIiTJ09y//59SpUqhbu7O9bW1kDSJKMZLbSjqCA4OzsbXOZ52d27dzN8fZkyZYCkTo9///2XUqVKKQkLwKJFiwgICKBw4cLMnj07xfM6nY6VK1dy4cIFbG1tGTp0aIqb4IQQIjuY+z4Ea2vr11pdTVFBWLt2rcHj+/fvM2/ePHr27KkoSFRUFEOHDmXLli1YW1sTExPDn3/+yZkzZzK8rbpZs2a0bduWn3/+OdXnL1y4wIMHD5g/fz5BQUEsW7aM7777TlFeQghhTObsVH769CmTJ09OdWVLJV/cQeF6CE2bNjX417NnT7Zv387KlSsVBfH29qZw4cLcuXNHvyiOu7u7fhGH9FStWpWCBQum+fy5c+do0qQJKpWKihUrEhMTw5MnTxTlJYQQecXQoUMJCAjgm2++4fHjxyxYsAAXFxe++OILxfvIch+Cra2t4r6BQ4cOER4ejrW1tf7SU/HixXn48GFWw+s9fvwYR0dH/WMHBwceP36c6loNvr6++juuZ8yYYfC61KT1vLRnY/uzSPPsJ63tTd2e0/Ix9XGZut0M+ZvLgQMH+Oeff3BwcMDS0pLOnTtTv359OnbsqLgoKCoIr94SHRsby549e2jXrp2iIIULFyYiIsKg7+Du3buZ6ktIS1oLVKemVatWtGrVSv84IiIi3X2n9by0Z2d7yvcye/aT1vambs9p+Zj6uEzdnn35m3tuNq1WS+HChQEoWLAgUVFRlCpVKlPLFSi6ZBQaGmrw78WLF3z55ZesXr1aUZBBgwbxwQcfcOTIEbRaLSdPnuTjjz/G29tbcaJpcXBwMHhjIiMjFa/kJoQQOdXRo0fx8PDA29ubo0ePZrh9rVq19Msae3h4MGzYMD799FMqVqyoOKaiMwSlfQVpGTt2LPny5WPYsGEkJCTg5eXFkCFDGDFixGvtF6B+/frs27ePd955h6CgIOzs7KQgCCFyJC8vL3bt2kWJEiW4cuWKvn3fvn2MGDECjUbDoEGD8PHxQaVSUbBgQV68eIFarc5w30uXLtVfMZk/fz7jxo0jKiqK3377TXF+6RaEu3fvYmFhoU8mNjaW6dOnc+XKFdzd3Rk9ejSWlpbpBrh27RqBgYG0b9+ekSNHKk4s2dy5cwkMDOTZs2d4e3vTo0cPEhMTAWjdujV16tQhICCAzz//HBsbG4YOHZrpGEIIkR369+/P8OHD6devn75No9EwbNgwDh48iFqtpkGDBnTq1AkPDw+aNm3Kv//+y5dffsnvv/+e7r5fHm5fvHhx/WzVmZFuQRg4cCBDhw7VF4Rhw4bh7+/P+++/z8qVK3n69Gm6QzxXrVrF4MGDKVq0KP/99x9r167N9AR3GRURlUqlvxNaCCHMycfHR//zq32WkHRz2KsLgyUvc5n8gd6zZ0927NhB1apVAShatChxcXGK4h87dowLFy4QHR1t0D5+/HhFr0+3IFy6dInWrVsDEBMTw8aNGzl27Bj16tVj4MCBeHp6plsQZs6cyZYtW+jcuTPbt29n2rRpMuOpECLPKleuXKZvTAsLC8PZ2Vn/WK1Wc/r0abZt28b+/fuJiorKcFI6gM8++4xNmzbh4eFB/vz59e1pDbJJTboFIT4+Xj8Z0tmzZylUqJB+ObbKlStnOEonPDyczp07A9ClSxcGDx6sODEhhHgTpDVS8v333+f9999XvJ/ff/+dK1euvNZop3RHGZUrV07fu/3nn3/SvHlz/XOPHj3Czs4u3Z2/fKAqlQqtVpvlRIUQIqcbMmRIpqetUKvVhIaG6h/fu3cvSx/qzs7OWV7ELFm6ZwiTJ0+mS5culC9fnmvXrhkMfdqxYwcNGzZMd+cxMTG4uLjoH//3338Gj0H5LdVCCJHTZWUuowYNGhAUFERISAhOTk5s2LCBdevWZTr28uXLGTx4ML169UqxVKbS+Y3SLQidO3fm/PnzXLx4kbp161KuXDn9c1WqVMHd3T3dnR8+fFhREkIIkRdkNJdRr169OHr0KBEREajVaqZMmcLAgQNZuHAhbdq0QaPR4OXlRbVq1TId+/z58+zduxd/f/8UfQhKv3hneB+Cq6urfoGFl73zzjsZ7rxp06aKkhBCiLwgozOEV1c0S+bp6Ymnp+drxR4/fjw7d+5MMbIpM0y2XNn8+fMzHCoVFxfH/PnzTZWCEEJkq6z0IRhLgQIFXmvqazBhQXjw4AFubm4MGTKEdevWcf78eW7cuMH58+dZv349Q4YMoUKFCkaZ4E4IIXKCJUuWcO7cObPEnjp1KiNHjuTBgwdotVqDf0qZbMW07777ji+//JJVq1axfPly/v77b6KioihatCg1a9bU38Pg4OBgqhSEECJbmXM9BC8vLyCpKCXL7FLFWSoIR44cwdLSMsPTE0dHR0aNGsWoUaOyEkYIIXIVc66YZoylihUVhKZNm/Ldd9/xzjvvMHPmTH766SesrKwYNmyY4luihRAirzPnGULyUsWvQ1FBuHLlCm+//TaQNKPe0aNHKViwIO+8844UBCGEyCH+/PPPVJfQVDrjqaJOZa1Wi0ql4ubNm+h0OqpUqYKzs7MsVSlyPQu0tCeQCRykPYFYIHfT51avvpcovG6eV0yZMoUhQ4ag1WrZvHkzDg4O7N+/nyJFiijeh6IzhHfffZfhw4dz//59unbtCsDNmzdz3BJyQmSKRsN+ltKIu9gRTyw2nMYFNL0w4XgLYQIWaFO8lzYfBcGS9ZDBFP3GZM4+hBUrVnDw4EGqV6/OypUrmTNnDr169eLbb79VvA9FZwirVq2iSJEi1KxZk8mTJwNJ6xwYY4EbIczF9vBhGnGXQsRjCRQinkbcxVbusM912nEtxXtpHRBASf/sfS/NeR9CVFQU1atXB8DGxoaEhAQaNmyoX0VNCUVfgxwcHFJMc92+fftMpCpEzmN95Qp2xBu02RGP7upVeLu1mbISWVGHsBTvper5cwpfu8qD5u+ZKavs5erqytWrV6lWrRrVq1dn8eLFFC1aNFMrSCoqCN98802az02dOlVxMCFykoTq1YnFhkIvfZDEYkNCFuaREeZ1AacU76Uuf37+q/zmvJfffvstkZGRAHz//ff07t2b6OhoFi1apHgfigrCy1OzQtJdyH5+fvr+BCFyo7gWLbiAS4o+hKotWkCsubMTmbGXypx+5b20qVuXB01aZGse5uxDeHkupEaNGhEcHJzpfSgqCCtXrkzRtm/fvjQnahIiV7C0pA2Dacc1ahPORUqzl8qEZmMnpEibBVracY06hHEBp3RHDWmxSPFe/rpuCsRm73tpzvsQAgMDOXbsGI8fP6ZYsWJ4eHjol+FUKstDKVq3bs2HH36Y1ZcLkSNosWA3VdlN5v7jCNNKa9SQBR3RpjEWJsV7+YYUdp1Ox8CBA1m9ejVqtZrSpUsTFhZGeHg4ffv2ZcWKFYqX0VRUEG7dumXwODY2lnXr1hmsAyqEEMby8qghSBo1pA0IoB0VpHi/4tdff+Xo0aOcOnWKBg0a6NvPnj1Lr169WLJkCd7e3or2pWjYqZubGxUqVMDNzQ03Nzfefvttjh07xurVq7N2BEIIkY60Rg3VJtxMGeVca9asYf78+QbFAJJWYps7dy5r1qxRvC/FdyprNBr9VKrR0dEcO3aMevXqZS5zIYRQIHnU0Mt0+fNzkdJv1B3JMTEx1KtXj127dqW5TWBgYJqLkTVt2pTAwEDF8RT3ISQmJnLixAnCwsJQq9W4u7tjZZU9d3NevHiRlStXotVqadmyJV26dDF4/urVq8yaNYsSJUoAST3s3bp1y5bchBDGl9aoof1/Vcx030JO4uXlxa5duyhRogRXrlzRt+/bt48RI0ag0WgYNGgQPj4+AMycOZMePXqku0+NRkOhQoVSfa5QoULGXw/h2rVrdOzYkefPn+Ps7ExoaCj58uVj586dVKlSRXGwrNBqtSxfvpwJEybg4ODAuHHjqF+/Pmq12mC7KlWq6H+JQojcLa1RQ21cJuXqvoX+/fszfPhw+vXrp2/TaDQMGzaMgwcPolaradCgAZ06dSI8PJyqVavy4sWLdPeZkJDAkSNHDCaze1liYqLi/BQVhKFDh/LJJ58watQofW/1jz/+yNChQzly5IjiYFkRHBxMyZIleeuttwBo3LgxZ8+eTVEQhBB5S2qjhtLrW8gJBeHlL6WtWrVKsb5xkyZNuH37tkHbmTNncHNzo3z58gD07NmTHTt2EB0dTUxMDIGBgeTPnx9PT08sLFKeBZUoUUK/OE5qkq+cKKGoIFy8eJGDBw8aDF0aOXIk06dPVxwoqx4/fmywqpqDgwNBQUEptrtx4wajR4+maNGi9O3bN9URUL6+vvj6+gIwY8aMDCfnS+t5ac/G9meR5tlPWtubuj2n5WPq48pke2p3JGNnx8WY0jki/3LlymX6xrSwsDCDzyu1Ws3p06dZuHAhkDSXnKOjY6rFAEhRYF6HooJQunRp/Pz8aNHi/+/6O3bsGKVLp3wTjC2106BXx9SWK1eORYsWkS9fPgICAvjhhx+YP39+ite9WrEjIiLSjZ3W89Kene0px09nz37S2t7U7TktH1MfV+baU+1bqFOHvX9VNnv+pUuXztKNaRl9xvXv3z/T+8wqRQXhu+++o1OnTnTo0IEyZcpw584ddu/ezdq1a02dHw4ODvr5OQAiIyNTTNZkZ2en/7lu3bosX76cp0+fYm9vb/L8hBDZJ62+Ba3LRnOnBmRt6gq1Wm0wPdC9e/ey5ct2ahR1y3fq1ImAgACqV6/Os2fPqF69OufPn6dz586mzg9XV1fu37/Pw4cP9SOdXv1lR0VF6atscHAwWq02zV53IUTulty3MJ1WSf0GufyO5AYNGhAUFERISAjx8fFs2LCBTp06mSUXxeNGK1asyIQJE0yZS6osLS3x8vJi+vTpaLVamjdvjrOzMwcOHACSptA4deoUBw4cwNLSEhsbG0aOHKn4Vm0hhDCWjC4Z9erVi6NHjxIREYFarWbKlCkMHDiQhQsX0qZNGzQaDV5eXlQz04y7aRaEvn37KvpQVbpW5+uoW7cudevWNWhr3fr/56tv27Ytbdu2NXkeQgiRnowuGaU1Iainp6fBbKXmkmZBcHNzy848hBAi1zPnbKfGkGZBmDRpUnbmIYQQuZ4510MwBsV9CPHx8Vy/fp2IiAiDYVIvD0UVQog3WZ49Q3jZX3/9Rffu3YmLi9MP53z27BnOzs4ppsYWQgiROykadvrFF18wZswYHj9+TKFChXj8+DETJ05k6NChps5PCCFyjSVLlnDu3Dlzp5Flis4Qbty4wYgRIwzafHx8KFeuHKNGjTJJYkIIkdvk9ktGis4QChcuzNOnTwEoVaoUgYGBPHnyhOjoaJMmJ4QQIvsoKgjvv/8+e/bsAWDgwIE0b96cevXq0b17d5MmJ4QQuckbcclo7ty5+p+/+uorGjZsSHR0NG3atDFVXkIIkevk9ktGWVryzMPDw9h5CCGEMDNFBeHu3btMmTKFCxcupOg3uHHjhkkSE0IIkb0UFYTu3btTuXJlpk6dSv78+U2dkxBCCDNQvKbyyZMn01yxRwghRO6fukLRJ3zHjh3x8/MzdS5CCJGrDRkyJNcWA1B4hjB//nwaN26Mq6urfrH7ZCtWrDBJYkIIIbKXooIwYMAALC0tqVKlivQhCCFEHqWoIBw+fJjw8HBZllIIIbLJP//8w7x584iIiKBly5Z8+umnJo+pqA+hZs2aBgvdCyFEZlmgpT2BTOAg7QkEjSbd9rzIy8uLEiVKUL16dYP2ffv2UalSJdzc3JgxYwYAVapU4ZdffmHTpk3ZdvezojOEFi1a0Lp1awYMGJCiD8HLy8skiQkh8g4LtOxnKY24ix3xxGKDzUdBWNGevSxP0W5BR7TKvq/mKv3792f48OH069dP36bRaBg2bBgHDx5ErVbToEEDOnXqRNWqVfnzzz+ZMWMGw4cPz5b8FK+H4OTkpF/YPplKpZKCIITIUDuu0Yi7FCIegELEow0IYAL5Um1vRwV2U9WcKZtEkyZNuH37tkHbmTNncHNzo3z58gD07NmTHTt2ULVqVTp16kSnTp1o3749H330kcnzU1QQjhw5Yuo8hBB5WB3CsPvfh34y1fPnvMPtVNtrE54rC4KPj4/+51atWtGqVasMXxMWFoazs7P+sVqt5vTp0xw9epRt27YRFxeHp6enSfJ9leK5jJ48ecLOnTsJCwvDycmJjh07UrRoUVPmJoTIIy7gRCw2+jMBAF3+/ByPLWtwhpDcfjG2tDnSfG3lypXL9I1pLy9JnEylUtGsWTOaNWtmxOwypugi3cmTJ3F1deWXX37h8uXLLFmyBFdXV06ePGnq/IQQecBeKnMaF55hgwZ4hg0JdevyLS1Tbd9LZXOnnG3UajWhoaH6x/fu3aN0afMUREVnCCNHjmTRokX07NlT37Zx40Y+//xzzp49a7Lkkl28eJGVK1ei1Wpp2bIlXbp0MXhep9OxcuVKLly4gK2tLUOHDtVfjxNCmJ8WC9owmHZcozbhXKQ0v66bQqLLxlTbtS4bzZ1ylmRl+usGDRoQFBRESEgITk5ObNiwgXXr1pkgu4wpOkO4ceMGPXr0MGjr1q0bwcHBJknqZVqtluXLlzN+/HjmzJnD8ePHuXfvnsE2Fy5c4MGDB8yfP59PPvmEZcuWmTwvIUTmaLFgN1WZTquk/gFLy3Tbc6OMFsjp1asX7u7uXL9+HbVazfLly7GysmLhwoW0adOGKlWq0KNHD6pVq5aNWf8/RWcIFSpUYMOGDQa93Js3b8bV1dVkiSULDg6mZMmS+uGujRs35uzZs6jVav02586do0mTJqhUKipWrEhMTAxPnjyRPg4hRI6yfv36VNs9PT2zreM4PYpXTOvQoQPz58+nTJky3L59m6CgIHbt2mXq/Hj8+DEODg76xw4ODgQFBaXYxtHR0WCbx48fS0EQQmSr3L5imkqXWhd3Kp48ecLu3bsJDw+ndOnSeHp6UqxYMVPnx8mTJ7l06RLe3t4A+Pv7ExwcbHD/w/fff0/Xrl2pXDmpI2rq1Kn06dMnRT+Cr68vvr6+APq7AYUQQiRRfCtg0aJF6dOnD2PGjKFPnz7ZUgwg6dv+y9NmREZGpvjm7+DgQERERLrbQNK44BkzZjBjxgyD8cLZSeLm7bjmjC1x83bc7KCoIISEhPDRRx9RtWpVXFxcDP6ZmqurK/fv3+fhw4ckJiZy4sSJFGN869evj7+/Pzqdjhs3bmBnZyeXi4QQIpMU9SF89NFHuLq6Mnv2bOzs7EydkwFLS0u8vLyYPn06Wq2W5s2b4+zsrJ9Go3Xr1tSpU4eAgAA+//xzbGxsGDp0aLbmKIQQeYGignD16lWOHz9utiU069atS926dQ3aWrdurf9ZpVIxaNCgTO1TyS3lpiBx83Zcc8aWuHk7bnZQ1KncoUMHpkyZQr169bIjJyGEEGag6AyhbNmytGnThvfff5+SJUsaPDd16lSTJCaEECJ7KSoIMTExdOzYkYSEBIM5N4QQQuQdiu9DSI1WqzVbv4IQ5qLT6VCpVOZOI9vI8b45svRp/vfffzN69GiD6SNyqteodxJX4hp4+PAhT58+JTY2NttiAty+fZt79+4RHh6erXHleN88itdDePToEevWrWP16tVcunSJd999l3nz5pkytywJDAwkODgYR0dHateunW3DZIODg3n06BElS5bEyckJGxubbIkrx5s9xxsQEMDWrVspVqwYzs7OtG3blkKFCpn8m+T58+dZv349ZcuWxcbGhm7dulGkSBGTn5nL8WbP8eY06RaEhIQE/vzzT1atWsX+/ftxc3OjV69e3Llzh82bN1OiRInsylORc+fOsXHjRmrWrMnt27cpWLAgNWvWBEx7GnjhwgXWrl2Lm5sbT58+pUuXLlSqVMkksV4mx5s9x3vp0iU2bNjAJ598gpWVFZs2bSIxMdHkHxaPHj1i/fr1DBs2jKJFi/L7778DEB8fT758+UwWV443e443J0q37L711lsMGTKESpUqcerUKQIDA5k4cWK2fRvMjPj4ePz8/BgyZAh9+/alePHi3Lt3j7CwMKKjo1GpVGi1WqPHvXHjBmvXrsXb25tPP/2UYsWK8e+//xIfH098fNIqUKa4rBEXF4e/v78cr4mPF5KWOOzVqxdubm4UK1aMsLAw1q5dy65du/j7779NEhPg+fPn2NvbU65cOWxsbLh27RorVqxg2bJl/PXXXyaLGx4ebpbjjYuLM8vxmuv9zYnSPUOoWbMmf/31F6dPn6ZChQqUK1cuR04JkfztMDo6msDAQBwcHDhx4gQuLi6EhIQQHx+Pl5cXhQsXNnpsR0dHBg8eTIUKFYiKiiIgIICIiAj++ecf7O3t6dy5s0kua6hUKp49e5btx1u8eHEGDRr0xhwvoJ+WOC4ujl9++YXmzZtTvXp1goKCOH36NK6uruTPn9/o3yhdXFzQarVMnDiRBw8e0KFDB9555x1u3bqFv78/lSpVonjx4kaNCdCuXTsg6UtWdh6vWq1Gp9Nl+/Emv7/Zfbw5ki4Dt2/f1k2dOlXn5uams7Gx0XXo0EFXqFAh3b179zJ6aba4dOmS7vDhwzqdTqe7fv26zsfHR/ftt9/q1q1bp9PpdLp///1Xt2TJEt25c+eMGvfJkye6qKgogzZfX1/d7t27dTqdThcYGKibO3euLiQkxOhx//vvP51Op9MFBQVl2/EmJCTo4uLiDNoOHjxo8uNNSEjQJSQk6HS6/z/eqVOnmvx4L126pNu7d6/++LRarU6n0+kiIyP129y9e1c3c+ZMXWxsrNHj7tq1S992//593ZIlS/Q5PH36VDd79mzdw4cPjRY3PDxc9+TJE51O9//HqtPp9G06nWmO9+W4yUJDQ3VLly416fFeuXLF4O83mamPN6fLsKemTJkyTJw4kaCgIA4dOkSpUqWwsLCgVq1ajBkzJjtqVpouX77M3LlzWb58OWFhYVSsWJEpU6ZQo0YN/YI6JUqUQKPR8OTJE6PFPXv2LHPnzmXWrFkcP35c396yZUv9t40qVaqg0WiIiooyetwZM2bg7++Pm5sbkydPpkaNGvr+HFMd76JFi1iwYAGXL1/Wj8Jo1aqVyY930aJFzJs3j4CAANzc3Pj222+pXbu2SY/32rVrzJs3DxsbG06ePMmKFSu4du0aiYmJFCtWTH9JLDw8nOfPn6PRaIwe99SpUyxbtozr169TsmRJbGxs9CsBXrt2jcjISCyNtLLY2bNnmTp1Klu2bCEiIsLg8tvLZ13GPt5X4yZTq9WoVCqWL18OGP94AwICWLZsGRqNhj179hissGjK480NMtV1/+677/Lrr7/y4MEDFixYYNbrawEBAaxZswYfHx969OjBpUuX0Gq12NjYUL16dY4fP46vry8nTpzg9u3bVK9e3Shxg4OD2bRpE15eXnzwwQfs27eP58+fp9ju9OnTPHjwwGhDc1+O261bNw4ePEhsbCy2trbUqFGDkydPmuR4b968yaZNm+jQoQPNmjVj69at7Nq1i4cPHxpsZ+zjfTXujh072LhxIy9evKBGjRqcOHHCJMcLSb/rtm3b0qJFCyZOnIidnR2nT5/m1q1bQNIlyt27d7Nt2zYGDBhAwYIFTRK3YMGCHD9+nJs3b9K6dWsePnzIpEmT2LRpE0OGDDHKFPSxsbEcOnQId3d3ihcvzp49e4iIiMDCwkJ/KVan07Fr1y6jHm9acZN5enry6NEjox9vREQE27Ztw8vLizZt2jBr1izu3r2rH+KafFnI2MebWygedvqyfPny0atXL3r16mXsfBSJiori5MmT9O/fn4oVK3L37l1OnjxJ27ZtgaRvGO3atWPfvn0UKFCAoUOHpphyI6siIiIoU6YMLi4uODo6otVqWbVqFZUrV6ZatWo4Ojpy+PBhdu/ezZdffmmwkpux465evZpKlSpRsWJFPD092bNnj9GP99GjR5QtW1a/2FBkZCS+vr44ODjQvHlzNBoNfn5+Rj/e9OJ6eHjQpk0bDhw4YPTjhaRr9zt37tQvBvXBBx+wefNmjh07RsWKFdFqtURHRzNs2DCjTgGfVlx/f38GDBjA119/zb1797C3t8fe3t4oMe3s7PDy8qJgwYKEh4dz9uxZdu/eTbt27fRnYXFxccTGxhr1eDOKW6pUKcaNG0dYWBiFChUy2vFaWVnx/vvvU6NGDRITE7GxscHa2pro6Gj9NnFxccTExBj9/c0NXutOZXPRaDTExcUZdF5Onz6dMmXK0KdPH31bYmIiOp0Oa2tro8UODQ1lw4YNFChQgMDAQJo3b46LiwtnzpzB1dWVtm3bEhISQr58+ShVqpRJ4zo7O3P27Fl93Pj4eFQqlVGP9+bNm+zfv593332XmjVrsm3bNv7991/Cw8Pp06cPlSpV4vbt29ja2hr1eFOL++DBA+7fv0/v3r2pXLkyCQkJAEY53oiICAoXLoxGo8HKyoo1a9ZQpkwZ6tSpQ9GiRYmPj2fSpEm0adOGZs2avXa8zMZ97733aNGihUnivjykMygoiLNnz5KQkEDfvn0JCQmhZMmSFChQIFvj3rp1i5IlSxrt23lERARFihQBkorCy3799VeaN29OhQoVCA4Oxs3NzSgxc6MsnSGYS3BwMAkJCdja2uq/OcbHx2NjY0Pnzp05ceIE0dHRFCxYEK1Wm+KNf524iYmJ5MuXj7Jly/Lhhx9y+/ZtoqOj+eCDDwDInz8/W7dupUmTJpQrVy7b4trZ2enjGmt0z8txXV1dUavV+Pr6snv3bgDGjRuHr68vly5dolKlSpQtWzbb4h48eJDLly9TuXJloxW+gIAAfv/9dypWrMjz58/p06cPTZo0Yd++feh0OipXroyTkxP16tUz6kiTzMQ11vXzV+PGxMTQs2dPSpcuDUCFChUA+Oeff/jmm28IDQ1l9uzZRikImYl77949fvzxR6MUhOS4FSpU4Pnz5/To0QMnJycSExOxsrIiNjaWuLg4/vrrLzZu3MjUqVNz5GjK7JBrCkJAQADr16+nTp06PH78GBsbGz755BP9PRGlS5fm9u3bnDhxgtatWxvtzsaX40ZERGBnZ8egQYNwcXEhODiYv//+mxo1avDixQtsbGzyXNwiRYrQr18/nj9/zr///qvvI4iOjjbq3aNK48bExBgtrk6nIzIykt9//x0vLy+cnJzw8/Nj4sSJTJs2jbZt23L8+HH8/f0pW7YsJ06cYNKkSXkq7rFjx5gyZQoTJkzA2dkZSPpwPn78OE+ePGHatGmvfRkwK3GnTp1qsrhTp041iFusWDG2b99OYmIiY8aMeWOLAZDxsNOc4MWLF7rp06frLl++rNPpdLpHjx7pBg4cqPv5558Ntjt58qTu22+/1b148cJg6Jyp4v7xxx+6OXPm6GbOnKkbM2aM0YZc5qS4AwYM0C1atMhgu927d+u+/PJLXWhoaK6Oq9PpdBqNRvfLL7/oIiMj9X8zO3fu1Hl7e+siIiJ0Ol3SkFpfX1/d/fv382Tc3bt36z755BNdWFiYTqfT6Z49e6YbNWqU7tatW29E3B07dug+/fTTHDOU3pxyxVSlOp2OfPny6YeEOTo60rRpU4KDg/ntt9/029WsWZMRI0Zga2trlFP7tOIGBQWxadMmOnbsSMeOHWnWrBlfffWV0S6d5KS4zZs3JygoSP97jo+PJzo6ms8++8xoI4rMEffBgwcEBwcTExNDbGwsx44d0//NdOjQgTZt2rBx40bi4+OpUqUKLVu2NErHdU6M6+npiaenJ9u3bycuLo6CBQsyffp0o1z6zMlxt23bRmJiIlWrVmXq1Kk4OTm9dtzcLkdfMkoebZEvXz7UajXfffcdH3/8MX///TeWlpaMHTuW7du3ExMTQ4ECBYx2DT2juD4+PmzdupX4+HhcXV1xdXXN03Ff/T13797dKAXXXHGTJ1ArUKAALi4ueHh4sHLlSrRaLV27dgWgcePGbN++3agd9Dk5rru7O3/88Yf+Eqwx4uf0uPfv38fKyuqN7kR+VY4tCOfPn2fOnDnUq1ePL774gh49euiHqNnZ2dGzZ0+srKx4+vQp8fHxRhsFoTRudHQ0z58/N9qkWzk97su/Z2N9KJsj7vXr11mzZg0jRoygXLlyLFmyhODgYKZNm8bXX3+NVqvlnXfe4dq1a4SEhBATE2OUjs3cEPfWrVv6uK/7u84NcZMHaLxJ9xlkJEcOO33x4gWzZ8+mUaNGXL9+nYSEBEaOHJliO39/fw4ePMjo0aONMk5Z4ubtuJD0gXH//n390NGnT5/y888/M27cOP7991+2bt2KjY0NwcHBDB061Gjj0CVu3o6bZ5ix/yJdkZGRuufPn+v+++8/3Y8//qibO3eu/rnExETdhQsXdD4+PkafO0fi5u24Go1GFxMTo/85IiJCN3r0aN3jx491Op1O9/DhQ11iYqJ+G4krcd8kObZTuVixYuTLlw97e3s++eQTEhISmD9/PpA0XW1iYiJjx441WoeqxH0z4lpYWOj7mnQ6HQUKFKBgwYIULVoUf39/tm/fjkajMfqMrRI3b8fNK3LkJaPUPH36lLVr13L9+nV0Oh2TJ082ytwmEvfNjgvw888/U7RoUS5dupSt0xVI3LwdNzfKsZ3Kr7K3t8fFxYULFy4wceLEbPuwkLh5N65Op0Oj0ehnMv3mm2+MOv2GxH0z4+ZmuaYgREdHc+HCBSZMmJCtFV7i5t24KpUKKysrPvjgA1xdXbPtw0Li5u24uVmuuWQE/z9vkcSVuMakM+F6zBL3zY2bG+WqgiCEEMJ0cuwoIyGEENlLCoIQQghACoIQQoj/kYIghBACkIIg3jDz589n0aJFBm2BgYF4eXnx5MkTM2UlRM4gBUG8UQYMGMCFCxe4fPkykDTUdcmSJfTr188oK2VpNJrX3ocQ5iLDTsUb5+TJk6xdu5bZs2ezbds2bt++Tbdu3fjtt9+4d+8exYsXp3///lSrVg2AI0eO8OeffxIZGYm9vT2dO3fmvffeA+Dq1assWLCAtm3bsnv3bmrWrMlnn31mzsMTIstyzZ3KQhiLu7s7J06cYN68eVy/fp2ZM2cyduxYhg8fTu3atbly5QqzZ89m7ty52NvbU7hwYcaOHctbb73FP//8w3fffYerqyvly5cHICoqiujoaBYtWoR8vxK5mVwyEm+kgQMHcuXKFbp168bx48epU6cOdevWxcLCgpo1a+Lq6kpAQAAAdevWpWTJkqhUKqpWrUrNmjW5du2afl8qlYoePXpgbW1tljuthTAWOUMQb6QiRYpgb2+PWq3mzJkznDp1ivPnz+uf12g0+ktGFy5cYMuWLYSHh6PT6YiLizOYb8ne3l4KgcgTpCCIN56DgwMeHh54e3uneC4hIYHZs2czfPhw6tevj5WVFbNmzTLYRubJEXmFXDISbzwPDw/Onz/PxYsX0Wq1xMfHc/XqVSIjI0lMTCQhIQF7e3ssLS0NRigJkdfIGYJ44zk6OjJmzBjWrl3LvHnzsLCwwM3NjcGDB5M/f34GDBjAnDlzSEhIoF69etSvX9/cKQthEjLsVAghBCCXjIQQQvyPFAQhhBCAFAQhhBD/IwVBCCEEIAVBCCHE/0hBEEIIAUhBEEII8T9SEIQQQgBSEIQQQvzP/wGdqbv0kdQZHgAAAABJRU5ErkJggg==\n", "text/plain": [ - "
" + "
" ] }, "metadata": {}, @@ -108,6 +111,59 @@ "plt.close()" ] }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Percentage of maximum water area estimates above the median 1.80% when damage is above zero: 90.0%\n" + ] + } + ], + "source": [ + "# Calculate the median of the maximum water area from the df_annualmax dataframe\n", + "median_max_water_area = df_annualmax['ensemble_water_area_km2_anomaly_pct'].median()\n", + "\n", + "# Merge the annual max water area dataframe with the damage dataframe on the 'water_year' column\n", + "merged_df = df_annualmax.merge(df_sac_damage, left_on='water_year', right_on='year', how='inner')\n", + "\n", + "# Now, filter for entries where damage is above zero\n", + "filtered_df = merged_df[merged_df['damage'] > 0]\n", + "\n", + "# From the filtered dataframe, calculate the percentage of records where the max water area is above the median\n", + "# This is done by first checking if the water area is above the median for each row, then calculating the mean of this boolean series\n", + "percentage_above_median = (filtered_df['ensemble_water_area_km2_anomaly_pct'] > median_max_water_area).mean() * 100\n", + "\n", + "print(f\"Percentage of maximum water area estimates above the median {median_max_water_area:.2f}% when damage is above zero: {percentage_above_median}%\")" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Number of years where the maximum water area is above the median and corresponds to a year with flood damages: 9\n" + ] + } + ], + "source": [ + "# From the filtered dataframe, find the instances where the max water area is above the median\n", + "years_above_median_with_damage = filtered_df[filtered_df['ensemble_water_area_km2_anomaly_pct'] > median_max_water_area]\n", + "\n", + "# Count how many years satisfy the criteria\n", + "number_of_years = years_above_median_with_damage.shape[0]\n", + "\n", + "print(f\"Number of years where the maximum water area is above the median and corresponds to a year with flood damages: {number_of_years}\")\n" + ] + }, { "cell_type": "code", "execution_count": null, @@ -118,7 +174,7 @@ ], "metadata": { "kernelspec": { - "display_name": "Python 3", + "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, @@ -132,7 +188,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.8.18" + "version": "3.9.7" } }, "nbformat": 4, diff --git a/table1.ipynb b/table1.ipynb index eab2725..bdb9f7e 100644 --- a/table1.ipynb +++ b/table1.ipynb @@ -13,7 +13,12 @@ "from rasterio.warp import Resampling\n", "import matplotlib.pyplot as plt\n", "import matplotlib.lines as mlines\n", - "LOC_DATA_DIR = \"/mnt/d/2023_nwm_cnn_california/\"\n", + "\n", + "#### User needs to modify this line to correspond to their local data, \n", + "#### which can be downloaded from Hydroshare:\n", + "#### http://www.hydroshare.org/resource/dbf8e4c2a39a4c228db867b04f9c21ed\n", + "LOC_DATA_DIR = \"../2023_nwm_cnn_california/\"\n", + "\n", "STATS_DIR = f\"{LOC_DATA_DIR}images_for_sacramento_stats/\"" ] }, @@ -912,7 +917,7 @@ ], "metadata": { "kernelspec": { - "display_name": "Python 3", + "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, @@ -926,7 +931,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.8.18" + "version": "3.9.7" } }, "nbformat": 4,