{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# DataArray - Arithmetic\n", "\n", "We can basic arithmetic operations (plus, minus, multiplication and division) with DataArrays, both with scalars, numpy arrays and other DataArrays. The output will in all cases be a new DataArray." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import mikeio" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "fn1 = \"../tests/testdata/oresundHD_run1.dfsu\"\n", "da1 = mikeio.read(fn1, items=\"Surface elevation\")[0]\n", "fn2 = \"../tests/testdata/oresundHD_run2.dfsu\"\n", "da2 = mikeio.read(fn2, items=\"Surface elevation\")[0]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Multiply or add scalar\n", "\n", "We can scale a DataArray or add a constant using *, +, / and -" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.18681715078618902" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "da1.values.mean()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(1.186817150786189, -0.813182849213811)" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "da1_A = da1 + 1\n", "da1_B = da1 - 1\n", "da1_A.values.mean(), da1_B.values.mean()" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(0.37363430157237804, 0.09340857539309451)" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "da1_C = da1 * 2\n", "da1_D = da1 / 2\n", "da1_C.values.mean(), da1_D.values.mean()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Difference between two DataArrays\n", "\n", "Assume that we have two calibration runs and we wan't to find the difference..." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAOEAAAEICAYAAACpj5OEAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABpOUlEQVR4nO2dd5wkVbn3v09Vh0k7s3kXdmEXlgySRJIKqFdRDHhfs15FBLmoF1HEHPGa4zUSFBUVIxhQESMYl5wkSFyW3WVzmBy6q573j3NO1amenpne2Zmd3tn+fbZ2uqvqVJ3qPk8/+XlEVWmggQamDsFUT6CBBnZ3NIiwgQamGA0ibKCBKUaDCBtoYIrRIMIGGphiNIiwgQamGA0i3AGIyCUi8kHv/ZtEZL2I9IjIHBF5qog8ZN+/eAqn2kAdQxp+wpEhIo8BC4AyEAH3Ad8FLlPVuOLcPNAFHK+qd9l9fwKuUdUv7cx5N7BrocEJx8YLVXUGsAT4FPBu4PIq5y0AmoB7vX1LKt7XDBHJjWdcA7seGkRYI1S1U1WvAV4BnCEih4nId0TkYyJyAPCAPXWbiPxZRB4B9gV+ZcXRooh0iMjlIrJWRNbYsSGAiLxeRP4hIl8UkS3AR+yYz4nI41bMvUREmu35p4jIahF5h4hssNc8081XRJpF5PMislJEOkXk797Y40XknyKyTUTuEpFTdtoH2cAwNIhwO6GqNwOrgad7+x4EDrVvZ6rqM1V1GfA4hpO2qeogcAVGtN0POAp4DnC2d/njgEeB+cDHgU8DBwBH2jGLgA955y8EOuz+s4Cvicgse+xzwJOBE4HZwLuAWEQWAb8BPmb3XwhcLSLzduiDaWDcqFsiFJFv2V/4e2o8/+Uicp+I3CsiP5jk6T2BWcA1Q0QWAM8D3qaqvaq6Afgi8Er/uqr6FVUtAwPAG4G3q+oWVe0GPlFxfgn4qKqWVPVaoAc4UEQC4A3A+aq6RlUjVf2n/SH4L+BaVb1WVWNV/QNwK3DaOD6HBiYA9ax3fAf4KsYQMipEZH/gvcBTVXWriMyf5LktArZs55glQB5YKyJuXwCs8s7xX88DWoDbvPMFCL1zNluCdegD2oC5GP30kRHm8TIReaG3Lw9cvz0P08DEoW6JUFX/KiJL/X0isgz4GmaB9gFvVNV/YzjG11R1qx27YbLmJSJPwRDh3zHiY61YBQwCcysIx4dvqt4E9AOHquqa7ZzmJgwnXQbcVWUe31PVN27nNRuYJNStODoCLgPOU9UnY3SZr9v9BwAHWMPGjSLy3Im+sYi0i8gLgB8B31fVf23PeFVdC/we+Ly9ViAiy0Tk5BHOj4FvAF90nF1EFonIqTXcKwa+BXxBRPYUkVBEThCRIvB94IUicqrd32SNPIu353kamDjsMkQoIm0YI8NPReRO4FJgD3s4B+wPnAK8CvimiMycoFv/SkS6MRzk/cAXgDNHHzIiXgcUMP7GrcBVpM9QDe8GHgZuFJEu4I/AgTXe60LgX8AtGNH500CgqquA04H3ARsxz/VOdqG1MN1Q1856K47+WlUPE5F24AFVHbZoReQS4EZV/Y59/yfgPap6y86cbwMNjAe7zK+fqnYBK0TkZQBicIQ9/AvgGXb/XIx4+uhUzLOBBrYXdUuEIvJDYDnG5L5aRM4CXgOcJSJ3YSJRTren/w7YLCL3Yax871TVzVMx7wYa2F7UtTjaQAO7A+qWEzbQwO6CuvQTzp07V5cuXTrV09jp6OnrBYG25tapnsqU4bbbbtukqkkI3anPaNXNW6Lh5909+DtVnXBX1FSgLolw6dKl3HrrrVM9jZ0CVeWy63/ELX9eTrhxkLgY8JXPfpGmYtNUT21KICIr/febtkTc9LvhLsz8Ho/M3WmTmmQ0xNEpQhzHXPK7Kzn3fedz61U30DUz5tHTWolC5au/uGKqp1dHUCKNh23TCXXJCaczyuUyl//+B9x8/XKkP6Zz/zzbTmqhVBDKccj6Q5sI/nYH257Xycz2jqme7pRDgZjpbTxsEOFOwuDgID/4w9Usv/7vxKWIrgNzbNuvSFwQVAVViFXZsleBOQ8O8t5vfpaLL/jYVE+7LhAzvThfJRpEOMno7e/je7/9Mbf99SZilK6Dc3QtayLOGeIDUAUhIBdAKQ5YdXQL+/1lIx+7+lI+8JL/nuInmFooSmmaiZ+VqIkIbRzmN4HDMBLCG1R1uXf8NZg4RzA5bW9S1btEZC9MKtJCIMbUZtkt6q10dnfxnd9cyT3/uJO4AN1PytG9b0AkAYoAggig5gNVlEgDhJiB2SErj2qBP97Oj/f+Pa94ynOm+GmmDgpEDXEUgC8B16nqS0WkgMlz87ECONnm8j0Pk+1wHCaL/B2qeruIzMDkxv1BVe+bqAeoN2zatpnLf/ldHr7l38QtQteTc/QtDdDA5gTaH3UTJCFECoIQq0kWRCCOha17NdPSGfP77/+cE5Ycxt7z95yiJ5p67PY6oQ2cPgl4PYCqDgFD/jmq+k/v7Y3AYrt/LbDWvu4WkfsxuXjTjgjXblrHl666mC13r6PcAT0nhPTvJcQERCqIAmK4n6CU1RimBSFCUcKEQwZiFt4ThzbR3FXmw1/8FJde9DkKhcKUPuNUQIFomkd11eKi2BeT8vJtEblDRL4pIqN5k88Cflu502ZEHAXcVG2QiJwjIreKyK0bN26sYVr1g3/ecxMf/uCH2bhiHdtODth8WsjgEiN2xoa/EWOMLwapLhiree+I0+0TUSQIWHFMG9pf4vM//MZUPFpdIK6yTSfUQoQ54GjgYlU9CugF3lPtRBF5BoYI312xvw24GlNfpavaWFW9TFWPUdVj5s3btWoOHbr0INqWziTshNZ/xXQsj2l+OIJBs1xceQrH6QACUUuAASCoKuqI0W6qSjkfsvLoNlYsv5tf3nrDVDzelEJRoirbdEItRLgaWK2qjoNdhSHKDETkcIzx5nQ/g0FMUdyrgStV9Wc7PuX6Q0dbB59/16d46zvfypKDlzKr2MGMu5U9fh4x89YyQY9ZNKqG00UxqIrjhyhGP8Tnlgpqv56uOXnW7tPMr77/U7Z2bdvpzzeVUIVSlW06YUydUFXXicgqETlQVR8AnkWFTiciewM/A15ry/+5/YIplHu/qn5hYqdeXxARDlt2KIctM5UP4zjmmuXXce11v2XGrwbpPCjHtkNyxHlzrqqCBIYQ1RVyUjs2ILZEGitECI8f1MaMrdu44Asf4/IPfYYg2F2CnYQIGfu0XRi1fpPnAVeKyN2YGpifEJFzReRce/xDwBzg6yJyp4i4wM+nAq8Fnmn33ykiu0VpvSAIePFTT+Oy//0KLznzlbStilj8mwFaV5ZBYwKBgBjD8jQjYMVquGIUY3XKAELh/iPbkU09fOFH356ah5oCKOaHqHKbTqjJRaGqdwLHVOy+xDt+Ntkitm7/32Ga/4zVgFOPfQbPOOqpXPSjLxDcuILmjTFbjskTBoLGMa6KoYgSa2A5JSDmw9PYujNaQx48Yia5v97Mwyf9B/stXjKFT7Xz0OCEDUwICvkCH3/te3jL28+jbWXErNtLoJoYYQJxWqGx/xmjqbWsGmYJQOeCIp1zi3z9p9+boifZuTDOehm2TSc0iHAn48j9DuP/nf0qZjwaMevusokb9RaVoTsjnqpCHENitMEYdR5b1kbPA6tYfu+dU/AEOxcKlDQYtk0nTK+n2UXwvCNO4cxzz6L9gTLtD5UAQ1xxLKgGqIbEsRCroASo5YRRZIixa2aBNXu1ctnll5tE4GkMRYgIhm3TCdPraXYhnHjYsZz6yhcx/84SLesijAJoCK8Upz5FyRhOxRKksOKADuIo4hu/nZZenwxilWHbdEKDCKcQLz3p+cw5bhmLb+yj0F2mFAmRpg59jVPuF6lYQ6o5FgWwds8W7rnl9ql7gJ2Ahk7YwKTjE2e8k2DvDva9vo+W9RFRbCJoQBJXRZpvaL4utfrhhgUt5Lb20dPTM4VPMNkQIg2GbdMJ0+tpdkGICJe841MccsLR7L+8i1mrBtHYGWRMhI0TSU1sjXHiiwiDTebrW791extE7TowmfXBsK0WiMhzReQBEXlYRIaFWorIQSKyXEQGReTC7Rk7kWgQYR0gCALe8epzOOnFz2PfO7qZ99gAqgGxmiDwKDZEl7owDK+MQqEcCreveGDU6+/KUBWGNBy2jQUxHZC/hukJeQjwKhE5pOK0LcBbMQ1Vt3fshKFBhHWE1z73xSx7/onsc183+/6rk7CUDQA3kWppBgaBsHF+E39f/s/qF5wmiG02ir/VgGOBh1X1UZt+9yPSiu2AaaFn+5WUtnfsRKJBhHWG97zwDM4471xmb4t5yvUbmbOuP3NcRBFxCT3K+kUt9K14glKpch1NDxjDTFUXxVyX+ma3cyqGLiLbdHW13VcLdmTsdqNRY6YO8bRDjmLB+bP4zMc+iUTq5SGqiazxGMFQS0CgMDQ0RD6fn4rpTjJkJEPMJlWtDKXMDhyOWqNOd2TsdqPBCesQA4MDfOyrX2DrgiKb9iwSBDFhEBOI4pInXLhbnDPrZc2mSWtOPKXYAcPMamAv7/1i4Ikab7sjY7cbDSKsQ1z41U+iQ2UeOWIGQWAdFmJEUZdvEdj3pWJA18w8V/zxmqmd9CQiUhm21YBbgP1FZB9bF+mVQK0f0o6M3W40xNE6w/rNGyk9uI6HT5pJlE8Tf8GKoarYgjWgxl+4bnELTXf9m6GhoWlXh0YRSrr9y1RVyyLyP5i2eSHwLVW916XfqeolIrIQuBVoB2IReRtwiKp2VRs7MU80HA0irDN85fc/YqA1YGhWQE5j65gXYuemEJdRYQolCrBxURNLH+zmi9f+kHe/+IypfYAJhjPMjGus6rXAtRX7/BS8ddiiZLWMnSw0iLCOMDg0yNqb72fz0iZExCw9MQJoYGvOKCZ8zdVHDABy0DszR8/WbVM298mCUrP4ucuiQYR1gjiOefelnyQmZuN+RqQUSU1yziIqACFEkdtvDohCfzw93RS1RsjsqmgQYR2gb6CfC778YUqrulj1tBYIQTxRtBrCACIF1RgwIWxB71DVc3dlqDLtYkUrMb2fbhfAzQ/cwds+9E4GNnax8lktDM4JyYdx4oIwAVpeCTYbtmZrCROGShDE9M7Ks/W+lbzlk+/n2zf8ckqeZTJgDDPhsG06ocEJpxC33H8n3/jKJXQtzbH+iFbiXJC4I/LECReUyCT3OvFUNSZ2NfPtn41LigzMyDFndQ9DP76Wu26+nXed/Rb2nD1/ip5u4jDdkngrMb2fro5x/+qHuOziS+jaP8fGYwo4K3zgEne9AI0gMBzQOSxENHFdSHIOdM/JsfKIGdx7yiw2btnEBz/z8V0+814ZntDbSOptYEJw6ZXfon+OsPmIvM2QMBpgnLQ/CQg8fTAQUCKwyUyhmC8vQ6xi3BYDM0IeOKGDoajEeZ/9EHG8axeOb5S3aGDC8Zd7bqRvxRa2HFlAXX80TCmnEFfSyZGi437Y6Bn1jqYloARTuS2058R5ePj4dtjcy0Xf+OJOfsKJg6k7GgzbphOm19PsIvjxz35Kz9KQUkfamSlN3IVAhMCGqDnNMBE7xeUUKmFSJpGkoYygtlyNMNQS8vAJHay940HueHDSAj4mGcNLWzTKWzSwQ7jmpj9QXtvD1sNcxoOkvvfkf7VEmYqRTux0VtGcGJ7oRFAf4i3Svpk5uufkuXb59ZPzQJMMU/KwYR1tYALxuxt+T8/SkHJbtv+ErbVN0hRGlDAIiIwb0LbUViBIOvsmDnz1gmjsNcWeDdC5IM8jDzw06c82GVCVaSd+VmJ6P12dYXBwkNLKLnqXuF/ylABd0V8HV2s0mzVRwRFJM+6tEyMd72mOA20hQe+uG03TKPTUwIThVzf+njiE/nnuY3dlDNNmoWCd8JJtEy2+qJpx3qcl8mOv/kyQXE0Zag4JByLK5fLkPdwkwTlnxlHeYpdBQxzdifjbP/5O79IQDVNTS7qgfPHUiKRBYoKxldY0jScVJel3X5nynQq2BrH9lkulErncrvaVj5hZP20wvZ+ujvDImsfof3wbPfvmElEzbfGVckBIg9Nia+nMipmeIkjWJONI29+nQByYMQMDAxP6TDsDxkUxvZ31u9rP4i6Lb179XQYWhgzMCpNS9pUEpMmW5g0KEGkawuafL7YOKRid0e/b5xrKKJIQ4WBp1wvwdrGj0xkNTrgTsGLt42y+bw2bD80lTS6NndPX3FKicQhFiNWvOSqpwUYl2Q/mmj6ZKmlFqDg0X/MTWzdN9qNOCsZb/HdXQYMT7gRc/rPv0T8/YGiuiYdxXC9OXAlZnuiOl2NDqmAIsuw6hybOCiHCVWOrNNe4q5liUF1z8/zyL7/j6AMPncxHnXCYVKb6FT9F5Ms1nNalqh8Y6WCDCCcZQ0NDbPjX42w7pZhJ0A0QynE29tORmG8TdXvLmq0z4wLd0qg3a9Dxrqbe4l23TzOtdzxIT18vbS2tE/+gk4g61wFPx7SLHw3vARpEOFXo6e1BFIbaA6JYrKElSIgqIEhcEeqc7gIkC8/qfEgSP+MILVJX+iIgTq6X6pYJFLYuLFDOC9/73S9403++ZrIed8JhsijqWvz8oqpeMdoJIjJrtON1/XTTAeu2bgSgJKHR01Qy/j/3OtK021Js3Q++ludMM+rJnIoQEKCa5aY+AQY2wgYR1u3bwq1/+ht/u/vWiX/QSYIJW6vfTr2q+n8iEorI20c7Z7Rr1M/TTFNc8csf0TcvRPMeQeD8fbYFdmKQCYwIWSF+OSupc+hLYnTxU5+y5zs90UsNZv2+Laxf0sx3L/4G19x8w6Q988RC6j6LQlUjdqBXRU1PIyIzReQqEfm3iNwvIidUHH+NiNxtt3+KyBHesZ3WYqre8Jd7b2brA0+w/qgiagOuEytowvlcHRlJHPJmoUnGkhqrZHQ89QnS7cO3ikrii4zt+TGw6uA21u3TwjU/uZqhoV3DZbGLRMz8Q0S+KiJPF5Gj3VbLwFp/Ur4EXKeqBwFHAPdXHF8BnKyqhwP/C1wG7PQWU/WG+x76N4OzAoZmhomYGatYR7y/kJyX3hFlagXFEp/jolknvyY6pAjDlmaaBJU9snr/VsqlMh/+0cUT/MQTD2cdHUcF7lr6E4qIfNkev9snGhF5u4jcKyL3iMgPRaRpjNudCBwKfBT4vN0+N+oIizENMyLSDpwEvB7AtorK/ISqqt+b60bSgqpJiyl7Lddi6r5aJrer44lN6xlqCVACUxcmIbbAc7NbZBaXpzNajldOQraFoL9Mc3fEvEcH2LhPE93z84lPUFE0TrnF8CtCnBdWH9BKeOsDDLxigKbiWOtrajEe8dNjAM/G9Ja4RUSuUVV/7T0P2N9uxwEXA8eJyCJM38JDVLVfRH6CKYX/nZHup6rP2O5JWtTydPsCG4Fvi8gdIvJNERnNxn0W8Fv7uuYWUyJyjmtztXHjxhqmVf9o72in2BmbsoSSFRuT/Hk1XXkjP9rFEaO6P5IYZCKN2f9vXex7Uw8da0vs989uDvn9NjpWD9rRzk2R/u84ZxL+prBxcRORKF/85ZUT/NQTix2oMVNLj8HTge+qwY3ATBHZwx7LAc0ikgNaGKMhjIgsEJHLReS39v0hInJWLROthQhzwNHAxap6FNCL8XtUm8gzMET4brerymmVdgSzU/UyVT1GVY+ZN29eDdOqf5x96iso9Ma0bIyMbqdp/3mAODbGGPG/huSYb4gxUTMohP1KU2/MBz/wAS699FI+9slPsOSog1h6Ww9Lbu1GItO5KbDXdPcyOmFKmBoErFvawiO33J2xrtYbFChrMGxjYvoTVj1HVddgRMnHgbVAp6r+foypfgfTu2JP+/5B4G21PGMtRLgaWK2qN9n3V2GIMgMRORz4JnC6qm72xu60FlP1hlntM2GfDmY9OGS5myShZ8bgElCOXf68kIatSeIfc8YZsLpRGFAqCD/6428AmDd7Dhe97nzOesdb6NimHPDPTqSsid/Qhx/ypgobFjcRdg3wyOqVO+cDGSdGsI5ucj/adrusYlgtDKDqOdavdzqwD4aoWkXkv8aY5lxV/Qn261TVMhCNMQaogQht04xVInKg3fUsKnQ6Edkb+BnwWlV90Du0U1tM1SPe/vI3MuOJMsXOtGgvCHEcWFdEyhl9OLeEs5oaTqpEuYCHn9zOQ/+4nR9e/+vk/OP3O5xPf/CjFPuUOY8OZERYc5VsrCoI5YIJjO4Z6Jvgp55AVBFFaxRHa2EAI53zH8AKVd2oqiXM2j5xjPv1isgccG5ZOR7orGWitWq85wFXisjdwJHAJ0TkXNdmChO2Mwf4uojcKSK3QvJr4FpM3Q/8ZDJbTNUjDl66P6WFTcx8fCjhRM7p7ggSBI2HE42DIpRjscHL0Dm3wGOHzeDPP/kVP73pT8l5c2fOZukxhzJ/1YDnK3TXyIbICWLcJgJb+3om49EnBMq4XRS1MIBrgNdZK+nxGLFzLUYMPV5EWsQ0+3gWwz0ClbjAXm+ZiPwD+C7GuDMmagpbU9U7gcrWxH6LqbOBs0cYu9NaTNUrDjvwEG6/7+40wNq6IhLLpRU5AwmtfuabVkhEWT/vcMPSZpr6In77o5/xvCcdn8SDdm7aSm+7+1qFONYkvtSFvRlOYq4eCwT1HZs5rtjRWvoTYtblacDDQB9wpj12k4hcBdwOlIE7sG63UXAvcDJwIOaLeoAamVx9hR5MUxx2wMG0bClT6IkzBBVrqh+a90qkUFafAM2xYT5AhccPaKUUKO+49NOACRbf9sgats5vwoi5av/6zn+foIUoF/DIhjWT/yGMEy5QYTxJvap6raoeoKrLVPXjdt8lrkehtYq+xR5/kqre6o39sKoepKqHqeprVXVwpPtYLFfVsqreq6r3WDF2eS3zbBDhTsBznvx0ZixdwN439yJx6ooHLEdMfXoCic8w0uFB2xmEwiNHdBA/sJ5LfvEj3n3xZ1Fgy4KiPcF3ixitsNIWsXlBE7fcfPOEPOdkwIjiwbCtXiAiC0XkyRh3xlFetMwpGNfGmGhkUewEiAgfftOFvOM976JtQ5meBUWvoqjH9UjD2RxhZrLlVQnEdO115/fMyrHqgDb099fT257niSe1E4eSXLRaDRoT/G1+CtYvambhTZsZGBigqak+nfZ1GqbmcComkGUx8AVvfxfwvlou0CDCnYSOGe2EM5vJ98ZE1jqTBlkbuPfgE6bPKTWJJU2JS3liv1bWLGvNBASQjPcNNGpSp0QRW16jr6NALHDTv//FyUc+ZRKefAeh49MJdxZsGtMVIvISVb16PNeoH76+O2BGkeZtURJg7ecKGl0xyOxPHO24v5LZ58ZZm/iw2wXJ+dmUKLxdUQC97XmuvfvGCXrIicWO6IQ7Gf+YzIiZBiYIrzj1dGY/PkihJ/LC0IbnDsJwx3ocQ2Wib5rElDXcuPC0su1xL0jiI0wqvcUpkZcKIe1x/RZT2kWI8NtMYsRMAxOEZx51AtEebSy83zjH04K/2WDrKHYRNdkYUjJnGQQV4qsgVTIMbRZ+FeuOIJTyworN68bzSJMORYjiYNhWh5i8iJkGJhZnnP5yZq0ZQvoV9autKZRjQ4CV6UdJlI1z9mtgU5wYRnB+doYM45gGgT0SYMLihoohuf76rc69i+QTTnrETAMThJOOOJZye4F5j2cL8VZL0q0Ow+n8aJvK60RJbKo7I3tVN94R50BzyGBnfUbNqO4y4mi1iJnzahnYsI7uZIgIBx51GKU77uaJ/U2Ui58tMVI+gzGoVnE3VJzjB6bhrqfZoopU6KClfICUhktOG7dt4X+/+xXCklJsb+UdLzuLeTNnj/WIEw6tT6LLQFVvF5FMxIx12I+JBhFOAV76tFNZccPt5AZiSk2pQUTdf04VTOI/A1ttTZP9JpPe+RL9MLfhnNFcwSuvgWlEmowLJJPQuG7zRi7+3U9Ys/we+lsD+mblmPHgBj74uY/x+fd/jNbmmnzQE4S65XwZ2CTi04ClGLp6joigql8YdSANIpwSLFu8hFJbnpnrhtiwpDnZn7gfMi3mjcYwzNiiAZEXgaqZEdWNM77xxs8hHCoG5EoRD6x4hI9d/lVaNvYx0Bqy7qAWNi4tEIiQL5tk4gv+76Nc+t5PjfPJx4ddgRMCvwIGgH9BJhZjTDSIcAogIiw+cF96Vq1goyVCx8VS2kgXniEqE3fqRNa0T2GW2ATDGStJMNmjw0Xe3hl5YoHPf/ozDMwv8uDJsxlqD0DS7H/Jw4pj2jj4z1t5ePVj7Ld46Y59CDVCFVOvtf6x2NZY2m40DDNThKX77ktLVzkRD1O/3XBjS4wm+l5ccbza8qxW9MnBiaL+HglgzT6tPHJoO/cfPZO+thAVL55VABX623N0z83xye9/nTjerh/7HcIuYh39rYg8ZzwDG0Q4Rdh37p7kSpq4Ggw8f6G66my+zzDrT6xmyFFvk4r9DsOy7gNYfWAH6/dutcRnXSDYDsEYh5cirDqsDX2ii7d95kP09PfuwCdQGwwnlmFbHeJG4Oci0i8iXSLSLSJdtQxsEOEUoZDPE0RKHLueEelXYWqFmjgXP+3IIRs7OhxJ4rB9nyQRa1p7Jr1q6kusJFQ3zg8i75+Z4/6TZtK5YTMf+MaYNocJwLgz63c2Pg+cALSoaruqzlDV9loGNohwijBUKplsB7K5fv6CH22pDdf6vHGSiqTZaziDjXqW0uqRNCYhWIjitD2bQ6k1x6NHz6DvvtVcf/fkp0GlIXy+xbju8BBwj46jalbDMDNF2Ni9jSh0+k0qPGYaxIzqFRwZgcfd0veaGHgcBK2gUk3OzXBahbKt4maqNyq9swtsXlzkJ9f+gmccfuy45lkr6lT8rMRa4AYbwJ0kADdcFHWMvIQeAUjiVsg0i/FoMJTqbodqeqHre+igqO3kG6fJFpVKY3JBvyCV/8pwRmd9FYWtC4rM+tc24jgmCCZHqDLW0V1CYFtht4LdakaDCKcIfaXBTK5f5Me6yHCRK/L7E5LV5wSva5M9X9HEoR9XjK2OCrJTY5SpZvgB0xexa26BYLCLux66n6MmsflonYqfGajqReMdu0v8xExH7DV/DwqDMWFpuKl/rEXnREXHOVOf4XAkxOjpifYiqZXVvg7ExZm6SJwR7u9yEQshPXML/PYf148+4R1EPVtHReQjO3pOgxNOEZ568BFckQ/o2DLElgXVy0r4kmMCTXU9Rzyj0qwVO5246wgojTNVj6tW00JdlocmInHStRtlw6IiD91936SJpCaRuX6IrgrOHsMVIZhyix8Z6YQGJ5wiBEFAbul8Fq/orcr6HEeqLg6O/D6gYsHat3EFtww8lmgS/TUzxHBGPyggOxdnpdyyR5HcYMTN991V7TEnBFplqyN8A5gxytZmzxkRDU44hfjQ69/Chz78IeY/0c/GRS2pPqdYosgSlB8T6ghFyXLDasYbh6zO6Ahx5BjTwNNNnSsjdX2YuUT5gM55RX54/W84/rCjtuPpa4SC1nHY2o7ogg4NTjiF2HPufGYcvowFq/uH+8B0JOEwXZBxQiDDMYwj4gjW8BKpgafE/mQq5uATe9ecPP0bto56rR3BeHXCHexPOGpj3IlEgwinGKcfezIztpU8K4jhMwk38onOM8I4VC9ZUT2bHlLOmXC0UdazHygnTmytcu3+tpB4ay/j8FPXhPE462tsUOv3JzwH05/Q4UuM3hh3wtAgwinGcQc9iTBWmnsjspkTQaZ8ReXCG20dVjvmXBg+MY0VB+3zysSdUuXi5UJAUNZMVM1Ewdx3XJxw3P0Jvca4lwOo6pCqbpuoZ6pEQyecYjQ3NxMHQr7CVRE74dE57UZAbes+NcrU8qvrdEet8BUa44xaH6IkHDeI1IbgTQIUGJ91tFrvweNqOGcRpv+Ea4x7BHAbcL6qjhixLiLzgDeSJvWa6au+YayJNjhhHaCcE8LyaLyt+iIcmwAN63KZEGIv5Q+rdgmt8jo16mTFY4CwFBPlJs94MoI4OlaT0LEebbRzam6M6+GXQAfwR+A33jYmGpxwihHHMflSTKmQZjf4K8V/n7GOysjH3N7AW2I+wVYjsvSdjR0Vj/PhuKj1NRJk0qEK/TH5tmYmBzKSdXSTqlZ2CvOxI/0JleGNccciwhZVffcY51RFgxNOMbZ1diJqSkzASH5A8//w8obpT3llvKgjwOGZFOn+SgTectAq1ll3/9jWMHUdppp6IwaaJrF48PgchePuT1hLY9wq+LWInFbTzCrQ4IRTjPYZM4hyQktPmcHmkb+OWnQ/vxBirCkh+iKle1+N28a4TIt0VMJhdbijAkv4szaVeOrJx489wfEgIfbtHLYD/QktXGPcAvBoxbFqOB94n4gMAa7KmtaSU9ggwilGLpcj3HMOczb0sXVeNnxttJA05zivJCZwgWimeYzvhhju/k/3O1RyWxemlhpUs+6TOWv6aeopc9pxJ436nDuEcXo+qjWodb0J7WsF3jLC2DsZ3hh3tHvNGN8sG0RYF2hdMJvmR7uH7fdFUYcAAameaVi5L6wg0pGqsPnwDTC+38/V9fa56p4P9bD4oR6Oev4zWLxgTyYP9Rsx40NEXoRxbQDcoKq/rmVcgwjrAHu2zWaNPubtSRe/n7rkRMYa40WIVRMXQ7UIGh+OQIeLnHZGmp3P/JV9LHq4h9Pf8GpecOzJNc1o3Nh5NaXGDRH5FPAU4Eq763wReZqqjmXQaRBhPeDuRx+gnM8SiST/GWQsmlqrjpiKoGNxQP94pRjsV/YOEJq7htjn3i5OfeV/Tj4Bjt9PuLNxGnCkqqkaKyJXYHrdj0mENVlHx4qjE5GDRGS5iAyKyIUVx94uIveKyD0i8kMRqc92sFOE2x+6j6aVW3l83zYyJhPx9TAvH9BtVa7l8S3wOOb2qlSVumNq7jHcdf6qfkoLZ/DSU567nVceH3aRGjMAM73XHbUOqpUTuji6l1prUWUd9C3AW4EX+ztFZJHdf4iq9ovITzCm4u/UOsHpjouvvpLNC5vonVkYRljJWqtCcdV1wuH+weyI2jmK8xH6XX6dnjhn7QDPOO3Umq+1w6hfovPxSeAOEbke80GfBLy3loFjcsJa4uhUdYOq3kJqmvWRA5pFJIch3kqH6W6Lbd1d6OObWb94IhzdfmXS0c9LOeXIzresW8LFbCpx2TSQ+ct1f+CH//hDZszj657gr/ffMQHPUjllGb7VGVT1h8DxwM/sdoKq/qiWsbVwwn3Zzjg6b2JrRORzwONAP/B7Vf19tXNt2NE5AHvvvXctc9/lcdGPv8FAc8jWucUKHS9bL6ZayFgWWiUjYvjokQw+MIKIZwyxacqUChoIt58wj30e7OaPV17Nc590LN/54zXc+bcbaeozFcW7X/8qnn/CKaM9+nZB6pgTishBqvpvLw1qtf27p4jsqaq3j3WNWojQxdGdp6o3iciXMMrmB2uY4CxMpPo+wDbgpyLyX6r6/cpzVfUy4DKAY445po4/9onDpjVr2TKvWEE9mka7+BJklQgW32E+ikdx2J5Ki2s1AnTGmMpjIkAIjx7YxpyN/Vz4vvcQA6uXtbF5QRPz1vXzs+//iAOXLmO/PfYafuHthQrUcVIvpi/hOZjiv5VQ4JljXaAWIlzN9sfROfwHsEJVNwKIyM+AE4FhRLg7Iq+Chi61qAqvq4j9FE+v88PKZARWMVYNGkdoI50j9iJVCTEQHjm4g44tQ6zet81ad4WNC5vY94Fu+gb6R3rs7Ucd/ySrqgscf56qZjq/1mqEHFMnHGccncPjwPEi0iIm2exZTGJy5K6GYEYzxf60OWcgaglqpC2FiDHAjGWEGXX9WktrSuzZsxOBVrIWWrdv67wijx3YTjmftv2e/8QA0awWDt/ngNHuvH0Y/aOoF/yzxn3DUKt1dFgcnR+DJyILgVuBdiAWkbdhLKI3ichVwO2YHK07sCJnAzDQ10e5EJAfKLNwVT/FwYgNi5vonV29dmx103ztYuhI16zGMTOjtfZ139pTpjxQ4rYH7uXJthbplq5tXLX8z7z25NNobhqHh6o+iQ4Au/YXYYyPR5H+pLUz3ItQFTUR4QhxdH4M3jpMGki1sR8GPlzLfXY3FLb0U2jPcdTfN6MteXIzWliwfCt3P20WfR15Kg0xpihwdR0vdeDbYDelYnyloUYyY9Njw1+PTQPp9R48tINl/+7i0i98mWhmM1Fznvw6UxHw/rv/xeff8cHtK41Y/876U4HXY9a/X/K+G3hfLRdoRMxMIRYduAx99DFmH7yUz577ToIg4Iz3nc/8NQM8PitPnKm+mxINUJU7JVxtBPEy3Tfyoh6p4O/IqHBt5AIeOmwma5a2MGvTEMXBiE3HzaFcDDjin2t5x+Vf4ItvvHCEa1VHPVtHVfUK4AoReYmqXj2eazSIcArxv+dekHn/63/+meLmAbYc0gFo0pzT8jaylpoKvjWKcaYSqtlrVSu7X81yOpo1tRJ9bXn62vIZHfKBI2dxyK0PcfPT7+bYg7ajqW0dE6GDql4tIs8HDgWavP0fHWtsI6m3DnDXyn/z3xe9g2u++2NWPamFvtl5rAESRKnu/8vucxww8BZ9lisOH1OJOBFhpSqxxVopuvqBdWPfZuvcIuv2auFrl17Ktu6a+meaq+vwrd4gIpcAr8DYTwR4GbCklrENIqwD/OTPv6a/t49HntnG5v2KhBLb1NpK5/pwp7tDhmMlBKjJex8jiavOnlotPnNk7jcSAWadkW78ioPaKeWF937xk5RK1QKsqmAXiJgBTlTV1wFbbUHgE8iWzhgRDSKsAwz09jLQEdLfERLa1ZpytKzYaPZLQih+AxdzEgy3c6auj+riqRkoiYlUqx539/XvP1ZgtfFFpudHgXDvUbPo37SNN3/kPTy4asXoH87Y3pp6gXOM9onInpgQzn1qGdggwnpAGNO2ocx+f+7m8F9s4/Bfb6XQE2VJyS7iuGKxVykHnAiUgaSbIcDhjncfcZzqfEn5+2HnSxWR1KNs9QnUxZyS1KMBYagp5PanzWNzrsxnPvEp/nrXLaN/PrsGEf5aRGYCn8W45B4DfljLwAYRTjGiKKLrgXXkB5SWGU0c89yTCWKlZZspBqxa2c46q4dViqHub6Uo6ml62WvZt1Gcku/wa2ffu9Bv/6h6d6hKnBkIpXzA/UfOYs3SVq745re4f9WjI5wLEg/f6g2q+r+qus1aSJcAB6nqh2oZ2yDCKcavb/wVMqBoqCycM4dbr/sLQQxzHx8gGIqpzmWyHG0k/S1LlKk4my1b4bgUGd0NmzkRq+lZH6vPhPz/08R3X2P1tFnICszeflh5wAy2zcrzyUu+MvKHtAtwQhG5S0TeJyLLVHVQVTtrHdsgwinGH2+4DlFBmmHPPfdL9ndsLDFnTdL63BAeTseqFBsrnfqV5OCRj2YNM0JqFXVcLEvkKVdLS9B73BKq/iAk7gyy862EIjx6UDvFzX384bblw45Xs4zWo3UUeBEmKuwnInKLiFwoIjWlAzWIcApx7c0/Yuhxs6Ke+cKTOOf0V9Izp5AY/9o3lTK6meNMjhhiDdCMqOjemZTcAPMF58QZb7KGGVWIKjMUPJEz4WDVjC4KGpstnd/I9tske9HjZI4DDzTn2LSwiSv/OEJdpF3AOqqqK1X1M6r6ZODVwOGYHvZjouGsnwKoKldc8yGWX7seEFr2LvHyk18DQG5mK2wegj06aO3qSa2RVXQsR5gSKCBEEQSBEkq2NlqszigTZMYZQ49n8bRhOGotpRmCUhKvfpZLqplDRcjcMLcIjhDVu2hKiIPFkA6pHjNbj+JnNYjIUuDlGH9hBLyrlnENTriToap86fJzWX7tBkAotJc478x3Jq6Is05/pTlxbSfNvRHF3ogoHl0HjGKIY8slkcSC6hz+YQChaKILahwSxU7fMxw1VgECQxQKceQ2U4ZePWJ1jhEzF0m4ohNZHees7r6odPKb11EuYEt3dTVqVxBHReQmTEZ9CLxMVY9V1Wo5hsPQIMKdjF9c+1Hutxb5BfvFfOQDn2DfPQ9Ojp904JF0zks5QinvviIx4p86MdIQQBRZ4qlY3JGNawtKES0bh2DQcKxSWShFEMUB5SjMhLDFEfj6n1kenpshhiiyrowqeqPa+IJEd/SI1qcbzYTgmSPdHXmKm/uGf2A6fuvojjQJtcdDEblDRGqpH3qGqh6tqp9U1ZFNvVXQEEd3Iu67/1f87ldrAGHh0k4+/I4fVM0oOPqYY7j/z8u55ZR5aOAIwfkIfae5mvITsXHqly336hgos+cjfbT0lGnujpAY4gBue/o8SkWnKToisISoQmz1Q6nSqnuYKGz3ZbM6Ur1Tqo3xidZznQBsm1UcWewcB+fzmoQ+G5OYfouIXKOqfi6s3yT0OEyTUL992vmY/NcxS9kDW0XkcmBPVX2ebUh6gqpePtbABifcSSiXy3zvez9EVWhpHeQt/33RiCk9Z538QoqDSrE/dVGoBlanMyKjxoLGAbEGRFFgCSgg0pA5jw0xoyvm4KMP54WvfwUXf+3rlIo5a+gx4+PEoKJGlI0tJ1XQOEiDAxLullpi3f5qltnkB0Irxnlc0dw3PR+EOAxGbq82PhfFuJuEAojIYuD5wDdrupupIPg7wJUifxB4Wy0DG5xwJ+GPf/gyWzabdgWvPeNZzJ994Ijnzpk5iyiA/IDSZwuxxRojBESR4VJxbAhGLB1HkUCkSBDS3Ftmxj6LeP8rzk2uKTNaKPZFhgAjxwlNiEwcBwQhiSgpIonVEzE6pkhs9MzAnoMhYGfAyUS8ifEdqkIQpCai1FjjDDqkHFGV/FB1OXMEHXCuiNzqvb/M1ily2JEmoWuB/8MYVmrtMTFXVX8iIu8FsA1porEGQYMIdwq6ujZz3W/vAfIcc3QXRx/x2lHPFxGGWvPM2Faic1ZTYmiJLYdCMcWPRNAotoQBaIBEyozOEk86+eDMNdv3XEDfoyuJlgRorEhgdT4FYj+DwlxfRXyDqHF4qGaIUwInjmrGVZIQJIZjOoZfLYjAYdHjvUT57RLMxupPWI2tVpJz1XNE5AXABlW9TUROqXE+vSIyx93DtVqrZWBDHN0J+NNvP07/YB4R5fkvqK2P5P4nPJmlD3fR1FWy3MsRiCUeGwcalwPiKDCiaTlg9tYBglg54+TnZa43f+Y8wpISD4ZoOU9cEuKhgHjIjEWdZdRYWDUyW1z2fXMCiUicWk/jKKjiyHcisye2epZTSP829ZXZ58FOnvuSF1f/MMYnju5Ik9CnAi8SkccwYuwzRWSs4mQXYPodLhORfwDfxaQ1jYkGEe4EPLrSpOycenLMnouOrGnMB19yJl3zWtj7sZ7EyOGc4XFkOJJGGEIpGWIigvbNQ/TPbKOpWMxcr7enl1I+ZzioCloOjbUmDiAKzH5nfY09J6HPK7zgbKc/DtPxnHU0kSzNa6dbJuM8X+LsTQOU24q87BlVqnqP3zq6I01C36uqi1V1qR33Z1X9r9FuZuuLnoypJvjfwKGqenctE22Io5MAVWXTps1s3tzJ3/55Gw8+MsRpz30Kp//n2dt1nde/8KX8/PLv8theM+hpLaBlY1XUKIfGsSUeTA+0GIiFsARBfrjTW3ziiW2Re009fhoBgU88eNEy6v1cW4e9fS0BOMUwsZRaOdbXAVU101EqvaQyc8sgMmcU1Wsc1tEJaBJaE0Tk/41w6AATIKE/G+saDSKcQHR19/C+D3yOwb6NBEE5c+yZz37Vdl/vtKc8lV/8+Xccddcm/nrcniQ6HEA5AEuUBFgihJk9A8w7YL9h11q85yLW3vEvGAzQnLWIiBop03E/MNRq+6lp6LXgjsgSIkb5M+NSZc8YitT9s/qlGgONWu1RjN5LrOz/YCezNw3y9Nc8v+pnIIzfOb8jTUK9c24AbhjllBeONhzjwB8VDSKcQHzmM5fQ072JTdtm019uplgYYtGsdex/0LHMaGsd1zWffMTR3PLYdey/opOhMGDt3DYGykVEA1QUKQsMKITCjL4hOvoGOfu5Lxh2nf937Incdc21FEsRg2L6yysCNuRNEoOKoJElxDJoiCVEwfk1NLD801lhfMd7jNFXYzC3UQJr3WntKxHnoJwTjrxzE+1dQ5QLAW9+21t58gGHjPwhjJMIdwZUdbu5ZyUaRDhB+P73f86G9Q+xfsNC+kvN5PMl9py1jtnz9+edF4z/ezrjGc/j78uXM6uvRG5ggH1XdfGP/fZmMAwIMJbOQI2hpHVwiIFcjkP32XfYdR5duxYFhoLQ6mSW5AIxXE4NwakzqgtGzByycmOIMeAY34Sxz6ikGRuWGP0UYyed7v14J0tXdVEoGSKOBcqzWnnj287hkCXLaG0ZpTyn1meYWiVEZAHwCcbhrG8Q4QTgttvu5a9/vY71a+cx2N9EEEQsmreeOXP34lMfuyBTomJ70VQs8r2LPg0Ynes1F7ydxZu7eXTuHEQEiQUpAShlDchHMXEcDwsEeGzDekphgEaBSbhxomYiZipEjqcFUHKiY4DEoBrbIHBN/RYBOH1PVDKiLACxki/HLFuxjWWnHM+bn/9SBocGueG+Oznl8Cczv2NWbR9CHSbxVsF3gG8D77fvHwR+jO1mNhoa1tEJQGyCLunubkZiZcHczYSB8KEPvH2HCLASIkLY3kHzUIkgMuIiJdM9N9SA3rAFQbn0l8NDHU867EkUopjm/hLEGDG2jOfUE6NjOt0wBomtSBrD8PQha01NjknKYSNrKEKYv6mfoUKOd7/8DDpmzGD+nLm8/On/UTsBwi4RwI111mOfXFXLmJ+4MdHghBVQVX5y6dV09nbS19NH+8wOjjzuSRx61CEUK8z+DrNmtiMCuahM66wB2tr6ePvbL6S1taYq6NuFQq5AEPURDLnUIzFuwyGI45BVLbMJr/sNH3zgYT584VvJ5cxXvHjhQvpzOWZ397N2dsGIsWEAJSUWy8mwbsDIrPLAJRJHqQdeXKC2eD4GZy9SSTmifR8HQhgE21d1uxL1SXSVGLezvkGEHlY+soovffrLdEedlLZGaKQETQE33vN34j7h69/+CoXCcPP/osULAGhuGWTugi085z9exCGHDLdQTgQGunsoW8NK4LhTydhXROHx1rns1buFDSse4Mzz3s3XPnMR7TPaAOjYa28Wrl3HuhkdBJEYHTAnhLF5DYoMYnKfrOXUuCnE6IqxM8SABEFi7XSWUmeQIYCmoYjDHtlIe+8gYdXAlBpRu3N+qlHprJ8HvLSWgQ0iBOI45ltfvoKb7l7OwIaITf/oo9wVJZYFycNeL5/Jt778Xc69cLiv787b7kVjKMgAWi7wsldUN7fvKPr6+2nq2saWWXuagnp2gQYxhBoxb6CbzriVm3MHUCxHHBKt5LMXf4//fdebAHj2iU/lmiu/T65H0YIRPbVkjSyQuiZc3HgAmjP3MFGQ1iLrQk89o2gmWC0W9l21lZYY/vPVr+T4gw/boeeuU/EzA1W9XUROBg7EfBgPqGpNhVV3eyL81+338OXPfQ3NRWy+qZ/eFZH1dwWoKkFTExrHbL0r4o6WW3jjGbdR0BbaZszgwEMO5NTTT+Hbl32P/u4coZbRwclRs1ev38C7P/Ep8hKwLWwltBZNx4QW921lae8mcy5z2RDNYj2zCB++h0dWPMbKNZv4yw3LiSQgPyCUEENgoXVAlA1X09BJmQKRZ+sMXLa9JTbL8ZIyio4QAyiWyize1M2rzzmbZxz9lB1/+F2ACCHRA+/d3nG7LRH29/dz0QWfZkv5Cfo2xmy5dYioP4ZcaM8QgkKeIAxQDelbqazZNETT4iL5Gb1sK/SyuWctN919A0ODAVsfaqV14RDF5sn5SD9/xZVQLnPnrL1NlIsLA7MhYauZTRgqe0WbWcwmOoJeHhpcxPxCJ5/51CdRhX6KrCosABUK3VBuhii03K5go2fKIFYnjENDVHGoUCAlBicGl42FVAIbyWOJc9HGLgZbmyeGAKk5TG2XxW5FhENDQ9zyj9v56x+W88jqB4jLMZtviRncEKFlkHwewtBYNNWUgzARWObXvjyUo/fxwKQFhCEEUGhXhvrzIAHRYMBgaQI71HpYuGABW1Y8zGCuYKySkdEDtQx5BTRgVWEem0odHFl+lBnSz8LcVv49tJSAiEgCNBBKOQiGDPfM9ZrH0BxEZYgLmEAaF60aKxqIMdAMxmgu9c1r6E0uBkLzOeVLZRZt7uaYp544MQ++6+iE48ZuQ4QP3vcwn/30F5B8xOBWGFgHXQ8LOmSJDCCXS00INlUHCVARyOVc0U78TitDPTk0DBARBjpzzGnuZ+0T69ljzwUTOv8Xnngcj/7zLxQHSwyq1+3IckM37wEpcK8s4TBdyZ65zRSlxMrSQtwZ4RBEobGqhja0LAohNwSlVAgwxJiImqZQFGWPCMUQKEC+HDG7u5f523pY0NnLUHORVz67SjD2OOCrnvUM24n6NcC+qvpRW+5woarePNbYaU2E657YwJ+u/Qt33XYPWwfX0b8RttxBEghNEEBTzhBWPkaiiCAMiMvWvSOCNBXTJpyhreUShqkdIhDIhahAuVyg1Bfyncuu4r0fGTUkcbtxyP7LGAzyzBroY2O+w4iBJTKVdx2B9Ekzj+kClrKeOWEXc8Iu7i0tpZ+iWdQDgO1BKjbmFIGwD9TvTSqgQVpdLbB6oOaAIWXWUB9Lt25lTk8f5XyOobmz+K9XvJqTjjp6Qv2juwgn/Drmk3wm8FFMk9CrgTFl8mlJhL09vVz25Su475G7iIdgcGtA/8YcfWss8YVxYngRy9XCfJ44jonLEUHBEmI+B7nQJK+WI1QCKOTTfn650KyPfGCz0ZXONS08WriHvr5+WlqaJ+yZ1m7YRC6OKBEaS2U51ZUS66GnJ27UmSjCYjaSl4h9cmu5L1pqrKlifIwupyHfD+UCBCHEeRIRMA4gN2TOiRxxxqDlmCdtXsfc3h5m7beMFz77OTz98MMnlvA87ArWUeA4VT1aRO4AUNWtNoVqTEw7Irzp77dx+Tcvp1yK2PbvHP3rjEgZhNa6F4LkQrQcGUIE49/C+r4KARrHSC5niBCbdR6EkM+lY0MhFkmCmU2tJKG7u5WOqI9vX/5j3nLe63foWeI45nd/vZGHVq7i9hv/zoA00aWtqQjqwjZ9bpgsWGEzM9kctTOLHspWiRPBOONtrCg2HDRXtnaWQXONGMgFEBWtA99x2hiWdm1hTl8f5114IUfsNzn+0Ax2DSIs2eJSzlk/jxoD7qYVEX7zKz/g5rv/QveakG0PFGwJCHMsjpxMZQnP/pVADIcLTbY6sY1/RNMWSCJQzJOUEQsDw3SaLJHGCsWcKQmRC9i6uYM77rqJ2+94Ckcfdeh2PYOq8qNr/8jja9bywL/uIBjqZ0hybJM2Hs/NM0YTF27mOKF/gcgSos0AlljYxowkHjvJOgJCmxxMTpAoDSFVbEkKhdyAEVEjNWKoxDH7dG7i1Je9fKcR4C5iHf0y8HNgvoh8HOOo/0AtA6cNEV7/u79z891/Ycs9OXpXhcNPqBSV8jmIrR8sDIzhBUy2QCCJaKVDQ0iYw2QIYFZoGKCOaAHJCQQmRIswoLurjUJziYu//hU0aELjEvMWLuENZ7yKA/YbuW+kqvKBz36J9Y/8mz4p0hM0s7q4iIgwERGT31ZrjPHbTogLwE6sJzbYwB33CEzjlI6NYCqEEcRF0iRfMX8lIrGMamCI9fBly8b+UiYKuwAnVNUrReQ24FmYj+rFqnp/LWNr8iyLyEwRuUpE/i0i94vICRXHDxKR5SIyKCIXbs/YHUVXZzef+NAX+OFV36PzkRy9a/JmlYkxmBAGhuCs+wGwJamNeEkYpvvTSScvw+Ym4nxI1JInntFE3JQnLuQIcqHhfGLE0thy0BiIA2Fd51ye2DKf9dvaWNc9hyfWrOZLXxs9oP7SH/yUdY88wD1NS7insJSVuQWGAP2pVWwJJSXOe/X8eSkBuveBE2WjdEwi1qoiQy7qhVTkLZuFko/L7Ne3HgVamydO3x0Lu0IAt40VXaOqX1PVrwKrRaSyultV1MoJvwRcp6ovtcpmZWTyFuCtwIvHMXbc+OcNN/Gt73ybaAi23l9gcLM1osSaipcCSV3AXAhYwrQLVvMBEgYQx2iMcc4HgfnVj2PKKDTl7fVAi+YjK0PimkBNEHSQC4kChZY8AD0DbaaMSzEgGIwpDFTv066qfO9n13DrX//Mo8WF9AZNZIrlubAxKyImxOUf90TNRBT1jkOWIBPjriW2ONQ0Sb+kBDZQO86Z6hlSKnPUwKPkZrTzhrP/h0ULJ9YFMyrqkOiq4GLAr+DdW2VfVYxJhCLSDpwEvB7AFlId8s9R1Q3ABhF5/vaOHS++9plvc9dDN9KzOqTz4bypm+JVfyCXM+6GWBNidBZMd5KaJg1GvFRDfCYDR2w+iiU8SMO21OqRgOQCozfa6wwFAnmPcwUm4kRFiAjROFvyAuCBR1fyua9eTNzbycr8AjZLe1YHSgokkW0EaKkowxUcF5QskSbE59Ri7xLuejm1hBibaJw4VMISUDABMm3BIAHKJZ/91KRZQUdCPXK+KhD1ypGraiwiNTG5Wk7aF9gIfFtEjgBuA85X1d6JHCsi5wDnAOy99+ht3f72p+Xc9dCNbL4rT9+69BHUGafyYWLx9B3riaQWCFLIpVkB1QphCmnj+EAIkhqboLGiBeMrjEMhzgdoKMbA4RzcziBkq86LxAQ6xJlvfjuSb6KlYzYnPeVJ/O4319AlTaxo2ochyWfFRycO+pzMP+50NrfTnRd553vjM/tIF7eqzbAfMppvoKYykvMjak5okogol9vpBJjRg+sbj4rIWzHcD+DNQE09KWrRCXMYlnqxqh6FYbPDmmvs6FhVvUxVj1HVY+bNmzfiBdet3cAV3/8OXY/m6Fvrc53AiJu5XFJtVmNNW9qBCVYu5qCQM453SdeoWj0ydpKqtYhqwUTExLEmRKyBEIQBUSDETaGV18QjfHNeOYepYAZ0lmewanABa4fa2NQPm9Y/wR9+/Qs2SjsP5BYzJPmM0cUntOSPt9+5KRJPeux1mBCP2JwBxx1TNZW6YzU6oKoRSctKLoLQEn6gEEQQDphz+uMCOS1x7/0PjfjdTAac8anedULgXEy5wzWk1b7PqWVgLZxwNbBaVW+y76+idiLckbHDEMcxn7zo85S6AzofyqVczhlirDtBIQnEllxoAovVio2QMg2xoihiEspzgRVNk1yBlItiDBcA5eaciSvN+5lymnCJWK1zOxATIJ2HmJBN0mE8HzGQs1wyTEsOSoXbIbHGevAXoAviFtUsl/S4h7OImmtrqj+SvZc/NnFxYIgx7FY0ztPd3spXL/k+F3/pojG+qQlGfRJdBlYle+V4xo5JhKq6TkRWiciBqvoAxgR731jjdnRsNVz5zavpG9rG5rubAEt8lusl+mBg0pDAcDON1RCXXVSRlRXDQCCfR7FiZi4w4qcjusisUAmFODTXjFHjhshZD7fNIHA+SCeGxpYDpizIXDOIDOElETd2+buKhRnR0RlWKhZg4CyxHkcYRoAVMFxU09f+PSQdJ57bIrC/W64KmwCd/e205jekkUY7CaL1T4Ui0gScBRwKNLn9qvqGscbWmvx2HnCliNwNHAl8QkTOdYVURWShiKzGZBd/QERWW6NM1bE13jODUqnE35ZfT+cjOSKXqCDggos1Bg3DtPGBO47ReTQwnX9oKkBz3oiSGP+e0yEz1orARMNETtR0BwPxVEgxTuxQiPNCFChRXohzQpQzgdGSr1isYs5HJInbdPl4Cc16xJXZBGLHocyDkfzxi/ZWGmecWKpV3ruxkSaLPYgVidT8tQELqkrfUJEwiHhi7foav7UJgI6w1YDx9icUkb1E5HrrUrtXRM6v4XbfAxYCpwJ/wZTU765lnjVZb1T1TqCy+YZfRHWdvWmtY7cb3/r6D4jiiJ7HijaVyBJNIFa3k0SUDIKUWzhCDfI5U3UniXqxHM/rGmQsmfaGIkjOEFmIC11TNLB9GAJjPCUUE4tpa3CK1UHd9eMYxMaIx6EjwOQWhgjD4VzPEZ1izkle2/cVKmM2ftTBuTk8cdV3fST3jLwbJ+4N8SxZdv6ElKOQ3//pZs587YtG+KYmHuPRAWXH+hOWgXfYbPkZwG0i8oeKsZXYT1VfJiKnq+oVIvIDTPXvMbFLRMz09fVx613L6Xo4hxDgChxJIElGg2J9fKpmTRUsV7RsJorVJs+Rsha14qQqWkgjYDQ2hBxZYipDIl6qiC1spK4vC7El/AhNiNnjncSREhfs+AqLbEJcStpdCY87OuORHx8qyW9LgkSHJDXiCBhO5o5H6YJW95weF05uHFdcMDkuDJSLbFi/djxf47gxzrC1pD8hgIi4/oQ+IZ2O7U8I3GgDS/ZQ1bWY9mioareI3I9pmTYaEbpSFttE5DBgHbC0lonuEkT4pU98g2gIetYY3x9lWwka0l/5MLRRKxijjNPDwoDIGlF8ThMEng8QPGuFESEjxHCtZIDjnFguZyjBVvajlMSEmR3OGJK4K5wNCWuYCVNCS8VmNzfSIryV4pcjQJ97eZwxc65mz8kYdTwdMD3fex8bAlaxRGAPDJYKPLpiJTsV1TnhZPcnBEBElgJHATcxOi4TkVmYeNFrgDbgg2OMAXYBIty8aQsr1t5H54N5wDrfwwCKBdO4HSsqWr3OLHy1upvlcsUQkERERNX2T7HcKQhSkTIw+2Ort6WtT0iMMQlBosQ5dx5AbIjeWjOiJEouMLqcI8bQEq9VX32uqZrYUDLMKAOt+CuGU0pc0S7CEnMYQFz2BADH7TLXtM/pDDLmcdxvSnKvfFBm4YKFI35fEw4dURydtP6EyUGRNkxO4NtUtWq4k4icr6pfAu5X1a3AXzH+8ZoxOVWJJghxHHPRez5FqVfoe8IZTmzcp1iDSi4EG0pmmlYa7oclxrgYmrC0wPu8RYhDICfG0pkPiG2Jvyhn/I1iuZnzAarVATUhOPPaSZYSALmAshjRM86LjZgRY5E1ErQhcCR16klKHD4nSqbrxEpHCHFKP+KfZ+eQWWUqRopWSfKQnVXT3F6z3LJsdMaclYEl6YOmzGjuZo9ZG2hr6mXJkiXj/EbHifEZZnakPyEikscQ4JVjdFY60/79Sk2zqoK65oSf/MCX6Is62XR7ISU+wHgnbGFa1WS/xo4fYQgnFJOy5HOM5LVkRFAg4X7uWxbnbLe6XwRJbZXUkGMWbFI20NfXPAJRQEIvmkazxxyMyGoNSx5c1flkuL1GwjWdzzB5GGu7iq1FtZxywcCNcfJyrCnXw6ZmWW4ZCLQ1d7Ng1ibmzV/Kkw4/hVe8NNuAdDIhjMgJx0LSnxDjQH8l8OqKc64B/sfqi8dh+xPaUhWXY7jbF8a4z/1imonOtx4Af+qqqoePNdG6JsLHV6+g+/GQcr8VP53hJEgXqflh1MQymrgoxISmOb0xOTmAKAxS4wZmYWoYIDlJFrViqrkH1gpb9iwigRNzLXdTgTBn5yQk6T4ZOsrbKDlP9hhpbTk3RCICWw7o9EwgzZrAlr4JzL5UD1WiyFh2HeG6EjmOlTo/IFgxNNLknoFnCOpo6WbWzIU87fhjGegvsWnjFuYvmLvTfIUSbz8V7mB/wqcCrwX+JSJ32n3vs63WKu/zKhFZaO8zLpNxXRPh7Nlz6dm8LmNoUZuipKqEYUgURYbbOfFUAAkSqUUEczyOE8MfOTGuhoDEh+gaZDqRLiqYwen3nxpbYjX+QIk1EWvLGE4Xp0E5CZczhOrlI4pHYGQZoy8dJroZ6QnqE546o60X16rpj4R/fhogkF5L3AHL+RKR2PwyuScmDGK2blnPz39xFaD8/k8/IxrK8YlPX8T8+XNr/TrHh9rFz+FDx9mfUFX/TnV9cSRsBP6lquOyWNU1EYrmCAqG8yWFzgo5yBkRMwLUlaUPjMtCRIw7woatuaRVzYeGWzgWiNXpKvyEMcafRxBkrJsaYiNx1MSRhoLmJCt6eqKlLzamFa5J7+2d5MYEXvRN1pJrOW8SB5oWX4oTgkkl6YSTakp8am8olnu7SmvEml1tAkRxWr8milj96EJrqFEoRxSkzMIDN/Pzn17Hf79l1C7SE4J6z6xX1UhE5opIwWYKbRfq2jATSZmoHEI+b36sbRSL/+uOTaxVGxeaiKlWTIxDsQRojTChJMcgXaAKEHohZyRqkbFySkpwGooRUwHnplBswHYOa0HF3CsQxDr0K53xWXeCV4LDnuBoNiEqp486AvGIWjD3TNwe2PnayRnfanpttbqgH3njdETnughEkgB0sa9FISqF9G5u5qYbx6zmNzHQKlv9YSXwDxH5oIhc4LZaBtYtJyyXy2zuXE9poGgMLy5LnpTLxGrLyORCRDXV5wRw+p0lqCBng7SdIQWj70VqezHkA3u+WfSCDbAupOdHqklEjMbGeRYrhng9gnCIrbyolgicCuVE1Qq6xO42c3MuE48LJOKpPcmWmLJ6bYX4arM+JM4ahsmMt29VIIpthzMhRo0umci+oOUYhsomuiaK6dlQYOaePTz+2Gr2Xlo1WGrCUKdZE5V4wm4BMGN7BtYtEV7xjatAIno3mDIKpqqZVWFsnU+iGLFGFmwen4RBQoCeTcPCxX0K5Ew7Ps0Zwg5sRr4qNpjbinPWCKQYsVJzRocSW8rCrz6hQBgYcdhx2iBIkhcycodvPc3oaxbuhyC7M3uvJJAbwyVdIxf33EYfJvkREe867rWqImW1C90+v8eRw0QvVkvQAmFAuT9ksDvHV79wBZ/58vuZNLhfmDqHql403rF1S4T/vOlvdG1sNW3A8hiflZj6LeRDs3ALuYxBQgPQQmqUqYz2V6sQJYYSK1/6emGcIyVUlChSApdjmLPBABUEmEqAJtZUcs4Io8PzUS2RuECCZCde+hTpCxFDXL4hUgKBkiYj3fMHVilMqgc4bltOL+lHySTE7ymyjm79OYq9UypImD3bHmuh0PY4f7/hFp52ysT0naiGetcJAUTkeqoIyqr6zLHG1iURdnf3EubLdG9uTcpFaBwZq6cTSRVMDz2xomSVX8u8LU2Yty4JR6jWiOJzSkXRotkRe2kMKlBOArttXTJPrBTf8IIhdLWN/QJbYNg/nozz37sfkiqfRaL/kYqpYLsnRdlBqkbHU88gIzhxLkPiWUuprRogseGKgbWCxVbGDRTjvnB1ENUQ+sC2PL3ri/zoB1dNGhGm8697+AXOmoCXYALBx0RdEuG2bZ30djYTaw7EcLM4yCE5E7ZmaYg4xgZemwWfGDDEhqRZsTOxEJKKtUnGRGCc+s4J7/x4Tnx0XFIkjRNV3Fgb6+xF0fjkpVZfc/7J5MLipSRJhT8R77xhV0zHOfYrGBE4jm0NUU0/HxXzGSUSsw3rQ21MqH8PK46bPEKPy1Y5h3KcuDd6NjTROn8bnZ1ddHS0M+FITbt1DVW9rWLXP0TkL7WMrUvraByX6etqMtzKWUNtQSbJB2hO0i6ylsCc3ha77IfAphuRdcCDdUHkgkSPdKs0sPqPb4TLWE/tvtheP2Ew/vnOMgqZawUuybfyWdVaHj39MJMb6O7rJhObLRSPI3tE4+7nfjz8e0o5FVVdsSrnI1Rbldvp1w6BS70KrevH+iRDGxQx2JUjGgr45U+vG/5wE4RdobyFiMz2trkiciomv3BM1CUnBGWgp2gsdVam0mQRYH10QaLPxbEpXZhWRLOW0LyNoslZP1vO1JURMZZOJ5a67zTyrYHeYkyIzgveTgKy1SM+LCe13C5IOKXhioG1PA5/XF+mtNchvXaii1nuJWQNNxqbLUz0NSUuq3sUc29rKcr86lpiTDLX7fhUVzTXCVwkTZhWKogGyskcezfkuXH5Tbzu7JeP+c2OC3VIdFVwG+nXVQZWYDLtx0SdEiGEzRGlwULW1GfFOLFlKIIwMJbIgiRPn17A/AlCQ5AiQeLzAhInepyzhGNTh2JbG0Zjc44LCTPiaSqa4hx59qUjFCey4qQoSWlMR1hNiciqiYqGRxeJE941cvGJPwgkCTerlNoSLdBQbXIvozfa60SY0h2xYTGxmMx6wOiIpCtL7MQkFMiH6GCZQGP61hdp36uLu++4n8OPOrjqM+4I6pHzVUJV9xnv2LoUR1UDCi2RMbpYE19CBKFxzovguQ6CJOk2LUiWkmSsUFYlsr46xTQ6iYvp9TQnkDel7GOwraSto96WM3Smw1gwOqRkbpPV7aqInuBJelqxuctVqEBpxIumgQjuuJCKl/bcRJz07h+43ETjmU/9f/bWSbC4CISG2N0G5vnd76CLYBAgKOZM49G+gL4NeS772neqP/SOQLG+yYqtTiAiT7Gxo+7960Tkl7ZsxuxarlGXRGjCGVOiigMrggKB9eGlZSi8tWxXoYaSLJw4Bs0HaF7QvBAFRiesNHsnzCLEZE848c9P7HWiq7fGRTJ0hPvdCGytUreOk4BxreI8x/jjgGHEqd7rYTVIrQ7n4l0TG4amYW2OQDMPiiE0dw2nawb2oFqJI5PJEUBUtuFsVhpxUTdEShBqSrQTjDrXCS/FFrQWkZOATwHfBTqBy0YZl6BOxdGYUtlOzQZXE5A4000Qt42zVCWwul9sf6E1R2K9dGlIQc5US0s4mSMOLPescKQ7You9RZyIdo4YM9Tn7cfol5WWT/fS7UsCrdVG43ifgG/RdcYa36/nJOK0+K83+Yp9UQBiMypEhFis/9MbIoqJAnKFp8BYktU0i3GJ0tjPXK2yqwLtew/SNLvMBe9984jf6A6hvq2joapusa9fgcnwvxq42svAGBV1SoRQivJpwDWYhRcK5cALmra6YfKLHZBkugNmkQZu0WuSMBDnzXX9mFFj+EiVMQ1c4LTlXi4u00nIkBiFUKt7xpqhyUQ/JNWrEojnbnDwfYJuQMZSmmh56Y9A5XURQssqXLa9yxV0TC/5PCFTbYBICJzrJPBYvGOZ6n68TI5lmFf2PnYrxMLzT30pS/fdm8lAnXG+SoQiklPVMqakp1/wd8LK4E8JymnoSoJhhgeRhIOogBTCjL5jSkhIwqFiS5Qu5Exc2RmRNAFQxGZMWOILrWXV43LJ/e1/5lgaqgbp2lV/RzVousarnSpoJpU++R0QkjC1JLvCPScp8QEQCnHZctpYySEmqKBskhR9Q5DPrR1Bmh+EAElaaZsHzTUPIgKnnvpCTn/ps0d4wB2EL23UJ34I/EVENgH9wN8ARGQ/jEg6JupSJwRoae5PdD/3PcQBhDYzIRZrNMm70hNiSks47pc3hhfNW/E0NGJtlHdV1DxJMvAJSm3fCauPVjGGJMYdz97p14VxqLZ2pOLvMJG22n57yJ+LcUlIRlRWWz8USDM5rGIp3o1jdVXhJO36a6/iPgsnNfjOcheU4L6Tgb4ifd1Fbrjer7c0sRBAIh221QtU9ePAO4DvAE/zmsIEmJq7Y6IuOWGx2MK8uZtZuaYJTcqSASJEkdrmK6Y4E2L9e4n10v511ktSnS+xaFpXRRCkQdCu/ou7XqLPOYMKDNP7gCzhUJ3wMsc9juNzSvfa6X/Og4AThePsiYInzorPoGxmhcuyh7Q5qKUuI+4aru9aZ1c+heOu6kkjUaymSp1N81cRure20tSynnK5TC43Ocup3itwq+qNVfY9WOv4uuSECxbMJZCQPeZtRETT0DJnUElCyzRt1OIWp/2rZZLFrKHhhpUE41KdVAyXFeuWQKzb0MsLHPGTGsbaKg5X7HeW0YBUjDTPYv9amgg9E6rTSf17OkbpdMMMpyONKBIxRYwdB/eJVyHpR+98O4nLI6mPmt3EWaADgRB6u5uBmKt/8ocRPqAdhI6wTSPUJRGKCB/84Ltobh5kwZyN5ELDBlzmvKu3koSF5VNOpc53aAnLEFfKvpxeBKmZXm3aQPLdeoSecI1kcHpOJrMhM38Sc71LqnVcLZPzaJGY9r3FFVX8+seVmQRuju4z0Ir9ko5JgnvsueLiZy3HdAcT3dQXgyVtHe7r11ijmRIy0Ffk0cceYXKgqUjsb9MIdUmEAEuXLubMN5xDx4yYJQtX0drcaxcQSXymhCSLwUGxHZGagtQRrxm6MoYZUcrESTLusB9ZOyDZXyGGSsW5mWWhhksnxFcF6hGVbyV1InYqp0pqpPSIOzMRyc7RSJqSzF8EU//GUpnigr2FIJcSlYoye3Yne+yxgdYZfWaemjr3BaMHJ0qjRX9/kRUrHhvhSXccde4n3GHULRECPPWEI/naVz9Lx5x9WTx3PXPbN6NeRmqlMUSCtD6oK03vELhgbbWiac6IUy5X0K8D6noUAokzfBjn81TQavCJLHHEV3JS+9efp0+QhoCyYqmffZEQrCdqZhBKwuUTd4mIDVUzupamcittrX3MnrWNllaltbV/mCytim094LNx0DhIqlFOChqccGqRy+X4wmcu5IzXn8Osth7232MlC9s3kQ9NqpYG6XcSo6Z2KGaRlWM1BEe6eKOiDcjGLKiy46yWGLRiJfuinp9N4UThylA1FzHjwyesJLRM0/PiWNNxbmzgicKGBabcWIanP2WmkREnLYe0k3eE5Dhnkg2CMrujk5mzl7HnHnuTnOwjkRgqfpECKJUmKfNWqWvr6ESg7onQ4aknPpmvf+2LvOCF/4/21phlc1Yyv20LKprofuoZUBQSv1bgG3W8he6CQCT0DBb2/KS4sD+JkWRLD8kPtScmJhPyzkl2aXZfBp44DClXTH50XKC1I7REJHA/HpbwYhIiNlZNtb8BhgvOnNHNfotXUsiXedObXkkuVySfT/NRh0nA4r0QY41tKhbG/nDGi4Zhpn5QKBR48YuezeWXfo7/fMmrmNnczQGzVtLR1Gv8fk4+FMl0nI7U1AeNc3bhQUakdT7BlMtppneEEx0dt6rklgkqF0cVLuk4WSWRBhViKaQGFbdX3YleezV3nn8ftRYgUyMHJDJB3EHWl4EC7c3dLOjYxCFPOpGvf+3/2HfpXrzghafQ1DxIYA1i7kfFGcYq6UECdW1BJgWurbe/1TRunP0Jaxk7kdiliNDHaaeezMVf/QJHHHUMezWvZWnLGvKBV/IxrvjhDE3qU2JswycAyRBGtUBkXyzFG1tN3Usg2XOqrR3xdLYkTckn/sBsSWFxeyHnwxw2SUgqrCXEW21yCq3FPubP2MKyA4/mgvNfS8HWcD3iiIOJygVmz+nMnG+4vKeIiyHAGR097Lf//tU+gYnBOHRCrz/h84BDgFeJyCEVp/n9Cc/B9CesdeyEYZclQoBiscDb3nImH/7wh9lz/iwOalnJvsVVLMhtYm5+G/OCbcwMugmC2BocNSXMHGgB0zseTQwzItWJAUYnuGRZVDlYGcTtn5L5oagiaiWc2v4KCLZ8ReCda8VMn+jcIWcgUktFAuSlxF4z17L3rCfYe+kBnP+W12XuGQQBZ519Jh0zu2hr7/XEgOzkc7mY+Qs3EQTwpjdVtnmYICg2dKliGxtJf0JbkNf1J/RxOrY/oXW4zxSRPWocO2Goy4iZ7cXiRXvw+U9+iF9e9zfu//d9PPLIo2jUb6I9ygMExHTFLfTSRG/QxNawDcXkyzl/nEDa4dfHSLocznDinVthcKmkuNSA5O2WdFxVTmk587DwOfUuX0XHdHNzVbgNq42Z17KVuc1bKbbM4ew3/DdHHn7g8JsCJ554NDff9FRU/0nHzG7KZSPvDwwUKA3laW7pZ0Z7L4VCC/99znm0tDRXvc6OQqhd/KzAjvQnrGXshGFaEKHD6c99Oqc/9+mZfXEcc8ONt/Lz664n7u9moGstpVKOR/ML6ZKWVE/DK1kBCWepYDYJ3LquZHxGr8tG8SQDMGJi7HpeVKH3zPXs3KJYs76xarqn3Wey5tP5Gx1WmZnrZo+WTYRhjhf/56s47TknZdwf1fC2t7+OlSufxWWX/ZjWljxPPLGFctRJFG2jpXkOy5YdyZvf8spJC1dLMCxSARi7SeioAssY59QydsIwrYiwGoIg4JknHsszTzwWgJ7ePj53yeUUH7yHRwp7sCnXkTnfp5tR173jXrgCueaMqqUX8bigpu/9xi1o6nPUyht7CqgklYS9Z6ycpyanU5RBFjVvoCkY4NDDj+ct57w60f1qwZIli/j4xy+o+fwJhxNHh2OsJqE70p+wUMPYCUNNOqHt5X2ViPxbRO4XkRMqjh8kIstFZFBELqwyPhSRO0Tk1xM18fGirbWFj7zjPJ7+nOezbGgd7VEvUEF4NfzmeUbGpDx+5XHxXlRjOrEvn2Ist8Nu7ZexrxZs7XSkil+OkIg9iuvZv3Uley6cz0c/+lHe/j+v3y4CrBeM0zqa9CcUkQKmP+E1FedcA7zOWkmPx/YnrHHshKFWTvgl4DpVfamdVEvF8S3AW4EXjzD+fOB+YBIKU44Pr3vJi7jzwRUc+Ni/ubdpb/qDpqwYB8MWfKVICp5YKSkHCyrO8wnWcbzMOrKJCa6jr3gE5et8vnuucq6KG6fMyW9jj+ImIoqcefabOOHYI2v8VOoU49AJd6Q/4UhjJ+JRqmFMIhSRduAk4PV2gkPYmhoOqroB2CAiz68yfjHwfODjwBTKNcPx+fe8lQ98+ovEjz3Cv4pLKAX5Ub/vauJpRsz0RNQMEuWswrjiXSxJ0nWEbEVVxdR+UcdtPQOODaQx146V1qCPxfn1FMKYp59yKq99xYts6tGujBEsVrWMHGd/wpHGThZq+Yb2xTRB/LYVKb8pIq3bcY//A97FGIZlETlHRG4VkVs3bty4HZcfP0SEj77zfKJCC/sNrU17OFS4CsT7O1KsqEPghZc50dZVxnbIrCnP7O475t3iSxt2prsTPUmNzSKnQywtrmW/4iraZi/k85/9FGe86sXTgACxlimt22prE4FavqUccDRwsaoeBfQCNUUQiMgLgA1VSoQPg6pepqrHqOox8+bNq+XyE4IwDPnMB99Dqw6xtLShqr7hMblhMZu+Cw2qJNq6a4wiYibvvc1xxdB53SOzz6VfCRCKsiDcxMFNK5jXUeBtF7yLr3/m/XS0b1dnrrrHeCNmdhXUQoSrgdWqepN9fxWGKGvBU4EXichjGIfnM0Xk+9s9y0nGwnlzeM2ZZzE/7ma/8hMjEmIlEn3MWjp9IqqmP/oWUMfNfBHUv6irtK2xZrik2oiRjqCbgworWFDs4dDjnsVXPvcxDjlo2biev+4xzbMoxtQJVXWdiKwSkQNV9QFMRan7arm4qr4XeC+AiJwCXKiqk99feRx45nFHs9f8uXzi059mj2gzT+SyvdhF0r7wDuq9cGXpk5ZnzsBS4XVSS3yOg2boz9P5nPHFNXRx5xUZZHF+A23Sz8GHH8ubz3oNTU3FCfsc6g7VxI9phlqto+cBV1rL6KPAmb6VyVYgvhVj/YxF5G3AIaraNQlznjTsv8/enP6fL+E3P/sJW4J2BoLUnO9VQ0zDxKoRpmpG3MR/rVZnJDXQSIDJbidr8TSGHklMpaFG7BFuZl6wlbZ5iznvnDewz5JFE/r89Ynpx/kqURMRquqdQKVj1LcyrcM4NEe7xg3ADds1uynA6c95Br/90184qHMV9xf3ZlDygBdN460HP5h5WESLxwED9Zz0kWZcDFpRE9QQsHqErszSThblN6FBnv86442cdMJoPuppiGlOhNPAfDaxEBG+8rH3M2PWbA4aepxAK4y61ayjVq8LvNf+6apeBUKya8p5GIJKh75CG/0cFKxkr9wmTjz5mXzjq5/fDQkQiOLh2zRCgwiroFAo8NkPv4+cCAujLaOe69ONVwjc7LcVsDN0W2kFhUyWBECBMvvIWg4MHqdj9gI++6lPcPZrXjr5MZp1CSsuVG7TCLvjt1oTmpubOPV5p/H731zD2nC26fxUzcleZayvOjoJNhCIo4oTfUMMJmNgAVvYI9hEXGjjLW+6gCMOrZ7lsFthmoujDSIcBS969in8+Te/oEiZISkmFlDfKpoQWLV14tkU4jgrbmZFUqWDXvaS9YQCz33h/+Mlz3/2mFkOuwUa1tHdG8Vi0RAZUZJ3mNhbPDbnV31z+p/LvvGtqsOyI4AmBlkSbKSNXpYdejTnnf1ftLZUhubu5mhwwt0XQRAwGBZpj/voDW3SqtPn8HL3HGKSkvMmxM34DP2UJVNJQxAtsyebmc8WWucs4m3nXsDSvXcHl8M40CDC3RtPOe54bl6+nLXhbNC0AQsRSS5hFKexnZK4Fkic9pk+9arMopO9ZAMSFnjGqS/l1adPUkej6QBViCqV6emFBhGOgbaWFvJqyv/5GpozK8cuu6EiCmYYhxRo0X6WyAaaGGSvg5/M+//n9bupxXM70eCEuy/6Bwb4y5/+wBP5OQxrPujDhahRfb3kKLOXbmSOdJKfvRcfufAtzJsza9LmPe3QIMLdFx+/+NsosC43Ky0oo6nvTyEJxM4gMdQoC9jKnrKJfEs7Z595Pk8+fNIq501TaMM6ujtj9UP3syY3B9Ugk+3uXBOZeFK840A7PeytG8hLzH8870W84vTnNlwO44GCTjPnfCUaRDgCuru7KUaDdOdaTOtpAIXQ+CwyUS8+aRUZYi82GL/fQYfz1rNfy8xplt+30zHNwtQq0SDCEZAURIpImrP4zvdAhDhKraKBxuzBZhayhbBtLue84RyObkS77DhcPtc0RoMIR0CxWKRESDEuMUCar2cSeE0zlkBAY2U2XezFRlRCTviPF/GGl542hTOfhmgYZnZfaK5IMS4Niw91mfTNcT9LWE8zgxxx3NM497UvJ5/PT8lcpzO0wQl3YwQhQRxnk+MF8lpmUbSJebqN/OzFfOiCN7Ng3pwpnOh0RiOpd7dGVI5Q2yTQEKKyIN7KIt2EFFo56w3/w3FHPWmqpzm90Qjg3r2RC5Q4EiSAGXEvS+INFCXi6f/xXF77khdOj5KCdQ4FtBG2tvtC45g8ZfYrr2EW3SxcdigXnPt6ZrbXTSHx6Q/VaZfEW4kGEY6CMAjYq7yJoeJM3nDW+Rx/RCPaZSqgDXF098WJTz+Z+XNn8/xnjd1GrIFJxDTnhFJZsq8eICIbgZVTPY9RMBfYNNWT2E7sKnNeoqpJCXYRuQ4z90psUtXn7rxpTR7qkgjrHSJy6xi98eoOu+Kcdxc0zHsNNDDFaBBhAw1MMRpEOD5cNvYpdYddcc67BRo6YQMNTDEanLCBBqYYDSJsoIEpxm5FhCLSJCI3i8hdInKviFxUcfxCEVERmevte6+IPCwiD4jIqd7+J4vIv+yxL4v15otIUUR+bPffJCJLvTFniMhDdjtjR+csIufZed0rIp+plzk3sJ1Q1d1mw2QktdnXeeAm4Hj7fi/gd5gggbl23yHAXUAR2Ad4BAjtsZuBE+w1fws8z+5/M3CJff1K4Mf29WxMb8fZwCz7etZ45ww8A/gjULTH5tfLnBvb9m27FSdUgx77Nm83Z5n6IvAu7z3A6cCPVHVQVVcADwPHisgeQLuqLlezWr8LvNgbc4V9fRXwLMtxTgX+oKpbVHUr8AdgzIiPUeb8JuBTqjpoz9tQL3NuYPuwWxEhgIiEInInsAGzwG4SkRcBa1T1rorTFwGrvPer7b5F9nXl/swYVS0DncCcUa41rjkDBwBPt+LjX0TkKfU05wZqx24XwK2qEXCkiMwEfi4ihwPvB55T5fSROp+N1RFte8eMiipzPgzz3c3CiKZPAX4iIvvWy5wbqB27HSd0UNVtmPbdp2N0p7tE5DFM2+/bRWQh5pd/L2/YYuAJu39xlf34Y0QkB3QAW0a51njm/Fx7vZ9ZcfVmTAniufU25wZqwFQrpTtzA+YBM+3rZuBvwAsqznmM1DBzKFkjx6OkRo5bMFzIGTlOs/vfQtbI8RP7ejawAsO9ZtnXs8c7Z+Bc4KN2/wEYsVHqYc6NbTvX5VRPYKc+LBwO3AHcDdwDfKjKOQkR2vfvx1gYH8BaE+3+Y+w1HgG+Shp91AT8FGMQuRnY1xvzBrv/YeDMHZkzUAC+b/fdDjyzXubc2LZva4StNdDAFGO31QkbaKBe0CDCBhqYYjSIsIEGphgNImyggSlGgwgbaGCK0SDCBhqYYjSIsIEGphj/H0lPWDx9BGllAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "da_diff = da1-da2\n", "da_diff.plot(title=\"Difference\");" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAANoAAAEICAYAAAA5nOpLAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABuwElEQVR4nO2dd4AkVdX2f+dWd0/YmdnZ3dkcWJYMkhEQFcQsIqCivkYEkdcMmBP6YQbF9JKVYEARwRxQVDARlCxI3sDmOLs7eaarzvfHvbeqOsxM7+zOdO/SD9ROd1XdW7e66/Q55znnniuqSh111DG+MNUeQB11PBNQF7Q66pgA1AWtjjomAHVBq6OOCUBd0OqoYwJQF7Q66pgA1AWtjjomAHVB2wEQkWtF5AsTcJ2lIvLi8b5OHTsedUGbQIjI60XkdhHpFZHbhjnnTSLyo+28zgwR+bGIrBKRLSLyTxE5anv6rGP7kKn2AJ5h2AR8E9gXeOEw55wA/G47r9MC/Bv4ILAOeAfwWxFZqKrd29n3Mw4i8sEKTutR1SuGO1jXaGOAiBwqIveKSJeI/ARodPuniMhvRGS9iHS61/N8O1X9k6reAKwapl8DvAS42b1/q4gsE5GNIvKponOPFJE7RGSziKwWkYtFJOeus1hVv66qq1U1VNUrgRywz3h8Hs8AfAT749U6wvahEXtQ1ZrcgKuxv8YPVXj+64H/Ag8DPxrHceWAZcC5QBY4FRgCvgBMA14LNLsP/6fAL8r0cSZwW5n9RwN3uNf7A93AsUAD8HUgD7zYHT/cnZ8BFgKPAOcMM+ZDgH5gcrW/151xAy7c3nOqfhMjDPxY4LBKBA3YC7gPmOLezxjnca0CJLXvduALZc49BOgss384Qfs8cJ57/Rng+tSxScCgF7Qybc8Bfl5mfxvwH+AT1f5On8lbzfpoqvo3EVmY3iciewCXANOBXuCdqvoo8E7gElXtdG3XjePQ5gAr1T3FDsvc+JqBbwAvB6a4Y60iEqhqWEHfJwBnpa6z3B9Q1R4R2ejfi8jeWC13BFaDZoB70p2JSBPwa+BOVf1yxXdYRwFG89FU9euj9bGz+WhXAu9X1cOBDwOXuv17A3s7du1OEXn5OI5hNTBXRCS1b4H7+yGsH3SUqrZhtR9A+tyyEJFZwGzg3tR15qeON2NNU4/LgEeBvdy1Ppm+jog0AL8AVgL/W+G91VEe3g87Ang3MNdt78Ka+KOiZjVaMUSkBTgG+GnqGW9wfzNY8/EFwDzg7yLyLFXdPA5DuQPrK31ARC4BTgKOBG7Ffhl9wGYRmQp8tugeAqxflwGMiDQCoaoOYbXZzSlNeSNwl4g8D/gX8DkKfxhbga1At4jsi30A1rvrZF37PuBtqhrt2I/gmQVVPR9ARP4IHKaqXe79/8P64aNiZ9JoBtisqoektv3csRXAL1V1SFWXAI9hBW+HQ1UHgdcAbwc6gTcAP3OHvwk0ARuAO3HsYQpvxT78lwHPd6+/444V0Pqq+jDwXuBHWO3Wib1Pjw8DbwK6XB8/SR07BjgReClW6Lvd9vyx3XUdDguwfrLHIJaIGhVS6GrUFpyP9htVfZZ7fzvwDVX9qTPdDlLVB5yp+EZVPU1EOrDEyCGqunHYzmsIIpIB1gB7qOqWao+njvJwIZbXAz8HFHg1cIOqfmnUtrUqaCLyY6wp2AGsxZphf8Fqg9lYE+x6Vf2cE7qLsCRECHxRVa+vxrjHAhGZAbxWVS+r9ljqGBkichjWGgH4m6reV1G7WhW0OuqoRTifeS9VvUZEpgMtzl0ZuV1d0OqoozKIyGexzOM+qrq3iMwBfqqqzx2tbU2yjh0dHbpw4cJqD2PC0dXbAwitzc3VHkrVcM8992xQ1en+/cuOn6QbN5WGIO95cOAPqjqeYZxyeDVwKC4Eo6qrRKS1koY1KWgLFy7k7rvvrvYwJgSqyjd/9zPu++s/adrSw1A24NKLvkFjQ8PojXdBiMiy9PsNm0Lu+sO8kvOys5/qmLBBJRhUVRURBRCRSZU23Jno/V0KURTxjV/dwDs++kEe/s0trG0O+Nth88mLcNFPfljt4dUQlFCjkq0SiMjVIrJORB4a5vgL3DSi+932mVG6vEFErgDaReSdwJ9IwjMjoiY12q6MfD7Pxb+6iXv+9g+yg0M8PbuNZft2MJQJIBKemD+FzJ1303nKqUxpm1zt4VYdCkSMmUe4FrgY+P4I5/xdVU+saCyqXxORl2ATBfYBPqOqt1TSti5oE4SBgQGu/u0v+fff/o4O5Vk5ZzJPz24lnwlIP0erprWycPUWPnTJN7n6E58dvsNnECLGlthSLl92e+EEqyLhSqMuaOOMnt5evvPrn/HAP+4gQlk6ZzIrZrUSBsVWu9j/jfDQoukc9fAqzrvuaj7/5jOqMu5agaIMlTcVO0Qk7chfqXbe3bbiOSLyAHZGxoddRk5ZiMhrgAuAGdi8UgHU5ZqOiIoETUTage8Cz8L+/p6hqnekjr8Z+Jh72w2822VszMeq7VlAhP0wvlXJNXd2bO7q4vKf/4TH7rqHfGBYMn8yK2e0oiadXyygSpwLrPafrZMaeGj36cjf7+KHu+3OW553/MTfQI1AgbC86bhBVY/Yzu7vBXZT1W4ROQGbhD1S6t6FwKtU9ZFtvVClGu1b2ITXU90s3mL+eQlwnKp2isgrsFn2R2GTbz+kqvc6GvQeEblFVf+7rQPdWbC+cxMX3/hjVtz3H/pzGRbvPpU1M1rQcgn8xUImyd/V01tp7Rtg6Pqf8vy99mW3mbMn7iZqDNvho40IVd2aev07EblURDpUdcMwTdaORcigAkETET/d4+1uQIMUJlaiqren3t6JzaBHVVdjE2JR1S4ReQQ7vWCXE7QV69bw5euuYeDxZfQ0ZVm8dwdrpzYzIrErQskzJCBqdz8+fyqtvYN8+msXctUXLyCXy43jHdQmFAjHKanCTU1a6yj7I7Ff1kj5sXe70hW/AAbiMar+bNgWDpVotEXY6RfXiMjB2MmFZ6tqzzDnvwP4ffFO55QeCtxVrpGInIWb9LhgwYJyp9Qsbrv/bq67/Dv0NGV5Yr/pbJjSZIVoWB++jHbz2sw9UyIKIjy41wyOeXAlF3z/as47813jcwM1jrHO8Unny4rICmy+bBZAVS/HlqF4t4jksTMp/qdoQm8x2rATjl+a2qckszeGRSWClsGWFHi/qt4lIt8CPg6cV3yiiByPFbTnFe1vAW7C1rTYWtwOwDmyVwIcccQRO1Ve2MF77M31c2bQsnIduy/fwsz1vXS2NbBuajND2QzWF6NUvrTMayncPxRkeGiP6TT8+z5+dsDfec1znlkzXRQdzkcbva3qG0c5fjGW/q+0v9PHNBAqE7QVwApV9ZroRqygFUBEDsISJq9IT09xkxBvAq6rRMXujJjS2sZl532OB554jGtvu5ktGzuZtnwD+y7exMqZrSyd005/Q6ZQsIabc+3OUU1O2NjazLKZk/nNj3/CC/Y/iKmTnznxNVUYqvLProh8VFUvFJH/o9TYR1U/MFofowqaqq4RkeUiso+qPga8iCIfS0QWYNXnW1X18dR+Aa4CHqmkrsLODBHhkL335Zt77wvYzI+f/PXP/P7m3zP/3qdZOqedxfPay9D6pfBCZrljQOHxeR20d6/knAu/xLWf/zLGPFOSeoRw9EoQ4w3/vI85L7DSb+v9wHUi8iC2stOXRORdIuKdhs9g61lc6lJZ/ICei51V/MJUmssJYx3szgRjDG88/iV8/4Kvc/Jb38KsTT08777lzNrQDRpR4JClfiNVkx3xaQhq4IFFMwk2b+VrP7x2Au+kulAg0tJtgvEG97ddVb9XvFXSQUX0vqrej50ekMblqeNnYkuoFbf7BxUUptnVceJzn89Ljzyaj1x9KQfc/whTtrbwyKJpjnUsct4E0GI2UkGF/oYs/1k4k8w/7+KJF7yYvXYy0misqAGNdriI7AacISLfp+iZVtVNo3XwTLE/qo5cNsu3/vds3nPO2cze0MM+SzdZIfP5BR4pjRbvTwndhvZJbGpr5pLrnxmJxzZgLSXbBONybP2XfbGse3qryJysC9oE4/B99uNVp72FeWu72evpTrvTC1I6aA0FxyT1/sk5U+l5ahm3P/jARAy5qlBgSE3JNqFjUP22KwR1taouUtXdU9uiSvqoC1oV8Oojn8vbzjqT3VZtZf5qF+1IC5sPB6T2aer9luZGlk+fzJVXXUV3z3DhzF0DihBiSraqjEX13SLyPBE5HUBEOkRk90ra1gWtSjjukMN54amnsO+SjUzd3Of2ivXP4gitWD8OQdRmjHhhe3xuB1EUcsWvfzHRQ59wRColWzXgShl8DPiE25UDKrLh64JWRbzpRS9n0sH7ccjja2nuHURDCjVbRJm0CCt4kRhWTm3j4XvuLT5hl0KN+Gger8YWzO0BW8oAW8h2VNQFrcr4xrvOJj9rGkc9vIppW/oQdRoMiXmS5LEqfMDWtLeQ3dpNd/euvOSZEKop2aqEQZeiVS9lsLNBRLj2059jn6OezeFPrGLWxq5Ek5XmICR/jdpsE2DtplHZ5Z0WVrGbkq1KqJcy2JlhjOHjp53BtdM64Ne/JZNXVkyfTIGkxWxkEhIIjZA3wt1PPc4eu2hMTVUY1KDawwDqpQx2Gbz9xJNY3tuN/uVvtPb188TcaeSz9iHzZGT8RwEjrG1v4R933M4bjt9115CPqh+wjjHWUgZ107HGcN7r38Tb3vsuZvQNcexDy5jR6f0vTc0RVcf+Kys7Wul/eiVDQ0PVGvK4wpIh1aX3xS6hvHW4rZI+6oJWgzj2wEP41DnnEkQRJoqcWIEVNrXfmlHIQF9TBqMwODg4Upc7MapPhqhqq6sL8k3szJW52MnNH8MuqTwq6oJWg+gfGODz3/4G69ubWd0xCQK1m8FuKUsq72YDrFy/vipjHW/UGBnyMlW9VFW7VHWrW5TktZU0rAtaDeKcr1+ADuV5aM/pELjqWJKegq1IYLXbUNbQ2dLANX/8dZVHPX4IVUq2ag1FRN4sIoGIGFeUqpIlk+uCVmtYu2E9unQlD+49nXzGID45BLWvDYjxaSIWK2e0svo/j+yS5qMiDGmmZKsS3oRdH22t217n9o2KuqDVGC767U30NGbY2pa1fpiJ7F8pMh3F1hWRQFk9cxKo8tVf3lDl0e941AIZEo9FdamqnqyqHao6XVVPUdWllbSt0/s1hIHBQdbf+xArZ7ciRojLxKhCQDKtxsOt5R2JYWtLA92bN0/wiMcfSlVNxQK49dDeiV1ON5YdVR21ym1d0GoEURRxzsUXEmrE03Nt+pyfFyrp0o9Om9nivcmkUVHoDXdNir+K5Ecxfgn8HZsRUpFv5lEXtBpAb38f773o87C6k/8c1EHktRexGKX+tRDjE0Vs1nF/Y0Cmb4BdDapUM7exGM2q+rHRTytFzdzBMxW3P/IA7//0xxjasIV7jphBV3uOIFDrf4mm0kBA0/NkHDliMiABbG3LsfWxJbz7C+dx5Z9/U6W72fGwZEhQslUJvxlrzZu6Rqsi7vrvA1x98WWsnd3ME3tNIwoMBkEFDIoiCIpG1mRM0h19CFtjVnLVnBZ6J+WYtbaboRt/zX/+fTefOOu9zJk6faQh7BSoFvlRBmcDnxSRAWCIbVjkombu4JmG/yx7iisvvZxV8yfx+P7tRJlkekw6yQoUE9iwbXxMUvE036EIm9sbeGyfqfzryFms37iRT33li3T39lbj9nYYlNJJn9Wa+OkyRIyqNqlqWypjZFTUBa1KuPyH17KlPctTe7fFZAei8VpgxggZI8TlGwUi538bAWN86f70TFErhL2Tstx36HQGwzzv+cpniaKxFtWuDdQKvS8ix5bbKmlbF7Qq4M8P/Juh5etYvPfkmFIUFxdL6qJaATKkjhtAEl3n6xuIRPavUQJj+wkzwoOHTMN0dvPZK3belbJsXUdTslUJH0lt5wG/Bv5fJQ3rglYFXP+zn7J2dhN9LVkEF4NOZeanX0dFJqIVSDszNJBEAI0BEwuefT/QlOHBQztY88Cj3PvYzrqAT2kZg0pLGYy2hnXqvGeLSCgip450nqq+KrW9BLte4NpKxlIXtAnGTbf/BdZuYemerYWFryAmNsQJm9ddxSaiCGQDQT0zGU+aKUV3W44t7Q389vbbxvfGxgm23NyYWcdrgZePdIKIBNhVPP8whuGtwArbqKizjhOM3996C+vmNjPYnIGCFYKSX2l1RfcDkYQSUSuARqyApTUfmg5dJ7353jd2NPLk40+M302NI1RlzKZihWtYvx+7CMuzR+uvaJELgy2PX1FxzbqgTSAGBgZgZScbDp0CpCqCx3xiIdJRM++7JcR+Spgc9W+QeHXMdF+9zRmCp3feAj7DBKy3ew1rEZmLrWz1QioQNAqrEueBH6vqPyu5Vl3QJhA/++efiQxsntow4nlxmlVK0tJCNByGOz7YFJAdCMnn82QyO9dX7gMbZbAj1rD+JvAxVQ1FRvf7Kl3QohzqPtoE4m///AfrZjcTVfSlOlPQnVogRF4QSYQyjeJdebdA/c5Z7mBcZ1gfAVwvIkuxq39eKiKnlIxAZFRNOdo5O9fP206MJ5YvI1yxkTXH2EyNxD1LRMJQWC9VnYFoZHhtVs7LKz4zCuyR/v5+mpqaxnYDVYKl98cnQK2qcTlvEbkW+I2q/qLMqaeISP8IXQlw/EjXqgvaBOHyG39IZ0cDPW3ZwpJWDsU+Wpy5T3mTUHCkSdH55RA5jTYwtPNNDPW5jmNBBWtYV4qPVHDO30c6WBe0CcCTK5+m67GnWXZkB5piCIvJi2I5CUQIhxGeEuJECzVaWnBDp9FWdm5g1oyZ23k3E4+xTpMZbQ3ronPfPsKxMftmHnUfbQJw5U0/YvPUBra2N7jQs5uw6dKGi+H3hCkVZcqcF+/T5E9xbA4gCgydUxv4xV+2uRxh1WGnydRMzZAxoy5o44zBwUE6/7uEp3dvKfKnxDFqWiAgpTS/fTca4zgaVsxvYdV/HqW7d+db5qlWkoq3B3VBG2d09/RgFHqbs36OJgaJtVVaU6mm/Kz4YRJ3Xim88KX7KPsIKqzvaCSfNXzv5l+O/WaqAJu9XzO5jmNG3UcbZ6zq3ABAPok4E0mSuxj7aJpoM1GvxxKxiSin7Urja8MeF+HpBa0Ef/k7f9tzH4496PAdcXvjDpuCVRuCJSJ7Y4mR3SisGfLC0drWBW2ccdVNP6FzSgNhRgjEEGpUGJCGwmRH/zadU+V9MPdeUudUHsQWVs5vITcY8oPLrmTTaW/ilKOP2+77G3+MPQVrHPBT7HrW32Eba4ZUdAci0i4iN4rIoyLyiIg8p+j4m0XkQbfdLiIHp469XEQeE5EnReTj2zK4nR1/fuDfdD+1nMf3aUdFiDTRXj71Kr35TBBJC2Bajsq8LjYV0+/9dXyuZISweM92ls9v5dc/vWmnqQMZISVblZBX1ctU9V+qeo/fKmlY6U/Ft4CbVXVf4GDgkaLjS4DjVPUg4PPAlRBnRl8CvALYH3ijiOxf4TV3ejz85ON0t+boackRRRCpTZIdjdcodzjNMKbjZeXyI6GcICf7lyxqIz+U59PXXbEtt1MV1Bjr+GsReY+IzBaRqX6rpOGopqOItAHHAm8HUNVBoOCnUFVvT729E7sAAMCRwJOqutj1dT1wMrCzTo7aJqxcv5b+hiDlczmoQZ1TNtoj44UqdNn6BiHTn2dSzxBzl3ezcn4LnVMbkgltFApiuWyR0BiWLGojc+8j9L9pgMaGkXMvq40aMh1Pc3/TAWwFFo3WsJI7WASsB64RkftE5LujLCn6DuD37vVcYHnq2Aq3rwQicpaI3C0id6/fRRZsaJ08mUk9Q85klFibxUSHpsy7tJYqCYYlQhSqcujd6znwgY1MX9/PIfdu4Oh/rmH62t6yUW+bwmRfp9nJ1bMmkRe46OfX7cA73vGosZohu5fZRhUyqEzQMsBhwGWqeih2oeyyvpaIHI8VNF/7rtwnUtZwUtUrVfUIVT1i+vSdv3ITwHtOOJWmvjztm5N6i0pSgVjVFuPxApAWOm9q+ik0ABpBpi+kuS/PeeedxxVXXMEXvvwl5h+8P/s91Ml+D21CQt9GSgS4gJ00wop5LSy++wGrXWsUCuTVlGzVgIhkReQDjq+4UUTeJyLZStpWMuIVwApVvcu9vxEreMWDOAj4LnCyqm5MtZ2fOm0esKqSge0KmDp5MoPzp7Hg6a54Xyw07tmO1CcSW42nmiZHpEDgVIUwCBjMGn78x98BMH3qNL54+vs449z3Ma0rz8H3rcfko1iLpVGsPVfPnkSmq58nVzw9Hre/w1BDcbTLgMOBS912uNs3KkYdsaquAZaLyD5u14so8rFEZAHwM+Ctqvp46tC/gb1EZHcRyQH/A/yqkoHtKvjom86gY0M/k7qtW1tOdxSX1Pf7TIqN9OxhmDH854AOnrrjbn74p9/F5x+z97O48DPn09wfMmdFT9kEY1M0p2YoZ7/+7r4aLklXxmysYmbIs1X1NFX9i9tOp7IJoxWzju8HrhORB7HTt78kIu8SkXe5458BpmHn89zvZ76qah54H7YewyPADar6cMW3tQvggN33pLejhZlregu0SbEgFGug4uC0RsmOzimNPLbXFP5606+4/vZb43M62qey4LADmbO6B4pMx2IzUhAbchDY1Jto3FqDTVOrGXo/FJE9/BsRWUSF8bSKAtaqej92klwal6eOnwmcOUzb3wG/K3fsmYJn7bsf9zz8YPy+WMh8/CwwLovDC0d8QurBcjtXzm2lqS/kDzfcxImHHEVLczMAWzZsoqvFTcVxc2dKuBVSMT0RghrPHayh3MaPALeKyGLsb+FuwOmVNKwZ3nRXxoF778vkLYM09ebxREWahbSCZAgjW/5bNS1kUBzYthCeWtTOoMDZl3wVsAnMW5csZ8O05sK2auw1CgTcjiMfCE+tq1232bKmtWE6quqfgb2AD7htH1W9deRWFnVBmwCccORzaV4wi2c9vBGJ/LxpC40FDyReabCcYKThzEIR/rvfNMxTq7js5z/hw5d8DQXWTW+mtKBBqv5BSmOun97Mv+7614694R0IRchHpmSbSIjIC93f1wCvBPYE9gBe6faNirqgTQBEhM+970NM6h1i6qZ+RJN0K1RKTDtLjniKXgrYQiHldwFbWht4amE79/zxVtasXsOj+0xDTeLhlc53S81hU8s86qoN9PePNFO/uqgBH80nhb6qzHZiJR3Uk4onCO2trUjbJJr68qiqq91YkNAxbGpVmhjxpEl639LdJrN0QRuIpPor/zAWZvsL3a0NRCLc+ch/eMGhFRFoEwutvo+mqp91Lz+nqkvSx0Rk9zJNSlDXaBOIqKWR1u4hECmg7qGUIEnnNGrRfs+XBGkp9Zn/qZONm1yafkyjeI/dG4rQ1ZLltw/cRS2ilnw0bKHVYtxYScO6RptAvPEVJ3HDd65m6W5t9Da6hAIvUEXPTrk4WAHlrzb/sdwxP7HU50d6L7B0So298FAuYHL1FvcbFdXWaCKyL3AAMLnIJ2sDGivpo67RJhAvOfwoBqa3sfuSLXZHQaDMyVqZfMVhoUn+YmIMFtYa8TuHm7cmAkMZw+L1Fa3VMOFQhDAyJdsEYx+sL9ZOoX92GHbx+FFR12gTjLe/5nX8+IqreGLRFAZziRZJB5SHq6+a1lh+qcKoKK3EK0fFrqMWn5OC99MMhkgjBhoCsv21W1y1igFqAFT1l8AvReQ5qnrHWPqoC9oE4/hDns21LT9kzupulu42ueioq8GvWkKSpN+XaKci0zPNYhpcInOqjW8fEYFAX2OGgVW1WZtfa4AMSeE+EXkv1oyMTUZVPWO0hnVBm2CICHsffCBD9z/IkgVpQUsvdCFxRn28VpqmzhhuTkQpN+ISlstmWMavhrIGGSrNJFq/eROf//63yeShsa2Zc089k+ntFc1z3KHQ2hG0HwCPAi8DPge8mdJJ0GVR99GqgDcc9xLaugdpGAhJM4CF4pBmEQWDid03bwmaYaJk5UzPwnMVYyR2CtWApJbfXbNxPZ/90SV84rxPsWn1OpblN7LysSV85qIv0DPhCcg1lVS8p6qeB/S4oqqvBA6spGFdo1UBe87fjYHmHB0b+1g5u3WYswodtqiA7tDYz3Jh7SJzsdBULOrRpUDad0ZgqCEgMxTy6JKn+OJVF9O8vpf+SQFr9m1mw8IcuWxINq/sfms3H/7W+Vz28Qt2wKdQOWpIo3lHdrOIPAtYAyyspGFdo1UBIsLcvfZkxqbekgzG1Fn4lKlydRtDlxDpSRF/wBIg5YQs1ofuP7s8LwLdrVkiEb5+wYX0Sch9x07lwRdOYd3CJkIC+oeyDAWG5Uc3Ey7ZzFMrlu7AT2NkqEIYSclWJVwpIlOw61f/CjtdrKJfnbqgVQkL99id1u7B+PGPUqxjMcpm9JNyy4qeu8K37gqp0lqBEbvmtbHvJYDlCyfx5H7tPHzoFHpbMqhfxB77w6Aa0DMpR3dHhi//8FKiKGKiMNYUrNHWsBaRk13ltvtdGY3njdLlNaraqap/VdVFqjpDVSuqcFQXtCphj45ZZMLCR8b7Vj6eVlo7BPzyuiIJfR+3oTAMl+47PQkgcrUl483A8r3aWT1/Era0q4lzML0rl48se7nywElEq7Zy7oWfobtv/MuLW59USrYKcS0jr2H9Z+BgVT0EOANbIWAkLBGRK0XkRVLJyoUp1AWtSmjINRCEShRpUhrOSUik9qE2Zb5Lb/z52v3p/f6F97+8CSmpSHiQ1lQkNEyUnA1AGGnMgmqK6exrz/DEcW1sXreRT3/3ou39GCrA2MkQVf0bsGmE492aFEyZxOipAvsAfwLeCywVkYsr0IJAXdCqhoHBAbecUjIlJl24B4ZLC3bB7WGeCcWaevGstVQnxoCKjdGlPbbC2d3+vc3AyId2TInQK0MtGZYe3kLvwyv563/Gf4pNevZCKrDf4aumue2ssfQtIq8WkUeB32K12gjj0D5VvUFVX4OtNNAG/LWS69QFrUpY37WF0KQ+fk9mYOI3ZQPTRcnIyX7XjdhkY4l1VRL8NiV60HZXQPvHGs+baYYoMgyFVtPmQyEfKr1Ts2ya18BPfvuLsdz+NmEY03GDr5rmtm1aKD7pW3/uCgOfgi3+OyJE5DgRuRS4Fxu0fn0l16nT+1VCzgQFFQo8XV9QEi4lE3GZgyKIlNo7oWrRL6hixGaCSFqwYiYlLalR0R4lUmuq5iPLf+ZVMRqxeVaW9gc7iaIIY8bnN9uyjuOvD1T1byKyh4h0qOqGcueIyBLgfuAG4COqWrGTWtdoVUL30EDB8x3FRfLLnS22zEGBWSmOqJCUFkr8M88JGpHEzyqLNOOiyauCBJS0N+eWUhTDlo4GzEDEA09WlBwxZgxjOm43RGRPT2qIyGFADtg4QpODVfXVqvrjbREyqAta1TB/xmwahkIy+TJFlOJnfvgImwKhuoUM0+SHJz5EHcPocyVL+4qNRAFEMcaneYk1E6PkzCgi3uL0rlxA9/Qcv//Hrdv1WYyGsbKObg3rO4B9RGSFiLyjqHrba4GHROR+7BoRb0iRI+UwS0T+7MMFInKQiHy6krHUTccq4dgDDuIHmYApW/pZP618hXVxgTLPBkpK4yR7C5VgKVGprgqJWNORRIbT/Rh3xPdlCu1LVK3vFxERRRrH4NbPzfHEg/8dN/NR2SY6v7DtKGtYq+oFVBhwdvgOthLWFa79gyLyI+ALozWsa7QqwRiDmT+Thau2lrWFDD5QnPCDUPiynH9Wbq1rTGrBQ/evSb2DxOT01wYKwgtpGiXSZNs0u5HMQMi//vvA6Dc9RmiZrUpoVtVimjVfScO6oFUR57/jPbT1DDJnfU9JerCtl5o4SjZAXfjgxxR+ql1Yho5UtealSZ2oYt/HswOKrh1poXZUR4oAGCMxOxpmDVtmNPDj2367jXdfIRQ0kpKtStjgCqhao0DkVGB1JQ3rpmMVMWf6dJoP2JO5S59m1YxWig3BEuNQvWZxhEdkszqSFhKzGIEkhEhhK/teNW0mlqdK8i7AZgqea9tx6H4ERGDr1CztazvH9iFUgBpKKn4vdu2/fUVkJXZdwLdU0rCu0aqMVx/9fCZ3DyBx1LjwoTI+P98THy4uJpIIn0ecUkVpDC4wpqB0XUyFOIG0ZioFx/1/kRPwZI6bxsykKvS3BkSdPYzMI4wd48U6bvs4dLGqvhiYDuyrqs9T1aWVtK1rtCrjOfsfyE8iZVLfED2TcrHeMTg/KC1QMaPoM0CSfjy1b/02GyqIcyfFmpRxEEBST6s3HWNm0l+s0JxUf2nREsJlKGsweU2133GwAl1djSYiHxxmPwCq+vXR+qgLWpXR1NREaIRsPix4sGOzr4BSTO0rctBiwUsFoVV9hr4/Ryn/yBaarHFNEk0lLmuanLH+otecJlSiYJyEQYHqm47DTRqsGHVBqwGEgSETetEq1CoFu9ybAoFxQjdctC0p+JNoMK/ZVAoFT51EpZcf86RIwnAa/M+AN0+DoYgwM37CUO11ElX1/O3to+6jVRlRFJEdChnM2opYPv/CvylICnbCk55OY7M0UmufpZwnKRLQ+JqaiHHyOlGRRVZlvOiGiUXUxIFwVcj1R2RbmrbnYxgBpYxjtVhHEdl7rAHruqBVGZu3bEGAgaw1LuIfb2fNpb0lP4/Mw5/riRKJRc+ykb6btPBqSkMlus3CU/Z+pdECH01TsTgtTAdr7AnpbxrHR6l2AmnfAT6BK2mgqg9iF9ccFXVBqzLaWlvJB4aW3sFCGynth5H2v4ZHMqNM7MKFBV3Z3EhP1ac5Tt+trfVYynqms0O8dotnBIgwZf0QzzuoePm8HQRleyZ+7mjUA9Y7KzKZDGZ2B9M7e0vyp8o9TiWLGLp/fBKwMy7xJqD3uyItHytTbPlw32/kV9ROUfn+OrGJGZuWwtSV/TR25znhqGO35ba3DbWj0eoB650ZzdOn0rxkmXtXHLBOO0ziNFJKaMqYkjFtISnhcKalr2wsBecmrKPPiSxFYoAqkEGY+3gXc5/o4dBXHs+8mXO26Z63DVVnHT3KBazfXEnDuqDVAOa1TmFVtLRoryYyIIkoRMUkR5GJGUMs0eGFzWBcCK5QHSSTQg1KlFptJqH3CylOi+lLe5jzZA8nn/EmTjzyuDHdd8WYuDpAI0JVFwMvFpFJgFHVihf/rgtaDeCBp55gKGOLd1toAeNoV/ckPl5cInwkeLGxZElRWMBdB9JkSqG2i7xJ6jUqQuPWAXZ/uIuX/c+rx1/IaiOOVoBtnYsGFfpoItIuIjeKyKMi8oiIPKfo+L4icoeIDIjIh4uOnSsiD4vIQyLyYxGpaJmbZwruefwRmlduYPG89oL96ZnQvrhOQeWq9LkF56eEteRVGcO0xN/RVJwuJXKxD6fMXNHH0KxWTn3BSAWmdhxqJQVre1ApGfIt4GZXW+FgSuuNb8Iunv219E4Rmev2H6GqzwICKqRDnym45IYfs3baJLa2NiQ7i5jBpCZjgnLCk/hjAlo8H8BTJKmLOCpTtIzgUtze9RJFdKzu54VHV1T8acegdsiQMWNUQRORNuBY4CoAVR1U1c3pc1R1nar+m6RkchoZoElEMkAzsGp7B72rYHNXF7JyHStntrg96vIQt6GTlHlZDE/Dx/mJqXh2TFWibkpOocCGmhSgi/tTRVQYyhluu/kWrr/9loLjy9es5O+P3rcNg68QKqVbFSAizSJynoh8x73fS0R22BrWi4D1wDUicjBwD3B2JXaqqq4Uka8BTwN9wB9V9Y/D3MRZwFkACxYsqGTsOz0+c9136W/IsKHdZ1U4dtCk3hblbYwOz1Jq8laS6Hfs20lhlkkCT7vEeSEFx9TAA8d0sPDxLm657iZe9qwj+f6ff8V9f7+Dhh6br9l1+hs54egXVDzi0VBGoVcL12Cff+86rQB+CvxmtIaVmI4Z7MqGl6nqoUAP8PFKRuXqlJ8M7A7MASaJSNn5O6p6pS8dNn369Eq63+mxYdVq1k9tLmU2Yhsv8cuGYxYTpOyqAso/0VXxK+fPFRdoLTAfi8iXeFgCGghL9m1hsMHwkU9+nAf+8k9WLWzi3hdN4+n9W/jZD67nyTXLK/wURoEKRGW26mAPVb2QJDOkjwp/AysRtBXAClX1q4nfiBW8SvBiYImqrlfVIeBnwDEVtt3lkVMhSk9NKXCWhvkZl9Tmd6VMwwItKC5LxEtX6rBd4E+3iVhIzE5ADEsPaGXNgkbue8E01uzezGBzwMbZDQR5pa+/r/KOK7pw0VYdDIpIkx+BC14PVNJwVNNRVdeIyHIR2UdVHwNehF1FoxI8DRwtIs1Y0/FFwN0Vtt3lIZOaaOxJQjEyBhspyVlM2sa7CrysVN+SvEhf0xqMRXE2iuZmpwS8s6OBLdMbCoS7Y2U/0dQmDly49zbfy7CoHdPxs8DNwHwRuQ54LvD2ShpWGkd7P3CdiOSAxcDpvmSXql4uIrOwAtQGRCJyDrC/qt4lIjdiq7rmgfuwkfU6gL6+PoYyhtzgEPPXdJMbyrN6VgtbJzsGsiKjpMgnK2oXT3HR4jaFwlnYuNDchCR2l47hpRRljOauPPn+PPc+/jCH7X0AAJ1bN3PTnX/izce+kqbGMWT514igqeotInIvcDT2zs8erthqMSoSNFW9HyjOGr08dXwNMG+Ytp/F/hLUUYSmzT30TcpxzP2riBpzBK2TmHf/Gu46fBbdLQ0l5ycz0YpIjDJmZPGkz4QE8XtKn14p86rknDKH0kL41IFtLPxvF5df9G2iKU1ETQGZ1XZ97EcefIivfvAz21aWroYC1iLyauAvqvpb975dRE5R1V+M1raeVFxFzNl7T9oHQ9r3XsS1F36Daz77JfqnNDN3XQ8SFDsjiqap+mGQDmqXh8bB7fLn2OsO57sVL8KRvLedRRlhycFtPPiCKayYa9jYFPLo89p4+EXtdC5bw0ev+Vpxl6NCtHSrEj6rqlv8GxfmqkiJ1FOwqogvvvecgve/+uetNHX2smaPGXiBSBPtJcx7gV9mwwKF883Km5OF5b7Lo3TSqJYxNiXurDgtrK81w0Bbxu2LEBGWHN5K5s6n+PdzH+TZ+x40yghSqBHTkfKKqSIZqmu0GsA9Sx7lzM9+hN/84Hqe2HsyXSkfrUSDFf2ci4AxiglADASB/8l3zGOyOE0Byj276cySqKBiFm6fOy8VoivMN7Gv43QtxdUYsce3zsyyYWEDl155OVu6to7+waRuuUY02t0i8nW3GMYiEfkGNq42KuqCVgP4yZ9/R19PL/cePZ2Vu7VgTOQEzBMShcITP+JFWso/4X4VTzGp92W2cig0WNMcpOMki4StRHP6dQCcX5WqFIUAK541iaEsfPKbX2JoqFwiUblB1UZmCJYUHAR+gg1U92OnzoyKuqDVAPp7euhpzdLdmsGkhMX/TW+mSFCMFGUlFj/4Ra+dzmE4eyytxbyGi1L0YmFyrxBGEEUp7VcUXbc1hxJNlzfCk0e10rthC+87/2M8sWLJyB9OuRhahRqtgjWs3+zWsH5QRG53mU/DD0W1R1U/7hIrDlfVT1SayV8XtBpAZJT2jQMcdsd6jvvTKp536yqaesr/2hc/Y+lFcf0eL0yGCCOKkdIl1tPCWaylCvdp8t4JXGEpASv18WqlKZMzLtjqhDFSQUQYasrw3+Pb6cwOceGXvsLfH/z3yB/QGAWN0dewXgIcp6oHYRchHDH0JCLTReSrIvI7EfmL3yoZSF3QqowwDOl9chUNAxFNbU0c/tIXYCKlrWuI0Wh2b1r6ZZfEuDln/phIgYloz4/ifiXVV5I9krZHk30ai2pKqFItfAFWjftMUlWiIh2qCkM5YfGRrazbo4nvffdqHl2+ePh7jUq3SlDBGta3q2qne3snw4SoUrgOeBSbUng+sBQY5VfCoi5oVcYvbv8dmQFFA5g9bSr3/PE2gghmr+ohMzT8E5X2kRRHiIh/wBUjSmASX86YtMkZlfYVb1qouZxGKpwLZuN5vlKWJ0miKOkvihLhtcSkFbq4fKUTvlX7NdE1NcMFV35r+A+pvEbbIWtYp/AO4PejnDNNVa8ChlT1r6p6BjZ4PSrq9H6V8cfb/oAohI3CvLl7sP7BZQBM3TTAjLU9rJo3cpHc2C1K7UnSqnyZ7siZe4VrX3s+MJ0LiZ9LE0OSQ+pied5fIzEtvTYLI3tNIzaX0v+S+2Kt4mIBfoaCirDiwGb2v2Uzf77v9pL7G4Fl3KCqO6T0logcjxW00SbZeXt+tYi8EjvlazQtCNQ1WlXxizt+Civsd/eCE4/jXae8gS3tjXFy+tTOkfNVvcaxryGREEGMYoySCSJymUQ4fG6kPz2M0v25v6nXXqOV5FvFZIjBJyl7vwwoIFD8wvNR5Au2SnyuAv3NAZvn5vjxn389zI2OH+soIgcB3wVOVtWRltUF+IKITAY+BHzYtTu3kuvUNVoVoKpc+YvPce/NKwEhM1958/F24nkweRKyuZ9oZjstXd1FNHpJT9Ysi000K1w280MxxpXuVsEEShTZashRlCoqrkl8zI4tuV4xOSLizEhJBNabkRKfZAUoqXditZdfyFDVTjQ1Xp06DTfUaGjV3DAfWGWf67ZCRBZgZ5S8VVUfH+G8C1T1Y0CTywzZAhy/Ldeqa7QJhqryte++l3tvXgUIQVvIuW//YBxvOus1bwDArN3MpN48Tb15144i7QVpyh2sprAhgAhjIoxAYCATWM2WMIiBFZgo0UKxaRkZNDJEoRRsXmNZ08+UsJbWXzNEvsqxOqFSIQxT2s5fRxwvqgZFCDNCZ/eWsp/ZWAPWFaxh/RlgGnCpiNwvIsPNLDlBRLLYKsVjQl2jTTBu/O3nefJuu0D81D2ED/7vF5g+eVZ8/Lh9D+bbU5vo2GTncw1lAnfEreSSMu/8C5PxZxCfG0WGTBASDEU0bQ7pmpRDjZAPhSgENIjNTTFWM0VhOsCcvAa1QukuYjWbmz5TEP22asqvRAOJ75aPlMBYYUTDgr4BuqdkmPVkd+kHppWzjCVNR1/D+kzgzAq6uhnYgJ24nE5psS6uattoHdQFbQLx8CO/5s+/WQEIHQu7+fyHf1A2k/2IIw7n0dtu559Hz0GDxAcqpd5dUdRIMcZrEmVy3xDzFvfQ3JOnsSvERBAZ4c6jZzKQDQCTetC9ALsFJNxlCtc6i0n8xDz0+9OTRyVp68crbsV7b4raYSempWVNha6OXKxxSzBOpmOlUNWPAB8RkV+q6slj6aNuOk4Q8vk8P/jB9agKTZMGOPuszw47XeTMF76SxsGIxv4wRUIkm3qhiAyqhig0jloXVA3Tnx6gtSti38MO4qS3v4HLLrmUwVzA5A1DqJqYmADigLJdocVJUIrOj0pqvEkJfV/MWHozNO7fmb2aOq+4vQYy/Bpr5en9CYWIBMCksbava7QJwi23fIvOjbba1VtPezEzpg0/A7mjfQqhEXIDEb1NKWYR4nizFQy1BfwUwnxAJBEmE9DUk6d197l8+g3vivuUlkk09uatEEQ+8OzsQTVIhlhIREBDdTJgbHFVcZWLA3cOVnv6n2o3GlQFYyzhgUocv/NVk9PBgXRcDlUyg+UlqBaK86hqKCK9IjI5PVWmUtQFbQKwdetGbv79w0CWww/byuEHj7y+uIgw0JRl8tYBNrc3pqh2ibWDD6BpaH0s+7waGIpo3TLEgcftV9Bn6+xZTF2ylCeHAqvCAsCXCVdFQxydqLF2U3GCpL7GcUSkzgY0Yq8beYYiMSO9IIKtdOwVlfX7EqlJFLoy/al+wtyOo+3HCf3Af0TkFmyRKgBU9QOjNaybjhOAP//+i/QPZBFRTjzxYxW12fM5R7DHsi00dw06UzF91DMSVjtpZF9rKEztHMBEymnHvaKgv1lTOsjkFQYF8oENvQ6J3SJjN8coxoHtSCCfjl8ZJ4TGmo8hlpGMWcXCOTnp2Jk3eaMouQdvPjb0KPP/28MJp55S/sOoAdPR4bfAecDfsNNj/DYq6hptArB4mQ1Kv+y4iDlzD6mozWdfdxqn3v8AC1ds5eF9psfBXZxGif2dUOJ9iNK6cYD+KS00NhSWQujp7mYoCCDCaqKhIMnOiBTJYLM34viZUFLC2JMfOLIknkEqcZwNPI2frNamkRKJTRMDiTNGAkeotq8bIN+a47XHlcn/3Q7WcUdDVb831rZ1QRsHqCobNmxk48Yt/P2Oe3j8qUFOePmzOfnVlTDJCU4/+VR+cfX3WDJ3Mj1NDUmKogZWTUTOfwuwrwUyecVkSwO/4ueHYbWgoJ6Nt0s15bEztGO+wgmzdbasUBrsP6qxDxcHriMrqD7rRApyK91yUZFijBXE0NnDgVHaNg4RdIyQalYDPhqAiCyhzGhUddFobeuCtgOxtaubT573Nfp71xNI4UKQL3zJiCGdsnjlkcfwiz/9kcMfXsvfDp9P2tKXvDMZDQlvrtDePUDH3vNL+po/Zy6rH3gQMyBEGawEGB/AxpEsBkxkNZ4CgZNErFvntZt46kMkiadh/TbrdyWUv6rrVsQxnpZYMQJGlYUPddO+fpBj31p+IUOhNsgQh3RuZSPwOmBqJQ3rPtoOxAUXXkH31g2s7pzG4+sXsmzzHMJIWLTPUbS2jI0ZPvyQQ8nmI/Z6upPdVm4m159He62fJQoyJEifIENC65YhJvcMcObLSsvBv/ro55ALIxryeSTEFv/LC4Q4DecWu4gMhE6YQs94OvM0jz1fAUxRxWBvN9rXGnp/LUnPau4eoqEvTzAQcdAdmzj69+vpWDXAe855P287eoRpYzXio6nqxtS2UlW/CbywkrZ1jbaD8MMf/oL1ax9nzYbZ9EXNZLNDzGtbzdSZe/PRc08fc79vf9HL+eeddzClb5BMfz97rNzMP/fcjYFAEJyGcH9bBoboz2Y4YPdSS2bx6tUoMGidMauTPOceWa2hAcl0NR94zjvm0WAFMgAitZyJ+gxG58957eaWevNabcHSLhYs6yLnpv1EAvmpzbzznLPYf7c9mNQ8wo+Q1o5GE5F0hW6D1XAjT69wqAvaDsA99zzM3/72e9asm05/fyMS5Jk3fTXTpi/gK58/tyjLYtvQ2NDADz73FcD6fm8+91zmbdzK4o5piAgSOa0E5DFkw4goikqC4UvXrGXIGAgNJrRC5ZOExWAFJCYdDOQdlSFWKCO1law0dHx95GNonj0xLmzgtKDY/jL5iN0Xb2HP447m3a88lcHBAW579H6Of9bhdEyeUtmHUCNkCHBR6nUeO/Hz9ZU0rJuOOwBRZJ/0ru5mjCqzp2wkMMJnPnXOdglZMUSEoH0yTYNDGC9goVM4KnRlmxBVrvhV6eImxx54ILkoomFgyAqXMx/TLKOEEpMmdi0new2NEjIFcOabI0DimiD+mM9asa9nrOtjMJfho68/jcmtrUyf1sHrnvviyoUMaqYKlqoen9peoqrvdGXyR0VdoxVBVbnhipvY0rOF3u5e2tonc8hRB3LAofvT0FBaPRhgSnsbIpDJ55nU3kfLpB7OPffDTJrUvMPHlwtymKgXGfK5GIABGQSNMjzdMpXg97/l0489yf/70AfIZOxXPG/2LPqCDNN6elk9eTIaggkEBm3NEu+GYXAaUpwFKEmWiDpzMiIJUgtgHOXvA9viaXwIjSEQs23ViYtRZdNRRF4FPKiqy9z7zwCvBZZhy4IvGa2PuqClsOyp5Xzrgm/TFW5hqDNEQ8U0Gu586B9EvXDpNReTy5VS53PnzQSgqamf6dM38ZIXn8T+++85LmPs7+4mLzYAZbwmGUpCXstaO1jQtYn1ix/j9A98jEsuOJ+2Vpv6NXneAmauWcuaSe0YtfPSCMBEQhS4rI1BG2cjUO99YckNZ2f6lCu3Cn28EK8XvAgw0DAY8qwn1tPWPUjAdqC6AWqPL+JKFriFB98CvBE4FFsa/2WjdVA3HYEoivjuN6/hC1/8POtXbWLFz7pY+csuVv26ixU/3cLS73eiRrn62+Xjlfff8zAaQY5+NJ/j9W945biMs7evj8Ytm9nYOMlmduSxghFCEIbM6tlMU1eeu3L7cI/ZE5Mf4MLLfxC3f8lzn0v7QC9BnyIRBKFgBkHyNigs6kzH0PZLHj9JwB4PDeSlsMSkOx6nh7n3ezzdSXMkvPpNb+DLn//Cdt13DZiOqqq97vVrgKtU9R5V/S5Q0WJ+z3iN9p97H+LbX7sEzYZsvKufnsXWcRFjmTzT0IBGEZ0PhNzXfDfvPO1ectpMS2sL++y/Ly87+QVcc+UP6OvKkIny6MD4/HatWLuOj33pK2TFsDkzCZMiCASY393J7l12YZPl0sG6aCprdQrBkw/x1JKlLF21gdv+egehGDJDVkY1CwSW7DB5G3BWg6tuLBDaGiFxUNxfDGJBs1NifPBb0EBpHAiZu66bN511Jscf9uztv/nqazQRkRagF7v02KWpY42VdPCMFbS+vj7O/+AFbMqvond9xKZ7hgj7IkhNtDTZDCYIUFV6lyorNwzSODdHtrWbzbluNnav4a4Hb2NwwND5xCQmzRqkoWl8PtKLvn8d5PPc17HAji0VU9IIlmemEuSUBYMbma8bmEwPT+TnMj2zhQu/8mUihT5pYHluJqiQ7YOhEDQDZAFXV0RCILRPdmQEiUCNrdLlp9C4dBKXXEySe+WYy7nruxhoadoxQkZNpGB9E7gf2Ao8oqp3A4jIocDqSjp4Rgna4OAg/779Pv52y+08tfwxonzExn9HDKwL0RAkmwVj4jw+jWxxUv+rnR/I0PO0sWnnmQw8Ark2ZbAvC2IIBwwDQztwpcsUZs2YyaYlTzKQyVnTKXT8ROiUjRqWNU1nXTCZw/sW0yZ9zA46eTS/EJGQ0BhUDHm1pqYAmUHQvO0jzDoNJ/j5LDarX1zIeVBdSMCxjUEqudH5eoKQHYyYs76bI47ZQQu71oCPpqpXi8gfgBnAA6lDa4CKgqTPGEF7/L9P8dULLkKyIQOd0L8Gtj4p6CBJ4c8gSIhqESSwwSYVgWwGl6iXFEkEBrszaGCFs39LhmlNfaxetZbZc2bu0PG/6rlHsfj2v9IwMMSAZH0ur/0bJhZdv8nxn4bdOGhgGXOCjTQwxNJoFv4Mk4coa881LqgcZSCTh7xP4PcxaCPxrGcR57s5hhKTrATakI+YurWHGZt7mNnZw0BTA//zklH5gYogJPdWTajqSmBl0b6KtBns4oK2ZtU6/vy7v/LAPQ/RObCGvg3CpvsE9UWAjYHGjH2ywggJQ0xgiPIuAiwgTY1xrEkCW0iGIPBTt+yvf8a+z0c5hnoDrv3OjXzisxWtfVAx9t9zDwaCLFP6e1nXMNlS+kOkEo2T+HGPaWJxZiaL8muZFmxlWrCVh8KF9NJoH9x+0ByJoLpsEDPgtFrsf4H6RGN/v5HYc/IR0/p72G3jZjq6eslnMwxOn8JbznoTxx562A6NH1Zbo+0I7JKC1tPdw5Xf/h7/feoBokEY6DT0rc/Qu9J9+UEUkx3iUsyDbJYoiojyISbnhC2bgUxgf9nzISoGcllcMWA064Q0a1x8Sdm8upnFuYfo7e2juXkMy8gOg9XrNpAJQ4bETnWxOYpOUDT56/2oddIOgTA/XE+WkN2D1Twsu5Nk3CfPbzAAYQYksALoWcXIQDBkTw69Wamg+ZCD1q1helc37XvtwUlveynPP/CgHStcKdRKCtb2YJcTtLv+cQ9Xffcq8kMhmx/N0LdGQAwmcFP0A5BMgOZDK2xgzUNAjIGcQaMIyWSsoOEmLpoAchkkm0GH8kggRCJoYJIEehG6ulpoD/u45uobeO/7Ttuue4miiJv/cSdPPL2c+27/B31BI5tlEr6Evr9uyfLVYu9pfdDOep3MVO0ir5n4uEbY8ICt02P9tdDmGBtXs1WBIONm4vhgNfa6u2/upKOnl/d/+MMcvOf4xAsLUP2A9d3AP7Elw29T1f5t7WOXErTv/t+P+NeDf6VrZcDmx3LxHC2AKPRZFE643F8xYjVVUFLcwj6RkXP6GyzhYZ/KwM6naszYPMAogoYsUWAFb9Omydx3/53ce98RHHboAdt0D6rK9X/4E8tWrOaxB+8jGOhjwGTYHLSwrGG6y6r3Jxf5L54kSWk2QeikLU7iSLcRl15Fxp5rXGa+AlFg+woia06KOoZSQxZt2sDLXv/6CROyGmAdj8aWC385cL6IbAT+APx+pMKraewygnbrH/7Bvx78K5seytCzwhc61OTXsNisyWZsFjq4KqOuTUBMboCiA0OWFME5/yr2/CAhRCQTWHLAGMgYurpbaGga5LLL/o/INEI0RMfs3XjH297I3nuUzhXzUFU+/fVvsfaJR+kJGujKNPF0+1zCKEi0VvFfD7/SCiS+V/rWo9L3/ndFHZliIoga3OvQ3qriWMqcDWwjNjXroEV7DP9l7GhUn3XMA7e5DRGZDbwCWyJ8T+BOVX3PSH1UJGgi0o6tM/4s7G2foap3pI7vC1wDHAZ8SlW/Vmnb7cXWLV1cfNF3WLr6MbY8laFnZTZ5ygLjNJSLwkaR3YLACQuOmi4KMqdqHgbNjeRRNJdBMwGStyXgjAihWxbJmlYKgdjn2QiremfQJS1IJiQUw+DKFXzzsqu49Gv/b9h7ueL6n7Lmicd4cPJu9IiLg2oZoRlJyNKCqGX2Ffdp8KmJNuti0BEl6fau71yYZ7fe9SgwaQf6n6Oh1nw0xzZeDVwtIgZ4zmhtKtVo3wJuVtVTRSQHFGfLbgI+AJwyhrZjxu233cXV115DOAidj+QY2BgklZmM9c0wbh4IuGB0BjLGPUiKZgQJDIShTeXLGNRlS2io5FFozKJBYKnvBsuN51E0k8TcIlFMxpAPQCfZj3XrQAs6BPlmIchH5HrLr9usqnz/57/i7tv+wpOts+jONtq0KHDqhTjTvuzkRy0SwNhsLH+u13iKM8sUX5MnjtEZIBqy+42BYDDPEVsWk2lt44zT38fcmTs2fDEiakzQ0lDVCOu/jYhRBU1E2oBjgbe7jgexmXDpi60D1rmlbLap7VhxyYXX8MATd9K9IsOWJ63XLib1nWQylqr3PpZYskOjxJzUIDEDo8A+UaF7ukNsEBa3miUmqbbry7tJxiQaMSMMGoFcSjuKZe5UhLwYiArLGwA8tmQZX73kMrR7C4tbZrIu15bUbtREEOyAk9fxWIqFKDaVKRUwTQ4Rx8ew/pkKUaBIZGOLmoFgCMKc1ditDGJUufzCr4wbuzgcxqrRRORq4ERgnao+q8zxYS2xHY1KNNoiYD1wjVvj9x7s1IBK1u6tuK1bSO4sgAULFozY6d//fAcPPHEnGx/M0rs6uQW11T0hG8RMYro2vHqfzAiSy9g1i9yi0ILL64sLY0iSXmQEkwreaqSoE6goI0RZQ5QRqx39yoCehMkABrvEbTTI6e87F7KNNE2ZynGHH8gff/0rtmQbeXLK7gyYbAFhES+XVPygea2U9r3SghQlfcTH/OecOte/9nU9GLSeqKT7V9CM0BjlCSUz4UKWpH2NCdcCFwPfH+b4SJbYDkUlGbAZrMRfpqqHYgtHfrzC/ituq6pXukW4j5g+ffiE6DWr1/G9H17L1sWZAiGL06JcGhU4gTASO/VkDNKQsTR9YE3E2BILjCM03D4j0JghygVETjPGNXCMYAJDGAhhcwbNSCzQVrjtefkGXAYFdGorS3QmK7SVdYOwcc0q/vSrX7A218Z/J81j0GQTdi0lKJ4tLPHNUiZkWjgLyJCUCRkfU7W5jE4bilrrmrySiWw9HoP9KyEE/SCq9GqOjA7x0KNPjPiF72gIjDl7v4Klddep6r9JFhgceSx2DetPisiVbiH6q53WHBWVaLQVwApVvcu9v5HKBW172pYgiiK+fP5FDHUZtjzhmUVJUqJccUEV4qKBkgniVKGY9IitR4GsNQ1DR9vbJU/UWV4SP6FK8gUPNQX2GW9wM8I8++i1mELYCBhDPlA0ByEB67Pt4GNTro5iXMY3LWRpaOFLU+ybpYXNn+Q2X3YxVrDObDYkx6G0raT2G4WgC6Koga0Nk/i/K3/IFV8/v9zXM36oHR/tl8DfgT+RBFkqwqiCpqprRGS5iOzjpm2/CPhvJZ1vT9tyuO67N9E7uJmNDzYCTsCc9or9MyOIIz80cMRINoifqtAl8gVGIJfBrj6pTshSZmZo69JLIEQZY/+qEnq/zgAZWxbbx+h8jmCUwwk/BX2aPG6CJfZB9/tJLeGe9sdSf/1bu2qMFghXiZClG+CFUQvP9ftTmtIzlxEuD9KZkb7A96ahyUzqXptk1EwQpPwqMx1F65ldqapXjvNQmt2ChNuMSlnH9wPXOdZwMXC6X8xNVS8XkVnA3UAbEInIOcD+qrq1XNuxDHRoaIi/33ErWxZnCPuwvpXPMnfpT5ozpQ+AwS3GYM1CslbThW7VSzECxiYTe0Gx7azJF2W8tiT21zylLyLWCgvE5gQa0KzYCtsZBRVbATg9HsVOOXEMn0aO8UtfO0Vo+FIdPo8xLjwqxLUWY1LECY3NwUy0nVfhwwpkysREwETqJoL6cdjGfUMNZJpCVq5ey7w5yZpu44pyPqrFDlvDehvwGxE5QVV/t60NKxI0Vb2fwuKRYKdw++NrGGbR7GHabjOuvvRHhFFI97LGJA6mWIHw63EFllWMf/WTUWCymUTXi8Qmm3p/DicogdVwooJkDCHWX4ncw+t9uAgrUATiioMq6pa1VbcfrDBIljhL3k9F8b6cOEEr0GapYapTd+nYO+nXxcKZ/hul/sY3merfvwhTO51DWzweDSEkQz4K+OOt/+KMN5/ERKGG4mhnA58UkUESv27XWYiwt7eXux+4g61PZhEx8Y+0rW0RxA+uydhE4VBxCcFuWosqYaT2vYdvY2x56qgxiKfLaKQ2IO3Mw7zXQo6qjwPUTrtFaoU7Jl1MoW8U5ZWw0fp7fg3omIaSgiEVCpO9RLy8ke9PioWLwuulhUSilMmY8r9QbE2UtI/n7ksiCgcTa0OhL2xg3dqKZ4fsEIw1BcstrfsCrJm5AvgsdpprJZZYCVS1ohqO5bBTCNq3vvQdwkHoXhlgjLj8Q28VueUeMkESM8sGsVlpAkMYKVEuKIwfucyQ0O/yjCEgGUuOaFow3cRrdVrNP3q+WO+Q8elcYEnyhJ/xT7+vwh05LVb8Q534YXallgKtIokwiILm4934Zc787SUmYrL0n9eKkjrPpN6X+Ieh84281nUn9ucbWLx4GROKMWq0CpbWHdYSGw4ichI2Ngw2wbi0tl8Z1LygbdzQyZLV/2XLE1nABaBNYDM0wshqDMWSGT4lyi3KJWIJjKjR5Sr6ZOHIxsxCT8kbE5t/PsgcZby2ItZQsc/k+rE+HIk/J0oQGEJj+1fnnwV+wqhv7/224swv36cTzvIcAInqiiXTmnYSpbp0WjOKrCKP3OKh5TL+42s54Yv9tahQuAGyZohZM7fp2dw+FI+1ihCRrwDPBq5zu84Wkeep6qhMek0LWhRFnP/xLzPUI/SuMtbvcRkc4jSX9c2cYePSr0xgLJOtStgY2JhZ+ssyzr/KOAkybrUTL2CBQdSuWqlZQdyyQ4jTZl6r+LxivwVCHrUL6qVqrIWOJYycEMdC5uTey4ym3oi3evOJuegXbE+X7TBihclePznf92xEQYXA2AXfPXEUhwlioVWb5qUpZtNZDgi05bppzfXQmuth4cKFO+Db3QbUiKABJwCHuLQrROR7wH1UELKqaUH78qe/RW+4hQ335VxA2j29Yh8GxQqTj5mp80fcp2BTrEyhkMUvU9n3HlEgTlu5szKCCcQWqREIxeZGxv14v0zdtBJwbKQVili9OOGSYk3mNVfy0sW9U2SOYyQjkqF5EhRnhsY35TLyfYf+/IiERYxIhCxe48xNJPX7ffiA0H5MLQ1bmd2yno5Zu3PQgcfzhteOsCDFDoZQOxrNoZ0kCD650kY1LWhPr1hC19MZ8n0BNORiwdAiVlFRyzjG9TywicEN7vbSkiYQZmwWSPxMKjZ9Kif4blUhDLDVfLFJxD5dy3iTFGJNZzJ2TGqs0PmJAV6AyDpyL/2JF8p5jChe84jYrBOnzeLfhjDxx2J2MiTxXVHCPGRSJqgXUI00VmI+xiZ+1jaWOZWU+TilsYv2KbN4/lHPpr9viA3rNzFjZseExdIkqhlJ+zJwn4jciv14jgU+UUnDmha0qVM76N641qZWeSETsTU6VG0puDC0WsuTIAKISeZSCfZ4FNkfb7FC5afJRJE6zSduwQerKfNNzl8jLdDELONQA0jofLSMuOXFJNZsEcRjUOOTi0lJXvJXUv37yxXMak49z15ojPfLBAS/ljQxkxpFmmhWUsSM7y7dd2QPFvhvMZEiBBKyeeNafvbzGxGUP/7pZ4RDGb70lfOZMaNjDN/sNiC2qasPVf2xiNyG9dME+JgjVEZFTQuaaBaTsxrM09ySs7S9qvsB9iW63eLlIkLeZ4PgQ0mK5gLrl8X0HUQiSNoUdIIVOpNT0vuzYIxBUcsaBoJmEvcm4dVJ+oM4kO3JiWKhEUl8rsBpyrTAgRMcZ+L5y8RZKK4jcT4aYaIRJUoELBEgK8SB07ZxBS3xH6Og+cjOvgYkDFm+dI7Tbgr5kAYdZNbem/j5jX/gf9/z5m36TseCas+wFpF9VfXR1LJNK9zfOSIyR1XvHa2Pmha0UIYIQwO5DFE+SrJB0qaQuMmWgRUM/0vuH94oECuYkASRfeFP3EPotV8mYSIhIR6iLLGZpAoapB5+d46IK2CTZhLda+Om24Qp7r2EAMGRJiltF/tYrp0JiBdmF9dv5Hwp3yauzoUgklp5ExJ632u3uO6IOt9PbL0UL9DikqxFMdjsF1EIhzL0bGrkrjvumhBBqwGN9kHszJKLyhxTKliMsGYFLZ/Ps3HLWob6G5LSAZ70IPlVjwRLkjifiVhovAA6AyhwgWaXDm/dOTvvLApwRXmw5qR7IG0mh8RkRIRdEjYS8PlPlpl0WquIrlfFpYB50iZlwglJfqZ/kExyX0Z86QQSDZy8jCXG6tjU9UgIkUgtwVGkSL39G+8zkUDkjGQX5rDjcN5eADoYwUAeCSPIh3SvzdE+u5unl65gwcLxpfurTYao6lnu5SuKC/OI+KnwI6OSaTJVwfe+cyNIRM96O2Xep0qJQJQLCBsCIudriXjfyvpvNGYKYkURbpVYI8lkzqyQNxA2miRH0gtFVtAGieeheXtNM5BvFPvz5PwxbSBmFVUSNjQ09pp+crdCvMa7P9c/P1IiCY75S5ujKX8L3I9MpDG17xdz94u02wAehT5j3HnKL1RFw8h9Vs7GjJItQGzRHhQTWc1GEJDvCxjoCrj4omu3/cvdFsRsVdFWHdxe4b4S1KxGu/1f/2DruknkTQaySVpQ5BKDRUAbMrHZ45NztcHEJpnPMvfmmbqYkjqhi9x8EUlpuTDrLUf7tIahWuLBkSji2vqANP5arkmInRbjyYzQmXLxOn6er/EBcpJrpX3CAvM4JRixJhxKWvoxeC0YuQGJC+BHYZKDKak4W6xdNdWJptazcGEGN+EoYTjd6y1LJ5FrWc4/bv03zzt+x9TZL4ca8NFmAXOBJldv33+EbVRYmqMmBa2rq4cgM0TXpmkxqaFhCBibTY//UVNUBMlaar/4d06zxv6SNzitpyn2z7/HkW6iVjvhmEZ3TmQsza+2wD0gsfmokFDr/vyUpvKzsrVIY/l2/k1sopX5LAqoeSFeE00DywbFlqXz5Yy9AavhPOOY6jh9fuynpRgTM2TzPEFdCEXdLANNSjeEVgv2bwzoWZPj+h/9dNwETYrGXyW8DFuOYx7w9dT+LuCTlXRQk4K2efMWurc2E5IBXDa+ZCyNnmLeohDrW4llBP0UEvGaT6yZl47LxjUbA0ccBGLLEDgTz2d6KFhTLJBYi8QEC05YXdHgglolKcTC73y1GJIwg8alSVUCJRVjIxEaI2J/HFxdRv/5qMGu1pm6rs84CQNbNDU+ZOzC8JZtVH96gWMoYokR8nm7ywg9axuZNHMLW7ZsZfLkUZPYtx3VNRXdEPR7wPdE5LWqetNY+qhJQYuiPL1bm0BcwR1Hp4sxtmy12l9VgTgzw2uOyFN1LnYFkiTaSuIrGffrnM7gMOKK8vg+XVeaEhSfRCye8Uv5YJblTJER/obEBr4jTVhAr9BChSDFYKZ4lpQT52l8YvYwJktIBINUkzh3Mz1hIU8S5A6EKK+YlOmNep9S8JF7Y8QO0gj5fETgPsTAGCKNGOjKEA4afvnTm3nbmRWtm77NqAGNBoCq3uQKUB1Aal00Vf3caG1rlAxR+nty8fdthcPS474OSNiUIWwwcUJvPmdshr4jSKJAIGusyeeL5zQGhA2C5mxOYuSyOPz3GNPvAUQ565NpYOn9yP1Vt4XOR7PJyfZ1FFgBBKxwZYEsSMbek4kLgBTIUMpXK/4UUucJMSvpg+ZeWjXCrknt+zbWL7MZ+xIXbBYtZCBVcBNfiYsw+/INKtZcjPJR/KOGr7GSMYSh+xlRpXdthjtv99UqxgFaZqsCRORy4A3YycwCvA7YrZK2NSpoEDSGVjhyhqjBuFxHK0DqZj2bbIbICGGDlN6J11IuJUuMiRlHsPGhKCOEjRILmwqEDRBm7WTSyAmdBlaYorSJ6J5YO8nTvjaefXR+YJhKAg41mR0NxA8zFDKCPogdnybJe2/1poU1Ln1X9AAW8BvOX4vjhfGFcNkwhkjs5xAZ4pxRE2nBhNN4wffAoLmMLYUeRvSsyZHPdPPgfY+M9JWOGTWwtK7HMar6NqBTVc/HFk4dvvR0CjUpaKqGhuZ8LBRIEnDWjMSJu1FcPs5x2p4d9Hai7c0mNKgSujgY2ApVYaNNPNZAiHICOZtAHGE1FAG2IKqru4qj0KOM02YBsZBBbG3ZqxYJReznFMfaio7b+09eJ/mJiS+XJis1ZXIiyUeW7i/IE2tqNVrQRiHJGxWJFxuMN/+5401ijcMDpiGLijDUE9C7NsOVF1/DDofif6UKt+rArzLZKyJzsLOsd6+kYU36aNY0Ejd9hWQ+CBBkTKwpIC1YJNkbBvtwqH04o0bHBmYlflhNmPhT3gqLlIKgcgTOJ/QDS0iSgrGmx+60WiB+xrVVCZ41VMoTID79qqBGvtcmTlv55zzNWHrWU6PkPrxf6V+nZC7WBj7BLNZu6pJuvLaMNE4+dg0IB0Jbmi5Sq0lDd6EwwgSKKf4V2UGoFR8NWzOkHfgqcC/2m/lOJQ1rUtAgYih0Q8tIMtFSrXZSwWZ+uPqExtH78S++85289hHBhgBSJEmUmqqiWBbOQzURqLSQRW5feppKQv1R4HuF6rJIkubxsXgGjBcuZ1amBcLndsaNKbyuYO8xKprsqUpBrAxcLC+0X7YYy1BGefUzcGx/aoXLEyZgiaVAbdTMkk1OOJ2v7K2I1t0GaZwW8sFPjLjOw9hRZdbRQ1U/717eJCK/ARpVdUslbWtU0GAoytmsBpM8MmKE0NXFt5rOOu7e9FFDUswUEmJArWbx2ifKYc0kJa6xGMfDnPBErv5PHvve+GM+LIBjKX2tRJNMnRGsZgvxb1KC6f7EhEYarrHPorc3mhyWlLiWPHr+ltWVr/TOocZWo701PwZP7fu6kqqQtylduHvzgovXzO4zRuznHGSV+UdugUg48WWvZeGiBWW+ye1HrWg0EXkA+AnwE1V9ChiotG1N+miA1WimcHjFWfK2VIF9HQlILiDNJIRO0/j6HKE4DecyO8g6X8wTGq7/KAeStVrPZMWeF9jzY6bPbV6YVSyLibHn+kTjAhT5Tr4jofBhKpl+FRWcHndFLAQWxv+wuM8qFqqMLTRkBU0JIrGJztZ5jcsXpO/LzwpIPneDyQU2QTtnF2nMNEeIwMtf/ipOOvWljAt0mK06OAn723uDiPxbRD4sIhX9utSsoDU3Wr8z/dlGAQSBxOlTYQC4OopqhKEgsg9UYNnDqNH6ZZoVNwdNCLMSnxM/tCYhEyJ8RofTfqmHzQuWF6BQNfZpYhMxRYKUex6k6G8hh5+cl9Y+/lDBWNT6dbG75sxI8q5tOptFE5NZXazQm8YFVbDSswvExQwlkVgRSxZ5X7O/t4HerkZuu/WeMne6YyCAhFqyVQOqukxVL1TVw4E3AQcBSyppW5OmY0NDMzOmbWDJ6nlomo0QIQydr5AxLlXIzX7OuifE+2cxle+E1BRpLywhELogsroHSwPnX5FkgngSwpuNaJIR4hEHt0e5t1jLeD/L7fev41Qx5zfhTMGYPDHEc8jSpmc8Hw3ipZfiEIAbl3EaXpx9qoLNBvEXNuATJeP8SpFY24eR2vor3qRU6NoyicbmNeTzeTKZ8XmchqlUXBWIyELg9dh4Wgh8tJJ2NanRZs7swJiAOR3rEJsGYVlEz4rF2Ria/Mq7JzVm3fLYzA+loHBpmoWL0DjmFQXWXPQ0ft4/uAGQEbesrGvnTVJ33fgxSKsrKdzl4ScDxBNZU6YoEGvSwJ9ob6PQig7sGKJymtT43E0ncEbsj5DgKzHESQCRD1WkMvZ9DK2EQnf3ZMMHEn8n3T1NIBE3/vSW4q9xx6CGTEcRuQv4GfapeJ2qHqmq5eaolaAmBU1EOO/TH6W5cYBZU9aTCexcKXEl5IiSIHAERA3uwTKJiadCbF5KNpEGISECjNhlYqMMiePvT8oQVxCOH74iQUqbicWvTUqbxkj7Tand8VSc1M6w6Ffcl6CLx+D8Qd9Xwdi935hc1mphr3GDlI8XYokhVwYciGdm++8irg1i0jfq2kmGvt5GlixdzPhAE1u5+tNkTlPVw1T1y6q6TTdck4IGsHDhPE4/4ywmt0YsnPE0LU09OJ/cspEBNrXJFSxNkxRhFmg0BTOq/TPqraQIGELJu4BzhDMv0w+0ZyxTD643mYrkIvbf/BuNWY7C+/JtotT70CcfewFNt/HX9UMqeNCJBS4eo/9N8Z+Tkw3JSRzSUE00l2SdNZCxmSHT2jYzZ/oaWpp77OeihQHuKPZLkx+S/v4GliypyFUZE2ooM6RTRK4Skd8DiMj+IvKOShrWrKABPPc5h3DJ/13I5GmLmDdlDdNbNqCpyVmphAb71/jFJiQu0+3PNXEQyuUxZoEgSZsqyPxIB6S95kg91P6aIs6kK7YPKfzRTU8ejf0mkvf+uFUwiW2r7p7SmivO+k8LpRemeNDur0upSscACZzGMt4XdDMXRGhr6qGjrZNJzcqk5t64o1gzRtaq8JkiHmFk/OT38UHtaLRrgT8Ac9z7x4FzKmlY04IGkMlk+PoFH+a0t59Fe3M3e89YyuyW9WQDO/MxnRVhyxLYh8MYIa+uTIEmpljYmMrAN3YmtBckXz8/rYkKzDIP7+doKRWf+rGP4YPS6eNadFwg0U7utc1QcYIXpLSqlF43/Tb9i+8zUnBa0wtLvIZALIwR0yZ10j59T+bMWUAsSinNlZij1keL5/UZGBwap9mZSs2wjkCHqt5AHP3QPBWuk1bzgubx3GMO57JLvsGJr3otbU0Re01exszmjWigjt7HOehYQRAgZ399jWPK/K97LESOZPD0vq9WbHwJA4pMwgpgtU6ZAykfzvepqX4LzFA/fp/s7JsWC3DKd/TkiBfSuFO1Pqm6C8dhC68BRZnasIV9pi6hIZPn3We9gUzQSCabPD9iii4c34tVpYLS2Jgb+YPZHoyRDHErcq4TkYeGOS4i8m0ReVJEHkxVuRoOPSIyzY9ARI4GtlQylp1G0AByuRynnPRirr7ia5zy2jfS3tDFvi1LmZzrTqh3J0yhaGyihU6zRQEl5lWqKG/iP6km01D8rzZlyI/081f8WlLPREqDauq4f684U9A/0E6zFWSW+MC7j4+lQwvp63ttZ6sJWVkIwUTEpQ2AmPKfktnKnKb17H/IMVx68TdZtHA+J554HE0N/RhbLMT5rfbHR1IaLv4sAq148upYIKolW4W4FhiprPIrgL3cdhZw2Sj9fRD4FbCHiPwTuzb2+ysZyE4laGm88mXHcfnFX+egQ45gt9xqFuVWkjWDyQl+iorXHL7qcEhcwz5e5M+nGzkywU6qLCU7YqKE5OFOZXuVkB/x5AN3rPjxiON9xs1pk1TZA/ejIYEVSF9lz8fu4mRrf73UgI06wfXXTWs4v0+hJehlVm4je+x3GB9671vJuRqZBx+8H/kwx7T2zoLx2giAJuMwIEZpa+1mzz33Lv6KdhzG6KONtoY1cDLwfbW4E2gXkdkj9HcvcBxwDPC/wAGq+mAlY6nJgHWlaGjIce57T2fFypfz9Uu+Q8uGpXTRRLc2MSRZVJU8AZ2mBXWF2WJ2Lkf8M5NPERACRMPQWt63EkpNPUiELx2rgziFMT7HrzxF6nwo//z4YlhKUb8BLhGz8JpemXlZTSpnaXxOToeY17CONtPNnN0P4Jx3va3gmsYYzjzzdK695gr6Bxro6mmh4ANy2jCbieiYugEjyrvfNeIKSWOHv6FS7IildecCy1PvV7h9BQvAichrhmm/tyuy9LPRLrRTC5rHvLmz+fqXPsMv//h3/vvIIyxe/BRRvs9KxlA/e4Sr2Rw202Ma6c40sjHbgmJr6KfjVels+hKkNFJaRtI0fYHQ+NdpP8y992RMmh3EETpadMimSSXjioemqe41eZ8euydofWkFoxGzgk3Mymwi1zKNM08/i0MP3KfsZ3rMcw7jrn89F334dtrbtpLP2/W++wdyDA7laG7spa2lm1xDM/971vtpbm4q28/2QhjWVNwRS+uW877LXexVI/Sh2CD2iNglBM3j5Jc+n5Nf+vyCfVEUceudd/OzW25Fe7sY2LKa3QYzPNE2i82mOQ5yQ2k2fSw3UqrF7PmUGN/GTSXxaV2x4LiM+ngdRSnSYGlyxJuBrn0Y+TSp1LG01KemaJtUe5/OBcoU2cr8zHoymQynvPqNnPDiY5NA9DA49+y3sezpF3PllT+huTnLqtUbCfNbCcNOmps72GPRwbznPf8zbqlXMcbPAVxB4QzpecCq4pNU9fTtvdAuJWjlYIzhRcccyYuOORKA7p5evnrlVTQ++hCPtc9mXdNke2KhZVRWe5T8/nlNkS4PXpwJ7OBNwDQzXbDWtiSnp33BWHP59wWVUO15gRbKnTcRDdCo/eym62g2fex/yNG87x1vin2xSrDbgjl88QvnVnz+DsfwpuOOwK+A94nI9cBRwBZVHXbdYBGZCXwJmKOqrxCR/YHnqOpVo12oIjJERNpF5EYReVREHhGR5xQd31dE7hCRARH5cJn2gYjc5ybLVRUtk5o5/9z387xXvJJ9Nq+hfaCn0LyL/xkeaXbTZ0+UyGBKcFzYqVBgoKDAjm0jxdYmEiZTxsiTCJofo6vFX2zXZqOQ+eEaDtClzJ09g8+d/zk++O63b5OQ1QrGyjq6NazvAPYRkRUi8g4ReZeIvMud8jtgMfAkdqb0aDNXr2WMAetKNdq3gJtV9VQRyVFanXUT8AHglGHanw08gq3sWhM47ZSTuO/JJRzw5KPcP20BPVlXPWyE7zD2yVLSELtiqQz71DocyTkpwqRg9nSqk9gP8/KX+gvEi1TEZmcR6UIEEikzw83Mj9aTNw28/R3v5pgjDqn8g6lFjDETpII1rBV47zZ02aGqN4jIJ1z7vIjsmIC1iLRhF1y7ynU+qKqbiwa8TlX/jS1WUtx+HvBK4LuVDGgi8Y0PfYCZe+zFgZtWkAuHRo1Lp0NrMenhjkUpqi+tYHwIQFJsSrnnxhMWceEbR3HGMh2RLFifGodPxRKgLerlwKGlLNCNPP9FL+OaS76+8wuZ/+WpjRSscQ1YLwLWA9c48++7IjJpGwb3TeycnREtbRE5S0TuFpG7169fvw3djx0iwuc/eDZhYzP7bV5tE22HEwInMCaWNAoode/DxWlcpIRQSVKpKLqE90HSDKE/4I7Fiwmm27gIe6SQDQfZe3AV+w8+Tcu0mVx04Vd4++tPScrD7cxQaqkK1rgGrDPAYcBlqnoo0EMFi2MDiMiJwDpVHXUKrqpeqapHqOoR06dPr6T7HYIgCPjqpz7OpPwge3StSzj2NE1OKkZV/B2LSyyOA8kJuVGSveFeFPhgafMwvUX2WOCdvTyuSKrEuYwByrzBDRw8sJgZLTnO/tBHufyLn2JyW+t2fSa1hu3IDNmh2J6AdSWCtgJYoaq+FO2NWMGrBM8FThKRpcD1wAtF5IcVtp0wzOqYxlvOeAez+rvYb8sqjCc3ioStYJNE+0RaNH0l3QhS0ePUB+46UqfNiil7vxiguuTEuAuXyTw138UhvYuZo13s99wXc/EFX+CAvffYQZ9IjaF2TEdUNa+qD6vqQ6pa4ioNh1HJEFVdIyLLRWQfVX0MeBHw3woH9QncYtoi8gLgw6r6lkoHN5F44bMPY/70Dr50wQXM693I8kmFazPHwexUvCv9dSfTVxKqv4RC9G2iAtmLISk71BMeGhX+GjaGAyzqX0tr2Me+hx3F+057E42NDdt387WMsmbEzodKWcf3A9c5xnExcLqnSFX1crd+1N1YVjESkXOA/VV16ziMedyw18IFnPSa1/LbG29gfUMbfUEuEZQCe8/9kdQaZ06oPNUf1wRRCgTOkF6cojBwHS80QWpGtQ90a8hu/RuYNdhJy6z5vP8dp7Nowdzx+SBqCtXVYDsKFQmaqt4PFKe7XJ46vgYbVR+pj9uA27ZpdFXAKS8+nt/d+lcO3LycB6YtYECysQBBoQYqLq0dZ3OkhEuwvlZc3BUtoOdjgXTt00FrO5tYmdG/hYV964kyWd5yxjs57qjtzTzayVAjgiY2lebNwCJV/ZwrNTdLVf81WttdgJbasRARLvl/n6J1ylQO2vg0xqX6p+n6EjihCKA0cwN8tbdkmljquO/OTXKO+0OhdaiXg7csY1H/ep7zghfy3W9e9AwUMuxqIcVbdXApdmELH5/rAi6ppGFd0Mogl8vxtU9/kowIc3s2xTmL8WqhKaTf+jSrePZKlJAaMTSl8bxvT+G+hnCIfbpWcfDWp2mbPoOvfvlLvPN/Th3/nMKahHNUi7fq4ChVfS/QD6CqnUBFqTbPxG+uIjQ1NfLSV57An371K5a3TkUDU1oUZhiLpoDgcO8NqRonFJ3gM0qiiLl9nSzo2UC+uZX3vudDHLzfOM7z2llQI6YjMCQivqY1IjKdCjMx64I2Ak45/gXc9stf0BDlGcg0xGtEF/hpFM0vS1P6XntJak5aen/ciTJ1sIc9utcSAC895TW87uUvGTW7/hmB2mIdvw38HJghIl8ETgU+XUnDuqCNgIaGBhTIRCF9PjufUjctranU+WphVLgPUoRJCk35AfbqWs/kgR4WHnIY55z2FiY1F6eSPsNRIxpNVa8TkXuwIS4BTlHVilZfrAvaCDDG0J9roH2wl65GN7GxiMjQlOSpy+bwtL2PqaVjazbNSpAwz27dG5nbu4mmmXM596wPsvu8ZwJdPwbUiKC53MaHVfUS975VRI5KJXMMi7qgjYIjjjoaveMOlrdMxaUGW4RWoAI3S1tSJqNnGf1cs6jI2Jzeu4VFXesgm+PYU07lra94yUTdzs4HVQgrSpCfCFxGYVZUT5l9ZVEXtFHQMqmZbJgvzeJwf4vnonk/rCQ+JtAy2MdeXetoyg8w9+DD+cw73/4MZRK3ETWi0QBRTQajqpGIVPQF1r/lEdDX38/fbrmF5ZOnUSBOPmsjTYCk/TAKE0qyYZ7du9czs38LZvZ8Pv+B9zJ96pQJuINdBLUjaItF5AMkZeneg82UGhV1QRsBn7vqGiJgZduUmAWJMzYgKcBTRNvHwqbKnN5OduvZQKa1jXe+82yOeNb+E3wXOzu0lljHd2GZx09jv+Y/Y+tBjoq6oI2AVY8+wvLWaUTFS7OkCA8tEjKfRjVloJs9uteRJeJFJ57EG1/58jpdPxYodr2FGoCqrgP+Zyxt64I2DLq6umgcHKBzanOcbW+Te1MnpaoBezTmB1nUvY6pgz3MPeAgzjn9rbS37lrzwyYc1Uu5KoCINALvAA4AGv1+VT1jtLZ1QRsGvoiNhNhScUWBZiNCFCbz1kwUsaB3I/N6NyFTOjjrPWdxxH7laybWsQ3QVEZ29fED4FHgZcDnsAnG9Tja9qChoYEhE9CYH6I3aEjmiOEyQSK1RXYiZfrAVhZ1r0dNwFEnnMSZJ59Q7eHvWqgdMmRPVX2diJysqt8TkR9hq2KNirqgjYAo20CjK9qTnmOG2vhZ02Afe3avpSU/wEHPeR7vftPryWaz1R30LgitHY3mZ1RvFpFnAWuAhZU0rAvaCNAgIPArYpAQHrkoz27dG5jVt5lg1jw++/73MKtjWtXGuWujpiZ+XikiU7Cs46+AFuC8ShrWBW0EhPmQKDAJy4gyp6+T3bo3IE2TeMd73sdRBx9Y7WHu2qiBpGIROVtVvwU84qbG/A1bHa5i1AVtBGREUbVpV+2DPezRtY4GQp7/0pfztlNetWuUc6txKKDVT8E6HVtE+P+ovDBVAeqCNhLCiGyQZ7/NK5k22MXMfQ7gQ2e+nSltNVNwedeHr1BUXTziKrnNEJF0eTlXmEwPGq2DuqCNgExg2L13AwMt7Zxx1tk858B6Vkc1oFU2HVX1ja4A1R+Ak8bSR13QRsBznn8cM6ZN5cTjR1/iqI5xRPU1Gthq3f9R1WVjaSxaO4xODBFZD4zphiYIHcCGag9iG7GzjHk3VY1LVYvIzdixF2ODqo60PvUOhxvLSao6OOrJxW1rUdBqHSJy9w5YbXJCsTOOudYgIldgyZBfYeeiAaCqXx+tbd10rKOOyrHKbQbYpgTWuqDVUUeFUNXzx9q2Lmhjw5XVHsAYsDOOuaYgIrdSpsigqr5w1LZ1H62OOiqDiByeetsIvBbIq+pHR21bF7Q66hg7ROSvqnrcaOfVTcc66qgQIjI19dYAhwOzKmn7jErWE5FGEfmXiDwgIg+LyPlFxz8sIioiHal9nxCRJ0XkMRF5WWr/4SLyH3fs226lEUSkQUR+4vbfJSILU21OE5En3Hba9o5ZRN7vxvWwiFxYK2PehXEPdnmye4A7gA9hZ1yPDlV9xmzY3LQW9zoL3AUc7d7Px6bYLAM63L79gQeABmB34CkgcMf+hV1ZRIDfA69w+98DXO5e/w/wE/d6KrZi0lRgins9ZaxjBo4H/gQ0uGMzamXM9a10e0ZpNLXodm+zbvNO6jewi9qnndaTgetVdUBVlwBPAkeKyGygTVXvUPtEfh84JdXme+71jcCLnOZ4GXCLqm5SO9XiFmDUzIYRxvxu4CuqOuDOW1crY97VICLPdrmO/v3bROSXziqYOlJbj2eUoAGISCAi9wPrsA/RXSJyErBSVR8oOn0usDz1foXbN9e9Lt5f0EZV88AWYNoIfY1pzMDewPOdqfdXEXl2LY15F8MVwCCAiBwLfAX7Q7WFCsMmzzgyRFVD4BARaQd+LiIHAZ8CXlrm9GGWHRx2/1jbjIgyY34W9rubgjUjnw3cICKLamXMuxgCVd3kXr8BuFJVbwJucj+Ao+IZp9E8VHUzdqnfk7G+zANuztE84F5nKqzA+m4e87ApOCsoXErY7yfdxpWLngxsGqGvsYz55a6/nznT8l/Y4ncdtTbmXQRBqvT3i4C/pI5Vpqyq7SRO5AZMB9rd6ybg78CJRecsJSFDDqCQWFhMQiz8G6tNPLFwgtv/XgqJhRvc66nAEqwWmuJeTx3rmLFVcz/n9u+NNfGkFsa8q21Yi+efwC+B+0jiz3sC/6yoj2rfxAR/YAe5D+pB4CHgM2XOiQUt9SE/BTyGY+nc/iNcH08BF6c+/Ebgp1gS4l/YhcV9mzPc/ieB07dnzNglXX/o9t0LvLBWxrwrbu4H6tXApNS+vYHDKmlfzwypo44JwDPWR6ujjolEXdDqqGMCUBe0OuqYANQFrY46JgB1QaujjglAXdDqqGMCUBe0OuqYAPx/3TYh+HqMtLQAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "da_diff = da1/da2\n", "da_diff.plot(title=\"da1/da2\");" ] } ], "metadata": { "interpreter": { "hash": "fa576ebcd40e010bdc0ae86b06ce09151f3424f9e9aed6893ff04f39a9299d89" }, "kernelspec": { "display_name": "Python 3.8.10 ('base')", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.10" }, "orig_nbformat": 4 }, "nbformat": 4, "nbformat_minor": 2 }