{ "cells": [ { "cell_type": "markdown", "id": "ea8884ce", "metadata": {}, "source": [ "# Thermogram Fitting\n", "\n", "This notebook provides a example how to fit a thermogram for a multi-stage decomposition. Polyvinylchloride pyrolysis shows a characteristic 2-stage decomposition:" ] }, { "cell_type": "code", "execution_count": 1, "id": "57314c14", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnwAAAGwCAYAAAAzL7gEAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAtphJREFUeJzs3Xd4k+X6wPFvkjZN994ttMxSRgsFShmyKkVx4EQcIPLT40bxqIBAUfSgIoooHpygx6MgDlzYIzIEoYIUyiqbFlrohO7dJL8/QqO1uyRNx/25rly2b573fe8gJHeecT8KvV6vRwghhBBCdFhKSwcghBBCCCHMSxI+IYQQQogOThI+IYQQQogOThI+IYQQQogOThI+IYQQQogOThI+IYQQQogOThI+IYQQQogOzsrSAbS2qqoq9u/fj7e3N0ql5LtCCCFEe6DT6cjMzGTgwIFYWXW69OWKdbo/sf379zN06FBLhyGEEEKIFtizZw9DhgyxdBjtTqdL+Ly9vQHDXxhfX18LRyOEEEKIpkhPT2fo0KHGz3HRPJ0u4asexvX19SUgIMDC0QghhBCiOWQ6VsvIn5oQQgghRAcnCZ8QQgghRAcnCZ8QQgghRAcnCZ8QQgghRAcnCZ8QQgghRAcnCZ8QQgghRAcnCZ8QQgghRAcnCZ8QQgghRAcnCZ8QQgghRAcnCZ8QQgghRAcnCZ8QQgghRAdn0YRv+/btXH/99fj5+aFQKNiwYUOj52zbto1BgwZhY2NDjx49WLNmjdnjFEIIITqTe++9l8mTJzfrnKCgIJYvX26WeMSVs2jCV1xcTFhYGCtXrmxS++TkZCZNmsTYsWNJTEzkiSee4P/+7//43//+Z+ZIhRBCdGT33nsvCoWCl19+ucbxDRs2oFAoLBRV44YNG8aDDz5Y49iqVatQKBS1OkTuvfdeRo0a1aTrvvnmmybvUElJSUGhUJCYmGjS64qmsWjCd8011/Diiy9y0003Nan9qlWrCA4OZtmyZfTp04dHH32UW2+9lTfeeMPMkbZder2eSq2OwrJK8ksr0en0dbapqNJRUlFFfmklOUXlZBWWkVVQRmZBGen5pVzIKyUtt4TUS4bH2YvFpOQUcya7iNPZRZzKKuJkZiEnMgs5llHA0fQCDqXlcyA1jzPZReSX1H1vIYRoLzQaDa+88gq5ubkmvW5FRYVJr/dXY8eOZdu2bTWObd26lcDAwFrHt23bxrhx45p0XWdnZ1xcXEwTpIWtXLmSoKAgNBoNkZGR7Nmzp8H269evJyQkBI1GQ//+/dm4cWON57/++msmTJiAu7t7gwlsfHw848aNw97eHicnJ6666ipKS0tN9bKazcpid26B+Ph4oqOjaxyLiYnhiSeeqPec8vJyysvLjb8XFhbWalOl1VGh1VGp1VOl1VGlMyRIF4sruFhUTqVWR4VWT2WVoV1eSSU6vR4nW2sqqgzJVnmVzniuVqenSme4ll4PVioF5ZU6Siu1VOn02Fqr0Or0lFdpKa/SUVKhpVKrQ6fXo9OBTq9Hrzf8t+bPoNXp0ev1lFZqKa00nK//S56lUipwsbWmUquj7HJMrZmHWasUWKuUxoetWomdtRVOtlaorQzHrJRK1FYKbKxUuNmrcXdQ46SxxsZKidpKiYONFT7OGvycbXGxs27T365b26LvjpCcU8yKOwbibGdt6XCEaFT1+5Ul2FqrmvX+ER0dzalTp1iyZAmvvvpqve2++uorFi5cyKlTp/D19eWxxx7jqaeeMj4fFBTEzJkzOXnyJBs2bODmm29mzJgxPPHEE3z66ac89dRTpKamcu211/LJJ5+wfv16YmNjyc/P55577uGNN95ApVI1KeaxY8fy8ssvk5GRgY+PDwC//vorCxcurPEakpOTOXv2LGPHjgUgNTWVp556ip9//hmlUsmoUaN48803CQoKAgy9gXl5ecapVoWFhTz44INs2LABJycnnnnmGb799lvCw8NrDOOWlJRw3333sX79elxdXZk/fz4PPPAAAMHBwQAMHDgQgNGjR9dKSk1t3bp1zJ49m1WrVhEZGcny5cuJiYnh+PHjeHl51Wq/a9cupk6dypIlS7juuuv47LPPmDx5Mvv27aNfv36AYXRy5MiR3H777dx///113jc+Pp6JEycyd+5c3nrrLaysrDhw4ABKpeX62dpVwpeRkYG3t3eNY97e3hQUFFBaWoqtrW2tc5YsWcLzzz9f6/jjn++nTJ1MVmEZOUXm+/bV2rQ6PReLG389CgUoAKVCYfhZoUBx+bjy8s/Ky40UgFJZ/bwCpQJAgZXS8HNhWRWF5VUAVGr1VGq1gGne4G2slDjbWuOoscLVTo23swZPBxucNFY426kJ9rCju6cDAa52qJQdOzHMKihjza4UAL7Zn8a9I4ItG5AQTVBaqSV0oWWm3SS9EIOduukfcyqVin/961/ceeedPP744wQEBNRqk5CQwO23386iRYuYMmUKu3bt4uGHH8bd3Z17773X2O61115j4cKFxMbGArBjxw5KSkpYsWIFa9eupbCwkJtvvpmbbroJFxcXNm7cyJkzZ7jlllsYMWIEU6ZMaVLMI0aMwNramq1btzJ16lSSkpIoLS1l5syZPPvssyQnJxMcHMzWrVvRaDRERUVRWVlJTEwMUVFR7NixAysrK1588UUmTpzIwYMHUavVte4ze/Zsdu7cyXfffYe3tzcLFy5k3759hIeH12i3bNkyFi9ezLx58/jyyy956KGHGD16NL1792bPnj0MHTqUX375hb59+9Z5H1N7/fXXuf/++5kxYwZgGCn88ccf+eijj5gzZ06t9m+++SYTJ07k6aefBmDx4sVs2rSJt99+m1WrVgFwzz33AIYh6vo8+eSTPP744zXu0bt3b1O9rBZpVwlfS8ydO5fZs2cbfz9//jyhoaHsSb6ElVPdmbZKaUhm3OzVeDjYoLE29EpZWylRqxQ42VpjpVSQV1KJjbUKR40VNsbeK8O5ViolKqUhmarS6tFYK7G1VqFUKiit0GKlVKCxVmFz+biNlcqYbCkvJ1XVyZVS+dffDT/bWKmwUxvOt1EZ/qtQQG5xJXmlFVgpldiqVVgrFYbXo1JirVIYXodKYfJes4oqHQVllVRqdVRW6S/3mBp6NUvKtcbnKqp0l5NCHWWVWi4VV5BdVE5xeRXlVTrKKw3Xycgv42JxBeVVOrIKy8kqLAeK672/WqUk2MOe7l72dPd0MD66edpjb9Mx/pqnXCwx/nzwfL4FIxGi47rpppsIDw8nNjaWDz/8sNbzr7/+OuPHj2fBggUA9OrVi6SkJJYuXVoj4Rs3blyNXr8dO3ZQWVnJv//9b7p37w7Arbfeyn/+8x8yMzNxcHAgNDSUsWPHsnXr1iYnfPb29gwdOpRt27YxdepUtm3bxsiRI7GxsWH48OFs27aN4OBgtm3bRlRUFDY2Nnz66afodDo++OAD42fB6tWrcXFxYdu2bUyYMKHGPQoLC/n444/57LPPGD9+vLG9n59frXiuvfZaHn74YQCeffZZ3njjDbZu3Urv3r3x9PQEwN3d3dgb2RKFhYUUFBQYf7exscHGxqZWu4qKChISEpg7d67xmFKpJDo6mvj4+DqvHR8fXyNnAMNIYlMWlVbLyspi9+7d3HXXXQwfPpzTp08TEhLCSy+9xMiRI5t8HVNrV5+EPj4+ZGZm1jiWmZmJk5NTnb17UPsvQvVfkvmT+tC7RxDejhq8nGywU6uwUhoSNmU77inycVbh46xp9fuqrZR4ONT+B3clyiq1ZBeWU1BWSUFpFbklFaTnl3GpuJyC0iouFpdzJruY5Jxiyqt0HM8s5Hhm7SF7P2cN3b2qk0B7uns50NPLEU9H08ZrTjtOZrP7zCXevCOcWWsTOZQmCZ9oH2ytVSS9EGOxe7fEK6+8wrhx4/jnP/9Z67mjR49y44031jg2YsQIli9fjlarNQ7FDh48uNa5dnZ2xmQPDCNUQUFBODg41DiWlZXVrHjHjBnD+vXrAcM8vTFjxgB/DpnOmDGDbdu2GYcfDxw4wKlTp3B0dKxxnbKyMk6fPl3r+mfOnKGyspKhQ4cajzk7O9fZYzVgwADjzwqFAh8fn2a/nsaEhobW+D02NpZFixbVapeTk4NWq61zZPDYsWN1Xru+kcSMjIwmx3fmzBkAFi1axGuvvUZ4eDiffPIJ48eP5/Dhw/Ts2bPJ1zKldpXwRUVF1Zo8uWnTJqKiopp9revC/AgIqD1+L9oOjbWKQDe7RttpdXou5JVyKruI01lFnM4u5nR2EWeyi8gpquBCfhkX8svYcTKnxnneTjb083MmIsiVqG7uDAhwaZPDwpeKK7jnQ8Mk4zfvCAfgVHYRZZVaNC38QBOitSgUimYNq7YFV111FTExMcydO7dGr11z2Nvb1zpmbV1z3q1CoajzmE6na9a9xo4dy0svvcT58+fZtm2bMVEdPXo07777LqdPnyY1NdW4YKOoqIiIiAj++9//1rpWdS9cS5ni9TQmKSkJf39/4+919e5ZUvXr/cc//mEcSh44cCCbN2/mo48+YsmSJRaJy6L/CouKijh16pTx9+TkZBITE3Fzc6NLly7MnTuX8+fP88knnwDw4IMP8vbbb/PMM89w3333sWXLFr744gt+/PFHS70E0QaolAoC3ewIdLNjbO+aSXxeSYUhAcwyrDY+nV3EicwiUnNLyCwoJ7Mgi83HDN8+XeysGdHDg2HBbkSHeuPrXHevcWtLTDWsGPRytOGGMD+eXn+QCq2OnKJyAlwbT4iFEM338ssvEx4eXqsXq0+fPuzcubPGsZ07d9KrV68mL7QwteHDh6NWq3nnnXcoKysjIiICgCFDhpCdnc1HH31kHPoFGDRoEOvWrcPLywsnJ6dGr9+tWzesra35448/6NKlCwD5+fmcOHGCq666qslxVs/Z02qvbI63o6Njk+L28PBApVLVOTJY35ByfSOJzRmC9vX1BWr3RPbp04dz5841+TqmZtGEb+/evcYVQ4Bx3Hz69OmsWbOG9PT0Gn84wcHB/Pjjjzz55JO8+eabBAQE8MEHHxATY5nhAtH2udipieiqJqKra43jJRVVHLlQwMG0fPYkXyT+9EXySir58WA6Px5MJ/a7I4zu5cntgwMZ38cbtZXlVlYdOW+YhjC8u6EEQFR3d+OKbiGEefTv35+77rqLFStW1Dj+1FNPMWTIEBYvXsyUKVOIj4/n7bff5p133jFLHH/v+KiLra0tw4YN46233mLEiBHGxFOtVtc4Xt37dtddd7F06VJuvPFGXnjhBQICAjh79ixff/01zzzzTK3FKo6OjkyfPp2nn34aNzc3vLy8iI2NRalUNms+uJeXF7a2tsTFxREQEIBGo8HZ2bkFfypNo1ariYiIYPPmzcYi0jqdjs2bN/Poo4/WeU5UVBSbN2+uUf2juSOJQUFB+Pn5cfz48RrHT5w4wTXXXNPs12EqFk34xowZg15ff82Quoo+jhkzhv3795sxKtEZ2KmtGBLkxpAgN2aODKZKqyMxNY/fTuWw81QOf6TksvV4NluPZ+PpaMM9w7pyV2QX3E08T7EpquclhvoZvtF+fN/QhpoLIUzkhRdeYN26dTWODRo0iC+++IKFCxeyePFifH19eeGFF1o89NuYv3d81Gfs2LFs377dOH+v2ujRo9m6dWuNzhU7Ozu2b9/Os88+y80330xhYSH+/v6MHz++3p6z119/nQcffJDrrrvOWJYlNTUVjabpc8atrKxYsWIFL7zwAgsXLmTUqFFmL8sye/Zspk+fzuDBgxk6dCjLly+nuLjYONQ6bdo0/P39jcOss2bNYvTo0SxbtoxJkyaxdu1a9u7dy3vvvWe85qVLlzh37hwXLlwAMCZ2Pj4++Pj4oFAoePrpp4mNjSUsLIzw8HA+/vhjjh07xpdffmnW19sgfSeTmpqqB/SpqamWDkW0YaezCvWv/HRUP+TFTfquz/6g7/rsD/qQ+T/pl8Yd0xeUVrRqLFPe3aXv+uwP+m8Tz7fqfYUQoj5FRUV6Z2dn/QcffNBq92zp5/dbb72l79Kli16tVuuHDh2q//33343PjR49Wj99+vQa7b/44gt9r1699Gq1Wt+3b1/9jz/+WOP51atX64Faj9jY2BrtlixZog8ICNDb2dnpo6Ki9Dt27GhW3Kam0Osb6GLrgNLS0ggMDCQ1NbXOGktC/FWlVsfGQ+l8sCOZQ5dLobjZq3lsXA/uHtYVa5X5h3pj3tjO8cxCPp0ZycieHsbjer1eilILIVrF/v37OXbsGEOHDiU/P58XXniBbdu2cerUKTw8PBq/gAnI5/eVsejWakK0ddYqJTeG+/PdoyN4954Iunnac6m4gue/T+Lmd3Zxoo4yMKZWXUjbzd4w4fntLScZsOh/vPxT3WUFhBDCHF577TXCwsKIjo6muLiYHTt2tFqyJ65c+1orL4SFKBQKYvr6MD7Eiy/2pvFK3DEOnc/n+rd+48XJ/bhtcKDZ7v3+tAhyiiro6m5YkauxVlFQVkVqbkkjZwohhGkMHDiQhIQES4chroD08AnRDFYqJXdGduHnJ6/iql6elFfpePrLg8z9+hDlVebZL3RgF1euDvU27hjS1d1Q3+vsRUn4hBBCNI0kfEK0gLeThjX3DuHJ6F4oFPD5nnPcviqenKJys9876HJP37mLJQ2uchdCCCGqScInRAsplQpmRfdk9b1DcLGz5kBaPvd8uIe8kgqT3WPfuVze/fU0CWdzjceqdx8pLK8it6TSZPcSQgjRcUnCJ8QVGtPbi68fGo6Hgw1H0wuYvvoPCstMk4htPZbFkp+O8WVCmvGYxlqF7+X9klMuFpvkPkIIITo2SfiEMIFung789/8iDT19qXnMXLOX0oorn9N3MrMIgB5eDjWOd3H7c1hXCCGEaIwkfEKYSG8fR/5zXySONlbsSbnEo5/tQ6e7sjl2J7MMZV96/i3hGxDgzNjenvg4N73KvRBCiM5LEj4hTKh/gDOrZwzBxkrJ5mNZvLPtVIuvVVapJeVyD15P75oJ33OTQlk9YyjDurlfUbxCCIN7770XhUKBQqHA2toab29vrr76aj766CN0Oh3btm0zPl/fo3qbsIqKCpYuXcqgQYOwt7fH2dmZsLAw5s+fb9yO66/i4+NRqVRMmjSplV+16Ewk4RPCxAYHubF4cj8Alv9ykkNp+S26zs5TOWh1erwcbfBxkp48Icxt4sSJpKenk5KSwk8//cTYsWOZNWsW1113HcOHDyc9Pd34uP32243tqx/Dhw+nvLycq6++mn/961/ce++9bN++nUOHDrFixQpycnJ46623at33ww8/5LHHHmP79u11JoRCmIIUXhbCDG6LCGDb8Sw2HsrgqfWJ/PDYKNRWzft+lZ5fBkBYoEudW6jp9XpyiirQWCtx1FibJG4hOjMbGxt8fHwA8Pf3Z9CgQQwbNozx48fzySef8H//93/Gtra2tpSXlxvbV3v55Zf57bff2Lt3LwMHDjQe79KlC6NHj65VSqmoqIh169axd+9eMjIyWLNmDfPmzTPjqxSdlfTwCWEGCoWCFyf3x8NBzYnMIt7bfrrZ16gu7+JqV3cy98B/Ehjy0i/8dDjjimIVojWUVFTV+yir1Jq8ramMGzeOsLAwvv766ya1//zzz7n66qtrJHt/9fcvb1988QUhISH07t2bu+++m48++kjqawqzkB4+IczEzV7NgutCmbU2kRVbTjFpgB/BHvZNPn/a8CBi+vqgsVbV+by/iy0Ap7KKTBKvEOYUuvB/9T43trcnq2cMNf4esfgXSivrXuUeGezGun9EGX8f+cpWLhXXrn2Z8rLp5sOFhIRw8ODBJrU9ceIEY8aMqXHspptuYtOmTQAMGDCAXbt2GZ/78MMPufvuuwHDkHJ+fj6//vprrWsIcaWkh08IM7ohzI9RPT2oqNLxyk/HmnWuk8aant6OxkLLfxfgakj4Mi4P/QohzEOv19c5raKp3nnnHRITE7nvvvsoKfmzlNLx48fZs2cPU6dOBcDKyoopU6bw4YcfXnHMQvyd9PAJYUYKhYIF14USs3w7cUcyOHw+n37+zia5truDGoCLxebfzk2IK5X0Qky9zyn/lkwlLIhuctvfnh17ZYE1wdGjRwkODm5S2549e3L8+PEax3x9fQFwc3OrcfzDDz+kqqoKPz8/4zG9Xo+NjQ1vv/02zs6mea8QAqSHTwiz6+XtyA1hhjf01zedaPJ5H+w4w/JfTtRbXNnd3gaAi0Wm28pNCHOxU1vV+/j7tAVTtDWVLVu2cOjQIW655ZYmtZ86dSqbNm1i//79Dbarqqrik08+YdmyZSQmJhofBw4cwM/Pj88//9wU4QthJD18QrSCWeN78v2BC2w5lsX+c7kM7OLa6Dmf7T7HmZxiorq508W99rDunz18kvAJYQrl5eVkZGSg1WrJzMwkLi6OJUuWcN111zFt2rQmXePJJ5/kxx9/ZPz48cTGxjJq1ChcXV05ceIEP/30EyqVIWH94YcfyM3NZebMmbV68m655RY+/PBDHnzwQZO/RtF5SQ+fEK2gm6cDNw0MAOCD35KbdE5u9Spde3Wdz1f38F0qrrjiHT2EEBAXF4evry9BQUFMnDiRrVu3smLFCr799ltjotYYjUbD5s2befbZZ1m9ejUjR46kT58+PPHEE4wYMYINGzYAhuHc6OjoOodtb7nlFvbu3dvkhSJCNIVC38nWf6elpREYGEhqaioBAQGWDkd0Ikcu5DNpxW9YKRXsmjMOrwaKKev1errP24hOD7vnjce7jrYVVTqe/CIRD3s18yb1wcaqaR9IQgjRHsnn95WRHj4hWklfP2ciurpSpdOz9o/UBtsWV2ip7rRzqqeostpKyco7B/H8jf0k2RNCCNEgSfiEaEX3DOsKGObnVWl19bYrKK0EwFqlQGMt/0yFEEJcGfkkEaIVXdPfB3d7NRkFZfxyNLPedoVlhp0CnDTWDdb/0ur0XCwqp6jcdDsLCCGE6Hgk4ROiFdlYqZgyJBCA//x+tt52BWWGHj4n24b3yH34vwlEvPgL3+w/b7oghRBCdDiS8AnRyu6M7IJCATtPXeR0dt3bovX3d2bTk1ex8s5BDV7L7fIK3ktSi08IIUQDJOETopUFuNoxupcnAD8eTK+zjcZaRU9vR0L9nBq81p+lWWS3DSGEEPWThE8IC7imnw8A/zuScUXXqa7Rl1tSecUxCSGE6Lgk4RPCAqL7eKNUwJELBaReqr112o6T2byx6QQ7TmY3eB1XO8Mcv+oizUIIIURdJOETwgLcHWwYEmTYSL2uXr4dJ3N4c/NJtp9oJOGrnsMn26sJIYRogCR8QljIxMvDuj8fqV2epboOn2M9RZerudkZEr48GdIVQgjRAEn4hLCQCX0NCd8fZy+RU1Rz0cWfdfisGryGj7OG68P8uLa/j3mCFEII0SFIwieEhfi72NLf3xm9Hn5JqtnLV12Hr7EePm8nDW9NHchzk0LNFqcQQoj2TxI+ISwouo83YJiz91cF1T18jRReFkIIIZpCEj4hLGhkT3cAdp3OQafTG48XGnv4Gh7SBajS6rhYVE5ZpdY8QQohhGj3JOETwoIGBLhgr1aRW1JJUnqB8XhBqaGHrykJ36QVvxHx4i/sTck1W5xCCNFerVy5kqCgIDQaDZGRkezZs6fB9uvXryckJASNRkP//v3ZuHFjjee//vprJkyYgLu7OwqFgsTExHqvpdfrueaaa1AoFGzYsMEEr6blJOETwoKsVUqGdfuzl6/a2gci+eqhKLp5ODR6DZfLtfguym4bQghRw7p165g9ezaxsbHs27ePsLAwYmJiyMrKqrP9rl27mDp1KjNnzmT//v1MnjyZyZMnc/jwYWOb4uJiRo4cySuvvNLo/ZcvX45CoTDZ67kSkvAJYWHDe3gA8Nupi8ZjPbwciejqhq1a1ej5no6G7dWyCyXhE0J0fIWFhRQUFBgf5eX1v/e9/vrr3H///cyYMYPQ0FBWrVqFnZ0dH330UZ3t33zzTSZOnMjTTz9Nnz59WLx4MYMGDeLtt982trnnnntYuHAh0dHRDcaZmJjIsmXL6r1Xa5OETwgLG9HD0MP3R/IlKqp0zT7f20kDQJYkfEKITiA0NBRnZ2fjY8mSJXW2q6ioICEhoUZiplQqiY6OJj4+vs5z4uPjayVyMTEx9bavT0lJCXfeeScrV67Ex6dtlM1qfIKQEMKsens74uGgJqeogv3ncunu5cCnv5/Fw8GGu4d1bfR8bydDD19mQZm5QxVCCItLSkrC39/f+LuNjU2d7XJyctBqtXh7e9c47u3tzbFjx+o8JyMjo872GRnN2/f8ySefZPjw4dx4443NOs+cJOETwsIUCgXDu3vw3YEL7DyVg8ZaxfJfTuLnrGlSwufleLmHr0B6+IQQHZ+joyNOTk6WDqNe3333HVu2bGH//v2WDqUGGdIVog2oHtbdefqicZeNxoouV/Oq7uErlB4+IYSo5uHhgUqlIjOzZmH7zMzMeodZfXx8mtW+Llu2bOH06dO4uLhgZWWFlZWhb+2WW25hzJgxzXsRJiQJnxBtwPDuhoUbB1LzyLqcuDnZNq0DvoubHZMG+HJNv7YxT0QIIdoCtVpNREQEmzdvNh7T6XRs3ryZqKioOs+Jioqq0R5g06ZN9bavy5w5czh48CCJiYnGB8Abb7zB6tWrm/9CTESGdIVoAwLd7OjqbsfZiyUcTMsDmt7DF+Bqx8o7B5kxOiGEaJ9mz57N9OnTGTx4MEOHDmX58uUUFxczY8YMAKZNm4a/v79x4cesWbMYPXo0y5YtY9KkSaxdu5a9e/fy3nvvGa956dIlzp07x4ULFwA4fvw4YOgd/Ovj77p06UJwcLC5X3K9pIdPiDYiMtgNgBOZRUDTii4LIYSo35QpU3jttddYuHAh4eHhJCYmEhcXZ1yYce7cOdLT043thw8fzmeffcZ7771HWFgYX375JRs2bKBfv37GNt999x0DBw5k0qRJANxxxx0MHDiQVatWte6LayaFXq/XN96s40hLSyMwMJDU1FQCAgIsHY4QRmv3nGPO14fwd9FwPq+Me4Z1ZfHkfo2fCOh0enKKy7FTW+FgI4miEKLjkc/vKyM9fEK0ERFdXQHIuLzatjk9fPeu+YOhL20m7nDzSgcIIYToHKQrQIg2orunA04aKwrKqnj1lgHGLdeawttRavEJIYSon/TwCdFGKJUKBnYx9PKVVmrp4m7X5HN9nQ21+NLzS80SmxBCiPZNEj4h2pCBXVwASEzNa9Z5Ps62AGTkSw+fEEKI2iThE6INCQ90AWDb8axmJW9/9vBJwieEEKI2SfiEaEOqE77ckkqOZxQ0+Tyfywmf9PAJIYSoiyR8QrQhLnZqlArDz+fzmj4fr7qH72JxBWWVWnOEJoQQoh2TVbpCtCF6vR7d5cqYqblNT/icba25MdwPTwcbKrU6NNYqM0UohBCiPZKET4g2pLjiz96501lFTT5PoVDw5h0DzRGSEEKIDkCGdIVoQwrLKo0/J6XnWzASIYQQHYkkfEK0IQWlVcaf03LLyC+pbKB1TTqdnqyCMrIKZeGGEEKImiThE6INqe7hs7q8cuPIhab38r2+6QRD/7WZtzafMktsQggh2i9J+IRoQ3r7OPLlg1EMDjLsuHHofNMTPh/ZbUMIIUQ9JOETog1x1FgzOMiNUT09ATh8oem1+PxcDAnfhTwZ0hVCCFGTxRO+lStXEhQUhEajITIykj179jTYfvny5fTu3RtbW1sCAwN58sknKSuTDzjRsfTzdwbgSDN6+PxcDNurNad+nxBCiM7BognfunXrmD17NrGxsezbt4+wsDBiYmLIysqqs/1nn33GnDlziI2N5ejRo3z44YesW7eOefPmtXLkQpjHHymX+Oi3ZCq1OgDO5BTXWLnbkC5udgDkl1aSV1JhthiFEEK0PxZN+F5//XXuv/9+ZsyYQWhoKKtWrcLOzo6PPvqozva7du1ixIgR3HnnnQQFBTFhwgSmTp3aYK9geXk5BQUFxkdhYaG5Xo4QV+yXo5m88EMSv5++aNw942h60/7O2qmt8HayASDlYonZYhRCCNH+WCzhq6ioICEhgejo6D+DUSqJjo4mPj6+znOGDx9OQkKCMcE7c+YMGzdu5Nprr633PkuWLMHZ2dn4CA0NNe0LEcKECssMZVkcNdb09TMM6x5uxrBuV3d7AM5eLDZ9cEIIIdotiyV8OTk5aLVavL29axz39vYmIyOjznPuvPNOXnjhBUaOHIm1tTXdu3dnzJgxDQ7pzp07l/z8fOMjKSnJpK9DCFMqKDUM3zrZWtHP3wmAw80ozXLdAF/+b2QwwR72ZolPCCFE+9Sutlbbtm0b//rXv3jnnXeIjIzk1KlTzJo1i8WLF7NgwYI6z7GxscHGxsb4e0FB01c9CtHa/trDF+hqmJN35HzT/85OiwoyR1hCCCHaOYslfB4eHqhUKjIzM2scz8zMxMfHp85zFixYwD333MP//d//AdC/f3+Ki4t54IEHeO6551AqLb7oWIgrUnB5gYajxsq4UvdkViGlFVps1SpLhiaEEKIds1iGpFariYiIYPPmzcZjOp2OzZs3ExUVVec5JSUltZI6lcrwIajX680XrBCtpKRcC4CDjWEBhoeDGp0ejmU0rZevSqsjI7+MM9lF5gxTCCFEO2PRLrHZs2fz/vvv8/HHH3P06FEeeughiouLmTFjBgDTpk1j7ty5xvbXX389//73v1m7di3Jycls2rSJBQsWcP311xsTPyHas5JKw5CurVqFQqEw9vI1tQDz72cuMWzJZh76dJ/ZYhRCCNH+WHQO35QpU8jOzmbhwoVkZGQQHh5OXFyccSHHuXPnavTozZ8/H4VCwfz58zl//jyenp5cf/31vPTSS5Z6CUKY1Jt3DCS/tJLung4A9PNzZtvx7CYXYHaxswYgV+rwCSGE+AuFvpONhaalpREYGEhqaioBAQGWDkeIBsUdTufBT/fRz9+JHx4b1Wj783mljHh5C2qVkuMvTkShULRClEIIYX7y+X1lZJWDEG1YdS2+4xmFVFTpGm3vYmvo4avQ6iip0Jo1NiGEEO2HJHxCtBGVWh0f70rhiz9Sqbq8tVqAqy3OttZUavWcyGx8xw07tQq1yvDPWoZ1hRBCVJOET4g2oqisitjvjvDMVweNxwwLNwwFmI80oQCzQqHA+fI8vvzSpu3BK4QQouOThE+INqK00jAEq1YpsVL9+U+zn3GLtaat1HW+PKybXyIJnxBCCIN2tdOGEB1Z9Zw7jXXN72F9L5dmOdTElbo3DfQnr6QCLyeNaQMUQgjRbknCJ0QbUXa5h89OXfOfZaivYUj3WEYBWp0elbLhlbePjO1hngCFEEK0WzKkK0QbUd3D9/ct1II97LG1VlFWqSM5R3bQEEII0XyS8AnRRlTP4bO1rpnwqZQK+vg6AnCkCTtulFVquZBXSlZhmemDFEII0S5JwidEG1Fa8ee2an8X6mcY1k1qQsL37q9nGP7yFt7YdNK0AQohRDu0cuVKgoKC0Gg0REZGsmfPngbbr1+/npCQEDQaDf3792fjxo01nv/666+ZMGEC7u7uKBQKEhMTazx/6dIlHnvsMXr37o2trS1dunTh8ccfJz+/afOwzUUSPiHaiEFdXfno3sE8HdO71nPVBZib0sPnbGuYA1ggZVmEEJ3cunXrmD17NrGxsezbt4+wsDBiYmLIysqqs/2uXbuYOnUqM2fOZP/+/UyePJnJkydz+PBhY5vi4mJGjhzJK6+8Uuc1Lly4wIULF3jttdc4fPgwa9asIS4ujpkzZ5rlNTaVbK0mRDtwMC2PG97eiZu9moT50Q1umfbN/jSeXHeAkT08+PT/IlsxSiGEMJ/qz++kpCT8/f2Nx21sbLCxsanznMjISIYMGcLbb78NgE6nIzAwkMcee4w5c+bUaj9lyhSKi4v54YcfjMeGDRtGeHg4q1atqtE2JSWF4OBg9u/fT3h4eIOxr1+/nrvvvpvi4mKsrCyzXlZ6+IRoB3p5O6JSKrhUXEFGQcNz81xs1YAUXhZCdEyhoaE4OzsbH0uWLKmzXUVFBQkJCURHRxuPKZVKoqOjiY+Pr/Oc+Pj4Gu0BYmJi6m3fVPn5+Tg5OVks2QMpyyJEm3EoLZ/jmYX08nZgQIBLjec01ip6eDpwPLOQI+cL8HW2rfc6TpcLL+eVytZqQoiOp64evrrk5OSg1Wrx9vaucdzb25tjx47VeU5GRkad7TMyMlocb05ODosXL+aBBx5o8TVMQXr4hGgjfjqczj/XH+Cb/efrfN64cCO94Xl8stOGEKIjc3R0xMnJyfioL+FrCwoKCpg0aRKhoaEsWrTIorFIwidEG2Gsw2dde5UuQF+/pu2p63J5L93C8iq0uk41RVcIIYw8PDxQqVRkZmbWOJ6ZmYmPj0+d5/j4+DSrfUMKCwuZOHEijo6OfPPNN1hbWzf7GqYkCZ8QbcSfO23UnfCFGhO+hnv4XO3U3BDmx92RXXk17hi7TuWYNlAhhGgH1Go1ERERbN682XhMp9OxefNmoqKi6jwnKiqqRnuATZs21du+PgUFBUyYMAG1Ws13332HRmP5rS4l4ROijfhzL916Er7LW6yl5ZY2uCBDpVSwYupAtHo9724/wz8+TaCiSmf6gIUQoo2bPXs277//Ph9//DFHjx7loYceori4mBkzZgAwbdo05s6da2w/a9Ys4uLiWLZsGceOHWPRokXs3buXRx991Njm0qVLJCYmkpSUBMDx48dJTEw0zvOrTvaKi4v58MMPKSgoICMjg4yMDLRabSu++ppk0YYQbURpPXvpVnOxU+PvYsv5vFKSLhQQ1d29wet5OBjmtRSWVXH2YjE9vR1NG7AQQrRxU6ZMITs7m4ULF5KRkUF4eDhxcXHGhRnnzp1Dqfyz72v48OF89tlnzJ8/n3nz5tGzZ082bNhAv379jG2+++47Y8IIcMcddwAQGxvLokWL2LdvH7t37wagR4+ae5snJycTFBRkrpfbIKnDJ0QbcfcHu/ntVA7Lp4QzeaB/nW0e+GQvPydlMn9SH/5vVLdGr3ndWzs4fL6Ad++JIKZv8+egCCFEWyGf31dGhnSFaCOqe/jqG9IF6O9v2HHj8PmmbdHj42Qo35JTVH6F0QkhhGjPZEhXiDbi6ZjeZBWWExboXG+bfpcTvkNNTPjc7Q1FmC8VSU0+IYTozCThE6KNGNat4Tl58GfCdyanmOLyKuxtGv4n7O5gSPguFkvCJ4QQnZkkfEK0I56ONng72ZBZUE5SegFDgtwabH9VL0/sbawYGOjSOgEKIYRok1qU8J07d46zZ89SUlKCp6cnffv2bdOVroVoD747cAFrpYLRvT3rXakLhnl8mQVZHD6f32jCN6ybe5N6DoUQQnRsTU74UlJS+Pe//83atWtJS0vjr4t71Wo1o0aN4oEHHuCWW26pscRZCNE0T32RSKVWT/zccQ0mfH39nPnlaFaT5/EJIYQQTcrMHn/8ccLCwkhOTubFF18kKSmJ/Px8KioqyMjIYOPGjYwcOZKFCxcyYMAA/vjjD3PHLUSHUqnVUak1fImqb2u1atUrdY+cb3jHDYCKKh1JFwrYfebilQcphBCi3WpSD5+9vT1nzpzB3b320JCXlxfjxo1j3LhxxMbGEhcXR2pqKkOGDDF5sEJ0VNXbqgHY1rO1WrXqhRsnswoprdA22D67qJxrV+zAWqXgxIvXoFAoTBOwEEKIdqVJCd+SJUuafMGJEye2OBghOqvSy9uqKRWgVjXc8e7tZIOHgw05ReUczShgUBfXetu62RlW6VZq9RSVV+Gosezm3UIIISxDJtsJ0QYUlVcBYG9j1WgvnEKhoJ+/YV/dxgow26pVxiHi3OL6998VQgjRsZks4Zs3bx733XefqS4nRKdSWGZI+BwbqatXrXoe36G0xhduuFUXXy6RWnxCCNFZmSzhO3/+PCkpKaa6nBCdSnUPn4OmeQnfwSYkfK72hmHcXCm+LIQQnZbJCi9//PHHprqUEJ1OT28H3rwjHBurhhdsVAu7XEj5ZFYhJRVVDZZxcb08j++SJHxCCNFpyRw+IdoAL0cNN4b7M7GfT5Paeztp8HHSoNPD4UbKs1Tvp3uxuPyK4xRCCNE+NbuH74UXXmjw+YULF7Y4GCFE0w0IcCYjqYwDqXkMDa5/x42J/Xzo6e3I8O4erRidEEKItqTZCd8333xT4/fKykqSk5OxsrKie/fukvAJ0QJH0ws4d6mE7p4O9PByaNI5YYEu/JyUyYG0vAbbTezna4IIhRBCtGfNTvj2799f61hBQQH33nsvN910k0mCEqKz2bD/PO9uP8P9o4J5blJok84ZEND0hRtCCCE6N5PM4XNycuL5559nwYIFpricEJ1OweWyLA42TS+MPMDfBYBzl0oaXZBRUFbJ2j3nOJlZ2OIYhRBCtF8mW7SRn59Pfr70NAjREtVlWRybWJYFwNnOmu6e9gDsO5vbYNsFGw4z5+tDvL7pRMuDFEII0W41e0h3xYoVNX7X6/Wkp6fzn//8h2uuucZkgQnRmRSWGXbBcGhi4eVqg7u6cTq7mL1nc4kO9a633bSoIL5NvMDPSZnkl1bibCtbrAkhRGfS7ITvjTfeqPG7UqnE09OT6dOnM3fuXJMFJkRnkldiSPic7ZqXiEUEubJubyoJZy813K6rKz28HDiVVcTOUzlc218WcgghRGfS7IQvOTnZHHEI0anlXd72rLpIclMN7uoKwIG0fMqrtA0Wbh7Zw4NTWUXEn74oCZ8QQnQyUnhZiDYg93IPn5t983r4gj3scbdXU1Gla7QA8/Du7gDsOp3TsiCFEEK0WyZL+N55551GizILIer2r5v6s/jGvvg42zbrPIVCwaDLvXyNDetGBhsSvpSLJSTnFLcsUCGEEO2SyRK+r776ijVr1pjqckJ0KpMG+HJPVFCzF23An8O6e1MaXqnrbGfNmN6eaHV6fj2e1aI4hRBCtE/N/3Spx+bNm011KSFEMwwOqu7hy0Wv16NQKOpt+949gzmTU0Rvb0cqqnRYqxQNthdCCNExyBw+ISwsq7CMuMPpJKbmtej8fv7OqK2UXCyuIOViSYNt1VZKenk5cv3bv9E3No7MgvIW3VMIIUT70uIevqSkJM6dO0dFRc0K/zfccMMVByVEZ3IwNZ8HP91HWIAz3z46stnn21ipGODvzN6zuexNuUSwh32D7ZVKBeWVOiq1eo6mF+DjrGlp6EIIIdqJZid8Z86c4aabbuLQoUMoFAr0ej2AcVhIq9WaNkIhOri8UsMKXacrKIYcEeTK3rO5JJzN5bbBgY22D/F14mRWEUczChgb4tXi+wohhGgfmj2kO2vWLIKDg8nKysLOzo4jR46wfft2Bg8ezLZt28wQohAdW8HlhO9Kdr8Y3NUNgL2NbLFWrY+vIwBJFxou5SKEEKJjaHYPX3x8PFu2bMHDwwOlUolSqWTkyJEsWbKExx9/nP3795sjTiE6rHwTJHwRl1fqnsoqIq+kApdGCjiHB7oAsONkDgVllThpZKs1IYRoC/Ly8vjmm2/YsWMHZ8+epaSkBE9PTwYOHEhMTAzDhw9v0XWb3cOn1WpxdDT0Dnh4eHDhwgUAunbtyvHjx1sUhBCdWUHZlQ/putmr6eZpmLuX0IRevqFBbvg6a8gvreSpLw60+L5CCCFM48KFC/zf//0fvr6+vPjii5SWlhIeHs748eMJCAhg69atXH311YSGhrJu3bpmX7/ZCV+/fv04cMDwAREZGcmrr77Kzp07eeGFF+jWrVuzAxCiszNFDx/8pR5fExI+K5WSVXdHoFDApqRMzjWyulcIIdqrlStXEhQUhEajITIykj179jTYfv369YSEhKDRaOjfvz8bN26s8fzXX3/NhAkTcHd3R6FQkJiYWOsaZWVlPPLII7i7u+Pg4MAtt9xCZmZmg/cdOHAgbm5uJCQkcPr0aT7//HOWLVvGiy++yDvvvMP27dvJyclhwYIFLF++nNdee61Zfw7NTvjmz5+PTqcD4IUXXiA5OZlRo0axceNGVqxY0dzLCdHpFZRWAVzxsGr1PL6ERgowVwsLdOH128P47P5IPBybt4evEEK0B+vWrWP27NnExsayb98+wsLCiImJISur7uLzu3btYurUqcycOZP9+/czefJkJk+ezOHDh41tiouLGTlyJK+88kq9933yySf5/vvvWb9+Pb/++isXLlzg5ptvbjDWpKQkXn31VUJDQ+ttY2try9SpU4mPj2fGjBmNvPqaFPrqZbZX4NKlS7i6uraLAq5paWkEBgaSmppKQECApcMRgp2ncjiTU0xUNzd6eDm2+Dqns4sYv+xXbKyUHFoUg9pKymwKITqOlnx+R0ZGMmTIEN5++20AdDodgYGBPPbYY8yZM6dW+ylTplBcXMwPP/xgPDZs2DDCw8NZtWpVjbYpKSkEBwezf/9+wsPDjcfz8/Px9PTks88+49ZbbwXg2LFj9OnTh/j4eIYNG9bcl24SJvlEcHNzaxfJnhBt0YgeHtwzrOsVJXsA3TzscbNXU16l4/CFfBNFJ4QQbUthYSEFBQXGR3l53QXkKyoqSEhIIDo62nhMqVQSHR1NfHx8nefEx8fXaA8QExNTb/u6JCQkUFlZWeM6ISEhdOnSpdHrJCQkMHbsWAoKaldQyM/PZ+zYscZpdc3VpITvwQcfJC0trUkXXLduHf/9739bFIwQouUUCoVxte6e5EtNOken0/Nt4nle/CGJ0gqpoSmEaPtCQ0NxdnY2PpYsWVJnu5ycHLRaLd7e3jWOe3t7k5GRUec5GRkZzWpf3zXUajUuLi7Nvs6yZcsYN24cTk5OtZ5zdnbm6quvZunSpU2O5a+aVJbF09OTvn37MmLECK6//noGDx6Mn58fGo2G3NxckpKS+O2331i7di1+fn689957LQpGiM5oU1ImDjZWDOrqgo2V6oquNby7O5uSMtl1+iIPju7eaHuFAhb/cJSconKuD/Mj7HK5FiGEaKuSkpLw9/c3/m5jY2PBaExr9+7ddQ41V7v++uv54IMPWnTtJiV8ixcv5tFHH+WDDz7gnXfeISkpqcbzjo6OREdH89577zFx4sQWBSJEZ1Sp1XH/J3sB2LfgahMkfB4A/JF8iYoqXaPz+BQKBV3d7cgpKud8XqkkfEKINs/R0bHOHrC/8/DwQKVS1Vodm5mZiY+PT53n+Pj4NKt9fdeoqKggLy+vRi9fU65z/vx5Y+m7ujg4OJCent7kWP6qyXP4vL29ee655zh06BA5OTns27ePnTt3cvz4cXJzc/nyyy8l2ROimYrLq4w/O2pavLW1US9vB9zt1ZRWajmQltekc6r30k3PL7vi+wshRFuhVquJiIhg8+bNxmM6nY7NmzcTFRVV5zlRUVE12gNs2rSp3vZ1iYiIwNrausZ1jh8/zrlz5xq9jqenZ4M1jY8dO4aHh0eTY/mrFi3acHV1JSwsjGHDhtGjR48rWrDR3Po4eXl5PPLII/j6+mJjY0OvXr1q1cgRor0ovjxvTm2lxFp15WuoFAoFw7q7A4bVv03hV53w5ZVe8f2FEKItmT17Nu+//z4ff/wxR48e5aGHHqK4uNhY0mTatGnMnTvX2H7WrFnExcWxbNkyjh07xqJFi9i7dy+PPvqosc2lS5dITEw0jnYeP36cxMRE4/w8Z2dnZs6cyezZs9m6dSsJCQnMmDGDqKioRlfoRkdH89JLL9X5nF6v56WXXqq1qKSprrxL4QpU18dZtWoVkZGRLF++nJiYGI4fP46XV+0N3SsqKrj66qvx8vLiyy+/xN/fn7Nnz9aaGClEe1FyuYfPXn1lQ7l/NaK7Bz8eTGfXqYs80YT3BR9nW0B6+IQQHc+UKVPIzs5m4cKFZGRkEB4eTlxcnHFhxrlz51Aq//yyPXz4cD777DPmz5/PvHnz6NmzJxs2bKBfv37GNt99912NGnh33HEHALGxsSxatAiAN954A6VSyS233EJ5eTkxMTG88847jcY7f/58IiIiiIyM5KmnnqJ3796AoWdv2bJlnDhxgjVr1rToz8Ikdfhaqrn1cVatWsXSpUs5duwY1tYtK1IrdfhEW5KYmsfklTvxd7Fl55xxJrlmSk4xY17bhrVKwYHYCdipG/5e99OhdB767z4GdXHh64dHmCQGIYQwtc7y+b13717uvfdekpKSjCOoer2e0NBQVq9ezZAhQ1p0XYv18FXXx/lrV2pj9XG+++47oqKieOSRR/j222/x9PTkzjvv5Nlnn0WlqruHpLy8vEaNnsLCQtO+ECGugLGHz8Z0PXxd3e3wd7HlfF4pf6TkMrqXZ4Ptq+fwZUgPnxBCWNzgwYM5fPgwiYmJnDx5Er1eT69evWoUd24JiyV8DdXHOXbsWJ3nnDlzhi1btnDXXXexceNGTp06xcMPP0xlZSWxsbF1nrNkyRKef/55k8cvhClUz+FrrBeuORQKBcO7u7M+IY1dp3IaTfh6eTvy+f3DjImfEEIIywsPD7/iJO+vmj1LfNy4ceTl5dU6XlBQwLhxphmSqo9Op8PLy4v33nuPiIgIpkyZwnPPPVdru5O/mjt3Lvn5+cbH30vKCGFJIT6OLL6xL/eNDDbpdUf0MKzi2nm68YUb9jZWRHV3J9jDniqtjk1JmeQU1V25XgghhHmtXbuWVatWUVJSYtLrNrtbYdu2bVRUVNQ6XlZWxo4dO5p8nZbUx/H19cXa2rrG8G2fPn3IyMigoqICtbr2BvA2NjY1ijLWtV2JEJYS6GbHPVFBJr/u8MsrdY9cKCCvpAIXu9r/NupyNL2Qxz7fR4iPExsekfl8QgjRmh5//HFOnz5Nt27dmDBhAr/99pvJrt3khO/gwYPGn5OSkmpsD6LVaomLi6tR+boxf62PM3nyZODP+jh/Xf78VyNGjOCzzz5Dp9MZV9WcOHECX1/fOpM9ITorLycNPb0cOJlVRPzpi1zT37fB9hsPpXMso5Cobu4oUJCYmkdZpRaNtenmFgohhGjYV199RVxcHP369cPGxoasrKw6q5a0RJMTvvDwcBQKBQqFos6hW1tbW956661m3Xz27NlMnz6dwYMHM3ToUJYvX16rPo6/v79xn7yHHnqIt99+m1mzZvHYY49x8uRJ/vWvf/H44483675CtBXJOcVk5JcR4GpLoJudSa89oocHJ7OK2Hk6p9GEb+0fqWw/kY2fswYrlQIqIfVSCT2966/4LoQQwrTCwsL48ssvOX78OG5ubi0uslyXJid8ycnJ6PV6unXrxp49e/D0/HMiuFqtxsvLq96VsvVpbn2cwMBA/ve///Hkk08yYMAA/P39mTVrFs8++2yz7itEW7F2zzne3X6G+0cF89ykUJNeO6q7O2t2pbDr1MVG2wa527EdSLlYgp+zLcfLCknPL5OETwghWtHq1atZsGABR48e5ccff6yRA12pJid8Xbt2BQzDrqb06KOP1juEu23btlrHoqKi+P33300agxCWUlxhKMtiylW61YZ1c0epgDM5xZzPK8XfxbbetkHu9gCs+vW08Vh6vuy8IYQQrcnb25v33nvPLNdu0afMyZMn2bp1K1lZWbUSwIULF5okMCE6g5JyQ1kWU9bhq+Zsa83ALq4knM3l1+PZ3BnZpd62wZ72tY7JzhtCCNG6Ro0axeTJk7nhhhvo2bOnSa/d7ITv/fff56GHHsLDwwMfH58a++gqFApJ+IRoBnP28AGM7uVpSPhOZDWY8PXwdDD+HN3Hm1+OZpKeJwmfEEK0pvvvv59vv/2WRYsWERAQwA033MANN9zA8OHDa+RbLdHsT5kXX3yRl156SebNCWECJRXm6+EDGNPbk9c3nWDnqYtUanVYq+qeDxLoZsfzN/RFqVTg56xBoYB+/k5miUkIIUTdpk2bxrRp0ygvL2fz5s18++233HbbbWi1WiZNmsQNN9xATEwMtrb1T9GpT7NnA+bm5nLbbbc1+0ZCiNqKy83bw9fPzxk3ezVF5VUknM1tsO304UHcM6wr4/t48/60wWapDyiEEKJxNjY2XHvttbz77rtcuHCB7777Dl9fXxYsWIC7uzvXXXcdO3fubNY1m53w3Xbbbfz888/NPU0IUYcS49Zq5unhUyoVXNXTsKz/1xPZZrmHEEII84qMjOSll17i0KFDHDp0iPHjx5Oent6sazSpW2HFihXGn3v06MGCBQv4/fff6d+/P9bW1jXaSk08IZrugau6kVFQRre/zKEztTG9vdiQeIFfj2fz7MSQJp2j1+vJLanESWOFVT3DwEIIIVpf9+7defLJJ5t9nkKv1+sbaxQc3LR9PhUKBWfOnGl2EK0pLS2NwMBAUlNTCQgIsHQ4QpjdxaJyBr/0C3o97Jk3Hi8nTaPnRC3ZTHp+GXFPjCLER+byCSEsTz6/4cCBAwwaNAitVtvsc5vUw5ecnNzsCwsh2gZ3Bxv6+ztzMC2fX09kc9vgwEbPcba1Jj2/jPS8Mkn4hBCiDWlCP12dzDNTXAjRKJ1Ozx8pl7C3saKPrxMq5ZUtuW/ImF6ezUr4unnacyyjkOOZhYwNMc0+jkIIIRp28803N/h8fn5+i8uzNDvhmz17dp3HFQoFGo2GHj16cOONN+Lm5taigIToLIorqpjynmHXmGOLJ6JSmmfhBsDo3p6s2HKKHSdzqNLqGp2XFx7owsZDGexrZGWvEEII0/n++++5+uqrjVvM/l1LhnKrNTvh279/P/v27UOr1dK7d28ATpw4gUqlIiQkhHfeeYennnqK3377jdBQ0+4NKkRHUnp5ha5SATZW5l0YERbggpPGivzSSg6k5RHRteEvZGEBLgAkpuah1+uvuOCnEEKIxvXp04dbbrmFmTNn1vl8YmIiP/zwQ4uu3exPmRtvvJHo6GguXLhAQkICCQkJpKWlcfXVVzN16lTOnz/PVVdd1aIVJEJ0JsXVRZfVVmZPqKxUSkb18gRg2/HGy7MMCHDB1lpFVmE5B9PyzRqbEEIIg4iICPbt21fv8zY2NnTpUv+uSQ1pdsK3dOlSFi9ejJPTnxO5nZ2dWbRoEa+++ip2dnYsXLiQhISEFgUkRGdRXXTZ1kw1+P5ubG/DXLymJHy2ahVjehsSxK3Hs8walxBCCINVq1axdOnSep/v06dPixfSNntINz8/n6ysrFrDtdnZ2RQUFADg4uJCRUVFiwISorOoLrrsYNM6a6dGX+7hO3Q+n6zCMrwcGy7PctNAfwLd7IyJohBCCPOysbEx27VbNKR733338c0335CWlkZaWhrffPMNM2fOZPLkyQDs2bOHXr16mTpWITqU4orL26qZaR/dv/N0tGFAgDMA20/kNNp+Ql8f5l3bh7BAFzNHJoQQori42Kztm53wvfvuu4wfP5477riDrl270rVrV+644w7Gjx/PqlWrAAgJCeGDDz5o7qWF6FRKyqu3VWu96khjeskwrRBCtEU9evTg5ZdfbnDLNL1ez6ZNm7jmmmtq7ILWFE3aaaMuRUVFxl01unXrhoOD+baGMiWp1C3ailNZhfzvSCbeThpujWidv4v7zuVy8zu7cNRYsW/B1Vg3Up6lrFLL7uRL5JVUcGO4f6vEKIQQdenon9/Hjx9n3rx5/Pjjj4SFhTF48GD8/PzQaDTk5uaSlJREfHw8VlZWzJ07l3/84x+oVE0fIWpx14KDgwMDBgxo6elCdHo9vBzp4eXYqvcMC3DB3V7NxeIK9qbkEtXdvcH2+87mMv2jPXg4qLlugJ9Zi0MLIURn1rt3b7766ivOnTvH+vXr2bFjB7t27aK0tBQPDw8GDhzI+++/zzXXXNOsRK9akxK+m2++mTVr1uDk5NRoFeivv/662UEIIVqHSqlgTG8vvtqXxpZjmY0mfIOD3HC2tSanqILfz1xkRA+PVopUCCE6py5duvDUU0/x1FNPmfS6TZrD5+zsbKwT5uzs3OBDCNE0abklHLmQT05Reaved3wfw6rbzUcbn8entlJybX9fADbsP2/WuIQQQphPk3r4Vq9eXefPQoiW+/C3ZFbvTOHhMd15ZmJIq913VE8PrFUKzuQUcya7iG6eDc+/vWmgP5/vOcdPhzP41839G533J4QQou1p0Tt3VVUVv/zyC++++y6FhYUAXLhwgaKiIpMGJ0RHVr1K176V6vBVc9RYExlsGMrdcqzxXr7BXV1xsLGiqLyK5JzmlQEQQgjRNjQ74Tt79iz9+/fnxhtv5JFHHiE721C1/5VXXuGf//ynyQMUoqMy1uFrpZ02/mpcSNOHdZVKBb28Db2AxzIKzRqXEEII82h2wjdr1iwGDx5Mbm4utra2xuM33XQTmzdvNmlwQnRkJX/ZS7e1Vc/j+yPlEvmllY227+1jWE18MlMSPiGEsKTDhw+36Lxmf9JULxNWq9U1jgcFBXH+vEzqFqKpqvfSba2dNv6qq7s9PbwcOJVVxK8nsrkhzK/B9o+M7cFDo3vg59LwdmxCCCFMr7CwkM8//5wPPviAhIQEtFpts6/R7B4+nU5X543S0tJwdGzdmmJCtGeW7OEDiO7jDcD/jmQ02jbA1Y4u7nZYyYINIUQ7s3LlSoKCgtBoNERGRrJnz54G269fv56QkBA0Gg39+/dn48aNNZ7X6/UsXLgQX19fbG1tiY6O5uTJkzXanDhxghtvvBEPDw+cnJwYOXIkW7dubXbs27dvZ/r06fj6+vLaa68xbtw4fv/992ZfB1qQ8E2YMIHly5cbf1coFBQVFREbG8u1117boiCE6IwsOYcP4Jp+PgBsPZZFWWXzvy0KIURbt27dOmbPnk1sbCz79u0jLCyMmJgYsrLqnr+8a9cupk6dysyZM9m/fz+TJ09m8uTJNYZRX331VVasWMGqVavYvXs39vb2xMTEUFZWZmxz3XXXUVVVxZYtW0hISCAsLIzrrruOjIzGv2BnZGTw8ssv07NnT2677TacnJwoLy9nw4YNvPzyywwZMqRlfxj6ZkpNTdWHhobq+/Tpo7eystIPGzZM7+7uru/du7c+MzOzuZdrdampqXpAn5qaaulQRCe3+rcz+pd/OqpPyy2xyP11Op1++JLN+q7P/qCPO5zeaNs3Nh3Xz/p8nz6vuKKVIhRCiD+15PN76NCh+kceecT4u1ar1fv5+emXLFlSZ/vbb79dP2nSpBrHIiMj9f/4xz/0er3hvdDHx0e/dOlS4/N5eXl6Gxsb/eeff67X6/X67OxsPaDfvn27sU1BQYEe0G/atKnBeK+77jq9k5OTfurUqfoffvhBX1VVpdfr9XorKyv9kSNHmvy669LsHr6AgAAOHDjAvHnzePLJJxk4cCAvv/wy+/fvx8vLq2VZpxCd0L0jgnl2Ygj+LraNNzYDhULBxMu9fHGHG/7WqVAo+Gz3OTYkXiD5opRmEUJYTmFhIQUFBcZHeXndxesrKipISEggOjraeEypVBIdHU18fHyd58THx9doDxATE2Nsn5ycTEZGRo02zs7OREZGGtu4u7vTu3dvPvnkE4qLi6mqquLdd9/Fy8uLiIiIBl/bTz/9xMyZM3n++eeZNGlSi7ZQq0+TJw+NHj2a8ePHM2bMGKKiorj77rtNFoQQwjKu7e/Dh78l80tSJuVVWmys6n9zCfKwJ6uwnJScYsIDXVovSCGE+IvQ0NAav8fGxrJo0aJa7XJyctBqtXh7e9c47u3tzbFjx+q8dkZGRp3tq4diq//bUBuFQsEvv/zC5MmTcXR0RKlU4uXlRVxcHK6urg2+tt9++40PP/yQiIgI+vTpwz333MMdd9zR4DlN1eQevuDgYFavXs2YMWNwcXEhOjqal156id9//71Fq0WE6My0Oj1JFwo4e7EYnU5vsTgGBrri7WRDYXkVO0/lNNi2m4c9gBRfFkJYVFJSEvn5+cbH3LlzLR1SDXq9nkceeQQvLy927NjBnj17mDx5Mtdffz3p6ekNnjts2DDef/990tPT+cc//sHatWvx8/NDp9OxadMm42YXLdHkhG/NmjUkJydz5swZ3nrrLfz9/XnvvfcYPnw4rq6uXHPNNSxdurTFgQjRmRSUVnLtih2MXroNnd5yCZ9SqWBiX8Ow7k+HGh7WDb6c8J3Klh11hBCW4+joiJOTk/FhY2NTZzsPDw9UKhWZmZk1jmdmZuLj41PnOT4+Pg22r/5vQ222bNnCDz/8wNq1axkxYgSDBg3inXfewdbWlo8//rhJr9He3p777ruP3377jUOHDvHUU0/x8ssv4+XlxQ033NCka/xds+fwBQUFcd999/Hxxx9z9uxZTp06xeOPP86uXbuYM2dOi4IQorOpXqFrY6W0eKmTif18Afg5KZNKra7edn18nQBIulDQKnEJIcSVUKvVRERE1NgUQqfTsXnzZqKiouo8JyoqqtYmEps2bTK2Dw4OxsfHp0abgoICdu/ebWxTUlICGOYL/pVSqUSnq/89tj69e/fm1VdfJS0tjc8//7zZ51drUQGws2fPsm3bNuMjKyuLYcOGMXr06BYHIkRnUl2Dz1IlWf5qaLAb7vZqLhZX8PuZi4zq6Vlnu75+hoQvOaeYovIqHFp5D2AhhGiu2bNnM336dAYPHszQoUNZvnw5xcXFzJgxA4Bp06bh7+/PkiVLAMNuYqNHj2bZsmVMmjSJtWvXsnfvXt577z3AMD/viSee4MUXX6Rnz54EBwezYMEC/Pz8mDx5MmBIGl1dXZk+fToLFy7E1taW999/n+TkZCZNmtTi16JSqYxlYlqiye/Yn3zyiTHBy8nJYfjw4YwePZr777+fIUOGYG1t3aIAhOiMjLtsWKjo8l+plAom9PXh8z3n2Hgoo96Ez93BBh8nDZeKK0jJKaafv3MrRyqEEM0zZcoUsrOzWbhwIRkZGYSHhxMXF2dcdHHu3LkaPXHDhw/ns88+Y/78+cybN4+ePXuyYcMG+vXrZ2zzzDPPUFxczAMPPEBeXh4jR44kLi4OjcawE5GHhwdxcXE899xzjBs3jsrKSvr27cu3335LWFhYg/Gmp6fz9ttv89JLLwEwcuRIY48hgJWVFd988w3+/v7N/rNQ6PVNm0CkVCrp0qULc+bMYebMme02wUtLSyMwMJDU1FQCAgIsHY7opHaeyuGuD3bTy9uBn5+0fM/4jpPZ3PPhHtzt1ex5LhqVUlFnu9RLJfg4a7CWHTeEEK2sM3x+L1iwgIsXL/LOO+8AhvmK9913H25uboChbMvIkSN57bXXmn3tJr9rv/POOwwbNoznn38eLy8vrr/+epYtW8bevXtpYs4ohLisLfXwAQzr5o6LnTUXiyvYk3yp3naBbnaS7AkhhJn88MMPTJ06tcaxWbNmERsbS2xsLM8//zw//fRTi67d5HfuBx98kLVr15Kens7OnTu59tpr2bNnD5MmTcLV1ZVJkya1KOMUojMqupzwtZV5cNYqJVdf3lv3p8MNlw0QQghhHikpKQQHBxt/v/rqq7G3tzf+3rt3b5KTk1t07RZ9VQ8NDeWhhx5i3bp17N+/n0cffZTffvuNZ599tkVBCNHZ9PRy5MHR3bmmf92lASzh2v6G1bpxhzPqrQ1YVqnl6fUHuP6t32T/XSGEMLHKykqys7ONv3/99dc1ijzn5ubWWv3bVM3uXsjKymLr1q3GBRwnTpzA2tqaYcOGMXbs2BYFIURn0z/Amf4BbWvRw/Ae7jjaWJFVWM6+c7kMDnKr1cbGSsnW49nkFJVz+Hx+nW2EEEK0TO/evdm1axcDBw6s8/kdO3bQq1evFl27yWniww8/TGhoKL6+vkybNo3Dhw9z6623smnTJvLy8ti2bRuxsbEtCkIIYXk2ViqiQ6uHdesuwqxQKIjo6gJAwtnc1gpNCCE6hTvuuIOFCxdy8ODBWs8dOHCAF154odYcv6Zqcg/f/v37mTx5MmPHjmXEiBHY2dm16IZCCDh3sQSFAjwdbdBYW74WX7WJ/Xz4Zv954g5nMH9SHxSK2qt1I7q68r8jmZLwCSGEiT3xxBP88MMPREREcPXVV9O7d28Ajh8/biwA/cQTT7To2k1O+OLj41t0AyFEbc9+dZD4Mxd5845wbgxvfj0lcxndyxM7tYrzeaUcTMsnLNClVpsBAYZjSemy44YQQpiStbU1mzZt4vXXX2ft2rVs27YNgJ49e7J48WKefPLJFpfFk/oKQlhAbkkFAK52agtHUpPGWsXYEC8ANtazWjfExxGAtNxSCsoqWy02IYTo6LZv345SqWTOnDkkJiZSUlJCSUkJBw4cYM6cOfXuG9wUkvAJYQF5JYZEqa0lfADXXt5b96dDGXXW2HSxU+PrbKgofyKjsFVjE0KIjmzs2LFculR/LdQrIQmfEK1Mr9dz6XIPn4td29uxZkxvTzTWSs5dKql32DbEx5EgdzvjnsBCCCGunDk3smgbVV+F6ETySiqpqNIB4OXU8u55c7G3sWJ0L0/+dySTnw5l0NevdvmYj+4dUueCDiGEEFfGXO+tkvAJ0coyCsoAcLNXY2PVdlbo/tW1/X3535FMNh5O56kJvWq9AUmyJ4QQ5nHvvfc2Olfv66+/bvZ1TZbwTZ8+ndTUVLZs2WKqSwrRIVUnfN5OGgtHUr9xIV6oVUrOZBeTlF5QZy+fEEII03N0dMTW1tbk1zVZwufv79/i7T6E6EwCXGx5eEx33Ozb3oKNao4aa64O9ebHQ+ms35tG3xtqJnwpOcXMWrufSq2ejbNGWShKIYToeFasWIGXl5fJr2uyhO9f//qXqS4lRIfW09uRZyaGWDqMRt06OIAfD6WzIfE8c68NqTH8bG9jxYG0fBQKKK/SttmhaSGEaE/MOV1GuuSEEHW6qqcnPk4a8koq2Xw0q8ZzHg5qbK1V6PVwIa/MQhEKIUTH0qZW6er1er788ku2bt1KVlYWOp2uxvMtmUgoRGdyMrMQjbUKH2cN1qq2+51LpVRwS4Q/K7ee5ou9qVzb39f4nEKhIMDVlpNZRaTllhDsYW/BSIUQomPYunUrbm5uZrl2sz9tnnjiCe655x6Sk5NxcHDA2dm5xkMI0bAnv0hk1Ktb2XEy29KhNOrWiEAAtp/IJiO/Zk9egKthUnFabmmrxyWEEB3R6NGjUSqVfPTRR1x33XX069eP/v37c8MNN/DJJ59cUQ9gs3v4/vOf//D1119z7bXXtvimQnRmF4sMRZc9HNpeDb6/C/awZ2iQG3tSLvHVvjQeGdvD+Fygmx0AabkllgpPCCE6FL1ezw033MDGjRsJCwujf//+6PV6jh49yr333svXX3/Nhg0bWnTtZvfwOTs7061btxbdTIjOTq/XGxO+trxK969uGxwAwPq9qTW+XUoPnxBCmNaaNWvYvn07mzdvZv/+/Xz++eesXbuWAwcO8Msvv7BlyxY++eSTFl272QnfokWLeP755yktlTd5IZqrsLyKCq1h3qu7fdvv4QNDEWY7tYqUiyX8kZJrPN7V3Z4gd7t2k7gKIURb9/nnnzNv3jzGjh1b67lx48YxZ84c/vvf/7bo2s0e0r399tv5/PPP8fLyIigoCGvrmnuB7tu3r0WBCNEZVPfu2atV2KrbRykTexsrrhvgyxd701i/N5WhwYYJxTF9fYjp62Ph6IQQouM4ePAgr776ar3PX3PNNaxYsaJF1252wjd9+nQSEhK4++678fb2li2WhGiGi0XlALi3g/l7f3X74EC+2JvGj4fSib2hLw42siujEEKY2qVLl/D29q73eW9vb3Jzc+t9viHNftf+8ccf+d///sfIkSNbdMO6rFy5kqVLl5KRkUFYWBhvvfUWQ4cObfS8tWvXMnXqVG688cYWT2IUojXlXO7hc3doX8OgEV1d6eZhz5mcYjYeTOf2IYHG57Q6PWWVWuwlCRRCiCui1Wqxsqr/vVSlUlFVVdWiazf7HTowMBAnJ6cW3awu69atY/bs2axatYrIyEiWL19OTEwMx48fb3BrkZSUFP75z38yapRs6yTajyAPOx4Z2x0vx7a7j25dFAoFtw4O4NW446xPSDUmfKt3JvPGphPcGdmVOde0/d1DhBCiLdPr9dx7773Y2NQ9ClReXt7iazd70cayZct45plnSElJafFN/+r111/n/vvvZ8aMGYSGhrJq1Srs7Oz46KOP6j1Hq9Vy11138fzzz8uKYdGuhPg48XRMCNOHB1k6lGa7ZVAASgX8kZLLmewiAJw01hSUVfH7mYsWjk4IIdq/6dOn4+XlVavGcfXDy8uLadOmtejaze7hu/vuuykpKaF79+7Y2dnVWrRx6dKlJl+roqKChIQE5s6dazymVCqJjo4mPj6+3vNeeOEFvLy8mDlzJjt27GjwHuXl5TUy4sLCwibHJ4T4k7eThtG9PNl6PJv1CWk8OzGEYd3dATh0Pp+i8iqZ2yeEEFdg9erVZrt2s9+dly9fbrKb5+TkoNVqa01Q9Pb25tixY3We89tvv/Hhhx+SmJjYpHssWbKE559//kpDFcIk0nJL0Or0eDraYKduf8nR7YMD2Xo8m6/3pfHU1b3wd7Gli5sd5y6V8Pvpi0SH1j/ZWAghhOU06xOnsrKSX3/9lQULFhAcHGyumOpVWFjIPffcw/vvv4+Hh0eTzpk7dy6zZ882/n7+/HlCQ0PNFaIQDXr++yQ2JWXy4uR+3D2sq6XDabbxfbxxs1eTWVDOlmNZTOjrw9jennwcf5afkzIk4RNCiDaqWXP4rK2t+eqrr0x2cw8PD1QqFZmZmTWOZ2Zm4uNTu77X6dOnSUlJ4frrr8fKygorKys++eQTvvvuO6ysrDh9+nStc2xsbHBycjI+HB0dTRa/EM1VWFYJgKOm/fXuAaitlNw+2LBg473tZwCMtfh+OZqFVtfyfR6FEEKYT7MXbUyePNlkJVDUajURERFs3rzZeEyn07F582aioqJqtQ8JCeHQoUMkJiYaHzfccANjx44lMTGRwMDAWucI0ZYUlhmW0ztprBtp2XbdNyIItUrJ3rO57E25xJBgN2yslFwqruC8bLMmhBBtUrO7GXr27MkLL7zAzp07iYiIwN7evsbzjz/+eLOuN3v2bKZPn87gwYMZOnQoy5cvp7i4mBkzZgAwbdo0/P39WbJkCRqNhn79+tU438XFBaDWcSHaouqEr7328AF4OWm4eZA/a/9I5d3tZ3h/2mDuiuyKxlqJ2qrZ3yGFEEK0gmZ/6nz44Ye4uLiQkJBAQkJCjecUCkWzE74pU6aQnZ3NwoULycjIIDw8nLi4OONCjnPnzqFUyoeI6Bj+HNJtvz18APdf1Y11e1PZlJTJqaxCFl4v82KFEKIta3YmlZycXO/jzJkzLQri0Ucf5ezZs5SXl7N7924iIyONz23bto01a9bUe+6aNWtklw3RLuj1+g7RwwfQ3dOBq/sYvpRVz+UTQoi2aOXKlQQFBaHRaIiMjGTPnj0Ntl+/fj0hISFoNBr69+/Pxo0bazyv1+tZuHAhvr6+2NraEh0dzcmTJ2td58cffyQyMhJbW1tcXV2ZPHmyKV9Ws11R15ler0evl0naQjRFWaWOqsuLGtp7wgfw4JjuAHyz/zypl0pIyy3h3MUSC0clhBB/qt7NKzY2ln379hEWFkZMTAxZWVl1tt+1axdTp05l5syZ7N+/n8mTJzN58mQOHz5sbPPqq6+yYsUKVq1axe7du7G3tycmJoaysjJjm6+++op77rmHGTNmcODAAXbu3Mmdd95p9tfbEIW+BRnbJ598wtKlS40Zba9evXj66ae55557TB6gqaWlpREYGEhqaioBAQGWDkd0IiUVVby15RSFZZW8cEM/lEqFpUO6YrevimdPyiXG9PJk24lsRvbw4NP/i2z8RCGEaKbqz++kpCT8/f2Nx21sbOrdiiwyMpIhQ4bw9ttvA4aFoYGBgTz22GPMmTOnVvspU6ZQXFzMDz/8YDw2bNgwwsPDWbVqFXq9Hj8/P5566in++c9/ApCfn4+3tzdr1qzhjjvuoKqqiqCgIJ5//nlmzpxpyj+CK9LsHr7XX3+dhx56iGuvvZYvvviCL774gokTJ/Lggw/yxhtvmCNGIToEO7UVz04M4cXJ/TtEsgfwj9GGrQ1/TzZsrZZT1PJ9HoUQoilCQ0NrbDe2ZMmSOttV7+YVHR1tPNbYbl7x8fE12gPExMQY2ycnJ5ORkVGjjbOzM5GRkcY2+/bt4/z58yiVSgYOHIivry/XXHNNjV5CS2j2uNJbb73Fv//97xp7ud1www307duXRYsW8eSTT5o0QCFE2zW2txe9vB04kWnYW1cSPiGEudXVw1eXluzmlZGRUWf7jIwM4/PVx+prU72eYdGiRbz++usEBQWxbNkyxowZw4kTJ3Bzc2vqSzWpZvfwpaenM3z48FrHhw8fTnp6ukmCEqIjKiqv4uzFYvJLKi0diskolQoeGdvD+Pul4gopviyEMCtHR8caGyrUl/BZik6nA+C5557jlltuISIigtWrV6NQKFi/fr3F4mp2wtejRw+++OKLWsfXrVtHz549TRKUEB3Rr8ezGb10G/d/stfSoZjU9QP8CPU17GCj00NuSYWFIxJCiObv5gXg4+PTYPvq/zbUxtfXF6DGNq42NjZ069aNc+fOXcErujLNTvief/55Fi5cyMSJE1m8eDGLFy9m4sSJPP/887zwwgvmiFGIDqG9b6tWH6VSwXOT/nxjO5iaZ7lghBDisubu5gUQFRVVoz3Apk2bjO2Dg4Px8fGp0aagoIDdu3cb20RERGBjY8Px48eNbSorK0lJSaFrV8vtod7sT55bbrmF3bt388Ybbxjr3/Xp04c9e/YwcOBAU8cnRIfRUWrw1WVEDw/s1CpKKrSs3HqKcX28Gz9JCCHMrDm7eQHMmjWL0aNHs2zZMiZNmsTatWvZu3cv7733HmDYYOKJJ57gxRdfpGfPngQHB7NgwQL8/PyMdfacnJx48MEHiY2NJTAwkK5du7J06VIAbrvtttb/Q7isRZ88ERERfPrpp6aORYgOraPsslGfqUO78NFvySScy+PXE9mM7uVp6ZCEEJ1cc3fzGj58OJ999hnz589n3rx59OzZkw0bNtTYvvWZZ56huLiYBx54gLy8PEaOHElcXBwajcbYZunSpVhZWXHPPfdQWlpKZGQkW7ZswdXVtfVe/N+0qA5feyZ1+ISlLPruCGt2pfDwmO48MzHE0uGYxeIfkvjwt2QC3Wz53xNXYafueL2ZQgjLkM/vK9PkOXxKpRKVStXgw8pK3tyFqE9BB+/hA3jy6l74u9iSeqmUpf873vgJQgghWkWTM7Rvvvmm3ufi4+NZsWKFcSmyEKK2glJDwudq1zETvkqtjtziCh4d14O5Xx9iza4Urunny9Bgy9ScEkII8acmJ3w33nhjrWPHjx9nzpw5fP/999x1112ySleIBowL8cbbSUNvH0dLh2IWW49l8cB/EggLcOa2iADWJ6TxxNr9/Pj4KFzt1ZYOTwghOrVml2UBuHDhAvfffz/9+/enqqqKxMREPv74Y4suNxairbszsgsv3dSfgV0sN2nXnHycDROWMwrKiL2hL8Ee9lzIL+ORz/ZRqZXefyGEsKRmJXz5+fk8++yz9OjRgyNHjrB582a+//77GqtXhBCdk4+TIeHLLixHY6Vk1d0R2KlV7Dp9kdjvjtDJ1ocJIUSb0uSE79VXX6Vbt2788MMPfP755+zatYtRo0aZMzYhOgy9Xk9KTjF5JRUdNvFxd7DBSqlAp4ecogp6+zjyxpRwAD7bfY7nNhxGJ9uuCSGERTR5Dt+cOXOwtbWlR48efPzxx3z88cd1tvv6669NFpwQHUVxhZYxr20DIOmFmA5ZrkSlVODlaMOF/DLS80vxcdYQ09eHV28dwLNfHeSz3ecor9Tx6q0DUCkVlg5XCCE6lSZ/6kybNg2FQt6khWiJ7MJyAOzUqg6Z7FULcLXjQn4ZZy+WGOcq3j44EBsrJbO/OMBX+9Ior9Ly+u3hqK1aNIVYCCFECzT5k2fNmjVmDEOIji2roAwAL0cbC0diXr18HNiTconjmYU1jt8Y7o+NlZLHPt/PDwfTSblYzOu3h9PLu2OuWBZCiLZGvmIL0QqyLvfweTlqGmnZvo0P8eahMd25qmftbdUm9vPlvWmDcbGz5vD5Aq5b8RuvxB0jp6jcApGKvyssqyQlp9jSYQghzEQSPiFaQXXC5+nUsXv4xoZ48ezEEKK6u9f9fG8vfn7iKsaFeFGh1fHvbacZ8fIWYr89TFpuSStHK6rp9XqmfbSHD347Y+lQhBBm0nEnEwnRhmQVdo4h3abwctLw4fTBbErKZOW20xxIzePj+LP8d/c5xoZ4MaK7O1HdPejh5SCLO8zs3V9P8/6OZMIDXdh/Lg8XW2vKq7TYWKksHZoQwsQk4ROiFWQXdI4hXYCLReUcOp+Ps611vUWmFQoFE/r6cHWoN7tOX2Tl1lPsOn2RTUmZbErKBMBapcDfxZZANzt8nDT4OmsIcLXDy8kGR401Xo42+DprsFLJQEVLHEzLY8lPxwD45ajhz9zbSVMr2dPr9VTp9FjLn7MQ7ZokfEK0glG9PNCoVYQFOls6FLNb+0cqS/93nMnhfo3uKqJQKBjRw4MRPTw4ciGfbcez+f3MRfam5FJaqSXlYgkpF+sf6lUpFfg6a3Cxs0atUqK2UqK2UqGxUmJvY4WjxgonjTXOttZ4Otrg7qDG1U6Nn4strnbWnbryQB9fJ1bdPYgHP91nPDYtKog/Ui6x7o9U9p/LJdTPmfjTORSUVnFPVFceHN0dT+mlFqJdkoRPiFZw08AAbhoYYOkwWkXvyytvj2UUNtKypr5+zvT1c+aRsT3Q6vSk55eSeqmU1NwSsgrKuJBfxvncUrILyykqryKjoIyKKh1puaWk5ZY2O04bKyWOGmucbK3wcdLg46zB20lj/Lm7pz1d3e07bM+WtUrJxH6+/Pj4SJ5cl8ionp6E+jmx/JcTfJmQBsDp7D8XcXz4WzIZ+WWsvGuQpUIWQlwBSfiEECbV28eQ8J3OLqJKq2vRkKtKqSDA1Y4AVzuiqHsBiE6nJ7uonHOXSigqr6KiSmd8lFZqKS6voqi8ivzSSvJKKskuLCe3pIKcogpyisopr9JRXlROTlE5Z7LrXp1qpVQQ5GFPTy8Heno50MPbkR6eDnT3su8w89z6+jnz85Ojjb8/Pq4nZZU6Vv16mmHd3Jg6tAtpuaV8tS+NJ6/uacFIhRBXQhI+IcysSqsj5WIxno4anDRWHX4Y0d/FFmuVgkqtnszCcvxdbM1yH6VSgbeToVeuucqrtGQVGHoKc0sqyCwoIz2/jKyCcjLyy7iQX8rprCKKK7ScyiriVFYRP/3lfGuVglBfJ/r5O9PP35n+/s6E+jqhbCeLTJ5efwB/V1umRwXhaq+u8ZxSqWDONSHMuSakxvFHxvZozRCFECYmCZ8QZpZRUEb069tRq5Qcf3GipcMxO6VSgY+zhtRLpaTnlZot4bsSNlYqAt3sGmyj1+tJzy/jZFYRJzMLOZVVZPy5oKyKA2n5HEjLN7b3crThql6ejOrpQWSwOz7ObWuBjk6nR6lUoNfr+Tkpk/zSSqL7eNdK+Bpz5EI+W49lMaKHR6NzNIUQbYckfEKYmbEGn6NNh+/dq+bnbEvqpVLO55Uy2NLBtJBCocDPxRY/F1tG9/qzkLReryctt5QDaXkcPl/A4fP5JKbmkVVYzpcJacb5b0OCXJk+PIgJoT4W30Zu6f+O8cGOZF6/PZzBQa7kl1aiUiro4eXQ7Gv9J/4sa/9IJT2/TBI+IdoRSfiEMLOsgj8Tvs7C73Kv3oW8MgtHYnoKhYJANzsC3ey4boAfYBgi/iM5l+0ns9l1OoekCwX8kZLLHym5eDiouT7Mj1sGBdDP3zKrtLu42dHP35l//3qKpzWGodogdzs01s2fh3jdAD/W/pHK1/vOc32YH/39nbG3kY8SIdo6+VcqhJll5BtWkPq0YK5Ze3X74EBG9fQgPNDF0qG0ChsrFSN7ejCypwcAmQVl/Pf3s3y2J5WconJW70xh9c4U+vs78+i4HkwI9W7V3t5h3dx59qtDqFVKDp83DEOH+Di16FpR3d3xdrIhs6CcO977nRAfR75+eDh2avk4EaIt65j1BoRoQ85dMiR8gW5tby6buUR1d+fmQQF082z+kGFH4O2kYfaE3sTPHcdH9w7mugG+qFVKDp3P5x//SeDmf+/i9zMXWy2eLm52OGmsqNDq+P7ABeDP1dTNpbq8qKN6mPpYRiFf7TtvsliFEOYhCZ8QZpZ6eY/YxhYJiI7HWqVkXIg3b985iN/njeeRsd2xtVax/1wed7z3OzPX/EFmgXmHvX87mcMfKbkEX06+q+sj9vJuWcIHhrqSe+dHc/tgQ23Jw39ZvCKEaJukD14IM0u91PkSvrJKLfFnLnKxqIJbIzpHwenGuNmreTomhOnDg3hr8yk+33OOzceyiFm+nRcn9zPOBzSFvJIK/r3tNLZqFXGHMziWUUiAa80e5pAW9vBVc9JY88+Y3jw2rqdxJXZWYRl7ki9xVS9PnDTWV3R9IYRpSQ+fEGZ2++BApg4NpGcLVkS2V2WVWmas/oN/rj9AWaXW0uG0KV6OGhZP7sdPs0bRz9+JvJJKHv1sP7PW7ie/pNIk98gtqeR0djHLfzlp7NGbOrQLAIO6uLDjmbF0McEXEC9HDYFudsb6gwdT81m/N43F3ydd8bWFEKYlCZ8QZnbfyGCW3DyAANfO08PnbGuNxtrw9pKR3/FW6ppCT29Hvn5oBI+N64FSAd8mXmDim9s5mJZ3xdcO9rAn1PfPHjwHGysmhHrjbGuNm726RpJmSl3c7bBTq9h2Itvk1xZCXBlJ+IQQJqdQKIwJbvUcRlGb2krJUxN6s/7B4QS525GeX8aUd39nU1LmFV/7wTHdsVcbyq4M7+5ODy8HEhdezQfTh1zxtf9q7Z5zPPLffew6lUOAqy0/Hc4gu9CwZZ0Qou2QhE8IM8osKON4RiElFVWWDqXVdb08ZHjukiR8jYno6sr3j43kql6elFZqeeA/e1m9M7nZ11mx+SSTV+7kh4MXsFNb8e49g4nu480T0b1QKBRmKQWzJ/kSPx5K56OdyRSUVhmHik9kFpr8XkKIlpOETwgz+mb/eWKWb+e5bw5bOpRW18X9csJ3URK+pnDUWPPh9MFMHdoFvR6e/z6J578/glanb/I1jqYXkJiaR+blYt8je3rwwfTBhPq1rOZeU/TxNVz7l6NZDFuymbTLPbqnsorMdk8hRPNJwieEGVXPX2tr+6q2huoevrOS8DWZtUrJv27qx5xrDLthrN6ZwkOfJlBRpWvS+SmX/6yDPVpvvujfk8nJA/0B+f8uRFsjCZ8QZmRM+DrRLhvVurrbA3BWhnSbRaFQ8ODo7qy8cxBqKyU/J2Xyz/UH0DXS06fX6zl7sRiAoMt/9q1hQIAz1irDUHF4oItxdxVJ+IRoW6QOnxBmlHG5qK53J0z4+gc4s+y2MLp5tl7y0ZFMGuCLg8aKmWv+4LsDF3B3UBN7fd9622cVllNSoUWpoFVXhDtqrHl/2mA+iT/Lo+N6UFBaadyFQwjRdkjCJ4QZVe+i0BmHdD0cbLhFii5fkdG9PFl2exiz1iayemcKfs623H9VtzrbJucYevcCXO1aPeEa09uLMb29ANDq9Bx7YaJZyr4IIVpOvoYJYSZanZ6sQsPk+c44pCtM48Zwf+Zda5jT99LGo3ybWHPf2qzCMtbsTOZAah4AQR6W7VFVKRWS7Ik2ZeXKlQQFBaHRaIiMjGTPnj0Ntl+/fj0hISFoNBr69+/Pxo0bazyv1+tZuHAhvr6+2NraEh0dzcmTJ+u8Vnl5OeHh4SgUChITE031klpEEj4hzORiUTlanR6VUoGno42lw7GIg2l5fPr7WQ7JXqtX5P5R3ZgxIgiAf64/wB8pl4zPPfrf/Sz6PoklPx3D38WW7jKELoTRunXrmD17NrGxsezbt4+wsDBiYmLIysqqs/2uXbuYOnUqM2fOZP/+/UyePJnJkydz+PCflRZeffVVVqxYwapVq9i9ezf29vbExMRQVla7yPwzzzyDn5/ptk28EpLwCWEm1iol/5zQi/tHdUPVSXs8Pv39LPM3HGbb8brfXEXTKBQKFkwK5dr+PlRq9Tz62T6yC8spq9Sy5y/J37anx7DwulALRmrw1uaTXPfWDr4/cMHSoYgOqLCwkIKCAuOjvLz+It+vv/46999/PzNmzCA0NJRVq1ZhZ2fHRx99VGf7N998k4kTJ/L000/Tp08fFi9ezKBBg3j77bcBQ+/e8uXLmT9/PjfeeCMDBgzgk08+4cKFC2zYsKHGtX766Sd+/vlnXnvtNZO99ishCZ8QZuJqr+bRcT2NJTY6Iz8XWwDO55VaOJL2T6lUsPTWMLp72pNZUM7D/03gaHqB8fnunvZk5JeZpbhyc+WWVHL4fAEv/pjEmWypxydMKzQ0FGdnZ+NjyZIldbarqKggISGB6Oho4zGlUkl0dDTx8fF1nhMfH1+jPUBMTIyxfXJyMhkZGTXaODs7ExkZWeOamZmZ3H///fznP//Bzq5tbKspCZ8Qwmz8JeEzKXsbK96bNhhHGyv+SMnllbhjAAS42vLL7NEEurWND5ZRvTwAyCwo5/Z3f6e0QmvhiERHkpSURH5+vvExd+7cOtvl5OSg1Wrx9vaucdzb25uMjIw6z8nIyGiwffV/G2qj1+u59957efDBBxk8eHDzX6CZSMInhJmk5BRzPKOQ4vLOt61aNX/XywlfriR8ptLd04GpkYEA/H7GMJwb3ce7TfTsVRvZw4ObBxkKMOcUlZNwNtfCEYmOxNHREScnJ+PDxqZtzZF+6623KCwsrDcRtRRJ+IQwkxVbThKzfDsfx6dYOhSLCXAx9DidzytFr2/6FmGifpVaHZ/EnzX+bqVU0NfP0YIR1WatUvL67eFMCDX0gpzKkn11Revz8PBApVKRmZlZ43hmZiY+Pj51nuPj49Ng++r/NtRmy5YtxMfHY2Njg5WVFT169ABg8ODBTJ8+/cpfWAtJwieEmRhr8HXikiw+zhoUCiiv0nGxuMLS4XQIqZdKKKs0bLXm4aCmSqfnw99SKKtse8OmwZdXDFfXCBSiNanVaiIiIti8ebPxmE6nY/PmzURFRdV5TlRUVI32AJs2bTK2Dw4OxsfHp0abgoICdu/ebWyzYsUKDhw4QGJiIomJicayLuvWreOll14y6WtsDim8LISZdOZt1aqprZS429uQU1RORn4ZHg5ta+ilPUq5vH2av4st3zwynGuW7+BYRiEv/3SMRTfUvxOHJQwMdCG6jxd9fJ0abyyEGcyePZvp06czePBghg4dyvLlyykuLmbGjBkATJs2DX9/f+PCj1mzZjF69GiWLVvGpEmTWLt2LXv37uW9994DDCvmn3jiCV588UV69uxJcHAwCxYswM/Pj8mTJwPQpUuXGjE4ODgA0L17dwICLFeMXhI+Icwks8BQKsC7E+6y8VdLbxuAxkpFsIULAncUZ7INCV9YoDNejhpeuy2MGWv+YM2uFK7q5cG4EO9GrtB6JvbzZWI/X0uHITqxKVOmkJ2dzcKFC8nIyCA8PJy4uDjjootz586hVP452Dl8+HA+++wz5s+fz7x58+jZsycbNmygX79+xjbPPPMMxcXFPPDAA+Tl5TFy5Eji4uLQaNr2e71C38km1qSlpREYGEhqaqpFM23RsRWWVdJ/0c8AHHk+Bnsb+W4lTGP+hkN8+vs5HhnbnadjDCV/nv/+CKt3puBmryZu1ii8OnGvsui45PP7ysgcPiHMoHr+nqONlSR7wqQKSqtQKCDI/c8e0znXhNDH14lLxRXM/uIAOl3b+R6v1+vJLiynuLxKFu4IYUHySSSEGaRdLkNSXXi4MzuVVcjvZy7h66xhfJ+2M9zYXq2YOpBXbx1Q45iNlYq3poZz3Vu/8dupHFbvSmHmyGALRVjTtI/2sONkDmqVEmc7a965axBDgtwsHZYQnY708AlhBl3d7Xk6pjd3RnZpvHEHF3/6IvM3HGbdH6mWDqXD0Fir0Firahzr4eXI/EmGbdWW/u8YZy+2jZWxXS4Xg67Q6sguLOeDHWcsHJEQnZMkfEKYQbCHPY+M7cH04UGWDsXiPB0N88myCuvf71I0rqCsEm0jQ7V3RXZheHd3yip1PPPlwTYxtHtDWM2N43eduihDu0JYgCR8Qgiz8nYylGLJujyvUTTfgdQ8Il/aTPd5G7n5nZ2s3HqqznYKhYJXbhmArbWK3cmX+O+ec60caW2R3dxZfe8QVs8YgtpKiZuDmtySSkuHJUSnIwmfEGaw9XgWp7OLGu2R6QyqV4xmF5VLz04LfbP/PKWXCyvvO5dnLM1Sl0A3O56d2BuAlzceJS23pFVibMjYEC/G9vbiYOwEfn16LG72akuHJESn0yYSvpUrVxIUFIRGoyEyMpI9e/bU2/b9999n1KhRuLq64urqSnR0dIPthWhtZZVaZq75g/HLfiWnSIYxPS8XW67U6qVnp4UKymr+uYUFOjfYflpUEEOCXCmu0DL360NtJtH++7xDIUTrsXjCt27dOmbPnk1sbCz79u0jLCyMmJgYsrKy6my/bds2pk6dytatW4mPjycwMJAJEyZw/vz5Vo5ciLrtO5uLTg9ejjZ4OcrOEmorJa521gBkFcqwbku8fns4BxZOwMNBjbOtNdGNrHZWKg1DuzZWSnaczGH93rRWilQI0VZZPOF7/fXXuf/++5kxYwahoaGsWrUKOzs7Pvroozrb//e//+Xhhx8mPDyckJAQPvjgA+PeeHUpLy+noKDA+CgslE28hXn9nnwJgOHd3VEoFBaOpm3wql64USA9ni3lbGfNjmfGsWvOuCaV++nm6cBTE3oBsPjHJONWf5Z0OruI21fFM3nlTkuHIkSnY9GEr6KigoSEBKKjo43HlEol0dHRxMfHN+kaJSUlVFZW4uZWd12nJUuW4OzsbHyEhoaaJHYh6rMn+SIAQ4PdLRxJ27Hohr6sfWAYYYEulg6lXbNVq5pVyHvmyG6EBbpQWFbFc99YfmjXSWPNnpRLHEzLo7RCa9FYhOhsLJrw5eTkoNVqjXvaVfP29iYjI6NJ13j22Wfx8/OrkTT+1dy5c8nPzzc+kpKSrjhuIeqj0+k5lJYPQERXVwtH03ZEdXdnWDd3nG2tLR1Ku7N+byoT3viVFZtPNvtclVLB0lsHoFYp2Xwsiw2Jlp364ulog7+LLTo93PTOTvZc7g0XQpifxYd0r8TLL7/M2rVr+eabb+rdtNjGxgYnJyfjw9HRsZWjFJ3J+bxSiiu0qFVKunnaN35CJ7PvXC4H0/IsHUa7ciq7iBOZRVwqrmjR+b28HXl8fA8AFn2XZPF5lA+N6Q7AsYxCpr7/Owlncy0ajxCdhUUTPg8PD1QqFZmZmTWOZ2Zm4uPj0+C5r732Gi+//DI///wzAwYMaLCtEK3leIZhjmg3T3usVe36+5TJFZdX8c8vDvDsV4csHUq7cu6ioaxK9Y4VLfGP0d3p6+dEfmklsd8eMVVoLXJXZBe+emg4w7q5odXp+bwN1AoUojOw6CeSWq0mIiKixoKL6gUYUVFR9Z736quvsnjxYuLi4hg8eHBrhCpEkwzq6sqquyN4IrqXpUNpc9LzSzmTU8zR9AIpwtwM5y5decJnrVKy9NYwrJQKfjqcwY8H000VXrMpFAoiurry6NieAOw4mW3xuYVCdAYW74KYPXs277//Ph9//DFHjx7loYceori4mBkzZgAwbdo05s6da2z/yiuvsGDBAj766COCgoLIyMggIyODoqIiS70EIYzc7NVM7OfDxH4N91B3Rj28HOntbZhSceRCgYWjaR/0er2xh6+re8sTPoBQPyceHmsY2l347eEWDxGbyuAgV+6K7MKi6/si9cmFMD+LJ3xTpkzhtddeY+HChYSHh5OYmEhcXJxxIce5c+dIT//z2+i///1vKioquPXWW/H19TU+XnvtNUu9BCFEE/XwdgDgVJZ8QWuKvJJKCsurAMMOGlfq0bE96O3tyMXiChb/YNkFbBprFS/d1J9r+vuiUkr5IiHMrenr+83o0Ucf5dFHH63zuW3bttX4PSUlxfwBCdECpRVaPthxhlA/J8aFeEkNvjp095SErzmqh3O9nWxMskuF2krJK7cO4KZ3dvLN/vNMHdqFocF1l7QSQnQsFu/hE6KjOJZRwLJNJ3j2q0OS7NWj++WVy8k59e8FK/5kq1YRHuhyRfP3/i480IU7hnQBDEO7VVqdya7dXOVVWn4/c5Ev9qZaLAYhOos20cMnREdwMtPQaxXiI6V/6uN/eYeIC/mlFo6kfejl7ciGR0bU2kv3Sj0d05ufDqdzLKOQT38/y70jgk16/abKL6nkjvd+R6GAif18cNJInUYhzEV6+IQwkVPZhoSvh5eDhSNpu7q62zN1aCB3RnaxdCjtiqkTITd7Nf+c0BuAZZtOkFNkmS3vvJw0dHW3Q6+HP6QIsxBmJQmfECZyMtNQg08Svvp5Otqw5OYBPDymh6VDafOOZRSQV2K+lbRTh3ahn78ThWVVvPLTMbPdpzEjengAsP1EtsViEKIzkIRPCBOp7uHrKQmfMIGn1x9k6Eub2XIss/HGLaBSKnj+hn4ArE9Is9iOF1f19ARg+8kci9xfiM5CEj4hTKC0QktarmFemvTwNay4vIpTWYVkF1pmGLE9OJiWx6Hz+ej0esICXMx2n4iurtwaEQDAou+OoLNAQbzhPdyxUipIzikm9fKqZCGE6UnCJ4QJnM4uQq8HVztr3B1sLB1Om/bMVweJfn073yaet3Qobdb7O5IBuD7Mz+x/n56dGIKjjRWHzuezPqH1V8s6aawZ1MUVgF9lWFcIs5GETwgT6OXtyE+zRvHW1EGWDqXN83PWAJCeL9ur1aWgrJKfj2QAcF8rrJ71dLRhVrRhm7NX446TX2raFcFNMbq3YVj3QGpeq99biM5CyrIIYQJqKyV9fJ0sHUa74OtsKM2SLqVZ6hR/+iLlVTqCPezp5986f6emRQXx+Z5znM4uZtnPx3nhxn6tct9qN4T5MSTIjcFdXVv1vkJ0JtLDJ4RoVX4uhh6+C3nSw1eX+NMXARjRw73VCnirrZQsvpzk/ef3syS2ck9boJsdQ4PdUMoWa0KYjSR8QpjA0v8d4/3tZyy+IX17ID18Dfv9jCHhG9bNvVXvO7yHBzcN9Eevh3lfH7LYDhx6fesvHBGiM5CET4grVFJRxTvbTvPSxqNU6Sy3TVV74Xu5hy+rsJxKC27r1VZ9eO8QXr89jOHdPVr93s9N6oOzrTVJ6QV8HH+2Ve9dqdXx7JcHGfnKVvJLWn8eoRAdnSR8Qlyho+mF6PXg5WiDl6PG0uG0eR72NlirFOj1kFkgw7p/5+9iy82DAnCzV7f6vT0cbJhzTQgAr/98vFV7Ya1VShLO5XI+r5Qdp2S1rhCmJgmfEFfoyIV8APr6yaKNplAqFfzfqG48dXUvNNYqS4cj/mbK4EAiurpSXKFl0XdHWvXeYy+v1t2UZJ5i00J0ZpLwCXEFjmcUsuznEwD093e2cDTtx7MTQ3hsfE88pGahUVmllllr9/Pe9tPkWnAuqFKp4KWb+mGlVPC/I5n80orJ17X9fQFDwldaoW21+wrRGUjCJ0QL6fV6Hv5vAvmllaiUCm4I97N0SKKdOpSWz/2f7OXbxAt8vOsstmrL9nyG+Dgxc5ShBmDsd0coqahqlfuGB7rQxc2OkgotGw+lt8o9Rce3cuVKgoKC0Gg0REZGsmfPngbbr1+/npCQEDQaDf3792fjxo01ntfr9SxcuBBfX19sbW2Jjo7m5MmTxudTUlKYOXMmwcHB2Nra0r17d2JjY6mosOyiPkn4hGihU1lFnM4uRqGA7x4dQQ8vR0uH1G5UaXWcyipib8olS4dicQdS85j8zk52XN5Ldv6kPm1iqHvW+J74u9hyPq+U5b+cbPwEE1AoFNw+2LDV2+d7zrXKPUXHtm7dOmbPnk1sbCz79u0jLCyMmJgYsrKy6my/a9cupk6dysyZM9m/fz+TJ09m8uTJHD582Njm1VdfZcWKFaxatYrdu3djb29PTEwMZWWGOcnHjh1Dp9Px7rvvcuTIEd544w1WrVrFvHnzWuU110eh72Rr4NPS0ggMDCQ1NZWAgABLhyPascKySnafuUReaaVxP1LRNH+kXOK2VfH4u9iyc844S4djUf/38R/8cjQLDwc1r98ezlW9PC0dktGWY5nct2YvKqWC7x8dSWgrzFPNLChj+Mtb6OPryJcPDm8Tya9oG1ry+R0ZGcmQIUN4++23AdDpdAQGBvLYY48xZ86cWu2nTJlCcXExP/zwg/HYsGHDCA8PZ9WqVej1evz8/Hjqqaf45z//CUB+fj7e3t6sWbOGO+64o844li5dyr///W/OnDnT3JdtMtLDJ0QLOWqsiQ71lmSvBXp4OgBwPq+U4vLWGS5siwrKKtlyzNDTsO4fUW0q2QMYF+LNNf180Or0PLfhEDqd+fsHvJ00TBkSiFaHJHuiToWFhRQUFBgf5eXldbarqKggISGB6Oho4zGlUkl0dDTx8fF1nhMfH1+jPUBMTIyxfXJyMhkZGTXaODs7ExkZWe81wZAUurm5Nfk1moMkfEKIVudqr8bT0bBg40RmoYWjsZw/ki+h00Owhz3dLyfBbU3s9X2xV6vYfy6Pz/9onWHWF2/sxzt3yb7Uom6hoaE4OzsbH0uWLKmzXU5ODlqtFm9v7xrHvb29ycjIqPOcjIyMBttX/7c51zx16hRvvfUW//jHPxp/cWYkCZ8QLfRJfApxh9M7dQ/VlQjxMcx5PJ7ReRM+bycNU4d24Yawtrvgx8dZw1MTegPwyk/HyC6suzfFlJRKBcEe9hSWVbJmZzLPfnnQ7PcU7UdSUhL5+fnGx9y5cy0dUr3Onz/PxIkTue2227j//vstGoskfEK0QEWVjue/T+LBT/dRUCa7ArREL29DwnesEyd8/fydWXJzf568upelQ2nQ9OFB9PN3oqCsihd/TGq1++qBF388yrq9qaTkFLfafUXb5ujoiJOTk/FhY1N3eScPDw9UKhWZmTVLC2VmZuLj41PnOT4+Pg22r/5vU6554cIFxo4dy/Dhw3nvvfea/gLNRBI+IVrg7MVitDo9DjZW+DjJ7hotMSDAULdwd7Ks1G3rVEoF/7qpP0oFfJt4gR0nW2cnDCeNNUOCDPOequc6CtFUarWaiIgINm/ebDym0+nYvHkzUVFRdZ4TFRVVoz3Apk2bjO2Dg4Px8fGp0aagoIDdu3fXuOb58+cZM2YMERERrF69GqXS8umW5SMQoh06nV0EQHdPexQKhYWjaZ9G9fREqYCj6QWdsvemvEpL0oWCVqtxd6UGBLgwLSoIgAUbDlNW2TqFkceGGBaybG+lJFN0LLNnz+b999/n448/5ujRozz00EMUFxczY8YMAKZNm1ZjSHjWrFnExcWxbNkyjh07xqJFi9i7dy+PPvooYCgd9MQTT/Diiy/y3XffcejQIaZNm4afnx+TJ08G/kz2unTpwmuvvUZ2djYZGRn1zvFrLVYWvbsQ7UxZpRaNtYpTWdUJX9ucaN8euNmreWpCbwYGutDV3c7S4bS6k5lFXPfWb3g4qNk7/2pLh9MkT03oxU+H00m5WMI7W08x+/LcPnMa1dMTOMbvZy4a//0J0VRTpkwhOzubhQsXkpGRQXh4OHFxccZFF+fOnavR+zZ8+HA+++wz5s+fz7x58+jZsycbNmygX79+xjbPPPMMxcXFPPDAA+Tl5TFy5Eji4uLQaAyjPZs2beLUqVOcOnWqVvkYS1bCkzp8QjSRXq/ntlXxPDMxhM/3nOOb/ed5OqY3j4ztYenQRDv03YELPP75fiK6uvLVQ8MtHU6TbTyUzsP/3Ye1SsFPs66ih5d5v/To9Xoi/7WZrMJy/vt/kYzo4WHW+4m2Sz6/r4wM6QrRRGcvlrD3bC6z1u4n6UIBYBjSFaIlTl/uJe7RznqJr+nnw5jenlRq9czfcMjsPRYKheJyLx9sPyHDukK0lCR8QjTR8cv14pw01nS5PATZz9/ZkiF1CL+dzOGVuGMcPp9v6VBalXFagFf7+tKgUChYfGM/NNZKfj9zia/3nTf7Pa/q5YFapaRQSiAJ0WIyh0+IJjp5OeHr6+fEvEl9qNLqCHDtfHPPTO3zPef48VA6jhqrTpNA/3DwAj8eSgfa5zzQQDc7Hh/fk1fjjvPSxqOMC/HC1V5ttvvF9PXh6lBv7NTykSVES0kPnxBNdLJ6CM7bAQ8HGz6YPsTCEXUM4YEuACSey7NoHK1p1+mLxp/DLr/+9ub+Ud3o5e3ApeIKnjPz0K7GWiXJnhBXSBI+IZqoekeI3pcLBquUUo7FFMK7uACQmJpn0RVsreFkZiHzvjlEkLsdE0K9eXZiCB4OdReNbeusVUpeuy0MK6WCjYcy+DIhzez31Ov1LPruCIfSOtfwvxCmIAmfEE1QUaXjTLahVlz1DhHCNPr5OaNSKsgqLCc9v8zS4ZhVYmoen+0+x7bj2bw3bTAPjelu6ZCuyIAAF+MuIYu+O8LZi+atp7jw2yOs2ZXC018e6PBfDoQwNUn4hGiCE5mFVGh1uNpZE+Bqa+lwOhRbtcrYa5qYmmfZYEzsj5RLxB3+s9hq9bSAnmYuZdKaHhzdnaFBbhRXaHlyXSJVWp3Z7nXXsC6AYTu+D3Ykm+0+QnREkvAJ0QT9/J1JmB/Nh/cOkZ01zKB6WPdAB0r4corKufP933nw0wR+P2OYs1c9LaBnB+olVikVvD4lDEcbK/ady+Pf206b7V4hPk48Pr4n/9/efYdHUe0NHP/upod00gNJKCmUEHpuRAQFKSoCFhDQVxT1egVRFK/gFQER4SoqRYqKUrwKigIWioQAAUIIkBAIIQmQQnohve2m7Hn/2DCwAgpKSDuf55mH3ZkzM2c4u9nfnDkFYHVY0h2b7UOSWgIZ8EnSTWprZUZvT/vGzkaLdLnjxrn6ntAtQVxWKTV1+seOe8/mIoQgLkvf9qyru01jZu22a2dvyfzR3QBYFnqe0xnFDXau6fd1xsnajMKKaqZvOkmdTj7alaSbIQM+SfoDyfnlrD6QRElVTWNnpUUb3tWV/TMH89XkltPzOb2wUnl9/GIRGUVVXCqvRq2CLq4tK+ADGNvLgwd7uFGrE7z6XQxV1Q1T+2ZspOb9sQEA7Dmby/RNJ9HJoE+S/pQM+CTpBo5cuMTUb0/y390JvLU1trGz06LZWprQwbFNi3pcnl50JeA7k1nCGz+cAqB/BwcsTFvefLAqlYqFY7rjamNOcn4F83+Ja7Bz3d/VhTkPdQWgpk5Hja7h2g1KUkshBzaSpOs4kVrIxLWRyvvn7+nYiLlpXXQ6gboFDHkz/T4fRgd6sDEilcF+zkQkXeJUegkzhvo2dtYajJ2lKR+NC+TJLyPZfDydvt4OPNanYeY8nXJ3B3p72tGjnR1GahUF5VpCE/IY17d9g5xPkpo7GfBJ0nUcTy1SXk8f4qO0MZMaTlZxFYt2JZB6qYKfpw1o9rV9bcyM6epuw+JHewAworsrr93vh62lSSPnrGEN6OzIjKG+fBxyjv9si8XH2arBBpfuVd+mtkxTwz0f7Keiuo4BnR3xsJM96SXp9+QjXUm6jqR8/fAZM4b68tr9LbdGpimxMjfmtzM5xGaWkFpQ+ec7NEMtPdi7bNq9nRni74y2VsfzG0+QXVLVoOezNjehc/1QN8dTChv0XJLUXMmAT5Ku4/LE9p1b0HhpTZ2NuQmB7fVz6R5LKfiT1E1bmaaGt7bFsiYsqVUOEKxWq1j6RE98XazIK9My+avjDd7x6XJt3+We0JIkGZIBnyRdh6OVGY5WpjLgu8P6d3AAILKZ1dJkFVcxdlU4U7+JJqdEQ1phJd9GpvHFweRm/2j6r7I2N+Gryf1wtjYjMbeMf359Am1tw42b5++qH9swoX6sw8KKap5df5zNx9Ia7JyS1JzIgE+S6qUVVPL29liOJhew9um+nHj7fnxdZMB3JwV1aAtAZHLzCvjWHkrhZFoxO2KzUamuDMnSzsGykXPWuNrZW7LumX5YmRlzNLmQmVtON9gQKv5u+qFuLgd83x1P50jSJbbHZMqx+iQJGfBJkmLp3nP872gaR5KuPE5srbUzjaWPlz1GahWZxVVkFDWfdnwR9TNpjO3lgYuNudIG0bOVB3wA3dxtWfNkH4zVKn45lcWiXfENch5fFytUKsgv03KpXEtofC6aGh0PBLhh1AJ6fUvS3yUDPklC32avt5e+DdDlabCkO6+NmTEBHvp2fFcH3k1ZqaaGhJxSAGY/4A/AufpaJl/ZJACAu30c+fBxfW/lLw6lsOrAhdt+DktTY7zqA+yIpAKi0/Q97Yd0cbnt55Kk5kgOyyK1eumFlQz9OEx5fyylkF7v7mHWSH/G9/NsxJy1Tvf5O6NWgblJ8xicOCW/AiHA2doMZ2tzskuq2HoyE4Du7WwbOXdNx9he7cgr1bJoVwIf7E4E4KXBnW/rOWbc74uRWkVOiQadgC5uNtiYG7MvIZfAdnYUVdbIdrlSqyUDPqlVqayuZc72OHq2t+WpYG8AQs7mAtDf24HM4ioyi6soqqzBztK0EXPaek0f4sP0IT6NnY2bdrG+vZ5XW33t0o9RGcq2oPpOKJLePwd1orpWx0ch5/hgdyJCwNR7b1/QN7qnBwAvfRMFwGA/J+5avI8yTa2SZv0z/Rjs53zbzilJzYV8pCu1Kqv2J/FjdAZzfoqjXFuLTif4MVr/Az0ywJVHeut/MDzsLBjk69SYWZWaibSCCgDa1z9OnBTkxWA/J94fG4Clqbyn/r2Xh/jwev3Ylh/+lsin+87f1qFrLpVrSS/Uj/s3vJsr9/4uuPv8YPJtO5ckNSfyr5HUapRqalgXngLAigm9qK3T8fKPp4nLKsXS1IiHA91pY2aMs405g3ycms0jxZaqTFODpkaHk7VZY2flGkUV1XwTeZFh3VwZGeCGo5UZ7ez1AZ99G1PWP9O/kXPYtL08xAeVCpbsOceSPefIL9Py9kNdMTH6+3UQZZpaunvYYmFiRGA7W56+y4sdsdlKT90jSQWkXqrA27HN3z6XJDUnMuCTmrUTqYXM3HKKGff7Ko9zbmTn6WwqquvwcbZiRHdXxn0Wwcm0YgDmjepGWyt9YPHUP7waOtvSn1i5/wIf7Unk/4K9mfdwt8bOzjWiLhaxfN8FNh1LZ9/MQXRykm09b9W0+3wwMzZi4c54NkRcJD6njE8n9sLZ2vxvHfdcbhmamjqWTeiJSqWij5cDO6cPxEgNP0Zn0sXNBnc59ZrUCslHulKztuDXs6QWVPLK5pg/fSx08Hw+AA/1cCevTMuZzBJMjFSsebI34/rJCdebEq+2lugEbI/JVKa5A8goqqSq+srgvUUV1ZzJbNiZFbKKqxi14jDl2ivtwO7zd8bB0pTM4qpm05u4KXr+no589lQfrMyMOZZSyKgVh4m6+PfGYBzezZVPxvfEzfZKUOfnak1nZ2veHOHPw4HumBrLnz6p9ZGfeqnZEkIQn12mvE++pG9LVaGt5e3tsby1LZaHPz3M2FXhJOWXE35B/8M80Fc/uXpvT3s+GteTEd3dGiX/0o0N7eKCn4s1xZU1PPH5UYorq9mXkMvd/93PC1+fUNK99E00D604zM7Y7AbLy87YbGIzS9gWnaFM26VWqxjSRd827PkNJziR2rwGim5Khndz5adpA/BxtiK3VMvjayJYvCsBTU3DzcpxPUKIVjkNntR6yIBPatYOvDFYeX0spZDCimqe+jKS/x1N49vINE5nlJBeWMn53DJKqmqwMTemR/04b1883ZeHA90bKefSHzE3MWLNU30wMVKRX6Zl15kcvjmqnyLr0PlLFFZUA/Bk/eP3LSfSb+v5Y9KL2XwsjYyiSrS1OgDm/BTHg8sP88JGfcA5tKt+fLdanWD6ppO39fytTScnK7ZPHcAjvTzQCVgTlsQDyw8R2QBjYp7NKuXzg0kcPn+JqIuFRKcVEZtRwoDF+xjyURjH7uC0fiWVNWQVV92x80mtmwz4pCYrp0TD29tj2ZeQe93tKpUKdzsLXr5PP6zDsZRCnll3jOj6dnkAahUsHBtAO3tLHghwZWR3N4zrG4bbmJs0+DVIf10HxzZMu1c/PMvBc/lE1Q+kC3A6oxgAdzvz+vd/77Fuda2O2jqd8v7HqAxmbY3li4PJ19wUXB5+JbhjW2VdJzm229/WxsyYj8f35POn+uBsbUZyfgXjPz/K7K2nKamquW3n2RGbxfs7E3hvx1keXxPBkt8Smf9LHFklGpIvVTDuswge/vSwcs6rPxe3U22djoc+PcTAD/azJy6nQc4hSVeTAZ/UZG07mcmmY+n88+soZW5SgLBz+RRXVivv+9ePdbbtZCanM0toY2rEb6/eQ8KCESS+N5Lh3Vzp7mHLqkl9+O9jPe74dUh/3d0+jgDsOpNDceWVH/3Y+gDPx8UagIKKakoq/3pQsOdsDoHz9zDv5zgAjiRdAuCuzo60d7Ckt6edkvapf3gD+lrIFRN60cvTjn8P9//L55YMDevmSshrg5gYpO8Is+lYOkM+CuPzg0kG7Sj/qgAPO0A/5663YxtG93TnxEX9zUTP9vptpzNK2HwsjYJyLY+tibjhTedle+Jy6PnuHmUUgJtxLKWQ9MIq6nSCT/ae/0vXIkm3QvbSbWFySzWsPpDE0C4uyo/lnVBQrqWmTuBqe+s97HafySGvTMOYXh4GtW4vDupIZEoBBxLzWXsomfmjuwPw08lMlu09x7i+7cksrmKQrxMv39eZoA5tcbczJzGnDD9X69t2bVLjCWxni72lCUWVNThZm2Fuombjs0G42pgzed0xurrZYGthQklVDakFFQRa2t3wWBsjUkkrqMSrrSW1OsHI7m7K5/VocgEV1XWoVPpgMim/AiO1in900NfivTHcn1lbT/NkkBeeba/Mjzsq0J1RslnAbWdrYcL7YwMY09OD2VtPk5Rfwfs7E/h03wWevsubyXd5K73qb1VvLzvldXJ+BX6uNqz9v76czS5l+hAf1h5K5r0d8Ww7mUmFtpaY9GL+dzSN+/xvPEVbZnEVxZU1zP/lLGN7edzUoO3Hrmr3GZ9dSrm2Fisz+ZMsNRz56WphZm+N5fCFS6hU3LGAL6OokodWHKZSW8f6Z/txV6ebP+/es7m8+L8oPOws8G7bhnuuGuxYpVLx/MCOHEjM57sT6cwc7odF/dh40WnFnMkqpbpWh1fbNrw+zE/Zr6OTfLzWUhgbqYmYPUQZE7FOJzBSq4jNKOFAYj4n04rxc7HmWGqhPuCrr6G5LOVSBZlFVZiZqHnnpzjDY6tVDPJ1pqBCy/4EfQ/ufQl5rAtPBeAfHR2wtdTfgAR3akvYG/c27MVK1+jfwYFdr9zD9phM1oQlkZxfwYp9F/jsYDIPBrgxob8n/bztUalUN31MZ2tzhndz4be4XO7q1Fap1bvcJnNikCcJOWW8dr8vxZU1LN93gYikAqprdTfs3Tv5Lm++jUzjfF45B89fUpoBCCHYGp2JTgge69POIJ+J9fMt21qYMP/hbhj97hqEEAbpS6pq2BOXw33+zn8Y7FZoa5my4ThV1XXMfbgbvT3trznu6rAkCsqrmTXS/7aMfSg1D7Kk/4b9CXm8vT1WaUDe2IQQHE0uwNXGnCfu4BywXx+9SHFlDdV1Omb9GKv0dMsv0yqDnV4tLquE5zee4EjSJb48nIKRWsWbI/2VYC+9sJLj9Xe/d3VqSyenNmhqdOyKzVHa34G+3RWAv6zNa9GuHgDbSK3/AYytH4olwMNWaVOXVlBp0MsyIaeUYZ+E8eSXkczeGnvNcatq6rjnw/2MXXWEzOIqrMyMeWlwJ2X7iG6uDXI90q0xNVYzrm97QmYMYs2TvQlsZ0t1rY5tJzMZ91kEQz4OY01YEnllmps+5sKxASwY051lT/S6ZpulqTFLHg/E3c4Cf1drHNqYUlVTx6n6dqPXo1KplL9fR6/qaLI9JpPXt5zijR9Oszc+z2CfhPqA79OJvRjTywML0yuf83XhKfi9vZvnN55QHmPP/ekMb/xwmhf/F/WH1xZ2Lp+jyYWcyijhra2x1/Q8PppcyAe7E/nycArf3+bOTqCfvrKp/CZKhmTA9xdpaur41zdRfHc8nZ9jMm95/4Jy7W0fdiCnVENldR1ZxVV0uIOjyIdfuKS8TiusJC6rlJCzufRbuJde7+65ZsiKMk0tWcVVTPwikojkAoQQ9PHS34Vqaup4edNJHl8TwcaIVFQqFY/2aQfAv388zZnMEqVtD0DbNqZ0dbO5A1cpNRWrDyTx1jZ9ANfdw5YBnR0Z17cd3o5t6PveXiavO0ZVdR0/nMigpk4wb1RXfnv1HlZN6q18Vgb5OvFo73aYXnUD8a/BnXi8T3vu8XXC18VKPqptYozUKkZ0d2P71AH8NHUA4/u2x8LEiOT8ChbvSqD/wlBGLjvEvJ/j+PlUFgXl2hsey9HKjKf+4fWns7io1Sqlc87jayLYn5iHEIKIpAK+PnqR1EsVxGeXUlOn4x/16a4O+DYfuxJQbY3OMDj2myP8eWWID93cbQ3W19TpWPJbItV1OkLO5vLV4RTKtbVsj8kC4KXfzT2sqalDW3vlt+RE6pXOTQk5ZZzJLDVIb2thgrmJ/nP/Q/28zz+fymLJb4lKBxUhBAfqr/Xq8yzbe57VB5JuOHxNSVUNQz8Ko//CvRxINAxwdTrBh78lsD48RQ5/00iaRMC3cuVKvL29MTc3JygoiGPHjv1h+i1btuDv74+5uTkBAQHs3LnzDuX0irisEjQ1Ojo5WfF/wd4G2+p0gp9PZRkMGPt7u87k8MiqI7fUA6yooprotKIb7nM+V38+r7aWmBqr0V2ndu2yC3nlzP3pDAk5pTdMA/ov/ryf43hw+SEu5JWjra3jhY0neGub/s6xpk7HuRz9eWeN9OeXaXfT1c2GDUdSASjV6NvARF0s5N1fzrIzNpvenvb4XNWr8T5/FzzsLNDW1vHhb4nEpBdjY27MPT76O+ZJQV7YmOtbH6w/kkpfbwf++2gA7ewtmPNQV9Tqm3+cIzV/Vw/BEtTBgTG9PPjgsUCcrM0oqKjmQGI+4z+PYO1hfQN6J2tzjNQqHghw47t//oOl43vy6cRetLUyY8GYbtzj68S3zwUx9d7OqNUqNj7bnz0zBt1UOyzpzlOpVAS2t+O/j/Xg2H+GsPiRAHrVd6qJzy5l/ZFUpm86Sd+Fexm9MpxPQs5x+PwlyjR/rVPPIL8rzUzWHkpmTVgyL30TxZztZ7jvowOMXHaI3gtCcLczR6XStwvMKKqkuLJa6QwCEJqQR+lVeRjR3ZUZ9/tiaWrEkQuX+LE++IrNLKHiqsHFvzuefqVXuq25wdzAFwsqGPThfgLm7uHFr6OIuljIWw/489PUAcrNzZYow1q8ru42SvOEmPRiiiqq2XE6i0/3X2B3fW/hlEsVTF53nOc2nFCCs/m/xPHJ3nPklFRRXadDpxNKTZ4QgqiLhXwSco6sEg21OsFX9U0jLlOrVXRzt+W9HfGsPXQl6NsZm81DKw7R9729jF0V3iDD4tzuGEMIwTvvvIObmxsWFhYMHTqU8+cNO94UFhYyadIkbGxssLOzY8qUKZSX3zgmuBMaPeD77rvveO2115g7dy7R0dEEBgYyfPhw8vLyrpv+yJEjTJgwgSlTpnDy5EnGjBnDmDFjOHPmTIPlsU4nSMovN7grScrTD/LrZG12TcCxMSKV+T/HGcwIcLXLd4hns0sJTdBfp04nSMwpM3gE+vb2WKW7/qVyLUM/DuORVUd4f2fCdY97Pk//YUq5VMHAD/YxbVP0Da9p1o+n2RBxkTd/OG2wvkJbq1xnTomGzcfTWX8klbisUj4LS2L3mRz2nM1Vxrg7n1tOdZ0OG3Nj/nlPRwLa2VJUWU1E/V2uo5UZE4M8OXKhgK/CU9hxOhtTYzUzh/vRyakNfi7W/HuEH5HJBfSYt4cv63+kXx/mp8x1aWthwuYXglkxoReLHgkAYHw/Tw6/eR9jev3xdGpSyzO4/gfP2syY4E5XhkaJNhi25cowLUEdHZTX1uYmjOnlgXV956Dx/TzZ+Gx/7up85zo4SbePtbkJT/T3ZNtLAzj+n6Ese6InTwd70cXNBiHgVHoxy0LP8+SXkfSYv4eBH+zjybWRvLUtls8PJhFyNpeEnFIq/qD37+ie7vT1ssfVxpwFo7tjpAbX+lk8Lv+5tjIzpqubDf289Z+1tYdSUKlUTOzvSXcPGzo7W1Fdq2Pv2Wt7++aWapi4Vt/soKZOp9QQDvJ1wt/Vmkd6exBdHzh29zCsDfz8YDK5pVqq63TsjsvhqS+PUa6tJbC9HW+O1Pcc33Yy85phbVxszPFzsUYICE+6hFdb/d/asER9W9aY9GIAZb8LeWX8GKV/kvVUsDdmxkZsPZnJP7/WB4Tv74zn0dURrK+/0QcY2f3aJhGmRmosTPTT6a3Yd4FybS2vfR/DmcxSLpVrKamqUYZaup6a2lsfIqchYowPPviA5cuXs2bNGiIjI2nTpg3Dhw9Ho7nSrGDSpEnExcUREhLCr7/+ysGDB3nhhRduOf+3U6N32vj44495/vnneeaZZwBYs2YNO3bs4KuvvmLWrFnXpF+2bBkjRozgjTfeAGDBggWEhITw6aefsmbNmr+VF51OXLe2KCa9iGnfnqR/BwelzUfSJX1w1fE6j053xeZgaqwms7iKjk5tUKtUBu2QVCoV7Rz0fzC2RmcwvJsrW6LSefPHWB7r044ljwdSpqnhh6gM/nc0jYVju1OhraWg/m7qx+gM5jzU5ZqGyhfy9G1CAjxsOZVRgrmx/pxnMkv47ng6c0d1xUitorSqVrnzPJVRQnFltVKbMX3TSbJLNCwY042LBZUGbZ/2JeQpbahA/+jB08GSjx4PpFxbq+Rnw5FU6nSCAA9bfnn5bgCC6h937IjNZm6Zhnb2loS+Plg5Vl6pRhngVq2CR3obBnJd3W3o6i4f3Uow434fzE3U3N3ZUfleaWrq2HJCX0Py7xF+rA9PJa9My0AfRxz/Ym9OqXlxsjZjdE8PZU7tnBINBxLziEguIOpiERlFVaQX6hcuXLu/vaUJ7R0scbM1x6GNKe62FrR3sMS+jSmzH/DHxtwEO0tT/i/Ym+cHdiQhp4yRyw4B8ECAGyqVince6sriXQm8OcIfC1MjFozpjqamjoikAsxNjDBSq3huwwkyiir59wg/enva4+lgibW5MWWaWs7nlivzew/2c+KZAR0AeHWzfmDvAA9bfjmVRfiFS8we2YUHA9wwMVKjVqnYHpNJYUU1e+JyGdevPXd3dsTH2QpLUyMKyrXYWpiQnF/OxcJKenvaM9DHkcTcMo4kFTCimyufH0xWeg5fDvgsTI0YsfQQibn635bLwau2to5PQs6RWVzFyfRitkZfadY00MeRjc/2v25HmsF+TkwI8uTzg8l8cTCZDo769tlWZsZ8PE7fbrKdvb5NbqmmhrjMUvp522NspEZTU8cT62Nu+XNxu2MMIQRLly7l7bffZvTo0QBs3LgRFxcXtm/fzhNPPEF8fDy7d+/m+PHj9O3bF4AVK1bwwAMPsGTJEtzdG6e5SKMGfNXV1URFRTF79mxlnVqtZujQoURERFx3n4iICF577TWDdcOHD2f79u3XTa/VatFqr7TlKCnR3/lnZxtOxfRFeBqbo7OZPawTgR42OFldeZzTFh3/CmrLf0OSiellh2MbE84mZ6HTVrLuQDyRienMG+lDO3tzIlOLOZqo/+GxF+V8GZLJJ/tSeKi7E6MDXSmpqsHB0pT+zipWayvZH3uR5FRnNh2KR6et5Psj53jEzxJ3WzOGdbRk++k83t8exc//7IPveB/+b+MpirRwPC4JN1szDicV4evcBhcbM2KT9Hn6h5sjJ5MqScqsIj4plXFroynX1PF9RCIOliY8N6A9Oq1+XLtpg7zIzsqizMyInNJqQk+nohOgKXEl6lw+Om0lZsZqtLU68rWV/HS8FJ1Wx3ujfHA1rqSiqJIgFwBjMjIyOJJcxL5Tmei0lQxs50RGhv7/wkmlw7hOQ3Wtjr7v/MLLg7yY0NfwQ9/BWpB0qYpJQR6UXMqlYWdIlZqzid2tAI3y+Rr2aSTlWv0NQ0fLGubf70F4chHjezsraaTW5253NXe7O8HdThRWVHOxQN9rO7O4ioyiKtIKK8gu0VBaVUt+KeTf5PjHKhWYm6gxVhmhVhtx8GQCEy+kYGFiRBsjNTM37MdUrcLSzAh7S1PMTYwwVqs4l1/BnlP6Ns+TP8vluWB3xgS64WmpI7akkh/CzzApwIERHUzxtteRkZGBTghikjLRaatwN9Py4c8JpBZq6G4vGOzblud662v91NpSPg9PZ0t4PDaU093NmiUPeuJgaYxKW0xGRjEbwtNYdzSTEV0dGd3DBaor6W6nw0FVgU5bSUpWJeeSU4lMSKv/LTFmaWIGuhr9d2tkJ1fl+9TZFtJzK9l5LJEtT/lxMKkQFXBPJwcyMw3btcfnlHM2u5xe7W14MsCK7w5XU1RaycrfTqHTVnKfnzNdbWtBlJGRUUa5tpanvz5NdomWRwJdeHmwF+/tTiIzT19RUVJSgo3NlQoAMzMzzMyuvbFriBgjJSWFnJwchg4dqmy3tbUlKCiIiIgInnjiCSIiIrCzs1OCPYChQ4eiVquJjIxk7NixN/5wNSTRiDIzMwUgjhw5YrD+jTfeEP3797/uPiYmJuLbb781WLdy5Urh7Ox83fRz584VgFzkIhe5yEUucmmBy9y5c+9YjBEeHi4AkZWVZZDm8ccfF+PGjRNCCLFw4ULh6+t7zbGdnJzEqlWrrnveO6HRH+k2tNmzZxtE67W1tcTHx9O+fXvU6kZvwthklZWV0bVrV86ePYu1tRz2pKmS5dQ8yHJqHmQ5NW06nY60tDS6du2KsfGV8OV6tXvStRo14HN0dMTIyIjcXMOGrLm5ubi6Xn8MLFdX11tKf72q3gEDBvyNXLcOpaX63rseHh4GVedS0yLLqXmQ5dQ8yHJq+jw9b36M2YaIMS7/m5ubi5ubm0Ganj17Kml+3ymktraWwsLCG573TmjUKi5TU1P69OlDaGiosk6n0xEaGkpwcPB19wkODjZIDxASEnLD9JIkSZIktT4NEWN06NABV1dXgzSlpaVERkYqaYKDgykuLiYq6sog2fv27UOn0xEUFHTbru+WNdrD5HqbN28WZmZmYv369eLs2bPihRdeEHZ2diInJ0cIIcRTTz0lZs2apaQPDw8XxsbGYsmSJSI+Pl7MnTtXmJiYiNjY2Ma6hBappKREAKKkpKSxsyL9AVlOzYMsp+ZBllPL0xAxxuLFi4WdnZ346aefxOnTp8Xo0aNFhw4dRFVVlZJmxIgRolevXiIyMlIcPnxY+Pj4iAkTJty5C7+ORg/4hBBixYoVwtPTU5iamor+/fuLo0ePKtsGDRoknn76aYP033//vfD19RWmpqaiW7duYseOHXc4xy2fRqMRc+fOFRqNprGzIv0BWU7Ngyyn5kGWU8t0u2MMnU4n5syZI1xcXISZmZkYMmSISExMNEhTUFAgJkyYIKysrISNjY145plnRFlZWYNd481QCSHnOJEkSZIkSWrJZDdVSZIkSZKkFk4GfJIkSZIkSS2cDPgkSZIkSZJaOBnwSZIkSZIktXAy4GtFFi1aRL9+/bC2tsbZ2ZkxY8aQmJhokEaj0TB16lTatm2LlZUVjz766DWDUKalpfHggw9iaWmJs7Mzb7zxBrW1tXfyUlq01atX06NHD2xsbLCxsSE4OJhdu3Yp22UZNU2LFy9GpVLx6quvKutkWTW+efPmoVKpDBZ/f39luywjqbWQAV8rEhYWxtSpUzl69CghISHU1NQwbNgwKioqlDQzZszgl19+YcuWLYSFhZGVlcUjjzyibK+rq+PBBx+kurqaI0eOsGHDBtavX88777zTGJfUIrVr147FixcTFRXFiRMnuO+++xg9ejRxcXGALKOm6Pjx43z22Wf06NHDYL0sq6ahW7duZGdnK8vhw4eVbbKMpFajUQeFkRpVXl6eAERYWJgQQoji4mJhYmIitmzZoqSJj48XgIiIiBBCCLFz506hVquVQSuFEGL16tXCxsZGaLXaO3sBrYi9vb1Yu3atLKMmqKysTPj4+IiQkBAxaNAg8corrwgh5PepqZg7d64IDAy87jZZRlJrImv4WrGSkhIAHBwcAIiKiqKmpoahQ4cqafz9/fH09CQiIgKAiIgIAgICcHFxUdIMHz6c0tJSpQZKun3q6urYvHkzFRUVBAcHyzJqgqZOncqDDz5oUCYgv09Nyfnz53F3d6djx45MmjSJtLQ0QJaR1LoYN3YGpMah0+l49dVXGTBgAN27dwcgJycHU1NT7OzsDNK6uLiQk5OjpLn6D9/l7Ze3SbdHbGwswcHBaDQarKys2LZtG127diUmJkaWUROyefNmoqOjOX78+DXb5PepaQgKCmL9+vX4+fmRnZ3N/PnzGThwIGfOnJFlJLUqMuBrpaZOncqZM2cM2rJITYefnx8xMTGUlJTwww8/8PTTTxMWFtbY2ZKukp6eziuvvEJISAjm5uaNnR3pBkaOHKm87tGjB0FBQXh5efH9999jYWHRiDmTpDtLPtJthaZNm8avv/7K/v37adeunbLe1dWV6upqiouLDdLn5ubi6uqqpPl9D7bL7y+nkf4+U1NTOnfuTJ8+fVi0aBGBgYEsW7ZMllETEhUVRV5eHr1798bY2BhjY2PCwsJYvnw5xsbGuLi4yLJqguzs7PD19eXChQvy+yS1KjLga0WEEEybNo1t27axb98+OnToYLC9T58+mJiYEBoaqqxLTEwkLS2N4OBgAIKDg4mNjSUvL09JExISgo2NDV27dr0zF9IK6XQ6tFqtLKMmZMiQIcTGxhITE6Msffv2ZdKkScprWVZNT3l5OUlJSbi5ucnvk9S6NHavEenO+de//iVsbW3FgQMHRHZ2trJUVlYqaV588UXh6ekp9u3bJ06cOCGCg4NFcHCwsr22tlZ0795dDBs2TMTExIjdu3cLJycnMXv27Ma4pBZp1qxZIiwsTKSkpIjTp0+LWbNmCZVKJfbs2SOEkGXUlF3dS1cIWVZNweuvvy4OHDggUlJSRHh4uBg6dKhwdHQUeXl5QghZRlLrIQO+VgS47rJu3TolTVVVlXjppZeEvb29sLS0FGPHjhXZ2dkGx0lNTRUjR44UFhYWwtHRUbz++uuipqbmDl9Ny/Xss88KLy8vYWpqKpycnMSQIUOUYE8IWUZN2e8DPllWjW/8+PHCzc1NmJqaCg8PDzF+/Hhx4cIFZbssI6m1UAkhRGPWMEqSJEmSJEkNS7bhkyRJkiRJauFkwCdJkiRJktTCyYBPkiRJkiSphZMBnyRJkiRJUgsnAz5JkiRJkqQWTgZ8kiRJkiRJLZwM+CRJkiRJklo4GfBJkiRJkiS1cDLgkyRJaiBz5szhhRdeuG3Hq66uxtvbmxMnTty2Y0qS1DrIgE+SWjiVSvWHy7x58xo7i7edt7c3S5cubdQ85OTksGzZMv7zn/8o6yoqKnjiiSdwc3NjwoQJVFZWXrPPyy+/TMeOHTEzM6N9+/aMGjWK0NBQAExNTZk5cyZvvvnmHb0WSZKaPxnwSVILl52drSxLly7FxsbGYN3MmTMbO4s3RQhBbW3tHT1ndXX1X9537dq13HXXXXh5eSnrli5dipWVFXv27MHCwsIgKE1NTaVPnz7s27ePDz/8kNjYWHbv3s29997L1KlTlXSTJk3i8OHDxMXF/eW8SZLU+siAT5JaOFdXV2WxtbVFpVIZrNu8eTNdunTB3Nwcf39/Vq1apeybmpqKSqXi+++/Z+DAgVhYWNCvXz/OnTvH8ePH6du3L1ZWVowcOZL8/Hxlv8mTJzNmzBjmz5+Pk5MTNjY2vPjiiwYBlE6nY9GiRXTo0AELCwsCAwP54YcflO0HDhxApVKxa9cu+vTpg5mZGYcPHyYpKYnRo0fj4uKClZUV/fr1Y+/evcp+gwcP5uLFi8yYMUOpxQSYN28ePXv2NPi/Wbp0Kd7e3tfke+HChbi7u+Pn5wdAeno648aNw87ODgcHB0aPHk1qauof/r9v3ryZUaNGGawrKirC19eXgIAA/P39KS4uVra99NJLqFQqjh07xqOPPoqvry/dunXjtdde4+jRo0o6e3t7BgwYwObNm//w/JIkSVeTAZ8ktWLffPMN77zzDgsXLiQ+Pp7333+fOXPmsGHDBoN0c+fO5e233yY6OhpjY2MmTpzIv//9b5YtW8ahQ4e4cOEC77zzjsE+oaGhxMfHc+DAATZt2sTWrVuZP3++sn3RokVs3LiRNWvWEBcXx4wZM3jyyScJCwszOM6sWbNYvHgx8fHx9OjRg/Lych544AFCQ0M5efIkI0aMYNSoUaSlpQGwdetW2rVrx7vvvqvUYt6K0NBQEhMTCQkJ4ddff6Wmpobhw4djbW3NoUOHCA8Px8rKihEjRtywBrCwsJCzZ8/St29fg/XTpk3js88+w8TEhHXr1vHKK68o6Xfv3s3UqVNp06bNNcezs7MzeN+/f38OHTp0S9clSVIrJyRJajXWrVsnbG1tlfedOnUS3377rUGaBQsWiODgYCGEECkpKQIQa9euVbZv2rRJACI0NFRZt2jRIuHn56e8f/rpp4WDg4OoqKhQ1q1evVpYWVmJuro6odFohKWlpThy5IjBuadMmSImTJgghBBi//79AhDbt2//0+vq1q2bWLFihfLey8tLfPLJJwZp5s6dKwIDAw3WffLJJ8LLy8sg3y4uLkKr1Srrvv76a+Hn5yd0Op2yTqvVCgsLC/Hbb79dNz8nT54UgEhLS7tmW11dncjOzjY4XmRkpADE1q1b//RahRBi2bJlwtvb+6bSSpIkCSGEcaNGm5IkNZqKigqSkpKYMmUKzz//vLK+trYWW1tbg7Q9evRQXru4uAAQEBBgsC4vL89gn8DAQCwtLZX3wcHBlJeXk56eTnl5OZWVldx///0G+1RXV9OrVy+Ddb+vJSsvL2fevHns2LGD7OxsamtrqaqqUmr4/q6AgABMTU2V96dOneLChQtYW1sbpNNoNCQlJV33GFVVVQCYm5tfs02tVuPq6mqwTghxS3m0sLC4psOHJEnSH5EBnyS1UuXl5QB88cUXBAUFGWwzMjIyeG9iYqK8vtwm7vfrdDrdLZ97x44deHh4GGwzMzMzeP/7R5wzZ84kJCSEJUuW0LlzZywsLHjsscf+tIOFWq2+JrCqqam5Jt3vz1deXk6fPn345ptvrknr5OR03XM5OjoC+jZ7N0pzNR8fH1QqFQkJCX+aFvSPgG/muJIkSZfJgE+SWikXFxfc3d1JTk5m0qRJt/34p06doqqqCgsLCwCOHj2KlZUV7du3x8HBATMzM9LS0hg0aNAtHTc8PJzJkyczduxYQB+Q/b4DhampKXV1dQbrnJycyMnJQQihBK0xMTF/er7evXvz3Xff4ezsjI2NzU3lsVOnTtjY2HD27Fl8fX3/NL2DgwPDhw9n5cqVTJ8+/Zqgs7i42KAd35kzZ66pCZUkSfojstOGJLVi8+fPZ9GiRSxfvpxz584RGxvLunXr+Pjjj//2saurq5kyZQpnz55l586dzJ07l2nTpqFWq7G2tmbmzJnMmDGDDRs2kJSURHR0NCtWrLimw8jv+fj4sHXrVmJiYjh16hQTJ068pnbR29ubgwcPkpmZyaVLlwB97938/Hw++OADkpKSWLlyJbt27frT65g0aRKOjo6MHj2aQ4cOkZKSwoEDB5g+fToZGRnX3UetVjN06FAOHz58k/9bsHLlSurq6ujfvz8//vgj58+fJz4+nuXLlxMcHGyQ9tChQwwbNuymjy1JkiQDPklqxZ577jnWrl3LunXrCAgIYNCgQaxfv54OHTr87WMPGTIEHx8f7rnnHsaPH8/DDz9sMMjzggULmDNnDosWLaJLly6MGDGCHTt2/Om5P/74Y+zt7bnrrrsYNWoUw4cPp3fv3gZp3n33XVJTU+nUqZPy6LNLly6sWrWKlStXEhgYyLFjx25qDEJLS0sOHjyIp6cnjzzyCF26dGHKlCloNJo/rPF77rnn2Lx5800/6u7YsSPR0dHce++9vP7663Tv3p3777+f0NBQVq9eraSLiIigpKSExx577KaOK0mSBKASt9paWJIk6U9MnjyZ4uJitm/f3thZaTRCCIKCgpgxYwYTJky4bccdP348gYGBvPXWW7ftmJIktXyyhk+SJKkBqFQqPv/889s6O0h1dTUBAQHMmDHjth1TkqTWQdbwSZJ028kaPkmSpKZFBnySJEmSJEktnHykK0mSJEmS1MLJgE+SJEmSJKmFkwGfJEmSJElSCycDPkmSJEmSpBZOBnySJEmSJEktnAz4JEmSJEmSWjgZ8EmSJEmSJLVwMuCTJEmSJElq4f4fzglUE6wAZVEAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import pyTGA as tga\n", "import matplotlib.pyplot as plt\n", "import pandas as pd\n", "import numpy as np\n", "from scipy.special import erf\n", "import os\n", "#laoding the data\n", "data_dir = os.path.abspath(os.path.join(os.getcwd(), '..', '..', '..', 'example_data'))\n", "pvc_exp = tga.parse_txt(os.path.join(data_dir, 'PVC.txt'), exp_type='general')\n", "pvc_exp = tga.calc_DTGA_stage(pvc_exp, stage_name='stage4', x='Temp', y='relative')\n", "data = pvc_exp.stages['stage4']\n", "data['norm. weight'] = data['Unsubtracted weight']/data['Unsubtracted weight'].max()\n", "\n", "#plotting\n", "fig, ax = plt.subplots()\n", "ax2 = ax.twinx()\n", "ax.plot(data['Sample Temp.'], data['norm. weight'], label='Norm. Weight')\n", "ax2.plot(data['Sample Temp.'], data['DTGA_Temp_relative'], label='DTGA', linestyle='--')\n", "ax.set_xlabel('Temperature (°C)')\n", "ax.set_ylabel('Norm. Weight (a.u.)')\n", "ax2.set_ylabel('DTGA (1/°C)')\n", "\n", "ax.set_ylim(0,)\n", "ax2.set_ylim(0,)\n", "ax.set_xlim(data['Sample Temp.'].min(), data['Sample Temp.'].max())\n", "fig.legend(loc='upper right', bbox_to_anchor=(1, 1), bbox_transform=ax.transAxes, frameon=False)\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "bea1878b", "metadata": {}, "source": [ "To fit the thermogram, we first set up a model. For this we are going to use the sum of two modified error functions with a step location ($x_0$) and a width ($\\sigma$) (Eq.1): \n", "\n", "\\begin{equation}\n", "0.5 \\cdot \\left(1 - \\text{erf}\\left(\\frac{x - x_0}{\\sigma \\sqrt{2}}\\right)\\right)\n", "\\end{equation}" ] }, { "cell_type": "code", "execution_count": 2, "id": "d65c5790", "metadata": {}, "outputs": [], "source": [ "from scipy.special import erf\n", "\n", "def error_function_desc_full(x, center=0, sigma=1):\n", " \"\"\" Error function for a descending curve form 1 to 0 with center and sigma (width).\n", " \"\"\"\n", " return 0.5 * (1 - erf((x - center) / (sigma * np.sqrt(2))))\n", "\n", "def double_error_function(T, A1, T1, S1, A2, T2, S2, R):\n", " \"\"\" Double error function for a descending curve with two steps and a residual.\n", " \n", " parameters:\n", " ----------\n", " T: Temperature array\n", " A1: Amplitude of the first error function\n", " T1: Center of the first error function\n", " S1: Sigma (width) of the first error function\n", " A2: Amplitude of the second error function\n", " T2: Center of the second error function\n", " S2: Sigma (width) of the second error function\n", " R: Residual value\n", " \"\"\"\n", " return A1 * error_function_desc_full(T, center=T1, sigma=S1) + A2 * error_function_desc_full(T, center=T2, sigma=S2) + R" ] }, { "cell_type": "markdown", "id": "0db171e8", "metadata": {}, "source": [ "For fitting:" ] }, { "cell_type": "code", "execution_count": null, "id": "246a71fa", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "fitted parameters:\n", "A1 = 0.69\n", "T1 = 302.73\n", "S1 = 19.84\n", "A2 = 0.32\n", "T2 = 439.18\n", "S2 = 25.33\n", "R = -0.00\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAGwCAYAAABVdURTAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAbRJJREFUeJzt3Xd4FOXexvHv7qZ3ICSBEAi9SO+hdxRFAVFUlKJHD4gFkFdFBQQLNgRFmg2wggVsIIqRAAJK70iH0EJCSe+78/4RWc2hJZgwKffnuuYiO/vM7G8ybPbeZ56ZsRiGYSAiIiJSQljNLkBERESkICnciIiISImicCMiIiIlisKNiIiIlCgKNyIiIlKiKNyIiIhIiaJwIyIiIiWKi9kFXG8Oh4OTJ0/i6+uLxWIxuxwRERHJA8MwSEpKomLFilitV+6bKXXh5uTJk4SFhZldhoiIiFyDY8eOUalSpSu2KXXhxtfXF8j55fj5+ZlcjYiIiORFYmIiYWFhzs/xKyl14ebCoSg/Pz+FGxERkWImL0NKNKBYREREShSFGxERESlRFG5ERESkRFG4ERERkRJF4UZERERKFIUbERERKVEUbkRERKREUbgRERGREkXhRkREREoUhRsREREpUUwNN6tWraJ3795UrFgRi8XCN998c9VloqKiaNq0Ke7u7tSoUYN58+YVep0iIiJSfJgablJSUmjUqBEzZszIU/vDhw9z880307lzZ7Zu3crIkSP5z3/+w08//VTIlYqIiJQeQ4YMoU+fPvlaJjw8nGnTphVKPfllari56aabePHFF+nbt2+e2s+ePZuqVasyZcoU6tatyyOPPEL//v2ZOnVqIVdaMhmGgd1hkG13kJntICPbTnqWnbRMOykZ2SRnZJOYnkVCWhbxqZmcT8nkbHIGZ5IziE1K53RiOnFJGaRl2jEMw+zNERG5ZkOGDMFisfDKK6/kmv/NN9/k6UaNZmndujXDhg3LNW/27NlYLJaLjmwMGTKE9u3b52m9b731VoEfGTly5AgWi4WtW7cW6HovpVjdFXzdunV069Yt17yePXsycuTIyy6TkZFBRkaG83FiYmLOD1OmgIeHc74DsDsMHFhIGfkEWXYHCWlZGIsW43LoANkOA4fDINthkGk3MAwDTzcbuwc8wLmMnHAQ9nsUZY4ewAEYBjgMI9e/v3W/gzQXdwzDoMb23wmJ3o/xV1vDMDD+qgMDVrW/lVQPLxwG1Nq7mbCj+3Je/0IgcUC2w0G23WDxDZ044+qNwzDocP4QzeIOggEOAxyGA7thwe7Iqfm7Wm057uqDYcANMQdofGrfZX93P9dsTZxPWQBqxR2h2Yk/cz3/zzjzW40WpAVXwN3FSrVzx2kcvQub1YKL1YqbzYK7mwsuVgtWCxxq0IrkCpVwtVkJT4ihzv6teLracLFZcLVZcXWx4eNuI8DLDY/OHaF69ZwXOXYMoqL+flFXV/D0zNmPnp5QqxaEhOQ8Z7fn/GJditV/cRExkYeHB6+++ir//e9/KVOmTIGtNzMzEzc3twJb3z917tyZxYsX55q3YsUKwsLCiIqKYsiQIc75UVFRDB48OE/r9ff3L8gyr7ti9Zc/JiaG4ODgXPOCg4NJTEwkLS0NT0/Pi5aZPHkyEydOvHhlkyblemj9a8q2WGma1NA5f/bi2dy4b91la+pra06miysAU7//jI67oy7b9im/ZiR4+gLw8rJv6bdt2WXbTgtoRIxfIAAdfl3G4I3fXrbt0goNyCiXs+2Nt/7Gw+sWXrZt1ODqGCE1AOh4eDNPrvrosm33lq/iDDdtjm7n+ch3L9t2kN9EVvmUA6D9zk2MWTb9sm3/2+cZfqrdBoBbd6/kge9fv2zbJ2/7P9a1uRE/D1c6/LmOp2Y9ddm2zJwJw4fn/LxuHbRvDwEBUL48BAZCaChUq5YTljp2hNq1L78uESkQhmGQlmU35bU9XW356nXp1q0bBw4cYPLkybz22muXbff1118zfvx4Dhw4QIUKFXj00Ud54oknnM+Hh4fzwAMPsH//fr755hv69etHp06dGDlyJJ988glPPPEEx44do1evXnz00Ud8+eWXTJgwgYSEBO677z6mTp2KzWbLU82dO3fmlVdeISYmhpC/vtytXLmS8ePH59qGw4cPc/ToUTp37gzAsWPHeOKJJ/j555+xWq20b9+et956i/DwcCCnlyc+Pt45FjYpKYlhw4bxzTff4Ofnx5NPPsm3335L48aNcx2KSk1N5f777+fLL7+kTJkyPPfcczz00EMAVK1aFYAmTZoA0LFjR6L++YW1ABWrcHMtxo4dy+jRo52PExMTCQsL46sbuuBly9l8yz8OqRh/vRGsFvD3dGV/vebYygRgs+T0OtgsYLNasGAh0+6gfe0g/P28cHe1Yo1vx9YgXyyAxWLBgvHXzzmPH+5RB3x8sVosVDI6cjAoJ8lbLBYshsGF96AFeLpfE7LLlMVqgcrenTleLqcm6191XPjZZrXwwaNdcA0LxWKxkB0WQ6x/2t+v+9e2WLBgscDsR7riXrsGVosF9+9SyAxIhL/aYRjOPwQW4NMne0G9uliwYF1ix/A+7Szwwu/MMMBuGEwZ0Y2zteuTlW3gviKdc1n7MBw5vVb2v3qYDHJ6sG7s1pgGdWuRke3A3yeWHdFtsTsMZ++S3WGQme0g2+HgmGcAx86lAWm4J1lZWbWpc5+5OrJxz87EIzsT9+wsZq44wb7s36hczotue7fSFyA+Pmfavz/3f4xp0/4ON8eOweLFcNNNULPmv/nvJiL/Iy3LTr3x5oyL3D2pJ15uef+Ys9lsvPzyy9xzzz089thjVKpU6aI2mzZt4s477+T5559nwIABrF27locffphy5crl6iV54403GD9+PBMmTABg9erVpKam8vbbb7NgwQKSkpLo168fffv2JSAggKVLl3Lo0CFuv/122rZty4ABA/JUc9u2bXF1dWXFihXcfffd7N69m7S0NB544AGeeuopDh8+TNWqVVmxYgUeHh5ERESQlZVFz549iYiIYPXq1bi4uPDiiy9y4403sn379kv2Mo0ePZo1a9bw3XffERwczPjx49m8eTONGzfO1W7KlCm88MILPPPMM3z11VcMHz6cjh07Urt2bdavX0/Lli355ZdfuOGGGwqtNwuKWbgJCQnh9OnTueadPn0aPz+/S/baALi7u+Pu7n7RfMv0tykfEkgFf08Cfdz+OhxixcVq4Q6bFZv1QtrvccWaOv7zQb+Gl2sGQL1/Pujw9BXbVvnng6aPA49ftm2uzsMRD+RMeWl7z50502Xk+t5w26050/+wkPOfqPxfEwD33Z4zXUblfz7oURueG3rJdikZ2VRNyuBsSgaJ6dmkZDQh+uH+xCVnkpiWxdmUTI6dS+X4+VTOJGfmLHQigR0nEvjRUZkXHv2UMmmJlElLJDg9iYaOBOqnnyE8IYa0ijUIzbTj6WaDVavg8cdzpho14K67YMiQvw+HiUip0bdvXxo3bsyECRP44IMPLnr+zTffpGvXrowbNw6AWrVqsXv3bl5//fVc4aZLly65enNWr15NVlYWs2bNovpff1v69+/Pxx9/zOnTp/Hx8aFevXp07tyZFStW5DnceHt707JlS6Kiorj77ruJioqiXbt2uLu706ZNG6KioqhatSpRUVFERETg7u7OJ598gsPh4P3333d+oZ07dy4BAQFERUXRo0fuz72kpCTmz5/PZ599RteuXZ3tK1aseFE9vXr14uGHHwbgqaeeYurUqaxYsYLatWtTvnzOp0S5cuWcvUyFpViFm4iICJYuXZpr3vLly4mIiMj3um5vFoafn19BlSaFwNvdBW93F8IDva/aNjUzm2Pn0jh6NoXoc6kcOZvC0bOpHD2bypb4NOwOgyX/XGATWDcvIzzQmz5nEunfNIIKOzZiOXAAXnwxZ+raFZ55Bjp3hiI8oFCkKPN0tbF7Uk/TXvtavPrqq3Tp0oUxY8Zc9NyePXu47bbbcs1r27Yt06ZNw263Ow8nNW/e/KJlvby8nMEGcoZVhIeH4+Pjk2tebGxsvurt1KkTX375JZAzrqZTp07A34d9hg4dSlRUFA8++CAA27Zt48CBA/j6+uZaT3p6OgcPHrxo/YcOHSIrK4uWLVs65/n7+1P7Eof2Gzb8+0u+xWIhJCQk39tTEEwNN8nJyRw4cMD5+PDhw2zdupWyZctSuXJlxo4dy4kTJ/joo5xxIcOGDeOdd97hySef5P777+fXX3/liy++YMmSJZd7CSklvNxcqB3iS+0Q34uey7I7OBWfzu5Tiew4Ec/24wnsOZXImeRMDsWl8CaVebP7s3h3SOW2E1u5589fuWHnH1giIyEyMueQVo0aJmyVSPFnsVjydWioKOjQoQM9e/Zk7NixuXpj8sPb++IvZa6urrkeWyyWS85zOBz5eq3OnTvz0ksvceLECaKiopyhrGPHjsyZM4eDBw9y7NgxunTpAuR89jZr1oxPP/30onVd6F25VgWxPQXB1P9xGzdudA5uApxjYwYPHsy8efM4deoU0dHRzuerVq3KkiVLGDVqFG+99RaVKlXi/fffp2dPc74VSPHgarNSuZwXlct5cWP9v7tCY5PS+fNUEjtPJrDu4FnWHz7HZ9Xa8Fm1NoS2jeW/G7+hjlsWGYY/7S6MRzIM9eKIlAKvvPIKjRs3vqh3om7duqxZsybXvDVr1lCrVq08DwIuaG3atMHNzY2ZM2eSnp5Os2bNAGjRogVxcXF8+OGHzsNXAE2bNmXhwoUEBQXl6QhGtWrVcHV1ZcOGDVSunDOwICEhgX379tGhQ4c813lhjI3dXvgDzE0NN506dbri9VEudY59p06d2LJlSyFWJaVFkK8HQb4edKhVnoc71SA9y87OEzlBZ/FWb8b7P5QTZj5YT50QXx6t7spNzw7DOnMGtGtndvkiUogaNGjAwIEDefvtt3PNf+KJJ2jRogUvvPACAwYMYN26dbzzzjvMnDmzUOr43yMYl+Lp6Unr1q2ZPn06bdu2dYYsNze3XPMv9KoMHDiQ119/ndtuu41JkyZRqVIljh49yqJFi3jyyScvGkjt6+vL4MGD+b//+z/Kli1LUFAQEyZMwGq15utstKCgIDw9PVm2bBmVKlXCw8Oj0E45172lRP7i4WqjeXhZHu1ak8jRHfnh0XYMbhOOl5uNP2OSyBj7LNYd23F07IhjypSc4CMiJdakSZMuOqTStGlTvvjiCxYsWED9+vUZP348kyZNuubDV1fzv0cwLqdz584kJSU5x9tc0LFjR5KSknIdJfHy8mLVqlVUrlyZfv36UbduXR544AHS09Mv25Pz5ptvEhERwS233EK3bt1o27YtdevWxeMf14u7GhcXF95++23mzJlDxYoVLxq7VJAsRim7tGxiYiL+/v4kJCRoQLHkSUJqFp9viOaLX3czYtE0bt/5KwCx995P0Nw5ulCgiJQ6KSkphIaGMmXKFB544PJn6Bak/Hx+q+dG5Cr8vVwZ1rE6y8bdzPkZ7/J69wdxYCHokw/Z2PFWTp1LNrtEEZFCtWXLFj7//HMOHjzI5s2bGThwIECh9r78Gwo3Innk5mLlPx2qM3TRdD4f8zpZVhvN1/7Ihk63sXTrcbPLExEpVG+88QaNGjWiW7dupKSksHr1agIDA80u65J0WErkGkW/9zGhw4ayv2wl7rn7Ze7r3ZzHu9bEatXZVCIiBS0/n98aLCByjSo/eB/28Ap8d96Pc5vjeCtyP3tjkphyZyO83fXWEhExiw5LifwLtu7dePLOlrzevyFuNisrtkVz7wd/kJyRbXZpIiKllsKNSAG4o2kov2av49cPhnN4bzTDP9lElv36X5VTREQUbkQKRloalb5ZSGj8aaYse5vV++J4bvHOK16kUkRECofCjUhB8PaGhQvBzY2u+37nrh3LWbjxGDNWHLj6siIiUqAUbkQKSuPG8NJLAEz8bR7lk8/zxs/7WPHn9b8jrogUjE6dOjFy5Eizy5B8UrgRKUgjR0LTprgnJTJ3y8cAjPlyG7FJ6ebWJSJXNGTIECwWy0XTa6+9xgsvvOBsFx4ezrRp08wrVPJE4UakILm4wPvvg81G/bU/c0/8n5xNyWTcNzvNrkxEruLGG2/k1KlTuaZmzZrh6+trdmmSTwo3IgWtSRN49FHw8OCxqlZcrBZ+2nWaZTtPmV2ZiFyBu7s7ISEhuaauXbs6D0t16tSJo0ePMmrUKGfPjhRNCjcihWH8eNi7l5CnRjGsY3UAXlq6h4xsu8mFiZgkJeXyU3p63tumpeWtbSFYtGgRlSpVYtKkSc6eHSmaFG5ECkOZMlC5MgAPd65OkK87x86l8cnv0SYXJmISH5/LT7ffnrttUNDl2950U+624eGXbncNfvjhB3x8fJzTHXfckev5smXLYrPZ8PX1dfbsSNGkcCNSyLw2bWCqdT8A7/y6n4S0LJMrEpFL6dy5M1u3bnVOb7/9ttklyTXSDXBECtPPP0PPnrQpV46Goz9me2IWc9ccZmS3WmZXJnJ9JSdf/jmbLffj2CtcPsH6P9/Jjxy55pL+l7e3NzVq1Ciw9Yl5FG5EClOXLlCzJpb9+3n19Gpu8mzLJ78fZVjH6ni42q6+vEhJ4e1tftsC4Obmht2usXNFnQ5LiRQmFxcYNw6AOp++S3UPgzPJmXy37aTJhYnItQgPD2fVqlWcOHGCM2fOmF2OXIbCjUhhu/tuqFEDy9mzPH9+AwAf/nZY950SKYYmTZrEkSNHqF69OuXLlze7HLkMi1HK/sImJibi7+9PQkICfn5+ZpcjpcXMmTBiBPZq1Wl411uk2OGz/7SiTY1AsysTESkW8vP5rZ4bketh8GAoUwbboYM868i5meYHvx02uSgRkZJJ4UbkevD2hmHDIDycrrVzurIj/4zlUNwVziAREZFronAjcr08+yzs30/wkHvoUicIgK83Hze5KBGRkkfhRuR68fbOOXsK6NMkFIAl209pYLGISAFTuBG53jIy6LF9BeUc6Rw5m8quk4lmVyQiUqIo3Ihcb1274jHoXkaf2wLA99t1zRsRkYKkcCNyvfXtC8DNv/8A6NCUiEhBU7gRud4GDwY3NwJ2baPZ2cMcP5/G1mPxZlclIlJiKNyIXG+BgdCnDwCPHF8LwLJdMSYWJCJSsijciJhh0CAAItYvx+aws3JvnMkFiUhhOXLkCBaLha1bt162TVRUFBaLhfj4+AJ9bYvFwjfffFOg6ywOFG5EzNCjBwQG4nHuDO2ObuXPmCROJ6abXZVIqTVkyBAsFgsWiwVXV1eqVq3Kk08+SXr6v39fhoWFcerUKerXr18AlUpeKNyImMHVFe66C4BeCYcAWLlPvTciZrrxxhs5deoUhw4dYurUqcyZM4cJEyb86/XabDZCQkJw+es6V1L4FG5EzDJmDOzdy4nRYwFYpXAjYip3d3dCQkIICwujT58+dOvWjeXLlwPgcDiYPHkyVatWxdPTk0aNGvHVV185lz1//jwDBw6kfPnyeHp6UrNmTebOnQtc+rDU0qVLqVWrFp6ennTu3JkjR47kquX555+ncePGueZNmzaN8PBw5+MNGzbQvXt3AgMD8ff3p2PHjmzevLlAfyfFlWKkiFmqVAGgo/s53v71AKv3n8HuMLBZLSYXJlJwDMMgNSvVlNf2cvXCYrm299POnTtZu3YtVf56n06ePJlPPvmE2bNnU7NmTVatWsW9995L+fLl6dixI+PGjWP37t38+OOPBAYGcuDAAdLS0i657mPHjtGvXz9GjBjBQw89xMaNG3niiSfyXWNSUhKDBw9m+vTpGIbBlClT6NWrF/v378fX1/eatrukULgRMVmjSgGUdTU4l5bFtuPxNK1cxuySRApMalYqPpN9THnt5LHJeLt557n9Dz/8gI+PD9nZ2WRkZGC1WnnnnXfIyMjg5Zdf5pdffiEiIgKAatWq8dtvvzFnzhw6duxIdHQ0TZo0oXnz5gC5elj+16xZs6hevTpTpkwBoHbt2uzYsYNXX301X9vXpUuXXI/fffddAgICWLlyJbfccku+1lXSKNyImCktDZd77mHNT8tp9dAHrNwbp3AjYpLOnTsza9YsUlJSmDp1Ki4uLtx+++3s2rWL1NRUunfvnqt9ZmYmTZo0AWD48OHcfvvtbN68mR49etCnTx/atGlzydfZs2cPrVq1yjXvQmjKj9OnT/Pcc88RFRVFbGwsdrud1NRUoqOj872ukkbhRsRMnp6wfz+eaSl0O/AHK/dVYlT3WmZXJVJgvFy9SB6bbNpr54e3tzc1atQA4MMPP6RRo0Z88MEHzrOclixZQmhoaK5l3N3dAbjppps4evQoS5cuZfny5XTt2pURI0bwxhtvXFPtVqv1oiuXZ2Vl5Xo8ePBgzp49y1tvvUWVKlVwd3cnIiKCzMzMa3rNkkThRsRs/fvDrl3ctHctDzXoyvmUTMp4u5ldlUiBsFgs+To0VFRYrVaeeeYZRo8ezb59+3B3dyc6OpqOHTtedpny5cszePBgBg8eTPv27fm///u/S4abunXr8t133+Wa9/vvv1+0rpiYGAzDcI4b+t/r5KxZs4aZM2fSq1cvIGcsz5kzZ65lc0scnS0lYrbbbweg45EteKWnsuag/jiJFAV33HEHNpuNOXPmMGbMGEaNGsX8+fM5ePAgmzdvZvr06cyfPx+A8ePH8+2333LgwAF27drFDz/8QN26dS+53mHDhrF//37+7//+j7179/LZZ58xb968XG06depEXFwcr732GgcPHmTGjBn8+OOPudrUrFmTjz/+mD179vDHH38wcOBAPD09C+V3Udwo3IiYrX59qFkTt+xMuhzcwPrD58yuSEQAFxcXHnnkEV577TXGjh3LuHHjmDx5MnXr1uXGG29kyZIlVK1aFQA3NzfGjh1Lw4YN6dChAzabjQULFlxyvZUrV+brr7/mm2++oVGjRsyePZuXX345V5u6desyc+ZMZsyYQaNGjVi/fj1jxozJ1eaDDz7g/PnzNG3alPvuu4/HHnuMoKCgwvllFDMWo5TdjjgxMRF/f38SEhLw8/MzuxyRHM88A5Mns6R2W6YPe5llIzuYXZGISJGSn89v9dyIFAV/HZrqfGgjR4/HEZ+qAYEiItdKA4pFioKmTeGOO/gorRw2u4ONR87TrV6w2VWJiBRL6rkRKQosFvjiC44MGU6yuxcbjmrcjYjItVK4ESlCmlXJuYDfpiPnTa5ERKT4UrgRKUJa+jrou/NX4rfvIT3LbnY5IiLFksKNSBFS+emRTF3yJjftXMGukwlmlyMiUiwp3IgUIZa/bnbX9eAGth5TuBERuRYKNyJFyV+XUW98ah+Hdh40uRgRkeJJ4UakKKlQgaT6jQHw+/Vnc2sRESmmFG5EihiX23oD0GDbGl3MT0TkGijciBQxnr1vBqDd0W1sP3LW5GpERIofhRuRoqZ5c1K8/fDLSOHEqvVmVyMiUuyYHm5mzJhBeHg4Hh4etGrVivXrr/zHfNq0adSuXRtPT0/CwsIYNWoU6enp16lakevAZiPqlXdpMeJjIj1Dza5GRKTYMTXcLFy4kNGjRzNhwgQ2b95Mo0aN6NmzJ7GxsZds/9lnn/H0008zYcIE9uzZwwcffMDChQt55plnrnPlIoUr+JbuxPmUYdvxeAzDMLscEZFixdRw8+abb/Lggw8ydOhQ6tWrx+zZs/Hy8uLDDz+8ZPu1a9fStm1b7rnnHsLDw+nRowd33333VXt7RIqbGyr6Y7NaiEvKICZRPZMiIvlhWrjJzMxk06ZNdOvW7e9irFa6devGunXrLrlMmzZt2LRpkzPMHDp0iKVLl9Lrr2uDXEpGRgaJiYm5JpGiztPNxsNHf+OTBc9y8qOFZpcjIlKsmBZuzpw5g91uJzg4ONf84OBgYmJiLrnMPffcw6RJk2jXrh2urq5Ur16dTp06XfGw1OTJk/H393dOYWFhBbodIoWldUI07Y5uw/bjj2aXIiJSrJg+oDg/oqKiePnll5k5cyabN29m0aJFLFmyhBdeeOGyy4wdO5aEhATndOzYsetYsci1y+qa06sZun4VaNyNiEieuZj1woGBgdhsNk6fPp1r/unTpwkJCbnkMuPGjeO+++7jP//5DwANGjQgJSWFhx56iGeffRar9eKs5u7ujru7e8FvgEghC+7dk4z/c6H82Rgce/dhrVPb7JJERIoF03pu3NzcaNasGZGRkc55DoeDyMhIIiIiLrlMamrqRQHGZrMB6IwSKXFqVg1mc9gNAJxd/IPJ1YiIFB+mHpYaPXo07733HvPnz2fPnj0MHz6clJQUhg4dCsCgQYMYO3ass33v3r2ZNWsWCxYs4PDhwyxfvpxx48bRu3dvZ8gRKSlcbFYONWwNQOYvkVdpLSIiF5h2WApgwIABxMXFMX78eGJiYmjcuDHLli1zDjKOjo7O1VPz3HPPYbFYeO655zhx4gTly5end+/evPTSS2ZtgkihSm/XEb6bQ5kNa8HhgEscehURkdwsRik7npOYmIi/vz8JCQn4+fmZXY7IFX238SitOzXheOWaNF21BAIDzS5JRMQU+fn8NrXnRkSurFHVQFqNmI+riws7A8riZnZBIiLFgPq4RYqwymW98Pd2J9PuYG9MktnliIgUCwo3IkWYxWKhYaUAAPZu2QvZ2eYWJCJSDCjciBRxjSr58/XHY+h/ayvYtMnsckREijyFG5EirmGlAM54B+Q8WLHC1FpERIoDhRuRIq5hJX/WVW4IgF3XuxERuSqFG5EiLtjPg703NM95sGYNZGaaW5CISBGncCNSDHg3acwZL39s6Wnwxx9mlyMiUqQp3IgUAw3DAvj9r0NT/PqrucWIiBRxCjcixUCDSv6sq9wg58GqVeYWIyJSxCnciBQDDUL9WVOlEZ81vpG0+/9jdjkiIkWabr8gUgwE+riTVa0Gz5R9hGoRrWltdkEiIkWYem5EiokGof4A7DieYHIlIiJFm8KNSDHRoJI/LvZsElb+Bl99ZXY5IiJFlg5LiRQTDSv5UyfuCGPeGAn+/tC3L9hsZpclIlLkqOdGpJhoEOrPnqCqJLl5QkIC7NxpdkkiIkWSwo1IMRHg5UZooC+bQuvlzNAp4SIil6RwI1KMNKzkz/qwG3IeKNyIiFySwo1IMdI4LCB3uDEMcwsSESmCFG5EipGGlQLYHlKLDBdXiI2FffvMLklEpMhRuBEpRuqH+pHt6sqWCrVzZujQlIjIRRRuRIoRLzcXagX78mb7e/njo29h0CCzSxIRKXIUbkSKmZxxN/VZWb4WuLubXY6ISJGjcCNSzDQKCwBg2/F4U+sQESmqFG5EiplGlQIA8Fq5AuPhh2HpUnMLEhEpYnT7BZFiplawDx6uVprv24hl/SLIzoZevcwuS0SkyFDPjUgx42Kz0iDUnw2V/rrezerV5hYkIlLEKNyIFEONKgWwsVLdnAd//glxceYWJCJShCjciBRDjSsHEO/px9EKVXNm/PabuQWJiBQh1xRuoqOjWb16NT/99BObN28mIyOjoOsSkStoWrkMAL+F/NV7o0NTIiJOeQ43R44c4amnnqJKlSpUrVqVjh07ctNNN9G8eXP8/f3p3r07X375JQ6HozDrFRGgYoAnIX4e/FHprzuEK9yIiDjlKdw89thjNGrUiMOHD/Piiy+ye/duEhISyMzMJCYmhqVLl9KuXTvGjx9Pw4YN2bBhQ2HXLVLqNa0S8Peg4vh4yMoytR4RkaIiT6eCe3t7c+jQIcqVK3fRc0FBQXTp0oUuXbowYcIEli1bxrFjx2jRokWBFysif2tauQxLd5RnzNQlvDFSp4KLiFyQp3AzefLkPK/wxhtvvOZiRCTvmvw17mZFkiuGYWCxWEyuSESkaNDZUiLFVP1QP9xsVs6mZHLkbCoYhtkliYgUCQUWbp555hnuv//+glqdiFyFu4uNRmH++KUn496nN4SGQmam2WWJiJiuwMLNiRMnOHLkSEGtTkTyoFXVciS6e+O3fQucOgWbNpldkoiI6Qos3MyfP59ff/21oFYnInnQqlpZsFjYFKZbMYiIXKAxNyLFWLMqZXCxWlgVXCdnhsKNiEj+7wo+adKkKz4/fvz4ay5GRPLHy82FBpX82XjhYn5r1oDDAVZ9bxGR0ivf4Wbx4sW5HmdlZXH48GFcXFyoXr26wo3IddaqajneP1ydTHcP3M6fh127oEEDs8sSETFNvsPNli1bLpqXmJjIkCFD6Nu3b4EUJSJ517paWWavdGF7WD2aH9icc2hK4UZESrEC6bv28/Nj4sSJjBs3riBWJyL50Dy8LDarhZ8rNSKjfUcICjK7JBERU+W75+ZyEhISSEhIKKjViUge+bi7UL+iH++2up26AxrRt0kls0sSETFVvsPN22+/neuxYRicOnWKjz/+mJtuuqnAChORvGtVrRzbjifwx6FzCjciUurlO9xMnTo112Or1Ur58uUZPHgwY8eOLbDCRCTvWlUty7urDvHH4XNw9iykpUElhRwRKZ3yHW4OHz5cGHWIyL/QPLwsFgv0WDIf/m8ePPggvPuu2WWJiJhCF8MQKQH8PV2pV8GPA+XCcmboYn4iUooVWLiZOXPmVS/wJyKFp1XVcmwKrZvz4M8/IS7O3IJERExSYOHm66+/Zt68eQW1OhHJp1bVyhLv6cfhkPCcGb/9Zmo9IiJmKbBwExkZyaFDhwpqdSKSTy3DywKwJuSv3hsdmhKRUkpjbkRKiDLebtQJ8WW97hAuIqXcNV/Eb/fu3URHR5OZmZlr/q233vqvixKRa9OmeiA/7vsr3GzZAsnJ4ONjblEiItdZvsPNoUOH6Nu3Lzt27MBisWAYBgAWiwUAu91esBWKSJ51qBXIh2vK83m7/tw19Cbn+1JEpDTJ92Gpxx9/nKpVqxIbG4uXlxe7du1i1apVNG/enKioqEIoUUTyqlXVcri5WBnbdggHb7kTvL3NLklE5LrLd7hZt24dkyZNIjAwEKvVitVqpV27dkyePJnHHnss3wXMmDGD8PBwPDw8aNWqFevXr79i+/j4eEaMGEGFChVwd3enVq1aLF26NN+vK1ISebrZnAOLV+3TqeAiUjrlO9zY7XZ8fX0BCAwM5OTJkwBUqVKFvXv35mtdCxcuZPTo0UyYMIHNmzfTqFEjevbsSWxs7CXbZ2Zm0r17d44cOcJXX33F3r17ee+99wgNDc3vZoiUWO1rBoJhcPSX1fD22/A/4+JEREq6fI+5qV+/Ptu2baNq1aq0atWK1157DTc3N959912qVauWr3W9+eabPPjggwwdOhSA2bNns2TJEj788EOefvrpi9p/+OGHnDt3jrVr1+Lq6gpAeHh4fjdBpETrWLs8k5fuYeRrj0JaIrRoARERZpclInLd5Lvn5rnnnsPhcAAwadIkDh8+TPv27Vm6dOlFdwy/kszMTDZt2kS3bt3+LsZqpVu3bqxbt+6Sy3z33XdEREQwYsQIgoODqV+/Pi+//PIVBzFnZGSQmJiYaxIpyWoH+xLs78GGSvVyZuiUcBEpZfIdbnr27Em/fv0AqFGjBn/++SdnzpwhNjaWLl265Hk9Z86cwW63ExwcnGt+cHAwMTExl1zm0KFDfPXVV9jtdpYuXcq4ceOYMmUKL7744mVfZ/Lkyfj7+zunsLCwPNcoUhxZLBY61CzPH5X+OiVcVyoWkVKmQC7iV7Zs2etyyqnD4SAoKIh3332XZs2aMWDAAJ599llmz5592WXGjh1LQkKCczp27Fih1ylitg61yrMh7B/h5q/eVhGR0iBP4WbYsGEcP348TytcuHAhn3766VXbBQYGYrPZOH36dK75p0+fJiQk5JLLVKhQgVq1amGz2Zzz6tatS0xMzEUXE7zA3d0dPz+/XJNISdeuRiB7QqqT4uoB58/D7t1mlyQict3kKdyUL1+eG264gV69ejFr1iw2bNjAiRMnOHv2LAcOHOC7777jySefpHLlykydOpUGDRpcdZ1ubm40a9aMyMhI5zyHw0FkZCQRlxn82LZtWw4cOOAc8wOwb98+KlSogJubW142RaRUKOPtxg2Vy7G5Yp2cGRp3IyKlSJ7CzQsvvMC+ffto27YtM2fOpHXr1lSuXJmgoCBq167NoEGDOHToEO+++y6///47DRs2zNOLjx49mvfee4/58+ezZ88ehg8fTkpKivPsqUGDBjF27Fhn++HDh3Pu3Dkef/xx9u3bx5IlS3j55ZcZMWLENWy6SMl20aEpEZFSIs+nggcHB/Pss8/y7LPPcv78eaKjo0lLSyMwMJDq1atf05ibAQMGEBcXx/jx44mJiaFx48YsW7bMOcg4Ojoaq/Xv/BUWFsZPP/3EqFGjaNiwIaGhoTz++OM89dRT+X5tkZKuY63yjLyhM7uqNWTOOw9f+43kRESKGYtx4eZQpURiYiL+/v4kJCRo/I2UaHaHQbMXlxOfmsUX/42gZdWyZpckInLN8vP5XSBnS4lI0WOzWuhcOwiAyD2nr9JaRKTkULgRKcG61g2iTuxhar0yDiZPNrscEZHrQuFGpATrUKs8VZJiuX3112TOnWd2OSIi14XCjUgJ5ufhiqVtOwDc9u+DON0pXERKvnyHmy5duhAfH3/R/MTExHzdfkFEro9WLWqyN7ByzgOdEi4ipUC+w01UVNQlrwacnp7Oal0oTKTI6VY3mA1/3Wcq/dcoc4sREbkO8nzpi+3btzt/3r17d66bW9rtdpYtW0ZoaGjBVici/1pYWS+O39AMtv5I2q8r8TC7IBGRQpbncNO4cWMsFgsWi+WSh588PT2ZPn16gRYnIgUj4MYu8OmL+P25E5KTwcfH7JJERApNnsPN4cOHMQyDatWqsX79esqXL+98zs3NjaCgoFw3tBSRoqNlh8Yc8w/GzWGn7P6DuDZpZHZJIiKFJs/hpkqVKgC5blopIsVD40oBdHt4Bofs7nziE0o7swsSESlE13S7mf3797NixQpiY2MvCjvjx48vkMJEpOBYrRZaNK3JoY3H+GXPadrVDDS7JBGRQpPvcPPee+8xfPhwAgMDCQkJyXXDTIvFonAjUkR1qxfMwo3HiNx5kgm9amNx0a00RaRkyveNM6tUqcLDDz9cbO/ErRtnSmmVlmlnUds+3Lwzivi5HxN+Vx+zSxIRybNCvXHm+fPnueOOO665OBExh6ebjXB3g4D0ZM58s9TsckRECk2+w80dd9zBzz//XBi1iEghc+neDQD/tatMrkREpPDk6bDU22+/7fw5JSWFN998k5tvvpkGDRrg6uqaq+1jjz1W8FUWIB2WktLs3P7DlK1VDQcWTh2MJrRaJbNLEhHJk/x8fucp3FStWjVPL2yxWDh06FDeqjSJwo2UdscrVKVSzBEiX5hB1+ceNrscEZE8yc/nd55Olzh8+HCBFCYi5kto25FKXx8he/kvoHAjIiVQvsfciEjxFtKvFwC1dvxBXFKGydWIiBS8fF/oYvTo0Zecb7FY8PDwoEaNGtx2222ULVv2XxcnIgWv3M092V21Pr+E3EDQjhPc1aaa2SWJiBSofF/npnPnzmzevBm73U7t2rUB2LdvHzabjTp16rB3714sFgu//fYb9erVK5Si/w2NuRGBGSsO8PpPe+lUuzzzhrY0uxwRkasq1Ovc3HbbbXTr1o2TJ0+yadMmNm3axPHjx+nevTt33303J06coEOHDowaNeqaN0BEClfPG0IAWHPgDInpWSZXIyJSsPLdcxMaGsry5csv6pXZtWsXPXr04MSJE2zevJkePXpw5syZAi22IKjnRiTHrS/+QODWDdz6zIP0aapTwkWkaCvUnpuEhARiY2Mvmh8XF0diYiIAAQEBZGZm5nfVInK9ZGXx5eS7+PDrSWxZutrsakRECtQ1HZa6//77Wbx4McePH+f48eMsXryYBx54gD59+gCwfv16atWqVdC1ikhBcXUlK6ItAN6/Lic+VV9GRKTkyHe4mTNnDl27duWuu+6iSpUqVKlShbvuuouuXbsye/ZsAOrUqcP7779f4MWKSMHx6dMbgPYHN7JsZ4zJ1YiIFJx8j7m5IDk52Xk14mrVquHj41OghRUWjbkR+cvBg1CjBllWG/998yc+fLyr2RWJiFxWoY65ucDHx4eGDRvSsGHDYhNsROQfqlcnq3pNXB123Fb+yunEdLMrEhEpEHm6iF+/fv2YN28efn5+9OvX74ptFy1aVCCFiUjhc72lF7z1Fh0PbuSH7ad4oF3e7iMnIlKU5annxt/fH4vF4vz5SpOIFCM33QRAp0Ob+G7rCZOLEREpGNc85qa40pgbkX9ITydp5rvcvN+baP8QosZ0IjzQ2+yqREQuUuhjbrKzs/nll1+YM2cOSUlJAJw8eZLk5ORrWZ2ImMXDA9/Rj1GleX0Avt920uSCRET+vXyHm6NHj9KgQQNuu+02RowYQVxcHACvvvoqY8aMKfACRaTw3dqoIgDfbjtJKevMFZESKN/h5vHHH6d58+acP38eT09P5/y+ffsSGRlZoMWJyHXgcHDL+iXM/XoicUdPsedUktkViYj8K3k6W+qfVq9ezdq1a3Fzc8s1Pzw8nBMnNCBRpNixWvGcNYPOB3bS+dBGvtvWmHoVNR5NRIqvfPfcOBwO7Hb7RfOPHz+Or69vgRQlItdZ374A9Ny3ju91aEpEirl8h5sePXowbdo052OLxUJycjITJkygV69eBVmbiFwvf90XruPhTZyJi2dz9Hlz6xER+RfyHW6mTJnCmjVrqFevHunp6dxzzz3OQ1KvvvpqYdQoIoWtSRMIC8MrK4N2R7fy7VadNSUixVe+w02lSpXYtm0bzzzzDKNGjaJJkya88sorbNmyhaCgoMKoUUQKm8Xi7L3pse93lu44RbbdYW5NIiLXKM8Dijt27EjXrl3p1KkTERER3HvvvYVZl4hcb336wPTp9Di4nrGJaaw9eJYOtcqbXZWISL7lueematWqzJ07l06dOhEQEEC3bt146aWX+P333y85wFhEipkOHSAoiPPValEuNVGHpkSk2MpzuJk3bx6HDx/m0KFDTJ8+ndDQUN59913atGlDmTJluOmmm3j99dcLs1YRKUwuLnD0KGe/X0acTxl+3hVDepa+uIhI8fOv7y116NAhPvzwQ6ZPn05ycnKR78XRvaVErszhMGj36q+cTEhn9r1NubF+BbNLEhHJ1+d3vi/iBzm3YIiKinJOsbGxtG7dmo4dO15TwSJSdFitFu4Kc2XBsWi+23ZS4UZEip08h5uPPvrIGWbOnDlDmzZt6NixIw8++CAtWrTA1dW1MOsUketlyhQeffJJAht05/myISSmZ+Hnofe3iBQfeQ43Q4YMoXLlyjz99NM88MADCjMiJVWTJlgcDm7Zv44JGcNYvPkEg9uEm12ViEie5XlA8cyZM2ndujUTJ04kKCiI3r17M2XKFDZu3KhLtYuUJB07QnAwfqmJtD2ylU//OKr3uIgUK3kON8OGDWPBggWcOnWKNWvW0KtXL9avX8/NN99MmTJluPnmm3njjTcKs1YRuR5sNrjzTgD6/bmafaeT2XBEt2MQkeLjX58tdfLkSWbOnKmzpURKknXroE0bMjw8aTpsPt1a1eCtu5qYXZWIlGKFerZUbGwsK1ascA4u3rdvH66urrRu3ZrOnTtfc9EiUoS0bg21auG+bx+99v7Gt14+jL8lg3I+7mZXJiJyVXk+LPXwww9Tr149KlSowKBBg9i5cyf9+/dn+fLlxMfHExUVxYQJEwqzVhG5XiwWuP9+AO49+BuZdgdfbjpuclEiInmT556bLVu20KdPHzp37kzbtm3x8vIqzLpExGyDBkGZMhyq1RaWHeGzP6J5qH01rFaL2ZWJiFzRvx5zU9xozI1I/qRl2mn18i8kpmczb2gLOtUOMrskESmF8vP5nefDUiJSOnm62bi9aSgWw8Gnf0SbXY6IyFUp3IjIlS1YwNNP30XPfeuI3HOak/FpZlckInJFRSLczJgxg/DwcDw8PGjVqhXr16/P03ILFizAYrHQp0+fwi1QpDTbuRP3Pbt4bNePOAxYsOGY2RWJiFyR6eFm4cKFjB49mgkTJrB582YaNWpEz549iY2NveJyR44cYcyYMbRv3/46VSpSSg0fDjYb9fZvpW7sIRasjybL7jC7KhGRyzI93Lz55ps8+OCDDB06lHr16jF79my8vLz48MMPL7uM3W5n4MCBTJw4kWrVql1x/RkZGSQmJuaaRCQfQkOhf38A/rt9KbFJGfyy+7TJRYmIXF6BhZvBgwfTpUuXfC2TmZnJpk2b6Nat298FWa1069aNdevWXXa5SZMmERQUxAMPPHDV15g8eTL+/v7OKSwsLF81igjw6KMA3LxjBWVSEzSwWESKtAILN6GhoVSpUiVfy5w5cwa73U5wcHCu+cHBwcTExFxymd9++40PPviA9957L0+vMXbsWBISEpzTsWMaLyCSb23aQNOmuGZmMGjLEn47cIZDcclmVyUickkFFm5efvll5s6dW1Cru6SkpCTuu+8+3nvvPQIDA/O0jLu7O35+frkmEckniwWefBKAh7b8gFdmGnPXHDG3JhGRy8j3vaUKUmBgIDabjdOncx+/P336NCEhIRe1P3jwIEeOHKF3797OeQ5HzsBGFxcX9u7dS/Xq1Qu3aJHSqn9/GDSIo11vI3WXG19sPMZjXWtS3lf3mxKRoiXf4cYwDL766itWrFhBbGysM1xcsGjRojyvy83NjWbNmhEZGek8ndvhcBAZGckjjzxyUfs6deqwY8eOXPOee+45kpKSeOuttzSeRqQw2Wwwfz51DYPGM9ey9Vg8H645zFM31jG7MhGRXPIdbkaOHMmcOXPo3LkzwcHBWCz/7j4zo0ePZvDgwTRv3pyWLVsybdo0UlJSGDp0KACDBg0iNDSUyZMn4+HhQf369XMtHxAQAHDRfBEpHBaLheGdqjNi3h98vO4o/+1QjQAvN7PLEhFxyne4+fjjj1m0aBG9evUqkAIGDBhAXFwc48ePJyYmhsaNG7Ns2TLnIOPo6GisVtPPWBeRC9LT6fHFLH778F26D3qL91cfZkzP2mZXJSLilO8bZ1atWpUff/yROnWKZ1e0bpwp8i9lZ0P9+rB3L+9E3Mns7vez+snOlPFW742IFJ5CvXHm888/z8SJE0lL0/1lREolFxd49VUAHtz4DT5xMby3+pDJRYmI/C3f4ebOO+/k/PnzBAUF0aBBA5o2bZprEpFS4NZboV073LMyeSbqQ+auOUJMQrrZVYmIANcw5mbw4MFs2rSJe++9t0AGFItIMWSxwLRpGC1bcuueVXy5rxuvLQvhzQGNza5MRCT/4WbJkiX89NNPtGvXrjDqEZHiolkzLI8+Cm+9xYs/z6RHpRu4N6IKTSuXMbsyESnl8n1YKiwsTANxRSTHCy9AaCgVUuO5IfYQE7/bhcORr3MUREQKXL7DzZQpU3jyySc5cuRIIZQjIsWKry988glJa39nX7UGbDuewNebj5tdlYiUcvk+FbxMmTKkpqaSnZ2Nl5cXrq6uuZ4/d+5cgRZY0HQquEjhmLPyIJN//JOy3m78NLKDbssgIgUqP5/f+R5zM23atGutS0RKsKFtq7L/6x/p8N2XPBv6CnOGttYJByJiinyFm6ysLFauXMm4ceOoWrVqYdUkIsWQW1oKr348Dlt8PGdnTubLhjO4s7nu9yYi11++xty4urry9ddfF1YtIlKc+fpimzsXgKGbvufQsy9x/HyqyUWJSGmU7wHFffr04ZtvvimEUkSk2OvTB8fLLwPw9M9z+H7E82RmO0wuSkRKm3yPualZsyaTJk1izZo1NGvWDG9v71zPP/bYYwVWnIgUP9annybxVBx+06fy38/f4FtfN/rMfkHjb0TkurmmG2dedmUWC4cOFe17zOhsKZHrwDA4Nvghwj5+H4CFz03nzkkjFHBE5JoV6tlShw8fvubCRKSUsFgIm/8uOzz9iPl1DWMzKnPwxz8Ze1MdBRwRKXT5Djf/dKHTR3+sROQiFgsN5kxh+5pDOL7fw7urDpEen8CEwCRsPbqbXZ2IlGD5HlAM8NFHH9GgQQM8PT3x9PSkYcOGfPzxxwVdm4iUAAPbVuOVfg2wWKDeqxOw9exB+q19oIgfwhaR4ivf4ebNN99k+PDh9OrViy+++IIvvviCG2+8kWHDhjF16tTCqFFEirm7WlZmxoBGpHj6kG2x4vH9tzhq1cIYOhT27jW7PBEpYa5pQPHEiRMZNGhQrvnz58/n+eefL/JjcjSgWMQ8u08mMnP6Yu5c8BYdjmxxzjfat8fy1FNw880mViciRVl+Pr/z3XNz6tQp2rRpc9H8Nm3acOrUqfyuTkRKkXoV/Zj24n3snvcldwx5k+U1WmG3WLGsXs2mb1eQkJqV0/DMGYiMhMREcwsWkWIp3wOKa9SowRdffMEzzzyTa/7ChQupWbNmgRUmIiWTi83KsI7V6dtkOB+v68XrP2+k+4afWOaowaEXfqZ2sC8PRK/jjqlPY1gsZFWrjqV2bVzq1MZSuzaEhkLr1lCunNmbIiJFVL4PS3399dcMGDCAbt260bZtWwDWrFlDZGQkX3zxBX379i2UQguKDkuJFC1pmXYWbTnOR2uPsvd0EgD9d/zC42s+Jyzh9CWXeePJd4ht0Y7QAC/arv6Ouh9OJ9vXF/z8MQICsPj7Yw3wx1YmANvgwbjXrZ2zYEwMREeDvz8EBOT86+FxnbZURP6NQr3Oze23384ff/zB1KlTnbdhqFu3LuvXr6dJkybXVLCIlF6ebjYGtqrCwFZViE1KZ8Ph82yOrsrzd9xD2rETeB/YS3BMNNXOHqdK/CmCk8+x5LwrhzceByBj3XaaHz962fXfFe3L1uqNKePlxsBNP/DI19NyPW93dSPbxxdHuXIkvT4V/1t64u5iyzmba9s2qFkTatRQCBIpRvLdc1PcqedGpPhJy7RzOjGdmMR04pIySErP5kxyBsfPp5Jy9ARlTx/HLSURa0IibilJuKUk4ZGSjGdaMu+26sdx/2AABmz7iUfXLsAvIxWfjFSs5P7zd+c9r7A+rD5Bvu4M2/Mz93/6OgCG1YqjShVsdepAw4bQsiV07gxlylz334VIaZWfz2+FGxEpsewOg+SMbBLTsjiXkkl8WhYJaVmcTkgnJj6V9PMJOM7Hk3n2HFmxsaz2Dee8S04Pza27V3L/xm+odu4kfhkpF617+duf4NPrRmoG+1AuPg4LQKVK13cDRUqRQjksZbVar3olYovFQnZ2dl5XKSJSqGxWC/6ervh7uhJW1uuq7Q3D4HxqFkfPprA/tiFLY4eyPyaRs4eO4XnoANXPHad+zEEandrH6EOuJL33OwATf5vP4DVfklCtFtxxB35D7sVSp05hb56IXEaee26+/fbbyz63bt063n77bRwOB+np6QVWXGFQz42IXIuUjGwOxiWz/3Qy+2OTORCbxP7YZKLPpfLGD2/SZ1cUNsPhbH+yah2yh95P2GMPYfH3N7FykZLhuh2W2rt3L08//TTff/89AwcOZNKkSVSpUuVaV3ddKNyISEFKz7Kz40QCf2zcj+P772n420+0PbIFV4cdgFj/8nz99WpubVaZ0ABPk6sVKb4KPdycPHmSCRMmMH/+fHr27MnkyZOpX7/+NRd8PSnciEhhSsu0s3XrAU7N+pDG33/GyvAmTOz2XwAiqpVjVF1vWrSrrxsOi+RToV2hOCEhgaeeeooaNWqwa9cuIiMj+f7774tNsBERKWyebjYiWtam39xXKX/sIL5TXyeiWjksFsheuZKmHRvzU8fb+X39PkrZ+Rwi102ew81rr71GtWrV+OGHH/j8889Zu3Yt7du3L8zaRESKNV9PN/q3r83nD7Xmt6e68Ij9CC6GgxtXL6ZOx2Z8cM//se9kvNllipQ4eT4sZbVa8fT0pFu3bthstsu2W7RoUYEVVxh0WEpEzHRuyc9kjXiE4KP7AdgSWocdL73FXQO74eaS79v9iZQahXIq+KBBg3SMWETkXyp7cw/ouZv4qdNxmzCOJif+pO5/ejNryaN0eXMcDSrpzCqRf0sX8RMRMYkRHU3cnfcS9MdqxvQayeJG3RndvRbDO1bHatWXSZF/KrQBxSIiUnAslSsTtG4liQu/Jv2ee7E7DF7/aS8PfrCO+NRMs8sTKbYUbkREzGSx4HdnP94Z2IzX+jckODOZsU8PYNp/X+ZgXLLZ1YkUSwo3IiJFxJ3Nw/jetoMa547z/CcTWTxwFGv2x5ldlkixo3AjIlKEBL0ykbRhDwMwZvn7HL77fiJ3njS5KpHiReFGRKQosdnwnPkO2a+/gcNi4d5NP5B09338tifG7MpEig2FGxGRosZiwWXMExiffordaqPPzl9J7D+AzQdjza5MpFhQuBERKaJsd9+NY+FCsm0u1D+xlydn/8ohDTIWuSqFGxGRIsy1/+1kf/MtL4yazgGbL4Pnric2Kd3sskSKNIUbEZEizuOWXkx+4laqlPPi2Lk0Jkz+kuT0LLPLEimyFG5ERIqBQB935g9tyd2H1vD2q0P56a4RZNkdZpclUiQp3IiIFBPhgd6MqB+Aq8PO7d++x9ePvUgpu4OOSJ4o3IiIFCOVxj/Jkf+OBKDf7BdY+s5CcwsSKYIUbkREipnwmVM43LkXbo5s2jz5ENtWbDS7JJEiReFGRKS4sVoJ/+FLjlSvT5n0JALu7Evc0VNmVyVSZCjciIgUQxYvL4J+/ZHYgCCqnDnOz49P1ABjkb8o3IiIFFNelSuR9fVipnUZynO1b2by0j/NLkmkSFC4EREpxkK7tKHe2y9jWKx8uOYw323TTTZFFG5ERIq5HjeE8HCn6nhmppPywEMc2bTb7JJETKVwIyJSAjzRozZzfv+AuzctIf32O8hI0y0apPRSuBERKQFsVgv13p1KoocPdY7uZuOQx8wuScQ0CjciIiVEYP3aHJ48DYCIL95n5yeLzS1IxCQKNyIiJUijkQ+wvtvtWDEIfvhBEo/p+jdS+hSJcDNjxgzCw8Px8PCgVatWrF+//rJt33vvPdq3b0+ZMmUoU6YM3bp1u2J7EZHS5oaFH3A0qDLlk85ytM9doPtPSSljerhZuHAho0ePZsKECWzevJlGjRrRs2dPYmNjL9k+KiqKu+++mxUrVrBu3TrCwsLo0aMHJ06cuM6Vi4gUTd5l/Umb/wkZNhcq7tnKzz/+YXZJIteVxTD5lrKtWrWiRYsWvPPOOwA4HA7CwsJ49NFHefrpp6+6vN1up0yZMrzzzjsMGjToouczMjLIyMhwPk5MTCQsLIyEhAT8/PwKbkNERIqY7ybN5IVYX9IDg1g2sgOhAZ5mlyRyzRITE/H398/T57epPTeZmZls2rSJbt26OedZrVa6devGunXr8rSO1NRUsrKyKFu27CWfnzx5Mv7+/s4pLCysQGoXESnqej07jNC61UhKz+aJL7bicOjwlJQOpoabM2fOYLfbCQ4OzjU/ODiYmJiYPK3jqaeeomLFirkC0j+NHTuWhIQE53Ts2LF/XbeISHHgYrMydUBjvNxsBPz4PWvGvGB2SSLXhYvZBfwbr7zyCgsWLCAqKgoPD49LtnF3d8fd3f06VyYiUjRUDfRmeoUEur4wmUybCwdv6UH1Lq3NLkukUJnacxMYGIjNZuP06dO55p8+fZqQkJArLvvGG2/wyiuv8PPPP9OwYcPCLFNEpFjr8t872NaoHW72bByDBpGekmZ2SSKFytRw4+bmRrNmzYiMjHTOczgcREZGEhERcdnlXnvtNV544QWWLVtG8+bNr0epIiLFlsVqpdJXHxPv6UvNE/v54z9PmF2SSKEy/VTw0aNH89577zF//nz27NnD8OHDSUlJYejQoQAMGjSIsWPHOtu/+uqrjBs3jg8//JDw8HBiYmKIiYkhOTnZrE0QESnyytUI5/hLbwDQduEctnwTeZUlRIov08PNgAEDeOONNxg/fjyNGzdm69atLFu2zDnIODo6mlOn/r7C5qxZs8jMzKR///5UqFDBOb3xxhtmbYKISLFQf9RD7GjTAxfDgffDw0hJ1uEpKZlMv87N9Zaf8+RFREqalOgTZNepi39aEgufm86AFx4xuySRPMnP53exPltKRETyx7tyKDtfe5vHVx5lVXZV6hyLp1FYgNlliRQo0w9LiYjI9VX/kSH4334rDgOe+no7WXaH2SWJFCiFGxGRUmjcLfUI8HIlft9hfn7tfbPLESlQOiwlIlIKBfq481o9V1q/9DAuhp3jPdpSqdkNZpclUiDUcyMiUkp179uBE+G18crK4Px9QzEcOjwlJYPCjYhIKWWx2fD7eC7pLm402LOBLRPfNLskkQKhcCMiUoqFtmzEpgdGAVDzteeJP3DE3IJECoDCjYhIKddi2iT+DKuDb3oKx+4eCqXr8mdSAinciIiUcm4ebjg+eJ9MqwsNNkaxbfo8s0sS+Vd0tpSIiFCve1tWDnyYfX9G89HZ8vyQmoW/l6vZZYlcE/XciIgIAK0+eJPP73ycYxkWJv2w2+xyRK6Zwo2IiADg4Wrj9TsaYrHA4o1HWffT72aXJHJNFG5ERMSpWZWyPF7Hmy8+e5qa/W8m8cRps0sSyTeFGxERyeW/fVtQPiuFwORz7Bv4oNnliOSbwo2IiOTi6edNyqx3cWCh+crv2f7e52aXJJIvCjciInKRev168sdtgwAI/r/HSTx9xuSKRPJO4UZERC6p0by3OV6uIsEJcey5b5jZ5YjkmcKNiIhckleAH4nTZwHQavnXbJ//tckVieSNwo2IiFxWvbtvZf1Nd7EzuDpvbD5HUnqW2SWJXJXCjYiIXFH9z+bwyOMzWeVZkck//ml2OSJXpXAjIiJX5BXgx8t3NgXgsz+iWbcj2uSKRK5M4UZERK6qTfVABjevyKjVn1C1fQuST+riflJ0KdyIiEiePNm1OrftX0tIQixH+t8HhmF2SSKXpHAjIiJ54l3Gj4R3PyTLaqP+uuXsnzLL7JJELknhRkRE8qxRn66suGs4ABWee5LUfQdNrkjkYgo3IiKSLxFzXmN7lRvwyUghtt9dYLebXZJILgo3IiKSL74+nmR8MJdkN0/Cd23k6BPPmV2SSC4KNyIikm8turZg6X+fI93FjXnRWcQmpZtdkoiTi9kFiIhI8dT7jad4oEwt1mR48ufnW/nkP62wWS1mlyWinhsREbk2nm42Jj1yI95uNtYdOsusxRs0/kaKBIUbERG5ZtXL+/ByvwY0PbGHvvffQvTIsWaXJKJwIyIi/85tjUMZUC6b0MQ4Ks14g/gvFpldkpRyCjciIvKv3fb6//FDm1uxGgZug+/Dvmu32SVJKaZwIyIi/5qHq426X8xlY+X6eKWnktjjJjh/3uyypJRSuBERkQJRPbQscXM/4bhfecqcjCb25r6QnW12WVIKKdyIiEiBualLI356cTapru4ErVvJ8edeNLskKYUUbkREpEANGdGXeQ9NJKpqM+6wNmb/6SSzS5JSRuFGREQKlM1qYcibTzBt9DROOVwZ+P4fHD2bYnZZUooo3IiISIHzcnNh7tCW1A72JTYxnaUDR3JOp4jLdaJwIyIihaKMtxsf/6clD0avZfhP7+Mz8C4SvvnB7LKkFFC4ERGRQhPk68Hgt55iZb22uGVn4XHn7SR/v9TssqSEU7gREZFCVSnIn8rLv2dlnQjcszJxvb0v579dYnZZUoIp3IiISKGrWrEMIT99z6o6rXHPysT79r6c+fBjs8uSEkrhRkRErovalctR+ZclrGjQATd7FgEPDuHw+u1mlyUlkIvZBYiISOkRHloWj6ilfNf7XjZ4h/D9shhmlq1ImxqBZpcmJYh6bkRE5LoKKetN+58XsrPPvcSnZnHvB38w9+NI7GfOml2alBAKNyIict2V8Xbj8wdbc3vTSnilp9J+5BDi6jYkbvXvZpcmJYDCjYiImMLD1caUOxsxpWMI7o5sQs6cxL9zBzY/9iyObLvZ5UkxpnAjIiKm6tmnPY4//uCPBu1ws2fRdPrL7L6hJYc37DS7NCmmFG5ERMR0VWpVpvmWlax58mVSXd2pv28zIW2aEzloJAnJGWaXJ8WMwo2IiBQJNpuVtq+OJXHNH+yt0xTP7Awc69bR/o0oJn2/mwOxyWaXKMWExTAMw+wirqfExET8/f1JSEjAz8/P7HKkFEnLSmPtsbVsOLmBjSc3EpsSS5WAKnSv1p3u1bpTwbeC2SWKFB2Gwd43ZvJifDlW230BqJAYRy/3JBoO6suNDSrg7mIzuUi5nvLz+a1wI1LIjiceZ8b6Gby7+V3OpZ27ZBsLFh5s+iCTu02mrGfZ61yhSNFldxis2h/Hp79H0+PNZ7hz+3K2VqjJwrb9CbjvLga0qU54oLfZZcp1oHBzBQo3cj0YhsG64+t464+3+Hr319iNnDM/KvpWpF3ldrSo2ILK/pXZfno7Px38iY0nNwJQ3qs88/rMo1fNXte95jOpZ1h2YBk7Y3ey/fR2dsbu5HTKaawWK2F+YVQtU5Vw/3DqBNahfZX2NAlpgs2qb85ynRgGyY88jsd7c3DJygTglE85vm7QlSM39aVJjwjaVg+kSjkvLBaLycVKYVC4uQKFGylMmfZMvtj1BW/98ZYzsAB0rNKRx1s9zq21b71kIFh9dDXDlwxnV9wurBYrs26exUPNHrouNW8/vZ2pv0/l8x2fk2HP+8DNAI8AOlbpSN86feldu7d6nOT6iIvDPnMW2e/MwP1MrHP2D7Xb8UifpwnydadJ5QBahJeldbVy1K3gh82qsFMSFLtwM2PGDF5//XViYmJo1KgR06dPp2XLlpdt/+WXXzJu3DiOHDlCzZo1efXVV+nVK2/fdBVupDCcTj7NnE1zmLVxFjHJMQC429y5p8E9PNbqMRqHNL7qOjLtmQz/YTgfbv0QgJm9ZjK8xfBCq/nPM3/ywqoX+HzH5xjk/BloGNyQtmFtaRjckAZBDQjzDyPbkU10QjSHzx/mcPxhtsRsYdXRVSRmJDrX5WJ1oXu17txa+1b61OlDiE9IodUtAkBGBnz7LWkfzsPtl+UsunEQzzS6nSy7gVdmGs/9+j6bQ+vyZ7X6lG9SnyZVytIoLIBqgd6U93XHw1W9jsVNsQo3CxcuZNCgQcyePZtWrVoxbdo0vvzyS/bu3UtQUNBF7deuXUuHDh2YPHkyt9xyC5999hmvvvoqmzdvpn79+ld9PTPDjWEYOAzHNU12w06WPYtMeyaZ9kwy7BmkZ6eTlpVGWnYa6dnp2B12gn2CKe9VHgODbEc22Y5ssuxZzp8TMxJJzkzGYTgwyKnnQl35fewwHNisNnzdfPH38KeMRxk8XDxwtbniZnPDzeaGu80dL1cvAjwCcLW54mp1xdXmis1iK/Zdx+nZ6UQeimThroUs3LWQTHtOV3kFnwo83OJh/tvsv5T3Lp+vdRqGwZifx/Dm728CsOD2BQyoP6BA6z6Xdo4xP49h3tZ5zlDTv15/RrceTetKrfO0X7Id2Ww5tYUl+5ew+M/FbD/9980PLVhoW7ktvWv1pnu17jQOaVzs93VBMQyDDHsGKZkpxKfHE58eT1JmEhnZGaRkpZCSmUJyZjIpWSnO96wFC0HeQfi4+eBmc8PX3ZcKPhXwdffF180XT1dP5/vR29W7xB0qdBgOUjJTSMxIJDEjkYSMhL9/Ts/52Z6UABYLpxzJ7I49SvzhvVTYvR2bA2wGOHAlycOXc55exHt4cMbblywPV9xcHLjaHLi72vB288DHzR0/D0/83D3xdvfAw8Xd+bfMzeaGq9UVm9WGBQs2q41sRzbert6U9SyLp6tnzt+/v/7GuVhdcLW64uvui4eLBzaLDavFis1qw2ax4WJ1IcQnpMTtr8JUrMJNq1ataNGiBe+88w4ADoeDsLAwHn30UZ5++umL2g8YMICUlBR++OEH57zWrVvTuHFjZs+efVH7jIwMMjL+7mpPTEwkLCyMvcf3kmJN4VjiMc6mniXLkZUrPGTaMzmffp641DjSstLIsGfkhIrsDJIzk7Ebdsp6liUjO4MzqWfIsGdgd9ixG/bL/puRneH8MBFwtbpSzqscwd7BBHgE4O7i7gxDgV6BVPKrRIhPCEHeQfi6+TpDUqBXID5uPqb8UUjKSGLp/qUs/nMxS/YvITnz71NTW4a25PFWj9O/Xn/cbG7X/BqGYTD6p9FM+2Ma7jZ3Fg1YVCBjcAzDYMHOBYz8aSSxKTnd+bfVvo0JHSfQpEKTf7XuPXF7+Hbvtyz+czHrT6zP9Zy3qzd1Aus4p7qBdakTWIdqZarh4eJx2eCTac8kKSOJ5MxkkjJz/r0QCs6nn8fusONidcFmzfnQsDvsZDuyne+5bEc2BgY2i80ZuF2sLrhYXfBx88HHzQcLFrIcWQDOMG6xWLBarM7H7i7uuFhdSM9OJzUrNdeUkply0bzUrFQy7BnOLxbn089zLPEYxxKOcTrlNA7D8a9+11fianXF280bf3d/3F3cMQwDA8P5b0Z2BmnZaaRlpeEwHHi4eDjfdxd+9nDxcP4eLnxAu7u44+3qnetD3sXqgoeLBwYGaVlpGBgEegXmvJdt7jgMB1aL1fl6Fz7Ys+xZZDmynIEsLSuN5MxkEjNzB5Z/TiX176bNYqOcVzkchoPUrFTn78nP3Y8AjwBnoHX/K2S529zJdmQ7A7GL1QVPV088XTxxsbpgYODv7k9Zz7L4uvnm+jKcac8k25HtXLeni6dzn1gt1jxNNovNGc7+GdL+Oc/V6orVYs3Znr/2fZYjCw8XD7xcvfB08bzmv90pSSn0a9Kv6IebzMxMvLy8+Oqrr+jTp49z/uDBg4mPj+fbb7+9aJnKlSszevRoRo4c6Zw3YcIEvvnmG7Zt23ZR++eff56JEyde/OJPAx4FsBGFxILF+R/qwh9bd9vf3yJcba54ung6vy14unhisVg4nXyaM6lnsFltzj/k/5x83Xzxdff9e91/vc6F18jz47/+vfBGi0+PJyEjIVc4zMjOcH4gpGSlFPjvyNvVm4q+FankV4lQv1CCvYMJ9ArMNZX3Kk85r3L4uvni7uJ+1XUahkFKVgonEk9wOP4wR+KPOA/HHI4/zPbT2509NAChvqH0qdOHexveS+tKrQts2+wOO3d+dSeL9izCy9WL1UNX07RC02teX1JGEg/98BALdi4AoG5gXd6/9X3ahLUpqJKdjiUcY/Gfi4k8HEnkocgr7ntXqys+bj7O3rwLH37x6fH5Gv9THHm7ehPgEYCfux9uNjd83HzwdvPO+dfV29lLYHfYiU2NJS0rjUx7JufSzhGXGkdSRhJp2Wlmb8Z142J1wd/dHz93P/zc/fD3yPnZ180XV5srhmFQ1rMsFX0rEuARkNPj7bBjz8og89gR7KdO4BYTh8upWIw+/UkuW5HkdPCOWkn1JV+QaeOS06J67ThYLgTDkkXVc0doHb2dbKuNDJsNLFZSXCHewyDdBQ6WCyLBwx1w4JmVgn/6eVJcHWRZDRwWA4cFHBYDw5KzbsNSeEG3REoHXiFP4cbl+lR0aWfOnMFutxMcHJxrfnBwMH/++ecll4mJiblk+5iYmEu2Hzt2LKNHj3Y+vtBzY7FYKO9dnjC/MIK8g3KFhgvfTMp4lCHQKxBvN+9cwcLHzQerxcr59PO429wJ9ArM6Xb8n0T7v/9eSN9XS8cWLCWuG9/usDt7xy78m2HP4GzqWU6nnCYpI4n07HQy7Dk9Y3EpcRxLPEZsSiyxKbGkZOV8Qz6fdt75YZmSlcL+c/vZf25/nmq48EHq6+6Lt6s32Y5s5+G9C9M/g8vl1Cxbk351+9G3Tl9ahLbAain4a2HarDYW9l/ILZ/dwk8Hf6L3571Ze/9aqgRUyfe6dsftpu/Cvuw7uw+rxcpz7Z/j2Q7P/qvepSsJ8w/jsVaP8Virx8iyZ3Hw/EH+PPOnc9pzZg9/nvmTxIxEshw5PRtX4m5zdx6CudB7V8azDC5Wl5yeGocdh+G46BvkhcMAdsPuDNwX2idnJpOcmYyBgavVFcD5DRdy/r9m2DPIyM5w9sJ4unji5ep12cnbzRsvFy88XT1xt7k7/5b4uvkS5h9GJb9KVPCpgI+bD16uXrjaXP/17zrbkU16drrz78aZ1DOkZKWQkJ7zRcNisTj/nliw4O7ijqdLzhciq8Xq3MYL770L25uenY5h5BzaznJkOb+k/PPLi91hJy07p6fB08UTgNiU2JxeBXsGFiw4DIezZ+FCT82FwzYXar7wu/N393cGFmeA+UeYuVIv379SrRvU7AHHj8PZs3DmDEbcGbLj4jDiznDzA6M52bg1KRnZVFw4n3Y/bAOy/5pymzDofpbXaUdGtoPOOyJ5Y/Grl33ZUb0e58sGTXFYEmkdvZMpS+fgkQ0OCyR4QII7pLrmhKAPm99IVLVGGGRS88wJxq78AvdsyLZCugukueb8DPDlDa1YE14NB2mEJCdy97ZfcbeDmx1sjpx1x3tAhg22VqjBzpBqOT0+6Yl0O/DHXwHs4ulg2YrsKx+KgQOP7HSandiF3ZLzunYr2C05/2ZZ4ZyXL7E+IVgMd9ztFjoc3kGGS872pLqSqx8u3sOHE/45Q1CshoO6sUcu+zuLt3hxhNQ87VZTw8314O7ujrv7xd/Yzz55ljIBZUyoqHSyWXM+cDxccneXVfavnO91pWenk5KZwvn085xIPMHxxOMcTzxOXGocZ9POEpcSx5nUM5xJPUNcapxz4OuFD9KrfZhCzrfqqmWqUjUgZwoPCKdqmarUDaxLrXK1rkv4dLG6sLD/Qtp82Ibdcbvp8lEXogZHEeYflud1fLztY4YtGUZqViqV/Cqx4PYFtK3cthCrzs3V5uo8HPVPhmGQmJFIUmYSSRlJuQ4nuVhdnD0aF3p15NIuHGK7ID//N+QvtWrlTP9gAS78r2vy1wRA/Sdg1FBIToaUFMjMhOxsyMqC7GwmNmzIxAtfvo/UxOhTj6z0TLKy7dgdDhx2B3a7gcPhYHjLNvQKDONcSgY+J5qQVDmMdAsYDgfpWXbItOOabcdqd9CzYStuqFyLjGwH5c7G4G2vjIvFgsMAN4eBp93BhWMwHWo0wr9CLVIy7QQln6XWiXrYrBZsFgtWi4X0bDsZWQ7sDgM3jxvwCrgBwwBf63laZS5zrifnUCY4DACDMt51cFjr4zDAKyuJ7sd/xDDAcYmDP7uCq7MmvDFWC5Qxsugc9wPZhkG2/eK2+wKrEOXRHABXexZ37P3usrvqT9/yHOG1POxUk8NNYGAgNpuN06dP55p/+vRpQkIufbZFSEhIvtpfjgZxFV8eLh54uHhQzqscNcrWuGr7LHsWKVk5AzUvjOFIyUpxjhm4MMbgwpgDL1cvvF29i0Tvmb+HPz/f+zMd5nXg0PlDtJvbjl8H/Ur1stWvuNyxhGOM+mkUX+/5GoCuVbvy+e2f53uAc2GxWCz4e+R8UxcpNjw9ITQ0b23Dw7GEh+MGXKqPtDzgjFQtKkOfy3/p6JzrUQN4sPtl29580Zxbr1rq3/rmo21/508Oh4HdMLA7DLIdBg7DwMPFhqvtwlGIPs526dl2rjwY5pbLPpOYmMjHQcUg3Li5udGsWTMiIyOdY24cDgeRkZE88sgjl1wmIiKCyMjIXGNuli9fTkRExHWoWIojV5srAbYAAjwCzC7lmoT6hbJi8Aq6fdSN/ef20/bDtrx/6/vcUuviPwKGYfDe5vcY/dNoUrJSsFlsPNP+GZ7v9HyhHD4TEbFaLVixcLWz661WC15u1x477O55X9b0w1KjR49m8ODBNG/enJYtWzJt2jRSUlIYOnQoAIMGDSI0NJTJkycD8Pjjj9OxY0emTJnCzTffzIIFC9i4cSPvvvuumZshUqgq+1cmakgUN316E9tPb6f3573pWb0nw5sPp3ZgbZIyklh5dCWf7/yczac2AxBRKYI5t8yhQXADk6sXEbm+TA83AwYMIC4ujvHjxxMTE0Pjxo1ZtmyZc9BwdHQ0Vuvf3zjbtGnDZ599xnPPPcczzzxDzZo1+eabb/J0jRuR4qyib0V+f+B3xq8Yz7Q/pvHTwZ/46eBPF7XzcvXipS4v8WjLR3X4VURKJdOvc3O96QrFUhLsP7ufd9a/wy+Hf+F44nG8XL1oFNyIm2vezID6AwjyvvgCmCIixVmxuojf9aZwIyIiUvzk5/NbIwxFRESkRFG4ERERkRJF4UZERERKFIUbERERKVEUbkRERKREUbgRERGREkXhRkREREoUhRsREREpURRuREREpERRuBEREZESReFGREREShSFGxERESlRFG5ERESkRFG4ERERkRLFxewCrjfDMICcW6eLiIhI8XDhc/vC5/iVlLpwk5SUBEBYWJjJlYiIiEh+JSUl4e/vf8U2FiMvEagEcTgcnDx5El9fXywWi9nlFFmJiYmEhYVx7Ngx/Pz8zC5HLkP7qXjQfioetJ+KNsMwSEpKomLFilitVx5VU+p6bqxWK5UqVTK7jGLDz89Pb/JiQPupeNB+Kh60n4quq/XYXKABxSIiIlKiKNyIiIhIiaJwI5fk7u7OhAkTcHd3N7sUuQLtp+JB+6l40H4qOUrdgGIREREp2dRzIyIiIiWKwo2IiIiUKAo3IiIiUqIo3IiIiEiJonBTikyePJkWLVrg6+tLUFAQffr0Ye/evbnapKenM2LECMqVK4ePjw+33347p0+fztUmOjqam2++GS8vL4KCgvi///s/srOzr+emlGizZs2iYcOGzguJRURE8OOPPzqf1z4qml555RUsFgsjR450ztO+Mt/zzz+PxWLJNdWpU8f5vPZRyaRwU4qsXLmSESNG8Pvvv7N8+XKysrLo0aMHKSkpzjajRo3i+++/58svv2TlypWcPHmSfv36OZ+32+3cfPPNZGZmsnbtWubPn8+8efMYP368GZtUIlWqVIlXXnmFTZs2sXHjRrp06cJtt93Grl27AO2jomjDhg3MmTOHhg0b5pqvfVU03HDDDZw6dco5/fbbb87ntI9KKENKrdjYWAMwVq5caRiGYcTHxxuurq7Gl19+6WyzZ88eAzDWrVtnGIZhLF261LBarUZMTIyzzaxZsww/Pz8jIyPj+m5AKVKmTBnj/fff1z4qgpKSkoyaNWsay5cvNzp27Gg8/vjjhmHo/VRUTJgwwWjUqNEln9M+KrnUc1OKJSQkAFC2bFkANm3aRFZWFt26dXO2qVOnDpUrV2bdunUArFu3jgYNGhAcHOxs07NnTxITE509C1Jw7HY7CxYsICUlhYiICO2jImjEiBHcfPPNufYJ6P1UlOzfv5+KFStSrVo1Bg4cSHR0NKB9VJKVuhtnSg6Hw8HIkSNp27Yt9evXByAmJgY3NzcCAgJytQ0ODiYmJsbZ5p9v8gvPX3hOCsaOHTuIiIggPT0dHx8fFi9eTL169di6dav2URGyYMECNm/ezIYNGy56Tu+noqFVq1bMmzeP2rVrc+rUKSZOnEj79u3ZuXOn9lEJpnBTSo0YMYKdO3fmOvYsRUft2rXZunUrCQkJfPXVVwwePJiVK1eaXZb8w7Fjx3j88cdZvnw5Hh4eZpcjl3HTTTc5f27YsCGtWrWiSpUqfPHFF3h6eppYmRQmHZYqhR555BF++OEHVqxYQaVKlZzzQ0JCyMzMJD4+Plf706dPExIS4mzzv2cSXHh8oY38e25ubtSoUYNmzZoxefJkGjVqxFtvvaV9VIRs2rSJ2NhYmjZtiouLCy4uLqxcuZK3334bFxcXgoODta+KoICAAGrVqsWBAwf0firBFG5KEcMweOSRR1i8eDG//vorVatWzfV8s2bNcHV1JTIy0jlv7969REdHExERAUBERAQ7duwgNjbW2Wb58uX4+flRr16967MhpZDD4SAjI0P7qAjp2rUrO3bsYOvWrc6pefPmDBw40Pmz9lXRk5yczMGDB6lQoYLeTyWZ2SOa5foZPny44e/vb0RFRRmnTp1yTqmpqc42w4YNMypXrmz8+uuvxsaNG42IiAgjIiLC+Xx2drZRv359o0ePHsbWrVuNZcuWGeXLlzfGjh1rxiaVSE8//bSxcuVK4/Dhw8b27duNp59+2rBYLMbPP/9sGIb2UVH2z7OlDEP7qih44oknjKioKOPw4cPGmjVrjG7duhmBgYFGbGysYRjaRyWVwk0pAlxymjt3rrNNWlqa8fDDDxtlypQxvLy8jL59+xqnTp3KtZ4jR44YN910k+Hp6WkEBgYaTzzxhJGVlXWdt6bkuv/++40qVaoYbm5uRvny5Y2uXbs6g41haB8VZf8bbrSvzDdgwACjQoUKhpubmxEaGmoMGDDAOHDggPN57aOSyWIYhmFmz5GIiIhIQdKYGxERESlRFG5ERESkRFG4ERERkRJF4UZERERKFIUbERERKVEUbkRERKREUbgRERGREkXhRkREREoUhRsRkUIybtw4HnrooQJbX2ZmJuHh4WzcuLHA1ilSEinciJRwFovlitPzzz9vdokFLjw8nGnTpplaQ0xMDG+99RbPPvusc15KSgp33XUXFSpU4O677yY1NfWiZR599FGqVauGu7s7YWFh9O7d23ljRzc3N8aMGcNTTz11XbdFpLhRuBEp4U6dOuWcpk2bhp+fX655Y8aMMbvEPDEMg+zs7Ov6mpmZmde87Pvvv0+bNm2oUqWKc960adPw8fHh559/xtPTM1cAO3LkCM2aNePXX3/l9ddfZ8eOHSxbtozOnTszYsQIZ7uBAwfy22+/sWvXrmuuTaSkU7gRKeFCQkKck7+/PxaLJde8BQsWULduXTw8PKhTpw4zZ850LnvkyBEsFgtffPEF7du3x9PTkxYtWrBv3z42bNhA8+bN8fHx4aabbiIuLs653JAhQ+jTpw8TJ06kfPny+Pn5MWzYsFxhweFwMHnyZKpWrYqnpyeNGjXiq6++cj4fFRWFxWLhxx9/pFmzZri7u/Pbb79x8OBBbrvtNoKDg/Hx8aFFixb88ssvzuU6derE0aNHGTVqlLN3CuD555+ncePGuX4306ZNIzw8/KK6X3rpJSpWrEjt2rUBOHbsGHfeeScBAQGULVuW2267jSNHjlzx975gwQJ69+6da9758+epVasWDRo0oE6dOsTHxzufe/jhh7FYLKxfv57bb7+dWrVqccMNNzB69Gh+//13Z7syZcrQtm1bFixYcMXXFynNFG5ESrFPP/2U8ePH89JLL7Fnzx5efvllxo0bx/z583O1mzBhAs899xybN2/GxcWFe+65hyeffJK33nqL1atXc+DAAcaPH59rmcjISPbs2UNUVBSff/45ixYtYuLEic7nJ0+ezEcffcTs2bPZtWsXo0aN4t5772XlypW51vP000/zyiuvsGfPHho2bEhycjK9evUiMjKSLVu2cOONN9K7d2+io6MBWLRoEZUqVWLSpEnO3qn8iIyMZO/evSxfvpwffviBrKwsevbsia+vL6tXr2bNmjX4+Phw4403XrZn59y5c+zevZvmzZvnmv/II48wZ84cXF1dmTt3Lo8//riz/bJlyxgxYgTe3t4XrS8gICDX45YtW7J69ep8bZdIqWLyXclF5DqaO3eu4e/v73xcvXp147PPPsvV5oUXXjAiIiIMwzCMw4cPG4Dx/vvvO5///PPPDcCIjIx0zps8ebJRu3Zt5+PBgwcbZcuWNVJSUpzzZs2aZfj4+Bh2u91IT083vLy8jLVr1+Z67QceeMC4++67DcMwjBUrVhiA8c0331x1u2644QZj+vTpzsdVqlQxpk6dmqvNhAkTjEaNGuWaN3XqVKNKlSq56g4ODjYyMjKc8z7++GOjdu3ahsPhcM7LyMgwPD09jZ9++umS9WzZssUAjOjo6Iues9vtxqlTp3Kt748//jAAY9GiRVfdVsMwjLfeessIDw/PU1uR0sjF1GQlIqZJSUnh4MGDPPDAAzz44IPO+dnZ2fj7++dq27BhQ+fPwcHBADRo0CDXvNjY2FzLNGrUCC8vL+fjiIgIkpOTOXbsGMnJyaSmptK9e/dcy2RmZtKkSZNc8/639yM5OZnnn3+eJUuWcOrUKbKzs0lLS3P23PxbDRo0wM3Nzfl427ZtHDhwAF9f31zt0tPTOXjw4CXXkZaWBoCHh8dFz1mtVkJCQnLNMwwjXzV6enpeNBhZRP6mcCNSSiUnJwPw3nvv0apVq1zP2Wy2XI9dXV2dP18Yw/K/8xwOR75fe8mSJYSGhuZ6zt3dPdfj/z1MM2bMGJYvX84bb7xBjRo18PT0pH///lcd/Gu1Wi8KEVlZWRe1+9/XS05OplmzZnz66acXtS1fvvwlXyswMBDIGWNzuTb/VLNmTSwWC3/++edV20LOYay8rFektFK4ESmlgoODqVixIocOHWLgwIEFvv5t27aRlpaGp6cnAL///js+Pj6EhYVRtmxZ3N3diY6OpmPHjvla75o1axgyZAh9+/YFcsLH/w7udXNzw26355pXvnx5YmJiMAzDGdC2bt161ddr2rQpCxcuJCgoCD8/vzzVWL16dfz8/Ni9eze1atW6avuyZcvSs2dPZsyYwWOPPXZRwIqPj8817mbnzp0X9XCJyN80oFikFJs4cSKTJ0/m7bffZt++fezYsYO5c+fy5ptv/ut1Z2Zm8sADD7B7926WLl3KhAkTeOSRR7Barfj6+jJmzBhGjRrF/PnzOXjwIJs3b2b69OkXDWb+XzVr1mTRokVs3bqVbdu2cc8991zUaxQeHs6qVas4ceIEZ86cAXLOooqLi+O1117j4MGDzJgxgx9//PGq2zFw4EACAwO57bbbWL16NYcPHyYqKorHHnuM48ePX3IZq9VKt27d+O233/L424IZM2Zgt9tp2bIlX3/9Nfv372fPnj28/fbbRERE5Gq7evVqevToked1i5Q2Cjcipdh//vMf3n//febOnUuDBg3o2LEj8+bNo2rVqv963V27dqVmzZp06NCBAQMGcOutt+a6YOALL7zAuHHjmDx5MnXr1uXGG29kyZIlV33tN998kzJlytCmTRt69+5Nz549adq0aa42kyZN4siRI1SvXt15+KZu3brMnDmTGTNm0KhRI9avX5+na/x4eXmxatUqKleuTL9+/ahbty4PPPAA6enpV+zJ+c9//sOCBQvyfLiuWrVqbN68mc6dO/PEE09Qv359unfvTmRkJLNmzXK2W7duHQkJCfTv3z9P6xUpjSxGfkeyiYhcxZAhQ4iPj+ebb74xuxTTGIZBq1atGDVqFHfffXeBrXfAgAE0atSIZ555psDWKVLSqOdGRKQQWCwW3n333QK9qnJmZiYNGjRg1KhRBbZOkZJIPTciUuDUcyMiZlK4ERERkRJFh6VERESkRFG4ERERkRJF4UZERERKFIUbERERKVEUbkRERKREUbgRERGREkXhRkREREoUhRsREREpUf4fY6Az0b3ln0wAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from scipy.optimize import curve_fit\n", "from scipy.optimize import Bounds\n", "\n", "def fit_model_erf(T, A1, T1, S1, A2, T2, S2, R):\n", " return double_error_function(T, A1, T1, S1, A2, T2, S2, R)\n", "\n", "# setting \n", "initial_guess = [ 0.5, 300, 1,\n", " 0.5, 450, 1,\n", " data['norm. weight'].min()]\n", "\n", "bounds = Bounds([0, -np.inf, 0, 0, -np.inf, -np.inf, -np.inf], # lower\n", " [1, np.inf, np.inf, 1, np.inf, np.inf, 1]) #upper\n", "\n", "\n", "\n", "popt_erf, pcov_erf = curve_fit(fit_model_erf, data['Sample Temp.'], data['norm. weight'], p0=initial_guess, bounds=bounds)\n", "# Extract fitted parameters\n", "A1_erf, T1_erf, S1_erf, A2_erf, T2_erf, S2_erf, R_erf = popt_erf\n", "\n", "digits = 2\n", "print(\"fitted parameters:\")\n", "for param, value in zip(['A1', 'T1', 'S1', 'A2', 'T2', 'S2', 'R'], popt_erf):\n", " print(f\"{param} = {value:.{digits}f}\")\n", "\n", "# Plot the original data and the fitted curve with residuals\n", "fig, ax = plt.subplots()\n", "\n", "ax.plot(data['Sample Temp.'], data['norm. weight'], label='Norm. Weight')\n", "\n", "# Plot fitted curve\n", "ax.plot(data['Sample Temp.'], fit_model_erf(data['Sample Temp.'], *popt_erf), 'r--', label='Fit')\n", "\n", "# Calculate and plot residuals on secondary axis\n", "residuals_erf = data['norm. weight'] - fit_model_erf(data['Sample Temp.'], *popt_erf)\n", "ax.plot(data['Sample Temp.'], residuals_erf, 'g-', label='Residual')\n", "\n", "\n", "ax.set_xlim(data['Sample Temp.'].min(), data['Sample Temp.'].max())\n", "ax.set_xlabel('Temperature (°C)')\n", "ax.set_ylabel('Norm. Weight (a.u.)')\n", "ax.legend(frameon=False)\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "ae421f7c", "metadata": {}, "source": [ "If you would like to have thermogram fitting implemented in the package itself, please bring it up as a suggestion in the issues." ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.12" } }, "nbformat": 4, "nbformat_minor": 5 }