{ "cells": [ { "cell_type": "markdown", "id": "9ce08d3d-75c1-463b-8b25-ec4665fc6dd1", "metadata": {}, "source": [ "## Initial structure" ] }, { "cell_type": "markdown", "id": "0bb5270a-7dbf-4905-8d9c-3a17b0741c62", "metadata": {}, "source": [ "To get the initial diamond structure in AMSinput:\n", "\n", "1. Builders -> Crystal Structure -> Cubic -> Diamond\n", "2. Click \"Ok\"\n", "3. File -> \"Export Coordinates (System)\" -> \".in...\"\n", "4. Copy the contents of the resulting .in file" ] }, { "cell_type": "code", "execution_count": 1, "id": "5886c299-e977-4ffa-aca8-6c8916b1375b", "metadata": {}, "outputs": [], "source": [ "from scm.base import ChemicalSystem\n", "from scm.plams import view\n", "\n", "diamond = ChemicalSystem(\"\"\"\n", "System\n", " Atoms\n", " C -0.44625 -0.44625 -0.44625 \n", " C 0.44625 0.44625 0.44625 \n", " End\n", " Lattice\n", " 0.0 1.785 1.785\n", " 1.785 0.0 1.785\n", " 1.785 1.785 0.0\n", " End\n", "End\n", "\"\"\")" ] }, { "cell_type": "code", "execution_count": 2, "id": "b88507c4-f90a-4901-aadd-80cc736aecb9", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAyAAAAGQCAYAAABWJQQ0AADJjElEQVR4Ae29CYBcVZn+fWrr6n1L0lkIIUDYQlgkBAICCWtAcRtNHEdmVBTCyIwBZlFwSYIfuP4FdXQkjDozjqMSHJ0RkU1JQMQAYQuEHbKRtdP7Xsv9nudUncrtSndS3emlqvu5enPueu65v1Mk73Pf9z3HGC0iIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiIAIiMJoEli9fHly8eHFoNNugZ4uACIiACIiACIiACIiACIiACIiACIiACIiACIiACIiACIjAkBEIsKYzzjjjyE984hPnYDM4ZDWrIhEQAREQgYIioH8ACqq71FgREAERKBQCXmDxYi+0YIEXvvrqpyKe5wWmT59f8t73vveRurq6h//2b//2OL4JQ7IK5Y3UThEQAREQgaEhEB6aalSLCIiACIiACDgCFBWB5OrVJsEja9cas2qVPRfr6rr4ZYiRno6Ojk53tUoREAEREIHxRcC6xMfXK+ttRUAEREAEhpEA/13x5s//U21t7bxTysvD55WUmOZw2Nw/fbp5ZeXKQHLGjJNqtmzZ0DiMbVDVIiACIiACeUxAAiSPO0dNEwEREIHCIkDPx8rkSSd9///V1Jz98erqU2qrqoyJRo1JJk0sHE6+EYsFf/HjHwdWMDxr9eqA9ZAU1juOidZakciwOPc2K1as8G977nggEHDb9h53XKUIiIAIHAqBzF84h1KJ7hUBERABERjvBDi07upEXd2ii6ZM+YsHa2vPMeXls5JlZUXJoiJjIhEThheEQsTs3Rtf/KtfRe6WCBmS34wTBgEIikyFfQkKiAme33dR5ur+N3hPMplEkREi/V+sMyIgAiKQIwEJkBxB6TIREAEREIEDEUgJkMmT3/trrO+BAElUVs6KlJaGDAUIxQfkRyKRCAYaG+Mte/Y8d9bjj899ZflyE2BY1oFqHmfn9hMUfjFBFti3ImIwogBJ/+gNE0R8nNfW1lYdCoWisVjMgzqcGI/HS1An+2IS9qPoLFNfXz/h1FNP3dTe3v7Epz/96UZ6TQbzXLZbiwiIgAg4AkpCdyRUioAIiIAIDJYAjObVMFznw4CNHB8KFQUDgZAJBIJYkY6e+dQVCMEDAhUSK0smO2AI86v6vjCgwT48z++zgmIA4U7OQ0FL372aO2b3V65c6Y5zFLFiiomurq4iPKOGYiISifDGKViD2A9isdsoA/Bm1PI4xAf6JVCCe6BBQqw/ZF1U2KDAiEajPevWrXvHnj17Dsf5p0844YTShx9++Hc41yURAkhaREAEDomABMgh4dPNIiACIiACaQJeRUVTaTA4YypEiAkGw7BVAz7xAamBBeZtMJHobmhtff5N3gdbupdxnec0rZhAGwcS7mTfj0a9790y234xQcP+hhtuKJ45c2YS3ogyjBRWATEBR0SiEiKgGscYCxWFkKhDXUlcT29GDcQExUoIawnFBOuB0AjhOno2+Fhey4XiIo7SKhuUzMFJxVd5Xpx1WBdVINCzfv36E7Zu3Tqzs7MzvGHDhhNRtlx55ZX3sjJcx0KLCIiACAyagATIoNHpRhEQAREQAT+B1tZd8aqqCD65R2Gk0gOSOutK7jEHxPNikWTSQ2q6aUtdMWp/7icoDhLu5ISD3zvBxrvjEFQrMy+TS7gTRMJECgfclPzqV79aNmnSpHKEO1EwUGhEISQYEgXRlghSXKCkkLPJ+7gPm14C19kNXNrDa3GMIiGG0v4bj20uyNEh/2QY1yd4I46FUXRijafvaywuLm5/4403zsa55PTp01956623ToTwCRYVFcVwjRs0IPO+rFOLCIiACAyUgATIQInpehEQAREQgSwCy2HIr/QmTDj+2GCwtCoYjMC4jcCG5dfyXpfiOL0H5tXJk+c0v/QSbeBenoFeFw9ih0+ztrW71y8oXO4Ez6Wf6wxpNsTd4o7Z/SxBMeThTniI9U6kH06xYT0SaA+3u7FyokYnEKAdGNZmG8vj3KaISKSPUVjY+VVQJnFtAwQHeaBI7oQ3JQkRw3p3QlTYUC1c0whPC8VFEO/accsttyw47LDD9iKs6/Vnn332+MmTJ+865ZRT1k+YMGEX+NFmYFu0iIAIiMAhEZAAOSR8ulkEREAERMCYjc56r8RHe3xhL4IRzRCsbDb24zwUQqJl7drz48ilphHdy+D33WHFBPbHbLgThIElBHEAXRCy3gm8rxUhOBfBdjcgUlAgpybQhrKd2zjXiXvqUQYRYtUJr0g98z5wrLu8vLypoaGB9SYhKOgRyXn5+te/flFPT8+ZfFZra2sZ6g8dddRRz1dVVXVgvw31Kf8jZ5q6UARE4EAEJEAOREfnREAEREAEcibAfOdQqBhf5YuQA5IKwaKDI6VEqDM8GNkBk0jQjvUCK1aYIMKUbP30TvA6HOf1qYtTT/Z7J3gkI1iyvBPMhzjg6E742j/q4U4QCdZbgfdsxHvakCmUuyEiuiEmGGrVhHa2IHk8VFpa2opj7Zs2bQp+61vf6kpzSVHJ7c+MeAPfjBwka97OY9xGvUHkd7wbwuWU6urqRnhHStDORuSiPIbnHw5REkEIFoWRFhEQAREYEgISIEOCUZWIgAiIwLgkYA3c009fGnzqKS85adIF80OhEgiICIzVELwbGa1gRQg0B4xdhg2V/yltTGfCeZyYwHELsq/RneAlYN7IJBjHFDCw0wc2uhPuHfVwJ+RY9NBDQW/CQH4xt912G8WZX0T4tzOg01x5jscy4m0xxNldixebwOrViTRr9p0VH//8z/98b11d3Xm7du36LQTQ78G5B6d+hbreCTESAWfGcA2ovQN5N10rAiIw/ghIgIy/Ptcbi4AIiMChEHCGr9/AjQUCq0xt7YUYdakCc35wLTWhUAArrP5k3K7QDAYJ1vCAPF35b//2tWlbtwYmlpQky2HkMl+hEgZvFRpGkdDn6E44b0d3goHsDG4rKHA9Bc4BR3eirc2XhiGN21PhTnimTcbGuREPd2JbXJvoieC+W+iV4DbeiQW3uUGx4N7bHWfZK/HdHkjdk960O6wpAfHBygIrli+nAGKuibnmmmsuAYd3wsNRXFJSUtPd3d0FL8ivMRrXXuSDFOEaej44otZuVphuq78dmedoQwREQARyJSABkispXScCIiAC45wADWZnBH/729+OVlZWBt5+++0Jra2xyMyZlbGbb/7VkR0dW01Ly3OBzs43MPt5HJMQBk083gTRwXzqULClpcccdtj2o2OxmveHw7FKDOoEkRLiF3YatTbMh8/A0mt0JxjAzH2g24NCw47oxMZg6TW6E7sofV07vubTiGdljSh6hTuhHh7fg5Jf9oPIoWgbznAntgttYMH3zBYUvQz6tIeC17ql13l3MIcy4EFsQG0kPztr1v/77JFHvhF48MHvQ7F4P/7Yx4r3nHjieajjRIR4fQPziBwGXpvnzJnzm6VLl9ZfffXVEfApRZvRNckw+kj2Qg7AdYkIiEBuBPQXSm6cdJUIiIAIjGsCCIliInTyH/7hH+qOPvrok+G1mETxAcMUQsQLxWKRrlCo+4L29rcgKMKBaLQO4iNsuqE7wmFqhAC8DykvxKRJpZt7EOSDYz20yWHosm5bUhgQNPZ7je4EI7gTx2yiNa5pwj4FBI3jXqM7cSI+XoaJ9BrwRT85GuFOaH5GMPQhJvh6mfPcGaYlcBeEFcRH4iPTpq06oqTkqnJMSe9dfvmCNRMmfPax009/R8+OHafurK+PzZ49uwUC5MUzzzzzvvPPP78N7QlMnTqV84pUYNuFrTWznfTO9PNOPK1FBERABHIiIAGSEyZdJAIiIALjlwDFB4zOJLwek15++eUbHn/88cNmzZr1IpKTwxh1qTmR8Iowuismsgt2YQQsiAf+0xKEoEhSVMD7Eccxhj1RhCThEenhSFnwPHgYKitIw7Ydx2nwNuErfAsFCbwR/Y3uFEdbMrkjufYK6rTCJh1CZG+jMe3ux7O5yX1uUMRkzqWPozh4uJO9aPT/sOJjCcKuPjJ9+qoFtbVXzauqiu2EGjy8rGzJO1pa5v9x27YfPvnaa9HTjzuOs6gj/WPXbyE+Ol1f4xUoPDjvB8PhODFiE0otIiACIjAkBCRAhgSjKhEBERCBsUnAGaTIB5iM/IC/wHCsk2CsVsOzMA8T1b2FMKx2DN7U2dDQGoBXo4aej5KSI01x8URTWVmHkpMSeignwaCPmsbGrs6dO//wQ4wYuyEeLw594QtX07PhN/YPCtKJCV44CEGReVY/X/Iz5w/akPy8ICM+PgrxcS7ExxlVVfGJ0WikCrkvMajBqmRyxvnr1n06MH/+9wJTpryEzrsHSe52iF28kvVwICekBv1SShGCMgBPlw11y89XVqtEQAQKjYAESKH1mNorAiIgAiNEwC8+YIMuhpeiCF/LQyg5LGt7RUXFYxgqdl15eX3D1772h5pA4PBbJ00639TUnB2oqJhuysrw+Rzfz2m6MuAqFjOBjo7OjhdeOPKlZcsWMSkE6QhL6SWhjWu9JQcQFBQGA/VO8BGFLij4DrkuvcTHeSnPR3wiYuEqw2FTjhXha2Fk3ifeWVEx+ZTXXrvyF4nEyZ/6xje6vHSInRN38G5F0ceczIV5OR0oG9kIlOOJZ67cdZ0IiMAACUiADBCYLhcBERCB8UAgW3zAcKX48DBrdg3yA944/fTTXweHF/7+7/9+e4rHhLbDDjsCYVgRCAnarJzzg+FXKfEBwxahV0xIbw3PmPFq8a5dnFRvOYbqtaMxWaOWIsTemKrwoKM7pS9TkSKQER9XwPNB8XF6yvMRrtonPqggDPRgCJ6QRHksNvOTGzd+5ypjPpENEfkz6Ks4+4WdkkC/W8GYfZ32RUAERGAwBCRABkNN94iACIjAGCaQLT4gHjj/Rjds1krMF7EZScuvdXZ2FkOQdKa+mAfMxInHHo85QDATegQjVXH01gC/ltPetStTCuxBE3i1tHROM+/DYke9GsMoR+rVeomPBRMmXHV6ZWV8AjwY2eKDrijm4XQh52ZXd3d8emnpx9suu+ypwMqV3/MWLw6l82KY3T8V/cO8nBi2vdraWqsOR+qF9BwREIGxTQBfn7SIgAiIgAiIQIpAX+KDRigEQwkm0Xv9hBNOeBlfxhmGFcBXco5eBSVhw5wgPqJMLsd+2Cc8HFmOoktvSKJl7drz45gFXQatQ3NoZUZ8/A0Tzik+4PmYEI2GqyIRF3aVUoFp8YHxiU0z4uEaMGLYT7duNf+6adPRbMISzBPiFgxLzHAtekAwHkC4ftmyZd0psemuUCkCIiACgycgATJ4drpTBERABMYUASc+vva1rx0DYbEYBmeU4gMvWQbB8WfsP4phbpHZQXs24PICrJCIxYLBUKgYx4tguKYESOo6/pla8OEdIVldNtzKHVN5SAR6iY/zID4w2pXN+ejP8+HER2NPT/y5lpbI7+vr7/zHl166YTmGLYP84KhXdkF/F8Pj5XZVioAIiMCQEpAAGVKcqkwEREAECpOAEx8Y7WoOjM+/gNjg7OCc9K8YYuPez3/+84/iS7gdCQnHOAeHywuwAgQJy/NDoVKIjyhm2A6lw656sUgLj/CfeXTNGozTq+VQCGTEx8f9no/eCed9ez4gPp5pbQ2vhfj48datV3scZjmdrJ8OwWKYVh0ax6F4g9jmZI1aREAERGDICCgHZMhQqiIREAERKEwCfvGBN3iXz+Bk7se9N9544wsoAwi9olFqw65wTVZeAGb0CHLIXXo/UgKE+R9uYfgVPSCwZeUBcVAGX/YSH/R8nMawq4OIjyaEXdHz8WxLS3jN3r1WfGCywhDyP1KjBvRuT5whWBCf7O/dPJUeoUz915uT9kRABAZBQF+gBgFNt4iACIjAWCGQLT7wXm60Iys+4Pl4AdfwY5UHL0c0PR1EJi9gwQKOZMUP7eHaYLAEZcQKEB7LWgKpEbFiDVnHtTswAr3Fx8SJV8114qOfnI+2eNw48fFMlvhYAk2Ix/cSFexvhF9V4Dh0JmaPRB7IwJqoq0VABETgwAQkQA7MR2dFQAREYMwSyFF88N8JGql9LnV1G63xCgFyWihkZ0GHJyQIEZKK/nFeEHg/0gIktJ4V1dX1Nnr7rFwHswlkxMcnDj981QKKD4x2VcvRrig+MOEK2Vv4LuEc4qMZKz0fT0N8uLArej76Eh/wdrA/Of+HFSAIweJ+MxsCD0gvocJjWkRABERgMAT0VWMw1HSPCIiACBQ4gVzFB2YLZx4AonRWYiLB2FRsJ+EFCbgwrd27Z9tAKwiQrn0hWCkBsg+Rh9GvONEgk9Dbi/cd19YACPQWHwi7eodPfJT1IT7o+aD4aEiLj0eQ8/Hv27Zd3Z/48LWFfc4M9CKsSYTeNfnOaVMEREAEDpmAPCCHjFAViIAIiEBBEQgsxnwPFBZMOEfL34W1r7CrIK/xvxnDr5gXgGMogjYvYOFCekfml2D+j+nBYDFO8OP5Pg9I6n77VR05JF2d8Xj3Nh6bPVsekBSbnP7sJT4WDlB8rG9uDuciPig02RqE2tWgf0uxiQEFIC3Tgw/k1FJdJAIiIAI5EJAAyQGSLhEBERCBAicQWLBgQZheD7yHt3r16sTNN998EgzOd2M/J/HBe+H1iOIeF4ZjPeg337wyWVHRVIrE8+nM/0gJkNQEhI5ZamZ0JrF3dTQ1vWgFCBwqrh53mcq+CWTExycRdrUQYVe5ej72wvNB8fEoEs5z9HzYFvT09LCfqSS5uOGWEdllhWTfrdRRERABERgAAQmQAcDSpSIgAiJQaATo7UCbvbVr18bp0fjsZz9bhUnlzsSkghdjIsFWnOO/A/6E8/08H2nDk+E4NdhGfjKjc1J5Adxobd2FEZOKerDaBHSX9+FKXsMRsDwvBqO2k/VoyY1AL/HBnI93VFTEazDJYDVyPg4UduXEx9oBig82C78L9JVd/MMt59ZiXSUCIiACORBQDkgOkHSJCIiACBQiAYoPejtQFkFwfDgajX7ohRdeOGvnzp2VCLPZfPLJJ//Pcccdt7Otre13X/rSlzja1X7iw733tGnTEnv27LEhOTiWyQugP2TChOOPxQhYVQjDwrCtEWiUVB60uxclZ0uHMWtenTx5TvNLL3F0V31N9/HpazMjPq6C5+M8iI9TKD6QcF4dDlvxgWz/Xgnnrcj3aMFQu3tjsfhTac/Hf+aW82GfjyRzhmB58HRNRf+EsMawnTXccl9N1TEREAERGBgBCZCB8dLVIiACIlAQBBhyBfERf9/73rcI8fz/D+LjxMbGRrNr1y7b/vb29mOffPLJa6uqqi75+te/vuGuu+4KLVmyxLo2/C+I7+BWLECkVCAHJAqD1IoQbNtJCdPXViIBPQwPCM4xcidbgARRDT0giZa1a8+PYzAlGwrmf462exHYT3ycmhYfHO2Kno8+xQcECD0fT2K0q0fg+fivAYgP/9MjeAYSz+kCyQy37H4H/uu0LQIiIAKDJaAQrMGS030iIAIikKcE6PlgyNVFF110FQzH+6AITty2bVv8jTfeSGDbgyDhGp8yZUoFPCJf4mt873vf49dvrn0uHR0dFdAeUaoRLC4vwF4fiwWDoVAxwq+KmANiR4FlJdkhWImEJiHsE27vgxnxcXV6qN1cxAdHu6rv7o4/yYTzQxAfbAr6uTgdZte7ZdoTAREQgSEiIAEyRCBVjQiIgAjkAwEXdjVv3ryr8BV7FRKKEzt27Ei8/fbb8FAEQ/ySTSfG4YcfHp44cWKisrLysksvvfQHFCzwmvi9Gr1eB6LDfhLn/VgTECRMXrcCBOFc80OhUoiPKDwgqVnQe90M50dqP/xnlmvW2LyTrEu0CwL7iY9TMMkgw6768nwkkFjDsCuKj71p8bG2oWHQng83zwd+H3VoC4fi5cADXeoZERABERhqAhIgQ01U9YmACIjA6BGwOR/HHHPMed3d3au6urqYtxGEAKHw4JdthtYY5HMgb2MC5uRIhCBKKEKWXnzxxUspQihg/M1P5wUYRFzRKA1CiKAqLysvAAPvBqPweKRmQc8OwWL4FZPQYcumhYj/CdpOE8iIj2uOOGLVgkmTrjqZ4gNBUP2JDzfPBz0fT9DzAfHx00GGXWX1AgYVsF4y9GtquGX3O8i6TrsiIAIiMCgCEiCDwqabREAERCAvCdDADyFfYyXEh9fQ0GAgQOjxoNjgRIJm8uTJprq62goRGJfwWiBoKhxOlJWV/euFF154VjppvZcI4ZtCgBSn0z4yeQELFthhfSE8wrVIQmeJdb9beXt6FvRYA3e07Eegl/g4D/N8nIJJBg8qPtCfLuyKM5wPhfjAQARh/FbsLOj43dDDpVzR/bpLB0RABA6VgATIoRLU/SIgAiKQHwT493kSCcSnIsF8YVNTEwVIyHk9EIplamtrTUVFhYF3hMOs2lYzoYMLh16FCPn1WWedVUcRgpMH/fehrm6jrQTC47RQiEPwhvnFHGUq/8PlgED/pAVIaD0fWlfnQrJsE8b7H0HMTB5cgkni/5aeD4iPkyE+MMxuZqhdf8I5w64yng8knNPzQfHxs+3bc5nh/ICs8TNgf3IUAStA6OnCfjNvggck9YPhjhYREAEROEQCB/0H5hDr1+0iIAIiIAIjQ4B/nwfg5TgP+RleS0sLh8q1Xg96PhBmZTAUr8E56/2gR8SJEHpBaHxCgNQhNOs3c+fOLcWXcLY64AxP1DEV99gRsCBqbF7A7t2zbQ4IhEcXZ0FPeUBSAmTfK3t4TioEK5FoL953XFsgEATAZJ/iA0PtlmaNdpXJ+eBQu93dsXVNTXaG86EQH77eYO5HAiv7NjPcsu+8NkVABETgkAlIgBwyQlUgAiIgAqNPAAnkbISH5PIuCA8KESs+6Pmg8EDclOns7MwIEp73ixCEV1GExOEhOaOuru77nLSQSekQIFZkMPwKOoVfwalXdvNhCxcaZHbML8H8H9OZA5KaBX2fB4TXQMPwHs6C3hmPd9tZ0GfPlgdkecrDlDy8qqrm6hkz7mPY1UnO83Fg8eHt7u5O/LmpKcJJBodKfKQFh8GAAjXo31L0mRWb+N30GVPHntUiAiIgAoMlIAEyWHK6TwREQATyhAAnEGQC+S233DLn2muvLYG3YxPERRihVjQibSspPhh65VYnULJECOP9Y5gb5GOXXHLJdazznnvuCbF+eD2iMFIpJrjYvICbb16ZrKhoKkXex3QmoKcECGYcTD3SXuh50CgQIBiCt6Op6UUrQFauHLcCJMAk/+WYo2VlIJCcAfFxSVXV/RdOmrTIhV1VHUR87Onp8R5vago9unfv7b9A2BWFDDwohOz6xnIf7B8QrOxnhmFxccMto081ceRgmeo+ERCB/QlIgOzPREdEQAREoGAIUBzQW0HxAcPx3fhiHfu7v/u775eXl7dBXAQpQpATYr0fTnzQK8LVeUGYJ+IWeDpofNITctsFF1xw8fr162M4V4pranAcVTI9JJUXwI1weBdGTCrqSY2C1f8cIJ4XQ72dRbxnPC7sJ7y3x/yalRB253hezaUQHxdMmjTvpIqKGHM+DiY+6nt6ko83NgafbGy8Gp6P65E7ElqZEh5DIj7YL8wFotDE4h9ueTx2md5ZBERgGAlIgAwjXFUtAiIgAsNJwC8+8Jx3QTx0QWxEMMxu58knn/wpGJFtEB0BHPMwKpYVISyzhQhFhRMhEBlMJA8hFAfejYqfIQxrDgSOFTM0SvGcTF4A/SHB4PHHYgSsKng/MGU6tQvcHT4PCK7HcfzhmVcnT57TnKpj3H1ND1EkwvtRdO3SpQs/+clPXnzM3Ln3X1hXN29+SQkTziOVfXk+0C92no9YzIPnI/lYY2PozxAfP3n77Tuvnjs3wtwR8sV6yIsLtcPvYCo6MYQVXZU93PIhP0YViIAIiIAlYN3oYiECIiACIlBYBLLFB1rfTRHBMXK3b9/+hz/96U/Pwq4twbEfI/yK8zpwIkImcPS58u1xqztnhQaS0icgn+QnH/vYxz6AL+MePSYwTAMQOv68gEp4P1B3EcK9IjibLUCCMGMpQBIta9eeH8dgStYTUFi0B99aCLgwQ9muuuqqeRB1PwagE0siEePNn2+e6uxMJltbw+9ta7MJ557DCmAUhRQfTfG4h6F2k39C2NU6iI//hviwIVxr19IzNeQLfjMcpIDagz+l+mXLlnVjm90+JEJnyBusCkVABAqSgDwgBdltarQIiMB4JtCX+EjziKK899Zbb3121qxZUQiGf8f+7TQmMfpVjKFY9IC41e8J4YhZNHpxra2KXhBsJOAFORVD+v4XNEcc+0xUd3kB1s8Ri+HCEEfAKoJ4cbOge728IIzwQg7IuDNgKT4eeeSR+JQpU+Yh5O1+GPQnokxgTaJvko1FRcHfTppkvjtliungiFdgT/5WfKAvKD72QHzQ87Fu716/+GBfDMsCwVrM52sRAREQgeEkIAEynHRVtwiIgAgMMYGDiY/Pf/7zL/Ca119/vQePjsCYvB5G7S9hWEYgQOIchtcvQJgLwpXeDa7ZIgRGcwLH3/nb3/72A5jEsAXnPdTRjbqtAMFX/fmhUCnERxQeECdA7Cn35mnhEf4zD6xZc/D5RdyNBV5a8QH280pLS++HkKuBJyoBAcfwJgo5BK15phLC7wWMUnYbZqZv5wvD+G+m+GDYVVp8/Lmh4er/3rHDeT6GRXwgBMv2E34ndWgFh+LlwAN2uOUC7wc1XwREIA8JKAQrDztFTRIBERCBvgjAKITdGkh+4xvfmAPR8C5cQyHAxXo+nPhgvkHqsLFeCxiSV8GDcSy8HCdhlvQEvRv+UCzUydAqu/I4t12UFUpOZpjYsWPHIoyI1Xj55Zff+8QTT/DfjrQhzJyRKO5JzYKeHYJFhwo9ILBlx5MHxPLB8MfzIN7o+ahB37FP/KFrNnkDisRUAdBL0ai5ExNFfmrHDtMCzweG2k3+iZ4PiI+fD7P4SP9WXMFwPf7QKJJ28yDECRXleOo/x0KlCIjAMBGQB2SYwKpaERABERhKAmnx4SG86iwYtZeibooP/h3en/jg42k00nhshFH5UZYQLjYpncPy9ucJceFYuN6KEX4NhxGd3LRp0+KHHnqoZtWqVR2zZl0aSZ0P14ZCJbgugrWXfc3TXOws6IlErCG1O+b/tOIDvObBgL+f82rAawQNl+z331squSp4PZ4uLTW/R65NS0dH8tFREB/wnIXRVjsLOtrL340+Uo75n6teUARGh0C/fyGOTnP0VBEQAREQgWwCTnx85StfuQjnLsA+xQeNwx4Yir+g54PX+Dwf/ioY0B+GaNmA8mqsQQiPOPNBnAhxuSAsXTgWQ7FQt60HhjQTz5mgHnrppZc+c+GFF054xzsqrAcEwuMdSECnAOHoWSgpWlIrb0YVVoCghvXcr6sb01/SM+IDwsN6PlAmGzF0Lmamt0n+ZNDXQtJF6MR1mLF+TXPziHs+0MUUqxzGzAoQNIX7zWwrPCDyfhCEFhEQgSEjIAEyZChVkQiIgAgMPQEKCxqHFB/Ynocn0CgswjGKkLu/+MUvbmbOB685wNMpFpgPcjfKlagn4pLSKUKcEKH4cAIEgsUKEAoRekogQJgTkKipqTkaOQ0/53wWqeeGuoPB4nQCekqA7GsHk6pTIViJRHvxvuNjcoseoTg9H058QJAlyY0sN2/ebAWIE3XZBNDLBjkhXnNRUWJdMrl0NcKuONQu5wzJvnYY95n7kcAKCblvuOVhfJ6qFgERGKcEJEDGacfrtUVABPKfAA3BLPHRiv0oVg6Nuhqej10UAf14PrJfkMO2MsRmBcr/gCHcKymd3g8KDecNodHMlQIEz4JXI2A9IDCu45gp/aJLLln0LT43GIwckc4BwSX7PCCph1tRFIjHuzrj8e5tPDZ79pj0gNDzEcsWHxAg9t/YKPI7Xn75ZYM8GsPtfkSIBxUZRAJG8/tfeunfl8NTtWr9+hERH2nBYdC3NWhyKd4FAwogqycc7jOmjv2oRQREQAQOhYAEyKHQ070iIAIiMEwE+hMfOD4Y8eFayXAshEQlPo3yKZRhJKUnnRfEiY9sTwiFiBMhtEuxJiory68/77wL/jkeD2ESQuZ/8DAT2d2jOPeHDeHC87o6mppetAJk5coxJ0Bs2BXmz1ji93xQfKR4WOFm+WFuFoMhjQ2S0+0+mWYtHCWr6PFPfap8JURA1rlh30W/R9Em25FohxtumeJzv4YOe2P0ABEQgTFNQAJkTHevXk4ERKAQCcAIhM3XK+zKej5w/FDEB1HQkKRE6ID4wETaZi/FhV+EZHtBuM+VxjJXfsFH20KRSJGpqSn7WnV1ZKLncQ4Q2q2oEYsruc00Es+L4WRnEffH2OLEx1UQIL+Ax4AeBHiF9okPQiE3nLPepd///vfm7bffthw56Z87j1viCG0L4dqf3PbDHzZi1nR6H0bU8Mdkk2wrF0xJ4iXSwy2PsS7T64iACOQDAf7lqUUEREAERCBPCND4w+LP+Rgq8eHekF/W+Xf/W1j/EuuDEBgcepWJ5jaRHNvWMGbiOdqTCRmisYyv/PYcw61KSoqTs2bFQ3v28Hoa073FB+pmagOMWfPq5Mlzml96KfVuOD4WFr/4WAWBweGNyS/j+SAvLizJkVzpXVqzZo058sgjzdFHH20wPwjFSBzHI8gTuX/Pnj2f5bW4Z8Q8IEgyZ0M9hIZNxXM5T0kM215tbW1aUtrX0B8iIAIiMGQEJECGDKUqEgEREIFDIzAC4sM1kLkFzAd5CEbxNTA4f4AwrDgER5hGslspQDgkL70k3MZ1NieElcDMxp/FwYqKJpx/CQb2hTiWLUCCtKVxLtGydu35cQymRK/7iH7VZ1uHYelLfAAbtce+OVX4XO5zsSBQki23X331VfPaa68xHCteVlYWxihZd+7FbOe81lYyCpzokUF/o3ke9FS4ftmyZfS4sTljoc+IVosIiECeEFAIVp50hJohAiIwvgk4Q8832tVQez6yAVOEcJLBO1CuotGJkJsEh+d1YVfMCcGxTDI6jWkKESZRc5uCIxyuMNXV60wkQocKF55Lbdk9fMdHDshYMmAHLT7A2LIjv3T4VRzhb+G33nrLio902BXpjQovtKsYonRf52lLBERABIaJgATIMIFVtSIgAiKQK4FREB+uaQzzCcHovAbl4yhDECCM/e81NC+9IFxpOHPlNo1pCg0moBcVeQjHugtVUHz0sp3TO+E/84GIPCr0f3MOSXyQmVvBMA7eYfC8E2KOno8QhjZmf/QCSG7DvaxIz/OBttThWRyKl0Mudw33c1W/CIjA+CVQ6P8YjN+e05uLgAiMCQIw9mB/DmvOx4E40di1Bi+M4fdjezfCrUIuKZ2eEIoRCg6uuMZ6Q7hNIcKFqQrBYDVEyLMoH8UR5k6nRo+FvY3ruI4JD8iQiQ+wi2MNg+OdYHw1fgOENiriA8/1L8wF8hhKhnU3T0Cc+PxZ/ku1LQIiIAKDJyABMnh2ulMEREAEDonAKIsP13YavjQyd0MIvQdlB0OvIEI8N0GhEyA4nxEhFCP7RAhHeaqEGPkxjr2Mkra6DeWxs6AnErEGHCjkZUjFBz0fSDq34gNQ8kJ8YD4Z5gRVoD3o1iR/D3xnLSIgAiIwLAQkQIYFqyoVAREQgQMTyBPx4RpJEcIv8k+gXddiOwjxkYAIsUPHOgHCkqM4cWVi+j4RwnCsMERIN87dhHsexL7NFbECBNvr+aC6upEPL+JzD3EZMvEBXjbsiuIDqw27QttG3fNBDxzawXGUrQDBb4D7zeTmwrO4rUUEREAEhoqAvnAMFUnVIwIiIAI5Esgz8eFazbipCL5+/zsSzU/B9nUIv4pBdCBfOmKvoQChd4QLR3OiR4QLy1Q+SBTJ1UEksd9ucx0ikUuslySRaC+2FxbeHxnxgeGHV4ELh9rlglc+8GhXtOH9K8UH2Nqwq3wSH74uYe4H3VZFWJPo6ybfOW2KgAiIwJASkAAZUpyqTAREQAQOTCBPxYdrtBue93oY2XNgNF8EwUGjlGFC1qCmAKHxjfMZEcL91DC9SZRhzGtRjST2b8FLEoQH5CIkM2/dzvtnzy4oD8iwiA/mfABFXoRdsU/c7xECqwb9WEoRgpJzwtg+5zVaREAERGCoCUiADDVR1ScCIiAC/RBwxt4IDrXbT0v6PczQGwoOjoL0lxAZD+NL+EnYZ5hQ0A3Pi+2012OfFwA2KwQJ5wrhhHshiJAar7v79kB3d2tbY+MHXuMtK1cWjAAZF+KD/egWeGWiECEMw2Jfd6Bs5DmU/E1oEQEREIEhJaAckCHFqcpEQAREoG8CBSA+XMNpcDK2am9JSckVECGtDL2Cgeox74MihF4QrjiWyQfhNckktQuNVnziD0UC0WhZsrT0v8snTvyruTjsLV68pBD+zRl34oN9hr6mN4QLZq73OBRzKtaOJ7WIgAiIwBATKIR/DIb4lVWdCIiACIwsARp1/JKcx56PbCDWCwLB8fz8+fN/gpNBCIwkxQeH5fWLECdEeielM0fEihBTVhYprqpq+a958+ZNuOuuu/h1PZU4kv3E/NgfEvEB7xET9DM5H+CVV2FXftRIMrf9gTZPxW80xN8pl9ra2nzuJ/8raFsERKAACUiAFGCnqckiIAKFQ6AAxYfNCwDh5JVXXlm3dOnSp48//vgHYKCGMDJWnALEDc/rxAc9IRQgToTgnW0HwXuCf2OCiYqKqqMmTpz4cxq3CxYsYG5BPhq3QyY+wCoz2lU+iw//f0UcaADdxY4LY6lftmxZN3+7/mu0LQIiIAJDRUACZKhIqh4REAERyCJQiOLD/wqHHXZYyZYtWyZ+5jOf+Z8pU6ZsgMAIQ4D0mind7w2hAGEoFr7+24R11hUKISEEMxNWVVVdtGjRotvWrl0bnzt3Lo39fFqGXHyARV4NtXsw2BBNxew3LSIgAiIwEgQkQEaCsp4hAiIw7ggUuvhgh8FjEcd7eJgPJApvyL8hSfl5eDtC7e3tSXpBOjs5TwjzQJgjEsvkgzhPiOt0fFkPY41XVlZed/HFF398/fr1MXhC8kWEjGvxgRAs666CAKlDf3EoXg5AgJHLtIiACIjA8BGQABk+tqpZBERgnBIodPHh8gIwEmsdwnEC8HoUzZgxo2HatGlL0KX1FBttba0Y5aoIOSEtpr19J4RIhxUgDMtyXhD3RR1ChmImhLqS5eXl37vgggvOoCcEs2+P9r9BvcQH2jeoeT5gsHO+Ext2VWieD99/onGGYGENYN3N4+534LtGmyIgAiIwJARG+y//IXkJVSICIiAC+UKg0MWHnyMEiJ1AEJ4Pb/v27YlNmza9AiHxVxjQCkPydnrl5R/xJkz4sikquhDigxMQUoTs84RQgNA455I2bDnaUmlZWdlvzjrrrLr01/fR+ndoP/GBd0MzBzbJ4FgQHxCCYfRVBboJr5Nk3ke+eKf409EiAiIwBgmM1l/8YxClXkkERGCcE+DwQbBhC2q0q5y6jO9VWloaw8UBGKgPYpTWG0pKjgwWFU3BMLtnmBkz7jTTp/8QXpAY1k7rCXGJ6RQhuN8+J23dJ+AFqaupqfk1WJnFixfT4B3pZGcrPqLR6FUQV3aG8/EqPvh7BX/O/2EFCPqK+83ssLRA5KYWERABERhSAhIgQ4pTlYmACIxTAjSgqT4KaajdfrvKGZ4IJ5oKAZHEe1F4MC/AO++8L4Xh1LittPSEdcFgUcjzGHnUbSorZ5vy8nchJKseAsTOG9LfyFgc6jWBpPSzkJT+g9WrVycgQkby36KM+EDIlRUfFEZY+JoMFbMr4XCbixNQLN0KHmMh7Mq+H/5g7kcCK184iRC6JndCpQiIgAgMBwG5WYeDquoUAREYNwRotGHx7rjjjkhjY+MFsMbfgf1WrCWA0IHzd3/+85/fxXyHlStXpuKRCoQOQrBoZPOLeAgGus0LWLjQJN94Y35JIlFSFgxGEVpVhFflQFeemThxiamv/5mJRDrtPk70Muhh8Ns3R10h1MuRsa5GUvrTECF3MCmdeSH2gr7/SKmB3udSrpXexw60lxEfGHZ2VSQUStjYMDQ0u63c54L+y5Tc5jpWxAfehe/twQtUg7IU+wkeQL+zQ7WIgAiIwLARkAAZNrSqWAREYKwTcAbct7/97SjExwdhmB4B+43iowzntiEh+7cQHU3p6wpKfFAw4X2iaDutVHal/ffi5ptXIpH8+NKKihl4V84dEUbiRBhD7yZNVdWJprT0LCSo/wnHp2AN0L63IsSVtG1ZH0UIFo6M9YMLL7zwjd///vcPZYmQIPatZ2ThwoXJvsQb27hx48bA7t27A3V1dR6EDNVCf5yt+ChF2FUQ4gMNSEAdBsvRGLYHjbLt4ovafZR4de7akttcx4r4sC+W/gPhclGIEIZhkV0HykaeQpkCkL5OhQiIgAgMFQEJkKEiqXpEQATGFYGUXR7wKD4wLO2HYKxNx7EWGNblKLch0frum266yU7mVmiGHNsL4z6KDq3BNpw6CTgvws2ug8PhXRgxqagnFIrCSA1j5RkPc36EkAvyafP88w+baLQLx739BAjqc8dsxjc8EUmw+vm55547Dx6Qt2bPnl0EUcF8kyT2rZhAyQeEzjnnnEo8w0N7AhgGmKIk0yZewMWJEogR3ksDmq3jv3Uxig8PYVdJiI+qcDg4u6wssAHDCadfIHVx6mXGjfgAFw4MwJHLLCv8dpOzZs3q5nEtIiACIjBcBCRAhous6hUBERizBPoSH/gy3uYXH24maRrzhQgCQ+4m9uzZY/M/0P5MXgCdAsHg8ccGgyVV8IDAO8IP5zwWQs5H0kyadJapqDjHtLQ8Crv+MBzvRpkRHXabPMCKa5AGLwTIBIiKuzBB4QLMEdLB8xAbs/FV/gwIjlNw2clYp0IETXF14ZLEzJkzN4J7F849gfJpfMl/AqLkbd7PJe1RIf9YOcRHDOKjNhJJvKemJrho4sTAO6urzT9u2mT+a/duM7moyLaNF6NNvH3Mez6Q60Nx5oHdVHBlmB0TenrgzUsBsBT0hwiIgAgMPQEJkKFnqhpFQATGMIGxLj7c+2HywQoY/1EYp1aEYNufFwBPRBTKowjnIlaAUIRQa7E84ojPmKefXmuH5WVUjxMNLLNXGL0UIQmMinU6kt6/N3/+/PX4In8FHnc6yhA8JNazgsvsvWyGqw/b5/Gnhv1LYDhz9K2Wj3zkI2tQzy+fffbZX8BzYr/klxUVXQfxcRvER/KvJkwIvnvSpMCM0lJThDr/39FHmyo84+76+ozwYJ0UIW4FAoZexfCMSAHP88HX6nMhY3pAsHCAgIalS5fGsI3NwhTPfb6kDoqACOQVAQmQvOoONUYERCCfCTijzB92BcN0THk+HH9MPlgBoz8KI7QHx1xegP1iHotRNhTDg8EcEAoC3sWP5kEY6gYjYp0C470aI2I1O0+HLVFXpuQdTEqn8YtwqhDyODgS08eRmP5xHqPggMBIYEXUFPNM7AR5LDN1YNuGaPHhPI/7KlG+F0LhvRAyn0NY11d/9l//dXYiEllai0kGKT4ug/iYjpCjSjy7BGKmBuu/nHiiOXPrVvOJV14xk/HsHrwE2mJX9C89BEmKD+T03AlD/Wo8j2LMhXhhs7AXvF4xBVyarV9oFvaLqfUiIAJ5S2Akhz7MWwhqmAiIgAgcjMB4Eh9kAbFgP4nzvbEmIEjoTbBSA6FR80OhUhisUXhAnABJeTdgzGJiwhIIi6NMa2sDvCDdWLswP0h3ZuUcIbwOxrx5++23zebNmzlkbxBCI4k1TuEBMcHRmegBCWNliXlHigL+lcd5DdZwuqQQYcIKPSonPPP00/8BrbS0Bp6Pj06cGKL4OAKeD+wb5ICYYogPvICJo20fqaszFyEkawe2+Q8jBQjaSPHBnJMQ2nx7Wnzw9JgQHytWrLChVnjHOnLDO3Ms4ia8nxYREAERGFYC8oAMK15VLgIiMBYI0AjH0ivhHEbbmPR8wCi1Xg6IgDoY3jRIrRFeW1trxUeqPzllOBPQI1hTAiTlBaHhnrSejYqKM8z27fea4mLmgKSS0dNf2O15hHiZvXv3YvSshEEOCJLWo/SGUIRwtZ4RtMF+lee+23Yl60LbMivbhT5hG0MhXP/7Bx5IPvfii960ysoAxAdzPgzFRzXOVUKAMPwKlVuxgUlOrNfjx7NmmZu3bDH/tWuXKcKroD56PkIUHgy9Yt18DNaxliNhhz8GU/JrxmpWpH8H3NYiAiIgAkNNQB6QoSaq+kRABMYUgfEkPvwdB0MfDgIbjQP7P1zPpPoFC5bbfzOQeF4bCpXA+OcIWNkROwEbhjVhwsXwagS8trZmhFh1WS8IPSH0fjQ1NZlt27ZZLwhmRYfHpMiKEgoNt/LZXCFK7MrjvI77/pLbbp8ipgz1rfnDHwzER/Cw6urQFRAfl/YjPuz7UnxABDWgXXD1mFmooxteDyxWfKC9VzP0CtdGsCawjinxgVHDwhBZdhb0tIDTh0n7w9AfIiACw0lAAmQ46apuERCBgiYwXsVHf51WV7fRGt8QHqdhGF4KEJuTQe/HvjWE0CrmgcyFZ2NWAAIkifAt5nnYEKzm5mYrQIqLi+EdKbYiw4mN7NKJEZYUHq7kNsWGEx8srfiAJ+Wh++83Tz/3nIH4MFfU1ppLJ00yTDivxj3M+3CeD/uOEB/dEB/NaHALyhfhlfnlnj0ep32H18Pv+aBRzqGBx9RC7xZeiMOYOQFCp4+dcBIekDEltMZUx+llRGAMEJAAGQOdqFcQAREYegLjVXw4wxMhR1PxZdyOgAWjtIuEd++ezRAdCg8Mfcsk9JQHxIVf8RwdBLDrcU00UVJy+LZEwgu2t3daEYL5UmzIFeedcIIiW3T497MFiBMaTnywpIhhXUheN3967DHzxFNPmek1NVZ8MOxqBp5F8VEFb0oRPTrB9D97PvHRFIuZTWjb/+3e7b3a3p7ENO2hHoRdpT0fFB8HmqHdvnUB/8E+Zv4HioCHd7Z9XcDvo6aLgAgUAAEJkALoJDVRBERgZAk4Y2w8jHbVH1kKAeQE8Cs4iqD9Kr5wIfMf5pdg9KvpqRwQfjxnHkbKA5KqKwBLFu6CWEdHdfXCS7D9TCzWHUT4VZICwgkGv9Bw23gOH2a9IhQfPO5ECEVGtseDdVGEVFRUmFZ4Vh7985/NFAgR6/mg+IDng0PsMuE8cgDx8SbEx73wfPxi795kOzwfiL0a8+KDv3H2F5jWoCjFivkZbX/zsBYREAERGFYCEiDDileVi4AIFBoBiY/UbOLwekTBwoXh2LyAm29emayoaCqFB2Q6E9DhAckIENfPTELnXYlEh/fqq198BdtX4KLN1QiJgmDgkLoZkeEEh790YsSJDwoP/+r3glB8WAGCvI9GhE+VIJTqysmTDT0fh6fFB8OuIvR6+DwfXemwq0Z4PvziA56QUHAciA/XVyyR4xKFp4v9G0Sfd2Byx3oeR5+5vueuFhEQAREYUgISIEOKU5WJgAgUMgGJj4zhiUGgTA2MUNimzLtOjYzEjXB4VxwekB5MRMhQLOv94HF6PdxCAYK155hjbpqAYxunTp3yFJLNgxAVSScwULcVIk58ZO87AcKyP++H9YDAq1JSWWmYxPA3EybYnI/D02FXfeV89CU+Vjc2Big+QhAfHamE87EeduW6yjAcDuwzYqOlpSWznblIGyIgAiIwxAQkQIYYqKoTAREoTAISH/v6bdq0aQl8DedQUJQVTMhu4lkKi2Dw+GODwZIqeEAQakUPSG/xgQF6OUwt8kDir7z++lf3TJlyxGUQEB+EkKD44EzbGeHBbf/qxAhLXNpnCJbzelB8FME7Ug1vR8OmTabooYfMe6dNM4fhOHM+DiQ+mPPxBsKufrdnj7mrocHb3d7eAfGxdDyJjxWpYXY5EtlU9EGIzFF2g6vtP/ahFhEQAREYLgISIMNFVvWKgAgUDAGJj1RXpQWHwRwdFRAAURy1SegIm/KPtVuJU1AeRRAgnAektwBJez9ozGLUKA+jZCW/g4k5MqIjbej2Eh4UIe64Ex+upBChB4ShWy4Ui2FYYRyvhMejCTOY1//Hf5h3IhzrCOSCMN/jQGFXfvGBnI8EwrCCs4855g8QH6vwvuPG8+H+4yRTsEbXe+zjxqVLl8bc78Bdo1IEREAEhpqABMhQE1V9IiACBUWAxhaWcTHJYK4dg2FzK/BlPEowWDqwNuJeG2QFex0LR8Ci4epmQfesEHH1pzwl1W8iH/wDmMh8Fua3S/Au1LOf8PAfwyX2PEuuFB8UHk6EOAESxrkqzFzeiBnUd9x5pzkPw+8ejRyTShw/mPh4Pe35WI1JEPd0dATPO+ssM//ss8/FMhXtZ7yZfU/3LmO9RD8Xp8PsyNwvNMf6q+v9REAERpGABMgowtejRUAERpeAxEff/CEK4mDDBXPzeQkIkm5caQ1zeB/mh0KlMFaj8I44AbLPZof4wC0M1+p+CZOafwo2reeEhRMbfT2V57i4a1m6lQLErfR8lGN+j7cef9y8fccd5lyIj6OgdCooPvA1/0AJ568jUZ1hV6sbGswuCJFFF1wQmHPSSQk8u2rq1KlX4/HeXXfdNS7+XVyRnucDAqQOnDkMLwViE/tAiwiIgAgMN4Fx8RftcENU/SIgAoVHgMY1Fnk+fF0Ho9SqABj7dThMgxSYPK+2tnafwoCPIjUEL8OvUgIkrR1sTdArSKcwpqHh15ciVOscHORkhTb3w/eoXpt4jt33l9zm6kQIS1RiyuDpeOXhh83OH/3InIfRr/oUH1RAmOejk6NdId+jCaNjvYXJEH8H4XE3xQeEyOWXXmpOnTuXM7TjMfY5VyxevLhoyZIl4y0Hws5xAr7shGZ2hPsd2E7RHyIgAiIwDAQkQIYBqqoUARHIbwISHwfuHwiQYnocsCACKly/bNmy7gULltt/L5B4XhsKceSklADZvybEWyVMsq3tT6fh2goa9wdb0B+ZS9z1LN1K8cFaiuHpePkPfzAtd99tLpsxw8xEzgc9H5zro5fnA/U58dEMIfJKU5P5TySqr0LI1m6Ijw9cfrk5/YwzTCdECd6TCRCco+QoTL44F4/xIETGRSjS8uXLwwi/4gBicIQkidgOt4xSiwiIgAgMKwEJkGHFq8pFQATyjYDEx+B6pK5uo1UJEBWnBYNFEAdhejZQwsWRXlM1BwxGzcLEgxurmSeSOkfbNrfFL0Yyd0BERCE2XluzxnT8+tfm0unT7SSDyJS34oM5IXCVpC73eT5asP0CQq7WoHwV95x9zjlm6Sc/aeadeaYTH2ifFToUIBwm+FJWsnv37twbnGlkYW3gvdmfHMbMCRCKEDvh5Ip0eFZhvZFaKwIiUEgE9LWjkHpLbRUBETgkAhIfB8ZHw3PlypUGnoCp9ArAIA/AKO3iXbt3z7ZGOYRHVyoJneFXLgfE1csEkAAmt9uF9W3YtpxOJJ0Ukr6EAiN75Sl3LH2ZLewxiIcw8jzeePRRE7j/frPosMPsJINWfCBBPVt82Hk+EHbViuc8t2uXeQ5zhBz3sY+ZkzByVk9Xl2lD7gfmurDCw/es9LsFzsWx4MKFC5Nr1671nR6zmxzljPkfRRQk3d3dtq/H7NvqxURABPKGgDwgedMVaogIiMBwEoCRBRtLOR+5MGb4Fbwb/EKOImi/ii9caJAbMb8EXo3pqRyQfbOg08vBhbOgc7urawu8IA3YDlthkTq3T4hQWLjFigyfKIHgsffYErFcIUyUtwUJ58EHHkiJD+xTfHC0q4z4YH24j2FXHGa3BfvP7NhhnkPI1jFXX21KMSt6B0KvKD66IULwO3CPdyUmQEfiijEnnX322WUQYcwD2e8id3Ghl/xvge8Ar08NilKsyXR/87AWERABERh2AhIgw45YDxABERhtAhIfufcA8gKCMP6jYOZUgvWU33zzymRFRVMpRMV05n8Eg/sEyL7aKTIoQN5EiWlADP+J6S08KCy4sLQiAze4ko/kas9BEAQwqeCOJ580JfBGXALPx3SIDwoPzvWRyfngA7G6nA8nPp6vqTGzr73WRFFHd2fnvhAtiA/3nKySo3UVQ6zQKOcyZgVI6vUMvFQ9UQgv9i/7vGPmzJn19sVT4VnuMpUiIAIiMOQEJECGHKkqFAERyCcCMDLl+cixQ+ghwqWMm6rBNmxT6xVodreHw7vi8ID0cI5CejecI8GVvI56oLt7C7bsvRlj3y8y/NtOBFjRkRYlmHqdn+fNnmefNdXr1pmLMMP5dAgJio/9PB9+8YH7n4bn43mMlDUH4gNf+E0XEs35UqzT/wz33HTJUDMKkHLMBH4M3wOJ6GNegJRA0KX7nK/M0DSi0iICIiACw05AAmTYEesBIiACo0XAiY877rijtLOz84MwtqajLW1Yq3BuW1lZ2d0c4cldN1rtzKfnTps2LQFjPEkmaBfs9ngT20dhEQwef2wwWFIFDwhUHT0gqQR0nncLr8OUEnYXdWQESF/b2YLA7kP0eBAOTS++aKY884y5aOpUc3g/4oN1ZjwfFB87d5oNFB9/93cmjJGxGHaFl2FSvImjXgoqrnyOK7PbhXNjfhhe5PpYcYVXnYpODEF4UYh0Q3yN+Xd3v1OVIiACo0tASeijy19PFwERGCYCDCWCUZW85ZZbJtfX1/8FDM1yGFrtOFaO9XmIjwclPvbBBx9gCXhtbW1IsQhFaYjzALb9Q9JW4hSURxHORXC2twCh+KD5jrlAaNDaymnPc6Wh77Zd6Yx/7lMQcJ6PJL7Kt7/yijnm9dfN+RAfh0F8VPTh+eC9fvHxDMQHcz5OvOYaO3t6e2urFR9IqGdSvRUhToCw9LfHbbMcTwtnloc4w2t77OPGpUuXxrDNbh9fIMZTp+tdRSBPCEiA5ElHqBkiIAJDR4Dig4nEFB8wphaj5mKUSAQwFVifvOmmmx5KP03GVhqEKzDrOabWCEYBpgfHOlA2oqSa8GIxnAkWwxPCHBA3ApaXERu8LC1AaMDiHoqOlBfECgxoGRr5XJ3owNi3TEAwHr0TzNd4800zZ8sWc/6UKWYaxEgVvs5TgPgTznl/tvh4BsP0HochdllfG0a5wtBOJtbTY8VHfyLECSGWrk24NqWc8AJjfcF7F7Nf0Klc/UJzrL+63k8ERGCUCSgEa5Q7QI8XAREYMgI0HEMLFiwI+8UHDOgojnN4UU6K9+SNN974EIxNXstVX3oBwb+Akf0kTkZYExAk3ThvjXLkVMwPhUphrEbhAXECxJ7KVAG7H8Z8PHPQCQ2/sU+j1+5TiNAbgRXaxiQhPE7autWcP3mymQYxUgUxUYGv9PuJD4RUcYZzTjL4DIbaXY9hemf+zd9wynXT2txsuiE8ONoVkqwzqxMhLPl8t7r2cR/eAIacIQFlbC8r0vN84N3rIDw4DC+9hU1j+631diIgAvlEQAIkn3pDbREBERgMgUB65mqKiQTmb4jDAzIF20soPmBc0YDuJT5wnNdKfACCW2CUWtGAiKs6HKNBCnSeV1tbmxET0CGw8YvxJycYTCWhpyOtbDXQEx49IIFA6XZsw5BPjThFIx9VZbweVnzgmBMfPagkuH27OQ1hVOfX1dmwK4501afnw4kPCIZncf0T8JJM/fCHbWc68dEF8cEV81rYlUKEwoPaIi00MgKE+2gPXzWI82247iW+zOrVq8dDPoQVWxAh7ONmvrf7HXBbiwiIgAgMFwGFYA0XWdUrAiIwEgT4ESUJYzFx3HHHwV4Nz62pqQlXVlYeCWMyDAM0jmN9iY+RaFtBPgMCpJgNh0EOdOF65sksWLA8vHbtSng9wrU8nRIgfUXshGDYG4wmFb0XHpK/RjX23xjnbWA6iRMfPMacjx54OCp27zbzETa1gJ6PaNSOdHVA8QHx8izu+TOunfDe91rxwbArekCc0OFz/IKD21ydEHHnnCDh9TjX3NDQwFC9Mb9ApIfx7gy3w6snw+gb2QNjvtf1giKQPwTkAcmfvlBLREAEBkAg7fVIQnhMmzNnzp0wljfCmHq4ubn5wV/+8pc3IJk6UlFRwbFkn3FhV2nPxwCeoktJoK5uo/UWQYCc5iYh5Edzej/cyuvg5IDAMIFIZPo9OLMr7WSCVySVZE5j360Mu+rGzSV795rzKD7g+chFfDShrmf37DF/goekZNGijPjogpcDI51Zz4fzgPhDsPxeECdG2BZu0wLHQgHyysaNG9uYQ8TX4TuNxSX93wGHMWNOVOrtk0k74SQ8IGP2vcdiX+qdRKBQCUiAFGrPqd0iMI4JUHzQ63HiiSfOw0g+T0F4fArG1HRnSeEr9vGPP/74Z1999dWuL37xi/ddffXVqSGbxjGzg726MzxhhE+FYW5HwAJP5s6Y3btn2zAsCJAuChBMQgjh4XJAXM1MLrcOCFNWdsEepI88iz7hScY3ZYQHjX6KD5uU09hoLoRoOG8A4uM5iI9HUW9gwQL7MI52xZyPLp/4cKFXrnTiw3k/nAhyXhCUSLCP0RK3gxOsWbNmPPzbyD5m/geKgAdWtq9db6oUAREQgeEkMB7+kh1OfqpbBERg5AkEKT5mz559Bh59P4ynqTAgOXwov9wyd4F/r8HGTNTu2rVrxYUXXvjhVatWxfBVm9awtYhRaumHAMOkIOgsS5T2q/jChQb5EPNLMPrV9JQHBMFTFAGgmdIYrCwAS9YEYrHWzokTa5+HQHnYeUB41hn9FB+d8DbUQDhciu1zJ00yUw8QdtUBD0UTxEEj7nmuvt6sQTfHzjzTig/O80GR4c/38G/znF98OAHiSgqQ9BpK1/Mg21pXVzdmvQAUHHxHDChQg6IUazLd3zysRQREQARGhIAEyIhg1kNEQASGiAD/zvJOPvnkOhjAv4HhVEOlge0IjlvDKv2cEEKykpjrY3J1dfXPL7300us4MhbvTYfXpC9T4SdANvACMHHfGeA2L+Dmm1cmKyqaSuEBmU7U8IBkBIi73/OYaG4Q8tYdq6kxSPwP/h/q4nToNlkE23ZY3HaIjonwVrwHauXciRPNVI525eb5gPhB5eilVDY7xQdHu2rCPRQfD2O//dRTrfjobG83PThH4eBWv/joK/yKwoMrxZATH9hHLnwigNG+nsb2y2RAgevea6yW4BPFe7N/2ecdM2fOrOe7omNd34/VV9d7iYAI5AEBCZA86AQ1QQREIDcCCL2iyGBex7/gS30dyjgMpl7Z0LSdccwgEd2KFXzpTUCE3Paud73rh0xUpxBJ54/k9tBxclXa8CzC69ZgG2itDd7sXr+1dVccHpAeTEQIvhwFK3XGldyjbIEQif3xj99CXnnzm9h7DCLRVlYEkVGDJPOjIDI+gMTzcyA+prihdnGOnhdbqRMfeH5GfCBP5PcIs9o7e7ZNzGC4VQxixHk4nABxosO/70QHS4oOt+8ECN4zyWMo/7R+/fqONanwK7+YdQjGVFmC0cPSfW7fq6WlRcJjTPWwXkYE8puABEh+949aJwIikCbg8j4QenU6jNrFtJBhQPUauQf79us2xAfyEMp4p53JG8fjVVVVV86aNev355577lR+4eZ8IYLbm8C0adMS+BrOnA0a4JwTo4lXUFhMmHD8scFgSRXEB0KtGIKV0gs8zwWRb+lha4MvbNzY0mTM8iC8UF9AHd3RSCRQO22aOQI3vR/ig2FXUxB25Yba7U98NNLz0dBgHoSXY9esWfY53Vniwy863DZLJzRcScHB7ewSYiWM5HUPYVdLL7vsshUcxhkPGrOeMuT6WHGFbp6KTgzhvyUKke7i4uJ0/1nM+kMEREAEhpWABMiw4lXlIiACQ0Vg9+7dNJy4fgQGE6OE9vtiy0O0jGFM2jAdtw8DlzFDcQiTeRimd/0FF1xwBg3NuXPnMnRr3C/gZI1SjBxWAVacuNEmoUNA+L1LlTiFOQGLIEBSHhCKELewN1IrJ1DnkL1fTm7duvUFeD4aJh52WGBaLOZd0NFhzq6tTXk+OPxuf54PiIUGrM9DfNwPwbFtxgzr+ejB3B70fPgFhhMdFBf+bSc8/NdSfLjViRGEXpnS0tIARkwLQ6Quv/zyy8eFpwyDN9hcH/Q9+7hx6dKlzKPy9ajrWZUiIAIiMPQEJECGnqlqFAERGAYCEAwUHLCRvFPThtJ+xhKcIgYCw0TxdZ1hOPjKa4UIRQlFCL72JuAZmQpD8yEkp1+JkJtYOhxrv7qG4RXyvkoY4xVgFgUvLh1YG9FoyyaGqcrdJITBoBsBi4Jv32sBNwz8zhYe8bwpE46YPv2BSdOnT50K8XFRZ2eAYVcu58OKD+Z7sAIoF85g2JEOu2qAmNiAEbLuQ57HW1OnWvERQ39SfDhhkS04/ELDnXMiw3k9su/lsL1cJsEjg4XvHINIHReeMvRzMf974YKO9QtNe+wAf9BukO1wAEA6JQIicHAC+kvk4Ix0hQiIwCgTSCeOJ+67775aJpbDeMJXeL/pa5WJDbsqLy+380E449OJEL4CDS3clkS4SQXyQn64aNGiFemE4zEbcjOQrgObOBUeFrqYEhAknEXeSgzk0swPhUrBMArvhhMg+9QHNQQWE4lMXot7imceEX1g4vTDT4f4SEJ8BF3ORyW8HgcUHxAZFB/3YoSrVyBYWGdf4sP1L8WFEx+udOf8gsO/zeu40vsxYcIEtDnlCINI5cZ+nrKxFK6HECzrOcR/F3X474HD8HLkuCa8d98LNCYU4L6OhncMFybT/032fY+OioAIiMBBCEiAHASQTouACIwuARo6TBxHOeW11177C3o3YDzRILIx6zCibNgVSxiQVnxkJyHzSy8NWS64BpcGPSan19bW2pCbs88+e1wnp8MotQYm2NQBEQ1S4PI88PEZnsjywBwgqfArzgOScl5YqPgDyee4KQCj/vXDD5tS8muEXZ02tac7flFHR/AceBgmp0e76k98NEFI2LAriI97MDHhCxUV1itixQf6z4kKfwiVExXOw+FKd9yVfmHC3wb3EW5mBStCr6y3jO+B97aeMpQJCFnrKbvkkkuuZLjeWPGUrTDoawgK/DeUSHiJIKRmKB6Jty7wFoSXvLgkwnKxtzi03Fse5JqWHx6PkdGpp556DIrD06PK+X4fPKtFBERABHIjoCTM3DjpKhEQgVEg4MTHLbfcMhmPXwzjMThv3rxfYKSiG2BcVsJG5lf6EA1P9xXbeTxciWus+GBJgUKhghU2ZiAEYWJDbnDuJCSnvw/ekB0u2X0UXnfUHwk+xWwEeCD9I1y/bNmy7gULlofXrmVOR7g2FOLIScz/6B2xwyF4wTTU0dFodu28YVndEcebKV2d9HyEz6b4SCeclzPnA/ytemG/IOyKw/JytKtG9OHzEB+/aWoyz5WWmgjOxyAUguizELbZb+zTdP9ZsUDBwJXHuHDbLexvLu53wJJClCuH6+V5DlbAbbxrph7WAQ78bVhPGTbpKTsJv43rWV+e/T4oJuyy3CwPbDQbMwB2G5szZc8tNAutWF9pVnorA+hLqIrordGSSCDSHgwEIzW7a9rXBmzyfaoy358nf+Pksj1te7zVgdUdi+9aHJr1zKxPXHPNNVf88Y9/vPEnP/nJf69YscJ+IPDdok0REAEROCgBCZCDItIFIiACo0HALz5gFC7GWgRjMTZjxoz2hQsX/tW//du//SOMxIUQInEMKRpm6JUzNmlc+re57wxSvgsN1rShGcF1NuQGhuZ61LsYhuZjDLnhV+/ReO98e2Zd3UZr4kKAnJbygIThCfFPQpiEAR807W2N5vXXP2pKq3aa6T2lyYu6OoJnI7ypLgfxQc/HcxAfv4X4WIe+iUAkwE1hkmnxwb50wsNfsg/d6rhx3/W163fez5Xiw3lDmCtEb4i7hvdTiOB34H4bnB+DSszDtde9+93vnrZly5Zr8PtoHLbfhy/UafHqxUGzOPVWfjFRZ+q82Wa27ROKCWgJu1BcpLb2/3OtWZs5eMqvTqmeUT8junP7zsN2xHfMwH8Z4baz2lrOe+O894TCoSKsIS/hnREoCkSCYUxI6Zljp0emJ8/87Jl/s7p09dpbvnJLEnOGHP72229fCdY/ffjhh0Npb0jmGdoQAREQgYMRkAA5GCGdFwERGHECfYgPTo7XDeOzor29ff3tt9/+ABp1Hwygf8Wxq5BMzOFjaTAytMQamyydcenfpqeERiYNWRqr2A7jvE1Ox/YahNx8+oEHHrgz/aWbX477NexGHMwwPRBfsT0YkZjFPDYVzDgML11E3SgDS5euCnreXckpU97XRQGCSQjBLZUDQs8HxUd3V5N5/vkrTKTkTXNkZ5G5INkWPGsSTOUcxcezmOeDQ+2+BFHQgwkH0XEmCTHAfqNocP3Vl/jgMS7sS//i73vUgwnYEwF6yujxYBgfhQgX/2/D3eN7Hqq1k6LEMHDBkiOxQIx8FOL0tZw8IU5QrDCBxSvsHDb2mQcUFPYKY1abnCZDDJ3z6DmVrfWtXvGU4qlFNUVTAzFMJFhiqvEOp6DpHoRENf48xYSNFwwFKyDZjulOdHsP/PcDxejvEgYyHnH6EUunHTYNP3R4mvA/ekjc/9wx+KF+/wnvE3/9yEWP/BZex/jLG1++h01ds2aN9a6km61CBERABHIi0Ptv7Jxu0UUiIAIiMHwE+hMfeGIFDKonb7zxxodo/N11110cKpZfqK+D4XgbW8S8DnhCQhhW1SDRPLPS4OSK83alCOHXbidCeC8MVDv/BYyyQGtr6+3333//9Tyek6HJC/Nv4d/vFE9M6Mi0DmIj8/c+hYc7QZZf//rXPwAOR4EPvBzBP9xwww1Pp88Hpk794AtTprx/dnX1/GRl5dHBkhIPfMmtyTz7zF8D/ptmVneROT/RY87CMMh1YM15PvoNu4IYaIAIeAZD7d4HofE6cjEa9uyBJ6XNYBp72zfOI8F+8okCKzZc36Hd/YoPvjfFJRaG6XmYbM+ycL8PTsbn/21w2/023O/D8UEdTNAPQwA3QPC+98EHH7SeMrMCVyx0V3FzX7iTyQRI7Tt/oK25d8yt6jmsJxgsC3qlZaXHhCvDFRAPeKx3SlFFUTU3sJwSLg5Xs7/w5lMD4cBUzDXvQThURSojIJX6HwRDeist0NAY97/Oxk7z1n++hWQfaMsezxz5kSO9isMqMG19wt6DNtrfCK+HHKW/KxkOhkOYR37vj82P63BYouNAHalzIiACByWQ+YfooFfqAhEQAREYZgK5iA8YYLC98JU3ZSTx77AkjNElOPYDbNfAcIxjpKwwJyKkgemEiN/QpBBxhiYNWa5c0l/KWXcQIuQuhtxs2LBh+EJu7FNz/oPvSgM08/f2iiwxAQa8gAZ5RljkWjvYl8LI/wsY2tPx7pW47wnMp/IUtNq0117bUfOb3+z8n8mT319aUzPfq6g4IgCNBzHnmWee/iBCpTaaWT0l5oJEd0Z8uNGuLFt6J9Auf87HXgiQZ+Ht+B3OvI5cjIadO5MQHwHkfKQ6Aw1wfUPxcSABwnfku3Ph+6dXesXoEQsi1CrR2NgYgnjg78GDSA3wt5H9+3AClSUFiBMhrm6KGawheFHiqMt6yuxD+/njnHvOqWmNtXqR2khFSVXJMYEQxETUiwQTwTOCxbDoi2y405nwUoThnYCMMLNRhrHv4X1rwiXwNuF/aMl+YoIywy28Jr3Y/A7sU5zY30D6XOoCyAYwDbZubjVv/QICJBIMhKIhM+sTs+C9iuDx/f9s2HvohWQyljzrrx/762f/8dV/DKxfuj7mHqxSBERABAZCQCFYA6Gla0VABIaNwADFB9tBa4krQ6jugqB4DeUP8aX7HTCg49jmcRtik136jFRrZLKytLELWzMVcoNhepfACBxYyA0rym2hQci25+ydSF+fLS4yFiNCqOyTaSzj/QLwXhQjVp+zmZfBq1OOk8CQZJL5JJRUcVG8M0e9oicpAqN68tNPP31W2ltQilybIoRaHYn9SE1NaXE4zPArjoLF8KtUUvjrr30BGQQvmWO6i634mJ/2fPQrPiA6mhBSxXk+nkHY1X1g8EZlZaBx507T0d6OYDjECTGuK80G7WRYnQ3DYv/w3ShEWLp9bmcW++q2d/Fp3wsl4vymn1jjlXqfa9nZcnI4FF7V3dMdSLYnMWYXRu1CF7B+/2rv3idirFB1z0LJ7Hsmp6Op4VWXXnTp8R0/7nguUhQ52mtHlcWpcCebOxGCxw7hTuXhcrgQgsXhinBJiOzSgoIlRQRL9z/3HjxOGYV9Cgn+zwoLnse274Wxy0MgFoCfAss+VcK9rAXOFN7PO1KrO4+7yIL/c0v2NvaTlGSJSGLB+eef/xRGyApDgLjLVYqACIjAgAhIgAwIly4WAREYDgKDEB/+ZjBZPAwj+xmUF8JI/BkMykUYZhW2cyKTF0LDkoamMzD92zzm/9qNbeRBJ+hJmYc6/3zxxRe/l8np6XAsZyD720CTjtbbQQUFDGZrWKZvphDw15OxAJ2g4Enw4d/VYQyL6+G9GN8fxfsy/GwiBEYJ6qAXiKOCVaNMfPWrXy3D5Hrl+ErPsDLObM7JBSk06BHA7cgzhhDAOTsTXfr+AI7hE3fSzhqPd2/FdrK0tKh9y5amw5LJomLkf6C9GDspEjRbt95uWtr/28yKVSLnI2YOJD5QjZ1k0IoPhF09TfEBofFmVVUQno+mpubmD8HrcCba8zdo/3F8Z2sO4wA3YTgHMGQs5AVCvtBktNeuvC5teOMELkzgSqzxrniHF/N+19nR+Y1W07qO12FZh3nQN6L873gsfnh3W7cXK0fCxARUSlmBbuDjuLK9ruSN7reBtlH40Mj3wpGwh7yQG8I/CpviLxWbCP7HhQa++5/btyco9HAGNdh38p2zPwCeY7hT5ri/koMIi2zZ4RcP6Wc7ngy7Mu072q2og2PFwAuDhiVN3ItD7Vjq9hZXR+oIm2zfzes0nUdzew3+p0UEREAEBktAAmSw5HSfCIjAkBA4RPHh2kARQjOyEcbjpTCwfwADcimSZTNhODQqnWHpDE1nZLrjztCkgYs6XHI6hp8NrYEI+TuIkDsoQvhQbNO4pyGcMShxOFdBUUwxAfEQRD5BLeoJMu8C90/BGoS4oJ1rt1EicilZjedEIT5YfwmuR5NCfC6/xKOw7g2es23iLlfsc+E2JxSkV8BDfcxl4GQobLsVJCg4nHE5BEgRrqPxzUrjuKQbYqOzsbGzIRAoRnRUGe4vNps3fdF0J39lZhRPMxd1N5ozMdTuJF/OB+tG5azIMndD7e6F5+NphF3dD2Zv1VYHG3bsaG2IN767sanxT7CQfz/9g9NvM0+YSxM9iY+auJkLw3omhEUAYT9sFupLF6jXLV7cbnuhMuQ8VAbfDk0Mra27pO6lCX87YU10WjTY9UrXF4NRhDVBOSGn4vj2V9uTXW91Bfbcv8dreLjB7H17rymJlJiyqjJTWVdpjXD3ewAHK0RYOs58N6x4AyiEkkAi+mjUxP8+7pXdWmaQu8FmUWVYEYK7MXYwtumdYH63lRqu5fuXGTGVdcqJAf9h/zG3nV3yensMiJy4wH8FBrkcJhFIGASUmaIJRSZRhNHB8D///dz2/w9MAj3BnkAsEXvK3w5ti4AIiMBgCEiADIaa7hEBERgSAtniA4Yev9bTWM4knOMY7eR9Fmf/T+bXfFqA/JJ/DYzEZ3DfDzjqEYwnZiLbcB4al/7VGZmsltvMDaGRiXtZUmzYkBsc/8Gll176DgiPa3gtF17D9g0k3AnX08ivoZiAiOAzylCP9ZygXRQG9GyweooJLryGM1ZTPFAw8D2tuMChOG7FNBmWEduCS+w1rMCKEBzjp/luXNeJkgw4mlgTt3GuE8f3Yp/XTjr22GP5vgns92CCvjt27tzZXl19fPzRR1+4bMqUS/DVPJ7csuWKUGnxa6Zi0pHmDCSMn11ba2qRvF2Fd2HCOe4lGMLsW3wEA96bZZWm/uVtm9/as/kiPO91rHbZtnob2/grrrM/O3tG25q2ecGK4AfC5eFzvG5vBiUmjPSUgZ8y6k3xMcWm9MTSQMmxyPc5rrgCI0EtgCfho0iatiKg6rgq/Cj2hTpNOnOSCZ6JI38ZDHS2dpqmJ5pM2zNtpvWpVtO8ptl0R7tNRV2F/SWBj30Hlm71i9QgxpUqnlZsut/oNm1fbjNVX6tKvQh/rUBAJ4b9RfJoWnzQqM9e/Mfctit5rdu2ZbpubvuP+7fdPe4aJz48OF/IoqOrA/oOmh11dXV3mbZ4m0nE4BFL4I4kVpYQdunSg1DxGLGG304s0Z14lPUvRLK9f3hfHtMiAiIgArkSSP+VmOvluk4EREAEhoaAEx9f+9rXpkEc/AUMPIqPHhjElVjXcbQrHMPmfuKDf29ZM4xGIZcVK1b0+rtszZo1QQyVGoeBvwjX/BSXTICA4FC7doQsf/Ixk9SZdOxGQHIJyDSkrTGNm2F40fiiIAg1Nzfff8opp/zL1KlTuR1FvcU4Xs7LUPIdeoU7oQ4b7oSSwoGiBZfBtMO7pfe5TcPfnkNh1QdfHNv2kzoezVGprAjhcVzfyTrS1zSi5JiyFBTNuLYFRnIQTJvw/i3wpvCdW3GsvaGhgXXGEd4F67P38s1vfpPhXH/D86iv64wzzvjR+eefb6+bNOnSW6qqTrspGHw8XlG+PVw54UgzDWFYyyBApoIdcz4OJj7Ww/PxQDhgXg+VJ5u6dntFV0Q/N+3KaS3JtuQM5E7EAonAKfBU1MBoZxunhkpDE0IVoTb0QimO1eAYBvzd9z8a99zjaE/W0McFLPk/mM4skzjPH4j9kfA4LwF5jupkd5Dwnbmf2SJ71uwxW+/YalofbzXFRcWmbFKZHcSgKJwaPc39NvwihLUnQ/Aq1MdM8N1BU/K3JalwMP5s0XtsC5fsMnMMpykQMvvpa3l99j3umL/kffh1pq6mePCvFBJuhaBIxjFtCH4+O361w3RtwWSM+EWVHVNmas6pwYye2HHCIwa90Y3ALFwPZ0kQ2waJ8aanoeex57c+fx6faVZagWs39YcIiIAIDJSA/YduoDfpehEQARE4FAJOfGCG8zkwmCkS+EWfBnYJDO2HP/e5z/0Jw+yGICxwKGWc8Xk0nlC4A1QnrhnumNvnF33Gud+PkZwuh+G9CiMhnYS6cMjmOfT6so022C/cruQzaWRyQXu44lEB5kckkJy+6LXXXpuL+v4T00K8BRHC0aMonLjwuZlwJ+zznXC79aiwsVZQ4DjDu5iMweeiCCFKyeZkMJejAadw2GqenayTz8f+HlyDkCg7RG4jhFMPBQXERBfbOZAFdWXAgTHFTxJC5QgIkCI8i88OPv/88yFMuGco5NAe5IPcZ6qr46am6ggTKo2aaXinycjnqAan0oN4Pp6C+HgQguW1YKlp7twdnPGjI0zt/NpvUBSEJqcEBLf5P4oKt2B7gttmybCnXiWuz/waUnQZGMVjqars1b3/wEhUmQN4TxND/oqtAxTqFtaZSQsnmfpH6s32H2437evaTTwUN6UTSzO/F7CwvxX+Pth/6B9Y6EFTNLHIxH4TMx3RDhO+Esn0rBS/Bnod7Dae6i/dNvMvuKSuSv+JdrFtGTHB7bSQYNgUt+nAgEBIiYaYvdaz+ynPBXxjOMa1B2IihqcgjI3beFCg/bV207WjyyDUDUhT18Vb4ybZA+8gVogRQ9GBPBqT6Eo0oZ4WKLUf4zfyr+ZnrME2WX+IgAiIwKAJ6K+RQaPTjSIgAoMhAMMKNlvAo/iAsXsZ6qDxjsMeQ4l+e9NNN73aX7281x/uhLyFCno2sFTCGKyGAU0hkxndCQZTCXItwpi1uQozp//tnj17joPhmIQHJABvCIaS3X++EH7ldiuNTK5oZ6ZJeBYeEw9DgHRMmTLlx8cff/yLaEcV2uY8Egx3siIEZRtubMdqw51gvNbjuiDq7EQd9fTEIESM7WlguyAoOCFgaoa8zBMPvkEuWKzRCkGR+Xsd21aY8RwWbnPDHuMBLk4MInF9PupZgJUhWJsgAv8Hp2kd1x599JH3T5gwYW5FRbVXVlYaLEO7DWaePwWG8T/s3o2wKEgFPAO9aA11f87HkxQf8Hy8miwxLT17zZz/PtFMPG2SiSfiBh6IjCcGd2fajWdy2z+6k/8cTuW2wLTe70L/MbedKWnY40kcEZe5Ei1Pt5jGXzaa9vswfG9ZsYkW7ZtPBr+7zG/D/T5ozMfb4iZxKvwpn4QRPy0tLuiVoIjgim0rIngsLSqsUOA2snPodbD7cGvYEiKCIVDwRFgPBkRBIC0qKBCCFAwUClYs9OCZ3MfcHgiVsiKCQoKCguftNkYk7nqzqzjQHegqDZfuTZQmWuJl8c3JzuTzfGcokhiuewJlHGuia2/XxrbWtg6zw3TsB1MHREAERGCQBAb1l/ogn6XbREAExjkBGGA0lL3bbrvtBBjw74bxbS1EHI9ADDwJwbAR20UwxmthoLvRnSbiHli81hguw/lybqPMhDthmwoBduD+4U6ok/M9xLh++9vfXvLGG28spEEOD0ISIiToRAjFgJsrhGVfIoT3cYEIsTklnE0bddxx9tlnfx2T3MF89IK4jyFQ7Zs2bQp+61vf6uL72pty/4NqzF7dl5jgiXQ7eBEbNND6WUVmcQIEgnABnjsfDKmwtn3pS1/6eWVlZS1G07oP6zy8Z4YXWZUhdK0bc61ciyF0T8Fkgh4MckAx7RhqtxkrE86f2gPxgdCdV2LFJnlYtzn2G8eaiSdPND3xHoMJ9DJtyHUDZvZ+l/qPue3skjfZY7g9l3An8ud1DO+iUf765a+bwOaAKanBbySy77fB3wj62npC/CIk0ZwwGGHLdFzcYTrP6bR1weCnOEtYIcA5L/0eip6UhwJeDSsonKhgmVm7UtvwSBiIBbissHYmOyE0kOQE3RD3dkB47MDRAOrZDCGzCb+TALwcmxHgtomuJSSQt8w0MzdffM7Ff82hhKfUTOmeOX3m3e9Z8p639wObfeAuwFjMJ/XRCdnXal8EREAEDkJg4P8CHKRCnRYBERCBvgjAqKM95H3+858/Gl/6r2f+BBKxq3A8gvkqXjj88MN3wqAvpwGM62xORLoeG4/ktlHaBGyW6WM2ZgjbVoTwGB+U3ueugUGNx3gBPKfxO9/5zjsx38UVuIRioc+8ECdAXF6I+9rtQm5YJ4xJGyMFIRXAe/zLAw88sIzHsxc+1x3LQVC4S0esZPuweF/5ylcWY5thWJHJkyf/eenSpS/MmjXr3okTJ56OifviEGphijWKD67FEGkBGOBTITQW7NplTm1uNl2c5wNrIxL/10GUPICwq1c6MX/IcXFz8k9OMcVVxQZGsGEYlBMJ/hfNPub2+ytpCvsFhf+67O2+9nkM75zxTNiQJ+5THMDpkOjEzOAY2ar+9nrT+q+tpnxauYliPpSiaCpnqC+Rat8Hv0QY/ibRlDCt01vN7rm7TSfGAAjGMSgWhQO8EnaFwOgV7tSZaIbgwEHm+yc2QqzQQ0Rx8Ry2m/CuAa/Lew5FE+RAINmQfA05NK3mcBNo/mNzC661eUZ+ptnbeN/wN772jU/i91vSFe+KdiY6f/a1jV97e27N3OD6V9anFF6dT2Sszvx3ljqXXaH2RUAERGAQBJQDMghoukUERGBwBGjs/vCHP5yJGcYrEA41kbXAyN1wxBFHUHxEID66eQ0OZ0Z34jUwlmD742s05qjgLldc12t0JxjR/YY7QUDsgSckhFnNY08++eQPsX8n7v9fPLOWIVVY+5y0ENfsM1CxzYVfutMrDfcQviR7ECqfefe73308hMhHkTNRP3fu3Mj69alZomnc2xtTf2S2EWrlO5wXm3EwCiOsrPnXv/71cYcddti/w/MxA6IjkS0+KMzC8HiEoPO2o7xz2jRzJI791ZYtpr672zyACQD/XFluNsMTEDnTmOO/M8eEq8IGxi7DrjLiA3Ttix+o9J/L3qb42O+Y67N0iBN+M6lQJxfyRHFxoHAn5FUwv4IiBGKA3qhg4ERk+df0mPYmzJ9RDg8Ek77xHFu3r6Q3hL8NThsYCkMhTAqYyoZKU3RvkXnDeyO+Pbl9TTgZboF34gXcSz2xL9wphHCnXV0bMTN5IlgUTDaub2zO5Vdhc142o73472bhioXhuo113u7ZuwNrV64lXPd7c2LdrPjWikpIyWIKHSTXY1oVvMxqk3jKe4rC3l2fy6N1jQiIgAgMmoAEyKDR6UYREIGBEKBxc8cdd0QwOd6Uo4466pW9e/dW1NTUNMyePfvtpqamMhhvsKEwz1nqKy4Nv0bcY0d3QrkbIqEb19jRnWDk9RrdaSDhTt/4BnKfjfkjvl7PR50/xToPwiEGgzBCgzLbqEQ7Msf4vjQyuVhDEwY4Sjs0LialuwQi6U/nnXfeRx955JEn/SLE3pDHfzDhf+PGjdXTpk1rffDBB4944YUXrodHimFxGQ8RvR8QW5nRwvj+mGQFoyTFTBQej+fBoqGiwtRDkLyOD/EtbzaY8g9GzVG3HA3FmDSIFrIztaQkQ8rO7WvbHXMl5aY1+Ckg3JoSEFZI2HwKOrh8wgLXefQ02ERtf7hTN3Im6IHACo+CP38iZD0SPMewp3S4E49h5CfTvqE9GWzEBCoBTNbHpHK2w/e74LZb+ftwnjIM02tMqTEl8ZLkCT0nBI/sPPKF//vd/13f30/Bion0yQXLF1gxwd3duyEo1q51IsKV9ko82ZZp8ZDJqUFoHYd7DrrJK/G7LobwTtTvqJ+zq23XETjeiD5tRJjdXlaQvt/WpT9EQAREYLgJSIAMN2HVLwIiYL/O0sDZvn37BBhn0xDmY2bMmLEJRvvDCGF6lohwnkPNDnp0J+SV2Oc43AcId6LFFob34zWUi9CeH8A4W8JZwyFGKCa4ZkSHf9sZnS4ky4kQ1BHCvXGEKh2DY/dfdNFF1zz00EN3pfMr+LyUlegal0cl++WKK64o5ohXDz/88FFo999Onz49AsPUzn9C4dGX+CAX5teg/0wnBIjp6DAbMEJTF6z+jtfqTfQDYTPtlsNMu9duBQgD5JyRnhESNOQpJlCmhUTGO+EEBY5nkrFxbUpoMKEaIztBRLjRnYLYdsIiYEUEw5woKHxiwoY7MUkb18Y7483ImdgX7hSnLwAhTzHzBDwUyMaA8GlJvt36fOuKeHd8GrwF8XhFPEyD3/0mWLp38m/zGH4S/E0jZcKK1KBX5HmlRaXXvevd75q/t3Pv+9f9Yd0u/AzoxUN6eWph3W6BByMjJngMvyU7eSXyplCTVwP2SLvpd/JK/oZrcVsQ4prtYD5VBL/P9ueee+7U1tbWSngd34YIb4HQPGjYlmuTShEQAREYKgL8EqhFBERABIaVAAwm2D8BD0bUDBhNH4Zx1sP1ggsu+LfTTz89Y4BlN4L38ZhfTKT3raWGOrnLbW7ss9549OAL80ys8YW23Ibt63gLwouSMNQyyen86u++/DP0iHH/FCDZIoT3wujjzOshiBsDQXP9fffddzsOs21ce325xn6+LMyFSb7vfe/7MgTIP+Jdi+GhSuILeZDvfSDxwUkeuWIUMLt2dXeYxq2Y2G9Gu5nxkxl46YCd4A5vboWDHeHJP7qTP9yJgiI1ZOw+7wRGc+I9zjthR3diEnbaS8EcDYoJCAtMkZdshmjh3BWt0ECvoSc450cTzj9nf3/Im4h1xZ4LhzApugnEcw13QsbLYfGi+K8hTk6HuEwg/yWEhPxUHkz6t+F+H/xtuN8IPSHOG8KOpkCBSLX5Tch/eg3rxx7H8pnPfKYOQzuX4PfE032O5ob228krUfI/JES+BUrS2xQadvLK9I/Jny9FvxBFvX022tj9+uuvz8DK4agDCLHbsXDhwp98/OMfX8t7+d9nug4VIiACIjDsBOQBGXbEeoAIiIAjwORlej9gNDGUKvCb3/yGQ+9mvvRmCwqfUdTLOOojf6LXefe8g5QUHxQGbMv1MNQ24Hm344t+BdpHY5Dzftgv3Cz922hz5su3MzLZdtRBUePRoIdBedu73vWu2ffee++ncSy+ePHiEGZRz7evzfw3IA5P1Nlr1qy5AWIDueXRJFYrwGhMO+OagssZs/R8+MVHZ0en6ejsSDbvak62TGk1tf9Ua9pebgsl2lKza1tvBcUEQ5/iXijjlXDhTjhOMWEFBXMv0qM7QUBAwZjXkNDdinb2PboTcuG9Dm9HYmdih03G3tTcZbYZhvIdcKE4cosLd2KoEzwGAeTvZPqp03S+bXrMObj2dqzXpH8bnFHeesr4W8Bvxf6u+dumR4jeD4oR7nPhPo9jn793jjJ2DETbQ8gb+glmn38K+VAT6NEgX9QVxPV28kpey/v53wG3cd5u4Nk2NBHHeC6G0v5bjm0udgAHXBPG9fZ+1Nf51ltvTd+6deux+G+wHX3aQUGC0MUtuJ5174PBB2oRAREQgWEmIAEyzIBVvQiIQMaDwVGjptJAwkqjiV/Z+UW3L/HQ17HhQMnn0EhjEvqPICZeQLt+DWNwKsJUOAFfJjkdxzNihNs4lxEhbBiNx/SKVwrYkCyEMV11+eWXz66vr38vxEcDwpzCiOXPCC7eN4qLFR9o4zwYy7/Bu5c68QFBYr/wOwFCkcV34ztni4+09yPZ2tQabIg0BEPvDJmGRxtM997u1LwVEBQYgrYNIiIG0kz53wgRkgl3whZntjjQ6E6NuTCygmKzNdbNecvP25eMvWatMWsz3qeMFyrXcCe83wSMBNaD/Jibb7zxxkmYhPKD8G6Rg9fZ2YmBwIqswPB7xCg4uJKZK93vA78dihBOLlmKIaGX3n333UdAqP4fREkEoi7M63Eep/efvBLH3VwzDEazIitdVwP6BpupySvRFptQjn7bCeHD39vZ8LJUIS/pj9gOQ4TEkX91L0Y6exP38Pc6Uv+95dKVukYERGAcEJAAGQedrFcUgXwhQCMNhhYNJdhG4fply5bZUa/ywACikRZG255AORdG4P/iS7RNTkeZSU6nAc4V7c+U3ObKd+OSNhxpeFK8JCoqKt6J7T+fe+65V0B8PJEnyelWfMAQnYeRru6HCKlhDgyMaRt2xRAjej5oXOcgPuKtba3hPXv3PNM6q/XXgfWBYKA90HLWrLMmhYMY+wpLVUnV//xizS/WUdhs/+P2AQmKxXctDu1+cXfgYKM7ZQQFTGl//gT6ZkCTV+J6G+4ET4gNdwKHEgySEMBqMDfK7773ve91rlu37kMw7Ivx20hAiITICZuZ0Dz7Q8Af+F1nRKr73fAYfyZYPfwuPEySeemvfvWrmgsvvPA/0R+dFAxA1sJr0JZO3FeP0k5eiWcgxz9Cz0s3wgSbMJw1K0seaPJKDK98Hq6ZieubURbh+W0ofw3xsZ1ssEh8AIgWERCBkSUgATKyvPU0ERjXBGA40WizRnoegqAIYfjKDrTxQhiHt6O9V+ILOPM6MiE32O5TgMCYswYo7rNfvWloYpueEI4kxeT0h2BkXvf73//+R+lwLH6Nz8X4CyB3JrBmzRoaraaurs7DyGH73ec/D2/Lgeq24oOeD3zZvx/D7tYgZyUJo5e5IFZ0wBC2uQw4tp/nA1/9Dde05yMOT1G4fnf9k29temuR2WQa2cYQooAu/sDFV7K/wSESKg69eud372xsDbQaCgoM+2pHdsoKd8p4J1iHExSrl/QOWwOLXqM7AXGU3gi0dSIE5H6TV2KG9zKIrHIKLPRRFP0QRbv4LIbe5RTuhHsse+RP1F511VWPcz6Z//3f/70KHIohGOiNoFcjIzbYfvY/6rc8yZVVMJ+G21xQBigm0Gb+7s7E8NAUBx9EWOJb9oKB/cE4KnvHihUrAieeeGJgyZIliVtvvfUStGEe6m3BWowLunDdLzEXzy5wxOT1gV7MB/ZIXS0CIiACgydgP8UM/nbdKQIiIAIHJwCjB7ZOarI7GEQzuYNl40033XSPO3fwWkbsChqb1jCDcbsC28v5ZIYnQUjY3Aga6PQQcGWYEo1KlvSCcHWGO96Rt9JbkklOh9G98v7771/B4zQC8fW6TyOQIiU9/CoN1AEvvJ83ZYmRjPig5wNrDbwSSYiJIISEbfuECRMMDNhM6BBtdSbVw9C2wsOJDxj0Caxs41MYUnkR7m/wTsPcLOtN4mPLP1Z0XPS4jyPYLoomhLs6un6C99zd10uAwaBGd2Jd4EvBQePfeit89fuTscl3wJNXokob7sQSq8spSeL9mzFyVMdTTz1V9aMf/eg68DkN4iKO34SdqJG/CZekzt+FS+Ln78L9RvC7sgKF7UXd9vcBHRJCP7wNph+ASH0S/VcEodmr7yEuUioD9+G9uc0fWOYYtlkfTtn/1i7C9jwcasV+FNv0Nq524qO/3x3r0CICIiACw01AAmS4Cat+ERABv1H0ARhBR8MgYh7Iw5/73OeeOJARPoro+HcjhQgTf69EW2/HdgWMyMy8GNkixAkRGp1OhPgNTRiqNkAf7x9qaWn5EeL/r3vllVda096QTNJzep9GszMsw5hb5BiImjNR99Go5wzUH0Hd1tLkdWgfv7A/j+1NOP84jNkX7rnnng7s2yVdJ98pTs8HxQc8KTUwjpMQIEGGF1GAUDhRgMyZMyfzNZ85HxQfTDqn+ICBTO9HAveEduzY8ShCiC5Dve18EEbRCsMTg2rCk3DNx5jTgMMUbY8gFK2N16B9lWhvFTfBwoY7obQCAsdzHt2JdWGxwgL3pTKpU2LE/ruG5wCRnbySI0FRRHBhn3bj8m6U9ACwTWz7fuFOYGjDnSCYM+FOCHnyvvvd7/Jeu4BVBTjcjh3+RhIQH0GEOgX423C/D/4usn8b7veBdlkhgnaRt/WkgFkrROp1DzzwwEA9Zf7/ziQ+Ul2kP0VABPKUAP9x0CICIiACw04AQiMMY64CRhdswySTbfP57x8a/xQFzP/ITk7Ha6RGyMK5jKGeNiL326cRTEMTK2xUm5weQ0LwlZgB/iSIgPfBQ7EjnZyeBCODL9NWjGB44rNhuH8Axvy7sXL0ohANV1cfS9Rn62bnof4LWNJbgXZt+cu//Mvfo50/Q/1/wGrr9IsPfKW3Hh0ayo2NjdaLg+dYw5nPYT0UHU54YNjY1FC7XV1JJNUztGzLX/3VX63EO5yHaypxb/yxxx6rQb2hnTt3ToC4OoNM0MYY3gUThRg7LCyOka31+uAcd3uN7oT34A+E78PrrRcH13HZb3Qn3ozr2tEWPoeVNaKwI0ShtJNX4jyP70HJNgTxbm1oa/tAJq/Eff6F9fF5QXh+WnHik6hvA47dRoGGtnCxI6ihZF/Y34Qr/RXhPtt/qIv1kSlHION/Iz9ctGjRDPTbCl6fLVL9dbhtPB9VyPPheKgUARHIbwLygOR3/6h1IjBmCDDUBl/wPwkjiUZ9FMbXPfCAvOgMpzx+UQolhsJMhWG4Gu1/J9Y4xEHYhdpkf+12X7xpyPtFA9+Rxiven8Z1GN6EHfja/f4//OEPTH63C5LUL8b9/4Sv6RezXp/gSIAZE5ft3A5+AYL28F7nMaF3yV5DEYEv+Ruw+1XMSbIJX+fvmTp1qvV8oO1B1g/PhNm2bZuhF4SeDogU6wWh8EDbrBBhPTSuUdqwKzxr68c+9rHvIH8kgWuKwMV5cJhLEoNXhKNFnZI+7s2bN+9x9L39wm8bmhrlKfPvD1j0Cndi+3GsG/c3gRWvQ5HcCZY2nAr7O5n3gcuSeIcGPDPJhGyIN4qMAS14TqYdzJ9wN+cS7oRreb1tH9q6BO35AfZr0J44WNuQLPfbcOF67rdxME8Z6gnBK/UjcOzTU4bzmYXvQGRIOJfnI0NFGyIgAvlMIJ+/QOYzN7VNBERg4AQYckNDlWE3SRjGTQOvYlTuoPjg1/cdMIIXwvD/Pt7jKoYhYd8mp2PffunGfsYD4t/2ixAYihQVdoQsiICpMFYfwtfua5HgvBEhPV+FoXoRjVVcQ7GRwMqv9gGUXKwgYQmD1//1nF/Q/XDImqNahVDXSfiy/tNzzjnHejporMM4tqNdUXxg38CjYeulUYzwMPsOFB38gs+SYVjorwRESAjXvw3x8T0ksCdwbRH2bX4L24Tr0dRIEqKkmttoE9vQhrZ2Y78Z17SDCwVSE/Zb0MYDje4Uh6DolQPhf8H+tlGnBXEAMcFbKdZ4HS13J9zccZb0RNnS94f/OneYx7hSVN8FFm/h+T8Fq2MgHvCK+zxl/D2434QrcW3m9+LrUzQp5SmDGLwS84ScNHny5I/CG/JaX54Qvi/fQeLDdYlKERCBQiAgAVIIvaQ2ikABE3AGEgzRGhhKpdi3X7thcNmQmgJ5NQonmxMCo/JqGJob8R6ZkBsYlJlJC3F8PzHCYxQhuC8jGmCU2/kgGHKDuv8Tngl6OGjgM3nBw/UMubLzQvA+rkTGleLDlWCaqZPbfBbbipIrRYENd4IRG0TolwdxgNuD1tMBQ9ls2bKF19gcEHo/KEKc6KDw4ErPBzwibO+bCLv6DsKuOuBx6EDbbVgTnmvDnVCyT9/G/C5zUBcnvOuEB+BlPOf/cLznC1/4gm0LG5jLgudlVNUgBEVGMPQhJvj4zPlc2nKQa2I4z2Gcn0Q5H33zf2D6TrDmXDL8bWQmLWSfOAHiSh7jyj5m37AfUQdFDT0pHMXqzxdffDHnknnMP5cM+WCR+DhI5+i0CIhA/hGQAMm/PlGLRGBMEoARG4VxG4bBRCO0A2UjX5QGVIG8MNtNgzgIQ/N2GIrb0fYf4L1q4AmIw+DMTFqYbVj6952QwL00MulB8RiWg3AmhtxQnNl5JSgwuDrx4e5zx1nSWHUr6/MvztDFfTanAm32IAhsaBbbw3Arhl5xm4nnvJ7tgMFrq6EI4TV4rzi8PQwXexKT170Xw9ruxRwSkVWrVmWS3P3P5fY3v/nNaaivHveyb5MIteti+9IGs33WAQQF7+HLDMQ7gcuHVFCwvoEuzlPWgPdeiP6xnjJww27CesrI2q2uf1j6V/azT4Q4T1kt6ltzySWXfBrJ6XdShEAEkk9Cno+BdpOuFwERyAcCEiD50AtqgwiMAwKMhacRjFelcZmAUZsZTaiAXp/t57pfyA2+dttQJGdgurI/Q9MnIGBH2iRkDyFRIeZecN8JDpZ9bfN+Xx32Ht7Hhc/k4toAA5jXBujNYF3M90CSeGYoYT7z6KOPpgiyiebM90D/sM44xFUY1z+J/UUY9aox/Qx+8beCgqUTE5x/4sUXX/QgXKajDfQYwZYOtqAMYBI/3JqadyJdR6qRONmHhyJzDqcLaenTU0a+6AsKEXqRMv3C/uG+O+a26S1zIgQlvUpJiEd0YWgVwvVmYxjn6wnly1/+8iW4Zy42W8FUQ+0SihYREIGCICABUhDdpEaKQOESSBunHIJ2Kowkxv/HsO0hTCdlLRfmq/UXckMjM/O12xmWztB0pTM0KQZ8hmaA4U80ziHU7HEnPCg0uM2VxqnbZ8nV1cF7ubj6WaI9dmWdvJfiY/PmzfaetIfDIDyLExzaMCwKDwoQfLm3kwxino8nEaa1CNU2YqUxTCObD/J7KKxgwPM4AR5DgmLY5nm+TxOvTd+DYswvvTxl6J+X8MY/o6cM/c+QLOspY9/0t5KQ+22wb7FSzHGyxST+u7nusssum3TSSSf9FAMhzEFfteMZEh+EpkUERKBgCEiAFExXqaEiUNgEaPzSA4IFdlS4ftmyZZwYDfZpwYRgZXdAdsjNv+B9lsKApyeEMf829Allv4Ymrs8YmjTWKSZY0t7ktltpjHIlQ3/pzrPkPbyXC+vlc91KzwfFDYfbxfwj9hqGW1FowIi1AoTCw4kPlH2JDxrBFB9c+vVQFNhwy6m3Gfo/yYcrh56+H/12Ifrih/j9vwNhdhkR4vrHX7q+c78NNo19ixXdazs4Do/KRzFy2kUQIt9EyFwI9XLoYU0ySFhaREAECoIA/0HRIgIiIALDTgBGVjG/xo+xhS9kDXO82zUwGq/BPsOogghd4hC11qj3G/cMx+FK74NL8qYBSoOTC41Nioy00ZkRIX7xwW0KCooIhOZkZmR32/6S13K4YD7z1Vdftc/hvfSIIPHdhl7BKLYTDKZDr2zY1a5du/yeD77jQRPIYR9boxtlBa9PL3YG9BW+Wbz5nuNkoUgNo5+fwe/jQvTx/WASZrgefh8evEyZ3wfZ+38b7vfBfnL/3VBoYo6VMAYAiMGjMvnRRx/9CARIC7YlPsbJD0qvKQJjhYA8IGOlJ/UeIpCnBGh4MsYfhlcdmsjhYekZGPB8DXn6emyWC7nhaEd3wEjchGP7hdzgXC9PiNtnycUvOrLFh/N0UEw48eG2eZ87n/pAvs8DQuOVYoTD7L70EiOBMAYyxAeN3WnTppnDDz/cGr0UMhRJMIw5hG6Y4mPr1q0u7Con8WErT/3B2ck550cRBQkM6bHU177XzHnTecoa0deXoq9+ADYZTxnEhfWU4Vgvz5V/n33N3wTmVzEQH+xDTlrooR9LMav9r6699toGeJ6C+O/soCIx51brQhEQAREYRgISIMMIV1WLgAj0IhCHEcWYK8aSuK/iNpeg11WFucMv//SGHDDkhkYlBYdbnZHpSooJLjQ2KSpcyeNu9QsQignuu2t5PRfWx6/mDK/as2eP2bBhgz3O6yk0pk+fbmbOnGmFCO9lHUiAT8LADSLR/EmUAxYfeCa6NeDhGTV4WClWznCecuvYp4/rP5ynLIB+uQZcngGrH9Djgd/CfoMXuN+DK3GtgSi0+TvsQ9zjYVCHAMTHDuTcNHB+EIgPPkOLCIiACBQEAQmQgugmNVIECpvAOMoLyITcoMcuhKH5MxiLi+BVsILACQ9XOgOT+26bJY1MGp0UFFwpEvwCxIkQXsdtnnMihL8U1sdjSCA3zzzzjP3x8FqKjyOOOMIceeSRGfHB6/g13Rm069atuxTPbkQ7XMK5vT/XPxAOxFnuw2g3Bx7ogNCp572oc7yLkX49ZeDfaxhn/gbcb4TsKCL5G6KnimF76GsPfRhA3z9PtLNnz04l//BiLSIgAiJQAARSn8sKoKFqogiIQGESSBuenP/D5gXAuKIh2sy3WTE28wJ6hdzgNVfBmORcGknE/VtDkkKAYVBu5ZdwxvyzJB4nJpwA6U+E+IUIBQYNVB7jDOc0WJ966ilitmFYDMei8JgzZ44VN3yGew7KAPsFo2HVfvjDHz4Rm0F8Vbf3DvQPDrfsFxvMdRhoHWP4erLIeMrQrxdin/khYfQXE/9tLg5zQ/jb4G9i+/btZu/evRYJxB1Lb8aMGZxBvhu/ozu4j1WhVySjRQREoGAIyANSMF2lhopAQRNg7gcNryKsSRjDTQX9NgdvPN+VH3gYcrMUSeDFEBd/A6Nyv3AbfukGG+shqaqqsgKC1cOIt6sTH66kx4LCgULDv/IY76EI4Rfzxx9/3H5F5z7FyKxZs8xxxx1nt3lt9or6k6gvijquxeMfxUoPSM4LxKQNp8P72OGWUR9He+pGDoqM4/0p9vKUoS9+AFZLIDxQ2BHUKAjtsMgUjhSXFCMsMVxyHL+nCPr053ffffdLCr/aH66OiIAI5D8BeUDyv4/UQhEoWAIwomxoCAynGhikzAtggnIARvSAjNsCBUDD24PhH4WwmIt4fRr9AX7ZpiBwIyCx5JdvCguOVuUWigmKDq7+be77xQPvc6KE9zc3N5vHHnvMig8moPOr+bHHHmtOO+00W7W7lnW4utP1M4meYyS/+13veteU1atXU0QNOLSHogj1outtCFcjZk3nnCADrsdxGMNlxlMGkfphvOft4BTE7yOA34fHZHP2JQUI+xB9Y6qrqxPgG8FcLrvh/fgcE8/RT/IwjeEfiV5NBMYqAQmQsdqzei8RyCMCzAuAccUwLC4dWBvZPJRj2XiiyErCiJyLdz8RoUlJ5F8EMWyqDbWigcmVX7Zramrsiutsr4GLLd0f3OdKweBfnRDhOSacc56PtWvXWqOVng/WTa/HvHnzrJeF9/YhPGzdqIMPTeK+ctR1GZ+9YMGCAQtFGMqZ4ZbxvAHf7955nJRO5HFUq+vR/0shODogSilCksz3YB+iZI5ODF0UwvFGiI/33HvvvTvTjORhGic/Fr2mCIwlAhIgY6k39S4ikKcEmBcA44oLQ0sSMKC687SpQ9ksqyJghC+i4Q/DnMa9Oeqoo2wiOIfAZUL4McccYzC7tRUI7uFg5Db7LKkV3MpryZd5Ag8//LAVHwzVYT4JxccZZ5xhw7t4PdvhSrft6nFl+vj7+GCE+xy4Ib7WrUjn8+A161AHh+HlaGdNvku02TcBMqaIiIDdqosvvvhGODm6IESC6MMEBEgSopETWzLsagcE7SJMQvhEOvRK4qNvpjoqAiKQ5wSUA5LnHaTmiUAhE4BRSiOcoyHZvAAYpDFsezC4e3/iL+SX7L/tNA5hiwfP5SV4dxqRNm+Dc3BQJKTDlaxAoOHPxS8+/Nv2ZNYfrI+ihjkfjzzyiK2H9ULgWfFx+umn223exvrRhF41ZO/zMtaJ43NhCJchvKcdx2wf9rrxwDsMLeLzeF8zt93vgNta9ifAfgbzGIbSXQTRF8BgAV8F+08g1OoIhmChP1vB826sX8D8LNspPnCe3hMtIiACIlCQBOQBKchuU6NFoLAIpA1tfulFikG4ftmyZd0wunpbw4X1SgdrLd8tieFRmfcyJy0kmBhhBQZDaxjXz5IGPxd3zpU8nr3tjrFEyI4VMBxql54P1sX8DuaUMOeDng+XO+Cv3z4s/Qfrz1ookpgHMhXCZjbPwdjN+d8JN9wybkM1STLQRy5C7H+hRxDaI+B95StfuQj/nZzG+T4wUlnjdddd982JEydeBo4XwusxG3OzXEnxgaqY9yHx0T9TnREBESgAAjn/w1IA76ImioAI5CkBGFGZvIA8beJQN8uKK+R4HI6KK2Bk0tK3x/xGP7e5gk+fpf8cBQevY8mv4szloOeDOR/ME+A+xQdDupjzwWO81tXBe91z3HP9pQOAYx6+tKO6UK07lktJIxrXZYZbxrO47Oa98IDsp3RyqXOMX8PfA9WHFR/APg/7rRB+YQgODwMXrEb/3oc+/AOOb8PKfBreo7ArQNAiAiJQ2AQkQAq7/9R6EchrAs7whCFah4ZyKF6G+HTldaOHpnFWbOBdp8PALEaVtOutGDhYiXv2Ew1+8UFPBwSCFR+PPvqo9XJQfCA/wBx99NFm7ty5dg4JXufECgULt/2re052e9hPaDM9KqcTBTws9l24ncPCe5n/gSLgQQSNh77OAUvvS8gHR7xvf/vb0VtvvfVybM/F2gpmJTjXjf69Cx6QXRgEgB4kCg/+W02vh4QcIGgRAREofAJyjxd+H+oNRKAQCMRhVNEqZWKy+ypujbBCaPxg24h3jbl7sw19/z7FgFt53G37BQO3KSpczseGDRusUOFQuxQfDLs69dRTbc6Hq5v3UHzwPidI/HXyOdx3z3T3sc1oO+dsyWnBfbg84HG4ZdzAsLME+zunm8fZRY4VxQc8Vh/C/hFg14y1HNvbINp+i1yQpvR1Np9mnCHS64qACIwDAhIg46CT9YoiMJoEXF4ADFLYu8kwvtaPp793Mt4DJyr6KmFsZkQABQE9Gn6hwG3mczDBnKNdbdy40XapEx+cZJAznDMEiwuMWVuf83w4AeJKf92uPWyD22aJOgYc6sPhlpnkQ/GBOjpmzpxZn26PxAhApEWF9XxQfIDxdBxrxikrPjCPy9033XSTzY/COTHjj0eLCIjAmCSgEKwx2a16KRHIDwJpIyqTFwBji0YVDa7xkheQMSL56jTs/aXb5nGKApb+bScUKBy4cN6QZ555xooRsLWCY/r06eaEE06w2xx6l7kfLLlStLjViQ+KEidM3DNdyfa4NuHYgP994HDAfsO5paUl8/72BcbxH+AKNIFe4gN93QaxlhEfbnAGP8NxjEyvLgIiMIYJjKcvkWO4G/VqIpDXBJj7wfh1hvQkYfw25XVrh7BxeO/M37FOXLiSRj89HW7fGf7c5zmKBBinNnSKo4hhSFbzxhtvWO8Gw7Do7YCHwSadc5vXsj4Yr/YNXH2uLidA/KU7557Jfd92T64oVmQNt4y2xFFPNzw0A/ai5PrMQroOfSHxUUgdpraKgAgMO4HMP47D/iQ9QAREYFwRcEYX8wJgfZViP0ErDEbyeJgd2375hyHejPem+8L+XeuMeyLAcSs0KBxo+LP0bztxwKF1KT42b95sfz9ung9OZIj5IqwQYX1cXR28kPW757m6/OIj2xPixAeOB3kO++tZz0AmI6RQwr14tMc+bly6dGkM2+z2cesJce/vcj7AYjr6RZ4P/ri0iIAIjFsCEiDjtuv14iIwMgSYFwCjmWFY/BregbKRTx7jRqk1uOGpeA0hUe1412q8sgej3k5G6IQBBQO3/fu41goJMqL4QBiT2bZtm/Vs0MDv7Ow0hx12mJkxY4ZNOKfw4HVOgPB+t7h6KS64OtFBIeK23Tl3HvcE0Gc4Hbf5G5jLJGfxgHvtcMtpITQehKZD3Wcp8dEnFh0UAREQATu0nzCIgAiIwLARYF4ADDEunHQtgVmdu4ftYflTMY32ACaSo/h4Iy0KmJhthYDf6Oc2j/vPURzwntbWVvPWW2/ZMCyws3kdEyZMMMz74Gzn2TkfLveDJYWKO8+8EK7+fBAnQJwYSbeJbaQAacT9rxInRmQ6qABBCJa9BvfWQXxwGF6OdtbE+8frwt87FuV8jNcfgN5bBETggAQGnGR4wNp0UgREQATSBGCU2k/xMEqnwhAL0RjjUltbu+8T/dimFVq/fn0M7/0YVr6pzYdw4oOCg6vb95fAZBobG82WLVusp4I3UzzU1NSYqVOnWnHhRIUr/aLDL0T8woN1ZIdh+Z8LUYImJXnNxieeeKINI5jx34iDChC2L73YYWMhQvjCzTzmfgfp8+OiQP+hyyU+xkVn6yVFQAQGRUAhWIPCpptEQARyJeDLC0CkULh+HI30Yw13GPgPIjzqM+BlP/g4g58eCIYqUZy4kkwZSsWcD4545UKrKBqmTZvGfAzr1eBxrk7A8H63psWO7R4KGa7uOn/J57Mt7hifgX2Px7H9O1SQXLNmDf+NyCmRfJwPt2x58w+JjwwKbYiACIhAvwTkAekXjU6IgAgMBQEYuDYvYCjqKrA6rOEOg34tGOxA2+kVsKLEGf+upBDgNpfm5mbOPm63KQ4oDBDKZSorKzMznNOT4TwfzsPB0nk+3Lb/GtbjX50Q4jHXDpQhhF7Fcd+v2IC1a9fmJD74tR+XZ4ZbRru52JeAB8S+s32hMf6HxMcY72C9ngiIwJARkAAZMpSqSAREwE/AGZ4wROtwnEPxcibCLv81Y3zbGuWYOLAVBvrP4KGgAOFwxBmvgxMBFAD0VDDhvL7e5n7ba3gcIWumvLzcigvu+0UEhYhfjPj33TZL3uP23f3uONuQPobNBJ+z9sknn3w9HX6VkwBJ92MS78n8DxQBDyJmPPW1PB/pH4EKERABEciFgEKwcqGka0RABA6FQBzGN61SJia7r+IZb8ChVFwA99KA58hXP8K7L8V2KVYKkwDFhFsYPsWEcyaWu9Aq3GMmTZpkMDt2JuyKIoUrz1EsMFyLJe9H/bZkndzmwmu5zevdtjvuxAz33TaESADC4Vs4lMRs6zmNYoV68YiAx+GWcR/fL8H+RjluFsdAQ+2Omy7Xi4qACBwiAXlADhGgbhcBEeifgMsLwBWwgZO0isfbRw8KkNCuXbtehJH/MAxzMrBeBQoCigeKg7a2Nht65Y7xeFVVleEIYgyjokDwiQTnsch4NZw3g9c6TwdL3kNRwxGxWCfr50LhwvlEWKavT+J8CO1Yj7Y+TO/H6tWrbTvtDTn8gXqiqIP9S09XByZJtK4cipMcbi/YSyQ+Crbr1HAREIFRJCABMorw9WgRGMsE0oZnJi8AhhoN0Wa+swvPGsvvn363CMo45gO5GiLgAqow7Gc8C0RCgUABwm2KBIoGhlxxtnMnPnicYVLuvAuj8pcUEm7fiQ+GdDGnxO3zfrfQa8JnYLZyWy/bgev+EXOOdML7QaE0IOFAseQXG3j2gO537SqkEn2GV9ZoV4XUZ2qrCIhAfhAYb18j84O6WiEC44cAcz9o9RZhTcK4bho/r269PTEY+VchrOoOGKoUCAxVyoRDUTAwcZyeCJ6nCGGyOUcOo2jgcR7jytAslhQOXN09bp/3c+XCa534KC0ttdfyOneNu5clBAj7J4QQsLufe+65NYsXLw7B+7FPqdga+/9jRWq4Zc4fYodbxnM4lm836h2QB6X/J+TnGfQFcEp85GfvqFUiIAL5TkACJN97SO0TgQIk4Iwz5gXASCvFfoLWGgzezNf/AnytgTSZf7fS80HxsQor8yIYmhSgF4IYKD7o4QCWTI4GPR8UDxQfLCk43MJtigjey5LOFJa835Vum6KGs6fTw8F6/Od5jbsuXQYhFjwM83vqmWeeORnig3k69I4PSED4hltmHzcuXbo05n4H7h3GSuneSzkfY6VH9R4iIAIjTUAhWCNNXM8TgXFEAIZ0FMYaw7C4dGBt5Ouj3GdZjz0e+4kPvC8T8LlY4eDEh//VacDzvAu74jXZeR/c7+sYr+VK4UKBQ/HBuURYF8UIS648z+t4DQUNF4gTtstDsvsszLL+Kx6DF4SulJQ7hQdyWCCIMsMto84xKzQlPnL4MegSERABETgIAQmQgwDSaREQgcETYF4ADDYunAY9gTyD7sHXVhB3upyPjOcDxj3scWvkW08EBQSTwt0CLtZTwdAsCgMnJA4kNHhN9nknLNrb283WrVut4OBznACh+HACxF3LZ3OBV4X/FiSQ+H7WZZdd9gOGYEGE5PTvA0KwbCUQIHV4Tw7Dy5dtYr1jbcG7WbEmz8dY61m9jwiIwEgTUAjWSBPX80RgHBCAUWqTmGGU2rwAWG0xbHuY02JAX9ULCBXfiwa7y/mwYVfZ4sMJDL4XeFiPB5PA6f1wYoAlxQVLhlsxhIrXcoWBb0UKS3c/npEJseJxigt6O1i68Cse58pr/SvbwWek6wihzjhyUK6++OKLn4YIuWPBggXhtWvX7hsvmDf0v9jr8ByyaOZlK9K/g/5vKZwzEh+F01dqqQiIQP4TyOkLV/6/hlooAiKQjwRoWMMg5Rdy2MLh+mXLlnXTkMvHth5Cm/g+XBN439vwnhQfnJQPr77P80HxQQ+EWygwyAfXZsSH/5zzcDjPhd97QXGRve+Eh6uX5+kN4ehW9IL05QlhmyhkuDgRgjZThPzgwgsvvIjigyLEtau/cqwPtyzx0V/P67gIiIAIDI6ABMjguOkuERCBHAjAuM3kBeRweSFeQuHBv0eTCKFahfU6eBRswv3BxAdDruh9oFBw3gzUw9GnuLqwpkxIFq9zeRzc9q/+HA8KFySUW0FB0UER4g/F8l9L0ZIlQqxqgjBKIifk5+eee+6RFCGcFwRt6nOBcLECE2UFLsCr2IWJ7PSA2Pfo88YCOSjxUSAdpWaKgAgUFIF+/1EpqLdQY0VABPKKgDM8YYrWoWEcipcjQHXlVSMPvTFOfCQoPmC0XwVBEcMaOpD44GMpPuj5ABPr/YAo8CASgMnjvUzgth4VlFaM8Drn9XACwpXZQoQeD7TFq6urS+C+OOpN+kUIRYn/nmwRgudz4VwkE5ATctfcuXNL0/3JNvW30OPD/A8UAQ9tGxN97d5HOR/9dbuOi4AIiMDgCBzUtT64anWXCIiACFgCcRiztEoZi+S+itOQLfQv4/uJDwgKvmsE7+nCmax3gca+f3FhV058oOQcGh48F0GIg2XIkynF9dejrjoYwC48ywoRGsS4Do+wi62WYgF75OlWti2EIXhDnFOE84FAlHjQRfRu2FyQ9D195oWgLnpmgvCMMBTrdJTfx7GP95UP4gx0CKoaPJPt5nDDhd63eA3AVMK55aA/REAERGA4CMgDMhxUVacIiIAZw3kB/YkPDjc8GPGRhJFPg5+zpX/n7bff/iq2T4YBfCV+Rg+i5PDF9Iywfpb8e5ujitkV11KU8H4mkYdZYj8BL8Q6TEK4DPc8xGsw23qc3pHscCznDXGeENRrf718HrZj1dXVH0NS+gqGYsEbwlG+9ltQRxTP4Actero6Zs6cWc+L8OyCFCNky7bL87FfV+uACIiACAwJAQmQIcGoSkRABPwE0oYnDWabFwCDjoZoM69ZUdh5Ab3EB778c7hdej4ORXxQMFwNwXAn8NDAD+/evXvX9u3bf4z1EjA8Afj+Agb+t1H+AetLuKYRpV3T21uxz3O3Y/8TKE/euXPn/B07dnwHguAvcewtPCMMEZJ0AqSvpHRckwkLw3PpCWF7YvCELL/gggs+tH79+lhfSekcbjnd57jc0OtSkMKDbQc7vErAg4CuRujah7A9HcfaUVZi3YrcmLvdYAq8jvdoEQEREAERGBgBhWANjJeuFgERyJ0Acz8YOlSENQnjtin3W/Pyyv3EBwz0QxIfYBKC5+FqrBQf/Ps4ln5zPou5IAmIiC0ouf4KK5fQjBkzaAxb45cGM0RGFyYf3De5SOo61kEBsRfXLsF1D8BTUQ0R4kEwZcKxcC7jtXHbzE/Bu7njYQitJETIqvPOO+8VeEI2YI6QEOcKgZjkMxhCZodbRr1xbHdjaOEBzaKeau7o/8lkezBI3nLLLZPxLh9Ci0opPrCWg8fTNTU1fxjLM7yPfg+oBSIgAuOFgATIeOlpvacIjBABGsRYPOYFoKQBZ0eFggFXyLNjj4T48M+3QXHh9m3IFfbZBhr2iS1btjSizF54HVcuvI4rE1DCEDpPAf/V2F4Nz0cc22EY2Jl8ELeN/rKiw5Xp4/a5COeqgRfmp/PmzTt/9uzZfD6P2wXXoZvtwj5uLEQjneJj5cqVVnzg/RfjPaJ4I4q6Crzfk5/73Ocesi9rI8vk+UizUCECIiACgyLg/rEa1M26SQREQAT6I8C8ABhwDE3iwjwGazSjpHFdSMtIi49sNlZ04CAFCbfZnr5WnuM17jps2oX7EYiHu1GuZJ8gFyTWXz4Ic0H6yAehsIhXVFScNHHixJ/TUEcoll9QVsHzYfsVxrr/uG1Avv+RLT7wG+VvtxvtrsD2kzfeeOND2HfMC+33m+/41T4REIFxSEACZBx2ul5ZBEaCAPMCYLRxYbJ0AgYvDbpCW4ZFfEAM+MOuKBAGstAA7ms9UB0M7WKC+gr0xS8hFiLIb9gvKR2iMTPBIfNBcD370NZLrwk2Yhia96JFixbdxqR0hIdZsYFrqmGoM+SOIUxN9oYC+SNX8YH3cswL5M3UTBEQARHIXwIKwcrfvlHLRKAgCWTnBcBwi/HrOIaXpTFfSMuwiY90wjn//h2o+DgUfszHYb4I5ys5AgLjdIiQJLbhtEiFY6GvMiFYPOYW5oRwwTHmlNATch1Gxnpu1apV/459crLvgfPcbsZq3O+A2/m6DFB85OtrqF0iIAIiUHAE9v0LU3BNV4NFQATymQDnu4BByq/GsF/D9W7koHxus69tY0188NXYF3yvRng2lqDcS48Hk9I5MhZXCKPMbOvcdl4QekK4pAVKCH2bxESF30MY1uk4zDlMOGcJJz5k/QXxYUvigz2qRQREQARGh4AEyOhw11NFYMwTgDFa7AzXAnvZA4oPega48t1owPuXfiYZ5AhgIVzvhtodac+Hv4nME+Hz38LK4Xk5+WECnhCPw/K6IXopPvzhWHxX9Ccut14Q6JCAQYhdKUTIPZdeeul0jHpFDwjnAOGym9fBA0LBk5eLxEdedosaJQIiMI4ISICMo87Wq4rASBBwhicM0To8z+YFYLtrJJ49BM84qPig8U3PwEDEB4faHaWwq76QUCwwH+QhvMv13IbwiEOEGIoQt/YlQnz5IMz1SGBo3snwbv0f6oImC1qFgvvyuq8lPvr6SeiYCIiACIwsAQmQkeWtp4nAeCLAOTLcnBPuqzgN/HxdchIfHCHKLz4oSKLRKOPMMpP4QXBxbgzr+aD4wOrm+RjJnI8DcbYiBELqdlz0H2hqBALEJqX7BYgTIXxnrvSEOBGCvg1RhGB43nc888wzV8IL0oZjB3rmqJ+T+Bj1LlADREAERMASyO9/LdRJIiACBUkAhh6/sNtZ0GHc0rDP97yAnMUHDXGKDrdgvhMbkoX3tMY5ynwXH67pTOwIop8+jfIJlHamdA7P6xch/lAsJ0JcBRQhEF4JHJ//3HPPLUJy+u49e/Y08TwY5VUIlsSH6zWVIiACIjD6BCRARr8P1AIRGFME0oYn5/+wAgRfzGmINvMlXXhWnr1wL/EBbwZHicrMcM6v+hQcNL6zxQc9HzxfgOKDXeAEQgfEx3uwvxvvF4AnJOlyQbKFCN8f19qVFZAL3j/ImdKbmpo++NJLLx13xx13tPEUz+fLIvGRLz2hdoiACIhAioAEiH4JIiACw0GAuR8JrDREGYrUNBwPGYI6e4kPGNJXwZ4eD+LDoWPeBufy2A0x8X4ehOjw3MhYFCAMw3KhWPSGOBFC0cUF9wUg2AIIwUpiXpBPX3DBBUfisLd48eK8mJBQ4sN2k/4QAREQgbwiIAGSV92hxohAYRNICw4DQ74Ghnwp3iZJAxVhOnlhjGbRzYgPGM+rKD7GiecjC4NhKFYIIvFx9N/fchsekCST0p0n5GAiBH0dwMqRsSrB8a65c+eWzp49mx4WMh61ReJj1NDrwSIgAiJwQAISIAfEo5MiIAKDIYAv5VEYswzD4tKBtZH1oHRhP4Opdijv6SU+MHzueBUfjilFSBhejTtQ3oG+CyEXpNdM6fSGuHwQ5wlxoWeshBMaon8TmCn99Lq6uu+vXLkyiXlCRk14SnywV7SIgAiIQH4SkADJz35Rq0SgoAngSzgTsrkEsCZgzHbn0QvtJz4GEnYFI7tQcz4O1gVWhCDH4xr02UPppPSMCKEXxIVkOSHiwrFwva0bHEPYjkGEfOySSy65bu3atXF4Qzh7+oguEh8jilsPEwEREIEBE5AAGTAy3SACItAfASSZ25AbfBmfCkOdw7TCHvW82traUQ3F8bX3kMUHXqfQRrvyvf4BN6kimNjBCQU5SeFbCMsKu6T07HAsig+uuCYzSSHFGTwh9HzFMCLWbcgH+dD69etj8ISM2ChoEh8H7GOdFAEREIG8ICABkhfdoEaIwNgiwBnB8TWcBi3SP8L1y5Yt64bhPtoiROLj4D+zVGa5MXshIpbg8g54PgL+pHQXiuUS050I4ehYXHAf80HCyAVJYqLCVeedd95J9ISMRFK6xIftAv0hAiIgAnlPQAIk77tIDRSBwiOAL+jFziDNk9bz7zquCSacM+djoGFXY9jzkd1FFCFheDaewjtfi+0AREeCSekUH271CxAXisWcEC5gSxHiYZLCGoRj/fScc86pGe6kdIkPi15/iIAIiEBBEJAAKYhuUiNFoDAIIATLJgPAEK1DizkUL8N5uka59fx7jpaxxEfuHcGZ0iPou39HuRL9GEYeT8wvQihA3OrCsSg6KdS4QIAwAT2OUKyT4Am5k0npyAdhKNaQe8IkPkhciwiIgAgUDgEJkMLpK7VUBAqJAOfS8NJfwnez4S4/ZIRfwomPGng+7jtEzwfs8STnNAnB4L4a6514FxrUNNbH4hLDS0UgKlZAVPwSrx6BAIlThGTng7hRsZwnxIkQ5oNgOwYB8sGLL774NuaDpEXIkPGS+BgylKpIBERABEaMgATIiKHWg0RgfBCAQRiG0VqBt6W9zq/dI5aAnEWYz6bnowYzlt8P8bEIgijnSQaRymC/5tOY5otgCY4j8eFQUlzRi3UVyg14fyalJ5wAOVg+CCuBCImQO0KxroMI+fhQJqVLfJCwFhEQAREoPAISIIXXZ2qxCOQtARjtNvEcpRUgMN6538wGu/CsEWo8xQfXCWnxMQ9GcAwrR2hieBCTpe0oTvxqz2234Hq7z6ZniY9GGNwfHgeeD4eCJfuPcBrB7KMo98LbsV9SenYoFoSKcTlAZIs1BAGYLCsr+x5GxDp9KJLSJT7QG1pEQAREoEAJjNaXyQLFpWaLgAjkQIC5HxwSqQgrQ5aacrhnqC/hx5UExMS3MBqTEx8RGsO5ig989beeD7xLEKKjEesiHHsS9Y7lsKu++oF9Gcb7b4A3g8PzPgjBEcfIWCGwBFIrMCxXx9bPmcd4HTgyKb0UDO+aP3/+6UhKb0JdLkQOm7kvEh+5s9KVIiACIpCPBOQBycdeUZtEoAAJwL60bgQY/DUwOEvxCklapxiGd6Rnw2Y7Ekh+noCv7u9HW5JYw1gPSXxASFF8cFK9sZrzgVfrd7FJ6fBqPIQuvR59zaR0mw/CMCx6QFw41gHyQawoLC8vP7Kmpua3SEr3MDQv+8r+bvp9ctaJbPGBtkSxcqLLCrTtyRtvvPEh/haxpLLhs+7XrgiIgAiIwOgTkAdk9PtALRCBMUUABmgUIoShTsy/6EDZyBccQYPQGroQHCfhmeVY7Rd4llzxJd+u3HaLP+yqL89HWnzw70smZo/Xhe/O4XlvhyfkVBj5H2NSOratuEtztowda1fiGnscZQgiJl5dXT1/0aJFP1i9evVSTlLIkKxcoKaFRfLLX/7yEbj+cuxHUfZgrcKz1kl85EJR14iACIjA6BOQB2T0+0AtEIExRaCkpIS5E1wYdpPAl3J+nR7JxX75hkH6FtZOCBEUAX4Sz0l8oM3QIMlM2JVPfORkJI/ki47CsxiOFYSI+DTKp1DamdLRx73mB6FHJNsT4tpKbxT6gknpVyMpfSnFRy6TFNLzwX685ZZb5qCOD6M+hvjFcYwi5A8SH46wShEQARHIfwISIPnfR2qhCBQEASSZW5cCjPepMApDNBa51NbW7nM1jMyb0PMSbGpq2oo2vARjlV/fkwfyfKS9HjbnA0Y1R7uyOR8SH/t1mAtr6kBo3RJ0bz2FBjwhXq4jY6Ev6JEK4f4EwrH+9cILLzwLnpDEgUSIL+xqDlr0LqwUteznKMLs7rvpppseR1tQtcKuwESLCIiACOQ9AQmQvO8iNVAECosADEIamDRUYWOG65ctW9ZN43CE34LPo4H6fawBGMkeBQiNX7e4sCuKD7TPP9SuxIeD1HdpBR6E2lvvec97yDeI/A8vW4T4vSAQcnZULLLmQqXABXOzGIyM9euzzjqrjiKEddkLfH/0Iz54BT0f9/7TP/3TC8474rtNmyIgAiIgAnlMYL+/7PO4rWqaCIhAARCAkVkM43S0W0pLN/D3f//3/wE79yWID+YpZBol8XHI3ZO89NJLS9///ve/NW/evF/SawQPSBIiZL9JCv1ChL8LJ0LQHzakCgKkbsKECb/BBIWlEBJsWEYlHkx8fP7zn7fig7OsH/IbqQIREAEREIERIyABMmKo9SARGNsEEIJlw3NgYNbhTTkULyew6xqNt6ZHA4tXV1d37oc+9KH/Qe5zB77CMywsDvFh80FoCKc9HwlsK+xqgB11/vnne9u2bSv71Kc+9dDMmTMfh8gLIReEOT+9RAhDtPz5II47H4d+oQiJY8SyM9BX36eQQFI6R00LSHwMsEN0uQiIgAgUEAEJkALqLDVVBAqEAGcb9/iBG+tuttnlh4xA+5l4Aps24H3lK1+5aO/eveeccMIJ2z74wQ/+C5Ked3B0LrYLRrCdqwQFxUdIOR+59wz58uqWlpYJ4Fy6Z8+eKDxN/wUR8TJFCLwgdqZ0f04IvSBcGQbHlZ6QtEikCLGji6F/PnbJJZdcl05Kj1CMMOEcj3I5H3ysDbuS54MotIiACIhA4RKQACncvlPLRSDvCOCrdRjGZQUaRhufhupIDvXN51F9WPEBA3ceclCaGxoaKhHe8xqSnM9F25bj+Nu4LogylBYfm2EUL1LCOagMYgE79nnstNNOOx+3v4H9ECYpTFKAcHXio79QLD6SIgT9FkdS+m0XXHDBxcgH6fnCF77wDpyS+CAgLSIgAiIwxghIgIyxDtXriMBoEaDhj2fTkLQCBAY+95vZHheexe3hWPAoKz7uuOOOyFe/+lXOWH46ntOKtpTA69G+Y8eO3373u999AyNj3YxrT8D5BRAcH0B5PozkOWnxwb8PNdRujh2EULZicudoVgi7Kn7sscd2Yp/D43ZQbLikdOcJYRiWEyHcBvNMPgjuobcshInSPXhSfnrVVVd9EJMVvhP94/pDno8c+0WXiYAIiEAhEJAAKYReUhtFoHAI2NCmtCBIwshsGu6m81kUP9/+9rejjY2Ni2G0nob9NqxlOLcdIzT99NZbb93BnAK0JYSwrFYIkUfwlf7XMJLX4FgbVp5TIjMgHGxx4XQQDHVYmcMRwMhnHUxKx73rcehalExKT/iT0t2M6S4fhAKEK/qIw6XxsQH0iQdvyqS33357FbaLIXJiOF+Mc/cq7IqItIiACIjA2CAwkuERY4OY3kIERGA/Ak4EwNtQA3u0FPsJGqYIrWFC8bAt7rkUHzB2P4RHTsexFhjB5Si3YYSluzFHhB0GGOecwKC3xP/xhcfduWFr6xis2Ob6gDP7uPGee+7pQKhbZP369f+Obj8Fx6+DZySG7Qh/BuiTzIq+cDOjW68IPFQGwpDJ6gyN42hlta2trVd84AMf+DGEy71f+tKXNNrVGPwB6ZVEQATGLwH/P8Ljl4LeXAREYEgIwICMwoBkGBaXDqyNrBilHZZqSB6SrgTP4TOs58OJD3g/2vziw81BkvV8toVGrluHvG1D+Z75VpcvnK4KvC07MLdC86mnnmLIVARejOvRP7/E+Qj6Js6RsegBcV4QhmIxGR35OeaVV14xGE3LjpSFfrKTFKJMwEt1yoMPPrgZ4mMDRsYKa6jdfPslqD0iIAIiMHgCEiCDZ6c7RUAEsgiUlJQwpIYLR6NibgBnrB7yhfVTVPg9HzmKjyFvy3itEH1QjT6wwy2jbPJxoAjhEMxXodyAMKswwqns8LxOhDAMCwLDbN++3YZhcfJK1JGpAts2p+fVV19tZl9jiF6JxAwdbYiACIhA4ROQACn8PtQbiMCoE3B5ATA6p8J25HwbsBs9r7a2dp9VOUStpEHK+iU+hgjoIKtBX9ODRI8F+7iZ2+nfAcUCjzXi3EdZQnAE4AnxXEI6PSIIsbIhWejLzJC8uDazsJ8xKlY7+zpzUBsiIAIiIAJjgoAEyJjoRr2ECOQHAX7JhtFJgxF5xeF6FwI1VK2T+BgqkoOvh4KAo42hhhpscx4Vio3sfEKKkzDCrDagvBprEN6POJPMGX6Fe3p5PHA+s6CPE/gN0YOy9rnnnnsKJ4IYlteKncxF2hABERABEShoAhIgBd19arwI5BcBGI3FiP8flkZJfAwL1kFVCjHBZPEohQRWLm7CSb+3wuWD3A2hwpyQCDwfTFy3no9+Hpxyb+FiJK7/A65JYvSyfi7VYREQAREQgUIlIAFSqD2ndotAHhFA6I01PGGI1qFZNi8A211D1USJj6EiOTT1VFZW0uthF3pE4NXor69juCiMPJDb8Xv4D4iPSHFxcQL9mb57X4FjVK7JtPdj6YYNG9Zj8siQks/3MdKWCIiACIwVAhIgY6Un9R4ikB8E7NCsMCI5sZz7Kp4xVgfTRImPwVAbnnvYF6x506ZNExFiV8oQLJT7q4nej6ewoMfkWgiRx7HNEbN4jPfhsEdBkoDHI4TfTAgetKs3btx4J0e+UugVCGkRAREQgTFIQAJkDHaqXkkERoMAQmXCMB4r8Gx87O4zL2DAzYJhChtXCecDBjcCN0BMcKSBEMqOmTNn1vOR7Ks+Hu1BTPDfmnZMVvh5CJYe/D5C6b61c8XgGPdfxPFLnfhYu3YtQ7i0iIAIiIAIjEECEiBjsFP1SiIw0gTShifn/7ACBMYkDdFmtsOFZw20TWkDVaNdDRTcCFyPGcqL2T/uUS0tLX0JD3uaM9BTTKA8GQLkpPnz53+3tLR0M0QHJ4jcg/X3EB6fRJ3vePHFF+/HTSGJD0dWpQiIgAiMTQLZI5eMzbfUW4mACIwEAeZ+MLSmCGsSX8abBvtQGrcUNRpqd7AEh+c+iEmKDg+hUnXoIxRBuru6kdfR50zyFB/M4bjlllvm4L5LMfdH7IQTTnjzmGOO+e7mzZv/uG7dupcxc7oVqukWu/Cs4XkB1SoCIiACIpAXBOQByYtuUCNEoHAJUCyw9UVFRTUwSEuxmaR6YFjNYN5K4mMw1Eb8Hpvrg75iHzcuXbo05n4HriVZ4uNdON6NNI8k5gKpQNL6A9///vfXpcVHkMnmOM/fEQWsFhEQAREQgTFOQAJkjHewXk8ERooAJpuLwghlGBaXDqyNfDbKfsNzstsm8ZFNJL/2feF0VfB82H6F6NxPaPYlPtJvEsXv4bdf+tKXNvhERzKdbJ7z7yS/qKg1IiACIiACAyUgATJQYrpeBESgTwIlJSV2SCOKCKwJzPnQ3eeF/RyU+OgHTB4eRl9VQ0jY4ZZRNrGJ6fAscyDxgcvu/fznP/8Cr5HoIDUtIiACIjA+CSgHZHz2u95aBIaMQNrw9PBFfCqM0RDWGL+O19bW2tCsXB4k8ZELpfy5Bv2bQD9zQkH2cbNrWbofXc6HDbtKn4uizIgPze3hiKkUAREQgfFJQB6Q8dnvemsRGHICkUiEBinDaJD+Ea5ftmwZRzk6qAiR+Bjyrhi2CiE6vDvuuCOCB9RgOwEhwv61H7JOPPFEHAp4SDg/Gcec+MBPIliMfYkPQNAiAiIgAiKQIiAPiH4JIiACQ0IAxmgx5gGhCMm5PomPnFHlzYVdXV2cVDCK/k5N+BII7EbjAkuWLEl87WtfuxCjn82DEOnANRGU3fg93HPjjTe+5kKz8uZF1BAREAEREIFRI5C7pTBqTdSDRUAE8pkAQrBs8jAM0jq00+YFYLvrYG2W+DgYofw8X1lZ2curBYHBvvZuvfXWSyBAz8B2uxMfKFd/9rOflfjIz65Uq0RABERg1AhIgIwaej1YBMYcATs0KwzSAFZ+FTfp/JD9XlTiYz8keX+AfcZGbtq0aSL6txSbyfSwuj1f+cpXzoO343Qca8W5Ino+KD6QcL5Lng9S0yICIiACIuAnIAHip6FtERCBQRGAkRnG1287Czq8H5m8gL4qk/joi0phHUv3MbRGsBm5P2ei9WfhWCuEB4dilvgorO5Ua0VABERgxAlIgIw4cj1QBMYWARidNvEcpRUgMEC538y3dOFZ7o0lPhyJgisD9Gax/+D1KEFfBzDxZOyll146GfO/nIzddno+cF7io+C6Vg0WAREQgZEnIAEy8sz1RBEYiwSY+5GggYqXSyIRuSn7JXkOi/ftb3872t7e/iFsT8dX8zYYruU4t62srOxuN3IWr8u+X/tDRsCGUqE2KyjYL1wZKuVWd4xl+qkeh85lv5SXl1ehrxLPPffcCXv37p2GEc/acJ08H0PWPapIBERABMY+AY2CNfb7WG8oAsNGgAYqjVJ8Da9BWYp9zg/Br+S9Zsd210l8DEtXUCRQsFFIZB7gz79xnih0DS9wF7HzMtf7jhuIjcxxhtdhJ9zS0hKYM2dO+csvvzztjTfeOL25ubl2woQJW4uLiwOxWKwH1yjnI0NNGyIgAiIgAgciIAFyIDo6JwIikBMBhOFEIULCMGiTuKEDZSNvpMHrxAcM2Wp4Pt6NY9NxrB1lJdbNpaWlv5TnoxdmKyjIzR3NFhPgxlnnLV9cM1BBUYxJIj14qYKdnZ21qCeIPA4+awrWIMQEczvsNkpOMMlZz6MTJ070Ghsbi3bt2nU8+rEM5+LQmWHcT/HxSyWcg4IWERABERCBnAhIgOSESReJgAgciEBJSYmBQUtDmIZsoqOjo5vXM6QHxitnxp4Mg/VDOEQvCYdpLYfx+nRNTc0fli5dGsM+LhuTYVdWTOC9D+qdIK80Aysosng4kZHxTuA8RUjghhtuKJ45cybD3sogHspRDTRDshjlJJQEGwV7N0RyEY7XtLW1eXRS4XgZztu24doQ9j2EVOESO5wyqvd4DUPrPJQBCM3E1KlTdx199NGb4Qk5tbq6egc8IKshIPewrzXDOdFpEQEREAEROBgBCZCDEdJ5ERCBfgmkv8zzK/lU2LEhrDEavfjCTqOWtq0VHygXoxLmCXSirIAx++TnPve5h9IV87qMgd3vw/LjRM6CIv1O7r0GFO5EDwVEQjVEQhSigqJgIgQGk7+TYFcJjtUUBl/96lfLJk2aVA6PBHNwokAU5TVYAxQUWDxODolzCeLDYWza7QCPY7+H16XPUQjafxN4PxZ7HOfDfBZvxLEw5gHZW1FRUX/OOec8vW3btmckPkhPiwiIgAiIwEAISIAMhJauFQER6JMAQnisB4QGLIzU+s985jM9KD16PlAuxmqTlHFzBbafxMzYD6XsWdrzmRCiPuse5oNWULAt7jn5EO5EDwU4laBdVkSgbaG0Z4LuDZ6j4GCTGfJmBQePYZueJ85UTo9FnPejpNDg+9lBR7AZRh1UJxQmrL89LVIobhpYPxZespN14hhVyx5c081QLbSjEfs91157rfVyYZvV89laREAEREAERCAnAhIgOWHSRSIgAgciAGO1mEZs2tilJyQX8eG8AweqeiDnrJjADTakyN2YLSjccbYR27YN6W13KtMul4yN8yMe7gSm9FpYcQEjP442hNLGPtsC1EF6Ofh3uBUhOBbBdjeuo5cpiHPduKaJ2zjXif163gcB0YmynqKxq6srGY1GGxBCl2xoaOBIV8znyHWxoi2LXa736joREAEREIFxTCDz1W8cM9Cri4AIDJIADFnYnwEPM2EvxvYR2Kb4eBrJzeuRlP5RGMB2bghUn+35yBj5B3j0foLCLyZ4H/b7EhAHqHLfKTe6U67hTrizDO9ocyxQ7hfuhHe14U4oe4U7kRGfCi4MY3LhTmx3JtwJ29Y7AZHAWeRduBOFFMWCHVkM1zSiDgoECopmXNsCMRGE8GuCF6MFoVohJPS34lg7xASfGYegiKMcyJIRb37WjjMrQhtY5NJ/vE6LCIiACIiACOxHwP5Lst9RHRABERCBHAjAEIY9agXIB2AQHwNDmAbvKzCKp8CQ5ghLNucD5VMYJelBJir7jdkDGLkDNnBRtx3dCV/1KXpqmDuRy+hOuJYvwfwKehi4nRlCGO/EfRdeZD0N2Offm+6YDXfCPgWEC3HiebuN0oCFC3ei28KGO+FwJtyJ21gy4U7Yz4Q74VGNSPLuoaCAmOjCuQEtfB13wwFY85IB83b1qhQBERABERCBgRLI/OM00Bt1vQiIgAjQwF29enUQIyL9NSeo27x5cw28HzUYJelNJE0TUBlECRPOHxwILdbrH90Jo2pVQEzAlk9Uor5q1M3QpP1Gd8J9VkDgHPMnuG2TsWH4OwPbJjiwLTiVgOFvw8aoAHgtF1cHb+d1OIVHhlhyn9exDgoMhjvZnAuUbdhvx7pfuBPaWs9RwhjuhLIBHorBhDuhaqiE1CuxNH0JCrTDXoY/uOHemce0iIAIiIAIiEDeELD/WuVNa9QQERCBgiNw/fXXl8yYMeOKPXv21G3YsGEutl864YQTNjU1NVXDIH4Cno+19E7wxQ4S7jQR15fgMnoX+gx3Shv+0BN9hzvhfisgcB2Nb+eR4Pagw50gPHajuu7+wp02bdoU/Na3vtXFR+M5A1kGEu4kQTEQsrpWBERABEQgrwlIgOR196hxIpD/BGCchzEU6z/DC3IOREflUUcd9QK+9MewtmJ+inXwWjAkagJW+/cNRQbW/cKd8KYZ7wS3ueI63pNzuBOeYfMlWGJl+BeX/cKd4E2x9ePcTheqhUcp3Mni0h8iIAIiIAIiMLwE7FfB4X2EahcBERiLBCgosHj3339/JcKLqiA0gt3d3UVbt26dfeSRR26cNWvWLkxcd2w6/InxWFaAUBxg24oLbGdGdyIjX7gT/26yIgHX9BrdCc/sN9wJbeDoTgy76kZIWFM6GTs5wNGd2BQb7sTSH+qU3reeDrSDu9zmBmH4PSCZbTyb12kRAREQAREQARFIE7D/goqGCIiACAyUgBMgDz/8cPGbb755+XPPPffOF1988aS6urptZ5xxxrPI2yiC+LB/xzB/AtsHHd1J4U4D7QVdLwIiIAIiIAKFR0ACpPD6TC0Wgbwh4ETIv/zLv0zB7NinP/vss8dNnDhxK8UGxARns9s5iuFOGe9E3gBTQ0RABERABERABFIhEeIgAiIgAodKgGIEdQQRimTnwRhIfel7cw53GkjdulYEREAEREAERCC/CMgDkl/9odaIQKES4N8lLu9BozsVai+q3SIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIgAiIw5gksN0HDVYsIiIAIiIAIiIAIiIAIiIAIDDOBXrPIL1++POh5nsTIMENX9SIgAiIgAiIgAiIgAiIw/gjcZUL2pR8xl5u15t2mlxQZfzj0xiIgAiIgAiIgAiIgAiKQKwF9sc+V1L7rnNyImLj5gYmZHyz2Focuf9/lS++4445PwQvC8+6afXdpSwREQAREQAREQAREQAREQPkLA/4NeGC2xCTM35sbzOHmsOARwemrP736ypPqTjp+xtEzQoFAwLvrrrsk7AYMVjeIgAiIgAiIgAiIgAiIgAj0JsCEc/o2FphTzb8Zz2w3XmBHwDN3mO5P3PiJK+D9KFVSem9k2hMBERABERABERABERABPwF9qffTyGXbw0VTzHvNTGR/JE09REe9Ocb8+ccNPy6D96PDnKjwq1ww6hoREAEREAEREAEREAEREIFcCaRS0A3yP75nV9y33LOjYCn3I1eGuk4EREAEREAEREAEREAERCBHAk+ZCAKwgli/i/U76e1wjnfrMhEQAREQAREQAREQAREYtwQUgjWYrp8LyRFAAFZqCfi2B1Ob7hEBERABERABERABERCBcUNAAmTcdLVeVAREQAREQAREQAREQARGn4AEyOj3gVogAiIgAiIgAiIgAiIgAuOGgATIuOlqvagIiIAIiIAIiIAIiIAIjD4BCZDR7wO1QAREQAREQAREQAREQATGDQEJkHHT1XpRERABERABERABERABERh9AhIgo98HaoEIiIAIiIAIiIAIiIAIjBsCEiDjpqv1oiIgAiIgAiIgAiIgAiIw+gQkQEa/D9QCERABERABERABERABERg3BCRADrGrQ3ET9IwJL19hwigDh1idbhcBERABERABERABERABERCBLAKeCdsjMfM9zIn+TbudOoJdiZAsWtoVAREQAREQAREQAREQgQwBeUAyKAa2YcGFTfuCP5nTIDq+2BU2N6CsCAQC3l2LF4cGVhuvtvfIgzJwcLpDBERABERABERABERABMY4AXhA6PAIeGblLd81HoSH54WNlzBmHbbL+fYoByEmeMtgxAufqEUEREAEREAEREAEREAE8p+APCCD66NkzJjgOY+YBe1eyNtiwt3bsUKXnPGvhx116zd/9Kuzl5i7Isa8UAQpAq1yoPUpXMf1r9+H60oDgdXQMR76Zbn6ZnB9o7tEQAREQAREQAREQAREYIwR8Oj8MNHFvzEvLlplYv+fMbE9ERPDse5lF53c6XnPvoNvHByQhPjPx4350WuI5lo0xmjpdURABERABERABERABEQgQyCdOp3Z10ZuBEJBE+iec3pw7c5/vXK2Z042HYm4aastNpsix5ujj57wj8Y0N3heIInq4vb/Bo4Nw91urF1Y4UOx292Mu2ox5tjKkpJ1syZN6rxvy5ZHf2HMy98y5lNP4hwW6h0tIiACIiACIiACIiACIlD4BCRABtGHywMmucJ4gUfOm/rzGzYt755knlz0RLK09aHGSx4wvzNeKGCqITYCnrf5SGNajzCmDAIiDKFBlwhXpokwT53bAXhKAuFkMjn7qKPeSBx33FmhLVtKPxwOT30qHg88kQrHYkqJFhEQAREQAREQAREQAREQgXFPgO4LqouJKKnmAjwwwCV1zwPPnX/+17xFi376U2POOyEVvuUNorYBPlyXi4AIiIAIiIAIiIAIiIAIFAaB5cvpwmCCOaUH1/2SzeHmYEL5gdaHceNduO76fy0qWvYX+/JGeI8WERABERABERABERABERABEcgiAA/IoXoqfPdTeCyX+MhirF0REAEREAEREAEREAEREIEhJUANojlAhhSpKhMBERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERABERhrBP5/WKMEJu4Z3SkAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "view(diamond, show_lattice_vectors=True, direction=\"tilt_x\", guess_bonds=True)" ] }, { "cell_type": "markdown", "id": "d6f61175-422b-403f-9e6c-fe926df73c38", "metadata": {}, "source": [ "## Run a volume scan for the diamond structure" ] }, { "cell_type": "code", "execution_count": 3, "id": "5566346f-182f-44e5-a003-fbc70d7fe1ac", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[22.06|14:05:56] JOB DiamondVolumeScan STARTED\n", "[22.06|14:05:56] JOB DiamondVolumeScan RUNNING\n", "[22.06|14:05:57] JOB DiamondVolumeScan FINISHED\n", "[22.06|14:05:57] JOB DiamondVolumeScan SUCCESSFUL\n" ] } ], "source": [ "from scm.plams import Settings, AMSJob\n", "\n", "sett = Settings()\n", "sett.input.ReaxFF.ForceField = \"C.ff\"\n", "sett.input.ams.Task = \"PESScan\"\n", "sett.input.ams.PESScan.ScanCoordinate.CellVolumeScalingRange = \"0.8 1.3\"\n", "sett.input.ams.PESScan.ScanCoordinate.nPoints = 11\n", "\n", "job = AMSJob(name=\"DiamondVolumeScan\", settings=sett, molecule=diamond)\n", "job.run();" ] }, { "cell_type": "markdown", "id": "c66ac952-7ab6-461c-9dca-8dbd5eea69e0", "metadata": {}, "source": [ "## Extract the results\n", "\n", " - Obtain the volumes in Å3\n", " - Obtain the energies in eV\n", " - Plot the results and fit with the Birch-Murnaghan EoS" ] }, { "cell_type": "code", "execution_count": 4, "id": "3e4cc7b8-81e1-4f56-bfff-15d2c2e96eb7", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaoAAAEgCAYAAADlpDdIAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAAsTAAALEwEAmpwYAAA+KklEQVR4nO3deXwU9f3H8dc7d8IVjnCFy1bF/rzA4lVra62K2loQxap41Nt6VUFUPBEPUDzBoqIWvGtVQFEUj9ZqPYsCCrV4gUq4whEgECDH5/fHDnYJm2QTksxm83k+HvPI7sx8Zz6z2d3Pzvf7nfnKzHDOOecSVUrYATjnnHPV8UTlnHMuoXmics45l9A8UTnnnEtonqicc84lNE9UzjnnEponKueccwnNE5VzzrmE5onKOdfkSeog6T1JX0var9KyTsGyf0r6u6QuYcXp6kZ+ZwrnXFMn6VKgHHgbuNHMBkYtSwXMzCok/QHoZmY3hxGnq5smf0YlaZGkw6pYNl/SIfW93aYsWY/LNXtfABnB9J/oBWZWbmYVwdNWwPxGjs3toCafqKpjZrub2Vthx+EigiRZIqk4arpvB7Z3kaRZkjZLmlxp2VuSNkXtZ0FDbUtSO0lTJW2Q9K2kk+t4PK9KGhVj/gBJyySlNWT5YN3qXocql9WwzV2C1++JSvNPlPR58Lp9LengqGVPSFoqaZ2kLySdXcNuXgGGAx8Bj8eIoY+kD4GLgE/ijDv6/bpG0suSusdTNihfXGkqlzQ+avkOv9aSMiU9Erzv1kuaI+moWsZR28/KiZI+DP5vK4LHF0hSjNdtuaTJklrG+7rFktSJqirxfGDDksix1ZNjzKxl1HTRDmxrCXAz8Jcqll8UtZ/eDbitPwNbgE7AEOB+SbvHdwjbeBQ4ZesHPsqpwJNmVtbA5aH616Gm16gqfwb+HT1D0uHAbcAZRM5yfgF8E7XKaKCXmbUGfgfcLOmn1ezjYCKv/2rg+MoLzWyOme0PXAeMqEXsx5hZS6ALsBwYX8P60fv84X0OdAZKgGejVqmP1zoN+B74JdAGuBb4m6RetYgD4vysSBoG3AuMDbbVCTgfOIjI2exWW1+3fYB+QVx1liyJal9J/wl+9UySlAXbVnMFj6+U9CmwQVKapO6SpkgqlLRK2/+67yPpU0lrJT2zdbtR2xseLN8Q/KrpJOmV4JfNG5LaBuuapJ2jyk6WdHOlbVWObZGky2PtX9JVwS/Q9cFxHxsdtKR9JM0Olj8blI2uk6/uuKrcdnUxhcHMppjZNGBVWNuS1AI4DrjOzIrN7F/Ai0SSQ6z1u0p6PnjPLZR0SdTiaUB7Il+6W9dvC/wWeCyOcHa0fLWvQ11eI0knAkXAm5UW3QiMMrMPzKzCzArMrCBqX/PNbPPWp8H042p2NYTI2dQjwDZntJKiv0DXAhvjjT8qnk3Ac8D/1bZs4DhgBfBO1DZ3+LU2sw1mNtLMFgWv40vAQqCqpL5dHPGS1AYYBVxgZs+Z2XqLmG1mQ6L+X9HxFRA5290jajvVfn/FkiyJagjQn8gbeVeqzt4nAb8Bcom88V8CvgV6AfnAXyutfwJwJLATsBfwh0rLjwMOD/Z5DJF/yNVAHpHX9hLi90NsUb98q9r/10S+jNoQ+cA/oaAnU/ChnApMBtoBTwOV3wjVHVeV266prKQJkibU4ph/sCNlqzFa0kpJ76qObZVxbGtXoMzMvoiaNxfY7oxKUgowPVieD/wauFRSfwAzKwH+BpwWVewE4L9mNremAOMp30Cvc0ySWhP5YhtaaX4qkV/ZeZK+krRY0n2SsiutN0HSRuC/wFJgRhX7yQAGA08G026S9olapY+ktyX9A7iUyNlAbY8lB/g98EGl+OJ9LU8HHrMG7r0mqROR92RV7XBVxRHPZ+VAIBN4oRbxdAeOBmZHza7pO2Z7ZtakJ2ARcH7U86OBr6OWHRb1+Myo9Q4ECoG0arZ7StTz24EHKi0fEvX8eeD+qOcXA9OCxwbsHLVsMnBzpW2dWZv9V1p3DjAgePwLoICgR2cw719b91eb7cbYdq3KVvGaFhP5hb11Oqce3gM3A5MrzdufSJVSJpEP53rgx/W9LSIfuGWV1j8HeCvGtvcHvqs0bwQwKer5z4PXJSt4/i5wWS1eix0qX93rEM+ySuvdC1wZPB4JPBE87hp8JmYRqVLrEMR5S4xtpAbHdC2QXsV+BgBlQKfg+Xzgjnp4X0W/X0uJVMftWYft9CTSI3Gnhnqtg3XTgTeAB2sTR7yfFeCUGO/194LXpwT4RYzX7VtgApBdTdxzCL5jqpqS5Yzq+6jH3xL5INS0XnfgW6u+3n5Z1OONQOUGweVRj0tiPK9NA+L3MebF3L+k0xRpNC2SVETktLpDsF5XoMCCd0AV267yuGrYdrVl4zTQzHKjpodqKiBpiP7X0PtKPDsxsw8tUjWx2cweJfJFeHQtY41nW8VA60pFWhP5sFfWE+i69bUNXt+ridTzb93Xv4CVwEBJPwb2A56qRaw7VL6+SOoDHAbcHWNxSfB3vJktNbOVwF3E+P9YpMfev4BuwB+r2N0Q4E0z2/r5exI4MTiD3VEDzSwXyCLSEeOfkjrXchunAv8ys4X1EE9MwbE+TqSttKp235hx1OKzsgrooKh2dDP7WfD6rGLbGrqtn/OeZnaBRc72t8Za03fMdpKl4T66J04PIr98Yqn85d1DUloNyao+bARyop53BhZXE1uVJPUEHiJSbfS+mZVLmgNsbUBfCuRLUlSy6k7kdHtHtx0KM9tapbNDm6H+jiN6W18AaZJ2MbMvg3l7E7vq5XtgoZntUsP2HyNSfdcbmBn1BRyvHS1fHw4hUqX+nSJ9O1oCqZL+z8z2kbSYbd/zNb3/04jRRhVULx4DlEva+iMqA2hLpIPBP3bgGP4XnFk5MEXSg0TO8J6rRfHTgDH1EUcsirzAjxD5wXO0mZXuYBxVfVbeBzYTOYN9vg6h1vk7JlnOqC6U1E1SO+Aa4Jk4ynxE5Et9jKQWkrIkHdRA8c0BTpaUKulIIh+gumpB5I1UCCDpDKIaKom8mcqBixTplDGAyK/q+th2QgmOL4tI9VBq8D9Mk5QrqX/U8yFEqkRfre9tmdkGYAowKngfHUTkg7xdF2ki77n1inScyQ7eD3tI2rfSeo8RORs5h0hPvug4J6vmruFVlq9JVa9DTctimEgksfQJpgeAl4m0JQNMAi6W1FGRDh+XEWkzJph3oqSWwWvUn0gbbuUOGQCDiJyh/V/Uvv6PyP9nSBzHG8/riSIGEEmAn9e0flS5nxFpj6zcy64+X+v7gZ8Q6WlXEmuFquKozWfFzIqItClNkHS8pFaSUoKz5xY1vxpAXb9j6lJ3m0gTkfrQEUQu8isi8sHMiVp2WOXHUWV7EOkptYpIdcm4Sts9LOr5SII69iqWPwGMjHp+NvBG8LgfkV/Y64l8gT3N9m1UlWOrcv/ALUS64W6tMvkncHbUuv2IJMdiIm/MKUR6pcVzXFVuO46yD1B9e9ciIl8qxVHT1HjKVrG9kfyvR9jWaSSRziz/Dl7vIiIN4IdXKvsKcHU9basdkffRBuA74ORqYu4a/P+XAWuC7R0WY723guWZlea/SRztetWUr+l/FPN1qGlZrNc0xnaj3yvpRNouioLXYhz/a1fLC953RcA64LOqjhl4HbgzxvyjYx1/jPWqfD0rvV/XA/PYtl26xvcs8CDweEO91kSqkw3YxLafqyHxxEEc7+8YZbb2sNxIJOF8CJwLZMT6nohRvtrvr1iT30KpGVDkQscHzGxS2LG4ulOkd9tcYC+runrHxclfz6bDE1USkvRLYAGRXyxDiPzy+5GZLQ01MOecq4Nk6UzhttWbyPU0LYhc7X+8JynnXFPlZ1TOOecSWrL0+nPOOZekPFE555xLaN5GBXTo0MF69eoVdhjOOdekfPzxxyvNLK+h9+OJCujVqxezZs0KOwznnGtSJH3bGPvxqj/nnHMJzROVc865hOaJyjnnXELzROWccy6heWcK55xrhqbNLmDszAUsKSqha242w/v3ZmDf/LDDiskTlXPONTPTZhcwYspnlJSWA1BQVMKIKZ8BJGSy8qo/55xrZsbOXPBDktqqpLScsTMXhBRR9TxROedcM7OkKOb4ilXOD5snKueca2a65mbXan7YPFE551wzM7x/bzKtbJt52empDO/fO6SIqueJyjnnmpn+aWu4ZNYzdE4tQ0B+bjajB+2ZkB0pwHv9Oedcs1M4fhyHF3/D+SN+TWrLlmGHUyM/o3LOuWZk4+zZbHj7HdqddVaTSFLgico555qVwnHjSG3XjnZDTg47lLiFkqgkDZY0X1KFpH5R8/eTNCeY5ko6torykyUtjFq3TzBfksZJ+krSp5L2aaRDcs65hLfhw4/Y+P4HtD/nHFJatAg7nLiF1UY1DxgEPBhjfj8zK5PUBZgrabpZpe4pEcPN7LlK844Cdgmm/YH7g7/OOdesmRmF995LWseOtD3pxLDDqZVQzqjM7HMz2+4SaDPbGJWUsgCr5aYHAI9ZxAdAbpDwnHOuWdvwr39R8skndPjj+aRkZYUdTq0kXBuVpP0lzQc+A86v4mwK4Jageu9uSZnBvHzg+6h1FgfznHOu2TIzCu+5l/T8fHKPOy7scGqtwRKVpDckzYsxDaiunJl9aGa7A/sCIyTFSv0jgN2CddoBV9YhvnMlzZI0q7CwsLbFnXOuyVj/xhtsmj+fDhdeiDIywg6n1hqsjcrMDtvB8p9LKgb2AGZVWrY0eLhZ0iTg8uB5AdA9atVuwbxY258ITATo169fbasYnXOuSbCKClaOG0/GTjvR5nfHhB1OnSRU1Z+knSSlBY97EjlrWhRjvS7BXwEDiXTCAHgROC3o/XcAsDYqqTnnXLOz7uUZbP7yS/IuvgilNc17PITVPf1YSYuBA4GXJc0MFv2cSE+/OcBU4AIzWxmUmSGpa7Dek5I+I9KO1QG4OZg/A/gG+Ap4CLigMY7HOecSkZWWUjh+PJm77UarI48MO5w6CyW9mtlUIomo8vzHgcerKHN01ONDq1jHgAvrKUznnGvSip6fQul339HtgftRSkJVoNVK043cOedclSo2bWLlhAlk9+1Ly1/+MuxwdognKuecS0JrnnqashUryLvsUiLN+U2XJyrnnEsy5cXFrJo4kRYHHUSL/fYLO5wd5onKOeeSzOpJkykvKiLv0kvDDqVeeKJyzrkkUrZqFasnTaJV//5k77lH2OHUC09UzjmXRFY++CAVmzeT96c/hR1KvfFE5ZxzSaK0oICip/9Km2MHkvmjncIOp954onLOuSRReN+fQSLvwuS6nNQTlXPOJYHNX33F2hdeoO2QIaR3Sa7RjTxROedcElhx9z2k5OTQ/txzwg6l3nmics65Jm7jJ7MpfvNN2p99Fmlt24YdTr3zROWcc02YmbHizjtJzetAu9NOCzucBuGJyjnnmrDit96i5OOPybvwQlJycsIOp0F4onLOuSbKysspvOtuMnr2bJJDzMfLE5VzzjVRa1+cHhkU8dI/ofT0sMNpMJ6onHOuCarYtInCe+8la889m/SgiPHwROWcc03QmieeoGzZMjpefnmTH8ajJp6onHOuiSkvKmLlxIdo8ctf0GL/pj+MR008UTnnXBOz8sGJVKxfT8ehw8IOpVF4onLOuSaktKCANU88QZtjjyWr965hh9MoQklUkgZLmi+pQlK/qPn7SZoTTHMlHVtF+Xei1lsiaVow/xBJa6OWXd9Ih+Scc41ixT33QkoKeRdfFHYojSYtpP3OAwYBD8aY38/MyiR1AeZKmm5mZdErmdnBWx9Leh54IWrxO2b22waK2znnQlMybz7rpk+n/TnnJN2NZ6sTSqIys8+B7XqqmNnGqKdZgFW3HUmtgUOBM+o5ROecSyhmxoqxY0lt2zYpbzxbnYRro5K0v6T5wGfA+ZXPpioZCLxpZuui5h0YVBu+Imn3avZzrqRZkmYVFhbWT/DOOddANrz9Nhs//JAOF1xAaqtWYYfTqBosUUl6Q9K8GNOA6sqZ2YdmtjuwLzBCUlY1q58EPB31/BOgp5ntDYwHplWzn4lm1s/M+uXl5cV9XM4519isrIwVd9xBes8etP39CWGH0+garOrPzA7bwfKfSyoG9gBmVV4uqQOwH3BsVJl1UY9nSJogqYOZrdyRWJxzLkxFU6ey+cuvyL/3XpSREXY4jS6hqv4k7SQpLXjcE9gNWFTF6scDL5nZpqjynRU0fEnaj8jxrWrQoJ1zrgFVbNhA4bhxZPftS6sjDg87nFCE1T39WEmLgQOBlyXNDBb9nEhPvznAVOCCrWdDkmZI6hq1mRPZttoPIslrnqS5wDjgRDOrtkOGc84lslWPPEJ54Uo6XXlF0t8qqSry73Ho16+fzZq1Xe2ic86FqnTZMr4+8ihaHXoo+XfdGXY425H0sZn1q3nNHZNQVX/OOef+p/Due6CigryhQ8MOJVSeqJxzLgGVzJ/P2hdeoN1pp5LRLT/scELlico55xKMmbFizG2Ri3vPOy/scELnico55xLM+jfeYOO//02Hiy9qdhf3xuKJyjnnEoht2cKKsXeQsfOPaXtC87u4N5awbkrrnHMuhtVPPkXpd9/R/aGJKM2/osHPqJxzLmGUrVnDygkTaHHwwbQ8+OCaCzQTnqiccy5BrBx/HxUbN9LpyivCDiWheKJyzrkEsPnLL1nzzDO0/f0JZO68c9jhJBRPVM45FzIzY/noMaS0aEGHiy8OO5yE4y11zjkXkmmzCxg7cwFLikrIa3kIl5xyOL3btg07rITjZ1TOOReCabMLGDHlMwqKSjBgRU47bilsw7TZBWGHlnA8UTnnXAjGzlxASWn5NvNKSisYO3NBSBElLk9UzjkXgiVFJbWa35x5onLOuRB0zc2u1fzmzBNVHU2bXcBBY/7OTle9zEFj/u71ys65WvnT/2WTWbZlm3nZ6akM7987pIgSlyeqOqjcCFpQVMKIKZ95snLOxcUqKtjvyXu59KtX6NomEwH5udmMHrQnA/s27yE9YvHu6XUQuxG0nLEzF/ibzDlXo3XTp1Mydy4n3Xorfxx0WNjhJDw/o6oDbwR1ztVVeXExy++4g6y99qLNwAFhh9MkhJaoJA2WNF9ShaR+MZb3kFQs6fIqyu8k6UNJX0l6RlJGMD8zeP5VsLxXfcfujaDOubpaOeF+ygtX0vnaa1CKnyvEI8xXaR4wCHi7iuV3Aa9UU/424G4z2xlYA5wVzD8LWBPMvztYr14N79+b7PTUbeZ5I6hzriabv/6a1Y89RpvjjyN7r73CDqfJCC1RmdnnZhbzyjZJA4GFwPwqlgs4FHgumPUoMDB4PCB4TrD818H69WZg33xGD9qT/NxshNFx42pu3K+tt08556pkZiy/5RZSsrPpeNllYYfTpCRcZwpJLYErgcOBmNV+QHugyMzKgueLga1ZIh/4HsDMyiStDdZfWWk/5wLnAvTo0aPWcQ7sm8/AvvmUr1/PV4f+mhYcCMfsX+vtOOeah/Wvv86G996n0zXXkNa+fdjhNCkNekYl6Q1J82JM1bUgjiRSpVfckLGZ2UQz62dm/fLy8uq8ndRWrWg75GTWv/46m7/+uh4jdM4li4qSElaMuY3MXXel7Uknhh1Ok9OgZ1RmVpd+l/sDx0u6HcgFKiRtMrP7otZZBeRKSgvOqroBWy9iKgC6A4slpQFtgvUbTLvTTmP15EdZ9dDDdB0zuiF35ZxrglY++CClS5bQ8/HHfHj5Oki4LidmdrCZ9TKzXsA9wK2VkhRmZsA/gOODWacDLwSPXwyeEyz/e7B+g0lr147cEwazdvp0tiz2i36dc/+zeeFCVj/yF1r/7hhy9t037HCapCoTlaR2cUy5dd2xpGMlLQYOBF6WNDOOMjMkdQ2eXgkMlfQVkTaoR4L5jwDtg/lDgavqGmNttD/zTEhJYdXDDzXG7pxzTYCZsfzmW1BmJp2GDw87nCarunPQJcFUXY+5VKD2PREAM5sKTK1hnZGVnh8d9fgbYL8YZTYBg+sS045I79yZ3GOPZe3zU+jwxz+S3qlTY4fgnEsw6197nQ3vvkunq68mbQfawpu76qr+PjezH5nZTlVNNHDbT1PT/txzsIoKVj3ySM0rO+eSWsWGDSwfM4bM3r1pe/JJYYfTpFWXqA6Io/yB9RVIMsjo1o02v/sdRc/8jbKVK2su4JxLWoUTJlC2dCmdb7jeO1DsoOoS1RRJpwTXNcUUVLO5KB3OOxcrLWXVpElhh+KcC8mmL75g9aOP0ea4QeTss0/Y4TR51SWqicBvgG8k/S3o/JDRSHE1WRm9etH66KNZ8/RfKVuzJuxwnHONzMxYNmoUqS1a0PHyqu5Z4GqjykRlZi+Y2UlAL+B54DTgO0mTJB3eSPE1SR3OPw8rKWH15EdrXtk5l1TWTnuBklkfk3f5MNLatg07nKRQ43VUZrbRzJ4xs2OBI4A+wKsNHVhTlrnzzrQ+6kjWPP64n1U514yUFxWxYuxYsvv0Ife448IOJ2nUmKgkdZJ0saR3gWnATMArXWvQ4Y9/pMLPqpxrVlbceSfla9fSeeQNPoRHParugt9zJP0d+ATYBRgedFe/yszmNlqETVTmLrvQ6sj+flblXDOx8eOPKXr2OdqdfjpZu+0WdjhJpbqUfyAwGuhuZpeY2XuNFFPS8LMq55oH27KFpTfcQFrXLuRddGHY4SSd6jpTnGlmrwMWdFO/Hn4YeXe7O0K47WXtuqufVTnXDKyaNJktX31N5+uuIyUnJ+xwkk48lagTiJxdbb20ej3w5waLKMnkXXBB5Kxq0uSwQ3HONYAt333HygkTaHX44bT61a/CDicpxZOo9jezC4FNAGa2BvDrqeKUucsutD7qKFY/8QRlq1eHHY5zrh6ZGctG3ojS0uh07TVhh5O04klUpZJSAQOQlAdUNGhUSabDRRdhmzax6qGHww7FOVeP1k2fzob33iNv2FC/EXUDiidRjSNyl/OOkm4B/gXc2qBRJZnMH+1Em2OOYc1TT1G6YkXY4Tjn6kHZmjUsHz2G7D59aHuij9rbkOK54PdJ4AoiPQCXAgPN7NmGDizZdLjwAqysjFUTfbwq55LBittup3z9ejqPutGvmWpgcd3S18z+C/y3gWNJahk9epA76FiKnnmG9medSXqXLmGH5Jyrow3vvcfaadNof955ZO26a9jhJL3qLvj9pKbC8azj/qfD+edjwMoJ94cdinOujio2bmTp9TeQ0asXHS74Y9jhNAvVnVH9RNKn1SwX0Kae40lq6fn5tD3hBNb89a+0P+tMMnr1Cjsk51wtFY4bT+nixfR8/DFSMjPDDqdZqC5RxXMPkPL6CqS56HD+eRRNmULh+PvIv/OOsMNxztVCyaefsvqxx8g98ffk7Ltv2OE0G1UmKjP7tjEDaS7S8vJod+qprJo4kfbnnO33BHOuibAtW1h6zbWkdehAx2HDwg6nWQmlq4qkwZLmS6qQ1C/G8h6SiiXFHHVM0pOSFkiaJ+kvktKD+YdIWitpTjBd39DHUhftzzqTlNatKbzn3rBDcc7FaeVDD7H5yy/pPPIGUlu1CjucZiWsPpXzgEHA21Usvwt4pZryTxKpmtwTyAbOjlr2jpn1CaZR9RFsfUtt04b2Z51F8VtvsfET74/iXKLbtOALVj7wIK1/8xtaHXpo2OE0O/GMR3WxpHodptLMPjezBVXsbyCwEJhfTfkZFgA+ArrVZ3yNod2pp5Ca14EVd91F5DCcc4nIyspYes01pLZq5bdJCkk8Z1SdgH9L+pukIyWpoYKR1BK4ErgxzvXTgVPZdsThAyXNlfSKpN2rKXuupFmSZhUWFu5Q3HWRkpND3oUXUjLrY4rfeqvR9++ci8/qyZPZNG8ena+71oeWD0k8d6a4lsjAiY8AfwC+lHSrpB9XV07SG0EbUuVpQDXFRgJ3m1lxnPFPAN42s3eC558APc1sb2A8kRGJqzquiWbWz8z65eXlxbm7+pV73HFk9OxJ4V13YeXegdK5RLP5m4UUjhsfuTP6kUeGHU6zFVcbVVDFtiyYyoC2wHOSbq+mzGFmtkeM6YVqdrU/cLukRcClwNWSLoq1oqQbgDxgaNQ+121NcmY2A0iX1CGeYwyD0tPJu+wyNn/5FWtfeDHscJxzUay8nKVXX42ys+l8/XU0YGWSq0GNt1CS9CfgNGAl8DCRIelLJaUAXxK5D2C9MLODo/Y7Eig2s/tixHQ20B/4tZlVRM3vDCw3MwsGd0wBVtVXfA2hVf8jyNprLwrHjaP10UeRkpUVdkjOOWD15EcpmTOHrmPHkhZSrYuLiOeMqh0wyMz6m9mzZlYKECSI39Zlp5KOlbSYyICML0uaGUeZGZK6Bk8fINJ29n6lbujHA/MkzSVy1/cTLcF7Kkii47BhlC1bxponnww7HOccsPmbbyi8915aHvZrWv/2N2GH0+yppu9xSe1izF6/NWElg379+tmsWbNCjeG7886j5JPZ/Pi1md5g61yIrLycb08ewpZFi/jRS9P9bKoakj42s+2uha1v8ZxRfQIUAl8QqeorBBZJ+kTSTxsyuOak0+WXU7FhAyvv9xvWOhem1ZMmUTJ3Lp2uu86TVIKIJ1G9DhxtZh3MrD1wFPAScAGRXneuHmTusgu5xx3HmqeeZsuiRWGH41yztGnBFxTeO45WRxxB698cHXY4LhBPojrAzH5oQzKz14ADzewDwG8dXI86XHwRyshgxV13hx2Kc83KtNkFHDT6TX7yly847ddX8dEJF3gvvwQST6JaKulKST2D6QpguaRUoKKmwi5+6R070v6sM1n/2mt+ayXnGsm02QWMmPIZBWs3YRIrstpwzcyFTJtdEHZoLhBPojqZyC2KpgFTge7BvFTghAaLrJlqf8YZpOXlsXzMbViF/w5wrqGNnbmAktJtL7gvKS1n7MyYd3lzIag2UQVnTfea2cVm1tfM9gkeF5rZFjP7qpHibDZScnLIu+wyNn36KetefjnscJxLekuKSmo13zW+ahOVmZUDPSVlNFI8DmgzcABZu+/OijvupGLjxrDDcS6pddLmmPO75mY3ciSuKvFU/X0DvCvpOklDt04NHVhzppQUOl09grLly1n1yF/CDse5pFX89tuc+tFzZFVqbs9OT2V4/94hReUqiydRfU2kO3oK0Cpqcg0o56c/pdVRR7LqkUcoXbo07HCcSzpla9aw5JprOCqnmNGD9iQ/NxsB+bnZjB60JwP75ocdogvUeGeKH1aUcswsKeuhEuHOFLFsWVzAN0cfTasjjiD/jrFhh+Nc0jAzCi65hPVv/ZOdnv0bWbvtFnZITVLC3JlC0oGS/gP8N3i+tyS/0LcRZHTLp92ZZ7DupZe8u7pz9Wjt88+z/vU36HjppZ6kmoB4qv7uIXKn8lUAZjYX+EUDxuSidDj3XNI6d2bZTTf7mFXO1YPNCxey7JZbyTngANqd8Yeww3FxiHc8qu8rzfJvzEaSkpNDpyuGs/nzzyl69rmww3GuSbMtW1hy+XCUkUHX28aglLi+Al3I4vkvfS/pZ4BJSpd0OfB5A8florQ66ihy9tuPwnvuobyoKOxwnGuyCsffx6b58+kyahTpnTqFHY6LUzyJ6nzgQiAfKAD6BM9dI5FEp2uuoXz9egrHjQs7HOeapA3vv8+qhx8md/DxtO5/RNjhuFqoMVGZ2UozG2Jmncyso5mdYmYJPWpuMsrqvSttTzqJNX99hpL588MOx7kmpWz1apZccSUZO+1EpxEjwg7H1VI8vf7yJF0taaKkv2ydGiM4t628P11Cart2LLtxlN8H0Lk4mRlLR1xNeVER+XfdSUpOTtghuVqKp+rvBaAN8AbwctTkGllqq1Z0umI4mz791DtWOBenNY8/QfE//0nHK67wruhNVFoc6+SY2ZUNHomLS+tjjuHZlz5i4nulFM5+ma652Qzv39uvoncuhpJ581kxdiwtDzmEtqcMCTscV0fxnFG9JKleh7qUNFjSfEkVkra7qllSD0nFQQ/DWOUnS1ooaU4w9QnmS9I4SV9J+lTSPvUZdyJ4Yc4S7sj7GSuyczGgoKiEEVM+87FznKukfP16Ci67jNQOHegy+lYfCLEJiydR/YlIstokaZ2k9ZLW7eB+5wGDgLerWH4X8EoN2xhuZn2CaU4w7yhgl2A6F7h/B+NMOGNnLmBT+ba3vfKxc5zblpmx9PrrKV2yhPw77yCtbduwQ3I7oMaqPzOr9xvQmtnnQMxfOJIGAguBDXXY9ADgMYvcwPADSbmSuphZ0tzV1cfOca5mRc88w/pXXiVv6FBy9km6ipVmJ55ef5J0iqTrgufdJe3XEMFIaglcCdwYx+q3BNV7d0vKDOblA9F30VgczIu1r3MlzZI0q7CwcIfibkxVjZHjY+c4F1Eyfz7Lbx1Ni5//nPZnnxV2OK4exFP1NwE4kMjw8wDFwJ9rKiTpDUnzYkwDqik2ErjbzIpr2PwIYDdgX6AdkeRWK2Y20cz6mVm/vLy82hYPzfD+vclOT91mXmbZFi7bx6s2nCtft46CSy8jtW1but5+m98iKUnE0+tvfzPbR9JsADNbE8+Iv2Z2WB3i2R84XtLtQC5QIWmTmd1Xadtbq/I2S5oEbO10UQB0j1q1WzAvaWzt3Td25gKWFJXQpVUGp77/PPsWVmCHT/IGY9dsmRlLrr6a0qVL6fnYY6S1axd2SK6exJOoSiWlAgaRC4CBBrna1MwO3vpY0kiguHKSCpZ1MbOlinwrDyTSOQPgReAiSX8lkvTWJlP71FYD++Zv0x19zTNrWHbDDax9/nlyjz8+xMicC8/qyY9S/MabdLzySnL26Rt2OK4exXNePA6YCnSUdAvwL+DWHdmppGMlLSZSpfiypJlxlJkhqWvw9ElJnwGfAR2Am4P5M4BvgK+Ah4ALdiTOpiJ38PHk7Lsvy2+7ndLlK8IOx7lGt3HWLFbccQetDj+Mdn84PexwXD2La4RfSbsBvwYEvLm1116ySNQRfmtjy6JFfDNgIC1+/nO63TfeqwBds1G6YgULjzuO1JwW9HruWVJb1XtHZVeFhBnhF8DM/mtmfzaz+5ItSSWLjF69yLvkYorffJP1r74adjjONQorLaXg0suoKN5A/vhxnqSSlHeJSSLtTj+drN13Z9lNN1O2Zk3Y4TjX4JaPHUvJJ5/Q5aabyNp117DDcQ3EE1USUVoaXW69hfL161l+001hh+Ncg1r7wguseexx2p56Km1++5uww3ENyBNVksnq3Zu8Cy9g3YxXWOdVgC5Jlcyfz9LrbyBn333pdMXwsMNxDcwTVRJqf/bZkSrAG0dRtsrHuHTJpWz1ahZffDGp7dqRf8/dKD097JBcA/NElYSUlkbXMaOpKC5m2cgbiadnp3NNgZWWUnDZUMpXrabb+PGktW8fdkiuEXiiSlKZu+xC3p8uYf3rr7Nu+vSww3GuXiwfPYaNH35Il1E3kr3H7mGH4xqJJ6ok1u6MM8jeZx+WjbqJ0oKkupOUa4bW/PUZ1jz1FO3OPJM2A6q7ZahLNp6okphSU+l6+21QUcGSq0ZgFQ1y5yvnGtyGjz5i2c030+KXv6DjsKFhh+MamSeqJJfRrRudrrmGjf/+N6snTQ47HOdqbcv331Pwp0vJ6NGD/DvuQKmpNRdyScUTVTPQZtCxtDr8MArvuYdNn/uNRVzTUb5uHd+f/0eoqKD7hD/7nSeaKU9UzYAkOo8aRWpuLgXDLqdi48awQ3KuRltvj7Tlu+/IHz+OjF69wg7JhcQTVTORFgwkt2XhQpaPHh12OM5Vy8xYdsstbHjvPbqMHEmL/RpkUHHXRHiiakZaHHgg7c8+m6Jnn/O7VriEtnrSZIr++gztzz6L3OMGhR2OC1k8Aye6JJJ3ycVs+PBDll53PVl77ElGt/yaCznXCKbNLvhh5Oq8jWmcf/TpnDHUe/g5P6NqdpSeTv6dd4AZBUOHYlu2hB2Sc0ybXcCIKZ9RUFSCASty2jG2xd68MDfpBuh2deCJqhnK6N6dLrfcwqZPP2XFnXeGHY5zjJ25gJLS8m3mlZRWMHbmgpAiconEE1Uz1br/EbQ99VRWP/oY6157LexwXDO3pKikVvNd8+KJqhnrNPxysvbck6XXXMuW774LOxzXTJUXbyCvdH3MZV1zsxs5GpeIQklUkgZLmi+pQlK/GMt7SCqWdHkV5d+RNCeYlkiaFsw/RNLaqGXXN/ChNGnKyCD/7rshJYXFl/yJihL/9eoal23ZQsElF3P6py+RVenbKDs9leH9e4cTmEsoYZ1RzQMGAW9Xsfwu4JWqCpvZwWbWx8z6AO8DU6IWv7N1mZmNqq+Ak1VGt3zy7xjL5gULWHrDDT4kiGs0Vl7OkquuYsN77zPk/EGMGdyH/NxsBOTnZjN60J4M7Ou9Ul1I3dPN7HOI3DGhMkkDgYXAhpq2I6k1cChwRv1G2Ly0PPhgOlx8ESvHjSd7r71pd8qQsENySc7MWHbTTayb8QodLx9G7nGDGAiemFxMCdVGJaklcCVwY5xFBgJvmtm6qHkHSpor6RVJPmBNnDqcfz4tf/Urlo8Zw8aPPw47HJfkCu+994cLetuffXbY4bgE12CJStIbkubFmKobSGYkcLeZFce5m5OAp6OefwL0NLO9gfHAtGriO1fSLEmzCgsL49xd8lJKCl1vG0NGfj6LL/kTpUuWhB2SS1Kr/jKJVQ88SO7g48kbNizscFwToDDbJCS9BVxuZrOC5+8A3YPFuUAFcL2Z3RejbAdgAZBvZpuq2P4ioJ+Zrawujn79+tmsWbPqeBTJZfM337Do9yeS3q0bvZ58gpScnLBDcklkzdNPs+zGUbQ66kgfsiMJSPrYzLbrEFffEqrqL+gk0cvMegH3ALfGSlKB44GXopOUpM4KGr4k7Ufk+FY1bNTJJfNHPyL/zjvYvGABS0Zc7YMtunpTNHUay24cRctf/Yr822/3JOXiFlb39GMlLQYOBF6WNDOOMjMkdY2adSLbVvtBJHnNkzQXGAecaN6NrdZa/uIXdLz8ctbPnMnK+/4cdjguCaybMYOl11xDi58dSP49d6P09LBDck1IqFV/icKr/rZnZiy95lrWTplClzGjyR04MOyQXBO17tVXKRh2Odl9+9Bj4kSvTk4ijVX153dPdzFJosvIGyhdsoSl111PepeutNjfxwRytbNu5muRJLX33vR48EFPUq5OEqqNyiUWZWTQbdy9ZPToweKLL2bzN9+EHZJrQta9/joFw4aRvddedJ84kZQWLcIOyTVRnqhctVJbt6b7gw+g9HS+P+dcSlesCDsk1wSsmzGDgksvI3v33en+0ERSW3qScnXnicrVKKNbN7o/8ACvZXbjoNv+zk5XvcxBY/7OtNkFYYfmEtDaF1+k4PLhZPftQ/dHHiG1ZcuwQ3JNnLdRubjMLGvLuL6D2RQMGVRQVMKIKZ8Bftsb9z9Fzz3H0uuuJ2f//ek+4c/eJuXqhZ9RubiMnbnghyS1VUlpuQ9s536wavJkll57HS0OOojuD9zvScrVGz+jcnHxge1cVcyMlePvY+WECbTq35/8sbejjIyww3JJxM+oXFyqGsCukzb70CDNmFVUsPyWW1k5YQJtBg0i/847PEm5eueJysVleP/eZKdve8ubLCo49aPnWHHHHZ6smqGKLVsoGDaMNU88QbvTT6fLzTehNK+kcfXP31UuLls7TIyduYAlRSV0zc3m8iN25YDsT1j9yF9QSip5Qy+LOcaYSz7lxcUsvuhiNn7wAR2HD6f9WWeGHZJLYp6oXNwG9s3froef9bkOzFj10EMAnqyagdJly/j+vPPZ/PXXdL1tDG0GVDdyj3M7zhOV2yFKSaHz9dcDRJKVVZA3bJgnqyS16fPP+f6886nYsIHuDzxAy58fFHZIrhnwROV22A/JSmLVw49QUbKJTtdcjVK8CTSZFP/znxRcNpSUNm3o+dRTZPXeNeyQXDPhicrVi63JKiUrm9WTJlGxcSNdbhrljetJwMxYPflRVowdS9Zuu9Ht/vtJ79Qx7LBcM+LfIq7eSKLjFcNJadmClePvo2LDBrqOvZ2UzMywQ3N1VLFlC8tuGMnaqVNpdcQRdB0z2i/kdY3OE5WrV5LIu/BCUlu2ZPnoMXy/Zg3d/nwfqa1bhx2aq6XS5SsouPRSSmbPpsOFF9Lhwgu8OteFwgdOxAdObChrp7/EkquvJrNXL7o/NJH0zp3DDsnFYdrsAm6b/inLNpSTt2ktl+3TlpNOPzrssFwCaqyBE/3nkWswbY75LT0mPkjpkiUsOvEkNn3+edghuRpM/WQxV/1tNks3VmASK7JzGfVVit8p34XKE5VrUC0OPJCeTz4BwKIhp7D+zTdDjshVpXz9esb89QM22baXFvjNh13YPFG5Bpe12270+tszZO68M4svupiVDz3kt1xKMCXz5rNw0HGsUFbM5X7zYRemUBKVpMGS5kuqkNQvan4vSSWS5gTTA1WUbyfpdUlfBn/bBvMlaZykryR9KmmfxjomV730jh3p+dijtD7qSArvvIuCoUOp2LAh7LCaPauoYNUjf2HRSSdhZWV0bhG7f1VVNyV2rjGEdUY1DxgEvB1j2ddm1ieYzq+i/FXAm2a2C/Bm8BzgKGCXYDoXuL9+w3Y7IiUri6533knesKGsn/kai048kc0LF4YdVrNVumwZ3515FivGjqXVIb9kpynPc+Uxe2538+Hs9FSG9+8dUpTOhZSozOxzM9uRSu8BwKPB40eBgVHzH7OID4BcSV12YD+unkmiwznn0OPhhygrXMmiwSew7pVXwg6rWTEz1k6fzjcDBlIydy5dbr6J/HHjSGvbloF98xk9aE/yc7MRkJ+bzehBe/oozi5UiXgd1U6SZgPrgGvN7J0Y63Qys6XB42VAp+BxPvB91HqLg3lLqUTSuUTOuujRo0c9he7i1eJnP2On55+jYOgwCi4byob33qfT1SNIyfYqpoZUtnIlS0eOpPiNN8nu04euY0aT0avXNuvEuvmwc2FqsDMqSW9Imhdjqu5Wy0uBHmbWFxgKPCWp2itFLdIqX+uWeTObaGb9zKxfXl5ebYu7epCen0/PJx6n/TnnUPTccywcPNi7sDcQM6No6jS++c1v2fD2O3QcPpyeTz6xXZJyLhE12BmVmR1WhzKbgc3B448lfQ3sClS+Gne5pC5mtjSo2lsRzC8Auket1y2Y5xKU0tPpOGwoOQfsz9KrRrBw8AnkXXgB7c85x+8TWE+2LFrE0pE3svGDD8ju25cut9xM5o9+FHZYzsUtobqnS8qTlBo8/hGRThHfxFj1ReD04PHpwAtR808Lev8dAKyNqiJ0CazlQQfxo+kv0rp/fwrvHceik4ew6Ysvflg+bXYBB435Oztd9TIHjfm7X4Aah4qSEgrHjeOb3w1g0/z5dB45kp5PPuFJyjU5odxCSdKxwHggDygC5phZf0nHAaOAUqACuMHMpgdlHgYeMLNZktoDfwN6AN8CJ5jZakUGQboPOBLYCJxhZjXeG8lvoZRY1r36KstuHEX5+vW0P/NM3v35QK6evoCS0vIf1slOT/VG/iqYGetnvsby22+jbMlSWv/mN3S88grSO/odz139aqxbKPm9/vBElYjK1qxhxe1jWTt1KqcfdT0rMrdvqszPzebdqw4NIbrEVTJ3LsvHjqVk1sdk9u5N52uvIWfffcMOyyWpxkpU3gjgElJa27Z0HX0rbQYMoHBqYcx1/G4J/7N54UIKx41j/SuvktqhA51H3kDu8cd7O59LCv4udgmtxQH70/Ufb1KwdtN2y5L9bgnTZhcwduYClhSV0DU3m+H9e29X1bnl229ZOeF+1k6fjjIz6XDBBbQ780xSW7YIKWrn6p8nKpfwhh+5GyOmfLZNG1Vm2RbOKyugbE1f0tq2DTG6hjFtdsE2x1xQVMKIKZ8BkeucNi34glWPPMy6l2egtDTanX467c8+i7T27cMM27kG4W1UeBtVUxB9dtGlZTrnrJvLAdMnoaws2p4wmHZnnJFU410dNObvFMSo2uySncIzS6ax4Z9vo5wc2g4eHElQfi2gC4F3pmhEnqiaps1ffcWqhx5m7UsvQUoKbY4+mrannEL2nnuEHdoO2+mql2NexS4zZr59K21PPYV2J59Mam5uY4fm3A+8M4VzNcjceWe63jaGDhdfzOpJk1g7dSprX3iBrL33ou3vT6R1/yNIadE022q65mZRULR9u1znTNj5rX+QkpkZQlTOhcPPqPAzqmRRXlzM2mkvsObJJ9mycCHKzqb1EYfT+rfH0OKA/VF6eq23GU+HhvpiZmz67DPWvfIqU/+9iLt7Hc7mtIwflmenpzB60F5+7ZhLGF7114g8USUXM6Nk9hzWTpvGuhkzqCguJqV1a1r96le0PPRQWhx4AKmtq72FJLB9hwao/wuNKzZsYMOHH1L8z7cpfvttypYuhfR0Wv7sZ7yz32+4b0U2S9dubvAk6VxdeKJqRJ6oklfF5s1sePdd1r/2Ouv/8Q8q1q6FlBSy99yTnAMPILtPH7L33jtmz8GqOjTsyIXGZStXUjJvHiUff8zGj/5Nyfz5UFZGSosWtPjZz2h56KG0+vWhcSVS58LmbVTO1YOUzExaHXoorQ49FCstpeTTT9nw7rsUv/suqx56GMojZ0vpPXuQteuuZO6yCxk//jEZ3bpVeUFxTRcaW1kZZYWFlC5dxpZFi9iy8Bs2f/0Nm/7zH8qWLYuslJZG9h570P7MM2lx4AHk/PSnKCOj2u0611z5GRV+RtVcVWzcSMln8yiZM4dN8+ez+csv2fLtt1BRAcDpR1zNipx225XrVFHClLL3sfJyrLwMKy2lYn0xFevXUV60lrJVq37YBgDp6WT07EHWbj8ha/fdyd5jd7J2352UnJzGOlTnGoSfUTnXwFJycmix/3602H+/H+ZVbNrElm+/o7SggEvmr+Dmggo2RQ0ykFlRxlkF77Jp5eeQloZSU1F6OimtWpHRaydS2rQmvWNH0jp3Jr1zZzJ69iQ9P99vZeTcDvBPj3NRUrKyyOq9K1m9d+WUQ6FlzF5/1Y396Zyrb56onKuGD8vuXPgSauBE55xzrjJPVM455xKaJyrnnHMJzROVc865hOaJyjnnXELzROWccy6heaJyzjmX0PwWSoCkQuDbOFbtAKxs4HDCkIzHlYzHBMl5XMl4TJCcx1X5mHqaWYMPL+2JqhYkzWqM+1o1tmQ8rmQ8JkjO40rGY4LkPK6wjsmr/pxzziU0T1TOOecSmieq2pkYdgANJBmPKxmPCZLzuJLxmCA5jyuUY/I2KueccwnNz6icc84lNE9UzjnnEponqjhJ+pOkeZLmS7o07HjqQtJfJK2QNC9qXjtJr0v6MvjbNswY66KK4xoc/K8qJDW5LsJVHNNYSf+V9KmkqZJyQwyxTqo4rpuCY5oj6TVJXcOMsbZiHVPUsmGSTFKHMGLbEVX8r0ZKKgj+V3MkHd0YsXiiioOkPYBzgP2AvYHfSto53KjqZDJwZKV5VwFvmtkuwJvB86ZmMtsf1zxgEPB2o0dTPyaz/TG9DuxhZnsBXwAjGjuoejCZ7Y9rrJntZWZ9gJeA6xs7qB00me2PCUndgSOA7xo7oHoymRjHBdxtZn2CaUZjBOKJKj4/AT40s41mVgb8k8iXYJNiZm8DqyvNHgA8Gjx+FBjYmDHVh1jHZWafm9mCkELaYVUc02vB+w/gA6Bbowe2g6o4rnVRT1sATaqHVxWfK4C7gStoYsezVTXH1eg8UcVnHnCwpPaScoCjge4hx1RfOpnZ0uDxMqBTmMG4uJ0JvBJ2EPVF0i2SvgeG0PTOqLYjaQBQYGZzw46lAVwUVNX+pbGaCjxRxcHMPgduA14DXgXmAOVhxtQQLHKtQpP89decSLoGKAOeDDuW+mJm15hZdyLHdFHY8eyI4Mfs1SRBwo3hfuDHQB9gKXBnY+zUE1WczOwRM/upmf0CWEOkjSAZLJfUBSD4uyLkeFw1JP0B+C0wxJLzIsgngePCDmIH/RjYCZgraRGRKtpPJHUONap6YGbLzazczCqAh4i02zc4T1RxktQx+NuDSPvUU+FGVG9eBE4PHp8OvBBiLK4ako4k0ubxOzPbGHY89UXSLlFPBwD/DSuW+mBmn5lZRzPrZWa9gMXAPma2LOTQdtjWH7WBY4k0izT8fpPzR1n9k/QO0B4oBYaa2Zshh1Rrkp4GDiFyq/7lwA3ANOBvQA8iQ52cYGYJ0YAaryqOazUwHsgDioA5ZtY/pBBrrYpjGgFkAquC1T4ws/NDCbCOqjiuo4HeQAWR9+D5ZlYQVoy1FeuYzOyRqOWLgH5m1qSG/Kjif3UIkWo/AxYB50W1cTdcLJ6onHPOJTKv+nPOOZfQPFE555xLaJ6onHPOJTRPVM455xKaJyrnnHMJzROVc865hOaJyrkkJ+lsSZ9JOiNq3k8kPSDpOUl/DDM+52riicq55HcccCgweOuM4O7y5wMnAAeFFZhz8fBE5VwDklTciPvqJalE0pxKiz4kcg/HDyut/zvgZWCGpOxgILwtTXGQP5fcPFE5l1y+DgYgjNYSeAdoEz3TzF40s6OI3OC2JCi3pFGidK4WPFE5FydJYyRdGPV8pKTLJQ2VNC+YLo1Rrlel4bwvlzQyatl/JU2W9IWkJyUdJuldSV9K2i9Y7xRJHwVnPQ9KSo0z5hQiN1E+DTh2azlJh0gaJ+lBoFFGaXWurjxRORe/Z4i06Wx1AvBv4Axgf+AA4BxJfWu53Z2JjOuzWzCdDPwcuBy4WtJPgN8DBwVnPeVEBhiMx6HAXDNbBMwNnmNmb5nZJWZ2npn9uZbxOteo0sIOwLmmwsxmS+ooqSuRu7KvIXIn6almtgFA0hTgYGB2LTa90Mw+C8rPB940M5P0GdAL+DXwU+DfkgCyiX/csCHA08Hjp4Pnr9ciNudC54nKudp5Fjge6EzkDEtxlClj29qLrErLN0c9roh6XkHkMyrgUTMbUZtAJWUTGd/p15JuD2JoJSnbzEpqsy3nwuRVf87VzjPAiUSS1bNEOikMlJQjqQWRweTeqVRmOdBRUntJmURG6K2NN4HjowbvbCepZxzlfge8YmY9gkH8egDTgWNquX/nQuVnVM7VgpnNl9QKKAgGjFsqaTLwUbDKw2Y2u1KZUkmjgnUKqOUItmb2H0nXAq8FnSNKgQuJDDJYnSHAY5XmTSXSpva32sTgXJh84ETnkoSkXsBLZrbHDmxjEU1wNFqX3Lzqz7nkUQ60iXHBb422XvALpBNpG3MuYfgZlXPOuYTmZ1TOOecSmicq55xzCc0TlXPOuYTmico551xC80TlnHMuoXmics45l9A8UTnnnEtonqicc84ltP8H4rRIyzf2/AQAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "from scm.utils.conversions import plams_molecule_to_chemsys\n", "from scm.base import Units\n", "\n", "from ase.eos import EquationOfState\n", "import matplotlib.pyplot as plt\n", "\n", "res = job.results.get_pesscan_results(molecules=True)\n", "\n", "# Obtain volumes\n", "volumes = [plams_molecule_to_chemsys(x).lattice.get_volume() for x in res[\"Molecules\"]]\n", "\n", "# Obtain energies\n", "energy_unit_conversion = Units.conversion_factor(\"Ha\", \"eV\")\n", "energies = [energy * energy_unit_conversion for energy in res[\"PES\"]]\n", "\n", "# Fit and plot\n", "eos = EquationOfState(volumes, energies, eos=\"birchmurnaghan\")\n", "v0, e0, B = eos.fit()\n", "eos.plot(show=True);" ] }, { "cell_type": "code", "execution_count": 5, "id": "727ef08f-0a70-428f-a877-1c8849711ffe", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "e0= -15.150 eV\n", "v0= 11.463 ų\n", " B= 4.439 GPa\n" ] } ], "source": [ "print(f\"{e0=:8.3f} eV\\n{v0=:8.3f} Å\\u00b3\\n{ B=:8.3f} GPa\")" ] }, { "cell_type": "markdown", "id": "17da7e05-a3e4-4159-b8d9-1424e5cd5b4f", "metadata": {}, "source": [ "## Scale cell to match minimum volume\n", "\n", "`set_density` applies uniform scaling." ] }, { "cell_type": "code", "execution_count": 6, "id": "c2464d4f-04da-4c2e-91f5-85871366092b", "metadata": {}, "outputs": [], "source": [ "diamond_scaled = diamond.copy()\n", "diamond_scaled.set_density(diamond.total_mass() / v0)" ] }, { "cell_type": "code", "execution_count": 9, "id": "2a7d9d11-78ba-458d-b55b-9f6614171357", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "System\n", " Atoms\n", " C -0.4473943191070795 -0.4473943191070795 -0.4473943191070795\n", " C 0.4473943191070795 0.4473943191070795 0.4473943191070795\n", " End\n", " Lattice\n", " 0 1.7895772764283184 1.7895772764283184\n", " 1.7895772764283184 0 1.7895772764283184\n", " 1.7895772764283184 1.7895772764283184 0\n", " End\n", "End\n" ] } ], "source": [ "print(diamond_scaled)" ] }, { "cell_type": "code", "execution_count": 14, "id": "2a2b6f7a-f4ab-4912-8594-c5953c9e3ec9", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "diamond.lattice.get_volume()=11.375\n", "diamond_scaled.lattice.get_volume()=11.463\n" ] } ], "source": [ "print(f\"{diamond.lattice.get_volume()=:.3f}\")\n", "print(f\"{diamond_scaled.lattice.get_volume()=:.3f}\")" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.12" } }, "nbformat": 4, "nbformat_minor": 5 }