{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "*This notebook contains course material from [CBE20255](https://jckantor.github.io/CBE20255)\n", "by Jeffrey Kantor (jeff at nd.edu); the content is available [on Github](https://github.com/jckantor/CBE20255.git).\n", "The text is released under the [CC-BY-NC-ND-4.0 license](https://creativecommons.org/licenses/by-nc-nd/4.0/legalcode),\n", "and code is released under the [MIT license](https://opensource.org/licenses/MIT).*" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "< [West Virginia Chemical Spill](http://nbviewer.jupyter.org/github/jckantor/CBE20255/blob/master/notebooks/B.01-West-Virginia-Chemical-Spill.ipynb) | [Contents](toc.ipynb) |
"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Ajka Alumina Plant Spill"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Summary\n",
"\n",
"This notebook demonstrates the analysis of the 'Red Sludge' tragedy that took place in Ajka, Hungary, in 2010."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Background\n",
"\n",
"On Monday, October 4th, 2010, an earthen dam failed at an alumina ore processing facility near [Ajka, Hungary](http://www.theatlantic.com/infocus/2011/09/a-flood-of-red-sludge-one-year-later/100158/), resulting in the release of red mud over a large area in western Hungary. According to the [New York Times](http://www.nytimes.com/2010/10/06/world/europe/06hungary.html), a flood of 700,000 cubic meters of sludge swept cars off bridges, damaged roads, and flooded businesses and homes throughout the affected region. There were reports of up to eight deaths and hundreds of displaced families. The red mud entered creeks about 45km upstream of the Danube raising the potential for devastating environmental damage. At the time there were fears that other earthen dams at the same processing facility were showing signs of impending failure."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wCEAAYEBQkJCQkJCQkJCQkJCQcHCQkJCQgHBwcHDAgMDAwICgoOERkTDg8XEAoKFCAUFxobHR4dDRMtJCEcJBkcHRwBBwcHCwkLDQoKDRwODA4cHBwcHBwcHBwcHBwcHBwcHBwcHBwcHBwcHBwcHBwcHBwcHBwcHBwoHBwcHBwcHCgcHP/AABEIAWgB4AMBIgACEQEDEQH/xAAcAAABBAMBAAAAAAAAAAAAAAAEAgMFBgABBwj/xABSEAABAwIEAgUEDAsGBQQDAQACAAMSBCIBBTJCE1IGB2JykhEjgqIUM0NTk7LBwtHS0/AVFyExVIGRobHh8hZBUWHi8QhjcYOzJCVzozREwxj/xAAbAQADAQEBAQEAAAAAAAAAAAAAAQIDBAUGB//EACwRAQEAAgIBBQACAgIABwAAAAABAhEDEiEEExQxUQVSIkEGkSMyM0JhcdH/2gAMAwEAAhEDEQA/APPmPQyr98p/E59RZj0Nq/fKb4Rz6i6CWCRiuSc+Tf24oP8AY2r98pvG59RRz2SOtnwyNqWHaKP8F0zGSrGZ+2ld2k/eyOcc2iKDojUvlFtxj0ic+QFM0fVfmLxRF6j/AFuVHyNqa6LEIueiukdHzGQ9oS8S0487knl45jNxy0eprNsRlx8v+FqvsVrDqazX9Jy/4Wq+xXe2HSiIxS1tfDluVcC/Ezmv6Tl/wtV9is/Ezmv6Tl3w1V9iu/eVKCMtQqdl3rzdn/VlmGXtC+8/QkJFEeG7UEUv1t4Krlk7uBeSTXlw/wA3PoXo/rOZJyhGOkXZS7K4m6xIii4JeiXyI20wu/tWfwO7zN/tL6Fn4JPna8RfQp9yncErbkyTZd0ku1a9Yhzykx/ObXiL6E0GXmW4P2l9CmTAkjAS5SHtRtU96cwgBvJ3cfySD9pfQpeg6E1j+h2mw7xu4fwBFU1K5aQkXhV16L0xCJSKJeqp9zJdwx0rdD1S5nUEQg/Q4RGWMnagR/c0isepbNv0nLvhaz7Bda6NVrbAmRERWoj+0BHa0Mi3SERJaTKuS27cbx6ms1/Sst+FrPsFmHUzmv6Tl3w1V9iuxY1lTtGKQdXWR2j6Ke6c7VyHDqXzX9Jy74Ws+wWfiYzX9Jy34Wr+xXWCr6uKb/CNSJXEIqeyphnf9OVfiYzX9Jy34Ws+xTdV1PZq01i6VRl+Ij/cLtVL9mLK6zhmT/NKMkNUZk+Q8IpRldbtT7NMeLP/AHHDnOhVWJeTiU3je+om3uh9WGo2PE59RdXqWGhc7XrJt9tiJD82SO1Zz8rl39jKyQDNjzgyEpORj4FmHQ2pnw+PSy7731F1HBhshbEdvekmsMubJ1xy6JWiOkpD2f7kdqrw5tW9CqpgMXCfpCEY6DeIv3treHQis4XHwcpib8glrdl8RdNfo2hCLsiEubUg2aWJcNt3zRN6ZCjdGooR9Aq6LZ8WlLi4SHDBx2XxE010HriKPlY/WTkfiLq9FgLTQtlHiD2pRRJvgTgxIRjq2o3Rpymu6uswYHAnHKXyFym98rabPoBXCTY8SlxxdGYxcext8C6t0yIn2mhEo3alA4kYcIpSJq3mtS7Vcwlc/qOhlU2RCTlMRD+cRN3Eu9o/Mh3ei1SO5rTL8hOfUXVcGxKpbMtLolLs/Qt1NC04+2I3NiMXOb/dReSjpHHmMldcKIk3y3EX0KXb6DVhDIXKYu6b31FaKfL+EZRHzZEXejL6FZKKmINN0hEolyqPeyP23MMOhVVd5XqUI8xvXd3yNoWn6MPuFEXGPE59RdTzBwfOW7YqCySkGV1pCSPep+0rOHQKt99pPG99RZ/YKulHiU2HpvR+IupRISlq7Kbq61uBStLu/Qj36XtVzY+ruvEcT4tJEe299mhy6DVnPTF3Se+ouj/hhuJCJXFzDJOZJjIiIrhK3tCX38i1xzuRZYac3Y6BVzmk6bVGOJuy+Ip+j6mM3dESF7Lxl747VB//ACXWujmSEPnHBtIpDLUrvRlERGOm2K2xlsc2XJq+HnjDqKzn9Jyz4es+wTw9Qmd4/wD7WV/D1n2C9INYyR1PttkrmCPcrzPh/wAPeefpWVfD1v2C3/8A56z39Kyn4et+wXqIcIreOKrpB7leW8f+HvPf0vKfh637BIx/4fs8/Ssq+HrPsF6jLFMuEjpBeWvL34hc6/Ssr+HrPsE3+IrOf0rK/h6z7Bem8be0hDxISuHUn7cKcteW6/qlzNjXUUGPddqvlaQbHVpmDn5n6P8AW5UfZruXTB7b2lAUlQQ2yjtR0gvLXMvxW5j7/Q/C1H2aCLoBWYC4XFpfNCRl5x6RCIyt82u1C9ISGWnSq/mFrT8iG5p//wARJdImc12jcUg0vFMOOd5eVI9UotyrWbYNzIt3hU8VSKgc5ISKQ6opniN6MFJ3vCW5dL6On5xrukK5h0bcHjjHliuiZU4LZSL/AFK8Lo88dzS9BbHtJaro9IBESHtW27UoMzfc9qEi7o2rbuyx9HnknsSjqtHmUfmGcNNbhktU9BmVT7mQjq0iKcDoBVukJOXelaqcvNw3C6V7OM2CsYcaLTIS+/kVZp+jjBFLhiXduXccg6EtMD59tsiKMhIRJTzHR6jGMWGh7o6lUicfDgFL0TG4RYIpaZDcpij6CyGRU2rmFdwxyxsdIiPhTg0I9lPqrtXGm+r6Q+0N+kK0fV2W1hrwrs40w81q2VEN1xau8l0HuVxnDq5dj7QOq3ai2+gtS1dwtQxIZCuwk2SaxYIh1JdNq7uN1PR6spGjccaHhiJSLdFRXRkeOX/phkQlpJdxr6An2ibK4XBISlHSoLoz0OaoScJsRKVw3ff8qcx0zsVhqmqxKTlMMe6iDccjL2Np0+bGKv8A7GKXtfiFbwoeZtvuimcyyjnYVjelymIebzacwqsvL3CQ7iIdKv2NIPvTXpCkfg5jcw16IxS0v3soo7buVF7m2MrbkR7AyriCPm4ldq1K2OZNSOaqZrl0x/gqJ1i9H3W76ESGIja2RW/t/N+RKRXyc01/YjL3/OCIxLtbUhzq4o+Xu3fyVKyDpBmVIIi4LpNj2Zetj/cus9Gs6brGh2uCN0uZPqz7XJVvxcMD4ltzq6YiUR9b+S6BisRobcuzPq3EmpDcQ6ZFd/BVXor0JN06kS9wdjLsyJd6xwSGKdsSMRbbGUSKIxkXa8n50aHZyDNert+ROAI3bZaVFOdXlYO2XpLuxuDplpimnKxsdRfFTmI7uHD0IrOGQut3bbhuFRNT0NrmxIuFaXqiu+PZ1SbjElHVWfZfuIfSEUdFTl088VlK+3ESEhIdyw2HWh4kZStkuy5jV5NUyEuAJEMRjbcuZZ5kbrRlFwnGykTZSt7v7FzcuNxbcNnJVYedIo7Ue3VaSlpGJIqn6PvmIuRG7uqUaydpvWQiXKueV6OPDJFZcORXSKSbaovOCQjFXNjK2CutJSDdA0PufxVNV7UnlUSYdLSPpJNLk5uFdd3tKupA1yx7sUxIR1DbzaU8cd0r1insZZSC642IiRNXOFK4RVq6OZKBENTw4sDaI++9olXM0aEX+KwQ3e2DzfSpIM/cabFtsrR22xXo8eGnl+o5fOnTWGGxEbe6nhHcufUXTE4jxZDcO21WOk6UUzjdzgiXaW81HFasrJI9jHxKsUvSKm0zGSLPpLRtDJx0YkqJYpS3LXlUTlvSGhfIRbdC7tCpdp9gpecHxI2mtEdqYxJFatOn0U3i1cmAuGpC5qQiMh1CpR1rsqs9KHyAezdJGxpzjpTWlIi7SgmamQlK0tQklZ9VE46QjK0u8iXqZsm2xk22XDuIktnMdhgribK60SG1Qeb5sRC4I7m3xLwkjXKQGy86+LktI8TSSkKvK2m6ZwhbEi4D5f8A1FpQc4/Ln+OdOaRld3Vpx6pc0y2q5tdG2uYUe3lVM3aQivM6vV251hTVJc3iSTy58iiRejJdRCmph/pSipRK8Rt9HSjqO0UXIacmiHzZF2hirzR4EQ6CkjsuoDdcEW2xKRdldFyro7TUzfHqbYjIh1EPorTDj2j5HVWej/Q918hMhERjKUl0fK+i9NTCPEbGRbiGXy/kwVWzDpsLUmqQSiOkhiKrlX0kzCpIZPkPpRW84i5fXZWaxdfdqqamEtMdNooQuk1IO7dyrlYE+ftrpF6S05GI2lq5lfRzTO5ecnY6HN2Kn2u7mRr7wAMit5uyqX0AoSHiGWmIxuuSus7N+BTC0255wyjqjsR/9Jyukz/aaju856pJ1vpHSFv9VcRo25CJEWoRlcjm8BHcUu8Svoy7Oxt51TFuLwp38JsbSXIGH4yjKXNKSZOtduiXrI6n2dobzRgvdIlKNwp1uva2kJdqMblwTiujIhMhKUtSf/CdYAl50ox5t3ZR1Pu78ZtxEpS7qG4jY7lWervMXanL5PyJwRIZFq0rmnSLpTWMVhNtukQjq2jq0/lS0qV3XiDzLYOdr0l59/trmG0yLwj/AARTfTasEZERWlzCp0XZ3vEh7y1auJ03WXUjrl6MVNU/WS2UZC5Hup9T7R1HEU1UsC4MSES5lSqHrDpiIW4kMuZWiizhh/SQ3afvipsG23clpCtJofDqS6DKaamKTTYgRd5GgUvSGSXipqoTjgh6ivBu0tq3mDhNtOENpCMhXBs9z+rN10eKWoh9ZOQsq7JmPSNhoZSt7tyrGYdPWhHzbly5gIPuCJOO2lpFHUNM0VpCKvTPsmq/pu+egiHwqFe6QVLvurnrInCmb5Uv2IPKpXijvZlS57o4SYxF0vbDLxKdbY7q09SSEtPiuRsqrjjJFuKXaJTLFWVS17GIouiPm3Cu2aP8kxjaSacEpCY2kJakZyZTyOHO4UIVc+wPDKUhtjqiShqysfIpcTVarnXMtuttuiNxDEuYiFCt0LRe2CNq8/PxX0XDlOTGaVZiufHS4UUT+EKnURF4vWVxboqS0uH6qj87fo2mi4YtkUS0rOY20uX/ABiCZzYhK4iLsyTmY52RN2iQiQ8w2qtuV7ZRiMeaSSJGW0o8q6MePTgvL2E4PERXERJ5yot3IcadwtIopjKn3NI3d6K6OznuFyNcUo3ERJOLpcxKR/s/VltHtXDpSv7NVctviFL3B8XYWnzAhHV6ydxrActdEi9K1GNdFqkrpDy6rlhdFKkd4il7pfEplnM22Ci3IbbbkS10pfHS654pJhzow/uISjquWH0cJOcyb6WrVl3T50RES4heH5VOUPTkSG4SkqK1lkSHujpUs3SNjqa9Za+7Ki+mq+f2pbIZDLTaKrufZ1xRK0u9aop0IjIVFvuFG6VyV5FY+nU3Ps0cYdJ1sdUtSrONVXVzhRfJsR2yVjzFibpSEreZKpsW2JDbcovJpr8dE5Nl74FJ90nB2yLcp7Ns2k04IkUeE6OrsEKjK6vKVpKKfdIm3LtpfEU+8r41dWBhsdQ3ClOONiWkVGv5kPMKCdzdsZXdnSsdvVw9JN+U2dQ2W0RR+U0DtS4Lbeki7qB6LZY7XEJDa3K4o7VbavMmMracCkIeKNpORld2fKtccduX1V4+KJUHqbKmyExEn46RuGX61Ta7M6uscJwiIWiK0ZRG3s4f3IGoMny4jjhOEVxERJ1uoiMRtHSunDB4Vtyy2fabER7RLaDJwkrE1ouXQ0D5fWRVAHHcFsSuJQ+JR1Kz9XNALtXxBu4YlaVyKO7quTU4sMNjpi2MlxfrCzHj1hCNwtkS7Hn9WNNSOu6bbV5zxqHHXTcJy5wiItKWMRl5SVDiQoocRLVtQpYRiIlcm8XC5lonaRdqBEdyjzqZFuUPm2YEMRkXoqIwzMRIiIiVaDoOLdolptUFnFURONtiUYkMi5hUQ30idK2XZHuoTF0n3BkUtXZU3wcei+rLCVJISERIikMeUFxvpU8JVz//AMnqyXZ+r2nFvKhIRja6Uuaxea+ldc4Nc7Fz3Qo+NQ0WfLxEpaUZwuYRiqnkOZHEic5lM09eJXSkhFSRNN8ooqnomiErVGYVbfMjcvqWyG0pERJlKIqssAtNtvKo9yqfpLhMrS2kQkrSyI280bu6oTN6cSac7JKMoqV0bqx6UDVti0fEIxlqL7/3LoGBLzV0EzIqbMGIyuMmiHsxXo6mOTbZcyjJrjTtSMm3B5hJebOkuMayrHTF0l6XPCVq839YDPDzWsbH3ySnGnkVlrkmm+yj6fVJQmVFaIqTwciUluw1pKz8QrROR7I+qh23Lbk/ASEVjm0wrY4iW5IHAeb0f91sxEd3hTThjHclFhnSGUSQ9SMpD85Kfxut5u8hscZlK4Y22qrNoF5fVCJCJSjpGOmSZr5CQkNpCSG8pNlaUiEtRKQdeaEWnSlArTj76uXPDVel6Pm/0jHa90dJXF2lDY5TUuyLiDdtkpQ81Ybd2xt2qwt5uwLYxEfDqU4ZYx183bJTsp6MOuOjKMZXXSV3YySmbGMRlttWqbMBKKJ9liP9Sq5M8OMptlseXsxEVvyDt1eFJwqx3aU42+BEs7k1x4zJCZdr0oph8XNpR9KSksHG+X1UjEWylEUSruKMYacEpESfwFzmtRgCMZcqT5BK5CdaDYtDzFLvLXAiiPIKZddEf6kSA2bS06+KEqa0R5lD1WZjzI3oSS+EhmFeIqCzHMhEVH1+ZiNxSVYzHNgKUdSO9K8UG5nmpEo114SuUa8UubxIjL2ZKbV446EBItKx6kciXLwy+KrVkeVDG7VHl+MpTPqJsKZwmyEiFp3TErYKYMp4QQtDKJSl2RkrB0T6NFXODARIZRKS3n9ONG7wIiRQkXMPd/xVQpOmr9HViDUhEiuKUdX/AEw/wW+HHuufm9brxHY+mPSBjIaH2IwI8UhuIY6iP/Jc/oa8nGhfMpEdxXfSqV0+z4q2pGJEUoyIpf7KfyV7zDTcrR1Lsxx08jktzu7VhbzMtItyHtInjEVxDH0lVqyp4UiGRCox7pE4JW6dwyWkumVjofHGN1skBVZm2FokPrfIq3UZyJNDd94qAfd4soulIrYpXJUxXvHOB3aV2XqybEaQX9PEG2I9teb6Nq5seIRDIRjFerOjNHwMvo2xtGIlGKntsdVZ66MzcYy8RkQk4QiMS+MuLZHWykJRkIj2lef+ICvFxximlpIi5tID+xc0oMBaHVKUfRTnhF8Lh7NHcWlNP1zcZCMlUnXS5reW1Cv17g2q+yP9pyqwFwpSj3lF1QCW4VHY1zhavjJJPS3J7VpJYRjaSlMuMR8JXKrtOx3bpI1qqtLw9pGRyPUfQkv/AGYIuS82UuzavNPSF3/1jso+2FqLVcvR/QQY5MIjbJsi0x2Ly90ydjWO8pEQ2lK5ZtKl8rfFthwoiMiUc5mhC4UZfN/cm6E/MEO7aop0ymRaRJOJsS9Pn0d1yl8sz6MpW+lcqc2A6tyJAoy9VFRp03Lek4EO6Q9oU+7njZCY8wl2tq5SdQ4MYo2lekMiIpJHJVhpa2NS0Q2xdEpal6m6OvToaRzmbG70V40dfi5q3CS9bdAKjiZRQlzNj9/yYLPNrItmGMV5/wCtxggzVwo2uiPxF6AwwXGeu1rz7B2xIY9rQSjG+VZfTntC7t5VJC8oOidKRD9VSOBRW+2Qwai7ctvVpCO7UhDOJRu08qFqHB0iWm5O+RodjmpDqTIZkV3a7Sjni5i+KmzK21TYaRdrSWqStKUYqJdfiMeygG8yujIhId21TarrauruMhIiGJDp7S3ktOVW0/SFKVr49n74qEZfM9va1SV/6vsodcdcfIYiTESHtLPP6acHjJQKrKCGRDIuHIh7ya9jPi3IdXLyq2VZxMrS9sIdNv8AJO09Q1qiMfjLgvh7+N3FaCrf4cokUdVyXTZq4UZNkPpaf2qaqXGhIdIjuGJKuVTxce3QXL/PBK5VUm0iOZDKJF6SJZzKO6UdXKhaKsbcEW+GJOSiUhUq3QsFK0RjEZR1ftVynozjnQ23D2Yl/P8AMiKbMpCXaLm2pNTlLUm9IxLaKQWXSIilp0qonQ1vMRESlJNY5uIjbIlGVVMUSGVyinKV8dOrTajZdVjczgY3WqKrM5Hn1bbVBVhOiJfFUJUVBJbHVYK/ORjq1KDqMxG67dzKJfkQ6k086ItxjIuZRs9QxmFeRWiW5RmGKkG8vIhIiKJaoodmnkS0jPOW1IMUzjoiIjLdy2q25DlvC9sEU30WoBiUiGQjzRVpYysjuErRuUWKgfMWCIRFqN2rmH0lC1uVVbTTkXLSafkMtQxVvZouHdqlH0UVmbQ+xn9pcCp0xL3Ik5PKcvpVendcQ5o+4MvONx07YrmWYDJ0juH0Y/xXtXM+h2X1JCRNNkWn2sdKqXSDq1y0pRERIuwMl34yx87nve3kk8CJwdRK15a7FoZAXoyL+C7rlXVbQ2iQjq08MVbKTqyy0RG3/wCtv6FV3FSvMVe+RDERcG3TElWHm3ZkUS9IV7MLq5y3lHwCsLq3y3lHwN/Is+2R6eQKdg+GNpCpKjaiMibGUdW5eqvxcUPKMe4KUHVvl9pRH4Nv6ES2nNPOfQpsn8ypgcAuHKREQyXqanqWhERI5CLYjqEdv7lGU/QCkacEgiJD/wAtv6EYHReJe2+qP0KpsV5t62K4n81ciJRFwhGN3Kqs+RiIxEl6lrOrqkdc4jhDL/4xL+KaPqzoyKVtv/LGKradR5cp6lwdpIkBccjbHwr00fVfRx0jd/y0n8WFJtiP/bS7F0eZH6HmGX3/AMkjCnIbYlHbqXp/Dqzpox82XeFaLq1pOUSLu2p9j6vKzrRCSQGDshjLWO0l6ld6q6YuUf8AtpGPVYwJam/RHcjsXXSV6POCOVNiW1gpF/2v8l5W6RNE5UukI+6F9/yr19T5GTTBNSERII2y0qoO9VjBFLiBdIitL5MUbDzfl4kIxJKeppaR09leig6q2B0uNfBl9KQfVS2Ol1q7/lufSjsrq8540kR0xWNtl4l6Dc6rJRi616Tbn0oeo6o3C0usesP8cUtjTz84BSinCaJduLqgfuKTFva/mmR6on4yk1LvCjY04dUtFLUUor1f1SvS6P0Oq1vd3lzU+qOpIhHiNDdyy/gur9E8lcyrKxpnSFwmi26btKWWRyLm3pHurl/Xgx/6Zp3cLojL0V0LL69twREbSEfWVY6ycvKry98REScEZDLmUYHZZ9vP7DjYkRbkS7USHtKuVjNSw4TZNlId0pfwW6cnzIYtldLcruRdRNdUvyttLmGRer/goeqxq2/OEZFIuUVZss6PVlScbhHVKX0q0UvV0+VxOCQ7Zfy/Mp76VMXJK6qqYjcQ8pR/kjcpqak4jdp5iXaMv6tmruKQ+kMlM/2dy+jacPhBIB5Yy/ejvs+vlyJzKXSbFxy3m7SFp6QZjIRiMZaVY+kmfi5EeHwx2iJKpY13eiS4+Tnetw+j7Y7d06I0OXuMCQtNyGMvNlLSrTgTTAWjEd0V57yPPzpiGJFGQ2yVrznpaT7TfDIhL3S4R/iieon1UZeisy8Bekplx33KQrRISjzKkFmFY2+W5twhiMitVpoKsZOcS6QjFSQsNORskUbSIRtUW7deONk1ULSU3HGLlpFbHV6yar8icuHi2jtkMVYyZb06e6KWzSjzSjtLclY0cydyx+mIiEjltIZF/FO5TnNXSOALnFcEnLpStXSHmRiNo+FI4FMI3A2Uk4SHzbpS2LZcMXJbbUJl3SXiMSISIiLUVqlMaGmcMSFvduG3xYKZHLmhEYiPdEVYVluvcddEokIx+/51Y6JwYRKMiLdFPtUYDLT3ULUUQk5LagguYutXDwxl2RJRmORMHdECldu+TH8il36AiIolGXN/JLp6Em9USt7SRbVap6JtxI5BEfS+VRzXRsuz6Su9dSk4MZRuHVJOYZcQjEnB0jtRINqRU5DbzFpjEiQTGTGLnN6OlXnChJspcTwyUa9mDYO6SiVsrudMtq+FDUgRCOki0kp/KHHWyKdwkI7rVMng3baifNRjEZDuQQTCsKVwpOciRUz5CUfMP+HhEl5iyJCJDESuj9/2KvZm+XAdkWlpwbbtpIn2V+npEa0kBXVRSLuoMHJSitPaV6ssv0+bv2coK+JCPaU6OZ+FVlvCPKigKSLFLGNeKcxrhVfDElv2TGSXUrVhwrh3aU5jWN6RlHmkuE9dPSl+jJhpgoi4LZFcXa+qqj0N6b1blTTN8UiInbhIiu/kl1Ldensa0ZRTwPy3KtE+RFpuiJF4UXg6jqvdTxHpWYOd1QHELmSuKi4msWBrMSVd4y3i6SXQLE2SWJEqd0hzj2HTcXlluXP6Hrek5GI28xJXAbdwwJYOKovQ/pe3mYmTYiJBaV0lYcKhLonaXIh8Sy2W5Q2D2lOey1UkG0ot44qL9kJbdSKejuQ7DHspQoLGsiswrR3fFS6l2HeTurPKgfZA9rwpfsoe0l1OZixiPKq50zqeGxbbIh3KabqRIh7pfFXMesvO+FwG+YrvWU5ySO70HHOTlkpvHOCaIYlH0leKKtF+mEyISEhjHUvPeadJxkQirl1UdIOKRNEUoiRXERLDgzm9Pf8A5X+NmHH2wi/VPRqjdLiE23q5VjXRykArWmx9G1adzJuRDdL4qewftG4pRXT4fJ408FMw3cINjHltWpDK1Ck5atAJFcpuO1yluuiREMpRu2qjdYtdEWm2+bm+hXOpm2JFEbdS5T0hk/U8UnBERIojJYZ3q6PT4XLJUuluEX24jaTUvSko5nLXHSkIlG0RttVjqswpmhj7YRS1RKKiq3pNtaERXn5+a+l4/E0YzDLiYERu9ZBcaJcvKmajMjdKRFL0kPi5Jc9l2rKJzL6wRckUhGQxjH5f7leaIxK4StIeyub0DZO2jar3kIONDFweJLSS6sJ4c+eOhLzJXXEMdUkhtwh0lLsin6gSckPMo3DKHBcEhct3LWshpm47pKP3/ilDRlbKRRuSqZgh7XoooakQlIh26lMAV6RRiJeqjKOW62Mi9FYdQMrvmikE1cJWlO7VK1VE6opvCVw7k5gPpIUHIjEdqcaqZWo2R/DBKwwGOqW0k0eMij/SmTetjEbSjaSAU7iJboxJIcOW5aB4o3CmX3C4bhCIyigH3qaW7V4lH1OUtEVzY29lCUNS6RcMhjHdKSj6t2padJxwxiRW7hQEq7TkJS2yGIxTpMkRWoeizAiHzhCQiMpDFF0tc077WRW6tMvD/gjYJffMRFuIqr58wQiXLw3SLwkrdg+Mv6Uzm7LZUzpS0sP/APiJCcvpfsszFt8RJuN3dRnH1CuG9AelpUzgtOXCRFcUiXYsvqm32+K2QkJbR1Lr4ObfivK9T6frdxIYYpQkQoZt6UtvZW+LIYj4l2WuIew72tSdxcUa26ngcSORHdLejLGaCJH7YIiI2iURVcybq5apqlupIrmylHhCPyK9C5/Stk6g9HweEikPKnsCUU3UiRF2U/g8g9D8HCWcVBC6lcVOAZg6KVgaB4i2LiCqA6z6Z1+hFtoS1F7XLk/yXnh7Iq4C9rcGRcpL1XiUhiUYoWopGCEZD8VTU7U3qPyl2mYqeIJCRxISiQyt/wA10V2pEUK2QtjEbRWHiKUTodTuyGScw7yAYNOiSoaGYktYEhscO0lSip2oRisxxQ+Dq3NMdT8lvBxDYmswxTh9RQPXeiQrkHWldWNN7uGJdnSS6tJcZ60Xv/d2oy9qEVz+qusNvU/ibJzeXN64C4jmrUStfVhWcCubEoxckJS1XKOeYkRFcN2ntIqiYJp1p8drglyyXgYeq1k+29ZnhycNn/w7RUMDxLS1EixtTVEfEabc3EIl6qzE/EvoOK7kr835JrI/JF0dW0wRE4Q2jpJRwlcKo/WpVuNNNk3bdd2lpldQTaV6edK2xpHOFqIoiIkuF5lnLhf1Kbqq/iN3FstHtKpFhLVzLzubPb2/Q8fjZs3CLUSyKeFkSJFDSjGV3oxIVyXJ6kBtEnAxWnMbiit4KaVqTyhzhPtFtIhXSscJjaURlLUqBRsDw2yL+n6FYsvqRiIi7IuWa048mWd2nniut0pJv7Vv2REJWyiWokMdaItE6Q3N7eZa2smPuEIuRlKO0ex/kqDmbxMOOuGbpSGI7f4q6nmYuDIW4l9/8FCZiw2/GRd4bfW8qUNT6npI+QiMtO6RSJTnRzpW5Y24IlGUSlH+Kj8yyljSxIiEroxIf3KIbpDbdEdMS1FIZKtzXgOmYZ4REJcMRl2lJNVDbo6olqKKqNCDvDEuGRCPugiRD3Zf4oD8IVLT5NN2iQy03D/JTBpb8vrHCqya026pakc1UkL4tEIxH3SRSL0f8VVxzuLZFESdFu6NooSnze2ZODKNwyEv4flVp6uhC4V0lHVlfwuW3m3KuUle4+EhcEY9pB1jT7giRHK+67b+vFGxpZMM0Fy0RHxJtwBNuLkbuW0rUHSYNiJCIjIu7JOutFESHaNw6VPYjuVsCx5pspC7zbVp3KhFwiEiEpEMRHV9CCxbccGQui0UhISPb2E9l9bIyk2RW3FzJ7hHaNh1gi1FLmRuccQaR0huk05LwEgM2ry0sDqIhItRD9CLcqn/AGM62cY+xn5EUtXCL96cLL6QlR0fFohISlzafCp/o/nTtGPC4XFEpFKRIKor24xbKTkhkMSSmHx2+kOpZ4240cmHZZcelzcouNC12rSR9H0oYL3Rv4o/yVMrm+KP+kZKGqOj4kJRJxdPyLPDiy9Fv6djo82YcH2wPRIS+VGjWtc4rgo0D7UYulbbFbCuzACiLhFy2yT+UznoMtu+YVbfvg+qsdrmBH29sfSXAqnOswEouFItOmKZwzl/SUpbhU31TSehruWX5q2++62JCXD0x1EpocLdQrz3Q9JnKZ0SGIkWqSnS6bGWkhl3oq8fUJvorHaW8C7ywxJcdb6XVI6SEfS+hOtdLqwiGUSHvfzVT1KZ6OuvYAXiSykPpLmdT06fARtGXeUM/wBYdTKPD095VPUxPwsnZcDScDu1LmOX9Y7gj52ml6UVPZN01pnyITi1ujxEe/Kzy9JlFvcdjq3JxkSc0igMuJp/S4JDKWpTx5rQ0Y3ONyjzirmbPLi6n6fLjIZJ7DLTVad60aFu2Yd2Yp3Lus7L3yiJN/CCtNs+qafpiFMt4ju+cpmir6asEYkN20SEv4JrM8uERkMkbGtIcSuWYkmccJFIpS8K15eXSKdoFSSSIhTGJrDJEpCJWri3WRIs5CNxCLRRXZMCtLukuJ9YFfws8aIokAtjLbt3eVZ89lx8u3+O/wDUOgAkNw6fv/gkG2PDK3SMu1qWN5/QiREJDLvDqUXUdKWouCLYyiUbhXi4+mx7dntcnJzbsk/xdeyB+VJTF2BFHeQlDdDnhcyqkc3EMo6lKv1ANCRObV7XHZMXg8mO8tFt4xJU7rPY49NGUY3CUooaj6aC++TUW4yKJCQ3Ivpy+0NIJEWqMbtSnPkli8eCyxx+sGIt26hHchHG+USUjUOtCVw+sgirxHSI+kvLyvl7vDhrEtmmkMit9FLccbAYiUvSQFVXG52UOLktXxlHV0SHSO6SksvouLHagqZqV21XDL2BERIR5VOVZ5UJS5aZWiVrZaokj6SjJp/bG3al45pw4iIyIiKUR0qXpXBIZFGRR2j8i0wY5GsG3HTISiIxtQlUyQge4uGUYqXwwG0o7u0hMweFsRiNw7pfNWiUJltU+VIJE1EpEMeyiWcvEhLVcSNp2CdGUu8IiOpFg3G1MkSzkwyIikI8wxUdneRiVwkW31Vb6ohARiWpDkAOjFGw5zUZu6w37GiUdxSK5Q3slxx23UQx9FdNrsoaIh/0qDq+jbcvNykJff8AMn2Ur1biLQA2PthDI4kolhgiOI6i+KrZjkvnZFIo7SkphvL2xEXBGJRG7UjsPCuZY1whi7IR2y3LdbmYCURbKO0uIrg7lnEGRS8Kh3cukXDIbR0lH7+VTcjiEoa+RhEt3MrNWZq22MiLUMSjqH9qCLKhbEibauG4bvm4ISoyBwwk6JARRLcX+ycm0X7PvONk2RjEh1RjqQFFmJNOSlaW0kuloxYGBOSCWmMbe9it5q9TcUeHp7275E7Afp6oRd4kic4pSK2MVN41QuMVJEPuTke9ElX8kNt8iGMhEh0qyP5IAtPmMg8w+XNK0i/UjXlOU8Csae0f9KeBjw7dKJxG3TpEk00OlSo2NOIjFL9jp7ELuyW5MP1AyIZWosPZAtN825O8ANKzARjKSJxjbDlS0NgKqnbjdzII6Nort3dTWfVZCQhEoy5Uiobg0JSlP1Uah7oSqyZt30iKJf7KHrcgIfaylFSFTmosNiA3Od5A5fWC4JCZRK64i1dn86WlVEO0To8xeL5E04JjGJF3fOKzMN223DK4pJb7BORHhkQp7G1VKqdtkXxkS3nDgcpeJTVdSN2jwy03IfDJWyCQyG7SjZ7AlnLkpRH79lMjVkZkWktQx5U5mOVOCVokSBZZMCuEk5YepftcQ6VFRMREpEW2RcqomcdI6l8rjIfSJazXEpXKEJdfF5jyvW+KdOrMtRF4iW2qtwSkJF4iQ2IpYCS224JNuhdEusmpoSb3AJSK76cV6M6EdZVJmQNjIRIreHIbf24LxkIKbyDNnaN0HG9pDzD8qXZftWva+YUhOecbkQ+j8ihniiJD2lUegHWKL7QtuxEiIhuIVb3KhpyRNkrxzmmfSsZKQ9pLHFDUeHEcIRIZIzClJXE5QmXxSJcC62zEs1K6RcMdMl3upx4AuOOaRErl566ftE/mTjjYyEhEhIdK4/VZ68PT/icJ7naqniSRL5yMxy5+MuGS2GTv8pCS4ZY+m5rJPDs/Vk/xMqYj7nIS7siUl00f4VC4fZFVvqmx4TZUzm4SKRc0kL1t53ECpBIboy1SFdk5P8HzHJxW8rmlLW8N0TG64iJTnTLpD7Lbaa22kUblVAKK04creVc/a162HDLGzdkkyWhCXN6yKYoJc3hU3TbroI4JJ6lpDK6Nqm6TKU5VNixaN3NFTc/BEs0LgjtEdUiKKlaeuaKLTZSK37/kUdgBVIjdER1CpbK8rGJ2qZNsskrR0wDdzff+9FtYDKRFHsxQwt2jtEUcUYCS1kZ5FiV0UC4152MpXDt+/wCVFBUNkMhtEUAyc3LdJbiTSmWWBEbbhTuLIiMkw3xI23RTrco3apJwEONSG5RTtL51uMYx3EposYj3kCdPIpXW6UUjjNLEpFclOxErt1qzFyJRJMvP6hjaXiSBx1touWUU4NKIiMfmptpndzDcnal63s6RQCsQjqj9/wDoku0wltUc3XXELsdJENyVR5wJcQXIj3UrBsSNMIlIRiRbhWVtLxBjqTlNVi5cO21OcQi/qVyaCr1WQkW3tRUDV9HLdRXFy6V1AStQp04uF62lEoVDonloscQZSIiG6KttZEaZ8Ze4VMfgiS6bLmxKUtXhWq8R9jPiWoWKmPwRKp9oyqPN8hkIyL1vkUfV1b4jIWiLtRUq/gIylG7xJOBDEo6VmtC0eYVbjgt8Moj7oSlSoplIrdJJnFyJxERGO7s/r/Mi8XyuIdw6Y/OTBY4xbISt26dS1SUxSG60pJMy3CN13L/BENVrbYy0x2oAZzL/ADgkV0S1LVVlkhkJDGUiFZhnDZHGNxEUZEpBx846hKOq2SQVWlyBsqkXCEt1pIXPMmblIW93aVuEyIZW+iMVtxoSGRD6KZ7QGUUXCbERGVsu6iazE4WjEhK6Kk8GBaK0oyEu1FLpKASuIpeHUp0cqm1LL7oyESlHSVqcy99rhxIxlyy3KwV7TYlEdRCRcyhmKMboSEiu+/8AgjSj7ES5f9Khs5JuUeFdpkKJzOldLhxcIY6o9lSNFQNk2JOFIxG63cg3OK/Ln3HCi2RIBzJ3RK4C5l006MeJEZF2ko8sbcltKMVc5LI5uTimd8uYUmVOObYipZrJWxuJwe7uJW/HK+EOmPa1KKraMi0kWnlReW1pjwYxGfglsxKI+ElF1OX8MiEpdlWvK6Fxsd121bzCi4hFbcs/dsbezgquX1LjDgk2RapErvSdLnyaEeJGPLH4qqtRQE3K3amKY42q7ybhX02Dr3QHppwyi64O7UK6Xl+etVMibcGXKK8wMPk3cJFJWDJuk9SxpMpbrlXH6i4sOf0eNnh3bpZmIsUL5lcXDivPz2Zm4RFKNxW8qmc86ZO1bBUxiN2opFpVRko5+T3Ltp6P03tpxqpcIdReqnsSIi1bVHU7zcRu+Ki2XBItS5fLuyTnRtwmnWykQiJXKtdMKg6uudIRkMu7apZutECHTErfv/im3aikEiKWqVsVvjldOXPi3ltXafLzLaSMp8vjqtT9ZnTQ2tCJd61QlRm7hW7e8iY2uiXUTWLrDQ7Skgn8z5Ij6KhReIt2lONYjuT66ZZZeUm3VOlpIi7oil4AZFKJImjJhsbR1bpSL9yNpnGy0lIuaJLKwbuiMoZInNJRV1YZFsRESlzKAp24xtJTAvDaJfGRImiaelHVLtdlNvYx1bluqrybERFsYkofMakiG0roraVGkm9FsCiQxK2Uhj+5Ch7XIZR7KgsmrCbB1pwiKeknLol8iPy9w9JERCOnaPyo2LE7ltYPDjKPZS6uubGN+4bSUW4x5uQlGNygKmTjoi4ZCMtQiJfKjadLm7mbEhGQyt1JwHxLTp8KrLGXNOlHilIC1RHT9KmgoyDzYlbHVuRshpgUpFyyTA4iVzhJbuLhCIyiI26bopbdEQiMilFMiMHiESjcO1IqKiIAUZFKUY7f2rKkhbbIi0ih3czpCaEiiNttyAZJoi4hOadsbvv+RBUeVcSLktRdkk5VVzbhSbck2O0tqLoKxsRLzgtg3pLaRfrQBtM0TVsSEdyjqusdlYMtVtycqqsiblxBLtIamdEW7olzST2SWoa1yF2rcNqNbeIhly8oqHpAldLVtRwVUbY3drlRKWhLThOWl6v8kLmBeaqbS9of5rfNEn28I3SESLchswOLD8i1MVMbdRcIleKaQyQvi5LbpT1C0I23EMlVejebk4QtuWuXauVWRtwmy2xLsioaFVTbZEUZD85QbuZEJE0IyJvcn85rXytbGUo7dqbZYbFvjmRC6RRjtSESbXEcj3dKIdoJCKboCEbS0kJRIVuqzAWmyIiLuxQYXDKxbdGQyjdIf5ImtfIRi2Po7lHtdIuKNsuyMSFIqswfBxgyEpOFErZWxQcTNNVkLYytLlTlTXW2iN3a0qDqXyYudkREUhttH9ih6ytGRE4RREuVAWIcwbEpO3DdHvfqS6PN24uXDaQxVTrs0aICiQjGMtSrrFY6RlBwhl2dSZ9XTnTbdkRbbRt1JFKTDhWxkIxLvKAyqoMmhbcHzg+siadlxsi2kkek04BEPZElomSFu26Q3Dyp2kwH3QroiSewuIuHdt9FI0d7XIiuIhHxJpqUhctiRDuRNczLUWlC07VyD0kqjAXO6g3aRuQ2lEd1yk2BEoiWlPVQiXDiWkruyp0AR0DdtqAzPLiGJDqU1iQjK5Bk6MpakukKWq3U5YRELfMoiryQhMhEdIq011SPE0xEU3CQkXZU6sbS7UyooibihDEhVtp6LiOzjaNqFzfKbiiPatSl/Wn2rRuCmcTRLtGYlpJDHTkP5YktZpUum+NpTrdaQoMsEjyYqusTchTtaRbkyTxEmcQJKAC5VXWRl3rZEkDgRfmRlLTkX1lLt5YPDHUJSuSuUiLki6WgcLbbFSlPQCIXRIi2yUi1TtxEW5SlEiiSc/BlwkIyisrlsvvyDayyV0hGIyG5E1FaIR4QjIdURtJaqmCtjJMOgTESuGW3akpZcvdJ9sTIYly6YlFakJEUiiqw/nzgNcNsoiWrT8ZaoK4yG64h091Gv9hYOkFaY0wxKRXd5RNI8QjJ0olylqRrkiaKXKWlDZfk5u8N0pRHtDqVbKeRlGLcboylbItqlKUdoiPeFV+kyt1x1wjIhFsSIdqYyTOHW3RA7pSG0bkbKrRVYah9ZMY0AvyG63d2lJs4CVzlpEOkrfvikey225DqISRPJGclyXgEThF3fvgptoJCUh5o7kBS5kRFpt3I5uolIto26U2eRQsiPZFJceEtUYjqKSiMyzQolcq9V1rpMFErSLl+cns4suY1VMTRSIe7JVLpAyxwh4erVaUiQOOHFGROjKPd9BRr9WQlE7o6bk9noTlLxNleW3cn6qviLgyIpEJaf5oNs+IIlq7KewoScERG0t0hikNJbJcwYMYuEVukS3KdpWmykXufLpiqp+CDbEXCLT6qOo85IfNFKPNEkFpZMKkW9sRHTuSal4SIHJR7SreZ5p5uLZbu0iqDM2ybESdGURGJIiU7WvuWi0UnPVQ+fPuFTSGROcMhIdojC5IZcuGJCXdJF1pRYflq4T8vgiVS+U5TwpuSYFxRcckV0RIoq2i7EolGMfv+ZQlLTk2IiQ3SlqGPiTtbWCLgjuIe8Mkl6TwmRRiLcea5AZhiUYkIkIlaQ/OQFFnQxcEi06bU0OciREMhOPZSVjE9R1QiMCLUO7b3UDmlZLzRRIdMi1Ku5nmIyEhLlTr1S46I8NsiLcQjJBp1tjhA25EYxS6Or4ok6RSFoto2/fyKOp3HSai7btiVsf2JLNI6LXDaIhkUtQxL9yAedqPZlSQXcIR9b6FuvytuMbiKUu8icuYFpseJ7YRXFzfsRRYMSIiduiURkmNKfW0IkJRG6Q2ltRn9nuBEuIMijEZKZby0ikd0S5k3mFE6UYlIe19/zIMVlFCWpzV6pCmKusEXf9PKiPZBMNjLVHvF2lCvVok5ErZbo3KdmnKdwrT08RStOQtkIuFEdUiiKgG6xtgREyERHmGSYzDMyftaIYj6KexpLV9YNxiXmxIvv/0WNvNi1xBIpFzaVGtZZUuNw4YxK6RI0sucFq4hiPzUpQkqWr4mqPo/zTkxESu1F2VXG33RHTqLu2ogHXBGToxlpuH5EHEw9qjK2N0v5Jkm24iVpXW7h/eozNqkhalxBldKMVEZTnTZEIuORESLs3bf+qD0slXTiRDb6QxFMYEIjwyUhR1LTlolIuaKHqLZDEdXetUU41SttxIRlLVKSZJqJERXWx7KUzU8O4rRT2Fa2RSGN33/AL0j2CxZbcES4Qy9VR9VlwlbD4qmHajdaI7UI2ThERXRLchW1aqsnIfc0ikyWe0hU44bsinaO0huReW4iZcMSGUdxCKfaixXs3yttrhiI6h1XEUk5S0LYtW3FpKSsvBbIuG4QEQ+kn8KJspcO4tyNosQlPlLbbYx4npRinSpBuHl7qmMWI2x3eqmK0iEZC3K6MiGKRSI6iauKN0vVU3S0lpDEfnJrLKUikQjp1KWYARtLcjRq7VNC24IkOrmQFdhISEtMojzKd6QMSJsh5SVRfdiUZbpR1IBmppGiEW4lIR2l8ZSGU5eMbdpbtRLeL7WkoiW77/3ogSERk3pIrSFVN6LQp0W4EUiAREuUdqGyXMHDaIRK0bpFqWjAnWybl6w6fpSX8rJprzGkhuKSegPrcyEaYikPEISGUSt7f8A1VeyaicI23xuiUZFpTr1U0VNwDEeLcIlH537VIZDiXCIYxEYyK75Uh4JriqePcZCI7RK3wpprM2AlxSiRFuReaYjwC85p5ZESpwUbrpebkV1pEnilcDzsBGLcXHS0iMRu/b+VTOUG+bY8UeHK6PKqAxQVLbg2jzDu/2Uq10iqaRwROMBlpul/NUmxPZ5TjdErtIj/qVbq8H22+HGQ6bbh8SVmWfFWebaERlG4h3J/K6N9wYmXduEbv8AZGtFA+X5RxBJ0ijEpRttaQONCL7pF7nK2PKpZ2gdbJ0SfkJCVolt/UougfhaUR4dumX+6FJLLsqEbSEo9pSY0bQyK4e7FA02ZWj/AKlL4vNOCJFbJFAb2I262WrtbfvgmxyYSGN2nVbJSjLIDKJXF9/70S1aN1opJqtfgoRIhISIdsub9S1V9FxiJDxJEUiES/kpb2U2LpERDHtCWpPeyHRcGQlw7YxlFNKPoKTgRE5R2kWpEuGT7D5XWtO6d3miRWZME6LRN3COqKSJcNh0SHUw/wD+Ik59lYqblSRXaoy3INnMylFyNsoyHSmKcyES4mr76vIhXmiItO1JpThPNkRERFLUMbRkmHXSIpDbbqEtST7DdIpQKJbo2p+lpo3FIo8sVXiFDLTD7twjKQq6ZGy3RsScIhKOooxl9CHpItNiWmQ3SEdK3T48eTZFxJaR1R7SVp62Qw+4+ZkRCQSKMVM0hCLYluHVyqNpmSaGJDERkUo3JujqyJxwSIiH1UjkHVDhOy22yTdKIyk5dLSKDfrBbIhErVEs5k5Io7e8imuTdYRFHS2OkURgIkQyLSJWqtUtW+7G2OnbFS+Z1bbBNyltu5vpQWialyRXXCO61AVuDfEAgFTD8SESbGQldK2P38qbBxgYi6TYlulwxU6FqCzCjddiLctUtyHy6lJoiFwtJaZKyuvtDcJj2SEhTDQtx4hE2RFzRIlPbQmSWxfkLYiW0e8kug6VgkUdw8qDbqWpBGNsRR2MhkV0i09pPatkZiLbcRK2Q2qE6Svxai2V0Y2o7OKR0okfEIRHmjEUANIJiUiJwtsiknKast1zul2RD2kTl+XiLrRGNpFzfNUy7lgkQk5dH7/rUoVDTEOoScEZDdGKLTNZxWexICIytlaO1RtLm5PuSICEdxS5kjOHx4dxScG24tqbo8Ci3aJERCRR5UgnXsfNWxKRDpTlO2QiIiMkyLxFpEbdsdyLoyIhu1EXdH7+RIEVFCToiMo90opzCkcFrUQ6RG7SXOiqh4WhGZalvGuYiIyu5UaMDhQlEhcKRS28q3S5aI3CPnOzqTj9aOkSu9H5EUxUCXDISHtdpLR7R+GTuNuDUkTmrcledEijaPo3Cic2rhbbtKMS5lCDnhcIo6iK2WkhkmSxsecbIiEhc3IauCQCIykJXRWma9txiRODIhu4dtyj2M1AiNsSIiG3vFzI0SRy6qJoXBtlEfuSOpSIRIyjIlCtlIhkWgZF2u8ttZkT5EI2tj6IoA6tqmyKJKn1VLF3iah02qdpmBd4o8TSURuuinmsubjEpFdqFBItikacKZRuiI6Y/wC6KfohbER3bRSwbEXCGO4U5WARDIRuHvSQas0lfw3HRcLu8yFr87dK1sojpikV7oG4ZFqjG4R2/wDRC0FCTjgW2EVxbVpComheEzEam5sSHbH+Cu9WQ8PzQkIi33RUR+Dx2i3Ee7L0lKlVk2NvLpiJS9FTslRrakhO6W0StlFNZpVkBALUhiPLEkXmFTJxsibkQl5wY/G/zUfUNlxSIh4crhlbEU4BFLmbpFFwe0JR0pWFKLguuOERFK3spll0ZCNokVspKRKicGUSttlyoJGNALDm7m0pfsx8ok2RXFpuRAUgTk7d4pR+hJiIuFEIgJDGSBYk8kYdN8RfIol6yMeycW3TccHzcpWju+lRo58LZkI6uaKJezg3TBspXXR26f4ppSJNtxGz1Uh1sdqJaIStIhjGSAdcIrtIpA0wbolIikIop3Oxlw4pWA+bkcdvLpUYTLbhETbYyEtv80ELffFwZRHmu0qKLpE6NhQIZd6P708YOiNol4UHjTXDxWxER5tSBpbuj1cLjFwiMttwqv5vnBtE6FsSF0bhkVwxtQGGbAwcm7pahiSjc3ruOUxttKXeVyeSyngbUiMiEiGPLpTbjwjHhjpHlld+tacyqpdInNu4pI2j6POEJSG4e0ouooHS1BuEIjKXoxRboRKO4tojpVgy/KxabtERIhuKMkzV5ddJZZZfhopyojESGQiO7chmq9tohIZCUtvKpl7IXHRGJW26lpzotLsx3CW5Vjl+mBrc44jbgyK4bSUc1IbZEJRlJSRdGXJEIlHvFK1Cnk5tFcREWmIqxtG1b0S1ERfOQDTpCVtslMfgV9yRCNvoyTTGSOmURG4S3Wq5qQkzklURsFIrhiKdzMW3Sk7cA6RkswyOrFqIiI80UNS5W75xtwpesShawfhBomBFoYxERHtKGcFw5C433S7KAr8rdYGQuadoptjOyG0ht9FPSanvwX5qUuXd9P5kqmYciIj3dQxkm2+kIttebbKRDEpDtULU5m+6Xm5S7JRU6Gkk2RNuuW8SMpdkh+RNN9IHDKMiHvKKD2SV0ikWq7UpChyghIXC7Pa+X8iPB6STvSEjbjqIpCW1ZlNWIiQEJS23D9CRR4OzIeAEdumQ/fBLN8WiGTd13KlDGk6LjkYkMtUopuooicck04I2xKXMgsKwnXB23R5VYDEGxIglON2q5ByKrU5Y6LjZOE2Qlbb9/wAqk2KSIyFwYiMbfmoCqrHRK4S1bij/ABTZVzpNEI2yJBjuI63cJW7vvinyzArZW6dKhKN5+JadW5TAgRMXxEpcu79qAdzHz7YuCZSEruZR+YMO8VuNokI3cS2UfpWMtEJCIlISu5bljucuCUYy4VulAZFxspEVw6rvi/kS6bM+AUnCIhO0QHUJcxf5KWy95ipYkRXaZFao5/Lm5CRQIhKQjqkkZVFgVSRFEhAi0lcj66gbIeE22PEG1AC+bbgiLcRK61E45gYucUhLdyp7Ko38B1heaEiEm9QytTdA6LBEJiMhKMlJDn7huRaiLjno2/InMsoi47jjoibbhdkruZBFtOyciI2kJalGsGLUhcIh+LH9SkHqgW3XCPSMuHElUcyqyfcKOmWnT9/yo0pY2XeI4PCKIlaXiU0BONEUrh291c9ZqibK0o6e1crnleZibBC6I8QdNupGk0TgYuEXCIhLUUhuko1/MhFsikUhkJCm8czFgmyjKVpWqB6QviTkmxERK4k5PIbi0Qk4Y80U7lT8uG3pECIrRuLs/kWZbUN+xjEhuIi7yhHXCEijaPZ5VfXYXPKa0SErdWki0yS6/F0SErSkNpCQqm0jp6RLtcqlMc4dgNxFHmU3FKx5NRsCDhO+dddlcV0S2qs5vTG05uKRW3SRFFmrnEEikPdRVbXNuOCJbeylAiKejLiSIStJTNdWcMRiUQIbpcyKx4fClw7SuH74KsZtUk4URtEdqr7IfldSLjlxSiVsuVF1cdqrYSEZCWkvWTx1rmmVu5PRUVTi2LpE4QxlyipPiN8QHBH0tKrmDki3eJYLhS7MkaJfgujt7RbuymHqwRKMbVXGa522JcqmKcBdKRd65IJEzEmokQyKMRJJp6LhecLdyqHraZz2SIiVolJWalbcKIlGO1Kp2bcakMpDHuqFrGeIJSKPL3VaSAoxiMUDhl4lLmigbUwcicK6VveSX8sEQIuUS3CW1WippyG1sh5UzUU0WnJe9ul6pK8b5LKrG64MrRtLsrMQG4Rjan8GRK71tyQ4z/VzLHKnCMEksVjgkK1i0RXRUyKOYVAjbp+skG7LcXiWgpvWRAU4j6SvRbOM3DpW8KJuUjbEh5rZJxhoRJFGqKVH+xWxEo7tMbVunpxEpfNTpEOmIp0RQq0jgCW3takrFhv0ttqXhilGQoEyA1+WtOhG2XLElXHOjgy0jy6VcRwWYEI7dqDV0OjzZahGPdRI5U2A2xlHlUjK7lTuONqiqiADKm9Vsk8xTi5tFS4YW95MutkJDH6qVNFxISIR8SDeopCPFIS26VYCiMiIZd1JAmnBiOoht2/7JyhSqml9jCTojJsS7qAfzdxz2yUY2CJDqVsrsncOIuCPDEhIhlqjtQdRkIkdoiIjpERtFPZ7VlrMitEhKKynJ/iCTZWiYy0l/FWPHJxEh82NvZ0rdS0I2jaJXW7kW6DKp9oYkUdMtNv7lHZvm7RNELUpFaOoRTFVlN0iIyHl/wB8UC/lwjcBR0xEtSJTS/RvE4kLvnAIZDzSQtaJcRwhG2RFp2pOWE6Lg3EVpCQp562W64kbOI4KkgK0oxuIRtUzl1YLoxIrtqbfZbiMh1CNyZo6Rsi9K0dKRpCsxITEmylZq3J/LXjfIG5DEtQko6oYclIit06tqapR4TgkOoSkMrpJEtGe5dTF7U3FxuIyG2Rbv1KGoq0RImilaVylMpB2pIi0iRSl8xRfSDInG3eK3G8okmV8Gs5zAWyhwrS7ur9ihsKbilYMbpEpbP6HhONXCREO1IGkJhonClcNvZVQtoiupBGRCUo2+kgctxk4IlLl1RRXswOBEriUZjjK7tLTELTncRFsRutu7Kq1Qd2r5yJZfIZSIiEuZDwuRPBMafIY8qW4+TkiLl2oulpxKJEMRSH6YWyt0kjtAHoiiVyN9huEQkIy8KaYYAtW0lbMmMY8ONqnKhE09K4LYyb901aoqVECIhG3VyiKk6kmB83plH0iW8QbbkIjpHUVxKBaiK/C2263Ttiqw7TXco8qttPEicIiIdvZVeraYZOkMpSJVEtvNNsUwiQyIilpQODAu6Yy7SHddMoiUrbRFGUNM6JDaQ7o7lV8QDMoy7zkXIp7OMti3IRuIkRT1PCujIrbbSTr58fs3SiplNFUeXibjYy73eTdRSPsHujtiVyk2aEWylKRS7SfaBx0huHl7qNlW8hN110iICjHUUbVYiZ97u7pJ9ukIW4y27bU2wyTZaiJG0aKHCI7vjIWrqOFcSNqCGMh5kJiyJWl3kFoxT5q2W26W5LzFoeE6XM07/4iW6hhoRGIxkJaVHOPSF1sSIvNOFdbbwiVY/abit7DVur1ljjHa9GSGnbuTYuXLMz5U/MQpIuRt5Uh2RaUtsPEno9l+VONpgjuitsmntNGgNy1VYR09nlWA4m3XJFpQcNgJSReA2pDRCn8UxszEk5hgP3uSTJbHFAZihqnEh0onAkPUFcPdSPYYCLcSJBxLFkUnBjteqkqUyTpJ3AiLcl8IUrhI0vYeqbIhIdMhig6WihG6XeipPHArk3FIi2Qlb/SnBbFNBilY6kaONu04lyxURW08bYjbptU5hpUfUEPeRYe0IDBOCRFzaU7gI6iGPZijTbKPL3U3gyfZU1SOfYArhbES7v0oU6Ps6lO4N7Y/f8AUtE3cNqW1RUqpktMtKVTU8buJpIZKeqqSUrVAVIEEhtS2eh7lYREI2x3ERRl+1M1OXjc58XtJlghIm2yukXKPL/mpXHHzgCOjdtVwj+TYcICkMRIdVqKzWvb4QkEijqlKIoPMHR/7Yjqlu7qjX3It3FIS0j9/wA6RWC8XhqZGTftYlES+aozpEbrjTYt9mQo7LcBKUrRiPZ/imaurbcLzYx4ZeL9qqUtKYeXuyEY6tIp+koCIiE7YqwhmLbRiTo6rh++CeewEyIrRb1dpV3thaV5yjLUI28qcao7RkMbhiRKVbwu4YjIR0ykPpfnUrV5bxBbuHulaP71FyVpA1jURERtlyonHLx4XMRDqjpJIIx4gjH3RFO1UbR0o2QWly62RRt7yfJsWyEhlFB19eXDiMu19/InKcjdaIhcjGNttyYGsm2TjfEGQiVpEVqX0mzARAeEMfOXRLaovF8iG/aQj3f2JvNnBiFuktwkKJ9oDMZiVxbSJEMuTErbvEh6k2jG60pStG2KVTk1FwgIpCMhFWbKCjbcIuIURG6Q3XKRqGXRNspRbjGVopTzDg0w1cR4ZRkPa5k1WZyxwxEW+IXDEZFtU+aTdQy6RlEJCPLzfQpLLKXSRCMtw3CQqP6N1RQc5iKXLuVkYASu3eJING03IRjcnqOpa02CUtRXKGqcwjU8MRLVqRANSckXzfkQFhCpkMYj4rkM88kUzkhIR+LFJNohK7UmDToEReknakOGMtRREoiSRhqEU+40fo7kJyBPZkJC3IY26iiIkSHfaHhvlESk07dqiPCL/BSr9M1GJDLlt0qLqq1tph9qMZNOiPatL96rH7Z1OODpTXDulci3orY4lFToNNGK07jFJI5JoTIijuQZomylIkQGCWDXMsLCIpWlI3gXxUjipnEiTzTZFFG2khbJ3IjBwkzwCHlSscO0n2k+2HJy443Vy6nfKlB2UhoR5hWCXdS7T9KcuH9p/wBw/hhEUyWqSSTnaHxJOJJ7ipyY36yh7B1Zi7JDFhJJxEvqob6FEt4Ykh8ESKoNJGOKdcHdJDjjJAbHBKx7SzydpMv1jbQlIvVRsQ7xCHTpQhDcJKN4hP8Ateke1FSdG0QxkpuUbahRFtTFQZRiKKdxTGOKyuQ0DpwKRERIpsrZS095axbSY7VG1RhSc5YqvZ7RlISEe1LmVh4Ue6scZEh1S7yNnpVadshiRCNtyIcqScKIqVfpx4Zjujaolqk4YiRDKIqpQFxfi25KRFIeVOu4iIyJBVrjZkLcYiOqMrUy8yZaZRVbhCcM7CJDG7tKPazHhlplK7upY5WW6PoxJPDSiPo8win3iATtWLo3Dut730JFwkJSKIxtuUvT5SQlKIxIvVU6zl7ZagEu9FFzCv1eYOSkAxEhiVpfIhzzh+0ZS27lZ84baYbblGJFbwxkQqLcbbIbLSLcQxTl2EKOJNkMtxFFGNjw4nIe0mCpSmIxulutFaecuJsbfrJELYwF8SHT4flUR5TadER5rvvgn2RJvml2U8zTxLicu0tyueAaqqmIx3SG7ak1L5EMtUkupxJwhKMY6kQZC42MR3XClske+0MQkUiIVmXNSKWkR1DzJ+sbGQFtHVbqR1fVtexOGLdxEMSj9CsCcwzYCo2mxGMbY+lyqrjdL0k81TERCJeipinyjuxuS3pOhORU48MeYlMY5gLcmxHTuiX0JDFIItCXKP3/ADKCdaIiIpRuIvRUbLSRYAXDJwi9VSYuiJXXSjG0lAs5O6PnPFd9C21ROC5xCujHUSZrXTMk3dJP44S2qvVGYkJNiREI7vvgjanMhIR4UijtjGSCHCHnOVEOHbEhl9VVahzDiEQ6oyLUUZI0c8bG1wSlGNskFUyLHDES1CQkonOMoadbcdjpafL1SWMZs1cJEVpR0p3MMzb4DsSKJNOCNpbhJXj9xll9JrAB0kUvCnMMNqZZw9JLLEZXeFSZuF2pJwC7vEicBHuprhDLUgm8MPvJaIVo2i2plxspKaez44CnsMUK0BS1IrDGKR7azHCTUdMiG4lEYt9pTwtk+QtNiJERWiXdR1P0dDh8WqdpaYSIhHju8CRDsHy/3rzPX8nXOR8P/wAmyy97CY/io8NNvNK3HkuVSiWZ0wf8wKmlfb/jhjgiT6PULTZOfhXL3WhtkFSzcXJHy6vIuPvlPOr/ANPmv/G1vy5u+2jOjAedc7nzv81bavo/lpShmtDHtVNPL+Kj6egYpnS4VSxUkQXcJwXYjzF5F2+m5LcpLK9n+G5M/k8cuxGLPa+Kti2nMCS/IvWfoxIgtGHaL1UsFs9XMg6ZPG0hQzIxRpghyw1CMUwa8lpXXIN9niDdtTgU5TlJEm2pyOItoOHpHVzKRHEiH+pJwZEkQ0AiMZEsbFygSAhK5acLso3EfSSIdn1UtNIBmSVjIkQ5h2UzjiQqVkYtkt4B2k5g4XZWzO3akZnVb8ZNvDbH4qUtJbLSIOhk7KPpCllTiMiEZFp3KVbHmituNDdd81PY0hKNkhK4bSHtKQClEtvqoyYjaI+JIwciSmxNhLIC3t226k7hcWlJF9ONODzJp0iMzy8nRttES8SZxoYgIy9X7/lU+RCkRHcrlVpCv05ORjISHuqFq6eR3D3oq3uBIUEdIMt3hVbQqFSJCMbrbttwrZOnGQ6R9JW38HiZcT1YrTeUt3CIld4VXYKydcZNR4fpbiQzDjjfNdzK6jlQbhitO5O32i7qNjUVZ/HiRiRRjzDqSacSIhlpH1lYyywR0y5dNqbwpx0pyppuioWyJuWovCpbGkJuIx1d5CBQEJSlHajAEoxIrtSaRVTTiTUbVEjl0R+ttR+BRHctOYi4MdMkGigfIbR+NqSWHCIStEZbil98Vs6ARkXEu5SJMv0/DbG+661EJvOuELcrpej8qh8tqxHSlV+J1LYiAxESu7SAbpX2ytVSEncBAS4sRa3FEYy/aiSbYdk6NwxjG3VzfkUHUU9S4IiQlzblM9HKMmxIS72pBVD+w3OLIbh7IyUpno8OmG7aXLyqYHHhlGIlqQWZ0guNuFKUWnbeWxPG+YzyWoTj/pTLuN0opb46YpvEVIKAxiNqUBXJrAU7whinojmGK3hgK0Ip/ARijRmcBW8W06lRQBeRNSrKYR3OfNJAf8RIcPLMv7VY/wD+AVJZD/8AmMf/ACfMJWvpV0Upc8pWGKo3wFh0ng4GLYkREMYlLDH8i8b1XqMOH1nFln/5ZK+O/m+XHj9ZxXP8rgfRutBro9WE7lxVOBv8EauTeAsOkAxl+SW3FFdSuWYvFmFTjiB4Cz7E4VpOiThCXsiOP5hw4eOHl/xLyLr2TdEKPL6Z3L2hdfpnycmNRwyIpBH+7Dyfk/IXooXIOgWX5K/7Mpn6wnYk3wjJgmHhIdJDgHlxHDHyY94U+X+W4csOfHHeNv1f3/8AHkcn8jxZY8uMvW36rh3Wlhi3nVSP5BjwAtERGQsNiVuH5PL5fKuxZ9lrTBsEACHEpGBKICIl5oStLD86RnHVvl9dVcd92q9kulxXRDFkW2h/Xh5fL5EZ0hfIxG4iaB0m2pcrdOI/wgieu4+e+nxwv+WP3/1HX/Hesw5eb0+ON/yx+0QOMdqXg4m8MFvy91ey+/OTiK21juSRLVpuS8ECmn6iO1C8a6UUus1JjFK1WJ3A7dMZLfliKawSscdqVpxsCTk0xgW1LxxUrkO4HbKJLBMdSZcxt9JD4ulplaotXoU6YkXzu0mSAUtl0YjctkMlCyMA7SbJq7V4RtRGDJc0STDsuaSSjeLfaW8Wu0sx73hW8S7JJFtvCn7SQ5T9pZxSTmGPEtuSgDYtjpScQT7jabxFUZrAUrH0U+2Eh2peFP3UyDinRwWEN2mSR5OylsqfKPKl4Rut9VMYkl4FK1XInRIEMS0pbQjqW8WdK2YxuTqdlxWcMZJImsF3mRoq2VPbGSjKinJsxG4rZd1SuDwoasOQkQ6k5SC1DxFaMhKKQ04Q2kMilqHSkOPRGREMkzhUkWkYq0jcT5vCo6ozRtsiHlLmu/3TjrxGQt8PTu3aUzWZY2TbjhEPEkPLJAFUzjcuIWnaNsiUfUsi6ZFKIbR1fKhWW3NJFLvff86O9jlHuilABEm2xIRIZW2pVO6RREh9IUw1QuOGJlddpEUbm4+xmvNlKQ7SG3veT8yrYbfzIScFsboxuUrGIi7HVaqRkoydIi3aVZmcwu4ZaeblQmwRWEIxIiIZeFA1Alw3b5Dw3C9UkdWUjbojcNt1xDcgatqInGNrTunTpJOfacp4WlzGS35VnkSmhuSSwRSsRT+AxSccFcBxrBOYprEloSQWyztW+JakuFLakHGI9rakaXyPD/1lMO0nPmkuiUbRFbcPaiS5LIoyHSOlOC+RWkXrEvL9b/Gz1Oczt11fNfzP8Nn6zlw5MbrrNOy1jODQTiUiiJYCN0u1/hgoioYICF1y50tAco/IuZE6W0il3iWG6XMXiXLn/C45f+55nJ/xbky8zKOhHSkMhEZOu6yjpHlJVrpY3HgCOgOK2JcxbiVe4hS3eIlmJq/S/wATODkmfbem38d/xzl9N6jDlyynXEvDDupJjas8topw9K9p9mQGKcwxTOGP3kkuupgjjCRLToR3LfAFbxEVNXCASCwuRMRjah3AKX9KSm22pJeDUVrBtYOEVNqpTD+FyaK1GWkmnGxlap02CeS5FtOEIimsQiRJ0OWSk23CMk2bRIkcB5iS4paCPwwJbxcLSjsQTLrUkaK0imak4Mk8LEf6VpkSEvRRjQkQp3FOwDwEgzUw60o19lTpcpkCinwdGKZglYYo0Vp8buytlj9xSME7gtJEZUjAE5hFKwaWQtRUGje07uZDVOZtCN2kdVpfQiTbQr2UtOykTl2q6IihNFUb7boybuFYdPIrVlDSNtDw29KKMREZFajake82Q8qGckI80kY7UbbUFVHp7qsILMHiF0bdQ9lO4EXDl8VJq5EekrVo8CERES7W5MhNC8Uh3bdUk/mGIuW6dOm3+CjssJ/jiMR3XR+lSr1GJSKREfLbFAVphwxcjcV24pKxg7Z2rUOzl5CV0R1apJ9qlKQkJae8QoJuR8JyOqNsR3KDpaV/hGJiRcRWOnfi7ErbpStTtSAkXDEtRJhRjYKmuHuqPCqKYykQyV6JkSImyG0dMoqBzLKhlIRL1UShuudIeEQCUdV1ykHH4sHzE278VBccm2AEolG0RKXyIcKgiaOXK7HVyon2mukYCtFglN4+ikP2kqYleXmSiIYofDV3VsnIpg8IrWAfe1M4mnQJBnMMdKYwwkRd6SU6RJA3eklQfwPUPMtYY3JOGNqTNJRw4kkEkzkscw0ig5G5cqS6cdSRhhciOHxB1DakZAOSj2U6HLJBt1DcybEht5UTiJfeKYOYtikQuWxxSsbkJZiIyWOxWxxWFjdpSXs0IrMATmIpl61TSlaKUkkkvAZdlIdwSafRylFPuNW/1JhnG30kW3hamvYDELitSeGjSwl/SmTC0lFPZjHvJeJDFMONlHupo8LUtDYyVtpSSTIhQrJFaiMUxttp1Etux3Ch8QSoXI2Wz4uDdIhTT+AlpJNvNxTTeCWjmRl3UttrbzREkYtkke9nScTmLqaitY4KtpFi8tYmmW4peJI2VOYkkYmIpJYoZ1FpCiqCTNSZGMe0hwJLxxS2ZLLBDcnsSHSXdimeIkcGRS1Ft7KNhlRTjLmEtqjswkLkY2jHlUuGBRimzp5K5kQBq7sEiwwiQxIS3Fama2lKI3fFFCsNOXXR5kwlMXBld2lEvVZNjG66UU6/U8O23l9JRbhOFqGXooTsaLQyEpWrWNYTdzd11sUUy95vTdaIiofMSIZEJR9FMxdJXjOTlvZJFykUmxkJFpUa1SiYyJz1YqVYNsR1S0jyoSX+DhcjMYj8VQua0otCYjylHwqSr3SGMdyh81YcdAiutEpR5YJ4/YROHWVX+80fgqPtEk+sivL3Kk8D32ixYvU9rH8fPe/n+tYdY1d71SeB77Rax6xa73qk8D32ixYj2sfwe/n+s/GJXe9Uvge+0Sx6x67D3Cj+De+0W1iPax/B8jP9ax6yK/3ql8D32iwOsivH3Kl8D32i0sS9rH8HyM/0rHrJr/eaP4N77RIx6xa73ql8D32ixYj2sfwfI5P1mHWLXe9Uvge+0W/xjV3vNL4HvtFixHtY/g+RyfrWPWNXe80vge+0TeHWBW3eap7uy99dYsR7WP4Pkcn6Fp+mNS2Uxappf5i9j89SH4yK73mk+De+0WliPax/D+RyfrePWRXe80vgqPtFr8Y9d71S/BvfaLFiPax/B8jP9bHrIrsPcqXwPfaJWPWVX+80fge+0WLEe1j+D5Gf6z8Zlf7xR/B1H2i1j1k13vFH8HUfaLFifs4fg+Rn+tfjKrveKPwVH2iTj1i13vVL4HvtFixHs4fh/J5f71sesauw9xpfA99onh6zswwwjwKPwVH2ixYj2cPwvlcv9mses7MPeKPwVH2iRj1k13vFJ8G99osWJezh+K+Vy/3pH4xa73mk+De+0WvxiVvvVL4HvtFixHs4fg+Vy/3rWHWFWe80vwb32i1j1hV3vVN4HfrrFiPZw/B8rl/vS8Osau96pfA99olYdZFd71S/BvfaLSxHs4fhfJ5f71hdY9dj7lS/BvfaJOHWLXSlwqXwPfaLFiPZw/D+Vy/3rPxiV3vVL4HvtFr8Ylb71S+B77RbWI9nD8HyuX+9Jx6wq33ml8D311mPWDWe9Uvge+usWI9jD8HyuX+9a/GDWe9Uvge+ulYdYld71S/BvfXWliPZw/C+Vy/3rf4xK33ql8D32iTj1g1nvVL4HvrrFiPZw/B8rl/vWv7fVnvVN4HfrrC6fVnvVN4HfrrFiXsYfh/K5f71r+3tZ71TeB7662PT6sH3Km8Dv11ixHsYfg+Vy/3rf4waz3ql8D311vHrBrPeqXwPfXWLE/Zw/B8rl/vTZ9PawhjwqbwvfXTGPTOq5GPC59daWJ+zh+F8rl/vSD6WVBY+UgZx/L5fzOfWTp9M6ohjw2PA55fjLFiPYw/B8nk/TYdLqkRjw2PAX0pT3S+oMYkzTl2ovS+OtLEezh+F8rk/sZw6T1H9wtYeif0rGulFQO1ou8JfStrEezh+H8nk/T/9sar+8GC7wufWW3OmdUQkMGMMCHEcYgeFviWLEvax/B8nk/X/2Q==\n",
"text/html": [
"\n",
" \n",
" "
],
"text/plain": [
" "
]
}
],
"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.7.3"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
\n",
"\n",
"
\n",
"\n",
"Following digestion and separation, assume the solid components of the red mud have a density of 4.0 g/cc, and that 35 vol% of the red mud is comprised of liquid retained in the void spaces. Based on this information, sketch an input/output diagram for the production of Al(OH)3(s) and estimate the minimum amount of water required to process 1000 kg/hr of bauxite ore. \n",
"\n",
"From this information, can you estimate the amount of caustic soda (NaOH) required to process the bauxite ore?"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Free Energy of Reaction\n",
"\n",
"The dissolution process involves reaction with hydroxide ions (OH-). Hydroxide ions are formed by the dissociation of water\n",
"\n",
"H2O (l) → H+ (aq) + OH- (aq)\n",
"\n",
"The dissociation constant at 25 °C is\n",
"\n",
"$K_{H_{2}O} = \\frac{a_{H^{+}}a_{OH^{-}}}{a_{H_{2}O}}=1.2\\times10^{-14}$\n",
"\n",
"The molar enthalpy and molar Gibbs free energy of formation of H+(aq) is normally taken to be zero at STP (1 atm, 25 °C). Using this information and additional data from Table B.3 of your textbook, estimate ΔĜ°f,OH- for OH- (aq)."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Estimating Equilibrium Constants\n",
"\n",
"The standard free energies of formation for selected compounds of aluminum at 25 °C. are given in the following table. Data is from [Parks, George, American Mineralogist, Vol. 57, pp. 1163-1189 (1972)](http://www.minsocam.org.proxy.library.nd.edu/ammin/AM57/AM57_1163.pdf). \n",
"\n",
"\n",
" \n",
"\n",
" Component wt % \n",
" \n",
"Al2O3 55% \n",
"Fe2O3 35% \n",
"SiO2 10% \n",
"\n",
"
\n",
"\n",
"In the Bayer process, alumina is dissolved by digesting Bauxite ore with sodium hydroxide (NaOH) solution at 175 °C. The digestion selectively dissolves the alumina by forming Al(OH)-4(aq) according to the reaction\n",
" \n",
"R1: Al2O3 (s) + 2 OH- (aq) + 3 H2O (l) → 2 Al(OH)-4(aq)\n",
"\n",
"(a) Using the results of problem 1 and the other available data, estimate the equilibrium constant KR1 of this reaction at standard pressure and temperature.\n",
"\n",
"(b) Because we don't have enough information to compute ΔĤ°rxn, use the Gibbs free energy as a crude estimate for ΔĤ°rxn. Estimate the equilibrium constant KR1 at 175 °C.\n",
"\n",
"(c) Write an equation for the relationship between equilibrium concentrations of OH- and Al(OH)-4 in the digester."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Estimate Equilibrium Yield\n",
"\n",
"Following digestion and separation of the undissolved solids, the liquid from the digester is cooled and seeded to precipitate Al(OH)3 by the reaction\n",
" \n",
"R2: Al(OH)4- (aq) → Al(OH)3 (s) + OH- (aq)\n",
" \n",
"The Al(OH)3 precipitate is a white, fluffy powder. \n",
"\n",
"(a) Estimate the standard molar free energy of reaction ΔĜ°R2 and the associated equilibrium constant KR2.\n",
"\n",
"(b) Write an equation for the equilibrium relationship between the liquid phase concentrations of Al(OH)4- and OH- in the crystallizer."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Flowsheet Analysis\n",
"\n",
"Prepare a mass balance of the Bayer process flowsheet shown in the accompanying figure assuming a basis of 1000 kg/hr of bauxite ore. \n",
"\n",
"![Ajka_alumina_plant.pn](https://raw.github.com/jckantor/CBE20255/master/images/Ajka_alumina_plant.png)\n",
"\n",
"Assume that the dissolution of Al2O3 reaches equilibrium in the digester at an operating temperature of 175 °C, and that the precipitation reaction reaches equilibrium in the crystallizer at an operating temperature of 25 °C. The spent solution from the crystallizer is mixed with makeup NaOH and water before recycling to the digester for reuse. The portion of the liquor from the digester lost through entrainment in the red mud is equal to 35% by volume of the red mud. The density of solids is 4 g/cc and the density of liquids is 1 g/cc. Adjust the flow of caustic soda to maximize the recovery of Al(OH)3(s). \n",
"\n",
"(a) What is the maximum recovery of Al(OH)3(s)? \n",
"\n",
"(b) What is the pH of the red mud when recovery of Al(OH)3(s) is maximized?"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### How much acid should be shipped to the accident site?\n",
"\n",
"There have been several proposals for remediating the environmental impacts of spilled red mud. At the time of the spill, one village poured vinegar into a river contaminated by the red mud. Suppose 10,000 cubic meters of sludge entered a creek. Assuming vinegar is approximately 5 wt% acetic acid, how much vinegar would be needed to completely neutralize the effects of the red mud? Roughly how many railroad tank cars would be required to transport the vinegar?"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Long Term Remediation\n",
"\n",
"There have been several proposals for a long term strategy to remediate the red mud from alumina facilities through [CO2 sequestration](http://www.epa.gov/climatechange/ccs/). Global alumina production is estimated to be 80,000,000 metric tons per year of Al2O3.\n",
"\n",
"Use the results of your calculations above to estimate how much CO2 would be required to neutralize the red mud produced as a by-product of alumina production through the reaction\n",
"\n",
"CO2 + 2 OH- → H2O + CO3-2 \n",
"\n",
"[Laboratory measurements](http://www.ncbi.nlm.nih.gov/pubmed/20036053) show that about 5.3 g of CO2 can be sequestered per 100 g of red mud. How does this compare to your estimate? Where would the necessary CO2 come from?"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"\n",
"< [West Virginia Chemical Spill](http://nbviewer.jupyter.org/github/jckantor/CBE20255/blob/master/notebooks/B.01-West-Virginia-Chemical-Spill.ipynb) | [Contents](toc.ipynb) |\n",
" \n",
"\n",
" Component ΔĜ°f kcal/gmol \n",
" \n",
" αAl2O3(s) -378.2 \n",
" Al(OH)-4(aq) -311.0 \n",
" αAl(OH)3 -275.3