{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Lab 5-1: Calculating the Water Balance for the area draining to Cedar Lake, Washington\n", "---" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We are going to use observational data to test the conservation of mass in a watershed. Monthly climate data are from the Cedar Lake station at 1560 ft elevation, at a latitude of 47.25 N and longitude of 121.44 W: [monthly climate from 7/1/1898 to 6/10/2016](https://wrcc.dri.edu/cgi-bin/cliMAIN.pl?wa1233)). Streamflow is from USGS gauge 12116400 for years 2001-2019. The gauge is listed as having an area of 217 $km^2$. Monthly ET is from climatological [maps](http://www.cpc.ncep.noaa.gov/soilmst/e.shtml). Pan evaporation is from the [Western Regional Climate Center WRCC](https://wrcc.dri.edu/Climate/comp_table_show.php?stype=pan_evap_avg). The Maple Leaf Reservoir site in Seattle, at a lower elevation, operated 1941 to 1960, and the Lake Kachess site, at a higher elevation, operated 1931-1977. Note the different units and time periods. This is often what you get from observational data. Your job is to make sense of it. " ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "import scipy.stats as stats\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Load the csv file with these various measurements from the Cedar River Watershed or nearby vicinity. " ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
monthAvg_max_T_degFAvg_min_T_degFAvg_total_P_inAvg_streamflow_cfsAvg_ET_mm_per_monthPan_evap_MapleLeaf_inches_per_monthPan_evap_LakeKachess_inches_per_month
0140.030.113.5221830.610.00
1243.831.510.0614170.820.00
2347.633.010.36132171.800.00
3453.736.28.01130303.262.37
4560.841.46.52142504.644.82
\n", "
" ], "text/plain": [ " month Avg_max_T_degF Avg_min_T_degF Avg_total_P_in Avg_streamflow_cfs \\\n", "0 1 40.0 30.1 13.52 218 \n", "1 2 43.8 31.5 10.06 141 \n", "2 3 47.6 33.0 10.36 132 \n", "3 4 53.7 36.2 8.01 130 \n", "4 5 60.8 41.4 6.52 142 \n", "\n", " Avg_ET_mm_per_month Pan_evap_MapleLeaf_inches_per_month \\\n", "0 3 0.61 \n", "1 7 0.82 \n", "2 17 1.80 \n", "3 30 3.26 \n", "4 50 4.64 \n", "\n", " Pan_evap_LakeKachess_inches_per_month \n", "0 0.00 \n", "1 0.00 \n", "2 0.00 \n", "3 2.37 \n", "4 4.82 " ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Read in a .xlsx file\n", "data = pd.read_csv('../data/Cedar_average_monthly_waterbalance.csv')\n", "data.head()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "These data are month of year (1=January), monthly average values of maximum temperature ($^\\circ$F), minimum temperature ($^\\circ$F), total precipitation (inches), average streamflow (cubic feet per second, cfs), evapotranspiration, ET (mm per month), pan evaporation at Maple Leaf Reservoir in Seattle (inches per month), and pan evaporation at Lake Kachess near Snoqualmie Pass (inches per month)." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Plot the temperature data." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0sAAAI3CAYAAABK08B4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAACft0lEQVR4nOzdd1hT5xcH8G/YYBE3iOCquCe17r333ttaR7WuDketC1fV1qq1zlpH1brqqtY6WsU96qpbq7i3Iqgg8/39cX6JIEEJEm4Svp/nyQPcm4QDWffc933P0SmlFIiIiIiIiCgeO60DICIiIiIiskRMloiIiIiIiIxgskRERERERGQEkyUiIiIiIiIjmCwREREREREZwWSJiIiIiIjICCZLRERERERERjBZIiIiIiIiMoLJEhERERERkRFMlojiOHz4MHQ6HXQ6HSZNmqR1OESUim7evImhQ4fC398fmTJlgpOTEzw9PVGnTh3Mnj0bz58/N+vv79atG3Q6HXbv3m3W3/M2uXPnhk6nw7Vr18xy/9WqVTPr/ZtqwIAB0Ol0GD16tNH9y5cvN3wunDhxwuh1WrduDZ1Oh1mzZpkzVADJe54sXrwYOp0OY8aMMVtcRLaKyRJRHL/88ovR7yllXLt2DTqdDtWqVdM6FJtlKQfcqSklnlfz5s1D/vz5MWXKFNy+fRsVKlRAq1atUKRIEezbtw/9+vVDvnz58OjRo5QLnCxC5cqVAQB79+41un/fvn2G7xO7zv79++PdF8XHZI2sGZMlov+LiorCqlWroNPp4OXlhfPnz+P48eNah0VEZvbTTz+hT58+cHBwwJIlS3Dv3j1s3rwZK1aswN9//41Hjx5hypQpCAsLM/voEqU+fYJz+PBhREVFJdi/b98+5MuXD+nSpTOaLF25cgV3796Fh4cHihUrZvZ4k6N58+Y4f/48Pv30U61DIbI6TJaI/m/r1q149OgRqlSpgl69egHg6BKRrbt16xb69+8PnU6HTZs2oUuXLtDpdPGu4+bmhi+//BKHDx+Gh4eHRpGSuXh5eSFfvnwICwtLcILs6dOnOHfuHKpWrYoyZcrEG2XS02+rWLEi7Ows87DKw8MDBQsWRJYsWbQOhcjqWOarmkgD+sSoU6dO6NSpEwDg119/RUxMjOE6x44dg06nQ7ly5RK9nylTpkCn02HEiBHxtkdGRmLGjBn48MMP4e7ujnTp0qFMmTJYuHAhlFIJ7ken0yF37tyIjIxEQEAAChYsCGdnZzRr1gwA8PLlSyxcuBBNmzZF3rx54erqigwZMqBKlSpYuXJlovE9ePAAPXv2hKenJ9zc3ODv748VK1a8dSrT77//jrp16yJz5sxwcXFB/vz5MXLkyCSfaR8zZgzy5MkDAAgMDDSsAdDpdOjWrVu86z58+BBffPEFChQoABcXF2TMmBH169fHnj17Etzv7t27Dffx4MED9OjRA15eXnjvvfdQqVIlHDhwwHDduXPnonjx4nB1dYWvry/Gjh2L2NjYBPcZ938/evRovP/++3BxcUHevHkxatQovHz50ujfqPVjrNPpsGTJEgBA9erV4/2P9etD3jZNTx9TXHGn0Fy6dAnt2rWDp6cn7OzssGHDBsP1Tp8+jY4dOyJHjhxwdnaGt7c3unfvbvLalJMnT2LIkCH44IMPkDVrVjg7OyNv3rzo27cv7ty5E++6pjyvjJk1axZevnyJNm3aoHr16m+8bqFChZAxY8Z4254/f46AgAAUK1YMbm5uSJ8+PapWrRrv//K63377DWXKlIGrqys8PT3RpUuXBH9XXHv37sWnn36K4sWLI2PGjHB1dUXBggUxbNgwPH36NMH1474m7t27h48//hg+Pj5wcHDA9OnT3/g3murp06f44YcfULduXeTKlQvOzs7InDkz6tWrhx07dph0XyEhIahSpQp0Oh0GDRoU7zXzru8/b5PYVLz9+/cjNjYWFStWRMWKFXHv3j38999/8a6jT5b096GUwq+//op27dohf/78SJcuHdzd3VGmTBnMnj3b6HuOUgorV65ElSpV4OXlBRcXF/j6+qJWrVr48ccfE417z549qFGjBtzd3ZE+fXo0bNgQ586dS3C9xKbBxX0/SOp9AfK8/+KLL+Dr6wtXV1cULlwYM2fOhFIqwXtItWrV0L17dwDA2LFj471GFy9eHO9+//jjD9SuXRsZM2aEi4sLChQokOjzfMyYMYb7OH36NJo0aYKMGTMiXbp0qFq1arz3fqJ3oohIPX36VLm4uChnZ2cVHByslFKqTJkyCoDaunVrvOsWLFhQAVD//fef0fsqUaKEAqDOnDlj2Pb8+XNVuXJlBUBlyZJF1atXTzVo0EBlzJhRAVC9e/dOcD8AlK+vr6pfv75Kly6datCggWrdurXq06ePUkqp8+fPKwDK09NTVa1aVbVt21ZVrVpVOTo6KgBq9OjRCe7z4cOH6v333zfcd9u2bVW1atWUnZ2dGjBggAKgqlatmuB2n332mQKgXFxcVJUqVVSLFi1Urly5FAD1wQcfqOfPn7/1f7x+/XrVsmVLQ8xdu3Y1XBYsWGC43vnz51WOHDkUAPX++++r5s2bqypVqignJydlZ2enli9fHu9+d+3apQCoJk2aqLx58ypvb2/VrFkzVbZsWQVAubm5qTNnzqgBAwYoFxcXVbVqVdWoUSPl7u6uAKivvvrK6P8+Z86cqnHjxsrV1VU1atRItWjRQnl4eCgAqmbNmio6OjrebSzhMe7atavh8a1bt268//HDhw8N1wGgdu3aZfRxAqBy5coVb9uiRYsUANWuXTuVPn16lSdPHtW2bVtVp04dtXnzZqWUUmvXrlVOTk6G50SrVq1UqVKlFACVOXPmeK+Ht2nbtq2yt7dXJUqUUE2bNlXNmjVTuXPnVgBU9uzZ1e3btw3XTerzKjFFixZVANTGjRuTHJ/evXv3VOHChRUAlSNHDtWkSRNVq1YtlS5dOgVATZo0KcFtfvjhBwVA2dvbqxo1aqg2bdqo7NmzK19fX9WoUSOjj03ZsmWVs7Oz+uCDD1SLFi1Uw4YNVfbs2RUAVaRIEfXs2bN419e/Jho0aKB8fHyUl5eXatWqlWrUqJGaN2/eW/8u/Ws7KCjordfdunWr4Xlcs2ZN1bZtW1W+fHml0+mUTqdTCxcuTHCbqlWrJrj/e/fuqZIlSyoAauzYsfGunxLvP2/z888/G95H4ho2bJgCoC5dumT4W3/++ed419F/Juzbt08ppVR4eLgCoDJmzKgqVqyo2rZtq2rWrKnc3NwUANW1a9cEv3/o0KEKgHJ3d1f169dX7du3V9WqVVNZsmRJ8HrUv4Y/++wzw+ukZcuWKn/+/IbX2927d+PdRv8aNvaeYep9hYeHGz4fs2bNqlq1aqXq1aunnJycDJ8jcWOeNGmSqlixogKgSpQoEe81unfvXsP1Jk6cqAAoBwcHw3PJx8dHAVD58+dX9+7dixfH6NGjFQDVr18/5ebmpvLnz69atmxp+Ax2cXFRp0+fTvC/JjIVkyUipdT8+fMVANWyZUvDtpkzZyoAqmPHjvGuGxAQoACogICABPdz7tw5BUCVLFky3vZPPvlEAVCdO3eOd2Dz4MEDw0G9/qBTD4ACoPLly6du3bqV4Hc9evRIbdu2TcXExMTbfvXqVZU7d25lZ2eX4GCnR48eCoBq3ry5evnypWH7zp07DQe6rydLq1atUgBUqVKl4t1fZGSk6tWrlwKgvvjiiwTxGRMUFJRoQqaUUtHR0YaD1xkzZqjY2FjDvuPHj6vMmTOrdOnSqfv37xu26w8MAai2bduq8PBwwz79h2nhwoVVjhw54h2wnz17Vjk5OSk3N7cEB5v6+/Px8VFXrlwxbH/w4EG8+OKylMf4bcnQuyRLANSnn36aIFG8evWqcnNzUx4eHiowMDDeviVLligA6sMPPzT6+4z566+/1J07d+Jti4mJUWPHjlUAVPfu3ePte9vzKjERERFKp9MpAOrmzZsm3VYpperXr68AqCFDhqjIyEjD9itXrqj3339f2dvbq1OnTsWL09nZWTk7O8f7/7948ULVrl3b8D9+/bHZsmWLevLkSbxtL1++NLz+Xk8u4r4mmjdvHu81kRSmJEtXr15V+/fvT7D9+PHjKkOGDCp9+vQJXl+vJ0tBQUEqX758SqfTqR9++CHedVPy/edNLl++bEgO4r7vVKpUSWXLlk0pJSfV7Ozs1EcffWTY/+jRI6XT6ZSLi4vhPTUqKkr99ttvKiIiIt7vePDggSpdurQCEO91Eh4erpydnVXu3LnV48eP490mKioqwWtK/xq2s7NTK1asMGyPjo42nDgYOXJkvNu8LVky5b7GjRunAKjy5curkJAQw/ZTp04ZTg4l9h5i7CSeUkodOXJE2dnZKXd3d3X48GHD9pcvX6rWrVsrAKp169bxbqN/fwegJk+eHG/foEGDDO/HRO+KyRKRUoYRgfXr1xu2PXjwQDk4OCQ4mL5y5YoCoAoUKJDgfkaMGKEAqKlTpxq23b9/Xzk6Oqo8efLES1D0Tp48qQCoxo0bx9uu/xBYs2aNyX/PggULFAA1c+ZMw7Znz54pFxcX5eDgoG7cuJHgNl26dDF6wKk/S3fhwoUEtwkPD1deXl4qQ4YMCQ7ojXnbQe369esVANW+fXuj+6dPn64AqO+++86wTX9g6OHhYRgV1AsJCTEcDL9+NlgppZo3b2704FT/v58/f36C2+jPLufPn9+wzVIeY6XMmyxlzZpVvXjxIsFtBg4cqAAkOmrRrFkzBUAdO3bsrX/X2+TIkUNlypQp3rbkJkv37t0zPAbGHrc3OXHihAKgKlSoEO/gWm/Dhg0KgOrfv79h28iRIxUA1bNnzwTXv3DhguG5mthj87qwsDDl4OCg/P39423XvyacnZ2NJuFvY0qy9Cb698NNmzbF2x43WTpz5ozy9vZWDg4OatmyZQnuIyXff95GP1p39uxZpZQk0y4uLqp58+aG6xQrVkz5+fkZftY/zlWqVEnS79ixY4dhJEfv/v37CoBq2rRpku5D/xru1KlTgn3Hjh0z+lp4W7Jkyn3pR3sOHjyY4DajRo1KVrKk//x5PTFTSv4/rq6uys7OLt7zWZ8sVapUKcFtHj16ZDQOouRwAFEad+3aNezbtw+ZMmVCgwYNDNuzZs2KunXrYsuWLVi/fj06d+4MAMibNy/KlSuHQ4cO4fjx4/D39zfcZuXKlbCzs0O7du0M2wIDAxEVFYV69erB2dk5we8vUaIE3N3dcfTo0QT7dDodGjdu/Mb49+3bh927d+P27dt4+fIllFK4e/cuAODy5cuG6x0/fhwvX75EpUqV4Ovrm+B+WrdujaVLl8bb9uDBA5w6dQqFChVCgQIFEtzGxcUFpUuXxubNm3H58mWj1zGFfo2Dfs3O6ypVqgQARv9XpUuXRoYMGeJtS58+PTJnzoxHjx6hdu3aCW7z/vvvA4Dh//W6uI+jXr169ZAxY0ZcunQJDx8+RNasWS3mMTa3WrVqwc3NLcF2/ePWtGlTo7erVKkSNmzYgKNHj8Z7vbzJ48ePsWnTJpw5cwZPnz41rB2MiorCkydP8OTJE2TKlCmZf4lQRtaRJVXcv/n1ghCA8eeqfm1LmzZtEly/QIECKFWqVKIVOG/fvo3ff/8dFy5cQGhoqGHdi5OTU6LPAX9/f+TIkcOEvyp5YmJi8Ndff+HAgQO4d++eYU2fPq7E4jt06BD69euH8PBwbNiwAQ0bNoy3P7XffypVqoQ1a9Zg3759KFy4MI4ePWp4z9SrWLEi5s6di/v378PT0zPBeqW4Tp48ie3bt+P69esICwuDUgrPnj0DEP9/ki1bNvj4+GDLli2YOnUqOnbsCG9v77fGW6dOnQTb8ufPDyDx97R3va8bN27g1q1b8PHxMbp2t3Xr1ggICDDpdwOv1op17Ngxwb5s2bKhTp062LhxIw4cOIDWrVu/NfbMmTMjc+bMJv8fiIxhskRp3rJly6CUQps2beDk5BRvX6dOnbBlyxb88ssvhmQJkDf0Q4cOYfny5YaDv0OHDuHKlSuoXr06fHx8DNfVL26fM2cO5syZk2gc4eHhCbZly5bN6ME3IIuhW7Rogb///jvR+9R/MAMwLCA3ligBQM6cORNsu379OgDg/PnzRg8I43r06NE7H6zo/1dt27ZF27Zt3/i7XpfYQWG6dOnw6NEjo/vTpUsHAIiIiEiwL2PGjHB3dzd6n7ly5UJwcDDu3LmDrFmzWsxjbG7GniPAq8fNy8vrjbdPao+iX3/9Fb169Xrj4v1nz569c7KUKVMm6HQ6KKXw8OHDeK/bt9H/zUOHDsXQoUMTvV7cv1n/Gkzs/5gzZ06jydK0adMwfPhwREZGJjm+N/2elHTr1i00atQIp06dSvQ6iT1HO3fujOjoaKxatSpBogSk/vtP5cqVsWbNGuzduxe9evWKV+VOT58s7du3Dy1btjSaLEVGRqJbt2749ddfE/1dr/9PlixZgnbt2mHIkCEYMmQI8uTJgypVqqBDhw5GkwEARp+v7733HgDj72lvktT7Ss7nSFLcuXMHOp0OuXLlMrpfXzDCWCGUxF637733Hh4/fpyseIjiYrJEad6yZcsAAH/99Ve8M4jAqw+Jv/76C3fv3kX27NkByMH84MGDsXLlSkydOhV2dnZYsWIFgIRnxvRnxEuVKoXixYubFJuLi0ui+4YOHYq///4bVapUQUBAAIoWLYoMGTLA3t4e27dvR926dROtwGaMsevqY8+ePXuiH9h6mTNnfuP+pND/vvr16yNbtmyJXq9gwYIJtr3tYOpt+03x+v/K0h7j5DJWpSspscbExECn06FLly5vvH2RIkXeGsP169fRrVs3KKUwffp0NGzYEDly5ICrqysAoEKFCjh48GCK/N1OTk4oXLgwzp49i+PHj5uULOkf88qVKyNv3ryJXi9uqWZ9zKY8Fw8dOoTPP/8cHh4emD9/PqpVqwYvLy9Dgu3t7Z3o2fM3PbdSyscff4xTp06hRYsWGDp0KAoUKAB3d3fY2dlh/vz56N27d6KPVfv27fHLL79g1KhRhipwcaX2+8/rFfH27dsHV1fXeKOh+sRp7969aNiwIY4fPw57e3uUL1/ecJ1p06bh119/RdGiRTF16lT4+/sjY8aMcHR0xKVLl1CgQIEE/5MaNWrgv//+w+bNm/Hnn38iMDAQS5YswZIlS9CmTRusWrUqQbwp+Z5m6n2l5O9+19+rVSyUdjBZojTtyJEjuHjxIgCZFpHYdJHY2FisWLECn3/+OQCZole7dm1s3boVu3fvRtWqVbF69Wo4OzujZcuW8W6rPwCrVq0apk2blmKxr1+/Hvb29ti0aVOC3i9Xr15NcH19onfjxg2j93fz5s0E2/Sxe3l5JSjxag7639enTx80adLE7L/vTYKDg/Hs2TOjo0v6/6H+f2opj3FS6EdPjY3aGHsOJIWPjw+uXLmCmTNnIn369Mm6D70//vgDkZGR+PzzzzFw4MAE+5P7dyemYcOGOHv2LFasWGHSc07/mLdq1QoDBgxI0m28vb1x6dIlXL9+HX5+fgn2G3ttrl+/HgAwfvx4dO3aNd6+8PBw3Lt3L8kxp7QXL15gx44d8PT0xOrVq2Fvbx9v/9seq4CAAHh7e2Py5MmoWbMmdu3aFe8kSWq//xQvXhweHh64fv06bty4gQMHDqBMmTJwdHQ0XCdPnjzInj079u3bh8OHDyMyMhL+/v7xnvf6x0yfMMX1pv9J+vTp0aFDB3To0AGAJMqtW7fG6tWr0a1bN9SvXz8l/9xkedvnSGLb38bb2xtBQUG4fv260RFC/Sij/vcTpSb2WaI0Td9b6csvv4SSgicJLtu3bwfwagRKTz+CtGLFCvz111+4f/8+GjZsmGDdTPXq1WFvb4/NmzfH69n0roKDg+Hu7m60Sebq1asTbPP394ezszMOHTqEW7duJdi/du3aBNt8fHxQoEAB/PvvvwgKCnrnmPUH6tHR0Ub316pVCwDe2KMmNRk7m7tt2zYEBwfDz8/PcGBnKY8x8Pb/sf5g49KlSwn26Z/rpkrJxy04OBiA8Wk+e/bswf379xNsf9vf/CaffvopnJ2dsXr1auzateuN171w4YIhvuT8zfqR6zVr1iTYd+nSJZw8eTLB9jf9P9asWZOiI4umCgkJQWxsLLJnz54gUYqOjjYkDW/yzTff4Msvv8S5c+dQs2bNeNMWU/r9523s7OxQoUIFAMD8+fPx5MmTBLMNABldOnnyJP78808ACdcrvekxS+x1a0y5cuUM079Pnz6d5NuZU65cueDt7Y1bt27h8OHDCfYb+xwB3v4a1f8Ply9fnmDfw4cPsX379niPD1FqYrJEaZZ+rjwg00ESU6NGDWTLlg0nT57EmTNnDNubNWuGdOnS4bfffsOiRYsAGF+cmiNHDnTr1g2XL19G586dja7bOHDgAP744w+T4s+fPz+ePn2a4ID++++/N3rQ5+7ujg4dOiA6OhqDBw+Ot/5h9+7dic6v//rrrxETE4OWLVvG+/v1rly5gp9//jlJMWfJkgWOjo64cuWK0aSiVatWKFiwIBYvXozJkycjKioq3v7IyEisW7cu1Q4cAgIC4jVUffToEYYMGQIA6Nu3r2G7pTzGAAwLw/Ujpq+rWrUqAFlfFXc+//HjxzFy5EiT4tP7/PPP4erqisGDB+P3339PsP/JkyeYPXu20TVbr9MvKl+2bBlevHhh2H779m306dPH6G3e9rx6E19fX0yfPh1KKTRp0gS//PJLggQkPDwc06dPR9myZRESEgJADmT1oyGDBw9OMFIXGxuL7du3G9a0AED37t3h5OSEpUuXxmt+Gh4ejoEDBxqdBqn/fyxcuDDe6+HcuXNvXCuVGrJlywYPDw+cOXMG+/fvN2yPiYnBkCFDjCbkxkyZMgWfffYZzpw5g5o1a8Z7Xibn/UffYDs507P0B+36RrCJJUsxMTGYO3duvNvo6R8z/X69tWvXJiiiA8hozOLFixEWFhZve0REhOF1nhrrz5Kqd+/eAOR1H3ft1ZkzZ/DDDz8Yvc3b3pf69esHOzs7zJgxA//8849he2RkJPr374+wsDC0aNEiVQqWECWQeoX3iCzLpk2bEi0B/rq+ffsqAGro0KHxtnfo0MFQetjDwyPR8sMvXrxQ1atXV4A0HaxcubKhwai+AevAgQPj3QZvKXu6bNkyw++uXLmyat++vSpcuLCys7NTgwcPVkDC5of3799XefLkUYA0XW3Xrp2qUaOGsrOzU59++qkCoGrXrp3gdw0ZMkQB0kizdOnSqnXr1qpu3bqGZowlSpR46/9Qr3HjxgqQZpqdO3dWPXr0iFfW+/z58ypnzpwKkAakdevWVa1bt1blypVTGTJkSFDiXV8m2VijR6VelUE2Rl96dtGiRfG26/8/jRo1Um5ubqpx48aqRYsWht9fvXp1FRUVFe82lvIY//PPP0qn0ylnZ2fVtGlT1aNHD9WjRw/16NEjpZRSsbGxhtLN2bJlU82bN1eVKlVSjo6O6osvvkhW2V+llPrtt9+Uq6ur4TXVrFkz1bRpU1WyZElDD6/XS7sbExERoYoUKaIAKC8vL9WyZUvVsGFD5ebmpipUqKAqVKhgtKz1255Xb/Pjjz8qZ2dnBUhz20aNGqkOHTqoGjVqGJqJent7G5r7KiWlx4sXL64AqEyZMqkaNWqotm3bqkqVKqmsWbMqAOr777+P93u+//57w2tJ33jT29tb+fj4GG1K++jRI+Xl5aUAqDx58qg2bdqoWrVqKUdHR9W6dWujz++3vSbeRn+fJUuWVGXLljV60f9dEyZMMPw9tWvXVm3btlW5c+dWrq6uql+/fkafN8aa0ir1qjdOyZIl4/WVMvX9R9/ewcHBweS/fe/evYbXnJ2dnXr69GmC6xw5csRwHQAJmqUGBgYqe3t7BUjT3Pbt2xv6K+lfY3HLcevL0Lu5uakqVaqoDh06qKZNmxqeQ2XKlInXs+ldyv8nVjrclPsKCwsz/D1Zs2ZVrVu3VvXr11fOzs6Gz5G45dWVkjLv2bJlM/zt3bt3Vz169IjXo0v/XHJwcFC1atVS7dq1U76+vob7S6wp7evv33pveu8nMgWfRZRm6RvdvekAUE//Aerj4xOvn8eWLVsMH5hxGxUaExUVpX766SdVtWpVlTFjRuXk5KR8fHxUlSpV1JQpUxI0xXzbgbT+95crV065u7urDBkyqFq1aqndu3e/8WDp3r176qOPPlJZs2ZVLi4uqkSJEmrJkiVq3759CoBq166d0d/1119/qebNmysvLy/l6OiosmXLpvz9/dWXX35pUv+c+/fvq86dOysvLy/DAcXrcT558kSNGTNGlShRQqVLl065ubmp999/XzVp0kQtWrQoXt8rcyVLuXLlUi9fvlRfffWVyp07t3JyclK5cuVSI0aMUGFhYUbvz1Ie4+XLlyt/f39D8vL6genTp09Vnz59lKenp3J2dlZFihRRc+bMSTSmpCRLSil16dIl1bt3b5U3b17l7OysPDw8VKFChVT37t3V5s2bjfYjMubJkyfqk08+Ublz51bOzs4qb968aujQoerFixeJHmgn5Xn1NtevX1dffvmlKlGihPLw8FAODg4qW7Zsqk6dOmrOnDnq+fPnCW4TFhampk2bpsqWLavc3d0NzUXr1Kmjfvzxx3jJld7q1avVBx98oJydnVWWLFlUhw4d1K1btxI9cL1586bq0KGDypEjh3JxcVGFChVSkyZNUtHR0WZNlt50iZv4L1myRJUqVUq5ubmpzJkzq6ZNm6pTp04l+rxJ7DFUSqn+/fsbkoy4ybUp7z/r1q1L0nuyMS9fvlQuLi4KgCpevLjR60RFRRkS6NeTAr2DBw+qGjVqqIwZMyp3d3dVoUIF9dtvvxntCRYaGqq+/fZb1aBBA5U7d27l4uKismTJoj788EM1c+bMBO83WidLSkkPu8GDB6scOXIoJycnVaBAAfXdd9+pmzdvKgCqXLlyCW5z9OhRVbt2beXh4WHoKfb6e+/mzZtVzZo1lYeHh3JyclL58uVTQ4YMSdCUWSkmS5R6dEppOOGZiCzG5MmTMWzYMHzzzTeaT+/Rmr6EbdwpeERkHQYMGIC5c+fi0qVLhpLTlDpWrVqFdu3aoU+fPm9so0BkTbhmiSiNMdbHZc+ePZg4cSIcHByMNswkIrIWu3btQvfu3ZkomdHJkycTrLE7ffq0YU2nvqIfkS1g6XCiNKZChQrw9vZGoUKFkC5dOvz33384ceIEAKlMlSdPHo0jJCJKPkupHGfL2rVrh9DQUBQrVgwZM2bEtWvX8M8//yAmJgZ9+vRJUPSCyJoxWSJKY4YPH44tW7bg8OHDCAkJQfr06VGnTh18+umnaNy4sdbhERGRhevfvz9WrlyJkydPIjg4GG5ubqhQoQJ69OiRoB8YkbXjmiUiIiIiIiIjuGaJiIiIiIjIiDQxDS82NhZ37tyBu7t7sprUERERERGRbVBK4dmzZ/D29oad3ZvHjtJEsnTnzh34+vpqHQYREREREVmImzdvwsfH543XSRPJkru7OwD5h6RPn17jaIiIiIiISCuhoaHw9fU15AhvkiaSJf3Uu/Tp0zNZIiIiIiKiJC3PYYEHIiIiIiIiI5gsERERERERGcFkiYiIiIiIyAgmS0REREREREYwWSIiIiIiIjKCyRIREREREZERTJaIiIiIiIiMYLJERERERERkBJMlIiIiIiIiI5gsERERERERGcFkiYiIiIiIyAgmS0REROY0ZgwwbpzxfePGyX4iIrJITJaIiIjMyd4eGDUqYcI0bpxst7fXJq7UwmSRiKyYg9YBEBER2bSRI+XrqFGvftYnSgEBr/bbKn2yCMT/W+P+D4iILBSTJSIiInMbORK4f1+Sg9GjAaWAXLmAv/8G9uyRhMLeHnBwePv3Sb2eJXxvbw98/bX8D9JqskhEVk2nlFJaB2FuoaGh8PDwQEhICNKnT691OERElNZs3Qq0bQs8fy6Jkp0d0LEjEBMjl+johN8b2/au32tFp5NLbKwkUDExTJSISDOm5AYcWSIiIjKnH34ABg0C/PyAixcBJycgMlJ+Ts1kQSlJVlI6ATPl+/795atOB3Trlnp/OxFRMjFZIiIiMofoaEmSfvwRqFABOHDg1WiKfhoakHoJk073amqcFsaNk2TN0RGIigLy5wd++QVo1UqbeIiIkoDV8IiIiFJaaCjQuDEwdy7QpEn8RAmQrwEBxqvk2aK4a5QiI4Hhw4GXL4HWrYGPPwZevNA6QiIioziyRERElJKuXQMaNQJu3QL+/BPYtw8oXTrhCJL+Zy3XEqUGY8UcJk4EXFyk2MXSpcDevcCKFcAHH2gbKxHRa1jggYiIKKUcOgQ0bQq89x6weTNQqJDWEWlvzBiZ+mdsuuG4cVIl8OBB4PRpSaI++0wKYBARmYkpuQGTJSIiopSwcqUULfjwQ2D9eiBLFq0jsh6RkZJMTZkC1KoFLFkCeHtrHRUR2ShTcgOeuiEiInoXSskUs/btZQ3Ozp1MlEzl5ARMngzs2AGcPQsULw5s2qR1VERETJaIiIiS7eVLoHNnWXszbpysv3F21joq61WrFvDvv0ClSjKdsW9fICxM66iIKA1jskRERJQcDx/Kwf1vvwGrVgFffy3luendZMki0xjnzAEWL5Zpjf/+q3VURJRGMVkiIiIy1blzQNmywOXLwO7dQJs2WkdkW3Q6oE8f4J9/pC/Thx8CM2bIlEciolTEZImIiMgUO3ZIk9l06YAjRyRpIvMoXFgqDPbtKw1+GzSQ6nlERKmEyRIREVFSzZ0L1K8vydL+/UCuXFpHZPtcXIDvvwe2bgWOH5fiD3/8oXVURJRGaJ4s5c6dGzqdLsGlX79+AAClFMaMGQNvb2+4urqiWrVqOHv2rMZRExFRmhITAwweDHzyiYxybNoEsBVF6qpXT9YuffAB0LAhMHCgFNggIjIjzZOlo0eP4u7du4bLjh07AACtW7cGAEyZMgXTpk3DrFmzcPToUXh5eaF27dp49uyZlmETEVFa8ewZ0KwZMHMmMGuWfHVw0DqqtMnTE9iyRR6DefOAMmWk1DgRkZlonixlzZoVXl5ehsvmzZvx/vvvo2rVqlBKYfr06RgxYgRatGiBokWLYsmSJQgLC8OKFSu0Dp2IiGzdzZtA5cpAYKAcpP9/1gNpSKcD+veX9WIxMUDp0sDs2Sz+QERmoXmyFFdkZCSWLVuGjz76CDqdDkFBQbh37x7q1KljuI6zszOqVq2KAwcOJHo/ERERCA0NjXchIiIyydGjMnLx9Clw4IBMAyPLUby4VMvr0UOS2GbNgEePtI6KiGyMRSVLGzZswNOnT9GtWzcAwL179wAAnp6e8a7n6elp2GfMpEmT4OHhYbj4+vqaLWYiIrJBv/0GVK0K5M4NHD4MFC2qdURkjKurTI3ctEkS2uLFgZ07tY6KiGyIRSVLCxcuRP369eHt7R1vu+61Jn9KqQTb4ho+fDhCQkIMl5s3b5olXiIisjFKAZMmAa1aAU2bAn//LetkyLI1bizFH4oWBWrXBr78EoiM1DoqIrIBFpMsXb9+HTt37sTHH39s2Obl5QUACUaRHjx4kGC0KS5nZ2ekT58+3oWIiOiNIiOBjz4CvvoKGDUKWLFCRi7IOmTPDvz5J/Ddd9LAtlw54OJFraMiIitnMcnSokWLkC1bNjRs2NCwLU+ePPDy8jJUyANkXVNgYCAqVKigRZhERGSLHj+WEYkVK4Bly4CxY6WQAFkXOzvgs89k6mRYGODvD/z0E4s/EFGyWUSyFBsbi0WLFqFr165wiFOOVafTYdCgQZg4cSLWr1+PM2fOoFu3bnBzc0OHDh00jJiIiGzGpUsyCnHunEy769hR64joXZUqBRw7BnTqBPTsKdMqnzzROioiskIW0Shi586duHHjBj766KME+4YMGYLw8HD07dsXwcHBKFu2LLZv3w53d3cNIiUiIpuyaxfQsiXg5SWjEXnzah0RpZR06aQXU926wMcfS/GHZcuAatW0joyIrIhOKdsfmw4NDYWHhwdCQkK4fomIiMTPPwO9ewPVqwOrVwMZMmgdEZnLrVtA587SL2vYMJlm6eiodVREpBFTcgOLmIZHRESUamJjgaFDpT/Pxx9Ls1kmSrbNx0dKik+cCEydClSsCPz3n9ZREZEVYLJERERpx4sXsn5l6lTg+++B2bM5wpBW2NvLqNKBA0BwsKxrWrKExR+I6I2YLBERUdpw5w5QpQqwfbs0MR00iBXv0qIPPwSOH5ekuVs3oEMH4OlTraMiIgvFZImIiGzfiRNAmTLAgwfA/v1Ao0ZaR0RacncHFi0CVq4Etm4FSpYE9u3TOioiskBMloiIyLZt2gRUqiRNS48cAUqU0DoishRt2wKnTgG+vkDVqsCYMUB0tNZREZEFYbJERES2SSngu++AZs2A+vWlElr27FpHRZYmVy4pIT9mDDB+vCRNQUFaR0VEFoLJEhER2Z6oKCkL/sUXsqh/9WrAzU3rqMhSOTgAI0cCe/fK2raSJYEVK7SOiogsAJMlIiKyLcHBMpK0eLH0Upo4EbDjxx0lQfnywMmTQOPGQMeOQJcuQGio1lERkYb46UFERLbjyhU54D1xAtixA+jeXeuIyNp4eADLlgG//AJs2CAlxg8f1joqItIIkyUiIrINe/cCZctK09lDh2TtCVFydeoko0xZs0oT2wkTgJgYraMiolTGZImIiKzf0qVAzZpAsWKSKPn5aR0R2YK8eSUJHz5c1jTVqAHcvKl1VESUipgsERGR9YqNBb7+GujaFejcGdi2DciUSeuoyJY4OgLjxgG7d0uVvOLFgbVrtY6KiFIJkyUiIrJO4eFAu3ZSwGHyZOCnnwAnJ62jIltVpYr0ZKpdG2jdGvj4Y+D5c62jIiIzY7JERETW5949oFo1YPNm4LffgCFDAJ1O66jI1mXMCKxaJVUWV64E/P2BY8e0joqIzIjJEhERWZd//5VCDjdvynqS5s21jojSEp1OqiyeOCGV88qXB6ZMkSmhRGRzmCwREZH1+OMPqUyWKRNw5AjwwQdaR0RplZ8fsH8/8Nln0vi4Th3g9m2toyKiFMZkiYiILJ9SwMyZ0iy0Rg0ZUfLx0ToqSuucnIBvvgF27gTOnwdKlAA2btQ6KiJKQUyWiIjIskVHA59+CgwcCAweDKxbB7z3ntZREb1So4ZMD61cGWjWDPjkEyAsTOuoiCgFMFkiIiLLFRICNGoEzJ8vl2+/BezttY6KKKHMmSWRnzcPWLIEKF1aqucRkVVjskRERJYpKAioUEGazP75J9Czp9YREb2ZTgf06iUV8pycgDJlgHr1gIAA49cfNw4YMyZVQyQi0zBZIiIiy3PwoFS8e/lSkqWaNbWOiCjpChUCDh8G+vWTRsmjRwNDh8a/zrhxwKhRHCklsnBMloiIyLL8+itQvTpQoIAccBYsqHVERKZzdgamTZNR0XTppLx4586yT58oBQQAI0dqGycRvRGTJSIisgxKAWPHAh06AG3aSIWxLFm0joro3dStC1y9CuTPDyxbBjg6MlEisiI6pZTSOghzCw0NhYeHB0JCQpA+fXqtwyEiote9fAn06AGsWAGMHw989ZWs/yCyFUoBDg7SvNbBAYiK0joiojTLlNyAI0tERKStBw9kTdK6dcCqVcCIEUyUyPaMHy+Jkp3dq3L4RGTxmCwREZF2zp2TQg5XrgC7d8v0OyJbE3eNUkgI4OUF/PgjMGSI1pER0VswWSIiIvMaM0YOFl+3fTvg7y/NOw8flqSJyNa8XszhvfeAo0fl69Spso+ILBaTJSIiMi97ezkgjJswzZkj/WciIqR/Uq5c2sVHZE4xMQmLOfj4yEiqoyOwdq1MzyMii+SgdQBERGTj9AeJo0bJQWFwMDBjhmwbM0Z60BDZqsSazn7wAbB6NdCiBfD118DEiakaFhElDZMlIiIyv5EjgRcv4h84snQypXXNmslUvC++APz8gO7dtY6IiF7DZImIiMwrNhZYuBCYP//VNicnJkpEAPDZZ8ClS0CvXkDu3NKQmYgsBtcsERGR+Zw6BVSqJAeC+nVJTk5AZKTxog9EaY1OB8yaJUlSy5bAxYtaR0REcTBZIiKilPfsmZwx/+ADKZX80UfAyZMy9S4iQr6+XvSBKK1ydJT1S9mzAw0bAo8eaR0REf0fp+EREVHKUQr47Tdg4EDg6VNgwgQgPBwYOzb+GqW4RR/i/kyUVmXIAGzeLCX0W7QAduwAnJ21jooozePIEhERpYwrV4AGDYDWrYHSpaXh7NChss9YMYeRI2V7TEzqx0pkifLkATZuBI4ckZL6SmkdEVGap1PK9l+JoaGh8PDwQEhICNKnT691OEREtiUiApg8WUofe3oCP/wANGmidVRE1mvVKqBdO5mm+vXXWkdDZHNMyQ04DY+IiJJv506gb18gKEjKH3/9NZAundZREVm3tm2By5dl9PX994H27bWOiCjN4jQ8IiIy3d27cgBXu7YsSj91Cpg0iYkSUUoZMQLo3Fl6Lx04oHU0RGkWkyUiIkq6mBiZZlewIPDXX8CSJcDu3UDhwlpHRmRbdDpgwQKgTBlpXnv1qtYREaVJTJaIiChpjh6VA7eBA2VU6cIFoEsXOagjopTn7AysXw94eACNGkmFSSJKVUyWiIjozYKDZV1S2bJAbCxw8CAwdy6QKZPWkRHZvsyZgS1bgHv3pNJkVJTWERGlKUyWiIjIOKWAX36RKXfLlgHffy+jS2XLah0ZUdqSP7+MMAUGAv36saQ4USpiskRERAmdPw/UqCHT7KpXlyl3AwcCDiyiSqSJqlWB+fNlHdO0aVpHQ5Rm8FOPiIheCQsDxo8Hvv0WyJUL2LYNqFNH66iICAC6dZOS4l9+KSXFmzXTOiIim8eRJSIiEps3A0WKyFnrESOA06eZKBFZmnHjgFatgI4dgWPHtI6GyOYxWSIiSutu3ACaNwcaN5a1EadPA6NHAy4uWkdGRK+zs5OS/cWKyWv25k2tIyKyaUyWiIjSqqgoYOpUoFAh4PBhYNUq4M8/AT8/rSMjojdxdQU2bgScnCRhevZM64iIbBaTJSKitGjvXqBUKWDYMKBnTyng0KYNeyYRWQtPT5k6e/Wq9D2LidE6IiKbxGSJiCgtefgQ6N4dqFIFeO894J9/gOnTgfTptY6MiExVtCiwZo2MCH/+udbRENkkJktERGlBbKyUHC5YENiwAZg3DzhwQEaXiMh61a0L/PADMGMG8OOPWkdDZHNYOpyIyNadOgV88glw8CDQtSswZQqQLZvWURFRSvnkE+DSJWDAACkpXq+e1hER2QyOLBER2apnz4DPPgM++AAICQECA4HFi5koEdmib78FGjaUtYenT2sdDZHNYLJERGRrlALWrpUqd3PnAhMmACdOyDolIrJN9vbAihUystSoEXDvntYREdkEJktERLbkyhWgQQOgdWsZUTp/Hhg6VEoME5Fte+894PffgehooGlTICxM64iIrB6TJSIiWxARAQQEAEWKAOfOSQ+WjRuBXLm0joyIUpOPjyRMZ87IGsXYWK0jIrJqTJaIiKzdzp1AsWLAuHHAoEGSLDVponVURKQVf3+Zkvfbb8CIEVpHQ2TVmCwREVmru3elGWXt2kD27MDJk8A33wDp0mkdGRFprWlTKfrwzTfAzz9rHQ2R1WLpcCIiaxMTA8yeDXz9NeDsDCxZAnTuDOh0WkdGRJZk8GApKd67N5AnD1C9utYREVkdjiwREVmTo0eBMmWAgQNlVOnCBaBLFyZKRJSQTicNa6tXB1q0AC5e1DoiIqvDZImIyBoEBwN9+wJly8qC7QMHpCx4pkxaR0ZElszREVi9GvD2lj5Mjx5pHRGRVWGyRERkyZQCli0DChaUr99/L6NL5cppHRkRWYsMGYDNm4HQUKB5c6meSURJwmSJiMhSnT8P1Kgh65GqVZMpdwMHAg5cbkpEJsqTB9i0SU62fPyxnIghordiskREZGnCwoCvvgJKlABu3QK2bQNWrZJpNEREyVWunBSEWbYMGD9e62iIrAJPTxIRWZLNm4H+/aUs+IgRwNChgIuL1lERka1o2xb47z+pppkvnxSKIaJEMVkiIrIEN27IFLsNG4A6dYDt2wE/P62jIiJb9NVXUlK8e3cgVy6gQgWtIyKyWEyWiIjMacwYwN4eGDky4b5x44DISCB9ermeh4dMt2vdmqXAich8dDpg/nzg2jVpXnv4MJA3r9ZREVkkJktEROZkbw+MGiXfx02Yxo2T7dmySSnf/v2BgABJnIiIzM3ZGVi3TtYxNWwIHDwoVfOIKB4mS0RE5qRPkOImTMOGAZMny8958gB//gmUKqVNfESUdmXODGzZIglTq1bA1q3Sl4mIDHRK2X7tyNDQUHh4eCAkJATpedaWiLSgH0lycACio6Vow/TpQM+egB0LkxKRhgIDgdq1ga5dZXoepwGTjTMlN0hWshQaGopDhw7h9u3bCA8PR5YsWVC4cGEULVo02UGbE5MlItJcdLRMe4mNleTo7l2ZgkdEZAkWL5aCD1OnAl98oXU0RGZlSm6Q5Gl40dHRWLt2LebOnYv9+/cjNjYWcfMsnU6HzJkzo2PHjujbty/8WMWJiEi8fAn4+0uipB9ZmjfPeNEHIiItdOsGXL4MDBkCvP8+0Ly51hERWYQkzf3YtGkTChcujC5duiBdunSYOHEitm/fjlOnTuHixYs4ePAgli1bhnbt2mHDhg0oXLgw+vTpg0ePHpk7fiIiy/bsGVC4MHD+PNChAxAVJYUcRo2SqXlERJZi3DipxtmxI3DsmNbREFmEJE3Dy5gxIwYPHow+ffogWxKmjfz111+YMGECqlWrhlH6Rc0a4jQ8ItLE48dAiRLA7dvARx8BCxe+2qdfwxQQwBEmIrIc4eFA9erS++3wYcDXV+uIiFKcKblBkkaWgoKCMGrUqCQlSgBQs2ZN/P333xgwYECSrn/79m106tQJmTNnhpubG0qWLIljcc5oKKUwZswYeHt7w9XVFdWqVcPZs2eTdN9ERJq4dQuoXBkIDgb69ImfKAGSIAUEADEx2sRHRGSMqyuwcSPg5AQ0biyj40RpWJKSpQzJrLuflNsFBwejYsWKcHR0xNatW3Hu3Dl899138W47ZcoUTJs2DbNmzcLRo0fh5eWF2rVr4xlfwERkiS5dAipWBF68AE6cAObMMX69kSOlGS0RkSXx9JSS4kFBQPv2PKlDaVqSpuHNnj0brVu3RtasWVM8gGHDhmH//v3Yu3ev0f1KKXh7e2PQoEEYOnQoACAiIgKenp6YPHkyevfuneA2ERERiIiIMPwcGhoKX19fTsMjIvM7cQKoW1f6l+zYAfj4aB0REVHybNsmDWv79QNmzNA6GqIUk+LT8Pr374+goCDDz7GxsciZMyfOnDnzbpFCikeULl0arVu3RrZs2VCqVCksWLDAsD8oKAj37t1DnTp1DNucnZ1RtWpVHDhwwOh9Tpo0CR4eHoaLL+fbElFq2LMHqFYNyJUL2LuXiRIRWbe6dYEffgBmzgRmzdI6GiJNJClZen3wSSmFW7duITIy8p0DuHr1KubMmQM/Pz9s27YNffr0wYABA7B06VIAwL179wAAnp6e8W7n6elp2Pe64cOHIyQkxHC5efPmO8dJRPRGv/8uBxalSwN//w1kyaJ1RERE7+6TT4DBg4GBA4E//tA6GqJUl+Q+S+YSGxuL0qVLY+LEiQCAUqVK4ezZs5gzZw66dOliuJ7utW7SSqkE2/ScnZ3h7OxsvqCJiOJatkx6lDRpAqxYAbi4aB0REVHKmToV+O8/oG1bYP9+oHhxrSMiSjVJGlkyp+zZs6Nw4cLxthUqVAg3btwAAHh5eQFAglGkBw8eJBhtIiJKdTNnAp07A127AqtXM1EiIttjby8ngvLlAxo1AhKZ2UNki5KcLIWGhuLJkyd48uQJHj9+nGBb3IspKlasiIsXL8bbdunSJeTKlQsAkCdPHnh5eWHHjh2G/ZGRkQgMDESFChVM+l1ERClGKWD0aJma8sUXwE8/AQ6aD9YTEZnHe+/JdOOYGBlFDwvTOiKiVJHkT/a6desm2FazZk2j140xocTk4MGDUaFCBUycOBFt2rTBkSNHMH/+fMyfPx+ATL8bNGgQJk6cCD8/P/j5+WHixIlwc3NDhw4dkvx7iIhSTGysJEmzZgGTJgFDhwKJTAsmIrIZPj7A5s1ApUpAly4ymm6n+SQlIrNKUrI0evRoswXw4YcfYv369Rg+fDgCAgKQJ08eTJ8+HR07djRcZ8iQIQgPD0ffvn0RHByMsmXLYvv27XB3dzdbXERERkVFyfqkX38F5s0DevXSOiIiotRTqpS8/zVrBowYISeMiGxYkvosWTtTaqkTESUqLAxo0wbYvh1Yvhxo3VrriIiItPH998BnnwELFwIffaR1NEQmMSU34AR7IqKkePoUaNwYOH5cpqHE6f1GRJTmDBoEXLoE9O4N5MkDVK+udUREZsGJpkREb3P/vjSbPXsW+OsvJkpERDqdVAOtUQNo0QK4cEHriIjMgskSEdGbXLsmi5kfPAD27AHKldM6IiIiy+DoKEUevL2Bhg2BR4+0jogoxTFZIiJKzNmzQMWKUiZ8/36gaFGtIyIisiweHsCWLcDz51L0ISJC64iIUhSTJSIiYw4fBqpUAbJkAfbtkzn5RESUUO7cwMaNwD//AD16yAkmIhvBZImI6HU7dgA1awKFCgGBgYCXl9YRERFZtnLlgKVLpVLouHFaR0OUYkyuhnfjxo1E99nZ2cHDw4P9j4jIeq1dC3ToANSqJd+7uWkdERGRdWjTBvjvP+m/lC+fvJcSWTmTk6XcuXND95ZO9X5+fhg+fDi6du2a7MCIiFLdggVAnz5A27bA4sWAk5PWERERWZfhw6WkePfuQK5csu6TyIqZnCzNnz8fEydOhJubG9q0aQNPT0/cvXsXa9asQXh4OD755BPs2LEDH330EZycnNC+fXtzxE1ElLImTwaGDQP69gV++AGw4yxlIiKT6XTA/PlSSbRZM1n/mTev1lERJZtOKdNW4Y0aNQonTpzApk2b4o0wKaXQuHFjFCtWDJMmTULLli1x8+ZNHDlyJMWDNpUpXXqJKI1RChg6FJg6FRg5Ehg7Vj7siYgo+R4/BsqXB+ztgYMHgQwZtI6IyMCU3MDkU6eLFi1Cnz59EkzF0+l06N27N5YuXQoA6NixI86dO2fq3RMRpZ7oaKBnT0mUpk8HAgKYKBERpYTMmaWk+LVrwAcfAFFRCa8zbhwwZkxqR0ZkEpOTpUePHiE8PNzovpcvXyI4OBgAkDlzZpg4aEVElHoiIl6tTVqyBBg4UOuIiIhsi5+fFHm4elVGmeIeF44bB4waJSNPRBbM5GSpZMmSmDhxoiEp0nvy5AkmTJiAkiVLAgBu3rwJL5bbJSJL9OyZdJvfsgVYtw7o0kXriIiIbNPChUDz5sCxY0C9erJNnygFBMj0ZyILZnKBh6lTp6JOnTrIlSsXatSoAU9PT9y/fx9///03oqOjsXPnTgDAiRMn0Lhx4xQPmIjonTx+DDRoAJw/D2zbBlStqnVERES2bd06ea/dvh1wdJQp0EyUyEqYXOABAP7991+MHz8ee/bswePHj5E5c2ZUrVoVI0aMQPHixc0R5zthgQciAgDcugXUqQM8egT8+Sfg7691REREaUNsrCRKsbHSliEiQuuIKA0zJTcweWQJAIoXL47Vq1cnKzgiIk1cuiSJklLAvn1A/vxaR0RElHZMmCCJEgBERkrl0dGjtY2JKAneqZHIxYsXsX//frx48SKl4iEiSnknTgCVKgGurkyUiIhSW9w1Sn//LdvGjJHtRBYuWcnS0qVL4ePjg8KFC6NKlSq4ePEiAKBNmzZYsGBBigZIRPRO9uwBqlWTTvJ79wK+vlpHRESUdrxezKF6dWD4cGn8PWoUEyayeCYnS2vWrEG3bt3g7++PWbNmxSsP7u/vz+l5RGQ5Nm8G6tYFSpeWs5lZsmgdERFR2hITk7CYw9ixwIcfSqPasDDNQiNKCpOTpUmTJqF79+7YtGkTevXqFW9foUKF2IiWiCzDsmVAs2ZSqnbLFsDdXeuIiIjSnjFjEla9c3QEVqyQROrmTU3CIkoqk5Ol8+fPo127dkb3ZcqUCY8fP37noIiI3snMmUDnztI/ac0awMVF64iIiCiuvHmBuXOB5cuBX37ROhqiRJmcLLm5uSEkJMTovtu3byNjxozvHBQRUbIoJWcxBw4EPv9cmiE6JKvoJxERmVuHDnJSq29f4L//tI6GyCiTk6WKFSsmWKukt3jxYlSrVi0l4iIiMk1sLDBggMyFnzQJmDoV0Om0joqIiN5k1izAy0sSp8hIraMhSsDkZGnUqFE4dOgQypQpg5kzZ0Kn02HdunVo3Lgx9uzZgxEjRpgjTiKixEVFydnJH38E5s0Dhg1jokREZA3c3WX90okTCdc2EVkAnTI2RPQWu3btQt++fQ0lwwHAz88P8+bNs8iRJVO69BKRlQkLA9q0AbZvl6IObdpoHREREZlq6lRgyBBgxw6gVi2toyEbZ0pukKxkSe/KlSu4f/8+smTJgvwW3OSRyRKRjXr6FGjcGDh+HFi3TsqEExGR9YmNleqlp08D//4LZM2qdURkw1ItWbIWTJaIbND9+5Ic3bghpcHLl9c6IiIiehd37wLFiwNlywK//87p1GQ2puQGSSoTtXTpUpMC6NKli0nXJyIyybVrQO3awIsXwJ49QNGiWkdERETvKnt2YPFioFEjKfzQv7/WERElbWTJzi5+HQjd/zP9uDfVxcn+Y2JiUiq+FMGRJSIbcvYsUKcO4Ooqc9vz5NE6IiIiSkkDB0qxnsOHgRIltI6GbJApuUGSquEFBQUZLgcPHoSvry8+/vhj7Nq1C+fPn8euXbvQo0cP+Pr64sCBAynyRxARJXD4MFClCpAlC7BvHxMlIiJbNHkyULAg0K6dFPEh0pDJa5bat28PLy8vfP/99wn2DR48GHfu3MGqVatSLMCUwJElIhuwYwfQvDlQsiSweTOQIYPWERERkbmcPw988AHQubOMMhGloBQfWYpr69ataNiwodF9DRo0wLZt20y9SyKiN1u7FmjYUEaVtm9nokREZOsKFQJmzADmz5dqp0QaMTlZio2NxeXLl43uu3z5MtJAcT0iSk0LFgBt2wKtWgEbNgBublpHREREqeHjj4GWLeXrzZtaR0NplMnJUr169TBixAhs2bIl3vbNmzfj66+/Rl32OSGilDJ5MtCrF9CnjzScdXLSOiIiIkotOp2cMHvvPaBjR8DCCohR2mBysjRjxgx4eXmhSZMmyJAhAwoUKIAMGTKgadOmyJYtG2bMmGGOOIkoLVFKOrkPGwaMHCklZO1MfrsiIiJrlzEjsHw5sH8/MGGC1tFQGpSsprQvX77E4sWLsXv3bjx+/BiZM2dG9erV0aVLF7i6upojznfCAg9EViQmBujdG1i4EPj+e2DQIK0jIiIirY0eDYwfL731KlbUOhqycqbkBslKlqwNkyUiKxERIVMtNmwAfv4ZYINrIiICgOhooFo1Wbt06hQL/dA7MWs1PCIis3j+XLq2b94slY+YKBERkZ6Dg0zHCwmR2Qe2f66fLESSkqWiRYti/fr1Sb7Tu3fvYsCAAfjmm2+SHRgRpSGPHwM1a0rT2W3bgCZNtI6IiIgsTa5cUvBh9Wpg0SKto6E0IknJUps2bdClSxfkzJkTw4cPx7Zt2/Dw4UNDmfDw8HCcOXMGP/30Exo3boxcuXLh2LFjaMIDHiJ6m9u3pX9SUBCwezdQtarWERERkaVq3Rro0QPo3x+4eFHraCgNSPKapbt372L69On4+eef8fjxY+h0Ouh0Ojg6OiIyMhIAoJRC5cqVMXDgQLRo0cKsgZuCa5aILNTly0Dt2jKdYvt2oEABrSMiIiJL9+IF8MEHgKsrcOgQ4OysdURkZcxa4CEqKgqHDx/GwYMHcefOHYSHhyNLliwoWLAgqlWrBh8fn3cK3hyYLBFpaMwYwN5eSoDHdeIEULky4OIi3/v6ahIeERFZoRMngHLlgH79gGnTtI6GrIwpuYGDqXfu6OiISpUqoVKlSskOkIjSEHt7YNQo+V6fMO3dKyNKEREylYKJEhERmaJUKWlcPngwUKcOUK+e1hGRjTI5WSIiMok+QdInTKVKAc2bSxnYESOkbwYREZGpBg6UKdxduwL//gt4emodEdkgJktEZH4jR8q6JH3CBMj3Y8dqFxMREVk3nQ5YvBgoXlwSpj/+AOzYFYdSFp9RRGR+L18C16+/+tnJiYkSERG9u2zZgKVLpe3E999rHQ3ZICZLRGReN29KafAlS+RnJycgMhIYN07buIiIyDbUqQN8/jkwfDhw7JjW0ZCNYbJEROYTGCjlXS9dAmJigIAAKeoQECDT8JgwERFRSpg4UabjtW8PPH+udTRkQ5gsEVHKUwqYOROoWRN47z0gJEQSJH2xh5EjmTAREVHKcXICfv0VuHMHGDBA62jIhiSrwMPDhw8xbdo07N69G48fP8b69etRpEgRzJs3D2XKlEGpUqVSOk4ishZhYUDv3sCyZcBnnwHp0gGOjgn7LOl/jolJ/RiJiMj2+PkBs2YB3bvL1Lx27bSOiGyAyclSUFAQKlasiJCQEJQoUQJXrlxBREQEAODff//FoUOHsGjRohQPlIiswLVrQIsWwIULwPLlQIcOb77+6wkUERHRu+jaVYo99O4NlC0L5MmjdURk5UyehjdkyBBkyJABly9fxp49e6CUMuyrVKkS9u/fn6IBEpGV+OsvoHRp4OlT4ODBtydKREREKU2nA+bOBTJlAjp2lJ5+RO/A5GTpr7/+wujRo+Ht7Q2dThdvX/bs2XHnzp0UC46IrIBSwLffypSHDz4A/vkHKFFC66iIiCit8vAAVqwAjhxhmwp6ZyYnSy9fvkSmTJmM7nvx4gXs2AyMKO148UIqD335JTBkiDQETOT9gYiIKNWULy+J0oQJUpmVKJlMzmwKFCiAnTt3Gt23Z88eFC1a9J2DIiIrcOWKfBht3gysXg1MmgTY22sdFRERkRg2DKhaFejUCXjyROtoyEqZnCz17NkTM2bMwIwZMxAcHAwAiIyMxNq1azF79mz07t07xYMkIgvz55+yPik8HDh0CGjdWuuIiIiI4rO3B375Raq0fvyxTBsnMpFOKdOfOb169cJPP/0EOzs7xMbGws7ODkop9OzZE3PnzjVHnO8kNDQUHh4eCAkJQfr06bUOh8h6KSUjSF9/DdSvLxXvMmTQOioiIqLEbdgANG8OzJkD9OmjdTRkAUzJDZKVLAHAoUOHsGXLFty/fx9ZsmRBo0aNUKFChWQFbG5MlohSwLNnQLduwLp10kx29GiAaxSJiMgafPIJsHixFCEqUkTraEhjZkuWwsPDkS9fPsydOxeNGzd+50BTC5Mlond08aKclbt1S6Y0NG2qdURERERJFxYGfPihTM07cgRwcdE6ItKQKbmBSaeFXV1dER4ejnTp0r1TgERkRX7/HShTBoiNlQ8YJkpERGRt3NyAlSuBS5ekgitREpk8h6ZmzZqJVsMjIhsSGytlV5s0AWrUkESpYEGtoyIiIkqeYsWA774DZs2SE4FESeBg6g2++uortGzZEi4uLmjRogWyZ8+eoDltYn2YiMhKhIQAnTtLWfBx44CvvuL6JCIisn59+wLbtwPduwOnTgE5cmgdEVk4kws8xG06+3qSpBcTE/NuUaUwrlkiMsG5c7I+6f596YDeoIHWEREREaWcR4+AEiVktsT27ewRmAaZkhuYPLI0atSoRJMkIrJy69YBXbsCuXNLxaB8+bSOiIiIKGVlySLFimrVAqZOlea1RIlIdulwa8KRJaK3iImRcuATJ0qD2Z9/Bt57T+uoiIiIzGf4cODbb4H9+6WQEaUZZquGR0Q2KDgYaNQI+OYbYPJkYNUqJkpERGT7AgIAf3+gfXsgNFTraMhCmTwNLyAg4I37dTodRo4cmeyAiCgV/fuvrE96+hT480+gdm2tIyIiIkodjo7Ar78CJUsC/frJ1Dyi17xTgQejd6jTscADkTVYtQr46CPAzw9Yvx7Ik0friIiIiFLf8uVAp07A0qVSCZZsnlmn4cXGxia4PHr0CD/99BOKFi2Ka9euJTduIkoN0dHSkK9dOxlVOnCAiRIREaVdHTtKktS3L/Dff1pHQxYmRQs8fPPNNzh8+DDWr1+fUneZIjiyRPR/jx5JkrR7tyxqHTgQYHVLIkolt24Bly/LgLaPj9bREMXx7BlQqhSQKROwbx/g5KR1RGRGmhV4KFOmDP7666+UvEsiSiknTgClS0sTvh07gEGDmCgRUapZuBDIlQuoUUO+LlyodUREcbi7y/qlEyekOizR/6VosnTq1Cm8xypaRJbnl1+AChWkt8SxY0D16lpHREQ2TCkgPBy4cwc4e1aWRfbsCcTGyv7YWKB3bxlpIrIYH34ITJgATJkC7NypdTRkIUyuhrd06dIE2yIiIvDvv//i559/RqdOnUy6vzFjxmDs2LHxtnl6euLevXsAAKUUxo4di/nz5yM4OBhly5bFjz/+iCJFipgaOlHaExUFfPEFMHMm0K0bMHs24OqqdVREZAViY2VmUnCwFMwMDn51iftzYt9HRr75/mNiZHkIp+ORRfniC5l90bmzVIzNmlXriEhjJidL3bp1M7rdxcUFnTp1wrfffmtyEEWKFMHOOBm8vb294fspU6Zg2rRpWLx4MfLnz4/x48ejdu3auHjxItzd3U3+XURpxoMH0mD2wAHgxx+BTz7htDuiNCYqShKY5CQ7ISGvRoJe5+4OZMwIZMggXzNmBAoWfPV93O0ZMwIvX0pngtfvb9o0ua8PPjDjP4HIFHZ2UhWveHGge3fg99/52ZnGmZwsBQUFJdjm4uICT0/P5Afh4AAvL68E25VSmD59OkaMGIEWLVoAAJYsWQJPT0+sWLECvXv3TvbvJLJpR48CLVrIkdKuXUClSlpHREQwvcCBfjpbYgnN2xKf58+N36+dXcKkJlMm4P33E25/PfHx8AAcTD56AObPl6l3MTGAvT3Qti1w6JAspaxeXYp01qvH41KyANmzA4sWAY0bA7NmAf37ax0RacjktzudTofs2bPD0dExwb7o6GjcuXMHOXPmNOk+L1++DG9vbzg7O6Ns2bKYOHEi8ubNi6CgINy7dw916tQxXNfZ2RlVq1bFgQMHEk2WIiIiEBERYfg5lF2ZKS35+Wcpf1qyJPDbb0COHFpHREQAfvpJkoXYWEkIunQBihV7e+KT2HQ2Z+f4SUyGDICvr5wQTyzZ0X/v7p76SUmPHkDdujL1Ll8+SRZjYoB164CpU4EGDYAiRWQWVIcOLEZGGmvUSJKkL78EqlaVFxalSSaXDre3t8fBgwdRpkyZBPuOHTuGMmXKmNSUduvWrQgLC0P+/Plx//59jB8/HhcuXMDZs2dx8eJFVKxYEbdv34a3t7fhNr169cL169exbds2o/dpbB0UAJYOJ9sWGSkV7ubMAXr1knVKzs5aR0WU5ikFLFsmydHr0qUDMmdOOHrztmQnY0bAxSV1/w5zUgrYu1eSps2bAW9v6WzQu7eMZBFp4uVLoGxZmaXxzz+Am5vWEVEKMaV0uMkjS2/KrWJiYqAz8VRV/fr1Dd8XK1YM5cuXx/vvv48lS5agXLlyAJDgPpVSb/w9w4cPx2effWb4OTQ0FL6+vibFRWRV7t4FWrWSN/P586XsFBFpKjYW2LhRimsdO2b8Ops3A9WqpWpYFkmnA6pUkcu5c8B33wEjRwLjx8u5n4EDZdSMKFW5uEg58dKlgcGDgXnztI6INJCs0uHGEpWIiAhs3boVWbJkeaeA0qVLh2LFiuHy5cuGdUz6ynh6Dx48eOMaKWdnZ6RPnz7ehchmHTggq6OvXQMCA5koEWksOhpYsUJm7bRoAbz3HrB8uawTisveXqajUXyFC0sPpmvXgH79ZPpi3rwyMvfvv1pHR2lO4cLA9OlyInLdOq2jIQ0kKVkaO3Ys7O3tYW9vD51Oh3Llyhl+1l/c3NwQEBCApk2bvlNAEREROH/+PLJnz448efLAy8sLO3bsMOyPjIxEYGAgKlSo8E6/h8jqKSVnuapVk1XZx44B/x+NJaLUFxkpSwYLFQI6dpSRkL17gd27ZQ3O/PmSIAHydd48ls1+k+zZgUmTgJs3ZXpeYCBQooSse9q5U94CiVJFz55y5uPjj+UJSWlKktYsbd26FX/88QeUUpg9ezZatWqVYGTH2dkZxYoVQ4cOHYwWf0jMF198gcaNGyNnzpx48OABxo8fj8DAQJw+fRq5cuXC5MmTMWnSJCxatAh+fn6YOHEidu/ebVLpcFPmJRJZhZcvgU8/ldOv/fpJ/V2uhibSRHi4JEmTJ8txVPPmwIgRxsth37oVv8ABJV1UFLBmjSROJ09KDZsvv5QOCSYcdhAlz5Mnkq3nzQv8/ferMx9klUzJDUwu8NC9e3eMGjUKefLkeacg9dq1a4c9e/bg0aNHyJo1K8qVK4dx48ahcOHCAF41pZ03b168prRFixZN8u9gskQ25dYtoGVL4NQpYO5caTZLRKnu+XN5CX77LfDwIdCuHTB8OGDCxxMlg1JyrDp1KrBtG5Azp9S2+fhjqfJHZDZ79kid+zFjZFEdWS2zJkvWiMkS2Yw9e+Q0qrOzzJ0uXVrriIjSnKdPgR9+kGUMoaFA167AsGFcf6SFf/+VZPXXX2VtWJ8+wIABMoWPyCxGj5aqLYGBQMWKWkdDyZQqydKZM2dw/vx5hIeHJ9jXxVh9VA0xWSKrp5Q0xvvsM2kwu3o1kDWr1lERpSkPHwLffy8vxchIWcbw5ZcyskHaunkTmDFD1oVFRACdOgGffy5r84lSVHS0rBW+eVNmeGTIoHVElAxmTZbCwsLQpEkT/P3339DpdIZS4nEr5JnSZyk1MFkiqxYeLqdLly6V0qVTpgAOJlf9J6Jkun1bRi/mzZOKdn37ynmL/xdsJQsSEiKP04wZwJ07QMOGktBWqZL6TXjJhl2//qrayMqVfHJZIVNyA5NLh48bNw7Xrl1DYGAglFJYt24dduzYgRYtWsDPzw/Hjx9PduBE9Jrr12Ukac0aqT08bRoTJaJUEhQk5yny5gUWL5aD7uvX5XwFEyXL5OEBDBkij93ixfJ4VasmfUXXrJFBAaJ3liuXDGOuXg0sWqR1NGRmJidLGzduxNChQw2lu3PmzImaNWtizZo18Pf3x5w5c1I8SKI06e+/pZzWkyfSS6lDB60jIkoTLlyQdUh+fsBvvwFjx8pB99ixQObMWkdHSeHkJI/hv/8Cf/wh65natAHy5wd+/BEIC9M6QrJ6bdoAH30E9O8PXLyodTRkRiYnS9euXUPBggUNPZfC4rzjdOzYERs2bEjJ+IjSHqWkfX3t2oC/P/DPP1Ijl4jM6tQpOf4pXFj6+Hz7rTRGHTYM4Axu66TTAfXry7mnf/6REaYBA2Sd2ahRwIMHWkdIVm3mTOkB0L69LJYjm2RyspQhQwa8ePECAJAtWzZcvnzZsC8qKsqwj4iS4cULGUH64guZ87N1K09lE5nZ4cNAkyZyTuLoUWDOHODqVSlHnS6d1tFRSvngA6mad+WKNA3+7juZTdWnDxDnUIYo6dKlkzVLZ89K3wCySSYnS8WKFcOlS5cAANWrV8fEiROxb98+HDlyBAEBAShRokSKB0mUJly9ClSoAPz+u8yD/uYbNr0jMhOlgN27ZQC3XDk5WF6yBLh0CejdW6rzk23KnVsKQNy8CXz9NbBhA1CgANCiBXDwoNbRkdUpVUo+r7//HvjzT62jITMwOVnq0aMHnj17BgCYMGECwsLCULVqVZQvXx7Xr1/Hd999l+JBEtm8bdukZ9KLF8ChQ9JLiYhSnFIyYFu5svSWfPhQzk2cOQN06QI4OmodIaWWTJmAESNkquX8+cC5c3K+qmJFSaBiY7WOkKzGwIFAvXqyUO7+fa2joRT2zk1pX7x4YSgjXqFCBWTKlCmlYksxLB1OFkspOSM1YoRMrF+2DMiYUeuoiGxObKwcAE+YABw/LmtXvv5aSkuz6i8B8hzZvBmYOhXYt0+KQXz+uSTRLi5aR0cW7/59KSdesqRUFbEzeTyCUpHZSoeHh4ejQ4cO2Ldvn2FbunTp0LhxYzRq1MgiEyUizY0ZA4wbl3D7s2dA0aLAV1/JUdvvvzNRIkph0dFSdb9YMaBlSyktvXOnTLdq1IiJEr1iZydr1/buledH0aKynilXLnkLf/xY6wjJonl6ylzebduA6dO1joZSkEnJkqurKzZu3IhYjk0TJZ29vZRdipswXb4MvP++zPto3x4ICOBZKKIUFBkJ/PQTULAg0KmTHPDu2ydV0WrWZJJEb1aunJSNv3hRkuyJE6WCXv/+0sOJyKi6daVj9bBhMoRNNsHko7OSJUvizJkz5oiFyDaNHCnJkD5h2rxZTnM/fCifvCtWaB0hkc0IDwdmzQLy5QN69pRZMceOyayYihW1jo6sjZ8fMHs2cOOGFCj99Vd5brVtK6XIiRKYOFE+49u1A54/1zoaSgEmr1nau3cvOnfujCVLlqBq1armiitFcc0SWYSBA6Ung96IEcD48drFQ2RDnj0D5s6VctAPH8qA7fDhQJEiWkdGtiQsDFi8GJg2TUqQV6smSVS9epwcQHFcuiR9Etu0AX7+WetoyAhTcgOTk6VixYrh7t27CA4ORsaMGZE9e3bo4sxn0Ol0OHXqVPIiNxMmS6Sp8+clKVq58lV5JScnNrAjSgHBwcAPP8gSgefPpRjV0KFy9p/IXGJipGDI1KnSp6twYWmP16EDy87T/y1aBHz0kXz2t22rdTT0GrMVeACAzJkzo1ixYqhSpQqKFSuGLFmyIHPmzIYLizwQ/d/ZszIMX6SIrBhu0EC2OznJggpjRR+IKEkePJCRo1y5gEmTZF3SlSvAggVMlMj87O1lLdPBg8CePbIE9aOPgDx5gMmTgadPtY6QNNetmyRJvXpJfXqyWu9cOtwacGSJUtXp05IIrV0L+PpKtbs7d2TdUkCArGEaN07WMOl/JqIkuXUL+PZb6Ytjbw/07SvrqT09tY6M0rrz52Ua6C+/yDmxXr1k9nXOnFpHRpp5+lS6IKdPL43nHRzi7x83ToYpx4zRILi0zawjS0SUiH//BVq1AooXB44eBebNk6p3Dx7ET5SAhEUfiOiNrl4FeveWM/hLlgBDhgDXr8tZfCZKZAkKFZIKjNevAwMGyFKV99+XUU/96oRbt4Bdu+QrpQEZMkiT+Zs3gVq14u/TnzS1t9ckNEq6ZCVLDx8+xPDhw1G+fHn4+fnh7NmzAIB58+bhxIkTKRogkcU7cQJo3lzKbp04ASxcKIs7e/aU04sxMcZHkPQJU0yMNnETWYHz56UpaP78wPr18pK5fl1OxHLWN1kiLy9pfnzzpoyC7tsnfUqLFJFpozVqyNeFC7WOlFLFggXyoAcGAj16yDbOLrEqJk/DCwoKQsWKFRESEoISJUrg8OHDOHr0KPz9/dGvXz+EhYVh0aJF5oo3WTgNj8zi2DF5o9u0SU4ffv010LEj4OiodWREVu/kSTng/O03wNtbRpI+/hhwc9M6MiLTREdLpcb+/eNvt7eXpSw+PpqERakpJkYWU1679mrdMhMlTZl1Gt6QIUOQIUMGXL58GXv27EHcXKtSpUrYv3+/6RETWZN//gEaNwZKl5bT3kuWABcuyGJOJkpE7+TgQaBRI6BUKTkfMW+eFG4YMICJElknBwfjJexjYoD//kv9eEgD9vZS6AmQRMnRkYmSFTE5Wfrrr78wevRoeHt7xysZDgDZs2fHnTt3Uiw4Ioty5AjQsCHw4YeyFumXX4Bz52SO0OuLNokoyZSSdRw1awIVKsj6pF9+eTWblaWYydr5+Rnvw7R6tRw7Uxqgn3Wl0wFRUUDTpvLmRxbP5GTp5cuXiZYHf/HiBezYlY1szcGD0nGwbFk5ilu+XMqCd+rEJInIRHEXuCsF/PEHUKmSTOl/8kSKSJ45w5cX2RYfn1cVHIFXpccXLAAqVpTzb2TD4q5RCguTmSmbNslitufPtY6O3sLkzKZAgQLYuXOn0X179uxB0aJF3zkoIouwfz9Qp46c6r51SxrLnTkjXQdZvYbIZAsXvlrgnjOnVNRt2FCSpi1bgOPH5QCS59zIFvXoIUtWdu2Sr2vXAgcOSHXpUqVk4IEDDTbo9WIOLi5SMbdVK6mimyePnIAli2XyR1LPnj0xY8YMzJgxA8HBwQCAyMhIrF27FrNnz0bv3r1TPEiiVLV3r5T4rFQJuHtX5kn8+680l2OSRJQst25J35nYWPlZKeDGDWDVKjkv0aCBzE4hsmU+PkC1aq+KOnz4oRRRbdtWmtq2awf8/9CKbEViFXHXrJGqHzodUKYMsHSpNvHRWyWrKW2vXr3w008/wc7ODrGxsbCzs4NSCj179sTcuXPNEec7YTU8SpLAQGDsWDntV7y4nAlq3pynuYlSwMKFUs3udbt2ycEjUVq3erWcUEifXmZ7V66sdUSUKsLCgH79gMWLZfjxhx8AV1eto7J5puQGyUqWAODQoUPYvHkzHjx4gCxZsqBRo0aoUKFCsgI2NyZLlCilgN27JUkKDJT5w6NHA02aMEkiSgFhYXJS9dtvE7YUY+lkovhu3JD1evv3AyNGyDk7rt1LIxYtkqTJz09GnfLn1zoim5YqyZI1YbJECSgF/PWXHMXt3Qv4+0uS1Lgx5wIRpZA//wT69gXu3JE2ZFmyAJ9+KkmTvb2UBdf3aCQiERMDTJokjZfLlJFRpjx5tI6KUsXp00Dr1sDt28BPP8n8TDILU3KDZJ+vOHjwIHbt2oXHjx8jc+bMqFatmsWOLBEZKAXs2CEjSQcOSEWa33+XVeZMkohSxL17wODBUhOlRg1g2zY5WQpID6X//pP+jBxRIkrI3l5OLtSqJfWESpQA5syRnudk44oVk+IPffrIArbAQGDaNCkKQZoxOVkKDw9Hu3btsHnz5ngNaXU6HRo0aIDVq1fDlXMtydIoJUdsY8cChw7J6botW4D69ZkkEaWQ2FgphTx0qPRcXLpUphTFfYn5+DBJIkqKcuWAkydlZlanTjJS++OPsqaJbJi7O7BsGVClCjBwIHD4sEzLy5tX68jSLJMXZQwZMgRbt27F+PHjERQUhPDwcAQFBWHcuHHYtm0bhgwZYo44iZJHX5O4XLlXidGff0rCxPJbRCnmzBlZkN6nj5T/vnAB6NyZLzGid5E+vTRoXrYM2LhRltUePKh1VGR2Oh3Qu7c82CEhslRg/Xqto0qzTE6WVq1ahZEjR2L48OHIlSsXnJ2dkStXLnz11VcYMWIEVq5caY44iUyjlEyvK1NG5v04OgLbt8uq2bp1eQRHlELCw4GvvpI+MY8fS72UhQuBzJm1jozIdnTsKKNMnp5yUmLcuIQFU8gGlSoFHDsmczJbtJD5zZGRWkeV5picLIWFhSW6NqlixYoIDw9/56CIkk0pOf1WurRUtHN1BXbulCIOtWszSSJKQdu3A0WLAt99Jy1ETp0CqlbVOioi25Q3r3yUffWVFH+oXl2q55GN8/CQaXgzZ8o8zCpVgOvXtY4qTTE5WSpXrhyOHj1qdN/Ro0dRpkyZdw6KyGSxsTJE7e8PNGsmc3537QL27AFq1mSSRJSC7t+XM9116wK5cknP5lGjAGdnrSMjsm0ODlLEdfduOV4uXlz6M5GN0+mkge2+fVJBp1QpYPNmraNKM0xOlmbOnIl58+bhxx9/RPD/20wHBwdj1qxZmD9/Pn744YcUD5IoUbGxwNq18sbRogWQKZN8iuzezU6XRClMX8ChYEGpl7J4sVTgL1BA68iI0pbKlWUkt25dqS790UfA8+daR0VmV6YMcPw4UKmStDoZNgyIjtY6Kptncp8ld3d3REVFISoqCgDg4OCA6P8/UI6OjnBycnp15zodQkJCUjDc5GGfJRukT5LGjZOV5bVqSZ+kSpW0jozIJp07J+uN9+0DunUDpk6VvklEpB2lgCVLpH+ZtzewYoXMQicbp5TMfx42DChfXvo05MihdVRWxax9llq2bAkdpzSRVmJiZO7uuHFy9FanjnS2ZI8vIrMIDwcmTACmTJHGmH//LWsliEh7Op2cvKhYUXoylS8PjB8PfPklYGfy3CGyGjod8MUX8oC3bStlEpcvl2MiSnEmjyxZI44s2YCYGGDVKkmSLlyQMuCjRklJcCIyi507pRT4zZvA8OFyEpO9EYksU2SkTLCYPFlOaCxdysGGNOHhQ+nTsH27dDMePVo6G9MbmZIb8LwDWbboaGkwUaSIrCh//31p0PbHH0yUiMzkwQP57K1dWxrInjol1beYKBFZLicnYNIkWUd48aIUf2BrnjQga1Y5Jho3TqYB1K4tRSAoxSRrZCkiIgI7d+7E9evX8fLly/h3qNNh8ODBKRZgSuDIkhWKjpYh5QkTgMuXZSHjqFGcjE1kRkoBixbJFB4A+PZbmeLDmddE1uXxY6BnT0mWevUCpk0D0qXTOioyu927gfbt5fsVKzhn+g1MyQ1MTpaOHTuGxo0b4/79+zB2U51OhxgL65TGZMmKREXJSNKECcCVK0DTppIk+ftrHRmRTTt/Xqbc7dkjo0rffScnLInIOikl1SsHDZIS/ytWSOFYsnH37skCtsBAqTM/fDgXsBlh1ml4ffv2Rfr06bF+/XpcvHgRQUFB8S5Xr15NduBko8aMkeFhY8aNk/1RUcDChVKD+KOPgBIlgBMngA0bmCgRmdHLl3I+okQJ4M4dWae0dCkTJSJrp9PJqNLx4zKFtmxZGWGKjdU6MjIrLy9gxw5gxAjpFt6gAfDokdZRWTWTq+GdPXsWK1asQJMmTcwRD9kie3s5GgPkhas3bpxsb9IE8POTDnutWkmCVLy4JqESpSV//y2jSdeuSfGGr77iuiQiW1OwIHDokBw7f/659EhbskSOqclG2dvLqFKlSrLeu2RJKZJVsaLWkVklk0eWcubMaY44yJaNHCkv2lGjXo0wjR4tP3t4AL//Lqe8Tp+WsuBMlIjM6uFDoGtXoGZNwNMTOHlSXqJMlIhsk7OzrEHctg3491/5mN2yReuoyOzq1JE3+Dx5gKpV5Ulg+0WwU5zJydKQIUPw7bffIiIiwhzxkK2KmzA5Osr3gAwPnzkjZzyKFtU2RiIbpxSweDFQqJCco/jpJ5nWXriw1pERUWqoU0eSpbJlgUaNgP79pZca2bAcOYBdu6Qv05dfylrwJ0+0jsqqmDwNr1u3brh27Rref/99VKtWDZkyZYq3X6fTYcaMGSkWINmIly+BjBnl++hoWWx49qzMDyAis7t4EejdW5Kjjh1l7UK2bFpHRUSpLWtWYNMmYPZsOX7evRv49Veer7RpDg7AN9/ItLwuXWQt+OrVQJkyWkdmFUyuhrdlyxa0aNECUVFRxu+Q1fAorvBwYP586ZJ3756c2nZ0lIIOAQHx1zARUYqLiJDeK5MmAb6+wJw50oaDiOjMGak0ffmyzNDq14+tAmze9etA27ZS+ePbb2V4MQ0+6Gathvfll1/C398fJ0+eREREBGJjY+NdLC1RIo2EhQHffw/kzSsrSj09JVEKCJA246+vYSKiFLd7t6xNmDBBZl+cPs1EiYheKVoUOHpUqub17y8tDR880DoqMqtcuaRHRL9+wMCBQOvWQEiI1lFZNJOTpWvXrmHMmDEoXrw4HB0dzRETWbMXL6RBS968cnTWoIG8A+tXkOtHkowVfSCiFPHoEdC9u/QjzJpVXn7jxwOurlpHRkSWxsUFmDkT2LwZOHJETrBs26Z1VGRWTk5yQvu336RfxAcfSLsWMsrkZKlgwYIIDQ01RyxkzZ4/B6ZOlYorw4bJ6anLl6V3koeH8Sl3+oSJo5FEKUIpKQlcsKBU4J8/X04gFimidWREZOkaNpTiDyVLAvXqAZ99JtN4yYa1aAEcOybHaeXLA3PnslqeEclas/TVV19h27Zt8LKSIv1cs2RGz54BP/4oo0khIXI6e/hwIHdurSMjSlMuXZKeSbt2SfP2adNk9isRkSliY2WkaehQqZz566/ylWzYy5eSHc+ZI4vY5s0D3N21jsqsTMkNTE6WmjRpgpMnT+LJkycoWbKk0Wp4GzduND1qM2KyZAahocCsWZIkPXsG9OghSRL7cBGlqogIqZ8yYYJUiJ0zB6hbV+uoiMjanTwpJ16uXZMZW716pck6AGnLypVAz57yYbJmDVCsmNYRmY1Zk6XcuXND94ZXi06nw9WrV025S7NjspSCQkKAH36Q09YvXsiLauhQKbNFRKlqzx4pB/7ff7JE8OuvATc3raMiIlsRFiY1mubOBZo1k95smTNrHRWZ1cWLUvThv/9k5lD37lpHZBZmTZasEZOlFPD0KTBjBjB9upQD79VLkqQcObSOjCjNefwYGDIE+PlnoEIFmTHBHilEZC4bNsgEEhcX4JdfgBo1tI6IzCosDBgwQNadd+smSZONnYkza+lwSmOCg4HRo2UN0jffAF27AlevyoRmJkpEqUopYNkyKeDw229ytnfvXiZKRGRezZpJ8YeCBYFataSOU2Sk1lGR2bi5yTDikiWvmteeP691VJpJVrIUERGBefPmoX379qhduzYuX74MANi4caPFTcGjZHryRKrV5c4tVe569JAkafp0wNtb6+iI0pzLl6VHUufOcrBy4YJMwbPjKS8iSgU5cgA7dsgaye++k1HtS5e0jorMqksXqScfGwt8+CGwYoXWEWnC5I/ZR48eoXTp0vjkk08QGBiIv//+G8+ePQMAbNiwAd9++22KB0mp6PFjYMQISZKmTZPpdkFB8s6YPbvW0RGlOZGRUryhWDHgyhXgjz+kOpWVFCMlIhtiZyfrIw8elDpP/v4yHdj2F3SkYUWKSMLUvDnQsaOcpXv5UuuoUpXJydKQIUPw9OlT/PPPP7hx4wbiLnmqXr06AgMDUzRASiWPHr0q+T1jBvDJJ5IkTZ3K+sNEGtm3T3qejB4NDBoEnD0L1K+vdVRElNaVLg0cPw60bSsTT9q2lVn7ZKPeew9YuhRYsECm5pUvLwUg0giTk6XNmzcjICAA/v7+Cari+fj44NatWykWHKWCBw9kpXju3FIK/NNPpU7o5MlAtmxaR0eUJj15IoUmK1eWXoHHj8uSQRtbX0tEVuy992T9/+rVMj2vRAlZQ0k2SqcDPv4YOHQIeP5chhXXrtU6qlRhcrIUGhqKXLlyGd0XFRWF6Ojodw6KUsH9+8AXXwB58sgq8YEDJUmaNAnIkkXr6IjSJKVkSnihQnIAMns2sH8/ULy41pERERnXujVw6pQcTlSrJsudo6K0jorMpmRJ4NgxoF49efAHDJCGfzbM5GQpT548OHjwoNF9R44cQYECBd45KDKje/ekS3OePDKc+vnnkiRNmMDmCUQaunJFmsl27CgHHBcuyGxYFnAgIkuXMyfw999AQICcc61SRWpCkY1Knx5YtUpmJM2dK9MggoK0jspskvQxvGfPHjx//hwA0LFjR0yePBkbN240rFfS6XQ4evQoZsyYgc6dO5svWkq+O3dk0UOePLIac8gQSZICAoBMmbSOjijNioyUg4uiRaWy1JYt8hnEeipEZE3s7aU+1L59MsO/ZElpdUA2SqcD+vUDDhwAHj6UaXmbNmkdlVkkqSmtvb09Dh48iDJlyiAqKgpNmjTBtm3bkDFjRgQHByNLlix4/Pgx6tWrh99//x12FnYqNE03pb19WxY7LFgAuLoCgwfLkGmGDFpHRpTm7d8vhYUuXJAB39GjgXTptI6KiOjdhIbKEuhffpHR8h9/lPWXZKOCg4Hu3YGNG2WJx8SJgKOj1lG9kSm5gUNS7jBuPuXo6Ig//vgDq1atwpYtW3D//n1kyZIFjRo1Qrt27SwuUUqzbt6UJOmnn+To6+uvgf79+W5FpKFbt6RfUtaswA8/APPnS6+/Y8dkcTQRkS1In16Kp9WrJ9OJ9++X9Zi+vvIe6OcH+PhoHSWlmIwZgfXrge+/B4YOldGmlSvlAbcBSRpZsrOzw6FDh1CmTJnUiCnFpamRpRs3ZE7PwoXybvX55zJMaut/N5GFW7hQ2pbFxsrPzs7SvqxPH5m+QkRki4KCgE6dpDcTIIVs7OzkZFGPHtrGRmZw8CDQpg0QHi7zMOvV0zoio0zJDZI8DPR6mXCyMNeuyXyefPmklOO4cfIONXw4EyUijd26FT9RAoDoaKBpUyZKRGTb8uQBli+X7/Wn52Nj5ZCF3WZsUPnywIkTMm2ifn0p/jB2rPHrjhsHjBmTquElR5Km4QHScDYpU+x0Oh1CQkLeKSgywdWrMjd0yRIZBp0wQca833tP68iI6P8WLIifKAFATIz09ONUFCKydUFBrxIlPb4H2rAsWYDNm6Vnp77qx4sXwJQpr64zbhwwapQUGrNwSU6WqlWrhqxZs5ozFjLFlSuSGC1dKiW/v/lG5vNwdTiRxXj8WBY5r1yZcJ+9vQwEExHZOj8/mXr3+kmj1auBihUtvhYAJYedncxuqlABaNQImDpVquYtWhQ/URo5UutI34prlqzN5cuSJC1bJqvEhw6V+T1ublpHRkRx/P67vDQjIqQS1IsXcj4jJkYSpXnzOF+fiNKOhQtl6p3+PbB5c2DDBuDDD+WEUs6cWkdIZnP/viRNV68CDg4yD13jRMmU3IDJkrW4eBEYP17KyXh6AsOGAT17SjlwIrIYISHS0mzxYqBhQ1nE7O0t+27dkmkn+fJx6gkRpT2vvwcePAi0awc8eybvmU2aaB0hmU1MjFQ2iokBnJzkTKKGzFLggTRy/rw0KShcGNi1C5gxQzLzAQOYKBFZmJ07gWLFgN9+k97Pv//+KlEC5OCgWjUmSkSUNr3+Hli+PHDyJFC1qhS8GTxYGnWTDZo48VWiFBkpU/GsBJMlS3XuHNC+PVCkCLB3LzBrlqxT+vRTwMVF6+iIKI7nz4G+fYHatWVu/unT0p+PRUSJiN4sY0Zg3To5F/zjj7KG6epVraOiFBV3jVJEhHwdNcpqEqYkFXiIfX1FHpnPmTPy5FmzRpp5zZkDdOsmQ5dEZHH27pWX6L17ck7jk09kXSsRESWNTicTZipUANq2BUqVkjVOrVppHRm9M2PFHPRfR42K/7OFSnI1PHoHY8bIakZjT4Zx42RYskULeSL99huQO7es/u7aVYYricjihIcDX38tDcsrVAC2bWN1OyKid1G6NHD8uBTHad1aTj5Nm8YJNVYtJsZ4MQf9zzExqR+TiSzq/OekSZOg0+kwaNAgwzalFMaMGQNvb2+4urqiWrVqOHv2rHZBJoe9vfHhRn22vXo1UKKEvEP89BNw6ZIUb2CiRGSRjhwB/P1lysiUKUBgIBMlIqKU4OEh1fHmzpW1n+XKyWERWakxYxIfORo50iqa0lpMsnT06FHMnz8fxYsXj7d9ypQpmDZtGmbNmoWjR4/Cy8sLtWvXxrNnzzSKNBlGjkw4P/OTT14NP0ZGSt35ixelljAbDhBZpMhIGU2qUEH6Ph8/DnzxhZwPISKilKHTSZnxI0eAly/l5NTy5VpHRWmVRSRLz58/R8eOHbFgwQJkzJjRsF0phenTp2PEiBFo0aIFihYtiiVLliAsLAwrVqxI9P4iIiIQGhoa76K5uAmTvb2cMsmUCViyBLhwQRY9MEkislinTkk/kMmT5UTYgQNSpJKIiMyjeHHgn39kpUKnTsDHHwNhYVpHRWmNRSRL/fr1Q8OGDVGrVq1424OCgnDv3j3UqVPHsM3Z2RlVq1bFgQMHEr2/SZMmwcPDw3Dx9fU1W+wmGTlSmnHFxsrX+/eBLl3keyKySNHR0gf6ww8BpYCjR2V0iec2iIjM77335LzyokXAr7/Ke7G1rcYg66Z5srRy5UocP34ckyZNSrDv3r17AABPT8942z09PQ37jBk+fDhCQkIMl5s3b6Zs0Mk1bpwceTk5yVcjfzMRWY7z52XK3ahRwJdfSqJUsqTWURERpS06nUzAOXpUvv/wQ1nPpJTWkVFaoGmydPPmTQwcOBDLli2DyxtKnehea1ailEqwLS5nZ2ekT58+3kVzVl5jnigtiYkBvvtOyteGhsqUuwkTWMGfiEhLhQvLOqaOHWWJd+fOgDUtYSfrpGmydOzYMTx48AAffPABHBwc4ODggMDAQMycORMODg6GEaXXR5EePHiQYLTJoiVWY54JE5HFuXJFOsx/+aU0mj1xAihbVuuoiIgIANzcgAULpODDxo1SbvzUKa2jIlumabJUs2ZNnD59GidPnjRcSpcujY4dO+LkyZPImzcvvLy8sGPHDsNtIiMjERgYiAoVKmgYuYneVGM+IMAqaswT2brYWGD2bFlQfPs2sHu39PdwddU6MiIiel2HDsCxY5I8lS0rdbM4LY/MQdPKAu7u7ihatGi8benSpUPmzJkN2wcNGoSJEyfCz88Pfn5+mDhxItzc3NChQwctQk6eN9WQt/CuxURpwY0bMqVj506gTx9g6lRZVExERJYrf37g4EHg88+lI8vff8uok4eH1pGRLbH4MmxDhgxBeHg4+vbti+DgYJQtWxbbt2+Hu7u71qERkZVTCli8GBg0CEifHti2DYhTfJOIiCyci4s0CK9eXU56+fsDq1bJ9DyilKBTyvYHLUNDQ+Hh4YGQkBDLKPZARJq7exfo1QvYvBno2hWYPh3IkEHrqIiIKLmuXgXatQNOnpQZAgMGSPU8oteZkhtoXjqciCg1KQWsXAkULSplaDdulNElJkpERNYtb15g3z6gf3+ZMdCsGfDkidZRkbVjskREacajR0DbtkD79kDt2sCZM0CTJlpHRUREKcXJSVo/bNoE7N0rLSAOHtQ6KrJmTJaIKE3YsAEoUgT46y8ZWVq5EsiSReuoiIjIHBo3lul4Pj5A5crAlClS9ZTIVEyWiMimPX0KdOkCNG8OlCsHnD0ro0tERGTbcuaUNhBffgkMHQo0agQ8fKh1VGRtmCwRkc3atk3WJunXJW3YAHh5aR0VERGlFkdHYNIk4M8/gX/+AUqWBAIDtY6KrAmTJSKyOc+eAb17A/XqAYULy9qkrl1ZFYmIKK2qW1em5eXPD9SoAYwbB8TEaB0VWQMmS0RkUwIDgeLFgeXLpaP7tm2Ar6/WURERkda8vaX5+MiRwOjRkkDdu6d1VGTpmCwRkU0IC5NSsdWqyTz1f/+V0SWOJhERkZ69PTBmjCRNZ88CJUrI90SJYbJERFbv0CEpDztvHvD998CuXdJvg4iIyJgaNWRaXokSQJ06wNdfA9HRWkdFlojJEhFZrYgIYPhwoGJFaSp74oSMLtnxnY2IiN7C01MKP0yYAHzzjSRQt25pHRVZGh5SEJFVOnECKF1amg+OHw/s3w8ULKh1VEREZE3s7OSk2+7dwNWrUi3vjz+0joosCZMlIrIqUVFAQABQpgzg4CClYIcPl++JiIiSo1IlmZZXrhzQsKH0ZoqK0joqsgRMlojIapw9C5QvL8nS8OHA4cNS+Y6IiOhdZckC/P67zFiYPh2oXBm4dk3rqEhrTJaIyOLFxABTpwL+/lL17uBBSZicnLSOjIiIbIlOB3z2GbBvn5QVL1VKGppT2sVkiYgs2uXLQJUqwNChwIABwPHjwIcfah0VERHZsrJlZW1sjRpA8+bAwIFSVIjSHiZLRGSRYmOBWbOkrOv9+8CePTK65OKidWRERJQWZMwIrF0L/PCDNDmvUAH47z+to6LUxmSJiCzO9etArVpA//7ARx8Bp07J4lsiIqLUpNMBn34q079DQmQ6+KpVWkdFqYnJEhFZDKWAn34CihWTs3c7d8roUrp0WkdGRERpmb+/TANv2BBo1w7o3RsID9c6KkoNTJaIyCLcuSMfQj17Aq1bA6dPAzVrah0VERGRSJ8eWLECmD8fWLpU1jVduKB1VGRuTJaISFNKAcuXA0WLymLazZuBhQsBDw+tIyMiIopPp5OTekeOSB+m0qWBX37ROioyJyZLRKSZBw+AVq2ATp2AevWAM2dkdImIiMiSFSsmTdFbtQK6dAG6dwdevNA6KjIHJktElGpu3QJ27ZKv69bJaNKePcCaNTK1IXNmrSMkIiJKmnTpgMWL5bJ6tbS1OH1a66gopTFZIqJUsXAhkCuX9Kzw9QVatgQqVpTRpFattI6OiIgoebp2lVEmBwegTBkpVKSU1lFRSmGyRERmd+sW0KuX9E7Ss7MDZs4EPD21i4uIiCglFCoEHD4sU/J69gQ6dgSePdM6KkoJTJaIyKxOnQI++SR+ogTIz1euaBMTERFRSnN1BebNA379Ffj9dyk3fuKE1lHRu2KyREQpLjJSPiwqVQJKlgSOHpUKQnHZ2wP58mkSHhERkdm0ayc9mdzdgXLlgB9/lGl5cdftkvVgskREKebWLWDkSCBnTqBDB8DJCVi7Frh5E1iwQBIkQL7Omwf4+GgbLxERkTn4+QEHD0rz2k8/lRLj+nW7uXLJOl6yDjqlbH8JWmhoKDw8PBASEoL06dNrHQ6RTVEK+PtvYPZsYONGwM1N5mz37QsULhz/urduAf/9JyNKTJSIiCgtWLBA1u3GZW8PXLvGz0KtmJIbOKRSTERkY0JCpIP57NnSwbxwYSnY0LmzTD0wxseHHwxERJS2GJtyHhMjJw/5mWj5mCwRkUlOn5YE6ZdfgIgIoHlzYO5coEqVhOuSiIiI0jo/P6kA+3qho4AAwMMDKFVKm7goabhmiYjeKjISWLVKEqLixWW63RdfANevSyO+qlWZKBERERnj4wPMnx9/3W6vXjI13d8faNNGZmiQZWKyRESJun0bGD1aFqO2aydnxlatkiRpzBjA21vrCImIiCxfjx6yRmnXLvk6bx5w7pwUejh0CChSBOjeXfaRZWGBByKKRylg926Zard+vfSN6NxZCjYULap1dERERLYlIkJGnsaPB4KDZdRpxAgge3atI7NdpuQGHFkiIgBAaKj0gihaVEqbnj0LTJ8uo0uzZzNRIiIiMgdnZ6B/f+DqVVnHtHw58P77wNChwOPHWkdHTJaI0rizZ4F+/YAcOYCBA4FChaQU+Nmz0huCg7FERETmly4dMGwYEBQEfP65nMDMmxcYNw549kzr6NIuJktEaVBUFLBmDVCtmowY/fYbMHiwzJVeuxaoXp0FG4iIiLSQIYMkSFevAh99BEyYIEnTtGlAeLjW0aU9TJaI0pC7d4GxY4HcuaX6TmwssHIlcOOGDP2z3wMREZFlyJYN+P574PJladMxZIj0bJo3T056UupgskRk45QCAgOBtm2BnDmBqVOBxo2BU6eAPXtku5OT1lESERGRMb6+UgDiwgWZEfLJJ0DBgtLvMCZG6+hsH5MlIhv17BkwZ470RapWTZKj776Tgg1z58p2IiIisg758knxh1OngGLFgC5d5LN8/Xo5MUrmwWSJyMacPy9VdXLkkAINfn7Azp2yfcAA6RZORERE1qlYMWDDBunP5O0NtGgBlCkDbN/OpMkcmCwR2YDoaCnSUKMGULgwsHq1JEbXrgHr1gE1a7JgAxERkS0pWxbYsUMq2Do6AnXrykySffu0jsy2MFkismL37knFnNy5gVatgMhIGaK/cUOa2/n6ah0hERERmVP16sD+/cDmzUBICFC5MtCwIXDihNaR2QYmS0RWRik5a9S+vRRs+OYboEEDeVPctw/o0EEa3BEREVHaoNNJgnT8uFS5/e8/wN9fKt9euKB1dNaNyRKRlXj+XMqFliwpZ42OHQOmTJGCDfPny3YiIiJKu+zspMrt2bPAwoXA4cNAkSJA9+4yNZ9Mx2SJyMJdvAgMHCgFG/r2BfLkkUWcFy4AgwZJ8zoiIiIiPQcHaWh76RIwfTrwxx9A/vxS+OnuXa2jsy5MlogsUHS0lAKtVUt6Kfz6K9Cvn3Tz3rABqF1bzh4RERERJcbZWSrkXr0qzeeXLwfefx8YOhR4/Fjr6KwDD7eILMj9+8CECUDevFIKNCxMms7dvAlMnAjkyqV1hERERGRt0qUDhg0DgoKAzz8HfvxRjjUCAqQvIyWOyRKRxpQCDhwAOnaU6nUTJgB16siapAMHgE6dWLCBiIiI3l2GDFJF9+pVoEcPORGbN680rQ8P1zo6y8RkiSgV3boF7NolX1+8ABYskGo1FSsCR45IZbvbt4GffpLtRERERCktWzZg2jSpmteihUzLy5cPmDsXiIrSOjrLwmSJKJUsXCjT6GrUkJLfWbMCvXvLaNKff0ohh88+AzJm1DpSIiIiSgt8fKTS7oUL0tC2b19ZK/3LL0BMjNbRWQYmS0Rm9OCBJEJDhwIffwzExsp2pYCXL6WJ3KZN0nWbBRuIiIhIC/nySfGHU6eA4sWBLl3k67p1csySljloHQCRLVAKuHNHmsHpL8eOyZQ6QBZWGrtNRETqxklERESUmGLFpBrvkSPAiBFAy5ZA6dKynrp2bWl+m9YwWSIykVLS2C1uYnT8uIwiAUCWLLLeqEsX+ervDzg6ArlzvxpZAgB7ezmTQ0RERGRJypQBduyQddYjRsgMmCpVJGmqVEnr6FIXkyWiN4iNlcWPx47FT4yePpX93t6SDPXp8yox8vExfuZl/nxZoxQTI4nSvHlyXSIiIiJLVL26LBn44w9JmipXBurXB8aPTzuFqHRK2f5MxNDQUHh4eCAkJATp06fXOhyyUNHRssAxblJ04gTw/Lnsz537VULk7w+UKgV4eZn2O27dkuQrXz4mSkRERGQ9YmOBtWuBkSOBS5eA1q2lT1PBglpHZjpTcgMmS5QmRUQAZ8/GT4xOnZKiCwDg5/cqKfrgA0mMMmXSNmYiIiIirUVHS7W8MWPkJHCXLsDo0XJS2VowWXoNk6W0LTwc+Pff+InR6dPSR8DODihUKP6IUcmSAJ8mRERERImLiJAlBhMmAE+eAL16yVS97Nm1juztmCy9hslS2vHsGXDyZPzE6Px5WSfk4AAULfpqtMjfX8piurlpHTURERGRdXrxAvjhB2DKFJmh078/MGQIkDmz1pEljsnSa5gs2abgYFlTFLf4wuXLUq3O2RkoUSL+iFHRorKdiIiIiFLW06fAd98B338vhaw+/xwYPBhwd9c6soSYLL2GyZL1e/AgYanuoCDZly6drCmKmxgVLCjluomIiIgo9Tx4AHzzDTB7tiRKw4YBffsCjx/LSW0/P+2LXDFZeg2TJctx69abXyhKSSPX1xMjfXNXD4/4SZG/v9yXvX3q/h1ERERElLhbt4Bx44CFC4H33gNCQ+U4z85O1jr16KFdbEyWXsNkyTIsXCiL/2Jj5YUybx5Qs+abm7vq1xbpL3nypM3u0URERETWaO9eaWgbl709cO2adiNMpuQGbEpLqSIo6FWiBMjXnj1f7Y/b3FWfIOXIwcSIiIiIyJpFRyfcFhMjfSe1no6XFEyWKMU9eiQ9i06dkpLdp04BZ868SpTi+uYboGtX05u7EhEREZHl8/OTGUVxjwPt7YF8+bSLyRRMlijZoqOBixdfJUT65OjOHdnv6goUKwaULg20aCHNy15/oXTsyESJiIiIyFb5+Mgapd69ZUTJ3l6WYljDqBLAZImS6MmThKNFZ89KQzIA8PWVUt3du0vvohIl5IxB3MIL3t7W+0IhIiIiouTp0QOoW1em3uXLZ13HfyzwQPHExEi1On1ipE+Obt2S/S4u0q9InxCVKCGjR5kyJe3+b92yzhcKEREREdkGFnigJAkOfjVKFHdt0cuXsj9HDkmGOnd+lRz5+QEO7/Cs8fFhkkRERERE1oHJUhoQEwNcuZJwtOjGDdnv5AQUKSLJUMeO8rV4cSBzZm3jJiIiIiLSEpMlGxMSYny0KCxM9mfPLslQu3avptHlzw84OmobNxERERGRpWGyZKViY4GrVxOOFl27JvsdHYHChSUZatv21TS6rFk1DZuIiIiIyGowWbICz55JIhS3RPfp08CLF7Lf01OSoVatXo0WFSgg0+uIiIiIiCh5NE+W5syZgzlz5uDa/4dEihQpglGjRqF+/foAAKUUxo4di/nz5yM4OBhly5bFjz/+iCJFimgYdfLduiXV5vz8EhY6iI2VkaHXR4uuXpX9Dg5AoUKSDLVo8Sox8vRM9T+DiIiIiMjmaZ4s+fj44JtvvkG+/7fxXbJkCZo2bYoTJ06gSJEimDJlCqZNm4bFixcjf/78GD9+PGrXro2LFy/C3d1d4+hNs3Ah0KuXJEV2dsDQoUCuXPFHi549k+tmySKJULNmrwouFCoEODtr+icQEREREaUZFtlnKVOmTJg6dSo++ugjeHt7Y9CgQRg6dCgAICIiAp6enpg8eTJ69+6dpPuzhD5Lt25JYhQbG3+7nd2r0aK4vYu8vACdTpNQiYiIiIhsltX2WYqJicGaNWvw4sULlC9fHkFBQbh37x7q1KljuI6zszOqVq2KAwcOJJosRUREICIiwvBzaGio2WN/m8uXEyZKALB1KxDnzyMiIiIiIgthp3UAAHD69Gm89957cHZ2Rp8+fbB+/XoULlwY9+7dAwB4vrYox9PT07DPmEmTJsHDw8Nw8fX1NWv8SeHnJ6NIcdnbS8U6IiIiIiKyPBaRLBUoUAAnT57EoUOH8Mknn6Br1644d+6cYb/utfloSqkE2+IaPnw4QkJCDJebN2+aLfak8vEB5s+XBAmQr/PmJSzyQERERERElsEipuE5OTkZCjyULl0aR48exYwZMwzrlO7du4fs2bMbrv/gwYMEo01xOTs7w9kCKyH06AHUrQv89x+QLx8TJSIiIiIiS2YRI0uvU0ohIiICefLkgZeXF3bs2GHYFxkZicDAQFSoUEHDCJPPxweoVo2JEhERERGRpdN8ZOmrr75C/fr14evri2fPnmHlypXYvXs3/vzzT+h0OgwaNAgTJ06En58f/Pz8MHHiRLi5uaFDhw5ah05ERERERDZM82Tp/v376Ny5M+7evQsPDw8UL14cf/75J2rXrg0AGDJkCMLDw9G3b19DU9rt27dbXY8lIiIiIiKyLhbZZymlWUKfJSIiIiIi0p4puYFFrlkiIiIiIiLSGpMlIiIiIiIiI5gsERERERERGcFkiYiIiIiIyAgmS0REREREREYwWSIiIiIiIjKCyRIREREREZERTJaIiIiIiIiMcNA6gNSg77sbGhqqcSRERERERKQlfU6gzxHeJE0kS8+ePQMA+Pr6ahwJERERERFZgmfPnsHDw+ON19GppKRUVi42NhZ37tyBu7s7dDqd1uGkaaGhofD19cXNmzeRPn16rcOhVMbHn/gcSNv4+KdtfPzTNkt6/JVSePbsGby9vWFn9+ZVSWliZMnOzg4+Pj5ah0FxpE+fXvMXCmmHjz/xOZC28fFP2/j4p22W8vi/bURJjwUeiIiIiIiIjGCyREREREREZASTJUpVzs7OGD16NJydnbUOhTTAx5/4HEjb+PinbXz80zZrffzTRIEHIiIiIiIiU3FkiYiIiIiIyAgmS0REREREREYwWSIiIiIiIjKCyRIREREREZERTJaIiIiIiIiMYLJEqWLSpEn48MMP4e7ujmzZsqFZs2a4ePGi1mGRRiZNmgSdTodBgwZpHQqlktu3b6NTp07InDkz3NzcULJkSRw7dkzrsCgVREdH4+uvv0aePHng6uqKvHnzIiAgALGxsVqHRmawZ88eNG7cGN7e3tDpdNiwYUO8/UopjBkzBt7e3nB1dUW1atVw9uxZbYKlFPemxz8qKgpDhw5FsWLFkC5dOnh7e6NLly64c+eOdgEnAZMlShWBgYHo168fDh06hB07diA6Ohp16tTBixcvtA6NUtnRo0cxf/58FC9eXOtQKJUEBwejYsWKcHR0xNatW3Hu3Dl89913yJAhg9ahUSqYPHky5s6di1mzZuH8+fOYMmUKpk6dih9++EHr0MgMXrx4gRIlSmDWrFlG90+ZMgXTpk3DrFmzcPToUXh5eaF27dp49uxZKkdK5vCmxz8sLAzHjx/HyJEjcfz4caxbtw6XLl1CkyZNNIg06dhniTTx8OFDZMuWDYGBgahSpYrW4VAqef78Ofz9/TF79myMHz8eJUuWxPTp07UOi8xs2LBh2L9/P/bu3at1KKSBRo0awdPTEwsXLjRsa9myJdzc3PDLL79oGBmZm06nw/r169GsWTMAMqrk7e2NQYMGYejQoQCAiIgIeHp6YvLkyejdu7eG0VJKe/3xN+bo0aMoU6YMrl+/jpw5c6ZecCbgyBJpIiQkBACQKVMmjSOh1NSvXz80bNgQtWrV0joUSkWbNm1C6dKl0bp1a2TLlg2lSpXCggULtA6LUkmlSpXw119/4dKlSwCAU6dOYd++fWjQoIHGkVFqCwoKwr1791CnTh3DNmdnZ1StWhUHDhzQMDLSSkhICHQ6nUXPNHDQOgBKe5RS+Oyzz1CpUiUULVpU63AolaxcuRLHjx/H0aNHtQ6FUtnVq1cxZ84cfPbZZ/jqq69w5MgRDBgwAM7OzujSpYvW4ZGZDR06FCEhIShYsCDs7e0RExODCRMmoH37/7V370FV1A0Yx5/1AAcxLgMpiI6BomGmJlpppWnqvF6wSRO81IiXLuYlE2WkmopuZHazhtTXV/CSYepAWTI1qZhmMVOjoY4zWhZS77yjlIaOIoiHff9gPHliNUU4e5DvZ4ah/Z3d33mO+0c889vdM8HuaPCyo0ePSpIiIyM9xiMjI1VaWmpHJNiosrJS6enpmjhxokJCQuyOc0mUJXjdrFmztG/fPu3atcvuKPCS3377TXPmzNGXX36pwMBAu+PAy2pqatSnTx9lZmZKknr16qUDBw5o6dKllKVmYP369Vq7dq1yc3PVrVs3FRcX66mnnlJ0dLRSUlLsjgcbGIbhsW2aZp0xXN+qq6s1fvx41dTUaMmSJXbHuSzKErxq9uzZ+vTTT7Vz5061b9/e7jjwkt27d6usrEy9e/d2j7lcLu3cuVNZWVmqqqqSw+GwMSEaU9u2bXXLLbd4jHXt2lV5eXk2JYI3paWlKT09XePHj5ckde/eXaWlpXrttdcoS81MVFSUpNoVprZt27rHy8rK6qw24fpVXV2t5ORklZSUqLCw0KdXlSTuWYKXmKapWbNmKT8/X4WFhYqNjbU7Erxo8ODB2r9/v4qLi90/ffr00UMPPaTi4mKK0nXu7rvvrvNVAT/++KNuuukmmxLBmyoqKtSiheefGw6Hg0eHN0OxsbGKiorSli1b3GPnzp3Tjh07dNddd9mYDN5yoSj99NNP2rp1qyIiIuyO9I9YWYJXzJw5U7m5udq0aZOCg4Pd1y2HhoaqZcuWNqdDYwsODq5zf1qrVq0UERHBfWvNwNy5c3XXXXcpMzNTycnJ+u6777R8+XItX77c7mjwglGjRunVV19Vhw4d1K1bN/3www96++23NXXqVLujoRGcPn1ahw8fdm+XlJSouLhY4eHh6tChg5566illZmaqc+fO6ty5szIzMxUUFKSJEyfamBoN5XLnPzo6WmPHjtWePXu0efNmuVwu99+D4eHhCggIsCv25ZmAF0iy/Fm5cqXd0WCTe++915wzZ47dMeAln332mXnrrbeaTqfTjI+PN5cvX253JHjJqVOnzDlz5pgdOnQwAwMDzY4dO5rPPvusWVVVZXc0NILt27db/v8+JSXFNE3TrKmpMV944QUzKirKdDqd5oABA8z9+/fbGxoN5nLnv6Sk5JJ/D27fvt3u6JfE9ywBAAAAgAXuWQIAAAAAC5QlAAAAALBAWQIAAAAAC5QlAAAAALBAWQIAAAAAC5QlAAAAALBAWQIAAAAAC5QlAAAAALBAWQIA+IxVq1bJMAwZhqGvvvqqzuumaSouLk6GYWjgwIGNluPbb79VRkaGysvL67wWExOjxMTERntvAIDvoCwBAHxOcHCwsrOz64zv2LFDP//8s4KDgxv1/b/99lu9+OKLlmUJANB8UJYAAD5n3LhxysvL06lTpzzGs7Oz1a9fP3Xo0MGmZACA5oSyBADwORMmTJAkrVu3zj128uRJ5eXlaerUqXX2P3HihGbMmKF27dopICBAHTt21LPPPquqqiqP/QzD0KxZs/TBBx+oa9euCgoKUs+ePbV582b3PhkZGUpLS5MkxcbGXvKywC+++EIJCQlq2bKl4uPjlZOT01AfHwDgIyhLAACfExISorFjx3oUkHXr1qlFixYaN26cx76VlZUaNGiQ1qxZo9TUVBUUFOjhhx/WokWLNGbMmDpzFxQUKCsrSy+99JLy8vIUHh6u0aNH65dffpEkPfLII5o9e7YkKT8/X0VFRSoqKlJCQoJ7jr1792revHmaO3euNm3apB49emjatGnauXNnY/xzAABs4md3AAAArEydOlWDBg3SgQMH1K1bN+Xk5CgpKanO/UqrV6/Wvn37tGHDBiUlJUmShg4dqhtuuEELFizQli1bNHToUPf+Z8+e1datW93zJCQkKDo6Whs2bFB6errat2/vvsyvV69eiomJqZPtjz/+0DfffOPeb8CAAdq2bZtyc3M1YMCAxvjnAADYgJUlAIBPuvfee9WpUyfl5ORo//79+v777y0vwSssLFSrVq00duxYj/HJkydLkrZt2+YxPmjQII/CFRkZqTZt2qi0tPSKs912220e900FBgaqS5cuVzUHAMD3sbIEAPBJhmFoypQpeu+991RZWakuXbqof//+dfY7fvy4oqKiZBiGx3ibNm3k5+en48ePe4xHRETUmcPpdOrs2bNXnK0h5gAA+D5WlgAAPmvy5Mn6448/tGzZMk2ZMsVyn4iICB07dkymaXqMl5WV6fz587rxxhu9ERUAcB2iLAEAfFa7du2UlpamUaNGKSUlxXKfwYMH6/Tp0/rkk088xtesWeN+/Wo5nU5JYqUIAJo5LsMDAPi0hQsXXvb1SZMm6f3331dKSoqOHDmi7t27a9euXcrMzNSIESM0ZMiQq37P7t27S5LeffddpaSkyN/fXzfffHOjfxkuAMC3sLIEAGjSAgMDtX37dj300EN64403NHz4cK1atUrz589Xfn5+veYcOHCgnn76aX322We65557dPvtt2v37t0NnBwA4OsM8+8XeQMAAAAAWFkCAAAAACuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACwQFkCAAAAAAuUJQAAAACw4Gd3AOB65XK5VF1dbXcMXCF/f385HA67YwAAAB9CWQIamGmaOnr0qE6ePCnTNO2OgytkGIZCQ0MVFRUlwzDsjgMAAHwAZQloYCdPnlR5eblat26tVq1a8Yd3E2Caps6cOaPff/9dLVu2VFhYmN2RAACAD6AsAQ3INE2VlZUpJCREN954o91xcBVatmypqqoqlZWVKTQ0lJILAAB4wAPQkFwul1wul0JCQuyOgnoICQlxn0MAAADKEtCAzp8/L0ny82PRtim6cN4unEcAANC8UZaARsAlXE0T5w0AAFyMsgQAAAAAFihLAAAAAGCBsgTgurBq1SoZhvGPPzExMXZHBQAATQR3oQO+LCNDcjik556r+9rLL0suV+0+0MiRI1VUVOQx1q9fP40dO1bz5s1zjzmdTm9HAwAATRRlCfBlDof0/PO1/31xYXr55drxl16yJ5cPat26tVq3bl1nPDIyUn379rUhEQAAaOq4DA/wZc89V1uInn++tiBJnkXJasWpkWRkZMgwDO3bt09JSUkKDQ1VeHi4UlNTdf78eR06dEjDhg1TcHCwYmJitGjRIo/jKysrNW/ePN12223uY/v166dNmzZ57PfRRx/JMAxlZWV5jL/wwgtyOBzasmVLo39WAAAAiZUlwHsqKqSDB6/+uJEjpf/976/CVF0tTZ9eO75nT/2yxMdLQUH1OjQ5OVkPP/ywHn/8cW3ZskWLFi1SdXW1tm7dqhkzZmj+/PnKzc3VggULFBcXpzFjxkiSqqqqdOLECc2fP1/t2rXTuXPntHXrVo0ZM0YrV67UpEmTJEnjx4/Xjh07NG/ePPXt21d9+vRRYWGhXnnlFT3zzDMaOnRo/T4zAADAVTJM0zTtDgFcLyorK1VSUqLY2FgFBgZ6vrhnj9S7tz3B/m73bikh4aoOycjI0Isvvqi33npLqamp7vFevXqpuLhY+fn5Gj16tKTaL3WNjo5W//79lZeXZzmfy+WSaZqaPn269uzZoz0XFb+qqir169dP5eXlKigo0KBBgxQfH69t27bJ4XBccWbDMDRz5sw6q1SXctnzBwAAmh1WlgBviY+vLSn18Z//SMuWSf7+f60sPfrotWWpp8TERI/trl27au/evRo+fLh7zM/PT3FxcSotLfXYd+PGjVq8eLH27t2rM2fOuMf/XkycTqc2bNig3r17KyEhQSEhIVq3bt1VFSUAAIBrRVkCvCUo6KpXcyTVXnq3bNlf9yhduGcpOtqr9yxdEB4e7rEdEBCgoKCgOoUnICBAp06dcm/n5+crOTlZSUlJSktLU1RUlPz8/LR06VLl5OTUeZ+4uDj1799fBQUFeuKJJ9S2bdvG+UAAAACXQFkCfJnVwxwu/LZ6Sp4PW7t2rWJjY7V+/XoZhuEer6qqstx/xYoVKigo0B133KGsrCyNGzdOd955p7fiAgAA8DQ8wKe5XNZPvbvwlDyXy55c9WAYhgICAjyK0tGjR+s8DU+S9u/fryeffFKTJk3S119/rR49emjcuHH6888/vRkZAAA0c5QlwJdlZFx65ei555rUF9ImJibq0KFDmjFjhgoLC7V69Wrdc889dS6vO3PmjJKTkxUbG6slS5YoICBAGzZsUHl5uaZMmWJTegAA0BxRlgB4xZQpU7Rw4UJ9/vnnGjFihF5//XWlp6dr4sSJHvtNnz5dv/76qzZu3KhWrVpJkjp27KgVK1Zo06ZNWrx4sQ3pAQBAc8Sjw4EGxKOnmzbOHwAAuBgrSwAAAABggbIEAAAAABYoSwAAAABggbIEAAAAABYoSwAAAABggbIEAAAAABYoSwAAAABggbIEAAAAABYoSwAAAABggbIEAAAAABYoSwC84siRIzIMQ6tWrbrmub766isZhnFFPwAAAPXlZ3cAAM1D27ZtVVRUpE6dOl3zXAkJCSoqKvIYGz16tDp16qQ333zzmucHAACQKEsAvMTpdKpv374NMldISEiduZxOp8LCwhrsPQAAALgMD2gi/vtfafv22t92yMjIkGEY2rdvn5KSkhQaGqrw8HClpqbq/PnzOnTokIYNG6bg4GDFxMRo0aJFHsdbXYZ3Yc4DBw5owoQJCg0NVWRkpKZOnaqTJ096+RMCAAB4oiwBTUB2tnTTTdJ999X+zs62L0tycrJ69uypvLw8Pfroo3rnnXc0d+5cPfDAAxo5cqQ+/vhj3XfffVqwYIHy8/OvaM4HH3xQXbp0UV5entLT05Wbm6u5c+c28icBAAC4PC7DA7ykokI6ePDqjzt2THrsMammpna7pqZ2OzpaioysX5b4eCkoqH7HPvbYY0pNTZUkDRkyRF9++aWysrKUn5+v0aNHS5IGDhyozZs368MPP9SYMWP+cc5p06YpLS3NPefhw4eVk5Oj7OxsHtIAAABsQ1kCvOTgQal374aZq6ZGGjGi/sfv3i0lJNTv2MTERI/trl27au/evRo+fLh7zM/PT3FxcSotLb2iOe+//36P7R49eqiyslJlZWWKrG8jBAAAuEaUJcBL4uNrS8rVOnZMSkz8a2VJklq0kDZvvraVpfoKDw/32A4ICFBQUJACAwPrjJ86deqK5oyIiPDYdjqdkqSzZ8/WPygAAMA1oiwBXhIUVP/VnOXLpccfl1wuyeGQ/v1v6aKFHAAAADQCyhLQBEybJv3rX9Lhw1JcnNS+vd2JAAAArn+UJaCJaN+ekgQAAOBNPDocAAAAACwYpmmadocArheVlZUqKSlRbGxsnQcewPdx/gAAwMVYWQIAAAAAC5QlAAAAALBAWQIAAAAAC5QlAAAAALBAWQIAAAAAC5QloBHwkMmmifMGAAAuRlkCGpC/v78kqaKiwuYkqI8L5+3CeQQAAM2bn90BgOuJw+FQWFiYysrKJElBQUEyDMPmVPgnpmmqoqJCZWVlCgsLk8PhsDsSAADwAXwpLdDATNPU0aNHVV5ebncUXKWwsDBFRUVRcAEAgCTKEtBoXC6Xqqur7Y6BK+Tv78+KEgAA8EBZAgAAAAALPOABAAAAACxQlgAAAADAAmUJAAAAACxQlgAAAADAAmUJAAAAACxQlgAAAADAAmUJAAAAACz8H/9BBYqUnS68AAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(10,5))\n", "\n", "data.plot(x='month', y='Avg_max_T_degF', c='r', linestyle='-', linewidth=1, marker='x', ax=ax, label='max T')\n", "data.plot(x='month', y='Avg_min_T_degF', c='b', linestyle='-', linewidth=1, marker='.', ax=ax, label='min T')\n", "\n", "ax.set_title('Average temperature at Cedar Lake, Washington', fontsize=15)\n", "ax.set_xlabel('Month', fontsize=12)\n", "ax.set_ylabel('Temperature (deg F)', fontsize=12);\n", "plt.legend(loc=(0.2,-0.36), fontsize=12);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**How is temperature related to the water balance?** What do you know about the total elevation range that drains to Cedar Lake? What do you know about how temperature changes with elevation? Do you think there are times when precipitation may fall as snow or as a mix of rain and snow? How do you think this will affect the monthly water balance? " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "---\n", "## **Basic Water Balance**\n", "\n", "Over a watershed, we can define the monthly change in storage as the sum of water input minus water output.\n", "\n", "$$ \\Delta S = P - ET - Q $$\n", "\n", "where $\\Delta$S is change in storage, $P$ is area-average precipitation, $ET$ is area-average evapotranspiration, and $Q$ is basin runoff, normalized by basin area.\n", "\n", "**Let's calculate the change in storage each month, in mm.**\n", "First, convert everything to the same units. Then combine.\n", "\n", "1 inch = 25.4 mm\n", "\n", "1 foot = 0.3048 m\n", "\n", "Basin area = 217 km$^2$" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYEAAAG2CAYAAACZJjoxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAACO3UlEQVR4nO3dd1hTZxsH4F9YYe8tMtwLR904UNwD96jbav0cdWurtiquiqOOVutobR11YF111Q3iniDiXqioICKyd/J8f7xNIAJKIJO893Wdy3DOyTlPYnKenHcKiIjAcRzH6SQ9dQfAcRzHqQ9PAhzHcTqMJwGO4zgdxpMAx3GcDuNJgOM4TofxJMBxHKfDeBLgOI7TYQbqDkATiMVivHnzBhYWFhAIBOoOh+M4rtSICCkpKXB1dYWeXtG/93kSAPDmzRuUL19e3WFwHMcpXHR0NNzc3IrczpMAAAsLCwDszbK0tFRzNBzHcaWXnJyM8uXLS69vReFJAJAWAVlaWvIkwHFcmfK5Im5eMcxxHKfDeBLgOI7TYTwJcBzH6TCeBDiO43QYTwIcx3E6jCcBjuM4HcaTAMdxnA7jSYDjOE6H8STAcRynw3gS4DiO03A5OTlKOzYfNoLjOE5DvHv3Dvfv38f9+/dx79496WMPDw9cuHBBKefkSYDjOE6FiAivXr0qcKG/f/8+4uPjC31OWloaiEgpQ93zJMBxHKcEIpEIz549K/Rin5qaWuTzPD09UaNGDVSvXh3Vq1eXPlbWXCc8CXAcx5VCVlYWHj16VOBi//DhQ2RnZxf6HAMDA1SqVKnAxb5q1aowNTVVafxqTQLr16/H+vXr8fz5cwBAzZo1MXfuXHTq1AkAMHz4cGzdulXmOY0bN8aVK1ekf2dlZWH69OnYtWsXMjIy0KZNG6xbt+6TkyhwHMeV1oYNG7By5Uo8ffoUYrG40H2MjY1RrVq1Ahf7ihUrwsjISMURF05ARKSukx8+fBj6+vqoVKkSAGDr1q1Yvnw5wsPDUbNmTQwfPhxv377F5s2bpc8xMjKCra2t9O+xY8fi8OHD2LJlC+zs7DBt2jQkJCTg5s2b0NfXL1YcycnJsLKyQlJSEp9PgOO4z8rJyYGdnR1SUlIAAFZWVgWKb6pXrw4PD49iX4cUrbjXNbXeCfj7+8v8/eOPP2L9+vW4cuUKatasCQAQCoVwdnYu9PlJSUn4448/8Ndff6Ft27YAgO3bt6N8+fI4ffo0OnToUOjzsrKykJWVJf07OTlZES+H4zgdcfnyZaSkpMDe3h4RERFwcXHR2vnJNaafgEgkQlBQENLS0tC0aVPp+rNnz8LR0RFVqlTBqFGjEBcXJ9128+ZN5OTkoH379tJ1rq6uqFWrFi5dulTkuQIDA2FlZSVd+PzCHMfJ4/jx4wCA9u3bw9XVVWsTAKABSSAyMhLm5uYQCoUYM2YMDhw4gBo1agAAOnXqhB07diA4OBgrVqzA9evX4efnJ/0VHxsbCyMjI9jY2Mgc08nJCbGxsUWec9asWUhKSpIu0dHRynuBHMeVOSdOnAAAdOzYUc2RlJ7aWwdVrVoVt27dQmJiIvbt24dhw4YhNDQUNWrUQP/+/aX71apVCw0aNICHhweOHj2KXr16FXnMz7WnFQqFEAqFCn0dHMfphrdv3yIsLAwAZEohtJXa7wSMjIxQqVIlNGjQAIGBgahTpw5+/vnnQvd1cXGBh4cHHj9+DABwdnZGdnY2Pnz4ILNfXFwcnJyclB47x3G65+TJkwCAevXqlYnrjNqTwMeISKbSNr/3798jOjoaLi4uAID69evD0NAQp06dku4TExODO3fuwMfHR6lxfvjwAatXr8b48eOVeh6O4zSLpD6gLBQFAWouDvr+++/RqVMnlC9fHikpKQgKCsLZs2dx/PhxpKamYt68eejduzdcXFzw/PlzfP/997C3t0fPnj0BsGZZI0eOxLRp02BnZwdbW1tMnz4d3t7e0tZCypKcnIwpU6ZAIBBgxowZvHKZ43SAWCyW3gmUlSSg1juBt2/fYsiQIahatSratGmDq1ev4vjx42jXrh309fURGRmJ7t27o0qVKhg2bBiqVKmCy5cvw8LCQnqMVatWoUePHujXrx+aNWsGU1NTaf8DZfLw8ECLFi1ARNi1a5dSz8VxnGYICwtDfHw8LCwsZFoxajO1dhbTFCXtLPbbb79h9OjRqF27NiIiIpQYIcdxmmDRokWYM2cOevTogQMHDqg7nE8q7nVN4+oEtEmfPn1gaGiI27dv4/bt2+oOh+M4JStr9QEATwKlYmtriy5dugAAduzYoeZoOI5TpsTEROm4ZUWNRqCNeBIopcGDBwMAdu7cWeQgUhzHab8zZ85AJBKhatWq8PT0VHc4CsOTQCl16dIFVlZWePXqFc6dO6fucDiOU5KyWBQElDAJREdH4/z58zhx4gTCwsKKbNevC4yNjdGnTx8AbPA6juPKHiLiSeDFixeYNWsWPD094enpCV9fX3Tq1AkNGjSAlZUV2rVrhz179uhkkYikSGjv3r3IzMxUczQcxyna/fv38erVKxgbG8PX11fd4ShUsZLApEmT4O3tjcePH2PBggW4e/cukpKSkJ2djdjYWPz7779o3rw55syZg9q1a+P69evKjlujtGzZEm5ubkhKSsLRo0fVHQ7HcQomuQto2bIlTExM1ByNYhUrCRgZGeHp06fYu3cvhg4dimrVqsHCwgIGBgZwdHSEn58fAgIC8ODBAyxbtgwvXrxQdtwaRU9PDwMHDgTAWwlxXFlUVouCAN5ZDIBiZhaLjIxE7dq1YWRkhJiYGJnZzziO017p6emwtbVFVlYW7t27h+rVq6s7pGLhncVUzNvbG97e3sjOzsbevXvVHQ7HcQoSGhqKrKwsuLu7o1q1auoOR+HkTgKS8X5cXV1hYGAAfX19mUWXSSqIeZEQx5UdkqKgDh06aPUMYkWRexTR4cOH4+XLl5gzZ45Wz6upDAMGDMDMmTNx7tw5vHjxAh4eHuoOieO4UirL9QFACZLAhQsXcP78edStW1cJ4Wi38uXLw9fXF2fPnsXOnTsxa9YsdYfEcVwpREVF4dGjR9DX10ebNm3UHY5SyF0cVL58efC65KJJioS2b9/O3yeO03KSuYR9fHxgZWWl5miUQ+4ksHr1asycORPPnz9XQjjar3fv3jAyMsK9e/f48NIcp+Xy1weUVcUqDrKxsZEp+09LS0PFihVhamoKQ0NDmX0TEhIUG6GWsba2hr+/P/bt24ft27fzYjOO01LZ2dk4c+YMgLJbHwAUMwmsXr1ayWGULYMHD8a+ffuwa9cuLF26VOdbTXGcNrp8+TJSU1Ph4OCAevXqqTscpSlWEhg2bJiy4yhTOnXqBGtra7x58wZnz54tsxVKHFeW5S8K0tMru12q5H5l+vr6iIuLK7D+/fv3/Bfvf4RCIfr16weA9xngOG2lC/UBQAmSQFEtXrKysmBkZFTqgMqK/COLZmRkqDkajuPkERsbi1u3bgEA2rdvr95glKzY/QR++eUXAIBAIMCmTZtgbm4u3SYSiXDu3Lky2aW6pJo1awZ3d3e8fPkShw8flt4ZcByn+U6ePAkAqF+/PhwdHdUcjXIVOwmsWrUKALsT2LBhg0zRj5GRETw9PbFhwwbFR6il9PT0MGjQIAQGBmLHjh08CXCcFinrvYTzk3sU0datW2P//v2wsbFRVkwqp4hRRAtz79491KxZEwYGBoiJiYG9vb3Cjs1xnHKIRCI4OTnh/fv3OHfuHFq0aKHukEpEaaOIhoSElKkEoEw1atRA3bp1kZubiz179qg7HI7jiuHmzZt4//49LC0t0aRJE3WHo3Ryjx0kEomwZcsWnDlzBnFxcQWmkwwODlZYcGXB4MGDcevWLezYsQNjx45Vdzgcx32GZKiItm3bFugMWxbJfScwadIkTJo0CSKRCLVq1UKdOnVkFk7WgAEDIBAIcPHiRTx79kzd4XAc9xm6VB8AlOBOICgoCH///Tc6d+6sjHjKHFdXV/j5+eHMmTPYuXMnZs+ere6QOI4rwocPH3DlyhUAGtQ/QCwGsrMBY2OlHF7uOwEjIyNUqlRJGbGUWfknm+Eji3Kc5jp9+jTEYjGqV68Od3d3dYcDREYCzZsDM2Yo7RRyJ4Fp06bh559/5hczOfTq1QvGxsZ48OABwsLC1B0Ox3FFkNQHqL0oKD0dmDkT+OIL4PJlYOtWQEmDc8qdBC5cuIAdO3agYsWK8Pf3R69evWQWeaxfvx61a9eGpaUlLC0t0bRpUxw7dky6nYgwb948uLq6wsTEBK1atcLdu3dljpGVlYUJEybA3t4eZmZm6NatG169eiXvy1IqS0tLdOvWDQAfRoLjNBURaUZ9wLFjQM2awNKlQG4u0LMncOcOYGurlNPJnQSsra3Rs2dP+Pr6wt7eHlZWVjKLPNzc3LBkyRLcuHEDN27cgJ+fH7p37y690C9btgwrV67E2rVrcf36dTg7O6Ndu3ZISUmRHmPy5Mk4cOAAgoKCcOHCBaSmpqJr164QiUTyvjSlkhQJ7dq1C7m5uWqOhuO4j929exevX7+GsbGxevoGxMQA/fsDnTsDz58D5csDBw8C+/cDbm7KOy9pGBsbG9q0aROJxWJydnamJUuWSLdlZmaSlZUVbdiwgYiIEhMTydDQkIKCgqT7vH79mvT09Oj48ePFPmdSUhIBoKSkJMW9kI9kZWWRnZ0dAaATJ04o7Twcx5XM8uXLCQB17NhRtScWiYjWrSOytCQCiPT0iKZOJUpJKdVhi3tdK/H4qO/evcOFCxdw8eJFvHv3rtTJSCQSISgoCGlpaWjatCmioqIQGxsrM3iTUCiEr68vLl26BIB16sjJyZHZx9XVFbVq1ZLuU5isrCwkJyfLLMpmZGTERxblOA2mlvqAiAjAxwcYNw5ITgYaNgRu3ABWrADyjc+mTHIngbS0NIwYMQIuLi5o2bIlWrRoAVdXV4wcORLp6elyBxAZGQlzc3MIhUKMGTMGBw4cQI0aNRAbGwsAcHJyktnfyclJui02NhZGRkYFejDn36cwgYGBMkVY5cuXlzvukpAUCe3fvx9paWkqOSfHcZ+XlpaGc+fOAVBR09C0NOC774D69YGrVwELC2DNGlYJrOIJbOROAlOnTkVoaCgOHz6MxMREJCYm4uDBgwgNDcW0adPkDqBq1aq4desWrly5grFjx2LYsGG4d++edHv+aS0BVnnz8bqPfW6fWbNmISkpSbpER0fLHXdJNG3aFF5eXkhNTcWhQ4dUck6O4z7v7NmzyM7OhoeHB6pWrarckx09yip+ly8HRCKgd2/g/n1g/HhADXOyyJ0E9u3bhz/++AOdOnWSturp3Lkzfv/9d+zdu1fuACT9Dho0aIDAwEDUqVMHP//8M5ydnQGgwC/6uLg46d2Bs7MzsrOz8eHDhyL3KYxQKJTGLllUQSAQYNCgQQB4kRDHaZL8RUGf+5FZYm/eAH37Al27Ai9eAO7uwOHDwN69QLlyyjlnMcidBNLT0wu9wDo6OpaoOOhjRISsrCx4eXnB2dkZp06dkm7Lzs5GaGgofHx8ALCxvg0NDWX2iYmJwZ07d6T7aBpJEjh+/LhC6lI4jis9pTYNFYmAtWuBatXYBV9fH5g+Hbh3jyUEdZO3xtnPz4/69u1LGRkZ0nXp6enUt29fatOmjVzHmjVrFp07d46ioqLo9u3b9P3335Oenh6dPHmSiIiWLFlCVlZWtH//foqMjKQBAwaQi4sLJScnS48xZswYcnNzo9OnT1NYWBj5+flRnTp1KDc3t9hxqKJ1UH4NGjQgALRmzRqVnI/juKI9efKEAJCBgYHirwHh4UQNG7JWPwBR48ZEt24p9hxFKO51Te4kEBkZSeXKlSM7Ozvy8/OjNm3akJ2dHZUrV47u3Lkj17FGjBhBHh4eZGRkRA4ODtSmTRtpAiAiEovFFBAQQM7OziQUCqlly5YUGRkpc4yMjAwaP3482drakomJCXXt2pVevnwpVxyqTgKrVq0iANSkSROVnI/juKL9+uuvBIBatmypuIOmpLBmnvr67OJvaUn0669Ecvw4La3iXtfknlQGADIyMrB9+3Y8ePAARIQaNWpg0KBBMDExUeRNisooa1KZosTGxqJcuXIQi8V4/PgxH4uJ49Soe/fuOHToEBYvXoxZs2aV/oCHDwPffANIGpz06wesXg24uJT+2HIo7nVN7lFEAcDExASjRo0qcXC6TtLz+cSJE9i5cyfmzp2r7pA4TidlZ2fjzJkzABRQH/DqFTBxInDgAPvb0xNYtw7o1Kl0x1WyEiWB169f4+LFi4VOKjNx4kSFBFbWDRo0CCdOnMD27dsxZ84c5bVI4DiuSBcvXkRaWhocHR1LPh+KSAT8+ivwww9AaipgYABMmwbMnQuYmio2YCWQOwls3rwZY8aMgZGREezs7GQuXgKBgCeBYurZsydMTU3x+PFjXL9+HY0aNVJ3SByncyStgjp06AA9vRIMoHDzJjB6NPsXAJo2BTZuBLy9FRilcsn9qufOnYu5c+ciKSkJz58/R1RUlHThM2cVn7m5OXr06AGA9xngOHUp1VARK1YAjRqxBGBlBWzYAFy4oFUJAChhP4Evv/yyZFmTkyHpM7Br1y7k5OSoORqO0y1v3rxBREQEBAIB2rVrJ++TgVmz2KxfX34JPHjA7gi08Lood8QjR47Enj17lBGLzmnXrh0cHBzw7t07nD59Wt3hcJxOOXnyJADW6dTBwUG+J69ZA+TkAM2aAbt2Af+NcKCN5K4TCAwMRNeuXXH8+HF4e3vD0NBQZvvKlSsVFlxZZ2hoiP79+2Pt2rXYsWMHOml4KwKOK0tK3Es4JYUV/QDAt98qOCrVkzsJLF68GCdOnJAOsvRxxTAnn8GDB2Pt2rU4cOAAUlNTYa6i4WM5TpeJRCLpcDNyJ4E//wQSE4HKlQF/f8UHp2JyJ4GVK1fizz//xPDhw5UQju5p1KgRKlWqhCdPnuCff/6RDjfNcZzy3LhxAwkJCbCyskLjxo2L/8TcXGDVKvZ42jStrAP4mNyvQCgUolmzZsqIRSfxkUU5TvUkRUFt27aFgYEcv4X37mUjgDo4AEOHKik61ZI7CUyaNAlr1qxRRiw6S5IETp48ibdv36o5Go4r+0pUH0DE5gAA2Nj/WjpMzsfkLg66du0agoODceTIEdSsWbNAxfD+/fsVFpyuqFy5Mho3boyrV69i9+7dvMMdxylRQkICrl27BkDOWcRCQ4GwMMDYmE0HWUbIfSdgbW2NXr16wdfXF/b29jLTNFpZWSkjRp0guRvYvn27miPhuLLt9OnTEIvFqFmzpnxTy/70E/v3q68Ae3vlBKcGJRo2glO8/v37Y8qUKbh+/ToePXqEKlWqqDskjiuT8g8VUWz37rFpIQUCYMoUJUWmHtpftV1GODo6Sj+UvIKY45SDiEo2VMSKFezfHj1Y09AyhCcBDZK/SKgE0zxwHPcZd+7cwZs3b2BiYoIWLVoU70kxMYCkmLYMdA77GE8CJZWbC2zeDHz/vcIO2b17d5iZmeHZs2e4cuWKwo7LcRwjKQpq3bo1jI2Ni/ektWuB7GzAx4eNElrG8CRQUpGRwIgRwNKlbPAoBTAzM0OvXr0A8CIhjlMGuesDUlOB9evZ4+nTlRSVesmVBHJyctC6dWs8evRIWfFoj3r1gO7d2SiCCxYo7LCSIqGgoCA+sijHKVBqaiouXLgAQI76gM2bgQ8fgEqVgG7dlBid+siVBAwNDXHnzh0+RpDE/Pns36Ag4O5dhRyyTZs2cHJywvv376UVWBzHld7Zs2eRnZ0NLy8vVC5O5W5uLiAZEHPqVEBfX7kBqoncxUFDhw7FH3/8oYxYtE+dOkCfPqwnoSQhlJKBgQEGDBgAgBcJcZwi5S8KKtYP2f37gefPWZ+AYcOUG5wayd1PIDs7G5s2bcKpU6fQoEEDmJmZyWzXuaGkAwKAffuAPXuAiAiWGEpp0KBBWL16NQ4ePIiUlBRYWFgoINBPS05ORkZGBpycnJR+Lo5TB7mGiiDK6xz2zTdaMVdwSQlIzraIrVu3LvpgAgGCg4NLHZSqJScnw8rKCklJSbC0tJT/AAMGsCKhHj2AAwdKHQ8RoXr16nj48CG2bt2KoUoYqCo+Ph4XLlzAuXPncO7cOYSHh0NfXx/nz5+Xb1RFjtMCT548QeXKlWFgYICEhITP/7A6dw7w9WVDRLx4ATg6qiZQBSr2dY04SkpKIgCUlJRUsgPcv0+kp0cEEN24oZCYFi5cSACoXbt2Cjneq1evaNeuXTR27FiqWbMmASh08fX1JbFYrJBzcpymWLt2LQGgVq1aFe8J/v7s+zx6tHIDU6LiXtdKnAQeP35Mx48fp/T0dCIirb5wlDoJEBENHsw+NF27KiSmp0+fEgDS09OjN2/eyPVcsVhMT58+pc2bN9NXX31FFStWLPSCX7NmTRo7dizt2rWLrl27RkKhkADQ8ePHFfIaOE5TdO3alQBQYGDg53e+d499lwUCoocPlR+ckigtCcTHx5Ofnx8JBALS09Ojp0+fEhHRiBEjaOrUqSWLVs0UkgQePSLS12cfnitXFBKXj48PAaCVK1d+cj+xWEx3796l9evX04ABA6hcuXIFLvh6enpUv359mjJlCh04cIDevXtX4DhTp04lAPTFF1+QSCRSyGvgOHXLzMwkU1NTAkDh4eGff8LXX7PvcY8eSo9NmZSWBIYMGUIdOnSg6OhoMjc3lyaBEydOUI0aNUoWrZopJAkQEX31FfvwdOigkLjWrVsnvSjnl5ubSzdv3qRVq1ZRz549yd7evsBF39DQkJo1a0azZs2iY8eOFeu1xcXFkbm5OQGgv//+WyGvgePU7cyZMwSAnJ2dP19iERNDZGTEvscXLqgmQCVRWhJwcnKiW7duERHJJIFnz56RmZlZCUJVP4UlgadPiQwMFPYBevfuHRkYGBAA2rVrFwUGBlKnTp3IwsKiwEXfxMSE2rRpQ/Pnz6eQkBBKS0sr0TkDAgIIAFWtWpVycnJK/Ro4Tt2+/fZbAkDDhg37/M6zZ7Pvb5MmRFpcxE2kxCRgbm5Ojx49kj6WJIFr166Rra1tCUJVP4UlASKiUaPYh6hNm9Ifi4j8/f0LLc+3tLSkzp0705IlS+jSpUuUlZWlkPMlJSWRnZ0dAaBNmzYp5Jgcp07e3t4EgHbu3PnpHVNTiWxt2fd3717VBKdESksCnTt3ptmzZxMRSwLPnj0jkUhEffv2pd69e5csWjVTaBJ48YLI0JB9kM6eLfXhTp8+TXp6emRvb0+9evWi1atXU1hYGOXm5pY+1iKsWLGCAJCbmxtlZGQo7Twcp2yvXr0iACQQCAqtB5OxZg373lasSKTE75eqKC0J3L17lxwcHKhjx45kZGREffr0oerVq5OTkxM9efJErmMtXryYGjRoQObm5uTg4EDdu3enBw8eyOwzbNiwAr+CGzduLLNPZmYmjR8/nuzs7MjU1JT8/f0pOjq62HEoNAkQEY0bxz5MLVsq5JYyIyNDpa2vMjIyyM3NrViV0hynyf78808CQI0aNfr0jrm5RF5e7Hv766+qCU7Jintdk3vYiBo1auD27dto1KgR2rVrh7S0NPTq1Qvh4eGoWLGiXMcKDQ3FN998gytXruDUqVPIzc1F+/btkZaWJrNfx44dERMTI13+/fdfme2TJ0/GgQMHEBQUhAsXLiA1NRVdu3aFSCSS9+UpxqxZgFDIOpwooPOcsbGxSsdrMjY2RkBAAABg8eLFSElJUdm5OU6Rit1L+MABICoKsLMDhg9XfmCaREVJqVji4uIIAIWGhkrXDRs2jLp3717kcxITE8nQ0JCCgoKk616/fk16enrFbu+u8DsBIqKJE9mvCh8fraxgysnJoSpVqhAAmjdvnrrD4Ti55eTkkI2NDQGgixcvFr2jWEzUqBH7vs6Zo7oAlUxpdwIA8OHDB/z0008YOXIkvv76a6xYsQIJCQmlTkhJSUkAAFtbW5n1Z8+ehaOjI6pUqYJRo0YhLi5Ouu3mzZvIyclB+/btpetcXV1Rq1YtXLp0qdDzZGVlITk5WWZRuFmzABMT4NIlQAtHAzUwMMDChQsBACtWrEB8fLyaI+I4+Vy/fh0fPnyAtbU1GjVqVPSOFy8C166xu/dvvlFdgBpC7iQQGhoKLy8v/PLLL/jw4QMSEhLwyy+/wMvLC6GhoSUOhIgwdepUNG/eHLVq1ZKu79SpE3bs2IHg4GCsWLEC169fh5+fH7KysgAAsbGxMDIygo2NjczxnJycEBsbW+i5AgMDYWVlJV3Kly9f4riL5OwMjBvHHs+dywak0jJ9+vRBvXr1kJKSgiVLlqg7HI6Ti2Qo9nbt2sHA4BNjZUoGihs2DNDFARTlvcWoWbMmjRo1SqZ1Sm5uLv3vf/+jmjVryn3LIjFu3Djy8PD4bIXumzdvyNDQkPbt20dERDt27CAjI6MC+7Vt25ZGFzHuR2ZmJiUlJUmX6OhoxRcHERG9fUtkaspuMw8fVuyxVeTYsWMEgIRCoVyV7Rynbo0bNyYA9McffxS904MHbHgIgD0uQ5RWHPT06VNMmzYN+vkmWNDX18fUqVPx9OnTEiWiCRMm4NChQwgJCYGbm9sn93VxcYGHhwceP34MAHB2dkZ2djY+fPggs19cXFyRwyILhUJYWlrKLErh6AhMmMAea+ndQIcOHdCyZUtkZWVhgQJnUOM4ZXr//j2uXbsGADJFxQWsXMm+l926AVWrqig6zSJ3Evjiiy9w//79Auvv37+PunXrynUsIsL48eOxf/9+BAcHw8vL67PPef/+PaKjo+Hi4gIAqF+/PgwNDXHq1CnpPjExMbhz5w58fHzkikcppk8HzM2B8HDgn3/UHY3cBAIBAgMDAQB//vknn1qU0wqnTp0CEaFWrVpF/7CMiwO2bmWPy+j8wcUi7y1GUFAQubu70/Lly+n8+fN0/vx5Wr58OXl6elJQUBBFRERIl88ZO3YsWVlZ0dmzZykmJka6SEYmTUlJoWnTptGlS5coKiqKQkJCqGnTplSuXDlKTk6WHmfMmDHk5uZGp0+fprCwMPLz86M6deoUu0OVUloH5Sfpiu7tTaSlA7NJRmHs37+/ukPhuM8aPnw4AaDp06cXvdPcuex72aiRVrbg+xyldRYTCASfXPT09KT/fg6KGNN+8+bNRESUnp5O7du3JwcHBzI0NCR3d3caNmwYvXz5UuY4GRkZNH78eLK1tSUTExPq2rVrgX0+RelJICGByNKSfeC0dGC2iIgIEggEBIDCwsLUHQ7HFUksFpOzszMBoNOnTxe+U1oakZ0d+07u2aPaAFWkuNc1uWcWe/HiRbH39fDwkOfQalPqmcWKY/58YN48oEYN4PZtrZy0etCgQdi5cyc6depUoMMex2mKiIgI1K1bF6ampkhISIBQKCy407p1rDmolxfw+LFWfh8/p7jXNbnnGNaWC7vGmTwZ+Pln4N49YPduYOBAdUckt/nz5+Pvv//GsWPHcP78ebRo0ULdIXFcAZJewq1atSo8AYhErEIYAKZOLZMJQB4l6izGlYCVVV7l0/z5QG6ueuMpgUqVKmHkyJEAgFmzZkHOm0iOU4mjR48CADp37lz4DgcPAk+fAjY2wFdfqTAyzcSTgCpNmMDGJnn0CNi5U93RlMicOXNgbGyMixcv8iIhTuN8+PBBOlJAkUlA0jls3DjAzExFkWkungRUycIC+O479njBAiAnR73xlEC5cuUw4b++Dz/88APEYrGaI+K4PCdPnoRIJEL16tULb3J+6RJw+TJgZJTXh0fH8SSgat98wzqRPX0KbNum7mhKZMaMGbC0tERERAR2796t7nA4Tkpyd9qlS5fCd5DcBQwdqptDRBSiVEkgNTVV+QOxlTVmZsDMmezxwoVAdrZ64ykBOzs7fPvttwBY8VCOFt7RcGWPWCzGsWPHABSRBB49yuuwOXWq6gLTcHIngaioKHTp0gVmZmawsrKCjY0NbGxsYG1tXWAQN64IY8awAeZevAA2b1Z3NCUyefJkODo64unTp/jzzz/VHQ7H4fr163j37h0sLS3RrFmzgjusWsWGiOjaFaheXfUBaii5k8CgQYPw4cMH/Pnnnzhz5gyCg4MRHByMkJAQBCtgAhWdYGICfP89e7xoEfDfiKjaxNzcHD/88AMAYMGCBcjIyFBzRJyukxQFtW/fHoaGhrIb370Dtmxhj3V5iIjCyNsLzczMrMAUkNpO6T2GC5ORQVSuHOuxuGaN6s6rQJmZmeTh4UEAaPny5eoOh9Nx9evXJwD0559/Ftw4bx77rjVsWCaHiCiM0kYRbdiwIaKjoxWdi3SPsTHw3y9pLF4MaOEvaaFQiHnz5gFgczRIJgXiOFWLjY3FzZs3AbA5SGRkZABr17LH06cDKpyqVRvInQQ2bdqEpUuXYuvWrbh58yZu374ts3ByGDkScHcHYmKAjRvVHU2JDBkyBNWrV0dCQgJWrFih7nA4HSWpEG7QoAGcnZ1lN27dCsTHA56eQK9eqg9Ow8mdBN69e4enT5/iq6++QsOGDVG3bl3Uq1dP+i8nByMjYM4c9jgwEEhLU288JaCvr49FixYBAFauXCkz9SfHqUqRvYTzDxExZQrwqRnGdJTcSWDEiBGoV68eLl++jGfPniEqKkrmX05Ow4axQazi4tigVlqoZ8+eaNiwIdLS0rB48WJ1h8PpmJycHJw8eRJAIU1DDx9mA8TZ2AAjRqghOi0gb2WDqakpPX78uMSVFZpILRXD+W3ezCqt7OyI8s2ToE1OnTpFAMjIyIieP3+u7nA4HRIcHEwAyMHBgUQfz9fRrBn7bn3/vXqCUyOlVQz7+fkhIiJC0blItw0eDFSuDLx/n1eBpWXatm0LPz8/ZGdnY/78+eoOh9MhkqKgTp06QU8v3yXt8mXg4kVW7Dp+vJqi03xyF5D5+/tjypQpiIyMhLe3d4H2uN26dVNYcDrDwAAICGDJYPlyNrSEsuY1UKIff/wRTZs2xdatW/Htt9+iOu+Qw6lAkUNFSBoqDB4M/DcdLVeQ3JPKyGTajw8mEEAkEpU6KFVTyaQynyMSAbVqAQ8esKGm585VTxyl1KNHDxw8eBC9e/fG3r171R0OV8ZFRUWhQoUK0NfXR3x8PKytrdmGJ0+AKlVYD+G7d9lkTjqmuNc1uYuDxGJxkYs2JgCNoa/PZh4DWGuGDx9Uc96sLODQIfZryd+fFUmVwqJFiyAQCLBv3z7cuHFDQUFyXOEkdwHNmjXLSwBA3hARXbroZAKQBx9FVJP07cvuBpKS2IdYWbKzgWPHgOHD2UiK3bsDO3YAR46wIa5LoVatWhg8eDAA4HvJ0BgcpySFNg2Nj88bk4sPEfFZchcHAcC1a9dw9uxZxMXFFRhPfqWkTa4W0YjiIIn9+4HevdncA1FRbBIaRcjNBUJC2NSW+/fL3mm4ugKtWrGJbgwN2WiLnp4lPlVUVBSqVq2KnJwcBAcHo3Xr1qUOn+M+lp6eDjs7O2RmZiIyMhK1atViG5YvZ/N21K8PXL+usz2ElVYctHjxYjRp0gSbN2/GjRs3EB4eLl1u3bpVmpg5AOjRA6hbF0hJyRv7vKREInbhHzuWXejbtwf++IMlACcnVgF97hwQHc3uBNq2ZRPdlLI+wsvLC//73/8AsLuBEvzO4LjPCgkJQWZmJtzd3VGzZs28DQcOsH9HjtTZBCAXedueOjo60ubNm0vQalVzqb2fwMcOHmRtm83MiN6+le+5IhHR+fNE48cTOTuz40gWe3ui0aOJzpwhys0t+Nzr19l+AgHR7dulegkxMTFkampKAOjgwYOlOhbHFWbs2LEEgMaMGZO3MiaGfX4Botev1RecBlBaPwE9Pb3Cx+rmFMffH2jQgA0jsXz55/cnAq5cYd3i3d2BFi1Yf4PY2LyekidOsDGKNmwA/PxYRfTHGjRg9RJEeUNdl5CzszMmTZoEgE1DyRsNcIpERIU3DT18mH1+GzVid7/c58mbXZYuXUqTJk0qYW7STBp3J0BE9O+/7NeMiQn7dfMxsZj9cp8+ncjdXfYXv6Ul0dChREePEmVlyXfehw+J9PXZcc6fL9VLSEhIIGtrawJAf/31V6mOxXH53blzhwCQUCiktLS0vA2dO7PP7uLF6gtOQxT3uiZ3EhCJRNSxY0eqUKECde3alXr27CmzaCONTAJiMVGTJuwDLUm6YjHRrVtEs2YRVawoe+E3NycaOJAVJWVmlu7c//sfO2azZqUeez0wMJAAkJeXF2XJm5A4rghLly4lANSxY8e8lcnJREZG7LN79676gtMQSisOmjBhAkJCQlClShXY2dnByspKZuEURCDIa665YQMwaxabEq9uXTbi6NOnbIayfv2AvXvZAHQ7dgDdugFCYenOPXcum+/g4kXgvyZ4JTVx4kQ4OzsjKioKmzZtKl1cHPefQouCTpxgzZ8rV+bTR8pD3uxibm5OR44cKXF20kQaeSdAxH6FN28u+4tfKCTq2ZMoKIgoNVV5554xg52vVq3CK5Hl8OuvvxIAcnZ2plRlxszphA8fPpC+vj4BoKdPn+ZtGDSIfWanT1dfcBpEaXcCtra2qFixoqJzEVcYgQBYvZp1f/f3B7ZvZ3Ol7t8P9O8PmJkp79wzZgDW1sCdO6z/QCl8/fXX8PLyQmxsLNasWaOY+DiddfLkSYhEIlSrVg0VKlRgK3Ny8u5ae/RQW2zaSO4kMG/ePAQEBCA9PV0Z8XAfq18fePiQDe0waBDrRKYKNjbAzJns8Zw5bHiJEjIyMsKC/4q2li5dig+qGhKDK5MKLQo6dw5ITAQcHYEmTdQTmJaSOwn88ssvOHbsGJycnODt7Y0vvvhCZuHKkAkTWDO7Fy9KPf3lgAEDULNmTSQmJmLZsmUKCpDTNWKxWDqVpMxQEf/8w/719y+8+TNXJLmTQI8ePTBt2jRMnz4dffr0Qffu3WUWeQQGBqJhw4awsLCAo6MjevTogYcPH8rsQ0SYN28eXF1dYWJiglatWuHu3bsy+2RlZWHChAmwt7eHmZkZunXrhlevXsn70riPmZqyIa4BYNEi1ou5hPT19fHjjz8CAJYvX47g4GBFRMjpmJs3byIuLg4WFhZo3rw5W0kEHDzIHvOiIPmppoqicB06dKDNmzfTnTt36NatW9SlSxdyd3eXqTxcsmQJWVhY0L59+ygyMpL69+9PLi4ulJxvBq4xY8ZQuXLl6NSpUxQWFkatW7emOnXqUG4xKzQ1tmJYE2RnE1WuzCrc5s0r1aHEYjENGjSIAJCNjQ09evRIQUFyuiIgIIAAUK9evfJW3rzJPp+mpkTp6eoLTsMorZ9AcYhL2LY8Li6OAFBoaKj0OM7OzrRkyRLpPpmZmWRlZUUbNmwgIqLExEQyNDSkoKAg6T6vX78mPT09On78eKHnyczMpKSkJOkSHR3Nk8Cn/P13Xl8EeYex+EhGRgY1btyYAFDVqlXpw4cPiomR0wkNGjQgAPTHH3/krZw7l30+8ycGTrGtg6pXr46dO3ciOzv7k/s9fvwYY8eOxdKlS0t0V5KUlASAtUAC2GiUsbGxaN++vXQfoVAIX19fXLp0CQC7PczJyZHZx9XVFbVq1ZLu87HAwECZvg3ly5cvUbw6o3dvVkGdmgqUciJ5Y2Nj/PPPP3Bzc8PDhw/x5ZdfIjc3V0GBcmXZ27dvpXNUdOrUKW+DpD6AFwWVSLGSwK+//opVq1bByckJ/fv3x/Lly7Fjxw7s27cPmzZtwtSpU9GoUSPUq1cPVlZWGDdunNyBEBGmTp2K5s2bS4eEjY2NBQA4OTnJ7Ovk5CTdFhsbCyMjI9jY2BS5z8dmzZqFpKQk6RIdHS13vDpFTw9YsoQ9Xr8eeP68VIdzdnbGoUOHYGpqihMnTmA6H/OdKwZJhfAXX3wBF8l0kVFRwO3brDL44+kluWIp1hzDfn5+uH79Oi5duoTdu3dj586deP78OTIyMmBvb4969eph6NChGDx4sOzsPnIYP348bt++jQsXLhTYJvhoOFgiKrDuY5/aRygUQljaXrW6pm1boE0b4MwZVlm8dWupDlevXj1s27YNffr0wc8//4yaNWti1KhRCgqWK4sKbRoqqRBu2RL4rwSBk5MqyqY+Z/z48eTm5kbPnj2TWf/06VMCQGFhYTLru3XrRkOHDiUiojNnzhAASkhIkNmndu3aNHfu3GKdn1cMF9O1awobalpiwYIFBIAMDAwoJCREIcfkyp7s7GyytLQkAHT58uW8Db6+7DO5erXaYtNUSusxrOAEhPHjx2P//v0IDg6Gl5eXzHYvLy84Ozvj1KlT0nXZ2dkIDQ2Fj48PAKB+/fowNDSU2ScmJgZ37tyR7sMpSMOGQJ8+rEneDz8o5JCzZ8+W1gv07t0bT58+VchxubLl4sWLSE5Ohr29PRo2bMhWxscD58+zx3I2T+fyUU1OKtzYsWPJysqKzp49SzExMdIlPV8zryVLlpCVlRXt37+fIiMjacCAAYU2EXVzc6PTp09TWFgY+fn58SaiyvLgQd5Q0xcuKOSQ6enp1LBhQwJA1atX5/8PXAHffvstAaAhQ4bkrdyyhX0O69ZVX2AaTK1NRIsLQKFL/pnLxGIxBQQEkLOzMwmFQmrZsiVFRkbKHCcjI4PGjx9Ptra2ZGJiQl27dqWXL18WOw6eBOQ0ahT78jVvXuqhpiVev35N5cqVIwDUuXPnYidwTjfUqFGDANCuXbvyVvbowT6HAQFqi0uTFfe6VqKJ5ssajZpoXhu8fg1UqgRkZgJHjiisVcbNmzfRokULZGRkYOrUqVixYoVCjstpt+fPn8PLywv6+vp49+4dawmYng7Y2wMZGUB4OBtinZOhtInmOQ7lygETJ7LHs2axCe0VoH79+tiyZQsAYOXKlfjjjz8UclxOu0laBfn4+OQ1BT99miUADw+gTh01Rqf95E4C+vr6iIuLK7D+/fv30OcDN+kOyVDTkZHArl0KO2y/fv0Q8N94RWPHjsW5c+cUdmxOOx39b4joQgeM696dDbnOlZjcSaCo0qOsrCwYGRmVOiBOS9jaskQAlHqo6Y/NnTsXffv2RU5ODnr37o2oqCiFHZvTLhkZGdLBBqX9A0QiNqE8wHsJK0CxOosBbAhpgHXc2rRpE8zNzaXbRCIRzp07h2rVqik+Qk5zTZwI/PIL60H8229s6GkF0NPTw5YtW/Ds2TPcvHkT3bp1w6VLl2ChqrkUOI0REhKCzMxMuLm5SUcSwKVLrHmojQ3QooV6AywDip0EVq1aBYDdCWzYsEGm6MfIyAienp7YsGGD4iPkNJdkqOkxY4CFC4HhwxU26Y2pqSkOHjyIhg0b4s6dOxg4cCD++ecfXuSoYyRFQV26dMkbAUDSS7hrV8Cg2JcwrijyNjtq1apVgd652o43ES2F/ENNz5+v8MNfvXqVjI2NCQB9++23Cj8+p7nEYjF5enoSADp06JBkJVHFiuzztnevegPUcLyJqBx4E9FS+vtvNuexuTnw7Bng4KDQw+/atQsDBw4EAGzevBnDhw9X6PE5zXTv3j3UrFkTQqEQ79+/h5mZGZvz2tsbEApZkVC+YmlOVnGva3LfS40YMeKT2//88095D8lpuz59gC++AMLC2FDT/xUdKsqAAQNw7949LFq0CKNHj0blypXRrFkzhZ6D0zySpqGtWrViCQDIKwpq144nAAWRu3XQhw8fZJa4uDgEBwdj//79SExMVEKInMbLP9T0unVsTmIFmz9/Pnr37o3s7Gz07NkTL5RwDk6zfLZpKKcQCikOEovFGDduHCpUqIDvvvtOEXGpFC8OUgAiNtx0cDAwbBjwX6cvRUpLS0OLFi0QHh6O2rVr4+LFizKt1LiyIykpCfb29sjNzcWTJ09QsWJF4NUroHx51i8gJgb4aJ4RTpZKewzr6elhypQp0hZEnA4SCPLuBrZtY2W3CmZmZoaDBw/CyckJt2/fxuDBgyEWixV+Hk79Tp06hdzcXFStWpUlAAA4dIj927QpTwAKpLBhI54+fcqnCdR1DRuyqSgVONT0x8qXL4+DBw9CKBTi4MGD+EFJ5+HU65NFQbyDmELJXRw0depUmb+JCDExMTh69CiGDRuGtWvXKjRAVeDFQQr08CFQsybr1XnhAqCkCtwdO3Zg8ODBAIBt27ZhyJAhSjkPp3pisRiurq54+/YtTp8+jTZt2gCJiazVWW4u8OgRULmyusPUeEorDgoPD5dZbt++DQBYsWIFVq9eXeKAuTKialXgq6/Y45kz2V2BEgwaNAizZs0CAHz99de4cuWKUs7DqV5YWBjevn0Lc3NztJD0CD52jCWA6tV5AlAwuZuIhoSEKCMOriwJCAC2b2d3Av/+q7QJwBctWoT79+/jn3/+QY8ePXDt2jW4u7sr5Vyc6kiKgtq1a5c3HpmkaSgvClK4EtcJxMXF4fz587hw4UKho4pyOszNLW8coVmzACVV3urp6eGvv/5CnTp18PbtW3Tv3h1paWlKORenOgUmlM/KYj8mAJ4ElEDuJJCUlIQhQ4bA1dUVvr6+aNmyJcqVK4fBgwcjKSlJGTFy2mjmTMDKig01vXOn0k5jbm6OQ4cOwdHREbdu3cKQIUN4iyEtFhcXh+vXrwMAOnXqxFaGhAApKYCLC9CggRqjK5vkTgKjRo3C1atXcfToUSQmJiIpKQlHjhzBjRs3MGrUKGXEyGmjj4eazs5W2qnc3d1x4MABGBkZ4cCBA5g7d67SzsUp17Fjx0BEqFevHlxdXdlKSVFQ9+6sYyKnWPIOSmRqakrnz58vsP7cuXNkamoq7+E0Ah9ATknS0ohcXNhgX7/8ovTTbd26VTpP9Y4dO5R+Pk7x+vXrRwBo9uzZbIVIlPcZOnZMvcFpmeJe1+ROq3Z2drCysiqw3srKKm/qN44D2FDTkl/lCxeyW3olGjp0qLTH+ogRI3iLIS2Tk5ODEydOAMjXP+DGDdY72MICaN1ajdGVXXIngdmzZ2Pq1KmIiYmRrouNjcW3336LOXPmKDQ4rgwYOZJNSv/uncIHlivM4sWL0a1bN2RlZaFr1664f/++0s/JKcbly5eRlJQEOzs7NGrUiK2UdBDr3JmNHMopnNxJYP369bhy5Qo8PDxQqVIlVKpUCe7u7rh06RI2btyIL774QrpwHAwNgUWL2OOffmLJQIn09fWxY8cONGrUCO/fv0eHDh0QHR2t1HNyiiFpGtqxY8e8yYP4gHFKJ3c/ge7du+fN8MNxxdG3L7BsmdKGmv6Yubk5jh49iubNm+Phw4fo2LEjzp8/D1tbW6Welyud/LOIAWA9g+/fZz8k8g8fwSkUn1QGfNgIlTh5EujQATAyYl9uDw+ln/Lly5fw8fHB69ev0bRpU5w+fRqmpqZKPy8nvxcvXsDT0xN6enp49+4dS9jLlwPffcfmDjh5Ut0hah2lDRtRoUIFvH//vsD6xMREVKhQQd7DcbqiXTvAz481FVVR3ZG7uztOnDgBGxsbXL58GX379kVOTo5Kzs3JR9JBrGnTpnl3bLyXsErInQSeP38OkUhUYH1WVhZevXqlkKC4MkggAAID2eO//gL27FHJaWvWrIkjR47AxMQE//77L0aOHMk7k2mgAr2E374FLl1ij7t1U1NUuqHYdQKHJGN5Azhx4oRMM1GRSIQzZ87Ay8tLsdFxZUujRqwn8ZIlrNVQvXqs5ZCS+fj4YM+ePejevTv++usvODo64qefflL6ebniycjIwJkzZwDkaxp6+DAbfLBBAzYMCac8xe14IBAISCAQkJ6envSxZDEyMqIqVarQ4cOHS9m9QT14ZzEVyskhataMdf6pV48oI0Nlp87fmWzZsmUqOy/3af/++y8BIDc3NxKLxWxl167sM7JokXqD02IK7ywmFoshFovh7u6OuLg46d9isRhZWVl4+PAhunbtqqRUxZUZBgZAUBBgZweEhwPTpqns1EOHDsXy5csBAN999x22bt2qsnNzRZMUBXXu3Jm1PExNBU6dYht501Clk7tOICoqCvb29go5+blz5+Dv7w9XV1cIBAL8I2kT/J/hw4dDIBDILE2aNJHZJysrCxMmTIC9vT3MzMzQrVs3Xjeh6dzcWL0AwCam//tvlZ16+vTpmD59OgBg5MiROHLkiMrOzRVERAVnETt5ko0cWrEim6CIUyq5+wksWLDgk9vlGbwrLS0NderUwVdffYXevXsXuk/Hjh2xefNm6d/S8cX/M3nyZBw+fBhBQUGws7PDtGnT0LVrV9y8eTOvwwmneTp1YsNMBwYCX38NfPGFSuoHAGDp0qWIi4vDtm3b0K9fP5w+fRo+Pj4qOTcn68GDB4iKioKRkRGbQQyQnUaS90lSPnnLmerWrSuz1KxZk0xNTcnS0pLq1atX0uIrAkAHDhyQWTds2DDq3r17kc9JTEwkQ0NDCgoKkq57/fo16enp0fHjx4t9bl4noCY5OUTNm6ulfiA7O5u6dOlCAMja2poiIyNVdm4uz08//UQAqF27dmxFdjaRjQ37TJw7p97gtJzSBpD7eHrJO3fuICYmBm3atMGUKVMUnKKAs2fPwtHREVWqVMGoUaNkJrC5efMmcnJy0L59e+k6V1dX1KpVC5ckzcsKkZWVheTkZJmFUwMDA2DXLsDentUPfDR/tTIZGhri77//ho+PDxITE9GhQwe8ePFCZefnmAK9hC9cAD58YJ8JfnemEgoZnNvS0hILFixQ+ABynTp1wo4dOxAcHIwVK1bg+vXr8PPzQ1ZWFgA2cJ2RkVGB0UudnJwQGxtb5HEDAwNhZWUlXcqXL6/QuDk55K8fWL8e2L1bZac2NTXF4cOHUbNmTbx58wYdOnRAfHy8ys6v65KTk3H+/HkA+ZKApCjI3x/gxbkqobAZGiQTzChS//790aVLF9SqVQv+/v44duwYHj16JP31UBQi+uT4RrNmzUJSUpJ04QOMqVnHjsD337PHo0YBjx+r7NS2trY4fvw4ypcvj4cPH6Jz585ITU1V2fl12alTp5Cbm4vKlSujUqVKrF8A7yWscnJXDP/yyy8yfxMRYmJi8Ndff6Fjx44KC6wwLi4u8PDwwOP/LhLOzs7Izs7Ghw8fZO4G4uLiPlnRJxQKIeTD0mqW+fOB8+fZ0q8fcPkyYGysklO7ubnh5MmTaN68Oa5fv47evXvj8OHDBRohcIpVoCgoIgJ48YLNQ9GunRoj0zHyVjZ4enrKLBUqVKDGjRvTrFmzKDk5uYRVGIVXDH8sPj6ehEIhbd26lYjyKoZ3794t3efNmze8YlhbvXpFZG/PKgXHjlX56a9cuUKmpqYEgAYMGEAikUjlMegKkUhEzs7OBIBOnTrFVgYEsP/7Hj3UGltZUdzrmtxJQJFSUlIoPDycwsPDCQCtXLmSwsPD6cWLF5SSkkLTpk2jS5cuUVRUFIWEhFDTpk2pXLlyMslmzJgx5ObmRqdPn6awsDDy8/OjOnXqUG5ubrHj4ElAgxw/zi4EAFG+Vl+qO/1xMjAwIAA0ceLEvB6snELduHGDAJCZmRllZmaylXXrsv/3LVvUG1wZoZIkEB0dTa9evSrx80NCQqTd+PMvw4YNo/T0dGrfvj05ODiQoaEhubu707Bhw+jly5cyx8jIyKDx48eTra0tmZiYUNeuXQvs8zk8CWiY779nFwNzc6JHj1R++u3bt0s/iz/++KPKz68LFixYQACoh+RXf1QU+z/X0yN6906tsZUVSksCIpGI5s+fT5aWlqSnp0d6enpkZWVFCxYs0NrbZ54ENExODlHLluyiUKeOSvsPSKxevVqaCH7//XeVn7+sa9y4MQGg3377ja34+Wf2/+3rq9a4yhKlJYGZM2eSg4MDrVu3jiIiIujWrVv066+/koODA33//fclDlideBLQQK9fEzk4sAvDmDFqCWHWrFkEgPT09D5bX8UVX1xcHAkEAgKQV5LQujX7v161Sq2xlSVKSwIuLi508ODBAuv/+ecfcnV1lfdwGoEnAQ114gSRQMAuDrt2qfz0YrGYRo4cSQBIKBRSaGioymMoi7Zt20YAqG7dumzF+/dE+vrs//nZM/UGV4YorcdwQkICqlWrVmB9tWrVkJCQIO/hOK5o7dvL9h949EilpxcIBNiwYQO6d++OrKws+Pv7IyIiQqUxlEUFBow7ehQQiYDatQE+J4nKyZ0E6tSpg7Vr1xZYv3btWtSpU0chQXGc1Lx5gK8vG164Xz8gI0OlpzcwMMCuXbvQokULJCcno2PHjnj27JlKYyhLcnNzceLECQCF9BLmw0arhdwTzYeGhqJLly5wd3dH06ZNIRAIcOnSJURHR+Pff/9FixYtlBWr0vCJ5jXcmzdA3brAu3fA6NHAhg0qDyExMREtW7ZEZGQkKlWqhIsXL8LR0VHlcWi78+fPo2XLlrC1tUVcXBz0s7MBBwcgLQ24eZONJssphNImmvf19cXDhw/Rs2dPJCYmIiEhAb169cLDhw+1MgFwWsDVFdi+nQ0rvHEjG3ROxaytrXH8+HF4enriyZMn6NSpEx94sAQkRUEdO3ZkQ72fOcMSQPnybLpRTvWKW8mwceNGiomJKWVVhWbiFcNaYvbsvP4DDx6oJYRHjx6Rg4MDASA/P7+8jk5csXh7exMA2rFjB1sxciT7Px0/Xr2BlUEKrxjetWsXPD090bhxYyxevBj37t1TVl7iuMLNmwe0aqW2+gEAqFy5Mo4dOwZzc3MEBwdjyJAhIPlKVHXW7du3ERkZCT09PXTo0IFVBh8+zDbyAePUpthJICQkBDExMZgwYQJu3bqFpk2bomLFipg6dSrOnj0LsViszDg5jg0tvHMn4OgI3L4NTJ6sljDq16+PAwcOwNDQEHv27MGZM2fUEoe2WbRoEQCgd+/esLOzA65cAeLiAGtroGVL9Qanw+SqE7CxscHgwYPx999/4927d/j111+RmZmJIUOGwMHBAUOHDsXevXuRlpamrHg5Xefiklc/8NtvLCmoQdu2bTFq1CgAwK+//qqWGLTJvXv3sHfvXgDIm3dEMmx0ly6AoaGaIuMUNoDcjRs3aM6cOVS7dm1asGCBog6rErxOQAvNmaP2+oG7d+9KexQ/f/5cLTFoi4EDBxIA6tWrF1shFhNVqsT+D/fsUW9wZZTSOostWLAA6enpBdbXqFEDBgYGiIiIwMyZM0udnDjukwIC1F4/UKNGDfj5+UEsFmODGpqtaouHDx8iKCgIADB79my28v594MkTwMgI6NBBjdFxcieB+fPnFzrzUnp6OubPnw+Azd/KcUr1cf3ApElqCWP8+PEAgE2bNiEzM1MtMWi6xYsXQywWw9/fH/UkzUAlRUFt2wIWFuoLjpM/CVARUzdGRETA1tZWIUFxXLG4uAA7drD6gd9/Z49VzN/fH+XLl0d8fDz+/vtvlZ9f0z158gQ7/vt/kZmDnPcS1hjFTgI2NjawtbWFQCBAlSpVYGtrK12srKzQrl079OvXT5mxclxBbdsCkovL6NHAw4cqPb2BgQHGjh0LAIUOp6LrAgMDIRKJ0KlTJzRs2JCtfPMGuHaNJe9u3dQbIFf8YSO2bt0KIsKIESOwevVqWFlZSbcZGRnB09MTTZs2VVqgysSHjdByIhGbkzYkBPD2Bq5eBUxMVHb6d+/ewc3NDdnZ2bh69SoaNWqksnNrsufPn6Ny5crIzc3FpUuX8q4PGzYAY8cCTZsCly6pN8gyrNjXNXlrnM+ePUvZ2dklqq3WVLx1UBkQE0Pk5MRam3z9tcpPP2TIEAJAQ4YMUfm5NdXo0aMJALVr1052Q4cO7P9pyRL1BKYjintdk3sAuZcvX35yu7u7uzyH0wj8TqCMOHOG3REQAX/9BQwerLJTX7t2DY0bN4aRkRGio6N1fnC56OhoVKxYETk5OTh//jyaN2/ONiQnA/b2QE4O8OABULWqegMtw5Q2gJynpye8vLyKXDhObdq0AebOZY/HjAFu3FDZqRs1aoSGDRsiOzsbmzZtUtl5NdXSpUuRk5ODVq1a5SUAADh2jCWAatV4AtAQcieB8PBwhIWFSZerV69iw4YNqFKlCvbs2aOMGDmu+ObMAVq3ZiNT+vgAS5eyOgMVkDQX3bBhA3Jzc1VyTk30+vVr/P777wCAuZKkLCFpGspbBWkORZU/HTlyhHy1dJJoXidQxrx/T9SjByt3Btik9Sro0ZuRkUH29vYEgPbv36/082mqSZMmEQBq3rw5icXivA3p6USWluz/5PJl9QWoI5TWY7goVapUwfXr1xV1OI4rOVtbYP9+YNMmwMwMOHeOTV34118sLSiJsbGxdDwhXW0uGhsbi40bNwJgdwEyfYo2bmR1Ah4eAG9BpTHkTgLJyckyS1JSEh48eIA5c+agcuXKyoiR4+QnEAAjRwIREawpYnIyMHQo0L8/oMS5sMeMGQM9PT0EBwfr5HDrK1asQGZmJpo0aYK2bdvmbUhNBRYvZo9nzwb0FPb7kysteW8xBAIB6enpySwCgYDc3d3p0qVLJb51USdeHFTG5eQQLVxIZGDAiiJcXYlOnVLa6Xr06EEAaNy4cUo7hyaKi4sjU1NTAkD//vuv7MbFi9l7X7EiURlrYq6plNZENDQ0VOZvPT09ODg4oFKlSjAwMFBYclIl3kRUR1y/zpqNPnrE/p40CQgMVHjHsjNnzqBt27YwNzfH69evdeYzNWvWLCxZsgQNGjTAtWvX8oqCEhMBLy/27/btwKBB6gxTZyits1hZxO8EdEhqKtHYsXmVxjVqEIWHK/QUYrGYqlWrRgBozZo1Cj22poqPjydzc3MCQAcPHpTdOHdu3nudm6ueAHWQUiuGHz58iPHjx6NNmzZo27Ytxo8fjwcPHpTkUBynWmZmwLp1wNGjgJMTcO8eq6RctkxhTUkFAoG0uejatWt1YvrJn3/+GampqahTpw78/f3zNsTHA6tWsccLFrDRXzmNIncS2Lt3L2rVqoWbN2+iTp06qF27NsLCwuDt7c37CXDao3NnIDKStVfPyQFmzAD8/IAXLxRy+KFDh8LCwgIPHz4s89NPJiYm4ueffwZQSIug5cuBlBSgXj2gZ081Rch9kry3GF5eXjRnzpwC6+fOnUteXl7yHk4j8OIgHSYWE23aRGRmxoosLC2Jtm1j60tp/PjxBIC6d+9e+jg12Pz58wkA1apVi0QiUd6GmBgiExP2vh45or4AdVRxr2tyJwETExN6/PhxgfWPHj0iExMTeQ+nEXgS4OjJE6KmTfPqCvr2ZZ3OSuH+/fvS6SejoqIUE6eGSUpKImtrawJAu3fvlt04cSJ7L5s0UUhS5eSjtDqBVq1a4fz58wXWX7hwAS1atJDrWOfOnYO/vz9cXV0hEAjwj2Siiby7FMybNw+urq4wMTFBq1atcPfuXZl9srKyMGHCBNjb28PMzAzdunXDq1ev5H1ZnK6rWJF1Klu4kJVb79nDhqU+fbrEh6xWrRratGlTpqefXLt2LRITE1GtWjX07t07b8PLl2zIaABYtIj12+A0k7zZZf369eTg4EDffPMN/fXXX/TXX3/RN998Q46OjrR+/Xo6ePCgdPmcf//9l3744Qfat28fAaADBw7IbF+yZAlZWFjQvn37KDIykvr3708uLi6UnJws3WfMmDFUrlw5OnXqFIWFhVHr1q2pTp06lCtHKwR+J8DJuHaNqEqVvLuCSZPYkAclcODAAQJAdnZ2lJGRodg41SwlJYXs7OwIAG3fvl1246hR7L1r1YrfBaiJ0oqDBAJBsRY9PT25jvtxEhCLxeTs7ExL8o05npmZSVZWVrRhwwYiIkpMTCRDQ0MKCgqS7vP69WvS09Oj48ePF/vcPAlwBXzclLRmTaJbt+Q+TE5ODrm7uxMA2rJlixICVZ9ly5YRAKpcuTLl5OTkbXj8mEhfn71vFy6oL0Adp7TiILFYXKxFVMrmdlFRUYiNjUX79u2l64RCIXx9fXHpv9mIbt68iZycHJl9XF1dUatWLek+hcnKyiow/AXHyfi4Kendu0DDhnI3Jc0//eSaNWvKTHPR9PR0LF++HADwww8/yHYUXbCAvUedOgHNmqkpQq64NHYAj9jYWACAk5OTzHonJyfpttjYWBgZGcHGxqbIfQoTGBgIKysr6VK+fHkFR8+VGQpoSjpy5EgIhULcvHkT165dU2KwqrNx40a8e/cOFSpUwMCBA/M23LvHegUDrH6F03glGufhzJkzOHPmDOLi4iAWi2W2/fnnnwoJTELwUYUSERVY97HP7TNr1ixMnTpV+ndycjJPBFzRHByAAweAP/9kQ01IRiX99Vc2BMJnPo8ODg748ssvsXXrVqxduxaNGzdWUeDKkZGRgWXLlgEAvv/+exgaGuZtnDePFaD17AnUr6+eADm5yH0nMH/+fLRv3x5nzpxBfHw8Pnz4ILMoirOzMwAU+EUfFxcnvTtwdnZGdnZ2gfPm36cwQqEQlpaWMgvHfVJho5IOGcKWYhQPSXoQ//3333j79q2yo1WqTZs2ITY2Fu7u7hgyZEjehlu3WKsqgQCYP19t8XFykreywdnZmbZt21aymopPQBEVw0uXLpWuy8rKKrRiOH/75Ddv3vCKYU65JKOSSio/i1nh26hRIwJAixYtUnKAypOZmUnlypUjALR+/XrZjf7+7P0YMEA9wXEylNY6yNbWlp48eVLiwPJLSUmh8PBwCg8PJwC0cuVKCg8PpxcvXhARayJqZWVF+/fvp8jISBowYEChTUTd3Nzo9OnTFBYWRn5+fryJKKcay5axi1758sVqQrpt2zYCQG5ubrKtabTI+vXrCQCVK1eOMjMz8zZcucLeCz09oocP1RcgJ6W0JPDdd9/RggULShxYfiEhIQSgwDJs2DAiYncDAQEB5OzsTEKhkFq2bEmRkZEyx8jIyKDx48eTra0tmZiYUNeuXenly5dyxcGTAFciGRksAQBE+e5Yi949gxwcHAgA7du3TwUBKlZWVpa0uWuB0VHbtmXvw1dfqSc4rgClzScwadIkbNu2DbVr10bt2rVlK4UArFy5suRlU2rC5xPgSmzbNmDYMMDKCnj6FLCz++TuP/zwAxYvXozWrVsjODhYRUEqxqZNmzBq1Ci4uLjg2bNnMDY2ZhvOngVatwYMDdlcDZ6e6gyT+09xr2tyJ4HWrVsXfTCBQOs+2ABPAlwpiETAF18At28DU6cCK1Z8cvfo6Gh4enpCLBbjzp07qFmzpooCLZ2cnBxUrVoVUVFRWLVqFSZPnsw2EAEtWwIXLgDjxrEWU5xG4JPKyIEXB3Glcvw4KwoxMiIqxkBxvXr1IgA0duxY5cemIFu2bCEA5OjoSGlpaXkbJK/d2Jjo9Wv1BcgVoNRJZTiOy6d9e6BNGyA7G5gz57O7S5qLbtu2DUlJScqOrtRyc3Px448/AgCmT58OU1NTtoGITRoPsLsAV1c1RciVRrGKg3r16oUtW7bA0tISvXr1+uS++/fvV1hwqsKLg7hSu3kTaNCAPQ4LY5OoFIGIUKtWLdy7dw8///wzJk6cqKIgS2bHjh0YPHgw7Ozs8Pz5c5ibm7MNBw8CPXqwITaePQMcHdUaJyeruNe1Yt0JWFlZSXvg5h9uobCF43RS/fqAZPiEGTM+uatAIMA333wDAPj1118L9LrXJCKRCIsWLQIATJs2LS8BiMV5dz2TJvEEoMXkrhgui/idAKcQUVFAtWqsWOjECVZMVISUlBSUK1cOKSkpOHnyJNq1a6fCQItv9+7d+PLLL2FjY4Pnz5/nfT927wa+/JK1ioqKAj4av4tTP4XeCXAcVwxeXsB/v/Dx3Xfs13IRLCwsMHz4cABsYhZNJBaLpXcBU6ZMybuQ5OYCc+eyx9Om8QSg5XgS4DhF+uEH9us4IgLYseOTu44bNw4AcPjwYTx//lwFwcnnn3/+wZ07d2BlZYUJEybkbdi+nfUHsLNjRUGcVuNJgOMUyc4OmDmTPZ49G8jMLHLXatWqoV27diAirF+/XkUBFg8RYcGCBQCAiRMnwtramm3Izs4bHG7mTIAXn2o9ngQ4TtEmTQLc3Ng8u5/pPCVpLrpp0yZkZGSoIrpiOXz4MCIiImBubp7XMQxgw2k/fw44O7NmoZzWkysJ5OTkoHXr1nj06JGy4uE47WdiwmbXAoAffwQ+McR6ly5d4OHhgYSEBOzevVtFAX5a/ruACRMmwNbWlm3IyMibKOaHHwBJfwFOq8mVBAwNDXHnzp3PTurCcTpv6FCgVi2WAAIDi9xNX19f46afPHbsGG7evAkzMzOZyZewcSPw5g1QvjwwapT6AuQUSu7ioKFDh+KPP/5QRixa5/VrIC1N3VFwGklfH1i6lD3+5RdWNFQEyfSTYWFhuHr1qooCLFz+u4Bx48bB3t6ebUhNzUtmc+cCQqGaIuQUTe7pJbOzs7Fp0yacOnUKDRo0gJmZmcx2bRxFtCRSUoCOHdnAiQcPsh9HHCejUyc2umZICOtYtXVrobvZ29tjwIAB2LJlC9auXYsmTZqoONA8p0+fxtWrV2FiYoJp06blbVi7FoiLAypWZKOmcmUGH0UUJessdvs2Gy4mPh5wcgL++QdQ43eX01Q3bgANG7IpF8PDgTp1Ct3t5s2baNCgAQwNDREdHf3J6VGVhYjQokULXLx4EZMnT8aqVavYhqQk1gfiwwfgr7+AwYNVHhsnPz6KqBxKOopoVBSRtzcbRFEoJFLCrJtcWdC/P/uQdOjwyd2aNGlCAGjhwoUqCkxWcHAwASChUEiv848IOncui79GDSI5Zuzj1Evpo4g+efIEJ06ckDZrIw2o0FI1T0/g0iWge3cgK4vVBc6c+cmOopwu+vFHVm544gRw+nSRu0mai27YsAG5ubmqik5KUhcwatQouEpGBI2PByR3BAsWsLoOrmyRN7vEx8eTn58fCQQC0tPTo6dPnxIR0YgRI2jq1KklyljqVtr5BEQiou+/Zz+WADbfdr5pkDmOaOJE9uGoV499YAqRmZlJjo6OBID27t2r0vAkdwFGRkYUHR2dt+G77z4bN6eZlHYnMGXKFBgaGuLly5d544oD6N+/P44fP66w5KRN9PTYj70dO1ijicOHAR8fNq4WxwFgvYctLVm9QFBQobsIhUKM+q/ppSrGE0pPT8eOHTvQvn17tGnTBgAwYsQIuLm5sR1iY4E1a9jjhQvZB50re+TNLk5OTnTr1i0iIjI3N5feCTx79ozMzMxKkK/UT5Ezi129SuTiwn482dkRhYYqIECubPjxR/bB8PQkyswsdJeXL1+Svr4+AaDIyEiFhyASiSg0NJRGjBhBFhYWBEC6tG3bluLi4vJ2lty9NGlCJBYrPBZt9/33RI0bE/32G1FGhrqjKai41zW5k4C5uTk9evRI+liSBK5du0a2trYlCFX9FD295KtXRPXrs++PgQHR778r5LCctktLI3J1ZR+MlSuL3K13794EgMaMGaOwUz99+pQCAgLIy8tL5sLv5eVF8+bNk36PpV6+ZNNlAkSnTyssjrJi+/a84l+AyNmZKDCQ6MMHdUeWR2lJoHPnzjR79mwiYkng2bNnJBKJqG/fvtS7d++SRatmyphjOC0tr1EIQDRpElFOjsIOz2mrTZvYB8LWtsgrRkhICAEgU1NT+lCKq0pSUhJt2rSJWrRoIXPht7CwoBEjRlBoaCiJiirn/9//WJytWvG7gI88fkxkbs7enu7didzc8r7nFhZE06ezH4LqprQkcPfuXXJwcKCOHTuSkZER9enTh6pXr05OTk705MmTEgesTsqaaF4sJlqwIO8D0r69Zv1S4NQgJ4c1tQSIZswodBexWEw1a9YkALR69Wq5Dp+bm0snTpyggQMHkomJifTCLxAIqF27drR9+3bZieIL8+QJu4UFiM6fl+v8ZV1WVt5dvq8vazGblUW0dStRzZp533VDQ6KvviK6e1d9sSotCRARxcTE0Ny5c6lLly7UqVMn+uGHH+jNmzclClQTKCsJSOzdS2Rqyj4cVasSPXyolNNw2uLwYfZhMDZmxS6FWL9+PQGgKlWqFP1rPZ979+7RjBkzqFy5cjK/+qtVq0aBgYGyLX4+Z8gQFl/HjsV/jo6YOjXvRu7jt1QkIjpyhKhlS9miIn9/ogsXVB+r0pLAixcvSFzE7eGLFy/kPZxGUHYSICIKDycqX559KKytiU6eVNqpOE0nFuddKYYPL3SXlJQUsrS0JAB04sSJQveJj4+ntWvXUsOGDWUu/DY2NjRu3Di6evVqkd/VIt29SyQQsNiuX5f3lZVpR4/mXdgPHfr0vpcvE/XsmfdWAkQ+PkT//KO6lrZKSwJ6enr09u3bAuvj4+NJT09P3sNpBFUkASKi2Fiipk3ZB0Jfn2jNGl7cqrOuXmUfBIGA6PbtQneZOHEiASB/f3/puuzsbDp48CD16tWLDA0NpRd+fX198vf3p71791JmES2PiqVvXxZXz54lP0YZ9Po1kb09e2smTiz+8x48IBo1Kq+OHSCqVo3ojz+KbCCmMEpLAgKBQLYZ2X+eP39Opqam8h5OI6gqCRCx//ihQ/M+EKNHE2VnK/20nCaSXHA7dSp084MHD6Tl+YcOHaJJkyaRg4ODzK/+unXr0qpVqwr9YSa38PDPJiZdlJtL1Lp1Xp+5kly837whmjmTyMoq77vv6kq0bBlRYqLCQyai4l/Xij2AnGRc8Z9//hmjRo2S6SgmEolw9epV6Ovr4+LFi8XtoqAxSjKAXGkQAStWsLnIiYBWrYC9e9nMhJwOefIEqF6dTdweHMxGHP1Ihw4dcPLkSZl1Tk5OGDRoEIYNG4batWsrLp5u3VhPxwEDgJ07FXdcLffjj6yvn5kZEBYGVKlS8mMlJwO//cZG4njzhq2ztATGjAEmTwZcXBQS8n/nUvAAcq1ataJWrVqRQCAgHx8f6d+tWrWi9u3b0//+9z9p/wFto8o7gfwOH2ZNygCiChWI7txR6ek5TTB+PPsANGhQaGHxyZMnSSAQkJGREfXt25eOHDlCOcpoa3zlCotDT4+3XMjnwgVWdAuwFkCKkpVF9OefRNWr590ZGBkRjRxJdP++Ys6htOKg4cOHq/xiqWzqSgJE7MJfoUJeG+PDh1UeAqdOb9/mNTrftavQXR49ekQJCQnKOX9KCtHChUSWliyGr75Sznm00Pv3RO7u7G0ZPFg55xCJWCVzs2Z5yUAgIOrRg+jSpdIdW6lNRFUlICBApvwTADk5OUm3i8ViCggIIBcXFzI2NiZfX1+6U4Kf0+pMAkRE796xNseSD8CyZbzCWKcsXMj+8728lF9bKJGVRbR2LZGTU97V54svWOE1R2IxqxsHiCpVUs2AkBcusM5n+ZuXtmjBfhiWpEWRUpPAtWvX6Ntvv6X+/ftTz549ZRZFCggIoJo1a1JMTIx0yV8pvWTJErKwsKB9+/ZRZGQk9e/fn1xcXChZzv8xdScBIvadlHTSBIiGDVPd9YBTs9TUvAGn5OwcJjeRiI154OWV92GrVIkoKIiPEprPunV5nb5u3FDtue/dIxoxgp1b8l9Us6b8iUhpSWDXrl1kaGhIXbp0ISMjI+ratStVrVqVrKysaHgRbZ5LKiAggOrUqVPoNrFYTM7OzrRkyRLpuszMTLKysqINGzZ88riZmZmUlJQkXaKjo9WeBIjYr481a/LKIJs2JYqJUWtInKps3Jg36qAymouIxayhe+3asgPerF/Pm6d9JCKCTRIFEK1apb44Xr9mI3lbWhK1ayf/85WWBLy9vWnt2rVElDeAnFgsplGjRtHcuXPlj/QTAgICyNTUlFxcXMjT05P69+8vHejq6dOnBIDCwsJkntOtWzcaOnToZ4/7cTGTJiQBiZMnWYcygHUw++glcsWgdcVpOTmsATlANGuWYo998SIrV5Bc/K2siBYvZncgnIzU1Lz/hi5dNONzlJhI9PH4fsWh8CaiEmZmZrh79y48PT1hb2+PkJAQeHt74/79+/Dz80NMTIw8h/ukY8eOIT09HVWqVMHbt2+xaNEiPHjwAHfv3sXDhw/RrFkzvH79Om8WJAD/+9//8OLFC5w4caLI42ZlZSErK0v6d3JyMsqXL6+yJqLF8egR4O/P/jU1Bbp2BZydWRMyySL5285O94Z6z80FoqPZnA2FLQkJbJ73//0P6NBBSybEOngQ6NEDMDEBHj8GypUr3fHu3AF++AE4dIj9bWwMTJjApr+ztS11uGXR118Df/wBuLoCt24BDg7qjqjkittE1EDeA9va2iIlJQUAUK5cOdy5cwfe3t5ITExEenp6ySMuRKdOnaSPvb290bRpU1SsWBFbt25Fk/9mdRcIBDLPIaIC6z4mFAohFAoVGquiVakCXL0K9O8PnDwJ/P130fsaGLCEUFSSkCxOToCRkepeQ2kQsTlNirrIR0cDItGnj3HwIFvKlwdGjgRGjGCPNVa3bkDz5sCFC0BAALBpU8mO8+IFe/62beyN1NNjLz4gAJBMGMMVEBTEEoBAAGzfrt0JQB5yJ4EWLVrg1KlT8Pb2Rr9+/TBp0iQEBwfj1KlT0tmJlMXMzAze3t54/PgxevToAQCIjY2FS74eFnFxcXByclJqHKpibQ0cPQocPw48fQrExOQtsbHs3/h49qv41Su2fI6dXdFJws6OzYxmaMiShZFR3uPC/jU0ZF+YkkpMlL2wP3uW9/j5cyAz89PPNzJi8zx7eRVcDAzYNXDrVpYw5s1jU+R27szuDjp1YvtoFIEAWLaMTUu3eTMwZQpQs2bxn//uHevZtH49kJ3N1vXuDSxaBFSrppyYy4hnz9jnAmAdwwrpt1dmyV0clJCQgMzMTLi6ukIsFuOnn37ChQsXUKlSJcyZMwc2NjbKihVZWVmoWLEi/ve//2HOnDlwdXXFlClT8N133wEAsrOz4ejoiKVLl2L06NHFPq6qewwrUnY28PZtXlL4OEnk/1sZc5cbGHw6WRSWON69Yxf6xMRPH1tPj/1wLewi7+XFEtfnisEyM4H9+4HffwfOns1bX64cuzsYORJwdy/tu6BgffoA+/axMsDDhz+/f0oKsHIl8NNPQGoqW+fnByxZAjRsqNxYy4DsbHYDdv06+zckRAN/IJRAca9rcicBVZo+fTr8/f3h7u6OuLg4LFq0CKGhoYiMjISHhweWLl2KwMBAbN68GZUrV8bixYtx9uxZPHz4EBYWFsU+jzYngeISi1k5+aeSxPv3QE4O+1IU9m9OjuLjcnQs+iJfvrxii68ePmQlLFu2sDsogP347tiR/Qrs0oUlKbV79AioUYOVd509C/j6Fr5fVhawcSP7pf/uHVv3xRfs4t+uncrC1XbffQcsXw7Y2LB6AI37UVBCSksCgwYNQqtWreDr64sqpRlEoxi+/PJLnDt3DvHx8XBwcECTJk2wcOFC1KhRAwAr/58/fz42btyIDx8+oHHjxvj1119Rq1Ytuc6jC0lAEYjY3UT+pFBUwvjUvzY27CLv6cnGY1G1rCzgn3/YGC7BwXnrXVxY0fnIkSw+tRo3jhXrNGzIKofyl7uJRGxsn7lzWbkZAFSuzIqCevfWvVYCpXD8OCsaBIADB1i9fFmh8LGDJP73v/9R1apVSSAQkIuLC3355Ze0fv16uq+oAS/UQBM6i3Hq8fgxm+DL0VG2236HDmwyILU1oY+NJTIzYwHt3s3WicVsjIFatfKCdXFhfQx4W3+5vXlD5ODA3sZvvlF3NIqn9GEjYmJiaNeuXTR69GiqVq0a6enpkbOzc0kPp1Y8CXBZWUR79rBOOfm77Ts5sSGA1TJz6rx5LIiKFYmCg2UHmLG2JlqyhE1mzclNJCJq25a9lbVrE2VkqDsixSvuda3E940WFhawsbGBjY0NrK2tYWBgAGdn55IejuPUysiI1ceePMlaYs2axVpPvX3LitgrVWLF7Hv25DW8Ubpp01i73qdPWUXvxYusrf+MGaw5y4wZrBMJJ7dly4DTp9nbFxTE3lZdJXedwIwZMxAaGoqIiAjUqlULLVu2hK+vL1q2bAlra2slhalcvE6AK0xODnDkCKs7OHGC/QQHWPvx4cOBUaNYUbxS/fYbMHo06+329dfAnDml70Sm4y5fBlq0YFUrf/zB6oHKIqVVDOvp6cHBwQFTpkxB9+7dUb169VIHq248CXCf8/w5u2D88QdrTSXRrh2wcCHQuLGSTkzEMlClSmzhSiUxEahbl/WnGzAA2LGjdH1dNJnSkkBERARCQ0Nx9uxZnD9/Hvr6+vD19UWrVq3QqlUrrUwKPAlwxZWbyzrw/f478O+/eXcHvXuzxjlVq6o3Pq5oREC/fmwWvwoVgPBwNqtXWaWyfgIRERFYvXo1tm/fDrFYDNHn+vJrIJ4EuJJ4/pz1Qt66lfXD0NdnzUsDAtjYM5xmkZSsGRgAly6V/X50xb2ulahiODw8HKtWrUL37t3RunVr/PXXX6hTp450HmKO0wWensCffwIREWywP5GIXWgqVWLjtiUlqTtCTuLOHWDSJPY4MLDsJwB5yH0nYGNjg9TUVNSpU0daBNSyZUut/gXN7wQ4RbhwgTXYuXSJ/W1rC3z/PfDNN7rd+kTd0tOBRo2Au3dZ7/CjR3WjP53SioOOHDmi9Rf9j/EkwCkKERu5edYs4P59tq58eVZ5PHiwlgxpXcaMHs3u0Jyd2V2bo6O6I1INpRUHde3alV8oOa4IAgHQvTtw+zZrSeTmxkYxHT6ctUo5ciSvMplTvj17WAIQCIC//tKdBCAPHbgp4jjVMzBg7c8fPWIdk6ytWbm0vz8bD+7yZXVHWPY9f876cgBsHp22bdUajsbiSYDjlMjEBPj227wOvsbGwPnzbMqAnj3ziow4xcrJYf0AkpKApk2B+fPVHZHm4kmA41TAxoYNP/H4Mev4q6fHRjKtVYv9XZwJgbjimzsXuHIFsLJiA65qxBDhGoonAY5TITc31tHszh12JyAWs7qDypXZncKHD+qOUDslJrJxnxYuZHPxLF3K1m/axJryckXT6EllVIW3DuLU5fJldvE/f579bW3NWhZNmMCKkriCcnJYxfvVq3nLw4cF95swAfjlF9XHpynKxMxiqsKTAKdORGwIipkz2R0CwMaImz8fGDasbEx1WFJEbJyf/Bf8sLDC55+uWJGN4dS4MasHaNCg7I4LVBw8CciBJwFOE4hEwPbtrDz75Uu2rnp11sO1WzfduKAlJbG5fvNf9OPiCu5nY8M6gEku+o0aAfb2qo9Xk/EkIAeeBDhNkpkJrFvHBqRLSGDrLC1ZQqhenU0/LHns6am9HdByc4HISNkL/oMHBftRGBiwPhaSC37jxqwORReSYmnwJCAHngQ4TZSUxPoY/PwzkJZW+D7Gxmzk0o8TROXKbKIcdSNiiSwmJm+RlOffvAlkZBR8jqcnu9A3acL+rVePD7tREjwJyIEnAU6TZWezpqX37rF+BZLlwQMgK6vw5+jrszLy/HcNNWoA1aoBZmaljyknh826Fhsre4Ev7O+cnKKPY2kpW6zTuDHv1asoPAnIgScBThuJRKxX7P37BRNEcnLRz3N3L1isVL06YGfH7jg+dUGXPI6Pl2/4C1tbwMWFLZUq5V3wq1bVjcHc1IEnATnwJMCVJUTAmzd5CSF/giisklXCxKTw4pmi6OuzQdmcnfMu8C4uBf92cgKEwtK/Lk4+xb2u6XDjM44rmwQC1sS0XLmC4+W8fy97xyBJEC9f5iUAU9NPX9Qlf9vb81/xZQFPAhynQ+zsgObN2ZJfaior43d0BMzNecsbXcKTAMdxMDdnC6d7+M0cx3GcDuNJgOM4TofxJMBxHKfDeBLgOI7TYTwJcBzH6bAykwTWrVsHLy8vGBsbo379+jgvGaCd4ziOK1KZSAK7d+/G5MmT8cMPPyA8PBwtWrRAp06d8FIyHi/HcRxXqDIxbETjxo3xxRdfYP369dJ11atXR48ePRAYGFhg/6ysLGTlG3krOTkZ5cuX58NGcBxXZhR32AitvxPIzs7GzZs30b59e5n17du3x6VLlwp9TmBgIKysrKRL+fLlVREqx3GcxtH6JBAfHw+RSAQnJyeZ9U5OToiNjS30ObNmzUJSUpJ0iY6OVkWoHMdxGqfMDBsh+GiwEyIqsE5CKBRCmG9YQ0mJWPKnxt/lOI7TIpLr2edK/LU+Cdjb20NfX7/Ar/64uLgCdwdFSUlJAQBeLMRxXJmTkpICKyurIrdrfRIwMjJC/fr1cerUKfTs2VO6/tSpU+jevXuxjuHq6oro6GhYWFgUefegTpKK6+joaK2vuOavRTOVpdcClK3XU9LXQkRISUmBq6vrJ/fT+iQAAFOnTsWQIUPQoEEDNG3aFL/99htevnyJMWPGFOv5enp6cHNzU3KUpWdpaan1H2gJ/lo0U1l6LUDZej0leS2fugOQKBNJoH///nj//j0WLFiAmJgY1KpVC//++y88PDzUHRrHcZxGKxNJAADGjRuHcePGqTsMjuM4raL1TUR1gVAoREBAgEyLJm3FX4tmKkuvBShbr0fZr6VM9BjmOI7jSobfCXAcx+kwngQ4juN0GE8CHMdxOownAY7jOB3GkwDHcZwO40lAQwUGBqJhw4awsLCAo6MjevTogYcPH6o7LIUIDAyEQCDA5MmT1R1Kib1+/RqDBw+GnZ0dTE1NUbduXdy8eVPdYcktNzcXs2fPhpeXF0xMTFChQgUsWLAAYrFY3aF91rlz5+Dv7w9XV1cIBAL8888/MtuJCPPmzYOrqytMTEzQqlUr3L17Vz3BfsanXktOTg5mzJgBb29vmJmZwdXVFUOHDsWbN28Ucm6eBDRUaGgovvnmG1y5cgWnTp1Cbm4u2rdvj7S0NHWHVirXr1/Hb7/9htq1a6s7lBL78OEDmjVrBkNDQxw7dgz37t3DihUrYG1tre7Q5LZ06VJs2LABa9euxf3797Fs2TIsX74ca9asUXdon5WWloY6depg7dq1hW5ftmwZVq5cibVr1+L69etwdnZGu3btpANGapJPvZb09HSEhYVhzpw5CAsLw/79+/Ho0SN069ZNMScnTivExcURAAoNDVV3KCWWkpJClStXplOnTpGvry9NmjRJ3SGVyIwZM6h58+bqDkMhunTpQiNGjJBZ16tXLxo8eLCaIioZAHTgwAHp32KxmJydnWnJkiXSdZmZmWRlZUUbNmxQQ4TF9/FrKcy1a9cIAL148aLU5+N3AloiKSkJAGBra6vmSErum2++QZcuXdC2bVt1h1Iqhw4dQoMGDdC3b184OjqiXr16+P3339UdVok0b94cZ86cwaNHjwAAERERuHDhAjp37qzmyEonKioKsbGxMjMOCoVC+Pr6FjnjoDZJSkqCQCBQyN1nmRk7qCwjIkydOhXNmzdHrVq11B1OiQQFBSEsLAzXr19Xdyil9uzZM6xfvx5Tp07F999/j2vXrmHixIkQCoUYOnSousOTy4wZM5CUlIRq1apBX18fIpEIP/74IwYMGKDu0EpFMr9IYTMOvnjxQh0hKUxmZiZmzpyJgQMHKmSEVJ4EtMD48eNx+/ZtXLhwQd2hlEh0dDQmTZqEkydPwtjYWN3hlJpYLEaDBg2wePFiAEC9evVw9+5drF+/XuuSwO7du7F9+3bs3LkTNWvWxK1btzB58mS4urpi2LBh6g6v1OSZcVAb5OTk4Msvv4RYLMa6desUckyeBDTchAkTcOjQIZw7d04r5jwozM2bNxEXF4f69etL14lEIpw7dw5r165FVlYW9PX11RihfFxcXFCjRg2ZddWrV8e+ffvUFFHJffvtt5g5cya+/PJLAIC3tzdevHiBwMBArU4Czs7OANgdgYuLi3S9PDMOapqcnBz069cPUVFRCA4OVtg8CbxOQEMREcaPH4/9+/cjODgYXl5e6g6pxNq0aYPIyEjcunVLujRo0ACDBg3CrVu3tCoBAECzZs0KNNd99OiRVs5fkZ6eDj092cuAvr6+VjQR/RQvLy84Ozvj1KlT0nXZ2dkIDQ2Fj4+PGiMrGUkCePz4MU6fPg07OzuFHZvfCWiob775Bjt37sTBgwdhYWEhLeO0srKCiYmJmqOTj4WFRYG6DDMzM9jZ2WllHceUKVPg4+ODxYsXo1+/frh27Rp+++03/Pbbb+oOTW7+/v748ccf4e7ujpo1ayI8PBwrV67EiBEj1B3aZ6WmpuLJkyfSv6OionDr1i3Y2trC3d0dkydPxuLFi1G5cmVUrlwZixcvhqmpKQYOHKjGqAv3qdfi6uqKPn36ICwsDEeOHIFIJJJeD2xtbWFkZFS6k5e6fRGnFAAKXTZv3qzu0BRCm5uIEhEdPnyYatWqRUKhkKpVq0a//fabukMqkeTkZJo0aRK5u7uTsbExVahQgX744QfKyspSd2ifFRISUuh3ZNiwYUTEmokGBASQs7MzCYVCatmyJUVGRqo36CJ86rVERUUVeT0ICQkp9bn5fAIcx3E6jNcJcBzH6TCeBDiO43QYTwIcx3E6jCcBjuM4HcaTAMdxnA7jSYDjOE6H8STAcRynw3gS4DgtMW/ePNStW1fdYXBlDE8CHKeBCpsukeOUgScBjuM4HcaTAMd9QqtWrTBhwgRMnjwZNjY2cHJywm+//Ya0tDR89dVXsLCwQMWKFXHs2DHpc0JDQ9GoUSMIhUK4uLhg5syZyM3NlTnmxIkT8d1338HW1hbOzs6YN2+edLunpycAoGfPnhAIBNK/Jf766y94enrCysoKX375pUbOmctpD54EOO4ztm7dCnt7e1y7dg0TJkzA2LFj0bdvX/j4+CAsLAwdOnTAkCFDkJ6ejtevX6Nz585o2LAhIiIisH79evzxxx9YtGhRgWOamZnh6tWrWLZsGRYsWCAd9lgy+9rmzZsRExMjMxvb06dP8c8//+DIkSM4cuQIQkNDsWTJEtW9GVyZw4eS5gCwSV5ycnLUHYbGcXR0RLt27TB9+nQAbBjpnTt3onLlyhgyZAgAYNasWfj3339x+/ZthISEoFGjRvjpp5+kv+KXLl2KFStW4LvvvoOenh4cHR3Rtm1bzJgxAwBQvnx57NmzB9euXUOLFi1gYWEBDw8P2NvbS+eQzczMhLm5Odzc3PDbb7/BzMwMADBu3DjcuHEDmZmZqn9zOBgaGmrdfBgf46OI6jgiQmxsLBITE9UdikaKjY2FkZERbG1tpetevXoFCwsLWFlZSde9ePECDg4OSEtLg0AggL29vXRbdnY2YmJiUK5cORgYGBR6zLi4OOjr60snC5Ecz9TUVLpPYmIi0tPT4erqKl2XnJyMlJQUlCtXTimvn/s8a2trODs7a+20lfxOQMdJEoCjoyNMTU219oOsTEKhUObCm52dDTs7O5kLfVpaGtzc3JCYmAgDAwOZi3JGRgZycnLg4eEhnQDk42Pq6elBX19fOoWo5Hj5E01cXBySkpJkZpmLj4/H+/fvtXrmOW1FREhPT0dcXBwAyExjqU14EtBhIpFImgAUOV1dWaKnpwcDAwMYGxtL1wkEAhgaGsqsAwAjIyOYm5vjw4cPEAqF0oSanJwMPT09WFhYQCAQFHpMfX196OvrS9cJBAIYGRnJ7GNgYAA9PT2ZdYaGhhAIBAVi4VRDMstfXFwcHB0dtbJoiFcM6zBJHUD+IgeudBwcHJCdnY3o6GhkZGQgMTERb968kbu4wMjICMnJycjJyZFpWcRpHsn3R1vr1PidAMeLgBTIyMgIlStXxqtXr/Du3TsYGBjA3t5e7qKC8uXLIzo6GvHx8TA0NETt2rWVFDFXWtr+/eEVwzosMzMTUVFR8PLy4sUJHFdC2v494sVBHMdxOownAY7jOB3GkwDHqYinpyeGDx9eoufOmzevQNnzunXrsGXLllLFtHjx4kIHqjt79iwEAgHOnj1bquNzmo/XCegwbS/L1Dbh4eGwtLRExYoV5X7uq1ev8OrVKzRp0kS6rlatWrC3ty/Vhdrc3Bx9+vQpkEySk5Nx79491KhRA5aWliU+vi7Q9u8Rbx3EcYXIyMiQtgFXlHr16pX4uW5ubtKOZKpgaWkpk3C4sosXB3EFEBHS0tI0ZinpzaqkCCU8PBy9evWCpaUlrKysMHjwYLx79066n6enJ7p27Yr9+/ejXr16MDY2xvz58wGwHtWjR4+Gm5sbjIyM4OXlhfnz5xdou5+VlYUFCxagevXqMDY2hp2dHVq3bo1Lly7JnCd/cZCkyGX79u2YOnUqnJ2dYWJiAl9fX4SHhxf6WvIf6+7duwgNDYVAIJAZbTQzMxPTpk1D3bp1YWVlBVtbWzRt2hQHDx6UOaZAIEBaWhq2bt0qPUarVq1kYvv4LuPQoUNo2rQpTE1NYWFhgXbt2uHy5cuFxnr37l0MGDAAVlZWcHJywogRI5CUlPT5/zhOpfidAFdAeno6zM3N1R2GVGpqqnTAtJLo2bMn+vXrhzFjxuDu3buYM2cO7t27h6tXr8LQ0BAAEBYWhvv372P27Nnw8vKCmZkZYmNj0ahRI+jp6WHu3LmoWLEiLl++jEWLFuH58+fYvHkzACA3NxedOnXC+fPnMXnyZPj5+SE3NxdXrlzBy5cv4ePj88n4vv/+e3zxxRfYtGkTkpKSMG/ePLRq1Qrh4eGoUKFCoc85cOAA+vTpAysrK6xbtw4AG4oCYAkpISEB06dPR7ly5ZCdnY3Tp0+jV69e2Lx5M4YOHQoAuHz5Mvz8/NC6dWvMmTMHAD5Z9LNz504MGjQI7du3x65du5CVlYVly5ahVatWOHPmDJo3by6zf+/evdG/f3+MHDkSkZGRmDVrFgDgzz///OT7wakYcTorIyOD7t27RxkZGTLrU1NTCYDGLKmpqSV6fQEBAQSApkyZIrN+x44dBIC2b99OREQeHh6kr69PDx8+lNlv9OjRZG5uTi9evJBZ/9NPPxEAunv3LhERbdu2jQDQ77///sl4PDw8aNiwYdK/Q0JCCAB98cUXJBaLpeufP39OhoaG9PXXXxd4LfnVrFmTfH19P/0mEFFubi7l5OTQyJEjqV69ejLbzMzMZGL6OLaQkBAiIhKJROTq6kre3t4kEomk+6WkpJCjoyP5+PgUiHXZsmUyxxw3bhwZGxvLvNayoKjvkbbgdwJcAaampkhNTVV3GFKlHdZi0KBBMn/369cPw4YNQ0hIiHRb7dq1UaVKFZn9jhw5gtatW8PV1VWm+KdTp06YPn06QkNDUaNGDRw7dgzGxsYYMWJEieIbOHCgTFGPh4cHfHx8EBISUqLjAcCePXuwevVqREREIC0tTbq+pBWXDx8+xJs3bzB58mTo6eWVIpubm6N3797YuHEj0tPTZf6vunXrJnOM2rVrIzMzE3FxcXBycipRHJzi8STAFSAQCEpV/KJpnJ2dZf42MDCAnZ0d3r9/L11X2LAOb9++xeHDh6VFRh+Lj48HALx79w6urq4yF8fSxCdZFxERUaLj7d+/H/369UPfvn3x7bffwtnZGQYGBli/fn2Ji2Ik71Vh75OrqyvEYjE+fPggkwQ+HpRQUlyVkZFRohg45eBJgCvzYmNjZYZ2zs3Nxfv372UuUoWN/2Jvb4/atWvjxx9/LPS4kqGgHRwccOHCBYjF4hIlgtjY2ELXlXRk1+3bt8PLywu7d++WeV1ZWVklOh6Qd0GPiYkpsO3NmzfQ09ODjY1NiY/PqQ9vHcSVeTt27JD5+++//0Zubq60JUxRunbtijt37qBixYpo0KBBgUWSBDp16oTMzMwSd9zatWuXTAuoFy9e4NKlS5+NTygUFvqrWjIMdf4EEBsbW6B10KeO8bGqVauiXLly2Llzp0ysaWlp2Ldvn7TFEKd9+J0AV+bt378fBgYGaNeunbR1UJ06ddCvX79PPk8y76+Pjw8mTpyIqlWrIjMzE8+fP8e///6LDRs2wM3NDQMGDMDmzZsxZswYPHz4EK1bt4ZYLMbVq1dRvXp1fPnll588T1xcHHr27IlRo0YhKSkJAQEBMDY2lramKYq3tzeCgoKwe/duVKhQAcbGxvD29pY2dx03bhz69OmD6OhoLFy4EC4uLnj8+HGBY5w9exaHDx+Gi4sLLCwsULVq1QLn0tPTw7JlyzBo0CB07doVo0ePRlZWFpYvX47ExEQ+z7E2U3fNNKc+2t6q4XMkrVRu3rxJ/v7+ZG5uThYWFjRgwAB6+/atdD8PDw/q0qVLocd49+4dTZw4kby8vMjQ0JBsbW2pfv369MMPP8i0WsrIyKC5c+dS5cqVycjIiOzs7MjPz48uXbokc57CWgf99ddfNHHiRHJwcCChUEgtWrSgGzduFPpa8nv+/Dm1b9+eLCwsCAB5eHhIty1ZsoQ8PT1JKBRS9erV6ffffy/0GLdu3aJmzZqRqakpAZC2Nvq4dZDEP//8Q40bNyZjY2MyMzOjNm3a0MWLFwuN9d27dzLrN2/eTAAoKiqq0PdaW2n794gPG6HDtL27++fMmzcP8+fPx7t372SmgtQUZ8+eRevWrbFnzx706dNH3eFwJaTt3yNeJ8BxHKfDeBLgOI7TYbw4SIdp+20sx2kCbf8e8TsBjuM4HcaTAMdxnA7jSYDjOE6H8STAcRynw3gS4DiO02E8CXAcx+kwngQ4juN0GE8CnM4pbL7e/HP/KsqWLVsgEAhw48YNhR+b4xSFjyLK6bwDBw58cm5djivLeBLgdF69evXUHUKpfDytI8fJgxcHcQUQAWlpmrOUZmCTo0ePom7duhAKhfDy8sJPP/1UYJ+Pi4PEYjEWLVqEqlWrwsTEBNbW1qhduzZ+/vlnmec9ePAAAwYMgJOTE4RCIdzd3TF06NACM3ilpKRg7NixsLe3h52dHXr16oU3b97I7LN79260b98eLi4uMDExQfXq1TFz5kyZ+YEBYPjw4TA3N0dkZCTat28PCwsLtGnTBgCQmJiIkSNHwtbWFubm5ujSpQuePXsGgUCAefPmyRzn8ePHGDhwIBwdHSEUClG9enX8+uuv8r69XBnA7wS4AtLTAXNzdUeRJzUVKMmUx2fOnEH37t3RtGlTBAUFQSQSYdmyZXj79u0nn7ds2TLMmzcPs2fPRsuWLZGTk4MHDx4gMTFRuk9ERASaN28Oe3t7LFiwAJUrV0ZMTAwOHTqE7Oxs6Xy6APD111+jS5cu2LlzJ6Kjo/Htt99i8ODBCA4Olu7z+PFjdO7cGZMnT4aZmRkePHiApUuX4tq1azL7AUB2dja6deuG0aNHY+bMmcjNzYVYLIa/vz9u3LiBefPm4YsvvsDly5fRsWPHAq/v3r178PHxgbu7O1asWAFnZ2ecOHECEydORHx8PAICAuR/szntpd7pDDh1KmoyjNRUIvb7WzOWfHO3yKVx48bk6uoq8/qSk5PJ1tZWZnKVjyd76dq1K9WtW/eTx/bz8yNra2uKi4srch/JJCrjxo2TWb9s2TICQDExMYU+TywWU05ODoWGhhIAioiIkG4bNmwYAaA///xT5jlHjx4lALR+/XqZ9YGBgQSAAgICpOs6dOhAbm5ulJSUJLPv+PHjydjYmBISEj752jlZ2j6pDC8O4gowNWW/vjVlKUlxd1paGq5fv45evXrJjOxoYWEBf3//Tz63UaNGiIiIwLhx43DixAkkJyfLbE9PT0doaCj69esHBweHz8bSrVs3mb9r164NgM0lLPHs2TMMHDgQzs7O0NfXh6GhIXx9fQEA9+/fL3DM3r17y/wdGhoKAAWmzBwwYIDM35mZmThz5gx69uwJU1NT5ObmSpfOnTsjMzMTV65c+exr4soOXhzEFSAQlKz4RZN8+PABYrEYzs7OBbYVti6/WbNmwczMDNu3b8eGDRugr6+Pli1bYunSpWjQoAE+fPgAkUgENze3YsViZ2cn87ekqEgywXtqaipatGgBY2NjLFq0CFWqVIGpqSmio6PRq1evAhPBm5qaFmjN9P79exgYGMDW1lZmvZOTU4H9cnNzsWbNGqxZs6bQeOPj44v1uriygScBrkyysbGBQCBAbGxsgW2FrcvPwMAAU6dOxdSpU5GYmIjTp0/j+++/R4cOHRAdHQ1bW1vo6+vj1atXCok1ODgYb968wdmzZ6W//gHI1EHkl7+Pg4SdnR1yc3ORkJAgkwg+fq02NjbQ19fHkCFD8M033xR6fC8vrxK8Ck5b8eIgrkwyMzNDo0aNsH//fmRmZkrXp6Sk4PDhw8U+jrW1Nfr06YNvvvkGCQkJeP78OUxMTODr64s9e/Yo5Fez5KKevzIZADZu3FjsY0iSx+7du2XWBwUFyfxtamqK1q1bIzw8HLVr10aDBg0KLB/fuXBlG78T4MqshQsXomPHjmjXrh2mTZsGkUiEpUuXwszMDAkJCUU+z9/fH7Vq1UKDBg3g4OCAFy9eYPXq1fDw8EDlypUBACtXrkTz5s3RuHFjzJw5E5UqVcLbt29x6NAhbNy4ERYWFsWO08fHBzY2NhgzZgwCAgJgaGiIHTt2ICIiotjH6NixI5o1a4Zp06YhOTkZ9evXx+XLl7Ft2zYAgJ5e3u+9n3/+Gc2bN0eLFi0wduxYeHp6IiUlBU+ePMHhw4cLtEbiyjZ+J8CVWe3atcM///yD5ORk9O/fH1OnTkXv3r0xYsSITz6vdevWOHfuHMaMGYN27dph9uzZaNOmDUJDQ2FoaAgAqFOnDq5du4b69etj1qxZ6NixI2bMmAGhUAgjIyO54rSzs8PRo0dhamqKwYMHY8SIETA3Ny/wq/5T9PT0cPjwYXz55ZdYsmQJunfvjvPnz2P79u0A2B2NRI0aNRAWFoZatWph9uzZaN++PUaOHIm9e/dK+xxwuoPPMazDtH1uVO7zdu7ciUGDBuHixYvw8fFRdzhlkrZ/j3hxEMeVEbt27cLr16/h7e0NPT09XLlyBcuXL0fLli15AuCKxJMAx5URFhYWCAoKwqJFi5CWlgYXFxcMHz4cixYtUndonAbjSYDjyoiuXbuia9eu6g6D0zK8YpjjOE6H8STAcRynw3gS4DiO02E8CXAcx+kwngQ4juN0GE8CHMdxOownAY7jOB3GkwCns7Zs2QKBQIAbN26oOxSlKmyOYWVbt24dtmzZotJzciXDO4txXBl3+fLlYk+Aoyjr1q2Dvb09hg8frtLzcvLjSYDjyrgmTZqoOwSFEIlEyM3NLTDvAlc6vDiIK7MePHiAAQMGwMnJCUKhEO7u7hg6dCiysrJk9ktJScHYsWNhb28POzs79OrVC2/evJHZZ/fu3Wjfvj1cXFxgYmKC6tWrY+bMmUhLS5PZb/jw4TA3N8eTJ0/QuXNnmJubo3z58pg2bVqB87569Qp9+vSBhYUFrK2tMWjQIFy/fh0CgaBAUcqNGzfQrVs32NrawtjYGPXq1cPff/9drPfh4+IgSTFYSEjIZ193YZ49e4Yvv/wSrq6uEAqFcHJyQps2bXDr1i0AgKenJ+7evYvQ0FAIBAIIBAJ4enpKn//y5UsMHjwYjo6OEAqFqF69OlasWAGxWCzd5/nz5xAIBFi2bBkWLVoELy8vCIVChISEIDMzE9OmTUPdunVhZWUFW1tbNG3aFAcPHiwQa2JiIkaOHAlbW1uYm5ujS5cuePbsWaFFZI8fP8bAgQNl4vr111+L9R5rM34nwBVEBKSnqzuKPKambOJjOURERKB58+awt7fHggULULlyZcTExODQoUPIzs6W+TX59ddfo0uXLti5cyeio6Px7bffYvDgwTKTqzx+/BidO3fG5MmTYWZmhgcPHmDp0qW4du1agUlYcnJy0K1bN4wcORLTpk3DuXPnsHDhQlhZWWHu3LkAgLS0NLRu3RoJCQlYunQpKlWqhOPHj6N///4FXktISAg6duyIxo0bY8OGDbCyskJQUBD69++P9PT0Ehe5FOd1F6Zz584QiURYtmwZ3N3dER8fj0uXLkmnwzxw4AD69OkDKysrrFu3DkDerGnv3r2Dj48PsrOzsXDhQnh6euLIkSOYPn06nj59Kt1f4pdffkGVKlXw008/wdLSEpUrV0ZWVhYSEhIwffp0lCtXDtnZ2Th9+jR69eqFzZs3Y+jQoQAAsVgMf39/3LhxA/PmzcMXX3yBy5cvo2PHjgVe07179+Dj4wN3d3esWLECzs7OOHHiBCZOnIj4+HgEBASU6D3WCsTprIyMDLp37x5lZGTIbkhNJWKpQDOW1FS5X5ufnx9ZW1tTXFxckfts3ryZANC4ceNk1i9btowAUExMTKHPE4vFlJOTQ6GhoQSAIiIipNuGDRtGAOjvv/+WeU7nzp2patWq0r9//fVXAkDHjh2T2W/06NEEgDZv3ixdV61aNapXrx7l5OTI7Nu1a1dycXEhkUhU5GskIgJAAQEBpX7dRETx8fEEgFavXv3Jc9asWZN8fX0LrJ85cyYBoKtXr8qsHzt2LAkEAnr48CEREUVFRREAqlixImVnZ3/yXLm5uZSTk0MjR46kevXqSdcfPXqUAND69etl9g8MDCzwnnTo0IHc3NwoKSlJZt/x48eTsbExJSQkFHn+Ir9HWoIXB3FlTnp6OkJDQ9GvXz84ODh8dv9u3brJ/F27dm0AwIsXL6Trnj17hoEDB8LZ2Rn6+vowNDSUzut7//59mecLBAL4+/sXOGb+44WGhsLCwqLAr9IBAwbI/P3kyRM8ePAAgwYNAgDk5uZKl86dOyMmJgYPHz787GssTHFe98dsbW1RsWJFLF++HCtXrkR4eLhMMc7nBAcHo0aNGmjUqJHM+uHDh4OICtyFdOvWTTqbW3579uxBs2bNYG5uDgMDAxgaGuKPP/6Q+b8IDQ0FAPTr10/muR+/x5mZmThz5gx69uwJU1PTAu9xZmYmrly5UuzXqG14cRBXkKkpkJqq7ijymJrKtfuHDx8gEomK3SLm44nVJUUXGRkZAIDU1FS0aNECxsbGWLRoEapUqQJTU1NER0ejV69e0v3ywjUtMMOUUCiUmfD+/fv3cHJyKhDLx+vevn0LAJg+fTqmT59eaPwlnez+c6+7MAKBAGfOnMGCBQuwbNkyTJs2Dba2thg0aBB+/PHHz86t/P79e5n6AQlXV1fp9vxcXFwK7Lt//37069cPffv2xbfffgtnZ2cYGBhg/fr1+PPPP2XOZWBgAFtbW5nnf/wev3//Hrm5uVizZg3WrFlTaNwlfY+1AU8CXEECAWBmpu4oSszW1hb6+vp49eqVQo4XHByMN2/e4OzZs9Jf/wCkZeAlYWdnh2vXrhVYHxsbK/O3vb09AGDWrFno1atXoceqWrVqieMoCQ8PD/zxxx8AgEePHuHvv//GvHnzkJ2djQ0bNnzyuXZ2doiJiSmwXlIhLXm9EoJC6oK2b98OLy8v7N69W2b7xxXvdnZ2yM3NRUJCgkwi+Pg9trGxgb6+PoYMGYJvvvmm0Li9vLw++bq0GS8O4socExMT+Pr6Ys+ePQr5BSe50HzcNHHjxo0lPqavry9SUlJw7NgxmfVBQUEyf1etWhWVK1dGREQEGjRoUOjyuV/fylSlShXMnj0b3t7eCAsLk64XCoWF3lG0adMG9+7dk9kXALZt2waBQIDWrVt/9pwCgQBGRkYyCSA2NrZA6yBJwt69e7fM+o/fY1NTU7Ru3Rrh4eGoXbt2oe/xx3dNZQm/E+DKpJUrV6J58+Zo3LgxZs6ciUqVKuHt27c4dOgQNm7cKNeF08fHBzY2NhgzZgwCAgJgaGiIHTt2ICIiosTxDRs2DKtWrcLgwYOxaNEiVKpUCceOHcOJEycAAHp6eb/PNm7ciE6dOqFDhw4YPnw4ypUrh4SEBNy/fx9hYWHYs2dPieOQ1+3btzF+/Hj07dsXlStXhpGREYKDg3H79m3MnDlTup+3tzeCgoKwe/duVKhQAcbGxvD29saUKVOwbds2dOnSBQsWLICHhweOHj2KdevWYezYsahSpcpnY+jatSv279+PcePGoU+fPoiOjsbChQvh4uKCx48fS/fr2LEjmjVrhmnTpiE5ORn169fH5cuXsW3bNgCy7/HPP/+M5s2bo0WLFhg7diw8PT2RkpKCJ0+e4PDhw59tMaXV1F0zzamPtrdq+Jx79+5R3759yc7OjoyMjMjd3Z2GDx9OmZmZRJTXSub69esyzwsJCSEAFBISIl136dIlatq0KZmampKDgwN9/fXXFBYWVqAlz7Bhw8jMzKxALAEBAfTx1+3ly5fUq1cvMjc3JwsLC+rduzf9+++/BIAOHjwos29ERAT169ePHB0dydDQkJydncnPz482bNjw2fcBRbQOKs7r/tjbt29p+PDhVK1aNTIzMyNzc3OqXbs2rVq1inJzc6X7PX/+nNq3b08WFhYEgDw8PKTbXrx4QQMHDiQ7OzsyNDSkqlWr0vLly2VaOUlaBy1fvrzQOJYsWUKenp4kFAqpevXq9Pvvvxf6HickJNBXX31F1tbWZGpqSu3ataMrV64QAPr5559l9o2KiqIRI0ZQuXLlyNDQkBwcHMjHx4cWLVpU5PtBpP3fIwERkXrSD6dumZmZiIqKgpeXV4GKTE49Fi9ejNmzZ+Ply5cqH+pBV+zcuRODBg3CxYsX4ePjU+rjafv3iBcHcZyarF27FgBQrVo15OTkIDg4GL/88gsGDx7ME4CC7Nq1C69fv4a3tzf09PRw5coVLF++HC1btlRIAigLeBLgODUxNTXFqlWr8Pz5c2RlZcHd3R0zZszA7Nmz1R1amWFhYYGgoCAsWrQIaWlpcHFxwfDhw7Fo0SJ1h6YxeHGQDtP221iO0wTa/j3iTUQ5juN0GE8CHMdxOownAQ68RJDjSk7bvz88CegwycBc6Zo0bDTHaRnJ96ewge60AW8dpMP09fVhbW2NuLg4AKy1SmFjtXAcVxARIT09HXFxcbC2toa+vr66QyoR3jpIxxERYmNjSzUYGsfpMmtrazg7O2vtDyieBDgAbP7WnJwcdYfBcVrF0NBQa+8AJHgS4DiO02G8YpjjOE6H8STAcRynw3gS4DiO02E8CXAcx+kwngQ4juN0GE8CHMdxOownAY7jOB32f7ewFEhgL3bHAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# As an example, I'll set up a plot where I ignore ET and calculate delta_S based only on \n", "# precipitation and discharge. Your job is to add in ET to the calculation and the plot. \n", "# See how your answer differs if you assume ET is equal to either of the pan evaporation measurements.\n", "\n", "P_mm = 25.4*data.Avg_total_P_in\n", "# The above is total precip per month, on average, in mm\n", "\n", "Q_mm_per_sec = 1000*(0.3048)**3/(217*1000*1000)*data.Avg_streamflow_cfs\n", "# check my unit conversions here -- I converted km to m and feet to meter and then multiplied by 1000 to get mm\n", "\n", "# Now we also need mm per month. I will assume 30 days per month to scale up. You may want to be more precise.\n", "Q_mm_per_month = Q_mm_per_sec*3600*24*30\n", "\n", "Delta_S = P_mm - Q_mm_per_month\n", "\n", "fig, ax = plt.subplots(figsize=(4,4))\n", "\n", "ax.plot(data.month,P_mm,c='k',label='precipitation');\n", "ax.plot(data.month,Q_mm_per_month,c='b',label='discharge');\n", "ax.plot(data.month,Delta_S,c='r',label='change in storage');\n", "\n", "# Scatterplot\n", "\n", "ax.set_xlabel('month')\n", "ax.set_ylabel('water input/output (mm per month)');\n", "plt.legend(loc=(0.2,-0.36), fontsize=12);\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Think about your plot**: Does the seasonality of the values make sense?\n", "\n", "In the western U.S., changes in storage in a mountain watershed over a year tend to positive sometimes and negative other times, but on average, sum to 0 over the year. Is that true here?\n", "\n", "Calculate the annual sums of each of the values in your monthly water balance above." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2568.448\n" ] } ], "source": [ "# I'll sum the precipitation, and then you do the rest. Print out the values.\n", "P_annual = np.sum(P_mm)\n", "print(P_annual)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Now sum all of those sums to see what that annual residual is for the watershed. Again, we expect this to be close to 0." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "If your value is not close to 0, think about what else might be happening here. You may want to look a the web links with historical information about the watershed or look at your lecture notes. You may want to discuss with another classmate." ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.13" } }, "nbformat": 4, "nbformat_minor": 4 }