{
"cells": [
{
"cell_type": "markdown",
"id": "f43cd6af-d772-43f7-a4da-30bcc41c7e50",
"metadata": {},
"source": [
"## Diffusion of a bell-shaped initial concentration along an initial gradient\n",
"A single chemical, whose initial concentration is a mix of a bell shape and a gradient. \n",
"Contrary to perhaps an intuition of a \"pile sliding down a sand dune as a unit\", the concentration peak\n",
"remains in place, and simply spreads out from there: reminiscent of an ice cream dropped on the floor!"
]
},
{
"cell_type": "markdown",
"id": "6398f935-72f9-458e-90ab-f02a05132302",
"metadata": {},
"source": [
"### TAGS : \"diffusion 1D\""
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "2b6f0de7-d66f-4f4e-8850-8a092fca20b8",
"metadata": {},
"outputs": [],
"source": [
"LAST_REVISED = \"May 3, 2025\"\n",
"LIFE123_VERSION = \"1.0.0rc3\" # Library version this experiment is based on"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "806647d0-4c0b-4abf-8d86-53460378d93e",
"metadata": {},
"outputs": [],
"source": [
"#import set_path # Using MyBinder? Uncomment this before running the next cell!"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "f6e0f263",
"metadata": {},
"outputs": [],
"source": [
"#import sys\n",
"#sys.path.append(\"C:/some_path/my_env_or_install\") # CHANGE to the folder containing your venv or libraries installation!\n",
"# NOTE: If any of the imports below can't find a module, uncomment the lines above, or try: import set_path \n",
"\n",
"from life123 import BioSim1D, ChemData, check_version"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "fe26e828-4428-4754-af9b-b54739b3a109",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"OK\n"
]
}
],
"source": [
"check_version(LIFE123_VERSION)"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "79045ff7-ed3f-4c3a-be79-9f8a4476f3de",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": 5,
"id": "a7c31adc-b1b2-419e-8a09-e049e6a95f5d",
"metadata": {},
"outputs": [],
"source": [
"# Initialize the system\n",
"chem_data = ChemData(names=\"A\", diffusion_rates=10.)\n",
"\n",
"bio = BioSim1D(n_bins=200, chem_data=chem_data)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "c984f746-d373-4a8a-89a4-9c3f0a3fca69",
"metadata": {},
"outputs": [],
"source": [
"# Set up the initial bell-shape concentration, with the peak close to one end of the system\n",
"bio.inject_bell_curve(chem_label=\"A\", mean=0.25, sd=0.1, amplitude=20., bias=0)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "b0064638-528f-40a6-9d07-8c4175af16e0",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
" \n",
" "
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"application/vnd.plotly.v1+json": {
"config": {
"plotlyServerURL": "https://plot.ly"
},
"data": [
{
"hovertemplate": "Bin number=%{x}
[A]=%{y}
System snapshot at time t=0"
},
"xaxis": {
"anchor": "y",
"autorange": true,
"domain": [
0,
1
],
"range": [
0,
199
],
"title": {
"text": "Bin number"
},
"type": "linear"
},
"yaxis": {
"anchor": "x",
"autorange": true,
"domain": [
0,
1
],
"range": [
-4.432342224959555,
84.21450227520528
],
"title": {
"text": "[A]"
},
"type": "linear"
}
}
},
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA2kAAAFoCAYAAADNS7MqAAAgAElEQVR4Xu3dCZgeVZ3v8X96zd5JZ09IGFxuHgUcVGTGkTE4LBIXRPQaGAUZdAyMKEQEXLiTyVx0WERAUIg6GkEGoleMiKIIQmZw5rncjHpFxuuCSEI2knTS6SSdXnP/501XW115l3rfc963TlV963nyJP121alTn3O68/7eU3XOuEO6CRsCCCCAAAIIIIAAAggggIAXAuMIaV60A5VAAAEEEEAAAQQQQAABBAoChDQ6AgIIIIAAAggggAACCCDgkQAhzaPGoCoIIIAAAggggAACCCCAACGNPoAAAggggAACCCCAAAIIeCRASPOoMagKAggggAACCCCAAAIIIEBIow8ggAACCCCAAAIIIIAAAh4JENI8agyqggACCCCAAAIIIIAAAggQ0ugDCCCAAAIIIIAAAggggIBHAoQ0jxqDqiCAAAIIIIAAAggggAAChDT6AAIIIIAAAggggAACCCDgkQAhzaPGoCoIIIAAAggggAACCCCAACGNPoAAAggggAACCCCAAAIIeCRASPOoMagKAggggAACCCCAAAIIIEBIow8ggAACCCCAAAIIIIAAAh4JENI8agyqggACCCCAAAIIIIAAAggQ0ugDCCCAAAIIIIAAAggggIBHAoQ0jxqDqiCAAAIIIIAAAggggAAChDT6AAIIIIAAAggggAACCCDgkQAhzaPGoCoIIIAAAggggAACCCCAACGNPoAAAggggAACCCCAAAIIeCRASPOoMagKAggggAACCCCAAAIIIEBIow8ggAACCCCAAAIIIIAAAh4JENI8agyqggACCCCAAAIIIIAAAggQ0ugDCCCAAAIIIIAAAggggIBHAoQ0jxqDqiCAAAIIIIAAAggggAAChDT6AAIIIIAAAggggAACCCDgkQAhzaPGoCoIIIAAAggggAACCCCAACGNPoAAAggggAACCCCAAAIIeCRASPOoMagKAggggAACCCCAAAIIIEBIow8ggAACCCCAAAIIIIAAAh4JENI8agyqggACCCCAAAIIIIAAAggQ0ugDCCCAAAIIIIAAAggggIBHAoQ0jxqDqiCAAAIIIIAAAggggAAChDT6AAIIIIAAAggggAACCCDgkQAhzaPGoCoIIIAAAggggAACCCCAACGNPoAAAggggAACCCCAAAIIeCRASPOoMagKAggggAACCCCAAAIIIJDKkLZi5e3y8PoNR7TehcvOlCsvOTcXrfrEk0/J8qtukmsuP1/OO/vUMdcc+Ky+4Qo5+aTjvfPYuHm7LH331XLc4mNk7eqV3tWPCiGAAAIIIIAAAgggkKRAqkPaQ/dcL4sWzCn4BaElL2/8CWlJ/thwbgQQQAABBBBAAAEE6ieQmZBmiG684z5Zs/YHRUeX6keYTMnlQloyNeKsCCCAAAIIIIAAAggg4EIgUyEtCC7BbY9BaDMjblf+453yy18/WzAL3yIYvXUyOhIXDkNP/uxXY26zjN5qeO+6R+XaW+4+ol2i+5Wrlzm4UhlBnaInOmPJiXLzqktHw2p4pNHsW+y46D7hup3/oU/Lzq7u0dNE9zXfWHLOZWP2mdnZIevvv7Vi3zz2lAslqG+4bsaqknOpwqupe3D+d7x5SeG2UbOF6x7HKqiHKSu8mT5k+lrYa9nyVbJtR1fBJrx/eJ+oZbHbd4vd6httlzj7VGwgdkAAAQQQQAABBBBITCAXIc3oFns+K3hT/PTja0YbwLxmtiBohN+sh8NWEAjC5ZqQtu6hJ8Y8ZxUEt2LHFqtX3DLKjaSFw0pwO2ix+pZ7zdStVMgIh5NwqA2eNQuHr1I9u1RIM/sXsyoWEKNlB9dTrO4mNIXbOQhKxW6PjWtV6npNICsW0oIPCcL1MHUt9oxesbKDcsv11zj7JPbbhhMjgAACCCCAAAIIxBLIVEgLRhCCN/TFwkqgUup7wZvjICiUC0Mm0M2d1Vlx8guz36uOf2lhlMts5epVqtWiZVQb0qKhKDhPeISnXN2CsBmE0krnrzSBS7mRtPBEKEF7xJkUplKbhsso5WEM4lpF7cr1rVL7mmNKfS8wD/pzqXqZ6w684+wT6zcDOyGAAAIIIIAAAggkJpCZkBa8oQ2P4pQLQ+XeNJs3usEb+jgTdERHRqK3v5nWDY/YxAlplcqoFJLMs3nBm/tq9i1Vt2hIC8KTubZaZpFsZEgzdYwG6lJhphqrcoEo7F8uiJULhUFdAt9g5LdcYI2zT2K/bTgxAggggAACCCCAQCyBVIe06BVGw0K5MBR9/idaVhD2qglpwUhe9Ha/aEAoV6+4ZVQTJqIBK3yt0Vv74oY0U0Y4qIXLjIbWYj3R15AW12rRgtmFZQSKBaZihqU+FChlGDYL3/5ZrN9G+32cfWL9dmAnBBBAAAEEEEAAgUQEUh3SKj2nVGkkzYhXWqerXBiKvvEuNbJSTUiLW0Y1Ia2afasJadEeGwScOMsgNDqkRc9nO5Jmrt1FSDPlhEduq/ktEA54pYJxnH2qOSf7IoAAAggggAACCNRfILchLc4th4a/XMAJv9Ev9+xU3JBWTRnRZ+fCXaXYtcV9zipuSDOBbOH82Ucsll3uNtJwHRsZ0oq1oYtn0ko9k1jNSJoxiWNm2nvtA48dsVh7+FyFDx0q7BNMJFP/Xy2cAQEEEEAAAQQQQKBWgdyGNANWbHZH87q55fCkV75MzAQW0Wn9A+i4s+gF54j7TFp0dslwPaMjVNWEhGJryJWbxTA6Shm9DTD4Ony7n4vZHaPLFdQycUj09j8TyKJ25UJaXKtiM3eGZwOtNDtm0JeKze4YfM+08d23faLwpRm5K9YHzPfMbKSlyinWp2r9hcFxCCCAAAIIIIAAAvUXyHVICwLZw+s3HCEd3D4WvOk2a2iF1wwrtR5Y9Hkg80bdrDcWngWy0ihenDJMhaPPM7leJy086lLsWa3wlPcBYJxZGM2+9RxJi7ZVsSUByoU0U7+466RF18YzIerEExYXFlUP34IYZ7Ss2LNk4VBWrE6l1vULd+g4t5/W/1cNZ0AAAQQQQAABBBCIK5DKkBb34lzsV+52RxflU4Y7gUrh192ZypdkRmJ/+tRvYy3q3ag6cR4EEEAAAQQQQACB9AgQ0iq0FSEtPZ250SHNjGSaUdJg4XMjRX9JT3+hpggggAACCCCAgK8ChDRCmq99s+p6JRHSzHNi0a2WdeOqvlgOQAABBBBAAAEEEMisACEts03LhSGAAAIIIIAAAggggEAaBQhpaWw16owAAggggAACCCCAAAKZFSCkZbZpuTAEEEAAAQQQQAABBBBIowAhLY2tRp0RQAABBBBAAAEEEEAgswKEtMw2LReGAAIIIIAAAggggAACaRQgpKWx1agzAggggAACCCCAAAIIZFaAkJbZpk3nhR17yoVyxpIT5eZVl6bzAqg1AggggAACCCCAAAKWAqkOaeYNfXQ7bvExsnb1SkuW4oc3eh2uulyE54XWI6TVo8yAsdqyyy12vWLl7fLw+g3i0zpr9VycO/rze83l58t5Z5/qeQ+leggggAACCCCAQP0FUhnS7l33qFx7y91HjLgEbyhndnbI+vtvda5HSHNOekSB1YaeODWqR5mEtDjyxffZuHm7mEXAwyOm9QyDtdeUIxFAAAEEEEAAgWQEUhnSzJvuciNmZkSiHrfLEdLq30nrEajqUWY9Qlr9das/Qz3CUzBi+PTja8ZUyLz+06d+W5cPWKq/co5AAAEEEEAAAQSSE0htSIvz3FIw4lbsNqrgzeeFy86UKy85t9ACwZvHcHM8dM/1smjBHAkCWrSpwmUHIwThfaK3ri1bvkq27eiST33s/bL8qptGdw1CZ1Dn4BtxrjPYN3qseT1cv/Ab7id/9qvCrXXBFjUqVla0PPN1OLie/6FPy86u7tEyA7uwx5JzLhuzT3TUMwhUJ73yZYXR0koOxvOXv3625H7FbomtNNJarB3NCaJtUW3ZgX+0DwXlFvsQIPCYP3emrFn7g9FDg34b7bOlbhmMuof7falfP5XqW+uvLVOXubM6j7gtOehzPt3uWes1chwCCCCAAAIIIGAjkMqQFrzhjPNmzuxrtujtj0FYCl4P3uyHP92PHltuJC14gxl+81ssJAbnCY8EhkNBsdfjvKEO6hY2Cd7AB2Ep/KY7/Ga+2HWZuq976Ikxb6SLXU84vIZDWTHP6AhosdveguAT51a4YiOqpV6rJuwG9Qr3hWJ1Nf2q2lG6ciNTpUKaOU+xflXq9XAfCOpdrF/FMak0klYsqBb7hRS2LGVW6Vw2v+g4FgEEEEAAAQQQSJNAKkNaqZGOYmGm2KfzwfHhoFLqjaN54xyMtJULaaWOj97CFQ2HQWcpNbpQ6vVoJytVrrn+173muMJoYKk3wYFHnDBo6vOq4186ejtpKZNomZXCSWBcyjHqUOq8xc5TbZAq9QNc7Ha8asuu5GBGy8Jht1j5pdqr2Ovl+oUZqSw22hm+ftfBqVTYNed0fa40/SKmrggggAACCCCAQFgglSEtfAHR293M94rdlhZ+LXiDHx01M7fqlQsq1QSDoI7Rc5V601zt69FuHIyalXtWr5aQVmkGzXLBNRyswsG63Aho3JBWysu4RMuoNkgFtsX6lvlenFGhUr9mXIS0YtcYnC/utQf1qDQa7To4EdL4DwgBBBBAAAEEEKgskPqQFr3E4I11+M1nNEiYN7LFwlj0uR1TdrlygnOXen4rXLdgxKLaMFYujJQKauHXw9dZTUgLQl808MYd0TJ1iO5bagQ0TuiJllVuhDF6y2O1IS1oz+iza8UmvKi27EaGtFLe4f5Radr7SiGN2x0r/5JlDwQQQAABBBBAoFqBzIW0UpOFBMHsqHmzChNSVLrNK/wGNwgRlUbSKo1KmMapZ0iLNn702b1qQlrcEa1abgGNhtvw6F/c89ZzJK1U2WkLaca51AcScX9RVAppccsJ78fEIbWocQwCCCCAAAII5EkgdSHNhKe1Dzw2+pxYtLFKhYbgeSIzq9z8uTPGTNFfqsxoWUEALBbw4o6o1CukFVt2IPoGO25IK/eMWtyRtOi5jN3C+bPl5JOOH9NkUY+4Ia1SYA6PEEWfo6v0A14qRBQLadWWXex5yKA+5WZ3jC4pUcop+no1I7HFXMrVt5Jjqe8zBX+tchyHAAIIIIAAAnkRSGVIMwvhmi06clVqkWuzb3hkLBqyis2AZ46Jzu5YLrwUm90xOK+Zmj48i6SZgr/SbJNBB4z7JjsYNQvfOlhqdsfoLW7FrqvYrJjBOcIjX8VmlTR1j95yWMyn1OyOxWYdLBaczDmKTeEffS6vVCioFCLC/SQ81X0x4+iaX+V+gZQKgfUIaaX6dtC/777tE4VJZWqpb62/JFnMulY5jkMAAQQQQACBvAikLqQFDVPs+THzvXK3MZZ7jqnYmlDFJuGIPn9WaZ00U6fwc2H1GkkLglG044bDQ9yRtFLGxtYEzvAaV0GwMGEpvEZasaBVbK256LOBcUfSwiG23Dppxa6l0jpp5phi648Fa8tFA1m4L8YpO/qsWJx10modSSv381JukplwPypVX9tfktHn2So9H2d7Po5HAAEEEEAAAQTSIpDakFYtMAvlVisWb/9yz6TFK4G9EEAAAQQQQAABBBBAICyQm5BWalFruoOdACHNzo+jEUAAAQQQQAABBBCICuQipJWa8ZHuYC9ASLM3pAQEEEAAAQQQQAABBHI5kkazI4AAAggggAACCCCAAAJpEMjFSFoaGoI6IoAAAggggAACCCCAAAJGgJBGP0AAAQQQQAABBBBAAAEEPBIgpHnUGFQFAQQQQAABBBBAAAEEECCk0QcQQAABBBBAAAEEEEAAAY8ECGkeNQZVQQABBBBAAAEEEEAAAQQIafQBBBBAAAEEEEAAAQQQQMAjAUKaR41BVRBAAAEEEEAAAQQQQAABQhp9AAEEEEAAAQQQQAABBBDwSICQ5lFjUBUEEEAAAQQQQAABBBBAgJBGH0AAAQQQQAABBBBAAAEEPBIgpHnUGFQFAQQQQAABBBBAAAEEECCk0QcQQAABBBBAAAEEEEAAAY8ECGkeNQZVQQABBBBAAAEEEEAAAQQIafQBBBBAAAEEEEAAAQQQQMAjAUKaR41BVRBAAAEEEEAAAQQQQAABQhp9AAEEEEAAAQQQQAABBBDwSICQ5lFjUBUEEEAAAQQQQAABBBBAgJBGH0AAAQQQQAABBBBAAAEEPBIgpHnUGFQFAQQQQAABBBBAAAEEECCk0QcQQAABBBBAAAEEEEAAAY8ECGkeNQZVQQABBBBAAAEEEEAAAQQIafQBBBBAAAEEEEAAAQQQQMAjAUKaR41BVRBAAAEEEEAAAQQQQAABQhp9AAEEEEAAAQQQQAABBBDwSICQ5qAxtuzqdVAKRSCQnMCE9mYZ39osu/f1J1cJzoyAA4HWliaZNqlVdnT3OSiNIhBITmDcOJG50yfI1i7eYyTXCpzZRmD+jAk2h+f+WEKagy5ASHOASBGJChDSEuXn5A4FCGkOMSkqUQFCWqL8nNyBACHNDpGQZudXOJqQ5gCRIhIVIKQlys/JHQoQ0hxiUlSiAoS0RPk5uQMBQpodIiHNzo+Q5sCPIpIXIKQl3wbUwI0AIc2NI6UkL0BIS74NqIGdACHNzo+Qpn5PPPmULL/qplHJC5edKVdecu7o1ytW3i4Pr99Q+Pq4xcfI2tUrx6gzkmbXCTk6eQFCWvJtQA3cCBDS3DhSSvIChLTk24Aa2AkQ0uz8ch/SNm7eLkvffbWsvuEKOfmk4yX69b3rHpU773pA1t9/a0F62fJVcuIJi8eEOEKaXSfk6OQFCGnJtwE1cCNASHPjSCnJCxDSkm8DamAnQEiz88t9SAtG0Z5+fM2o5JJzLpOLLzhLzjv71CNCWTS0mYMIaXadkKOTFyCkJd8G1MCNACHNjSOlJC9ASEu+DaiBnQAhzc4v9yEtGB375a+flYfuuV5H0l6QT1735dGRs3BgM/sWC3WENLtOyNHJCxDSkm8DauBGgJDmxpFSkhcgpCXfBtTAToCQZudHSFO/G++4Tx780X/Izq7ugmb4mbRjT7lQrrn8/MKoWjikmUC3aMGcwmvDhw7ZtQJHI1BngSEt/xe9vfKfB3plQ+HPAfll70EZCJ13UlOTvKitTV7c3iYvam+XF7e16t/6dVu7HK1/t9S5jmksXpcxYvNIoG9gWNpbmzyqEVVBwFbA/JbhPUYcRZTiKDV2nybzSQNbzQK5D2lmZCw8cmYkTTALglqckbStLGZdcwfkwPoJ/Htfr3z/wD75mf79s367hX3H6y/av5owSd40cbIsHT9JJmigY+Otk499gJE0H1uFOtUiwEhadWrEgeq8GrH3PBaztmLOfUgzo2gbfv7rMTM2mslB5s+dITevupRn0qy6Fwc3WqB7eEi+ub9H7urplmcG+8ecfm5zs7yibby8qn28vFL/PkH/njzucNgytzv2NY2TDd375LnBAdk4MCC/17+f0zKeHRyUnUODY8p616SpcmnHdHlxS1ujL5HzIVBWgJBGB8mKACEtKy2Z3+vgdke7ts99SAueMYvO7hjc4sjsjnYdjKMbI7Bbw9nnu3fLV/ftkYMjt99Ob2qWs3Tk6xQdAXuV3rI4s7n0DYuVnkkzwe2R3v3yr/rnJzoy1ztyjrdo+R+e2inHavlsCPggQEjzoRWogwsBQpoLRcpIUoCQZqef+5Bm+EwQu/aWu0clw8+gmRdZJ82uk3F0/QR6Dg3L6r17ZHXPbjkwPFw40TkTp8g7J0+RJXpbYtytUkgLl2MC4Zd69shX9Lzm/GY7VYPgZRrWXq2jc2wIJClASEtSn3O7FCCkudSkrCQECGl26oQ0O7/C0czu6ACRIqoS6NVA9rX93XKbjp7t0dBktqUalK6eNlNe2lr9LYjVhLSgovs0oK3RsHanhjUT3Mz2uvYJsqJjhrx2/ISqroedEXAlQEhzJUk5SQsQ0pJuAc5vK0BIsxMkpNn5EdIc+FFEdQJf0lB0696u0WB0io6YfXxapxynz5nVutUS0oJzmRE8c5vlF/U5uODZtbdNmiIrO2bKnBbmhKy1TTiuNgFCWm1uHOWfACHNvzahRtUJENKq84ruTUiz8yOkOfCjiHgCvx3olw/v3Ca/GDg8U+OrNJR9QkfOXIxa2YS0oPbmWbgv622X/7RnV+ElMynJVdNmyPumTIt3geyFgAMBQpoDRIrwQoCQ5kUzUAkLAUKaBZ4eSkiz8yOkOfCjiMoCd+zdLdfu2VnYcY5OCHLdjDlyht7e6GpzEdKCuvxBZ4a8qusFnWDkQOGl1+hzanfOnC9mdkk2BOotQEirtzDlN0qAkNYoac5TLwFCmp0sIc3Oj5DmwI8iSguYZ70+sGOrmDXPzHaRjkpdraNTwdT5ruxchrSgTt/cv1c+290lG3VmyE4Nll+YOVf+cvxEV1WmHASKChDS6BhZESCkZaUl83sdhDS7tiek2fkR0hz4UURxgd/o7Y0X7Ngsm3SdsmNaW+W2GXML65vVY6tHSDP1NFP1f1xH1UxgM9sVHZ3yEZ1YhA2BegkQ0uolS7mNFiCkNVqc87kWIKTZiRLS7PwIaQ78KOJIgX87eEDepyNo+3UGRfPM2VdnzZcpIwtP18OrXiEtqOs9+7rl73fvKKzhZq5ntQbOGWXWbavHNVJmPgQIaflo5zxcJSEtD62c7WskpNm1LyHNzo+Q5sCPIsYK3KmTb3xq904xK5AtmzRVbtTnz+r9NFe9Q5q5wl/phCfvfWGLbB4alNka0L40c56cyLpqdH/HAoQ0x6AUl5gAIS0xek7sSICQZgdJSLPzI6Q58KOIPwpcvmv76K2BK6fNkg9MbczMiI0IaeYq9+p0/R/auVUe0ZFCs/29zk65fOp0ugACzgQIac4oKShhAUJawg3A6a0FCGl2hIQ0Oz9CmgM/ihDZpaNLF2l42dB3sHBb4xdnzZPXN3CSjUaFtKCtv6CzVX5qZLbK9+tkKKumz6IbIOBEgJDmhJFCPBAgpHnQCFTBSoCQZsXHFPx2fIeP3rLr8Mx7bAjUImDWP/vrFzbLFg1qR+niz1+ftUBe2tpWS1E1H9PokGYq+t0D++RiDaZmO2viZLlDb39kQ8BWgJBmK8jxvggQ0nxpCepRqwAhrVa5w8cxkmbnVziakOYAMadFmJEzM4Njt94GaJ7PWqMThEzX6eobvSUR0sw1mrXU/uaFwxOk/EW7TpAye77z5QUabcn5khUgpCXrz9ndCRDS3FlSUjIChDQ7d0KanR8hzYFfXov4Ye9+uWjHlsLlnzJ+knxZb3GcYP5XTmBLKqSZS32q/6C8Wx12DQ3JSRrU/mX2gsQcEqDnlI4FCGmOQSkuMQFCWmL0nNiRACHNDpKQZudHSHPgl8cizLphZpIQs50zcYrcpgs9J7klGdLMdZs14d62bZPs1RG1N0yYpLd8zk+Sg3OnWICQluLGo+pjBAhpdIi0CxDS7FqQkGbnR0hz4Je3In6kI2gXjoygmdkbzSyOSW9JhzRz/T/TEbV36bN5B/TWT55RS7pHpPf8hLT0th01HytASKNHpF2AkGbXgoQ0Oz9CmgO/PBVhgsg7tz9fWNT5XboG2s26BpoPmw8hzTj8e1+vvFuDWr9nPj60EXWIJ0BIi+fEXv4LENL8byNqWF6AkGbXQwhpdn6ENAd+eSnimcF+eave0mcmCTlNp9f/qj571eTJxfsS0gzHIzrro1mOYEj/zfT8nnSQFFWDkJaixqKqZQUIaXSQtAsQ0uxakJBm50dIc+CXhyK26aQYb972nJi/zSyO35x9lLQlNElIMW+fQpqp37r9PXLprm1ySP99VccMuayjMw/dhGt0IEBIc4BIEV4IENK8aAYqYSFASLPA00MJaXZ+hDQHflkvYp9OhvGmrZvEjKS9vK1dvj3nKO+mmfctpJk+cde+bvl41wuF7vFPnbPlgskdWe8qXJ8DAUKaA0SK8EKAkOZFM1AJCwFCmgUeIc0OLziaddLcOGaxFPNs1TJ9xupJfdbq6JZW+d7chYmsg1bJ1seQZup8c/cu+Ux3l5iFCW6fMVfOnjSl0qXw/ZwLENJy3gEydPmEtAw1Zk4vhZBm1/CMpNn5FY4mpDlAzGARw3pN79NZHB/W2RxnNrfIQxrQ5uvfPm6+hjRjtXL3Dvlyzx4xS3x/ZeY8OW3iZB8JqZMnAoQ0TxqCalgLENKsCSkgYQFCml0DENLs/AhpDvyyWsRH9Va9e/WWvY6mpsItjotb2729VJ9DmkH74M5tsu5AT+E5vnt0wpW/0EWv2RAoJkBIo19kRYCQlpWWzO91ENLs2p6QZudHSHPgl8UibtHb9G7U2/TG6/+y35qzUE7QZ9F83nwPaWamx4v0ttFHDh6QCWr6v1Jg6nN7Z7luhLQst26+ro2Qlq/2zuLVEtLsWpWQZudHSHPgl7Ui7tPRsytGJrz42qz5ctqESd5fou8hLQD8aw1q6zWoTR3XJN/V20df0trmvS0VbKwAIa2x3pytfgKEtPrZUnJjBAhpds6ENDs/QpoDvywVYdb4eq+u8WW2W3Wh6nfqgtVp2NIS0np1jbn37tgqP+k7UJiI5UEdUetsNk+rsSFwWICQRk/IigAhLSstmd/rIKTZtT0hzc6PkObALytF/Kz/oLxFF6s228pps+QDU6el5tLSEtIC0LO3b5L/03ewsObcdzSosSEQCBDS6AtZESCkZaUl83sdhDS7tiek2fkR0hz4ZaGI5wcH5U3bN8ouXaz6kinT5ZrpM1N1WWkLabuHzeLgm+S5wQF5+8QpcvvMuanyprL1EyCk1c+WkhsrQEhrrDdncy9ASLMzJaTZ+RHSHPilvYgeXaz6LA0Mvxnol1P1+bO79Dm0tG1pC2nG93fqbYKaWSz8qo4ZcllHZ9rYqW8dBAhpdUClyEQECGmJsHNShwKENDtMQpqdH69+85sAACAASURBVCHNgV/aiwgms3i5zuD4ndlHyUSdcj9tWxpDmjH+N51E5FydTMRs9+nU/H85fmLa6KmvYwFCmmNQiktMgJCWGD0ndiRASLODJKTZ+RHSHPiluYiP6yyOd+lsjrN1kervebxYdSXjtIY0c1236lIHN+iSB51NzfLwvEUyz9MFwyu1Ad93I0BIc+NIKckLENKSbwNqYCdASLPzI6TZ+RHSHPiltYh1+3vkg7u2FdZC+7aOoL1CJ7FI65bmkGbMz9+xRX7cu19erW3wABOJpLUbOqk3Ic0JI4V4IEBI86ARqIKVACHNik8IaSN+9657VK695e7CVzM7O2T9/beOyq5Yebs8vH5D4evjFh8ja1evHKO+ZVevXStwdOoEnh0YkNO3PSe9hw7JP3XOlgsmd6TuGsIVTntI69aJRM7YtlHMBC7vnzJNVk2fler2oPK1CxDSarfjSL8ECGl+tQe1qV6AkFa9WfgIQppq3HjHfbJm7Q/k6cfXHKFpwtuddz0wGtqWLV8lJ56wWK685NzRfQlpdp0wbUcf1GC2VAOBmSjkjTpRyFdSOFFI1DztIc1cz3/198lbdGr+Pm2fNdomp6dgEfG09f001JeQloZWoo5xBAhpcZTYx2cBQppd6xDS1O/YUy6Uh+65XhYtmHOEZjSURUObOYCQZtcJ03b0lfoc2r/oc2gLW1rkkXlHy+Rx6ZsoJIshzVyTaRfTPpO0TX44d5Ec09qatu5FfS0FCGmWgBzujQAhzZumoCI1ChDSaoQbOSz3Ie2JJ5+ST173ZdnZ1T0qecaSE+XmVZcWvl5yzmVy8QVnyXlnn1r42uy//Kqbxoy6EdLsOmGajv7ugX1y8c6tYt76f09DwLE6o2MWtiyMpAXt8OGd2+RbB3rkpa1thaDWbt7psOVGgJCWm6bO/IUS0jLfxJm/QEKaXRPnPqQFz6KFb3U0I2sXLjuzcEuj+fc1l59/REgLj7zt7um3awWOToXA7/v75fWbnpUDejvddTNnywemZWddrrbWJmltbpL9BwdT0RblKmludzxt0x/kab398V2Tp8qdc9O3bl3qGyGhCxjWtjchbVJ7i+w9MJBQLTgtAo4E9POlaZPaZM8+3mM4EqWYBgtMn9LW4DNm63SEtMgzZ6Z5zUQhZjOjaXFG0nr7h7LVK7iaIwQG9M3fnz/zjPyqT597mjJFvrFoUaaUmpvGifnTPzicieva2D8gf/7738meoWH5/Pz58jfTp2fiuriI8gJDQ4dDWmvLOOkbyEZfps3zK2DuARjf1iy8x8hvH0j7lU/Q/stWu0DuQ1qx2xfDIY1n0mrvXFk68mP6nNPdGXsOLdw+WbrdMbiu9Qf3y1+/sKVwu+ODOi2/WWycLfsC3O6Y/TbOyxVyu2NeWjq718ntjnZtm/uQZvjMaNlbTn9t4fbGjZu3y9J3Xy2rb7hCTj7peGF2R7sOloWjf6jrb12k63C16MV8P0PPoWU9pJnru14Xuf6cLnZ9lE7y8uO5R8ukpvRP8pKFn6l6XgMhrZ66lN1IAUJaI7U5Vz0ECGl2qoS0ET/z7FmwhZ9BM6+xTppdJ0vz0ZsGB+S0rRtl36FhWTltlnxg6rQ0X07JumdxJM1crLnhbdn25+Xf+3oLU/KbqfnZsi1ASMt2++bp6ghpeWrtbF4rIc2uXQlpdn6Fo5nd0QGih0WY59DepOtumfW33qBv8L+e4Tf4WQ1pplvt1oWuT9/6nGwdGpL/MX2mXDyF59M8/HFzViVCmjNKCkpYgJCWcANwemsBQpodISHNzo+Q5sDP1yKu2f2CfLWnW+Y1N8ujuh5aR1N2H4DNckgz/esXfQflLA3ch/Tf9+vzaa9uH+9rt6NelgKENEtADvdGgJDmTVNQkRoFCGk1wo0cRkiz8yOkOfDzsYjwc2jf0efQTsj4pBNZD2mmj63RwP1JDd6zmlvksXmLZHqGQ7ePP1ONqhMhrVHSnKfeAoS0egtTfr0FCGl2woQ0Oz9CmgM/34oIP4d2zbSZcsnU7N8el4eQZvrZcl2I/EFdkPwv2ifI2jlHCdOI+PbTZ18fQpq9ISX4IUBI86MdqEXtAoS02u3MkYQ0Oz9CmgM/n4owSzkv3bYxF8+hhd3zEtJ6h4dlqd72+NuBfvlwR6dc3THDp+5HXRwIENIcIFKEFwKENC+agUpYCBDSLPAIaXZ4wdFMHOLG0YdS/mHPDvnS3j2FW+LW6y1xWX4OLY8hzVzzswMD8kYN4gd0xs57Zs+XJeMn+dD1qIMjAUKaI0iKSVyAkJZ4E1ABSwFCmh0gI2l2foWjCWkOED0o4hG9De69ejuc2b6tt8KdpLfE5WXLy0ha0J7f17b+W23rafpc2o80jM/XUM6WDQFCWjbakavQW53GicydPkG2dvXCgUAqBQhpds1GSLPzI6Q58POhiM2Dg3L6tuekW2+H+9i0GfKhqZ0+VKthdchbSDOwK3fvkC/37JFX6qQwD+rkMGzZECCkZaMduQpCGn0g/QKENLs2JKTZ+RHSHPglXYR5Du1tevvbz3U9tNeNTCihH2DmastjSOvXdfCWbtsk/2+gr7BIuVmsnC39AoS09LchV3BYgJE0ekLaBQhpdi1ISLPzI6Q58Eu6iGt375Q7enbnemr2PIY00+9+rxOInK4B/aAGNrNYuVm0nC3dAoS0dLcftf+jACGN3pB2AUKaXQsS0uz8CGkO/JIs4rHe/fKeHVsKU7Hfr8+hvSZHz6GF3fMa0oxBsH5aR1NTYdHyeTyfluSPpPW5CWnWhBTgiQAhzZOGoBo1CxDSaqYrHEhIs/MjpDnwS6qIrUODcurWw8+hfVSnY1+R4+nY8xzSTP8794XN8m8HD8jL9fm0789ZKK3m3RFbKgUIaalsNipdRICQRrdIuwAhza4FCWl2foQ0B35JFfH27c/Lk329+hzaRPnGnAVJVcOL8+Y9pHUNDRUmjtmmf/+tPp/2Dzyf5kW/rKUShLRa1DjGRwFCmo+tQp2qESCkVaN15L6ENDs/QpoDvySKuG7PLrltb5dM1ynYH9Mp2M26aHne8h7STNv/Hw3s52hwH9Z/83xaen8aCGnpbTtqPlaAkEaPSLsAIc2uBQlpdn6ENAd+jS7iJ30HZNn2zYXTrtURNDOSlveNkHa4B9zcvUs+092li5jzfFpafyYIaWltOeodFSCk0SfSLkBIs2tBQpqdHyHNgV8ji9ihz6G9YetG2T08JJfpc2hX5fg5tLA7Ie2wxiH98/btm3RU7aCcoM+nfUfXT8v3GGsjfzrdnIuQ5saRUpIXIKQl3wbUwE6AkGbnR0iz8yOkOfBrVBHmNrZzRt6Av6Z9vM7muLAwqyObCCHtj70gHOT/bup0+eS0mXSRFAkQ0lLUWFS1rAAhjQ6SdgFCml0LEtLs/AhpDvwaVcSNeivbLXorG8+hHSlOSBtrEtwSa0bWzKQy3BLbqJ9S+/MQ0uwNKcEPAUKaH+1ALWoXIKTVbmeOJKTZ+RHSHPg1ogieQyuvTEg70ucGDfW3Euob8ePp9ByENKecFJagACEtQXxO7USAkGbHSEiz8yOkOfCrdxHh29cundopH582o96nTF35hLQjm4zbY1PXjQsVJqSls92o9ZEChDR6RdoFCGl2LUhIs/MjpDnwq2cRTAQRT5eQVtwpHPA/rBPNXM1EM/E6VIJ7EdISxOfUTgUIaU45KSwBAUKaHTohzc6PkObAr55FMKV6PF1CWmknc6vsu3TJhnG6C0s2xOtPSe5FSEtSn3O7FCCkudSkrCQECGl26oQ0Oz9CmgO/ehXB4sTxZQlp5a0+vWenfH7vbiadid+lEtuTkJYYPSd2LEBIcwxKcQ0XIKTZkRPS7PwIaQ786lGEWQfNrIdmbldbPmW6/P10plEv50xIK98LB/Xbb9u2UX7e3ydm+YZv6/INZmSNzT8BQpp/bUKNahMgpNXmxlH+CBDS7NqCkGbnR0hz4Oe6CPMc2rLtz8tP+npZkDgmLiGtMtRWDfynbn1OuoeH5aP6fNoKnk+rjJbAHoS0BNA5ZV0ECGl1YaXQBgoQ0uywCWl2foQ0B36ui/icTpt+vU6f3tHUJI/OO1rmNbe4PkXmyiOkxWvSx3r3y3t2bCksgn7/nKN0VG1CvAPZq2EChLSGUXOiOgsQ0uoMTPF1FyCk2RET0uz8CGkO/FwWwXNotWkS0uK7/cOeHfKlvXtklob/x+YtKjynxuaPACHNn7agJnYChDQ7P45OXoCQZtcGiYe0Y0+5sOorOG7xMbJ29cqqj6vXAVt29daraMqtQiD8HNr7p0yTVdNnVXF0vnclpMVvf/N82lJ9Pu2/9Pm01+lI2lodUeP5tPh+9d6TkFZvYcpvlAAhrVHSnKdeAoQ0O1kvQtrTj6+JfRX3rntU1j30BCEttlh+dnzXyHNoL29rl+/rxA6t5n84tlgChLRYTKM7bRockNN0Ypp9h4YLi6ObRdLZ/BAgpPnRDtTCXoCQZm9ICckKENLs/Alpdn6FoxlJc4BoWcQXdHr0T+k06ZPHNckjegvawpZWyxLzdTghrfr2/qE+n3aRPp9mnnj8ztxFhUlq2JIXIKQl3wbUwI0AIc2NI6UkJ0BIs7MnpNn5EdIc+NkWYaZFN9Ojm9vQvjJrvrxxwiTbInN3PCGttia/ZvcL8tWebp2cprkwSU0Hz6fVBunwKEKaQ0yKSlSAkJYoPyd3IEBIs0NMPKTZVd+PoxlJS64dunU9NDMt+tahIXnv5A75dOfs5CqT4jMT0mprvIFDh+RN2zcVnk97g3448HX9kIAtWQFCWrL+nN2dACHNnSUlJSNASLNzT0VICyYXqebZtVpZlpxzmcyd1TnmmbcVK2+Xh9dvKBRZbNISQlqt2vbHmenQzbToPIdmZ0lIq90v/Hzaymkz5QNTp9deGEdaCxDSrAkpwBMBQponDUE1ahYgpNVMVzjQ25BmJgi59pa7R6/uoXuul0UL5thdbYWjTUAzWzikmXrcedcDsv7+WwvfW7Z8lZx4wmK58pJzR0sjpNW1WUoWvkZvM/uk3m7Gc2j2/oQ0O0OeT7Pzc3k0Ic2lJmUlKUBIS1Kfc7sQIKTZKXoX0sKjVubSVt9whZx80vF2VxnjaBO+zl56sjy/dYds+PmvR0fSoqEsGtpM0YS0GMCOd/nNQH9hGvSDervZbTPnyjkTpzg+Q76KI6TZt/fHurbL3fv26qQ1LTp5zdGFDw/YGi9ASGu8OWesjwAhrT6ulNo4AUKanbU3IS28Xto1l58v5519qpjXGhHSwkHsxjvuGxPSzOjaxRecVaiP2Z548ilZftVNEr71kpBm1wmrPbpfg9nSbZvk/w30yZsnTpYvzpxXbRHsHxEgpNl3CfOBgfngwHyAYCavMZPYsDVegJDWeHPOWB8BQlp9XCm1cQKENDtrL0JaqWfOGhHSzMid2W5edWnh72hIM3UIQmM4pIVvv+wbGLJrBY6uSuCKrdvk87t2ycLWVvnpS18iU5oYsagKsMjOTfpuoKlpnAwODdsWlevjn+nrl9f87ndyQAPbzfPmyiUzZuTao9EXPzh0SExIa2keJwOD9OVG+3M+9wJtLc3SP8h7DPeylNgIgfbW5kacJrPn8CakzezsGH3uK9BuREgzI2U7u7qPaOCgPnFG0nZ292W2g/h2Yet7D8g7t24SE8u+v2CRvLp9gm9VTGV92tuapK25SXp6zUIGbDYC39nfI+/fvkXMSn0PL/gTOa6d9dNsPKs59pDu3KoBbfKEFtmzb6CaQ9kXAe8EzEha55R22bWX9xjeNQ4ViiUws4P//2JBldjJi5Bm6ha+3TG4xbERIS3qEh1J45k0m+7l9tgtQ4OybPvz8vvBAbmqY4Zc1tHp9gQ5Lo3bHd02/hX6fNp9+nzan+ii6j/UxdV5Ps2tb7nSuN2xcdacqb4C3O5YX19Kr78AtzvaGXsT0oLLSGrikOD80ZDG7I52Hczl0Uu3bpRf6HNoJ4+fKGtnL3BZdO7LIqS57QK9I8+n/VafT+O5Sbe2lUojpFUS4vtpESCkpaWlqGcpAUKaXd/wLqQFlxOegr/YrZB2l1366GhIM3uyTlq9tOOX+5Fd22Xt/r2yoLmlMDIxvYn7nOPrVd6TkFbZqNo9wjOQXj99trxnSke1RbB/DQKEtBrQOMRLAUKal81CpaoQIKRVgVVkV29DWriujVzMuhZOZnesRS3+MV/b1y2f6HpB2vV/rO/NXSgva+Ue5/h68fYkpMVzqnav+7TvXqF9d7z23YfmLpL/1tpWbRHsX6UAIa1KMHb3VoCQ5m3TULGYAoS0mFAldks8pJkAFp7OvtLlmBG2dQ89MbqOWaX9G/F9Qlr9lP+z76C8Y/smMVMA3K7rob2d9dDqgk1IqwtrodCLd26V7x7YJ8fobKQ/mrNIJjAbaf2wtWRCWl15KbyBAoS0BmJzqroIENLsWAlpdn6FowlpDhCLFLFtaEjeqOtO7dQJQ94/ZZqsmj6rPieiVCGk1a8T9A4Py+nbN8qzAwPyVl3X707W9asfNiGtrrYU3lgBQlpjvTmbewFCmp0pIc3Oj5DmwK9YEQM68cI5OpPjT/sPyp+PnyDfmH2U8BRanbC1WEJa/WxNyeHn0z7bOUeWTZ5a3xPmuHRG0nLc+Bm7dEJaxho0h5dDSLNrdC9CWrWXcNziY7jdsVq0lO1/lT7Hc48+zzO3uVkemXc0E4XUuf0IaXUG1uLv1v78sZHn0x6Ze3Th9kc29wKENPemlJiMACEtGXfO6k6AkGZnmXhIs6u+H0dzu6Pbdvh6T7dcvfvwRCHfmXOUHN823u0JKO0IAUJaYzrFRTu2yA979xcmEDETiZgJRdjcChDS3HpSWnIChLTk7DmzGwFCmp0jIc3Or3A0Ic0B4kgRTBTizrKakghp1WjVvu++Q8Ny2tbnZNPgoJw3uUM+0zm79sI4sqgAIY2OkRUBQlpWWjK/10FIs2t7QpqdHyHNgV9QhJkgxEwUYiYMuUgnCvmfTBTiULd8UYS0hlHL0/198mbt52bGUjOJiJlMhM2dACHNnSUlJStASEvWn7PbCxDS7AwJaXZ+hDQHfkERb9eJQp7s65U/a58g9+ttjmyNEyCkNc7anOmrPXvkmt07ZOq4Jnlcn7mc09LS2Apk+GyEtAw3bs4ujZCWswbP4OUS0uwalZBm50dIc+BnijCLVZtFq+fpRCEP66QKnfo3W+MECGmNsw7O9H59Pu2hkefTvj9nIeunOWoCQpojSIpJXICQlngTUAFLAUKaHSAhzc6PkObAb+2+vfKRru1i5rr77tyFTBTiwLTaIghp1YrZ779f109bum2TPDPYL2/SWx6/xPpp9qhaAiHNCSOFeCBASPOgEaiClQAhzYpPCGl2foQ0S7+ndB20t+ltjn26LhrrR1liWhxOSLPAszj0ucEBOX3rRtmvE4p8ctpM+bup0y1K41AjQEijH2RFgJCWlZbM73UQ0uzanpBm50dIs/ALTxTyXp3p7tPMdGehaXcoIc3Oz+boH+stjxforY9mMv5v6qLtZvF2ttoFCGm123GkXwKENL/ag9pUL0BIq94sfAQhzc6PkGbhd46OoP1vnSjkJJ0o5NtMFGIhaX8oIc3e0KaEW7u75IbuXYUFru+auUBexELXNXMS0mqm40DPBAhpnjUI1alagJBWNdmYAwhpdn6EtBr9/l5ntvtnneFurk4Q8gNd1HdWM7Pb1Ujp5DBCmhNGq0Iu2blVHjiwT/6kpVW+qxOJMHlObZyEtNrcOMo/AUKaf21CjaoTIKRV5xXdm5Bm50dIq8Fv3f4e+ahOFNKrz6E9qBOFvLJtfA2lcIhLAUKaS83ay3qXji7/REeX/7StXb6ltz5OaGqqvbCcHklIy2nDZ/CyCWkZbNScXRIhza7BCWl2foS0Kv029JmJQjYVjrplxhz575OmVlkCu9dDgJBWD9Xqy9yrMz6an4/fDPTLaeMnytdmL6i+kJwfQUjLeQfI0OUT0jLUmDm9FEKaXcMT0uz8CGlV+G3SmezO1CnH9wwPyQd1FrtP6Gx2bH4IENL8aAdTi61Dg4Wp+Xfo3+/RDzGu1w8z2OILENLiW7Gn3wKENL/bh9pVFiCkVTYqtwchzc6PkBbTr0enGF+6baM8OzBQGCFYoyMEZjY7Nj8ECGl+tENQi18N9Mnbtz0v5ufm6o4Z8uGOTr8q6HFtCGkeNw5Vq0qAkFYVFzt7KEBIs2sUQpqdHyEtpt/FOinCd3VShMKzNjopwgTzvw+bNwKENG+aYrQiTxw8IMte2Fz4+raZc+WciVP8q6SHNSKkedgoVKkmAUJaTWwc5JEAIc2uMQhpdn6EtBh+QUA7qqWlMGvdbGZyjKHW2F0IaY31jnu2+w/0yId2bivs/o05C+R17RPjHprb/QhpuW36zF04IS1zTZq7CyKk2TU5Ic3Oj5BWwW/Fru3yjf17ZaYGswd0LbSjdXpxNv8ECGn+tUlQo9v2dsl1e3bJlHFN8u25R8nLWtv9rawHNSOkedAIVMGJACHNCSOFJChASLPDJ6TZ+RHSyvit1LXQvqxroU3WN5cP6JvLxby5dNDb6lMEIa0+rq5KvVo/7Pi6fthh1hN8SJetmMdodElaQpqrXkc5SQsQ0pJuAc5vK0BIsxMkpNn5EdJK+N3a3SU3dO+S8fq/zFpd7+nEdtZCc9DV6lYEIa1utE4KHtZSLtyxRR7t3S8vaWmTB+ctLIyssR0pQEijV2RFgJCWlZbM73UQ0uzanpBm50dIK+J3175u+XjXC9Ki37tr9nxZMn6SA2WKqKcAIa2eum7KNou/v10Xu36q/6D8+fgJct+sBdLKBDxH4BLS3PQ3SklegJCWfBtQAzsBQpqdHyHNzo+QFvFbt79HLt11eKKD1TPnyZsnTnYgTBH1FiCk1VvYTfldQ0OydPtGeX5wUM7Sn6079GeMbawAIY0ekRUBQlpWWjK/10FIs2t7QpqdHyEt5PeITrF/kU61P6Sv/VPnbLlgcocDXYpohAAhrRHKbs5h1ho8a/sm6WJR+KKghDQ3/YxSkhcgpCXfBtTAToCQZudHSLPzI6SN+P17X6+8W9d06tdbsj6qC++u0AV42dIjQEhLT1uZmppbHs/WWx8P6s/b9dNny3um8IFI0IKEtHT1ZWpbWoCQRu9IuwAhza4FCWl2foQ0FdjQd1DOfeF5Mc/MnK+jZ9fpKBpbugQIaelqL1NbM4mImUzEbGtmzZdTJ/Dsp7EgpKWvL1Pj4gKENHpG2gUIaXYtSEiz88t9SHu6v0/eqZ/o7z00LO+cNFVunTHHgShFNFqAkNZocTfnu0cn6blKJ+kxm/lwxHxIkveNkJb3HpCd6yekZact83olhDS7liek2fnlOqT9fqBf3qoBbY8+G3Pa+InytdkLHGhSRBIChLQk1N2c8zO63MXNutyF2T6nH5K8Qz8syfNGSMtz62fr2glp2WrPPF4NIc2u1Qlpdn65DWlbhgblrds2yjadbc6sgfZNXQutjenAHfSmZIogpCXj7uqs1+x+Qb7a0y3NWuAX9dbHM3N86yMhzVWvopykBQhpSbcA57cVIKTZCRLS1G/JOZfJzq7uUcnVN1whJ590/OjXK1beLg+v31D4+rjFx8ja1SvHqG/Z1WvXCik7eqcGtLfpCNofBgfk+Lbx8r/mLJDJLKybslYcW11CWqqbr1B5szahWaPQbOvmHCWvaZ+Q/ouq4QoIaTWgcYiXAoQ0L5uFSlUhQEirAqvIrrkPaRs3b5ebv/hNuXnVpQWee9c9Ktfecrc8/fia0a/vvOsBWX//rYWvly1fJSeesFiuvOTcUc48hbRuvbXxHds3y68G+uRFLa3ynTkLpbPZfH7PlmYBQlqaW++Pdf/Iru2ydv/ewgtmncK35HCdQkJaNvoyVyFCSKMXpF2AkGbXgrkPaVE+E9qWvvtqeeie62XRgjlHhDIT4sKhzRyfl5C2f3hYlutC1Y/prHLzmlvkgbkLZb7+zZZ+AUJa+tswuAIzkYiZUMRseVyvkJCWnb6c9yshpOW9B6T/+glpdm1ISIv4PfHkU7L8qptGR9LMrZAXX3CWnHf2qYU9o983r23ffdCuFVJw9G4dQTt322b5v7o+0/SmZnlw3kJ5cWtbCmpOFeMIjG9rkvaWZuk+MBBnd/bxXOCGPTvls3u6CrVcMa1Trp420/Mau6neITkkrc1NMnViq3T19LsplFIQSEjAhLRZHePlhT3Zf4+REDGnrbPAnOnj63yGbBdPSIu0rwllbzn9taO3Mx57yoVyzeXnHxHSgpE2c/jQ8KFM95Ln+gfkjGd+L8/098vclhZ59CUvksXt7Zm+5rxd3Dh9N2DeEAxnvC/nqV2/2tUlH9i0WWOLyPs7O+ULCxdIU8YBBgaHpa21WZq0L2f993LGm5LLGxFo1s5MX6Y7pFXA9F+22gUIaSE7E9BedfxLR59PM9+KM5KW5dsdzSyOS7dtEjNZyEta2uRenSSEWxxr/4Hz9Uhud/S1Zezq9b0D++QDO7cWClmqMz5+WWd+zPrG7Y5Zb+H8XB+3O+anrbN6pdzuaNeyhLQRv2IBzXwrOlFInp5Je+LgAfnbHVsLC1W/RqfZv2vWApnalPXP4u1+oNJ6NCEtrS1Xud7/cbBX3rtji+zXn+M/0xkf75o9P9OzsRLSKvcJ9kiHACEtHe1ELUsLENLsegchTf3MLY0XLjtzzIyNAWs0lOVldkczlbeZ0ttsb5s0Rb4wY65dT+NorwUIaV43j3Xl/qu/T5a9sFm69NnSl7e2y9d1RG2O3rqcxY2QlsVWzec1EdLy2e5ZumpCml1r5j6kBVPuRxnPWHLi6G2PeVonbUgh/l7D2RoNaeZO4is6OmVFxwy7XsbR3gsQ0rxvIusKPj84qEHt8PqG5pblb+gC9Me0tlqX61sBhDTfWoT6i6B85wAAGd1JREFU1CpASKtVjuN8ESCk2bVE7kOaHd/ho7PyTFqP3g71N3pblLk9arz+7/AFXWfpjfocC1v2BQhp2W9jc4VmllYzova0jqyZWVrv01sfj9MF6bO0EdKy1Jr5vhZCWr7bPwtXT0iza0VCmp1fZkLaHwYGCs+t/G6wX+bom7ev6Zu34zP25s1BU2e2CEJaZpv2iAvr1fUO36eTiazXZ04n6LvA+3VE7RX6zGlWNkJaVlqS6yCk0QfSLkBIs2tBQpqdXyZCmhk5e7++adujn7K/Qp9XWZPh51UcNHcmiyCkZbJZy17UpTu3ybcP9BT2+UznbDlvckcmEAhpmWhGLkIFCGl0g7QLENLsWpCQZueX+pC2dt9euapruwzqlZhbG80tjuZWR7Z8CRDS8tXewdXe1N0lN3fvKqyl9oGp0+R/TJuV+rXUCGn57MtZvGpCWhZbNV/XREiza29Cmp1fqkPaP+7eKat7dheuwUwQ8hEmCHHQG9JZBCEtne3motZmLbXLdm2T3kOH5IS2drlOR9XSfKszIc1Fr6AMHwQIaT60AnWwESCk2ejpaPoh3eyK4Og0ThxygT5/9mjv/sKo2ed0ev03T5xMQ+ZYgJCW48bXSzcTiVywY7NsGzLzu4pcM22mXDJ1eipRCGmpbDYqXUSAkEa3SLsAIc2uBQlpdn6pG0n7vU4Q8r6dW+Q3A4cnCDEL22ZtdjcHTZq7IghpuWvyIy64SwPaCh1Re0QnFDHbX46fKLfNmCOzdLr+NG2EtDS1FnUtJ0BIo3+kXYCQZteChDQ7v1SFtAf298hHdQ20/TrV/kntE+ROHUHL6oK2Dpo1V0UQ0nLV3GUv9mu6RuInRhay79QPcm7W2x9PS9FIOyGNvpwVAUJaVloyv9dBSLNre0KanV8qQtozOq3+/9Tnz36ktzea7YN6G9Mn9HYmNgQCAUIafSEs8Fsdab9EZ3/81UBf4eX3TZkm/zh9ViqQCGmpaCYqGUOAkBYDiV28FiCk2TUPIc3Oz/uQdovO3HajzuBmtmn6qbi5femvWKDaQatnqwhCWrba09XVrNy9Q77cs6dQ3Mt0eY7Vs+bKi1vaXBVfl3IIaXVhpdAEBAhpCaBzSqcChDQ7TkKanZ+3Ic18Av6hndtHPwk/VYPZZ/W2pZkpe77EQfNQRAwBQloMpJzu8oQ+o3axjqrt1nUUzURDK3VE7QKP11QjpOW0o2bwsglpGWzUnF0SIc2uwQlpdn7ehbR+nazzszpydsfersLaZyaUXatvqt6aomdKHDQJRVQpQEirEixnu+8cGpQP7touJrCZzXzo8zkdlTej875thDTfWoT61CpASKtVjuN8ESCk2bUEIc3Oz6uQtqHvoFzetU2e1RkczfbOSVNl1fSZXr6RcsBOEQ4FCGkOMTNc1Bf37pFVe3YUrvC/tbbJx/XZ1jM8u32akJbhDpizSyOk5azBM3i5hDS7RiWk2fl5EdIODA/Ldfrs2T+PPDsyX0fPbtJPuV+vU2izIRBHgJAWR4l9jIBZU+2DOlW/mVzEbGdqSPtU5xyZ2+zHqBohjX6aFQFCWlZaMr/XQUiza3tCmp1f4iHt3/T2oyv0NqTNejuS2f5mSod8vGOmTGpqcnBlFJEXAUJaXlra3XXeprdUX7dnV6FAM1X/R6fNkHfp6P0E884ywY2QliA+p3YqQEhzyklhCQgQ0uzQCWl2fomFtB0ayu7QW49W9+wu1OFFLa1ys657dmL7eAdXRBF5EyCk5a3F3VzvxsEB+YTOAPnYyPIeZuHrD+sSHxfplP1JbYS0pOQ5r2sBQpprUcprtAAhzU6ckGbn1/CQ1qsTg3xx7275vP4xi1K3aA3+bmqnrOjolLaEP8F2QEkRCQkQ0hKCz8hpf6whzdxybW6FNNs8E9b0d9J5OrLW2uDfS4S0jHQqLkMIaXSCtAsQ0uxakJBm59fQkHbvvm75oj539pvQsyCf0IlBfF+3yAExRdRZgJBWZ+CcFH//gR65SW+B/IOOsJntqJYWWaEfIp3bwCn7CWk56Ww5uExCWg4aOeOXSEiza2BCmp1f3UNaj46WrdFg9s893WJucTTbn7a1y7W65tmr2ri10UHzUYQKENLoBq4EzG+pe/QDpVt0KZAXRn5nHdPaWghrb9eRtXo/LUtIc9WSlJO0ACEt6Rbg/LYChDQ7QUKanV/dQtr2QX3mTJ83u2f/XjGzN5rtde0T5BJ93uMNnk157YCQIhIWIKQl3AAZPL25Nfsr+jvs9u7dslc/bDKbGfX/dOcsObmOM88S0jLYmXJ6SYS0nDZ8hi6bkGbXmIQ0Oz/nIe2Zwf7ChCDm1sZge4suRG2eOzMjaGwI1EOAkFYPVco0AruHhwpB7c6RSY7May9vbZfL9Jk187vN9UZIcy1KeUkJENKSkue8rgQIaXaShDQ7P2ch7Wf9B+VOnQzkwQP7Rmt0/uSpOnLWKUfrzI1sCNRTgJBWT13KNgLbhobkJp1c5F9CH0CZWWlXTZ8lf6Z3CbhaNoSQRn/LigAhLSstmd/rIKTZtT0hzc7PKqTt1dsYv68P2j+ma50F4WzKuCZ5r05f/be63tlMnSGNDYFGCBDSGqHMOYxAl4a1f9nfrc+t7RUzhX+wmVG1sydOkaWWo2uENPpZVgQIaVlpyfxeByHNru0JaXZ+VYe0/RrMfqHTVJv1zX40sraQKcRMWb18ynR5j46eTWAhagetQhHVCBDSqtFiX1cC/6ofUH1Tw9r3e/fJQX2GzWyT9IOqN5vANmmyLBk/qepTEdKqJuMATwUIaZ42DNWKLUBIi01VdEdCmp1frJC2RWc4+97+ffLIwf3yhL4pCW9mMpDzdeTsrZafHju4DIrIsQAhLceN78Glm8mRvqdB7Zv7euQnfX/8HTm9qVneVghsU+Q1+rsyzkZIi6PEPmkQIKSloZWoYzkBQppd/yCk2fmVDGn/2XdQHtGRsh/qG49fj6xrFpzKzHD2xomTdKHXDnmRTk3NhkDSAoS0pFuA8wcCZtr++3R07RsH9sqzA3+8HXKB3m1wjk7hbz7QOrbMJEqENPpSVgQIaVlpyfxeByHNru0JaXZ+oyGtVz8J/te+3kIwe0gn/zAzmoW31+uU03+lt+6YcLaIiUAcqFOESwFCmktNynIl8Av9sOvbvT3ybb0TIVgn0pT9Ev2g6+Kp0+R4DWvHRdaLJKS50qecpAUIaUm3AOe3FSCk2QkS0uz85Eu7uuS+Hbvlx6Hny0yRnXqbzum6ntkZ+qnv6/U2nYk8Z2YpzeH1FCCk1VOXsl0ImN+x39rfI+t0sqXwZn7XnqEffr1DR9mm6fNsfzpxgkyb1Co7uvtcnJYyEEhMgJCWGD0ndiRASLODJKTZ+cm4n/1itITj9RPdUydMlNM0nL0y8umu5Wk4HIG6ChDS6spL4Q4FzORLZjZccyv5Tw72yr6RhbKDU0zUd7avmThRjm1uk1fo7+ETdLRtIXcvOGwBimqUACGtUdKcp14ChDQ7WUKanZ98cNNmedFwcyGYzWLKfEtNDk9KgJCWlDzntRX43xrUzGQj/9nXJz/r75VuDXHRbZqOtr26/XBgMx+gvVL/bV5jQ8BnAUKaz61D3eIIENLiKJXeh5Bm51c4esuuXgelUAQCyQkQ0pKz58xuBTbLkPzm0IA8sWefhraD8n/1TzC9f/hM5tlgE9rO1tskp+jt6OZWyWNa22SCeWfMhoAHAoQ0DxqBKlgJENKs+ISQFsNvxcrb5eH1Gwp7Hrf4GFm7euWYowhpMRDZxWsBQprXzUPlqhAoNnHI07o25S80rD2lf/+8ENyKP682VYPase3thbNN0X+b9Svnt7QU/l7Q3Cpz9e8/YUbeKlqDXW0ECGk2ehzrgwAhza4VCGkV/O5d96jcedcDsv7+Wwt7Llu+Sk48YbFcecm5o0cS0uw6IUcnL0BIS74NqIEbgbizO5plUswo2290iZQ/DA7Ic/pno/6Js83UsFYIcM3NfwxyOjI3v+lwqGMG3ziK7FNJgJBWSYjv+y5ASLNrIUJaBb9oKIuGNnM4Ic2uE3J08gKEtOTbgBq4EYgb0kqdbYuu07Z1cFA2a2Dbqv8ufD3y2pahAdk2NHZ5lWLlmOeTXxIacfvT1vFFZ/g1q2SamX8n6bvxiTpyZ/49sWmc3nKpf+trk0b+Nt+bwAzBbjpIikohpKWosahqUQFCml3HIKRV8FtyzmVy8QVnyXlnn1rY84knn5LlV90kTz++ZvRIQppdJ+To5AUIacm3ATVwI2Ab0uLUYrOGuG0a3EyQ2zJ8ONQdDnf6tf57e2SdzDhlxtmnEOI0sAWhzjw/d/i1wyFv0si/zd9NMvbZuiY5VHjNvKo5sPDvwt+HpPDq4dfM90b2Ma8Hr418//Cxuq9+z+w7zvzbvBa8bl4Lvg6/VjjPyL4j5z5ctilL66X/HlOH4HyRYw7vd7jOR9Q3fD1aeBaeLSSkxfmpYB+fBQhpdq1DSKvgd+wpF8o1l59/REh76J7rZdGCOXb6HI0AAgggkEmBXTri1jWofzS87Rsalr36tVk+oNgf8/1yr3fFGL3LJCIXhQACqRY49MpXpLr+SVeekFahBeKMpCXdiJwfAQQQQCD7AntNyNNAt2807I0Ev1DIM987oH+GdLQpvA3pSNqwvmYWKBge/Xex1/T7Oro1dr/wayKHRr8/9vXR4wrnGSkj2HfMa5E6mDqV/H419T1cn179w4YAAskLENLs2oCQVsGPZ9LsOhhHp0OA2x3T0U7UsrJAI253rFwL9kDAXoDbHe0NKSFZAW53tPMnpFXwY3ZHuw7G0ekQIKSlo52oZWUBQlplI/ZIhwAhLR3tRC1LCxDS7HoHIS2GH+ukxUBil1QLENJS3XxUPiRASKM7ZEWAkJaVlszvdRDS7NqekGbnVzia2R0dIFJEogKEtET5OblDAUKaQ0yKSlSAkJYoPyd3IEBIs0MkpNn5EdIc+FFE8gKEtOTbgBq4ESCkuXGklOQFCGnJtwE1sBMgpNn5EdLs/AhpDvwoInkBQlrybUAN3AgQ0tw4UkryAoS05NuAGtgJENLs/Ahpdn6ENAd+FJG8ACEt+TagBm4ECGluHCkleQFCWvJtQA3sBAhpdn6ENDs/QpoDP4pIXoCQlnwbUAM3AoQ0N46UkrwAIS35NqAGdgKENDs/QpqdH0cjgAACCCCAAAIIIIAAAk4FCGlOOSkMAQQQQAABBBBAAAEEELATIKTZ+XE0AggggAACCCCAAAIIIOBUgJDmlJPCEEAAAQQQQAABBBBAAAE7AUJajX4rVt4uD6/fUDj6uMXHyNrVK2ssicMQaKzAsadceMQJV99whZx80vGF1+nbjW0Pzla9wLLlq+TEExbLlZecO+bge9c9Ktfecvfoa08/vmbM95ecc5ns7OouvHbhsjOPOL76mnAEAnYCN95xn2z4+a+PeA/B72k7V45GIAsChLQaWtG8Ebjzrgdk/f23Fo4u9YahhqI5BIG6C5j//MOhLHxC+nbd+TmBhUD4A4RoyNq4ebssfffV8tA918uiBXMk+ubXHGu2m1ddWvi73M+BRRU5FIFYAuEPFIp90Mvv6ViM7IRApgUIaTU0bzSURd/Y1lAkhyDQMIFy//nTtxvWDJzIQsCMiL3l9NeOGQmLhrJoaIv2+2hos6gOhyJQs4Dph1u27So6klbqwzR+T9fMzYEIpEqAkFZDc5k3CBdfcJacd/aphaOfePIpWX7VTRK9taaGojkEgboLRG+jCb8RoG/XnZ8TOBAoFtKKha4gmC1aMHvMKJupQqnbzBxUjyIQiC1QLqSFC+H3dGxSdkQgMwKEtBqa0vzHf83l5x8R0oLbbGookkMQSETAvFFds/YHox8w0LcTaQZOWqVAsZBmRhfmz50xejujKTLozwvnzy58kBb+HW36/oM/+o/R29arrAK7I+BEoFRICxfO72kn1BSCQOoECGk1NBmjDTWgcYi3AuHbwOjb3jYTFQsJMJJGd8iKQJyQFnzgEIym8Xs6K63PdSBQXoCQVkMP4X7wGtA4xFuBcEijb3vbTFSsQkjjmTS6SBoFaglp/J5OY0tTZwSqFyCkVW8mzIBXAxqHeCFg+q7ZgucpzRuEnz7129FbvujbXjQTlaggUGwkjdkd6TZpFCgW0vg9ncaWpM4IuBcgpNVoylpSNcJxWKICwSQ34UpEJ7yhbyfaRJy8jEC4bwa7hfsv66TRfdIiEO2rpt7B7Yz8nk5LK1JPBOorQEirry+lI4AAAggggAACCCCAAAJVCRDSquJiZwQQQAABBBBAAAEEEECgvgKEtPr6UjoCCCCAAAIIIIAAAgggUJUAIa0qLnZGAAEEEEAAAQQQQAABBOorQEirry+lI4AAAggggAACCCCAAAJVCRDSquJiZwQQQAABBBBAAAEEEECgvgKEtPr6UjoCCCCAAAIIIIAAAgggUJUAIa0qLnZGAAEEEEAAAQQQQAABBOorQEirry+lI4AAAggggAACCCCAAAJVCRDSquJiZwQQQAABBBBAAAEEEECgvgKEtPr6UjoCCCCAAAIIIIAAAgggUJUAIa0qLnZGAAEEEEAAAQQQQAABBOorQEirry+lI4AAAggggAACCCCAAAJVCRDSquJiZwQQQAABBBBAAAEEEECgvgKEtPr6UjoCCCCAAAIIIIAAAgggUJUAIa0qLnZGAAEEEEAAAQQQQAABBOorQEirry+lI4AAAggggAACCCCAAAJVCRDSquJiZwQQQACBPAjceMd98uCP/kPW339rHi6Xa0QAAQQQ8EyAkOZZg1AdBBBAwGeBJedcJju7uo+o4tOPrxl9bcXK2+WnT/021QGHkOZzL6RuCCCAQPYFCGnZb2OuEAEEEHAmYELaW05/rVx5ybmjZS5bvkq27ehKdSiLAhHSnHUZCkIAAQQQqEGAkFYDGocggAACeRUoFtKigabU1xdfcJZce8vdo3Th0bdinib8zZ87o/Cth9dvKPx93OJjZO3qlYV/b9y8XZa++2pZfcMVcvJJx48WcewpF8o1l58v5519auG14OvwuS9cdqYsO+sNheODLXxMcA0mkK5Z+4OSdTb7hb9frIzwdYe/n9c+xHUjgAACCFQWIKRVNmIPBBBAAIERgWIhLfpasZBmgswZS06Um1ddWijJBDCzBYGrGLDZ55e/fvaIwGUClhnJqyakmfKDUHjvukcLYXFmZ8fo6F/wWrBPEL6CcxWrc/Q6o/UJyghfNx0JAQQQQACBOAKEtDhK7IMAAgggUBAo9UxasRGkYNKNYrcOmlB0510PlL1FMhhJC4KdOb953s1s5rVqQlq4fsWOC1576J7rZdGCOVKszk88+ZQsv+omCfYxI3TBv4PuEa4ft0zyQ4MAAgggUKsAIa1WOY5DAAEEcihQbCQtCDjBqFOl2x8NW3TkqhhlqZC2ZduuwghcPUJacOtksYAVPp+prwlsxbbglkxCWg5/QLhkBBBAwJEAIc0RJMUggAACeRAoFtLMdZsRpCA85SmkRUfSwn2AkJaHnwiuEQEEEKiPACGtPq6UigACCGRSwKeQZoDNLYdxJw4JJhIpd7tjuZG04HbH4Lm16AQl0QYnpGXyR4CLQgABBBoiQEhrCDMnQQABBLIhUCykBeEleO6rUSNpRtTU51XHv3R0QhIzomdmgiw2u6NtSCt1rvBomrH41vfWF+pDSMtGn+cqEEAAgSQECGlJqHNOBBBAIKUCpSYOCY9mNTKkBaNiAacJZ2bmRlchLbxwd7FZGqNT8Jt6lBuNS2mzU20EEEAAgQYLENIaDM7pEEAAAQQQQAABBBBAAIFyAoQ0+gcCCCCAAAIIIIAAAggg4JEAIc2jxqAqCCCAAAIIIIAAAggggAAhjT6AAAIIIIAAAggggAACCHgkQEjzqDGoCgIIIIAAAggggAACCCBASKMPIIAAAggggAACCCCAAAIeCRDSPGoMqoIAAggggAACCCCAAAIIENLoAwgggAACCCCAAAIIIICARwKENI8ag6oggAACCCCAAAIIIIAAAoQ0+gACCCCAAAIIIIAAAggg4JEAIc2jxqAqCCCAAAIIIIAAAggggAAhjT6AAAIIIIAAAggggAACCHgkQEjzqDGoCgIIIIAAAggggAACCCBASKMPIIAAAggggAACCCCAAAIeCRDSPGoMqoIAAggggAACCCCAAAIIENLoAwgggAACCCCAAAIIIICARwKENI8ag6oggAACCCCAAAIIIIAAAoQ0+gACCCCAAAIIIIAAAggg4JEAIc2jxqAqCCCAAAIIIIAAAggggAAhjT6AAAIIIIAAAggggAACCHgkQEjzqDGoCgIIIIAAAggggAACCCBASKMPIIAAAggggAACCCCAAAIeCRDSPGoMqoIAAggggAACCCCAAAIIENLoAwgggAACCCCAAAIIIICARwKENI8ag6oggAACCCCAAAIIIIAAAoQ0+gACCCCAAAIIIIAAAggg4JEAIc2jxqAqCCCAAAIIIIAAAggggAAhjT6AAAIIIIAAAggggAACCHgkQEjzqDGoCgIIIIAAAggggAACCCBASKMPIIAAAggggAACCCCAAAIeCRDSPGoMqoIAAggggAACCCCAAAIIENLoAwgggAACCCCAAAIIIICARwKENI8ag6oggAACCCCAAAIIIIAAAoQ0+gACCCCAAAIIIIAAAggg4JEAIc2jxqAqCCCAAAIIIIAAAggggAAhjT6AAAIIIIAAAggggAACCHgkQEjzqDGoCgIIIIAAAggggAACCCBASKMPIIAAAggggAACCCCAAAIeCRDSPGoMqoIAAggggAACCCCAAAII/H92MOxsoujzsgAAAABJRU5ErkJggg==",
"text/html": [
"