{
"cells": [
{
"cell_type": "markdown",
"id": "29473a1b-493c-40a0-80b6-4873d04b9407",
"metadata": {},
"source": [
"## Diffusion of a bell-shaped initial concentration along a persistent gradient\n",
"Similar to the experiment `diffusion_along_gradient_1`, but now the gradient is NOT an initial condition; rather, it's a persistent dynamic condition (explored in experiment *gradient_1*).\n",
"\n",
"The one-chemical system starts out with a uniform concentration. \n",
"The persistent concentration gradient is attained by continuosly injecting and draining, at opposite ends.\n",
"\n",
"After a stable gradient is established, a one-time injection is performed, to\n",
"add a bell-shape concentration near one end of the system, on the \"uphill\" side of the gradient.\n",
"\n",
"Just as seen in the case of `diffusion_along_gradient_1`, the concentration peak\n",
"remains in place, and simply spreads out from there"
]
},
{
"cell_type": "markdown",
"id": "ac60b210-df99-4660-a428-7405ef79f553",
"metadata": {},
"source": [
"### TAGS : \"diffusion 1D\""
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "eae00b7a-49a4-4c4d-91c8-527c471d24f0",
"metadata": {},
"outputs": [],
"source": [
"LAST_REVISED = \"Aug. 8, 2025\"\n",
"LIFE123_VERSION = \"1.0.0rc6\" # Library version this experiment is based on"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "d73e2c5a-503f-4a6a-823d-e2f90e7efe5d",
"metadata": {},
"outputs": [],
"source": [
"#import set_path # Using MyBinder? Uncomment this before running the next cell!"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "9f96dbc0",
"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": "77d9a293-b76b-41a6-b056-f03374ece489",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"OK\n"
]
}
],
"source": [
"check_version(LIFE123_VERSION)"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "dc5fef25-fc4b-4497-81e0-7276922b5c16",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": 5,
"id": "99ab76d6",
"metadata": {},
"outputs": [],
"source": [
"# Initialize the system with a uniform concentration (of the only species)\n",
"chem_data = ChemData(names=\"A\", diffusion_rates=3.)\n",
"bio = BioSim1D(n_bins=200, chem_data=chem_data)\n",
"\n",
"bio.set_uniform_concentration(chem_label=\"A\", conc=150.)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "d1c3fcac-410b-4f3f-85dd-972964b7e607",
"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": [
149,
151
],
"title": {
"text": "[A]"
},
"type": "linear"
}
}
},
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA3EAAAFoCAYAAAD5MnI2AAAgAElEQVR4Xu3df6gl5Z3n8ad/2bZj1lk1qygr+EeQRV0ywRWGFXQIExSGIIGhDcHQf01HcElENLAbaHqRZVVEwwqmw/4hEdEmrEiYxZBBVsFhwJEkrMrQ5A9BRkc3aTMxP7S1W/fWma7bdevWqafO8z3n3q5zXuefxHOfp36861PV9T7fp57a8enaJ/kggAACCCCAAAIIIIAAAgiMgsAOEjeK42QjEUAAAQQQQAABBBBAAIEJARInCAgggAACCCCAAAIIIIDAiAiQuBEdLJuKAAIIIIAAAggggAACCJA4GUAAAQQQQAABBBBAAAEERkSAxI3oYNlUBBBAAAEEEEAAAQQQQIDEyQACCCCAAAIIIIAAAgggMCICJG5EB8umIoAAAggggAACCCCAAAIkTgYQQAABBBBAAAEEEEAAgRERIHEjOlg2FQEEEEAAAQQQQAABBBAgcTKAAAIIIIAAAggggAACCIyIAIkb0cGyqQgggAACCCCAAAIIIIAAiZMBBBBAAAEEEEAAAQQQQGBEBEjciA6WTUUAAQQQQAABBBBAAAEESJwMIIAAAggggAACCCCAAAIjIkDiRnSwbCoCCCCAAAIIIIAAAgggQOJkAAEEEEAAAQQQQAABBBAYEQESN6KDZVMRQAABBBBAAAEEEEAAARInAwgggAACCCCAAAIIIIDAiAiQuBEdLJuKAAIIIIAAAggggAACCJA4GUAAAQQQQAABBBBAAAEERkSAxI3oYNlUBBBAAAEEEEAAAQQQQIDEyQACCCCAAAIIIIAAAgggMCICJG5EB8umIoAAAggggAACCCCAAAIkTgYQQAABBBBAAAEEEEAAgRERIHEjOlg2FQEEEEAAAQQQQAABBBAgcTKAAAIIIIAAAggggAACCIyIAIkb0cGyqQgggAACCCCAAAIIIIAAiZMBBBBAAAEEEEAAAQQQQGBEBEjciA6WTUUAAQQQQAABBBBAAAEESJwMIIAAAggggAACCCCAAAIjIkDiRnSwbCoCCCCAAAIIIIAAAgggQOJkAAEEEEAAAQQQQAABBBAYEQESN6KDZVMRQAABBBBAAAEEEEAAARInAwgggAACCCCAAAIIIIDAiAiQuBEdLJuKAAIIIIAAAggggAACCJA4GUAAAQQQQAABBBBAAAEERkRg5SXupZdfTQfvfWj9kH3nW7enr976xTTt+xEdW5uKAAIIIIAAAggggAACS0hglBL31LPPp/seeaLzcBzYf3O6547bNv2t7tP8+5tvvZtu+dq3Uy1udadp38/z+NfruOaqK9PRI4fmuWjLQgABBBBAAAEEEEAAgSUmMGqJmyZf1fE68sDd6Ybrr10/dF0S9+BjT6fHj/44vf7C4xsO8bTv55kDEjdPmpaFAAIIIIAAAggggMDqEFgqiasP2/6Dh9Nrx95Izz15f7ri8kumHs27Dj2afvrqL9KLz3x3Q5tp369OLOwpAggggAACCCCAAAIInK0EllLi6irXl268Lj18+M4J+/oZt7p6V4te88BUQxurTyWA7e+rIY9Vn3d++V6n9P3kxVc2VPTqbWgup7k91fdX33Qgtb+rv2/2a1ccm/vy8s/+IVXrrj/ttmdr8GwXAggggAACCCCAAAIIlBFYSomrUNSSVg+VbEtc1aaquLXlq+/7oRJXr6v5/F09nLMpWW2J6+rX9121rc3l1cNAcxXIsqjohQACCCCAAAIIIIAAAmcDgaWVuFrQaqHZSombJlNVde5v//61yeyX1actcdMksb28rn2plldX/6ZN7nI2BM42IIAAAggggAACCCCAQIwAiWsNg6xwTqvQDa3ENWfPbE+a0jxcbYmbNryyLW0kLhZ6vRFAAAEEEEAAAQQQGDOBpZW47RxOWQWi6zUI7dcJNKWtr4rWHlJJ4sZ8ytl2BBBAAAEEEEAAAQRiBJZS4oZMbNJXcYtW4roOSb3M5lBHlbhYePVGAAEEEEAAAQQQQGAVCSylxHW9YmArn4mrhK2eFbMZqpy0eSZuFU9B+4wAAggggAACCCCAwGwElkrimtP6t2donIfE1UMkmy8SrytsFfb6+bf6u2a7WWan7Ho1QrOCZzjlbCHXGgEEEEAAAQQQQACBZSIwaonrOhDTZmach8RV66tniqzXXU3xX7+rrTmJSdd76JpSV/Wfx3vi6pkuq+WZnXKZTk37ggACCCCAAAIIIIBAN4FRSpyDiQACCCCAAAIIIIAAAgisKgESt6pH3n4jgAACCCCAAAIIIIDAKAmQuFEeNhuNAAIIIIAAAggggAACq0qAxK3qkbffCCCAAAIIIIAAAgggMEoCJG6Uh81GI4AAAggggAACCCCAwKoSIHGreuTtNwIIIIAAAggggAACCIySAIkb5WGz0QgggAACCCCAAAIIILCqBEjcqh55+40AAggggAACCCCAAAKjJEDiRnnYbDQCCCCAAAIIIIAAAgisKgESt6pH/izd7wcfezo9fvTH6bkn709XXH7JWbqVNgsBBBBAAAEEEEAAge0jMGqJ23/wcHrt2Bub6L3+wuMLIfrSy6+mg/c+lL7zrdvTV2/94kLWseoLXYTELWKZ9XEqWfbVNx1IX7rxuvTw4Ts3HO6nnn0+3ffIE+nA/pvTPXfcdtZEYdr2RjewZlcv5+ILL0gvPvPd6GL1RwABBBBAAAEElp7AKCXuzbfeTbd87dup66avuuGsPouo5JC4xZ8PJVKU26pFLJPE5aj3//2uQ4+mn7z4yobztPpR5p1fvkfkYmj1RgABBBBAAIEVIDBKiasrcNMqbtVN+/4v/9nch+ORuMWfEYsQrkUscxESt3i6ZWuYdyWu/hGmXXGsv1fpLjtOeiGAAAIIIIDA6hAYrcQN/cX+xq98c3I0u4ZptWWwvolsHv562FstcO1otIfFVev71Xu/WW/WvlFtiuCzz720YTjokQfuTjdcf21qL2OWqmK7b7taWd+QX/8n/24ydK/+tPeji0XVtmsY4NBlVv3rCkyTY3P/msJ1+3/6bxtYdnHoOi5dy2sft5wodG1ntYzosutKcXN76mPU9SNBk0dVfW5+qh8x2sdp2pDEeqhmu3/uUte3vbm+0/5eb0ud92a7Kr+XfvbCdPTIodLF64cAAggggAACCCw9gVFKXH2DPeTZofqGsX3T3r5hrv+7ucx2375KXH0zfc1VV67fgNbfNcWnKR3Nm9imNHR9P+Q5v+qGO7f++qa8a5uajOptb663a3+qM6RvmU2eXRXUtmQ3n5NqClPXULu6bZNX33ezyHB1PC679OINz6Z1DQEsqfJNq2z1SVzFuXksmrLe/r4tQfV2d+VqCJO+Sty051LbV86uY9S1bkMql/7fHDuIAAIIIIAAAnMgMEqJq/a7XXGqvmsKTO7X/erG9qev/mK9QjftZrwSl7/9+9cmE5n0Sdy0m89aBOsb1mnLmFad6KtaNPexb9uqfasnyph2Qz60AtLmVktcV4WulodaMqatu7l9045DF4dpy2sfixLR6jq3uob7lSy7ROLawtNmW29v+/u+XFTH/AvXfm7TBCvtfZ/3cMouGa7XSeLmcFW3CAQQQAABBBBYegKjlbj6yLRnuKu+bw8p67rRrm5Mu6puVf9pVa++G+LcjXldiZi2jFm/byezOayua5ha3X5WiRsyA2hutsX2MNG+CupQictJa/M1BSWiVfGaNoS2uf0ly85lpVkRnbb8od/3bV91bKtPbugiiVv6fwfsIAIIIIAAAgiMjMDoJa7Nu75pbYtC80a0btOWta7nhprVvWniMO35sea21Tfms8raLJOpTNuO5n4OlbiaRVuIuypAQyWu4tFVQR0y1K5dieurULaHVJaIVr2dXUNMxyRx057tq7M5ZFp/wylHdlW3uQgggAACCCCw9ASWTuKmPbfVFLehw8jaz97lKnFDntFbpMS101qLTlNEh0rctGFts0hcTp6a0lmL5tlQicvJ+pgkbtoPFrNc2eZdiTOxySz0tUUAAQQQQAABBDYTGKXEVSLRflFyvWu5G/Dq2a32+6mqvtOW2byB7ZsCfeizPIuSuOrG+N9e9m8ms1s2P+3tGipx056Rm0Ximuuu2B390f/Z9BLrtrQNlbhqH4c+E9d+LjF3IZgmGV1T48+67Loa2fUsWm52yisuv2R904cOp5ylkjuNy9AfPXJc6797xcBQUtohgAACCCCAAALdBEYpcfVQt2nvmZo2RKwWkL5JOJpD+7pmtpwmN12zU9bIqz5P/I//PHlv3SIlrnplQJNJV1VyqMR1TT7RHJrXHqLZnlSmPaRxGp9ps1O2J/LoEqt6HV3PkDWP4zRpmHZR6OLWrBp2MR5Sha3XN21SkkVIXP0DRfXDRftZyYrpyz/7h+zEJtO2N3JR9bLvCD19EUAAAQQQQGDVCYxS4po3pu0D2Pf+r9xMj12TeLRvfNvPneXeE1dt35Dn6maVu67gdk3y0vdsYHMZXXLafp6qYlvd9FdC0Ja4Spyb78erlt0Wsa6JQqbJ3xCJq9aRe09cvY/t5x1z74nrev9aJeLVu9raTGdddrVNzWcDh74nrqQSN23/6+/7JsFp56M+vkOeoxtyYW3ndV7LHbJubRBAAAEEEEAAgTETGK3ElUAfOo1+ybJXuc+8n5laZZb2HQEEEEAAAQQQQACBHIGVkbhpL/3OAfL3PAESl2ekBQIIIIAAAggggAAC8yKwMhLXfvZqXgAtZ/oEI9gggAACCCCAAAIIIIDA/AmsjMTNH50lIoAAAggggAACCCCAAAJbT4DEbT1za0QAAQQQQAABBBBAAAEEigmQuGJ0OiKAAAIIIIAAAggggAACW0+AxG09c2tEAAEEEEAAAQQQQAABBIoJkLhidDoigAACCCCAAAIIIIAAAltPgMRtPXNrRAABBBBAAAEEEEAAAQSKCZC4YnQ6IoAAAggggAACCCCAAAJbT4DEbT1za0QAAQQQQAABBBBAAAEEigmQuGJ0OiKAAAIIIIAAAggggAACW0+AxG09c2tEAAEEEEAAAQQQQAABBIoJkLhidDoigAACCCCAAAIIIIAAAltPgMRtPXNrRAABBBBAAAEEEEAAAQSKCZC4YnQ6IoAAAggggAACCCCAAAJbT4DEbT1za0QAAQQQQAABBBBAAAEEigmQuGJ0OiKAAAIIIIAAAggggAACW0+AxG09c2tEAAEEEEAAAQQQQAABBIoJkLhidDoigAACCCCAAAIIIIAAAltPgMRtPXNrRAABBBBAAAEEEEAAAQSKCZC4YnQ6IoAAAggggAACCCCAAAJbT4DEbT1za0QAAQQQQAABBBBAAAEEigmslMQ9+NjT6ZWfH0tHjxzaAGz/wcPptWNvbPjuwP6b0z133Lb+3Usvv5oO3vtQev2Fx4th64gAAggggAACCCCAAAIIRAmshMQ99ezz6b5HnpiwuuaqKzsl7rrPX7VB2mqwb771brrla99e50ziopHTHwEEEEAAAQQQQAABBCIEVkLiakB3HXo0vf3O8Zkkru5biyCJi8RNXwQQQAABBBBAAAEEEIgSIHFrBNvDKdtDKSvIJC4aNf0RQAABBBBAAAEEEEBgHgRIXIti/ezbkQfuTjdcf+36X0ncPOJmGQgggAACCCCAAAIIIBAlQOI6CFaVufYzcn0S9/bxD6LHYSn771jbq0sv2pf+CZ+px/eCP9qTPj71afrDhyeXMgPRndq1c0e6+IK96d1ffxhd1NL2/9efOSd9eOJU+uCjU0u7j5Ed27NrR/rj889Jv/zNichilrrvRf9qb/rdBx+nEx9/stT7Wbpze/fsTOfv25OOvy9D0xh+9o/PTf/82xOTf898NhPYd86udO7eXenXv/0IngaBy9buEX3KCZA4EleenkxPEpdHS+L6GZG4fIZIXD8jEpfPEInrZ0Ti8hkicf2MSFw3HxKXP7f6Wqy8xFWzTz78/R+mhw/fOeFUV9yee/L+dMXll6yzU4mbPWgkLs+MxJG4fEr6W5A4EhfNEIkjcdEMkTgSV5IhEldC7UyflZC45isG6l1vPvN29U0HNlBs/q39ioGqYXviE8Mpu0NI4vInJ4kjcfmUkLgII5W4PD0SR+LyKelvQeJIXEmGSFwJtRWTuBiifG8SR+LyKeluQeJIXGl26n4qcf0ESVw+YSSOxOVTQuIijAyn7KZH4iKpSmklKnExRPneJI7E5VNC4koYeSYuT43Ekbh8SvpbkDgSF82QSpxKXEmGSFwJNZW4GLVWbxJH4koDpRKnEleaHZW4YeRU4vKcSByJy6dEJS7CSCVOJS6Sn2l9VeLmQJXEkbjSGJE4EleaHRI3jByJy3MicSQunxISF2FE4khcJD8kbhH0Ti+TxJG40niROBJXmh0SN4wcictzInEkLp8SEhdhROJIXCQ/JG4R9EhcL1WzU+ZDR+JIXD4l/S08E9fPh8TlE0biSFw+JSQuwojEkbhIfkjcIuiROBIXzBWJI3HBCCUSR+KiGSJxJC6aIROb9BMkcSQueo519fdM3ByoGk7ZDVElLh8uEkfi8ilRiYswUonL0yNxJC6fEpW4CCMSR+Ii+VGJWwQ9lTiVuGCuSByJC0ZIJS4DkMTlE0biSFw+JSQuwojEkbhIfkjcIuiROBIXzBWJI3HBCJE4EheNUCJxJC4aIsMp+wmSOBIXPce6+htOOQeqhlN2QzScMh8uEkfi8inpb+GZuH4+KnH5hJE4EpdPiUpchBGJI3GR/KjELYKeSpxKXDBXJI7EBSOkEqcSF42QSlyG4N49O9P5+/ak4++fCLNe1gWoxKnElWT7sov2lXTT5zQBlbg5REElTiWuNEYkjsSVZqfupxKnEhfNkEqcSlw0QySOxJVkiMSVUDvTh8TF+E16kzgSVxojEkfiSrND4oaRM5wyz4nEkbh8SvpbkDgSV5IhEldCjcTFqLV6kzgSVxooEkfiSrND4oaRI3F5TiSOxOVTQuIijDwT102PxEVSlZJKXIzfpDeJI3GlMSJxJK40OyRuGDkSl+dE4khcPiUkLsKIxJG4SH6m9SVxc6BK4khcaYxIHIkrzQ6JG0aOxOU5kTgSl08JiYswInEkLpIfErcIeqeXSeJIXGm8SByJK80OiRtGjsTlOZE4EpdPCYmLMCJxJC6SHxK3CHokrpeq98TlQ0fiSFw+Jf0tzE7Zz4fE5RNG4khcPiUkLsKIxJG4SH5I3CLokTgSF8wViSNxwQh5T1wGIInLJ4zEkbh8SkhchBGJI3GR/JC4RdAjcSQumCsSR+KCESJxJC4aIS/7zhD0su98xLxioJ8RiSNx+bNo9hYmNpmd2aYenonrhmg4ZT5cJI7E5VPS38Jwyn4+KnH5hKnEqcTlU6ISF2FE4khcJD8qcYugpxKnEhfMFYkjccEIqcSpxEUjpBKnEhfOkEqcSlxJiLwnroTamT4qcTF+k94qcSpxpTEicSSuNDt1P5U4lbhohlTiVOKiGSJxJK4kQySuhBqJi1Fr9SZxJK40UCSOxJVmh8QNI2c4ZZ4TiSNx+ZT0tyBxJK4kQySuhBqJi1EjcYP4eSYuj4nEkbh8SvpbqMSpxEUzROJIXDRDJI7ElWSIxJVQI3ExaiRuED8Sl8dE4khcPiUkLsJIJS5Pj8SRuHxKVOIijExs0k2PxEVSlZJn4mL8Jr0Np+yGSOLy4SJxJC6fEhIXYUTi8vRIHInLp4TERRiROBIXyc+0viRuDlRJHIkrjRGJI3Gl2an7GU7ZT5DE5RNG4khcPiUkLsKIxJG4SH5I3CLonV4miSNxpfEicSSuNDskbhg5EpfnROJIXD4lJC7CiMSRuEh+SNwi6JG4XqqGU+ZDR+JIXD4l/S1U4lTiohkicSQumiETm/QTJHEkLnqOdfU3nHIOVFXiVOJKY0TiSFxpdlTihpFTictzInEkLp8SlbgIIxJH4iL5UYlbBD2VOJW4YK5IHIkLRiipxKnERTNE4khcNEMqcSpxJRkyO2UJtTN9VOJi/Ca9VeJU4kpjROJIXGl2VOKGkVOJy3MicSQunxKVuAgjlTiVuEh+VOIWQU8lTiUumCsSR+KCEVKJywAkcfmEkTgSl08JiYswInEkLpIfErcIeiSOxAVzReJIXDBCJI7ERSOUSByJi4bIcMp+giSOxEXPsa7+ox5O+eBjT6dXfn4sHT1yaMO+7T94OL127I0N3x3Yf3O6547bJt/ddejR9JMXX5n8/2uuunJT/2bHl15+NR2896FN7F5/4fH17wyn7I6m2SnzpyyJI3H5lPS38ExcPx+VuHzCSByJy6ekvwWJI3ElGfJMXAm1M31GKXFPPft8uu+RJ6ZKWCVx133+qnVpayKq+n7vBz9KLz7z3cnXfW2rv9cS15S2NnISR+JKT0MSR+JKs1P3I3EkLpohEkfiohkicSSuJEMkroTayCWu3vyqovb2O8c7K3HTJK4tbW2pa+MkceUBU4nLsyNxJC6fkv4WJI7ERTNE4khcNEMkjsSVZIjElVBbAYlrDqdsDqW88SvfTN/4+pfTV2/94oRCTtK6hlO2q3IqcSpxpachiSNxpdlRiRtGznDKPCcSR+LyKelvQeJIXEmGSFwJtSWXuCaSWsKOPHB3uuH6a9PVNx1I3/nW7Zsk7rkn709XXH5JlmZVyas+zefwfvuHj7P9VrXB+eftSb/DZ+rh33vOrvTJJ5+mj09+sqoR6d3vHTt2pPPO3ZV+/8FJfKYQOHfvrnTy5Kfp5CkZaiOqzq3du3emc9fOsz98KEPTTqJ9e3enj06eSqdOfeo86yCwa9eOdM7uXemDEzI0LSDnnbs7ffjRqcm/Zz6bCezetXPtWrQjfXjiFDwNAp9Zu0f0KScwymfi6t2dNpyyjaM5hHLWSlx7WV2Vu9+6wexMYDWc8vx9uxM+00/Qc/fsTNV9E4nrZrRzLUTnrd1g/s4N+NQQVbOeVQL3sRvwTYyqG8qqEldJ3O/dPE3N0HlrPwR89PEn6aQb8E5Gu9cuROesXav/IENTM/RHlcStSa7LUDei6jpUidwHa6Lrc4bAZ9buEX3KCaycxM36TNwQiTOcsjuAnonLn5iGU/Yz2rV283TxBXvTu7/+MA9zRVt4Jq7/wBtOmT8xDKfsZ7R3TeDO37cnHX//RB7mirYwnLL/wHvFQDcfwyljF4ylk7g333o3Pfz9H6aHD985IVPPZFkPl8zNTlm9tuDxoz9O9XNvVbXvrr/6y/WhllUl7wvXfm59+dU6SByJKz0NSRyJK81O3Y/EkbhohkgciYtmiMSRuJIMkbgSamf6jFLimq8YqHelfuat+u/qubfmp/m36vu+98S1Ja7+73p5Xe+VI3EkrvQ0JHEkrjQ7JG4YOZW4PCcSR+LyKelvQeJIXEmGSFwJtZFLXGyX59+bxJG40lSROBJXmh0SN4wcictzInEkLp8SEhdhZDhlNz0SF0lVSqOsxMV2ef69SRyJK00ViSNxpdkhccPIkbg8JxJH4vIpIXERRiSOxEXyM60viZsDVRJH4kpjROJIXGl2SNwwciQuz4nEkbh8SkhchBGJI3GR/JC4RdA7vUwSR+JK40XiSFxpdkjcMHIkLs+JxJG4fEpIXIQRiSNxkfyQuEXQI3G9VL1iIB86Ekfi8inpb2F2yn4+JC6fMBJH4vIpIXERRiSOxEXyQ+IWQY/EkbhgrkgciQtGKJE4EhfNEIkjcdEMmZ2ynyCJI3HRc6yrv2fi5kDVcMpuiCpx+XCROBKXT4lKXISRSlyeHokjcfmUqMRFGJE4EhfJj0rcIuipxKnEBXNF4khcMEIqcRmAJC6fMBJH4vIpIXERRiSOxEXyQ+IWQY/EkbhgrkgciQtGiMSRuGiEEokjcdEQGU7ZT5DEkbjoOdbV33DKOVA1nLIbouGU+XCROBKXT0l/C8/E9fNRicsnjMSRuHxKVOIijEgciYvkRyVuEfRU4lTigrkicSQuGCGVOJW4aIRU4jIE9+7Zmc7ftycdf/9EmPWyLkAlTiWuJNuXXbSvpJs+pwmoxM0hCipxKnGlMSJxJK40O3U/lTiVuGiGVOJU4qIZInEkriRDJK6E2pk+JC7Gb9KbxJG40hiROBJXmh0SN4yc4ZR5TiSOxOVT0t+CxJG4kgyRuBJqJC5GrdWbxJG40kCROBJXmh0SN4wcictzInEkLp8SEhdh5Jm4bnokLpKqlFTiYvwmvUkciSuNEYkjcaXZIXHDyJG4PCcSR+LyKSFxEUYkjsRF8jOtL4mbA1USR+JKY0TiSFxpdkjcMHIkLs+JxJG4fEpIXIQRiSNxkfyQuEXQO71MEkfiSuNF4khcaXZI3DByJC7PicSRuHxKSFyEEYkjcZH8kLhF0CNxvVS9Jy4fOhJH4vIp6W9hdsp+PiQunzASR+LyKSFxEUYkjsRF8kPiFkGPxJG4YK5IHIkLRsh74jIASVw+YSSOxOVTQuIijEgciYvkh8Qtgh6JI3HBXJE4EheMEIkjcdEIedl3hqCXfecj5hUD/YxIHInLn0WztzCxyezMNvXwTFw3RMMp8+EicSQun5L+FoZT9vNRicsnTCVOJS6fEpW4CCMSR+Ii+VGJWwQ9lTiVuGCuSByJC0ZIJU4lLhohlTiVuHCGVOJU4kpC5D1xJdTO9FGJi/Gb9FaJU4krjRGJI3Gl2an7qcSpxEUzpBKnEhfNEIkjcSUZInEl1EhcjFqrN4kjcaWBInEkrjQ7JG4YOcMp85xIHInLp6S/BYkjcSUZInEl1EhcjBqJG8TPM3F5TCSOxOVT0t9CJU4lLpohEkfiohkicSSuJEMkroQaiYtRI3GD+JG4PCYSR+LyKSFxEUYqcXl6JI7E5VOiEhdhZGKTbnokLpKqlDwTF+M36W04ZTdEEpcPF4kjcfmUkLgIIxKXp0fiSFw+JSQuwojEkbhIfqb1JXFzoEriSFxpjEgciSvNTt3PcMp+giQunzASR+LyKSFxEUYkjmtkdqAAACAASURBVMRF8kPiFkHv9DJJHIkrjReJI3Gl2SFxw8iRuDwnEkfi8ikhcRFGJI7ERfJD4hZBj8T1UjWcMh86Ekfi8inpb6ESpxIXzRCJI3HRDJnYpJ8giSNx0XOsq7/hlHOgqhKnElcaIxJH4kqzoxI3jJxKXJ4TiSNx+ZSoxEUYkTgSF8nPXCtxV990YOZtueaqK9PRI4dm7jeGDiSOxJXmlMSRuNLskLhh5EhcnhOJI3H5lJC4CCMSR+Ii+Zm7xL3+wuODt+epZ59Pzz73EokbTGw5GhpOmT+OJI7E5VPS38Jwyn4+JC6fMBJH4vIpIXERRiSOxEXyQ+IWQe/0MlXiuuGSuHzoSByJy6eExEUYkbg8PRJH4vIpIXERRiSOxEXyQ+IWQY/E9VIlcfnQkTgSl08JiYswInF5eiSOxOVTQuIijEgciYvkZ64St4gNGfMyVeJU4krzS+JIXGl26n6GU/YTJHH5hJE4EpdPCYmLMCJxJC6Sny2XuHryk1menVvEDjaX+eBjT6dXfn5s07N5+w8eTq8de2PD6g/svzndc8dtk+/uOvRo+smLr0z+f9cELSSOxJVml8SRuNLskLhh5EhcnhOJI3H5lJC4CCMSR+Ii+dkSiasmMLnvkSfW1/Xck/enKy6/ZBHbPdMym9vVJWGVxF33+avWpa258Krv937wo/TiM9+dfN3VlsSRuJkC2WhM4khcaXZI3DByJC7PicSRuHxKSFyEEYkjcZH8LFTimpWqakVHHrg73XD9tYvY3tAyq+18+53jnZW4aRLXlra21FUbROJIXGkwSRyJK80OiRtGjsTlOZE4EpdPCYmLMCJxJC6Sn4VIXPN9cd/51u3pq7d+MVXfjVHimsMpm0Mpb/zKN9M3vv7lyb5Vn5defjUdvPeh1BwmSuJIXOnJSeJIXGl2SNwwciQuz4nEkbh8SkhchBGJI3GR/Mxd4qY98zZGiWvCqSWtFtFqf2pBbUpcPVR0x8/+7yKOi2UigAACCCCAAAIIILC0BD79k3+/tPu2FTu249O1T8mKKrm5+MIL1p8Vq5cxdomr9qM5hDJXiSNxJenRBwEEEEAAAQQQQGCVCZC42NEvlrhqtc3hlM3K1diGU7YRNiXOM3HlAfOeuDw7wyn7Ge3auSNdfMHe9O6vP8zDXNEWXjHQf+ANp8yfGIZT9jPau2dnOn/fnnT8/RN5mCva4rN/fG7659+eSB+fKqoLLD01wym7D/FlF+1b+mO/yB0MSVy9YWOe2OTNt95ND3//h+nhw3dOdqeeybIeLml2yvL4kbg8OxJH4vIp6W9B4khcNEMkjsRFM0Ti+gmSOBIXPce6+s9F4uoFN6fy7xpquYgdGLLM9qsPqj7NamGzotj+W/Xf3hM3hPLmNiQuz43Ekbh8SkhchJFKXJ4eiSNx+ZT0tyBxJK4kQypxJdTO9JmrxDU35Wx82XcM1fTeZqfsZkPi8okjcSQunxISF2FE4vL0SByJy6eExEUYqcR10yNxkVSlVCRxlaA1p9jPbUJVCXv2uZc2vZ8t128sfydxJK40qySOxJVmp+5nOGU/QRKXTxiJI3H5lJC4CCMSR+Ii+ZnWl8TNgSqJI3GlMSJxJK40OyRuGDkSl+dE4khcPiUkLsKIxJG4SH5I3CLonV4miSNxpfEicSSuNDskbhg5EpfnROJIXD4lJC7CiMSRuEh+5i5xs27MNVddaTjlrNBG3t4zcfkDSOJIXD4l/S0Mp+znQ+LyCSNxJC6fEhIXYUTiSFwkP3OVuEVsyJiXqRKnEleaXxJH4kqzoxI3jByJy3MicSQunxISF2FE4khcJD8kbhH0Ti+TxJG40niROBJXmh0SN4wcictzInEkLp8SEhdhROJIXCQ/JG4R9EhcL1XDKfOhI3EkLp+S/haGU/bzIXH5hJE4EpdPCYmLMCJxJC6SHxK3CHokjsQFc0XiSFwwQonEkbhohkgciYtmyMu++wmSOBIXPce6+he9YmARGzLmZRpO2X30VOLyqSZxJC6fEpW4CCOVuDw9Ekfi8ilRiYswInEkLpIflbhF0FOJU4kL5orEkbhghFTiMgBJXD5hJI7E5VNC4iKMSByJi+SHxC2CHokjccFckTgSF4wQiSNx0QglEkfioiEynLKfIIkjcdFzrKu/4ZRzoGo4ZTdEwynz4SJxJC6fkv4Wnonr56MSl08YiSNx+ZSoxEUYkTgSF8mPStwi6KnEqcQFc0XiSFwwQipxKnHRCKnEZQju3bMznb9vTzr+/okw62VdgEqcSlxJti+7aF9JN31OE1CJm0MUVOJU4kpjROJIXGl26n4qcSpx0QypxKnERTNE4khcSYZIXAm1M31IXIzfpDeJI3GlMSJxJK40OyRuGDnDKfOcSByJy6ekvwWJI3ElGSJxJdRIXIxaqzeJI3GlgSJxJK40OyRuGDkSl+dE4khcPiUkLsLIM3Hd9EhcJFUpqcTF+E16kzgSVxojEkfiSrND4oaRI3F5TiSOxOVTQuIijEgciYvkZ1pfEjcHqiSOxJXGiMSRuNLskLhh5EhcnhOJI3H5lJC4CCMSR+Ii+SFxi6B3epkkjsSVxovEkbjS7JC4YeRIXJ4TiSNx+ZSQuAgjEkfiIvkhcYugR+J6qXpPXD50JI7E5VPS38LslP18SFw+YSSOxOVTQuIijEgciYvkh8Qtgh6JI3HBXJE4EheMkPfEZQCSuHzCSByJy6eExEUYkTgSF8kPiVsEPRJH4oK5InEkLhghEkfiohHysu8MQS/7zkfMKwb6GZE4Epc/i2ZvYWKT2Zlt6uGZuG6IhlPmw0XiSFw+Jf0tDKfs56MSl0+YSpxKXD4lKnERRiSOxEXyoxK3CHoqcSpxwVyROBIXjJBKnEpcNEIqcSpx4QypxKnElYTIe+JKqJ3poxIX4zfprRKnElcaIxJH4kqzU/dTiVOJi2ZIJU4lLpohEkfiSjJE4kqokbgYtVZvEkfiSgNF4khcaXZI3DByhlPmOZE4EpdPSX8LEkfiSjJE4kqokbgYNRI3iJ9n4vKYSByJy6ekv4VKnEpcNEMkjsRFM0TiSFxJhkhcCTUSF6NG4gbxI3F5TCSOxOVTQuIijFTi8vRIHInLp0QlLsLIxCbd9EhcJFUpeSYuxm/S23DKbogkLh8uEkfi8ikhcRFGJC5Pj8SRuHxKSFyEEYkjcZH8TOtL4uZAlcSRuNIYkTgSV5qdup/hlP0ESVw+YSSOxOVTQuIijEgciYvkh8Qtgt7pZZI4ElcaLxJH4kqzQ+KGkSNxeU4kjsTlU0LiIoxIHImL5IfELYIeieulajhlPnQkjsTlU9LfQiVOJS6aIRJH4qIZMrFJP0ESR+Ki51hXf8Mp50BVJU4lrjRGJI7ElWZHJW4YOZW4PCcSR+LyKVGJizAicSQukh+VuEXQU4lTiQvmisSRuGCEkkqcSlw0QySOxEUzpBKnEleSIbNTllA700clLsZv0lslTiWuNEYkjsSVZkclbhg5lbg8JxJH4vIpUYmLMFKJU4mL5GcpK3EPPvZ0euXnx9LRI4emsrnxK99Ml372wg1t9h88nF479sakzzVXXdnb/6WXX00H731o0/Jff+Hx9e9IHIkrPTlJHIkrzQ6JG0aOxOU5kTgSl08JiYswInEkLpKfpZK4p559Pt33yBNZCasErvo0Je6uQ4+mt985vi5uVZsvXPu59PDhOzsZ1RLXlLZ2QxJH4kpPThJH4kqzQ+KGkSNxeU4kjsTlU0LiIoxIHImL5GepJK7embaQNXeyqrbdessN6R//6ZcbqnVX33QgHXng7nTD9ddOmldC+L0f/Ci9+Mx3SdycE2Z2yjxQEkfi8inpb+GZuH4+JC6fMBJH4vIpIXERRiSOxEXys1ISVwncdZ+/Kt1zx22pPeSyLXG5SlvXcMp2VU4lrjteJC5/ypI4EpdPCYmLMCJxeXokjsTlU0LiIoxIHImL5GdlJK6qzlWfenhkW+Iqwbvs0ovW/56TuDa4qn/1aT6H99s/fLyIY7MUyzz/vD3pd/hMPZZ7z9mVPvnk0/TxyU+W4njPeyd27NiRzjt3V/r9ByfnveilWd65e3elkyc/TSdPyVD7oFbn1u7dO9O5a+fZHz6UoWmh37d3d/ro5Kl06tSnS3NezHNHdu3akc7ZvSt9cEKGpnE979zd6cOPTk3+PfPZTGD3rp1r16Id6cMTp+BpEPjM2j2iTzmBUc9O2TWcsnrG7Vfv/WYTkYsvvGB9yGRVjWt+mn/LoeySvt+6wezEVlXizt+3O+EzPVXn7tmZqvsmEtfNaOdaiM5bu8H8nRvwqSGqfuGtBO5jN+CbGFU3lFUlrpK437t5mpqh89Z+CPjo40/SSTfgnYx2r12Izlm7Vv9BhqZm6I8qiVuTXJehbkTVdagSuQ/WRNfnDIHPrN0j+pQTWDqJa6PIzWBZ/f3td341dWKT9vK6JM5wyu4AGk6ZPzENp+xntGvt5uniC/amd3/9YR7mirbwTFz/gTecMn9iGE7Zz2jvmsCdv29POv7+iTzMFW3hPXH9B95wym4+3hMXu2CstMTVQvbck/enKy6/ZEKykrrHj/441c+9VdW+u/7qL9f/3jWbJYkjcaWnIYkjcaXZqfuROBIXzRCJI3HRDJE4EleSIRJXQu1Mn1FKXPMVA/WuNGecbCJpV+LaE5W0JylpS1z93/Uyu94rR+JIXOlpSOJIXGl2SNwwcipxeU4kjsTlU9LfgsSRuJIMkbgSaiOXuNguz783iSNxpakicSSuNDskbhg5EpfnROJIXD4lJC7CyHDKbnokLpKqlEZZiYvt8vx7kzgSV5oqEkfiSrND4oaRI3F5TiSOxOVTQuIijEgciYvkZ1pfEjcHqiSOxJXGiMSRuNLskLhh5EhcnhOJI3H5lJC4CCMSR+Ii+SFxi6B3epkkjsSVxovEkbjS7JC4YeRIXJ4TiSNx+ZSQuAgjEkfiIvkhcYugR+J6qXrFQD50JI7E5VPS38LslP18SFw+YSSOxOVTQuIijEgciYvkh8Qtgh6JI3HBXJE4EheMUCJxJC6aIRJH4qIZMjtlP0ESR+Ki51hXf8/EzYGq4ZTdEFXi8uEicSQunxKVuAgjlbg8PRJH4vIpUYmLMCJxJC6SH5W4RdBTiVOJC+aKxJG4YIRU4jIASVw+YSSOxOVTQuIijEgciYvkh8Qtgh6JI3HBXJE4EheMEIkjcdEIJRJH4qIhMpyynyCJI3HRc6yrv+GUc6BqOGU3RMMp8+EicSQun5L+Fp6J6+ejEpdPGIkjcfmUqMRFGJE4EhfJj0rcIuipxKnEBXNF4khcMEIqcSpx0QipxGUI7t2zM52/b086/v6JMOtlXYBKnEpcSbYvu2hfSTd9ThNQiZtDFFTiVOJKY0TiSFxpdup+KnEqcdEMqcSpxEUzROJIXEmGSFwJtTN9SFyM36Q3iSNxpTEicSSuNDskbhg5wynznEgcicunpL8FiSNxJRkicSXUSFyMWqs3iSNxpYEicSSuNDskbhg5EpfnROJIXD4lJC7CyDNx3fRIXCRVKanExfhNepM4ElcaIxJH4kqzQ+KGkSNxeU4kjsTlU0LiIoxIHImL5GdaXxI3B6okjsSVxojEkbjS7JC4YeRIXJ4TiSNx+ZSQuAgjEkfiIvkhcYugd3qZJI7ElcaLxJG40uyQuGHkSFyeE4kjcfmUkLgIIxJH4iL5IXGLoEfieql6T1w+dCSOxOVT0t/C7JT9fEhcPmEkjsTlU0LiIoxIHImL5IfELYIeiSNxwVyROBIXjJD3xGUAkrh8wkgcicunhMRFGJE4EhfJD4lbBD0SR+KCuSJxJC4YIRJH4qIR8rLvDEEv+85HzCsG+hmROBKXP4tmb2Fik9mZberhmbhuiIZT5sNF4khcPiX9LQyn7OejEpdPmEqcSlw+JSpxEUYkjsRF8qMStwh6KnEqccFckTgSF4yQSpxKXDRCKnEqceEMqcSpxJWEyHviSqid6aMSF+M36a0SpxJXGiMSR+JKs1P3U4lTiYtmSCVOJS6aIRJH4koyROJKqJG4GLVWbxJH4koDReJIXGl2SNwwcoZT5jmROBKXT0l/CxJH4koyROJKqJG4GDUSN4ifZ+LymEgcicunpL+FSpxKXDRDJI7ERTNE4khcSYZIXAk1EhejRuIG8SNxeUwkjsTlU0LiIoxU4vL0SByJy6dEJS7CyMQm3fRIXCRVKXkmLsZv0ttwym6IJC4fLhJH4vIpIXERRiQuT4/Ekbh8SkhchBGJI3GR/EzrS+LmQJXEkbjSGJE4Eleanbqf4ZT9BElcPmEkjsTlU0LiIoxIHImL5IfELYLe6WWSOBJXGi8SR+JKs0PihpEjcXlOJI7E5VNC4iKMSByJi+SHxC2CHonrpWo4ZT50JI7E5VPS30IlTiUumiESR+KiGTKxST9BEkfioudYV3/DKedAVSVOJa40RiSOxJVmRyVuGDmVuDwnEkfi8ilRiYswInEkLpIflbhF0FOJU4kL5orEkbhghJJKnEpcNEMkjsRFM6QSpxJXkiGzU5ZQO9NHJS7Gb9JbJU4lrjRGJI7ElWZHJW4YOZW4PCcSR+LyKVGJizBSiVOJi+RHJW4R9FTiVOKCuSJxJC4YIZW4DEASl08YiSNx+ZSQuAgjEkfiIvkhcYugR+JIXDBXJI7EBSNE4khcNEKJxJG4aIgMp+wnSOJIXPQc6+pvOOUcqBpO2Q3R7JT5cJE4EpdPSX8Lz8T181GJyyeMxJG4fEpU4iKMSByJi+RHJW6NwIOPPZ1e+fmxdPTIoaksb/zKN9Oln71wQ5v9Bw+n1469MelzzVVXbupP4khc6clJ4khcaXbqfiSOxEUzROJIXDRDKnH9BEkciYueY139V6IS99Szz6f7HnliqoTVYCqBqz5Nibvr0KPp7XeOr4tb1eYL134uPXz4znWeJI7ElZ6cJI7ElWaHxA0jpxKX50TiSFw+Jf0tSByJK8mQ2SlLqJ3psxISV+9uW8ia6Kpq26233JD+8Z9+uaFad/VNB9KRB+5ON1x/7aR5JYTf+8GP0ovPfJfEZbJnOGX+5CRxJC6fkv4WKnH9fEhcPmEkjsTlU0LiIoxU4rrpkbhIqlIicWv8KoG77vNXpXvuuG3TkMu2xL308qvp4L0PpddfeJzEkbjY2bfWm8SRuGiISByJi2aIxJG4aIZU4voJkjgSFz3HuvqvvMRV1bnqUw+PbD83VwneZZdetP53Ejc8hipxeVYkjsTlU6ISF2GkEpenR+JIXD4lKnERRiSOxEXyM63vyktc9Yzbr977zSY+F194wfqQyaoa1/w0/7aIg2KZCCCAAAIIIIAAAggggACJWyPQ90xcDSg3g2X197ff+dWGiU3ECwEEEEAAAQQQQAABBBDYKgIrX4lrg+6TuHoo5XNP3p+uuPySrTpG1oMAAggggAACCCCAAAIIrBNYCYlrvmKg3vPmjJPNPLQlrha3uk1zQhM5QgABBBBAAAEEEEAAAQS2msBKSNwioFZDM3/y4iuTRXe9AHwR6zybl9l+trAtyW2RXmUZrlm0GTUZHth/82S21FX7NHPSfvZUhv5lJt3Xjr2xHovmefTmW++mW7727fW/Tfuhahkz1ZxhuLl/ucysynW8a0KuilNz/6v/bl93cvyWKUvTMlTvY31+tRmtynV7WoYqPu1rT3O0kgz1X7crfquSoWW6Xpwt+0LiCo5E+11xuYt/wSpG1aW6gD/8/R+uPydYX7TrG8z6Al9f2HPPHY5q52fc2Do71WQ6zZvs9iyp7VdbzLiaUTavcvH40R9veH1H+wZqlTPUfqa3/d/VjcA3vv7l9NVbv5j6brhGGY4pG92UkPbNde66swrX8fbNdfvHs+rfrqNHDk3o1m3r61KO37LkqC9D7etP9cPSX/z5n67/wLYK1+1chuprTdePRjK0eS6Grut4lbN6hvRV/Ld/Wa4l27EfJK6Aelvaul4AXrDYpemSu3C3/740O57Zkeofu//y3//nZNbT9oW6/d/tm4NVYFQxmPa8aVv8VzFDfdedLmlrSt2y56fa1+bNdbW/ucys0nW8/cPatDw0M5Pjt2yZ6spQvY/1teme//q99XfKVn9bpev2tAxV59Gtt9ww+fGo/ZGhje8hrvi07xdXKUPLds04G/aHxBUchfbN0ar86j0UVZtHl5Cs2q9NbSbN/e8SklWrVtaC23zdx5duvG7910kZSuvVtZpL8+ap64ekVRoh0HUDnsvMKl3Hh0pc87qU4zf034OxtJsmcU0mzXNq1a7b0zJU8akqlPW1uzkMXob6r9urlqGxXAvGtJ0kruBoVRet73zr9vVfnsxauRFi+x/D9gvT618wmwwLDsNounRdqJs3Bl35qSTur//m79bfVTianS3c0K4bhIpRPURu1TPUrAg0b5jq4XFdeeliVnh4zvpuXTfgucys0nV8iMS1h3nl+J31oZhxA7sy1Bb9psSt2nW7K0PtIbgV8opR9amG6crQv4SwLbr1dXvVMjTjKan5AAIkbgCkdpNV+gV3VjwVmy9c+7kN79FbtV/j2sy6Zket21Qi+x//wzWTCSmaQwlXrRLXVUlq5mbVM1TfHDWHLTWfIVSJ2zycMpeZVbqO5ySuYvXTV3+x4UejHL9Z/20429u3Ja79LFhz+6tq+F1/9Zcrdd3uk7j2RCbf+8GPJlla9QzlrtsqcWf7VeHs3z4SV3CMVulZilnwdAlc1X/VxsUPYeaZuI2UuoYkN28AZGjz8zfNX3HffOv/pYP3PrRhUhjPxD2dXvn5sU0Td9Q3nKt0He+TuC6BW8Xrdt8zcfXVqp2ZVXqeqW84ZXNSk+YPSqt23e7KUDsj7erbKmVoyL2RNrMRIHGz8Zq0XoVZzWbF0hz61u6bm+hk1nUtQ/vchXvVnhmsjmnzH8BVnSmvL9vNYUpVu/bN9yrOTlnz6rp5yl13Vuk63jcpRcWwnqGymb8cv2W4Djf3oUTiVmF2yprRtAx1zZJbTzIkQxuHl3Zdt1cpQ8t2zTgb9ofEFR6FVXm/0BA804YLNiemWKV3xQxh1iVp3hXzL9Wm+tN+ZlKGNr5PqP0evVV8T1z7PWdVdprT6Ocys+zX8a4hgfVzptOGCzZzleM35Fp3trfJZai5/V2TBS37dbsvQzWb5vsrm//uV3+Xof7rdv0DZj0xzKq+I/Zsv06crdtH4s7WI2O7EEAAAQQQQAABBBBAAIEOAiROLBBAAAEEEEAAAQQQQACBEREgcSM6WDYVAQQQQAABBBBAAAEEECBxMoAAAggggAACCCCAAAIIjIgAiRvRwbKpCCCAAAIIIIAAAggggACJkwEEEEAAAQQQQAABBBBAYEQESNyIDpZNRQABBBBAAAEEEEAAAQRInAwggAACCCCAAAIIIIAAAiMiQOJGdLBsKgIIIIAAAggggAACCCBA4mQAAQQQQAABBBBAAAEEEBgRARI3ooNlUxFAAAEEEEAAAQQQQAABEicDCCCAAAIIIIAAAggggMCICJC4ER0sm4oAAggggAACCCCAAAIIkDgZQAABBBBAAAEEEEAAAQRGRIDEjehg2VQEEEAAAQQQQAABBBBAgMTJAAIIIIAAAggggAACCCAwIgIkbkQHy6YigAACCCCAAAIIIIAAAiROBhBAAAEEEMgQePCxp9Nf/83fpRef+S5WCCCAAAIIbDsBErfth8AGIIAAAuMlcONXvpl+9d5vNu3A6y88vv7dXYceTT999RejFiASN96M2nIEEEBgGQmQuGU8qvYJAQQQ2CIClcT9xZ//abrnjtvW17j/4OH0zi/fG7W0tfGRuC0KlNUggAACCAwiQOIGYdIIAQQQQKCLQJfEtYVn2n9/4+tfTvc98sT6YpvVu651VXJ42aUXTf70kxdfmfzvNVddmY4eOTT5/2++9W665WvfTkceuDvdcP2164u4+qYD6Tvfuj199dYvTr6r/7u57gP7b077v/xnk/71p9mn3odKWB8/+uOp21y1a/69axnN/W7+XcIQQAABBBAYSoDEDSWlHQIIIIDAJgJdEtf+rkviKtH50o3XpYcP3zlZZiVo1acWsi7UVZvXjr2xScgqAasqgbNIXLX8Whqfevb5iUxefOEF69XD+ru6TS1n9bq6trm9n+3tqZfR3G+RQgABBBBAoIQAiSuhpg8CCCCAwITAtGfiuipQ9aQgXUMTK2n63g9+1DsEs67E1eJXrb963q76VN/NInHN7evqV3/33JP3pysuvyR1bfNLL7+aDt77UKrbVBW++v/X8WhunyGZThoEEEAAgXkRIHHzImk5CCCAwAoS6KrE1QJUV61ywysrbO3KVxfKaRL39jvHJxW8RUhcPTSzS8Ca66u2txK6rk895JPEreAJYpcRQACBBREgcQsCa7EIIIDAKhDokrhqv6sKVC1XqyRx7UpcMwMkbhXOCPuIAAIIbA0BErc1nK0FAQQQWEoCZ5PEVYCrIY1DJzapJzrpG07ZV4mrh1PWz821J1BpH3ASt5SngJ1CAAEEtoUAidsW7FaKAAIILAeBLomr5aZ+7myrKnEV0Wp7vnDt59YnTKkqgtVMll2zU0Ylbtq6mtW4isX/+t8vTraHxC1H5u0FAgggcDYQIHFnw1GwDQgggMBICUyb2KRZDdtKiaurajXOSt6qmSfnJXHNF5t3zTLZfsVAtR191byRHnabjQACCCCwzQRI3DYfAKtHAAEEEEAAAQQQQAABBGYhQOJmoaUtAggggAACCCCAAAIIILDNBEjcNh8Aq0cAAQQQQAABBBBAAAEEZiFA4mahpS0CCCCAAAIIIIAAAgggsM0ESNw2HwCrRwABBBBAAAEEEEAAAQRmIUDiZqGlLQIIIIAAAggggAACCCCwzQRI3DYfAKtHAAEEEEAAAQQQQAABBGYhQOJmoaUtAggggAACCCCAAAIIILDNBEjckdwrxwAAAq9JREFUNh8Aq0cAAQQQQAABBBBAAAEEZiFA4mahpS0CCCCAAAIIIIAAAgggsM0ESNw2HwCrRwABBBBAAAEEEEAAAQRmIUDiZqGlLQIIIIAAAggggAACCCCwzQRI3DYfAKtHAAEEEEAAAQQQQAABBGYhQOJmoaUtAggggAACCCCAAAIIILDNBEjcNh8Aq0cAAQQQQAABBBBAAAEEZiFA4mahpS0CCCCAAAIIIIAAAgggsM0ESNw2HwCrRwABBBBAAAEEEEAAAQRmIUDiZqGlLQIIIIAAAggggAACCCCwzQRI3DYfAKtHAAEEEEAAAQQQQAABBGYhQOJmoaUtAggggAACCCCAAAIIILDNBEjcNh8Aq0cAAQQQQAABBBBAAAEEZiFA4mahpS0CCCCAAAIIIIAAAgggsM0ESNw2HwCrRwABBBBAAAEEEEAAAQRmIUDiZqGlLQIIIIAAAggggAACCCCwzQRI3DYfAKtHAAEEEEAAAQQQQAABBGYhQOJmoaUtAggggAACCCCAAAIIILDNBEjcNh8Aq0cAAQQQQAABBBBAAAEEZiFA4mahpS0CCCCAAAIIIIAAAgggsM0ESNw2HwCrRwABBBBAAAEEEEAAAQRmIUDiZqGlLQIIIIAAAggggAACCCCwzQRI3DYfAKtHAAEEEEAAAQQQQAABBGYhQOJmoaUtAggggAACCCCAAAIIILDNBEjcNh8Aq0cAAQQQQAABBBBAAAEEZiFA4mahpS0CCCCAAAIIIIAAAgggsM0ESNw2HwCrRwABBBBAAAEEEEAAAQRmIUDiZqGlLQIIIIAAAggggAACCCCwzQRI3DYfAKtHAAEEEEAAAQQQQAABBGYhQOJmoaUtAggggAACCCCAAAIIILDNBEjcNh8Aq0cAAQQQQAABBBBAAAEEZiHw/wEaEfogcTtFLgAAAABJRU5ErkJggg==",
"text/html": [
"