Created
August 9, 2025 13:05
-
-
Save julian-klode/320dbeff578897b62be6a09d82ee0880 to your computer and use it in GitHub Desktop.
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| { | |
| "cells": [ | |
| { | |
| "cell_type": "code", | |
| "execution_count": 1, | |
| "id": "01cbd225", | |
| "metadata": {}, | |
| "outputs": [], | |
| "source": [ | |
| "import numpy as np\n", | |
| "import pandas as pd\n", | |
| "import matplotlib.pyplot as plt\n", | |
| "from scipy.interpolate import UnivariateSpline\n", | |
| "\n", | |
| "# Extracted from .shot file by hand\n", | |
| "raw_stamps = pd.to_datetime(list(map(float, \"-8.673 -8.480 -8.223 -7.585 -7.360 -7.173 -6.910 -6.721 -6.047 -5.784 -5.598 -5.410 -5.145 -4.510 -4.283 -4.059 -3.834 -3.648 -2.333 -2.073 -1.886 -1.622 0.154 0.327 0.568 0.741 0.964 1.189 1.415 1.640 1.865 2.282 2.546 2.727 2.952 3.186 3.413 3.628 3.815 4.039 4.264 4.490 4.715 4.911 5.145 5.360 5.577 5.820 6.064 6.253 6.477 6.664 6.900 7.115 7.344 7.565 7.790 8.014 8.202 8.428 8.652 8.876 9.102 9.327 9.580 9.783 9.965 10.189 10.414 10.639 10.877 11.114 11.315 11.502 11.727 11.953 12.178 12.406 12.664 12.852 13.078 13.311 13.529 13.715 13.941 14.202 14.427 14.613 14.839 15.064 15.262 15.477 15.702 15.926 16.153 16.380 16.602 16.837 17.027 17.240 17.465 17.728 17.940 18.140 18.364 18.576 18.806 19.001 19.263 19.460 19.677 19.901 20.127 20.324 20.552 20.779 20.989 21.214 21.439 21.686 21.852 22.076 22.304 22.527 22.752 22.993 23.249 23.464 23.727 23.914 24.140 24.364 24.590 24.776 25.043 25.226 25.489 25.676 25.902 26.126 26.353 26.550 26.803 27.020 27.213 27.446 27.664 27.891\".split())), unit=\"s\")\n", | |
| "raw_weights = pd.Series(map(float, \"0.3 0.2 0.0 0.0 0.4 0.4 0.0 0.0 0.0 0.0 0.2 0.6 0.4 0.0 0.0 0.5 1.0 0.7 0.0 0.0 0.4 0.0 0.2 0.0 0.0 0.0 0.0 0.0 0.3 0.5 0.0 0.0 0.0 0.2 0.3 0.8 0.6 0.3 0.3 0.5 1.2 0.7 0.0 0.0 0.3 0.8 1.2 2.3 2.3 2.1 2.6 3.1 3.2 3.0 3.5 4.2 4.7 5.0 5.0 5.1 5.2 5.3 5.5 7.0 7.6 8.1 8.0 7.9 8.0 8.6 9.6 10.2 10.4 10.6 11.2 12.1 13.1 13.7 13.5 14.1 14.8 15.4 15.4 15.3 15.5 16.8 17.2 17.5 18.3 18.6 18.7 18.9 19.8 19.8 20.0 21.4 22.2 21.9 22.1 22.2 22.4 23.3 24.8 25.2 25.4 25.3 26.0 26.7 26.9 27.0 27.4 28.4 29.4 28.9 29.0 30.0 30.4 31.6 31.3 31.7 32.9 33.2 32.9 32.5 33.1 34.6 35.1 35.6 35.9 36.5 36.6 36.1 37.0 37.0 37.1 36.9 36.7 36.4 36.2 36.1 36.0 36.2 36.4 36.5 36.5 36.4 36.3 36.3\".split()), index=raw_stamps)\n" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 2, | |
| "id": "31af7e09", | |
| "metadata": {}, | |
| "outputs": [ | |
| { | |
| "data": { | |
| "image/png": "iVBORw0KGgoAAAANSUhEUgAAA90AAAHqCAYAAAAZLi26AAAAQHRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjErZGZzZzEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvzRIYmAAAAAlwSFlzAAAPYQAAD2EBqD+naQAApKlJREFUeJzs3XmcjeX/x/HXmTE7M5YZjHUGIfuaNQyRVERFaUGlfd8pIZWllMqvpE3fUlEpkWQLCSGElLIOGWksY5sZM3Pu3x+3OeOYM+Occc6cZd7Px+M85tz33MvnPtd9ls99Xfd1WQzDMBARERERERERtwvydgAiIiIiIiIigUpJt4iIiIiIiIiHKOkWERERERER8RAl3SIiIiIiIiIeoqRbRERERERExEOUdIuIiIiIiIh4iJJuEREREREREQ9R0i0iIiIiIiLiIUq6RURERERERDxESbeISACYNm0aFouFdevWubzuqVOnGDVqFEuXLnV/YB4yatQoLBZLkdadN28eo0aNcm9AZ3Tp0oUuXbqcd7mEhAQsFovtERUVRYsWLZg8eTKGYRRp3ytXrmTUqFEcPXq0SOu7w4YNG+jcuTMxMTFYLBYmTZqUb5mcnBzKli3LFVdcke9/r732GhaLhRtvvDHf/8aMGYPFYmHTpk1Ox7N06VIsFkuRzu3du3djsVh45ZVXzrtsUc+pOXPmcPXVV1OpUiVCQ0MpX7483bp1Y/r06WRlZbkch4iI+CYl3SIiJdypU6cYPXq0XyXdF2LevHmMHj3a22HQoUMHVq1axapVq/j444+JjIzkgQceYOzYsUXa3sqVKxk9erRXk+7bbruNlJQUPv/8c1atWsUNN9yQb5ng4GAuvfRSVqxYQXZ2tt3/li5dSlRUFD/++GO+9ZYuXUqFChVo3Lix0/G0aNGCVatW0aJFC9cPxgWunlOGYTBkyBB69+6N1Wrl1VdfZdGiRXz00Uc0bdqUe++9l7feesuDEYuISHFS0i0iIh6RlZWVL6mSPGXLlqVt27a0bduWfv36MW/ePGJiYnjnnXe8HVqRbdmyhcsuu4wrrriCtm3bUrlyZYfLJSUlceLECbuWGVarlZ9++ol77rmHf//9lz/++MP2v9OnT7Nq1Sq6dOniUguH6Oho2rZtS3R0dNEPygNefvllpk2bxujRo/nuu++46aab6NSpE1dffTWvvvoq27Zto2XLlt4OU0RE3ERJt4hIgBo8eDClS5dm+/bt9OrVi9KlS1O9enUee+wxMjMzAbPpalxcHACjR4+2NXcePHiwbTt///03AwcOpGLFioSFhXHxxRfzf//3f3b7ym3G+/HHH/PYY49RtWpVwsLC2L59u63p+8KFCxkyZAjly5cnKiqKq6++mp07d+aL+4MPPqBp06aEh4dTvnx5+vbta5eAFWTGjBn06NGD+Ph4IiIiuPjii3n66ac5efKk3WuSG/vZzbt3794NmDWQb731Fs2aNSMiIoJy5cpx3XXX5YvTMAwmTJhAzZo1CQ8Pp0WLFnz//ffnL5RCREdHU7duXf7991+7+QsXLqRPnz5Uq1aN8PBw6tSpw1133UVqaqptmVGjRvHEE08AkJiYaDuus1svzJgxg3bt2hEVFUXp0qW5/PLL2bBhg1OxbdmyhT59+lCuXDnCw8Np1qwZH330ke3/uWWcnZ3N22+/bdt/QZKSkgDs4vvtt984cuQId955J/Hx8Xa13b/88gvp6em29QDWrVtH7969KV++POHh4TRv3pyZM2fa7aeg5uXvvvsudevWJSwsjAYNGvDpp58yePBgEhISHMb76quvkpiYSOnSpWnXrh2rV6+2/e9859S5srKyGD9+PPXr12fEiBEOl6lcuTIdO3Z0KY7c1+SGG24gISGBiIgIEhISuPHGG9mzZ4/dcrnl9eOPP3LPPfcQGxtLhQoV6NevH/v377dbNjMzk8cee4zKlSsTGRlJp06d+PXXX0lISLD7nAA4cOAAd911F9WqVSM0NJTExERGjx6ti28iUuIp6RYRCWBZWVn07t2bbt26MXv2bG677TZee+01xo8fD0B8fDzz588H4Pbbb7c1d85NBrZu3Urr1q3ZsmULEydOZO7cuVx55ZU8+OCDDpvTDhs2jOTkZKZMmcKcOXOoWLGi7X+33347QUFBfPrpp0yaNIk1a9bQpUsXu+bQY8eO5fbbb6dhw4bMmjWL119/nU2bNtGuXTv+/vvvQo/177//plevXrz//vvMnz+fhx9+mJkzZ3L11VfblhkxYgTXXXcdgO1YV61aRXx8PAB33XUXDz/8MJdddhnffPMNb731Fr///jvt27e3S4ZHjx7NU089Rffu3fnmm2+45557GDp0KNu2bXOleOxkZ2ezd+9e6tatazd/x44dtGvXjrfffpsFCxbw3HPP8csvv9CxY0fbfb933HEHDzzwAACzZs2yHVdus+qXXnqJG2+8kQYNGjBz5kw+/vhjjh8/zqWXXsrWrVsLjWvbtm20b9+e33//nTfeeINZs2bRoEEDBg8ezIQJEwC48sorWbVqFQDXXXedbf8Fadq0KeXKlbNLrH/88Ufi4+O56KKL6NSpk12inLtcbtL9448/0qFDB44ePcqUKVOYPXs2zZo1Y8CAAUybNq3Q45k6dSp33nknTZo0YdasWTz77LOF3l7xf//3fyxcuJBJkyYxffp0Tp48Sa9evUhLSwPOf06da926dRw+fJg+ffq4VGt/vjjAvIhWr149Jk2axA8//MD48eNJSUmhdevWdhdpct1xxx2EhITw6aefMmHCBJYuXcrNN99st8yQIUOYNGkSQ4YMYfbs2Vx77bX07ds3320MBw4c4JJLLuGHH37gueee4/vvv+f2229n7NixDB061OnjFBEJSIaIiPi9Dz/80ACMtWvX2uYNGjTIAIyZM2faLdurVy+jXr16tun//vvPAIyRI0fm2+7ll19uVKtWzUhLS7Obf//99xvh4eHG4cOHDcMwjB9//NEAjE6dOhUYW9++fe3m//zzzwZgvPDCC4ZhGMaRI0eMiIgIo1evXnbLJScnG2FhYcbAgQNt80aOHGkU9hVmtVqNrKwsY9myZQZg/Pbbb7b/3XfffQ7XXbVqlQEYEydOtJu/d+9eIyIiwnjyySdtcYaHhxd4PJ07dy4wrlw1a9Y0evXqZWRlZRlZWVnGnj17jKFDhxohISHG3Llzz3tce/bsMQBj9uzZtv+9/PLLBmDs2rXLbp3k5GSjVKlSxgMPPGA3//jx40blypWN/v37FxrrDTfcYISFhRnJycl286+44gojMjLSOHr0qG0eYNx3333nO3zDMAzjmmuuMaKiooysrCzDMAzj6quvNm644QbDMAzjrbfeMuLi4gyr1WoYhmEkJSUZFStWtK1bv359o3nz5rZ1c1111VVGfHy8kZOTYxhG3nn5448/GoZhGDk5OUblypWNNm3a2K23Z88eIyQkxKhZs6Zt3q5duwzAaNy4sZGdnW2bv2bNGgMwPvvsM9u8gs4pRz7//HMDMKZMmeLU8q7Eca7s7GzjxIkTRlRUlPH666/b5ue+J++991675SdMmGAARkpKimEYhvH7778bgPHUU0/ZLffZZ58ZgDFo0CDbvLvuussoXbq0sWfPHrtlX3nlFQMwfv/9d6eOV0QkEKmmW0QkgFksFruaXoAmTZrka27qSEZGBosXL6Zv375ERkaSnZ1te/Tq1YuMjIx8zVuvvfbaArd300032U23b9+emjVr2moxV61aRXp6er4mq9WrV6dr164sXry40Hh37tzJwIEDqVy5MsHBwYSEhNC5c2cAp5qnz507F4vFws0332x3rJUrV6Zp06a2mtBVq1aRkZFR4PE4a968eYSEhBASEkLNmjV59913efPNN7nyyivtljt48CB333031atXp1SpUrblnT2uH374gezsbG699Va74woPD6dz587n7UBvyZIldOvWjerVq9vNHzx4MKdOnSq0RrswSUlJnDx5krVr19ru587t+b1z5878999//P7772RmZrJ69WpbLff27dv5888/ba//uedlSkpKgS0Otm3bxoEDB+jfv7/d/Bo1atChQweH61x55ZUEBwfbpps0aQLg1HvInZyJ48SJEzz11FPUqVOHUqVKUapUKUqXLs3Jkycdniu9e/e2mz53m8uWLQPI93pdd911lCpVym7e3LlzSUpKokqVKnZlkttLfe62RERKolLnX0RERPxVZGQk4eHhdvPCwsLIyMg477qHDh0iOzubN998kzfffNPhMuc2WS2oSS3gsFOtypUrc+jQIdv+CtpGlSpVWLhwYYHbPnHiBJdeeinh4eG88MIL1K1bl8jISPbu3Uu/fv1IT08vcN1c//77L4ZhUKlSJYf/r1Wrll2cBR2Pszp27Mhrr71GTk4Of//9NyNGjOD++++nYcOGtvt5rVYrPXr0YP/+/YwYMYLGjRsTFRWF1Wqlbdu2Th8XQOvWrR3+Pyio8Ovvhw4dKrBMcv9fFGc3FQ8NDeXo0aO2iyQNGjQgLi6OpUuXcujQIbv7uXOP5/HHH+fxxx93uG1HTanPjtVRGVeqVIldu3blm1+hQgW76bCwMACnXntHatSoAeBwX4VxJo6BAweyePFiRowYQevWrYmOjsZisdCrVy+H8Z5vmwW9XqVKlcq37r///sucOXMICQlxGH9BZSIiUhIo6RYREYfKlStHcHAwt9xyC/fdd5/DZRITE+2mC7tH9cCBAw7n1alTB8hLAFJSUvItt3//fmJjYwvc9pIlS9i/fz9Lly61JW6AS8NnxcbGYrFY+Omnn2zJx9ly5+XGWdDxFNQZ17liYmJo1aoVAG3atKFNmza24aI2btxIUFAQW7Zs4bfffmPatGkMGjTItu727dtdOi6AL7/80qWa+FwVKlQosEzO3r6rGjVqZEusw8LCqFSpEvXr17f9v1OnTvz444+2xC836c7d37Bhw+jXr5/DbderV6/AYwHydVYHjsvTE1q1akX58uWZPXs2Y8eOLfJ48+dKS0tj7ty5jBw5kqeffto2PzMzk8OHDxdpm2e/XlWrVrXNz87OznexJTY2liZNmvDiiy863FbuRRoRkZJISbeISAlXUM1dZGQkSUlJbNiwgSZNmhAaGnpB+5k+fbpd8/OVK1eyZ88e7rjjDgDatWtHREQEn3zyCddff71tuX379rFkyRJbZ1WO5CYu5ybLjobfOvt4IyIibPOvuuoqxo0bxz///JOvOe3Z2rZtS3h4eIHH42zSfa6LLrqIJ598ktGjRzNjxgxuvPHGIh/X2S6//HJKlSrFjh07Cm3+X5Bu3brx9ddfs3//frvE6X//+x+RkZG0bdvW5W2CWWadO3fm+++/JygoyO5iCZhNzEePHs2hQ4eoUqWKrYO5evXqcdFFF/Hbb7/x0ksvubTPevXqUblyZWbOnMmjjz5qm5+cnMzKlSuLnBgWdE45EhISwlNPPcVTTz3FmDFjeO655/Itc/DgQf7+++8Cm7w7YrFYMAwj37ny3nvvkZOT4/R2ztapUyfA7Pn+7LHOv/zyy3w9kl911VXMmzeP2rVrU65cuSLtT0QkUCnpFhEp4cqUKUPNmjWZPXs23bp1o3z58sTGxpKQkMDrr79Ox44dufTSS7nnnntISEjg+PHjbN++nTlz5rBkyRKn97Nu3TruuOMOrr/+evbu3cszzzxD1apVuffeewFz3OoRI0YwfPhwbr31Vm688UYOHTrE6NGjCQ8PZ+TIkQVuu3379pQrV467776bkSNHEhISwvTp0/ntt9/yLdu4cWMAxo8fzxVXXEFwcDBNmjShQ4cO3HnnnQwZMoR169bRqVMnoqKiSElJYcWKFTRu3Jh77rmHcuXK8fjjj/PCCy/YHc+oUaNcal7uyOOPP86UKVMYPXo0/fv3p379+tSuXZunn34awzAoX748c+bMcdjUPve4Xn/9dQYNGkRISAj16tUjISGB559/nmeeeYadO3fSs2dPypUrx7///suaNWuIiopy2BN9rpEjR9ru133uuecoX74806dP57vvvmPChAnExMQU+XiTkpL48ssvWbBgAZMnT7b7X+fOnTl06BDLly9n4MCBdv975513uOKKK7j88ssZPHgwVatW5fDhw/zxxx+sX7+eL774wuH+goKCGD16NHfddRfXXXcdt912G0ePHmX06NHEx8eft6l9QQo6pwq6UPXEE0/wxx9/MHLkSNasWcPAgQOpXr06aWlpLF++nKlTpzJ69GiXku7o6Gg6derEyy+/bHv/Llu2jPfff5+yZcsW6bgaNmzIjTfeyMSJEwkODqZr1678/vvvTJw4kZiYGLvX6/nnn2fhwoW0b9+eBx98kHr16pGRkcHu3buZN28eU6ZMoVq1akWKQ0TE73m5IzcREXGDgnovj4qKyreso56/Fy1aZDRv3twICwvL1yvxrl27jNtuu82oWrWqERISYsTFxRnt27e39TpuGHm9RH/xxRcFxrZgwQLjlltuMcqWLWvrpfzvv//Ot/x7771nNGnSxAgNDTViYmKMPn365Ov52NExrFy50mjXrp0RGRlpxMXFGXfccYexfv16AzA+/PBD23KZmZnGHXfcYcTFxRkWiyVfj98ffPCB0aZNGyMqKsqIiIgwateubdx6663GunXrbMtYrVZj7NixRvXq1Y3Q0FCjSZMmxpw5c4zOnTs73Xv5lVde6fB///d//2cAxkcffWQYhmFs3brV6N69u1GmTBmjXLlyxvXXX28kJyc77HF+2LBhRpUqVYygoCC7XrsNwzC++eYbIykpyYiOjjbCwsKMmjVrGtddd52xaNGi88a7efNm4+qrrzZiYmKM0NBQo2nTpnavaS5c6L0899gAAzC2bNli9z+r1WqUL1/eAIx3330337q//fab0b9/f6NixYpGSEiIUblyZaNr1652vYKf23t5rqlTpxp16tQxQkNDjbp16xoffPCB0adPH6N58+a2ZXJ7DX/55ZcdHufZr/35zqmCzJ4927jyyiuNuLg4o1SpUka5cuWMpKQkY8qUKUZmZqbLcezbt8+49tprjXLlyhllypQxevbsaWzZssWoWbOm3Xva0edFQa9XRkaG8eijjxoVK1Y0wsPDjbZt2xqrVq0yYmJijEceecRu/f/++8948MEHjcTERCMkJMQoX7680bJlS+OZZ54xTpw4cd7XQ0QkUFkMwzCKMccXEZESZtq0aQwZMoS1a9fa7mEW8SVHjx6lbt26XHPNNUydOtXb4fi8lStX0qFDB6ZPn56vFYKIiOSn5uUiIiJSYhw4cIAXX3yRpKQkKlSowJ49e3jttdc4fvw4Dz30kLfD8zkLFy5k1apVtGzZkoiICH777TfGjRvHRRddVGBHdiIiYk9Jt4iIiJQYYWFh7N69m3vvvZfDhw/bOoObMmUKDRs29HZ4Pic6OpoFCxYwadIkjh8/TmxsLFdccQVjx47NNxyhiIg4publIiIiIiIiIh5StG46RUREREREROS8lHSLiIiIiIiIeIiSbhEREREREREPKVEdqVmtVvbv30+ZMmWwWCzeDkdERERERET8lGEYHD9+nCpVqhAUVHB9dolKuvfv30/16tW9HYaIiIiIiIgEiL1791KtWrUC/1+iku4yZcoA5osSHR1dpG1kZWWxYMECevToQUhIiDvDk2KiMgwMKsfAoHL0fyrDwKByDAwqR/+nMvQvx44do3r16rY8syAlKunObVIeHR19QUl3ZGQk0dHReiP4KZVhYFA5BgaVo/9TGQYGlWNgUDn6P5WhfzrfrcvqSE1ERERERETEQ5R0i4iIiIiIiHiIkm4RERERERERDylR93Q7Kycnh6ysLIf/y8rKolSpUmRkZJCTk1PMkYk7lIQyDA0NLXTYAhERERERKR5Kus9iGAYHDhzg6NGjhS5TuXJl9u7dq7G+/VRJKMOgoCASExMJDQ31digiIiIiIiWaku6z5CbcFStWJDIy0mFCZrVaOXHiBKVLl1ZNop8K9DK0Wq3s37+flJQUatSoEbAXFkRERERE/IGS7jNycnJsCXeFChUKXM5qtXL69GnCw8MDMmErCUpCGcbFxbF//36ys7M13ISIiIiIiBcFZsZRBLn3cEdGRno5EpELl9usPFDvWRcRERER8RdKus+hprgSCHQei4iIiIj4BiXdIiIiIiIiIh6ipFsCyvvvv0+PHj08tv3MzExq1KjBr7/+6rF9iIiIiIhI4FDS7ecGDx6MxWLBYrFQqlQpatSowT333MORI0e8HZpTRo0aRbNmzdyyrczMTJ577jlGjBhRpPUHDx7M008/XegyYWFhPP744zz11FNF2oeIiIiIiJQsSroDQM+ePUlJSWH37t289957zJkzh3vvvderMZ0+fbrY9/nVV19RunRpLr30UpfXtVqtfPfdd/Tp0+e8y95000389NNP/PHHH0UJU0REREREShAl3R6QkpbOyh2ppKSlF8v+wsLCqFy5MtWqVaNHjx4MGDCABQsW2P6fk5PD7bffTmJiIhEREdSrV4/XX3/d9v/NmzcTFBREamoqAEeOHCEoKIjrr7/etszYsWNp165dgTEkJCTwwgsvMHjwYGJiYhg6dCgATz31FHXr1iUyMpJatWoxYsQIW0/x06ZNY/To0fz222+22vpp06YBkJaWxp133knFihWJjo6ma9eu/Pbbb4W+Dp9//jm9e/e2m5ednc2DDz5I2bJlqVChAk899RSDBw/mpptuslvu559/JigoiDZt2nD69Gnuv/9+4uPjCQ8PJyEhgbFjx9qWrVChAu3bt+ezzz4rNB4RERERya+4fyuLeJvG6XazGWuTGTZrM1YDgiwwtl9jBrSuUWz737lzJ/Pnz7cbm9lqtVKtWjVmzpxJbGwsK1eu5M477yQ+Pp7+/fvTqFEjKlSowLJly7j22mtZvnw5FSpUYPny5bZtLF26lM6dOxe675dffpkRI0bw7LPP2uaVKVOGadOmUaVKFTZv3szQoUMpU6YMTz75JAMGDGDLli3Mnz+fRYsWARATE4NhGFx55ZWUL1+eefPmERMTwzvvvEO3bt3466+/KF++vMP9//TTT/mS6fHjxzN9+nQ+/PBDLr74Yl5//XVmz55Nx44d7Zb79ttvufrqqwkKCuKNN97g22+/ZebMmdSoUYO9e/eyd+9eu+UvueQSfvrpp0JfDxERERGx5+3fyiLeoKTbjVLS0m0fIgBWA4bP2kKnunHEx0R4bL9z586ldOnS5OTkkJGRAcCrr75q+39ISAijR4+2TScmJrJy5UpmzpxJ//79sVgsdOrUiaVLl3LttdeydOlSBg0axEcffcTWrVupW7cuK1eu5JFHHik0jq5du/L444/bzTs7AU9ISOCxxx5jxowZPPnkk0RERFC6dGlKlSpF5cqVbcstWbKEzZs3c/DgQcLCwgB45ZVX+Oabb/jyyy+588478+376NGjHD16lCpVqtjNf/PNNxk2bBh9+/YFYPLkycybNy/f+t9++y2vvPIKAMnJyVx00UV07NgRi8VCzZo18y1ftWpVdu/eXejrISIiIiJ5vPVbWcTblHS70a7Uk7YPkVw5hsHu1FMe/SBJSkri7bff5tSpU7z33nv89ddfPPDAA3bLTJkyhffee489e/aQnp7O6dOn7Tow69KlC1OnTgVg2bJljBkzhl27drFs2TLS0tJIT0+nQ4cOhcbRqlWrfPO+/PJLJk2axPbt2zlx4gTZ2dlER0cXup1ff/2VEydOUKFCBbv56enp7Nixw+E66elm86Tw8HDbvLS0NP79918uueQS27zg4GBatGhhd8/5H3/8wb59+7jssssAs0O17t27U69ePXr27MlVV12Vr0f0iIgITp06VehxiIiIiEgeb/1WFvE23dPtRomxUQRZ7OcFWywkxEZ6dL9RUVHUqVOHJk2a8MYbb5CZmWlXsz1z5kweeeQRbrvtNhYsWMDGjRsZMmSIXeLZpUsXfv/9d7Zv386WLVu49NJL6dy5M8uWLWPp0qW0bNmSMmXKnDeOs61evZobbriBK664grlz57JhwwaeeeaZ83ayZrVaiY+PZ+PGjXaPbdu28cQTTzhcp0KFClgsFoe9tlss9oViGPaf9t9++y3du3cnIsL8sG/RogW7du1izJgxpKen079/f6677jq7dQ4fPkxcXFyhxyEiIiIiebz1W1nE25R0u1F8TARj+zUm+EySF2yx8FK/RsV+5W7kyJG88sor7N+/HzDvdW7fvj333nsvzZs3p06dOvlqjHPv637hhRdo2rQp0dHRdkn3+e7nduTnn3+mZs2aPPPMM7Rq1YqLLrqIPXv22C0TGhpKTk6O3bwWLVpw4MABSpUqRZ06dewesbGxDvcVGhpKgwYN2Lp1q21eTEwMlSpVYs2aNbZ5OTk5bNy40W7d2bNn5+uALTo6mgEDBvDuu+8yY8YMvvrqKw4fPmz7/5YtW2jevLlLr4eIiIj4PnXy5Tm+8ltZpLipebmbDWhdg05149ideoqE2EivfIh06dKFhg0b8tJLLzF58mTq1KnD//73P3744QcSExP5+OOPWbt2LYmJibZ1cu/r/uSTT2z3bjdp0oTTp0+zePFiHnroIZfjqFOnDsnJyXz++ee0bt2a7777jq+//tpumYSEBHbt2sXGjRupVq0aZcqU4bLLLqNdu3Zcc801jB8/nnr16rF//37mzZvHNddc47AZO8Dll1/OihUrePjhh23zHnjgAcaOHUudOnWoX78+b775JkeOHLHVfh88eJC1a9fyzTff2NZ57bXXiI+Pp1mzZgQFBfHFF19QuXJlypYta1vmp59+YsyYMS6/JiIiIuK71MmX5/nCb2WR4qaabg+Ij4mgXe0KXv0QefTRR3n33XfZu3cvd999N/369WPAgAG0adOGQ4cOORzHOykpiZycHLp06QKYiXjumNfn9vbtjD59+vDII49w//3306xZM1auXMmIESPslrn22mvp2bMnSUlJxMXF8dlnn2GxWJg3bx6dOnXitttuo27dutxwww3s3r2bSpUqFbi/oUOHMm/ePNLS0mzznnrqKW688UZuvfVW2rVrR+nSpenRo4etg7Y5c+bQpk0bKlasaFundOnSjB8/nlatWtG6dWt2797NvHnzCAoy3y6rVq0iLS0tX5NzERER8V8FdfKlGm/384XfyiLFyWKce4NrADt27BgxMTGkpaXl68wrIyODXbt2kZiYaNcZ17msVivHjh0jOjraloSJ7+jfvz/Nmzdn2LBhDv9vtVq5+OKL6d27N+PHj+eaa66hY8eOPPnkk07v4/rrr6d58+YMHz7cXWG7nbPnsz/Lyspi3rx59OrVy26IPPEvKkf/pzIMDCpHWLkjlYHv/pJv/mdD29KudgUHa/gelaP/Uxn6l8Lyy7Mpa5SA8vLLL1O6dGnb9J49e3j33Xf566+/2Lx5M/fccw+7du2y1VJ37NiRG2+80entZ2Zm0rRp0/MOnyYiIiL+RZ18iYinKOmWgFKzZk274dKCgoKYNm0arVu3pkOHDmzevJkFCxZQr149AJ588kmqV6/u9PbDwsJ49tlnbT2di4iISGBQJ18i4inqSE0CWvXq1fn555/t5uXeIiAiIiJyNnXyJSKeoKRbREREROSM+JgIJdsi4lZqXi4iIiIiIiLiIUq6RURERERERDxESbeIiIiIiIiIhyjpFhEREREREfEQJd0iIiIiIiIiHqKkW3xCQkICkyZN8sq+p02bRtmyZT22/d27d2OxWNi4cSMAS5cuxWKxcPToUY/tU0REREREfIOSbj83ePBgrrnmGpfWsVgsfPPNNx6Jp7h4OlH2pPbt25OSkkJMTIy3QxEREREREQ9T0i1FlpWV5e0Q/FJoaCiVK1fGYrF4OxQREQBS0tJZuSOVlLR0b4ciIiIScJR0B5guXbrw4IMP8uSTT1K+fHkqV67MqFGjbP9PSEgAoG/fvlgsFts0wJw5c2jZsiXh4eHUqlWL0aNHk52dbfu/xWJhypQp9OnTh6ioKF544QVbU+nvvvuOpk2bEh4eTps2bdi8ebNdXF999RUNGzYkLCyMhIQEJk6cWOhxvPrqqzRu3JioqCiqV6/Ovffey4kTJwCzefaQIUNIS0vDYrFgsVhsx3j69GmefPJJqlatSlRUFG3atGHp0qV22542bRqNGjWidOnS9O3bl0OHDhUay+nTp7n//vuJj48nPDychIQExo4da/e6vP3221xxxRVERESQmJjIF198UeD2zm1enltr/8MPP3DxxRdTunRpevbsSUpKit16H374IRdffDHh4eHUr1+ft956q9C4RUScMWNtMh3GLWHgu7/QYdwSZqxN9nZIIiIiAUVJdwD66KOPiIqK4pdffmHChAk8//zzLFy4EIC1a9cCZgKXkpJim/7hhx+4+eabefDBB9m6dSvvvPMO06ZN48UXX7Tb9siRI+nTpw+bN2/mtttus81/4okneOWVV1i7di0VK1akd+/etprwX3/9lf79+3PDDTewefNmRo0axYgRI5g2bVqBxxAUFMQbb7zBli1b+Oijj1iyZAlPPvkkYDbPnjRpEtHR0aSkpJCSksLjjz8OwJAhQ/j555/5/PPP2bRpE9dffz09e/bk77//BuCXX37hjjvu4Pbbb2f9+vUkJSXxwgsvFPp6vvHGG3z77bfMnDmTbdu28cknn9hdrAAYMWIE1157Lb/99hs333wzN954I3/88Ueh2z3bqVOneOWVV/j4449Zvnw5ycnJtmMCePfdd3nmmWd48cUX+eOPP3jppZcYMWIEH330kdP7EBE5V0paOsNmbcZqmNNWA4bP2qIabxERETcq5e0AfNo7neHEQbtZFiDasGKxePB6RemKcNeyIq/epEkTRo4cCcBFF13E5MmTWbx4Md27dycuLg6AsmXLUrlyZds6L774Ik8//TSDBg0CoFatWowZM4Ynn3zSti2AgQMH2iXbu3btAsxkvHv37oCZ9FerVo2vv/6a/v378+qrr9KtWzdGjBgBQN26ddm6dSsvv/wygwcPdngMDz/8sO15YmIiY8aM4Z577uGtt94iNDSUmJgYLBaL3THs2LGDzz77jH379lGlShUAHn/8cebPn8+HH37ISy+9xOuvv06PHj145JFHiI6Opn79+qxcuZL58+cX+HomJydz0UUX0bFjRywWCzVr1sy3zPXXX88dd9wBwJgxY1i4cCFvvvmm07XRWVlZTJkyhdq1awNw//338/zzz9v+P2bMGCZOnEi/fv1sr0nuxZHcMhMRcdWu1JO2hDtXjmGwO/UU8TER3glKREQkwCjpLsyJg3B8v90sy5mHL2vSpInddHx8PAcPHixgadOvv/7K2rVr7Wq2c3JyyMjI4NSpU0RGRgLQqlUrh+u3a9fO9rx8+fLUq1fPVtP7xx9/0KdPH7vlO3TowKRJk8jJySE4ODjf9n788Udeeukltm7dyrFjx8jOziYjI4OTJ08SFRXlMIb169djGAZ169a1m5+ZmUmFChVssZzb8Vy7du0KTboHDx5M9+7dqVevHj179uSqq66iR48eBR5/7nRub+XOiIyMtCXcYF9m//33H3v37uX2229n6NChtmWys7PVGZuIXJDE2CiCLNgl3sEWCwmxkd4LSkREJMAo6S5M6Yr5ZhmAcaam22PJt4P9uiIkJMRu2mKxYLVaC13HarUyevRoW03q2cLDw23PC0p4HcntKMwwjHydhhmG4WgVAPbs2UOvXr24++67GTNmDOXLl2fFihXcfvvthXbeZrVaCQ4O5tdff82XyJcuXfq8+y1IixYt2LVrF99//z2LFi2if//+XHbZZXz55ZeFrudKR2mOyiw31tyye/fdd2nTpo3dco4uWIiIOCs+JoKx/RozfNYWcgyDYIuFl/o1Ui23iIiIGynpLoyDJt6G1cqxY8eIjo7GEuSft8SHhISQk5NjN69FixZs27aNOnXqFGmbq1evpkaNGgAcOXKEv/76i/r16wPQoEEDVqxYYbf8ypUrqVu3rsOkcd26dWRnZzNx4kSCzrzGM2fOtFsmNDQ03zE0b96cnJwcDh48yKWXXuowzgYNGrB69Woeeughu9jPJzo6mgEDBjBgwACuu+46evbsyeHDhylfvrxtG7feeqvdNps3b37e7TqjUqVKVK1alZ07d3LTTTe5ZZsiIrkGtK5Bp7px7E49RUJspBJuERERN1PSXQIlJCSwePFiOnToQFhYGOXKleO5557jqquuonr16lx//fUEBQWxadMmNm/efN6OxgCef/55KlSoQKVKlXjmmWeIjY21NeN+7LHHaN26NWPGjGHAgAGsWrWKyZMnF3i/c+3atcnOzubNN9/k6quv5ueff2bKlCn5juHEiRMsXryYpk2bEhkZSd26dbnpppu49dZbmThxIs2bNyc1NZUlS5bQuHFjevXqxYMPPkj79u15/fXXGTBgAIsWLSq0aTnAa6+9Rnx8PM2aNSMoKIgvvviCypUr240T/sUXX9CqVSs6duzI9OnTWbNmDe+///55XzdnjRo1igcffJDo6GiuuOIKMjMzWbduHUeOHOHRRx91235EpGSKj4lQsi0iIuIh/llVKxdk4sSJLFy4kOrVq9tqYy+//HLmzp3LwoULad26NW3btuXVV1912GmYI+PGjeOhhx6iZcuWpKSk8O233xIaGgqYtegzZ87k888/p1GjRjz33HM8//zzBXai1qxZM1599VXGjx9Po0aNmD59ut0QXWD2YH733XczYMAA4uLimDBhAmD2yn7rrbfy2GOPUa9ePXr37s0vv/xC9erVAWjbti1Tp05l6tSptGjRggULFvDss88WemylS5dm/PjxtGrVitatW7N7927mzZtnq4UHGD16NJ9//jlNmjTho48+Yvr06TRo0MCp184Zd9xxB++99x7Tpk2jcePGdO7cmWnTppGYmOi2fYiIiIiIiPtZjKLc5Oqnjh07RkxMDGlpaURHR9v9LyMjg127dpGYmGh3D/O5rGc1Lw/y0+bl7rR06VKSkpI4cuSIXc2vL3N3GVosFr7++ut8HbR5k7Pnsz/Lyspi3rx59OrVK9898eI/VI7+T2UYGFSOgUHl6P9Uhv6lsPzybMoaRURERERERDxESbeIiIiIiIiIh6gjNbkgXbp0KdIwXIGkpB+/iIiIiIgUTDXdIiIiIiIiIh6ipFtERERERETEQ5R0n8NqtXo7BJELpibvIiIiIiK+Qfd0nxEaGkpQUBD79+8nLi6O0NBQLBZLvuWsViunT58mIyNDQ4b5qUAvQ8Mw+O+//7BYLBpqQkREREqslLR0dqWeJDE2iviYCG+HIyWYku4zgoKCSExMJCUlhf379xe4nGEYpKenExER4TApF99XEsrQYrFQrVo1goODvR2KiIiISLGbsTaZYbM2YzUgyAJj+zVmQOsa3g5LSigl3WcJDQ2lRo0aZGdnk5OT43CZrKwsli9fTqdOnVSL6KdKQhmGhIQo4RYREb+g2khxt5S0dFvCDWA1YPisLXSqG6dzTLxCSfc5cpvkFpSMBQcHk52dTXh4eMAmbIFOZSgiIuIbVBspnrAr9aQt4c6VYxjsTj2lpFu8IvBuaBURERERn1dQbWRKWrp3AxO/lxgbRdA5dxAGWywkxEZ6JyAp8ZR0i4iIiEixK6w2UuRCxMdEMLZfY4LP9N0TbLHwUr9GquUWr/Gr5uX//PMPTz31FN9//z3p6enUrVuX999/n5YtW3o7NBERERFxQW5t5NmJt2ojxV0GtK5Bp7px7E49RUJspBJu8Sq/qek+cuQIHTp0ICQkhO+//56tW7cyceJEypYt6+3QRERERMRFqo0UT4uPiaBd7Qo6p8Tr/Kame/z48VSvXp0PP/zQNi8hIcF7AYmIiIjIBVFtpIiUBH5T0/3tt9/SqlUrrr/+eipWrEjz5s159913vR2WiIiIiFwA1UaKSKDzm5runTt38vbbb/Poo48yfPhw1qxZw4MPPkhYWBi33nqrw3UyMzPJzMy0TR87dgwwx2nOysoqUhy56xV1ffE+lWFgUDkGBpWj/1MZBgaVY2BQOfo/laF/cbacLIZhGOdfzPtCQ0Np1aoVK1eutM178MEHWbt2LatWrXK4zqhRoxg9enS++Z9++imRkeqkQ0RERERERIrm1KlTDBw4kLS0NKKjowtczm9quuPj42nQoIHdvIsvvpivvvqqwHWGDRvGo48+aps+duwY1atXp0ePHoW+KIXJyspi4cKFdO/enZCQkCJtQ7xLZRgYVI6BQeXo/1SGgUHlGBhUjv5PZehfcltSn4/fJN0dOnRg27ZtdvP++usvatasWeA6YWFhhIWF5ZsfEhJywSexO7Yh3qUyDAwqx8CgcvR/KsPAoHIMDCpH/6cy9A/OlpHfdKT2yCOPsHr1al566SW2b9/Op59+ytSpU7nvvvu8HZqIiIiIiIiIQ36TdLdu3Zqvv/6azz77jEaNGjFmzBgmTZrETTfd5O3QRERERERERBzym+blAFdddRVXXXWVt8MQERERERERcYrf1HSLiIiIiIiI+Bsl3SIiIiIiIiIeoqRbRERERERExEOUdIuIiIiIiIh4iJJuEREREREREQ9R0i0iIiIiIiLiIUq6RURERERERDxESbeIiIiIiIiIhyjpFhEREREREfEQJd0iIiIiIiIiHqKkW0RERERERMRDlHSLiIiIiIiIeIiSbhEREREREREPUdItIiIiIiIi4iFKukVEREREREQ8REm3iIiIiIi4JCUtnZU7UklJS/d2KCI+r5S3AxAREREREf8xY20yw2ZtxmpAkAXG9mvMgNY1vB2WiM9STbeIiIiIiDglJS3dlnADWA0YPmuLarxFCqGkW0REREREnLIr9aQt4c6VYxjsTj3lnYBE/ICSbhERERERcUpibBRBFvt5wRYLCbGR3glIxA8o6RYREREREafEx0Qwtl9jgi1m5h1ssfBSv0bEx0R4OTIR36WO1ERERERExGkDWtegU904dqeeIiE2Ugm3yHko6RYREREREZfEx0Qo2RZxkpqXi4iIiIiIiHiIkm4RERERERERD1HSLSIiIiIiIuIhSrpFREREREREPERJt4iIiIifSklLZ+WOVFLS0r0dioiIFEC9l4uIiIj4oRlrkxk2azNWA4IsMLZfYwa0ruHtsERE5Byq6RYRERHxMylp6baEG8BqwPBZW1TjLSLig5R0i4iIiPiZXaknbQl3rhzDYHfqKe8EJCIiBVLSLSIiIuJnEmOjCLLYzwu2WEiIjfROQCIiUiAl3SIiIiJ+Jj4mgrH9GhNsMTPvYIuFl/o1Ij4mwsuRiYjIudSRmoiIiIgfGtC6Bp3qxrE79RQJsZFKuEVEfJSSbhERERE/FR8ToWRbRMTHqXm5iIiIiIiIiIco6RYRERERERHxECXdIiIiIiIiIh6ipFtERERERETEQ5R0i4iIiIiIiHiIy72XZ2ZmsmbNGnbv3s2pU6eIi4ujefPmJCYmeiI+EREREREREb/ldNK9cuVK3nzzTb755htOnz5N2bJliYiI4PDhw2RmZlKrVi3uvPNO7r77bsqUKePJmEVERERERET8glPNy/v06cN1111H1apV+eGHHzh+/DiHDh1i3759nDp1ir///ptnn32WxYsXU7duXRYuXOjpuEVERERERER8nlM13T169OCLL74gNDTU4f9r1apFrVq1GDRoEL///jv79+93a5AiIiKFSUlLZ1fqSRJjo4iPifB2OCIiIiI2TiXd9913n9MbbNiwIQ0bNixyQCIiIq6YsTaZYbM2YzUgyAJj+zVmQOsa3g5LREREBChC7+V79+5l3759tuk1a9bw8MMPM3XqVLcGJiIicj4paem2hBvAasDwWVtISUv3bmAiIiIiZ7icdA8cOJAff/wRgAMHDtC9e3fWrFnD8OHDef75590eoIiISEF2pZ60Jdy5cgyD3amnvBOQ+K2UtHRW7kjVBRvxS3nnb4a3QxERB1weMmzLli1ccsklAMycOZNGjRrx888/s2DBAu6++26ee+45twcpIiLiSGJsFEEW7BLvYIuFhNhI7wUlfke3KIg/O/f87Z9ooRfq60LEl7hc052VlUVYWBgAixYtonfv3gDUr1+flJQU90YnIiJSiPiYCMb2a0ywxQKYCfdL/RrpB6Y4TbcoiD9zdP7O2BnEuyt20WHcEga++wsdxi1hxtpkl7apVh8i7uVyTXfDhg2ZMmUKV155JQsXLmTMmDEA7N+/nwoVKrg9QBERkcIMaF2DTnXj2J16ioTYSCXc4pLCblHQuSS+ztH5a2Dh5QV/Y5xzIalT3bjzntNq9SHiGS7XdI8fP5533nmHLl26cOONN9K0aVMAvv32W1uzcxERkeIUHxNBu9oVlCSJy3JvUTibblEQf+Ho/AXDlnDncqavC7X6EPEcp5PuEydOANClSxdSU1NJTU3lgw8+sP3/zjvvZMqUKe6PUERERMRDdIuC+LNzz98gC/SuYS3ShSR1TCniOU43L4+NjSUpKYnevXvTp08fqlSpYvf/hIQEd8cmIiIi4nG6RUH82dnnb9WYUDb8vIS2LRowYvYf5BiG0xeS1DGliOc4nXRv27aNb7/9lq+++oqHH36YJk2a2BLwJk2aeDJGEREREY+Kj4lQsi1+K/f8zcrKYgNwfctqJF1c2aULSbm15sNnbXEpWReR83M66a5ZsyYPPPAADzzwAGlpacybN4/Zs2czceJEypUrZ0vAO3fuTHBwsCdjFhERERGRQhTlQpJafYh4hssdqQHExMRw44038vnnn5Oamso777yD1WplyJAhxMXFMX36dHfHKSIiIiWAhisS8S51TCnifi4PGZZvA6VK0aNHD3r06MGbb77J+vXrycnJcUdsIiIiUoJouCIREQlELifdmzZtcjjfYrEQHh5Ow4YNCQsLu+DAREREpOQoaLgiZ8YWFhER8WUuJ93NmjXDYsk3IKBNSEgIAwYM4J133iE8PPyCghMREZGSobDhipR0i4iIP3P5nu6vv/6aiy66iKlTp7Jx40Y2bNjA1KlTqVevHp9++invv/8+S5Ys4dlnn/VEvCIiIhKAcocrOpuGKxIRkUDgck33iy++yOuvv87ll19um9ekSROqVavGiBEjWLNmDVFRUTz22GO88sorbg1WREREApOGKxIRkUDlck335s2bqVmzZr75NWvWZPPmzYDZBD0lJeXCoyvE2LFjsVgsPPzwwx7dj4iIiBSPAa1rsOLpJD4b2pYVTyepEzUREQkILifd9evXZ9y4cZw+fdo2Lysri3HjxlG/fn0A/vnnHypVquS+KM+xdu1apk6dSpMmTTy2DxERESl+Gq7IP2moNxGRgrncvPz//u//6N27N9WqVaNJkyZYLBY2bdpETk4Oc+fOBWDnzp3ce++9bg8W4MSJE9x00028++67vPDCCx7Zh4iIiASulLR0dqWeJDE2Ssm9G2ioNxHn6LOn5HI56W7fvj27d+/mk08+4a+//sIwDK677joGDhxImTJlALjlllvcHmiu++67jyuvvJLLLrvsvEl3ZmYmmZmZtuljx44BZs18VlZWkfafu15R1xfvUxkGBpVjYFA5+j9/K8Mvft3Hs7O32hLEF/o04PqW1bwdltcVtRxT0jLyDfU2bNZm2iWWIz5Go9gUN397P5Ykzn72qAz9i7PlZDEMwzj/Yr7h888/58UXX2Tt2rWEh4fTpUsXmjVrxqRJkxwuP2rUKEaPHp1v/qeffkpkpHpDFRERKUmOZsKo9cEY5HWTbsFgVIscyoZ5MTA/9neahclbg/PNv79BDhfF+M1PTBGP0mdP4Dp16hQDBw4kLS2N6OjoApdzqqZ71apVtGvXzqkdnzx5kt27d9OwYUPnInXS3r17eeihh1iwYIHT438PGzaMRx991DZ97NgxqlevTo8ePQp9UQqTlZXFwoUL6d69OyEhIUXahniXyjAwqBwDg8rR//lTGa7eeRhj/Tq7eQYWajdrS5vE8l6KyjcUtRxT0jJ464/ldmOsB1mgf68k1XR7gT+9H0sSVz57VIb+Jbcl9fk4lXTfeuutJCQkMHToUHr16kXp0qXzLbN161Y++eQTPvzwQyZMmOD2pPvXX3/l4MGDtGzZ0jYvJyeH5cuXM3nyZDIzMwkOtr/SGhYWRlhY/stHISEhF3wSu2Mb4l0qw8CgcgwMKkf/5w9lWKdyNEEW7BLEYIuF2pWifT724uJqOdaIDXE41FuN2DIejFLOxx/ejyVJUT57VIb+wdkycirp3rp1K++88w7PPfccN910E3Xr1qVKlSqEh4dz5MgR/vzzT06ePEm/fv1YuHAhjRo1uqDgHenWrZttSLJcQ4YMoX79+jz11FP5Em4RERGRs2kscM8Y0LoGnerGsTv1FAmxkXo9Rc6hzx5xKukOCQnh/vvv5/7772f9+vX89NNP7N69m/T0dJo2bcojjzxCUlIS5ct7rmlWmTJl8iXzUVFRVKhQwSNJvoiIiAQeJYieER8ToddSpBD67CnZXO69vEWLFrRo0cITsYiIiIh4nBJEEfEGffaUXC4n3b5k6dKl3g5BREREREREpEBB3g5AREREREREJFAp6RYRERERERHxECXdIiIiIiIiIh5yQUl3RkaGu+IQERERERERCTguJ91Wq5UxY8ZQtWpVSpcuzc6dOwEYMWIE77//vtsDFBEREZGSJyUtnZU7UklJS/d2KCIiF8TlpPuFF15g2rRpTJgwgdDQUNv8xo0b895777k1OBEREREpeWasTabDuCUMfPcXOoxbwoy1yd4OyW10MUGk5HE56f7f//7H1KlTuemmmwgODrbNb9KkCX/++adbgxMRERGRkiUlLZ1hszZjNcxpqwHDZ20JiCQ1kC8miEjBXE66//nnH+rUqZNvvtVqJSsryy1BiYiIiIjn+HJt667Uk7aEO1eOYbA79ZR3AnKTQL6YICKFK+XqCg0bNuSnn36iZs2advO/+OILmjdv7rbARERERMT9ZqxNtiV/QRYY268x/ZrFezssm8TYKIIs2CXewRYLCbGR3gvKDQq7mBAfE+GdoESkWLicdI8cOZJbbrmFf/75B6vVyqxZs9i2bRv/+9//mDt3ridiFBERERE3KKi2tV1iOe8Gdpb4mAjG9mvM8FlbyDEMgi0WXurXyO8T00C9mCAi5+dy0n311VczY8YMXnrpJSwWC8899xwtWrRgzpw5dO/e3RMxioiIiIgbFFTbmnzYt5puD2hdg05149ideoqE2Ei/T7ghcC8miMj5uZx0A1x++eVcfvnl7o5FRERERDyooNrWGuUjOeS9sByKj4nweEKakpbOrtSTJMZGFUvyG4gXE0Tk/IqUdIuIiIiI/ym4tjWcDd4Orpg5urd9QOsaHt9vcVxMEBHf4nLSXa5cOSwWS775FouF8PBw6tSpw+DBgxkyZIhbAhQRERER93FU21rSRqAp6N72TnXjlBCLiNu5nHQ/99xzvPjii1xxxRVccsklGIbB2rVrmT9/Pvfddx+7du3innvuITs7m6FDh3oiZhERERG5ACW9tlU9iYtIcXI56V6xYgUvvPACd999t938d955hwULFvDVV1/RpEkT3njjDSXdIiIiIuJz1JO4iBQnl5PuH374gfHjx+eb361bNx577DEAevXqxdNPP33h0YmIiIiIuJlLPYlbcyArHbIzzEdWRt5zh9NOLlvY/wwrVG0JdbpBncugUiNwcHuniPgHl5Pu8uXLM2fOHB555BG7+XPmzKF8+fIAnDx5kjJlyrgnQhERERERW/KbCdln/tpNn5vUFrZsJgOy0ulT7ySZ6elEWE4Tui4LVjtY1prtnePd/ZP5WDQKSlc+k4B3g1pJEFneOzGJSJG4nHSPGDGCe+65hx9//JFLLrkEi8XCmjVrmDdvHlOmTAFg4cKFdO7c2e3BioiIiIiXWXPOU1Ob7kRtsKPk+DzLWt3f2Vv4mYdXlQrPe4Sc+Xv6FKQl5y1z4gBsnG4+LEFnasEvMx9VmkNQsPfiL2GKe5g5CQwuJ91Dhw6lQYMGTJ48mVmzZmEYBvXr12fZsmW0b98ewNbMXERERCTQBeSP8OMHYP8G8/HPevj3d8g87rHk1ycEh+UlvaXCISQCSoVBqTN/z54+e7mzk2W7aQfrOZouqNn4oR2wfTFsX2TWeGedMucbVti31nwsHQsR5aB2VzMBr9mp+F6vEshbw8yJ/yvSON0dOnSgQ4cO7o5FREQkIAVkUiZAgPwIP3kIy9611D0wm+AvPoOUjXA8xXvxBIedlbiem6wWluQ6kSwXtJ3gMAgK8t4xO1Khtvloc6dZ25+8ykzAty+G//7IWy79CGz5CrZ8RQjQOaIGQRG/Qt0eUL0NBId47RACiYaZkwtRpKTbarWyfft2Dh48iNVqtftfp066wiYiIpIrIJIyccgvf4SnHzWT6txa7P0b4GgypYCLAQrKtcPLQumK7ktyC9qOLya/viAkHGonmY/LX4S0fXm14DuXQWaabdGy6cmw8nXzEVoGanU27wWv3Q3K1fTiQfg3DTMnF8LlpHv16tUMHDiQPXv2YBj2Z57FYiEnJ8dtwYmIiO9S7e35+WVSJk4r7Ed47v+9+v7IPAEHNpnNw3MT7MM7zr9eWDRUaWbeK5z7KFtTvWf7kphq0HKQ+cjJhn/WwfZFWP9eSFDKxrzlTh+HP+eaD4DYumfuBe8GNTuYFzrEKRpmTi6Ey0n33XffTatWrfjuu++Ij4/Hog9gEZESR7W3zlHNSGAr6Ef4pn+OctN7q4v3/ZGVDge2nEmuzyTZ/20DjMLXC4mE+KbkVG7Chn8tNL1iMCEV66m22Z8El4IabaFGW3IufYofZn9O91pBlNq1FHYshlOH8pZN/ct8rH7LbF2Q0DGvQ7YKdXRhpRAuDTMncg6Xk+6///6bL7/8kjp16ngiHhER8XGqvXWeakYCm6Mf4U/2rMf47//07Psj+zQc/P2sjs42wMGtYJyntWFwGFRubNZcV21h/o2tC0HBWLOy+GfePJpWqKOE282Ku1XQ6ZBojMa9oMVNYLWatxPkNkXft8bsiA3MTvG2LzIfAGVr5CXgiZ0gTMP/nmtA6xp0qhvH7tRTJMRG6jtPnOZy0t2mTRu2b9+upFtEpIRS7a3zVDMS+M79Ee7290dONvz3p/092P9ugZzTha8XVAoqNTyriXgLqHixOtUqZl5vFRQUZF5gqdoCOj9hdrq2c1leh2zH9+ctezQZ1n1gPoJKQY12Z8YGvwwqNVIt+BnxMRH6DBeXuZx0P/DAAzz22GMcOHCAxo0bExJi/+HdpEkTtwUnIiK+R7W3rlHNSOA790d4kd8fVisc2p7XPHz/BkjZZI57XRhLEMTVNxPrKs3Mv5Uamp1vidf4ZKugiHLQ8BrzYRhw8I+82u7kVXkXc6zZ5jBlu3+CRaOgdKW8e8FrJUFkee/EL+KnXE66r732WgBuu+022zyLxYJhGOpITUSkBFDtretUM1JyOP3+MAw4sitvHOz9GyHlN7Pjq/OpcJF9E/HKjSE0yiPHI0Xn862CLBao1MB8dHjQ7Hhv94ozSfhCOLI7b9kT/8LG6ebDEgRVW+Y1Ra/SHIKCvXYYIv7A5aR7165dnohDRET8iGpvRQqW7/0RHQ5H957VRPxMkp1x9PwbK5dg30Q8vgmEx3j4CMQd/K5VUFhpqNfTfAAc2pF3L/junyDL7JUfwwr71pqPpWPN2vPaXc0EvHZXKFPZe8cg4qNcTrpr1tT4fiIiotpbkQIdP0D8gQ3E798Aq88k2if/O/960dXyD9WlZrx+y+9bBVWobT7a3AlZGWbz89x7wf/7I2+59COw5SvzAWbLi9xa8GqXQKlQ78Qv4kNcTrpzbd26leTkZE6ftu/Io3fv3hcclIiIiIhfOHnIvpOz/RvsO6cqSFTFvObhuY/SFT0frxSrgGkVFBIOtZPMx+UvQtq+vFrwnUsh81jesgc2m48Vr0FoGajV2bwXvHY3KKfKOymZXE66d+7cSd++fdm8ebPtXm7ANl637ukWERGRgJR+1Lzv+uyOzo4mn3+9iHJ5zcNzE+zoKuoNuoQIyFZBMdWg5SDzkZMF+9bldciWsjFvudPH4c+55gPMIepyO2Sr2QFCAux1ESmAy0n3Qw89RGJiIosWLaJWrVqsWbOGQ4cO8dhjj/HKK694IkYRERGR4pV5Ag5sOqujsw1weMf51wuLhvim9h2dla2pBFsCV3AI1GxnPrqNgBP/wY4lZgK+YzGcOpS3bOpf5mP1W1AqHBI65jVFr1BH7xMJWC4n3atWrWLJkiXExcURFBREUFAQHTt2ZOzYsTz44INs2LDBE3GKiIiIeEZWOhzYYt/R2X/bAKPw9UIi8xLs3Jrs8rXMsZFFSqrScdB0gPmwWs2a79ym6PvWmB2xAWRn5NWOA5StkZeAJ3aCsDJeOwQRd3M56c7JyaF06dIAxMbGsn//furVq0fNmjXZtm2b2wMUERERcZvs03Dwd/t7sA/+YY5LXJjgMLODqLPvwY6rp6GSRAoTFGS2+KjaAjo/YXa6tnNZXodsZ/d/cDQZ1n1gPoJKQY12ZjP0OpdBpUaqBRe/5nLS3ahRIzZt2kStWrVo06YNEyZMIDQ0lKlTp1KrVi1PxCgiIiLiupxsSN2W1zx8/wb4dwvknC58vaBSULGBfUdncRerF2aRCxVRDhpeYz4Mw7zglVvbnbwq771pzTaHKdv9EywaBaUr5d0LXitJvfqL33E56X722Wc5efIkAC+88AJXXXUVl156KRUqVGDGjBluD1BERETkvKxWOLT9rHGwN0DKJshOL3w9SxDE1bdvIl6podlbs4h4jsUClRqYjw4Pmv0o7F5xJglfCEd25y174l/YON18WIKgasu8puhVmqvFifg8l5Puyy+/3Pa8Vq1abN26lcOHD1OuXDlbD+YiIiIiLrNazRouazYYOWeeO5qXY9aI/ffnmSR7o/k4ffz8+6hwUV6CXbWF2WQ8NMrTRyYi5xNWGur1NB8Ah3bk3Qu++yfIOmXON6ywb635WDrWrD2v3dVMwGt3hTKVvXcMIgUo8jjdZytfXk08REQkcKSkpbMr9SSJsVHuH+rHarVPHq3Z5o/I3MSysHnWnPzr2s3LnV/AvPMmswUkuHbr2u8/OCeLDof+I/ijyWfF7ExMVqzWLAxrDkFGDpbz3VNdFGVrkl6xKf+WvpgytVpToc4lEB7j/v2IiPtVqG0+2twJWRlm8/Pce8H/+yNvufQjsOUr8wHmhbTcWvBql+i2EPEJLifdJ0+eZNy4cSxevJiDBw9itVrt/r9z5063BSciIlLcZqxNZtiszVgNCLLA2H6NGdC6RuErnTgIX98FqX+flZSenaSeNe98PWL7mSAgFuBE0dZ1m+iq9p2cVWnOjN9P5JXlqtOM7ZfGgNbOJ90evfgiIs4LCYfaSebj8hchbZ+tFvz030sIzT7rA+jAZvOx4jUILQO1Opv3gtfuBuVqeu8YpERzOem+4447WLZsGbfccgvx8fFqUi4iIgEjJS3dlqQBWA0YPmsLnerGFZx0GQbMedgcl1bAEmx2RBZ01l/bvFIQFEQ2wew5nEE2weQQTA4WcgimQdVyhIaEnln3zDp22ztrO5ZgiKlmNhGPbwZlKtmFYZblKtfK8ixFuvgiIsUjphq0HERKnf503riAJmync/AmOgf9RpOgXXnLnT4Of841HwCxdfM6ZKvZAUJ0MU2Kh8tJ9/fff893331Hhw4dPBGPiIiI1+xKPWlL0nLlGAa7U08VnKj9Pgu2fWc+Dw6DqFgHiWawk/OC7JNKh/McJJ8Fzjsr+b2QeXZx2s/LshrM/2EhPXtdRUhomFPD+qzZkcrAd3/JN/+zy9rSrnYFV4vNoSKV5RlFuvjiJNWei7jPrtSTnDZKsY76rMuuz0T6U4E0Pu+WzkXHVsOOxXDqUN4KqX+Zj9VvQalwSOho1oBf1ANi63jvQCTguZx0lytXTvdwi4hIQEqMjSLIgl2yFmyxkBAb6XiFk4dg3pN50/3egYZ9PRukr8nKwhoUciYJd671m8uvcxFcyD4uJGEvjGrPRdzL0fv8qKUspS/pCzG3m/1UpGzM65Bt3xqz7wmA7Iy84cp+GAaXPgZdR2g8cPEIl2+nGjNmDM899xynTp3yRDwiIiJeEx8Twdh+jQk+86Mr2GLhpX6NCk605j8Fp1LN5/WvggbXFE+gfs7l1xmzhnjljlRS0s4zBNgF7CNX7g/5s13oRYGCas+dPR4Rye+87/OgIPMWlM5PwO0/wJM74fqPoPktUKaK/cZ+mgiLR5u3DIm4mVM13c2bN7e7d3v79u1UqlSJhIQEQkJC7JZdv369eyMUEREpRgNa16BT3Th2p54iITay4CRt23zY/IX5PDwGrpyoGhIXOP06U/QaYlf2cbbcH/LDZ20hxzBcStgL4qnac5GSzqX3eUQ5aHiN+TAMOPgHbPnSTLjB7HzNEgxdn9XnubiVU0n3Nddc4+EwREREfEd8TEThP9wy0mDuI3nTl48tEWPDuvt+5PO+zlz4/dXO7MORoibsBSmOJvUiJVWR3ucWC1RqAJWeg+gq8N1j5vyfXjFvl0ka7v5ApcRyKukeOXKkp+MQERHxHwtGwPH95vPa3aDZQO/GUwy8dT+yN2uIi5qwF7Qtd9eei4ibtL7DrPme97g5vWy8WePd5SnvxiUBw+WO1NauXYvVaqVNmzZ283/55ReCg4Np1aqV24ITERHxOTuXwfqPzOehpeHqSQHfDLGw2ubYSJd/SrgkkGqI3V17LiJudMlQsOaYfXUALH3JvCe80xPejUsCgssdqd13333s3bs33/x//vmH++67zy1BiYiI+KTTJ+HbB/KmLxsFZQO39+nczst+3XOkwNpmT7uQDtF8UXxMBO1qV/Db+EUCWtu7zduFci15Ie9+b5EL4PLl6a1bt9KiRYt885s3b87WrVvdEpSIiIhPWvICHN1jPq/RHlrd7t14POjs5uQWzMfZeXdx1jarhlhEik27e8HIgQXPmtOLnzebmnd82KthiX9zuaY7LCyMf//9N9/8lJQUSpXybBMzEREpWVwdJsqj9q6F1W+bz0uFQ+83zaaHAejc5uS5yXbuMFreqG1WDbGIeJLd9037B+Cy0Xn/XDQSVr7pveDE77mcJXfv3p1hw4Yxe/ZsYmJiADh69CjDhw+ne/fubg9QRERKJm913OVQdibMvg9b+tllGMTW8U4sxcBR52UG8OYNzalQOky1zSISUBx+33R82KzxXvy8udCCZ8ESBO10O624zuWke+LEiXTq1ImaNWvSvHlzADZu3EilSpX4+OOP3R6giIiUPBc6TJTbLX8ZUreZz+ObQbv7iz+GYlRQ52UtE8op2RaRgFLo982lj4FhNW8tAvhhuNnUvO3d3gtY/JLL7eKqVq3Kpk2bmDBhAg0aNKBly5a8/vrrbN68merVq3siRhG/4FPNYEX8XGHDRLnbed+7BzbDitfM50GloM//QXBg304VaJ2XiYgU5LzfN52egC5njdk9/ylY827xBSgBoUi/GqKiorjzzjvdHYuI3yqoGWxKWjq7Uk+SGBulH6siLiiuYaLO24Q9J9tsVm7NNqcvfQwqN3LLvn3980Gdl4lISeDU902Xp8ym5svGm9PzHjeHimx9R/EGK34rsC/VixSDgpolHU3PYvz3f/rG/agifia3pnX4rC3kGIZHalqdasK+8g1I+c18HncxXPq4U9s9XzLtU/erFyI+JkLJtogENKe/b7oMM8fx/ukVc/q7x8ym5q2GFH/Q4neUdItcoIKaJY37/k8MX7kfVcQPebqmtbAmhfExEZD6NywdZ/7DEmQ2Ky8VWug2nUmmfe5+dRGREs6p7xuLBbo+a9Z4595yNPdh8/uh5aBijVf8T2COdSIlljfuq85tlnS2ILAl3Lk8dT+qO+m+9MARKGXpyWGiHL13bU0KrVaYfT/kZJr/aHsvVGtZ6PYKSqbPLYPivF9dRESc49T3jcUC3UZCh4fy5s15CDZ84vkAxa8p6ZaAMWNtMh3GLWHgu7/QYdwSZqxNLpb9Oupw6Kkr6hf8Y95Heev1E/dTWTqn0M7C1r4Le1ebC5ZLhKRnzrs9Z5PpQpN9ERHxbRaLOYa3bRQLw7xIu/FTr4Ylvs3l5uW1atVi7dq1VKhQwW7+0aNHadGiBTt37nRbcCLO8nZzTUfNkspGhnj0flR38vbrJ+6jsnSNwyaFR/bAotF5C/V+E0LPnxA72/lbcdyvLiIiHmSxQI8XzOHEVr8FGPDNvWZT86Y3eDs68UEuJ927d+8mJycn3/zMzEz++ecftwTlyNixY5k1axZ//vknERERtG/fnvHjx1OvXj2P7VP8x3nvzSwG53Y45E89//rC6yfuobJ0nd171zBgzoOQddKcbjkEEi91ejvOJtP+9PkgIiIOWCxw+Utm4v3LFMzE+x6zc7Um13s7OvExTifd3377re35Dz/8QExMjG06JyeHxYsXk5CQ4NbgzrZs2TLuu+8+WrduTXZ2Ns888ww9evRg69atREVFeWy/4h+Ka3ghV/lLz7+++vqJ61SWF2jDJ7Bzqfk8uip0f96l1V1Jpv3l80FERApgsUDPcWav5mvfNRPwr+805ze+ztvRiQ9xOum+5pprALBYLAwaZN9DX0hICAkJCUycONGtwZ1t/vz5dtMffvghFStW5Ndff6VTp04e26/4BzXXvDB6/QKHyvICHEuBH866d/uqSRAe7fJmlEyLBBZnhgGUEsxigV4vm72ar/vATLxn3QlBwdCwr7ejEx/hdNJttVoBSExMZO3atcTGxnosKGekpaUBUL58ea/GIb5DzTUvjF6/wKGyLALDMMdczTS/W2gyAOr28G5MIuJ1zgwDKP7H7RdSLBboNdFMuH+dZibgX95u3uPdoM+Fb1/8nsv3dO/atcsTcbjEMAweffRROnbsSKNGjQpcLjMzk8zMTNv0sWPHAMjKyiIrK6tI+85dr6jri2fFRpYitoZZM1VQGakMC+bM6+crVI6F85ey9JVytGz9mlLbvgPAiIoju9sY8OHXzZf4ShnKhVE55peSlpGvY8phszbTLrEc8THh3g2uACrH8/vi1308O3ur7ULKC30acH3Lau7Z+OUTCM7JJmjjJ2DkYHx5Gzl938eof6XTm1AZ+hdny8liGOeOJnx+ixcvZvHixRw8eNBWA57rgw8+cHVzLrvvvvv47rvvWLFiBdWqFfwmGTVqFKNHj843/9NPPyUyUvc3iogIhGYfp+sfTxOWfRyAtQn3s7/cJV6OSkS87e80C5O3Buebf3+DHC6Kcfnns187mgn/ZViICzcoG+btaIruaCaMWh+MQd64jRYMRrXIcd9xGVaaJb9PzcM/AWAlmLWJ93OgbEs37UB8yalTpxg4cCBpaWlERxd8S5rLSffo0aN5/vnnadWqFfHx8Vgs9oONfv3110WL2EkPPPAA33zzDcuXLycxMbHQZR3VdFevXp3U1NRCX5TCZGVlsXDhQrp3705ISEiRtiHepTIMDCrHwOAL5Rg8+26CtnwJgLXeleRcO81sKihO8YUylAuncswvJS2DLhOX23VMGWSBpY918umabneXo0drhovZ6p2HueXDdfnmf3JbK9okuvGWVWsOwXMfJGjzDACMoBByrv0Qo27P866q96J/OXbsGLGxsedNul1uXj5lyhSmTZvGLbfcckEBusowDB544AG+/vprli5det6EGyAsLIywsPyXrUJCQi74JHbHNsS7VIaBQeUYGFwpR7fei/fXD3Am4SY8hqCrXiUoNPTCtllC6b0YGFSOeWrEhjjsmLJGbBmvxeTs55+7yjElLd2WcIPZxH7E7D9IuriyX/YXUqdytMMRPmpXinbzeR8Cfd8GC7BpBhZrFqVm3QYDpjvdX4jei/7B2TJyOek+ffo07du3dzmgC3Xffffx6aefMnv2bMqUKcOBAwcAiImJISLC/970IiLiOrd2apSRBnMezpu+fCyUqeyWOEUkMPhSx5Te6NRtV+pJuwQVIMcw2J16yi+T7mId4SMoGK552+xcbfMXkHMaZtwEN3wGF13m/v2JTwtydYU77riDTz/91BOxFOrtt98mLS2NLl26EB8fb3vMmDGj2GMREZHil5KWnq9To+GztpCSll60DS58Do7vN5/X7gbNBronUBEJKPExEbSrXcGrSabbP/+clBgbRdA5d9sEWywkxPpv30gDWtdgxdNJfDa0LSueTvLshYugYLhmCjS61pzOOQ2fD4Ttizy3T/FJTtV0P/roo7bnVquVqVOnsmjRIpo0aZKvSv3VV191b4RnFKG/N3EDjU0pIr7CrTUuu5abw7oAhJaGqyfpPm4JGPruDjzeqnEu1prhYhQfE1F8xxBcCvpOBWsObP0GcjLh85vgxs+hdlLxxCBe51TSvWHDBrvpZs2aAbBlyxa7+ed2qib+TWNTiogvya1xOfdePJdrXE6fhG8fyJu+bBSU1WebBAZ9dwcmt33+FYEvNbH3W8Gl4Nr3zKbmf3wL2Rnw2Q0wcCbU6uzt6KQYOJV0//jjj56OQ3xMQc2YOtWN04etiHiF22pclrwIR3abz2u0g1a3uz1WEW/Qd3fg8naNc7HWDAeq4BC47gP4YjD8OddMvD8dADd9AYmXejs68TCXO1KTkiHQOs4QkcBwwTUue9fC6rfM58Fh0PtNCHK5exMRn6Tv7sCmGucAEBwC130IXwyCbfMgOx0+7Q83fQkJHbwdnXiQy0l33759HTYjt1gshIeHU6dOHQYOHEi9evXcEqB4hzebMYmIFKbINS7ZmTD7PuDMB1vSMIi9yK2xiXiTvrsDn2qcA0CpULj+I5h5C/w1H7JOwfTr4eYvoWbxjxAlxcPly/sxMTEsWbKE9evX25LvDRs2sGTJErKzs5kxYwZNmzbl559/dnuwUnxymzEFnynjQOk4Q0RKsOUvQ+o283l8M2j3QKGLi/ibgr67AVbuSPV4T9ciBUlJS9c5eLZSodD/f3DRmTG7s06aiXfyau/GJR7jck135cqVGThwIJMnTyboTJM8q9XKQw89RJkyZfj888+5++67eeqpp1ixYoXbA5bio2ZMIhIwDmyGFa+Zz4NKQZ/JZsc2IgHm3O/u5X/9R4dxS9SxmniNOvcrQKkw6P+xOXb39kVw+gR8ch2WgV94OzLxAJdrut9//30efvhhW8INEBQUxAMPPMDUqVOxWCzcf//9+Xo2F//kC2NTiohckJxss1m5Nduc7vgoVG7s3ZhEPCj3uxvwytjOIrm8Nb643wgJhwHTodaZocNOHyf4s+spe3KHd+MSt3M56c7OzubPP//MN//PP/8kJycHgPDwcA0fJiIivmHlG5Dym/k8rj50ety78YgUk8I6VhMpDjoHnRASDjd+Bonm0GGWzOO03/Eylv3rvRyYuJPLSfctt9zC7bffzmuvvcaKFSv4+eefee2117j99tu59dZbAVi2bBkNGzZ0e7AiIiIuSf0blo4zn1uCoM//mU36REqA3I7VzqaO1aQ46Rx0UkgE3Pg5JJhDh4XknCL4s+th/wYvBybu4vINba+99hqVKlViwoQJ/PvvvwBUqlSJRx55hKeeegqAHj160LNnT/dGKiIi4gqrFb59AHIyzem290K1Vt6NSaQYeXtsZxGdgy4IjYSBM7BOv56gPT9jyUiD/10Dg76F+Kbejk4ukMtJd3BwMM888wzPPPMMx44dAyA6OtpumRo11DmCiIh42dr3IHmV+bxcIiQ94914RLxAnaKKt+kcdEFoFDn9P+XwWz2IPbkNMo7C//rAoDnqi8TPudy8/GzR0dH5Em4RERGvO7IHFo3Km+79plmLIFICqVNU8Tadgy4IjWJ17cewVmtjTqcfgY96wwF1Uu3PnKrpbtGiBYsXL6ZcuXI0b9680E7S1q/XTf8iIuJFhgFzHjLHPQVoOQQSL/VuTCIiIk7KCQ4n54bPCfp8AOxbA+mH4X+9YdBcqNTA2+FJETiVdPfp04ewMLPjmWuuucaT8YiIiFyYjdNh54/m8+iq0P1578YjIiLiqrAycPNX8Ek/2LcWTh2Cj66GwXOh4sXejk5c5FTSPXLkSIfPRUREfMrxA/DD8LzpqyZBuP/cBpWSls6u1JMkxkapGaZ4RUpaBvvS0nQOSrHT558D4dFm4v1xX/jnVziVeibx/g7i6nk7OnGByx2pARw9epQvv/ySHTt28MQTT1C+fHnWr19PpUqVqFq1qrtjFBHxOfpx4IMMA757DDLSzOkmA6BuD+/G5IIZa5MZNmszVgOCLDC2X2MGtFbHpFJ8Vv1r4ZGJy3UOSrHT518hwmPg5lnw8TXmEGIn/8tLvGMv8nZ04iSXO1LbtGkTdevWZfz48bzyyiscPXoUgK+//pphw4a5Oz4REZ8zY20yHcYtYeC7v9Bh3BJmrE32dkgBISUtnZU7UklJSy/aBn7/Gv6caz6PjIWe49wXnIelpKXbfnACWA0YPmtL0V8LERelpGUwY2eQzkEpdvr8c0JEWbjl67yhw078C9OugtTtXg1LnOdy0v3oo48yePBg/v77b8LDw23zr7jiCpYvX+7W4EREfI1+HHjGF7/uu7ALGScPwbwn8qZ7vQyR5d0bpAftSj1pO6dy5RgGu1NPeScgKXH2HDqFgX1HuToHpTjo889JEeXglm/yhg47cQA+ugoO7fBqWOIcl5PutWvXctddd+WbX7VqVQ4cOOCWoEREfJV+HLjf0Ux4dvbWC7uQMf9p8143gPpXQcO+7g/UgxJjowg6Z2CQYIuFhFgNcybFo2aFSCzYf7jpHJTioM8/F0SWh1u/hUpnEu/jKWaN9+Gd3o1LzsvlpDs8PJxjx47lm79t2zbi4uLcEpSIiK/SjwP3+y/DcmEXMv76ATbPNJ+Hx8CVE6GQoS19UXxMBGP7NSb4TNzBFgsv9Wuk/gKk2MTHhDOgltX2+aZzUIqLPv9cFFkebp0NFRua08f3w7Sr4fAu78YlhXK5I7U+ffrw/PPPM3Om+QPHYrGQnJzM008/zbXXXuv2AEVEfEnuj4Phs7aQYxj6ceAGceEGQRbsEm+nL2RkpMGch/OmL38JylR2e4zFYUDrGnSqG8fu1FMkxEbqnJJi166Swb39OvFP2mmdg1Ks9PnnoqgKMOhbs5b7vz/g2L68ztXK1fR2dOKAy0n3K6+8Qq9evahYsSLp6el07tyZAwcO0K5dO1588UVPxCgi4lP048C9yobBC30aMGL2H65fyFj4nHmVH6B2V2h2k2eD9bD4mAidT+JV8THh1Igt4+0wpATyxc8/nx6pJCo2L/FO3QZpe817vAd/B2XV87uvcTnpjo6OZsWKFSxZsoT169djtVpp0aIFl112mSfiExHxSb7448CfXd+yGkkXV3btQsau5fDrNPN5SJQ5JrefNSsXERHf5BfDmJWuCIPmmMl26l9wNNlMwofMg5hq3o5OzuJy0n3q1CkiIyPp2rUrXbt29URMIiJSAJ++6n6BXLqQcfoUfPtg3vRlo9SkTkRE3KKgkUo61Y3zve/eMpXMxHvalXBoOxzdYz4fPA9iqno7OjnD5Y7UypYtS/v27Rk+fDgLFizg5MmTnohLxCkXPK6viB/R+OBn+fFFOHKm05ga7aD1Hd6NR3yCvhNKHpW5eILfjVRSpjIMmgvla5vTR3abtd/H9ns1LMnjck33smXLWLZsGUuXLmXy5MlkZGTQokULunTpQufOnbniiis8EadIPn7R7EfETfzqqrun7VsHq98ynweHQe83Icjla8gSYPSdUPKozMVTckcqKVIHn94SHQ+D58KHvcyL0od35nWu5qcdjAYSl3+ltGvXjqeffpr58+dz5MgRli9fTv369Zk4cSJXXXWVJ2IUyaegBERXuiVQ+d1Vd0/JzoTZ94FhNaeThkHsRd6NSbxO3wklj8pcPMlvhzGLrmIm3uUSzOlD2817vI//W+RNqjWJe7hc0w3w559/snTpUluNd1ZWFldffTWdO3d2d3wiDhWWgPj8B6JIEfjlVXdPWP4K/Pen+Ty+GbR7wKvhiPec3b+BvhNKHpW5eJrfjlQSU81saj6tl9mx2qG/z9R4zzU7XnOBWpO4j8tJd+XKlcnKyqJr16506dKF4cOH07hxY0/EJlIgJSBS0mh8cODAZljxqvk8qBT0mQzBRbp2LH7u3B+CT/Wsr++EEka/A6Q4+O1IJWWrn0m8rzSHEkvdZibeg+ZC6TinNqHb2tzL5ebllStX5sSJEyQnJ5OcnMy+ffs4ceKEJ2ITKZDfNvsRuQADWtdgxdNJfDa0LSueTipZV5tzsmH2/WDNNqc7PgqVdcG3JHL0Q3DC/G08dUV9fScUE19obqrfASLnUa6mWbsdfWbosP/+hP/1hpOpTq2u29rcy+Uqgo0bN3L06FGWL1/OsmXLGDFiBL///jtNmjQhKSmJcePGeSJOkXz8ttmPBKziGM7LW1fdvT5U2ao3IWWj+TyuPnR6vPhjCBBeL8sLVNAPwSZVy7Li6aSA/U7wlXLzpeamrvwO8MTr5ytlIlKgcglm4j3tSjj2DxzcCv/rA7d+C1EVCl1VrUncq0jt8sqWLUvv3r3p2LEjHTp0YPbs2Xz66aesW7dOSbcUK79t9iMBx5d+iLqb148t9W/4ceyZCQv0+T8oFVZ8+w8gXi9LNyjsh2Cgfif4Srm5o7mpuxNVZ8rcE6+fr5SJyHmVTzwzjvdVcHw//LsFPj6TeEeWL3A13dbmXi43L//666956KGHaNq0KRUrVuSee+7h5MmTvPbaa2zatMkTMYqI+LRA7kXX68dmtcK3D0BOpjnd9l6o1qp49h1gvF6WblLSmhX7UrldaHPTGWuT6TBuCQPf/YUO45YwY22yB6K054nXz5fKRMQpFWqf6UjtzNBhBzabNd7pRwpdrUTf1uZmLtd033XXXXTq1ImhQ4fSpUsXGjVq5Im4RET8RiD3ouv1Y1v7HiSvMp+XS4Cuz3p+nwHK62XpRiXp9iJfKrcLaW7qrU6ZPPH6+VKZiDgtN/GediWc+BcObIL/XQO3zoaIsgWuFqgtiIqby0n3wYMHPRGHiIjfCuT7nrx6bEf2wKJRedO934RQ/39NvSXQztOS8kPQl8rtQpqbeitR9cTr50tlIuKS2IvympqfPGj2lfJxX7jl60IT73OpPwPXudy8XERE7AVyc1evHZthwJyHIOukOd1yMCR28uw+A1wgn6eBzNfKrajNTXMT1bMVR6LqidfPU2XiC73CSwkQV89MvKPODB22fz18ci1kHHNqdW/cJhIINMCpiIgbBHJzV68c28ZPYeeP5vMyVaD7857fZwkQyOdpIPO1citKKwNvdsrkidfP3dtUx2xSrCrWNztS++gqOHUI/llnJt63zIKwMgWuprG7i05Jt4iImwRyc9diPbbjB+CHYXnTV0+C8Jji2XcJEMjnaSALhHLz5sUDT7x+7tqmEhnxikoN8pqapx+GfWvgk+vg5i8LTLzVn0HROdW8fNOmTVitVk/HIiIiJZ1hwHePQUaaOd24P9S93LsxiYjbxMdE0K52Bf1AP8uF9govUmSVGsKgbyGinDm9dzVM7w+ZJxwu7q3bRAKBU0l38+bNSU1NBaBWrVocOnTIo0GJiEgJtfUb+HOu+TwyFnqO82o4IiKepkRGvKpyY7MH8/Cy5nTySvi0P5w+mW9RX+tjwp84lXSXLVuWXbt2AbB7927VeouIiPudOgzznsib7vUyRFXwXjwiIsVAiYx4XXzTM4n3mVu59vwMnw6A0/lbW2js7qJx6p7ua6+9ls6dOxMfH4/FYqFVq1YEBwc7XHbnzp1uDVBEREqI+U/Dyf/M5/WvgoZ9vRuPiEgx8bXO8qQEqtIMbvnGHLs7Mw12/wSf3QADZ0CI/fkYCH1MFDenku6pU6fSr18/tm/fzoMPPsjQoUMpU6bgnu1ERERcYdm+EDbNMCfCY+DKiWCxFL6SiEgAUSIjXle1hTlm98fXQOYx2LUMPrsRbvwsX+ItrnG69/KePXsC8Ouvv/LQQw8p6RYREbcolZNO8Lyn82Zc/hKUqey9gEREREqqai3h5lnwcV84fdwcvvPzm+CGTyEk3NvR+S2n7uk+24cffmhLuPft28c///zj9qBERKTkaPDP51iO7zcnaneFZjd5NyARsUlJS2fljlRS0tK9HUqJo9devKZ6a7j5KwgtbU7vWAwzbobsTO/G5cdcTrqtVivPP/88MTEx1KxZkxo1alC2bFnGjBmjDtZERMQllj0rSDz0ozkREgVXTVKzchEfMWNtMh3GLWHgu7/QYdwSZqxN9nZITgmEZNVfX3sJIDXawE1fmt/NANsXwoxblHgXkctJ9zPPPMPkyZMZN24cGzZsYP369bz00ku8+eabjBgxwhMxiohIIfz2B+bpUwR/90je9GWjoFxNr4XjTxyVubPzvMWXYpHzS0lLZ9iszbbxo60GDJ+1xefLLxCS1eJ87fW+dI+AfR1rtoObvoCQM8PX/f0DzBwE2afzLRqwr4GbOH1Pd66PPvqI9957j969e9vmNW3alKpVq3Lvvffy4osvujVAEREp2Iy1ybYfZ0EWGNuvsf8M3/Hji1iOmMNRWqu1Iaj1HV4OqHApaensSj1JYmyUVzs7clTmgFPzvHVu+PV5WkLtSj1pS/py5RgGu1NP+WxnXwUlq53qxvlszI4U12uv96V7BPzrmNABBs6E6ddDdjr89T18OQSunwbBIUAJeA3cwOWa7sOHD1O/fv188+vXr8/hw4fdEpSIiJyfv9ZEAbBvHax+C4AcSwg5V02CIJe/koqNr9SeOSrzYV9tdmqet84Nvz5PS7DE2CiCzrnTI9hiISE20jsBOaGwZNWfFMdrr/ele5SY1zHxUnPosFJnOlL7c66ZeOdklZzX4AK5/AunadOmTJ48Od/8yZMn07RpU7cEJSIi5+e3PzCzM2H2fWCY/YBsq3wNVLjIuzEVwpd+UDgqc+uZmM43z1vnht+epyVcfEwEY/s1JvhMHwvBFgsv9Wvk0zXG/nihwJHieO31vnQPT72ORzNh9c7DvpW41uoMN36el3j/MQe+up3dB4/qXHKCy83LJ0yYwJVXXsmiRYto164dFouFlStXsnfvXubNm+eJGEVExIHcH5hnf9n5xQ/MnybCf38CYFRuwvZKvfDdlNu3mtk6KvMgACfmeevc8NvzVBjQugad6saxO/UUCbGRPp1wQ16yOnzWFnIMwy8uFBTE06+93pfu4YnX8Ytf9zFqfTDG+nW+11S7dpI5dNhnN0JOJmydTbNsgxDL9WQZwbbFdC7l53JNd+fOnfnrr7/o27cvR48e5fDhw/Tr149t27Zx6aWXeiJGERFxwB9rojiwxUy6AYJKkX3VGxiW4MLX8TJfqj1zVOZjr23s1DxvnRt+eZ6KTXxMBO1qV/Cb8hrQugYrnk7is6FtWfF0ku8kK0Xgydde70v3cPfrmJKWzrOzt2Jgbs8nm2rX6WYm3sGhAET89S2LEj8l1GK2XtO55JjLNd0AVapUUYdpIiI+wK9qonKyzWbl1mxzuuMjUKkR4Nu9C/ta7VlBZe7sPF+KWcQT4mMidI45Qe9L93Dn6+hLLasKddFlMGA6zLgJck5Tc//3/NY0ko0tx5FQsYxvxeojipR0i4iI7/CbH5irJkPKRvN5XH3o9AQYha7hM3ztx6mjMnd2nrf4UiwiYtL70j3c9Tr6VbP/uj2g/8cw42awZhHx51e0Cw2Ba97ydmQ+yXe7ihURkcCRuh2Wjj0zYYHek6FUmFdDcpW/NbMVERH/Eh8TwQt9GmA5c0Xa2y2rzqteT+j/EQSdqcfd9Dl8+wBYrd6NyweppltERDzLaoVv74fsDHO67b1QvbV3YxIREfFB17esRlbyJmo3a0vtStG+m3Dnqn+lOWb3F4PN28c2TgeLBa5+06eHAi1ueiVERMSz1r0PyavM5+USoOuzXg1HRETEl5UNgzaJ5X0/4c518dVw3QeQ2zHqhk9g7sOq8T5LkZLu7OxsFi1axDvvvMPx48cB2L9/PydOnHBrcCISYI7tN5sZ60O45DiaDItG5U33fhNCffDeNBERESm6Bn3guvfzEu/1H8G8x8Dwk85bPMzl5uV79uyhZ8+eJCcnk5mZSffu3SlTpgwTJkwgIyODKVOmeCJOEfFHWRmQvBK2L4YdS+DgVnN+lRZw+UtQs53bdpWSls6u1JMkxkb5z5XhQGcYMOchOH3mgmzLwZDYyashiYiIiIc07AuGFb66w/y77gOwBEGvV8wm5yWYy0n3Qw89RKtWrfjtt9+oUKGCbX7fvn2544473BqciPgZw4DUv2HHYjPR3r0Csh2MLbl/PXzY02yOdNloqFD7gnY7Y20yw2ZtxmpAkAXG9mvs12OzBoyNn5oXWwDKVIHuz3s3HhEREfGsRteaLRq/vtNMvNe+ZybeV0wo0Ym3y0n3ihUr+PnnnwkNDbWbX7NmTf755x+3BSYifiIjDXYuO5NoL4G0gsZctkDVlnD6JPz3hznrjzmwbT5cMtQcPiqyvMu7T0nLsCXcYA6zMXzWFjrVjVONtzcdPwA/DMubvnoShMd4LRwREREpJk2uNxPur+8CDFgz1Wx23nNsiU28XU66rVYrOTk5+ebv27ePMmXKuCWowrz11lu8/PLLpKSk0LBhQyZNmsSll17q8f2KyBlWqznW8vbFZqK9dw0Y+T8TACgTD7W7QZ1uUKuLmVRbc8wONn58EU78C9YsWP2W2dtlpyfNBNyFoaT2HDplN54lQI5hsDv1lJJubzEM+O4x84IMQOP+UPdy78YkIiIixafpADPx/uYewIBf3oagYOjxQolMvF3uSK179+5MmjTJNm2xWDhx4gQjR46kV69e7owtnxkzZvDwww/zzDPPsGHDBi699FKuuOIKkpMLqlkTV6SkpbNyRyopaQ6aA7uwjCvLubr+hW7XXXwljgvl9HEc/xc2fgZf3g6v1IF3k+DHF8weqc9OuINDzeS6+xi4ZxU8+gdc83/QqF9eLXZQMLQcBA+sN5PsUmcS44w0WPAM/N8l8Ps3Tne8UbNCJEHnfHYHWywkxObvrOvc43V0/M7OczdfOucveJ9bv4E/55rPI2Oh57hC9pFR9ECdEAjv1UA4BilcoJZxoB6XiLsE+nskJfEatrcbh8GZH2qrJsPC55z6jRdor43LNd2vvfYaSUlJNGjQgIyMDAYOHMjff/9NbGwsn332mSditHn11Ve5/fbbbfeOT5o0iR9++IG3336bsWPHenTfgc6Ze2KdvW/2Qu+vLWh9X7lv11fiuFCFHkf2adi7GrYvMpuM/7u54A1VuMisya7dDRI6QGiUcwGElYauz5iday15AX77DDDgyG74YhBUbwuXvwjVWhW6mfiYcMb2a8zwWVvIMQyCLRZe6tcoXy33ucfbt3lVvt7wj93xA/leE0fz3F3evnTOX/A+Tx2GeU/kTfeaAFEV7BY5dx/9Ey144pJtILxXA+EYpHCBWsaBelwi7hLo75G846vODcF3MC7kXfMfK98w7/G+bFSBNd6B+NpYDMP1ftzT09P5/PPP+fXXX7FarbRo0YKbbrqJiAjPNeU8ffo0kZGRfPHFF/Tt29c2/6GHHmLjxo0sW7Ys3zqZmZlkZmbapo8dO0b16tVJTU0lOjq6SHFkZWWxcOFCunfvTkhISJG24WtS0jLoMnG5XRPdIAssfawT8THhTi/jynKuxjLzzjb0n/pLkbd7tgspwws9Pl/h6DgSLQf4qvspyu5fjmXPz1iyTjpc1wgtjZHQCaN2V6y1ukJZN30IHthE8KLnCNqzwm62tcE15CSNgLI17eafW44paRkkHz5FjfKR+crC0fGey4L52X/2Mo7mubu8i+Ocv9BYXNln8Lf3ErR5JgDWuleQc93/7L5UHe3DgsHih9tTvYL7blEKhPeqvxxDIH4vFhdfKmN3lqMvHVdJo/ejfyjsPRIbGez3Zejo+AYGL+alkPdt0zntH8HaZXi+xNvfPj+OHTtGbGwsaWlpheaXLtd0L1++nPbt2zNkyBCGDBlim5+dnc3y5cvp1Mkzw8GkpqaSk5NDpUqV7OZXqlSJAwcOOFxn7NixjB49Ot/8BQsWEBl5YePELly48ILW9yV/p1mwGsF286wGzJz3IxfFGE4v48pyrsbyyfyVF7RdR4pShhd6fL7i7zQLEcZp2gVtpVPQJjoH/UbNoIOw3PHyRyMS+De6CQejG3MkqjaGpRQcAA5sAba4L7ByQ6kU3JqG/3xOmcwUAIK2foPxx1x2xvXgr0pXk13Kvib93HI8BGw4Z7OOyu1cBvlbOzma5+7yLs5zvqixOLvPimm/0W6nmXBnBUeyJLQnGd9/f959GFiYtWCFW48rEN6r/nYMgfS9WFx8sYzdUY6+eFwljd6Pvs2Z94g/l6Gj4/s0pxtJ8dl0T/0IgP+2/MiaU00wLMHnXdeXPz9OnTrl1HIuJ91JSUmkpKRQsWJFu/lpaWkkJSU57GTNnSznXA0xDCPfvFzDhg3j0UcftU3n1nT36NFDNd1nSUnL4K0/8l9R6t8rya6m+3zLuLKcq7Hc3LM93zqo9XN2u2e70JruCzk+rzKs8O8Wgnb+SNbBRTwYtoZQi+P3qxFVEaNWF6y1umIkdiEqKpZaQK1iCfRKyHmSnA0fE/TTeCynDhFsZHPRwXnUOb4a66VPYm0xiCwrTpejo3I7lys13e4s7+I45y80Fqf2mXmcUu88bZu09BxL12Y3ObUPCwb9enR0e023375Xz/CXYwjE78Xi4ktl7O6abl85rpJG70f/UNh7JFBquh0dX70bXiTn7wZYdi2lQt/3uMJBx7n+9vlx7Ngxp5ZzOekuKMk9dOgQUVFO3stZBLGxsQQHB+er1T548GC+2u9cYWFhhIXlL8yQkJALPondsQ1fUSM2xOE9sTViy7i0jCvLuRpLq8TYC9quI0Upwws9vmJ3MtUcJ3n7YvPvyYMABAOc9TY+bQRzpEILKjXvBXUuw1KpEZagINd7WnSXkBBodxc0vwFWvAar3oKcTCzphwle8DTBv76PpavZEYcz5eio3K5pXoVvNuy3K0cgX9k6mufO8i7Oc76osTi1zx/GwPH95vNaSZRqNcjhvVrn7sO8p9tK9Qpl3PqZ6nfvVQf87RgC6XuxuPhiGbujHH3xuEoavR99W2HvkaysLMC/y7DQz4DYu6HtXQQVUGnqb58fzpaR0/d09+vXD4DZs2fTs2dPu2Q2JyeHTZs2Ua9ePebPn1+EcJ3Tpk0bWrZsyVtvvWWb16BBA/r06eNUR2rHjh0jJibmvG3uC5OVlcW8efPo1auX374RCpKSls7u1FMkxEYWONSSM8u4spyr61/odsE9ZeiOODwiJwv2rT3TAdpiSPkNs5G0A+USOVmjM3vLtadsw65Ujosr1lBdcjQZFj8Pm7+wm/1f6Ysp2/9NQmq0dmoz55abo3J0dp67efKcd1csBdr1E3x0lfk8JAruXQXlaha6Su4+qsaEsuHnJR77TPXZ96oLfP0YAvl7sbj4Qhl7ohx94bhKGr0f/Yuj90ggleGFfAb4y+eHs/ml0zXdMTExgFnTXaZMGbtO00JDQ2nbti1Dhw69gJDP79FHH+WWW26hVatWtGvXjqlTp5KcnMzdd9/t0f2WFPExEec9qZ1ZxpXlXF3/QrfrLr4SBwBH9pjjZW9fDLuWQ2YBzVxCoiCx05mexrtChdpEAfWLNdgiKlsDrn0P2txjDiuWvAqAuBN/wAeXQZMboNsIiKlW6GbOLTdH5ejsPHfzpXPepX2ePgXfPpA3fdnI8ybcZ+8jKysr3z347uRT79UiCoRjkMIFahkH6nGJuEugv0cu5PgC7bVxOun+8MMPAUhISODxxx/3aFPyggwYMIBDhw7x/PPPk5KSQqNGjZg3bx41a57/B55IwDh9Enb/fCbRXgSHthe8bKXGZpJdp5s5BFep0OKL0xOqtYQh38MfczAWPoflyC5z/qbPzbGh290HHR+BMN9sghSQfnwRcsuhelto7dmLryIiIiL+xuV7ukeOHOmJOJx27733cu+993o1BpFiZRhwcOuZ+7IXw56VkHPa8bKRFaBWEtS5zKzNLuO4vwO/ZrFAg95k1+rGnx8/QaND32HJOArZGfDTRFj/P0gaDs1vhWCXP+LEFft+hdVnbvcJDoM+kyHIaz0BiIiIiPikIv0i/fLLL5k5cybJycmcPm3/43/9+vVuCUykRDt1GHb+CNuXmIn28RTHy1mCofolUPtMbXZ8s5KT9ASHsrPi5dQfMJqQVZPgl3fAmgUn/4O5j5jTPV4wL0AU0FmHXIDsTJh9n9krPkCXpyH2Iu/GJCIiIuKDXE6633jjDZ555hkGDRrE7NmzGTJkCDt27GDt2rXcd999nohRJPDlZMM/v+bdm/3PrxTYAVpMDajT1UwmEztBeEyxhupzIsrC5S9C69th0SjYOtuc/9+fMP06qNXFTL4rN/ZikAHop4nw3x/m8/im0P5B78YjIiIi4qNcTrrfeustpk6dyo033shHH33Ek08+Sa1atXjuuec4fPiwJ2IUCUxp+/KajO9cChlpjpcrFQEJHc/cm30ZVKijmltHyteC/v+D5F/Mztb2rTXn71wKUy6F5jdB0rMQHe/VMAPCgS1m0g0QVAr6/J+a8ouIiIgUwOVfScnJybRv3x6AiIgIjh8/DsAtt9xC27ZtmTx5snsjFAkUWenm/di5ifZ/fxa8bMUG5j3ZdbpBjfYQEl58cfq7Gm3g9oXw+yyz5vtoMmDAhk9gyyzo8BC0fwBCi78zyICQk202K7dmm9MdH1ErAhEREZFCuJx0V65cmUOHDlGzZk1q1qzJ6tWradq0Kbt27cLJIb9FSgbDgNS/8sbM3vOz2dmXI+FloXaSeW927a4QU7VYQw04Fgs0uhbqXQlr3oHlEyEzDbJOwdKx8Os06PosNL0RgoK9Ha1/WTUZUjaaz+PqQ6cnvBqOiIiIiK9zOenu2rUrc+bMoUWLFtx+++088sgjfPnll6xbt45+/fp5IkYR/5F+FHYtO5NoL4Fj+xwvZwmCqq3OjJndDaq2UPLnCSHhZs12s5th2ThY94FZQ3s8xaytXT0FeowxL3jI+aVuNy9aAGCB3pOhVJhXQxIRERH/kZKWzq7UkyTGRgXUONzn43LSPXXqVKxWs7fau+++m/Lly7NixQquvvpq7r77brcHKOLTrDmwf2PemNn71oGR43jZMlXO6gCtM0SWL9ZQS7SoCtDrZbjkTlg4ErZ9Z87/dzN8fA1c1AO6j4GK9b0apk+zWuHbB/Jaa7S9B6q39m5M4rNK6o8qERF/4ehz2tOf3TPWJjNs1masBgRZYGy/xgxoXaNEfGe4nHQHBQURdNaQRP3796d///4A/PPPP1StqmaxEuCOpcCOM0N57VgC6UccLxccBjXb53WAFldfHaB5W+xFcOOnsOsns7O1lN/M+X8vMG8BaDkIugyD0hW9G6cvWvc+JK80n5dLMJvnizhQ0I8qkUBXEhIHCQyOPqcBj352p6Sl27YPYDVg+KwtHE3PYvz3fwb8d4Zbups9cOAAL774Iu+99x7p6enu2KSI78jOhORVZzpAWwL/bil42di6eWNm1+wAoZHFF6c4L/FSGLoUNn8Bi0fDsX/MFgrrPoBNX0DHh6HdfRCiH02A2RndolF501e/oY7oxKGCflR1qhunJEQCmi42ib9w9Dk97KvNYMGjn927Uk/atp8rxzAY9/2fGCXgOyPo/IuYjh49yk033URcXBxVqlThjTfewGq18txzz1GrVi1Wr17NBx984MlY5Swpaems3JFKSpoucridYZj3rv7yDkzvD+MT4H99YOUb+RPusGi4+Gq4ahI8vBnuXwtXjIOLuvtUwq3zxYGgIGg6AB74FbqOgNDS5vzTx2HJGHizFfw2w2xWXZIZBsx5GE6fMKdbDIJanb0akviugn5U7U495Z2A/JQ+s/1LQRebVH7iixx9TlvB45/dibFRBJ3T4DMIOLcf7kD9znC6pnv48OEsX76cQYMGMX/+fB555BHmz59PRkYG33//PZ0760eYW3xzHxzcanaqZQk+8zfIbvqfY6fZsv8k2Vg4RBBNqpenZmy0w2Xtt1Gq0O06nl/KhWXP+uvKsrnzg0qd8z+nrwlduIxjsGv5mXuzF8PRPQUsaIEqzfM6QKvWCoJDii/OItDV9/MIiYBOj0OLW+HHl2D9R2BYzU7wvr4TfnkberwICR28Hal3/PaZ+b4As1+CHmO8G4/4nJS0DP5Os5CSlmH7UXX2j7dgi4WEWN+5COnr9Jntfwq72BRotXXi/xx9TgeBXU03uP+zOz4mgrH9GjN81hZyDINgi4Une9Zj/Pw/S8R3htNJ93fffceHH37IZZddxr333kudOnWoW7cukyZN8mB4JdB/f8D+9YUuUhWoenZH1/vPPAJRvgT9TDJuS+5L5Z93niQ/GAttDx0h+PP/mQlzxlHYtzZv3OFzla6U12S8VpLZKZefUFNPF5SuCFdPgjZ3wYIRsH2hOX//BpjWC+pfBZeNhtg6Xg2zWB3/F+YPy5u+6jUIj/FePOJz8hLEYN76Yzlj+zXO96PqpX6N9HnjJH1m+yddbBJ/4ij5falfIwCPf3YPaF2DTnXj2J16ioTYSOJjIigbGVIivjOcTrr3799PgwYNAKhVqxbh4eHccccdHgus5LKceWjMc8C8zzangN7AiygIqARwfHMBC4RAjbZm52d1ukGlRn7bAZquvhdBxYvh5i/N+/cXjMi7peDPufDXfGh1O3R+yq8uvhTZvMfMi1IAja+Hej29Go74loISxBVPJ7Hi6SS7H1XiHH1m+6eCkhiVmfgqR8kv4HCeu8XHRNhtu6BYAo3TSbfVaiUkJK8ZbXBwMFFR6kjH7YaeacZpGOZwVEaO3d8DaSfo/fpyMKwEYyXYYiUEgxl3taZiVIiDdaxmDe65887Zbr75dutYC1/WOLO808ue9dfpZa3n7KegZR3s33Dintzytc7UZl8GCR0hzLy/NyUtnV07D/ltT6S6+n4BaneFuzrDxk9hyQtw4oB5/q15B3773GyS3uauwB2n+vdv4I855vPIWOg53qvhiO8pLEFsV7uCX35meps+s/1XSUkcJHCcm/wWNM9bsQQap5NuwzAYPHgwYWHmD8yMjAzuvvvufIn3rFmz3BthSWWxQHApzi2iyhFleaxfx7yrqZhXUysm6H6vAp1zASPrdAYLfphPj8u6ERJ8pgm6gzGzA+G+Ol19v0BBwdDiFmjYF1a+aXaml3UKMtNg4QhY+y5cNgoa9vPb1hAOnToM8x7Pm+41oWTU7ItLlCC6nz6z/VtJSBxEpGicTroHDRpkN33zzTe7PRhxjq6muujcCxiWELJLRZmJdojjTtAC6b46nS9uEFYakoZBy8Hw4wuwYTpgmENpfXkbrD7T2VqNNt6O1D3mD4OT/5nP611pXlQQOUdugnj2xUkliBdOn9kiIoHH6aT7ww8/9GQc4iJdTfWsQLuvTueLm0THQ5//gzZ3w4JnYedSc/6+tfBBD2hwjVnzXT7Ri0FeoL8WwKbPzedhMXDlxMCqxRe3GtC6Bu0SyzFz3o/075VEjdgy3g4pIOgzW0QksBTjmEwi3uPqmKeOxhJUs0mxqdwYbvkGBn4BsfXy5m/9Bv7vEvjhGUg/4q3oii7jGMx9OG/68hfNCw0ihYiPCeeiGIP4mHBvhyIiUuK4+htXvMPpmm7xDylp6exKPem3HX95wrn3Zr/QpwHn6wJQ99XJeVksULeH2eHa+o/MMb5PpULOaVg1GTZON3s5b3U7lAr1drTOWTQSjv1jPq+VBM11G5GII/quFRFfEAj9D5UUSroDiN54+Tm6N/vZ2VsZ2fz86+q+OnFKcClofbs5pNbPk2DV/0F2hlnTPf9pWDMVuj9vjvPty820d/0E6z4wn4dEwdWv+3a8Il6i71oR8QWB1P9QSaDm5QGioDdeSW9q4ujebKsB/2U4l0zEx0Ro6BtxTng0dHsO7l8HTQbkzT+8E2bcDB/2gn9+9V58hTl9CuY8mDd92UgoV9N78Yj4KH3XioivKKz/IfE9SroDhN54jjm6NzvIAnHhhuMVRC5U2erQbyoM/RFqdsibn7wS3u0KX91h9nruS5a+ZF4cAKjeBloP9W48Ij5K37Ui4ivU/5B/UdIdIPTGcyz33uzgM81kgy0WXujTgLJhXg5MAl/VFjD4O7jhU6hQJ2/+5i/gzVawaBRkpBV7WPk6XPnnV7NJPEBwGPSeDEH6ahBxRN+1IuIrHP3GLaj/IXW25n26pztAqOOvgp17b3ZsZCnmzdvk7bCkJLBYoP6VcFEP837ppeMg/TDkZMKK12D9x9DlaWg55MxY8p517r2o46+pz/W/3g+G1Vygy1MQV9fjcQQSdahVsui7VkR8iTP9D6kfCt+gpDuAqOOvgp095mlWVpaXo5ESJzgE2txl3uv90yvwyztmL+enUmHe42c6WxsDdS/3WOdlju5F3T/nRSi11ZxRuQm0f7DgDUg++iFTMum7VkR8ydm/cc+lztZ8h9oQBhh1/CXiwyLKQo8X4P610LBf3vzUv+CzAfC/3pDimVYY596LWs+SzL3B35gTQaWgz/+ZFwfEKepQq2TTd62I+AP1Q+E7lHSLiBS3cglw/Ydw+yKodkne/F3L4Z1O8M29cGy/W3d59r2oweQwIWQqIZYcc0aHhyG+iVv3F+j0Q0ZERHyd+qHwHUq6RUS8pXpruH0BXD8NyuYO0WXAxunwRgtY8iJknnDLrs7ucOX24Hk0DTrTW3lsPej8pFv2UZLoh4yIiPg6VzpbE8/SPd0iIt5ksUDDvlCvF6x5F5ZPMHs1z043n6//CJKegeY3Q1DwBe1qQOsaJMUdI/bjWZADYIE+k6FUyenO310dn6lDLRER8Qfqh8I3KOkWEfEFpcKg/f3QbCAsmwBr3wVrNpz4F+Y8CL9MgR5joM5lRd+H1UrFJU+YvacDtL0Hql9S+DoBxN0dn+mHjIiIFJcLuWhcWGdrUjzUvFxExJdElocrxsF9a6D+VXnzD26FT66Fj/vBv1uLtu1170PySvN5uQTo+uwFh+svPNXxmTrUEhERT5uxNpkO45Yw8N1f6DBuCTPWJns7JHGRkm4RCTgpaems3JHq3z1JV6gNN0yHwfOgSvO8+TsWw5QO8O2DcPxf57d3NBkWjcqbvvoNCI1yW7i+Th2fiYiIP9JoGYFBSbeIBJSAuxqc0AHuWAL93oOY6uY8w2re6/1mC1j2Mpw+T+JoGDDnYTh9plO2FoOgVmePhu1r1PGZiIj4I100DgxKukUkYATs1eCgIGhyvTm+d7eREFrGnH/6BPz4ArzZEjZ+Clar4/V/+8ysIQcoU8W8N7yEUQ+uIiLij3TRODAo6RaRgBHwV4NDIuDSR+HBDdD6DrCc6c38+H745h6Y2tkc6/tsx/+F+cPypq96DcJjii9mHzKgdQ1WPJ3EZ0PbsuLppAvqRE1ERKQ46KJxYFDv5SISMHKvBp+deAfk1eDScXDlRLjkTlj4HPw135x/YBN8dDXUvQK6Pw9xdWHe45Bx1Px/4+uhXk+vhe0L1IOriIj4G42W4f9U0y0iAaPEXQ2OqwcDZ8Cts6Fy47z5f30Pb7WFT2+AP74150VWgJ7jvBOniIibBERHmSJFoNEy/JtqukUkoJTIq8G1usCdy+C3z2HJGDieAkaOmXznumICRMV6LUQRkQs1Y22yrd+OIAuM7ddYt4kEuAsZm1rElyjpFpGAUyKbEAcFQ/OboOE1sOr/YMUkyDpp/q9eL2h0rTejExG5IAV1lNmpblzJ+7z3Ye5MknWRRQKJkm4RkUASGgWdn4QWt8Kq/2/v/mOrLO8+jn9OS1va0p5RjvSHgMXJCrMotiAgG4U9G5BIgLFloEgkGiJLwM7FLeCyAfrMMscw2RSd/NwMGyzKErIlBhwCOrrCCs7WBrZiawXaYRm0SKEt9vv8wdMzTn8i9uac+z7vV9KEc9/Xuc91zudcLd/71/WC1PyJ9D8/kXy+3p8LABGqpxtlUnRHhr4sktnJAq+h6AYAL0rJkKb973U15fQ9AJEuam6U6VJ9XSSzkwVew43UACCK7Thco0lr9urBDSWatGavdhyuCXeXAKCTqLtRpsv09ZSdzE0Nr+FINwBEKU7fA+AmUXmjTJfo6zMR2neyPLWzXJ+asZMFrkfRDQBRitP3ALhNVN4o0wWcKJLZyQIvoehGn+CaUG8jX2/iGkkAQF9xokhmJwu8gqIbnxtTOngb+XoXp+8BAPoSRTLQNYpufC5cE+pt5Ot9nL4HAADgLO5ejs+lr+9WichCvtEh05+oiV8cRMENAIDL1TZc0sET9aptuBTuruAaHOnG58I1oZ+N266NJl8AAAB34JLAyMWRbnwuzJt5/dw4HzL5AgAARL7uLgnkiHdk4Eg3PjeuCe2dm6+NJl8AAIDIxjSgkY2iG32Cu1X2zO2/CMkXAAAgcnFJYGTj9HLgJmj/RXgtfhECAACgL3BJYGTjSDdwEzAfMgAAAJzEJYGRi6IbuEn4RQgAAPD5uG0mmJuNSwIjE0U3cBPxixAAAODGMCUW3IprugEAAABENKbEgptRdAMAAACIaD3NBANEOopuAAAAABGNmWDgZhTdAAAAACIaU2LBzVxxI7Xq6mo988wz2rt3r+rq6pSVlaWHHnpIP/rRjxQfHx/u7gEAAABwGDPBwK1cUXQfO3ZMbW1t+vWvf6077rhD5eXlWrx4sS5evKi1a9eGu3sAAAAAbgJmgoEbuaLonjFjhmbMmBF8fPvtt+v48eN66aWXKLoBAAAAABHLtdd0NzQ0KC0tLdzdAAAAAACgW6440t3RiRMn9Ktf/Uq/+MUvemzX3Nys5ubm4OPGxkZJUmtrq1pbW2/otdufd6PPR/iRoTeQozeQo/uRoTeQozeQo/uRobtcb04+M7Pemzlj1apVWr16dY9tDh8+rLFjxwYfnz59WgUFBSooKNDGjRtvaPu/+93vlJTE9AIAAAAAgBvT1NSkBx98UA0NDUpNTe22XViL7vr6etXX1/fYJjs7W/3795d0teCeOnWqxo8fr61btyompuez47s60j106FDV19f3+KH0pLW1VXv27NE3vvENxcXF3dA2EF5k6A3k6A3k6H5k6A3k6A3k6H5k6C6NjY0KBAK9Ft1hPb08EAgoEAhcV9tTp05p6tSpys/P15YtW3otuCUpISFBCQkJnZbHxcV97i9xX2wD4UWG3kCO3kCO7keG3kCO3kCO7keG7nC9Gbnimu7Tp09rypQpGjZsmNauXauPP/44uC4jIyOMPQMAAAAAoHuuKLp3796tyspKVVZWasiQISHrwnh2PAD0qrbhkqrqL2p4IDkq5hWNtvcLAADQG1cU3YsWLdKiRYvC3Q0A+Ex2HK7Rip1lajMpxicVzR2teeOGhbtbjom29wsAQHfYCY1ruXaebgCIZLUNl4IFqCS1mfTUznLVNlwKb8ccEm3vFwCA7uw4XKNJa/bqwQ0lmrRmr3Ycrgl3lxBmFN0A4ICq+ovBArTdp2aqrm8KT4ccFm3vFwCArrATGl2h6AYABwwPJCvGF7os1udTdiApPB1yWLS9XwAAusJOaHSFohsAHJDpT1TR3NGK9V2tRGN9Pj07N9ez13VF2/sFAKAr7IRGV1xxIzUAcKN544Zp8pduUXV9k7IDSZ4vQKPt/QIA0FH7TuindpbrUzN2QkMSRTcAOCrTnxhVf2ij7f0CANARO6HREUU3AAAAAPQhdkLjWlzTDQAAAACAQyi6AQAAAABwCEU3AAAAAAAOoegGAAAAPK624ZIOnqhXbcOlcHcFiDrcSA0AAADwsB2Ha7RiZ5naTIrxSUVzR2veuGHh7hYQNTjSDQAAAHhUbcOlYMEtSW0mPbWznCPewE1E0Q0AAAB4VFX9xWDB3e5TM1XXN4WnQ0AUougGAAAAPGp4IFkxvtBlsT6fsgNJ4ekQEIUougEAAACPyvQnqmjuaMX6rlbesT6fnp2bq0x/Yph7BkQPbqQGAAAAeNi8ccM0+Uu3qLq+SdmBJApu4Caj6AYAAAA8LtOfSLENhAmnlwMAAAAA4BCKbgAAAAAAHELRDQAAAACAQyi6AQAAAABwCEU3AAAAAAAOoegGAAAAAMAhFN0AAAAAADiEohsAAAAAAIdQdAMAAAAA4BCKbgAAAAAAHELRDQAAAACAQyi6AQAAAABwCEU3AMBRtQ2XdPBEvWobLoW7KwAAADddv3B3AADgXTsO12jFzjK1mRTjk4rmjta8ccPC3S0AAICbhiPdAABH1DZcChbcktRm0lM7yzniDQAAogpFNwDAEVX1F4MFd7tPzVRd3xSeDgEAAIQBRTcAwBHDA8mK8YUui/X5lB1ICk+HAAAAwoCiGwDgiEx/oormjlas72rlHevz6dm5ucr0J4a5ZwAAADcPN1IDADhm3rhhmvylW1Rd36TsQBIFNwAAiDoU3QAAR2X6Eym2AQBA1OL0cgAAAAAAHELRDQAAAACAQyi6AQAAAABwCEU3AAAAAAAOoegGAAAAAMAhFN0AAAAAADiEohsAAAAAAIdQdAMAAAAA4BCKbgAAAAAAHELRDQAAAACAQyi6AQAAAABwCEU3AAAAAAAOoegGAAAAAMAh/cLdgZvJzCRJjY2NN7yN1tZWNTU1qbGxUXFxcX3VNdxEZOgN5OgN5Oh+ZOgN5OgN5Oh+ZOgu7XVle53Znagqui9cuCBJGjp0aJh7AgAAAADwggsXLsjv93e73me9leUe0tbWptOnTyslJUU+n++GttHY2KihQ4fqo48+Umpqah/3EDcDGXoDOXoDObofGXoDOXoDObofGbqLmenChQvKyspSTEz3V25H1ZHumJgYDRkypE+2lZqaykBwOTL0BnL0BnJ0PzL0BnL0BnJ0PzJ0j56OcLfjRmoAAAAAADiEohsAAAAAAIdQdH9GCQkJWrlypRISEsLdFdwgMvQGcvQGcnQ/MvQGcvQGcnQ/MvSmqLqRGgAAAAAANxNHugEAAAAAcAhFNwAAAAAADqHoBgAAAADAIa4tuouKijRu3DilpKRo8ODBmjNnjo4fPx7SZtWqVRo5cqSSk5M1cOBAff3rX1dJSUmP262urpbP5+v088Ybb4S0e/HFFzVq1CglJiYqJydHv/3tb3vt85QpUzptd/78+SFtzp07p4ULF8rv98vv92vhwoU6f/789X0oLuTVHLOzszu1Wb58+XV+Ku61fv16DR8+XP3791d+fr7efvvt4Doz06pVq5SVlaXExERNmTJF77//fq/bLCsrU0FBgRITE3Xrrbfq6aefVsdbUezfv1/5+fnq37+/br/9dr388su9bre5uVnLli1TIBBQcnKyZs2apZMnT4a0ibbxKHkvQ8Zi5Of4yiuvaMqUKUpNTZXP5+tyjEXjWJS8l2M0jke3ZPif//xHy5YtU05OjpKSkjRs2DA9/vjjamhoCGnHWPRGjtE4FsPOXGr69Om2ZcsWKy8vt3fffdfuv/9+GzZsmH3yySfBNtu2bbM9e/bYiRMnrLy83B599FFLTU21M2fOdLvdqqoqk2Rvvvmm1dbWBn+am5uDbdavX28pKSm2fft2O3HihP3+97+3AQMG2K5du3rsc0FBgS1evDhku+fPnw9pM2PGDMvNzbWDBw/awYMHLTc312bOnHmDn1Lk82qOt912mz399NMhbS5cuHCDn5I7bN++3eLi4mzDhg1WUVFhhYWFlpycbB9++KGZma1Zs8ZSUlLs9ddft7KyMps3b55lZmZaY2Njt9tsaGiw9PR0mz9/vpWVldnrr79uKSkptnbt2mCbDz74wJKSkqywsNAqKipsw4YNFhcXZ6+99lqP/V2yZIndeuuttmfPHjty5IhNnTrV7r77brty5UqwTbSNRy9myFiM/Byff/55KyoqsqKiIpNk586d69Qm2saimTdzjLbx6KYMy8rKbO7cubZr1y6rrKy0v/zlLzZixAj71re+FdKOseiNHKNtLEYC1xbdHZ05c8Yk2f79+7tt09DQECzEutNerB09erTbNhMnTrQnn3wyZFlhYaFNmjSpxz4WFBRYYWFht+srKipMkv3tb38LLisuLjZJduzYsR637RVeyNHs6i+z559/vsc2XnPvvffakiVLQpaNHDnSli9fbm1tbZaRkWFr1qwJrrt8+bL5/X57+eWXu93m+vXrze/32+XLl4PLioqKLCsry9ra2szM7Ic//KGNHDky5HmPPfaYTZgwodvtnj9/3uLi4mz79u3BZadOnbKYmBh74403zCw6x6PXMjRjLLaL1Byv9dZbb3VZrEXjWDTzXo5m0Tce3Zphuz/84Q8WHx9vra2tZsZYvJabczSLvrEYCVx7enlH7adNpKWldbm+paVFr7zyivx+v+6+++7g8kWLFmnKlCmd2s+aNUuDBw/WpEmT9Nprr4Wsa25uVv/+/UOWJSYm6tChQ2ptbZUk7du3Tz6fT9XV1SHttm3bpkAgoDvvvFNPPvmkLly4EFxXXFwsv9+v8ePHB5dNmDBBfr9fBw8e7P1D8AAv5NjuZz/7mQYNGqQxY8bopz/9qVpaWnp9/27V0tKi0tJSTZs2LWT5tGnTdPDgQVVVVamuri5kfUJCggoKCkK+2x1zLC4uVkFBQchcldOnT9fp06eDmRQXF3d63enTp+vvf/97tzmWlpaqtbU15HlZWVnKzc0N9ifaxqMXM2zHWIzcHK9HtI1FyZs5touW8eiFDBsaGpSamqp+/foFt8tYvMrNObaLlrEYKfr13iTymZm+//3v6ytf+Ypyc3ND1v3pT3/S/Pnz1dTUpMzMTO3Zs0eBQCC4PjMzU21tbcHHAwYM0Lp16zRp0iTFxMRo165dmjdvnn7zm9/ooYceknT1C79x40bNmTNHeXl5Ki0t1ebNm9Xa2qr6+nplZmYqKSlJOTk5iouLC257wYIFGj58uDIyMlReXq4VK1boH//4h/bs2SNJqqur0+DBgzu9v8GDB6uurq5PP7NI5JUcJamwsFB5eXkaOHCgDh06pBUrVqiqqkobN2506uMLq/r6en366adKT08PWZ6enq66urrg97er9R9++GHwcccc6+rqlJ2d3ek57euGDx+uurq6Lrd75cqVbnOsq6tTfHy8Bg4c2GV/29tE03j0YoYSY7FdpOZ4PaJtLErezFGKrvHo9gzPnj2rZ555Ro899ljIazMWr3JzjlJ0jcVI4Ymie+nSpXrvvff0zjvvdFo3depUvfvuu6qvr9eGDRv0ne98RyUlJcFfGkVFRSHtA4GAnnjiieDjsWPH6ty5c3ruueeCxdqPf/xj1dXVacKECTIzpaena9GiRXruuecUGxsrSbr33nt17NixkG0vXrw4+O/c3FyNGDFCY8eO1ZEjR5SXlydJ8vl8nd6DmXW53Gu8lOO1r33XXXdp4MCB+va3vx3cq+hVHb+nHb+7va3vmGN3z+m4vLc2XeXYld7621Ubr/FahozFq9yWY2+v3VWfvchrOUbjeHRjho2Njbr//vv15S9/WStXruzxtbvqsxd5LcdoHIvh5vrTy5ctW6Zdu3bprbfe0pAhQzqtT05O1h133KEJEyZo06ZN6tevnzZt2vSZXmPChAn617/+FXycmJiozZs3q6mpSdXV1aqpqVF2drZSUlJCjr72Ji8vT3FxccFtZ2Rk6N///nendh9//HGnPV1e46Ucu3ttSaqsrPxMfXaLQCCg2NjYTnu6z5w5o/T0dGVkZEhSt+u7k5GR0eVzpP/uEe6uTb9+/br9w5GRkaGWlhadO3eu2/5E23j0YoZdYSxGVo7XI9rGouTNHLvi5fHo1gwvXLigGTNmaMCAAfrjH/8YcgSVsfhfbs6xK14ei5HCtUW3mWnp0qXauXOn9u7dq+HDh1/385qbmz/Tax09elSZmZmdlsfFxWnIkCGKjY3V9u3bNXPmTMXEXP9H+v7776u1tTW47YkTJ6qhoUGHDh0KtikpKVFDQ4Puu+++z9Rnt/Bijt29tqQe27hZfHy88vPzQ06xl6Q9e/bovvvuC56Of+36lpYW7d+/v8fv9sSJE3XgwIGQ64x2796trKys4GlZEydO7PS6u3fv1tixY7v9I5Ofn6+4uLiQ59XW1qq8vDzYn2gbj17MsCuMxcjK8XpE21iUvJljV7w8Ht2YYWNjo6ZNm6b4+Hjt2rWr031vGIv/5eYcu+LlsRgxnLxLm5O++93vmt/vt3379oXc7r6pqcnMzD755BNbsWKFFRcXW3V1tZWWltqjjz5qCQkJVl5eHtzO8uXLbeHChcHHW7dutW3btllFRYUdO3bMfv7zn1tcXJytW7cu2Ob48eP26quv2j//+U8rKSmxefPmWVpamlVVVQXblJSUWE5Ojp08edLMzCorK2316tV2+PBhq6qqsj//+c82cuRIu+eeezpNUXTXXXdZcXGxFRcX2+jRoz09FYMXczx48KCtW7fOjh49ah988IHt2LHDsrKybNasWU5+lGHXPqXGpk2brKKiwr73ve9ZcnKyVVdXm9nVKTX8fr/t3LnTysrK7IEHHug0pUbHHM+fP2/p6en2wAMPWFlZme3cudNSU1O7nFLjiSeesIqKCtu0aVOnKTU65mh2dbqpIUOG2JtvvmlHjhyxr33ta11OGRZN49FrGTIW3ZFjbW2tHT161DZs2GCS7MCBA3b06FE7e/ZssE20jUUz7+UYjePRTRk2Njba+PHjbfTo0VZZWRnyf7Jo/rto5r0co3EsRgLXFt2SuvzZsmWLmZldunTJvvnNb1pWVpbFx8dbZmamzZo1yw4dOhSynYcfftgKCgqCj7du3WqjRo2ypKQkS0lJsfz8fHv11VdDnlNRUWFjxoyxxMRES01NtdmzZ3eaKqF9yoz2Aq6mpsYmT55saWlpFh8fb1/84hft8ccfD/lPhZnZ2bNnbcGCBZaSkmIpKSm2YMGCLqfd8Aov5lhaWmrjx483v99v/fv3t5ycHFu5cqVdvHix7z64CPXiiy/abbfdZvHx8ZaXlxcy9VtbW5utXLnSMjIyLCEhwSZPnmxlZWUhz++Yo5nZe++9Z1/96lctISHBMjIybNWqVcHpNNrt27fP7rnnHouPj7fs7Gx76aWXQtZ3zNHs6ndr6dKllpaWZomJiTZz5kyrqakJeV60jUczb2XIWHRHjitXruzx74BZdI5FM2/lGK3j0S0Ztj/u6ufanBmL7s8xWsdiuPnM/v+KfAAAAAAA0Kdce003AAAAAACRjqIbAAAAAACHUHQDAAAAAOAQim4AAAAAABxC0Q0AAAAAgEMougEAAAAAcAhFNwAAAAAADqHoBgAAAADAIRTdAABEoVWrVmnMmDHh7gYAAJ7nMzMLdycAAEDf8fl8Pa5/+OGH9cILL6i5uVmDBg26Sb0CACA6UXQDAOAxdXV1wX/v2LFDP/nJT3T8+PHgssTERPn9/nB0DQCAqMPp5QAAeExGRkbwx+/3y+fzdVrW8fTyRYsWac6cOXr22WeVnp6uL3zhC1q9erWuXLmiH/zgB0pLS9OQIUO0efPmkNc6deqU5s2bp4EDB2rQoEGaPXu2qqurb+4bBgAgglF0AwAASdLevXt1+vRpHThwQOvWrdOqVas0c+ZMDRw4UCUlJVqyZImWLFmijz76SJLU1NSkqVOnasCAATpw4IDeeecdDRgwQDNmzFBLS0uY3w0AAJGBohsAAEiS0tLS9Mtf/lI5OTl65JFHlJOTo6amJj311FMaMWKEVqxYofj4eP31r3+VJG3fvl0xMTHauHGjRo8erVGjRmnLli2qqanRvn37wvtmAACIEP3C3QEAABAZ7rzzTsXE/Hd/fHp6unJzc4OPY2NjNWjQIJ05c0aSVFpaqsrKSqWkpIRs5/Llyzpx4sTN6TQAABGOohsAAEiS4uLiQh77fL4ul7W1tUmS2tralJ+fr23btnXa1i233OJcRwEAcBGKbgAAcEPy8vK0Y8cODR48WKmpqeHuDgAAEYlrugEAwA1ZsGCBAoGAZs+erbfffltVVVXav3+/CgsLdfLkyXB3DwCAiEDRDQAAbkhSUpIOHDigYcOGae7cuRo1apQeeeQRXbp0iSPfAAD8P5+ZWbg7AQAAAACAF3GkGwAAAAAAh1B0AwAAAADgEIpuAAAAAAAcQtENAAAAAIBDKLoBAAAAAHAIRTcAAAAAAA6h6AYAAAAAwCEU3QAAAAAAOISiGwAAAAAAh1B0AwAAAADgEIpuAAAAAAAcQtENAAAAAIBD/g8xRveMFZo8iQAAAABJRU5ErkJggg==", | |
| "text/plain": [ | |
| "<Figure size 1000x500 with 1 Axes>" | |
| ] | |
| }, | |
| "metadata": {}, | |
| "output_type": "display_data" | |
| } | |
| ], | |
| "source": [ | |
| "import numpy as np\n", | |
| "import pandas as pd\n", | |
| "import matplotlib.pyplot as plt\n", | |
| "from scipy.interpolate import UnivariateSpline\n", | |
| "\n", | |
| "# Assume raw_weights is a Pandas Series with datetime index\n", | |
| "timestamps = raw_weights.index\n", | |
| "weights = raw_weights.values\n", | |
| "\n", | |
| "# Step 1: Convert time to seconds relative to the first timestamp\n", | |
| "time_seconds = np.array([(t - timestamps[0]).total_seconds() for t in timestamps])\n", | |
| "\n", | |
| "# Step 2: Compute rate of change (g/s)\n", | |
| "dt = np.diff(time_seconds)\n", | |
| "dw = np.diff(weights)\n", | |
| "rate = dw / dt\n", | |
| "rate_time = (time_seconds[1:] + time_seconds[:-1]) / 2 # midpoint times for rate\n", | |
| "\n", | |
| "# Step 3: Fit a smooth spline to the rate\n", | |
| "spline = UnivariateSpline(rate_time, rate, s=0.001) # s controls smoothing\n", | |
| "\n", | |
| "# Step 4: Plot\n", | |
| "t_fine = np.linspace(rate_time.min(), rate_time.max(), 10)\n", | |
| "rate_spline = spline(t_fine)\n", | |
| "\n", | |
| "# Convert back to datetime for plotting\n", | |
| "t_fine_datetime = pd.to_datetime(t_fine, unit='s', origin=timestamps[0])\n", | |
| "\n", | |
| "plt.figure(figsize=(10, 5))\n", | |
| "plt.plot(timestamps[1:], rate, 'o', label='Raw rate (g/s)', markersize=3)\n", | |
| "plt.plot(t_fine_datetime, rate_spline, '-', label='Interpolated spline', linewidth=2)\n", | |
| "plt.xlabel(\"Time\")\n", | |
| "plt.ylabel(\"Rate of weight change (g/s)\")\n", | |
| "plt.title(\"Interpolated Rate of Weight Change\")\n", | |
| "plt.legend()\n", | |
| "plt.grid(True)\n", | |
| "plt.tight_layout()\n", | |
| "plt.show()" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": null, | |
| "id": "f0529fcd", | |
| "metadata": {}, | |
| "outputs": [], | |
| "source": [] | |
| } | |
| ], | |
| "metadata": { | |
| "kernelspec": { | |
| "display_name": "Python 3 (ipykernel)", | |
| "language": "python", | |
| "name": "python3" | |
| }, | |
| "language_info": { | |
| "codemirror_mode": { | |
| "name": "ipython", | |
| "version": 3 | |
| }, | |
| "file_extension": ".py", | |
| "mimetype": "text/x-python", | |
| "name": "python", | |
| "nbconvert_exporter": "python", | |
| "pygments_lexer": "ipython3", | |
| "version": "3.13.5" | |
| } | |
| }, | |
| "nbformat": 4, | |
| "nbformat_minor": 5 | |
| } |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment