diff --git a/about_initial_conditions.ipynb b/about_initial_conditions.ipynb deleted file mode 100644 index dea3767..0000000 --- a/about_initial_conditions.ipynb +++ /dev/null @@ -1,3865 +0,0 @@ -{ - "cells": [ - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "# About initial conditions\n", - "\n", - "When creating spatio-temporal trajectories of the semi-linear PDEs available in\n", - "`Exponax`, we need to decide on two things:\n", - "\n", - "1. The concrete dynamics (this include the PDE, the boundary conditions, the\n", - " constitutive and numerical arguments)\n", - "2. The initial condition\n", - "\n", - "The second point affects how the state evolves in time. Some dynamics are highly\n", - "sensititive to it, some mostly ignore it. There are many options to define\n", - "initial conditions. This notebook introduces you to the ones available in\n", - "`Exponax`." - ] - }, - { - "cell_type": "code", - "execution_count": 1, - "metadata": {}, - "outputs": [], - "source": [ - "import jax\n", - "import jax.numpy as jnp\n", - "import matplotlib.pyplot as plt\n", - "from IPython.display import HTML" - ] - }, - { - "cell_type": "code", - "execution_count": 2, - "metadata": {}, - "outputs": [], - "source": [ - "import exponax as ex" - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "Let's again start with the advection equation in 1d" - ] - }, - { - "cell_type": "code", - "execution_count": 3, - "metadata": {}, - "outputs": [], - "source": [ - "DOMAIN_EXTENT = 1.0\n", - "NUM_POINTS = 20\n", - "DT = 0.01" - ] - }, - { - "cell_type": "code", - "execution_count": 4, - "metadata": {}, - "outputs": [ - { - "name": "stderr", - "output_type": "stream", - "text": [ - "2024-03-21 17:17:14.932545: W external/xla/xla/service/gpu/nvptx_compiler.cc:679] The NVIDIA driver's CUDA version is 12.2 which is older than the ptxas CUDA version (12.3.52). Because the driver is older than the ptxas version, XLA is disabling parallel compilation, which may slow down compilation. You should update your NVIDIA driver or use the NVIDIA-provided CUDA forward compatibility packages.\n" - ] - } - ], - "source": [ - "adv_stepper_1d = ex.stepper.Advection(1, DOMAIN_EXTENT, NUM_POINTS, DT)" - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "Let us define two initial condtions for it\n", - "\n", - "$$ u_0(x) = \\sin(\\frac{2\\pi}{L} x) $$\n", - "\n", - "and\n", - "\n", - "$$ u_0(x) = 2 \\cos(\\frac{2\\pi}{L} x) $$\n", - "\n", - "where $L$ is the length of the domain." - ] - }, - { - "cell_type": "code", - "execution_count": 5, - "metadata": {}, - "outputs": [], - "source": [ - "grid = ex.make_grid(1, DOMAIN_EXTENT, NUM_POINTS)\n", - "u_0_first = jnp.sin(2 * jnp.pi / DOMAIN_EXTENT * grid[0:1])\n", - "u_0_second = 2 * jnp.cos(2 * jnp.pi / DOMAIN_EXTENT * grid[0:1])" - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "and produce trajectories for both of them." - ] - }, - { - "cell_type": "code", - "execution_count": 6, - "metadata": {}, - "outputs": [], - "source": [ - "trj_first = ex.rollout(adv_stepper_1d, 100, include_init=True)(u_0_first)\n", - "trj_second = ex.rollout(adv_stepper_1d, 100, include_init=True)(u_0_second)" - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "visualizing them in a spatio-temporal" - ] - }, - { - "cell_type": "code", - "execution_count": 7, - "metadata": {}, - "outputs": [ - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAA1UAAAE8CAYAAAAom5t9AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/SrBM8AAAACXBIWXMAAA9hAAAPYQGoP6dpAABcOklEQVR4nO3de3RU5b0//s/cZ89MMjPJ5EIgJEhRg0FQrBy0Hlu/VI568NibVK0g9fK1Xo6a1aNQlYhWUL+Waj1YlrRUv6f2aG3Vr6dSvKAcW4XSAuGYEi1GCEgIySRzSea2Z2bv3x/83M/+7MxQyJBMEt6vtVjr8+z37Ev2PLOf/WSGiUlVVZUAAAAAAABgSMzFPgAAAAAAAICxDJMqAAAAAACAAmBSBQAAAAAAUABMqgAAAAAAAAqASRUAAAAAAEABMKkCAAAAAAAoACZVAAAAAAAABcCkCgAAAAAAoACYVAEAAAAAABQAkyqAUWbfvn1kMpno2WefLfahAAAAjDomk4keeOCBYh8GAINJFUARPPvss2QymXL+W7p06Qnf38qVK+nVV1894dsFAIDR6cMPP6RvfvObVFdXR06nkyZOnEhf/epX6amnnir2oQGMS9ZiHwDAyezBBx+kKVOmsGVnnHEGPfvss2Sz2U7YflauXEnf/OY36Yorrjhh2wQAgNHpgw8+oK985Ss0efJkuvHGG6m6upoOHDhAW7dupSeffJJuv/32Yh8iwLiDSRVAEV1yySV0zjnnDGndWCxGbrf7BB8RAACMdQ8//DB5vV7685//TD6fj2Xd3d3FOSiAcQ4f/wMYZXL9n6rrrruOPB4Ptbe306WXXkolJSV0zTXXEBHRnj176Bvf+AZVV1eT0+mkSZMm0be//W2KRCJEdOSz57FYjJ577jntI4bXXXddEX4yAAAYCe3t7XTGGWcMmlAREVVWVrL2L3/5S5o9ezZJkkRlZWX07W9/mw4cODBovT/96U906aWXkt/vJ7fbTWeeeSY9+eST7DHvvPMOXXDBBeR2u8nn89G//Mu/UFtbG3vMAw88QCaTiT755BO67rrryOfzkdfrpSVLllA8HmePTaVSdNddd1FFRQWVlJTQ5ZdfTp999tkQzwrA8MI7VQBFFIlEKBgMHtNjM5kMzZ8/n770pS/R448/Ti6Xi2RZpvnz51MqlaLbb7+dqqur6eDBg/S73/2OwuEweb1e+o//+A+64YYb6Nxzz6WbbrqJiIimTp06nD8WAAAUUV1dHW3ZsoVaW1upsbEx7+Mefvhhuv/+++nKK6+kG264gXp6euipp56if/zHf6SdO3dqk7K33nqL/vmf/5kmTJhAd9xxB1VXV1NbWxv97ne/ozvuuIOIiN5++2265JJL6JRTTqEHHniAEokEPfXUU3T++efTjh07qL6+nu37yiuvpClTptCqVatox44d9LOf/YwqKyvp0Ucf1R5zww030C9/+Uu6+uqr6bzzzqN33nmHLrvsshN+vgBOCBUARtwvfvELlYhy/tu7d69KROovfvEL7fGLFy9WiUhdunQp287OnTtVIlJfeumlo+7P7XarixcvHoafBAAARps333xTtVgsqsViUefOnavefffd6htvvKHKsqw9Zt++farFYlEffvhhtu6HH36oWq1WbXkmk1GnTJmi1tXVqaFQiD1WURStnjVrllpZWan29vZqy3bt2qWazWZ10aJF2rLm5maViNTvfve7bFtf+9rX1PLycq3d0tKiEpF6yy23sMddffXVKhGpzc3Nx3dSAIYZPv4HUERr1qyht956i/07mu9973us7fV6iYjojTfeGPSxCQAAODl99atfpS1bttDll19Ou3btoscee4zmz59PEydOpNdee42IiF5++WVSFIWuvPJKCgaD2r/q6mqaNm0avfvuu0REtHPnTtq7dy/deeedgz5OaDKZiIjo0KFD1NLSQtdddx2VlZVp+Zlnnklf/epXacOGDYOO8eabb2btCy64gHp7eykajRIRaev867/+K3vcnXfeOfQTAzCM8PE/gCI699xzB31Rxb59+3I+1mq10qRJk9iyKVOmUFNTE61evZqef/55uuCCC+jyyy+n73znO9qECwAATj5f/OIX6eWXXyZZlmnXrl30yiuv0I9//GP65je/SS0tLbRnzx5SVZWmTZuWc/3Pv4G2vb2diOioHyPs6OggIqLTTjttUNbQ0EBvvPHGoC9Xmjx5Mnuc3+8nIqJQKESlpaXU0dFBZrN50MfVc+0DYDTApApgjHA4HGQ2D35z+Uc/+hFdd9119P/+3/+jN998k/71X/+VVq1aRVu3bh00CQMAgJOL3W6nL37xi/TFL36RTj31VFqyZAm99NJLpCgKmUwm+v3vf08Wi2XQeh6PZ1iPK9c+iYhUVR3W/QIMF0yqAMaBGTNm0IwZM+i+++6jDz74gM4//3xau3Yt/fCHPyQi8RENAAA4eX3+yYhDhw7R1KlTSVVVmjJlCp166ql51/n8naLW1laaN29ezsfU1dUREdHHH388KPvoo48oEAgc958AqaurI0VRqL29nb07lWsfAKMB/k8VwBgWjUYpk8mwZTNmzCCz2UypVEpb5na7KRwOj/DRAQBAMbz77rs53/H5/P8pnXbaafT1r3+dLBYLrVixYtBjVVWl3t5eIiI6++yzacqUKfTEE08MGkc+X2/ChAk0a9Yseu6559hjWltb6c0336RLL730uH+GSy65hIiIfvKTn7DlTzzxxHFvC2Ak4J0qgDHsnXfeodtuu42+9a1v0amnnkqZTIb+4z/+gywWC33jG9/QHjd79mx6++23afXq1VRTU0NTpkyhOXPmFPHIAQBguNx+++0Uj8fpa1/7Gp1++ukkyzJ98MEH9OKLL1J9fT0tWbKEfD4f/fCHP6Rly5bRvn376IorrqCSkhLau3cvvfLKK3TTTTfR97//fTKbzfTTn/6UFixYQLNmzaIlS5bQhAkT6KOPPqK//vWv9MYbbxAR0f/5P/+HLrnkEpo7dy5df/312leqe71eeuCBB477Z5g1axZdddVV9PTTT1MkEqHzzjuPNm3aRJ988skJPlsAJwYmVQBj2MyZM2n+/Pn0X//1X3Tw4EFyuVw0c+ZM+v3vf0//8A//oD1u9erVdNNNN9F9991HiUSCFi9ejEkVAMA49fjjj9NLL71EGzZsoGeeeYZkWabJkyfTLbfcQvfdd5/2LX5Lly6lU089lX784x/TihUriIiotraWLr74Yrr88su17c2fP5/effddWrFiBf3oRz8iRVFo6tSpdOONN2qPmTdvHm3cuJGam5tp+fLlZLPZ6MILL6RHH32UpkyZMqSfY/369VRRUUHPP/88vfrqq3TRRRfR66+/TrW1tUM/OQDDxKTifwQCAAAAAAAMGf5PFQAAAAAAQAEwqQIAAAAAACgAJlUAAAAAAAAFKOqk6r333qMFCxZQTU0NmUwmevXVV//uOps3b6azzz6bHA4HfeELX6Bnn3122I8TAAAAAAAgn6JOqmKxGM2cOZPWrFlzTI/fu3cvXXbZZfSVr3yFWlpa6M4776QbbrhB+zpPAAAAAACAkTZqvv3PZDLRK6+8QldccUXex9xzzz30+uuvU2trq7bs29/+NoXDYdq4ceMIHCUAAAAAAAA3pv5O1ZYtW2jevHls2fz58+nOO+/Mu04qlaJUKqW1FUWhvr4+Ki8vJ5PJNFyHCgAAOaiqSv39/VRTU0Nm88n733oxNgEAjA4nalwaU5Oqrq4uqqqqYsuqqqooGo1SIpEgSZIGrbNq1SrtD9oBAMDocODAAZo0aVKxD6NoMDYBAIwuhY5LY2pSNRTLli2jpqYmrR2JRGjy5Mn0zNvbyeX2UG9cZo/vi6e1uqc/ybJgv/itYjDK14snRDul2wYRUVLXTiUzLEvHxT7kRJRlmURcqxU5wddLxkSW4ceiZMT+lCzPVCWrqxU6VibdzN1ktrDMbLGL2mrjmVVkNqebZ3YxCbZKLpbZpVKxnsvJModTdFuni+/PoWu7JDvLAqWiHShxsKyiROyjzLDNcpdYr1zimU8Sx+Ky8d9ueHRttyGz6J5rUzLCMnOqX6uVaB/LlHBQq7PhHpbJ4bBWp/r6WZbsi+TNEmHRB1Nh3ueTEdHnU9EUyxIx0c+imSzLYlnxqeJkln/COJYR/U5WeKZvpw2Zfg/H3nP5fxy1GDKbWbwjYDfzdwf0bbeVP39Oiy6z8PVKrWIvktvQP0tFv3N6eR90+EQflHyGPl9WItYr8+bN7D4fyyy+Cq02+wIsM5eWabXiKGGZ6hT7UHSvRSKigbSSsyYiiuva4YS41vUmRF+JxwbopnmzqaSE7/Nkk29suoYmkn2UfzGvS/famF7Kr7OTasTzGji9nGXlDZO12nvaFJbZJp8mGlU8y/prtborzsfQAxExxrX3xVjWelBc93Yf4NfZ3kPiOhg6dIhlkY42rc6k+DZhbLM6xH2It66BZf4JE7S6fAK/Pk2vFdfExon8Gjy1TGyz1stfD1UucY9gDR3gB3N4r1am93/MosjHIutt28+y4Ee9Wv1ZJx/Pd+vuTePKqPjfPWOCTAo9TwcLHpfG1KSqurqaDh8+zJYdPnyYSktLc75LRUTkcDjI4XAMWu5ye8jlKaG4mU86nCbRdmT56bFlxA2SNc1vMK2qyDJZPqmy6CY2FoUPCNmMGDzNGb6eWXfzafyvb6aM2I6J+E2dvj3oQyQm3a2p6cRMqky6iZO+HpTZ+I2iWdc22/jzZ7aLSZbFztezOKy6mu/P6tQ9R4ZJlU0S/cDh4n3C6dbd0Lr5ei7dpMptmHB5dJMs48TJYxdtj3FSZRPPp8nGJyRmu8gUhU9ylLQ4T1mZnxc5JX6mpJP3T5vuPNns/Gew2ERfslp5n7dYRGYxPu9m0X/Sho6mmHT91dDPMrrHmkyGfq1rGzP93GzIkyrDcdp1H7OyGz5ypW87Tfz5k3RtybCeS3eeJMM5c1pE22k4106baEuG58ipe/6cTt4/nbp+bTf8AsLiFv3F7OG/uDCXeLRacXhYpkpiYDFOqky6iZNJ5s+EWZelreJ6lrDwaxsRnfQfccs3NtnJPOonVQ7dcyeZeB936/q4x8b7cYmu75Ya+qpN3z9LeH/Mloo+GLPwMdStiGudJPPzZneJa6vVydezOERmHH/YuJUZ3Hdh7NI/t8bn3eIQfdBq+CWw3SX6pOThN97uEvHYEsMvbEvduklVhvdrion9pQ2vB0X3WkkZXkcJ3WvM+Ppz6MamDGFSdbwKHZdG95XbYO7cubRp0ya27K233qK5c+cW6YgAAAAAAOBkV9RJ1cDAALW0tFBLSwsRHfnK9JaWFtq//8hbncuWLaNFixZpj7/55pvp008/pbvvvps++ugjevrpp+nXv/413XXXXcU4fAAAAAAAgOJOqv7yl7/QWWedRWeddRYRETU1NdFZZ51Fy5cvJyKiQ4cOaRMsIqIpU6bQ66+/Tm+99RbNnDmTfvSjH9HPfvYzmj9/flGOHwAAAAAAoKj/p+rLX/7yoP8rpPfss8/mXGfnzp3DeFQAAAAAAADHbkz9nyoAAAAAAIDRBpMqAAAAAACAAmBSBQAAAAAAUABMqgAAAAAAAAqASRUAAAAAAEABivrtf8W0/o/7yCq5KR6X2fJkTPz19KQhSyXEX2SX40mWyfGIVmcSMZZl5YRWpxMDLFMycs7a2FaVLMtURaFjYTKbDW3x17ctdv7XxM36vzSuq4mIbJL4S+DG9ayS+GvidpeXZXbdXwl3SLy7OV1iH043/4vhLl1WUcr/Qnmlrl1Zyv8KeblbrBdw8Z8h4BL78Et8f26bOE8eOz9nHl1mSURYZkoEtdoc7WeZEhGZEuphWSbUrdVyKMSyZG9UV0fyZolQgmch0SeTYd4/k5GUWE/Xx4mIImnRt2JZ/m2ciazoZwMZ3udkRc1ZG9uGTdKx9dzBv/Gx6P7QuWTmf/XcrmvbDZnHKrYkWfhW3bqNem38L9NLuj7p9PI+6PSJfuf08z4o+cXrw1leytcr9+bN7H6/Vlv8lSwz+ytE7Q2wTHGUaLUq+Xgmif2F0/zMD8iiHTNkobDoI8HOMMuCuutib4xfs7qjSV0t+lyPrjZeH2Hk6Ps7EVGj7lo6eVIJyyqmiz4XaKxjma9hqlbb6hv4TqpFli3j6x2KiTF0m+6aRES0p1f0i11bwyxr7dim1cHOKMv6DnZqdXhfK8sySfS18cLqdLO2r75Rq8sm1rAsUCOurY11PpbNrBXtaeV8m3W663y1m9+vWPs6RKOL97N0S5tWh9vaWbanVazXs5vfB+z/TNwztEb568E4FsPYgHeqAAAAAAAACoBJFQAAAAAAQAEwqQIAAAAAACgAJlUAAAAAAAAFwKQKAAAAAACgAJhUAQAAAAAAFACTKgAAAAAAgAJgUgUAAAAAAFAATKoAAAAAAAAKgEkVAAAAAABAAazFPoBi2bPrMzLbXSTHI2x5JhHT6qycYFk6MaDVSkZmmb5tzFQlq6uVYz5Gk9msqy0ss9glrTZb7SzTt22SJ+96VsnNMrvLq6udLHNIoqs4XXx/TrdNq12GrKLUodWVuvpIW+yj3M3XC+i2E3DZWOaXRNtt478X8NhF22PILAnxXJsSQZaZo/1arUR4poR6tDoT6maZHAppdbI3yrJkbyRvlgiJvpUMJfl6YdFORlJ8vVhaqyPpLMtiWVU8Lsv72UBGtGVFZZm+bcx0m6Rj77n8tzUWE88ks1hgN/NQ3/ZY+fMnWUTbbdio1yZeH5Kb9xenV/Q7p4/3a6dftCW/xLPyUl3tzZvZ/X6WWfyVWm32V7DM7A1oteIoYZkq+UQm8f2F0+LsD8j8mYjpslA4zbJgZ1jUcX5d6o2JdneU98HuqOh3PVHeB+O67SRjfH9JXZZKZLRajovtK3Kc4Pjp+3yj4Vo6eZLoSxXTeZ8LNNZpta9hKsts9Q2iUc2zbJlY71Asw7JtuuvSnt4Yy3ZtDWt1a8c2lgU7xXWw72Any8L7WrU6k+TbhLHN6hT3Gr76RpaVTazR6kBNKcsa63xaPbPWx7Jp5WKbdV7+eqh2i/sVa18HP5gu0c/SLW0sCre1a/WeVr5ez25xH7D/s36WtequkfpxGE4+eKcKAAAAAACgAJhUAQAAAAAAFACTKgAAAAAAgAJgUgUAAAAAAFAATKoAAAAAAAAKgEkVAAAAAABAATCpAgAAAAAAKAAmVQAAAAAAAAXApAoAAAAAAKAAmFQBAAAAAAAUAJMqAAAAAACAAliLfQDF0vvJDjJZHaRkZLZc31aVLMtURTmmbZvMZkPbotUWu8Qys9WesyYiskmevOtZJbdW211eltldTq12SPwpdrrEPpxuG8tcuqyi1MGySl27stTJsnK3WC/g4j9DwCX24Zf4/tw2cZ48dn7OPLrMkoiwzJQIarU52s8yJSIyJdTDskyoW6vlUIhlyd6oro7kzRKhBM9CSVGHkzyLpMR6sTTLImnRt2JZlWWJrOhnAxne52RFzVkb24ZN0rH13MG/ZbGYRC2ZTSyz69p2Q+axii1JFr5Vt26jXpuFZZKuTzq9vA86faLfOf28D0p+8fpwlpfy9cq9eTO736/VFn8ly8z+ClF7AyxTHCVarUo+nklif+E0P/MDsmjHDFkoLPpIsDPMsmBcXJd6Y/ya1R1N6uoUy3p07Xicr5fU9cmkIUslMlotx3m/luPi9ZFJxFiWlcXrI50Y0Gp2Xc3wYwTuhq+fRh6bjQKNdWy5r2GqVtvqG/hK1SLLlvH1DsXEc7ktws/9nl7x/O3aGmZZa8c2rQ52RlnWd7BTq8P7WlmWSfI+AWOX1elmbV99o1aXTaxhWaBGXFsb63wsm1kr2tPK+TbrdNf5aje/X7H2dYhGF+9n6ZY2rQ63tbNsT6tYr2c3vw/Y/5m4Z2g1XC+NYzGMXfr7DCKihhLRz6bUivG7YvqRcX4gnaZfvHyg4P3inSoAAAAAAIACYFIFAAAAAABQAEyqAAAAAAAACoBJFQAAAAAAQAEwqQIAAAAAACgAJlUAAAAAAAAFKPqkas2aNVRfX09Op5PmzJlD27ZtO+rjn3jiCTrttNNIkiSqra2lu+66i5LJ5FHXAQAAAAAAGC5FnVS9+OKL1NTURM3NzbRjxw6aOXMmzZ8/n7q7u3M+/le/+hUtXbqUmpubqa2tjX7+85/Tiy++SD/4wQ9G+MgBAAAAAACOKOqkavXq1XTjjTfSkiVLaPr06bR27VpyuVy0fv36nI//4IMP6Pzzz6err76a6uvr6eKLL6arrrrq7767BQAAAAAAMFyKNqmSZZm2b99O8+bNEwdjNtO8efNoy5YtOdc577zzaPv27dok6tNPP6UNGzbQpZdemnc/qVSKotEo+wcAAFBMGJsAAMYXa7F2HAwGKZvNUlVVFVteVVVFH330Uc51rr76agoGg/SlL32JVFWlTCZDN99881E//rdq1SpasWLFoOVyLEwmi33QcpPZrKstLLPYJa02W/m6+rZN8uRdzyq5WWZ3eXW1k2UOSTw9Thffn9Nt02qXIasodWh1pa4+0hb7KHfz9QK67QRcNpb5JdF22/hc3GMXbY8hsyQiWm1KBFlmjvZrtRLhmRLq0epMiH8cVA6FtDrZy29Ekr2RvFkilBBZiP8/vGRYtJORFF8vltbqSDrLslhWFY/LKiwbyIi2rKgs07eNmW6TxLd4dPozbzHxTDKLBXYzD/Vtj5U/f5JFtN2GjXpt4vUhuXl/cXpFv3P6eL92+kVb8ks8Ky/V1d68md3vZ5nFX6nVZn8Fy8zegFYrjhKWqZJPq7MS3180Lc7+gMyfiZguC4XTLAt2hkUdl1nWGxPt7ijvg91R0e96orwPxnXbScb4/pK6LJXIsEyOJ3V1hGWZREyrs3KCZenEgFYrGf4z6NvGTFWyujp371Wzcs7lJ5t8Y9PUp39JpaWl1BXjz+U23XVpT2+MZbu2hrW6tYN/ciPYKa6DfQc7WRbe16rVmSTfJoxtVqe41/DVN7KsbGKNVgdqSlnWWOfT6pm1PpZNKxfbrPPye4tqt7hfsfZ18IPpEv0s3dLGonBbu1bvaeXr9ewW9wH7P+tnWavuGqkfh2Hs099rNJTwfjalVozhFdP5WB9orNNqX8NUltnqG0SjWmTZsiPrRKNRoperh37Q/7+if1HF8di8eTOtXLmSnn76adqxYwe9/PLL9Prrr9NDDz2Ud51ly5ZRJBLR/h04cGAEjxgAAGAwjE0AAONL0d6pCgQCZLFY6PDhw2z54cOHqbo692zx/vvvp2uvvZZuuOEGIiKaMWMGxWIxuummm+jee+8ls3nwHNHhcJDD4Ri0HAAAoFgwNgEAjC9Fe6fKbrfT7NmzadOmTdoyRVFo06ZNNHfu3JzrxOPxQRMni+XIR5BUFW//AgAAAADAyCvaO1VERE1NTbR48WI655xz6Nxzz6UnnniCYrEYLVmyhIiIFi1aRBMnTqRVq1YREdGCBQto9erVdNZZZ9GcOXPok08+ofvvv58WLFigTa4AAAAAAABGUlEnVQsXLqSenh5avnw5dXV10axZs2jjxo3al1fs37+fvTN13333kclkovvuu48OHjxIFRUVtGDBAnr44YeL9SMAAAAAAMBJrqiTKiKi2267jW677bac2ebNm1nbarVSc3MzNTc3j8CRAQAAAAAA/H1j6tv/AAAAAAAARhtMqgAAAAAAAAqASRUAAAAAAEABMKkCAAAAAAAoACZVAAAAAAAABSj6t/8Vi93tI5PVQWarnS23SR6tttglllklt1jf5eXbczm12iHx0+p0iX043TaWuXRZRamDZZW6dmWpk2XlbrFewMV/hoBL7MMv8f25bWIe7bHzObVHl1kSEZaZEkGtNkf7WaZERKaEeliWCXVrtRwKsSzZG9XVkbxZIpTgWSgp6nCSZ5GUWC+WZlkkndXqWJb/sehEVtHqgYzCMllRc9bGtmGTxLeSn/E3GxaTqCWziWV2XdtuyDxWsSXJwrfq1m3Ua+N/003S9Umnl/dBp0/0O6ef90HJL14fzvJSvl65N29m9/u12uKvZJnZXyFqb4BliqNEq1XJx7KsJPYXTfMzPyCLdsyQhcKijwQ7wywLxmWt7o3JLOuOJnV1imU9unY8ztdL6vpk0pClEhmtluO8X8tx8frIJGIsy8ri9ZFODLBMycg5a2NbVbIsU5Vj670mwx9jN5lF39JfP/XXWTWTogxBPt/4+XayOt0U7Iyy5X0HO7U6vK+VZZkk7xMwdlmdbtb21TdqddnEGpYFasS1tbHOx7KZtaI9rZxvs053na928/sVa1+HaHTxfpZuadPqcFs7y/a0ivV6dvP7gP2fiXuGVsP10jgWw9ilv88gImrU3cNOnlTCsorpYqwPNNaxzNcwVatt9Q18J9Uiy5bx9Q7FxMiyLcL72Z5ecY3ctTWs1a0d24joxF1D8U4VAAAAAABAATCpAgAAAAAAKAAmVQAAAAAAAAXApAoAAAAAAKAAmFQBAAAAAAAUAJMqAAAAAACAAmBSBQAAAAAAUABMqgAAAAAAAAqASRUAAAAAAEABMKkCAAAAAAAogLXYB1As5V84m8x2F9ldTrbcIYlT4nTZWeZ027TaZcgqSh1aXamrj7TFPsrdfL2AbjsBl41lfkm03TY+//XYRdtjyCyJiFabEkGWmaP9Wq1EeKaEerQ6E+pmmRwKaXWyN8qyZG8kb5YIJUQWSvL1wqKdjKT4erG0VkfSWZbFsqp4XFZh2UBGtGVFZZm+bcx0myS+xaPTn3mLiWeSWSywm3mob3us/PmTLKLtNmzUa7OIx7l5f3F6Rb9z+ni/dvpFW/JLPCsv1dXevJnd72eZxV+p1WZ/BcvM3oBWK44SlqmST6uzEt9fNC3O/oDMn4mYLguF0ywLdoZFHZdZ1hsT7e4o74PdUdHveqK8D8Z120nG+P6SuiyVyLBMjid1dYRlmURMq7NygmXpxIBWKxn+M+jbxkxVsrr62HuvyWzW1RaWWeyij5it/Jqlb9skT971rJJbq+0u8TwrcpziW4/5ME86f3rhV2Sy2P/+A2HUsTpFn/fVN7KsbGKNVgdqSlnWWOfT6pm1PpZNKxfbrPPye4tqt7hfsfZ18IPpatXKdEsbi8Jt7Vq9p5Wv17Nb3Afs/6yfZa26a6R+HIaxT3+v0Wi4h508SYzhFdP5WB9orNNqX8NUltnqG0SjmmfZMrHeoRgfQ7fp7gf39MZYtmtrWKtbO7axLNgp7j/7DnayLLxPvB4ySb5NIiI1Kw9aNhR4pwoAAAAAAKAAmFQBAAAAAAAUAJMqAAAAAACAAmBSBQAAAAAAUABMqgAAAAAAAAqASRUAAAAAAEABMKkCAAAAAAAoACZVAAAAAAAABcCkCgAAAAAAoACYVAEAAAAAABQAkyoAAAAAAIACWIt9AMUybeYkskpuqih1sOWVunZlqZNl5W67VgdcdpYFXDat9ks2lrltYu7qsfN5bIkuMyciLDMlgiKL9rNMiYhMCfWwLBPq1mo5FGJZsjeqqyN5s0QowbNQUtThJM8iKbFeLM2ySDqr1bGsyrJEVtHqgYzCMllRc9bGtmGTxLeSn/G3CRaTqCWziWV2XdtuyDxWsSXJwrfq1m3Ua7OwTHKLPuL08j7o9Il+5/TzPij5JZGVl/L1yr15M7vfr9UWfyXLzP4KUXsDLFMcJVqtSj6WZSWxv2ian/kBWbRjhiwUFn0k2BlmWTAua3VvTGZZdzSpq1Ms69G143G+XlLXJ5OGLJXIaLUc5/1ajovXRyYRY1lWFq+PdGKAZUpGzlkb26qSZZmqHFvvNZnNhrboWxa7xDKz1Z6zJiKySZ6861klt1bbXV6W2V2iTzokPoQ4dddFp66Pu3TLM4kYHfglAYwaVqebtX31jVpdNrGGZYEacW1trPOxbGataE8r59us013nq938dWPt6xCNrlaWpVvatDrc1s6yPa1ivZ7d/D5g/2finqHVcL00jsUwdunvM4iIGnX3sJMnlbCsYroY6wONdSzzNUzValt9A99JtciyZXy9QzExhm6L8H62p1eMm7u2hlnW2rFNq4OdUZb1Hjig1ZH9bSzLJPlYPNrgnSoAAAAAAIACYFIFAAAAAABQAEyqAAAAAAAACoBJFQAAAAAAQAEwqQIAAAAAACgAJlUAAAAAAAAFKPqkas2aNVRfX09Op5PmzJlD27ZtO+rjw+Ew3XrrrTRhwgRyOBx06qmn0oYNG0boaAEAAAAAALii/p2qF198kZqammjt2rU0Z84ceuKJJ2j+/Pn08ccfU2Vl5aDHy7JMX/3qV6myspJ+85vf0MSJE6mjo4N8Pt/IHzwAAAAAAAAVeVK1evVquvHGG2nJkiVERLR27Vp6/fXXaf369bR06dJBj1+/fj319fXRBx98QDbbkT8sWV9fP5KHDAAAAAAAwBTt43+yLNP27dtp3rx54mDMZpo3bx5t2bIl5zqvvfYazZ07l2699VaqqqqixsZGWrlyJWWz2bz7SaVSFI1G2T8AAIBiwtgEADC+FO2dqmAwSNlslqqqqtjyqqoq+uijj3Ku8+mnn9I777xD11xzDW3YsIE++eQTuuWWWyidTlNzc3POdVatWkUrVqwYtPy7X6onl6eEAi4bW+6XRNtt43NOj120SwyZORHRalMiyLNov1YrEZ4poR6tzoS6WSaHQlqd7OUDbrI3kjdLhBIiCyX5emHRTkZSfL1YWqsjaT5RjWVV8biswrKBjGjLisoyfduY6TZJfItHpz/zFhPPJLNYYDfzUN/2WPnzJ1lE223YqNdmEY9z8/7i9DpE7XPyzC/akl/iWXmprvbmzex+P8ssfvGxWLO/gmVmb0CrFUcJy1TJp9VZie8vmhZnf0Dmz0RMl4XCaZYFO8Oijsss642JdneU98HuqOh3PVHeB+O67SRjfH9JXZZKZFgmx5O6OsKyTCKm1Vk5wbJ0YkCrlQz/GfRtY6YqWV197L3XZDbragvLLHbRR8xWO8v0bZvkybueVXKzzO7y6mrePx2SuPw7XXx/Tl0/dxmyilLR5yt19ZG22Ee5W6wX0G0jPtBP7xDkG5vg+Fidos/76htZVjaxRqsDNaUsa6zzafXMWh/LppWLbdZ5eR+vdovXjbWvgx9MV6tWplvaWBRua9fqPa18vZ7d4j5g/2f9LGvVXSP14zCMffp7jUbDtXTyJDGGV0znY32gsU6rfQ1TWWarbxCNap5ly8R6h2J8DN2mux/c0xtj2a6tYa1u7eDfexDsFPeffQc7WRbeJ14PmSTf5nhV9C+qOB6KolBlZSU988wzNHv2bFq4cCHde++9tHbt2rzrLFu2jCKRiPbvwIEDI3jEAAAAg2FsAgAYX4r2TlUgECCLxUKHDx9myw8fPkzV1dU515kwYQLZbDayWMRveBsaGqirq4tkWSa73T5oHYfDQQ6HY9ByAACAYsHYBAAwvhTtnSq73U6zZ8+mTZs2acsURaFNmzbR3Llzc65z/vnn0yeffEKK7uM2f/vb32jChAk5J1QAAAAAAADDragf/2tqaqJ169bRc889R21tbfS9732PYrGY9m2AixYtomXLlmmP/973vkd9fX10xx130N/+9jd6/fXXaeXKlXTrrbcW60cAAAAAAICT3JA//tfe3k6/+MUvqL29nZ588kmqrKyk3//+9zR58mQ644wzjmkbCxcupJ6eHlq+fDl1dXXRrFmzaOPGjdqXV+zfv5/Muv/YXVtbS2+88QbddddddOaZZ9LEiRPpjjvuoHvuuWeoPwYAAMBxOxFjIAAAjB9Deqfqv//7v2nGjBn0pz/9iV5++WUaGDjyDVq7du3K+y18+dx2223U0dFBqVSK/vSnP9GcOXO0bPPmzfTss8+yx8+dO5e2bt1KyWSS2tvb6Qc/+AH7P1YAAADD6USOgQAAMD4MaVK1dOlS+uEPf0hvvfUW+79MF110EW3duvWEHRwAAMBogzEQAACMhjSp+vDDD+lrX/vaoOWVlZUUDAZzrAEAADA+YAwEAACjIU2qfD4fHTp0aNDynTt30sSJEws+KAAAgNEKYyAAABgNaVL17W9/m+655x7q6uoik8lEiqLQ+++/T9///vdp0aJFJ/oYAQAARg2MgQAAYDSkSdXKlSvp9NNPp9raWhoYGKDp06fTP/7jP9J5551H991334k+RgAAgFEDYyAAABiZVFVVh7rygQMH6MMPP6SBgQE666yzaNq0aSfy2IZFNBolr9dLXV1dVFpaSpZEhOWmRFirzal+likR8Vl5JdTDsmyoW6vlUIhlyd6oro7kzRKhBM9CSVGHkzyLpMR6sTTLIumsVsey/OlNZMUfTh7IKCyTFTVnbWwbNkl8K/kZZ/AWk6jtZhPL9G1j5rGKLUkWvlW3bqNeG/9WSMlt02qn18Eyp88par+TZZJfEll5KV+v3Js3s/v9Wm3xV7LM7K8QtTfAMsVRotWq5GNZVhL7G0jzMz8gi3bMkIUSoo8E47y/BOOyVvfGZJZ1R5O6OsWyHl07HufrJXV9MmnIUomMVstx3q/luHh9ZBIxlmVl8fpIJwZYpmTknLWxrSpZlqnKsfVek9lsaIu+ZbbyPzyubxszm+TRaotdYplVcmu13eVlmd0l+qRD4n8Jw+kS+3Dq+jgRkUuXVZTyPl+pa1eW8j5f7hbrBVz8Zwi4xD78Et+f2ybOk8euq3XLo9EoVVdXUyQSodJS/poZa07kGPj52GSdcQ2ZLOP/j9lbnW7W9tU3anXZxBqWBWpEP2ms87FsZq1oTyvn26zTXeer3fx1Y+3rEI2udpal97VpdbiNZ8FWsV7Pbn4fsP8zcc/QarheGsdiGLv09xlERI26a+nkSSUsq5guxvpAYx3LfA1TtdpW38B3Ui2ybBlf71BMjKEdEd7P9vSKcXPXgTDLWjtEO9gZZVnfwU6tDu9rZVkmycfi8UjNypT58PmCx6Uh/50qoiN/N6q2traQTQAAAIxJGAMBAOBzQ/r43ze+8Q169NFHBy1/7LHH6Fvf+lbBBwUAADBaYQwEAACjIU2q3nvvPbr00ksHLb/kkkvovffeK/igAAAARiuMgQAAYDSkSdXAwAD7g4efs9lsFI1Gc6wBAAAwPmAMBAAAoyFNqmbMmEEvvvjioOUvvPACTZ8+veCDAgAAGK0wBgIAgNGQvqji/vvvp69//evU3t5OF110ERERbdq0if7zP/+TXnrppRN6gAAAAKMJxkAAADAa0qRqwYIF9Oqrr9LKlSvpN7/5DUmSRGeeeSa9/fbbdOGFF57oYwQAABg1MAYCAIDRkL9S/bLLLqPLLrvsRB4LAADAmIAxEAAA9Ib0f6oAAAAAAADgiCG9U5XNZunHP/4x/frXv6b9+/eTLMss7+vrOyEHBwAAMNpgDAQAAKMhvVO1YsUKWr16NS1cuJAikQg1NTXR17/+dTKbzfTAAw+c4EMEAAAYPTAGAgCAkUlVVfV4V5o6dSr95Cc/ocsuu4xKSkqopaVFW7Z161b61a9+NRzHekJEo1Hyer3U/dpaKnVLlA11s1wOhbQ62cv/3kiyN5I3S4QSIgsl+Xph0U5GUny9WFqrI+ksy2JZ8dQksgrLBjKiLSv8KdS3jZluk8S3eHT62bfFxDO72ZSzNrY9Vj6Hlyyi7TZs1GuziMe5bSxzeh2i9jl55hdtyS/xrLxUV3vzZna/n2UWf6VWm/0VLDN7A1qtOEpYpko+kUl8f/1pcfYHZP5MxHRZKJFmWTCe1tX8t+O9MdHujvI+2B0V/a4nyvtgXLedZIzvL6nLUokMy+R4UldHWJZJxLQ6KydYlk4MaLWS4T+Dvm3MVCWrq4+995rMZl1tYZnZas9ZG9s2ycMyi130LavkZpnd5dXVvH86JPHhAKeL78+p6+cuQ1ZRKvp8pa4+0hb7KHfz9QK67QRc/HXkl0TbbeOvTY9dtD2GzJIQz7UpEWaZOdWv1UokKOpQj1ZHYwmqvPxmikQiVFpaSmPRcIyBn49N1hnXkMky+G9gjSZWp+jzvvpGlpVNrNHqQA1/fhvrfFo9s9bHsmnlYpt1Xt7Hq93idWPt6+AH09Wulel9bSwKt4ks2MrX69kt+uT+z/pZ1qq7RurHYRj79PcajYZr6eRJYgyvmM7H+kBjnVb7GqayzFbfIBrVPMuWifUOxfgY2qG7H9zTG2PZrgNhrW7tCLMs2CnuP/sOdrIsvK9VqzNJvk3IT83KlPnw+YLHpSG9U9XV1UUzZswgIiKPx0ORyJFB9p//+Z/p9ddfH/LBAAAAjHYYAwEAwGhIk6pJkybRoUOHiOjIb+zefPNNIiL685//TA6H42irAgAAjGkYAwEAwGhIk6qvfe1rtGnTJiIiuv322+n++++nadOm0aJFi+i73/3uCT1AAACA0QRjIAAAGA3p2/8eeeQRrV64cCFNnjyZtmzZQtOmTaMFCxacsIMDAAAYbTAGAgCA0ZD/+K/e3Llzae7cuSdiUwAAAGMKxkAAABjypOrjjz+mp556itrajnzbTkNDA91+++102mmnnbCDAwAAGI0wBgIAgN6Q/k/Vb3/7W2psbKTt27fTzJkzaebMmbRjxw5qbGyk3/72tyf6GAEAAEYNjIEAAGA0pHeq7r77blq2bBk9+OCDbHlzczPdfffd9I1vfOOEHBwAAMBogzEQAACMhvRO1aFDh2jRokWDln/nO9/RvmYWAABgPMIYCAAARkOaVH35y1+mP/zhD4OW//GPf6QLLrig4IMCAAAYrTAGAgCA0ZA+/nf55ZfTPffcQ9u3b6d/+Id/ICKirVu30ksvvUQrVqyg1157jT12NOr49WvksdsoEUqw5clQUtThJM8iKa1OxNIsi6SzWh3LqixLZBWtHsgoLJMVNWdtbBs2SXwr+RlnzRaTqCWziWV2XdtuyDxWsSXJwrfq1m3Ua7OwTHLbtNrp5X8U0+lzitrvZJnkl0RWXsrXK/fmzex+v1Zb/JUsM/srRO0NsExxlGi1Kvl4Jon9hdP8zA/Ioh0zZKGw6CPBzjDLgnFZq3tjMsu6o0ldnWJZj64dj/P1kro+mTRkqURGq+U479dyPKLVmUSMZVlZvD7SiQGWKRk5Z21sq0qWZapybL3XZDYb2qJvWewSy8xWe86aiMgmefKuZ5XcWm13eVlmd4k+6ZD4pdLpEvtw6vo4EZFLl1WU8j5fqWtXlvI+X+4W6wVc/GcIuMQ+/BLfn9smzpPHzs+ZR5dZEhGWmRJBrTZH+1mmRESmhHpYlgl1a7UcCrEs2RvV1ZGcywdkfu0ci8bqGGh1ulnbV9+o1WUTa1gWqBHX1sY6H8tm1or2tHK+zTrddb7azV831r4O0ehqZVm6pU2rw23tLNvTKtbr2c374/7PRN9tNVwvjWMxjF36+wwiokbdtXTypBKWVUwXY32gsY5lvoapWm2rb+A7qRZZtoyvdygmxtBtEd7P9vSKcXPX1jDLWju2aXWwM8qyvoOdWh3ex18PmSQfi2FsGNKk6pZbbiEioqeffpqefvrpnBkRkclkomyW31ABAACMZRgDAQDAaEiTKuUYf9MMAAAw3mAMBAAAo+P6P1Vbtmyh3/3ud2zZ//2//5emTJlClZWVdNNNN1EqlcqzNgAAwNiFMRAAAPI5rknVgw8+SH/961+19ocffkjXX389zZs3j5YuXUr/9V//RatWrTrhBwkAAFBsGAMBACCf45pUtbS00P/6X/9La7/wwgs0Z84cWrduHTU1NdFPfvIT+vWvf33CDxIAAKDYMAYCAEA+xzWpCoVCVFVVpbX/+7//my655BKt/cUvfpEOHDhw4o4OAABglMAYCAAA+RzXpKqqqor27t1LRESyLNOOHTu0r5MlIurv7yebzZZv9bzWrFlD9fX15HQ6ac6cObRt27a/vxId+S2hyWSiK6644rj3CQAAcDyGawwEAICx77gmVZdeeiktXbqU/vCHP9CyZcvI5XKxP3T4P//zPzR16tSjbGGwF198kZqamqi5uZl27NhBM2fOpPnz51N3d/dR19u3bx99//vfxx9aBACAETEcYyAAAIwPxzWpeuihh8hqtdKFF15I69ato3Xr1pHdLv5I5fr16+niiy8+rgNYvXo13XjjjbRkyRKaPn06rV27llwuF61fvz7vOtlslq655hpasWIFnXLKKce1PwAAgKEYjjEQAADGh+P6O1WBQIDee+89ikQi5PF4yGKxsPyll14ij8dzzNuTZZm2b99Oy5Yt05aZzWaaN28ebdmyJe96Dz74IFVWVtL1119Pf/jDH466j1Qqxb7iNhqNHuXRAAAAuZ3IMRBjEwDA+HJc71R9zuv1DhpMiIjKysrYb+3+nmAwSNlslv3HX6Ijn1vv6urKuc4f//hH+vnPf07r1q07pn2sWrWKvF6v9q+2tvaYjw8AAMDoRIyBGJsAAMaX43qnqtj6+/vp2muvpXXr1lEgEDimdZYtW0ZNTU1aOxqNUm1tLX20sZ1cZgtF0ln2+FhW1epEVmHZQEa0ZUVlmb5tzHSbJL7Fo9PPeC0mnklmscBu5qG+7bHyebNkEW23YaNem7hJkNz8P1s7vQ5R+5w884u25Jd4Vl6qq715M7vfzzKLv1Krzf4Klpm94rlXHCUsUyWfyCS+v3BanP0BmT8TMV0WCqdZFuwMizous6w3Jtrd0STLuqPit9A9Uf4HQeO67SRjfH9JXZZKZFgmx5O6OsKyTCKm1Vk5wbJ0YkCrlQz/GfRtY6YqWV197L3XZDbran7zabGLPmK28htQfdsm8d/469ezSm6W2V1eXc37p0MSlzmni+/PqevnLkNWUSr6fKWuPtIW+yh38/UCuu0EXPx15JdE223jr02PXbQ9hsySEM+1KRFkmTnar9VKhGdKqEerMyH+/1TlUEirk738XZJkbyRvlgiJvpUM8T6fDIt2MiL6fELXx+MKv+aerPKNTflYnaLP++obWVY2sUarAzWlLGus82n1zFofy6aVi23WeXkfr3aL1421r4MfTFerVqZb2lgUbmvX6j2tfL2e3aI/7v+sn2WtumukfhyGsU9/r9FouJZOniTG8IrpfKwPNNZpta+B/19FW32DaFTzLFsm1jsU42PoNt11aU9vjGW7toa1urWDf2FasFNcB/sOdrIsvE+8HjJJvk04uRR1UhUIBMhisdDhw4fZ8sOHD1N1dfWgx7e3t9O+fftowYIF2jLl/7/Rs1qt9PHHHw/6T8IOh4McDv4iBgAAKCaMTQAA48uQPv53otjtdpo9ezZt2rRJW6YoCm3atInmzp076PGnn346ffjhh9TS0qL9u/zyy+krX/kKtbS04OMTAAAAAAAw4or+8b+mpiZavHgxnXPOOXTuuefSE088QbFYjJYsWUJERIsWLaKJEyfSqlWryOl0UmMj/9iDz+cjIhq0HAAAAAAAYCQUfVK1cOFC6unpoeXLl1NXVxfNmjWLNm7cqH15xf79+8lsLuobagAAAAAAAHkVfVJFRHTbbbfRbbfdljPbvHnzUdd99tlnT/wBAQAAAAAAHCO8BQQAAAAAAFAATKoAAAAAAAAKgEkVAAAAAABAATCpAgAAAAAAKAAmVQAAAAAAAAUYFd/+Vww7wklymMwkKypbrm9neUTKMW7bOFO1mEQtmU0ss+vadkPmsYotSRa+Vbduo16bhWWS26bVTq+DZU6fU9R+J8skvySy8lK+Xrk3b2b3+7Xa4q9kmdlfIWpvgGWKo0SrVcnHM0nsL5zmZ35AFu2YIQuF01od7AyzLBiXtbo3JrOsO5rU1SmW9eja8ThfLxkT+0saslQio9VyPMkyOR7R6kwixrKsnNDqdGKAZUpGzlkb26qSZZmqHFvvNRn+fIHJLPqWxS6xzGy156yJiGySJ+96Vsmt1XaXl2V2l+iTDolfnpwusQ+nro8TEbl0WUUp7/OVunZlKe/z5W6xXsDFf4aAS+zDL/H9uW3iPHns/Jx5dJklEWGZKRHUanO0n2VKRGRKqIdlmVC3VsuhEMuSvVFdHcmbJUIJnoVEn0yGef9MRkSfT+j6OBFRJC36VsxwkUxkRT8byIhaf11Nqcd6JT05zfn21WR1uqmxzseWz6wV7WnlbpbV6a7z1W7+urH2dYhGVyvL0i1tWh1ua2fZnlaxXs9u3h/3fyb6bqvhemnsEzB26e8ziIgaddfSyZNKWFYxXYz1gcY6lvkapmq1rb6B76RaZNkyvt6hmBhDt0V4P9vTK8bNXVvDLGvt2KbVwc4oy/oOdmp1eB9/PWSSfCyGscvq5NdI72TR78pra7U6UHPkfjaTjNH7Hz5f8H7xThUAAAAAAEABMKkCAAAAAAAoACZVAAAAAAAABcCkCgAAAAAAoACYVAEAAAAAABQAkyoAAAAAAIACYFIFAAAAAABQAEyqAAAAAAAACoBJFQAAAAAAQAEwqQIAAAAAACgAJlUAAAAAAAAFsBb7AIolnFHInmO5fpZpMfFMMosFdjMP9W2Plc9VJYtouw0b9dos4nFuG8ucXoeofU6e+UVb8ks8Ky/V1d68md3vZ5nFX6nVZn8Fy8zegFYrjhKWqZJPZBLfXzitaPWArLAspstC4TTLgp1hUcdllvXGRLs7mmRZdzSl1T26mogorttOMsb3l9RlqUSGZXI8qasjLMskYlqdlRMsSycGtFrJ8J9B3zZmqpLV1fycHY3JbNbVFpZZ7KKPmK285+vbNsmTdz2r5GaZ3eXV1bx/OiRxaXG6+P6cun7uMmQVpaLPV+rqI22xj3I3Xy+g207AxV9Hfkm03Tb+2vTYRdtjyCwJ8VybEkGWmaP9Wq1EeKaEerQ6E+pmmRwKaXWyN8qyZG8kb5YIib6VDPE+nwyLdjLC+3xC188j6SzLYllVPC7L+9lARrRlRWWZvm3MdJukfD1XzpsAEdFvr59NpaWlZO3r4EFXq1amW9pYFG5r1+o9rXy9nt2iP+7/rJ9lrbprpL4/wNinv9doNFxLJ08SY3jFdD7WBxrrtNrXMJVltvoG0ajmWbZMrHcoxsfQbbrr0p7eGMt2bQ1rdWvHNpYFO8V1sO9gJ8vC+8TrIZPk24SxzeoU9xq++kaWlU2s0epATSnLGut8Wj2z1seyaeVim3W6e+tq95F7lWg0StX3DvmQNXinCgAAAAAAoACYVAEAAAAAABQAkyoAAAAAAIACYFIFAAAAAABQAEyqAAAAAAAACoBJFQAAAAAAQAEwqQIAAAAAACgAJlUAAAAAAAAFwKQKAAAAAACgAJhUAQAAAAAAFMBa7AMoFp/VTA6TmexmE1vusYp5pmThc063RTzWa7OwTHLbtNrpdbDM6XOK2u9kmeSXRFZeytcr9+bN7H6/Vlv8lSwz+ytE7Q2wTHGUaLUq+Xgmif2F0wrLBmTRjhmyUDit1cHOMMuCcVmre2Myy7qjSV2dYlmPrh2P8/WSMbG/pCFLJTJaLceTLJPjEa3OJGIsy8oJrU4nBlimZOSctbGtKlmWqQo/T/mYzGZDW/Qti11imdlqz1kTEdkkT971rJJbq+0uL8vsLtEnHRK/JDhdYh9OXR8nInLpsopS3ucrde3KUt7ny91ivYCL/wwBl9iHX+L7c9vEefLY+Tnz6DJLIsIyUyKo1eZoP8uUiMiUUA/LMqFurZZDIZYle6O6OpI3S4QSPAuJPpkM8/6ZjIg+n9D1cSKiSFr0rVhWZVkiK/rZQIb3OVlRc9bGtmGTdGw9d/Bv5XSXSJJ011b9dTalElGGII/2W75DHpuNenbz/rj/M9F3Ww3XS2OfgLFLf59BRNSou5ZOnlTCsorpYqwPNNaxzNcwVatt9Q18J9Uiy5bx9Q7FxItzW4T3sz29YtzctTXMstaObVod7IyyrO9gp1aH97WyLJPkYzGMXVanm7V99Y1aXTaxhmWBGnFP21jnY9nMWtGeVs63Wae7v6528/sVa1+HaHTxfpZuadPqcFu7Vu9pPbLOQJqPuUOFd6oAAAAAAAAKgEkVAAAAAABAATCpAgAAAAAAKAAmVQAAAAAAAAXApAoAAAAAAKAAmFQBAAAAAAAUAJMqAAAAAACAAoyKSdWaNWuovr6enE4nzZkzh7Zt25b3sevWraMLLriA/H4/+f1+mjdv3lEfDwAAAAAAMJyKPql68cUXqampiZqbm2nHjh00c+ZMmj9/PnV3d+d8/ObNm+mqq66id999l7Zs2UK1tbV08cUX08GDB0f4yAEAAAAAAEbBpGr16tV044030pIlS2j69Om0du1acrlctH79+pyPf/755+mWW26hWbNm0emnn04/+9nPSFEU2rRp0wgfOQAAAAAAAJG1mDuXZZm2b99Oy5Yt05aZzWaaN28ebdmy5Zi2EY/HKZ1OU1lZWc48lUpRKpXS2tFotLCDBgAAKBDGJgCA8aWok6pgMEjZbJaqqqrY8qqqKvroo4+OaRv33HMP1dTU0Lx583Lmq1atohUrVgxafrbPSS6zhSS3jS13eh2i9jl55hdtyS/xrLxUV3vzZna/n2UWf6VWm/0VLDN7A1qtOEpYpko+kUl8f+G0otUDssKymC4LhdMsC3aGRR2XWdYbE+3uaJJl3VFxY9Cjq4mI4rrtJGN8f0ldlkpkWCbHk7o6wrJMIqbVWTnBsnRiQKuVDP8Z9G1jpipZXc3P2dGYzGZdbWGZxS76iNlqZ5m+bZM8edezSm6W2V1eXc37p0MSL2eni+/PqevnLkNWUSr6fKWuPtIW+yh38/UCuu0EXPx15JdE223jb4h77KLtMWSWhHiuTYkgy8zRfq1WIjxTQj1anQnxjw7LoZBWJ3v5jWuyN5I3S4RE30qGeJ9PhkU7GeF9PqHr55F0lmWxrCoel+X9bCAj2rKiskzfNma6TdKx91z+MQWLiWeSWSywm3mob3us/PmTLKLt1m3UaxOvjbiSJeJP30kp39j0s5c/JnvxP0QCQ6Dv842Ga+nkSWIMr5jOx/pAY51W+xqmssxW3yAa1TzLlon1DsX4GLpNd13a0xtj2a6tYa1u7eD/Jz3YKa6DfQc7WRbe16rVmSTfJoxtVqe41/DVN7KsbGKNVgdqSlnWWOfT6pm1PpZNKxfbrPPy10O1W9yvWPs6+MF0iX6WbmljUbitXav3tPL1enaL+4C9B/pZ1tYvXg/6cfhz8nGNnvmN6Sv3I488Qi+88AK98sor5HQ6cz5m2bJlFIlEtH8HDhwY4aMEAADgMDYBAIwvRX2nKhAIkMViocOHD7Plhw8fpurq6qOu+/jjj9MjjzxCb7/9Np155pl5H+dwOMjhcOTNAQAARhrGJgCA8aWo71TZ7XaaPXs2+5KJz790Yu7cuXnXe+yxx+ihhx6ijRs30jnnnDMShwoAAAAAAJBTUd+pIiJqamqixYsX0znnnEPnnnsuPfHEExSLxWjJkiVERLRo0SKaOHEirVq1ioiIHn30UVq+fDn96le/ovr6eurq6iIiIo/HQx6PJ+9+AAAAAAAAhkPRJ1ULFy6knp4eWr58OXV1ddGsWbNo48aN2pdX7N+/n8y6LwP46U9/SrIs0ze/+U22nebmZnrggQdG8tABAAAAAACKP6kiIrrtttvotttuy5lt3ryZtfft2zf8BwQAAAAAAHCMxvS3/wEAAAAAABQbJlUAAAAAAAAFwKQKAAAAAACgAJhUAQAAAAAAFGBUfFFFMZz+T1PJY7eRs7yULXeWe3U1z+x+v1Zb/JUsM/srRO0NsExxlGi1Kvl4Jon9hdMKywZk0Y4ZslA4rdXBzjDLgnFZq3tjMsu6o0ldnWJZj64dj/P1kjGxv6QhSyUyWi3HkyyT4xGtziRiLMvKCa1OJwZYpmTknLWxrSpZlqkKP0/5mMxmQ9ui1Ra7xDKz1Z6zJiKySeJr/I3rWSW3VttdXpbZXU6tdkj8Zeh0iX043TaWuXRZRSn/w6GVunZlqZNl5W6xXsDFf4aAS+zDL/H9uW3iPHns/Jx5dJklEWGZKRHUanO0n2VKRGRKqIdlmVC3VsuhEMuSvVFdHcmbJUIJnoVEn0yGef9MRkSfT+j6OBFRJC36ViyrsiyRFf1sIMP7nKyoOWtj27BJOraeO/g3YRaTqCWziWV2XdtuyDxWsSXJwrfq1m3Ua7OwTNL1SaeX90GnT/Q7p1/Ukl+8NgbkNNEvPyaA0ULf34mIGnXX0smTSlhWMV2M9YHGOpb5GqZqta2+ge+kWmTZMr7eoZgYQ7dF+Li8p1eMm7u2hlnW2rFNq4OdUZb1HezU6vC+VpZlknwshrHL6nSztq++UavLJtawLFAj7mkb63wsm1kr2tPK+TbrdNf5aje/X7H2dYhGF+9n6ZY2rQ63tbNsT6tYr2c3vw/Ye0DcM7T189eDcSwebfBOFQAAAAAAQAEwqQIAAAAAACgAJlUAAAAAAAAFwKQKAAAAAACgAJhUAQAAAAAAFACTKgAAAAAAgAJgUgUAAAAAAFAATKoAAAAAAAAKgEkVAAAAAABAATCpAgAAAAAAKAAmVQAAAAAAAAWwFvsAiqXuysup1C2R2V/Blpu9Aa1WHCUsUyWfyCQvy8JpRasHZIVlMV0WCqdZFuwMizous6w3Jtrd0STLuqMpre7R1UREcd12kjG+v6QuSyUyLJPjSV0dYVkmEdPqrJxgWToxoNVKhv8M+rYxU5Wsrubn7GhMZrOutrDMYpe02my1s0zftkmevOtZJTfL7C6vrnayzCGJl5DTxffndNu02mXIKkodWl2pq4+0xT7K3Xy9gG47AZeNZX5JtN02/vsSj120PYbMkhDPtSkRZJk52q/VSoRnSqhHqzOhbpbJoZBWJ3ujLEv2RvJmiZDoW8kQ7/PJsGgnI7zPJ3T9PJLOsiyWVcXjsryfDWREW1ZUlunbxky3STr2nst/i2Ux8UwyiwV2Mw/1bY+VP3+SRbTdho16beL1Ibl5f3F6Rb9z+ni/dvpFW/JLPCsv1dXevJnd79dqi79Sq6OxBNEvNxDAiaDv842Ga+nkSWIMr5jOx/pAY51W+xqmssxW3yAa1TzLlon1DsX4GLpNd13a0xtj2a6tYa1u7djGsmCnuA72HexkWXhfq1ZnknybMLZZneJew1ffyLKyiTVaHagpZVljnU+rZ9b6WDatXGyzzstfD9Vucb9i7evgB9Ml+lm6pY1F4bZ2rd7Tytfr2S3uA/Ye6GdZW794PejH4fEM71QBAAAAAAAUAJMqAAAAAACAAmBSBQAAAAAAUABMqgAAAAAAAAqASRUAAAAAAEABMKkCAAAAAAAoACZVAAAAAAAABcCkCgAAAAAAoACYVAEAAAAAABQAkyoAAAAAAIACWIt9AMWinnsFqaWlFE4rbPmALNoxQxYKp7U62BlmWTAua3VvTGZZdzSpq1Ms69G143G+XjIm9pc0ZKlERqvleJJlcjyi1ZlEjGVZOaHV6cQAy5SMnLM2tlUlyzJV4ecpH5PZbGhbtNpil1hmttpz1kRENsmTdz2r5NZqu8vLMrvLqdUOiXd9p0vsw+m2scylyypKHSyr1LUrS50sK3eL9QIu/jMEXGIffonvz20T58lj5+fMo8ssiQjLTImgVpuj/SxTIiJTQj0sy4S6tVoOhViW7I3q6kjeLBFK8Cwk+mQyzPtnMiL6fELXx4mIImnRt2JZlWWJrOhnAxne52RFzVkb24ZN0rH13MG/fbKYRC2ZTSyz69p2Q+axii1JFr5Vt26jXpuFZZKuTzq9vA86faLfOf28D0p+8fpwlpfy9cq9eTO736/VFn8ly8z+ClF7AyxTHCVarUo+sVwS+1KjUSK6meDkpu/vRESNumvp5EklLKuYLvpcoLGOZb6GqVptq2/gO6kWWbaMr3coJsbQbRE+Lu/pFePmrq1hlrV2bNPqYGeUZX0HO7U6vK+VZZkkH4th7LI63aztq2/U6rKJNSwL1Ihra2Odj2Uza0V7WjnfZp3uOl/t5vcr1r4O0eji/Szd0qbV4bZ2lu1pFev17Ob3Afs/E/cMrYb7VONYDPnhnSoAAAAAAIACYFIFAAAAAABQAEyqAAAAAAAACoBJFQAAAAAAQAEwqQIAAAAAACgAJlUAAAAAAAAFwKQKAAAAAACgAKNiUrVmzRqqr68np9NJc+bMoW3bth318S+99BKdfvrp5HQ6acaMGbRhw4YROlIAAAAAAACu6JOqF198kZqamqi5uZl27NhBM2fOpPnz51N3d3fOx3/wwQd01VVX0fXXX087d+6kK664gq644gpqbW3N+XgAAAAAAIDhVPRJ1erVq+nGG2+kJUuW0PTp02nt2rXkcrlo/fr1OR//5JNP0j/90z/Rv/3bv1FDQwM99NBDdPbZZ9O///u/j/CRAwAAAAAAEFmLuXNZlmn79u20bNkybZnZbKZ58+bRli1bcq6zZcsWampqYsvmz59Pr776as7Hp1IpSqVSWjsSiRARUX9/PxERDaQV9nh9O27IYomMyBJpliXislYn4zxLxZNanU6kWJZJiPUySZll2VQ6Z01ElJXFsShykmVKOqGreabq2mqGH4uakXU135+a1WVKlmcKP095qYY5vGoRpfGhuiWqIVXTotsqJhPLFKvYhyLbWJa1iuPMmnnXz1rEz5ux8PUyJpGlbYbn1iraSUuGZQlVbCeu2FkWy4rMluHHotjEz6DazHkzS6KfZabkgFabUwMsUwbioo4lWJbV9U/Z0D9Tuj6ZTPH+mZLFz55I8589lRHtZJb3l5Su/yQMfSmh6jKVP+9JXTul8j4n6zLZsF5a1+Z7IzrGnjvot08W3S5UlfdBfduYWVX2QL4P3WNthgNTFXEEWcP5zOrOdcbwPGR1z1Ha8Pyldc+t8Xm3O0SfsNh5fzHbRF8ymw39TNb97Gnx+lbSYvnn115VNb7qTy75xib5mHvl2GY1vDb0r/2YoY9LadGPHYZx0qy7ftl01zkiIuoX/TNrjfIoLl4rsX6+zcRATKvlOO/jmaTIsim+P/3Yqx9PifgYCmObmuH3CPrn3dgnMklxHZTjfKxPDIgsZuev+36zuDa4snw9q65fk6HPp3Wvh37Da2VA9zoyvsb0r79B4+ugO7Tx5/PrbqHjUlEnVcFgkLLZLFVVVbHlVVVV9NFHH+Vcp6urK+fju7q6cj5+1apVtGLFikHLvzBt2hCPGgBAx3gPfLR74qHeVwWHuN4o1t/fT16vt9iHUTT5xqbn6WARjqYIjK+TcJ6aiGj3sB4JwHHJGNrJv4j68FHWe384DgZOqELHpaJOqkbCsmXL2Dtb4XCY6urqaP/+/Sf1gG4UjUaptraWDhw4QKWlpcU+nFED5yU3nJfBcE5yM54XVVWpv7+fampqin1oRYWx6djgdZUbzstgOCe54bzkpj8vJSUlJ2RcKuqkKhAIkMViocOH+dz+8OHDVF1dnXOd6urq43q8w+Egh8MxaLnX60XnyqG0tBTnJQecl9xwXgbDOclNf14wacDYdLzwusoN52UwnJPccF5y+/y8nIhxqahfVGG322n27Nm0adMmbZmiKLRp0yaaO3duznXmzp3LHk9E9NZbb+V9PAAAAAAAwHAq+sf/mpqaaPHixXTOOefQueeeS0888QTFYjFasmQJEREtWrSIJk6cSKtWrSIiojvuuIMuvPBC+tGPfkSXXXYZvfDCC/SXv/yFnnnmmWL+GAAAAAAAcJIq+qRq4cKF1NPTQ8uXL6euri6aNWsWbdy4Ufsyiv3795PZLN5QO++88+hXv/oV3XffffSDH/yApk2bRq+++io1NjYe0/4cDgc1Nzfn/NjFyQznJTecl9xwXgbDOckN5+XY4DzlhvOSG87LYDgnueG85DYc58WknuzfawsAAAAAAFCAov/xXwAAAAAAgLEMkyoAAAAAAIACYFIFAAAAAABQAEyqAAAAAAAACjAuJ1Vr1qyh+vp6cjqdNGfOHNq2bdtRH//SSy/R6aefTk6nk2bMmEEbNmwYoSMdWcdzXtatW0cXXHAB+f1+8vv9NG/evL97Hseq4+0vn3vhhRfIZDLRFVdcMbwHWATHe07C4TDdeuutNGHCBHI4HHTqqaeOy9fR8Z6XJ554gk477TSSJIlqa2vprrvuomQyOUJHOzLee+89WrBgAdXU1JDJZKJXX331766zefNmOvvss8nhcNAXvvAFevbZZ4f9OEcDjE25YWwaDONSbhibcsPYxBVtXFLHmRdeeEG12+3q+vXr1b/+9a/qjTfeqPp8PvXw4cM5H//++++rFotFfeyxx9Tdu3er9913n2qz2dQPP/xwhI98eB3vebn66qvVNWvWqDt37lTb2trU6667TvV6vepnn302wkc+vI73vHxu79696sSJE9ULLrhA/Zd/+ZeROdgRcrznJJVKqeecc4566aWXqn/84x/VvXv3qps3b1ZbWlpG+MiH1/Gel+eff151OBzq888/r+7du1d944031AkTJqh33XXXCB/58NqwYYN67733qi+//LJKROorr7xy1Md/+umnqsvlUpuamtTdu3erTz31lGqxWNSNGzeOzAEXCcam3DA2DYZxKTeMTblhbBqsWOPSuJtUnXvuueqtt96qtbPZrFpTU6OuWrUq5+OvvPJK9bLLLmPL5syZo/7v//2/h/U4R9rxnhejTCajlpSUqM8999xwHWJRDOW8ZDIZ9bzzzlN/9rOfqYsXLx53g9fxnpOf/vSn6imnnKLKsjxSh1gUx3tebr31VvWiiy5iy5qamtTzzz9/WI+zmI5l8Lr77rvVM844gy1buHChOn/+/GE8suLD2JQbxqbBMC7lhrEpN4xNRzeS49K4+vifLMu0fft2mjdvnrbMbDbTvHnzaMuWLTnX2bJlC3s8EdH8+fPzPn4sGsp5MYrH45ROp6msrGy4DnPEDfW8PPjgg1RZWUnXX3/9SBzmiBrKOXnttddo7ty5dOutt1JVVRU1NjbSypUrKZvNjtRhD7uhnJfzzjuPtm/frn0M49NPP6UNGzbQpZdeOiLHPFqdDNdcI4xNuWFsGgzjUm4Ym3LD2HRinKjrrfVEHlSxBYNBymazVFVVxZZXVVXRRx99lHOdrq6unI/v6uoatuMcaUM5L0b33HMP1dTUDOp0Y9lQzssf//hH+vnPf04tLS0jcIQjbyjn5NNPP6V33nmHrrnmGtqwYQN98skndMstt1A6nabm5uaROOxhN5TzcvXVV1MwGKQvfelLpKoqZTIZuvnmm+kHP/jBSBzyqJXvmhuNRimRSJAkSUU6suGDsSk3jE2DYVzKDWNTbhibTowTNS6Nq3eqYHg88sgj9MILL9Arr7xCTqez2IdTNP39/XTttdfSunXrKBAIFPtwRg1FUaiyspKeeeYZmj17Ni1cuJDuvfdeWrt2bbEPrag2b95MK1eupKeffpp27NhBL7/8Mr3++uv00EMPFfvQAMYFjE0Yl44GY1NuGJuGz7h6pyoQCJDFYqHDhw+z5YcPH6bq6uqc61RXVx/X48eioZyXzz3++OP0yCOP0Ntvv01nnnnmcB7miDve89Le3k779u2jBQsWaMsURSEiIqvVSh9//DFNnTp1eA96mA2lr0yYMIFsNhtZLBZtWUNDA3V1dZEsy2S324f1mEfCUM7L/fffT9deey3dcMMNREQ0Y8YMisVidNNNN9G9995LZvPJ+TutfNfc0tLScfkuFRHGpnwwNg2GcSk3jE25YWw6MU7UuDSuzpzdbqfZs2fTpk2btGWKotCmTZto7ty5OdeZO3cuezwR0VtvvZX38WPRUM4LEdFjjz1GDz30EG3cuJHOOeeckTjUEXW85+X000+nDz/8kFpaWrR/l19+OX3lK1+hlpYWqq2tHcnDHxZD6Svnn38+ffLJJ9pATkT0t7/9jSZMmDAuBi2ioZ2XeDw+aHD6fHA/8n9nT04nwzXXCGNTbhibBsO4lBvGptwwNp0YJ+x6e1xfazEGvPDCC6rD4VCfffZZdffu3epNN92k+nw+taurS1VVVb322mvVpUuXao9///33VavVqj7++ONqW1ub2tzcPG6/tvZ4zssjjzyi2u129Te/+Y166NAh7V9/f3+xfoRhcbznxWg8fsvS8Z6T/fv3qyUlJeptt92mfvzxx+rvfvc7tbKyUv3hD39YrB9hWBzveWlublZLSkrU//zP/1Q//fRT9c0331SnTp2qXnnllcX6EYZFf3+/unPnTnXnzp0qEamrV69Wd+7cqXZ0dKiqqqpLly5Vr732Wu3xn3917b/927+pbW1t6po1a06ar1TH2DQYxqbBMC7lhrEpN4xNgxVrXBp3kypVVdWnnnpKnTx5smq329Vzzz1X3bp1q5ZdeOGF6uLFi9njf/3rX6unnnqqarfb1TPOOEN9/fXXR/iIR8bxnJe6ujqViAb9a25uHvkDH2bH21/0xuvgdbzn5IMPPlDnzJmjOhwO9ZRTTlEffvhhNZPJjPBRD7/jOS/pdFp94IEH1KlTp6pOp1Otra1Vb7nlFjUUCo38gQ+jd999N+e14vNzsXjxYvXCCy8ctM6sWbNUu92unnLKKeovfvGLET/uYsDYlBvGpsEwLuWGsSk3jE1cscYlk6qepO/1AQAAAAAAnADj6v9UAQAAAAAAjDRMqgAAAAAAAAqASRUAAAAAAEABMKkCAAAAAAAoACZVAAAAAAAABcCkCgAAAAAAoACYVAEAAAAAABQAkyoAAAAAAIACYFIFMMpcd911dMUVVxT7MAAAADQYmwCOzlrsAwA4mZhMpqPmzc3N9OSTT5KqqiN0RAAAcLLD2ARQOJOKVwjAiOnq6tLqF198kZYvX04ff/yxtszj8ZDH4ynGoQEAwEkKYxNA4fDxP4ARVF1drf3zer1kMpnYMo/HM+gjFl/+8pfp9ttvpzvvvJP8fj9VVVXRunXrKBaL0ZIlS6ikpIS+8IUv0O9//3u2r9bWVrrkkkvI4/FQVVUVXXvttRQMBkf4JwYAgNEOYxNA4TCpAhgDnnvuOQoEArRt2za6/fbb6Xvf+x5961vfovPOO4927NhBF198MV177bUUj8eJiCgcDtNFF11EZ511Fv3lL3+hjRs30uHDh+nKK68s8k8CAADjBcYmAAGTKoAxYObMmXTffffRtGnTaNmyZeR0OikQCNCNN95I06ZNo+XLl1Nvby/9z//8DxER/fu//zudddZZtHLlSjr99NPprLPOovXr19O7775Lf/vb34r80wAAwHiAsQlAwBdVAIwBZ555plZbLBYqLy+nGTNmaMuqqqqIiKi7u5uIiHbt2kXvvvtuzs/At7e306mnnjrMRwwAAOMdxiYAAZMqgDHAZrOxtslkYss+/+YmRVGIiGhgYIAWLFhAjz766KBtTZgwYRiPFAAAThYYmwAETKoAxqGzzz6bfvvb31J9fT1ZrXiZAwBA8WFsgvEM/6cKYBy69dZbqa+vj6666ir685//TO3t7fTGG2/QkiVLKJvNFvvwAADgJISxCcYzTKoAxqGamhp6//33KZvN0sUXX0wzZsygO++8k3w+H5nNeNkDAMDIw9gE4xn++C8AAAAAAEAB8GsBAAAAAACAAmBSBQAAAAAAUABMqgAAAAAAAAqASRUAAAAAAEABMKkCAAAAAAAoACZVAAAAAAAABcCkCgAAAAAAoACYVAEAAAAAABQAkyoAAAAAAIACYFIFAAAAAABQAEyqAAAAAAAACvD/Ac4pC90/dpsiAAAAAElFTkSuQmCC", - "text/plain": [ - "
" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], - "source": [ - "ex.plot_multiple_spatio_temporal(\n", - " jnp.stack([trj_first, trj_second])[:, :, 0, :],\n", - " grid=(1, 2),\n", - " figsize=(10, 3),\n", - " domain_extent=DOMAIN_EXTENT,\n", - " dt=DT,\n", - " include_init=True,\n", - " titles=[\"First\", \"Second\"],\n", - ");" - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "p = ex.plot" - ] - }, - { - "cell_type": "code", - "execution_count": 8, - "metadata": {}, - "outputs": [ - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjgAAAGiCAYAAADqYLxOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/SrBM8AAAACXBIWXMAAA9hAAAPYQGoP6dpAACHXklEQVR4nO3dd3xUVfr48c+0THrvjdCLlECACGJDOqugrooNQcS1sLsu7q6yu8qqu7qWdf3p+rUDsgooFkRFJKLY6L2FXpKQBgnpJJnM3N8fkxmIpGdm7szkeb9eeWUyc+fMc3KnPHPvc87RKIqiIIQQQgjhRbRqByCEEEII4WiS4AghhBDC60iCI4QQQgivIwmOEEIIIbyOJDhCCCGE8DqS4AghhBDC60iCI4QQQgivIwmOEEIIIbyOJDhCCCGE8DqS4AghhBDC6zg1wfnhhx+49tpriY+PR6PRsGLFihbvs27dOoYMGYLRaKRHjx4sWrToom1effVVUlJS8PX1JT09nc2bNzs+eCGEEEJ4LKcmOJWVlQwaNIhXX321VdsfP36cyZMnc/XVV7Nz504eeugh7rnnHr7++mv7Nh988AFz585l/vz5bN++nUGDBjF+/HgKCwud1Q0hhBBCeBiNqxbb1Gg0fPrpp0ydOrXJbR555BG+/PJL9u7da79u2rRplJSUsHr1agDS09MZNmwY//3vfwGwWCwkJSXx29/+lkcffdSpfRBCCCGEZ9CrHcCFNmzYwJgxYxpcN378eB566CEAamtr2bZtG/PmzbPfrtVqGTNmDBs2bGiy3ZqaGmpqaux/WywWiouLiYiIQKPROLYTQgghhHAKRVEoLy8nPj4erbb5k1BuleDk5+cTExPT4LqYmBjKyso4d+4cZ8+exWw2N7rNgQMHmmz3mWee4YknnnBKzEIIIYRwrezsbBITE5vdxq0SHGeZN28ec+fOtf9dWlpKcnIyx48fJygoyKGPZTKZ+O6777j66qsxGAwObdsdOKp/mmPfof/kbpSwbtTNWuvACJtntij8fLSIT3ac4sfDRVjqT9CG+Oq5dlAcUwbEcHTXRoK7D+az3QV8e/A0JrN1I3+jlkmXxHLjkAT6xDr2eeNK8hz1bN7eP3DPPuo+/y3ag19gvuxhLCPmdKgtd+yfozmrj+Xl5XTt2rVVn91uleDExsZSUFDQ4LqCggKCg4Px8/NDp9Oh0+ka3SY2NrbJdo1GI0aj8aLrw8PDCQ4Odkzw9UwmE/7+/kRERHjlE9dh/fMbDV9qoOo4+GnAP9xxQTYiv7SaD7dm88GWbE6VnLNe6ePPpV3DuT09mfGXxOJr0GEymcg/4s/4Id35VXofiipq+GhbDks3Z3GiqIpP9pXwyb4SBiWFcvvwZH41KA5/H7d6GbVInqOezdv7B27ax5I9YNRAvyshIqJDTbll/xzMWX20tdWa8hK3emceMWIEq1atanBdRkYGI0aMAMDHx4e0tDTWrl1rL1a2WCysXbuWOXM6llELF/MPh/DuUHwUTm2HnmNavk8bmS0KPxw+zZJNWXx7oBBz/eGaUH8DNw5J5NbhSfSIbv5bQESgkd9c2Z3Zl3dj47Ei3t+cxZp9+ezKLmFXdglPfbGfqYMTuC09mb5xjk2WhRBuorwASrIADcQPUTsa0UpOTXAqKio4cuSI/e/jx4+zc+dOwsPDSU5OZt68eZw6dYrFixcDcN999/Hf//6XP//5z9x99918++23fPjhh3z55Zf2NubOnctdd93F0KFDGT58OC+99BKVlZXMnDnTmV0RzpA4zJrg5GxxaIJTUFbNh1uyWXbh0RpgeEo4t6YnMbF/HL4GXZva1Go1jOwRycgekZy54KjOyaIq/rfxJP/beJLUpFBuS0/m2oHx+Pm0rX0hhBs7tdX6O7ov+MoXGU/h1ARn69atXH311fa/bXUwd911F4sWLSIvL4+srCz77V27duXLL7/kD3/4A//v//0/EhMTefvttxk/frx9m1tuuYXTp0/z+OOPk5+fT2pqKqtXr76o8Fh4gMShsHuZNcHpIMsFR2vWXnC0JsTPwA1DErhteDI9YxxTNxMZaOS+K7tz7+XdWH+0iKWbs/h6Xz47s0vYWX9U5/r6ozp9YuXNUAiPZ3uPShyqbhyiTZya4Fx11VU0N81OY7MUX3XVVezYsaPZdufMmSOnpLxB4jDr71NbwWKBFob8NaawrJrl9UdTcs6eP1oztEsYt6UnM2lA24/WtJZWq2FUz0hG9YzkdHkNy7dls2xzNlnFVSzecJLFG04yJDmUW4cn8ys5qiOE58qpP4Jje88SHsGtanBEJxNzCej9oLoUio5AVK9W33XjsSIW/XyCbzILqKs/WhPsq+eGIYnclp5MLwcdrWmtqCAjD1zVg/uu6M7PR8+wZFMWGfsL2J5VwvYs61GdG4YkMvOyFLpEBLg0NiFEB5jrrHWCIAmOh5EER6hHZ4D4wZC13noIuJUJzifbc5j74S7732ldwrhtuPVojdpHSbRaDZf3jOLynlEUllezfGsOy7ZkkV18jkXrT/Dx9hw+uX+kw06XCSGc7HQmmCrBGAyRvdWORrSBrCYu1GU7p93KOpyNx4p45OPdAExJjefrh67g4/tHcmNaourJzS9FB/ny4NU9+P6PV7P47uEMTAyhvLqOmYu2cLq8puUGhBDqs703JQxp12l0oR7ZW0JdtkO+tnPczTh6uoLf/G8bJrPC5AFx/OfmVHp7wIR7Wq2GK3pFsWjmcFIi/Mk5e457Fm/lXK1Z7dCEEC2R+huPJQmOUJftTaNwH9RUNLlZUUUNMxduofScicHJofz75kFotZ61jlh4gA8LZw4n1N/AruwS5n64E4vFJWvdCiHayz6CShIcTyMJjlBXcBwEJ4JigdzGR89Vm8zc+79tZBVXkRTux1vThzptZJSzdY0M4I070jDoNHy1N5/nvj6odkhCiKacOwtnDlkvJ8gQcU8jCY5QXzN1OBaLwp8+2s22k2cJ8tWzcMYwIgMvXnbDk6R3i+C5Xw8E4PXvj7J0c1YL9xBCqOLUNuvv8G4Q0LHlGYTrSYIj1NdMHc5/vjnE57ty0Ws1vHFHWotLK3iK6wcn8tCYngD8bcVefjx8WuWIhBAXkfobjyYJjlCfPcHZAhdMDLl8azavfGtd6uPpGwYwskekGtE5ze+v6cn1gxMwWxQeeG87B/PL1Q5JCHEhqb/xaJLgCPXFDQStASoL6xe0g/VHzjDvkz0AzLm6BzcPTVIzQqfQaDT868YBDO8aTnlNHXcv2kJhebXaYQkhwDq7uv0IjtTfeCJJcIT6DH4QO8B6OWcLRwrLue+9bdRZFK4dFM/csa2f4djTGPU63rwzjW6RAZwqOcc978rwcSHcQvFRqC4BvS/E9Fc7GtEOkuAI91B/CPjc8U3MXLSFsuo60rqE8fyvB3rccPC2CvX3YcGMYYT5G9idU8pDH+yQ4eNCqM12eip+sHXWdeFxJMER7qE+wcna/T3ZxedIDvfnzTvTPHY4eFulRAbw5vSh+Oi0fL2vgH+tPqB2SEJ0brKCuMeTBEe4BUv9HBMppqNE+sLCmcOI8PDh4G01LCWc52+yDh9/84djvLfxpMoRCdGJSYGxx5MER7iF5zdXc0YJxqipY+EEI92jAtUOSRVTUhN4uL7maP7Kfaw7WKhyREJ0QrWVULDPelkSHI8lCY5Q3Qdbsnjt+2PssPQAYIBySOWI1DVndA9uHJKI2aIwZ8kOMvPK1A5JiM4ld4d1dvXgBAiOVzsa0U6S4AhV/XT4DH/9dC8Afl0vtV7ZypXFvZVGo+GZGwZwabdwKuqHjxeUyfBxIVxG6m+8giQ4QjWHCsq5v344+JTUeC67arz1hlasLO7tfPRa3rhjKN2iAsgrrWbWu1uoqq1TOywhOgeZwdgrSIIjVHG63Lo6eHlNHcNSwnju1wPRJKQBGijNgvJ8tUNUXYi/gUUzhhMe4MPeU2X8ftlOzDJ8XAjnUhQpMPYSkuAIlztXa+aexVs5VXKOlAh/3rhzKEa9DoxBEN3PupEcxQEgOcKft6an4aPXkrG/gKdXZaodkhDerTQbKgpAq4e4QWpHIzpAEhzhUhaLwh8+2Mmu7BJC/Q0snGk9QmHXzMrinVVal3D+fZP1jfadn46zeMMJdQMSwpvZ3ntiB1hnWRceSxIc4VLPrj7A6n35+Oi0vHnnULpGBjTcoJmVxTuzawfF86fxvQH4+8p9fHdAho8L4RRSf+M1JMERLvP+ppO88cMxAJ779UCGdw2/eCPbm0rudjBLUe2FHriqOzcPTcSiwJwl29mXW6p2SEJ4H6m/8RqS4AiX+P7QaR7/zDpx1h/G9GLq4ITGN4zsBcZgMFVB4X4XRuj+NBoN/7x+AJf1iKCy1sysRVvJL5Xh40I4TF0N5O2yXpYh4h5PEhzhdAfzy3nw/e2YLQo3DE7gd9f0aHpjrRYS0qyXpQ7nIgadlv+7PY0e0YHkl1mHj1fWyJEuIRwifw+Ya8E/AsK6qh2N6CBJcIRTFZZVc/eiLVTU1DG8azjP3DgAjaaF1cGlDqdZIX4GFs4YRkSAD/tyy/jd0h0yfFwIR7jw9FRL71PC7UmCI5ymqrbOPhy8W2QAb96ZZh0O3hJ7giNHcJqSFO7PW3cNxajXsvZAIU99IafzhOgwmcHYq0iCI5zmb5/uZXdOKWH+BhbMGEaov0/Ld4Lzby5Fh6Gq2HkBerghyWH855ZUABatP8HH23LUDUgITycFxl5FEhzhFDuzS/hkxyk0GnjjzqGk/HI4eHP8wyG8u/Xyqe3OCdBLTBoQxx/GWFcff3b1AVnOQYj2Ki+AkixAA/FD1I5GOIAkOMLhFEXh6S+tM+5ePzih8eHgLZHTVK1231XdSAr3o7C8hrd+OK52OEJ4plP1NX/RfcE3WN1YhENIgiMcLmN/AZtPFGPUa/njuN7ta0RmNG41o17Hn8f3AeCNH45SWC5Dx4VoM6m/8TqS4AiHMpkt/Gv1AQBmjepKfGg7pzq3HcE5tRUsFgdF571+NTCO1KRQqmrNvPTNYbXDEcLzyAzGXkcSHOFQy7Zkc+x0JeEBPtx3Vff2NxRzCej9oLoUio44LkAvpdFo+OvkvgB8sCWbI4XlKkckhAcx152v95MEx2tIgiMcprzaxP/75hAAv7+mJ8G+hvY3pjNA/GDrZTlN1SrDUsIZ1y8Gs0XhX18dUDscITzH6UwwVVpnUY9s52l14XYkwREO88b3xzhTUUvXyABuS0/ueINSh9Nmj0zsg06r4ZvMQjYcLVI7HCE8g+09JmGIdTZ14RVcsidfffVVUlJS8PX1JT09nc2bNze57VVXXYVGo7noZ/LkyfZtZsyYcdHtEyZMcEVXRBPyS6t5+yfrQpqPTOiDQeeAp5bMaNxm3aMCuW24Nbl8elUmFpnhWIiWSf2NV3J6gvPBBx8wd+5c5s+fz/bt2xk0aBDjx4+nsLCw0e0/+eQT8vLy7D979+5Fp9Nx0003NdhuwoQJDbZbunSps7simvHvNQepNlkY2iWM8ZfEOKZR25tN4T6oqXBMm53A78f0JNCoZ8+pUj7fnat2OEK4P5ngzys5PcF58cUXmT17NjNnzqRfv368/vrr+Pv7s2DBgka3Dw8PJzY21v6TkZGBv7//RQmO0WhssF1YWJizuyKakJlXxkfbrbPo/mVy35bXmmqt4DgITgTFArk7HNNmJxAZaOT++gLv51YfpNpkVjkiIdzYubNwxlo7SIIMEfcmemc2Xltby7Zt25g3b579Oq1Wy5gxY9iwYUOr2njnnXeYNm0aAQENZ8Jdt24d0dHRhIWFMXr0aP7xj38QERHRaBs1NTXU1NTY/y4rKwPAZDJhMpna2q1m2dpzdLvuorH+Pf3lfhQFJl4Sw4C4QIf2XZeQhrYsB3PWJiyJlzqs3aZ4y/67c3giizec4FTJORb+dIx7RqXYb/OWPjZF+uf5XNlHzcnN6AElrCt1PsHggseUfdjxdltDoyiK007S5+bmkpCQwPr16xkxYoT9+j//+c98//33bNq0qdn7b968mfT0dDZt2sTw4cPt1y9btgx/f3+6du3K0aNH+ctf/kJgYCAbNmxAp7t4Mce///3vPPHEExddv2TJEvz9/TvQQ3GgRMNrmTp0GoW/pJqJ9HVs+90Lv6L/qaXkhQxhc7eHHNu4l9tUqGHJUR1+OoXHBpsJ6MCgNiG8Ve+8T+mT/ynZYSPZnnKf2uGIFlRVVXHbbbdRWlpKcHDzM0479QhOR73zzjsMGDCgQXIDMG3aNPvlAQMGMHDgQLp37866deu45pprLmpn3rx5zJ071/53WVkZSUlJjBs3rsV/UFuZTCYyMjIYO3YsBoP3faJc2D+tTs9r/7cBqOCOS7swfVIfhz+eJicS3l1KbF0OkyZOBEed/mqCN+2/8RaF7f+3gQMFFRwydOOv9fvHm/rYGOmf53NlH3VL3wUgfvgUYodOcupj2cg+bD/bGZjWcGqCExkZiU6no6CgoMH1BQUFxMbGNnvfyspKli1bxpNPPtni43Tr1o3IyEiOHDnSaIJjNBoxGo0XXW8wGJz25HJm2+7AYDCwYlc+BwoqCPLV89CY3s7pb+IQ0BrQVBZiqMyDsC6Of4xGeMP+MwB/mdyP6Qs28/7mbGaO6kaXiPOner2hj82R/nk+p/fRYoHcbQDouqSjc/H/U/Zh+9prLacWGfv4+JCWlsbatWvt11ksFtauXdvglFVjli9fTk1NDXfccUeLj5OTk0NRURFxcXEdjlm0zrlaM/9eYy3Mm3N1D8ICfJzzQAY/iB1gvSzz4bTZFb2iuLxnJCazwnNfH1Q7HCHcS/FRqC4BvS/E9Fc7GuFgTh9FNXfuXN566y3effddMjMzuf/++6msrGTmzJkATJ8+vUERss0777zD1KlTLyocrqio4E9/+hMbN27kxIkTrF27lilTptCjRw/Gjx/v7O6Ieos2nCS/rJqEUD/uGpni3AeT+XA65C+T+qLRwJe789iRdVbtcIRwH7YvTfGDrbOnC6/i9ATnlltu4YUXXuDxxx8nNTWVnTt3snr1amJirHOlZGVlkZeX1+A+Bw8e5KeffmLWrFkXtafT6di9ezfXXXcdvXr1YtasWaSlpfHjjz82ehpKOF65Cd748TgAfxrfG1/DxYXdDmVPcOQITnv0jQvm10MSAevkf04cVyCEZ5EVxL2aS4qM58yZw5w5cxq9bd26dRdd17t37ybfhP38/Pj6668dGZ5oo9XZWiprzPRPCOa6QfHOf0Dbm0/+bqirAb0ksm318LjefL47ly0nzpKR2fgkm0J0OjLBn1eTRTdEmxw7Xcn6AutIpr9M6otW69xRTQCEpYB/JJhrIW+38x/PC8WG+HLPqG4APL/mMGaLygEJobbaSijYZ70sCY5XkgRHtMkLGYexoOHq3pGM7B7pmgfVaCCpfqoAOU3Vbr+5shsRAT6cKKpifaELElMh3FnuDuss6cGJEOyCI9HC5STBEa22+XgxGZmFaFD407hern1wWVm8w4J8DTw0pidgPc1YXl2nckRCqEjqb7yeJDiiVRRF4Z+rMgEYEa3QMzrQtQHISCqHmDY8mW6R/lTUaXirvlBciE5JVhD3epLgiFb5Ynceu7JL8PfRMTFJhQKO+MGg0UJpFpTnu/7xvYRBp7UffVuw/iS5JedUjkgIFSiKFBh3ApLgiBbV1Jl57usDANwzKoVgJ83p1yxjEET3s16Wozgdck2fKLoHKdTUWeyTNQrRqZRmQ0UBaA0QN1DtaISTSIIjWvS/DSfJLj5HdJCRWZe5ZqmERkkdjkNoNBqmdDED8MmOHPbntn5tFyG8gu09JHaAdbZ04ZUkwRHNKq0y8cq3RwCYO7YX/j4qrs8qdTgO0yUIJvePRVHgma8y1Q5HCNeS+ptOQRIc0axX1x2h9JyJXjGB3DQ0Sd1gbG9GudvBLCOAOmru2B4YdBp+PHyG7w+dVjscIVxH6m86BUlwRJOyi6tY9PMJAOZN7IvOFZP6NSeiJxhDwFQFhfvVjcULJIf7c9eIFACeWZWJ2SJLOIhOoK4G8nZZL8sQca8mCY5o0vNfH6TWbOGyHhFc1TtK7XBAq4XENOtlqcNxiDmjexDsq+dAfjkfb8tROxwhnC9/j3VWdP9I6yzpwmtJgiMatTunhJW7ctForEdvNBo3mflW6nAcKtTfh9+Otk7+9++Mg1TVyqk/4eUuPD3lLu9rwikkwREXURSFf35pLTy9PjWB/gkhKkd0AVlZ3OGmj+xCYpgfBWU1vCOT/wlvJzMYdxqS4IiLrM0sZNPxYnz0Wh4e31vtcBpKqD9FVXQYqorVjcVLGPU6/lS/n1///iiny2tUjkgIJ5IC405DEhzRQJ3ZYh82fPdlXUkIdbM5IvzDIaKH9fKp7erG4kWuHRjPwMQQKmvN/L+1Mvmf8FLlBVCSBWggYYja0QgnkwRHNLBsSzZHT1cS5m/ggau7qx1O4+Q0lcNptRr+MqkvAEs3Z3OksELliIRwglP1tXvR/ayzowuvJgmOsKuoqeOlb6zf3n9/TU+CfQ0qR9QEmdHYKS7tFsGYvjGYLQr/+uqA2uEI4XhSf9OpSIIj7N78/ihnKmpJifDntnQVl2RoSUL9m9OpbdZF84TDPDqxDzqthm8yC9h0rEjtcIRwrFPbrL9ttXzCq0mCIwAoKKvmrfoRNI9M6IOP3o2fGtH9rIvkVZdAyUm1o/EqPaIDmTbMOmP106syscjkf8JbKMr5Cf7iB6sbi3AJN/4UE6704ppDnDOZSesSxoT+sWqH0zy9D8TUryxue8MSDvPQmF4E+OjYlVPKF3vy1A5HCMc4ewKqS0HnA1F91I5GuIAkOIID+WUs35YNwF8m9XGfSf2aE5dq/Z27U80ovFJUkJH7rrQWmD+3+gA1dWaVIxLCAfJ2Wn/HXGL9kiS8niQ4gmdWHcCiwKQBsaR1CVc7nNaJT7X+liM4TnHP5d2ICTaSc/Yci9fLaUDhBWzvFbYvR8LrSYLTyW06VsT3h05j0Gn483gPOmwbN8j6O2+nFBo7gZ+PjofHWif/e+Xbw1TUyBIOwsPZjvba3juE15MEp5N7+ydrYfFNQ5NIiQxQOZo2iL4EtHqoKoKyU2pH45VuTEukW2QAZdV1fLQ1W+1whGg/RTl/isp29Fd4PUlwOrGTRZV8k1kAWGct9igGX4iyTkwndTjOodNqmHlZCgAL15+QEVXCc5Vmw7mz1tGX0f3Ujka4iCQ4ndii9SdQFLiqdxQ9ogPVDqft4i84TSWc4oYhiQT76jlZVMW3BwrVDkeI9rF9CYruC3qjqqEI15EEp5MqrzaxfGsO4IFHb2xsxYJSaOw0AUY9tw5PBmDBz7LSuPBQ9gJjqb/pTCTB6aQ+3JpDRU0dPaIDubxnpNrhtM+FQ8Wl0Nhppo9MQafVsP5oEZl5ZWqHI0TbSf1NpyQJTidktigsWm/9Nn73ZV09Y96bxsRcAhotVBZCeb7a0XithFA/JlxinfxxoRzFEZ5GUS4YQZWqZiTCxSTB6YS+ySwgu/gcof4Grh+coHY47efjf35GUqnDcaq7R6UAsGJnLmcqatQNRoi2KMuFqjOg0Vm/FIlOQxKcTmhB/dDw24Yn4+ejUzmaDrLPhyN1OM40JDmMQYkh1NZZWLIpS+1whGg923tDVB8w+Kkbi3ApSXA6mb2nStl0vBi9VsOdI9x4xfDWkiUbXEKj0XD3KGsx+v82npTlG4TnkPqbTksSnE5m4c8nAJg0II64EC/4NiNHcFxmYv84YoKNnC6v4cvdsgin8BAygqrTkgSnEyksr+bzXbkA9m/jHi92AKCB8lyokHlanMlHr2X6iBQA3vnpOIqMXBOeQAqMOy1JcDqR9zdmUWu2MCQ5lNSkULXDcQxjIET2tF6WozhOd+vwZIx6Lftyy9hy4qza4QjRvPJ8qMi3jraM7a92NMLFXJLgvPrqq6SkpODr60t6ejqbN29ucttFixah0Wga/Pj6+jbYRlEUHn/8ceLi4vDz82PMmDEcPnzY2d3waNUmM+9vsq4K7TVHb2ykDsdlwgN8uGGIdeSdrVhdCLdl+9IT2Qt8PGitPeEQTk9wPvjgA+bOncv8+fPZvn07gwYNYvz48RQWNn06ITg4mLy8PPvPyZMnG9z+3HPP8fLLL/P666+zadMmAgICGD9+PNXV1c7ujsf6fFcuZypqiQ/xtc9p4jVsxYMyVNwlbDNfr9mfT3ZxlcrRCNEMOT3VqTk9wXnxxReZPXs2M2fOpF+/frz++uv4+/uzYMGCJu+j0WiIjY21/8TExNhvUxSFl156ib/97W9MmTKFgQMHsnjxYnJzc1mxYoWzu+ORFEVhQX1x8V0jU9DrvOzMpBQau1TPmCAu7xmJRYF3159QOxwhmiYFxp2a3pmN19bWsm3bNubNm2e/TqvVMmbMGDZs2NDk/SoqKujSpQsWi4UhQ4bw9NNPc8kl1gmajh8/Tn5+PmPGjLFvHxISQnp6Ohs2bGDatGkXtVdTU0NNzfnJycrKrNPNm0wmTCZTh/t5IVt7jm63IzYdLyYzrww/g5YbB8d1KDZ37B+R/TAAlGZjKs0H/4h2N+WW/XMwR/TxrkuT+PHwGZZtyebBq7oSaHTqW0mbePs+9Pb+geP6qM/dgQaoi+6P4kb/L9mHHW+3NZz6rnTmzBnMZnODIzAAMTExHDhwoNH79O7dmwULFjBw4EBKS0t54YUXGDlyJPv27SMxMZH8/Hx7G79s03bbLz3zzDM88cQTF12/Zs0a/P3929O1FmVkZDil3fZ4+4AW0JIWXsfP3zkmLnfqH8A1xhgCawrYsvJtTgcP6HB77tY/Z+hIHy0KRPvqKKyu4x/vZXBFnPuNqPL2fejt/YOO9dHHVMbE8lwUNHy9K5e6vascGJljyD5su6qq1p8Wd5+vXfVGjBjBiBEj7H+PHDmSvn378sYbb/DUU0+1q8158+Yxd+5c+99lZWUkJSUxbtw4goODOxzzhUwmExkZGYwdOxaDweDQttvjZHEVezf+BMBjt1xOt6iOFdq5W/9sdDWfwv5PSU8yYrlsUrvbcdf+OZKj+lgalc3fP89kS1kgT88chVbrHmuaefs+9Pb+gWP6qDm6FvYCEd0Zd+2Njg2wg2Qftp/tDExrODXBiYyMRKfTUVBQ0OD6goICYmNbV+hqMBgYPHgwR44cAbDfr6CggLi4uAZtpqamNtqG0WjEaDQ22raznlzObLst3t+cg6LA1b2j6B0f6rB23aV/dgmDYf+n6Ar3oHNAXG7XPyfoaB9vHpbMixmHySo+x49HzzKmX0zLd3Ihb9+H3t4/6GAfC/cCoIlLddv/k+zD9rXXWk6tNvXx8SEtLY21a9far7NYLKxdu7bBUZrmmM1m9uzZY09munbtSmxsbIM2y8rK2LRpU6vb7CzKq00s35oDeOHQ8F+SoeIu5++j59b0ZAAWyCrjwt3IEg2dntOH08ydO5e33nqLd999l8zMTO6//34qKyuZOXMmANOnT29QhPzkk0+yZs0ajh07xvbt27njjjs4efIk99xzD2AdYfXQQw/xj3/8g5UrV7Jnzx6mT59OfHw8U6dOdXZ3PMqHW3OoqKmjZ3Qgo3pEqh2Oc8UNtP4uOQnnZAI6V5k+IgWdVsP6o0Vk5rX+0LEQTicjqDo9p9fg3HLLLZw+fZrHH3+c/Px8UlNTWb16tb1IOCsrC632fJ519uxZZs+eTX5+PmFhYaSlpbF+/Xr69etn3+bPf/4zlZWV3HvvvZSUlDBq1ChWr1590YSAnZnZorBovfVb9d2juqLRuEd9hNP4hUFYCpw9YX1j63aVygF1DgmhfkzoH8uXu/NY+PNxnvu1fJgIN1BVDCX1q97HDlQ3FqEalxQZz5kzhzlz5jR627p16xr8/Z///If//Oc/zban0Wh48sknefLJJx0Votf5JrOA7OJzhPkbuH5wgtrhuEbcIElwVHD3ZV35cnceK3bm8ucJfYgMvLjeTQiXsh29CesKfqGqhiLU42Uzvgkb2zT6t6Un42vQqRyNi0gdjiqGJIcyKCmU2joLSzZlqR2OEFJ/IwBJcLzS3lOlbDpejF6r4c5LU9QOx3VkRmNVaDQa7r4sBYDFG05SU2dWNyAh7Es0yCnTzkwSHC9kG9EyeWAcsSGdqC7JdgSn+ChUl6oaSmczaUAcMcFGzlTU8MWuPLXDEZ2dvcA4VdUwhLokwfEyheXVfL4rF4CZl3n50PBfCoiAkCTr5fw96sbSyRh0WqaPSAGsCbaiuN/MxqKTOFcCZ+unLZAjOJ2aJDhe5r2NWZjMCmldwkhNClU7HNezvaFJHY7L3TY8GaNey77cMjYfL1Y7HNFZ5e+2/g5NBv9wdWMRqpIEx4tUm8y8v/EkYB3Z0inZigqlDsflwgJ8uGFIIiAT/wkV2etvUtWMQrgBSXC8yMpduRRV1hIf4sv4S9xr2nyXsb2p2UZRCJeyFRuv2V9AVlHrF8UTwmFkgj9RTxIcL6Eoin1o+F0jU9DrOumutb2pnTkMNRXqxtIJ9YwJ4vKekSgKvLvhhNrhiM5IhoiLep30U9D7bDhWxIH8cvwMOqYNS1Y7HPUERkNQPKBIobFKbOuefbAlm/Jqk8rRiE6lugyKrAszyykqIQmOl1jw0wkAfp2WSIi/d69O2yJ7Hc5ONaPotK7sGUW3qAAqaur4aFuO2uGIzsT2pSY4EQK8fP090SJJcLzAiTOVrD1QAMCM+hqITk0m/FOVVquxT1GwaP0JzBYZMi5cROpvxAUkwfECi9afQFHg6t5RdI8KVDsc9cmSDaq7cUgCwb56ThZV8e2BQrXDEZ2F1N+IC0iC4+HKqk0s35oNnK996PTshcYHoVZG8qjB30fPrenWWjBb8bsQTidHcMQFJMHxcB9uyaay1kzP6EBG9ZBzzgAEx0FgDCgWKNirdjSd1vQRKei0GjYcK2J/bpna4QhvV1sJZw5ZL0uBsUASHI9mtigsWn8CsB690Wg06gbkTqQOR3UJoX5M6B8LwEKZ+E84W/4e65eawFgI6qTzgIkGJMHxYBn7C8g5e44wfwPXD05QOxz3InU4bsE2o/ZnO3M5U1GjcjTCq9m+zEj9jagnCY4Hs02Hf1t6Mr4GncrRuBk5guMWhiSHMigplFqzhfc3ZqkdjvBm9iUapP5GWEmC46H2nipl8/Fi9FoNd16aonY47sf2Le50JpiqVQ2lM9NoNPblG/638SQ1dWZ1AxLey15gnKpqGMJ9SILjoWxHbyYPjCM2xFflaNxQcAL4R4ClDgr3qR1NpzZpQByxwb6cqajhi115aocjvJHpHJw+YL0sR3BEPUlwPFBheTWf78oFYJYMDW+cRiN1OG7CoNMyfWQXwJqYK4pM/CccrGAfKGYIiILgeLWjEW5CEhwP9N7GLExmhaFdwhiYGKp2OO7LvmSD1OGo7dZhyfgatOzLLWPz8WK1wxHeJneH9XdcqvXLjRBIguNxqk1m3t94EpCJ/VpkLzTeqWoYAsICfLhhSCJw/vSqEA4jE/yJRkiC42FW7sqlqLKWhFA/xvWTuR6aZTtFVbAf6mpVDUXAzJEpAKzZX0BWkcwwLRxIlmgQjZAEx4MoimKf9v6ukV3Q62T3NSs0GXxDwWKCwv1qR9Pp9YwJ4opeUSgKvLvhhNrhCG9RVwOFmdbLcgRHXEA+IT3IhmNFHMgvx99Hxy1Dk9UOx/1pNBfU4exUMxJRzzZk/IMt2ZRXm9QNRniHgn3W0ZJ+4RCSpHY0wo1IguNBFvx0AoBfpyUS4m9QNxhPIRP+uZUrekbRPSqAipo6PtqWo3Y4whvYvrzEDZICY9GAJDge4sSZStYeKABgRn0tg2gFGSruVrRaDTPrl29Y+PMJzBYZMi46SJZoEE2QBMdDLFp/AkWB0X2i6RYVqHY4nsN2BKdgH5jllIg7uGFIAiF+BrKKq/gms0DtcISnkyUaRBMkwfEAVbXnD+fPrK9hEK0U3g2MIWCuOT/TqVCVv4+eacOttRLv1U95IES71NWeH0AgSzSIX5AExwN8sTuPipo6ukT4c1n3SLXD8SwaDcQNtF6WOhy3cdtwa5H8T0fOkF0sQ8ZFO53OBHMt+IZAWIra0Qg3IwmOB1i22boK87RhyWi1UkTXZrZD11KH4za6RAQwqkckigIfbs1WOxzhqS6c4E8KjMUvSILj5g7ml7M9qwS9VsOv0xLVDscz2Q5dyxEct2I7TfXh1mzqzBaVoxEeSepvRDMkwXFzS+uP3ozpG0NUkFHlaDyUbXRF/h4w16kaijhvbL8YwgN8KCir4buDp9UOR3gi+xGcVFXDEO5JEhw3Vm0y8+mOU8D5b7uiHcK7g08g1J2DosNqRyPqGfU6bhySAJw/DStEq5nroGCv9bIkOKIRkuC4sdV78yk9ZyIh1I/Le0apHY7n0mohtr7QWOpw3Motw6zFxt8dLCSv9JzK0QiPcuYg1FWDT5B1tKQQvyAJjhuznZ66eWgSOiku7hj7kg1Sh+NOekQHMjwlHIsCy7fKzMaiDS6sv9HKR5m4mEueFa+++iopKSn4+vqSnp7O5s2bm9z2rbfe4vLLLycsLIywsDDGjBlz0fYzZsxAo9E0+JkwYYKzu+FSx05XsOl4MVoN3DxMios7zL5kw05VwxAXs51+/WBLNhaZ2Vi01oUjqIRohNMTnA8++IC5c+cyf/58tm/fzqBBgxg/fjyFhYWNbr9u3TpuvfVWvvvuOzZs2EBSUhLjxo3j1KlTDbabMGECeXl59p+lS5c6uysu9cEW69DZq3pHExfip3I0XsA+kmo3WGTEjjuZNCCOYF89p0rO8eORM2qHIzyF7cuKLNEgmqB39gO8+OKLzJ49m5kzZwLw+uuv8+WXX7JgwQIeffTRi7Z///33G/z99ttv8/HHH7N27VqmT59uv95oNBIbG9uqGGpqaqipqbH/XVZWBoDJZMJkcuz0/bb2OtJubZ2F5dusCc5NQ+IdHmNHOKJ/qghJQW/wR2OqxFRwACJ7NrqZx/avDdytjzpgyqA4/rcpmyUbTzCya2iH2nO3/jmat/cPWtFHixl9/h40gCnqEvCw/4Xsw4632xoaRVGcdky4trYWf39/PvroI6ZOnWq//q677qKkpITPPvusxTbKy8uJjo5m+fLl/OpXvwKsp6hWrFiBj48PYWFhjB49mn/84x9EREQ02sbf//53nnjiiYuuX7JkCf7+/u3rnBPtLNKw8JCOYIPC39PM6KT8xiFGHXqKiMrDbOtyHznhI9UOR1zgVCU8t1uPVqPwxBAzwT5qRyTcWdC5U4w+MI86rZEvB74BGqnB6Syqqqq47bbbKC0tJTg4uNltnXoE58yZM5jNZmJiYhpcHxMTw4EDrVsX6JFHHiE+Pp4xY8bYr5swYQI33HADXbt25ejRo/zlL39h4sSJbNiwAZ1Od1Eb8+bNY+7cufa/y8rK7Ke+WvoHtZXJZCIjI4OxY8diMBja1cbyd7cBRdw+ohvXjm38SINaHNE/tWh1P8DWw6TGahg4ZlKj23hy/1rLXfu4ungju3PKKIvoy7TLu7a7HXftn6N4e/+g5T5q9nwIB0Abn8qkyb9SIcKOkX3YfrYzMK3h9FNUHfGvf/2LZcuWsW7dOnx9fe3XT5s2zX55wIABDBw4kO7du7Nu3Tquueaai9oxGo0YjRdPkmcwGJz25Gpv29nFVfx8tAiAW9NT3PbJ78z/ndMkDIGtoMvfg66F2D2yf23kbn28bXgXdufs4aPtuTxwdU80HZx6393652je3j9opo+F1vlvtAmD0Xrw/6BT78MOtNdaTj2uFxkZiU6no6CgoMH1BQUFLdbPvPDCC/zrX/9izZo1DBw4sNltu3XrRmRkJEeOHOlwzGr7cGs2igKjekSSHOF+p888mm20Rb4UGrujawfFE+Cj4/iZSjYeK1Y7HOHOZIkG0QpOTXB8fHxIS0tj7dq19ussFgtr165lxIgRTd7vueee46mnnmL16tUMHTq0xcfJycmhqKiIuLg4h8Stljqzxb7woMxc7ARRfUDvCzVlcPa42tGIXwgw6rkutX5m4y0ys7FogsVi/ZICMoOxaJbTK7Pmzp3LW2+9xbvvvktmZib3338/lZWV9lFV06dPZ968efbtn332WR577DEWLFhASkoK+fn55OfnU1FRAUBFRQV/+tOf2LhxIydOnGDt2rVMmTKFHj16MH78eGd3x6nWHTxNQVkN4QE+jO0X0/IdRNvo9BBzifWyzIfjlm6tT+y/2ptPSVWtytEIt1R8FGorQO8Hkb3Ujka4MacnOLfccgsvvPACjz/+OKmpqezcuZPVq1fbC4+zsrLIy8uzb//aa69RW1vLr3/9a+Li4uw/L7zwAgA6nY7du3dz3XXX0atXL2bNmkVaWho//vhjo3U2nsT2rfXGIQkY9RcXSwsHsH3jkyUb3NKAhBD6xQVTW2fhk+2nWr6D6HxsE/zF9rd+aRGiCS55dsyZM4c5c+Y0etu6desa/H3ixIlm2/Lz8+Prr792UGTuI7+0mm8PWCc/tK3PI5zAPqOxLNngjjQaDbcOT+Kxz/axbEsWMy9L6XCxsfAyuTusv6X+RrRAJg9wEx9uzcaiwPCu4fSIDlQ7HO914ZpUzpsCSnTAlMEJ+Bq0HCqoYHvWWbXDEe7GvkRDqqphCPcnCY4bsFgU+9IMt0pxsXNF9QWdD1SXQMlJtaMRjQj2NTB5QDwASzdnqxyNcCuKYl1uBeQIjmiRJDhu4McjZzhVco5gXz0T+3v2SDC3p/eB6H7Wy1KH47Zsif4Xu3Mpq/be6exFGxUfg5pS0Bkhuq/a0Qg3JwmOG1i22VpcfMOQRHwNUlzsdBeephJuKa1LGD2iA6k2WfhsZ67a4Qh3YXvNxlwCOu+eIE90nCQ4KjtdXkPGfutEiDL3jYvYC413qhqGaJpGo2HaMOvrwfYFQAj7a1ZOT4lWkARHZR9vz6HOopCaFEqfWMeuiyWaYCtOlEJjt3bDkER8dFr25ZaxJ6dU7XCEO7AdwbEdhRWiGZLgqEhRpLhYFdH9QKuHqiIozVE7GtGE8AAfxve3LumyVGY2FooiSzSINpEER0UbjxVz/EwlAT46fjUwXu1wOg+D7/kCRanDcWu31p+mWrkzl8qaOpWjEaoqybKOftQazg8UEKIZkuCoyDZz8XWpCQQYZUZOl5I6HI9wabcIukT4U1FTx5e781q+g/BettdqdF/Qe/as9cI1JMFRSUlVLV/tzQfk9JQqZMkGj6DVaphWP7O3nKbq5KT+RrSRJDgq+WT7KWrrLPSLC2ZAQoja4XQ+9kLjnVJo7OZ+nZaIXqthR1YJB/PL1Q5HqEXqb0QbSYKjAkVR7Kenbh2eJGvtqCG2P2h0UHkayuXUhzuLCjIypq91cd6lMmS8c1KUC5ZoGKxuLMJjSIKjgu1ZJRwqqMDXoGXK4AS1w+mcDH4Q1dt6WQqN3Z5tjqhPd5yi2mRWORrhcmWnoOqM9UtJjBQYi9aRBEcFtonLJg+IJ9hXZuNUjdTheIzLe0aREOpH6TkTq+tr10QnYvsSEt3X+uVEiFaQBMfFyqpNfFE/GkSKi1VmH0klR3DcnU6r4eah1teLnKbqhKT+RrSDJDgu9tnOXM6ZzPSIDiStS5ja4XRu9jWpdqoZhWilm4clotXApuPFHD1doXY4wpXsSzSkqhmF8DCS4LiY7fTUtGFSXKy6mP6AxlpkXF6gdjSiBXEhflzVOxrAPgO46CTsBcZyBEe0niQ4LrQnp5R9uWX46LTcOCRR7XCEMRAie1kvy2kqj2BbgPPjbTnU1llUjka4RFkeVBSARmsd/ShEK0mC40K2icom9I8lLMBH5WgEIHU4HmZ0n2iig4wUVdaSsV+OunUKttdmZC/wCVA3FuFRJMFxkcqaOlbuzAXOD3kVbkDqcDyKXqflpqHWo5/LZGbjzkHqb0Q7SYLjIl/uzqOipo6UCH9GdItQOxxhY5/RWI7geIpbhlqXbvjx8Bmyi6tUjkY4nSzRINpJEhwXsZ2eumVYshQXu5PYAdbfpdlQWaRuLKJVkiP8GdUjEpBi405BhoiLdpIExwUO5pezI6sEvVbDr9OkuNit+AZDRA/rZTlN5TFuHW49irN8WzZ1Zik29loVhVCeC2jOfxkRopUkwXEB28RkY/rGEBVkVDkacRF7ofFOVcMQrTe2XwwRAT4UlNXw3cHTaocjnESTv9t6IaIHGIPUDUZ4HElwnKzaZObTHacAKS52W7Jkg8fx0Wu5sf5o6DKZ2dhr2RMcqb8R7SAJjpOt3ptP6TkTCaF+XN4zSu1wRGNkqLhHuqV+TpzvDhaSV3pO5WiEM2jyZYI/0X6S4DiZ7fTUzUOT0GmluNgt2d48S07CubPqxiJarXtUIMO7hmNRYPnWHLXDEU6gsc9gnKpqHMIzSYLjRMdOV7DpeDFajXUdHeGm/EIhLAUATf4eVUMRbWNbsPaDLdlYLIrK0QhH8qkrR1NWn7jGDVQ3GOGRJMFxItsQ1qt6RxMX4qdyNKJZ9d8QNfk7VQ1DtM3E/nEE++o5VXKOH4+cUTsc4UAhVSesF8K7gW+IqrEIzyQJjpPU1ln4aJv124dt/RzhxupPU9mLGoVH8DXouGGIFBt7o1BbgiP1N6KdJMFxkrUHCimqrCU6yMjoPtFqhyNaUj9KQyOFxh7HNjoxY38Bp8trVI5GOErIuRPWC1J/I9pJEhwn+WCrdWj4TUMT0evk3+z2bKeozh5Hb5bp/z1Jn9hgUpNCqbMo9qOmwvOFVp20XpAjOKKd5JPXCYqq4eej1mn/bevmCDfnHw4h1n0VYntjFR7jfLFxFooixcYe71wJAbWF1suS4Ih2kgTHCTYWWv+tl/eMJDnCX+VoRKvVj9QItR0aFx7jVwPjCfDRcaKoik3HZai/p9MUWGvhlJBk65cPIdpBEhwHqzNb2FRone9m2jA5euNR6utw7KM3hMcIMOq5LjUBgA9kThyPZ6uFU+TojegAlyQ4r776KikpKfj6+pKens7mzZub3X758uX06dMHX19fBgwYwKpVqxrcrigKjz/+OHFxcfj5+TFmzBgOHz7szC602veHzlBq0hAeYGBsvxi1wxFtETcYuGD0hvAottNUX+8voNKkcjCiQ2yjGZVYSXBE+zk9wfnggw+YO3cu8+fPZ/v27QwaNIjx48dTWFjY6Pbr16/n1ltvZdasWezYsYOpU6cydepU9u7da9/mueee4+WXX+b1119n06ZNBAQEMH78eKqrq53dnRZ9UF/keH1qPD56OUDmUeq/LQbW5ENNucrBiLYakBDCJfHBmMwKW87IrOGezLZEgxIrE/yJ9tM7+wFefPFFZs+ezcyZMwF4/fXX+fLLL1mwYAGPPvroRdv/v//3/5gwYQJ/+tOfAHjqqafIyMjgv//9L6+//jqKovDSSy/xt7/9jSlTpgCwePFiYmJiWLFiBdOmTbuozZqaGmpqzg8fLSsrA8BkMmEyOe6rXl5pNd8fsk42dsOgWIe27S5sffLGvmEMRRcUj7Y8F3PuLug6Su2InMKb9+FNaQnsyy1jQ4GW2tpatcNxCm/efwDUlGMoPgZAbWRfFC/sp7fvw8y8cv773RF6aDSMdXAf2/I/c2qCU1tby7Zt25g3b579Oq1Wy5gxY9iwYUOj99mwYQNz585tcN348eNZsWIFAMePHyc/P58xY8bYbw8JCSE9PZ0NGzY0muA888wzPPHEExddv2bNGvz9HVcE/HWOBouio3uQwpEdP3Nkh8OadjsZGRlqh+AUw7UxxJHL4e+XcyyzTO1wnMob96FvHfhodeSf0/DWp9/QNUjtiJzHG/cfQET5AUYBVYZwMn7eAXjvG6m37sPlx7T8VKDldISG3g7uY1VV66fxcGqCc+bMGcxmMzExDWtRYmJiOHDgQKP3yc/Pb3T7/Px8++2265ra5pfmzZvXIGkqKysjKSmJcePGERwc3LZONaNfUSWxm7Kg6ARjx47FYDA4rG13YTKZyMjI8Nr+Kev2wM876BtaQ59Jk9QOxym8fR9urN3Nil35nNAl8uAk7zvF4e37T7vpBByBUv+uXttHb96H52rN/O3574E6RkQrDu+j7QxMazj9FJU7MBqNGI3Gi643GAwO/cf3jA1l3qQAVq067vC23Y239q8uIRUAXeFeNF7Yvwt56z68dXgyK3bls3p/IU+aIdjX+/oI3rv/KLDWW5b4dyHSW/tYzxv34We7CyivriMxzI+eIeUO72Nb2nJqFWxkZCQ6nY6CgoIG1xcUFBAbG9vofWJjY5vd3va7LW0K0Vr2URtnDkFtpbrBiHYZnBRCrJ9CtcnCZztz1Q5HtFXeTgBK/LqqG4doF9uacDenJaBVudbfqQmOj48PaWlprF271n6dxWJh7dq1jBgxotH7jBgxosH2YD1Padu+a9euxMbGNtimrKyMTZs2NdmmEK0WFEu1PgSNYoGCfWpHI9pBo9EwIsYCyAKcHqemAs5Yp/wo9e+icjCirQ4XlLP15Fl0Wg03DklQOxznDxOfO3cub731Fu+++y6ZmZncf//9VFZW2kdVTZ8+vUER8u9//3tWr17Nv//9bw4cOMDf//53tm7dypw5cwDrm9dDDz3EP/7xD1auXMmePXuYPn068fHxTJ061dndEZ1AiX+K9ULuTjXDEB0wLFLBoNOwL7eMPTmlaocjWit/D6CgBMZSYwhVOxrRRsu2ZAMwuk800UEXl4W4mtNrcG655RZOnz7N448/Tn5+PqmpqaxevdpeJJyVlYVWez7PGjlyJEuWLOFvf/sbf/nLX+jZsycrVqygf//+9m3+/Oc/U1lZyb333ktJSQmjRo1i9erV+Pr6Ors7ohMo9U8htmwXyMriHivAAOP7xfDFnnyWbsliQOIAtUMSrZEn8994qpo6M59st84DZ5t0U20uKTKeM2eO/QjML61bt+6i62666SZuuummJtvTaDQ8+eSTPPnkk44KUQi7Er8U64X6WgDhmW4ZmsgXe/L5bMcp/jqpLwHGTjGmwrPVv+aUuEFQoW4oom2+3lfA2SoTcSG+XNkrGou5Tu2QZC0qIX7JfoqqMBNM6s+OLdonvWsYKRH+VNaa+WK3FBt7BDmC47Fs9W43DU1Cp3Z1cT1JcIT4hWpDOIp/JChmKTT2YBqNhlvqF7xdujlb5WhEi2qr4LR1fjQlNlXdWESbnDhTyfqjRWg0cPPQRLXDsZMER4hf0mjOf4OU01Qe7ddpiei1GnZml3Ag37tnpvZ4BftAsUBAFATJlB+exFZcfEXPKBLDHLc6QEdJgiNEI+zz4UiC49GigoyM6Wsd0LBMjuK4N9trLS4VNO5xikO0zGS28NE29youtpEER4hGKHG2BEdGUnm6afVvup9sz6HaZFY5GtEkW4ITn6pmFKKN1mYWcKaihshAI9f0jWn5Di4kCY4QjbCfoirYD3U1zW8s3NrlPaNICPWjrLqOr/bmqR2OaEpu/ZcJ25cL4RFs9W2/TkvEoHOvlMK9ohHCXYQkgV8YWExQuF/taEQH6LQabhlmPYojxcZuylQNpzOtl+NSVQ1FtF7O2Sp+OHwagGnD3Ov0FEiCI0TjNJrz3yTlNJXHu2loIloNbD5ezNHTMsGK2yncB5Y68AuHEPcZhSOa9+HWHBQFRnaPICUyQO1wLiIJjhBNsX2TlCUbPF5ciB9X944G4IMtchTH7dheY/GpUmDsIcwWheVbra+lacOTVY6mcZLgCNEUOYLjVWxvwh9vy6G2zqJyNKKBPKm/8TTfHyokr7SaMH8D4y9xr+JiG0lwhGiKbTRHwT4wm1QNRXTc1b2jiAk2UlRZS8b+ArXDERe6cIi48Ai2erYbhiRi1OtUjqZxkuAI0ZSwrmAMAXONfYZV4bn0Oi03pVkLIZdtyVI5GmFXV2sdrQhyBMdDFJZV8+2BQsD95r65kCQ4QjRFo4G4+uHiUofjFWyjqX48fIbs4iqVoxGAdZSixQS+IRCWonY0ohWWb8vBbFEY2iWMHtFBaofTJElwhGiO1OF4laRwfy7vGQlIsbHbuLD+RgqM3Z7FotiPgLprcbGNJDhCNCd+sPW3LNngNabVL8D54dZs6sxSbKw6qb/xKOuPFpFdfI4gXz2TB8SpHU6zJMERojm2Izj5e8Fcp24swiHG9oshIsCHwvIaex2BUJGMoPIoSzdbj95MTU3Az8c9i4ttJMERojnh3cEnEOrOwZlDakcjHMBHr+XGNOtkcsvkNJW6zCbrlwc4f7RUuK2iihrW7M8Hzq/x5s4kwRGiOVot2Nalkjocr2ErNl53sJC80nMqR9OJnT5oHaXoE2QdtSjc2sfbczCZFQYmhnBJfIja4bRIEhwhWmKbD0fqcLxG96hAhncNx6LAh1ty1A6n87LX3wyyfpkQbktRFPsRT1sdm7uTZ5QQLbEVP8oRHK9im7/jw63ZmC2KytF0UrbXlO1LhHBbm48Xc+x0Jf4+Oq5LjVc7nFaRBEeIltiHiu8Gi1ndWITDTOwfR7CvnlMl5/ixfkVk4WK2+aWkwNjt2Y7eXDswnkCjXuVoWkcSHCFaEtkTDP5gqoSiI2pHIxzE16DjhiH1xcabpdjY5cx1kL/HelmGiLu10ioTq/bkAZ5RXGwjCY4QLdHqIHaA9bKcpvIqtjfrbzILOF1eo3I0nUzRYevoREMARHRXOxrRjE935FBTZ6FPbBCpSaFqh9NqkuAI0Rq2b5iyZINX6RMbzODkUOosCh9tk2Jjl7Kfnhpo/RIh3FLD4uIkNB4027QkOEK0hizZ4LVurR8R8sGWLBRFio1dRib48wg7s0s4kF+OUa/l+sGJaofTJpLgCNEa9qHiu8Ai0/t7k18NiiPQqOdEURUbjhWpHU7nIUs0eARbfdrkAXGE+BtUjqZtJMERojUie4PeF2rL4exxtaMRDuTvo7cPe5ViYxexWKyjEkGO4Lixipo6Pt+dC7j/wpqNkQRHiNbQ6SGmv/Vy7g51YxEOZztNtXpvPmcra1WOphMoOmIdlaj3g8heakcjmrByZy5VtWa6RwUwLCVM7XDaTBIcIVpL6nC81oDEEC6JD6bWbOHj7VJs7HS211Bsf+uXB+GWlm2xLqw5bViyRxUX20iCI0RryZINXs12CH7ZlmwpNnY2qb9xe3tPlbI7pxSDTsMNQxLUDqddJMERorUuPIIjH4BeZ0pqPH4GHUcKK9h28qza4Xg3GUHl9mxHb8ZdEktEoFHlaNpHEhwhWiuqL+h8oLoUzp5QOxrhYMG+BiYPjANgqRQbO4/FImtQubmq2jo+22EtLr7VQxbWbIwkOEK0lt4HovtZL0sdjleyLcD55Z5cSs+ZVI7GS509DjVloDNCVB+1oxGN+HJ3HuU1dSSF+zGye4Ta4bSbJDhCtIXU4Xi1Iclh9IwOpNpkYeXOU2qH451sr52YS0DnWfOqdBbnZy5ORqv1vOJiG0lwhGgLW1GkHMHxShqNxl5svHSzFBs7hZyecmuHCsrZdvIsOq2Gm9I8a+biX3JqglNcXMztt99OcHAwoaGhzJo1i4qKima3/+1vf0vv3r3x8/MjOTmZ3/3ud5SWljbYTqPRXPSzbNkyZ3ZFCCtbUWTuTik09lI3DE7AR6dlf14Ze06VtnwH0Tb2NaikwNgd2Sa7HN0nmuhgX5Wj6RinJji33347+/btIyMjgy+++IIffviBe++9t8ntc3Nzyc3N5YUXXmDv3r0sWrSI1atXM2vWrIu2XbhwIXl5efafqVOnOrEnQtSLuQS0ejhXDKVSiOqNwgJ8mNA/FpBiY4dTlAtGUKWqGoq4WLXJzCc7rPNA2erRPJnTZljKzMxk9erVbNmyhaFDhwLwyiuvMGnSJF544QXi4+Mvuk///v35+OOP7X93796df/7zn9xxxx3U1dWh158PNzQ0lNjYWGeFL0Tj9EaI7gv5e6xv1KGeO8JANG3a8CRW7spl5c5T/G1yXwKMMhmdQ5SchOoS0BqsryPhVr7el09JlYm4EF+u7BWtdjgd5rRX7YYNGwgNDbUnNwBjxoxBq9WyadMmrr/++la1U1paSnBwcIPkBuDBBx/knnvuoVu3btx3333MnDmzyZkWa2pqqKmpsf9dVlYGgMlkwmRy7EgJW3uObtddSP9AFzMQbf4ezDnbsfSY4KrQHEb2YcuGJgXTJdyfk8VVfLYj261qETx5/2myt6EHlOi+1ClaaKIPntzH1nDX/i3ddBKAGwfHYzHXYTG3vy1n9bEt7TktwcnPzyc6umEGqNfrCQ8PJz8/v1VtnDlzhqeeeuqi01pPPvkko0ePxt/fnzVr1vDAAw9QUVHB7373u0bbeeaZZ3jiiScuun7NmjX4+/u3skdtk5GR4ZR23UVn7l9KkY5BwJk937CxKtVlMTlaZ96HrTEwUMPJYh1vfLOPgILdDorKcTxx//XN/ZRewElTGLtWrWpxe0/sY1u4U/9On4ONx/VoUIgoO8SqVYcc0q6j+1hVVdXqbduc4Dz66KM8++yzzW6TmZnZ1mYvUlZWxuTJk+nXrx9///vfG9z22GOP2S8PHjyYyspKnn/++SYTnHnz5jF37twGbSclJTFu3DiCg4M7HOuFTCYTGRkZjB07FoPB+4ZASv9AcyoGFi0m2pzHpIkTwcPWaJF92DrDK2r46vkfOFkB3YdcTu/YIAdG2X6evP90SxcCkDTsVySkTWpyO0/uY2u4Y/+eX3MIOMEVPaO44/ohHW7PWX20nYFpjTYnOA8//DAzZsxodptu3boRGxtLYWFhg+vr6uooLi5usXamvLycCRMmEBQUxKefftriPyc9PZ2nnnqKmpoajMaLp5Q2Go2NXm8wGJz25HJm2+6gU/cvYRBodGgqT2OoPgPBF9eTeYJOvQ9bIS7MwNh+MXy1N5+PduTx9+vCHRhdx3nc/lMUyLceCdMlpqFrRewe18c2cpf+mcwWPtmRB8Ct6V0cGpOj+9iWttqc4ERFRREVFdXidiNGjKCkpIRt27aRlpYGwLfffovFYiE9Pb3J+5WVlTF+/HiMRiMrV67E17flYWo7d+4kLCys0SRGCIcz+FlnYC3cZx3y6qEJjmjZtOHJfLU3n0+25/DoxD74GnRqh+S5SnOgqgg0OutoROE21mYWcKaihshAI9f09fziYhunDRPv27cvEyZMYPbs2WzevJmff/6ZOXPmMG3aNPsIqlOnTtGnTx82b94MWJObcePGUVlZyTvvvENZWRn5+fnk5+djNlurnT7//HPefvtt9u7dy5EjR3jttdd4+umn+e1vf+usrghxsQsX3hRe6/IekSSE+lFWXceqPXlqh+PZbK+V6L5g8Oz5VbzNkvrpEG4amohB5z3z/zq1J++//z59+vThmmuuYdKkSYwaNYo333zTfrvJZOLgwYP2oqHt27ezadMm9uzZQ48ePYiLi7P/ZGdbd4DBYODVV19lxIgRpKam8sYbb/Diiy8yf/58Z3ZFiIZkyYZOQavVcMsw63wgy2ROnI6xvVZk/hu3kl1cxY+HTwMwbZjnz31zIadO7hAeHs6SJUuavD0lJaXBVOhXXXVVi1OjT5gwgQkTPG9orvAycgSn07hpaCIvfXOIzSeKOVJYQY/oQLVD8kz2Cf5kBmN3snxrNooCI7tH0CUiQO1wHMp7jkUJ4UqxAwANlOdBeYHa0Qgnigvx4+re1rqED7ZkqRyNh1KU80s0yBpUbqPObOHDrdaZi21rsHkTSXCEaA+fAIjsZb0sR3G8nu3N/+Ptp6ip68DsZ51VeT5UFoJGCzH91Y5G1Pv+0Gnyy6oJ8zcw/pIYtcNxOElwhGgvqcPpNK7uHUVMsJHiyloy9ssRuzazvUYie4OPcyZXFW1nW2vthiGJGPXeN0JQEhwh2stWLClHcLyeXqflpjRrAeZ7G0+qHI0Hsr1G5PSU2zhVco7vDlrnqvO24mIbSXCEaC9bsaSttkB4tVvTk9FpNWw8Vsz+3NbPpio4/xqRAmO3sXjDCcwWhRHdIugZ4x6zdDuaJDhCtFfcQOvvshyoPKNuLMLpEkL9mNDfOgv7wp+PqxyNh5Eh4m6lqraOpZusBfN3j+qqcjTOIwmOEO1lDIKIHtbLUofTKdx9mfXD4LOduZypqFE5Gg9RXmAdbYimfvShUNvH209RVl1Hlwh/RvfxnpmLf0kSHCE6wvaNVE5TdQpDkkMZlBRKrdnC+xtlyHir2OpvInuCUeYQUpvFotiPQM4YmYJO61mLBbeFJDhCdIRM+NepaDQa7r4sBYD/bTwpQ8ZbQyb4cyvfHz7NsdOVBBn13DTUO4uLbSTBEaIjZKh4pzNpQBwxwUbOVNTwxS5Zn6pFUn/jVhb8ZD16c/OwJAKNTl3MQHWS4AjREbH1hcYlWVBVrG4swiUMOi3TR6QAsODn4y0uL9PpyREct3GooJwfD59Bq7GenvJ2kuAI0RF+oRBWPwpBTlN1GrcNT8bXoGVfbhmbj0ti26TKIiitX6TUNupQqGbhzycAGNcvlqRw759wURIcITpK6nA6nbAAH24YkghYj+KIJthOT4V3A98QVUPp7M5W1vLJduu6U948NPxCkuAI0VFSh9Mpzaw/xL9mfwFZRVXqBuOupP7GbSzZnEVNnYX+CcEMSwlTOxyXkARHiI6SIzidUs+YIK7oFYWiwLsbTqgdjnuS+hu3YDJbWFz/HL37sq5oNN47NPxCkuAI0VG2b6fFx6C6VNVQhGvZhox/sCWb8mqTusG4I9v8ULIGlapW7cmjoKyGqCAjkwfGqR2Oy0iCI0RH+YdDSLL1ct5udWMRLnVFzyi6RwVQUVPHR9ty1A7HvZw7CyX1C5PGSoGxWhRFsQ8Nv/PSLl65anhTJMERwhHibaepdqoahnAtrVbDzPrlGxatty5eKOrZTk+FdrF+CRCq2J5Vwq6cUnz0Wm5LT1Y7HJeSBEcIR5AlGzqtG4YkEOJn4GRRFd8eKFQ7HPchp6fcgm2U39TUeCIDjSpH41qS4AjhCLYERwqNOx1/Hz23Drd+M37np2MqR+NGpMBYdadKzrF6bz6A/UhjZyIJjhCOYHsTLzoCNeXqxiJcbvqILui0GjYeK2ZfrhSaAzJE3A0s3mA9bTqyewR944LVDsflJMERwhECoyA4AVAgf4/a0QgXiw/1Y2L/WOD8bLGdWnWpdVQhSIKjkqraOpZusq54f3cnPHoDkuAI4ThSh9Op2WaHXbkzl9PlNSpHozLbaMKQJAiIUDeWTurj7acoq66jS4Q/o/tEqx2OKiTBEcJRZMK/Tm1IchipSaHUmi28v+mk2uGoS+pvVGWxKCysLy6eOTIFrbZzTOz3S5LgCOEosmRDp2c7ivPexpPU1JlVjkZFUn+jqu8Pn+bY6UqCjHp+PTRJ7XBUIwmOEI5i+7Z65hDUVqobi1DFxP6xxAb7cqails935akdjnrkCI6qbBP73TIsiUCjXuVo1CMJjhCOEhQLgbGgWCB/r9rRCBUYdFqmj+wCWD9kFKUTTvxXUw5nDlsvyxw4LneooJwfD59Bq4G76heE7awkwRHCkaQOp9O7dVgyvgYt+/PK2HS8WO1wXC9/L6BAUBwEds7iVjXZRvGN6xdLUri/usGoTBIcIRxJ6nA6vbAAH24YkgicP1XQqUj9jWrOVtbyyXbrmmi2erDOTBIcIRxJjuAIrCNXADIyC8gqqlI3GFeT+hvVLNmcRU2dhf4JwQxLCVM7HNVJgiOEI9m+tRZmgumcqqEI9fSMCeKKXlEoinURzk5F1qBShclsYfGGE4B1Yj+NpnMODb+QJDhCOFJwPPhHgmKGgv1qRyNUdPdlKQB8uDWb8mqTusG4Sm0VnDlovSxHcFxq1Z48CspqiAoyMnlgnNrhuAVJcIRwJI3mgjqcHaqGItR1Rc8oukcFUFFTx/KtOWqH4xoFe62jCAOirUXGwiUURbHXe915aReMep3KEbkHSXCEcDRZskEAWq3GvoLzovXWRQ+93oWnp+QUictszyphV04pPnott6Unqx2O25AERwhHk0JjUe+GIQmE+BnIKq5ibWaB2uE4nxQYq2JB/bIMU1PjiQw0qhyN+3BqglNcXMztt99OcHAwoaGhzJo1i4qKimbvc9VVV6HRaBr83HfffQ22ycrKYvLkyfj7+xMdHc2f/vQn6urqnNkVIVrPdoqqMBPqOvmii52cv4+eW4dbv1HbPoS8mgwRd7lTJedYvTcfkKHhv+TUBOf2229n3759ZGRk8MUXX/DDDz9w7733tni/2bNnk5eXZ/957rnn7LeZzWYmT55MbW0t69ev591332XRokU8/vjjzuyKEK0XkgR+YWAxQaEUGnd2d43sgk6rYeOxYvbllqodjvOYqq1JPcgRHBdavMF6+vOyHhH0iQ1WOxy34rQEJzMzk9WrV/P222+Tnp7OqFGjeOWVV1i2bBm5ubnN3tff35/Y2Fj7T3Dw+Z22Zs0a9u/fz3vvvUdqaioTJ07kqaee4tVXX6W2ttZZ3RGi9TQaqcMRdnEhfkwaYC24tc0y65UK9llHD/pHQEii2tF0ClW1dSzdlAVYh4aLhpy2CteGDRsIDQ1l6NCh9uvGjBmDVqtl06ZNXH/99U3e9/333+e9994jNjaWa6+9lsceewx/f397uwMGDCAmJsa+/fjx47n//vvZt28fgwcPvqi9mpoaamrOnyooKysDwGQyYTI5dvimrT1Ht+supH+to40ZgO7Yd5hP7cAy6A5HhOYwsg9db3p6Ip/vyuWznad4eEz3DtVJuGP/ALQ529ABltiBmDtYMuCufXQUR/Xvw81ZlFXX0SXcn1Hdwtzq/+WsfdiW9pyW4OTn5xMd3XAdEr1eT3h4OPn5+U3e77bbbqNLly7Ex8eze/duHnnkEQ4ePMgnn3xib/fC5Aaw/91Uu8888wxPPPHERdevWbPGnjg5WkZGhlPadRfSv+bFn7UwDCg7+AM/sMoxQTmY7EPXSgnUcaICnnj/WyYmdXxElbv1b1DWF6QARyoDyVzlmOe8u/XR0TrSP4sC/7dTB2gYGlzO6tVfOS4wB3L0Pqyqav3M4G1OcB599FGeffbZZrfJzMxsa7N2F9boDBgwgLi4OK655hqOHj1K9+7d29XmvHnzmDt3rv3vsrIykpKSGDduXIPTX45gMpnIyMhg7NixGAwGh7btDqR/rXS2L/zffwmtOcWk8WNB5z7/K9mH6lCS8nnow91sOevHC7OuwKhvX4WAu/ZP/84LAHS77Hq69p3UobbctY+O4oj+rTt0msKNOwjy1fPYHaMJMDrteEW7OGsf2s7AtEab/yMPP/wwM2bMaHabbt26ERsbS2FhYYPr6+rqKC4uJjY2ttWPl56eDsCRI0fo3r07sbGxbN68ucE2BQXW4ZdNtWs0GjEaLz4kbDAYnPbicWbb7kD614KonmAMQVNTiuHsEYgb6LjgHET2oWtNHpTAs18fIq+0mtX7T/PrtI7VqbhV/+pq7AXG+qQ0cFBcbtVHJ+hI/xZvzAZg2rAkQgP9HBmWQzl6H7alrTZ/hYiKiqJPnz7N/vj4+DBixAhKSkrYtm2b/b7ffvstFovFnrS0xs6dOwGIi7MW6Y0YMYI9e/Y0SJ4yMjIIDg6mX79+be2OEM6h0ZxPamQ+HAEYdFqmj0gBrKuMK4oXTfxXmGkdNegbCqFd1I7G6x0qKOfHw2fQarA/p8TFnDaKqm/fvkyYMIHZs2ezefNmfv75Z+bMmcO0adOIj48H4NSpU/Tp08d+RObo0aM89dRTbNu2jRMnTrBy5UqmT5/OFVdcwcCB1g+LcePG0a9fP+6880527drF119/zd/+9jcefPDBRo/SCKEa+5INO9WMQriRW4cn4WvQsj+vjE3Hi9UOx3Hs898MkhmMXWBh/ZxK4y+JJSncOXWk3sCp8+C8//779OnTh2uuuYZJkyYxatQo3nzzTfvtJpOJgwcP2ouGfHx8+Oabbxg3bhx9+vTh4Ycf5sYbb+Tzzz+330en0/HFF1+g0+kYMWIEd9xxB9OnT+fJJ590ZleEaDvbUHE5giPqhfr7cOMQ66kp29pBXkFmMHaZ4spaPtl+CpCJ/Vri1Kqk8PBwlixZ0uTtKSkpDQ7TJiUl8f3337fYbpcuXVjloCp9IZzGluDk7wVzHejcqwhQqGPmZSm8vymLjMwCsoqqSI7wgm/gF65BJZxq6eYsauosDEgIYWiXMLXDcWuyFpUQzhLeDXyCoO4cnDmkdjTCTfSIDuLKXlEoinURTo9nNlkn+QNZosHJaussLN5wAoC7R6WgkdOBzZIERwhn0WovKDTeqWoowr3YTi18uDWb8mr3mZytXU4fAHMNGIMhTE6ZONNXe/MoKKshKsjI5AHxaofj9iTBEcKZZMkG0YgrekbSIzqQipo6Ptyao3Y4HWN7bscNsib1wikUReGd+rqt6Zd2waed8yh1JvIfEsKZbEWXUmgsLqDRaJh5WQoAi9Yfx2zx4CHjUmDsEtuzzrI7pxQfvZbb0pPVDscjSIIjhDPZii7zd4PFrGoowr3cMDiRED8D2cXn+CazQO1w2s8+RDxVzSi83oKfTgBwfWoCER1Yy6wzkQRHCGeK6AGGADBVQdERtaMRbsTPR2f/Ju6xQ8bNddZRgiBHcJwo52wVX+3NA2DmqBR1g/EgkuAI4UxaHcQOsF6WOhzxC9NHdEGn1bDpeDF7T5WqHU7bnTlkHSXoE2hN5oVT/G/DSSwKXNYjgj6xjl0/0ZtJgiOEs0kdjmhCXIgfkwZYl6FZ+PMJdYNpD9tzOnaAFBg7SWVNHUs3ZwFw92UySq0t5BkphLPJkg2iGXfXFxt/viuXgrJqdYNpK6m/cbqPtuVQVl1HSoQ/V/eOVjscjyIJjhDOZj+CsxssFnVjEW5ncHIYQ7uEUWu28NI3h9UOp21kBJVTVdbU8cq31tq9WZd3Q6uVif3aQhIcIZwtsjfofaG2HIqPqR2NcEOPTOwDwAdbsjhcUK5yNK1kMVuTdpAlGpzkzR+OcaaihpQIf24ZmqR2OB5HEhwhnE2nh5j+1stymko0YlhKOOMvicGiwDNfHVA7nNYpOgqmStD7QURPtaPxOgVl1bz5g/UL0SMT+sjEfu0g/zEhXEHqcEQLHpnQB71Ww7cHCll/9Iza4bTM9lyOHSALyTrBfzIOcc5kJq1LGBP6x6odjkeSBEcIV5CRVKIF3aIC7fPiPL0qE4u7z24s9TdOc6ignA+3ZgPwl0l9ZFHNdpIERwhXsK9JtUsKjUWTfn9NTwKNevaeKmPlrly1w2le7g7rb0lwHO6ZVZlYFJjYP5a0LuFqh+OxJMERwhWi+1prFWpKZUZj0aSIQCP3X9UdgOe/Pki1yU2X9zDXnU9wEoeqG4uX+fnIGb47eBq9VsMjE/qoHY5HkwRHCFfQGSB+sPVyzhZ1YxFu7e7LuhIb7MupknO8u/6E2uE0rnC/dfkRY7B1lKBwCItF4elVmQDccWkXUiIDVI7Is0mCI4Sr2L7pSoIjmuHno+Phcb0A+O93RzhbWatyRI2wPYcThsgMxg60Yucp9uWWEWTU87trZGRaR8kzUwhXSRxm/Z2zVd04hNu7YUgifeOCKa+u4+Vv3XDyP9tz2PacFh1WbTLzwtcHAXjg6h6EB/ioHJHnkwRHCFexfRgU7oOaCnVjEW5Np9Xwl0nW+ov3Np7kZFGlyhH9gu0IjiQ4DrPw5xPkllYTH+LLzPrlO0THSIIjhKsEx0FwIiiW8wWaQjTh8p5RXNErCpNZ4bnVB9UO57yqYiiqP6qUIAXGjlBcWcv/fWcdfPDH8b3xNehUjsg7SIIjhCtJHY5og3kT+6DRwJd78tiedVbtcKxObbf+Du8GARHqxuIlXl57mPKaOi6JD2ZqaoLa4XgNSXCEcCWpwxFt0DcumF8PSQTg6S8zURQ3mPxPTk851PEzlby38SQAf53UVxbUdCBJcIRwJXuCswXc4cNKuL2Hx/XG16Bl68mzfL0vX+1wJMFxsOdWH6DOonB17yhG9ohUOxyvIgmOEK4UNxC0BqgshJIstaMRHiA2xJfZl3cD4NnVBzGZVZwJ22KBU7YRVFJ/01Hbs0r4am8+Wg3Mm9RX7XC8jiQ4QriSwc+6OCFIHY5otd9c2Z3IQB+On6lk2ZYc9QIpOgLVpaD3hZj+6sXhBRQF/lVfPH7z0CR6xQSpHJH3kQRHCFeTOhzRRoFGPb8fY53875XvjnKuTqVAbEl5/GDr7Nyi3XYVa9iRXYqfQcfcsb3UDscrSYIjhKtdWIcjRCtNG5ZEt6gAzlaZ+CZXpbdue/2NnJ7qiNo6C5+ftO7D2Vd0IzrYV+WIvJMkOEK4mu3DIX831NWoG4vwGAadlkfrF1/8PldDXmm164OQGYwdYumWbM7UaIgM9OE3V3RTOxyvJQmOEK4WlgL+kWCuhbzdakcjPMjYfjEM7RKKSdHwn29cvIRDTYV1Fm6QBKcDSs+ZeHXdMQB+P7oHAUa9yhF5L0lwhHA1jUZOU4l20Wg0PDrBunr3il157Mstdd2D5+6wzsIdnADB8a57XC/z2rqjnK0yEeOn8Osh8n90Jkkdm2E2mzGZTG26j8lkQq/XU11djdlsdlJk6nF1/wwGAzqdF05bnjgUDn0lCY5os0GJIQyJsLC9SMszqw7wv1nD0WhcMDmc1N902KmScyz4+TgA13WxoNfJMQZnkgSnEYqikJ+fT0lJSbvuGxsbS3Z2tmvedFxMjf6FhoYSGxvrXf9PGUklOuBXyRb2lOj46cgZvj90mqt6Rzv/QaX+psNe+PogtXUW0ruGcUnoabXD8XqS4DTCltxER0fj7+/fpg9Wi8VCRUUFgYGBaLXel527sn+KolBVVUVhYSEAcXFxTn08l0oYAmigNAvK8yEoVu2IhAeJ8IU705NZsP4kz6w6wOU9o9A5c4p/RZEZjDto76lSPt1xCoBHx/cma5ckOM4mCc4vmM1me3ITEdH2heQsFgu1tbX4+vp6bYLjyv75+fkBUFhYSHR0tPecrjIGQXQ/a9Fmzlbo+yu1IxIe5v4ru/HR9lMcLCjn42053DwsyXkPVpJlnX1bq4e4Qc57HC+lKApPr8oEYEpqPP0TgsnapXJQnYBTP6GKi4u5/fbbCQ4OJjQ0lFmzZlFRUdHk9idOnECj0TT6s3z5cvt2jd2+bNkyh8Rsq7nx9/d3SHui42z7oq31UG5PVhYXHRDqb+C3o3sC8MKag1TVOnH2P9tzNHaAdTZu0SbrDp5m/dEifPRa/jiut9rhdBpOTXBuv/129u3bR0ZGBl988QU//PAD9957b5PbJyUlkZeX1+DniSeeIDAwkIkTJzbYduHChQ22mzp1qkNj96p6Dw/ntftC6nBEB00f2YXEMD8Ky2t4+8fjznsgqb9ptzqzxX70ZubIFJLC5cuzqzgtwcnMzGT16tW8/fbbpKenM2rUKF555RWWLVtGbm5uo/fR6XTExsY2+Pn000+5+eabCQwMbLCtrfDU9uPrKzNBCg9j+7DI3Q5mtebeF57MqNfx5/rJ/974/iiny500caTU37Tb8m05HC6sINTfwANX91A7nE7FaTU4GzZsIDQ0lKFDzw8pHDNmDFqtlk2bNnH99de32Ma2bdvYuXMnr7766kW3Pfjgg9xzzz1069aN++67j5kzZzb5Tb+mpoaamvMv/LKyMsB6yuOXpz1MJhOKomCxWLBY2r5qr6Io9t/tub+7U6N/FosFRVEwmUxOr8GxPR9ccjostCt6YxCamnJMubvPL8LpZC7towo6W/8m9I1kYEIwu0+V8eKaAzx5XT/HPmBdDfr83WgAU2wquOD/6i37sLKmjhfXWBfUfODKbvjrG37ueHr/muOsPralPaclOPn5+URHNxy6qNfrCQ8PJz8/v1VtvPPOO/Tt25eRI0c2uP7JJ59k9OjR+Pv7s2bNGh544AEqKir43e9+12g7zzzzDE888cRF169Zs+aiWhu9Xk9sbCwVFRXU1ta2Ks7GlJeXt/u+jvarX/2KAQMG8MwzzzisTVf2r7a2lnPnzvHDDz9QV+eaIx0ZGRkueZwRPl2IrtnL/jXvciJytEse08ZVfVRLZ+rflaGw+5SeD7Zk09V0ghgHlsmEVR7hCnMtNfogVq/fD5pMxzXeAk/fh19lazhdoSPCqBBRvI9Vq/Y1uN3T+9caju5jVVVVq7dtc4Lz6KOP8uyzzza7TWZmx18A586dY8mSJTz22GMX3XbhdYMHD6ayspLnn3++yQRn3rx5zJ071/53WVkZSUlJjBs3juDg4AbbVldXk52dTWBgYLtOeymKQnl5OUFBQS6vHZk5cyaLFy++6PoNGzbQt29fgoKC2t22Tqfj448/ZsqUKS7vX3V1NX5+flxxxRVOPxVpMpnIyMhg7NixGAzOXy1Z+/1u+GkvA8Kq6TdpktMfD1zfR1frrP3b//4O1h44zabqOF6/cbDDHk+7+XU4BIauI5k0ebLD2m2ON+zD0+U1zHvpJ8DM41MGMWnA+akgvKF/LXFWH21nYFqjzQnOww8/zIwZM5rdplu3bsTGxtrnL7Gpq6ujuLiY2NiW5/z46KOPqKqqYvr06S1um56ezlNPPUVNTQ1Go/Gi241GY6PXGwyGi/7xZrMZjUaDVqtt1zBo22kbWxuupNFomDBhAgsXLmxwfVRUVLOndmpra/Hx8Wmxfa1Wa09qXNk/2+M2tr+cxWWP1eVS+Am0p7ahdfEbnSv/n2robP2bN6kf6w79wNoDp9mWXcal3do+zUWjcrcDoE0aLs/RNnhl3QGqas2kJoVy3eDERr8QenL/WsvRfWxLW21OcKKiooiKimpxuxEjRlBSUsK2bdtIS0sD4Ntvv8VisZCent7i/d955x2uu+66Vj3Wzp07CQsLazSJcQRFUThnat2yBBaLhXO1ZvS1dR1OAPwMujYfJTEajRclkFdddRWpqam89NJLAKSkpDBr1iwOHz7MihUruOGGG3jzzTeZO3cuH3/8MWfPniUmJob77ruPefPmkZKSAmCvm0pKSuLEiRMd6puol2B9bVB0GKqKwT9c3XiEx+oRHci0YUm8vymLp1dlsuKBy9A6YvI/GUHVZocKyvlgSxYAf53c13tHgro5p9Xg9O3blwkTJjB79mxef/11TCYTc+bMYdq0acTHWxcYO3XqFNdccw2LFy9m+PDh9vseOXKEH374gVWrVl3U7ueff05BQQGXXnopvr6+ZGRk8PTTT/PHP/7RWV3hnMlMv8e/dlr7Tdn/5Hj8fZyzi1544QUef/xx5s+fD8DLL7/MypUr+fDDD0lOTiY7O5vs7GwAtmzZQnR0NAsXLmTcuHFtOgcqWuAfDhE9oOgInNoGPceqHZHwYA+N6cWKHafYnVPK57tzmZKa0LEGy/Ots22jqZ99W7TGv746gEWB8ZfEMCxFvrSoxakzGb///vvMmTOHa665Bq1Wy4033sjLL79sv91kMnHw4MGLPjAXLFhAYmIi48aNu6hNg8HAq6++yh/+8AcURaFHjx68+OKLzJ4925ld8RhffPFFgyH1v5w/yGb06NE8/PDD9r+zsrLo2bMno0aNQqPR0KVLF/tttqNotqH5bTkHKlohcZg1wcnZIgmO6JCoICP3Xdmdf2cc4vmvDzKhfyxGfQdGHtqO3kT3s86+LVq0/ugZvj1QiF6r4ZH6IfxCHU5NcMLDw1myZEmTt6ekpNiHHV/o6aef5umnn270PhMmTGDChAkOi7E1/Aw69j85vlXbWiwWysvKCQoOcsgpqra6+uqree211+x/BwQEcOutt1603YXD9wFmzJjB2LFj6d27NxMmTOBXv/pVowmmcILEobBrqcxoLBzinsu78d6mk+ScPcfi9SeZfUW39jcmK4i3icVyfkmG29KT6RYV2MI9hDPJWlStoNFoWn2qyGKxUOejw99Hr8paVAEBAfTo0fJkUgEBAQ3+HjJkCMePH+err77im2++4eabb2bMmDF89NFHzgpV2NhnNN4GFgt44RpmwnX8fHQ8PLY3f/54N698e5ibhiYS6t/yIIJGSf1Nm6zclcveU2UEGvX8/pqeaofT6ck7qbALDg7mlltu4a233uKDDz7g448/pri4GLCeGjSbW1doLdoo+hLQ+0FNqbXYWIgOujEtkT6xQZRV1/HKt0fa14i5zj6CShKcllWbzDz/tXVSv/uv6k5EoHMGvYjWkwRHAPDiiy+ydOlSDhw4wKFDh1i+fDmxsbGEhoYC1tOJa9euJT8/n5KSElVj9To6/fkCTjlNJRxAp9Xw6ERr/cfiDSfIKmrHwIDC/WCqAmMwRPZycITeZ9H6E5wqOUdssC93X9ZV7XAEkuCIekFBQTz33HMMHTqUYcOGceLECVatWmU/zfbvf/+bjIwMunTpwhVXXKFytF5IVhYXDnZlrygu7xmJyawwf+VezJaL6x2bZXsuJqTJadMWZBVV8ep31iNlfxzfGz8f5y4pI1pHanC8yKJFixq9ft26dQ3+bmwOm9mzZzc7Eu3aa6/l2muvxWKxyCgqZ5CVxYWDaTQa/jKpL1Ne/ZnvDp7mX19l8tfJbVinSupvWqW0ysTMRZspr65jUFIo1w/u4NB84TCSlgvhDhLqj+AU7oca91nHTHi2vnHBPP/rgQC89eNx/rfxZOvvLCuIt6i2zsJ9723j6OlK4kJ8efPONHSOmFxROIQkOEK4g+A4CEkCxQK5O9SORniRKakJ/HGctYZm/md7+e5gYQv3wDqrtq3gXYaIN0pRFP7y6R42HCsiwEfHO3cNIybYuWvlibaRBEcIdyF1OMJJHry6B79OS8SiwJz3t7M/t4XTzKfqR0+Fd5flQ5rw6ndH+GhbDloN/Pf2IfSLD275TsKlJMERwl1IHY5wEo1Gw9PXD2BEtwgqa83MencLBWXVTd9BTk8167Odp3hhzSEAnpjSn6t7R6sckWiMJDhCuAt7grMFGpnhW4iO8NFref2ONLpHBZBXWs3di7ZQWVPX+MYyg3GTtpwo5k/LdwNwz6iu3HlplxbuIdQiCY4Q7iJ2IGgNUHkaStpQDCpEK4X4G1g4YzgRAT7syy3j98t2XDx83GKBUzKCqjEnzlRy7+Kt1JotjOsXw7xJfdUOSTRDEhwh3IXBF+KsI17kNJVwluQIf966ayhGvZZvMgv5x5f7G25QdASqS62za8dcok6Qbqikqpa7F23hbJWJgYkhvDQtVUZMuTlJcIRwJxeephLCSYYkh/HizakALPz5BO+uP3H+RttzL34w6Awuj80d1dSZufd/2zh2ppKEUD/evmtoq9cnFOqRBEcIdyIJjnCRyQPj+POE3gA88fk+vj1QYL1B6m8aUBSFeR/vYfPxYoKMehbMGEZ0kAwH9wSS4AjVzZgxg6lTp6odhnuwfajk7QZTM6NchHCA+6/szi1Dk6zDx5fsYF9uqcxg/Asvrz3CJztOodNqePX2IfSODVI7JNFKkuB4kdOnT3P//feTnJyM0WgkNjaW8ePH8/PPP6sdmmit0C4QEAUWE+TvVjsa4eU0Gg3/uL4/l/WIoKrWzIMLf0Ap3Ge9URIcPt2Rw3++sQ4Hf2pKf67oFaVyRKItJMHxIjfeeCM7duzg3Xff5dChQ6xcuZKrrrqKoqIitUMTraXRyGkq4VIGnZb/uz2NntGBxFYeQKNYsAQnWGfX7sQ2HSvikY/2APCbK7pxW3qyyhGJtpIEpzUUBWorW/9jqmrb9k39tGEulJKSEn788UeeffZZrr76arp06cLw4cOZN28e1113nX2be+65h6ioKIKDgxk9ejS7du1q0M7nn3/OsGHD8PX1JTIykuuvv95+29mzZ7nrrrtISUkhMDCQiRMncvjwYfvtixYtIjQ0lK+//pq+ffsSGBjIhAkTyMvLs29jNpuZO3cuoaGhRERE8Oc//xlF5nxpSGY0Fi4W4mdgwYxhXOZ7HICtdd2pM1tUjko9x05X8Jv3tlFrtjCxfyyPTOijdkiiHaQMvDVMVfB0fKs21QKhjnrcv+SCT0CrNg0MDCQwMJAVK1Zw6aWXYjQaL9rmpptuws/Pj6+++oqQkBDeeOMNrrnmGg4dOkR4eDhffvkl119/PX/9619ZvHgxtbW1rFq1yn7/GTNmcPjwYZYsWUJcXBzz5s1j0qRJ7N+/H4PBOtqiqqqKF154gf/9739otVruuOMO/vjHP/L+++8D8O9//5tFixaxYMEC+vbty7///W8+/fRTRo8e7YB/mJeQGY2FCpLC/ZmeeBqy4OvSJL74Yj9PXHcJGk3nGgpdXGkdDl5SZWJQUigv3pyKVoaDeyRJcLyEXq9n0aJFzJ49m9dff50hQ4Zw5ZVXMm3aNAYOHMhPP/3E5s2bKSwstCc/L7zwAitWrOCjjz7i3nvv5Z///CfTpk3jiSeesLc7aNAgAA4fPszKlSv58ccf6d+/P8HBwbz//vskJSWxYsUKbrrpJgBMJhOvv/463bt3B2DOnDk8+eST9vZeeukl5s2bxw033ADA66+/ztdff+2S/5HHiB8MGi2UZkNZXqc/VSBcRFEIKdoJwE6lJ9s2nCQlIoC7R3VVNy4XqjaZuXfxVk4UVVmHg08fip+PTu2wRDtJgtMaBn/r0ZRWsFgslJWXExwUhFbbwTOABv82bX7jjTcyefJkfvzxRzZu3MhXX33Fc889x9tvv01lZSUVFRVEREQ0uM+5c+c4evQoADt37mT27NmNtp2ZmYleryc9PZ3KykoAIiIi6N27N5mZmfbt/P397ckNQFxcHIWF1tWLS0tLycvLIz093X67Xq9n6NChcprqQsYgiO4HBXutM8oGX6t2RKIzKMmCykLQGpg4Zhzb1pzgqS/3kxTuz9h+MWpH53SKovDIx7vZevIsQb56Fs0cRlTQxUfCheeQBKc1NJpWnyrCYgGD2bp9RxOcdvD19WXs2LGMHTuWxx57jHvuuYf58+fzwAMPEBcXx7p16y66T2hoKAB+fn4dfnzbqSobjUYjyUt7JA61Jjg5W6CvJDjCBWw1X7EDmHV1P46WWFi6OYvfLd3Bh78ZwYDEEHXjc7L/fHOYz3bmotdqeO32NHrGyHBwTydFxl6uX79+VFZWMmTIEPLz89Hr9fTo0aPBT2RkJAADBw5k7dq1jbbTt29f6urq2LRpk/26oqIiDh48SL9+/VoVS0hICHFxcQ3aqKurY9u2bR3ooZeSOhzhahfMf6PRaHhyyiVc3jOScybr6uO5JefUjc+JPt6Ww8trrQMm/nl9f0b1jFQ5IuEIkuB4iaKiIkaPHs17773H7t27OX78OMuXL+e5555jypQpjBkzhhEjRjB16lTWrFnDiRMnWL9+PX/961/ZutX6xjZ//nyWLl3K/PnzyczMZM+ePTz77LMA9OzZkylTpvCb3/yGDRs2sGvXLu644w4SEhKYMmVKq+P8/e9/z7/+9S9WrFjBgQMHeOCBBygpKXHGv8Sz2RKcU9vB3MSKz0I4kn0GY+tzz6DTWie2iwmisLyGuxdtobzapGKAzrHhaBGPfmKdc+qBq7pzyzAZDu4tJMHxEoGBgaSnp/Of//yHK664gv79+/PYY48xe/Zs/vvf/6LRaFi1ahVXXHEFM2fOpFevXkybNo2TJ08SE2M9v37VVVexfPlyVq5cSWpqKqNHj2bz5s32x1i4cCFDhgxh2rRpXHbZZSiKwqpVqy46LdWchx9+mDvvvJO77rqLESNGEBQU1GAouqgX0ROMIVB3DmwTrwnhLHU15yeWvGCJhmBfAwvqa1EO5JczZ8kOrxo+fqSwgt/8bysms8LkgXH8cVxvtUMSDiQ1OF7CaDTyzDPP8MwzzzS5TVBQEC+//DIvv/xyk9vccMMN9hFOvxQWFsa7775LWVkZwcHBFxVRz5gxgxkzZjS4burUqQ1qcPR6PS+99BIvvfRSy53qzLRaSEyDo99av1nHDVI7IuHN8naDuRb8IyEspcFNCaF+vHPXUG5+YwPfHzrN/JX7+MfU/h4/fLyownpUqqy6jiHJofz7pkEyHNzLyBEcIdyV1OEIV7nw9FQjicvAxFD+37TBaDTw/qYs3vnpuIsDdKxqk5nZi7eSVVxFcrg/b00fiq9BhoN7G0lwhHBXsmSDcJVWrCA+/pJY/jqpLwD/XJXJ6r35rojM4SwWhT8u38X2rBKCfa2rg0cEynBwbyQJjhDuKiHN+rvoCFQVqxuL8G6tXEF81qiu3HFpMooCD32wg13ZJc6PzcH+nXGQL3bnoddqeP3ONHpEB6odknASSXCEcFf+4RDRw3r5lAylF05Sng+lWYAGEoY0u6lGo+Hv117Clb2iqDZZmPXuFt7fdJKKGvcf6XfsdAV/X7mPV7+zTmz6zA0DGNldhoN7M0lwmmCxeM9IAU/XqfeFnKYSzmY7ehPdzzqLdgv0Oi3/vW0wfeOCOVNRy18/3cvwf37DvE92syen1MnBtk1NnZmVu3K59c2NjP739yxafwKA347uwU1Dk9QNTjidjKL6BR8fH7RaLbm5uURFReHj49Om0QIWi4Xa2lqqq6s7vlSDG3Jl/xRFoba2ltOnT6PVavHx8XHq47mlxKGwa6kkOMJ5WlF/80tBvgY+/M2lfLAlmyWbszh2upKlm7NZujmbAQkh3Do8metS4wk0qvMRc+JMJUs3Z7F8Ww7FlbWAtXb66t7R3HFpMlf3jlYlLuFakuD8glarpWvXruTl5ZGb27r1py6kKArnzp3Dz8/P44dRNkaN/vn7+5OcnOyVCWOL7EdwtlmXAemM/wPhXK2sv/mlIF8D91zejVmjurLpeDFLNmWxem8+e06VsufTPfzzy/1MGZzAbcOT6Z/g/GUeaussrNmfz5JNWaw/WmS/PibYyC3DkrllWBIJoR1fjkZ4DklwGuHj40NycjJ1dXWYzeY23ddkMvHDDz9wxRVXtGkCPE/h6v7pdDr0er1XJoutEn0J6P2gphSKDkOUTEQmHMhcB7nbrZfbmODYaDQaLu0WwaXdIiiurOXjbTks3ZzFsTOVLNmUxZJNWQxMDOG24clcOyieAAcf1TlZZD169NG2bM5UnD9ac2WvKG4bnszoPtHodfLFoDOSBKcJGo0Gg8HQ5g9xnU5HXV0dvr6+XpngeHv/3I5Oby38PPmz9VSCJDjCkQr3g6kKjMEQ2avDzYUH+DD7im7cc3lXNh4rZsnmLFbvzWN3Tim7c/bwjy8zmZIaz23pyVwS3/6jOrV1Fr7JLGDJpix+OnLGfn10kJFbhiVxy7AkEsP8O9wf4dkkwRHC3SUOPZ/gDL5D7WiEN7HV3ySkOfT0p0ajYUT3CEZ0j6Cooh8f1R/VOVFUxfubsnh/UxaDkkK5fXgyvxoUh79P6z6KsoqqWLoli+VbczhTUVP/WHBFzyhuHZ7MNX2jMcjRGlHPac+Ef/7zn4wcORJ/f39CQ0NbdR9FUXj88ceJi4vDz8+PMWPGcPjw4QbbFBcXc/vttxMcHExoaCizZs2ioqLCCT0Qwk3IjMbCWdpZf9MWEYFGfnNld759+CqW3JPO5IFxGHQadmWX8OePd5P+z7U8tmIvmXlljd7fZLbw1Z487nxnE1c8/x2vrTvKmYoaooKMPHh1d37409W8e/dwJvSPleRGNOC0Izi1tbXcdNNNjBgxgnfeeadV93nuued4+eWXeffdd+natSuPPfYY48ePZ//+/fj6+gJw++23k5eXR0ZGBiaTiZkzZ3LvvfeyZMkSZ3VFCHUl1I9uKdwPNeWtGsorRKv8YgVxZ9JqNYzsEcnIHpGcqaixH9U5WVTF/zae5H8bTzI4OZRbhyczoW8URdXwYsZhPtqRy+nyGns7l/eM5Pb0ZK7pGyMJjWiW0xKcJ554AoBFixa1antFUXjppZf429/+xpQpUwBYvHgxMTExrFixgmnTppGZmcnq1avZsmULQ4da3/RfeeUVJk2axAsvvEB8fHyjbdfU1FBTc/4FUlpqnauhuLgYk8nU3i42ymQyUVVVRVFRkVfWqEj/1OCD3icOTXkudfu/Q0m+rEOtuWcfHUf610rnSjDkHrK26dcNiopauIPjaICb+odyY78Qtpw8y8fbT/HdwTNsO1zFtsO5POajpbrWjMI+ACICDExJjWdqahxJ9bU1ZSVnXRavo3n7cxSc18fy8nKABos4N0lxsoULFyohISEtbnf06FEFUHbs2NHg+iuuuEL53e9+pyiKorzzzjtKaGhog9tNJpOi0+mUTz75pMm258+frwDyIz/yIz/yIz/y4wU/2dnZLeYVblNknJ9vXbgtJiamwfUxMTH22/Lz84mObjhBk16vJzw83L5NY+bNm8fcuXPtf1ssFoqLi4mIiHD48OOysjKSkpLIzs4mODjYoW27A+mf5/P2Pkr/PJ+399Hb+wfO66OiKJSXlzd5xuZCbUpwHn30UZ599tlmt8nMzKRPnz5tadbpjEYjRmPD1WJbW/jcXsHBwV77xAXpnzfw9j5K/zyft/fR2/sHzuljSEhIq7ZrU4Lz8MMPM2PGjGa36datW1uatIuNjQWgoKCAuLg4+/UFBQWkpqbatyksLGxwv7q6OoqLi+33F0IIIYRoU4ITFRVFVFSUUwLp2rUrsbGxrF271p7QlJWVsWnTJu6//34ARowYQUlJCdu2bSMtLQ2Ab7/9FovFQnp6ulPiEkIIIYTncdoYu6ysLHbu3ElWVhZms5mdO3eyc+fOBnPW9OnTh08//RSwTgz10EMP8Y9//IOVK1eyZ88epk+fTnx8PFOnTgWgb9++TJgwgdmzZ7N582Z+/vln5syZw7Rp01p1Ps4VjEYj8+fPv+iUmLeQ/nk+b++j9M/zeXsfvb1/4B591ChKa8Zatd2MGTN49913L7r+u+++46qrrrI+uEbDwoUL7ae9FEVh/vz5vPnmm5SUlDBq1Cj+7//+j169zk8hXlxczJw5c/j888/RarXceOONvPzyywQGBjqjG0IIIYTwQE5LcIQQQggh1CLTQAohhBDC60iCI4QQQgivIwmOEEIIIbyOJDhCCCGE8DqS4LTDq6++SkpKCr6+vqSnp7N58+Zmt1++fDl9+vTB19eXAQMGsGrVKhdF2jbPPPMMw4YNIygoiOjoaKZOncrBgwebvc+iRYvQaDQNfmwrv7ubv//97xfF2tKs256y72xSUlIu6qNGo+HBBx9sdHt3338//PAD1157LfHx8Wg0GlasWNHgdkVRePzxx4mLi8PPz48xY8Zw+PDhFttt62vYmZrro8lk4pFHHmHAgAEEBAQQHx/P9OnTyc3NbbbN9jzXnaWlfThjxoyLYp0wYUKL7brLPmypf429HjUaDc8//3yTbbrT/mvN50J1dTUPPvggERERBAYGcuONN1JQUNBsu+197baFJDht9MEHHzB37lzmz5/P9u3bGTRoEOPHj79ohmWb9evXc+uttzJr1ix27NjB1KlTmTp1Knv37nVx5C37/vvvefDBB9m4cSMZGRmYTCbGjRtHZWVls/cLDg4mLy/P/nPy5EkXRdx2l1xySYNYf/rppya39aR9Z7Nly5YG/cvIyADgpptuavI+7rz/KisrGTRoEK+++mqjtz/33HO8/PLLvP7662zatImAgADGjx9PdXV1k2229TXsbM31saqqiu3bt/PYY4+xfft2PvnkEw4ePMh1113XYrttea47U0v7EGDChAkNYl26dGmzbbrTPmypfxf2Ky8vjwULFqDRaLjxxhubbddd9l9rPhf+8Ic/8Pnnn7N8+XK+//57cnNzueGGG5pttz2v3TZrcTlO0cDw4cOVBx980P632WxW4uPjlWeeeabR7W+++WZl8uTJDa5LT09XfvOb3zg1TkcoLCxUAOX7779vcpvWrhbvDubPn68MGjSo1dt78r6z+f3vf690795dsVgsjd7uSfsPUD799FP73xaLRYmNjVWef/55+3UlJSWK0WhUli5d2mQ7bX0Nu9Iv+9iYzZs3K4By8uTJJrdp63PdVRrr31133aVMmTKlTe246z5szf6bMmWKMnr06Ga3cdf9pygXfy6UlJQoBoNBWb58uX2bzMxMBVA2bNjQaBvtfe22lRzBaYPa2lq2bdvGmDFj7NdptVrGjBnDhg0bGr3Phg0bGmwPMH78+Ca3dyelpaUAhIeHN7tdRUUFXbp0ISkpiSlTprBv3z5XhNcuhw8fJj4+nm7dunH77beTlZXV5LaevO/A+nx97733uPvuu9FoNE1u50n770LHjx8nPz+/wT4KCQkhPT29yX3UntewuyktLUWj0bS4YHBbnutqW7duHdHR0fTu3Zv777+foqKiJrf15H1YUFDAl19+yaxZs1rc1l333y8/F7Zt24bJZGqwP/r06UNycnKT+6M9r932kASnDc6cOYPZbCYmJqbB9TExMeTn5zd6n/z8/DZt7y4sFgsPPfQQl112Gf37929yu969e7NgwQI+++wz3nvvPSwWCyNHjiQnJ8eF0bZOeno6ixYtYvXq1bz22mscP36cyy+/nPLy8ka399R9Z7NixQpKSkqaXSDXk/bfL9n2Q1v2UXtew+6kurqaRx55hFtvvbXZFZrb+lxX04QJE1i8eDFr167l2Wef5fvvv2fixImYzeZGt/fkffjuu+8SFBTU4ukbd91/jX0u5Ofn4+Pjc1HC3dLnom2b1t6nPdq02KboPB588EH27t3b4nnfESNGMGLECPvfI0eOpG/fvrzxxhs89dRTzg6zTSZOnGi/PHDgQNLT0+nSpQsffvhhq75ReZp33nmHiRMnNrtOmyftv87OZDJx8803oygKr732WrPbetJzfdq0afbLAwYMYODAgXTv3p1169ZxzTXXqBiZ4y1YsIDbb7+9xUJ+d91/rf1ccBdyBKcNIiMj0el0F1WHFxQUEBsb2+h9YmNj27S9O5gzZw5ffPEF3333HYmJiW26r8FgYPDgwRw5csRJ0TlOaGgovXr1ajJWT9x3NidPnuSbb77hnnvuadP9PGn/2fZDW/ZRe17D7sCW3Jw8eZKMjIxmj940pqXnujvp1q0bkZGRTcbqqfvwxx9/5ODBg21+TYJ77L+mPhdiY2Opra2lpKSkwfYtfS7atmntfdpDEpw28PHxIS0tjbVr19qvs1gsrF27tsG34AuNGDGiwfYAGRkZTW6vJkVRmDNnDp9++inffvstXbt2bXMbZrOZPXv2EBcX54QIHauiooKjR482Gasn7btfWrhwIdHR0UyePLlN9/Ok/de1a1diY2Mb7KOysjI2bdrU5D5qz2tYbbbk5vDhw3zzzTdERES0uY2WnuvuJCcnh6KioiZj9cR9CNYjqmlpaQwaNKjN91Vz/7X0uZCWlobBYGiwPw4ePEhWVlaT+6M9r932Bi/aYNmyZYrRaFQWLVqk7N+/X7n33nuV0NBQJT8/X1EURbnzzjuVRx991L79zz//rOj1euWFF15QMjMzlfnz5ysGg0HZs2ePWl1o0v3336+EhIQo69atU/Ly8uw/VVVV9m1+2b8nnnhC+frrr5WjR48q27ZtU6ZNm6b4+voq+/btU6MLzXr44YeVdevWKcePH1d+/vlnZcyYMUpkZKRSWFioKIpn77sLmc1mJTk5WXnkkUcuus3T9l95ebmyY8cOZceOHQqgvPjii8qOHTvsI4j+9a9/KaGhocpnn32m7N69W5kyZYrStWtX5dy5c/Y2Ro8erbzyyiv2v1t6Dbtac32sra1VrrvuOiUxMVHZuXNng9dlTU2NvY1f9rGl57q79K+8vFz54x//qGzYsEE5fvy48s033yhDhgxRevbsqVRXVzfZP3fahy09RxVFUUpLSxV/f3/ltddea7QNd95/rflcuO+++5Tk5GTl22+/VbZu3aqMGDFCGTFiRIN2evfurXzyySf2v1vz2u0oSXDa4ZVXXlGSk5MVHx8fZfjw4crGjRvtt1155ZXKXXfd1WD7Dz/8UOnVq5fi4+OjXHLJJcqXX37p4ohbB2j0Z+HChfZtftm/hx56yP6/iImJUSZNmqRs377d9cG3wi233KLExcUpPj4+SkJCgnLLLbcoR44csd/uyfvuQl9//bUCKAcPHrzoNk/bf999912jz0lbHywWi/LYY48pMTExitFoVK655pqL+t2lSxdl/vz5Da5r7jXsas318fjx402+Lr/77jt7G7/sY0vPdVdqrn9VVVXKuHHjlKioKMVgMChdunRRZs+efVGi4s77sKXnqKIoyhtvvKH4+fkpJSUljbbhzvuvNZ8L586dUx544AElLCxM8ff3V66//nolLy/vonYuvE9rXrsdpal/YCGEEEIIryE1OEIIIYTwOpLgCCGEEMLrSIIjhBBCCK8jCY4QQgghvI4kOEIIIYTwOpLgCCGEEMLrSIIjhBBCCK8jCY4QQgghvI4kOEIIIYTwOpLgCCGEEMLrSIIjhBBCCK/z/wFmqo5EQpv8AAAAAABJRU5ErkJggg==", - "text/plain": [ - "
" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], - "source": [ - "p = ex.plot_state_1d(\n", - " jnp.concatenate((u_0_first, u_0_second)), labels=[\"First\", \"Second\"]\n", - ")" - ] - }, - { - "cell_type": "code", - "execution_count": 12, - "metadata": {}, - "outputs": [ - { - "data": { - "text/html": [ - "" - ], - "text/plain": [ - "" - ] - }, - "execution_count": 12, - "metadata": {}, - "output_type": "execute_result" - } - ], - "source": [ - "HTML(\n", - " ex.make_animation_1d(\n", - " jnp.concatenate([trj_first, trj_second], axis=1),\n", - " domain_extent=DOMAIN_EXTENT,\n", - " dt=DT,\n", - " include_init=True,\n", - " ).to_html5_video()\n", - ")" - ] - }, - { - "cell_type": "code", - "execution_count": null, - "metadata": {}, - "outputs": [], - "source": [] - } - ], - "metadata": { - "kernelspec": { - "display_name": "jax_fresh", - "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.10.13" - } - }, - "nbformat": 4, - "nbformat_minor": 2 -}