From 70b90f47a58be4ae626da247b71be6e449cd0c74 Mon Sep 17 00:00:00 2001 From: jonhood11 Date: Thu, 28 Jan 2021 20:38:47 -0500 Subject: [PATCH] Update twoDGaussianTesting.ipynb --- .../twoDGaussianTesting.ipynb | 181 +++++++++++++++--- 1 file changed, 151 insertions(+), 30 deletions(-) diff --git a/custom_2dgaussianfit/twoDGaussianTesting.ipynb b/custom_2dgaussianfit/twoDGaussianTesting.ipynb index 7d3ea91..3d3aa28 100644 --- a/custom_2dgaussianfit/twoDGaussianTesting.ipynb +++ b/custom_2dgaussianfit/twoDGaussianTesting.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "code", - "execution_count": 9, + "execution_count": 68, "metadata": {}, "outputs": [], "source": [ @@ -39,7 +39,7 @@ }, { "cell_type": "code", - "execution_count": 42, + "execution_count": 78, "metadata": { "scrolled": false }, @@ -50,19 +50,10 @@ "text": [ "(160, 150)\n", "(160, 150)\n", - "(160, 150)\n" + "(160, 150)\n", + "(50, 80)\n" ] }, - { - "data": { - "text/plain": [ - "" - ] - }, - "execution_count": 42, - "metadata": {}, - "output_type": "execute_result" - }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPQAAAD4CAYAAADb7cuFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAABLS0lEQVR4nO29b6x131Hf9521zr3PYwMuuEBqGbc2EUFyo7ZYFlRNi5DcJECpnapt5CipnEAVVUpSaIqCHb9I3iA5pU3Li7aRm9C6rRuaOKD4RZLaoKZRpEJiOza2cRzMn4DBYEJSgWL/nnvOXtMXM7PWrLX/nP/33HvumkfnuWf/Oefsf58935k1ey1iZnTr1u06LFx6A7p163Y660B363ZF1oHu1u2KrAPdrdsVWQe6W7crstWlNwAAbukZP8eXXHozunV7FPZb+Gf/hJm/amrZgwD6Ob4E30RvuvRmdOv2KOzH+H3/eG5Zl9zdul2RdaC7dbsi60B363ZF1oHu1u2KrAPdrdsVWQe6W7crsg50t25XZB3obt2uyDrQ3bpdkXWgu3W7IutAd+t2RdaB7tbtiqwD3a3bFVkHulu3K7IOdLduV2Qd6G7drsg60N26XZF1oLt1uyLrQHfrdkXWge7W7YqsA92t2xVZB7pbtyuyrUAT0Q8R0eeJ6BMTy76XiJiIvtLNewcRfYaIPk1Ev/fUG9ytW7d528VD/y8AvrWdSUSvAfC7Afyim/d6AG8F8K/qZ/4HIoon2dJu3bptta1AM/PfAfBPJxb9twD+FAA/wPRbAPwwM79g5p8H8BkA33iKDe3Wrdt2OyiGJqI3A/hlZv5Ys+jVAH7JTX9W5019xx8log8R0YfWeHHIZnTr1q2xvYfCIaKXA3gngN8ztXhiHk/MAzO/G8C7AeAV9MrJdbp167afHTK21W8H8DoAHyMiAPgaAB8hom+EeOTXuHW/BsCvHLuR3bp12832ltzM/HFm/mpmfi0zvxYC8RuY+VcBvB/AW4noGRG9DsDXAfh7J93ibt26zdouzVZ/GcD/C+DrieizRPRdc+sy8ycB/BUAPw3gbwH4Y8w8nGpju3XrtmxbJTcz/4Ety1/bTH8/gO8/brO6det2iPVKsW7drsg60N26XZF1oLt1uyLrQHfrdkXWge7W7YqsA92t2xVZB7pbtyuyDnS3bldkHehu3a7IOtDdul2RdaC7dbsi60B363ZF1oHu1u2KrAPdrdsVWQe6W7crsg50t25XZB3obt2uyDrQ3bpdkR00FA4R/QAR/UMi+iki+lEi+nK3rA+F063bhezQoXA+COB3MvO/BuAfAXgH0IfC6dbt0nbQUDjM/AFm3ujkT0D63wb6UDjdul3UThFDfyeAv6nv+1A43bpd0I4CmojeCWAD4L02a2K12aFwmPmNzPzGGzw7ZjO6deumdshQOAAAInobgO8A8CZmNmj7UDjdul3QDh198lsBfB+ANzPzF9yiPhROt24XtK0eWofC+RYAX0lEnwXwZyBZ7WcAPqgD1v0EM/9nzPxJIrKhcDboQ+Fcn9FEVMV98NCHYsQP4GS8gl7J30RvuvRmPG2bAvXU9gCutWuwH+P3fZiZ3zi17OAYutsjtPuA9tDf77CfxDrQ12yXBngfa7e1A36QdaCvzR4TxEvm96PDvbN1oK/BrgXiOetw72wd6Mdo9wUwHVF3xOl02+Gtw71oHejHZOcA+RhoD/3eU8He4R5ZB/ox2ClBPhfAx2zDKQC3Y/TEwe5AP2Q7BcgPAeBtdkrAnzjYHeiHaMeAfCTAFE4v6zntCdcpAH+iYHegH5IdCvKBEJ8D3l1/Zy/I/f7tC/cTA7sD/VDsEJj3BPm+AN7F2m3ZGXDb50PAfgJQd6AvbfuCvAfEewN8znh7C4B7A36I134C3roDfSk7A8g7A3yJRNncb87AuBfg+3rtKwa7A30J2wfmU4C8D8DnkOW7wGi2BfCTg31lUHeg79MeGsj3FVNP/c4cmFugPDnYV+atO9D3ZbvCvAXCRYiXPrsFXrqHctLq2fttkG/x3P44zMK9L9hXAHUH+tx2bpAPgHgneMMJ4uzUQDjxu7OQt5AuJMG2eu1dwb4Cb92BPqddAuZ9IT4FuHM29d0LkO8F9wScFGi7FL9yb33oUDivJKIPEtHP6N+vcMv6UDjASWCmQGOYKZSXt0AjmIkov8p6oX6123yu19zvu22Y3N6ZfZs7DpPHrP3cLvZIH0k9dCictwP4cWb+OgA/rtN9KBxgfAHPrjcBZV60ALI3u9DduiMopuCZg24KtlO9lkCf2Mbpm9F4f+eOzSLYC8e+Xm/Hc/mA7KChcCBD3rxH378HwO9z85/uUDhHeuW9Qc4/uwCx37YliCZAOtVr6fe2evGp/Zs4BnPHaivYu9gjgvrQGPq3MfPnAICZP0dEX63zXw0Z68pscSgcAH8UAJ7j5QduxgOyXb3y7KKJz0+BXP1kK00nJPTS8qnvmPmtg0zj2fY3cqzstyelenvbdTT2tu8q30HVb8lK41h5Nr6+srj61EmxqatgdigcAO8GpBvfE2/H/doRMJ8c5C0QT8an237/UIuYTl41s5h5O9z7gj2R2Z7Nhu+TBX/gUB8K9K8R0avUO78KwOd1/tMbCudAmI8GeQ7iJYCnJOrC7zVfNL/MbPJi1xSKh6gBnUDVcm6VhM+MN+CfHOxHDvWht+L3A3ibvn8bgL/u5j+doXDuCebJWLTdhrmYs00k+exwaJbPJsti/bm5V4gL39H+VpOldssm4++pYz6xz1PH76jzMGUPOKY+dCicdwH4K0T0XQB+EcB/DABPaiicbSf1RBJ7H688e0H776/mb7/o94+l40xZZ+P5/E8lN8Gp8rLZA3vJbdvt42wnxStvfUhs/Yg9dR8K5xA7FczHgNy037afrzxf+9n2t+d+t7WlZTPX0ej68vB4aPx6to5frvOq7zMpPjVvtG6zHQ2wswUp28C+AD99KJxT2jlgPhXISxBPLKt+a4eM+N7mYlwAcvFHH/O6+JpQ4Mk/7bx9UA+cGR4nytrkWRVfT8XWTVx9UBb8gXnqDvQ+dgmYl6T1HMhTnliXTUK8lCXfNt9be2HHOJ5n0th/zH4+S2+fBIMCb/sxA3YrxXeR4VcIdQd6V7skzK1XbqV165H3gXjRS+8bP7v1MzRNvOshb+GOVDx3BiuN420XZ1Oa8NYG2BOEugO9ix0A8ylB3tkjO7jzZ+YAnoy7lxJme3ro0ExPxbMt7CaTVZpXcFdyvJbivl2bW2+9iwRvmraOata6sHWgt9l9wLwQK0/CvA3kJYhHn52GmEf7sL/kpjQDNHE9L7FCHCfh5griJs5OIctwy4qPMuITsfXB3noJ6vYGdQHrQB9jO7VZ7pFNXoLZe+WJRFflkds6af/bfl5uEjqBl2aulzODo7/ATSJzgSzDBlnul3kA7SsDxnG2nxdSgbqNrdubzY5QT9oD9tQd6CXbs4DgYM/cyO6tXnlfkGkMM7fzR9txvIcuzU/lr4e8gtsy4MzC6hzYJsU92BPeOifNtkjwbVDv1OVRaxeMpzvQc7an1D4aZu+Vp2LlRl6PQG5ldeOJec5DL0hwrgCfPgyVOadF2fOiAlrW4yKlDWao97bP2fxWjk+CPeOtIWDPSvBj5PcDTZJ1oKfswIz25PKpTPMEzDt55TZO9s8Z22enQF6C2Dw2UQ3tlNfeZg5ghr+gCUgKOSAxdBs/A2CdTz5pRg5sYAT2KL723lqz4ZPe+gCoR/YApXcHel/blgSbK7PM8w6HedIre4+c129AngG+gridn7e33dmJfRpJbb3xNImwPJloDLeT3yOw83Y42AxGuKy1eWs/bU1cPmG2D9R+1w9Jkt2zl+5At7bkkQ6EeUlmLya+2lh5yivHWm6PQJ6D2MDPv9/sf16/2eGpm1RzkedWoEpqOy8ddBlzgXtCZldg2zJT2LZPzKCUioRfkuDNzWAnqCeatB4y1B1ob/cF81wWe84rtwmvWH9XFR+H+ns4hJEXbqEGFFxLAGWodzguQK7gzNI5/+e8dGK57m09TX5R4AK3QU+pkuTspmlITrITMKRyPFNSJWAedpwwq+Q3IGDvKb8fMtQd6ANtl2FnlrLZe8M85XG3eGQmksyxQRvdeigQM9EIXqZ6e3dLirn1mYU7i6WJwFDAE0BMAp8up6FIcEaQ9XJiLeSbAMdQe+wYShzuZbj31l6Ct01bts9zUM/Y1h5GL2QdaLM9vfPs8jYJtovMnsti7wKzLW9BVg+cPXRehhriFmBTmK2X9vOmDkGV1dYPZi9LNeBB4mlSeU1MWV4jEQhJPmpy3MfZzGDU1WFZhtvx3gfqbfL7VEmye/LSHWjg5FJ7b5m9LV6ODurWK7tl3iNnqEOYhNgDbPNbiL0Mz9Nzh8kvTN47A4gCb4FYL26egNuY8HJ84BpsYBxf2zKT4C3UBE2Q4WRQH/ws9RmtA71kp4LZVl+C2Wexo+v5o/XKHm710DW8Hmro+piEWG4MDtRAldSW937n5w9V5XyYcnVn9swW+irkGcQpuAflwmAmBqdUZ8dd4owGB5BJcEuYDUPzNNc9QD1n9+ClO9DbEj5+1bm4eaqt2X/3NpndNkm1MFvia8oru7jY3osUNy9dIJXl3iOjVG5l4OE8NNWSeyHyqCQwo1y4TAJpEDCIkYFjJlBQmAcW0Blgi5mTOlYGCAHW9GXvzSNzDAqvz26Xc2C7wAHAgBJabIN66jzv0px1QendgT7UluJqX7kVGi8+FTMD9fptTDwBc46NWw8di7S2z3mQOaCGPGfAkbPlHFBLcF1u5tV1fs4iAr5/V2KoJ0YBnDUxlgrcnHQdQgGbUJqmgiTMGMVD5/cpyfFQkDmGMdQGKTfedQLqfPwN6qUk2QMsKgGOBJqI/gsA/ynklH0cwB8B8HIA/yeA1wL4BQC/n5n/2VFbeS7bI3beW2q797mccy5mts7upjxwm/hSaCuvrHEyx1BgVbCzrI4qwXXaAK7gNfntZbiX33n/3fsWYkBlNuXlArfKZYU8t0QxgTQxRgODg0yLx9abAHGR4SkVbx1CyY7nhJmT80J0kd8p6c0ga3DdwJCbtKqHOuwcHiK9L+SlDwaaiF4N4D8H8Hpm/qJ2DvhWAK+HDJPzLiJ6O2SYnO87ydbel+3yFBWwexKsbZpqPzuTyc4wtxLbxcocqI6TVxMgG+QWO0fUHjuDTDXQWXrbcVk4FqXVSb/AeWn10MRUgLZlg25XEigEZlQeW95abKxyWxNljFRLcKBkwYO6aZcBJytAAeR3fAHJFNTb4uklu4AXP1ZyrwC8jIjWEM/8KwDeAeklFJBhcv42HiLQh8bOc7AvJcHa3yUH+1TM3HrmVmK7WHnklUMB2GAu8lukcQacmveEWoK7GHqU4SZU3hlQpgCF2LLXloSmPG3FHyK1GRj0dxKDiWSXByfFEwS+fHyt5BMgBAHYLw+hQJ2ldVkun1WIc5Ks7MxkG/WcPbAE2cFAM/MvE9F/DenG94sAPsDMHyCiuWFyKqOHOhTOoVJbJnTZQhJsrp25TYC5JqsRzK3EDiVOTpHKvKDAmke2eeaNA9y089iwaezmpSvvDAW69dIl6UWJ9CkoldpBZfIgn8nt4+atYQ5bPDgrocSUPXSOq4HcdJW9d0poY2oC6nbqpSQZsOylD5XeZ7BjJPdXQAanex2A/w/AXyWiP7Tr5y86FM6cd972SKTZDnFzXm8JZmuecgC3MptDUM9boE6rULyu98pxDHIVJ4ca4jqGLsmw6Vi67JZPllUPYHiwHcwwiF2iLMfNCfmFWN6bhzYZToO2OZu3VgmeYaYIoiRVZkAtvy1rbjE1sBPU+XyeS3qfwUsfI7n/XQA/z8y/DgBE9CMA/i3MD5PzMGwPqV1/bktcPRc327Ilz7wLzCazzQOrxE4xCKAmrydgrkC2aSrTJUlWe2ssQG3G2YXrKhbKNl7ZirUK0KTemsu0SXAiJ8U5/zCTHNuApMVoJr81SYYgVWYIJQMOjOU307ynnkqS2W89Aul9DNC/CODfJKKXQyT3mwB8CMA/hwyP8y7Uw+Rchy15Zz+/Wb+eF8qNZappysfMmvwymM0bw4PqYE4ragCeew9NkhWoq0ITL7lb2T1hVTdhOXamEkOnkuFGADDofIuZifUvgUjQZJL1A0psDWj8DeVHN4oG6JhZBAQuUOfjiSK//bkbmoFdprzurtL7AdgxMfRPEtH7AHwEMuzNP4BI6C/FxDA5D952iZ1bOCey2nm9JandNE9NxswtxJrJ5hhqab0KtcSOtVdOsQE5x9MFcj99LNBwQFd/LW42rxwEVs6Sm7IUDwMhEXIXRQkKLEGnUwa/ynglAJEEXDj5DQApCdRMVdt1TpLNSO+lrHd1vi1TfuFY+qgsNzP/GchYV95eQLz1w7ND5PZU+edUIszm7xo3N/XZ1TwXMwvAoSS6zANn0BuvHAu8KZr0RgW0FaKUaUxDDffeHzr/nuv3OW4GmpjZv1d5nUSW53hZl3EAwsAFcj2MlICwUYh985Y1W+W7SompYQADtXc2qHdIks1J712eytrJTii7e6UYsOydzbYkwiabqOz7fNxs32VSMIY6braikaC/EajKZNtfyWY7mNUTJwM4Fg8NBzjCGPBx/FzmoQV7bh8V3txsZd6YoZ6x9cgKeRJpHbRCjMl8LmUZHnIbGSGtRIJXd5WYN9HFEAIzR9K2bkL+ZK4msx90UE8PZz6b9S7XwMPw0h3oY2wqXp7wzrJs7JFz3Kx/OUNOuYlqlMF2MJdlDuYR2GVennYxdFIYfEyd/xLvJLt9+3MVP3uonadFC/ZQA5xz0ZocE+dFsASZNK1JXB3sd80sprb2adZYPDjIgyXgqJxDi6cDLWe9bbtO7aVPZE8H6EOaqtplS81UrVefi5tDAzM5mLPUnoA5Fpg5OohXjcxe1eDm6QruGuqSFGP3Xv8COwEtC+UlnpJHMtueuELjpSmYpBbQArE88TjIPEAz4UPGNMtvzXFLFqfZQELQZFgCD6jORY6nA42LTgxkPxJyIClHXUqQ+XO/r5c+kex+OkAfY1MSvG2mAmrvDHcBTTxwkee7JFj+62FehQK0ZbdX3iNTSXxFKOBFdicPceOdZZpHWW+Du02M7SK520RYabbiKk7G4Lx0Bhqa/JK6biKV5RDZTNqkFaDJMgLCWrscWkGLUAKwSW67qG7OsgUuAz7/IAcVL80+s92s56+TC2e8nwbQJ/DO2eaaqWZqtfPv+7i5kdq5eUr/5lpse0rKZbXrZqkCdc5oxxIrp2jwT0PNkV1dt0xX3lmhrj30+IIlfYrKgGaFmdVLyzPOejOyohIqibGc3NK/IrlLbJ24SHBmRoqU4+i0AgKXmJoSy01Rt0cSlQxrzoJ5ZZPe+S7F1Y04P2pZXQ/D2Bu3sbRf/wJe+mkAfYy1yTCZmFwnL6Mw9s552YTUbuJmL8vLqwHby24Ha57vXkmhLd6bnYd2Ujs4iP37vG/NtBpbBwV5BoTK5OAOpYAEWkDClvEmIGyUtUE+byFuAiEw57g6RZtWqLnE1IDeCI1NPWajeLo9D3YOnaQmQB4I8bXeCx74obRLXz/QhzRVzdlSZnspETaT1a7LOpu4WT1xiYn1fW6mkos7rYoXTpHAKwU3YAJszhDDeegCuHlkdYkGsM4jW+ZNQbZHIYuXVnBZPfVARWIPqobMWwc5fjRA2p/tMGaPTRJXZ/5IPLeiygypKgMjrSh34Y2NblNiUNB4mgf4pixTEbnZ0JqyphJkVlfuM95T7dLH2JFe+vqBnrM9Ckm2eucpee69sv6dldr2rLJ55Tj+W0HsYmbzwJYoy17Y5q8U3AAnvzn/hQFNACKLLA0Kb2CNW1HB7I9AvvT0aaoCM4GteyFtmjLPzS5mhkrxRCrDCRnI5A6TZcED9CuL3wXYNWWx7DuAUjmmoFKUeCI/nRXcNOu2UvO3vWZ2yXjv0oR1Jnu6QB9jYQZgn9mek9o6naV25bmd1CaUJFhbHGIJsOx9aeyR8zRX89JK4+bAAnSQ9wYz5Re0GYmzdybnndvrnTPQpO/FI7OBTVTAttdAYGKEDcmVqDePBELY+O8GLN3OkbXks8hxDnoTgcXslCUy2R1B1YOUkDrpPTo/nEOl/Ox0m/E+hZ2pTboD3douTVU237c71wvL35lmKp8Iq55HrhJhmIihC9jwsbQDPHmZXclvBq/UK5s39iBH+RtI/6qHpgZoDzZnyS1AG8xJQU4sfzko2IkUZGQZK6WdZJWdCBZDDy581Z+UJiv3+5Yc1Bw22weidldkI1oGk95AlSDTc5UTZNaM1Xppn/G26rGlQhN/PT2GxycfhR2S3Tbb0lQ1+Z3mEfK6DcQap7F+T0mEhZzVznGze++bp9KKKlmdY2YP8MrFyya5I8ArBdpk9YpBMcmmh4QQE4KCHEJCpAJ0cEAHB3RyQCcH9CaFDHZKhDQEpBTEaw8BvFHANSOPDRRwhVPpTlly6ykxNZtZtJ5ARe6yyWxAM+AAWdY7Bi2CUcGeEsji5vZmawUqlvFe8tJLsru5zs6d7b5uoM9kkzXbZeFIbk/FzkVyo8lq63eQl9pU/c0yvJXaOdPNDnC9yLN3lhfFhLAScCkwYkyIMSEQI4aEGOSvwCwPSQZ931pi9cT6npmwSgFDIgwpIDFhCIxhYPHaAUhE4CFIe7R+j8hhRmB9IMOeukoAVgo1I3tdKyu1ZqwcT1vvokyQvsfKo5kmtVlP32Qzlnnn9mmsLMk1lm4LTR6APXmg964MW/quVhG0oLuYuUhuKt0F2bST2vWTUXYBO8irdmXnjV0WGyuNlyMDqwSKjBAZIQ6IUb1xTFg5kKN66BhSBjlgXnInCNQF4IQhBQyJMSQr62SkFMQbrwISJe0soVRnM9SLZc/PmlyDhAqsTs3Gh7Z9Z+Qui/KNL0tuQDrtJznfDNX1LmYeJT6L1JYpnXaed65Neqly7Nz29IDepQPAfeW2/8xsVVgDMrnEF3lwLbM97Z1zXfZEAixV8pvBK4xgDjcqq9Ujr1YDIjFWMWEVB9yE4qVvwlBgVqADjS/MxCEDba91ihhSwCYFDEyIQ8SGGAMnEEVsNuLpmEgrLxVqKgkukd/yNgDARvfNkmLaJs2aCEuRtTqtSG/W5J55abZY2t1gyc5fcmCr167O71xebI/KsXM/sHG9QM/Fz/vatsowzW6Tv8s7kH0yzBeRlGec22YqB63zxj7R5eNla5qS2LnEyslgXiVgpQmvVUJciVeOMeEmDrhZDQqvTocBq5CwooSVvheQBeY4IbkHJoGaCRv7myI2LECvU8Q6JKxjwGaIiMQIIWIYGIOpEBLJbR4b/mEMZqcEkJuFgeKxRY4L1MGazjQjngcTiJAilxhAnKTohKmq8+YgZaflRhxmCk0ko7dz5di+dmAcfb1An9OmPLjZlOxuZXb22GgeWZR5tbduJTfVMnvUPOVi5hwvK8wxYbUq8vpmNeA2CrQ3ccCzuKlAvg0bgQ+MVRgyzAHlQkvaKj0oxAmEgQl3aYVNCthQlJsCMYLCTBT10JiOTpDuteWbWbsukgOhcXfUSuyk+8tAeeJM1q+OUS4Y8V4ZObFdVYw5eZ1fPru9dB08oN5KgCcO9M7FJNvMy22bNu/sl482wHlui6UtOaZg59rqEdgOZnvAgjzcdSa7ktkO5mcrgfgmDnge19k734QhA70ikd4RaZQUsyz3gIANpQx0JMaaIjYh4W4QgAMKzOZxsUJu7kpIYA5A0vZfk83avRCzNEexFqfovUBktc1zN8hyHFmlrkGMGkafHFNZnpuwLPFllWM7NkcvFZmc05400DubFRp4yH2pp63T/nUy3BeS5DiZULU310UkTezs26EnYmepBGPX1switdUjh5hwczPgdrXBjUL9spt19sq3cYPncY1nQWU3Dbgh8dQG8o12fh1dHD1IITXWHJECCdgpYB02eJFW2KSI2xDw0pBwRyvEJL8NADEE0CaiWASQ5CZheTGr/GKCDbfDA8r7rFrMS5vHlliVI/Lz2VJMo96XrMDEFZrYebMsdzvd9uU9Vwq67XI6Y9XYsUPhfDmAvwjgd0J0y3cC+DQuPRTOnHfddUQMs6kk2ESpJ03B7LelymrD3QyKt2i989gLz7zPZZ0lfs5NUwpz1L83ccDtasCzOGSJ/TyucRvk/cviGrdhgxuF+VnYIFLCDQ0SPy8kxQYQ1mklHjoEvEgrPAsbrDniLq3U0yesUsRLww0A4MVmlTPm2WMjglk99UbyA5blZllRbmxFrWtSTJ2gvirYnZfOvYrqDZK03jxnvb0Mr/4quL43halSULtu9omjT9gefayH/kEAf4uZ/yMiuoWMnvGn8diHwtnV2gcxgDw9anv2TVaA6xnEVU1liU2oml/c+7ZSrJbjmg7Wqq+glV8xSjPUKtbJr+dxjedxg2dhI0CHO5HjNOB5WGewg0I9JbkBkd0DAtY0IHHAmiMiJbyUbvTRRj00KO3YyZ5h3qwwxIQhJY2NrdKMSylpYJXcdkOsoS3Hx0IPyplt/4CJ3TTJStGA3NrQJjVz2/URsvsSdkxH+68A8M0A/jAAMPMdgDsieguAb9HV3oOHOhTOlM01V+1qUx66+n7U2e6cBFMpHlF32JdjQRcfBrfcnpIykLVG2wpHiDzM4pU9zLdhyDDfBgVbPXQBeoNAjFvaIGAhy42AO15JkxWvsOaIABa4UXufxIRNECpSlJh7SKQJMakqk2o29ZzRYluNhXWZdQnsY+h8vFJJftlyubGyHvciu+G9rgFs53HqQQ0K2KVq7BJ2jIf+WgC/DuB/JqJ/HcCHAXw3gJ2GwnlQNiehW6uSaBOy2s/3TVfVe9RgkoOYvBd2VWFtEcnsQxguqx0ZcTVgtRKZvYqSBHvZSmT187jBl6zu8LJ4l0F+ebjLID8La9zSJnvmGxoQKSFgQnIjYFDPPED+3vEKNzRgzRE3dINAjBc05AQbAIRBLj9ru87fp3BLSXWQeNMexdQsdwZYXylKsxVyDywaQ/tjmaxSDA5u5Ioxi6fH57MBvJ0/Y5dIjB0D9ArAGwD8Ce2j+wch8nonu/TYVrPD3EyZT4g1D2OMEmUys/k8kIePCQVyy277XkIMeu+pS4cD3mNz+euz2lYF5qq/VlHiZvPEzzVefhbk7/OwxjPa4OXxBZ7T2nnoDW7VS0ewemhNaIEx5CargEQBdxzVU0fcsNwI1hxHcXfigE0MbloqzCQ2JgwxIWmdd0oM3rCroOOS+HL5htJEpdxlVaPg6nxYtttiZi0oyWD7TTWAXZFJFUcHmn9Y40J2DNCfBfBZZv5JnX4fBOidhsK56NhWp7Jt7dFTdd7I10b20HWM56Y9wG5+3TSDyuPbM8z2kEUMnEs4VyFpu7O8nrk4+XlY4zmts4cWsBugKWX5XAFNAYkDAhISQvHiAYha9bEOEc91fZHlUkG24pDbwaVkNORtZ3vMKpjc9vkFBdDf7FyYQs0x9dOkaghaoz7ywm2V2NR5PzKOPlem+5iRM36ViH6JiL6emT8N6Vz/p/X1NlxqKJx9Yt4l27F+e+r3uQXdsttwnjrD67y1l94NuB5g3zunJYmklxF7llnanc07r0JpV15RwjOLncNGvLO+DGoB/A7PaY0IzkDnZisnNQ3qNSJuaIM1r2R9TngpSds1grRTmwnQUavLAjZhgzVFqf+OA8ImIgRGCqnEzLmFwCcUC8REmLg5Uk6K2TGVw6+JMj0fZLWlVtutVj2ssa9dqGLs2Cz3nwDwXs1w/xyAPwI5NA9vKJy5Jqul+Hnp5lANLRvquMs+a6DqdLng6mKSKQ9dvDNNFpXkLoRcRwUU1StraWdwJZ23FkPH9ShufqYgvzyI5L5VD/2c1giUcGvt0OadXWJsUE+2Jmk/vtNMd4SogzuOCJykvTqoZ9d4O5eLxoC7NOR68JV6a8t2UxRv7avD5Ekrk99eassNIEvwGY/tE2I5yw3UcXRzLqunr06dGDtRfH3sUDgfBfDGiUVvOuZ7H5Xtogh8/AwAzb2lSGoayfBqGm6eey/fyU5u2wMW5SGLFQ25AszamL28fh4E5Fsnt0WOJ9xodvtGvVd0uzzYPJXgAYy1k+OAJM2eh7U8iUXyXoAmJBBepBVWlLChlB/dDPk5bEayYW5yuz2DByqSW4/FqFKsuoFqHG3nwboS9advh0TXvnbfT171SrFt5qX3LtnwNkPaCoOR5C6LuJlXZ8BRA+6kp8lK0p5GiJAfgfRPS2Wwta3ZmqWmXwXmW0oIAG50u3xtV44uWGJ1ew5y0CBzAGEA4YYibikikbyPlPITXDca0684YJ3cI5uuSyTOYUWRzu0xIUx75Pr4N39tJ6biWTuXFZAHSvB7sg70vtYmuvz8mWluJBy3n58Ae8pDTwHOrodO6wcsaMcEuTqLUqnRziCXl/fMOYbW1w0YkYAbAJEo358iCAMYN5DEcARjDUgmvPFyAwfc0oBEawwUqt9ehwErLQddJdnGaPuQAqyTwgKzPfpYHxNqjpmHuVoOD79rh9YwqOr0YMaR5ky33bD1ySs575fNdHegzbY1Y+3yhNUW+c0V1HASmyYh9vPLchtvyiXG1DPDvBqQZatkt4vczt45iCeus9pDhvkGCTdg3BBwqyDfqEKJSsYNiuReIyECuNOe9gcZ9BUgiMROUWHe4IYinoV1jqVvdPtWYUAMMYcLIy+tstsSg5QoS/GqNQBjyOskmd4Y3HmhbVJ7p8z35T13B3oH2+vJqxnzF9rSOj4LPt4Q93nzWu6ilXsEVzK7vOS55qhtyVHLMG1aPKskwCIYt5RGMEcQgvroSISBWfv7SgAFrDnhVj1egnvyCgk3GLCmhFsa8BKxJM2sTlybxKa2mxx4IyUDd6OrjmMBlLfExeatF+0UMvueHrV8kkDvVVQyZzt47FZql+XlLU949xrasj430/772t45S+d+Ci45aFybcjCYHeQBXCpJIR45ahwMIEMtqtSeW1SPDM6fC2DcUMKaBd7AqfnduuMED7EcEtfbqN9vnyCj+jhMxcge2sobz6xzqD2EkSifJNDntFEbNFASY36RNVm18+ds5IXgvHS5iFqQBU7O8EQq8NrfWxpyO7PMl7g5koEcYP9kd6odQSBgzZo044REElNHMBLr9yrU+S/KtkxtJzmwbb/yFLWQNvfNKTVkctsd/5FE1rbo/N1LyuyBFpd0oM9lh8j0GW886bEnvLRdsKVfhVrKSmGmzSueWdYtPXtF9dABQCCR2ZIQky+OPtvP1jZN2re99PMlnp1zjB0czME8dC4jrZ/i8r2Lkr9h5fcLYcnUMW1sH2+cuyR6JHZgOVS3bIdWlD0gs1LNNjs9ZbFpuvPTYeZyWvre0Cybet4awDyBjWw+SDdP5itOEJZdwB7/1fiUbSmM3+HjU49CAtMXxRysuy7fxdpnrU+R6nhq1oF+zLbgUHcRiVa2mfQysOquKR+ZdO7QVDnZdJr8FKrvnbLklsljlL7m2624ww4R77beTvaAi0eWrAN9rM0VERxyQbgLMl+cPuflL9j26920gWCA+I78pF46ZHCsTzCbTg7qBCAxZ1iT/sjAKb/8/ISk68tnB1D5vonfy/NQtrFsd6ulgVGb3y4A73IaFtahB9ar5zbrQN+HGfMn8x7+PVUXfO6/urn4DZLBe8DGM1uZZutRBzAGtrZldgCzvhIGLgmw+rM0eiVIH2QDS8cIKSsFqjw2u/2p97dOWI+OyZTlroH1s1OgJlt+PxA/qMcnn7zpCIetUdKKJsiFwXNyM2l3OoB6mpK9JZ64PrmWlMS1CMhDuDIqqA0SD8rADij7S+q5kRQ8ATRoXUayLkJy0io4Gc7ZO68V/jVbB/zyvQXeULw0qIbZbXcZ1dKBPQXxlIJp5/vzsw3W+qDOrzcB46XboIEnCjQrTDuvzww6tOzARiW05xYaEKkBOU83F2dZ389X2tj90Yt/0NEfNzYcTSBsOGI9eq2kDpsH3JGIcmu6ilrGmbR9ORA7z18aYQdtb16DcaeSewBhjYA7BNxpg9maV7mLojVHvaEE2a4UdcSNUMbH0n0oNz2331ZrskV2kw1cN3lu6haxs9o9SfcnCfSk6QDhi8vjzDIj1PfnvIspvHLRuWpiLnCT90x28UHAJR3XCYmy9GYm7QlE4DCpvUkCzSYFrClinVYYgtZS00a7DYq4xYB1BnoASDwwCAjMlac3GyAQ36lnTgDWHPQVtY+xchMZELBO0pGgDZVj25hDA9dpIKqXQUw13NweN1THrdws5bh7gHf22nt67EtYB/pUNiPBAZPeRYJXAFfxMMbeJl+MyN3a5qFUGfI9zCpNi0RNzmNvOCJhyN3tmhQ3L2leeuA17igCDC0uYZRIW0o6kxWtoKQGJH5GhvnOwTzkDgTFO99pj6B5tEq4MbFU/nulkfcpGbg1wO1xGknuKZnuzkt7nqrz6W0O2AcCslkHepulVIpHOGmncBBVzDz/4EbSK8pGLM9SupbYXk4St/COs93tRZtfiWQMqPwS6brhgJhsNIsgnpqKt3yJbxBZ6qvtPUiLTVhi2xtKWsIJrLUKLO+m/l2zymyF+SVe4Q4RL6UbvMQ38lvpZiT5NylizSGHBia5/WDxeThZOw5pfLyoPUZqU8es9t6N903uPZr3eVbzGd+U94g7CewGjE+4uZQpM1mvF6SsLxeoQMxOSsp8tou0PAOR+6MmHTUCrN+ZCJwCUmJsEmOVAtZDxIoS7tIKq5TwIg25HNRqtwH5vqAd+j2nNRIFeX4ZA9YccEMJd9aJAdpeP83j6wgaeqN4iW9w596/xDd4keT1UrrBi7TCi7SSge00jh7sRpQCUgpgGR2+hBUOZhUQY1DtmOZlnI9hPg9AuRvNedmp87hn4muUKHvA3fg+LpO+W7fOX/S6W3+DSwydk2GsoJqHRpmftEDfeYzW47SeJcPN0LGQkbvT4USVZxuSxNPeS29YklAvEuOGbuTxRdZRLlDGafVNWgLv4GR4DbRIZ5XwkI72K5gVYPPSHuYXGtffDTGPJT04pcGsHlobtymR7HcqN8IW7iVVQ+74y3my5azg8/imvG/2+qAahNOAfjTQJMMJfgjALzPzdxDRK3Hpsa1OZUtJLpZ+OgCIzPJdYsx8vkhsgTt75RkJLhebeO0sw9Ujk3Yib2MzsUpFGeupZLrzi0RuvxhW+ammdRDAAOkQYc3yVJQ36T5owMDy1+LrtVtuY1uVDPYqwyyjaKyyd/ZS278GDtX2ine2/bJ9phpM77Ez3FyB62+gdk4K4LyYEKPkwGau5fQ9jCR5iJ3CQ383gE8BeIVOvx2XHNtq30zznPnYec/fN9iqbTKZrID6GNmyssykUlF64zBgKXun8mK7gBOk88kAYNA4mgnDIP1br4cgozxSlJEfUyyPVg6ln22zAUH7/gpIYY07jrjV5qxb7cbXP0Bh1V6+k/2EkGPnO15VMvsL6RZfHG7wxeFG5XbEXYpYDxHrFLAeDGjKUGPwHlphbY5HjqG9BE9Oarvj7NVQdX6WwiVvuyTCThVL7+ntjx198msA/HsAvh/An9TZbwHwLfr+PXigY1vt1RadEjgEaYu2xFgk2DjGWaL7TDdzlRijoNI7QXrnVE9iMTAN0G5oJQNeS2sF2EFMiYBB9oEigEHD7yGAA0s1VmBsQsRd9ThlfYEE4lwSCmhiiyLuaJVhfsl6GmF93NJ9RSkSUZg5qHe+xZpj9spfSLf4wnCLF2mFLw4C9Z2ORPliWGGdAjaDgD0MAcMQkAYCDxJD0+BfJr3L8ci5hcpb63lw6yCxgs3lHEFvtD6mnsl4cztfPTW3HnyLPcjhZAH8dwD+FIAvc/N2Gtvq0kPhHGRzbdV2cVhXNXPxl2+uyh4blbwmrmWlz+jmBJDz0mhfg4xVnChhGAhDDNgMMtC6dRpYdUs03FS10wMH3IQNntMGCMCak3SgT6vc0wjQxtBBm6fKQHUis1dYpxVe8CrDbHGzeOaVeGZ9bQaR2sNgcbP3zuSgdIDmG5wpGztuRXrbzRNo4m046NtzufR+6rp4IHbM6JPfAeDzzPxhIvqWfT//oIbCscSYArtTYkyLLRaXu8IRuYNbE5YtL7EzmHN5o12cHGxaPjuW4NqmPZBsS5A7v8EgcGgPmoM0W8UU8wP71jWRWSAGAkTuUsjxs4xPJSNQgsvz00BJntnIkxZH+3i5zWjn16ByOwVsBm2yUqh5IPBAZT8NanbeOTU3vEZ2V/J7VGxSIDe5TQs34snzO3fe81t9b8DfQ9x9jIf+XQDeTETfDuA5gFcQ0f+OHce2ejS2LTGW5bcmxhIA4lLTbbJbLxjOco9y7Gbw+niwig9Ncg9y36HoJPhQ8nEcSLZlA0l7EWMgiZ85Uu6aSHarZKjXIWITog5LE3LPoDYk7E3wY1uNL2SfGLNB3y2rvUkBL1IdM7803ODFZpWl9t1GPfQmYtgYzEG880Cgjcrs9pXK+zqm5vo4DnAxNYMGdp67yWyr3M4JsdzEpTA6mf0Q7Zixrd4B4B0AoB76e5n5DxHRD+CSY1vJxk1DONd0tff3uww3sOzRfcGCFZjYIGwD9CZQuqTNF2WgnByzixRENcSDrhcYGEhqra2vsgFIFEAB2GzE+wvQ+vCDq2rzlVprDrhJ0qWuje0c0q0bTrYGOsfeeXgb8dIv0ioXjdyllcTMQ5Qms80KX9zcVHGzwTwMAWkTgI3CbDHzYPtOdaysrzCUY+dVToZY4+eijgq8OSnmz9fEOZxsUz6F3J67QRzQ/HWOduh34SGObXWI+Ux3SkCMrjiEx6Webbzl42qbn6QzdysYyXG1NU/5C5Lqi1aAJzftlmk7NBR4BgAipIFg46SGEHA3tZuRcMsbAFIksqGElT7WaH2RrUMZuN0n1vyz1psUcilnC7TIa4mZfRLsbmNJMMpSW+Jmaaai/H6sWrxHbgEfZ72LV/ZVZVWSf+r87RA3jxJij72jfWb+25BsNpj5N/BIxraqMt1LcbSCOfLEuWE4LsdUSTuK1wvKZLd8vNQms7sIS7a7QIwR4CYp1btDYbYm8YHAG414V8BmIzeJu03tWKyuOnFAigXoTYrS+T0xQiodDraWPwstZHEgSzlnxEvDSrz4EPFCQd4MUhm22USkIZastnpmmJQeyn7m/R9aqA1kzpC3We9SfMIjuU2t3J45j9LQf1z8fK4MN/CUKsXOaeatzaNr3Jzbo7Oso0JultB6AblYWCQkIUFuIBQ0DiRI/BcoPxxRuv6WgeKwIQFYY3SmBOYg15W2f69WlJ85HlLATQq4CRHrOOBO26lvwwYrHWNKMuMpJ9C87C6SW5+W0pvD3bCSqrQUMshrLUXdDKHI7IEE5rV65w2BVHLTAIQ1IQwAbQhhA9BGjkXIsTSPY+rkoVZIB3bTGMvt1iN7ST7ncR9QdtvseoFeSmbtY02BCVs8Fpxnbh/UqC6OMq94ZVXB7L22dhdLciH7jDYl/W5LiqkXJjcNyE9RkH6zE3TIHM2Uc+6QIGLDrqySBykRjUm8aQwIxFhTxE0cqo75V6Fc2PWz0dBHNUMGej1EbLjUZW+0aCQZ1JuApNNpQ+BNkdoYIPAOLmcwkQgbe2mu3+f1OWe67TxIy0E5LwYyedAd3DwCPeXrYXS9nMIOiJ+BawZ6znZJjM21N7fVY7YeO7gt1gayfKuy3W2RiWbGKbB6UFfdRC4x5ry3zatiaJXYMsaVyXBIbA3kzHeygFtmAQCCeusUdegaYgwhYcMBK0ql4/40LkyRw0D5xZpgM5AHVQHWLJVSwGYTMGyi3FQ0LMhx88YVkGR4yWWq5wB283KMbAArqEPtnXMxiffGet6qc+6vncnr5QgPfuKM+dMDurGdK8a2efypKrFIozt9bpMe1OvaunpexTvrBUwquRO0w3kaeWUK9UCJFETWiyRXqBlgu8dAnsbiCNgzxyEEpJgwpDI+80aHo7UhanyTl/1tO/aznkYYcHXZrqZ80JJOVwVmMjtntDf6amCWMMR750ZqT3hnqHfOmW7GvHduvfJCtdjs9ZHfzq97zvgZ6EDXNvfk1ZLsbp+PBkqbtJNxbJ6RDGYGpySRL2kMDSjA4mWxEc9MA2cZTSTzgsl3IqsmdS1jhJRY55VKK5OOPDA4EniVkFQxDAMhxqBD0QIhlGFdZQTdZjSLfCxKe7Z1rJDUM6cUwAwt41TYp7yyxc4DEKzNeSN5hBwzu79hw7XsHlhhV6g3JX7OWfBBbpwivxXkweUzDHQPdkrzcjuZytohu32P7dbXDfScVz2l7PbfaUUm1XeglnLeS1vbM6vcVulN1j6d26VVXrcS2+Jo/cnsqanE0YHME5dKNI6QvxxEoUQpbDG1EgKDKGAI0tRlI0B6uU3EGWbZ1dJlkD1YIXE69BltBdnKOR3MpXCEMpAhtz8XmNs4OnjPXclwzrGzyWpfSLLonQ3afAnVMfPONuGJR975hO3PZtcN9I62JLunmq+aD1cgZy/druulXFU55uYZsNpsJhe83ntUfgdNnBWJzZoRR4Y+jz0OQiKWs8yivZn0Ucyo85jlQouMITBCVO+uRSrtiJbA9G61vY7aU19sPY4kK+U0mE1Om0cmB2rxzL6ZKlgG280vBSWc16uftDKJ7TxxC28LUCu/zebA3lFu34d1oFub895NkQmHoLLbwQ4H91QpaNWEBYggVhmumpmhkG6AtEKGHNoslYicBDev7KahX83yCiAkZilsYxlIjoOCFYts50GhHgTgZAOtk24/nNT2iTG2G5n1kWZeWSFmOIi1zd1KOg3MjTRNITUy2zzyBllSh2qZg3lQKW4eerDmwAKxzUdKWXpn1bSt1LOJpx+i3AaeMtBtvHysl54qBW2ldkLOeBMJhtlbD/UFA4jUzNCuAIMaJqdB+XmPxIzAJZ7mWK7PoB0FcBSZzAHgQb10JPmADW1K1tzF1ZC1BvbkAykMlGe8Kd9MStdBGjZoLF/Ja/O+G/PabfLLJDfn+Rlug9jDPHB53wBeYmYsx8762qnUc8o771JMcibQrx/obdnpyc8o7HNxdP5qLl7a4mXz0sTFq5uktu2pYmm5OAV0yLYO4hFZvXHOekP2hZJ54nKjseng5DiQQ/Eshc1bW4cKrB47934SSLZFs+XqmnXwAJoHGqhh1t5EyP6azOYCcwHYPXyRCsRINbyjJFgTN1dZ7Qwpild2x302dvae13vnNq72dspSzyMl+/UDfUozQO2kL/RoMlloohls3y4tRSKcIdZGaYGZymfEKH8G0Gz3IDIcUKhJPLcPAS0JBrY2Z5X/2i8vBYCDZL6t7VqzYLJtCvVW098wiLNnZoPVQ23wOpjZe2RbXsMcPMwGb37aqrRB5+nBgERRQkveWc/X1lg4paPhO4c9baCXZPdcbfcuTVi+vtt7abgbAjMwJCCG0oylTVeVR4eCms+UxtJAbpYKEC+bk+qJkSIpwOpk7GcTJNscVJZHAZYDgTfyY2zNcK43FjbJLZvgjqH9paazBsCeX84Qm1dOLn52srutAguurblKiGWoWeNtAVPkN4M2yWW79WaZUpHcBvmSd7ZrwGe50wTouv5DkNvAUwH6ENm9y3dOeWmjSpNgVbu0JcgGHe4tBGBIICJpkw6hkuDywL18SdgIgHIpKMSWuWZCiiWGli55tcNAbQ7joNdn1PdU4IZ56OAgJpQYGvV73YQCs66Q70EGtEHMBq56aAOXG4ht2YYryMPGZbN9m3OGm523buJmVphNamePnZz3djCbd26TZNU5nvHOxxaNnMDjPw2gj7V9vLR5fG3OyoDYTaWV3vIFkjxCKTShIYFjKFVNarlz0Q0BK87wBd2OqkmLRYLn2DZvMKz6UxJ1AUBQ6U0CtzWBGcj2AMjUJUf+e0cwOy+dCuBhqOe1EjvHz8k3TbmYOT+E0UCcGFV5Z5sI81lt1PPsXI6ap7YUkly6qcrb0wF6zkvvIruXrPXSUwUmdiE4L53nA+KlgZL1dvE0I+QqMgCl6ESfrDL5bYkoIvHUJqkDnPfVSq4Mrnpri3E52DIqkhuyvRns2eNrIE9JbwXRgT6W3FzD7Us6M7wK96ZMV0mwXA2WqiYr+f00L7XtPI7kdAP2rnaI3D7RTeHpAL1kuzRhzXlpfb/ULp3BcAAjhuyFi9coLcryF8hJMsBlt5PWiQNAaVvmyACoAJyc3Nb+xqxW3JquLGZmat7LVyvIlKfzMSL4ML8oAOehAQfyKI52wLaAVwB72GuZjaRye1OKRzzMFje3MGepnf8OjdRu4G69s3s/ip2ry+r+YmezDvQ2m/DSi+3S0vBbJ8h8sUkeJ0tOdo6nAZHiwBjqjSbPtGok16fE4vAt6SWlnb7dGUhRmomyd87xsjZXGbgeaH2N4mjb1OYwTcXP9j5DXHlorpaFYQpw55V1Gnm6kdm7wOzj5jazbeduIhHWnPzqOpi8Xi5oTwvoQ5Njc+3SU9Vjo6ITiUerBzcs653j2Dqezu3T5rE19uaESn4DCiNIo1t9zBL23dqcpR7cCk7seewMNpVukZBKnM4+jsaWQ8e7Ql2k9agJK5VlwZ5hdnCXpqgZmc0oEE+0N1dNVFOJMDtf7a7NxM5l2j57IMwnjMGP6cb3NQD+VwD/EsRRvJuZf/DRDoWzR+XY+LPlRjFuxtJ0spfePp725aCAeBPr5CAKAPZUllyxViuGzK09CgmWJ6x4pXI7uMqw5B6Z1PgYOWbmHDNb/Oy9s5fg88dAV/Gxc37PJdttNw0XV2cp7b33RPycYR5cBntwntlntLWYZGsTVSu1gflmqj3Au4TcBo7z0BsA/yUzf4SIvgzAh4nogwD+MHDBoXC22Ym89CiWds1Ys23TU/G0amZSkDlKUxYCScY7JZBCSZGQ7wKMDLgUkElhSLIqMEt62XuLqYODNzi4s+yW7crzASe554+bH/wNcLGz7mMlt7lAWpq0uCwzic21xM4PXWwMYAlZaEgu8VWgLpnuIrmrroem4uYpmH3CbId2573sxBnyY7rx/RwAGyHjt4joUwBejUcyFM6k7eqll5qx5rLeLdR+PSswgY+pJWDObdOpSZQxifffJFAkeeAC5vAlY0V2nZpnTiV2luaqkjwzoA1i32RVgJ6/+HwirGq6Gv11ktuBTSbHvbzOAHPpmMDL7fwXY5mdeBwzV5I7zcNcneux1F5KhMkpv4x3Bk4UQxPRawF8A4CfxGMYCudYL73D91Z13u13tFVk9ln/Pcn0uLZNxwZqvzpQz0sMRCqxNRMQ5SYkNdzaDKXdBjNxSY45uHPeb8cYGpiOo0eS20B13jmXhyrUWU47qW1xsJ+uYmaT2anI6cmX97hTsO0rtQ/t0fPE3hk4zXCyXwrgrwH4Hmb+zV3HVn5QQ+F4O9ZLt6C2WW/tWL+Kp7VirPLSMeiFJ78tBSfI2W8ElrgaQH4ggljludZyqwzPg8Nr27RA7OS3Nk2xxdPyg02Ge/t5HY+5XAAXoLmR3g5kgzo5kL1XZjmWxCK5wc4z544M0nLz1FTcDBwntbfZY3p8kohuIDC/l5l/RGc/jqFwlrz0rsUmW6BmhXprPO2hJspt1AJZ2UYyHRwZeawrJikdjZaShia7uAIbbduzJbyiwuqKRwAHeU6I7XABey9t8Op8k88AkDshMChtWWpATm4+L3hlD20qCbKtMB8bN9v14K6T+ngswHwG7wwcl+UmAH8JwKeY+c+7Re/HpYfCOYPtXEF2SDzdrutOdn6GuoIdlQTP3ppYJDlrt8Cq7kkT7QisBSkEezSSTARUyTHSGFqb2tx9jyfugVMFJl5uVzXdOl1B7kDOjzsa1IzSJJVc8msK5vyZeZl9aNw8affUef4+duxgdf8JgI8T0Ud13p/GYxoKZw8vPWutl26/u60i2wVqJ78BJ8EBV2WGLMHBBEbS30+wvrhzk5g+Fpmfiya5OeUqsApuKR8tZZ+6O2TSf96mMt1TEGfvneNoLus18jonvloP3Xrg1jM3y7bCPFUN5i+HU0rtM3ln4Lgs99/F/Pl9FEPh7GsHSe+pJNmuUJNmsoOWiQaXqQqUJThHkmeQibQirZHhSbaLdWC7IrEljib9Pl9IUnlmzbjt02w1Kb8Tjz13joUNRNTyeuG9T36VLoTSKJu9M8z5RB8ptZfsjDADT61SbMr29NKHxtN1+/Q81GADzL5X/6MCLQPI9Fl7NUGkN8uA7/aeLA4OSqoNx0MKt2a1MXAdSwM5Q5+z3Jg5Vj5EsLcWL3uA4d5XSTLXDKXTlRfOsfYY2FEpp0EOiDqa631kCuZTxM3N8vu2DjRwnPReqvU+AGrtd7eJxZHLRAE03lqbtsxbexlukFuyLLE8kglAn63McOdeSmBeWjy7f8KK1VNP2XQczfX7BuK8zGLkFt5JsCe8sl82JZuPhbk930t2Ialt1oHexfZsyqrnjZNkO8vvQTvdjtazAZemrQTx6Na0Zd6aqCTKHOQ2TZyyxy5tzQ3ckCw5gEpm79IODThvbNe2h1jXzfHxAsiyn6mGdSr5NTgY9ckp+dkdZfa2JNgDLCCZsw602ZKXnlp9V+kNFKid156FmgeFMKichlyk9nktC60kuHlrtvkGKOUHL0bTA7RzhNI0JYmyrK/ljzsmUxlus9pDF3BlGTcSnJ03xLS0br3uFMjt8rl42c6Le3pqyTPXy/eQ2kt2D94Z6EDXdoZ4GnBQ+AtnFmqV3FPJMkAlfPHIS2D7xBk0qVXDDQCsHfmX5BmAURwNIPdE2hq1F6uH1++3yW7MQGzrHAiyfHwB5rw5EzC3y9x+yIIj4uZ7ghnoQO9n+8TTrfxu26QPgbr9nlFszbNgZ8/dxtZABbw9nFHURX2Tm+6EyB8DHr9vgbW/3hvneQsgT62zD8xzMntyehrm6X2+vNQ260C3doz0lhnLmW9gXn4DzWOXDmpwkeA+2dZAC2j7L3gMNlC8dkKW4NDnMP1808vs3svnmwMwcS3TFNRT3jnPb6YNZDc9AtmOA7dlnC5etukpmZ1/e3eZLZPNDW0r7PfnnYEO9LRtk97AcpJs15ja3gPjMtFBf8PiaqIabNvGKbCBSorL9zu4dbnshs5PqPY5F5K0AA9bjp27vufBnvDWcxDbMdwV5Lz+gldupneR2TL5sGEGOtDzdiZPDQCTxSe7SnAAVZv1HNgKPjTxVclx2Qjdztp7m5WbD+r1tx2zuenWC7u/NOGh0cIN1CDn7/Te1sGty88G8za7AMxAB/pw25Yka9dZatKagxooEhwoT2txmTdq4qridJTyTyK9ARTg2W4UVSyt5gHON5QdrIJixkMDY2/s12kBn4qTgXmYlyS2n57d7h1gfkBxs7cO9JJt89IHQj0pv6egtuKOUDxz9rytt8YE2OaxgSLH3T5RI7+5veHs46WXvDMagKf+TnnwFuS83rJXrtc/g2d+gFLbrAO9zc4INeDkN7AcVwPzYGu8bfF1luLyA9NwA7ViIH06q93XPcKOsr8zcE9BtCvEwElArtbx3+2/K08+LpiBDvRudiao5av39NYA8ggdiV3W2cOOaa8N1N8/+P3isbSmPaT26Jg0n5uS4jPwLYJs8+dAnvi+pwIz0IHe3S4BtU3DeetJGY7dwHbfJ59xv2OA23v5tFt3B089CcFMTO23Yx+Q/fxtXrn5nWuHGehA72fngBpYluDeW2NBhlszV/a06pW9HNffyr8zaBtUlQAL44vTw77Npi7sCXjlbQMqsB1iN2+ySGQfr9x+Lx43zEAHen87FdTAbt4aKLH1EtjA2Gu3cAPTgPubydA0NJ8ihsaO3nEJ8DR1I9gD5KXfzbMeN8xAB/owOxBqAPt7a2Asw6fAzlKc6htG/l6/Nd7buQIZ8Fhat9uyYLOjMC55xW2g7wpyM/9orzyx3nj5w4IZ6ECfz2bqvhclOLCbt/Y3FJc4A2Y8NlCBOwt3u57f7rl4elvBxRQUu3jrKZCBrV55vP6ZYH6gdjagiehbAfwggAjgLzLzu871Wxcx7y1n12mkdZ69gwQH5r01MH7Qw77Kza/gxlBAzPE2gCa2lm5BnU3t37byz0M8dbNsK8TN/EV53f5WnnWAxG634YHZWYAmogjgvwfwuwF8FsDfJ6L3M/NPn+P3Lmrb5DewnwSXmfJ3V7AnpDjQwO2bvYDaewMN5NWW5m3YakveegqUOS8MzEPcLms/eyjIc9s4Wufhwgycz0N/I4DPMPPPAQAR/TBkiJzrAxo4GGqZPdH7yRYZDuwAdrPMX/QjuIE69m7Nx+LbbAmKie9YBFlWmF12Mq88s+54nYcNM3A+oF8N4Jfc9GcBfJNf4aJD4ZzD7gNqYD+wZYVR23MLwiTg3lpvvmRbwJ9MnC1BPLH8EJBl9nXDDJwP6KkruzoifigcIvr1H+P3/XMA/+RM27OPfSUO3Y5jzvk4Lj18O05rfTtqewjb8a/MLTgX0J8F8Bo3/TUAfmVuZWb+KiL6EDO/8Uzbs7P17ejb8Ri2Y84W+tM5yv4+gK8jotcR0S2At0KGyOnWrdsZ7Swempk3RPTHAfxfkGarH2LmT57jt7p161bsbO3QzPw3APyNPT7y7nNty57Wt6O2vh21PZTtmDSaLdfr1q3bo7NzxdDdunW7gHWgu3W7Irs40ET0rUT0aSL6DBG9/R5/9zVE9H8T0aeI6JNE9N06/88S0S8T0Uf19e33sC2/QEQf19/7kM57JRF9kIh+Rv9+xZm34evdPn+UiH6TiL7nPo4HEf0QEX2eiD7h5s3uPxG9Q6+XTxPR7z3zdvwAEf1DIvopIvpRIvpynf9aIvqiOy5/4VTbcZQx88VekAz4zwL4WgC3AD4G4PX39NuvAvAGff9lAP4RgNcD+LMAvveej8MvAPjKZt5/BeDt+v7tAP7cPZ+XX4UUMJz9eAD4ZgBvAPCJbfuv5+hjAJ4BeJ1eP/GM2/F7AKz0/Z9z2/Fav95DeV3aQ+eab2a+A2A132c3Zv4cM39E3/8WgE9BSlYfir0FwHv0/XsA/L57/O03AfhZZv7H9/FjzPx3APzTZvbc/r8FwA8z8wtm/nkAn4FcR2fZDmb+ADNvdPInIEVSD9YuDfRUzfe9Q0VErwXwDQB+Umf9cZVYP3RuqavGAD5ARB/WGncA+G3M/DlAbj4AvvoetsPsrQD+spu+7+MBzO//Ja+Z7wTwN93064joHxDR/0NE/849bcOiXRrorTXfZ98Aoi8F8NcAfA8z/yaA/xHAbwfwbwD4HID/5h4243cx8xsAfBuAP0ZE33wPvzlpWtn3ZgB/VWdd4ngs2UWuGSJ6J4ANgPfqrM8B+JeZ+RsA/EkA/wcRveLc27HNLg30XjXfpzYiuoHA/F5m/hEAYOZfY+aBmROA/wknknNLxsy/on8/D+BH9Td/jYhepdv5KgCfP/d2qH0bgI8w86/pNt378VCb2/97v2aI6G0AvgPAH2QNoFXy/4a+/zAklv8d59yOXezSQF+s5pvkucO/BOBTzPzn3fxXudX+AwCfaD974u34EiL6MnsPScJ8AnIc3qarvQ3AXz/ndjj7A3By+76Ph7O5/X8/gLcS0TMieh2ArwPw9861EdrzzvcBeDMzf8HN/yrtyANE9LW6HT93ru3Y2S6dlQPw7ZAM888CeOc9/u6/DZFqPwXgo/r6dgD/G4CP6/z3A3jVmbfjayFZ248B+KQdAwD/IoAfB/Az+veV93BMXg7gNwD8C27e2Y8H5AbyOQBriAf+rqX9B/BOvV4+DeDbzrwdn4HE7HaN/AVd9z/U8/UxAB8B8O/f17W79Oqln926XZFdWnJ369bthNaB7tbtiqwD3a3bFVkHulu3K7IOdLduV2Qd6G7drsg60N26XZH9/829XDslevYyAAAAAElFTkSuQmCC\n", @@ -89,7 +80,10 @@ "img = img.reshape(x.shape) # undo the ravel \n", "\n", "print(img.shape)\n", - "plt.imshow(img, origin='lower')" + "plt.imshow(img, origin='lower')\n", + "\n", + "min_index = np.unravel_index(img.argmax(), img.shape)\n", + "print(min_index)" ] }, { @@ -101,7 +95,7 @@ }, { "cell_type": "code", - "execution_count": 43, + "execution_count": 71, "metadata": {}, "outputs": [], "source": [ @@ -123,7 +117,7 @@ }, { "cell_type": "code", - "execution_count": 44, + "execution_count": 79, "metadata": { "scrolled": false }, @@ -151,7 +145,7 @@ }, { "cell_type": "code", - "execution_count": 45, + "execution_count": 73, "metadata": { "scrolled": true }, @@ -188,7 +182,7 @@ }, { "cell_type": "code", - "execution_count": 54, + "execution_count": 65, "metadata": {}, "outputs": [ { @@ -235,6 +229,13 @@ "needs_background": "light" }, "output_type": "display_data" + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "(22, 25)\n" + ] } ], "source": [ @@ -256,24 +257,31 @@ "ratio = (img+1)/(bcg+1)\n", "plt.imshow(ratio, cmap=plt.cm.jet, origin='lower') # set axis manually: vmin=0, vmax=2\n", "plt.colorbar()\n", - "plt.show()" + "plt.show()\n", + "\n", + "min_index = np.unravel_index(ratio.argmin(), ratio.shape)\n", + "print(min_index)" ] }, { "cell_type": "code", - "execution_count": 56, + "execution_count": 63, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ + "(22, 25)\n", "[-0.22368031 24.56356777 24.5694857 8.58063247 10.80474537 0.99137914]\n", "[0.00427531 0.16302866 0.20599005 0.20294417 0.25014993 0.00165347]\n" ] } ], "source": [ + "\n", + "\n", + "\n", "initial_guess = (-0.2, 30, 30, 10, 10, 1) # amplitude, xo, yo, sigma_x, sigma_y, offset\n", "optimum, cov = fit_twoD_Gaussian(ratio, initial_guess)\n", "print(optimum)\n", @@ -282,21 +290,20 @@ }, { "cell_type": "code", - "execution_count": 108, + "execution_count": 96, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ - "[-0.22368031 24.56356777 24.5694857 8.58063247 10.80474537 0.99137914]\n", - "[1. 0.8824969 0.60653066 0.32465247 0.13533528]\n", - "[0.76769883 0.79398196 0.85571018 0.91876078 0.96110731]\n" + "[0.135 0.325 0.606 0.882 1. ]\n", + "[10.135 10.325 10.606 10.882 11. ]\n" ] }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAS4AAAD4CAYAAABSUAvFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAABwjUlEQVR4nO2dd3yUVfbGv3cmvZAQCB0JJRIpUowiIoIoiEqxrqLY666uuurP3ruuujbWulZsa1uxFxAQFCQoSO+9SiAJpGfm/v6YCZlzIpnADMkkvM/nM5/Mk/d979xM3rlz73PPeY6x1uLAgQMHDQmu+u6AAwcOHOwtnIHLgQMHDQ7OwOXAgYMGB2fgcuDAQYODM3A5cOCgwSGqLl+sWXNj22eYqhdfq3Y0W6oLPIqXKJ4naVmnmv8cg3w9t5UvUGgSBU/eWij49hape2xvG83FsVZslm1tkm0RL2lualPBm1XsELwkKmaPrw0Qu7FcNpgkaUUT+R0VVeiVJ8jmKY6Old0tLpXH4+UFsbZM9k81b1T3iuP03yMRt1O250mWZ7jL1b1ToRoItlke8OcUN5V9iS+Sr71tkby0eUfJd6UlCJ60tUieIN/K6n+svm31fV6quLqVbKeq52vWwrZtVr/CXqGLMbYo+GkAbIJvrLXDQ3m9fUGdDlztMwyTcqr+i2lXq//QDeqCXMUXKz5B0jX/bUpNiEXekMmlOwWfGXuo4EPG/Sz4+KuOU+1V3VGvcKk4djOPyrYekm1xiKRvnDpU8PO3/1fwJWmtBXerUT3z7vWywQGS/jFMfrjSZ+6SJ7SXdG6bgwTvNX+Z4At6tBM8o3S14DElcuSK3iTbX5Alr9d/T9akNYIXDJK3apMtaiTcLmm1gUwj4M9ZcKbsS/eclYK/fLi89LIHJZ82Jkvwo8f9Kk9QA121T10LxecrvlLxHElL3q16PmAgIaMIuKKW596D+sauI9TpwOXAgYPIhyHyB4ZI758DBw7qGC6qKRkRhzoduKJWWdLOCVgequVJtem91rjeVbyHpB2+/UPwGcN6Cb6KDMGzY2cLPuRmuZxb82i64P2YKfgceu9+3pkV4tj3HC940W3yVhjxwSTB9VLp7bTTBU9GLmvd6s3KPFUtFa+TNL2FXBqWy9UN0fKtoHkbtU5Xq/ru8+X6ZUYP+V73K5wr+KSs/oInIFWUbTQTPHVInuCtvs0XvLyf7E+0khFsb8k3pqUJ/lvvqhOi1Hsfky0lhcu+UO/tZZIe3VMuDT+9apjgo9/9Vl6gNCpWSVp+s+TRv6nz5Soed+CtEIZEGANEh97MfoUz43LgwIGAs1R04MBBg4Mz43LgwEGDQ0OYcZm6dIfI7mJszpMBv0ja46k+qG3iT3tI7UDrPm3YKLjWTcpUQM1G2gjeki2CD505TfDv+h0teBFVIQbdWCiOfcIpNfZV6ypJ6vhqtYfek98F36qC3lqovo/YKDW0khRByUuUv2g1XmpI1fREFb7x6ePDqAk5HCb4ubwj+PU8Kfh5vCn4mOmfygZfqbk/nCrpG5l/Ebw3Uij6hSqRLFUFBB7MEsHn0EfwC16UoSo6Tmvlha0E77RRxvRpTaukt+QzE48QvKvqz/ucJXgGq3c/vyF7Bstz8kOK4+pkjL2/lueOhdnW2uxQXm9fEOkDqwMHDuoYzq6iAwcOGhwcjcuBAwcNEpE+MNRt/ywi76pEZtAQd53k21+OE1zrSLfysODfl8rYqSdjrxc8npozsKYi8yXi+8nzh74rNa/yk6qev58i465O5X+C67irF7lS8Pasq5H34xfBA2PIoLou82sbKQIVqcl/dqFKS1GxRU9/cbngJ/CN4Md6fhD8d3dPwY9nouBtPFJ/vMt9n+DrVFDfrwNk/7scuVzwJpeqlJ9Bkuq4uMkcK/gZfLj7udZCe62S6U2ejupjonMJZYornaYrTUtKqcwdkCl4NQ0tR2po47PlvfW3/JcEj95a9fx+nde4D2g0My5jzGpgJ76Q0AprbbYxJg14H8gAVgN/sdbu2FMbDhw4aBhoCLuKe2Nrc6y1tnfADsItwERrbSYw0c8dOHDQwFEpztfmUV8IxY9rNPCG//kboPb/HThw0CBRuVSszaO+UNsZoQW+NcZY4EVr7UtAS2vtJgBr7SZjjDbnAMAYczlwOcBBLRGxWXMSZX5b+5elrvOD0iV0fttAfhS8RexWwX9TOlDHgHgXqK6r9GSe4B9xhuBfjjlZ8IcWV0W7/Jgi9bEuSE1mJjK5Tv8tzdgm+MmeLwXf5pY6TD+PzJts8q7UfFaPlXFggyZIjUxJaNU0Iq1p5SBDdcau/Ujwo+cpzexgSWdlycTSIx+SuYxHDpD8tkF3Cn6Ve5zgO15zC95hmcxTzVDBUnmkCr6Erruf6/i/goPkR7LvdGXIpWLiWCvp+LFSkxr7inyvml8q80AHI/VC/b859wR5vXlCvX5g7mIZIaMhLBVr278B1tqN/sHpO2OMdsbaI/yD3EsA2V2NUwvNgYMIR0MQ52u1VLTWbvT/3Ap8AhwBbDHGtAbw/9y65xYcOHDQUFA546rNI2hbxrxqjNlqjNH2iJXHs4wxPxtjSo0xN9a2j0EHLmNMojEmufI5MAyfR+ME4AL/aRcAn/55Cw4cOGhICLPG9TpQk7XzduAa4PG96WNtBs2WwCfGmMrz37HWfm2MmQX81xhzCb5V/plBWyqEQEurbYOkbhOoOwCc6vlE8CbLVOyOip/Z0F56LrWdL/18F/ToJPgOpXvo/EDtI3883wv+btbo3c91LtztPCD4B+rt0flmbZDexsn58m/9PU0GA7ndMk6pidJ4Br2rNC2d2zdVcZU/l3WTtE7uHCt5Ne+0OMWXSto5S2p+yFAmNgyS/zv9Xredo7yZlcUVyl9sU6Z8v3R7WYur/p6Xs8aKY5etHS8bkyF3cJriSh8cmyM1qZIx8njbSepv0RblUs7DXKqO609toDRbTMgwhG/H0Fo71RiTUcPxrcBWY8zJezrnzxB04LLWrgR6/cnvc4Hjql/hwIGDhgwDRNdW/a6guTEm0AX/Jb+uvV8R6ZsHDhw4qGMYA1G1H7i2Oe4QDhw4qHcYA9Hu4OfVJ+p04PK0MhTcUPWSJ+RLz6joh9QF8jDlUqZgZ4oUVpoVSu3g1x5S2FlBZ8FP8MhYpXfc5wreR+lWOlfyS6qSFQ9TNaOu4t+Cv85Fgm9Rflo6ly4vLVXw0VOkqGNlamC1fLhqflrKx1yFrMF5km5uIYOVdBxU1jqleek7SabfkTZJJvhZJTLo92PIRun/P6m39KxXIXoMuV2ef+bUz+X1j8vrs7ZU9X9w1mTZmNKYuE/xOYqrMnnfPSV9244q/Ume8K4qOqn0Oa5X/APFR0laMKZKJve8FawuW3Ds1YyrnhDh3XPgwEFdwxhQ9YBDaMu8CwwGmhtj1gN349+QtNa+YIxpha9SZBPAa4y5DuhmrS2oqV1n4HLgwIFEGEPnrbVjghzfDLSr6Zw/gzNwOXDgQKIB5PzUaffW047r3Tft5s/z95oveEfRFOkjrj3jixJl9MmPHCP4Wbwn+N3uewXfSbLgvZWYoXPaRvLZ7uffq8iQY1QepYbO/Vutaj7q9tYNknmV5/C24GnDlEmUipNC6YclH0v+Y6LUZYZuld5jP7U4Svanh+zP0Mfk+eqtqtYf85rkfQfJfMCHs68TXNe0HJgvNS0d0zftwb6CD1mlzg/I85iHFAwzr1B1FHVep64HqiIYda0CJX/CAMWV3X/ByTK0M36QjOmbmSL/tkM9WrAMA5yBy4EDBw0KBnB2FR04cNCg4CwVHdQ5XOkQ3RXadIOEgyGmFcSkwzPpkNwM3FFgDLHxvgpWlgKwufQGysmjmE0UsRaiCsGzAmre3HHQGGGoVnIt0lCnA1dztnFpQIG8jSnp4njCozLRKtYjDbS1D7vWuLQOsohugl9U+Ibg8xKlR1QpMYLr3EWtQwXqVNrbS/tZnbbuK8krJH+443WC5yofdO0v9SJXEk8ibehAWzrQrWMKyXQjqlKnawOWIjxsxcs2Yko2QNFc8JYDloK0eMDgogkuUknPS4SY7hB/AhgXpN3ja6diFZTPYiBLKWUuJfyGpbRaXBeqzGK5fOuYmCIT+nr3myN4q+lS6ElFuoAvUQZfQ75UmpUKX4pVxlRvdJT6aLeOVTF5p30g/xeTzpQxXy0zZc3KdUPk/3pgqdQzEx9ScVpTJC38RnoblMbK+86jPpa/pciguH6FMg+1LC7MJjTOjMtBuJFMBi05ir4cTxs6AODFQyHL2Mp3FLGKYtbRkol4WY/PA7L6wLC1k9yBzpzqF6RdsZDYBXp1hahMiO4L0YeTzl/8r1VEEdOI5mfymEq542bU+OAMXA7CgSYkk8l5tGMYif5Z5kbWMJUvWctKNrOOw5ksrkmvthVWS3hLYecCKFkgfr2yzaHE0ZcEBpPIEDL8U6ydzOYPPgKzAmwYrAkcRAYccd7BvsAAmXSiP33pxsG4cfEHOazgXbYwg5y9j9kLCR62U8j3FPI9fwAlHEZTjqU5o+nEQ9DjZtjxNmwbByULg7bnIILhzLgkXHhFbcMOX0oPKZ2buPlxmS83ZKbUNdz9pCeV1pkuRAYLDU2U+X4/rpLCzHcdZSzT2GXSV+nlTOnb1JIq7UPHeGmdQucGLjtJDjzH+n3HDVE05yTachZxZFJOLrm8wiImUxzwGkUqMKq9Mj5v9YEKDlJxTpnfqlglaQdGifJV13UhJ9ESWAgspA2ZnOw+iPjmF2Ga/40SviXu9UdgVZWXetK9O8X1PzBYcM8A+X7p+gDDlymh6EhJVaojhz8tDTcTrpUe/wfnr9z9vOQkcYheUfI++0/F1YLrGL3EL6WmteZeqd12mC7v88QZ8vzE6TIG7/nbLhD8XI8MaIxTuZFx06vivNzK/36f4AxcDmoLF3G05C+04lxiaUkxi1nDjeTzFZZyiqtlUUcONrKMHTxMPreRyCUkcilcOAk2/gqT74ElnwVtw0EEoQHsKoZSnsxBmJDESHrxCR24nhJWs4i/spRR5DEBS3nwBiIEXnawk8fZTG+YcBnEJMI5E+CC70mgS313z0FtEU7T+f0EZ+CqR8RwCG35gNb8mwryWcBFLOIK8vk5+MURjRKY/QqM6wFfXA2tetOHd+jC7UTTtL475yAYGsDAZaytu4ph2d2MzQlcrsswLQqyZTxKkxlqtqG8uZ8bdYng53veFPwa99OCj0IuWR5HFhX5lNGCp8/cJfjn/YYI/j5n735+I/8UxyZyvOATAkyUYnFxB8kMoB/FlPA1E7mdg/CK82Vfum6X/ldmJRJJin+j+BzFr1Jce7grTezTa6UeqL3Nrv/geXnB9IDn8alsePg5WnMWHnaxkidIV7URZqvczQylcWlkjlcanUoPnPWCjNHrnS81r+jAe0/57aPyOLXn/IKOsnbBVGRNzb++K+MFdRxXNcNzVbaxWn2A3opfpHhA7mP2G5Cz2RpCQHaasTnHBz8PwHzAbMcB9QBAT5pyMz1oRyI/k8M3TKSYEryiqmcjQ3Eeq3iSzXxCJnfSlfsp42iKuAavrqbqoP7hiPMOKhGLi8s4mNPowGaKeZE3WBlkVtHYUMwqfudSWnMmnfgbKcymiP+jlNfru2sOAuGI8w4AUknn3/TnDDL4mDVczPQDbtCqgpdNvE8+h1HBLBJ5nkRexlWtvpmDeoOjcUn0znbb73ISdvP0F6WG9NwVUrPS+XlRyLitFirdROcSFpEgeDIylkj7b2kPrPN5S/BU8gSPDRDpdN7k3fi8vgbSnDvphqWC5/mBefi0mYVKyHiBvwqu45y6qkKFXZB1Crt/q0QvpZM81l7GIgXmjAIk58tYorwUKZp5VCh1cql8LxOPU/l5L0j6cI/rBF9BZwyGEWQzkiPIZy1TeZyd/oCygarwY4fXZSzUHxfK/s3kCMFH3KyCApUP5+beVYFq1WLe1Kp9Sj/Z9qDFMlfw+SwZd9VNiVaDnpDn36PqNV9pZdBcq3WyP2vaq7iwxSr+MeDWyL4ecpaHqHG1MDbn9Nqda16oH43LmXHtJ7iAy+jIoxzKOoq4k092D1oOfLBYPmMWTzGBeJpyIo/QlsPqu1sOIOJnXM7AtR/gIob76cFFdOQzNvJXfiWXXcEvPECxkHV8yU3sZBODuIlMbTXhoG5RaSRYm0c9wRHnw4wokujNwzSlBU+zjPf3Ndn5AEMR2/iWuxnI9fTjcvLZRR6P1He3Dkw4u4rqxcq9pG8MmHko7+5LC/9T4/V5iVIL0LmJ2ldJ56DFqfp4Ux6V2sUupXm9wqWC63zEBwrv3P38tcSLSCaRyziHZJrzIt8wj2W7HcG0f5XWz7RGdhSyFp/OfdR625PDpEb2u/JR1/5UWt97K+V8wd1KT7x6q/rfKDebX6dJUa3v7VLnOfbBHwTX+mXle/81zzCACzmE60nZlQo7Lgc81epGan8uXfNS2ZlVqyvp/qbq7yscJRcenig5lRh0gtSotGf8X1vKuK25VyiD/TmS3jNR9e00pbHdJem89vJ/+WWWvO9Pyvpy9/Oy+6SP2T6hAewqRvi42nCQShOuYCwpJPEf3mMWYbiBDkBYvEzjVQ7JnwIp94ArGXLHgBpIHexHODOuAwSmDVdyHonE8yJvs4YNoB1CHewdCu4FuxNSnwBbjq/UtjfYVQ7CAWfgOhDQgpi4iSSRwEu8zdpqdbkc7DN2PglEQeqj0KMc5l+MM3jVAZyBS6EQWWNOxRrF6ZwtlT/X6l2pBbS6WWkD0tacuO2qPRnKJPy0oHrc11RVl/EYFVv018S3uYfhtCSZ25jBAppRKa5ojUprUj+q/DYdJ7WVFoKfVzhe8JMTZd7lLTxa4+uN4xqaEk0UBheGf/qNCHdRSjHlPMb/ifNzaS54tfFYaUh9z1P/vBskLVWiyZgXPxW8mi60W7Z5DPrGwOH3U5FeTEWZT8s7euav4vRX+8lArW9ukj7uOn8w/ZMqrfXJU6U+OFi5ybb8Rt4nba9UN5byLuvVdpng5Up+i1b7NYs/7iB41gkyL/Xkd2VM2sa0NNmfrVX9iVHe+/sMxwG1cSIaF//HENqSwiNMZAGRYVscQ1OakEUzjqAZHUikGQk05UJScREYl1hVgKEMDy6G4GEr5SyljHmUso4iFuOlsO7/CI1fH6CiVwJRMbeC3UxF+b3Br3Gw72hMMy5jjBvffGmDtXaEMSYNeB/IAFYDf7HWHhCKtAGu5ii60YqnmcI8NlFfmlYUhl40I4vraU5/4mkF+ETufDaxkz/IZQ2TyWQb5ZT5l1ot2IQBkoklhXiGsBA3rYjneJIYQ+V3eiG/k89kiH8NiufUx58IQEX5bWBaEhVzD17vclCVvB2EEY1sV/FafAYcTfz8FmCitfYRY8wtfn7zni4GKE2NYuWoqiVIp0mb5QlPqQtOkLRgjrK9uUvZ3uil5nmSbsiSU+yZ9BP8ZrXc0uEWlSlBp3A4A8jgjuIinijrCfTkrpT7xbnvcI7gXVkiuF46Vai5+RblRXxHomz/aGI5ni70pS2JxFBOKev5nS1M5A9WsYMVVASkJK2kq7h+qkqPOoQ7dj9305JOG1pA3OEkJpxIYtzV0PUaytnITj4hn3foeJ1cCuvwiK97y3X78KuVt4v635TpT8p/JY1LBFxXwpAMYtL/w/J+XgoDdAddDq7tdLmcK58s24sK2KS8frG05CnsWHN4BCqsRuO7a6UF+NAvleeOSilaauTSMEulJ5kPJM+7IlXwtoUBf2s4Nl8by4zLGNMOOBl4ELje/+vRsDuh7g1gMkEGrsaAw+nEKWQzlUU8UdYy+AVhRCwxHEZPjiKbVqSzk1JmsJZfWEdL3sIT4Jaq47D2Bh62QNFsKPoKtt8H7hZs6TSWRIbTlCtpyl/him/h55dg/gTw1lGogrccpp4OJ/xMRso4lnIq5c5mSPjRWAYufHOhm0BELba01m4CsNZuMsa0+LMLjTGXA5cDtDmoYWcYtSKVizmW5WzmDaYCZ9bJ68YQyzCO4Rj6EUcs69jIOH5iOmso8w9QJ+xPi2fPVgr4LwX8lyha0YQxNGv9F7j4Y9i6BL68ExZ/CHWRsF+WB5NHYkb/SgbPspwxWBVc6yBEVKb8RDCCjiTGmBHAVmvt7H15AWvtS9babGttdlp6wx243MRwFcOowMO/+RZPHWzLu3CTzQCu4laGcQxLWMHTvMrTvMoPrNw9aNUlKtjMdv4F92XAq6eDpxwu/C88mwPZw+umEzuXs5abSKAnbbi9bl7zQEIDsLWpzUsPAEYZY04C4oAmxpjxwBZjTGv/bKs11VSO6thFMj8GhBi06SeFjLiT5fkfXDFCcB2ucNSDUmfJnCLdF6ZkyZSe17hQ8ObKC1qn5VxR+NLu51Exr+CObsoCruE0vyd8VErVwPEiV4hrz+J9wc9H2kofWThD8OREaRNzJgUcyzmk0oJ1LOEaCljETnx7IRnVbKav50nBdQqM1oC0BdDFvCp4TKb01dbhFpkXrwc+htz/QcUYSjMfJvbBryhgKuu5mxO2TxbnK2fmalbLLfrJkINq3oLKtrt14QtURHekecxNNCn5H16PvJemDFBWNONl2s6nY6sSuYvS4sWxNqpWmy55j07ZUaXShq5Sf5y2bpYu0oyarI7rokjScYjuq6SF0ZSOVX/rzljV+L7AQKTbowWdAllrb7XWtrPWZgBnA5OstWOBCUClEdEFwKd7aKLBw+X+C1HRl7COV8nbz4Us4onmUgZwKtfhxcunPMvHPOkftCIRXsh/myUMZwMPkEhfuvI5xP4NRPhF+FFRfjtez09Ex74MphFbX9c1GoA7RChrt0eAocaYZcBQP2+EaEF07Di8npms5aXgp4eAbrTmEU5lMJnk8DXvcD+r9ddzhMJSzjbeZAkjKGQ2JnEcJE8CV8Z+fNUKykvPAYxv8HIQHjSApeJeDVzW2snW2hH+57nW2uOstZn+nzpOvVEgOvZ5IIny0gvZX4m+BjiTvtzBiVTg4V6+YDqf4CFcYdB1h3I2sopLsbsuBncfaJIDUbUsGbMPsHYNFWU343YPwx118X57nQMOET5w1al1c8/sGPtxTlUIQeYnUpNadqosS9++UB6flyhLTum0lv/jMcHfLh1bY382xkqvFK37ZNGX9jzNZh4ll5e4RU0qA19fl+vSGtf3qlxZZYxYIm4eoDPHksYH7OAeNlGCZZyqH6bjvuaomlVa/3ts9t2CDznsc8F/yj9K8BNStAbWUXBt9axtcW5d9ZTgzACSOsPAT6BJNzylt+AtqypJFl0gT5/bXqb89LpRps0ohyHmZgWeb+i16wWI7wvLukPFBrb3kCKNtgUKLD03t5967fnqtVXExa/DZK6aTt+6dqaamevvH/2BT1P8bsV7K659FgPSq7JnQU5BiNbNBxub81ztzjUnONbNEYVoUmjN3RTxG7nqQxsutCGWN+jOQJpyL5u4hY2UUHdfJPsdu1bAd/1h/ce44/+JO348+yck28LGS8FEQ5vng5/uoGaEcalojHnVGLPVGPOnmofx4RljzHJjzO/GmL616aIzcO0BXbgCN8ls5Fb2hyNBVxJ4k+60JIa/sYg3aZQrbagohOl/wVNyG66Yc3EnfkW1arPhQNlK2HInNBkJyaOCn+9gz6hM+anNIzheB2qKkzkRyPQ/Lgdq9c3jDFx/ghS60Y6R5PI6pSwLfsFeojOpvEQ3yrGcz3x+oSD4RQ0c3tKHqSg6F+M+xjd4meTgF+0tcp+BkvnQ+mlQS2cHe4EwzristVOhxm/l0cCb1ocZQKo/vKrmLtalxtU1O8k+n3Pobj7kExVaoOJd1jwlyzLpfDZdbkzHYel8vwSKBM/eLmevW9JSABfN+AEX6YxgAcUBgnwZ0iolJiBiO7BUmT4G8PmKMwA4LBa+aw2l7mJu4Gc2+10l9N82M1fmUd7bTAofOi7rsA3S5oUS9XVYomQPNXHPPGuu4DqW6W3OFbztY+pevFDSapuhlTJS/BnQ7B3yzRLmcj0V/iIiOpfzNyXsjJikyo2pamy7bXaaDYRjplJR9hAV5VXBqd8kSmEo8P81/EV14zWRtFrZAGVjo1JemdJbxZBdJOPAtr8m9be0V1Sgls6FVAPEypWtBO/UpyrnN3sJ5BSFqHF1MzbnreDnAZhs1gDbAn71krVWiHzGmAzgc2utFKl9xz4HHrHWTvPzicDN1tocfW4gnBmXQjznE00vdnK7GLTCgV4x8H1ryPPC9QGD1gGF4g8h90yS6UpvnsYd7plR7o+w9g3c0TdiTJfwtn2gYO9mXNsqM2P8j72NGfqzQTbobMoZuARiSeImyviZEj4Ja8tdouHbNlDghcEbYcuBOGhVovhT5nEbSWTSk4cw4Y5knH8zUEFUjN6ec1Br1F0A6npk2Zx2VLetrAZn4ApAAhfhpg27eCis7SYSzeetfG/20E2wtuGFZ4UduUxnMQ+TxhEczI3hbbx0C57yZ3G5z8GYbsHPdyBRtwGoE4Dz/buLRwL5leYNNaFOQ8iS/yhkyItVulbhhXLcTBwgd+86TFKlxmVoERs6ygCYrBzpa/RTtoxV2qb8hmemVYkTUcRwATezhd+ZhBsYXC02a/0G6c91Q9uquKSfkK81c0s/f7vwRWocnWO9PMf7DDtoPcOAFchlzC+FUihp30wKK7rE/N25ygV0l0ouk2FbcLTiahXVNbCOO8H1wnk3yZJZw9dJnWjZEBmTp2Odri/15Vbe4d7FjVGjmcGzzKOqvNxlH0iram3zjVZLlKYWte4xyPgrsd77YNMZjNY+4IFvl3IAX3aF7PvdyPf6YW4RvMM58j5NfUf5aaq4qyIjNa00HZf1qOLjJO30kPKxeyHg+YWEjjAaCRpj3sVnf9XcGLMeX5RaNIC19gV8Jt0nAcuBIuCi2rQb4a47dYeeDCaepkxTQayh4p/JMRwf62Y8X7KC9cEvOMDwoGcXXUwUo91nk8dG1jE3+EW1gXc75D0Bze6F2L5Q+mvwaxz4EEY/LmvtmCDHLaho61rAWSoCbqLpzQlsZg7bqtmo7jtOi3VzTUI0TxWWM7OB5BzWNSxwZUUeuaxlMFeSpIt0hIK8f4FnB6TdFr42DwQ0tlzFxopuHE0iKczng+An1xId3YZXmsQyo8zDLbvq1+jOZaBTLAxLgVFN4dSmcHqq73FkIjSv53l3CfA9z+LCxRD+hgnXbendCXnPQeKpEJ0VnjYPBDSAgatuX7oQXw6bH4lZKiJd90bF0xQcJD3nf+BYwU/IlrFNGrpk2FK64sLNYQxnM4v92lZVnNmWXKnrjG0rg1sCcxUH8uPu5y4MVzbPJIZ0Jse8wLUt83nsC7nDlTZ8g+C75siYtRsHPC649scv/1C9OalVT5vFwMltYXhX6NEKMptDnHzrqqHAns5Kr4dp3nK+rShje/xaygN2pf/N38T5Oo7s5fYyL1T39zzke3derPQna8HPLONRunMvQzgaCuXxallX+l45TPHKuLGSp8FeT37G/WwK2ATIui5AD1VeX5N3O5L7EKNi9IpVCMfD71wn+K2PPSUbnCbpO1bWybvpCZUYqDa0Z02Tgt7hE+Ts3QamWoZJm7IR7oB6wGtcnTmaJJoxnfDZohzDADrQnnf4iB1a+d1PSI6C8zPg9HZwTDq4XbAhH37dAF8tgcUuWFYAhRXgseBtDcZA+1joHA+HdCqhqyuKS6LiuTo6gUKS+YldfEAuU+rIC2wL39GM/mRwETR7CXJrjEGsHby5UPgSTZL/zh88QQVBN6wOeFgXlEW4keABP3AdwlC2s5YN/E61rbZ9QBqpHM8xzGMhc+tA12oaC9f0gGsPhaYxMD8fHloEn34PswP3AvRKyT9hm+c3SUhq56ufGA8c447h1LhyjieFoaSwlGKm8QczWbrfLauX8ARNOYzY/i/AF0eADcPr7XwKkq8llbPYVq2UlAMNa6DCXdvlev1UFj+gNa5mdKQ5nVjCpOAn1xKjOQkvXibwddja/DMkR8Mj/WDNOXBPNvywFbK/hZ5fw13z1aC1FygGvvGUcS8bOJaF3MgavMDFHM/DnM9R1UbA8MJDIct4GpofBl3/FvyCWjW6lkJ+JIUzOcBv+VrBGoMnKqpWj/pCneYqZrcyNueCKr790Zrno8n5Mt4l+kV5/OubZO2+HSr2qCOrBT9yjtxq/7H3XPpxCLfyEsWUVsuX24j068pVcWATqHIh6MZCjqEl99GH51jE3+eqwKM7JGW54hmK63vi7Kqnw9vAS0dB2xh4bxs8tB4WPKPO1+OLCnrNvF2+F8te7FVjfx484TvOoBcH04Lf2ch67qKQKp/49zlLnL9N7Q6+wfmCd5gvY5829JAxea3L38FE9adiZxbYTUSrsC7mSLryZZm/l6AyE1rNOw56fgRzT4Lcr0TU94JBnagJ3afIxMjy3vJ4tCx5CUMU17mHapPz81HyAl0LoWfpPMFLY2XOrHjp7FLm5HhDylXsk+2yk3Jqt1ZMM8WOH1edwpXI4WSRwxKKlfi6L4jDzd85hGUU8Alrw9DB6kiMglf6wVfHQoEH+v8O5y6FBUXBrw0Vc9nIXXzFK8ygC80Zxji6MJL95SvvKbkKiMEd98/wNLjtMyjbAm0uC097jRgWgwd3rR71hQN34EoZTRwx/MyCsDR3Bh1IJ46nWYhnP5gBHtIEfjsRLuoMDy+AvnPgl11BLwsrLPAdS7iBT9nK7/ThcvpzMy72PAPYZ3hX4C19AlfMueDqHXp7thw2vQHNRkB0evDzD2BYDBW4a/WoLxy4A1fT88ilgJVsCH5uEKTi4mw6Mo0tzFd20uFA/yyYNgySomHw93DbHCirR6PU7RQxnfuYwyu0pT+DuJ+k/TB4eUsfw3pzcceFKXd00+vgioZWNVt6H+iwGMqIrdWjvlC3Gtchxua8XsULe8txMy82VfAW+dLzaWlKzVpEe4/M79vhlu3NoQ8AcTRjGG9RXvEQ5Z57dh9fHitzE29VHvO6VuETC3zC1SMt4f+aW27gc9b5wx8+el99OGQZRXhK1RHkY0mb/xWAEcPg/ZdhfT6c8E9YXSkNSUmnmic7IxS/Th+Xla+TmucJvut1NStR57PLFxh2ehMY3xa2ugq5mRw2+bWlNirBX8d9DWSq4PpD0Mbjuz7G/IM416Ms5QJ2MWv38QxWifMD63UCjH5X5SZWSo6ZM8CVQGENs7jECWqnrKc6QXnEL2shcxszF6udER2DpsdNpSyUHCf5i4mXC35tjnKOCchlzP4McraF5sfVMzvGTsipXQZDJ7PJ0bjqCu05DoObCs+bwU8OgpZRcE0a/Miq3YNWuHDuGfDJGzB/MQy4P2DQiiB8VADHr4FUohnHkXQkKaztl9nn8dr1tKk28u4jtr8G8T1xmT7haa+RwtG4IhBtOIbtLMRWs9Hce1yVBrEGPmBe8JP3AqNPhDfHwZSfYMhpsC1S68EC04vgambiwfIIh5EW1iVECaX2MZLoQ2K1cjf7gLz/gq3A7To99LYaKRyNKwIRTwtSyWQT00NuK87AlU1hwk7YHMbo8kMPbcb452HWbzByLBQWhq3p/Ya1FHIrs2lCNA/Sh5gw3tTl9i0qyKeFCqnYJ3h2wK4pRLlGh95WI4VvVzGqVo/6Qt2+cmVpbz+2xLYQh2cjl8pnvidNpTqPlTOkBxNvFfy+/IcF35Imv/nLiKGDXwtZwyxeiX1AHF+H9Nv6B/8SXHtuPdYth3SymdnkB5bQVRxzDZajTfuzpP62ZpsKtOrh07RaJMCE8yFvF5xyCxRXho7pIPwHlOb0ikpG1H5cKo6s11WzBZ878kh5gvb2u0W1f4akR3X35YF+wGrGMobHaM987qTShVfHeem80azF0kvtgywp0rXhN/ozjKmcRz55nMGH4vjom5WmdZLq/+0Bz0f+D9eVz5I4rQsULmXKKOURHys94tG7t8oiPq9FquAPZMm/1f24tAC/dd1TsgH1xRSnMp06D1JBf3Mk3fxalQl+eXboW80+cX4/7BSHEQfcjKst/cljFYVhyFkbTSYr2MF8UStg3+E28OEpkJ4Ao/8Bm8PTbJ1iMUv4mm9pwRA6Er7K0r/wM2A5nP6hNzbjU9/P1s6s689gwVkqRhKiiKcZWWwK2J3aV7SjAweRwufVQuD3HXcOgIHt4dKv4Nfw2YLVOabxE5v4ggwuIpVa1fcMigLyWMxCepONK9QPzB/roOB3SB8alr41PkT+UvGAGria0x0XUWwNg8tmH7IppoJp1WpX7Rv6tYE7+sPr8+DdBjxoVWIJT1DMeg7hVlxhEut/YxaJJHFwOPIl/5gEaUeDK7KXRPWBhhA5X69xXLrWavlIyaO1fq51Hh09opb3b4z6i+CDOJ32nMJERuCljOXKDULnIurjlT7ssbh5k5G8VxDFJZUhCutl6EzKkdIXPH+GCrwKyB2MNjD7YEiJhx53wM4SqhtVtFM8T/HeiqswnCO6S0/4X56QeZ7VPOn1e63Gik4DZMbBygXdBW/VfSX9ieNj2vA8eeQp33btoT8O6VHVdqaM4Rvf73QMLk7heXawmnzk/1bnE+pi2SUqdTTuk1Ew9lN4+Rho+aM4tmyUisu6WsVlDZAUFV74cb8TBT9t2VfyeKY6PkEer/beXyjptDZyFtuSrVVtZW9hXk5ZSHFcWdmJ9qWcaiUQ/xSDzC9OHNf+Rhp9yWMeXkJzJD2KdiQQzWth2ki8/iDo2Q6uess/aDUS/EwJb1HA5aTQstrIu/eweFnJZNrQmyhaBL+gJqyaAl4PdNIZ0Q4awozrgBm4ooknmU5s5/eQ2zqadmyhkGlhGGTSouDWDJjwG3wepjoRkYT7ySUPL8dWC+XfN6xiCi7cJBOiPlWSD5t/hw56+uTAYigltlaP+kLQgcsYE2eM+cUYM9cYs8AYc6//92nGmO+MMcv8P5vu/+7uO5rRCYOL/BCLYcQTRW9aMiMMOY4AN2dAshtu+ygszUUcdmJ5ljw60pWDwmDUWMBG8tlAk1AHLoC1P0H7I8FEuE9xHaMhzLhqsy1QCgyx1u4yxkQD04wxXwGnAROttY8YY24BbgFurqmhbYlpvNrvhN38qH4ylqdUxY70OkiJYFoH1+l+Kh7mgin/rSLtb4ZOcC8nsNNvmPSqKuF2tzJO2qqWI99OH83pzSAmC56cd7A0f8yTr53/odK0MiRNOXIzrXHx95h0/ustYUHHeFk3Uvt3Xam4DtfRmTbKf+uXHKVpaWiX5MGS9u8uzRZ/3iBj2nR/BgTEaS3BxR+cSl/G8Kbfm1/XbWw7R2paC/pJ4agPv+1+XsZnNLeXkFG2g8qiiCVKZZmXqHzaJynh6FIgcSrEXgU/9IblVXFtmeOUpqVyB/lS0nJlGX/aeVKzeu6tSwTXMWjV7mulz01pI/VAPWAEchsGm6HKgSuSEXTg8tc9q7wto/0PC4ym6vZ+A5hMkIGrXtHkKNZSyE7Kg59bA05Jgz/KYXoB1QeLvcQNUUlEAQ9X7MJnmrz/0DYBhreFlGiIdkF0iS8e+JfNMHkdYXAk2zPK8fIWi7mePhxJK2awOfhFNWAn35NursTtOgmP9919b6jEL8p3P0YMXA6o1xit2qBWgRjGGDcwG99e1zhr7UxjTMvKUtnW2k3GmD9VS40xlwOXAzQ7KOHPTqkbNOnHwjBYzgxNhW92gCfomTWjKYZzXfG84y1mTcit/TlaxcAZLeCslnB06p7P21UG362Hz1fDu8uguGLP5+4rvmMd53Awp9Ep5IGrmDlYm4vbdVxoA5dnE5SvhszDQ+pPY0Nlyk8ko1a9s9Z6gN7GmFTgE2NM7fZKfde+BLwEkJHdrH5cpGLaQExLlihbmr1FVjy0jIEfwmACca47nnhjeMkTfvvSaBfc0Qlu7eB7Pm8X3PEbfLQGNhRBuRfKt0K0Gwa3gxGdYWRnOLUT3JENf52CMqEJHV4sX7CGS+hGe5JC/AqxeLxTcbuCLH9rg7JfoVN4gmQbCxrFUjEQ1to8Y8xkYDiwxRjT2j/bag0BwSR7gBcXRQE16bSne6q6nZ/rIbWBq6P+IxtUWsDm9imCr/YLS00ZRFfg8rLLudT+vPv4b7G9xfkJyEFE+28d3zcL6EpB5s/0zyxmC1V1F+/tJesmfq+EkTde/Kvg53VowvR8+OlGf8CV3qHUcVU6QF8Xw3jB96NPL3j9BTi0I7z1Czz8HSzaTHX/rs2+WePXy+HryXD1chh8KPz7Kvh6JLy3HK77Cbb4rdtjussQkrgk+V6d2u8dwfX/NpmdTGY5F5DFCNqzTm2SfNdb/sFDx8tihAvGSs3rV3c+R9CJr2KvoJBcsmOlSHf4OqlpLR7SQfBKPbUFa2mdcAo8ngRevyKiJ3HaT+suSaPVOnvSWzItqR8zBW+5XX3zqbqOBbfLvFB9Xx6eo/S6J6uexkmbsn2Cb1cxsgNza7OrmO6faWGMiQeOBxYDE4DK0hcXAJ/upz6GjEQOweLBa0OLN+hDKlsoYYMqxLC3OKY1dG0HL34V/Ny9wS03wMwfIL05jHwRzn/LP2jVEpN/h95Xw11vwakdYfFZMDT08KvdyKOUn9nAcWQQHWIkzmb/l0pLuoXUTjELwbggtlfwkw8QNAR3iNrcPa2BH4wxvwOzgO+stZ8DjwBDjTHLgKF+HpGIpzMlrAdCW5Z1pwnzwmAW+JfOUFgCH04Lfm5tcdM/4OF74KNPofsR8Pk+lnQsK4f734FDP4A1u+DDodA9jIEu37OKZGLoS2i+7ztYSznFpIcYYlFSWdkpJrQBsLGhwYdDWGt/B6rZRVprc6m+URyRiCeDElaHNPlNII7WxPNxGOK3RneAb2ZDcZi28k7s5xu03vsQxlwU9PRaYWk+nPwV/HIqfH4i/IMYdoSYcQDwO1sppJz+tGJmtXiWvYEll1Wk6XybvUQ5m8FbDDEHh9ROY0Kj07hCRTIFHMsPu/k3nCCOj+Qzwa+eIDWtlaOkUNOyVMpqrebL2VCrRXMBA6e1J2H5/5jbVXrKv6KM2v+PxwV/M8C4rgNtAVjPVpL9poGB2sVrKqFsa4D+5euM78dhqdAuCW57C+lDr+yw+FBH+KsB80hfvtvBB8E7t8PcVXDJeKq0sdfV5Q8onqH48YqP973iyHnw431wq2cgJ5fk7ZbiuqSsEKdrffLd2dLS5sTDfJ76HiCHzRxNcyYxH+v368pQNTB1PmD3rTIXMb5FEXFMpQkXMJCfWKc0Nb6R52eVSr+vzVdV6aE7XRspaXYkc5r53ryhw9RUWH+GVRzXmrFy9pisTCUP/0BNf3WwvrpVNDprgVN+TNj+TlUNxIrs0L8NG8LA1fhTfhIOAncc7FwS/Nwa0NY/8iwPcak4uhVUeOGLb4OfGwwJcfDpY77l3Sn3QdF+CMb6dSWMfQaOcEfzbGxyWNr8mU0kk8RB/i+DfUUpc3ERT4wycdxbFLGOBA4KqY3GhEaR8tPgkegPR98Vmr98S1qwjWIKQlwuDWkOv+TB9h0hNQPA7RdCVgc4605YG3RPd9/xyS/wQFkhZ0fFcYI79N2mX/1LxC7Vpn17hzJ/TcyYEG1uitlIXLBpzwGEhpDy0/gHrnj/tlhxaL5Z6TRjfbU8m71DjMu3VJyeG/TUoOjUFm48B974Eib/Gnp7wfBEeRErvB5ujQ49iHgX5WxmKxkhzrjKWY3FQzSdg59cA0r4AzexRJMS/OQDBJE+cNWpxhW/sYzud1fNfL65t4aT/wSdrpR7+2tekNpCYqGs31VwWjQxJoM4oGDYVrqUyoHnydjrBb8JWe69Z0Dlnpak8VVRDN/vrBKD2resGgxfQMZpLUGKvSOjJtInFeLc8HMBvszOQMxR/JZDJX9c8htvB2vg1lXAKcDX+nrFdVyYDpPQm3OnSJo3OJ1HToGXb3fT76Z0lozLE8dXq9nTfw+T5mq3qE3n7SyiE4fRmRWA5VGVLdavo4x9ilEz3UpNrQVbKKIvq1XhyjaXyrqOmTNl4FurnKolf6s3p8Iz1zL41PWweC4rF0kttdO76s1SdRE7jJP3XdMr8wRfcKbcQOj+rpz9b7hCFmpsu0rmbVIq09Q+vXeY4L0D8ji9Okl1H1BZ5SeS0ehnXC7a4rW5EELsVTxJxJHIEo83+Mk1oL8/rODnEJeJ6elwUU94Yz5sCr02Qq3x5hewYSvcdmHobW1mGXEk0ZTWIbWzi3Ukher1tcU/qLUMbQbYWNBY4rgaNIxpiQ1p2x2a+J1RV3lCy1jqlQwbSmBziCL6aadCXBQ8W8d5wWXl8O8PYcjh0CrEm/YP/y5iM1VZaW9RxBbiQzUV3OafUTXX6QUHLsK5VDTGDDfGLDHGLPc7yejjTY0xnxhjfvdbaAVNKWz8AxdNsTqnYi+R4Nc+NoU4cGUmwtIw1EgceTKs2AHz66EK0Mf+aJZhOt9qL1HgzxBLCVEUL2U7MSTjJjr4yXtCnn9p1iSiLeXqDJXlyWrzCAa/QcM44ESgGzDGGKOjfW8D5lhrDwXOB54O1m6dzvV2tEnhg3sH7uY6F9Ct1+cq/GXGCzIto5tHJU1PlDQ+qxx3UlOsdxXxReUUJ8mbe6naRn8YWafxM3w6TTP/tu/oZs8zOCBGZydV4QGZG6WGclqbj2VnsiyZKfC/QiALDjn5N3F40a0q0Xe1pJV+WQkuOK47vPAj1XWqQKRKGp1RIPjAw6YKnqCW0p/HnSIb8Gtki4Gl5VEcU9Gch7ZWfdCf7CD1Qh2jl4pcH1++ZhzD21oWl5zC33JHck8HaUB2kgqWmspAwbP9b0gCB+MzZW9FUUCZuMxJKplThXl9nVWVoP3M7GuYgJePbrqMN246ii83qirXMtyPZb2VJ/0U+Vpb3HIw1nmb3ROlxtV2vtS0/ugh/ZLuVn79scqEKHAAKUT65+8LwqxxHQEst9auBDDGvIfPEivww9sNeBjAWrvYGJPhd5/Z41Kp0c+4MGlgQ51xNcXiZad2KtwLNHFBC7dhWWh2YAxNgjgXTNjHlJ5w4NMiL8fGGZqE6Fm3rMKSGcJECcDjlwESQtwR3EkZyaHM2hoRwqxxtUVaJa73/y4Qc/EZk2KMOQLoQPXyMAIHwMCVgrUFwc+rAbEkUUYRXvZdnG/n/wJbG+LAdUQClFuYFgYXgH3F9yWWGGPoExPayLWmAtq7Q2vDi2+2EkdowbFFVJDgDFy7sRcaV3NjTE7A43LV1J/9g7Xm8gjQ1BgzB/g78BvVPHwlItstLCyIAxuam0MUcZSH6AiR5h+4ckPbmCQzFlaVQfn+8R6sFRaX++67zGiYEsJGww6vpakrtIHL+mfB0SFGcZfhCdmxorFgL1N+tgUpT7YexA5MO0DEqljfzOIiAGOMAVb5H3tEnQ5cTXflc+b0z3fzrwdIIzitixxzm1yvH7lV2dJ8oF5Ahu7wU8qRDCKGdXHNWRN3BD1VhZ8fOFbw71XCXmX+XCwxVFBSzRfpsblVHlwLe0m9MUbpEN3brgY6ktxmGZkUs+hFpWkNV3+Ldo7wx2llXgjLcqketzVY0rQrZW7jGW7pc/6j0oy26UKM76nZR1U6HOuA0mu9DGy2mV+brQbgPmVSpfXLhfny/bmuw+NkcRSJDOKGDv+qFqf1NucIPlsV0Ty10O+iZNZDAkSpgeu4ITKh73glgAZ+MD8uPJ2YuJkcZLfxcenpoCboK3vL3cZ1aic080upca0elLHH1wKqfySV/pZ+uIxxeWSG3Ij70H2G4IOZvPt5cogOKBB2jWsWkGmM6Ygv/fVskP9cv21WkbW2DJ+iONUGWSY16hmXyy9aekJ0VI8ijopqTn97hxT/jVAQok1zl6Y+j/j6hAXWUcxBIfrkF/nf04TAUXGvO+P7kIc644ISMCH0oxHBt6sYnjxEa22FMeZqfKa6buBVa+0CY8yV/uMvAIcAbxpjPPhE+0v22KAfjXrgMv4/z4YYTewmKuSI5Hj/MqQoBJ0s1g1JMbAlDCEVoWIbZTQNURMq8xcuiQmpHd+XUkjhEIClHFOPScORhHC7Q1hrv0R5avgHrMrnPwOZe9Nmox64IgmVamQoio7X30iIslBY4Mbgraax7h0q/wwbUjtVrYSO+imJEImIdFubuh24XBD4pdZC2dQPnzRFnq+ovVZyc5JqX/01R385G06CLvPX0mXlL5Qo28Nb5z8l+KR+0iv8S3wv0I8kooiqpoEFSAtMSZL6XIfO0kanxP+HF5LETjykX7FWHO+iPJd2DpK7ZPO/OXz3ItPVmur+XWdL2tJdc7bAokeVxqY1Nb1qUvJr4apUAH591R/gpWOdz5L03hTpyb+ajlT4Z0mlRLOInuJ4Gy1YKlyf+AQAsUTzL6CcaIoD6hl0Zak4/1/8Q/CvA0TFzxJHMIjmGKKYnDiCM7/5XJzbyS0D5jo9pALoHpS0n0fmWU51Sz2x2nurNa/XJG0yQ25FX1yhTPEDPidRmwgZDSFXsZHPuMLzDRrO7/KQZlz+nyFGEIQFbuOrFhQKjP/diJy5kjPjgkZUnqzBwvp1KVdo+ocHD/Eh6h8l/g9FXIiVhvPLoXkEFGBpngiL/wh+Xk2I8n+rV+Chus1o7eD2X+cNcdPD4MaGoD82JlSm/EQyGnfgiscfe+UKbferhFJiQ/xHbvd/sNJCnIIvKoRDQksTDBkxbshMgwUhmhcm+XclC0OIkYvzLzfLQtz1dRMb8u5zY0HlUrE2j/pC3c64tgPvVdG+BbK2XjUTB2k7hFmmjkuZiLlXyI2JXouX+QohpPu0jzhtl6zaG/Llz4IX3ev7YDWjD01oyVm8L45PqxhaReLkTb8lXzoW3JjyD+AN/slj7GBatfwznSuYk6tEJX+o0MKOcFI7quUiohKu1xXKWKPXSlQVjTnq+u8V1x71eVVPD24JUW5Y8CJV/0/V3fntZHXo9wfI96cNm7Cks5My1tBO5H1CdZ917Umf4ReG0vyZITtwsyPgTdG1DHUM3sKAsmZjZ34Ehz4KRas5c/nn1fa3cg+RvNmjkvOQ4ipFeMTiSfIXT6jzn1Jc2nPxWI+rBb9pwnPyhMBspzCNJZG+VGzcMy4AWwSu0Fw7yykiilDb8CUZRxOaA8HCHdAqAdJCNyLdZ/Twj8nzQ8yXbEIs+SHOcqL9SndpqLMlEw/e8FcVb4hwrJsjAd6d4GoSUhPlFBBFAq4QvoXK/dY6sSF6R83254sfE5pbcUgY3BF27YIlodUfIZ0E8kJc4sXhc1IoDjElC3eTqkrWBzicgSsS4MkFd7OQmijxDzqJ1dZntYeXEkrYTDwdgp9cA6ZthtwSOO3Q4OfuD7hdcFo3+OwzKA8xYbwtyaxXpbz2Fon+GWyBztPZK7ggqimUh6EYQCOBo3EFohXY26qoUffJ10Nk7uLw01Qgl9a4VCxTr0HqhInAyFyIae573lFdP0Rxlfs44gm/NtHBDWfcxFTOYU7AB2TIDVXxPmcgcwF1bl0eqexkPTF0Jo9UTuAbasI5zd6WPO5/gC9l/tMNcFoWRC+G8sqAfhUbtKtEVYrWy7r37lG/OFHSyf0k92tggwZDeiJ8UABcGXD8BRVKMF/unm4dID2qdpRmkhIby9SK5vzoGUjRLrn2nbJLBt2N7iDf30rNqjvt8OLlf/QX+4pNlR6pMXbmR1VkTRr0c8GmXN89Ju3r+blC3igjOknNyqqKdU3uUiP6PElVWieo71Wr/pd9lCD59ChpwHDN9peqiIoB2xd4cYUt5Wd/ofHPuEq2QVxoMy52+YIh00PcWdzFWpLCUL/vo7WQmgTHVasvvv9x5mlQWAhfLwitnYON7ztzqQ0tlSqZJhSyK7RgiFj//VHmzLgq4SwV6xvFWyEhRC/xnb7ty7ahJAMDO1lNFAkkaDuAvcT3m2FbPvxtZPBzw4m0NDj3bPhkAhSHuEw81PjCGBZ7Qxu4UmhKvqqivdeI8+uOJSEGpjUSOBpXJGDXOohJhpgQHDLLdsLODXQM0Q1hu3+9lqbSW/a6O1546mMY2R+yDw5+frjwj79DYiI88njobR3pimad9bA+xKDPdFqwjRADyhL9oSNF9Wy7ESGwRL7GZaytuzSHbtnxdnxOldBUpAaCo+eryqZKc1p4n2pPhbOg4m24GTjuDHjkAzi7JzwshZ7FHaVQnnXfGsGfu6vKXWM0N+AmgTcCgm4u5ZXdzycrQyztL+WLRDaMYjwbmcHZSkM4DFmyR/vva81swYghLH4DNuZCv7+CRxfq1nUU9YbZ54prtfNGSVtPgsXfwFdT4OzrqOb/pePIqsWFqTivDZ/BD+thbKXUd4qawuXJbIdWveQf2JxcknDzM335F+uYxXfiuI781h72SwLqDZzLQRzKRXzCX6iguFpp+bFzPhK8mtaqpFgdX8jtkm7op+ooPqTqKOrcRTWzXjBK1mkM7O/Y7FUszCkOKT0jPrub7ZjzdvATgUWm7+wgRoL7BY1/xrXF/y3aKjRtaQebaBZymXZLLotoTrfgpwZBQSFc8ywcdjBcc3rw80PF8/dCdBTc8a/Q2+rYGtokwbSa86iDIsO/dF8VYkhFPM0pYxcVoYZUNBI4S8VIwKbVvp9tO9V4WjBsZyOxxNEkxADSrcwjmXa0DlEvA/hwCkyYDo9cBsf0C37+vuKeITD6eLjtCVi+Jvj5wXCCP6h+8vqazwuGg/0z9hUhDjiJtKIIR9+qhMVQSkytHvWFoAOXMaa9MeYHY8wiY8wCY8y1/t+nGWO+M8Ys8/+MzKJ0uVtgZx5kZIXUzB/4PrGtQqyavMG/136stkreR5z/MKzcBB8/D51C37Cshov6wt3HwqsfwlOvh6fN04+BxdthcYgVvXuSSD4VrA0xaj6ZtuxkQ/ATDxA0hErWtXnlCuAGa+2vxphkYLYx5jt8xewmWmsf8VenvQWfqrRHJCwpoe/ggPxE5WOk6yIqWQKlBLBepnDRTvtznen/WbAIjjiEhztcJw5rHaToLhlLFOgJtYVNWCo4AktLfDUJcwMCcLTH+lWME/xYftj9vC8bOIUYNgYYej22RPpVdeoq4w10XULu8f3IB0Zsgplp8Pl7cMyTsG0X1TWtwYpfqbjy82IODOsDL42Cb3+DKx4FAsv9vaICPk9R2Qn3SBp9oe/8NGMY3DSJR18zMD7ghC5S03r2ZFnM8H1l8DVtylC6HAYzymD+vMOZXyFzIx86TvpvLVTL823+/10ULhJpxRwW7Na97pijkgmnS4oMN6RwlPz+T3xdbTioYt2ppXnyF1ohGiBpiTr+nnovAr3HSgiDIReRbyQYdMZlrd1krf3V/3wnsAhfXbTRwBv+094ATtlPfQwd2xdDs9BmXB7KKWIlSdV2APYei5hPB9rTJMSSWpVYUQKnvQQdm0POLXB4aMH5AJwzCD66FeavgTMehorQohZ2Y2RMFFHG8NGk4OfWhDgX9EiEWaEF3tOCFFy4Qt+ZbERodBqXMSYD6APMBFpaazeBb3CDP0/CM8ZcXllz7Y8QY3/2GbkLIakN8SGWjd/JApLpQajS4EJ/taE+hC9vZ+oyOPpxMAZ+vgkeOxvi90GCyOwM334Mb98Iv6+Gk+6FnWHUrMfGRrPc4+W3EPMcs5MhyoQ+cLXzWzFsq2ZNcuDCYvB43bV61Bdq/Qk0xiQBHwHXBSsdFAhr7UvW2mxrbXZ6fdXb3Owr1x6qPpVPDlE0IZHQgqdy2cYKVtOPviHaCkrMXgs974dXpsP/nQzrn4ZrT4C2tVAfD+kKD90F86bDEYfBVc/DwFtgk16fh4DubheDoqN4paQs+MlBMCTV58E/NT/oqTUig3QqqOAPZ+DaDes1lJbE1upRX6hVHJcxJhpf5M831ton/b9bAgy21m4yxrQGJltru9bUTnaWsTn/CfiFzp9rrbgOMFc+R6ibtkRZccdVtu9uAtn5sOl22FJlnjSrtzRK17UGdSzW1Wu/gIM2wfb/g/zHKQmYY85J7CXO7VkqE9SWx0o7hxXcyWmcx3v8h+UsYqP6Y7sipyRao3kz/zzBS66WsUEnvpTHTTEJDHb7pl0/lXj5uNiyyWMpt1C+MIooFwxuCSe3gQy/fvX2z3DDe7DlTSQeUfzKryQfoXIdL1TnPw7jboSLR0C70ZD7f/JwytHSx71ZrEy/0Xrk1zQhjihuxtcP/f5pz/rTVS7pmE/8dRkHfE9+i7b8HCD6DX9XBWZJay+2P1XzjrD2t79//sOCz+2hfOPuVoFhKqd284UyeLrVfHXjj6p6mr0BckptSN+Hrj59bNyU2q3li1PS6iWOK6g4768s+x9gUeWg5ccE4AJ8t/QFwKf7pYfhgKcAipdAwuHBz62xnc1QthDihkB+aOHji5lHAXn04xiWsyj4BXuJH73l/FiSTxfj5pSoWEa7Eni8acAE2y8wF1bA95vgoQ/hq3mwPowzrECkJMF5w+H9iZAb4iwpDsPBNOdLQlxvAqQeRr5/s8WBHxY8FZEtztdmV3EAcB4wzxgzx/+72/ANWP81xlyCL1b4zD+/PEJQ+AukDA1+XjAUfw/Jl4JJgBCqBnvx8gvTOJ4RtCODjYS+fPozLLceHi8v4p5tMbRyQ7KBaAPRc3wF5xfmQ6mX6tHfYcb1YyA5Ef71XvBzg+Eo4onGzbxQd9CSukJMKvksDr1TjQjWGirKG/jAZa2dxp6L0xy3h99HHnb+BM3Pg5gOUBZCFGXRJ5ByDcSfiE/y23fM5if6MZDjGMEvfBxSW7XBZg/sXpCFGEO1N2gVBzecBO99B3N1usw+4ESSKKSM+aHqUs19084dzA29U40KBq8nsq2b67R3xUmxzB1QFSUZO0DOMrLGBxlQzqv5cJxecQXuc8b4TZb+OBJW+V4no/dqdbr8IHRYJqOpv84chMHFYHawveVfaf1u1cB1pLr5V46RjhSpysHgG3x1GA0buIlu5JdfxBfeqkDK22JlkJv2UP8+RdZ4XDlCaly7Vis/rjipZaafKxPq/pivole1jKMnJT2UpnWpOv6/qqf3XADRsXC7h6ocSJXbmL9Yvl9H9fpJ8K8WnAb4LNWP7wo5UaspwVBpsr7EI+XVK9wvCj5mvlQy7CAgcRB4N9LCJQO1PhgzQvAlY2TbLdV9ctl94wW/P05qWoyRtNdWOXoX3qbiwNbKOLBW69Ta+i1JP19Z5ReWnz2LkGGBCF8qNv6Un0qUzYOKIkjXlVT3DhYPW5lGOv3BFfquypdsZDW7uNudHOEhf/uGrNZw6THw/EpYWRh6ewMToHkUzCAMTg5R/aFCR5c6wGugJKp2j3rCgTNw4YFtv0CLgcFPDYLNTPYVz2g7IvjJQXtleZHlHOyK4q/ueqyAsR9gDDx/PuwqgQfCtP9wTioUeuHXUPUt0xbj7gjlPwU/90BERS0f9YQDaOACNk6EZn2qHC/3Ebn8RjFboNMlwU+uBabxB196SrjNnUwX03jmXZcPgsFZcMP7sC0Mew9JLhjTBN7Lh9IQC8AS7ZdnK3SemQO/IVdED1x1OteLLyql15yq9f2M3jL2acFY+aHt/ooymVJxX+WyNCFlcUoreEvljLX7Dsz90HYIrP2A9Jkqoa9E8vGDpF/MSXwRwF6B1rfCQe2hdB1Mk011Wifjkja3l7E42qP+P+4oBjCAd2JiuIoZbFJxSTqOa+WC7oJnniU1tmUj5Xub9j8Z1/THEqVpyZA2yFFWMdOU6KUzqOSfS5eHS/lnixi+L/Hyn2HlDGkrDbp+yj9K8JKnpEb31VOnyQbvsJydDEluwyuFlhlFp4rDY9tK4ecWpM7k6SHvrYsXHQ8dtmJen8/j1z4jjmk9cSA/Cq5rND5813WC37rsKcELW8j7cl6sNJI8LF/+777LlGZqQz9RN5f61GYHeLnpvu0TKgeuCMaBNePKzYGyPGgdelhEKW+CcUHrC0NuCyCXUh5lHgeTwuXUGMcb8Yh1w/tp0ZQBF+WFL8/r0iYwv8wyIxwFpw86HtZNxPcpdSBggfJaPuoJB9bAZT2weRK0PYk9R3jUDl7WwvbvoM3lYMKTy/QTf/AxaziTDA6je/ALIhT/GQR9Y1xcuKOc9SGu6CpxWCz0izO8HEoVMj/S6ACJrWGdtml1APgGrtJaPuoJB9bABbDuY0hoC83D4Ly37l8Q1w5anhX83Fri3yzmV3I5m5PoSNuwtVtXuC8bzs2EW/PL+bwkND/5QNyQAvkey2thGLg6kA1eD6z8LPTGGiMcjUtibUI7/ta7Ko9r3PYbxHGjHXiV79GCS6WLafenpQYWrQsv6G+EHkDhZ+Atg4PPgLNmVD8egNRBMkrzIWUePrDXNI5lDbbbvaScJGN5tN9Vq6tkLM7Ok6SlTfOAwKZnmMjjDOESTudFXmEbubxSKAOlDuku/fkXLegrX1Blj23frDYktPuo5jlK01IaVjW/rx5wQXu4sy+8sgYemRYNVM1EJ+2SO7ApF8oGS8aqJdvrVTPig1LgzGR4m0106uSLPzsfmUz5igokO1QVMwzMXexCH5a6NvLU5bcAUITczf0EqZ89zC2CL1VL+T78Jni58kl5J/YcwbX/fbRKCxz6ltK0TpAU9Z3bcnvVvRUdjsHE0bgiEBUFsO1baBUOo3bLCj4khU7Q5Pjgp9cSuyjjDcbjxculXEx6mNxS9yfOaguv9IbvtsJfwxyIfm1/sFjeqTZ67j0SaEcynZjD8jD0rJGiAcy4DryBC2DTB5CQAYf2D7mp9fxACbnQWpcnDg25bOcVf1niy7mEI1yRm4Jx3RHwzmHw03Y4bRZUhFHvbpkEl2fDt+SyJQz5nK0ZgsXL3GqlehzshjNwRSg2fwwVhTDyopCb8lLOEt6B5IGQcnIYOleFP9jGS7xCMSV8HpfKye76K07wZ3AZeGYY/GsofLwJTvgZdoX5Zr5rMMRGwQvV1rL7AkNbTiSXX9lRba3rQCDCB646rauYmd3E/iunylpmxHy1uNflqnIk3X6b1F1+QsYCjVin2tMe9IFW4t1fw9v2NHLpDP7yVunj1c2sNLan218u+Gv4Br4oDDkFaeApgS97g/WCzgZS/vj+ydRunHPSfwT/xlMlbDTDxQR3HJmk8h5LeZslrFamiLo459ZSWUpN5wJWGwd0ncU5iisPxmbHwusnwIjO8NgsuGWpCizQLhA602qw4moVmH76WjoRxTRa8wa7eBbpwb/OI43cty+WGxlpWbL4xUD3VA6lJfcxhCf4ibyAP1jXrHyAOwQPrC0AcPx2qUEZ5awx61QpliYjbVqzJqicXCWPcoPi2rhX2wIFfG6yb4GcFaH5cZlO2ZYHcoKfCHCuCerHZYwZjs9Nzw28Yq19RB1PwfcuHIRPd3/cWvtatYYCcGDOuAA2vYmLJsRyYvBzg6ACC3Nuh6Y9oWOQTPB9QC5ebmI6E1nHuXTlXwwko9rIWHcYlQ0LLoBhGfDX7+HmqfsnGup2UinB8kS1T+q+4Tg6sYsyZoZl9taIEcalojHGDYwDTgS6AWOMMbqw6FXAQmttL3xfaU8YY2pcXhy4A9f2KXjYQBxjw9Pe2g/gj5+hzyMQkxqeNgNQiocn+I37+YU2JPIJh3AJLYkOq/lzzWjdFP57PXx6M2wshOzx8MJ+coQZSCwjSWAcBWzTu8X7gESiOZJ2/MgaykJNF2rsCK/GdQSw3Fq70lpbhm8uPvpPXjHZb1qahK+gV42tH7gDF15KeINojsdFRnia/OVvENvcN3jtJ0xjE5cziakUcANt+YbujCWduP04gLVrAY+fD0uehhGHwa1vQ7+3Yd624NfuC+Jc8DhprKSc5whD4BYwmAxiieJbVoSlvUaNvRu4mlcWw/E/LlettQVh5bHe/7tAPAccgm/ROw+41lpb47dVnWpc2cnG5hwW8Iub1AmrFFf16DaMkvlsbedIn2HtIX/4eJXc2FvS3B6ZpLGQYp6lkOq7gjM5QvCdqpxYHqm7nzfH55F+KBfTlVOZwxXkBwhFR+X/Iq79X4qMawpsC+AHjhW8s9q+n0dPDqUlf6E7PWjJNip4hTw+ooDNeGiG9Gyf21OKTNGT5YDQrZmsCzn37SM5NBWuPwTOyQBj4b35cM8kWLEd6IKEXn31VlyFzFW7PmDT9P5ecEcPeJb3WIovbuvFXBkYV16iVhKTpf7Z4VxpIPZlRSaFwIB8n7dOarO83ce2zpf13D7tMUzwgcrauel2GXBo9AblHMW11qr1visUV2ma6PjGJEk/z6ry4/pH9iyW5RSEpnF1yLbcXEuN66qaNS5jzJnACdbaS/38POAIa+3fA845A5/T8vVAZ+A7oFdNRXkid4+9DuBlE2V8QRznU8iDhCOHYQFv046jyOQWZnMedj8mdP3OFn5nC4fQnGEcxS005xaa8xsl/EgsU9jBSoprrT9FYziMJgwklX6joEuyz5N+3FL41+ewNjxSU43olgI3d4OZzN89aIWKbOLpFuXmsl1hrLPWmBHeANT1yClIO6pvw10EPGJ9s6jlxphV+NL4f2EPOKAHLoBiXiaW0cRxNiW769vuOzyU8ivPM5B7yOAyVvHvMPSyZixiG+PZSGeiGU4Sw0niGg7iGg6iGA9rKOG3R2HpGigqAY8HTFwMLqCdy9DJ7aIbvWhLLDG4KMHLxAJ4YhG8vwZ2lFF9J2s/INoFb/SH/DL4JO6H4BfUEhfTlDyv5YPSeswKbkjwAuEb42cBmcaYjsAGfDXTz1HnrMVnA/+jMaYl0BVQ1jASB/zAVc4UyplNAtdTwngIg3C7mdls5H+05zzyyGHHnr84wooVlDOOHYxjB4ewk2ySySKRTsTT79BEzhoOrt2qpm9ptctaVni8rKCIyezgVwqYRQEzJx+xx9fZX7j/UMhuBqdOgXaDwvPJySSGk2jCQyWl4TB8OTBgCcfHwNeUtRXGmKuBb/CFQ7xqrV1gjLnSf/wF4H7gdWPMPHzuBzdba2tUUOt24GqOrLfXSR1X6XHfDVG+RONk/MyCq2QDh3+iNC2ZrgZDJE2f7o/bSrsXDplAIVeynQm7jxerHLbJSne6gipfc+2XtZp/kkIPsriXOZzD7BTpj3XmOhk4taa99IjvrETk1WoD4WSV7+ZWd9qbnMdOfF93s4DYdmW4cOHGhQsXz027CQPke3y+7UcMku/dBec+r15fFvub8p/hggd6zP8Z+t8vY+x+ni7/GSe3g5s7wIvb4X/NIalQBgeX5zQRvO8J8l6YN1xWBa+sS3kLh1FEBakJd/BIQpV3tI6tCkQ3pN4XWyoj9nemSDeQin4yhi6toxKlRknKtYpLqzY29JBa7gRGCt5PfRGOGF/13t4TrvJyYQwutdZ+CfKG9Q9Ylc83AsP0dTXhAN5VDMD2z2DXHFpxBeF6S7yUsJibcRNHVx7GRMDk1ouXcioopYwCD+RHSFRARiy81Q5+K4brQk9H3I22JHIMbfmcVZQQBsP7AwVOyk8DwvoHiaMjTattAe07ilnNch4khT504P+CX3AAItkNnx7iWx+csQ5KwrjJfQ5dKcfDR04IxN6hAQxcdTsNcAMBs+AZWXL55M6SUwCdWqHtPGJU0u2kU2XSdNmpcst8+LcyN+P5YRcEMMNFLKItf6OUj4EKTiqVy7Ez8+Xy7rkWVZ7zJ/CNOFa5tFvLVGJ5lwzG4GIB+ZW5Pmo10eEhWQqtw6WSf9dCLkd0io9eKr7jOVfw7ctl6Mxbg86Q5yu9dCsyZWhOaW/ZYbmKJ/0SVe7sUWkN/fMTcmnY44ZZuIHnOJhDSObEb12sDKh/sWuaXDoPuUK+91epTY+uzWRV65+4nOMYyQx+oDff0F5VBQq0Y36jx1/EMW3d7K6QIUWJiyR/t7f8shtzuyrqriSRkqskj3tF8pnqRu/IatkfPWIEnq7Tg/YFXqqHYEQYnBnXblj+4J/EkEEqfwl++l5gGS+wi29I526SCG8idkPG7XTgaFJ4gDVMDLFoj8YxnEgppcxgUvCTHVRHhM+4nIErAIVMoojZNOMaTFhzAb1s5u+UMJtWPENCtQzjAw//oB1n0oJX2MTHujpsiIjnULpyKDP5gRJnL3Hv0QCWis7ApbCNfxJNK9K4LKztWorZyEWUsoTWvAwJe7WJ0qjwyEC4mNa8x1aeDnvCs6Elt1HITnJUxLuDWqIBFMuoU43LJkN5gFXMkctkhq5V7sJGZ+FcKKlOwdGlmXTZe1RFruORNfXexOfscBbz6Mq1cMQ4WLO66oSP5fVXt6myonk1S9ZZ1zrJv/kbAIlM5ErSadtmAuWlp+P1+EqexamCqV+3kJ46w6dLfe6BAdL7RNsPv+2WGtfbXSV/C+liod8rnXJ0UqzU+4q6xguuNbGeN/8u+Aq6YIC7aMFFpPFCSRnXFsZhK22Qj1eiynoZG/NLodR9+iXKkIDK/rfmZBI5nFXcyTEB5eR+ULPcQE2wjQrkbl8q9bA4ndKj4r7HzFea1oXUiDi1wbnhWhn+sJzOgrdRxW+/R7rtdux03+7nntgwTIPCGMe1v+DMuP4EX/EFXry4H3027G0XUszzvIv1ziM69mNcUeG3wYlEuIGHaMVFpPEy27mmsCTsVjhRNKEzV5HHXLbxafALHPw5nKViw0QBBUxiIq7hIzAn6ejB0FFECWUlx+P1/khM7JtERd/vq1ffSJFENK/TnjGk8gzbeJCt++V1uvA3okhiCf/EqZcYAiy+lJ/aPOoJzsC1B8zgJ+yCebj/OQ5SUvfDK+RTXjKcivKXiYq5A/7+IcQlBb+sgSGDZP7N0fQjgf9jE0+GWYivRFMOow2jWMd7FNac5uYgGCqXirV51BOCalzGmFeBEcBWa20P/+/SgPeBDGA18Bdr7Y49tbG7rQqIDkxJuFsd1/eb0qRmZMu4r9ZKbOiwWMY+ocLA9BudmS+F4fvdsmy7zfwamswgaukzUHg+y5QGl7ml6nlPlV+kY3H+tfU2wd9tUeml9gVdcdHr8AvZdvgK3ucDNrOlWvmtGQPk337HuicEn9s+U3Cti+i4rxXKV+Yij3TKneoeKPjryBScdcpz6FzeEXwbzTmUHpzGMEoo5Srm8DsFu5WwNeuzxPkn9pIC4levS2+XE/rJOLmeVGloUcTTlafZyQbm8BEeUnkQ+X5fGZCeBVLjmor8W2NiZXzgUR2lnhZ9naDMnSzf+173KVHsXUlRaTltx8hfXPTU64KXIeMRTwrQ7gDGuasCw/5Q/4d9QiMpT/Y6oBLTuAWYaK3NBCb6eeODZzaUPICJPQ+iTw1+/j5iCZ/xOm+RQDxXcSXHMoiGnP8eTxJncyZn8xc2sZlxPM/vYTIE/DP04SLiacZsnsBTn+WVGwsag8ZlrZ1Kte8IRsNuD5g3gFPC260IQvGD2IocSHwRd4sWwc/fR6xgJU/xHPNZwFCOI4NPiFWFRxsCupLNBdxJNw7hOybyMq+ycz9W1GlNX7owjKV8xHYWB7/AQXA0gHCIfdW4WlprNwH4f+6/T3S9owJ2nQ8mmZavvbZfRfRiinmfD3mbd4miNRl8Tgvuxk3T/faa4UIzWnMKf+NkLqaAXJ7jeX5gCt4w+MXvCXGkcARXk8caFlUrleMgJDR0jStU+D2oLwc4qI06qKyZte3MrEulFfNypcscOUXGgc0apKybJyqbG20NvUVxJU38uvuzsIjmV/6DDs8/T+Y3t8Lkh3y/DggDO/xo+Vrt28tYIP3a/VrIOK9YkXeZwxxepD1/pwXn04QzaMFTFPIitnIrR4ag8cuFUlM7ju8F1xY9OrdyiVvO7nTJrju5T/B5+GxkEknnUM6iI3dQhIcnWcl7bGQXA8T5Wo+MbieXjl9Nl5rWdbdLvXFcrkzw69ZsERdwEdEk8yrjufrqz8Txb56Tdet1zF9pgG50/xz5WgU9pW1NtI7jkmFX9LpdnaBi8rRWu1BNDLupz0X6Q3KG+utthwj+ODcK3ixgw6MM2fd9QiPOVdxijGkN4P+5x/1ta+1L1tpsa212etqezop8bHvhBZgzHo67HzofH/yCEFFBHqu4n7mcQQE5pHAnLZlFEv/ARXrwBvYzUmjH4VzKKJ6hA/15iw2MIoe32VgnX8THMZTOdOFzJrBFF2V0EBoa8VJxAlBprXABHCDRfv+7Av5YCGe9Aynt6uQlS1jFUv7BH5xMOctI4Q5a8Tsc8w60GgR1WJ7MEEtTTqEL7zKSp+jCcazkBz7l7zzDavLrSK0dFR3FIAYzi1/4ldl18poHFBpJOMS7+Io0NjfGrMcXxPAI8F9jzCX4/KLP3J+djBiUF8E7p8NfZ8G5n8D0QeCpmyTeMmaQy6lE0YVELiSp3RjoNAYKN8Da/9GGnWxiMTbMd1Mc8WTSlQ4MIpmBuEmmlFXM5k1WMpnS/bhb+Gfo6XbxalI861nHF3wW/AIH+4YID4eo2/Jk3YzNCQhPssqnSFviNtmo5qI6LusQxacrrqQHq3IfjazgVb30uXa3qTSMaHUSHD0B1n8K088E662mY7w8QBaavWyxEo+V+cT2jjI3L+1FKTIUXignxzNijyWdo2nJsTSnH27iwLsDymdCxSxmJxZTwBJK/WXKfqenuP4s3hf8cW7EhaElzWhPS4aynUS6k8ShGKLIo4g5rOUnlrOITRyvNDSdJ6pjo56bK2vRpWTJ5d0ZsR8KrjWpnswjiSZczDWAocm/u8OuDVUnSEmtOlSUxIwhVXFxRy5WVW2l/Fdt6+nrMSqPdJTMI+Vqdb3SI6t5VerzdfF5devc+dStgt+/qkqjyx4FOfNsaOXJ4rMtHWtZnmxRzeXJ9hcabrBQfWLzlzDnH9DnGej7DMzWd97+h5dStjCRLUzERSzH5SVB7MkQdQQk3Mph/n+tlzJKyKUHxRSSj4cKLNCMgYDBTVNcpHM/nUgmAbc/UNVDMcUsZRP/IY/JPMKQek2iiSGWs7iYOOJ5g3FcFjhoOQgvGoA47wxc+4rlz0JsW+h2M5RsgZ3311tXvJRC2UQoq5Qa45nR4mKakEkc6cTSHA8ZNKOVf2AyxNICA3jYgYctLMJQQCGbyWU9WxjGMxAQymD1lm8dIhoXf+EiWtCa//IaW6qV5XMQVjSAyHln4AoFc2+BuBZw6H2wMh82PVPfPfKjmDzmk0dViEawpeI7aoud/Rh/tTeIwsU/GEwG7fmEt1nhBJnWDSJ84KpbjaursTkvBfxC6ULaL3tGC5mfd+RDSotQuYPVNDBt7KBiqWbdJOO+MpS3d/p4GU9jlTZRlOgC3MRGvUuU+3TY+jfI95f1UiFka06VIQwHbZd5lUZ5fXGc4koTe67NJYI/w98Ff5A7BNe5lFpD0n+7bk/HdWmN6yh+EjxHna9f7xUuFVznZi6kG25c/INBHM5BrOVetgUMtn03ymCpxW06CJ41c43geunz66AqgTQVmWbbaavU3za0kHE8q1SpuF9UXuqlHmki/5u7j+B5pAo+evq3gq8ZIO8V7XWm9cQWAdFIQ7OLmJPjCU3jis62pNZS49rmaFwNFB5KK84lqjgKWvwbTCzkPVXfnWrwiMbFNf5B6z/MoI+aITrYj3CMBA8UlMOmM2HnfyH9X5B2X/BLHOwRMcRyC8dzBAfxGjP5liXBL3IQPjSAJGtnxhU2lMPmMeDNh2Z3Qu8WMPcqsBH+1RVhSCKZc7iIdFryLD8yzfHWqnt4qVeTwNqgbgeuGISutaGj1A60LtJSJxPKKvCgJza6fp3SpLYkSu3g8MVSiJqUJesyNhsrTe963SgDwxKvrxKwCzpVxqBdRaw3l9iOt0Baa1hzLnh3sQmZkLY0TeYG9rxUalDbTb7ghVbqcVff+B/Jb5L80xayGMdPHCW41pzG+T3xK6H9vL6cf7rgX/eQsUz6f3csPwj+mvLz+lo5JS2kG8kcxJFcTwxNmMC/KWfhbvVIa2JPpsigu21K8Lyi3z8EP0lWgGf0qgBd6QNxiKdvulxw/V7dsUx6oR298VdqQuwgGUQ2epzUtHQM2o8cI7jOK03PUW4bAf6TUeEKY4jw71tnxrUfUGrvInbdWmj3DBz8C6w6LfhFBzha0Z++/B8eSpnGjax2rJfrFxH+9jsa1/5C7vOwYii4m8HBv5DG/k/Mbphw04zb6Mc97GI9U7mWfJbXd6ccRDicgWt/YtdkWNoXSuZzMI9zENdjwmE70kjgpiXteJ80/soqPmMaN1C8nwppOGhcqNM4rh7ZsfbDnNa7eTI7xfG2G5XRqs49VLFMn4+S0dzHlk4WPPFjGUS5YYyuXyf9vbQO1FXtZg3kR8HTpwRoDbqvgfltrhjyL36eFC6mjEVs5VoeUhrPjfxT8E6PKasWGSrEskHSnaKlR+qBpW75Zv1PmdR+g/Sr+juyFJuObdJ1Fz3Kw77TRtnfBW2k4PgbMpapJ4dyCNfgIpYFPEGUStDrvkyJ8qoWYbnSO6Pfk3zDFfJ/retEjp30URXRvmzKH2v7IJVHeqMSklSqItcr/qLiul6o8qS3j0r+e5rytJ8gtdbvRh29+/lV2XNYmrMztDguk22hlnFc1E8clzPjqgt4y8jldjYxFhdNacuXnEx/3Afg2x9PGoO5hV7cRREb+Zkr2KSCWR3UNyq3FSO3PpkjztchipnIegbTjPs5mTPpRRfeZxIraPwJwwYXXTiOPozFRRSLeI41fEykpBY5CESlk2Dkwhm46hhe8vmDaxhPNGdzHDdwFrNZQhTvUsH64A00QCTSn5N4gKZksJl5zOAF4lhQ391ysEdEfpZ1nQ5ccYVlModMJ/kra+cFZ0qdJAZZ727EJ5ME//RUGbs0+kgZL9N2ldTQzu74gODXIJOkc1VsUPonKn6md8Dz8+QhrWO8z9mC37rsVDDx0PT/OKzpzXhdP7KDF8jjZbwUgAyTYlY/Gcf1L2Sc0lFumSuoa+8dwkLBV6mguEEfyNqB1XQb9fec0kLW77u3jSySuYIuJNKGnlxCG/pTzno2cBU7+YIsquub96kXOD5TLh97M0e131kev0Ieb/ui/F+PHfKR4IUDqpbpv8X2FseiVBDTka/LHNnFj8u8yBhl9tXpCan3rZzcSvAyJdZmZcq8SqNcknq5lbGcOn789qok3SZhGW+cGZeDmmCLYft9UPAqhR2fpRnXkcol5DMeou+H8k3B24hAxNODw7mWdhxNBaXM5zWieQTr1DxsIIj8gevAU4cjERXr2czfWcMJFDGJplwOvVdBx5cgoW99966WiKIJQ+jMWxzMx7Qim6V8zLdcxlI+cAatBgVLOMV5Y8xwY8wSY8xyY0y14tHGmP8zxszxP+YbYzzGmBpL6zgzrghCGQvZzNVE808y/jgH0i+CFpdB4W+k8xm5fIGXuvG4ry1akcTxdKYbk4mmBWVsZAMPMZvZVER6wpuDPSB8Gpcxxg2MA4YC64FZxpgJ1trd2oW19p/giwcyxowE/mGt1UWoZbt1GcfVPTvOvpNTlayotYTuh6vYHZVruOBeqXltUT5FbZRolnWl1A6ee0F6WF1aKPP74lSsULX63Ur3WXxFldaRNV/5P6naetvPlLFATbfLWCBzrXqtg4C4FDh0DBxxBbTuDSWFkPMVzPiUwr9/BQGxVgtju4nLW6jgJF1XMV4NgEtV1eyh9ylzs8BQoqTOsHkkHHkK9BwEHg8F7h/I4wN2MhmoYCGyPzoObOgq2b7OW52pAteWqP7d+thTsn9rJa3mxabDko4MeC7TRNl+rfxf6fdO63HNkMULtDb64gfXCf7pmVKL7aIyBbqfpj4HcyRVaZfCayz7HMhZGKLnvOluqwWX7RG9aozjMsb0B+6x1p7g57cCWGsf3sP57wA/WGtfrulVnRlXJKMkH355wfcoPgKOuxCOHA1Hn0GCrcBrf8TjnYrXOwM3+XiU4B02NGkP7ftB8/7Q5iRokuX7/ZoF8Obt8P3rrHszruY2HDQghHVXsS0QWB15PdXCqX0wxiQAw6lePqQanIGroWDJL77H81dBZjbl/zwFt+skot13YKLc9MFLCasoZgWlrCWaeVSwCg/b8LIDX5mbPcVMxRJDGtE0IZ5WxNMWhp0OaV2gTTYk+XfFPCWwdSosHQfPfwGbAy1lO/1pyw4aIvZKnG9ujAmcz75krQ30Of6z2d+elnkjgenBlongDFwND9bC0lmUe2ZT7rkTSMJljiA3ZhSJ9CCeTFI5FpfKibR48VKEpRiLF+MfxFw0wUUiHfTr9N0F21fA8q9h4y/g/gXyfgev/4Z2ikc3YuzVjGtbkJSf9UD7AN6O6oFQlTibWq5R63TgKiWWFQH5gSPzv5InqDApsiRNUkshrXFpXWXjCzLp7OrpUtPaPCBF8B8SBws+pkIW6C64VA4GWd/WEJOm/O3TslV+m9YptL6WqbgqueuuqJw9FQDf0+7BgLgnlxtWdYDWnSGlOSQ3w5zdHLerCbgSAENBis/g38tOvOxgFSmUs4si/mAXmxn1iCrmpzcF75W0+wSpy3T/Uuk06tb+/TLJDz1Tfsnu/K/0wLp18VPyAlWf4NXnxgh+cY68/2fcpuoXrAuIzVL7VzrO6nZkvF+C2nR4n7MEvxTpOc8Lko6eovy4grke3az4UsUDJbWwJCKE1UlwFpBpjOkIbMA3OJ2jTzLGpOCLXhyrj/0ZnBlXY4TXA5tW+h6VkHow21JkUORKPbI4OIARvjgua22FMeZqfGV23cCr1toFxpgr/ccrh/VTgW+ttfor/E/hDFwOHDj4E4Qv5cda+yVqjREwYFXy14HXa9umM3A5cOBAIfIj5+u2rmIXY3MCbaeU1zcjFVeTxo8vPVHwbir/Tscu3a/ibY5gpuAdVS1BjQkqGOihD1S16kAJ7Tp1sdK8flf80Asl3/BazXFM7cWOMhx+tCrcqNLZeDQIlxbuvHyDlBZG8pngrRZLD/xZWTJ38vBxsj8vXyXbu2yS0sy0B5b6Cl18ptwu0J7yR69TPu+6pubTimudKEDenDZEZifoPEqNroXyzf4x8WjBh16pOvMX1YD2hdd63aCa9Tr9uVnzaFUthZHZO/g9pzzEOK7OFh6q5dlnO3UVHThwEAlw3CEcOHDQ4BD59cmcgcuBAwcKka9x1e3AFQ8ESiN641PJINoD6rQpMu7r00Fyjz/rA5kv+Lcz/y145hRp1Kd92z9TIttDq5SmpfPdAvIRp6ljR6vwgy5SIqoW51Wk8uHOzpUm6mX3y5izacrj/mglq2wYIDWztoNknNTLHaUGlaH0vjxSBW+VJv+Awx9SGpuUvDjL877gDw+5TnBdK7DvVpncqes06vy/o/+rNC4d99ZCca3CBIRaRQ2RObO9tirBUJaErKYXVsvrVDm244fImpQ6p1a/1xe/LjWt5y+8QF6fLe2OAus+7mIioSPyl4oh2doEs6tw4MBBQ0TljKs2j/rBPs+4amNX4cCBg4aIyJ9xhbJUPAJYbq1dCWCMeQ8YDTgDlwMHDRqRL87vcxyXMeYMYLi19lI/Pw/oZ629Wp13OXC5n/YAlDgSUWgObKvvTtSASO8fRH4fG3v/Olhr04OftmcYY77296M22GatHR78tPAilBlXrewq/BYXLwEYY3LqI1ittnD6FzoivY9O/4KjPgaivUUo4vze2FU4cODAQdgQysC1267CGBODz65iQni65cCBAwd7xj4vFfdkVxHkspeCHK9vOP0LHZHeR6d/jQB1mmTtwIEDB+GAU1fRgQMHDQ7OwOXAgYMGhzoZuCIxNcgY86oxZqsxZn7A79KMMd8ZY5b5fzatx/61N8b8YIxZZIxZYIyv8mKk9NEYE2eM+cUYM9ffv3sjqX8B/XQbY34zxnweof1bbYyZ56/inBOJfYxE7PeBKyA16ESgGzDGGNOt5qvqBK/jq+EWiFuAidbaTGCin9cXKoAbrLWH4CtfepX/fYuUPpYCQ6y1vYDewHBjzJER1L9KXIsszxtp/QM41lrbOyB+KxL7GFmw1u7XB9Af+CaA3wrcur9ft5Z9ywDmB/AlQGv/89bAkvruY0DfPsWXFxpxfQQSgF/xFfqMmP7hiy2cCAwBPo/E/zGwGmiufhdRfYzER10sFf+skm3bOnjdfUFLa+0mAP9PbY5SLzDGZAB9gJlEUB/9y7A5wFbgO2ttRPUPeAq4CVm0K5L6B75sk2+NMbP96XEQeX2MONSFH9feVLJ1oGCMSQI+Aq6z1hYYE5KdeFhhrfUAvY0xqcAnxpgeQS6pMxhjRgBbrbWzjTGD67k7NWGAtXajMaYF8J0xZnF9d6ghoC5mXA0pNWiLMaY1gP/n1vrsjDEmGt+g9ba19mP/ryOqjwDW2jxgMj7NMFL6NwAYZYxZDbwHDDHGjI+g/gFgrd3o/7kV+ASf60pE9TESURcDV0NKDZoAVNpNXoBPV6oXGN/U6j/AImvtkwGHIqKPxph0/0wLY0w8cDywOFL6Z6291Vrbzlqbge+em2StHRsp/QMwxiQaY5Irn+Mr2zufCOpjxKKOBMiT8BUOXwHcXt/Cnr9P7wKb8Nk4rgcuwVfMfCK+Yl8TgbR67N/R+JbUvwNz/I+TIqWPwKHAb/7+zQfu8v8+Ivqn+jqYKnE+YvoHdALm+h8LKj8bkdTHSH04KT8OHDhocHAi5x04cNDg4AxcDhw4aHBwBi4HDhw0ODgDlwMHDhocnIHLgQMHDQ7OwOXAgYMGB2fgcuDAQYPD/wNUscwofyaqKQAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAS4AAAD4CAYAAABSUAvFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAABMM0lEQVR4nO1deXgUVfY9N01IIKyBsO+QIQKaqEEQRBDZRhHEERVFwRHFEUYZdVDUcRc3dNSfjICoMCKoqAyoKCIICgoSNGwCBllkk0AigSQkZHm/P7qT7nMb0oEsXS3vfF9/6dP1quqluvLy3ql7zxVjDCwsLCxCCWHB7oCFhYXFqcIOXBYWFiEHO3BZWFiEHOzAZWFhEXKwA5eFhUXIoUplnqxefTHNW4n35L+qJ5oN1Q4Fiucofpjp8TYl/zoCPp/L8AmyJIp4zdQs4ukN6pz0eIdQn7Y1wm98rP18LFRjmlanLvF6+b8Tz6lS9aTnBoCIfXl8wBpM82vx/6gqWYXcgA+PY+ER3N1juby9Gu8QYY5z/9ThRXXvWKT+fRiRR/l4BTW5hStP3Tv56gAnelguYUDVlkD+QSA709uXutyXatl87kOb+TD1WzPPjK5OvEZqNjfgS+n/y/retq6aAOoB2buBovszV7VXt5Jp432/61fg0CGjz3BKaCdisgM3AwDsBxYZYwaU5Xyng0oduJq3EixN8n6L0WPVSHSP2iFN8S2KL2C66/26KAkR4BuyZu5R4qsjziHee/J3xGeNuVQdz3tHTcco2nYfnuVjTeRj4SymM4f0JX5T+vvEt0Y3Ju5So3rsI3v4gN2ZHuzHf1wxqzO5QXOm65q0IB6/MYX4pk7NiLfK3Um8ag6PXOH7+fib4nh//fvELd1F/EhPvlVrHVAjYTpTv4EMACLaAWdtBXbdAHw/29uXodyXjknbib/emQ9z61PMVwyLI37R5B+4gRro/P7qGvi8b3gn0PhJYHYkUOi5JttV+ySmOXO877v3QJmRDWB0Kds+CvUfu5JQqQOXhUVQUcUzQuSlBrcfJSG8BXDsAFCop1mVB4HzBwan98/CovxQJcb9M9/BA1fVFkD2r0HtQhj8lAzHoVIHrio7DKKv91kequWJ3/Rea1xzFO/EtOUXB4mv6hdPfAdaEU+MWEu89328nNv1bAzxLlhNPBkJxe/b4hfa9iX6EM9+gG+FgXOXEtdLpXei/0K8JnhZ61IXK3aIWiqOYxrTgJeGeby6QThfCtRvotbpalXfcSOvX1Z14mvdJWsd8aVxFxKvDlZRDqEe8Tq9DxNv9EUG8bwu3J9wJSOYBOb7oqNRHa1QF8D+uAKsjetdvK2KuvZVE1lSuPVTdW1vZXrR2bw0nD+mH/HBc77gHZRGhR0+71s0Q+F3KSgY6/0o/EfVnlfxcPneCuWQCCMAwst+mAqFfapoccYgDLUBAEY/1XESqjWAORjcGWHRUrE0r2DBLhUtzhiEoRYM8mBwLNhdOQkEiKgHHDoYuGnF9sLxMy47cFmcMRDUQiEyAjcMFiKjgTAXcMgZMy4nQyrTHSKxnZikF30+qHHSpm40YDq/E2sHWvdpgn3EtW5yXAXU7EMT4g1xgHjf1SuIL+5yEfFseEMMOuAn2jYPV5bYV62r1FDbd6pn6GdjPfFUFfTWQPV94D7W0HJqE8XhKP6g0Sz1B631RBW+MX9SP5SEJJxP/AbMJn43XiR+I/5LfNjK+XzA6SX3B0OYzoy9hngCfkRzPI8onIst6IPv4RXJ6qil45+wlXgyziU+YiqHqug4re0jGxFvs49j+kjTApCT4P4p8idEVN+KzXgYqVhUvL296s97uJZ4K+wsfn9P4ipsS8ooUxxXGxHzRCnbDgfWGmMSy3K+04HTB1YLi3JDGCJR6NhlIgCpCQAo8FPvKxf2qaKFhYMQhmowfmHozoHAPXDlo7Rx6xXVD6txWVg4Bu4Zl84bcxDErZ0UBHngApw/MFRu/wwo7yqHM2gQOY55+uuRxLWONAFPE/8yl2OnXoy4m3i1ADfE1+B8iWpduH3fOax55V3mff9ebY67GoL/EddxV1NxO/Hm2F0i74LvifvGkAH+uswPTVgEylaT/8QslZaiVicvf3ob8f4+mgsAXFLwFfH1rrOJ98ES4k0KWH982PU48d0qqO+H7tz/dl23Ea81SqX89GSq4+KW4RI0RAxykI1luARX44PibVoLjd/B6U0FrdWfiR77OMUVbVYqTYulVKzrHku8SENrhW7oCeD8jd8DOd7AtFmJfG/dkTGNeLiPlv9EOUwo/zAzLhHZCeAo3CGh+caYRBGJBvAegFYAdgK4xhjz+8mOYWERbAjCYPyimp2DMLjcb1TCemUjFJ4qnkoA6iXGmASfJwj3A1hijIkFsMTDLSwcizCEoRCFgRsGCd6B60QZ4pXZD7c4X5pXsFCWyPnBAGZ63s8E1PN/CwuHwT3jcu7AJQ4ZuIqWiqV5BQulnREaAF+IiAEw1RgzDUBDY8x+ADDG7BeRBifaUURuA3AbALRoCIrNSo7i/Lbmr7Ou8xUuIa7z23rgG+INIjhw70elA7X2iXcB/HWVs7GB+Ie4mvjCYZcTn7jFG+3yTW3Wx9qBNZnV4OQ6/bvUwyHilxcsJH7IxTpMlwLOm6w1hzWfncM5DqznAtbIlITmpxFpTSsJHKoz/NcPiV+0QWlmf2K6Jo4TS7tO5FzGrt2ZP9DzX8THuCYT//0tF/GWKRxt3koFSx1GHYRBABTChQJsRfvibTr+70gL/pM8b6Uy5FIxcVA50bOGsyY1fDpfq/qjOA+0F9x6YQ20dH+wt4C04Bv68/7ygjq/b+5iOawyQ2GpWNr+dTfG7PMMTotFRDtjnRSeQW4aACS2F1sLzSJoMDAQPxc/J8Ez0xJXyc0qGKEgzpdqqWiM2ef5mQpgHoALABwQkcYA4PnpYK8QCwugEAUQB/sKmOKBK7jznfJMshaRN0UkVUQ2nmR7nIh8JyK5InJvafsY8FsUkSgRd0iviEQB6AdgI9z+oyM8zUYAmH/iI1hYOAMGBmEOHriKfZzCgjvfKWeNawaAkqyd0wHcCWDSqfSxNINmQwDzRKSo/WxjzOcisgbA+yJyC9yr/KEBj5QF+FpaHerJuo2v7gAAQwrmEa+VomJ3VPzM3ubRxJtuZD/fTZ3aEP8ddYjr/EDtI98HXxKfEze4+H0C2DTpQTxJfK66PDrfrAnY27hmBv+u66M5GMjl4sf6tZTG03OO0rR0bt/Xiqv8ubjxbJ3cNoK5n3dapOI/M20bx5ofOJQJe3vyd6evddNk5c2sLK6g/MX2x/L16oMvURNjUR0Z6IMvEbfF+/u8Hjec2t766yw+GIfcAVcprvTB4UmsSeUM4+1Nl6rfpUjyap4GdAXwTlXA53ILu4L7/9X6SrPlkNEkKL8nhsaYr0WkVQnbUwGkisjlJ2tzIgQcuIwx2wHEn+DzNACX+u9hYeFMGOQgzG+EdRDyPWaPkVElt6tgCIDw0q5W81FfRHxd8Kd5dO0KhdMfHlhYlBsKcQwulFxQJagoGriqBbJNqViIAFVKP3Adsu4QFhYVCPeMy8G+B3kea6PqNYPaDREgPLgPNgOiUgeugkaCI/d4T9k/gz2jwieqHXgz8lj2wNHaPO2vl8XawQ+dWNj5BW2J9y/gWKXZrhuIn6t0K50ruRDeZMXzVc2oMfgP8Rm4mfgB5ae1TMWsHY6uQ3zwchZ1DKcG+uXD+flp3af4BsVvZPpbAw5WOqz0wLjdSvPSdxJbWCF6KSf4GSUy6OvRex/7/y9NYM96FaKH3g9y+6Fff8L7T7oQdSCoghrYhyaIO+Dtf6+4ZXwwDhkDHlc8WXFVJm/xS+zb1i33W24wRwXBFulzRz0D19W13FpXEeaq8w1iemSYVyYveLvswaunNOMKEhzePQuL8kM+MuHST3SchGOeIOSqQSlVWAwRQNUDLsOxZA6AXgDqi8geAI/A80DSGDNFRBrBXSmyFoBCERkHoIMx5khJx7UDl8UZg3xkoUpA290gIucwUJgPVI0J2LRCUY6h88aYYQG2/wagWUltTgQ7cFmcMchHJsJQxcFPFg1wLO0PNXBVFCq1e3vQDHe7xhfz1/D3kneYrWht9hHXnvHZUSy8foOLiV+Ld4k/4nqM+FGwKJqgxAyd03YFPi5+/6WKDLlY5VFq6Ny/narmoz7e7p6cV3k93iEe3U+ZRKk4KSj9MOcj5t9EsS7TN5W9x75t0I3704n70/c5bq8ulV9/5C3m5/XkfMCnE8cR1zUte2SwpqVXgCueOo947x3fATU7AfWBXr/+THkeG8CCYexoVUdR53XqeqAqglHXKlDyJ9BdcV+7/2OpyGvZGMeae3Wraj05pm91bf7dzinQgmU5wOEDl5PDiC0syhcFniBdV5BnNCUhczfC0DS4fRAArlK+ggSHj6sWFuWIAs80y3VCIxNn4OivEHQObh/sUtHCwkEIiYFrN8IkBu6kmyBVJBL4lVxzGip14KqPQxjlUyBvX22esld/lr+oiAI20NY+7Frj0jrIZnQgfnPWTOIbotgjKhdVievcRa1D+epU2ttL+1ldtfsz5vnMn249jnia8kHX/lLas75ba44Vat+aPegbpXHdxN1R/CCn7wKly6hYaO1VpuO6oMos5vGlw5LanNCX0CWZ+7eS+1cH7AK+VRl89V6oNC4VvhShjKlmtr4G4aiG6wEkxfSDGeo9/1Vz+btYOpRjxhrGcs3K3b35u+6Ry3pm1EQVp7WcadYiVmhyI7z3XVXsRw0A+a52KPAkfP5Ym4PiumRxHurxyHJOyg6BGZfVuCzOGOThGPJwDNURHbhxkFCAXwAALhUsXakoT1+bCoIduCzOKGQhHdXVbNZJKPTMslx+j4UrGVact7BwDrKRhijUw8HATYMCg8MoxEG40C54nQiBpWKldi8MhVTbsOVCdfuo3MTfJnG+XO/VrGu4urAnldaZRoKDhfpGcb7fNztYmFncmmOZhqewr9Lrsezb1BBe7UPHeBXoS6tCbVIuY43pEnCdQl0Dcr+OWVPJws2V8XmjuawZ6Tin2C9UrBLbgSFH+arrupDfguO6miTw79/0Ec4brfHYUeJfoRfxgu58vXR9gAEpSijqylSlOqLzy2y4Wf0u9/Wshq2ogYuRkOHdnnMZNUV8Fb7P3sgfS1zH6EUtZE1r12Os3bZcyfd51CpuH7WSY/D2P5AOoBv+B7d3/Q0FHNAYqXIjI1d647xcyv/+tGAHLgsLZyEPexCOhshDJPwruzoDadiDWFwYuGFFIQSeKlqNy+KMQm6RtWhYm5IbBhGH8CsiUB01lQNvpcGK8xYWzkKeZ+CSsCA+tQuANE/YT32qO1aJCIGBq1JPXS07F/HJKd4P1MOdI89yPIqO7QGXo/PLMbup4L/E73S9TPwf+Dfxbq2XEJ+PwcShrMEbKx3ra59cyHvxPG1bgj7E777sReI65myeqqe7QPWlfzprPAO2K81Hmx7oXMFkxccorvLpIqcz33DXOcQL1COlpnPVxVJfnYaOc1urAsdaKY0rJZY1wdhZSqNTYWhrpnCMXrGmJYeBWkCVwnbFtQvDVbG9yHuYj9/xKvFNrXm29tqQEcT/NofjBXUcl5/huYpB+8tHzwBDHsSAzRnAppl+3mN+fmG+uY/lMRUpSvlxMKzGZXFmwaTDFB6EVNHVQxyEgmNAxgYgukvgthUBK85bWDgPpnATxNUhcMNgIu07oMUNcE+hCgO1Ll9Ycd7CwnkwBT8BVRw+cB36FgivBdQKQj9DQOMSY0ylnSwh0WUWJ1Uv5jFTM2n7q6NvIa7z86qA47YaqOLZOpcwG9WJ1wTHEmn/Le2BdRPeJl4Hh4lHwJtLqfMmHwF7fbX1pHIU4SdV6HAK/kZcxzm1V4UK24HrFHb8YjtxXUfxueYci+SbMwoANTM4NOBwbRbNtKZVM5evZdSlalYwhenTncYR1/7/vbCMeA9V+LHlDI6FOjiS+7caFxAfeJ8KCvT14aw/Bmj2KlIRh0L85h/zpjTx5V342D23cK7ga3GscXUAe4v1fIHbP6rqNd9uOGiu0e4MoEpboPE2IP127Ipm87SWW1T8o8+tkXg3kLTNCMqAxAZikv5SurYyBWuDUeXHzrgszjwcWwcAqAJdccRByP8FyN8DRF4SuG1FwOEzLjtwWZx5OJYMAAj3r3PsLOQuBSJ6w712q0SEgJGgHbgszjwUZiIf25w/cOUsAVwxCNfr/opGCGhclXrqKnmFiNnno2sp7+5RWW+UuP/hKNYCdG6i9lXSOWiRqj7e8mdZu8hUmtd0jCKu8xGfzPpX8fu3orhuoq7BqP2rtH6mNbJuYH8tnfuo9bYX+7FGtl4tg7Q/ldb33q59E3GX0hPHpqrvhuVF/LCC/7jOe5B1nkue4lxMrV/qa99yt9JxVN1I7c+lr7efAYSqK5m/KBnhuAAuFCBrEP//LqjCU4me/Vmj0p7xf2vIcVvrRitnh2Smj3L4IHCV0tge9vw8shSoB/yOa7EDHxRvXhjH9/1lcQuL3x9/nH3MTgv2qaKFhTORj7VwoQXCdHa2k5C3Bzi2CQ3Ug4cKRwjMuOzAZXFGIg/uWVQVBCnIs7Q48inq4Ry4lBtIhcIOXBYWzkQ+kmGQi/DKns2cKo58ijCEIwbnBW5bXgiBgatyT50FzolTmmMkyyIA22eh0RzWAhrdp7QBtjVHpEqfA4cykZ8W4B/39bWqy3ixii0aETWj+L3WsLRGpTWpb9CDuI6TSgUXdLgxaxbxy6M+Jn4/ni3xfK/gTuLjVW7lc/gn8TTtTKBzH5WGdN6N6stT+X65SjQZNnU+cT9daCHTHLZCw0WrfyD+ZhcumLxoPNcP0PmDMfPSgYvXojoS8WK90bRNx5Q1XMT3SdPb1Y2lvMvim6YQz1PyW7iq07jlo5bE4/rv8hLXtzCfHMb5ea2ArCcAAPui2Xq6aaq3P1VV3uNpw+G5inbGZXHmIn0lUCcRLidnvhXkA3kfA+FXotLmGSEw4yr1wCUiLhH5UUQ+8fBoEVksIimen3UrrpsWFhWAQ98Argg0UhkXjsPx9yFh0UC4tpWoIBQ9VSzNK0g4lTHzLgCbAdTy8PsBLDHGPCMi93v4fSfbGQBy61TB9kHeJUibpb9xg5fUDv2ZHklm25taD3Npcuil5o1M98bxFHu1EmbvU8stHW6hU4Ley7i2+P3DtZ+gbbNxPfH24HJheumUr+bmB9TTroei+Ph11VJQpzvpcAZtI6PbH1f9OS+FL2ZKgrKVGadsZVR4xOcJvG4fMFZ5u6jvRp8f7zONVNbT64bz0lL/fk1X8nIubxnvX6UAgKyEALj2YBqQ9lrxtqzWJYdHQIXVaCy+iy3A+y5Unjsqpehn2UU8jle9kHe/AG7KgOy9Bli+CIdH16HtTbN8flf+2k8PIeAOUaoZl4g0A3A5QAlugwEUCQczAWUoZWHhdJh0mPyNQLWLA7cNJgqPA7vmA62HAGFVA7cvK/5AS8WXAIwH+2s0NMbsBwDPzxOWBxaR20QkSUSS0g9Wsj2HhUUg5H8FVO8BSCUMCGVByjtARF2g5RUVf64/QsqPiAwEkGqMWXs6JzDGTDPGJBpjEqNj7LMAC4ch7wsgrDpQrVvgtsHE3i+BzD1A+5sDty0rQmDGVZpTdwcwSEQuAxAJoJaIzAJwQEQaG2P2i0hj+Kkc/shETXzjE2LQpAsLGZGXc/u5owcS1+EK3Z7ikIPY5ay7LI/jGJ23MJJ4feUFrUMaRmdNI94qStns1PYKClPBj9SvxXvEbwLbSnfNWkW8ZhTbxGhr54Xg/CVtM3032Bpap8AE0rj+ijeJV43NJa7DLWJH8bVe10lpTspqWjkz+1ktN+jCIQeYodor2+72WSkoCcu7KyuaWZy2M3+4uzRdFbjwZ+Rjc4sxWO+Ja2iiarXpkvfQKTuqVFrfHeqX09bNG5kOWqa2f6x4DgAUApv/C5x/Hzr+mg0UePXh5a29v+vRCHXw04HA/ZfuYAScAhljJhhjmhljWgG4DsBSY8xwAAsAFBkRjQAw/ySHsLBwLPJxDIfwMxr5Gbs7ED/NBMJcQI0bA7ctC/4IS8US8AyAviKSAqCvh1tYhBz2Yy3qoR2qweERPYd/BvZ+DdQcjQq1ugmBpeIpDVzGmGXGmIGe92nGmEuNMbGenzpO3cIiJLAXawAATdA5yD0pBdb/BwhvC1TrF7htWeDwgatSrZvPTqxqPkryxifFzmOdJGUIxwo1z+LtG6K45JROa/knniP+Tq7KE1HYF8FeKVr30TY296tJpe/5tRWx1ri+VOXKdIyYtpmZrOqH6bivZLW00frfc2sfId77/E+If5vBYnT/2loDa01cWz3r/k7Y8RJxsISHPGUxFH6E+brmrJHF36s0LHYYwro41X4jt0/vxCKNtgWKWe21V1rXJRZxWIIcpGAnbvc7lk53+qEf56rp9K27VrM2qsuP+f3BRyv+iOIJPu/DwoG7dgPZq4EdHp3TJ70qcQ2QdKSM1s1/EpP0aunaSn9r3WxhETQcwVLURDeEqX8AjkNhHpA2Hah1OVC1VcWcoxyXiiLypoikisgJnxqIG6+IyDYRWS8ipcomtwOXhQWADHyBMESips7UdyLSXgNMARDzj4o5fvmm/MwAMKCE7X8GEOt53QbgtRLaFsMOXBYWALKwFnlIRR0VduJI5O0Ffn8HqDcKcGmr13JAOc64jDFfw68mPGEwgP8aN1YBqOMJryq5i5WpcbVPrGFeS/KWcu897ztuoOJddr0UQ1zns+lyYzoOS+f7VUc28cR0nr0eiGZ/kuvwrjo/R1dX9bFD9i1VprcBwCe/XE28d9tP1bH5d1udxnmUj9Vj4UPHZZ2/l21ekKP+HeYo2UNN3GOvXUdcxzK9gxuIN31O3YsjmerjQ7nWLG1+IXGdy/mj0vAGLlXlxlQ1Nm2zk6O060VR/IHv9zVgqufG6/4K0P4WYF4MUOBzrygbGm1jo70IlyeoGLKbOQ4s/S3W36Knc2k4v1xINUBs394I4fgTmmM50vE8os8dX7wtcSuQlF1GjauDmKS3A7cDAEnELgCHfD6aZowhkU9EWgH4xBjDIrV72ycAnjHGrPDwJQDuM8Yk6ba+sDMuC4si7PgICK8ONPlzsHsSEHn4GVlYhNq4BagWFXiHU8GpzbgOFWXGeF7TTnjMks+mEXA2ZQcuC4si7P8GyP4NaDkscFsH4DD+Dy5EA0P/FrjxqaLyAlD3gMvmNIO/baUf7MBlYVEEUwD88h7Q5HIgXK8HnYdcrEU2lgEjx5fvrKtyA1AXALjJ83SxK4CMIvOGklCpIWQ1D2ah91SvrpU1ksfNqO7sHtFyqSpRxaFF2NuaA2DiktjX6NtEjlU6pISQ1dEsTtSjpbp/bNaevRx7dU/TSd5zgc+1+gAfe2xbtkr+Be2If5/F7ZvXY2FFl5h/JO0x4shUyWUctgVcpDifHu1967gjsF64YTyXPxuwmwXKlN4ck6djne7O5dzK5yPYOvrWubO4g7q0oFZLlKYWqf5nD9Y+4L6Xy9cB/Md3gLPvwoGrR+GIJxbvEfC1fhr3E295Pd+ndWarEmFKb8sW1rSidSzps4pPZtpmoo+PXdMJwIjvgHfvAPY/7681ng7KsTyZiMwB0AtAfRHZA3eUWjgAGGOmwG3SfRmAbQCyAZQqi9zhdmEWFpWM39bgOHagJoYUD1yOxt5VwOHPgcb/BFKnAOqB1WmhHI0EjTElrruN++ngmJLanAh2qWhhoXAUH6I6uiMcLQM3dgL2PASExwCNxwduWxr80XIVLSzOBBzBezAoQC1cF+yulA5Za4FD7wCN7wnctjQIgYGrck+dBcphi4pTjqi6N7WYHmnBnvNf4RLi/RM5tklDlwz7Ge2Jb8A5xA+ksa4zvCkHt/jmKvbAN7StR0Pmz33KcVjRA/YSz0zmmLV7u08irv3x8z5QF6cOUyg7f7+4qlZMl2X1It44ikWi/+AO4jqO7PXmnBeq+3sj+NrdGMH+ZG3xC3coS/V3uuL6XjlfcRU3tqW5KgE2zkcPVV5fXyIO/bAeMRiO5UhCVRWjd0ylBT09exzxCc+9xAdU9lyzDdfJG/+CSgycx3TNChb0Oi/gL9PEAgh7AJC/lJs2ZRxensxqXBYWJ8AWfIWWOBctkaDHHWei8Fcg57nA7UoBEwYcd7iRoB24LCxOgN1Yh0wcQgf0AdRs2rE4pm0lTg9GgHxXaVWk4NSRsBqXhcUJYFCIzViKpuiExn7r8D82jAgKqlQp1StYqNRcxcRGYpJGeHn6syXPR2tmcLxL+FTe/vl4zuT/Xd1grbGTeNdkzse7I+EF4jpfbh/YrytNxYEtwKDi99rjffk6lRD/EFNsU7yV4vqe0Dqxbq81oDjFlSdU7IN8LVKmxpd4/Gv6cwn7q/EB8fdwLfFDqE98Jm4i3nIjxz7t7cQxeQ0yOBcyXIV1IZnp9tcbEa+OY8QbLc0g7hv1valnG5wILtTDn7ACYXv+A/zidWLIS+B24VzyEuituM49fIDpJ4N4B10L4ezcDcRzI05ekah3Yi6SkwrLlKt4bmKYWZpUurVitByzflwWFk5CAdJwFIuARiOBsHLOB3QwDAQFcJXqFSzYgcvCogSkYSZQpQ7QqBLKgjkEBoJ8uEr1ChbswGVhUQKO4UcgYyXQ7B8IalmbSoSB4DgiSvUKFipXXasL4Covjchlz6rDEXWI1wRrXJvGsxbRXBkldSvgOK3fXXy8+QmcFPZ8LufHbYvg3MQJymNe61gHN7Uofl+/49e0DVuY6txAfKLqCG75iHl9lfFfR+3fSPH7FR+o+DimKbs6EK8xnDWnzBkcV/Z+3PXMM0cQ79nxc+Lar3+Zirnr0Ymvl5/XWg1OXdkwhi9gK3CNSx2DN3iOyk1UuY5ZZ3n/Z7fK3UnbohaoJ2UXTgKazQPaDQGOfoBwDutCyiTOy4zdwrUS/PJEWzAduIC9xnIu5e1To24jfleSco7xyWWsshNlRtFS0cmwMy4Li0DIXAAcTwGi/xm47R8EVuOysAh5FAJpzwPVLgCq9w12ZyocVuOysPij4MhMIO9XoP7Dwe5JhcO9VKxSqlewULlnLirt7cGBiAa0eS04HGTou2wq1XY4G40/FTWB+OMZTxM/EM26ifaMfyjiSeK61uE/8G/i2nNrYMe5xe+3qrzHsF6cbNf8Wtbjdh1SgVadlKbFdvf+uYZP5jGfznmcfn5cKo4sfsxa4uuu6MoN7lX736+Ozxb66NaR9cU++JK4jvPSeaNxW9hLbW4ci3Q/gTU5HUc2+D6laemaFw8yjRrt1bGWD1Ie8RHsEY9MADgO7HsOaPkqsPdiIMur0R1uUIeaPxnHv6trUgHxCbtf4uOrvMxI5bbetqcK+ktm+ttbXtPDvMRMlBVucf7ksWJOgJ1xWViUFofeAPJ+AxrpaOI/Fgxgl4oWFn8YmBzgt+eAmn2BqIuD3ZsKhPOXinbgsrA4FaS+BhzfCzR+Ktg9qTCEQuR85eYqniUmaYbPBym8Pe8K5uEr1QG0zqMzpNTyfuaga4g3BMdObVPBVToXUW/XPuxUK3EPp4fV7sqGWBmrVOCVyh2EzsXTcV/NFD+seILinCqICzqyJ/z3L6iKzTrWSF9rJcm16b6J+PZNHYk36sh65J34P+LaQ38y2KOq6WrOVZzV5S/Ez8WPxDsuV4UWVYZOjorjipzoQ1RdxJRBKi5rrIrLGjEa6DwFWHYZsP8zQKU6ftSFy5tdlfIZb49V2xfwdr9rP5LpiiZcpb4hUr3HSjyADUnHy5SrGJcYZaYl+ZVAPCF6yvc2V9HCIiSw/U0gczsQ/xROXBYwtBEKMy47cFlYnCoK84D1/wLqngu0uiFw+xCDgSAXEaV6BQsBBy4RiRSR70VknYhsEpHHPJ9Hi8hiEUnx/Kxb8d21sHAIds0B0r4H4p8GwqoFuzflilCYcZXmsUAugN7GmEwRCQewQkQ+gzvrcIkx5hkRuR/ubLn7SjrQoahovNmlfzHv1oVjeXJV7Eh8CyWCcSgUlGTlFw8zYvn7xHv1ZC3hTVXC7RFlnJQKjjP7YuVgPkENn/eHeVPGB0rTasXUTwPbqdrrJ+63K67DdWoorjS075OUpqWhYofQi+mFHTmf7ru9HNOm+9NdxWktUwfUemHTZNa0NnVh4UhrWjq/MEepLBuilE/7UiUcjfJ5fxdvip2sNC2VO4iFAGCAbfcA//oGBWf/A4W5XtHsqhv5Pnv17VuI6xg0v/ta6XPLm7AeqAcMX27KYen6h8hVNG4U3ZbhnpcBMBhAkbvcTABXVkQHLSwci59XAEnzEBZxPyANA7cPIfwh4rhExCUiyQBSASw2xqwG0LCoVLbnZ4OT7HubiCSJSFLmwZwTNbGwCF28Ox5ABFyRuvx06CIUUn5KNXAZYwqMMQlwP5S/QERK96zUve80Y0yiMSaxRozDS4dYWJwqDmxDYe4khFUdAXHpmJLQxB9F4yqGMeawiCwDMADAARFpbIzZLyKNAZ9gkpOgEGHI9qlJpz3d6yih6NVOrA2MrfIGH1BpAb81r018pxKWPs29nPiPEQnEqyObuPbfOtq9JvED8C4PHovnCitfKmFk5lTORcx4UmlaejKq/wa0R72SYTBF8ZGKa/8uXXdRH1/dGVU7sndaZA2+VkO6zCauv9uaqjS8/q4XJ/Av3HcWFwXbNJw1rw8iOFkyMYJFus67WdPa0pvrKvrqqfHPKS11DlM/P3+VZ+068hQQPhxVIicD6edh6dusSXXBauIN05X/varreORBzgvV92XnJKXXveh9G8k2ZacF91PFEM9VFJEYEanjeV8NQB+4bfIWAChykxsBYH4F9dHCwuHIBjLHAVXOAaqNCXZnyoxQWCqW5syNAcwUERfcA937xphPROQ7AO+LyC0AfgUwtAL7aWHhbOTOA3I/A6KeRARuRG7gBYij4fSnigEHLmPMegDnnuDzNPg/KLawOHNx9A6g3kb8CfdgQ8mRQY5GKIRDVOpcryaO4BJ8VcwXoT9tvwIfEx+7gDWt7YNYqGmYy//VGm1k7aDRZq4duG5oLPHpFMwD/BOTiP9X1QLUOo2vdvGWEpVSoR6Pa41pleLKDgsfrFcf7FXtOd/NTxObofiTirdSvI/iKndy7fDziber/QtxrVnNWftX4n8+nz31zwbXCmylamCiO9OOqZyLWK0B6z67tZH7Im4fl8t+X7+N8eqhi+OUvtaP9TW/v+GFTHcNL/Lnz0ItPI8YPIoW6ILfPZ5knecqTUr9bvpW0WirBUj+M0H6bO9Dr/xEZYh/GrADl4XFGYYjmIqquBYtMAFHsBoF6p9dKKAo5cfJsLmKFhbligLsxKMIRzRa+JVeCg2EQjiEHbgsLMoZ2fgJ+/A66mMQ6qB3sLtzWnD6wFWpS8Vq+46j4yNe7WHRY6e2f5vbOfho1xSu/ReVxbUBj1zF8TDtclmXeTHibuLj8TxxrcNMOsBG7M0bepPMpoDjtLbiT8SvqLKEuN8/42TF7+c6gZik+JWq/eeK6+NrDUzHcWn/L3X8zF58rQ+vOUxcx8y9fz6bq92valRegQXEn1VidpfWHPtUFRxHpjU1nfvYZBTXdYxdzYFvjZK8emijG1nT2r6ZBck2c9TFGs605WS+7+refhjAvSgI64q2eBBbhn6LAnhzMTvOYf1t7+ho4k13cN4mcrm+wPzHuD5ogk8eZ6Gf0dupo6jKj5NhZ1wWFhWCPBwr/CsEtdHE78mIsxEKcVx24LKwqCAUYhNyzcOohf6oi+uC3Z1TQnkuFUVkgIhsFZFtHicZvb2uiMwTkfUeC62AKYV24LKwqEAcNy8hE9+gEf6FCMQG3sEBKCpPVppXIHgC1ycD+DOADgCGiUgH1ewBAMnGmHMA3ATg5UDHrdS53u9NamPuYz2Kuc4FdOn1uQp/WTUlnniHAt4fSkaqFsfawLEarHn9rGohPg2u0/gxWKe5tyHHeR2FN3cxdh9rKFc14bglxLG3/1mXs7/U5gnsI67Dmvz8smagZNRhGt7qCPEe539NvDqOEf8k8ko+gNLI9uxqRfzFlqwX6hi9Ovid+B27Xif+aEs2ILtMBUt9jR7EE9UFWaC+q9ilKpmTUyfxeZzXn+yVzXfStoX72N9ehfshJUF50i/ncx1wcWDWPkxCV7yFhpiM73EbOkaxxtV0I2taBzuxudojYDE4Ahyr5TuAZOEblBXlrHFdAGCbMWY7AIjIu3BbYvn+8XYA8DQAGGO2iEgrEWlojNGOe8WwMy4LiwrGcaRjE55CTbTFn/D3YHcnIMpZ42oKtkrc4/nMF+vgNiaFiFwAoCX8y8MQ7MBlYVEJSMP32InZaI4rgabDgt2dgDgFjat+kd+e53WbOtSJLFl1abFnANT1eP79HcCP8K+DRbCR8xYWlYRtmIba6Ii68a8DR9YBR38KvFMQcIopP4cClCfbA6C5D28GgGJVjDFHALePuogIgB2e10lRqQNX3cwMDF35STH/vDv7oGtd5OIHeL3eNZVzDzFXnYBDd/BtbfZFOhuc//cVLiH+pUrY0/lz2hfpuXVeD66f4llvrKp0iNi2fO7NU5WmNYApVLqcX5yW5r2YRt/OuY1Xu9jn/BulGR3ShRjfZT0QygPyvLu4g48rkyqtX/6UwddnXEvWC3Wc1ju4nvhaVURzSJZyUVLebJf25oS+PkoA9f3D/ChLaVosB2J7Asd17aa/QyB2IWtcO3u2Oum5kvE4LsmdDMR/CLzbGcjL9NPfYjqzgf8zq/hB3Acu9iLrhWXF72uqe/R0UM4a1xoAsSLSGu6E2+sA/nI9tlnZxpjjcCuKX3sGs5PCLhUtLCoRuUgDPrsOqBML9Hkj8A5BgPupYkSpXgGPZUw+gLEAFgHYDLct1iYRuV1EikrAnAVgk4hsgfvp410nPpoXdqloYVHZ2LscWHk/0ON54OCPgMoqCDbK2x3CGLMQylPDGDPF5/13wKnFitiBy8IiGPhhEtDgXKDbU8C2TUDGx4H3qURYWxtfhAG+s8sGyiVywNLl3F5RoyaQcpk6vvptei74nniOsj2csPEl4ku7XEh8IfgEq9GFD7DMp6s1WJ9r2XYrce3hHTP6V+LtlOfS0Z7sb79xUWc+t/bvUoHZDV0nDYEBAGx+VmlsWlPTdU2U/PrDv1Rgl451vpbpY7XZk38nWhPfgLOJN9GCpcLdUS8Q19evPX4m/m/8g/jnPqLix1EDadvQRZ8Qb+PiXMU2E1Xu4lNMuxRwnuXXLtYTi6/tN7cAdWOB1rOBpd2ADE9u7FvcvNYqjkf8a74yxff5O6myH2WGzVW0sLA4OQpygEWDgbwMoPvHQMQJK/xVOmyuooWFRcnI3g+sHAxExADdFwCuasHuUbmm/FQU7MBlYRFs/L4WWD0MiO4MdJmNYP9ZFi0VnVzJunLneukA3vXS845s5u1almHbIYgqfweWibBuND+YiN/CO0R+ofZXx+u98Dvi2Y/xfz8d17Uiv6/PwTlu60AGT/sX1+5LXOef6VzBpDQlKuk6inUUP8R0dxbHGr2VczM3SFb7f6m4dmI5HGC76u7GZqzJvdedr08TsBjjm/cJ+Pus62vfSsUn6jqOupahjsH7Cd64suGrP6Rt+vlW2lnM6+mi1RMVVynCA7cs5Q9YngNeApCxAPjpTqDjq0Cd/wMyvGXOnus0lpqPX/Aq7+9bTrScxpJgLgNLA2f3zsLiTMKuyUBkc6DtfUDBPiDzqcD7VABssQwLC4tTw9YJQOPGQK0ngcI0IFuXKK942IHLwsLiFGGAw7cAYXWA2pMBkxFwj4qA08MhKnfgagSYB7xU0njz5705d3HAVSqQS2tcKpYpvqdqoPy5VOgQ/OoYqNzHgS+wNnH+PZw72fseb7zP1eBcQJ1bpz3R+2MRSsL19d5hHvk/bqC8ynTcVWYOe8T7tX/3UfWBqtO4TMWsaQ3sdsWnqIT/jWwKkNqdPap+yWWT++zM6sSXZ3LQ3eCWfH21ZvWFynOti/dQEkjX0veVqnn5XT7fKAPb8H2hx5ZaD3PclSpdAJXWCdRjaiLzgexrANdnQJ3/oh8exkEfn62XB7EBw53p07xExYCdDgoRVqp0nmDCPlW0sHAkcoGjg4CCJJyDx1BPBz9XMJxe5ccOXBYWjkUmcHQAMrEd8XgK9XBB4F3KAbauooWFRdlgMrAW/0AWdiEeEytl8DKA4+O4xBhtRlhx6JBYzcxK8gpN2eA4qYs2/sA7KM3pp8fV8VQ4C1S8jSrVB7zPdEvrlsTjHt9F/NWHbyGuY41GYXrx+2XKEEv7S+ko40ngGo3nYy1x7b+vNbPvLlEC3U6mfnUUMxX/RHGtdt6r+LuK91JcxZH5aWLaau5KzZUudJj9wBrFs097fbBAqmtg6uutPey3+tQb0P71uvz88GQV56U1MSXF6vhCPMh0bxdVR3GiqqOoLfSuABBeF7jwS6BmR+xyjUEmvjphf4cn7sBPScdO5DpaalRL7GBaJ70TuCGAzXLe2gBGghUCO+OysAgF5P0OfNcHOLIeLfAaavk5T5Yf7FLRwsKi/OAZvI5hPZrhFdTBXwLvcxowEOSiaqlewULAgUtEmovIVyKyWUQ2ichdns+jRWSxiKR4ftat+O5aWJzhyD+CnRiBLKxEUzyHerg58D6niFBwhyjNmfMB3GOM+UFEagJYKyKLAYwEsMQY84ynOu398FeVCNW35uC8Xj75iTqjQcddsSwBpQRgD6dwoZn25xrK9OnW44hrHST7YY4l0p5QPcC1CNN8AnC0x/oYTCZ+iY8mAbBPOAA8t5X9qtq030Rc1yXEo0z9NCitafVSXMdh6ULLyYrXUHy6sgS/shbzR5mGj+T2eQNU+3asaf3f5VzM8D1l8LViOed+bszn3MiJl7L/lm9uIgAc8vnutqr6mg8lq2TClUzB4YbIGsT//6NmFHIDThtFndzD/IFWiLozzVHb52AQXFiMoYhGRzyEQ+iIZLg1qRyUgyEXnG8kGHDGZYzZb4z5wfP+KNy+0U3hLuo409NsJvzlVgsLiwpCAQrwHubge6xGJwxBN/wdYeU0AwoFjeuUflMRaQXgXACrATQ0xuwH3IObiJzQBc1TZ+02AGjh7GBcC4uQgoHBx5iPqtiOBAxDNdQFAmRklO64goJCZ8+4Sj1wiUgNAB8CGGeMOeIufxYYxphpAKYBQGJNqbzYCwuLMwQb8SGycQhd8bdyOZ4pFOTmOHuWUaqBS0TC4R603jHGfOT5+ICINPbMthoDykD+RGgK1rV0/lyC4kqzukj5HEHliOUoK+5IdfwJyS8RX5PARum61qALBXz+HRxnluMzx0yOiqdtn+ZeTnxbRFviH2MQ8b+1f5F4e7BnvdZofkpgnjOdY4NqTDlIPPNQHeJYoeomXskU/1NcF6K5XQk/+SrXsRXTvIFK01LHq30R+7hrj3itR/6l5yzi2o9rLc7n9iqXdNg8b13Gz4co0UrZxOm4rfQx2pCf8cxo7vsTG58mvq0J3wvxK9UJVE7t4Z61ievjuW+lD4GuHyF6w+ES+1YaGCMoyHf2jKs0TxUFwBsANhtjfP+6FgAY4Xk/AsB8va+FhUUlYtWy8jmOAQryXaV6BQulmXF1B3AjgA0ikuz57AG4/2e+LyK3wB0rPPTEu1tYWIQSjBHk5zl7xhVw4DLGrABwMkHr0pN8bmFhEbIQFBY426qvUnt3rEYE1nVvUcwjunM+X9ysXXoXxo0lb47U2oR+zqlyylol7FTN2fS+ZQrrRJ/HKr+wOd4DdsU62rZ9WCPidZRp+yLlH7Uhl+sKPhDBQW7aQ/3L2n34fANZ48rcqfy4Ivm5SMwNnFB3cGML4n51Fbco3klpWqPU9v8prnMfVW5jxha+Xt3ivyX+2aariJ/fkfMLtxZwLNZo11TiwzaykmF8vkrtZz93GNdZ3DqMj91Q3Se3Ps562xORSoMaxjQ+lTWtrAdUHNivHAfWaLcSc99m+sl2b95qRuIalBkGgMM1LmcPqxYWFpWPQgFynD00OLt3FhYWwUF+4CbBhB24LCwsGG5DLkejUgeuatm5iE/2ru9XJXDs06bhvK7uOJ09mHTcVx6XJsTxSKUVvK1yxpRuE7NaJfTlMJ/Vk7PvL8On3N7XRn0Fb2qzm+OSfmvOsTjao94VwTFj+1Vcko7j2r6pI/HYa1ljS7mCr230/zjv8uBWpWlxSBuQlMN8hbp4car9b4o/yfv3bsoGXd9mdCOe8xJrdJ+9xJoWHmKNbv7eIcSHN2Xh536wzlTQie+tv77sDfqbdBcLcFpP7OHj9w74+90//fA44hNSXiKe1YDvyw0RrGeen8Hf3eJYNlPrO0/dXOqvNtHHy0337bRgBy4LC4uQgwGQF7BVUGEHLgsLC4YBkBuwVVBhBy4LCwuGXSoyfq3eDHckePO4JqffQ9tFySra92jTqDbEO77MGlg4lKal/2toHWdsydvr9GQPrInKPLxHb6/2MbjWF7yz8rtqNIZjcY5exv719VVgUzUcIz49iwOlzurIeZObN53HJ1QeTum/qeJ9e1AyT1KaltawtN+XvrbLeP+lmRwbVXskHzBnuMq/n8Exz/FtWXe6Cf8lPl0Fkp2jPOi1t9odd3k9t7LBPmzzwPrZ07if+M/Kv+tc/Eg8T8UPzo64nrj2vw/nMo3o+7bStDjkD7pSWcN0770VXh4Djh24LCwsQg524LKwsAg52IHLwsIiJOHwgatS6yrGJtYy/07yeoMP3KgW9/vUDpyOhvQHWDf5FhwLNHC3Op72oFdW4gf7sZF6zCwl3CiN7eXmtxF/y6dQQfKsC7mx9mHTetpbTK+/7A3iiwpY2Eh08cXYrYzMdXHO1NyGxHUuoJ+mpessJiveTPGuiqvD+3ng6/a9FFcaWsxfOJdSa1S7C/j3T9/SlHh03F7iPVxcL8C3DqOuWfkkHiLuW1sAAPqkswYlKgd2zRAW/GriKPG4BSonl1MdgXsUj1JcpS76XprE+4GkX0yZ6ipKm0SDJ5MCNwSAGyRgXUURGQDgZQAuANONMc+o7bXhvgot4J5MTTLGvOV3IB/YGZeFhQWjHJeKIuICMBlAX7j/Xa4RkQXGGN/qMmMA/GSMuUJEYgBsFZF3jDHHT3BIALauooWFhUbRwFWaV2BcAGCbMWa7ZyB6F+5CO/qMNT2mpTXgLuhV4tHtjMvCwoJxajOu+iLiu66c5qkzUYSmAHb78D3wC+jAq3A7Ku8DUBPAtcYYFdvEqNSBq/bWoxjYy0eHGq8a7FBcxQYdU/E2A5NZ09Ie8p3vU8mNnP7nh4PDWfNajQuIa7+u0ZhS/H7u8IGqLVvwd+v3PfH/1eb2uu6irvzUFtuIV1M5aVqnaRWxk/i64SxChS/jOocdLue6kOveUaKU9uOqr/hOxXXdxlWKb1Nc3YnXKpFsahofMC9HVVFOZlqzI+tKC9O4BsC39bzFC1M3tqRt8zv1I67rafpB3Vedp6r7TmutuqbFaMW1pqXjG1VI3ifdffy4ni0HPy7gVAauQwE0rhPpbVpY7w/3N9gbQFsAi0XkG2PMEb1jEexS0cLCglG+S8U94JK4zeD/GO5mAB8ZN7bBPYXRafwEO3BZWFgwCgEcK+UrMNYAiBWR1iJSFe6a6QtUm1/hsYEXkYYA2gNQ1jAMq3FZWFgwDKAq853+oYzJF5GxcFeqdQF40xizSURu92yfAuAJADNEZAPcS8v7jDGHTnpQVPbAVR/ASB/eRm1X6XGLeytfoskcP7NpDB+g8zylLXC6mnsF7YOYlRy39UP3s4hrTW0ZLiE+Gl5fc+2X5VLz6LW12R9r6G4OnNrVnD3i2+IX4jtVocLLVb6brgH5X2XQ32cD+2G9sJxjlSJ6cmLniBteU+fnYn/L3xhA3M9jXuHCJ1iP/G6l+jLqMJ2RdTPxvCSuy3hef74XNgw4h7iuS3lnvVeI69gqX3QA630RufxU/mhtrkmZ34UFyejWSpTiEprAXYqzVRv2dmJvsgW4gngXsF46cJb32j6ajvJBOQagGmMWAnzDegasovf7APTT+5UEO+OysLBg2JQfCwuLkIMduBRcAHxmwaviePnkiuPljk6t0NEfVcFT+KVDOO3m+BB+ZD7gC87NeK3fCOKXYBnxy3J5OTY0g5d3rza4pfh9fyyibXppp0tg6UfcLSdyKbSWo5gvbsDLEZ3io5eKswtuIJ6+jVNi3u55NbcHW6+kglOGknMTuMO8ikfMLarc2bNsDf3dC7w07HQPP7bfOLMz8cwVvHTuPZqv/Rj8h3j7erw0XKDWZ80plIjtmGd2uoa2aetmVz6HFEVtZj4ngeMdhj2oirorSSRnDPPI6cxXqxu9tbp3tAxBzXUoxemgEP4hGA6DnXFZWFj4w864LCwsQgp2qWhhYRFysMUyGKYmkOdjFdM1hcsyGZXKIA+rA4xkehRsf6xLM+my91AVufpgCXEdQvBwQ1VK/SOmY5t4rWjejOM661on+Q/uIP50kweIR27mY3/egD11Bqxkfe7J7ux9ou2H33GxxvVOe+Zvq99VX6uvVOjHZRGs92W3r0Zca2Jn37ee+C9oR3xrGtsfo48SVfZwbMz3Waz7dInikADd/z7g8I+vlI+OryaoLXOa57IeFpkChor7HrZRaVojUSIis5jvvYvDH7ahLfEm2E/8S/Qh3rrN48XvCyLKYapUjnFcFQU747KwsGDYpaKFhUXIwaC06TxBgx24LCwsGH+EpaKIvAlgIIBUY0wnz2fRAN4D0ApuQ5NrjDG/n+wYxcfKB8J9UxIeUdt1WqXSpFYlctxXYyU2tNzCsU9QYWD6y4jNYP/iJ1ysaRnVnxSlwcX6uNycrfKLdCzOv1NZ05rTgL3UNrzNKSu6/Naq7vy7P7SbfajXNY8lrnURHfelNaebC9gp92tXD+IzwCk42jr6Bswmfkj53mQqPXLXHk7+/3M8C4ifzbiKeP8uHCd3NlhDO6xyhp4CX+/bfdKzANa4vgb/rlUjOD6wW2vW08LHEcW6ZXzt4x9XotgcplBpOU2H8Qc3vzSD+HFwPOJl+JT4ZJc3MOyg+h5OCyGwVCyNO8QMACoxDfcDWGKMiQWwxMMtLCz+CChfW5sKQcCByxjzNfz+R2AwgJme9zMBXFm+3bKwsAgaisIhSvMKEk5X42pojNkPAMaY/SLSINAOFhYWIYRQ17jKChG5DcBtANBCWyc3V1w5nawZxVbM25Qu03U5x4Gt6amsm5comxttDX1AcSVN/KDKRp3Psgl8w8A6X8Tnat6cY4H0ubs04DivCJV3qTUbnWunQtDw/UjW1C5VcUzaokfnVm51cVyVtoL+Fx4nvgGsyf0b40o8n9Yjw5uxK+9nK1nTGvcg642T0zjB76x6bD3zxFhuv+hVLu+mY/5yfXSjJ5J53yNns21NuI7j4rArxD+oGqiYPK3V/qRssDuov4uYicpu6QG2W5qEe4nXg9e66ji476eFEMhVPF0H1AMi0hgAPD9TT9bQGDPNGJNojEmMiT5ZKwsLC8cgBJaKpztwLQBQZK0wAsD8EtpaWFiEEorCIUrzChJKEw4xB+66w/VFZA/cQQzPAHhfRG6B2y96aEV20sLCopLh8HAIMUZXCqo4JHYQk+QTnmSUT5G2xK21T81FdVzWWYqvVFxJD0blPkqaaq9Ln1+jeITivl+u0jFe7z6c+K1blGCmjpXemnPzoqeyyJA1kifH30Z0I943lS+OznVcj7OJX4v3iGvd5Ca8TfwV3Elc5wLqPFEdG/XqOq5FVzvuN+JXR3xAXGtSOk7uoec4jg3dUTLYmRqrenvj4rpuYa1UyX+AevT0+TCVRzqI80gxVu2v9Ei/cmW6vS4+r26df700gfgTO7waXeIgIGmDOVFJsFJDqiUatE4K3BAANsvaAOXJKgQ2ct7CwoIRAuK8HbgsLCwYIRA5bwcuCwsLf9iBywcFAHy8iCSDN9fKZ01rVXPOz+u6Q2kRukC3LvOuykKJ8vZeM57jvlot2Ek8ZhbH05hLef/sKK/uFPUr+5DfOo+FiV1D2EO9RTrnVUZPV3NzVawpKo2Pv7UJx12NafAq8afA5ccG4WPi2j/rYTxB/BX8nXiaqvuuuS4H1h4/E58Y/w/i0zGKuI5T0+XeLlPl2MASIrY0aUk8bvUubqD+EH3rFWyPa0Tb2kSz/ra3Acfx1FClzV5c8Dfiowr4Rvux37nEdYze4ClfEN+VyPdKaiJ/V9cpffJg6xrF7/MjWGs8LVgjQQsLi5DDH8EdwsLC4gyD1bgsLCxCDoWwRoKEqqB4p72tWTtIUvlxDXUyIVeBh0qf869fp+JlDkSxdtB5C+cXLo3juoz1hh8iHn8vB4ZF3e3VnY60UTFoWSwS7AcnpP0czRrV2aM4TildCYBZhvW4sfe+wXw88/kNWCT7Fhz3peOkJitPfO3ntXDjX4h/3oljmfR3dwm+Iv6W8vP6XDklaU2rptKRtCb2Ym0OujukNLfRXVhT0xrZ4B0+utJc2oSXx99GXF+rh1I4huyifT+gJET05CCywZNZ09IxaN/gYuI6rzQmibVXeCUuVCmvMAa7VLSwsAg5VF5c+mnhdHMVLSwsLIIGO3BZWFiEHCo1V7FTYoT5IKlxMdc6RtN9ymhV5x6q/L5PBrGB1yW5y4hHfcSxT3uH6fp17O+ldSAdm9QD3xCPWe6jNei+qvy27aM4VmgS/kn8XjxPvM1zHEukLOyR0rMZ8YYFrAfmuvhi/U+Z1C4C+1X9Hf9HvA5+V/ww8QLlYd9mH/d3UxMWHH8ExzLVB+uHOo6rY4oy/Fe1CPOU3hn+LvO9o/m71nUihy/90Eu0L5vyx0rvqfJI71VCkkpVxN2KT1Vc1wtVnvTmWebro5Wn/QLWWhcPuqj4/ZjEZPycdLRsuYqSaIBS5irC5ipaWFg4As5/rGgHLgsLCwXnh87bgcvCwkLB+RGolTpwRWYd5xyyfaqBsnbeNJR1kqrKl33gvKXE5w/h2KXBXTlepukO1tCua/0k8TvxCnGdjxczT8XPJPi8v5E3aR3jPVxH/D8pHIeUEsuaFThMCmu6cBzXv8FxSt1c3xLXtffOAnu071BBcT3ncu1AP91G/T5XNuD6fY814SKZum5jotJMtL75uDpBn1j2+0pAsjp+W94+mrc3ncrf9fDeHxLP6u59LvVjRAJtq6KCmLrO4BzZLZM4L7KqMvtq8wLrfduXsb55XIm1cbGcVynKnyvepYzl1PY+6V4vtlrlMt7YGZeFhUXIwfkDlw2HsLCwUDBwi/OleQWGiAwQka0isk1E/IpHi8g/RSTZ89ooIgUiUmJpHTtwWVhYKJRfKWsRcQGYDODPADoAGCYilN9ljHneGJNgjEkAMAHAcmOMLkLNx63MOK6OiZFmdpI3WVFrCR07q9gdlWu46THWvA4oT6kmSjSLu521g1en3EJ8VBbn90WqWCG/+t1K99ky2qt1xG1U/k+qtl76UI4FqpvOsUBylzqX8rBXEg+yPuL/OT9FcK5fAxWcpOscVlMe8T+Dcyf7Pq4M/jmUyM8HfdOn/N3o3EMdB9Z3Bx9f562uVoFrW1X/Jjz3EnfgV9U/5cXmF5bU1ec9p4ki/S7+rvS103pcPXDxAq2NTp07jvj8oazFtsM24h2vUn8HyUy1NZmvzXLi9UDST2X0nJeOxi+47KSILzGOS0QuBPCoMaa/h08AAGPM0ydpPxvAV8aY10s6q51xWVhYKJTfjAtAU4Cii/d4PvODiFQHMADAhyfa7gsrzltYWCickjhfX0R857PTjDHTfPiJZn8nW+ZdAWBloGUiYAcuCwsLP5xSHNehACk/ewA09+HN4B8IVYTrUMo1aqUOXLmIoPieKzI+4wYqTApxTLXXt9a4tK6ybwonnY1dyZrWb91rE/8qqhfxYflcoPvIKPbcivuihJi0HUyjE1V+m9YptL6mNSVVcteVz3mYnSeyt5if/z5b0mN7LMcWaR90P6jD4zGmHRewLtNxodJp1K29/lbm5wzlf7JH32cPrAlbXuIdopi++eow4n9N4vt/1QOqfsFun9gs9fxKx1k9CI73q66epr2Ha4mPgipuMIXp4OXKj+sqlIz7FP9ZcV9JrRDlgHJN+VkDIFZEWgPYC/fgdL1uJCK14Y5eHK63nQh2xmVhYaFQfnFcxph8ERkLd5ldF4A3jTGbROR2z/aiYX0IgC+MMfpf+AlhBy4LC4sToPxSfowxC6HWGD4DVhGfAWBGaY9pBy4LCwsF50fOV2ocV2I7MUm+tlPK6xtXKK4mjR+N+jPxDir/TscuPaHibS7AauKtsbOE3gILVDDQxLlce5B8m8apnZXmtV7xc0Yy3/tWyXFM2q+q80VKdFLpbHg2AGcLd7x+D0sLV6g6jI22sAf+mjjOnew8mfvz+hg+3q1LVeCX9sBS/0K3DOV8QO0pf9Fu5fOuws7wsuJaJ/KRN1f0Po826TxKjfZZfLG/ibqIeN/bVWeuUQfQvvBar+tZsl6n/252PeutpXBF4u9Yn5RXxjiutgaYWMrW11k/LgsLCyfAukNYWFiEHKyRoIWFRcjB+RpX5Q5c1QD4SiP6waeSQbQH1FXLOe5rfk/O+Yqby/mCdwz9D/HY5XuIa9/2j5XINnGH0rR0vptPPuIKte0i7hrasUTkF+eVrfLhrktjE/XjT3DM2QrlcX+RklX2dmfNrGlPjpN6vTVrUK2U3qfjuhpF8y/gFzfGkheuLXiP+NO9xxHXtQLPS+XkTl2nUef/XfS+0rh03Jvy/NdxZL6hVlV6c85sfKoSDLkkpJ9e6JfXqXJsZ/XmmpQ6p1Zf67/OYE3rtZEjeP9ErnnpW/cxE0tQdjh/qVimXMVAdhUWFhahiKIZV2lewcFpz7h87Cr6wh3Wv0ZEFhhjfip5TwsLC2fD+TOusiwVLwCwzRizHQBE5F0AgwHYgcvCIqThfHH+tOO4RORqAAOMMaM8/EYAXYwxY1W72wDc5qGd4J/15iTUB1TBP2fB6f0DnN/HP3r/WhpjYgI3OzlE5HNPP0qDQ8aYAWU53+mgLDOuUtlVeCwupgGAiCQFI1ittLD9Kzuc3kfbv8AIxkB0qiiLOH8qdhUWFhYW5YayDFzFdhUiUhVuu4oF5dMtCwsLi5PjtJeKJ7OrCLDbtADbgw3bv7LD6X20/fsDoFKTrC0sLCzKA7ZYhoWFRcjBDlwWFhYhh0oZuJyYGiQib4pIqohs9PksWkQWi0iK52fdIPavuYh8JSKbRWSTiLvyolP6KCKRIvK9iKzz9O8xJ/XPp58uEflRRD5xaP92isgGTxXnJCf20Ymo8IGrNJVsg4QZcNdw88X9AJYYY2IBLPHwYCEfwD3GmLPgLl86xnPdnNLHXAC9jTHxABIADBCRrg7qXxHuApfndVr/AOASTyXnovgtJ/bRWTDGVOgLwIUAFvnwCQAmVPR5S9m3VgA2+vCtABp73jcGsDXYffTp23y480Id10cA1QH8AKCLk/oHd2zhEgC9AXzixO8YwE4A9dVnjuqjE1+VsVQsdSVbB6ChMWY/AHh+anOUoEBEWgE4F8BqOKiPnmVYMoBUAIuNMY7qH4CXAIwHF+1yUv8Ad7bJFyKy1pMeBzivj45DZfhxnUolWwsFEakBd0nyccaYIyJlshMvVxhjCgAkiEgdAPNEpFOAXSoNIjIQQKoxZq2I9Apyd0pCd2PMPhFpAGCxiGwJdodCAZUx4wql1KADItIYADw/U4PZGREJh3vQescY85HnY0f1EQCMMYcBLINbM3RK/7oDGCQiOwG8C6C3iMxyUP8AAMaYfZ6fqQDmwe264qg+OhGVMXCFUmrQAgBFdpMj4NaVggJxT63eALDZGPOizyZH9FFEYjwzLYhINQB9AGxxSv+MMROMMc2MMa3gvueWGmOGO6V/ACAiUSJSs+g9gH5wu6c4po+ORSUJkJfBXTj8FwAPBlvY8/RpDoD9cNs47gFwC9zFzJfAXexrCYDoIPbvIriX1OsBJHtelzmljwDOAfCjp38bATzs+dwR/VN97QWvOO+Y/gFoA2Cd57Wp6G/DSX106sum/FhYWIQcbOS8hYVFyMEOXBYWFiEHO3BZWFiEHOzAZWFhEXKwA5eFhUXIwQ5cFhYWIQc7cFlYWIQc/h+TJByvfuSgWAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] @@ -319,24 +326,138 @@ "\n", "\n", "# Find height levels for contours at different standard deviations\n", - "print(optimum)\n", - "sigma = 1\n", - "xnorm = np.array([0, 0.5, 1, 1.5, 2]) # standard deviations for heights\n", - "heights = np.exp( - (xnorm)**2 / (2 * sigma**2) ) \n", + "#print(optimum)\n", + "#sigma = 1\n", + "#xnorm = np.array([0, 0.5, 1, 1.5, 2]) # standard deviations for heights\n", + "#heights = np.exp( - (xnorm)**2 / (2 * sigma**2) ) \n", + "heights = np.array([0.135, 0.325, 0.606, 0.882, 1])\n", "print(heights)\n", "\n", "os = optimum[-1]\n", "amp = optimum[0]\n", "contour_levels = os + heights*amp \n", + "contour_levels= np.sort(contour_levels)\n", "print(contour_levels)\n", "\n", "# Final plot\n", "#fig, ax = plt.subplots(1, 1)\n", "plt.imshow(img, origin='lower', cmap=plt.cm.jet, extent=(x.min(), x.max(), y.min(), y.max())) #origin='bottom' didn't work\n", "plt.colorbar()\n", - "plt.contour(x, y, img_fitted, 5, colors='w', levels=contour_levels)\n", + "plt.contour(x, y, img_fitted, colors='w', levels=contour_levels)\n", "plt.show()" ] + }, + { + "cell_type": "code", + "execution_count": 83, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[1. 0.8824969 0.60653066 0.32465247 0.13533528]\n" + ] + } + ], + "source": [ + " sigma = 1\n", + " xnorm = np.array([0, 0.5, 1, 1.5, 2]) # standard deviations for heights\n", + " heights = np.exp( - (xnorm)**2 / (2 * sigma**2) ) \n", + " print(heights)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [ + "[1.,0.88, 0.606, 0.325, 0.135]" + ] + }, + { + "cell_type": "code", + "execution_count": 90, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "array([10.135, 10.325, 10.606, 10.88 , 11. ])" + ] + }, + "execution_count": 90, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [] + }, + { + "cell_type": "code", + "execution_count": 100, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[[1.0200708 1.064459 1.001646 1.0051079 0.937599 1.0175885\n", + " 0.96301293 1.0054469 0.96290076]\n", + " [1.05853 1.0182664 1.0239749 0.93696785 0.9614606 0.9738683\n", + " 0.9017772 0.9327952 0.9521524 ]\n", + " [0.98196447 1.0818717 1.0192714 0.96637696 0.9451842 0.9370096\n", + " 0.9464692 1.0042784 0.9692586 ]\n", + " [0.95845973 1.0763768 1.0382882 0.9582531 0.9462724 0.96113765\n", + " 1.0498253 0.9213169 0.9789404 ]\n", + " [1.0155486 1.0051879 0.9592701 0.96256775 1.0141811 0.93864226\n", + " 0.9912149 0.98979616 0.8719947 ]\n", + " [0.96749 0.9859212 1.015961 1.0136578 0.94090223 0.9777953\n", + " 1.018927 0.9675724 0.94836783]\n", + " [0.96860135 0.9384735 0.952281 1.0193214 1.0817233 1.0644275\n", + " 0.9791999 0.966643 0.99475867]\n", + " [0.99147207 0.9844064 0.9521766 1.0357393 1.0254681 1.0550414\n", + " 1.0165464 1.0327921 0.9531771 ]\n", + " [0.9835919 0.9764188 0.98863244 1.0181607 0.99789697 1.0286107\n", + " 1.0462499 0.9830032 0.9930277 ]\n", + " [1.0419786 1.0235844 0.95347303 0.9305422 1.0051019 0.9533537\n", + " 1.0270752 1.0441222 1.0376717 ]]\n" + ] + }, + { + "ename": "IndexError", + "evalue": "shape mismatch: indexing arrays could not be broadcast together with shapes (10,) (9,) ", + "output_type": "error", + "traceback": [ + "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m", + "\u001b[1;31mIndexError\u001b[0m Traceback (most recent call last)", + "\u001b[1;32m\u001b[0m in \u001b[0;36m\u001b[1;34m\u001b[0m\n\u001b[0;32m 1\u001b[0m \u001b[0mprint\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mimg\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;36m10\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m1\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;36m10\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 2\u001b[0m \u001b[0mind\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m(\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m1\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m2\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m3\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m4\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m5\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m6\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m7\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m8\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m9\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;36m1\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m2\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m3\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m4\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m5\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m6\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m7\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m8\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m9\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 3\u001b[1;33m \u001b[0mprint\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mimg\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mind\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m", + "\u001b[1;31mIndexError\u001b[0m: shape mismatch: indexing arrays could not be broadcast together with shapes (10,) (9,) " + ] + } + ], + "source": [ + "print(img[0:10,1:10])\n", + "\n", + "print(img[ind])" + ] + }, + { + "cell_type": "code", + "execution_count": 102, + "metadata": {}, + "outputs": [ + { + "ename": "SyntaxError", + "evalue": "invalid syntax (, line 1)", + "output_type": "error", + "traceback": [ + "\u001b[1;36m File \u001b[1;32m\"\"\u001b[1;36m, line \u001b[1;32m1\u001b[0m\n\u001b[1;33m [0:10]\u001b[0m\n\u001b[1;37m ^\u001b[0m\n\u001b[1;31mSyntaxError\u001b[0m\u001b[1;31m:\u001b[0m invalid syntax\n" + ] + } + ], + "source": [] } ], "metadata": {