{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Shan-Chen Two-Phase Single-Component Lattice Boltzmann" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from lbmpy.session import *\n", "from lbmpy.updatekernels import create_stream_pull_with_output_kernel\n", "from lbmpy.macroscopic_value_kernels import macroscopic_values_getter, macroscopic_values_setter\n", "from lbmpy.maxwellian_equilibrium import get_weights" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This is based on section 9.3.2 of Krüger et al.'s \"The Lattice Boltzmann Method\", Springer 2017 (http://www.lbmbook.com).\n", "Sample code is available at [https://github.com/lbm-principles-practice/code/](https://github.com/lbm-principles-practice/code/blob/master/chapter9/shanchen.cpp)." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Parameters" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "N = 64\n", "omega_a = 1.\n", "g_aa = -4.7\n", "rho0 = 1.\n", "\n", "stencil = LBStencil(Stencil.D2Q9)\n", "weights = get_weights(stencil, c_s_sq=sp.Rational(1,3))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Data structures" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "dh = ps.create_data_handling((N,) * stencil.D, periodicity=True, default_target=ps.Target.CPU)\n", "\n", "src = dh.add_array('src', values_per_cell=stencil.Q)\n", "dst = dh.add_array_like('dst', 'src')\n", "\n", "ρ = dh.add_array('rho')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Force & combined velocity" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The force on the fluid is\n", "$\\vec{F}_A(\\vec{x})=-\\psi(\\rho_A(\\vec{x}))g_{AA}\\sum\\limits_{i=1}^{q}w_i\\psi(\\rho_A(\\vec{x}+\\vec{c}_i))\\vec{c}_i$\n", "with \n", "$\\psi(\\rho)=\\rho_0\\left[1-\\exp(-\\rho/\\rho_0)\\right]$." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "def psi(dens):\n", " return rho0 * (1. - sp.exp(-dens / rho0));" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "zero_vec = sp.Matrix([0] * stencil.D) \n", "\n", "force = sum((psi(ρ[d]) * w_d * sp.Matrix(d)\n", " for d, w_d in zip(stencil, weights)), zero_vec) * psi(ρ.center) * -1 * g_aa" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Kernels" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "lbm_config = LBMConfig(stencil=stencil, relaxation_rate=omega_a, compressible=True,\n", " force_model=ForceModel.GUO, force=force, kernel_type='collide_only')\n", "\n", "collision = create_lb_update_rule(lbm_config=lbm_config,\n", " optimization={'symbolic_field': src})\n", "\n", "stream = create_stream_pull_with_output_kernel(collision.method, src, dst, {'density': ρ})\n", "\n", "\n", "config = ps.CreateKernelConfig(target=dh.default_target, cpu_openmp=False)\n", "\n", "stream_kernel = ps.create_kernel(stream, config=config).compile()\n", "collision_kernel = ps.create_kernel(collision, config=config).compile()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Initialization" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "method_without_force = create_lb_method(LBMConfig(stencil=stencil, relaxation_rate=omega_a, compressible=True))\n", "init_assignments = macroscopic_values_setter(method_without_force, velocity=(0, 0), \n", " pdfs=src.center_vector, density=ρ.center)\n", "\n", "\n", "init_kernel = ps.create_kernel(init_assignments, ghost_layers=0, config=config).compile()" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "def init():\n", " for x in range(N):\n", " for y in range(N):\n", " if (x-N/2)**2 + (y-N/2)**2 <= 15**2:\n", " dh.fill(ρ.name, 2.1, slice_obj=[x,y])\n", " else:\n", " dh.fill(ρ.name, 0.15, slice_obj=[x,y])\n", "\n", " dh.run_kernel(init_kernel)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Timeloop" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "sync_pdfs = dh.synchronization_function([src.name])\n", "sync_ρs = dh.synchronization_function([ρ.name])\n", "\n", "def time_loop(steps):\n", " dh.all_to_gpu()\n", " for i in range(steps):\n", " sync_ρs()\n", " dh.run_kernel(collision_kernel)\n", " \n", " sync_pdfs()\n", " dh.run_kernel(stream_kernel)\n", " \n", " dh.swap(src.name, dst.name)\n", " dh.all_to_cpu()" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "def plot_ρs():\n", " plt.figure(dpi=200)\n", " plt.title(\"$\\\\rho$\")\n", " plt.scalar_field(dh.gather_array(ρ.name), vmin=0, vmax=2.5)\n", " plt.colorbar()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Run the simulation\n", "### Initial state" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAACQcAAAQKCAYAAAD3zl6bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAewgAAHsIBbtB1PgAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzde7BlZ1kn4N97cki4ChjA3FAuwQmRqEHIBCloRUFQQJJCS1AUAQVnJqWUpaQEB7zMjChThRMZhCGCMirKJZIZCYgXWkQYwpSiYwggiiYESAJyD+mw9zt/nN32tu1zduf02Weftfp5qnadb6317e97z+rOP51fvV91dwAAAAAAAAAAgPFZW3UBAAAAAAAAAADAcggHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAHtCVd2xqp5bVe+squur6ubZ58NVdWlVfc2qawQAAACAoanuXnUNAAAAwHGuqh6T5H8kOX2LaV9M8pTuft3uVAUAAAAAwyccBAAAAKxUVb0gyfPnbn0wyZVJbk7ytUm+Ye7ZF5I8sLvfv2sFAgAAAMCACQcBAAAAK1NV/zHJz8wur03yQ9395sPmPDbJ65KcNLv1qu7+wd2rEgAAAACGSzgIAAAAWImqekSSP0xSSf4xyTd290c2mfvTSX52dvnZJHfp7umuFAoAAAAAAyYcBAAAAOy6qjohyVVJvjpJZyMY9K4t5p+R5Jq5W2c5WgwAAAAAFltbdQEAAADAcem7sxEMSpLf2SoYlCTdfW2Sz8zduvOyCgMAAACAMREOAgAAAFbhB+fGLznK79wyN9YKGQAAAACOgmPFAAAAgF1VVXdK8k9JTkhyY5J79IJ/oKiq9SQ3JVmf3Tq1uz+21EIBAAAAYAR0DgIAAAB22zdmIxiUJH+2KBg0c2YOBYM+m+TjyygMAAAAAMZGOAgAAADYbf92bvy+o/zOw+bGf3qUgSIAAAAAOO4JBwEAAAC77dy58fVH+Z0L5sZv2cFaAAAAAGDUhIMAAACA3TYfDrrNoslVdWaSb5tdfjHJby6jKAAAAAAYI+EgAAAAYNdU1Zcn+aq5W2cfxdf+aw79G8Yru/uTO14YAAAAAIyUcBAAAACwm8497PqCqrrrZpOr6uIkj59dfjLJTy+rMAAAAAAYI+EgAAAAYDfNh4MOJLlzktdU1V3mJ1XVnavqJUn+y+zWJMlTuvsTu1MmAAAAAIzD+qoLAAAAAI4r8+Gg5yb5hSSPSvIPVfVHSW5MckaSfUluP5s3SfKs7n7TbhYKAAAAAGNQ3b3qGgAAAIDjRFVdleT+s8t7ZuPIsEuyeXfja5P8UHe/eRfKAwAAAIDREQ4CAAAAdkVV3T7JZ5KckOST3X3y7P5Dk/xokvOSnJrk00muSvL6JJd29xdWUzEAAAAADJ9jxQAAAIDd8nXZCAYlyXsP3uzudyR5x0oqAgAAAICR26xlNwAAAMBOO3du/JcrqwIAAAAAjiPCQQAAAMBumQ8HvXfTWQAAAADAjhEOAgAAAHaLcBAAAAAA7LLq7lXXAAAAAIxcVa0n+VySk5LckuSO3X1gtVUBAAAAwPjpHAQAAADshrOzEQxKkqsFgwAAAABgdwgHAQAAALvBkWIAAAAAsAKOFQMAAAAAAAAA4LhQVQ9M8ugkD0vygCT3SHJLkuuS/HmSS7v77TuwzwuSPP8op39zd7/tWPfczPqyFgYAAAAAAAAAgL2iqvYnefgRHp2Y5H6zzw9U1auTPKO7D+xmfcsiHAQAAAAAAAAAwPHg9NnP65K8Nsnbk/xjkhOSPCTJj8/mPCUbmZon79C+5yx4/vc7tM8ROVYMAAAAAAAAAIDRq6r/neQ3kry+uydHeH63JO9I8tWzWw/f7hFj88eKdXdtq+AdsrbKzQEAAAAAAAAAYDd092O7+3ePFAyaPb8xG92DDnri7lS2XMJBAAAAAAAAAACw4W1z4/uuqoidJBwEAAAAAAAAAAAbTpwbT1dWxQ4SDgIAAAAAAAAAgA375sZX78SCVfXWqvpEVR2oquur6m1VdXFV3XUn1l+4f3fvxj6DVlUnJTlndnlDkiOePQcAAAAAAAAAe9wJSe4+G/91d9+8ymJYrqpaT3LKquvYplNyFBmN7r52pzasqrUk70xy3uzWg7v7Pdtc6wVJnr9g2qeSPLW737idPY7W+jIXH5Fzkly56iIAAAAAAAAAYAc9OMm2gg8MxilJrll1EUtWO7jWs3MoGHTZdoNBc/46ye8leXeS65LcJsm/SfK9SR6V5C5JXl9Vj+vuK45xr03pHHQUqupBEQ4CAAAAAAAAYFy23RWFYaiqMzLycFB370g4qKr2JfnDbDTauT7J13b3x49hvbt096e2eP7MJL86u7wuyZndfdN299uKzkFH54aDg3e/+9059dRTV1kLADBS3/f1P7XqEgAAANhD/udf/udVlwAAjNBHP/rRnHfewcYoh/5fOOP3rivumVPvccKqy1joo9dPcv5j/jnP9OAkH1v2nlX1NUkuy0aO5uYk330swaAk2SoYNHv+slmzmmckOS3JhUl+81j23Ixw0NH55/PrTj311JxxxhmrrAUAGKnbnnDHVZcAAADAHuLfogGAXTBZPIWxOPUeJ+SM026z6jJurY9197XL3KCq7p3kD5LcNRv/TTypu/cvc885L8tGOChJ9mVJ4aC1ZSwKAAAAAAAAAAB7WVWdlo2jxE5L0kme1t2X7WIJV82NT1/WJjoHAQAAAAAAAACM3DSdaaarLmOhaXpX9qmquyV5a5L7zG5d1N2/sSubz5WxG5voHAQAAAAAAAAAwHGjqu6c5C1Jzp7duri7X7KCUs6eG1+3rE2EgwAAAAAAAAAAOC5U1e2T/H6SB85u/afufuGKynnm3Hj/sjZxrBgAAAAAAAAAwMhNeppJ7/1jxZZZY1WdmOSyJA+d3frl7n7eNtZ5apJXzi5/prtfcNjzc5Lc1N1/u8Uaz0zy9Nnlx2Z1LYVwEAAAAAAAAAAAx4PfTvKo2fiPk1xaVQ/YYv6B7v7ANvb5hiSvqKo/SXJFkr9O8ols5HTOSvJ9SR45mztJ8szu/vw29jkqwkEAAAAAAAAAABwPLpwbPyLJXy2Y/w9J7rXNvU5I8q2zz2Y+keTp3X35Nvc4KsJBAAAAAAAAAACwc96UjSPDHpLk3CRfkeTkJJXkk0nem+TNSV7V3Z9ZdjHCQQAAAAAAAAAAIzdNZ5pedRkLLbPG7q4dWudVSV61xfPrk/za7LNya6suAAAAAAAAAAAAWA7hIAAAAAAAAAAAGCnHigEAAAAAAAAAjFxnmmmmqy5joR5AjUOjcxAAAAAAAAAAAIyUcBAAAAAAAAAAAIyUcBAAAAAAAAAAAIzU+qoLAAAAAAAAAABguSbdmXSvuoyFhlDj0OgcBAAAAAAAAAAAIyUcBAAAAAAAAAAAI+VYMQAAAAAAAACAkZumM83eP7JrCDUOjc5BAAAAAAAAAAAwUsJBAAAAAAAAAAAwUsJBAAAAAAAAAAAwUuurLgAAAAAAAAAAgOWaJpmkV13GQtNVFzBCOgcBAAAAAAAAAMBICQcBAAAAAAAAAMBIOVYMAAAAAAAAAGDkpulMB3Gs2N6vcWh0DgIAAAAAAAAAgJESDgIAAAAAAAAAgJESDgIAAAAAAAAAgJFaX3UBAAAAAAAAAAAs16Q7k+5Vl7HQEGocGp2DAAAAAAAAAABgpISDAAAAAAAAAABgpBwrBgAAAAAAAAAwctPZZ68bQo1Do3MQAAAAAAAAAACMlHAQAAAAAAAAAACMlHAQAAAAAAAAAACM1PpubVRVd0vytCTfmeS+Se6a5BNJrknyp0ne0N3vXLDGo5P8cJLzktw9yQ1J3p3k5d395uVVDwAAAAAAAAAwXNN0JulVl7HQdAA1Ds2uhIOq6ruSvDTJyYc9OnX2OS/J/ZI8YZPvV5JfzUYwaN7pSS5IckFVvTzJs7rb3xIAAAAAAAAAAMguhIOq6vuTvDIbR5hdn42Q0J8l+WSSU7LRRehxSW7ZYpmfz6Fg0F8k+cUkH5p99yeTnDt7fkOS5+34LwEAAAAAAAAAAAO01HBQVd0/ycuzEQx6e5LHdfenjzD1kqo6cZM1zsxGAChJ3pPk4d190+z6yqq6PMn+JA9K8pyqemV3f2gnfw8AAAAAAAAAABiitSWvf0mSk5LcmOTCTYJBSZLuPrDJo2fnUIjporlg0MHvfSHJRbPL9SQ/dkwVAwAAAAAAAACMzCTJpAfwWfWLGqGlhYOq6qwk3zK7/JXuvnEba1SS75xdXt3d7zrSvNn9988unzD7HgAAAAAAAAAAHNeW2Tnou+bGrz04qKq7VtX9qurko1jj3klOn433L5h78PkZSe51tEUCAAAAAAAAAMBYLTMcdP7s56eTvK+qvreq3pvkk0k+kOTGqvq7qnp+Vd1xkzXuPze+esF+88/vv+ksAAAAAAAAAIDjzHRAH3bW+hLXPnv288NJLkny748w595JXpDkiVX1bd193WHP7zk3vnbBftds8r2FquqMBVNOuTXrAQAAAAAAAADAXrDMcNCXz36eleTrknwqycVJ3pDkM0nOSfKzSR6T5AFJXltVD+vu+RDYnebGn1uw3+fnxpt1ItrMNYunAAAAAAAAAADAsCzzWLE7zH6elGSS5DHd/bLuvqG7b+7u9yR5bJIrZvO+McmFh61x27nxgQX73Tw3vt02awYAAAAAAAAAgNFYZuegL+ZQQOi13f2uwyd097SqfiIb3YOS5ElJXnfYGgeduGC/k+bGN93KWhcdQ3ZKkitv5ZoAAAAAAAAAAHvCNJVJatVlLDQdQI1Ds8xw0GdzKBx0xWaTuvtvquojSU5P8uAjrHHQoqPC7jA3XnQE2eE1XLvV8yp/8QAAAAAAAAAAGJ5lHit2zdx4y/DN3Nx7HHZ//ntnLFhjvvvPNZvOAgAAAAAAAACA48QyOwf9TQ51AjphwdyDz7902P2r5sZnLVhj/vn7FswFAAAAAAAAADhuTHvjs9cNocahWWbnoD+dG993wdz7zH5+5LD7f5/kutl434I1Hj63xocXFQcAAAAAAAAAAGO3zHDQ5UlumY0v3GxSVe1LcvLs8u3zz7q7k7xxdnlWVZ2/yRrn51DnoDfOvgcAAAAAAAAAAMe1pYWDuvsTSV4xu3xkVX3P4XOq6k5JXjx362VHWOrFOXTc2CVVdbvD1rhdkktml186bD0AAAAAAAAAADhurS95/ecn+Y4kX5nk1VX10CRvSPKZJOckeU4Odfx5aXdfefgC3f2BqnpRkouTPCjJO6rqhUk+lI3jyp6T5NzZ9F/q7g8u8fcBAAAAAAAAABicSSqT1KrLWGgINQ7NUsNB3X1DVT06G0eMnZnkP8w+h/u1JD+6xVLPTXKPJE/LRhDoNUeYc2mS5x1TwQAAAAAAAAAAMCJLO1bsoO5+X5KvT/ITSf5Pkk8mOZDk2iS/k+QR3f307r5lizWm3f30bHQhemOS62ZrXDe7/vbufkZ3T5f6ywAAAAAAAAAAwIAs+1ixJEl3fz7Ji2afY1nnTUnetCNFAQAAAAAAAAAcJxwrdvxaeucgAAAAAAAAAABgNYSDAAAAAAAAAABgpISDAAAAAAAAAABgpNZXXQAAAAAAAAAAAMvVnUy7Vl3GQt2rrmB8dA4CAAAAAAAAAICREg4CAAAAAAAAAICRcqwYAAAAAAAAAMDITVKZZO8fKzaEGodG5yAAAAAAAAAAABgp4SAAAAAAAAAAABgp4SAAAAAAAAAAABip9VUXAAAAAAAAAADAck2ylskAesgMocah8UYBAAAAAAAAAGCkhIMAAAAAAAAAAGCkHCsGAAAAAAAAADBy3ZVp16rLWKgHUOPQ6BwEAAAAAAAAAAAjJRwEAAAAAAAAAAAjJRwEAAAAAAAAAAAjtb7qAgAAAAAAAAAAWK5JKpPUqstYaAg1Do3OQQAAAAAAAAAAMFLCQQAAAAAAAAAAMFKOFQMAAAAAAAAAGLlJr2XSe7+HzBBqHBpvFAAAAAAAAAAARko4CAAAAAAAAAAARko4CAAAAAAAAAAARmp91QUAAAAAAAAAALBc01SmA+ghM02tuoTR2ft/6gAAAAAAAAAAwLYIBwEAAAAAAAAAwEgJBwEAAAAAAAAAwEitr7oAAAAAAAAAAACWa5rKJLXqMhaaDqDGodE5CAAAAAAAAAAARko4CAAAAAAAAAAARsqxYgAAAAAAAAAAIzfptUx67/eQGUKNQ+ONAgAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASK2vugAAAAAAAAAAAJZrmso0teoyFhpCjUOjcxAAAAAAAAAAAIyUcBAAAAAAAAAAAIyUY8UAAAAAAAAAAEZumrVMBtBDZjqAGofGGwUAAAAAAAAAgJESDgIAAAAAAAAAgJESDgIAAAAAAAAAgJFaX3UBAAAAAAAAAAAs16TXMum930NmCDUOjTcKAAAAAAAAAAAjJRwEAAAAAAAAAAAj5VgxAAAAAAAAAICRm6YyHUAPmWlq1SWMzt7/UwcAAAAAAAAAALZFOAgAAAAAAAAAAEZKOAgAAAAAAAAAAEZqfdUFAAAAAAAAAACwXNOuTLpWXcZC0wHUODQ6BwEAAAAAAAAAwEgJBwEAAAAAAAAAwEg5VgwAAAAAAAAAYOQmWctkAD1khlDj0HijAAAAAAAAAAAwUsJBAAAAAAAAAAAwUsJBAAAAAAAAAAAwUuurLgAAAAAAAAAAgOWa9lqmvfd7yAyhxqHxRgEAAAAAAAAAYKSEgwAAAAAAAAAAYKQcKwYAAAAAAAAAMHKTrGUygB4yQ6hxaLxRAAAAAAAAAAAYKeEgAAAAAAAAAAAYKeEgAAAAAAAAAAAYqfVVFwAAAAAAAAAAwHJNk0y6Vl3GQtNVFzBCOgcBAAAAAAAAAMBICQcBAAAAAAAAAMBIOVYMAAAAAAAAAGDkplnLdAA9ZIZQ49B4owAAAAAAAAAAMFLCQQAAAAAAAAAAMFLCQQAAAAAAAAAAMFLrqy4AAAAAAAAAAIDlmvRaJr33e8gMocah8UYBAAAAAAAAAGCkhIMAAAAAAAAAAGCkhIMAAAAAAAAAAGCk1lddAAAAAAAAAAAAyzVNZZpadRkLDaHGodE5CAAAAAAAAAAARko4CAAAAAAAAAAARsqxYgAAAAAAAAAAIzfttUx67/eQmQ6gxqHxRgEAAAAAAAAAYKSEgwAAAAAAAAAAYKSEgwAAAAAAAAAAYKTWV10AAAAAAAAAAADLNUllMoAeMpPUqksYnb3/pw4AAAAAAAAAAGyLcBAAAAAAAAAAAIyUY8UAAAAAAAAAAEZu2pVp7/0ju4ZQ49DoHAQAAAAAAAAAACMlHAQAAAAAAAAAACMlHAQAAAAAAAAAACO1vuoCAAAAAAAAAABYrmnWMhlAD5npAGocGm8UAAAAAAAAAABGSjgIAAAAAAAAAABGyrFiAAAAAAAAAAAjN+21THvv95AZQo1D440CAAAAAAAAAMBICQcBAAAAAAAAAMBICQcBAAAAAAAAAMBIra+6AAAAAAAAAAAAlmuSyiS16jIWGkKNQ6NzEAAAAAAAAAAAjJRwEAAAAAAAAAAAjJRjxQAAAAAAAAAARm7aa5n23u8hM4Qah8YbBQAAAAAAAACAkRIOAgAAAAAAAACAkRIOAgAAAAAAAACAkVpfdQEAAAAAAAAAACzXNMkkteoyFpquuoAR0jkIAAAAAAAAAABGSjgIAAAAAAAAAABGyrFiAAAAAAAAAAAjN+21THvv95AZQo1D440CAAAAAAAAAMBICQcBAAAAAAAAAMBICQcBAAAAAAAAAMBIra+6AAAAAAAAAAAAlmvSa5n03u8hM4Qah8YbBQAAAAAAAADguFBVD6yqn6qqK6rqmqq6uao+V1UfqKpXVdXDlrDn91TVW6rqo1X1xar6cFW9uqrO3+m9jkTnIAAAAAAAAAAARq+q9id5+BEenZjkfrPPD1TVq5M8o7sPHON+t03y2iSPPezRV80+T66qF3T3zx3LPosIBwEAAAAAAAAAjFynMk2tuoyFerk1nj77eV02QjtvT/KPSU5I8pAkPz6b85RsZGqefIz7XZpDwaA/SfLLs73PSfJTSe6b5Ger6qPd/Ypj3GtTwkEAAAAAAAAAABwPrs5GKOf13T057Nm7Zh2D3pHkq5M8qape2t1v385GVbUvh8JF/yvJBXN7XllVlyf5v0m+MskvVtXruvtT29lrkbVlLAoAAAAAAAAAAHtJdz+2u3/3CMGgg89vzEb3oIOeeAzb/eTs5yTJvzt8z9lez5ld3jXJ049hry0JBwEAAAAAAAAAwIa3zY3vu50FquqOSb5ldvnW7r52k6lvSPKZ2fjC7ex1NBwrBgAAAAAAAAAwcpNey6T3fg+ZPVDjiXPj6TbXOC/JSbPx/s0mdfeBqnpXkkclOa+qbtPdt2xzz02t/I0CAAAAAAAAAMAesW9ufPU217j/rVjj4PP1JPfb5n5b0jkIAAAAAAAAAIC96JSq2nLCFkd23WpVtZbk4rlbv7vNpe45N15U3zWHfe+qbe65KeEgAAAAAAAAAAD2oiuPYs7W6aFb59nZOBIsSS7r7vdsc507zY0/t2Du5+fGd9zmfltaajioqvoop+7v7m9asNajk/xwNv4Q7p7khiTvTvLy7n7zsdQJAAAAAAAAADBm01SmvZM5muWY7mjW5+hV1b4kvzC7vD7JjxzDcredGx9YMPfmufHtjmHPTe35zkG10R/qV7MRDJp3epILklxQVS9P8qzuPtowEgAAAAAAAAAAe9uDk3xs2ZtU1dckuSwbOZqbk3x3d3/8GJb84tz4xAVzT5ob33QMe25qt8JBL03y37d4/vktnv18DgWD/iLJLyb5UJL7JvnJJOfOnt+Q5HnHXCkAAAAAAAAAAHvBx7r72mVuUFX3TvIHSe6aZJLkSd29/xiX/ezceNFRYXeYGy86gmxbdiscdH13/79b+6WqOjMbAaAkeU+Sh3f3wZTUlVV1eZL9SR6U5DlV9cru/tCOVAwAAAAAAAAAMBKTVCZZW3UZC0128VixqjotyR8mOS1JJ3lad1+2A0vPB5rOyEbmZTP3nBtfswN7/yt7/U/92TkUYLpoLhiUJOnuLyS5aHa5nuTHdrE2AAAAAAAAAAAGqKruluStSe4zu3VRd//GDi1/1dz4rAVzDz7/UpK/3aH9/4U9Gw6qqkrynbPLq7v7XUeaN7v//tnlE2bfAwAAAAAAAACAf6Wq7pzkLUnOnt26uLtfsoNbXJnkwGy8b4s6Tkxy/sHvdPeBzeYeiz0bDkpy7ySnz8aLznI7+PyMJPdaVkEAAAAAAAAAAAxXVd0+ye8neeDs1n/q7hfu5B7d/dkkfzS7/NaqOmOTqRcm+bLZeCeOMzui3QoHfVdVvb+qbqqqz1bVB6vq16vqm7f4zv3nxlcvWH/++f03nQUAAAAAAAAAcBzqrkwH8Ole3oFRs049lyV56OzWL3f387axzlOrqmefF2wy7UWzn+tJXlJVJxy2xt2SHAwlfSrJK25tHUdrfVkLH+bsw67PnH2+v6p+L8lTu/vTh82559z42gXrX7PJ947KFgmtg065tWsCAAAAAAAAALCn/HaSR83Gf5zk0qp6wBbzD3T3B7azUXf/cVW9Jsn3JHl8krdW1YuTXJfknCTPTfKVs+kXd/c/bWefo7HscNAXklyejVZJVyf5XJK7Z+M8tWclOTnJE5K8saoe2d23zH33TnPjzy3Y5/Nz4ztuo85rFk8BAAAAAAAAAGDALpwbPyLJXy2Y/w9J7nUM+z0tG8eGfXuSb5595k2T/Fx3v+wY9lho2eGg07v7U0e4/9aquiTJFUnOzUZY6EeS/Le5ObedGx9YsM/Nc+PbbadQAAAAAAAAAICxmmYt06ytuoyFhlDj0erum5J8R1U9OclTk3xdkrsk+XiStyf5le5+57LrWGo4aJNg0MFnH6+qJyZ5X5ITk1yUfxkO+uLc+MQFW500N77p1taZxUeRnZLkym2sCwAAAAAAAADAHtDdtUPrvCrJq27F/N9K8ls7sfd2LLtz0Ja6+++q6q1JviPJmVV1WndfN3v82bmpi44Ku8PceNERZEeq49qtnlftyN8NAAAAAAAAAADYVXuhF9NVc+PT58bzgZ0zFqwx3/nnmmOuCAAAAAAAAAAARmClnYNmNmvLMx8aOmvBGvPP33ds5QAAAAAAAAAAjMukK5OdOVVrqYZQ49Dshc5BZ8+Nr5sb//3c9b4Fazx89vMjST68M2UBAAAAAAAAAMCwrTQcVFX3SfLI2eXfdfdHDj7r7k7yxtnlWVV1/iZrnJ9DnYPeOPseAAAAAAAAAAAc95YWDqqqx1XVpseWVdVXJHldktvMbr3kCNNenORLs/ElVXW7w9a4XZJLZpdfms0HAAAAAAAAAGDOtGswH3bWpuGdHXBJkttU1euTvDMbx33dlORuSb4pybOSnDyb+2c5Qjiouz9QVS9KcnGSByV5R1W9MMmHktw3yXOSnDub/kvd/cFl/TIAAAAAAAAAADA0ywwHJclpSS6afTbz+iTP6O6bN3n+3CT3SPK0bASBXnOEOZcmed4x1AkAAAAAAAAAAKOzzHDQDyTZl+QhSe6TjY5BX5bkc0muSfLnSX69u9+51SLdPU3y9FkHoh9O8uDZWjcmuTLJy7r7imX9EgAAAAAAAAAAMFRLCwd19/4k+3dwvTcledNOrQcAAAAAAAAAcLzoXsu011Zdxo1MB5cAACAASURBVEI9gBqHxhsFAAAAAAAAAICREg4CAAAAAAAAAICRWtqxYgAAAAAAAAAA7A2TVCapVZex0BBqHBqdgwAAAAAAAAAAYKSEgwAAAAAAAAAAYKSEgwAAAAAAAAAAYKTWV10AAAAAAAAAAADLNe1k2rXqMhaa9qorGB+dgwAAAAAAAAAAYKSEgwAAAAAAAAAAYKQcKwYAAAAAAAAAMHLTXsu0934PmSHUODTeKAAAAAAAAAAAjJRwEAAAAAAAAAAAjJRwEAAAAAAAAAAAjNT6qgsAAAAAAAAAAGC5OpVpatVlLNQDqHFodA4CAAAAAAAAAICREg4CAAAAAAAAAICRcqwYAAAAAAAAAMDITboy6b1/ZNcQahwanYMAAAAAAAAAAGCkhIMAAAAAAAAAAGCkhIMAAAAAAAAAAGCk1lddAAAAAAAAAAAAyzXttUx77/eQGUKNQ+ONAgAAAAAAAADASAkHAQAAAAAAAADASAkHAQAAAAAAAADASK2vugAAAACAeY/ff9W2vnf5vrN3uBIAAACA8ZimMu1adRkLTbP3axwanYMAAAAAAAAAAGCkhIMAAAAAAAAAAGCkHCsGAAAAAAAAADBynRrEkV09gBqHRucgAAAAAAAAAAAYKeEgAAAAAAAAAAAYKeEgAAAAAAAAAAAYqfVVFwAAAAAAAAAAwHJNO5l2rbqMhaa96grGR+cgAAAAAAAAAAAYKeEgAAAAAAAAAAAYKceKAQAAAAAAAACM3LTXMu2930NmCDUOjTcKAAAAAAAAAAAjJRwEAAAAAAAAAAAjJRwEAAAAAAAAAAAjtb7qAgAAAOB49Pj9V626hNHxTjd3+b6zV10CAAAAsGLTrky7Vl3GQkOocWh0DgIAAAAAAAAAgJESDgIAAAAAAAAAgJFyrBgAAAAAAAAAwMhNU5lm7x/ZNYQah0bnIAAAAAAAAAAAGCnhIAAAAAAAAAAAGCnhIAAAAAAAAAAAGKn1VRcAAAAAAAAAAMBydVemXasuY6EeQI1Do3MQAAAAAAAAAACMlHAQAAAAAAAAAACMlGPFAAAAAAAAAABGbjqQY8WGUOPQ6BwEAAAAAAAAAAAjJRwEAAAAAAAAAAAjJRwEAAAAAAAAAAAjtb7qAgAAAGAvePz+q1ZdAizNbv/9vnzf2bu6HwAAALDYtCvTrlWXsdAQahwanYMAAAAAAAAAAGCkhIMAAAAAAAAAAGCkHCsGAAAAAAAAADByjhU7fukcBAAAAAAAAAAAIyUcBAAAAAAAAAAAIyUcBAAAAAAAAAAAI7W+6gIAAAAAAAAAAPj/7N1/sKZnWR/w7xWOIRFwtCJmmqBQYpvslGlpgcGBslqpQxFXQq0j/hgzEAE7kw5pR0LVjhH8oyKd4tAfhEqNdVTGVjCrQ7S0xZ3WFl0cbOmsUZuG6cbIT2WEEMjAXv3jPDt5u+7uu3v2PPu+z30+n5l73vt5n/u5n+/JbM7555rrnlcnOZXadIy1etMBBqRzEAAAAAAAAAAADEpxEAAAAAAAAAAADMqxYgAAAAAAAAAAgzvVlVO9/ceKLSHj0ugcBAAAAAAAAAAAg1IcBAAAAAAAAAAAg1IcBAAAAAAAAAAAg9rZdAAAAAA4myPHTmw6ArBHm/j/9+jhQ5f9nQAAALAkp1I51bXpGGudyvZnXBqdgwAAAAAAAAAAYFCKgwAAAAAAAAAAYFCKgwAAAAAAAAAAYFA7mw4AAAAAAAAAAMC8TnXlVNemY6y1hIxLo3MQAAAAAAAAAAAMSnEQAAAAAAAAAAAMyrFiAAAAAAAAAACDc6zYwaVzEAAAAAAAAAAADEpxEAAAAAAAAAAADEpxEAAAAAAAAAAADGpn0wEAAAAAAAAAAJhZV7pr0ynWW0LGhdE5CAAAAAAAAAAABqU4CAAAAAAAAAAABuVYMQAAAGZz5NiJTUcADoi9/r45evjQPicBAACA7XQqlVPZ/iO7lpBxaXQOAgAAAAAAAACAQSkOAgAAAAAAAACAQSkOAgAAAAAAAACAQe1sOgAAAAAAAAAAAPM61ZVTXZuOsdYSMi6NzkEAAAAAAAAAADAoxUEAAAAAAAAAADAox4oBAAAAAAAAAAyuu9ILOLJrCRmXRucgAAAAAAAAAAAYlOIgAAAAAAAAAAAYlOIgAAAAAAAAAAAY1M6mAwAAAAAAAAAAMK9TnZzq2nSMtU71phOMR+cgAAAAAAAAAAAYlOIgAAAAAAAAAAAYlGPFAAAAWOvIsRObjgAwi73+fjt6+NA+JwEAAIB5dVd6AceKLSHj0ugcBAAAAAAAAAAAg1IcBAAAAAAAAAAAg1IcBAAAAAAAAAAAg9rZdAAAAAAAAAAAAObVXTnVtekYa/UCMi6NzkEAAAAAAAAAADAoxUEAAAAAAAAAADAox4oBAAAAAAAAAAyuk3RvOsV6C4i4ODoHAQAAAAAAAADAoBQHAQAAAAAAAADAoBQHAQAAAAAAAADAoHY2HQAAAAAAAAAAgHmdSuVUatMx1lpCxqXROQgAAAAAAAAAAAalOAgAAAAAAAAAAAblWDEAAIAD5MixE5uOADCEvf4+PXr40D4nAQAAgAvTXene/iO7lpBxaXQOAgAAAAAAAACAQSkOAgAAAAAAAACAQSkOAgAAAAAAAACAQe1sOgAAAAAAAAAAAPM61ZVTXZuOsdYSMi6NzkEAAAAAAAAAADAoxUEAAAAAAAAAADAoxUEAAAAAAAAAADConU0HAAAAAAAAAABgXt27Y9stIePS6BwEAAAAAAAAAACDUhwEAAAAAAAAAACD2khxUFW9sap6ZXzdBTzzwqp6Z1U9UFWfmz7fWVUvvAyRAQAAAAAAAACWqyu9gJGuTf+XGs5lLw6qqr+S5LaLWF9VdWeSe5LclOTaJFdOnzcluaeq7qwq/zoAAAAAAAAAAGDFZS0OqqorkvzrJDtJPnqBj/1okldO8w8keVmSZ0+fH5i+f2WSN+xfUgAAAAAAAAAAWL7L3Tno7yd5VpJ7k7x93eKquj7Ja6fL9yd5bne/o7uPd/c7kjxv+j5Jbq+qp82QGQAAAAAAAAAAFumyFQdV1ZPzaHef70vyyAU8dlt2uwwlya3d/fDqze7+TJJbp8udJK/Zh6gAAAAAAAAAAEPprsUM9tfl7Bz0L5M8PslPd/evr1tcVZXkW6bLe7v7fWdbN33/e9PlS6bnAAAAAAAAAADgwLssxUFV9W1JXpzkj5N8/wU+9tQk107zY2vWnr5/XZKnXGw+AAAAAAAAAAAY0c76JZemqr40yU9Ml7d398cu8NEbV+b3rlm7ev/GJPdf4DsAAAAAAAAAAIZ3qiunFnBk1xIyLs3sxUFJ3pjkmiT/LcnbL+K5J6/MH1iz9uQ5nrsgVXXdmiXXXOyeAAAAAAAAAACwabMWB1XV85LckuTzSV7d3X0Rjz9hZf7pNWsfWpk//iLecdrJ9UsAAAAAAAAAAGBZrphr46q6MsnbklSSf9bdH7zILa5amT+yZu3nVuZXX+R7AAAAAAAAAABgSHN2DvqBJDcm+b9JfmQPz392ZX7lmrWPXZk/vId3rTuK7Jokx/ewLwAAAAAAAADAxnXvjm23hIxLM0txUFXdkOQfTZe3dvdD51t/Dp9ama87KuxxK/N1R5D9Gd39wPnuV9XFbgkAAAAAAAAAABs3V+eg27Lb7ef/JPniqvr2s6z5yyvzv1lV10zzX56KiVYLdq5b877Vzj8nLzYsAAAAAAAAAACMaK7ioNPHfP2FJD9/Aev/8cr8qUkeSnJi5bsb1jy/ev93L+B9AAAAAAAAAAAHxu6xYtt/cpJjxfbfFZsOcB73J3lwmh9es/b50+cfJvnQXIEAAAAAAAAAAGBJZikO6u6bu7vON5L8yMojX79y70PTHp3k7un+DVX1nLO9a/r+dOegu6fnAAAAAAAAAADgwNvmzkFJ8uYkn5/mb6mqq1dvTtdvmS4/P60HAAAAAAAAAACS7Gw6wPl09+9X1ZuSvC7JM5P8RlX9WJL7kjwtye1JnjEt//Hu/oPNJAUAALi8jhw7sekIAOzBXn9/Hz18aJ+TAAAAcNB0KrsHPW23zvZnXJqtLg6a/GCSJyV5eXYLgd5xljVvT/JDlzMUAAAAAAAAAABsu20/Vizdfaq7X5Hkm5LcneTBJI9Mn3cneVF339LdpzYYEwAAAAAAAAAAts7GOgd19x1J7riI9e9O8u658gAAAAAAAAAAjKqnse2WkHFptr5zEAAAAAAAAAAAsDeKgwAAAAAAAAAAYFCKgwAAAAAAAAAAYFA7mw4AAAAAAAAAAMC8uivdtekYay0h49LoHAQAAAAAAAAAAINSHAQAAAAAAAAAAINyrBgAAAAAAAAAwOh6GttuCRkXRucgAAAAAAAAAAAYlOIgAAAAAAAAAAAOhKp6UlW9uKpeX1X3VNXHq6qncdc+vueOlX3Xja/br/eejWPFAAAAAAAAAAA4KD6y6QCXm+IgAAAAAAAAAIDRdaW7Np1ivcub8WSS303yjTO/5+lr7t8/58sVBwEAAAAAAAAAcFC8PsnxJMe7+yNV9ZTMXJzT3f9rzv3XURwEAAAAAAAAAMCB0N0/vOkMl5viIAAAAAAAAACAwXXvjm23hIxLc8WmAwAAAAAAAAAAAPNQHAQAAAAAAAAAADOpqvdU1Seq6pGq+mhV/XpVva6qvuxyvN+xYgAAAAAAAAAAbKNrquq8C7r7gcuU5VK8YGX+FUkOT+P2qrq5u++e8+WKgwAAAAAAAAAABtdd6T5/oc02OCPj8Qt4ZJt/qA8m+aUkv5XkwSRflOQvJfnOJN+Y5EuT/GJVfXN33zNXCMVBAAAAAAAAAACwv97c3Xec5fvfTPJvq+pVSd6a5DFJfrKqru/uh+cIojgIAAAAAAAAAIBt9KwkH950iL3o7k+uuX9nVT0zyS1J/nySlyb52TmyKA4CAAAAAAAAAGAbfbi7H9h0iBndmd3ioCQ5HMVBAAAAAAAAAADsSdfu2HZLyLh/TqzMr53rJVfMtTEAAAAAAAAAAHBOl6USSnEQAAAAAAAAAABcfodW5g/O9RLHigEAAGzQkWMn1i8C4MDb69+Lo4cPrV8EAADAgdC9O7bdEjLuo1etzI/N9RKdgwAAAAAAAAAA4AJV1c1V1dO44yz3n15V16/Z41VJXjFdfjjJu/Y/6S6dgwAAAAAAAAAAOBCq6nlJVgt3nrgyv76qbl5d39137eE1fz3JT1bVe5Pck+SDST6R3TqdG5J8V5K/Na39QpJXdfdDe3jPBVEcBAAAAAAAAADAQXFLku85x73nTmPVXXt8z2OSvGAa5/KJJK/o7qN7fMcFURwEAAAAAAAAADC6nsa2W0LG9d6d3SPDvjbJM5J8ZZIvT1JJ/jjJ/0jyq0nu6u4/nTuM4iAAAAAAAAAAAA6E7r45yc2XuMddOU9Hoe7+aJJ/M42Nu2LTAQAAAAAAAAAAgHnoHAQAAAAAAAAAMLjuSndtOsZaS8i4NDoHAQAAAAAAAADAoBQHAQAAAAAAAADAoBQHAQAAAAAAAADAoHY2HQAAAAAAAAAAgMugNx2ATdA5CAAAAAAAAAAABqU4CAAAAAAAAAAABuVYMQAAAAAAAACAwXVXumvTMdZaQsal0TkIAAAAAAAAAAAGpTgIAAAAAAAAAAAGpTgIAAAAAAAAAAAGtbPpAAAAAAAAAAAAzKynse2WkHFhdA4CAAAAAAAAAIBBKQ4CAAAAAAAAAIBBOVYMAAAAAAAAAGB4NY1tt4SMy6JzEAAAAAAAAAAADEpxEAAAAAAAAAAADEpxEAAAAAAAAAAADGpn0wEAAAAAAAAAAJhZT2PbLSHjwugcBAAAAAAAAAAAg1IcBAAAAAAAAAAAg3KsGAAAAAAAAADA6BwrdmDpHAQAAAAAAAAAAINSHAQAAAAAAAAAAINSHAQAAAAAAAAAAIPa2XQAAAAAAAAAAABm1km6Np1ivd50gPHoHAQAAAAAAAAAAINSHAQAAAAAAAAAAINSHAQAAAAAAAAAAIPa2XQAAAAAAAAAAADm1b07tt0SMi6NzkEAAAAAAAAAADAoxUEAAAAAAAAAADAox4oBAAAAAAAAAIyup7HtlpBxYXQOAgAAAAAAAACAQSkOAgAAAAAAAACAQSkOAgAAAAAAAACAQe1sOgAAAAAAAAAAADPr2h3bbgkZF0bnIAAAAAAAAAAAGJTiIAAAAAAAAAAAGJRjxQAAAAAAAAAABldJqjedYj2Hiu0/nYMAAAAAAAAAAGBQioMAAAAAAAAAAGBQioMAAAAAAAAAAGBQO5sOAAAAAAAAAADAzHoa224JGRdG5yAAAAAAAAAAABiU4iAAAAAAAAAAABiUY8UAAAAAAAAAAEbXtTu23RIyLozOQQAAAAAAAAAAMCjFQQAAAAAAAAAAMCjFQQAAAAAAAAAAMKidTQcAAAAAAAAAAGBmPY1tt4SMC6NzEAAAAAAAAAAADEpxEAAAAAAAAAAADMqxYgAAAAAAAAAAo3Os2IGlcxAAAAAAAAAAAAxKcRAAAAAAAAAAAAxKcRAAAAAAAAAAAAxqZ9MBAAAAAAAAAACYWU9j2y0h48LoHAQAAAAAAAAAAINSHAQAAAAAAAAAAINyrBgAAAAAAAAAwOi6dse2W0LGhdE5CAAAAAAAAAAABqU4CAAAAAAAAAAABqU4CAAAAAAAAAAABrWz6QAAAAAjOHLsxKYjAMCfsde/T0cPH9rnJAAAAGxa9e7YdkvIuDQ6BwEAAAAAAAAAwKAUBwEAAAAAAAAAwKAcKwYAAAAAAAAAMLqexrZbQsaF0TkIAAAAAAAAAAAGpTgIAAAAAAAAAAAGpTgIAAAAAAAAAAAGpTgIAAAAAAAAAAAGpTgIAAAAAAAAAAAGpTgIAAAAAAAAAAAGpTgIAAAAAAAAAAAGtbPpAAAAAAAAAAAAzKuSVG86xXq16QAD0jkIAAAAAAAAAAAGpTgIAAAAAAAAAAAG5VgxAACAfXD08KE9PXfk2Il9TgIAj9rr3ycAAAAG1LU7tt0SMi6MzkEAAAAAAAAAADAoxUEAAAAAAAAAADAoxUEAAAAAAAAAADConU0HAAAAAAAAAABgZj2NbbeEjAujcxAAAAAAAAAAAAxKcRAAAAAAAAAAAAzKsWIAAAAAAAAAAKNzrNiBpXMQAAAAAAAAAAAMSnEQAAAAAAAAAAAMSnEQAAAAAAAAAAAMameujavqS5K8KMmzkjwzybVJviLJ1Uk+meREkncneXt3f+IC9nthklcmefa0z8eS/FaSt3X3r87xMwAAAAAAAAAAjKB6d2y7JWRcmtmKg7JbxPPz57j3FUkOT+P7q+q7uvvXzrawqirJW7NbGLTq2iQ3Jbmpqt6W5NXd7Z8IAAAAAAAAAABM5iwOSpKTSd6b5Len+R9l9yiz65J8a5KXJnlikqNV9azu/p9n2eNH82hh0AeSvDHJfUmeluS1SZ4x3f9Ykh+a7ScBAAAAAAAAAICFmbM46L3d/VXnuf8LVfWSJO9KcmWSH07yd1YXVNX12S0ASpL3J3l+dz88XR+vqqNJjmX32LLbq+qnuvu+/fwhAAAAAAAAAAAWr6ex7ZaQcWGumGvj7v7CBaz5pST3TpfPP8uS2/JoAdOtK4VBp5//TJJbp8udJK/ZW1oAAAAAAAAAABjPbMVBF+Gh6fOq1S+rqpJ8y3R5b3e/72wPT9//3nT5kuk5AAAAAAAAAAA48DZaHFRVNyb5q9PlvWfcfmqSa6f5sTVbnb5/XZKn7Es4AAAAAAAAAABYuJ31S/ZXVX1xdot+vjnJa5M8Zrr1E2csvXFlfmbh0JlW79+Y5P5LyQgAAAAAAAAAMJSexrZbQsaFuSzFQVV1c5KfOs+SNyX52TO+e/LK/IE1rzh5jucuSFVdt2bJNRe7JwAAAAAAAAAAbNpl7xx0ht9J8uru/s2z3HvCyvzTa/Z5aGX++D3kOLl+CQAAAAAAAAAALMvlKg76pSTvn+ZXJ3lakm9LclOSn62q13T3r5zxzFUr80fW7P+5lfnVlxIUAAAAAAAAAGA01btj2y0h49JcluKg7v5kkk+ufHU8yTuq6ruT/HSSu6vqFd1918qaz67Mr1zziseuzB/eQ8R1R5Fdk93MAAAAAAAAAACwGBs9Vqy7f6aqXpzdLkL/vKru7u4/mW5/amXpuqPCHrcyX3cE2dlyPHC++1V1sVsCAAAAAAAAAMDGXbHpAEnunj4fl+Rvr3y/WrBz3Zo9Vjv/nNyPUAAAAAAAAAAAsHQb7Rw0+djK/KtX5idW5jes2WP1/u9eciIAAAAAAAAAgJF07Y5tt4SMC7MNnYOuXZmvHgl2f5IHp/nhNXs8f/r8wyQf2p9YAAAAAAAAAACwbNtQHPR3V+YfPD3p7s6jR47dUFXPOdvD0/enOwfdPT0HAAAAAAAAAAAH3mzFQVV1c1VdtWbNbUleNF1+KMl/PWPJm5N8fpq/paquPuP5q5O8Zbr8/LQeAAAAAAAAAIAz9QIG+25nxr3vSPJPq+oXs1v0c192jw17QpKnJ/nOJM+d1j6S5Hu7+/OrG3T371fVm5K8Lskzk/xGVf3YtNfTktye5BnT8h/v7j+Y8ecBAAAAAAAAAIBFmbM4KEn+XJLvnca5PJDk5d39H89x/weTPCnJy7NbCPSOs6x5e5IfuoScAAAAAAAAAAAwnDmLg74hyQuSfH2SG5N8ZZIvT/LZJB9J8jtJfiXJL3T3Z861SXefSvKKqQPRK5M8K8kTk3w8yfEkd3b3PTP+HAAAAAAAAAAAsEizFQd1933ZPf7rzn3a791J3r0fewEAAAAAAAAAHCTVu2PbLSHj0lyx6QAAAAAAAAAAAMA8FAcBAAAAAAAAAMCgZjtWDAAAAAAAAACALdHT2HZLyLgwOgcBAAAAAAAAAMCgFAcBAAAAAAAAAMCgFAcBAAAAAAAAAMCgdjYdAAAAAAAAAACAmXVSvekQF2AJGRdG5yAAAAAAAAAAABiU4iAAAAAAAAAAABiU4iAAAAAAAAAAABjUzqYDAAAAAAAAAAAws57GtltCxoXROQgAAAAAAAAAAAalOAgAAAAAAAAAAAblWDEAAAAAAAAAgNE5VuzA0jkIAAAAAAAAAAAGpTgIAAAAAAAAAAAGpTgIAAAAAAAAAAAGtbPpAAAAAAAAAAAAzKt6d2y7JWRcGp2DAAAAAAAAAABgUIqDAAAAAAAAAABgUIqDAAAAAAAAAABgUIqDAAAAAAAAAABgUIqDAAAAAAAAAABgUIqDAAAAAAAAAABgUDubDgAAAAAAAAAAwMx6GttuCRkXRucgAAAAAAAAAAAYlOIgAAAAAAAAAAAYlGPFAAAAAAAAAAAGV707tt0SMi6NzkEAAAAAAAAAADAoxUEAAAAAAAAAADAoxUEAAAAAAAAAADConU0HAAAAAAAAAADgMuhNB2ATdA4CAAAAAAAAAIBBKQ4CAAAAAAAAAIBBOVYMAAAAAAAAAGB0nWUcK7aEjAujcxAAAAAAAAAAAAxKcRAAAAAAAAAAAAxKcRAAAAAAAAAAAAxqZ9MBAAAADrKjhw/t6bkjx07scxIAttle/14AAADAadW7Y9stIePS6BwEAAAAAAAAAACDUhwEAAAAAAAAAACDcqwYAAAAAAAAAMDoehrbbgkZF0bnIAAAAAAAAAAAGJTiIAAAAAAAAAAAGJTiIAAAAAAAAAAAGNTOpgMAAAAAAAAAADCv6t2x7ZaQcWl0DgIAAAAAAAAAgEEpDgIAAAAAAAAA4ECoqidV1Yur6vVVdU9Vfbyqehp3zfTOb6+qX6uqP6qqz1bVh6rqZ6rqOXO870yOFQMAAAAAAAAAGF1PY9vNn/Ejs79hUlVXJfl3SV58xq2vnsZ3VNUd3f2GOXPoHAQAAAAAAAAAwEF0Msl/mHH/t+fRwqD3JnlJkmcneUWS+7Jbt/P6qrplxgw6BwEAAAAAAAAAcGC8PsnxJMe7+yNV9ZQk9+/3S6rqcJLvmC5/OclN3f2F6fp4VR1N8ttJvirJG6vq33f3J/c7R6JzEAAAAAAAAAAAB0R3/3B3/0p3z3282Gunzy8k+XsrhUGnc3w8ye3T5Zdlt5vQLBQHAQAAAAAAAACMrhc0Fq6qHp/kG6bL93T3A+dY+s4kfzrNXzpXHsVBAAAAAAAAAACwf56d5LHT/Ni5FnX3I0ned/qZqvqiOcIoDgIAAAAAAAAAgP1z48r83jVrT9/fSfI1c4TZmWNTAAAAAAAAAAC4RNdU1XkXnOfIrk168sp8Xb6TZzx3Yr/DKA4CAAAAAAAAABhc9e7YdmdkPH4hj8yT5JI8YWX+6TVrH1qZP36GLI4VAwAAAAAAAACAfXTVyvyRNWs/tzK/eoYsOgcBAAAAAAAAALCVnpXkw5sOsQefXZlfuWbtY1fmD8+QRXEQAAAAAAAAAMDwehrb7v/P+OHufmBDSS7Fp1bm644Ke9zKfN0RZHviWDEAAAAAAAAAANg/qwVN161Z++SV+ckZsigOAgAAAAAAAACAfXRiZX7DmrWn738+yf+eI4ziIAAAAAAAAAAA2D/HkzwyzQ+fa1FVXZnkOaef6e5HzrX2UigOAgAAAAAAAAA4CHoBYwDd/akk/2m6fEFVnetosZcm+ZJp/q658igOAgAAAAAAAACAC1RVN1dVT+OOcyx70/S5k+RfVNVjztjjiUl+bLr8ZJKfnCXsFAAAAAAAAAAAAIZXVc9Lcv3KV09cmV9fVTevru/uu/bynu7+z1X1jiTfnuRIkvdU1ZuTPJjk6Ul+MMlXTctf191/spf3XAjFQQAAAAt09PChPT135NiJfU4CwMXY6+9vAAAAuFTVu2PbXYaMtyT5nnPce+40Vt11Ce96eXaPDXtRkq+fxqpTSd7Q3XdewjvWcqwYAAAAAAAAAADss+5+uLu/Kcl3JnlPko8meSTJySQ/l+R53X3H3Dl0DgIAD4vzrQAAIABJREFUAAAAAAAA4EDo7puT3HyJe9yVi+go1N0/l91ioI3QOQgAAAAAAAAAAAalcxAAAAAAAAAAwOh6GttuCRkXRucgAAAAAAAAAAAYlOIgAAAAAAAAAAAYlGPFAAAAAAAAAAAGV707tt0SMi6NzkEAAAAAAAAAADAoxUEAAAAAAAAAADAoxUEAAAAAAAAAADConU0HAAAAAAAAAABgZj2NbbeEjAujcxAAAAAAAAAAAAxKcRAAAAAAAAAAAAzKsWIAAAAAAAAAAKNzrNiBpXMQAAAAAAAAAAAMSnEQAAAAAAAAAAAMSnEQAAAAAAAAAAAMamfTAQAAAAAAAAAAmFdNY9stIePS6BwEAAAAAAAAAACDUhwEAAAAAAAAAACDcqwYAAAAAAAAAMDoehrbbgkZF0bnIAAAAAAAAAAAGJTiIAAAAAAAAAAAGJTiIAAAAAAAAAAAGNTOpgMAAAAAAAAAADCzTqo3HeICLCHjwugcBAAAAAAAAAAAg1IcBAAAAAAAAAAAg3KsGAAAwAFy9PChPT135NiJfU4CsGx7/X0KAAAAG9NZxpFdS8i4MDoHAQAAAAAAAADAoBQHAQAAAAAAAADAoBQHAQAAAAAAAADAoHY2HQAAAAAAAAAAgMugNx2ATdA5CAAAAAAAAAAABqU4CAAAAAAAAAAABqU4CAAAAAAAAAAABrWz6QAAAAAAAAAAAMyrendsuyVkXBqdgwAAAAAAAAAAYFCKgwAAAAAAAAAAYFCOFQMAAAAAAAAAGF1PY9stIePC6BwEAAAAAAAAAACDUhwEAAAAAAAAAACDUhwEAAAAAAAAAACD2tl0AAAAALbf0cOH9vTckWMn9jkJwP7a6+83AAAAWJrq3bHtlpBxaXQOAgAAAAAAAACAQSkOAgAAAAAAAACAQTlWDAAAAAAAAABgdD2NbbeEjAujcxAAAAAAAAAAAAxKcRAAAAAAAAAAAAxKcRAAAAAAAAAAAAxqZ9MBAAAAAAAAAACYV/Xu2HZLyLg0OgcBAAAAAAAAAMCgFAcBAAAAAAAAAMCgZi0Oqqq/VlU/UFX3VNXJqvpcVX26qn6/qu6qqr9xkfu9sKreWVUPTHs9MF2/cK6fAQAAAAAAAABg8XpBg321M9fGVXUsyfPPcuvKJF8zje+pqp9Jckt3P3KevSrJW5O88oxb1ya5KclNVfW2JK/ubv9MAAAAAAAAAAAg83YOunb6fDDJTyT51iTPTvK1Sf5Bkj+c7n93krvW7PWjebQw6ANJXjbt9bLpOtP9N+xDbgAAAAAAAAAAGMJsnYOS3JvkB5L8Ynd/4Yx775s6Bv1Gkr+Y5GVV9a+6+7+cuUlVXZ/ktdPl+5M8v7sfnq6PV9XRJMeSPDPJ7VX1U9193ww/DwAAAAAAAAAALMpsxUHd/eI19z9eVf8wyS9PX31rkj9THJTktjya89aVwqDT+3ymqm5N8t+nda9JcuulZAcAAGB/HD18aM/PHjl2Yh+TAKO7lN83AAAAcCD0NLbdEjIuzJzHil2IX1+ZP+3Mm1VVSb5lury3u993tk2m739vunzJ9BwAAAAAAAAAABxomy4OunJlfuos95+a5NppfmzNXqfvX5fkKZcWCwAAAAAAAAAAlm+2Y8Uu0OGV+b1nuX/jmvs5x/0bk9y/11AAAAAAAAAAACOp3h3bbgkZl2ZjxUFVdUWS16189QtnWfbklfkDa7Y8eY7nLiTLdWuWXHMx+wEAAAAAAAAAwDbYZOeg25I8e5q/q7vff5Y1T1iZf3rNfg+tzB9/kVlOrl8CAAAAAAAAAADLcsUmXlpVh5P8k+nyo0m+7xxLr1qZP7Jm28+tzK/eYzQAAP4fe/cfbN1V1gf8+4RLQhIUMyLNaKBaYgmoMAihoISgKLUEMaJVsVUDSUGdQRFUsDgVQ4vDrxEG8RdGIkprUYGAoAULRmT4ER2m2iEQAmgTMWj4oSTkh+E+/ePsd97D67333Pu+99xz9r6fz8yas/bZa6/1XIas3HPy3GcBAAAAAAAwGQdeOaiqvirJ64a1b0vyXd398W2G3zrXP3nB1KfM9W/ZY1iLjiE7M8lVe5wTAAAAAAAAAGB99KoDYBUONDmoqr4iyVuSnJHkc0me0N1X7vDIZ+b6i44KO32uv+gIss/T3dfvdL+q9jIdAAAAAAAAAACshQM7VqyqvjTJHyX50sxy0Z7U3a9b8Nh80s5ZC8bOV/+5bu8RAgAAAAAAAADAtBxI5aCqunuStyb5V8NbT+3uV+3i0ffP9c9ZMHb+/tV7CA8AAAAAAAAAYNKqO9Xrf67YGGIcm6VXDqqquyX5X0nuN7z1rO5++S4f/2iSjw398xeMfcTw+jdJ/movMQIAAAAAAAAAwBQtNTmoqk5L8qYkXzu89d+6+/m7fb67O8kVw+U5VfXQbdZ5aI5WDrpieA4AAAAAAAAAAA61pSUHVdXJSV6X5OuHt17a3T99HFO9JMkdQ/9lVXXqMeucmuRlw+Udw3gAAAAAAAAAADj0NpY49/9I8uih/7Ykl1XVV+8w/vbuvubYN7v7mqp6UZJnJXlwkndW1fOTfDjJvZM8M8kDh+Ev7O4P7dcPAAAAwOq84fz7LR60jx535fsPdD2YsoP+5xcAAADYhR7auhtDjCOzzOSgx8/1vzHJXywY/9dJvnybe89Oco8kT8osEei3txhzWZLjqUwEAAAAAAAAAACTtLRjxfZTd29298VJLkhyRZKPJbl9eL0iyWO6+5Lu3lxhmAAAAAAAAAAAsFaWVjmou2sJc745yZv3e14AAAAAAAAAgCmrnrV1N4YYx2YUlYMAAAAAAAAAAIC9kxwEAAAAAAAAAAATJTkIAAAAAAAAAAAmamPVAQAAAAAAAAAAsGQ9tHU3hhhHRuUgAAAAAAAAAACYKMlBAAAAAAAAAAAwUZKDAAAAAAAAAABgojZWHQAAAAAAAAAAAMtVPWvrbgwxjo3KQQAAAAAAAAAAMFGSgwAAAAAAAAAAYKIcKwYAAABJ3nD+/Q50vcdd+f4DXY/D7aD//w0AAACsoR7auhtDjCOjchAAAAAAAAAAAEyU5CAAAAAAAAAAAJgoyUEAAAAAAAAAADBRG6sOAAAAAAAAAACA5aqetXU3hhjHRuUgAAAAAAAAAACYKMlBAAAAAAAAAAAwUY4VAwAAAAAAAACYuh7auhtDjCOjchAAAAAAAAAAAEyU5CAAAAAAAAAAAJgoyUEAAAAAAAAAADBRG6sOAAAAAAAAAACA5atedQSsgspBAAAAAAAAAAAwUZKDAAAAAAAAAABgohwrBgAAACvwhvPvt+oQ1tbjrnz/cT3nf1MAAACAHXTP2robQ4wjo3IQAAAAAAAAAABMlOQgAAAAAAAAAACYKMlBAAAAAAAAAAAwURurDgAAAAAAAAAAgOWqnrV1N4YYx0blIAAAAAAAAAAAmCjJQQAAAAAAAAAAMFGOFQMAAAAAAAAAmLoe2robQ4wjo3IQAAAAAAAAAABMlOQgAAAAAAAAAACYKMlBAAAAAAAAAAAwURurDgAAAAAAAAAAgOWqzVlbd2OIcWxUDgIAAAAAAAAAgImSHAQAAAAAAAAAABPlWDEAAABgrbzh/PutOgQAAACA6emhrbsxxDgyKgcBAAAAAAAAAMBESQ4CAAAAAAAAAICJkhwEAAAAAAAAAAATtbHqAAAAAAAAAAAAWK7qWVt3Y4hxbFQOAgAAAAAAAACAiZIcBAAAAAAAAAAAE+VYMQAAAAAAAACAqeuetXU3hhhHRuUgAAAAAAAAAACYKMlBAAAAAAAAAAAwUZKDAAAAAAAAAABgojZWHQAAAAAAAAAAAMtVPWvrbgwxjo3KQQAAAAAAAAAAMFGSgwAAAAAAAAAAYKIkBwEAAAAAAAAAwERtrDoAAAAAAAAAAAAOQK86AFZB5SAAAAAAAAAAAJgoyUEAAAAAAAAAADBRjhUDAAAAAAAAAJi46llbd2OIcWxUDgIAAAAAAAAAgImSHAQAAAAAAAAAABMlOQgAAAAAAAAAACZqY9UBAAAAAAAAAACwZN2ztu7GEOPIqBwEAAAAAAAAAAATJTkIAAAAAAAAAAAmyrFiAAAAAAAAAAATVz1r624MMY6NykEAAAAAAAAAADBRkoMAAAAAAAAAAGCiJAcBAAAAAAAAAMBEbaw6AAAAAAAAAAAAlqyHtu7GEOPIqBwEAAAAAAAAAAATJTkIAAAAAAAAAAAmyrFiAAAAAAAAAAATVz1r624MMY6NykEAAAAAAAAAADBRkoMAAAAAAAAAAGCiJAcBAAAAAAAAAMBEbaw6AAAAAAAAAAAAlmyzZ23djSHGkVE5CAAAAAAAAAAAJkpyEAAAAAAAAAAATJRjxQAAAAAAAAAApq6Htu7GEOPIqBwEAAAAAAAAAAATJTkIAAAAAAAAAAAmSnIQAAAAAAAAAABM1MaqAwAAAAAAAAAAYLkqSfWqo1isVh3ABKkcBAAAAAAAAAAAEyU5CAAAAAAAAAAAJsqxYgAAAAAAAAAAU9dJegTnio0gxLFROQgAAAAAAAAAACZKchAAAAAAAAAAAEyU5CAAAAAAAAAAAA6dqrpXVb2oqq6uqpur6pNV9d6q+vGqOu0E535OVfUu2yP36Ufa0sYyJwcAAAAAAAAAYA10Ur3qIHbhgGKsqguSvDrJ3ebePi3JuUO7pKoe090fOZiIlkdyEAAAAAAAAAAAh0ZVPSDJazJLBropyc8leXuSU5N8T5L/lOQ+Sd5UVed2900nuOTXLLj/0ROcf0eSgwAAAAAAAAAAOExeklli0B1JHt3d75q797aq+lCSFyQ5J8nTk1x6Iot19/89kedP1EmrXBwAAAAAAAAAgAPQI2pLVFXnJnnkcHnZMYlBR7w4ydVD/2lVdeflRrVckoMAAAAAAAAAADgsLpzrv3KrAd29meRVw+UZOZpMNEqSgwAAAAAAAAAAOCzOG15vTvLnO4y7cq7/8OWFs3ySgwAAAAAAAAAAOCzuO7xe29137DDuA1s8c1yq6q1V9Ymqur2q/q6q/riqnlVVZ5zIvLu1cRCLAAAAAAAAAACwOtWd6l51GAsdE+OZVbXj+O6+ftdzV90lyd2Hyx2f6+5PVdXNSU5Pcs/drrGNb5rrf0mS84f2zKq6qLuvOMH5dyQ5CAAAAAAAAACAdXTVLsbsnD30+b5grn/TLsYfSQ666x7WmPeXSV6f5L1JPpbkzknuk+Q/JHl0ki9K8ntV9a3d/QfHucZCkoMAAAAAAAAAADgM7jLXv30X428bXk89jrVe0t3P2eL99yR5VVU9JckvJ7lTkl+rqrO7+5bjWGchyUEAAAAAAAAAAKyjc5PcsI/z3TrXP3kX408ZXvectNPdn15w/1eq6sFJLknypUken+TVe11nNyQHAQAAAAAAAABM3ebQ1t3nx3hDd1+/j7N/Zq6/m6PCTh9ed3ME2fH4lcySg5Lk/CwpOeikZUwKAAAAAAAAAADrpLtvTXLjcHnWTmOr6owcTQ66bkkhvX+u/2VLWkNyEAAAAAAAAAAAh8bVw+vZVbXTiVvnbPHMfqslzft5HCsGAAAAAAAAADBx1Z3qXnUYCx1AjH+a5LzMqgI9KMl7thl3/lz/nUuK5X5z/Y8taQ2VgwAAAAAAAAAAODReP9d/4lYDquqkJN8/XH46yduXFMtT5vpXLmkNyUEAAAAAAAAAABwO3f3eJO8YLi+uqodtMewZSe479F/a3f80f7OqLqqqHtpzjn24qr6mqs7eKY6qekqSi4fLG5K8bg8/xp44VgwAAAAAAAAAgMPkRzM7KuzUJG+pqudlVh3o1CTfk+TJw7hrkrz4OOZ/UJJfq6q3J/mDJH+Z5BOZ5emck+Q/JvnmYeznkjylu28+vh9lMclBAAAAAAAAAABT10NbdwcQY3e/r6q+O8lvJfnCJM/bYtg1SS7o7s8c5zJ3SvJNQ9vOJ5Jc3N1vOM41dkVyEAAAAAAAAAAAh0p3v7Gq7p9ZFaELkpyV5PYk1yb5nSS/0N2fPc7p35zZkWEPS/LAJP8iyRcnqSSfTPJ/kvxhksu7+x9P5OfYDclBAAAAAAAAAAAcOt3910mePrS9PHd5kst3uP93SX59aCsnOQgAAAAAAAAAYPI6aeeKHUYnrToAAAAAAAAAAABgOSQHAQAAAAAAAADAREkOAgAAAAAAAACAidpYdQAAAAAAAAAAACxX9aytuzHEODYqBwEAAAAAAAAAwERJDgIAAAAAAAAAgIlyrBgAAAAAAAAAwNR1z9q6G0OMI6NyEAAAAAAAAAAATJTkIAAAAAAAAAAAmCjJQQAAAAAAAAAAMFFLTQ6qqntU1WOr6tKq+oOqurGqemiXH8d831JVr62q66vqtuH1tVX1LUsIHwAAAAAAAABgEmpzPI39tbHk+T++H5NUVSX55SRPPubWlyX59iTfXlW/muQHu7v3Y00AAAAAAAAAABi7gzxW7LokbznOZ/9rjiYGvS/JE5I8ZHh93/D+k5M890QCBAAAAAAAAACAKVl25aBLk1yV5Kru/nhVfXmSj+5lgqo6O8lPDpd/luQR3X3LcH1VVb0hyZVJHpzkmVX1yu7+8H4EDwAAAAAAAAAwCd2ztu7GEOPILLVyUHf/THf/fnefyPFiP5ajSUxPnUsMOrLGZ5M8dbjcSPK0E1gLAAAAAAAAAAAm4yCPFduzqqok3zZcfqC7373VuOH9Dw6XFw7PAQAAAAAAAADAobbWyUFJviLJlw39KxeMPXL/rCRfvqyAAAAAAAAAAABgLDYWD1mp+871P7Bg7Pz9+yb56P6HAwAAAAAAAAAwQj20dTeGGEdm3ZOD7jnXv37B2Ou2eW6hqjprwZAz9zIfAAAAAAAAAACsg3VPDvqCuf5NC8bePNe/6x7XuW7xEAAAAAAAAAAAGJd1Tw66y1z/9gVjb5vrn7qEWAAAAAAAAAAARqm6U73+Z3aNIcaxWffkoFvn+icvGHvKXP+WPa6z6BiyM5Nctcc5AQAAAAAAAABgpdY9Oegzc/1FR4WdPtdfdATZ5+nu63e6X1V7mQ4AAAAAAAAAANbCSasOYIH5pJ2zFoydr/5z3RJiAQAAAAAAAACAUVn3ykHvn+ufs2Ds/P2rlxALAAAAAAAAAMA4dc/auhtDjCOz7pWDPprkY0P//AVjHzG8/k2Sv1pWQAAAAAAAAAAAMBZrnRzU3Z3kiuHynKp66FbjhvePVA66YngOAAAAAAAAAAAOtbVODhq8JMkdQ/9lVXXq/M3h+mXD5R3DeAAAAAAAAAAAjugkmyNoysHsu41lTl5VD09y9txbd5/rn11VF82P7+7Lj52ju6+pqhcleVaSByd5Z1U9P8mHk9w7yTOTPHAY/sLu/tC+/QAAAAAAAAAAADBiS00OSnJJkh/Y5t7XD23e5duMfXaSeyR5UmaJQL+9xZjLkvz03kMEAAAAAAAAAIBpGsOxYunuze6+OMkFSa5I8rEktw+vVyR5THdf0t2bKwwTAAAAAAAAAADWylIrB3X3RUku2sf53pzkzfs1HwAAAAAAAADAYVDdqe5Vh7HQGGIcm1FUDgIAAAAAAAAAAPZOchAAAAAAAAAAAEyU5CAAAAAAAAAAAJiojVUHAAAAAAAAAADAknWS7lVHsdgIQhwblYMAAAAAAAAAAGCiJAcBAAAAAAAAAMBEOVYMAAAAAAAAAGDqukdyrNgIYhwZlYMAAAAAAAAAAGCiJAcBAAAAAAAAAMBESQ4CAAAAAAAAAICJ2lh1AAAAAAAAAAAALNnm0NbdGGIcGZWDAAAAAAAAAABgoiQHAQAAAAAAAADARDlWDAAAAAAAAABg4qo71b3qMBYaQ4xjo3IQAAAAAAAAAABMlOQgAAAAAAAAAACYKMlBAAAAAAAAAAAwURurDgAAAAAAAAAAgCXrnrV1N4YYR0blIAAAAAAAAAAAmCjJQQAAAAAAAAAAMFGOFQMAAAAAAAAAmLyRHCuWMcQ4LioHAQAAAAAAAADAREkOAgAAAAAAAACAiZIcBAAAAAAAAAAAE7Wx6gAAAAAAAAAAAFiy7llbd2OIcWRUDgIAAAAAAAAAgImSHAQAAAAAAAAAABPlWDEAAAAAAAAAgKnbHNq6G0OMI6NyEAAAAAAAAAAATJTkIAAAAAAAAAAAmCjJQQAAAAAAAAAAMFEbqw4AAAAAAAAAAIDlqu5U96rDWGgMMY6NykEAAAAAAAAAADBRkoMAAAAAAAAAAGCiHCsGAAAAAAAAADB13bO27sYQ48ioHAQAAAAAAAAAABMlOQgAAAAAAAAAACZKchAAAAAAAAAAAEzUxqoDAAAAAAAAAABgyTZ71tbdGGIcGZWDAAAAAAAAAABgoiQHAQAAAAAAAADARDlWDAAAAAAAAABg6jpJj+DIrhGEODYqBwEAAAAAAAAAwERJDgIAAAAAAAAAgImSHAQAAAAAAAAAABO1seoAAAAAAAAAAABYtk66Vx3ELowhxnFROQgAAAAAAAAAACZKchAAAAAAAAAAAEyU5CAAAAAAAAAAAJiojVUHAAAAAAAAAADAknXP2robQ4wjo3IQAAAAAAAAAABMlOQgAAAAAAAAAACYKMeKAQAAAAAAAABM3WbP2robQ4wjo3IQAAAAAAAAAABMlOQgAAAAAAAAAACYKMlBAAAAAAAAAAAwURurDgAAAAAAAAAAgCXrzVlbd2OIcWRUDgIAAAAAAAAAgImSHAQAAAAAAAAAABPlWDEAAAAAAAAAgKnrnrV1N4YYR0blIAAAAAAAAAAAmCjJQQAAAAAAAAAAMFGSgwAAAAAAAAAAYKI2Vh0AAAAAAAAAAABLttmztu7GEOPIqBwEAAAAAAAAAAATJTkIAAAAAAAAAAAmyrFiAAAAAAAAAABT1z1r624MMY6MykEAAAAAAAAAADBRkoMAAAAAAAAAAGCiJAcBAAAAAAAAAMBEbaw6AAAAAAAAAAAADkD3qiNgBVQOAgAAAAAAAACAiZIcBAAAAAAAAAAAE+VYMQAAAAAAAACAqesex7FiY4hxZFQOAgAAAAAAAACAiZIcBAAAAAAAAAAAEyU5CAAAAAAAAAAAJmpj1QEAAAAAAAAAALBkm5tJba46isU2RxDjyKgcBAAAAAAAAAAAEyU5CAAAAAAAAAAAJsqxYgAAAAAAAAAAU9c9a+tuDDGOjMpBAAAAAAAAAAAwUZKDAAAAAAAAAABgoiQHAQAAAAAAAADARG2sOgAAAAAAAAAAAJase9bW3RhiHBmVgwAAAAAAAAAAYKIkBwEAAAAAAAAAwERJDgIAAAAAAAAAgInaWHUAAAAAAAAAAAAs2WYn1auOYrHNEcQ4MioHAQAAAAAAAADAREkOAgAAAAAAAACAiXKsGAAAAAAAAADAxHVvpntz1WEsNIYYx0blIAAAAAAAAAAAmCjJQQAAAAAAAAAAMFGSgwAAAAAAAAAAYKI2Vh0AAAAAAAAAAABL1kk2e9VRLDaCEMdG5SAAAAAAAAAAAJgoyUEAAAAAAAAAADBRjhUDAAAAAAAAAJi67llbd2OIcWRUDgIAAAAAAAAAgImSHAQAAAAAAAAAABMlOQgAAAAAAAAAACZqY9UBAAAAAAAAAACwZJubSTZXHcVimyOIcWRUDgIAAAAAAAAAgImSHAQAAAAAAAAAABPlWDEAAAAAAAAAgKnrnrV1N4YYR0blIAAAAAAAAAAAmCjJQQAAAAAAAAAAMFGSgwAAAAAAAAAAYKI2Vh0AAAAAAAAAAADL1Zub6WyuOoyFenP9YxwblYMAAAAAAAAAAGCiJAcBAAAAAAAAAMBEOVYMAAAAAAAAAGDqumdt3Y0hxpFROQgAAAAAAAAAACZKchAAAAAAAAAAAEyU5CAAAAAAAAAAAJiojVUHAAAAAAAAAADAknUnm73qKBbrEcQ4MioHAQAAAAAAAADAREkOAgAAAAAAAACAiXKsGAAAAAAAAADA1HUn2Vx1FIs5Vmzfja5yUFXdq6peVFVXV9XNVfXJqnpvVf14VZ226vgAAAAAAAAAAGBdjKpyUFVdkOTVSe429/ZpSc4d2iVV9Zju/sgq4gMAAAAAAAAAgHUymspBVfWAJK/JLDHopiTPTvJ1SR6V5BXDsPskeVNV3XUlQQIAAAAAAAAAwBoZU+Wgl2RWJeiOJI/u7nfN3XtbVX0oyQuSnJPk6UkuPfgQAQAAAAAAAADWT292unrVYSzUvf4xjs0oKgdV1blJHjlcXnZMYtARL05y9dB/WlXd+SBiAwAAAAAAAACAdTWK5KAkF871X7nVgO7eTPKq4fKMHE0mAgAAAAAAAACAQ2ksx4qdN7zenOTPdxh35Vz/4UneurSIAAAAAAAAAADGojeTbK46isV6BDGOzFiSg+47vF7b3XfsMO4DWzyzUFWdtWDImbudCwAAAAAAAACA9VdV90ryI0kuSHKvJLcluTbJa5L8Ynd/dp/W+Z4kT0xy/8xOw7ohyTuSvLy7370fa+xk7ZODquouSe4+XF6/09ju/lRV3Zzk9CT33MMy1x1neAAAAAAAAAAAjExVXZDk1UnuNvf2aUnOHdolVfWY7v7ICaxxlyS/k+Sxx9z6l0P73qp6Tnc/93jX2I2Tljn5PvmCuf5Nuxh/8/B61yXEAgAAAAAAAADAiFXVAzKrDnS3zHJRnp3k65I8KskrhmH3SfKmqjqR/JPLcjQx6O1JLkzykCQXJ/lwZnk7l1bVJSewxkJrXzkoyV3m+rfvYvxtw+upe1hjUZWhM5NctYf5AAAAAAAAAADWRm92unrVYSzUfSAxviSzKkF3JHl0d79r7t7bqupDSV6Q5JwkT09y6V4XqKrzk3zvcPnGJN/e3Z8brq+qqjck+fPMjjN7QVX9bnd/+rh+mgXGUDno1rn+ybsYf8rwestuF+ju63dqmZ31BgAAAAAAAADAiFXVuUkeOVxedkxi0BEvTnL10H9aVd35OJYDYvLXAAAXS0lEQVT6yeH1c0l+eC4xKEnS3TcmeeZweUZm1YSWYgzJQZ+Z6++mVNPpw+tujiADAAAAAAAAAODwuHCu/8qtBnT3ZpJXDZdn5Ggy0a4MR5E9arh861CYZiuvTfKPQ//xe1ljL9Y+Oai7b01y43B51k5jq+qMHE0Oum6ZcQEAAAAAAAAAMDrnDa83Z3as13aunOs/fI9rPCRHT766crtB3X17kncfeeY4KxQttPbJQYMjpZrOrqqNHcads8UzAAAAAAAAAACHW2+Opy3XfYfXa7v7jh3GfWCLZ/a6xrHz7LTORpKv3OM6u7JTos06+dPMMrdOT/KgJO/ZZtz5c/137uP6dzrS+du//dt9nBYA4KhbP+dUVAAAAI66/vrtTh4AADh+x/w37zttN47puS23Jr3qKBa7LbfOX55ZVTuO3+HIrn+mqu6S5O7D5Y7PdfenqurmzHJV7rnbNQbz4xfFN38y1j2TvH+Pay00luSg1yf5qaH/xGyRHFRVJyX5/uHy00nevo/rf8mRzkMe8pB9nBYAAAAAALZ2z3v+5qpDAACm70uS/PWqg+BgXJW3rTqE43HVLsbsnD30+b5grr+bv9o+khx01z2ssdd1bp7r73WdXRnFsWLd/d4k7xguL66qh20x7Bk5Wpbppd39TwcSHAAAAAAAAAAAY3CXuf7tuxh/2/B66hLXuW2uv9d1dmUslYOS5EczOyrs1CRvqarnZVYd6NQk35PkycO4a5K8eJ/X/ssk5w79v0/yuX2eHxinM3M0U/XcJDesMBZguuw1wEGw1wAHwV4DHAR7DXAQ7DXAQVjmXnOnHD095y/3cV7W0w3Z+5FY6+LM7H+OxvyZZSfvYvwpw+stS1znlLn+XtfZldEkB3X3+6rqu5P8VpIvTPK8LYZdk+SC7v7MPq99W5I/2885gfE75mzLG/ZyliXAbtlrgINgrwEOgr0GOAj2GuAg2GuAg3AAe42jxA6J7r4jyVj/XbWMuOfzSXZzhNfpw+tujiA73nVOn+vvdZ1dGcWxYkd09xuT3D/Jz2eWCPTZJJ/OLHHnmUke2N3Xri5CAAAAAAAAAADWUXffmuTG4fKsncZW1Rk5mrhz3R6Xmk9s2nGdfH5lp72usyujqRx0RHf/dZKnDw0AAAAAAAAAAHbr6iTnJTm7qjaG6kpbOeeYZ/bi/dvMs9M6dyRZSkGcUVUOAgAAAAAAAACAE/Cnw+vpSR60w7jz5/rv3OMaVyW5fYt5Pk9VnZzkoUee6e7btxt7IiQHAQAAAAAAAABwWLx+rv/ErQZU1UlJvn+4/HSSt+9lge7+TJL/PVx+U1Vtd7TY45N84dB/3V7W2AvJQQAAAAAAAAAAHArd/d4k7xguL66qh20x7BlJ7jv0X9rd/zR/s6ouqqoe2nO2WepFw+tGkpdX1Z2OmePuSZ4/XH46ya/t7SfZPclBAAAAAAAAAAAcJj+a5JbMEnfeUlU/VVUPrapvqKpfSfKCYdw1SV58PAt099uS/PZw+bgkb62qx1XVg6vqiUneneRew/1ndfenjveHWWRjWRMDAAAAAAAAAMC66e73VdV3J/mtzI71et4Ww65JcsFwRNjxetIw/2OSfMPQ5m0meW53/8oJrLGQ5CCA49Td1yepVccBTJu9BjgI9hrgINhrgINgrwEOgr0GOAj2Gli+7n5jVd0/sypCFyQ5K8ntSa5N8jtJfqG7P3uCa9yS5IKq+t4kFyV5QJIvSvLxzI42+4XufteJrLEb1d3LXgMAAAAAAAAAAFiBk1YdAAAAAAAAAAAAsBySgwAAAAAAAAAAYKIkBwEAAAAAAAAAwERJDgIAAAAAAAAAgImSHAQAAAAAAAAAABMlOQgAAAAAAAAAACZKchAAAAAAAAAAAEyU5CAAAAAAAAAAAJgoyUEAe1RV96qqH6qq/1lVH6yqm6vq1qq6vqquqKonVNXGHub7qqr65aq6tqpuqaq/r6o/qaqn7GUeYHqG/eZFVXX1sNd8sqreW1U/XlWnrTo+YD1V1ddW1X+uqj+oquuq6raquqmqrqmqy6vqvD3O9y1V9drhd53bhtfXVtW3LOtnAMatql5QVT3XHrmLZ+w1wEJVdfeq+smqemdV3TDsFx+rqvdU1Qur6mG7mMN+A2yrqk6uqour6g+r6m/nPk99sKp+vaoeust57DVwyFTVParqsVV16fCdzI1zn4kuP475TngfqarTquonhu+UPznsZ1cP3znfa68xAeNW3b3qGABGo6ouTfLTSWrB0D9L8h3d/f8WzHdxkpcnOWWbIe9O8tju/sReYwXGraouSPLqJHfbZsgHkzymuz9ycFEB666qrkzyiF0M/c0kl3T37TvMVUl+OcmTd5jnV5P8YPtgCQyq6gGZfR6a/0OHb+juP95mvL0G2JWq+vdJfinJF+8w7IruvnCb5+03wI6q6p5J3pTkaxYM/fkkz9hqr7DXwOFVVTv9M/0b3X3RLufZl32kqu6d2Z52n22G/EOS7+3uN+8mLmD8VA4C2JsvzSwx6OYkv5XkiUkenuTBSb4vyVXDuAcn+aOquut2E1XVv83sF7hTknw8yY8k+TdJ/l2S1w7DHprktVVlv4ZDZPiPaq/JLDHopiTPTvJ1SR6V5BXDsPskedNO+wxwKH3Z8PqxJC9N8p1JHpLkYUmenuRvhvvfl+TyBXP91xz9Iup9SZ4wzPWE4TrD/efuQ9zABAyfW16RWWLQ3+3yMXsNsFBVfX+S384sMejvkvxskm9O8qAkF2T2ncpbk/zTDtPYb4BtDRXc5xOD/iLJRZl9lnp0kksz+044SX4syY9vM5W9BkiS65K85TifPeF9ZPjO+PdzNDHoFZl9t/x1mX3XfFNm3z3/TlXd/zjjBEZG5SCAPaiq5yf5RJJf6u7PbHH/Tkn+e5LvGt76L939z35BGz5sXp3k7CT/mORru/vDx4x5eZIfHi5/oLtftW8/CLDWqurtSR6Z5I4kj+judx1z/yeSvGC4/JnuvvRgIwTWVVX9fpJXJfm97v7cFvfvnuSdSf718NYjuvsdW4w7O7PfVTYyqwDyiO6+Ze7+aUmuzCwh+o4k5xz7uwxw+FTV0zL7S/oPJHldkp8abm1ZOcheA+xGVd03s/8QdkqSdyT51u7+h23GnrxVZUT7DbBIVX1Hkt8dLt+V5LxjP1NV1YOGe3dO8qkk9+juO+bu22vgEKuqn83sD8iv6u6PV9WXJ/nocHtXlYP2ax+pquck+Znh8ie7+4XH3H9Ykj8Z1nl7d3/jrn9QYLRUogDYg+5+Zne/YKvEoOH+5zJL6DnyRdR3bjPVt2eWGJQkP7fNh8CfyOxD5pE+cAhU1bmZJQYlyWXHJgYNXpzZh8QkeVpV3fkgYgPWX3c/trtfs1Vi0HD/xiTPmHtru99VfixHjwR66vwXUcM8n03y1OFyI8nTjj9qYAqGYziO/GHED+XoZ6Kd2GuA3XhZZolBNyZ5/HaJQUmyw5Gp9htgka+f6//cVp+puvvPM6vEkSRnJDnnmCH2GjjEuvtnuvv3u/vjJzDNCe8jw3fFPzpcXp3Zd8nHxvquJJcNl98wJD8CEyc5CGCfdfcnMis7myT33mbYhXP9y7eZ57OZHSuUJF9dVV+5LwEC625+f3jlVgO6ezOzyiDJ7MuoRy45JmBa/niu/89+VxnOtv+24fID3f3urSYZ3v/gcHnh8BxweP1ikrtm9hexf7xosL0G2I2qOiezIzCS5BeGROe9zmG/AXbj5Ln+R3YYN/9Hnqcc6dhrgBO1j/vII5N80dD/jeG75K1cPtd//J4DBkZHchDAchz5YLjdL13nDa8f7O4bdpjnyrn+w084Kvj/7d15rDZnWQfg311K6QJCaVlSllAo2g+xoQJSNktJNZElGkSkEqFQQNAQQSSFIgESFdQERNzYDFWSsqi1FiTKYtECCiVgDC1UKGhbFMtSCi20FG7/mDmc8fhu5zsH6Xm/60pO3mdmnnnOfH98v8w7557nYS/YyIdrk3xkQT/5AOyv6UPvWfcqxya509h+34zjUxvH75zkbju7LGCvqqrHJXlUki9l9VlPZQ2wip+ZtN+20aiqI6vqnlV11ApjyBtgFZdO2ndf0G/jBYtO8m+T/bIG2KndypGHzug3y0UZnkEnni/DAUFxEMAuq6rbJ9k3bn5ixvFbZrhhm3l8i+nxfXN7Aetk4//6p6br1s8gH4D9dfKkPeteZN+S45lzXBbBAaiqbpPkVePmmd191YqnyhpgFSeNn19JcklVPaGq/iVDMeKlSb5QVZdV1YvH5y2zyBtgFeckuWZsn1lVN9vaoapOTPLIcfPN3X3N5LCsAXZqt3JkpXHGZ88bs6HJIjgAKA4C2H3Py+aasG+dcfzOSTamebxiyViXT9p32eF1ATdxVXVokqPHzYX50N1fzuabHfIBWElVHZTk+ZNds+5VppniXgVY5reT3DHJB5K8YRvnyRpgFfcaPz+b5NVJ3pTkhC19jk3ykiQfrKpjZowhb4ClxgLn05N8PcmDk3y4qp5YVSdV1alV9eIMM3AckuRjSX5lyxCyBtip3cqRje1ru/vqFce5XVXdYmFPYM9THASwi6rqAUmePW5ekeQPZ3S71aT9tSVDXjtpz3sDDlgf28mHZDMj5AOwquck+ZGxfW53XzSjj3sVYCVV9ZAkT01yY5JndHdv43RZA6zituPn8Ul+KcnVSZ6R5PZJDk1y/yTvHPvcO8nbxmLoKXkDrKS7z01yvwwFz/dJcnaSDyZ5V4YixOsyFAU9pLv/a8vpsgbYqd3KkY1xtvN8edY4wJpRHASwS6rqDkn+PMOsQZ3kSd193Yyuh07aNywZ9vpJ+7CdXSGwB2wnH5LNjJAPwFJVdXKSl4+b/53kmXO6ulcBlqqqQ5K8NsOsqK/s7n/d5hCyBljFEePnLZJ8K8lPdPdruvuq7r5+LHR+VDYLhB6U5DFbxpA3wEqq6uZJfi7Jo7M58/vUHZKcluRhM47JGmCnditHNsbZzvPlWeMAa0ZxELCWqurgqupd+Dl9xd93qyTvyLBkWJKc1d3vndP9G5P2IUuGnk7j+PVVrgXY07aTD8lmRsgHYKGq+sEk52YoYr4+yeO6+/NzurtXAVZxVpJ9Sf4jyUv343xZA6ximhVv6+5/2tqhu7+dYYn3DactGEPeADNV1RFJ3p3khUmOyrB06r4MuXDrJD+e5MIMM5adX1W/vGUIWQPs1G7lyMY423m+PGscYM0oDgLYoao6NMl5Se477npFd798wSlfnbSXTdN4xKS9yhSQwN62nXxINjNCPgBzVdWxSf4uyZEZ3rg/rbvft+AU9yrAQlV1fJIXjJvP6u5rF/WfQ9YAq5hmxTvnderujye5cty8/4Ix5A0wz0uT/OjYPqO7z+zuT3T3Dd19TXe/K8kpSf4+w6xCr6iqEybnyxpgp3YrRzbG2c7z5VnjAGvm4O/1BQB8N3T3jVW1bxeG+s9FB6vq4CRvzfDFMEle393PXTLmFZP2nef2Gtxl0r58SV9gj+vub1TVF5IcnSX5UFVHZvPLm3wAZqqqYzK8/XpMhmVPn9Ld5y45zb0KsMxzMryFelmSw6vq8TP63HvSfnhV3XFsnz8WE8kaYBWXJ9nIjysWdRz73inJ7bfslzfAQlVVSZ48bl7a3WfP6jc+c35RhhmEDhrPec54WNYAO7VbOXJFkgckOaKqbtPdV68wzlXdff2CfsAaUBwErK3u/sR3c/yqOijJn2VYgzpJ3pLkF1a4rq9V1eUZbrqOX9J9evyS/blOYM+5JMlDkxxXVQd3941z+skHYKGqOjrJu5Lcfdz1rO7+0xVOvXjSdq8CzLIx9fzdk5yzQv8XTdrHJrk2sgZYzcezORPQzZb03Ti+9TuUvAGWuUOS247tjy7p+5FJe5oZsgbYqd3KkYuT/PSk3/9ZljX5zsvv95gzBrCGLCsGsP9ek2TjDdm3J/n5cZ37VVw4fv7A5A3aWU6etN+/zesD9qaNfDgim8sVziIfgLmq6tZJ/jbJvcZdz+/uP1jx9M8k+dzYPnlRx2xOu39lks9u5xqBA56sAVbxD5P2Peb2GmwURF+5Zb+8AZaZFhUue6n+5nPOkzXATu1Wjlw4aS8a537ZnJne82U4ACgOAtgPVfWKJE8dN9+T5LHd/c1tDPFXk/bpc37H4UkeN25e3N2Xbvc6gT1pmg9PntVhnLnsiePm1RnWuwdI8p17iHck+eFx129092+ten53d5Lzxs3jq+qkOb/npGy+qXbeeB5wAOju07u7Fv0keenklFMmxz47jiFrgFX8dZKN5y2Pmdepqk5OctS4+Y/TY/IGWMGXklwzth84zqYxz/QP7Z/ZaMgaYKd2MUcuSPKVsf2kcenEWU6ftJctQQ+sAcVBANtUVS/J5lrSH0jyk/uxFuu5ST49tl9QVbPefvudJEdO2sABoLs/lM2H2WdU1QNndHtukn1j+1XbLE4E1lhVHZLhPuPB465Xdfev7cdQv5vNt2BfXVWHbfk9hyV59bh549gfYLtkDbBQd38xyevHzR+rqsdv7VNVt8r/zofXzBhK3gBzjbPBv2PcPCbJC2f1q6ojk0xfvHj7li6yBtipHedId9+Q5PfGzX1JfnVrn/GZ8xnj5vu6+8M7v3Tgpq4UJQOsrqqelc2bqiuT/Gw2K7Dn+eSsP9xX1SOSnJ+hUPPzSX49yYcyFAQ9LZtrwl6Y5GHd/a0d/wOAPaGqTswwlethSb6W5DczzA50WIblDJ8+dr00yf26+6vfi+sEbnqq6i+y+Vb9e5M8O8miL303zJudsKpeluT54+ZHMzwE/3SGJT3OTHLieOxl3X3WDi8dWDPjSxUvHjdP6e4L5vSTNcBCVXW7JBcluWuGP4L9cZK/zDDLxw9lyIqNt+f/qLt/cc448gaYq6qOT/KRJIePu85PcnaSy5IcmuSkDN+v7joef093nzpjHFkDB6iqekiS4ya7js7mi9/vz2bBc5Kku984Z5wd58hYPH1Rku8fd702yZuTfD3JKUnOSnLLcftB3f2xVf6NwN6mOAhgG6rqgixf63WrYzemzp8x3tOS/H6SQ+ac+6Ekj+zuL2zzdwJ7XFU9OsmbknzfnC6XZsiHT/3/XRVwU1dV2/2C9+/dfbc5Yx2U5HVJnrLg/Dckefr4pi3Ad2yjOEjWAEtV1b4MS4wdt6DbnyR5xryZVeUNsExVnZrknAx/0F/kvUke291fnjGGrIEDVFW9McmTVu0/Lsc8a5xdyZGqOi7J3yS555wu1yR5QndvnQUNWFOWFQP4Huru1yW5b4YbvcuSfCPJFzPMFvTMJA9WGAQHpu4+P8kJSV6ZoRDouiRXZ3jj48wkJyoMAr6buvvb3X1GkkdmWPP+c0luGD/PS/KI7n6qB9rATsgaYBXdfUmS+yR5XpJ/TvKlDFlxRZK3JHl4d5+xaMlleQMs093vzjAT2ZlJLkhyVZJvZphZ4zNJ3prkp5KcOqswaBxD1gA7sls5Mj47PjFDpl2U4dnydUk+meGZ8wkKg+DAYuYgAAAAAAAAAABYU2YOAgAAAAAAAACANaU4CAAAAAAAAAAA1pTiIAAAAAAAAAAAWFOKgwAAAAAAAAAAYE0pDgIAAAAAAAAAgDWlOAgAAAAAAAAAANaU4iAAAAAAAAAAAFhTioMAAAAAAAAAAGBNKQ4CAAAAAAAAAIA1pTgIAAAAAAAAAADWlOIgAAAAAAAAAABYU4qDAAAAAAAAAABgTSkOAgAAAAAAAACANaU4CAAAAAAAAAAA1pTiIAAAAAAAAAAAWFOKgwAAAAAAAAAAYE0pDgIAAAAAAAAAgDWlOAgAAAAAAAAAANaU4iAAAAAAAAAAAFhTioMAAAAAAAAAAGBNKQ4CAAAAAAAAAIA1pTgIAAAAAAAAAADWlOIgAAAAAAAAAABYU4qDAAAAAAAAAABgTSkOAgAAAAAAAACANfU/Fk9Xo30BqhkAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "init()\n", "plot_ρs()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Run the simulation until converged" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAACQcAAAQKCAYAAAD3zl6bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAewgAAHsIBbtB1PgAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdfZDtd10f8Pdn7+Ym4flJTEiw8mQDijYUKOpI6hPVilUYddQWRUClD0zrdCqZgoVqO63WzmhTqlApKGNrBaTQlmCx1ogIJXRaawsRi6IJAUKgQCDJvbnnfPrHntt7env3npvf7tmzv9++XjNn9vs7v+/DZ8/ZvRkO7/1+q7sDAAAAAAAAAABMz9amCwAAAAAAAAAAANZDOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAA4FKrqflX1kqp6V1XdVlUnFo8PVdWrq+qLN10jAAAAAIxNdfemawAAAACOuKr6xiT/PMkV5+l2d5LndPcbDqYqAAAAABg/4SAAAABgo6rq5UletvTU7ye5McmJJF+a5E8v3bszyZO6+/cOrEAAAAAAGDHhIAAAAGBjqurvJPm7i8tbknx/d7/trD7PTPKGJBcvnnptd3/fwVUJAAAAAOMlHAQAAABsRFV9TZJfS1JJ/jjJV3T3h3fp+yNJfnRxeUeSB3X3/EAKBQAAAIAREw4CAAAADlxVHUvyviRflKSzEwx693n6X5nk5qWnrnK0GAAAAACstrXpAgAAAIAj6TuyEwxKkn99vmBQknT3LUk+s/TUA9dVGAAAAABMiXAQAAAAsAnft9R+xQWOuWepbStkAAAAALgAjhUDAAAADlRV3T/J/05yLMntSR7eKz6gqKrtJHcl2V48dXl3f3SthQIAAADABNg5CAAAADhoX5GdYFCS/NaqYNDCY3MmGHRHko+tozAAAAAAmBrhIAAAAOCg/Zml9vsvcMxXLbV/8wIDRQAAAABw5AkHAQAAAAft6qX2bRc45llL7V/dx1oAAAAAYNKEgwAAAICDthwOumhV56p6bJI/t7i8O8kvrqMoAAAAAJgi4SAAAADgwFTVQ5L8iaWnnnABw/5xznyG8Zru/uS+FwYAAAAAEyUcBAAAABykq8+6flZVPXi3zlV1bZK/sLj8ZJIfWVdhAAAAADBFwkEAAADAQVoOB51M8sAkv1RVD1ruVFUPrKpXJPkHi6dmSZ7T3Z84mDIBAAAAYBq2N10AAAAAcKQsh4NekuQfJnlGkj+qqv+Y5PYkVya5Jsl9Fv1mSV7Y3W89yEIBAAAAYAqquzddAwAAAHBEVNX7kjx+cfnI7BwZdl123934liTf391vO4DyAAAAAGByhIMAAACAA1FV90nymSTHknyyux+6eP4rk/z1JE9NcnmSTyd5X5I3Jnl1d9+5mYoBAAAAYPwcKwYAAAAclC/LTjAoSX7n9JPd/c4k79xIRQAAAAAwcbtt2Q0AAACw365eav+3jVUBAAAAAEeIcBAAAABwUJbDQb+zay8AAAAAYN8IBwEAAAAHRTgIAAAAAA5YdfemawAAAAAmrqq2k3w2ycVJ7klyv+4+udmqAAAAAGD67BwEAAAAHIQnZCcYlCQ3CQYBAAAAwMEQDgIAAAAOgiPFAAAAAGADHCsGAAAAAAAAAMCRUFVPSvINSb4qyZckeXiSe5LcmuS3k7y6u9+xD+u8PMnLLrD7V3f3b+x1zd1sr2tiAAAAAAAAAAA4LKrqhiRPP8et40ket3h8b1W9LskLuvvkQda3LsJBAAAAAAAAAAAcBVcsvt6a5PVJ3pHkj5McS/LlSf7mos9zspOp+e59WveJK+7/4T6tc06OFQMAAAAAAAAAYPKq6t8l+YUkb+zu2TnuPyzJO5N80eKppw89Ymz5WLHurkEF75OtTS4OAAAAAAAAAAAHobuf2d2/fK5g0OL+7dnZPei0bzuYytZLOAgAAAAAAAAAAHb8xlL7MZsqYj8JBwEAAAAAAAAAwI7jS+35xqrYR8JBAAAAAAAAAACw45ql9k37MWFVvb2qPlFVJ6vqtqr6jaq6tqoevB/zr1y/uw9inVGrqouTPHFx+fEk5zx7DgAAAAAAAAAOuWNJPm/R/t3uPrHJYlivqtpOctmm6xjoslxARqO7b9mvBatqK8m7kjx18dRTuvu9A+d6eZKXrej2qSTP7e43D1njQm2vc/IJeWKSGzddBAAAAAAAAADso6ckGRR8YDQuS3LzpotYs9rHuX4oZ4JBbxoaDFryu0n+TZL3JLk1yUVJ/mSSv5jkGUkelOSNVfXN3X39HtfalZ2DLkBVPTnCQQAAAAAAAABMy+BdURiHqroyEw8Hdfe+hIOq6pokv5adjXZuS/Kl3f2xPcz3oO7+1Hnu/2CSn11c3prksd1919D1zsfOQRfm46cb73nPe3L55ZdvshYAYKKe8+Qf2XQJAAxR+/mHSXDI+KMygI163Xt/bNMlAAAT9JGPfCRPferpjVHO/H/hTN+7r39kLn/4sU2XsdJHbpvlad/4f/NMT0ny0XWvWVVfnORN2cnRnEjyHXsJBiXJ+YJBi/uvXGxW84Ikj0jy7CS/uJc1dyMcdGH+7/l1l19+ea688spN1gIATNQlx+636RIAGGJra9MVwPrM55uuAOBI81k0AHAAZqu7MBWXP/xYrnzERZsu4976aHffss4FqupRSf5Dkgdn53fiu7r7hnWuueSV2QkHJck1WVM4yCeYAAAAAAAAAAAcOVX1iOwcJfaIJJ3ked39pgMs4X1L7SvWtYidgwAAAAAAAAAAJm6ezjyHf5fgeQ7mmPOqeliStyd59OKpF3X3LxzI4ktlHMQidg4CAAAAAAAAAODIqKoHJvnVJE9YPHVtd79iA6U8Yal967oWEQ4CAAAAAAAAAOBIqKr7JPn3SZ60eOrvd/ePb6icH1xq37CuRRwrBgAAAAAAAAAwcbOeZ9aH/1ixddZYVceTvCnJVy6e+unufumAeZ6b5DWLy7/b3S8/6/4Tk9zV3f/rPHP8YJLnLy4/uqhrLYSDAAAAAAAAAAA4Cv5Vkmcs2r+e5NVV9SXn6X+yuz8wYJ0/neTnquo/Jbk+ye8m+UR2cjpXJflLSb5+0XeW5Ae7+3MD1rkgwkEAAAAAAAAAABwFz15qf02S/76i/x8l+cKBax1L8nWLx24+keT53f2WgWtcEOEgAAAAAAAAAADYP2/NzpFhX57k6iSfn+ShSSrJJ5P8TpK3JXltd39m3cUIBwEAAAAAAAAATNw8nXl602WstM4au7v2aZ7XJnntee7fluRfLB4bt7XpAgAAAAAAAAAAgPUQDgIAAAAAAAAAgIlyrBgAAADsxXw+bNyWv9fhAA39OQUAAAAmozPPPIf/M4IeQY1j45NIAAAAAAAAAACYKOEgAAAAAAAAAACYKOEgAAAAAAAAAACYqO1NFwAAAAAAAAAAwHrNujPr3nQZK42hxrGxcxAAAAAAAAAAAEyUcBAAAAAAAAAAAEyUY8UAAAAAAAAAACZuns48h//IrjHUODZ2DgIAAAAAAAAAgIkSDgIAAAAAAAAAgIkSDgIAAAAAAAAAgIna3nQBAAAAAAAAAACs1zzJLL3pMlaab7qACbJzEAAAAAAAAAAATJRwEAAAAAAAAAAATJRjxQAAAAAAAAAAJm6eznwUx4od/hrHxs5BAAAAAAAAAAAwUcJBAAAAAAAAAAAwUcJBAAAAAAAAAAAwUdubLgAAAAAAAAAAgPWadWfWvekyVhpDjWNj5yAAAAAAAAAAAJgo4SAAAAAAAAAAAJgox4oBAAAAAAAAAEzcfPE47MZQ49jYOQgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZq+6AWqqqHJXlekm9J8pgkD07yiSQ3J/nNJL/S3e9aMcc3JPmBJE9N8nlJPp7kPUle1d1vW1/1AAAAAAAAAADjNU9nlt50GSvNR1Dj2BxIOKiqvj3JzyR56Fm3Ll88nprkcUm+dZfxleRnsxMMWnZFkmcleVZVvSrJC7vbTwkAAAAAAAAAAOQAwkFV9T1JXpOdI8xuy05I6LeSfDLJZdnZReibk9xznmn+Xs4Eg/5rkp9I8sHF2B9OcvXi/seTvHTfvwkAAAAAAAAAABihtYaDqurxSV6VnWDQO5J8c3d/+hxdr6uq47vM8djsBICS5L1Jnt7ddy2ub6yqtyS5IcmTk7y4ql7T3R/cz+8DAAAAAAAAAADGaGvN81+X5OIktyd59i7BoCRJd5/c5dYP5UyI6UVLwaDT4+5M8qLF5XaSv7GnigEAAAAAAAAAJmaWZNYjeGz6hZqgtYWDquqqJF+7uPyn3X37gDkqybcsLm/q7nefq9/i+d9bXH7rYhwAAAAAAAAAABxp69w56NuX2q8/3aiqB1fV46rqoRcwx6OSXLFo37Ci7+n7Vyb5wgstEgAAAAAAAAAApmqd4aCnLb5+Osn7q+ovVtXvJPlkkg8kub2q/qCqXlZV99tljscvtW9asd7y/cfv2gsAAAAAAAAA4IiZj+jB/tpe49xPWHz9UJLrkvzVc/R5VJKXJ/m2qvpz3X3rWfcfudS+ZcV6N+8ybqWqunJFl8vuzXwAAAAAAAAAAHAYrDMc9JDF16uSfFmSTyW5NsmvJPlMkicm+dEk35jkS5K8vqq+qruXQ2D3X2p/dsV6n1tq77YT0W5uXt0FAAAAAAAAAADGZZ3Hit138fXiJLMk39jdr+zuj3f3ie5+b5JnJrl+0e8rkjz7rDkuWWqfXLHeiaX2pQNrBgAAAAAAAACAyVjnzkF350xA6PXd/e6zO3T3vKr+VnZ2D0qS70ryhrPmOO34ivUuXmrfdS9rXXUM2WVJbryXcwIAAAAAAAAAHArzVGapTZex0nwENY7NOsNBd+RMOOj63Tp19/+sqg8nuSLJU84xx2mrjgq771J71RFkZ9dwy/nuV/nBAwAAAAAAAABgfNZ5rNjNS+3zhm+W+j78rOeXx125Yo7l3X9u3rUXAAAAAAAAAAAcEevcOeh/5sxOQMdW9D19/9RZz79vqX3VijmW779/RV8AAAAAAAAAgCNj3juPw24MNY7NOncO+s2l9mNW9H304uuHz3r+D5Pcumhfs2KOpy/N8aFVxQEAAAAAAAAAwNStMxz0liT3LNrP3q1TVV2T5KGLy3cs3+vuTvLmxeVVVfW0XeZ4Ws7sHPTmxTgAAAAAAAAAADjS1hYO6u5PJPm5xeXXV9V3nt2nqu6f5KeWnnrlOab6qZw5buy6qrr0rDkuTXLd4vLUWfMBAAAAAAAAAMCRtb3m+V+W5JuSfEGS11XVVyb5lSSfSfLEJC/OmR1/fqa7bzx7gu7+QFX9ZJJrkzw5yTur6seTfDA7x5W9OMnVi+7/qLt/f43fDwAAAAAAAADA6MxSmaU2XcZKY6hxbNYaDuruj1fVN2TniLHHJvlri8fZ/kWSv36eqV6S5OFJnpedINAvnaPPq5O8dE8FAwAAAAAAAADAhKztWLHTuvv9Sf5Ukr+V5D8n+WSSk0luSfKvk3xNdz+/u+85zxzz7n5+dnYhenOSWxdz3Lq4/vPd/YLunq/1mwEAAAAAAAAAgBFZ97FiSZLu/lySn1w89jLPW5O8dV+KAgAAAAAAAAA4IhwrdnStfecgAAAAAAAAAABgM4SDAAAAAAAAAABgooSDAAAAAAAAAABgorY3XQAAAAAAAAAAAOvVncy7Nl3GSt2brmB67BwEAAAAAAAAAAATJRwEAAAAAAAAAAAT5VgxAAAAAAAAAICJm6Uyy+E/VmwMNY6NnYMAAAAAAAAAAGCihIMAAAAAAAAAAGCihIMAAAAAAAAAAGCitjddAAAAABxJ8/mwcVv+zudIG/pzAwAAABx5s2xlNoI9ZMZQ49h4RQEAAAAAAAAAYKKEgwAAAAAAAAAAYKIcKwYAAAAAAAAAMHHdlXnXpstYqUdQ49jYOQgAAAAAAAAAACZKOAgAAAAAAAAAACZKOAgAAAAAAAAAACZqe9MFAAAAAAAAAACwXrNUZqlNl7HSGGocGzsHAQAAAAAAAADARAkHAQAAAAAAAADARDlWDAAAAAAAAABg4ma9lVkf/j1kxlDj2HhFAQAAAAAAAABgooSDAAAAAAAAAABgooSDAAAAAAAAAABgorY3XQAAAABwL8znw8du+RuhQ2Mv7yMAAADAAPNU5iPYQ2ae2nQJk3P433UAAAAAAAAAAGAQ4SAAAAAAAAAAAJgo4SAAAAAAAAAAAJio7U0XAAAAAAAAAADAes1TmaU2XcZK8xHUODZ2DgIAAAAAAAAAgIkSDgIAAAAAAAAAgIlyrBgAAAAAAAAAwMTNeiuzPvx7yIyhxrHxigIAAAAAAAAAwEQJBwEAAAAAAAAAwEQJBwEAAAAAAAAAwERtb7oAAAAAAAAAAADWa57KPLXpMlYaQ41jY+cgAAAAAAAAAACYKOEgAAAAAAAAAACYKMeKAQAAAAAAAABM3DxbmY1gD5n5CGocG68oAAAAAAAAAABMlHAQAAAAAAAAAABMlHAQAAAAAAAAAABM1PamCwAAAAAAAAAAYL1mvZVZH/49ZMZQ49h4RQEAAAAAAAAAYKKEgwAAAAAAAAAAYKIcKwYAAAAAAAAAMHHzVOYj2ENmntp0CZNz+N91AAAAAAAAAABgEOEgAAAAAAAAAACYKOEgAAAAAAAAAACYqO1NFwAAAAAAAAAAwHrNuzLr2nQZK81HUOPY2DkIAAAAAAAAAAAmSjgIAAAAAAAAAAAmyrFiAAAAMCZb/s5nEoa+j/P5/tYBAAAAHBmzbGU2gj1kxlDj2HhFAQAAAAAAAABgooSDAAAAAAAAAABgooSDAAAAAAAAAABgorY3XQAAAAAAAAAAAOs1763M+/DvITOGGsfGKwoAAAAAAAAAABMlHAQAAAAAAAAAABPlWDEAAAAAAAAAgImbZSuzEewhM4Yax8YrCgAAAAAAAAAAEyUcBAAAAAAAAAAAEyUcBAAAAAAAAAAAE7W96QIAAAAAAAAAAFiveZJZ16bLWGm+6QImyM5BAAAAAAAAAAAwUcJBAAAAAAAAAAAwUY4VAwAAgCTZ8vczjMBYfk7nNgAHAACAw2aercxHsIfMGGocG68oAAAAAAAAAABMlHAQAAAAAAAAAABMlHAQAAAAAAAAAABM1PamCwAAAAAAAAAAYL1mvZVZH/49ZMZQ49h4RQEAAAAAAAAAYKKEgwAAAAAAAAAAYKKEgwAAAAAAAAAAYKK2N10AAAAAAAAAAADrNU9lntp0GSuNocaxsXMQAAAAAAAAAABMlHAQAAAAAAAAAABMlGPFAAAAAAAAAAAmbt5bmfXh30NmPoIax8YrCgAAAAAAAAAAEyUcBAAAAAAAAAAAEyUcBAAAAAAAAAAAE7W96QIAAADgnLb8PQuM1iZ+f+fzg18TAAAARmSWymwEe8jMUpsuYXIO/7sOAAAAAAAAAAAMIhwEAAAAAAAAAAAT5VgxAAAAAAAAAICJm3dl3of/yK4x1Dg2dg4CAAAAAAAAAICJEg4CAAAAAAAAAICJEg4CAAAAAAAAAICJ2t50AQAAAAAAAAAArNc8W5mNYA+Z+QhqHBuvKAAAAAAAAAAATJRwEAAAAAAAAAAATJRjxQAAAAAAAAAAJm7eW5n34d9DZgw1jo1XFAAAAAAAAAAAJko4CAAAAAAAAAAAJko4CAAAAAAAAAAAJmp70wUAAAAwYVv+JuVQqdp0BZzWvekKpmfovzfz+f7WAQAAAIfULJVZDv/nQ2OocWx8SgsAAAAAAAAAABMlHAQAAAAAAAAAABPlWDEAAAAAAAAAgImb91bmffj3kBlDjWPjFQUAAAAAAAAAgIkSDgIAAAAAAAAAgIkSDgIAAAAAAAAAgIna3nQBAAAAAAAAAACs1zzJLLXpMlaab7qACbJzEAAAAAAAAAAATJRwEAAAAAAAAAAATJRjxQAAAAAAAAAAJm7eW5n34d9DZgw1jo1XFAAAAAAAAAAAJko4CAAAAAAAAAAAJko4CAAAAAAAAAAAJmp70wUAAABwgLb8jci+qzrYcVsDx+3FQX+PB637YMftxfyAa93E93jQDvrfxfn8YNcDAACAhVlvZdaH//PBMdQ4Nl5RAAAAAAAAAACOhKp6UlX97aq6vqpurqoTVfXZqvpAVb22qr5qDWt+Z1X9alV9pKrurqoPVdXrqupp+73Wudg5CAAAAAAAAACAyauqG5I8/Ry3jid53OLxvVX1uiQv6O6Te1zvkiSvT/LMs279icXju6vq5d39Y3tZZxXhIAAAAAAAAACAietU5jn8x8D3emu8YvH11uyEdt6R5I+THEvy5Un+5qLPc7KTqfnuPa736pwJBv2nJD+9WPuJSf52ksck+dGq+kh3/9we19qVcBAAAAAAAAAAAEfBTdkJ5byxu2dn3Xv3Ysegdyb5oiTfVVU/093vGLJQVV2TM+Gif5vkWUtr3lhVb0nyX5J8QZKfqKo3dPenhqy1ytY6JgUAAAAAAAAAgMOku5/Z3b98jmDQ6fu3Z2f3oNO+bQ/L/fDi6yzJXzl7zcVaL15cPjjJ8/ew1nkJBwEAAAAAAAAAwI7fWGo/ZsgEVXW/JF+7uHx7d9+yS9dfSfKZRfvZQ9a6EI4VAwAAAAAAAACYuFlvZdaHfw+ZQ1Dj8aX2fOAcT01y8aJ9w26duvtkVb07yTOSPLWqLuruewauuauNv6IAAAAAAAAAAHBIXLPUvmngHI+/F3Ocvr+d5HED1zsvOwcBAAAAAAAAAHAYXVZV5+1wniO77rWq2kpy7dJTvzxwqkcutVfVd/NZ4943cM1dCQcBAAAAAAAAAHAY3XgBfc6fHrp3fig7R4IlyZu6+70D57n/UvuzK/p+bql9v4Hrnddaw0FV1RfY9Ybu/rMr5vqGJD+QnTfh85J8PMl7kryqu9+2lzoBAAAAAAAAAKZsnsq89zNHsx7zfc36XLiquibJP1xc3pbkL+9hukuW2idX9D2x1L50D2vu6tDvHFQ7+0P9bHaCQcuuSPKsJM+qqlcleWF3X2gYCQAAAAAAAACAw+0pST667kWq6ouTvCk7OZoTSb6juz+2hynvXmofX9H34qX2XXtYc1cHFQ76mST/7Dz3P3eee38vZ4JB/zXJTyT5YJLHJPnhJFcv7n88yUv3XCkAAAAAAAAAAIfBR7v7lnUuUFWPSvIfkjw4ySzJd3X3DXuc9o6l9qqjwu671F51BNkgBxUOuq27/8e9HVRVj81OAChJ3pvk6d19OiV1Y1W9JckNSZ6c5MVV9Zru/uC+VAwAAHCYbW1tuoL1q4HbBw8dtzV03LD3oo8NfA+PHRs2Lkm2h43t7YG1DnwveuC4Grqh8MBxdWo+bL0kOTUbNm42bFzNBtY6HzruYN+LPY89SEP//R76XgAAAMDCLJVZDv/nirMDPFasqh6R5NeSPCJJJ3led79pH6ZeDjRdmZ3My24eudS+eR/W/v8c9nf9h3ImwPSipWBQkqS770zyosXldpK/cYC1AQAAAAAAAAAwQlX1sCRvT/LoxVMv6u5f2Kfp37fUvmpF39P3TyX5X/u0/v/j0IaDqqqSfMvi8qbufve5+i2e/73F5bcuxgEAAAAAAAAAwP+nqh6Y5FeTPGHx1LXd/Yp9XOLGJCcX7WvOU8fxJE87Paa7T+7Wdy8ObTgoyaOSXLForzrL7fT9K5N84boKAgAAAAAAAABgvKrqPkn+fZInLZ76+9394/u5RnffkeQ/Li6/rqqu3KXrs5M8YNHej+PMzumgwkHfXlW/V1V3VdUdVfX7VfXzVfXV5xnz+KX2TSvmX77/+F17AQAAAAAAAAAcQd2V+Qge3es7MGqxU8+bknzl4qmf7u6XDpjnuVXVi8fLd+n2k4uv20leUVXHzprjYUlOh5I+leTn7m0dF2p7XROf5QlnXT928fieqvo3SZ7b3Z8+q88jl9q3rJj/5l3GXZDzJLROu+zezgkAAAAAAAAAwKHyr5I8Y9H+9SSvrqovOU//k939gSELdfevV9UvJfnOJH8hydur6qeS3JrkiUlekuQLFt2v7e7/PWSdC7HucNCdSd6Sna2Sbkry2SSfl53z1F6Y5KFJvjXJm6vq67v7nqWx919qf3bFOp9bat9vQJ03r+4CAAAAAAAAAMCIPXup/TVJ/vuK/n+U5Av3sN7zsnNs2J9P8tWLx7J5kh/r7lfuYY2V1h0OuqK7P3WO599eVdcluT7J1dkJC/3lJP9kqc8lS+2TK9Y5sdS+dEihAAAAAAAAAABTNc9W5tnadBkrjaHGC9XddyX5pqr67iTPTfJlSR6U5GNJ3pHkn3b3u9Zdx1rDQbsEg07f+1hVfVuS9yc5nuRF+X/DQXcvtY+vWOripfZd97bOrD6K7LIkNw6YFwAAAAAAAACAQ6C7a5/meW2S196L/v8yyb/cj7WHWPfOQefV3X9QVW9P8k1JHltVj+juWxe371jquuqosPsutVcdQXauOm453/2qffnZAAAAAAAAAACAA3UY9mJ631L7iqX2cmDnyhVzLO/8c/OeKwIAAAAAAAAAgAnY6M5BC7tty7McGrpqxRzL99+/t3IAAAAAAAAAAKZl1pXZ/pyqtVZjqHFsDsPOQU9Yat+61P7DpetrVszx9MXXDyf50P6UBQAAAAAAAAAA47bRcFBVPTrJ1y8u/6C7P3z6Xnd3kjcvLq+qqqftMsfTcmbnoDcvxgEAAAAAAAAAwJG3tmPFquqbk1zf3ad2uf/5Sd6Q5KLFU684Ry764j4AACAASURBVLefSvL92anzuqp6enfftTTHpUmuW1yeWvQHAAA4eFuHYWPWQ6j2sAXwsYGv6bFjg4b1xRet7rSP42aXDhs3v2TY95ckpwaOnR8f9j7Otwe+/0N/bAb+udDWqWEDt04O//uk7btnw9YcOO7YXfcMGlcnDnZcZsO+v52x82HjxvJ3Zgf935n5wNcTAACAQ2velfkIjuwaQ41js7ZwUHZCOxdV1RuTvCs7x33dleRhSf5skhcmeeii72/lHOGg7v5AVf1kkmuTPDnJO6vqx5N8MMljkrw4ydWL7v+ou39/Xd8MAAAAAAAAAACMzTrDQUnyiCQvWjx288YkL+juE7vcf0mShyd5XnaCQL90jj6vTvLSPdQJAAAAAAAAAACTs85w0PcmuSbJlyd5dHZ2DHpAks8muTnJbyf5+e5+1/km6e55kucvdiD6gSRPWcx1e5Ibk7yyu69f1zcBAAAAAAAAAABjtbZwUHffkOSGfZzvrUneul/zAQAAAAAAAAAcFd1bmffWpstYqUdQ49h4RQEAAAAAAAAAYKKEgwAAAAAAAAAAYKLWdqwYAAAAAAAAAACHwyyVWWrTZaw0hhrHxs5BAAAAAAAAAAAwUcJBAAAAAAAAAAAwUcJBAAAAAAAAAAAwUdubLgAAAAAAAAAAgPWadzLv2nQZK8170xVMj52DAAAAAAAAAABgooSDAAAAAAAAAABgohwrBgAAcNrWEfj7iRq4bfDQcdvHho1L0pccHzRufp9h407d/+JB404+cNj/tD55/2E/byfvN3zr51P3GTZ2PuwlzXzo2z/0V3E+cLnZwHEnh41Lku07h/3cHP/ssH21j98x7Of7+KdPDRq3fceJQeO27hz+otbdA8eeGvgD0AP3OB867qDt5b+J84G/jAAAAKzVvLcy78P/GegYahwbrygAAAAAAAAAAEyUcBAAAAAAAAAAAEyUcBAAAAAAAAAAAEzUsAPuAQAAAAAAAAAYjU5lntp0GSv1CGocGzsHAQAAAAAAAADARAkHAQAAAAAAAADARDlWDAAAAAAAAABg4mZdmfXhP7JrDDWOjZ2DAAAAAAAAAABgooSDAAAAAAAAAABgooSDAAAAAAAAAABgorY3XQAAAAAAAAAAAOs1763M+/DvITOGGsfGKwoAAAAAAAAAABMlHAQAAAAAAAAAABMlHAQAAAAAAAAAABO1vekCAAAA9t3WEfg7iKph444Ne236+EXDxt3n4kHjkuTUAy4ZNO7uhx0fNO6uhw57bU48ZNh7cfL+PWjc7D7DxiXJ/JL5oHG9PWxcjg2sdeCPd4a+NLNhC9ap4f/WbN09bM1jdw4bd/yOYR8BXfzJY4PGXfqJYf9mXHL7yUHjkmT7M3cPGld3nhg27uQ9g8ZlNvD3qYf/7h+4of8dng98bQAAALgg81TmPfSDl4MzH/zhELs5Ap+YAwAAAAAAAADA0SQcBAAAAAAAAAAAE+VYMQAAAAAAAACAievUKI7s6hHUODZ2DgIAAAAAAAAAgIkSDgIAAAAAAAAAgIkSDgIAAAAAAAAAgIna3nQBAAAAAAAAAACs17yTedemy1hp3puuYHrsHAQAAAAAAAAAABMlHAQAAAAAAAAAABPlWDEAAAAAAAAAgImb91bmffj3kBlDjWPjFQUAAAAAAAAAgIkSDgIAAAAAAAAAgIkSDgIAAAAAAAAAgIna3nQBAAAAR1rVsHHbxwYN6+MXDRo3f8Clg8adeNiwcUnyuc8f9j9Z7/r8Ya/p3Q/pQeNOPeDUoHG5dDZo2LGLh41LkouOzQeN29oa9tpUHey47mHv/dBx8/nA398ks9mwv9c6eWLY7/7Ju4aNu/shw8adGDju0gdcMmhcktz3Y8P+zbj49mG1bn3mrkHj6uQ9g8bl1MDf/R72+wQAAMD0zLsyH/g5yEEaQ41jY+cgAAAAAAAAAACYKOEgAAAAAAAAAACYKMeKAQAAAAAAAABM3DyVeQ7/kV1jqHFs7BwEAAAAAAAAAAATJRwEAAAAAAAAAAATJRwEAAAAAAAAAAATtb3pAgAAAAAAAAAAWK/uyrxr02Ws1COocWzsHAQAAAAAAAAAABMlHAQAAAAAAAAAABPlWDEAAAAAAAAA/g979x9reVrXB/z9nHPmzszuLBV2hcVdEAtN2aakJQVCAwGt1FjECrQ1Vm3cAIJtQiN/CLTauKL/iJhiaBuxUtcaldj6Y9GAlrZIUlt1NdrSIGoppLsiCywguzu7c+fe8/SPezacTmf2znzuPT++3/t6JU/me875fr7P5/z6Pufc+ZznAUZuPpBlxYaQ49CYOQgAAAAAAAAAAEZKcRAAAAAAAAAAAIyU4iAAAAAAAAAAABip2aYTAAAAuKLJgH7P0IrrYM+mpbB+eqcUt//460pxDz/pTCnuwS+r3b8kOf/kWtzujfu1wHMXS2GnTu/V4k7V8pxO56W4JJm2Xoprxbiqan99zevRH6W//WLs/unaefHi2dp7ce/62p+OHrrhVCnu4rn6eX/vulqu507XzotnZ7Vcp587X4pr2S3FZa94Tuzrfd8nqY/78/p5EQAA4CSZ95b5mv9+UjGEHIdmQH9pBwAAAAAAAAAAroXiIAAAAAAAAAAAGCnLigEAAAAAAAAAjJxlxU4uMwcBAAAAAAAAAMBIKQ4CAAAAAAAAAICRUhwEAAAAAAAAAAAjNdt0AgAAAAAAAAAArFZPMk/bdBqH6ptOYITMHAQAAAAAAAAAACOlOAgAAAAAAAAAAEbKsmIAAAAAAAAAACM37y3zvv3Lig0hx6ExcxAAAAAAAAAAAIyU4iAAAAAAAAAAABgpxUEAAAAAAAAAADBSs00nAAAAnACTgfwuoR1hLetp7T72nVOluP3HX1eKO//ks6W4Lzx1Wuyvl+KSZO/GvVLc7PqLpbid07W42XReiptOqnH1x7SqtfX3WbPePHuvnzOqf5DZL75uqq/TvVP7pbjdU7X+LhTPiUmyv1N7VKtxfVI7n9bO3sn0s7XXd+vF98V+7TlMklT7rKp+zpgf4T4CAAAM0Dwt8yP8PWNd5tn+HIdmIH+hBwAAAAAAAAAArpXiIAAAAAAAAAAAGCnFQQAAAAAAAAAAMFLVJe4BAAAAAAAAABiIeW+Z97bpNA41hByHxsxBAAAAAAAAAAAwUoqDAAAAAAAAAABgpCwrBgAAAAAAAAAwcpYVO7nMHAQAAAAAAAAAACOlOAgAAAAAAAAAAEZKcRAAAAAAAAAAAIzUbNMJAAAAAAAAAACwYr2l97bpLA43hBwHxsxBAAAAAAAAAAAwUoqDAAAAAAAAAABgpCwrBgAA8KhWn66275wqxc0fd7YU9/CTzpTivvDUaSnu/C29FLd308VSXJKcPnehFLezs1eKm03mpbjppPbYtFaLG5LJQO7jvDhV9Saew9m01mf1dTqtvi+mtbjd2X4pLkkuTE+X4s5PaufvpHY+bfPa+fu6vdpjOvmz4jnqkd1SXJKkD+O9DwAAcNLM0zLP9i/ZNYQch8bMQQAAAAAAAAAAMFKKgwAAAAAAAAAAYKQUBwEAAAAAAAAAwEjNNp0AAAAAAAAAAACrNe8t8942ncahhpDj0Jg5CAAAAAAAAAAARkpxEAAAAAAAAAAAjJRlxQAAAAAAAAAARq73lj6AJbuGkOPQmDkIAAAAAAAAAABGSnEQAAAAAAAAAACMlOIgAAAAAAAAAAAYqdmmEwAAAAAAAAAAYLXmPZn3tuk0DjXvm85gfMwcBAAAAAAAAAAAI6U4CAAAAAAAAAAARsqyYgAAwNWbDOT3Ba04Ne5sWu6yX3e6FHfhprOluAe/rJbr+SfX5uTdu+liKe70uQuluCQ5vbNXiptN90txk+LLprX1znM8WXN/J8GQHtPq1N/V1+m0eNpvrfY+bDu1/pIk52ph1bPU+fmpUtx0t3b+nl6ojRdnLhafi71aXJLkYvE91df8Xqx+rpnPjzcPAACANem9pQ9gWbEh5Dg0A/nLPgAAAAAAAAAAcK0UBwEAAAAAAAAAwEgpDgIAAAAAAAAAgJGabToBAAAAAAAAAABWq/eWeW+bTuNQfQA5Do2ZgwAAAAAAAAAAYKQUBwEAAAAAAAAAwEhZVgwAAAAAAAAAYOR6kt43ncXhBpDi4Jg5CAAAAAAAAAAARkpxEAAAAAAAAAAAjJTiIAAAAAAAAAAAGKnZphMAAAAAAAAAAGC15mmZp206jUMNIcehMXMQAAAAAAAAAACMlOIgAAAAAAAAAAAYKcuKAQDASTQZyO8EWnH62Gnt/vUzO7X+kuw97kwp7qEn1b6WnX9yKSx7N+6V4k6fu1CL26n1lySz6X4pbjrp5T4rJm29/XGyrfv1Ni/GTcuzf9fe90mS6in8XC3swn7tTp7frZ33Z+eLcQ/VxqdTu/Xzd5sXXzn7xbi+5vNw9XNU9XEBAAA4Jr239L79S3YNIcehGcj/CAAAAAAAAAAAANdKcRAAAAAAAAAAAIyU4iAAAAAAAAAAABip2mLlAAAAAAAAAAAMxry3zHvbdBqHGkKOQ2PmIAAAAAAAAAAAGCnFQQAAAAAAAAAAMFKKgwAAAAAAAAAAYKRmm04AAAAAAAAAAIDV6v2gbbsh5Dg0Zg4CAAAAAAAAAICRUhwEAAAAAAAAAAAjtZHioNbaW1trfal95VXEfG1r7Rdaa/e21i4s/v2F1trXriFlAAAAAAAAAIDh6i19AC29bfqRGp21Fwe11v5Kkjdcw/6ttfbOJO9L8ooktyTZWfz7iiTva629s7Xm1QEAAAAAAAAAAEvWWhzUWpsk+ddJZkk+dZVhP5DktYvt30vy95M8b/Hv7y2uf22S7z++TAEAAAAAAAAAYPjWPXPQP07y3CQfSfKuw3ZurT0jyRsXF38nyQt67+/uvd/de393khcurk+SN7XWnr6CnAEAAAAAAAAAYJBm6+qotfaUfHF2n3+Y5CuvIuwN+WKOr++9P7x8Y+/9fGvt9Un+22K/70zy+mNJGAAA2Lzq6sHTaSlsft1Orb8kj9xUi334SbX7uHvjfiludv3FUtzOzl6tv2ktzySZTno5tmLS1tsfDEH1fTHvtXPb0d73tfNNL57694vn090ba7/Ve/jB2th2+gu1Ozg9Xx8Tpxdqj03mxee/O38DAABcjd5bevE7+zoNIcehWefMQf8qybkkP9l7//XDdm6ttSTfsLj4kd77b15uv8X1f7i4+PJFHAAAAAAAAAAAnHhrKQ5qrX1jkpcl+WyS77rKsK9Icsti+4OH7Pvo7bcmedq15gcAAAAAAAAAAGO08mXFWmtfkuRHFhff1Hv/9FWG3ra0/ZFD9l2+/bYkH7vKPgAAAAAAAAAARm/eW3kp8HUaQo5Ds/LioCRvTXJzkv+a5F3XEPeUpe17D9n3nivEXZXW2q2H7HLztR4TAAAAAAAAAAA2baXFQa21FyZ5TZK9JN/Re+/XEH7D0vaDh+z70NL2uWvo41H3HL4LAAAAAAAAAAAMy2RVB26t7ST5sSQtyT/vvX/oGg9xZml795B9Lyxtn73GfgAAAAAAAAAAYJRWOXPQP01yW5L/k+T7CvGPLG3vHLLv6aXthwt9HbYU2c1J7i4cFwAAAAAAAABg43o/aNtuCDkOzUqKg1prz0zyTxYXX997f+ix9r+CB5a2D1sq7Pql7cOWIPv/9N7vfazbW2vXekgAAAAAAAAAANi4Vc0c9IYczPbzv5Nc11r7psvs85eXtv9Ga+3mxfYvL4qJlgt2bj2kv+WZf+651mQBAAAAAAAAAGCMVlUc9OgyX38+yc9exf7/bGn7K5I8lOTDS9c985D45dv/4Cr6AwAAAAAAAAA4MQ6WFdv+lZMsK3b8JptO4DF8LMknFtsvPmTfFy3+/ZMkH19VQgAAAAAAAAAAMCQrKQ7qvd/ee2+P1ZJ831LIVy3d9vHFMXqSuxa3P7O19vzL9bW4/tGZg+5axAEAAAAAAAAAwIm3zTMHJcnbk+wttt/RWju7fOPi8jsWF/cW+wMAAAAAAAAAAElmm07gsfTe/6i19rYkb07ynCS/0Vr7wSQfTfL0JG9K8uzF7j/Ue//jzWQKAAAbMtn2ev8jmtTWv+6nT5Xi9m44XYpLkodvrD0XjzyhOPnpuYulsJ3TtbjZZF6KKz6FRzJpJpSFTau+D+e9ftKonm+q57fq+fThc9XxohZXHZ9Of64+Jk7O75bi2v5+rcPaU7h+R/ncNh/KnQQAALZZT0s/wnfvdenZ/hyHZquLgxa+O8kTk7wqB4VA777MPu9K8j3rTAoAAAAAAAAAALbd1v/MuPc+772/OsnXJbkrySeS7C7+vSvJS3vvr+m9+/kMAAAAAAAAAAAs2djMQb33O5LccQ37vzfJe1eVDwAAAAAAAADAWPVF23ZDyHFotn7mIAAAAAAAAAAAoEZxEAAAAAAAAAAAjJTiIAAAAAAAAAAAGKnZphMAAAAAAAAAAGC1em/pvW06jUMNIcehMXMQAAAAAAAAAACMlOIgAAAAAAAAAAAYKcuKAQAAAAAAAACMXV+0bTeEHAfGzEEAAAAAAAAAADBSioMAAAAAAAAAADgRWmtPbK29rLX2ltba+1prn2mt9UW78xj7uWPpuIe1rzyufi/HsmIAAAAAAAAAAJwU9206gXVTHAQAsC0ma57UcT5fb3+cbK3V4orvi376VClu98/VvyJdeELtPu49bq8Ud+p0LW42rb33p5PaQt+t1RcInxwhluM178X38Mh5jR6/ozym1U820+JHsOr5dFY8f198XG2MuvCEWtxRxsRTn6+Nw+2R3VqHrfjsd+9h1mjd3/cAAOByeksfwt951pvjPUn+IMnXrLifZx1y+8dW2bniIAAAAAAAAAAAToq3JLk7yd299/taa0/Liotzeu//c5XHP4ziIAAAAAAAAAAAToTe+/duOod1UxwEAAAAAAAAADByvQ9jheUh5Dg0FjoGAAAAAAAAAICRUhwEAAAAAAAAAAAr0lp7f2vt/tbabmvtU621X2+tvbm19vh19G9ZMQAAAAAAAAAAttHNrbXH3KH3fu+acjmKlyxtf2mSFy/am1prt/fe71pl54qDAAAAAAAAAABGrveW3h+70GYbXJLj3VcRss136kNJfinJbyf5RJJTSf5ikm9J8jVJviTJz7fWvr73/r5VJaE4CAAAAAAAAAAAjtfbe+93XOb630ryb1trr0vyo0mmSX68tfaM3vvDq0hEcRAAAAAAAAAAANvouUk+uekkKnrvnz/k9ne21p6T5DVJvizJK5P89CpyURwEAAAAAAAAAMA2+mTv/d5NJ7FC78xBcVCSvDiKgwAAAAAAAAAAKOntoG27IeR4fD68tH3LqjqZrOrAAAAAAAAAAADAFa2lEkpxEAAAAAAAAAAArN9fWtr+xKo6sawYAMBJNTkBdeLz+aYz4FGt9uOHPq29TvfPnirF7d5Qf1/s3tBrgWf3S2GnTtXippPa+6K14v3j2M1P1rTKW20Tz8XEe/HYVc9v1fNp9fx9sThe7N4wLcbVx8QzxXG4PVjrs+0V34vd+2lrnITvJgAAsAV6H8ZXoSHkeIxet7T9wVV14lsXAAAAAAAAAABcpdba7a21vmh3XOb2Z7XWnnHIMV6X5NWLi59M8ovHn+kBMwcBAAAAAAAAAHAitNZemGS5cOempe1ntNZuX96/935noZu/luTHW2sfSPK+JB9Kcn8O6nSemeRbk/zNxb77SV7Xe3+o0M9VURwEAAAAAAAAAMBJ8Zok33aF216waMvuLPYzTfKSRbuS+5O8uvf+nmIfV0VxEAAAAAAAAADA2PVF23ZDyPFw783BkmF/PcmzkzwpyY1JWpLPJvnvSX41yZ299y+sOhnFQQAAAAAAAAAAnAi999uT3H7EY9yZx5hRqPf+qST/ZtE2brLpBAAAAAAAAAAAgNUwcxAAAAAAAAAAwMj13tJ723QahxpCjkNj5iAAAAAAAAAAABgpxUEAAAAAAAAAADBSioMAAAAAAAAAAGCkZptOAAAAAAAAAACANeibToBNMHMQAAAAAAAAAACMlOIgAAAAAAAAAAAYKcuKAQAwXpNiLfx8frx5HKaa55BMWi1uOi2Fzc/U4nbPFfNMsn9dbT7e6en9Wty09jqdTtY7b/Ckmaf4Sua9/nrj5Kq+bk7Ce7F6H6uPafV8Wj5/F8eL6vh0lDGxOg5Xx/1M9mpxa/7IdyQ+1wIAAMeg95Y+gL9JDSHHofFtDQAAAAAAAAAARkpxEAAAAAAAAAAAjJTiIAAAAAAAAAAAGKnZphMAAAAAAAAAAGDF+qJtuyHkODBmDgIAAAAAAAAAgJFSHAQAAAAAAAAAACNlWTEAAAAAAAAAgNFri7bthpDjsJg5CAAAAAAAAAAARkpxEAAAAAAAAAAAjJTiIAAAAAAAAAAAGKnZphMAAAAAAAAAAGDF+qJtuyHkODBmDgIAAAAAAAAAgJFSHAQAAAAAAAAAACNlWTEAALjURA39FbW23v5m01LY3pli3HX1+zc/My/FnZrW4qbN3LoAQ1Q9f0+L48XFM7X+9q6rfx6qjsM7xXG/rPq5pg9oDPa5FgAAWGZZsRPLt0MAAAAAAAAAABgpxUEAAAAAAAAAADBSioMAAAAAAAAAAGCkZptOAAAAAAAAAACAFetJett0Fofrm05gfMwcBAAAAAAAAAAAI6U4CAAAAAAAAAAARkpxEAAAAAAAAAAAjNRs0wkAAAAAAAAAALBavR+0bTeEHIfGzEEAAAAAAAAAADBSioMAAAAAAAAAAGCkLCsGAAAAAAAAADB2fdG23RByHBgzBwEAAAAAAAAAwEgpDgIAAAAAAAAAgJFSHAQAAAAAAAAAACM123QCAADACdBaKazPar9nmO/U+pvvlMKSJH02L8VNJrUFtFtbbxxXNu+11xusU/V1OnHOuKJ1n4er40V1fJrvTEtxB7HrHfdb8XMGAADAidPbQdt2Q8hxYMwcBAAAAAAAAAAAI6U4CAAAAAAAAAAARsqyYgAAAAAAAAAAI9eSDGEFcYuKHT8zBwEAAAAAAAAAwEgpDgIAAAAAAAAAgJFSHAQAAAAAAAAAACM123QCAAAAAAAAAACsWF+0bTeEHAfGzEEAAAAAAAAAADBSioMAAAAAAAAAAGCkLCsGAAAAAAAAADB2vR20bTeEHAfGzEEAAAAAAAAAADBSioMAAAAAAAAAAGCkFAcBAAAAAAAAAMBIzTadAAAAcAK04hrRxbj5rBg3LYUdmPZSWGu1uHWbDCRPgGtVPb/Ne3FsKyqPF8Xx6ShjYnUcXvfnBQAAgBOnL9q2G0KOA2PmIAAAAAAAAAAAGCnFQQAAAAAAAAAAMFKWFQMAAAAAAAAAGDvLip1YZg4CAAAAAAAAAICRUhwEAAAAAAAAAAAjpTgIAAAAAAAAAABGarbpBAAAAAAAAAAAWLG+aNtuCDkOjJmDAAAAAAAAAABgpBQHAQAAAAAAAADASFlWDAAAAAAAAABg7Ho7aNtuCDkOjJmDAAAAAAAAAABgpBQHAQAAAAAAAADASCkOAgAAAAAAAACAkZptOgEAAIAr6a24tnR1Seqj/Hyimmrra40DYLPWft4f0JhYHvcBAAC4Kq0ftG03hByHxsxBAAAAAAAAAAAwUoqDAAAAAAAAAABgpCwrBgAAAAAAAAAwdn3Rtt0QchwYMwcBAAAAAAAAAMBIKQ4CAAAAAAAAAICRUhwEAAAAAAAAAAAjpTgIAAAAAAAAAABGSnEQAAAAAAAAAACMlOIgAAAAAAAAAAAYKcVBAAAAAAAAAAAwUrNNJwAAAAAAAAAAwGq1JK1vOovDtU0nMEJmDgIAAAAAAAAAgJFSHAQAAAAAAAAAACNlWTEAAGBrtV6c47Y6Ne68GHeEPnuvTZJbjas/OAAch7Wf9wc0JpbHfQAAAK5Obwdt2w0hx4ExcxAAAAAAAAAAAIyU4iAAAAAAAAAAABgpxUEAAAAAAAAAADBSs00nAAAAAAAAAADAivVF23ZDyHFgzBwEAAAAAAAAAAAjpTgIAAAAAAAAAABGyrJiAAAAAAAAAABjZ1mxE8vMQQAAAAAAAAAAMFKKgwAAAAAAAAAAYKQUBwEAAAAAAAAAwEjNVnXg1trjkrw0yXOTPCfJLUm+NMnZJJ9P8uEk703yrt77/VdxvK9N8tokz1sc59NJfjvJj/Xef3UV9wEAAAAAAAAAYAxaP2jbbgg5Ds3KioNyUMTzs1e47UuTvHjRvqu19q2991+73I6ttZbkR3NQGLTsliSvSPKK1tqPJfmO3ruXCAAAAAAAAAAALKyyOChJ7knygSS/u9j+0xwsZXZrkr+b5JVJbkryntbac3vv/+Myx/iBfLEw6PeSvDXJR5M8Pckbkzx7cfunk3zPyu4JAAAAAAAAAAAMzCqLgz7Qe3/qY9z+c621lyf5xSQ7Sb43yd9Z3qG19owcFAAlye8keVHv/eHF5btba+9J8sEcLFv2ptbaT/TeP3qcdwIAADgG1Uk+i3GTvWLcfinswH4rhfVei1u3eTHPiTmAgS1XPb+tW3m8KI5PRxkTq+Pwuj8vAAAAnDh90bbdEHIcmMmqDtx7P/RPCL33X0rykcXFF11mlzfkiwVMr18qDHo0/nyS1y8uzpJ8Zy1bAAAAAAAAAAAYn5UVB12Dhxb/nlm+srXWknzD4uJHeu+/ebngxfV/uLj48kUcAAAAAAAAAACceBstDmqt3Zbkry4ufuSSm78iyS2L7Q8ecqhHb781ydOOJTkAAAAAAAAAABi42eG7HK/W2nU5KPr5+iRvTDJd3PQjl+x629L2pYVDl1q+/bYkHztKjgAAAAAAAAAAo9IXbdsNIceBWUtxUGvt9iQ/8Ri7vC3JT19y3VOWtu89pIt7rhB3VVprtx6yy83XekwAAAAAAAAAANi0tc8cdInfT/IdvfffusxtNyxtP3jIXI0vvwAAIABJREFUcR5a2j5XyOOew3cBAAAAAAAAAIBhWVdx0C8l+Z3F9tkkT0/yjUlekeSnW2vf2Xv/lUtizixt7x5y/AtL22ePkigAAAAAAAAAwNi0ftC23RByHJq1FAf13j+f5PNLV92d5N2ttX+Q5CeT3NVae3Xv/c6lfR5Z2t45pIvTS9sPF1I8bCmym3OQMwAAAAAAAAAADMZGlxXrvf9Ua+1lOZhF6F+01u7qvX9ucfMDS7setlTY9Uvbhy1Bdrk87n2s21tr13pIAAAAAAAAAADYuMmmE0hy1+Lf65P8raXrlwt2bj3kGMsz/9xzHEkBAAAAAAAAAMDQbXTmoIVPL21/+dL2h5e2n3nIMZZv/4MjZwQAAAAAAAAAMCa9HbRtN4QcB2YbZg66ZWl7eUmwjyX5xGL7xYcc40WLf/8kycePJy0AAAAAAAAAABi2bSgO+ntL2x96dKP33vPFJcee2Vp7/uWCF9c/OnPQXYs4AAAAAAAAAAA48Va2rFhr7fYk7+69P/IY+7whyUsXFz+e5L9cssvbk3x7DvJ8R2vtRb33h5fizyZ5x+Li3mJ/AABg2xRr+NvevBQ32a31N9kthSVJ2l7ttxfzeW2K3F6cWrca15rfYVzJpPjYzE2PzBpVX6dc2brPw9Xxojo+HWVMrI7D1XG/+jkDAADgRPIV6kRaWXFQkjuS/HBr7edzUPTz0RwsG3ZDkmcl+ZYkL1jsu5vk23vve8sH6L3/UWvtbUnenOQ5SX6jtfaDi2M9Pcmbkjx7sfsP9d7/eIX3BwAAAAAAAAAABmWVxUFJ8oQczPzz7Y+xz71JXtV7/49XuP27kzwxyatyUAj07svs864k33OEPAEAAAAAAAAAYHRWWRz01UlekuSrktyW5ElJbkzySJL7kvx+kl9J8nO99/NXOkjvfZ7k1YsZiF6b5LlJbkrymSR3J3ln7/19K7wfAAAAAAAAAAAwSCsrDuq9fzQHy3+985iO994k7z2OYwEAAAAAAAAAnCStH7RtN4Qch2ay6QQAAAAAAAAAAIDVUBwEAAAAAAAAAAAjtbJlxQAAAAAAAAAA2BJ90bbdEHIcGDMHAQAAAAAAAADASCkOAgAAAAAAAACAkVIcBAAAAAAAAAAAIzXbdAIAAAAAAAAAAKxYT1rfdBJXYQg5DoyZgwAAAAAAAAAAYKQUBwEAAAAAAAAAwEhZVgwAAC41n6+3v8mAavb7mudz3dsvhc0eKcadr39FmjzSSnH7+7Xnf7/X+vMlEGCzqufv6nhRHZ9m5+tjfnUcro77Zev+XLMJPtcCAAAQfxcGAAAAAAAAABi/vmjbbgg5DoyfcgAAAAAAAAAAwEgpDgIAAAAAAAAAgJGyrBgAAAAAAAAAwNhZVuzEMnMQAAAAAAAAAACMlOIgAAAAAAAAAAAYKcVBAAAAAAAAAAAwUrNNJwAAAAAAAAAAwGq1ftC23RByHBozBwEAAAAAAAAAwEgpDgIAAAAAAAAAgJFSHAQAAAAAAAAAACOlOAgAAAAAAAAAAEZKcRAAAAAAAAAAAIyU4iAAAAAAAAAAABip2aYTAACAlZnPN53B1TlKnpOB1PvPey1uf78UNnmkFrfzYDHPJNPzrRS3e2Faits/XXvu9ye119tsWnts5r32uCTJpNWfjyGo3r+jPKYM39jfF0ex7vfG/rzW3/5+8fxdHC92iuPTUcbE6jhcHffLnzOGZOyfa4fymRYAAIauL9q2G0KOA+NbFwAAAAAAAAAAjJTiIAAAAAAAAAAAGCnLigEAAAAAAAAAjFzrB23bDSHHoTFzEAAAAAAAAAAAjJTiIAAAAAAAAAAAGCnFQQAAAAAAAAAAMFKzTScAAAAAAAAAAMAa9E0nwCaYOQgAAAAAAAAAAEZKcRAAAAAAAAAAAIyUZcUAAAAAAAAAAMauZxjLig0hx4ExcxAAAAAAAAAAAIyU4iAAAAAAAAAAABgpxUEAAAAAAAAAADBSs00nAADAhsznm86Ak6TXFolu+7XX6fThi6W4nQdOl+IOYmtfr3YfnpbiLp6txc2mxcd0UnwOmwXCj9tkA4/pvLe19zkEm3guOH69+Pren9d+c3fxYu38neJ4sfNA7f7tPLBfikvq43B13K9+zmCLnITvJhO/0wUAYPNaP2jbbgg5Do1vJAAAAAAAAAAAMFKKgwAAAAAAAAAAYKQsKwYAAAAAAAAAMHZ90bbdEHIcGDMHAQAAAAAAAADASCkOAgAAAAAAAACAkVIcBAAAAAAAAAAAIzXbdAIAAAAAAAAAAKxW6wdt2w0hx6ExcxAAAAAAAAAAAIyU4iAAAAAAAAAAAE6E1toTW2sva629pbX2vtbaZ1prfdHuXFGf39Ra+7XW2p+21h5prX28tfZTrbXnr6K/S1lWDAAAAAAAAABg7PqibbvV53jfyntYaK2dSfLvkrzskpu+fNG+ubV2R+/9+1eZh5mDAAAAAAAAAAA4ie5J8h9WePx35YuFQR9I8vIkz0vy6iQfzUHdzltaa69ZYQ5mDgIAAAAAAAAA4MR4S5K7k9zde7+vtfa0JB877k5aay9O8s2Li7+c5BW99/3F5btba+9J8rtJnprkra21f997//xx55GYOQgAAAAAAAAAgBOi9/69vfdf6b2venmxNy7+3U/yj5YKgx7N4zNJ3rS4+PgczCa0EmYOAgDYFvP5pjOA1enFRaKL74t24WIpbufP9kpxSXL6s9NS3CNPqMXtXV/7Ord3av/wnS5jOqk9F9Mj/CSlelactCEsnD4sHlPWZd5bObYXY/fntbi9/doJbu9C7fw9+0JtvDj92dr79yhjYnUcLn8ern7OgHXyfQ8AgG3QF23bDSHHQ7TWziX56sXF9/fe773Crr+Q5AtJHpfklUl+eBX5mDkIAAAAAAAAAACOz/OSnF5sf/BKO/Xed5P85qMxrbVTq0hGcRAAAAAAAAAAAByf25a2P3LIvo/ePkvyF1aRjGXFAAAAAAAAAADYRje39thLgj/Gkl2b9JSl7cPyu+eSuA8fdzKKgwAAAAAAAAAARq71g7btLsnx7qsJWU0mR3LD0vaDh+z70NL2uRXkYlkxAAAAAAAAAAA4RmeWtncP2ffC0vbZFeRi5iAAAAAAAAAAALbSc5N8ctNJFDyytL1zyL6nl7YfXkEuioMAAAAAAAAAAEavL9q2+39z/GTv/d4NZXIUDyxtH7ZU2PVL24ctQVZiWTEAAAAAAAAAADg+ywVNtx6y71OWtu9ZQS6KgwAAAAAAAAAA4Bh9eGn7mYfs++jte0n+1yqSURwEAAAAAAAAAADH5+4ku4vtF19pp9baTpLnPxrTe9+90r5HoTgIAAAAAAAAAOAk6ANoI9B7fyDJf1pcfElr7UpLi70yyeMW27+4qnwUBwEAAAAAAAAAwFVqrd3eWuuLdscVdnvb4t9Zkn/ZWptecoybkvzg4uLnk/z4SpJdJAAAAAAAAAAAAKPXWnthkmcsXXXT0vYzWmu3L+/fe7+z0k/v/T+31t6d5JuS/O0k72+tvT3JJ5I8K8l3J3nqYvc3994/V+nnaigOAgCAIZvPa3GTgUwiOq/NIdsuXCzFzR64UIpLkrP3nyrFXXjC9PCdLuOhG2r97Z6qvWam01pca/uluCSZtlrcvNcCJ20kcxbDFqi+D4/WZy1ub14bE3cv1M7DebAWd+aztcf07P218/BRxsTqOFx+Eoei+rkNAADgmLR+0LbdGnJ8TZJvu8JtL1i0ZXceoa9X5WDZsJcm+apFWzZP8v2993ceoY9DDeR/BAAAAAAAAAAAYDh67w/33r8uybckeX+STyXZTXJPkp9J8sLe+x2rzsPMQQAAAAAAAAAAnAi999uT3H7EY9yZa5hRqPf+MzkoBtoIMwcBAAAAAAAAAMBImTkIAAAAAAAAAGDs+qJtuyHkODBmDgIAAAAAAAAAgJFSHAQAAAAAAAAAACNlWTEAAAAAAAAAgJFr/aBtuyHkODRmDgIAAAAAAAAAgJFSHAQAAAAAAAAAACOlOAgAAAAAAAAAAEZqtukEAAAAAAAAAABYsb5o224IOQ6MmYMAAAAAAAAAAGCkFAcBAAAAAAAAAMBIWVYMAADYXr04f+z+filscn631l+SM5+pxZ593JlS3MVztd96XNg5VYrbndUe07ZTCluo9Tmd1F43895KcZNmnmPGq/q+qNqf1/vb25+W4nZ3a38e23uodj49fX/t/H32vtq5pjo+HWVMrI7D5XEfAACAq2NZsRPLzEEAAAAAAAAAADBSioMAAAAAAAAAAGCkFAcBAAAAAAAAAMBI1RZVBwAAAAAAAABgMNqibbsh5Dg0Zg4CAAAAAAAAAICRUhwEAAAAAAAAAAAjZVkxAAAAAAAAAICx64u27YaQ48CYOQgAAAAAAAAAAEZKcRAAAAAAAAAAAIyU4iAAAAAAAAAAABip2aYTAAAAAAAAAABgxXrS+qaTuApDyHFgzBwEAAAAAAAAAAAjpTgIAAAAAAAAAABGyrJiAABwEs3ntbjJmn9f0Ivzx+7X7l97ZLfWX5LZFx4pxV1/X+1r2d51tbj9nVrchenpUlzO1cKSJDvVwP1S1KTVeiu+m8omg5j7mVWZ9+ILtagX+5sXX6Z7+9NaYJILu8Xz24O189vs/lp/1/1pKSzX37dXiquOT0cZE6vjcHncX7fq5ygAAIBN6xnGkl1DyHFgzBwEAAAAAAAAAP+XvbuPlfS66wP+/d293hfbITUEahGHQjGNQwsIiCPeEqNCIxRTkSBUEtpCICbQSikhIJICUlGgtLxEDaK0vKWxUlARlARDE1qoEiyIEuKgqKWKQ14I1MZJqPNCbK93r/fO6R/zrHay7O7ce+6dO/M8+/lIR3OemfN7zm/unZe9s785B2CiFAcBAAAAAAAAAMBEKQ4CAAAAAAAAAICJ6tscHQAAAAAAAACAcWnrToB1sHIQAAAAAAAAAABMlOIgAAAAAAAAAACYKMVBAAAAAAAAAAAwUdvrTgAAAAAAAAAAgNWqNm+bbgw5jo2VgwAAAAAAAAAAYKIUBwEAAAAAAAAAwETZVgwAAAAAAAAAYOra0DbdGHIcGSsHAQAAAAAAAADARCkOAgAAAAAAAACAiVIcBAAAAAAAAAAAE7W97gQAAIARmc364raO+HsJrXNT6nO73VPW6bNdcScePNYVd/2Ja7vido/3/Rl4euuarri+n8rg+r6wdrwvbnur7/F9rPPhXdX3OJ216ptwDbY67+NRG9PPtFfrvI+7s764c7O+J8bOTv9HVWcfPtEVt/1g3+vbtR/o+9lc/8C5rrgTDz7aFdf7/nSQ98Tu9+Gj1vvvGgAAgJGqNm+bbgw5jo2VgwAAAAAAAAAAYKIUBwEAAAAAAAAAwETZVgwAAAAAAAAAYOra0DbdGHIcGSsHAQAAAAAAAADARCkOAgAAAAAAAACAiVIcBAAAAAAAAAAAE7W97gQAAAAAAAAAAFitavO26caQ49hYOQgAAAAAAAAAACZKcRAAAAAAAAAAAEzUSouDquqLqur7q+q3q+q+qjpbVQ9X1bur6s6qevo+z/c1VfXaqrp/ONf9w/HXrOo+AAAAAAAAAACMXhtR41Btr+rEVXV3kmdc4qbjST5naN9SVf85yR2ttZ0rnKuS/GySF1500xOTPCfJc6rq55N8Z2vNwwQAAAAAAAAAALLalYOeOFw+kOSnknxDkqcl+dIkL0nyF8Pt/zTJnUvO9SO5UBj0jiTPG871vOE4w+0/fAh5AwAAAAAAAADAJKxs5aAk70ry/Ul+vbW2e9Ftbx1WDHpzkr+T5HlV9R9ba79/8Umq6uYk3zccvj3JM1prjw7H91TVbya5O8lTk7y0ql7dWnvfCu4PAAAAAAAAAACMysqKg1prX7vk9ger6nuS/NZw1Tck+WvFQUm+OxfyfNFCYdD585yuqhclecsw7sVJXnSQ3AEAgKvUAXYprp3HuuK2Pv7o8kGXcGq7byHYtnWqKy451hV1enZN53zJ2d3qitu9ru93cfxEX9z2sVlX3LGt3rij3027qm/OWev7HY5FW8P92511Pi9mfa8Z53b74nbO9j33zz3S/5qx/eG+j7mu/UDfz/ST/u/F34Xbm1MfOtMV1/t+0fv+dJD3RAAAADZUG9qmG0OOI7PKbcX24vcW+p998Y1VVUm+bjh8V2vtrZc6yXD9nwyHzx7iAAAAAAAAAADgqrbu4qDjC/1LfWXys5I8cejfveRc52+/KclnHiwtAAAAAAAAAAAYv5VtK7ZHty3033WJ25+y5PZc5vanJHl/b1IAAAAAAAAAAFNSbd423RhyHJu1FQdV1VaSly1c9auXGPakhf79S05532Xi9pLLTUuG3Lif8wEAAAAAAAAAwCZY58pB353kaUP/da21t19izOMW+g8vOd8jC/3r95nLfcuHAAAAAAAAAADAuGytY9Kqui3Jvx0O/zLJP7vM0JML/Z0lpz270D/VmRoAAAAAAAAAAEzGka8cVFV/N8nrhrnPJvlHrbUPXWb4mYX+8SWnPrHQf3SfaS3bhuzGJPfs85wAAAAAAAAAAJujrTsB1uFIi4Oq6rOS/E6SG5LsJnlea+3uK4Q8tNBftlXYdQv9ZVuQfYLW2v1Xur2q9nM6AAAAAAAAAADYCEe2rVhVfXqS/5nk0zOvRfu21trrloQtFu3ctGTs4uo/9+0/QwAAAAAAAAAAmJYjWTmoqp6Q5HeT/O3hqhe11l6zh9B3LvRvWTJ28fZ795EeAAAAAAAAAMCkVWuptvn7io0hx7FZ+cpBVfX4JP8jyecOV72stfYzewx/f5IHhv5tS8Y+Y7j8iyR/tp8cAQAAAAAAAABgilZaHFRV1yZ5fZIvGq761621H9trfGutJblrOLylqr7kMvN8SS6sHHTXEAcAAAAAAAAAAFe1lRUHVdXxJK9L8uXDVT/VWvvBjlO9Msm5of/TVXXqonlOJfnp4fDcMB4AAAAAAAAAAK562ys8939J8syh/8Ykr6qqv3eF8TuttXdffGVr7d1V9ZNJXpbkqUneXFU/luR9ST47yUuTfOEw/Cdaa+85rDsAAAAcktmsL25r5Tshf6KDLEK623cfa+exrrhjHz3dFXdtV1RSs5Ndccd2jnXOmJze6fuTdedT+h43j17fF7d94tzyQZdwzTW7XXHHjnU+n5Icq77HeHXG9eqdr7U65ExWN99uZ+zubt/j9LHH+p6L5852fnT08DVdYSc+3P+6f+0H+uKuf6DvOXzqQ2e64npfv3vfL3rfnw70nnjUev+dAQAAcLVpQ9t0Y8hxZFZZHPT1C/2/n+R/Lxn/50k+8zK3/UCST0vybZkXAv3KJca8KknPykQAAAAAAAAAADBJR/w13D6ttVlr7QVJbk9yV5IHkuwMl3cleVZr7Y7Wmq8JAQAAAAAAAADAYGUrB7UVrOXdWntDkjcc9nkBAAAAAAAAAKas2rxtujHkODajWDkIAAAAAAAAAADYP8VBAAAAAAAAAAAwUYqDAAAAAAAAAABgorbXnQAAAAAAAAAAACvWhrbpxpDjyFg5CAAAAAAAAAAAJkpxEAAAAAAAAAAATJTiIAAAAAAAAAAAmKjtdScAAAAAAAAAAMBqVZu3TTeGHMfGykEAAAAAAAAAADBRioMAAAAAAAAAAGCibCsGAABsrtmsL25rDd+DaJ1r3Z7b7Qqr7HTFHftIX57Xnuv7XRw7e6orLkm2T/f9yfrow8e64s58ct/j5rFP6svzsVN9v/tjJ/rikuTYsb7f49ZW3+OmOteA7o1rrY40bjbri0uS3d2+x9vu2b7Hdx7ti9v+eF/cyY/0/WxOfah/3fDrPnSuK+7Eg492xW19vC+udh7riut9v+h+f1qH3vd9AAAA9qYNbdONIceRsXIQAAAAAAAAAABMlOIgAAAAAAAAAACYKMVBAAAAAAAAAAAwUdvrTgAAAAAAAAAAgNWqNm+bbgw5jo2VgwAAAAAAAAAAYKIUBwEAAAAAAAAAwETZVgwAAAAAAAAAYOra0DbdGHIcGSsHAQAAAAAAAADARCkOAgAAAAAAAACAiVIcBAAAAAAAAAAAE7W97gQAAAAAAAAAAFi9auvOgHWwchAAAAAAAAAAAEyU4iAAAAAAAAAAAJgo24oBAACsU+tcx/fcbldYdc639Vd9cScf68szSbYfOdkVd+Ljx7viHv2Uvu/PnP3kvj+tdx53rCtu99r+tZ8fO9kX27ZnfRMe68y1+sLS+6PZ7ZuwzvV/52rrTN+cx093xj3UF3fiI30/1FMf7nvun3xwpysuSbY/fqYrrk6f7YvbeawrLrudz6fe9wsAAAA4r7Vx/H05hhxHxspBAAAAAAAAAAAwUYqDAAAAAAAAAABgohQHAQAAAAAAAADARG2vOwEAAAAAAAAAAFar2rxtujHkODZWDgIAAAAAAAAAgIlSHAQAAAAAAAAAABNlWzEAAAAAAAAAgKlrQ9t0Y8hxZKwcBAAAAAAAAAAAE6U4CAAAAAAAAAAAJkpxEAAAAAAAAAAATNT2uhMAAAAAAAAAAGC1ajZvm24MOY6NlYMAAAAAAAAAAGCiFAcBAAAAAAAAAMBE2VYMAACYnlnnurNbI/r+RGt9cbt9P5s6s9MXd263Ky5Jrtk51xV37PTxrrgTHz3RFbfz+L4/rXce1/d427m+uuKS5Ny1fXPOjh/ri+sL6/8qU+9Tv/NhutX3tEiSbJ/uew4ff7gz7qG+O3n8r/qeh9sPne2K2zrd/0PtfZ1K7+tU7+twb9yY9L4PAwAAsFptaJtuDDmOzIg++QYAAAAAAAAAAPZDcRAAAAAAAAAAAEyU4iAAAAAAAAAAAJio7XUnAAAAAAAAAADAalWbt003hhzHxspBAAAAAAAAAAAwUYqDAAAAAAAAAABgomwrBgAAAAAAAAAwda3N26YbQ44jY+UgAAAAAAAAAACYKMVBAAAAAAAAAAAwUYqDAAAAAAAAAABgorbXnQAAAAAAAAAAAKtVbd423RhyHBsrBwEAAAAAAAAAwEQpDgIAAAAAAAAAgIlSHAQAAAAAAAAAABO1ve4EAAAANsZs1h+7NZLvXrTODbt74x7r3yC8On8fx84+1hW3dXqnK+6aj13TFXfyVF/c7OSxrrgkOdcZOztefXHbfXHpDEvnw23rXF/g1k7/43v7zG7fnJ1xxx7te15U5/OpNy67ffdvHtv5Gt77+jZ1B3lPBAAAYHP5M/iqNJJPrwEAAAAAAAAAgP1SHAQAAAAAAAAAABNlWzEAAAAAAAAAgImrNm+bbgw5jo2VgwAAAAAAAAAAYKIUBwEAAAAAAAAAwEQpDgIAAAAAAAAAgInaXncCAAAAAAAAAACsWGvztunGkOPIWDkIAAAAAAAAAAAmSnEQAAAAAAAAAABMlG3FAAAAAAAAAAAmrtq8bbox5Dg2Vg4CAAAAAAAAAICJUhwEAAAAAAAAAAATpTgIAAAAAAAAAAAmanvdCQAAAEzCbHa0822N5Lse7QAbhO92/kxnfXPW7m5f3JmdvriHO3+Hx471xSU5vt0X27Y7c63qm68zrnofb51xde4Az/tzfY+39D5Ou59PR/s8PNBrxkFix+Co32cAAACYnja0TTeGHEdmJJ8mAwAAAAAAAAAA+6U4CAAAAAAAAAAAJsq2YgAAAAAAAAAAE1dt3jbdGHIcGysHAQAAAAAAAADARCkOAgAAAAAAAACAiVIcBAAAAAAAAAAAE7W97gQAAAAAAAAAAFixWZu3TTeGHEfGykEAAAAAAAAAADBRioMAAAAAAAAAAGCibCsGAAAAAAAAADB1bWibbgw5joyVgwAAAAAAAAAAYKIUBwEAAAAAAAAAwEQpDgIAAAAAAAAAgInaXncCAAAAdJjN+uK2RvQdkda5uXhvXOePNNUXWOeqb76tc31xB1DVmWtv3FE76sfaQcyOONd13Mep6339BgAAgAOqJDWCP/VH8onSqIzoU2EAAAAAAAAAAGA/FAcBAAAAAAAAAMBE2VYMAAAAAAAAAGDqWsaxhfgIUhwbKwcBAAAAAAAAAMBEKQ4CAAAAAAAAAICJUhwEAAAAAAAAAMBVp6o+o6p+sqrurapHquojVfW2qvreqrr2gOf+oapqe2xfeUh36ZK2V3lyAAAAAAAAAAA2QEuqrTuJPTiiHKvq9iS/nOTxC1dfm+TWod1RVc9qrf3p0WS0OoqDAAAAAAAAAAC4alTVFyT51cyLgR5O8m+SvCnJqSTPTfLtSZ6c5PVVdWtr7eEDTvl5S25//wHPf0WKgwAAAAAAAAAAuJq8MvPCoHNJntlae8vCbW+sqvck+fEktyR5SZKXH2Sy1tr/OUj8QW2tc3IAAAAAAAAAAI5AG1Fboaq6NclXDoevuqgw6LxXJLl36L+4qq5ZbVarpTgIAAAAAAAAAICrxbMX+q++1IDW2izJa4bDG3KhmGiUFAcBAAAAAAAAAHC1ePpw+UiSP7rCuLsX+l+xunRWT3EQAAAAAAAAAABXi6cMl+9trZ27wrh3XSKmS1X9blV9uKp2quovq+r3quplVXXDQc67V9tHMQkAAAAbYjY72vm2roLvpLTOTdB74w7yK6w6QDCHqvf3z+E76tdFAAAAWJNqLTWCzyQuyvHGWvKZVmvt/j2fu+pkkicMh1eMa619tKoeSXJdkiftdY7L+OqF/qcmuW1oL62q57fW7jrg+a9IcRAAAAAAAAAAAJvonj2M2c834h630H94D+PPFwddv485Fv1xkt9I8rYkDyS5JsmTk/zjJM9M8jeS/HpV/cPW2m93zrGU4iAAAAAAAAAAAK4GJxf6O3sYf3a4PNUx1ytbaz90iev/MMlrquo7kvxskmNJfrGqbm6tPdoxz1KKgwAAAAAAAAAA2ES3JvngIZ7vzEL/+B7Gnxgu912001r72JLbf66qnprkjiSfnuTrk/zyfufZC8VBAABtsuvPAAAgAElEQVQAAAAAAABTNxvapvvEHD/YWrv/EM/+0EJ/L1uFXTdc7mULsh4/l3lxUJLclhUVB22t4qQAAAAAAAAAALBJWmtnkjw4HN50pbFVdUMuFAfdt6KU3rnQf+KK5lAcBAAAAAAAAADAVePe4fLmqrrSjlu3XCLmsNWKzvsJbCsGAAAAAAAAADBx1VqqtXWnsdQR5PgHSZ6e+apAX5zkDy8z7raF/ptXlMvnLvQfWNEcVg4CAAAAAAAAAOCq8RsL/W+91ICq2kryzcPhx5K8aUW5fMdC/+4VzaE4CAAAAAAAAACAq0Nr7W1Jfn84fEFVfeklhn1PkqcM/Z9qrT22eGNVPb+q2tB+6OLgqvq8qrr5SnlU1XckecFw+MEkr9vH3dgX24oBAAAAAAAAAHA1+a7Mtwo7leR3qupHM18d6FSS5yZ54TDu3Ule0XH+L07yi1X1piS/neSPk3w48zqdW5L8kyT/YBi7m+Q7WmuP9N2V5RQHAQAAAAAAAABMXRvapjuCHFtr76iqb0zyS0k+KcmPXmLYu5Pc3lp7qHOaY0m+emiX8+EkL2it/WbnHHuiOAgAAAAAAAAAgKtKa+23qurzM19F6PYkNyXZSfLeJL+W5N+31k53nv4NmW8Z9qVJvjDJ30zyKUkqyUeS/K8k/z3Jna21jx/kfuyF4iAAAAAAAAAAAK46rbU/T/KSoe0n7s4kd17h9r9M8p+GtnaKgwAAAFid2aw/dmvr8PJgro1h3WjodJDXGwAAALgqtJF8PjSGHMfFJ60AAAAAAAAAADBRioMAAAAAAAAAAGCiFAcBAAAAAAAAAMBEba87AQAAAAAAAAAAVqvavG26MeQ4NlYOAgAAAAAAAACAiVIcBAAAAAAAAAAAE2VbMQAAAAAAAACAqWtt3jbdGHIcGSsHAQAAAAAAAADARCkOAgAAAAAAAACAiVIcBAAAAAAAAAAAE7XS4qCq+rSq+tqqenlV/XZVPVhVbWh3dpzva6rqtVV1f1WdHS5fW1Vfs4L0AQAAAAAAAAAmoWbjaRyu7RWf/0OHcZKqqiQ/m+SFF930xCTPSfKcqvr5JN/ZWmuHMScAAAAAAAAAAIzdUW4rdl+S3+mM/ZFcKAx6R5LnJXnacPmO4foXJvnhgyQIAAAAAAAAAABTsuqVg16e5J4k97TWPlRVn5nk/fs5QVXdnOT7hsO3J3lGa+3R4fieqvrNJHcneWqSl1bVq1tr7zuM5AEAAFij2RGvH7x1lN+fgYk76ucvAAAAsFxr87bpxpDjyKz0k8/W2r9qrf231tpBthf77lwoYnrRQmHQ+TlOJ3nRcLid5MUHmAsAAAAAAAAAACZjo78WWVWV5OuGw3e11t56qXHD9X8yHD57iAMAAAAAAAAAgKvaRhcHJfmsJE8c+ncvGXv+9puSfOaqEgIAAAAAAAAAgLHYXj5krZ6y0H/XkrGLtz8lyfsPPx0AAAAAAAAAgBFqQ9t0Y8hxZDa9OOhJC/37l4y97zJxS1XVTUuG3Lif8wEAAAAAAAAAwCbY9OKgxy30H14y9pGF/vX7nOe+5UMAAAAAAAAAAGBcNr046ORCf2fJ2LML/VMryAUAAAAAAAAAYJSqtVTb/D27xpDj2Gx6cdCZhf7xJWNPLPQf3ec8y7YhuzHJPfs8JwAAAAAAAAAArNWmFwc9tNBftlXYdQv9ZVuQfYLW2v1Xur2q9nM6AAAAAAAAAADYCFvrTmCJxaKdm5aMXVz9574V5AIAAAAAAAAAAKOy6SsHvXOhf8uSsYu337uCXAAAAAAAAAAAxqm1edt0Y8hxZDZ95aD3J3lg6N+2ZOwzhsu/SPJnq0oIAAAAAAAAAADGYqOLg1prLcldw+EtVfUllxo3XH9+5aC7hjgAAAAAAAAAALiqbfq2YknyyiTfnnmuP11Vz2itPXr+xqo6leSnh8Nzw3gAAADYn9ls3RnszdZGf8+HVRvL4xQAAADYPC3JGD5asBzMoVtpcVBVfUWSmxeuesJC/+aqev7i+NbanRefo7X27qr6ySQvS/LUJG+uqh9L8r4kn53kpUm+cBj+E6219xzaHQAAAAAAAAAAgBFb9cpBdyT5lsvc9uVDW3TnZcb+QJJPS/JtmRcC/colxrwqyQ/uP0UAAAAAAAAAAJimUaxF3lqbtdZekOT2JHcleSDJznB5V5JntdbuaK2NYQEsAAAAAAAAAAA4EitdOai19vwkzz/E870hyRsO63wAAAAAAAAAAFeDai3V2rrTWGoMOY7NKFYOAgAAAAAAAAAA9k9xEAAAAAAAAAAATJTiIAAAAAAAAAAAmKjtdScAAAAAAAAAAMCKtSStrTuL5UaQ4thYOQgAAAAAAAAAACZKcRAAAAAAAAAAAEyUbcUAAAAAAAAAAKautZFsKzaCHEfGykEAAAAAAAAAADBRioMAAAAAAAAAAGCiFAcBAAAAAAAAAMBEba87AQAAAGAfZrP+2C3fEdoYB/k9AgAAAPSYDW3TjSHHkfGpIAAAAAAAAAAATJTiIAAAAAAAAAAAmCjbigEAAAAAAAAATFy1lmpt3WksNYYcx8bKQQAAAAAAAAAAMFGKgwAAAAAAAAAAYKIUBwEAAAAAAAAAwERtrzsBAAAAAAAAAABWrLV523RjyHFkrBwEAAAAAAAAAAATpTgIAAAAAAAAAAAmyrZiAAAAAAAAAACTN5JtxTKGHMfFykEAAAAAAAAAADBRioMAAAAAAAAAAGCiFAcBAAAAAAAAAMBEba87AQAAAAAAAAAAVqy1edt0Y8hxZKwcBAAAAAAAAAAAE6U4CAAAAAAAAAAAJsq2YgAAAAAAAAAAUzcb2qYbQ44jY+UgAAAAAAAAAACYKMVBAAAAAAAAAAAwUYqDAAAAAAAAAABgorbXnQAAAAAAAAAAAKtVraVaW3caS40hx7GxchAAAAAAAAAAAEyU4iAAAAAAAAAAAJgo24oBAAAAAAAAAExda/O26caQ48hYOQgAAAAAAAAAACZKcRAAAAAAAAAAAEyU4iAAAAAAAAAAAJio7XUnAAAAAOzDlu/5TELv73E2O9w8AAAAgKvHrM3bphtDjiPjE0UAAAAAAAAAAJgoxUEAAAAAAAAAADBRthUDAAAAAAAAAJi6lqSNYMuuEaQ4NlYOAgAAAAAAAACAiVIcBAAAAAAAAAAAE6U4CAAAAAAAAAAAJmp73QkAAAAAAAAAALBqLWlt3UnswRhyHBcrBwEAAAAAAAAAwEQpDgIAAAAAAAAAgIlSHAQAAAAAAAAAABO1ve4EAAAAAAAAAABYsdbmbdONIceRsXIQAAAAAAAAAABMlOIgAAAAAAAAAACYKNuKAQAAwDps+b4OHXofN7PZ4eYBAAAAjM+szdumG0OOI+OTSAAAAAAAAAAAmCjFQQAAAAAAAAAAMFGKgwAAAAAAAAAAYKK2150AAAAAAAAAAAAr1mbztunGkOPIWDkIAAAAAAAAAAAmSnEQAAAAAAAAAABMlG3FAAAAAAAAAACmrrV523RjyHFkrBwEAAAAAAAAAAATpTgIAAAAAAAAAAAmSnEQAAAAAAAAAABM1Pa6EwAAAAAAAAAAYMVmbd423RhyHBkrBwEAAAAAAAAAwEQpDgIAAAAAAAAAgImyrRgAAAAAAAAAwNS1Nm+bbgw5joyVgwAAAAAAAAAAYKIUBwEAAAAAAAAAwEQpDgIAAAAAAAAAgInaXncCAAAAAAAAAAAcgdbWnQFrYOUgAAAAAAAAAACYKMVBAAAAAAAAAAAwUbYVAwAAAAAAAACYutbGsa3YGHIcGSsHAQAAAAAAAADARCkOAgAAAAAAAACAiVIcBAAAAAAAAAAAE7W97gQAAAAAAAAAAFix2Syp2bqzWG42ghxHxspBAAAAAAAAAAAwUYqDAAAAAAAAAABgomwrBgAAAAAAAAAwda3N26YbQ44jY+UgAAAAAAAAAACYKMVBAAAAAAAAAAAwUYqDAAAAAAAAAABgorbXnQAAAAAAAAAAACvW2rxtujHkODJWDgIAAAAAAAAAgIlSHAQAAAAAAAAAABOlOAgAAAAAAAAAACZqe90JAAAAAAAAAACwYrOWVFt3FsvNRpDjyFg5CAAAAAAAAAAAJkpxEAAAAAAAAAAATJRtxQAAAAAAAAAAJq61WVqbrTuNpcaQ49hYOQgAAAAAAAAAACZKcRAAAAAAAAAAAEyU4iAAAAAAAAAAAJio7XUnAAAAAAAAAADAirUks7buLJYbQYpjY+UgAAAAAAAAAACYKMVBAAAAAAAAAAAwUbYVAwAAAAAAAACYutbmbdONIceRsXIQAAAAAAAAAABMlOIgAAAAAAAAAACYKMVBAAAAAAAAAAAwUdvrTgAAAAAAAAAAgBWbzZLM1p3FcrMR5DgyVg4CAAAAAAAAAICJUhwEAAAAAAAAAAATZVsxAAAAAAAAAICpa23eNt0YchwZKwcBAAAAAAAAAMBEKQ4CAAAAAAAAAICJUhwEAAAAAAAAAAATtb3uBAAAAGDUtnzvhhHofZzOZoebBwAAALA2bTZLy+b/rd98HnHofIIJAAAAAAAAAAATpTgIAAAAAAAAAAAmyrZiAAAAAAAAAABT19q8bbox5DgyVg4CAAAAAAAAAICJUhwEAAAAAAAAAAATpTgIAAAAAAAAAAAmanvdCQAAAAAAAAAAsGKtJbO27iyWayPIcWSsHAQAAAAAAAAAABOlOAgAAAAAAAAAACbKtmIAAAAAAAAAAFPXWpLZurNYzrZih250KwdV1WdU1U9W1b1V9UhVfaSq3lZV31tV1647PwAAAAAAAAAA2BSjWjmoqm5P8stJHr9w9bVJbh3aHVX1rNban64jPwAAAAAAAAAA2CSjWTmoqr4gya9mXhj0cJIfSPJlSb4qyS8Mw56c5PVVdf1akgQAAAAAAAAAgA0yppWDXpn5KkHnkjyztfaWhdveWFXvSfLjSW5J8pIkLz/6FAEAAAAAAAAANk+btbRq605jqdY2P8exGcXKQVV1a5KvHA5fdVFh0HmvSHLv0H9xVV1zFLkBAAAAAAAAAMCmGkVxUJJnL/RffakBrbVZktcMhzfkQjERAAAAAAAAAABclcayrdjTh8tHkvzRFcbdvdD/iiS/u7KMAAAAAAAAAADGos2SzNadxXJtBDmOzFiKg54yXL63tXbuCuPedYmYparqpiVDbtzruQAAAAAAAAAA2HxV9RlJ/kWS25N8RpKzSd6b5FeT/IfW2ulDmue5Sb41yednvhvWB5P8fpKfaa299TDmuJKNLw6qqpNJnjAc3n+lsa21j1bVI0muS/KkfUxzX2d6AAAAAAAAAACMTFXdnuSXkzx+4eprk9w6tDuq6lmttT89wBwnk/xakq+96Ka/NbRvqqofaq39cO8ce7G1ypMfksct9B/ew/hHhsvrV5ALAAAAAAAAAAAjVlVfkPnqQI/PvBblB5J8WZKvSvILw7AnJ3l9VR2k/uRVuVAY9KYkz07ytCQvSPK+zOt2Xl5VdxxgjqU2fuWgJCcX+jt7GH92uDy1jzmWrTJ0Y5J79nE+AAAAAAAAAICN0WYtrdq601iqtSPJ8ZWZrxJ0LskzW2tvWbjtjVX1niQ/nuSWJC9J8vL9TlBVtyX5puHwt5I8p7W2OxzfU1W/meSPMt/O7Mer6r+21j7WdW+WGMPKQWcW+sf3MP7EcPnoXidord1/pZb5Xm8AAAAAAAAAAIxYVd2a5CuHw1ddVBh03iuS3Dv0X1xV13RM9X3D5W6Sf75QGJQkaa09mOSlw+ENma8mtBJjKA56aKG/l6Warhsu97IFGQAAAAAAAAAAV49nL/RffakBrbVZktcMhzfkQjHRngxbkX3VcPi7w8I0l/LaJB8f+l+/nzn2Y+OLg1prZ5I8OBzedKWxVXVDLhQH3bfKvAAAAAAAAAAAGJ2nD5ePZL6t1+XcvdD/in3O8bRc2Pnq7ssNaq3tJHnr+ZjOFYqW2vjioMH5pZpurqrtK4y75RIxAAAAAAAAAABXtzYbT1utpwyX722tnbvCuHddIma/c1x8nivNs53kc/Y5z55cqdBmk/xB5pVb1yX54iR/eJlxty3033yI8x873/nABz5wiKcFALjgzK5dUQFGaVbrzgBWp7V1ZwBwVbv//svtPAAA0O+i//M+drlxTM/ZnElG8Kf+2ZxZPLyx6sqfv11hy66/pqpOJnnCcHjFuNbaR6vqkcxrVZ601zkGi+OX5be4M9aTkrxzn3MtNZbioN9I8i+H/rfmEsVBVbWV5JuHw48ledMhzv+p5ztPe9rTDvG0AAAAAABwaU960p3rTgEAmL5PTfLn606Co3FP3rjuFHrcs4cx+/n23uMW+nv51vb54qDr9zHHfud5ZKG/33n2ZBTbirXW3pbk94fDF1TVl15i2PfkwrJMP9Vae+xIkgMAAAAAAAAAYAxOLvR39jD+7HB5aoXznF3o73eePRnLykFJ8l2ZbxV2KsnvVNWPZr460Kkkz03ywmHcu5O84pDn/uMktw79/5dk95DPD4zTjblQqXprkg+uMRdgurzWAEfBaw1wFLzWAEfBaw1wFLzWAEdhla81x3Jh95w/PsTzspk+mP1vibUpbszh12gs7ll2fA/jTwyXj65wnhML/f3OsyejKQ5qrb2jqr4xyS8l+aQkP3qJYe9Ocntr7aFDnvtskrcf5jmB8btob8sP7mcvS4C98loDHAWvNcBR8FoDHAWvNcBR8FoDHIUjeK2xldhVorV2LslY36tWkfdiPcletvC6brjcyxZkvfNct9Df7zx7Moptxc5rrf1Wks9P8u8yLwQ6neRjmRfuvDTJF7bW3ru+DAEAAAAAAAAA2ESttTNJHhwOb7rS2Kq6IRcKd+7b51SLhU1XnCefuLLTfufZk9GsHHRea+3Pk7xkaAAAAAAAAAAAsFf3Jnl6kpurantYXelSbrkoZj/eeZnzXGmec0lWsiDOqFYOAgAAAID/397dR/tW1fUef3/weHjSlCRsoDJEMTllFIJ5UB691C3BMjQLG+pBzLQ7uEJIKNZF7MHE0stVU0QLzHtVKIiALPEBFMIAh3XvSOBcRW+AiTyIyIMc0e/9Y83tXm1/j+e3z+Hs336/xjhjz/lb3zV/a/9xvmOtub9rTkmSJEmSpBlc0X7uDOw3Iu6QXvvKKb/jGmDTgHH+gyRrgfUL51TVpmGxs7A4SJIkSZIkSZIkSZIkSavF3/TaxwwKSLId8NLWvQv41DRfUFXfAj7RuocnGba12FHAD7X2BdN8xzQsDpIkSZIkSZIkSZIkSdKqUFVXA59p3WOTHDAg7ERgXWufUVXf6R9MsiFJtX9vHPJVf9J+rgHeleRhS8bYFXhL694FvG+632RyFgdJkiRJkiRJkiRJkiRpNXkNcD9d4c7Hkrw+yfokhyU5Ezi9xW0E/nRzvqCqPgl8uHV/Ebg0yS8m2T/JMcBngT3a8ddV1Tc295cZZ82WGliSJEmSJEmSJEmSJEna1lTV55P8KvBBum29/mhA2EbgiLZF2OZ6eRv/ucBh7V/f94Dfr6ozZ/iOsSwOkqTNVFU3A3mor0PSfDPXSNoazDWStgZzjaStwVwjaWsw10jaGsw10pZXVRcl2YduFaEjgMcDm4AvAucB76yq+2b8jvuBI5K8GNgA/BTwaOBWuq3N3llVV83yHZNIVW3p75AkSZIkSZIkSZIkSZL0ENjuob4ASZIkSZIkSZIkSZIkSVuGxUGSJEmSJEmSJEmSJEnSnLI4SJIkSZIkSZIkSZIkSZpTFgdJkiRJkiRJkiRJkiRJc8riIEmSJEmSJEmSJEmSJGlOWRwkSZIkSZIkSZIkSZIkzSmLgyRJkiRJkiRJkiRJkqQ5ZXGQJEmSJEmSJEmSJEmSNKcsDpKkKSXZI8mrk3wkyQ1J7k3y7SQ3J7kwydFJ1kwx3k8keU+SLya5P8ltST6d5DenGUfS/Gn55k+SXNdyzZ1Jrk7y2iQ7PdTXJ2nblOTpSU5J8tEkNyV5IMk9STYmOTvJQVOO9/NJzm/3Og+0n+cn+fkt9TtIWtmSnJ6kev8OneAcc42ksZLsmuR3klyZ5GstX3w1yT8leWuSAyYYw3wjaagka5Mcm+Tvk/x773nqhiR/nmT9hOOYa6RVJsluSY5M8qY2J3N775no7M0Yb+Y8kmSnJCe1OeU7Wz67rs057zHtNUla2VJVD/U1SNKKkeRNwO8CGRN6LfCCqvq3MeMdC7wL2H5IyGeBI6vqjmmvVdLKluQI4H8CjxoScgPw3Kq6cetdlaRtXZLLgYMnCP1L4BVVtWnEWAHeA7xyxDjvBV5VPlhKapL8FN3zUP9Fh8Oq6rIh8eYaSRNJ8ivAu4HHjAi7sKqeP+R8842kkZI8AbgE+MkxoW8HThyUK8w10uqVZNT/6XOqasOE4yxLHknyZLqc9tQhId8EXlxVfzfJdUla+Vw5SJKmsztdYdC9wAeBY4ADgf2BlwDXtLj9gY8necSwgZL8Z7obuO2BW4H/CjwT+AXg/Ba2Hjg/iflaWkXaH9XOpSsMugd4A/As4D8BZ7WwpwKXjMozklalx7WfXwXOAF4I/AxwAPDbwC3t+EuAs8eM9QcsTkR9Hji6jXV069OO//4yXLekOdCeW86iKwz6+oSnmWskjZXkpcCH6QqDvg6cBvwssB9wBN2cyqXAd0YMY76RNFRbwb1fGPS/gQ10z1I/B7yJbk4Y4ATgtUOGMtdIArgJ+NhmnjtzHmlzxhezWBh0Ft3c8rPo5prvoZt7Pi/JPpt5nZJWGFcOkqQpJHkLcAfw7qr61oDjDwP+F/Ci9tF/q6ofuEFrD5vXAXsBdwNPr6ovLYl5F/BbrfuyqvrAsv0ikrZpST4FHAo8CBxcVVctOX4ScHrrnlpVb9q6VyhpW5XkYuADwF9X1XcHHN8VuBL4sfbRwVX1mQFxe9Hdq6yhWwHk4Kq6v3d8J+ByuoLoB4G9l97LSFp9khxP9yb99cAFwOvboYErB5lrJE0iyTq6P4RtD3wGeF5VfXNI7NpBKyOabySNk+QFwF+17lXAQUufqZLs1449HPgGsFtVPdg7bq6RVrEkp9G9QH5NVd2a5InAl9vhiVYOWq48kuSNwKmt+ztV9dYlxw8APt2+51NV9ZyJf1FJK5YrUUjSFKrq5Ko6fVBhUDv+XbqCnoWJqBcOGeqX6QqDAN485CHwJLqHzIW2pFUgyTPoCoMA3r+0MKj5U7qHRIDjkzx8a1ybpG1fVR1ZVecOKgxqx28HTux9NOxe5QQWtwQ6rj8R1ca5DziuddcAx2/+VUuaB20bjoUXI17N4jPRKOYaSZN4B11h0O3AUcMKgwBGbJlqvpE0zrN77TcPeqaqqs/RrcQBsAuw95IQc420ilXVqVV1cVXdOsMwM+eRNlf8mta9jm4ueem1XgW8v3UPa8WPkuacxUGStMyq6g66ZWcBnjwk7Pm99tlDxrmPblshgKclecqyXKCkbV0/P/zFoICq+h7dyiDQTUYduoWvSdJ8uazX/oF7lba3/S+17vVV9dlBg7TPb2jd57fzJK1efwY8gu6N2MvGBZtrJE0iyd50W2AAvLMVOk87hvlG0iTW9to3jojrv+S5/ULDXCNpVsuYRw4FHt3a57S55EHO7rWPmvqCJa04FgdJ0pax8GA47KbroPbzhqr62ohxLu+1D5z5qiStBAv54V7gcyPizA+SNld/0nvQvcqewONa+/IBx/sWjj8eeOJslyVppUryIuBI4E4mX/XUXCNpEr/Sa5+30EiyS5KnJHnMBGOYbyRNYmOv/aQRcQsvWBTwf3ufm2skzWq58shBA+IGuZZuDhqcX5ZWBYuDJGmZJdkNWNe61w84/gi6G7aBx5foH183NErSPFn4v/7F/r71A5gfJG2uQ3rtQfci68YcZ8hxc5G0CiV5NHBG655cVbdNeKq5RtIk1ref3wSuS/LrSf6FrhhxI3B7khuTnNrmWwYx30iaxIeAu1v75CQPWxqQZF/giNb9cFXd3TtsrpE0q+XKIxON0+aeF1ZDMxdJq4DFQZK0/E5icU/YcwccfzywsMzjzWPGuqnXfsKM1yVpG5dkB2DX1h2ZH6rqGyy+2WF+kDSRJNsBr+t9NOhepZ9TvFeRNM7pwI8C/wi8f4rzzDWSJvHj7edXgHcAHwT2WRKzJ/BG4Kokuw8Yw3wjaaxW4LwBuB94NnBNkpcmWZ/k8CSn0q3AsRb4Z+C3lwxhrpE0q+XKIwv9e6vqrgnH+ZEk24+MlLTiWRwkScsoyTOB41v3ZuDPBoQ9ste+Z8yQ9/baw96AkzQ/pskPsJgjzA+SJnUC8DOtfUFVXTsgxnsVSRNJciDwCuBB4FVVVVOcbq6RNIkfbj/3Bv4LcBfwKmA3YAfgGcBHW8zTgPNaMXSf+UbSRKrqAmB/uoLnnwbOAa4CLqUrQryPrijowKr62pLTzTWSZrVceWRhnGnmlweNI2nOWBwkScskyWOBv6JbNaiAl1XVfQNCd+i1N40Z9oFee8fZrlDSCjBNfoDFHGF+kDRWkkOAP27drwOvHhLqvYqksZKsBd5Ltyrq26vq/0w5hLlG0iR2bj+3B74L/EJVnVlVt1XVA63Q+UgWC4SeBRy1ZAzzjaSJJHk48GLgeSyu/N73WOBo4NABx8w1kma1XHlkYZxp5pcHjSNpzlgcJGkuJVmTpJbh34YJv++RwCV0W4YBnFJVnxwS/u1ee+2YofvLON4/ybVIWtGmyQ+wmCPMD5JGSvITwAV0RcwPAC+qqluHhHuvImkSpwDrgH8DTtuM8801kibRzxXnVdVnlwZU1ffotnhfcPSIMcw3kgZKsjPwceANwGPotk5dR5cXHgX8HHAF3YplFyV5zZIhzDWSZrVceWRhnGnmlweNI2nOWBwkSTNKsgNwIbBf++htVfXHI075Vq89bpnGnXvtSZaAlLSyTZMfYDFHmB8kDRPQNvMAAAgBSURBVJVkT+BjwC50b9wfXVWXjzjFexVJIyXZG3h96x5XVfeOih/CXCNpEv1c8dFhQVX1r8AtrfuMEWOYbyQNcxpwcGsfW1UnV9X1VbWpqu6uqkuBw4BP0a0q9LYk+/TON9dImtVy5ZGFcaaZXx40jqQ5s+ahvgBJ2hKq6sEk65ZhqH8fdTDJGuBcugdDgPdV1Yljxry513780KjOE3rtm8bESlrhqurbSW4HdmVMfkiyC4sPb+YHSQMl2Z3u7dfd6bY9fXlVXTDmNO9VJI1zAt1bqDcCOyX5tQExT+u1n5PkR1v7olZMZK6RNImbgIX8cfOowBb7OGC3JZ+bbySNlCTAMa27sarOGRTX5px/j24Foe3aOSe0w+YaSbNarjxyM/BMYOckj66quyYY57aqemBEnKQ5YHGQpLlVVddvyfGTbAf8Jd0e1AAfAX5zguu6J8lNdDdde48J7x+/bnOuU9KKcx1wELBXkjVV9eCQOPODpJGS7ApcCjypfXRcVX1gglO/0Gt7ryJpkIWl558EfGiC+N/rtfcE7sVcI2ky/8riSkAPGxO7cHzpM5T5RtI4jwV+uLU/Pyb2c712P2eYayTNarnyyBeAF/TifmBbVvj+y+9PHjKGpDnktmKStPnOBBbekL0YeEnb534SV7SfT+29QTvIIb32lVNen6SVaSE/7MzidoWDmB8kDZXkUcA/AD/ePnpdVb1rwtO/DHy1tQ8ZFcjisvu3AF+Z5holrXrmGkmT+HSv/eShUZ2FguhblnxuvpE0Tr+ocNxL9Q8fcp65RtKsliuPXNFrjxpnfxZXpnd+WVoFLA6SpM2Q5G3AK1r3E8ALq+o7UwzxN732hiHfsRPwotb9QlVtnPY6Ja1I/fxwzKCAtnLZS1v3Lrr97iUJ+P49xCXA09tHf1hVb5n0/Koq4MLW3TvJ+iHfs57FN9UubOdJWgWqakNVZdQ/4LTeKYf1jn2ljWGukTSJvwUW5luOGhaU5BDgMa37mf4x842kCdwJ3N3aB7TVNIbp/6H9ywsNc42kWS1jHrkM+GZrv6xtnTjIhl573Bb0kuaAxUGSNKUkb2RxL+l/BH5pM/ZivQD4Umu/Psmgt9/eCuzSa0taBarqahYns49NcsCAsBOBda19xpTFiZLmWJK1dPcZz24fnVFVv7sZQ/13Ft+CfUeSHZd8z47AO1r3wRYvSdMy10gaqaruAN7Xuj+b5NeWxiR5JP8xP5w5YCjzjaSh2mrwl7Tu7sAbBsUl2QXov3hx8ZIQc42kWc2cR6pqE/A/Wncd8NqlMW3O+djWvbyqrpn90iVt62JRsiRNLslxLN5U3QL8KosV2MPcMOgP90meC1xEV6h5K/AHwNV0BUG/weKesFcAh1bVd2f+BSStCEn2pVvKdUfgHuCP6FYH2pFuO8NXttCNwP5V9a2H4jolbXuS/DWLb9V/EjgeGPXQt2nY6oRJ3gy8rnU/TzcJ/iW6LT1OBvZtx95cVafMeOmS5kx7qeLU1j2sqi4bEmeukTRSkh8BrgX2oPsj2HuA8+lW+fhJulyx8Pb8u6vqt4aMY76RNFSSvYHPATu1jy4CzgFuBHYA1tM9X+3Rjn+iqg4fMI65RlqlkhwI7NX7aFcWX/y+ksWCZwCq6uwh48ycR1rx9LXAj7WP3gt8GLgfOAw4BXhE6z+rqv55kt9R0spmcZAkTSHJZYzf63WpPReWzh8w3m8A7wTWDjn3auCIqrp9yu+UtMIleR7wQeCHhoRspMsPX9x6VyVpW5dk2ge8/1dVTxwy1nbAWcDLR5z/fuCV7U1bSfq+KYqDzDWSxkqyjm6Lsb1GhP058KphK6uabySNk+Rw4EN0f9Af5ZPAC6vqGwPGMNdIq1SSs4GXTRrftmMeNM6y5JEkewF/BzxlSMjdwK9X1dJV0CTNKbcVk6SHUFWdBexHd6N3I/Bt4A661YJeDTzbwiBpdaqqi4B9gLfTFQLdB9xF98bHycC+FgZJ2pKq6ntVdSxwBN2e918FNrWfFwLPrapXOKEtaRbmGkmTqKrrgJ8GTgL+CbiTLlfcDHwEeE5VHTtqy2XzjaRxqurjdCuRnQxcBtwGfIduZY0vA+cCzwcOH1QY1MYw10iayXLlkTZ3vC9dTruWbm75PuAGujnnfSwMklYXVw6SJEmSJEmSJEmSJEmS5pQrB0mSJEmSJEmSJEmSJElzyuIgSZIkSZIkSZIkSZIkaU5ZHCRJkiRJkiRJkiRJkiTNKYuDJEmSJEmSJEmSJEmSpDllcZAkSZIkSZIkSZIkSZI0pywOkiRJkiRJkiRJkiRJkuaUxUGSJEmSJEmSJEmSJEnSnLI4SJIkSZIkSZIkSZIkSZpTFgdJkiRJkiRJkiRJkiRJc8riIEmSJEmSJEmSJEmSJGlOWRwkSZIkSZIkSZIkSZIkzSmLgyRJkiRJkiRJkiRJkqQ5ZXGQJEmSJEmSJEmSJEmSNKcsDpIkSZIkSZIkSZIkSZLmlMVBkiRJkiRJkiRJkiRJ0pyyOEiSJEmSJEmSJEmSJEmaUxYHSZIkSZIkSZIkSZIkSXPK4iBJkiRJkiRJkiRJkiRpTlkcJEmSJEmSJEmSJEmSJM0pi4MkSZIkSZIkSZIkSZKkOWVxkCRJkiRJkiRJkiRJkjSnLA6SJEmSJEmSJEmSJEmS5pTFQZIkSZIkSZIkSZIkSdKcsjhIkiRJkiRJkiRJkiRJmlMWB0mSJEmSJEmSJEmSJElz6v8DAKI153JiN4gAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "init()\n", "time_loop(1000)\n", "plot_ρs()" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [], "source": [ "assert np.isfinite(dh.gather_array(ρ.name)).all()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.9.7" } }, "nbformat": 4, "nbformat_minor": 2 }