{ "cells": [ { "cell_type": "markdown", "id": "d146d425", "metadata": {}, "source": [ "# Estimating the Start and Stop of Motion\n", "In this notebook we show how to detect the start and stop of the probe using the \n", "accelerometer data. \n" ] }, { "cell_type": "code", "execution_count": 11, "id": "99a7023d", "metadata": {}, "outputs": [], "source": [ "# Import the function to read the data \n", "from study_lyte.io import read_csv \n", "\n", "# Import the function to the events\n", "from study_lyte.detect import get_acceleration_start, get_acceleration_stop\n", "\n", "# Impor plotting functions\n", "import matplotlib.pyplot as plt \n", "\n", "# Open the file \n", "df, meta = read_csv(\"./data/acceleration_example.csv\")\n", "\n", "# Rolling mean to smooth out the signal\n", "df = df.rolling(window=100).mean()" ] }, { "cell_type": "markdown", "id": "9096a529", "metadata": {}, "source": [ "## Raw Data\n", "Below is the raw NIR profile. In it you can see lag time before the probe enters the snow this can make interpretation difficult.\n", "A good starting place is to trim the data only to when the probe is moving. " ] }, { "cell_type": "code", "execution_count": 12, "id": "97f277f7", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAJOCAYAAAC6OADbAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABDtUlEQVR4nO3deZxcVZ3//9fpTu/7vqRXQhKSdMJiwLBHkFVWkRFGDco4OCrqiKMizgN/8HNckK/jODoyOPJ1R0AJoKIERyEsCZBA9o3urL3v+959vn90VVG9d6CrT3Xu+/l43EffuvdW1afuI6l33XPOvddYaxEREe+KcF2AiIi4pSAQEfE4BYGIiMcpCEREPE5BICLicQoCERGPUxCIhIgx5kPGmA2u6xCZjoJAwpYx5rAxpscY02mMqTXG/NQYkzhH7/3/GWOsMebdM9y+xLf9Av8ya+2vrLWXhqC2n/re66ygZScbY2zQ4+eMMR/3za81xgz79mOHMWa/MeZjs12XzF8KAgl3V1trE4HTgNOBr4T6DY0xBlgHNPv+hqNm4OvHsX21bz8mA58HfmyMWRqSymTeURDIvGCtrQWeYSQQADDG3GmMqfD9yt1jjLk+aN0RY8y7fPMf8v2CXuF7/A/GmCemeLvzgTzgs8BNxpjooNeNM8b8H9/rtxljXjTGxAEbfZu0+n55n22M+agx5kXf835kjLk/+E2MMU8aY+7wzecbY35njGkwxhwyxnx2ml3yM2CVMebCabYbxY54mpEgWXU8z5UTl4JA5gVjTAFwBVAetLiCkS/tFOAe4JfGmDzfuueBtb75C4GDwAVBj5+f4u1uAX4PPOp7fHXQuvuBdwHnAOnAl4DhoNdOtdYmWms3jXnNh4EP+o42MMakAZcCvzHGRPjebzuwELgY+GdjzGVT1NgNfAP4tym2GccYE2GMuQbIZPS+FA9TEEi4e8IY0wEcA+qBr/lXWGsfs9ZWW2uHrbWPAG8C/nbz5xn5woeRsPhm0ONJg8AYEw/cCPzaWjsA/BZf85DvC/tW4HPW2ipr7ZC19mVrbd8MPscLgPXVAvABYJO1tho4E8iy1t5rre231h4EfgzcNM1r/jdQZIy5Ygbvn2+MaQV6gPXAHdbaN2bwPPEABYGEu+ustUmM/Lo/hZFfsgAYY9YZY7YZY1p9X3JlQeufB873HSFEMvLr/lxjTAkjRxDbJnm/64FB4Gnf418BVxhjsnyvHcvIkchxsSNXd/wNcLNv0d/7XhugGN8XddBnuQvImeY1+4D/3zdNp9pam8pIH8H3gYuO9zPIiUtBIPOCtfZ54KeMNM1gjClm5Ffz7UCG70tuF2B825cz0nzyGWCjtbYdqAVuA1601g5P8la3AInAUWNMLfAYEMXIF3cj0AssmqjEGXyMh4EP+Gp/N/A73/JjwCFrbWrQlGStvXIGr/l/gVTg/TPY1h8eXwZWGmOum8lz5MSnIJD55HvAJcaYU4EERr58GwB8wyHLxmz/PCNB4W8Gem7M41GMMf72+asY6ZQ+DTgV+DawzhceDwHf9XXuRvo6hWN8dQwDJ01WvK8pphH4H+AZa22rb9WrQIcx5su+zuhIY0yZMebM6XaItXaQkeayL0+3bdBz+oH/A9w90+fIiU1BIPOGtbYB+Dlwt7V2DyNfZpuAOmAl8NKYpzwPJPHWiJ6xj8f6CLDNWrvBWlvrnxhpSllljCkD/gXYCbzGyMibbwMR1tpuRjpuX/I176yZ5D1+DbzX99f/uYZ4K3wO8VZYpEy7U0Y8DNTMcFu/hxjpX7h62i3lhGd0YxoREW/TEYGIiMcpCEREPE5BICLicQoCERGPWzD9JuElMzPTlpSUuC5DRGRe2bp1a6O1NmuidfMuCEpKStiyZYvrMkRE5hVjzJHJ1qlpSETE4xQEIiIepyAQEfE4BYGIiMcpCEREPE5BICLicQoCERGPUxCIiHicgkBExOMUBCIiHqcgEBHxOAWBiIjHKQhERDxOQSAi4nEKAhERj1MQiIh4nIJARMTjFAQiIh6nIBAR8biQBoEx5nJjzH5jTLkx5s4J1scYYx7xrX/FGFMSynpERGS8kAWBMSYS+CFwBbAcuNkYs3zMZv8AtFhrTwb+Hfh2qOoREZGJhfKI4Cyg3Fp70FrbD/wGuHbMNtcCP/PN/xa42BhjQliTiIiMEcogWAgcC3pc6Vs24TbW2kGgDcgY+0LGmNuMMVuMMVsaGhpCVK7Mpdtvv53rrrvOdRkiwjzpLLbWPmitXW2tXZ2VleW6HJkF1dXVVFRUuC5DRAhtEFQBhUGPC3zLJtzGGLMASAGaQliThInIyEiGhoZclyEihDYIXgMWG2NKjTHRwE3AU2O2eQq4xTf/AeCv1lobwpokTCgIRMLHglC9sLV20BhzO/AMEAk8ZK3dbYy5F9hirX0K+AnwC2NMOdDMSFiIB0RGRjI8POy6DBEhhEEAYK19Gnh6zLK7g+Z7gRtDWYOEp4iICB0RiISJedFZLCceNQ2JhA8FgTihIBAJHwoCcUJBIBI+FATihIJAJHwoCMQZjRQWCQ8KAnGip6eHuLg412WICAoCcaS/v5+YmBjXZYgICgJxxBijpiGRMKEgECcUBCLhQ0EgTigIRMKHgkCcUBCIhA8FgTihG9GJhA8FgTijIwKR8KAgECfUNCQSPhQE4oSCQCR8KAjECQWBSPhQEIgTCgKR8KEgECcUBCLhQ0EgTigIRMKHgkCcUBCIhA8FgTihIBAJHwoCcUJBIBI+FATihIJAJHwoCMQJBYFI+FAQiBMKApHwoSAQJxQEIuFDQSBOKAhEwoeCQJxQEIiEDwWBOKEgEAkfCgJxQkEgEj4UBOKEgkAkfCgIxAnds1gkfCgIxBkdEYiEBwWBOKGmIZHwoSAQJxQEIuFDQSBOGGMYHh52XYaIoCAQR9ra2khNTXVdhoigIBBHuru7iY+Pd12GiKAgEEestRpCKhImFATihIJAJHwoCMQZBYFIeFAQiBMaOioSPhQE4oSahkTCh4JAnFEQiIQHBYE4oaYhkfChIBAn1DQkEj4UBOKMgkAkPCgIxAk1DYmEDwWBOKGmIZHwoSAQZxQEIuFBQSBOqGlIJHwoCMQJNQ2JhA8FgTihIBAJHwoCERGPUxCIEzoiEAkfCgJxQkEgEj4UBOKMgkAkPCgIxAkNHxUJHwoCcUJNQyLhQ0EgzigIRMKDgkCcUNOQSPhQEIgTahoSCR8KAnFGQSASHhQE4oSahkTCh4JAnFDTkEj4UBCIMwoCkfCgIBAn1DQkEj4UBOKEmoZEwoeCQJxQEIiEDwWBiIjHKQjECR0RiIQPBYE4oSAQCR8KAnFGQSASHhQE4oSGj4qEDwWBOKGmIZHwoSAQZxQEIuFBQSBOqGlIJHwoCMSJ2NhYenp6XJchIigIxJH4+Hi6u7tdlyEiKAjEka6uLhITE12XISIoCMSRvr4+YmNjXZchIigIxJH+/n6ioqJclyEiKAjEkYGBAaKjo12XISIoCMSRgYEBHRGIhAkFgTihpiGR8KEgECeGhoaIjIx0XYaIoCAQRyIjIxkcHHRdhoigIBBHmpqaSE9Pd12GiKAgEEeam5vJzMx0XYaIoCAQB4aHhzV8VCSMKAhkzrW1tQGQlpbmuBIRAQWBOFBZWQlAbm6u40pEBBQE4kBVVRUARUVFjisREVAQiAMNDQ0AZGVlOa5EREBBIA60tLQAkJqa6rYQEQEUBOJAZWUlMTExGj4qEiYUBDLnjh07RkFBgW5eLxImFAQy544ePUphYaHrMkTER0Egc6qvr49t27axcuVK16WIiI+CQObUK6+8Qnd3NxdffLHrUkTER0Egc+q5557DGMMFF1zguhQR8VEQyJx68cUXWbVqlS4vIRJGFAQyZ6y1bN26ldWrV7suRUSCKAhkzmzfvp3m5mbOPvts16WISBAFgcyZJ554goiICK6++mrXpYhIEAWBzJn169dz7rnnkp2d7boUEQmiIJA5cfDgQXbs2MF1113nuhQRGUNBIHNi/fr1AFx//fWOKxGRsRQEMicee+wxTj31VEpLS12XIiJjhDQIjDGXG2P2G2PKjTF3TrD+o8aYBmPMNt/08VDWI24cPHiQV155hZtvvtl1KSIygQWhemFjTCTwQ+ASoBJ4zRjzlLV2z5hNH7HW3h6qOsS9733ve0RGRioIRMJUKI8IzgLKrbUHrbX9wG+Aa0P4fhKGKioqeOCBB7j11lt1a0qRMBXKIFgIHAt6XOlbNtYNxpgdxpjfGmMmvDaxMeY2Y8wWY8wW/20OJfxZa/mnf/onYmNj+drXvua6HBGZhOvO4t8DJdbaVcCzwM8m2sha+6C1drW1drXuczt//PKXv+Qvf/kL3/zmN1m4cKLfACISDkIZBFVA8C/8At+yAGttk7W2z/fwf4B3hbAemUO7du3is5/9LGeffTaf/OQnXZcjIlMIZRC8Biw2xpQaY6KBm4CngjcwxuQFPbwG2BvCemSOvPrqq7z3ve8lLi6OX//610REuD7wFJGphOx/qLV2ELgdeIaRL/hHrbW7jTH3GmOu8W32WWPMbmPMduCzwEdDVY+E3uDgIPfddx8XXHABcXFx/PWvf6WkpMR1WSIyDWOtdV3DcVm9erXdsmWL6zIkSH9/P08++ST33HMPu3fv5vrrr+e///u/UX+OSPgwxmy11k54DfiQnUcgJ7bm5mY2btzIhg0bePTRR2lqamLp0qX87ne/4/rrr8cY47pEEZkhBYHMSF9fHxs3buTPf/4z//u//8uOHTuw1hIXF8c111zDunXruPTSS1mwQP+kROYb/a+VKQ0PD/PhD3+Yhx9+GIDIyEguvPBC7r33Xt7znvdw5plnEh0d7bhKEXknFAQypcHBQYL7ZIaGhqioqGB4eJjy8nL+/Oc/U1RUNGpKSEhwWLGIHC8FgUwpOjqaXbt2sXXrVvbs2cOBAwfYvHkzL730Es8999yEz3n88cd1uWmReURBIFOy1vKTn/yE9evX88Ybb9DY2Djl9meccQZlZWVzVJ2IzAYFgUyop6eHgwcP8vzzz/PpT38agDPPPJPPfOYz5OXljZqys7PVSSwyj+l/rwAjzTnr16+nsrKS8vJyKisrx23z4x//mFNPPdVBdSISSgoCAeCee+5hx44dAKxbt46TTz45MC1atIj09HTHFYpIqCgIPGx4eDhwhdCjR48CsGDBAv7rv/5LI39EPERB4GGHDh3illtuAeCaa67hzDPPZO3atQoBEY/RZSE9rK+vLzC/aNEikpKSaGhoYOvWrRw5coTOzk7m27WoROT46YjAwxYvXswdd9zBE088wQMPPEBPT8+4baKjo0lPTyc9PZ2MjAyysrLIysoiMzMzMO+fiouLSUtLc/BJROSd0NVHBRg5X6CxsZGKigoOHjxITU0NtbW11NbWUldXR21tLTU1NdOeR+D3wgsvcN5554W4ahGZKV19NMSstQwODjIwMDBu6u/vn3D5XG071Xb+df6/g4ODs7ZPhoaGZu21RCS0PBMEfX193H///Tz88MPExcWRnJxMYmIi0dHRDA0NvaMv4dn8Ap1KREQEUVFRgSk6OnrU44km/2ccu9y/7O38jY6OJj4+ftQUFxdHamoqMTExc7IvRGT2eCYI3njjDf71X/911l83IiKClJQUkpOTSUpKIikpKTA/dtnY5TExMRN+eU/2Ba9bPopIKHgmCNasWcOBAweoqamhubmZ5uZmWlpa6O7uHjV1dXXR19dHT08Pvb29gb8Tzff39zM8PExbWxttbW3HXVNUVBSxsbHExcURGxs7an6iZcezbXJycqBTV8NBRWQqngkCGBkls3jx4ll7veHh4VGhMV1wHO98c3PzpOuPpzkqLi6OjIwMoqKi6OrqoqOjg9tvv51vfetbOsoQEY0amq96e3upqamhqqpq3FRZWRkY+TOVuro6srOz56hiEXFJo4bCmLWWjo4OmpqaRk3Nzc1TPp6qKSoyMpKMjAxWrFgRGOOfm5tLSUkJxcXFlJSUsGTJEpKSkubwk4pIuFIQOFJRUcGyZcsYGBiYcruUlBQyMjICJ3QtXrx41GP/FPw4OTlZN48XkRlTEDgyPDw8bQhERkYSHR2NMYahoSF6e3vp6OgAYGBggJ6eHtrb22lsbAyMSvIHR2ZmJhkZGcTFxc3FxxGReUxB4MjixYux1tLe3k5VVRXV1dW0tbXR3t5OR0cH7e3toyb/ssrKylHLg68XNJH4+HgyMzMDU3BIJCUlkZCQMGqKj48ftywxMVE3nhE5gel/t2P+X/LLli17W8/v6+sLhERHRwetra00NTXR2NhIY2NjYL6pqYmGhgYqKipobGw87uGuMTExJCUlkZiYGJgmeuwftpqTk0N2dnZgSkxMVHOVSJhSEMxzMTExxMTEkJmZeVzPGxgYoKurKzD5z6EYu6yzs5Ouri46Ozvp7Oyko6MjMN/Z2UlNTc2ox729vRO+X2xs7KhwSE9PJyUlZcrJf4TiP3NZQ11FQkNB4FFRUVGkpqaSmpo6q6/b19dHQ0MD9fX11NfXU1dXF5j3T9XV1ezatSvQFDY8PDyj1/afMDf28hb+oJjJ8pNOOok1a9bM6mcWme8UBDKrYmJiKCgooKCgYEbbW2vp7OwMnJ0dPI0967unp2fcMv/y5uZmdu/eTX9//7Tv+dRTT3H11Ve/048qcsJQEIhTxpjA9ZhmGh5TvdZUsrKy+NCHPsSFF174jt5H5ESjM4vlhLFhwwYuu+wyzj77bJqamqisrKS7u3vcdikpKeTl5ZGXlxe4Cm3wCKmxo6ammmJjY9UJLvOCziwWT7j00ktH3VrTWktra2vgshs1NTWjptraWg4dOhToDPd3kh/Pj6OIiIgpg8J/tOOf/KOrJpr86xQuMtcUBHLCMsaQlpZGWloaZWVlM3qOtZaenp5xI6j8U3BgTDfV1dXR0dERGGk13TkffpGRkdOGx2SBMtHl0KOiot7JbhQPUBCIBDHGBEYYZWVlzepr9/f3B4bgBk8TLZtofX19/dsKltjYWDIyMsad2zHRlJWVRWxs7Kx+bgl/CgKRORIdHU16ejrp6emz8nqTBUvwmegdHR20tbXR1NQUGL67d+9e6urqJj3nIycnh+uuu45zzjmHG2+8UZcp8QB1Fot4kLWWrq6uced41NfXs3XrVh5//HEAbrjhBr71rW+Rk5Ojs8Pnuak6ixUEIjLKwMAAS5cu5dChQ6OWx8XFkZOTQ05ODmlpaYFbtPrPBA+en2id7mftlkYNiciMRUVFUVFRwZEjR9i3bx91dXWByX+2eFNTEwcPHgycHd7T0zPt6yYkJLBo0SL279/PQw89xN///d/PwaeRmVAQiMg4xhhKSkooKSmZ0fb9/f2BK+L6zwwfO19XV8f3v/99AD70oQ+xZMkSli5dqhskhQEFgYi8Y9HR0YFLnU9l586d/O1vfwPgzDPPBKCwsJBly5axfPlyVq5cyXXXXTdrHeoyM+ojEJE51dXVRXl5OQcPHmTv3r3s3buXPXv2sHfv3kATU3FxMbfddhtlZWWUlZVRWlqqjup3SJ3FIhJ2rLX09vYGmo9aWlp46qmn+OY3vzlu27i4uAkvFyIzp85iEZl1wV/iwVNra+uEyyeaZnK1WIAf/vCHIf403qYgEJFJNTc38+1vf5sNGzaQmJgYuNvdTL/Ek5KSRg0pzc7OZvHixYHHqamp425I5L9PRmpqqs5dmCMKAhEPs9YG7lY39qzk6upqbr311lHbf+ADHyAzM3Pau8ulpKSQlJREZGSko08mx0NBIBLment7R11zaOwtQ6e7ac90y4eGhqatIScnh//8z//kxhtvnINPLHNNQSASBn75y1/ykY98JPA4OTmZqKgoOjo6ZtyODiMng012u87s7Owpb+M59uqlWVlZFBcX64xgD1AQiISBffv2jXrc3t4+bpu8vDyKi4spLi6mqKiI4uJiCgsLA1NqaioLFui/tBw/DR8VCRNDQ0McPHiQffv2UVFRwcGDBzl69ChHjhyhurqa5uZmBgcHJ31+fHx84OqmK1eu5Nxzz6WoqCgwpaSkzOGnkXCj8whETgDWWjo7O2lubh43NTU1BeYbGxt5+eWXaWpqGvX85ORkioqKKCwsJC0tbdwNbiaaD16WkJBARESEo08v75TOIxA5ARhjAl/KxcXFU247PDxMbW0tR48eHTcdO3aMAwcOBO5fMJMLxvn57+vsD4ngaapl/iuRJicnB+YVLOFDQSByAoqIiCA/P5/8/HzWrFkz5bZDQ0PjRiRNNe8fauoftdTU1MThw4cDjzs6OmY0EskfbGMDwt9hPdW9oOPj4yddpyGrx09BIOJxkZGRgbH/s8FaS39//6gA8V+Z1H9F0snmm5qaOHToEO3t7YF7P88kVILFxMQEgiL4/s3+KfjxRPOpqank5uYSHx8/K/tjPlAQiMisMsYQExNDTEzMtFcjnY4/VLq7uwPB4J9msqyzszMQNLW1taMCaXh4eMr3TkxMJDc3l5ycHBYuXMjKlSu5/PLLOf3000+4ow51FouI51hr6e7uHnUmtX9qbm4edTOeQ4cOsWXLlnFHJt///vf5zGc+4+gTHD91FouI5wwMDNDS0kJLSwvNzc20traO6/OYqB/E/7e1tZWmpiZ6e3snfP1nn312XgXBVBQEIjJvbd++nY9//ONERkYSExMTuJx1c3MznZ2d0z4/JiZm3LDZ5ORk8vPzSU1NJSMjg7S0NLKzs8nOziYrKyswn5iYOAefcG4oCERk3vrOd77DdE3FGRkZgfMn/COp8vLyyM3NJT09naSkJHJzc8nMzPTslU7VRyAi81Z/fz8vv/wy+/fvD5wn0dTUNOp+yf6mnqnOyoaRo4OCggKio6Pp7e3lgx/8IHfffTdxcXFz9GlCS2cWi4inWWvp6+ubsE+gvb2dmpoaKisr+c53vjPqeZdddhlXXXUVmZmZZGRkkJGREZiPj4+fV0cQCgIRkRmoqqriyiuvZMeOHcDIUNjJviNjY2PJzMwkMzOTrKyswN9zzz2XU089lZNPPjmshpkqCERE3obBwcHAtZz8d2drbGwcNd/Q0EB5eTn79+8f9/zXX3+d008/3UHl42n4qIjI27BgwQLS09Pp7e2lqamJrq4uGhoaAleF9fdLdHR0jHvuokWLpr0mVLhQEIiIBDl27BhFRUUAFBYWUlVVNe4s5IyMDIqLi1m8eDEXX3xx4B4R/vtEZGVlzasL6ikIRER86urqAiEAI6GwePFi/uVf/mXUl31CQoLDKmff/IksEZEQy8zM5LOf/eyoZW+++SannXYal112GcuWLTvhQgAUBCIiAZGRkfzHf/wHTz755KjlX/ziF7nkkku49957J73kxHympiERkSCtra288sornHLKKYF7SW/cuBGAv/zlL1x11VWcccYZLkucdQoCEZEgzz33HN/4xjcmXHfuueeyceNGysvLOe2001iyZMkcVxcaOo9ARCSItZbnn3+el19+maNHj1JbW0tdXR21tbXU1NTQ19c3avsvfvGLrF69mpUrV5KZmUlaWhoLFoTfb2ydUCYiMgustbS0tLB7924+//nPs3fvXrq7u8dtl5qaSnp6euCyFP759PR0UlNTx03p6emkp6eTkJAQsstWKAhEREKko6ODnTt3cvjwYZqamgJnIk8039raOuVrRUVFBULBHxqRkZH09fWxatUqvvKVr5CWlva26lQQiIg41N3dTUNDA3V1dVRUVFBeXs7hw4c5dOgQBw8e5MiRIzN6nR/84Ad8+tOffls16BITIiIO1NfXk5OTM+U2iYmJFBYWBpqJ0tLSxv3NyckhPz+fc889NyR1KghEREIkLi6Os846i1dffXXcuoSEBAoKCvjXf/1XPvzhDzuo7i1qGhIRCaHe3l6ee+45nn76aTZu3Mj27dvHbTMX38NqGhIRceSrX/0q3/3ud0ctS09P5+KLL6akpITPfe5zjip7i44IRERCqLKykq9+9atUVFTw0ksvTbjNG2+8wWmnnRbSOnREICISAtZaampqOHLkCJs2beLo0aMANDQ00NDQQH19PQ0NDTQ2NtLf3z/p67S3t89VyRNSEIiIvA07d+5k1apV45YnJycHbl1ZWFjIGWecQVZW1qjJPwooKysrLM5Cdl+BiMg8VFBQwN/93d/x6KOPBpYtX76c3bt3O6zq7VEQiIi8DWlpaTzyyCOUlpby7W9/G4A9e/awZMkScnJyyMnJYeHChRQWFlJYWEhBQQGlpaXk5eWF7DISb5c6i0VE3qG2tjaeeeYZNm/eTHV1deAidZWVlXR2do7aNjU1leXLl5OWlsadd97JeeedNyc16hITIiKOtLW1UVlZybFjx9ixYwcPP/ww27ZtC6z/9a9/zc033xzyOjRqSETEkZSUFFJSUjhw4ABf/vKXR63LzMzk/PPPd1TZW3SrShGROVBUVDTq8dVXX82xY8coKChwVNFbdEQgIhIC7e3tfO1rX+Pw4cM0NDQEbnvp9/vf/56mpiYWLlzoqMK3KAhERGZRV1cXf/rTn7jxxhsDyy688EJuuOEGli5dyrJly1i6dClFRUVhcQ4BKAhERGbFd7/7Xb7whS+MW37gwAEWL17soKKZUxCIiMyCJ554YtTj7OxsioqKuOOOO8jLyyM3N5e8vDwKCgooKiqirKyMyMhIN8WOoSAQEZkFGzdu5NVXX2Xz5s2BG937zyXYsmUL9fX1DA8Pj3pOb28vMTExjip+i4JARGSWnHXWWZx11lnjlg8NDXH48GE2bdrERz7ykcDy2NhY6urqyM7Onssyx1EQiIgch+HhYZqamqiurqayspLW1tYpp+bmZo4dO8bAwMC417r66qtJTU2d+w8xhoJARCSItZa6ujp+8pOfsH//fowxgctJ+y8dMTg4OOFz4+PjA/ceTk1NJTs7myVLlnDjjTeyaNGiwFRQUBA2/QOgIBARD2ppaeHQoUOjpsOHDwf+9vT0BLYtLi4mMzOTrKwsysrKyMvLC3T+Lly4kOzsbFJTU0lJSSE6Otrhp3r7FAQickJ79dVX+cUvfkFDQwNvvvkmFRUVtLW1jdomJSWF0tJSTjnlFK644gpKS0s56aSTOP/880lKSnJU+dxREIjICesXv/gF69atA6CkpISlS5eyZs0aSktLR01paWmOK3VLQSAiJ5TnnnuOrVu3sm/fPv7nf/4nsPzQoUMOqwpvCgIROWHs3r2b97znPeOWf+xjH+Ob3/wmaWlppKamBv76p4SEBOLj44mI8OZ1OBUEInLCWL58OV/60pfYtm0bHR0dtLS00NLSwq9+9aspbx7vFxcXR3x8PAkJCYFwCP47dj4hIYH09PTAvYizs7MpLS0NqxFBM6Eb04iIJ/T29tLS0hIY3x8839XVRVdXF93d3YH5sY/Hznd3d0/6XkVFRWzfvj0szhHw041pRMTzYmNjA0M/Z8Pw8DDd3d00NTXx5ptvcskllwTWHT16lAceeIA777xzVt4r1BQEIiJvQ0REBImJiSQmJvKnP/0psPy2227joosu4rrrrnNX3HFSEIiIzEBfXx9vvvkmx44d4+jRoxw7dowjR47w+uuvj7rpzM0338zatWvdFfo2KAhERGbgmmuuYcOGDZOuP++887j22mtJSkqirq6OrKyseTMKSUEgIjIDd999N/v27aOnp4fBwUF6enro7e0NrH/xxRd58cUXRz2nrKyMTZs2kZiYONflHhcFgYicUKy19Pb20tnZGRjlM1vzwV/8M1FQUDAvrj80ZRAYY3bM4DUarLUXz1I9IiJTOnDgAPfddx9DQ0PjLvnc1tZGR0fHpFcHnUhkZCSJiYkkJCQE/iYkJJCRkUFRUdGo5cnJyYEpJSVlwsdxcXEYY0K4B2bfdEcEkcCVU6w3wFOzV46IyIjy8nK+/vWv09bWRldXF3V1dYHJb+XKlaSmplJYWEhZWVngyzgpKSkwosf/xT7ZfHR09Lz74p5t0wXBJ6y1R6bawBjzqUmWPwRcBdRba8smWG+A/2AkaLqBj1prX59R1SJywnvsscf42c9+NmpZTk4Oa9euZdmyZdx3331h3/Y+Xxz3mcXGmDSg0Fo7ZbORMeYCoBP4+SRBcCXwGUaC4N3Af1hr3z3d++vMYhFvsNby/ve/f9xN4YNlZWVRXFw8aioqKqK4uJjS0tKwOrPXtXd8ZrEx5jngGt/2W4F6Y8xL1to7JnuOtXajMaZkipe9lpGQsMBmY0yqMSbPWlszk5pE5MRmjGH9+vUA9Pf3U1NTQ319fWCqqqri6NGjHDlyhF27dvH000+PuqGMX2ZmJs8//zzLly+f648wb8x01FCKtbbdGPNxRr68vzbDjuSpLASOBT2u9C0bFwTGmNuA22DkGh4i4i3R0dGBX/zBXnrpJc477zze+973kpyczO7du8ddA6ixsZHt27crCKYw0yBYYIzJA/4O+GoI65mQtfZB4EEYaRqa6/cXkfD00Y9+FIC//OUv49YtWrSIU089lbVr15KVlcX27dtZuHAhmZmZc1xl+JtpENwLPAO8aK19zRhzEvDmO3zvKqAw6HGBb5mIyIzs3r2bX//61xw4cIC9e/fS0NBAQ0MDjY2NVFRUUFFRweOPPz7qOZmZmSxbtoxly5ZRUlLCmjVrJryHgZeE9DLUvj6CP0zSWfw+4Hbe6iz+vrX2rOleU53FIjITg4ODNDc309jYSGNjIw0NDRw5coS9e/eyd+9eNm/ezNDQUGD7u+++m3vuucdhxaH1tjuLjTG3+ZpljnsbY8zDwFog0xhTCXwNiAKw1j4APM1ICJQzMnz0Y9N/FBGRmVmwYAHZ2dlkZ2cDMDAwwP333z/q9pXBgs9P8JrpmobuNMY0TrHeAJ/D134fzFp781Qv7Bst9OlpKxQReZsGBgaora2lpqaGG2+8kaNHjwKQm5vLF77wBYqLiznppJNYuXLlvLgURKhMFwTPA1dPs82zs1SLiMisWL9+Pe9///snXV9To1HqwaYMAmutmmtEZF7p7u7m4YcfDjyOjY3lscceIy8vj/z8fHJychxWF5509VERmbc+8YlP8NOf/pRTTz2VhoYG6uvrx51H8MEPfpCrrrrKUYXzg4JAROaVzs5O9u/fz+rVbw2AaWlp4bzzziMrK4usrCyys7MpKChg0aJFnHTSSQ6rnR8UBCIS9tra2vjtb3/LE088wR//+EeCh73/9a9/9fx5AO/UTK81lAN8A8i31l5hjFkOnG2t/UlIqxMRAb70pS/x4IOjByeuWrWKiy66iK1bt9LX18eyZcvIzs4mLi7OUZXz10yPCH4K/F/eurzEAeARQEEgIiF3zz33kJaWRn9/P42NjVRWVlJVVcWPf/xjurq6Rm2bkJBATk4OeXl55ObmkpeXR1FREUuWLGHJkiWcdNJJxMTEOPok4WlGZxYbY16z1p5pjHnDWnu6b9k2a+1poS5wLJ1ZLCJ+1lpaW1vZtm0bFRUVgUtM1NbWBqaamhpaW1vHPfePf/wjV1451X23Tizv+DLUQJcxJgOwvhdcA7TNUn0iIsdteHiYmpoajh07RmtrKwMDA/T39zMwMACMnFkcGxtLbGzshM8fO7rIy2YaBHcwckvKRcaYl4As4AMhq0pEJIi1lubmZvbu3ctrr73Giy++OO5icn5paWlkZGSQmZlJfn4+Z5xxxrib1yxcuJCoqKg5/hTha0ZBYK193RhzIbCUkctK7LfWDoS0MhHxrPb2du666y7+8Ic/MDg4SEdHB+3t7YH1paWlvP/972fZsmWcffbZ5Ofnk5ubS1ZWFgsWaDDk8ZrpqCH/TexLfM+51BiDtfa7IaxNRDxq27Zt/PCHPxy3/Mknn6SsrEznBsyymUbn74FeYCcwHLpyRMQrtm7dyuHDh+ns7KSxsZGmpqbAJaMnuhbQ5ZdfzjXXXOOg0hPfTIOgwFq7KqSViIhnVFdXjzoz2C83N5eMjAxycnL4h3/4B4qKioiPj+fzn/88kZGRDir1hpkGwZ+MMZdaazeEtBoROeFZa0lMTOTf/u3fuP/++2lpaQms8w/5jI+P50c/+pE6dOfITINgM7DeGBMBDDDSYWyttckhq0xEwlpXVxfV1dW0trbS0tJCS0vLqPnJHre2tjI8PHULc2Nj47TbyOyZaRB8Fzgb2GlDeW9LEQl7u3fvpqxs3N1nR4mOjiYtLS0w5eTksHTp0sDj1NRUUlJSSEpKIjk5meTkZJKSkgIjf4wxc/RpBGYeBMeAXQoBEZmorf7FF18MfMGnpaXpej/zzEyD4CDwnDHmT0Cff6GGj4qcuLq6ugKXaBg7XXrppWzY8FaX4bnnnuuwUnmnZhoEh3xTtG8SkTDnv0BbU1MT1dXVNDY20t7eTltbG+3t7VNObW1t9Pb2jnvNBQsWBC7kdtVVV1FSUsJdd93l4NPJbJrpmcX3hLoQETl+fX19tLW1sXnzZvbs2UNVVRVVVVWsX79+yudFRkaSkpISaJ9PTk4mJyeHJUuWBB6npaWRl5cXuIpnfn4+GRkZREREzNGnk7kyZRAYY35grb3dGPN7fBecC2at1dkdIpPo7++no6ODzs5Oenp66Onpobe3d8q//qm7u3vUFLzM/4u9ra2Nvr6+Ue+ZmprKwoULgZFr7nzjG98IjMvPzc0lOTmZlJQUYmNj1SErAdMdEawDbgfun4NaRMKCtZampiaOHj0aaCrp6OiYcJpsXXt7O/39/W/r/aOjo4mPjyc+Pp64uLhR89nZ2Zx88smkpKSMmpYuXco555xDQkLCLO8N8YLpgqACwFr7/BzUIhJSQ0NDtLS0jLqcQV1dXaBDtLa2lmPHjlFeXk5b29RXWU9KSho3lZSUBIZDBi9PTEwMfJHHxsYSFxc36XxsbKzOoJU5N10QZBlj7phspUYNyVyy1tLZ2Rk4Kck/tbW1jVs2dvKf0DTZCOjMzExyc3NZuHAha9asYdGiRZSUlJCamjruCz8hIUHt5HJCmS4IIoFERs4kFnnbrLX09vYGmk2Cp+P5Up/ubNP4+HhSU1MDk/9EptTUVDIzMwPXqc/IyCAjI4Pc3Fyys7OJjtZgOPGu6YKgxlp775xUIvPG0NAQR48e5ciRI6OaVVpaWgJt5hMNSRwcHJz2tRMSEkZ9kefl5bFs2bJRyyaa/G3l+kIXOX7TBYGOBARrLU8//TSPPPIIr7/+OuXl5eNGq0RHR5Oenj7qcgGlpaWjhicGX04geFnwl7kuMiYy96YLgovnpAoJa48++ig33XQTAGeeeSa33347p5xyCqWlpYFx5qmpqRqOKDJPTRkE1trmuSpEwlfwzb9fe+01tm3bRl5eHgsXLiQ/P5+FCxcG2tqzs7PJyckJzMfHxzusXERmQjf3lGlde+21dHR08Prrr7Nz587A2atVVVXs3r2bDRs20NHRMeFzExMTycnJGTfl5eVRWFgYmFJSUnREIeKImW8XFF29erXdsmWL6zJkjJ6eHhoaGqirq6O+vp76+nrq6urGTbW1tTQ1NY17fmJiIgUFBRQVFbFixQrOOOMMTj/9dJYuXaqbkYvMAmPMVmvt+NvCoSMCmSVxcXEUFRVRVFQ07bYDAwPU1dVx7NixcdORI0f40Y9+FLjgWVxcHO9+97s5/fTTKSoqIj8/P9AvkZmZqSMJkVmgIwIJO4ODg+zfv5833niDrVu38sILL7Bnzx56enrGbRsZGRk4J2DseQJj/2ZlZZGbm0tiYqLCQzxnqiMCBYHMC9Zampubqampobq6mpqamsBlIoL/Bs8PDAxM+Frx8fHk5uYGpsLCQpYvX84555wz7Z23ROYrNQ3JvGeMCfzyn8mXtbWWjo6OUcHg77fw3yC9traWffv28ec//5nu7m4WLFjAE088wfve9745+EQi4UNBICckY0zgpLXS0tIptx0eHubll1/mwx/+MFdddRWZmZk88MAD3HDDDXNUrYhbahoSzxgeHqahoSFwu8Xq6mrKy8vZuXMne/bs4dChQ4FtFy1aRHl5ucNqRWaXmobEM+rr69mzZw/79+/nwIEDlJeXU1VVRU1NDXV1dQwNDY3aPioqimXLlvHud7+bdevWUVxcTGFhIe9617scfQKRuacgkHmtp6eHDRs28Pjjj/Pss89SU1MTWBcbG8uiRYsoKChg1apVgWGn/iGo+fn55Ofn6/pG4nkKAplXWltbefbZZ9myZQuvvfYar7zyCt3d3aSmpnLFFVdw5plnsmLFCpYuXUphYaHuGyAyAwoCmVcuueQS/H1ESUlJ3HLLLVx//fWsXbtWv+xF3iYFgcwr//iP/xgIgo6ODh566CE2b97MypUrWblyJatWrWLlypXk5ubqpDGRGdKoIZl3+vv72b17Nzt27GDnzp2BKbh/ICMjg5UrV3LOOefwgQ98gNNPP91hxSLu6cxi8YTGxkZ27drFzp072bFjBzt27GDr1q0MDQ3xxS9+kfvuu891iSLOaPioeEJmZiZr165lzZo1bNq0iezsbCorK6murmb//v2uyxMJWwoCOSEcPXqUZ555hj/84Q88++yz9PT0EBkZyZo1a7j77rtZt26d6xJFwpaCQOat7u5u/v3f/51f/vKX7Nu3D4CioiJuvfVWLrvsMi688EKSk5MdVykS/hQEMq9Yazl06BBPP/00X//616mrq+Piiy/mtttu45JLLmHFihUaLSRynBQEEtaamprYvHkzmzZtYtOmTbz++uu0trYCcNZZZ/Hb3/6W8847z22RIvOcgkDC0h//+Eduv/12Dh8+DIzcgObUU0/lgx/8IKeffjrnnnuufv2LzBIFgYSlBx98MBACH/3oR/ne975HSkqK26JETlC6EIuEpYceeoibbroJgJ/+9KfcddddjisSOXHpiEDCQk9PDy+99BLPP/88mzZt4pVXXqGzsxMYOT/guuuuc1ugyAlMQSBO1dTUcOedd/LII4/Q19cX6Au45ZZbOOecczj77LMpKSlRX4BICCkIxKmLLrqIffv28b73vY9PfepTnH/++SQlJbkuS8RTFATiRFtbG3v27CE1NRWA+Ph4rrzySrdFiXiUgkBCqqOjg71797J79+7AtGvXLiorKwPbLFmyhE984hMOqxTxNgWBvGPWWqqqqti/f3/gXsH79u1jz549HDt2LLBdTEwMy5Yt48ILL6SsrIyysjJWrFihPgARxxQEMmPt7e2BL/rgL/0DBw7Q3d0d2C4hIYElS5ZwwQUXsHz5cpYtW8aKFStYtGgRkZGRDj+BiExEQSDjNDc3B5pwdu3axe7du9m/fz+1tbWBbSIiIigtLWXJkiWsXbuWpUuXsmTJEpYuXUp+fr5+4YvMIwoCoa+vj2effTYw7d27N7AuOTmZFStWcMUVV7B06dLAF/6iRYuIiYlxWLWIzBYFgce98MILfOxjH6OiooLY2FguuOAC1q1bx2mnnUZZWRkLFy7Ur3uRE5yCwONuvvlmIiIiePLJJ7n00kuJjY11XZKIzDFda8jjmpubueGGG7jmmmsUAiIepSDwuISEBLq6ulyXISIOKQg8bHBwkMbGRvLz812XIiIOKQg8rLGxESBwmQcR8SYFgYf5b/mYmZnpthARcUpB4GH+JqGqqirHlYiISwoCDxscHATQZR9EPE5B4GHl5eUAlJaWOq5ERFxSEHhYcnIygIaPinicgsDD/EcChw8fdluIiDilIPCwQ4cOAVBcXOy4EhFxSUHgYe3t7QBkZGQ4rkREXFIQeNjw8DAwcm8BEfEufQN4mD8ArLWOKxERlxQEwsDAgOsSRMQhBYGH+fsI4uLiHFciIi4pCDzMfx5BZ2en40pExCUFgYdFR0cDOqFMxOsUBB7mv8TEihUrHFciIi4pCDzMfySg+xGIeJuCwMNyc3MBqKiocFyJiLikIPAw/83qm5ubHVciIi4pCDysra0N0LWGRLxOQeBh/j6ClJQUx5WIiEsKAg/zn0fgPzIQEW9SEHiYgkBEQEHgaTU1NcBbo4dExJsUBB7mv+roggULHFciIi4pCDxs0aJFAOzbt89xJSLikoLAw/z3I/DfoEZEvElB4GF1dXUA5OfnO65ERFxSEHiY+ghEBBQEnuY/EqiqqnJciYi4pCDwMN2zWERAQeBpHR0dACQlJTmuRERcUhB4WGRkJAB9fX2OKxERlxQEHrZy5UoiIyN59dVXXZciIg4pCDwsISGBpKQkGhoaXJciIg4pCDxueHhYJ5SJeJyCwOPOOOMMdu3a5boMEXFIQeBxCxcupLa21nUZIuKQgsDjsrOzqa+vd12GiDikIPC46Oho+vr6dFKZiIcpCDwuLi6OoaEhhoaGXJciIo6ELAiMMQ8ZY+qNMRP2RBpj1hpj2owx23zT3aGqRSaXlZUFoCGkIh4WystO/hT4AfDzKbZ5wVp7VQhrkGn4jwR6e3sdVyIiroTsiMBauxFoDtXry+xoaWkBoKioyHElIuKK6z6Cs40x240xfzLGrJhsI2PMbcaYLcaYLWrCmF0JCQmAmoZEvMxlELwOFFtrTwX+E3hisg2ttQ9aa1dba1f727RldpSWlgJQXV3tuBIRccVZEFhr2621nb75p4EoY0ymq3q8Kjo6GoDBwUHHlYiIK86CwBiTa4wxvvmzfLU0uarHq9ra2gBISUlxXImIuBKyUUPGmIeBtUCmMaYS+BoQBWCtfQD4APBJY8wg0APcZHVWkzM9PT2uSxARR0IWBNbam6dZ/wNGhpeKQ2VlZQDs3r2b0047zW0xIuKE61FD4pj/EtT+vgIR8R4Fgce1trYCkJaW5rYQEXFGQeBxhYWFABw6dMhxJSLiioLA4/znERw9etRxJSLiioLA47q7uwFITk52XImIuKIg8Dh/Z7HvlA4R8SAFgcclJiZijKG9vd11KSLiiILA44wxWGt1RCDiYQoCj4uIiCAmJkb3IxDxMAWBkJCQQFdXl+syRMQRBYGQnp5ObW2t6zJExBEFgXDyySdTXl7uugwRcURBIBw6dIiCggLXZYiIIwoCYXh4mKioKNdliIgjCgKhqKiIqqoq12WIiCMKAmFoaEi3qhTxMAWB0Nraqj4CEQ9TEAitra26Z7GIhykIhJ6eHmJiYlyXISKOKAiEvLw8dRaLeJiCQBgaGtI9i0U8TEEgdHV1kZCQ4LoMEXFEQSAkJibS1tbmugwRcURBIMTExOg8AhEPUxAIjY2NZGRkuC5DRBxREIj6CEQ8TkEg5OXlUV9f77oMEXFEQSCkpKTQ1NTkugwRcURBIBw+fJjCwkLXZYiIIwoCobu7W5eYEPEwBYGQl5en8whEPExBICQkJNDe3u66DBFxREEgNDc3k5qa6roMEXFEQSD09vYSFxfnugwRcURBIERFRTE0NOS6DBFxREEgxMXFqY9AxMMUBB5nrWX//v0sWbLEdSki4oiCwOOMMaSkpOgOZSIepiAQLrzwQrZt2+a6DBFxREEg7N+/n+LiYtdliIgjCgJhz549rF692nUZIuKIgkDIzMyksbHRdRki4oiCQACIjIx0XYKIOKIgEKKioujq6nJdhog4oiAQYmNjGRgYcF2GiDiiIBCGh4eJjo52XYaIOKIgEPr7+4mKinJdhog4oiAQ+vv7dUQg4mEKAqGnp0e3qhTxMAWBYK11XYKIOKQgELKysqipqXFdhog4oiAQIiIiiIjQPwURr9L/fqGvr0+dxSIepiAQrLU6IhDxMP3vFzo6OkhISHBdhog4oiAQ4uPj6enpcV2GiDiiIPC4gYEB6uvrKSwsdF2KiDiiIPC4/v5+rLXExsa6LkVEHFEQeFxCQgIJCQnU19e7LkVEHFEQeFxPTw8DAwMsWLDAdSki4oiCwONaWlro7+9XH4GIhykIPC43N5ekpCT27dvnuhQRcURB4HERERHExcXR39/vuhQRcURBIHR0dJCcnOy6DBFxREEgDA4O6g5lIh6mIBDS09M1fFTEwxQEQnp6Oh0dHa7LEBFHFARCXFycrjUk4mEKAiE2Npbe3l7XZYiIIwoCYcGCBQwNDbkuQ0QcURAICxYsYGBgwHUZIuKIgkCIiopSEIh4mIJAdEQg4nEKAiE6OlpBIOJhCgJR05CIxykIhMjISI0aEvEwBYEQERHB8PCw6zJExBEFgRAZGakgEPEwBYEQERGhpiERD1MQiI4IRDxOQSA6IhDxOAWB6IhAxOMUBKJRQyIepyAQNQ2JeJyCQNQ0JOJxCgLREYGIxykIREcEIh6nIBAdEYh4nIJAdEQg4nEKAtHwURGPUxCImoZEPE5BIEREjPwzsNY6rkREXFAQCMYYQEEg4lUKAlEQiHicgkACTUPqMBbxJgWB6IhAxOMUBKLOYhGPUxBI4IhATUMi3qQgEB0RiHicgkB0RCDicQoCUWexiMcpCERNQyIepyAQNQ2JeJyCQHREIOJxCgLREYGIxykIRJ3FIh6nIBA1DYl4nIJA1DQk4nEKAtERgYjHhSwIjDGFxpi/GWP2GGN2G2M+N8E2xhjzfWNMuTFmhzHmjFDVI5PTEYGIty0I4WsPAl+w1r5ujEkCthpjnrXW7gna5gpgsW96N/Aj31+ZQzoiEPG2kB0RWGtrrLWv++Y7gL3AwjGbXQv83I7YDKQaY/JCVZNMTEcEIt42J30ExpgS4HTglTGrFgLHgh5XMj4sMMbcZozZYozZ0tDQELI6vUrDR0W8LeRBYIxJBH4H/LO1tv3tvIa19kFr7Wpr7eqsrKzZLVDUNCTicSENAmNMFCMh8Ctr7eMTbFIFFAY9LvAtkzmkpiERbwvlqCED/ATYa6397iSbPQWs840eWgO0WWtrQlWTTExHBCLeFspRQ+cCHwF2GmO2+ZbdBRQBWGsfAJ4GrgTKgW7gYyGsRyahIwIRbwtZEFhrXwTMNNtY4NOhqkFmRp3FIt6mM4tFTUMiHqcgEDUNiXicgkB0RCDicQoC0RGBiMcpCERHBCIepyAQHRGIeJyCQDR8VMTjFASipiERj1MQiJqGRDxOQSA6IhDxOAWB6IhAxOMUBKLOYhGPUxCImoZEPE5BIGoaEvE4BYHoiEDE4xQEoiMCEY9TEIg6i0U8TkEgahoS8TgFgahpSMTjFARCQkICp5xyCnFxca5LEREHQnbzepk/zjrrLPbu3eu6DBFxREcEIiIepyAQEfE4BYGIiMcpCEREPE5BICLicQoCERGPUxCIiHicgkBExOMUBCIiHqcgEBHxOAWBiIjHKQhERDxOQSAi4nEKAhERj1MQiIh4nIJARMTjFAQiIh6nIBAR8TgFgYiIxykIREQ8TkEgIuJxCgIREY9TEIiIeJyx1rqu4bgYYxqAI2/jqZlA4yyXEyqqNTRUa2io1tkXijqLrbVZE62Yd0HwdhljtlhrV7uuYyZUa2io1tBQrbNvrutU05CIiMcpCEREPM5LQfCg6wKOg2oNDdUaGqp19s1pnZ7pIxARkYl56YhAREQmoCAQEfG4eRsExphCY8zfjDF7jDG7jTGf8y1PN8Y8a4x50/c3zbfcGGO+b4wpN8bsMMacEfRat/i2f9MYc8sc1vodY8w+Xz3rjTGpQc/5iq/W/caYy4KWX+5bVm6MuXOuag1a/wVjjDXGZPoeh91+9a37jG/f7jbG3Be0PKz2qzHmNGPMZmPMNmPMFmPMWb7lLvdrrDHmVWPMdl+t9/iWlxpjXvHV9IgxJtq3PMb3uNy3viTotSbc33NQ669877nLGPOQMSbKtzzs9mvQ+u8bYzqDHs/dfrXWzssJyAPO8M0nAQeA5cB9wJ2+5XcC3/bNXwn8CTDAGuAV3/J04KDvb5pvPm2Oar0UWOBb/u2gWpcD24EYoBSoACJ9UwVwEhDt22b5XNTqe1wIPMPICX2ZYbxf3wP8BYjxrcsO1/0KbACuCNqXz4XBfjVAom8+CnjFV8OjwE2+5Q8An/TNfwp4wDd/E/DIVPt7jmq90rfOAA8H1Rp2+9X3eDXwC6AzaPs526/z9ojAWltjrX3dN98B7AUWAtcCP/Nt9jPgOt/8tcDP7YjNQKoxJg+4DHjWWttsrW0BngUun4tarbUbrLWDvs02AwVBtf7GWttnrT0ElANn+aZya+1Ba20/8BvftiGv1bf634EvAcEjDMJuvwKfBL5lre3zrasPqjXc9qsFkn2bpQDVQbW62q/WWuv/ZRrlmyxwEfBb3/Kx/7f8/+d+C1xsjDFMvr9DXqu19mnfOgu8yuj/W2G1X40xkcB3GPm/FWzO9uu8DYJgvkOm0xlJ2BxrbY1vVS2Q45tfCBwLelqlb9lky+ei1mC3MvJLhSlqclarMeZaoMpau33MZmFXK7AEON93OP28MebMMK71n4HvGGOOAfcDXwmHWo0xkcaYbUA9I1+KFUBr0A+X4PcN1ORb3wZkuKrVWvtK0Loo4CPAn8fWOqYml7XeDjwV9L3lN2f7dd4HgTEmEfgd8M/W2vbgdb5fA2EzPnayWo0xXwUGgV+5qm2s4FoZqe0u4G6XNU1mgv26gJFD/DXAF4FHfb+knJug1k8Cn7fWFgKfB37isj4/a+2QtfY0Rn5JnwWc4raiyY2t1RhTFrT6v4CN1toXnBQ3xgS1XgDcCPyny7rmdRD40v53wK+stY/7Ftf5DvXw/fU3C1Qx0sbtV+BbNtnyuagVY8xHgauAD/mCKxxrXcRIW+R2Y8xh3/u+bozJDcNaYeQX0uO+Q/FXgWFGLuIVjrXeAvjnH+OtQ3yntfpZa1uBvwFnM9KMsmCC9w3U5FufAjQ5rPVyXy1fA7KAO4I2C7f9+h7gZKDc938r3hhTPrbWkO/XmXQkhOPESMfLz4HvjVn+HUZ3Ft/nm38fozuJXrVvdRIdYqSDKM03nz5HtV4O7AGyxixfwejOoIOMdGgu8M2X8lan5oq5qHXMNod5q7M4HPfrPwH3+uaXMHIYbcJxvzLSV7DWN38xsDUM9msWkOqbjwNeYOTHymOM7iz+lG/+04zu1Hx0qn/Hc1Trx4GXgbgx24fdfh2zTXBn8Zzt11n7kHM9Aecx0uyzA9jmm65kpA3tf4E3GRk5ku7b3gA/ZKStcyewOui1bmWkw6Uc+Ngc1lrOyJeUf9kDQc/5qq/W/fhGlfiWX8nIiJMK4KtzVeuYbQ7zVhCE436NBn4J7AJeBy4K1/3qW77V9x/7FeBdYbBfVwFv+GrdBdztW34SIx2v5YyEgn9UVqzvcblv/UnT7e85qHXQ977+fe1fHnb7dcw2wUEwZ/tVl5gQEfG4ed1HICIi75yCQETE4xQEIiIepyAQEfE4BYGIiMcpCEREPE5BICLicf8PoStW/CE8ZM8AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Plot up the raw stuff\n", "fig, ax = plt.subplots(figsize=(6,10))\n", "ax.plot(df['Sensor3'], df['time'], color='k')\n", "ax.set_title(\"Raw Active NIR\")\n", "ax.set_ylabel('Time [s]')\n", "a = ax.invert_yaxis()\n" ] }, { "cell_type": "markdown", "id": "a335f213", "metadata": {}, "source": [] }, { "cell_type": "code", "execution_count": 15, "id": "dac2911f", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmkAAAFzCAYAAABl1J6yAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABmr0lEQVR4nO3dd3hc5Zk3/u89RVPUrWJLltyxAVewMdVg0xNCWyDAD/KG7AZSIcm+SyCb3WySd9M22YSELAlsQggEHIgTCAmEFjDNYFywcQEbd8uSrWZ1aerz++PoTJGmnClnZiR9P9fly9LMmZnHluX56r6fIkopEBEREVFhseR7AEREREQ0GkMaERERUQFiSCMiIiIqQAxpRERERAWIIY2IiIioADGkERERERUgW74HYIbq6mo1Y8aMfA+DMrCrYxcAYF7VvDyPhIhG2aV9f2Ievz+JsmHTpk3tSqmakbePy5A2Y8YMbNy4Md/DoAysfGglAGDtLWvzOg4iimHlSu33tWvzOQqicUNEDsa6ne1OIiIiogLEkEZERERUgBjSiIiIiArQuJyTRkRERGE+nw9NTU0YGhrK91AmNKfTiYaGBtjtdkPXM6QRERGNc01NTSgtLcWMGTMgIvkezoSklEJHRweampowc+ZMQ49hu5OIiGicGxoaQlVVFQNaHokIqqqqUqpmMqQRERFNAAxo+Zfq14AhjYiIiPLinnvuwcDAQMqPe+ihh9Dc3GzCiAoLQxoRERHlRTohLRAIMKQRERERZUt/fz8uu+wyLF68GAsWLMC3vvUtNDc3Y9WqVVi1ahUA4HOf+xyWLVuG+fPn4z/+4z9Cj50xYwbuuusunHrqqVi9ejU2btyIm266CUuWLMHg4GC+/kim4+pOIiKiCeTLz30ZW45uyepzLpmyBPdcek/Ca5577jnU19fjmWeeAQB0d3fjN7/5DV555RVUV1cDAL7zne9g0qRJCAQCuOCCC/Dee+9h0aJFAICqqips3rwZAPCrX/0KP/rRj7Bs2bKs/jkKDStpRAWqv9+LvXu78j0MIqKsWLhwIV588UXcddddeP3111FeXj7qmieeeAKnnnoqTjnlFOzYsQM7d+4M3Xf99dfncrgFgZU0ogJ1660vYPXqD+D1fgV2uzXfwyGicSJZxcssc+fOxebNm/Hss8/i3/7t33DBBRdE3b9//3786Ec/woYNG1BZWYlbbrklaruK4uLiXA8571hJIypQq1d/AABoaxu/8y2IaOJobm6G2+3GzTffjDvvvBObN29GaWkpent7AQA9PT0oLi5GeXk5jh07hr/97W9xnyvyceMZK2lEBa6vz5vvIRARZWzbtm248847YbFYYLfb8Ytf/AJvvfUWLr30UtTX1+OVV17BKaecghNPPBGNjY04++yz4z7XLbfcgs9+9rNwuVx466234HK5cvgnyR2GNKICNzjoz/cQiIgydskll+CSSy6Jum3ZsmW4/fbbQ58/9NBDMR974MCBqM+vueYaXHPNNdkeYsFhu5OowDGkERFNTAxpRAWOIY2IaGJiSCMqcAxpREQTE0MaUYHzeAL5HgIREeUBQxpRgfP5GNKIiCYihjSiAuf3q3wPgYiI8oAhjajAsZJGROPFd77zHcyfPx+LFi3CkiVLsH79etxzzz0YGBjI99AKEvdJIypwPl8w30MgIsrYW2+9hb/+9a/YvHkzHA4H2tvb4fV6cf311+Pmm2+G2+3O9xALDitpRAXO72dII6Kxr6WlBdXV1XA4HACA6upqrFmzBs3NzVi1ahVWrVoFAFi9ejUWLlyIBQsW4K677go9vqSkBF/5ylcwf/58XHDBBWhra8vLnyOXWEkjKnCspBFRVn35y8CWLdl9ziVLgHvuSXjJxRdfjG9/+9uYO3cuLrzwQlx//fW444478OMf/xivvPIKqqur0dzcjLvuugubNm1CZWUlLr74Yjz11FO46qqr0N/fj2XLluEnP/kJvv3tb+Nb3/oWfv7zn2f3z1FgWEkjKnAMaUQ0HpSUlGDTpk144IEHUFNTg+uvv37UMVAbNmzAypUrUVNTA5vNhptuugmvvfYaAMBiseD6668HANx888144403cv1HyDlW0ogKHNudRJRVSSpeZrJarVi5ciVWrlyJhQsX4re//W3azyUiWRxZYWIljajAsZJGROPBrl278OGHH4Y+37JlC6ZPn47S0lL09vYCAJYvX45XX30V7e3tCAQCWL16Nc477zwAQDAYxJo1awAAjz32GM4555zc/yFyjJU0ogIlAijFLTiIaHzo6+vD7bffjq6uLthsNsyZMwcPPPAAVq9ejUsvvRT19fV45ZVX8P3vfx+rVq2CUgqXXXYZrrzySgBAcXEx3nnnHfznf/4namtr8fjjj+f5T2Q+hjSiAqSUghrew5aVNCIaD5YuXYp169aNuv3222/H7bffHvr8xhtvxI033hjzOX784x+bNr5CxHYnUQFSEYcMcE4aEdHExJBGVIACgXAwYyWNiEhrl040DGlEBSgQCJfSGNKIiCYmhjSiAhQMhkMa251ERBMTQxpRAWIljYiIGNKIClDknDRW0oiIJiaGNKICFNnuZCWNiMaLp556CiKCDz74IGvPeeDAASxYsCBrz/fd73436vOzzjora8+dKoY0ogIU2e5kJY2IxovVq1fjnHPOwerVq/M2Br/fn/D+kSEt1t5uucKQRlSAOCeNiMabvr4+vPHGG/j1r3+N3//+9wCAQCCAf/mXf8GCBQuwaNEi3HvvvQC0g9bPOussLF68GMuXL0dvby8CgQDuvPNOnHbaaVi0aBHuv//+Ua8R75q1a9dixYoVuOKKK3DyyScDAK666iosXboU8+fPxwMPPAAAuPvuuzE4OIglS5bgpptuAqAdDA9om4zfeeedWLBgARYuXBg68WDt2rVYuXIlrr32Wpx44om46aaboCI3u8wATxwgKkCck0ZEZvnyl1/Gli2tWX3OJUtqcc895ye85s9//jMuvfRSzJ07F1VVVdi0aRPeeecdHDhwAFu2bIHNZkNnZye8Xi+uv/56PP744zjttNPQ09MDl8uFX//61ygvL8eGDRvg8Xhw9tln4+KLL446aD3eNQCwefNmbN++HTNnzgQAPPjgg5g0aRIGBwdx2mmn4ZprrsH3v/99/PznP8eWLVtGjf9Pf/oTtmzZgq1bt6K9vR2nnXYazj33XADAu+++ix07dqC+vh5nn3023nzzzaycLcqQRlSAouek8exOIhr7Vq9ejS996UsAgBtuuAGrV6/G/v378dnPfhY2mxZHJk2ahG3btqGurg6nnXYaAKCsrAwA8MILL+C9994LHbLe3d2NDz/8EHPnzg29RrxrioqKsHz58lBAA4Cf/exnePLJJwEAhw8fxocffoiqqqq443/jjTdw4403wmq1YvLkyTjvvPOwYcMGlJWVYfny5WhoaAAALFmyBAcOHGBIIxqvouekZadsTkQEIGnFywydnZ14+eWXsW3bNogIAoEARCQUxIxQSuHee+/FJZdcEnX7gQMHkl6zdu1aFBcXR33+0ksv4a233oLb7cbKlSsxNDSU3h8OgMPhCH1stVqTznszinPSiApQ9Jw0VtKIaGxbs2YNPvGJT+DgwYM4cOAADh8+jJkzZ2Lx4sW4//77Q6Gms7MT8+bNQ0tLCzZs2AAA6O3thd/vxyWXXIJf/OIX8Pl8AIDdu3ejv78/6nWMXANoFbbKykq43W588MEHePvtt0P32e320OMjrVixAo8//jgCgQDa2trw2muvYfny5dn5C4ojryFNRC4VkV0iskdE7o5x/y0i0iYiW4Z/fTof4yTKtegTB1hJI6KxbfXq1bj66qujbrvmmmvQ0tKCadOmYdGiRVi8eDEee+wxFBUV4fHHH8ftt9+OxYsX46KLLsLQ0BA+/elP4+STT8app56KBQsW4DOf+cyoipWRawDg0ksvhd/vx0knnYS7774bZ5xxRui+2267DYsWLQotHNBdffXVoXGef/75+K//+i9MmTIli39Lo0m2ViCk/MIiVgC7AVwEoAnABgA3KqV2RlxzC4BlSqkvpvLcy5YtUxs3bsziaCnXVj60EgCw9pa1eR1HvnzwQQdOOuk3AIDly6dg/fqb8zwioggrV2q/r12bz1FQCt5//32cdNJJ+R4GIfbXQkQ2KaWWjbw2n5W05QD2KKX2KaW8AH4P4Mo8joeoYHBOGhER5TOkTQVwOOLzpuHbRrpGRN4TkTUi0hjvyUTkNhHZKCIb29rasj1Wopzi6k4iIir0hQN/ATBDKbUIwIsAfhvvQqXUA0qpZUqpZTU1NTkbIJEZ9EqazWbhPmlERBNUPkPaEQCRlbGG4dtClFIdSinP8Ke/ArA0R2Mjyit9M1un08oTB4goK/I1B53CUv0a5DOkbQBwgojMFJEiADcAeDryAhGpi/j0CgDv53B8RHmjtzsdDhsraUSUMafTiY6ODga1PFJKoaOjA06n0/Bj8raZrVLKLyJfBPA8ACuAB5VSO0Tk2wA2KqWeBnCHiFwBwA+gE8At+RovUS7p7U6Hg5U0IspcQ0MDmpqawDnb+eV0OkMnExiR1xMHlFLPAnh2xG3fiPj4awC+lutxEeVbZEgbGBi9qSIRUSrsdnvUkUg0NhT6wgGiCUmfk8ZKGhHRxMWQRlSAwnPSrJyTRkQ0QTGkERUgzkkjIiKGNKICpIc0p5OrO4mIJiqGNKICFNnu9PmCXDZPRDQBMaQRFaDIhQNA9DFRREQ0MTCkERWgyDlpADgvjYhoAmJIIypAkScOAOC8NCKiCYghjagAjaykMaQREU08DGlEBSjygHWA7U4ioomIIY2oAIUraWx3EhFNVAxpRAUocgsOgJU0IqKJiCGNqABxThoRETGkERWg8Jw0rd3JShoR0cTDkEZUgILDmSxeJa29fQC33fYCjh7tz/XQiIgoRxjSiArQyBMHRlbSHn30ffzv/76HBx7YmvOxERFRbjCkERWgZHPSmpp6h3/vy+3AiIgoZxjSiEw0MODD6af/Dn/84+6UHpdsdefBgz0AgNbWgSyMkoiIChFDGpGJduxoxzvvHMXtt/89pcclq6S1tQ0O/86QRkQ0XjGkEZmos3MIQOqrM5Ot7uzr8wJgJY2IaDxjSCMykR7SlErtcSPbnSMraX19PgAMaURE4xlDGpGJ9JBWVJTat9rIdqfPF4i6v79fC2k9PV54PP5Mh0lERAWIIY3IRF1dHgCA3Z5uSNPP7owuxfX1+UIBTp+fRkRE4wtDGpGJvF6tApbunLR4lbS+Pi9mziwHwMUDRETjFUMakYn0kDYwkFpLcvSctHAlzesNwOcLYtYsLaRxXhoR0fjEkEZkIj2k6XPIjEo0J01f2TlrVgUAhjQiovGKIY3IRF6v1rb0+4OhwGZEIKBgsUhoLlvk6k59ZWe43ck5aURE4xFDGpGJIueipVJNCwYVrFaBzWYZ9Tx6Ja2+vgR2u4WVNCKicYohjchEkdWzwUHj89ICgSAslnBIi1VJKy21o6bGzYUDRETjFEMakYkiQ9rAgPFKWiCgVdL0dmd0JU17npKSItTWullJIyIapxjSiEwUHdKMV9K0dqclTiVNa3eWlNgZ0oiIxjGGNCIT6QsHgNQraZELB+JV0mpqXFw4QEQ0TjGkEZkokzlpkQsHWEkjIpp4GNKITOTzBVFUpO11lkolLfHqzshKmhv9/b6UnpuIiMYGhjQiE3m9AZSXFwFIbU6a3u4UEVitErOSVlysVdIAHg1FRDQeMaQRmUgLaQ4A6azu1L497XbrqC04nE4bbDZLKKQdO8aQRkQ03jCkEZnI6w2iokIPaanPSQMAm01GbYpbUmIHAMycWQYA2LDhKG6//e/48MPj2Ro6ERHlmS3fAyAaz7zeAKZMKQaQ2sIBfU4aEKuS5g2FtHnzJsFut+CLX/w7AKCry4NHHvlotoZPRER5xEoakYmi56SlvgUHMLqS1tfnQ3GxFtLsditOPHFS6L7Nm49lY9hERFQAWEkjMpHXG0BxsR12uyXlhQOJ5qTpIQ0AfvnLi/Dcc/vh8QTw3/+9ET5fAHa7NXt/CCIiygtW0ohMpG/B4XLZ0tqCA4g9Jy0ypJ111lR8+9vn4IQTKhEIKDQ392XvD0BERHnDkEZkIq83ALvdArfbnvLCAb3dObKSNjKk6RobSwEAhw/3ZjhqIiIqBAxpRCbyerVKmtudWiVNP2AdAGw2S8JKmo4hjYhofGFIIzKR1xsYDmn2NFZ36nPSLIYqafX1JQDAdicR0TiR15AmIpeKyC4R2SMid8e43yEijw/fv15EZuRhmERpC4e01Ctp4dWd0ZW0gYHYIa283AGbzYKOjqHMB05ERHmXt5AmIlYA/wPgIwBOBnCjiJw84rJ/AnBcKTUHwE8A/CC3oyRKXzCoEAgoFBWlNyctvE+asUqaiKC62sUjooiIxol8VtKWA9ijlNqnlPIC+D2AK0dccyWA3w5/vAbABSIiORwjUdp8vgAAbeJ/qqs7481J8/kC8PmCMUMaAFRXu9DePpjhyImIqBDkM6RNBXA44vOm4dtiXqOU8gPoBlCVk9ERpUgpBaVU6HOvVwtWWiXNllIlLRhUEas7w5W0/n4t6LndsUNaTQ1DGhHReDFuFg6IyG0islFENra1teV7ODQB3XLL33DppX8Mfe71apW0eAsHlFKha0aK3Mw2spKmh7RElbS2NoY0IqLxIJ8h7QiAxojPG4Zvi3mNiNgAlAPoiPVkSqkHlFLLlFLLampqTBguUXzbt7fh4Yd34oUXDqCvzwtgdEgb2e78l39Zi9ra+3D4cM+o54tsd8aqpMULaTU17qhKmsfjj5rPRkREY0c+Q9oGACeIyEwRKQJwA4CnR1zzNIBPDn98LYCXVWQ/iahA/OAHG0If6/uUhUPa6HbnSy8dxI9/vAnd3R788Y8fjnq+6BMHUqukdXYOIhAIwu8P4pRTHsZll/0pC39CIiLKtbyFtOE5Zl8E8DyA9wE8oZTaISLfFpErhi/7NYAqEdkD4J8BjNqmgyjf9u/vwurV72P58ikAIkOaPictXEnTf8b4+tdfx6xZ5aitdePtt5tHPWf0iQOpVNJcUAro7BzCzp0deP/9TrzwwgF0drIFSkQ01uR1TppS6lml1Fyl1Gyl1HeGb/uGUurp4Y+HlFLXKaXmKKWWK6X25XO8RLH86EcbYbEIvvvdFQCAjg4tEIVXd1rgctkQCCj4fEEcPdqPd945ittuW4SlSydj9+7jo55z9OpO7bmSV9LcAIC2tgHs3dsVuv299zhPk4horBk3CweI8uHYsX48+OB2fPKT8zF/fjUA4PhxbTPZ6EqaDYC2Ee1bb2mVsxUrGjB7dgX27u3CyC5+5IkDNpsFfr8afrzWMo0X0iZPdg+PawB79oTD34EDo+e9ERFRYWNII8rAT3+6GR6PH3feeRoqKx0AgOPHPQBGLxwAgMFBPzZvPgarVXDqqZMxe3YFenq8oeqbLvLEAbvdeCUtMqTt3duNkhI7RIADB7qz+ccmIqIcYEgjSlNfnxf33bcF11wzF3PnToLDYYPLZUNXl15JiwxpeiXNj/37u9HYWAqn04aGBv28zf6o5448cSCykpY8pBUD0Cp8e/d24aSTqjB1aikraUREYxBDGlGaHnlkJ7q7Pfjnf14Wuq2iwhFRSYvczFYLVQMDPhw40IMZM8oBAHV1WkhraYk+FD1ydWesSpoe+kaqrHTCZrMMV9K6MHt2BWbMKGMljYhoDGJII0rTb3+7A0uW1OKMM+pCt5WWFsXZJy1cSTtwoBszZpQBAOrqtMpXc3N0SBt5wLpeSevp8YReJxaLRVBb60ZTUy8OHuwZDmnl2L+fIY2IaKxhSCNKw9Gj/Vi/vgXXXHMCIo+TLS62h6pd+t5m+upOAOjqGkJzc19EJU0LaS0tI9ud4YUDkZW07m4vSkrsoftiqa114513jiIYVJgzpwIzZ5ahqakv9BxERDQ2MKQRpeGZZ7TdYK64Yk7U7W53OKTFWjiwa9dxKIVQJc3lsqOiwhEjpI2ck6YFvp4eD8rLHQnHNnmyG7t2dQIAZs+uwMyZ5QgGVWj/NiIiGhsY0ojS8NJLB1FXV4yFC6ujbi8utoe2yYjV7ty5UzvVTK+kAVo1bWS7c+QB63pVrrvbg7Ky2K1Onb7CE0Co3QmALU8iojGGIY0oRUopvPLKIaxaNS2q1QlEtzujj4XSKmnhkFYWekxdXUmcdqf23A6HFT5fEMGgQne310AlTWuhOp021NUVY+ZMLaRxhaf5Dh7sxsqVv8emTUfzPRQiGgcY0ohStGtXJ44dG8CqVY2j7tMqaXpIiz4WCtBCmtUqmDq1NPSYWJW0yDlp+ny2wUGfoXbn/PlVAACn0woRQUNDKaxWYSUtBx55ZCdefbUJ//RPz4/aoJiIKFUMaUQpWrv2MABg5crRIc3tto2qpNntllC7s6NjEI2NpbDZwt96U6YU49ixgag39cgtOFyu8Ea4RippF188A5WVTtx993IA2py2xsZShrQc2Lq1LfT7mjW78zwaIhrrYm+2RERxvfLKYTQ0lGL27IpR98Va3VlUZEVxsR0WiyAYVFHz0QDtUPShIT/6+30oKdHmm0UesB6upPkNzUmrqytBa+vno4LgzJnl3CstB44c6cN55zWgrW0QP/zhBlx33bx8D4mIxjBW0ohS4PcH8fe/H8KqVY2j5qMB4YUDSqmohQNWqwW1tdqE/unTy6IeU1MTPhRdFzknLTKkdXV5UFGRuJIGICqgAdpChT17uvDnP+/BDTf8JdSSpexqbu5DY2MZPvGJk7Fhw1EcPsx5gESUvrghTUR6kvzqFRHW82lCef31JnR0DOLKK+fEvN/ttiMYVPB4AlEhDQAmTXICQGgivy4c0sLnd0YesK6HtNbWAQwN+UMLA1Jx+ul1aG0dwFVXPYXHH9+FJ5/8MOXnoMSCQYXm5j5MnVqCyy6bBUCruhIRpStRJW2vUqoswa9SAP0JHk807jz00Ha4XDZceumMmPfrZ2r29/vg9QYgglBFTG9/Rq7sBLR2JzC6kjay3amvztQrcqm47LJZKCqyorJSC4rr1jWn/ByUWFfXEHy+IKZMKcb8+dUoL3fgzTeP5HtYRDSGJQpp1xh4vJFriMaFPXuO49FH38dnP7sYxcWx54XpCwQGB/3weoMoKrKG2qLnnz8NADB/fvTearEqaZGb2eoh7eDB9ENaQ0MpNmy4Ge+++wmcc85UbNvWnvJzUGI9PdpxYGVlRbBYBGeeWccwTEQZiRvSlFL7AEBEikXEMvzxXBG5QkTskdcQTQTf+c7bsNut+OpXl8e9JnL+mNcbgN0e/hb7/vdX4A9/uBxLl06OekysSlowiKyGNABYtKgG06eXY9YsLiIwQ2+vFtL0c1WXL6/Dzp0dnP9HRGkzsnDgNQBOEZkK4AUAnwDwkJmDIioUwaDCoUM9+POf9+Dhh3fic59bjClT4s8J0/dDGxjwwecLhuajAUBFhRPXXjtv1IKD0tIiiIQrMUDs1Z1bt7YCAGbNip7TlqqZM8vR1NQbmjNH2TEypC1ZUotgULFqSURpMxLSRCk1AOAfANynlLoOwHxzh0WUfz09Hixf/jtMn/4ArrrqKTQ0lOAb3zgz4WNGVtIiQ1o8IhJ15icQezPbd945imnTSlFR4Uz3jwRAW12qFNDUxLM8s2lkSDvllFoAwJYtrXkbExGNbUb2SRMRORPATQD+afi25O88RGPcL36xBZs2HcP3vrcC06aV4YILpiUNSHqgGhjQQ5qxXW5KSsIhTSk1YjPb8LfpwoU16fxRouiVwGPHBjBrVkXGz0eakSFt+vQyVFQ48O67DGlElB4jIe3LAL4G4Eml1A4RmQXgFVNHRZRnwaDC/fe/h/POa8Ddd59u+HHhhQO+0MIBIyI3wdUPHtBDWllZeF+0RYuyF9KOHuXi7GwKhzSt5S0iWLKkFu++eyyfwyKiMSzpj/lKqVeVUlcopX4w/Pk+pdQd5g+NKH9eeOEA9u/vxuc+tySlx42upBkPaX19WkgLBLStOvQ5aZEnDGQjpOn7rB07xpCWTb292tdPr6QBWsvzvffa4fcH8zUsIhrDEm1m+xsReVBEfpLLAREVgl/8Ygtqa924+uoTUnqcvnAg1urORCIraYGAVkrT56RFLjRYtKh69INTVFPjgggradk2st0JAIsX12BoyI89e47na1hENIYlanc+NPy7N8E1ROPO4cM9+Otf9+Guu5YbroTpwpW00as7E4kMacGgHtJGHzt1wgmVKY0nFrvdiqoqF44dG0h+MRnW2+uF3W6BwxH+b1Wf83fwYA9OPLEqTyMjorEqbkhTSr2ay4EQFYpf/nIrlFK49daFKT92ZCUtlZCmb2arV9L0dicAPPXUVTh0qAd2e3bW7EyZUsxKWpb19nqjqmgA0NhYCgA4fJgraYkodUkXDojINgBqxM3dADYC+E+lVIcZAyPKh5aWPvz0p5txzTVzMXNmRcqPj96CI2h4dWd0u1ObvxRZSYt3Vmi6Jk92M6RlWayQNnVqCUQY0ogoPUbeQf4G4BloW3DcBOAv0ALaUXBTWxqjtmxpxec//yJeeukglAr/DPKNb7wJrzeA731vRVrPa7dbYLVKWgsH9J3pw+1OYwEvHaykZV+skGa3W1FXV8KQRkRpMbIFx4VKqVMjPt8mIpuVUqeKyM1mDYzILIFAENdd9zT27OnCL36xFWedVY+vfe102O0W/OpX2/B//+8yzJmT3twvEYHLZRvegsP4wgGHwwqPJzA8Pr3dmdYQDKmtdUedFUqZixXSAKChoYQbBxNRWoyENKuILFdKvQMAInIawpvZ+k0bGZFJ3nqrGXv2dOHBBy/B0FAA3/3uelx++ZMAtIn53/zmWRk9v8tlS7mSFiukmVlJq652ob/fh8FBH1wuu2mvM5H09npjbnY8eXIxDh3qycOIiGisMxLSPg3gQREpGf68F8A/iUgxgO+ZNjIik/z1r/tgt1twzTVzUVbmwD/+4wK89NJB7NjRgf/zf+ajpGR0NSQVbrcdg4P+lFZ3Op02DA1pP/PEmpOWbTU12iHtbW2DmDaNIS0benu9oYUCkWpr3di0iRvaElHqkoY0pdQGAAtFpHz48+6Iu58wa2BEZnnrrWaceurk0E7+DocNl102G5ddNjsrz69V0nwpHQvlcFgRCCj4/cHQnLTI1Z3ZVlPjAgC0tQ1g2rQy015nIunt9cVsd9bUuNDaOgClVNSed0REySTazPZjkZ8rpbpHBLRR1xAVOr8/iI0bj+L00+tMew29kpbKsVBOp/bzksfjj2h3ml9Ja2/nvLRsiTcnrbbWDb8/iK4uTx5GRURjWaJK2g9F5AiARO8U3wXw1+wOicg827e3Y2DAb2pIS3dOGgB4PIGczEkLV9IY0rJBKRU3pIVbywOorBw9Z42IKJ5EIe0YgB8nefyHWRwLkenWr28BAJxxhpmVNBv6+lJf3QloIS3RiQPZEhkcKHODg34EgypuJQ0AWlsHMHfupFwPjYjGsEQnDqzM4TiIcuLtt5tRXe3CzJnlpr2Gy2VDa+sAPB7jlTS93Tk05B91wLoZKiocsFqFlbQsiXVup45VSyJKl4k7MREVnvXrW3D66XWmTuB2u+2hdqdeIUsmdrvTvDGKCKqrXaykZUmikKZX0o4d4+bBRJQahjSaMLq6hvD++52mtjoBrZLW2+uFUjAc0sKVtEBOThwAtJYnFw5kR+JKWni7EyKiVDCk0YSxYcNRADB10QCghTR9JV8mlTQz252A1oZjcMiORCGtqMiK8nIHq5ZElDIjm9lCRM4CMCPyeqXUwyaNicgU69e3QARYvtzckOZ2hzemdTgMfYtFhDQ/bDYtnJnZ7gS0Cs+WLa2mvsZE0dennbsaK6QBWsuztZUhjYhSk/QdREQeATAbwBYAgeGbFQCGNBpT3n67BSeeOAnl5Q5TX8flCn9bpdPu1B9vdkirqHCgu5t7d2VDuJIW+/QGhjQiSoeRH/OXAThZKaXMHgyRWZRSWL++BR/72CzTX8vtDr9RF3K7s7zcwQ1WsyRRuxPQWst79nTlcERENB4YmZO2HcAUswdCZKb9+7vR3j6IM86oN/21Iitpxrfg0K7TtuDIzcKBigoHPJ4APB6/qa8zEeghLd65r6ykEVE6jFTSqgHsFJF3AIR+7FZKXWHaqIiy7O23tU1szV40AKRbSdOPhQrk5IB1AKG2b3e3B7W1xubOUWzJKmm1tdpK2mBQmV4hJaLxw8j/zN80exBEZlu/vgVutw0LFlSb/lrpzEnL9YkDQGRI86K2ttjU1xrvenq8cLlssNliVz9ra90IBhU6OwdRXe3O8eiIaKxK2k9RSr0K4AMApcO/3h++jWjMePvtZixbNiXum2g2ud2phzT9+CifL7dz0gBt/zjKTG+vF2VlsatoQHivNLY8iSgVSd+xROTjAN4BcB2AjwNYLyLXZvKiIjJJRF4UkQ+Hf6+Mc11ARLYM/3o6k9ekicvj8WPLlractDqB9Cpp+tw1rzeY0zlpgFZJo8z09MQ+XF0XeX4nEZFRRtqdXwdwmlKqFQBEpAbASwDWZPC6dwP4u1Lq+yJy9/Dnd8W4blAptSSD1yHCli1t8HoDeQppxuZ6RVbSct/u5ArPTGmVtPhbu9TWaud3MqQRUSqM/Khu0QPasA6Dj0vkSgC/Hf74twCuyvD5iOJ6661mADD9OChd5MKByNZnInolzecL5uSAdQAoL9cqPwxpmdMqabH3SAPClTSe8EBEqTAStp4TkedF5BYRuQXAMwCezfB1JyulWoY/PgpgcpzrnCKyUUTeFpGrEj2hiNw2fO3Gtra2DIdH48nrrzdhxowyTJ1ampPXi6yklZTEf+OOpFfSvN7cHLAOABUVTgDgXmlZkKySVlXlgkjhVdJefvkQfvnLLfkeBhHFkfTHfKXUnSJyDYCzh296QCn1ZLLHichLiL2/2tdHPL8SkXgb5U5XSh0RkVkAXhaRbUqpvXHG+QCABwBg2bJl3HiXAGib2L7xxhFccsmMnL1mZCUt3r5ZI1mtFlgsAp8vmLMD1ktLiyDCSlo29PR4UFpaFfd+q9WCqipXwYW0G2/8K1pbBxAIKHzhC6fkezhENIKhXoxS6o8A/pjKEyulLox3n4gcE5E6pVSLiNQBiHmAoFLqyPDv+0RkLYBTAMQMaUSxfPjhcbS2DuCcc6bm7DXTqaQBWjUtspJmdrvTYhGUlhYxpGVBstWdgNbyLKRD1tvbB0Kh8Z//eS0uvHA65s2blOdREVGkuD+qi8gbw7/3ikhPxK9eEenJ8HWfBvDJ4Y8/CeDPMV6/UkQcwx9XQ6vk7czwdWmCeeONIwCAFSsacvaakfPQ9DM5jSgqskbNSTO73QnwaKhs6elJHtJqalxobS2cOWmHD/cCAO67T/t5+t57N+dzOEQUQ9yQppQ6Z/j3UqVUWcSvUqVUWYav+30AF4nIhwAuHP4cIrJMRH41fM1JADaKyFYArwD4vlKKIY1S8vrrTaiqcuHEE3NXIYispIkYD1p2u2U4pOVmThqghbSeHm7BkQmfLwCPJ5BwCw6g8I6GamrqAwAsXToZV101B2vW7A612omoMBjZJ+0RI7elQinVoZS6QCl1glLqQqVU5/DtG5VSnx7+eJ1SaqFSavHw77/O5DVpYnrjjSM455ypKYWlTBk9r3Mkvd2pv1Hm4vigsrIi9PSwkpYJPeQmWjgAFGJI0yppDQ2l+NjHZuHYsQG8++6xPI+KiCIZmZk8P/ITEbEBWGrOcIiy5+jRfuzZ04UVK3I3Hw3QqmcrVjTg2mvnpvS4cLszNwsHAFbSsqGjQ2thVlU5E15XW+vG8eND8PkCuRhWUkeP9gMAJk924+KLZwDQVnsSUeGIO2FGRL4G4F8BuCLmoAkAL4ZXURIVsjfeaAIAnHNO7uaj6V5++eMptyvDCwdyNyetrKwIe/d2mf4641l7ux7SXAmv04+Gam8fRF1dienjSkY/JcFqtWDy5GLMmVOBdeua8z0sIoqQaE7a95RSpQB+OGI+WpVS6ms5HCNRWl57rQkulw2nnlqb89e22Swpt1j1SlpQy2ickzZGdHRoZ58aqaQBhbNXWk+PJ2qxw1ln1ePNN49AKc5LIyoURg5Y/9rwSsvlInKu/isXgyPKxEsvHcS55zbAbk9vjliujayk5WpOGrfgyEx7uxa6qqsTV9ImT9ZCmt5mzLeRK1LPPnsq2toGWVklKiBGFg58GsBrAJ4H8K3h379p7rCIMnP4cA/ef78zNNdmLMjHnLSyMgcGB/0FM09qLApX0hKHtIYG7cSLI0f6TB+TEVpICy92OPtsbe7mm28eydeQiGgEI+8CXwJwGoCDSqlV0DaU7TJzUESZevHFgwCAiy+enueRGKdvwZGrA9aB8CHrvb1seaarvX0QNpsl6RYc9fXaPDR9f7J8G9nuPOmkKrhcNrz3Ho/VIyoURkLakFJqCABExKGU+gDAPHOHRZSZF144gLq6YsyfX53voRhWVGTN6YkDAEJv0uOx5XnwYDe2bIl5mElWHTnSh/r64qRzEIuKrJg82R3a+iLfRrY7LRZBfX0JWloKox1LRMaOhWoSkQoATwF4UUSOAzho5qCIMhEIBPHiiwdx+eWzc7o/WqbsdguGhnK7ulOvpI23xQPBoMIZZzyGo0f70dt7h+EzVNNx+HAPpk0ztr93Y2NpgYW06L3d6uqKGdKICoiRhQNXK6W6lFLfBPDvAH4N4CqTx0WUtnffbUVn59CYmo8GjN7MNjdz0sZnJW337s7QBP233jJ3W4lDh3oNh7SGhtKCbXcCDGlEhSbhu4CIWEXkA/1zpdSrSqmnlVLj68duGlf+9rf9EAEuvHBavoeSkpELByzmZ7RxW0nbvft46GMz51gNDvpw8GAPZs8uN3R9Q0Np6DimfAoGVcxD4auqXDh+fChPoyKikRK+DSilAgB2icjYerejCe3pp/fg9NPrUFtbnO+hpCS8BQcraZnSQ5rLZcO2be2mvc6OHR0IBhUWLza2F19jYym6uz15X6jR3++DUqOPsqqocKCry8O90ogKhJF3gUoAO0Tk7yLytP7L7IERpaO5uQ8bNx7DFVfMyfdQUhaupOXyxIHxWUnbv78blZVOrFjRYGpI27pVq9ItXlxj6Hp9G458z0vTQ+LIFakVFQ54vQEMDvrzMSwiGsHIwoF/N30URFny1FMfAgAuv3xWnkeSOm0LjtxW0srLtTfpsXTI+q5dnaitdaOyMv4O/83N2orLBQuqcN99WxEIBE35+9y6tRXFxXbMmlVh6PqGBm0bjqamXpx0UlXWx2OU/vUe2e6sqND+Tru6PHC77TkfFxFFM7Jw4FUABwDYhz/eAGCzyeMiSplSCvfdtwVLltSOqa03dNoWHMGcnjjgdNpgs1nQ3T02KmlKKZx00oO44IInEl7X0tKP+voSLFxYg6Ehv2m76G/d2oaFC6sNf60aGwujkqZXTkeHNK2y2tXFeWlEhcDIiQO3AlgD4P7hm6ZC246DqKCsXXsYO3Z04I47ThlTW2/oIitpuWh1AoCIDJ/fOTYqab29XiilreBNpKWlD3V1xaFq1a5dxxNenw6lFLZubTM8Hw0Ib2ib78UD4ZAWPSdNr052dY2Nfw9E452R+v8XAJwNoAcAlFIfAsj9idVESfzsZ5tRXe3CjTeelO+hpCVcSVM5aXXqtPM7x0Yl7dix5IeTK6XQ0tKPurqS0KpLMypphw/3orvbY3g+GgA4HDZUVDhw7Fh+t7mI1+7UV/sypBEVBiPvBJ7ILTdExAaAS3+ooBw40I2nn96LW29dBKfTyFTLwqNX0vz+IGy23FUCx1Ilzcjh5B0dg/D5gqivL0FVlQtlZUXYsyf7lbSNG48CAJYsSe1n1smTiw2FTTPFa3eWlGjz0Pr6xkZoJxrvjIS0V0XkXwG4ROQiAH8A8Bdzh0WUmvvu2wIR4HOfW5zvoaQtcp+03FfSxkZIi6xAxTsUXt+Mta5OO6pp+vQyU9qL69Y1w+Gw4tRTUw1p7gIKadHtzuJiLaT19/tyPiYiGs3IO8HdANoAbAPwGQDPAvg3MwdFlIqBAR9+9attuPrqE9DYaGzn90KkH7CurUTMXSWtrKxozGzBEVlJi7fpanOzFsjq6rR98mpr3aa0F9eta8ayZVPgcKRWudVCWmG0O0duwcGQRlRYjKzuDCql/lcpdZ1S6trhj9nupILx6KPv4/jxIdxxx6n5HkpG7HYrAGBoKJDTSprW7hwbIS2yAhVvHl24kqZN0jejvdja2o933mnBypWNKT+2psaNjo78rp7s6fHC4bCiqMgadTtDGlFhifsjoIhsQ4K5Z0qpRaaMiCgFSin87GebsWRJLc45Z2q+h5ORoiItmA0N+XNeSRsr7c7ISlq8DVcj252AVrlqbc1uSPvDH3YjEFC44YZ5KT+2rKwo7ycO9PZ6R1XRAO2EBhGGNKJCkahO/7GcjYIoTWvXHsb27e349a8vGZPbbkTSK2mDg7kNaWO1kjYwEDtINDf3obzcEdqMdfJkN/r7fejv96K4eHQwScdjj72PhQursWCB8ZWdurIybVd/j8efcqs0W+KFNBGB221nSCMqEHF7Kkqpg/qv4ZtOGP64FUBnTkZHlMR///fG4W03Tsz3UDIWXUnL5cIBLTQMDRX+UUBHj/aHVu/Gr6T1hapogDYnDUDWqmkHD3Zj3brmtLd60VdU5jMYxwtpgNby7OtjSCMqBOlsZtsAbmZLBWD79jY888w+3HHHqXC5xv4RNvmrpI2do6GOHevHzJna4pB4lbS2tsFQMAOyv0HrX/6yFwBw3XVz03q8Ho7yGdL6+nwJQxoraUSFgZvZ0pj1wx9ugNttwxe+sCTfQ8kKvZI2OOiHzZbbShpQ+IesK6Vw7NgAZs7UNqgdGIhdSevp8YY2ZQUijzrKTkh79dUmTJ9ehjlzKtN6fLiSlr9Q3NvrDe2JNlJJCUMaUaHgZrY0Jh0+3IPHHvsAt966CJMmufI9nKzQV9rlupKmh4ZCXzzQ3e2BxxMIhbR47c7ubs+IkJbdStru3cexaFHqc9F0hRCKk7U7GdKICgM3s6Ux6Qc/eAcA8JWvLM3zSLInXyFNDzSFfjSUvmhg1qwKAPHbnVpICweQbB4arpTCvn1dmDWrPO3nGAtz0hjSiAoDN7OlMaepqRf/+7/b8KlPLcD06em/WRaa6JCWu3bnlCnaJPvm5j54PIW7eEDffiNRu1Mphe5uT9RO+tlsd7a1DaCvzxcKiunQ24z53IaDIY1obDDyTuAC8KC+mS2AB4dvI8o5ny+Az3/+JYgA//qvp+d7OFkVOSctl5W06dPLIKIdUD9p0s+xdu0hw499+eVDeOGFA+YNLoK+S384pI0OEgMDPgQCKqqSVlbmgEh2QtqhQ70AgGnTStN+Dn1rkHjtWrMppZIuHODqTqLCYCSk/R3RocwF4CVzhkMUX1vbAC666A/4y1/24gc/OBczZoyfKhqQv0qa02nDpEkubNhwFAMDfvzpTx8afuwFFzyBSy5ZY+LowpqbtZDW2FgKi0Xg8Yw+u1Nv2UbOSbNYBGVljqyEtM5OrWVaXZ3+z6kul76FSH6C0OCgH8GgirtwwO225y1AElE0I+8ETqVU6HTi4Y/dCa4nyrrt29uwfPnvsH79UTz88EfwpS+Nn7lounzNSQOAs86qD328d2+Xocf4/cGYH5vlyJFeOBxWTJrkRFGRNU5I04JYZEgDtJZnNuak6SFt0iRn2s/hdmshLd7qVLPpbVa2O4kKn5GQ1i8ioUMRRWQpgEHzhkQU7d13j+G88x6HxxPAa69dj098Yn6+h2QKvd3p9+f2gHUA+OUvL8Jzz12Da645AXv3dht6TF9feE5VLg4MP3KkDw0NpRAROBxWeL2jQ5o+GT92SMtGJU37ry+TFcX6nn75qlYlC2luty3uogwiyi0jZ5J8GcAfRKQZgACYAuB6MwdFpGtp6cOFF/4BxcV2rF17fUYTtgtd5GHXudwnDQDq60tQX1+CZ57ZhxdfPJj8AYie+H78+BCmTk1/npYRTU19mDpVOzTd4UitklZWlp2jr/RKWmWlI8mV8VksWsjMVxBKHtLs8PmC8PkCoQ2WiSg/koY0pdQGETkRgH6S8C6lFH/Mopy4667X0Nfnw5tv3jiuAxoQHdJyXUnT1dS40dPjNXSuZGRI6+0157+EgQEf7rtvC84/fxr27DmOiy+eAcBISIsOIKWl9qhzP9PV2TmE4mJ7xmdu5nPel74oIFG7E9DaseXlDGlE+WTkWKgvAChWSm1XSm0HUCIinzd/aDTRdX9YiUce2Yl/+ZdlOPHEqnwPx3TRIS23lTSdfpxSe3vyGQ2RKwDN2k7i//2/t3Dnna9i6dJH0NLSjwULqgEg6Zy0yC04AC2QZGOMnZ1DqKpKfz6azuWyFeyctPCcOf4sTpRvRt4JblVKdemfKKWOA7jVtBERAVBBYM/vFqChoXTcbbURjz4nDchnJU2ba9XWljykRVfSzAlpa9bsjto49uyzpwJA3Dlp8dqd2Qxp2Tjhwu225X1OWqLVnUD+FjYQUZiRmr1VREQppQBARKwAYv8IRpQlLa9OQ9+hcvzq9+ehuHhi/HMrhHanvrVEIVTS9uw5jj17unDvvedj8eJarFt3BGecUQcgcbtTZHSVKLshLfNKmtttL9g5aeF2JytpRPlmJKQ9B+BxEbl/+PPPDN9GZIqmpl7s+8NJqDixHR//+LzkDxgnHI78tzvDR0QlXwlpdiXtuef2AwA+8pFZmD27AitWNITuix/StJ30LZbokFtaWoS+Ph+CQTXqvlR0dg5i/vzqtB+vGwvtTm7DQZR/RkLaXQBuA/C54c9fBPC/po2IJjSlFD7zmRegAhbM/dR7EMlPRSkfCqGSpoe0np7UQtrQUPYDx3PPHcCcORWYPbti1H1aSBv9mj09nlGtTiAcSPr74++0b0RHR7YqaflvdyZa3Qmw3UlUCJL+uK6UCiqlfqmUunb4WKidAO41f2g0ET3yyE48++x+zLzmA7hqM1+NN5YUwsKBVA5bj9wnLduBIxAI4tVXD4dWc45UVGSF1zt6A93ubm/CkJZJxU8pNS7anX19PjidtrjbvLDdSVQ4DL0TiMgpIvJfInIAwLcBfGDqqGhC2rTpKD7/+ZdwzjlTMfXC/fkeTs5ZrQK9cGiz5aeSpocZo+1OEW3cQ0OjW4+Z2LGjA319vqiTECIlmpNmVkjr6/PB7w9mJaS5XPmtpMVbNACw3UlUSOKGNBGZKyL/ISIfQKucHQYgSqlVSilW0ihrlFK4//6tOPvs1aisdOLxxy+H5KeQlFciEqqm5auSZrNZUFxsN9ju9KG42G7Knl/btrUBAE45pTbm/fHanVpIG93Gy0ZIy8ZpAzptTlr+Fg4kavmG250MaUT5luid4AMA5wP4mFLqnOFglt0fl2nCW7++Beed9zg++9kXcd55jdi06WbU15fke1h5o4e0yO04cq2srMhwu7O0tAguly3pnLQjR3pTOt/z8OFeAMD06WUx73c4bHEraSP3SAO0zWyBTENa5ud26pzO2OPPheQhLb9nixJRWKKFA/8A4AYAr4jIcwB+D+1YKKKUKaXwwgsHcP/972HXrs7QsTMHDvSgutqFX/7yItx666KMVt6NB3pIy+dxPOXlDsPtzpKSIng8/oQhratrCDNn/i++8pWl+MEPzjM0hsOHe1FZ6Yy7/UpRkSXuPmlmtTuzG9KsWW8RG5UspOlz0tjuJMq/uD+uK6WeUkrdAOBEAK9AO8OzVkR+ISIXZ/KiInKdiOwQkaCILEtw3aUisktE9ojI3Zm8JuXH0JAfDzywFSef/BtceukfsW7dEZx44iQsXToZZ55Zjx//eCX27bsVn/nM4gkf0IBwBc1uz18lrbzc2DmXfX3aSkmnM/H8qi1bWuHzBfHzn79reAyHD/eisTH+WaDxK2leE9ud2Qtp8dq1uaB93eLPSdMPgGe7kyj/jJzd2Q/gMQCPiUglgOugbcvxQgavux1ape7+eBcMb5r7PwAuAtAEYIOIPK2U2pnB65rK5wugubkP9fUlWa2EKKWwdWsbgkGFyZPdmDKlOC9zlpRSCAaVodc+fnwIDzywFffcsxlHj/bj1FMn45FHPoqPf3xe1CpGihZud+bv70hrdxqrpJWW2hEIBBNWhfQzMwMBZXgMR46ED1OPJdbCAY/HD683kLCSlskh69lud+azkjZjRuw2MqAdAO905m8fNyIKS+mU4OEjoR4Y/pU2pdT7AJLtgbUcwB6l1L7ha38P4EpoW4Dk1WuvHUZTUx/a2wfR1jaAlpZ+vPtuK7Zvb4fXG0BxsR2XXTYL118/D5dfPjutwKaUwvbt7fjDH3bh0Uffx7593aH7bDYLGhpKMG1aGaZNK8OiRdW47rp5mDGjPMEzas+5Y0c7XnjhIFpbB+B0WjF5cjFmzSrHCSdUwuMJ4MMPj2PXrk7s2nUcu3cfx5EjvQgEFIaG/Dh+XNvNfdWqRlx66UwsXFiDqionyssdwxUVK958sxm/+91O/OEPuzE05MdFF03Ho49ehlWrGifUnmfpCrc781tJ0+eEJdLb60V9fQmGhgIJ25169crnMz4n7fjxIZx00qS498cKafGOhAIKs5IWDCr4/cG4W2GYRW9TJ1JcbHyLkE2bjuLdd1tx5ZVzIAK89loTNm48hlmzyrFs2RScdNKkjA+kT5dSCj/72WZ85ztvo7bWjSuumIMrr5yDJUtqMhpTV9cQmpv70NExhI6OQQSDCqecUouZMyuyN3gipBjScmwqtBWluiYABXGI4x13vIytW7XVZxaLoLrahUWLavClL52K2bMrsGVLK/70pw/xxBO7MHVqCW69dREqK53o7fWirKwI8+dXY8oUNwYH/ThypA+vvHIYr7/ehJaWflRUOFBSYse+fd1obx+ECHDBBdPx9a+fgaoqF5qb+3D4cC8OHerBoUO9eO21w/jd73biq199DeefPw1nn12PY8cGcPhwL1pa+uHzaW9kRUVWtLcPht58tX2m4v8kX1vrxrx5k3DmmfWw2SwoKrJi0iQnBgZ8eP75A/jyl1+J+9iysiJ86lPzcdtti7FkSezVeRRbobQ7jVTS9HZnX583YbtTD0bBoPFKWk9P7P3OdLH+/eqLHWI9rrjYDpHMV3c6nbZQOzAT+ukSHo8fNltujz1LNicN0BYPGJ2TdtttL2Lz5mO49dZwc0UEUCr8cX19CaZPL4v6NWmSEyUlRSgpsaO83IGKCgccDm3Va3v7INrbBxEIKBQX29Hb68X+/d3Yv78bR4/2o61tAG1tg/B4ArDZtMqf9rWxoaTEjpNPrsI550zFQw/twJo1u0M/JP7Xf72D731vPaxWwZw5lZg/vwozZpShsbEMDQ0lmDKlGC0t/di5swM2mwVnnFGHWbPKEQwC+/d345VXDuH55w9g8+ZjoT9fpNNPr8PKlY2oqXEhEFDYu7cLu3Z14tAhbeGMUlpwBLS/H6tVUFxsR3GxHQ6HFUVFVjgc1tDH+v8Dg4N+DAz4MTjoh8fjh8cTCP0aGtLuCwYVSkuLUFZWFFrQoz+X9ssGp9OKQEChv9+Hvj4vhoYCENG+N7WxGfqS51wqP9+ndm3yix0OK5577lrjT5plpoU0EXkJwJQYd31dKfVnE17vNmgnI2DatGnZfvoojz56Gex2C6qrXaiocMacS/Xzn1+AZ5/dj3vv3YxvfnNdwudzOm0444w6LFpUg+5uD3p7vbjyyjk4/fQ6XHbZrKSrHQ8c6MbvfrcTDz64HS+/fAjV1S40NpaioaEEDocVSgFebwBz5lTg3//9THzkIzPR0FAKny+A1tYB7N3bhQ8/7ILDYcUJJ1Ri7txKVFYmrhYcPKj9h9nRMYSeHg96e33o7fVi4cJqXHTR9Ky8kU1EetU1vyHNeLuzpMQOp9OG48eHEl6nM3Isk1JqeJVm/CDhcGhvNoFAMNR+T1RJExGUlGR2fmdn5xCqqjKvogHa9zwADA0FUFyclac0RCkVCteJaJvtJm93KqWwZ89xrFzZiPPPnwan04qzz56KpUsn48CBHrz7bivef78DBw/24NChHqxf34I1a3anVFWNVFJiR319CWpq3JgzpwIulw0+XxAeT2A4yPjQ2jqAZ5/dj+9//x1YrYIf/OBc3HnnaRARdHYO4qWXDmLbtnZs3679evbZ/aMqwZEhM5LVKjjzzHp861tnY86cClRXu1BV5YJSCn//+yGsWbMb//3fG0MrmauqXJg7txJnnlmPoiILRLS9EPVwEAgE0d/vQ3+/Dx5PAF5vEMePD8HrDcLj8Yf+nlwuW+hXZaUTTmc4gOnhVP8hRP81NBRAX58WeCNDnYj291hSos0nVUr7nrRYJKWAkyupBEeVwsWRgTmRVH64NINpIU0pdWGGT3EEQGPE5w3Dt8V7vVAbdtmyZab+rRo5u89qteDyy2fj8stno7W1HxaLoLS0CB0dQ/jgg060tQ3A5bKhutqFU0+dHPpPOx0zZpTj3/7tTHz962cgEFCG2yd2uxVTp5Zi6tRSnHtuY/IHRJg+vRzTpydur1Lq9OOg8jknraLCiYEBP3y+QMJWvV6RcblsaGmJX5WNPIi9r88bc4uMSENDfvj9wYSVtHAlKgC3OzqkxQt3mR6ynq3TBoDo8Y/k9wfR2+tN+oNSOgYH9YpL4h+ijLY7OzoG0dOj/VD55S8vjbpv3rxJmDdvdMs6EAji6NF+dHV50Nen/XDX3e1BV5cHXm8AbrctFH5sNgv6+31wu22YObMcVVUuQ9WP7m4PtmxpxYwZZVH/T02a5MLHP34iPv7x8LX6SRJNTb1obu5DXV0J5s2rxOCgH5s3t+Lw4R5YrRbU1RVj+fK6uP8uly6dgq9+dTkCgWAo4GZyBBkRUNjtzg0AThCRmdDC2Q0A/r/8Dik9tbXhH5Xr60tM2wdMRPK2Uz1lh/7+k+92J6C90VVXu2NeEwyq0DmYyVZ3RgajePuYRdLblskqaYAe0uyh544c/0iZhrTjx4dQUZHtStrov7cvfvHvePjhHejtvSPrC4SSndupM9ru7OjQKqi1tbH/ncRitVpCPxyapbzcgfPOM/aDp4igqkoLhYsXh6dnuFx2XHjh9JRf22q1MJxR1uTlnUBErhaRJgBnAnhGRJ4fvr1eRJ4FAKWUH8AXATwP4H0ATyilduRjvES5olcJ8rlPWkVF8vM79Tdwrd1pNbRwADC2Qap+2kGyOWlAdCVKX7kZawsOQNvQNpOQ1tvrSxgcU5Goknb//VtD81WzzXhIM9buDAdjhhIiM+SlkqaUehLAkzFubwbw0YjPnwXwbA6HRlQQ8nnigB6OurqSzzPT251GK2lGjo9KVhEDwiEncvGAsUpa+nt/9fZ6ccIJFWk/PlKiSpqutXUA06bF3yojHfrXwsjqzpaW5CGxq0v7O89WhZGIok3AExKJCle43VnYlbS+vvCbvbbnV/ZCml4RM9ru1IXnpMUOaWVljgwPWE++KtKoRJU0XWvrQFZeK1K22516kGcljcgcDGlEBURfbKS/iedD5Jy0ePTFAEYrafpRQ0Ymo6dSSYvctb+724PiYnvchTOZzkkzsr+YUU6nNv6R4TZyJVmiFbPpivy6JWK83an9fbKSRmQOhjSiAqIvuddDTT7oVRGj7U6n04ZAQMU9QL231xeaWJ69SprWLoyupCXeWy2TOWnhhRLZ+brEGj8Q3b41I6SFv27ZWd0ZrqQlXgxCROlhSCMqIHrQyWdI06sixtqddrhcWuCIF8B6e72YPNl4SDNSSdMrUSPbnYmCXSaVtIEBH5RKPpfLqHAlbeTRVpEhLfledalKpd1ptJJmsQhKSrgvIpEZGNKICoheSXO787c7jh50ErU79Qn40SEtduWlt9ebUiVNf91EQUJvd0aGnO5uT5JKWhE8nkDoFI5UGG0TGhWrXQtEV9L0Y6iyyejCAbfbDr8/mPBUEkD7mrvdNh75RmQShjSiAhIIaCEtWxWbdOj7POkr92LRK2l6uxMYXRUC9B3uwyHNSAutp8cLt9uWcPFEuF0YDjnJjpLK5PzOyMphNsT7O4uupJnZ7kxeSQOSf72GhgJ5O5eTaCJgSCMqIHfccSoAYOpUczY8NirZ+Z16ZamkpChhu1NvE6ZaSUu24W2sdqFWSUvc7gTSC2lGw41R8VZ35iKkOZ22pKeShBd6JP56DQ35Q18LIso+hjSiAvKFL5yCtrbPo64u3yEt8fmd4bZZ4jlpels01YUDySaix9uCw0glTV+YkIrIUJoN8fZJi144kP05adq5ncmrgfopDsm24RgaCmR0pB0RJcaQRlRALBaJexRTLlVUOJO2Ox0OK+x2a9ztJIBwmNPOXDR24kCyihgQO+QYDWmZVNKy1e7MZyXNSDXQaLvT42EljchMDGlENEqySppWkdHe7BNX0sJtQrfbHndxQaSenuSHsI8MOT5fAAMDftPnpGW73Tky2Op/nilTivMa0oy3O1lJIzITQxoRjVJR4Uza7tSrSi6X9nuykOZyGd3WIfVKmpG91TKrpIVXs2aDzWaBxSJx90kzM6QZadkab3f687rxMtF4x5BGRKOUlydb3Tm6kpao3VlaaofbnfhkAl06lTQje6vpc7EKoZImInA4Rh9MH66kuTEw4B+1RUemtEqakTlpxld3spJGZB6GNCIaRV/dqZSKeX9X11AoSOlzkoy1O41W0lILaXolzVi7M/VD1o3uL5YKp9MWY06a9vczZUoxgOwvHogM14lwdSdRYWBII6JRKiud8PmCcdtdx497UFWlnUyQaE7ayDM+k1VmgkGF3l5vwrYloO3lZrNZQpUoY5W0TOak+VBUZEVRUfYCiVZJG9nu1PbJmzxZD2nZbXkaXzhg7KxVj4eVNCIzMaQR0SgNDaUAgKam3pj3d3YOorIyOqQlbncWGTq0W7/eyFmQDoc1VInSW7OJ5rI5HDbY7Za0253ZPvrI6bSOamdGLhwA8hnStK+pkTlpDGlE5mFII6JRpk3TQtqhQ/FC2hAmTdJCmv4mnajdWVxsHz4PMvGbvpGKmM7ptIWCoX7Qtx4c40n3/E6j4SYVDoctxokDI9ud2Qtp2ukPPkMt23C708iJA2x3EpmFIY2IRmls1ENaz6j7fL4A+vp8oZCWbAsOt9sGq9ViqJLW06OFtGTtTiB2Ja2iInG4SzekaeHGjEpa7HbnlCnaXnnZnJM2MOBDMKgMLRzQgzfnpBHlF0MaEY1SX18Cu92CXbs6R92nV3f0kGa1WmC3W2Ke3RlZgYq3unPr1la0tQ0AALq7jbc7nc7wnC59TMkeV1iVtFirO82rpKVySLyIwO22sd1JlGcMaUQ0it1uxZln1uPHP96Eq69+KhSiAKCtbRAAMGmSK3Sb0xk7gEWGm1gLB/bt68IppzyMyy9/EkA6lTS93elBWVkRrNbE/6VlVknLbkiLtbpTr6Tpx2h1dmYvpKV6/mhxsZ0LB4jyjCGNiGK6555VuOGGE/HUU3vwH//xZuj2I0f6AEQfAu9yJQ9psdqda9bshlLA+vUt6OwcTKmS5nDYotqdyVqdwFiopGl/nuJiO0pLi7JaSUs1pCVrTwcU4PMFOSeNyEQMaUQU0ymnTMajj16Gm28+GY8++j78fq3Kc+SItpigoSE6pMVa3dnTE93uHFmZ+dOfPgx9/N57bSlV0iLbnVpIS7xoANA2tC2c1Z3x90mz262orHTkOaQlbnd6gtrbB+ekEZmHIY2IErryytno6fHinXdaAAAHD2qLCerrwyHNSLvT7bbD4wkgGNQ2yD1+fAgbNhzFpz+9EACwc2dHSqs7I9udx48PobLSaCUtnc1ss79wIFYlzesNwm7XjoyqrHRmdeFAqhvyapW0+H9XQ0EtnLHdSWQehjQiSuj886dBBPj73w8BAN5+uwULFlTD4Qi/OSdqd+pVMX3vLf2Q9ZdfPoRgUOETnzgZRUVWHDzYg54eLywWCW0BkYi2BUdqlbSyMkfalbRstzvjVdL09qEW0rJZSdMXDhgLm9qctPjtznAljSGNyCwMaUSU0KRJLixeXIu1aw/D7w9i3bpmrFgxNeqaeO3OkQsHgPC2DuvWHYHTacOZZ9Zj2rRSHDzYg64u7UgoEUk6rugtOIZSmpMW77irWHy+ADyeQNYXDsQ6ccDjCYRONZg0KbshLdXzR5PtazfEdieR6RjSiCiplSsbsW5dM955pwW9vV6ce25j1P1GFw4A4Q1St2/vwMknV8Fut2L69DIcPNhjuG0JRG9m29lpNKTZEQyqpKsWI6WydUUq4u2TFl1Jy367M5WFA4nmpOntzsiKKhFlF0MaESW1cmUjhob8+NGPNgLAqEqa02kdFdICgSAGBvwxQpp23Y4d7ViwoAoAMGNGOQ4c0ENa8rYlEK6k9fd70dfnC+0tlkg6h6zrFahczEmLbHdWVTnR3j4YmsOXqVTnpCVrd7KSRmQ+hjQiSurccxsgAjz55IeYNascU6eWRt0fq5KmnwKgh67wnDQ/eno8OHKkDyefrIe0Mhw92o/m5v6UQ9rRo9oebnV1qYQ04/PSUq1AGRW5hYgust05c2Y5vN5AaDVtpnp7vXA6bbDZjP23n3zhAOekEZmNIY2IkqqsdGLx4loAwOWXzx51v8tlH1UV6ujQNr2tqtJDWrjdqR/cPmNGOQBg9uwKAMCWLa0ptztbWrR921KrpBkPaXq704xjofz+IAKBYOi2yHbnnDmVAIC9e7uy8nq9vT7DiwYAA1twKH11JytpRGZhSCMiQ37yk5W47rq5uOuu5aPui9XubG/XQ5p2MkHkwoHm5n4A4eqXHtKA5Iek68KVNP25SpI8otAqaVq4iaymRbY758ypAAB88MHoo7nSkeoKVb3dGW+RBStpROZjSCMiQ1aunIYnnrgiZhiK1e7s6NBWJoYraXpI84WqX/pea7NmlYceZ6QiBkRW0vqHH+dO+pjMKmnZ34IDQFQFMrLdOX16GRoaSvHMM/uy8nqpnprgdmuLLLze0WeyAuGQxhMHiMzDkEZEGdO24Ih+Mw+3O7VKmt7uHBz0o7lZC2l6JU2/BgAaG6Pnu8XjcFgRCCg0NfXCahVUV5sV0nJXSYtsd4oIbrllPv7613147722jF8v9ZCmhch4LU9uZktkPoY0IsqYXtWKbI3pbUj9sPBwJU1rd5aVFaG4WAsNkfui6fPUkr+mFhIOHuzB5MnFsFiS762WSbvTjGOhgJGVNH9UZeqf/3kZyssd+OY312X8eqmGNH1D4XgrPLm6k8h8DGlElDF9vllk4DhwoAdVVa5QmzBy4UBLS1/UsVIA8IUvLEFpaRFOO22KodfU9+c6cKDH0MpOILN2pxknDgCIqkBGtjsBbX7epz61AM88sy90BFa6Uj3aauS+diPxxAEi8zGkEVHGwiEtHDj27evCjBllo64ZGPChubl/VEj76U/Px759nzZ0bicQruDs399teB6bHlLSqaTpoSVb9PFHzuXzegOj5nidfXY9vN4Atm1rz+j1Ul04EFn5jIXtTiLzMaQRUcb0AKYHjvfea8O6dc1YunRy6Br9zVxrd/aNqn5ZrRZD88p0euBobR0wXEmzWi1wu20pz0krKbEbaqemYuTfGTC6kgYAJ52k7SW3Z09XRq+Xbrsz/pw0LhwgMhtDGhFlTK8y9ff7sG9fF1atehwVFQ7cfXd4uw6LReBy2YbbnaMraamKPAbKaCUN0M/vNH7igNYmzG6rE4gf0kaGnpkztTl6+/Z1pf1aSin09flSXt0JxG93MqQRmY91aiLK2KRJ2jYbnZ1D+OlPN8HjCWD9+pswc2ZF1HVutx3Nzf0YGvJnIaSF91MzWkkDwoesG6VX0rIt1jy+WO1Ot9uOKVOKM9rUdmDAh2BQZb3dabdbYLXyZ30is/C7i4gypoe0pqZePPnkHnzqUwtCO+ZHcrtt2LPnOIDUglUsmVXSUpuTlu1FA4DxdiegbUty5Ehf2q+lVw7TqaTFa3d6ghbORyMyGUMaEWVMD2lr1uzG0JAfV1wx+ugoQAsme/d2A0Ce252pre40p90Z3jdOF6vdCQCTJ7tx7NhA2q+VziHx4S044rQ7lYXbbxCZjCGNiDKmh7Tf//4DlJTYce65DTGvc7vtoU1uM62kRR4fZXQDXCC9dmcqZ14aNbKSppSK2e4EtBB67Fh/2q+VztFWydudltA2KERkDoY0IspYZGC65JIZcd+89Td+IPOQFtlqM3Jupy71dmduFg74/UEohZjtzsmTi9HaOoBgMPY5msmkF9KSre60ctEAkcn4YxARZSxy8vjHPha71QkgtAdaebkjdNpAJj75yfk4dqwfNpvxnzfTq6SZH9L046HitTsDAYWOjkHU1BjfpkSXTkhzOKywWCThZrZsdxKZiyGNiLLijDPq8PbbLfjoR2fGvWbyZC1gZFpF0/3mN5em/JiystQXDpixutPhsEIEGBzUQpB+kHmsKqQ+5+7o0f40Q1rqCwdEBG63LeEWHFw4QGQufocRUVY8/PBHsWtXJ2pr4wcwPWxMnZrZogFd5JmfRpWWFmFgwI9AIJh0+4h09hczSkTgdNpGVdKKikaPST//tK1tMK3XSvf8UbfbHndOmidoYbuTyGQMaUSUFSecUIkTThi97UakyZO1kKa3+vJBD1x9fb6kR1ANDvoRDCpTKmmA9vdgpN2ph7TW1vRWeOqrO1MNm8XF9oRbcDCkEZkrLwsHROQ6EdkhIkERWZbgugMisk1EtojIxlyOkYiyb/HiGgDAiSdOytsYUjlkPd1wY1RkSEvU7sw0pIUraan9ORK1Oz2KCweIzJavH2e3A/gHAPcbuHaVUiqzk4WJqCCsXNmIF1+8DmeeWZe3MaQS0vS5XGas7gRiV9JitTsrK52wWiWjkOZy2VJaYAEkbndyThqR+fLyHaaUeh9Ibz4JEY1dIoILL5ye1zHoe54VWiUt3O4c/d+yxSKornahrS3dkJbevDq2O4nyq9D3SVMAXhCRTSJyW6ILReQ2EdkoIhvb2tpyNDwiGmtSa3fqlTTz56SF252xg09trTujSlo6IS1hu5Mhjch0plXSROQlAFNi3PV1pdSfDT7NOUqpIyJSC+BFEflAKfVarAuVUg8AeAAAli1blt6Oj0Q07qXW7kxvLpdRRtudQGYhLd1D4hOu7lQ8cYDIbKZ9hymlLszCcxwZ/r1VRJ4EsBxAzJBGRGREegsHzKukdXQMAUjc7gS0kLZhw9G0XifdSlqiducQN7MlMl3BtjtFpFhESvWPAVwMbcEBEVHa9LDS01MIlTS74XZnTU1+2p3x56RxdSeR2fK1BcfVItIE4EwAz4jI88O314vIs8OXTQbwhohsBfAOgGeUUs/lY7xENH4U6pw0I+3Onh4vPJ7Y7cdEMlk4kHhOGtudRGbK1+rOJwE8GeP2ZgAfHf54H4DFOR4aEY1zLpcNFoukNCct36s7gehTBxoaSlN6nXQrafFOZ/ArQRDCShqRyQq23UlEZAYRMXzIek+PFw6H1bSKUaqrO4H0NrRNd+GA3uYd2fIcCmpvHZyTRmQuhjQimnCMhzSPaVU0IF67M96cNBeA1ENaJuePxmsNe4ZDGtudROZiSCOiCaey0oHjxz1Jr+vt9aGszLyQVlJix9CQH35/MOHZnUBkuzO1kDYw4EMwqNIKaXr1LX5IYyWNyEwMaUQ04dTXl+DIkd6k1/X0eFBWlvgQ9kyED3v3htqdiRYOAKlX0vSjrTKppOkLKHQMaUS5wZBGRBPO1KklOHKkL+l1PT3pTbg3KjIEJVs4UFpaBIfDmkZIS3/xQ7xK2lBQC2c8u5PIXAxpRDThTJ1agqNH+3HDDX/BP/5j/J19enq8Jrc7w3O+PJ4ALBaJewi6iKS1V1omIS1uJU2xkkaUCwxpRDThzJhRDqWAxx/fhd/8ZnvceV69veaGtMjD3r3eQNxFA7p0jobST03IZHUn56QR5QdDGhFNOCtXNsJikdDn773XFvM6s1d3Rq6e9HgCSUNPba0LbW2DKb1GdippI9udDGlEucCQRkQTzqxZFVi//ibs2HELAOD99ztiXperdqc+Jy15SEun3Zn+woFkqzs5J43IXPwOI6IJadmyKQgGFWw2S8xFBH5/EIOD/pys7ky13amUgogkvFaXSSWtuFgLaVzdSZQfrKQR0YRlsQjq6orR3Dw6pJl9JBQQXakyUkmrqXFjcNAf99DzWDL5c1itFrjdttGVNKWNkyGNyFwMaUQ0oU2dWoLm5v5Rt/f0aJvdmrtwIDzny+PxG2p3AqltaKvPJ9OrYumMcWQlbZDtTqKcYEgjogmtrq4kZiWto2MIAFBV5TLttd1uO0T0dmfQULsTAI4dMx7Senu9cLlscbf2SKakZPQRWgMBbZzpBj8iMoYhjYgmtOpqFzo6Rq+YbG8fDN1vFotFUFxsR2+vsYUDdXXFABAzVMbT25veuZ06rZI2IqQFGdKIcoEhjYgmtEmTnOjsHIJSKur2XIQ0IHzY+9CQH05n4pDW0FAKAGhqSiWkZXZqQkmJfVQlrX+4kuZ2s91JZCaGNCKa0KqqnPD5gqPmXenzvnIR0rQ5aYGkc7yqq11wOKxoakp+7qgu05Cmh8hIAwErbBKE3c6FA0RmYkgjoglNn3M2suXZ3j4Ii0VQWek09fXLyorQ3e3B0FDyhQMigoaG0pRCWl9fZiGtqsoVmp+nGwha4bYE0n5OIjKGIY2IJrRJk7QQ1tkZHUTa2wdRVeWMOpnADBUVTnR3ezE0lLySBgANDSU4fNh4SNMOiU9/7lh1tSvU+tX1B2wotjKkEZmNIY2IJrRElTSzW50AUFHhQFfXkKE5aQBSrqR1d3tQXp7+hrzV1a7hfdz8odtYSSPKDYY0IprQqqq0StrIlp4W0tymv74W0jyG5qQBWkg7cqQPwaBKei2gH22VWUgDov9+BgJWuFlJIzIdQxoRTWh6Ja2zM7qS1tY2kJNKWnm5FtKMzEkDtHanzxc0vKGtVklLf06a/ncQ+XoDQQsraUQ5wJBGRBOavjAgdiUtN+3OwUE/enq8hippjY1lAGBoXprH44fHE8io3VlTo/0dRM5L45w0otxgSCOiCa2oyIrS0qKoOWlKqZyGNAAIBpXBOWklAGBoXlpPj7Z1RiZHW9XUjD7loC9gRbHVH+8hRJQlDGlENOGNXMHY3e1BIKByFNLCW3wYnZMGGAtp3d3a+aOZVNKmTdNeb//+7tBtXX47KmwMaURmY0gjoglvZEjTP9ZbfWbSK2kADM1Jq6lxw263GDp1IBshrbi4CHV1xdi7tyv8vAEbKmy++A8ioqxgSCOiCa+mxoW2ttEhLZftTsBYJc1i0Ta0PXiwJ+m12Wh3AsDs2RWhkBYMKnT77aiwMqQRmY0hjYgmvJoad9TqRT2w5TqkGT0ZYNas8qj2YzzHj2uLITKppAHA3LmV+OCDTgDaMVMKwkoaUQ4wpBHRhDe63ZmbczuB6DlpxkNaBfbt60p6nR42a2sz2+9t/vxqtLYOoL19IBT8OCeNyHwMaUQ04dXUuDA46Ed/v9YeDLc7c7OZrc5oSJs9uxxtbYOjDj4fKVuHxM+fXwUA2LGjAy0t/QCAKUWejJ6TiJJjSCOiCU/fZkKvPLW3D6KoyIqSkvTPvDSquDj8Gkbnjs2aVQEASatpra0DKC93wOFIPtctkQULqgEAO3a048gRbVXpVMdgoocQURYwpBHRhKdXmvQKWnv7IGpqXBAx93B1AFGvkcqcNADYty/xvLS2tsGsrFCtry9BebkD27e34/XXj8BpCWCOy9iJB0SUvsx+vCIiGgfClTQteLS25uZIqJEmTzbWXp09uwIAsHt3Z8LrWlsHQn+2TIgI5s+vwrZt7ThwoAeXVLbxxAGiHGAljYgmPL3apLc7W1r6UV9fkrPXt1q1app+RFUyFRVOTJ9ehi1b2hJed/BgD6ZPL8t4fIDW8nzjjSNoaurFNTUtWXlOIkqMIY2IJryR7c7m5r6chrQXXrgODz/8kZTaq0uW1GLLlta49wcCQRw82IOZM8uzMUSsWNEQ+vjyqmNZeU4iSoztTiKa8MrLHbDbLWhtHYDfH8SxYwOory/O2euff/60lB+zZEkNnn56D/r7vSguHj2XbefODvj9QcybV5mNIeLaa+fiL3/Zi7PPnoqKP/01K89JRIkxpBHRhCcimDatDPv3d6O1dQDBoMppJS0dS5dOgVLApk3HcO65jaPu/9vf9gMALrpoRlZez+m04fHHL9c++VNWnpKIkmC7k4gIwLx5ldi1qxN79hwHAMyYkZ02oVnOOqseAPD660di3v/ss/uwZEltwYdNIoqPIY2ICMC8eZOwe/dxbN/eDgA46aRJeR5RYlVVLsyfX4U33mgadV9X1xDeeOMIPvrRmXkYGRFlC0MaERG0ifGDg3584Qt/R1lZERobs7Mq0kwrVjRg3bpmBALBqNuffXY/AgGFj3yEIY1oLGNIIyICcNlls3DddXMBAP/wDyfAYjF/I9tMrVjRgJ4eL957L3orjl//ehumTSvFWWdNzdPIiCgbuHCAiAhAUZEVTzxxBQ4c6MaUKblb2ZmJFSu0EPbSSwcxc2Y5XnrpIJ58cg9efvkQfvSj88ZE0CSi+BjSiIgiFPqCgUiNjWU455yp+OpXX8PXvvY6AgGF8nIHPve5xbjjjlPzPTwiyhBDGhHRGParX12Cb31rHWbNqsBHPjITp59eB5uNM1mIxgOGNCKiMWzevEl47LGP5XsYRGSCvPy4JSI/FJEPROQ9EXlSRCriXHepiOwSkT0icneOh0lERESUN/mqib8IYIFSahGA3QC+NvICEbEC+B8AHwFwMoAbReTknI6SiIiIKE/yEtKUUi8opfzDn74NoCHGZcsB7FFK7VNKeQH8HsCVuRojERERUT4VwuzSfwTwtxi3TwVwOOLzpuHbiIiIiMY90xYOiMhLAKbEuOvrSqk/D1/zdQB+AI9m4fVuA3AbAEybNi3TpyMiIiLKK9NCmlLqwkT3i8gtAD4G4AKllIpxyREAjRGfNwzfFu/1HgDwAAAsW7Ys1vMRERERjRn5Wt15KYCvArhCKTUQ57INAE4QkZkiUgTgBgBP52qMRERERPmUrzlpPwdQCuBFEdkiIr8EABGpF5FnAWB4YcEXATwP4H0ATyilduRpvEREREQ5lZfNbJVSc+Lc3gzgoxGfPwvg2VyNi4iIiKhQFMLqTiIiIiIagSGNiIiIqAAxpBEREREVIIY0IiIiogIksbcoG9tEpA3AQRNfohpAu4nPT/nHr/H4xq/v+Mev8fg23r6+05VSNSNvHJchzWwislEptSzf4yDz8Gs8vvHrO/7xazy+TZSvL9udRERERAWIIY2IiIioADGkpeeBfA+ATMev8fjGr+/4x6/x+DYhvr6ck0ZERERUgFhJIyIiIipADGkJiMilIrJLRPaIyN0x7neIyOPD968XkRl5GCalycDX9xYRaRORLcO/Pp2PcVJ6RORBEWkVke1x7hcR+dnw1/89ETk112OkzBj4Gq8Uke6I7+Fv5HqMlD4RaRSRV0Rkp4jsEJEvxbhmXH8fM6TFISJWAP8D4CMATgZwo4icPOKyfwJwfPjA+J8A+EFuR0npMvj1BYDHlVJLhn/9KqeDpEw9BODSBPd/BMAJw79uA/CLHIyJsushJP4aA8DrEd/D387BmCh7/AD+r1LqZABnAPhCjP+nx/X3MUNafMsB7FFK7VNKeQH8HsCVI665EsBvhz9eA+ACEZEcjpHSZ+TrS2OYUuo1AJ0JLrkSwMNK8zaAChGpy83oKBsMfI1pDFNKtSilNg9/3AvgfQBTR1w2rr+PGdLimwrgcMTnTRj9jyN0jVLKD6AbQFVORkeZMvL1BYBrhkvoa0SkMTdDoxwx+m+AxrYzRWSriPxNRObnezCUnuHpRKcAWD/irnH9fcyQRhTfXwDMUEotAvAiwlVTIhobNkM7bmcxgHsBPJXf4VA6RKQEwB8BfFkp1ZPv8eQSQ1p8RwBEVk4ahm+LeY2I2ACUA+jIyegoU0m/vkqpDqWUZ/jTXwFYmqOxUW4Y+R6nMUwp1aOU6hv++FkAdhGpzvOwKAUiYocW0B5VSv0pxiXj+vuYIS2+DQBOEJGZIlIE4AYAT4+45mkAnxz++FoALytuPDdWJP36jpjXcAW0+RA0fjwN4P8Mrw47A0C3Uqol34Oi7BGRKfo8YRFZDu09jz9IjxHDX7tfA3hfKfXjOJeN6+9jW74HUKiUUn4R+SKA5wFYATyolNohIt8GsFEp9TS0fzyPiMgeaJNXb8jfiCkVBr++d4jIFdBWGHUCuCVvA6aUichqACsBVItIE4D/AGAHAKXULwE8C+CjAPYAGADwqfyMlNJl4Gt8LYDPiYgfwCCAG/iD9JhyNoBPANgmIluGb/tXANOAifF9zBMHiIiIiAoQ251EREREBYghjYiIiKgAMaQRERERFSCGNCIiIqICxJBGREREVIAY0oiIiIgKEEMaEY1LIlIlIluGfx0VkSPDH/eJyH0mvN5DIrJfRD6b4JoVIrJTRLZn+/WJaPzhPmlENO6JyDcB9CmlfmTiazwE4K9KqTVJrpsxfN0Cs8ZCROMDK2lENKGIyEoR+evwx98Ukd+KyOsiclBE/kFE/ktEtonIc8PnBkJElorIqyKySUSeH3FkWLzXuU5EtovIVhF5zew/FxGNPwxpRDTRzQZwPrTzWX8H4BWl1EJoxwhdNhzU7gVwrVJqKYAHAXzHwPN+A8AlSqnFw89NRJQSnt1JRBPd35RSPhHZBu0c1+eGb98GYAaAeQAWAHhx+KxuKwAjBzi/CeAhEXkCwJ+yPWgiGv8Y0ohoovMAgFIqKCK+iAO4g9D+jxQAO5RSZ6bypEqpz4rI6QAuA7BJRJYqpTqyOXAiGt/Y7iQiSmwXgBoRORMARMQuIvOTPUhEZiul1iulvgGgDUCjyeMkonGGlTQiogSUUl4RuRbAz0SkHNr/m/cA2JHkoT8UkROgVeL+DmCrqQMlonGHW3AQEWUBt+Agomxju5OIKDu6Afy/ZJvZAvgLgPacjYqIxixW0oiIiIgKECtpRERERAWIIY2IiIioADGkERERERUghjQiIiKiAsSQRkRERFSA/n/aboCim2s4UwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "n_points = 200\n", "start_threshold = 0.1\n", "stop_threshold = 0.8\n", "\n", "# Use simple thresholds to get the start and stop of the signals\n", "start = get_acceleration_start(df['acceleration'], n_points_for_basis=n_points, threshold=start_threshold)\n", "stop = get_acceleration_stop(df['acceleration'], n_points_for_basis=n_points, threshold=stop_threshold)\n", "\n", "# Plot the acceleration with starts and stops\n", "fig, ax = plt.subplots(figsize=(10, 6), ncols=1)\n", "ax.axvline(df['time'].iloc[start], -2,2, color='green', label='start')\n", "ax.axvline(df['time'].iloc[stop], -2,2, color='red', label='Stop')\n", "ax.set_ylabel(\"Acceleration [g's]\")\n", "ax.set_xlabel(\"Time [s]\")\n", "ax.plot(df['time'], df['acceleration'], color='darkblue', label='Acceleration')\n", "a = ax.legend()\n" ] }, { "cell_type": "code", "execution_count": 16, "id": "dab31f1f", "metadata": { "tags": [ "nbsphinx-thumbnail", "nbsphinx-gallery" ] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Index(['time', 'Sensor3', 'acceleration'], dtype='object')\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ4AAAJOCAYAAACHqmG9AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABxvklEQVR4nO3dd3hUZfr/8fc9k94DCSW9kAChhl6UqogVO1jRdXXtveva3XXXuurvq4uKYq9rFxELIkrvHUIPNUB6z8zz+2MOGCMlQGZOyv3ymiszzylzn3HCJ+ec5zxHjDEopZRSvuKwuwCllFItiwaPUkopn9LgUUop5VMaPEoppXxKg0cppZRPafAopZTyKQ0epRopEblIRL6zuw6lGpoGj2q2RGSjiJSLSImI7BCRN0QkzEfv/ZCIGBHpX8/5U6z5/fa1GWPeMcaM8kJtb1jv1a9WWwcRMbVeTxORv1rPh4mI2/oci0VktYhc3tB1qZZDg0c1d6cbY8KAnkA2cI+331BEBLgU2Gv9bIz2Ao8dwfzbrM8xArgFeEVEOnqlMtXsafCoFsEYswOYgieAABCRu0VknfVX/AoROavWtE0i0tt6fpG1h9DFen2FiHx2iLc7HmgP3AiME5GAWusNFpGnrfUXisgMEQkGpluzFFh7FgNF5DIRmWEt95KIPFX7TUTkcxG51XoeJyKfiEieiGwQkRsP85FMArqLyNDDzPcHxuMbPMHV/UiWVWofDR7VIohIAnAykFOreR2ekIgEHgbeFpH21rSfgWHW86HAemBIrdc/H+LtxgNfAh9ar0+vNe0poDcwCGgF3Am4a607yhgTZoyZWWed7wFjrb0pRCQaGAW8LyIO6/0WA/HASOBmETnpEDWWAf8AHj/EPH8iIg4ROQOI4Y+fpVL1psGjmrvPRKQY2ALsAh7cN8EY85ExZpsxxm2M+QBYC+w77/EznoABTzj9s9brgwaPiIQA5wHvGmOqgY+xDrdZAfEX4CZjzFZjjMsY85sxprIe2/ELYKxaAM4FZhpjtgF9gVhjzCPGmCpjzHrgFWDcYdb5XyBJRE6ux/vHiUgBUA58CtxqjFlYj+WU+hMNHtXcnWmMCcez99IJz1/qAIjIpSKySEQKrH9Uu9aa/jNwvLUH5MSz9zJYRFLw7CEtOsj7nQXUAN9Yr98BThaRWGvdQXj2tI6I8Yzm+z5wgdV0obVugGSsYKi1LfcCbQ+zzkrgUetxONuMMVF4zvE8D4w40m1Qah8NHtUiGGN+Bt7Ac6gLEUnGs1dwPdDa+kd1GSDW/Dl4DkfdAEw3xhQBO4CrgBnGGPdB3mo8EAZsFpEdwEeAP56g2A1UAOkHKrEem/EecK5Ve3/gE6t9C7DBGBNV6xFujDmlHut8HYgCzq7HvPvC6i6gm4icWZ9llKpLg0e1JM8BJ4pIDyAUzz/2eQBW9+Cudeb/GU8w7TusNq3O6z8QkX3nV07D04mhJ9AD+BdwqRVWE4FnrM4ATqsTQaBVhxtIO1jx1qGt3cCrwBRjTIE1aQ5QLCJ3WZ0XnCLSVUT6Hu4DMcbU4Dn8eNfh5q21TBXwNPBAfZdRqjYNHtViGGPygDeBB4wxK/D84zkT2Al0A36ts8jPQDi/9zir+7quS4BFxpjvjDE79j3wHJrqLiJdgduBpcBcPD3D/gU4jDFleE70/2odLhtwkPd4FzjB+rlvu1z8HnYb+D2cIg/7oXi8B2yv57z7TMRzfuj0w86pVB2iN4JTSinlS7rHo5RSyqc0eJRSSvmUBo9SSimf0uBRSinlU36Hn6V5iYmJMSkpKXaXoVqyygJwVULIIa/vVKpRmT9//m5jTGxDrKvFBU9KSgrz5s2zuwzVkuV8AQVroc9tdleiVL2JyKaGWpcealPK5wzWAAlKtUgaPEoppXxKg0cppZRPafAo5WvGgOihNtVyafAopZTyKQ0epZRSPqXBo5TPaa821bJp8CillPIpDR6llFI+pcGjlM9przbVsmnwKKWU8ikNHqVsoXs8quXS4FHK1/R286qF0+BRSinlUxo8StlCD7WplkuDRymf00NtqmXT4FFKKeVTGjxK2UGv41EtmAaPUr6mvdpUC6fBo5RSyqc0eJTyOR2dWrVsGjxKKaV8SoNHKaWUT2nwKOVrRkenVi2bBo9SSimfsiV4RCRKRD4WkVUislJEBopIKxGZKiJrrZ/R1rwiIs+LSI6ILBGRXrXWM96af62IjLdjW5RSSh0Zu/Z4/gN8a4zpBPQAVgJ3Az8YYzKAH6zXACcDGdbjKuAlABFpBTwI9Af6AQ/uCyulGjft1aZaNp8Hj4hEAkOA1wCMMVXGmAJgDDDJmm0ScKb1fAzwpvGYBUSJSHvgJGCqMWavMSYfmAqM9tmGKKWUOip27PGkAnnA6yKyUEReFZFQoK0xZrs1zw6grfU8HthSa/lcq+1g7X8iIleJyDwRmZeXl9eAm6LU0RB0oFDVktkRPH5AL+AlY0w2UMrvh9UAMMYYGvA30xgzwRjTxxjTJzY2tqFWq9TREQcYt91VKGUbPxveMxfINcbMtl5/jCd4dopIe2PMdutQ2i5r+lYgsdbyCVbbVmBYnfZph3vz1XmrGfbSsMPNppTXFOzZSLW7ktjFH9tdilK28HnwGGN2iMgWEelojFkNjARWWI/xwBPWz8+tRb4ArheR9/F0JCi0wmkK8I9aHQpGAfcctoAyYH5DbpFSR2Yb+ZRIObFJdleilD3E2DBSroj0BF4FAoD1wOV4Dvt9CCQBm4DzjTF7RUSAF/F0HCgDLjfGzLPW8xfgXmu1jxtjXj/ce/eRPmYe8xp2g5Q6AudHDmSZcyMr9mw//MxKNRIiMt8Y06ch1mXHoTaMMYuAA23AyAPMa4DrDrKeicDEBi1OKS9z4cahIxeoFkxHLlDKx9y4cYr+6qmWS7/9SvmY2xjd41EtmgaPUj7mwoXDob96quXSb79SPubG4NTgUS2YfvuV8jGXceFw6KE21XJp8CjlY26Mdi5QLZp++5XyMTduHPqrp1ow/fYr5WMu49ZDbapF0+BRysdcuHE6nXaXoZRtNHiU8rEyU0lIUKDdZShlGw0epXysBhf+fraMVqVUo6DBo5SP1RgX/k4NHtVyafAo5WPVxoWfQ4NHtVwaPEr5WI1x4+ennQtUy6XBo5SP1RgXfnqoTbVgGjxK+VilqcE/QINHtVwaPEr5WKmrirCQULvLUMo2GjxK+ZDBTZmritBwDR7VcmnwKOVD5VKMAULDNHhUy6XBo5QPlTkLAQgJD7G5EqXso8GjlA+VO4sBCAoJtrkSpeyjwaOUD1VSDkBQmAaPark0eJTyIZfUAOjo1KpF0+BRyodKHZ5zPKGR4TZXopR9NHiU8qEC524Aotq1trkSpeyjwaOUD+U78gCIbtfW5kqUso8Gj1I+VCB7AYhqHWtzJUrZR4NHKR/KxxM80THtbK5EKfto8CjlQ/mSj1MchEW2srsUpWyjwaOUDxVQSJR/EOLQXz3Vcum3XykfyncXEuWvw+Wolk2DRykfKjDFRAfpAKGqZdPgUcqH8l3FRIWG2V2GUrbS4FHKhwrcpUTrqAWqhdPgUcqH8l2lREVG2F2GUrbS4FHKR9y42F1dRpt2evGoatk0eJTykT3+W3AZN23jdbgc1bJp8CjlI7v9tgIQk9DG5kqUspcGj1I+stu5DYDY1DibK1HKXho8SvnIbucOAGLi2ttciVL20uBRykf2yC4AWsdq8KiWTYNHKR/Jt26JENlKz/Golk2DRykf2W62E+LnT0Qr7U6tWjYNHqV8ZKfJo22IXjyqlAaPUj6y1ZVHXGS03WUoZTsNHqV8ZEPNLtLax9tdhlK20+BRygeqqSC3ooCU5CS7S1HKdho8SvlAbuAq3BhSOqbaXYpSttPgUcoHNgasACClW6bNlShlPw0epXxgo3MNACndsmyuRCn7afAo5QMbZQMCJKR2srsUpWynwaOUD2w0m4kPjSQgKNjuUpSynS3BIyK3iMhyEVkmIu+JSJCIpIrIbBHJEZEPRCTAmjfQep1jTU+ptZ57rPbVInKSHduiVH1sqNlGig6VoxRgQ/CISDxwI9DHGNMVcALjgH8BzxpjOgD5wBXWIlcA+Vb7s9Z8iEiWtVwXYDTwfyLi9OW2KFVfOVU7yEjSrtRKgX2H2vyAYBHxA0KA7cAI4GNr+iTgTOv5GOs11vSRIiJW+/vGmEpjzAYgB+jnm/KVqr8KKWF7VRHJyYl2l6JUo+Dz4DHGbAWeAjbjCZxCYD5QYIypsWbLBfZd4h0PbLGWrbHmb127/QDL/IGIXCUi80RkXh55DbtBSh1GTtB8ADJ6drS5EqUaBzsOtUXj2VtJBeKAUDyHyrzGGDPBGNPHGNMnFh0ZWPnWqoAFAHTs39PeQpRqJOw41HYCsMEYk2eMqQb+BwwGoqxDbwAJwFbr+VYgEcCaHgnsqd1+gGWUajQWOxbgECGr1yC7S1GqUbAjeDYDA0QkxDpXMxJYAfwEnGvNMx743Hr+hfUaa/qPxhhjtY+zer2lAhnAHB9tg1L1tti1go7RbQkO01siKAWek/w+ZYyZLSIfAwuAGmAhMAH4GnhfRB6z2l6zFnkNeEtEcoC9eHqyYYxZLiIf4gmtGuA6Y4zLpxuj1GEY3MyuWMPJvfvaXYpSjYbPgwfAGPMg8GCd5vUcoFeaMaYCOO8g63kceLzBC1SqgWwMXMKuyhL6H6fBo9Q+OnKBUl40J+h7APqfNsLmSpRqPDR4lPKi2fIbgQ4/uvUfancpSjUaGjxKedEPVXMYmNwB/8Agu0tRqtHQ4FHKS3b4r2NJ2VZOGq57O0rVpsGjlJd8EToRgFP/Os7mSpRqXDR4lPKS912fkhEVS9f+Q+wuRalGRYNHKS9YH7iQacUrufjU0xCH/popVZv+RijlBc+FPISfOPnro3fbXYpSjY4Gj1INbHvAWl4rmsxFA4cSl5ppdzlKNToaPEo1sDuC/0qNcXPvi/+0uxSlGiUNHqUa0M9h7/NO4XTuPOsCMrL1voRKHYgGj1INZHnwdM6vupKU8Fbc8/qLdpejVKOlwaNUA1gRPIMRNafidDj4dsrXhIRH2l2SUo2WBo9Sx+i30E8Z7joZp9PBT99OoePAAXaXpFSjpsGj1FGqlFLui76I40vPISQogJ9+nELHoRo6Sh2OBo9SR2F26JdkB3fgH/nvculxw1m0foXu6ShVTxo8Sh0Bg5t7oi9gUOkYip3lTH7zNV7/5QciW7e1uzSlmgxb7kCqVFP1a9inPJH/PgMz0/l25kwiWsXaXZJSTY7u8Sh1BDIqs4kPjGTWmnWcO3AEk/7xCLu2brC7LKWaFA0epY5A2+o0FrmW8PfoS1m7eRuX3fcgbRPS6Ns+hQfOvZRZX32Bq6ba7jKVatTEGGN3DT7VR/qYecyzuwzVDLhxsSjkeyYHfszk6mnMLFmHG0OYXwAv3HUflz32gN0lKtVgRGS+MaZPQ6xL93iUOkoOnPQqO4kbCv/B3Y7buSRqGAAlNVW8+8FH9hanVCOmnQuUOgozQz/jq4CP+LFmFnNLNuIqchPo9GNEemdOHjmCyx680+4SlWq0NHiUOgI1VHFF1Em8WTANZ5mDfqmp3H3OBYw8+xQGjjyDoJAwu0tUqtHT4FHqCHwVMYE3C6ZxYo+ufPzjj9qdWqmjoOd4lDoCadVdAFi3YbuGjlJHSfd4lDqMnf7r+TT0FT5yfc604lUA9MvsZHNVSjVdusej1EFsCVzO5VEjia/J4JqCJ8j128M9Yy9k0YzveW/uDLvLU6rJ0j0epQ7iseDbeKPgR6468SSuf+RuuvYbgjj0bzWljpX+Fil1EBfWXAZAu5C2dBswTENHqQaiIxcodRAGN52D25NPGduK9uL087e7JKVsoyMXKOVFBjcfRT5D39BUVpfvIiY0DLfLZXdZSjUbGjxK1fFtxETOL7yNHY5CJtx/Hws25+AfGGR3WUo1G9q5QKk6SqUYgImvPMeosZfZW4xSzZDu8ShVh5/xnMup2a23N1DKGzR4lKpjp2MbAItz5um5HaW8QINHqTrGlF5Bz9AE7n1uAqnRsdx8ypl89cZ/2bF5nd2lKdUsaHdqpQ6ginI+jHyWD8zHfF+ylAp3DQDtQsLJTkoju0sXsof0IXvkENI6Z+s1PqrZa8ju1Bo8Sh1GhZQwJ/QbFvj9xkIWs7BqDSvKd+AybgCiA4M5PjWLYf0HM/TiU+gxbIRe86OaHQ2eY6DBoxpChZSwLOQXFvr/yixm8XPFYtZV7AagdWAIp/cawCOvPk9iVhebK1WqYegFpErZxOCmxLGXvIBNOIyDJFc6w80Irgu+kPMjBwKwp7KMN2b+yPN3PmRvsUo1Unodj2r23LgodeRT7LeXYsdeip35lEghxY58iqWQYimiRIrZ/58pociUUGRKKXKXUuguo8hV7nnUVOJ2G6j88/sEO/3p0ro9nROSGDX6BC689zbfb6xSTYAGj2p0XFRT4txLiTPfEw7OvRQ7CiiRIoodBRRT5AkLSig2xRRTQrEppYQyik0Zxe5ySlwVFLsqKHZVUuqqAjdQdfj3DvULIDwgkIiAYCKCg4kMDqVNaDQRIWFEhocTERVBRHQEkTFRtEluT2xye2LjE4htl0hUTDvtZKBUPWjwqAZhcLPLfyN7/bZ7wsLam/CERSGeeCix9ilKKKaUEiskit1lFLsqrLCopMxVDS48j8MI8w8gPDCI8MAgwoKCCA8KJj4khvDQUMJCQggPDyM8IpzwyHDCIsMJbxVBeOtIwltHExYRSXhENOFRrQiPbEVoRDQOp9Prn5VSLZ0GjzpqJY69XB9xLh8Vz8AYKK+uhkNc7C9AWEAg4UGBhAUGER4cRHhICImhbQgLCyE8PJTwiDDCo8IJiwojPDqC8PAIwiOiCAuPJDwyivDIaMIjWxEe1ZqQsEgNCqWaIA0eddSWBP/MpIKf/tAW5O9Hx7h2ZCYnkN4hhfTOHUjL7ERax24kpHXCzz/ApmqVUo2FdqdWx2R58HR+DvyapSxlac1aNlbnsaOqeP81LrU5RGgdGEKbkAhiwyNpExlNbKvWtGkbS2xcLG1S2tMmLZ7Y+ATaxCURFdNO92iUaiT0Op5joMHjfW5c5PlvIs9vC7v8cslzbGeXYzt55LHL5JFn9rLLlU+eq5Bd1cXk15QfcD1OcRAbFEpscDhtwqKIjYimTXRrYlvH0qZtW2Lj2tEmJY7YDgm0yUgmsnUbPbmvlJc0ZPDooTbV4Bw4aVudRtvqtHrNX00FuwO2kOfMtYJqG7tkB3liBZVrL7v25jN/1xp21RRTWFNxwPX4i5PYwFBiAyNoExJFbEQUbaJjiI2JpU1cW2LT2tGmezzp3XsQ2z5ZQ0opm2jwKNv5E0T7qgzak1Gv+SullN3+W8jzy2WXM5c853ZPUO3bo3Lns6son/X5O9iVU0yx688X3bQKDKFz+zg6p6SS1aUz3Yb0oe/IUUS2btvQm6eUqkMPtalmr0JKyAvYRJ4zlx1+m1nrWMFKVrHStZ4VFVvZXV26f96usXF89fNkkjt3t7FipRqfJnGoTUQmAqcBu4wxXa22VsAHQAqwETjfGJMvIgL8BzgFKAMuM8YssJYZD9xvrfYxY8wkq7038AYQDHwD3GRaWoqqegkyYSRWdiGRP4+blu/cxrtRz/J42QS2VxWxLG8bs7+cqsGjlBd58yD3G8DoOm13Az8YYzKAH6zXACcDGdbjKuAl2B9UDwL9gX7AgyISbS3zEnBlreXqvpdSf+CimlVBM3kv8l/cGX0+x4VlEONK4PqCp6iUGi7sdxwfPfcM591+i92lKtWseW2PxxgzXURS6jSPAYZZzycB04C7rPY3rT2WWSISJSLtrXmnGmP2AojIVGC0iEwDIowxs6z2N4Ezgcne2h7VtFRK6f7RoxfKQhbWrGJx2WbKKqqhAvwdTnq2TeTeERdxyvhz6HfGqXorA6V8xNedC9oaY7Zbz3cA+87kxgNbas2Xa7Udqj33AO0HJCJX4dmTIomkYyhfNTYGN3n+m1gZOJuVfotZKkuZU7WMxWW5VJd6xtwJDwikZ2IyV3Q6meyB2WQPP46sXoMJCAq2uXqlWibberUZY4yI+OScjDFmAjABPJ0LfPGeyjsKnTv5Knwi85nLItcqFldsYm912f6hesL8AuiTkMYt3c+hz6j+ZA/13CFUL0RVqvHwdfDsFJH2xpjt1qG0XVb7ViCx1nwJVttWfj80t699mtWecID5VTNW5Miji19HthYUEuT0o3u7BM7JPo6srp3pPLAHnfv3JSGts4aMUo2cr6+g+wIYbz0fD3xeq/1S8RgAFFqH5KYAo0Qk2upUMAqYYk0rEpEBVo+4S2utSzVTRX67KbBGOeidkMKQXn3oe3xfep40gK6DBmroKNVEeO06HhF5D8/eSgywE0/vtM+AD4EkYBOe7tR7rfB4EU/PtDLgcmPMPGs9fwHutVb7uDHmdau9D793p54M3FCf7tR6HU/Tti1gDf8v9GGmVv/G4rItVLl/v3dCiJ8/6VGxpLeJo0NyCh2yMujQtzPpPXuQmN5ZOw8odQx0rLZjoMHTfLioZkvASnICF7PWbxlrzVrWuTeTU7OVdRW7qXTX7J83wOEkM7ot3ZLS6JrVmW5Ds+k2ZDDJGV116Byl6kGD5xho8LQMblxsDVzFOv+l5PgtZ62sYYVrLUurNrKpIn//fBH+gfSIS6Zn105kj8imZ/9BZPUaTGBwqI3VK9X4aPAcAw0eVeTIY1nwDJb6z2WJLGFR9SoWl2/x3CIbzzU+WW3a0yUxhcT4eBJSEonvmERCxzQS0jJpE5+ih+1Ui6PBcww0eNSBuKhmXdACFgXOYJFjPgtdy1lVlcvWykKqzR/vwe0nDuJCI0mIakVCbBuSExJI79jBcz6ph55PUs2TBs8x0OBRR8KNi93+m8n1X0Ou33q2OjeSSy65Zhu5rp1sqdnNpsq9f+jkEOBw0rN9Iv17dqf/SQPoP/xE0rN66bkk1aRp8BwDDR7V0PadT8oJWMw650pWyyrm1ixlbtl6ylyeK1tbB4XQLz2D/gOz6T/qePoNOYlWbQ862IZSjY4GzzHQ4FG+UkMVy0N+YXbAj8xmNrOrVrCibDv7fuMyYmLp37kz/fr2ISkrjdikONq0TyS2fSIR0bG6h6QaFQ2eY6DBo+xU5MhjXuh3zPabzmyzgNllq9lRVfyn+QIcTmKDw4gJCScmPILW4ZHERLUipnVrWrdpTUxCG2LS25GS1ZmUzO467pzyOg2eY6DBoxoTg5sdAevY7reeXc6t5Dm2s8uxg13sIs/sYY+7gN3uAna7itldU0J+dRl1f2MdCImhUaS3akfnjHRue/ZRUrv3tGNzVDOmwXMMNHhUU+aimny/7ez220qeXy4bnKvJcazhp+pZzChZC0DftinM2bHB5kpVc9Mk7kCqlGp4TvxpVRPPV2Fv8L35gfXurSwryaXYVbl/nqH9GuTfBqW8RoNHqSZmbug3XF7wIADHp2VyaYcT6ZbdjW7D+5I9eCTBYRE2V6jUoWm3GaWamJiauP3PN+7ehZ/DSduEODK799bQUU2CnuNRqgkqcezli/AJvGs+4IeSZVRYA6J2bxPPiD59GX7OKIaccRZRMe1srlQ1F9q54Bho8KjmplJKmRf6LT/7T+Gnmt+YUbp6fxDFhUaQFtOW1HZxJCcmkpiWRFJWOolZHUjqkEV4dIzN1aumQoPnGGjwqOauUkqZFfolv/h/T45ZzwbXVjZU7WJrZSHuOp2xI/2DSAprRWJ0LEnt4khMTSAxK5nEbunEp6QRn9KRkPBIm7ZENSYaPMdAg0e1VDVUsT1wLVv81rDZL4ctjo1sNpvZ4t7B5pqdbKney+7q0j8t1yowhKTIViTFtiWpfXuSUpJI6phKcrdMkjp2ol1iut75tQXQ4DkGGjxKHVyZFJIbuIpc/xxyHRvIdWxmi8lli3s7m2t2sblqL4U1FX9Yxt/hJCEkiqSoGJLataNDxzS6j+hF94GDScnsrqHUTGjwHAMNHqWOTaFzJ1sCVrHZfw2bHevZLJvY7N7KZtcONlXnsaUyf/8BvXD/QLrFJdA9PZOOnTJJ69qR9D5dSO3UQw/hNTEaPMdAg0cp7yp15LM8+FcW+89iCYtZ4lrDkvLNFNTZU2ofFEFaRBvSYxNIS0whvVMG6cdlkdavO23iU3SQ1EZGg+cYaPAo5XsGN3v9trIucDHr/VayzrGG9WYD62pyWVe1g61VhX/o9pAZ3YbhPXvRZ2Af+pwylC59jsM/MMi2+pUGzzHR4FGq8amQEjYGLmV9wDJWO5cxueZH5pTn7D+fFOjwo3u7BPpkZtFvaH/Ouela7QruYxo8x0CDR6mmwY2LdUHzmR/4M/Mdc5hXs4wF5Rsoqqmkc6t2fDX9G9K6ZNtdZovRkMGjB1GVUo2SAycZFf0YV3gHT+Z/xE/FK8mvKeWziBfZVlRAZrc+vP+PJ+0uUx0FDR6lVJPhwMmYouv4Z/iNuIyblYuW2V2SOgoaPEqpJud79zQANu/KZdWCmbhdLnsLUkdEb4uglGpynq14jdCo63lr+jTe6D2IYD9/Osa2pXNiCp2zMuncvytZffvSoUtvvS14I6SdC5RSTdaWwOV8F/wBy1nGSvc6VlZtYVNF/v7pDhFSwluR0TaOzKRUMjtkktEji8zjs0nq1AWnn7+N1Tct2qvtGGjwKNW8lTryWRU8m5V+C1jjWMUa9zrW1uSypmIHJa6q/fO1DQnjhnPHcc0TD9GqfbyNFTcNGjzHQINHqZbJ4Gan/wbWBM5njd9SPnF/wbdFS0gJb8X6gjwdKeEwNHiOgQaPUgqginICCQHAVVOjg5kehl7Ho5RSx+i16If2P5/y7uvUVFcdfGbVoDR4lFIt0gVFt/C3qNHE+odyyqVXEhfRijtOP49lM3+2u7RmT4NHKdUiRbna8XLBZHKr8/gk4j8MDuzMc1//j26DhnHr6LOprqw4/ErUUdHgUUq1aAEEc3bRjXxaOJftzg1cH3Uaz075lOtHnolxu+0ur1nS4FFKKUtMTRIvFHzJkPBMJvw6haFpnfnwmacpLymyu7RmRYNHKaXq+LZkLi9G3U7O9h2Mve12eiR0oGhvnt1lNRsaPEopVUullLIuaCExpi0XhJ4AwNrCPCorSm2urPnQsdqUUi1OkSOPjYHL2OS/ii2OjWyWTax1b2B51UZyKnbjKndDOfg7nIzO7M4dD95JbFyK3WU3Gxo8SqlmyY2L9YELWRL4G0ucC1jqXsmGmu1srNpNfnU5lON54AmY1MjWdOmQzHlZo8jq35Uug/qR2a0vQSFhtm5Hc6TBo5Rq8ty4WBD6HXP8p7GYxSxxrWVpWS6llVVQ6RksND0qhg5JcQxM7E5KajIpXdJJ7ppJYlpH2iak6cgFPqTBo5RqslYFzeShoNv5qWwJu0pLAIgOCqZ7YhJXdB5N94E96D6gP116H0dIeKTN1ap9NHiUUk3WhKAn+aDgNwBevPt2zrjiMhLSOuuAn42cBo9SqklaFTSTZws+BcApDkZdcB6JHbrYXJWqDw0epVSjYnBT4NzBLv/N5PltJc+5jV2ynTzZxTaznRzXZnKqt7OhYg8AN503hsdemUhYZCubK1f1pcGjlPKpQudOlgbPYJXfIjbKBja6t7DDvZs8VyF5NcXkVZdQ43KD68/LRgcE06FVWwbEdebKblmcNP5seg0/yfcboY6JBo9SqkGVSSFbAley2X81m5w5bJZNbDa5bHBtZX3VTnIrC8HTDwCnOEgIiyQushUpEe3oG9WZNjGxxLaNJTahLW3S44hNjiO2fSKxcUkEBofau3GqQWjwKKWO2fzQb7nUfTmbqvZS6qqCCjwPPF2Z40IiSGnVhuHte5LVpRPdh/Sm84C+JKZn4ecfYGvtyvc0eJRSx2xawNesyN+x//WN487g3PEXk5TakbiUTPwDg2ysTjU2eutrpdQxM7j5NexTnnM+ySeFswHoHhvPaUOG0m/0IDr360Na52zdu2nCGvLW1xo8SqkGtS1gDR+G/B+fuiYzo3gtbjz/xvg7nHSIjKVT20Q6J2fQqXMWnYb2pMtIvbizKdDgOQYaPEr5ToljL8uDf2WV/0JWyUpWunJYVb3ZMxCn8dxkrVVACP+65Vb++sSjNlerDkWD5xho8ChlvyrKWR+0iBUBc3nRvMxPxSvp0qY9p40Ywunjz2PACWfg9PO3u0xVS5MIHhGZCJwG7DLGdLXangROB6qAdcDlxpgCa9o9wBV4eu/faIyZYrWPBv4DOIFXjTFPWO2pwPtAa2A+cIkxpupwdWnwKNW4uKjmv1H384n7S6YXr6bGuGkdFMIpffsy+szRdO7Xi/SsbCJaxdpdaovWVIJnCJ7e+m/WCp5RwI/GmBoR+ReAMeYuEckC3gP6AXHA90Cmtao1wIlALjAXuMAYs0JEPgT+Z4x5X0ReBhYbY146XF0aPEo1XoXOnUwJe5uv5Eu+KZ3Lnuqy/dNaB4aQFt2G9HbxpCUmk5aZTnrfLNJ6dSU+taPuIXlZkwgeABFJAb7aFzx1pp0FnGuMucja28EY809r2hTgIWvWh4wxJ1nt91htTwB5QDsrxAbWnu9QNHiUahpcVLMsZDrr/Jazzrma9WYD61xbWF+9g02Ve6mxzhGBp+NCSngr0mPbkxaXSFp6KuldO5LWrwtp3XrocDoNoCGDx87reP4CfGA9jwdm1ZqWa7UBbKnT3h/P4bUCY0zNAeb/ExG5CrgKIImkYy5cKeV9TvzpUTaSHoz807QaqtgSuIL1AUtZ71zNOslhvXsT67dvY9bGNRRMr/jD/G2CwkiLakNa2zjSk1NJ69aBDoO60LXfYKJi2vlqk5TFluARkfuAGuAdX7yfMWYCMAE8ezy+eE+llPf4EUBqZU9SK3seIJYg37mN9UFLWO+3gnWyhvVsYH1pLr+tXsH7i3/F/cXv/wy0D4kguVUMKe3bk5ycQHKnZJK7diClQ2eSM7oQGhHtuw1rIXwePCJyGZ5OByPN78f5tgKJtWZLsNo4SPseIEpE/Ky9ntrzK6VauGhXHL1L4+jN6D9Nq6aCzYHLWR24gCWO+awhh02F25i7ZzWfzJtFtfnj6KStg0JIjo4hpU07kuPjSUxKJKFjMgld0khIzaB9coZeGHuEfBo8Vg+1O4GhxpiyWpO+AN4VkWfwdC7IAOYAAmRYPdi2AuOAC40xRkR+As7F07NtPPC577ZEKdVU+RNEemVv0it7cwpX/mGaGxfbA9ayyX8lm/zWsEk2sInNbCzdysq1m/l22ULKXNV/WMYhQrvgcBKiWpHStj0ZqWl0yMqkQ98udOjZk7YJqXpjujq82avtPWAYEAPsBB4E7gEC8eyxAMwyxlxtzX8fnvM+NcDNxpjJVvspwHN4ulNPNMY8brWn4QmdVsBC4GJjTOXh6tLOBUqpo7XvXkG5AWvI9c8h17GRXNlCrtnKFtdONtTsYEPF3v0XxwKEOgPoENGGDm3i6ZCcQoesDDoc14WM7J7EpWQ2mVBqMr3aGiMNHqWUN1VTwabAZeQELCHHuZIcWUuOaxM51dtYX7HnD4fy2gWHMzgti8H9+nPc+aPofeJoHE6njdUfnAbPMdDgUUrZxUU1WwJWkhO4mFXOJcwyc/i1chkbK/YC0Kd9Mq9Nep3uJw63udI/0+A5Bho8SqnGZlvAGiaHvM19Zc9S5KrkznEXctvzTxHeKsbu0vbzWfBYF2ZeDBwPtAfKgWXA18DbxpjChijClzR4lFKN1baANYxwHs/q8l04xcH6FYtJ6vSn6+9t4ZMLSEVkMrANT2+xx4FdQBCeoWyGA5+LyDPGmC8aohCllGopDG7y/DexMnA2K/wXssKsYLlrHfPL11NU5ekjFeznR3BYmM2VesdB93hEJMYYs/uQC9djnsZG93iUUr6w73qhddboCutlHevcm1hfs431lXkU1fzeCTfcL5Cs2DiyO3Vi0OjBDDp1NGmdsxtVjzef7PHUDhQRSQYyjDHfi0gw4GeMKW5qoaOUUg2tUkqZH/IdC/1/YzFLrPHkdrK5Ih93pQErXwKdfqRGtiatfTuOT+xFWmYqWf16kjWoP/FpHRtVyHjbYS8gFZEr8Yxz1gpIxzNKwMtwwJEqlFKqRdgcuIxbg/7K18XzqSj1DBsZHRhMx5j2DGqfxcXJKaR2SSO9T2fSs3oQl5LZaLtK+1p9Ri64Ds/tCmYDGGPWikgbr1allFKN3MPBN/NJwWwAbrn4LG647wFSMru3qD2Xo3XY7tQiMtsY019EFhpjskXED1hgjOnumxIblp7jUUo1hM2By7gm8GKmFi+j2rgIdPrROaYdGe0TSE1MIi09nbSuHUnt04Wkjl0ICAq2u+Rj4uvbIvwsIvcCwSJyInAt8GVDvLlSSjVVSZVd+bpyEcWO3XwX8TazHb+xtHw1i1av47PFc/8wQoEDISEkirRWbenYIYWeg3rSf/Rwegwa2SIPv9Vnj8eB55bUo/AM2jkFzy2om+SVp7rHo5TyNhfVbAtcwwb/5az3W8UGWcd690bWu7ayomIrBTWe+wWlRrTmX/ffy3l33GpzxYenIxccAw0epZSdDG42BSzjjdBneDh/EgAzPv2YwWeeY3Nlh9aQwXPYs2AicpqILBSRvSJSJCLFIlLUEG+ulFItjeAgpao702vm7G+rcVRTU11lY1W+VZ9DbTnA2cDSpnp4rTbd41FKNQaLQ37gNueN/FyyihrjJtjpT3b7JPp27UrfYf3pMWIQGV37EBgcanepgI8PtVk3XBtpTK0bTDRhGjxKqcZkr3Mr34W9y2zHb8ytXsaCsk2Uuz03m3OKg06t2tI9KZ0Bg/sy/sG7iIxpa0udvg6evsCjwM/svwYXjDHPNEQBvqbBo5RqzGqoYkXwrywPmMMyxxKWuFaypHIjmyvz6dSqLQ/ddzsnXzaeiFaxPq3L18HzHVACLAX27/UYYx5uiAJ8TYNHKdUUfR/+FhdWXENedSmjOnRlytqlPn1/X1/HE2eMaRzjciulVAtTQxW/hH/M584PcePZUTj1pBNtrurY1Cd4vhGRUcaY77xejVJKtXAGNxsCFzM/6CemOX7k09Jf2F5cRLCfP6O79GTspedz/s2N/7qfQ6lP8FwD3C4ilUA1notIjTEmwquVKaVUC7HbbzM/hX7M1/IV35TOIa+yFCohxM+fUZ17cOHfLuCUSy4nNCLa7lIbxGGDxxgT7otClFKqpdkasIrLg85jatEyKITogGBO7tKbIWceT+/hx9Ot75BG0526IR3qDqSdjDGrRKTXgaYbYxZ4ryyllGr+Eqo6QxXcf9HFnDz+bPoNOxU//wC7y/K6Q+3x3IrnPjxPH2CaAUZ4pSKllGoB/hN1CxR4nt/z0ouEhEfaWo8vHSp4lgAYY4b7qBbbVEdXk/tQLhUdKuoxiJCqFzcE5QSR8FAC/vn+dlejVKMTStj+5/3Tsnj99Vfpc9rJNlbkOwe9jkdEFhhjDniYrSk70HU8G/6zgfB+4bT2a40gNlXWvBgMe2r2UDynmNSbUu0uR6lG6fvwtzix+FIAusXGsWTXVpsrOjifDhLaElR0qNDQaWCC0NqvtWcvUil1QBuda/Y/f/PDSTZW4luHCp7u1mjUdR/Nb3RqBxo6XiCI/mmj1CH8teDR/c/nff8zZcWFNlbjO4f6Z2GpMSbiAI9wvYan6Zq3Yh43PnVjg6/3oQkP8dRbTzX4epVq7j6J+A/dQuK48vHHiImO5eqhJ1NVUW53WV6lf482MsYY3G7vDQTeJ6sPz9/+vNfWr5Q6MmcX3ci3NT9xSdRQyl3V/Hf6t6xdOtfusrzqUMHzkc+qaOE2bttIx3M6cumDl9J1bFe27NzCNU9cQ59L+9Dl/C48+N8HAZi7fC5n33E2AJ///DnBxwVTVV1FRWUFaWPS/rTej77/iK5ju9Ljwh4MuWoIANPmT+O0W04DIC8/jxOvO5Eu53fhr4/9leTTk9ldsJuN2zbS+bzOXPnYlXQ5vwujrh9FufUX2CufvkLfS/vS48IenHPnOZRVlPniI1Kq2TG4eSb6RrJDE4mv6shbBT+TEtGKT//9Al36DrG7PK86aHdqY8w/DjZNRE4zxnzlnZJstvdbqNrRsOsMaAetRh9ylrVb1jLpoUkM6DYAgMeveZxWka1wuVyMvHYkS9YuIbtjNovWLALgl4W/0DWtK3NXzKWmpob+Xfv/aZ2PvPoIU16YQnybeAqKC/40/eFXHmZEnxHcc/k9fPvbt7z2+Wt/qOe9x97jlftf4fx7zueTHz/h4lMu5uzhZ3PlWVcCcP9L9/Pa569xw9gbjvKDUarlclHD7fkvYICeqQk8+/yTDD3lfMTR/A9EHe0W9m3QKhTJ7ZP3hw7Ah99/SK+Le5F9cTbL1y9nxYYV+Pn5kZ6QzsoNK5mzYg63XnQr0xdM55dFv3B8z+P/tM7BPQZz2cOX8cqnr+Byuf40fcaiGYwbNQ6A0YNGE11rHKjUuFR6duwJQO9Ovdm4fSMAy9Yt4/grj6fbuG688+07LF+/vAE/BaVaDj8C+CjyWWL9Q1m0IZcRp19An/hUbh99Ll+9+F8K9+y0u0Svqc8goX9ijHmwoQtpNA6zZ+ItoUG/j8e0YesGnnr7KeZOmkt0RDSXPXQZFZWebslDsocw+bfJ+Dv9OaHfCVz28GW4XC6evOnJP63z5XteZvay2Xw942t6X9qb+W/Or3c9gf6B+587HU7KXZ5DbZc9chmfPfkZPTJ78MaXbzBt/rSj3GKl1DmFN3Mm1zEv9FumBnzO96UzeGHqZzw95ROCbrqeCQ88xCUP3md3mQ3usHs8InL2AR4jRaSNLwpsiYpKiwgNDiUyLJKde3Yyeebk/dOOzz6e5957joHdBxIbHcuewj2s3ryarul/vmXSutx19O/an0eufoTYqFi27Nzyh+mDewzmw+8/BOC7Wd+RX5R/2NqKS4tpH9Oe6ppq3vn2nWPcUqWUE3/6l57O/fmvMq14FQVmLz+Gv4MBLn3ofh467zKK9ubZXWaDqs8ezxXAQOAn6/UwYD6QKiKPGGPe8lJtLVaPzB5kZ2bT6bxOJLZJZHD3wfun9e/Sn517dzIk23PysXuH7uzYswORP1+HdMd/7mDtlrUYYxjZdyQ9Mnvw84Kf909/8MoHueC+C3jrm7cY2G0g7Vq3IzwknJKykoPW9ujVj9L/8v7ERsXSv2t/ikuLG3DLlVLBJoLhxReyKqArNwRdxsMfT+K1bz/j7GFDGXz6EAafcCrxaZ3sLvOY1OfW11OAS40xO63XbYE3gQuA6U3t7qQHGjJn5eSVdI7pbFNF9qmsqsTpcOLn58fMJTO55olrWPTuogZ9j5W7V9L55Jb32SrVUGaFfsH9jnv4rXQt5e5qAFIiWjEwM4usTh3p1KMLHQdmk5Hdh6CQsMOs7ej5+tbXiftCx7LLatsrItUNUYSyx+Ydmzn/nvNxGzcBfgG8ct8rdpeklKpjQOkZfM8ZVFPBotAf+NV/Kr+aWfy6bDnvzZsBb3vmEyA1tDWd2ifQqWM6Xfp0Jeu43mT1GkhEq1hbt6Gu+gTPNBH5it+v6znHagtl/6DeqinKSMpg4TsL7S5DKVUP/gTRt/RU+nIqN1ttpY581gbNZ5X/QlY5lrPKvZZV2zbz4/rlVHz9v/3LpkfE0Dspg97detL39OMYct55OP3sGzW+PofaBDgbOM5q+hX4xBxuwUZKD7X5lh5qU8r3XFSzIXAxKwLnstS5gOk1s/iueNn+6Y9c8Bf+/u5rh1jDn/n0UJsVMJ9YD6WUUo1UNRWsC1rIqoD5LHMuYrF7GYsr1pNT/nuvuNjgUAaOGmZfkRzldTwiMsEYc1VDF6OUUqp+KqSE78Pf4Rfnj6xyrWN19RbWVeympsIN1t1I0iJa07NjOpf0PJ3sIf3oPWIE7ZMz7C2cowwe4L8NWoVSSql6q6GKvsEdWVa0jQCHk4zoNnRNTeXcjifSsVdnOvbrQVavQYRFtrK71AOqd/CISIgxpgzAGFP/S+BVk/TQhIcICw7j9ktut7sUpVQdgrCsbBsA3RMSuO7KvzBy7DkkZnSxubL6qc/IBYNEZAWwynrdQ0T+z+uVqf1qamrsLkEp1Yg48afEsZe7o8execceLv/7gyRldqVv+xR+fKPxjyhSnz2eZ4GTgC8AjDGLRaR5j9ltgze/fpOn3n4KEaF7h+44HU6CAoNYuHohg3sM5tJTLuXqJ66mrKKM9IR0Jv59ItER0Qz72zB6ZHhGJKhx1TDxgYn069KPhyY8xLrcdeTk5rC7YDd3XnLn/lGln3zrST6c+iGV1ZWcNewsHv7bwwA8PvFxJn09iTbRbUhsm0jvTr3t/EiUUocQ6o7mn/nv8Tgulgf/wndBn3D7jhcZefnFfFVRyqlXN97T8PU61GaM2VJnSJY/D3XcTNz89M37bz3QUHpm9uS525476PTl65bz2MTH+O2134iJimFv4V5uffZWcnfl8ttrv+F0Oul+QXdeuP0FhvYeygMvP8DDrzy8f51lFWUsencR0xdM5y+P/IVlH3i6TS7JWcKsibMorSgl+6JsTj3uVJatW8bazWuZM2kOxhjOuO0Mpi+YTmhwKO9/9z6L3llETU0NvS7ppcGjVBMgCNVSSTG/D1+V2jnLxooOrz7Bs0VEBgFGRPyBm4CV3i2rZflx3o+cN/I8YqJiAGhlnRA8b+R5OJ1OCksKKSguYGjvoQCMP20859193v7lLzjpAgCG9BpCUWnR/nvvjBkyhuCgYIKDghneZzhzls9hxqIZfDf7O7IvygagpLyEtVvWUlxazFnDziIkKASAM4ac4ZNtV0odnd1+m3ky/A7eK/uOLWUFSBkMT+vMXffdQdbQ4w6/AhvVJ3iuBv4DxANbge+A67xZlJ0OtWfia6HBoYefCf40QOi+1wdqN8Zwz2X38Lez//aHac+9+9zRF6qU8rnHwm/kP/mfA3DJyKE8+fpE2ib++U7EjdFhOxcYY3YbYy4yxrQ1xrQxxlxsjNnji+JaihF9RvDRDx+xp8Dzse4t3PuH6ZFhkURHRPPLwl8AeOubtxjaa+j+6R9M/QDw3NgtMiySyLBIwHN77IrKCvYU7GHa/Gn0zerLSQNPYuIXE/ePQL1111Z27d3FkF5D+OznzyivKKe4tJgvf/nS69utlDp6d5Q+wXVRp5EYGMX7P85g6sfvY9xuu8uql4Pu8YjI/cD/GWP2HmT6CCCk2d4C24e6pHfhvsvvY+jfhuJ0OsnOzP7TPJMenLS/c0FafBqvP/D6/mlBAUFkX5RNdU01Ex+YuL+9e0Z3hl8znN0Fu/n7FX8nLjaOuNg4Vm5YycC/DAQgLCSMtx95m16dejH2xLH0uKgHbaLb0DdLbzKrVGMWX9WJF6u+5M7AZXR39uWSW+/j7w89zZmDh3DPhGdok5Bqd4kHddCx2kRkDHAnnmtgFwB5QBCQAfQEvgf+YYxpUncoam5jtQ372zCeuukp+mT9cQilxnIdjo7VppR3GdwsD/6FO/1vYnLRYgDO7NKHT5fNbdD3acix2g56qM0Y87kxZjCeczzLASdQhGcQ7n7GmFsOFToiMlFEdonIsgNMu01EjIjEWK9FRJ4XkRwRWSIivWrNO15E1lqP8bXae4vIUmuZ5+VAd0JTSqlmplyKmBH2CU9FXc85kf2JD4imW/kwJhctJsI/kDO69ObWh+60u8xDqlevNmPMG7UbrMAoP8xybwAv4rlpXO1lE4FRwOZazSfj2ZPKAPoDLwH9RaQV8CDQBzDAfBH5whiTb81zJTAb+AYYDUymhZn232kHbH/oqod8WodSyjsMbmaGfs4HAW/wW/ViFpVuoabEcy4nLao1I7J7MvC4fgw4fSTdBwzDPzDI5ooPrz7BM1dErjTGzAIQkXOAfwKZh1rIGDNdRFIOMOlZPIfwPq/VNgZ40xoJe5aIRIlIezy32Z667zyTiEwFRovINCCiVk1vAmfSAoNHKdU8Gdz8LepkXin4DkqBUugVn8Ttp57PwFOG0H/4SU2mF1td9QmeC4GJ1j/2cUBrYMTRvJl13mirNfpB7UnxwJZar3OttkO15x6g/WDvexVwFUASSUdTulJK+VS+czuTCn/Y//rZW2/i5qefs6+gBlSf7tRLgcfxnOsZDlxvjMk99FJ/JiIhwL3AA0e67LEyxkwwxvQxxvSJpXHdAlYppQ6klSuelQFz6R3q+WP5tmeep3DPTpurahj1GST0NeBmoDtwOfCViBzNBaTpQCqwWEQ2AgnAAhFph+fC1MRa8yZYbYdqTzhAu1JKNRvtq9KZX+o5HX7liScRbo1u0tQdNniApcBwY8wGY8wUPCf/ex1mmT8xxiy1LkBNMcak4Dk81ssYswPPAKSXWr3bBgCFxpjtwBRglIhEi0g0nk4JU6xpRSIywOrNdil/PGfU5D337nOUVZQd8XJvfPkG2/K2eaEipZSvBZsI7oj2DI+1ZtMmVi2c2WQuEj2U+hxqew5IEpETrKYqPHtAhyQi7wEzgY4ikisiVxxi9m+A9UAO8ApwrfXee4FHgbnW45FaF7ReC7xqLbOOZtax4Ln3jzx4XC4Xb3ylwaNUc/JE/nvcFD2G6WtX06Xv8bQOCWNMlz68cN1trJz3a5MMooNeQLp/BpEr8ZyYb2WMSReRDOBlY8xIXxTY0BrjBaSl5aWcf8/55O7KxeVycd4J5/H4xMfpmNyRmKgYfnr5J6554hrmrphLeUU55448d/+tDFLOSGHsiWOZOnsqt154K1c/cTXxsfEEBwYzc+JMgoOCbdsu0AtIlWooWwKXMzX4Q37jN36qWMj6Cs8QW91i47j6/Iu56IFbiWzT1mvv35AXkNanV9t1QD8818tgjFkrIm0a4s0bo5uX38yiokUNus6eET15rstzB53+7cxviYuJ4+vnvgagsKSQ1798nZ9e/mn/iNWPX/M4rSJb4XK5GHntSJasXUL3jO4AtI5szYK3FwDw6uevHnAkA6VU05ZY2YW/VD7MX6zXGwIX8W3Ie7xa9j7X/b9/c8fLz3LxoBE8/eW7jfaW1/vU5xxPpTGmat8LEfHDczGnaiDd0rsxdc5U7nrhLn5Z+Mv+QT5r+/D7D+l1cS+yL85m+frlrNiwYv+0sSeO9WW5SqlGILWyJ9fk/4v5pZuYG/oN0c4QJvwyhX/edIPdpR1WffZ4fhaRe4FgETkRz7mVZjt08aH2TLwlMzmTBW8t4Jtfv+H+l+5nZN8/HsXcsHUDT739FHMnzSU6IprLHrqMisqK/dPre/sEpVTz1Kf0ZF6MeIyzqm7gmbc/5LGJbyGO+uxX2KM+ld2NZ4DQpcDf8HQEuN+bRbU02/K2ERIUwsWnXMwdl9zBgtULCA8Jp7jUc0fBotIiQoNDiQyLZOeenUyeefB+FOEh4RSXFR90ulKqeVkTNJu/Rp3IucU3Eez0519X39CoQwfqscdjjHHj6Wn2ivfLaZmW5izljufvwCEO/P38eenul5i5dCajbxxNXGwcP738E9mZ2XQ6rxOJbRIZ3H3wQdd12emXcfU/r240nQuUUt51lpzJioIdXD/6NP7+6vO0iW+8t0PY51C3RVjKIc7lGGO6e6sob2qMvdqaM+3VppR3nRs5gE8KZ/PCTbdz3TP/8trejq96tZ1m/dw3SsFb1s+L0c4FSinVKLxQ/haF4Wdxw3+eAidc//STdpd0WIe6H88mY8wm4ERjzJ3WyANLjTF34RlBQCmllM2CXeGc4XcSAPe88HyTuKC0PvtkIiKDa70YVM/llFJKeYHBzYywTxgfNZw4k8SN+c/QNz6Fd//zXKPvWAD16059BZ7bIkQCAuTD/muYlFJK+djx4R35tTiHcP9ALh06kr/+/UZ6Dz2pSYQO1K9X23yghxU8GGMKvV6VUkqpA1oc8gO/FucAnksxGvsoBQdy2OARkUDgHCAF8Nt3AzdjzCNerUwppdQfrAucT88yz3jN//fInU0ydKB+52o+x3Nr6hr234CVUm8W1RI9PvFxupzfhe4XdKfnhT2ZvWz2Ud8aQSnVvFRKKS9E38ZQM5Igpx9vPvkYV9/3T7vLOmr1OceTYIwZ7fVKWrCZS2by1YyvWPD2AgIDAtldsJuq6irGvj+Wi0+5mJCgELtLVErZ6D9Rd3BX/ksATHvvPYaOG2dzRcemPns8v4lIN69X0oJt372dmMgYAgMCAYiJiuHjHz5mW942hl89nOFXDwfgvSnv0W1cN7qO7cpdL9y1f/mwIWHc8swtdDm/CyOvGUlefp4t26GU8o5zyq4g0OHZT3ju8af5vzvuYuEv31FTXXWYJRun+tyPZwXQAdgAVOLp2Waa7cgFT98MaxY17Jtm9oTbnjvo5JKyEo7763GUVZZxQt8TGHviWIb2HkrKGSnMe3MeMVExbMvbxoDLBzD/rflEh0cz6oZR3Dj2Rs4cdibSV3j7kbe56OSLeOSVR9iVv4sX73yxYbfhKOnIBUo1jC2By3ky5F4+LfuZ3EpPH69QvwAGp2Yy5oRTGHPr5cR36OS192/IkQvqs8dzMpCB56LR0/GMaHB6Q7y58ggLCWP+W/OZcO8EYqNjGXvvWN748o0/zDN3xVyG9R5GbHQsfn5+XDT6IqYvnA6Aw+HYf2uEi0++mBmLZvh6E5RSXpZY2YXn8z9nc+VeNgYu5r3If3N52Cg25u7iupf+TUJGZ/7fTXfYXWa9HPQcj4js6y7RsoY6PsSeiTc5nU6G9R7GsN7D6NahG5O+mnTU69rX81Ap1fwIDpIru5Nc2Z1xeILmoejLeDh/Etc//xTXPuu98doayqGqmw/Ms37Wfcw7xHLqCK3euJq1m9fuf71ozSKS2yf/4dYI/br04+cFP7O7YDcul4v3przH0F5DAXC73Xz848cAvDvlXY7reZzvN0IpZZsH8ycS7PAHYMqEN+wtph4OusdjjGn8Y2s3EyXlJdzw5A0UlBTg5/SjQ0IHJtw3gfemvPeHWyM8cf0TDL96OMYYTj3uVMYMHQN4bgQ3Z/kcHnvtMdq0asMH//jA5i1SSvmKi2puiT6X8vxqAF797wRGX924B5c5bOeC5qY53hYhbEgYJdNL7C7jgLRzgVINq4YqNgUuJSdgCQv8ZvFK+SdsqNgDwCcT/sOpF19BoBfuSuyr2yIopZSyicHN5IjXWOZcyHqzgfWuXNbX7GBTxV5qKt2ePsZA15g4/u/6O/nbE//A4XTaW3Q9afA0A411b0cpdWTy/DcyNeR9ZsgvfF85l7VFnmvyYoJDSWsVS7/2nRibnESHLhl06NeFrN4DaN0u0eaqj1x9xmp7GphojFnug3qUUqpF+TLiJd6UScyqWOW5PqcQwv0D6ZeUzo2nj+fSv99JRKtYu8tsUPXZ41kJTBARP+B14L1mN0K1GwwGQbshNySDgcZ/TyqlbLEqaCaTgp/nifz3ARjXZzB9emfT75QhDDrlTJx+/jZX6D31uS3Cq8CrItIRuBxYIiK/Aq8YY37ydoG+EJQTxJ5We2jt11rDp4EYDHtq9hCUE2R3KUo1GlWU81HUczxfM4E5JRtxVjo4pVMPnnj5OboNHWZ3eT5Tr3M8IuIEOlmP3cBi4FYR+ZsxpmmPVgckPJRA7kO55HXI03urNhS3J9ATHkqwuxKlbLXDfx0zQ77hJ8f3fFj6IzsLSmgTHMZT117PRXfcTLuUdLtL9Ln6nON5Fs8wOT8C/zDGzLEm/UtEVnuzOF/xz/cn9Sa9bEkp1TDKpJD3Ip/mherXWVyaC4UQ7PRnZGZXrr7rak4cewkBQcF2l2mb+uzxLAHuN8Yc6B48/Rq4HqWUavIujTiJTwpm0zUmjqcuu56Bp4+gz9DRLTpsaqvPOZ7XRSRaRLoAQbXapze7TgZKKXWM5oZ+zSeFswFYsnNLox83zQ71OdT2V+AmIAFYBAwAZgIjvFqZUko1Qa1r4vY/19A5sPp8KjcBfYFNxpjhQDZQ4M2ilFKqKVoVNJN7gq4DoHdKks3VNF71CZ4KY0wFgIgEGmNWAR29W5ZSSjUdJY69nBHZm84Vg/i8eC5XDT2JaYsW2l1Wo1WfzgW5IhIFfAZMFZF8YJM3i1JKqaak0G8XXxYuAOD6s8cw8pzRFO7ZRWh4lB5uO4AjGp1aRIYCkcC3xpgmebPvA41OrZRSx+p/Ec/zYM2/WFa2bX9bq8AQusUl0r1jR7r3705G7y4kZ2YRn5KJf2DTuri6IUenPmTwWBeOLjfGeO9G3j6mwaOU8qZ85zaWBf/KEr85LGEpS2rWsrQ8l1LX73+rO0SID4kkKSKG5Nh2JCckktwhhaSuaST37ERyxy6ERkTbuBV/5rPgsd7sc+AGY8zmhnhDu2nwKKV8zY2LjYFLWBewhE3OHDbLJjaZLWyq2c6m6jxyqwqoMX8c2LB1QAjJYTEkxbYhPT2Z408ZzvALzrdtwFBfB890PD3Z5gD7LyI1xpzREAX4mgaPUqqxcVHN9oAcNgWsYJNzLZscG9hkNrPZvZ1N1TtZX7mbCncNgQ4/Zk79muwRo3xeo6+DZ+iB2o0xPzdEAb6mwaOUamoqpZTPI15mbOHtAFRXVeLnH+DTGnx2B1LrHM9/m9M5HqWUamr8TAB5smP/681rl5OWlW1jRcfmkP38jDEuYLWI6JVQSillk/GRJ3B9wVMM69CJeT9PbtKhA/W7jicaWC4izeIcj1JKNTVTyuYD8P3KJc3iBnH1ubLp73hui/AI8HSth1JKKR+4JuxsACY9+JjNlTSM+oxO/YdOBCJyHHAB0CQ7FyilVFOxPWAtj4fcyhtF3wEQFNg8bqtQ3zuQZgMXAucBG4BPvFmUUkq1VBsCF/F2yAtMq5nJ9JLV1BS4uaD3cdz77GN0Pf6AnYybnIMGj4hk4tmzuQDP7a4/wNP9eriPalNKqRblhejbuDH/GaiEQKcffx16EuP+dilDx42zu7QGdag9nlXAL8BpxpgcABG5xSdVKaVUC+Q0v/+TXOmq4c0Z37N8/TqGvfcNw848iUFjxxAUEmZjhQ3joBeQisiZwDhgMPAt8D7wqjEm1WfVeYFeQKqUasxcVLMy+DcWBPzCAsc8fqtazPzSTbgxhPj5MzAxg77dunPbS/8iJs53V7r4euSCUGAMnkNuI4A3gU+NMd81RAG+psGjlGpK9vhtYX7wD0x0vMIHhb/tb3/h+tu5/oUnfVaHz0YuADDGlALvAu+KSDSeDgZ3AU0yeJRSqrEwuCl27GFr4Bq2+q1nm2MTGx3rWetex9qaLayt3MnemjIo9szv73DSuXU7zj/hZP7yj/vtLf4Y1KtX2z7GmHxggvVQSil1EFWUsz1wLdv81rPVbwNbZTPbZBtb3TvYanaxrWYvWysLPLdLKP/jsomhUWS2iuP8+KFkZHYgo19nMgb1JD0ru8ndx+dAjih4joSITMRz4ekuY0zXWu03ANcBLuBrY8ydVvs9wBVW+43GmClW+2jgP4ATzzmmJ6z2VDznnVoD84FLmurN6ZRSTVellJITNJ8lATNZ4JjHgprlLKvcwq6qEqjE87AEOJ3EhUcSH92Knq3TOaVNW+IT44lLSyA+K4X41A7Ep2QSHBZh2/b4gteCB3gDeBHPOSEARGQ4nvNFPYwxlSLSxmrPwtORoQsQB3xvdecG+H/AiUAuMFdEvjDGrAD+BTxrjHlfRF7GE1oveXF7lFKKfOc2/hNxL7NcC1hTtY1NFXtxlxsoh0CHH93bJnB69gASMxKI75hAfHIycUmpxCdn0Lptgt4KGy8GjzFmuoik1Gm+BnjCGFNpzbPLah8DvG+1bxCRHKCfNS3HGLMeQETeB8aIyEo8HR0utOaZBDyEBo9Syos2BS4hpbIH5EN2u0QGdOrMpZnpZPboRJfBfcnqNbhZHArzNm/u8RxIJnC8iDwOVAC3G2PmAvHArFrz5VptAFvqtPfHc3itwBhTc4D5/0RErgKuAkhCB9pWStVPqSOfFyLuYYaZzfzydeyo9JzlH5zWgRnr1tpcXdPl6+DxA1oBA4C+wIcikubtNzXG7O8Q0Uf6HLr/uFJKWV6MvJd78v+Lv9PJBQOOp3fvbEZdeh4ds/vbXVqT5uvgyQX+ZzwXD80RETcQA2wFEmvNl2C1cZD2PUCUiPhZez2151dKqQZxesVF3M3LVLtcXHP3LQw4Te8G0xB8fZbrM2A47B8LLgDPOHBfAONEJNDqrZYBzAHmAhkikioiAXg6IHxhBddPwLnWescDn/tyQ5RSzV9W+XHE+ocCMPD0MQyK68C1Q0/ly5dfoqZaO9EeLa8Fj4i8B8wEOopIrohcAUwE0kRkGZ6u0OONx3LgQ2AFnuF5rjPGuKy9meuBKcBK4ENrXvBcxHqr1RGhNfCat7ZFKdVybarZxttR/+DG6DPwL3Hyzm8/cMY115LVNoG/X3IJi2Z+j3G77S6zSTnskDnNjQ6Zo5Q6FtVU8EnkC7zsfo1fitfgxpAW1ZoTu/elz3F96HvaMLr3H97suk37dKy25kaDRynVUHb5b+Dz0Nf41P0Vv5WtprCmAoBRGV2Z+NVHxGd2srnChtOQwdO8IlkppXyoTXUqVxY8xjdFi9hbU8K34RNpFxDOd2uX8d4/X7S7vEbL173alFKqWTC4WRzyI78GTOEX8xu/VCxjW3ERAN1j4zn/9qttrrDx0uBRSqmjcH30Gfxf/tdQBvGhkQzp2Y3jhw5m+Pln0Cl7YLM7x9OQNHiUUuoohOO5E2h8ZCS3XHQJnQd2p9uAQSSkddbQOQztXKCUUkfBRTWPRl/Jf0v/x46q4v3tbYLCGJzRmVGnj+TkSy8guWN3G6tsONqr7Rho8CilGtpuv82sCprLYr9ZzGU+P1UsYnNFPgCndO7BrfffxvCxF+JwOm2u9Ohp8BwDDR6llLcZ3KwJmsMHwRN4pvg9CmsqiA4MZmSn7jz93qskde56+JU0Mho8x0CDRynlSxVSwv8iXuRt8z6TixYDkLNoHuk9ettc2ZHR4DkGGjxKKTtUU0FqQDu2VhWSGtma9QW77S7piOgFpEop1USsDZrDLdFn0Tkoka1VhQD84767bK7KXtqdWimlvGSn/3oyK/pDBYzu2J07x1zFhffeRlhkK7tLs5UGj1JKeUmyqyMAj13xV+579RWbq2k89FCbUkp5gRsXle4aAE4bP87mahoXDR6llGog1VTwStT9nBc1gDb+Efvba8qrbayq8dFDbUop1UAuiRzJBwW/0S40nNMG9mfESSM4+bKLiY1Lsbu0RkWDRymlGoDBzZqaLQCc1nsAr/z8nc0VNV56qE0ppRqA4ODigDMByN27095iGjnd41FKqaNkcLMlcAXLA2ex0rmE/1VPBmDmmjU2V9a4afAopdRReDr6Bv5d+ga7Kkug0tMWExTKxX2HcOfTD9tbXCOnwaOUUkeo0LmT2/M9t7Z+8obr6D9qKJ179yemfZLNlTUNGjxKKXWEbg+/BArgrO59ueWZZ3H6+dtdUpOiwaOUUvVUKaW8E/lvJhX+iL/DyUfzf9XQOQoaPEopVQ8GNyeEZTOjYC0dW7fhn3ffr6FzlDR4lFLqMEod+ZwVOoQZxWv51zXXcMeLLyIOvRrlaOknp5RSh7EucBFTi5cB8Mp7H3HV8aN544FH2L19s82VNU0aPEopdRjdy4ezLnABz0XdTIY7no/n/MLljz5IXHwad48bx/qVC+0usUnRO5AqpdQRcuNiccgPPBfwT94smAZAXEgEvZM70KdHD44/5wSGnjUWh9Npb6ENSG99fQw0eJRSDWlt0By+Df6A2e55zK9azerynRggLaI1fzn7LMbfdysJHTrbXeYx0+A5Bho8SilvKnLk8VX4a7zinsS04lU4xcFFfYbwr/dfpV1aut3lHbWGDB49x6OUUg0owh3LhYV382Pxcj6P+D9cxs2bc6fxwu0P2F1ao6HdqZVSqgHUUMWq4JksDJjBLMdMvi77jU1F+QCclpXN1f+41+YKGw8NHqWUOgouqpkbOpnvAj5jWs0sZpXmUF5eDeUQ4ufPyNRu3Df6Wk699VLiUjLtLrdR0eBRSqkj4MbFU9E38O+SSewpLUNKIbtdElcddwp9TuxHryFDyOzeDz//ALtLbbQ0eJRSqp6+DZ/IDTV3kZO/m0Fp6Vx/05WMOv9CWrdLtLu0JkWDRymlDiEncB5vhbzAR5XfsbJ4B51bt+Pdx55g3M136LA5R0mDRymlDqCaCm6MOpuXCybjqBKGpHbk+nMvZfx9dxMaEW13eU2aBo9SSh1AXsBmXi7w3Mo6o3Usl194EZc+er/NVTUPup+olFIHEFeVycqg37gq6iRW797F+Mf+zo7N6+wuq1nQ4FFKqQPI89/I8sDZTCiYAsDfThhNm/gUe4tqJvRQm1Kqxct3buPz8FdZLStZ697AiupNrCzbAYWe6beeeTZPf/qJvUU2Ixo8SqkWbU3QbDpWDIACz+uOrdqQGZ/AJd3OYOj5o+lzymgCgoJtrbG50eBRSrVobaqTSQqMZnNlPmf16MvH82c2q9sZNEZ6jkcp1aJFudqxsXI3UX5BfLp4LjedeSYzv/6MsuJCu0trtjR4lFItnuBgWsBX9AlN5sWvvmLQaWcRGhHFd6+/YXdpzZIealNKtUh7/LawNnAha/yWsNaxmrXu9bhq3H+YJzUry6bqmjcNHqVUs1UhJawPWsRa/8Usdy5mpXs1a11bWFO5g/zqcqjxzOcQITWiNRkJcdyU0ose/XpyxjVX0Lq9jsHmDRo8SqlmZX3gQj4NmcinNd8ys3gd7nID5Z5piWFRZLaLY2zyMDI6dyCzbxcyuvcgtVMP7bnmQxo8Sqlm48Ho8TyS/yZUQo92Cdx90gV0GdiNDtld6dxzAOHRMXaXqNDgUUo1I37WP2mxIaF889MPxHXSG7A1RtqrTSnVbPw9/zU+iHyKvLJSLh59pt3lqIPQ4FFKNSvbZAsAa/K32VyJOhivBY+ITBSRXSKyrFZbTxGZJSKLRGSeiPSz2kVEnheRHBFZIiK9ai0zXkTWWo/xtdp7i8hSa5nnRUS8tS1KqcbN4GZN0Gxejfo737h+AODcAcPsLUodlDfP8bwBvAi8Wavt38DDxpjJInKK9XoYcDKQYT36Ay8B/UWkFfAg0AcwwHwR+cIYk2/NcyUwG/gGGA1M9uL2KKUakSrKWRs8lylBH/N21WcsLN0CFRAdEMxpnbO59bnH7C5RHYTXgscYM11EUuo2AxHW80hg377wGOBNY4wBZolIlIi0xxNKU40xewFEZCowWkSmARHGmFlW+5vAmWjwKNWsFDh3sD5wMev9V7DBkcN6s4Ec12ZyqrezuSJ/f1fpXu2T+M9VtzD6srFkdO2rt6Ru5Hzdq+1mYIqIPIXnMN8gqz0e2FJrvlyr7VDtuQdoPyARuQq4CiCJpGPaAKWUd9RQxdKQabwV8DI/V89nfeUuCmoqoOz3eVoFhtChdRsGxWVxaWoaHbpnMujM0aR37W1f4eqI+Tp4rgFuMcZ8IiLnA68BJ3j7TY0xE4AJAH2kj/H2+yml6qdcipge9glXVN7C1qpCKAMpg6Hpnbiww0jSOqeS1iOT1I6dSc3sRmTrtnaXrBqAr4NnPHCT9fwj4FXr+Vag9tgUCVbbVjyH22q3T7PaEw4wv1KqCSh07uSG8LF8XPQb5cXVAAT7+/PaPx/l+DPPICG9s80VKm/ydfBsA4biCY8RwFqr/QvgehF5H0/ngkJjzHYRmQL8Q0SirflGAfcYY/aKSJGIDMDTueBS4AUfbodS6jCKHHnkBq4m1y+HXOcGtkoum8xmllWvY17ZJlwFbq4acRJnXnYux512po4q0IJ4LXhE5D08eysxIpKLp3falcB/RMQPqMA674KnV9opQA6eI7qXA1gB8ygw15rvkX0dDYBr8fScC8bTqUA7FijlZcWO3WwJXMUuZy55zm3kyU7yZBd5Zjd5Zg957gJ2uPLJrcyn2FW5f4y0fdoEh5HVJoFrM07jjEvP5YRLLrFnQ5StxNORrOXoI33MPObZXYZSjVKhcycbA5ex0W8Vm5zr2Gg2stHksqlmJxurdrO3uuyAy0UHBhMbEk5seARtI1uR2LY9CQkJJHRIJqFrGvEd04lLziAwONTHW6QaiojMN8b0aYh16VhtSrUwFVLCjLBPWeo3j41sZJPZxsaaHWys2E1hnV5kIX7+JEe2JqV9W/rHdyU5NYmkjim0zUggtl08se0Tad02Hv/AIPs2SDU5GjxKtQDbAtbwRcjrfOGezE8ly6ko9tyIJsw/gJTWMaTEteX45J4kd0gmpVM6yWmZpGR2IaZdkl4ToxqcBo9SzVAV5SwIncoPAV/yedX3zC3dCFXQITKGv40+lVEXnEr/E0bTqk28BovyOQ0epZq4GqpYGfwbiwJ+ZaHMZ75rOXNK11NRWgOl0L99Kv845SrGXHsJnYcM0qBRttPgUaoJKXXkszT4FxYG/MpCs4iFNatZWraVyvIaKIcgpx/d2yZw9cDTOP704Rx31um0iU+1u2yl/kCDR6lGqEJKyAmaz+qAhaxyLGeZeyWLqnJYU74Ld6mBUk9PsuzEFK7vMobsob3JHnI8md374ecfYHf5Sh2SBo9SNqmQEtYHLSLHfyk5zpXkmHX7B8DcVLF3/wCYAElh0WSnpTO222iyT+hL9pChJKZn6WEz1SRp8CjlRaWOfNYF7guXVVa4bCGneju5FQWYcvaHS1RAEBmt29E/rhOXpKfTqVcWHQf2JLNHX8IiW9m6HUo1JA0epRpAuRQxK+wrZvn9zNp94VK1g+2VRZ5gscIlNjiUDm3aMSy+Bx06eEZX7tCrKx2ysmnV9qADrCvVrGjwKHWUKqSELyMm8KZ5m+9KllBV7AKgfWgEHdq2Y3RSXzp0SqNDdkc6ZHUlvXNPHV1ZKTR4lDoqr0b9nTtLnia/sJy40AiuP2UMI8aexODRZxAV087u8pRq1DR4lDpCr0U9wJUFjzEiozN3PXYnI8++CKefv91lKdVkaPAodYReqXmbrm3i+G7FYg0cpY6C9sVU6ggl+7UjZ88u/nHV1WzJWW53OUo1ORo8Sh2h58omcWJYdx54fSJJGV3JimnPLaedxeS3XqU4f7fd5SnV6On9eJQ6SmuD5vBF8Ft855rG9JJVVLhrcIiQ3T6R43v1YshpIzju9NOJjUuxu1SljllD3o9Hg0epBlAuRcwI+5Rf/H5ges0cZpeuo8LtufVA5+h2HN+lB0NGDWXIuDNIzOhic7VKHTkNnmOgwaN8oVJKmRf6Lb/4TWW6exa/lq2iqKYSgOTwaIZ07s7Qk4Zw+t8u10E8VZOgwXMMNHiUHVxUsyTkJ34J/JZf3L8xvXw5u6pKcIqD4YmdGTv2LC64/1ZCI6LtLlWpA9LgOQYaPKoxMLhZEvITHwZO5IPy71hXsZvWgSFcd/a53PD0o8S0T7K7RKX+oCGDR3u1KWUDwUGPspE8nv8Oayt2MiPsY44LzOKR996kQ3JHXrznDlw11XaXqZRX6B6PUo3I8uDp3Ox/Ld8XLSctsjVnDDqewacNYdAppxKXkml3eaoF00Ntx0CDRzV2BjefRrzIf3llfzdtgOTIaAZldWHQiAEMGn0C3fsP15u+KZ/R4DkGGjyqKaminEWhP/Cb//f8Zmbza8UKtlUWARDi50/3tolkZ3Qiu1822aMH07X/8QSFhNlctWqONHiOgQaPasoMbrYEruC3oMnMcvzKwpqVLCrftL+rtlMcdG7djsFZ3TlhzCiGX3gurdsl2ly1ag40eI6BBo9qbty42Bi4hIVBP7PIMZ/5rqXMKF1FsasSP3Ewtsdx3P7Pv9Nz9Al2l6qaMA2eY6DBo1qCaiqYF/otHwZM4tWib6h0u1i1fC5pnbPtLk01UdqdWil1SP4EMbD0TJ7N/5SZAd9RbVy8+c//2F2WUoAGj1LNXpfy4zkzsg9Pvvsu2zettbscpTR4lGruBAcPVD9OldtFj07ZbFu3xu6SVAunwaNUM7cieAYfBb5BjXGTV1HK9I8+t7sk1cLpra+VaoaqKGdm2OdMcP4/3iucgaPCwbDUTow77yzOueUGu8tTLZwGj1LNgMHNyuDfmBr4CVPNNKaVrKC0pIoQpz93njOW255/Qm9IpxoNDR6lmiiDm5/DPuQtv1eZUj6PreWFUA4ZrWMZP+pETjx3NMPHnEdk67Z2l6rUH2jwKNXElDj28t/IB3m5/H1ySnYT4R/ISb2yOfGskZx45vmkdOxud4lKHZIGj1JNyG6/zZwYOJhF+bkcl5jBA1fczDm3Xk9IeKTdpSlVbxo8SjURi0N+4Bz3OLaWF/DVK//l1L9eZXdJSh0V7U6tVCO3x28Lt0WfQ7+KkyhzVPHDZx9q6KgmTYNHqUbs08gXSCOD5wo+5cLBQ1m8ejGDTj/L7rKUOiZ6qE2pRuxp9/MU1VSydPY0uvYbanc5SjUI3eNRqhE73+90AG695HpmfP4Jxu22uSKljp0Gj1KN2A35T/Gf6FtYuGEDx595Lv0T03n3ySeoqii3uzSljpoGj1KNmODgxvxn2OjawotRt1OYX8ZFd95Daut2/GP8VeTv2mZ3iUodMQ0epZqAUHc01xU8ycrybXwdMYGuziTue/MVMpIyeOm2u3DVVNtdolL1psGjVBPiwMkpRVcypXgpC0Om0i0giWuf+TcDkzPZvW2z3eUpVS8aPEo1UT3LTuDH4uW8F/lvlu7I5eYzL7G7JKXqRYNHqSZMcDCu8A5GhfVg1soVdpejVL1o8CjVxP0Y/g7fly4hMy7e7lKUqhcNHqWaKBfVPBF1NScWX0JKVGsmfvuJ3SUpVS8aPEo1QVsDVnFieA/uKfgv5/YZwKy1S2mXmm53WUrViwaPUk3MVxH/pYfpzeyyHF578AHenz2D8OgYu8tSqt68FjwikigiP4nIChFZLiI3We2tRGSqiKy1fkZb7SIiz4tIjogsEZFetdY13pp/rYiMr9XeW0SWWss8LyLire1Rym7VVHBD1OmcXnQ1CRGtmD/zZ/7y0MOIQ/9+VE2LN7+xNcBtxpgsYABwnYhkAXcDPxhjMoAfrNcAJwMZ1uMq4CXwBBXwINAf6Ac8uC+srHmurLXcaC9uj1K2uibqdF4s+IqbTh3DrM0r6dR3oN0lKXVUvBY8xpjtxpgF1vNiYCUQD4wBJlmzTQLOtJ6PAd40HrOAKBFpD5wETDXG7DXG5ANTgdHWtAhjzCxjjAHerLUupZqVDyOf5rWC77n33At57qvPCAoJs7skpY6aT/bRRSQFyAZmA22NMdutSTuAttbzeGBLrcVyrbZDteceoP1A73+ViMwTkXl55B3bxijlY1sCl/O30vvpH5fKQ2+/Znc5Sh0zrwePiIQBnwA3G2OKak+z9lSMt2swxkwwxvQxxvSJJdbbb6dUg6ihiteiHqCPGUANbt7+8EP8A4PsLkupY+bV4BERfzyh844x5n9W807rMBnWz11W+1YgsdbiCVbbodoTDtCuVJNWTQWvRz1Ep6B4/lrwKOmxbZk++Us6DO5jd2lKNQhv9moT4DVgpTHmmVqTvgD29UwbD3xeq/1Sq3fbAKDQOiQ3BRglItFWp4JRwBRrWpGIDLDe69Ja61KqSVoePJ3uISn8peBhIqND+PSlF/h18xqyTxhld2lKNRhv3vp6MHAJsFREFllt9wJPAB+KyBXAJuB8a9o3wClADlAGXA5gjNkrIo8Cc635HjHG7LWeXwu8AQQDk62HUk3SxoAlDKg6kdDAAD6f8H+cfsXftKu0apa8FjzGmBnAwa6rGXmA+Q1w3UHWNRGYeID2eUDXYyhTqUbjx5BPKCmo4tV/P84ZV15jdzlKeY3+OaVUI3F66eUkBkZx4+0Ps3buHLvLUcprNHiUaiRiq1OYKpNxuwwnDB/Njg3r7C5JKa/Q4FGqEelYMYDrws5lc2k+37/7gd3lKOUV3uxcoJQ6AsWO3dwRcQn/zf+WrjFxjLn2KrtLUsordI9HKZsZ3Lwf+SRd/TOZUPAtd4w+j7kbV+uI06rZ0uBRykaLQ36gf2g6FxTeSXRkKDM++x//nvwhQaE6FptqvvRQm1I2WRTyPcOqTic4wJ83HnqYi++7B6efv91lKeV1GjxK2aCaCsa4zic8KJBf588gKVMvR1Mthx5qU8oG64MWsbkyn2svGauho1ocDR6lbJBe0Yvs0ESee+Ndqisr7C5HKZ/S4FHKBn4EcEnAWewqL2H1ktl2l6OUT+k5HqV8rFyKeCPqcR4unkBW63Zkdutnd0lK+ZTu8SjlI5sCl3BX9FgS/Npzbf6/6RGfxFc/TSYgKNju0pTyKQ0epbxsYch3nBPZn7TKnjxV8BHDs7rx8/8+YNr6laR262l3eUr5nB5qU8oLDG6+C5/Es/IsU4qWEh0QzJ1jx3LNo/eRlKG92FTLpsGjVAMqlyLejvw3z1W9xoriHbQLCefR8Vdw/VOPERXTzu7ylGoUNHiUagArg3/lzaAXebXkS3YXlNKzTQKT7n2UsbfeQmBwqN3lKdWoaPAodZQKnTv5IPx5Xqt5lzklG3FWODg1K5tb/n47Q887X29brdRBaPAodQQMbqaHfchrfi/zcdFvlBdU0yWmPU+Nv56L7ryZdknpdpeoVKOnwaNUPVRICZMi/8FzlRNZVbKTCP9ALh02kivuvp4+I0/WvRuljoAGj1KH8U7UP7m97J/sKCimT/tkJt33KOfeeAMh4ZF2l6ZUk6TBo9RBlDj2cl3kObyZP42BSem8+/TLDDt7nO7dKHWMNHiUOoANgYsY7TiJnII8HrxqPPe/OAE//wC7y1KqWdDgUeoArgi4kF2VRfzwv3cZduY4u8tRqlnRYwZK1VHkyOOn4pXcPG6cho5SXqDBo1QduYGrAcjs28nmSpRqnjR4lKqj0LEbgOhWsTZXolTzpMGjVB0BBAFQVV5pcyVKNU8aPErVIUYAMG63zZUo1Txp8ChVR41UA+AU7fSplDdo8ChVR5EjH4DI+FY2V6JU86TBo1Qd1eI5t+Pn729zJUo1Txo8StVRJiUAhOpYbEp5hQaPUnVUSDkAgYHBNleiVPOkwaNUHZVUABAYHGJzJUo1Txo8StXhxtON2qGjUCvlFfqbpVQdUcbTm61gb57NlSjVPGnwKFVHiAkDoLy0xOZKlGqeNHiUqqOaKgD8A/T+O0p5gwaPUnWUSxkAIaERNleiVPOkwaNUHfuu4wkOCbW5EqWaJw0epeqolH3dqTV4lPKGFjgK4mpgmN1FqEbs9OJ1dAZaj7sYAgLtLkepZqflBY8DrNutKHVA4nJADWD0tghKeUPLC57sjjBvmt1VqEYsfs8qwua/gf+gmyCsvd3lKNU4iDTYqvQcj1J1OQOIDAsCd5XdlSjVLGnwKFWX23MjOBx6WwSlvEGDR6m6XJ778eDQC0iV8gYNHqXq2rfH49TgUcobNHiUqsvt8vx0OO2tQ6lmSoNHqT8x1k/99VDKG/Q3S6m69l2/04DdR5VSv/Na8IhIooj8JCIrRGS5iNxktT8pIqtEZImIfCoiUbWWuUdEckRktYicVKt9tNWWIyJ312pPFZHZVvsHIqIH5dWxM/v2eDR4lPIGb+7x1AC3GWOygAHAdSKSBUwFuhpjugNrgHsArGnjgC7AaOD/RMQpIk7g/wEnA1nABda8AP8CnjXGdADygSu8uD2qxbCCR/SAgFLe4LXfLGPMdmPMAut5MbASiDfGfGeMqbFmmwUkWM/HAO8bYyqNMRuAHKCf9cgxxqw3xlQB7wNjRESAEcDH1vKTgDO9tT2qBdFDbUp5lU/+pBORFCAbmF1n0l+AydbzeGBLrWm5VtvB2lsDBbVCbF/7gd7/KhGZJyLz8vL0dsbqcPRQm1Le5PXgEZEw4BPgZmNMUa32+/AcjnvH2zUYYyYYY/oYY/rExsZ6++1UU7fvHI/u8SjlFV4dJFRE/PGEzjvGmP/Var8MOA0Yacz+M7lbgcRaiydYbRykfQ8QJSJ+1l5P7fmVOga6x6OUN3mzV5sArwErjTHP1GofDdwJnGGMKau1yBfAOBEJFJFUIAOYA8wFMqwebAF4OiB8YQXWT8C51vLjgc+9tT2qBTFuQHSPRykv8eYez2DgEmCpiCyy2u4FngcCgamebGKWMeZqY8xyEfkQWIHnENx1xhgXgIhcD0wBnMBEY8xya313Ae+LyGPAQjxBp9SxcQYCBlxVOmyOUl7gteAxxszgwMcqvjnEMo8Djx+g/ZsDLWeMWY+n15tSDccv2POzukyDRykv0AsVlKpr3+0Q3DWHnk8pdVQ0eJSqa//9eFreDXqV8gUNHqXqqi71/PQPtbcOpZopDR6l6nJVeYbLceodSJXyBg0epepyV+thNqW8SINHqT8x6K+GUt6jv11K1eXwB3dVrdsjKKUakgaPUnU5gzyjF2h3aqW8QoNHqbr2nd/R4FHKKzR4lKpr32gF7ip761CqmdLgUaqu/TeCc9pbh1LNlAaPUnXpyAVKeZUGj1J17Tu349ALSJXyBg0epepyV+O5H4/+eijlDfqbpVRd7mrPcDl6IzilvEKDR6m63C70V0Mp79HfLqXq8gsGV+XvvduUUg1Kg0epuhx+gNHgUcpLNHiUqmtfpwIdq00pr9DgUaou4/L81M4FSnmFBo9SddVUeA636QWkSnmFBo9SdRkXiIaOUt6iwaNUXeL4/XCbUqrBafAoVZczyHMRqVvDRylvENPCeu6ISDGw2u466iEG2G13EYfRFGoErbOhaZ0Nq6nU2dEYE94QK2qJB7JXG2P62F3E4YjIvMZeZ1OoEbTOhqZ1NqymVGdDrUsPtSmllPIpDR6llFI+1RKDZ4LdBdRTU6izKdQIWmdD0zobVours8V1LlBKKWWvlrjHo5RSykYaPEoppXyqyQePiEwUkV0isqxW20MislVEFlmPU2pNu0dEckRktYicVKt9tNWWIyJ3+6jOD2rVuFFEFlntKSJSXmvay7WW6S0iS606nxdp2JEsRSRRRH4SkRUislxEbrLaW4nIVBFZa/2MttrFqiNHRJaISK9a6xpvzb9WRMb7oMYnRWSVVcenIhJltdvyeR6izkb1/TxEnY3q+ykiQSIyR0QWW3U+bLWnishs6z0/EJEAqz3Qep1jTU+pta4Dfs5ervMd6/2WieffA3+rfZiIFNb6PB+otS5v/n8/WJ1viMiGWvX0tNob7nfdGNOkH8AQoBewrFbbQ8DtB5g3C1gMBAKpwDrAaT3WAWlAgDVPlrfrrDP9aeAB63nKIeabAwwABJgMnNzAdbYHelnPw4E11uf2b+Buq/1u4F/W81OsOsSqa7bV3gpYb/2Mtp5He7nGUYCf1f6vWjXa8nkeos5G9f08WJ2N7ftprTPMeu4PzLbe60NgnNX+MnCN9fxa4GXr+Tjgg0N9zj6o8xRrmgDv1apzGPDVAdbj7f/vB6vzDeDcA8zfYL/rTX6PxxgzHdhbz9nHAO8bYyqNMRuAHKCf9cgxxqw3xlQB71vz+qRO66/C8/F8GQ9KRNoDEcaYWcbzf/xN4MwGrnO7MWaB9bwYWAnE4/k8JlmzTar1vmOAN43HLCDKqvMkYKoxZq8xJh+YCoz2Zo3GmO+MMTXWbLOAhEOtx9uf5yE+y4Ox5ft5uDoby/fT+o6VWC/9rYcBRgAfW+11v5v7vrMfAyOtbTnY5+zVOo0x31jTDJ6APuT3E+//fz/Y53kwDfa73uSD5xCut3YHJ4p1WAjPL9OWWvPkWm0Ha/eV44Gdxpi1tdpSRWShiPwsIsdbbfFWbft4tU7r0EQ2nr+E2hpjtluTdgBta9Vk22dap8ba/oLnr7N9bP08D1Bno/x+HuTzbDTfTxFxWof8duH5B24dUFDrD47a77n/c7OmFwKt8cHnWbdOY8zsWtP8gUuAb2stMtA65DVZRLrUrd+GOh+3vp/PikjgYeo54jqba/C8BKQDPYHteA4TNGYX8Me/JrcDScaYbOBW4F0RifBlQSISBnwC3GyMKao9zfqLzfZ++AerUUTuA2qAd6wmWz/PA9TZKL+fh/h/3mi+n8YYlzGmJ569hX5AJ1+875GqW6eIdK01+f+A6caYX6zXC4BkY0wP4AXgM5vrvAfP59oXz+Gzuxr6fZtl8BhjdlofqBt4hd93o7cCibVmTbDaDtbudSLiB5wNfLCvzToEsMd6Ph/PX3WZVk21d8+9Uqf1F9knwDvGmP9ZzTut3ep9h1R2We22fKYHqRERuQw4DbjICkhbP88D1dkYv5+H+Dwb3ffTet8C4CdgIJ5DPvvGnaz9nvs/N2t6JLAHH/6+16pztFXHg0AsnsDeN0/RvkNexphvAH8RibGrTuvQqzHGVAKv443vp2mgE1V2PqhzshNoX+v5LXiO5wJ04Y8nFdfjOYHnZz1P5feTeF28XafVNhr4uU5bLNbJTjwnFrcCrazXdU/entLANQqeY/PP1Wl/kj92Lvi39fxU/njCcY75/YTjBjwnG6Ot5628XONoYAUQ2xg+z0PU2ai+nwers7F9P633jbKeBwO/4Pkj4yP+2LngWuv5dfyxc8GHh/qcfVDnX4HfgOA687fj94v5+wGbrc/P2//fD1Zn+1rfi+eAJ6zXDfa73iAbYOcDzyGA7UA1nmOLVwBvAUuBJcAXdX7R78PzF9pqavW4wdNjY4017T5f1Gm1vwFcXWfec4DlwCI8u+Gn15rWB1hm1fnivi9sA9Z5HJ7DaEus919kfTatgR+AtcD3/P4PjQD/z6pnKdCn1rr+gufEbQ5wuQ9qzMFzrHlf275/dGz5PA9RZ6P6fh6szsb2/QS6AwutOpfxey+7NDyBl4MnhAKt9iDrdY41Pe1wn7OX66yx3nPfZ7yv/Xrr81yMp1PMIB/9fz9YnT9a389lwNv83vOtwX7XdcgcpZRSPtUsz/EopZRqvDR4lFJK+ZQGj1JKKZ/S4FFKKeVTGjxKKaV8SoNHKaWUT2nwKKWU8qn/D8JhjOIAzOxUAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Plot up the raw stuff\n", "fig, ax = plt.subplots(figsize=(6,10))\n", "print(df.columns)\n", "cropped = df.iloc[start:stop]\n", "\n", "ax.plot(df['Sensor3'], df.index, color='darkorange', alpha=0.5, label='raw signal')\n", "ax.fill_betweenx(df.index[start:stop], df['Sensor3'].iloc[start:stop], color='magenta')\n", "\n", "ax.plot(cropped['Sensor3'], cropped.index, color='black', label='cropped')\n", "ax.axhline(start, color='green', label='start')\n", "ax.axhline(stop, color='red', label='Stop')\n", "\n", "ax.set_title(\"Raw Active NIR\")\n", "ax.set_ylabel('Array Index(e.g.Time)')\n", "ax.set_xlim(1500, 3500)\n", "ax.set_ylim(5000, 23000)\n", "ax.legend()\n", "a = ax.invert_yaxis()\n" ] }, { "cell_type": "markdown", "id": "102dd5a9", "metadata": {}, "source": [ "# " ] }, { "cell_type": "code", "execution_count": null, "id": "2cedaf6f", "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.8.10" } }, "nbformat": 4, "nbformat_minor": 5 }