{ "cells": [ { "cell_type": "markdown", "id": "7545f9ad", "metadata": {}, "source": [ "# Armadillo real world example - Frequentist fit\n", "\n", "In progress" ] }, { "cell_type": "code", "execution_count": 1, "id": "d175860d", "metadata": {}, "outputs": [], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import pandas as pd\n", "from pysip.regressors import Regressor\n", "from pysip.statespace import TwTi_RoRi" ] }, { "cell_type": "code", "execution_count": 2, "id": "1b96ade6", "metadata": {}, "outputs": [], "source": [ "# Load and prepare the data\n", "df = pd.read_csv(\"../../../../data/armadillo/armadillo_data_H2.csv\").set_index(\"Time\")\n", "df.drop(df.index[-1], axis=0, inplace=True)\n", "inputs = [\"T_ext\", \"P_hea\"]\n", "outputs = \"T_int\"\n", "sT = 3600.0 * 24.0\n", "df.index /= sT # Change time scale to days" ] }, { "cell_type": "code", "execution_count": 3, "id": "c4e18d00", "metadata": {}, "outputs": [], "source": [ "# Parameter settings for second order dynamic thermal model\n", "parameters = [\n", " dict(name=\"Ro\", scale=1e-2, transform=\"log\"),\n", " dict(name=\"Ri\", scale=1e-3, transform=\"log\"),\n", " dict(name=\"Cw\", scale=1e7 / sT, transform=\"log\"),\n", " dict(name=\"Ci\", scale=1e6 / sT, transform=\"log\"),\n", " dict(name=\"sigw_w\", scale=1e-3 * sT**0.5, transform=\"log\"),\n", " dict(name=\"sigw_i\", value=0.0, transform=\"fixed\"),\n", " dict(name=\"sigv\", scale=1e-2, transform=\"log\"),\n", " dict(name=\"x0_w\", loc=25.0, scale=7.0, transform=\"none\"),\n", " dict(name=\"x0_i\", value=26.7, transform=\"fixed\"),\n", " dict(name=\"sigx0_w\", value=0.1, transform=\"fixed\"),\n", " dict(name=\"sigx0_i\", value=0.1, transform=\"fixed\"),\n", "]" ] }, { "cell_type": "code", "execution_count": 4, "id": "721b1df5", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Optimization terminated successfully.\n", " Current function value: -331.057569\n", " Iterations: 25\n", " Function evaluations: 510\n", " Gradient evaluations: 34\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
θσ(θ)pvalue|g(η)||dpen(θ)|
Ro0.0175930.0009270.00.0000223.230694e-17
Ri0.0019840.0000750.00.0000492.539865e-17
Cw169.5971127.6625830.00.0000234.657315e-17
Ci18.9463500.7715800.00.0000223.731826e-17
sigw_w0.5213440.0469860.00.0000073.178810e-17
sigv0.0343250.0023330.00.0000298.487470e-18
x0_w26.5945390.1305170.00.0000110.000000e+00
\n", "
" ], "text/plain": [ " θ σ(θ) pvalue |g(η)| |dpen(θ)|\n", "Ro 0.017593 0.000927 0.0 0.000022 3.230694e-17\n", "Ri 0.001984 0.000075 0.0 0.000049 2.539865e-17\n", "Cw 169.597112 7.662583 0.0 0.000023 4.657315e-17\n", "Ci 18.946350 0.771580 0.0 0.000022 3.731826e-17\n", "sigw_w 0.521344 0.046986 0.0 0.000007 3.178810e-17\n", "sigv 0.034325 0.002333 0.0 0.000029 8.487470e-18\n", "x0_w 26.594539 0.130517 0.0 0.000011 0.000000e+00" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Instantiate the model and use the first order hold approximation\n", "model = TwTi_RoRi(parameters, hold_order=1)\n", "reg = Regressor(model, inputs=inputs, outputs=outputs)\n", "fit_summary, corr_matrix, opt_summary = reg.fit(df=df)\n", "fit_summary" ] }, { "cell_type": "code", "execution_count": 5, "id": "8cf43718", "metadata": {}, "outputs": [], "source": [ "# Predict the indoor temperature each minute\n", "dt = 60 / 3600\n", "tnew = np.arange(df.index[0], df.index[-1], dt)\n", "ds = reg.predict(df=df, tnew=tnew)\n", "ym = ds[\"y_mean\"].sel(outputs=\"T_int\")\n", "ysd = ds[\"y_std\"].sel(outputs=\"T_int\")" ] }, { "cell_type": "code", "execution_count": 6, "id": "80d9366b", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAHVCAYAAAB8NLYkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAACSAElEQVR4nOzdd3iTZfv/8XeSNunI6Er3optSyhZRnGwHiDgffRBBFEQUUVEUVARkiAMcqCgIAqKgKC5Q0SIoG8qmrBYKXXQl3WmT/P7oY39fFKGFtuk4X8eR47DJnbtnbEk+ve7rOi+F3W63I4QQQgghmj2lowsQQgghhBD1Q4KdEEIIIUQLIcFOCCGEEKKFkGAnhBBCCNFCSLATQgghhGghJNgJIYQQQrQQEuyEEEIIIVoICXZCCCGEEC2EBDshhBBCiBZCgp0QQgghRAvRZILdzJkzUSgUjBs3DoD8/HzGjh1LbGwsrq6uhIaG8vjjj2MymS54nmHDhqFQKM659e/fvxFegRBCCCGEYzk5ugCA7du388EHH5CYmFhzX0ZGBhkZGcyZM4f4+HhOnjzJqFGjyMjIYNWqVRc8X//+/Vm0aFHN1xqNpk712Gw2MjIy0Ol0KBSKur0YIYQQQoh6ZLfbKSoqIjAwEKXywmNyDg92xcXF3HfffSxYsIBp06bV3J+QkMCXX35Z83VkZCTTp0/n/vvvp6qqCienfy9do9Hg7+9/yTVlZGQQEhJyyc8XQgghhKhv6enpBAcHX/AYhwe7MWPGcPPNN9O7d+9zgt35mEwm9Hr9BUMdQFJSEr6+vnh6enLjjTcybdo0vL29//X4iooKKioqar622+1A9f9AvV5fh1cjhBBCCFG/zGYzISEh6HS6ix7r0GC3YsUKdu3axfbt2y96bG5uLlOnTuXhhx++4HH9+/fn9ttvp02bNhw/fpznn3+eAQMGsHnzZlQq1XmfM2PGDKZMmfKP+/V6vQQ7IYQQQjQJtZkeprD/NTzVyNLT0+natSs///xzzdy666+/no4dO/LWW2+dc6zZbKZPnz54eXmxZs0anJ2da/19Tpw4QWRkJL/88gu9evU67zF/H7H7Kxn/NUIohBBCCOEoZrMZg8FQq1zisFWxO3fuJCcnh86dO+Pk5ISTkxMbNmxg3rx5ODk5YbVaASgqKqJ///7odDpWr15dp1AHEBERgY+PD8eOHfvXYzQaTc3onIzSCSGEEKK5ctil2F69erFv375z7nvwwQeJi4vj2WefRaVSYTab6devHxqNhjVr1uDi4lLn73P69Gny8vIICAior9KFEEIIIZokhwU7nU5HQkLCOfe5u7vj7e1NQkICZrOZvn37UlpaytKlSzGbzZjNZgCMRmPNfLm4uDhmzJjB4MGDKS4uZsqUKQwZMgR/f3+OHz/OhAkTiIqKol+/fo3+GoUQQojmwmaz1VwtE41LpVJdtI1JbTl8Vey/2bVrF1u3bgUgKirqnMdSU1MJDw8HICUlpaZpsUqlYu/evSxevJjCwkICAwPp27cvU6dOrXMvOyGEEKI1sNvtmEwmSktLHV1Kq+bm5obBYLjs/rkOWzzRlNVlkqIQQgjRnBUWFlJaWoper0etVktj/kZmt9uxWCyYzWbc3Nzw8PD4xzF1ySVNdsROCCGEEA3LZrPVhDqtVuvoclottVoNVAc4vV5/WZdlm8xesUIIIYRoXH/NqfsrWAjH+etncLnzHCXYCSGEEK2cXH51vPr6GUiwE0IIIYRoISTYCSGEEKJFuP766xk3bpyjy3AoCXZCCCGEaHWSkpJQKBQUFhY6upR6JcFOCCGEEKKFkGAnhBBCiGanpKSEoUOHotVqCQgI4PXXXz/n8U8//ZSuXbui0+nw9/fnP//5Dzk5OQCkpaVxww03AODp6YlCoWDYsGEArF27lp49e+Lh4YG3tze33HILx48fb9TXdjmkj50QQgghatjtdiodtAuFs5tbrVeHPvPMM2zYsIFvvvkGX19fnn/+eXbt2kXHjh0BqKysZOrUqcTGxpKTk8P48eMZNmwYP/zwAyEhIXz55ZcMGTKElJQU9Ho9rq6uQHVgHD9+PImJiRQXF/Piiy8yePBgkpOT623br4YkwU4IIVqZM2eKqKqy4ePjSmWljYoKK15eLjg7qxxdmmgCKktLmeegZsWPFxejdne/6HHFxcV8/PHHLF26lF69egGwePFigoODa44ZPnx4zX9HREQwb948unXrRnFxMVqtFi8vLwB8fX3P2e1hyJAh53yvhQsXYjQaOXjw4D/2uG+Kmn70FEIIUW+Kiy0cPJjHjh1ZbNqUwaaN6Xw3dgKLb7yZM2l5ji5PiFo5fvw4FouF7t2719zn5eVFbGxszdc7d+7k1ltvJTQ0FJ1Ox3XXXQfAqVOnLnjuo0ePcu+99xIREYFer6/Zm/5iz2sqZMROCCFakaNHC9i+/QzXXBNCeVkVv898jfI/fyCEM+z96juCxj/g6BKFgzm7ufF4cbHDvnd9KCkpoV+/fvTr149ly5ZhNBo5deoU/fr1w2KxXPC5t956K2FhYSxYsIDAwEBsNhsJCQkXfV5TIcFOCCFaidzcUkaN+o5t206j0ahwslkoqQwCRuJNPqN+TKbfuKEolbILQWumUChqdTnUkSIjI3F2dmbr1q2EhoYCUFBQwJEjR7juuus4fPgweXl5zJw5k5CQEAB27NhxzjnOt4VXXl4eKSkpLFiwgGuuuQaATZs2NcZLqjdyKVYIIVqBoiILr776B9u2nQagosJKSaUKHUVolFby8OK7beVUVFQ5uFIhLk6r1TJixAieeeYZfv31V/bv38+wYcNqFjeEhoaiVqt5++23OXHiBGvWrGHq1KnnnCMsLAyFQsF3333H2bNnKS4uxtPTE29vbz788EOOHTvGr7/+yvjx4x3xEi+ZBDshhGjhysur+PrL/bz79p8ADPFJ5knmM0a1hGXTopj6bEcAsopUlJVUOLBSIWrvtdde45prruHWW2+ld+/e9OzZky5dugBgNBr55JNPWLlyJfHx8cycOZM5c+ac8/ygoCCmTJnCc889h5+fH4899hhKpZIVK1awc+dOEhISePLJJ3nttdcc8fIumcJut9sdXURTYzabMRgMmEwm9Hq9o8sRQohLZrfbSfrpCHcOXESexZUoTjCST3HW6mg/ZxEu7a7g6JGzDB+xBieqSF7Xn3Z9r3Z02aKRVFZWcvbsWYxGI87Ozo4up1W70M+iLrlE5tgJIUQLduJQBsMHzyfP4om3ooCnB7uhMb6Mtuv1mL38sZgtRER6o1LYqLI7seenPyXYCdGMSbATQogWyG63k5FRzLO3TSWtLAA3ypk9pQce7ToRFqYnJERHWVkVBoMGg0FDoCek58P+Pw44unQhxGWQYCeEEC1MZaWVI0cK+GXxOtYc9QHgseExeCZ0JjhYS0KCDxrNuW//ERGepOebOHYoyxElCyHqiSyeEEKIFqSioop9+3L5c9Mp3nh7G5U4E+9TQs/beuLl5UJ8vPc/Qh1AQpdwALJMUJyV3chVCyHqiwQ7IYRoIex2O0eOFLBu7VGeefIbTpV74UI5T790M87OTiQk+ODurj7vc9t1qN6KKRdvTm38ozHLFkLUI7kUK4QQLURGRjHffPAj099Lodyuxp9sJj3oj0dQIFFRHnh7u/7rc2NjvQE4izeFJ9Mbq2QhRD2TYCeEEM1ccbGF7OwS1iz9k+nvHqIcFyJVp5n+eAS63neh02kICdFd8BzR0dUboufjgTmnoDHKFkI0AAl2QgjRTFVWWjl1ysyJEyaSd53m5Vf+rA51bvl8+MVY7C56ysqqiIryOO+8uv8rKEiPWmXHYlVxMr2wcV6AEKLeSbATQohm6vDhfHbsyOT77w6z5puDVNrUhCszmL9oGJVOWmwWG4mJvvj7X3zfT6VSQZAXpJ6FtDNljVC9EKIhSLATQohmKDe3lOXL9/He239iLrECCqI4wcwXrwbPAOxWOx07GgkI0Nb6nCG+zqSereJ0tqXhChfNhsViparK1mjfz8lJiVqtarTv11JJsBNCiGbGarXx2vRfmf3WTgACyOJWzSYGjL0H1yv7YLFY6djRt06hDiDAzxUOFJFnbrwPc9E0WSxWtm3LpLi4stG+p1brzBVXBNQ63BUVFTF58mRWr15NTk4OnTp1Yu7cuXTr1q3mmGHDhrF48eJzntevXz/Wrl0LQEVFBQ899BDffPMN/v7+vPfee/Tu3bvm2Ndee41Tp07x9ttvX7Qes9nMrFmz+PLLL0lLS8PDw4OEhAQeffRRBg8ejEKh4Prrr6djx4689dZbtXqNl0KCnRBCNDNbf9jOu3M3Aa70ZAsjb/cj9MGlFFa54eSkJCHBh6CgCy+WOJ/gIB1QRH6pEpvNjlKpqPfaRfNQVWWjuLgStVqJRtPwo2gVFVaKiyupqrLVOtg99NBD7N+/n08//ZTAwECWLl1K7969OXjwIEFBQTXH9e/fn0WLFtV8rdFoav77ww8/ZOfOnWzevJkff/yR//znP2RnZ6NQKEhNTWXBggXs2LHjorUUFhbSs2dPTCYT06ZNo1u3bjg5ObFhwwYmTJjAjTfeiIeHR+3/h1wGCXZCCNGM5OaWMm3MQkrsgfg5mZj07miqgtpRZLMTHq4jKsoTrfb8veouJjjUA8igoNyZykrrRRdciJZPo1Hh4tI4vwcWS+1HisvKyvjyyy/55ptvuPbaawF4+eWX+fbbb5k/fz7Tpk2rOVaj0eDv73/e8xw6dIiBAwfSrl07IiIieOaZZ8jNzcVoNDJ69GhmzZqFXq+/aD3PP/88aWlpHDlyhMDAwJr7Y2JiuPfee3Fxcan1a7tc0qBYCCGaiYqKKn5YksQv6UYAHh3TnVJjHHq9hiuuCKBDB99LDnUAIeHVvexMVRoqK+VyrGi6qqqqsFqt/whMrq6ubNq06Zz7kpKS8PX1JTY2ltGjR5OXl1fzWIcOHdi0aRNlZWWsW7eOgIAAfHx8WLZsGS4uLgwePPiitdhsNlasWMF99913Tqj7i1arxcmp8f5Ikj/HhBCimTh+vJCpU36lEh2xHmYSru1IeLiBuDgvXF2dL/v8YVG+AJhs7lRUVF1WSBSiIel0Onr06MHUqVNp27Ytfn5+fPbZZ2zevJmoqKia4/r378/tt99OmzZtOH78OM8//zwDBgxg8+bNqFQqhg8fzt69e4mPj8fHx4cvvviCgoICXnzxRZKSkpg0aRIrVqwgMjKShQsXnnOJ9y+5ubkUFBQQFxfXmP8L/pUEOyGEaAby88uYNflbjpl1aKjgmQk9cXfXEBXlWS+hDqBNXPWHVjkunD2Tj7e3W72cV4iG8OmnnzJ8+HCCgoJQqVR07tyZe++9l507d9Ycc88999T8d/v27UlMTCQyMpKkpCR69eqFs7Mz77777jnnffDBB3n88cfZvXs3X3/9NXv27GH27Nk8/vjjfPnll/+ow263N9yLvARyKVYIIZq48vIqlry3gaVfpQFwb8xpvOISCAvTodPV36ial58HGioASEs5U2/nFaIhREZGsmHDBoqLi0lPT2fbtm1UVlYSERHxr8+JiIjAx8eHY8eOnffx3377jQMHDvDYY4+RlJTETTfdhLu7O3fddRdJSUnnfY7RaMTDw4PDhw/Xx8u6bBLshBCiCauqsvHDsj+Y9GISNpR0cT3OnS+ORKt1vug2YXWlUCgwqEoBOHk0p17PLURDcXd3JyAggIKCAtatW8egQYP+9djTp0+Tl5dHQEDAPx4rLy9nzJgxfPDBB6hUKqxWK5WV1e1eKisrsVqt5z2nUqnknnvuYdmyZWRkZPzj8eLiYqqqqi7x1dWdBDshhGjC9uw8w+OjV1NidyXEOY9XFz9KuasfkZEeuLvX/xw4T+fqEbvU47n1fm4h6tO6detYu3Ytqamp/Pzzz9xwww3ExcXx4IMPAtWB6plnnmHLli2kpaWxfv16Bg0aRFRUFP369fvH+aZOncpNN91Ep06dALj66qv56quv2Lt3L++88w5XX331v9Yyffp0QkJC6N69O0uWLOHgwYMcPXqUhQsX0qlTJ4qLixvmf8J5yBw7IYRoonJzSxl953ucqfREqyjlrY/uoVzjQ5CfG6GhF2/BcCm83aqgHLIyzA1yftG8VFScf5SqKXwfk8nExIkTOX36NF5eXgwZMoTp06fj7Fw951SlUrF3714WL15MYWEhgYGB9O3bl6lTp57Tyw5g//79fPHFFyQnJ9fcd8cdd5CUlMQ111xDbGwsy5cv/9davLy82LJlCzNnzmTatGmcPHkST09P2rdvz2uvvYbBYKjz67tUCntTm/XXBJjNZgwGAyaTqVb9a4QQor5VVdkY/+AS3l56EhVWZjziS/ytA1CrVXTt6o+nZ8P0xbqtzUi+SQtmQBc1329/DoVCmhS3ZJWVlZw9exaj0VgTiKB57DzR0vzbzwLqlktkxE4IIZqgXdtOsWDZUUDNXTEZxAwYiru7MwkJPg0W6gD8PFWQBrkFFVRW1n4XANGyqNUqrrgiQPaKbYYk2AkhRBNTXl7Fk/+ZT7ndjRBVDndPGYVW60ynTn4YDJqLn+Ay+PtUh8Z8sxWLxSoftK2YWq2Sn38zJIsnhBCiCbHb7Syb8zV/nnRDgZ2xI9uh0noQH+/d4KEOIMDfFYCCEhp1tEYIUT8k2AkhRBNRVGRh17ZTzJ76MwA9gwqJ7NObqChPAgK0jVJDUFD1/J3CclWjTZwXQtQfCXZCCNEEVFZa2b//LCuef4cjlkBUWHnw+bvw8nIhPLzxVtT5B3ugwI7NriAzs/FaNAjHknWUjldfPwMJdkII0QScOVPMb2t2sei36lWI/bu54RloJCrKExeXxpsOrfPxwp0SALKyJNi1dCpV9Rw6i8Xi4ErEXz+Dv34ml6rJBLuZM2eiUCgYN25czX1/dYH29vZGq9UyZMgQsrOzL3geu93Oiy++SEBAAK6urvTu3ZujR482cPVCCHHpSkoszJ27hUkzd5Jn98TgVM7dY2/Gz8+NgAD3Rq1F4+GBjupAd+aM9LJr6ZRKJW5ubpjNZoqLi7FYLFRWVsqtEW8Wi4Xi4mLMZjNubm4olZcXzZrEqtjt27fzwQcfkJiYeM79Tz75JN9//z0rV67EYDDw2GOPcfvtt/PHH3/867lmz57NvHnzWLx4MW3atGHy5Mn069ePgwcP4uLScC0ChBDiUtjtdlZ/eYi3527GjoJO7GXSW8NR6asvwapUjfv3t4uHB9r/jdilp0uwaw3+ap5rNsvP25Hc3NzqpZGxw4NdcXEx9913HwsWLGDatGk195tMJj7++GOWL1/OjTfeCMCiRYto27YtW7Zs4corr/zHuex2O2+99RaTJk2q2StuyZIl+Pn58fXXX3PPPfc0zosSQoha2r87nSce+gyr3ZVEDvDCmLZowuPw8nLB19et0evReHrWjNjl5JQ0+vcXjU+hUODh4YFer//X/VBFw1KpVJc9UvcXhwe7MWPGcPPNN9O7d+9zgt3OnTuprKykd+/eNffFxcURGhrK5s2bzxvsUlNTycrKOuc5BoOB7t27s3nz5n8NdhUVFVRUVNR8LX+1CCEaQ3l5FY/f9S75lW54KwqY/FJvPHr2o6Sk0iGjdXDupdjcsxLsWhOlUllv4UI4jkOD3YoVK9i1axfbt2//x2NZWVmo1Wo8PDzOud/Pz4+srKzznu+v+/38/Gr9HIAZM2YwZcqUOlYvhBCXZ8WCDSQdrx6Ve250HK5de1FZaaNDByN+fo07t+4vGoPh/we7HDNWq80hAVMIcWkc9q81PT2dJ554gmXLljl87tvEiRMxmUw1t/T0dIfWI4Ro+QoLy3lh4k8AXO1zhui+vXBzc6ZrVz/CwgwolY7Zo9VJo8HgXL06Ly+3RJoUC9HMOCzY7dy5k5ycHDp37oyTkxNOTk5s2LCBefPm4eTkhJ+fHxaLhcLCwnOel52djb+//3nP+df9f185e6HnAGg0GvR6/Tk3IYRoKFarjZnPfE5GiRtulPLQUzdhMGjo1MkXH5/Gn1f3d9766o+GgoIyKisl2AnRnDgs2PXq1Yt9+/aRnJxcc+vatSv33XdfzX87Ozuzfv36muekpKRw6tQpevTocd5ztmnTBn9//3OeYzab2bp1678+RwghGlNZWSU7dmTx4eJDANzatgiv+ATi4rzx8GgaK/d9DNV9tApNFgl2QjQzDptjp9PpSEhIOOc+d3d3vL29a+4fMWIE48ePx8vLC71ez9ixY+nRo8c5Cyfi4uKYMWMGgwcPrumDN23aNKKjo2vanQQGBnLbbbc15ssTQoh/sNvtHD6cz9svfUlBpSsGzNz62J2Ehurw83P8SN1fjN4aOAFFpTbKyirx9GwagVMIcXEOXxV7IW+++SZKpZIhQ4ZQUVFBv379eO+99845JiUlBZPJVPP1hAkTKCkp4eGHH6awsJCePXuydu1ah8/jE0KInJxSFi7czWfrcgElgzta8IkMJyLCA4XCMXPqzsfbxx0lNmwoycwsJjBQ5+iShBC1pLDLBnH/YDabMRgMmEwmmW8nhKgXlZVWXnllI9OmbQCgi3I/z7w/lg7XtCMuztvB1Z3r+/vv5z/L/DGj49tv7+WWW2IcXZIQrVpdcomsYRdCiEawP/kM8177DYBr2MzYwQY8QkMJCmp6o2HnbitW5OBqhBB1IcFOCCEaWF5uCS8OfB5zhRIv8rn/Khvutz1KWJgOnU7t6PL+wRAejvZ/wU62FROieWnSc+yEEKK5M5srWP3yu/ycFQrAiLsj8LnrUfz83AkNbZpTPfyvuAIdmwDIzS11cDVCiLqQYCeEEA3AbreTnl7Enz8fYNJ7p6nAmyg/6HH3Tfj5uZGYaMTNzdnRZZ6XX5cu6BSlYIfsU2cdXY4Qog7kUqwQQtQzq9VGSko+339/lAnjviHb7o2nUykz5g5BrVYSG+uFVtv0LsH+Re3ujn9A9dy/nJNZVFRUObgiIURtyYidEELUI5vNTkpKPosX72Hem79TXumKOyW88Up3nF00BAVp8fJq+u2XwuOCIaN6v9jS0io0Gvm4EKI5kBE7IYSoR2lpJj78cCevvbaR8kolbTjJWw/raHNVN5ydVYSFGZpUz7p/E901FoA8UxUlJZUOrkYIUVsS7IQQop5kZBSzYMEu5s3djN2uoBu7mP14KIGD/kNRkYXISI9ms4tDz0HV2zDmVbpz6mS+g6sRQtSWBDshhKgH+fllfP9tCm++loQdBd3ZwTNjO6K6/j+4uTnTrVsAsbFeji6z1iK7J+KjKABgw9fbHFyNEKK2ZNKEEEJcptLSSjb/nsrkcV9QYXUjklQmTLweut5KUJCW+HifJr1Y4nyUKhVR3hZycyH5j8OUl1fh4iIfGUI0dTJiJ4QQl8FqtbF/RxqTHphHdrkbBsxMmXw1dL2VsDA9HTv6NbtQ95cO7TwASDlaQGmpzLMTojmQYCeEEJfh6MFMXr1zMsnmAFRYeWVCR9w7X09wsJaEBB/UapWjS7xk1/ZrD0BqvhPFxRYHVyOEqA0JdkIIcYlOpeYx+aZn+SanLQCPDY0k5oar0evVxMf7NPsWIQNG9EeBDbPNneT1uxxdjhCiFiTYCSHEJcjNLeWZfs+z6nQ0AHf2MzJ4+A1UVtpo29a72V5+/b88fT0JcqveUizpiyQsFquDKxJCXIwEOyGEqKPKSitzH57DyqMBAPx3gBdPTRlEXl45UVGe+Pu7O7jC+pMQU72f7f7dpygsLHdwNUKIi5FgJ4QQdfTd/DW8vroUOwp6J9gZO2UI2dllBAS4Exnp0SwaENdW30EdAdib7c7ZjALHFiOEuCgJdkIIUQfZp84y6elvKcOVSEMJz775AKdPl6DTqYmL827WiyXO5/4xN6LCSjZGfv10HVVVNkeXJIS4AAl2QghRS2Vllcy6dwoHK8NQYeXlN++iuLiKmBhPuncPaDa7StSF0ehOgl/1Jdj13++Ty7FCNHES7IQQohYqK60kfbqOxX9Wr3Qd0tuIm0FPRIQHbdt64+bm7OAKG86APmEA7Dhu4+zZUgdXI4S4EAl2QghxEVVVNvbtPsPL41eSjydeLhbufKQPfn5uxMR4olS2nDl15/PAMwNRYuNMlTcbv9pEZaWsjhWiqZJgJ4QQF2Cz2Tl0KI93HpnNtpIIFNh46eUb0WiciI72bBXbbMUkBBPnYQbgh+V/kJ8vl2OFaKok2AkhxAWkpZlY9cZKPk32BOD+/l6ExYUSHKzDaHRzcHWNQ6lUcNstbQD4I8VGZmaxgysSQvwbCXZCCHEedrud9HQzaz76iTmL06nCiS6BpQx9aiAuLk60aWNoUW1NLmbk5DtRYyHXqmft0iTZYkyIJkqCnRBC/E15eRUHD+bx3vgPeX7GbkrtLoS75DF74XCKiqqIjPTAw6PlrYC9kNAof7r6mQD48cudHD9eiN1ud3BVQoi/k2AnhBD/Y7fbycgoZsuf6cy8fxqzV5kpw4VInYn5Kx6iqFxJUJA74eEGR5fa6JRKBXff1Q6AbSedOXgwl6ysEgdXJYT4Owl2QohWzW63U1xswWSq4ODBPDau3ccr905nabInNlRc3aacxd8/gcJFh7u7M7Gx3jg5tc63zvueuwtfRR7ldjVfL/qV48cLsVqlYbEQTUnLX84lhBAXkJpq4siRAirLy9i7cBHzvyvhjD0EgPt7ufPEzIcpKrJQWWmlfXt/DAaNgyt2HO9Ab+7qUsE7O+Db33IY+ICJ2FivVrOIRIjmoHX+2SmEEEB2dgmHDuXh5KTgtxenMu1bJ87Y/XBXVjB1XCyjptzD6dPFlJVZiYvzJjBQ6+iSHW7ky/fiTzZlVSq++mI3mZlyOVaIpkRG7IQQrVJ+fhkbNpzi3Xe3cmz/KTLzY7CjpF2oE6/N/w/uWhdyckpp08ZAWJihRW4Xdini+vTkVr9XWZDtx0/rjjHsoR5ERhpwd1c7ujQhBDJiJ4RohfLyylg6/TMeHf45v/9+kox8O3aU9Iq38dHnw1Br1BQUlBMf70P79kYJdf+HWu3E7aP64EEhxRUKfvnlGLm5ZY4uSwjxPxLshBCtSmFhOZ++8ikT30ghr0SBkVxGsoRZib8zY9HD5OWVU1lpo2NHX2JiPFGp5G3y7654eBjdVXsB+PnbPZw+XYTNJq1PhGgK5B1LCNFqmM0VfLVgPS/OO0opboS55DMh8g+6tNVy1fRZZGaV4erqRJcufoSFta4GxHXhFWhkYO9AlNg4lFpGcnIWBQWyzZgQTYEEOyFEq5CTU8JPq3fx7MRfKUJLsIuZj795jCsWfkPgrC/JsegxGDR07uwnqzxr4cYnHyKeFAC+/+YAZ8+WOrgiIQRIsBNCtHBVVTaOHy/gj1+P8Oyo5eRa9Xirinhn4T0UlUJJSSWhoXq6dfOna1e/VrejxKWK6XMN1wdmA7BhQxrHjxdSUVHl4KqEEBLshBAtlsViZffubDYmpfLsyMWcKPfBhQpef603NhcPAgLc6dEjkA4dfAkJ0cvKzjpQKpXcPvpmvMin1AI//niUvDy5HCuEo0mwE0K0SDabncOH89ny2yGeffQzjhZ7oqGCVyd2xCMiCqPRlYQEHxmhuwzdRg3nStUeANZ9u4eMjGLZP1YIB5NgJ4RocaqqbBw7VsA3izcyceLP5Fi0eCrMvP1iAvHXdUWrdaJ9e6OM0F0mrY8XgweEosLKsVPlbN58GrPZ4uiyhGjVJNgJIVqUoiILO3dmMXfyV0x5I5kSmwuhzjl8/EE/Ym7sSWWljfh4n1a9NVh96vv8WBI4BMA3K3fJIgohHEyCnRCixTCZKti2LYNpYxby3hdnsKKik+4Mi1aOwCsmjoKCCuLivAkIkK3B6ktoj270ia2eW7dp8xlSUvKwWKwOrkqI1kuCnRCiRTCZKvj5y+2Mum0u3+2sDhb9g0/z7rcTUHn4UFhYQXy8N5GRHo4ttAW654X/4stZKqqUfPPlPhm1E8KBJNgJIZq9oiILvyz7jScf/oxjxR5oqGB830peWfUSpVXOmEwW4uO9iYrylKbDDaDjvUO4xnACgB++3sPJkybZiUIIB5FgJ4Ro1kpLK/lt6TqefvwrTlt90SrL+fCNntw26RHOZJZRWWmjQwcjUVGeKJUS6hqCyknFA2N7ocbCmQIF635IIS9P9o8VwhEk2Akhmq2Kiip+X/Q1Lz22jDRrEG5KC++9P5CAhHgKCiqIivLgyisDZHuwRtDrmVH0cD0MwKpFGzl1ykxVlc3BVQnR+kiwE0I0S1arjQ3zP+PNx98n2dYWlcLGa2/chHeQLzYbdO7sR7t20qeusbjp3RnxSHeU2DiYoWD9uiMcP14ofe2EaGQODXbz588nMTERvV6PXq+nR48e/PjjjwCkpaWhUCjOe1u5cuW/nnPYsGH/OL5///6N9ZKEEI1k78/bWPzU6/xkuxaAZydcTXR8EEqlkk6djAQH62SUrpHd8vwoOmmOA7Dyo185ciSf7GxZSCFEY3Jy5DcPDg5m5syZREdHY7fbWbx4MYMGDWL37t3ExcWRmZl5zvEffvghr732GgMGDLjgefv378+iRYtqvtZopF+VEC1FRUUV6elFLB03gy9stwJw/30J9B0Qh9lsoVMnP3x93R1cZevkadTz0IOJ7H6/mB2pKrZvScNodMPX103mNwrRSBwa7G699dZzvp4+fTrz589ny5YttGvXDn9//3MeX716NXfddRda7YV7UGk0mn8890IqKiqoqKio+dpsNtf6uUKIxlNcbCE5OYev3/2Gt1PaU4UTV3fz4ZFRV3D2bBkJCT4EB+scXWardsfLo/hiyTB+K+3AB+9son3HYCIiDBK2hWgkTWaOndVqZcWKFZSUlNCjR49/PL5z506Sk5MZMWLERc+VlJSEr68vsbGxjB49mry8vAseP2PGDAwGQ80tJCTkkl+HEKJhVFXZ2L07mzfmbOTNFTlU4kzHgHJenXMrZ8+WER5uICLCw9FltnrevjpGj70aA2ZyzAq++2Y/6elFMtdOiEaisDv4X9u+ffvo0aMH5eXlaLVali9fzk033fSP4x599FGSkpI4ePDgBc+3YsUK3NzcaNOmDcePH+f5559Hq9WyefNmVCrVeZ9zvhG7kJAQTCYTer3+8l6gEOKy2e12dm4/w+D+H3G6oPqS3hXOB3j9u5cotGgwGFzo2tUPV1dnB1cqAM5mmRkXexvLzdehd7Wz8NN7ueGGULy8XB1dmhDNktlsxmAw1CqXOPRSLEBsbCzJycmYTCZWrVrFAw88wIYNG4iPj685pqysjOXLlzN58uSLnu+ee+6p+e/27duTmJhIZGQkSUlJ9OrV67zP0Wg0Mg9PiCbsdLqZsQNnc7rAGy3F3OPxJ3dOfYxS3HB2hvh4bwl1TYiPn467xt7Cj9NPUVDmyY/fHiAmxkuCnRCNwOGXYtVqNVFRUXTp0oUZM2bQoUMH5s6de84xq1atorS0lKFDh9b5/BEREfj4+HDs2LH6KlkI0UisVhtZWSVMu+tltmR7o8DGxHs9+O9Xn1IS1AWFQkF8vDfe3hIYmhKFQkGPMSPorT8AwOpVe0lNNWE2V1zkmUKIy+XwYPd3NpvtnMuiAB9//DEDBw7EaDTW+XynT58mLy+PgICA+ipRCNEITKYKtm3LZMmTc1i8tTq43Xu9G10euBdTsY2YGC969AgkJESmSzRFRn899469GT1m8ksU/LBmH5mZJY4uS4gWz6HBbuLEifz++++kpaWxb98+Jk6cSFJSEvfdd1/NMceOHeP333/noYceOu854uLiWL16NQDFxcU888wzbNmyhbS0NNavX8+gQYOIioqiX79+jfKahBCXr7LSyuHDeRz96TfeXpFBBRriA6z899k7USiUdOniT3y8N1qt2tGlin+hUCi4+rHh9NIfAuCrz5M5caKQkhKLgysTomVzaLDLyclh6NChxMbG0qtXL7Zv3866devo06dPzTELFy4kODiYvn37nvccKSkpmEwmAFQqFXv37mXgwIHExMQwYsQIunTpwsaNG2UOnRDNSGqqiS8W/cnTr+zgNIG4O1cx/b37KC+3EhfnRWCgVpoPNwNGPx33P94fd0o4W6zk+6/3cvp0saPLEqJFc/iq2KaoLqtPhBD1x263c/KkmUnjv2bZ6jQA3JUVzHj9JkLC/QgN1dOxo680u21GzmYXMSLufr4t7Iy/tpJFK4dzzTXBuLvLaKsQtVWXXNLk5tgJIVonm83OiROFTHnis5pQd6PbHr5cdjOx8UHodM5ER3tKqGtmfHy1jJx8N+6UkFXszOpP/5BROyEakAQ7IYTDlZdXsXdvDrOf/5pP1mQDcJPPIaZ++SLakDDKyqzExnqj08koT3OjUCi4ethgBoRUbxH5xapDHDmST3GxzLUToiFIsBNCOFRxsYXdu3NYMO93FnxxEoA+fieZ/OU0SlQenD1bRps2BoKCLryVoGi6PD1dGP3aQ3hgotCiYfm870lPL3J0WUK0SBLshBAOYzZXsGtXNm/N2cB7iw5iR0FPbQqTPhlPVr4dpRI6d/ajXTtvWSzRjCkUCjr17c7gTpUArPk1mwP7sjCZpK+dEPVNgp0QwiFKSizs2ZPDzOm/8sVXRwC4Ub2DSe+PIrdUTViYnu7dAwkN1aNSyVtVc+fp6cKYuWPxU+ZRatOweMZKDh/Oo6KiytGlCdGiyLulEKLRVVXZ2LfvLDOm/sqPP6WhwM6dLut59r3RlLoFEBfnRWKiUfrUtTCxncK5r783AOuSq0jecpxDh/KwWm0OrkyIlkOCnRCi0eXklPLOvC38+HN1qLvH9RdGzp9EsWcUEREGYmO9ZJSuBdJq1Tw8ZzTh6rNU4szC17/j5Ekz2dmlji5NiBZD3jmFEI2qqsrGwoW7WfbZfgBud/6FYe9OokgfTkSEgfh4bwl1LVh4hBdPPNoJJTa2pbux5addHD9eSGWl1dGlCdEiyLunEKJR7duXw6uvrAfgGsVWhs9+lBKPCKKjPUlIMOLsrHJwhaIhaTRODBl/J739TwPw/kd7SE83k5Ehve2EqA8S7IQQjaakxML4USspq1QQwhmeeOY6ykO7Eh5uoG1bb5yc5C2pNQgI0PLErPvxIZ9Ci4aVH/zE8eOFlJfLQgohLpe8iwohGkVlpZVPFmwnaVs+CuwM72ZCeeVggoK0tGvnI6GuFXFyUtK57xXcc0V1kPv293wO7s/i9GnpbSfE5ZJ3UiFEo9izJ5sXn/sBgKtdD9Ll6Wfx8HAhPt4HtVouv7Y2vr5u3D3tUeJUaVhRsejNtaSmmigpkR0phLgcEuyEEA2uqKiCUXd8QH6FGh/yGPPibSjdqxdKyDZhrZNSqSCmQzgPDAlHiZVdaZC07iAnT5odXZoQzZoEOyFEg5vxxBJ2nlSioooJ//HBJa4bcXFe+Pm5O7o04UBGoxt9n3mIa7UpACz64E+OHSugoKDcwZUJ0XxJsBNCNKgjh3N4f3H1B/ftcQVE3H434eF6wsMNDq5MOJpCoSA61sjwp29CSzGZxc588eF6TpwoxGazO7o8IZolCXZCiAZz5kwRU0cvpMCmw6Ao5u6XHsbHx00aEIsaOp2aq/87iMExeQB8+UM6+3ZnkJVV4uDKhGie5J1VCNEgCgvL+eWXVL76vXrO1F1XO+Oic6dtWy9cXZ0dXJ1oSkJCdAx/80lCVNmU29UsnLmKo0cLZB9ZIS6BBDshRL2rqrKxZcsZnhu/hlKbhiBFNjcOH0REhAe+vjKvTpzL2VlFfNcIht8dAcDG4xp+X/0naWmykEKIupJgJ4Sod2lpJp564juy8q14UMiEm214BfkREqJzdGmiiTIa3Rg8YSjdvbMAmP/+Ng7syyI9XcKdEHUhwU4IUa+KiizMHD2fg0cKUWPhEdevCBh0H2FherRaaW0izk+hUBAZ6cGTM+5BSymZFXqWz/6c5OQc0tJMji5PiGZDgp0Qot7Y7XZ+nL+aZb9U7yAwOOAIPd78AK82IQQHax1cnWjqtFo1Pfp14v5ebgCs3m7n6K4UUlLyKSyUFihC1IYEOyFEvTl9soCXXvqVclyI9CzjoU9m4hQcTXy8D+7uMlonLi4oSMvwVx+ikz4DGyqmvbqJ7Oxijh4twGq1Obo8IZo8CXZCiHqRlVXC03e9yeHyQJyp5KVZgygqsREX50VAgIzWidpRqZTExvnw3Mwh+JBHgcWFd6at4fTpIjIyih1dnhBNngQ7IcRly80t5cPZP7Jqe/VbyqhbdBj8jAQFuUsjYlFner2GK2++ksdu1qCiiu1Hq9jw035OnDBhsVgdXZ4QTZoEOyHEZTGbK1g88yumv7kbG0q6e5zm9vF34+SkJDLSUxoRi0sSHKzjpsmPc7PnAQA+WJBMSkqeNC4W4iLkHVcIccny8sr48q3VvPj6PiyoiXfPYdqiMeQXWAgP1+Pt7eroEkUzpVQqiIoxcv+LQwkjnXKrisVv/0RqaqE0LhbiApxqc9C8efPqfOIHH3wQnU56VgnRUp0+XcT2Tcd5ceoflOJDmHsRb6wYQ1GlE1FRnsTEeDm6RNHMeXq60L7/NQz9Zh1Tk+DP/WVs3XiMoCAdUVGeji5PiCZJYbfbL7rTslKpJDg4GJVKVauTpqenc+TIESIiIi67QEcwm80YDAZMJhN6vd7R5QjR5Jw5U8TWrRlMfXQBydk6tMoyPll6NwqNO1FRnsTHe+PkJBcExOWrqKhi66ZUxt8ylZ3lkbTzKeH1pY/So0cAer3G0eUJ0SjqkktqNWIHsGPHDnx9fWt1rIzUCdFy5eeX8ccfp5n01NcczdahoooXH2+L0kVLcLCWuDgvCXWi3mg0TrTrGMzwJ/qQPOsoB3LdWbvwO3x976RDB1+USoWjSxSiSanVu+9LL72EVlv7dgXPP/88Xl5yGUaIlujYsQImPfcjR09bcKOUF27MJ753T3x8XGnXzge1unYj+0LUlre3K/1G3kKfiOrtxRauOsXe7alkZkr7EyH+rlaXYlsbuRQrxPmZzRUM6LeEP7dk4EYpT0du4dp3P8Ric6ZbN398fNwcXaJooSwWK0k/HWHo4A/JrvKgR0Ah05c9SffuAbi5OTu6PCEaVF1ySa2vl5SVlbFmzRqKiorO+w3XrFlDRUVF3asVQjQbH320iz+3ZKCiihG6H7lm5msUlSmJifGUUCcalFqtIqFzCI8+1B4FdjZnevD9+99w5Eg+VVWyI4UQf6l1sPvwww+ZO3fueefP6fV65s2bx0cffVSvxQkhmo6iogpmT1kLQB/FRvpOf4liJ0/atvWmTRsPxxYnWgV/f3dufugmev/vkuxHq06xZ+sJUlLyZbsxIf6n1sFu2bJljBs37l8fHzduHIsXL66PmoQQTdD0sYvJNivRUsz9Y/tSFZJIu3bexMR4ygR20SiUSgUxMV48NutBjE5mTDYtC6Ys58iRfE6cKHR0eUI0CbUOdkePHqVDhw7/+nhiYiJHjx6tl6KEEE1L1skcPlhyBIDbEkpwuXYI4eF6IiI8UCgk1InGo9Op6dQ9jEcfSkSBnT8yfdj00WccOVJAdrbsSiFErYNdVVUVZ8+e/dfHz549S1WVdAMXoiV64f55FNp1eKpKGDRpFAEBWuLivGW7MOEQQUE6bnukP/3bV+8b+/4PpZzZupVDh/IoKbE4uDohHKvW78rt2rXjl19++dfHf/rpJ9q1a1cvRQkhmo7j+0/y+aZKAO4d4Iur1p3oaE9cXGrdBlOIevXXJdknX3+QEPcSSnFjxoyNZB49w+HDMt9OtG61DnbDhw9n6tSpfPfdd/947Ntvv2X69OkMHz68XosTQjiW3W7n+aHvUYIbvs5mrrv/JiIjPTAaZQWscCw3N2fat/dj4rSBuCospFX5M//5hZxMLSAtzeTo8oRwmFr/yf3www/z+++/M3DgQOLi4oiNjQXg8OHDHDlyhLvuuouHH364wQoVQjS+5N8P8M3u6obD998Wio/RnfBwg4OrEqKan58bV14XzbhR6cycn8LG3GCiXvsY15cfxdPTBS8vV0eXKESjq9MEmaVLl7JixQpiYmI4cuQIKSkpxMbG8tlnn/HZZ581VI1CCAew2ey88N93qEBDiIuJHnf1JzzcIJdgRZOhUCgIDzfQ85bu3HejOwBLtrmxbc3vHD1aIP3tRKtU53fou+66i7vuuqshahFCNCE/v/cFP6f7APDgQ93w9HQhMLD2WwsK0Rg8PV0ID9czeMxAjhz6gG2ZBt5ccIjQLu0xGt2IiPBwdIlCNCpZ0iaE+IeK4hJmPPclVTgTa7SQeH1n2rSR0TrRNMXEeNGmjQePTbkTX2UBJqsbC6at5MiRAvLyyhxdnhCNqtbBbtu2bbRr147IyEhWrFjRkDUJIRyooKCcZc+8zh8lMQAMe7yXjNaJJs3JSUl8vDcxiWE8en8ECmxsPunKj4vWcvhwHuXl0opLtB61DnajR49m6tSp/PLLLzz00EP1si/s/PnzSUxMRK/Xo9fr6dGjBz/++GPN49dffz0KheKc26hRoy54TrvdzosvvkhAQACurq707t1bGicLUUtlZZXs3nGGpYuTq0frAhVEtw0iIsJDRutEk6bROJGQ4MOND93OgPDqnqsffHGaret2s2dPjvS3E61GrYNdXl4ewcHB+Pn5UVFRQWlp6WV/8+DgYGbOnMnOnTvZsWMHN954I4MGDeLAgQM1x4wcOZLMzMya2+zZsy94ztmzZzNv3jzef/99tm7diru7O/369aO8vPyy6xWipcvIKOabuV+wqawtAP95qCfe3q4EBLg7uDIhLs7dXU379j6Meesxol3PUoGa2TOTOLLvNLt3S7gTrUOtg90LL7zA/fffzw033MADDzyAp6fnZX/zW2+9lZtuuono6GhiYmKYPn06Wq2WLVu21Bzj5uaGv79/zU2v1//r+ex2O2+99RaTJk1i0KBBJCYmsmTJEjIyMvj6668vu14hWrKyskomT/6Ned9bqERNQpCVdolBhIcb0GhktE40Dx4eLrTvHMxzrw7GoCwmu1LPWxOXkZ1VTEpKgTQvFi1erYPdyJEjWb9+PQsXLuSjjz6q90KsVisrVqygpKSEHj161Ny/bNkyfHx8SEhIYOLEiRccKUxNTSUrK4vevXvX3GcwGOjevTubN2/+1+dVVFRgNpvPuQnR2nzzzRE++2wfADc4bePVNwfj4aHB319G60TzEhyso1uvRMY/koCKKnblevHlzE84edJEaqo0LxYtW53+DA8JCan3Avbt20ePHj0oLy9Hq9WyevVq4uPjAfjPf/5DWFgYgYGB7N27l2effZaUlBS++uqr854rKysLAD8/v3Pu9/Pzq3nsfGbMmMGUKVPq6RUJ0fxkZhbx7JPfAtCdnYz+bxSVKndCQ/Uyt040OwqFgpgYT/oO7UNaymkW/VrFip0uhH3xA2r1rWi1avmDRbRYtXrHNpvNF7wE+ndFRUXodLpaHRsbG0tycjImk4lVq1bxwAMPsGHDBuLj48/ZyaJ9+/YEBATQq1cvjh8/TmRkZK3ruZiJEycyfvz4mq/NZnODhFghmqLy8ipeHL+aU1kVuFHKfd0rcb99DGoXNQEBshJWNE8ajRPt2vlw/wv3k3ribZLSdLyx4izG8J24uFyBq6sTBoPG0WUKUe9qdSnW09OTnJycWp80KCiIEydO1OpYtVpNVFQUXbp0YcaMGXTo0IG5c+ee99ju3bsDcOzYsfM+7u/vD0B2dvY592dnZ9c8dj4ajaZmZe5fNyFai3XfHeSTFdX/pu4IOkLw+NloXDW0b29Eq1U7uDohLp1OpyYuzpsxs4bTVp+PBTXTZm8n9cApDh7MpaJC2qCIlqdWI3Z2u52PPvoIrbZ2f71XVlZeckE2m+1fW6kkJycDEBAQcN7H27Rpg7+/P+vXr6djx45A9ejb1q1bGT169CXXJERLlZdXyuPDl1OFK23Vpxg86zk8/Qx07OgnoxmiRfD3dycy2ptnZ9/Ds4+tJLvKwPSJq5n+/jDc3Jxp396IUqlwdJlC1JtaBbvQ0FAWLFhQ65P6+/vj7Ox80eMmTpzIgAEDCA0NpaioiOXLl5OUlMS6des4fvw4y5cv56abbsLb25u9e/fy5JNPcu2115KYmFhzjri4OGbMmMHgwYNRKBSMGzeOadOmER0dTZs2bZg8eTKBgYHcdtttta5fiNbAbrcz6cH3OVXkiitlPPPcDWi8jMTH+0ioEy2GUqkgMtKDvLwAXn75Op6atJm0Ug/mPreUCfNG4ubmRFSUJwqFhDvRMtQq2KWlpTXIN8/JyWHo0KFkZmZiMBhITExk3bp19OnTh/T0dH755RfeeustSkpKCAkJYciQIUyaNOmcc6SkpGAy/f9VThMmTKCkpISHH36YwsJCevbsydq1a3FxcWmQ1yBEc3X0UDZLv8sBXLnvKjuenboTG+uF0ejm6NKEqFd6vYZ27XywWNrz3IgMpnycybYMLctmfYZ68n1oNE6EhsoUHNEyKOx2u93RRTQ1ZrMZg8GAyWSS+XaiRSotrWRE39ms+KMKo8rEe0sfJCzSh86d/VCpZAtp0TKdOmVmz54cfnh9Ee9vqJ4/+lA/VwY9djtt2njQpo0BN7eLX20SorHVJZfIO7gQrYzdbmfT72l8/Wd1T8h7btCi89QRHe0poU60aCEhOtq186HfuGEMis0H4ON1pfy65AdSUvLZtStbdqcQzZ68iwvRyhQWlvPcYyspt6sJVWXRc+htREYa8PCQ6QqiZVMoFLRpYyAhwYcHXh3NjSF52FEwd2Ue277+jdzcMvbty6W8XFbLiuZLgp0QrcyH7+9g9/FKVFQxaoArXgFeBAZKvzrROigU1Ysp2sYbGT3nEa72z8WGktmfZrLxqw1kZBRz4EAuFovV0aUKcUkk2AnRihQVVfDa9F8A6KfeStQddxMeLvOKROuiUCiIjvakbTs/nnjjYa71y8SOgjc+PUPSyg2cPGnmwIFcKisl3Inm55KC3caNG7n//vvp0aMHZ86cAeDTTz9l06ZN9VqcEKJ+fTznB/JKlLhTwj2P9ccr2J+gIBmtE62PSqUkLs6L+MQAHnvzUXr5nQLg7RUZrF3yM6mpJg4dyqeqyubgSoWomzoHuy+//JJ+/frh6urK7t27a5oJm0wmXn311XovUAhRPyrKKvhwzk8AXBtSiPvVtxAd7Ymrq4zWidbpr3AXFefHqDcfZ0BQOgALvsllzYLvOXasgMOH87BaJdyJ5qPOwW7atGm8//77LFiw4JwmxFdffTW7du2q1+KEEPXnyykfcKjUHwV2bhl7NyEhOhmtE62eSqWs7t8Y6MWoeU9yW5vqq1CL15pZOe9rjhzJJyUlX8KdaDbqHOxSUlK49tpr/3G/wWCgsLCwPmoSQtSziuIS3ptbPVWiUziEtg0jIsJD2psIAbi5OZOQ4IO7wZ3hbz7FXbFZAKzYUMbS177k8OE8DhzIIz+/DGn9Kpq6Or+r+/v7c+zYsX/cv2nTJiIiIuqlKCFE/frsmTlsLm8LwJ0jb8DHxw1PT2lvIsRfvL1d6dzZH72HK/fOGMf9ibkArN5SycdTV3DiRAGbN2dw5EgBNpuEO9F01TnYjRw5kieeeIKtW7eiUCjIyMhg2bJlPP3004wePbohahRCXIbSvDzeXHAUG0q6xTgT2zaQ4GCdbHwuxN94errQubM/fv5abnvxMR7sWgjAD8kw99klaDRKDh3K48iRfAl3osmq1V6x/9dzzz2HzWajV69elJaWcu2116LRaHj66acZO3ZsQ9QohLgMn06cz15rNApsjJowAJ3OGR8fV0eXJUSTpNOp6dTJD7Vaxc3PjULz1kI+3ORC0hFnMu96j1kLhnL4cD5Wq524OC+ZziCanDr9RlqtVjZu3MiYMWPIz89n//79bNmyhbNnzzJ16tSGqlEIcYnsNhsfLK2eOnF9ezUenu4EBmpxcanz33RCtBouLk4kJhqJjfWi7/iHeG6IBhfKSSnQ8tB/PqUwI5uUlHwOHsyTdiiiyalTsFOpVPTt25eCggLUajXx8fFcccUVaLWysk6IpmjdB1+TXBYGwB0jrsfLy4XwcIODqxKi6XN2VhEf70NCgjdXPnAXs5+KwUtRSK7FnUcf+4Gj2/Zx9GgB+/eflV0qRJNS5zHkhIQETpw40RC1CCHqkd1u541Xf8GOkvZ+ZYRE+BEb6yW7TAhRS0qlgqgoTxITjYRfdw3vzutHhDqbMruG52ftJenz9Zw4YWLfvrOyv6xoMi6pj93TTz/Nd999R2ZmJmaz+ZybEKJp2LXhIEmnPQEYeEdHQkJ0+Pu7O7gqIZoXhUJBaKietm29cQmOZP7nI+nueQYbSt7+IoclM7/gxIlC9uzJITe3VC7NCodT2OvYlEep/P9ZUKH4/6vq7HY7CoUCq7X5D0mbzWYMBgMmkwm9Xu/ocoSoM5vNzq3tJ/LDQVfCXQt567NHueqqIIxGN0eXJkSzZLPZOXAgl+PHC3HX2Hj/ibl8ddQPgC5BZTw9+37ctC54eGho29Zb/q2JelWXXFLnGdS//fbbJRcmhGgcOzYc5KeD1Zdc7709Aj8/d7y9ZSWsEJdKqVTQrp0POp2aw4fzGfb6k4TO+4R3flGy84wrzz28mNnv3UWRyoudO7OJj/cmJER3zgCIEI2hziN2rYGM2InmzG630z/uWX464k6UWy6zPn2M7t0DCQrSObo0IVqE3NzSmp0ojq1dxysfplOMOxqFhcf+G0H//95IaamV2FgvoqJkhxdx+Rp0xO7333+/4OPn225MCNF49v15iF+PqAG47552eHu74esrl4WEqC8+Pm506+ZMSko+igH9mRu2lxlTfuGYxZ/Xl5zm998X8cKb93HwYC4VFVbatvXC2Vnl6LJFK3FZc+xqTvJ/hppljp0QjnVPhyf5fK8Hoa6FvPHpKLp1CyA0VH6PhahvVquNtDQTKSkFVBSZ+fzlD1l12AcrKjzUFl6Z3ofANgGEhxuIj/fG1VVWpItLU5dcUufx4YKCgnNuOTk5rF27lm7duvHTTz9dctFCiMuXlpzCj3urRwZuGxiLj4+rrIQVooGoVEoiIz3p2tUPTz8v7p05njdGavEjl0KLmnHP/Ma6lZtJTTWxc2c2eXllVFRUyXZkokHV+VKswfDP5qZ9+vRBrVYzfvx4du7cWS+FCSHq7vVRb2PGiMG5nKsGX0d4uAG1Wi4BCdGQfH3d0WicSE7ORnnr7SzsdoSXn1jO9rJoFnx+kuRdZ3jipUGYTBU4O6swGNRER3vi5SULmkT9q7cZnX5+fqSkpNTX6YQQdZR1MIXVW6ubpN7SJwQfH1f8/GS0TojGYDBoSEz0Ra9XU+4dxasrnmFom8M4UcX2o1U8PnwZx5OP4uKiIienlK1bMzl1Snq/ivpX5zl2e/fuPedru91OZmYmM2fOpKqqik2bNtVrgY4gc+xEczTlhpG8nBSMWmllwaK76N49kNhYL0eXJUSrUlZWSWqqidRUE3arlUPLFjP7i2Jyqf63GB+s4JFn+tG2rS82G3Tt6i+tiMRFNeiq2I4dO6JQKPh7HrzyyitZuHBhXU8nhKgHpqxcVv5eCkCvq4z4+LjJ3DohHMDV1Zn4eB+MRjdSUvKJvPsBFvRJZ86zn7E5P5CDp1U88cRaunQ0MnZcTw4ccKJzZz+0WrWjSxctRJ2DXWpq6jlfK5VKjEYjLi4u9VaUEKJuVk5fyEFbNAC33NuDgAB39Hr5oBDCUYxGN/R6NceOFXLihJLnlzxH1refseTjrWyq6sTO5LNMe2kdr8y4GaVSQWKiEb1e4+iyRQtQ5zl2GzZswN/fn7CwMMLCwggJCcHFxQWLxcKSJUsaokYhxAXk5pYy66NT2FHQPQqCggwEBmql470QDqbROBEf703Xrv64azV49L2bSUsm8HzAd2gp5tjJUt58+WuysorZtSub3NxSR5csWoA6B7sHH3wQk8n0j/uLiop48MEH66UoIUTtnD1byvSnV3Ks3IgzlTwy7ka8vFxlzo4QTYRCoSAgQEv37gFER3tS7BrA1fMXM7bdkeqFFUcqef6RTzh5PIedO3NISzNRWFj+j+lOQtRWnYOd3W4/70jA6dOnz9sKRQjRMCwWK998c4QPlx0HYFBUDoYAP0JDdTg5yRZGQjQlLi5OtGvnQ2KiEaWbni4vzeTpvhW4UsbRfFeeeGQVB/7Yy+7d2fzxxxkOHsyjsrL5N/wXja/Wc+w6deqEQqFAoVDQq1cvnJz+/1OtViupqan079+/QYoUQvzTL9/uZdyjX1Fa5UQYp7jzwevQ6dQYjbJ9mBBNkVKpICLCA39/d86cKcJ5zMOEd97Iy7N3k2X1ZsLU7YwceIJ7xt9OSko+FouVtm29cXGp83R40YrV+rfltttuAyA5OZl+/fqh1WprHlOr1YSHhzNkyJB6L1AI8U+mswU8eu9HlFT6EsIZXhiiRh1/JTExnri5ybZFQjRlbm7OREV5AmC39+T9T2J5+bFPSTb78/4aE79u/pBpb99NaqoJi8VKQoIP7u6yGErUTp372C1evJi77767Ra+ClT52oimz2+1MHfo6Ly0tQU0lS+Z0wrlNe6KjPWnXzgelUhZNCNEc2Gx2TpwoJCUlH1ulha9fXcCy7S5UoMFVVcXEZ68ioWsURqMr7dsbMRhk1Wxr1aB7xT7wwAMtOtQJ0dSdOmVmwcozAAzopEAd0Z6QEB2xsV4S6oRoRpRKBVFRnnTt6o/e051bJo7inScDCVeepszqxIuvbmPxu7+QlVXCzp1Z5OSUOLpk0QzUOdhZrVbmzJnDFVdcgb+/P15eXufchBANp6jIwpsvf8vpCg80VDB4RG/8/d1p185H9oQVopny83OnW7cA2rQx4H1NX977+G766PcD8M36HJ4b+zknTuSzfXsWR4/mU1RkkYUV4l/VOdhNmTKFN954g7vvvhuTycT48eO5/fbbUSqVvPzyyw1QohDiLzt3ZvLB0iMA3BSejU94MLGxXjKvTohmzs3NmfbtjcTEeFGmD+X5L17hycQU3CjleJadx0evZtufx9m3L5eNG0+zeXMGZnOFo8sWTVCdg92yZctYsGABTz31FE5OTtx777189NFHvPjii2zZsqUhahRCABUVVYx97DvKq1SEcYrB9/cgKEiLl5f0rBOiJVCplMTFeREX54W5wpler77Cm8PUhJFOWZWKV17dzOJ31+PkBAUF5aSk5FNVZXN02aKJqXOwy8rKon379gBotdqaZsW33HIL33//ff1WJ4So8fmiLew/kIsaC8N8/sSz27WEhsriHiFaEpVKSdu23lxxRQBarQafW//LvLf60tstGYA1v+Zw+8BP+f33VNLTi0hLM0kzY3GOOge74OBgMjMzAYiMjOSnn34CYPv27Wg0smJHiIZgOnmSN59cAMAV6kPEPjWF4FADnp6ykEmIlkahUODv784VV/gTFeWJPbwTTy55iac6HcWPHEoq4I03NrNrxykOHsyTcCfOUedgN3jwYNavXw/A2LFjmTx5MtHR0QwdOpThw4fXe4FCCPjm6VfYVx4GwJ0vPIix6xVERXnKfrBCtGCurs4kJPjQpYs/Gg8Prpr0Im897MG1/AnAnFkbyD11hv37czl+vBCbTcKduIQ+dn+3ZcsW/vzzT6Kjo7n11lvrqy6Hkj52oikpycnhnsA7+M7ai5gwV2a9NZguXfwICZHfTSFai+JiCydOmDh1ykTRtiQmTd/CCVsIXqoiXn2xB36JiURGehAb6yUr5FugButjV1lZyfDhw0lNTa2578orr2T8+PEtJtQJ0dRsmP0Of1g7A9D31kT8/NwICtI5uCohRGPSatW0b+9D164BGK/tzbNTB2J0MpNv1fHkS9vZvXwlR47ks2dPDiUlFkeXKxyoTsHO2dmZL7/8sqFqEUL8TWV5OfPe20UBnnjrlVx5ZQhBQTppRCxEK/TX3LsuXfyJ7tGR19+7g0idmTJcmbrSwuIJc0lLyWTbtiyOHy+guFgCXmtU5zl2t912G19//XUDlCKE+LvV0z/ml7JOADw27mp8fNzx9pYFE0K0Zp6eLiQk+GAM8Wfu8kfom2DDjoJvT/jx/MiPOLljL3v35rJ1ayZZWbJbRWvjVNcnREdH88orr/DHH3/QpUsX3N3dz3n88ccfr7fihGjNrJWVzHx9B1bC6Rhqo3OXEPz8XGUzcCEERqMb8fHe7N2bw/g5D9Bt1QbmfHyMY+V+PDZhIxOGHuWKewexc2cWcXFehIcbUKnqPJYjmqE6L55o06bNv59MoeDEiRO1Ptf8+fOZP38+aWlpALRr144XX3yRAQMGkJ+fz0svvcRPP/3EqVOnMBqN3HbbbUydOhWDwfCv5xw2bBiLFy8+575+/fqxdu3aWtcliydEU7Bi0rv8Z/pZ7ChY9EF/PP28ufLKAPz83C/+ZCFEq5CZWczBg7kUFVVSkpPNxCe/Jr20eg7urTH5jJkzkjKrhrAwPW3beuPiUufxHNEE1CWX1Pkn/H8XTlyu4OBgZs6cSXR0NHa7ncWLFzNo0CB2796N3W4nIyODOXPmEB8fz8mTJxk1ahQZGRmsWrXqguft378/ixYtqvla+uuJ5sZqsfDWG39iJ4ZubezofTyJiDBgNLo5ujQhRBMSEKDF3d2Zo0cLSAfe++xh3pm4mHUH1Xx7xIsDQ97l5UlXkWrrjsViJTHRiKurbEHYkl1yuxOLxUJqaiqRkZE4OdXfXwBeXl689tprjBgx4h+PrVy5kvvvv5+SkpJ//Z7Dhg2jsLCwTvMAKyoqqKj4/3vumc1mQkJCZMROOMzXL81jyCu52FAx69XrufaGaDp18kWjkb+2hRD/ZLXaOHOmmJSUfEpLK9n27QbeWJhKuV2NG6WMu1VD1wfvx9/fnfbtjWi1MqWjOWmwdicApaWljBgxAjc3N9q1a8epU6eA6mbFM2fOvLSKAavVyooVKygpKaFHjx7nPeavF3SxIJmUlISvry+xsbGMHj2avLy8Cx4/Y8YMDAZDzS0kJOSSX4cQl6uqvJx35iRhQ0VCCLRNCCIqykNCnRDiX6lUSkJD9XTvHoDR6MYVt17Hp8vuINyjnFLcePVbFQuf/5D0UyaSk3MoKCh3dMmigdQ52E2cOJE9e/aQlJSEi8v/X53Xu3dvPv/88zoXsG/fPrRaLRqNhlGjRrF69Wri4+P/cVxubi5Tp07l4YcfvuD5+vfvz5IlS1i/fj2zZs1iw4YNDBgwAKvVesHXZDKZam7p6el1fh1C1Jc/3prPn6UxAAy67yr8/d3x9nZ1cFVCiOZAr9fQvr0RDw8Nzm5aPv32UW7qWH1h7rvDOiY98jFHD2SwdWsmycnZFBVJS5SWps6XYsPCwvj888+58sor0el07Nmzh4iICI4dO0bnzp0xm811KsBisXDq1ClMJhOrVq3io48+YsOGDeeEO7PZTJ8+ffDy8mLNmjU4O9d+fsCJEyeIjIzkl19+oVevXrV6jiyeEI5SWVrKUL+bWVF8PX4GeG/hPVx1VTD+/rJgQghRewUF5Rw6lEdWVgkGg5r1C77izc/PUo4LXk4lzJjRG9+oNnh6akhM9JV9p5u4Br0Ue/bsWXx9ff9xf0lJySXtW6lWq4mKiqJLly7MmDGDDh06MHfu3JrHi4qK6N+/PzqdjtWrV9cp1AFERETg4+PDsWPH6lybEI1t65vzSCpuC8Att3fA31+L0SijdUKIuvH0dKFbN386dvTFbodud9/KvJc74qMsJL/KnbHPbGDvL5sxmy3s2pVNbm6po0sW9aTOwa5r1658//33NV//FeY++uijf50bVxc2m61mIYPZbKZv376o1WrWrFlzzqXf2jp9+jR5eXkEBARcdm1CNKSy/Hw+mfElWfihcYZrr48iJEQnvaeEEJfE2VlFRIQHPXoEEh3tgW/HzrwzfxCxbtlYUPPyuydY+dYqysoq2b07R5oZtxB1no396quvMmDAAA4ePEhVVRVz587l4MGD/Pnnn2zYsKFO55o4cSIDBgwgNDSUoqIili9fTlJSEuvWrasJdaWlpSxduhSz2VxzmddoNKJSVW9yHBcXx4wZMxg8eDDFxcVMmTKFIUOG4O/vz/Hjx5kwYQJRUVH069evri9ViEa19ulJ/FDSFYDbb48nMFCLj4+M1gkhLo+7u5p27XzQ6zUcOuTMqx8/wrtPf8CvZ/z49NcKUlM/4unXh5GcnENiopHAQK2jSxaXoc5DAT179iQ5OZmqqirat2/PTz/9hK+vL5s3b6ZLly51OldOTg5Dhw4lNjaWXr16sX37dtatW0efPn3YtWsXW7duZd++fURFRREQEFBz+7+LG1JSUjCZTACoVCr27t3LwIEDiYmJYcSIEXTp0oWNGzdKLzvRpB37YydffbKRTPxxc1Ey4Oa2BAZqpd+UEKJeKBQKQkP1dOvmh2+wN2PmPsHwHqUosbEp1ZXxQxeQd+Ysyck5pKebucROaKIJuOQ+di2ZLJ4Qjamioor3evTlld0dKcTA8OGduP32dvToEYiXl4zYCSHqV0mJhV27qlueHFjzA1M/zqAMVzxUpbw8+RradIojLs6LiAgPnJxkKkhT0KA7T0B1z7nVq1dz6NAhAOLj4xk0aFC9NioWorX4c8FnLNgdQCEGggPd6Ns3Gn9/d1mlJoRoEO7uahISfNi1K4uEQTfxQcx+nnk+iWyrJ8+8vIUn7j5B1V39yc4uJSxMT1CQVub6NiN1HrE7cOAAAwcOJCsri9jYWACOHDmC0Wjk22+/JSEhoUEKbUwyYicaS9bpAoa3fYAfi7vgrLLz4Ue34+3txpVXBsr2YUKIBpWXV8ahQ3lkZ5eirshj4iOfcqDYD4Db2hby0IyRVNqdiYjwIC7OC7Va5eCKW68GbXfy0EMP0a5dO06fPs2uXbvYtWsX6enpJCYmXrR5sBDi/yspsTB96Bx+LK6em/rMUz3w9nbD399dFk0IIRqct7crXbv6ExPjic3dyMylY7m5bfXK2K8PeTDp/jdRFaRz7FgBe/bkUFFR5eCKRW3UOdglJyczY8YMPD09a+7z9PRk+vTp7N69u16LE6KlslptfPbh78z/rfqf4O1XudDv5ngUCggPN1xST0ghhKgrFxcnEhJ86NrVDxd3V0bOeJjx9xpRYSW5OJhRI1Zh3vk76elF7N+fi8Xy77s4iaahzsEuJiaG7Ozsf9yfk5NDVFRUvRQlREt38EAOE5/9CStOdNBn8fRr/yE3t4yoKE8ZrRNCNCqFQkFQkI6uXf0xGl25+o7evPFqT3ROFjLsfox9NYXj69Zx8qSZAwck3DV1dQ52M2bM4PHHH2fVqlWcPn2a06dPs2rVKsaNG8esWbNq+s3VdWsxIVqLkhILowa/TW6lO54KE9PevofsnHICAtyJiPCQ0TohhEN4eLjQpYsfcXFetEmI4P1FdxGqK6MYLRM/zCFp4Zekpppk5K6Jq/PiCaXy/2fBvz6A/jrF//1aoVBgtTbPH7wsnhAN6fWnl/L068dRYuPVBw1E33ILAQFa2rXzRq+XfotCCMey2+2kpxeRnJyDs5Odlx76gO2ZOgBub1/MvZOHExnlRUKCD87OsqCiMTRou5PffvvtkgsTorU7uXM/r725DfCmd3g+MbfeT3i4nvh4bzQaaRckhHA8hUJBSIiOsrIqDh7M5dUlY3h/wgK+3K3hq31azOPf5uHXHkWhgNhYL2mk3sRIg+LzkBE70RDKCwt5IOq/fJHXFa2qgo8W34VvoBdXXOGPu7va0eUJIcQ5rFYbqakmjh8vpLzcytoFX/P+99WrZnv4nuWJNx7GL8iD2FgvAgJkG7KG1OANisvLy9m7dy85OTnYbLZzHhs4cOClnFKIFi9p4iR+yqteYPTQyK646rTExnpJqBNCNEkqlZKoKE+8vV05dCiPPsMHovX4hTeWZbM5x0jZ2AW89P5D7C6twmazExSkc3TJgksIdmvXrmXo0KHk5ub+47HmPK9OiIZkKS5m9Se/U8gQtG5KrroxnpAQHQEB7o4uTQghLsjT04VOnXw5eDAPbu/NCy5JzPj4FMkFPkx56F1efm84e/dWZ4DAQBm5c7Q6r4odO3Ysd955J5mZmdhstnNuEuqEOL/khYv5o7wdAP1vikOrVdOmjfSrE0I0D66uznTo4Et8vDddB97IC6NjqnvdFQXx4sgPKD2dxt69Z8nKKnF0qa1enYNddnY248ePx8/PryHqEaLFqaqy8utbizhEDADXXBNOWJgeDw/ZC1YI0Xw4OSmJjfWiQwcjV9zSk+fHJeKElb3l4Uwcs4KCI4fZu/csubmlji61VatzsLvjjjtISkpqgFKEaJn2fvsbq1MDsaGkfTtvoqK8ZS6KEKJZUigUhIUZ6NTJlx59OzP9patwUVRytDKIZ8d/x5k9B9i3LxezucLRpbZadZ5j984773DnnXeyceNG2rdvj7PzucucH3/88XorTojmrshczscTP2AbXVBg5977OmE0uqLXy4IJIUTzFRCgpaysioqKSOa95cb48WtJt/oyceLPTJutQa1W0rGjrywOc4A6B7vPPvuMn376CRcXF5KSks6ZI6RQKCTYCfE/drud5SOfZUlKOAB33BREXJyRwECtzK0TQjR7YWF6TKYKbDY7779/C2NHfU221ZOXn/uOaW9qUCoVJCYaJdw1sjpfin3hhReYMmUKJpOJtLQ0UlNTa24nTpxoiBqFaJZ+nfUO078opRgt4X4qho+9Eb1eLXvBCiFaBJVKSVycF2Fhely9fHjjzX7olaWcqfTk5fGfc2jzQXbtyuHsWZlz15jqHOwsFgt33333OVuLCSHOVZxzlucn/UY6wWg1dt6YfyelpVaCg/Wyw4QQosX4v6tlvcLCePWVnrgryjlpMTL52W85tnkP27dnkZ4u+8c3ljqnswceeIDPP/+8IWoRosVYNP4ttlk7oMDOrDduwcXFGW9vV8LCZCcTIUTL8tdq2Y4dfWnTqR2zpl+PTlVGRpU3k19YS/GZdA4cyJORu0ZS56EDq9XK7NmzWbduHYmJif9YPPHGG2/UW3FCNEfF+SbeW5EBhHNDJ1fatfPFbLYQE+OJi4uM1gkhWh6FQkFoqB5nZyUqlYLX5rjwzNM/cqbKh+cf/4IZ7z/A3r1KOnf2w9NTWj01pDp/yuzbt49OnToBsH///nMekwnhQsB7j7/DYWs4SmwMe6IvJpOF2Fgv/P1llwkhRMsWEFC9OEypVPDy5GuZNGUTx8t9eWXMIl58Zyh7VAo6dvSVPp4NqM7B7rfffmuIOoRoEXIyTbzzeSZgpE8HJwyeOjp29CUkRCd/+AghWgV/f3ecnZU4OyuZ+EQRU9/ax/4Sf2aN+YAJrz/AHoWCTp180es1ji61RbrkFRDHjh1j3bp1lJWVAdWtHYRozex2O68Me4v0KiMuigrufrQ/ISE6CXVCiFbH29uVDh18uaJ/d14Y3wklNnaWtmHGuEVkHEhh796zlJRYHF1mi1TnYJeXl0evXr2IiYnhpptuIjMzE4ARI0bw1FNP1XuBQjQXKcmpLP6penLw3Tfo8A2qbgMgoU4I0Rp5eroQHe1BYs/2THn+CpywklwRyfTxn3Fy7xH27culvLzK0WW2OHUOdk8++STOzs6cOnUKNze3mvvvvvtu1q5dW6/FCdGcvPTAOxTjhp+zmZsfHkhQkFYmCQshWrWwMD2hoXriu0QzY8pVOCuqOFARyivjP+fIjhT27Tsr4a6e1TnY/fTTT8yaNYvg4OBz7o+OjubkyZP1VpgQzcmRP3bzwz4VAPffFY1W70pYmEFG64QQrZpKpSQhwYewMD0RCRHMfuVqNIpKjlQEMvXplRzacYxdu7IpKCh3dKktRp2DXUlJyTkjdX/Jz89Ho5GJkKJ1mjXiTYrR4qmuoPug6wkN1ctonRBCABqNE+3bG4mIMBDWLoLZr1yFi8LCsQo/po1fwcmDJ9m5M4vcXOlzVx/qHOyuueYalixZUvO1QqHAZrMxe/ZsbrjhhnotTojm4NCaH/kupbrx8ODbYvH0dCE8XBoRCyHEX9RqFe3bG0lI8CE8MYbpk7vjoqjgaLkvM59YhPlMJnv2nJWRu3pQ53Yns2fPplevXuzYsQOLxcKECRM4cOAA+fn5/PHHHw1RoxBNls1q5b0xb5BDT1ydbNxwa1fCwvSy6bUQQvyNSqUkKsoTjcYJJyclLz5TxYuzd7O3NJg3xr7L+LcfZa+Tgs6d/dHp5D30UtV5xC4hIYEjR47Qs2dPBg0aRElJCbfffju7d+8mMjKyIWoUoslKXriYH08HAXDLLTEYjW74+2sdXJUQQjRNf+1Q0amTL51v7Mz4UQkosPNHUTRvPvYOOcdOsm/fWUpLKx1darNV5xG7U6dOERISwgsvvHDex0JDQ+ulMCGaOpvVyqqX5nOcW1Aq7PS9uR0BAVr5S1MIIS4iIECLk5MSpbI7FZV25n18mE3FsdieeJ+xrz+MUqmgQwcjrq7OFz+ZOEedR+zatGnD2bNn/3F/Xl4ebdq0qZeihGgOdi1ayg+Z4QDceH0Yfn5agoJktE4IIWrDaHQjPt6bG27tyrNj26PAxp8lMbw1fgFpySns2XOWsjIZuaurOgc7u91+3hYOxcXFuLjIKkDROlirrHz4wqfsofoywq23JWA0uslKWCGEqIPAQC3h4XquuKE9k5/qjAIbW0ujeWP8x6TuOsT+/blYLFZHl9ms1PpS7Pjx44Hq6+OTJ08+p+WJ1Wpl69atdOzYsd4LFKIp+m7OEpbmdAXgv/fGExrqQUiIDqVS+tYJIURtKRQKoqM9KSqy0PHqeF5SKXll9g62l0Xz+tNLeGruSJydlbRt641GU+fZY61Srf8v7d69G6gesdu3bx9q9f+fR6RWq+nQoQNPP/10/VcoRBNjzjMx6eWNlBFGpLGKe4d2Ra1WYTT+s7+jEEKIC3Nzc6ZzZz8OHMjFdkUsLz2nYMrM7ewsi2DOU4t55u1HKCurom1bbzw85KrIxSjsdru9Lk948MEHmTt3Lnp9y+3TZTabMRgMmEymFv06xaV56eYJvPKDO0psLFs8ELVWR2KikagoT0eXJoQQzZbFYuXgwVxOnDCx/8/9vPJ6MlZUdPbIYuK7D+Md4E3Hjr6tcspLXXJJnefYLVq0SMKOaLUOb9rF+z9Ud0cf2FOPq0GPr68bgYGyaEIIIS6HWq0iIcFIbKwXHa5pzwuPxqCiil2F/rz72FsUnsliz54cCgulifGF1DnYCdFaVVXZmDJ0LjkY0TlbuHN0f7y8XOnQwRc3N1mSL4QQl8vJqXo+Xfv2Rq64qQdPDqvuj5tUEMFHY2eRd/IMe/acpajI4uBKmy4JdkLU0oav/uTr1EAAhg9NwGBwIT7eW/rWCSFEPVIqFUREeNCxoy+97r6O4XeFA/BDYSLvPT6Xs6np7NmTQ0mJhLvzkSUmQtRCVlYJrzyzmnL0hLoX07N/J8LDDXh7uzq6NCGEaJGCgqo7DdgeuAZntRMfLD3Gr+a2lI59nyfefBiVSkFiolG2cPwbGbET4iJsNjvLP9nOxlM6AB7+bwwGgwthYTLXVAghGlJAgJboaA/6Du7Gs2PiUWBjS3Eks55YROq+NPbule3H/k6CnRAXceqUienTkrCj4GrdMeJ6X0tEhAGtVv5KFEKIhhYebiAgwJ0efTrxyoROqLCSXBrK1Cc+5dju43JZ9m8k2AlxEeNHryS/RIEPeYy4Pw5vb1eCg3WOLksIIVoFZ2cVbdt64+XlQrsr2zFrclfUikoOlwfw4pNfcHj7EbZuzSQ93YzVanN0uQ7n0GA3f/58EhMT0ev16PV6evTowY8//ljzeHl5OWPGjMHb2xutVsuQIUPIzs6+4DntdjsvvvgiAQEBuLq60rt3b44ePdrQL0W0UKf+3MKPa1MBuC94P/prbyUiwkM6oAshRCPy8HChSxc/QkN1tOkYx+tTe+CmKCfNYuSFp7/h9MHj7NyZzeHD+dhsdWrP2+I4NNgFBwczc+ZMdu7cyY4dO7jxxhsZNGgQBw4cAODJJ5/k22+/ZeXKlWzYsIGMjAxuv/32C55z9uzZzJs3j/fff5+tW7fi7u5Ov379KC+XvjeibiwlJczp/wjlaNCqKrh25msEhvlIzzohhHAANzdnEhONREd7EhQfw5uzrkevLCGz0pOnx31HaWY6R48WcORI6w53dd55oqF5eXnx2muvcccdd2A0Glm+fDl33HEHAIcPH6Zt27Zs3ryZK6+88h/PtdvtBAYG8tRTT9Vsb2YymfDz8+OTTz7hnnvuqVUNsvOEAEhZuZKH7/qI37mKfr3DeWL8NVxxRYCshBVCCAey2eycOFHIoUN55BxL49ln1pJn1eGpLOKNOX0whEcQE+NJbKwXTk4tY8ZZg+480VCsVisrVqygpKSEHj16sHPnTiorK+ndu3fNMXFxcYSGhrJ58+bzniM1NZWsrKxznmMwGOjevfu/PgegoqICs9l8zk2II19+ySFiAGjfMYioKA8JdUII4WBKpYLISA86dfIjOD6Sma8NwNfJRIFNxxNPrSdr22ZSUvLZt+8s5eVVji630Tk82O3btw+tVotGo2HUqFGsXr2a+Ph4srKyUKvVeHh4nHO8n58fWVlZ5z3XX/f7+fnV+jkAM2bMwGAw1NxCQkIu70WJZq+yrIyt3/7BWXxQKeHKK4Px95dLsEII0RQoFAqCg3V06eJHVKdoXp0zkCCNCbPdnadmHmTnks84caKQ3btzWt0uFQ4PdrGxsSQnJ7N161ZGjx7NAw88wMGDBxu1hokTJ2IymWpu6enpjfr9RdNz8qef2F0aBkBCe3/8/LRotbJtmBBCNCU+Pm507OhLaHw4by4cSpRHKWW4MuMbO0snzuVMWi67dmWRn1/m6FIbjcODnVqtJioqii5dujBjxgw6dOjA3Llz8ff3x2KxUFhYeM7x2dnZ+Pv7n/dcf93/95WzF3oOgEajqVmZ+9dNtG57PlvFRqrncfbsGYa/vzsKhcLBVQkhhPg7Dw8X2rb1RqPV8t4Xo+jfAewo+PqoH1NGzufk3mPs2JHNsWMFrWL0zuHB7u9sNhsVFRV06dIFZ2dn1q9fX/NYSkoKp06dokePHud9bps2bfD39z/nOWazma1bt/7rc4T4O6vFwopvTlKEDl8vZ66/vg0Gg8bRZQkhhPgXQUFaYmI8KS218vTrw3hmaBDOVLG/xJ/x477jwE+/s29fLlu2ZHDmTJGjy21QDg12EydO5PfffyctLY19+/YxceJEkpKSuO+++zAYDIwYMYLx48fz22+/sXPnTh588EF69OhxzorYuLg4Vq9eDVRfcx83bhzTpk1jzZo17Nu3j6FDhxIYGMhtt93moFcpmpsjP6zjl/IuANxzXyc8PTUS7IQQoglTKBTExnrRsaORykob1955A/PfuBZv52LybAYmvH6EDR+vxGazkZycQ2pqIU2sKUi9cWiX1ZycHIYOHUpmZiYGg4HExETWrVtHnz59AHjzzTdRKpUMGTKEiooK+vXrx3vvvXfOOVJSUjCZTDVfT5gwgZKSEh5++GEKCwvp2bMna9euxcXFpVFfm2i+Fs/6mkJCMbhY6XlNBD4+bqjVKkeXJYQQ4gIUCgVhYdXbPR45kk9FRSgfLx/KS6MXsSfXwNzvK9l/8G3Gv/Uw+/blYrXaiYz0aHHTbJpcH7umQPrYtV4VZjNXeo4k2RbP4D5+jHiiD1deGShtToQQohmxWKykppo4dqwAq9XG8umf8sU2JXYUhLsV8uq8u3Dz8SEuzpvoaE+UyqYd7pplHzshmoLk5Ss5YIsCoHuvRMLDDRLqhBCimVGrVcTGetGtmz8Ggwt3PXs/L48MxJUy0ko9eOThr0nbdYBDh/I4fDivRe0xK8FOiP9j+ds/UYkao9ZGly4BtGljcHRJQgghLpGvrztXXOFPmzYGEvpfz4K51xHolE+RzY2nX9nJppW/kJJSwL59Z7FYrI4ut15IsBPif4rOnOG3/7VQvPrqUAIDdbi5Se86IYRozlxdnYmP9yYoSIt7aBSLvhhGB0M2VTjx+vIcPn11OUePFrBrVzaZmcVUVTXv0TsJdkL8z46FyzhENADX9muPj49cghVCiJZAo3GiXTsf/PzcKcGdOSvGcUtkLgCrd8DMMe9z8mgm27Zlsnt3NqWllQ6u+NJJsBPifz77YANVOONvgPbtffH0lJXUQgjRUmi1arp29SMhwYcKq4pRc5/giQEKVFSx66wnT41YQlHqUdLTi9i1K5uCgnJHl3xJJNgJAZzdu5c/zlSvNOp5fRSBgVqcnaXFiRBCtCTOziqio6sXVWi1aq4beS+vTWiLTlHKmUovHn7yd/av/pa8vDK2bcvkxInCZrewQoKdEMDWj5aSQvVq2GtujJGVsEII0YL5+rrTtas/oaE6wrp35f0PBhLhbqIMV15eWsTiR6dQcTKFPXtySE7OoaSk+WxFJsFOtHp2m40vlmzDihPBPiratZPLsEII0dK5u6vp2NGXxEQjHkEBzFk8kl5tq7CjYHVGNI8+8jUHV3zGyZNmtm/PIiurxNEl14oEO9HqpW/YwA5TEABX3xCLv7+77DQhhBCtgEqlJCrKk27d/NF7uPHU6w/y7GPtcVVVcYYAJn0B7zzxNjmn89ixI5MdO7I4caKQsrKmu7hCgp1o9XYt+ZwjRABwVc9wWQ0rhBCtjNHoRkKCDwC9B3Xmmx8fpE/76t0ofjvtwyPDVrB7w15yckrYvTubQ4fym2xbFAl2olWz2+18+/UBrDgR4O1EfLyPXIYVQohWyM/PnXbtfKiqslNWZmPqgoeY83QkPkoThVY3Xpyzj2nPrMRms3LypJnjxwtpiruySrATrVpOcjK7Cv0A6H51BP7+OrkMK4QQrVRoqJ7u3QPw83MjK6uU2OuuZNnn99DXeAIlNjYfrOC+ez7n+MFTHD6cx/79uZSXVzm67HNIsBOt2qHVazhEDABXXhUml2GFEKKV8/R0oWtXf7p3D8BodKPM2ZNxC59nzpASIkml0qbkxcm/kL72W1JS8klJyXd0yeeQYCdatW8/20oFGjzcoVOnADw9NY4uSQghhIOpVEr8/atbonTt6o+bzo2w/zzM7LcGEqHJoczuwgvvZ/DThOcxZWQ5utxzSLATrVZJdjZ/HKu+7Nr9imB8fd3QaJwcXJUQQoimQqVSEhSk48orA0hMNKKOTGTGJ48Q7GnHjI4FqYmMum0up/7c4uhSa0iwE63W0W+/4wBxAHS/Jhqj0c3BFQkhhGiKNBonoqI86dDBF61ey/wl9zN0SBgqrLhb8vjyhuvI3L7d0WUCIMMTotX6btGvFBODm7ONbt0C8fCQy7BCCCH+XUiIHoVCwZ49OQwdfQM9r4+mdHkyQfob8Ovc2dHlARLsRCtVWVbG+m2FAHRO9MbX1w13d7VjixJCCNHkBQfrsFpt7N+fi97oTfs35tO5vQdKVdPoqCDBTrRKaT//zN6q6r1hr7whHn9/rYMrEkII0VyEhupxc3MmLc2ERuOE2t3d0SXVkGAnWqXvPviefAJxUtro3j0Yg0EuwwohhKgdhUKB0eiGj48rNlvTalIsiydEq2OzWvlqfR4A3RN0+Plp0enkMqwQQoi6USgUqFRNK0o1rWqEaATHf9vEjopoAPrc1oWAAHeUSoWDqxJCCCEunwQ70ep8OPsHynHBx7WCxI7BshpWCCFEiyHBTrQqdrudb343A9D7Ck88PDQyv04IIUSLIcFOtCr7ftvF0Qp/AK67/SpCQ3U4OzeNJepCCCHE5ZJgJ1qVD2d9D0CUtpCIuCD8/ZvOEnUhhBDickmwE63K2k3/Ww3bwZPgYK00JRZCCNGiSLATrcaBbUc4XuqFAjs33tmTwEBpSiyEEKJlkWAnWo0PX/0GgEjXPOK6xaLXy6IJIYQQLYsEO9FqrEvKBKBbvFZ61wkhhGiRJNiJViE3s4CjJh0APQd1x9PTxcEVCSGEEPVPgp1oFT6f/zM2lPgqC2h/fRe5DCuEEKJFkmAnWoXv1hwCoJ1vGYGBWrkMK4QQokWSYCdaPLvdzpZDFgC6JHrj4SGXYYUQQrRMEuxEi7drVyaFFjXOWOjetxN6vfSuE0II0TJJsBMt3ldLtwIQxUkir7lCthATQgjRYkmwEy3ez2tTAGjrU4pvkJeDqxFCCCEajgQ70aJVVFSx52gZAB0SvGU1rBBCiBZNgp1o0TZuPInFqkRLMe2v7YBW6+zokoQQQogGI8FOtGhrfzgCQDQniLjhGhQKaXMihBCi5ZJgJ1q09T9W96+Lcc3FLy7SwdUIIYQQDUuCnWixSkos7E0pAiAxTifz64QQQrR4EuxEi3XqlAmbHVwoJ6prPG5uMr9OCCFEy+bQYDdjxgy6deuGTqfD19eX2267jZSUlJrH09LSUCgU572tXLnyX887bNiwfxzfv3//xnhJognJyKgerdNTRNCV3RxcjRBCCNHwHBrsNmzYwJgxY9iyZQs///wzlZWV9O3bl5KSEgBCQkLIzMw85zZlyhS0Wi0DBgy44Ln79+9/zvM+++yzxnhJoglJO5oNVAe7kB5XOLgaIYQQouE5OfKbr1279pyvP/nkE3x9fdm5cyfXXnstKpUKf3//c45ZvXo1d911F1qt9oLn1mg0/3iuaF2O7Kge/fV0qcIYHujgaoQQQoiG16Tm2JlMJgC8vM6/O8DOnTtJTk5mxIgRFz1XUlISvr6+xMbGMnr0aPLy8v712IqKCsxm8zk30fydOJgOgI9Ri1ot24gJIYRo+ZpMsLPZbIwbN46rr76ahISE8x7z8ccf07ZtW6666qoLnqt///4sWbKE9evXM2vWLDZs2MCAAQOwWq3nPX7GjBkYDIaaW0hIyGW/HuF4p08WAOAX5odK1WR+1YUQQogG49BLsf/XmDFj2L9/P5s2bTrv42VlZSxfvpzJkydf9Fz33HNPzX+3b9+exMREIiMjSUpKolevXv84fuLEiYwfP77ma7PZLOGuBcjOrQAgODbUwZUIIYQQjaNJDGM89thjfPfdd/z2228EBwef95hVq1ZRWlrK0KFD63z+iIgIfHx8OHbs2Hkf12g06PX6c26ieSs9e5YCS3XfupgeiQ6uRgghhGgcDh2xs9vtjB07ltWrV5OUlESbNm3+9diPP/6YgQMHYjQa6/x9Tp8+TV5eHgEBAZdTrmhGMrZtx0z1ApvItvJzF0II0To4dMRuzJgxLF26lOXLl6PT6cjKyiIrK4uysrJzjjt27Bi///47Dz300HnPExcXx+rVqwEoLi7mmWeeYcuWLaSlpbF+/XoGDRpEVFQU/fr1a/DXJJqGtH0nsKFCgZ3AQBmBFUII0To4NNjNnz8fk8nE9ddfT0BAQM3t888/P+e4hQsXEhwcTN++fc97npSUlJoVtSqVir179zJw4EBiYmIYMWIEXbp0YePGjWg0sqVUa3H6VPXCCZ26Cje3JjOVVAghhGhQDr8UWxuvvvoqr776aq3O4+rqyrp16y67NtG8nTlTBLjj6WbH2VlanQghhGgdmsTiCSHqW1ZO9eV8T50SZ2f5NRdCCNE6yCeeaJGy8iwA+HioJdgJIYRoNeQTT7RIuSYbAD4+bnIpVgghRKshwU60SLkl1b/avv46lEqFg6sRQgghGocEO9EiFZSrAfAP9nRwJUIIIUTjkWAnWhy73U6Rtbq1TUiEn4OrEUIIIRqPBDvR4lSYzZTgBkBotOw6IYQQovWQYOcgVquNjIxizOaKWvfzE7WTfTIbK9ULJoLb1H0LOiGEEKK5kpb8DlJcXMm+fbnYbDYCAtyJjfXC1dXZ0WW1CGeOZwKgUVjQ6dQOrkYIIYRoPDJi5wA2m407uk5izbzP0WiUnDhhYufObPLzyy7+ZHFRGSdzAdCqKnBykl9xIYQQrYd86jnAwpeX8csxLe/9UM5jd75DVXEhBQXlbN+ezcmTJmw2uTR7ObJO/2+fWGerNCcWQgjRqsinngPc/9xdjLjOhhOVHMrXMXTYN2xc/DUqpY3du3PYsyeHkhKLo8tstnIyzQDoXGwyYieEEKJVkU89B3Bx0/D6N8/z1vMJRGmyqLCree3LEp66fTYFe3eSlmZm27YssrNLHF1qs5STUwyAwU0pwU4IIUSrIp96DhR+VTc+/PZJ7uhmQ4mVvcUBjHzpIAuefIPsY6fYsSOLw4fzyM8vk8uzdZCbVw6Ah95Jgp0QQohWRT71HMTd3ZmgIC0FRfDEGyNZuqA/iX4lWFHx40l/Ro/6mr3frefw4Xz+/PMMe/fK5dnayiusBMDDw8XBlQghhBCNS4Kdgzg5KUlIMBIb60VeXhneoQEs/H48Mye0w+hcRKFNxwvvn+GTZ99GXZZHaqqZLVsyycgolr53F1FQZAXAy9vNwZUIIYQQjUuCnQM5OSlp29abLl38UakUnDxppkufrqz4diTXhRVhR8k3R7wZeffHHFvyPpaSEnbsyOLgwTwsFqujy2+yCksVAHgbdQ6uRAghhGhcEuwcTKlUEByso3v3QBITjTg5KTGVKpn4waNMeiQcV1Ulp+xBPPu1lhn/mUrJgW2kpOSze3c2ZnOFo8tvkkyW6r7bgaHeDq5ECCGEaFwS7JoIrVZNVJQnV10VRNeufri4ONGpbw8Wf/Yfrk10xY6SjSVxjHh+L2vnLCAjvYBt2zJJTzdjtdocXX6TYbfbKa7SABAULtuJCSGEaF0k2DUxarWKkBA93bsHEh/vjVrtzNS37+W9t3oRqi+nDFcWbNYyZfjbnN6RzK5d2WzZkkF6upmqKgl4ZnMFVf/bJzY0JtDB1QghhBCNS4JdE+Xi4kRsrFfN4or4jqGs/HksD92kQ0UVe4oCGD1xK1+9NJf8jFx27Mhi584sCgvLHV26Q505kQOAGgveIX4OrkYIIYRoXBLsmjCFQkF0tCexsV6UlVk5c6aY/zw9hI/evpFwfSnluLBkj5FJQ+ehOH2ArKxStm7N4sSJwlY7encmNQsAd0px0bk7uBohhBCicUmwa+JUKiXx8T5cdVUgbdt6U1lpxxAUxCdrxjDm7gCcsLK3PIzhj2/g93c/RqWoYs+eHPbvP0tlZetbOZt5MhcAnZNF9okVQgjR6sgnXzOh12uIi/PmqqsCiYrypKjIQr/7buTdt/vi7WKhAE/e+NmFR2+fR86RE5w4YWL//lwqKqocXXqjykzPB0DnXCm7TgghhGh15JOvmdFq1SQk+NC9eyC+vm74BvuxeNWD3Hu9K85UcrTEm9HPbWHVnM84fjSP3btzKCpqPTtW5GSZAdC72CXYCSGEaHXkk68ZUigUGI1udO3qT4cORqqqYNSUe/nsoz5c4ZmBHQWfbYZ5o1/jxObdbNuWQVqaqVU0Nc7OLgZA765EpZJfbyGEEK2LfPI1YyqVkvBwA7GxnuTllaELCWXeD5MZ2ccZBXY25ofz3FNr2DP/XXZuO822bZlkZhZjs7XcLcnO5pYB4Kl3cnAlQgghROOTYNfM/bVytlMnPzQaFRmZpQx+8h6mv9AVN6dK0u2BTP5SydanR5G9ayfbtmWxf38u5eUtc+7d2fzqy85eni4OrkQIIYRofBLsWoC/Ru6uuiqIbt0C0OnUxHaOYfFn9xEd5EQpbsxNv5bPn5xG+W8rOHo0nx07ssjMLG5xbVFyTdWXm/2DPR1ciRBCCNH4JNi1IM7Oqv/tOxtAbKwXzs7OvPbunXTr7IsFNStttzBqbh6/T3uVnP0H2bYtkz//zODYsQLy8spaxCheYVn1r3RwZICDKxFCCCEan0xEaoGcnVXEx3tjMGhIScnn6WdvYOPGVJYu3EphqQfzkj34bfTH9G1TSKd+V5LTpTfOxgDc3JwIDNQSGqpHp1M7+mXUWVWVjeKq6rrD40MdXI0QQgjR+BR2u73lzqS/RGazGYPBgMlkQq/XO7qcy1JcbGHfvrOcPVuGt7cL70z/js9/zMKOAgAtxUSRysArXbh6/OOUOHng7q4mJsaTkBBds1pZevpUASFh81Bg42jyf4nsEOXokoQQQojLVpdc0nw+tcUl0WrVtG3rjaurE2azhSdfHsjSZXcwsH8YrmooRksy7XllSzST7p3Nma+XoaSK3btzSE7OaVZ7z6YdSANASykeQXIpVgghROsjwa4V8PBwISHBBxcXJzIzSzAYXHn+5b6sT3qIDz+8lX49fQDYbm3P2CVKxg95nbMHD3D6dBFbt2Zy9Gh+s9jBIu1QOgA6pwrUGmcHVyOEEEI0Pplj10oEBGjx8XElL6+c1NRCMjNLUCoVRER4M3nmIO48fJYPZ//A9iOV7Cv2Z9Sk3Vwb9Csjn72FkpJwzpwpxmh0Q6dT4+vrhotL0/vVOX0sCwCDiw0nJ4WDqxFCCCEaX9P7dBYNxtlZhb+/Oz4+ruTklJKXV8bZs6WYzRY8Pd2Zt+RBTh44wZz/1969B0dV3/8ff+5usrnuJiSQGwlJuAVQwt18I4JYEARFUy+DlClBRJk2FPKj1EoLBaojMP06hValDF9/iTe+1B8/QylVFLEkX+RiCI0CCsr9lhvXbO7J7vn+Ed0SCWgwZJPl9ZjZGc4n53zOe/dMJi8+ez6fM389+aWdyD0Tzs7ZuYzt+Q8e/z9plMfG4XQa2O1WevQIpWvXYHx9LZ5+W25nTp4HIDTI3K7qEhERaSsKdrcgHx8zMTHBxMQE43S6qKys59ChC5w65aBrUgKr3n2WHRs/ZsUfd3DUYWfT4TDez9jKPYkOZsybgGFLYu/eEk6fdpCYGEJERGC7CFIlReWAP506+es5sSIickvSX79bnMVixm73Izm5CwkJIVy+XMfJk+X0uWsI//3hL3nuF91JDLpMPb58cCyMyRkf83z6f1J96igORx35+cXs2HGGI0cucvlyLRUVdR57Jm1JWeNEj85dgj1yfhEREU/TiJ0A4Ofnw6BBETgcdVy4UMOxY5c5fbqCoROGM+YnPyLv/+fxf/+rgEOXbOSe7UJe5k6GhW9k4v09GZQ2lk/P1+DnZ8FkMuHvbyE+3k54eABBQb5tNpp3/nLjBI/O0aFtcj4REZH2RsFO3EwmE3a7H3a7H9HRQZSUVHHyZDklJdX0GTGM7EdGkv/+Htb8aRv7LoTyyflIPnndgf31bFLjaxj/8CCGpI2hps7g00/LsFrNBAb6EhMTTEREIHa7H2azidraBnx8Wv8+uItVjQPQMYmRrdqviIhIR6EFipvhTQsU/1BOp4tz56o5fdrBiRPlREY2zoj9NPdT/nvN/7DjSydVLn/3/hGWS4zoa+LuYeEkDuqFf9JQKqqcWCwmgoJ8sVhM1NQ0zlq126106uRPUJAvNpuVoCBfTKbG2awul4HZ/P1ntroaGrD7/oZKgvh/2WN5ND211T8LERERT2hJLlGwa4aC3dXq650cPHiBo0cvfT2yZyU42Jf6eicf/TWPf7xTyJ5TVur59/pxJlxE+VziP24PYsgd8fQY0JPOvRLx8/elocFFdXUD9fUuDAP8/X2w26107hxAeXkdDkcdnTr5ERYWQGCgD8HB1ususXLh2Ak6d8/CwMSe/BkMGdq1LT4WERGRm07B7gdSsGue0+mipKSKsrLGV0VFvXskLijIl6ryCv6x5h98uL2UAyVWGlxXz82x4yA2sJLIMB+iY+zEdo8grm88Xfv1JCA4gJoaF76+Zvz9Le7gZzabCQz0ITTUj/DwAHx8zJSUVALQrVvjvXx7N35EyiMfY8LF6TPziImxtfXHIyIiclMo2P1ACnbfraamgfPnqzl3rprz56uprKynocHlXrzY5TIoPXuB3Dff43+2n+LoeR9K66//WQZbaukWWk9UZ6s79HXtGUt0j650CgukttZJXZ0LwzCwWi24XAYuF9jtVj5b/zcyX7mAzVLD8dLfERYW0EafhIiIyM3VYYLd0qVLeeeddzh48CABAQHceeedLF++nKSkJPc+o0aNIjc3t8lxM2fO5C9/+cs1+zUMg0WLFrFmzRouXbrE8OHDWbVqFb169fpedSnYtUx9vROHo44zZyo4evQyVquZ4GBffH3NmEwm95pylY4aDuz6gi//dZTTx0opLiqn5EID56p9uWTYMK6z+o4FJxEBNcSGmegcE0an7onckZrA0KExVFfXs/hnWXx82EKCrZJPTz+H3e7XVm9fRETkpuowwe6+++7j8ccfZ9iwYTQ0NPCb3/yG/fv38/nnnxMUFAQ0BrvevXvz+9//3n1cYGDgdd/Y8uXLWbp0Ka+99hqJiYksXLiQffv28fnnn+Pv73/N476hYHdjDMPgzJkKTp4sx+Gop6HBictl0NDgwmxuDHiBgT5XhS7DMCg/eYJ9uYUc3n+KorPllJyro9QBF2utXCSEBpp/9qvN6iQixMSRMjNmXGTedZHnP3iRgAA9K1ZERLxDS3KJR5c72bx5c5Pt7OxsIiIiKCgoYOTIke72wMBAoqKivlefhmGwYsUKFixYwEMPPQTA66+/TmRkJBs2bODxxx9vvTcgTZhMJmJjbXTtGkx1dQMNDS6cToOqqnoqK+upqqqntLSakycd+Pk1js4ZBphMYAmK5PYHxzP4UTMBAT7/nh3bUE910VlO7j/M0c9PcfSrMsoOHeaCw8k++uGoC8JR1nj+R9nI0AEPtounYIiIiHhCu1rH7vLlywCEhYU1aX/rrbd48803iYqKYuLEiSxcuJDAwMBm+zh27BjFxcWMGTPG3RYSEkJKSgo7d+5sNtjV1tZSW1vr3i4vL2+Nt3PLMplMBAb+e8SsU6d/j5I6HHWcPu2guroBk6kx1Llcjffs1dY6qapycu5cDYGBFnx8LJjNYAqNJm5ENIn3mBlrMWE2m6g4sJeLBw9w6OB5Dh+5gP+FY/QxXyJs+L1YLN9/mRQRERFv0m6CncvlIjMzk+HDh3P77be723/yk58QHx9PTEwMn332Gb/+9a85dOgQ77zzTrP9FBcXAxAZ2XSR2sjISPfPvm3p0qUsWbKkld6JXI/NZqVv3/Bmf+ZyGdTUNFBaWsWpUw4aGly4XAZg4HIZVFc7cTpdNDS4MIX3ofOEZLo96ss4s4nq6sbJHH0GRrhH+0RERG417SbYZWRksH//frZv396k/emnn3b/u3///kRHRzN69GiOHDlCjx49WuXc8+fPZ+7cue7t8vJy4uLiWqVv+f7M5saRvoSEEBISQnC5DAzDwDBwBzqn06C6uoGiogrKyqq5dKkCw2j8Cj4pKYxu3XRPpIiI3LraRbCbNWsWmzZtIi8vj9jY2Ovum5KSAsDhw4ebDXbf3ItXUlJCdHS0u72kpISBAwc226efnx9+fppF2d40PnmicfTNx8fMN5fIbvcjMjKImpoGHI466uqcGAZERQVptE5ERG5p115fog0YhsGsWbPIycnho48+IjEx8TuPKSwsBGgS2q6UmJhIVFQUW7dudbeVl5eze/duUlP1mClv4u/vQ5cugXTtaiM21uZeVkVERORW5dG/hBkZGbz55pusXbsWm81GcXExxcXFVFdXA3DkyBGee+45CgoKOH78OBs3bmTq1KmMHDmS5ORkdz99+vQhJycHaLxxPzMzk+eff56NGzeyb98+pk6dSkxMDGlpaZ54myIiIiJtwqNfxa5atQpoXKvuSllZWUybNg2r1cqHH37IihUrqKysJC4ujkceeYQFCxY02f/QoUPuGbUAzzzzDJWVlTz99NNcunSJu+66i82bN3+vNexEREREOio9UqwZWqBYRERE2ouW5BLdlCQiIiLiJRTsRERERLyEgp2IiIiIl1CwExEREfESCnYiIiIiXkLBTkRERMRLKNiJiIiIeAkFOxEREREvoWAnIiIi4iUU7ERERES8hIKdiIiIiJdQsBMRERHxEj6eLqA9MgwDaHzoroiIiIgnfZNHvskn16Ng1wyHwwFAXFychysRERERaeRwOAgJCbnuPibj+8S/W4zL5eLs2bPYbDZMJtNNOUd5eTlxcXGcOnUKu91+U84hN07Xp33T9WnfdH3aN12f9uta18YwDBwOBzExMZjN17+LTiN2zTCbzcTGxrbJuex2u36x2jFdn/ZN16d90/Vp33R92q/mrs13jdR9Q5MnRERERLyEgp2IiIiIl1Cw8xA/Pz8WLVqEn5+fp0uRZuj6tG+6Pu2brk/7puvTfrXGtdHkCREREREvoRE7ERERES+hYCciIiLiJRTsRERERLyEgp2IiIiIl1Cw84CXX36ZhIQE/P39SUlJ4ZNPPvF0SfK1vLw8Jk6cSExMDCaTiQ0bNni6JPna0qVLGTZsGDabjYiICNLS0jh06JCny5KvrVq1iuTkZPfCqqmpqbz33nueLkuuYdmyZZhMJjIzMz1digCLFy/GZDI1efXp0+eG+lKwa2N//etfmTt3LosWLWLv3r0MGDCAcePGUVpa6unSBKisrGTAgAG8/PLLni5FviU3N5eMjAx27drFli1bqK+vZ+zYsVRWVnq6NAFiY2NZtmwZBQUF7Nmzhx/96Ec89NBDHDhwwNOlybfk5+ezevVqkpOTPV2KXOG2226jqKjI/dq+ffsN9aPlTtpYSkoKw4YN46WXXgIan0sbFxfHL37xC5599lkPVydXMplM5OTkkJaW5ulSpBllZWVERESQm5vLyJEjPV2ONCMsLIw//OEPPPnkk54uRb5WUVHB4MGDeeWVV3j++ecZOHAgK1as8HRZt7zFixezYcMGCgsLf3BfGrFrQ3V1dRQUFDBmzBh3m9lsZsyYMezcudODlYl0PJcvXwYaw4O0L06nk3Xr1lFZWUlqaqqny5ErZGRkcP/99zf5OyTtw1dffUVMTAzdu3dnypQpnDx58ob68WnluuQ6zp07h9PpJDIyskl7ZGQkBw8e9FBVIh2Py+UiMzOT4cOHc/vtt3u6HPnavn37SE1NpaamhuDgYHJycujXr5+ny5KvrVu3jr1795Kfn+/pUuRbUlJSyM7OJikpiaKiIpYsWcKIESPYv38/NputRX0p2IlIh5ORkcH+/ftv+B4UuTmSkpIoLCzk8uXLrF+/nvT0dHJzcxXu2oFTp04xZ84ctmzZgr+/v6fLkW8ZP368+9/JycmkpKQQHx/P22+/3eJbGRTs2lDnzp2xWCyUlJQ0aS8pKSEqKspDVYl0LLNmzWLTpk3k5eURGxvr6XLkClarlZ49ewIwZMgQ8vPzWblyJatXr/ZwZVJQUEBpaSmDBw92tzmdTvLy8njppZeora3FYrF4sEK5UmhoKL179+bw4cMtPlb32LUhq9XKkCFD2Lp1q7vN5XKxdetW3Yci8h0Mw2DWrFnk5OTw0UcfkZiY6OmS5Du4XC5qa2s9XYYAo0ePZt++fRQWFrpfQ4cOZcqUKRQWFirUtTMVFRUcOXKE6OjoFh+rEbs2NnfuXNLT0xk6dCh33HEHK1asoLKykieeeMLTpQmNv0xX/g/p2LFjFBYWEhYWRrdu3TxYmWRkZLB27Vr+9re/YbPZKC4uBiAkJISAgAAPVyfz589n/PjxdOvWDYfDwdq1a9m2bRvvv/++p0sTwGazXXU/alBQEOHh4bpPtR2YN28eEydOJD4+nrNnz7Jo0SIsFguTJ09ucV8Kdm1s0qRJlJWV8bvf/Y7i4mIGDhzI5s2br5pQIZ6xZ88e7rnnHvf23LlzAUhPTyc7O9tDVQk0LoALMGrUqCbtWVlZTJs2re0LkiZKS0uZOnUqRUVFhISEkJyczPvvv8+9997r6dJE2r3Tp08zefJkzp8/T5cuXbjrrrvYtWsXXbp0aXFfWsdORERExEvoHjsRERERL6FgJyIiIuIlFOxEREREvISCnYiIiIiXULATERER8RIKdiIiIiJeQsFORERExEso2ImIiIh4CQU7EfFK27Ztw2QycenSpTY/t8lkwmQyERoaet39Fi9ezMCBA9ukpm/O901tK1asaLPzikjbUbATkQ5v1KhRZGZmNmm788473Y+38oSsrCy+/PJLj5z7WubNm0dRURGxsbGeLkVEbhI9K1ZEvJLVaiUqKspj5w8NDSUiIsJj529OcHAwwcHBWCwWT5ciIjeJRuxEpEObNm0aubm5rFy50v014/Hjx6/6KjY7O5vQ0FA2bdpEUlISgYGBPProo1RVVfHaa6+RkJBAp06dmD17Nk6n091/bW0t8+bNo2vXrgQFBZGSksK2bdtuqNZly5YRGRmJzWbjySefpKampsnP8/Pzuffee+ncuTMhISHcfffd7N271/3z6dOn88ADDzQ5pr6+noiICF599VUA1q9fT//+/QkICCA8PJwxY8ZQWVl5Q/WKSMejYCciHdrKlStJTU3lqaeeoqioiKKiIuLi4prdt6qqij/96U+sW7eOzZs3s23bNn784x/z7rvv8u677/LGG2+wevVq1q9f7z5m1qxZ7Ny5k3Xr1vHZZ5/x2GOPcd999/HVV1+1qM63336bxYsX88ILL7Bnzx6io6N55ZVXmuzjcDhIT09n+/bt7Nq1i169ejFhwgQcDgcAM2bMYPPmzRQVFbmP2bRpE1VVVUyaNImioiImT57M9OnT+eKLL9i2bRsPP/wwhmG0qFYR6cAMEZEO7u677zbmzJnTpO2f//ynARgXL140DMMwsrKyDMA4fPiwe5+ZM2cagYGBhsPhcLeNGzfOmDlzpmEYhnHixAnDYrEYZ86cadL36NGjjfnz51+zHsDIyclp0paammr8/Oc/b9KWkpJiDBgw4Jr9OJ1Ow2azGX//+9/dbf369TOWL1/u3p44caIxbdo0wzAMo6CgwACM48ePX7NPwzCM+Ph4449//ON19xGRjkkjdiJyywgMDKRHjx7u7cjISBISEggODm7SVlpaCsC+fftwOp307t3bfX9acHAwubm5HDlypEXn/uKLL0hJSWnSlpqa2mS7pKSEp556il69ehESEoLdbqeiooKTJ0+695kxYwZZWVnu/d977z2mT58OwIABAxg9ejT9+/fnscceY82aNVy8eLFFdYpIx6bJEyJyy/D19W2ybTKZmm1zuVwAVFRUYLFYKCgouGrCwZVhsLWkp6dz/vx5Vq5cSXx8PH5+fqSmplJXV+feZ+rUqTz77LPs3LmTHTt2kJiYyIgRIwCwWCxs2bKFHTt28MEHH/DnP/+Z3/72t+zevZvExMRWr1dE2h+N2IlIh2e1WptMeGgtgwYNwul0UlpaSs+ePZu8Wjrjtm/fvuzevbtJ265du5psf/zxx8yePZsJEyZw22234efnx7lz55rsEx4eTlpaGllZWWRnZ/PEE080+bnJZGL48OEsWbKEf/3rX1itVnJyclpUq4h0XBqxE5EOLyEhgd27d3P8+HGCg4MJCwtrlX579+7NlClTmDp1Ki+++CKDBg2irKyMrVu3kpyczP333/+9+5ozZw7Tpk1j6NChDB8+nLfeeosDBw7QvXt39z69evXijTfeYOjQoZSXl/OrX/2KgICAq/qaMWMGDzzwAE6nk/T0dHf77t272bp1K2PHjiUiIoLdu3dTVlZG3759f9gHISIdhkbsRKTDmzdvHhaLhX79+tGlS5cm96T9UFlZWUydOpVf/vKXJCUlkZaWRn5+Pt26dWtRP5MmTWLhwoU888wzDBkyhBMnTvCzn/2syT6vvvoqFy9eZPDgwfz0pz9l9uzZza6FN2bMGKKjoxk3bhwxMTHudrvdTl5eHhMmTKB3794sWLCAF198kfHjx9/YmxeRDsdkGJoHLyLSmkwmEzk5OaSlpd2U/isqKujatStZWVk8/PDDLT4+ISGBzMzMq57WISIdn0bsRERugsmTJ7f6o7tcLhelpaU899xzhIaG8uCDD7bo+BdeeIHg4OBWHdEUkfZFI3YiIq3s8OHDQOMs1dacjXr8+HESExOJjY0lOzub0aNHt+j4CxcucOHCBQC6dOnisefoisjNo2AnIiIi4iX0VayIiIiIl1CwExEREfESCnYiIiIiXkLBTkRERMRLKNiJiIiIeAkFOxEREREvoWAnIiIi4iUU7ERERES8xP8C3Fza4/gpwfcAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(df.index, df[\"T_int\"], color=\"darkred\", label=\"data\")\n", "plt.plot(tnew, ym, color=\"navy\")\n", "plt.fill_between(\n", " tnew, ym - 2 * ysd, ym + 2 * ysd, color=\"darkblue\", alpha=0.2, label=r\"95% CI\"\n", ")\n", "plt.xlabel(\"time [days]\")\n", "plt.ylabel(\"temperature [°C]\")\n", "plt.tight_layout()\n", "plt.legend(loc=\"best\", fancybox=True, framealpha=0.5)" ] } ], "metadata": { "jupytext": { "cell_metadata_filter": "-all", "main_language": "python", "notebook_metadata_filter": "-all" }, "kernelspec": { "display_name": "pysip-QZQT68sk-py3.11", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.3" } }, "nbformat": 4, "nbformat_minor": 5 }