From ff6c16b5097cc01d0fb186b4b8ffd2706aa92443 Mon Sep 17 00:00:00 2001 From: 8645139112c21501d9914ade8e4a034f <8645139112c21501d9914ade8e4a034f@app-learninglab.inria.fr> Date: Tue, 14 Apr 2020 07:56:13 +0000 Subject: [PATCH] New save --- module3/exo3/exercice.ipynb | 161 ++++++++++++++++++++++++++++-------- 1 file changed, 127 insertions(+), 34 deletions(-) diff --git a/module3/exo3/exercice.ipynb b/module3/exo3/exercice.ipynb index 8b1921b..c5b8608 100644 --- a/module3/exo3/exercice.ipynb +++ b/module3/exo3/exercice.ipynb @@ -638,7 +638,7 @@ "cell_type": "markdown", "metadata": {}, "source": [ - "On peut désormais faire un graphique des données de concentration en fonction du temps pour avoir une vision des donées. Nous remarquons l'oscillation mentionnée dans la consigne au début du document, nous allons désormais superposer la coubre de l'évolution plus lente. " + "On peut désormais faire un graphique des données de concentration en fonction du temps pour avoir une vision des données. Nous remarquons l'oscillation mentionnée dans la consigne au début du document, nous allons désormais superposer la coubre de l'évolution plus lente. " ] }, { @@ -649,7 +649,7 @@ { "data": { "text/plain": [ - "" + "" ] }, "execution_count": 7, @@ -686,20 +686,20 @@ "raw_data[-75:].plot(\"convert\",\"Concentration\")" ] }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Nous utilisons un filtre pour faire un fit des données. Nous importons la fonction *signal* de *scipy* et utilisons le filtre *savgol*. Nous superposons les données et le fit sur le même graphique. Cela permet de voir l'évolution lente." + ] + }, { "cell_type": "code", - "execution_count": 13, + "execution_count": 26, "metadata": { "scrolled": true }, "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "(array([ 7, 8, 9, ..., 3147, 3148, 3149]),)\n" - ] - }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD8CAYAAAB5Pm/hAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3Xd4VFX6wPHvmZKZ9F5JQui9I1IUFVEQFbuia13rqthdYXV1dZe1r/pzdZVVd1k7FhCxAoL0DqGXQAIJhPQ2KVPP74+ZTDKAAkJCEt7P8/CQc+bOnXN59J2T9577HqW1RgghRNtlONkDEEII0bQk0AshRBsngV4IIdo4CfRCCNHGSaAXQog2TgK9EEK0cUcd6JVSRqXUOqXUbF/7RaXUNqXUBqXUDKVUVKNjJyulspRS25VSY5pi4EIIIY7Osczo7we2NmrPAXprrfsCO4DJAEqpnsAEoBcwFnhTKWU8McMVQghxrI4q0CulUoELgXfq+7TWP2qtXb7mciDV9/MlwCdaa7vWOhvIAoacuCELIYQ4FqajPO5V4I9A+C+8/nvgU9/P7fAG/np5vr5fFBcXpzMyMo5yKEIIIQDWrFlTrLWOP9JxRwz0SqmLgEKt9Rql1NmHef1xwAV8WN91mNMcUmdBKXUHcAdAeno6q1evPtJQhBBCNKKU2nM0xx1N6mYEMF4plQN8AoxSSn3g+5CbgIuA3+mGojl5QFqj96cC+w8+qdZ6qtZ6sNZ6cHz8Eb+QhBBC/EZHDPRa68la61StdQbem6w/aa2vV0qNBR4Dxmutaxq9ZRYwQSllUUp1ALoAK5tg7EIIIY7C0eboD+efgAWYo5QCWK61vktrvVkpNR3Ygjelc4/W2n38QxVCCPFbHFOg11ovABb4fu78K8dNAaYcz8CEEEKcGPJkrBBCtHES6IUQoo2TQC+EEG2cBHohhGhGDpen2T9TAr0QQjQTp9tD1ye+a/bPlUAvhBDNJLu4GoDSagcALreHhmdNm44EeiGEaCZ7SrzPlu4t9f79rwW7WJxV3OSfK4FeCCGaSbHNzsD0KH+g/2RVLsmRwU3+uRLohRCimUz+ciNdE8Mpr/GmbtpFBdM5IazJP1cCvRBCNIN7P1oLQFpMCJW1TjwejTpcrd8mIIFeCCGaSMakb/B4vDdbZ2/IB2B0/r/53YpLqJkzhQ7BNb/29hNGAr0QQjSBHzcfAKDEt8IG4CLDMrrteIto+z7Clr3I37InwKz7oK6iSccigV4IIZrAzkIbAK/N2wEuOw8GfcXL5rcAmGG9jPykczBpB+xeAEFNm6eXQC+EECfAc99tC1gT73J7f04Is8Csidxv+BSLcuLpczUPll/J+xnPsebiOTD+dTAYm3RsEuiFEOI4rd1bxls/76LIZvf3vTJ3BwBd930OGz6lRlsoufxTDJdP5aK+KazZU0a7zn2h41lNPj4J9EIIcZyq6lwA7C0JvLn6xtmK0Tn/AGCS8zZi+44FpUiLCWHz/koig83NMj4J9EIIcQwKKutYvDPwadbvNnpX1OT4Ar3bowmmjjFbJ2PSTj72jGaWZ4T/+PSYEGx2F1Zz84RgCfRCCHEMfth8gOvfXYHb05CP/2RVLgC5pTXgcvDIm5/ytOUjTBU55Jg6sHvQEzx1cU//8X3aRQKgmmkh/fHsGSuEEKechHArAN9uzOfifikBr1UX58Fbv+OV4h2gABT/jryfaGsIPZMj/Mf1TI5gwSNnN9uYZUYvhBDHwOJLtzjdgXXlz0o3c/PuB6F4B9WEoGO7wIUvsSe4J9nF1SRGWP3HGgyKjLjQZhuzzOiFEOIYuN2aEZ1jya+oA6CqzokZF1OD3sDizKE2qgsjDjzC+okTAOicv5lZmftJirT+2mmblMzohRDiGLz043ZWZpeSV1YLWvPFt9/zqvmfWPKWUGaIZuuo9ygn3H98UqSV0moHVnPTrpX/NTKjF0KIY1DndNO7XSS66gBMvZeb89eDETCH8lbSs3RzxfDI+Q2lh1Ojm74M8ZFIoBdCiKNU43Cxt7SGl67sg/XTq6BuPVWGSKy9LsA84l7sKw1kFdoC0jQX9U1hXO/kkzhqSd0IIcSvKrbZ/aUNcktr8WjovOt/9K5biw6OZVTN3zFf8TYk9SEuLIgdBbZDHoQyGJqpHvEvkEAvhBC/wOHyMPhvc/1bACoFiZQSueJFALLPeJEiov3Hx4VZmLu1gF4pEYc938kigV4IIX5Blq8C5f6KWgAcDjvPmP+LctbwrXsIWVEjuO/cLv7jw63emXyHuKbfNepYSKAXQohf4PCtld9XVgv5G0j430jGGFeDKZitvR/hi7V5mBqlZU7L8M7ujSc5VXMwCfRCCOFT43AFlBqusbsY1T2BrNx8+HgCCc59lFrTYMIHWBM6saekhhGd4/zHx4dbeOHKvidj6L9KAr0QQvj0fPIHznlpgb+9Pq+cronhdMr5CCr3scHTAdstC6HzaJIjrWw7UBWwwkYpxdWD007CyH/dUQd6pZRRKbVOKTXb145RSs1RSu30/R3d6NjJSqkspdR2pdSYphi4EEKcSHVON9BQgRLghe+3s37LFsaUfwrAzLg7SE+MAfCXNIgPszTzSI/dsczo7we2NmpPAuZprbsA83xtlFI9gQlAL2As8KZS6uQ9EiaEEEdhfW45YRYT5/dMDOj/o/tdIlUNjg7n8nVlw43XLoneG65BppafGDmqESqlUoELgXcadV8CTPP9PA24tFH/J1pru9Y6G8gChpyY4QohxIlRVu2gosbpb0+Yuhyb3UVFrbdP15RxY9ACBtYsplaF8HnKoxTZGjb6Tgi3sv7J85p72L/J0X4VvQr8EWhcri1Ra50P4Ps7wdffDshtdFyery+AUuoOpdRqpdTqoqKiYx64EEIcj2v/vZynZm069AXtgW8eQb2QwTOGqQB8GXMbofHpPDqmW8ChUSFBzTHU43bEQK+Uuggo1FqvOcpzHm5dkT6kQ+upWuvBWuvB8fHxR3lqIYQ4MbKLq4ltlF8/u1s8943qzAU1X8GqfwNgD4qC027jY89o9pXXEhrUOrPQRzOjHwGMV0rlAJ8Ao5RSHwAFSqlkAN/fhb7j84DGt51Tgf0nbMRCCHEChAQZySmu9re7JIRxafg2bqz0Bvk7HQ9QcOdWuPBlNuXbeOH77QFfDK3JEQO91nqy1jpVa52B9ybrT1rr64FZwE2+w24CvvL9PAuYoJSyKKU6AF2AlSd85EII8RvtLamhrMbJvG3e+am7NIcdS2bSfv69GPBQMfh+ytqPJT02BIA3rhsIgKGZtv470Y6neuVzwHSl1K3AXuAqAK31ZqXUdGAL4ALu0Vq7j3ukQghxguwp9c7kh3SIgfnPYvz5OaYFAQ7YEX0WVb3vg307/cd3S/KusBnTK/Fwp2vxjinQa60XAAt8P5cA5/7CcVOAKcc5NiGEOCG+XJtHr5RIuiV5NwS54V1vkmG4Yyn8/BwA1aHphHY8nRXxD8GBalZml/rfbzF5c/MmY8tfSnk4rXPUQghxDB6ansmC7YUBfVenFHJ3sXc++kPSney/cSlc8Q7tEmJZuKOIG4e19x+bHGnlg1tPb9Yxn0gS6IUQpwSXx7v4r8bhIhIbz3peI0i5yekwgQ/NV5Ac5d0JqmdyJPO3FZIR27B5t8lo4IwucYc9b2sggV4I0aZpremVEsHy3SVQnIV6ZzSZ1jswlmdTEdGNJZ0eos7pIczizWRHhZhxeTTJJ3Ez7xNNAr0Qok3Jr6hlY16Fv11Z62Lz/kosdcXw/qUEF67DgwESe7PzrNfZXuJkZU5DPj7Il4cf3rn1zuAPJoFeCNGmvLMomwlTl/nbWUXezUNurXkPKnLZYezMvWlfwB+WEJrSgw2NvhTAu+2fxWQ4ZDvA1kw2BxdCtBkut4cftxyg2tGwovuD5Xu4JTmXYWVz8RgtfNbhGe4627suPiHcQlahjQdGdwk4z/a/XdCs425qMqMXQrRa1XYXm/c3StPUucgt9W77Z7O7qHO6mb1uD78reQ2AvT3/wNqqaLokeJdZRocEYbO76NMusvkH34xkRi+EaLV6PfUDADnPXQiAy9NQdzEnv5iUzW/zedBMOhv2UxWawdLk61mzagfBvpo1Bt+WfxLohRCiheqXFsW2/Ep/O6fYu2nIvWe3J3XuH4jK+4kYA2hTMPtGvsCeEhfDO8UGnKP+S6Itk9SNEKLVig8LokNcw3r3q99eBmgu2vMiUXk/UWWI4D7Hvaj71hLUcQTfbsr3L6M8lUigF0K0SnaXm7lbCwm3mvB4NDa7C4AHTV/QPX8mLoOVT7u8xM13PQIRKcSGWcgtrSUqpO2spjlap95XmxCiTfhxcwEA8eEWyvK2YVn0HF8HraePIQeNgZ96P8e6ui6M9z3xWr9ccthBqZtTgczohRCtwrq9ZWRM+sbfnvjxOgBOM+4i8sMxhO2cSR9DDigje858ntXWoZTY7AE15GdPPIOL+qY099BPOpnRCyFahcveXHpIXzuKuC7rz5hclRQknMH82GuYMHYUYcZ4cmduwqPBaGioId+7ja+u+SUyoxdCtCpae4uTKTy8ZH4bi6uS3NgRDN97J0s8fSAyldjQIHLLarCaW+fWfyeaBHohRItT53SzaV/FYV+rtFXDd5NYGvIIw4xbcFlj+Sz1cUb1TOGvl/QCQCnFpn2VxJyCN14PRwK9EKLF2bSvgoteX+xv1zhc/p+3TLsfVvyLZM8BtDmU6gvfIKvGSkWtM6A+zdPjezFhSHqzjrulkkAvhGhxnv56CwBLs4oBKK5ycH7PRJ7psothxZ/jMZi51fEw6pEdRPQeS3GVg9JqB6rRnq43Dc9gaMdTb4XN4UigF0K0OBt9aZv6ypMLdhTS2VzMtQdeACBn4GP0HXUtWMJQSpFTUt2m6sefaBLohRAn1cx1+zjnpQWHfS062AzbvyN/9rPctOcxzM4qtkefxaeGC6lulM4prLIT0YbKCp9osrxSCHFSvTp3BzklNbjcHkxGA3VOb4nhXikRpGx7F7a9zGNmoA4cUR35KGkSZRV2Lh3QsB5+2u+HkBEbcpKuoOWTGb0Q4qSqL12QXVwNwIjnfgJgnGchg7a9DMCS4LNg1BPUXDebaevKqLa7GNQ+xn+Os7rG077RHq8ikAR6IUSz+jpzf0D7jM5xjOgcS2WdE4CSagdnG9Zzd+U/AKg6+xk+TvsLjHyU8DjvLL78oBU24tdJ6kYI0WzyK2qZ+PE6Lu7XkHYpq3GSHBnM/K0FDMqfzidB0xiotqM8br4Ov4bE9jeQVlcIeJ9y7RgXSmm142RdQqskgV4I0WzqZ/PVdhehFhN2l5vdxTb+MLI9Yd/8AcrmMbQ+zzDkTj7Yeynn76sI2BikR3IEWw9UHubs4pdI6kYI0Wy0hpjQIHYUVAHw5dp95JbW0nPTy/Qum4e2hDM15hG4dzWMewGtFBvzygMCvcmocLn1ybqEVklm9EKIZmOzuxjaMYb8ijoGOGuJdpUwOWk1EevfxoWR7HPfpbAwDeK8m3WvzC4lyGjglWv6+88xOCOG0FNw85DjIf9aQogmU5+iqffNhnzO65kIhVvh+zsYW13kf+29iLtxVnciOaqhENnVg1OZvjov4InXG4a2b57BtyGSuhFCNIl5Wwvo9dQPFFbWAeD2aHYXVzMoppaRK++C6iJc2oDbaIWzJjE39EIOVNQxonND2YInL+7F7IlnnKxLaDMk0AshmsTt/1sNwHZfPn7a0hwsOBi+4h7C7AV4Uk+nl/09tt6yHc6ZzJb9lewsrCItuuHBpzCL6ZStIX8iHTHQK6WsSqmVSqlMpdRmpdTTvv7+SqnlSqn1SqnVSqkhjd4zWSmVpZTarpQa05QXIIRoGSpqnAHttBhvwN5RYINNXxDz8+MsSptKWNkWcjyJvJf+d5Jjo+idGgVA+9gQNuZVSP69CRzNv6gdGKW1timlzMBipdR3wDPA01rr75RS44AXgLOVUj2BCUAvIAWYq5TqqrV2N9E1CCFOMpvdRb9nfmTbX8f6N/sY2zuJ9JgQQjZ/Anv/zqUARYAy8GbsJPpHJvDA6GT/OdKiQ9i8X5ZNNoUjzui1l83XNPv+aN+fCF9/JFD/uNslwCdaa7vWOhvIAoYghGizVuwuAfBvFuJye1i+q4RRHUMZlfcvAFaZB8PAG+GaDyiM6E1ZjQOruSEEvXJNf5ZNHtX8gz8FHFWOXillVEqtBwqBOVrrFcADwItKqVzgJWCy7/B2QG6jt+f5+g4+5x2+lM/qoqKig18WQrQiNQ7vL+yvzdsJQEWtk5SoYJK3TSPSU4Y9cQD/af88jH8dul9IfJiFVTmlJEUG+88RHGQkuVFbnDhHFei11m6tdX8gFRiilOoN/AF4UGudBjwIvOs7XB3uFIc551St9WCt9eD4+PjfNnohRIsw8eN1xIQG+Tf6qKxzkWCuhSWvAbC770NkxIX5j48Lt7BgexH9UuVGa3M4plU3WutyYAEwFrgJ+NL30mc0pGfygLRGb0ulIa0jhGgDLvnnYopt9oC+Hsnh1FRXwewHiZp+ORP3PgB1FWwPHsBS3YfEiIaNQeZv89auabw+XjSdo1l1E6+UivL9HAyMBrbhDd5n+Q4bBez0/TwLmKCUsiilOgBdgJUneuBCiJOjzukmM6+CL9bk+fsMCu45pzPDdr0Gq98junA5cdU7QRlZ2H4im/ZVBAT6ab8fwse3Dz0Zwz8lHc2qm2RgmlLKiPeLYbrWerZSqhx4TSllAuqAOwC01puVUtOBLYALuEdW3AjRdpT4Kkc++9027jyrE3VONwnhVk4r/w5z2QwwmPlf6M1cNyAWU5dzce+OZeOaPG4c1vBEa2KENSDwi6Z1xECvtd4ADDhM/2Jg0C+8Zwow5bhHJ4RocYqr6lM2Gha9jFr4KtPtFsxfexdVuM+fwuzMvtw4ahgACUV5ZBXaSI2WHaBOFnkyVghxTJbvLuGGoe250zgb5j2DxVlJusEb5N9VV7Cj/bU4XB7/8Qnh3pl7fLjlpIxXSFEzIcQRlFY7CLeaMBu988Jnv9vGZ5eE0M80HYA/OW/l5qGpdE1LYsv2niyfs4P1ueX+94dZJcycbDKjF0L8qoF/ncPVby/zt4Opo/+y+wlSblyDb2dN3KV0vfgh6H8dkSFBxIQE8WqjssK9UiKYec+IkzF04SOBXgjxi7T2PgKzbq93hr499wCPmT7BXLmH/ZZOzG13t79oGUBUiJkdhVUkRzbcaDUbDfRPi2regYsA8juVEOIXFdnshFtNVNW5YPV/6DD7UbqZvMXLFnZ5jL0FgXu3llY7WLe3XFbUtDAyoxdC+F36xhJmrtvnb2cV2rhleAYXR2XDNw8RhBNnZAcY/09cqUPZW1rDTY2WTe4q8pbFSoiQG68tiQR6IYTf+txyFmcV+9uvz8siVlXyjOMl0B5mhl6N4f51MPAGkiKsrM4pY3BGjP/4P43rQdfEMEKCJFnQkkigF0IEWLzTG+i35leyfHcR1+RNIdpThjttOB+H3YjR4C1bEBsWxIHKOoZ1atgRqkdyBD8+eNZhzytOHgn0QggANuSV0z42hNiwIHA72bPsC142v4V1z3yqjRHsPec1LJaGlEznhDCuHZJOXJikaVo6+f1KiFNUxqRvmHBaGs9d0ReA8f9cQrjVREaYG96/jLE5i8AIoJjd4QmslaEkRdT63x9uNfPs5X1OzuDFMZEZvRCngIxJ31DraCg5Vb9s8pNVuQHHfX5jd54qnQw5i6g2RVHS8ya45Vvyk87h/k/WsyVfdoBqjSTQC9HG1Zcj6PHk9/6+Zbu8O0INbh8NgKckm/4qi65zb6KjcweuiHT+HP8aIZe+Au2HU+f0nkPWw7dOkroRoo0rqKw7pC+3rIaHz+vK8l2FuL99DOPKt5hpAfZDRVASeeM+I29BGcFB3v1fo0PMANjqXM05dHGCyIxeiDZGa01Jo01BymocnNXVu4vbvvJanG4Pj32xkX7JVh6reg7jyrdwY2S3JxnSh/P9wH+x3xMdsFfcSN/7/zi2e7NeizgxJNAL0cYs2F7EoL/N9bfX55bTz5dyqamto3rPWrqrvQxedjd9qxbiCorgOvuf2H/DYvj9d4Qkd2fz/gr/LB68yyafuLCHPPHaSknqRog25pb/rgpoP/nVZs7vmciTg92kfTwKa+VuvrcAuVAXFMOCoe/i3m7ljC5xAHSKD+Mvszb7vxzq3XZmx+a6BHGCSaAXoo3SWqNKd9NfZTG0fCU35vwLk8dOrSUeg7sOS1QyW4f/k3UFEVjMFf73dYwPpaTawfh+KSdx9OJEkkAvRBvi9niXTQ5Ij6Ly+78RueIl703WUu/ru9Mu54HK6/jzpQM5LSOGuLJaVi5fR9eEcP85rGbvDViDQTbubisk0AvRymmtUcoblL9c692we6BtAZErXgKgwppKpNnN5xE38Z7tDCJCTQxIjwalSIq0sm5vOaN7JAacc0yvRM47qE+0XnIzVohW7JU5O7j8X0v97f8uzSGFYiY7/gnA/A4Ps/mqn+HhbfS48G6GdIjB4wGTb7eo+l2j0mIC93N9+4bB/qWVovWTQC9EK/bzjiLszob9Wc/oEsf/Ur7E5KphW9RZfG25mDTfptzxYRb2lddiNh2akkmNDm62MYvmJ4FeiFaiosbJjHV5AX0Wk6Fh5q01pYvfo3PpAjzmMKYn3Mf+yjr/bk/x4RY25lWQEhkY1LOfHcfA9OhmuQZxckigF6KV+H5zPg9+monHd8MVQANmo6KuvADXfy/mRfNUAOwjHyfXFRWQplFKcaCyjnCrOeC89fl90XZJoBeilXjsi40AHPCVNNBa43B56BwThP74Wkx7FlFjjIDzp2AZfhdzthRgPUyePbvY1qzjFiefBHohWiCn28P01Q2VJV1ubx7+tIxo9pV7SwVf8sYS1ueWMaH0XwQXrKE2OInPhnwGw+/F4JvF920Xeci537phUDNcgWhJZHmlEC3QgGfmYLO7uHpwGgDTlu0BYHz/dhwoq4LKeQwqW89zvWrpuesz3MrMtNSnSU1ODzhPp4TQgPaWZ8ZgMclqmlONBHohWiCb3VslsqLWSWSwmao6JzcOa08HQyHdfpoIVVt5CmCX9/gvM55kT3BPRidF+M/x7k2DGdQ+8Car7OV6apLUjRAtVIe4UPaUVAOwLb+KkVFFnD7/WuKrtuIITmCJuxcEhVN77hS+18Moq3YSExrkf/+5PRKJCgn6pdOLU4h8vQtxkq3KKeWbDfn8ZXwvwDuLB/j9iAz2FJTRd9urTMr9hvTc/RgcNpbr3mwe/CYbijyMmDAAq9bYtixnRXYpb/xu4Mm8FNFCHXFGr5SyKqVWKqUylVKblVJPN3ptolJqu6//hUb9k5VSWb7XxjTV4IVoC1bsLuG/S3P87QMV3lU1XSI1p82/Dha/QoZjBwaHjbrO47jJ/ggeczgX9kkGvMsjV2SXkhxpxSj1acRhHM2M3g6M0lrblFJmYLFS6jsgGLgE6Ku1tiulEgCUUj2BCUAvIAWYq5TqqrV2/8L5hTil+RbUYHe5sZiMrN1bhhE3fZY/SKhtC57IdP5lvol7rhqHNaEH/aYuJ6+shmGdYgPOI9v8iV9yxBm99qpfeGv2/dHAH4DntNZ233GFvmMuAT7RWtu11tlAFjDkhI9ciFbq+ndWBLRfmbuDfmlRrNtbDsCu3bv5v/D/Ebp3PlXGSHZd8CEFaWMhsSf4Hm7aXlBFRlzgihpXowephGjsqG7GKqWMSqn1QCEwR2u9AugKnKmUWqGU+lkpdZrv8HZA463l83x9QrRpszL3H9J3wWuLyC2t8bcLK+tYnFVMnbPhF9whHWL43ZB0DhQVw/uX8cS2S7jQOQdtDOIfMU+x3hZNz+SG1TTpMSHsLqomzBL4C/nBbSHqHVWg11q7tdb9gVRgiFKqN960TzQwFHgUmK68z1IfLkl4yFRDKXWHUmq1Ump1UVHRb74AIVqK+z5eFxDA52wpYGt+JbllDYF+xrp9AGQVen9JLq9xsDK7lOQQN0OW3AG7fsKJCR3XFXXlf9hs6snOQhtdEhvqxXs8mqqDNul+8cq+XD+0fVNenmjFjml5pda6HFgAjMU7U//Sl9pZCXiAOF9/WqO3pQKHTHW01lO11oO11oPj4+N/4/CFaBlqHd4AX+lbMQNQ4CtVUG1vCP5v/exd+L6ryBvo95bWYMDDwOX3kVK5HldYCqPsL6LuXQU9LmJldim7Cm2kNyojHGIxclqHmIDPv2pw2iFr5oWod8Tf9ZRS8YBTa12ulAoGRgPPAzZgFLBAKdUVCAKKgVnAR0qpf+C9GdsFWNlE4xeiRagP6kU2Owm+DbTtLu9d1rJqh/+4shon4/uleB+IclTz4BvTecr0I6G5P1NpiGLrOf/jmorAm6rzthUSF9awHv4vF/dC0vHiWBxNUi8ZmKaUMuL9DWC61nq2UioIeE8ptQlwADdprTWwWSk1HdgCuIB7ZMWNaOsmf+ktOLZlfyW9Urz1ZXYWVPH4uB6U1jgCju2RHMGMr2fxuwX/YJ7Fu8efNph4JWoyGfZ4MuIOfcipcYXJ+mqUQhytIwZ6rfUGYMBh+h3A9b/wninAlOMenRAtlM3uCrj5uWx3CdCQe9dasyGvgov6prAoy3sPyun2MLxTLNekFHK96VmoraZSRRIRl4wa+Sj/+TCYGwptDGqfdugHCnEc5Da9EL9B76d+YMbdw717r/rcPDyDPSXeG6+5pbVsya8kNcpKUv4CePtOKgr2cYWzJ9Gfr0VRje5+EXeV38VHd50JwNndVrI1v5L02IZ8/Mx7RtAjORwhjof8DijEMaovOFY/e683eVx37NVl8POLRHz7B642zqfd97dwy95JkJ9JnKeYK4wLUQ4bi6znsP+8f5EQFeZ/f2GlnXW55UQ02hikf1qUVJsUx01m9EIcwcWvL+aDW08nMsQbgF/6YTsA2w5UAeD23Rm1FG/l2aJ74UA+UcALZmAX1KgQQs5/gpdXVPNwwhrIOIPXNwziwXJnwGqaM7vEsSW/slmvTZwaZEYvxK+oc7rZuK+Cl+ds9/fV16V5d3E2FO/E8/UDvBP5dy/mAAAgAElEQVT8OrxzLknufEjqw49R1+CM6wF9JzDG9TLu0+9mhfVMuO5TGD4Ro8FATkk1aY0CfWKEldDD7AglxPGSQC/EryiqsgOQmVcBeG+yAjw6phuXxubCu+dhXvdfRutl4Krj55DzcN38A/+23ozpnmVw+dvkOiNZu7eMiOCGlExMWBDr9pYFzOhvHNae1U+c14xXJ04VEuiFaOSJmRuZsS7P3y6orOOWERlk5nrr0NhdHsIsJu5J3MzztiegtoxVnq6sSPs9XPMBX6U/Tp4NjAblXxJ536jOrNhdQueEhnx8p/gwpq/OC7jxajIaCJYZvWgCkqMXopEPlu9l0c5iLhuQCkBhlZ30mBCGZkTCon+gN3/L5+EKpq/EosDW50YmrDqPjTeMgyATUVlbmL+9MGADkJjQIP7y9RZuPaODv89q9s6xkiODm/cCxSlJZvRCNKIUJEda/e27P1xLlNnNvcV/g3lPE3xgFd1t3ge9FyXeyLaBf2FQh3j/Fn0mo+Lpr7cElD0o8T0Ze9uZDYF+h+9GrhDNQWb0Qvi8tzibxHArHg/g8bDkh0+52LCJi9e/gsm1CndQBN+l3MuAeEW79I5sLO6HY1cJvVIaKkvafUXNJl3Q3d93Zpd4Xv8pK2D2/tgF3bl/dNdmuzZxapNAL4TPM7O3ANA92oP+9DpGbP+OEUHAfqgMSmTXmP8yPdPEOaMHgsVE3Kpc5mwt4PRGBcauH9qeacv2kBLVENRNxkMLukrKRjQnSd2IU1bGpG948qtNAX1dQ6p5vvR+1PbvcBhCKE45G/pdy4/DPiDXlEFFjYNQX+mDiGAzu4psxIdb/O+PDfP+HGFtmEMNTI9mx98uaPoLEuIXSKAXp6zQICM7Cxpq0xhxMzvlPyQ691EX04Nza6eQeeZUuOwtQuPSWLarxL/MEiA4yMjuomraNZq9x4QGMfOeEQFFyACCTPK/mjh5JHUjTkkOl4dqhxuH2wNaU5U5iy+i3yQobznVQbGsH/kuScsrGNU9AYC4cAvr9pYxuVHuvd7Be7XK3q2ipZFphjglfLY6l/nbCv3tIpv3QSiTQcHcp4iYeSP9a5eD0cKHqX/hrpl5KNWwFj4uzMK2A1WkRjcqWdA5jlWPj5aywaLFk/9CxSnh0c83cOcHa/ztwso6rh2STsSeH2HJa3iUkTUd74G7lzFq7GX0T4siKaJhmWX9xh+RjZ5uNRhUQH5eiJZKAr1oc+ZtLeA/S7L9bYdvp6ceSQ3lfi97cyl5e3bygnkqACs730/RwIkQ24nYUAuLdhazr7zWf3x97fkRnWOb4xKEOKEk0Is258mvNvP011v87R0F3oeTLGajv1ZNEE7+6niJaGVjV9Rwrt00iN7tvDtD1c/aL+id5D+HUorsZ8cdcpNViNZAbsaKNqfxTBzAbDRww9D2FFbVUZm9huDZ97DVsh1jrcYeksSPXZ+mnd3mz78bDN5gfu2Q9IDzSJAXrZXM6EWrtnRXMf9euNvfrt8UpGtiQwExp9uD2WhgpF5F+EcXE1S6DaPS6OgMSi76LzttQaRGBz7AtP7J8/zr5YVo7STQi1bt68z9TPl2q7+dW1rDjcPakxwZTI3DBbZCZq/aQerer7hu92QMrhry0i7mw5HzUfetJ7brEFbsLqVdVEjAeaNCDt2gW4jWSqYsolXzeOr/1hgMiryyWtKiQzDhpnrWHwnZ9C6TGh2f2ekuLtl8Jm+dngJKYTEZ2VdeK/uyijZNAr1oVbTWAbny3LIaxvRKpLDKTlKklY155fSNcTN+zUTiKzIBcGkDRoNib7+HmB92NcPrSjm3R0LAedvHhjbrdQjRnCTQi1Zja34lF7y2iJznLvT3Od0eOsWHkVdQTNL6j4lcuIZhkZsJrd5LtSWB95KeYEFFEl/cNZSaMsWqb7ZQ53RjbvSQ08iu8aTFSJEx0XZJoBetxvebDgBQVeck3Gr2/50caSXp50dh37fcagKqwR6eznud3+DlZVX0SA4BayTdEjVLskoY3ilwLfz/fj/kJFyNEM1HbsaKVmNfeS3j+6Ww0VdY7M0Fu4gNDWJw0QxS932Lx2hlfsj5MPj3lF09gx014ZzbPYFP7xwKNCybDDbLdn3i1CKBXrRIk7/cQMakbwL6Pl+Tx+CMaIoK98P3kwle9gpjI3Lotv7vAKzq9wy7hj0PF71CQrtOfJ25n8o6JxFWc8B5Nu6rQIhTiaRuRItUVOXdfs/ucmMxNczAU8KM9J57H1Ss4z4FLPsEgAUR47l5aSpvXR/40NPBDznNe/gszAaZ34hTi/wXL046h8vD/O2FAX1hFiOnZUSTVeitF+9ye+iVHE6f9U+TVLEOj8FMrfLeQNWdRvHvkNsY2jHGX1YY4PKB7egUH7iaplN8GOmxgWvmhWjrJNCLk+6yN5dwy39WUVHjBLyBf+b6/YzqnkhBZR0AvZ76gdMLp5O46zNqdRALR37EjPMXw4NbUDfMwKWC8HgCN/iIDDbTOUHWxwshqRtx0rnc3kJj2SXV9A+JYu3eMgakR5ESZSVvbw6se4BN5nmYtbeG/CPOuxhq7UZUsBkiUwDwaM2uouqA83aIC6VHcgRCnOqOOKNXSlmVUiuVUplKqc1KqacPev0RpZRWSsU16puslMpSSm1XSo1pioGL1qna7iJj0jeU+Db+AIgNC+Kxsd0p9M3eJ0xdzrq95Qz3rGX8yutg+7eYtR1tMMOYv/ONZyg5xdUB9WlW5ZRRWu0I+Kwbh2VwWkYMQpzqjmZGbwdGaa1tSikzsFgp9Z3WerlSKg04D9hbf7BSqicwAegFpABzlVJdtdbuJhi/aGWue2cFAHtKa4gNs+BweVi6q4QJQ9I5UGaDdR/wd+tskkI08bO+B0CnDeWMndcy76mrsQaHMGLzcnYUVNE5oUvAuZMjrYd8nhDiKAK99hbwtvmaZt8f7Wu/AvwR+KrRWy4BPtFa24FspVQWMARYdqIGLVqvzNxyABbuKGJgejSfr8kDIMNqI+6HG6BmG9cB1AAoPgy/mYFjn+SsmH1Yg703Uc1GA0VVdsIbLZtMjwnh7RsGNe/FCNFKHFWOXillBNYAnYE3tNYrlFLjgX1a68yDlrC1A5Y3auf5+sQp6ImZG3l0THf/Zh49kyN4dEw3vlq/D/atxZL5LcMNZrrM/hPBNTnoyFSmu87imoHJ0GkUX/1gwLDPFpCmWbC96JDPWfjHc5rtmoRobY4q0PvSLv2VUlHADKVUX+Bx4PzDHH643Rn0IQcpdQdwB0B6evohbxCtX2FVHR8s38tFfVMY2tFbdiDMaqJ3u0gK5r0O//4/rgCuCAKqYJMng4QJs1k4v4hrzh3oO8sy9pTU0D8t0n/eD287HbtLMoFCHK1jWl6ptS4HFuBNz3QAMpVSOUAqsFYplYR3Bp/W6G2pwP7DnGuq1nqw1npwfHz8bxu9aNHu/mAtAF+uzfP32Z1uYjPfZkLR/wGwxdQTT0g89L6Cl5JfZluVlW6JDUsiV+WUkplbTvekhtUzIzrHMap7YjNdhRCt39Gsuon3zeRRSgUDo4F1WusErXWG1joDb3AfqLU+AMwCJiilLEqpDkAXYGWTXYFoMV6du4M1e8r87dW+n90ewGXn7fc/YNSBdzDM/TMA+sJXeDL2ZQx/zIIr36PGEMruIhud4ht2h9Ialu0uIUlutArxmx1N6iYZmObL0xuA6Vrr2b90sNZ6s1JqOrAFcAH3yIqbti8zt5xX5+5kf3ktg9pH+/sfHdONPdk74K3bubN4u/+/uFcsdzE+42qSs3b6jy2rdrDtQBWnd2yoLnlu9wTmbSvEKoXIhPjNjmbVzQZgwBGOyTioPQWYclwjE63Ks995t/ObvjqPF67s5++/tLMJ4+KHwJ1PkY4gJLUPoaffQtamTnyxJo/eKQ0pmZ2FNrKKbPxlfC9/38F7uQohjp2UQBC/ySOfZVLW6AGl/Io6JpyWxrCOseDxsC5zLfGU0+7r60hy56OT+/FwwjuE3v4t9L2KPu0iWZxVTK+Uhpusf7m4J+2iggNm72kxIZzVVe7hCHE8JNCL3+TzNXn8sPmAvx1kNHDrGR2Id+bB22cyYMY5rLLeDYWbyTenkz32fQwhDSmdmJAgNuRV0KvRjD4pMjggPw9w25kdmSYbgwhxXKTWjThme0qqSYsJZkeB9zm61Tml7Cy0kVq5linFD4CuwmmwYPbYIb4H05OfJaXIQFp0Q9VIu9u7q3d0aJC/r0tiGCM6B+7+JIQ4fhLoxRH9acZGxvVO5owu3nJGOwtsXDUojRXZJVB1AOfit/i7KZPgjxeBduLuMpY/G+7n9pGd6JSSgJ63k635VQH59jDLoTdXO8WHHTKjF0IcPwn04ldV1jn5aMVejEr5A31eWQ2d4sPI2poJU29mWFU+w0yAB5bGXE7Cua+y48tNdExJAKUIs5hYsbuUwRnt/ecd368do7rJWnghmoPk6EWAGocr4Cbr4p3FALy/fI+/L7eslg6WCp4sfQyq8sn0dGRpu1vgsqls7Ps4y3PKSYq0BuzutGx3SUDqxmhQRIYEbvEnhGgaMqMXAXo++QPhVhMb/+KtLm13ubl6cCqFVQ1lhQtLSuia+zgmTzGudkOYYp/EJ7eOAoMiY/MBPlm5l76pUf7jLb5VNGkxslRSiJNBZvSnsGe/3cqqnFJ/u9bhfa5tcKMHnh78NJMZ6/ZRY3dD9kLs/x7D87uvwFSQSamlHbvOfQdMIf49WqNDgli9p4zOCQ259mG+B6CiQhpuvAohmo8E+lNEUZWdpVnFAX1vL9zNjHX7/O195TX0S43EZncFHPfSVf0YVrcQ/f7lWPYtJ0TZISKVZae/yZwcBxHBDb8YhllMVNW5AgJ954QwsqZc0ERXJoQ4Egn0p4hnv93q3/QDvDdZAT5a4d8zhryyWs7vleTNrRduRc//O+MMyxmT/xYPlv8d5XEywzyOj4d9DfevJ75DX96Yv4tuSQ1FyIKDvGmaDnGBm3KbjPKfmhAni+ToTxF2l4dgsxGtNUop3l+255Bjbv7PKq4fms7Awhnot95DeZy8GQSsAA8GcgdN4tvSUbxwRj8wmgkyGah1ugP2Zc2IDeGj20+X2jRCtCAS6E8R32zMByCnpIYOcaEkRVi5alAqmXnlVB3IIjxnDlcZdzPJYyfMPQ2A7dFnE+fKJ9ai+SH1AYwdz6V8326ifKtlwize/3y6NiorrJRieKc4hBAthwT6U0B5jYPuSeGM7Z1E8f4cOuxazMczbdz9u6uJX/8GlrdngHbwohnYAC5lZn3/p5nuPIM/nN2Z2LhQchbsojK3HIXyL5usz8Onx4T8yqcLIU42CfRt0H+WZLM1v9JfRTIzr4ILeiczwLGavl9PAmc5n1uAz59hlG8H4NoO57O9wEb/RBOZnSey0tmRvJ3F/qdZB6ZHcc3U5Vw2IHBXyMfH9ZA0jRAtnAT6Nujpr7cANAT63HLG1X5F5zV/Q6EhoSeVRXlE6Epqw9L4scNjfFjUGVO04qObhqL2lrFz2R6cbg9m303UNN+s3btXfIPbR3ZsxisTQvwWEuhbuWW7Snhn0W7evfk0f9+QDjENBxzYxPA1k+hSswiAr2Nuxnn6wyzcvIdXzw0hR7dn9eoDRIbU8cIVfQGwmAwByy4BUqK8M/tVOWUIIVoXWfPWimzIK2fbgcqAvs/W5DJvWyFOXzXIeiuzS6nb8i166tkMrlkERguei1/nQ+sEHvpsIzUEQ7tBpMRGsr+8ljlbCvw3Wetz7i9e2feQMcy8Z0QTXZ0QoqlIoG/BDg7e4/+5hLGvLgro+3Ktd+a9Zb/3CyC3tAaHy8MIw0ZMn92E8jhZFXUB3J+JYdCNuD3e1MvoHt6CYhFWE/vKazm7W7z/Jmu41Rvw48ItAZ/12oT+xB/UJ4Ro+STQt1CZueV0efy7w75WYmuoO5MWE8zvR3Qgr6wWtObMF34iKG8Z7wa9jEk72N3hOrYNeRYikgFv6iU6xMzVp6UB3uWQ2w5UBayFB1j46DmM7BK4s9Ml/QNvxAohWgfJ0bdQT361KaBdfxP08oHtyC2rJTbMQlahjdzSWtrHBBOxezb8+BLrLVVEqWoA1sVdxJyEiZyVFBjEy2qch3xeQWVdQDs9VpZMCtFWSKBvoTLzKgBYs6eMQe2jKa9xkhhhoW+7SA4UFYNtMW++v4jnRpzFkGW30tG2BoAoX2Vge48reLP2TqrzKrl9ZGf/efumRnLbmYErZWbeM4I+7SIRQrRNEuhbiAMVdSRFWv3tDnGhjOmVxPrccga1j+bLdfuotrsZULWAbvOfBFcl/wgC1rwFQKkOI7//A2ywDuTaPlFY2g2k8I0lZOZVBGzXN+veMw757P5pUYf0CSHaDgn0J0Gtw+0v/gWQU1zN2S8tYOGj5/hTJsU2O1cOascP38+GrImck1/CZYkRxCxfDoA7qT87Sx1013vwdDmPSSUTKMqP4IqBqZDq3cmp/rcCIcSpTQJ9MyutdjDwr3PY/fdx/hru//fTTgB2F9tIjw3BZncxMD2a1LxvuWP3RMBJR4BC0EYL74TcSp/zHmbhzmK6j+2OASh/exlup5szOkudGSFEIFl108RyS2sCnibdVWQDYEdhlb9Pa+iRHMEnK3PBUcPapXO5zj4d66w7MONED7yZGxyTcI55AfWHpcwJG8/SXSUMarRBCEBxld3/BCtAsNnIbWd0aOIrFEK0dDKjb0LVdhdnvjCfL/4w3B+Uc0triA0NwlbXsLlHWY2Dv1/Wm2WL5sBrExhZXeh/7f3w2xg5/Gni67IwD+sPgNlYzK7iaiYMSQ/4PLfWGA0N+7S+fu0AMuJk9YwQpzqZ0Teh5btLAChstHRx6a4SSqodzMrc732tqo4F24voVvQjt2ZNhOpC8nQ87tQhcNlUvg67km0FNro32txjSVYJ32zIJzGi4eZtWbUDiymwuNjonol0TghHCHFqkxl9E8opqeHqwansLvaua/9sdS6fr8nj9jM7YK4rhR8ex561kf+aywn5OhOAul7XcPaaC8m67RIAVn78DbY6F49d0P2Q8zeevZ/RJY6fdxQ1w1UJIVobCfRHYduBSsIsJlKjfzkNYne5GfvqIn58cKS/4uPW/ErG90vh67XZ8M3DJC5fQX91FRNjSzD9OBk8NtKANCNgMPFu6B0MOv2P3BKSH3DuLfmVATP67GfHYXcFlkd48qKeuDyBlSWFEAIkdXNUxr66iFfn7gzoe/nH7bwxP8vf/mlrIdnF1fT48/eA90nWBduLGB5bzc27H4JV7zDSuJGZlieJ+H4iIR4bjvZnMdFxL4x7Ce5awoKoS9lRaAv4Qpk4yvuwU0KjGjNKqUNqwCul/F8wQgjR2BEjg1LKqpRaqZTKVEptVko97et/USm1TSm1QSk1QykV1eg9k5VSWUqp7UqpMU15ASeaw+XB7nL72x6PJi7Mgmp0jNuj+XLtPtbsaSjZ6/Ro4sMtRAabweOmaNM8Uj37MX1wCb0cG9EhcWwP6uU//pPwW1g78j1SzrwehtwOCd1JjQ7mu435/s0+AC4b0I7+aVH+gmNCCHGsjiZ1YwdGaa1tSikzsFgp9R0wB5istXYppZ4HJgOPKaV6AhOAXkAKMFcp1VVr7f6lDzhZymschFpMATPhrk98x6juCbznq+9e7XDRIS6E8tqG+jA/bD7AlYNS/TdbAe77eB1/GtedzRvXwFtnklC4mZkAZZAb1AnjhM948psDfDo8D8ITWbw8goq8CjJiQ/3n+HhlLgCPjOnm7+sYHyalgYUQx+WIM3rtZfM1zb4/Wmv9o9a6fo3gciDV9/MlwCdaa7vWOhvIAoac4HH/Jo1n6gD9n5nD1W8v87cr67zBvL7kL8DOQhsGpahsFOjfX7aHG4a1Z0V2KQDu/E0kUcIZtT/xfPFEKNzc8CEhcXzU8TnWV4TQNTkS+l0DHc8mMtjM+txy2jda916/2XZCeMNqGiGEOF5HldRVShmVUuuBQmCO1nrFQYf8HqivqdsOyG30Wp6v7+Bz3qGUWq2UWl1U1PSrRW58byXdnvg+oK9bYjjRIQ11YKav8g67ss5JrcMNdZVs2V/JOd0TUB43OKrZvL+C1Ohg4sK8OXP7ivcwvj2C5daJ9Fz2CFZdh+59JQ91+obsq3+CiWuYtsXDEzM3kRjRkGcvrXbw3aYDAVUip1zWG4CYRrVphBDieB1VoNdau7XW/fHO2ocopXrXv6aUehxwAR/Wdx3uFIc551St9WCt9eD4+PjDvOXEWnjQ0kOtNdsLqqhoNFP/av1+Hh3Tjau7meDDK+G5NAYtuZOx8aW8Wnwb+rl07Kvep3+693bEhN6hmOc+0XBSk5VpsQ+yYsDzfLm5gnZdB0BwFDUON6XVDhIarXu/7nTvw07JkQ35+Iv6ppD51PkByyaFEOJ4HdMyDa11ObAAGAuglLoJuAj4nW54zj8PSGv0tlRg/3GP9BiszillaVbxIf0RVpN/h6XFvtfNRoXT5YbyXCKDNDek5PNo9q0E7/kJgB5Vy8iYPpokdz7K46Lf+qfo49oKHg+X2qZjcFazNfR0Mq9cBA9v4726s7j9fW/J4CCT95/3j2O9OffGK2e6JXqXSzYO6kaD8t7MFUKIE+hoVt3E16+oUUoFA6OBbUqpscBjwHitdU2jt8wCJiilLEqpDkAXYOWJH3qDxrVkAO76YA3XvbMi4PWuiWFU1rn4aZu3vEBlrYs7R3akV5gN17RL4dXe/Dd/PBEfX0Soq5yi+KEsGvYuNksCACXWdIq7XIVRu+g75xp4JpqhB7y/xPwn6Fr69uoDwdHsKamhqs7FK9f0839+fc69cZooIcLKzikXNM0/iBBCNHI0q26SgWlKKSPeL4bpWuvZSqkswALM8S39W661vktrvVkpNR3Ygjelc09Trrj5YfMBvlybx9s3DPb3FdscAcds3l/JjgIb1w9Nx1ixFxZ/TMUuC1dkRJOw4iGCVTVaGTBp70NIOzpcT2b3R8gqrSNy/Jf0rV7OvNpBBEfEEZ1bxAj7QpTWeIwWPg29jun5CbzgW/740Hld+cecHTjdDV8+F/dL5omZG+l3UN13WfcuhGgORwz0WusNwIDD9Hc+zOH1r00Bphzf0I7Onb40ycH6pUZ6677X7OeFaXPoFB7LzbXvkzFnGnjquA4gG1CwwDOAgnNexmpwc0mfOD5cVMO0GVs4s0sc95wzEKx9iNlSwNb8Sj6MmcQZd3wFwO7CKib/YyFDO8b4P/e+c7vwjzk7OLNLQ7lgi8nItr/K7F0IcXK0mRIIdpcbi8nIZa/OATQ9k8OxLXyD4CVP8T/fTJ3t3r/q2p9DWe5WklUp2X0eYJZ9HF/+sJ8Zdw+HmGgmj3OzPrccu9NDhNWbMw+1mHh5zg5uGZHh/8xw32t/GtfjkPE0vskqhBAnU+sO9I4aRhnWMi58N+7/vAG2vcyoyGVfxmiqCizEb/wGgEodQoSqoS5pMO+G3EZirzNx9nRz7aBkLDY3+oftJIRb6JvqTa1YzUYq61z0TGlY+lhf7rddVEMAD7d6//m6HFQhcta98oCTEKLlaN2Bft7TvBf0lvfZ3X0N3e0OzAXAZbSyuMdT3L+pA5kPD8BpimHtp5m031/JpQNSwBSExWRn5vp9aB24Aia7uJoRnWP97foZev36efBu7AEEbAsI+L8whBCiJWjddwO7nA8pAznQfyJ3OB6k5s6VTHA8gSe2K7WRnfhX2ov8aWc3nr9yAIQnEWY1U2V3saOgiq6+5Y0hQSa09ub0D9YtKSKgvfqJ0VzcL8XfVkrx44Mjm/YahRDiOLXuGX3nc6HzuSQBFQXLmF8UQbv+52O4+lEcNU5efuZHYkM9nNXV+0CWUgqtNQ6Xx1/9sX42/vyVfQNOvfLxc4kNtQT0NZ7N16v/whBCiJaqdc/oD/LF2jwSfGUGIkO8N0rTYkICUisVtU6sB6VaZtw9nO4Hzd4Twq3yhKoQok1o3TP6RkKCjIQEGbl8QENZnf5pUbg8gRt07Ciw+XPr9QakB26y/f/tnW2IFVUch5+frinrC6S5tZW5CooIadqiBaVWWGRUREi+gWFfAj/YC4VLfYuiJMIPBhGkBJVGL1BCZLWEH7K3XdNaV1dTzCTRxEDMMq3ThzkXx+u9u+tmO3OH3wPDnPufcw7z3Jn535kzc+81xpgiUZgz+kUzx/JJ52HGpH4NsmlUPUPK/kf1ihFDuGVSQ3+vnjHGZIbKfz4gC5qbm0NbW9t/7ufgbyfP+Xem43+eRpx93t0YY4qEpPYQQnNP9QozdAOc95+uI5zgjTGmOEM3xhhjKuNEb4wxBceJ3hhjCo4TvTHGFBwnemOMKThO9MYYU3Cc6I0xpuA40RtjTMHJxTdjJf0K/NRDtcuAo/2wOv83RfCwQz4oggMUwyMrh7EhhNE9VcpFou8Nktp681XfvFMEDzvkgyI4QDE88u7goRtjjCk4TvTGGFNwainRv5r1ClwkiuBhh3xQBAcohkeuHWpmjN4YY0zfqKUzemOMMX0g00Qvaa2kI5I6UrGpkr6U9IOkjZJGxHiTpD8kbYvTKzE+PBXbJumopNV5dIjLpsRlO+LyIbXkIGlx2br+I+m6GnMYJOn1GN8pqSXGM3Xog8clktbF+HZJc/LgIWmMpM/je7tD0ooYHynpU0l74vzSVJsWST9K6pJ0R9YeF+ogaVSsf0LSmlQ/me9TAIQQMpuAWcB0oCMV+xaYHcvLgGdiuSldr5s+24FZOXWoA74HpsbXo4CBteRQ1u5aYF8NbodFwIZYrgf2A01ZO/TBYzmwLpYb4voOyNoDaASmx/JwYDcwGVgFrIzxlcALsTwZ2A4MBsYBe7M+LvrgMBS4CXgYWNNNv/2+T4UQsk30UfycBA4c5+y9gzFAZ6V6VfqaAPxcap9Dh3nAG7XsUNbmOeDZWnMAFgIbSWwR4csAAAMJSURBVD54R8WDeGQeHC7Q42VgSapeKzAjLx6pdfgAmAt0AY0x1gh0xXIL0JKqvwm4MU8ePTmk6j1YLdFn6ZDHMfoO4J5Ynk+yY5cYJ+k7SZsl3Vyh7ULg7RDf1Qyp5jARCJI2Sdoq6ckKbfPukOYBYH2FeN4d3gV+Bw4BB4AXQwjHytrmxQGqe2wH7pVUJ2kccD3nb6dMPSQ1AdOAr4HLQwiHAOK8IVa7iiQBljgYY2ky8+ilQ2/IzCGPiX4ZsFxSO8kl018xfgi4JoQwDXgMeCs99h1ZQOXE099Uc6gjubxbHOf3SbqtrG3eHQCQNBM4GULoqNA27w4zgL+BK0mGCh6XNL6sbV4coLrHWpKk2AasBrYAZ8raZuYhaRjwHvBICOF4d1UrxMqTYSYeF+DQGzLbFrn7c/AQwi7gdgBJE4G7YvwUcCqW2yXtJTlDbot1pwJ1IYT2LNY7TTUHkoNycwjhaFz2Ecl4bGt8XQsOJSrutDXisAj4OIRwGjgi6QugGdgX6+bGAbo9Js4Aj5bqSdoC7Em9zsxD0iCSBPlmCOH9GD4sqTGEcEhSI3Akxg9y7pXI1cAvqb4y8bhAh576ynSfyt0ZvaSGOB8APA2Unq4ZLWlgLI8nGe/al2q6kJycgVVzIBl7nCKpXlIdMBvoTDWtBYdSbD6woULTWnA4ANyqhKHADcCuVNPcOEC3x0R9XH8kzQXOhBAy358kCXgN2BlCeCm16ENgaSwvJRn3LsUXSBoch6AmAN+k2vW7Rx8ceiLbfSqLGxupmxPrSYZkTpN8qj8ErCC5ObYbeJ6zN6HuB3aQjEtuBe4u62sfMCnPDrH+kujRAayqUYc5wFdV+sq9AzAMeCduh07giTw49MGjieTm4E7gM5JfMszcg2RYMpA8YbYtTvNIbny3klx1tJK6AQ48RfK0TRdwZ9YefXTYDxwDTsRtNznrbVGa/M1YY4wpOLkbujHGGHNxcaI3xpiC40RvjDEFx4neGGMKjhO9McYUHCd6Y4wpOE70xhhTcJzojTGm4PwLmsdtNLYdWqkAAAAASUVORK5CYII=\n", @@ -711,7 +711,32 @@ "needs_background": "light" }, "output_type": "display_data" - }, + } + ], + "source": [ + "from scipy import signal\n", + "\n", + "raw_data[\"smooth\"]=signal.savgol_filter(raw_data['Concentration'], 57, 1)\n", + "\n", + "fig, ax = plt.subplots()\n", + "ax.plot(raw_data[\"convert\"],raw_data[\"Concentration\"],linewidth=.75)\n", + "ax.plot(raw_data[\"convert\"],raw_data[\"smooth\"],linewidth=2)\n", + "\n", + "plt.show()\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Afin de pouvoir caractériser la période, nous ramenons les données autour de 0 en retirant aux données la courbe du fit." + ] + }, + { + "cell_type": "code", + "execution_count": 27, + "metadata": {}, + "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAD8CAYAAACcjGjIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xd8m9XZ8PHfkSx5770dO9uZzk6AsDekQEuhE9q+dDxv96STzqelLX37dDyUpy0tfVqgFFpCgYRAmCEhcYazncSO95b3tqzz/iHJsSNvy5IsXd/Px58k0u37vuzc0qVzznXOUVprhBBCiOEM3g5ACCGE75HkIIQQwoUkByGEEC4kOQghhHAhyUEIIYQLSQ5CCCFcSHIQQgjhQpKDEEIIF5IchBBCuAjyxkUTEhJ0Tk6ONy4thBBz1sGDB5u01omeuJZXkkNOTg6FhYXeuLQQQsxZSqlyT11LupWEEEK4kOQghBDChSQHIYQQLiQ5CCGEcCHJQQghhAtJDkIIIVxIchBCCOFCkoMQQggXkhyEEEK4kOQghBDChSQHIYQQLiQ5CCGEcCHJQQghhAtJDkIIIVxIchBCCOFCkoMQQggXkhyEEEK4kOQghBDChSQHIYQQLiQ5CCGEcCHJQQghhAtJDkIIIVxIchBCCOFCkoMQQggXkhyEEEK4kOQghBDChSQHIYQQLiQ5CCGEcCHJQQghhAtJDkIIIVy4LTkopYxKqcNKqX+765xCCCG8w50th88Cp9x4PiGEh/RbbXT3W70dhvAhQe44iVIqA7gJ+CHwBXecUwgxu3oHBnn9TCM7jtfx8ql6IoODeOMrVxBklN5m4abkAPw/4CtApJvOJ4SYZff95SBvnGkkJszEioxo9pyz8M75ZrbMT/B2aMIHzPgjglLqZqBBa31wguPuU0oVKqUKGxsbZ3pZIcQMVDZ388aZRj5+WS4HvnE1f/jwOsLMRl44Vuvt0ISPcEf7cQtwq1KqDHgCuFIp9b8XH6S1fkRrvVZrvTYxMdENlxVCTNf2ohoAPrAxG5PRQIjJyBWLk9h5oo5Bm/ZydMIXzDg5aK3v11pnaK1zgLuA3VrrD8w4MiHErHmuqIY12bFkxoUNPXbjslSaOvs5UNbsxciEr5CRJyECTHFdB6frOrh1ZdqIx69YnEiIycCL0rUkcHNy0Fq/prW+2Z3nFEK41/aiaowGxY3LU0c8HmYO4opFSbx4vA6bdC0FPGk5CBFAtNZsL6phc148iZHBLs/fsDyVho4+Dla0eCE64UskOQgRQA5XtlLZ3MO2VemjPn/l4iTMQQapWhKSHKZLa42ls2/C4+raeqlq6XZ5/EBZM5c+uJua1p7ZCE+IUW0/UoM5yMB1+cmjPh8RHMTWhYnskK6lgCfJYZpeLW5gw49eobLZ9Y1/uK//8xi3/fZt2nsHhh6z2TTffe4Elc097DnXNNuhCkFDey+/efUcTx+q4qrFSUSGmMY89sblKdS29VJU1erBCIWvkeQwTceq2rHaNEcqx38BVbf00NjRx0MvnRl67LmjNRyvbgfszXwhZkt77wCf+utBNv14Nz/dWczS1Ci+eO2icb/nykXJBBkUO0/UeyhK4YskOUxTeXMXACdq2sc9rqmzD5NR8djeMo5WtdJnHeTBHfYX6Zb58RyukOQgZs+Pnj/FjuN1fOySeez+4lae/Pgm5idFjPs90WEmNuXF89LJOg9FKXyRJIdpqrDYu5NO1LSNeYx10EZzdz8f2JhNfEQw3/jncR7dU0Z1aw9fv3EJa7JiKa5rp6tPVsMU7vfm2UaeOFDJfZflcf+NS8hNHD8pDHft0mRKG7s419AxixEKXybJYZrKHWMNJ2va0Xr0gbvmrn60htyEcL5981KOVbfxkx2n2bowkUsWJLA6KxabhqNVYycYIaajq8/K154+Rm5iOJ+7esGUv/+apSkA0rUUwCQ5TENXn5XGjj6So4KxdPVT3z561VKjo5opMTKYm1ekcukC+2qXX7thMQCrMmMAOFwpNeXCvR7ccZqath4evGMFISbjlL8/JTqElZkxvHRCupYClSSHaahwtBpuWGafYXqydvRP/o0d9uSQEBGMUorfvr+AZ/9jC0tSowCIDTeTmxAu4w7CrQ5VtPDnveXcszmHtTlx0z7PdfnJFFW1UdtmL7fuHRjky08V8fJJaU0EAkkO01DuGG+4Lt/e9D5RPfqgdFNnP8DQTNTIEBMrMmJGHLMqK4bDFa1jdk0JMVVP7K8gMjiIL183flXSRK51dC29dKIem03zxaeKeOpgFX/ZV+6OMIWPk+QwDRWOSqWlaVHkxIeNWbHU1Hmh5TCW1VmxNHX2US2T4YQbWAdt7DpZz1VLkggzz2wvr/lJEeQlhvPSyTp++MIpnj9aS3pMKAfLW7AO2twUsfBVkhymoczSTUyYiehQE/lp0ZwYp1sp1GQkPHjsF+lq57iDdC0JN3jnfDMt3QNcvyx14oMn4dr8FPacs/CHt85z75YcvnL9Ijr7rJysHb+EW8x9khymocLSTbZjHfylaVFUNvfQ1jPgclxTZ9+oi5sNtzglkhCTQZKDcIsXj9cSajKydaF7NtS60ZFkblqeyrduWsqGefEA7D8vez74O0kO01De3EV2fDgA+Wn2weVTo3ySauzoIyHCPO65gowGVmTESMWSmLFBm2bniXquWJxIqHnqFUqjWZ4RzY7PXcr/u2sVBoMiJTqE7Pgw3pHk4PckOUxRv9VGdUsP2fEXWg4w+kzpybQcAFZnxXCiup0+66B7gxUB5WB5C40dfW7rUnJanBKFyXjhrWLDvDgOlDXLwnx+TpLDFFW39mDTkOXoVkqKDCExMnjUmdJNnf3jDkY7rc6MpX/Qxm9fLeHXu8/ywPYTslWjmLIXj9diDjJw5eKkWb3O+nnxtHYPcEZmT/u1mZUzBKByi71SydmtBPaupZMXtRwGBm00d00uOazJjsVkVPzylbMAKAUljZ385aMb3Bi58Gdaa3Yer+OyBQlEjFMA4Q4b5tnnTuw/38zilKhZvZbwHmk5TJFzAlxO/IWN2fPTojjX0EnvwIVuoeaukXMcxpMYGczuL17Om1+5gtPfv56712dxpLJVmu1i0oqq2qhp63V7l9JoMmJDSYsOkXEHPyfJYYrKLd2Emowj3vSXpkZjtWnO1ncOPTZ8dvRkZMaFkRkXRojJSEFWLB29Vs41dk78jUIAO47XEWRQXLNk9E183Ekpxfp5cbxT2iyTN/2YJIcpKrd0kRUXhlJq6DFnxdLwZTSGr6s0VQVZ9rkPB8ulgklMzt6SJgqyYokOG3sTH3fakBtPU2cf55u6PHI94XmSHKao3NJN1rAuJbAPToeajJyuuzBA1+RoOSROsuUw3LyEcGLDTByS5CAmoavPyvGadtbNi/XYNdcPG3cQ/kmSwxTYbJqK5u4R4w0ABoNiYUokp2svJAdnyyEhcvx5DqNRSlGQFcuhCkkOYmKHK1oZtGnWzWCRvanKTQgnISKYfaUWj11TeJYkhylo6Oijz2oja1ilktPi5EiK6zuG+mCbOvoJNxunvb5NQXYsJY1dtHb3zyhm4f/2lzVjUPaqN09RSnHN0iSeP1ZLcZ2UtPojSQ5TUOYsY40Lc3luUUokzV39Qy2Gxs4+EqYx3uBUkGV/ocuyGoFHa83gFCrVDpxvZklqFJEhnhlvcPrStYuIDDHxlaePTileMTdIcpgC59ag2fGuyWFxSiTA0Keopo6+aY03OK3MjMZoUNK1FICeKqxiw49ept868cqn/VYbhytbPNql5BQfEcx3bllKUWUrf3zrvMevL2aXJIcpKLN0YTQo0mJCXZ5bdHFy6OybdBnraMLMQSxOiZTkEIC2F9XQ1NlPyyS6FI/XtNE7YBsaIPa0W1emcfWSJH6+q5gyqVzyK5IcpqCwvIWlqSPXmXGKjwgmMTKYU45BaXu30tQHo4dbkx3LEcdgI8D/vFHKlT9/bcRkO+FfuvutQxVAk0kOBxzHeqPlAPaxhx+8azkmg4H7nzkm8x78iCSHSerut3K4ooXN8+PHPGZxSiTF9e30W220dg+QGBEyo2sWZMXS1T9IcV0H24tq+OELpyht7KKoUsYh/NW+Ugv9jo10nLPsx3OgrJl5CeHTmk/jLinRIXzx2oXsLbXI3Bw/IskB++5Zj+45P25l0IGyFgYGNZvzEsY8ZlFyJGfrO2no6AWmV8Y6nHNQ+vdvlfKlvxexyrExkCzK579eL24c+ntrt+seIcPZbJoDZS2sy/FcldJY7lyXSVRIEI/uKfN2KMJNJDlg7y767nMn+cifDtDTP3qXzdslTZiMatwX4uLUKPqstqFPTzMZkAbIjAslISKYZw5VkxEXyp/uXcfC5Aj2l8mnM3/1+pnGoRn3E3UrnW3opK1nwGtdSsOFmYO4e30WO07UUSNb3vqFGScHpVSmUupVpdQppdQJpdRn3RHYbPjX4Wq+99xJl8dLHGsYHapo5TNPHB61LO/tcxZWZ8WOO2/BWbH01tkmgBmVsoK9P3djbhxx4WYevWcdMWFm1s+L42BZs+zh64fKmroos3TzrlXpALRM0K2039GC9NZg9MU+uCkbrTV/2Vfu7VCEG7ij5WAFvqi1XgJsBP5DKbXUDed1uycOVPCXfWUMXPTGWtLQRajJyAO3LGXXyXq+s/34iIG11u5+jte0sWWcLiWwb8huUPDWOXtymGnLAeDHd6xg1+cvG1oifF1OHF39g0MD38JzBm2an79UTJNjLou7vXHW3qV0zdJkwsxGWiboVjpwvpmkyOChvUW8LSM2jGuXpvD4/ooxW+Bi+ib6sOBuM04OWutarfUhx987gFNA+kzP625aa07WtDMwqF1K7kqbOslNDOeeLfP4+NZc/ndfBY/tvfDpZ1+pBa1hyziD0QAhJiM5CeHUtjnGHNyQHCKCg4gfdp6hNW1k3MHjzjZ08Kvd53jxeN2snP/14kay48PISQgnNsw8brdSUWUrO0/Uccn8hBGLQHrbvVtyaO0e4F9Hqr0dit/Z6+GlStw65qCUygFWA++M8tx9SqlCpVRhY2PjxU/PuqqWHtp7rQCcqR+5FHZJYye5iREAfPW6xWxdmMhPdxZT325/k99zzkKY2ciKjJgJr+PsWooIDnLbPr7DpUaHkhkXyv7zsqaNp7X32O+f2lnoU++zDvJ2iYWtCxMBiA03jTkgXdPaw8ceKyQxMpiv37TE7bHMxPp5cSxJjeLRPeelrNXNnD0SnuK25KCUigCeBj6ntXbZUFlr/YjWeq3Wem1iYqK7Ljtpw7fxLK6/0CXTOzBIVUsPeYn2bhuDQfG9bfn0D9r40QunAPtg9Pp5cZiDJv51OXfGSoiYWaXSeNbnxFNY1iIvPg9r77G/WTtbhu5UWNZCz8Agly9yJIcw86ilrF19Vj7250J6+gf54z3r3NI6dSelFPdszuZMfSdHq1y3zhXT9/ZcTA5KKRP2xPBXrfUz7jinu52sabfPbo4O4eyw5HC+qQutIc/RcgD7FqCf2JrHs0dqePZINSWNXROONzg5Z0rPZt35+nmxWLr6KWmUGame1OZIDrNRjfPCsVrMRgMbc+1dlzFh5lFLq7/0VBGn69r51ftWszA50u1xuMNljtaPzO53n6qWbsocy/d4ijuqlRTwB+CU1vqhmYc0O07UtJOXGM6y9GjODEsOpY432OHJAeBTl+eRERvKl586CjDu5LfhnN1Ks/mJzlm6KGvpe1Z7ryM5tI2fHN4408jnxqh6u5jWml+9cpa/vlPBbavTh6rh4sJMLgPSzV39vHi8jk9szeOKRUnT/ClmX2p0KMlRwRyRyZpu8/Y5z3cju6PlsAX4IHClUuqI4+tGN5zXrU7UtJOfFs3C5EjKLN30We3VFM4y1nkJI5fhDjEZeeAWe/dSbJiJJZPcSD0zNozoUBMZsa7rL7nLPMda+jIZzrOcYw51bb3j7u/9wrFa/nWkhh0TDFxrrfnxi6f5+a4z3F6Qzg9vWzb0XEyYmbaegREly84xsGXp0TP5MTxiVWaMJAc32lPSNKtd1aOZ3mYDw2it3wJ8p1xiFJbOPurae8lPiyIpKoRBm6a0sYslqVGUNHaSHhM66uDx1UuTuXNtBvERwRgMk/sRDQbF05/c7JYy1rHY9/CNlZaDhzlbDgODmqbOPpKiRl8exbl15sOvl3Dj8pQxq4l+9MIp/ufN83xwYzbfvTV/xD0W69jus61nYKharc6RHJKjfGucYTSrs2LZeaIeS2ffiGo7MXVaa/acs7A5L4GDHrxuQMyQPlFjHx9fmhbFwmR795Gza6mksZO8pIgxv/fBd6/kq9cvntL15idFzPpevutz4qhu7aFaZqN6jHPMAaBmnEHpMksXUSFBHKtu4+2S0bsDOvus/OntMm4vSOd72/JdPnzEhts/JQ7vWmoYSg4zW7PLE5xLvRRVSethps7Ud9LU2ccl8yc37ukugZUcUqOYlxCO0aA449i1rbSxi9wE153dfF2BY9evY/Li85j2ngGMjjfxsQalu/ut1Lf3cc/mHBIignn49ZJRj3un1MLAoObdBRmjtixiw+zJYfigdH27Y19yLy6yN1krMqIxKDgim1XN2B5HldJkxz3dJUCSQxvpMaHEhJkJDjIyLyGcM/Wd1LX30t0/OG7LwVfNd8R8rqFzgiOFu7T3DgyVPI+VHMqa7BUli1Ki+MglObx5tonj1a4lnW+ebSLEZGDNGGt1OZPD8HLWuvZe4sLt97CvCzMHsSglisMy7jBje841kR0fRkasZ2fCB0RyOFnTPrSYGcDC5AjO1ndQ0uCsVJp7LYcwcxDpMaGcleTgMW09VrLiwgg1GalpHb1bybmVbE5CGO/fkE1EcNCorYc3zjSyMTd+zDf6GEe3ZOtF3UpzoUvJyTkoPd7gvRjf+aYuXjvTyDVLkj1+bb9PDl19Vs5bushPu1DhsSApkvLm7qGJcfMT517LAWBBcgRn6yU5eEp7zwBRoSbSYkKoHaOc1TkYnRMfTnSoifdvyOKFY7WUWy7MSals7qa0qYvLFow9GfTCmMPIlsNcGIx2Wp0VQ0evldImuUen6xe7zmA2Gvj41jyPX9vvk8Op2na0ZkTLYVFKJFrDSyfriQgOmhN9uKNZkBRBSWOnbO7uIe29A0SFmEiLCR2nW6mLpMhgwoPthYAfuWQeRoPi929e2GPZuQzCZQvHHmAMNxsxGw00XzTmkDKHWg6rHYPSh2XcYVpO1rSzvaiGe7fkeOU9yu+Tg3MwOj99ZLcSwMHyFvISw31q4bKpmJ8UQZ/VRnWLVCzNNptN09lntbccokPHrFYqs3SRM6zAITkqhNtWp/PUwUosjtVc3zjTSGp0iMvEy+GUUsSEmWjtsncrWQdt45bP+qK8xAgig4NkvsM0PbSrmMiQID5+medbDRAQyaGNuHDziE9c2fHhmIz2hDDeC9TXzU+yz8Y+2yDLd8+2jl4rWkNUSBCpMSE0dvQNTaQc7nxTN/PiR45h3XdZLr0DNh7bW4510Maec01cumDi1VSHr8za2NmH1nNjjoOTwaBYmRkjLYdpOFjewsunGvjE1rxZL4sfi98mh47eAf77tRJePF7HsvToES9Ek9EwlBRy5+BgtJOzYkkGpWefcwJcdKi9Wwmgvm3kvg4dvQM0dfaNaDmAPYlfvSSZx/aW8c75Ztp7rUPrD40nJuzCyqzOMta51K0E9kHp4voO2d9hin62s5iECDP3bM7xWgx+lxy01vz2tXNs+fFufrLjNAVZsXznFte9hxY4Fi2byy2H6FATyVHBMijtAc4JcM5uJcBlAmK5Y2G0eQmuJYef2JpLS/cA9z9zDKWY1EKOceHmoTGHura5MwFuuNVZMQzaNMdGKecVo6ts7mZvqYWPXZo7NHblDd678iz5e2ElD+4o5srFSXzu6gVj7sGw0PGpey7OcRhuQVIk56RbadY5Ww5RIaahrp2LK5aGKpVGmVS5NieONdmxHCxvYWVG9FA10niGr8za0DFXk0OsfXfEs40+s52pr3Nu6nPlYu8uruhXLYdySxfffe4km3Lj+f2H1o67Oc8dazL4zFUL5mwZq9P8pAjONnTK3g6zzLmXQ3SoiVRHy+HifR2cySE7bvSuyo9flgvApeOUsA4X6+hW0lpT396L0aCIn0RS8SVx4WY25cWzvahG7tFJ2ldiIS7czAIvf3D1m+QwaNN88e9FGA2Kn925csKF8tJiQvnCNQsnvaCer5qfFEF3/+C4a/2ImXOuyBoVat/hLy7c7NKtVNbURWp0yJg7AF69JJnvbcvnw5PsR44NM2O1aTr6rNS19ZEUOfkFIH3JtpXplFm6ZfOfSdBas6/UwsbcOK9XUfpNcvjdGyUUlrfw/W3LSI+ZveWyfc0CWUbDI4a6lULtlSOp0SEu24Wet3SREz92gYPBoPjQpsnXrA9NhOvqp6Fjbs2OHu66ZSmYjQbZV3oSKpt7qGnrZVOuZ9dRGs2cTg5aaw6WN/PNfx3jF7vOcNOKVLatSvN2WB7lHFgfvrudcL/2ngGUggjHZjz2iXAjW2tlTV2jjjdMl3PZ7pbuAera5tbs6OGiQ01csTiR54pqZcLmBPY5xhs2SnKYvlO17Vz+s9e447/38o+DVdy8Io0fvWu515tinhYXbiY+3Cwth1nW1mOfHe3s1kmLDhmxI1xb9wAt3QOjVipNV0zYhSU06tt751wZ63DbVqXT1NnH3jGWMBd2+0otxIebh8rUvWnOVis9sb+ChvY+fvaelVyXn0xkiHcmiviCPMegtJg97b1WokIvvFzSYkLp6LXS0TtAZIiJ85YLayq5i7PlUNvaS3uvdU7Njr7YlYuTiAgO4tkj1VyywLP7EswVWmv2llrYmBvvEx9y52zL4UhVGysyonn3moyATgxgH3c469ifQsyOdkfLwSk1ZmTFUpmjUuni7WZnIs4x5lBcZ18CZq6OOYB9293r8lPYcbyO3gGZEDeaiuZuatt62ZjrGyW/Ppcc+qyD/ObVcyNWsbxYv9XGqZp2VmaOXaoaSBYkRdDea6Wxo2/ig8W0tF2UHNJj7G/Uzoql801dKAWZce7rVooKMWFQcKrOPp40l7uVALatSqOjz8prxQ3eDsUnOccbNuV5f7wBfCw5aK356j+O8tOdxdzz6IER2zIOd7qunf5BGyvHmccQSJyD0jLuMHvaeweIDh3WcnDMdShv6mJviYU3zjaSFh1KiMl9G/EYDIroUBPFjuQwVweknTbnxZMQYea5olpvh+KT9pU2kxBh9plVG3wqOfz8pTP860gNtxekU9nczeeeODxqdUORY5XHlZnRLs8FImc5a5HUkc+a9p6RYw5JkcEYDYoHnjvJ3f+zj6LKVq5fluL268aGmYc+JCVHz+2WQ5DRwDVLU3ituEG6li6itWZviYUNPjLeAD6UHB7fX8GvXz3H3esz+fl7VvLArfm8WtzIQ7uKXY4tqmojPtwcUPMZxpMYGcz6eXH8/s3SoXp84V7OvRycgowGPn3lfO7ZnMMjH1zD4W9fy7dudl3Da6accx1CTUYivbjOjrtcl59MV//g0L7Iwq7M0k1de69PlLA6eT05aK352zsVfPNfx9m6MJHvb1uGUor3b8ji7vWZ/ObVEnYcH9kMLapsZWVmjM9kWG9TSvGtm5bS3N3Pb3af83Y4fmdg0EZ3/+DQBDinz129kAduzefa/JQRXU7u5KxYSo4K9ov7fXNeApHBQew8UeftUHzK80drANg6yaVVPMGryaGnf5AvPXWUr//zGJvz4vnN+wsIMtpDUkrx3VuXsTglkl/sOjtUidPZZ+VcYycrMqRLabjlGdG8uyCDP+45P+5gvpi64esqeZpzrsNcrlQazhxk4IrFSbx8qgHroM3b4fgEm03zZGElG3PjyIp3X0HDTHktOdS09nDbb/fwzOEqPnvVAv5073oiLmo2m4MMfGhTDsX1HUO7SR2rakNrpFJpFF++bhEmo4H/fOG0t0PxK+29F9ZV8rQLLQf/SA4A1+Wn0NzVT2F5i7dD8Qn7Si1UNvfw3nWZ3g5lBK8lh0feKOV8UxeP3rOOz1+zEOMYC4rdsjKVUJORJw9UAnC0yjEYLZVKLpKiQviPK+az40SdzER1I2fLIcoL82mcYw4pc3wwerjLFyViDjJI15LDk4WVRIYEccOyVG+HMoLXkkNheTNrsmO5fNH4a5ZHhpi4eUUqzxXV0NVnpaiqlcy40KEJQmKkj14yj9ToEH73Rom3Q/Ebwzf68bRYR7dSkhc2mJ8t4cFBXDo/gZdO1Af8xM227gFePF7Hu1alu7UM2h28khxsWnOqtoM12bGTOv6u9Zl09Q/y/NFaiirbpNUwjhCTkauWJLH/fDMD0qfrFsO3CPU0f+xWAnvXUnVrDydq2r0dilc9W1RNv9Xmc11K4KXk0NM/yKBNUzDJ5FCQFcv8pAh+90YJ1a09khwmsCk3ge7+Qdma0U2G9nLwQrfSguRIQk1GlqZFefzas+mqJUkYFAHftfTkgUqWpkaxLN33Cmy8khy6HJuNF2ROLjkopbhrXSYljfYqHBmMHp9zbRYZd3CPC3s5eH5AOi8xglPfv95nZs26S3xEMGtz4ng1gJfSOF7dxomadu5a73utBnBTclBKXa+UKlZKnVNKfW2i47v7rcxPiiA6bPKfxG5bnY7JqDAoWJbuX5+i3C0+IphFyZFDa7WImWnrGcBkVIT6WJ/wXLc4JZIKS7e3w/Cav75TTnCQgW0r070dyqhmnByUUkbgN8ANwFLgbqXUuFNFu/sHWZM1uVaDU3xEMNtWpbM2O44w89yfKTrbNuXFU1jWQr9Vxh2m6lxDB88V1Qz927kiqz9MQvMlqdGhtPda6eyzejsUj2vu6ueZQ9XcXpA+pQ/JnuSOlsN64JzWulRr3Q88AWwb7xsGbXrSg9HD/eSOFTx+38bpRRlgNubG0zMwSJGj9FdM3sOvl/LZJw7T1m3vTrLv5eCbL+C5LM2xsu3F260Ggr+9U06f1cZHtszzdihjckdySAcqh/27yvHYCEqp+5RShUqpQoCC7KmPGxgNasz5EGIk+wblMu4wHecaOrFp2FNiX/98pfvTAAAbCElEQVTHvly3tFbdzbk2Wk1b7wRH+pd+q43H9pZz6YKEoRWVfZE7ksNo79Yuxcta60e01mu11muNSpGb4F8DbL4mJszM4pQoSQ5TpLWmpNG+9PnrxY2Ao1tJWg5u59wwqSbAWg7PH6uhoaOPj17iu60GcE9yqAKGD7dnADVjHAtAmNk4tBevmD2bcuM5WNEiyyNPQWNHHx29VoIMitfPNKK1tq/IKsnB7ZIjgzGowEoOWmv+8NZ58hLDucyHFtkbjTuSwwFggVJqnlLKDNwFbB/vG8L8YOnhuWBTXjz9VhuHK2TcYbLOOVoNNy5Ppa69lzP1nfa9HAJ8K9rZEGQ0kBwVQk1r4HQrHShr4Xh1O/dumefzH5BnnBy01lbg/wI7gVPA37XWJ8b7njCzlAR6wvp5cRgU7JWS1klzzqW5d0sOAK8VNzi6leQDzWxIiwkNqJbDn/eWER1q4o6CDG+HMiG33PFa6xeAFyZ7vCQHz4gONZGfFs2+Egtc4+1o5oaShk4igoNYlRnDouRIdp6oo3/Q5pWlMwJBanQIxwNkJn9H7wC7TtZz97pMQufAe6BXZkgbpF7cYy5dkMDBihYsnX3eDmVOKGnsJC8xHKUUWxclcsjRJSfdSrMjPSaUmrbegFiAb8fxOvqtNm5d5ZuT3i7m9Z3gxOy6aUUqgzbNzhP13g5lTihp6BxaqmLrwgsDhjIgPTvSYkLpt9qwdPV7O5RZt72ohqy4MAqy5sbyP5Ic/NzS1ChyE8L599FxC8gE0NVnpaatl7wke3JYmxM7tGSGzHOYHamOfSr8fdyhoaOXPeea2LYqbc7MtJfk4OeUUty0IpV9pRYaO6Rryamte4CvPX2U/eebhx4rdQxGO1sOwUFGNufZN3yXMYfZkRYgcx2eK6rFpmHbqjRvhzJpkhwCwM0r0rBpePF4rbdD8QmVzd3c/t97eOJAJb/afXbocefkt/lJ4UOPXbnEvhlVQoT/bLbjSy4kB/8uZ91+pJr8tCjmJ/nujOiLSXIIAItSIlmQFMG/j0pyOFrVym2/3UNjRx9XLk7i7RILzY7+7nMNnRgNiqy4C8nhvWszefqTm8iM852N3/1JbJiJEJPBr1sO55u6KKpq411zZCDaSZJDgLh5RRoHypqpb/fvT2jjqW7t4b2/20eIycgzn9rMF69d6Bist284U9LYSXZcGOagCy+LIKOBNdlx3grZ7ymlSIsJpdaP11d69kg1SsEtK+dOlxJIcggYN69MRWt4PoBbD7tP1dMzMMij96xjflIkS1OjyIkPG/qdlDR2Dg1GC89Jiw6l2k9bDv1WG08fqmLjvHhSoufWVq+SHAJEXmIES1KjeP5Y4CaHt841kR4TynxHAnAO1u8ttdDQ0cv5pi6/23FtLkiLCfHbbqXH9pZR2dzDfZflejuUKZPkEEBuXpHKwfKWgKxasg7aeLvEwqULEkaUEt60PI1Bm+b3b55nYFCTlxg+zlnEbEiNDqWxs8/vNqZq7Ojjly+f5YpFiVyxOMnb4UyZJIcAsmGeve+8qDLwFuI7Vt1GR6+VLfMTRjy+JDWSeQnh/GVvOcBQq0J4TnpMKFrjd+NhP3+pmJ6BQb5587gbY/osSQ4BJD8tGqNB+cXucIM2zacfP8yBsuaJDwbeOmvfuOfi5KCU4qblqfQ4ljXPlW4lj0uN8b+JcMer23iysJJ7t+TM2a5KSQ4BJNRsZGFyJEVVc3+hs7r2Xp4rquEfhVWTOv6tc03kp0URF252ee6mFakAJEYGy2Q3Lxia69DmH8mhu9/Kd7afIC7MzKevWuDtcKZN1gQIMCszotlxog6t9ZyZxj8a56fMQxUtEx7b3W/lUEULHxlj563Fjnkgzp3JhGelRY+cCGezafaWWmjtHsBqs49DXLk4iUgfX/zQZtP883A1P91ZTF17Lz9/z8o5vWCjJIcAsyIjhicOVFLR3E12/NwdfK1usSeHsw2dtPUMjPuJ/53zzQwMai65qEvJSSnFnz6yniAf33zFX4WajcSGmYYS/s93FfObV0tGHPORLfP49i2+23df397L/3mskKNVbazMiObX71vN2py5PT9GkkOAWZkZDUBRVdvcTg7D+qcPV7Rw+aKxq0H2nG3CHGRg3Tgv1nRpNXiVc9Of54/W8ptXS3jPmgw+dmkuRoPiwR2nefpQFV+5fhEhJt/cB+FnO4s5XdfBL967km0r031+l7fJkDGHALMwOZLgIANH53jFUlVLD5HBQRgUQ3sujOWtc02sy4n12TcWYS9nPVbdzpeeKqIgK4Yf3LaMRSmRzE+K4J7NObT1DPCCj87ROdfQydOHqvjQxmxuW53hF4kBJDkEHJPRQH5a1JyvWKpu7WFeYjiLUqI4PM64Q2NHH6frOlyqlIRvSY8Joamzj6jQIB7+wBqCgy4k8k158cxLCOdv71R4McKx/WLXGUJNRj55eZ63Q3ErSQ4BaGVmDMer27EOzt1JR9Ut3aRFh1KQFcORilZsttF3Enu1uAGAS+cnjvq88A0LU+wt2t99cC1JUSOXmVBKcff6TArLWzhT3+GlCEd3vLqN54/V8tFLc4n3s5V7JTkEoJUZMfQMDHK2odPboUyL1prq1h7SY0MpyIqlo8866s9iHbTx8OslLEqOJD8tyguRism6e10W+79xNasyR98l7d1rMjEbDT7XevjZS8XEhJn42KWjV8LNZZIcAtCKDPug9NE52rXU3NVP74CN9JhQCrJjgdFLWrcX1VDa2MXnr1ngN/3A/spgUONWnMWFm7l+WQpPH6qip3/Qg5GN7UBZM68VN/LJrXlzumR1LJIcAlBOfDhRIUFzdjKcs1IpPTaUnPgw4sLNHCofmRwGBm388pWz5KdFcV1+ijfCFG72vg1ZdPRafWbL2xeO1RJqMvKhTTneDmVWSHIIQAaDYkVGzJxdY8k5xyE9JhSlFAVZMRy8qOXwzKEqyi3dfP7qhXN6sp+4YMO8OHITwnn60ORmxc+2cks32fFhhJr9swpOkkOAWpERTXFdB70DvtFEnwpnyyEj1j43YXVWLKWNXbR223d067fa+K9XzrEyM4arlsy91TDF6JxLrO8/34yl0/srC9snkvrvDoGSHALUyswYrDbNa8WN3g5lyqpbewg3G4f6qAuy7OMOhytaqW3r4QfPn6S6tYcvXCOtBn9zXX4KNg27TtZ7NQ6bTc/5VQYmIskhQG3Oiyc3MZxPP36Ivx+o9HY4U1LdYq9Ucr7xr8y0rzZ7/zPH2Pzj3Ty2t5xtq9K4bIHMbfA3+WlRZMaF8uLxOq/GUd/RS7/VRpYf7y0uySFARYaY+OentrAxN56vPH2UH/z7JINjzBXwNdWtPSOWuwgzB3HZggRMQYpPX7mA1798Ob+8a7W0GvyQUooblqXydkkTbT0DXouj3NIN4NfdSrK2UgCLDjXx6D3r+MHzp/j9W+ftn75vXOLtsCZU3drD6qyR9fCP3rveS9EIT7suP4VH3ihl9+l6blud4ZUYKhzJQVoOwm8FGQ08cGs+d67N4A9vneecj0+M6+qz0to9MLQHgAg8qzNjSI4K5sVj3utaKm/uwmhQfn0fSnIQAHzl+sWEmox8/98n0dp3u5eG5jj48YtSjM9gUFyfn8LrZxrp7rd6JYaKZnvXpsnov2+hM/rJlFI/VUqdVkodVUr9Uyk1+tx34fMSIoL57NULeP1MI7tPN3g7nDE55zg4y1hFYLpuWQp9VpvXqu0qLF1+Pd4AM2857AKWaa1XAGeA+2cekvCWD2/OIS8xnO/9+yR9Vt+c/1A11HLw7xemGN/6nDjiws1eq1oqb+726/EGmGFy0Fq/pLV2tuv2Ad4ZHRJuYTIa+PYt+ZRbunl0T5m3wxlVTWsPJqMiKdK/VsAUUxNkNHDzilSeP1rDjuOe3eehrWeA1u4BaTlMwUeAF914PuEFWxcmsmV+PP+7r9wnxx6qW3pIjQ6VhfQEX7thMSszY/jM40d4+1yTx64bCJVKMInkoJR6WSl1fJSvbcOO+QZgBf46znnuU0oVKqUKGxvn3qzcQHLLijSqWno4Vetba+eD6xwHEbjCzEE8es86chLCHPs3e2atsPLmLgCy4vx3djRMIjlora/WWi8b5etZAKXUh4GbgffrcT5qaq0f0Vqv1VqvTUyUjVd82dVLk1EKdp7w7izU0ThnRwsBEBNm5rGPbCA23Mw9jx6go3f2J8ZVNDtaDtKtNDal1PXAV4Fbtdbd7glJeFtCRDBrs2N5ycvr11ys32qjvqNXWg5ihJToEB66cxXNXf28cmr2K+0qLN0kRJiJCPbvOcQzHXP4NRAJ7FJKHVFKPeyGmIQPuHZpCqdq26ls9p2cX9fWi9ZIy0G4WJsdS3JUMC8cm/3B6XKL/1cqwcyrleZrrTO11qscX59wV2DCu67NTwZ8q2vJ2debIS0HcRGDwb7m0mtnGunsm92Jcf6+GquT/07vEzOSHR/O4pRIn+paOl7dDsBS2Q9ajOKGZSn0W228OouTOPusg9S09ZApLQcRyK5dmkxhmW9srAJwrLqVrLgwYsLM3g5F+KC1OXEkRgbz4izOe6hq6UFryJbkIALZtY6NVTwxyDcZR6vaWJ4R7e0whI8yOtZc2n26YdbWXHJWKvn7BDiQ5CDGkZ8WRXpMKC+d9P64g6Wzj6qWHlZKchDjuGF5Cr0D7l1zqbK5e2gL2qEJcAGQHPy7FkvMiFKKa5Ym88SBCvqtNsxB3vsscay6DYDl6bK2oxjbhnnxxIebeeFYLTcuT53x+Q6UNfP+37+DAm5ekUZbzwBhZiOJEf6/fIu0HMS41uXE0Ttg41Rtu1fjOFZlTw7L0mUwWozNaFBct8zetdQ7MLPFI8/Ud/DRPx0gIzaUO9Zk8OLxWl4+VU9WXFhA7DIoLQcxroJs+yf1QxUtrMz03qf2oqo2chPDiQwxeS0GMTfctDyVv71TwUO7znD/DYun9UZe09rDh/+4nxCTkT/fu57MuDDuv2ExzxXVBsQcB5CWg5hAanQoadEhHCxv8Wocx6pbWZEu4w1iYpvz4vnAxiweeaOUX+w6M+Xv7x0Y5J5H99PZa+VPjsQA9n3X37chi0sWJLg7ZJ8kLQcxodXZsRyu8MyiZqOpb++lvr2PFRky3iAmppTie7cuY8Cq+a/d5wgyGvjMVQsm/f17Sy2cqe/kt+8vCOg5NdJyEBNakxVLdWsPtW09Xrm+c7xhhVQqiUkyGBT/eftybi9I56FdZ/jSU0WTvn/3lVgwGRVXLEqa5Sh9myQHMaGC7FgADpV7p/VwtKoVg5KZ0WJqDAbFT9+9ko9flsv2IzVc/tPX+M8XTtHWPf7KrXtLLazOjCXUbPRQpL5JkoOY0NLUKIKDDByq8M64w9HqNhYkRRJmll5QMTVGg+L+G5fwyhe3ctOKVB55s5Q7Hn6b9jGW9m7rGeB4dRsb8+I9HKnvkeQgJmQOMrAiI9org9Jaa47JzGgxQ5lxYTx05yr+96MbKGvq4jOPH2bQ5rr9zP7zzdi0fVA70ElyEJNSkB3LiZq2GdeOT1VNWy+Wrn6ZGS3cYsv8BL63bRmvFTfy4xdPuTy/t8RCcJCB1VlS/CDJQUzKmqxYBgY1xx0zlT3laKV9nGO5VCoJN3nfhizu2ZzD/7x5nr8XVo54bm+phTXZsQQHBfZ4A0hyEJPkHJT2ZNfSwKCNh98oJSbMxOKUSI9dV/i/b960hEvmJ/DNfx2n3GLfJ6Slq59Tte1sypUuJZDkICYpISKY7Pgwjw5K/3r3OYoqW/nhu5YTYpJPcsJ9gowGfvaelZgMiu9sP4HWmn2lFgA2yXgDIMlBTEFBViwHy1vR2nUgz90Olrfwq91nub0gnZtWzHwBNSEulhIdwheuXcRrxY3sPFHH3lILYWajTLZ0kOQgJq0gO5amzj5KGjvdfm7roI3ufis2m6azz8rnnzxCWkwo37013+3XEsLpw5uyWZIaxXefO8kbZxpZmxPn1dWHfYn8FsSkXbc0GaNB8VRhlVvPa7Npbvqvt1j67Z3kfv0FCr63i6qWbh66c5UstCdmVZDRwA/etYzatl7KLN0y3jCMzCoSk5YUFcLVS5J46mAVX7h2odsqOvadt1Bc38F712aSHB1CT7+VNdlxrJ8X55bzCzGeNdmx3L0+k8f3V8r8hmEkOYgpuXt9FjtP1LPrZD03r0hzyzmfKqwiMiSI727Ll4Fn4RXfunkpWxcmyfpdw0i3kpiSyxYkkh4Tyt/eqXDL+dp7B3jhWC3bVqVJYhBeE2YO4vplKQGxic9kSXIQU2IwKO5en8nbJRbON3XN+HzPFdXQZ7Vx59pMN0QnhHAXSQ5iyt6zNhOjQfHEgZm3Hv5+oJLFKZEsl418hPApkhzElCVHhXDV4iT+UVhFv9U27fOcrmunqKqN96zNlOa8ED5GkoOYlvdtyMLS1c+9f9rPW2ebRkyMm+wkuacKqzAZFbetTp+tMIUQ0yTVSmJati5M5Bs3LuGRN0v5wB/eYWlqFGkxIZRbuqlo7mZNdix/vGfdmIPMvQOD/PNwNdcsTSYu3Ozh6IUQE5GWg5gWpRT/57Jc3vrqFTx4xwoMBqhs7iEnIZzbC9LZW2rhS08VjdmKePj1Epq7+vnQphzPBi6EmBRpOYgZCQ4ycue6TO5cN7LaKDs+nB+/eJq8xAg+f83CEc+VW7r47Wsl3LIyjY0yI1UIn+SWloNS6ktKKa2USnDH+cTc9/HLcrmjIINfvnKWZ49UDz2uteaB7ScwGw1886YlXoxQCDGeGbcclFKZwDWAe2ZFCb+glOJHty+jsrmbzz95hL0lFr5w7UKOVLTyanEj37xpCclRId4OUwgxBnd0K/0C+ArwrBvOJfxIcJCR39+zll++fJY/v13Gc0U1hJiMLEqO5MObc7wdnhBiHDPqVlJK3QpUa62L3BSP8DNRISa+dfNSdn1hK5csSKC9d4Dvv2sZJqPUQgjhyyZsOSilXgZSRnnqG8DXgWsncyGl1H3AfQBZWVlTCFH4g3kJ4fzug2vpHRiUNZSEmAPUdHf1UkotB14Buh0PZQA1wHqtdd1437t27VpdWFg4resKIUSgUkod1Fqv9cS1pj3moLU+BiQ5/62UKgPWaq2b3BCXEEIIL5KOXyGEEC7cNglOa53jrnMJIYTwLmk5CCGEcCHJQQghhAtJDkIIIVxIchBCCOFCkoMQQggXkhyEEEK4kOQghBDChSQHIYQQLiQ5CCGEcCHJQQghhAtJDkIIIVxIchBCCOFCkoMQQggXkhyEEEK4kOQghBDChSQHIYQQLiQ5CCGEcCHJQQghhAtJDkIIIVxIchBCCOFCaa09f1GlOoBij194ahKAJm8HMQkSp3tJnO4zF2J0miuxLtJaR3riQkGeuMgoirXWa7107UlRShX6eowgcbqbxOk+cyFGp7kSq1Kq0FPXkm4lIYQQLiQ5CCGEcOGt5PCIl647FXMhRpA43U3idJ+5EKPTXInVY3F6ZUBaCCGEb5NuJSGEEK601hN+AZnAq8Ap4ATwWcfjccAu4Kzjz1jH4/GO4zuBX190rruBY8BRYAeQMMY11ziOOwf8FxdaOZcBhwAr8G5fjHHY8+8GNLDWF+MEfgEccXydAVq9HOcPgUqg86LHg4EnHfG/A+TMUpzvdcR4AnhwnNeDN+/NGcXowXtzpr9LT92bk41zyvfmNGO9Bjjo+J0cBK6c7P/pdO/PMX/miQ5wnDQVKHD8PdLxn7UUeBD4muPxrwE/cfw9HLgE+MTw/wzspbMNON4cHN//wBjX3A9sAhTwInCD4/EcYAXwGCNfgD4T47AY3gD2MfIF6FNxDjvm08AfvRznRsd1L34Bfgp42PH3u4AnZyHOeKACSHT8+8/AVT52b844Rg/dm26J0wP35lTinPK9Oc1YVwNpjr8vA6qn8ruazv051tekupW01rVa60OOv3dgz4LpwDbHL9T5i32X45gurfVbQO9Fp1KOr3CllAKigJqLr6eUSgWitNZ7tf2nemzYucu01kcBm6/G6PB97DfAiPP7YJxOdwOPeytOxzn2aa1rR3lq+DX/AVzlOJc748wFzmitGx3/fhm44+JAvHxvzjhGh9m+N90Vp9Ns3ZuTitNxjinfm9OM9bDW2vn6OAGEKKWCJ/u7ms79OZYpjzkopXKwZ7d3gGTnL8zxZ9J436u1HgA+ib3JU4M9g/5hlEPTgaph/65yPDYnYlRKrQYytdb/9uU4h8WRDcwDdnsxzvGkY2/So7W2Am3YP/W5LU7sTfDFSqkcpVQQ9hdU5hixeOXedEeMnrg33RHnsDhm7d6cQpzjmdS9Oc1Y7wAOa637mPx9N6P7c7gpJQelVATwNPA5rXX7VC+mlDJhf6NYDaRh7+u7f7RDR3lMz4UYlVIG7P2lX/TlOC/6913AP7TWg16Mc9zTjPLYiJ9hpnFqrVsccT4JvAmUYe+bnXIsY/F2jJ66N938u5y1e3MKcY5nUvfDVGNVSuUDPwE+PpXrTOG4CU06OThe5E8Df9VaP+N4uN7RjHE2ZxomOM0qAK11iaPJ83dgs1LKqJQ64vj6HvZslzHs+zIYoyvCB2OMxN5X+JpSqgx7X+V2pdTQ1HwfiXO4uxjWbPdSnOOpwvGJzvEJLxpodnOcaK2f01pv0Fpvwr7211kfuzdnGqOn7k13/i5n896cbJzjGffenE6sSqkM4J/Ah7TWJcOu4/K7ctf9OZpJJQdHH9ofgFNa64eGPbUd+LDj7x8Gnp3gVNXAUqVUouPf1zjOOai1XuX4+rajmdWhlNrouPaHJjq3r8SotW7TWidorXO01jnYB/1u1VoX+lKcw35vi4BYYO/wk3s6zgnOMfya7wZ2OxKNO+NEKZXk+DMW+0Dj733s3pxRjB68N93yu/TAvTmpOCc4xZj35nRiVUrFAM8D92ut9zgPHuf/dMb355j0BCPWjp/zEuxNk6NcKC+7EXvf2ivYy7FeAeKGfU8Z9gzaiT2bLXU8/gnsgzJHgeeA+DGuuRY4DpQAv+ZCOdY6x/m6AAtwwtdivOiY1xhZEeJTcQIPAD/2kf/zBx3fZ3P8+YDj8RDgKex9xPuB3FmK83HgpOPrrnFeD968N2cUowfvzRnHiWfuzcnGOeV7czqxAt903D9Hhn0lTfb/dDr351hfMkNaCCGEC5khLYQQwoUkByGEEC4kOQghhHAhyUEIIYQLSQ5CCCFcSHIQQgjhQpKDEEIIF5IchBBCuPj/VbYSepA4PugAAAAASUVORK5CYII=\n", @@ -723,10 +748,31 @@ "needs_background": "light" }, "output_type": "display_data" - }, + } + ], + "source": [ + "fig, ax = plt.subplots()\n", + "ax.plot(raw_data[\"convert\"],raw_data[\"Concentration\"]-raw_data[\"smooth\"])\n", + "ax.set_xlim(\"2018-01-01\",'2020-01-01')\n", + "\n", + "plt.show()\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Nous rappliquons un filtre à ces nouvelles données pour avoir un signal plus facile à visualiser." + ] + }, + { + "cell_type": "code", + "execution_count": 32, + "metadata": {}, + "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAD8CAYAAACcjGjIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xl8VNX9//HXyb5vkITsIYSw7wFERBTEDRQQpOBaq19r1VZrv99fN9vaVttvrdUutvql7pUqCC4U3EBAVLYs7AkJEBKyL2RPyDrn90eGFhmWLDdzZyaf5+ORB0mYnPNmuJPP3HvOPUdprRFCCCHO5mZ2ACGEEI5HioMQQggbUhyEEELYkOIghBDChhQHIYQQNqQ4CCGEsCHFQQghhA0pDkIIIWxIcRBCCGHDw4xOBw8erBMTE83oWgghnFZGRkaV1jrcHn2ZUhwSExNJT083o2shhHBaSqkCe/Ull5WEEELYkOIghBDChhQHIYQQNqQ4CCGEsCHFQQghhA0pDkIIIWxIcRBCCGHDsOKglHJXSu1VSm0wqk0hhBDmMPLM4REg28D2hBBCmMSQ4qCUigXmAy8Z0Z4QQghzGXXm8Efg/wEWg9oTQghhoj4XB6XUAqBCa51xicfdr5RKV0qlV1ZW9rVbIYQQ/ciIM4eZwM1KqXzgbWCOUurNcx+ktV6ptU7VWqeGh9tlUUEhhBC91OfioLX+sdY6VmudCCwHtmit7+hzMiGEEKaR+xyEEELYMHQ/B631NmCbkW0KIYSwPzlzEEIIYUOKgxBCCBtSHIQQQtiQ4iCEEMKGFAchhBA2pDgIIYSwIcVBCCGEDSkOQgghbEhxEEIIYUOKgxBCCBtSHIQQQtiQ4iCEEMKGFAchhBA2pDgIIYSwIcVBCCGEDSkOQgghbEhxEEIIYUOKgxBCCBtSHIQQQtiQ4iCEEMKGFAchhBA2+lwclFI+Sqk9Sqn9SqnDSqlfGhFMCCGEeTwMaKMVmKO1blRKeQJfKqU+0lrvMqBtIYQQJuhzcdBaa6DR+qWn9UP3tV0hhBDmMWTMQSnlrpTaB1QAm7TWu41oVwghhDkMKQ5a606t9UQgFpimlBp77mOUUvcrpdKVUumVlZVGdCuEEKKfGDpbSWtdC2wDrj/P363UWqdqrVPDw8ON7FYIIYTBjJitFK6UCrF+7gtcAxzpa7tCCCHMY8RspSjgdaWUO13FZo3WeoMB7QohhDCJEbOVDgCTDMgihBDCQcgd0kIIIWxIcRBCCGFDioMQQggbUhyEEELYkOIghBDChhQHIYQQNqQ4CCGEsCHFQQghhA0pDkIIIWxIcRBCCGFDioMQQggbUhyEEELYMGJVViGECSwWTd3pdmpPt1N3up0wPy/iB/mZHUu4CCkOQjiRxtYOtudWsjmrnC05FdQ2t//775SCmydE88jc4SSFB5iYUrgCKQ5COIFOi+b1Hfn84dMcmto6CfHzZM7ICMZGBxPi50mwrydp+TW8viOfDQdKWTo5lp/dNJoAb3mJi96RI0cIB3ewqI6fvHeQg8V1zE4J58GrhjElIRQP968PGc4dFcm9VwzlhW3HeX1nPvsKa3np7lTiwuRSk+g5pbW2e6epqak6PT3d7v0K4Ww+PFjKw//MZFCAN7+4aTTzx0WhlLrkz31xtJKHVmXi4e7GC7dPZnrSIDukFf1NKZWhtU61R18yW0kIB5WeX82jq/cxKT6UzY/NZsH46G4VBoBZw8N5/6GZhPh6cvtLu9l4oLSf0wpXI8XhIlraO8koqCGjoJrG1g6z44gBJK+ykfveSCcmxJe/35VKsK9nj9tICg/gvQdnMjEuhO+9vZcNB0r6IalwVTLmcI6apjb+b3seXx2rIru0ng7Lfy67xYf5MX1oGN+fl0J0iK+JKYUrO9XYyjdfTcNdKV67Zyph/l69bivYz5PXvjWNb72axiNv78Oiu2Y0CXEpUhysOjotrNp9kmc35dLY2sG0xDDuvzKJiXEhuLspskvrySqtZ/3+EjYcKOXhOcncN2so3h7uZkcXLkRrzQ/XHaCsvoW377+MhEH+fW4zwNuDV++Zyj2vpfHo23vxcFPcOC7KgLTClfW5OCil4oA3gCGABViptf5TX9u1p+La09z7WhpHyhq4fNggfnHTGEYMCfzaY+aOigSgsLqZJzdm8ftPclibUcTLd6fKnHJhmPf3FbM5u4LH549icnyoYe36e3vw2j1TuevlPTy6eh/hgd5MTQwzrH3hevo8W0kpFQVEaa0zlVKBQAawSGuddaGfcaTZSjVNbSx9cQcVDa38ful4rhszpFuDfp/nVvLY6n0opVh133SbYiJET1XUtzDvue0kRwSw5tszcHfr3uBzT9Q0tbHkhR1UN7ex7juXM0ze2DgVp5qtpLUu1VpnWj9vALKBmL62aw/NbR3c81oahTWneemuVK4f271pggCzU8JZ/e3LcHeD5St3cqi4rp/TClemteYn7x2ipb2Tp5eO75fCABDq78Wr90zFXSnueTWNqsbWfulHOD9DZysppRKBScDu8/zd/UqpdKVUemVlpZHd9kp7p4WHVmVyoKiWPy+f2Kt54MkRgaz59gz8vDxYsXIXWSX1/ZBUDATr95ewObuc/752RL+/m08Y5M9Ld6dS0dDCg6sy6bTY/14n4fgMKw5KqQBgHfCo1trmt6TWeqXWOlVrnRoeHm5Ut732zKc5bM2p5NeLxnL92N4PziUM8mfNAzPw83bnu29lcrqt08CUYiBoau3gyY3ZTIgL4VtXDLVLn5PiQ3lq0Tj2nKhm5fY8u/QpnIshxUEp5UlXYViltX7XiDb70+68U6zcnseKaXHcPj2hz+3FhPjy7LKJHK9s4smNFxxqEeK8Vm7Po7KhlZ8vGN1vl5PO55bJMdw4bgjPbsqRy6LCRp+Lg+q6SP8ykK21frbvkfpXQ0s7P3hnP/Fhfjw+f7Rh7c5MHsz9VyaxavdJNmWVG9aucG3l9S2s3J7H/PFRTEkwbnZSdyileGrROEL9vPj+6n20tMtZr/gPI84cZgJ3AnOUUvusHzca0G6/+PWGLEpqT/Pssgn4G7xi5Q+uTWFMdBA/XHeAivoWQ9sWrukPn+bQadH88LqRpvQf6u/FM7dO4GhFI7/7+IgpGYRjMmK20pdaa6W1Hq+1nmj9+NCIcEbblFXOmvQiHpg9jCkJxs/x9vZw50/LJ9Hc1sGjq/fJQJ+4qKySet7JKOLuyxNM3aTnypRw7p6RwGs78sk8WWNaDuFYBszaSs1tHfzig0OMHBLIo9ek9Fs/yREBPLloHDuOn+LZTTn91o9wfr/9KJtgX08evnq42VH4n+tHEhnow0/fO0RHp8XsOMIBDJji8OK245TUtfCrhWPx8ujff/bSKbGsmBbHX7ceZ7OMP4jzSM+v5oujVTx0VTLBfj1fVM9oAd4ePHHzaLJL63ltR77ZcYQDGBDF4eSpZl7cnsfCidFMG2qfJQN+cdMYxsYE8diafZw81WyXPoXz+POWYwwO8OKOy/o+W84o140ZwtyRETy7KZfi2tNmxxEmGxDF4cmNWXi4KX58wyi79enj6c4Lt08B4OG3MmmXU3VhtfdkDdtzK7lvVhK+Xo6zcKNSiiduHoNFa365/rDZcYTJXL44bM+t5NOsch66OpkhwT527TsuzI/fLRnPgaI6/rr1mF37Fo7rL1uOEernyZ0OdNZwRlyYH49ek8KnWeVySXSAc+niYLFontqYTcIgP+6bZZ87T891w7goFk2M5vktxzhYJDcaDXQHi+rYcqSC+2YlGT6V2ij3XjGU4REB/GL9YbnjfwBz6eKwKbucnPIGvn9Niqn7Lvzy5rEMCvDisTVyo9FA9+ctRwny8eCuGY531nCGp7sbTy4aS3HtaZ7fetTsOMKqryto95TLFgetNX/beoz4MD8WjDd3Y5NgP09+t2Q8RysaeXZTrqlZhHlyyhrYlFXOt64YSqCP+TOULmZ60iBumRzDyu15HKtoMDuOAMrsfGOtyxaHr46dYn9RHQ/MHoaHu/n/zKtGRHDb9Hj+/kUeh0vk8tJA9Pcv8vD1dOfuGYlmR+mWn9w4Cl9Pd372/mG7v2sVto6WN9q1P/N/a/aTv249RkSgN0umOM7WEj+8fiShfl78cn2WvNgGmLK6Fj7YV8w3psYR2oc9oe1pcIA3/+/6kezMO8W7mcVmxxnwjlZIceizjIIaduad4v4rkxxqj+dgX0/++9oR7MmvZsOBUrPjCDt6bUc+nRbNt2aaMzGit1ZMi2dKQii/2pAl64WZzN6X91yyOLyw7Rghfp6smBZvdhQb35gax5joIH77YbbMBBkgGls7WLW7gBvGRpm6hlJvuLspfrdkPKfbO3n8/UNyxmsiuazUR8crG9mcXcE3L090yKmC7m6KX9w0hpK6Fl74/LjZcYQdvL3nJA0tHaZNp+6r5IgAHpvXde+DnPGaQ2stl5X66o0d+Xi5uznUsgTnmjY0jJsmRPN/nx+nRJYpcGntnRZe/SqfaYlhTIq3734NRrrviqFMiA3mF+sPc0r2nba7ysZW6k6327VPlyoODS3trM0oYsH4KAYHeJsd56J+dMNItIbn5c5pl/bxoTKKa0/zX1cmmR2lTzzc3Xh66QQaWzr4zYey74O92fuSErhYcXhvbzFNbZ3cdXmi2VEuKSbEl+XT4liTVkhhtSzM56r+sauAuDBf5o6MMDtKn40YEsg9MxN5d2+RbCtqZ0fL7X+vicsUB601r+/IZ0JcCBPjQsyO0y0PXpWMm5vi+S1y9uCKcsoa2HOimjumJ+Bmx72h+9ODVycT4uvJUxuzZXDajo5WNBLkY98xVJcpDl8dO8XxyibuduBlCc41JNiH26bFszaziIJTTWbHEQZ7c1cBXh5u3JoaZ3YUwwT7evLoNSnszDvFliMVZscZMI5WNJISGWjXPl2mOLy2I59B/l7MN3mpjJ568KpheLgp/iJnDy6lsbWDdzO7xr/CnOSmt+66bXo8SYP9+c2H2bIUvZ0cq2hkeGSAXft0ieJQWN3MZ0fKWTEt3qFueuuOiCAf7rgsgXczi8irtP+gk+gfZ8a/HHFZ7r7ydHfjRzeM5HhlE2/vOWl2HJd3qrGV6qY2kiOc8MxBKfWKUqpCKXXIiPZ6ak16IYqudzTO6IHZw/D2cOe5zbICpivQWvPmzgLGxgQ5zfhXT80bHcn0oWH86bNjcjNnPztzf8PwCOc8c3gNuN6gtnqko9PCO+lFzE4JJzrE14wIfRYe6M29VwzlX/tLZBaIC0jLryGnvIE7L0tAKdcYiD6XUor/uW4EVY2tvLEz3+w4Lu3MTCWnvKyktd4OVBvRVk9tP1pJWX0L35jqnGcNZ9w/O4kQP09+/0mO2VFEH/1zdwGBPh7cNCHa7Cj9KjUxjNkp4bz4+XEaWzvMjuOyjlY0EuDtwZAg++5k6fRjDm/vKWRwgDdzRzn3PPIgH08evGoYn+dWsivvlNlxRC/VnW7no0NlLJwYjZ+X4y3fYrTH5qVQ09zOq1+eMDuKyzpa3khyRIDdz0LtVhyUUvcrpdKVUumVlZWGtFlR38JnRypYOiUWTwfYs6Gv7pqRyJAgH57++IjMIXdS6/eX0NphYZkLTV+9mAlxIVwzKpKVX+RR12zf5R0GAq01OeUNpNj5khLYsThorVdqrVO11qnh4eGGtLk2s4hOi+YbU13jhejj6c6j1wwn82Qtm2Rzd6f0TnohI4cEMi4m2OwodvPYvBQaWjp46cs8s6O4nIqGrplKo6KC7N63077d1lqzOq2Q6UPDGDrY3+w4hlk6JZZh4V1zyGW/aeeSXVrPgaI6vjE1zmUHos9ndHQQN44bwitfnpBF+QyWVVIPwGhnLQ5KqbeAncAIpVSRUupeI9q9mF151RScamb5NNc4azjDw92NJ24eQ/6pZlZul3dizmR1WiFe7m4smug4uw/ay/evSaG5vZMXZRl6Q2WVdhWHkc5aHLTWK7TWUVprT611rNb6ZSPavZjVaScJ9PHghrHOdUd0d8waHs78cVH8desxWZTPSbR2dPL+vmLmjYl0mm1AjTQ8MpDFE2N4Y2cB5bJjnGGyS+uJDfUl2NfT7n075WWluuZ2PjxUxuJJMfh4Otcd0d31+IJRuLspnlh/2Owoohs2ZZVT29zONwbIQPT5PHpNCp0WLQtJGiirtN6U8QZw0uLw/r5i2josLjMQfT5Rwb48es1wPjtSIYPTTmBNehHRwT7MTB5sdhTTxA/yY9nUON5OOylnvAY43dZJflWTFIfu0lrz1p6TjIsJZky0a88IuWfmUFIiA/jVhsO0dsjgtKMqr2/hy6OVLJkSi7uLLM3dW9+dk4xSij99JkvB9FVOeQMWbc5gNDhhcThYXMeRsgaXPms4w9Pdjcfnj6aw+jRv7pIFzhzVe3uLsWi4ZXKs2VFMFxXsy53WhSTN2KDGlZg5UwmcsDi8nVaIj6cbN0907aUJzrgyJZwrkgfz/Jaj1LfITUaORmvNuowipiSEutSU6r546Opk/L09eHJjttlRnFp2aT2B3h7EhpqzZpxTFYfmtg7W7yth/rhognzsP3pvlh/dMJKa5nZe3CbTBB3NweI6jlY0csvkgTd99ULC/L14ZO5wPs+tZGuObAjUW1ml9YyMCjRtF0GnKg4bD5TS2Nrhcvc2XMrYmGAWTYzmla9OUFYn0wQdybuZxXh5uLFg/MA4k+2uu2YkkjjIj6c2yoZAvWGxaI6YOFMJnKw4rE4rJGmwP6kJoWZHsbsfXDsCiwWe25RrdhRh1dZh4YN9xcwbHWnKPHRH5uXhxk/nj+ZYRSP/3C3jZT1VWNNMU1unaeMN4ETF4XhlI+kFNSwbYEsTnBEX5sedMxJ4J6OQYxUy0OcItuZUUNPczlIZiD6va0ZFMDN5EM9tzqW2uc3sOE7lzGC0nDl0w5r0Qtzd1IC+tvvQ1cn4errzrJw9OIR1GUWEB3oza/jAvbfhYpRS/GzBaOpPt8se6T2UXVqPm4IRQ+y7NejZnKI4tHdaWJdRzNUjIogItO+GF44kzN+L+2Yl8eHBMg4WyY5xZqpuamNrTgULJ0Tj4QLLxfeXkUOCuHVKHP/YWSA3xvVAVmk9SeEBpq4A4RRH9bacSqoaW1mWKqfv980a2rVj3KeyY5yZ1u8rpr1Ts2SKHJOX8v15Kbi5wTNyzHZbdmmDqZeUwEmKw5r0rt3erh7p3Lu9GSHQumPcdtkxzlTrMosZHRVk+gvYGQwJ9uHeK4bywb4SOePthtrmNoprT5s6GA1OUBwqGlrYcqSCJZNjXGK3NyPcNSORyCBvnvkkR3aMM0FueQMHi+vkrKEHvj17GGH+Xvz2o2w5Zi9hv7WATog1d3kgh/9t+25mMZ0Wza0DeLXLc/l4uvPdOcNJL6hhx3E5e7C3dZlFeLgpFg6Qu/SNEOTjyXfnJLPj+Cm25RqzTbCr2l9Yi1IwVorDhWmteSe9kCkJoSRH2H8PVUe2dEosoX6evLmrwOwoA0qnRfP+3mKuGhHO4ABvs+M4ldunJxAb6it3+l/C/sJahoUHmL4KhEMXh/1FdRyvbOJWOX234ePpztIpsWzKKqdCNlexmy+PVVFe38oSubehx7w83FgxLZ7dJ6rJr2oyO45D0lqzv6iWCbEhZkdx7OLwbmYRXh5u3Dje9XZ7M8KKafF0WDRr0gvNjjJgrMsoItjXkzmjZHJEbyyZHIubQo7ZCyipa6GqsY0JceZvR+CwxaGtw8L6/SVcOzrS9NMrR5UUHsDM5EG8taeQTosM8vW3hpZ2Pjlcxk0TovD2cM0dCPvbkGAfrh4RwdqMIjpkzSUb+wtrAeTM4WK25lRQ29wuM0Iu4fbpCRTXnma7DPL1u08Pl9PaYWHxpIF7l74Rlk2No6Khlc/lmLWxv7AWL3c3RkaZd2f0GQ5bHN7NLGJwgDezBvC2i90xb3Qk4YHerNotA9P9bcOBEmJCfJkcP/AWfjTSnJERDA7w5u00ubR0rn2FtYyKDnKIM1NDioNS6nqlVI5S6phS6kd9ba+mqY0tRypYNFGWJrgUT3c3lqXGsuVIBSW1p82O47Lqmtv54mgV88dHDciFH43k6e7GkikxbDlSQUWDTKY4o9OiOVhcx0STp7Ce0effvEopd+CvwA3AaGCFUmp0X9r814ESWZqgB5ZPjUcDr+/INzuKy/rkcBkdFs0CmRxhiGWpcXRaNO9mFpsdxWEcq2ikua2TCXHmjzeAMWcO04BjWus8rXUb8DawsC8NrsssZpQsTdBtcWF+LJoYw+s78+WdWD/514ES4sP8GBfjGO/qnN2w8ACmJoayJq1Q7pi2+vdgtAsVhxjg7IuHRdbv9crxykb2F9Zyiwz69cgjc4fT3ql5QW4wMtypxlZ2HD/FArmkZKjlU+PJq2pi94lqs6M4hP1FtQR6ezB0kGPsRW5EcTjfq8XmrYBS6n6lVLpSKr2y8sKzFN7LLMZNwcJJsjRBTyQO9mfJ5BhW7T5JaZ2MPRjp48NldFq0bAVqsBvHRRHo48Hbe2SnOOgqDuPjgk3bM/pcRhSHIuDshY9igZJzH6S1Xqm1TtVap4aHh5+3IYtF897eYmYNDx/Q+zb01nfnDEdrzV+3ysYqRtqwv5SkcH9GOcD0Qlfi6+XO4kkxfHiobMDvFNfS3smR0gaHuL/hDCOKQxowXCk1VCnlBSwH1vemod0nqimuPT2gd3vri7gwP5alxrE6rVA2VjFIRUMLu0+cYsH4aLmk1A+WT42nrcMy4AemD5fU0WHRjHel4qC17gAeBj4BsoE1WuvDvWnrvb1FBHh7cO3oIX2NNWA9PCcZpRR//uyo2VFcwieHyrBouElmKfWL0dFBTIgN5u20kwN6YDotvwaAKQmOcw+NITcRaK0/1FqnaK2Haa2f6k0bp9s6+fBgGTeMHYKvl/k3gDirqGBf7piewLrMIo5VNJodx+l9fLiMYeH+DI+US0r9Zfm0eHLLG8k8WWt2FNPsOVFNUrg/4YGOs9Kvw9xh9mlWGY2tHdwiq1322UNXD8PX053nNuWaHcWp1TS1sSuvmuvHyplsf7ppQjR+Xu68NUAHpjstmrT8aqYPDTM7ytc4THF4b28xMSG+DvcEOaNBAd7ce8VQNh4s5VCxbMvYW5uzy+m0aK4fI5eU+lOAtwc3T4hmw4ES6k63mx3H7o6U1dPQ0sE0B/vd5xDFoaKhhe25lSyaFO0w07ic3X1XJhHi58nTn8im7r31yeFyYkJ8GRsjN2P2t9unJ9DSbuH9vQNvYHqP9T6PaUMHmZzk6xyiOKzfV4JFw+JJcknJKEE+njx41TC251ayK0+2Eu2pptYOth+t5NoxkTJLyQ7GxQYzPjaYVbsLBtzAdFp+NTEhvsSE+Jod5WscojiszShiQlyIbAVqsLtmJBIZ5M3vP8kZcC+4vtqWU0lbh4Xrx8h4g73cPr1rYDq9oMbsKHajtWbPCccbbwAHKA5ZJfUcKWtgqdzbYDgfT3e+O2c4GQU1sql7D318uIxB/l6kJjrei9ZV3TQhmkAfD1YNoH3R86qaqGpsc7jxBnCA4rAuswhPdyVLE/STZalxxIb68uynuXL20E2tHZ1sPVLBvNGRuMsYmN34eXmwZHIsHx4so7ppYNwxfWa8YaoUh69r77Twwb5i5o6MJNTfy8woLsvLw41H5g7nYHEdn2aVmx3HKew4dorG1g6ukymsdnfb9HjaOi2szRgYGwHtOVHN4AAvkgY7xmJ7ZzO1OGzPraSqsU32behniyfFkDTYn2c/zcUie01f0oYDpQR6e3D5MMeaPTIQpEQGMi0xjH/uPjkgjtU9J6qZNjTMISc9mFoc3s0sJszfi6tGnH8hPmEMD3c3HrlmODnlDWw8WGp2HIfW1NrBR4dKmT8+yiG2ahyIbr8snvxTzXx5rMrsKP2qqKaZ4trTTHPQcS3TikNdczubssq5eUI0nrIVaL+7aXw0IyIDeW5zLh2dFrPjOKyPDpXR3NbJUjmbNc31Y4cwOMCLN3bmmx2lXznq/Q1nmPZb+V8HSmjrtMiL0E7c3BQ/uDaFvMqmAbtMQXeszSgkcZCfQy2ANtB4e7hz27R4PjtSwclTrru68JfHqgj182TEEMdct8u04vBORhEjhwQyJlruPrWXeaMjmZE0iD9syh3w6+efT2F1M7vyqlkyOdYhrwEPJLdNT8BdKf6xK9/sKP3CYtFsz61k1vBwh50RZ0pxaGm3sL+wlqVT5EVoT0opfn7TaOpPt/PHzbKk97nezSxGKbhFzmZNNyTYh+vGDmF1WiGn2zrNjmO4rNJ6qhrbmJ3iuOOtphSHmuY2PNwUi2WfaLsbFRXEimnx/GNXAUfLG8yO4zC01qzLLOLyYYMcbhmDgeqblydS39LB+/tcb72lz603pV4pxeHrapvbmDsqgkEBjrN2+UDy2LwU/Lzc+fXGbLkxziotv4aT1c0skSXjHUZqQiijooJ4fUe+yx2n23IqGBsT5FD7N5zLlOLQYdEsS4279ANFvxgU4M2j16SwPbdSboyzWptRiL+Xu+zd4ECUUnzz8gSOlDWw2zqzxxXUnW4n82StQ19SApOKg4ebcvgnxtXdNSOBUVFB/Oz9Q9Q1D7w19M/W1NrBxgNd9zb4eXmYHUecZeHEGEL9PHn5yxNmRzHMjmNVdFo0V42IMDvKRZlSHEL9vPCQextM5enuxu+XjudUUxtPbswyO46pPjxYSlNbp5zNOiAfT3fuuCyBzdnlnKhqMjuOIbblVBLo48GkuBCzo1yUOcVB1lFyCGNjgnlgdhLvZBT9e4BsIHonvYikwf5yb4ODunNGAp5ubrziAmcPWms+z61k1vDBDv8G2ZR03h6O/aQMJN+bO5zkiAB+vO4ADS0D7/JSXmUje/KrWZoq06odVUSgD4smRfNORiE1Tr5aa255I2X1LU5xWb1Pv6WVUrcqpQ4rpSxKqVSjQgn78fZw5/dLx1NW38IzA3BL0bUZRbgpZJaSg7tvVhIt7RZW7XbuvR625VQAjj2F9Yy+voU/BNwCbDcgizDJpPhQ7pqRyD92FXCwqM7sOHbT0WlhXWYRV42IIDLIx+w44iJSIgOZnRLO6ztrXwOSAAART0lEQVQLaO1w3pviPsuuYOSQQKKCHf9emj4VB611ttZ64L3ddEGPXZvCoABvHn//IJ0DYKlkgC+OVlFe38qyVDlrcAb/NSuJyoZWPthXYnaUXimvbyGtoJobxkaZHaVb7HbxXyl1v1IqXSmVXlk5cAc/HVWQjyePzx/F/qI6/jlAFuZ7J6OQMH8v5oyMNDuK6IaZyYMYFRXE/31+3Cn3evjoYClaw/zxznEvzSWLg1Jqs1Lq0Hk+FvakI631Sq11qtY6NTzc8a+3DUQ3T4hmZvIgnv74CJUNrWbH6VcVDS1syipn0cQYvGSChFNQSvHgVcM4XtnEx4fLzI7TYxsPljJySCDJEY65Cuu5Lvmq0Fpfo7Uee56PD+wRUNiPUopfLRxLS3snv/0w2+w4/eqt3YW0d2ruuCze7CiiB24cF0XSYH/+suWYUy2pUVbXQlp+DfPHOcclJTB5JzjheIaFB3D/lUm8u7eY9HzXWbLgbO2dXbNerkwJJyk8wOw4ogfc3RQPXp1Mdmk9W60zf5zBh9YdGG8cP0CKg1JqsVKqCJgBbFRKfWJMLGGmh65OJirYh59/cNglB6c/PlRGRUMrd89IMDuK6IWFE6OJDfXlz585z9nDxoOljIoKYpgTvRnp62yl97TWsVprb611pNb6OqOCCfP4eXnw0/mjyCqtd8nB6Td25hMf5ufwa9uI8/N0d+OB2cPYV1jLjuOnzI5zSSW1p8koqGGBE501gFxWEhcwf1wUM5IG8cwnOVQ7+V2pZztcUkdafg13zUhw2B24xKUtnRJLRKA3f9ni+JtWfXSoa/D8RicabwApDuIClFL8cuEYGls7+L0L3Tn9+o58fD3duXWKLLLnzHw83fnOVcPYlVfNl0erzI5zURsOlDA6Koihg/3NjtIjUhzEBaVEBvLNyxN5a89J0lxgcLqmqY0P9pWwaFIMwX6eZscRfXTb9HhiQnz53cdHHHbs4UhZPXtP1rJoUrTZUXpMioO4qMfmpRAb6ssP1x6gpd15ly0A+MeuAlo7LNwzM9HsKMIA3h7ufH9eCgeL6/596cbRrNp1Ei8PN6c8U5XiIC7K39uD/71lPHlVTfxxs+Nf372QlvZOXtuRz5yREaREOsdNSOLSFk+KYXhEAM98mkNHp8XsOF/T1NrBe3uLWTAuyim3KZDiIC7piuGD+UZqHH//Is9pF+Z7J6OI6qY2vn1lktlRhIHc3RT/fd0I8iqbWJdZZHacr3l/XzGNrR3cfplzTpmW4iC65SfzRzE4wIv/Wbuf5rYOs+P0SKdF89IXeUyMC2Ha0DCz4wiDXTs6kolxIfxx81FOtznGpU+tNW/uOsmoqCAmxzv2jm8XIsVBdEuwrye/vWUcueUN3Pb33U41vfWTw2UUnGrmgdlJsqGPC1JK8eMbRlJa18ILnx83Ow4AmSdryS6t547L4p32mJPiILptzshI/nb7FLJK61n64g4Kq5vNjnRJWmv+7/PjDB3sz7zRzrEapui56UmDuHlCNC9+fpyCU+bvNb1qdwEB3h4snBhjdpRek+IgeuT6sUNYdd90qhpaWfLCDo5VNJgd6aJ2Hj/F/qI67ps1VG56c3E/nT8KTzfFrzdkmZqjsqGVDQdKWTQpmgBvD1Oz9IUUB9FjUxPDWPudy7FouPuVNCoaWsyOdF4Wi+Z/Pz7CkCAf2QZ0AIgM8uF7c4ezObuCLUfKTcvxt23H6LRovjVzqGkZjCDFQfRKSmQgr3wzleqmNu59Ld0hB6nX7y/hQFEd/3PdCHw83c2OI+zgnplDGRbuzxPrs0y5L6ek9jSrdp1kyeQYp1/xV4qD6LXxsSE8f9skDpfU8d1/7nWoeeYt7Z08/fERxsYEsXiS8173FT3j5eHGrxaO5WR1M89tzrV7/3/ZcgyN5ntzh9u9b6NJcRB9MndUJL9cOJbPjlTwsw8OOcwyBi9/eYKSuhZ+euNo3GSsYUCZmTyYFdPiWbk9j9159lu1teBUE++kF3LbtHhiQ/3s1m9/keIg+uzOyxJ4+Opk3tpTyK83ZJteIKoaW3lh23GuGRXJjGGDTM0izPH4/FHEh/nx2Jr9NLS026XPP24+ioe74qGrk+3SX3+T4iAM8YNrU7hnZiKvfHWCP3xq/9P5sz3zSQ4t7Z38+MaRpuYQ5vH39uDZZRMprTvNL//V/7OXskvreX9fMXfPSCQiyKff+7MHKQ7CEEopfr5gNMunxvH81mP8afNRU84gduWd4u20Qu6ZmehUu24J401JCOWhq5NZm1HEhgMl/dZPW4eFx9bsJ8zPiwdmD+u3fuzNeSfhCoejlOKpxeNo67Tw3OZcCmua+c3icXh52Oc9SEt7Jz959yBxYb48Nm+EXfoUju17c4fz5bEq/vud/cSH+TE+1vilLP64OZfs0nr+fleqUy6wdyFy5iAM5e6m+MOtE3hk7nDWZhRxx8v2W2rjr1uPkVfVxG8Wj8PXS6auiq4tRVfemcogf2/ufT2dktrThrafUVDNi58fZ1lqLPNGRxrattmkOAjDKaX4/rwU/rR8IvsKa1n8t68Mf1Ge60hZPS9sO84tk2OYNTy8X/sSziU80JtX75lKS1sn33otjcZWY+7JaWrt4Pur9xMd4svPFow2pE1H0qfioJT6vVLqiFLqgFLqPaWUcy4/KPrFwokxvPVfl1Hd2MZtf99FWV3/3El9oqqJh/+5lyBfTx6f73ovUtF3KZGBPH/7ZI5WNPJfr6dT38cZTG0dFh55ex+FNc08u2wigT6ut7NgX88cNgFjtdbjgVzgx32PJFzJlIRQXr93GlXWAlFRb2yB2HCghJv+8iVVja38ZcUkwlzomq8w1uyUcJ65dTxp+dUse3Fnr9+stHZ08p03M9icXc6vbh7jssvA96k4aK0/1VqfOUfbBcgCNsLG5PhQXrtnKmX1Laz4+y5yyvq+WF9bh4VffHCIh/+5l5TIADZ+bxYzkwcbkFa4ssWTYnnlm1MprG5m8d++6vGx2NLeyQP/yOCzIxU8uWgsd85I7J+gDsDIMYdvAR8Z2J5wIamJYbz6zalUN7Ux/89f8NuPsnu9HlNVYyt3vLSb13cWcN8VQ1n97RnEhPganFi4qitTwlnzwAw6LZqFf/2SZzfl0tSNcYi9J2tYvnIXW3Mq+c3icdzhpDu8dZe61Fx0pdRm4HwL4f9Ua/2B9TE/BVKBW/QFGlRK3Q/cDxAfHz+loKCgL7mFk6ppauN/PzrC6vRCYkJ8eWTucBZOisbbo3uziw6X1HH/GxlUNbby9NLxTr1evjBXad1pntqYzYYDpUQEevPINcO5ekQEUcE+/96gR2tNXlUTf9p8lPX7SwgP9ObnC0Zz04RoUzIrpTK01ql26auvNyoppe4GHgDmaq27tftLamqqTk9P71O/wrml5VfzxPrDHC6pJzzQm7tnJHD79IQLzhO3WDRv7i7gtx8eIdjXk5V3TemXOeti4MkoqOHJjVnsPVkLQESgN2Oig6hqbON4ZSPNbZ34eLpx/6wkvj17GP4m7tHgNMVBKXU98CwwW2td2d2fk+IgoOtd2Y7jp1i5PY/Pcyvx8nBjwbgobr8snsnxof9+93a8spEfrTtAWn4Ns4YP5g/LJhAR6BpLFAjHoLXmQFEde0/WsK+wliNlDYQHejMsPIDkiADmjoogKtj8S5fOVByOAd7AmaUPd2mtH7jUz0lxEOfKLW/gzV0FvJtZTGNrByF+nvh6uuPl4UZpXQs+Hm78bMFolk6Jddo9eYXoK6cpDr0lxUFcSFNrB+v3l3C4pI62DgttHRaCfT15aE6ynC2IAc+exUHWVhIOxd/bgxXT4s2OIcSAJ8tnCCGEsCHFQQghhA0pDkIIIWxIcRBCCGFDioMQQggbUhyEEELYkOIghBDChhQHIYQQNqQ4CCGEsCHFQQghhA0pDkIIIWxIcRBCCGFDioMQQggbUhyEEELYkOIghBDChhQHIYQQNqQ4CCGEsCHFQQghhA0pDkIIIWxIcRBCCGFDioMQQggbSmtt/06VagBy7N5xzwwGqswO0Q2S01iS0zjOkPEMZ8k6QmsdaI+OPOzRyXnkaK1TTeq7W5RS6Y6eESSn0SSncZwh4xnOklUplW6vvuSykhBCCBtSHIQQQtgwqzisNKnfnnCGjCA5jSY5jeMMGc9wlqx2y2nKgLQQQgjHJpeVhBBC2NJaX/IDiAO2AtnAYeAR6/fDgE3AUeufodbvD7I+vhF4/py2VgAHgQPAx8DgC/Q5xfq4Y8Cf+c9ZzpVAJtABLHXEjGf9/VJAA6mOmBN4Dthn/cgFak3O+RRQCDSe831vYLU1/24gsZ9yfsOa8TDw9EVeD2Yem33KaMdjs6/Ppb2Oze7m7PGx2cus84AM63OSAczp7v9pb4/PC/6bL/UAa6NRwGTr54HW/6zRwNPAj6zf/xHwO+vn/sAVwANn/2fQNXW2AusvB+vPP3GBPvcAMwAFfATcYP1+IjAeeIOvvwAdJuNZGbYDu/j6C9Chcp71mO8Cr5ic8zJrv+e+AB8EXrR+vhxY3Q85BwEngXDr168Dcx3s2OxzRjsdm4bktMOx2ZOcPT42e5l1EhBt/XwsUNyT56o3x+eFPrp1WUlrXaq1zrR+3kBXFYwBFlqf0DNP7CLrY5q01l8CLec0pawf/kopBQQBJef2p5SKAoK01jt117/qjbPaztdaHwAsjprR6td0HQBfa98Bc56xAnjLrJzWNnZprUvP81dn97kWmGtty8icSUCu1rrS+vVmYMm5QUw+Nvuc0aq/j02jcp7RX8dmt3Ja2+jxsdnLrHu11mdeH4cBH6WUd3efq94cnxfS4zEHpVQiXdVtNxB55gmz/hlxsZ/VWrcD36HrlKeErgr68nkeGgMUnfV1kfV7TpFRKTUJiNNab3DknGflSACGAltMzHkxMXSd0qO17gDq6HrXZ1hOuk7BRyqlEpVSHnS9oOIukMWUY9OIjPY4No3IeVaOfjs2e5DzYrp1bPYy6xJgr9a6le4fd306Ps/Wo+KglAoA1gGPaq3re9qZUsqTrl8Uk4Bouq71/fh8Dz3P97QzZFRKudF1vfQHjpzznK+XA2u11p0m5rxoM+f53tf+DX3NqbWuseZcDXwB5NN1bbbHWS7E7Iz2OjYNfi777djsQc6L6dbx0NOsSqkxwO+Ab/eknx487pK6XRysL/J1wCqt9bvWb5dbT2POnM5UXKKZiQBa6+PWU541wOVKKXel1D7rx6/oqnaxZ/1cLBe4FOGAGQPpula4TSmVT9e1yvVKqX/fmu8gOc+2nLNO203KeTFFWN/RWd/hBQPVBudEa/0vrfV0rfUMutb+Oupgx2ZfM9rr2DTyuezPY7O7OS/mosdmb7IqpWKB94C7tNbHz+rH5rky6vg8n24VB+s1tJeBbK31s2f91XrgbuvndwMfXKKpYmC0Uirc+vU8a5udWuuJ1o+fW0+zGpRSl1n7vutSbTtKRq11ndZ6sNY6UWudSNeg381a63RHynnW8zYCCAV2nt24vXNeoo2z+1wKbLEWGiNzopSKsP4ZStdA40sOdmz2KaMdj01Dnks7HJvdynmJJi54bPYmq1IqBNgI/Fhr/dWZB1/k/7TPx+cF6UuMWFv/nVfQdWpygP9ML7uRrmtrn9E1HeszIOysn8mnq4I20lXNRlu//wBdgzIHgH8Bgy7QZypwCDgOPM9/pmNNtbbXBJwCDjtaxnMes42vzwhxqJzAE8D/Osj/+dPWn7NY/3zC+n0f4B26rhHvAZL6KedbQJb1Y/lFXg9mHpt9ymjHY7PPObHPsdndnD0+NnuTFXjcevzsO+sjorv/p705Pi/0IXdICyGEsCF3SAshhLAhxUEIIYQNKQ5CCCFsSHEQQghhQ4qDEEIIG1IchBBC2JDiIIQQwoYUByGEEDb+P0TYRdebPyxTAAAAAElFTkSuQmCC\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAD8CAYAAACcjGjIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xd8XNd14PHfnUHvANE7C8BOgiRIsYnqXRbVi21JcZNL7MSJ14mdeB3XeOPsOtmNY8eyFVsu6l2UZBWLYgULSLCBFQDRQZRBbwPMzN0/MJApAiRRBvPemznfzwcfgijvHYL34cxt5yqtNUIIIcT5bEYHIIQQwnwkOQghhBhDkoMQQogxJDkIIYQYQ5KDEEKIMSQ5CCGEGEOSgxBCiDEkOQghhBhDkoMQQogxQoy4aXJyss7Pzzfi1kIIYVkHDhxo01qn+ONehiSH/Px8SktLjbi1EEJYllKqxl/3kmElIYQQY0hyEEIIMYYkByGEEGMYMucghJg+p8tNY+cgLd2DNPc4yYiPYHV+ktFhiQAhyUEIi2nuHuTJ3dX8YW8tXQPDH/ncpzbk881bFhIWIoMCYnokOQhhEf1DLr772nFeKqvH5dHctCid6xelkR4XQUpsOM/sr+XXu6o5WNPBTz++kpykKKNDFhamjDgJrri4WMtSViEmrrN/iE/9Zj+H6zp5eG0en944m7xZ0WO+7o/HzvH1Fw5jU4oXv7iOeamxBkQrZopS6oDWutgf95K+pxAm19Q1wH3/VUJ5Yzc/+8Qqvrt5ybiJAeDmJels+cpGQu02PvNkKZ39Q36OVgQKSQ5CmFhz9yD3/ryEpq5BnvzUGm5ekn7Z78mbFc0vHl5FU+cgf/nUQYbdHj9EKgKNJAchTOwHb5ygtdfJ059by7q5syb8favyEvnnu5eyq8LB97ccn8EIRaCSCWkhTGp3ZRuvH27kr68rYGl2/KS//95V2Zxu7uHx7VWsyE3grhXZMxClCFTScxDChIbdHv7p1XJykiL54tVzp3ydv795AavzE/n2q+U0dg74MEIR6CQ5XITT5WbrqRb+8eWj/MPLR3njSBMdfTK5J/zjyd3VnGnp5du3LyYi1D7l69htiv9zXxFuj+brLxzG4/H/6kRhTTKsdIFep4sfbDnOliNN9DpdRIfZsSnFU3trUQo2zkvm3x4oIjkm3OhQRYBq6R7k3987wzXzU7h+Yeq0r5c7K4r/efsivvnSUX5bUs1fbJg9/SBFwJPkcJ6mrgE+/ZtSTjf3cM/KLG5eks76ucmE2BSH67vYdrqVx7dXctfPdvHrv1jDvNQYo0MWAeg/t1bgdLn5p48tRinlk2s+uDqHd8rP8aO3TrKxIEXarrgsGVbyOlrfxeaf7qKuvZ8nHi3mx/cu59oFaUSE2gmx21iVl8jf3lDIs4+tY2DIzd0/28WeKofRYYsA09E3xHOl9dxZlEV+8vh7GaZCKcW/3LOMiFA7//jyUYzY/CqsZdrJQSmVo5TaqpQ6oZQqV0r9tS8C86fTzT3c/4sSQu02Xvzieq6ef/Gu/PKcBF7+0gZS4yJ45Il9HK7r9GOkItD9fk8NA8NuPrdpjs+vnRoXwd/dPJ+9Z9t5/UiTz68vAosveg4u4Gta64XAWuAvlVKLfHBdv/B4NP/w0lEiQm28/KX1zE+/fLmBnKQonv/8OpJjwvjrZ8roc7r8EKkIdIPDbp4sqebq+SkUps1M2YsHV+eyJCuOH75xnF5pt+ISpp0ctNZNWuuD3vd7gBNA1nSv6y/PH6ijtKaDb966kNS4iAl/X2J0GD95oIia9n6+97psMhLT90pZA229Qzw2A72GUXab4rt3LKG528l/vH9mxu4jrM+ncw5KqXxgBbDXl9edKY5eJz966yRrZidx36rJbxBaO2cWX7p6Ls+W1vHWUemmi6nzeDS/3FHFkqw41s2Z+E7oqViVl8h9q7J5YsdZKlp6Z/Rewrp8lhyUUjHAi8BXtdbd43z+MaVUqVKqtLW11Ve3nZZ/fvMkfU4XP7xzyZRXhXz1+kKWZcfzjZeOcq5r0McRimCx9VQLla19fO7KOT5boXQpf3/LAiLD7Pzzmydm/F7CmnySHJRSoYwkhj9orV8a72u01o9rrYu11sUpKSm+uO207DvbzosH63ls0xwKpjG+G2q38X8fXMHAkJv/J910MUW/3lVNVkIkty3N8Mv9kmPC+fymObx/soVjDV1+uaewFl+sVlLAE8AJrfVPph+Sf/z7e6dJiQ3ny9cUTPtas5Ojubc4mxdK62nqkhIFYnLqO/rZVdnG/cU5hNj9t7r84XX5xIaH8PMPKv12T2EdvmiJG4CHgWuVUoe8b7f64Loz5lBdJ7srHXx242wiw6ZemuB8X7xqLh6t+cW2Kp9cTwSPFw80AHDPKv+u44iPDOXhdXm8eayJylaZexAf5YvVSju11kprvUxrXeR9e9MXwc2Un39QQVxECJ9Ym+eza+YkRXHXiiye3ldLa4/TZ9cVgc3j0Tx/oI4Nc5PJTvT/sZ6f3jib8BCb9B7EGEG3Q7qipYe3y5t5dH0+MeG+rR7ypWvmMez28Ksd0nsQE7PnrIP6jgHuKzamnHZyTDgPrs7llbIG6jv6DYlBmFPQJYeff1BFRKiNv1if7/Nrz06O5mPLM/ndnhqp4Com5PnSemIjQrhp8eVPeJspj22ag1Lwy+3yokb8WVAlh4bOAV491MCDq3OZNUNVVb98zTz6h9z8aqc8aOLSugeHeetYE3csz5xWWe7pykyI5O4V2Tyzv07OnBYfCqrk8MSOswAzUrdmVEFaLLcvy+DXu6pl7kFc0pbDTQwOe7ivOMfoUPjUxnycLg/PldYZHYowiaBJDoPDbl44UMctSzPISoic0Xt97cb5OF0e/nNrxYzeR1jb8wfqKEiNYfkUjgD1tQXpcVwxO4nfltTglgOBBEGUHLYcaaJ70MXH1+TO+L1mJ0dzf3E2T+2tlUk+Ma669n7Kaju5a2WWX3ZET8Sj6/Op7xhg68kWo0MRJhA0yeGpvTXMSYlm7Zwkv9zvr64rAAX//p7smhZjvemtxXX70kyDI/mzGxalkR4XwZMl1UaHIkwgKJLDyXPdHKzt5ONrcv32Ki0jPpJH1ubx0sF6zjT3+OWewjreONrE8ux4cmf5f2/DxYTabXxybS47zrTJpjgRHMnh6b21hNlt3LPSv2vJv3TNPKLCQvjJu6f9el9hbrWOfo7Ud3HbMv/UUZqMB9fkEma38buSGqNDEQYL+OQwMOTmpbIGblmaTmJ0mF/vnRQdxqc25PPH8nNSGll86A3vkNKtfiqyNxnJMeHctiyDFw7Uy2FAQS7gk8OWI430+GkiejyPrs8nzG7j8e1SnkCMeONoI0U5CYaUy5iIR9bl0et08fLBeqNDEQYK+OTwzP465qZEs2a2fyaiL5QcE84Dq3N4uaxBznsQVLf1cayhm9tNOKQ0qigngaVZ8fy2pAatZVlrsAro5NDUNcCBmg7uXplt6HLBz105B4+GJ2TXdNAbHVK6xYRDSqOUUjy8Lo8zLb3sPdtudDjCIAGdHP547BwAtywxrm4NjFRsvX1ZBk/traWrf9jQWISx3jjSxMrchBnfiDldH1uWSXxkqExMm4i/NycGdHJ469g55qfFMiclxuhQ+PymufQNufndnmqjQxEGqW7r43hTtyknoi8UGWbn/uJs3i4/R3O3DIeaQUuPf/8fAjY5tPQMsr+6nZsN7jWMWpQZx9XzU/j1rmqcLrfR4QgDvHu8GcDQCqyT8cm1ebg8mqf21hodigCa/DxnGbDJ4Z3yZrQ213LBT22YjaNviHfKm40ORRjgnePnWJgRR06SOVcpXShvVjRXFabw9L5aht0eo8MJek2dkhx84o/HzjEnOZrCNOOHlEZdOS+ZrIRInt4nr8SCTVuvkwM1Hdy4KM3oUCblkXV5tPQ4P5y/E8bx9/n0AZkcOvqGKKlycPOSdNMUNQOw2RQPrclhd6WDs219Rocj/Oj9Ey149Ej9Iiu5en4q+bOi+NWOKlnWarBG6TlM37vHm3F7tKmGlEbdV5yD3aZ4Zr/0HoLJO8ebyUqIZHFmnNGhTIrdpvjMxtkcru9if3WH0eEEtXPd0nOYtjePNZGdaM4HMS0ugusWpPJCaT1DLhnHDQb9Qy52nGnlhkVppurJTtS9q3JIiArll3I2uqGk5zBNPYPD7Kpo4xaTDSmd76ErcnH0DX24ekUEth1n2nC6PJYbUhoVGWbn4bV5vHeiWYZDDSRzDtO0q6KNYbfm+oXmfRA3FaSQlRApQ0tB4p3yZuIiQgwr4eILD6/LI9Rmk13+Bhl2e2jx87HDAZcctp5sJTYihJV5iUaHclF2m+KB1TnsONNGrUNOigtkLreH9082c+2CVELt1n3cUmMjuHNFJi8cqKe9b8jocIJOc/cg/l4PYN3WOg6tNVtPtbCpMMX0D+I9q0bOlnj9SKPBkYiZdKCmg47+YW60yMa3S/nslXMYHPbw1F4pqeFvRhTtNPdv0Ekqb+ympcfJNfNTjQ7lsrISIinOS+S1Q5IcAtnWU62E2BRXFiQbHcq0FabFsn7uLJ4trcPj5zo/wa5RksP0fHBq5GD0qwpTDI5kYj62PJNTzT2cOifHiAaqbadbWZWXSGxEqNGh+MT9xTnUtQ9ItVY/a+r072Q0BFhy2HqqlWXZ8aTEhhsdyoTcujQDm4LXD0vvIRA1dw9yoqmbqy3Qk52omxanExsewvOldUaHElSaugaJDQ/x6z0DJjl09A1RVtthqQcxJTac9XOTee1wo+w+DUDbTrcC1unJTkRkmJ2PFWXy5rEmegal/Ly/NHUNkB4f4dd7Bkxy2H6mFY+GaxdYJzkA3LE8k9r2kQPnRWDZdqqV1NhwFmbEGh2KT91fnMPgsIctR5qMDiVoNHUNkuHnM0ACJjlsPdnCrOgwlmXFGx3KpNy0JJ1Qu+I1GVoKKC63hx1nWrmqMMW0mzGnanl2PAWpMTK05EeNnYNkSs9h8twezbbTIw+izWatBzE+MpSrClPZcqRRVoAEkEN1nXQPurhqfuAMKY1SSnF/cQ4HazupaJHFFDPN6XLT1uskI156DpN2uL6Tjv5hyz6IdxRl0tztZF+1rAAJFNtOt2JTcOU8a7bJy7lzRRZ2m+L5A/VGhxLwmrtGdkZnJEjPYdJ2nmlDKbiywJoP4vULU4kKs/NKWYPRoQgf2Xa6lRW5icRHBcYS1gulxIZzzfxUXilr8PvZxsFmtKZShgwrTd7OijYWZ8aRFB1mdChTEhUWws2L03njaBODw3KEqNW19To5Ut/F1QG0Smk8d67w9nhlz8OMGj0e1JLDSkqp/1ZKtSiljvniepPR53RRVtvBRot33+9amUXPoIs/nWgxOhQxTTvOeJewWnSYc6KuW5BGdJid1w5Lj3cmNXp7DpkWHVb6DXCzj641KfvOtjPs1mycZ+3yBOvnJpMWF87LZTKGa3U7TreRFB3GkkxrrZybrMgwOzcuTufNo+fkbJIZ1NQ5SHxkKFFhFtwEp7XeDhjSt9xZ0UZ4iI3ifPNWYZ0Iu02xuSiLD061StVLC9Nas6uyjXVzZ1lu5dxU3FGUSdfAMNu9G/6E7zV1Dfp9vgECYM5h55k2VucnERFqNzqUabtrRRYuj2aLVGq1rMrWXpq7nZbvyU7UxnnJJEaF8qrs05kxTV0DgZ0clFKPKaVKlVKlra2+eZXR0jPIqeYeNgTIg7gwI44F6bG8dFDGcK1qV4UDgA1zA6NNXk6o3cZtyzJ49/g5+pwuo8MJSPUdA2Ql+ncyGvyYHLTWj2uti7XWxSkpvpmo21XRBhAQ5ZBH3b0yi0N1nVS19hodipiCnRVtZCdGkjsryuhQ/GZzURaDwx459nYGdA8O0zUwTE6i/9uTpYeVdp5xkBgVyqKMOKND8ZnNRVnYFLwsex4sx+X2sKfKETRDSqNW5SaSGR8hJWBmQEPHyEqlbKsmB6XU00AJMF8pVa+U+owvrnspWmt2VrSyfl5yQE38pcVFsKkwhWf31zHslhUgVnKssZueQRfrgyw52GyKjxVlsv20LKbwtbr2kWOEc5IsOqyktX5Ia52htQ7VWmdrrZ/wxXUvJZAn/h5em0dLj1O66RYzOsy5fu4sgyPxv83LRxZTvCGLKXyq3uo9ByPsPDPyIAZicrh6fipZCZH8fo+c1WsluyraWJAeS3KMNQ6b8qWFGbEUpsXwqhx761N1Hf1EhdlJNKAMi3WTQ0UbebOiyEkKvIk/u03x8Sty2V3pkKqXFjE47Ka0piNgVs5NllIj+3RKazo+HAoR01ffMUBOYpQhZd8tmRyG3R72VLUH9IP4wOocQu2K3++pNToUMQGl1R0MuTwB2ZOdqM1FmQAyMe1D9R0DZBuwjBUsmhyO1HfS63RxZQA/iMkx4dy6NIMXD9TTPyTrx81uV2UbITbFmtlJRodimOzEKFbnJ/JKWYMce+sDWmvq2/slOUzGDm+J7nUBPvH38No8epwuGce1gN0VbRTlJBDt50PgzWZzURZnWno53tRtdCiW1z3gosfpMmzo3JLJYVdFG0uz4kmIsmaJ7olalZfIgvRYmZg2ue7BYY42dAX8i5WJuHVpBiE2JS9ofKCuY2TuRnoOE9TrdFFW2xkUY7tKKR5cnUN5Yzcn5JWYaZVWt+PRsG6OJIek6DCuKkzhtUONcgjQNNV/mByk5zAhe6scuDzWL9E9UXcUZRFqV7woxzGaVkmlgzC7jZV51q4M7Ct3rsjiXPcgJZUOo0OxtNE9DkaUzgALJofREt3B8iAmRYeNHMd4qBGX7Jg2pZIqBytyEwKiMrAv3LAojfjIUJ4trTM6FEura+8nNjyEuEhj5rEslxx2VbSxZnZglOieqHtXZdPW62T7GamZbzZd/cOUN3bLfMN5IkLt3LUii7ePnaNDymlMWX3HANlJxuxxAIslh+buQU439wbNkNKoq+enkhQdxosHpBif2ew960DLfMMYD6zOYcjt4ZVD0manqq7DuGWsYLHkMFq7ZmMAleieiLAQG3csz+Td48109Q8bHY44T0mVg/AQG0W5CUaHYioLM+JYlh3Ps/vrZM/DFGitDd0ABxZLDjsrRs7mXZgeOCW6J+reVdkMuT28LoXNTGVPVTvF+YmEhwTPMOdEPbA6h5PnejhS32V0KJbT3jdE/5DbsMlosFBy0Fqzq6KN9UFyNu+FFmeOnBL3gqxaMo2OviFONHXLkNJF3LE8k8hQO8/sl4npyfpzNVbpOVxWZWsfzd3OgK6ndClKKe5ZmS2nxJnI3rMjSzXXSnIYV2xEKLcuzeD1w41SAmaSPlzGamBhUcskh9H5hmA5m3c8dxRlohS8IrtPTaGk0kFkqJ1l2TLfcDEPrsmh1+nijSNNRodiKaO7o404O3qUpZJDsJ3Ne6G0uAg2zE3m1UNS2MwMSqocFOcnEhZimcfI74rzEpmdHC3DoZNU39FPfGQocRH+P8dhlCVatdujKQnCs3nHs7kokxpHP2V1nUaHEtTaep2cbu6V/Q2XMTIcmsXes+1yzsMk1LUPGHI06PkskRyONnQF5dm847l5STrhITZeLZP140baUzUy3yCT0Zd318pslIIXD0rvYaLqO/rJTjB2lMQSySGYz+a9UGxEKNcvSuP1I00MSzkNw+ypchAdZmdpVrzRoZheVkIk6+fO4sWD9XikGN9lmWGPA1gkOeyuDN6zecdzV1EW7X1D7JByGoYpqXSwenYSIXZLPEKGu2dlNnXtA+yvbjc6FNNr6XHidHkMPwLZ9C17cNjN/urgPZt3PJsKU0iMCuWVMlm1ZISW7kEqW/tkSGkSbl6STnSYXYaWJqC6rQ+A/ORoQ+MwfXI4UDNyNu+GefIgjgoLsXHbsgzeOX6OXqesH/e3ktH5BhnmnLCosBBuXZrBm0fPyZ6Hy6h2jCSH2bMkOVzSrorRs3nlQTzfHcuzGBz2sO2UDC35254qB7HhISzOlPmGybh3VTa9Thdvl58zOhRTO9vWT6hdkZkQYWgc5k8OlQ6KchKICfKzeS+0MjeB2IgQtp+W5OBvJZUO1sxOwh6EZVymY3V+ElkJkbxxRJLDpVS39ZGTFGX4fJapk0PP4DBH6zul+z6OELuNjfOS2X6mVTbE+VFT1wDVjn5pk1Ngsyk2FSaPnOYoK+0uqtrRZ/iQEpg8OeyXs3kvaVNhCk1dg1S0SK0lfxnd3yD1lKZmw7xkepwujjRIpdbxaK2pdvSRJ8nh0vZUtcvZvJewqTAFgG0ytOQ3JZUO4iNDWZQRfGXjfWG9tzbarjNtBkdiTs3dTgaHPcxONr5MkMmTg4MiOZv3orISIpmbEs12edD8pqRqZL4hGMvG+0JSdBiLM+PYWSFtdjxnTbKMFUycHLoHhznW0CVDSpdxZUEKe6scDA67jQ4l4DV0DlDXPiBtcpo2zkvmYG2HLGkdx+gy1nwZVrq4/WdH5htkbPfSripMwenysO+s7DydaaXe3b1rZicZHIm1bZiXzLBbs7+6w+hQTKfa0UeY3UZmgrGlM8DEyaGk0kFYiI0VcjbvJV0xJ4kwu02WtPrB/up2osPsLEiPNToUS1udP9Jmd8nQ0hgjy1gjTbFM2rTJYc9ZBytlvuGyosJCWD07kR0y7zDjSqs7WJmXaPj6c6uLDLOzKi+RndJmx6hu62e2CeYbwKTJoWtgmPLGbhlSmqBNBSmcau7hXNeg0aEErK6BYU4197BKVs75xMaCZI43dePodRodiml4PCPLWM0w3wA+Sg5KqZuVUqeUUhVKqW9M93r7zrajZb5hwkaXtMrQ0swpq+1A65EhETF9o+X3d1c6DI7EPM51D+J0eUyxUgl8kByUUnbgP4FbgEXAQ0qpRdO55p4qB+EhNopyZL5hIhakx5KdGMkbR+Wc3plyoKYDu01Jm/SRpVnxxEaEyLzDeT4suBcoyQFYA1Rorau01kPAM8Dm6VxwT5WDlbmJMt8wQUop7lieyc6KNtqkmz4j9le3sygjjmip8eUTIXYb6+bMYseZNin/4lXdNnKMqll6Dr5o6VlA3Xl/rweuuOR3nDoFV1897qdcHs3/rG4nOzEKnjJ+OZdVfHnIzab6Tlx/jIY4Y6s5BhqPhr+tbic1LhxeNMeDGwi+1z1IdVsfg68nECkvBNnQ3s+zXYNk7jHH0KUveg7jrbka81JAKfWYUqpUKVU6PDx80Yt1D458Li4y1AehBY+oMDuRYSE4eoeMDiXg9A+58GhNbIS0SV+K9z7jXQMX/30QTAaH3USE2sb9hWoEX/Qc6oGc8/6eDYw5okxr/TjwOEBxcbHmgw/Gvdh/bjnO7/fUcOQ7N0KIvJqYjJ1bK/jXt0+x4++uMfyIwUDyzI4qfvDGCfb+w3XSK/OhCK352r9sZXFmHI8/Umx0OIb7yr9tI29WNL+81M9C+S91+KLnsB8oUErNVkqFAQ8Cr031YnuqHKzKSyRcEsOk3bE8E4DXj8jxob50oKaDnKRI0iQx+JRSiisLkimplBLeHo+mxmGePQ7gg+SgtXYBXwbeBk4Az2mty6dyra7+YY43yf6GqcpJimJlbgKvHZLk4Ctaj5R5KM4zxzhwoNlYMFLC+3B9cJfwbhpdxmqSPQ7go30OWus3tdaFWuu5WusfTvU6+6plf8N0bS7K4uS5Hk6d6zE6lIBQ295PW6+T4nzZ/DYTNsxNRimCfrf02dbRgnvmGQ421Q7p0f0Ny3PkbN6punVpBnab4tVDDUaHEhBKvcXhpOcwMxKjw1iaFc/OiuDewFnRMvJibl5ajMGR/JnpkoPMN0xPSmw4G+Yl89rhRjweWT8+XWV1HcSGh1CQap6HNtBsnJdMWW0nvc7gLeF9uqWX+MhQUmLCjQ7lQ6ZJDjLf4Dt3rcikvmOAA7VSEnm6ymo7WZ6TIIf7zKCNBcm4PJo9QVxKo6K5l8K0GJQfVyNdjmmSg8w3+M6Ni9KJDLXzcpkMLU1H/5CLk+d6pGz8DFuVl0hkqD1oT4fTWnO6pYd5qeYqBW+a5CDzDb4THR7CTYvTeONIE06XnBA3VUfru3B7tCSHGRYeYmfN7CR2nAnOeQdH3xCd/cOmG7o0VXKQ+QbfuXNFFl0Dw3xwKjgfOF84VNcJQFGOrFSaaVcWJFPZ2kdD54DRofjd6eaRyejCNOk5jCHzDb63cV4yyTFhvCJDS1NWVttJ/qwokqLDjA4l4F3lLTu/IwjLzle09AJQYKKVSmCS5DA633CFnM3rMyF2Gx9bnsmfTrRI7Zop0FpzsLZDSnT7ybzUGDLiI9gWhMnhTHMvsREhpMaaZ6USmCQ5lFR6z2+QsV2fumtFFkNuD2/JOQ+T1tQ1SEuPkxW5MqTkD0opNhWksLOiLehKaZxu7qEg1VwrlcAsyUHmG2bE0qx45qREy6qlKSirHZlvkMlo/9lUmELPoIvD9Z1Gh+JXFS29pptvABMkh46+IU40dbNO5ht8TinFnUVZ7KtupzEIJ/qmo6y2g/AQGwvS44wOJWhsnJeMTcG2IFpE4eh14ugbYp7JViqBCZLD3rMjG1/WzZXkMBM2F2WiNWyRSq2TUlbXydKseMJCDH9EgkZ8VChFOQlsC6I6S2c+nIyWnsMYe6raiQy1syxbuu8zIW9WNMtzEnhVKrVO2JDLw9GGLhlSMsCmwhSO1HfS0Rcch1aNJodCk61UAhMkh5JKB8X5ifIKbQZtXp5JeWP3h8W9xKWdaOpmyOWRyWgDbCpMQWuCZrd0RXMPMeEhpJvwrBBDfyM7ep2cau6R/Q0z7PZlGdgUcs7DBJV5a1JJz8H/lmcnEB8ZGjRLWk839zLPhCuVwODksKeqHZD5hpmWGhfB+rnJvHq4Ea2lUuvl7K/pID0ugoz4SKNDCTp2m2JjQTI7zrQGRVs909JrurIZowxNDiVVbUSH2VmaJfWUZtodRZnUOPqD/sSty9F6pDqovGAxzlUFKTR3OzkZ4AdWdfQN0dbrNOUyVjA6OVQ6WD07iVC7zDfMtJuXpBMWYpNDgC7jdHMvjr4hSQ4Gumr+SCmNQK8LVtE6MhltpgN+zmfYb+WWnkEqW/tkvsFP4iJCuXZ+Kq8fbgoy5MueAAAXt0lEQVS6HaiTsbtyZCJU9t0YJy0ugkUZcWw91WJ0KDNqtOCeDCtd4MP5BnkI/ebOFVm09TrZHqSlkSeipNJBTlIkOUnmOcs3GF27IJUDNR0BXRfsRFM3seEhZJp0bsuw5FBS2UZseAiLM2UHqr9ctzCV5Jgwnt5XZ3QopuT2aPaebZcXLCZwzYIU3B4d0Gc8lDd2sygzzrSnDBqWHHZXOrhizixCZL7Bb0LtNu5Zlc37J1to6R40OhzTOdHUTdfAsMw3mEBRTiIJUaFsPRmYycHt0Zxo6mZxpnkX4xjym3nI7aHG0c96eQj97sHVubg9mucP1BsdiumUeM8wXjcn2eBIhN02UqV12+kWPJ7AW9Ja1drL4LDH1CMnhiSHPqcLgA3z5CH0t9nJ0aydk8Sz++sC8qGbjpIqB3OSo0mPN99u1WB07YJU2nqHONYYeMuvR/9NS0y8jN+Q5NA76CI5JsyU9USCwUNrcqlt76ekymF0KKbhcnvYd7ZdhpRMZFNhCkrB+ycDb9VSeUM34SE25qZEGx3KRRmTHJwu1s1NNuWW8WBw0+J04iNDeXpfrdGhmMbRhi5vu5TkYBZJ0WEU5SSwNQD3O5Q3drMgI87Uc66GRObyaDbIQ2iYiFA7d63I4p3yZtqDpPrl5ez2zjfIvhtzuWZ+KkfqO2nrdRodis9orSlv7DL1fAMYuFpp/VyZbzDSx6/IZcjt4Zn90nsA2FPloDAthuQYc53jG+yumZ+K1oF1AFB9xwDdgy5JDuMJtdvInSWbjIxUmBbLlQXJPLm7miFXcO+Ydrrc7K9ulxcsJrQ4M460uHD+dLLZ6FB85liDdzLaxMtYwaDkEBMeYsRtxQU+s3E2zd1OXj8c3KW8y2o7GRz2yOo5E7LZFNctTGPbqVacLrfR4fhEeWM3dptifro5C+6NkuQQxK4qTKEwLYZf7qgKivLIF7O7og2bgjWzk4wORYzjhoVp9A25P9yHYnXljV0UpMYQEWo3OpRLMiY5REhyMAOlFJ/dOIeT53rYVREYD95U7Kp0sNR7yIwwn3VzZxEVZufd44ExtHTMWzbD7AxJDiEmrSUSjDavyCQ5Jpxf7qgyOhRD9DpdHK7rlNVzJhYRamdTQQrvnWi2/MbNlu5BWnucpi6bMcq8i2yFX4SH2Hl0XR7bTrd+WEI4mOw76xhZWi3zDaZ2w6I0mrudlt8tXd7YDcAS6TkIK/jk2jwiQm38etdZo0Pxu10VDsJCbKzKSzQ6FHEJ1yxIxaaw/NBSuTe5BfywklLqPqVUuVLKo5Qq9lVQwr8So8PYvDyLV8oaA7p+/nh2VzpYlZto+snBYJcUHUZxfpLlk8PRhi7yZkURG2H++a3p9hyOAXcD230QizDQw+vyGBh280IQVWt19Do50dTNhnky32AFNy5K4+S5Hura+40OZUq01hyo6WRlrjV6qdNKDlrrE1rrU74KRhhnSVY8K3MT+P2eGstP+k3UaOHB9TLfYAnXL0wDrDu0VNc+QFuv0zJDmDLnID70yLp8zrb1sbOizehQ/GJXhYOY8BCWmbhssviz/ORoClJjeLv8nNGhTElpzcjRyAGTHJRS7ymljo3ztnkyN1JKPaaUKlVKlba2Bk6dlEByy9J0ZkWH8duSGqND8YvdlW1cMTvJ1JUxxUfdsiSd/dXttPZYrxBfaU0HseEhFKaZe2f0qMs+FVrr67XWS8Z5e3UyN9JaP661LtZaF6ekpEw9YjFjwkPsPLgmh/dPNlPfYc1x3YmqcfRR4+hnY4EMKVnJrcsy8Ggs2Xs4UN3BirxE7BbZ5yUvmcRHfPyKPAD+sDewq7Vu9R4gc838VIMjEZMxPy2WOcnRvHWsyehQJqVrYJjTLT0UW2RICaa/lPUupVQ9sA54Qyn1tm/CEkbJSojkuoVpPLe/LqCrtW491cqc5Gjyk817EpcYSynFrUszKKl04LDQGQ9ltR1obZ35Bpj+aqWXtdbZWutwrXWa1vomXwUmjPOJK3Jx9A1Zsus+EQNDbkqqHFwtvQZLumVpOh4N71ho1dKBmg7sNkVRToLRoUyYDCuJMTYVpJCdGMkf9gbmxPTuyjaGXB6uWSBzX1a0KCOO/FlRvHnUOkNLpdUdLMyIJdpCFaklOYgxbDbFx6/IZU9VOxUtvUaH43NbT7UQFWaXEt0WpZTilqUZ7K500GGBY25dbg+H6jopzrNWe5PkIMZ136ocQmyKp/cF1sS01pqtJ1vZMC+Z8BApmWFVty3NwO3RvHPc/EOfJ5p6GBh2s9JC8w0gyUFcREpsODctSeeFA/UMDgfGCVwAFS29NHQOyColi1ucGUdOUiRvHDV/chjd/GallUogyUFcwifW5NI1MGypsd3LeX90CavMN1iaUorblmayq6KNNpOvWiqt6SAzPoLMhEijQ5kUSQ7iotbNncWc5Gh+vydwJqa3nmphQXosGfHWelDFWHetyMLt0aY+A11rTWl1O6vyrTXfAJIcxCUopfjE2jwO1nZyuK7T6HCmrXtwmNLqDq5ZIENKgWB+eiyLMuJ4pazB6FAu6nRzL83dTjZasPKvJAdxSfcXZxMbHsKvdlr/IKBtp1pxeTTXSnIIGHevzOJwfReVreZcVbft9Mgw5qZC6w1jSnIQlxQbEcqDa3J482gTDZ0DRoczLW+XnyM5Jswy9fTF5d2xPBObgpcPmrP3sO10K/PTrDmMKclBXNaj6/MBeHJ3taFxTIfT5eaDU63csCjNMoXPxOWlxkWwYV4yL5c1mO4ckv4hF/vPdrCp0JrFHSU5iMvKToziliXpPL23ll6ny+hwpmR3pYNep4sbF6UbHYrwsbtXZtHQOUBpTYfRoXzEnioHQ24PVxVacxhTkoOYkM9eOYcep4vn9tcZHcqUvFN+jugwO+stODEoLu3GRelEhtp5ucxcR9xuO9VKZKid4nxrDmNKchATUpSTQHFeIv+96ywut7Wqtbo9mnePN3P1glTZFR2AosNDuHlJOluONJlqw+b2M22smzuLiFBrtjlJDmLCvnj1XOo7BvjJu6eNDmVSymo7aOsd4qbFMqQUqO5dlU3PoMs05zzUOPo429bHJgsfJiXJQUzYdQvTeGhNDj/7oJL3T1qnXPLb5ecItSuumW+95YRiYtbNmcXs5Gh+v8cctcC2nx45CvkqC5dpkeQgJuWfPraYRRlx/M2zhy1xlKjWmneON7N+bjKxEaFGhyNmiM2m+MQVuRyo6eBEU7fR4bDtdCu5SVHkz4oyOpQpk+QgJiUi1M7PPrESj0fz5afKTH9a3KnmHmoc/TKkFATuWZlNWIiNpww+4tbpcrO70sFVhSkoZd1l05IcxKTlJ0fz43uXcaiuk59/UGl0OJe05XATNgXXL7Ju915MTGJ0GLcvzeDlsgb6DFxyvf10G/1Dbq5daO02J8lBTMktSzO4bVkGP/uggrp2cw4vuT2aFw7Us6kwhdTYCKPDEX7wibV59DpdvHrIuGJ8rxxqICk6jI3zrDsZDZIcxDR867aF2JTie1uOGx3KuHZWtHGue5D7i3OMDkX4ycrcBBakx/KHvTVo7f8d071OF+8db+b2ZRmE2q3969Xa0QtDZcRH8pXr5vHu8Wa2nmoxOpwxniutIzEqlOss3r0XE6eU4pNr8yhv7OZgrf8rCb997BxOl4fNRZl+v7evSXIQ0/LZjXOYkxzNd18rx+kyzwakjr4h3i1vZnNRlmx8CzJ3rcgiPjKUX2zz/3zYK4cayEmKDIjijpIcxLSEhdj4zh2LqXb08/i2KqPD+dCrhxoYcntkSCkIRYeH8Oj6fN453szp5h6/3be1x8muijY2L8+y9CqlUZIcxLRtKkzhtqUZ/Mf7FZzx48N4Kc8fqGdJVhyLMuOMDkUY4FPr84kKs/NfflxNt+VIIx5NQAwpgSQH4SPfuWMxUeF2/u7FI7gNLp18rKGL8sZu7lslvYZglRgdxkNrcnn1cKPfVtO9cqiRRRlxFKTF+uV+M02Sg/CJlNhw/uljiyir7eQ3Bp/78FxpHWF2W8C8ghNT87kr52BT8Pj2mR/uPNvWx+G6Tu5cEThtTpKD8Jk7i7K4Zn4K//vtU9Q6jNn70D04zIsH6rltWQYJUWGGxCDMIT0+gntWZvNsaR0tPYMzeq/f7DpLiE2xuShrRu/jT5IchM8opfjhXUux2xRfe/4QwwaU9n5ufx19Q24+vWG23+8tzOcLV83F5fbM6E7+tl4nz+yv464VWaTFBc5mS0kOwqcyEyL54V1L2F/dwff9vDnO7dH8Znc1q/MTWZod79d7C3PKT47mgdW5/K6khoqWmVks8etdZxlye/jC1XNn5PpGkeQgfG5zURaPbZrDb0tqeHqf/4qgvXeimfqOAek1iI/4HzcWEhlm57uvH/f5rumewWF+W1LDzYvTmZsS49NrG02Sg5gRf3/zAq4sSObbrx6jtLrdL/f8751nyUqI5IZFaX65n7CGWTHh/M31hew408Z7J3y7k//3e2rpGXTxpavn+fS6ZiDJQcwIu03x04dWkpUQyRd+f5D2vqEZvV95Yxd7z7bz6Po8Qixe00b43sPr8piXGsMP3jjus538g8Nunth5lisLkgNyGFOeIjFj4qNC+fknV9E1MMS3Xz02o/f65fYqosLsPFCcO6P3EdYUarfx7dsXUePo55c+Wtr61N5a2nqdAdlrAEkOYoYtzIjjq9cXsuVIE28cmZnzfX+z6yyvHGrkkXX5xEfJaW9ifJsKU7hlSTr//t4ZDtRMb6jzbFsf//r2Ka4sSGbtnCQfRWgu00oOSql/VUqdVEodUUq9rJRK8FVgInB8ftMclmfH861XjtLa4/TptV873Mh3txznpsVpfP2m+T69tgg8/+ueZWQmRPKlPxykrXdqbdHl9vC3zx0i1K7413uXB0QdpfFMt+fwLrBEa70MOA18c/ohiUATYrfxv+9bTt+Qm2+9ctRnK0Z2nGnla88dYnV+Ev/3wRXYbYH5kArfiY8M5eefXEln/zBfeaoM1xT24vxiexVltZ18/84lpMcHzr6GC00rOWit39Faj57HtwfInn5IIhAVpMXytRsKebu8me++fhzPNOsvnWjq5gu/O8DclBh++UgxEaFSlltMzOLMeH5w5xJKqhz8+O1Tk3qxcqyhi3979zS3LcvgjuWBUypjPCE+vNangWd9eD0RYD535Rxaepw8sfMsrb1OfnL/8imdtdDW6+SzT5YSExHCbz61hvhImWcQk3NfcQ6H6jp5fHsVXf3DfO/OxZdti4fqOnnst6UkRYfxg81LAnY4adRlk4NS6j0gfZxP/aPW+lXv1/wj4AL+cInrPAY8BpCbKytKgpHNpvjWbQtJjQ3nR2+dpKNviJ9/ctWkfrk7XW4+/7sDOPqcPPf5dQHdrRcz6/ubl5AUHcZ/vF9BZWsv//XwKpJjwsf92pfL6vn7F4+SGhvOrx4tJjE68Ot2qemO/yqlHgW+AFyntZ5QtbXi4mJdWlo6rfsKa3vpYD1/98IRZsWE8aO7l3LtgstvXNNa87XnD/PSwQZ++vEV3L4ssLv1wj9eP9zI1184THxkKPeszOamxeksy46ne9DFwdoO3ilv5ul9taydk8TPPrGKJAMTg1LqgNa62C/3mk5yUErdDPwEuEpr3TrR75PkIACO1Hfy9eePcKq5h7tXZvHt2xddtJKqx6P53pbj/GZ3NV+9voCvXl/o52hFIDvW0MU/v3mCvWfbcXs0CVGhdA0Mo/XIhs5PXpHLt25fRKjBGyytlBwqgHDA4f3QHq31Fy73fZIcxCiny81P36/gZx9UEh8ZyjduXsC9q7KxnbfyyO3RfOPFIzx/oJ7PbJzNt25bGPDjvcIYHX1DvH+yhd2VDvJmRVGcn0hRTgJRYb6cnp06yySHqZLkIC50vLF7pA5TTQdFOQl8ftMc4iNDiQyz88TOs2w50sRfXVfA31xfIIlBBC1JDiIoaa156WADP3rrBG29H63F9A+3LuCxTYFVElmIyfJncjBHX0kIRg4LumdVNjcvSaeipZeBYTcDQ24So8MoypHN90L4kyQHYTrR4SEsl2QghKGk8J4QQogxJDkIIYQYQ5KDEEKIMSQ5CCGEGEOSgxBCiDEkOQghhBhDkoMQQogxJDkIIYQYQ5KDEEKIMSQ5CCGEGEOSgxBCiDEkOQghhBhDkoMQQogxJDkIIYQYw5DDfpRSPcApv994cpKBNqODmACJ07ckTt+xQoyjrBLrfK11rD9uZNR5Dqf8dZrRVCmlSs0eI0icviZx+o4VYhxllViVUn47QlOGlYQQQowhyUEIIcQYRiWHxw2672RYIUaQOH1N4vQdK8Q4yiqx+i1OQyakhRBCmJsMKwkhhBhLa33ZNyAH2AqcAMqBv/Z+PAl4Fzjj/TPR+/FZ3q/vBX56wbUeAo4CR4A/AskXuecq79dVAP+PP/dyNgEHARdwrxljPO/z9wIaKDZjnMC/AYe8b6eBToPj/CFQB/Re8PFw4Flv/HuB/BmK8wFvjOXAjy/xPBjZNqcVox/b5nR/lv5qmxONc9Jtc4qx3gAc8P5MDgDXTvT/dKrt86L/5st9gfeiGcBK7/ux3v+sRcCPgW94P/4N4F+870cDG4EvnP+fwcjS2Ra8vxy83/+di9xzH7AOUMBbwC3ej+cDy4Df8tEH0DQxnhfDdmAPH30ATRXneV/zFeC/DY5zrfe+Fz6AXwL+y/v+g8CzMxDnLKAWSPH+/UngOpO1zWnH6Ke26ZM4/dA2JxPnpNvmFGNdAWR6318CNEzmZzWV9nmxtwkNK2mtm7TWB73v9zCSBbOAzd4f6OgP9k7v1/RprXcCgxdcSnnfopVSCogDGi+8n1IqA4jTWpfokX/Vb8+7drXW+gjgMWuMXt9npAF85PomjHPUQ8DTRsXpvcYerXXTOJ86/54vANd5r+XLOOcAp7XWrd6/vwfcc2EgBrfNacfoNdNt01dxjpqptjmhOL3XmHTbnGKsZVrr0eejHIhQSoVP9Gc1lfZ5MZOec1BK5TOS3fYCaaM/MO+fqZf6Xq31MPBFRro8jYxk0CfG+dIsoP68v9d7P2aJGJVSK4AcrfUWM8d5Xhx5wGzgfQPjvJQsRrr0aK1dQBcjr/p8FicjXfAFSql8pVQIIw9UzkViMaRt+iJGf7RNX8R5Xhwz1jYnEeelTKhtTjHWe4AyrbWTibe7abXP800qOSilYoAXga9qrbsnezOlVCgjvyhWAJmMjPV9c7wvHedj2goxKqVsjIyXfs3McV7w9weBF7TWbgPjvORlxvnYR/4N041Ta93hjfNZYAdQzcjY7KRjuRijY/RX2/Txz3LG2uYk4ryUCbWHycaqlFoM/Avw+cncZxJfd1kTTg7eh/xF4A9a65e8H272dmNGuzMtl7lMEYDWutLb5XkOWK+UsiulDnnfvsdItss+7/uyuchQhAljjGVkrPADpVQ1I2OVrymlPtyab5I4z/cg53XbDYrzUurxvqLzvsKLB9p9HCda69e11ldordcxUvvrjMna5nRj9Ffb9OXPcibb5kTjvJRLts2pxKqUygZeBh7RWleed58xPytftc/xTCg5eMfQngBOaK1/ct6nXgMe9b7/KPDqZS7VACxSSqV4/36D95purXWR9+3b3m5Wj1Jqrffej1zu2maJUWvdpbVO1lrna63zGZn0u0NrXWqmOM/7uc0HEoGS8y/u7zgvc43z73kv8L430fgyTpRSqd4/ExmZaPyVydrmtGL0Y9v0yc/SD21zQnFe5hIXbZtTiVUplQC8AXxTa71r9Isv8X867fZ5UfoyM9bef+dGRromR/jz8rJbGRlb+xMjy7H+BCSd9z3VjGTQXkay2SLvx7/AyKTMEeB1YNZF7lkMHAMqgZ/y5+VYq73X6wMcQLnZYrzgaz7goytCTBUn8B3gf5nk//zH3u/zeP/8jvfjEcDzjIwR7wPmzFCcTwPHvW8PXuJ5MLJtTitGP7bNaceJf9rmROOcdNucSqzAt7zt59B5b6kT/T+dSvu82JvskBZCCDGG7JAWQggxhiQHIYQQY0hyEEIIMYYkByGEEGNIchBCCDGGJAchhBBjSHIQQggxhiQHIYQQY/x/2fIKTqh/4VMAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] @@ -738,36 +784,83 @@ } ], "source": [ - "from scipy import signal\n", - "\n", - "raw_data[\"smooth\"]=signal.savgol_filter(raw_data['Concentration'], 57, 1)\n", - "\n", "raw_data[\"smooth_b\"]=signal.savgol_filter(raw_data[\"Concentration\"]-raw_data[\"smooth\"], 57, 3)\n", "\n", "fig, ax = plt.subplots()\n", - "ax.plot(raw_data[\"convert\"],raw_data[\"Concentration\"],linewidth=.75)\n", - "ax.plot(raw_data[\"convert\"],raw_data[\"smooth\"],linewidth=2)\n", - "#ax.set_xlim('2015-07-02','2020-02-01')\n", - "\n", - "fig, ax = plt.subplots()\n", - "ax.plot(raw_data[\"convert\"],raw_data[\"Concentration\"]-raw_data[\"smooth\"])\n", + "ax.plot(raw_data[\"convert\"],raw_data[\"smooth_b\"])\n", + "plt.axhline(y=0, color='r', linestyle='-')\n", "ax.set_xlim(\"2018-01-01\",'2020-01-01')\n", + "ax.set_ylim(-2.5,2.5)\n", + "\n", + "plt.show()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "On peut estimer la période plus facilement." + ] + }, + { + "cell_type": "code", + "execution_count": 42, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[22, 23, 26, 26, 26, 26, 26, 25, 46, 1, 26, 1, 1, 1, 1, 1, 1, 1, 1, 1, 47, 26, 24, 24, 24, 52, 27, 26, 26, 26, 26, 26, 51, 28, 27, 51, 26, 26, 51, 52, 52, 27, 25, 27, 26, 51, 53, 27, 26, 24, 24, 25, 26, 26, 26, 26, 25, 27, 26, 26, 26, 26, 25, 27, 79, 26, 25, 27, 25, 27, 26, 26, 26, 27, 24, 28, 26, 25, 26, 27, 25, 27, 25, 27, 25, 27, 25, 25, 25, 50, 24, 25, 51, 52, 52, 104, 27, 25, 51, 26, 26, 26, 25, 27, 78, 27, 24, 54]\n" + ] + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAD8CAYAAACcjGjIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xd8XNd14PHfncGgd6J3NrCTIAlSbKKoTlWqF9uS4iaXjddOvG6J193xxtl1shvHRbZiy0W9i5KiYrETLCDBBlYARAdRBr0NMDN3/8BApggQRBnMezNzvp8PPgRR3jsE78OZ285VWmuEEEKIi1mMDkAIIYT5SHIQQggxgiQHIYQQI0hyEEIIMYIkByGEECNIchBCCDGCJAchhBAjSHIQQggxgiQHIYQQI4QYcdOkpCSdl5dnxK0D15kzQ3/Om2dsHEKMl7TZCTt06FCL1jrZF/cyJDnk5eVRXFxsxK0D16ZNQ39u325kFEKMn7TZCVNKVfnqXjKsJIQQYgRJDkIIIUaQ5CCEEGIEQ+YchBBT53C6qG/vp6mzn8YuB+lx4azKSzQ6LBEgJDkI4WcaO/t5am8lf95fTUff4Ec+98n1eXzrlgWEhsiggJgaSQ5C+IneASfff/0kL5fU4nRrbl6Yxg0LU0mLDSc5JoxnD1bzuz2VHK5q4+cfW0F2YqTRIQs/JslBCD/Q3jvAJ39/kKM17TyyJpdPbZhJ7oyoj3zNd+9YxFUzZ/C1F49y+7/v5qUvrGVOSoxBEQt/J31PIUyuoaOP+39VRGl9J7/4+Eq+v2XxiMQwbPPiNLZ+aQM2q4VPP1VMe++Aj6MVgUKSgxAm1tjZz32/LKKho5+nPrmazYvTrvg9uTOi+PUjK2lo7+e/PX2YQZfbB5GKQCPJQQgT+9Gbp2judvDMZ9ewdvaMcX/fytwE/umeJewps/PDrSenMUIRqGTOQQiT2lvewhtH6/ny9XNZkhU34e+/b2UWZxu7eGJnBctz4rl7edY0RCkClfQchDChQZeb775WSnZiBF/YNHvS1/nG5vmsykvgO6+VUt/e58UIRaCT5HAZDqeLbWea+MdXjvMPrxznzWMNtPXI5J7wjaf2VnKuqZvv3L6IcJt10texWhT/5/4CXG7N1148itutvRilCGQyrHSJboeTH209ydZjDXQ7nESFWrEoxdP7q1EKNsxJ4l8fLCApOszoUEWAaurs59/eP8e185K5YUHKlK+XMyOS/3n7Qr718nH+UFTJ36yfOfUgRcCT5HCRho4+PvX7Ys42dnHvikw2L05j3ewkQiyKo7Ud7DjbzBM7y7n7F3v43d+sZk5KtNEhiwD0H9vKcDhdfPeORSilvHLNh1Zl827pBX7y9mk2zE2WtiuuSIaVPI7XdrDl53uoae3lyccK+el9y7hufirhNishVgsrcxP4+xvzee7xtfQNuLjnF3vYV2E3OmwRYNp6Bni+uJa7CjLJSxp9L8NkKKX453uXEm6z8o+vHEdrGV4SY5tyclBKZSultimlTimlSpVSX/ZGYL50trGLB35dhM1q4aUvrGPTvMt35Zdlx/PKF9eTEhvOo08e4GhNuw8jFYHuT/uq6Bt08dmNs7x+7ZTYcL6+eR77z7fyxrEGr19fBBZv9BycwFe11guANcB/U0ot9MJ1fcLt1vzDy8cJt1l45YvrmJd25XID2YmRvPC5tSRFh/LlZ0vocTh9EKkIdP2DLp4qqmTTvGTyU6en7MVDq3JYnBnLj988Sbe0WzGGKScHrXWD1vqw5/0u4BSQOdXr+soLh2oormrjW7cuICU2fNzflxAVys8eLKCqtZcfvCGbjMTUvVpSR0v3AI9PQ69hmNWi+P6di2nsdPDvH5ybtvsI/+fVOQelVB6wHNjvzetOF3u3g5+8fZrVMxO5f+XENwitmTWDL26azXPFNbx9XLrpYvLcbs1vdlWwODOWtbPGvxN6MlbmJnD/yiye3HWesqbuab2X8F9eSw5KqWjgJeArWuvOUT7/uFKqWClV3Nzc7K3bTsk/vXWaHoeTH9+1eNKrQr5yQz5Ls+L45svHudDR7+UIRbDYdqaJ8uYePnv1LK+tUBrLN26ZT0SolX9669S030v4J68kB6WUjaHE8Get9cujfY3W+gmtdaHWujA5Odkbt52SA+dbeelwLY9vnMXcKYzv2qwW/u9Dy+kbcPH/pJsuJul3eyrJjI/gtiXpPrlfUnQYn9s4iw9ON3GirsMn9xT+xRurlRTwJHBKa/2zqYfkG//2/lmSY8L422vnTvlaM5OiuK8wixeLa2nokBIFYmJq23rZU97CA4XZhFh9t7r8kbV5xISF8Mvt5T67p/Af3miJ64FHgOuUUkc8b7d64brT5khNO3vL7Xxmw0wiQidfmuBiX7hmNm6t+fWOCq9cTwSPlw7VAXDvSt+u44iLsPHI2lzeOtFAebPMPYiP8sZqpd1aa6W1Xqq1LvC8veWN4KbLL7eXERsewsfX5HrtmtmJkdy9PJNnDlTT3OXw2nVFYHO7NS8cqmH97CSyEnx/rOenNswkLMQivQcxQtDtkC5r6uKd0kYeW5dHdJh3q4d88do5DLrc/HaX9B7E+Ow7b6e2rY/7C40pp50UHcZDq3J4taSO2rZeQ2IQ5hR0yeGX2ysIt1n4m3V5Xr/2zKQo7liWwR/3VUkFVzEuLxTXEhMews2LrnzC23R5fOMslILf7JQXNeKvgio51LX38dqROh5alcOMaaqq+rfXzqF3wMVvd8uDJsbW2T/I2ycauHNZxpTKck9VRnwE9yzP4tmDNXLmtPhQUCWHJ3edB5iWujXD5qbGcPvSdH63p1LmHsSYth5toH/Qzf2F2UaHwic35OFwunm+uMboUIRJBE1y6B908eKhGm5Zkk5mfMS03uurN83D4XTzH9vKpvU+wr+9cKiGuSnRLJvEEaDeNj8tlqtmJvKHoipcciCQIIiSw9ZjDXT2O/nY6pxpv9fMpCgeKMzi6f3VMsknRlXT2ktJdTt3r8j0yY7o8XhsXR61bX1sO91kdCjCBIImOTy9v4pZyVGsmZXok/v99+vngoJ/e192TYuR3vLU4rp9SYbBkfzVjQtTSYsN56miSqNDESYQFMnh9IVODle387HVOT57lZYeF8Gja3J5+XAt5xq7fHJP4T/ePN7Asqw4cmb4fm/D5disFj6xJodd51pkU5wIjuTwzP5qQq0W7l3h27XkX7x2DpGhIfzsvbM+va8wt2p7L8dqO7htqW/qKE3EQ6tzCLVa+GNRldGhCIMFfHLoG3DxckkdtyxJIyEq1Kf3TowK5ZPr8/iv0gtSGll86E3PkNKtPiqyNxFJ0WHctjSdFw/VymFAQS7gk8PWY/V0+WgiejSPrcsj1GrhiZ1SnkAMefN4PQXZ8YaUyxiPR9fm0u1w8srhWqNDEQYK+OTw7MEaZidHsXqmbyaiL5UUHcaDq7J5paROznsQVLb0cKKuk9tNOKQ0rCA7niWZcfyhqAqtZVlrsAro5NDQ0cehqjbuWZFl6HLBz149C7eGJ2XXdNAbHlK6xYRDSsOUUjyyNpdzTd3sP99qdDjCIAGdHP7rxAUAbllsXN0aGKrYevvSdJ7eX01H76ChsQhjvXmsgRU58dO+EXOq7liaQVyETSamTcTXmxMDOjm8feIC81JjmJUcbXQofG7jbHoGXPxxX6XRoQiDVLb0cLKh05QT0ZeKCLXyQGEW75ReoLFThkPNoKnLt/8PAZscmrr6OVjZymaDew3DFmbEsmleMr/bU4nD6TI6HGGA9042AhhagXUiPrEmF6db8/T+aqNDEUCDj+csAzY5vFvaiNbmWi74yfUzsfcM8G5po9GhCAO8e/ICC9JjyU405yqlS+XOiOKa/GSeOVDNoMttdDhBr6FdkoNX/NeJC8xKiiI/1fghpWFXz0kiMz6CZw7IK7Fg09Lt4FBVGzctTDU6lAl5dG0uTV2OD+fvhHF8fT59QCaHtp4BiirsbF6cZpqiZgAWi+Lh1dnsLbdzvqXH6HCED31wqgm3Hqpf5E82zUshb0Ykv91VIctaDVYvPYepe+9kIy63NtWQ0rD7C7OxWhTPHpTeQzB592QjmfERLMqINTqUCbFaFJ/eMJOjtR0crGwzOpygdqFTeg5T9taJBrISzPkgpsaGc/38FF4srmXAKeO4waB3wMmuc83cuDDVVD3Z8bpvZTbxkTZ+I2ejG0p6DlPU1T/InrIWbjHZkNLFHr4qB3vPwIerV0Rg23WuBYfT7XdDSsMiQq08siaX9081ynCogWTOYYr2lLUw6NLcsMC8D+LGuclkxkfI0FKQeLe0kdjwEMNKuHjDI2tzsVksssvfIIMuN00+PnY44JLDttPNxISHsCI3wehQLstqUTy4Kptd51qotstJcYHM6XLzwelGrpufgs3qv49bSkw4dy3P4MVDtbT2DBgdTtBp7OzH1+sB/Le1jkJrzbYzTWzMTzb9g3jvyqGzJd44Vm9wJGI6Hapqo613kJv8ZOPbWD5z9Sz6B908vV9KaviaEUU7zf0bdIJK6ztp6nJw7bwUo0O5osz4CApzE3j9iCSHQLbtTDMhFsXVc5OMDmXK8lNjWDd7Bs8V1+D2cZ2fYFcvyWFqtp8ZOhj9mvxkgyMZnzuWZXCmsYszF+QY0UC142wzK3MTiAm3GR2KVzxQmE1Na59Ua/WxhnbfTkZDgCWHbWeaWZoVR3JMmNGhjMutS9KxKHjjqPQeAlFjZz+nGjrZ5Ac92fG6eVEaMWEhvFBcY3QoQaWho5+YsBCf3jNgkkNbzwAl1W1+9SAmx4SxbnYSrx+tl92nAWjH2WbAf3qy4xERauWOggzeOtFAV7+Un/eVho4+0uLCfXrPgEkOO88149Zw3Xz/SQ4Ady7LoLp16MB5EVh2nGkmJSaMBekxRofiVQ8UZtM/6GbrsQajQwkaDR39pPv4DJCASQ7bTjcxIyqUpZlxRocyITcvTsNmVbwuQ0sBxelys+tcM9fkJ5t2M+ZkLcuKY25KtAwt+VB9ez8Z0nOYOJdbs+Ps0INosfjXgxgXYeOa/BS2HquXFSAB5EhNO539Tq6ZFzhDSsOUUjxQmM3h6nbKmmQxxXRzOF20dDtIj5Oew4QdrW2nrXfQbx/EOwsyaOx0cKBSVoAEih1nm7EouHqOf7bJK7lreSZWi+KFQ7VGhxLwGjuGdkanx0vPYcJ2n2tBKbh6rn8+iDcsSCEy1MqrJXVGhyK8ZMfZZpbnJBAXGRhLWC+VHBPGtfNSeLWkzudnGweb4ZpK6TKsNHG7y1pYlBFLYlSo0aFMSmRoCJsXpfHm8Qb6B+UIUX/X0u3gWG0HmwJoldJo7lru6fHKnodpNXw8qF8OKyml/lMp1aSUOuGN601Ej8NJSXUbG/y8+373iky6+p385VST0aGIKdp1zrOE1U+HOcfr+vmpRIVaef2o9HinU72n55Dhp8NKvwc2e+laE3LgfCuDLs2GOf5dnmDd7CRSY8N4pUTGcP3drrMtJEaFsjjDv1bOTVREqJWbFqXx1vELcjbJNGpo7ycuwkZkqB9ugtNa7wQM6VvuLmshLMRCYZ55q7COh9Wi2FKQyfYzzVL10o9prdlT3sLa2TP8buXcZNxZkEFH3yA7PRv+hPc1dPT7fL4BAmDOYfe5FlblJRJusxodypTdvTwTp1uzVSq1+q3y5m4aOx1+35Mdrw1zkkiItPGa7NOZNg0dfYGdHJRSjyulipVSxc3N3nmV0dTVz5nGLtYHyIO4ID2W+WkxvHxYxnD91Z4yOwDrZwdGm7wSm9XCbUvTee/kBXocTqPDCUi1bX1kJvh2Mhp8mBy01k9orQu11oXJyd6ZqNtT1gIQEOWQh92zIpMjNe1UNHcbHYqYhN1lLWQlRJAzI9LoUHxmS0Em/YNuOfZ2GnT2D9LRN0h2gu/bk18PK+0+Zych0sbC9FijQ/GaLQWZWBS8Inse/I7T5WZfhT1ohpSGrcxJICMuXErATIO6tqGVSln+mhyUUs8ARcA8pVStUurT3rjuWLTW7C5rZt2cpICa+EuNDWdjfjLPHaxh0CUrQPzJifpOuvqdrAuy5GCxKO4oyGDnWVlM4W01rUPHCGcn+umwktb6Ya11utbaprXO0lo/6Y3rjiWQJ/4eWZNLU5dDuul+ZniYc93sGQZH4ntblg0tpnhTFlN4Va2/9xyMsPvc0IMYiMlh07wUMuMj+NM+OavXn+wpa2F+WgxJ0f5x2JQ3LUiPIT81mtfk2FuvqmnrJTLUSoIBZVj8NzmUtZA7I5LsxMCb+LNaFB+7Koe95Xapeukn+gddFFe1BczKuYlSamifTnFV24dDIWLqatv6yE6INKTsu18mh0GXm30VrQH9ID64KhubVfGnfdVGhyLGobiyjQGnOyB7suO1pSADQCamvai2rY8sA5axgp8mh2O17XQ7nFwdwA9iUnQYty5J56VDtfQOyPpxs9tT3kKIRbF6ZqLRoRgmKyGSVXkJvFpSJ8feeoHWmtrWXkkOE7HLU6J7bYBP/D2yJpcuh1PGcf3A3rIWCrLjifLxIfBms6Ugk3NN3Zxs6DQ6FL/X2eeky+E0bOjcL5PDnrIWlmTGER/pnyW6x2tlbgLz02JkYtrkOvsHOV7XEfAvVsbj1iXphFiUvKDxgpq2obkb6TmMU7fDSUl1e1CM7SqleGhVNqX1nZySV2KmVVzZilvD2lmSHBKjQrkmP5nXj9TLIUBTVPthcpCew7jsr7DjdPt/ie7xurMgE5tV8ZIcx2haReV2Qq0WVuT6d2Vgb7lreSYXOvspKrcbHYpfG97jYETpDPDD5DBcojtYHsTEqNCh4xiP1OOUHdOmVFRhZ3lOfEBUBvaGGxemEhdh47niGqND8Ws1rb3EhIUQG2HMPJbfJYc9ZS2snhkYJbrH676VWbR0O9h5Tmrmm01H7yCl9Z0y33CRcJuVu5dn8s6JC7RJOY1Jq23rIyvRmD0O4GfJobGzn7ON3UEzpDRs07wUEqNCeemQFOMzm/3n7WiZbxjhwVXZDLjcvHpE2uxk1bQZt4wV/Cw5DNeu2RBAJbrHIzTEwp3LMnjvZCMdvYNGhyMuUlRhJyzEQkFOvNGhmMqC9FiWZsXx3MEa2fMwCVprQzfAgZ8lh91lQ2fzLkgLnBLd43XfyiwGXG7ekMJmprKvopXCvATCQoJnmHO8HlyVzekLXRyr7TA6FL/T2jNA74DLsMlo8KPkoLVmT1kL64LkbN5LLcoYOiXuRVm1ZBptPQOcauiUIaXLuHNZBhE2K88elInpifprNVbpOVxReXMPjZ2OgK6nNBalFPeuyJJT4kxk//mhpZprJDmMKibcxq1L0nnjaL2UgJmgD5exGlhY1G+Sw/B8Q7CczTuaOwsyUApeld2nplBUbifCZmVplsw3XM5Dq7Ppdjh581iD0aH4leHd0UacHT3Mr5JDsJ3Ne6nU2HDWz07itSNS2MwMiirsFOYlEBriN4+RzxXmJjAzKUqGQyeotq2XuAgbseG+P8dhmF+0apdbUxSEZ/OOZktBBlX2Xkpq2o0OJai1dDs429gt+xuuYGg4NJP951vlnIcJqGntM+Ro0Iv5RXI4XtcRlGfzjmbz4jTCQiy8ViLrx420r2JovkEmo6/s7hVZKAUvHZbew3jVtvWSFW/sKIlfJIdgPpv3UjHhNm5YmMobxxoYlHIahtlXYScq1MqSzDijQzG9zPgI1s2ewUuHa3FLMb4rMsMeB/CT5LC3PHjP5h3N3QWZtPYMsEvKaRimqNzOqpmJhFj94hEy3L0rsqhp7eNgZavRoZheU5cDh9Nt+BHIpm/Z/YMuDlYG79m8o9mYn0xCpI1XS2TVkhGaOvspb+6RIaUJ2Lw4jahQqwwtjUNlSw8AeUlRhsZh+uRwqGrobN71c+RBHBYaYuG2pem8e/IC3Q5ZP+5rRcPzDTLMOW6RoSHcuiSdt45fkD0PV1BpH0oOM2dIchjTnrLhs3nlQbzYncsy6R90s+OMDC352r4KOzFhISzKkPmGibhvZRbdDifvlF4wOhRTO9/Si82qyIgPNzQO8yeHcjsF2fFEB/nZvJdakRNPTHgIO89KcvC1onI7q2cmYg3CMi5TsSovkcz4CN48JslhLJUtPWQnRho+n2Xq5NDVP8jx2nbpvo8ixGphw5wkdp5rlg1xPtTQ0UelvVfa5CRYLIqN+UlDpznKSrvLqrT3GD6kBCZPDgflbN4xbcxPpqGjn7ImqbXkK8P7G6Se0uSsn5NEl8PJsTqp1DoarTWV9h5yJTmMbV9Fq5zNO4aN+ckA7JChJZ8pKrcTF2FjYXrwlY33hnWe2mh7zrUYHIk5NXY66B90MzPJ+DJBJk8OdgrkbN7LyoyPYHZyFDvlQfOZooqh+YZgLBvvDYlRoSzKiGV3mbTZ0Zw3yTJWMHFy6Owf5ERdhwwpXcHVc5PZX2HHLfMO066uvY+a1j5pk1O0YU4Sh6vbcEmbHWF4GWueDCtd3sHzQ/MNMrY7tmvyk3E43XT2y9rx6Vbs2d27emaiwZH4t/Vzkhh0abqkzY5Qae8h1GohI97Y0hlg4uRQVG4nNMTCcjmbd0xXzUok1GqRs6V94GBlK1GhVuanxRgdil9bledps33SZi81tIw1whTLpE2bHPadt7NC5huuKDI0hFUzE2iXB23aFVe2sSI3wfD15/4uItTKytwEeUEzisqWXmaaYL4BTJocOvoGKa3vlCGlcdo4N5m+AScDTlk7Pl06+gY509jFSlk55xUb5ibRO+CUysIXcbuHlrGaYb4BvJQclFKblVJnlFJlSqlvTvV6B863omW+YdyGl7RK72H6lFS3ofXQkIiYuuHy+x19Mu8w7EJnPw6n2xQrlcALyUEpZQX+A7gFWAg8rJRaOJVr7quwExZioSBb5hvGY35aDGEhVuw9A0aHErAOVbVhtShpk16yJDMOq0XRKS9oPvRhwb1ASQ7AaqBMa12htR4AngW2TOWC+yrsrMhJkPmGcVJKMSM6lI6+QVq6HUaHE5AOVrayMD2WKKnx5RUhVguxETba+wal/ItHZcvQMapm6Tl4o6VnAjUX/b0WuGrM7zhzBjZtGvVTTrfmf1a2kpUQCU8bv5zLX2RWnSHe4aR34zUQa2w1x0Dj1vD3la2kxIbBS+Z4cAPBzNpz9A+66N+wkQh5Icj61l6e6+gnY585hi690XMYbc3ViJcCSqnHlVLFSqniwcHLdyU7+4c+Fxth80JowcOqFBalsHfL0JK39Q44cWtNTLi0SW8aXq4pS1qH9A+6CLdZRv2FagRv9BxqgeyL/p4FjDiiTGv9BPAEQGFhoWb79lEv9h9bT/KnfVUc+95NECKvJsZt0yY62vu4bfN32PX1aw0/YjCQPLurgh+9eYr9/3C99Mq8yLJpE2XV7fzqKz/niUcLjQ7HcF/61x3kzojiN2P9LJTvUoc3eg4HgblKqZlKqVDgIeD1yV5sX4WdlbkJhElimLAZnjO23zgmx4d606GqNrITI0iVxOBVCoiPsFFULiW83W5Nld08exzAC8lBa+0E/hZ4BzgFPK+1Lp3MtTp6BznZIPsbJis8xMKKnHhePyLJwVu01hysbKMw1xzjwIEmLtJGl8PJ0drgLuHdMLyM1SR7HMBL+xy01m9prfO11rO11j+e7HUOVMr+hqnaUpDJ6QtdnLnQZXQoAaG6tZeWbgeFebL5bTrEhttQCnYHeWXh883DBffMMxxsqh3Sw/sblmXL2byTdeuSdKwWxWtH6owOJSAUV7YBSM9hmtisiiWZcewuC+4zScqahl7MzUmNNjiSvzJdcpD5hqlJjglj/ZwkXj9aj9st68enqqSmjZiwEOammOehDTQb5iRRUt1OtyN4d0ufbeomLsJGsmfe0AxMkxxkvsF77l6eQW1bH4eq24wOxe+VVLezLDteDveZRhvmJuF0a/aV240OxTBljd3kp0ajfLga6UpMkxxkvsF7blqYRoTNyislMrQ0Fb0DTk5f6JKy8dNsZW4CETZr0J4Op7XmbFMXc1LMVQreNMlB5hu8JyoshJsXpfLmsQYcTpfR4fit47UduNxaksM0CwuxsnpmIrvOBee8g71ngPbeQdMNXZoqOch8g/fctTyTjr5Btp8JzgfOG47UtANQkC0rlabb1XOTKG/uoa69z+hQfO5s49BkdH6q9BxGkPkG79swJ4mk6FBelaGlSSupbidvRiSJUaFGhxLwrvGUnd91NvhezJQ1dQMw10QrlcAkyWF4vuEqOZvXa0KsFu5YlsFfTjVJ7ZpJ0FpzuLpNSnT7yJyUaNLjwtkRhMnhXGM3MeEhpMSYZ6USmCQ5FJV7zm+QsV2vunt5JgMuN28fbzA6FL/T0NFPU5eD5TkypOQLSik2zk1md1lL0JXSONvYxdwUc61UArMkB5lvmBZLMuOYlRwlq5YmoaR6aL5BJqN9Z2N+Ml39To7Wthsdik+VNXWbbr4BTJAc2noGONXQyVqZb/A6pRR3FWRyoLKV+iCc6JuKkuo2wkIszE+LNTqUoLFhThIWBTuCaBGFvduBvWeAOSZbqQQmSA77zw9tfFk7W5LDdNhSkIHWsFUqtU5ISU07SzLjCA0x/BEJGnGRNgqy49kRRHWWzn04GS09hxH2VbQSYbOyNEu679Mhd0YUy7LjeU0qtY7bgNPN8boOGVIywMb8ZI7VttMWJOehDyeHfJOtVAITJIeicjuFeQnyCm0abVmWQWl954fFvcTYTjV0MuB0y2S0ATbmJ6M1QbNbuqyxi+iwENJMeFaIob+R7d0OzjR2yf6GaXb70nQsCjnnYZxKPDWppOfge8uy4omLsAXNktazjd3MMeFKJTA4OeyraAVkvmG6pcSGs252Eq8drUdrqdR6JQer2kiLDSc9LsLoUIKO1aLYMDeJXeeag6KtnmvqNl3ZjGGGJoeiihaiQq0syZR6StPtzoIMquy9QX/i1pVoPVQdVF6wGOeauck0djo4HeAHVrX1DNDS7TDlMlYwOjmU21k1MxGbVeYbptvmxWmEhljkEKArONvYjb1nQJKDga6ZN1RKI9DrgpU1D01Gm+mAn4sZ9lu5qauf8uYemW/wkdhwG9fNS+GlL12eAAAX9klEQVSNow1BtwN1IvaWD02Eyr4b46TGhrMwPZZtZ5qMDmVaDRfck2GlS3w43yAPoc/ctTyTlm4HO4O0NPJ4FJXbyU6MIDvRPGf5BqPr5qdwqKotoOuCnWroJCYshAyTzm0ZlhyKyluICQthUYbsQPWV6xekkBQdyjMHaowOxZRcbs3+863ygsUErp2fjMutA/qMh9L6ThZmxJr2lEHDksPecjtXzZpBiMw3+IzNauHelVl8cLqJps5+o8MxnVMNnXT0Dcp8gwkUZCcQH2lj2+nATA4ut+ZUQyeLMsy7GMeQ38wDLjdV9l7WyUPocw+tysHl1rxwqNboUEynyHOG8dpZSQZHIqyWoSqtO8424XYH3pLWiuZu+gfdph45MSQ59DicAKyfIw+hr81MimLNrESeO1gTkA/dVBRV2JmVFEVanPl2qwaj6+an0NI9wIn6wFt+PfxvWmziZfyGJIfufidJ0aGmrCcSDB5enUN1ay9FFXajQzENp8vNgfOtMqRkIhvzk1EKPjgdeKuWSus6CQuxMDs5yuhQLsuY5OBwsnZ2kim3jAeDmxelERdh45kD1UaHYhrH6zo87VKSg1kkRoVSkB3PtgDc71Ba38n89FhTz7kaEpnTrVkvD6Fhwm1W7l6eybuljbQGSfXLK9nrmW+QfTfmcu28FI7VttPS7TA6FK/RWlNa32Hq+QYwcLXSutky32Ckj12Vw4DLzbMHpfcAsK/CTn5qNEnR5jrHN9hdOy8FrQPrAKDatj46+52SHEZjs1rImSGbjIyUnxrD1XOTeGpvJQPO4N4x7XC6OFjZKi9YTGhRRiypsWH85XSj0aF4zYk6z2S0iZexgkHJITosxIjbikt8esNMGjsdvHE0uEt5l1S30z/oltVzJmSxKK5fkMqOM804nC6jw/GK0vpOrBbFvDRzFtwbJskhiF2Tn0x+ajS/2VURFOWRL2dvWQsWBatnJhodihjFjQtS6RlwfbgPxd+V1ncwNyWacJvV6FDGZExyCJfkYAZKKT6zYRanL3SxpywwHrzJ2FNuZ4nnkBlhPmtnzyAy1Mp7JwNjaOmEp2yG2RmSHEJMWkskGG1ZnkFSdBi/2VVhdCiG6HY4OVrTLqvnTCzcZmXj3GTeP9Xo9xs3mzr7ae5ymLpsxjDzLrIVPhEWYuWxtbnsONv8YQnhYHLgvH1oabXMN5jajQtTaex0+P1u6dL6TgAWS89B+INPrMkl3Gbhd3vOGx2Kz+0psxMaYmFlboLRoYgxXDs/BYvC74eWSj3JLeCHlZRS9yulSpVSbqVUobeCEr6VEBXKlmWZvFpSH9D180ezt9zOypwE008OBrvEqFAK8xL9Pjkcr+sgd0YkMeHmn9+aas/hBHAPsNMLsQgDPbI2l75BFy8GUbVWe7eDUw2drJ8j8w3+4KaFqZy+0EVNa6/RoUyK1ppDVe2syPGPXuqUkoPW+pTW+oy3ghHGWZwZx4qceP60r8rvJ/3Ga7jw4DqZb/ALNyxIBfx3aKmmtY+WboffDGHKnIP40KNr8zjf0sPushajQ/GJPWV2osNCWGrissnir/KSopibEs07pReMDmVSiquGjkYOmOSglHpfKXVilLctE7mRUupxpVSxUqq4uTlw6qQEkluWpDEjKpQ/FFUZHYpP7C1v4aqZiaaujCk+6pbFaRysbKW5y/8K8RVXtRETFkJ+qrl3Rg+74lOhtb5Ba714lLfXJnIjrfUTWutCrXVhcnLy5CMW0yYsxMpDq7P54HQjtW3+Oa47XlX2HqrsvWyYK0NK/uTWpem4NX7ZezhU2cby3ASsfrLPS14yiY/42FW5APx5f2BXa93mOUDm2nkpBkciJmJeagyzkqJ4+0SD0aFMSEffIGebuij0kyElmPpS1ruVUrXAWuBNpdQ73glLGCUzPoLrF6Ty/MGagK7Wuu1MM7OSoshLMu9JXGIkpRS3LkmnqNyO3Y/OeCipbkNr/5lvgKmvVnpFa52ltQ7TWqdqrW/2VmDCOB+/Kgd7z4Bfdt3Ho2/ARVGFnU3Sa/BLtyxJw63hXT9atXSoqg2rRVGQHW90KOMmw0pihI1zk8lKiODP+wNzYnpveQsDTjfXzpe5L3+0MD2WvBmRvHXcf4aWiivbWJAeQ5QfVaSW5CBGsFgUH7sqh30VrZQ1dRsdjtdtO9NEZKhVSnT7KaUUtyxJZ2+5nTY/OObW6XJzpKadwlz/am+SHMSo7l+ZTYhF8cyBwJqY1lqz7XQz6+ckERYiJTP81W1L0nG5Ne+eNP/Q56mGLvoGXazwo/kGkOQgLiM5JoybF6fx4qFa+gcD4wQugLKmbura+2SVkp9blBFLdmIEbx43f3IY3vzmTyuVQJKDGMPHV+fQ0TfoV2O7V/LB8BJWmW/wa0opbluSwZ6yFlpMvmqpuKqNjLhwMuIjjA5lQiQ5iMtaO3sGs5Ki+NO+wJmY3namiflpMaTH+deDKka6e3kmLrc29RnoWmuKK1tZmedf8w0gyUGMQSnFx9fkcri6naM17UaHM2Wd/YMUV7Zx7XwZUgoE89JiWJgey6sldUaHcllnG7tp7HSwwQ8r/0pyEGN6oDCLmLAQfrvb/w8C2nGmGadbc50kh4Bxz4pMjtZ2UN5szlV1O84ODWNuzPe/YUxJDmJMMeE2HlqdzVvHG6hr7zM6nCl5p/QCSdGhflNPX1zZncsysCh45bA5ew87zjYzL9U/hzElOYgremxdHgBP7a00NI6pcDhdbD/TzI0LU/2m8Jm4spTYcNbPSeKVkjrTnUPSO+Dk4Pk2Nub7Z3FHSQ7iirISIrllcRrP7K+m2+E0OpxJ2Vtup9vh5KaFaUaHIrzsnhWZ1LX3UVzVZnQoH7Gvws6Ay801+f45jCnJQYzLZ66eRZfDyfMHa4wOZVLeLb1AVKiVdX44MSjGdtPCNCJsVl4pMdcRtzvONBNhs1KY55/DmJIcxLgUZMdTmJvAf+45j9PlX9VaXW7Neycb2TQ/RXZFB6CosBA2L05j67EGU23Y3HmuhbWzZxBu8882J8lBjNsXNs2mtq2Pn7131uhQJqSkuo2W7gFuXiRDSoHqvpVZdPU7TXPOQ5W9h/MtPWz048OkJDmIcbt+QSoPr87mF9vL+eC0/5RLfqf0Ajar4tp5/recUIzP2lkzmJkUxZ/2maMW2M6zQ0chX+PHZVokOYgJ+e4di1iYHsvfPXfUL44S1Vrz7slG1s1OIibcZnQ4YppYLIqPX5XDoao2TjV0Gh0OO842k5MYSd6MSKNDmTRJDmJCwm1WfvHxFbjdmr99usT0p8Wdaeyiyt4rQ0pB4N4VWYSGWHja4CNuHU4Xe8vtXJOfjFL+u2xakoOYsLykKH5631KO1LTzy+3lRoczpq1HG7AouGGh/3bvxfgkRIVy+5J0Ximpo8fAJdc7z7bQO+DiugX+3eYkOYhJuWVJOrctTecX28uoaTXn8JLLrXnxUC0b85NJiQk3OhzhAx9fk0u3w8lrR4wrxvfqkToSo0LZMMd/J6NBkoOYgm/ftgCLUvxg60mjQxnV7rIWLnT280BhttGhCB9ZkRPP/LQY/ry/Cq19v2O62+Hk/ZON3L40HZvVv3+9+nf0wlDpcRF86fo5vHeykW1nmowOZ4Tni2tIiLRxvZ9378X4KaX4xJpcSus7OVzt+0rC75y4gMPpZktBhs/v7W2SHMSUfGbDLGYlRfH910txOM2zAamtZ4D3ShvZUpApG9+CzN3LM4mLsPHrHb6fD3v1SB3ZiREBUdxRkoOYktAQC9+7cxGV9l6e2FFhdDgfeu1IHQMutwwpBaGosBAeW5fHuycbOdvY5bP7Nnc52FPWwpZlmX69SmmYJAcxZRvzk7ltSTr//kEZ53z4MI7lhUO1LM6MZWFGrNGhCAN8cl0ekaFWfuXD1XRbj9Xj1gTEkBJIchBe8r07FxEZZuXrLx3DZXDp5BN1HZTWd3L/Suk1BKuEqFAeXp3Da0frfbaa7tUj9SxMj2VuaoxP7jfdJDkIr0iOCeO7dyykpLqd3xt87sPzxTWEWi0B8wpOTM5nr56FRcETO6d/uPN8Sw9Ha9q5a3ngtDlJDsJr7irI5Np5yfzvd85QbTdm70Nn/yAvHarltqXpxEeGGhKDMIe0uHDuXZHFc8U1NHX1T+u9fr/nPCEWxZaCzGm9jy9JchBeo5Tix3cvwWpRfPWFIwwaUNr7+YM19Ay4+NT6mT6/tzCfz18zG6fLPa07+Vu6HTx7sIa7l2eSGhs4my0lOQivyoiP4Md3L+ZgZRs/9PHmOJdb8/u9lazKS2BJVpxP7y3MKS8pigdX5fDHoirKmqZnscTv9pxnwOXm85tmT8v1jSLJQXjdloJMHt84iz8UVfHMAd8VQXv/VCO1bX3SaxAf8T9uyici1Mr33zjp9V3TXf2D/KGois2L0pidHO3VaxtNkoOYFt/YPJ+r5ybxnddOUFzZ6pN7/ufu82TGR3DjwlSf3E/4hxnRYfzdDfnsOtfC+6e8u5P/T/uq6ep38sVNc7x6XTOQ5CCmhdWi+PnDK8iMj+DzfzpMa8/AtN6vtL6D/edbeWxdLiF+XtNGeN8ja3OZkxLNj9486bWd/P2DLp7cfZ6r5yYF5DCmPEVi2sRF2vjlJ1bS0TfAd147Ma33+s3OCiJDrTxYmDOt9xH+yWa18J3bF1Jl7+U3Xlra+vT+alq6HQHZawBJDmKaLUiP5Ss35LP1WANvHpue831/v+c8rx6p59G1ecRFymlvYnQb85O5ZXEa//b+OQ5VTW2o83xLD//yzhmunpvEmlmJXorQXKaUHJRS/6KUOq2UOqaUekUpFe+twETg+NzGWSzLiuPbrx6nucvh1Wu/frSe7289yc2LUvnazfO8em0ReP7XvUvJiI/gi38+TEv35Nqi0+Xm758/gs2q+Jf7lgVEHaXRTLXn8B6wWGu9FDgLfGvqIYlAE2K18L/vX0bPgItvv3rcaytGdp1r5qvPH2FVXiL/96HlWC2B+ZAK74mLsPHLT6ygvXeQLz1dgnMSe3F+vbOCkup2fnjXYtLiAmdfw6WmlBy01u9qrYfP49sHZE09JBGI5qbG8NUb83mntJHvv3ES9xTrL51q6OTzfzzE7ORofvNoIeE2KcstxmdRRhw/umsxRRV2fvrOmQm9WDlR18G/vneW25amc+eywCmVMZoQL17rU8BzXryeCDCfvXoWTV0Ontx9nuZuBz97YNmkzlpo6XbwmaeKiQ4P4fefXE1chMwziIm5vzCbIzXtPLGzgo7eQX5w16IrtsUjNe08/odiEqNC+dGWxQE7nDTsislBKfU+kDbKp/5Ra/2a52v+EXACfx7jOo8DjwPk5MiKkmBksSi+fdsCUmLC+Mnbp2nrGeCXn1g5oV/uDqeLz/3xEPYeB89/bm1Ad+vF9PrhlsUkRoXy7x+UUd7cza8eWUlSdNioX/tKSS3feOk4KTFh/PaxQhKiAr9ul5rq+K9S6jHg88D1WutxVVsrLCzUxcXFU7qvuMSmTUN/bt9uZBTj9vLhWr7+4jFmRIfyk3uWcN38K29c01rz1ReO8vLhOn7+seXcvjSwu/UBzyRt9o2j9XztxaPERdi4d0UWNy9KY2lWHJ39Tg5Xt/FuaSPPHKhmzaxEfvHxlSQamBiUUoe01oW+uNeUhpWUUpuBbwDXjDcxCAFwz4os5qRE87UXjvGp3xdzz4pMvnP7wstWUnW7NT/YepKXD9fxlRvmSmIQXnPHsgxmJkXxT2+d4tc7K/jF9nLiI2109A2i9dCGzsfW5vLt2xdiC6INllOdc/g5EAa85xl/26e1/vyUoxJBYWlWPK9/aT0//6CMX2wvZ/uZZr65eT73rczCctHKI5db882XjvHCoVo+vWEmX75+roFRi0C0ODOOpz+7hraeAT443cTecju5MyIpzEugIDueyFBvTs/6hyn9i7XWgbk1UPhMWIiVr940j1sWp/Od107w9ZeO8fSBaj63cRZxETYiQq08ufs8W4818N+vn8vf3TA34CcChXESokK5d2UW966UhZfBlw6FKS3MiOWFz6/l5cN1/OTtU3zhz4c/8vl/uHU+j28MrJLIQpiZJAdhGkop7l2ZxebFaZQ1ddM36KJvwEVCVCgF2bL5XghfkuQgTCcqLIRlkgyEMFTwTL0LIYQYN0kOQgghRpDkIIQQYgRJDkIIIUaQ5CCEEGIESQ5CCCFGkOQghBBiBEkOQgghRpDkIIQQYgRJDkIIIUaQ5CCEEGIESQ5CCCFGkOQghBBiBEkOQgghRlBaa9/fVKku4IzPbzwxSUCL0UGMg8TpXRKn9/hDjMP8JdZ5WusYX9zIqPMczmitCw2697gopYrNHiNInN4mcXqPP8Q4zF9iVUoV++peMqwkhBBiBEkOQgghRjAqOTxh0H0nwh9iBInT2yRO7/GHGIf5S6w+i9OQCWkhhBDmJsNKQgghRtJaX/ENyAa2AaeAUuDLno8nAu8B5zx/Jng+PsPz9d3Azy+51sPAceAY8F9A0mXuudLzdWXA/+OvvZyNwGHACdxnxhgv+vx9gAYKzRgn8K/AEc/bWaDd4Dh/DNQA3Zd8PAx4zhP/fiBvmuJ80BNjKfDTMZ4HI9vmlGL0Yduc6s/SV21zvHFOuG1OMtYbgUOen8kh4Lrx/p9Otn1e9t98pS/wXDQdWOF5P8bzn7UQ+CnwTc/Hvwn8s+f9KGAD8PmL/zMYWjrbhOeXg+f7v3eZex4A1gIKeBu4xfPxPGAp8Ac++gCaJsaLYtgJ7OOjD6Cp4rzoa74E/KfBca7x3PfSB/CLwK887z8EPDcNcc4AqoFkz9+fAq43Wduccow+apteidMHbXMicU64bU4y1uVAhuf9xUDdRH5Wk2mfl3sb17CS1rpBa33Y834XQ1kwE9ji+YEO/2Dv8nxNj9Z6N9B/yaWU5y1KKaWAWKD+0vsppdKBWK11kR76V/3homtXaq2PAW6zxujxQ4YawEeub8I4hz0MPGNUnJ5r7NNaN4zyqYvv+SJwveda3oxzFnBWa93s+fv7wL2XBmJw25xyjB7T3Ta9Feew6Wqb44rTc40Jt81JxlqitR5+PkqBcKVU2Hh/VpNpn5cz4TkHpVQeQ9ltP5A6/APz/Jky1vdqrQeBLzDU5alnKIM+OcqXZgK1F/291vMxv4hRKbUcyNZabzVznBfFkQvMBD4wMM6xZDLUpUdr7QQ6GHrV57U4GeqCz1dK5SmlQhh6oLIvE4shbdMbMfqibXojzovimLa2OYE4xzKutjnJWO8FSrTWDsbf7qbUPi82oeSglIoGXgK+orXunOjNlFI2hn5RLAcyGBrr+9ZoXzrKx7Q/xKiUsjA0XvpVM8d5yd8fAl7UWrsMjHPMy4zysY/8G6Yap9a6zRPnc8AuoJKhsdkJx3I5Rsfoq7bp5Z/ltLXNCcQ5lnG1h4nGqpRaBPwz8LmJ3GcCX3dF404Onof8JeDPWuuXPR9u9HRjhrszTVe4TAGA1rrc0+V5HlinlLIqpY543n7AULbLuuj7srjMUIQJY4xhaKxwu1KqkqGxyteVUh9uzTdJnBd7iIu67QbFOZZaPK/oPK/w4oBWL8eJ1voNrfVVWuu1DNX+OmeytjnVGH3VNr35s5zOtjneOMcyZtucTKxKqSzgFeBRrXX5RfcZ8bPyVvsczbiSg2cM7UnglNb6Zxd96nXgMc/7jwGvXeFSdcBCpVSy5+83eq7p0loXeN6+4+lmdSml1nju/eiVrm2WGLXWHVrrJK11ntY6j6FJvzu11sVmivOin9s8IAEouvjivo7zCte4+J73AR94Eo0340QpleL5M4GhicbfmqxtTilGH7ZNr/wsfdA2xxXnFS5x2bY5mViVUvHAm8C3tNZ7hr94jP/TKbfPy9JXmLH2/Ds3MNQ1OcZfl5fdytDY2l8YWo71FyDxou+pZCiDdjOUzRZ6Pv55hiZljgFvADMuc89C4ARQDvycvy7HWuW5Xg9gB0rNFuMlX7Odj64IMVWcwPeA/2WS//Ofer7P7fnze56PhwMvMDRGfACYNU1xPgOc9Lw9NMbzYGTbnFKMPmybU44T37TN8cY54bY5mViBb3vaz5GL3lLG+386mfZ5uTfZIS2EEGIE2SEthBBiBEkOQgghRpDkIIQQYgRJDkIIIUaQ5CCEEGIESQ5CCCFGkOQghBBiBEkOQgghRvj/433fobcxmmMAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "raw_data_period=[]\n", + "raw_data_period_mean=[]\n", + "raw_data_date=[]\n", + "raw_data_value=[]\n", + "\n", + "for i in range(len(raw_data[\"smooth_b\"])):\n", + " if abs(raw_data[\"smooth_b\"][i])<=0.1:\n", + " #print(i, raw_data[\"smooth_b\"][i])\n", + " #print(i, raw_data[\"convert\"][i])\n", + " raw_data_period.append(i)\n", + " raw_data_date.append(raw_data[\"convert\"][i])\n", + " raw_data_value.append(raw_data[\"smooth_b\"][i])\n", + " \n", + " \n", + "#print(raw_data_period) \n", + "#print(len(raw_data_period))\n", + "\n", + "for i in range(len(raw_data_period)-1):\n", + " #print(i)#,raw_data_period[i+1]-raw_data_period[i])\n", + " raw_data_period_mean.append(raw_data_period[i+1]-raw_data_period[i])\n", + " \n", + "print(raw_data_period_mean)\n", + "2.*np.mean(raw_data_period_mean)\n", "\n", "fig, ax = plt.subplots()\n", "ax.plot(raw_data[\"convert\"],raw_data[\"smooth_b\"])\n", + "plt.axvline(x=raw_data_date[-1], color='r', linestyle='-')\n", + "plt.axvline(x=raw_data_date[-2], color='r', linestyle='-')\n", + "plt.axhline(y=0, color='r', linestyle='-')\n", "ax.set_xlim(\"2018-01-01\",'2020-01-01')\n", + "ax.set_ylim(-2.5,2.5)\n", "\n", - "#plt.show()\n", - "raw_data_period=[]\n", - "raw_data_period.append(np.where(raw_data[\"smooth_b\"]== 0.0))\n", - " \n", - " \n", - "#for i in range(len(raw_data[\"smooth_b\"])):\n", - "# if raw_data[\"smooth_b\"][i]==0.:\n", - "# raw_data_period=raw_data[\"convert\"][i]\n", - "#print(raw_data[\"smooth_b\"]) \n", - "#print(raw_data_period)\n", - "print(np.where(raw_data[\"smooth_b\"]<= 0.01))" + "plt.show()" ] }, { -- 2.18.1