{
"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"
]
},
{
"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": "iVBORw0KGgoAAAANSUhEUgAAAwUAAAFoCAYAAAAGkFfUAAAgAElEQVR4Xu2dCZxeVX2w/8ms2Sf7JCGB2GKsQEVL8VORgCwCUoyoXwAbxKUGPlmCFGgRTWNRWURAsBC1lqUIaZEvAgqyCAGslVL1K1CLtEICSSaZzGTPTGb9/udO7nDn5l3u+55737s996e/YWbuOfec53/mzXnu2UYN6iVcEIAABCAAAQhAAAIQgEBuCYxCCnIbeyoOAQhAAAIQgAAEIAABhwBSQEOAAAQgAAEIQAACEIBAzgkgBTlvAFQfAhCAAAQgAAEIQAACSAFtAAIQgAAEIAABCEAAAjkngBTkvAFQfQhAAAIQgAAEIAABCCAFtAEIQAACEIAABCAAAQjknABSkPMGQPUhAAEIQAACEIAABCCAFNAGIAABCEAAAhCAAAQgkHMCSEHOGwDVhwAEIAABCEAAAhCAAFJAG4AABCAAAQhAAAIQgEDOCSAFOW8AVB8CEIAABCAAAQhAAAJIAW0AAhCAAAQgAAEIQAACOSeAFOS8AVB9CEAAAhCAAAQgAAEIIAW0AQhAAAIQgAAEIAABCOScAFKQ8wZA9SEAAQhAAAIQgAAEIIAU0AYgAAEIQAACEIAABCCQcwJIQc4bANWHAAQgAAEIQAACEIAAUkAbgAAEIAABCEAAAhCAQM4JIAU5bwBUHwIQgAAEIAABCEAAAkgBbQACEIAABCAAAQhAAAI5J4AU5LwBUH0IQAACEIAABCAAAQggBbQBCEAAAhCAAAQgAAEI5JwAUpDzBkD1IQABCEAAAhCAAAQggBTQBiAAAQhAAAIQgAAEIJBzAkhBzhsA1YcABCAAAQhAAAIQgABSQBuAAAQgAAEIQAACEIBAzgkgBTlvAFQfAhCAAAQgAAEIQAACSAFtAAIQgAAEIAABCEAAAjkngBTkvAFQfQhAAAIQgAAEIAABCCAFtAEIQAACEIAABCAAAQjknABSkPMGQPUhAAEIQAACEIAABCCAFNAGIAABCEAAAhCAAAQgkHMCSEHOGwDVhwAEIAABCEAAAhCAAFJAG4AABCAAAQhAAAIQgEDOCSAFOW8AVB8CEIAABCAAAQhAAAJIAW0AAhCAAAQgAAEIQAACOSeAFOS8AVB9CEAAAhCAAAQgAAEIJF4KLl5+izy65vn9InXO4pPk0vPOyEUEn33uBVl62fVy5bIlcuai40bU2eWz8tpL5KgjD0scj3XrN8nJn7hcDl0wX1atXJ648lEgCEAAAhCAAAQgAAGR1EjBw3dfI/PmzHRi5naS89LRRAr4U4UABCAAAQhAAAIQiJJAKqXAALnu1nvl9lWPFHx7HiWwOPIuJQVxlIdnQgACEIAABCAAAQhki0BqpcDtKLvTiFxJMCMKl37lNnnx5VedSHmn3PinIvlHGryd7+d+/dsR05b8U3fuWf2EXHXjXfu1Bv99pcplEpfLwy2T/0EnLjxCblhx/rAceUdSzL2F0vnv8ZZtyQVfky2d24cf47/X/GLh6ReNuGfalEmy5v6byv5FHHLMOeKW11s2w6oc52KZV1J29/kf/dBCZxqWubxlD8LKLYfJy3uZNmTampfX4qUrpK2902Hjvd97j59loelwhabO+eMS5J6yAeIGCEAAAhCAAARyTyBzUmAiWmh+vdsJe+mp24eDbn5mLrdj6+0cejv3bgfUm6+RgtUPPztinrwrCoXSFipX0DxKjRR4O8fu9KpC5S31M1O2Yp1ab2fYK1HuWgFvZ7/YX1MxKTD3F2JVSEj8ebv1KVR200n3xtntmBeabhaUVbH6GgEoJAWulHrLYcpaaI1FobzdfEu11yD35P4TDgAQgAAEIAABCAQikFopcN+Quh3IQp1jl0Cx37mdMbdjWqrzbQSidfqUsotlzX3vOuxg5y2+uUqVq1iE/HlUKgX+Trj7HO8b7FJlc+XGlaByzy+34LvUSIF34bQbjyCLyMvF1JtHMR6GQVBWfnal2laxe02aYr9zmbvtuVi5TL1d3kHuCfQpwE0QgAAEIAABCOSeQCqlwO1Aed9Sl+p8l+qkmY6V24EMsqDX/+bXP53EtCjvG+kgUlAuj3KdcrO2wu1MVnJvsbL5pcDtrJu6VbPLUS2lwJTRL3DFOs+VsCrVAffyL9XxLyUhbllcvu7IVilBCnJP7j/hAAABCEAAAhCAQCACqZECf238ndNSnW///G1/Xq5cVCIF7kiFf/qMv0NaqlxB86ik8+rv0Hvr6p8qE1QKTB5eMfDm6ZekQq0uqVIQlNW8OTOcbVULddALMSwmocUYepl5p1MVarf+dh/knkCfBNwEAQhAAAIQgECuCaRGCsrNMy83UmCiXG6f/FKdb39Hr9ib40qkIGgelUhBJfdWIgX+vxK3Qx1kW9haS4H/ebYjBabuYUiBO1IQZHqUn7dXKIqJWJB7cv1pR+UhAAEIQAACEChKIBdSEGQKjyFUqkPt7ViWmvseVAoqycO/9sEbzUJ1CzpPPqgUGAGYO3vGfoejlZqW5S1jLaWgUAzDWFNQbE1JJSMFhkkQZibeqx54cr/D+bzPciS3zD3uwnM+/yAAAQhAAAIQgEA5ArmQAgOh0O5D5udmCs+R7/wj56Rg/zanLrygu7y4zwi6psC/+5G3nP438JV0Sgud4VBqlx3/KIx/Wo37vfcNdxi7D/m3b61mobF/Oo0RAD+7UlIQlFWhnaW8u1WV273JbUulTng2Mb7r5iucWwudAu1tL8XyKdSmyn0I8HsIQAACEIAABCCQGylwBeDRNc/vF3V3OobbyTN72Hv37C+2H79/PrfpGJr9/r27FJUbpQiShymwfz562OcUeN8qF5pr790C1AUYdBpMlCMF/lgV2iK1lBSYugQ9p8B/NoWRjyMOX+Acoued0hNkNKDQWgCvzBQqU7FzNbwNOsh0Lj72IAABCEAAAhCAgJ9A4qWgliHj5OBa0rZ7VjnZsss9eGoz0vSrF14JdIhb8Fy5EwIQgAAEIAABCNSWAFLg4Y0U1Lbx2Tyt1lJgRmrMKJD3BGfai00ESQsBCEAAAhCAQJIIIAVIQZLaY+CyxCEFZp6//6rm3IbAleRGCEAAAhCAAAQgUCMCSEGNQPMYCEAAAhCAAAQgAAEIJJUAUpDUyFAuCEAAAhCAAAQgAAEI1IgAUlAj0DwGAhCAAAQgAAEIQAACSSWAFCQ1MpQLAhCAAAQgAAEIQAACNSKAFNQINI+BAAQgAAEIQAACEIBAUgkgBUmNDOWCAAQgAAEIQAACEIBAjQggBTUCzWP2J1DupGGYQQACEIAABCAAAQjUhkBqpMB0IP3XoQvmy6qVyyMhVet98COpRMIzjUIKosjTxVhp3qUONzMnIT+65nlJ0jkHUR7G5v/7vXLZEjlz0XEJb6EUDwIQgAAEIJAfAomXgntWPyFX3XiXnLjwCLlhxfnDkXE7MNOmTBpxymxYoUMKwiJZPJ9KO9lBShRFnkhBEPKF7zEnQZtD37x/v1HKR/UlJSUEIAABCEAg3wQSLwWmk1dqRMC8cfXKQljhRArCIokUpOWteBSddXdE5KWnbh/REMzPf/XCK5EIffQtlydAAAIQgAAEskcgFVLgHyUoFAZ3RKFQB8zt7Jyz+CS59LwznORuZ8Wb18N3XyPz5swUVwj8z/Hm7b4B9d7jnwqyeOkKaWvvlK/+1Wdl6WXXD9/qSo5bZvcXQerp3utPa37uLZ+3g/fcr3/rTFVxLz+jQnn58zPfe0VpyQVfky2d24fzdNl5eSw8/aIR9/hHddy3+ke+84+c0aByHAzPF19+teh9haaYlRtJKhRH8wB/LCrN2+Xvb0NuvoWk0+Uxu3Wa3L7qkeGkbrv1t9lisuHn7m33xT7CypW32o8+U5bW6VP2m+bntrkkTZ+qto6kgwAEIAABCGSBQOKlwO3gBOk8mHvNteb+m0bExu2cuz93O5fet5f+tKVGCtwOjbezVUhK3Od4Rzq8ndBCPw/SgXPL5mXidhjdzrm3k+ftPBaqlyn76oefHdFxK1Qfryx5JaAQT/8IT6FpJG5HO8jUkkIjRsV+VolcueXytoVCZTUNqtKpSaXevBeTAvOcQu2q2M+9bcAtd6F2FYRJuZGCQmJU6EPQy7IYs3LPysKHK3WAAAQgAAEIpIlA4qWg2JvcQp3nQm8f3fTejnGxjorpqLkjCaWkoFh6/5QIv4y4DaPY29NiP/c3qGL5mvq/708PdUY7inW6XB5B5MOU512HHTw8PasYE3+e5TrDLuNiHP0cij230HMq7bgX+2MtNL2l0rzLcTCjAV65KpR/sXgV+nmpdmFGYgqN5njrH3ZHvZhcmWeG/aw0fehSVghAAAIQgEASCSReCrzQ/NNHzO8KTfPw/sztUPpHBczUl1Id40o6om4Z/c8q1kmr9Of+huOOCpRaa1GNFJTb4amUKHk78l6RKzXCE1QKivEyXPx5VNpxd9kWalvmd0Heehf7ww5DCgrV0X1e0Lq75Sg32hZ2Rx0pSOJHPmWCAAQgAAEIFCaQKinwV8HtyHk7O/6Oq+k4Fer8++ddm7xL5eM+u9j8e2/Z3DeylXb+S3V+i4mB9+feelYiBa5k+AUr6Bt7Uwb/vcVGeIJ0sv15lRpB8U8hqlQK3Hj61x4UWiBbad61lIJivL3to9yC53JSwPQh/hmBAAQgAAEIZJdAqqWg2OJiVwQOmDXdWcBabtqEt0PldlrLjRSUe+tqmkyUUuBvkv61F5VIQdA39tVMqfLLlHd0I+hzoxwpKJZ32qTAHVEIMi2s2MdZOSmo5mOQhcbVUCMNBCAAAQhAoPYEEi0FprO+6oEnh+f5+/EU66S688HNriezW6eO2LK0WJ7+vFzhKCQUQd8YRyUFhbZh9XfogkpBqTUGQUcK/M8y7ObOniFHHXnYiJD5eQSVgnKC5n0D7l8HUe5PqlintZAUVJp3ofUsbnlK7T7k32K3GCf/zysZaSrEpVR5y3Es9nu2JK2WHOkgAAEIQAACtSWQeCkwBx+Zy/9mvtihZuZe75t/f6e+0A4tJo1/96FSneVCuw+5zzVbdXp3OTJbkpbbDckNedBOnTsq4J2KU2z3If+UkUL1KrRrk/sM75v9QrsembL7p/AU4lNs96FCu+IU6qibZxTa0tS/rqJYJ7Rcp9XbTrxbfxZi7N9zv9SfbDHpiEIKirVtt33fdfMVziL0aspb7ccSh5dVS450EIAABCAAgdoSSLQUuCgKzf83vys1LajUPPRCe7IXWrTrXz9Q7pwCUybv9I2oRgrcjri/qXg7q0FHCooxNmyN4Hj3mHc7sqZz7j2joFDHvtBZD/6pLUFHCrzSVOqcgkJ1KXdOgUlTaP9/92wHvwB422KQvP1z/YOcU1DtSEGpv5dSi9K97ahYeW0/lvzrEcqtb7B9HukhAAEIQAACEKiMQCqkoLIqiXAwUqXEgt3PKc/BOHEXBCAAAQhAAAIQSBuBTEpBsUPM0hacpJUXKUhaRCgPBCAAAQhAAAIQCIdA5qSg2I5E4eDKdy5IQb7jT+0hAAEIQAACEMgugcxJQXZDRc0gAAEIQAACEIAABCAQDQGkIBqu5AoBCEAAAhCAAAQgAIHUEEAKUhMqCgoBCEAAAhCAAAQgAIFoCCAF0XAlVwhAAAIQgAAEIAABCKSGAFKQmlBRUAhAAAIQgAAEIAABCERDACmIhiu5QgACEIAABCAAAQhAIDUEkILUhIqCQgACEIAABCAAAQhAIBoCSEE0XMkVAhCAAAQgAAEIQAACqSGAFKQmVBQUAhCAAAQgAAEIQAAC0RBACqLhSq4QgAAEIAABCEAAAhBIDQGkIDWhoqAQgAAEIAABCEAAAhCIhgBSEA1XcoUABCAAAQhAAAIQgEBqCCAFqQkVBYUABCAAAQhAAAIQgEA0BJCCaLiSKwQgAAEIQAACEIAABFJDAClITagoKAQgAAEIQAACEIAABKIhgBREw5VcIQABCEAAAhCAAAQgkBoCSEFqQkVBIQABCEAAAhCAAAQgEA0BpCAaruQKAQhAAAIQgAAEIACB1BBAClITKgoKAQhAAAIQgAAEIACBaAggBdFwJVcIQAACEIAABCAAAQikhgBSkJpQUVAIQAACEIAABCAAAQhEQwApiIYruUIAAhCAAAQgAAEIQCA1BJCC1ISKgkIAAhCAAAQgAAEIQCAaAkhBNFzJFQIQgAAEIAABCEAAAqkhgBSkJlQUFAIQgAAEIAABCEAAAtEQQApC4LqhoyuEXMgCAoUJjB9TL6NGjZKde3pBBIHICJh2Nlrb2Q7aWWSMyVhkfLO2s9G0M9pCNARmTx0TTcY5yRUpCCHQSEEIEMmiKAGkgMZRCwJIQS0o8wykgDYQJQGkwI4uUmDHz0mNFIQAkSyQAtpArASQgljx5+bhSEFuQh1LRZECO+xIgR0/pCAEfmRRmgAjBbSQWhBACmpBmWcgBbSBKAkgBXZ0kQLl9+xzL8jSy64fJnnO4pPk0vPOGP7+4uW3yKNrnne+P3TBfFm1cvkI6owU2DVCUiMFtIH4CSAF8ccgDyVACvIQ5fjqiBTYsc+9FKxbv0lO/sTlsvLaS+SoIw8T//f3rH5CbrvzAVlz/00O6cVLV8gRhy8YIQ1IgV0jJDVSQBuInwBSEH8M8lACpCAPUY6vjkiBHfvcS4E7SvDSU7cPk1x4+kVy7tmnyZmLjttPAvySYBIhBXaNkNRIAW0gfgJIQfwxyEMJkII8RDm+OiIFduxzLwXu2/8XX35VHr77Gh0p2CxfvPp7wyMDXkEw9xaSCKTArhGSGimgDcRPACmIPwZ5KAFSkIcox1dHpMCOPVKg/K679V556LFfyJbO7Q5N75qCQ445R65ctsQZNfBKgRGIeXNmOj8bGBy0iwKpIaAEnt/TJa/09Mg68/+9Q1/X9vbJ6/q1R9vYtPp6mVJfJ29paJBFk1vk1IkTZHJdXWB2owLfyY15I7C3d0CaGkbnrdrUN4UE+Nc2hUGrYZHNWStc1RPIvRSYN//ekQGD0oiAKwZBRgo2cnhZ9S0wxym3DQzIP+zaJs907ZF/3VvdAXhHNY+RU8ZMkJPHjpeZZQSBf0xz3NgCVJ2RggCQuMWagO1IAV0+6xBkOoNZHF5mFd/cS4EZJXj+Ny+P2FHILCae3TpVblhxPmsKrJoXiQsR2DkwKLfu6JDv7twuewYHhm85tnmsvK95nEzV0YApo0bLuNFDowBjmur0RGORPd39zvcv9nTLz7p2ydPdI0XiXY3N8tmJLfLhsRMAD4GKCSAFFSMjQRUEbKWgikeSJEcEmD5kF+zcS4G7RsC/+5A7ZYjdh+waGKnfJNClAvB9FYFv7+iU7TpKYK4/qG+QCydNlQ/qm/4JRYY9i51TYPL7uYrBk127VRL2yLr+XifPeXUNckHLFDlr3ETwQyAwAaQgMCputCCAFFjAI2lZAkhBWUQlb8i9FBg6puN/1Y13DYPyriEwP+ScArtGRmrRUYFtcvP2DunYJwN/VN8oF0+eKh8aM74snqCHl/1YRw9u2Nohv+3rcfL844YmuWX6LEc8uCBQjgBSUI4Qvw+DAFIQBkXyKEYAKbBrG0iBHT8nNbsPhQAxo1k819Mll27ZJP/dN/QW33TUl7XoyMCYcYFrHFQK3Awf2bNL/qpzs7QPDE03+vLkabJ0wuTAz+PGfBJACvIZ91rXGimoNfF8PQ8psIs3UmDHDykIgV8Ws9iluwUt79wk9+7e6VTvvU3N8nmdJnSMrhuo9KpUCkz+O3RE4m9UDFbtGXr+uxvHyI3TZzpTi7ggUIgAUkC7qAUBpKAWlPP7DKTALvZIgR0/pCAEflnLYou+oT9r0xvyUu/QNJ6/bpkm50+s/k19NVLgMn1G1xxc0tEm6/v7pElGyVemTJc/Hz8pa8ipTwgEkIIQIJJFWQJIQVlE3GBBACmwgKdJkQI7fkhBCPyylIVZ7HtG23pZq1/fqucJ3DC1VQ7XXYFsLhspMM/tlkG5eusWZ12DuT6ii5pvmTbLpkikzSABpCCDQU1glZCCBAYlQ0VCCuyCiRTY8UMKQuCXlSzMyIAZITAjBe9UEbh7xhyZNNr+QChbKXD5/vveblnSvt7Z+chMJ7pDyzdhNLt+Z6X92dYDKbAlSPogBJCCIJS4p1oCSEG15IbSIQV2/JCCEPhlIQuzoHjJpvVi1hIs1APF/n76bBmjZw2EcYUlBaYsr+mC58Wb35A3+vpkvu5KtGrGATJHT0rmggBSQBuoBQGkoBaU8/sMpMAu9kiBHT+kIAR+ac/iET0n4DPtG5xqmIPD/m5aa6hVClMKTMHMrkRmROM/dWRjrgrB6plzpbUOMQg1aCnMDClIYdBSWGSkIIVBS1GRkQK7YCEFdvyQghD4pTmLH+zeIZd2bHKq8KkJk+SqyTNCr07YUmAKaEY0/nJLmzyoZxu8RUcMVrfOlan7TlAOvQJkmAoCSEEqwpT6QiIFqQ9hoiuAFNiFBymw44cUhMAvrVnctGOrXLtti1P8v9SzBy6eOCWSqkQhBaagffr/P9epRGaHooN1UfSDMw9kjUEkEUxHpkhBOuKU9lIiBWmPYLLLjxTYxQcpsOOHFITAL21ZDGqBL9PRATNKYJbpXjt1ppw1bmJk1YhKCkyB9+rORGfqWohf7u1yDla7r/UAGRfSWojIgJBxJASQgkiwkqmPAFJAk4iSAFJgRxcpsOOHFITAL21ZXKj7/v9w36Fkt02fJX82ZnykVYhSCkzB9wwOyP9WMfh1T7ccqbsS3ds6xznTgCtfBJCCfMU7rtoiBXGRz8dzkQK7OCMFdvyQghD4pSkLs9//zTu3SotuNXrbtNnyft1pKOoraikw5TcnIH900+vO4uNjmsbK3TPnRF0t8k8YAaQgYQHJaHGQgowGNiHVQgrsAoEU2PFDCkLgl5Ys7tHpQn+p04aa9S366llz5TCdblOLqxZSYOrRobsSfbjtdXlVty39QPNYuV3PMairRQV5RiIIIAWJCEPmC4EUZD7EsVYQKbDDjxTY8UMKQuCXhix+smeX/MWWjU5R79LOsuk01+qqlRSY+mzs75PTVAw26NdFur3qt0PeXrVWzHhO5QSQgsqZkaJyAkhB5cxIEZwAUhCcVaE7kQI7fkhBCPySnsXLOqXmQxvXSZcuyv2GLio+M8JFxYVY1FIKzPPX6kiBGTEw5xmYupo6c2WfAFKQ/RgnoYZIQRKikN0yIAV2sUUK7PghBSHwS3IWOwcG5cSNa2Vdv54ErG/OvxnDm/NaS4GJx+9UhBbpGoPtutZg6YTJ8uXJ05IcJsoWAgGkIASIZFGWAFJQFhE3WBBACizgaVKkwI4fUhACvyRn8eeb18uT3Xuc7TofnjUvlqLGIQWmov/Rs1c+oiMG3TpCcl3E267GApaHjiCAFNAgakEAKagF5fw+Aymwiz1SYMcPKQiBX1KzuG57h9y4vdPZaehnsw6SmXXxLLuNSwpMXB7SE4+Xtm90Fhz/UM8w+FPdspQrmwSQgmzGNWm1QgqSFpFslQcpsIsnUmDHDykIgV8Ss3hCRwc+qaMEZrf++2YeIO9uiq8zHKcUmNhcs61DvrWjUyapHD2ioyXz6hqSGDLKZEkAKbAESPJABJCCQJi4qUoCSEGV4PYlQwrs+CEFIfBLWhZmS86TdR3BzsFB+dLk6XLuhJZYixi3FJjKf6p9gzzatVve2tCgYnAgh5vF2iKieThSEA1Xch1JACmgRURJACmwo4sU2PFDCkLgl6QsuvR035Pb1skrvb1yytjx8t1ps2IvXhKkwJx6fNyGdc6C609NmCRXTZ4ROxcKEC4BpCBcnuRWmABSQMuIkgBSYEcXKbDjhxSEwC9JWXxWzyJ4WM8k+MN680Z8nowZNTr24iVBCgyEF3v3Oluz9ul//0DPalhYw7MaYg9CDgqAFOQgyAmoIlKQgCBkuAhIgV1wkQI7fkhBCPySksW3dd7813T+/DgVgcdUCA5UMUjClRQpMCxu0oXX1+oC7Km6vuCZ2fOddQZc2SCAFGQjjkmvBVKQ9Ailu3xIgV38kAI7fkhBCPySkMUv93bJxza9IQNamDv1LfhxCXoLniQpMHzMice/7umWE8aMk9unz05C+ChDCASQghAgkkVZAkhBWUTcYEEAKbCAp0mRAjt+SEEI/OLOYn1fn5zYtla26UFdF0ycIn/VMjXuIo14fpKkwBRsQ3+fHLvhNdmlC7GvnjJDloyflCheFKY6AkhBddxIVRkBpKAyXtxdGQGkoDJe/ruRAjt+SEEI/OLO4oO609CLeoLv+5vHyL0zDoi7OPs9P2lSYAr4E1138Re6/sJcz84+SOYnZKpV4oKXogIhBSkKVoqLihSkOHgpKDpSYBckpMCOH1IQAr84s/ibre3y3Z3bZHZdvTyuW20mcY58EqXAxOwLW9pk1Z6d8o7GJvlJazynPcfZdrL2bKQgaxFNZn2QgmTGJSulQgrsIokU2PFDCkLgF1cWP9MDypbsO6DshzEfUFaKQVKlwLtN6YU67eryhE27iqtdpfW5SEFaI5euciMF6YpX2kqLFNhFDCmw44cUhMAvjizaB/rl6A2vyo6BQblwknZoJyVrHYGXSVKlwJTR3abULED+0cy58q6m5jjCyTNDIIAUhACRLMoSQArKIuIGCwJIgQU8TYoU2PFDCkLgF0cWZqehX+iOQ3+indgHtDOb5CvJUmC43bxjq1y9bYszBetJXV8wftSoJOOkbEUIIAU0jVoQQApqQTm/z0AK7GKPFNjxQwpC4FfrLO7ctV3+unOz03l9fPaBMrcuGecRFOOQdCkw5T697Q35ZU+XfHTcBPnW1NZah5TnhUAAKQgBIlmUJYAUlEXEDRYEkAILeJoUKdjH757VT8hVN97lfDdtyiRZc/9Nw2QvXn6LPLrmeef7QxfMl1Url4+gvqGjyy4KpK4ZgTbdTvPoDWtl9+CAfF230zw7BdtppkEKvNuUrpw+S04dM+WCLX0AACAASURBVL5mMeVB4RBACsLhSC6lCSAFtJAoCSAFdnSRAuV33a33yu2rHpGXnrp9P5pGFm6784FhSVi8dIUccfgCufS8M4bvRQrsGmEtU5+5+Q15urtL3tnYLA+1JnvakMslDVJgyupuUzpBR2Ce1tOOZ9TV1TK0PMuSAFJgCZDkgQggBYEwcVOVBJCCKsHtS4YUKIhDjjlHHr77Gpk3Z+Z+NP0S4JcEkwApsGuEtUr9oO6tf67urd+gndYndfvRtOytnxYpMHG8uGOT/NPuHfK/msbIfbqjE6sLatW67Z+DFNgzJIfyBJCC8oy4o3oCSEH17EzK3EvBs8+9IF+8+nuypXP7MMkTFx4hN6w43/l+4ekXyblnnyZnLjrO+d7cv/Sy60eMKiAFdo2wFql36i5D79nwe9mqpxZfqrsNLUvwbkN+HmmSAu82pVfoFqWf161KudJBAClIR5zSXkqkIO0RTHb5kQK7+OReCty1BN6pQ2bk4JzFJzlThMx/X7lsyX5S4B1Z2Lqrxy4KpI6cwAWbNsrdO7bL2xob5ekD3yL1kT8xvAc0N9bpG/dR0tXTF16mEeb0H3u75fh1rzlPeOrA+fJ2PdyMK+EEBkWam9LVzhJOlOIVIdDcoO1MR2vT8nlGINNFYPL4xnQVOGGlRQp8awZMfMzCYnOZ0YIgIwVde/sTFlaK4yWwZvduOXntUCf16be8RY5oHpMqQPV1OglH/9fXpz23lFzf6NgiX960SQ5ubJBf/sHB0sw2pYmO3MDgoDTWj9bOmkhvitpZoqFSuIIEzOcZ7YzGERWBMfpyg6t6ArmXgkLTgbxSwJqC6htXElJ2y6Acu36trOvvlc+Mb5GvTJmehGJVVIY0TR/yVszdpvQvJrTI30xOH/eKgpSBm5k+lIEgpqAKTB9KQZBSXESmD9kFL/dSYPCZ0YBTT3iPM11o3fpNcvInLpeV114iRx15mLD7kF0Dizv1VVu3yK07t0qr7oTztB6sNW7U6LiLVPHz0yoFr/b1yrEb10qvvoW+Tw+Iew+nHVcc+1omQApqSTu/z0IK8hv7WtQcKbCjjBTs42fWDriXdw2B+RnnFNg1srhSv9i7V07ZuE7M5K57ZhwgR6ds2pDLLa1SYMp/y45O+fq2DpkxWqVsznwx25VyJZMAUpDMuGStVEhB1iKarPogBXbxQArs+Dmp2X0oBIghZ2FE4PiNr8nventl0dgJ8u1p6T1lN81SYOJwio4WvNjbIx8bN1Fumrr/tr8hh57sqiSAFFQJjmQVEUAKKsLFzRUSQAoqBOa7HSmw44cUhMAviixu3bFNrtrW7ryZ/rm+oZ6qb6rTeqVZCgzz36oQHK9iYK47ps+W48eMS2soMl1upCDT4U1M5ZCCxIQikwVBCuzCihTY8UMKQuAXdhZmUfF71r/mZPutqa3y0XETwn5ETfNLuxQYWF/ftkWnEm2Vmbq2Y80snUY0mmlENW1EAR6GFASAxC3WBJACa4RkUIIAUmDXPJACO35IQQj8ws7i9E1vyC/3dsm79VTd+/VU3bRfWZACs9jYLDo2i4/PGj9JrpsyI+1hyVz5kYLMhTSRFUIKEhmWzBQKKbALJVJgxw8pCIFfmFncu2u7XNK5WZp1Y/+n5hwoc+sawsw+lryyIAUG3K97uuXP2l7XTWLFkTUjbVzJIYAUJCcWWS4JUpDl6MZfN6TALgZIgR0/pCAEfmFl0T7QL+9f/6rs1LfSX26ZJksnTg4r61jzyYoUGIhfVGG7XcVtbn29PKVbxBp540oGAaQgGXHIeimQgqxHON76IQV2/JECO35IQQj8wspi6ZaN8tCeXXJoQ6P8ZNaBkt6lxSOJZEkKdg8OyNEbXpO2/n753MQWWd7CoWZhtX/bfJACW4KkD0IAKQhCiXuqJYAUVEtuKB1SYMcPKQiBXxhZPN61Wz7ZvsERgcdVCN6qYpCVK0tSYGLydHeXnLn5DWeM4JFZ81TimrISqlTXAylIdfhSU3ikIDWhSmVBkQK7sCEFdvyQghD42WZh3j6/T3cbMtOHzp84Rf66ZaptlolKnzUpMHAv6tgk9+3eofLWII/NOkjqE0U8n4VBCvIZ91rXGimoNfF8PQ8psIs3UmDHDykIgZ9tFlfoPPU79s1Tf0Y7mA0ZOzU3i1Kwc2BQ3rPh97J1YEAumTRFvjApWyJn26bjSI8UxEE9f89ECvIX81rWGCmwo40U2PFDCkLgZ5NFHna0yaIUmJg/qOs/ztV1IGaU4LGMTfmyadNxpUUK4iKfr+ciBfmKd61rixTYEUcK7PghBSHwqzaLPk14jC5aNXvfn6l7338jo3vfZ1UKTNw/uXm9PN69x1kc/oiKAXsRVfvXYJ8OKbBnSA7lCSAF5RlxR/UEkILq2ZmUSIEdP6QgBH7VZvGN7Z1yw/YOmT66Tp6Znd1TcrMsBd5tZM1ORGZHIq54CCAF8XDP21ORgrxFvLb1RQrseCMFdvyQghD4VZOFGR0wowRmtGDltFly6tjx1WSTijRZlgITgB/oepBLM3bgXCoalq+QSEEao5a+MiMF6YtZmkqMFNhFCymw44cUhMCv0izMibjmZFyznuD45rFyx4w5lWaRqvuzLgUmGKdvekN+ubfLOeXYnHbMVXsCSEHtmefxiUhBHqNeuzojBXaskQI7fkhBCPwqzeL7O7fJl7a2y7hRo+Xncw5ypg9l+cqDFLzeryM/69dKtwzK9bo25AxdI8JVWwJIQW155/VpSEFeI1+beiMFdpyRAjt+SEEI/CrJoq2/T0/EXSvmbIKva+fx7Bx0HvMgBaYNfGfHVlmxbYtM0C1ln5kzP/OyV0m7r8W9SEEtKPMMpIA2ECUBpMCOLlJgxw8pCIFfJVmYk3DNibjvbGyWh1rnVpI0tffmRQrMtLCTNq6VF3t7cjEtLGkNEilIWkSyWR6kIJtxTUqtkAK7SCAFdvyQghD4Bc1i9e6d8vmONudwsid1+8r59Q1Bk6b6vrxIgQnS71QITlAxyMMC8qQ1SqQgaRHJZnmQgmzGNSm1QgrsIoEU2PFDCkLgFyQL7wm4l+oJuMtydAJunqTAtIXrdavZb+Zgq9kg7b6W9yAFtaSd32chBfmNfS1qjhTYUUYK7PghBSHwC5LFRR2b5L7dO+StDQ16+u1Bzim4ebnyJgVmlOCEja/pqEGvfHzcRLlx6sy8hDrWeiIFseLPzcORgtyEOpaKIgV22JECO35IQQj8ymVhtqo0W1aa024fmTVPT79tKpckU7/PmxSY4L3Yu1fXF6zTvYhE7plxgBzdPCZTMU1iZZCCJEYle2VCCrIX0yTVCCmwiwZSYMcPKQiBX6kszBaV5pCy1/v65DMTWuQrk6dH/MTkZZ9HKTBR+Jtt7fLdHdukta5Ot56dL82OFnJFRQApiIos+XoJIAW0hygJIAV2dJECO35IQQj8SmWxQjuG39nXMXx69kHO2QR5u/IqBV4h/LQK4d/mUAhr2daRglrSzu+zkIL8xr4WNUcK7CgjBXb8kIIQ+BXLwkwhOUWnkPTrDXmeQpJXKTDtwjt17EHdgtZsRcsVDQGkIBqu5DqSAFJAi4iSAFJgRxcpsOOHFITAr1AWRgSO3zi02HTR2Any7WmtET0p+dnmWQpMdP6yc7Pcs2u7swWt2YrWbEnLFT4BpCB8puS4PwGkgFYRJQGkwI4uUmDHDykIgV+hLG7e0SlXb+twTrf9xZy3yOTR+Zs25HLJuxSY7Wjfv+FVaR/olwt1O9rLc7QdbUR/XgWzRQpqSTu/z0IK8hv7WtQcKbCjHLsUHHLMORXX4NAF82XVyuUVp4sqwYaOrqiyzmW+r/f3aidwrfQODjrbUZptKfN85V0KTOwf79otn2zf4GxF+5iOFry1oTHPTSKSuiMFkWAlUx8BpIAmESUBpMCObiKk4KWnbg9ci3tWPyGrH34WKQhMLH03fnzT6/Ive7udbSjNWoK8X0jBUAs4t32jPNi1S7ekbZSfqhhwhUsAKQiXJ7kVJoAU0DKiJIAU2NFFCuz4OakZKQgB4r4sfrh7p1zY0SZNuv3kmjkHyty6hvAyT2lOSMFQ4Lbo9KGj1r8qO3UE6cqW6XLexJaURjSZxUYKkhmXrJUKKchaRJNVH6TALh5IgR0/pCAEfm4WOwYG5L06d3yrfl2unb7P0elz0CAFbzayf9ZTrZfp6dZmsfEzs5HGEP/8nHY2Wrnu2NMbZrbkBYERBJACGkSUBJACO7qxS4Fd8ZORmpGCcOJwme4yc7fuMmO2nTTbT7LHzBBXpGBk+zpz8xvydHeXvLtpjNw/k+ll4fz1DbUzpCAsmuRTjABSQNuIkgBSYEc3FVLgLkauZO1BtVgWnn6RtE6fMmLNwsXLb5FH1zzvZFlokTNSUC3tN9P9u64hOE3XEpiFpE/pIWVm+0kupKBQG2jr75OjdSH67sEB+fqUGXL2+Ek0lRAIIAUhQCSLsgSQgrKIuMGCAFJgAU+TJlYKzILiq268a7h2D999jcybM9OutmVSGyEwl1cKTDluu/MBWXP/Tc7vFi9dIUccvkAuPe+M4dyQAvuwuIuLL2+ZJhdOnGyfYYZyYKRg/2DeoSNKV+jIkjnh+mmdRtRaZ3SSy4YAUmBDj7RBCSAFQUlxXzUEkIJqqL2ZJnFS4H0rb4q58tpL5KgjD7OrZYDUprO/6OSj5I2N7fL8b14eHinwS4BfEkzWSEEAwCVuuUHPJPiGnknwtvpGeUI7eFwjCSAFhVvEorbX5d96uuUDzWPlrhlzaDaWBJACS4AkD0QAKQiEiZuqJIAUVAluX7LESIH3vIIrly2RMxcdJ+ZntZACb8f/ulvvHSEFZvTg3LNPc8pjrmefe0GWXna9eKcyIQXVN8L/7O2RkzeudTJ4WLeZfDv7z+8HEyko3L5e7euVY7XtmPMsbtETrz+iJ19zVU8AKaieHSmDE0AKgrPizsoJIAWVM/OmSIQUFFszUAspMCMT5rphxfnOV78UmDK4kuKVAu90pr29/XZRyGnqPq33Ea+8Iv+1t0eunDFd/z8jpyRKV7tu9NCS63492ZdrJIFr29vly5s2y5S6OnnxrQc7X7kqJ9DXPyhNDaN1cf8o6dPdv7ggEBUBPs+iIku+hkBTA/8G2LSExEjBtCmThuftuxWqhRSYkYAtndv3Y+iWJ8hIwZbte21ikNu012ztkG9s3SKHNDbJ4wcc5Cwy5tqfwJimOhmlW0Xu6TYaxeUlYHR84euvyss64nT6+ImycsYsAFVBwOjmWG1nZveh3bSzKgiSJCiBMY3azvRFB+0sKDHuq4TAtElNldzOvT4CiZACUybv9CF3ylAtpMDfIvwjBawpiOZvxjtt6PFZB8nBDew2VIw004dKt8EXe/fKKRvXiREEcwK2OQmbq3ICTB+qnBkpKifA9KHKmZEiOAGmDwVnVejOxEiBW7i4Fhq7z/dLAbsP2TWwQqnN++7jN74mr/T2ymWTpspFk6aE/5AM5YgUlA/m13Wh+i26YL1Vpw89rVvaml2JuCojgBRUxou7qyOAFFTHjVTBCCAFwTgVuytxUuAW1LslaaGpRXbVLp7aLwXmTs4pCJf2Nds75FvbO51FxT/VxcV030rzRQrKtz+z2NgsOjaLj5foNKKrp0S7fXH5EqXvDqQgfTFLY4mRgjRGLT1lRgrsYpVYKfBWq5aHl1WDk92HglMz04Y+qJ23Op27/FjrgUwbCoAOKQgASW/5tW5PeqpuU2ouc9KxOfGYKzgBpCA4K+6sngBSUD07UpYngBSUZ1TqjtilwHT4Kzmp2IwgrH742REnDtshsE+NFARjaN7mntC21pk29FctU+WCiUwbCkIOKQhCaegec6CZOdhsbn29PKNrVRpUPrmCEUAKgnHiLjsCSIEdP1KXJoAU2LUQpMCOn5MaKQgG8Wqd932zzvtm2lAwXu5dSEFwXrsHB+ToDa9JW3+//B89GfuLekI2VzACSEEwTtxlRwApsONHaqQgyjaAFIRAFykoD5FpQ+UZFbsDKaiM3dPdXXLm5jfE7Fb9k1nz5NAGtqgLQhApCEKJe2wJIAW2BElfigAjBXbtIxFSUGkVDl0wn+lDlUKL8X7vtCHz5ta8weUKTgApCM7KvfPCLW3ywz075a261e1jOo2IMzDKM0QKyjPiDnsCSIE9Q3IoTgApsGsdsUuBXfGTkZqRgtJx+Oq2LfJ3O7bKO/SQsoda57HbUIXNFimoEJjevlVP5X3/hledr2x7G4wfUhCME3fZEUAK7PiRujQBpMCuhSAFdvyc1EhBcYj/r2ev7gizztltaI1uP3pgPYeUVdrkkIJKiQ3d/8DunXJeR5uz2PhJbXvzaXslQSIF1bUzUlVGACmojBd3V0YAKaiMl/9upMCOH1JQgp+ZNrRww1pZ298rX5o8Tc6dwLShapobUlANtaE0n9y8Xh7v3iPvbGyWB1vnCnsRFWeJFFTfzkgZnABSEJwVd1ZOACmonJk3BVJgxw8pKMHvb7dukdt2Dk0b+rFOG6JDVl1jQwqq42ZStQ/0y/vWvyZmV6K/nTxdPj2hpfrMMp4SKch4gBNSPaQgIYHIaDGQArvAIgV2/JCCIvzMtKEP6bShRlWBJ2czbcimmSEFNvRE7tZzCy7T8wuatS3+fM5B0lrHsuNCRJECu3ZG6mAEkIJgnLirOgJIQXXc3FRIgR0/pKAAv70yKMeuH5o2tLxlunxuIm9nbZoZUmBDbyjt6ZvekF/u7ZKjm8fIPTMOsM8wgzkgBRkMagKrhBQkMCgZKhJSYBdMpMCOH1JQgN+Kbe3ynR3b5E91Hvf/ZR63dQtDCqwRyusqqMeoqHarsN44daZ8fNxE+0wzlgNSkLGAJrQ6SEFCA5ORYiEFdoFECuz4IQU+ft5pQ0/PPkgOqGeqhm0TQwpsCQ6lX6nb4n5Ft8edoLsRPTNnvkwfbY4343IJIAW0hVoQQApqQTm/z0AK7GKPFNjxQwo8/LzThr6iizo/w6LOEFqXCFIQCkbp12xO2bhWXuztkVPHjpeV02aFk3FGckEKMhLIhFcDKUh4gFJePKTALoBIgR0/pMDDb/nWdvneTqYNhdCkRmSBFIRH9HcqBCeoGPRplndMny3HjxkXXuYpzwkpSHkAU1J8pCAlgUppMZECu8AhBXb8kIJ9/My0oY+0ve6cVvzknANlbh2HlIXQtJwskIKwSA7l880dnXL9tg6ZodOH1ugUt4mjTavlQgpoA7UggBTUgnJ+n4EU2MUeKbDjhxQoAbMX/Al6SJn5ep0u4jyLRZwhtKo3s0AKQsXpZHbsxtfkd729cqa21W9om+Uaks/Rut5ix55ecEAgMgJIQWRoyVgJIAV2zQApsOOHFCiBj+p2j/+q2z2erPO0v8c87RBa1MgskILQkcpverr1HI3XnYwfnDlX3tXUHP5DUpYjUpCygKW0uEhBSgOXkmIjBXaBQgrs+OVeCr6q24/+nW4/OkcPhPqZTsUYr28aucIlgBSEy9PNzW27ZhrRU9p2J+V8GhFSEE07I1ffS45mHZEazYgU7SIaAkiBHVekwI5frqXAjA6YUQJzPaBvW/+Et60htKb9s0AKIsHqZOqOch3fPFbumDEnugelIGekIAVBykARGSnIQBATXAWkwC44SIEdv9xKwf/09cpJG9fJnsEBuaJlmnx+4uQQSJJFIQJIQXTtwqyDOVHXw2zWr1fpNrqfyvE2ukhBdO2MnN8kgBTQGqIkgBTY0UUK7PjlUgp2Dg7Kh3RbRyMGHxozXr4znf3eQ2hGRbNACqKkK/JvPV2yqG1oxOu7uibmFF0bk8cLKchj1GtfZ6Sg9szz9ESkwC7aSIEdv1xKwZLN6+Vn3XvkbfWN8uPZ86RZWEcQQjNCCqKEWCbv7+/aJl/qbJcmbcsPzZonb29ojLE08TwaKYiHe96eihTkLeK1rS9SYMcbKbDjlzspuEb3d/+W7vPeoosyH5t1oMzWBcZc0RJgpCBavm7uF2zZKPfv2eW06UdUDKbqAuQ8XUhBnqIdX12RgvjY5+HJSIFdlJECO365koJHunbLZ9o3iOkq3dd6gBzZOCYEemRRjgBSUI5QOL/v02lxp216XcxBfO9obHIWz9fnaDctpCCcdkQupQkgBbSQKAkgBXZ0kQI7frmRgpd7e+QUXVjcLYPyFV2Q+ZkcL8gMoclUlAVSUBEuq5s7dMGxWUC/ob9PTte1BTfn6NwNpMCq6ZA4IAGkICAobquKAFJQFbbhREiBHb9cSMG2gQE5QRcWm47Sx8dNkBuntoZAjSyCEkAKgpIK5z5z0vHJ2t7zJsBIQTjth1wYKaANxEcAKbBjjxTY8cu8FPRrDT+mO7M8pzu0HKqLL3/cOi9XUypCaB7WWSAF1ggrzuBRnSr3qX1T5e6deYC8tyn7U+WQgoqbCQmqIMBIQRXQSBKYAFIQGFXBG5ECO36Zl4LlW9vlezu36aLL0fL4rINkRl2+Fl+G0Dyss0AKrBFWlcH12zvlm9s7ZIKuK3h09oEyr66hqnzSkggpSEuk0l1OpCDd8Ut66ZECuwghBXb8Mi0FD+zZKedtaROzv9CPWufK4Y3NIdAii0oJIAWVEgvvfjNaYEYN3lLf4OxING7U6PAyT1hOSEHCApLR4iAFGQ1sQqqFFNgFAimw45dZKXixd6/82cbXpUcXFps1BGYtAVc8BJCCeLibp5p1BWZ9gVlncEzTWPnHmXMyeyoHUhBfO8vTk5GCPEW79nVFCuyYIwV2/DIpBWYHluM3rJXN+tXsMmR2G+KKjwBSEB9782SzwP4kFYMOXXD/+YmT5YqWafEWKKKnIwURgSXbEQSQAhpElASQAju6SIHyW3j6RbKlc/swyZXXXiJHHXnY8PcXL79FHl3zvPP9oQvmy6qVy0dQ39DRZReFBKU2e7V/WPdq/43u1W7OITDnEbCKIN4AIQXx8jdPN2cXnNa2Tvr0v1dOnyWnjhkff6FCLgFSEDJQsitIACmgYURJACmwo5t7KVi3fpPc8J1/lhtWnO+QvGf1E3LVjXfJS0/dPvz9bXc+IGvuv8n5fvHSFXLE4Qvk0vPOGCafJSlY1tEm/7x7p3Oqqzmx2JxczBUvAaQgXv7u0+/Xv4sL9O+jUScQPaTrCw7R3biydCEFWYpmcuuCFCQ3NlkoGVJgF8XcS4Efn5GEkz9xuTx89zUyb87M/STASINXEkz6rEjBal1Y/HldWDzGdHpmz5O31Wer02P3pxJfaqQgPvb+Jy/v1N24dm2TmboL109VmqePzs44GlKQnHaW5ZIgBVmObvx1QwrsYoAU+Pg9+9wLsvSy64dHCszUonPPPk3OXHScc6f/9+Znm7Z220UhAal/smeXfHrzBqck35sxW07V01y5kkFgXHO9jNJtMXd19SajQDkvxWI9t2NN9x55R2OTPKDy3JSJpcejZFxznYzWdraTdpbzFh5t9cc11cvo0bSzaCnnN/eZk9kl0Sb6SIGPnpGAU094z/D0oEOOOUeuXLZkPylwRxJM8v6BQZsYxJ72pzt3yYd+/6pTjutmz5KLp2dzIWXsoKssgPbTnEuXe3AlgMD2/gF59yuvyH/v7ZGlU6fItw+Yk4BS2RWhr39QGuqHGhrtzI4lqUsT4POMFhIlgToVTq7qCSAFHnZGCN512MHD6wvMr4KMFKR5+tDLvT3yAd1ZxVznjJ8kX50yo/rWRMpICDB9KBKsVpn+vq9HPqxb9nYODsgn9e/maxn4u2H6kFWTIHFAAkwfCgiK26oiwPShqrANJ0IK9qEoJATmV/6FxVlaU/BvPV3y55vWyy59NWh2UzG7qnAljwBSkLyYmBL9VoX6I7oj0U79+/k/E1vkiy3p3roXKUhmO8taqZCCrEU0WfVBCuzigRQoPzNF6JzFJ43YUcjF6peArOw+9K97u2SJriHYo286FzaPkR/MOMCuJZE6MgJIQWRorTM2W5V+fNMbslv/jq7W0YIlOmqQ1gspSGvk0lVupCBd8UpbaZECu4jlXgrcLUj9GE9ceMTwNKKsnVPw8+4uOXvzeue01hPHjJN/mD7brhWROlICSEGkeK0zf05H3M5qWy9d+vd0SctU+cLEKdZ5xpEBUhAH9fw9EynIX8xrWWOkwI527qXADt9Q6jStKXhahWDJ5jecQ5hOUiFYqUJQHwYE8oiMAFIQGdrQMn6+p1vO0ql4ZsTgs+NbZMWU9E0lQgpCaw5kVIIAUkDziJIAUmBHFymw45cqKXhCt1H8tI4QGCH48NgJcvO0Vk4rDiH+UWeBFERNOJz8zVSixXoauFljcMa4CfKNqa2p2qwUKQinHZBLaQJIAS0kSgJIgR1dpMCOX2qk4BE9h2Dplo2OEHxMOyw3pqzDEkKYUpsFUpCe0L2ki48/rmKwfWDAWbx/qy7eT8uZ4EhBetpZmkuKFKQ5eskvO1JgFyOkwI5fKqTgR3pS8QV6UnG/lvascRPl2qkzU/UGM4QQpToLpCBd4fudisFH24a2K/2gTtH7zrRZUu9uzp7gqiAFCQ5OhoqGFGQomAmsClJgFxSkwI5f4qXgvt07ZVlHmy6BFPnUhEly1WTOIQgh5DXNAimoKe5QHvY/fb2OGLQP9MtRurvXHTPmSHPCVRwpCCX0ZFKGAFJAE4mSAFJgRxcpsOOXaCl4sGuXnNu+0SljFvZRDyFUqcwCKUhl2GStisHpOpWorb9fjmhsln+ceYBMSPCIAVKQznaWtlIjBWmLWLrKixTYxQspsOOXWCn4ux3b5Kvb2p3yLZs4WS5tmRZCTckiDgJIQRzUw3nmG319zjkG6/p75Y8bmuQHKgaTRydzlQFSEE7MyaU0AaSAFhIlAaTAji5SYMcvkVLw5c52+ftd25yyfUW3RvyMbpHIlV4CSEF6Y2dK3tbfJx9TMXhVRw7MiIFZfDy7LnkbASMF6W5naSk9UpCWSKWzloYXxwAAFjZJREFUnEiBXdyQAjt+iZKCvbpy4Esdm+Xu3Tuccn1XFzieMnZ8CDUkizgJIAVx0g/n2WZtwWIVg5d1EXKLjhR8X88HeXfTmHAyDykXpCAkkGRTkgBSQAOJkgBSYEcXKbDjlxgpMDue/EX7BvlvfRtpOh1mYaN5K8mVfgJIQfpjaGpgDjY7X3cBe7Rrt1OhW1TaP5IgaUcKstHOkl4LpCDpEUp3+ZACu/ghBXb8EiEF39epQl/SKUPmekdjk55SPEvm1jWEUDOySAIBpCAJUQivDDfv2CpXb9viZGjOMrhOtwiemIB1BkhBeDEmp+IEkAJaR5QEkAI7ukiBHb9YpWDnwKBc3NkmD+vBZOb6tK4d+FtdQ8CVLQJIQbbiaWrzTHeXjhpslC06rai1rk5PF58l7415OhFSkL12lsQaIQVJjEp2yoQU2MUSKbDjF5sU/Gpvt3NC8QZdxDhetzk0nYoT9aAkruwRQAqyF1NTo0499fj8LRtkjQrCKP3+cxMmyxUtU2M76AwpyGY7S1qtkIKkRSRb5UEK7OKJFNjxq7kUDOgTv7WjU765rcM5ofjQhkb57ozZMo/pQiFEMplZIAXJjEtYpVq5c6t8Xf+eewcH5W31jXKbLkI+uKH20/+QgrAiSj6lCCAFtI8oCSAFdnSRAjt+NZWCzXoI0lI9jOy5ni7nuZ8cP0lWTJ4uDQk+ECkEvLnPAinIfhP4T90o4LObN8haPc+gSccNrpg8TT47obZbCSMF2W9nSaghUpCEKGS3DEiBXWyRAjt+NZOCJ7r3yJW63ag5BMlMF/rmtFb5kC5S5Mo+AaQg+zE2NTS7E32xc7P88+6dToWPbh4jN06dJTN1zUEtLqSgFpR5BlJAG4iSAFJgRxcpsOMXuRR06EJEs7PQj/YMdRQOMdOFdHrBgfW1n14QAiqyqIIAUlAFtBQnWa1/65d3bJJdOp1ojh5y9tUpM+SEGqwXQgpS3GhSVHSkIEXBSmFRkQK7oCEFdvwilYL79I3hddu3yBt9fc5zLtFFiF+YOCWEEpNFmgggBWmKVjhlNSOC5+pUwf/Xs9fJ0BxC+DWVg+mjoxs1QArCiR25lCaAFNBCoiSAFNjRRQrs+EUiBY/o4UZf39ruHERmLnMI2fVTW+UPY1h8GAIesrAkgBRYAkxx8n/ctV2u1TMNOnSnInOZcw0uaZkmb43gswApSHFDSVHRkYIUBSuFRUUK7IKGFNjxC1UK/l23GV2xrV3MV3OZHYX+evJUOW3shBBKSRZpJYAUpDVy4ZTbTCO6RXcn+u7ObdItg06mH9bPBDNy+AchTiNECsKJF7mUJoAU0EKiJIAU2NFFCuz4hSIFv+/rkau2bpGf6giBuabpFIFlk6bIEt1dqJ6dhUKIULqzQArSHb+wSm/ONbhle6fcsU8OzESiReMmyBcmTZWDQpADpCCsSJFPKQJIAe0jSgJIgR1dpMCOn5UUtOsi4mtUBv5p9w7nzAGzq9C5KgNLdSvCsaNGh1AyssgCAaQgC1EMrw7mFOSbt3WKmVpkRg6MHHxMRw6W6aiizXklSEF4MSKn4gSQAlpHlASQAju6SIEdv6qk4FVdK/Dgnl1yo04J2LtvOsCnx7fIspYpMjXChYQhVJUsYiCAFMQAPQWPNOeW3KwjB9/ftW24tGeNm6ijjFNlTn19xTVACipGRoIqCCAFVUAjSWACSEFgVAVvRArs+FUkBb/QtQLf1DUD/7JvzYBJ/BHdVeRyPahoLicShxCJbGaBFGQzrmHVysjBt3Z0yD/s3D6cpTnY8EIddWzVLU2DXkhBUFLcZ0MAKbChR9pyBJCCcoRK/x4psONXVgqe6e6S1bu368jAbudwIvc6XWXggknR7CISQpXIIkEEkIIEBSPBRWnr75ObdOTgTp1W5F4XTZys6w4m6m5FjWVLjhSURcQNIRBACkKASBZFCSAFdo0DKbDjV1AKzN7iP9IzBu7XtQJm3YB7zdCpQX+mc3+X6j/U1Qzvh1BUskghAaQghUGLschGDr6pUxPv1s8f9/qTpmZnQfLhur1xy+jC65WQghiDlqNHIwU5CnYMVUUK7KAjBXb8hqXAHDb0QxWB/7trh/zPvvMFzC8n6T/AZm/xD+vbuvc1jwnhaWSRNwJIQd4iHk59zQuJf9LPI7Nb0XoVBff6UxWDE/XlhDkl+WDPeQdIQTjcyaU0AaSAFhIlAaTAji5SYMdPbt3SIX+/uWP4bAGT3TjdOejkseOcYftjm8daPoHkeSeAFOS9BdjX/1/2dskju3fJY7rtsXmB4V7zdStTIwcn6XTG41omyGjdAW3Hnjd/b/9kcoDASAJIAS0iSgJIgR1dpMCOn4z69X84OYyRUXK8/uNq9g3/gApBo37PBYEwCCAFYVAkD5fA73p7VQ52yaO6A9qverrFXek0Wac3fnD8ePlA01jnZQbbItNmoiCAFERBlTxdAkiBXVtACuz4yRfWb5A/GGiQD47hH1FLlCQvQgApoGlERWCbHoj2aPduHUXYKWu690i3np7sXmfqSOcs3b1opm5vepy+8DD/zQUBWwJIgS1B0pcigBTYtQ+kwI6fk3pDR1cIuZAFBAoTQApoGbUgYNrZk7pL2uqt253T1Td61iG4z5+rgrCgvlEOaWySP9a1CX/Y0KT/b6hF8XhGRgggBRkJZEKrgRTYBQYpCMDv4uW3yKNrnnfuPHTBfFm1cvmIVEhBAIjcUjUBpKBqdCSsgIB/ofHvenvkeZ1e9K86gvDS3r3yX309RXN7h0rCAt329O0qCm9XUThUvzebLHBBwE8AKaBNREkAKbCjixSU4XfP6ifktjsfkDX33+TcuXjpCjni8AVy6XlnDKdECuwaIalLE0AKaCG1IBBk96EXevfKb1UWXtYDGF/S//6vnp4R2y57yzld1yj8UWOjvLd5nBzRNLTz2oTRo+RQlQau/BJACvIb+1rUHCmwo4wUlOHnlwC/JJjkSIFdIyQ1UkAbiJ9AECkoVEqzLuFFPZvlpd5ulYS98rJKgzmrJeg1WUcUWlQgptTVyUTduc2co9Ci/90yqk4m1w99NaMO7s9a6kbLNL2fK50EkIJ0xi0tpUYK7CKFFJTht/D0i+Tcs0+TMxcd59z57HMvyNLLrpeXnrp9OCVSYNcISY0U0AbiJ1CtFBQruTmv5RWVg//UaUftfX3S2d8vnQN9srV/QLbr6e5GJvZ4TnmvlIDZ+tkRCP3/BJ8kNOjmb/W6A1yDbrFap/9v0PXT9fqzBv1ZvX7v/N/5vfn5KGe3uDodxRj6vd5n0u1Lb37v/rzeyWcoXze9k/++fM3qCpP+zfyH7hv6flDGaJnzfiEFeW8B0dYfKbDjixSU4XfIMefIlcuW7CcFD999jcybM9OOPqkhAAEI5JxAmxGGvn6VhaGvjjzo1236tcOVCff3+35n7uGCAAQg4Ccw+M4/BooFAaSgDLwgIwUW/EkKAQhAAAJVENiqYmD+bwRil446eK8eHYHo1bf6vbrFqvl/376v7vdDP/P8Xt/i9wzoz/Sr9x7/ffvn53uGJ/3QM9/8vb+MVVSZJBCAQBkCSIFdE0EKyvBjTYFdAyO1PQEWGtszJIfyBMKePlT+idyRRwJMH8pj1GtXZ6YP2bFGCsrwY/chuwZGansCSIE9Q3IoTwApKM+IO+wJIAX2DMmhOAGkwK51IAUB+HFOQQBI3BIZAaQgMrRk7CGAFNAcakEAKagF5fw+Aymwiz1SYMfPSc3uQyFAJIuiBJACGkctCCAFtaDMM5AC2kCUBJACO7pIgR0/pCAEfmRRmgBSQAupBQGkoBaUeQZSQBuIkgBSYEcXKbDjhxSEwI8skALaQPwEkIL4Y5CHEiAFeYhyfHVECuzYIwV2/JCCEPiRBVJAG4ifAFIQfwzyUAKkIA9Rjq+OSIEde6TAjh9SEAI/skAKaAPxE0AK4o9BHkqAFOQhyvHVESmwY48U2PEjNQQgAAEIQAACEIAABFJPAClIfQipAAQgAAEIQAACEIAABOwIIAV2/EgNAQhAAAIQgAAEIACB1BNAClIfQioAAQhAAAIQgAAEIAABOwJIQZX8yp1yXGW2JMs5gUOOOWc/AiuvvUSOOvIw5+cLT79ItnRud/77nMUnyaXnnZFzYlS/EgKmfXnbk5u21OfZuvWb5ORPXD78mELpKykD92afQKF29uxzL8jSy67fr/IvPXX78M/4fMt+26CGySaAFFQRn3tWPyG33fmArLn/Jif14qUr5IjDF9BBq4IlSUYSKNZpM3eZjpu5blhxvvO11L1whYCXgLez5e/Ul/s8M2nPPfs0OXPRceJ27LwdOUhDwCVQqp2Vazt8vtGOIBA/AaSgihj4JcD/j2oVWZIEAmU7+n4J8P8jCkIIlCLgvvH3S0Gpz7NCHTmvJEAcAn4CxdpZOSng8422BIH4CSAFVcTA/49iuQ+7Kh5BkpwS8E8fcjtw7j+0D999jcybM9Ohc92t98rzv3lZVq1cnlNaVLsSAsU6a6U+zwq98GBktBLq+bu3nBR4ibgjTny+5a+dUONkEkAKqoiL6bhduWyJM5xuLlcKvB22KrIlCQRGEDCd/ttXPSLmH85Cbcz8/qHHfjE8jQ18EChFoFhnrdTn2aoHntyvjRkpmN06dXgaG9Qh4CVQrJ35KZl2ZC7zUoPPN9oQBJJBACmoIg6MFFQBjSRVEXCH1OfNmeEs9mSkoCqMJFICjBTQDGpBIKgUeEfYGSmoRWR4BgTKE0AKyjPa7w7WFFQBjSRVEfDOs2XObVUISbSPAGsKaAq1IFCNFJhy8flWi+jwDAiUJoAUVNFCyu3WUUWWJIGAmHZlLndamllI/KsXXhmeHsTuHDQSGwLFOmvlPs/YfciGev7SFmtn5vPr4s99fHhNlGlX7zrs4OFpaHy+5a+tUOPkEUAKqowJ5xRUCY5kRQkU2sfbv/Uj+3jTgKoh4G03Jv2hC+aPWKDOOQXVUCWNn0CpduaukXLT+Nug+Tmfb7QpCMRLACmIlz9PhwAEIAABCEAAAhCAQOwEkILYQ0ABIAABCEAAAhCAAAQgEC8BpCBe/jwdAhCAAAQgAAEIQAACsRNACmIPAQWAAAQgAAEIQAACEIBAvASQgnj583QIQAACEIAABCAAAQjETgApiD0EFAACEIAABCAAAQhAAALxEkAK4uXP0yEAAQhAAAIQgAAEIBA7AaQg9hBQAAhAAAIQgAAEIAABCMRLACmIlz9PhwAEIAABCEAAAhCAQOwEkILYQ0ABIAABCEAAAhCAAAQgEC8BpCBe/jwdAhCAAAQgAAEIQAACsRNACmIPAQWAAAQgAAEIQAACEIBAvASQgnj583QIQAACEIAABCAAAQjETgApiD0EFAACEIAABCAAAQhAAALxEkAK4uXP0yEAAQhAAAIQgAAEIBA7AaQg9hBQAAhAAAIQgAAEIAABCMRLACmIlz9PhwAEIAABCEAAAhCAQOwEkILYQ0ABIAABCGSHwHW33isPPfYLWXP/TdmpFDWBAAQgkAMCSEEOgkwVIQCB+AksPP0i2dK5fb+CvPTU7cM/u3j5LfKrF15JdYcaKYi/rVECCEAAAtUQQAqqoUYaCEAAAhUSMFJw6gnvkUvPO2M45eKlK6StvTPVEuDHgBRU2DC4HQIQgEBCCCAFCQkExYAABLJNoJAU+DvQxb4/9+zT5Kob7xoG5B1dKETNyMbs1qnOrx5d87zz9dAF82XVyuXOf69bv0lO/sTlsvLaS+SoIw8bzuKQY86RK5ctkTMXHef8zP3e++xzFp8ki0871knvXt40bh2MAN2+6pGiZTb3eX9fKA9vvb2/z3ZLoXYQgAAE4iGAFMTDnadCAAI5I1BICvw/KyQFpuN84sIj5IYV5zvETIffXG4HvxBGc8+LL7+6XwffdOjNSEUlUmDydyXkntVPOHIybcqk4dEN92fuPW5n331WoTL76+kvj5uHt945ay5UFwIQgEDNCSAFNUfOAyEAgTwSKLamoNAbcneRbqGpOKYTftudD5SccuSOFLgiYXib9QrmMj+rRAq85SuUzv3Zw3dfI/PmzJRCZX72uRdk6WXXi3uPGYFw/9ttC97yMQUpj38h1BkCEIibAFIQdwR4PgQgkAsChUYK3A61+1a93HQiA8r/Zr4QvGJSsKGtwxlhiEIK3KlIhTr03ueZ8hpBKHS5U5yQglz8SVBJCEAgYQSQgoQFhOJAAALZJFBICkxNzRtyt7OeJynwjxR4o44UZPNvgFpBAALJJoAUJDs+lA4CEMgIgSRJgUFqpvAEXWjsLjwuNX2o1EiBO33IXXfgX9DsDzFSkJFGTzUgAIFUEUAKUhUuCgsBCKSVQCEpcDvL7rz9Wo0UGIamPO867ODhBcxmxMLsVFRo9yFbKSj2LO9ogWHxwx+vccqDFKS1lVNuCEAgzQSQgjRHj7JDAAKpIVBsobH3bX0tpcB96+8CNDJgdhYKSwq8B7UV2kXIvyWpKUep0YbUBJqCQgACEEgpAaQgpYGj2BCAAAQgAAEIQAACEAiLAFIQFknygQAEIAABCEAAAhCAQEoJIAUpDRzFhgAEIAABCEAAAhCAQFgEkIKwSJIPBCAAAQhAAAIQgAAEUkoAKUhp4Cg2BCAAAQhAAAIQgAAEwiKAFIRFknwgAAEIQAACEIAABCCQUgJIQUoDR7EhAAEIQAACEIAABCAQFgGkICyS5AMBCEAAAhCAAAQgAIGUEkAKUho4ig0BCEAAAhCAAAQgAIGwCCAFYZEkHwhAAAIQgAAEIAABCKSUAFKQ0sBRbAhAAAIQgAAEIAABCIRFACkIiyT5QAACEIAABCAAAQhAIKUEkIKUBo5iQwACEIAABCAAAQhAICwCSEFYJMkHAhCAAAQgAAEIQAACKSWAFKQ0cBQbAhCAAAQgAAEIQAACYRFACsIiST4QgAAEIAABCEAAAhBIKQGkIKWBo9gQgAAEIAABCEAAAhAIiwBSEBZJ8oEABCAAAQhAAAIQgEBKCSAFKQ0cxYYABCAAAQhAAAIQgEBYBJCCsEiSDwQgAAEIQAACEIAABFJKAClIaeAoNgQgAAEIQAACEIAABMIigBSERZJ8IAABCEAAAhCAAAQgkFICSEFKA0exIQABCEAAAhCAAAQgEBYBpCAskuQDAQhAAAIQgAAEIACBlBJAClIaOIoNAQhAAAIQgAAEIACBsAggBWGRJB8IQAACEIAABCAAAQiklABSkNLAUWwIQAACEIAABCAAAQiERQApCIsk+UAAAhCAAAQgAAEIQCClBJCClAaOYkMAAhCAAAQgAAEIQCAsAkhBWCTJBwIQgAAEIAABCEAAAiklgBSkNHAUGwIQgAAEIAABCEAAAmERQArCIkk+EIAABCAAAQhAAAIQSCkBpCClgaPYEIAABCAAAQhAAAIQCIsAUhAWSfKBAAQgAAEIQAACEIBASgkgBSkNHMWGAAQgAAEIQAACEIBAWASQgrBIkg8EIAABCEAAAhCAAARSSgApSGngKDYEIAABCEAAAhCAAATCIoAUhEWSfCAAAQhAAAIQgAAEIJBSAkhBSgNHsSEAAQhAAAIQgAAEIBAWgf8PfSNPbLWuUcMAAAAASUVORK5CYII=",
"text/html": [
"