{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Sensitivity of the *ecal* parameter\n", "This example shows a simple sensitivity analysis, i.e. running SWIM sucessively with the basin-wide *ecal* parameter changed each time. The discharge changes between those runs will be plotted." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "nbsphinx": "hidden" }, "outputs": [], "source": [ "# hidden cell to setup\n", "import swimpy, os\n", "%matplotlib inline\n", "\n", "project_path = os.path.join(os.path.dirname(swimpy.__file__), '../tests/project')\n", "os.chdir(project_path)\n", "if not os.path.exists('swimpy'):\n", " p = swimpy.project.setup()" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "ecal=0.9\n", "ecal=1.0\n", "ecal=1.1\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEiCAYAAAD05tVnAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsvXmcZGV18P8999ba+zI9wzADzOAoIorLC4pxeRU1GmNQf0mMvu5bYn4matQ3Ro1bNjUhiUbUBEOMRsUkggGNoBEFRQUEZJFBBWGG2aeX6e7a6y7n/eM+tXVXd9+a6eqe7n6+n09/uuqup+reuuc5y3OOqCoWi8Vi2bg4qy2AxWKxWFYXqwgsFotlg2MVgcVisWxwrCKwWCyWDY5VBBaLxbLBsYrAYrFYNjhWEVgWRUTyInLmIuvvEZFnxDjOy0XkW8sq3MLneoqI3Gdkf9EKnfODIvKFlTjXAuffIyLPXobjvEZEblwOmbpF3HvOEh+rCNYQIvJUEfmhiMyIyJSI/EBEzu/mOVW1T1UfMOf/VxH5iznrz1HV62Mc54uq+qu19yKiIrJr2QWO+DPgEiP7f3XpHMuGiFwvIm9YZP0O833lzd8REfmUiCRXUs5u0em9EPees8THKoI1gogMAF8HPgGMANuADwGV1ZTrJOUM4J7VFqILDKlqH/AY4MnAm1dZHss6wSqCtcMjAFT1clUNVLWkqt9S1btqG4jI60TkXhE5JiLfFJEzmtapiLzJuEymReSTIiJm3S4RucFYGhMi8u9z9tslIr8LvBz4YzMq/ZpZv0dEni0ip4pISURGmvZ9vDlestnlICLfM5vcaY71OyLyUxH5jaZ9k2bfx7f7MkTkjSJyv7GMrhaRU83yXwJnAl8zx0632fdUEblCRMZF5EEReUvTOldE3iMivxSRnIjcJiKnmXUfF5F9IjJrlj8tzoUTkWER+bo53zHzertZ95fA04BLjLyXLHU8VT0K/A/wqAXO90QR+ZG5zodE5BIRSTWtX/BeaHOsvxGRG0Vk0Lxf1nus3b1glr9ARO4wx/mhiJzbdJ66G0wil9x/iMjnzfW6R0TOW+o7tMxBVe3fGvgDBoBJ4HPArwHDc9a/ELgfOBtIAH8K/LBpvRJZFEPA6cA48Dyz7nLgvUQDgwzw1Dn77TKv/xX4iznn3QM827z+DvDGpnV/A/yjef0a4MZ2xzXv/xj49zmf5+4FvosLgQngCUCayEr6XjuZ2uzrALcB7wdSRErjAeC5Zv3/Be4GzgIEeCwwata9Ahg13+87gMNAxqz7IPCFBc45Cvwm0AP0A/8J/FfT+uuBNyxy7XeY7yth3p8K3Am8boHr8L+AC4ycO4B7gbfFvBdeA9xovqfPAN8EelbqHjPvHw8cBZ4EuMCrzedLt/msHwTKwPPNth8Gblrt3+ta+7MWwRpBVWeBpxL9aD4DjJuR8BazyZuAD6vqvarqA38FPK55xAZ8RFWnVfUh4LvA48xyj8idcqqqllX1eIOFXwJeBmBGgi81y+LwBeD5ErnAAF4J/NsC274c+BdVvV1VK8C7gSeLyI4Y5zkfGFPVP1PVqkbxj88YWQHeAPypqv5cI+5U1UkAVf2Cqk6qqq+qf0ukhM5a6oRmnytUtaiqOeAvgf8dQ9a5TIjINHAAKABfWeB8t6nqTUbOPcA/tTnfQvcCQJLowT0C/IaqFs3ylbrHfhf4J1W9WSPr93NELtALFtj+RlX9hqoGRPfMYxc5tqUNVhGsIcwP8DWquh14NNHI8GNm9RnAx40pPQ1MEY1otzUd4nDT6yLQZ17/sdn2FmNav+44RbyC6IG8FXg6EALfj7Ojqh4EfgD8pogMEVk9X1xg81OBvU375omspW0LbN/MGcCpte/JfFfvAWoK9TTgl+12FJF3GrfIjNlvENi01AlFpEdE/klE9orILPA9YEhE3BjyNrNJVYeILIsfEI3W253vEcb9dNic76/ayLnQvQCwi2j0/yFVrTYtX6l77AzgHXOu0WlE170dc8+ZEZHEIse3zMEqgjWKqv6MyFXzaLNoH/B7qjrU9JdV1R/GONZhVX2jqp4K/B7wKWmfxbFoqVpVPQZ8C/gd4P8AX1bVTsrbfo7I/fLbwI9U9cAC2x0kelgAICK9RO6XhbZvZh/w4JzvqV9Vn9+0/mFzdzLxgD8GXkLklhsCZogebkvxDiLL4UmqOkCkJGnat6MSwKpaIrr2F4hIO0X0aeBnwMPN+d4TU84a9wKvBa4RkWaLZyXusdp5/nLOeXpU9fIOPoOlA6wiWCOIyCNF5B1NQcbTiNwwN5lN/hF4t4icY9YPishvxzz2b9eOCxwjejCFbTY9QuRTX4wvAa8CfovF3ULtjvVfRH7/twKfX2Tfy4HXisjjJAoG/xVws3GDLMUtQE5E3iUiWYmCw4+WRhruPwN/LiIPl4hzRWSUyLfvE/m9EyLyfqK4TRz6gRIwLVEw/QNz1sf5XuuYz/xKopHw5ALnmwXyIvJI4PfjHruGeei+B/i2iNQUY7fusbmf/zPAm0TkSeYa9IrIr4tIf6efwxIPqwjWDjmi4NnNIlIgUgA/JRptoqpfBT4KfNm4A35K5F6Jw/nmuHngauCtxnc+l8uARxlzfaH8/KuBhwOHVfXORc75QeBz5lgvMZ+hRORe2glcudCOqvpt4H1m20NEI/iXLrT9nH0D4AVEvusHiYLO/0zk5gH4O+A/iCybWaLPnCVyw1wL/ILILVUmGrnG4WPmGBNE1+3aOes/DvyWycT5h0WOM22u0RGi9NGLFrC43klkkeWIHqr/3mabJTG++T8DviMiO7p4j32QpntBVW8F3ghcQqQ07icKYlu6hHRmuVss3cWMtB+hqq9YbVkslo2CDahYThqM2+T1RG4Pi8WyQljXkOWkQETeSORquUZVv7fU9haLZfmwriGLxWLZ4FiLwGKxWDY4VhFYLBbLBmdNBIs3bdqkO3bsWG0xLBaLZU1x2223Tajq2FLbrQlFsGPHDm699dbVFsNisVjWFCKyd+mtrGvIYrFYNjxWEVgsFssGxyoCi8Vi2eCsiRiBxWKxLDee57F//37K5fJqi3LCZDIZtm/fTjJ5fG2srSKwWCwbkv3799Pf38+OHTuQ9p061wSqyuTkJPv372fnzp3HdQzrGrJYLBuScrnM6OjomlYCACLC6OjoCVk2VhGsI1SVMLQlQyyWuKx1JVDjRD+HVQTriBd/6oec+Z5vrLYYFoslJtdeey1nnXUWu3bt4iMf+ci89Xv37uVZz3oW5557Ls94xjPYv39/V+SwimAdcce+6dUWwWKxxCQIAt785jdzzTXXsHv3bi6//HJ2797dss073/lOXvWqV3HXXXfx/ve/n3e/+91dkcUqAovFYlkFbrnlFnbt2sWZZ55JKpXipS99KVdddVXLNrt37+bCCy8E4JnPfOa89cuFzRqyWCwbng997R52H5xd1mM+6tQBPvAb5yy4/sCBA5x22mn199u3b+fmm29u2eaxj30sV155JW9961v56le/Si6XY3JyktHR0WWV1VoEFovFcpJy8cUXc8MNN/D4xz+eG264gW3btuG67rKfx1oEFotlw7PYyL1bbNu2jX379tXf79+/n23btrVsc+qpp3LllVcCkM/nueKKKxgaGlp2WaxFYLFYLKvA+eefz3333ceDDz5ItVrly1/+MhdddFHLNhMTE4RhCMCHP/xhXve613VFFqsI1iG2/ajFcvKTSCS45JJLeO5zn8vZZ5/NS17yEs455xze//73c/XVVwNw/fXXc9ZZZ/GIRzyCI0eO8N73vrcrsqyJnsXnnXee2n4ES/PxP30ND5ODPO+D15JwrY63WBbj3nvv5eyzz15tMZaNdp9HRG5T1fOW2rerMQIR+SPgDYACdwOvBbYCXwZGgduAV6pqtZtybBTemvgqAOVQSSx/PMlisaxTujZsFJFtwFuA81T10YALvBT4KPD3qroLOAa8vlsybFSqQbjaIlgsljVEt/0HCSArIgmgBzgEXAh8xaz/HPCiLsuw4fB8qwgsFkt8uqYIVPUAcDHwEJECmCFyBU2rqm822w9sa7e/iPyuiNwqIreOj493S8x1xd8ND/HmLWN4wckf97FYLCcP3XQNDQMvBHYCpwK9wPPi7q+ql6rqeap63tjYWJekXF98dmiA7/Vk8axryGKxdEA3XUPPBh5U1XFV9YArgacAQ8ZVBLAdONBFGTYkNkZgsVg6oZuK4CHgAhHpkahY9rOA3cB3gd8y27wa6E4VpQ2M5werLYLFYonB6173OjZv3syjH/3otutVlbe85S3s2rWLc889l9tvv70rcnQzRnAzUVD4dqLUUQe4FHgX8HYRuZ8ohfSybsmwUfE8b7VFsFgsMXjNa17Dtddeu+D6a665hvvuu4/77ruPSy+9lN///d/vihxdnUegqh8APjBn8QPAE7t53o2OX62stggWiyUGT3/609mzZ8+C66+66ipe9apXISJccMEFTE9Pc+jQIbZu3bqsctiic+sQ3zv+3qUWy4bkmj+Bw3cv7zFPeQz82vyuY53QrlT1gQMHll0R2DoE65DAsxaBxWKJj7UI1iFetbTaIlgsa4sTHLl3izilqpcDaxGsQ6qV4mqLYLFYloGLLrqIz3/+86gqN910E4ODg8vuFgJrEaxLKn5htUWwWCwxeNnLXsb111/PxMQE27dv50Mf+lA96+9Nb3oTz3/+8/nGN77Brl276Onp4bOf/WxX5LCKYB3iedY1ZLGsBS6//PJF14sIn/zkJ7suh3UNrUOsIrBYLJ1gFcE6xPNtjMBiscTHKoJ1iJ1HYLFYOsEqgnWIF1hFYLFY4mMVwToksIrAYrF0gFUE6xDPtzOLLRZLfKwiWIcEoVUEFstaYKky1D/72c948pOfTDqd5uKLL+6aHFYRrEOCwCoCi2UtsFQZ6pGREf7hH/6Bd77znV2VwyqCdYJqo0+xH1ZXURKLxRKXpz/96YyMjCy4fvPmzZx//vkkk8muymFnFq8TmvQAgVUEFktHfPSWj/KzqZ8t6zEfOfJI3vXEdy3rMbuFtQjWCUGTJgitIrBYLB1gLYJ1QhA2FEGgtlWlxdIJa2Xk3i2sRbBOaHENWUVgsVg6wFoE64QW15BVBBbLmmCpMtSHDx/mvPPOY3Z2Fsdx+NjHPsbu3bsZGBhYVjmsIlgnNLuGQvVXURKLxRKXpcpQn3LKKezfv7/rcljX0DpBw7D+2ioCi8XSCVYRrBOCwEeMeyjAKgKLxRIfqwjWCUEY1C+mWkVgsVg6wCqCdYIGATXnUECwqrJYLGuF5hn5a5kT/RxWEawTAt9HRQAICZfY2mKxZDIZJicn17wyUFUmJyfJZDLHfQybNbROaK4vZC0Ci2Vptm/fzv79+xkfH19tUU6YTCbD9u3bj3t/qwjWCb7fiAtYi8BiWZpkMsnOnTtXW4yTAusaWif4QZNFIGvb1LVYLCvLghaBiFwdY/8pVX3N8oljOV6CoDGbOLAWgcVi6YDFXENnA29YZL0An1xecSzHS4siEKsILBZLfBZTBO9V1RsW21lEPrTM8liOEz9sjhFY15DFYonPgjECVf2PpXaOs41lZQj8hkUQ2hiBxWLpgOMKFovIpcstiOXE8MPmGIFVBBaLJT6LBYsXaqQpwPO7I47leAmCJteQtQgsFksHLBYjGAf2Ej34a6h5v7mbQlk6J2hOH7UWgcVi6YDFFMEDwLNU9aG5K0RkX/dEshwPgTZmEwcSTTsXkUX2sFgslojFYgQfA4YXWPfXcQ4uIkMi8hUR+ZmI3CsiTxaRERH5HxG5z/xf6ByWDmhNH9WWRjUWi8WyGItlDX1SVe9cYN0nYh7/48C1qvpI4LHAvcCfANep6sOB68x7ywnSHCMIAN8qAovFEpOOsoY6yRYSkUHg6cBlAKpaVdVp4IXA58xmnwNe1IkMlvY0WwS+WEVgsVji02n66HkdbLuTKOD8WRH5iYj8s4j0AltU9ZDZ5jCwpUMZLG0INZpNnNAoRuAHdnaxxWKJR6eK4GgH2yaAJwCfVtXHAwXmuIE0KgTedugqIr8rIreKyK3roUxst6nNI0gaReAF1iKwWCzx6EgRqOrzOth8P7BfVW82779CpBiOiMhWAPO/rXJR1UtV9TxVPW9sbKwTMTckoYkRpBB8hCCwPQksFks8llQEIvKY4zmwqh4G9onIWWbRs4DdwNXAq82yVwNXHc/xLa2EptZQUgVfwPOqS+xhsVgsEYs2phGRZwEfBJ52nMf/Q+CLIpIimpfwWiLl8x8i8nqiCWsvOc5jW5oIaooAh7JIS+0hi8ViWYzFSky8HHgH8NzjPbiq3kH7APOzjveYlvY0KwIfCLzK6gpksVjWDItZBJ8FHqmqNlK7Bgi1SRFYi8BisXTAYjGCDwCXiUh2pYSxHD9hGAWHk7j4IvhVaxFYLJZ4LDaz+MPAvwD/tXLiWI6XukUgLgBVv7Sa4lgsljXEosFiVf03ETm4UsJYjp+gZhFIdEk9v7ya4lgsljXEkumjqnrdSghiOTFqrqGUUQRVz1oEFoslHnHmEfy2iPSb138qIleKyBO6L5qlE0JThjrpJAGo2nkEFoslJnFmFr9PVXMi8lTg2URF5D7dXbEsnVKbUJZyUgD41jVksVhiEkcR1GoV/Dpwqar+N5DqnkiW46HWmCZlLAIbI7BYLHGJowgOiMg/Ab8DfENE0jH3s6wgaqqPptw0YBWBxWKJT5wH+kuAbwLPNf0ERoD/21WpLB1TDxa7kbHm+XYegcViicei6aMAqloErmx6fwg4tPAeltVAzTyClJsBwA+sIrBYLPFY0CIQkceIyE0isk9ELm3uLSwit6yMeJa41GIEmUSkCILAZg1ZLJZ4LOYa+jRR5dHHAL8AbhSRh5l1yS7LZemQWvpoOhkpAi+0FoHFYonHYq6hflW91ry+WERuA64VkVeyQFcxy+qhRMHiTDIqDWUtAovFEpel+hEMquoMgKp+V0R+E7iCKGBsOYmoBYvTyR6g0brSYrFYlmIx19BHgbObF6jqXUS9BK5su4dl1ag1r+9J9wLgB1YRWCyWeCxoEajqlxZY/hDwxq5JZDkuQjPvr9cogiC0riGLxRKPrvUstqwsGkZhm0yqpgisRWCxWOKxqCIwPYs/tUKyWE4ANRZBNtUHQKhWEVgslnh0tWexZeWop48aRVDrYWyxWCxLYXsWrxNq6aMp6xqyWCwdYnsWrxNqWUOppFEE1jVksVhiYnsWrxPUzPFL1xVBsNjmFovFUsf2LF4n1C2CRNq8tzECi8USD9uzeJ1QixGkTT+CwCoCi8USkyXLUIuIS9SdbEfz9qr6d90Ty9IpoXEN1foRWNeQxWKJy5KKAPgaUAbuBjPstJx0qIaIar15fWgVgcViiUkcRbBdVc/tuiSWE0JRBHDEwVGtl5ywWCyWpYjTqvIaEfnVrktiOSFqigDABasILBZLbOJYBDcBXxURB/AAAVRVB7oqmaUjVLWu1RMKgVUEFoslJnEUwd8BTwbuVlXbkOYkJSREzNVxtZFOarFYLEsRxzW0D/ipVQInO00WAWJdQxaLJTZxLIIHgOtF5Bqg3gjXpo+eXLTECBQC203UYrHEJI4ieND8pcyf5SSkOUbgIoRiXUMWiyUeSyoCVf3QSghiOTEUrccIEtgYgcViic+CMQIR+eBSO8fZxrIyaFOMwFUhsBaBxWKJyWIWwRtEZHaR9QK8FPjgskpkOS5a5xFIveSExWKxLMViiuAzQP8S+39mGWWxnACtisAhFJs1ZLFY4rGgIliu2IApWncrcEBVXyAiO4EvA6PAbcArVbW6HOfayCiKU5tHgBDYslAWiyUmceYRnChvBe5tev9R4O9VdRdwDHj9Csiw7plrEQRiXUMWiyUeXVUEIrKdqIT1P5v3AlwIfMVs8jngRd2UYaMwN0ZgHUMWiyUu3bYIPgb8MY3y1aPAtGq9a8p+YFuXZdgQtGQNWYvAYrF0wJKKQEQeISLXichPzftzReRPY+z3AuCoqt52PIKJyO+KyK0icuv4+PjxHGJDoTrHNWSzhiwWS0ziWASfAd5NVHkUVb2LKG10KZ4CXCQie4iCwxcCHweGRKQWpN4OHGi3s6peqqrnqep5Y2NjMU63sVFpTChzxSWQBbZTxQtsINlisTSIowh6VPWWOcuWbIirqu9W1e2quoNIcXxHVV8OfBf4LbPZq4GrOpDXsgCRRRA9/SOLoD1/+d/38vD3XkMYWovBYrFExFEEEyLyMIh8DSLyW8ChEzjnu4C3i8j9RDGDy07gWBaDSuNi1iyCdg/7L9y8F4B81Ta3t1gsEXGKzr0ZuBR4pIgcICpA94pOTqKq1wPXm9cPAE/sSErLkrTECIwi8EMl5bT6iHpSCcpelVzZZyCTXHlBLZY1Qq7skXQdMkl3tUXpOnGKzj0APFtEegFHVXPdF8vSKSrNwWIXH/DDkNQcoy9rbupc2QOyKyqjxbKWeMwHv8WuzX18++3/e7VF6TpLKgIRefuc9wAzwG2qekeX5LJ0iCqImhiB4+IreH44r3D4+c7P+KPkNcwWL1gFKS2WtcX9R/OrLcKKEMc1dJ75+5p5/wLgLuBNIvKfqvrX3RLOEp8Wi0AS+CIEvsdcTfDXxfeTcj2+W8gRhWgsFstGJ06weDvwBFV9h6q+A/hfwGbg6cBruiibpQNaGtNIAh/B9+aXcKpIJvqfm1pB6SyWtUupuv7n6cdRBJtpalFJNJ9gi6qW5iy3rCLNJSYS4uIL+L43b7uKE8UFqgWrCCyWOOw7VlxtEbpOHNfQF4GbRaSW7/8bwJdM8Hh31ySzdIQKjXkEToJQBL9anrfdoWSGO9JZvML0SotosaxJHpos8ogtS1XkX9vEyRr6c9O4/ilm0ZtU9Vbz+uVdk8zSEc0WQdKJ0kIrXmnedn+ySXgoOca7CxMrKJ3FsnZ5aGqDWwSml8A9qvpIop4ClpMUpWERJGRhRTBpUqJLFVu/yWJZiObJmDOl+S7W9caiMQJVDYCfi8jpKyRPV7lr/zSf+d4Dqy1GV2gpQ+1G+r1cnT+SyWh0ySsVaxFYLAvhhY16XCVv/QeL48QIhoF7ROQWoFBbqKoXdU2qLnHRJT8A4A1P21mbD7FuCGnMI0gY15DvzY8RpNUBAiq+jRFYLAvhBcoLnRvZr2MUKutiHLwocRTB+7ouxQpTrAb0puN89LVEo+hcwo3mDlT9NoogjCyCajC7cqJZLGsMPwj5eOpTALy9euEqS9N94gSLb1gJQVaCFzvf58XujeTKz1p3iiAEHKMIasHiahuLIGusBk83xoxJi+V48IJGjKBo5xGAiFwgIj8WkbyIVEUkEJE1N5wsVn3+PvVpnu7eTb68/qY/aJNFkDQWgR/Mn1CWMa7PKvMDyRaLJcJvihEUN0CMIM6EskuAlwH3EVUpewPwyW4K1Q32TjYCp/nc+vOPR/MIIhKusQjauYbMQMcXqwgsloXw/IZFEJTXv/Ucq2exqt4PuKoaqOpnged1V6zlZ+9kPc5NeT0qAphnEVz2vV8wW/bmbBfd4FVn/afEWSzHS3PWUF/58CpKsjLEcZQXRSQF3CEif03UlKbbTe+XndlSoxFLJX9sFSXpDi0Tytw0AI7j8dBkkUdvG6xvFxLd4GXHj3oYrLPsKYtlOfD9hiIYrK5/RRDngf5Ks90fEKWPngb8ZjeF6grVAp8YGuTxO06jug7LKzRbBMO9vQCc69xPxW/1b4bGIii5IRXf9i62bFzu2j+NavuWrZ7vcSDhMu04jHhWEaCqe1W1rKqzqvohVX27cRWtKdKlw1w6PIgvQqUwudriLDvNimCgpweAFyR+gEy0XqrAWAR5Z2NkQ1gs7fjB/RNcdMkP+PyP9rZdH/gezzttG8897VT6/fXnQZhLnKyhp4jI/4jIL0TkgdrfSgi3nKRLR+qvi6X1V16huR9BfUIZ4BdbrZ9QohFQwYGS7Vts2aAcmI6SJX7yUPuHfGBSr4uOQyosLmg5rBfixAguA/4IuA1Ys0PITOlo/XWhsv4sguZ5BImRMwHwRdA5ZSbqriFHKJWLQM9KimmxnFQsVD7C9xop5j1aohqEpBPrt3dxnBjBjKpeo6pHVXWy9td1yZaZsCmVsuyt7xhBIhn1HPiXwQG8SqFlu8CYDSVxqBbWf1qcxdKOY4Vojk3Jax8nC5uaOiWcwrpvTrOgRSAiTzAvvysifwNcSVMjGlW9vcuyLSuzYeOhV/LX3Hy4JWlWBCknSh+9NZvhOaX9LduFZsuiI1SK6+97sFjiMGUUwXRx/qRLgNBvWATlZJFCNWBoAeO55jZayxl4i7mG/nbO+/OaXiuwpgpwTIWNh141XH8j4WZFsHNwJ48dOps7p++l7LV+1ihGIAQiFIvrPwhmsbSjpggOHGs/sbLZNVRIlheNpz3iT6/hgjNH+bfXP2l5hVxBFlQEqvrMlRSk2xzTXP11RQuLbLk2CZuLzjkJ/v+zX83v/ehPqPpzXENNrwul1naVVT/k1j1TnLNtkMFsstsiWyyrRk0RTBaqlKoB2VSr/7+5cu9ssrJohp0XKN+/b22XdY+TNfRXIjLU9H5YRP6iu2ItP5WwoeE95pdeWOs0l5gA6M9Gl6watAaLA2lkPxRKrbGSN37+Vv7PP9/M53+4p1tiWuZQ9oIFM1cs3WOy0HAJHc3Nfx5Uvcbvpur6FCrrO0YQJ1j8a6paf2Ko6jHg+d0TqTuE2riQZVm9onMzRY/n/N0N3L1/ZlmPq7T6KLOZaDax57eavgFQ0wWlSqsieHAish6OFW35iZXiHf95Jy/+1A/rI1TLyjBTKPHG5LWkqbYd7Vf9hiLwHZ+r7zy4rlNI4ygCV0TStTcikgXSi2x/UlJTBEmFirN6+fN37p/mvqN53nfVT5f1uJHnv3E5M6kBAKpBQxGoKoFArylFXaq2BotfX/osn0leTE9+z7LKZlmY6+6N5rcsFLS0dIfzCt/jve7neVviiraKwG9q8+q5IZff8lBL4coa1abZ+X6wdmfqx5lH8EXgOhH5rHn/WuBz3ROpOwQaPfwHNEHZqaxanR0FUnjsPrS8GTugDejTAAAgAElEQVShNEb6AJlUlOLgNbnEglAJgV51yeNT8XItx3hR+G0G3QLhxH+zBo2+NUnZpC9aK2zlqPoh+CVIwiky1TY11GtKNy+JIoTkK/MHkKVqwLsSl7NHT2Gq8Cw2D2S6Knu3iNOY5qMicifwbLPoz1X1m90Va/mpWQQDpCg6ZSp+SCa58hNEwokH+EXm1byl+maq/nNJJZanfl9z1hBANhHNJfC1MdL0Q8UX6Jck4FP2GxlFYah8YqSPI8ksT6m2BpEt3WemZC2CE+EDV/2UqaLHJ172+CW3rQYhaSLFm8aj2CYjyA8aiiDvOPRSbms5lLyA3098DYCf5t63ZhVBnGBxL/AtVX0n8BkgLSJrLqWkpgj6JU3RUQpttPuKMPMQAG9NXMlDU8uXvaSA03Q506YCqacNi6CmCHqotbJsmLplP+DrA0mu7+3h1uT6L7J1MtB8Dx4rWIvgRLjpgSl+eH+8zB0/CBmTKEbXQ6Xt7GKvRREIvZQptFEYzUpkPL92G17FGY5+D8iIyDbgWqJqpP/aTaG6QU0R9LlZSo5QKK3ORauUClzXk+V059CyZiKEtGYNJZwESVV8Gg8Y3w8IELJONGppzigqlD0yoelVwPobnT44UWD/sfk+3tVkT1OPjOmSVQTHi6qy71iRyUKVmRguNi9QNhEpglNkqu1IPwii50NGXAri0Celti6k5n3Hc+tbEYiqFoH/D/i0qv42cE53xVp+ApNB3+f0UHCEUmF5s3bisrd4P2/bMsZfjQ5TnT64bMeNsoZaL2dKBU+bFYGHL9ArNWuhqexGsUDRiVRJs/JYLzzz4ut56ke/u9pitDCZbyhcGyw+fiYLjcyfByaWniwahMomYxFskwmKbbwDNYtg2M2Sdxz6KLX1IpTLjaDyzPTadanGUgQi8mTg5cB/m2VrrvpSSIirSjbZS8FxqKySIqhWo/Ne2d9HMHtkia3jE80jaA1+p1Xwady8ge8RIPQYRVANGzdxPj9FyYluh6rYqqQrQRA2ovvTNlh83Dw01bD0Hhhf2t3qBWFdEfRLibDNDHs/jBTzULKPvCNclX4/QZvtKqU8h12XGUeQ3PIN7FaaOIrgbcC7ga+q6j0iciZwcg2tYhBqgKvQl+rHEyGfXx3tXTZ1jkIRpgvLpwgi11CrIkip4EnDdPW9Kr5AUiK3UbkpoyiXb9QRbN5nPdDqiz95Rt5+kyI4Zi2C42LvZIFX/8st9fe1uTCLEYQN1xCAFOfHFgLz2xjpGaPsOHiAzB6Yt51XzPGc07fxa9u3kVjPikBVb1DVi1T1o+b9A6r6lu6LtryEhCSAvnSUX59vc/FXgorfSNmcLC63a2iuInDwmopK+H6FUISEkyCr0hJIzhca38d6UwSHZhqWzy+O5BbZcmUJTF/cdMJhxsYIjour7jhIrhwp+k19qVhxID8MGZUZrh86nRAIy/PviSCMrsdw7ylA1JcgbNPE3i9HA7uc65AptU+y2DNR4Mrb97ddd7KwWPXRj6nq20Tka8C8KXWqelFXJVtmQg1xFQbMjNt8cXUqaVe0caMWq8vnnormEcxRBOLiSeMBU6lGfk9XEmRwQcqUvYBM0iVXanwfVdF11c/4wHQjFvKLo3medOboKkrTwA+j/PSzekscK/attjhrkoOmwczLnng6d+2fjqVQg2qJvWnlD4fhE5UsWpn/gA/DKrgwnI3ulZwjhJX5CsMrNyyQTLl9w6sXf+oHHCt6vODcU5ctXXy5WWwewb+Z/xevhCDdJiTEAfp7hgEollfHNeSFTSmby6gI5s4sBkjh4jXNMvNMIa2EmyAradSZ5shsmTNGeymWI/+no1B1dNXmWXSDg9MlzpSDVEnw0OTJU3AwCJU/S/wrr6x8m+fwpdUWZ03ywESB83cM8+EXPpJXXfajeIrAq7A3GT36HkwlwGsz0jfzb4YzkSIoOA60UwSlxsRQz2tfM6o2WXA8X2HbUHZJ+VaDBdWTqt5m/t8A7AZ2GzfRDWbZmiIkJKEw3Btd2GJllYLFjZYOVLzleyg1dyirkSJB1YkUgapy5W1Rh1FXEowlB5hIwPhUVG+oWIlu6CF1qIiuq37GB6dLfCf9Tm5Mv41CFz/X7oOz3NeB6ykIlVcmvg1Ab/XoEltb2vHgRIGdm3rhqj/gXdN/HlMRVDmYiBTBQ4kETnW+Iqi7hoxF8FvbtlKtzm9oVWoaUM7q4sUDD8+cvMUuF7VTROSDIjIB/Bz4hYiMi8j74xxYRE4Tke+KyG4RuUdE3mqWj5geyPeZ/8Mn/jGWJrIIhKH+TQBUvFVSBM2ummD5+iIoMj9GIEmqovhByI9+Ocm3dkd+yoSbZGvPFg4kEswc2QNAySjGIdKUHaFYWT/By+YfYLncvR/j8//h+zzn778Xe/vM7IP11z3++uua121yZY/xXIWdm/rgwG08rHwPszEUQehXOWTaTu5LJnHaWASBSbve3LulvmzWnx9XbLbqc878sjFhqDxMDvAy9zqOzrbee7myd9IkLyyoCETk7cBTgPNVdURVh4EnAU8RkT+KcWwfeIeqPgq4AHiziDwK+BPgOlV9OHCded91Aglxgb5spAiayyssN0dny+ybah+0qojHYBh97ZVgeSY4qWoUI5hzOdNOiopA0QtwCXhr4isAJCTJGUM7ybkOx8Z/HslirJMhp4eSCOXi+mneE5Zm+cTQIJcN9jOQf6Ar5zieyWrZQiOA2B9Mr+vqlt2glja6YyQNx/aQCQtkyktn4gW+xyFjEexLJHD9+deuVpvsnNFzePlZLwXAC+dbe1Wv8fCfduYf52iuwivcb/Ph5GWMH2u1GN73Xz/l8X/+P7FnRHeTxSyCVwIvU9X6sEVVHwBeAbxqqQOr6qFaO0tVzQH3AtuAF9IoWvc54EXHJ3pnhCiuCn09YwBU/O75ii/48HU87a/nZ9iGoVJ1AoY1KvHghcszOlVtnzWUdtNUHKFYLDF46Eae6t4FRBbBzrFHAjA9+wsAKkYxDiUHKDkOpfz6GaH2lg9y6fAgHxsZZltxd1fO8cP7O08+0KAxet0ks1T8tVu9cjWouS83BZNgXDk7wn2UF2hIX6PZIjiUcMHLoaq856t3c/MD0XUMzPybtJvmeTujAox+Gwu+2pQFWJH5o/u9kwW2SxRELk881LLupgcit9Lnf7R3iU/afRZTBElVnaeqVHUc6KjWkIjsAB4P3AxsUdVDZtVhYMsC+/yuiNwqIreOj7ePxneCmmBxT6oXgErYvXIDtfTwuSO8O/ZPU3aUITdDUpWqLo8iCMOwbYwg7WYpi1AqzDLuVXnO6dsASDhJtm9+NAD5ajQqrX0fQ6mRaHlh/TRL6ak00nRHy3u6co479jcUZ+xG502KYITZdd8gfbnxjOLsKzYepA+XA0vGCXyvwsFEgl4njYpQYoaJfJUv3fwQf3D5T5gqVOsxgqSbrD8zAp3f1rJZEXiOhzenFPXh2TKnygRFEfRYqyJ4k/c5fpT+AzZN39HBp+4OiymCxZxXsR1bItIHXAG8TVVbnGgaPSnb2sOqeqmqnqeq542NjcU93YIEKC5CbyK6qMv1EF6MiXzr13Tt3YcoOspoupdMKPNq+hwrVLltb+cP4DAMCUXmlZjIJLJmdD/DvaU99eWum2Tb4A4ApsLoAVY11slIOgrZ5Eurk17bDVyv4YvP0R1Lp3kUGrfmjAaNiW6jMkvZt4qgE7xQ6aXE6XdE7dVDSbArhiKoeAWKjsPpyei5ElDgiPHfj+cqvOfKu/GMayjpJOlLRqm9PvMVQXMWYMUN5pWhqPohPxjK84LtW3FyDUVQ9gKeFNzOVpniuTP/2elHX3YWUwSPFZHZNn854DFxDm6qlF4BfFFVrzSLj4jIVrN+K7Ai6RI115DruKRDbZlM1Y4D06VYsxTb8Znkxfw4/fvsm+M3PnjvTRQcYTA7TEYFb05Nn3ddcRe/+ekfdly8KghrD5BWi6DH3MC5/ARuk9sh6aYYSA2QVKVgRjlVrZBUrc+zKKyjxvYBDYtgwunOhDI/aIxn2rU+bIeGjYfGJrEWQad4fshL3OvpHY9G1KWBnWySmSUVQa0N5XCyH4BACvVr9oyzxrj2nsOUfQ9XFUccepPR4NFn/u+yVqaljwRlJ6xPbqshlRkeSAvjiQRBZU99+aGZMg9my3x0ZIhxd/VrFC2WPuqq6kCbv35VXdI1JJHD+jLgXlX9u6ZVVwOvNq9fDVx1Ih8gLiFad51kdf5ofC5P+ch3eObF1x/XuZ7j3s6YzLBvTs76hcVryTkOfZsfRVadeTV9zpn8Fjem38I37+5sFmJgRpbOHIugNx3d6PniJOo3f94EIkImBM/kS3taJRtCXyaaeV2bV7AeqEjDuplIdCk25Je5LPk3/HHiyxyNq8gDj6/39vCC07bTz0y9SY0lHl4Q8mgnsva+9Lz3cll/mmHJLZk5VCsoN5KM7vWASr372K8+6hSzLCBpJmjWFUGbFreeGUhtkh4KDhTmZKVl8vs5aOIRPofqyw9N5fjqkPCFwQGuHSiueqJAN6e5PYUo4HyhiNxh/p4PfAR4jojcR9Ts5iNdlKFOKJFrCGj7EO4GU0f2tbx/lN5J2XEY6NlMGpfKnFIOv1P4VzLuFA/d3lnfn9CYsXNrDfWb0X2uNNUSHE+4kR7PNCnEinhkVBjIGougsrwd1FaTabdMQiEbCkWnO6Ucnj35JZ7l/oQ3uV+Lb9GFPu/evIm9CQdN5trWxbcsTDUIOUf2UjjjQj7883/jM8lZhsgtaRF4fnR9hmqKwK1yz75J7ky/gadORZl1KkE9EJpwEqQVAsebZ7X5xrOwKTFAznEo5Vpdj9niQQ6YDKVS08j/2NF9HDEKouAGbbufrSRdUwSqeqOqiqqeq6qPM3/fUNVJVX2Wqj5cVZ+tqitiFwVNFkEGl8oSiuDS5N9yTepdx3Wuj44M8dpTNuMd+UXL8nsz0Q36uM2PI02CqrSOAD8yOsYzz9jOw6b/p6Pz+QtYBH2ZIQAK5WNUm3qwnjEaLY8UQbSvR0BahcHeKFhcXkeKIOdUGQmFXhx8x+9Kb9kRLxrtHWUotmvoUNDIxdBEkYpVBB0RVEvskgN4Yw1PdS5VXDpYbNpQDqejQU/gVjm0/wEGpcjpN38o2khCks0d/3DxJGCy0KrkPaqIwqbUELOuQznfakmHlWmOGkUwkyzXY0nF8b0cNopg1lWmC6vby+DkLHzRBUIU13zcDEkqsvjD4Ffd2zjb2bfoNgvxhcEBbs1mSM3cX1/mBSE/7nHIqssTtjyBjKQoO9piEv4yET2sZzI/78hUDMKaRdB6OYfMLOpyZYZqU8cl1yiMtLp1yyhASajQ3xvNs6h6J09xthPFc3x61CGDi+8EFLvwwHWMVbZJZpguxFME09r4jvMJ31oEHZKdfYCEhIyP7Kgv+2UmYLY4P6jbTE0RDKUHcYGq41OdbPzW70q/noyUSDT9nnpI4LsBU3MmgAVSJaPCYHqAWcehkmtVBFNh4/14MqyXG/fye+pl36dch5ljqzuXYOMoAmm2CJJUnHDRh+1lg/18ZGT4hHx3YaVhERQrPndnEpzNCEknScZJUxZacse3etG5ftQbdFTiITDZJnOzhgb7otF9yZvFa+o9kKhGryP3VHR+X0ISCNmeSHk0p8WtZVQVT0LSuGRIUHVCisvYGa6GY9INE4RoIV7GVRA0Hip5Vym3aYVoWRgxpSEOOo3v8Z5MimpucSeDZ0pMpxMZ+iRJ1Qk4RSZ476YRfpJOMSAlUk6ZhDRKsfU60eBxskkRVPyAwPHJ4DCUGcYXoVhsTXWf9E3fA01wNAkzhei351UixXOGO8Ss6zIzNb/E9UqyYRRBAA2LwElTclh0As/HRob54mD/cU3ySZtdck3ZKqVSgZzjMOj2GhkylBzIFZvKH5jR+f6Uw9Rs/Jm9obEI5s4jGOozs6i9PF7QMD2Tm3ZFckqy7p4KCHFx6DO1VSrByVOc7USo+CFVJyRNgqyToupo296zJ4qjDeXixuwzUZu9CjDpOlTL6+M7XylCU0RxTzUaTQ+4GQ67LuESitg3CjidyDDgpCi7IcPJQ1zd38erTo2CxZ4ISadJEbhpKo4y1eT2y5d9PMcnS6JewyxfaFUEeY2u6Y7EEMccl/x0JGsxiJIlH963E4Cp6dWdVLZxFIFovbl71snE7ltcKHceXDzV7DKZyBGa2WXlwixFR8i6PZEMyR4KjkN+pnHTlpzoYTLpuhybOERcajEC5mYNZc0s6qBQHwVFK6LlaSL3FDRcQ32pKNOo2sUJdytJqRpQFSUjSbJOmrIoxXI3FIHHN3uyXJ/NkirFy4huVgQTrovfpt69ZWHUr3JjNsPFD/47KSfFYwZ2Mum6sESJed+4SVPJDP3JKNtnc2pPff3NmTQVEVJOqr6sP9VH0RHyM40Hfb7iU3WUrJNkpG8zAIVK67l9YyluSg0x4zoUZqJ7Y9ooiEeNnAXA5Ozq9ivYMIqgeeZtNtFDQYRifumAaLHU+QMxYebITaYq9abkpeIMJcchm4gUQV9yAF+E2Vx0YwShUjAPZU+EyYn4NXGCBSyCnkz0UK8ExXqaKEDR1FZJOynKEp3bF8XFJZvI4qhSWYEJdytB2Q+oOJBxUvS4GYqOUCwt/8jbDX3euWWMPzxljGwl3kz4miLISpJJ1yFsUxffsjDiV/haX2Rhf+LCT3BKzxYmXBentLhrqDZrOJPMMpAeIuc4jKYaMYI3bN3CvakUg8lGj4j+zAB5R6hON5rP5IoVio7S52QZ7IsKJJQqreeuXePNPZuj3/ts5CWYkej39ZhTzjXHij/w6wYbRhEEAq5EUfqeRC8VxyGfW9qXW46hLOadyyiC6UTIuJmxWGsFmTU5yX3pKHPn2Ex0A+SLFQqukAmjh/n49J7459P2WUM9RulUwzJ+UxP72kgn7aQpmUqjvigJcRARelTw2kyeWYuUqgFlidyB2USWouNQLS5/5dl8U52ZPi9ejKCW9rslNcik67btgGVZGN8r8f1slhdueya/su1XGO3dypTr4FaWiBEYN2k6maE/O8qs45BLRNfvRbui0mcTCZdRk3UH0JcejrqUNfUZL+ammXEdBhI9DPRGFkHZb723AtP7eOtAVN5lJhfFAqrqkQ6VLSORmzZfscHiFaE5RtCbijT9TH5pE34233mdI98MzI8kHKanohsnZ1pj1lwvg6ZBzrSRYXZmklnHYZtED+9jufimopqWh3ODxQknQTpUqlrBV4/eUPn4Mz/OeaecB0DG7aHkOBTzs6YER6QosyptC2itRUqVMiVHyLgZ+pK9FEWoFpc/EL432fi+0m3KFbejpsC3ZkaYSLiEVRsj6ITDwRFyrsNTTzkfgLGB0whFUH/xGE2txHQm1cNA7xZyjsPBZIJRJ80rzn5Ffbsxk0EH0NczSs5xoCn+UykcY9ZxGEgO0F+PrbXeW7VrfOrwGQDki4fry5MKI/2nAlAIVrfI44ZRBKE0Rsx9Jn94prD0Qz5X6FxT18KG467L7HgUBMqXojSyXjNzd9hUQc2ZLINj04fwRdiWikYWszHK6dbPZ4qXOW1aS0alLKp4+PSEwoWnX9hYZyyGmdlxfIGEUQQ9TWmla51KYZaS45Bxs/Sl+yg7DsX88s+afjDdCBZX3dl5xcfaEZo7ZTQzTM5x0IpVBJ1QK5Q4YkrLb+rbCkCwxNSkmt8+lcwymB1l2nXYn0iwrfcURrONNqajPY16mJv6T8MTISg3snvK+UlmHIehzCADqQEjU+s1DEwSweZBowhMYNsnmrA2kBrAUaUUrO68nQ2jCHwaFkFfvZ7O0nPZZo+jyX1g2kOGIoxP3gdAsRJp/H5jbo4MmOCS6XA0ORNZANt6TgcgV41f9M0Pa+mj81tLRpPGPHzx61Pma/QYN9WxQqQIaq6zaJ7F+shpr13jnmQPQz1mAJA78Wq2c5lwmxRBMj+v5kw7QvOQGEoPUxAhbNMgxbIwvnHxpMx9vMkohNCdobpItl9g4mXJRIadgzvxRbg5m+HUkbMYTg/XH4rDRrEAbB6Ifpde0BigFfJHCUUYyo7QZ7wM1TkVSmsWwdjADgBKpp1lQEBCo2ZS/aFQaVPQbiXZEIpAVQkEHPOgG8iaUssx6ukUSp2PHn2BYdN85lg+qoVSawVZcwmNDRmT0HQ4ms5FJuP2vtNwVDlWnpxXyXAhwgWKzgFk1aWMT5WA1Jz1WeOmyhcm8QQStRiKJCk7ShCubv2T5aBWRbUn2VsvuTETM8+/E3xRUkbRlpJlcjGyzWo17wczw6gIwTL2sN4IBGF7RYBbYHaR779eYjqR5eFDD68vP3vTObiOy4hEMbSRgdPq67aYTmWeNu6dvMkOG+0bI+kkyYZQkdYkixAfV5WRnki2knEd+c21jNSl3KaO0UqyIRRBEEaKoOYDHzT59bVR+lyaJ5EttM1CREpHGFPjdilHWQLlqlEEJt94aCAabZRNh6NcKRqljvRvYYgEWXeCb9wVL04QBpEicNpYBP1uEs/xqWpAUueUoKgVpSsdixQBUd50RjKUupRvv9IUay65dD+96UbtpeXGF2VTmCShkE9UO7MIzAOsusrugbVGrc5P2ozGa4qgmigvWmaiNkpPJrLsHNxZX/64sccBkasOGu5biLJ+AHx3tl4molCOvAWjJnW0Rx0qcyoKBxqQ0KhwnatKmch1FJgJnAB9pCg567TW0MmEH2oULDZTukf6owtcqrb/4flNI+FSh6O0IFR8YJOYEtAmcFgyRd+G+iJFUAsulUzXo4JROEO9m9nRv409acX75Y3xzqk119B8i2Cgbwx1PFJOiRStiqLXuKlyxenIdeZEZbZ63CzFmPMsTnZK5eh77U0P0GN+4OUuVFb1UVI4DEmSkuvH6p0bECCqDJh7wfOtIugE37h4Uslo0JVJZMiog+dWYyqCDD1mX4BHjT4KgE3DUSbPaKYRLxgzc3K8RKk+u7hi+kwP9kaT0HpwKc95oAcEJIh+mwOhQ5nIYvAJSRiLoM/JUHTCVS1DviEUQRAE+CL1EfPoQGTmFReop+P5jYtZ7tBv6wchvkCfpEmHSlEjRVI2/YkHzeghm8jiqtYDXqXa+p4RnrHrIu5Np0gd+36sc6oubBGMjp3N/kSSvOOQnKMIBoybKleawhOpT6nvSfSQF4dibm23q/zleJ4rbo/KfPRlBuk1P9jA64IiECWBy4CTpuQGzMawCAINSQB9xlXprZPZ3CtFzcWTTqTry4acFGXXX0IR1LqPRS6gp5z6FE7vP51MIgNQDxiPmOsSbZtkiCTFRJWpfE0RRL/tgb7ovuolRckJWjwKASEJ87ZfXUqmHEYgWq9l1Ov2knfhWG71rv8GUQS+SR81PvBa32Kv/QjMqzRSATttcu8Hfv2hOhy6FMU0fjEP+tooXEToDaFsJm5V/Wi7/p5BnrnjuQDs1ofmHr4ttcY07bKGnrvrIqZdh/tSKRxtVQSjJmCdL40TipAwFkFfqp+qI0zPrEjPoBbyFZ+fPLQ8D+qfH86RMt//QM8QPaaOEl1I1fPND3vAzVJwIFdceiJiiHEbpKOMEy+0iqATAuOGSbtNisDtIecq+fzCv9u6RWDKsf/jc/6Rr7/46/X1ZwycwUhmpD4Pp8Ymt5eZhDI5Gw0gy+Z6DRr/f4+kKTnaUicsNK4hgD7SFJ2AMNTIIjCP34HUIDOOw+yxlf+91dgQiqBaLaMiuKZ2SF+qD1Gl3KYZNUC5qeZLpzV3vEoZv6YISJN3PfwgrM/UzTaZoj0qVMzEraoJfA30DLFjcAepUJkmXr57uEjW0K+c+iuMSjTSEb81+Lt5dDvQmMySNIqgFlStTXZbSf7gS7fz4k/9MHagfDG8IMR1IkUw0j9adwMkyC3L8Wuoaj39dijVx4zrUI2RohoQ4mqj8YnXpifuSlD1Qy665Eb+/cfxBh4nC7WRfcptlIIYTg0w5bpUZhfODKu5UpuLyjW7VV99zqv5ym98ZZ6rdSw1xFHXpTgVJXaUTKZPLXW0181ScIR8k0u1JRbgZMm7Sq7s1wcOAIOZEUqOw/QqFp7bGIrAiy5MLT3SdVz6Q6FI+1Fbqdp4AFfCzkotVKvmXE6CIbePGTdkKl+hElbIhNoy+zerLmUxZWlNKlxPNrIY+tWhFLPMQ1i3COZfzoST4IKRcwDYJK0PmuF+owiMieuaGcdDZuRcy2RaSe7cF43WlyNQXfFDXCf6DnvSw/UfrOMWmMgvX/wjCBVPICkuQ+kBZpx4iiAkekjUO2CtkiI4Vqxy1/4Z3nXF3dyxb+24A2sWQXNNoE3ZYaZcBy+3cNp3QICriuvMHzhBZGGM9czvk35K72aOJlwq09EAqaJVkqp1l1Jvopec45CfbVz7QANcM/4aSPQx6woz+QI+Wp+3M2xqf61m4bkNogiih4HbNAIYUIeStH/QlposgnKHiqDsFcy5koykhplIOBydHKcQlOmZk43ZQ5Kyydev9VBOm5uqL3QpO/Fm99b7EbRRBADnn/EMAArpuemjvSRUKRJZRknzgxrtj2IoucLKK4JaQL+wDKWiK35InxspuWzPCEPpIRwFL1FmIr98M6f9sMkiSA8x6zgEMXo+hxL5j2vN0T0qq9KysNzUB+F1//pjStWAo7nySd9D2W/zQN+U3cQxxyXIj7NnosAvjsy3qgP16+6aTtg2sJ1jrkt5Jhq5e3hkmqYr9Kf6ybkOszMNJdTcGbHhAjrYMoFzbCCKMUzPWougq9S6c7lNZWX7NEVhgdmz5aap/s11/ONQMbX+E06CsZ7NlByHb9z0Yxwt0EOiZdseibIFwlDxtEoq1Lo52icpik68H2K4SIwA4PGnPx2AB/pGWpaLCH0h5I1CrAXPtoxElkItpXUlSbrRZ8gvR4XQqTs5MPQgaYS+vq24jsuwpKi41UXzzGfKrmgAACAASURBVDvF8wOqIiScBIPZEcqOg1dqehiEysXf/DkHp+dMNiLEbbIIQsc7rrLnJ0qtV/KvP2YrU4Uqt+yZ4ol/eR3v/erdKy5LJwTq1+du1Bju3ULVEaqlIzzj4uv51b//3vz9COf8EuNRmx2cK+wBwFOPdNP5h0yW0eRUI+27OTtoODNKKMLE5D6qTfN2Thkxxy3Fryaw3GwIReCZxu1ukw+9X7Lk3fbNaSrVhsuoKp0pgprSSTgJThmMZiMWHryaw2mPrf1bW7YdcPuYcSFXqpqbqrGuV7IUnDBWqYJG+mh7U3fnwE5+/cxf58+f9uF56/pCh5xTC561TqQpVpZ/4tVSjDo5nun8hFzlxB/Uh2d+yIOpJO951Bvr5vuI20MuEZIrLJ8bxveqxjWUYMDkspeaKpDeP57nku/ez5u/dHvLfrWMkpSbIqkQOvPnH6hq162EmkXw8C2RZfL2f78DgBt+sfIDgU4IJCA556sZ7osmak7MNB7Gc2cZ14L0nbJl+GEA5CvR3CBfWidpDvVFLp6ZpjphAY0YQa1U9dTs/sgiMB6KTUPmOVFevcJzG0IRVE1ruhaLwO1l1oXZ0nwXQdlrKIIZ128xnZeiYvZNOClOGzsTgE3Zm3ggleR5Z7+0Zduh1DDTrsvR8f346pNqujn73d7Inxij7aEuEiOAaOT/kad9hCdufeK8dT3qMmOcmEmTfTFkptYXV2GC00crf8VnU39DOX/is2w9U277cWPn1pdtSg0w4bpUFwkmdorvVSKLQJIM1magVhs/aj+Ivt+5FkHUPtWURsclcLx5M5Jf8Ikb+dT1vzxu2e4/mmPHn/w39x5a+FrW7u9dmyNFUMuTD1ZACZ0IUb2eVotgxMS9pouHea5zCy92vs+DE3Pq/5jyDp2yuS+qIFr0o3vHIyDVNElzi6kWMJNvNKQKm/qgbDYDwenC4frAAWDUpJ8W/NWLz2wIReD5tQBusr5sMDnItONwbGq+OVY17p3eECYScORY/GqVtcB0wkmyfeyRAHxzQHBU+dVdL2zZdsRMUjk0fh8erWZuf7KfGcdhZnrplLJwkXkES9FDikk3ug1qimAgPYio1qfDrySnh9Foqlw8cSVUa66TNVlQAJt7NjHpOgS55TPDfa9KVYSkJBg01SR9vzF7OV/xOVd+SbHUqgj8Jv9xjyTxHZ/8nGym+47m+c7Pjj+t8Lp7o32vuG3hWeplP0QI2TqYqS/b1Jdmuui1tGY82QgknOcaGjGB14Sb459SH+PvU5/m53PiBMftGqqVmjZzgzwJSDU9Qrea3smz5cYgI0DrFsEpo9HIf7p4xLiGoufRUHoIUaUQrl572A2iCIxrqFkRZEZREQ6PPzhv+9qofkuYZjLhcuRg/BGZZwLTCSfB2FA0fX1vMslw6DDUVN8cYMyYsUem9uDRWhRuIDMSFa2bWDqlLzBlqBeyCBYjK2nUxBZSxn3iOi59oVBeIKuqm4RGmXmFE7cI/CB68GaaFUHvFibdKJi4XARGESScJIPGNRRoQ/5gcg9Xp9/Hu/VfWvYLCXHNiLLHSVF1QnJNE6GCUKn6IT89MBPLRdiOgWx0zy82wapaLvBg5hVsv+uT9FAmic9n+z/NC50bue/IyVsIrzkXv0atzMQpib1MOA6HXZf9+1qzcYImv30n9Cf7SatQdfNU/ABPQpJNc3NGTAyh0FQwMlL2kYxjxgWUKx2N5hqZ51HCSTAQCoVV+L3V2BCKwPdro/TGOGDUmPBH2zSAqfn5T0lED5D9R++Jfa6KmRiWdNP0pfrJmHIVw5qat+0pozsAmMjtxydsMXNrdU4mpqPOSZP5yrzRYo2aRSBO55ez1+mtv64pAoA+dSitQnOa0NyS3jJYBJ7J+MpkGkHy0b6tBCKUigcX2q1jql4JX4Skk2TQ1DMKtCF/pRiZ/E9zWoOvUf0rkxzgZsm7QjHfcA8UTQptxQ/5+eHjGy26RskvFhxXoxS33HYxuzOv488Sn+Uxx77NOxP/yS/HT2JFICHJOY+w/9feeYfHVZ2J+/3uVI2aJUuWJUtyw4VmG2wwBBJMC81sEkIJ4UeWkE5gYRfCbgIbyIZNxYS6NIc0eg0d0+wQsDHuDeMmbMu2bElWb1PP749zRzOjmZFmRsW2fN/n0aM7t51z2/nO+c5XwmEhxNHG6WPLObtyDIHdq+OOsyUI0NgXIsJIXHTYvdS1evGJwhlliVhofrMdwahOABFBUJBrqpZ8dQTN9yVMvrLTbhy4kC6HhSDo8oetYiI3vmSEfij1jdVx+4fnFEo9usdeU78p5bL8ZmJsh82JiHRHIS2IanDDVBbrmCbNnfvMaISRxxE2KdtvmqrNvONdzrhzUcIyI1ZD6auGsu2RdHxOR7QgsB+QQFjhEUGwawDmCEKxQckAisIT4b6BM40NqxIdhpMCM56R3+jo1q/725tZ7XLiko6YBlmrDcz4V848Gg0bvqisedHmm6sytO8PC5PeRgQ9M6Ndbl8IQI50srvx4M1dHRAVF0jRYXOQGzL4OCvyLrub1sbsE87PnQmj7HnUOxS1Dc1aEEQpmTwOjxlWJnI/o0NJeJzZeEKKNvSztEtUTmTlojVFK8HB4LAQBG3esFdvxBV9XKnW3ze0xQsCv5ncujxfWwlU12+myx+kuqHvj8IfniMwLXDylS6z0FEQt++YYh0Ct9m3H5+omFhA5eZEc1PUxFNta+IeQ2SOIP3HmeOMqE1cjqzIenHT0SNuylAQND9Q1dV/fWlA+XCFFEa0nXmeHr53pZhFLBU6zYQyDpuTLHsWWUrw2bws3rafZ5dVU9e+nSvLRnN/kYvdjZF5gqCo7jzThe5CGmwGgbbI3EJ0qILVOzMTBO2+IG68NHcmF+ohbyvNhtDy5bsjK3NLKZA2uhKoTg8WdHrVBBF3Qw5WuiOCwBvaGB//J4MRAUBZVgl77Haa9+3Aa4BTHDHb80IGHUQMPAJRUY8BCkM2mmzxHdM8w0OLLdRrHoXBZFgLgltfXM1P/rqITjMBfZYr8nJUlhwLQFNXfBiF8IigLH8idgU++ZR/u+12FtxzCYE+dLVhIeIwA2GNMMKx0kfF7et255EdCtEWbNb6xqgXptScaG7t7FuFEeqeI0h/RJDviURYjBYEubYsWm3Q2DFw9vapUGMzeDEnm9AACQJ3D0FWZE74eUMDF3gurEoMhzrIx0GXzc+v39zIz19ZT4eZE2Gxxx1jwRKIUhuMyi6i2WbDF5URLywI7IawZldmgqBg/yo+c3+bCa3Lku7T2dXAqWMruKPpQ7xmoLWHJhzHYrebnIaD15cg0KPzFGa0Xb/HpTb97bW49sc4EAaj9PbpUjGikkabjf37NuOVWK9mgFxlpyMqzWt0QiyAAtzUOfT36jCiYiTZ82mwGTQ2HxjLoWEtCKbvvJYj991Ih5llLC8n0ujluPPIDimag/Gx6cPqnVxXHhcWzGBhrrBm0ovcM2En1bvjRxCxx8aOCAqdeoK4NL8i4f4FQaE51IZfVLc5GUBxXhmGUrQE6vvsJXTPEWQwIigrjMRjdzkj6qs8Zz4thrC3dmgDYb2VY+e24pH4fP3PGeDHH2dVEp5M7EoxjlMqhI0Lwo1Cvrhps4c4fu9zzOMu6kxT1RbDoGpP5H4GhW61QVHYhLAz0jHp9Ac4z1jKV0ob2VbXlpETXGGTbshP8n6UdMK5zbRyebPuE04tyeUfWW4eaFnPtaOLyW+rSrvMoUApZX4z8YKg3NDPfHb+EZTjYa+7i537I887GGWtlS4TzQ7anvpP8YrgjGrMAfLETZstMpIOSqz/UqEtjxqH/s4dtkjHdIRbOyLW1G7PqF79ZVgLgldzFS/nd9LSoU0FSwoqY7aPDNpoIT6oXMDs1bscHn58xm9xGQ46zInYdVsW9lpmwLRQCk+8TirSZU4Zc2TC/YtDbhqMDnwCjih9o8NwUBASWmiJsS1P5NMQniNIFjulNyaXT+9edjsjI4KinFI6DYN9uz9N+5z9ocX8iL0DYFMdIBDj+Qk6xLZLQWeajoK94fVHDAQACuzaUfAi2z8511jW/f51GgaNu7d0Hxckos4Lp0WM9i7t8AV50HkP8/Zfwzhq+HRP+hPotcrg5uKRuI2mGLVUNF3+yL3uEsX1JWbsfRHyvbsOykx1wZCK+2a6MUOkTBk9iwmeUqqcNmqrI5Z/gai5mXQZU3wMAC1tm+kyjLgRQZ7NQ6tNG3eAqRqK6uCNdBZ1L0cHyxudp/0ftu9en1G9+suwFgQlzmJqHEJLm55wHVUwLmZ7IVk02+LtpMMjArfTQ0luGS9d9AYXlp0JwJa9yYfYAIFgWBDoRuGIKXOxYTCu8tTEdbQXss+h8Ao4eugbR4ZctBpdtHYFmC5bGSc1Ce26VT/mCMaXRQuCyIhgXPFkAHbVro07ZjBpNy+hKzQA5qPEhyAQEUaE7HQYvgGb/wiHEA/HiRqTXch+m43x9u3sctjIIxJDJlgbEazR8WYKzYagNcojubMromu+zLaQnSnMUfVkjdrImznZfFzYkDDuDoDXTIjzwMl36DqKMMZ0ngpk1VCXZG7qQOIPKjPcuyNuW6uZ/a+05FiOLJpKtd1O056IiiucOyITxpgm4V1miHiXLStme4G7iAabwZ492mTV32OOYJTpZwKxVnpHjNGq6po6SxAMOOW543QvLKgf2oge6pkCI596u6Kzh1lmuDHPcumwxaU5pdw+59fYlGJPe+8+Bd1WQ6YgOG3c2bz59bcYbSZF6UlpdiVdhtBiM+IEQaGRQ6MtRH2bl5ddP2eR68aEH2V4jkAyGBG4nZGw2OHrBZhsCoja5q1pnzNTlFJ0GvpafKr/qhu/hGIcfsLki4t2eyBmMrY/hOeUws+82FNEk2Fw0ZjRzK0oY4Qt0rjneNd0j+qCUV6nhWFv7qiRUKC9kWVuF9scdr5orGNXBoLAE9CjyU6bly21iU1BvWY49mJTGAHcMPMGAOrdLexuOvgsh5o6uvAJuIx4QfAfM/+D88efz6ljTmVK2SyUCPWNqwCdOCoIOIxMXMr0pH52SGhy6ffTFdWYA5SPGI9fhB07VxAyMyNGT2hXjpzYveyKstKbOvYkAGrbDszk/LAWBBNK9DCuxtFKVkjFPbTirFE022xsr45t7PyhsCDI7V7ndGRR4RdqVO+OSIFQeESgewoiQmlOadL9J4w6KlJGj2HmSOdIau3Ctl0RU8f6RIKgHyOCaNxR11sxSvdQGr1DFxHRFwzRYeheenuwHW+gfw21DgEQLxwLbdk02xQ76gcmEUxkRKCf+ciiKSgRau26sdnhijwXZ9Z2qup0uQGku5EI278HVWP3SCXU3sDVpSV8tbwMX1YNjXXp+z6EzOQpXfYAW2oiE+TRo6GwIMjJKuaE0ScAcGbFmZQY2VS7QlTXHLhgaMnY19CCVwSPI94/Z1z+OH77pd/isrmYNHomAE1ePdfR1OEjICHsNnfccakgIlQaeWxy62fq7pG8ZorZgdpdtw5/MIg/Kg8KwJSxs7qXw20EQHHuGNwhRUNg6CP+wjAXBEdPOBmATW47eQnmySpHavPN9VsXxawPmAkvstw5MesryKfa0YXPn9wULzyaiJb2vXHshC90LzttsRNPZXnleA2Dnds+5O6CfF7Iyaa5IX7ytttqKIMRQTQed+SlzncX4AkpmkNDF3iuyxfsVg0FjC421vRvVOAzEluVFGcVst9m9DnxnyphA4FwysSRI6fEbP/E7SI/KOTjYJenmS37dK8/2rQw15mLocBma+lW/3V1RBrgBwvyyamNDVqXCl6le/N1NgPfHq12WL+7mfE/fYMVO/SEvM/Mg5DryuO+M+7jzYvexGFzMDW3gs1OB427Pku73MGmoaEevwgep6fX/SpyK3EoaBL93TQ31tNkM8hz5Pd6XG9M8FSw1xTyLkesf9DkCi1I69qr8Pm6CIrEzBGMH31c97IzykpPRCgJ2mniwOStHtaCoLJ4Cnaz55MTim8QTpii9f5V+2L1/gEzF2p2VqwgmJQ3lSabwZIN/0haZvhYp9OVdJ9oxo2Z2b3sMGKFx5GVuvfgq32PP47I5/bikXh3x+vswyMCWwbmo9FkuSLXKyKMCjpoljZCQzRZ6O1sp9WMe+S3+/qdslJ7fsarDsYUjqfZZqNpx4p+nb+7HFMQOM3e4bSiaUwpmEKW2eOrs9vxBA1OGTGZjz1O6qrWdGc1Cz8zQwyKxEWHo5O9ZnC69g5tQTQuaxRLsrJQXavSrluX6LpVO+wUNK2ivs3Lkq37mGss4aUVWmUaFgQ5zhyyHdmUm4HbphQfw3aHA3/N0M4TpULb3i0ERMjOLep1P5tho1y52e/Qznx1dVV0GUZ3nK9MmFYaCWLo6SEIRuePw6EUjf693ZkOo+cxon0H3I5YIVYsOdTb/RmHE+kPw1oQ2A07pQEtjbOJH0JOLj8RZ0ixpyNWNRROiu3oMXycPfk8AJZtejVpmd0JtR2991TC2OwOxphGQU57rPA47oizAfDZIz2y2s0fx1lxhFT/RgTPnf0YPxo7F5czduKrxDGCOoefLXsGwIS0pYaVT/+SeY89QVWSsAW+zjZajLAgCLAyQyeq7vOJwkG8IJhcPhuAuvqP+3X+MOGghmGrq5LsEp7/l+dZ+s2llDu1I2FxyM8plWfQaLPRsG8RwWCIQM/eorOYz50GdbVaBdTm1SOCK4+4SP92VKUVCRfoTm4UFKEoaz3vbdzH+F2vcL/zPo7Z8zwAXnT2PHsPvfmUspMIiuDrSH8kMth0NGpVT15e4rm3aMY6R1HthOo9tdQ06G99VE5ZH0cl56jxc7qX3c7cmG02w0ZJ0EYjzVGRiHvkITG/35CKfTfLssrY4zDYsHloLfVgmAsCgIk2/aJ4QvGTQ3abg8qgg32ie55vrtMWEt6AV2c+ssUec+IxF+AJhdjYsDSpY1k4W5gzRdUQQIXSL5O3R3rGopzRjAzC2pyI9cgUWcRHW2LnKboFQYYjgqllJ3DNnPhcBTNLZ7PLYWfFsuczOm80De/eydRN87ii+t+5488vJRxldHa2dJvp+m1+Vu/on1pKO/zEC4Jjy08BoN43MCqPSCiLWEEqIpxYrq3Fsp1ZzBj/ZQCau9biD/i1/jjqmU0pOILtDge1O3Vsq06/9n85etRxFCk7Te4m1u9Oz5qq0whQHtRGlp05u3l/fTVBMxd1lmkI0IUPT4KQC9NGaX13m317d6iKg4WmNh32pcI05+yNKYWTqbXbee2DBexv1qOg0oKxGZd91OjIKN7jyonbXmrks8/hZWedfn62Hu/gEQ7ttDd6ZGxn8YTxcwiIsGztsxnXLVOGvSA4qljr7FqTRPartI9ihzPIkk+3sPzpX3Lzwy/RHqwnLxT/Ydgdbk6zVbIsq5XHnv9lwvOFVUPuPnSX0RxZqJ1UgrZ4/eA4lc9WZ2Q043XvZ9OK92P26Z4s7uccQU/On3kVAOt3vd5vU8sNVYs4v7yciytGcWnXPN5YFx8WubktMvJolwClzaupbU0vVWiYQDCEVyTGezNMiaeEgpBBva2OXX3E0qlp7mT59t6d27qNCxKMAs83Q4+vdggV+WPJDxk0Ofayda8+Z/SI4JjS4wiIsGnnYgDag3pENDJvLNM9ZWzOUqzenLqDl1KKDiNEscpiRt5ElnqEI6v+SF17FeeUl9EW3MXe5i684scTim8KRmePpkw87Mjq4JOt8R74yQiGFLe9vJ4Ne2KFVqdv4EKWNPr1/M7E4qP73Hf6xDkA5NY8zI56bdZZMWpSxmVH2//nJhAEM4pnUu20s37DG0AkBWyYeV95mkvHnM60yefHrD99+qWIUmxtGJiRajoMe0EwZ9rXAHAkmRw6fdKFtNkM3l5wE6vHvsc5oVuosXdQGYoPEgdw26VPURAS3m16lpa2eKuToAqrhrLitiXjhxf8gQs8E7jm/N/FbZtRcnz3cpY9i/kj8vFU/T3mg1LmiCBZhrJMGTtyKuUhBxtdVby++JOMz+Nta+QvOS10GU4MVzYPjO5izWu/p7lH+Irmjkh4hf12O5fa/5FxjJ2Org58RrznJ+ie+tHZFWxxB3n/nVd6Pc9/PLWMqx96l68/uJhHP6hKqL+N9jvpyaySWeQ6c7n2uOsQEY51l7ItK8CKtXp+IlptMMkcqRjej2nu8NMR0h2DwuxiTqycQ43dTuPmJ1K8AzpqaZuhyDayuGT696lyOugo+gc7bdvY47BTnVfLL17dQJcEyUpgXQVwctE0VmY52b4suTq0J5/tbeGZJZu54N5/dr+nzZ1+Zv7879z77pY+jo7gC4RoSJIPoZF67AoqchN77Eczc/zZlNpzeaW4BXunvu9jTD+ZTJmVrcstKS2P23bOcd8EYFP9IiA2IRbA6Nwy/vuse2OijwLke4oYG3SwU3bT0JZZByhThr0gOLLyeG4/9ibmXfx0wu0XzP4+eUHF84VVbHY5+eNIJ9udBpXu+AcMkO0ewWWlF7LRbXD/n75Kly+2MQurhtzu1EcEHnc+v7nkZUYVTIzb9s0v3dy9fN8Z97HDYacmewUX3vsB7ab/Q1gQ2DIIQ90XN5xwI1tdDt5e9W3ufvAOnl20nMXb6tlY08K6Xc28sWIri7fW8/Dj9/DnF/7Gtro29rV00eYNEAiGWP75fh569g6WZrmZO/IM7j7z/9jtcLBl5AKuu/fPrNqpzSXrWr28uFLrRie5K6hy2Jnk/oTFKzLTT7e0a3Wf00isortk1o/Z47Czu/pXPPn+iu7wDaGQYmttK298vI7XV+9kcsvPOLPip8zt+CHeRT/i+rvm849NtTGqrXaf/mgTqQlsho3Fly/mW0d/C4Bzj7qY3Q47n2+Zr7dHjQgmFE6iwshmZX4DD734Fu2qnZygwmlzcu70q3Eo2Bt8nw83JTYjrW/z8j+vftodz6iusZlmm0GOLYfzJ5zP18tO49n8LD7M043rMo+gPn2FTgmRlSRVy5lHf5MOw2Bv3TPdYdCbOnz8fdXupB7HGzZv4abiG7kr+3+497kFhEKKpRs2cl7FT9iz9BYWb63nn+u2sWB5RDX3eX07L63aFdPBefCFt/n5r26NC/ZY39rFfnsHY0KuOP17IjwOD78+4w/sdth5vlDhDCnysgr7PK437pr7BDccfwOTy06K2za5bBYjg8IGtzYFtdvi5yeTcdqok1ifZfD0y3f3vfMAkplXRT8RkXOBewAbMF8p9ZvBLO/rx/9r0m0Ou5ML8k7gqfblAN1mYZOLpiU95nvn/IKPHl/CM7k1dD7yZb5+1r1MHj8Vj8tBQJmCII0RQW+MyivnOHcF7YF2ZpfO5vQRR/JCaANzdzzI3X+p5oTyLFr2boWs+J7HQHDOMVewqW4d83e8xiKepnLr44zeqDCUgYGizq4YEYQ1bjsun2LZs7/CHrJhD9kwQnaa3G187HFSHBCuP/u/yc3K5/bjbuCW1XdT6PwDj7/4IM+FcvAZbXyS7wUMriiby927nuBnJQH+Ze/N/M9dj4LhQsTAJjYMDGyGoZ2xxEBEUIgZa0lQIviD7ZAT7/kZ5vTx53Lq2j/yJJs4ZdM32Lg6Hwdj8NFMs2svex0Kl1KsKXEBHt7NUVC0h3L/Xbz11t28+8pE8jwTcLmL2N25BiNHUTSipM/7ee4xV/K7lffy0YjNgD1GEBhi8KPpP+Bnq+6itPknrM8yyDWt3QqzRnLWiKN4hw143riA5QvPoWLc6RSNqiQrdwQul5sXXn8Se9MLPLXRT2nRHLo63XTmGIwp1BOj1556O288cwbVDgfZ2NnsgoLKJ9nqdjPbn7ixOqXyNCrFw0f5+3A+eBNHn/ojHlq2ln3tb/HGikuYe8I0Tp9STK470rv9fPN8Hi/KATqZV/Vdnnj+9+xuXcTbOdmQ8xnbH3uUM0Y8CfZWnt78n3z5wn/lN3+aT27XYhav+Q6XzZnFjIoR7Nt7IysnwnOP7mbWebfxpWm6F79m5QdUueCo7MSdtUTMLJ3NmTlTeK9tEz5DEMks1lCYAncB3zn2Owm3iQiXlJ7FQ7XvAOCQ1AXB9868g2eeOY0lLX8l9OJxnPWFk5k8KhfD6F99+0KGOsywaP3FZuBsYBewDLhcKZV0qnzWrFlq+fLlg1anYDDALU9dzJHFM1jY/AkrOqt59ZwnGDc6uTDo8LVz3dNf5RO1F5tSlAZCeILCZtPRZN2/DlzUxpAKIeiXt6qpiitfvogWgjjNHpnPfEmeOPUhpk08ZcDKjWZbw2beWfcX1u9bS03XfnwhP0Gg1DmCpkA75e5i2gMd7PY10qmCdEqQThSjcZJjOPnPk25l9uQLus+3bPt7zF/ya5Z496EE7ArGG9lkY+eBi//OltYd3Pz+9dT6+hdq4uayy7jy7FsTbvMH/cz/+De8uPUV9kaFDi5WDsa5itnp20+r8vHX8/6K4cxmzZ6lvLfxGT5p246vx3d55chTuHnuQynVaeGWV7jho1sICfyk7GK+dfZt3duUUjyy5H95ePMz+AVmBD387eqlALT6WvnZa1exqHUzAA6lyAmFyAopnEqx02EnFNXA5QZDtNoM7vnSPM4wJ6rX7VvNg0t+yQ9P/Amr6tYyf/WDNBHg0qKT+e8LHklY38U73uf6hdfTJVDh99Nq2GiyGVT4/YzxKdwhOyFlJygGCFQ7vex02qjILqW6vYYTO7pY53bSmWDE6gwppnpD7LcrdjvsjAwEGeeD7ICTD/KC3dd5SocXT7AYu+TRSi0Lc/3c/4Vfc9qkuSnd8/D9u+bFiymyZfOHS15M+bhM8Aa6+NLjJ9IhilvKvsk3zv5pysc+v/RufvHZHxnr83Nsm4evzrmT2cedkVE9RGSFUmpWn/sdAEFwMnC7Uuoc8/dPAZRS8WYrJoMtCKIJqRDVLTsZmz+uz32VUry76kn+uelVajpr6FBdhJTiOq98MgAADtFJREFU6Nyp3HrZXwetjtubqnh9zaN0edvwhgSH3cmEkZO4aMb3+93TGUiUUn3WJxgK0uJr0Uk9ejjUBUIBdrbspDPYSSjoJxgKEAoFI/9VAFQIUFo9pvT/8LLD7mbm2DlxuthEtPpaqW7ZSbYju/vZK6XoDHTi6TEJ3BnoZFvDZvbt/4zWjnpKRkzgpAnnpnXvW32tbN6/iWNHTYuZfIzevrL6A8YXTqGy4IiYbdubqlhdtYBNe1bS2NFEV6CDzqCPcbmVXHDi93G58vnw07+zascHnDz+TK445cak9egMdLK2bi2TCyZ3J9VJxM6mz3l73V9YvnMJnYF2Ths9k4+bN1Hf2URTyEtIhTBQgCDAxSXn8t1z7+DPq+7nhQ1P0xDq4JLCLzJ7+td4fvl9NHqb+e4pt/L+mqdZXb+Kdgly+bjz2db4Geubq9lpdJGlhN99aR4f7nibBTsX0oC3O63q1KCHZ7/98UH1vvekw9vGKyvu4cLjf0x2jzS1fbHws+d5ePk9bAw08foFz3XnP0+Xg1kQXAycq5T6rvn7SmC2UuraHvt9H/g+QGVl5cwdO3bEncvCwmJ4opRCoWLCpgRDQToCHSgU2fbsjKLtHmq0+drIccbPPaVKqoLgoJ0sVko9opSapZSaVVycuReghYXFoYeIxMXOshk2cp255DnzDgshAPRLCKTDgRAEu4Fom69yc52FhYWFxQHgQAiCZcAkERkvIk7gG0DvxtwWFhYWFoPGkJuPKqUCInItsABtPvqYUmrDUNfDwsLCwkJzQPwIlFJvAG8ciLItLCwsLGI5aCeLLSwsLCyGBksQWFhYWBzmWILAwsLC4jBnyB3KMkFE6oAdQD6QKOZAovXhdUVAfR/7DsT6cDm91aU/ZRYB/gE4T1/7DsX96u1eDWSZPa+lt/0zLXMw7tdgPv++1juIv2cD/U5H37ND+X5VAjsH6dwDcb8q0W18345YSqlD5g94JNX14XXA8kzPkc76cDm91aWf17N8gM7T675Dcb96u1cDvH55qvtnWuZg3K/BfP59rU90zwb6nY4u4xC/X3WD+Sz6e78S1S/Z36GmGkoWFD3R+nT2Hez1Vl0OTF2Scbjfl+FWxwNVl0TJMg6m+5JyMo9DQjXUH0RkuUoh1sbBXs5wuY6hKmOoyhlO1zJU5VhlDE056Rx3qI0IMiFxbN1Dr5zhch1DVcZQlTOcrmWoyrHKGJpyUj5u2I8ILCwsLCx653AYEVhYWFhY9IIlCCwsLCwOc4aNIBCRtkE+f1BEVkf9jetl3zki8lqa51ci8njUb7uI1KV7nhTL+qpZXmZpj3o/95Bdh3n+QX3u6ZQlIotEJKPJw8F8Jj3KuUVENojIWvM9nj0IZZSLyMsiskVEtonIPWak4WT73yAinmTbkxyjRGRe1O+bROT2flQ7URnhb36DiKwRkRtFZNDazKF8l3sybATBENCplJoR9bd9gM/fDhwjIuFs62eTZp4GEUk1iODlwIfm/3TOn0o2kH5fx2FKRs8kHcw0sXOB45VS04CzgOoBLkOAF4G/K6UmAZOBHOB/eznsBiAtQQB4gYtEpCijiqZG+Js/Gv0enwfc1scxhyTDShCISI6IvCciK0VknYh8xVw/TkQ2isijpnR/O6qh6k95NhH5vYgsM3tYP4janCcir4vIJhF5KMWexBtAOMP75cBTUWWdKCJLRGSViCwWkSnm+qtE5BUReR94L4U65wCnAt9B54IIj2A+SFRfEWkTkXkisgY4OYVryPQ6PhCRGVH7fSgi01MprOcITETuF5GrzOXtIvKLqHeiXz3u3srqxzmTPZNk13S+iHwmIitE5N40RlulQL1SyguglKpXSu0RkZki8g/zfAtEpNQsZ5HZm18tIutF5MQUyjgD6FJK/cksIwj8O3C1iGSLyJ3mudaKyHUi8m9AGbBQRBameB0AAbRVzL/33GB+7++bZbwnIpUiki8iO6Le62wRqRaRvhNa6+uoRafOvVY0Sb99EflP811bIyK/SeOahrwNCzOsBAHQBXxNKXU8cDowz+yhAEwCHjClexPw9TTPnSURtdBL5rrvAM1KqROAE4Dvich4c9uJwHXAUcBE4KIUynga+IaIuIFpwNKobZ8BX1RKHQf8HPhV1LbjgYuVUqelUMZXgLeUUpuB/SIys4/6ZgNLlVLTlVIfpnD+TK/jj8BVACIyGXArpdakWF5f1JvvxIPATQN0zoEk2TOJw7ynDwPnKaVmAunkcX0bqBCRzSLyfyJymtkQ3od+f2YCjxHbe/copWYA15jb+uJoYEX0CqVUCzoUw3eBccAMc0TyhFLqXmAPcLpS6vQ0rgXgAeAKEcnvsf4+4C/hMoB7lVLNwGog/I3MBRYopfypFqaUqkLnUBlFkm9fRM5DP8/ZSqnpwO/SvKbBbMOSMtwEgQC/EpG1wLvAGKDE3Pa5Umq1ubwC/UKmQ7Rq6Gvmui8D3xKR1ejGbiT6YQF8opSqMntET6F7fL2ilFpr1uty4vM15APPich64A/oDy7MO0qphhSv43J0Q435P6yKSFbfIPBCiufuz3U8B8w1G6argT+nU2YfvGj+z+S5DwXJnkkipgJVSqnPzd9P9bJvDEqpNmAmumdbBzwD/AA4BnjHfI9vRaePDfOUeewH6FHuiFTLS8Ac4GGlVMA8Z6rvbEJMAfNX4N96bDoZeNJc/huRd/kZ4DJz+Rvm70xJ9u2fBfxJKdVh1jHdaxzMNiwpByQxzSByBbqHNFMp5ReR7YDb3OaN2i8IDMSwSoDrlFILYlaKzAF6Omik6rDxCnAn+qMZGbX+l8BCpdTXRE9UL4ra1p5SZUUK0UP3Y0VEoXs3Cni9l/p2mcIhXdK6DqVUh4i8g+5NXYpusFIlQGynxt1je/jZB+n/O99XWWnRyzN5eSDLCWM+y0XAIhFZB/wY2KCUSqb2S/c9/hS4OHqFiOShA6BtT7e+KXA3sBL4Uwr7voJuZAvR79f76RQkIhPQ71Atyb/9c9I5ZwKGug0Dht+IIB+oNW/g6cDYQS5vAfCjsJ5RRCaLSLa57URzqGigeyGpqlUeA36hlFrXY30+kUnXqzKs78XA35RSY5VS45RSFcDnwBf7Ud9kZHId84F7gWVKqcY0ytoBHCUiLrPHemYG9T1QZSV7JkaScjYBEyRitXZZzxMmQ0SmiMikqFUzgI1AseiJZETEISLRo83LzPWnolUhiSJfRvMe4BGRb5nH2YB56BHeAuAHYho1mA0yQCuQm+p1RGP2uJ9Fq2rCLMaca0E3rP80921D50y/B3gtnQ6OiBQDDwH3K+2Fm+zbfwf4tphWUFHXmCpD3YYBw0QQmC+WF60PnGX2dL6F1kcPJvPRPaCVpqrjYSI9zmXA/egP7XPgpYRn6IFSapepN+3J74Bfi8gqMu/VXp6gHi+Y6zOqbzIyuQ6l1AqghdR6d93PXSlVjW4M1pv/V/Wn7kNcVrJn8o1E5SilOtH6+rdEZAW6Ee2rcQ6TA/xFRD41VQ9HoedpLgZ+K9ogYDXwhahjusxn9RCxjW1CzEbya8AlIrIF2IzWe/8M/b3sBNaaZX3TPOwR83rSmSyOZh46VHOY69CN8VrgSuD6qG3PAP+P1NRC4XnBDWg1zdvAL8xtCb99pdRb6JHHclNtlNKc1AFsw3T5wyHEhGjrkkeVUqlYNVj0wFRl3aSUmnuA61GGVltMVUqFUth/yJ77wfSOiUiOUqrNnER8ANiilPrDIJSzCP1eLB/oc1vEcqDfr0N+RCAiP0RPaN16oOtikTmmKmEpcEuKQmDInvtB+I59z+xtbkCrEh4+wPWx6AcHw/s1LEYEFhYWFhaZc8iNCETkMRGpNfVy4XXTRTsprRORV00rBUTEKSJ/kohzxxxzvUe089Rnop0z0nL6sLCwsOgPA9GOmdveMtdtEO0Imor3fxyHnCBAWx+c22PdfOC/lFLHoifefmKu/x6Auf5stHNG+JrvVEpNBY4DThHtCGJhYWExFPyZgWnHLjUd145Bm51ekkllDjlBYDq29HTSmAx8YC6/Q8Tj7ihMW2HTRbwJmKWU6lBKLTTX+9B2yOVYWFhYDAED0Y6Zv1vMfeyAk9T9lWI45ARBEjagHZFAS8QKc3kN8C+iI2CORzuRVEQfaNpoX0gKcXosLCwsBpGM2jERWYB2cmsFns+k4OEiCK4GrjHtqnMBn7n+MWAXsBztgbgY7ZEHdNvuPoWORVI1pDW2sLCwiCWjdkwpdQ46oKAL7aWeNsMixIRS6jN07I9wwLILzPUBoqITishitINLmEfQNth3D11tLSwsLOLpRzuGUqpLRF5GjyjeSbfsYTEiEJFR5n8DbYv7kPnbEw75ICJnAwGl1Kfm7zvQNtg3HJBKW1hYWESRbjsmOmR1OGS4HS04MvJEPuRGBCLyFDqQWZGI7EInisgRkR+bu7xIJETBKGCBiITQ8W2uNM9RDtyCvmkrtYMm9yul5g/VdVhYWBy+DEQ7hg4R/4qIuNCd+oWYwiPt+lgOZRYWFhaHN8NCNWRhYWFhkTmWILCwsLA4zLEEgYWFhcVhjiUILCwsLA5zLEFgYWFhcZhjCQILix6IyAgRucZcLhORjNz2LSwOFSzzUQuLHojOB/yaUuqYA1wVC4sh4ZBzKLOwGAJ+A0w0s4BtAY5USh0jIlcBX0U78kwC7kRHfLwSnW/2fKVUg4hMRKeQLAY6gO+Z4QMsLA5KLNWQhUU8/wVsU0rNIBITPswxwEXACcD/Ah1KqeOAJehk46BjWF2nlJqJTl7+f0NSawuLDLFGBBYW6bFQKdUKtIpIM/CquX4dME1EcoAvAM+ZoUtAR4W0sDhosQSBhUV6eKOWQ1G/Q+jvyQCazNGEhcUhgaUasrCIpxUdDz5tzIxRn4vIJQCimT6QlbOwGGgsQWBh0QOl1H7gIzOx+O8zOMUVwHdEZA2xWacsLA5KLPNRCwsLi8Mca0RgYWFhcZhjCQILCwuLwxxLEFhYWFgc5liCwMLCwuIwxxIEFhYWFoc5liCwsLCwOMyxBIGFhYXFYY4lCCwsLCwOc/4/xlucxv5IOpgAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import os\n", "import pandas as pd\n", "from matplotlib import pyplot as plt\n", "\n", "import swimpy\n", "\n", "# load the project instance\n", "p = swimpy.Project()\n", "\n", "# run for two years and without subcatch\n", "p.config_parameters['nbyr'] = 2\n", "p.basin_parameters['subcatch'] = 0\n", "\n", "# run the model while changing the basin-wide ecal parameter\n", "ecal_values = [0.9, 1.0, 1.1]\n", "q = pd.DataFrame()\n", "for ecal in ecal_values:\n", " p.basin_parameters['ecal'] = ecal\n", " p.run(quiet=True, save=False)\n", " # keep daily discharge at Blankenstein station\n", " q[ecal] = p.station_daily_discharge['BLANKENSTEIN']\n", " print('ecal=%s' % ecal)\n", "\n", "# visualise results\n", "ax = q['1992':].plot()\n", "plt.title('Sensitivity of ecal at Blankenstein')\n", "yl = plt.ylabel('Discharge [m^2s^-1]')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Convert to a (general) project function\n", "If you are doing this sort of sensitivity testing more often, you can convert this little script to a project function. It is then accessible on the commandline and as method on the *project* instance. Put the below code into your `swimpy/settings.py` file." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "def sensitivity_basin_parameter(project, parameter, values, plot=True):\n", " \"\"\"Run the model with the parameter set to each value.\n", " \"\"\"\n", " q = pd.DataFrame()\n", " for val in values:\n", " project.basin_parameters[parameter] = val\n", " project.run(quiet=True, save=False)\n", " q[val] = project.station_daily_discharge['BLANKENSTEIN']\n", " print('%s=%s' % (parameter, val))\n", " if plot:\n", " # visualise results\n", " ax = q.plot()\n", " plt.title('Sensitivity of %s at Blankenstein' % parameter)\n", " yl = plt.ylabel('Discharge [m^2s^-1]')\n", " return q" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "nbsphinx": "hidden" }, "outputs": [], "source": [ "# hidden: like adding it to settings.py \n", "p.settings(sensitivity_basin_parameter)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The function `sensitivity_basin_parameter` is a generalisation of the above script and alows testing the sensitivity of any basin parameter (in the `.bsn` file). You can now use it in a python script, e.g.:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "ecal=0.9\n", "ecal=1.0\n", "ecal=1.1\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEiCAYAAAD05tVnAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsvXmYZGV5sH8/59TW26wMyywsigruGhCNy6diYsQticZoiLgb/fziihg1KmqiGI3BuENcUBH1pxggAZQgi4hAAEVgQIdlgBmYYfbeajnL8/vjvKfqVHV1d1XX1qf7va+rr66qs9RbdU69z/vsoqpYLBaLZfniDHoAFovFYhksVhBYLBbLMscKAovFYlnmWEFgsVgsyxwrCCwWi2WZYwWBxWKxLHOsILDMiYhMisgj5th+u4g8t4XznCwiP+/q4GZ/r2eKyBYz9j/v03ueLiLf68d7zfL+W0XkBV04z+tF5JpujKlXtHrPWVrHCoIUISLPEpFrReSAiOwVkV+JyPG9fE9VHVXVe8z7f1tE/qlh++NU9coWznOuqv5p/FxEVESO7vqAIz4BfMmM/T979B5dQ0SuFJE3z7H9SPN9TZq/nSLyFRHJ9nOcvaLde6HVe87SOlYQpAQRWQH8F/BFYA2wAfg4UB7kuBYpRwC3D3oQPWCVqo4CTwCeAbxjwOOxLBGsIEgPjwZQ1fNUNVDVoqr+XFV/F+8gIm8UkTtEZJ+I/ExEjkhsUxF5mzGZ7BeRL4uImG1Hi8hVRtPYLSI/bDjuaBF5K3AycJpZlV5ktm8VkReIyHoRKYrImsSxTzHnyyZNDiJytdnlFnOuvxaR20TkpYljs+bYpzT7MkTkLSJyl9GMLhSR9eb1u4FHABeZc+ebHLteRH4iIrtE5F4ReWdimysiHxKRu0VkQkRuEpFNZtsXROQBERk3rz+7lQsnIqtF5L/M++0zjzeabf8MPBv4khnvl+Y7n6o+DFwGPHaW93uaiPzaXOeHRORLIpJLbJ/1Xmhyrs+KyDUistI87+o91uxeMK+/RER+a85zrYg8MfE+VTOYRCa5H4nId8z1ul1EjpvvO7Q0oKr2LwV/wApgD3AO8CJgdcP2lwN3AccCGeAfgWsT25VIo1gFHA7sAv7MbDsP+DDRwqAAPKvhuKPN428D/9TwvluBF5jHvwDektj2WeBr5vHrgWuandc8Pw34YcPnuXWW7+L5wG7gqUCeSEu6utmYmhzrADcBHwVyRELjHuCFZvv7gVuBxwACPAlYa7b9LbDWfL/vA3YABbPtdOB7s7znWuAVwDAwBvx/wH8mtl8JvHmOa3+k+b4y5vl64BbgjbNchz8Cnm7GeSRwB/DuFu+F1wPXmO/pbOBnwHC/7jHz/CnAw8AJgAu8zny+fJPPejpQAk4y+34auG7Qv9e0/VmNICWo6jjwLKIfzdnALrMSPsTs8jbg06p6h6r6wKeAJydXbMAZqrpfVe8HrgCebF73iMwp61W1pKoLdRZ+H3gNgFkJvtq81grfA06SyAQG8Frgu7PsezLwTVW9WVXLwAeBZ4jIkS28z/HAOlX9hKpWNPJ/nG3GCvBm4B9V9fcacYuq7gFQ1e+p6h5V9VX1X4mE0GPme0NzzE9UdVpVJ4B/Bv5PC2NtZLeI7Ae2A1PAj2d5v5tU9Tozzq3A15u832z3AkCWaOJeA7xUVafN6/26x94KfF1Vr9dI+z2HyAT69Fn2v0ZVL1bVgOieedIc57Y0wQqCFGF+gK9X1Y3A44lWhmeazUcAXzCq9H5gL9GKdkPiFDsSj6eBUfP4NLPvDUa1fuMCh/gTogn5MOA5QAj8spUDVfVB4FfAK0RkFZHWc+4su68H7kscO0mkLW2YZf8kRwDr4+/JfFcfAmKBugm4u9mBInKqMYscMMetBA6a7w1FZFhEvi4i94nIOHA1sEpE3BbGm+QgVV1FpFn8imi13uz9Hm3MTzvM+32qyThnuxcAjiZa/X9cVSuJ1/t1jx0BvK/hGm0iuu7NaHzPgohk5ji/pQErCFKKqt5JZKp5vHnpAeDvVHVV4m9IVa9t4Vw7VPUtqroe+DvgK9I8imPOUrWqug/4OfDXwN8AP1DVdsrbnkNkfvkr4Nequn2W/R4kmiwAEJERIvPLbPsneQC4t+F7GlPVkxLbH9l4kPEHnAa8isgstwo4QDS5zcf7iDSHE1R1BZGQJHFsWyWAVbVIdO2fLiLNBNFXgTuBR5n3+1CL44y5A3gDcImIJDWeftxj8fv8c8P7DKvqeW18BksbWEGQEkTkGBF5X8LJuInIDHOd2eVrwAdF5HFm+0oR+asWz/1X8XmBfUQTU9hk151ENvW5+D5wCvBK5jYLNTvXfxLZ/d8FfGeOY88D3iAiT5bIGfwp4HpjBpmPG4AJEfmAiAxJ5Bx+vNTCcP8D+KSIPEoinigia4ls+z6R3TsjIh8l8tu0whhQBPZL5Ez/WMP2Vr7XKuYzv5ZoJbxnlvcbByZF5Bjg7a2eO8ZMuh8C/kdEYsHYq3us8fOfDbxNRE4w12BERF4sImPtfg5La1hBkB4miJxn14vIFJEAuI1otYmq/hT4DPADYw64jci80grHm/NOAhcC7zK280a+ATzWqOuzxedfCDwK2KGqt8zxnqcD55hzvcp8hiKReeko4PzZDlTV/wE+YvZ9iGgF/+rZ9m84NgBeQmS7vpfI6fwfRGYegM8DPyLSbMaJPvMQkRnmUuAPRGapEtHKtRXONOfYTXTdLm3Y/gXglSYS59/nOM9+c412EoWPvmwWjetUIo1sgmhS/WGTfebF2OY/AfxCRI7s4T12Ool7QVVvBN4CfIlIaNxF5MS29AhpT3O3WHqLWWk/WlX/dtBjsViWC9ahYlk0GLPJm4jMHhaLpU9Y05BlUSAibyEytVyiqlfPt7/FYuke1jRksVgsyxyrEVgsFssyxwoCi8ViWeakwll80EEH6ZFHHjnoYVgsFkuquOmmm3ar6rr59kuFIDjyyCO58cYbBz0Mi8ViSRUict/8e1nTkMVisSx7rCCwWCyWZY4VBBaLxbLMSYWPwGKxWLqN53ls27aNUqk06KF0TKFQYOPGjWSzC2tjbQWBxWJZlmzbto2xsTGOPPJIpHmnzlSgquzZs4dt27Zx1FFHLegc1jRksViWJaVSibVr16ZaCACICGvXru1Is7GCwGJJGapKGNrSMN0g7UIgptPPYQWBxZIyPn7RZh7xoYsHPQxLF7j00kt5zGMew9FHH80ZZ5wxY/t9993HiSeeyBOf+ESe+9znsm3btp6MwwoCiyVlfPvarUCkGVjSSxAEvOMd7+CSSy5h8+bNnHfeeWzevLlun1NPPZVTTjmF3/3ud3z0ox/lgx/8YE/GYgWBxZJSAmseSjU33HADRx99NI94xCPI5XK8+tWv5oILLqjbZ/PmzTz/+c8H4HnPe96M7d3CRg1ZLCnFyoHu8fGLbmfzg+NdPedj16/gYy993Kzbt2/fzqZNm6rPN27cyPXXX1+3z5Oe9CTOP/983vWud/HTn/6UiYkJ9uzZw9q1a7s6VqsRWCwpI/YLhtY0tOT53Oc+x1VXXcVTnvIUrrrqKjZs2IDrul1/H6sRWCwpQwDFmoa6yVwr916xYcMGHnjggerzbdu2sWHDhrp91q9fz/nnnw/A5OQkP/nJT1i1alXXx2I1AoslZcShglYjSDfHH388W7Zs4d5776VSqfCDH/yAl73sZXX77N69mzAMAfj0pz/NG9/4xp6MxQoCiyVlOLFpKBzsOCydkclk+NKXvsQLX/hCjj32WF71qlfxuMc9jo9+9KNceOGFAFx55ZU85jGP4dGPfjQ7d+7kwx/+cE/Gkoqexccdd5zafgTdY/90hRWFLI6zNJJplhvv+8cP8xbnAg7+wM2sGckNejip5Y477uDYY48d9DC6RrPPIyI3qepx8x3bU41ARN4jIreLyG0icp6IFETkKBG5XkTuEpEfioi9k/vI3qkKT/7EZfzrZb8f9FAsC+RfM1/iGOeBqsnAYumUngkCEdkAvBM4TlUfD7jAq4HPAP+mqkcD+4A39WoMlpnsn64AcPGtOwY8EkunhIE/6CFYlgi99hFkgCERyQDDwEPA84Efm+3nAH/e4zFYErjGHOTb1WTqCcJg0EOwLBF6JghUdTvwOeB+IgFwALgJ2K+q8VJmG7Ch2fEi8lYRuVFEbty1a1evhrnsiAVBECx+35ClOeePjPJnG9cTWkFg6RK9NA2tBl4OHAWsB0aAP2v1eFU9S1WPU9Xj1q1b16NRLj8yTnTJfRuDnlo+dvAatmczhJ41DVm6Qy9NQy8A7lXVXarqAecDzwRWGVMRwEZgew/HYGnAyAEbg74E8PzKoIdgWSL0UhDcDzxdRIYlyoA5EdgMXAG80uzzOqA3VZQsc2KzUtNPEFpBkHbe+MY3cvDBB/P4xz++6XZV5Z3vfCdHH300T3ziE7n55pt7Mo5e+giuJ3IK3wzcat7rLOADwHtF5C5gLfCNXo3B0gQz/1vTUHoRo835Nmoo9bz+9a/n0ksvnXX7JZdcwpYtW9iyZQtnnXUWb3/723syjp7WGlLVjwEfa3j5HuBpvXxfy/zYDlfpJa415AfeoIdi6ZDnPOc5bN26ddbtF1xwAaeccgoiwtOf/nT279/PQw89xGGHHdbVcdiic8uMePq3GkF6ifPBg8CahrrGJf8AO27t7jkPfQK8aGbXsXZoVqp6+/btXRcEttbQMsU6i9NLTRBY05ClO1iNYJmh1keQemJB4FuNoHt0uHLvFa2Uqu4GViNYZmhQ4ZzsGTyZLYMeimWBiJHhgc0OX/K87GUv4zvf+Q6qynXXXcfKlSu7bhYCqxEsO9wD9/F/3N+xSR4G3j3o4VgWgNUIlg6vec1ruPLKK9m9ezcbN27k4x//OJ4XBQG87W1v46STTuLiiy/m6KOPZnh4mG9961s9GYcVBMsMdbKEgIstT5BWYkFgi86ln/POO2/O7SLCl7/85Z6Pw5qGlhmTfoknHXU4/7nKXvq0UnUWhzZ81NId7GywzNgbTAFwyZhVBtOKFQSWbmMFwTJDzCUPbXOy1OLEzuLAmvcs3cEKgmWGI0YQDHgcls4JQusjsHQHKwiWGY7VCFKPNQ1Zuo0VBMsMqxGkHzGFQkKrEVi6hBUEy4w4GclqBOlFjE5gw0fTz3xlqO+8806e8YxnkM/n+dznPtezcVhBsOyIdAEfKwnSSvyj9dU6i9POfGWo16xZw7//+79z6qmn9nQcVhAsM9RqBEsGaxpKP895znNYs2bNrNsPPvhgjj/+eLLZbE/HYYPJlxlqNALrI0gvVfOeWkHQLT5zw2e4c++dXT3nMWuO4QNP+0BXz9krrEawzFA1piGxKkFasWWoLd3GagTLjKppaLDDsHRAtdaQ1Qi6RlpW7r3CagTLjkgSBFYjSD1haJ3Flu5gNYJlhu1Mln7i1Vtoo4ZSz3xlqHfs2MFxxx3H+Pg4juNw5plnsnnzZlasWNHVcVhBsMywesDSwWoE6We+MtSHHnoo27Zt6/k4rGlomRHarlapJy46Z8NHLd3CCoJlhmJNQ0sFaxqydAsrCJYdVhCknVrUkBUElu5gBcEywzqL00+11pD1EXSMLpHfQ6efwwqCZcZSufGXM1Yj6A6FQoE9e/ak/jehquzZs4dCobDgc9iooWWHdRanHSsIusPGjRvZtm0bu3btGvRQOqZQKLBx48YFH28FwbIj3asfS00QqBUEHZHNZjnqqKMGPYxFgTUNLTcSciAMrVBIM6Fa7c7SHWbVCETkwhaO36uqr+/ecCy9JkxIAj9Uco5NMUsbogKoNQ1ZusZcpqFjgTfPsV2AL3d3OJZeo4lVZGA1glRifQSWbjOXIPiwql4118Ei8vEuj8fSaxJzvx+GgDuwoVgWRk0QWNOQpTvM6iNQ1R/Nd3Ar+1gWF0kdwGoE6cQ6iy3dZkHOYhE5q9sDsfSH5OThW0GQakIbCmzpEnM5i2drpCnASb0ZjqXnJOZ+qxGkkyizWK1GYOkac/kIdgH3UV+5WM3zg3s5KEvvSGZRWo0g3VgfgaVbzCUI7gFOVNX7GzeIyAO9G5KllySrjwa+nUjSSNy8Xq1pyNIl5vIRnAmsnmXbv7RychFZJSI/FpE7ReQOEXmGiKwRkctEZIv5P9t7WHpAUiPwAm+AI7EslJqz2AoCS3eYK2roy6p6yyzbvtji+b8AXKqqxwBPAu4A/gG4XFUfBVxunlv6RHIVGfpWEKQRGz5q6TZtRQ21Ey0kIiuB5wDfAFDViqruB14OnGN2Owf483bGYOmQpI/Atx2u0khVEFjTkKVLtBs+elwb+x5F5HD+loj8RkT+Q0RGgENU9SGzzw7gkDbHYOmAZMXdMLCCII1UfQRWI7B0iXYFwcNt7JsBngp8VVWfAkzRYAbSyGDdNHRFRN4qIjeKyI1LoUzsYqHeWWxNQ2nGOost3aItQaCqf9bG7tuAbap6vXn+YyLBsFNEDgMw/5sKF1U9S1WPU9Xj1q1b184wLS0SWmdxKhEjzK1GYOkW8woCEXnCQk6sqjuAB0TkMealE4HNwIXA68xrrwMuWMj5LQujLo/As4IgjahxElgfgaVbzNmYRkROBE4Hnr3A8/89cK6I5IjyEt5AJHx+JCJvIkpYe9UCz21ZCAlBEIRWEKQRrf63CYGW7jBXiYmTgfcBL1zoyVX1tzR3MJ+40HNaOiNpVw788gBHYlkosQCw4aOWbjGXRvAt4BhVtZ7aJUSdszioDHAklk6xzmJLt5jLR/Ax4BsiMtSvwVh6j9aZhqxGkEasacjSbebKLP408E3gP/s3HEs/CUKbR5BG4unfOost3WJOZ7GqfldEHuzXYCy9JxlyGIbWNJRGqhqBWo3A0h3mDR9V1cv7MRBLf7A+gvQTX0MVKwgs3aGVPIK/EpEx8/gfReR8EXlq74dm6QnJxjRWEKQTm0dg6TKtZBZ/RFUnRORZwAuIish9tbfDsvSKOo3A5hGkEmsasnSbVgRB3A/vxcBZqvrfQK53Q7L0kqQgCK1GkEpqzmIrCCzdoRVBsF1Evg78NXCxiORbPM6yGEmGj9qEpFRS1Qisj8DSJVqZ0F8F/Ax4oeknsAZ4f09HZekLgdrw0TRiTUOWbjNn+CiAqk4D5yeePwQ8NPsRlsVMMnzUt3kEKUWJapBaQWDpDrNqBCLyBBG5TkQeEJGzkr2FReSG/gzP0nUSq8jQCoJUYn0Elm4zl2noq0SVR58A/AG4RkQeabZlezwuS4+ocxZb01AqictQWx+BpVvMZRoaU9VLzePPichNwKUi8lpm6SpmWfwkBYE1DaUTW2vI0m3m60ewUlUPAKjqFSLyCuAnRA5jSxpJ9iy2E0kqsaYhS7eZyzT0GeDY5Auq+juiXgLnNz3CsuhJTh6qwRx7WhYr2vDfYumUWTUCVf3+LK/fD7ylZyOy9JaEs1hDKwjSiNUILN2mZz2LLYuVpLPYJpSlkaqz2AoCS5eYUxCYnsVf6dNYUoOq8p+/2U7FT99EGqoVBEuFWCBYLJ0yVx7BycBngb/s33DSwe93TvDuH/6WX25JYxfPhGnIVq9si6v/sIuLbhl8ew5rGrJ0G9uzeAF4fvQDLHkpnEiTUUPWWdwWp3wzyqN86ZPWD3QcNUFgsXQH27N4AcS2WT9M308xqQVoCsdvqQkAm1Bm6Ra2Z/FC8Ct8JPNdpLh/0CNpm+TUYRubtM7DE6Xq48ny4kjEs1fP0i1sz+IFsPKei3hT5hLu3jwMz/jmoIfTFsmKlTaPoHVO/Nerqo93jpcYXTc6sLHYPAJLt7E9ixeCmUwlKA94IO0jNmpoQUyUalrAvqnBNvTRaqtKKwos3WHeMtQi4hJ1Jzsyub+qfr53w1rchI5rHqRvRZ2cOtQKggUx6LDhqrPYho9ausS8ggC4CCgBt2LNkgBorEilUhAknMVWECyIcjDY763qLB7oKCxLiVYEwUZVfWLPR5IiVCJBkEYbu9rw0Y4ZtEYQYwWBpVu00qryEhH5056PJEXsD8u8aONhPEBx0ENpG7EJZQvmkbKdx8s9AxcE1jRk6TataATXAT8VEQfwAAFUVVf0dGSLmP8t3s22bJafhvt47qAH0yZJB7Htedsel+ejVt0/8f9iwCOJsLWGLN2iFUHweeAZwK1qZw4AIv95+qM2rEawMHxvsNFicdRQgFUJLN2hFdPQA8BtVgjUcIz8TGNCVr2PIH3jXwwUxu8d9BAAm1ls6R6taAT3AFeKyCVAdSm0nMNHSbFGUKcFWNm+IJzSYDPKa7WGrEZg6Q6tCIJ7zV/O/C17HCfWCNI9kaZRo1kMhF5p/p16iO1ZbOk28woCVf14PwaSJtLsI0jmDljTUHv46pCRkNAfsI/A/A/FagSW7jBXP4LT5zu4lX2WIo5RycOUm1ass7g1VJUjZAcZMd/XYhEEAx2FZSkxl0bwZhEZn2O7AK8GTu/qiFKAY+RnOjWCZNE5O5W0gipclns/lw8Psdt1yQ44aoh4IUJ0PcVqBpYOmUsQnA2MzXP82V0cS2qI1ajUC4IUjn8QKJCTgHcfsg6A9/v7Bj4eiMJIQwXXygFLh8wqCLrlGzBF624EtqvqS0TkKOAHwFrgJuC1qjrYco5t4phA7jDl4XtWI2iNWHiKKirCnmDvYMdTzSOAIFRcx0oCS2e0kkfQKe8C7kg8/wzwb6p6NLAPeFMfxtBVYlU8jRpBskKNjRpqjdB8ZVnz3w8HZxpS1UTUEIS2y5ylC/RUEIjIRqIS1v9hngvwfODHZpdzgD/v5Rh6QVUjGPA4FkKdOSjlzu5+EX9nrvkf6iAFQe1xIGLLhFi6Qq81gjOB06jNmWuB/aoad/nYBmzo8Ri6TjwxpFEjqPcRpFGU9Z/4KzNdKAgHaMlU6juUBcHiaJtpSTfzCgIRebSIXC4it5nnTxSRf2zhuJcAD6vqTQsZmIi8VURuFJEbd+3atZBT9IzYtp5GQZCsPpr28Nd+45rvK2CAgiCs3XUBEKawJ4Zl8dGKRnA28EGiyqOo6u+Iwkbn45nAy0RkK5Fz+PnAF4BVIhI7qTcC25sdrKpnqepxqnrcunXrWni7/hFrBGmcRus6lKXc2d0vYoFZ0wi8gY1FE8Y9FSH0I0HgByFhaK+nZWG0IgiGVfWGhtfm1UdV9YOqulFVjyQSHL9Q1ZOBK4BXmt1eB1zQxngXBxqbhtJHnU3ZRg21hGp0wx9wop9LMEhBEEaRSzF+GI3l6A9fwtu+tyDl22JpSRDsFpFHYhaTIvJK4KEO3vMDwHtF5C4in8E3OjjXQIht60EqV9TJqKE0jr//qCqv2nAogZmAAwapEUT3XtVMlfAR/HzzzoGMyZJ+Wik69w7gLOAYEdlOVIDub9t5E1W9ErjSPL4HeFpbo1xkaJo1guRjqxG0RKghW3K1eovh/Apx78Zi+iU7RD4CPxicULLMTdkPKPshKwrZQQ9lXlopOncP8AIRGQEcVZ3o/bAWN7EgSON62mYWt482xOp7DM5BW9MIwBMIQisIFit/+x/X879b97H1jBcPeijzMq8gEJH3NjwHOADcpKq/7dG4FjVVjSCVCZ1WELRNgxM2GKAuGCeQxTbdwI+0k49mvsP14TFEaTuWxcD/bh1sKZJ2aMVHcBzwNqJ4/w3A3wF/BpwtIqf1cGyLljDF4aNWELRPY4imL4PTCOJ7z61mOUcawRszl/L13JmDGpZlDsr+4g/xbUUQbASeqqrvU9X3AX8EHAw8B3h9D8e2aKkllKWP0JqG2qbxe/IHeOVjbbSqEQSLf5JZ7hwoLn7zXSuC4GASLSqJ8gkOUdViw+vLhnSbhmrY8gStoTM0gkEKAmMaatAILIuPOMp3PAWCoJWooXOB60Ukjvd/KfB94zze3LORLWriDNM0SoKERjDACS1NNGZgD9IkGJuGqqXQjY/gmqECh3s+hw9oXJaZjOYzTJT8VGgErUQNfdI0rn+meeltqnqjeXxyz0a2iKmWmEihHKjLJ7MaQWs0hNkOsmpr7K9wzSLEN2W73n7owQDcOphhWZowtlQEgeklcLuqHkPUU8BCzWacynJfick/nc7u/qNBgyAYYCJhfO/FpqEw8G1piUVKLhPpbSVv8Wvec/oIVDUAfi8iVuNMEmd1prBFYL3j004grRA0+AiChu9t6+4p/u2yP/RFw0omlEHkI7DFAxcnccMgL1j8gqAVH8Fq4HYRuQGYil9U1Zf1bFSLnDChEaSvZ6zVCNol1HpB0KgRnPLNG7h/7zQnP/1wDh4r9HQstaghAZQwDBrTHCyLhMeGW3iScxde8KRBD2VeWhEEH+n5KFKHEQQiqe4Za8NHWyMI642AjQJ092T/gudi/0RsGgoCqxEsVr44dSrk4LzgHYMeyry04iy+qh8DSROaMA15QYjruPMcsXiweQTt01iTqbHY4NP8m/iH3HlUyr+CHmsEVDOLI40gCIO6dpXp01CXPv5SSCgTkaeLyP+KyKSIVEQkEJHxfgxusZK0Bfsp1sutaag1GpO2Gn/Wn81+nWOcBwgnd/d8LNXMYhM1FAmC2oiCFN+PS5XAKw56CPPSSkLZl4DXAFuAIeDNwJd7OajFTnKFmAZpX09idWtNCi3R2NKz0UeQM2WpK31wCsarfzGCIAz9OkGQ5oXJUsUpL/51c0s9i1X1LsBV1UBVv0VUayh1/PruPZxxyZ1dOFNCI/AH17ZwIdTHDCk//c02Sl7ahFl/aewL3PhtZc0rfqXU87HEPoLYGHnF7x+qC29NQ4TKcsOpLP6Cza0IgmkRyQG/FZF/EZH3tHjcouP1Z1/NBVc1Nltrn6RtPaikSxAkRYEXBLznh7d0STguXRrLUAdCXex+rBH0QxBUx6KRRnDzfXvqBFX6NNSlS0nyALiVpaERvNbs9/+Iwkc3Aa/o5aB6xXm5j3Nt4e87TsCpMw15vf/xd5VkQpkxcWx5eIKb709Pydx+E9AYPlpvgjl/xTBPOOpwiuXer/xiH0HOtM08ojBZV2/I85Zl+a9FwUMHijw8XpsPdrlDbMu4uEtBI1DV+1S1pKrjqvpxVX2vMRWljtM3lXjZhsO6assNvXRpBPH05Wqc5FGqAAAgAElEQVRNr/nVXXv4y69cax2NsxCaFfdrCifwrPBgfOqdsl9ZtRKA8fL+no8l9ldkjCB4Z/jdui5lvrf4yxksVZ7x6V/wtE9dXn3+pkNHedGmDTje1BxHLQ5aaUzzTOB04Ijk/qr6iN4Nqzfcm4taxpUrHoXswkM+kyGYgZ+uFVgc8eQqaIPTM22hsP0ivt6uOLjiEgj4YUhsqfdMtKZX6UN0iDENZQMPXGOmSkQ1+f0Yg6UlHsqadXaw+BeLrSSUfQN4D3ATM/1kqaFUqU3YpdIUK0c6iPdOCIJySn94Ls0FQScCcqmiJrPYwSUjLgH1GoFn4vbL3nTvx2L+V/IHge7g3sxRrEiYhkopvR+XMl7Q+/uiU1rxERxQ1UtU9WFV3RP/9XxkXWayWHPYFEuTHZ0rGU5YSZtNNqkRmGllDeM807kVL7CmoWbUMosFBxdfwE+YFyuxIPB7bwKoJo85kXZ7R87B92uCoJwCM8RyoxQs/msyq0YgIk81D68Qkc8C55NoRKOqN/d4bF1lYromCKbLnQqCGhU/Xc7iePJ3qZXRPjf3KY517men9/+A3MDGtlip9gl2ItOQL2Js8ZFWGRcfLPVDI6hqJ9F7nrtqir9M5BGUy+m6H5cijdnd5TDFggD414bnxyUeK/D87g+nd0wXDyQed+jFT5qG0ioIEhrBsc79AHiVMjA6qKEtWuIoMREh42QIAvCb2H29oPdmmdhf8Sj3EG4O7udxJbcufNRqBINB6+aEyMQ6HMK0AyWtLRDece7N/PetD7H1jBcPYpizMqsgUNXn9XMgvWYqkd1XrnT2Y6k3DaVLEJDQCBoNQV7KIqD6RVhdhdc0gtCfGZ3TD+2wGj4qGR7pZclpQJCojupVigShcu3du3naUWvIZ6zPpx8kzarjpSgYpVAVBLUFwn/f+tAghjcvrdQa+pSIrEo8Xy0i/9TbYXWfpF+gVO5MECSjhrwgXT6CeOSOyoziaelLjusPVbu8OLiSIRDBbyI0K0EfFgXm3hMEFyEgrPMRVPwin774Dl77jRu48ve7ej+eLnPXwxN9rebaLYqJ7PzJUqShZUzSX6Az75XF1h2wFWfxi1S1GiCtqvuAk3o3pN4wnZj8O9UI6rJzU+csjv455slHM9+hJMK2jIufslDYfhEaZ7EjDhknUqK9Jpqg14EgOP3C2/noBbfNP5ZqMqPgqkMgWh0fRFrJbQ9GZtDpSvp66L3g81fzwn+7etDDaJtiJeAEuYOnyh+qyYaxLjZZKc1IYi37i6sUSCuCwBUxudKAiAwB+Tn2X5Qksz5LHVYDTErzSgpihOuJWx0Krvi8MXMp71+3lhdt2kDFhh42peojQHCNIGhmEvTChQvSb1+7le/8+r4WxlLTCDLimOS22mq04pd51vTlfDf7KbS0+DNak8Sfbc9U2n5TkUbww/wnOT9/OhUzybvm8xS9Mlv31C8+pyuLKxK/lTyCc4HLReRb5vkbgHN6N6TekIzoKHcqCJJF51ImCJJRQ4EJbLhyZBiActkKgmbEq/BII4jCNpNhw1lVPBEq4cLvBYcQQeftJxDXGhJxcHAIResa5/h+kROnL+FY9zau3nkNcOyCx9RvFtvk2A7FyswKsFXvjPjVagYvd67hBOcOpsrPY83I4onQa6UxzWdE5BbgBealT6rqz3o7rO6TjKaodFMQhClL6TdDd3Vmz+VOv5elSpiIGoo1gmSEUMYIgk7uhavz78YhZLz4IlYOZ2cfC/FYIIOLL0oQJn0EZX45VOHza9bxyhQUO0uyb3rxL6q8IORRH76ED590LG95Tq24QjLBMM4xceOXJMDzoydfyH0FgN8vMqHXirN4BPi5qp4KnA3kRWT2O3WRUkpEdHh+hxNessRE6jSCWvNzT+DfV6+sbvOsIGhKHLsvCY0g6SzOmtvBa+IUbJWNspv1spdd8zhKk2ZJF4eAhhITQZmfj3hcNzTEHd7WBY9nEOyfXvyLqgPFaIz/9j9/qHs9qDMXx6XCo9fU8WfUN5sqL67P2oqP4GqgICIbgEuJqpF+u5eD6gXJSa7T6I46jaBJGGEayKiyLZvl7FU1QZC25Lh+ETd+STqLK4nFRPwjCnTh98Jd2Sz3ZDNU5ovcis1UiAllre+pnDRVVjrwWbTDzffvqzONLJQ0CYLGAo3JSELfr+8rjQR1megApanF5b9pRRCIqk4Dfwl8VVX/Cnhcb4fVfcqJH67X6YSXvOgd/PgHQTz0ZtHlzSJhLLXeA4KDG2sEiQk3vht8XXiUzl9sPIyXb1wP4w/OPZbquwmZWBAk8wgSi5ygA59Fq+yeLPOXX7mW9//4lo7PlYYop3EjCMKG8E9N3A+hCUyJF4wqwYzyLZXp3leqbYeWBIGIPAM4Gfhv81rqslTKiR+I33Hsf5pNQ7XM4kY6FpBLFDW1FsVxyLiRgy/5XcVfZUDnE1lQmtuuH2snIg4uUQE830+EjwYVYs+P3wdBEMfM33Dv3o7PlYYq6AeqgqD+dUlE3MXRWiaNgA3OTnQ6Ks82IcKDGZdwkUV0tSII3g18EPipqt4uIo8ArujtsLpP0i/gdTh5J5u+J9XyNBBPEu6MvGLwrSBoioa1zOKcW68RqGo1+iroQCOI2TP98HyjAYyz2MngSb15suKVaxpKHwIZJsvRZ+5Gi8zFlmTVjC/+ImrFMkMjSJb2qETJq/E34ongTu0E4G/WH8oLN21AEyVvFgOtRA1dBVyVeH4P8M5eDqoXeEG5Ogt6na6UEvdAsilIGqg6i5tqBOnSbvpFkAwfNRpBfN1VkxrBwuzkyWSjA/M0t6lmOSfKXVQSiYBRpruacfd+kTJeir6HShcSpNKgEdx0X/NOflKZ4oFMBgclNNcjXjBWRAiMsN5qeqKoKXnz8ESJ/9n8MH9zwuG9HvqczFV99ExVfbeIXMTMsjSo6st6OrIuUwnLVYNWp7bTup7Fffix9YKoTWV9+GjnJrOlSbLoXLYqCKLvSqmt/AIWNhkmI0oq/twVTKsJZSK4ko3CVhO+Hd+fRsw+QR80golSrBF0PouHqqxgkvIirYCb1FhGc/VTp+MXOWnTegA+bxZUsaZYFiFsqEDgGYH/tu/exM337+d5x6zjsJVDvRr6vMylEXzX/P9cPwbSa/ywUhUEHdtONb2moXjozdaufj9q5aSQah4BLtmM8REENadh1RSzQI2g7CUzg+cO4dVEiYmMZPBECBKCIAymqw2H/C74LOZjouRznNzJHeGRHZ8rVOV3hbdybfBY/ODlZNxWLNf9Y8pERuWp4DVcak1kd6tfWyRAZBpq7GQ45UUawUMHomvXDY2qE+aqPnqT+X+ViKwzj9NXxYpI9d49PcXIWMiU43Qc6ZPUCMLUaQTR2MMmyatBysxc/aJafVSkahoKjUYQhiGhScxbsEaQCBmdL6kvvvccopLYhHDvrr3VRU45nCbIRNpeP7RVb3wnP85/gp/pCcArWjpm53iJ+/ZM87Sj1tS9njGr5D92NzMdhItOEOybqpDF5/eF1/M1/yWonlTNAg8TpcDVaATxb+yykWFObGhgP20EQaxJxb6WQTHnNy0ip4vIbuD3wB9EZJeIfLSVE4vIJhG5QkQ2i8jtIvIu8/oaEblMRLaY/6s7/xhz88vrruVIZxujxr7aSbw3pDuzOB57c43AmoaaUTXHOA7ZTFRmywu96rZOTUNeJZHsOI9WptVKqFJNbrt9e219VtYSQR99BN5U5PR8mmxu+Zg/+fxVvOrrv57x+uh45Ij11K2WaVhM7J/2WEnkFD7Fvay+cJzWHmvVR1Bjr7e77lzFIBIMa4NdvMm9mKmST8kLeHD/YJI6ZxUEIvJe4JnA8aq6RlVXAycAzxSR97Rwbh94n6o+Fng68A4ReSzwD8Dlqvoo4HLzvKfkdlzHBWOjVETIqOJ3bM5JsUYQm4Zk5g+tH3HnaUQT/QhyxjQUa09hWJv+A2ldEEyUPH6/I5oMyqVaxIk3TwXYWomJWpbzcZnbq9tLUsaPo5j60GJcy/s45bCD2Vxo/bOPG79CY0XOrNEIpijgL8K2qQeKHislighyCeoEQSWxiPK9UrRASGjdcURewXzmadO+8rP6OT6S/R7B3vv4wuVb+OMzfsG9u/vfXGgujeC1wGtU9d74BRMx9LfAKfOdWFUfittZquoEcAewAXg5taJ15wB/vrCht86Nuh2Afa5LVjuP99Y6H8HiqhkyH3NpBKE1DTWl6iNw3KqzOBaaQRCgxjzgNwnJnY03n3MjLzzzasJQKSb6ac+rESSrj5qxPMHdUt1exuurIChWHuQ3hQIfX7dy/p0bj20wtMcCN483IxN3MVAJAlYRC4KQsl8bfzHh5P/11gc4+5f3EABrw8hm54fRdY3zd+JKtWs00qiKpWmzMFCu2dJ/C/xcgiCrqrsbXzR+grZqDYnIkcBTgOuBQ1Q1btOzAzhklmPeKiI3isiNu3Z19sVoYnWR1c4yQGFp+AiCZj6ClJm5+kVsjnFEyGWjPsWxVllX3qGJljUbN98fhSFOVnweHq+FJM5XyroWwQRZoxH847q11e0lJ6gJAu29IAj9aBW/kF/BVKNd3Hy2IaksStNQEMIqoxFkJKTs1YTVpF9bxT88uY9PXXwnIZDXaIr1/UhLiPN3PCMYFJgWwZ/ez3OnfsYd+TfgTGzrzwdKMJcgmMtO0LINQURGgZ8A71bVurRJjZY3Ta+4qp6lqsep6nHr1q1r9e1m4Pk+f9hzUfV5ViHscKVUJwj6sOrqCc5MSdCPTNRB8OUr7uLt37uJz1x6J9fePWNtMy+xOcYRh1zGNKw3QjPwFyYIRvNRnMZ40WPrrj21cwQVpso+19+zp+lxYVUjqOU0TDrRzzgfKlOO1gRBG6aqheIZQRCItJ0QNtVYnyisjXcxmoZCVU5yb6g+T2oE5UQ1WtcxZkOBnMQaQSTc4gnXNwUKv786wwlHbmJyaiePn7qWIamwfmf/G/PMJQieJCLjTf4mgCe0cnJTpfQnwLmqer55eaeIHGa2HwbMl0rZEdt33s1Vw7X43AzdUJmXqGloiWoEn/3Z77nkth189cq7+Zuzr2/7+GpooDjkspGzOA4SqNMIGo675YH9TJSaf6fvcn7E1sLfcGC6woE9D9TOEVb4xEWb+euzrmNrU1txzVmcc+v7Qw0rlEWJ37EfpiEviDUCmVGIrRlJYdGoEWiyZlK4+ExDYai8wv0lABV1KSU0gqSVwXGM2RAhbwIz/bBMxQ9xzOePBcHlo9H2ieJO7nUC/m31Sia85klrvWRWQaCqrqquaPI3pqrzmoYkiqv6BnCHqn4+selC4HXm8euACzr5APOxZ3+tWfSZzzuTrMqC471jkjdz2Af1u6uYofvOzHJR6TNzzU83qmLWGtO4FHLRouLB/ZGJIAiaawRlP+DlX/4Vb/nOjU3P+Vr/JwBMTOxnbKpW0tjXCoft/AWb82/gru07ZhzXzEcQk1EhEBI+gt6vqithZBv3pLWkMlV4u3sh38p+polpSHnF+kM5a+WKRakRBCYw4LR1a7mhMFSnESTzCEQ8IHIW50y1Wj8wgsDsE4ewO3H2sbeXK4YO8M1VK7klnL9TXbfpZaDuM4kczs8Xkd+av5OAM4A/EZEtRM1uzujhGNg7vrP6+Pmbnk/GNPzuFmkzDcUagVeY6dwLU1ZJtRVKjZk/C0ATpqGhXGQamipHtny/ThDUJup4rXDzfc1LRkw5YwCUD+xm2KtVHPVDj5dP/4jpjEdp600zx1JdhEg1gikmp4IvSqXDvIZ2KGpk685oa6v4QJUPZH/A89xbZnQkK4ZF/pDP8cU1q/CDxfe7ckv7OOA4XDI6wjsPW00pUS01Gap7cvbn/Cz3AQKgICYBMazUJQ7GgiBeO5S9fVSMLlcJ+5/Y2TNBoKrXqKqo6hNV9cnm72JV3aOqJ6rqo1T1BaraednCOdg/VbM8iQgZnI5tp2HiB5Y2jSAWBCLLQyNIlm84LfMDnuHcPsfezYknX8dxcU0ewTBR6GeycbwnNfNI4PtsLfwNr5P/anrOH6wY448P30h5fBd7qUWc+OrxnZUFnnf4Rqb2z4zNr45FhKxbqNtWQKiIVBPc2vFZLJSycRfmVPFayI5NFmubKtaX03jAj36rWVWCyuLLcs8VH2bc+GM8EbxybfxJE3FZhMc42whEyDm12lSVShnPaGs+sUYQMR1MEhsX+9VHIsniSt3rAePT9U63jHZXI9A+pPE/sHe6JftrK8Rzw8eePDN9I+zDZ+k3FT/k5c41fChzLv83cyHn5f657XNUSz/jkMlE/Z2P53cA+HFlUtXIPGImw7ikwKnuD5ue8yurHCZch+LUQ+x0kj0EPG7PRNdhV+XuGcfVSkzUkttickAlEQPQLFek28QrYU189rkIEyafoFifbTtlVsLDYTijJMNiIFvay/5EtrNfqgmCpH+jHGcbA3k3jjKr4FWm8aqhxnFuSkRRi3gSvearFQRdZ6oUOV4+9fgPAZDB6XilVO8j6K36PV7yePa/XMHJ/3FdV84XawTrhw+bsS1t2k0rlP2QL+S+wlsz/119rd0Y9XjydRyHTC4SBHG/59jBntdolVgxLSx93+MNhx7Mr4aau9PGwuinN1V8kL1OmRUBDIWKT8C6IJrgd4czzUpJjSCfrRcEeYVSog91P3wE8eLBQ/C8+Sew5O+lsSZ/bIIriRB4iy+CTfwyB5yEIChPVh8ngwZ+PDbKN1aOEYhQMA59Tz28crFqtvOlwSymJcpO9FqF/ptol7wgmKxECRuPPfiJAGTE7fwHEv8YVatNS3rFdDk6/3X3dNeC5jAzfDStlVTnIlnMa7/j4NF+XZe4nLDjOGQkcv5dNjwa9SIwGkFOwRehZLKES16RG4cKnHbImqbnHFMTLVLehUdAQaPQ5kB9hoz2N6Uzyw1Ue2GIVENZq5/VHakmtwF4/dAIzP3vCXiV+QVB0rmulfqoqFhIlB2nWq9nMaHqs9+tmVSDhGkouYi6fqjAmWuiyjmFTBRcEIQeXqmI36ARVKqmogpFJ37c/9/hkhcE+8u7yKiyft2RQJc0AvNjzGjvQ/SCNmOz56OWA9FEECxF05BXW109+4iNvPOQddXSya1SbV6PS96s8H5byBCEWp3YcqYdVdGscv1qKeImmXvAiAm8m/T3ExDiqkTJjgTViaCiMyfDqnYiQi5bX7bYb4giapY02G3iCbCxHPZsaBBwey7LeWOjUK7XCJKTqb8YBUEQ8LVVK6rPk4JAZ9Gmh40GqVqhlPi8vvFTViuUqseUye3xrCDoPjv8hznMg6HCCAAZXALRjrohxUdm6X3UUBgqo0yzmrlbGLZObFpwOHrV0XVbOg2rXYwEpQl+NjzEv6+OoqSuGR6qthtslZqz2MF1XJ4bbmBFGOKHWvUfZI0gKBlzQdy4ZrZvdMRoBFPBBIGEZBET0eZXY9K9JnmbVR+BCPkGQUCD0Gmn5MVCif1tgQheef6CaWEY8OoNh/Gpg9ZAg0agSbPRIjQNBaHHjkytYLOXmNhny+IeyY4hCkqFYmJ/zwiCahczfKZjQSBWEHSdh90ih+pI9XlUw73TRhrRsVmVugiiXhCEyrX5v+c3hbd154RVy4LwvZO+V/9eS1AQ+OUipx6yjrNX1cJlJ6faK+qVDB8FyDtZ4w/wqrbhHLEgiM4d9yvQWTSCuPRAOSzi12kEYXVF2GxlmMwjyOfqTUOjUv/cW0C2b7sk75liZe6mOlBzvAP4lXofSFIjWIymoT3+PjwRjs4dBEDFqy3OZvOvDedGyAKKR6mS1Ajqy8H7+NX8j0ofwn4bWfKC4ICjrJCx6vOsZPClPj28XeIfV2Qa6rEgUGWFFOvetxOSpqGR7AjrEnON34eSBP0mqMyc9IsT7flbqkXnTMhtrTNYJaERRD+lohEEwTwTWSYuPmbKRmcQMuoQECTsx01WhtWxOOQzNY3gMB3h2MyG6nNXo1ITva7Zk7z/K5XJOfaMSNbtr3j1gkAT5/LmadAzCPaYLOqjh6NOZKVKre9wOEsOxXBuzCwYPcqJezH239T6GtdEaqUPvp1GlrQgKJanmXQdRjI1QZBxXHzptCOQ0Qhw8HssCMJEYk2lGw3CzX/HqKHJ9epS1Agq5ZmCoNSmIIiLzrlGI8gajcArl6rho1nTGabkmUzbYL4VbXQlypQJRMngkEHwCROCYOb1jp3FAhTyw9XXn8qh1cqoEEcxgdfBgqcVkjk5pXIrGkFif7/e3JlcVZe9/pdino+4LMSaQuQI9vzaCn+2HJx1Y6sp4OCLX+cjqBjzdE0QhFV/UqVZQ/Ees6QFwfad9wAwmquZBbKSxZPOJtX4MmVUCHu9ii7uYWsmw7QIpUoX3qtqGoouvas1UbAUNYI4fDiJN92mRhBPviZ0MO4D8KZv/ZKi6S6WM4KgYlZ9yYzjuc5ZwcMXxcUxyY5aDS1saiJImobyNVOQKy5ZpxZOGoezepXehiImNYLpFgSBJr6XUlCvQSQ13vladg6C2A8wZuaTil8b/2ymoaH8GEO4VJyQKXNvjIRCxcxBsWnIS3yTZasRdJfLb/4+AGOFg6qvZZ0sFREqXicOmVgjkJ7Hak+U9vPSTet5ycbDKHVldVdbUUL9DdBr7WYQlBtaBAJMTe9ssufs1PIIosk+a7JFdx7Yz63bI6GSNWGlZdNqMpynkmu8FvTw8akJgshHYLQFJ5zh2NaEaWg4W/N9uU6GrFvLWciqROGsLdjtOyF5/7fyXn6iWmu5ITw2mWPgz9OycxDEWeRjQ6sA8IKa1qKz5BMNZYYZIhsJglIkOEbUMXNQbfKvOFrVCMoCXp/7MSxpQfC1/RcCMFaodcPMGLW+XF54Cns1fLQLoajzUTR2112ZDMWu9DWtlTEGcElqBIuv0FenlL2Zduvpcnv9LWohm0YjMCaYywqnIkbjqAkCYxry575W8T1UFp9AlCxu9X6KI0oqomy+v744b7IMdT4/Wn3dlQzZRF5BXuud170iGSxRbOG9khVuy2G94EhG4HmLsG1qvOpfMRT1f/CC5nkESQqjhzHkZCk5ynSsEZClLEKpXKlqBEktoOQIpR5rco0saUHwxtE/BeCYw59WfS2u/VEqze/Ymo2kIOi1RuAnuoaVW1C95yMerRofwVIXBJUmtuZSpb0yv/HqPXYW54wgqIiw5kBUuyhrejXFTk5/Hh9BmFj1e4LRCFwjCKJ99roOpYdum3GkGQyZRImJjGTIJspSx87rTu7zVkiWsShV5l/Fe4n72dP6xVjSNOQHi08jiBMuY0Hgh8XEtvoVfJxvMlRYyZDkKTlKyWg5o5Kn5Ajl0nRVjMaho44qJRGmp2dqsr1kSQuCv3/l5/nVK6/mhCf8SfW12KE2Xe7sByIaqfO9rueSnFDKxW4IgrjoXHzj1fsIeh1u2G+aOR1LwYEme85OvAp3qdcITjx8AwdMJEnO+A3iH3swT1/smkYQVavKSIYMbuTglVhIOEzsqy+SFwcBOUDGqcW0I5lq9zRIOK97qBFs2TlRl6vQStRQcmHTWFMnrIsaWnxF5+KksRWmcm9SkMWLhRdtiuaalflon5ybYzgzxJTjEBjn8og7hCdCsTgxQxAMK4QiTE02r1rbK5a0IABYMbK67nm8mit38ANRVQSqzr1eUicISp0LAqkKgnrTkJhww24Vt1ssTBRnrqyKfnuLgHgCcLPRZJ9LVP28jz11r1ViQTBP/+f4Wy46kSbmikvGhDaHwMF+dH12lBt7EsQagYMjDqvCSBgEPnUJZrHzurhALfLOHePsnpzbPHPt3XsIRMmZe8ZvIdInWUzObyh7XucjGEAFzvmII4NW5GJBUE5si+6R9zz5Hdz6uls58fATARjJjjDiDjElQmiipEbdyKQ3Ob2/6hcomkCEIWMrmpjub3OaJS8IGskZdbpYXrjqpUTOVrcPPoIgIQi8Uueru8YCE7EgKGg0IXWWaLf4mDCmkXc+9vXV18phe99jaCb1rCnyluwM5pmQwrimTKUaPjqfIIgmvWknqlGUwY00AqK+BqtM5vGk31iGoeYsBtgU58h4JYZyNedx3CKx2MIqvRl/duYvecHnr5pzn9UjOXxRhs004nktOIuTCWUziqslchIW6CPYPVnmeZ+7kpvu6/5EGpt/Rs1q309kfsfXJWNMz6cdfxqXvfIyVuZXMpIdZcpxUOMTGc1E12lq+sCM8IwhY9KbaDOyrVOWnSDIG4dauaNoCmMqkEzPI++TK8tiufMyE8nMVICciYTJGEHQjVyFxcSUuc7P2fAcbvzbqFtYo216PuJokaogSDhlK2ZVOGR+3GXTxDycRxDE3/KU41ARYxoyWe8hUDA+h0qjrVzrTXvrTSDEykyJoUItXyYnnfvC9k/P/Rn8ICQQGDWO8qTzdDaS+RWNgiCpEXjBwkxDf9gxwb27p2btDNcJsUZQyI3gar0giLXGrBt9Fxknw6EjhwIwmh8lEKGC8RFko+s0XdxPKEI2YY4tmNaWU0VrGuopsUbQiSBQjGlIIptu2ENzSrKhfKnUrXpDUd0cgNFcdOO5RM7Pfoet9ZIwVIrGFzRUWEXOyeGqUqG9SSbuT1wwK+6kLd4z5o2hXKTux03M5zcNGXOKCEXHIedmyTg5PNNqsioIGrpVVX08RpAfsTrKJnZHhJGkIDA5BaUONN/58IMozmeFWQW3IgiS30vQUIo56SOYz9k+G66xte+dqnDHQ+Nd/W3G48s4eQoKJfWqmke8zW0o/AcwZnwKFYmu5YgxLRVNxNlw4icXX/fpUnt+rE5ZdoKgkK1X4ReKo5AV15Sr6N3kmXSuFbsgCLQhyim76ojogZM1NZiWjiC4d88UBaIfVD6/AhFhWKVpMbe5iAXBcCHK5C1kaxm98apwyGSvV8xK1k84i5s54MOG67ByaIScm6ciUd/hvJg69nBU2rwAACAASURBVA0mkmT1UYCDD382ALsPeTQjJr4dIO9EwqrUJHy2a+z+JbszwtqhQwDwtQVBkAgf9fApeQGfvvgOJst+fdTQPHkYs5EsqfGiL/ySc2+4f0HnaUbsB8hkshRU8MXnFV+9ljDUqjbjOjP7T6woRNelLNFnGskbjaAc3ZuFRMDGkNHkuqH9t8OyEwR5UxbW6yCFveosNs69bvTFnY1k3HWl0r2bo5pHMBr9iB3JUO5DJmo/uevhSUacaEVcMOacoVDabvzRqBEMJUo7lE09oLjuT7yCT0YNlbzZS0XEHLTuEeTdIVSEKceprgy9sFEQxGFD0fU76ZEv4Znrn8lbnvIOxoZrJZKHTSe1hUYNnSB3sFHmzrd4cPIWAhFeeOTLgNY6ayW/F19CfnzTNr5+9T188fItdbWGggX6CBprK939cPcEYSwIXHHJq4NvGslMlPyagHZnCoJVw5H5ruRG+4/mjCAwyY5xAUKAgikcOL1A385CWXaCYMSoaWWvE2exySNwc3gilDpITpuPZCMPr9wNdTEulxCtQkZMdmo2rsraQix4WtgzWSHjTiHAsFnFR+n+7QnuOH48b0xCw0O1CbfkxJNDhlyoVedxcsKbmJr5o27UzI7Z9McMZWoCJtNwvupxcUlsYxoay43xtT/5GhtGNzCa0AgafRbtoKr8MP9Jrsm/a8794kCGjaNRETa/Bd9LXR5BoqTJeMmrE46Blnl4osSR//DfnH/ztlnP96qv/Zr/e+5N1eeN3efy2e5NcbVVv0sBl5wzzVvdi9g3XamZhpIhvYYVI1FzoqIb4qoyZDSCOOu9QK3ZzZAbLSiaJUL2kmUnCEaHTMGoDgQBxkeQNap8uYfJH3Vx15XOBUEtjyB6ftrxp/GGx7+Bp2ePoCxS14c17eydmKaSKbHWyZM1KnuBDGVpUxAYjSD2L60YroUkTzm1XgV5JSEIatftwIGZK+tGQfDoQ55MIVkyQty68804rkl566QgGI59FgsQBK2GEMefcTg/SkYVT+bXtJIC0nNC8hmHp8ofKFeCOo0gDCvsOBAJlm9cc++s57th614uvrUWYhtpBLXx753sXjnrOPPZFZc8GYoifCh7HvsmJqrO4jj7PMnoUFTiZspRXI1KUwOUTBhznprwGM7E2/pbdG/ZCYIx8yMutxlLniQKH1UKbiFKDJnunT0vOaEEfvcEjmNWISvzK3nvH72XgjFLlNtUSf/+vN9w7vX3dW1c3STcex+7MsIh+drEPUSWcpuJcz5+lEBoQjJXjdVqV02axZxjJu6KcR4n6+7vG39o5tigLlpk9fC6qu0YwHEyUfvLhp4E1VpDTX66mUzNUTmSi7SWhRRva7V0dTypF7KFyGbegu+lThCIcviuKzk/fzrH772gmrgHEGiFXCb6jHvamMxzB+5ja+FkXuv+nG9nP8PwvjtbPnY+QgJEFUcccmTY47pRQci9D9XamTYRBMPD0f1ywHXIAMOF6NqU/WjRVaB23YazkQDvd9G9ZScIVgxHalqlhQiH2YijhgpGek/2MOY3Wd620iAI9k9XuPS2HfMm/iSZbf6LU+Lbza94cMtv+O2dd7V1TL/IH9jCzozLIcZ0AVCQHEVH23LwBxqQ1VrI5tjIuuq2uJm567jkVarNZJIT3r7x+npBEN1DI4khZN0sw7lEuXRxyZnzbds3XV2lV/NAnJkagSS0hNHYBLqAMMyWBYERejkjCFrprJWMGqqIkpncDsBhxXofQUgF3+S07Jlq/f4uTETawyez3+a57i2csO/Clo+dj1C1asTJaIZ7c1leumk9xb3bCQlxVOuuQcyoKUmx33HIKIzEgsDkFRQSVWPje6Dc5xIby04QDA9HF6UTQVA9Vy7OEOydIEhqBKWwfrX+5E9cxlnnnsc7v/mLNs5YM2UkiQuWldqMVvhJ+B4+tvWUto7pF9l9d3JfNssRa4+pvlZwCxQdaauAX6ABSctvfN2h1pPYEZecOlVBkBTgB6ZmCoJQYDisvwax/yo6X4Y8Qll9nvWZK/iXn0Ur25omM3dD4iETouiH7QuCoEUhGQu7fHaYgjp4LWTVVI8JoeLAlO9yYyGPBGVClCEjhDz1qjkt7SQ5etTb6MeDmeGcCyUkwDVDcbX2PuMH7kE1nHUyXZGPJv5QBBcYGzbhpEZID7s1k+BQJrq3FppHsVCWnyAYMhrBAn4gMXHUUCy9p0u9S/5IrizLifC8XRNlsvicnz+d9+/9WPsnbli5xFE1pQX4CEa1v46tVpgs+5TKv8cT4bEHP6X6+pAzxJQjTLfxOX31qx3FgLpyzzG5TIY8bjWKKEiYhu7asY1ipX6SVJQhdeteG0vY+DPiklWXiplcr7jTCJM4OmWen+5wPg5nbT/6ptXG8XUaAS5eC074eGEzpEJZ4FsP3MkbDjuEa/K7gRBHIyHh4+G1IJAy+GSS5jPf40BCW1rh7W7ps7RCSFjVCAoJp/S+qfuMRtD8uLybr5oBXYXRISMITJTVSKa2sMi5eTKq1W39YtkJgpybx1VtO7u0HkUURo2kn+qTICglqh3+YecEq5jgv0aGWZ1p3Ubf6KSMGTZhtVOl1k1DqkoF8OlOG81usmuijF+I+g48du1jq68PZ0YoOg7TbXQpCwjrBEEzspkMBbKUTCRMkNAINu05h7O+9eW6/ZVaOYGYlSNrq48dJ0OOTLVJjV81DdU3yZmNQia+zxciCMpcNjzEluxMgZck1nqy2WHyZJp2VGskFpBDOBTFYZroftueKRISmV7yavo0JExUszmwb8u/if/Jvb/aevbByoM864hN/GhslE+vWU3gdFEQaFjVCNYWat///vLDKIo7y3EQNaOBKHFzuFBfq2goIQhcJ0s+ZEHXrROWnSAQEQpKRxI3viVjO2yxi/H9jYRJQSC1MU+WfZ6duZkPHnwQbzs0mkAuuuVBzr76nnnO2Nw0NFyI1NN2PosfKn901OG845B1MxqoDJqK57OnMMUIGTaNbaq+HkdsHJhovSdBSMDMoMB68rnhqv8B6p3F/7J2Ndu8LzacU+tKgAOsXVFzQmecDEOSo2RmntheXnMWz20acp1c5LxeQPG2wPN47yHr+MuNh829XzXBqkCeDGVH510QxAubIc1QcoTXONcAcYcuxQHyCJ4EVPyA17iXs4JJ9sziByuIx5HOTiZL0XkP+FFk3ScPWsP3V45xyeqJriVJhtTMP46bCA/29xLOYRqCqBkNRB0B8yaUOc5nSfoIXMc181P3op1aYdkJAoBCCJUFZi5CzVkcVzadbtIFq1skV5blRCJUsRLwgfy3AXjYRFf8/Xm/4Z8vvqPFM9dPJHHoYbmN7OWKadN47fAQD+9fXOYhb//9XDOc49HZw+oceHGbwb3jrXcp8wnI6NwT7+jQaoakwLQTaUfxhLfWieLCrxgp1O0fSi0XIObQtUdWH7uSZUjyVcFSdRZXo0fn/ulmMtmm4aetEHitCY9YI8i4efKSo9xCln1sGiqQwRPBNR+oQhDZ2TVKsPIkIHvgPj6d/QY/zP0TO8ebj+n7Y6NcNDJcXYhow+96WH2mutLQyQhv8/1Pj9aEdlEnqs7i2Rg2CYIu0WQfmX+iMTtS0yUybpacSnVbv1iWgmBYHUpt1pupJxIEB62Iokemu1gDqJFYlS6EUExEZRTLFfbHESttWGWqiakNl37ECDWvDWexlyhxvG/X7Ek/g2DH+FYezmQ4YfiYutdXmeSe/ZOtC4KAcN7veHR4FUPuMFOOUCxVqgL8I0/4EACHe/Xft1ITBOuHDgYgV6glqmXcDENu5M840bmJrGmL2NhPYjZcJ8tw6FBmARpByz4Cn4wq4jgUnDxFB8pNsqjrjzF9f01Oxm8L0f9QAgLje8vjUpEQNWVgjnXub6pxhqHy6YPW8KGDD2LvVDTmxpIcRTdgotQdQaAJH0Fp5Ybae1AkVJ1zMh02pSPiHuF5hYrMFASuk2VYXUptZr93yrIUBCs0x6R0ohFE6+mV5mZot+NVO8SrrjF1KCaccTq5kwNudPkaL+Lc6nl9P4KYsZFohVP2Wk9aS/ZHmNq9veXj+kFclymu9BizZmVUEXJyorHO/+wEhGQaVu8/ePEPePnYH1efj42sZjgzii/C3gM7q6ahNYU1PKM8yoTTWGBNEYRrXn0NP/2Li2a8pytZhtwRphyHs3L/ykf8LwG18tWN2kQjWTfLiGYotpDk1YiX6Bc8170UaC2KJu/kKTpCqTx32GOsKY2Zml8/WhFdn9Dx8Y15JRYEfmmKi0aGGXeE8eJMgeaFNaGze8J0hmvQCB7OOExMdkdjDwir3/v7j38/L37EizkkcJmSMqEGc06mo05kkowbQY2EwrRZ2MW9sCHSCEbIMu10R3i1yrIUBKMMMe52YjeM7v41Y5ENdcrvpSAw9Uk0y1TCBpuZeJAvr47MHE6D2aJZbZuYxszimJWj0Wcp+q0LgqRjubh3kQkC4+soJEI9AdavPQqAA6V2NIKZ9vzHHfQ4nrCq5oReMbyKYSN0du3bXlegbIUzwl4XwjAZJx/Z+VfmV1bLXwBk4uvrZhkx55tyhMdo5PvRmm2o6VgLxoTkulmGyFJ02r/Pi4nufePF2SekUAMy5n4azgwxLcL4+NyBE7GALOTX1L1eckL80JiGJIfnKNsnt/Ghgw/i5MMOxds38/7yE2Glk/v3mNfqBcYe12Fqf3ccxoHWTEObxjZxxrPPYK3mmXIClGDG7zDJymyk7WWNj2dEXSbcOAktqRFkGKbA9AKuWycsS0Gwwh1lvwvBAp1IsUaQd/OMBcqk9tBHEK+gKDDpCuWKx96pClu23s3NhcjunG1YtU2U51gF1jKS6l5ePRyZuYpB658laRLbPb64ksrKplH4UH6k7vVD1jwSgMly65ODT9jUR7BqpGYnHiqsZpXJUdm55/6aIHBzrM6u/v/bO/PwKMps4f9O9Zp0VrKQQNiMoAgIKO46go7b5Toyjs7oeF3GGWfmuo2zeNU7flede12+mcFt1M8FEXcdmcUFFRXhQ9kEZAfZAoGEkITsa3e6671/VKXTnQWSNit5f8+TJ1XV9Vadrj71nnc557wEDOFAZUsvREn7LqC+5orccIWdEa7NziLfbelBONdQR4bA/jwUMokXD3UxVCiNET2CA1Udt/BDhMLeVAnuBJQIhyoOb2DDeZsShkcdP+Q0ENUUfq8CYlJl5+Tf63ZhVrTNItoUbOllrdy0hdF3L6C8wfrd5854iu+nnEK5w0H1EWTqLKYKtWkQJEocNQ6TgArgOYwhSPVaumLanlXxykmVw/YkEkf4HXY63MQ74qgxeja9fWsGpSFI9qRTbxjk7evsxGprVFgdUk0HNfRMfp6VeWXhFyfREU+1YVBdVc4vX1vLvurd4fO8yoxyr+vMmGjreiTRnYihFPVm579LZMuxzN+/0kw0p1+O9yRFHU+1s63WBjvv8hsSFV6vOJK0pKzwttPhIiNlJAAlFXvChsDldDPEZ1V6u/ZtCJ9v9QjaEm+7GTodbtLsHmee28XDGZbRb4m+bd9ZMcFOX+32CPGOeOoNaRPDcCQaI9KMNOf7aY/IACufbQSr2qmwo8rYDZtJSdFzN+UOBxh+DCwX31rDoKGhxViX1e1sc61AoEW24qItzDRWcqjWMgTHpI5lVNJIQiJUVnScq6grCME2Q6pJziQqHAZNNEVlEW1NZoKlK825qXx4qLfn+MRwkNT8EykH8Y5Eagyhsqb31iQYlIYgN3MSAOu2fxpTeRWR1yoZL1WOQLf70a/MK+Oq51ey55DVQk/ypFFvGJQU5ZFXWgduyw/+BFKpdAgHS1pam4czBM1DQ61zojgMBwmmUKc6n+yqPiIdRUWo82PuvUFzHpd4b/QcgctwkRhSlAUqKanpnMNAEBNHOw6kI7LGRe2PzBwLQHltYXiS3+Vwk5VqDUcVFLfkvbHmCNq+fsn2GsRKGeRkHNsiQ3O68+YeQTspJgCeu/QNbsz5LsePOo0kdwo1hlBw4PCV8/Yv/8GcOU+HJ2T9ET2CosMYgpBqmTzNTLGMXVnl3sPeq3myOMOTzqIrF0V9Vu4KYSCkuFOocBiY/hYX34qGtg2NQISzwvWed3ja/STjxGogeVwJLcOANd1jCJQEMVrpQYo7lRqHQaPRFJVFtDXD021Z7F6Az9EyHOgwHCTaUeZ1DTWkeFNRIuwt3N4tcneGQWkITht/MQDbipYf9jzTVPzH/A1t1z+Vlh5Bjmc4e91CflH3LYABcKCyefLLqtTTk3IAOFi0lexkLzitSnhS0iQaDIPtu5aFy24sOHJrt71KKAEngS70biLXwy2X7m29BL7lYj/NlVnkYi3NpJpOdvlqeOyZxzp1rXrDDC8YEklmSk7U/phh1pxBQdUBymqbDZGX0dn28UM7wueqdtxHwRoCBKiuL2NM9viW7yPCjqLyI3oNjUw/nl+f/xgiwrCU0Zgi7Mxf1fGXU4qxn93A5QfuZf7ir4Do9CulFR3rkhkxZDYi3TKKFXX7O74X0Yu7pHnToj4rdhqIgrS4dEIi1AZbhnSqA23TdNQ3tOhcyG1t+wyrIeNxxZGdZslU3fDtPdqUUjQYJr5WepCRYHl8lTjBfZhok9zhJwLW0qQASe4WvXSIkwx7mDEjI57s1GMA2Fuwgd5iUBqCY4dPYkzA4DO1hTtfnMnDr/+EZ/5+F298/Cc+WDqXL9d9wLa8NazavoWDRb/gT28+EFVeqZZu/dRhZ9IkwodLn+pWGVsmwqz/x2VbilRQvo1s135KkvNxKcXpo6cDsG3/MlKp5hHn8+xeE93SCoVMisqro67XXjWSJUkcdAepqO1cS7kx0NJ7OORs7LagsvV5RTxw329ZsS324aYmOwq7OZw/El8oRJnTQWnyXKraSWj2zpI13DO7JRK4zoB4Ox4gktaVcUrCUNymwu3fTnrDVgCSEjOZNPYcUkImu6pXh8812ykPMNTOVOl0N5BlVwhgGYKyXWsj1pw+8qt73LApAOwp6rhCCdVX8pvMdM4dlUPtN/Ote0X0CKR4S8dlI9wpR2QeB0BZ/YGoc5bvPkRpTcszbp7zchhuHEZ0C7rOsGZNRtrLb1aZLQ2w6lDbSPCGCEPwUYKPSWNGss9tPRen4SRniNVDq2hqO0dQHwhyz983UlrjJ2SqIwadNQSCVDmsOYFIxuZYxrrM6cBDx5HYw9Oih8KGJY0Kbxti8Mj3XudHOedx2oQrODZ7IgAHKnqvR3CkgMmjEhHht9Me4PE19/GJKx8zuA9qgPbmSZOApEVMf3EC8abgVQbljpYx18u/cwsvvvIyL5sfsu7ZL0hzphPnTCTelYjbGcfOqo1sdZQwMZTJpZNuJjFhKLuK8shOSiA9KYuX15STnJTFXTNPxutqeTEO1flJNw4idoqBqWPOhE1wsHIzAVnKLo+DRFNx9oSZuNc/wPbK1dzuLOQq5xLyD22jrvE6fF5LMZ97+y6W1L7Hzdl3UFdfDAntpyg4LjmXt6rL+GL5Ar534Q86fH578/eyc/Un1MVbD2xUyEWex2Tdts1MP2lqh+U6y8qFd7N19EpyPjnAGePfBqxhkcjn05pA0MRUKnxO86LvifFD2pxrxKeAWcmqOC8fffwsV/0gegGWA2v/jc1JUFl9A16Pg1pDiCO+zXVaIyJkhxy8n2LNACSGFB7bTXKKymaVt4jSqkNkJKdj0v5k8d2z5mB+9HOuv+h/cEWkla5xGAT2L0fZUx4dTRZHMvnYc2Hjg+wrW9vhOfl7NrLIZ323UpZSVd9EINjSEEipWA9c3W5ZEzPsF5+RmIMoRZ2/pfXd2BTipheWUIeXBbefw4RhyeHFXVxOq+q59JhL2Vy2mb2VeSixzNuwzONgFzQ6qwEHbgV1RtueSX3EQk3N32GhLx63aWUBTfOm4TGhxmzrGLBkeynvfbWD6sYgZVX17C6pYvV9M6POCZmKyvoAaQkeamtrqHQYjHdGe6GdOHYmxvpHMUXwtNNrbCbOFa0/ozMmQMVCwHIMyEzI5t7znwBg/KhT4Gsoru6eIa3OMCgNAcC5U2dx7tRZ+AONFBTvoqB4J+U1xVTVl1Hrr6C+sYr6phpCBDBwUhIqwk+ABppICwkjDWvyx+OK4/EZL/DE53eyy1XBKkctSsRKwBMEr0uRGRQWew6xeMcfWgSISFHvrFQsezkECA4lJCgDU4H/OFhnn5OTPJq0kJPP3fuoNwRw8F/HX4/XHc+JZirLvBU05hzkNdcwLqmrx/fCVSTmTMPldPP/mj4Gj5v7S/5CWYJVUbaXN33W6b/kzU9X8+qe/8Oyp5/G50rDKW4azRrinakkeNPxuhLZWPo2eR7FqGIgwcWZSSfyZt1a/rnsFpatHE2cN5WrZz1MQrxVCbqdBh5n20r8rx/PYcywCYweOYn0ZF+4hfylbGSX2403tImfBIL88R8r2LltPU/c+QuG+Np/2a57YTl5ZfV8de+FAARCATAgOSG1zbkPzXyRVzY8x3v5C/ny4KtcRYshUEoxJ90yoNt3f8XwrFGERPA5EttcB+C1GXOpi0jLMZxk8rFasokR7qLfGXEhS4pe5e1PZnPrlQ+jpP00EckJWfzfK9tPnbypeClu3ynWzhEii8Fyb84Nutnm3s8r73/IlKlnkl9czq6v32DKtMs55YRxzF+7IHz+iuR6Xl+0gpKyQ+Euo9RvsJIstmN4InsELoeL0SqOAnc56/MOMCY7g5c/W8tlw+8hJ+Tnb3OvIO3WBwnYaajd9iLvD53zEAAXvziFQqflgpk7/HQAttiN7+NVEkXucnYXV5E7NJlDtX6e+nwXJyVbzzlXJbBbrGHKoEjY80pEGKa8HHLWUlLdyJe7DjFrynAMQ1CHtnJr1q+p3HUJY+JWc2ZiNbPf8XD5uaewuaiO8dlJLF2znq3L3+E/7nyEouICag2DVHf0UGOmbygTzWQ2OqpxmIevTh8791GG+qx6Y8rYs2DHo4CVVyqSlMQshoTgYFP3Djcfjj4xBCJyMfAEluvDHKXUI30hB4DH7SV3xERyR0yM+RoTjzmdF475ArCWlqysKaa0vJDymlJOHHsWCfEpLFz5GnuLt1Lrr7RWolIOapuqSXI1UeKvID9YgZhBFE3UErB7AkGGBhW5koRhOPj9aQ8we+V9BKSJq93TuPiMOwH4r395jrsX3ECxw8GQpJG87NwL7IDqHVFyqiMEIZ0w7FTuHH0Nr+1+iw+9pUBkPp49YeNGgpPUECy1G0cXTbyegrVFLEo6AGwCYP4/zyCu2SVOQbwJLmXgUBJ2wdsQp6AY0leHOMbvxSNuTBS7PQHAYHOcg2vmTWKHx8nwnCD/Pe9vDE84FsPhwSVOHOLgQN1OSkPFTPZXMibRwey5Z5M76kyK/fsQr8LtinYfBThmyDjunzGbsnlbWRG/j+f/eh/nnnItGenDKCreGz7v0Gc/ZsmwqwBIcLcdYgKYPPKUqP3jUyezvHYJAAciejCzZtzGE6++wqqyD/hRxV12j+DIrfofZ53PGwcX4VZCg3c/+3anQlbnegQAN029mbs3Pc5rJb9lxXsm+W5hv9vF2DWvMuoLD58lKwyleGjyr7h745Ns2Pfv5HtMcFlVw7aknTzw6A/xuLw4DSeGOC3PGXHSYDTiiPCUuWjUBTy7/33efP8SAsrNsuR6GpK8gJehiR/hm7OcQncTJMPQtMwoOcd5sigMFdJgBEmLT2eU6SHfYQ0pnZF9Ks8Vf8b8N29i2jn3s3TrbhyFj/K1KMiC6dlnsvvgJ+Fr1Tlans0E32gWyVZu+PMfKU0qIK/kVn5z4SS+KXyPeakpnOH9nBVxcUA8Nxf+hEXPHcNaczzPBcYwLeNdFo6tIfe9F8jJsPQoM3Vkm2c8M2cGG4vepaidXksk3x19QXg7O21sxCfRv6WIMI4h7PKUsmjdNs6fOp6eRno7a6SIOIAdwAVAAbAauFoptbWjMtOmTVNr1qzpJQkHPvur91NQtoOSg3nUNlQzLmcyU475DoYYVNQewOVwk5yQddhrBEIByutKCTTVkhSfSUVNAaUV+ZRVlzAu5yTGZJ3IwrWvkBSXylkTLgPg692fsmrbx9Q01VJUU4g/5MdAaFJB6lUjQRUkiEkQkxCKfc4mkk3BLwqfKTSJwlBCraG4KXMW35SvZ7t/P/tdJiNDBgVGCLOdCtAXUlEvfzOpIcXSGzd3+B2/2b+Gn316Q9iTI860Ilvr2hk2u33I97np0j+0Od6a2oZKbnt7FmukjJOMNF6+dkn4s9n//BXzqj5naDBEsdPBeU1ZPPGzw3uuNb+fd777QxZWtXgdvXH2c0zKPbOjYlEs27WAl1Y/SaG/HJ/hZruqJk25qCeIieLqzHP5zSV/4ZnP7+L5/R/hBI6XZDIzxvJZ6eHfu8nBeF77qTUZ7Q/5uefd6/i0xnqVM5SLWoL85+m/56W1z5AXtMb5jwvFMf/Gr6Kus37f/+faxbcyKZTAGzeuYMm2d7jtK+t5L/vRF1zx1nkUSROpoRCC7W7a/CwumseXO99jbdlmVlXtIFvi+OQ66/rr9i3lusW3hM89saGJ5KCHg64gO72dmyKd1BjAj5cdXpN5Zz/BybnnRX1e769lxltnc++Em7l02s87dU2A614+h3VU8tj43/DdU38S9dnCr+fyu02PMb0arjn7SU6fOqPT141ERNYqpaYd8bw+MARnAPcrpS6y9+8BUEo93FEZbQg0ADUNFRSX59EYqCcQCtAUCjAifSxpidlszvuSibnn8sXWBRQVb8Uf9DNj8uXkDj/pCNcsZ/7Sv5BftoP6phqaVBNxDh9Xnf0r/rrscRoCtQxxJ3HLZU+Tkph52GtFUtZQhtfpxRfRI1FKMe/TB1i2/3MaCXDj5Ns5b9qPO3W9ioZy3lr1Jyoa68jNOJ4fnvTvR8w31BlaD/v4Q35Mrs/YbgAACbpJREFUZRLnjEMpxZ6qPdQFqmkK+gkEGzHNAEqZKDOEMkMcl3UymYnDoq5ZUl+CP+hnRNIImswmXIaLkBliY9EqCko2c9a4WQzxtX2WBTUFuAwXQ31WrEdhVT7VDWWMzzqJan8V7615ik2F66gL1HLRhCvZUrgeUSHuuvSZ8DUKawsJmkFGRUzGfrHzXeatfIJNZjlZyk21acUrHOfKwulLYW3VNsb7RrCvsZQRRhqrQwWMxUeeqiVN4vGbfqoMk4yQ8On1X+NwdM9ASr2/lgXrnmHWKXfgckQPeSqleOSjXzK/ZDn/mPkmIzNiG7Hoz4bgCuBipdTP7P1rgdOUUre2Ou/nwM8BRo4ceXJ+fv8KWNJoNEc3pjIxxEApRV1TnbXATDuLEvUk5Q1lDIlLO/KJHdBZQ9Bv3UeVUs8rpaYppaZlZGQcuYBGo9F0I80OFSJCgjuh140A8K2MQFfoC0NQCIyI2M+xj2k0Go2mD+gLQ7AaGCsiY0TEDVwFtO8vp9FoNJoep9fdR5VSQRG5FViI5T46VynVcfiiRqPRaHqUPokjUEp9CHzYF/fWaDQaTTT9drJYo9FoNL2DNgQajUYzyNGGQKPRaAY5vR5QFgsiUgNsB5KBria+Hwl0NXtTV+/TG3LFco/BLFcs99Fy9fxvqeXqvd/SBfiUUkcOxFJK9fs/YI39//kYypbGUKZL9+kNuWK8x6CVq7/+jkeTXDF+Fy1XL/2WzfVmZ/4G2tDQ+zGU6fzitLHfpzfkiuUeg1muWO6j5er5MlqurtPjsg2UoaE1qhP5Mrq7bE+i5eoaWq6uoeXqGv1VLohdtq6UGyg9guf7qGxPouXqGlqurqHl6hr9VS6IXbZOlxsQPQKNRqPR9BwDpUeg0Wg0mh5CGwKNRqMZ5Bw1hkDEXr26HyEis0REicjxfS1LexzpmYnIEhHptQk0EckRkXdFZKeI7BaRJ+wMtR2df4eIxPeSbP1Ov0DrWBdl6bf6Zd+vz3TsqDEE/ZSrgS/t/53GXtd5UCHWeol/B/6plBoLjAMSgAcPU+wOoNde1H6K1rFOoPXr8BxVhkBEEkRkkYh8LSKbROQy+/hoEdkmIi+IyBYR+URE4npaFuBs4KdYay4gItNFZKmILBCR7SLyrIi1DJKI1IrIbBHZAJzRk7K1knO6iHwQsf+UiNzQW/eP4DygUSn1EoBSKgT8GrhRRHwi8mcR2SwiG0XkNhG5HRgGLBaRxb0hYH/Sr2Z50DrWWfq9fkHf6dhRZQiARuD7SqmTgBnAbLslADAWeFopNQErQOMHPSzLZcDHSqkdQJmInGwfPxW4DTgByAUut4/7gFVKqclKqS97WLb+yARgbeQBpVQ1Vmj9z4DRwBSl1InA60qpJ4EDwAyl1IxekrE/6RdoHesKA0G/oI907GgzBAI8JCIbgc+A4cBQ+7M9Sqn19vZarB++J7kaeMvefouWrvtXSqk8u0XyJlaLDiAE/K2HZRqoTAeeU0oFAZRS5X0kR3/SL9A61l1Mp3/oF/SRjvXJwjQ9yDVABnCyUqpJRPYCXvszf8R5IaDHuu4iMgSrKzpJRBTWSmwKWGD/j6R5v9F+cXubINENAm9HJ/YwW4ErIg+ISBJWwq29fSFQO/QL/QKtYzEwEPQL+kjHjrYeQTJQYj/AGcCoPpLjCuBVpdQopdRopdQIYA9wDnCqWOs1G8CPsCb6+pJ84AQR8YhICnB+H8mxCIgXkesgPJk5G5iHtazpL0TEaX82xC5TAyT2ooz9Rb9A61hXGQj6BX2kY0eFIbB/QD/wOjBNRDYB1wHf9JFIVwP/aHXsb/bx1cBTwDasF7f1eb1C8zNTSu0H/gpstv+v6wt5lBXi/n3gShHZCezAGi/9T2AO1ljuRnui88d2seeBj3t6Mq8f6hdoHesS/Vm/oO917KhIMSEik4EXlFKn9rUsh0NEpgO/U0r9az+QZUA8s/7AQHpWWscGJn39rAZ8j0BEfok1IXZvX8syUNDPrPPoZxUb+rl1nv7wrI6KHoFGo9FoYmfA9QhEZK6IlIjI5ohjk0VkhR2A8b7tDYCIuEXkJfv4Brvb3FzmQRHZL/00dYCm7+gOHROReDuo6xs7AOiRPvo6mn5GN9ZhH9vHtogVOBhztPiAMwRYs/wXtzo2B7hbKTUJa2LsTvv4TQD28QuwgjOav/P7WIE3Gk1r5tE9OvZnpdTxwFTgLBG5pKcF1wwI5tE9+vVDpdRkYCKWy+mVsQo04AyBUmop0DrgYxyw1N7+lJaIuxOAz+1yJVjReNPs/ZVKqaIeF1gz4OgOHVNK1SulFtvHA8DXQE4Pi64ZAHRjHVZtn+ME3LSNH+k0A84QdMAWrHB7sKziCHt7A/A9EXGKyBjg5IjPNJquELOO2b7zl2L5sms07RGTfonIQqAEK+Zhfqw3P1oMwY3AzSKyFisAJGAfnwsUAGuAx4HlWBF5Gk1XiUnHbP/wN4EnlVJ5vSqxZiARk34ppS4CsgEPVqR5TBwVKSaUUt8AFwKIyDhgpn08iJVhEPuz5ViBJBpNl/gWOvY8sFMp9XjvSasZaHybOkwp1Sgi72L1KD6N5f5HRY9ARDLt/waWL+6z9n68iPjs7QuAoFJqa58JqhmwxKJjIvI/WCkD7ugToTUDhq7ql1jpqrPt404swxFzFPKA6xGIyJtY2QLTRaQAuA9IEJFb7FP+Drxkb2cCC0XEBAqBayOu80esUPJ4+zpzlFL398qX0PRrukPHRCQH+D3Wy/m1WJmEn1JKzemt76Hpn3RTHeYD3hMRD1aDfjG28YhJJh1QptFoNIObo2JoSKPRaDSxow2BRqPRDHK0IdBoNJpBjjYEGo1GM8jRhkCj0WgGOdoQaDStEJEUEbnZ3h4mIjGH7ms0AwHtPqrRtEJERgMfKKUm9rEoGk2vMOACyjSaXuARIFdE1gM7gfFKqYkicgMwCyuYZyzwZ6ysj9dirTf7L0qpchHJBZ7GSg1cD9xkpxDQaPolemhIo2nL3cBupdQUWvLCNzMRuBw4BXgQqFdKTQVWYC02DlZ+oduUUicDvwOe6RWpNZoY0T0CjaZrLFZK1QA1IlKFtcARwCbgRBFJAM4E3rHTSoCVGVKj6bdoQ6DRdA1/xLYZsW9ivU8GUGn3JjSaAYEeGtJo2lKDlRO+y9irRu0RkSsBxGJydwqn0XQ32hBoNK1QSpUBy+zFxf8UwyWuAX4qIhuIXnlKo+mXaPdRjUajGeToHoFGo9EMcrQh0Gg0mkGONgQajUYzyNGGQKPRaAY52hBoNBrNIEcbAo1GoxnkaEOg0Wg0gxxtCDQajWaQ87/k7j58yRrpDgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "q = p.sensitivity_basin_parameter('ecal', [0.9, 1.0, 1.1])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Or on the commandline, e.g.:\n", "```console\n", "$ swimpy sensitivity_basin_parameter ecal \"[0.9, 1.0, 1.1]\"\n", "```" ] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.14" } }, "nbformat": 4, "nbformat_minor": 2 }