{ "cells": [ { "cell_type": "markdown", "id": "3979c593", "metadata": {}, "source": [ "# Article for Journal of Open Source Software\n", "Martin Vonk & Aaron Peche (2026)\n", "\n", "This notebook replicates the results presented in the article submitted to the Journal of Open Source Software (JOSS). The article can be found here: https://joss.theoj.org/papers/095927c33334e8ebb4612243ea40fdae doi.org/...\n", "\n", "JOSS is a developer-friendly, open-access academic journal (ISSN 2475-9066) dedicated to research software packages and features a formal peer-review process. The pre-review and review of the Pedon package are publicly available in issues https://github.com/openjournals/joss-reviews/issues/9769 and https://github.com/openjournals/joss-reviews/issues/10409, respectively.\n", "\n", "## Setup" ] }, { "cell_type": "code", "execution_count": 1, "id": "906c9e13", "metadata": {}, "outputs": [], "source": [ "from pathlib import Path\n", "\n", "import matplotlib as mpl\n", "import matplotlib.pyplot as plt\n", "import matplotlib.ticker as mplticker\n", "import numpy as np\n", "from cycler import cycler\n", "\n", "import pedon as pe\n", "\n", "rc_params = {\n", " \"axes.prop_cycle\": cycler(\n", " color=[\n", " \"#3f90da\",\n", " \"#ffa90e\",\n", " \"#bd1f01\",\n", " \"#94a4a2\",\n", " \"#832db6\",\n", " \"#a96b59\",\n", " \"#e76300\",\n", " \"#b9ac70\",\n", " \"#717581\",\n", " \"#92dadd\",\n", " ]\n", " ),\n", " \"axes.titlesize\": 10.0,\n", " \"axes.grid\": True,\n", " \"axes.labelsize\": 9.0,\n", " \"xtick.labelsize\": 8.0,\n", " \"ytick.labelsize\": 8.0,\n", " \"legend.fontsize\": 8.0,\n", " \"legend.framealpha\": 1.0,\n", " \"figure.figsize\": [7.0, 6.0],\n", "}\n", "\n", "mpl.rcParams.update(rc_params)\n", "fig_path = Path.cwd().parent.parent / \"paper/figures\"" ] }, { "cell_type": "markdown", "id": "adce1604", "metadata": {}, "source": [ "## Defining a soil model" ] }, { "cell_type": "code", "execution_count": 2, "id": "cc2b7f28", "metadata": {}, "outputs": [], "source": [ "# Mualem-van Genuchten parameters for Sandy Loam\n", "mg = pe.Genuchten(\n", " k_s=106.1, # saturated conductivity (cm/d)\n", " theta_r=0.065, # residual water content (-)\n", " theta_s=0.41, # saturated water content (-)\n", " alpha=0.075, # shape parameter (1/cm)\n", " n=1.89, # shape parameter (-)\n", ")\n", "\n", "h = np.logspace(-1, 6, 8) # pressure head (cm)\n", "theta = mg.theta(h) # water content (-) at pressure head values\n", "k = mg.k(h) # hydraulic conductivity (cm/d) at pressure head values" ] }, { "cell_type": "code", "execution_count": 3, "id": "c8e10c3c", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array([[1.00000000e-01, 4.09984348e-01, 1.03389137e+02],\n", " [1.00000000e+00, 4.08791540e-01, 8.59093947e+01],\n", " [1.00000000e+01, 3.43096726e-01, 1.34676049e+01],\n", " [1.00000000e+02, 1.21823289e-01, 4.55156715e-03],\n", " [1.00000000e+03, 7.23953055e-02, 2.81336545e-07],\n", " [1.00000000e+04, 6.59528265e-02, 1.67662224e-11],\n", " [1.00000000e+05, 6.51227480e-02, 9.98706595e-16],\n", " [1.00000000e+06, 6.50158130e-02, 5.94891489e-20]])" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "np.array(list(zip(h, theta, k)))" ] }, { "cell_type": "markdown", "id": "d4fce495", "metadata": {}, "source": [ "## Loading soil models from dataset" ] }, { "cell_type": "code", "execution_count": 4, "id": "952f797c", "metadata": {}, "outputs": [], "source": [ "hydrus = pe.Soil(\"Sand\").from_name(pe.Genuchten, source=\"HYDRUS\")\n", "staring = pe.Soil(\"B05\").from_name(pe.Genuchten, source=\"Staring_2018\")\n", "rawls = pe.Soil(\"Sand\").from_name(pe.Brooks, source=\"Rawls\")\n", "clapp = pe.Soil(\"Sand\").from_name(pe.Campbell, source=\"Clapp\")" ] }, { "cell_type": "code", "execution_count": 5, "id": "b258ae56", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPkAAAFECAYAAADlbJiPAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjksIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvJkbTWQAAAAlwSFlzAAAPYQAAD2EBqD+naQAAZJpJREFUeJztnXd4VGXah+8zJZM+CWmkh5ACSSAh9N4FCysWVtQVYsWFXcSy6rquYsfdFWyoLLuCruKnggUFC0hHRJDQWwKpkELqpGcm835/TDIkpJBJJskkOfd1zRVOe89zDvOb93nb80hCCIGMjEyPRdHVBsjIyHQssshlZHo4sshlZHo4sshlZHo4sshlZHo4sshlZHo4sshlZHo4sshlZHo4sshlZHo4sshlZHo4sshlZHo4PUrkq1atYurUqUyaNAm9Xt/V5sjI2ASqrjbAWqSnp3Ps2DF++umnrjZFRsam6DE1+Y8//khZWRlTpkxh6dKlXW2OjIzNYJM1+eLFi9m4cSNpaWkkJiYSFxdnPpaUlMT8+fPJy8tDq9Wydu1aoqOjyc3NRZIktm3bxvz58zl8+HCD6+qoqqqiqqrKvG00GikoKMDDwwNJkjrh6WRkrIMQgpKSEvz8/FAoWqivhQ2yc+dOkZGRIYKDg0ViYmKDY5MnTxZr1qwRQgjx+eefi2HDhgkhhHj77bfF2rVrhRBC/Oc//xH/93//12TZzz77rADkj/zpMZ+MjIwW9SQJYbtBI0JCQvjqq6/MNXJubi5hYWEUFBSgUqkQQuDr68uePXvQ6XSsW7eOf/3rXyxZsoTbb7+dkSNHNiqzriZfvXo1q1evxmAwcO7cObbv3sXBU6fwUuZwk+KfGEZ/hPAcy5Hbp1J2+ggD3/4U97FTm7QzMaOCpT9cIqSPijdu9rPa8+v1erZv387kyZNRq9VWK9cayLa1DWvaVlJSQr9+/SgqKkKr1TZ7nk26682RkZGBr68vKpXJbEmSCAoKIj09nSlTpvDhhx8yadIkIiIimhQ4gEajQaPR8Oijj/Loo4+i0+nQarV4eHji4OiAWumKq0ICN1fw8MBFowaFhJuLMx4eHk2WqS0tQ2VfgcbRrtlz2oJer8fR0REPDw+b/LLKtlmONW2ru/5qzcxuJfKr8frrr7f63JUrV7Jy5UpqamoAUNf+cOhF3YuvdXDq2jq26/DIyLRItxJ5YGAgWVlZGAwGs7uenp5OUFCQxWUtWrSIRYsWmWvyxiI3IUm1Ijca221/b6GmpqZB52Zno9frUavVlJeX22RNbqltGo0GpVLZ5nt2K5F7e3sTHx/PRx99REJCAhs2bCAgIICwsDCLy7qyJlfVvkSzyOtq7lpXSMgibxU6nY7k5GS6uqvHx8eH5OTkLrWhOSy1TZIkwsLCcHV1bdP9bFLkCxYsYNOmTWRnZzNjxgxcXFzML2XVqlUkJCTw8ssv4+rqypo1a9p0jytrcrvaX1UjSgxCharWXZdkd73V1NTUkJycjIuLC76+vi0P68i0CqPRSFZWFsnJycTGxrapRrdJka9atarZY5GRkezbt8/q96zvOlXhYBZ5XZtcCLkmvxpVVVXmEQ9nZ+euNqfH4Ovri06no6qqCkdHR4uv77U/tStXriQqKorhw4cDJpdIY2cHQDX1XmRdz6XsrrcauQa3Lu19n732f2PRokWcPHmSAwcOmPfVuexVwsHcuW7ueJPddZluSq8VeVPU1eRVOHDlEJrc8SbTXem1Ir/SXQeucNdrO97k3vVuTUhICIcPH26wb9KkSXz22WdERERQXl5u3v/ZZ58xZMgQ9Ho9kyZNol+/fsTFxREWFsb06dPZtGmT+dwdO3bg4OBAXFwcgwcPZuTIkfzyyy/m4wkJCY3mbSxdupQlS5YAUFRUxB/+8AdiYmIYPHgwMTExrFu3zurPDzba8dYZXNm7DmCnrq3JhcPlE83tIdldtwQhBJWGjn1n9iqpzYuK7OzsmDZtGk8++SRvvvkmubm5PPLII3z//ffmTtgVK1Ywe/ZswCTquXPn8s4773DzzTcDpk7guh+Qt99+m3vuuYeTJ0+26v5PP/00Xl5eHDt2DEmSKCkpITs7u03PcjV6rcibooG7LupqcnkyTFuoNAiufedch97ju4X9cVC3feXgP//5T4YMGcLOnTt54403WLx4MTExMU2eO2nSJJYuXcorr7xiFnl9pk6dyhNPPNHqe2dmZjJmzBjzj5SLiwsuLi5te5Cr0GtFfuVkGKjvrtdrk8vuerfntttuw8HhsndWN+fCycmJNWvW8Lvf/Y6BAwfy2GOPtVjOyJEjeeSRR5o8tn79eubOndtqmx566CFuvfVWPv30U0aPHs3MmTO54YYbWn29JfRakTflrl+uyeu1yeXJMG3CXiXx3cL+HX6P1vDpp582iC0wadIk87/Hjh1LVFQUS5YsuepQ1ZWz+M6cOUNcXBzZ2dkYDAb2799vPtZcM6Ju/+TJk0lPT2fnzp38/PPPLFiwgNmzZ7Ny5cpWPZMl9FqRN0XDIbQra/Ka5i6TaQJJktrlSncmSqWyVTPJDhw40MCdr2uT6/V6Fi5cyJ133sm+ffuQJAkvLy/y8/MbXJ+Xl4e/v79528nJieuuu47rrruOG264gWuuuaZDRN5re9eboqG7bkKq+8+X3fVeze7du1m6dGmT7W61Ws0bb7xBZmYmX331FQAzZszg888/p6CgAICsrCy+/vprpk+fDpjClRUWFprL+O233+jfv2M8H7kmr0edyCvru+tK0ysSNXJN3tt4+OGHWbp0KWVlZQQHB7N69epm282Ojo689NJLLF26lNmzZzN16lQWL17M5MmTkSTTKMBLL73EsGHDADh27BiPPvooQggUCgW+vr589NFHHfIcvVbkTXe8aQCoFE6YRa6qE7mh022UaT+pqamN9u3YsaPF7eb21WfSpEmNxt/nz5/P/PnzzdsLFy5k4cKFTV5fF7SkM+i17npT01od7GtFjvPlIbS6mtwgi1yme9JrRd4UDhp7wNQmrzHWdrzVtsnlmlymuyKLvB4aOzskTB1slXrTX7lNLtPdkUVeD0mSsFeYwhZV1IlcJdfkMt0bWeRX4FAn8mqTyBVym1ymmyOL/ArqavJKvdwml+kZ9FqRN7XUFOrV5Fe2yeWavNvyxRdfMHToUOLi4hgwYABTpkzBaDTy+uuvt2vlV1xcHCUlJVa0tGPotePkTc1dB+q1ya8cJ5c73ixCCKgpv/p57UHpeDk8VzNkZWXxwAMP8NtvvxEcHAzAoUOHkCSJ119/nUmTJtG3b1+LblsXEvzKcXJbpdeKvDkcFNXAZXf9cu+6XJNbRE05fGG9lFFNcvNFUDm1eEpOTg5KpZI+ffqY98XHx/P8889z8eJF8wq1tWvXkp+fz9NPP01lZSXV1dU88sgj3HvvvYApCIRCoSA5OZnc3FxOnz6NJEkUFhbi5uZGSEgI8+bNY8uWLWRnZ3Pvvffy9NNPA3D69GnuvvtudDodkZGRlJaWcscdd5CQkNBhr6Y+ssivwEHZTE0uu+vdksGDBzNu3DiCg4OZOHEiY8aM4Y477uCZZ57h/fffb7BCrbCwkD179qBUKikoKGDIkCHMmDGDgIAAwDS/fM+ePc2u+y4qKmLfvn3k5eXRv39/7r77bvz9/bnrrrtYuHAhd999N6dOnWLIkCHccccdnfUKZJFfSV2bvLy6bqmp3PHWJpSOppq2o+9xFRQKBRs2bOD06dPs3LmT7777jpdeeomDBw82Ojc/P597772Xs2fPolKpyM/P5/jx42aRz5kzp8XADnXC9fT0JDQ0lJSUFFxcXDh8+DDz5s0DYODAgYwbN64tT9tmZJFfgZOyEqgncrlN3jYk6aqudGcyYMAABgwYwIIFC5g5cyYbN25sdM6DDz7Iddddx4YNG5Akifj4eCorK83HrxZL3t7e3vxvpVKJoRnvr60hq9pKr+1dbw4nlakmL6sWCCHkNnk358KFC+zdu9e8XVhYSEpKCv3798fV1ZXi4uIGx4KDg5EkiV27dnHkyJF239/V1ZXY2FjzCrMzZ86wZ8+edpdrCT2mJk9NTWX06NFERkYSFBTEhx9+2KZyHJV6AAxGqNbr67XJ5Zq8O2IwGHj++edJSUnB0dERg8HA/PnzufHGG7l06RL3338/jo6OrF27lmXLlrFw4UJeeOEF4uLimk1/bSkffvgh99xzD//85z8JCwtj+PDhuLm5WaXs1tBjRA5w/fXX85///KddZaiVYEc51ThSXlFhDhoh1+Tdk+DgYH744Ycmj913333cd999DfYlJSU1ee7atWsb7asfDurKJa312/xBQUHmiDEpKSmMHj2aoUOHtvIJ2k+Pctd/+OEHxo8fz8cff9z2QiQFTphcuLKKcnkyjEy7+fnnn4mLiyM2NpZZs2axYsUKAgMDO+3+NlmTL168mI0bN5KWlkZiYmKDIHxJSUnMnz+fvLw8tFota9euJTo6Gl9fX86cOYMkScyYMYOZM2fi4eHRqOyqqqoGubN1Oh1gyhut1+tRCgknqZhC4YuutBSn2pDMNdVV6PX6Ju011LryQohmz2kLdWVZs0xr0ZRttminLXDNNddwzTXXtLucuu9o/e3WYJMiv/XWW3n88cebHGpYsGABDzzwAAkJCaxfv56EhAQOHDiARqMxnzN+/HjOnTvXpMhfeeUVnnvuuUb7t2/fjqOjI8P1l3DEE4BDiYm4nzmDA5CTmUHq5s1N2ptW4QQEU6IrYXMz57SHLVu2WL1Ma1HfNrVajY+PTxda07PZu3dvA2HXz/7SEpLo6mzxLRASEsJXX31lrslzc3MJCwujoKAAlUplTpO7Z88e+vbti7OzM0IIZs6cydq1a/H19W1UZl1Nvnr1alavXm3OqZ2VlYWHhwfKX+/ll3RINE4nJjyCsPQznHksAbcxUxj0wXdN2nkwvYK/fptDqIeaVbf5N3lOW9Dr9WzZsoXp06c3SK1sCzRlW3l5OcnJyQwcOLBNKXZlmqa8vJxTp04RFhbW4L3qdDo8PT0pLi7G1dW12ettsiZvjoyMDHx9fVHV9nhLkkRQUBDp6emcP3+ep556CrVaza233tqkwAE0Gg0ajcYcY6tu7rparTZ9WRUqnDCF0q2sqkJd+1KFQd+s0FSqarM9HSFGs202SH3bbNXGnsKV34PWvu9uJfKWsLTd01QgRwAkJU5SEQCl5eVItfnRRL12vIxMd6Jb9a4HBgaSlZVlnkkkhCA9PZ2goCCLy2oqkCMAkgIXTLGyS8pKUdRGcDXqq9tnvIxMF9GtRO7t7U18fLx59tCGDRsICAggLCzMejeRlLhIJpGXVVQgamtyY7Vck3dHQkJCiIyMJC4ujsjISJYtW2aVcuunIbZ1bNJdX7BgAZs2bSI7O5sZM2bg4uJiTlK3atUqEhISePnll3F1dWXNmjVtukfz7roCB0pRSQKDkKisHUKTRW4ZQgiMFR27nlzh4NiqeeB1K80uXLhAVFQUU6ZMYcSIER1qmy1hkyJftWpVs8ciIyPZt29fu+/RXNAIJAWSBC52NRRWqSinNheaLHKLMFaUszOqY1Lx1jHxZAlKx9YvgvH392fAgAGkpaVRUlLS5NrxsrIyAgMDycnJQa1WM2LECMLCwli3bh3p6elMnjyZc+capmT+5ZdfWLRoETU1NRgMBhYtWsQf//hHaz9um7FJkXcGLXW8AbjYGSisUlFWG39drsm7P6dPnyY/P59JkyahUqmaXTseExPDvn37GDRoENXV1Rw4cAAhBFu2bGHq1KmNyn3llVd47LHHuP322wEa5DizBXqtyJuvyetErgfsKaudzSaL3DIUDo5MPNmx8c8UDq0bi7/ttttQKBScOXOGFStW4OXlRXJycrNrx6dNm8bWrVu5dOkS11xzDadOneLYsWNs3bqV2bNnNyp/8uTJvPDCCyQlJTFlypROXy9+NXqtyJulTuRqUw9+md6AGlnkliJJkkWudEdS1ybfunUrs2bNYsqUKTz88MPNrh2fNm0ajz32GJcuXeKmm27C39+fLVu2sH37dt58881G5S9ZsoQbb7yRrVu38tRTTxETE8M777zT2Y/ZLL1W5M2667UDDi5q05BZaXU17oCorjatL+/kBf8y1mPatGn88Y9/5Omnn25x7fiIESM4c+YMubm5LF++HH9/f66//nr69u2Ll5dXo3LPnDlDZGQk999/P4GBgTz11FOd+VhXpdeKvKWONwAXtanmLqmowL32kNDrkWrTG8t0T/7+978TFhbG448/zpNPPtnk2nGVSsXYsWMpLS3FwcGB6Oho9Hp9k+1xgLfffptt27ZhZ2eHUqnktdde66zHaRW9VuTNUuuua9Um162sshKjSo3CoMdYVYlCFnm34sp13u7u7uTnm6YtP/HEE81ed2V4qAsXLjTYXrp0qfnfb731VvuM7GC61WSYTkFhErGDVIld7dxgvYc3ADUVZV1mloxMW+m1Im8ug0qdyCVRjZuLaWWPwS8EgJpyWeQy3Y9eK/Jm564ralf2GA1oXU2TOfR9TSF5ZZHLdEd6rcibpbYmx3i5Jq/2NC1bNcruukw3RBb5ldQXee1C/Cp307CJXJPLdEd6rcibb5PXuet6c01epTXl0ZJFLtMd6bUib7ZNrqxfk9e2ye0dqdE4yL3r3RSDwcBzzz3HgAEDiImJIS4ujgceeKBBaLGejDxOfiX13HU7tR2ODg6UV1RQ7e0n1+QWIITA0MGx6lVKVatmIN57770UFBSwb98+3N3dEUKwfv16CgoKOtQ+W0EW+ZXUc9cBPNzcKK+ooNInQBa5BRhqDPz70//r0Hs8cNtc1KqW45wlJyfz+eefk56ejru7ae6iJEnMmTOHHTt2mM8zGAxcf/315OfnU1FRQWxsLKtXr8bJyYkdO3bwpz/9ifj4eA4dOoRGo+G///0vcXFxLR6zFXqtu94s9WpygD5aNwCqfAJkd70bcujQIcLDw/H09GzxPKVSybp16zh48CDHjx9Hq9U2mMl24sQJ5s+fz/Hjx3niiSeYO3euOYNKS8dsAYtr8rKyMo4fP05hYSHu7u7ExMTg5GQbq42sgrkmN4ncw83061/lE0BNWWlXWdXtUClVPHDb3A6/h7UQQrBixQo2bdqEwWCguLiYMWPGmI+HhISY567//ve/54EHHiAjI6PFY22JPdgRtPotffPNN7z11lvs2LEDe3t7XF1d0el0VFVVMWHCBP70pz9x4403dqStnYO5Jr/srgNU9Q3AkH6ii4zqfkiSdFVXujOIj48nKSmJ/Pz8JpNt1LFu3Tq2bdvGzp07cXV15c0332Tbtm3Nni9JUrP9AS0d6wpa5a6PHz+e119/nTvvvJPz58+j0+nIzMxEp9Nx7tw57rrrLt566y0mTJjQ0fZajatNa8VoWqDirtWCENQ4uVJRXtHJVsq0l7CwMG655RbuvfdeioqKAFOtvWHDBs6fP28+r7CwEE9PT1xdXSkpKWmU4DA1NZXt27cDsH79enx8fAgICLjqMVugVTX5q6++2sB1qU9AQADz5s1j3rx5Vom91lk0u9RU6WD6W2MSuVqlwlmCUkCH7fw6y7Se999/nxdffJGRI0eiUqkwGo1MmDCBa6+91nzOvHnz+Prrr4mMjMTLy4vx48eTlpZmPh4dHc3atWtZvHgxdnZ2fPLJJ+bauqVjtkCrRN6cwK9k9OjR7TLGJlDV9i8YLkca1dqpKa02oFPJy0y7I2q1mueee67JHHh14Zy0Wi1bt25ttgyVSsUHH3xg8TFboE09F4cPH+bQoUOUljbsiFq8eLFVjOpSzDX5ZZH3cXLkQrWOUidtMxfJyNguFov86aef5rXXXmPw4MENkq9JktSzRG7Umz4KNd7ufaBQR1kf7661TaZLmDRpEocPH7b4mK1gscjfffddDh06xMCBAzvCnq5HVW840FAOdlq8vX3gfCpVfbyp1ptmwsk0jy2NEfcE2vs+LZ4M06dPH/r169eum9o0CjtznLc6l93Vywd14SUAsrOzu8oym6cu22yVnBzSqtS9z7r3aykWX/Wvf/2LRYsW8de//hVv74bua0s5kjuL119/nW+//bbFTpQWkSRQOoGhxNz5pnR2wf5CKnp3L3KzswgKtI1JDraGWq3G2dmZCxcuYGdnh0IhT6hsL0ajkQsXLuDs7Nzm1NAWi9zNzY0ff/yxwThiXajixuGNOxe9Xm+d9pHKwSTy2ppcUihwzM+mBMjJu9T+8nsokiQREhLCyZMnOXPmTFeb02NQKBRERES0eVjOYpHfc8893HPPPcydO7dBx5st8L///Y/bb7+df/7zn+0rSFn7XIbLk19cy0vIAfJ0uvaV3cPRaDTExsZSVVXVZW1zvV7P3r17GTt2bJtrv47CUtskSUKj0bTLK7JY5Pn5+SxdurRDB/sXL17Mxo0bSUtLIzExscGKnqSkJObPn09eXh5arZa1a9cSHR2N0Wjkhx9+4NNPP21R5FVVVQ3ajLpa0er1evR601RWldIRCTBUFSNq97kaDWA0UmqAYp0ORwcHcxmG2lRKQghzGdagrixrlmktrmZbW9uP1kKv16NWq21O5GC5bTU1NU16ya39Xlj8PzF37lw2btzYofPUb731Vh5//PEmc0otWLCABx54gISEBNavX09CQgIHDhzgiy++4He/+91Vy37llVeanBSxfft2s2cyrtqAB5C4fwcXlSaX3b68HE1OJlW+QWze8iOO9eZlp1U4AcGU6ErYvHlz2x66BbZs2WL1Mq2FbFvbsIZt5eWtSw0tCQt9qlmzZrFlyxaGDh2Kj49Pg2NffPGFJUVdlZCQkAbRO3JzcwkLC6OgoACVSoUQAl9fX/bs2cOnn37Kjh07UCqVHDx4kGXLlnHfffc1KrOpmjwwMJCsrCzzAgblvjtRZP+AIW45ot88AJKXLiZRV07hqGnEhEcwZki8uYyD6RX89dscQj3UrLrN32rPr9fr2bJlC9OnT7e5Gkm2rW1Y0zadToenpyfFxcUtdnpbXJMPGzaMYcOGtcu4tpKRkYGvr6/ZFZQkiaCgINLT0/nb3/7G3/72N8CU86opgYOpzajRaBrlQmvgPmlMMd1UxhKo3Wfv1ReHYz9QOGoaufn5Df6DVKpqsz0d8aWyVbcTZNvaijVsa+31Fov82WeftdiYzqY1w2fNLlABsKvdri6+vMvTB8e0swBcKiygWq83Z1iRkbFlLO6y+/DDDxtkgATTXPaPPvrIakY1R51bbTCYYocJIUhPT2/T4vxml5oCqN1Mf/VFl3d5eqMuLsCuvAQhBDn5eW14AhmZzsdikT/77LP4+zdsd/r7+/P3v//dakY1h7e3N/Hx8eYflA0bNhAQEEBYWJjFZTUbrRXAzs30t7ro8q7afGhOF1IBuJiTY/E9ZWS6AotFXlBQ0ChelqenpzlTpDVYsGABAQEBZGZmMmPGjAYiXrVqFatWrSIiIoJly5axZs2aNt2jxZq8TuT6y+66xjcQAPtThwDIlKe3ynQTLG6T9+vXj507dzJx4kTzvt27dxMcHGw1o1atWtXsscjISKsEp2ixTa6ua5MXmXdpfPxAknA6Y2qq5OTnUVVdjUZOZSxj41gs8ieeeIJbbrmFJUuWEB4eTnJyMm+88QbLly/vCPs6jCt71xugqfVUqi63uxVqtUno2RdwsVNTUq3nQk4OoYGBnWSxjEzbsFjkt99+O87Ozrz77rt88sknBAcHs3r16m4XxLHFmty+duFNZQ4IYVq0gsllr8q+gLckKAEys7NkkcvYPG2aezhr1ixmzZplbVtshzqR11SCXmceUrP3C0KX+At9Sgo5Z+dCRnZWFxopI9M6WtXxdvLkyVYVdurUqXYZ05m02PGmcgR17Qyiylzzbnt/U63tfDEFSZIo0ukoKZMTLsjYNq0S+e23386cOXPYvHkzZVd8qcvLy/nuu++YM2cOt99+e4cY2RG0OIQGDV32Wup62GsyUvDxMLXb0y5c6FA7ZWTaS6tEnpiYyPXXX8+zzz6Lm5sbISEhDB48mJCQELRaLc888wzXXXcdhw4d6mh7Ow/72nn59UTuEGIayqtITSakNq52yoWMTjdNRsYSWtUmVygUJCQkkJCQQFZWFocPHzanSYqLi8PX17ej7ex86kRecbnd7dgvAoDy1CQG+fnzy+FEMrOz8Qru2OydMjLtweKON19f3x4h6haH0ACcaqfKlqWbd9kHhCCpVBgrK3As06F1caG4pITiwhyQEy/I2Ci9NgjXVdvkTrWTe8ouZ9FQqFQ4BPUHoCI1iRB/k8teVCD3ssvYLr1W5FfFLPLUBrsdQ2td9pSz9AswdcQVFWQjIYchlrFNZJE3h1OI6W9ZumlCTC0O/cIBKD9/Fl8vL+w1GgyGarzUxcjhxmVskV4r8hbHyQEcAwDJFLG16nKEVqf+AwAoSzqJQqEgLMhU4wdprLdAR0bGmvRakV+1Ta7UgGPtktrSyylunQfGmnadOoIQgvDgEAD87QoAYwdaLCPTNlrVu96vXz+LorMuWbKkZ+RFcx0A5ZlQfAo8RwHgFBkDCgX6/EtU52bh6+2LnZ0DVFfgKgqBHpxdRqZb0iqR79ixw6JC3dzc2mCKDaIdCNlbTSKvRWnvgGNoJOXJpyg5eRhPHz88vALIupBEH3JbKExGpmtolcituVa8W6GNMv3VNZy77xIVR3nyKUpPHsFz8nV4eptE7ka+HPtNxuaweDLM888/32Db2dmZwYMHM23aNKsZZTO4mjrZKD7dYLdzVCw5Gz+h9ORh07azOyU19rgoKzmXnsbA/paHo5KR6SgsFnliYmKD7bKyMl544QVuvvlm/vvf/1rNMJvANRKQTL3rlbnmRSsuMaaY67ojpk47SSGRUunFYKcMTp1LlkUuY1NYLPIvv/yy0b7c3FxiY2M5e/YsERERVjGso7nqtFYw5Sp3jQTdaSj4DfyuBcA1dgQoFFRmplKVcxEJN9IqvRjkmEHWpUsUFhfjfmUgChmZLqJNCQ+bQq1W88ADDxAaGgrA+++/3z7LOpgWI8PUx2OYSeR5v5pFrnJxxTlyEKWnjlB8cC/EXk+lsKOIPrhTwMlzyYyNH9pJTyIj0zIWj5NrtdomP0qlEicnJ/N2j8FjhOlvQcPxdO2wMQAUHdxrHl7MNfYF4EzK+S5P4ywjU4fFNfmKFSsa7Tt9+jQffPABK1euJCQkxBp22Q4etTPiCg6B0QAK0yvTDhvLhf+9S/GhfdTFay2kD44ODpRXVHA+I53wEHnMXKbrsVjkQ4YMaTAxpqysjLS0NJ544omeJ3AAl0hTKCi9DoqOQh9Tp5t26FgASo8fwr2sLme5RHRYOAeOHeXImdOyyGVsAotFvmTJkgbbdUNo4eHh1rLJtlAowWscXNwMOdvNIncICMYhuD8VaefQH9oFxIGAmPAIfjtxnJy8PHLy8vC5IhGFjExnY7HI58+f3xF22DZ9p5hEnr0NBj5q3t1nwjVc+N+76H/9CYLjAHB0cCA8OIQzKec5cuYU13iO7yKjZWRMtErkH374oUWFxsXFMXjw4DYZ1FZ+/fVXHn74YYQQTJkyhRdffNF6hftMMf3N3w/6ElC7ANBnvEnk1fu3QvCj5hXlsQMGcCblPOfS0igbMhQnR0fr2SIjYyGtEvn27dstKtTNza3TRT5kyBD27t0LwNSpU9HpdC0mZrcIl/6m9eVlqZC7E/xvAMB99GRTOKgL53EvSgc3U9QYrz4e+Hp5k3Upl6NnTjN6SLx17JCRaQOtEnlbkwp2JnUJ2WtqavDz88PR2rWn33WQ9A5kfG0WucrFFe3QsRTt30nE+W1kBvc3nx43MIqsS7kcTzpLfHSMnDNNpstoUwaVjmbx4sVs3LiRtLQ0EhMTiYuLMx9LSkpi/vz55OXlodVqWbt2LdHR0QCsW7eOpUuXMmPGDFSqph+tqqqKqqoq87ZOZ+oZ1+v16PX6Zm2SfG9AlfQO4uJmDJUloLQHoM/031G0fycDk34kY9J95jICfHxwd9VSqCvm8KmTxEdFW/we6spqya6uQratbVjTttaWIQlhe0GLdu3aRWhoKOPGjeOrr75qIPIpU6Ywb948EhISWL9+Pa+++mqDwA9Go5E5c+awdOlSBg0a1KjspUuX8txzzzXav27dupZrf2HkmupHcKCA/aqHyFaaZrRJRXk4P3s/Aon//PEn5g64nO64XK+noKoCBRJ9nZxRWLAmX0bmapSXl3PHHXdQXFzcYtPUJkVeR0hISAOR5+bmEhYWRkFBASqVCiEEvr6+7Nmzh8DAQDQaDQAJCQn87W9/a3JYr6maPDAwkKysLDw8PFq0R3Hs7yiT38Xoex01oy53Ru67cSyGkwfZc91z/PWNJ837jUYjn32/GV1pKaNi4xgcOcCi59fr9WzZsoXp06ebmyO2gmxb27CmbTqdDk9Pz6uK3Cbd9ebIyMjA19fX7IpLkkRQUBDp6ekkJiaycuVKjEYjEyZMaHbcXqPRoNFoGi1QUavVV3/p/edD8rsosn9AYcgDB1P8efspN1F68iBhx79Brf57g0uGRsewff8vHDlzmkGRA9q01rxVtnURsm1twxq2tfb6HhPjbc6cOezYsYNdu3a1avjsqjHemkI70BQGStRAykfm3fbX3IZRUtA3/SDlKUkNLokM7Y+rswsVlZUcPd19EkLK9BxaVZM3t/LsSjp65VmdW20wGMzuenp6OkFBQRaX1aqlpk0Rejfk/QLn18KAJaBQo/L241zwOMJTd5G1fi39//KS+XSlQsGo2Fh+3LuHxFMniYmIxL62WSEj0xm0qiavv9pMoVCwbt06srOz0Wg05OTk8Mknn6BUKjvaVry9vYmPj+ejj0y16IYNGwgICCAszPIgDW2qyQECbgSNlynAY8YXgClB0pGomwDI3vAh4oofjrDgEDzd3anW6/ntxHGLbZWRaQ+tEvmKFSvMn6KiItavX8/mzZt599132bRpE+vXr6egoMBqRi1YsICAgAAyMzOZMWNGAxGvWrWKVatWERERwbJly9o8hn/VuOvNoXKA8AdN/z61AoQRhQRn+k+l0l5LVfYF8nd81+ASSZIYFTcEgGNnTss5zWU6FYvb5D/++CPXXXddg30zZ85ky5YtVjNq1apVZGZmYjAYyMnJITk52XwsMjKSffv2cfbsWQ4ePNjkMFlraHNNDhB2H6hcQHcKLnyLQiFRo9JwKu5WADLef6PRJUG+fvh5+1BjNLIvsQeleJaxeSwWeUhISKO299q1a7tdRNc21+QAdm4QvsD072PPoxSm1MVHhs4DhYLCvT9RcvJIg0skSWLc0GEAJKWlcjE3BxmZzsBikb/99tv85S9/ISYmhuuvv56YmBgee+wxVq5c2RH2dRjtqskBBjwEGg8oScI162MAilz98L72FqDp2tyrTx+iw0xDe7sOHsBolDOuyHQ8Fot8woQJpKSk8OSTTzJx4kSefPJJzp8/z4QJEzrCPttF7QpRTwDgnvIPnCUdRgGB9z4MQM7GdVRlX2h02cjYODR2duQXFnKyXjNERqajaNM4uZubG3/4wx94/PHH+cMf/oC7u7u17epw2uWu1xF6N7hEoKy+xALXtzAK0MaPQjtiPKK6mtSVLze6xMHenhGDTfnUfjmSSEVlZdvvLyPTCto0423fvn3s2LGDvLw86s+KXb58udUM62haHa21JZR2MOwN2H4tsxy/ZEf1dUB/Qh9+jsTbp3Dx//5D8INPYO/fcBw/JjyCk8nJ5BcVsue3g0wfO679DyQj0wwW1+QrV65k6tSp/Prrr6xcuZKUlBRWrVpFdnZ2R9hn+3iNoSLwLgAecX4O9CW4j56E2+jJCL2e1Lcbz75TKBRMHjUKSZI4m5pC6oXMzrZaphdhschff/11vvvuO7788kscHBz48ssv+fzzz82LQ3oj5QOWklPjg78qEw6ZwkOFPmJa6Zb12RpKzzSeAOPj4UnsgIEA7Px1P9U2uCxSpmdgschzcnKYOHGi6WKFAiEE1157LRs3brS6cR2JVdrktUh2brxY+BI1QgFpn0LqJ7gNH4fXzJsRNTUkPf8wTS32GzE4FldnZ0rLy/nlcGITJcvItB+LRd63b18uXrwImMbMd+zYwYkTJ1Aoutdal3YPodVDoZA4po9jbekDph2/PQKFhwl76h8o7DQU7v2JvC2NfwTVKhWTRprynh87e4bM7Kx22yIjcyUWK/OPf/wj+/fvB+CRRx7hmmuuYciQISxatMjqxnUXVLVv8ePSuzF6T4WactgzFwdPewLvfwSApBceoaa88XTWwL6+5rHzrT//TGW9te4yMtbAYpE//PDD3HSTaTHGnXfeyfnz5zl69ChLly61tm3dBjulKeKLESXlw/5rSnlckQV75xL8wGI0foFUZqRw7p9/a/L6sUOH4ubiSllFOTt+/aVJ115Gpq20ycc2GAzs3buXTz/9lMDAQIKCgijrxYsu1EqJusBOVZILjPs/02y4wiOoDj/IgJdMswEz175F0cG9ja9XqZk+dhwKSeJcejqnz5/vROtlejoWi/z06dMMHDiQ2267jXvvvReAn376ifvuu8/qxnUkVu14kyQ0KpPMq2sEOPeDsf8HSkfI2YaH+mN85ySAEJz6yz3UVJQ3KsPbw8M8SWb3wV8pLC5udI6MTFuwWOQLFy7koYceIjMz0xx+ZtKkSezZs8fqxnUk1ux4A7CrE7mh1tX2HGGq0RUauLiZsOml2Pn4UpGSRNLzDzdZxpCoaPx9fNAbDHy/eyd6g8Eqtsn0biwW+eHDh1m4cCGAOfGhq6srJSUl1rWsm2FX2/tWbai36MRnIoz9GBR2qIt+IGp+CEgSFz9ZTc7G/2tUhkKhYPrY8Tg6OFBQXMzugwfk9rlMu7FY5D4+PqSmpjbYd/bsWQICAqxlU7ekrvOtynCFKH2nw/jPQOVEH8/jhFxvmuJ6+q8PNIoHB+Dk4MCMceORJInk9DTKDPIkGZn20aYhtFtuuYVvvvmGmpoafvzxR+666y7+9Kc/dYR93Ya6NnkjkQP4TIaJ34CdOyFj83ELc6SmrJSj98/GoGvc9vbz9jGnViqqqiQ7L69DbZfp2Vgs8sWLF3P//ffz17/+lZqaGpYsWcJdd93Fgw8+2BH2dRsc1KZXWaFvZo24x1CY8gMKbSjRcyU0WgXlyac4/ufbG8WEA4gbMJB+td7Rj3t3oyst7TDbZXo2bRpCW7hwIcePH6e0tJSTJ092y1rcmr3rAC4a06ssqWohEIRrJEzdhiZsEoPn26NQQ8HO70l+8ZFGp0qSxKThI1ErFFRWVbFpx3aq9dVWsVWmd9GmpaZlZWWcPXu2UWdbdwocYZWlpvVwrhV5aeVVQjxr+sCEL3Bx/StRBe9w/KNKMta8hZ27E8GLG64/V6vVeNg7ojMaKCgu4sc9e7hu4qRuN4VYpmuxWOSfffYZ9913H9XV1Tg4OJj3S5Jk1Yit3Q0Xe1NI6hZr8joUaoj/F95eY+lffDfnvtFxbvky1KpC/P74DtTLmaZSKLhmzHi+2b6NtIsX2HvoN8YNHWYe2ZCRuRoWVwmPP/44K1eupLKyksLCQvOnNwsc6tXkVRYkawi8ieBXEgmaGQjA6X+uInfFRKjMbXCadx8Ppo0eC8DRM6c5JMdul7EAi0VeXFzMH/7wh46wpVvjWluTF1dYmJHFuR/93z6L34zhIODE23vIeTka0j6DemPkYcHBjI03ZVL95chhTiY3Hn6TkWkKi0X+hz/8gS+++KIjbOnWeDiZWj55ZRaKHJBUGiLf+Zm+s2YhjHDif/lkrUxAuecmXIyXo8bEDYwiPjoGgB2/7udcerp1jJfp0bSqTX7TTTeZ24AGg4E777yT1157jb59+zY4rzeL38vZ9CovlbZtKqqkVDLwjS9ROC/g4if/5dRnldTotzFp9D7EsQwY9BSoXRkVG0dlZSUnzyXz497dzLKbQkBfX2s+ikwPo1Uir8sPXsfQoUM7wpZ2sWfPHh577DEUCgW33HILjz76aKfe37NW5HmlBoQQbeoYkxQKIl/+NwqNI5lr3+Lsl1VUFgr6z3wHMtZD9BNIoQlMHDGSyuoqzmdksGnHdm6YPAV/n75Xv4FMr6RVIn/22Wc72o52Exoayq5du7Czs2Py5Mn88Y9/xNHRsdPu7+lkapNX1wiKKmpwd2xb6ndJkgh/9nXU7h6krFhK+o5qKnUODLwlF+Whx+DsOygGPcP0MbP4btdO0rMu8u32bbLQZZqlxwy4+vn5YWdnB4BSqez0sWQ7lQIfF5OwMwrbN99ckiT6PfQMEf/4L0KpIvdQMYfX+VBd5Q6l52FfAqod07k2wkiQrx+Gmhq+3b6NCzm9NGKuTIu0rbrpYBYvXszGjRtJS0sjMTGxQXMhKSmJ+fPnk5eXh1arZe3atURHR5uPb9myhf79+2Nvb99k2VVVVVTVC7Gk0+kA0Ov16NsZMTXIXU1OiYGUvAoGerf/1fa54TYq0i/g8sFrFJ84x69v9SXq0bm4q75BKjiE6udbmaEdyvfu95NRaODb7duYOX4Cft4+7b731ah7V+19Zx1Bb7GttWVIwgbXMu7atYvQ0FDGjRvHV1991UDkU6ZMYd68eSQkJLB+/XpeffVV85rwzMxM5s2bx8aNG3F2dm6y7KVLl/Lcc8812r9u3bp2u/c7CnxI1HkQ75rPxD7WS2go5VzA4b/LUOZkIhRKam6cS7+x1fQzbkNFNQahYmPNYrJEKAAe9g44qNRWu7+MbVJeXs4dd9xBcXExrq6uzZ5nkyKvIyQkpIHIc3NzCQsLo6CgAJVKhRACX19f9uzZQ2BgINdffz0rV64kMjKy2TKbqskDAwPJysrCw8OjXfZ+f6qE17bnM9hPw2uz29/jrdfr2bJlC9OnT0dRXcXZpxaQt3k9AJ7XzyH8b0uxy1mH4vx/qTFU8UPNvaSKQUgIJg4bQkTogHbb0Brb6oKH2Aq9xTadToenp+dVRd4mn3L79u2sW7eO7OxsvvnmGw4ePEhJSQmTJ09us8GtISMjA19fX1Qqk9mSJBEUFER6ejq7d+/m5MmTLFhgSin88ccf4+/v36gMjUaDRqNh5cqVrFy5kpraFWBqtbrdL31wgDOQz5ncaiSFCpXSOlNP1Wo1akdHBq38lMz33yD5lcfJ2/Q5JYm/EPXaB7jf8DCqs+8wM+k/bK8s54wYyY6DhzFc2ELsqNvAsfF7sBbWeG8dRU+3rbXXW9w79Z///Ie77roLHx8fdu3aZb7ZM888Y2lRVuXuu+/m4sWL7Nixgx07djQp8PpYO/wTQKC7GmeNgiqDIDnP+qGVJUki8N4lxH+6E/ugUKouZpB4x1SS/vUyxoi/oJx1nKnDYonV/AbAniwv9n39V8S++6HgN6vbI9M9sFjk//jHP/jxxx958cUXzT3YUVFRnDp1yurGXUmdW22ojX0mhCA9PZ2goKCrXNkYay81BVBIEoP8TIt2DqY1DtZoLbRDRzPiu8P4zr0XhCBj9XIO3DgC3clTSBH3M/amfzEyxGTHIeN0fkxxwrDlGtgyAc6tAX3vDtXV27BY5Pn5+URFRQGXY7xJktQpq6K8vb2Jj4/no48+AmDDhg0EBAQQFhZmcVkdUZMDjOnnBMDPKR0bolrl5MzAZasZtPor1B5elJ0+xsGbRnN26UPUlJczbOytTBk1BoUkkSyG8nXNQ1QUnIPflsA3A+DgEig80qE2ytgGFos8NjaWDRs2NNi3ceNG4uPjrWbUggULCAgIIDMzkxkzZjQQ8apVq1i1ahUREREsW7aMNWvWtOkeHVGTA4yuFfmp7Eryyzo+2qrX9N8x8odj+My+E4xGMte+xf7p0Vz6/ksG9u/PrCnT0NjZkS36sV75EoWOw8FQCufXmGr2rZPh/Idy7d6TERby22+/CQ8PD3HzzTcLjUYjbr/9duHj4yMOHz5saVE2QXFxsQBEXl6e1cpc+Gm6mPj6WbHuQH67yqmurhZfffWVqK6ubtX5+bt+FD9PCBM/BUvip2BJHLn3d6Ls/FlRUFQkPvzqC/H2Rx+Kf3/6iUg5/p0QPycI8bmHEJ+6mj7r+wrxy/1CZP0kRI3B6rZ1Jr3FtrrvbnFxcYvnWVyTx8fHc+LECUaPHs19991HbGwshw4dIjY21vq/QB1IR9XkANdHm4Yzvjmuw9iJI5R9xk9nxA9HCf7TU0gqFXlbvzHV6q8/y42jRuHr5U21Xs+mw5c44Pww4vqTMGgpuISZ8relfQq7boJNMXD0WdCd6TTbZToOi0RuMBjw9fVFq9Xy2GOP8fbbb/PEE0/g5+fXUfZ1GB3VJgeYHOGCk52Ci8V69p7r3PRRSnsH+j/2IiO+P4LH5OsQBgOZa97k8DUxxKefJKa/qenz69EjfPfrCarDFsHMgzB1K/S/B9RaqLgIp1+H70eY3Pmkf0NV7w4K0p2xSOQqlQpnZ2eqq+WAgi3hoFZwU6wpbtz7v+R3am1eh1PYQGLXfEvc/37AacAgDMWFnH/xUeyemMdQjYRSoSAlM5PPv/+O/OIi8BgOQ1fA75Jg9Afgdy1ISig4BIl/gW8iYM8dkL4BDL037113xGJ3/W9/+xt33303x44do6ioCJ1OZ/50JzrSXQe4Ld4dJzsFKfnV/Hiq6zq1+oyfzohNhxiw7N/YefpQkXaO8sfmEf71+zhIUKTT8fn333EyOcmUrUWpgcDZphRPs85A3CvgNhiMeri4CX65B74OQ3ngfvrW/AY1cqplW8fiaa31V3fVDZuJ2vXTNU3ED7d16qK15uXltXta65WsO1jAv/fm42qv4MN5Ibg5KC26Xq/Xs3nzZq677jqrzNwylJaQ+cFbpP/7NQzFhRicXMiZ/yg6/34AhAeHMGnkSOzUdo0vLjoBGRtMNXlZqnm3ULsi+c+CoJvBe6IpSGUXY+33Zk2saVvdd/dq01otrslTUlLMn/Pnz3P+/Hnzv2Ua8vsh7oR62KGrNPLPrTldntdM5exCyKKnGLP7PP2WPItGIeH37nN4f/8pGI0kpaXy6aZN5ObnN77YLRoGPQPXHYZp26gJ+yMVuCPpdZD6Mey6Bb6JhN8ehtw9IFoRtVamU7B47npwcHBH2NEjUSklnrzGh0WfZbL3fBmfHipi7lD3rjYLlauWfkueJeDuxWT8ZznqD1bikHqGC7ctRAes/24Tw6OiGBoX33hdviRBn6EYXQbzY/pIrh/pgeri15D5FVTlwbn3TR8HXwi4CYJugT5DG4SZlulcLBb5Pffc0+yx999/v13GdCZXLlDpKCK87fnTBE9WbL/Eqj15+LiomBzh0qH3bC1qrTuhj75A0ILHufjJv3H+3zukj7uOkuhh/HrqFGcPHWDq2HH0HTio6QIkBcJzNPhOgCGvQu4uk0uf+Q1UZEHSO6aPUzAE3mwSvDZGFnwnY7G7rtVqG3zKy8tZv369zbV9rkZHDqFdye8GablxkBYBvPRDNvs6eMqrpaicXQi6/1EmbDnO1NhYQnZ9i6KynCIHZ774ZT/f/30JeTu+RxhbcMEVKug7BYavNPXQj/s/CJoDKicoS4PTK+DHcaZhuROvQkly5z1gL8fimnzFihWN9m3fvp3Vq1dbxaCeiCRJLJ7kRXFlDTuSSnn624s8Od2H6QOa7yzpChQaDf633YPfnAQytn3HrlOnKO7jw7kBQ8ne/hOhrz1Dv9l30Pfmu8C5BduVGtMQnN+1YCiHrB9MHXZZP0LJWTjxsunjNthUuwfeDE6WLzKSaR1WCYQ2adIkNm3aZI2ieixKhcTTM/oyLdKFGiO89EMOa3/Jp8ZoezE7JIWCoGnXc+efHmVESDAKo5GyiMGcuOk+9m/ZzO5RAZxechfKM0dbrt0BVI4QeBOM/QhuTIYRq6DvdJBUUHTUNLNu0yD46RpIWgUV1ouoI2PC4pr8yvHw8vJyPvjgg0Yx2GUao1JKPDXDBzcHJesPF7F2fwEnsir52wwf3NoY3bUjkSSJ4WMnEBYTy7af95BdUEDOrLvQxY6m6sv3cdz0GQc3rsFv7r30vWU+9r4BLReodoWQuaZPVT5kbjTV8Jf2QP5+0+fwk+A13lTD+88yJYiUaRcW1+Rubm64u7ubP35+frz33nu89957HWFfh9HRk2GaQyFJ/GmiF3+9xgeNSuJAejl3f5TOtrMlXT7E1hzuWi03z7yOicNHoFapqQgKI+XPL5I74/eUZ2Vy/l9/5+cxwRy6fSoXP3sfg6746oVqPKD/3TD5W5h1CuKWmWbdCSPk7oSDi2FjOOz+fe0su4qOf9AeisWTYdLS0hpsOzs7W30SSWfSkZNhrsb5vCqe/z6b1HzTNOEx/Zz48yQvfF1NnZi2OKmjtLyMnQd+JTXTlL7JCUHgz98jbfo/8zkKOw0e02bRd/adeEy6FoVdE5Nrmr1BKmR8YfoUHbu8X601zcQLvh08R7XYQ2+L762OrpgM0+5AjomJiahUKgYNamaYxcbpSpEDVBuMfHywkI8PFGAwglopcdNgLXcO74OjymiTX1YhBGdTUtj2y8/meflBHh70TzmGbsOHlCdfjhKkcuuD14zZeF97K+5jplgm+OLTkP65aXVcecbl/U4hEHI7BM8F55BGl8kib4jF7vqNN97Inj17AJPLO2bMGEaPHt3t3HVbwU6l4O5RHqy+I4j4QAf0NYLPEou4Y20qHx0sorLG9vJfSJJEaGAgfR2dGRQRiSRJpOfns8s9EMXydQzZeIDA+x7GzqsvhqICsj59nyMJ17FnWF9OPno3eds2YaxqxZx37QAY9He4/ihM+hZC7gSVs2la7YlXYHMsbLsWUj+R3fkWsLgm9/b25sKFC6jVaqKioli9ejVarZabbrqJpKTul063q2vy+gghOJBezqo9eZzLM7nwKsnIDTFabhvmYXbjbYH6NZKurIzdB3/lQo6pZ9zV2Zmx8UMJ8fWj6NddXPpuA5e++4LqvMs950oXVzyn/Q6va26kz/hrUDm3coKQoQwubILUdZCzA6j9+tq5mVz5/nejdwiVa/J6WCzyukIvXLjA0KFDyc42peZxdXXtdivRwLZEXodRCLafLWXdwQKz2BUSjAt14voYLcOCHFEqunbW2JVfViEEyelp7P3tN8oqTEEsfb28GRMfT19PL0RNDUUH93Lpu/XkfvcF1TkXzWVJdna4j5yE57Qb8Jg6C4eAVk6dLr9oEnvKB1B2OY2z0WMUicVxDL7+GdQaJ6s+d3vpFiIfPXo0v/vd70hLS6OiooIPPviA3NxcYmNjycrKapfRnUn9aa1nz561KZHXUV1dzXtf7CJVHcmhzErzfh8XFddGuXJtlCs+XVS7N/dlrdbrOXTiOIdPnzJPGQ4LDmZU7BC0LqbaWhiN6BJ/Ife7DeRt/YaK1Iaz35wGDMJz6g14Tp2Fa+xwJOVVVu8ZayBnuylu3cXvQJjuK+x9kMIeMAXDsJGhuG4h8sTERP70pz9hZ2fH+++/T79+/fjf//7H1q1b+eCDD9pldFdgizV5HfW/EBnFRjad0PHjKR0lVZcnoAzys2dyuAsTwpzN6ZM727amvqyl5WXsP3KE0+fPAaYlyoMiIhkWMwh7jcZ8nhCC8nNnyPvpG/J++pbig3uh3gQbtbsH7mOn0WfCNXhMuAZN36skiii/SM25D6g+tQoHCk37lI7Q704I/yO49G//w7eDbiHynkZ3EXndF6LKYGTPuTK+PV5MYublziYJGOzvwKRwZ8aGOuHt0rE1fGu/rHmFBfx86BAZ2SYvT2NnR9zAKAZHDsCuiev0hfnk7/iOvK3fkr/re2pKGjYBnSKi6TPhGvpMmIHbiPEo7R2atO27TRu5blAlqnPv1huKkyDgRoh+ArRRbX/4dtAtRL5582YiIiIICwsjNTWVhx56CJVKxRtvvEFAwFVmPNkg3U3k9blUamBnUgnbk0o5kVXZ4Fh/TztGhTgxqp8TUX3trd6Gt/TLmn7xIj8n/kZ+UREA9hoN8VHRxEREolY17YEY9Xp0h/dTsOsHCnZtQXf0ANT7uio09mhHjMd99GTcR03CZfAwFCpVQ9tUKtPquLNvm+bOAyCZxtyjngDtwHa+CcvoFiIfMGAAW7duJSAggDlz5mBnZ4eTkxNZWVl888037TK6K+jOIq9Pjk7PzuRSdiaXcjKrkvr/qa72CkYEOzEsyJH4QAer1PJt+bIajUaS09L49dgRiktMIbEc7R2Ij44mOjwC1VXa3vrCfAr2/lQr+h+pyr7Q4LjSyRm34eNxGT6OEzVqZty/CDuHejV90Qk4+Spkfl27QzINyw162rT+vRPoFiJ3c3OjqKgIg8GAl5cX6enpaDQa/Pz8yMvLa5fRXUFPEXl9iipq+DWtjF9Syvg1rZzSqoaLSPy1auIDHRgS6MiQAAfc2zBvvj1fVqPRyJmU8xw4dpSSMtOyW2dHR4ZGD2Jg//4or9bRRm1bPvkUBXu2UrhvB0X7d2IoLmxwjtLZFbcR43EfPQm3ERNwjopDoVbXin2Zae48mNrskYthwEOmBTUdSFeI3OL/XQcHB3Jycjh27BgDBgzAxcUFvV5vkwnfeytuDkquGeDKNQNcMRgFJ7Mq2Z9axqHMCs7kVHKhWM+FYj3fHDe1d/t52DHYz4FoP3tifB3wdVV1aNorhULBwP5hRIT049T5cxw8dozS8nJ2HtjPweNHiRsYRXR4OOoWcqxLkoRTeBRO4VEE3r0YYTRSeuoohfu2U/DzNvL3baemVEf+tk3kbzOtkFQ4OOIaNxK3YWPQDr0T7Yh7UJ17GfJ/NYk+7RNTxNq+Uzvs2bsCi0U+b948hg8fTlVVFc899xwABw8eJDQ01OrGWUJBQQHTp0/nzJkzlJaWdqkttoRKITHY34HB/ia3tbSqhqMXKkjMrOBQRjnn8qpJyTd9vj5mWljSx1FJjK890X4OxPjaE+6lwU5l/Zl3SqWSmPAIBoT250TSWRJPnqCsooK9h37jtxPHGRw5gEERkQ1645tDUihwiY7DJToO3/l/YvO33zA+xJ/SA7sp/GUnxQf3YtAVUbRvO0X7ttdeJOEcOQjtwLFoXY+g9UvBvvQmpJC5pkg3dl0fqssaWCzyV199lWnTpqFWq5k0aRJgSl382muvWds2i3BxcWHLli38/ve/71I7bB1njZIxoc6MCXUGTK79kcxyTmRVcjyrkrO5lRSU17DrXBm7ahNDqJUSkd4aon3tieprT4SPPR721huUUSmVxA4YSEx4BGdSzvPbiRPoSkv49egREk+eICY8gtiBA3FysMCVVihxiYmnz5CRBD3wGMJopCz5FMUH91J88GeKDu6hMv08paePUnr6KHWteztXCdfAtbj2/xrXmU/jOn4+KhfbCu5hKW0eQrtw4QIZGRmMGjXK2ja1i2nTprF169ZWn98T2+Ttocpg5GxuFccvVnA8q5ITWZUUVTSOg+eiUeAm6Rg9wI/Ivo5Eemvwc1OjsIKbbzQaOZeexm8njpt745UKBZGh/YkdMJA+Wm2L17f2vVXlZlH8288UH9xL0cG9lJ5IRBiuSFIpSTiFDcQldgSuscPRDhmJU+QgU9u+DXSLNnlWVhZz585l79692NvbU1paymeffcaWLVusFgJq8eLFbNy4kbS0NBITE4mLizMfS0pKYv78+eTl5aHValm7di3R0dGtLruqqoqqeosj6qbi2mK/Qp09nWmXAhjgpWKAlwu3xroghOCizsDJrCpOZFdy9pLJtS+pMlKCMxmHdYDpHTqqJfp72RHuqSHMy45wLzsC3dRtGr4L8Q8g2M+f9KyLHD51kpz8fE4mJ3EyOYnAvr4MiojE38enyb6D1r43hbsn7tN+h/u03wFQU1FO6YlESg7tpnT3B5ScSaGyUFCWdJKypJNkr19ruk5jj1NUHC6Dh+EcPQTn6CE4hkYiNTMU2BbbWkNry7C4Jp89ezYRERE8//zz+Pr6UlhYSH5+PsOHD7da7PVdu3YRGhrKuHHj+OqrrxqIfMqUKcybN4+EhATWr1/Pq6++2iAY49Vq8qVLl5r7Euqzbt06HB07tme1p2AQEvnVGnKr7U2fKnsu6e2paSJ/pkoy4mVXibf5U4GHXRVKC3QvhKDaWENJdTWVNZdrWpVCgYvaDkeVukM6CoNqdjKw8H1KM/TkZLhyMdMbRfp5pIrGgTiF2g6jXzA1AaEYA0KpCeyP0TcIWug8bC/l5eXccccd1h9C8/LyIisrC5VKRZ8+fSgoMCXCq3MbrElISEgDkefm5hIWFkZBQQEqlQohBL6+vuzZs8ecw/xqIm+qJg8MDCQrK8sm3fUtW7Ywffp0m1xNVd+2GqMgvVBP0qVqki5VkZxXTfKlaioNjb9eaoWpRz/cy44wLw3hnnb081C3qnOvuKSE40lnOZOagqHWtbbXaIjqH0ZUWDiO9vZWfW9S3l6UvyQg6QsR2kHox3xORVYBJUcOUHr8EKUnDlN26jA1ZY07eyWVCsfwaJyj4nCOjsMpKg670Ei279tvFdt0Oh2enp7Wd9e1Wi0FBQV4e3ub92VmZuLj49M2Sy0gIyMDX19fVLVukSRJBAUFkZ6eTlhYGNOmTSMxMZFp06bx+uuvExMT06gMjUaDRqNpFHddrVbbnJDq6A62qYGIvnZE9L286qvGKLhQrCcpt4ozuZUk5VZxNreKsmojZy9Vc/ZSNWASh1IBIX3siPC2J9xbQ4SXhv5eGhzUDYXv2acPk0aOYvSQeE4kJ3HszGlKy8s5dPIEh0+fIjQwkAH9+iOEsM57850EkzfBzhuRio9ht+/32E3ejDYiGuYkAKYFNxVp5yg5fsj0OZFIyfFDGIoKKDt1hLJTR8jZcHldh5OHD2c3Dsc1KhbnAYNxihyEQ0gYila4+/Vp7bNZLPI77riDO++8k+XLl5sihJw9yyOPPEJCQoKlRVkdSzrcFi1axKJFi8ydFzLWR6mQCHK3I8jdjqmRtSvQhCBLZ+BsnegvmYRfXFHDubxqzuVV891J0/USENTHjnAvDZHeGsK9NYR5aXDWKNHY2REfFU3sgIGcT0/nyJlT5OTlkZyWRnJaGiqFguNJZ4kKC0djSTSapnCLNgl9+3WmCLP7Ekxx5WvzvkkKBY79wnHsF47PrNvMz1l1MaNW9IcoOXGYkuOHqM65iCI/h4KfvqXgp2/Nt1Bo7E3j/gMG4zxgUO1nMHae3k1ZZBEWi/yZZ57hqaeeYsyYMZSVlTF06FAefPBBHn/88XYbczXq3GqDwWB219PT0wkKsjxmd2dlUJFpiCRJ+GnV+GnVTAq/LPxLpQbO5lbVE34l+WU1pBVUk1ZQzdYzlzPD+mvVRNSKPsJLQ0TfQG4NCeFSQQHHk85wNiUFQ00NPyce4tdjR4kICSEmPBKvPu1YbuoaCeM+gx3XQ/ZWOPYCxD7f4nPa+wdh7x+E14zZ5v3ll3LY9uF/iHFzoiLpBGWnj1N65hjGinKzJ1Aftac32vgx9Bk/He+ZN2PnZbnHbFGb3GAwsGrVKu69917s7e3Nw04dNTvqyjY5mGK8JyQkmDveli1bxsGDB9t8D3kIrW10hm35ZQaz6JNyKzmbW0VOiaHJc/20aiK9NUT62BOsNXL20HZUDmqK6gUy8fH0JCosnLCg4CZXwLWKjK9g33zTv8evB9/pFl3e1HsTRiMV6ecpO32sdtz+GKVnjpnW2deTp71/MGP2ppi3O3zuekeyYMECNm3aRHZ2Nh4eHri4uJCcbAoscObMGRISEsjPz8fV1ZU1a9a0KYhkdwga0dtF3hTFFTUmwV+qqm3rV3GxuOmhpAA3JTF9KvEki+qSLHPIa5VKRVhQMAP798fXy9vySurQY5C8Ghz8YOZ+Uzz5VmLJe6spL6P01FHytm8i7e2XAZiSenkdQoeJ/LbbbuP+++9n2rRpllxms8g1eduwJdtKKms4W9u5dyanijM5lWRfUeNrpGr62V8izDEPB+nyOnytswsD+vdnQGh/nFs7hGqogB9GmQJKhj9omgLbStry3qov5bBnuGmVXFtEbnGb3MXFhdmzZ3PNNdcQFBTUILXt8uXLLS1ORqbduNgrGRrkyNAgk0j1ej0bvvmBoNgJnMs31IpfxelSf05X+OGhKiXEPpdAu3yKS0vYf+Qw+48cxs3dh7jIcAaEBLW8Ek7lAEOXw66b4dx/IWKhKXOrjWKxyGtqapgzZw6A1cfFOxO5461n46CsYXiQA2P6X64t88sMnM4xTdU9keXNjzlleCnz6Gd/CS91CUWFOez4JYetv6hQO/sS2S+E0QMCcbRrQvB9p4L3RFO2l5P/hOFvd+LTWYbFIl++fDnu7t1/dY48hNb78HBSMTbUmbG1i3MMNYLkvGBOZlVyMjOfkvx0PKVcHJXViNIMTh/L4NBhO8rVPvj7BhPf35vBfo6o6qbrxTwN23aa4r7HPA0OtpkPsNUiP3ToEDfeeCMXL14kJCSEjRs3WjRnXEbG1lApJQb42DPAx56b49yA/mQXV/Hz6UzSMlNRVOTgqKzG0ZhB1YUMvk9zYK3BCy/vQEaGeTIiZChunqMg7xdI+R9E/aWrH6lJWr1I+LHHHmPOnDkcPXqUWbNm8Ze/2OYDtZauSngoY9v01Wq4eWR/Hr5lKgvn/p4RQ8fiqPVFIKFVVRBpn04f3V727d/GYx/+zNpLpsUthnMfNhjusiVaXZMfPXqU77//Hjs7O15++WXzXPHuiuyuy1wNtUrF8AGhDB8QSlV1NUlpaRw5e46iokt4qkvwVJdQqndlg/Qw4TWH+PbT7wmPHMXUCBf6ONlOKupWW1JdXY1d7fRAR0dHKisrr3KFjEzPQWNnR0x4ODHh4ZSWl5OUlsrplFQKCvPJpj/Zxv5oyePkkR1sOeCBX98Arh3sYxPZbiwS+ZtvvmnerqysbLANpnXg3QW5d12mrTg7OjJkYBRDBkZRUlbGucNfci4tmWzRH091KZ7qUihJ44cdzqxXeDEqKpQb4vqi6YAQWq2h1SIfNWoUX375pXl75MiRDbYlSepWIpfddRlr4OLkRFz8VOIuLKIUd85Ff8WJtCwKCvPwUJfiQSnZZ1J45YQ7AwcO54ZBbp1uY6tFvmPHjg40Q0amG+PgC46BOJdnEOtZTGz0tZRVlHM6JZ1DZ89TXZaPl6qQfUePsv5Yf4bbu3JtJ3bS2V7yaxmZ7ohH7ShNvilKkZODI0OjBnD/7OuIjzatrXCzq+ZSaQ2b8wL467c5ZOs6J6yXLHIZGWvgHmv6qzvV6JCrsymQxrAgR+YNd0OJkd8yKrlvXTq7kjs+fHivFbk8Ti5jVVwiTH91Z5s9RSHBXcPduMv/PAN87CitMvLMpiw+3J9PR+Yd7bUiX7RoESdPnmwQBFJGps241oq8JBmEscVT3dXVrJjty+/j3QB4/5cC3tp5qcOE3mtFLiNjVZxCQFJCTTlUZF31dJVSYuF4Lx6e7AXAF0eKWftLQYeYJotcRsYaKFRgXxuaqSK71ZfdONiNR6aY4rh98GsBP5+3fhtdFrmMjLWoW4VW2XqRA/xukJabYk1zNV75MYfC8qZDXLUVWeQyMtbCvlbkFtTkdSwc70W4l4aSKiOrf863qlm9VuRy77qM1WljTQ6mpJIPTTK1zzef0JFWUG01s3qtyOXedRmrozGJlKq8Nl0e4+fA2FDTmPrniYXWsqr3ilxGxurURW3Vt73z7PfxpqhLW8+UUGVoeSiutcgil5GxFmpTsggMJS2f1wKD/ezxdlZRqRccSCu3ilmyyGVkrEWdyPVtF7kkSYwPM8Wg23u+cfbUtiCLXEbGWqjqRK5r+byrMCLYFFr6eFbFVc5sHbLIZWSshRXcdYABPvYAZBTq0VW2P6iJLHIZGWuhar+7DqB1UOLragr1kJLf/qG0HiXyP//5z4wfP56XX365q02R6Y0oTTUwxvYLM8DdFE8xs1AWuZmDBw+iUqnYvXs3hw4dIicnp6tNkult1OYrt4bI/bWmsi40k8zREnqMyPfv38+UKVMAmDhxIr/99lsXWyTT61CYal+M7RdmX1eTyHObSdVsCbYTHLoeixcvZuPGjaSlpZGYmNggP3lSUhLz588nLy8PrVbL2rVriY6OpqioiJiYGMCUlLG59MpVVVVUVVWZt3W1+av1ej16feeE42ktdfbYml0g29YkNaAGEDXoqytNS0+BGoOp88xoFK22zaX296KwXI/ecDmkc/3rWvt8NinyW2+9lccff5xx48Y1OrZgwQIeeOABEhISWL9+PQkJCRw4cAA3NzezYEtKSujfv3+TZb/yyis899xzjfZv374dx9amru1ktmzZ0tUmNIts22VUooLra//9/eZvMEompZbqTe57Tk6O2aar2ZZS7gwEkZFTyE9bE3Gu3b9582bzOeXlrZssY3F+8s4kJCSEr776ylyT5+bmEhYWRkFBASqVCiEEvr6+7Nmzh8LCQj755BOWL1/OnDlzePvtt/Hx8WlUZl1Nvnr1alavXk1NTQ3JyclkZWXZZH7yLVu2MH369C7PAX4lsm1NUFOBemOgyYYbUsxDaqfOnWP3bwcI9vNnyshRrbLtdE4Vf96QhbezkjXX2bF/dBAA45MaeqGenp7Wz0/elWRkZODr64tKZTJbkiSCgoJIT09nypQprF27lvHjxzNjxowmBQ6g0WjQaDQ8+uijPProo+a462q12ua+rHXItrWNTrdNebmLS60Eau+tVJncdoVCMttzNdvcnU11b0mVEbXq8nn1r2nts3UrkV+NlStXWnSunEFFxqoolCApTDHe2tn5Vpdtpbqm/Y52t+pdDwwMJCsrC4PB1OMohCA9PZ2goCCLy5KXmsp0CJJ1htHsanOg1xjBYGyf0LuVyL29vYmPj+ejjz4CYMOGDQQEBLQpw6ocNEKmQzAPo7VT5KrLPertrc1tUuQLFiwgICCAzMxMZsyY0UDEq1atYtWqVURERLBs2TLWrFnTpnvINblMhyDVSap9wqyryQH07RS5TbbJV61a1eyxyMhI9u3b14nWyMh0PkqFhEoBBiPo2xk8wiZr8s5AdtdlOhQrjEzbWanzrdeKXHbXZWydOpe9ve56rxW5jEyHIElXP6e1RVmpnF4rctldl+lYbGciaa8Vueyuy3QX2tu877Uil5HpGKznrlurqF4rctldl+lYZHe9y5HddRlbp64ib+/PRa8VuYxMx2BFd91KyCKXkekIbChMgyxyGRkbpW7IXXbX24jc8SbTIVhxMoy16LUilzveZDoW2V2XkZG5CmafoJ3te1nkMjJWxVoDX9ZDFrmMTA9HFrmMjK1S24kn967LyNgUdeNesrve5chDaDK2jryevJ3IQ2gy3QV5qamMjC0hyb3rMjIyrUR212Vkegly77qMjE0hu+syMjKtxFprXXqMyAsKChg6dCjOzs5XP1lGpjsh966bcHFxYcuWLYwaNaqrTZHp1ciTYToMtVpNnz59utoMGRmbo0sTHi5evJiNGzeSlpZGYmIicXFx5mNJSUnMnz+fvLw8tFota9euJTo6ut33rKqqoqqqyrxdXFwMmNx9W0Ov11NeXk5+fj5qtbqrzWmAbFvTqMoNSJUCfWEhGPMBKC4qpqK8gtKSUvLz81ttm75Ch6HSSFGxgrLaHOX5+fnm4yUlJQCIq3kNogvZuXOnyMjIEMHBwSIxMbHBscmTJ4s1a9YIIYT4/PPPxbBhw4QQQpw4cUJMnDixweeVV14xXzd16tQW7/nss88KTK0c+SN/esQnIyOjxe+8JETXNx5CQkL46quvzDV5bm4uYWFhFBQUoFKpEELg6+vLnj17GuQqb4pp06axdevWZo9fWZMXFRURHBxMeno6Wq3WKs8zfPhwq0yX1el0BAYGkpGRgaurqxUss13brGUX9B7bhBAMHTqUs2fPolA03/K2yfzkGRkZ+Pr6olKZzJMkiaCgINLT01sU+bRp00hMTGTatGm8/vrrxMTENDpHo9Gg0Wga7ddqtVYTklKptFpZAK6urj3eNmvbBb3DNjs7uxYFDjYq8rbSUg3emSxatKirTWgWW7XNVu2CHmCbZa3ojuHKNnlOTo5wcXERer1eCCGE0WgUPj4+Iikpyer3Li4uFoAoLi62etntRbatbci2NcQmh9C8vb2Jj4/no48+AmDDhg0EBARctT3eFjQaDc8++2yTLnxXI9vWNmTbGtKlHW8LFixg06ZNZGdn4+HhgYuLC8nJyQCcOXOGhIQE8vPzcXV1Zc2aNQwaNKirTJWR6bbYRO+6jIxMx2GT7rqMjIz1kEUuI9PD6TUiT0pKYsyYMURERDB8+HBOnDjR5Hn//e9/CQ8Pp3///tx///3o9XqbsC01NZVJkyah1WobTP+1Bdu2bdvGiBEjiIqKIjo6mscffxyj0WgTtu3bt4+4uDji4uKIjo5mwYIFDSZDdZVddQghmDJlCm5ubh1nUKf143cxzU2Trc/58+eFr6+vyMrKEkajUcyaNUu8/fbbNmFbfn6+2L17t/j2229FbGxsh9tkiW2HDh0S586dE0IIUVFRIcaOHWu+pqttKysrE9XV1UIIIWpqasTs2bPF8uXLu9yuOl577TVx3333Ca1W22H29AqRt3bc/R//+IdYsGCBeXvTpk1i7NixNmFbHdu3b+80kbd1vsKiRYvEs88+a3O2VVRUiBkzZogVK1bYhF3Hjx8X48ePF8nJyR0q8l7hrrc0TbY+6enpBAcHm7dDQkIandNVtnUFbbEtOzub9evXc8MNN9iMbampqcTGxuLp6YlWq2XhwoVdbpder+f+++9n1apVKJXKDrMHelGbXKbj0el0zJo1i8cff5xhw4Z1tTlmQkJCOHLkCNnZ2VRVVfHFF190tUk899xz3HzzzQwcOLDD79UrRB4YGEhWVhYGgwEwdXakp6cTFBTU4LygoCDS0tLM26mpqY3O6SrbugJLbCspKWHmzJnceOONPPLIIzZlWx3Ozs7MnTuXjz/+uMvt2rlzJ2+99RYhISGMGzcOnU5HSEgIly5dsrpNvULkrZ0me8stt7Bx40ays7MRQvDee+8xd+5cm7CtK2itbaWlpcycOZOZM2fy9NNP25RtycnJ5hGS6upqvvzySwYPHtzldu3evZu0tDRSU1PZs2cPrq6upKam4uXlZX2jOqy1b2OcPn1ajBo1SoSHh4uhQ4eKo0ePCiGEuPfee8XXX39tPu/f//63CA0NFaGhoeKee+4x98x2tW1lZWXC399feHp6CrVaLfz9/cWTTz5pE7a9+OKLQqVSidjYWPPnxRdftAnbVq1aJaKjo8XgwYNFVFSU+POf/ywqKiq63K76pKSkdGjHmzytVUamh9Mr3HUZmd6MLHIZmR6OLHIZmR6OLHIZmR6OLHIZmR6OLHIZmR6OLHIZmR6OLHIbYPfu3QQEBJi3J02axOuvv951Bsk04N///jcBAQG4u7vz5z//2TxltTmWLVvG448/3uzxl156ib/97W/WNrN5OmyaTS/j9OnT4oYbbhAeHh7CxcVFREZGimXLlrWprIkTJ3bockghTEtWO3KWVUcANEqn1R6effZZceONN7Z4zs8//yzc3NzE8ePHRUFBgYiPj2/x/7WoqEh4eXmJS5cuNXtOSUmJ8Pb2FllZWW013SLkmtxKXH/99cTGxpKenk5hYSEbNmwgNDS0q83qUDojak5X89577zFv3jyio6Nxd3fnmWee4b333ms2yeD//vc/JkyYgKenZ7NlOjs7c+211/Lf//63o8xuSKf8lPRwLl26JACRnp7e7DnZ2dlizpw5wtPTUwQGBoqnnnrKHFjgylq1pZp85syZ4t133xVCmGoNpVIpnnjiCSGEKUCBp6enOHjwoBBCiL/85S8iKChIODs7i4EDB4rPPvtMCCFEXl6esLe3F4BwcnISTk5OYteuXUIIIbZs2SKGDx8utFqtiIqKajDXev78+eKee+4Rc+bMES4uLuLNN99s0sZ169aJwYMHCxcXFxEUFGSOkmI0GsW//vUvERoaKtzd3cWMGTPMEWWEMCXZePXVV8XIkSOFs7OzmDBhgvmdDh8+XADCwcFBODk5iZdeekkIIURycrK44YYbhKenpwgKChIvvPCCqKmpEUIIsWbNGhEbGyuef/554eXlJby9vc3v9csvvxRqtVoolUrzO2iK8PBw8fnnn5u3L168KABx8eLFJs+fPn26WLlyZZPH6vPhhx+KESNGXPU8ayCL3AoYjUYRGRkppk6dKj799FORmpra6JwpU6aIO+64Q5SUlIjU1FQRFRVl/qJaIvJXX31V/P73vxdCCPHVV1+J/v37m78shw8fFu7u7uYv+UcffSRycnKEwWAQn3zyidBoNOL8+fNN3lMIIY4cOSLc3NzETz/9JGpqasTu3buFq6urOH36tBDCJHIHBwfx/fffi5qaGlFWVtbIvo0bN4o+ffqYy8jJyRGHDh0SQgjxwQcfCD8/P3H06FFRUVEhHnnkEREVFWX+sQsODhaDBg0S58+fFxUVFeLaa68V8+fPN5fNFe56WVmZCA4OFitWrBBVVVUiLS1NREdHi//85z9CCJPIVSqV+Ne//iWqq6vF9u3bhUqlEsnJyUKI1rnrzs7OYv/+/Q32aTQa8w/plXh5eYmffvqpxTKFEOLgwYPCzs7uqudZA1nkViIrK8v8pVUoFGLgwIHixx9/FEIIkZmZKQCRnZ1tPv/jjz8W4eHhQgjLRP7rr78Kb29vIYQQixcvFm+++abw8PAQxcXFYvny5WL27NnN2hgbGys++uijJu8phBALFy4US5YsabDvjjvuEM8//7wQwiTyq4li5syZ4rnnnmvy2LRp0xq0ZysrK4WLi4vYu3evEMIk8jovRQjTj1RMTIx5+0qRf/bZZyIuLq7BPf7973+LKVOmCCFMIu/bt2+D42FhYWL9+vVCiNaL3MXFRWi1WvNHkiTx22+/NXm+SqUSBw4caLFMIYQ4e/asAJr8obQ2cpvcSvTt25fXXnuNEydOcOnSJa699lpuuukmCgoKyMzMxN7eHh8fH/P5oaGhZGZmWnyf+Ph4KisrOXHiBNu2bWPatGmMGzeO3bt3s23bNqZMmWI+d8WKFURHR6PVanFzc+P48ePk5eU1W3Zqairvvfcebm5u5s/XX3/NxYsXzedcLZhFWloa4eHhTR7LzMwkJCTEvK3RaPDz82vwHvr27Wv+t5OTEyUlJS3ae/z48Qb2Pvroo2RnZ5vPqf/OW1Pmlfj6+vL9999TVFREUVERhYWFqFQqfH19mzzf3d0dnU5n3n7wwQdxdnbG2dmZBx980Lxfp9NhZ2eHo6Njq21pK7LIO4A+ffqwdOlSysrKSElJISAggMrKSnJycsznpKamNhg2ay1KpZIJEybw6aefUlBQwMCBA5kyZQpbtmxh165dTJ48GYA9e/awdOlSPvzwQwoLCykqKiImJsbcYdRUutvAwEAeeugh8xe6qKiI0tJS3n33XfM5V0uTGxwcbE51dSUBAQGkpqaat6urq7l48WKr34MkSY3sHTp0aAN7dTpdiyGQ63O1ZwEYPXo0Z8+eNW+fP38ef3//Bj9G9YmLi+P06dPm7ffee4/S0lJKS0t57733zPtPnjzZaaG1ZZFbgcLCQp5++mlOnz5NTU0N5eXlLF++nD59+jBgwAD8/f2ZPHkyjz32GGVlZaSnp/PSSy8xf/78Nt1v8uTJvPHGG0ycOBGAKVOmsGbNGuzt7c052XU6HUqlEi8vL4xGI++//z7Hjx83l+Hj40NJSQm5ubnmfQsWLGDNmjVs376dmpoaqqqq2LdvH6dOnWq1bQsWLOCNN95g586dGI1GcnNzSUxMBOAPf/gDb7/9NidPnqSqqoqnn34af39/RowY0aqyfXx8OHfunHn7hhtuICcnh3feeYfKykpqamo4c+YMO3bsaHV5aWlpLY57P/jgg/zjH/8gJSUFnU7H448/zoMPPtjoB6eOWbNmsX379qvee9u2bR0e7LIOWeRWwM7OjgsXLnDdddeh1WoJCgpi7969fPfddzg5OQGwbt06KioqCA4OZuzYsVx//fUtTphoicmTJ6PT6cyueUxMDA4ODuZaHGDmzJnceuutDBo0CD8/P06cOMHYsWPNxyMjI7n33nuJiorCzc2NPXv2MGTIED755BOefvppvLy88Pf35+9//7tFyQhmz57N8uXLWbRoEVqtluHDh3Ps2DEA5s2bx5///GduuOEG+vbty5EjR/jmm2/MkU2vxgsvvMDixYtxd3dn2bJlODs7s3XrVn766SdCQkLw8PDgjjvuaOCut8ScOXNwdXXFy8ur2eQGo0ePZsmSJYwdO5agoCD8/Px49NFHmy3zrrvuYufOneTn5zd7TllZGZs3b+a+++5rlZ3tRY4MIyNjZV555RWKiop49dVXmzz+8ssvU1ZWxksvvdQp9sgil5Hp4cjuuoxMD0cWuYxMD0cWuYxMD0cWuYxMD0cWuYxMD0cWuYxMD0cWuYxMD0cWuYxMD0cWuYxMD0cWuYxMD+f/AZqPY1Q+q8MFAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(2.6, 3.35), layout=\"tight\")\n", "pe.plot.swrc(hydrus.model, label=\"HYDRUS\", ax=ax)\n", "pe.plot.swrc(staring.model, label=\"Staring\", ax=ax)\n", "pe.plot.swrc(rawls.model, label=\"Rawls\", ax=ax)\n", "pe.plot.swrc(clapp.model, label=\"Clapp\", ax=ax)\n", "ax.set_ylim(1e-1, 1e6)\n", "ax.set_xlim(0.0, 0.45)\n", "ax.set_xticks(np.arange(0.0, 0.5, 0.05), minor=True)\n", "ax.legend(ncol=1, loc=\"upper right\", frameon=True)\n", "# ax.set_title(\"Sand SWRC from different sources\")\n", "ax.set_ylabel(\"Pressure head |\\N{GREEK SMALL LETTER PSI}| (cm)\")\n", "ax.set_xlabel(\"Soil water content \\N{GREEK SMALL LETTER THETA} (-)\")\n", "fig.savefig(fig_path / \"dataset_swrc.png\", dpi=300, bbox_inches=\"tight\")" ] }, { "cell_type": "markdown", "id": "8119f30b", "metadata": {}, "source": [ "## Obtaining a soil model from pedotransfer function and databases " ] }, { "cell_type": "code", "execution_count": 6, "id": "ef33492e", "metadata": {}, "outputs": [], "source": [ "# Create a soil sample with easily measured properties. Note that\n", "# not all pedotransfer functions require all of these properties,\n", "# but this is a common set that covers most cases.\n", "ss = pe.SoilSample(\n", " sand_p=60.0, # sand (%)\n", " silt_p=30.0, # silt (%)\n", " clay_p=10.0, # clay (%)\n", " om_p=2.5, # organic matter (%)\n", " rho=1.5, # bulk density (g/cm3)\n", ")\n", "\n", "# Estimate van Genuchten parameters using Wösten (HYPRES)\n", "wosten: pe.Genuchten = ss.wosten()\n", "\n", "# Estimate van Genuchten parameters using Rosetta v3\n", "# Optional dependency requiring installation of `httpx`\n", "rosetta: pe.Genuchten = ss.rosetta(version=3)\n", "\n", "# Estimate Brooks-Corey parameters using Saxton\n", "saxton: pe.Brooks = ss.saxton()" ] }, { "cell_type": "code", "execution_count": 7, "id": "5f4ae504", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQ8AAAFaCAYAAAAAbNh6AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjksIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvJkbTWQAAAAlwSFlzAAAPYQAAD2EBqD+naQAAXuhJREFUeJztnXlYVNX/x18zDIvIDqLIqqKoqCDu+y6W2qrlUkpmYprWr8Xq26KWmZWZLZS0aeWS5ZZlmftuKoorbiirsogIww7DnN8fI5MkKoMDM8h5Pc88Oveee+77jjNvP+dzNoUQQiCRSCQGojS1AIlEUjuR5iGRSKqENA+JRFIlpHlIJJIqIc1DIpFUCWkeEomkSkjzkEgkVUKah0QiqRLSPCQSSZWQ5iGRSKqENA+JRFIl7inziIyMZMCAAfTt25eSkhJTy5FI7mlUphZgLBITEzlx4gRbt241tRSJpE5wz0QemzZtIi8vj/79+zNr1ixTy5FI7nnMMvKYPn0669evJyEhgejoaIKDg/Xnzp8/z/jx48nIyMDR0ZElS5YQGBhIeno6CoWCbdu2MX78eI4ePVruujKKioooKirSv9dqtWRmZuLq6opCoaiBp5NIzBshBDk5OTRu3Bil8jbxhTBDdu7cKZKSkoSvr6+Ijo4ud65fv35i8eLFQgghfv31V9GxY0chhBBffPGFWLJkiRBCiG+//Vb8/PPPFdY9c+ZMAciXfMnXHV5JSUm3/Z0qhDDfxYD8/PxYt26dPoJIT0/H39+fzMxMVCoVQgg8PDzYs2cParWa5cuXM3/+fF544QVGjx5Nly5dbqrzv5FHdnY2Pj4+nDt3DhcXl0pri3luFFl7txAz/E1+93mEqT2dGdzSvlLXKo+/gUXCckpbPI824IVK37OkpITt27fTr18/LC0tK32dKagtWmuLTqg5rTk5OTRp0oSsrCwcHR1vWc4smy23IikpCQ8PD1QqnWyFQoGPjw+JiYn079+fH3/8kb59+9KiRYsKjQPA2toaa2trIiIiiIiIoLS0FAAXFxdcXV0rrcUjqAMl+7fiX5CMysaOq5r6lb/epQFcUYCtFgy4Z0lJCba2tri6utaKL3pt0FpbdELNaS2r+07N+FplHndi4cKFlS47depUpk6dilqtvq273gq7lm0BcEg9A0EQd7XoDlfcgOp6hKLJNfi+Eom5UKvMw9vbm5SUFDQajb7ZkpiYiI+Pj8F1/TfyMBS7lu0AUCWcBqElPrO48hdb2un+lOZxV2i1WgoLC++qjpKSEiwtLcnPz68VkYcxtdrY2Nw+IXoHapV5uLu7ExISwtKlSwkLC2P16tV4eXnh7+9vcF13G3nYNmuJ0soabZ4a5+wkrip8UReW4mBjceeLVdfNoyTH4PtKdBQVFRETE4NWq73ruho2bEhsbKwRVFU/xtSqVCpp3bo11tbWVbreLM0jPDycDRs2kJqaSmhoKPb29voPLDIykrCwMObOnYuDgwOLFy+u0j3uNvJQWlpi1yoI9bGDtFKfZp+TLxczigj2sr3zxZay2XI3CCGIj49HpVLRpEmTu/rfs66i1WqJi4sjPj6eFi1aVGmYglmaR2Rk5C3PBQQEsH///ru+x91GHgD2bUNQHztI86zT7PMZQuyVSppHWc5DRh5VoqSkhNzcXJo0aYKdnZ2p5dRaPD09iYuLo6SkBCsrK4Ovl5Z9F9i37QBAw9RTAJy/Usmkqcx53BUajQagyuG2REfZ51f2eRpKnTWPiIgIWrduTadOnapch31gCAA2CSdACM6nV9I8ZM7DKMgRwXfH3X5+ddY8pk6dSkxMDIcOHapyHfVbBKKwsoLcLJyzk4jPLKZIU4kEnkpGHpLaT501D2OgtLLSd9n6X4tBK+BiRiW6bC0ddH+WFoC2aiGjxLyYPXs2EydO1L/fs2cPCoWCHTt26I9NnjyZt956i7FjxxIVFcUrr7zC+vXrb1tvVlYW8+bNqy7Zd0WdNQ9jNFsA7Nvomi4B6jNAJfMeqhuSfBrZdLkbhBAUlGir/XWnWRz9+vUrZxTbt2+nS5cuNx3r378/y5Yto2PHjnz00Uc88MADt63XnM3DLHtbagJj9LbAv0nTRmknATifXgjcoT4LK7CwgdJCKFGDlXOV71/XKdQI7vvyQrXf568pzahneescQdeuXbl8+TLJycl4eXmxY8cO3n77bT788EMAUlJSSExMpFmzZjzyyCOcP38eIQTTpk0jPDwcrVbL9OnT2bp1K1ZWVqhUKvbu3cvkyZPJyckhODgYlUpFVFQUqampTJ8+nfj4eAoKCnjwwQeZM2cOoJsPNm7cODZv3kxqaipPP/00b775ZrV8JnU28jAWDm07AmAbdwyEtvI9LmWGUXytmpRJahIrKyu6d+/O9u3bKSoqIi4ujvvvv5/k5GQKCwvZvn073bp145VXXiEgIIATJ06wbds25syZwz///MOxY8fYunUrp06d4tixY2zbtg0rKysWLVqEvb09R48eJSoqCoDx48czdepUDh48SHR0NFFRUfz66696LVlZWezfv59Dhw7x0UcfcenSpWp55jobeRiL+i3boqxnizZfTYOrF7igakFJqcDS4g6ZbCtnKEiR5nGX2KgU/DWlWY3c506UNV18fX3p3LkzoItI9u/fz44dO+jXrx8LFy7k8OHDgG7E9COPPMKWLVuYNm0aGo2GCRMm0K9fP4YOHVrh4Le8vDy2bt1KWlqa/lhubi5nz57Vvx8zZgwAbm5uNG3alLi4ODw9Pe/q+SuizprH3Y4wLUOpUuEQ1Jmsf3bQPOMo+9yac+FKES0b2dz+Qhl5GAWFQnHb5kRN0q9fP7777jt8fHzo27cvAH369GH79u1s376dJUuW3DR5s6y71NHRkZMnT7Jz5062b9/O66+/zq5du/QzyMsoy738888/2NhU/B278biFhUWVx3HciTrbbDFGV20Zjh26AdAq8zgAMamVmKylN4+su76/xDzo1KkT6enpLFu2rJx5/Pzzz6SkpNC5c2cGDhzIN998A8CVK1dYs2YNgwYN4sqVK+Tl5TF48GDmzp2Ln58fMTExODg4UFBQQHGxrhfPzs6Ofv36lUuiluVaapo6ax7GxDFEZx7uSUeAypqHk+5PGXncM1haWtKzZ09ycnJo2bIlAC1atCAnJ4eePXtiaWnJZ599xunTp2nbti39+vXjjTfeoEuXLiQlJTFo0CDatWtHmzZtaNOmDffddx8uLi6MGzeOdu3a0bGjLr+2bNkyYmNjadOmDW3btuWRRx7h6tWrNf68dbbZYkwc2ncFwPJyLDaFWZxKqcTHKpst9yQbN2686VhKSor+7w0bNmTNmjU3lQkJCdHnQv5LWaRShru7O0uXLq2wbHx8fLn3ZUnW6kBGHkbAysUN26YtAPBOOUaKWkNm3h3amdI8JLWcOmsexhokVoZDSHcAAq8dA+B02h2aLjLnIanl1FnzMGbCFMAxRNd08UvTmUdMSmXNQ0YektpJnTUPY+PYQRd52McdQVlazKk7JU3LzKMos5qVSSTVgzQPI1G/eWssXdxQFOXTOO0kp1MLKb7dDFvr66umF2XUjECJxMhI8zASCqUSpy59AGiZeogijeBM2m2Gqts01P1ZfBW0dzdQTSIxBdI8jIhz1+vmkabrHjt2qeDWha1cQaEEoZXRxz2Cn58fAQEBBAcHExAQUC2zYY8ePcrPP/9c7tisWbPuehX5qiDNw4g4de0LgHNcFMrS4tubh9ICrN10fy9Kr35xkhph5cqVHD16lG3btvH+++9z8OBBo9ZfkXnMnj3bJOZRZweJGWtuy43UbxGIpYsbJZkZeKae4KRNRzSlAtWtJslZu0NhOhSmAW2NpqNOIQSU5lf/fSxswYBl+zw9PWnZsiUJCQm4uLgwefJk0tPTUSqVzJo1i4ceeoiCggLCwsI4ceIElpaWNGzYkE2bNgHw008/8cUXX1BSUoKdnR2ff/45Hh4evP3222RnZxMcHEzXrl319+vVqxcWFhZs2rSJLVu28Omnn1JcXIxWq2XOnDkMHz7c6B9JnTUPY63ncSMKhQKnrn258ucqWqQeIsmzA2fTCwn0qFfxBTbukA0UXjHK/eskpfmwpnH13+eRy6CqX+niZ86c4erVq/Tt25dhw4YxYcIEwsPDOX/+PF27dqV9+/YcOXKErKwsYmJiAMjM1PW87d27lxUrVrBr1y6sra3ZvXs3Y8aM4dSpU7zzzjusW7eOdevW6e8VGRnJ7t27cXJyAiA0NJTRo0ejUCiIj4+na9euJCQkGH3B6DprHtWFc9c+XPlzFa3SotiKLu9xW/OA65GH5F7g8ccfR6lUcvbsWT755BNsbGw4cuQIe/fuBaB58+b07NmT3bt30717d06fPs2UKVPo06cP999/PwC//fYbx44dK7ffcmZmJgUFt2kG30BcXBxjx44lOTkZlUpFZmYmcXFx+vk2xkKah5G5Me9hoSniaHIBYzreorDePGTkUWUsbHVRQU3cpxKsXLmS4OBgtmzZwvDhw+nfv/9NZcqm4Tdt2pSYmBi2bdvGli1bmDFjBkePHkUIwfjx45k7d26VpI4aNYp58+YxYsQIQLeJe3XkRGTC1MjUb94aK3cPFMUF+Fw6zLFLBbdeUV1GHnePQqFrTlT3y8BtCgYOHMizzz7Lm2++SUhIiH5nw9jYWPbs2UPv3r1JTk5GoVDwwAMPMH/+fIQQJCUl8cADD7B06VISExMB3e5uZRPcHBwcyM7OLncve3v7cseuXbtGkyZNAFi6dCnXrlXPKGZpHkZGoVDg0mswAG0u76VIIzh5+RauXzbWozC1htRJapK33nqLPXv28Nlnn7Fy5UqCgoIYMWIE3377LT4+Ppw4cYIePXoQFBRE+/btefLJJ2nXrh29evXiww8/5OGHHyYoKIjAwEB9D8uAAQMoKiqiXbt2TJ48GYCXXnqJQYMGERwcTHp6Op9++ikjRoygffv2REdHV2kj+Eoh6jjZ2dkCEBkZGUarM/W3FWKrr0L83r216LPwnPhq95WKC6btFmKlgxAb2leq3uLiYrFu3TpRXFxsNK3VRXVqzcvLE1FRUSIvL8/oddclbvU5lv0msrOzb3u9jDyqAZdeg0ChwPbSaexy0ziUkFdxQVsv3Z/5ybouR4mkFnHPmEd8fDweHh707duXcePGmVSLpbMr9u10WdJmCXu5kFHM1YrW96jXGFCAtkiOMpXUOu4Z8wAYOnQoO3bs4McffzS1FFx76/Iewam6LrqoxAoGMllY/Zv3yK/5NSglkrvhnjKPv//+m169erFs2TJTS8GlzxAAvC7uRaEt5VDCLUZB3th0kUhqEWY5zmP69OmsX7+ehIQEoqOjCQ4O1p87f/4848ePJyMjA0dHR5YsWUJgYCAeHh6cPXsWhUJBaGgoQ4YMwdXV9aa6i4qKKCr6d7arWq0GoKSkhJKSEqM9g21gCBb2jpBzDc/U4xyy7UBhUTEWyvJdfhb1PFESRWlOAto73L9MnzF1VhfVqbU2PH9t4r/f/cp+vmZpHiNGjGDGjBn07NnzpnPh4eFMmjSJsLAwVq1aRVhYGIcOHSo39LZXr15cuHChQvN4//33mT179k3Ht2/fjq1t5QYCVRab5m2xPLKHVnFb2Ny4Pd+v24GnTflRgoGaYvyBuJjdnDrnXal6N2/ebFSd1Ul1aC2bByIxDnv37i1nGPn5lZsrpBDCfNP8fn5+rFu3Th95pKen4+/vT2ZmJiqVCiEEHh4e7Nmzh0aNGmFnZ4cQgiFDhrBkyRI8PDxuqrOiyMPb25uUlJQKzeZuSP/jF87+35PkN2rO/NF/8FiwA890dylXRhkbicWJN9A2foDSLt/ftr6SkhI2b97MoEGDsLS0NKpWY1OdWvPz84mNjaVVq1ZGN/y6RH5+PqdPn8bf37/c56hWq3FzcyM7OxsHB4dbXm+WkcetSEpKwsPDQ7+LlkKhwMfHh8TERC5evMj//vc/LC0tGTFiRIXGAWBtbY21tfVNs2otLS2N/iV37z+UcyoVtqnncbkWz/6E5kzp8597ODQFQFmQgLKS968OrdVFdWg112dfs2YN7733HqWlpRQWFtK4cWO2bNlS4baRd+Lo0aOcOXOGUaNGVYPS8vz336iyn2+tMo/bMXjwYAYPHlzp8tUxq/a/WDo64dS1L9f2bKFl3Db2OfuReK0YH2erfwvZ6cyDnIu6sR4GDoOu6wgh0BZU/5R8ZT1b/ZyUikhJSWHSpEkcPnwYX19fAI4cOXLba27H0aNHWbduXY2YR1WpVeZR1rzQaDT6ZktiYmKVht9Wx3oeFdFg0ANc27OFkKTt7AuZwL6Lefh0uNE8/AAFaHKg6Mq/810klUJbkM/O1vbVfp8+MTlY2N56Sn5aWhoWFha4uPzbLA0JCQHg5ZdfZufOnZSUlODg4MA333xDQEAAZ8+eZcCAAezatYumTZsyf/58tmzZwo8//njTuh2LFi3i77//5vXXX0ej0eDs7MxXX31F69at2bFjB8899xy9e/dm7969aDQafvjhB/0Oc9VFreqqdXd3JyQkRL9b1urVq/Hy8sLf39/guoy99cKtcB2gW4TFOT6KegWZ7L2YW76AhQ3YXk+U5lyoVi2S6qNdu3b07NkTX19fHn74YT766CMuXboEwKuvvsqhQ4c4evQoU6ZM4fnnnwcgICCAjz76iMcee4wdO3YQERHBTz/9hLu7O++88w79+vXj6NGjLFq0iPT0dMaMGcMPP/zA8ePHmTRpEiNGjNBvfH3mzBnGjx/PsWPHmDZtGm+88Ua1P7NZRh7h4eFs2LCB1NRUQkNDsbe3JzY2FtAtfBIWFsbcuXNxcHDQz1Y0lJqKPOp5+WLXOpjcmKM0v7iTE/Ue5mqeBtf6N3z09k0hPxFyL0CDbtWq515DWc+WPjE5NXKf255XKlm9ejVnzpxh586d/PXXX7z33ntERUVx8OBBPv/8c3JyctBqtfpFfwBGjx7N9u3bCQ0NZevWrTRo0KDC+g8cOEDbtm1p21a34tzYsWOZOnWq3qD8/f31639069aN+fPnG+Oxb4tZmkdkZOQtzwUEBLB///67vkdN5DzKaDD4QXJjjtI5aTPHAx9mV2wuDwc5/VvArhmk7ZCRRxVQKBS3bU7UNC1btqRly5aEh4czZMgQli9fzsKFCzl06BDNmjXj+PHj9O7dW19eo9Fw8uRJXFxc9EZQFWxsbPR/t7CwQKO5w3anRqBWNVtqK+5DRwLgcX4X1oVqtp/7z/+Uds10f+ZK86itXLp0Sb9aGOjW1IiLi8PR0RFLS0s8PDwQQvDFF1+Uu+61114jICCA3bt38/LLL+sj7P+u29G1a1dOnDjByZMnAfj555/x9PTE09OzBp6uYswy8qgJaqrZAroFguoHtCHv7ElaXtjCcZtHuJKroYHd9Y/f/rp5yMij1qLRaHjnnXeIi4vD1tYWjUbD+PHjef7557l48SKBgYG4urry0EMP6a/5448/2LhxIwcPHsTW1pYFCxbw2GOPsW/fPgYMGMD8+fNp164d3bt3Z9GiRSxbtoxx48bpE6a//vprlXtzjEJ1rBNQm6iO9Twq4uKn74itvgqxbHA/0WfhOfHLkcx/T+Zc1K3r8aubEKUlt6xDruehQ67nYRzkeh61BPehjwHQMHYP9Qqusf3cDb0u9X11S91pi2XTRVJrqLPmERERQevWrenUqVON3K9+swDsWgWhKNXQ8sIWYlILSVVfn0+gUILD9ZWts0/ViB6J5G6ps+ZRU+M8bsR9mC5x2iXhbwC2nr0hcerYWvdnVkyN6antaLW32Uhcckfu9vOrswlTU+A+7HEufvQmDWL3Ypebzt+nLRnT0VmX9HIM1BVSnzatyFqAtbU1CoWClJQUPDw8qjR3pK6j1WpJSUlBoVBUeTMoaR41iK1vMxxCuqE+sp/25/5gt90EzqQV0aqRzb+Rh2y23BELCwv8/f2JjY3Vr8ciMRyFQoG/vz8WFhZVur7OmkdNdtXeiMej41Af2U/n87+xu/1TbIxRXzeP65FHbjyU5IKlXY3qqm04ODgQFBRUbnmFqlBSUsLevXvp0aOH2c7WLcPYWq2tratsHFCHzaMmR5jeiPuwxzk/+wXqXz5Doyun2XquDVN6u2Ft46ZbELngMmQdhwbda0xTbcXCwuKu1/MoW0XL1ta2VpiHOWmVjcUaxtLRCbdBDwLQ9fxv5BZp2R93fWsGF90sTDIPm0idRFJ5pHmYgEaPPglA4Jk/UJaW8Nep6+12vXkcMZEyiaTySPMwAS69Q7Fya4iFOgP/+F0cTMjXjfmQkYekFlFnzaOmB4ndiFKlouFDYwDoe2EtAvjjZDY4t9cVyEuAQrkJlMS8qbPmYYpBYjfSeNQzADQ6vRWHnBQ2nFRTYuEI9tcXNroWbRJdEkllqbPmYWrq+7fEqWtf0GrpfnY11wpK2X0hF1yuLx2X8Y9J9Ukkd0KahwnxHDMJgJCTv6LQavjteDa4XV9J7Mo+EyqTSO6MNA8T0iD0YSxdG6C6lkpA3A6OXSog2aqz7mRmFGgKbl+BRGJCpHmYEKW1NR4jwgAYcP4XAH4956rb/FpbrDMQicRMkeZhYhpfb7q4ntmFU1YSG8/kUOzSQ3dSNl0kZkydNQ9TdtXeiK1vM1x6DQYhGBy7kiKN4Ejh9fEeV3abVJtEcjvqrHmYuqv2RrzCngMg4PBKrIpz+Sn++iS5jAOgyTOhMonk1tRZ8zAnXPvdj23TFijy1fSMXcepPC/yLb11eY/0XaaWJ5FUiDQPM0ChVOL1lG4XsW5Hf0Kh1bK/8Pqs2pTNJlQmkdwag80jLy+PAwcOsHHjRg4cOEBengyrjYHHo+NQOTpjkRZP28QdbMnuqjuRulm3AbZEYmZU2jx+//13Bg8ejLOzM4MGDWLixIkMGjQIFxcXBg0axG+//VadOu95LGzr6weNDY75iejijpQIS8hLhJxzJlYnkdxMpcyjV69eLFy4kLFjx3Lx4kXUajXJycmo1WouXLjAk08+yeeff15uGz2J4XiNfw6FSoXt2f24p1/gaPH1XpeUTaYVJpFUQKVWEvvggw/o3r3ila28vLwYN24c48aNM8oesnUZ60aeNHxgNKlrfuKRM9+zz7c3nawPQPJ6CJhmankSSTkqFXncyjj+S7dupt/hfeHChQwcONDUMqqM77OvAuAQtYHTKU3RCgVcPQj5Vd8EWSKpDqq0hunRo0c5cuQIubm55Y5Pnz7dKKKqSklJCUePHjWphrulfvPWNBjyCFc2rmHQ6XWcbBpEO6ujuuijyURTy5NI9Bjc2/Lmm2/SrVs3IiMjWbt2rf61bt26apBnGD/99BOjR482tYy7xnfq6wA4H1zL/hRd3iPvwhpTSpJIbsLgyOOrr77iyJEjtGrVqjr0ALoIZv369SQkJBAdHU1wcLD+3Pnz5xk/fjwZGRk4OjqyZMkSAgMD0Wq1/P3336xcuZKPPvrolnUXFRWVW66/bN+PspWpzYF6Ldvh3Hsw13Ztwv14MvhC/ZyD5KsTAMxG5+0o02juWmuLTqg5rZWt32DzcHFxoUmTJgYLMoQRI0YwY8YMevbsedO58PBwJk2aRFhYGKtWrSIsLIxDhw6xZs0aHnjggTvW/f777zN79uybjm/fvv2ul/E3JhbBfbDdtQmHfb8R3S2E9m6nObz5E6g3mM2ba8/AsdqitbbohOrXmp+fX6lyCiEMG4H022+/sX79el5//XXc3d3LnXNwcDCkqjvi5+fHunXr9JFHeno6/v7+ZGZmolKpEELg4eHBnj17WLlyJTt27MDCwoKoqCjmzZvHxIk35wgqijy8vb1JSUnB1dXVqPrvlmNjBqI+tBv69KL/0KNc1vpxwHo2oaGDzWLfjttRUlLC5s2bGTRokFlrrS06oea0qtVq3NzcyM7Ovu1v2uDIw8nJiU2bNrFkyRL9MSEECoWi2ndfS0pKwsPDA5VKJ1uhUODj40NiYiJvvPEGb7zxBgADBw6s0DhAt0uWtbX1TTvGWVpamt2Xp9lL7xA9qh+KvfvJ6upAY9d41PkpZqn1VtQWrbVFJ1S/1srWbXDCdMKECUyYMIGTJ09y8eJFLl68SFxcHBcvXjRYZHWxZcuWO5Yxp1m1t8K5ax+cew1CaDSc2+kMQCPNPoo1cnd4iekxOPK4evUqs2bN0u3sXsOUNS80Go2+2ZKYmIiPj4/BdZlqr1pDafbyHKJ2byb3YBx5PWzo2WAr64+mMrqbr6mlSeo4Bkceo0aNYv369dWh5Y64u7sTEhLC0qVLAVi9ejVeXl74+/sbXFdtiDwAHII64Tb4IdBqOb9FRX1lHldiVnI1T2NqaZI6jsGRx6VLl3j88cfp0KEDDRs2LHduzRrjjEUIDw9nw4YNpKamEhoair29PbGxsQBERkYSFhbG3LlzcXBwYPHixUa5pznT9KV3yNj8G5nH1OT0sWV4wxV8s+dJXgv1MLU0SR3GYPPo2LEjHTt2rA4teiIjI295LiAgwChzaGpLswXALqANDR8cQ9q6ZZzfWELIxIukXdzO6dRHaNXIxtTyJHUUg81j5syZ1aGjxpk6dSpTp05FrVbj6Ohoajl3pMn/zSJ9wy9knSvh6lkVj/r+zGc7exHxmBdKE+SfJBKDcx4//vgjx44dK3fs6NGj+jxEbcFcFkCuLLa+zWg8bioAsX8U0VW1E/WVc2w5k2NiZZK6isHmMXPmTDw9Pcsd8/T05K233jKaqJqgtiRMb8Rnyuto69uTl6Yl9VAxY+yW8NXuDHIKzb/pJbn3MNg8MjMzcXNzK3fMzc2Nq1evGk2UpGJUDk4U3zcKgLhNxQxgA1bFyUTuzTCxMkldxGDzaNKkCTt37ix3bPfu3fj6ynEHNUFJ98HUaxpASZ4geXseo+r/yB8n1Ry/JLemlNQsBpvHq6++yqOPPsqcOXNYuXIl7733HiNGjODVV1+tDn3VRm3LeeixUNH0tXkAJO0pYUD+GlyVV1iwLZ2SUrlQsqTmMNg8Ro8ezeLFi9m3bx/vvPMOe/fu5ZtvvuGJJ56oDn3VRm3MeZTh3Pc+XPsOQZTChXU5THL8hvjMYlYcvmZqaZI6RJVWEhs+fDjDhw83thZJJVEoFDSf/TnXBgVy7XwxQWd/xcdjDD8dVNCrWX2auFqbWqKkDlCpyCMmJqZSlZ0+ffquxNQktbbZch1b32b4TtGtOBb7ewGvOH1GSang/U1paGTzRVIDVMo8Ro8ezciRI/nzzz9v2uQpPz+fv/76i5EjR9aqJQBrc7OlDJ/Jr1LP25titcBm0yY61z/OufQilkdlmlqapA5QKfOIjo5m6NChzJw5EycnJ/z8/GjXrh1+fn44Ojry9ttvc//993PkyJHq1iu5AQsbG1q8swiA5D0lvFg6FwVafjiYyfkrRXe4WiK5OyqV81AqlYSFhREWFkZKSgpHjx7l2rVrODs7ExwcjIeHnKBlKlz73UeDQfdzZfOfJC89xgv/t4FPkoYzb1Mqi0b5YGkhh65LqgeDE6YeHh7SLMyMFu9Gkrm3OeqkQoJ3vIt3295cyHDku/1XmdzT7c4VSCRVwOCu2nuF2p4wvRHrRp40f2shAHF/ZjOn/nwAfj58jagEuRG5pHqos+ZxLyRMb8Rj1DM4d+6IVgNXv1nN1OYnAZi7KY1r+XLhIInxqbPmca+hUCho+fFKlNYqsuJK6bbzOVq6lpCZX8q8TWloDVskXyK5I9I87iHqeTeh2Yy5AFxYm8K7dh9iZaHgQEI+q49mmVac5J5Dmsc9htdTL+IY1JbSYkj8dBkz258CIHJPBqdTC02sTnIvUaneliZNmhi0WvoLL7xg8k2v6yoKpZLWn6/l4OBWZMeV0PTnZxgyeBsbL1ozc0MKX4/2xsm2SrMSJJJyVOpbtGPHDoMqdXJyqoIUibGo59OUFrM/4/SrzxK3IYPwDs8T4/QNiVka3t2YyocPeWKhlOM/JHdHpcxDrtVR+2j02CQyNq/mypYtnI3YyYLPVvLEkZEcTipg8T9Xmdhdjv+Q3B0Gx6/vvPNOufd2dna0a9eOgQMHGk1UTVCbVk+vCgqFgpYfrSB7YHPyr2Rx9bN3ePeVbryy25Olh67RsqENPZvZmVqmpBZjcMI0Ojq63Gvjxo2MHDmSp59+ujr0VRv32jiPirB0dqX1Jz8DcGl/EU3+Gsu4tro5L+9vSiMhs9iU8iS1HIMjj7Vr1950LD09naCgIM6dO0eLFi2MIkxiHFx6D8b7qakkLY4gZmkKo5tN4oTHj0SnaHl9/WW+etwbx3oWppYpqYUYbB4TJkyo8LilpSWTJk2iadOmAHz//fd3p0xiNJq9/jHZUXtQnzjGqYiDvD9nLmF5b3I5u4S3N6Qw/2FPOYFOYjAGN1scHR0rfFlYWFC/fn39e4n5oLSyIvDLtagc7FEnaUn8ZhkR7dZga6Xk2KUCPtmejpAjUCUGYnDk8cknn9x07MyZM/zwww9ERETg5+dnDF0SI1PP24/WC5ZyfOKDJO8pwanJ23z6YEPC93Xlz1NqfF2seDzE2dQyJbUIgyOP9u3bExISon8FBAQQHBzMtGnTpHGYOW4Dh+Mz6WUATv9aiOfBacwK0S0xuWh3Bnsu5JpSnqSWYXDk8cILL5R7X9ZV27x5c2NpqhIHDx7k//7v/xBC0L9/f+bMmWNSPeZK01feI/vIPrKj9nHih1y6PT+Zp1su5bszPrzzVyofP+JJ28b1TC1TUgsw2DzGjx9fHTrumvbt27N3714ABgwYgFqtxsHBwcSqzA+lpSVtvvyVQ8M7kZd2mdPLrjJ24jOk+/zA74mN+N/6y3w+0gs/uQK75A5Uyjx+/PFHgyoNDg6mXbt2VRJUVSwtLQEoLS2lcePG2Nra1uj9axPW7h60XbSaI4/1IeNUMQl/XubFYU+T0+hbdqR68Mq6y0Q85oW7vaWppUrMmEqZx/bt2w2q1MnJ6a7MY/r06axfv56EhASio6MJDg7Wnzt//jzjx48nIyMDR0dHlixZQmBgIADLly9n1qxZhIaGolLJyV+3w7F9FwLe+4ozM54mbnMxdo2TeTvkGfJcIjmU6cmM3y7z+Qgv7G3kGBBJxVTqF7Z48eLq1lGOESNGMGPGDHr27HnTufDwcCZNmkRYWBirVq0iLCxMP0p0zJgxjBo1ipEjR3LixAnatm170/VFRUUUFf27srharQagpKSEkpKSanoi41Cmz1g6Gzz8BOoTh7n805fE/FxMB9dLzPMNZ3pxJKeuevL6+ku8P6whNpaGr9xgbK3VRW3RCTWntbL1K4QZd/D7+fmxbt06feSRnp6Ov78/mZmZqFQqhBB4eHiwZ88evL29sbbWtdPDwsJ44403Kkzizpo1i9mzZ990fPny5XWzqVOqod6Xs1HFnsTKWUWn56zR2LvxwpUIzmua42OTy4PuSaiUZvs1kRiZ/Px8xowZQ3Z29m3zhrUqtk9KSsLDw0PfJFEoFPj4+JCYmEh0dDQRERFotVp69+59y96f119/nRdffFH/Xq1W4+3tTb9+/XB1da2R56gqJSUlbN68mUGDBulzPEapt2cPjj3Wm4L4WI7/aEXIM1f5yuNZXrzyKccLAzmkCOHtUHeDRqFWl1ZjU1t0Qs1pLYvG70StMo/bMXLkSEaOHHnHctbW1lhbW980q9bS0tLsvzxlGFurpXsjghZv4PAj3clJuMqp1W60ffwan7g+y/+uzuefhM7M23qVt+9rhMrAdUBqy+daW3RC9WutbN21ahlCb29vUlJS0Gh0q4ELIUhMTMTHx8fguurCrFpDsG3SnLZfr0VhZUVGdAaxO9yx0Obzvsvz9Lfdyq7YXOZtSqNUK5svEh21yjzc3d0JCQlh6dKlAKxevRovLy/8/f0Nrute2rfFWDh16kmrj3TJ8aSNF0mOaYVSlPCW42s8breULWfVzN+aLldilwCVbLbcaibtfzHWTNrw8HA2bNhAamoqoaGh2NvbExsbC0BkZCRhYWHMnTsXBweHKvcETZ06lalTp6JWq+VEvhto9OBoCpMucnH+W5z74RBWL96He8PdPGu/EC+LBBbGvIpWCF4Z2NDgJozk3qJSkceNs2eVSiXLly8nNTUVa2tr0tLSWLFiBRYWxhsPEBkZSXJyMhqNhrS0NL1xAAQEBLB//37OnTtHVFRUhd2xlUFGHrfGd+r/aDxmEgjBqc82k6kZBygYbruWj1ymsffMJd77OxVNqYxA6jKVijxunEk7YsQIVq1axbBhw/THNmzYUOvW75CRx61RKBQEvBuBJvsa6Rt+5cScJbT/7F0c1B8SwiG+cgvjtQsLmfWX4O0hjbBS1arWr8RIGPyvvmnTJu6///5yx4YMGcLmzZuNJkpiehQWFrRe8CPOPQdSmp/H0ZfnkOf3Fdh64a1KJNJtHKXJf/HWHykUabSmlisxAQabh5+f301RxpIlS2rdCuuy2XJnlNbWtI1cg0NQZzRZmRyd8hwFLX8A187YKXN53+VF2lydz//WJVFQIg2krmGweXzxxRe88sortGnThqFDh9KmTRtefvllIiIiqkNftSG7aiuHqr4d7Rb/ga1/K4pSL3H06dEUtfwa/MMBeNL+e0YXTOTtNSfIKrg3V6KXVIzBg8R69+5NXFwcf/zxB5cvX6Zx48YMHToUZ2e5CtW9ipWLG8E//c2Rx/pQkHCBI08MIWTlDqxdO6E9NI2O1gfxLh3JwlXzmfzgfTRyqB2DrSR3R5VGmDo5OfHEE08YW0uNcq/v22JsbDy8aL9iG9GP96Ug7jzRowfQfsU2rAdto2T3EzTMv8CbyjBWrnuWbve9gbeTNJB7nSqZx/79+9mxYwcZGRnlFs5dsGCB0YRVN7K3xXDqefnS/udtHHm8L/kXznB07EDar9iGVegOCvc/h03qb4y1+YJjm/7hTKcvTS1XUs0YnPOIiIhgwIABHDx4kIiICOLi4oiMjCQ1NbU69EnMjHreTXQRRyNP8s7HED12IMXqImx6/UBB8OcUUY8gyyh8jwyC/FOmliupRgw2j4ULF/LXX3+xdu1a6tWrx9q1a/n111/10+El9z62vs10EUfDxuSdPUn06P4UXUmjXotxMHAXyYpAHJTZPGjxAYlbp0GJXFj5XsRg80hLS6NPnz66i5VKhBDcd999rF+/3ujiqhPZVXt32DZpTvvlW3UGcu4URx7rTUFyAtYuLWj08HYOWOu2H22mXkHWb13RpO0xsWKJsTHYPBo1asTly5cB3ZiPHTt2cOrUKZTK2jXKUHbV3j31mwXQ4ddd2Hj5URAfy5HHepMfdx6Vypr2983j++L3SNU0wkmbhHLnMEqiZoAm39SyJUbC4F/8s88+y4EDBwB48cUXGTx4MO3bt2fq1KlGFycxf+r5NCXk113YNg2g6HISR0b2JvfsSRQKBa723pwN2sJfBQ+hRGB5MRLNxh6QccDUsiVGwGDz+L//+z8efvhhAMaOHcvFixc5fvw4s2bNMrY2SS3BxsOLkF92YteyHcUZaRx5vC85Jw4D0L1FI5qERvJu3udcKXVHlX8RsW0IHH1D5kJqOVVqa2g0Gvbu3cvKlSvx9vbGx8eHvLw8Y2uT1CKs3Nxp//M2/VD2E0+GYnH2OAAtG9kw6dHRzCpdw8b8YSjQwrkvYGNnSF4Pcn2QWonB5nHmzBlatWrF448/ztNP65JiW7duZeLEiUYXV53IhKnxsXRyIXjZZpy69qU0L4d6ke+Svv5nABo6WPLRY4HscP6QVzMXclnjCQWXYN+TsOcxyI0zsXqJoRhsHlOmTOH5558nOTlZv9Zh37592bOndmXTZcK0elDZ2RP8w1+43T8CRamGsy+NJ/Hr+QghsLVS8t7wxjRrO5ynrqzkx5yn0aCClE3wd1eI+QhKi+58E4lZYLB5HD16lClTpgC6dR8AHBwcyMnJMa4ySa1FaW1Ny09+orjPcABi587g/Dv/h9BqsVAqmNTDjdfu82N54RSeSv+Zk6WdobQQTs6BTT0gbYdpH0BSKQw2j4YNGxIfH1/u2Llz5/Dy8jKWJsk9gEKppOiRCTR5/QMAkhd/xqlpoyktLASgXwt7vnjMm5L6/jyXHsE89XsUqdwh5zzsfBD+eRoK5Khlc6ZKXbWPPvoov//+O6WlpWzatIknn3yS5557rjr0SWo5XhNeIPCz5SgsLUnf8CtHx4VSnJkBgH8DayJH+RDiZcvGvFAeSf6Fk/WfRKCExFWwsROcjwStnLxojhhsHtOnT+eZZ57h9ddfp7S0lBdeeIEnn3ySyZMnV4c+yT1AwwdGEfzDRizsHcg+uJuoh7qSdz4GAMd6Fnz4sCcjgp3IE3Y8F/s8Hyt/RuMYDCVqiJ4BW/tB5mHTPoTkJqrUVTtlyhROnjxJbm4uMTExMuqQ3BHn7v3ouGYfNj5NKUy8SNQj3bm6YyMAKqWC5/o04O37GlHPUsEfl5oyKvFbkpvOA0tHuHYMtgyAwy9C8TUTP4mkjCqZR15eHtHR0ezatavcqzYhu2prnvrNW9Nx3T84de5NaY6aYxOGkbT4M/2yDv1b2LNolA9+rlZk5MO4fQNZ67EZ4fMYIODCd/BXJ4j/WY4NMQMMNo9ffvkFDw8PunXrxoMPPqh/PfTQQ9Ugr/qQXbWmwcrFjeClm/B47CnQajk/+wXOvvEs2us7s/u6WPHV494MammPVsCnB5S8eWU2+d1/A/sWUHQFDobDjmGQfcbET1O3Mdg8ZsyYQUREBIWFhVy7dk3/yszMrA59knsQpZUVLT/4Fv8354NCweXlXxM9ZgBF6brelXqWSv43uCEv9ddtrr33Yh5Pb/blXPut0PZtsKgHV/bounWPz5KT7UyEweaRnZ1d65cglJgehUKBz8QXaffdel0i9dAeDg3rQPbhffrzw9s6EvGYFx4OKlLUGqauSmdN0dOI0H+g8X0gNHDmE9jYBS7/ZeInqnsYbB5PPPEEa9asqQ4tkjqIW/+hdFp/iPotAilOT+HIqH4k//SlPg/Swt2Gr0f70KNpfUpKBZ/tvMIb26zI6rAMeiwHW2/IT4Q9o2DPaMhLNPET1R0qtYbpww8/rB9NqtFoGDt2LB9//DGNGjUqV06aiqQq2DZpToe1+zkzYyLpG37h3FvPoT4WRcCcCCxs6mFvY8GcYR6sO57Nl7sz2BeXx9PLEngjtB8hQw5AzIdw9gu4/KdudGrrGdBiKlhYmfrR7mkqZR7BwcHl3nfo0KE6tEjqMKr6dgR+sQKHoI7EznuN1FVLyD19jDYRK7H180ehUPBwkBPtGtfjnY2pJGQW89KaS4zp6MxTXWeh8h0FR16CK3vhxCxI+BlCPgb3nqZ+tHuWSpnHzJkzq1uHRKLLg0x6Gbs2IZx6bhS5p6I5NKwDLed9Q8NhjwHQrIE1kaO8idh1hd9PqlkWdY0jyQW8NcSfxn03QMIKOPYmqM/AjqHgOxqC5oCNm4mf7t6jdq0deBv27NlD165d6d69Ox9//LGp5UjuApfu/em04QiOnXpSmpvDqedGceZ/kyktLADAxlLJSwMaMvv+RthZKzmdWsjE5YlsPZcLfmNgSBQ0fQpQ6MxkY0eIWyrHhhiZe8Y8mjZtyq5du9i3bx9//PEH+fmy+642U7bJlO9z/9N35x5+qBt5F87qy/Rpbs93Y3xo29iG/GIt725MZd6mVPIVTtBxIQzYDI5tdKNSD03VjQ1RnzfZM91r3DPm0bhxY6ysdAkyCwuLWrcgs+RmlCoVzV6eQ/CPG7F0bUDumeNEDe9I6rpl+jINHSz55FEvwrq4oFTAxtM5TFqRyNm0QnDtBIN2QLt3bhgb0h1OzZPrhhiBKu0YV91Mnz6d9evXk5CQQHR0dLmE7fnz5xk/fjwZGRk4OjqyZMkSAgMD9ec3b95Ms2bNsLGxqbDuoqIiior+/eKo1WoASkpKKLk+ytFcKdNn7jrBuFrtu/YlZP0hzrwURvY/O4h54UmubPsT/5kLUTk4ATC2gwPtPKx4f8sVkrNKmPpLEhO6OvNokAPKZlOg0VAsjs1AmbYVTr2PSPiV0vYfU+LY2Wg6q5ua+vevbP0KIQxvCG7fvp3ly5eTmprK77//TlRUFDk5OfTr189goRWxa9cumjZtSs+ePVm3bl058+jfvz/jxo0jLCyMVatW8cEHH+iHmCcnJzNu3DjWr1+PnZ1dhXXPmjWL2bNn33R8+fLl2NraGkW/pJrQlmK1aRVWf/+CQqtF69yAwiemU+rfRl+ksFTJ5quNic13AMDXJpfQBpeob1EKQtBYe5C2mqXYkA1AnLIfMapRaBT1TPJI5kh+fj5jxowhOzsbBweHW5Yz2Dy+/fZbZs2aRVhYGJ9//jnZ2dkcO3aM5557jt27d9+18Bvx8/MrZx7p6en4+/uTmZmJSqVCCIGHhwd79uzB29uboUOHEhERQUBAwC3rrCjy8Pb2JiUlBVdXV6PqNzYlJSVs3ryZQYMG6ZeANFeqU6s6+gBnXw6jMPEiKBR4TXwR3+dnory+a6EQgg0xuXy1J5PiUoFTPSUzBjSgk891gyjORnnqHSzifwAgHzeUXb/CwsM4//lVFzX1769Wq3Fzc7ujeSAMpHnz5uLUqVNCCCGcnJyEEEIUFxcLV1dXQ6u6I76+viI6Olr/PioqSrRo0aJcmU6dOomtW7eK77//Xnh4eIg+ffqIPn36iOTk5NvW/cUXX4hWrVqJFi1aCEBkZGQYXb+xKS4uFuvWrRPFxcWmlnJHqltrSW6OiJnxtNjqqxBbfRXiwH3tRe65U+XKXMwoFE/9FC/6LDwn+iw8JyJ2pYtijfbfAmk7hfb3NkKsdNC9Dr8kRHFOteg1BjX175+dnS0AkZ2dfdtyBmcVr169SuvWrYF/1zBVKBT6v5uKp556isuXL7Njxw527NiBp6fnbcvLWbW1G1V9O1p98C1tF63G0tmV3JijHBrWkcRvP0GU6lYea+JqzVejvHmonSMAvxzJYuovSSRfK9ZV4t4bTf+dxCn7697HfqObbHdVficqg8HmERQUxOrVq8sdW79+PSEhIUYTdSvKmhcajQbQhaeJiYn4+PgYXJdcz+PeoMGQh+n893FceoeiLSokds5LHHmsD/kXzwFgrVLyQj935gzzwMFGybn0IiauSGRjjFo3f8bSnuOWYWh6rNbNk8mLh21D4PQCEFrTPpyZY7B5zJ8/n/DwcB599FF9YmXKlCl8+OGH1aGvHO7u7oSEhLB06VIAVq9ejZeXF/7+/gbXJSOPewdrdw+CfviTgLmLsKhvR/bhfRy8L5jEbz7WRyE9m9nx3Vhfgr3qUVgimLc5jff+TiOvWGcQwr0PDN4L3o/oZuuemA27HpaLMN8Gg80jJCSEU6dO0a1bNyZOnEhQUBBHjhwhKCjIaKLCw8Px8vIiOTmZ0NDQcuYQGRlJZGQkLVq0YN68eSxevLhK95CRx72FQqHAc8wkOv99Audeg3RRyHuvcHhEL/JidYsGNbBT8fHDnjzdzRWlAraczWHKL5dJLbrerW/lCF2/h45fgIWtbpLdpu6QutV0D2bGGNTbotFo8Pb2Ji4u7pbjKGobarUaR0dHMjIyakVvy59//sn9999fK3pbTKVVCEHKyu84P+clSnNzUFpZ0+Sld/B++v9QqnRDm05eLuDdjamk5WhQIpjY3YVRHV1RluXu1Gd12z9knQCFEtrOhIDnwYS5vZr6TMt+E3fqbTEo8lCpVNjZ2VFcXHzXAiWS6kKhUNB41ES6/H1ClwspLuLC+69y+KGuqK9vwN2mcT2+HeNDr6a2aFHw9b5rvPbbZa7l6/JpOATAgC3Q5Eld7uP4TPjnKdDIPZnLMLjZ8sYbb/DUU09x4sQJsrKyUKvV+ldtQjZb7n1sPH0I+uFPWn74LSoHJ3JOHiHqwS6cm/0Cmtwc7G0seCu0AQNcL2NloeBgQj5PL0skKvH6vCgLG+j4uW5qv0IFSWthW6jMg1zHYPOYMGECa9euJSgoCFdXV5ydnXFycsLZ2bk69FUbMmFaN1AoFDR+bAJdt56m4YNjQKslefFnHBgUyJVNv6FQKGhnn0XESA/8XK3IzC/llbWXiNyTgaZU6Jop/hOh7x9g3UDXjNk6SNesqeMYbB5xcXH618WLF7l48aL+7xKJuWLVoCGBny4l6MeN2Pg0pSglmROTHiZmykgU1zLwc7EicpQ3D7R1RAArDl/j+dXJpOdcn+fRoJtulq5dM92yh9sGQ8Y/Jn0mU2Owefj6+t7yVZuQzZa6iWvvwXTZdALfKa+hUKm4unk99d+fRvJ3n2ApSnmxvzuz729EfSslp1IKeWZFElEJ1/Mcdk10BuLaCYqzYOfDkL7HpM9jSgye2zJhwoRbnvv+++/vWlBNI3tbqofaoDX37EnOvB6O+sh+AGybtaTFrE9x6TWIy9klzPozhXPpRSiA8V1ceLKzCxZKhW6rh31P6LpwLWyh10pw713temt1bwuAo6NjuVd+fj6rVq0y2y+IRHIr7ALa0G7FNgpGT8XSpQH5F85w9MlQTjw7Euecy3w+0ovhbRwQwJIDmbz222WyCkpBZatbub3RQCjNh92PwZX9pn6cGsdg8/jkk0/KvX7++Wd+++03cnJyqkOfRFKtKJRKNF0H0nHzSbzCpqGwsODKX6s5MLA1lyPm8EIPB14f3BBrlYJDifk8szyRmNRCXU9Mj2XQaBCUFsDeUXUuiWqU5bb69u3Lhg0bjFFVjSFzHpIbUTk40WLWp3TacASnLn3QFhYQ98ksDgxqQ0jSdr58zAtvJ0uu5Gp4flUyG2PUOgPp/uO/OZBdI+pUN67B5nHjuA61Wk1qaioffvjhTXu4mDuyq1ZSEXYt29L+520Efr4C60aeFCbFcWLSw2RPu58F7TL0m0/N25xGxK4raJT1oOfKf3th9o6F0roxiNJg8ygb01H2aty4MYsWLWLRokXVoU8iqXEUCgUNhz9Ol62n8Z3yGkora7L+2cGJR7owZvvrTGhWCMCv0Vm89ttlcoQT9F4Flo6QGQXH3jDtA9QQVRrnUTa+4+LFi1y5coW4uDijLUEokZgLqvp2NJsxly7bzugGmAlB2pqf8Hm5K2+nf4uDNo+oxHwm/5xEUokXdPlad2Hs15C4+vaV3wPc9TiPxMRETpw4UR3aJBKzoJ6XL4GfLqXjun9w7NQTbWEBLPuIF3+6jz7nV3H5WiFTf0niBH2g1Uu6i468CAVpphVezRhsHg8++CB79ugGxkRERNC9e3e6detW65otMmEqMRSH4M6E/LKTNotWUc/PH21mGn3+eINpPz9Mo1NbeHF1MrusnwPnIF0C9ciL9/RGUwabx/79++nSpQug+wFu2bKFf/75p9bt0iYTppKqoFAocB/yCF02naT525+gcnTGKe0so9Y/yxPLR/H913+y2X4eKC3h0h9w6XdTS642DDaPoqIiLC0tuXTpEpmZmfTo0YM2bdqQlnZvh2gSyY0orazwnvA83XbF4vvsqyht6uGdEs24VeOIf/U1Dmc+oCt4/O17doMpg82jdevWvP/++7z77ruEhoYCui0R6tevb3RxEom5Y+noTLNX36fbrlg8n3gWYaGiWcIesucu4ehSLXkXY3UJ1HsQg83jyy+/5I8//uDs2bPMmjULgL///pvBgwcbW5tEUmuwdvcgYE4E3befoajPSAQKMo/nceDjfE6//T8KEk6bWqLRMdg82rdvz969e9m+fTtNmjQB4Mknn+SHH34wujiJpLZRz6cp9/2wkoLPd3OuaT8QkHIgj38GtOPM6+EUJCeYWqLRMNg8/vzzT2JjYwGIj4/nwQcf5NFHHyU5Odno4iSS2sqw4d1pErGWs+Mm4+xvgdCUcnnFN/zTtzlnXp9EQVK8qSXeNQabx4svvqhf/PiVV17Bzs4OV1dXnn32WaOLq05kV62kuglt5YD3o7PwfdqLkGfrQes2CI2Gyyu+5Z9+LTj92jMUJMWZWmaVURl6QWpqKl5eXmg0GrZs2UJiYiLW1tY0bty4OvRVG1OnTmXq1Kn6tQskkupgWJA7Z9JG0NLia6w9fMi2+QiP3xdwbfdmUn7+jtRVP9DokSfxe+4N6vk0NbVcgzA48qhXrx5paWns2LGDli1bYm9vj0KhoKSkpDr0SSS1npZdnwGgq/VellyxR/veWjqs3oNLr8EIjYaUXxbzT78AYl56irzY2pNYNdg8xo0bR6dOnRg7dizjx48HICoqiqZNa5drSiQ1hkMLhGtnLBRaelhtY+afKWQ16UjwTxvpsHovLn2GIEpLSV39AwcGBnJ80iNkHzH/9VENNo8PPviA7777jpUrVzJ58mQALC0ta90IU4mkJlF4Dgcg1HE3uUVaXl9/mZyiUhw7dCP4hz/puO4fGoQ+DEDGpnUcfqQ7R0b15+rOvzFwpdAao0qLAQ0aNIjmzZvzzz86d+zYsSP9+/c3qjCJ5J7CcygAAYoomjrkczm7hI+2pOuNwSG4M20jV9Nlyyk8RoahUKnI+mcHx8bfx6GhHUj7fSXi+gbv5oLB5pGSkkKfPn3w9fVl4MCBAPzyyy8888wzRhcnkdwz2DcD+xYohIZ3OsWiUsKu2Fx+O55drlh9/1a0+uh7uu26gPfTL6CsZ0tuzFFOTRtNVGhbLPf+TWlhgYkeojwGm8ezzz5Lly5dyM3N1S96PGDAALZulZsBSyS3xb0XAF4lUYT3cAMgYncGsVdunvti09ib5m8toMe+BJq8MBOVkwuFiRex+WURB3v7c3HB2xSlm3bJQ4PNY+/evcydOxcbGxsU1zf9dXV15erVq0YXZwiZmZl06NABOzs7k+qQSG5Jg+66PzP2M6K9E92b/Lukoaa04ryGpbMrTV6YSY99CTR9Yz5a5wZormUQ/9kc9vX0I+blCeSePl6DD/EvVdp6ITMzs9yx5ORkGjZsaDRRVcHe3p7NmzfTtWtXk+qQSG6J2/XvZtZxFJo8XhnojoONktgrRaw8cu22l1rY1sczbBp5b31Fy0+X4dC+K6K4mNRVSzh4XzDRYweRsW0DQqutgQfRYbB5jBkzhrFjx3LixAmEEJw7d47JkycTFhZWDfIqj6WlJS4uLibVIJHcFlsv3UuUQtYxnG1VTO3dANDtC5N8rRILJ1tY0OD+EXRcu48Oa/bhPvQxUCq5tncrxycM58DAQC4tXURpQX41P0wVRpi+/fbb/O9//6N79+7k5eXRoUMHJk+ezIwZM4wmavr06axfv56EhASio6MJDg7Wnzt//jzjx48nIyMDR0dHlixZQmBgYKXrLioqoqjo3zamWq0GdLtxmftAtzJ95q4Tao/WmtZpYd8KZX4ypZkn0Tp1pl8zG/72suFIciGf7khn7rBbR/D/1WrbtgMBC3/C95U5XP7pS1JXfk/+xbOcfXMKF+a/SaPHJuAx+hlsvPwM0ljZz8Kg7SY1Gg2RkZE8/fTT2NjY6LdoLMt9GItdu3bRtGlTevbsybp168qZR//+/Rk3bhxhYWGsWrWKDz74oNxqYAMHDmTLli23rHvWrFnMnj37puPLly/H1tbWqM8hkfyX1pqVNC/dQJyyP8ctwwDIKrHkh0v+aFHwaMMEfOrlVa3ywgIs/9mC1a4NKK/qFucSCgWlgR0p7nkfpQFBoLxzYyM/P58xY8bccbtJg/eqdXJyIisry5BLqoyfn18580hPT8ff35/MzExUKhVCCDw8PNizZw/+/v7Anc2josjD29ublJSUWrFX7ebNmxk0aJDZb+9ZW7TWtE5F4i+oDk9B69qV0t5/6I9H7L7KuhM5+LtZETHSA2UF/yFXVqsoLeXqtj9IWbqIrH3b9Mfr+fnjMWYS7o+Mw9LR+ZbXq9Vq3Nzc7mgeBjdbQkND2bJli36MR02SlJSEh4cHKpVOtkKhwMfHh8TERPz9/Rk4cCDR0dEMHDiQhQsX0qZNm5vqsLa2xtramoiICCIiIigtLQV0ORNz/pLfiNRqfGpMp4vuO6nMOYvyhvuFdXVj09k8YjOK2RtfRP8W9res4o5aLS3xuH8EHvePIO/CWS799CUpq3+gID6Wi3NnEL9gJo0eHovnk1Oxbx1UYf2VwWDzsLe356GHHmLw4MH4+PigvCEMWrBggaHVGZXbRRz/Rc6qlZgEu+tzwIqvQXE2WOm+e062Kh5r78SSA5ksj7pGv+Z2RkkH1G8WQItZn9L0lfdIW7eM5J++JO/MCS6v+JbLK77FsWMPvMZNxX3oSBQWFgbVbXBvS2lpKSNHjsTR0ZHs7GyuXbumf1U3Zc0LzfVhukIIEhMT8fHxMbguuZ6HxCRY2oOVk+7vBZfKnXo4yAkblYLYK0UcTjRub4mqvh2eY8Pp/NdRQn7ZifvQx1CoVGRH7eXU9DGcnPq4wXNoDI48FixYgLPzrdtL1Ym7uzshISEsXbqUsLAwVq9ejZeXlz7fYQgy8pCYDFtv3b4u+cng2Fp/2LGeBUPbOLL6aBYrDl+jo6/xFxVXKBQ4de6FU+deFKWncGnpIuI/e5crG9eQd+YEdq3aVbquSkceR44cwdvbGzc3N5o1a8apU6eqJL4yhIeH4+XlRXJyMqGhoeXMITIyksjISFq0aMG8efNYvHhxle4hIw+JybD11v2Zl3TTqRHtnVAAh5MKuJRVvRtmW7t70PTF2Vh76iJ3Q+fMVNo8Xn75ZUaOHMnx48cZPnw4r7zyimFKDSAyMpLk5GQ0Gg1paWn6NVMBAgIC2L9/P+fOnSMqKoq2bdtW6R5y0yeJybD10v2Zf7N5eDhY0tFXN2Tgrxh1jchRKA3LdZRR6WbL8ePH2bhxI1ZWVsydO7dKTQWJRALYXB8IVnSlwtNDAx04lJDPxhg1YV1dUSmNO47KWFQ68iguLsbKygoAW1tbCgsLq01UTSCbLRKTYX19GkVRxZ0M3ZvUx8FGSUZeKScum8f0+4qodORRXFzMZ599pn9fWFhY7j3ohpXXFmTCVGIyrK6bR3FmxadVSno0teOvGDW7Y3Np72WeI58rbR5du3Zl7dq1+vddunQp916hUNQq85BITIb17c0DoFez+vwVo2bXhTye6yMqHHFqaiptHjt27KhGGTXPf0eYSiQ1htX1oQ5FtzaPjj621LNUkJGrIfZKES3cbWpIXOWp0hqm9wKyt0ViMvTNlmtwi4FZViolwdebK0eSzDPvUWfNQyIxGWWRh9CAJveWxUK86wEQnVT9a3NUBWkeEklNY1EPuJ7DKL11VBFyPfI4drnglssUmpI6ax6yq1ZiMhQKUF3vQdHceu2OJm5WONgoKSwRFS6SbGrqrHnInIfEpFjomiS3izyUCgUBDXWJ0rPp5jeuqs6ah0RiUvSRx+3zGS0aWANwLl1GHhKJBMDiunmU3sE83KV5SCSSGylrttwp8rjebIm7WkSJmSVN66x5yISpxKSoKhd5NLJXUc9SgUYLKWrzWom+zpqHTJhKTIpF5XIeCoUCTyfdhNRL2bV8o2uJRGIElDpDQHvnBX+8nHQLEl/KkpGHRCJRXp9Wpr3z3Koy80iW5iGRSFBcNw9x56aIbLZIJJJ/0Uced44mGjvqIo/0XGkeEolEcX1jpUpEHq71dWuMZuaV3moSrkmos+Yhu2olJqVs0WFx55yHa31dlFKoERQL8/nJmo+SGkZ21UpMSlnOQ3vnyKOepZL6Vrqfap7G4K2Wqo06ax4SiUlRljVbKteDUtZ0ySuV5iGR1G0MiDzg36ZLrjQPiaSOo6h8zgPAqZ6ufKFWmodEUrdRVr63BcDOWvdTLdSaz0/WfJRIJHUJReXHeQDYWesijyJt1baGrA6keUgkpkBZ+RGm8G/kUVQqzaNamDZtGr169WLu3LmmliKR3J6ynEcl5rYA2NuURR7m85M1HyV3SVRUFCqVit27d3PkyBHS0tJMLUkiuQ2G7QBnr895yMjD6Bw4cID+/fsD0KdPHw4fPmxiRRKJ8dA3W8zIPMyn3+cGpk+fzvr160lISCA6Oprg4GD9ufPnzzN+/HgyMjJwdHRkyZIlBAYGkpWVRZs2bQCwt7cnKyurwrqLioooKvp3PUi1Wg1ASUkJJSXmNeX5v5TpM3edUHu0mkqnUqvFAtBqSymtxL2tlFoAioXS+FqvT5jRaDQG/Q7M0jxGjBjBjBkz6Nmz503nwsPDmTRpEmFhYaxatYqwsDAOHTqEk5OT3ghycnJo1qxZhXW///77zJ49+6bj27dvx9bWPHcj/y+bN282tYRKU1u01rROf81ZAoHk5GSi0/68Y/n0ImugGRqhMLrW+gUFKIF9+/ahvZxBfn7ldqhTCGFO8/TK4+fnx7p16/SRR3p6Ov7+/mRmZqJSqRBC4OHhwZ49e7h27RorVqxgwYIFjBw5ki+++IKGDRveVGdFkYe3tzcpKSm4urrW1KNViZKSEjZv3sygQYOwtLQ0tZzbUlu0mkqn8txnWJx6B63PaEo7fH7H8slZJTy1/BLWilLWTvQxqtZD/QIoTI4n6NfdOAR3Rq1W4+bmRnZ2Ng4ODre8ziwjj1uRlJSEh4cHKpVOtkKhwMfHh8TERPr378+SJUvo1asXoaGhFRoHgLW1NdbW1kRERBAREUFpqS7bbWlpadZf8huRWo1Pjeu00OUulEoFykrct75uEXU0QmF8rQpd8lalUhlUd60yjzsRERFR6bJTp05l6tSpqNVqHB0dq1GVRHL3WKl0CdNSlJRqBeZgx7Wqt6WseaHR6AbWCCFITEzEx8fH4Lrkeh4S02JYV62N6t/yxWayf0utMg93d3dCQkJYunQpAKtXr8bLywt/f3+D65LreUjMgkqmHK1uMI8ijTSPWxIeHo6XlxfJycmEhoaWM4fIyEgiIyNp0aIF8+bNY/HixSZUKpHUDEqFAsvrQzzMxTzMMucRGRl5y3MBAQHs37//ru/x34SpRFKzGNZsAbBRKSkp1VKk0VaDHsMxy8ijJpDNFklto6zpUmwmkUedNQ+ZMJWYB5U3gusdLpSaR+BRd81DRh6S2obF9fEYpWYyrrPOmodEYlIUhuc8lDLykEgk/1L5KEJGHmaCzHlIahsW13+tWhl5mBaZ85CYFsObLRZKGXlIJJIyDDACi+t+I3MeJkY2WyS1DWVZ5KGVkYdJkc0WiWmpSrNF96eZzIuru+YhkZgHVehtkZGHRCIxBKXMeUgkkqpQ1ttiJoGHNA+JxCToR5ga0GzRjzA1D/eos+Yhe1sktY1/R5iaWMh16qx5yN4WSW1DRh4SiYSqdNUq9b0txtZSNaR5SCSmxIARpkq938jIQyKR1GKkeUgktQzziDukeUgkJsLwrlrDsyTVS501D9lVK6mtmMmM/LprHrKrVlLrMLPQo86ah0RiUqowwrQMMwk8pHlIJLUFMws8pHlIJJKqIc1DIjElVcl+mkm7RZqHRGISDG+EVGGrl2rlnjGPzMxMOnTogJ2dnamlSCTVipkEHveOedjb27N582a6du1qaikSiQGYixUYzj1jHpaWlri4uJhahkRSSareBjEXuzGpeUyfPh0/Pz8UCgVHjx4td+78+fN0796dFi1a0KlTJ06dOmUakRKJmaAws6SHypQ3HzFiBDNmzKBnz543nQsPD2fSpEmEhYWxatUqwsLCOHToEDExMUyZMqVc2SFDhvDaa69V6p5FRUUUFRXp32dnZwO6nIm5U1JSQn5+PlevXsXS0tLUcm5LbdFqKp1KdS4W+QJtThGlV69W6pr83Gw0hQXkZiu4etV4RpJXUkqhVnAtKxvN1avk5OQAIO7UEyTMAF9fXxEdHa1/n5aWJuzt7UVJSYkQQgitVisaNmwozp8/f8e6BgwYcNvzM2fOFOgiP/mSL/m6zSspKem2vyWTRh63IikpCQ8PD1QqnTyFQoGPjw+JiYn4+/vf8rqBAwcSHR3NwIEDWbhwIW3atLmpzOuvv86LL76of5+VlYWvry+JiYk4Ojoa9Tk6depk1LkzarUab29vkpKScHBwMFq9xtYJtUdrdemE2qP1vzqFEOTk5NC4cePbXmeW5lFVtmzZcscy1tbWWFtb33Tc0dHR6F8eCwsLo9cJ4ODgYNR6q0sn1B6txtYJtUdrRTor8x+pWfa2eHt7k5KSgkajAXROmJiYiI+Pj4mVGcbUqVNNLaFS1BadILVWB1XVaZbm4e7uTkhICEuXLgVg9erVeHl53bbJYo7c618eUyC1Gp9aaR7h4eF4eXmRnJxMaGhoOXOIjIwkMjKSFi1aMG/ePBYvXlwtGqytrZk5c2aFTRlzQ2o1PrVFJ5ifVoUQ5rIukUQiqU2YZbNFIpGYP9I8JBJJlZDmIZFIqkSdMY/KzpX57rvvaN68Oc2aNeOZZ56hpKSkhpVWTmt8fDx9+/bF0dGR4ODgGtdYRmW0btu2jc6dO9O6dWsCAwOZMWMGWm3N7plYGZ379+8nODiY4OBgAgMDCQ8PLzeVwZy0liGEoH///jg5OdWcwBtuXifo16+fWLx4sRBCiF9//VV07NjxpjIXL14UHh4eIiUlRWi1WjF8+HDxxRdf1LDSymm9evWq2L17t/jjjz9EUFBQzQq8gcpoPXLkiLhw4YIQQoiCggLRo0cP/TU1RWV05uXlieLiYiGEEKWlpeKhhx4SCxYsqEmZQojKaS3j448/FhMnThSOjo41I+4G6oR5VHauzIcffijCw8P17zds2CB69OhhllrL2L59u8nMo6pzkKZOnSpmzpxZAwp1VEVnQUGBCA0NFZ988kkNqdRhiNaTJ0+KXr16idjYWJOYR51ottxursyNJCYm4uvrq3/v5+d3Uxlz0WoOVEVramoqq1atYtiwYTUl0yCd8fHxBAUF4ebmhqOj400zuM1Fa0lJCc888wyRkZFYWFjUqMYy6oR5SMwDtVrN8OHDmTFjBh07djS1nArx8/Pj2LFjpKamUlRUxJo1a0wtqUJmz57NI488QqtWrUymoU6YR2Xnyvj4+JCQkKB/Hx8fX+PzaWrTvB5DtObk5DBkyBAefPDBcrOazU1nGXZ2dowaNYply5bVlEyg8lp37tzJ559/jp+fHz179kStVuPn58eVK1dqTGudMI/KzpV59NFHWb9+PampqQghWLRoEaNGjTJLreZAZbXm5uYyZMgQhgwZwptvvmm2OmNjY/W9a8XFxaxdu5Z27dqZpdbdu3eTkJBAfHw8e/bswcHBgfj4eBo0aFBzYms8y2Iizpw5I7p27SqaN28uOnToII4fPy6EEOLpp58Wv/32m77c119/LZo2bSqaNm0qJkyYoM++m5vWvLw84enpKdzc3ISlpaXw9PQUr732mllqnTNnjlCpVCIoKEj/mjNnjtnpjIyMFIGBgaJdu3aidevWYtq0aaKgoKBGdVZW643ExcWZJGEq57ZIJJIqUSeaLRKJxPhI85BIJFVCmodEIqkS0jwkEkmVkOYhkUiqhDQPiURSJaR5SCSSKiHNoxaxe/duvLy89O/79u3LwoULTSdIUo6vv/4aLy8vnJ2dmTZtmn6I+a2YN28eM2bMuOX59957jzfeeMPYMo1HjQ9Lq+OcOXNGDBs2TLi6ugp7e3sREBAg5s2bV6W6+vTpU+1Txrdv326S0Yt3A1Bu+9K7ZebMmeLBBx+8bZl9+/YJJycncfLkSZGZmSlCQkJu+++alZUlGjRoIK5cuXLLMjk5OcLd3V2kpKRUVXq1IiOPGmbo0KEEBQWRmJjItWvXWL16NU2bNjW1rGrFFKux1TSLFi1i3LhxBAYG4uzszNtvv82iRYtuuVn0Tz/9RO/evXFzc7tlnXZ2dtx3331899131SX77jC1e9Ulrly5IgCRmJh4yzKpqali5MiRws3NTXh7e4v//e9/+oVh/hsF3C7yGDJkiPjqq6+EELr/5SwsLMSrr74qhNAtMOPm5iaioqKEEEK88sorwsfHR9jZ2YlWrVqJX375RQghREZGhrCxsRGAqF+/vqhfv77YtWuXEEKIzZs3i06dOglHR0fRunXrcnMuxo8fLyZMmCBGjhwp7O3txWeffVahxuXLl4t27doJe3t74ePjo189S6vVivnz54umTZsKZ2dnERoaql+JTAjdxugffPCB6NKli7CzsxO9e/fWf6adOnUSgKhXr56oX7++eO+994QQQsTGxophw4YJNzc34ePjI959911RWloqhBBi8eLFIigoSLzzzjuiQYMGwt3dXf+5rl27VlhaWgoLCwv9Z1ARzZs3F7/++qv+/eXLlwUgLl++XGH5QYMGiYiIiArP3ciPP/4oOnfufMdypkCaRw2i1WpFQECAGDBggFi5cqWIj4+/qUz//v3FmDFjRE5OjoiPjxetW7fW/wAMMY8PPvhAPPbYY0IIIdatWyeaNWum/xIePXpUODs76388S5cuFWlpaUKj0YgVK1YIa2trcfHixQrvKYQQx44dE05OTmLr1q2itLRU7N69Wzg4OIgzZ84IIXTmUa9ePbFx40ZRWloq8vLybtK3fv164eLioq8jLS1NHDlyRAghxA8//CAaN24sjh8/LgoKCsSLL74oWrdurTdRX19f0bZtW3Hx4kVRUFAg7rvvPjF+/Hh93fyn2ZKXlyd8fX3FJ598IoqKikRCQoIIDAwU3377rRBCZx4qlUrMnz9fFBcXi+3btwuVSiViY2OFEJVrttjZ2YkDBw6UO2Ztba036P/SoEEDsXXr1tvWKYQQUVFRwsrK6o7lTIE0jxomJSVF/2NQKpWiVatWYtOmTUIIIZKTkwUgUlNT9eWXLVsmmjdvLoQwzDwOHjwo3N3dhRBCTJ8+XXz22WfC1dVVZGdniwULFoiHHnrolhqDgoLE0qVLK7ynEEJMmTJFvPDCC+WOjRkzRrzzzjtCCJ153OnHNmTIEDF79uwKzw0cOLBcvqCwsFDY29uLvXv3CiF05lEWVQmhM782bdro3//XPH755RcRHBxc7h5ff/216N+/vxBCZx6NGjUqd97f31+sWrVKCFF587C3txeOjo76l0KhEIcPH66wvEqlEocOHbptnUIIce7cOQFUaMCmRuY8aphGjRrx8ccfc+rUKa5cucJ9993Hww8/TGZmJsnJydjY2NCwYUN9+aZNm5KcnGzwfUJCQigsLOTUqVNs27aNgQMH0rNnT3bv3s22bdvo37+/vuwnn3xCYGAgjo6OODk5cfLkSTIyMm5Zd3x8PIsWLcLJyUn/+u2337h8+bK+zJ0WL0pISKB58+YVnktOTsbPz0//3tramsaNG5f7HBo1aqT/e/369cnJybmt3pMnT5bT+9JLL5Gamqovc+NnXpk6/4uHhwcbN24kKyuLrKwsrl27hkqlwsPDo8Lyzs7OqNVq/fvJkydjZ2eHnZ0dkydP1h9Xq9VYWVlha2tbaS01hTQPE+Li4sKsWbPIy8sjLi4OLy8vCgsLSUtL05eJj48v1z1bWSwsLOjduzcrV64kMzOTVq1a0b9/fzZv3syuXbvo168fAHv27GHWrFn8+OOPXLt2jaysLNq0aaNP9CmVN39FvL29ef755/U/lKysLHJzc/nqq6/0ZSq67kZ8fX2JjY2t8JyXlxfx8fH698XFxVy+fLnSn4NCobhJb4cOHcrpVavVt93S4Ebu9CwA3bp149y5c/r3Fy9exNPTs5zJ3UhwcDBnzpzRv1+0aBG5ubnk5uayaNEi/fGYmBiTbq1xO6R51CDXrl3jzTff5MyZM5SWlpKfn8+CBQtwcXGhZcuWeHp60q9fP15++WXy8vJITEzkvffeY/z48VW6X79+/fj000/p06cPAP3792fx4sXY2NjQpk0bQPc/m4WFBQ0aNECr1fL9999z8uRJfR0NGzYkJyeH9PR0/bHw8HAWL17M9u3bKS0tpaioiP3793P69OlKawsPD+fTTz9l586daLVa0tPTiY6OBuCJJ57giy++ICYmhqKiIt588008PT3p3Llzpepu2LAhFy5c0L8fNmwYaWlpfPnllxQWFlJaWsrZs2fZsWNHpetLSEi47biNyZMn8+GHHxIXF4darWbGjBlMnjz5JiMrY/jw4Wzfvv2O9962bVuNLhZtCNI8ahArKysuXbrE/fffj6OjIz4+Puzdu5e//vqL+vXrA7B8+XIKCgrw9fWlR48eDB069LYDiW5Hv379UKvV+iZKmzZtqFevnj7qABgyZAgjRoygbdu2NG7cmFOnTtGjRw/9+YCAAJ5++mlat26Nk5MTe/bsoX379qxYsYI333yTBg0a4OnpyVtvvWXQBkkPPfQQCxYsYOrUqTg6OtKpUydOnDgBwLhx45g2bRrDhg2jUaNGHDt2jN9//12/ovidePfdd5k+fTrOzs7MmzcPOzs7tmzZwtatW/Hz88PV1ZUxY8aUa7bcjpEjR+Lg4ECDBg1uublSt27deOGFF+jRowc+Pj40btyYl1566ZZ1Pvnkk+zcuZOrV6/eskxeXh5//vknEydOrJTOmkauJCaRmIj333+frKwsPvjggwrPz507l7y8PN57770aVlY5pHlIJJIqIZstEomkSkjzkEgkVUKah0QiqRLSPCQSSZWQ5iGRSKqENA+JRFIlpHlIJJIqIc1DIpFUCWkeEomkSkjzkEgkVeL/AQaEjlrdlijtAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(2.6, 3.35), layout=\"constrained\")\n", "pe.plot.swrc(wosten, label=\"Wösten\", ax=ax)\n", "pe.plot.swrc(rosetta, label=\"Rosetta\", ax=ax)\n", "pe.plot.swrc(saxton, label=\"Saxton\", ax=ax)\n", "ax.set_ylim(1e-2, 1e6)\n", "ax.set_xlim(0.0, 0.46)\n", "ax.xaxis.set_minor_locator(mplticker.MultipleLocator(0.05))\n", "ax.xaxis.set_major_locator(mplticker.MultipleLocator(0.1))\n", "ax.legend(ncol=1, loc=\"upper right\", frameon=True)\n", "ax.set_ylim(1e-1, 1e6)\n", "# ax.set_title(\"SWRC from pedotransfer functions\")\n", "ax.set_ylabel(\"Pressure head |\\N{GREEK SMALL LETTER PSI}| (cm)\")\n", "ax.set_xlabel(\"Soil water content \\N{GREEK SMALL LETTER THETA} (-)\")\n", "fig.savefig(fig_path / \"ptf_swrc.png\", dpi=300, bbox_inches=\"tight\")" ] }, { "cell_type": "code", "execution_count": 8, "id": "95f42ef8", "metadata": {}, "outputs": [], "source": [ "# Estimate parameters using HYPAGS\n", "ks = 1e-3 # saturated hydraulic conductivity (m/s)\n", "hypags: pe.Genuchten = pe.SoilSample(k=ks).hypags()" ] }, { "cell_type": "code", "execution_count": 9, "id": "c1ee315c", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Genuchten(k_s=0.001, theta_r=0.105, theta_s=np.float64(0.30896365527138614), alpha=np.float64(9.439965924107144), n=np.float64(1.3453214358373247), l=0.5)" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# pe.plot.curves(hypags, label=\"HYPAGS\")\n", "hypags" ] }, { "cell_type": "markdown", "id": "732dd9ed", "metadata": {}, "source": [ "## Fitting a soil model" ] }, { "cell_type": "code", "execution_count": 10, "id": "6f010ea6", "metadata": {}, "outputs": [], "source": [ "# Fitting a Brooks-Corey soil model to existing Mualem-van Genuchten soil model\n", "bcf = pe.SoilSample(h=h, theta=theta, k=k).fit(pe.Brooks)" ] }, { "cell_type": "code", "execution_count": 11, "id": "11e43ca7", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAggAAAFRCAYAAAARlVXBAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjksIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvJkbTWQAAAAlwSFlzAAAPYQAAD2EBqD+naQAAeVRJREFUeJzt3Xd8U9X/x/FXVvduKRS6KGWWvXfZU5zgAJkqIFXwK4oDfwJORAFXlYoKIkMUZKiogOwhG2XIhg5WoaWbtml6f39EIqUtNCVp0vbzfDzy/Zp17zsJJ/3knHPPVSmKoiCEEEIIcRO1rQMIIYQQwv5IgSCEEEKIQqRAEEIIIUQhUiAIIYQQohApEIQQQghRiBQIQgghhChECgQhhBBCFCIFghBCCCEKkQJBCCGEEIVIgSCEEEKIQqRAEEIIIUQhUiAIIYQQohApEIQQQghRiBQI5UBoaCj+/v7o9XrTbRs3bkSlUvHcc8/d1bbnz5/P/ffff3cB7cC6devo3LkzYWFhtGzZktatW/PFF1/YJEuXLl1YuXJlkfetXLmSP//8s2wDlSOhoaHUrVuXpk2bUr9+fQYPHkxmZqZF9zF16tS7bjeJiYmMHDmSsLAwmjVrRvPmzXnnnXcsE7AUnnjiCRo0aMADDzzA6tWr+d///gfAuXPnmDNnjs1yifJNCoRyIjg4mNWrV5uuf/XVV7Rs2dKGiezH2rVrGT58ONOnT+fMmTPs3buXH3/8kfPnz9s6WiFSINzZ0qVLOXjwIEeOHCE1NZX58+cX+TiDwVC2wf51/fp1IiMjCQkJ4eTJkxw4cIBt27bh6upq1nby8vIskufy5ct89913HDp0iBUrVnDvvfcye/ZsQAoEcXekQCgnRo4cyddffw1Aamoqf/75J3369DHdf2tPwM8//0yXLl0AuHTpEl27dqVFixZERETwzDPPkJ+fX+R+vv32W9q0aUPz5s3p3Lkzf/31l2n7PXr04LHHHqNBgwa0b9+eo0eP8sADD1C/fn169epFRkZGoe1t376dRo0aFbitS5curFq1iry8PHr37k3Lli2JiIgo8Gtx06ZNNGzYkHHjxtGkSRMiIiLYu3dvkZnfeOMNXn/9ddq3b2+6LTAwkGnTppmu79mzh27dutGyZUuaNWvGDz/8ABi/QL28vJgyZQotWrQgPDycNWvWmJ6nUqlISUkxXffz8+PcuXMA/PPPP/Tu3ZvGjRvTuHHjAl/E27Zto1OnTtSqVYuxY8cCsGbNGlavXs37779P06ZN+fLLL0v8njdq1IiWLVty5syZIt+Du6YokJdp3YuimBUpNzeXrKwsvL29AeP70bVrVx566CEaNWrE7t27+f3332nevDmNGzcmMjKSo0ePmp7//vvvExERQaNGjRgyZAipqamF9nH06FEaNmzIr7/+yvXr13nkkUdo0KABTZo0oVevXkXmWrx4Me7u7kydOhWNRgOAi4sLEyZMACAjI4NRo0bRsGFDGjZsWODfYZcuXRg/fjzt2rUzbf+DDz6gdevWNG/enD59+hAbG0t2djbVqlUjPj7e9NxXX32Vl156qUCWlJQUunbtSnZ2Ni1atGD69OkFvgvGjh3L8ePHadq0Kffee69Z778QWlsHsHejl8SRnGW9Xyo+Lhq+eCz4jo/r0KEDn332GRcuXGD16tUMGjTI9OV0J15eXvz000+4ublhMBi47777+P7773n00UcLPG779u0sWbKELVu24OjoyNatWxk8eDBHjhwBjH9kDx06RHBwMEOHDmXAgAHs2LGDqlWrcs899/DNN98QFRVVKHdOTg579+41/YE7fvw4/fv3R6PRsHjxYnx9fVEUhXHjxvHJJ5/w8ssvA3Ds2DG++uorPvvsM+bMmcPkyZP5/fffC72+/fv388knnxT7+lNSUhg9ejRr1qwhICCAq1ev0rx5c1NBkZqaSuPGjZk2bRq//fYbEyZMoF+/frd9T/Py8rjvvvuYNm0ajz32GABXr1413X/69Gk2btyIXq+nQYMG7Ny5k379+nHvvffStGlTUxd3Sd7zgwcPUrNmTV5++WXee+89YmJibputVAxZ8GN1y2/3Zg9eAO2df2U/8sgjODs7c+7cOVq0aMHDDz9sum/Xrl0cOHCAunXrkpiYSP369dm0aRONGjVi0aJFDBw4kCNHjvDbb7/x9ddfs3PnTry8vBg9ejQvv/wyn3/+uWlbmzZtYty4cSxcuJDmzZuzYsUKUlJSTEVGcnJykfn27dtHu3btis3/5ptvkpOTw99//83169fp2LEj9erV45FHHgHgxIkTbNmyBZ1Ox+LFizl+/Dg7d+5Eo9Hw7bffMm7cOH755ReeeOIJPv/8c9555x1ycnKYN29eod4nLy8v1qxZQ9OmTTl48CBAgR6XOXPm8Nxzz5nuE8IcUiDcQXKWgasZlukKvFtDhw5l/vz5rFy5kkWLFrFo0aISPS8/P5+XXnqJbdu2oSgKiYmJNGzYsFCBsGrVKv766y/atGljui05OZnr168D0K5dO4KDjcVMy5Yt0ev1VK1aFYBWrVpx8uTJIvc/cuRI5s2bR8uWLfnmm28YMmQIWq2W/Px8Zs+ezS+//EJeXh6pqakFegHCw8NNWdq1a8cHH3xQotc7ZMgQjhw5wqVLlzhx4gQ7duzgzJkz9O3bt8Djjh8/TlhYGE5OTjz44IOm/Zw+ffqO+zh+/DjZ2dmm4gCMvQs3PPLII2i1WrRaLU2bNuX06dNF/lEpyXtes2ZN03/frhCqKJYuXUrTpk3Jy8tjzJgxvPTSS8ycOROA9u3bU7duXcBYLDRq1MjUQzVkyBCioqI4f/4869ev55FHHsHLywuAp59+mkGDBpn2sWHDBn777TfWrl1r+jfdpEkT/vnnH8aNG0dkZOQdi8TirF+/npkzZ6JWq3F1dWXYsGGsW7fOVCA8/vjj6HQ6wDjktGfPHlq0aAEUHDYZN24crVu3ZsqUKfzwww+0bt2akJCQUmUSojSkQLgDH5eS/Uovi+0PGzaM5s2bU6dOHWrXrl3gPq1WW+DLJTs72/Tfs2bNIjExkV27duHk5MTzzz9f4P4bFEVh+PDhxU62cnJyMv23RqMpdP3GmGr79u3JysrC0dGRXbt2MXz4cJo0acIHH3zAggUL+PnnnwFjV+2GDRvYvHkzHh4efPzxx2zYsKHY/RU3ZtusWTN2795Ns2bNAEyFk0qlIj8/H0VRiIiIYMeOHYWee+7cORwdHVGpVKb93Pw+3nq9qPetJO9VcdnNfc8tNW5diMbF+AvfmjQuZj1cq9Xy0EMP8eKLL5oKBDc3t1Lt+sbne0N4eDjHjh3jzz//NBUIYWFhHD16lA0bNrB+/XomTZrEwYMHmTJlClu2bAGMw0EtWrQwawLsrfu++TUoisIrr7zC6NGjCz2vRo0adO7cmaVLl/L555/zxhtvlHifQliCFAh3UJLu/7JSvXp13n33XerVq1fovvDwcFOX5o2uyxuuXbtGtWrVcHJy4tKlS/zwww889NBDhbZx7733MmTIEMaOHUtwcDD5+fns37/f7MmQt/4hrl69Oq1ateJ///sf/v7+REREmHL5+fnh4eFBeno68+fPN31Zm+P//u//eOKJJ2jSpAlt27YFKDDzvX379pw9e5b169fTo0cPAA4ePEiDBg3uuO3w8HB27dpFv379+PHHH03brVu3Li4uLixZsqTAEMPNvQhF8fDwKDAWbqn3/K6pVCXq/i9rGzZsMPUY3Kpt27YcOnSIw4cP07BhQ7777jtq1KhBjRo16NGjBxMnTuT555/Hw8ODmJiYAnMKgoODiY6Opnfv3mRmZjJy5EgSEhLw9vbm3nvvpU+fPqxcuZL4+Hg+/vjjAvutVasWM2bM4M033+TVV19Fo9Fw/fp15s6dy/jx4+nRowdfffUVkZGRZGVl8e233xaaO3DD/fffz8yZMxk4cCA+Pj7o9XoOHz5sKnYnTJjAoEGDcHNzM/3bNcet/96EMIcUCOXMyJEji7y9bdu29OvXj4YNGxIQEECHDh3YtWsXYPySGThwIBEREVSvXr3YL5pOnToxY8YMHnjgAfLy8sjNzaV///4W+WM1cuRIHn744QJjwMOGDWPVqlXUrVuXKlWq0KlTJ2JjY83edp8+ffjqq6948cUXuXDhAlWqVMHBwYFPPvkEd3d3NBoNv/zyCy+88AITJ05Er9cTHBxc7KGIN5s9ezbjx4/ntddeo3///vj6+gLGX7erVq3i2Wef5Z133kGtVjNu3DjGjBlz2+0NHTqUESNGsHLlSqKionjyySet9p6XVzfmIOTl5RESElLsLPwqVaqwaNEihg0bRl5eHt7e3vzwww+oVCr69u3L4cOHadeuHWq1msaNG/PZZ58VeH5AQAAbNmygT58+pKenU7t2bV555RUURSEvL4+hQ4fSuHHjQvt1cXFh8+bNvPzyy4SHh+Pm5oZKpWLw4MGAsWAdP368aehj0KBBBeZR3GzIkCEkJSXRtWtXwDi3ZdSoUaYCoW3btnh6ejJmzJhCPREl0bhxYyIiImjYsCFhYWEFjoQS4k5UimLm1GIhhBBl4vz587Rs2ZITJ07g7u5u6ziikpHDHIUQwg69/vrrtGnThunTp0txIGyiQvUgxMTE8P3332MwGFi3bp1pprAQQgghzFNh5iDExcVx6NAh/vjjD1tHEUIIIcq9CjPEsHbtWjIzM+nWrRtTp061dRwhhBCiXLPLAmH8+PGEhoaiUqkKrQB28uRJ2rdvT506dWjVqpVpxbnExERUKhUbNmzg7NmzsnKYEEIIcRfscohh4MCBTJo0iY4dOxa6b8yYMYwePZoRI0awbNkyRowYwZ49e/D09CQyMhKAzp07m9Yfv1VOTg45OTmm6/n5+SQnJ+Pr61uqw4iEKI8URSE9PZ3q1aujVpfsd4K0HSHsV2nadEk2ardCQkKUAwcOmK5fvnxZcXd3V/R6vaIoipKfn69UrVpVOXnypLJv3z5l4sSJiqIoyoQJE5Q///yzyG1OmTJFAeQiF7mAEh8fX+L2KG1HLnKx/4s5bfpO7LIHoTjx8fEEBASg1Rpjq1QqgoODiYuLo1u3bixYsIAuXbpQp06dAmvb3+yVV17h+eefZ+7cucydO5e8vDxOnz7NiRMn8PHxsV72vxYTFjeZPTltOFJzDkNaeJn1fL1ez8aNG+natavdHp0hGS2jLDKmp6dTs2ZNsw6fs1XbMZe9fsaSyzySyzzJycnUqVPHoofElqsC4U4+/PDDOz7G0dERR0dHJk6cyMSJE0lLS8PT0xMfHx/TKnnWkObhioeLCje1Fjd3L7P3pdfrcXFxwdfX167+Ud5MMlpGWWS8sV1zhgZs1XbMZa+fseQyj+QqHUsO99nlJMXiBAUFcfHiRdMJaxRFIS4urlTr90dHR9OgQQNatWpl6ZjFkDFaUTGUfdsRQthCuSoQ/P39ad68OQsXLgRg+fLlBAYGEh4ebva2oqKiOHr0KHv27LF0zNtSoZTp/oSwNFu1HSFE2bLLAmHMmDEEBgaSkJBA7969CxQAMTExxMTEUKdOHaZPn868efNKtQ/pQRCidKQHQYjKwS7nIMTExBR7X926ddm5c+dd7yMqKoqoqCjTOKoQomSk7QhROdhlD0JFpNzoQVApMsgghBDC7lXaAqHMu0llhEFUEDLEIETlUGkLhLKfaKX693+l/0CUbzJJUYjKodIWCPIrSIjSkbYjROVQaQuEMv8VJGvViwpCehCEqBwqbYFgKzLEIIQQojyQAkEIIYQQhVTaAsFWCyWpMC4RLUR5JXMQhKgcKm2BIOOoQpSOtB0hKodKWyCUPZmkKIQQovyQAqHM3BhiyLdxDiGEEOLOKm2BYLs5CDL/QJRvMgdBiMqh0hYIZT2OqqiMb7UUCKK8kzkIQlQOlbZAKHvSgyCEEKL8kAKhrPzbg6BWyRwEIYQQ9k8KhLIiQwxCCCHKESkQysy/PQhyFIMQQohyQAqEMqKo5DBHIYQQ5UelLRDK/lCtG0MMQpRvcpijEJVDpS0Qyv5QrX97EGSSoijn5DBHISqHSlsglLUb6yDIHAQhhBDlgRQIZea/oxjkZI5CCCHsnRQIZUX9Xw+C1AdCCCHsnRQIZUR100qK0oMghBDC3lWYAuHcuXMEBATQpUsXhg0bZus4hd20UJLUB0IIIeyd1tYBLKl///58+eWXto5RjP+WWlakC0EIIYSdqzA9CAC///47nTp1YtGiRbaOUpj0IIiKKuOcrRMIIazALnsQxo8fz+rVq4mNjeXAgQM0bdrUdN/JkycZPnw4V69exdPTk/nz5xMREUFAQADHjx9HpVLRu3dv+vTpg6+vb6Ft5+TkkJOTY7qelpYGgF6vR6/XW+01GfKNZYGafAyGfLP3dePx1sx4tySjZZRFxtJsu7i2o2zsS3rLN3GqMxJUtl8KzF4/Y8llHsllHmvkUSl22N+9ZcsWwsLC6NixIytXrixQIHTr1o1hw4YxYsQIli1bxnvvvVdowZbJkydz33330bp160Lbnjp1KtOmTSt0++LFi3FxcbH4a7nBkB3Pg6rJXMirwQc5XxPpc9lq+xLiTrKyshg8eDCpqal4eHiU6DnFtZ3Uee54uKg4k9+cE04jyFF5WTitEOJOStOm78QuC4QbQkNDCxQIiYmJhIeHk5ycjFarRVEUAgIC2LZtG9WqVcPNzQ1FUejTpw/z588nICCg0DaL+hUUFBTExYsXi+xxsJSzJ3dS5/AALuUFsDRgE+M6mrcvvV7PunXr6NmzJzqdzkop745ktIyyyJiWloafn59ZXybFtZ0ZH7/MBP8vcVDlkoUn+qYf4FLzAavkLgl7/Ywll3kkl3mSkpIICAiwaIFgl0MMxYmPjycgIACt1hhbpVIRHBxMXFwcZ86c4dVXX0Wn0zFw4MAiiwMAR0dHHB0diY6OJjo6GoPBAIBOp7Pqh63VOfz7XwpqtabU+7J2TkuQjJZhzYyl2W5xbedyleGMvhrJq16vU0d3HA4+xeXYNVSN/BAcvCwb3Az2+hlLLvNIrpKxRpYKM0mxV69e7N27l507dzJx4sQ7Pr6s15NXmY5iUOQoBlGu3dp2Xu7hx9AeHXg5YwHfpo/CoKipem0FqavakBH3h43TCiFKq1wVCDeGAvLy8gBQFIW4uDiCg4PN3laZn5HOdBSDrKQoyrdb245KpaJ7XXe+fDycE1Vf5NmkL0nIC8JTuYTbnw+SsPF/kHfdxqmFEOYqVwWCv78/zZs3Z+HChQAsX76cwMBAwsPDzd5WmZ+R7t/Z3WryZSVFUa4V13Z8XLW80T+Ah7r34n/p37Eq8yEAAq98zZUV7Um/IGd/FKI8scsCYcyYMQQGBpKQkEDv3r0LFAAxMTHExMRQp04dpk+fzrx580q1j7LuQVBJgSAqiNu1nRu9CTGP12NvlbeYlPwRVw1+VFHO4Ly1D+e2vAH5eTZILYQwl11OUoyJiSn2vrp167Jz58673kdUVBRRUVGkpaXh6el519u7E5VKA/y7kqLV9yaE9ZSk7dzoTdhw4iEmbG7Ek47v0NV5PaGXZpLw4zo8IufiUaVeGScXQpjDLnsQykJZ9yAoauNRDBoMMklRlGslbTs3ehM+HtKEP3w+5q1rb5KR70Zg/t84bojk5I5PkO40IexXpS0Qyn4OgrEHQYN0r4ryzdy24+uq5c17qtM2cgTj075nX05rHFXZ1E54jVM/9ic1Od7KiYUQpVFpC4QypzIeo6pRGWSIQVQ6KpWKHvU8eH9IG1Z4zuOj1BfIURwJN2xHvbY9R3cvtHVEIcQtKm2BUOaTFNXGHgQtedKrKsq1u2k7vq5a3hxQg4hOE3gubTHHc+vjrk6jwbkoDv/4GKmpV62QWAhRGpW2QCjzIQaNsQdBrVJQ8vPLZp9CWMHdtp0bvQlvDe7CQo+lzE9/CoOioWHeGvJ+bcuh/T9bOLEQojQqbYFQ5lT/HTCilnkIQuDrquWNAcHU6PA6L6bNIy4vGF/1FRqdGsLuFU+Tkp5m64hCVGqVtkAo+yGG/9bJVin2dZpQIcxhybajUqnoWc+DyYP787XbSn7MfBiA1vrFZP3cgQN/bbnrfQghSqfSFghlPcSg3NSDoFKkB0GUX9ZoO76uWqYMCMO93UxeT4vmisGf6po4Gh+7ny0rXyElM9ti+xJClEylLRDKmkotQwxC3M6N3oQJjz3KHNfVrL/eG43KQOfcz7i6qgt7Dx2wdUQhKhUpEMqIWqXGoPx7wibpQRCiWL6uWl4bUA+lzZe8l/EuafkehGv/odHR3vy26j1Ss2SIToiyIAVCWVGBAeOhjirFYOMwQti3G70JTz4yms9dVrM7ux2Oqhz65LzDuRX92HX0uK0jClHhVdoCocwnKQIGxTjMIEMMojwry7bj66pl0oCmpLRayueZL5OtONJEt5uIv7uxcnUMqVnSloSwFrMLhMzMTHbt2sVvv/3Grl27yMzMtEYuqyvrSYoatcrUg6DIUQyiHCvrtqNSqehV35OHH36RaMcVHM2NwE2dwf3Zkzjy4yPsPBZXJjmEqGxKXCD89NNP9OrVC29vb3r27MmTTz5Jz5498fHxoWfPnqxatcqaOcs9tQry/u1BIF+GGIQwl6+rlufva8/5Fj+zMOtp8hQN7R3WU+dAF5b8tJTU69KuhLCkEhUInTp14sMPP2TIkCGcOXOGtLQ0EhISSEtL4/Tp0wwdOpRPPvmEzp07WztvuXVzD4KsgyBE6ahUKnrW96HvoLeIcfyO2LxQfDVJPHZ9NH8ue4rtJy7bOqIQFYb2zg+B9957j/bt2xd5X2BgIMOGDWPYsGHs3LnTouEqEo1aRbZiPOWzKl8KBCHuhq+rlnH39eSPo2s5uG8K9zktobfjcs7v3cW3x2fgYtDYOqIQ5V6JehCKKw5u1a5du7sKU5FpVJDLjQJBFn0R4m6pVCp6RFSl40Of8pXuaxINVamhTWB41uOQspadp1NsHVGIcq1EPQi3OnjwIPv37ycjI6PA7ePHj7dIqLIQHR1NdHQ0BkPZjFtq1CpyFUfjfyu5ZbJPIayhrNvOnfi6ahl1/4NsPNIaDkyim+MaHnX7hhP7/mTO6Zk81q09ns7SoyCEucwuEF577TVmzpxJ48aNcXFxMd2uUqnKVYEQFRVFVFQUaWlpeHp6Wn1/ahXk/jvEoFGuW31/QlhLWbedklCpVHRrGMTV0G9ZuPYbBuROo47uOKEZD7Lo+2eo3eF/dAz3sHVMIcoVswuEzz//nP3791O/fn1r5KmwjD0I/xYI+Tk2TiNExeTnpuXhe4Yy5wdv2ud9SUuH7Yx0nsWBnZv55MT7DO/WFA8n6U0QoiTMXgfBx8eHmjVrWiNLhaZR898QA1IgCGEtKpWKEDcNNfovY4V2KtfznWjmuI+RaQOY/90nbDuVbuuIQpQLZhcIH3zwAVFRUZw6dYq0tLQCF1E8tUplmqQoPQhCWJ+fm477H3iOvfXW8o++MW7qTMa7TCF/x1BmrTlKWrZ9zKEQwl6ZXSB4eXmxdu1a6tati7e3N97e3nh5eeHt7W2NfBVKLjd6EGSSohBlQaVS0alpE6oMWMdvmgnkKRo6O21kROo9fLLkW7afybjzRoSopMwuEEaNGsWoUaM4fPgwZ86c4cyZM5w9e5YzZ85YI1+Fov93iEGryGGOQpQlPw8nej84jd21fyY2rxY+miQmu07g2pYo3v/tlPQmCFEEswuEpKQkpk6dSv369QkJCSlwsQcffvghPXr0sHWMIulvDDEoMsQgRFlTqVS0b94e1/5b2KIZQb6i4h6XlQxJvpfpi1dKb4IQtzC7QHj00UdZvXq1NbLcNb1ez8GDB20do1h5/w4xaGWSohA24+fpRqcHP2RvzaUkGqpRXXuet1yfIHbTq0z/LU56E4T4l9kFwvnz53nkkUfo0KEDDz74YIGLrX377bc89thjto5RLFOBID0IQtiUSqWideveqPvsZK/6ftQqhcFu3zAw6UFeX7SeHdKbIIT56yC0bNmSli1bWiOLyfjx41m9ejWxsbEcOHCApk2bmu47efIkw4cP5+rVq3h6ejJ//nwiIiLIz8/n999/Z+nSpbz//vvFbjsnJ4ecnP/+QN84+kKv16PXW/ccCbkqZwAclCyz93Xj8dbOeDcko2WURcbSbNuWbccc5rx/nm7uNL53Lnv39KF2/MuE604yQzuErzY8zYZjo3m6UxWLrZtgr//2JJd57D2XJakURVEsvtW7tGXLFsLCwujYsSMrV64sUCB069aNYcOGMWLECJYtW8Z7773Hnj17WLZsGTk5OQwZMoQePXqwfv36Irc9depUpk2bVuj2xYsXF1gZ0hoSk/bxlNtHbM7uQYrnMKvuS4jbycrKYvDgwaSmpuLhUbIVBm3ZdspCnj6TsJxvaabbAcBfOc34KP01GnqrqOUiPQrCvpWmTd+J2QXCggULaNKkCU2aNDHddvDgQQ4fPszjjz9ukVA3hIaGFigQEhMTCQ8PJzk5Ga1Wi6IoBAQEsG3bNpYuXcqmTZvQaDTs3buX6dOn8+STTxbaZlG/goKCgrh48SK+vr4WzX+r+Us/5SmHqezJ7UDTR1aZ9Vy9Xs+6devo2bMnOp3OSgnvjmS0jLLImJaWhp+fn1lfJrZsO+a4m/dPyc/n2O55hF+Yhosqi6x8Fz5Jm4g+8DHGdfLF/S56E+z1357kMo+95kpKSiIgIMCiBYLZQwxTpkxhz549BW6rUaMGDzzwgMULhFvFx8cTEBCAVmuMrVKpCA4OJi4ujsmTJzN58mQAevToUWRxAODo6Iijo2OhE87odDqrf9h5ajdjBrJKva+yyHm3JKNlWDNjabZry7ZTGqXN1bjTWJISe3B581PUVO/nJa832XZ1MxO/n8JTXevRPszNJrmsTXKZx95yWSOL2ZMUk5OT8fPzK3Cbn58fSUlJFgt1t4obXrA1vcYdAGdVJnY4siOE+JevfzihD63jWNVX0StaOjptYZbLQH5Zu4R3fr9EuhzpICoBswuEmjVrsnnz5gK3bd26tUzWQbjRnZmXlweAoijExcURHBxs9raioqI4evRood4QazKoXQFwUWWiN0iBIMonW7QdW1CptdSLfIn0Thu4qKqDt+Yab/u8QNOLL/D0wiNypIOo8MwuEF566SUeeugh3nrrLZYuXcrbb7/NwIEDeemll6yRrwB/f3+aN2/OwoULAVi+fDmBgYGEh4ebva3o6GgaNGhAq1atLB2zWPkaY9ekiypLCgRRbtmi7diST/UmVHtgK2f8niZfUdHP5Sc+cHmYpb//zLtrpTdBVFxmFwiPPfYY8+bNY8eOHbzxxhts376duXPnWnT+wZgxYwgMDCQhIYHevXsXKABiYmKIiYmhTp06TJ8+nXnz5pVqH7b4FWQqENRZ5EqBIMqpytKDcDOV1omwbtNJbbuaZGpQTXuR2T5jqZnwFk8tPCG9CaJCMnuSIsCAAQMYMGCApbOYxMTEFHtf3bp12blzp9X2bU35WjfIAUdVDim5ueBSqrdfCGEj3iGdUQL+JH7TiwSlfMejbt/SWr+Dd359g821W/JM5yp3daSDEPakRD0IR48eLdHG/vnnn7sKU5Zs0U2qaP6b/ZyXK+ekF+VTZRtiuJXKwYOgXjGkNP+WdHwJ053mc7/h+MZ9zBMLz7DzbKatIwphESUqEB577DEGDRrEmjVryMws+I8/KyuLX3/9lUGDBtn1Mse3skU3qVbnQFa+cUGZvOvJZbZfISypMg4xFMUr/F7c7t3FJffe6FR5jPaI5v+cRvHJL7tlboKoEEpUIBw4cID+/fszZcoUvLy8CA0NpXHjxoSGhuLp6cnrr79Ov3792L9/v7XzlmuOWjWp+V4AGLLt57BQIUTpqJyqUK3PUtIafUI2rjRy+Juv/B5De24BIxeek94EUa6VaBBcrVYzYsQIRowYwcWLFzl48CDXrl3D29ubpk2bEhAQYO2cFnfrYi9lwdXBWCAEcIG861fLbL9CWJIt2o5dU6nwqD8MJagzSZvG4Jv1Jy96vc3O7M28//NrtKpbk2c6V0GmJojyxuxZcgEBAeWyILhVVFQUUVFRpKWl4enpWSb7dHFQk5Zv3Fd+tgwxiPLJFm2nPFC5heLb/1fSD32M87G3aOe0jXkOjzDz7KuMjOvFc5H2sxy1ECVh9mGOovScdSpSFS8AlFwZYhCiwlGpcW/8HJpem0lzaoCnOpU3vF/iSe1rvPvrGX6/Wp2MHOl5EeVDpS0QbDET28XhvzkIqhzpQRDlU2U/iqEkVF4RePTfTGbYc+Sjpo/LL3xd5REcc0/y5HcX+FPmJohyoNIWCLaYiW0sEIxdshq99CCI8kmOYighjQOuLaeh6vormbpgqmouM9v3aR5Vz+D11WeYvvYS6dKbIOxYpS0QbMFF918PgjZPehCEqAxUVdries9OMmsMBWCQ2xJi/IZx+uQeRi6Mk94EYbekQChDzg5qkg3GM2E65122cRohRJnRueHQejY7tc+TrfGjpu4Mn/sNpy9zmLw6TnoThF0q0VEMNWvWRKVSlXijzz33HOPHjy91qIrKzUHNlXx/AFwMl2ycRghR1hI1TdF0H0P2vhdwuvwTT7jPoa3jdt45MY2R8WG82N2fNqGuto4pBFDCAmHTpk1mbdTLy6sUUcqWLY7l9nTWcMVgLBDclKuQnwdqOR+DKF9kHYS75OiLU+dvUWK/I2/vi0Q4HOJLv8F8nv4cL616iL4NPBnX2Q93R1k4QdhWif46hYSEWDtHmbPFsdyujmrSFG8MigaNygDZieBSvUz2LYSlyDoIFqBSoQp9DF2VjuT++TTOSVt53nM6HRy3MOPY/7EnLkt6E4TNmf3z9Y033ihw3c3NjcaNG9OjRw+Lhaqo1CoV7s4OXM33o6rmMly/KAWCEJWZaxAO3VajnJxD/l9TaeO0g68dHmV26su8tKonfRt4SG+CsBmzJykeOHCgwOW3335j0KBBPPHEE9bIV+F4Omu4+u8wg5KVYOM0QgibU6lR1RmHptcW9B6N8VSnMtX7FSZ7vcbWfxIYuTCOXefkSAdR9szuQVixYkWh2xITE2nSpAknTpygTp06FglWUXk5a7iSYSwQ9BnncbBxHiGEnfCsh67nHyhHZ6D8M4uezr/R1GE/01Om8NKqNvRr4MHT0psgypDZBcKoUaOKvF2n0zF69GjCwsIA+Prrr+8uWQXl5azhQmogAPrUM1IgCCH+o3FA1eg1VNX7kLfzKapknWGmbxTLMx/hi6PPslvmJogyZPYQg6enZ5EXjUaDq6ur6bq9s9VysV7OGs7nBQGQn366TPcthCXIUstlwLcl2j7bUGo9CcBDrkuZW2UIvjl/8dKqC8xYd1nO6SCszuwehNmzZxe67dixY3zzzTdER0cTGhpqiVxWZ6uZ2P7uWnYbjD0ImqwzZbZfISxFjmIoI1pXVC1mQvW+GHaPI5hYPvMdxYKMJ1h4dBR74rJ4QXoThBWZXSA0a9aswKJJmZmZxMbG8tJLL5Wb4sCW/N20ph4Ep5wEyNeDWmfjVEIIuxXQA02fP1H2PY8mYQUj3b+greM23kl5g5dW5dHv3yMd3GRugrAwswuE5557rsD1G4c51q5d21KZKjR/dx1X86uQozjiqMqBrHhwC7N1LCGEPXP0QdV+PsTdQ/7e56nPUeZWGUJM2nhWHh1kXDehhz+tQ6Q3QViO2QXC8OHDrZGj0vB316KgJiEviFq6U5B2XAoEIUTJBA9E7dcOZU8UTpc3MsHzfdo7beG9lClMWplHvwgPxnWS3gRhGSUqEBYsWGDWRps2bUrjxo1LFaii83PVogLO5NU2FggpR6B6X1vHEkKUFy41UHX+EU59ifLX67Ry3MW8Ko8wO/Ul1hzpw55Y6U0QllGiAmHjxo1mbdTLy6vMC4Tdu3fzv//9D0VR6NatG2+99VaZ7r+ktBoV/u5aTuvD6ekMSuoRSn4aLCGEAFRqqD0aVdWuKLtG435tP697v0bH65uZnfqy9CYIiyhRgTBv3jxr57hrzZo1Y/v27QB0796dtLQ0PDw8bJyqaME+Dpy5ZJyzYbh22PxxHiGEAPCojar7WvhnJsrRGXRzXkcjh4PMSHmdNUfaSW+CuCtmr4Ngr3Q645EABoOB6tWr4+LiYuNExQv2duC03lggaDJOgSHbxomEEOWWWgcRL6Pqvg7FvTZVNFd43/dZJni8R3pmOpNWXuD99bJugjCfXRYI48ePJzQ0FJVKxcGDBwvcd/LkSdq3b0+dOnVo1aoVR44cMd23ePFi6tevj5eXF1qt/f4uD/bWkZTvR4rBCxX5kHLY1pGEEOWdTwtUPbdA+BgAHnD9gbl+Q6ivO8wvR9IYuTCOPbFyTgdRcnb5V3TgwIFMmjSJjh07FrpvzJgxjB49mhEjRrBs2TJGjBjBnj17ABg8eDCPPvoogwYN4tChQzRq1KjQ83NycsjJyTFdT0tLA0Cv16PX6630igqq4aEGVBzVN6K9ZiuGxJ3kezS57XNuZCurjKUhGS2jLDKWZtv20HZKwl4/47LJpYNGb6Oq2hPNvmcJyo7jU98nWJgxkgUZT/Liygv0re/GmA4+uDqoyzCX+SSXeayRR6UoimLxrVpIaGgoK1eupGnTpoDxpFDh4eEkJyej1WpRFIWAgAC2bdtGUFAQjo6OAIwYMYLJkycXuTbD1KlTmTZtWqHbFy9eXGbDElkGDTHxdRni+jVPeXzGeXVr9uqeKZN9CwGQlZXF4MGDSU1NLfFcHXtoO6LkdEomjfIWEJS/E4DjufV5O/UN4vJq4q7R09PvAiHO0qNQUZSmTd9JuSoQ9u3bx+DBgzl+/LjpMa1bt2b69OkkJSURHR1Nfn4+nTt3LvYohhu/gubOncvcuXMxGAycOnWKixcv4uvrWxYvC4AhC+KpnruLD33HojjXIK/PX7d9vF6vZ926dfTs2dM038LeSEbLKIuMaWlp+Pn5mfVlYi9t507s9TO2VS5Vwko0B19ApU8hR3Hki7Rn+DHrERTU9K3vxhOt3dm++Q95v0rIXnMlJSUREBBg0QKhREMMxZ3B8Va2PIPjoEGDGDRo0B0f5+joiKOjIxMnTmTixImm9eR1Ol2Zftj1qjmz53QDDIoazfXz6PSXwSXwjs8r65ylIRktw5oZS7Nde2k7JSW5/lVzEFTrCHuicLz0B896zjQtrvTrP9XYG3+dTq6u8n6Zyd5yWSNLiSYp3nzWRrVazeLFi7l06RKOjo5cvnyZJUuWoNFY/1jboKAgLl68SF5eHgCKohAXF0dwcLDZ27L1GenqVXXkuuLCCX094w2JW2ySQwhz2brtiFJwDoBOy6H5TBSNCy0c9/B1lUfp6byGKxl5/Hg5hNmbrpIpRzqIm5SoQJg9e7bpkpKSwrJly1izZg2ff/45v/zyC8uWLSM5OdnaWfH396d58+YsXLgQgOXLlxMYGEh4eLjZ24qKiuLo0aOmCY5lrV5VJwD25rQ13nDpD5vkEMJctm47opRUKgh/ElWvbeDTEjd1BpO9Xmeq18t4qlJYczSDkQvj2CtHOoh/mX2Y49q1a+nXr1+B2/r06cO6dessFmrMmDEEBgaSkJBA7969CxQAMTExxMTEUKdOHaZPn17qRZxs/SuoXlUntGrYk/tvgXB5Iyj5NskihDls3XbEXXKvBd1+h4avoai0dHH+g6+rPEIbx20kZuTxwsoLzPzjsvQmCPMLhNDQ0EJzDebPn09ISIjFQsXExJCQkEBeXh6XL1/m1KlTpvvq1q3Lzp07OXHiBHv37i3yUMaSsPWvIBcHNQ2qOXE0txFZ+S6QkwQpf9skixDmsHXbERag1kKDF1F1X4/iXgdfTRLv+TzH8x7v4KzK4qfDaYxaJL0JlZ3ZBcKnn37Kiy++SMOGDenfvz8NGzbkhRdeIDo62hr5rMYefgW1CHYhDx37c1sab7jwu82yCFFS9tB2hIX4NCOv6x+c0vQG4F7XH/mqymAa6P7mcrr0JlR2ZhcInTt35uzZs7z88stERkby8ssvc+bMGTp37myNfFZjD7+CmgcZjx3flt3FeEPCKptlEaKk7KHtCAvSOHNEO4S8jj+CSyDVNQl86vckT7hHo0X/X29CXJatk4oyVqqVFL28vHj88cctnaXSqV/NCQ8nNduzI8lTNGhTj0DaSfAovMCTEEJYk1KlM/TaDgdeRh27hKFu82jntIO3r03jbHo4L6w4z4CGHozt6IernCGyUihVgbBz5042bdrE1atXuXmdpVmzZlksmLVFR0cTHR2NwWC7rjOtWkXHWm6sOZLPvpw2tHHaAQkrocGLNsskxJ3YQ9sRVuLgBW3mQI2+sPc5wjnOF1WG8UXaOJZlDuanw2nsjs3ixR5VaRksq2dWdGYPMURHR9O9e3d2795NdHQ0Z8+eJSYmhkuXLlkjn9XYSzdpZLgbAJuyuxtviP0e7HdxSyHspu0IKwq8D3r/CQG90JFLlMeHfOT3NFU1F41zE1acZ+Yfl8nKlSOvKjKzC4QPP/yQX3/9lRUrVuDs7MyKFSv44YcfTOdBEOZpHuSCh5Oazdndua44QfoJSNpt61hCiMrOuSp0/B5afAhaVxrr9vGN/6P0cf4JUPjpcBojF8bK3IQKzOwC4fLly0RGRhqfrFajKAp9+/Zl9erVFg9XGeg0KnrV8yBLcWPT9Z7GG88ssG0oIYQA4+JKtUZCr23g2wYnMnnZaxpv+0zCU33N1Jswa0Oi9CZUQGYXCNWqVePChQuAcU2ETZs2ceTIEdRqszdlU/Z0qNY9DY0n1vg5634AlPgfITfFdoGEuA17ajuijLiFQddfodEUUOvo4LiRhVUfpZ2jcYn41YdSGbUolgPx0ptQkZj9V/3pp59m165dADz//PP06tWLZs2aERUVZfFw1mRP46ihvo40DHDiiL4xp/XhqAxZcNp2J74S4nbsqe2IMqTWQP3nofsG8KiPO0m86/M8L3u9ibMqk0tpefzvx/N8uFF6EyoKswuE//3vfzzwwAMADBkyhDNnzvD3338zdepUS2erVAY28wJULM389/DRk3PAkGPLSEIIUZh3Y+i5Ceo+C6jo47yKb6sNoZHuIAAr/07liUWxHEyQ3oTyrlTjAnl5eWzfvp2lS5cSFBREcHAwmZmyJOfd6FTLjSAvHRuu9+aKwR+yLxuPaBBCCHujcYImb0GXn8ElGD8S+NjvKcZ5foKOXC6m5fHc8vN8svkK1/XSm1BemV0gHDt2jPr16/PII4/wxBNPAPDHH3/w5JNPWjycNdnbOKpGreKxlt7koWNZ5mMAKP98AIZcGycToiB7azvChvw7Qu/tEPo4KhQedvmG+QEjCNOeBGD5wRSeWBTH3+ev2zioKA2zC4Rx48YxYcIEEhIS0Ol0AHTp0oVt27ZZPJw12eM4aq96HgR56ViVNZBkgy+qzHNw5htbxxKiAHtsO8KGdB7QOho6LAJHP2pwgrn+w3ncfQFqDFxI1TNhWQLRW66QLb0J5YrZBcLBgwcZN24cACqVCgAPDw/S09Mtm6wS0mpUjOnoR7bizDcZxh4Z5egM0GfYOJkQQtxBjXuMiytV74dGyeVJt4/5ImAc1TTnUYAfDqTw5OI4Dl+Q3oTywuwCoWrVqpw7d67AbSdOnCAwMNBSmSq1DmGuNA105pes+zmfVwNVTiIc/8jWsYQQ4s6cqkCHxdDyU9C6Ec4+vq06mAGuqwCFhBQ9z/6QwOdbr5CTJ70J9q5Uhzk+9NBD/PTTTxgMBtauXcvQoUN55plnrJGv0lGpVDzXpQoqjQNfpD8LQP4/H0L6KdsGE0KIklCpIGyo8cRPfu3RKZlM9HiTj6u9iLc6CQVYut/Ym3DkovQm2DOzC4Tx48fz1FNP8corr2AwGHjuuecYOnQoY8eOtUa+SinU15ERbXzYnN2d3dntUCu5qA5OknM0CCHKD7dQ41EOjd8AtQONVZtYEjCYSOdNAMRfM/YmxGy7Kr0JdqpUhzmOGzeOw4cPk5GRwdGjR8tl74G9z8R+pLk3dfyd+ChtErmKA9qrWwgybLV1LCHsvu0IO6LWQL0J0GMjeDbEKT+JaV4v8E61N3FRZZCvwJJ91xi9JJ5jl7JtnVbcolQFQmZmJgcOHGDLli0FLuWJvc/E1mpUvN6nGimaEOaljwYgIm8RZMbaOJmo7Oy97Qg75NUQemyAes8BKtqrVvF94BBaOO0HIDY5l3HfxzN3+1VypTfBbmjNfcL333/Pk08+SW5uLs7OzqbbVSoVycnJFg1X2QV6O/Byz6pM+WUo7Zy20djhIGk7xqLr87uxMhdCiPJC4wiNp0FAH9g9BrfMWD7wHsO6/GF8cHk0uYoji/ZeY/vZTF7pWZUwH/mOszWzexAmTZpEdHQ02dnZXLt2zXSR4sA6OoW78XALP95JmUZWvgseGXtI2j3V1rGEEKJ0qrQzTmCsOQwVCr3U37A0eBR1HY4DcC4pl6eXxjN/1zUMMu3KpswuEFJTU3n88cetkUUU46n2voSH1OaD1MkA+MZ9TOI/y22cSgghSknnDq0+gY7fgWMVvPXHmeM3gvH+36LGQL4Ci/alsvhCGKeuyDlpbMXsAuHxxx/nxx9/tEYWUQyNWsXLPfw4pnRiacYQANz/fob4c3/ZOJkQQtyF6n2NiyvVuAeVoudBzUd8FzqOQF0CAFf1Tjyz/CLz/kxCL90JZa5EcxAeeOAB06qJeXl5DBkyhJkzZ1KtWrUCj5PCwXoctWru849jQ/ZE6uQco5njPpz+fJhT6l8JDw6zdTwhhCgdJz9ovxBil8D+Sfjn7GNB1cF8mzOReVfuxZCv4ptdyWw/k8nLPasSXsXR1okrjRIVCE2bNi1wvUWLFtbIcle2bdvGCy+8gFqt5qGHHmLixIm2jmRxTpp8pt9bg7d//hDfvMcJ1saSsW0gu1qupk0dWclSCFFOqVQQOhiqdITd41Bf2cpw7Zu09v+d15LeIMngx6krOYz5Lo7hrX0Y3MoHrVpl69QVXokKhClTplg7x10LCwtjy5YtODg40LVrV55++mlcXFxsHcvi3Bw1TH0gghmr5jLeMJia2tOk7x7C6szFDGha3dTTI4QQ5Y5rMHRZDSc+Qzn0BvXZzdLqjxFz/TV+uBqJIR++/jOZHWczeaVXNUJ8HGyduEIr1ToI9qh69eo4OBj/sWg0GtTqCvPSCnFz1DD5wVb84P4VGfluNHY4SOiRwcxae1rOliaEKN9Uaqj7DHld/yBFFYw27xpRuonMrfUW7mrjieuOXc7hqcVx/LD/GvmywqzVmL0OQlkYP348q1evJjY2lgMHDhQY4jh58iTDhw/n6tWreHp6Mn/+fCIiIkz3r1u3jlq1auHk5FTktnNycsjJ+W9WbFpaGgB6vR69Xm+dF2QBN7Ld+H81MKpXe37aMo8+V0fS2OEgmisjeP67z3mxTy0CvXQ2z2iPJGPBfZijvLQde/2MJZd59M612KKbSt/gg2hPfUrtrJUsD9rNzPSp/J7cnFyDQvTWq2w5nc6L3fwI8Cib7zy7fb+skEelKPZXfm3ZsoWwsDA6duzIypUrCxQI3bp1Y9iwYYwYMYJly5bx3nvvmVZ0S0hIYNiwYaxevRo3N7citz116lSmTZtW6PbFixeX2yGJlIxEHtC8gYc6jeO59XktZSaNPA00dr+GjDiIomRlZTF48GBSU1Px8PAo0XMqYtsR5YN3/kma62NwIxGAzfoHePvqRHIx/hDUqQxE+lymoVtKpf3OK02bvhO7LBBuCA0NLVAgJCYmEh4eTnJyMlqtFkVRCAgIYNu2bQQFBdG/f3+io6OpW7dusdss6ldQUFAQFy9exNfX19ovqdT0ej3r1q2jZ8+e6HSFK+WLcfvx2vsIHqprXMqrxsvXPsI3IIKJXf2o4lY2HUV3ymgPJKNRWloafn5+Zn2ZlJe2Y6+fseQyT6FceRmoD09Fc3Y+ANedajMt+Q3+TK1tek6rYGee7+qLn6v1vvPs9f1KSkoiICDAogVCqd7FjRs3snjxYi5dusRPP/3E3r17SU9Pp2vXrhYJVZz4+HgCAgLQao2xVSoVwcHBxMXFsXXrVo4ePcqYMWMAWLRoETVq1Ci0DUdHRxwdHYmOjiY6OhqDwQCATqezqw+7OMXlDK7Vhuve67m24SGqac/xqe8oplyewVPfteOp9r4MaOSJpoxm/ZaH97KyZyzNdstb25Fc5rH7XDpvaPURBPaHPc/gnH2Sd12Hss33GaaeeRQDWvbEXWf0dxd4rqs/3eq4WXXStr29X9bIYvZMvi+//JKhQ4dStWpV0wmadDodr7/+usXDmWPkyJFcuHCBTZs2sWnTpiKLg5tVxBPOOPuE4z1gI2nubXBTZ/Kez3j66r7lw02JPL00nn/kbGnCAipi2xHlSEAv6L0TAu9DpeTR6fqHrAwfR0N34+JK6Tn5vPnbJab9eomU6wYbhy3fzC4QZsyYwdq1a3nrrbdMRwo0aNCAf/75x+LhbnWjOzMvLw8ARVGIi4sjODjY7G1V2FPWOvrg0esn9DUGoVUZeMZjNtO8XiLhShLjlsbz/vrLXM3Is3VKUY5V2LYjyg9HX2j3DbSZCzpP3DP384nnYCaF/gQYR803ncxg5MJYdpzJsG3WcszsAiEpKYkGDRoAmLpvVCpVmRx/7+/vT/PmzVm4cCEAy5cvJzAwkPDwcLO3VaF/BWkc0bWfC83eJ1+lI9J5AzF+w6ijO8IvR9J4/JtzfLnjKhk5Ul0L81XotiPKD5UKQh6G3jvAvzMqQxb9cqaxvM6LhDonAXAty8CrP11k+tpLZMr3ndnMLhCaNGnC8uUFTxS0evVqmjdvbrFQY8aMITAwkISEBHr37l2gAIiJiSEmJoY6deowffp05s2bV6p9VPhfQSoV1B6NuttvKM6BBGnj+Mx3FMPc5qLP07NwzzUGzz/HDweukSPnXxdmqPBtR5QvLoEQuQqavgcaJ3zTN/F1lUcZHbzF9JDf/knnycVxHLpw3XY5yyGzC4QPPviAMWPG8NBDD5kOqxg3bhwzZsywWKiYmBgSEhLIy8vj8uXLnDp1ynRf3bp12blzJydOnGDv3r00atSoVPuoNL+CfFui6rUVgh5AozIwyj2GT32fJFh7lrTsfKK3XOXReef4bt81snKlUBB3Vmnajig/VGqoMxZ6bgHvpqj11xisf57F9d7E39E4xHAxLY8JyxL4amcSeXLipxIxu0Bo3rw5R44coV27djz55JM0adKE/fv306RJE2vks5pK9SvI0QfazoM2X4LOk/oOh5lXZTAj3GJwIIdrWQbmbLvKo/POsmB3MunSFSduo1K1HVG+eNSF7uuhwUug0lA9dRVLAh5jYMABAPIV+HZ3Ms/8EE/8tVwbh7V/ZhUIeXl5BAQE4OnpyQsvvMCnn37KSy+9RPXq1a2Vz2oq3a8glQpCBhnH6wJ6o0HPCPe5LKk+mOYOxvcgLTufr3cm8fBXZ/l4UyIJ0oBEESpd2xHli1oHDV+FbmvBPRxNzkWe4Snm1PkEF7XxSK4bSzX/dCgVO14KyObMKhC0Wi1ubm7k5sofjnLLJRA6LoX2C8CpGr5KLLN8n2ZO0KsEao2HCV3XK/z4VypDF8TyyuoL7IvLkkYkhChffFtCz60Q/hQA9dK/YUXoCDp7HwcgO09h5oZEJv90kWtZcmRXUcweYpg8eTIjR47k0KFDpKSkkJaWZrqUJ5W6m1SlgsD7oM9uCB8NKjX18tbyrf/DfBD2OT4645idAuw8m8nEFecZuSiOHw+mkJ4tww+VXaVuO6J80bpA8w+g84/gHIDj9VNMcx7Gu2HfoMFYFOw4m8mohXH8eTbTxmHtj9kFwqhRo1ixYgVNmjTB19cXb29vvLy88Pb2tkY+q5FuUsDBE5q/D712QLXuqJRcWl7/imXVH+TDiB8IdPvv5B/nknL5ePMVHvryLO+uvcShC9elV6GSkrYjyp1q3Y3fc0EPoVIMtLv+CStrjSXCNQ6Aa9cNvLz6Ap9uvoJeJjCamL3U8tmzZ62RQ9iSZ31jhX1xPfw1GXXaMZomv8e3fvM4VXccn124hwMXjQ/NNSj8/k86v/+TTqiPA/0iPOhe1x1fK659LoQQd83RB9p9DTX6wb7ncc86yKdeg/nJ53lmx9+PgpplB1P4+8J1pvStRg0vB1sntjmzexBCQkKKvZQn0k1ahIAe0Gs7tJ4DrqGochKpfWEqs50HsKzzrzzaWIeb43//ZM4l5/LZ1qsM+uosL648zx8nMtDnV9JTqVUi0nZEuRY8EHr/CVW7ocrP5t68d/iu1vME6IxnijyRmMNTS+L543i6jYPantk/+0aNGlXsfV9//fVdhSlLUVFRREVFkZaWhqenp63j2A+1FkIfMzaic4vh6PuQFY/fyf9jrMNMnmz3JNvUQ1h+TMehC8YZwfkK7InNYk9sFjpVHY7/cYXeDbxoFuhcZieIEmVH2o4o91yqG3tNT30Jf/8fVbO2sbDaY8zJeoUfknqQlWs8n8P++CyejayCk87s39IVgtmv2tPTs8AlKyuLZcuW2dVZrYQFqHUQNhz67oMWH4JbLchNQXv8A7qcaMcnYbNY/FAuQ1v7UM3jvzpTr2hYdzyTF1ac55Gvz/LJ5iscTMjCkC/jekIIO6JSQe2njEc6+LRAk5dKlMPLfBb8Ou6qVAB+OZLG2KXxxCZXziP3zO5BmD17dqHbNm7cyNy5cy0SSNgZjSPUGgk1h8GFX+DYh5C8D05/TfXTX/OEfyQjezzBYVU3fvsngz/+SSVH0QBwNdPA8oMpLD+Ygqezhg41XekU7kbzIGcctZWzIhdC2BmP2sY1E/6ZCUffo4F+DT8E7uXNq//H9uvtOJeUy9jv4ni5Z1Uia7vbOm2ZssjMsi5dunD//fdbYlPCXqk1EHgv1BgAV3fA8U/h4m+QuBl14mYaO1cnInQYDdJr4NLgXv44mcWuc5ncOM1D6nUDa46mseZoGs46FW1DXelYy422oS64Omps+9qEEJWbWgsRL0FAT9g1Bqf0E7zt9Swb3B9mxpVnua53ZsqaSzzaIocRrTxsnbbMmF0g3LreQVZWFt988w3VqlWzWKiyEB0dTXR0NAaDHNdvFpUKqnQwXjLj4PQ8OLsArl9A8890+qGGxJ/p1uxxMnr05M9YPVtPZ7IrNpNsvXGY4bpeYePJDDaezECnUdEiyJmOtdzoEOaKt4scDWHvpO2ICsunufF8Doemwsk5dNN8T5Pqu/m/xCkc1Tfiu33XOHbpOm3UleNHjdnfxl5eXgVO7awoCiEhIeVqgiLIRCuLcA2GxlMg4mVIWEX+ybmok3fDpd/h0u+4OfjQI+RherQfTE7vRuyNy2LrqQx2nM0kLdvYtaA3KPx5Los/z2UxawM0DHCiUy03Ooa7EeAh81rskbQdUaFpnaHZe1C9D+yOwvf6OaL9nmBhxkjmpz/FwfPZnNKE0Sgxh4Y1KvZ31F2vg+Dm5oavr6/FAolySOMIIQ9jqP4AG3+eS5fQi2jiv4fsS3ByDpycg6NnBB2CHqRD+wfJ6xHGofPX2XI6g22nM7mSYVzRLF+Bvy9k8/eFbKK3XqV2FUc61TIORdT0dShQmAohhFVV7Wo8d83+F1HFfc9Qt69o77yDN5KnEZsXxsSVl5jcuxqdw91sndRqzC4Qbl3v4MCBA1y4cKHUp10WFUuGugb5DZ9C02QKXN4I5xbB+V8g9YjxcvhNtN5NaRb0IM1aPcD4yFCOX85h6+kMtpzOIP7af6s3nrySw8krOXz9ZzI1PHV0DnejYy1X6ldzQi3FghDC2hy8oO3cfxdX+h+1cv/hS7/HiUl7huVZj/L6LxcZ3cGXx1p4V8gfMGZPJb/vvvvYtm0bYByLbN++Pe3atWPOnDkWDyfKMbXWOOGn3XwYcBxafmKsyFUauHYQ/n4dfmmEakNP6qXO5anGqXw7LJRvhobwZDtf6vo7Ftjc+VQ9S/ZdI+r7BB7+6iyzNySyNy5LzusuhLC+oAeMiytV64lOlcsznrP4wCcKf/UlvtiexPvrEyvkEs1m9yDs3LmTNm3aAMYCYf369Xh6evLAAw8wduxYiwcUFYCjD4QNM16yr0LCKohfAVe2QdIe4+Wv18CjLiHV+xES1o/HW7XkcrqBbWcy2Xo6g7/PX+fGUgpXMw2sOpTKqkOpuDuqaVfTlU613GgV4lJpFzQRQliZczXo9AOGk3NRDk6mheMevq7yCB+lTWLN0X5cycjjjXsCcK5A30FmFwg5OTnodDrOnz9PcnIyHTp0AODy5csWD2dNMhPbRpz8IPwJ4+X6ZUhYDRd+hsRtkHbceDk2GxyrULV6Hx6q3o+HIjqTkhfAzjMZbDmdyd64LFO1np6Tz9pj6aw9lo6jVkXrEBc61XKjXU1X3J0qx0zjsiZtR1RaKhX5NUey+YhCN7fvcbu2l8leU+jouJlZ8a/ywop83r23Oh4V5LvH7AKhQYMGvPvuu8TGxtK7d28AEhMTcXV1tXg4a5KZ2HbAuapxJbPaT0FuClxaD+fXwKV1kHMFzn5rvKh1ePm2pW+1bvTt1J0slwh2x11n66kMdp7LIivXeERETp7C1tOZbD2diUYNTWs4G4+IqOWGn5scPmkp0nZEZZeproah88+oT0ejHHmXSOcNNHT4i/eT/48Jy7rx/v01KsR3jtmv4LPPPuOZZ57BwcHBdGjj77//Tq9evSweTlQiDl7G8z8EDwRDrnExpvNrjIsxZcbCla3Gy6FpuDhWoUvVrnSp353cTl04cNWdracz2H46k2vXjb9qDfmwL/46++Kv8+GmK0QEONE53I32IY63zyGEECWh1kKDF1AF9CBn+2h8s44z3ec5Vmc+wMvLJ/HewDrl/iy3Zqdv1qwZ27dvL3Db0KFDGTp0qMVCiUpO4wBVuxgvzID003BpA1z+AxK3GHsX4r6HuO9xANp4NqBNlY78r08HjuW3YFOcE1tPZ3ApLc+0ySMXszlyMZvPt4K/Q01S96XQta4nQd5ySlchxF3wbopj3y2k752K67k53Ou6ghZ5e/h8xTtEPdi/XC/+ZnbyNWvWUKdOHcLDwzl37hwTJkxAq9Xy0UcfERgYaI2MorJzr2W81H7K2LuQtAsu/WG8pPwNqUch9SiaU18QAUR4NmBckw5ccmzLxmuNWXfOkbNJ/51sJTHXma93pfD1rhTCfB3oHO5G53BZa0EIUUoaJ9zbTCfZvzeGXU9TQ5vAK8oIfl01kk73vo2nq7OtE5aK2QXC888/z/r16wF48cUXcXNzw9XVlaeffpqffvrJ4gGFKEDjAP6djJfGU41HRVzdbpzkeGWbqVhQpR4lgLkMBgb71Sc9rA1/Zzfil0t12XGlKmAsBM4k5XImKZn5u5IJ9XWgRx13utd1J8CzYq+QJoSwPJ+aXbnkvp1Nv0+gi+4n7tF8Rfwv23Hu/hUOvg1tHc9sZhcIly5dIjAwkLy8PNavX09cXByOjo5Ur17dGvlKLDk5mZ49e3L8+HEyMjJsmkWUISc/CLzPeIGiC4a0f3BP+4cOQActXA/w4KpLa/ZmNmZDUn2O6+uTixPnknL5cmcSX+5MIiLAie513ela261cdxEKIcpWNT9f8vp/zXsrvmGM01sEaY6R90cX8hu/jrruM6AqP4dBmv3N5+zszOXLlzl06BD16tXD3d0dvV6PXq+/85OtyN3dnXXr1vHwww/bNIewsSILhh1wdRck7Ua5dhDn/DSCstYTpFrPA36Qj5Y46rM3syFHcxtyXN+AIxcDOXIxm083X6FVsAt9IzzoEOaGTiNDEEKI2wv0cuD+fkN5enkjnnV7i/ZOW+Hv/4OLv0Prz8A15M4bsQNmFwjDhg2jVatW5OTkMG3aNAD27t1LWFiYxcOZQ6fT4ePjY9MMwg45+RlPUx14LwB52ens/DWGDnW1aK7thaTdqLMvE8ohQl0Pwb9H66bnu3NCX59j+gYcv1yfzxMa8KFDDXrV96B/hCchPjK5UQhRvLpVnXi2T0Ne+3kWfZ1W8YzHTJyvbIPfOxhPBhU62Hh2XDtmdoHw3nvv0aNHD3Q6HV26dAGMf5xnzpxpsVDjx49n9erVxMbGcuDAAZo2bWq67+TJkwwfPpyrV6/i6enJ/PnziYiIsNi+RQWnceKaujb5tfuh0elAUYyHUSbtMU5+TN4PKYdxJ50Wjrtp4bjb9NQUgxfHzzZg04kGZLs3oU7dNrRvEI6jrmIsiiKEsKz2YW483bEK0VvvZ39uS17znkoEB2HPOLiwBlp8ZPwRY6dKNbjas2dPzp8/z59//knbtm1p2bKlRUMNHDiQSZMm0bFjx0L3jRkzhtGjRzNixAiWLVvGiBEj2LNnT4m3nZOTQ05Ojul6WloagF0Mk9zOjWyS8e4UmdGxBlSvAdXvN17Pz4W0Y6iuHUSVcpDsxAM4Zf2DlyaFNpodtHHaYXzccUj7x4urzvXxqtYYR99GKJ4R4FEX1KXvYSiL97E02y4vbcde/x1KLvNUlFz3NXTl7/NZbD0TyLNXY3im6hIe0ESjOv8zytVdGJrNQgnoa7FclqRSFMWsM0xcvHiRRx99lO3bt+Pk5ERGRgbff/8969atY+7cuRYNFxoaysqVK009CImJiYSHh5OcnIxWq0VRFAICAti2bRvh4eEA9OjRw3SURVGmTp1qGhq52eLFi3FxcbFoflFxqJVcHA0XyM2+iLMhllDNCUK1Z9GoCi83nI+GdFV10lRBpKqCSVMHk6oKJlflYYPkRcvKymLw4MGkpqbi4VGyXNJ2hCidnHw1iy6EkZpn/OHwkPdOhjnPwFOJByBW3ZnD2iHkqUp/OGRp2vSdmF0g3H///dSpU4c33niDgIAArl27RlJSEq1ateLMmTMWCXXDrQXCvn37GDx4MMePHzc9pnXr1kyfPp1u3brRo0cPDhw4QLNmzfjwww9p2LDwYSU3fgXNnTuXuXPnYjAYOHXqFBcvXsTX19ei+S1Jr9ezbt06evbsiU5nn4fgVZaMiqLwz+UcfvrrCvFxR6ipPUG47iS1tCcI153AXZ1e9POcqqK410PxqIviXhc86hn/38HL4hnvJC0tDT8/P7O+TMpL27HXf4eSyzwVLdexyzlM+PEi+QroNPDFID+CL3yA+mQ0KhQUlxAMLT5F8WtXqlxJSUkEBARYtEAwe4hh+/btLFu2DK1Wa1pUxtfXl6SkJIsEuhu36zm4wdHREUdHRyZOnMjEiRNN68nrdDq7+kdYnPKQszJkbBLkQJMgdxLTg1j5dzu+OZRKelo+oOCvvkwt3Ql6VYulrfdZnDOPQMYZVNmXUWVfhiubC27MqZpxWMKzPnjUReVSG52SYdX3sTTbLW9tR3KZR3KZx9xcjQJ1PNTUix8OpKA3wEfbMpn94FuoAvvB7rGoMmPRbr0X6o6HhpNBY96y8NZ4j8wuEDw9PUlOTsbf3990W0JCAlWrVrVosKIEBQVx8eJF8vLyTEMMcXFxBAcHm70tOSOdsAR/dx2jO/jxeCsfVv2dwvf7U0i8Xo3EnGrsjAVdgooHGnsytLMD7jknIPUYpN24HIeseMi+ZLwkGgsHLdAPUNZMBc964NsK6jxjPG22HZC2I0TpjGrny7bTGVxMy+NgwnV2nM2kQ1gH6LUdDr4CZxfC8Y+MJ67r+B24mv+3zZLMXrFh8ODBDBkyhEOHDqEoCidOnGDs2LGMGDHCCvEK8vf3p3nz5ixcuBCA5cuXExgYaJp/YI6oqCiOHj1q1gRHIYrj4qDmsZY+LBkZSlRnPzydjUc26A0K3x9IYdiSJDYk1UOp+Tg0fRs6L4d7DsMDCdB9A7SKhrrPQrWeKC5BAKhyEo3nnvhnJqxtD/qihy7KmrQdIUrHWadmXKcqputfbE8iL18BnYfxO6DDYnDwgdQjcOhNGyY1MrtAeP3112nWrBnt27cnNTWVFi1aUL9+fSZNmmSxUGPGjCEwMJCEhAR69+5doACIiYkhJiaGOnXqMH36dObNm1eqfURHR9OgQQNatWplqdhC4KRTM6iZN4uGh/B4K28c/l1Y6VqWgTd+vcTLqy5wJeO/k0ihcwffFlDzcWjyFnReRl7vA/ziEENel7XQ8lPj465fhOsXbPCKCpO2I0TpdazlSkSAEwCxybms/Sftvztr9DeukQDGw69tzKwCIS8vj5iYGN544w3S09NJTEwkLS2N999/H63WcsvRxsTEkJCQQF5eHpcvX+bUqVOm++rWrcvOnTs5ceIEe/fupVGjRqXah/wKEtbk5qjhyfZ+LBweQsdarqbbd8Vm8cSiWHacuf1y4HkqZxTv5hA2FBy8rR3XLNJ2hCg9lUrFmA7/rX3w3b5r5N98rIDWfo4IMqtA0Gq1TJ48GScnY/Xj5+cnZ78T4jb83XW8dU913uwfgJ+rcdghLTufV3+6yJxtVwt+MQghKoXGNZxpUsN4SGPcNT1/ns20caKimT3E0Lt37xIdLWDvpJtUlKVO4W58NSSEDmH/9SZ8t+8aU365SLY+34bJzCdtR4i790hzL9N/LzuQYrMct2P2uIC7uzv3338/vXr1Ijg4GLX6vxpj1qxZFg1nTVFRUURFRZkO1RLC2jydNbx1TwAr/krl0y1XyFdg6+lMXlhxnhn318DFoXyc5U3ajhB3r21NVwK9dCSk6NmfcJ2LaXoCPOzrcE6zv5EMBgODBg3C09OT1NRUrl27ZroIIW5PpVLxYFMv3ruvOq7/FgSHL2YzaeV5snLLV0+CEKL01CoVfRv8t6BRgcmKdsLsHoRZs2bh7W1fk6ZKQ47lFrbUKsSVDwcGMvHHBNKy8zl8MZs3fr3I2wOq2zraHUnbEcIyetZz58sdSSjA7/+kM6y1D/Y0q6/EPQj79+8nKCgIPz8/atWqxZEjR6yZy+pkJrawtdpVHJn5YCDujsZm+Oe5LGK2X7VxqjuTtiOEZfi762gaaJyseCFVT2xyro0TFVTiAuGFF15g0KBB/P333wwYMIAXX3zRmrmEqBRqV3Hkjf4BaP5tid/vT2F3bJZtQwkhykz7myYubz9jX0czlLhA+Pvvv5k+fToRERG88847HDx40IqxrE9mYgt70SzIhaibVlebuTGJbIP9TliUtiOE5bSv+V+BsMPODncs8bdQbm4uDg7GU1W6uLiQnZ1ttVBlQbpJhT15oIknbUKMC6QkZxnYnep3h2fYjrQdISynhpcDIT7Gv63/XMomR28/a6OUeJJibm4uH3/8sel6dnZ2gesA48ePt1wyISoRlUrF8938GbogllyDwoE0Xy6k6gnxs6/DnoQQlte0hjOxybnkKxB/LQfzzy5kHSUuENq2bcuKFStM19u0aVPgukqlkgJBiLtQ1UPHw829WLjnGvmo+G5/Ki/1sp9lV4UQ1tGouhOrDqUCEJusL38FwqZNm6wYo+zJoVrCHj3Wwpsf/0ohK1fhjxOZjOlkwMvWoW4hbUcIy2r877LLAOeu2c+RDPY7E8rKZBxV2CNXRw1967sDkGtQ+O2o/S2eIm1HCMvyd9fh72b8vX4hRW/jNP+ptAWCEPZqQEN3039vOJFuwyRCiLJSq4ojADl59jNJUQoEIexMDU8d/g7XATiRmEO+/XxfCCGspKavg60jFCIFghB2qLbLf0MLeQapEISo6MKkQBBClESw838LpuRJF4IQFV6Yn6OtIxQiBYIQdsjfIRsXB+NpWwxykkchKrxAL/tb86TSFgiyXKywZ2oV1K9q/EVhb/WBtB0hLM9Bq8bPVWPrGAVU2gJBDtUS9q6mj/2NSYK0HSGspZrHf70IBsX2Q4uVtkAQwt6F+tpfl6MQwnpuLhD0djA5WQoEIexUqLd99iAIIawjwPO/AiFXCgQhRHH83Uu8EroQogKo4vZfm7eHw5ulQBDCTnk6q9FpVLaOIYQoI94u/01SzLOD2ckVqkB49tln6dSpE++8846towhx19Qqld3NahZCWI9PgQJBehAsZu/evWi1WrZu3cr+/fu5fPmyrSMJcddu7nLM1tvBTwohhNV4Od9UIMgQg+Xs2rWLbt26ARAZGcm+fftsnEiIu+fl8l+BkJEjBYIQFZn3Te3dYAc9CHY5C2r8+PGsXr2a2NhYDhw4QNOmTU33nTx5kuHDh3P16lU8PT2ZP38+ERERpKSk0LBhQwDc3d1JSUkpcts5OTnk5OSYrqelGde81+v16PX2c5rNW93IJhnvTnnL6HLTkY4Z2Xo8LZS7NK+/vLQde/2MJZd5KmMunQp0/3Yi5OWbtw9r5LHLAmHgwIFMmjSJjh07FrpvzJgxjB49mhEjRrBs2TJGjBjBnj178PLyMn1hpaenU6tWrSK3/e677zJt2rRCt2/cuBEXFxfLvhArWLduna0j3JFktIx169aRmFwVnI3Xd+7Zj8uhWItsOysry+znlLe2Y6+fseQyT2XLpVX+LUDyDKxZs6bEzytNm74TlaLYwXJNxQgNDWXlypWmHoTExETCw8NJTk5Gq9WiKAoBAQFs27aNa9eusWTJEmbNmsWgQYP49NNPqVq1aqFtFvUrKCgoiIsXL+Lr61tWL81ser2edevW0bNnT3Q6+1xARzJaxs0ZvzuYyUNxLfBUp7Kv/h80rtfEIvtIS0vDz8+P1NRUPDw8SvSc8tJ27PUzllzmqay55nz/Lc/q/seR3CbUeeSPEj8vKSmJgIAAs9r0ndhlD0Jx4uPjCQgIQKs1xlapVAQHBxMXF0e3bt2YP38+nTp1onfv3kUWBwCOjo44OjoSHR1NdHQ0BoMBAJ1OZ1f/CItTHnJKRsvQ6XR4OP/XRHMMKotlLs12ylvbkVzmkVzmsVYuJ51xamA+Cqi1JT7U2RpZylWBcCfR0dElfmxUVBRRUVGkpaXh6elpxVRClJ6r43+zmrNy7WOSorQdIazHUauCPON/Z+Xm4+lsu0Ody9VRDDe6M/PyjO+eoijExcURHBxs9rbkjHSiPHDU/vfrwR4WTgFpO0JYk5P2vz/LmTb+UVCuCgR/f3+aN2/OwoULAVi+fDmBgYGEh4ebvS05I50oD3TqmwsE+5guJG1HCOtx1P3X5qVAKMKYMWMIDAwkISGB3r17FygAYmJiiImJoU6dOkyfPp158+aVah/yK0iUBzePP9pLgSBtRwjrcbypB8HWw4p2OQchJiam2Pvq1q3Lzp0773ofMo4qygPtzQWCHaysBtJ2hLCmAkMMOQYbJrHTHoSyIL+CRHngoJE5CEJUJg43/Wy3dQ9CpS0QZBxVlAdaOxxikLYjhPXoNP/9WdbbuNew0hYIQpQHuptaqL0MMQghrOemEQZypUCwDekmFeWBTFIUonLR3nTkUm6eFAg2Id2kojy4eYhBb9v5SibSdoSwnpvbvPQgCCGKpVHdtMyq/Z42RQhhIQWGGKQHQQghhBAA2psmKUoPgo3IOKoob+yl/0DajhDWI5MU7YCMo4pyoWQncitT0naEsJ6CkxRlHQQhhBBCUPDIJelBEEKUiL0MMQghrEcOc7QDMo4qygM7HGGQtiOEFWk1//13jhQItiHjqEKUjrQdIaxHfdOhzQYbL45WaQsEIYQQwt7cvPSJrRdPlQJBCDtmj0MMQgjrKdCDYOPF0aRAEEIIIezETXMUMdj4FO9SIAhRTshRDEJUfDf/Uc6XHgTbkJnYojxQ2eEYg7QdIaxHddNfZelBsBGZiS1E6UjbEcJ6bp6DID0IQogSkZM5ClG5yFEMQgghhChE1kEQQgghRCE2PhWDFAhClBcywiBE5ZIvPQhCiOKo7PEwBiFEmZCjGCwkOTmZFi1a4ObmZusoQgghxF2ToxgsxN3dnXXr1tG2bVtbRxHCOuQwBiEqFTmKwUJ0Oh0+Pj62jiGERckAgxCVl62PYtDacufjx49n9erVxMbGcuDAAZo2bWq67+TJkwwfPpyrV6/i6enJ/PnziYiIuOt95uTkkJOTY7qempoKGIco7JlerycrK4ukpCR0Op2t4xRJMlrGzRnTclU4Z+WjUitk5WWQlJRkkX2kp6cDoJjRK1Fe2o69fsaSyzyVNZfqWjraLIUMfR7ZWWklbvM32qE5bfqOFBvavHmzEh8fr4SEhCgHDhwocF/Xrl2VefPmKYqiKD/88IPSsmVLRVEU5ciRI0pkZGSBy7vvvmt6Xvfu3W+7zylTpigYJ4TLRS6V/hIfH1/i9iptRy5ysf/L6dOnS9ym70SlKLYf2AwNDWXlypWmHoTExETCw8NJTk5Gq9WiKAoBAQFs27aN8PDw226rR48erF+/vtj7b/0VlJKSQkhICHFxcXh6elrk9QC0atXKokvRpqWlERQURHx8PB4eHhbZpmS0TEawbM6yyKgoCunp6VSvXh21umQjjZZqO6V5r8x5TmnfP8ll3VzmPl5ymff41NRUgoODuXbtGl5eXiV+3u3YdIihOPHx8QQEBKDVGuOpVCqCg4OJi4u7bYHQo0cPDhw4QI8ePfjwww9p2LBhocc4Ojri6OhY6HZPT0+LfhlrNBqLbu8GDw8Pi21XMtr3523tjOYWxJZqO6V5r0rzHHPfP8ll3VylbSOSyzwlLfhLwi4LhNK6Xc9BWYuKirJ1hDuSjJZTHnLaS8bS5CiL7JLLuvsoq39/kstyKtwQg7nS0tLw9PQkNTXVKr9ULaU85JSMllEeMoL95pRc5pFc5qlMuezyMEd/f3+aN2/OwoULAVi+fDmBgYEWLw7A2G06ZcqUIrtO7Ul5yCkZLaM8ZAT7zSm5zCO5zFOZctm0B2HMmDH88ssvXLp0CV9fX9zd3Tl16hQAx48fZ8SIESQlJeHh4cG8efNo1KiRraIKIYQQlYpdDDEIIYQQwr7Y5RCDEEIIIWxLCgQhhBBCFFJpCoSTJ0/Svn176tSpQ6tWrThy5EiRj/vqq6+oXbs2tWrV4qmnnkKv19tVxnPnztGlSxc8PT0LLE1dlkqSc8OGDbRu3ZoGDRoQERHBpEmTyM8vu3OXliTjzp07adq0KU2bNiUiIoIxY8YUWAjIHjLeoCgK3bp1s9gCKEIIcUcWW5PRzhW3dPPNzpw5owQEBCgXL15U8vPzlQEDBiiffvqpXWVMSkpStm7dqvz8889KkyZNyizbzUqSc//+/aYlP69fv6506NDB9Bx7yZiZmank5uYqiqIoBoNBuf/++5VZs2bZVcYbZs6cqTz55JOKp6dn2YQTQlR6laJAuHz5suLu7q7o9XpFURQlPz9fqVq1qnLy5MkCj5sxY4YyZswY0/VffvlF6dChg11lvGHjxo02KRDMzXlDVFSUMmXKlDJIWLqM169fV3r37q3Mnj3b7jIePnxY6dSpk3Lq1KlyUSAkJSUpzZs3V1xdXU23nT17VqlWrZoSGRmpDB061G5ybd26VWnTpo3Srl075YMPPrBJrpvl5uYqgwYNUiIjI5UZM2bYOk4Bc+bMUbp166ZERkaaCmtb27Vrl9K+fXulXbt2yuTJk20dR3nmmWeUjh07Km+//bato5jczXtUKYYYbrd0883i4uIICQkxXQ8NDS30GFtntLXS5Lx06RLLli3jnnvusbuM586do0mTJvj5+eHp6cm4cePsKqNer+epp54iJiYGjUZTJtnulru7O+vWraNt27YFbu/fvz+bNm1iwYIFdpMrLCyMLVu2sGPHDn7++WeysrJsku2GH3/8kfbt27Np0yb27dvHlStXbJrnhri4OA4dOsQff/zBpk2b7Obsis2aNWP79u3s2LGDnTt3kpaWZrMse/fuRavVsnXrVvbv38/ly5dtluVmd/MeVYoCQdhOWloaAwYMYNKkSbRs2dLWcQoJDQ3lr7/+4tKlS+Tk5PDjjz/aOlIB06ZN48EHH6R+/fq2jlJiOp0OHx+fQrf//vvvdOrUiUWLFtkgVdG5qlevjoODA2Bc+96S69iXxtmzZ03rvdSvX9+iJyq7G2vXriUzM5Nu3boxdepUW8cxuVGoGAwGqlevjouLi82y7Nq1i27dugEQGRnJvn37bJblZnfzHlWKAiEoKIiLFy+Sl5cHGCd8xcXFERwcXOBxwcHBxMbGmq6fO3eu0GNsndHWzMmZnp5Onz59uO+++3j++eftMuMNbm5uPProo2X2x6ukGTdv3swnn3xCaGgoHTt2JC0tjdDQUIv/shw/fjyhoaGoVCoOHjxY4D5zJlMWJyAggOPHj7N27VpiYmJKfI57a+e6Yd26ddSqVQsnJyezn2vJjHXr1mXz5s0oisLWrVtJTU01O481ciUmJqJSqdiwYQNnz54ttD1b5QJYvHgx9evXx8vLy9Qjd7dKkzElJcW0xLG7uzspKSkWyXK3ueAu3iOLD3jYqcjIyAITwlq0aFHoMadPny40SfGTTz6xq4w32GoOgqKULGd6errSvn17Zdq0aWWczqgkGU+ePGkaS83JyVEefvhh5dVXX7WrjDc7e/as1eYgbN68WYmPj1dCQkKUAwcOFLivuMmUR44cUSIjIwtc3n33XdPzunfvXuS+Xn31VWXXrl12kys+Pl7p2rWrkp6eXqJM1syYl5enjB07Vunevbvy6KOPKuvXry9VJkvn+vTTT5X58+criqIoX375pfLdd9/ZRa4bDAaD8uCDDyp///13qXPdbcZPP/1UWb16taIoivLxxx8rv/zyi0Wy3G2uG0rzHlWaAuHYsWNK27Ztldq1aystWrQwvUlPPPGEsmrVKtPjvvjiCyUsLEwJCwtTRo0aVaaTcUqSMTMzU6lRo4bi5+en6HQ6pUaNGsrLL79cZhlLmvOtt95StFqt0qRJE9PlrbfesquMMTExSkREhNK4cWOlQYMGyrPPPqtcv37drjLezJoFwg23fvGUdlKqohT8Q3zjj29+fr7Sq1cv5cKFC3aRKzs7W+nevbty7Ngxs/JYO2N+fr4ydOhQJTMz0y5y7du3T5k4caKiKIoyYcIE5c8//7SLXNnZ2ab/Hj58uHLixIm7zlXajLt371b+97//KYqiKAMHDlQuXbpk0SylzXU371GlKRCEEHd26xfP3r17lTp16hR4TKtWrZQ//vjjttvp3r274uPjo3Tv3l05dOiQ8vvvvystWrRQ2rZtW6qjBayV6+uvv1YCAgJMv0oTEhLMzmbJjLGxsUpkZKTStWvXIotEW+VSFGNhEBkZqTz11FN2k+v7779XIiMjlU6dOlnlKAZzM44bN07p2LGj8uabb1o8S2lz3c17ZJkBGyGEuMn69esLXG/YsCG9evWyUZr/FJVr5MiRNkpTWHBwMJs2bbJ1jCJ9+OGHto5QyKBBgxg0aJCtY5hER0fbOkIhd/MeVYpJikKI0rHXybP2mutm9ppRcpWevWa0Vi4pEIQQxfL396d58+YsXLgQgOXLlxMYGEh4eLjkugN7zSi5Ss9eM1ot110NhAghKoTRo0crNWrUUDQajeLv76/UqlXLdF9xkykrc67ykFFyVbyMZZ1LpSiKYoECRgghhBAViAwxCCGEEKIQKRCEEEIIUYgUCEIIIYQoRAoEIYQQQhQiBYIQQgghCpECQQghhBCFSIEghBBCiEKkQBBCCCFEIVIgCCGEEKIQKRDs3NatWwkMDDRd79Kli12eVa2y+uKLLwgMDMTb25tnn33WdLKU4kyfPp1JkyYVe//bb7/N5MmTLR1TCCHMJgVCGTh+/DgDBgzAz88PDw8P6tWrx3vvvVei53bq1ImEhAQrJyxo06ZNeHl5lek+75ZKpeLgwYMW297UqVO5//77b/uYnTt38tJLL/H7779z5swZduzYwcyZM4t9fGpqKrNmzbptgTBhwgS+/PJLLl26VNroQghhEVIglIH+/fvTpEkT4uLiuHbtGsuXLycsLMzWsaxKr9fbOoLVzZkzh2HDhhEREYG3tzevv/46c+bMobjTm3z77bd07twZPz+/Yrfp5uZG3759+eqrr6wVu8JbuXIloaGhVt1HaGgoK1euBCAiIoKff/7Zqvu7W9Z+T8x5D+zh/fruu+94+OGHbbLvWz+LXr16sX79eptkuaO7Pt2TuK0rV64ogBIXF1fsYy5duqQMGjRI8fPzU4KCgpRXX31V0ev1iqIoysaNGxVPT0/TYyMjI5XZs2cXuZ0+ffoon3/+uaIoipKSkqJoNBrlpZdeUhRFUfLz8xU/Pz9l7969iqIoyosvvqgEBwcrbm5uSv369ZXvv/9eURRFuXr1quLk5KQAiqurq+Lq6qps2bJFURRFWbdundKqVSvF09NTadCggbJq1SrTvocPH66MGjVKGTRokOLu7q58/PHHRWZcvHix0rhxY8Xd3V0JDg5W5s2bZ8r3wQcfKGFhYYq3t7fSu3dv5fTp06bnhYSEKO+9957Spk0bxc3NTencubPpPW3VqpUCKM7Ozoqrq6vy9ttvK4qiKKdOnVLuuecexc/PTwkODlbefPNNxWAwKIqiKPPmzVOaNGmivPHGG0qVKlUUf39/0/u6YsUKRafTKRqNxvQeFKV27drKDz/8YLp+4cIFBVAuXLhQ5ON79uypREdHF3nfzRYsWKC0bt36jo+rKIr7Nw0oBw4cMHt7K1asUEJCQu461+2EhIQoK1assOo+LMmS74klX/vdbis3N1dxdHRU/vjjD9Ntq1evVqpUqaJ88sknRT7HYDAoNWvWtOpZGLOzsxVPT0/lypUrhe679bPYtGmT0rRpU6tluRvSg2Blvr6+1K1bl5EjR/L9998TGxtb6DGDBw9Gp9Nx9uxZtm7dysqVK5kxY4bZ++ratSsbN24EjMMEoaGhput///03BoOBZs2aAdCkSRP27NlDSkoKr7/+OkOHDuXs2bP4+vry66+/4unpSUZGBhkZGXTq1Im///6bQYMGMX36dJKTk4mJiWHo0KEcP37ctP8lS5bwxBNPkJKSwhNPPFEo308//cQzzzzD7NmzSUlJYc+ePTRp0gQw/rqeNWsWK1eu5MKFC0RERDBgwIACY/oLFy5kyZIlXLlyBVdXV/7v//4PgN27dwOwY8cOMjIyePXVV8nKyqJ79+50796d8+fPs3XrVr777jvmzZtn2t6RI0dwcXHh/PnzLF26lBdffJHTp09z//338+qrr3LPPfeY3oOiXLx4keDgYNP1gIAAHB0duXDhQpGPP3jwIPXq1bvDpwgNGjSw6HCJ+I+iKBgMBlvHEBZy5MgRcnNzadGiBYqiMHXqVMaOHcuKFSt45plninzOmjVr8PHxoVGjRlbLtXHjRho1anTb3sIbOnfuTEpKCtu3b7dantKSAsHKVCoVmzZtokmTJkybNo2wsDAaNGjAunXrADh//jwbNmxg1qxZuLm5ERISwuTJk5k/f77Z++ratSubNm0CYMOGDUyYMIHTp0+TlpbGhg0biIyMRK02fuRDhgzB398fjUbDo48+Sr169dixY0ex246JiWHEiBF069YNtVpNx44dueeee/j+++9Nj+nVqxe9e/dGrVbj4uJSaBufffYZEyZMMG3D39/fVLB8++23jB8/nkaNGuHk5MQ777xDfHy86Y8/wLhx46hZsyZOTk4MGTKEffv2FZv3l19+wdvbm+eeew4HBweCg4OZMGECixcvNj3Gz8+PiRMnotPp6NKlC6GhoWb/Ye7RowdeXl6mS25uLiqVqsjHXrt2DQ8Pjztu08PDg9zcXLKysszKUlF99NFHdOnSpcBt3333HQ0aNAAgISGBXr164eHhQYsWLTh69GiBx4aGhvLuu+/Stm1bXFxcOHr0KLNmzaJ27dq4u7tTq1YtPv300wLPuXVOy4cfflgow83bvzHcAJCWlsYzzzxDSEgIHh4etGrVivj4+CKfe7vHXr58mYcffpgqVaoQHBzM5MmTTQVzaGgoM2bMoG3btri7uxMZGVlgH3d6T+70+orLNWjQIOLi4njsscdwc3Nj7NixBd6D2bNn061btwL7Wrp0aYHC+MZji9rWnT7rW+3fv5/w8HAURWHAgAGsXbuWvXv30qFDhyIfD7B69eoCGW/3Gdz4t9OqVStcXV3p27cvycnJjBs3Di8vL2rXrl3k9+ZPP/3EvffeC5Tss+jWrRurV68uNrOtSIFQBqpVq8bMmTM5cuQIV65coW/fvjzwwAMkJyeTkJCAk5MTVatWNT0+LCysVBMTmzdvTnZ2NkeOHGHDhg306NGDjh07snXrVjZs2FCgUcyePZuIiAg8PT3x8vLi8OHDXL16tdhtnzt3jjlz5hT4Y7hq1aoCv5Zv/jVdlNjYWGrXrl3kfQkJCQXG5RwdHalevXqB96FatWqm/3Z1dSU9Pf22eQ8fPlwg78SJEwtM/rv5PS/JNm8VEBDAb7/9RkpKCikpKVy7dg2tVktAQECRj/f29iYtLc10fezYsbi5uRX4ogXjF5aDg0ORRVZl9Pjjj7Nr1y7Onj1rum3evHmMHDkSMPbABQQEcOnSJRYtWsTcuXMLbWP+/Pl88803ZGRkULduXUJCQtiwYQNpaWl8+eWXvPjiixb7BTdixAhOnTrFzp07SUlJ4YsvvsDZ2dnsx96pZ7G4HrWSvieleQ0//PADwcHBLFmyhIyMDObMmVPgeYMHD2bbtm0FipVvv/2WoUOHFtpHUdu602d9q/379+Pm5karVq0ICgpi06ZNxba/G27tybvT57V06VJ+/PFHLly4QHx8PG3btqVHjx4kJSUxePDgAm33hpsLhJJ8FvbaaygFQhnz8fFh6tSpZGZmcvbsWQIDA8nOzuby5cumx5w7d67AoY0lpdFo6Ny5M0uXLiU5OZn69evTrVs31q1bx5YtW+jatSsA27ZtY+rUqSxYsIBr166RkpJCw4YNTZPrbvQy3CwoKIgJEyaY/himpKSQkZHB559/bnpMUc+7WUhICKdOnSryvsDAQM6dO2e6npuby4ULF0r8Ptz6qz0oKIgWLVoUyJuWlsaRI0dKtL07vRaAdu3aceLECdP1M2fOUKNGjQKFzM2aNm3KsWPHTNfnzJljGsK4+Yv26NGjNG3atEQ5K4pXXnmlQDF381E0vr6+3HvvvXzzzTeAsddt8+bNDB06lPj4eLZu3cr777+Pi4sL9erVK/IL++mnn6Zu3bpoNBocHBx46KGHCAoKQqVS0bVrV3r37m3qfbsbly9fZsWKFXzxxRdUr14dtVpNs2bNiuxqvt1jS9KzWFyPWknfE0u8hltVrVqVHj16sGjRIgASExNZt25dkQVCUW73WRdl//79xMbGotFo+Pjjj3FwcLjjPm7uySvJa3366acJCgrC09OTfv364evry4MPPohGo+GRRx7h8OHD5Obmmh5/8OBBnJycqFu3bok/Cw8PD65du1ai96gsSYFgZdeuXeO1117j2LFjGAwGsrKymDVrFj4+PtSrV48aNWrQtWtXXnjhBTIzM4mLi+Ptt99m+PDhpdpf165d+eijj4iMjASgW7duzJs3DycnJxo2bAgYf6FqNBqqVKlCfn4+X3/9NYcPHzZto2rVqqSnp5OYmGi6bcyYMcybN4+NGzdiMBjIyclh586d/PPPPyXONmbMGD766CM2b95Mfn4+iYmJHDhwADD+Svz00085evQoOTk5vPbaa9SoUYPWrVuXaNtVq1bl9OnTpuv33HMPly9f5rPPPiM7OxuDwcDx48dL/EegatWqxMbG3nZdg7FjxzJjxgzOnj1LWloakyZNYuzYscUOMQwYMMA0J+R2NmzYwD333FOinBXFu+++W6CYS0lJKXD/qFGjWLBgAYqisGDBAnr16kW1atW4cOECTk5O+Pv7mx4bEhJSaPu39m4tWrSI5s2b4+Pjg5eXF2vWrLltD1pJxcbG4ujoeMfetDs9tiQ9i8X1qJX0PbHEayjKsGHD+PbbbwHjvKT27dubta3iPutbGQwG/vrrL+bPn4+bmxtjxowp0fZv7skryWu9+TNwcXEpdF1RlALDgTf3HpT0s0hLS8Pb27tE+cuSFAhW5uDgwPnz5+nXrx+enp4EBwezfft2fv31V1xdXQFYvHgx169fJyQkhA4dOtC/f//bHit/O127diUtLc00nNCwYUOcnZ1NvQcAffr0YeDAgTRq1Ijq1atz5MiRAmN2devW5YknnqBBgwZ4eXmxbds2mjVrxpIlS3jttdeoUqUKNWrU4P/+7//Iyckpcbb777+fWbNmERUVhaenJ61ateLQoUOA8Uvl2Wef5Z577qFatWr89ddf/PTTT2i12hJt+80332T8+PF4e3szffp03NzcWL9+PX/88QehoaH4+voyePDgEq8vMGjQIDw8PKhSpUqxa0K0a9eO5557jg4dOhAcHEz16tWZOHFisdscOnQomzdvJikpqdjHZGZmsmbNGp588skS5awsevbsSV5eHps3b+abb74xdTlXr16d7OzsAsVsXFxcoeff3CMUFxfH8OHDmTFjBomJiaSkpNCvX78Ch6e6uroW+NK/ePFiiXKGhISQk5NT7JyDkj72bnoWS/Ke3O713ek13Kl37b777iMhIYF9+/YVO7xwu20V91nf6tixY2RlZdGxY0dWrVrFb7/9VqL1ZW7uyTPn8yqp1atXmwqEkv77tNteQ9sdQCFE5fPOO+8okyZNKvb+t99+W3n11VfLMJHtlfQwx9dff11p0qSJUqVKFSU3N9d0e4cOHZSRI0cqWVlZyrFjx5SwsLACh5HdeijdkSNHFLVarfz111+KwWBQfvnlF8XZ2VmZMGGC6TGdOnVSxo4dq+j1euXAgQNKQECAEhkZWeQ2b93+fffdp/Tr10+5cOGCYjAYlP379ytXr14t8rXf7rFdu3ZVHn/8cSUjI0OJjY1VGjZsqLz11ltF7vPWQ+fu9J7c6fXdLlfbtm2VDz74oMDruDXPqFGjlN69eyvOzs5KWlpasY8taluKUvxnfbMFCxYooaGhpuu7d+9WXFxc7njY5OrVq5WWLVuW6LXe+rqmTJmi3HfffabrZ8+eVQDl2rVriqIYD3X29fVV8vLyTI+502ehKIpSs2ZN0+Hk9kR6EIQoQ6+88sptf+W8+uqrvP3222WYqPwYOXIkf//9N48//jg6nc50++LFi4mPj8ff35/BgwczatSo226nQYMGTJ48mW7duuHr68vSpUtNv/hu+OSTT9i5cydeXl689NJLZg35ffPNNwQFBdGyZUu8vLwYO3Ys169fN/uxd9OzeKf35E6v73a5Xn31VT799FO8vLwYN25ckfsfNmwYv//+O/fffz/u7u7F5ixuW8V91jfbv38/zZs3N11v1aoVX331FY8//jj79+8vdp/9+vXj6tWrpmFVcz6vO/n555/p27cvGo3GdNudPoutW7fi4eFBp06dSrVPa1IpSjHLvgkhhB3JysrC39+fP//80zSfRlRM1v6slyxZwsqVK1m6dKlFtztgwACGDRvGoEGDSvyc3r1788ILL9CzZ0+LZrEEKRCEEHZPURRmzJjBzz//zNatW20dR1hRef6sZ8yYwdNPP33bXpPyRAoEIYRdMxgMeHl54efnx/Llywt0K4uKRT5r+yIFghBCCCEKkUmKQgghhChECgQhhBBCFCIFghBCCCEKkQJBCCGEEIVIgSCEEEKIQqRAEEIIIUQhUiAIIYQQohApEIQQQghRiBQIQgghhChECgQhhBBCFPL/Jxg2TV0aIDgAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, axes = plt.subplots(\n", " ncols=2, figsize=(5.2, 3.35), layout=\"tight\", sharey=True, width_ratios=[1.0, 1.2]\n", ")\n", "axes = pe.plot.curves(mg, axes=axes, linewidth=2.0, label=\"Mualem-van Genuchten\")\n", "axes = pe.plot.curves(bcf, axes=axes, linewidth=1.5, label=\"Brooks-Corey fit\")\n", "\n", "axes[0].set_xlim(0.0, 0.46)\n", "axes[0].xaxis.set_minor_locator(mplticker.MultipleLocator(0.05))\n", "axes[0].xaxis.set_major_locator(mplticker.MultipleLocator(0.1))\n", "axes[1].set_ylim(h[0], h[-1])\n", "xt1 = np.logspace(-15, 3, 19)\n", "axes[1].set_xlim(xt1[0], xt1[-1])\n", "log_locator = mplticker.LogLocator(base=10.0, subs=(1.0,), numticks=20)\n", "axes[1].xaxis.set_minor_locator(log_locator)\n", "axes[1].xaxis.set_minor_formatter(mplticker.NullFormatter())\n", "axes[0].set_ylabel(\"Pressure head |\\N{GREEK SMALL LETTER PSI}| (cm)\")\n", "axes[0].set_xlabel(\"Soil water content \\N{GREEK SMALL LETTER THETA} (-)\")\n", "axes[1].set_ylabel(\"\")\n", "axes[1].set_xlabel(\"Hydraulic conductivity $K$ (cm/d)\")\n", "\n", "handles, labels = axes[0].get_legend_handles_labels()\n", "fig.legend(\n", " handles,\n", " labels,\n", " bbox_to_anchor=(0.13, 0.92),\n", " loc=\"lower left\",\n", " ncol=2,\n", " frameon=False,\n", ")\n", "fig.tight_layout(w_pad=0.1)\n", "fig.align_xlabels()\n", "fig.savefig(fig_path / \"fit_swrc_hcf.png\", dpi=300, bbox_inches=\"tight\")" ] } ], "metadata": { "kernelspec": { "display_name": "pedon (3.13.5)", "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.13.5" } }, "nbformat": 4, "nbformat_minor": 5 }