{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Periodic boundaries\n" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "This tutorial will show how to implement Periodic boundary conditions (where particles that leave the domain on one side enter again on the other side) can be implemented in Parcels\n" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "In case you have a global model on a C-grid, then there is typically not much that needs to be done. The only thing you have to be aware of, is that there is no 'gap' between the longitude values on the right side of the grid and on the left side of the grid (modulo 360). In other words, that\n", "\n", "```python\n", "fieldset.U.lon[:, 0] >= fieldset.U.lon[:, -1]\n", "```\n", "\n", "This is the case for most [curvilinear NEMO model outputs on a C-grid](https://docs.oceanparcels.org/en/latest/examples/tutorial_nemo_curvilinear.html), so that periodic boundary conditions work out of the box, without any special code required.\n", "\n", "For [Hycom models](https://github.com/OceanParcels/parcels/issues/1240), we have found that it can help to add an extra column of longitudes at the end of the longitude array:\n", "\n", "```python\n", "lon_c = np.expand_dims(fieldset.U.grid.lon[:, 0], axis=1)\n", "fieldset.U.grid.lon = np.hstack((fieldset.U.grid.lon, lon_c))\n", "```\n" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "In case of 'simple' models on an A-grid, the idea in Parcels is to do two things:\n", "\n", "1. Extend the fieldset with a small 'halo'\n", "2. Add a periodic boundary kernel to the `.execute`\n" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "We'll start by importing the relevant modules\n" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import math\n", "from datetime import timedelta as delta\n", "\n", "import matplotlib.pyplot as plt\n", "import xarray as xr\n", "\n", "from parcels import (\n", " AdvectionRK4,\n", " FieldSet,\n", " JITParticle,\n", " ParticleSet,\n", " download_example_dataset,\n", ")" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "We import the Peninsula fieldset; note that we need to set `allow_time_extrapolation` because the Peninsula fieldset has only one time snapshot.\n" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "example_dataset_folder = download_example_dataset(\"Peninsula_data\")\n", "fieldset = FieldSet.from_parcels(\n", " f\"{example_dataset_folder}/peninsula\", allow_time_extrapolation=True\n", ")" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "Extending the fieldset with a halo is very simply done using the `add_periodic_halo()` method. Halos can be added either in the zonal direction, the meridional direction, or both, by setting `zonal` and/or `meridional` to `True`.\n", "\n", "But before we apply the halo, we first define two new fieldset constants `halo_east` and `halo_west`. They store the original zonal extend of the grid (so _before_ adding the halo) and will be used later in the `periodicBC` kernel.\n", "\n", "**_Note that some hydrodynamic data, such as the global ORCA grid used in NEMO, already has a halo._** In these cases, **do not** extent the fieldset with the halo but only add the periodic boundary kernel, where you use the explicit values for halo_east and halo_west\n" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "fieldset.add_constant(\"halo_west\", fieldset.U.grid.lon[0])\n", "fieldset.add_constant(\"halo_east\", fieldset.U.grid.lon[-1])\n", "\n", "fieldset.add_periodic_halo(zonal=True)" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "The other item we need is a custom Kernel that can move the particle from one side of the domain to the other.\n" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "def periodicBC(particle, fieldset, time):\n", " if particle.lon < fieldset.halo_west:\n", " particle_dlon += fieldset.halo_east - fieldset.halo_west\n", " elif particle.lon > fieldset.halo_east:\n", " particle_dlon -= fieldset.halo_east - fieldset.halo_west" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "Now define a particle set and execute it as usual\n" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "INFO: Output files are stored in PeriodicParticle.zarr.\n", "100%|██████████| 86400.0/86400.0 [00:00<00:00, 202338.06it/s]\n" ] } ], "source": [ "pset = ParticleSet.from_line(\n", " fieldset, pclass=JITParticle, size=10, start=(20e3, 3e3), finish=(20e3, 45e3)\n", ")\n", "\n", "output_file = pset.ParticleFile(name=\"PeriodicParticle\", outputdt=delta(hours=1))\n", "\n", "pset.execute(\n", " [AdvectionRK4, periodicBC],\n", " runtime=delta(hours=24),\n", " dt=delta(minutes=5),\n", " output_file=output_file,\n", ")" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "And finally plot the particle trajectories\n" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAloAAAGwCAYAAABxbMuTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAADTm0lEQVR4nOydd3gj1dX/P6MuuUjuvW3z9t5hqYGlLiUJkE0BwgskIQFCCSk/kmwSygJJSGgheUMJSYCXUJaaUDaU3WX7sr3bu65ykywXdWl+f8iWLUuuK9uSfT/Po0fSzJ07d6TR6DvnnHuOJMuyjEAgEAgEAoEg6ihGewACgUAgEAgEYxUhtAQCgUAgEAiGCSG0BAKBQCAQCIYJIbQEAoFAIBAIhgkhtAQCgUAgEAiGCSG0BAKBQCAQCIYJIbQEAoFAIBAIhgnVaA9gLOH3+6mpqSEpKQlJkkZ7OAKBQCAQCAaALMu0traSm5uLQhFdG5QQWlGkpqaGgoKC0R6GQCAQCASCIVBZWUl+fn5U+xRCK4okJSUBgS8qOTl5lEcjEAgEAoFgILS0tFBQUBD8H48mQmhFkU53YXJyshBaAoFAIBDEGcMR9iOC4QUCgUAgEAiGCSG0BAKBQCAQCIYJIbQEAoFAIBAIhgkhtAQCgUAgEAiGCSG0BAKBQCAQCIYJIbQEAoFAIBAIhgkhtAQCgUAgEAiGCSG0BAKBQCAQCIYJIbQEAoFAIBAIhgkhtAQCgUAgEAiGCSG0BAKBQCAQCIYJIbQEAsG4o8bpZoO1lRqne8z2JRAIYgNRVFogEMQNNU43ZQ4XE/RacnWaIfXxz5om7jpciZ/AneYDU/K5OjsVPyAjI8sgA35ZDjwDsgx+ZAD83da/XmflvrLaYF9rJuWyOjcNjaRAJQ2uQG3PcT1SWsDq3LQhHSNE57MSCASnjiTLsjzagxgrtLS0YDQasdlsJCcnj/ZwBIKYYCh/+LIs4/DLtHh92Lw+Wrw+1tVb+WtVIzIgARelG5mWqMfp9+P0+3H5ZRw+f8d7GWf3134/Dp+fdp+fZq9vWI+3OxpJQq2Q0Cok1B2vNZKi41lC07FcRmaLzR6yrQR8pyCDHK2aRKWSBJWCJKWSRKWCRFXos0YR6pyItmgTCMY6w/n/LYRWFBFCSzCWOFWLiF+W+WtVA784VoOfgHBYnZPKjER9iIDq/tzq9Qffe2L00iQBCqnTyhUbaBUSCUoFiUolOoXEEbsrZL0E/LgkmxKDjlS1kjS1ihS1ihS1Eq2i7wgSYRkTjAeE0IoThNASjBUiWUS+lpNKu89Po8dLk9tLo8dLo7vj4fF0PIeuO1UhogCMKiUaSaLO4w1bf35aMsV6LTqFhE6pQKdQBF/rO18rFMF1rV4v1+wuCxmXAvhk8VRytWokSUKiS0xJSB3PHY9ursAap5uFnx8I6UsJfLJkKulqFR5Zxu2Xg89uvx+3LOPptqzzuc7l5t5jNXS/GEvAFVkm/DK0+fy0eX20+/y0+ny0ef20+3w4/Kd++U5UKkhRq0hVK0lVqzoegdfH7C5eq7Mid3xOv5mcx/V56YNyiQoE8YAQWnGCEFqC0WYo1gevX6bB48Hs8mJ2uTnc7mRtuZmeFwatBK4oXC2WGhOYaNCSpFJiVClJ7uM5QalAkqReRc22ZdMHbWX5Z00Tdx+uxNfRx8On4FYb7b68fpk2n482n5/WDiF20uHi+wcrwkTbmSlJOPx+LB4vFo8Pq2doQlgnSWRp1WRp1WRqVGRpIr1Wk6pWoughTIVlTBCrCKEVJwihJRhNelqhHi7N56IME2aXJ/BweyK+bhik5UmvkEjXqElXq0jXBB5palXI+3S1Cp8MF+44EhVx1Hl80RI1NU435Q4XJVH404/FvgbyWfnlQAycxePrEF/eoAizeLwcbHPwkaV1yGNQSZChCQgwj1/mQLsTCIi+G/PT+XpuOnlaNYkq5ZD3IRBECyG04gQhtASnwmDv+Fu8Pqqdbqqcbg60OXgwghVqoCglyNSoydaoMamUfGxtDelLAbwxbxIzkvQkKAf+xxhNcQTRFTVjnVP9rCJZERXA63MnIUtQ7/ZS7/ZQ5/JQ5/ZQ7/IGnt1emiK4eXsjWaUgV6shR6smT6shV6cmV6smN/hag0HZFUcmLGOC4UAIrThBCC3BUOlpjVo7JZ8zU5OodnmodrqpdnqodnV/dtPqG5gdKkWlJEurJkerJksTeM7u/tCoSdeoUHZz88Sq9Ugwsgz1PPD4ZRrcHurcXj6xtPBguTmsTYJCQbt/4OdwjlaNDBxqdwZnnt5SmMn1eelka9Uh569AMFiE0IoThNASDBS7z0+V002l082+VvuQrVEpKiV5Og1paiWfWtvCrFAblkxlgkE3pDEKgSSA4bGMdbqQk1VKalwealxuapyekNfVLjc1Lg/tA7ihUEmQq9WQr9NQoNOQr1MHXxfoAtaynikwhGVM0J3h/P8WCUsFglMk0gW7u5AKeTgCz40DcK2ogHy9hjythjxdwK2Sr9OQq1WTp9OQp1WT0C2+JZL1YagiCyBXpxF/QIJTPg9ydRoeKS0IOzc7+5yiUjIlIfJ5KssyrT4/1U43HzS1cH9ZbVgbBeCVocLppqKXjPoSkKMNiK98nQabx8t6S2twNuUDk/O5Nj99yMcoEPSFsGhFEWHRGj84OoTU32oa+d+OJJoAhVoNdr9/QEIqUamgUKchQ6OKaI3aunQa+XrtoMYlrFCCWGW4LGObl05DIUlUdcQrVjrdVDk9gfeuwDLnANJgZKhVlBi0FOg0FHZYwgr1gedcrQa1IrJrUljGxgbCdRgnCKEVX/R1gfT4ZWpcbioc7uCdcqXTTYXDRYXTTb27fyGVpFQEXBf6LhdG94dRpQzmI4p20LhAMBYZyu9ElmUaPd6gVXmDtY0XapoGtV8FAYtYoV5DoU4bFGGH2508VVEvMvCPAYTQihOE0IofugefSwQSX5rUKiqcLiocbmpdnn5THugUUsQ75Uem5HNppgmjenCeeWGNEgj6Z7hmUz43qwSH30+lo+vGqvPhGkRi2CsyTUxN0HeIssBNVYZG1WeSV2EVG32E0IoThNCKDyJdaCOhU0ghFqhCvTbEpeDw+li0+WDU8kQJBIKRYTCWMb8s0+D2hli1K51u9rTa2dvmHND+9AqJfF2HNaxDgBV2XEe2NLcFy1QJq9joIYLhBYIoUuZwRRRZX8tOYUVqcshdqKKvKeNqVZ9BvgKBIDZZnZvGWalJA7KMKbplwl9kTAguj3TD1pmMtdnro6Jj4kuNy4PDL3PU7uJojxqUPfEDdx6upNrlZmaivivMYADWcWEVi12ERSuKCItWfBDNci6d/QmXn0Aw/hiIZczt91Pj8nTFe3ZYxCqcbo7Zndi8/aevMKqUEeM8O+M/365vDqtNKqxig0O4DuMEIbTiBxF8LhAIosGp3Gj1ZhU7NzU5GMA/mCz73ftYOyWf2UkG8jvy7A2kEPh4tooJoRUnjDehFe8/SmGJEggEo01/N33tPh9VTk9YLr7Ox0CEmF6hIL9bLr7OhK75Og15Og05GjX/Z7aMa6uYEFpxwngSWj1Lxoy3H6VAIBBEi1O56Ttmd3LGlkNhVrGZiTrq3V7qBpCKRgFhcasS8PupBcxKMpCrDdRAHctWMSG04oR4EloD+THIsozN66Omo95ejctDrcvDUbuTdxpsIW3FbDuBQCAYHfqyirn8fmp6JHANJnTtuK57BiADDEoFeSHFvruKgHc+r6uL31gxMetQEFV6WqNuyE+nNEHfTUy5O8SVB8cAi776gHKHSwgtgUAgGGH6mkWpVSgoMWgpMUSuMuGXZfa2Orhgx5GweqtTDTrqPR4sHh92n39AMyeD/RKYQdns8TItUU92R2F74wAsY/FqFesNYdGKIiNl0ervJPT6ZercHupcHmrdASuUueNx0uFiW4t9UPtLVSsDdzFaNbk6DYlKBU9U1If8KIVFSyAQCOKXvqxiDp+f2o6C351WsO6Fv2ucbloHUPwbAjnFsrVqsjusY50CLFsTeN7c3MZ9ZbUjbhUTrsM4YSSEVndrFMC0BB0SYHYH7jqGikaCAp2WHK2anA5TcE7HD0CnUISIKhmZ9U0t/G9VY4hV7OzU5B7tOp47TrFIJ1pIe7mr/7B1PdqELOs2rp59RdpX6D7l3tdF3Gfv7ZHDl/UcW1/jD91n759ByC57fGaRxh15PN376P/7GehnEHmfPdr38d2EtotwTH3us/f2hIy/j/Yj9Jn1v88I7XvsbLg+MyK26719n7/Jfj6D3vbTW/vw8URoP+DzK3S7XvsdwGcQ+fyK0H7A51f4tr2NP9K++v9N9mjfscAny3hlGZUkheQQHMhn1uLx4oz05Z4iI3UTL1yHAiBgyeousgAOtg8sM3F/uGU47nBx3DEws3B3/MBfqhr5S1VjVMYiEAgEAgGMjbAUIbTiiN4ymvckSakgTaNCiUR3V3jnS5vXF1IUOUujwtSRebi751zq+RzSl9RreyLsM2J7Kax5ePuIfYXT1VfX2kjtI+8zdGyRwgcitu9r/IP8DCLta6CfQdc++z72nvsc6GcQeZ+9fwY9t+s+tr7HH2mf4cdExPFHaC/1XBdpn9H7zPrfZ+/fZ6T3PT+zSP32+5n12MFQPoOe++rvN9bz+xnoZxD599rHdabHuEL672P8kffZ+3kTsq8e4+p13GF9Dey6FHFfQ7wO9HXedH8z+Gtz6H56rmtwB+J783RqMjTqPr+fnvusd3n4yhfHw5JJl+gjx5fFC0JoxRET9NqwabhDNauKHFICgUAgiDbT0A9520kG3ZgsayaEVhyRq9NE7STM1Wni/uQVCAQCwdhiMHUo4wUhtOKMsXgSCgQCgUDQyVgzBAihFYeMtZNQIBAIBIKximK0ByAQCAQCgUAwVhFCSyAQCAQCgWCYEEJLIBAIBAKBYJgQQksgEMQNTmctFuvnOJ21MddfrPYlEAhGFxEMLxAI4oKamv/j4KGfQUfhp2lT7yM396pe28uyjCy78fns+HyOjueu1w2NH1FT8xKBIiIS2VmXYzItCJQXkWUChUYCzyHvI6xradlNfcN/gn3l5FxJWuoZKBQ6FAotCqUOpUKLQqFDqdR1LNehVGqRJE1I8sfBHqdAIIhtRK3DKDJSRaUFgnjD6azF7jiBQV+MTpfTb3tZ9uP12nC7rXg8Ftrtxzl06GeEVmyTSE1dAcjdRFSXkPL7Hcjy0Ot/jhxSUIBJkgq3uyFsfW7O1ej1eajUJtQqI2q1CbXahEoVeFYqDSFirTuD/ewFgvGIqHUoEAjillALjURB/rUkJc/C47HicTfh9lg7Xls6XlvweJqh34JTMhbLpwMagySpUSoNKJV6lEoDsuzH4TgR1i45eT5aTVpHDRIJCUW32idSoISIpOgoJSKBBC5XI1brhrC+EhKmoVRq8fud+HxO/H4nfr87+L7r+GT8fgd+v6PX46ypfanf41OrjajVKai6CTGX04zFupFOS1tJ8ffJz/8manUKkiQiRwSCkUBYtKKIsGgJxhIDsYTIsozP14bTWYvLZcblMuN0mXF1vLc7KnE4yoc8BpUqCbU6BaUykba2Az3WSkwouQOdLjsoohRKAyqlAaXSgEKpR6noWK5Qhx3bxk1nECrmFJy2/NNBW32G0lfArenB73cFRZjP78TpqGL3nhvpabnLzb0GZB8ebzMeT+Dh9djweK34/e5BjRcCwkyryUCrzUKjzUKrzUSrzUaryUSrzep4ZKJSJYUdq7COCcYiwqIlEAhGlJ5WqPz8b5KQMDkooJwdosrlMuPz2Qfdf1LiTAyGYtSaFNTqNDTqFNQdD40mreO1CYWiKzFvNGOXdLocpk29L6y/oYiHofQlSRKSpEGh0ISImcSEyUybev+Aj1OWZfx+Z1B8ebwdAsxjpaVlLzW1L/eynQenqwanq6bPY1MqE9BqM9FoMvH7XbS07CZoHSu5jaLCG1AqDX32IRCMd4RFK4oIi5ZgtBmsxcHns+NwVOJwVuFwVOBwVNLedhRr86ZB7VelMqHTZQesItpsdNpstNocJEnNgYN3EmqhGZrlCALH53CcRK8viopFJZr9xVpfvVnali1dj0KhxOWqw+WqDwhmd33H+8Ayt7sOr7d1QPtRq9PQ6/PR6fLR6wvQ6fLQ6wo6luWiUGgjjk1YxgSxhLBoCQSCfolk8cnJ+TIuV11ATDkqcTgrcToqA6LKWYnb3Tjg/pOT55OUNL1DRGWj1WUHXyuV+l63k2VXVCxHELAeRfOPOZr9xVpfvVnaDIaCjvW5fW7v9bbj7hBgFstGTpx8MmI7j6cJj6epw9rVEwmtNisgwnT56PT5uJy11JpfIyC+xaxKwdhHWLSiiLBoCUYDv9+DzbaLnbtWE2o5AlADnj63V6mS0OsL0esK0esLUKqSKSv7LdGyQkH0LVGCgTOc1rEli99Blj04HFUBEe+swuGoCj73HuAfSkrK6SQllqI3FGPQF2MwFKPVZouAfcGIISxaAsEYpz9Xiiz7A8Hl9nLs9hPYHeXB105nZR9pDDxIkqqbOyfw0OkDokqvK0CtNoZtpdWkRc0KBdG3RAkGznBaxxITpwCQlDQjbBtZlvF4mrpEmKOaZtt2mpr+G9bWat0QNnNTodCh1xdi6BBf3UWYRpOJJEnCBSmIC4RFK4oIi5ZgKPR0+RUV3ojBUBIQUo4T2O3lOBwn8ftdvfYhSTpk2dljqYIF8/8Po3E2kqQc9LiEFUrQk1M9J3qzjJWU3IbHY8HhONlx81CFLHt77Ueh0KNWG3G5zB1LJIqKvkNR4Y0RbxwEgv4Yzv9vIbSiiBBagoHg87lwdAgoW8sXVFT8ZUDbSZI6/A7fUILBUIJWk0Vt7Ssio7gg5hnI7FG/34PTWY3dcQKH/UTHDccJHPaTOJxV9JVjTa1OJcEwMfDbSJgQfK3TFaBQRHbiCMuYQAitOEEIrfFJpIu0LMu4PU3Y28uw24/Tbi/Dbi/D3l7W7x8FQFLSLIzJcwOiylCCXl+MTpfX6x9F97EIK5Qg1jmV89Tvd2M2v83BQ3cParvOG5UEwwQMnY+EElpa9nL06H2IG5TxjYjREghilKrqf3L48C/ozDdlTJ6HjB+7vQyvt6XX7VSqJAyGCWi1OTQEa+R1omD2rKeGnNNJCCxBrHMq56lCoSE1dRmgoKcLcumSf+Pz2TviF8totx/veF2O3+/Ebj+O3X68j979HDz0U0CBybQQnS6/35sbgaA/xBkkEAwAj6c5cNFu77x4l9HWdhins6pbKxlby85u7yV0unwSEiZg6HBfBNwYE9Bo0oO16SK5UoRYEgh6p7fg/ISEiQAkJ88KaS/L/oDludOybC+n3X6ctraDeDzWHr3LHDx0DwCSpMFgKMJgmBiwhCVMDFrEVKrEsHEJF6QgEsJ1GEWE6zC+6HlRlGUfDkdV151we5fLz+OxDLjf4uJbycq8AL2+GKUyPFljb2MRLj+BYHAMT3C+hMEwEaezss8JKFptdkeMZEB8OZxVVFY+h3BBxiciRitOEEIrfghYkX5Kp8tOo8nC47Eiy73XjdNqcwIWqY4AW5XKxP4DPyQa9fIEAsHo0FtwfsAKVtMVYxm88To+wES/EgX515GcPJuEhIBFW5Qril3GhdB64IEH+OlPf8ptt93Go48+CgQCitesWcOf//xnrFYrS5Ys4YknnmDGjK6cLS6Xi7vuuosXX3wRh8PBueeey5NPPkl+fn6wjdVq5dZbb+XNN98EYNWqVTz22GOYTKZgm4qKCm655RbWr1+PXq9n9erVPPLII2g0XbXW+kMIrfgg8l1sAIVC23GXOiHo5uvrIhnN+nsCgWB0GKxlzOOxdbN8l2Ft3kZLSNhAZAJWsK6ZkJ3Xl+7JWYX7cXQY88Hw27Zt489//jOzZ88OWf7QQw/xu9/9jueee44pU6bwm9/8hvPOO4/Dhw+TlBQoxHr77bfz1ltv8dJLL5GWlsadd97JJZdcwo4dO1AqA7mDVq9eTVVVFf/+978BuOmmm/jmN7/JW2+9BYDP5+Piiy8mIyODDRs20NTUxLXXXossyzz22GMj+EkIRgK74wSRRNaM6b8nK+viQeWcys29itTUFcLtJxDEMYMNzlerjRiN8zAa5wG9uyAzMy8KJhr2eCzBQuxWa2gtUYVCh8FQgiSpaG3dR2fh7smT76Ww4NpTPj7B6DLqFq22tjbmz5/Pk08+yW9+8xvmzp3Lo48+iizL5Obmcvvtt3PPPYHARJfLRVZWFmvXruXmm2/GZrORkZHBCy+8wNVXXw1ATU0NBQUFvPvuu6xcuZKDBw8yffp0Nm/ezJIlSwDYvHkzy5Yt49ChQ5SWlvLee+9xySWXUFlZSW5uoP7XSy+9xHXXXUd9fX2v6tblcuFydfnwW1paKCgoEBatGKe3pInC5ScQCIZKf9Ztj8fa4YIsDwblt9vLcDgqkOXey2SpVKbAhBp9CQZDMfqO3HkGfVGvNUaFVWzwjGmL1i233MLFF1/Ml770JX7zm98El5eXl2M2mzn//PODy7RaLWeeeSabNm3i5ptvZseOHXg8npA2ubm5zJw5k02bNrFy5Uo+//xzjEZjUGQBLF26FKPRyKZNmygtLeXzzz9n5syZQZEFsHLlSlwuFzt27ODss8+OOPYHHniANWvWRPPjEIwAvc1YEhckgUAwVPqzbqvVKZiMCzAZF4Qs9/u9OJ2V1Ne/x/Gy34b16/U2Y7PtxGYLd01qtdkdZYlKggmM21oPUlb+B0Q4Q+wwqkLrpZdeYufOnWzbti1sndkcKK2QlZUVsjwrK4uTJ08G22g0GlJSUsLadG5vNpvJzMwM6z8zMzOkTc/9pKSkoNFogm0i8ZOf/IQ77rgj+L7ToiWIfYTLTyAQRJuh5AdTKFQYDCVkZ1/B8bLf09PSPnfO/+L1tnWU5CrHbj+J3V6O19vc5Yps3txL74G8YA5nJUlJMwMVJfQFAwrKF1ax6DFqQquyspLbbruN999/H51O12u7zlxDnciyHLasJz3bRGo/lDY90Wq1aLUDm74/VonnH6NI7ikQCGKF3iztaWlnRmzv8VgDBea7FZlvbTmAw3miR0uZEyeeDFmi0WSg1xdh0Bei1xei1xehNwTeq1QmUc4ryoya0NqxYwf19fUsWNBlRvX5fHz66ac8/vjjHD58GAhYm3Jyuv4M6+vrg9an7Oxs3G43Vqs1xKpVX1/P8uXLg23q6urC9t/Q0BDSz5YtW0LWW61WPB5PmKVLEMDnc1BR8VfKyh8lELgpfowCgUBwKgzG0q5Wp2A0pgQD8qH3oPz0tHNwuxuwO07i9dpwuxtwuxuw2baH9atUJuDztXdb4ufgoZ+h1xdjNM5BoRjfxoWhMGpC69xzz2Xv3r0hy66//nqmTp3KPffcw4QJE8jOzuaDDz5g3rzAieR2u/nkk09Yu3YtAAsWLECtVvPBBx9w1VWBP/ja2lr27dvHQw89BMCyZcuw2Wxs3bqVxYsXA7BlyxZsNltQjC1btoz77ruP2traoKh7//330Wq1IUJwLNKbRcrvd+N01uB0VuNwVOJwVuF0VOJwVuN0VkbIIxP4MaamrhBWIoFAIBgip2Jp780qFhqUb8PhOInDUYG947nz4XKZe4isTvzs3PU1QEKrzUKny0OvK0Cn73zOR6/LR6vNCZYsimdvR7QZ9VmH3TnrrLOCsw4B1q5dywMPPMCzzz7L5MmTuf/++/n4449D0jt897vf5e233+a5554jNTWVu+66i6amppD0DhdeeCE1NTU8/fTTQCC9Q1FRUUh6h7lz55KVlcXDDz+MxWLhuuuu4/LLLx9Ueod4yqPl93uoqHiG42UP0zmVODl5LgpJhcNZhctVR3+FjyMxf94/SElZGu3hCgQCgWCADDVjvs/npNm2gy++uJbQ+qsgSTpk2dnn9pKk7MgJpsHhKO9cSmHhDeTlrkany45Zi9iYnnXYFz/60Y9wOBx873vfCyYsff/994MiC+D3v/89KpWKq666Kpiw9LnnnguKLIB//OMf3HrrrcHZiatWreLxxx8Prlcqlbzzzjt873vf47TTTgtJWBrL9HbHIMsyXm9zh0Wq4+HqfF2Ly1mDy93TnSrT0rIrZIlCoUOny0evzw+7e1FIGrZsvZiegZt6fdGwHa9AIBAI+meoVjGlUkda6mlMm3p/mFUsJ+ereDyWDu9GVTcvRxVOZxUORzWy7MbprO7Rq0xFxf9SUfG/AGg06ei0uWh1ueg6H9rcjjHnolanBWOjx4pVLKYsWvHOSFq0QnO2RPKrDw2TaSnG5Dno9Plo1GnQcy5At7PFYtlIdc2LdP4Y83K/Rmrq8l77lunrVOtjXb+naO/r+9znkE/9vrcb+j6Huq6/ffa95fDsc2j99n2O9MOw7LOf8Yjj7Gef0e+z/30Ox+/vFNae0j6H45o5sO/L623D67WhUhm7FdDuY1vZh9vdQGvrAdrbj/Sx/6Ew/DHA46IEz1hgpIRWXyVkBAKBQCAYewxvUulx6zoURKa3EjKdqNUpaLWdvvDe01N43E04nBXB93pdIRpNWuBNnyk0hrZO6mu7Ie+vb4Zjn332OeT99dVjP9uNge9q6Pvrb5/DcYxDG8twfaZDPx/77jX6+xuuY+xrw1g6/4f6mfbTbQz9xjv7bGs7jMW6kc7439SU00hMLO13f36/B5erjva2I9gdZT0a+XE4TsalC1EIrTjEoC8GFESjhMxQgyYFAoFAIOiNU/lv6a1MWrzGACtGewCCwdM5hbfr6xt6CRmdLoeUlKVCZAnGLTabjfLycmw2W0z2JxDEI6fy3xLN/7hYQMRoRZGRTu8grFGC8YjNZsNisZCamorRaDylvnbu3Mlbb70VrAJx6aWXMn/+/JjpL5rHKhDEGyP5HyeC4eOEeMqjJRDEI/0JGa/Xi8vlGtCjtbWVgwcPhu0jNzc3JD0MBFKm9KTnMq/XG7EKxaxZs0hMTAyW7NJoNCHPPV+rVKoBHatAIIgeQmjFCUJoCQSRGYplxu/3Y7fbaWlpobW1lfr6ej766KOwdkajEY/Hg8vlwufzRXvoI45CoUCj0eB0hieHPP3008nIyCAxMZGkpCSSkpLQ6XT91n/tRFjIBILIiFmHAoEgbulpmbnkkkuYMWNGUEB1f3Rf1tbWht/ffwqTSLFQarU6aCHq7SHLMhs2bAjZTpIkLr74YgwGQ0Tx0tcyu93OunXrwtYvXboUSZJwu924XK7gc/fXbrcbj8cDBARmJJEFhI0XQKVSBUVXUlJSiAjr/ti/fz9vv/22sJAJBCOMEFoCgSDqOBwOmpqaqKqq4t///ndwuSzLvPXWW8HyVwOhUzjodDrKy8tD1kmSxFVXXUVqamqIG66n6683UlNTo+qe6zy+ofTn8/lwu9243W6ampp44YUXwtyT06ZNC7o9W1tbcTqdeL1erFYrVqt1UON88803kWWZ/Px8UlJS0Gg0gzpWgUAwMITrMIoI16FgLNGfm8lut2OxWLBYLDQ1NQVfWywWHA5Hv/3rdLqgtSU5OTmiFSYxMTFENA1H3FK03WnR6m8gx+rxeGhrawuzDHZaBLsLsv5ISkoiNTU17JGSkoJOpxu24xQIYgERoxUnCKElGCv0/JOfO3cuycnJIaKqvz/vTgFVXR1a+0ySJG655RbS09OHNLbx9AcfrWNtamri8ccfD7OQZWVlYbPZ+v0uExISQoRXS0sLu3btEm5IwZhBCK04QQgtQbxit9upr6+nvr6eqqoq9uzZM6DtOq0gaWlpYZaQTleUmD0XG/T1PXS3Tlqt1hDrZHv7wGqozp8/n4KCAjIyMsjIyECr1Q7n4QgEUUUIrThBCC3BaNOfBcTpdNLQ0BAUVfX19TQ0NNDW1tZv35MnT6aoqCgoqgYT1zOerFCxzFC+B6fTGSK+Tp48ybFjx/rdzmQykZGRQWZmZvA5PT097JwR54YgFhBCK04QQkswmvS0WKxYsYLU1NQQQdVXtnKTyURmZibJycls3749ZJ0kSdx+++3ij1CAzWbj0UcfDXFDSpLEnDlzsNls1NfX92kFS0lJITMzk8zMTNrb24ULUhATCKEVJwihJRhpZFmmpaWFsrKyiKkFIpGUlBT8o+tuceju6hHuPkFf9Hd+dLqiu1tPGxoasNvt/fa9YsUKiouLyc7OJiEhYTgPQyAIIoRWnCCElmA48fl8NDQ0YDabg4+6uro+Z/hlZWVRVFQUIqr0ev2A9idcOoK+GMr50dbWFhRdx48f58iRI322T0pKIisri+zs7OAjNTUVhSK8TK84XwWngkhYKhCMcXr+STgcDurq6kJEVUNDQ8TM5wqFgtTUVBobG0OWS5LE6tWrh/ynYzQaxR+WoFeGcn4kJiaSmJjIhAkTmDp1KkePHg2bCTl58uTgzNbO9BTdY8LUajVZWVkhAqy2tpb33ntPWGAFMYkQWgLBKCLLMhs3buTDDz8MLtPr9b1aqbRabcjdfXZ2NhkZGahUqojuHCGUBLGK0Wjk0ksv7dUF6XK5qK+vD7PgejweqqqqqKqqithvZ9LYnJwccnKGtxCxQDAQhOswigjXoaAv/H4/FouF2tpazGYztbW11NTU9JrDyGg0hokqk8nUZ1074T4RxBuDOWf9fj9NTU0h1t6qqqpef0NJSUnk5OSQnZ0dFF5Go3HAtSEF4wcRoxUnCKE1Pon0R+H1eoPxVLW1tdTW1lJXV4fb7R5Qn1/72tcoLS0dzmELBGOCSLMg+0Kv1weFV+dzWlpaMO5L3KyMT0SMlkAQo3R31wEUFhbi8Xior6+PGE+lUqnIysoK3l0nJiby0ksvhU2Vz87OHrFjEAjimd5ckDNmzKCuri54o2M2m6mvr8fhcFBeXh5SN7Mz7kulUnHixAkAEesliBrCohVFhEVr7NPe3h60UlVWVnL48OFe22q12qCg6n7n3LPgsUilIBCcOgOxRHm93mDcV3dLs8fj6bXf6dOnU1hYGHTfR6r7KIh/hOswThBCK77o68IsyzLNzc3BO+HOC3Nra2u//Z511lnMmTOn33iqgY5FIBAMH51xX3v27OGzzz7rt31KSkqY6zEpKSmkjfg9xx/CdSgQRJmeVqQzzjiDlJSUoKAym824XK6I26amppKdnU1KSgobN24MWSdJEvPmzRv0xVWkUhAIRgeFQkFGRgYLFy5kw4YNYW78JUuWYLVaqa2tpaWlBavVitVq5eDBg8F2CQkJQdHlcDjYuXOnsFALggihJRh32Gy2kLgqWZb55JNPwtopFAoyMzODd66dj+4Z1NPS0kRKBYFgDNBfugkIZLzvfjNWW1tLU1MT7e3tHD9+nOPHj4f0Kcsyb775Jm1tbRQXF5OVlSWKbY9DhOswigjXYXxQXl7O888/H7Y8KysrWPojJyeH9PR0VKr+70WEm0AgGDsM9vfsdrupr6+ntraWo0eP9pvtPiUlJZhstTPxqslkEtnuRxnhOhQIokhqaiqSJIW5CIaaRV24/QSCscNgf88ajYb8/Hzy8/OZMmVKxGz3xcXFNDU10draGnQ9Hjp0KKSP7tnus7KyRLb7MYQQWoJxR28uAiGWBALBqdCf+7G9vZ26urrgo7O0ltvtprKyksrKyoj9dma7NxqNFBUVDcjSLogdhOswigjXYXwhzPICgWA4GMy1xefz0dTUFCwxVFdX12e2+87appmZmcFHRkYGqampYaljxDVu4Ij0DnGCEFoCgUAgOFV6y3av0Wh6rS6hVCpJT08nIyODzMxMWltb2b59u3A9DpCYEloul4utW7dy4sQJ7HY7GRkZzJs3j5KSkqgOLB4Zz0JL3DkJBAJB9IiUyHjevHm0tLTQ0NBAfX198NHQ0NBn0tVOli5dSl5eHunp6aSlpaHRaEbgSOKDmBBamzZt4rHHHuONN97A7XZjMpnQ6/VYLBZcLhcTJkzgpptu4jvf+U5Y8rbxwngVWiKzuUAgEESfgd7A+v1+bDZbUHSVlZVRVlbWb/9Go5G0tDTS09NDHklJScFky+PlJnrUhdZll13Gtm3bWL16NatWrWLhwoUYDIbg+rKyMj777DNefPFFdu/ezd/+9jfOO++8qA40HhhPQsvlcmGxWKiqquKdd94JWSdJErfffvuY/lEKBAJBrNKb63HmzJm0tLTQ2NiI3W7vdXuNRkN6ejoKhYKqqiogcF2/4IILWLJkybCOfbQYdaH1xBNPcOONNw7IzLh//35qamqE0IojodXbHYvb7cZisWCxWGhqagp53dbW1mef1157rXAnCwQCwSjRn6fBbrfT2NgYfDQ1NdHY2IjFYgkTaN3R6/Wkp6eTkpJCampqyKO7AaaTeLGIjbrQEgyMeBRaW7duDeZqASgsLEShUARzvvSFXq/HaDRiNptDlguLliDe8NpceBsdqNL1qIzRydw9HH0KBINhKCLH6/VitVrZv38/H3/88aD2p9PpQoRXS0sLu3fvjouwkpgVWm1tbfj9/pBl8SIwhoPRElp9/Zi8Xi82mw2r1Upzc3PIs8Vi6XUKcSd6vT74o0lLSwt51uv1gIjREow8QxUxsteP3+ENPOwe/HYvjgNN2LfXBdvopqeiKUgCJOioCR4IV+l6H3zusUzqaOyqaMHxRUNwedLZhSQsykKRoEZSKwZcbDwaxywQDIVI7sfOxM7dvR2dj/5uzDspLi4mIyMDk8kU8jAYDL3+LkbCKhZTQqu8vJzvf//7fPzxxyF/0p1/sj6fL6oDjCdGQ2h1FzkAU6ZMQa/XBwVVS0vLoPucP38+EydOJC0tLSim+qO1tZXm5mZMJtOQJkMM5Y8nWozmvkd7//F47PZd9TS/XQYyIEHi0hw0BUn47V78Ti9yp5ByePA7fIHXTi9+hw/Z3XV9GrUjVylQ6lUoEtQoDGoUemXg2aDuWK5C0qs71itR6NU49zXS/FbXMZsunYhhfuaQhxCP3/tY2Pdo73+w+961axfvvPNO8P/94osvZt68eRHbejye4A28xWLh5MmTHDt2bMD7UqlUYeLLZDJRV1cXLPY9nDfyMSW0li9fDsBtt91GVlZW2Bd35plnRm90ccZIC63eAh4FAoFAIBiLDFdoSkzVOtyzZw87duygtLQ0qgMRDJ7eghYVioG7Jvx+f5hpOFJx0/4YbbE32vuPlTEIBALBWEaWZSwWS1zFAA9aaC1atIjKykohtGKA3ooj33bbbYM6CeNlVohg+ImGWOzZR9uGamzvlQfcXoAq0wB+Ga/VCb7w/XVfokjUoErVoUrTokrVodCraH6nPLSRBNk/WhSMWxrsMdi312F942jQLZdy+WQMC7MG1UdP2rebsb5xrFufkzAsGHyfsizjbbRT/8cvQo8ZUOcn4rO58LdGTlQpd2wg6VSo0nSoUvWBzzFNj6fBQdtnVV3ju2xo4+tv7KONGMPo77+lpQWr1UpKSsqgw0paW1v505/+FPYfl5qaGu1hDiuDdh0eP36c73znO3zjG99g5syZqNXqkPWzZ8+O6gDjidGO0RKB6ILRRPbLeJsceGra8dS24a5px13Vimz39r6RQkKVokWZpu8QVJ2CQIcyVYdCowzbpH2bGetr3YTRlZNJWJR9SmOP9VmHfR2z3+3D2+TE1+TA2+TA2+jseHbga4lcriUSpssnopuUgjJVh6QY3TgmgaCTkfqPi6kYrc2bN7N69WpOnDjR1UmHVUUEw8ferEOBYKj0JRRkjw+P2Y67pg1PbTuemjY85nZkt7+X3kJJvqgYw8wMlEYtklLMwBsIQzlmv9uHz+LE2+jA2xQQYO7KVjy17b1uI2mVqHMS0OQmos5LRJ2biDpTj6QcfEiBQBANxt2sw+nTpzNt2jR+9KMfRQyGLyoqiuoA44l4zKMlEESipwUl8bRclMnagLCqacfbYA9zZQFIagXq7ATUuQmocxJRJKuxvHAw3NX348XjRiDFGl6bC/ODW8O+P1W2AW+jA7wRvliVhDq7Q3zlJgTEV3ZC0OI4HoWvYGwRU0IrISGB3bt3M2nSpKgOZCwghJYg3vG7vDgOWbC+eLjftooEdfBPV5MTeFal68PcTsPh6hOcGr19J7LPj7fBgbu6DU9NwP3rqWlDdkXwVEigyjCg0CpwV7YFl4nvVxCPxNSsw3POOUcILYFgDCDLMt4mJ+6TLbgrWnBXtOIxt0e0VAFoSpLRTUlBnZOIJjcBRZJmQLNbExZlo52SIiweMURv34mk7LBIZidAR3C87JfxWZ0Ba2Z1e4dVsw1/mwdvfY96eTJYXz2K60QLuikmNEXJAffwKOeuEghGk0ELrUsvvZQf/vCH7N27l1mzZoUFw69atSpqgxMIBNHD7/bhrmzFXdEaEFeVLfjbwwPVFcka/D2DqCVIvWbqkEWSyqgVAivGGOh3IikkVGl6VGl6mJUBBES6v9VN+456Wv5zImwb+4467DsC2fYVyRq0RcloCpPRFCWhyU1EUol4L8H4YdCuw75yLIlgeOE6FIwuwViZNB34ZFydoqqiJWCt6hmrrpLQ5CWhKUxCU5iMtigJZbJWuPsEAyJivJcEhvlZeOoCbsdez7miwPmmKUxGmaQJ6VNYPwUjTUy5DnvWNhQIBKOP7PVj+/AkbR9X9dlOadQELAv9WBeEu08wEFRGLSlXTu4z9YSnqg1XRUtA8J9swW/3Bl93RHahTNWhLUpGlmUcuxuEwBeMKU6pqLQgFGHREowUfrcPd0ULrvIW3OU2XBUtEWeLqXMS0E4woilKRlOULASTYFgYqBUqkIDVEXRfu062BOK8+vgXSv3GVPRT04S7UTCsjLpF66WXXuKaa64ZUIeVlZVUVFRw2mmnndLABAJBF367B9eJFlwnAsLKXd0G/v7vkYyXTEA30TT8AxSMawYc7yVJqDMMqDMMJHQE2/udXtwVrbTvqsexqz5sG8vfD4FKgbYwCU1xMtqSwI1DpGS2AkEsMiCh9dRTT/HLX/6S66+/nlWrVjFt2rSQ9TabjY0bN/L3v/+dDz/8kL/+9a/DMliBYKzS0yLga3HhKm/BdcKGu9yGpy78rl9p1KApMaItNqJK19H4131hsTKqdP2IHodAMFgUOhW6KSmosgw4vqgPO88lvRLZ4cNVZsNVZqOVSlBIaPIS0ZQkoy02oi1ORmEInZglYr0EscKAXYdvv/02jz32GB9++CEJCQlkZWWh0+mwWq2YzWYyMjK4/vrruf3228nMzBzuccckwnUoGAptW2tpfv1Y8A9GSlAhR5gNqMrQoy02Bv9clCmh0+ZFALsg3ol0DhsWZuFtcOAqt+E+0YKrzIbP5grdUAJ1liFw41FixNfs6qqvKX4LggEQUwlLm5qa2LBhAydOnMDhcJCens68efOYN29enzMSxwNCaAkGguyX8dbbcZ2w4TxsxXnQErGdOjehQ1gF7ti7z8zqDXEXL4h3BnIOe63OoBvdVW7D2+Dou1MJsu9ZhMqkG4YRC8YCMSW0BL0jhJYgErJPxlPbhqvcFgheP2HD31ehZSDt2unop6WN0AgFgvjG1+oOCi/HYQu+JmdYG4VBhXZyCtqSQJyXKsMgimcLgox6MLxAIIhMpLtv2ePHXdWK60SnsGpBdofml5PUCjRFyaizDbRtrAmLrVLnJo7gUQgE8Y0ySYNhVjqGWekk2vIj1nL02704djcE0kcQEF6aYmNQeKlzEsMKnAsLsSAaCKElEAyRnvEk2tIUZFcg+3rPVAuSThkI2i0JxFhp8hKRlAFXuzorISwuRVzUBYKhESm3l3HVRNSZhoCr8URXPi/ngSacB5oAkDRKNEVJaDvivDx1dprXHRNxXoJTRrgOo4hwHY4PfC1uHAcbaX79eK9tFInqwAW7OBlNiRF1dkKfbgpx5ywQRJe+flOyz4+7ug13eUvApX+iBdnZtzsfCbJ/vFj8PscoIkYrThBCa+wRDFw/GXABuk624LOEx390krgij4TF2ajS9aKQrkAQJ8h+GU+dPRhc7zxqRXaGl5NTpGjRT04J5PMqSkaZqhO/8zFCTMZoud1uysvLmThxIiqV8EAK4o/I8VU+3JVtuE52TCU/2Rp+pysFUi146x1hyxNPzxN3vAJBnCEpJDQ5CWhyEkhcnou32Yl57bbwOC+ri/atZtq3mgFQJKkDBbOLAtZrdW5CMCSgE2GtFgxaIdntdn7wgx/w/PPPA3DkyBEmTJjArbfeSm5uLj/+8Y+jPkiBINqExFcB2skmZKcPd00b+HrEV6kVgaLLRYH8VZrCJBQ6VcScP+JCKhDEPyqTLjzO66ISVCm6gHX7ZAvu6jb8rR4c+5pw7OuI81IrUOcnBUIGipLxNjmwvV0m4rzGOYMWWj/5yU/YvXs3H3/8MRdccEFw+Ze+9CV+8YtfCKEliHm8NleIyAJwHW0OvlYkaYIXSm1xMuqc8LtUEIWXBYKxTG+/b/3MdKDD+l3VFgwrcFd0FMwuD1RzCEMG62tH0U4yoUoR+bzGE4MWWm+88QYvv/wyS5cuDfFNT58+nePHew8OFghiBW+jI2IR28Sz8klclD2ouIuB1ngTCATxR1+/b0mtDM5QhI54zkZHoGzWiRacx6z4WzyhG8lg/u2OQN3Ggq6HUlxDxjSDFloNDQ0RS+y0t7eLoEBBXKBK14NEWO6qxGW5QjSNY1qbGrHW1pCSk0tSWnrM9imITSSFhDrTgDrTAItz8NpcEfN54fUH6zZ2okzWoClIQt0pvvITUWi7/p5FnFd8M2ihtWjRIt555x1+8IMfAATF1V/+8heWLVsW3dEJBMNApDw7Ir5qfLN3/ft88OfHkGUZSZI476YfMOuc82OuT0H8EOk6Y7piEtrCZNyVrcGHx9weSBmzvwnH/kCsFxKoMg1oCpKQfXJXsW0R5xWXDDq9w6ZNm7jgggv4+te/znPPPcfNN9/M/v37+fzzz/nkk09YsGDBcI015hHpHeILcZcYvwzGUuTzerG3NGNvbqbdZqW92drtdTMtDfWYjx0O206t06FQKpEUShQKBVLHQ9H9Wep6Tbd1fp+P+vIeoRSSxIprriWjuISktHSS0jLQGgzR/FgEMUh/1xm/24enqg13VZf48jW7IvTURfIFxegmmlBnG5DUyuEa+rgi5vJo7d27l0ceeYQdO3bg9/uZP38+99xzD7NmzRpUP0899RRPPfUUJ06cAGDGjBn8/Oc/58ILLwRAlmXWrFnDn//8Z6xWK0uWLOGJJ55gxowZwT5cLhd33XUXL774Ig6Hg3PPPZcnn3yS/Pz8YBur1cqtt97Km2++CcCqVat47LHHMJlMwTYVFRXccsstrF+/Hr1ez+rVq3nkkUfQaPov5NuJEFoCwfDT01I076LLyCqeQHuzlXZbM/bmgJhqb7ZitzXjaG0Z7SH3ikZvCIiu9AyS0zKCrwNCLCDGVBGuQcIlObbxtbhxV7Vi39OA44uG3hsqQJVhQJObiDo3EU1eAurcRBS6UGeVuKnsn5gTWtHirbfeQqlUMmnSJACef/55Hn74YXbt2sWMGTNYu3Yt9913H8899xxTpkzhN7/5DZ9++imHDx8mKSkJgO9+97u89dZbPPfcc6SlpXHnnXdisVjYsWMHSmVA6V944YVUVVXx5z//GYCbbrqJ4uJi3nrrLQB8Ph9z584lIyOD3/72tzQ1NXHttddy5ZVX8thjjw34eITQEgiij8veTsPJchpOllN96CCHP/900H1ICgUGo4kEYwoGU+A5wWTCYExBoVSy/rmnodulUJIUXPWL+9EnG5H9fmS/H7/fD7KM3+8Lvu++rvtru83KB395IqRPgMKZc3C02GhtasTZ3jagseuTjUHRlZSWTnuzhaNbPwfhkhzz9BbnpSlOxtvgwN/uibidMlWHJi8RdW4CvjYP7ZtqhOuxH2JKaL377rsolUpWrlwZsvw///kPfr8/aI0aKqmpqTz88MN8+9vfJjc3l9tvv5177rkHCFivsrKyWLt2LTfffDM2m42MjAxeeOEFrr76agBqamooKCjg3XffZeXKlRw8eJDp06ezefNmlixZAsDmzZtZtmwZhw4dorS0lPfee49LLrmEyspKcnNzAXjppZe47rrrqK+vH/CHPp6FlrhjEpwqsizT0lBH/clyGk6UBcWVrb6u320zSyaSll8YEFOmFBJMKSGv9YlJSIrwFB2d7F3/Ph/85XFkvx9JoeC8G78fnRitPvp0Ox20NjbS2tRAa1PguaWx83XgvdfVtwupk1nnnE/BjNlkTZhMSnZOn8cqiC8i5etLWJSNLMv4W9y4q9vw1LThrmnHU9PWr9sRwHhhMdoJJlSZBhRa4XqEGMsM/+Mf/5gHH3wwbLksy/z4xz8estDy+Xy88sortLe3s2zZMsrLyzGbzZx/fteFSavVcuaZZ7Jp0yZuvvlmduzYgcfjCWmTm5vLzJkz2bRpEytXruTzzz/HaDQGRRbA0qVLMRqNbNq0idLSUj7//HNmzpwZFFkAK1euxOVysWPHDs4+++yIY3a5XLi6XQhbWmLXRTGc9HYhEAggspvL43bRVFlBw8ly6jtEVWPFCVz29oh9JKVlkFFUTHJmNl/85+1Q65NCweV333tKLrRZ55xP8Zz5NJtrMGVHxx3XX58anZ60/ALS8gsibi/LMs72NlobGwJirLGR6kP7ObQp3KK3d/377F3/fqBfvYGskolkTZxM1oRJZE2YhCkrR8wKj1N6y+clSRJKoxa9UYt+elqwva/dg6e2DU91O44jFtzHw3N62d47EXytTNWhzjJ0PBJQZQVmTkoqkeE+WgxaaB09epTp06eHLZ86dSrHjh0b9AD27t3LsmXLcDqdJCYm8vrrrzN9+nQ2bdoEQFZWVkj7rKwsTp48CYDZbEaj0ZCSkhLWxmw2B9tESkeRmZkZ0qbnflJSUtBoNME2kXjggQdYs2bNII94bOB3evE2BHLG2N4p71rRmZRvSor4MQpC4qlAInvSZDxOJ5bqKmTZH9ZeoVSRVlBIZlEJGUUTyCgqIaOoGH1S1x1mZlFJmKUoGsKoMy4qmpxKn5IkoU9MQp+YRGbxBAAmLlzC4c8/I8QRIUlMP+McrDVVNJwox+2wU3lgL5UH9gabaBMSyCqZRNbEyWR3iK/kjKyg+BIxX7HNYPL1KRPUKCeloJuUgn5uRu+uxyYH/lYPPosTn8WJ86Clq4ECVGl61NkJqLMM+FrdgbJD4mZ6SAxaaBmNRsrKyiguLg5ZfuzYMRISEgY9gNLSUr744guam5t59dVXufbaa/nkk0+C63vehXUGwPZFzzaR2g+lTU9+8pOfcMcddwTft7S0UFAQ+e401olY96/DNO1psOOtd3Q82/E0OPC3uHvvTA4kBRVCa3zistupPrSfY9u3sPejf3dbI2M+diT4Tp+UTEZxQExlFpWQUTyB1Nx8lP3UTh0O61O8kJSWznk3/aBXl6TP68VSXYm57Ch1x49RV3aUhpPluNrbqdi3m4p9u4N96RKTyJowCYVSSfkXO0TM1xikt1Q2nSLJ1+7BW9eOp86Ox9z5bEfuuJH2Njhw7O3RqQzWV4/iqmoNxIGlG1Bl6FEkqvv8vxzPFrFBC61Vq1Zx++238/rrrzNx4kQgILLuvPNOVq1aNegBaDSaYDD8woUL2bZtG3/4wx+CcVlms5mcnJxg+/r6+qD1KTs7G7fbjdVqDbFq1dfXs3z58mCburrwGI+GhoaQfrZs2RKy3mq14vF4wixd3dFqtWi18X/CtG2upXndseBdj6YwKZDluMGB7AqvYN+JIkmDKkWLu6I1dIXUkRRUMC7wOJ1UH9pPxYG9VO7fQ13ZMWR/uLWqkxVfu5bpZ5xDQkrqkN1Zw2F9ihf6EppKlarDCljCrLM7xZeHxsoK6jrFV/kxGk6ewNnWysk9u0L6lmWZ959+DEmhYMqS5Wj0Iv1EvNNXqTBlghrlBBPaCabgMlmW8be68ZjteOracR5txnXEGtavfYsZe7f3klaJKkOPOl2PKsMQ2FeGHlW6HsfuhnEdXjJoofXwww9zwQUXMHXq1GAKhaqqKlasWMEjjzxyygOSZRmXy0VJSQnZ2dl88MEHzJs3DwC3280nn3zC2rVrAViwYAFqtZoPPviAq666CoDa2lr27dvHQw89BMCyZcuw2Wxs3bqVxYsXA7BlyxZsNltQjC1btoz77ruP2traoKh7//330Wq1cZMXrN9cLS4v3iZnx8OBz+LE2+jA0xhunQoRTlLAhKzK0KPKNKDOMKDK1KPOMKDQB04fUVx5fOFxu6g9cojK/Xuo2L8X87HD+H2hgtyUlUP2xMkc+vyzsHiqaSvOJjE1rWe3gkEwGKGpVKkDMVslE+HcwDKvx0NjxQkObviYne+u67GFzH+eepQP/vwY2ZNKKZo1h8JZc8mZVNqvtVEQmwzG9ShJEspkLcpkLbopKehnR3A/SmBYkIWvxY230YHP6kR2BfKBear6mU3bEV6iSNagLUwO/o+MZYaU3kGWZT744AN2796NXq9n9uzZnHHGGYPe+U9/+lMuvPBCCgoKaG1t5aWXXuLBBx/k3//+N+eddx5r167lgQce4Nlnn2Xy5Mncf//9fPzxx2HpHd5++22ee+45UlNTueuuu2hqagpL71BTU8PTTz8NBNI7FBUVhaV3yMrK4uGHH8ZisXDddddx+eWXx0V6h55CJ2FhFkqjNiCqLAFh5W+LPA24N5K+VIhhVjqqNH1YUGQkxrNZeCzSPWZHn2zEfOwwFfv2UHlgD7VHD+PzhJ5PSekZFM6YQ8GMWRTMmE1yegYwPLP5BNGjtamRv9xyPT3/BpLTM2lprA9ZptbpKZg+k6JZcymcNZe0/EIRYD9O6G/Ck+zx47UEXI2eRkeH29GOt9GB3+7ts29Jp0KVpkOVqkOZqkOVEnitStWhNGmD/z/D/R8TU+kdoskNN9zARx99RG1tLUajkdmzZ3PPPfdw3nnnAV0JS59++umQhKUzZ84M9uF0Orn77rv55z//GZKwtHuslMViCUtY+vjjj4clLP3e974XlrB0MK7B0RBavdbTOlV0yh4xbP20D1sv9bN+6Nv213xQ+460vq+D7Xcs/exsOD/HsN0P7XNwtLbQ2tiI3HFSSZLU8UfcdZIplCo0ej0avQGNXo9Spe71c/R7PXjdHlQadaDdQMYzzJ/joM/nQZwT4ev7btD/WE7hsxjAti0N9dSfON5heZTILJmIMSMTj8uJvcVGu82Go6UZn9dL93NAqdZgMBo78pOZUGm0IIHX5cLtdKLR6VD1vH4O4lii/9s6te3DROWgrmnDPJbBbD/YbTve+10+/A4vCr0qpA5jf5cZr83Vd9LVfsaiNAbElrfREVw2HK7HmBNaH330ER999BH19fWBJH7deOaZZ6I2uHhjNISW83gzjX/pGa0oEAgEAsEYRYLsHy+OqmUrpvJorVmzhl/96lcsXLiQnByRm2W0UaXrA7cQPfznGTfNRpnUVbojTE330Ne+Vje+ZhdKkzawXX/yuz993t9qubc3g+/r1McywvsP/zIG3FVf20ZeH7m5LMu0WRqpPXoY87EjNFSe7GXngd/3iq8FavQN+rMYxuMJH+7I7XtAzUfwswhvOrz79nm92OrrsNZUYampClhAQzLrB18x/6JVGDN7WB9G8LPouSDqv7FTvD70P55TuD4N8lj6a34qn0XAIuZBoVP1niQ1wvZemwtnZ7Htbu3iaWb7oIXWn/70J5577jm++c1vDsd4BIOkt+m72hLjoPpRZ4jZRWMdn9dD1YH9lO3cStnObTTX1YasT80rIH/aDPZ89J+wAPaUhUXo01JHesiCGMZIIYUsAsDZ1sauf7/Jplf+Gdbu6Es7SM3NZ8KCxUxcsJjcKdNQKEU2csHA8NpcmA80hRkT4mlm+6CFltvtDs7WE8QGfU3fFYwfIiWdbG+2Uv7FDsp2buXknl24HY5ge4VSRcGMWUyYv4gJ8xZhyg7MuM2eOGVYEoIKxi66xERmnn0+n//rxTALsaRQYOmwfG1/6zV0SclMmLuACQuWUDxnPlqDuMkT9E5vxoR4+p8bdIzWPffcQ2JiIvfee+9wjSluGc+1DgWjS0gGdkli4oLFtDdbMR8/GmKdMhhNAWE1fxFFs+b2mieptalxXCYEFZwakWaZTll6WkDs79hK+a7tIcW0FUoV+dNnMnHBEiYuWIwxs/e8hYLxzbiadXjbbbfxt7/9jdmzZzN79mzUanXI+t/97ndRHWA8IYSWYDSwmmt55vabeo3lyCyZyIT5i5k4fxFZEyaJgsOCYaUvke73+ag5fJBjO7ZQtmMr1trqkPXpBUUdLsYl5EyaQpvVIkoDCUaEmBJavRVYhsCU0/Xr15/yoOIVIbQEI4WzrY3yXds4tn0Lx3dsxecJL4m04JIrWHjx5SI5qCBmsdRUU7ZjC8d3bqX60IGQigIavT7o6halgQTDTUwJLUHvCKElGE5aGus5vn0Lx7ZtpurgvrBs7N2RFApufPwZYQUQxA2OtlZO7NrO8R1bKdu1DY/TGdbmzG98m+lnnosheXCTfQSC/hBCK04QQksQTWRZprHiBMe2bebY9s3Ulx8PWZ+WX8ikRUuZtHAp9SfL+PB/nxQZ2AVjghO7d/Lq/T+PuE6SFOSWTmPSwiVMXLSUlOzcER6dYCwSc0Jr27ZtvPLKK1RUVOB2h7osXnvttagNLt4QQkswFLrPFkwwpVB9aD/HOixXLQ3dCqJLEnml05i4cCmTFi4hJScvrB8RwC4YC0QsDSRJpOUV0lR1MqRt4IZjGZMWLhExiIIhE1MJS1966SW+9a1vcf755/PBBx9w/vnnc/ToUcxmM1dccUVUBycQjHVCZgsCKq0Wr8sVXK9SayicPZdJi5Yycf5iDEZTr30NptCwQBDLJKWlc95NP4hYJ7O7C73ywF6aqipoqqpgy+svk5iaxsQFS5i0cAkFM2cHSz5FSn0iEIwUg7ZozZ49m5tvvplbbrmFpKQkdu/eTUlJCTfffDM5OTmsWbNmuMYa8wiLlmCgtFoaOfDpf9nw4vNh67SGhICwWriE4tnzUet0ozBCgWD06c9K231SSPkXO/A4u/LEafQGSuYuQGMwsG/9+8iyLILqBb0SU67DhIQE9u/fT3FxMenp6fz3v/9l1qxZHDx4kHPOOYfa2tr+OxmjCKEl6A1Zlqk/Ucbx7Vso27mVurJjvbb9yv/7DUWz5o7c4ASCMYDX7aZy/x6ObdvM8R1baG+2RmwnSQpufEJMFBGEElOuw9TUVFpbWwHIy8tj3759zJo1i+bmZux2e1QHJxDEM50X/uM7tnJ851bamhq7VkoSmUUl1J8op3ttCUmhIDU3f+QHK6DN6qS53oEpU09iirAixhsqjYaSeQspmbeQL/m/R+2xI+x6700Obfo0pJ0s+3np5z+idPkKJixYTO7kqaIkkGBYGbTQWrFiBR988AGzZs3iqquu4rbbbmP9+vV88MEHnHvuucMxRoEgZukZ+2FvsVG2cxvHt2/h5J5deFxdU9RVWi3Fs+cxYcFiJsxbRIIpJWImbXGnPfIc2FjDx38/REdifc76xlSmnxad2WxCwI08kkJB7pSpJKWlc/jzz8LKArU01rPtzVfZ9uar6BKTKJm3kAnzF1EydwFaQ8IojVowVhm069BiseB0OsnNzcXv9/PII4+wYcMGJk2axL333ktKSspwjTXmEa7D8UVoILuEMSsbW705JEN7YkpqINP1wiUUzJiNWhNeOkLMFhx5fD4/bRYnLQ1O6itb2Px6WVibifMz0BrUKBQSCmX3hwKp+zJFYFno+8Cy6iNW9n1aHazRdtbqUmasyAsfkGDY6Hkzc9a3/gdDspGyndsilARSkjd1BhMXLGbC/EUhM3tFQP3YJqZitAS9I4TW+MDtdHBk8wb+89QfIq7PLJ7IxIWBMiKZJRORJGmERzh2Gah1SJZl7C1uWhqdtDQ6aG1yYGt00trowNbooN3q6q1i0bCTlpdAekESqTkJpOQkkJpjIClNj0IhzpPhorebGb/PR82Rg4EkqTu2YqmpCtkuJTefCfMXIUkSO95+XQTUj2FiSmgplUpqa2vJzMwMWd7U1ERmZia+PrJVj3WE0BqbdCYOLf9iByd276T60AH8Pm/Etpfcfg+ly1aM8AjHBz3de6dfNZncyaagmGppdNLS5KClwUFrkxOvx99nf0q1guQ0HQajhurDzaErJZh3fiEarRKfT0b2yfh9Mn5/x7PP3+114CH7O5Z3tHO0ummqbh/QsSnVCkxZBlI7hFdAgCWQnKFHqew9L5RwS0YXq7mGsh3bKNu5tc/qC6LywtgjpoSWQqHAbDaHCa2amhomTpyIw+HoZcuxjxBaYwdHWysVe7+g/IsdnNy9kzarJWR9UnoGrY0NIcvExXf4sNS28+KvtnSfN9A/EiSmaElO05OcriM5Xd/tocOQpEHqsCId2FjDx/84hOwHSQFnff3UY7TarE7+9tNN9Mi5yZmrS3G0urHU2rGa27Ga7fh6EYUKpYQpy0BKdqgAM2UaOLzVPGxxZQJw2ds5sXsXez/6Nyf3fhG2PiktncmLl1M0Zx4F02aJNCxxTkwIrT/+8Y8A/PCHP+TXv/41iYmJwXU+n49PP/2UEydOsGvXrqgOMJ4QQiu+CMnInpJCXdkxTnyxk/LdOzAfPYIsd/35qbRaCmfMpnjOfIrnLiAlOzdiILtwJ0QHr9tHbZmN6kNWqg5bqTvRElFkqXVKUrIMQfGUlKbHmK4nKV1HUqoOpWrgWcLbrE5s9Q6MUbQODUTA+f0yrU2OgPCqbcdS2x54NtvxunrxEEiEfR6SAr5133Jh2YoyEbPU90CpUpE3dTpFs+dTNHsemUUlIkN9nBETQqukpASAkydPkp+fj7LbdFiNRkNxcTG/+tWvWLJkSVQHGE8IoRU/9JeRHQKlPYrnLqBkzgLypk5HpdGE9SMC2aODz+en/kQr1YctVB22Yj7egs/bt+tPkuBb98e+sBiqgJP9Mm3Nri7hVduOtdaOpbYdtyOy67pkTjrTT88lrzQFtUakLIgWkQLqE1PTOLl7Fyf27KSloT6kvcFoomjWXIrnBIRXgilFBNPHODEhtDo5++yzee2118b17MLeEEIrPujtDlWj0wctVsVz5ouL4TAi+2Uaq9qoOmyl+rCVmqPNeHpYbxKMGvKnppI/NYW80hQqD1qi7t6LR2RZpqGilVce3N6rK1WpUpA7xUTRjDQKZ6RiyjKISRmnSG83VbIsY62t4eSenZzYs4vKfXtC0roAJKam0WZpAhDB9DFKTAmtnvh8Pvbu3UtRUdG4F19CaMUHFfv28Mqvfxq2XGRkjy7dA7UTTFqa6+xUHQoIq6ojVlztoVYZXYKavFIT+aUBYRVJHAyHey9eCXFLSjDzzDz8fqjY10SrJfSPPjldR+GMNIpmpAWsXVph7RoufF4PNUcOcXLPLk7s3tlrFYg5513IhPmLyZs6XeTuigFiSmjdfvvtzJo1ixtuuAGfz8cZZ5zB559/jsFg4O233+ass86K6gDjCSG04oNIFi0RyB5dus8QBNDolbgdoRYrtVZJ7mRT0GKVnpcYDE4XDIxIwlOWZaxmOxX7mzi5r4maY834vV3nukIlkTvJRNHMNApnpJGS3SVoxSzG6HN06ybe/O39va6XJAUZRSXkT59J/rQZ5E+biT5J/H+MNDEltPLy8li3bh0LFy7kjTfe4JZbbuG///0vf/vb3/jvf//Lxo0bozrAeEIIrfhBBLIPD26HlwMba9j4r/C7eIVSImeSkfzSFPKnppJRlNRn6gJBdHA7vVQfaaZiXxMn9zfR2hRq7UpKC1i7FArY90m1mMUYZSKGKkgSpUtXUFd+lGZzeH3g9IIi8qbNpGD6TPKnzRQxXiNATAktnU7HsWPHyM/P56abbsJgMPDoo49SXl7OnDlzaGlpieoA4wkhtOILEcgeHbweHxX7LBzZZubEnqZeg9gv/cEcCmekjfDoBN2RZZnmOjsn9zVRsb+J6qOh1q7uxMtkg3igrxu7NksTVQf3UXVwP1UH99FUVRG2vcFowm5rBgIxXl+68fvMPnflSB7CmCemhFZRURF/+ctfOPfccykpKeHJJ5/kkksuYf/+/Zx++ulYrZErpo8HhNASjBf8fpnqI1aObq3j+K6GkFlwyek6WhpDrSYi9UBs4nH5qD5sZf+GGk7saQxbnzfFxLzziyiYloJCWB9PiYHe2NlbbFQf3E/lwb1UHdxPw4nw8lAA+dNnkT9tJjmTp5A9cQqGZONwDX1cMJz/34MuKn399ddz1VVXkZOTE5g9cd55AGzZsoWpU6dGdXACgSB2kGWZ+hOtHN1Wx9Htddhb3MF1CSYtkxdlMWVRFukFiRzcVBs2Q1CIrNhDrVVSPDud9IJETu5tDCtLVH2kmeojzeiTNUxZnMXUpdmk5yeNzmDjnKS09AFZzg3JRiYvWc7kJcsBOLZ9C+se/nVYu6oDe6k6sDf43piVTc6kUrInTiFn8hQyiieE1FYVrsfRY0izDv/1r39RWVnJV7/6VfLz8wF4/vnnMZlMXHbZZVEfZLwgLFqCsYjV3M6RrXUc2VZHS0NX5QdtgopJ8zOZsjiLnImmsEB2MUMwvuiZXHXBBUW4HD6ObqvD2eYJtkvLS6B0SQ5TFmeRYAovki6ILhEn70gSy676Os21NdQeO4K1R41GCBTIzigqIXviFHweD/s/+VDUauyDmHIdCnpnvAotcac0Nug+4wzg6LZ6jmwz01jZFmyj0igomZPBlEVZFExPHVTmdUHsE0kc+3x+KvZbOLy5lvI9jcGYLkmCgmmplC7NpmRuhkiQOoz0N3nH2d6G+fhRzMeOUHvsMOZjR4IxXZGRWLH6OgpnzCKtsCjE8jVeGXWh9cc//pGbbroJnU4XLMXTG7feemvUBhdvjEeh1T3DurhTil96pmPojkIhUTAjlSmLsiienY5GN+iIA8EYwdnu4diOeg5vNmMuswWXq7VKJs7PoHRpDnmTu6ybIl1E9BjM5B1ZlmltbKD22GGObN7Ikc0bem0rSQpScnLJKCoJPIoDz4kpaSF57Mb6DfWoC62SkhK2b99OWlpasBRPxM4kibKyyIF744HxJLTstmYqD+zl7T88BCIfVdwiyzLlexp576m9Yesyi5OYtjyXifMz0CeGlx8SjG+a6+0c2WLm8BZzyOSHxFQtpYuzUWkUbH2rXKSLGGV6Sy+RWzoNa001jhZbxO30SclB8eVqb2ffJx/CGL6hHnWhJRgY8Sq0+rpTsduaaaysoKm6gqbKCpqqAg9Ha+9pPK76+f0UzJg93MMWnALOdg+Ht5g5sKEGS017xDaX/3AeeaXju9qDoH9kWab2uI3Dm80c21Hfax1GMfN09OjN9SjLMu1WCw0ny6k/WU5Dx8NaU40s911rdO75F5E9qZTUvHxSc/PjPru9EFpxQjwKrZDiypLE1NPOQKtP6F9QSRJJaem0NjaELhYWrZil8w/xwGc1HNtZj88TuJAqVRK+HrmUxJ+iYCh4PT5O7Gniiw8rqCsPv3acuXoKM8/IH4WRCQbjevS4XTRVVtBwspyynds4tu3zfvtPSEklLS+flNwCUnPzSc3LJy2vgMTUgAsy1l2Poy607rjjjgF3+Lvf/e6UBhTPjKbQ6u8klmUZR4sNq7mWZnMNzeYa6k+UU7Zza98dSxLGzCzS8gtJzy8kreORmpePWqsTGdbjAGdbwHq1/7NqrGZ7cHlaXgIzVuQxZXEWx3c1iILNgqjRZnXyt59uihjzl1GYxPTTc5myOEvE+8UBvbkep51+Fm2WJiw1VbRbLb1ur9bq0Ccl09JYH9x2yWVfZf7Fl6FPSo6ZYuejLrTOPvvskPc7duzA5/NRWloKwJEjR1AqlSxYsID169dHdYDxxGgJrZ4B6QsvvZLUvAKazbVYO0RVs7kWt8Pef2fA1NPOoGTuwhBB1Rciw3rsIcsyNUeb2f9ZDWW7GoLZ2lUaBZMXZjF9RS5ZxaEXOZGOQRBNeha9zixOprGyLXguqrVKJi/OYuaKPDIKRW6uWKa/G2qXvR1LTRWW6ios1ZXB1811tfh9vl771ej1GDOzMWZmY8rOCTxnZmHMziE5PQOlSh1sO9wWsVEXWt353e9+x8cff8zzzz9PSkogfsNqtXL99dezYsUK7rzzzqgOMJ4YDaEV8W7jFMksnoBSExr8LNHjrqPbH3T4DYnUo2nPbXu+7b3vgfQf3r7vO6Tw8fTYPnyDQfXff/uBH99APjuvx4fb5UOjVSIpFLQ2OWlpdOB2dl3gtAYVxgw9Sal6FMpunYz6Z9fP8Q32swsbz2D6P7XzcOyf54P77EHC7fTiavegTVCj1avwuv00VrXSUNEWkpsrwaQloyiZtLyEQP3LQX92p/qb7Nl89D+7wbQfzutz5wJnWxv2lmYSjCZ0iUm9tu8+Fr/PR+X+PRza9GnPTvtFkhQkpadjysrG5/VSfeggMHzB+DEltPLy8nj//feZMWNGyPJ9+/Zx/vnnU1NTE9UBxhOjIbQq9u3hlV//dET2JRAIBALBaDMcscAxVYKnpaWFurq6MKFVX19Pa2tr1AYmGBgpOblIkhSWNXjld25H3+1kCdfTXe8drS20WSwkpKSiT0ru3KBbyx7b9uwr7G3f7cOlfc/1/Wwftrkc+XV/bSPtK6x5f2Pp57Po43MfVNvwXdNmcbDj3yfpSYJJQ+5kExmFSajUXQlF+/vc+/6eovgdRdo+rHk/51xf7Qf4HTkP7Kd1/X/pzD+QePbZ6KZPDxv9UD8356GDtH3yaVf/Z5yBbtrUPo9t+D+3fo5tEOf3YM/X3sbidfux1rVjrW3H7eiywuqT1aRkGzBmGPB5/bjsHjR6FWqtMvJY+/iewj+X3ttGbj9y53d/fUf1GhSh+8F+FgNp29rUSGNl17UqraAoIJIG8Lk521qpKzsW2s7vp9lcEzehKoMWWldccQXXX389v/3tb1m6dCkAmzdv5u677+bKK6+M+gAFfZOUls55N/0gzH8+46xzR3togmFAlmVqjjSz/7Nqju2sR6ULT79w4XdFWobekL1eXMfLaN+4kfp/vkp295UvvoqmeAeSWh3wtygUgWepwx0SfC8F3C3dHwop0KZjO7/bjXPXrtCdv7KOCe/ehLa4eASPOH6Q/TJVR6zs/7SG8i8a8PpkGqrBUq8IzpAV+bjil6HG8kYsQaRQYMqOn3Ng0K5Du93OXXfdxTPPPIPHE/Cxq1QqbrjhBh5++GESEuI7l8apMNqzDkVA+tjF2ebh0OZa9n9WQ3Nd75MaRFqGLmS/H/fJkzj37cO5bx+OvftwHjyI7HD0v/FwoVCgnVqKfvZs9LPnoJ8zG01JCZJClDLqjr3FzcFNNez7pJo2qyt0pQSrf7mUlCzD6AxOMOKMxOz2mIrR6qS9vZ3jx48jyzKTJk0a1wKrk3jMoyWIXWRZpvZYM/s+reH4rvpgjTm1VsmUxVnMWJFHQ2WrSMtA4LPy1tQExNT+DlG1fz/+COEMCoMBzeTJOPfsCXVdKBTkPvIwKpMpcPcsA7I/0EaWA8v8MtDx3t+5jo5nf7Cdz2Kl7v77+3V5AiiSktDPmomuU3zNnoUqXdwsAVQdtLDuD1+ELVdpFExblsP0Fbmk54sZi+OB4TYmxKTQEoQjhJYgGjjbPRzeHJ73KqMwiRkrcpm8KDT/0FhOy+Axm3GfOImmuAh1dpejz1Nfj3Pffpz79uLYtw/nvv34LOG5fCStFt20aehmzgyImZkzgxak5n/9i9qf/wL8flAoyPnVGkxf+UrUxt6z/+w1a0g8bTmOPXtw7N6DY88enPv3IzudYduqc3PRzemyeummT0ehC/1ue/tsxhJ95ePqJKskmemn5zJ5YVZX/JZAMEiE0IoThNASDJbOorvGDB1tFhf7e2RtV2mVTFmYyYwz8sgsGl/nVIhQkSQSzz0HZHDu24e3ri58A5UK3ZQp6GbNQjdzBvpZs9BOnBiIueoFj9mM+2QFmqLCYREr/fUvezy4jh4NiK89e3Hs2Y37eFm4Jazz2DrEl8dcS+Njjw+bSIwlQvJxKeDMr5WSnK5n/2eBWC6/v8PSq1MyZXE2M07PFXm5BINGCK04QQgtwWA4sLGGj/9+KOLdenpBYiBr+6IsNPrxlz3bVVZG2cWX9O56UyjQTpyAbmY3UVVaikKrHdmBDgO+1tZATFmH1cuxZw++xsa+N1IomLT+ozFt2YpktbW3uDn0eS0HNtRga+iKvcssCmSf72n9FQh6QwitOEEILcFA8Ptljm4z8+GzB8PWTVqQydwvFZJZnBQzpSlGClmWcezaRfNrr9Hy1tvILldYG9PXvobx4ovQTZuGYpzEhXbFn+3FsXsPbRs34j5yJKxdwhkrSP/Od9HPmzv+zh2/TPURKwc21HB8VwN+X1c84+TFWcw4PRdDsobmegemMehiF5w6QmjFCUJoCfqiuc7Ooc9rObzFHD6TqoPLfzj+UjN46uuxrVuH7dXXcJ840XvDMW61GSges5lj55wbcBtGQDNxIqavfAXjZatQpaaO8OhGH0erm0ObzRzYEHmGrkgRIYjEqAutN998c8Adrlq16pQGFM8IoSXoicvh5dj2Og59bsZcZgsu1+iUISVyYHylZpDdblo/+QTbq6/R9tln0FEPTTIYSF65EtOXr8RVfgLzL4YvWD2e6Rlon3rDt/E1NNLy7393Bder1SSdcw6mr3yZhOXLkZTjK1C8s97nFx9VcmJ3uOv1rG9MZerSbJQqkVpDEANCSzHAHC+SJOHro4DkWEcILQEEXIPVh6wc/LyWsi8aQpItFs5MY+rSHEpmp3N4q3ncpWZwHjmC7bXXsb35ZsgsQf38+Zi+fCVJKy9AmdjlEhzuYPV4JtJn42ttpeWdd2n+179w7tsXbKvKycF05ZWYrrwCdV7eaA15VKg6bGXd73dFXKdLUDN5YSalS3PGpbte0MWoCy3BwBBCa3zROWOwM+ajuc7Ooc21HN4c6hpMyUlg6rJsSpdkk2DUhvUxVlMzdBL883/ttUDuqg6UGemYLr8c4xVXop1QMoojHJs4Dx2i+V+vYnvzTfwtLYGFkkTC8uWYvvoVEs85B0WP4vFjkd5SROgS1SGFrU1ZBkqXZDNlSRbJafoRHqVgtBFCK04QQmv80HPGoDFDHzLrSWtQMXlRFlOX5ZBZND7ulLvndVJlZmLfupXmV1+j9f33uwLbVSqSzj4b45VXkLhiBZJKzAgbbvxOJ60ffEjzq69i37w5uFxpMmG87DJMX/ky2smTR3GEw0/PFBFnfT3gNqw6ZOXwFjNluxrwerpi3vKmmChdms3EeZnjctbveCTmhFZ7ezuffPIJFRUVuN3ukHW33npr1AYXbwihNT6w1tn55y83R6j+CoUz0pi2PIfi2Wmo1OMnJqZnziul0YivuTm4Xjt5EsYvfxnjpZeiSksbvYGOc9wVFTS/9hq2117HW18fXK6fMwfTV79C8oUX4mttHZOJUPuyHrudXo7vbODwllqqDzcHl6vUCkrmZlC6NJuULD0tTS4xa3GMElNCa9euXVx00UXY7Xba29tJTU2lsbERg8FAZmYmZWVlUR1gPCGE1tjF4/Zxcm8Tx3fWU7a7IVgOpzsrb5zJpAWZozC60cVjNnPs7HPCcl5JBgPGVZdi+vKX0c2cOS6sevGC7PXS9tlnNL/6Km3//Tg4GQGNBjpvnsfpBISWJgdHttZxeLM5cl1RCc5aXcqMFeMr1m2sE1NC66yzzmLKlCk89dRTmEwmdu/ejVqt5hvf+Aa33XYbV155ZVQHGE8IoTW26BRXx3bUc3JfI1535On0ML5mDHanfctW6u6/H9fhw2HrCv7yZxJXrBiFUQkGg7ehAdu6dVhffAlPdXXoynGcUkOWZepPtrLv4yoObTaHrZ+6NJvSZTnkTjKiUIqZi/FOTAktk8nEli1bKC0txWQy8fnnnzNt2jS2bNnCtddey6FDh6I6wHhCCK34omcwO4DH5ePkvsjiKilNx6T5mUxamElDZSuf/PPwuJox2Iksy9i3bKHx8Sewb98eudE4/oOOV9o2b6byuuvDluvmzCb7Zz9DP3v2KIxq9Olr1iIEgupL5qQzYW4GBVNTUaqF6IpHhvP/e9BRfmq1OugCyMrKoqKigmnTpmE0GqmoqIjq4ASC4aJ7MLskwfTTc3G2e8PEVXK6jonzM5m0IJOMwq6g9syiZIpmpI35GYPdkWUZ++bNNDz+BI4dOwCQ1GpMX/0K6tw86n/3u5CcV0JkxRfa4mJQKMISoTp37+HEVVeTeOaZpH//FvSzZo3OAEcJU6YeSQr1jEsSTJifSdUhC842Dwc31nJwYy0anZKiWelMnJdB4Yw0UeRaAAzBonX++edz3XXXsXr1ar7zne+wa9cubr31Vl544QWsVitbtmwZrrHGPMKiFR/0Nt27k97E1XhFlmXaN22i8YkncezcCYCk0WD66ldJu/F/goJK5LyKf3omQs249VbcJ05ge/PNoAAbj4Ir0qzF6afl4vf5qT7aTNmuBsq+aMBu65ocplIrKJyRxoR5GRTPSsPj8okSQDFMTLkOt2/fTmtrK2effTYNDQ1ce+21bNiwgUmTJvHss88yZ86cqA4wnhBCKz7ozRUweVGgzqAQVwFkWaZ94yYaH38cxxdfAB0C66qrAgIrK2t0BxgjmNvNVLRUUJhcSHZC/AvMSILZfeIEjU/9Cdtbb41bwdVfzjvZL1N3oiUwYeaLBloancF13S1iogRQbBJTQkvQO0JoxQeRLFrjNZg9ErIs075hA42PP4Fj924AJK0W09VXkXbD/6DOGn8zK3vjtaOvsebzNfhlPwpJwS+W/YIrJ4/dCUERBddZZ5F+yy3oZ80c5dHFDrIs01jVRtmuBo5ur8NW7whrM3VZNhPnZ5I3JUW4GGMAIbTiBCG04ofeXAHjGVmWaf/sMxqeeALn7kAGd0mrJeWaq0m94QbUmfElsKJhaZJlGbvXjsVhocnZhNVpxeK0YHFaqGyt5PVjr4e0l5B46eKXmJ4+PRqHELMIwTVw+gumVyglciYZKZiWSuH0NNLzE5EUwqI+0sSU0Kqrq+Ouu+7io48+or6+np6bi1qHQmjFC+Oh/E1fdGZyVxcV4jp8mMYnnsS5dy8Akk5HyjXXkHbDt1FlZIzySAdPX5Ymj88TFEvdH03OJiyO8OUun6ufvYWTk5DD7IzZzE6fzeyM2UxLm4ZWqe1/wzijP8HVvVrAeI3bixgTKsGURVnUHrfR2uQMaa9PUpM/NZXC6akUTEslwaQN9iNivIaPmBJaF154IRUVFXz/+98nJycnLJblsssui+oA4wkhtATxQkjQczcknY6Ur30tILDS00dpdKeGud3Myn+txE/oseUn5mNz22h1tw66T51SR5o+jVRdavChVWp5+fDLyJFKBPRApVAxNWUqszNmMytjFnPS55CflD9mYgFd5eU0/elP2N56O3hOaUtLcR05EghOGqfJTzvpzYIuyzK2egeVBy1UHLBQfdiKxxVqrEjNTcBg1FB1yAqyiPEaLmJKaCUlJfHZZ58xd+7cqA5kLCCEliAe6C2Tu+maq8n4wQ/itkSOX/az3byd/937v3xe+3mfbZWSMkQ0pepTSdGmhImpzodBbYjYTyTL2crilexr3Mfexr3sbtjNnoY9WJyWsG1TtCnMypgVtHrNTJ9JkiYpKp/FaBEUXG++FXZ+jffcagOxoPt8furKWqg40ETlAQv1Fa0RS31BIDt98Zz0sEL1gqERU0Jr+vTp/OMf/2DevHlRHchYQAgtQazjPHyY2nt/jnPPnrB1hc8/T8KSxaMwqlOjsqWSN8ve5M1jb1LTXhOxjYTEI2c+wqSUSaTp0kjSJKGQopNY0txuprK1koKkgoixYLIsU9New56GPYFH4x4ONh3E4/eEjXGCcULA5djxmGiciFKhDO4nXmY3Nr/xBrU//knY8tzfPoLx4otHYUTxibPNw+71lWx/90SvbZLTdWRPNJIz0UTORCOpOQkixmsIxJTQev/99/ntb3/L008/TXFxcVQHE+8IoSWIVdxV1TQ+9sfIlgaIO2tDu6ed90+8z7rj69hRtyO4PEmdxAUlF2DUGnlm3zMxOxvQ7XNzyHIoKLz2NOyhuq06rJ1BZWBm+ky0Si0bqjcgI8fk8fTEYzZz7Jxzw1zTqFSkfv3rpH3nZlQpKaMzuDijt7x/pmxDoBZjj+UavYrskuQO8WUkq8SIWqsUMV79EFNCKyUlBbvdjtfrxWAwoFarQ9ZbLOEm8vGCEFqCWMNrsdD4pz/R/OJLyJ6ABSXpwgvQlU6l4Y9/DMnkHuvxM37Zz466Hbxx7A0+OPkBDm9gyryExLLcZVw+6XLOLjgbnSrwJ9KfpSnWaHI0sbdxb9DytbdxL3ZvhKLGgEJS8J8v/yemj6tn8lNNcTHusjIAFImJpP3PDaR+61soDJHdsoIueovxcju8mMtt1B63YT5uw1zegrdHjJekkEgwaWizdEzqkODM1aXMFEWxQ4gpofX888/3uf7aa689pQHFM0JoCWIFv92O5fnnafrfv+JvbwfAsGwpmXfcGZx+Hy+Z3Ktaq3jz+Ju8efzNEKtPcXIxl026jEsmXBLTgmOo+Pw+ymxlrDu2jucPhF93L594OXcvvptkTexea7qfY6qsLNo3bqL+t7/FdfAgAMqMdDJuuQXTl7+M1OOmXRDKQGK8/D4/TdXtHcKrmdrjNtqskWfNpuUlkFWcTHpBEhmFSaTlJfaaz2s8WMNiSmhFkwceeIDXXnuNQ4cOodfrWb58OWvXrqW0tDTYRpZl1qxZw5///GesVitLlizhiSeeYMaMGcE2LpeLu+66ixdffBGHw8G5557Lk08+SX5+frCN1Wrl1ltv5c033wRg1apVPPbYY5hMpmCbiooKbrnlFtavX49er2f16tU88sgjaDSaAR3PeBJa4+GHF4/IHg/N//oXDU88ia+xEQDd9Olk3HkHiaedNsqjGzh2j533T77PumPr2F7XVbg6UZ3IBSUXcNnEy5iTMWfMzNrrC3O7mZWvrsQv+8PWJaoT+drUr/HN6d8kRRcfrjjZ76fl3fdoePRRPFVVAGiKisi4/TaSLrhgXHynI8nRbXW8/9f9/baTJDBlGUgvSCK9IJGMgiQyCpIo290QUhd2rM54jDmh5fP5eOONNzh48CCSJDF9+nRWrVqFUjm47LYXXHAB11xzDYsWLcLr9fKzn/2MvXv3cuDAARISEgBYu3Yt9913H8899xxTpkzhN7/5DZ9++imHDx8mKSkwQ+e73/0ub731Fs899xxpaWnceeedWCwWduzYERzThRdeSFVVFX/+858BuOmmmyguLuatt94KHtPcuXPJyMjgt7/9LU1NTVx77bVceeWVPPbYYwM6nvEgtLxuHzv+czIQnCkDEpw9Rn948YTs99P6739T/4c/4DkZKO6uLiwk47ZbSb7wQiRFdAK/h5NO1+C6Y+t4/+T7Ia7BpTlLuWzSZZxTeA56lX6URzry9JzdeMWkK9jdsJtjzccA0Kv0XF16NdfOuJZ0fXyk5ZDdbqz/9wqNTz6JryPkRDdzJpl33UnC0qWjPLqxQ8RKGBKccU0pbVYnDZVtNFa2Ym9x995JNyQJvvbLpaRkjS2Xb0wJrWPHjnHRRRdRXV1NaWkpsixz5MgRCgoKeOedd5g4ceKQB9PQ0EBmZiaffPIJZ5xxBrIsk5uby+23384999wDBKxXWVlZrF27lptvvhmbzUZGRgYvvPACV199NQA1NTUUFBTw7rvvsnLlSg4ePMj06dPZvHkzS5YsAWDz5s0sW7aMQ4cOUVpaynvvvccll1xCZWUlubkB0fDSSy9x3XXXUV9fH/GDd7lcuFxdZtmWlhYKCgriRmj1ZZVyObxYa9uxmtux1Nqxmtux1raH1O/qRJSvGXm6J4J0l5dT/8hvce4P3LUq09JI/953SfnqV5EGaI0dDTpn0akVajbXbmbd8XUhrsGi5CIum3gZl068dEy6BgdLz5gzv+znvxX/5ek9T3PQEnDFaZVavjz5y1w/8/q4+cx8be1Ynn0Wy7PP4rcHYtISTjuNzDvvQDd9bGfYHykGUgmj3eaisbKNhspWGitbaahso6UhvHRQJ0mpOlKyDaRkJ5CSYwi+1ieFX3PiwQMSU0LroosuQpZl/vGPf5CamgpAU1MT3/jGN1AoFLzzzjtDHsyxY8eYPHkye/fuZebMmZSVlTFx4kR27twZkk7isssuw2Qy8fzzz7N+/XrOPfdcLBYLKd1mscyZM4fLL7+cNWvW8Mwzz3DHHXfQ3Nwcsj+TycTvf/97rr/+en7+85+zbt06dnfUdoOAuzE1NZX169dz9tlnh433l7/8JWvWrAlbHg9C68DGmqA5uDNLsTZBHRBXte202wZ2d9PJ5T+cR15pfLgu4p3eko0qDAZSb/g2adddh6LDIhyrvHrkVdZ8viYs2WeCOoELii/g8kmXjxvX4KkiyzKfVX/G03ueZk9DIG2HSqHi8kmX8+2Z36YgqWCURzgwvE1NND71J6wvvwwdEzeSL76YjNtuRdJoxn2G+VNlKJUwLLXtvPirLb3m8oqELkHdIboMmLITaLU42ftxVcwnWx1OoaUa7AaffPIJmzdvDoosgLS0NB588EFOO4UYEFmWueOOOzj99NOZOTMQrGs2mwHIysoKaZuVlcXJkyeDbTQaTYjI6mzTub3ZbCYzQp22zMzMkDY995OSkoJGowm26clPfvIT7rjjjuD7TovWSDGQuwS/X6bN4qSl0UFLoxNbgwNLTRsn9jZ1NZLhyNa6sG0TjBpSchJIyUkgteNuRaNX8coD28IKMhszx587ZzTwmM3U3vvzsBQNxiuvJPPOO+Ii2eg28zZ++fkvw5b/ePGPuXLylePSNXgqSJLEGflnsCJvBVvMW3h699Nsr9vOv478i9ePvs7FEy7mf2b9DyXGktEeap+o0tLI/n8/I/Vb36ThD3+k5Z13Ao/33uu6qYiTGbKxSGKKbtDWpNScBM7+xtQwa1jJnHSsnZ4Os73j0U6rxYmz3UPt8cBMyJ7IMvz3hUNUH7aSlpdIcrqe5HQdyel6tAZVxBureLCG9ceghZZWq6W1NbyERVtb24CDxiPx/e9/nz179rBhw4awdT0/fFmW+73T7dkmUvuhtOmOVqtFqx2drLzdLVKSBAsvKiY9Pwlbg4OWRge2RgctDQ5am5z4/QO7HZm0IJPCGakdpuAEtPrIp8dZEX548foDiCfclZWYf7kmYh4s42WXxbzIsnvs/GXvX3hm3zMR109JmSJE1ikgSYFYtqU5S9lRt4O/7PkLG2s28ubxN3nr+FtcUHwB/zP7f5iSMmW0h9onmsJC8n77CGk3fBvz/Q/g2N41GQK/n9qf/4KE008fVcvWoBxBA2w64B4HuO+B99f36tIl2eRNMWFrcGBM7xI7WcXJZBWHWn48bh/N9fagCKs51kxdWUtYn4Eb+9Cbe4VKwpiuJylNjzFdR3KGHluDg32fVse8Naw/Bi20LrnkEm666Sb++te/snhxIIv0li1b+M53vsOqVauGNIgf/OAHvPnmm3z66achMwWzO35IZrOZnJyc4PL6+vqg9Sk7Oxu3243Vag2xatXX17N8+fJgm7q6cItNQ0NDSD9btmwJWW+1WvF4PGGWrtGmzerscvsR+N1te+fEKfdbe9xGXXn4jyJSLTd9sga/T0ahlNj2Tjnb3invbBxVBnw9G+GLz9D6HNqHI/v8yA4HfpcLuBSWXxrWRvmyA1757JQGOJCafYPpr3t7j+zB6XUhyxO5lvsiNJLYtaudL6RP++lrdL7nQR3yQL/nqP8Bh76dyzXM5ip8fl9gxuJmeO+lcv4jnUQpKft3yw74tzfQAQ7ht5J0PfKZ4SmD1v9yP5J0YFh+z4LRwe+Vg9axSMgyfPyPQxROT427G/tBC60//vGPXHvttSxbtiyYrNTr9bJq1Sr+8Ic/DKovWZb5wQ9+wOuvv87HH39MSUmoabukpITs7Gw++OCDYIyW2+3mk08+Ye3atQAsWLAAtVrNBx98wFVXXQVAbW0t+/bt46GHHgJg2bJl2Gw2tm7dGiIObTZbUIwtW7aM++67j9ra2qCoe//999FqtSxYsGCwH9Ow0lzvGOp/dp+0N0fOtyKIBTTQh8XYY/cBvl7Xjz4SGvq+OHp8fiA8hYHg1FCgREG3GeFyQHcMWFiPNr2UShq9xESC0UL2g63eMfaFlslkYt26dRw9epRDhw4hyzLTp09n0qRJg975Lbfcwj//+U/WrVtHUlJSMBbKaDSi1+uRJInbb7+d+++/n8mTJzN58mTuv/9+DAYDq1evDra94YYbuPPOO0lLSyM1NZW77rqLWbNm8aUvfQmAadOmccEFF3DjjTfy9NNPA4H0DpdcckkwZ9f555/P9OnT+eY3v8nDDz+MxWLhrrvu4sYbb4y5wHZTph5JImy67gXfmYUhufc/Y4nQO1h7i5u2ZieJJl3odgOMP452nHLUA5+jfhwDa3iqh+Gz2Wh+9VVa3nkX2eMGGfSzZpHyja+jmxGYheVtbMRTU4s6NwdVenrH6AZ6wANsdgrH4fK5eOXIK7x25DU8she1pOKKyVdwVelVwcztDY4GatvM5CRmk6HPGNqO4vQ7HuJuT/k3crLlJC8depH1lf/tyMvVdRFRSAp+uOCHXFhy4ZAHOODhDeIwWt5+m7qHHqYzTsJ42Sqcu/fgKg9Y0ZVJSaSsXo3py1dGLct8tI8j2r/NwTDwYxngdzzgHXe9bGt20dLgIDlDT6JpYOE2bc0uXlqzZUzEA49qwtLeLhrPPvss1113HdCVsPTpp58OSVjaGTAP4HQ6ufvuu/nnP/8ZkrC0e2C6xWIJS1j6+OOPhyUs/d73vheWsHSgcVgjmUdrINN1BfGFr7mZpmefw/LCC8gd09z18+aRcdutcZNXSJZlPq78mLXb1gZTNSzPXc5PFv+EYmPxqI4tiK0aLMchdSIYR6gMyWjssw8qWip4bNdj/PvEv0OWx2ppn55VDDrzxjU8/kSwrI8yNZW0G28k5WvXoNDFl8VDEJmR/J8b9fQOd9xxB7/+9a9JSEgImWUXid/97ndRG1y8MdIJS4cyXVcQe/ja2rA8/zyWZ5/D39YGgG7GjIDAWrEiblIcVLZU8uC2B/m0KhBnlZ2QzT2L7uHcwnNj5xh2/g3euo3glfvSP8D8bw2tr+6XzuDrCMt2/R3e+WF09hlFttZu5Yb3bwhbfvv827lhVvjyWET2+Wh5+20anngST0UgUa8yI530m27GdPVVKGI4j5xgYIzU/9yoC62zzz6b119/HZPJFDGfVLAzSWL9+vVRHWA8MR4ywwuig8dsxnX4MO07dmB7+f/w2QJTobVTppBx260knnNO7IiTfnB6nTyz7xn+uvevuP1uVAoV1824jhtn3YhBHUPZoxuPweMLCY+MVnS4TSIJJyK0P0UkJdy+d9QtW32V9rl80uXcvSi26yh2R/Z4sK1bR+OTT+GpqQFAlZ1N+ne+g+nKK/BaLCIPl6BPRl1oCQaGEFqCgWD5xz+o+819IX/mmgkTyPj+LYFab3FQLqeTT6s+5YEtD1DVFqhZtzRnKT9d8tPYytlkKYdt/wvbnwVP+2iPJsDCb8M594Ihtf+2w0jP0j6LsxezpXYLMjKZhkzWLF/D6Xmnj+oYB4PsdtP82ms0PvUnvB0zzZUmU+BGRpZFHi5BrwihFScIoSXoC09dPY1/+hPNL74YukKSmPjhB2jyRj92Z6BUtVaxdttaPq78GIBMQyY/WvQjzi86PzYscX4/HF8PW/8MR9+nV6uUpID/WQ9JHVaO4Ni7HUPI8UgDXNZBqxn+dFrAbdgTlQ5mfQUW3Qi5c/s9pOGiZ2mfnXU7uXfjvVS0BlxxV06+krsW3kWSJmnUxjhY/C4XzS//Hw1PPYXfag1dqVAwaf1HwrIlCGHUhdaVV1454A5fe+21UxpQPCOEliASzsNHsDz7LLZ33gmWFulJ4fPPk7Bk8QiPbOB01iXMTsjm3fJ3+d+9/4vL50IlqfjmjG/yndnfiQ03odMGX/wTtv4lEHzeycRzYfFN0GaGt+8A2Rdw4V366PDHS+38G7x1e9c+Z18FdfvAvLerTf5iWHwjTL8MVKOTBLk7Dq+DP+78I/84+A9kZLIMWfxq+a9Ynrd8tIc2KNo+/ZTKm24OW5504YVk/ehu1N3yMwrGN6MutK6//vrga1mWef311zEajSxcuBCAHTt20NzczJVXXsmzzz4b1QHGE0JoCTqRZZn2TZuwPPsc7d2qHehmzcK5b19oDFCM32F3dy91Z3H2Yn625GdMME0YpZF1o+4AbPsL7H65yz2oTYa5X4dF/wPp3dLP2KrBUgapE0Z41mG3fcoyVG4JCMID68DfIcATMmD+tbDwejDm993nCLCjbgf3bryXytZKAL48+cvctfAuEjWJozyygeExmzl2zrlhdUEBUKlIvuACUq+/Dv2MGSM/OEFMMepCqzv33HMPFouFP/3pTyiVgSR4Pp+P733veyQnJ/Pwww9HdYDxhBBaAtntxvbuu1iefQ7X4cOBhQoFSeefT9r116GfMye0KHSMx4yY282s/NdK/D0Sif50yU+5pvSa0XUT+rxw+J2AWDnRLSt+xrSAdWj21aCNA0HQWhewem1/BloDgdxISph6UcAKV7xiGBJ3DRy7x84fdwWsWwA5CTn8cvkvWZ4bH9atnr+3lNVfw3X0GPZulUAMixeTet11JJ51ZlzFSAqiR0wJrYyMDDZs2BBM9NnJ4cOHWb58OU1NTb1sOfYRQmv84mtpwfryy1hf+Dve+noAJIMB05e/TOq130KTH2qd6JkXKFZ58dCL3L/l/rDlz6x8hkXZi0ZhREBbA+x8LhDc3hLI1RUQJhd3CJPTR1WYDBmfBw6/G0E4Tg1Y5eZcA9rRi5PaZt7Gzzf+PDjx4StTvsJdC+8iQZ0wamMaKJF+b479+7E893ygaLXXC4CmuJjU667FeNllKPTxlxhTMHRiSmilpKTw7LPPcvnll4csf+ONN7j++uux9gw8HEcIoTV+8JjNuE+cRNJqaf33ezS/8i/8HUlGVRkZpHzzm6RcfRVKo3GURzo0fH4fz+x7hsd3PR5mzRq1pJZVOwLB7ftfA587sMyQDguuixlXW9SoPxgQXLtf6nKFapICYmvxjZBR2vf2w4TdY+fRnY/y4qHAhI6chBzWLF/DstxlozKeaOAxm7G88ALN//cK/tZWIDBTMWX110hZvRrZ6xWpIcYBMSW07rjjDp577jl++tOfsrQjW/XmzZt58MEH+da3viUSlgqhNeaxvvIvzD//eVixNe3kyaR++9skX3xRXCdKNLeb+clnP2F73XYAZqbN5IDlQDAFwC+W/YIrJw98gswp4XHC/tcDAqtmZ9fyvAUB69WMK2IieHzYcLYExNa2v0Djka7lJWcEZiuWXgTKQVdSO2W2mbdx78Z7g9n/r5pyFXcsvCMurFu94Wtrx/baq1ie/xue6g5LqVIJvo4aojHu5hecGjEltPx+P4888gh/+MMfqK2tBSAnJ4fbbruNO++8Mxi3NR4RQmts46mpwfL3v2N5JnzCR87atRhXXRobqQ1OgQ9PfsgvNv2CFncLBpWBny39GZdOuJQ6e11ICoBhoXuZGtkfiFna+TzYO8IRlBqY+eWAwMiPrULvw44sQ/knASvX4Xe70kUk5wWsefOvg8SMES31Y/fY+d2O3/Hy4ZcByEvM47b5t5GmS6MwuTDmyvgMFNnrpfXDD2l8+mlcBw+FrpQkStatQzdl8ugMTjBsxJTQ6k5LSwuAEBUdCKE19vA7HLR++BG211+j/fPNYVasTmI9PUN/2D12Htr2EK8efRWAGWkzeOiMhyhMLhyZAXQvjQMEclN1fNbJ+bDo24HZeAnpIzOeWKa5EnY8CzueB3tjYJlSA9mzA1a/ES71s6V2C7/Y9IugdQtAgYJfLB9By+cw0L55CxUdNXdDUKlIOu9LmC6/nITTTkNSjbxFURB9YlZoCUIRQmtsIMsyji++wPb6G7S8+26w/iCAbu5cnLt3x1V6hv442HSQH336I060nEBC4tszv80tc29BrVSPzABs1fDozPCkngVLYPkPYMqFo+Iei3m8Ltj/RsCtWr09fP0Ilvopt5Wz6o1VIctitUD1QOkzNUQHyox0jKtWYbr8crSThZUrnhnO/+8BXb3mz5/PRx99REpKCvPmzevTPbJz585e1wkEsURnQHtnkKunrh7bunXYXn8dd3l5sJ06Lw/j5ZdjvOJyNPn5EdMzxKPI8st+/n7g7zy681E8fg+Z+kzuX3E/S3KWjOxA9r4SOXP6OfdCyYqRHUs8odLCnKsDj+3PwNs/DF0v+wK5u0ZAaDXYG8KW+WU/z+57lh8v/nFcutTV2dnk/GpNyG89e80a9DOm0/z6G7S89Ra+hkYsf30Gy1+fQTdzJsYrLif5ootQpaQA4dcYwfhkQELrsssuQ6sNBJz2nG0oEMQjIWJJktBMmoT7+PHg3auk15N8/vkYr7gCw+JFIbl1TF/5Cgmnnx4X6Rl6o9HRyP/b8P/YWLMRgHMKzmHN8jWYdKaRG4TfDxt+B+t/Hb5OUgaSew4ztTYH5Y3tlKQnkGOM4+n8k1cG3IU9BWvD4RERq4XJhSgkRVhS238e+ie17bX8avmvRvbcihK9/dazp08n6+67aPv0U5pff4O2Tz7BuW8fzn37qHtwLUlnn40qMwPrP1+Mi3x5guFlUK5Dn8/Hhg0bmD17Nikdil3QhXAdxgd9uQT0CxZguuJyki64AGViHCS7HAKfVn3KvRvvxeK0oFPquHvR3Xx1yldH1urgbIE3vguH3g68L1weyJQ+gqVx/rH5JP9v3T5kOZB2667zS/nGkiKSdCoUivizwISU+unOsu/Dl9YMu/u1Z4HqlcUr+fDkh3j8HrIMWTy44kEWZi8c1jGMFt6mJlreeYfm19/AdfBg5EZxHmIw1ompGC2dTsfBgwcpKSmJ6kDGAkJoxQe9Bbnmrl2L8bJV4RuMEVw+F7/f8ftghu8pKVN46IyHmGiaOLIDqT8EL38Dmo4GgrgvegQWXDtipXG8Pj/PbiznvncPRVyvkMCoV2MyaDAZ1Jj0alIMGoyGwLPJ0LFO3/29mkStalBidVisaZ2fYUpRwJ244feB5cUr4CvPBmYmDiM9C1R3j/9TSApunn0zN82+CZVi7MbcOQ8dovHJp2h9//2wdQmnLcf01a+ScPrpY/ZGLl6JKaG1aNEiHnzwQc4999yoDmQsIIRWfBDRojWG7zbN7WY+r/mcZ/Y/wwnbCQC+Me0b3L7gdrTKEc5BdWAdvPE9cLcFUhNc9cKIpWrw+2Xe2VvL7z88QllDe9T7VykkTAY1Rn0EQZagCVm+pbyJx9cfwy8HhN0DV87i6kXDMMPzwJsBy6G7DZJy4aq/QcHIZvS3e+zcv+V+1h1fB8D8zPmsPWNt3AbJD4R+A+nVahIWLybx7LNJOvss1HkjVHNT0CsxJbTef/997rnnHn7961+zYMECEhJCE9SNZ4EhhFb8EE/1Bk+FV4+8yprP1yB3pEpIUCXw8JkPsyJ/hIPMfV5Y/yvY+IfA+xGysEBgFulHB+v57QdHOFjbkZJGp6LV6aX7xU8pSfz3rjPRaZQ02z002z1Y7W5sHc/NDg/NdndweWebZocbp6f3mWkDQSlJbPjx2cMTJ9ZwBF7+eiDhqUINF66Fhd8e8TJFb5e9zW82/4Z2TzvJmmR+fdqvOafwnBEdw0jS8xqTev314PfTtn497pMnQ9pqp04l6ZyzSTz7HHQzpiMpFCKQfoSJKaGl6BYU3N1MLssykiTh8/kibTYuEEIrvoiXeoND5ZDlEF9966shy0Zlyn17E/zr+kDCTQikbDj3lyOSsmHjsUYeef8wuyqaAUjSqvifFRP49unFvLu3lp++tg+fLKOUJO6/cuaQrUpOj6+HAHNj7RBh3d/b7B6qm+1UNzvD+vjn/yxh+aRhyhPmag1YEg++GXg/9xtw8SOgHtkJAJUtldz96d3sb9oPwDWl13DXortG3rI6QvR2jXGVldP23//S+t/1OHbuCrF8qTIyUBcV4dixI5BGZgzfCMYSMSW0Pvnkkz7Xn3nmmac0oHhGCC1BrLCpehN3f3o3Le6WsHUjWhC6eif837fAVgnqBLjscZg5/Eksd5y08Mh/jvB5WSCrvE6t4LrlJdx8xgRSErrKI9XaHJxotFOcbhixWYe1NgenPbgef48r75x8I49eM4+S9GEqYyPLAYviR2sCsxNz5gRctylFw7O/XvD4PPxx1x95bv9zQCBW8OEzHmaCafhnmcYiXquV9k8/pXX9f2n/7LNgzdQQJInchx8i6ZxzUBgMIz/IcUBMCS1B7wihJRht3D43j+58lBcOvBBx/YhatHa+AO/cCT5XoCTMNf+AzGnDust91TZ++/5h/ns4kNdJo1Swekkh3zt7IplJumHd92B4eVtF0JomSaBWSLh9Mjq1grtXTuW65cUoh2vmY9nH8K9vB0ob6VPgy3+FSSMfc7uhegM/2/AzLE4LepWeHy/+MVdMuiIuc25FC7/bjeVvL9DwyCORG6hU6GfMwLBoIfqFCzEsWIAyKWlkBzlGiTmh9dlnn/H0009TVlbGK6+8Ql5eHi+88AIlJSWcfvrpUR1gPCGElmA0Od58nB99+iOOWAPFh68pvYZJpkncv/X+kS0I7XXBe/cEysRAoPDxFX8CnXHYdnmsvpXffXCEd/eaAVAqJL66IJ8fnDuZPFNs5sfqbk3z+mR+/NoeNh4LWOAWFqXw8FfnDJ91q7kyYGms2QlIcO69cPodIx631WBv4Kcbfsrm2s0AXFB8Ad+Z8x2aHE1xXS/xVIgYSC9JKDMz8NXVhzaWJLTTpmJYuBDDokUYFi5ElZIi4ruGQEwJrVdffZVvfvObfP3rX+eFF17gwIEDTJgwgSeffJK3336bd999N6oDjCeE0BKMBrIs8/Lhl3lk+yO4fC5StCn8+rRfc2ZBwI3fc8r9sGKrDvyBV28HJDjnZ3D6ndAttjOaVDTZefSjI7yxqxp/Rz6sVXNyuf1LU4ZPpAwTsizzz60V3P/OQdrdPrQqBXevLOX600qGx7rlccJ7dwfybwFMvQQufwp0I3vt6swg//iux/HK3uDyEbsxiEF6m6zjqa7Gvn077du24di2PSyoHggIsvqOTP2SRPaaX5Jy1VUjfATxR0wJrXnz5vHDH/6Qb33rWyQlJbF7924mTJjAF198wQUXXIDZbI7qAOMJIbQEI43FaeEXG3/Bx1UfA3Ba7mn85vTfkK4fheLLJzbAK9dBewPoTAGX1OQvDcuuzDYnf1x/lP/bVom3I9hp5YwsfnjeFKZmx/dvr8pq58ev7mXDsUDB6IVFKTz0ldlMyBimvEs7nod37wKfG9ImwdX/gMypw7OvPlhfsZ7b/ntbyLJ4r5d4Kgxkso6nrh7Hju3Yt2/Hvm0brqPHIrbTzZiBfu5cdNOno5sxHe3EiUjqEaplGifElNAyGAwcOHCA4uLiEKFVVlbG9OnTcTrDZ9OMF4TQEowkG6s38rMNP6PJ2YRaoeaOBXewetpqFNLwWI96RZZh85Pw/r2BrORZs+DqFyA1+kmNG9tcPPXxcV7YfBK3N+BaOWNKBnedP4XZ+aao72+0kGWZF7dWct87B0bGulW9A17+FrRUBSYtXP4EzLgi+vvpg621W7nh/RvClv9y2S/58pQvj+hY4pWWDz+k+vs/6LedpFajnTIlKLx006ejnTIFhU43bt2Oo15Uujs5OTkcO3aM4uLikOUbNmxgwoTxOWtkLDNef3SxjMvn4tEdj/L3g38HYKJxImvPWEtpaunID8bdDm/+APa9Gng/+2q45FHQRHdmlM3h4S+flvHMxnLs7kAKmcXFqdx5/hSWTEiL6r5iAUmSWL2kkDOmpPOT1/by2dFGfvPOQd7bZ+bh4bBu5S2Amz/pSMPxacAyWb1jxNJwQO/1En+9+ddYXVaum3HdmM4oHw30M2cG3PQ9kjFn/eQneGpqcB44gPPAAfytrTj378e5fz+80tFOqUSVno63ri7wXpLIuOMO0m74dkitV8HgGbRF66GHHuL555/nmWee4bzzzuPdd9/l5MmT/PCHP+TnP/853//+94drrDHPWLNojZeknvHEMesx7vnsnmDA+9emfo07FtyBTjXCM+ps1XDiM/j0kUApHYUKVj4Ai2+MakB1u8vLc5tO8PQnx2lxBuJ3ZucbufP8Us6YnD4uZqjJssxL2yq5752DtLm8w2vdipRYduX94GwOzBwdxtJIEF4vcbJpMoethwGYnTGb+067j2Jj8bCOId7p77otyzKeqiqc+wOiy3nwIM79+/FZLJE71GjQTpiAdkIJmpIJaCaUoJ0wAU1xMQp96ESTeL4xjynXIcDPfvYzfv/73wfdhFqtlrvuuotf//rXUR1cvBGrQmugJ78sy/iam/FUVODYs5e6++8PuIU6GcNlamIZc7uZky0n+aL+C/6y9y+4fC5Sdan8+rRfc0b+GSM/oJ1/gzdvhc686tpkWP1/ULTslLvurP+Xa9Tz4cE6nvr4OE3tbgBKs5K44/wpnD89a1wIrJ5UNzv48at7+OxoIHZrQUfs1sThiN3qXiqpE0kBl/5h2It9d5+8kWXIYt3xdazdupY2Txs6pY7bF9zO16Z+beRd5HHEYJMxy7JMy7//Q80Pfzio/ahzc9FMCIgvX3MzLW+9HbdJVmNOaAHY7XYOHDiA3+9n+vTpJIoCmTEptHre3WT/4hcknLYcT2Ul7opKPJUVuCsqcVdV4qmoxN/W1md/hc8/T8KSxSM0esFrR19jzaY1+OlyBZyWdxq/OW2UAt5t1fD7GdC9eI2kgNv3nbK14+VtFfzktb1hiTyL0wz88LwpXDI7d/hyS8UJkaxbd51fyrdPHwbrVvmn8PylocskJdy+d9gtWz2pbavl3k33sqV2CwBLspfwq9N+RW5i7oiOYyzTWw3Yohf+hq+1FXdZOe7yMlxl5biPH8fX3Nxvn/pFi9BOnIA6Lw9NXh7q/HzUeXkoU1Mj3iyNpkUsJoWWIJzh+qIGc/L5XS68ZjMecx2uw4eoe+DBUKvUAFBlZaHKysS5Z2/oCmHRGlFq22pZ+erKYJ1CAAmJ/3z5P+Qk5ozOoD5+AD5+MHz5tW9DydDrJ/aWLf0nF0zl2ytKUCuF9aI7Pa1b8wtNPPzVOdG1bkUSWnDK3/VQ8ct+Xj78Mr/f8XscXgcJ6gTuWXQPl0+6fMxYOIN/x7L8/9u77/AoqvWB49/d9LqkkEoCAekBRAhFUEQUvBSvYkEQBPXaQBERFctP0OsVrNdrAcUCVrAQBRQElCaGXgMBFEiBkJBKes/5/TFkyaZAQnazm+T9PM8+kJkzM2dnZ3ffPefMey48Kv2tLhSsef2FArWsr/QGq2F91qqfSXnzTeMPc78nZ+I5apTJ8So2KcvMoOjECYpPnCRvz24K9+6r9/O1DwzEISgIh+AgyrKzydu8xWotYjYRaN1333112uFnn33WoAo1ZZZ4oaq2SLWePh2XK6+kJDmJ0uSzlJxN1v5NTqY0OZmyzMxL79TODsfQUBxCQ3AMCcUxNASHNiHn/22D3tm5xmPXduFb+8PhwvEr7aPK8ep8jFr2Ua/1SlUUrPaclVIX9nGR9adzEvjowCLO5CaiA3SVTsGsPk/SqVWnavswfR2qH6PaeuNJqaWOVdcn7oUdi1AVg5UrVuv0cMPL4NKq0vY176O2OkQnniNyz2l0lV5rHYpJ/UNp6+N24ZzWVsca16sLdTQ+v7qtN9a16jGq7eNS6y+jDpXOmbrIeqUUp9LziEnKprSsHDudjs7+brTzcUNXdfta3qOKi9ShtAjio0wuJQD8u4Oz14XzxEWegxnWVz0PJWXFpOWnUlSqDV1xtXfB28kbO52d5eqg6vA5Vpf1tbwOF15zYdTIP+xtItDS6/W0bduW3r17c7FNfvzxR7NVrqkx9wtVY1OuGeicndHZ2dXtw0Ep7fgVvxjlw0EIIZofne7Co9LfuprWVyqjO/+vKilBFRdf2J2zMzonJ+P6ivKqtJTynJw6Vakxh6rYRHqHhx9+mGXLlnHy5Enuu+8+Jk6ciLe3t1krI0wVx8WbPcgCUIWFNJsQqYEfDrXt41Lr0WndeOZcn19WQHpBGiXlpSgduDm442zvTFphuvaLW6fD380fg1Mri9Wh2nqdDvLS4ewhoBw8ArQWjbJiKCnQ0jg4OGv7qeF1uLB7HdRQpkwptp/MIDO/BHs7HSXlWsOCDugV6kVbX7dK3UI117Hib5Nj1Hl9LXWsbf1lHaPK+ooylzpGfdYDexMyWXPoLIVl5djp9dzQzY9BV7Qmp6iUtLxifN2dMLg61v0YFc+hIPN8EloD7FmCLuNvLdfWdc+cn5DadHt0ukrH4MKxLrm+4m/qvP50biKLDi4iLieeSu2R6PQ67u/xL4aGDq37MUzOyWWur3St6iq97qbr67qPhhzj/HnmEusrHcsc6jsIv/J2NY0Pc2wbatb6WUu9xmgVFRURGRnJZ599RlRUFKNGjeL+++9n+PDhZn/BmqJGadHS6wn56EPsfX3r/MYtTUun5MwZrS/cr7V8ONiY1PxUXtv1Gmvj1gIQ4BbAc/2eM35JNOoUOlX9tRa+nagFVt3HwtiPzZZXqbSsnIe/2sNvR1LwdLbnh0euxsPZ3jj/X6DBNucotGVnzhUwOzKaLX9pU7CEertwOrOAcgV6Hcwb24NxEQ348irMhq/Gwuld2oTUk1dBQA8z1f7ylJSX8N/d/+XLI6YTqbfkrPJNkbXTCdlE12FV8fHxLFmyhC+++IKSkhJiYmJa/J2HjTFGq6ndMitqV1Zexvd/fc//9v6P3JJc9Do9E7tOZNqV03B1MG/Cz8vy93pYNkELsrrdok2pY6YgSynF7OXRfLv7FE72er76V38i2kkLuTkopfhu9yleXqVlla/MTqdj6+yhDQtiC7Pgy7HafJYu3ueDrfAG1rphassq/8RVT3Bfj7qNLxbWd7ktYuZgE12HVenOtywopSi3QPeW0LS6/XbcBg+22sUnLONoxlFe3vYy0WnanZ3hPuG8OPBFuvp0tXLNzjv+Gyy7Wwuyut4Mt31i1gzhb6//i293n0Kvg/fG95Ygy4x0Oh3jIkJxdbTnsaX7TNaVKUVcWn7DAi1nA0yKhC9ugTN74YubtTsR/bs1rOINUFtW+f/u/S+7z+5mVsQs2htk5hJb5xAQ0Cy/4+p1z3RRURFLly7lxhtvpHPnzkRHR/P++++TkJDQ4luzLMkhIAC3/v2a5QXYkiTnJbPl9BZeinqJu36+i+i0aNwc3Hi237N8NfIr2wmyTmw4H2QVQZfRcPtnYGe+CWi/3BbHexu0yW9fuaUHw7vLdW0Jfdt5UTW1ll4H7XzN0FrqbIBJP0JQb8hP19JApBxp+H4vU4BbAHMGzjEmMdXr9AwMHIi9zp4/Ev/gthW38drO18gqyrJaHUXLVeeuw6lTp7Js2TJCQ0O59957mThxIj4+zW+OsYawxYSlwjYs/2s5L217ySQn1vC2w3mm3zP4ufpZsWZVnNwE34yD0kLoPAruWAL2jmbb/ZroJKZ+sxelYMYNHZlxQyez7VtUVzUJrI+bIxuevA6Dq5kC54JM+OKfkHQA3FrDlF+gtRXm3Dyv6njGuKw43trzFptObQLA4GRg2pXTuKPTHTJvojBhE2O09Ho9oaGh9O7d+6KDkSMjI81WuaZGAi1RlVKKn47/xItRL5os16Nn7e02NlA3dgt8fSeUFkCnm+DOL80aZG0/mc49n+6kuKycCf1D+c8t4c3uxgZblJRVwMFT55izMobk7EKu69yaTydHmC+TfH6G1n2YHA1ufueDLdsKoKPORPHGrjc4fk5rSe1g6MBTEU8xKHiQlWsmbIVNBFpTpkyp04fi4sWLG1yppkoCLVFBKcUfiX+wYP8CDqcfrrHMZyM+IyIgopFrVou4rfD1HVCSDx2Hw7ivwN7JbLs/mpzNHR9uI6ewlOHd/Fk4sU+Ln06nsR1KzOL2D6MoLCln2tAOPDWii/l2np8Bn98MZ6PBPUALtnyvMN/+zaC0vJTIvyN5f9/7ZBZlAnBN8DXMipiFq70rCdkJhHqG2taPH9FobCLQEpcmgZaoCLAW7l/IofRDADjbOVNUVmTSbWhTt57HR8FXt2lB1hU3wLivwcHZbLs/nZnPbQujOJtdREQ7L768vz/ODnZm27+ouxX7E3l82X4AFt59Ff/oYcapnPIqxmodBo9ALdjy6WC+/ZtJdnE2iw4s4uujX1NaXooevXEuUb1Oz5yBcxjbcayVaykamyW/v2UCMSHMQCnF1sStTFw9kWm/T+NQ+iFc7F24t/u9rL19LXOvnmsyUHfOwDk2EmRtg69u14KsDtebPcjKzCvmns92cja7iE7+7nxyT4QEWVb0zyuD+dfgMACe/P4Ax5LrlqG7Ttx8YPJK8OsGOUmwZDSknzDf/s3E09GTWRGz+OmfPzEgcIDJhO3lqpyXtr1Ecl6yFWsomhtp0TIjadFqeZRSRJ2JYsGBBRxMPQhoLVh3dbmLKd2n4ONy4YYRqyYerSorEY7+DL/N1YKs9kNh/FJwMF+S0PziUu7+ZAf7Es4RZHBm+dSrJQmpDSgtK2fy4p38eTydtj6urJw22HyD4wFyU+Hz0ZB6FDyDtfxr5SXg3QEMweY7jhnUln8rwj+CmX1nEu5r3fxgovFI12ETIYFWy5Ccl0x8djxpBWksPbqUA6kHAC3AGtd5HFPCp+Dr4mvlWl7E3i9g1eNQkXPItyM8uEWbTsdMSsvKefDLPWw4moLBxYEfHh5IR38Ps+1fNExGXjFj3ttK4rkCru3UmsVTzDg4HiA3RWvRSjt2YZlOD2P+B1fdY77jNFByXjIjlo+oln+rQh//PkzuNpkhIUOMLdKieZJAq4mQQKv5+/7Y9/x7+79Nxls52TlxZ+c7uS/8PtsOsEBryXon/EKQBaCzgxnRZmttUErxzPKDfLf7NE72er55oD992kpCUltz+EwWty3UBsc/cl0HnrnJjIPjAc4cgEXXmi4z87VmDpF/R/LStpcoV+XodXoe6vkQp3NOsyZ2DaWqFIB2nu2Y1G0SN3e4GWd783WtC9shgVYTIYFW83U65zRfxnzJN0e/MVmuQ8e3o7+1nWSjl7L7M/j5ierLJ/8MYdeY5RBvrj3G+xuPo9fBR5P6cmM3f7PsV5hf5cHxH0y4ilE9zTg4PnaLNji+KjNea+ZSU7f+2byzfHP0G74/9j05JdpYNi8nL+7qchfjOo/Dx8WH5LxkuVuxmbDJKXiEaO6Kyor4Pf53Io9HsiNpR41lFIrcktxGrtllSjoI6+dUX66zA2/zTE/yeVQc72/UchW9emsPCbJs3D+vDOZQYhYf/xHLrO8P0MHPjS4BZvqS8e6gdRdW7ZZzdDPP/s0owC2gWqDk7+bPE32e4MGeD/Lj3z/y1ZGvSMxNZOGBhXwa/SnhvuHsS9mHQsndiuKipEXLjKRFq3k4lnGMH4//yKoTq8guzga0lqur/K9i79m9tpum4WLOHtbGzBRkgFcYnEsAVaYFWWPeMcu4mV8OJvHoUi3r+8wbOzF9WMeG11tYXGlZOVMW72Lr8TRCvV1Z+eggWrmaKVHt3i9g1QztWqvgc4WW+sHDxt8zVZSWl/J7wu98fvhz4xyllTWZzwJRI+k6bCIk0GpaKjf7uzu4syZuDZF/RRrzX4H2S/fWK27llituIcg9qNp4jibxKzblKCwZBflpENxHm6OuKBcyTmotWWYYL7PtRDqTP9Oyvt/dP5RXJOt7k5KZV8zNH2zlVEYB13T0Zcm9/cw3OD4rUbvW7J3gh/shKwF8O8OUn8HdhqafqiOlFN8c+Yb5u+ZXW9fLtxd3db2L60Oux9XBfDeXCMuTQKuJkECr6agcMAE46B0oKS8BwF5vz9CQodzW8TYGBA7ATm+a98mm0jRcSupfWpCVlwKBveCeleDSymy7T8oqYNPRVF75JYa84jJu6h7AB3dfJVnfm6CYM9mMXfgnhSXlPDSkPc/+wwLjDjNitesxO1HLtzX5Zy3/VhNzqbsVne2cubbNtYwMG8ngNoNxsjPfLAvCMiTQaiIk0GoaavuQDPUI5c7OdzKmwxi8nZvBXXLpJ2DxSMhNhoAeWpDlar7nVXXC4jBfN9Y8fo0kJG3CVh44w/Sl+wB4b3xvxvQKMv9B0k9owVZOEvj30JKcmvG6bCxVW7cf6fUIpeWlrIldQ0JOgrGcu4M714dezz/C/kH/wP446B1kEL0NkkCriZBAq2moLUnhp8M/pV9gPyvUyAIyTsLiUZBzBvy6w+RVZm05SMoqYND8DcYgC0Cvgz9nXy9JSZu4eauP8NGWk7g42LH8kavpFmSBz7K047BkJOSehYCeWrDl4mX+41hYTa3bSiliMmJYc3INv8b9ytn8s8byXk5edGjVgT1n98ggehsjU/AIYUahnqHVkg/qdXpCPUOtVCMzy4yDJWO0IKt1F7hnhdm7Z2LT8kyCLIByBXFp+WY9jmh8T9/UhWs6+lJQUsZDX+0mM6/Y/AfxveJ88N8akg/Cl2OhMMv8x7GwALcAIgIiTFqldDod3X26MytiFutuX8eSm5YwrvM4vJ29ySzKZPfZ3cYbaspVOXOj5rI7eTfS5tF8SaAlWpwAtwDmDJxjm3MPNtS5U1ruouzT4NNR6y50b232w6TlFlVbZqfT0c5XBgA3dXZ6He+N702ItwunMgqYvmwfpWU1j0VqkNadz3dn+8CZvdrE5oXZ5j+OFel1evr49+GFAS/w+x2/M/OqmdXKKBT3rr2XG3+4kRe2vsCqE6tIzU+1Qm2FpUjXoRlJ12HT0qQGtddFVqLWHZMZp+UwmvILeJoxAeV5hSVl3PTOFuLS89EBCi3IenVsOOMimkmroOBIUjZjF0RRUFLGg9e257mRFkrKmxyt/TgoyISQATBxOTi5W+ZYVlbb+FB7nb0xC32FDoYO9A/sz4DAAfQN6IuHo4eM7bIgGaPVREigJawmO0kbYJxxArzawZTVFpvm5K11x3hvw3H8PZ348v7+pOcW087XVcZmNUM/HzzDo99og+NfurkbHf09CPN1M/9rfWY/fHGz1n3YdhDc/b1NJjY1h5pSxPwj7B/sS9nH9qTt7EjawZH0Iyb5+ux0dgS6BZKYm6iN7ULPnKtlbJc5SaDVREigJRpdViKc3g3rX4RzcdAqVAuyWoVY5HB/n81h5Lt/UFKmWHj3Vfyjh/lbzIRtmb/mKB9uPmH8W6+DeWN7mL/1MnEPfHELFGVDSH+45inw72ZT8yKay6Va088VnmPX2V1sP7OdHck7iM+Or3E/g4IG0ce/D918utHNpxtezk3vhgJbIYFWEyGBlmhUe7+AVY9fmOLExQse3AxebS1yuPJyxbhF29gVl8mwLn58MrmvJCVtAU5n5jP4tY0my+x0OrbOHmr+lq1Tu7SW2bLzYwB1ehjzP7PMXNCUrTm5hqf/ePqS5YLcguju290YeHX36Y7BySBdjnUgcx0KIUxlJZoGWaB1u+gt95b+bvcpdsVl4upox8uS+b3FSMiofidpmVLEpeWbP9DyDILziYMB7fpeNQM6DGuWLVt11du/N3qd3mRslx49D/R8gIScBGLSY4jPjudM3hnO5J1hffx6Y7lWTq04V3QO0KYSm3rlVO4Pvx8HO4fGfhotlgRaQjRFR1ZVn6xXlWv5syzwhZSaU8Srq48A2jyGwa1kPFZLEebrhl5HtZxpFrnDNONEDdd1GcT+AVfeZf7jNREVd0pfbPqvnOIcjqQfISY9hsPphzmcfphTOaeMQRZodzh+sP8DFu5fSIhnCGGeYYQZwmhnaEeYIYwwzzBaObcyOba0hjWcdB2akXQdikZx8Hv46RHTX/6gTRA9I9oigdaMZfv4af8Zugd5smLaIOztJDNMS/LtrgSei4ym7Py3RWd/d36dca35WzWzEuGd8OrBlrMB7loK7QaZ93hNTH3vlN6YsJHpG6fX6xheTl7G4Cu3OJf18etbxAB8GaPVRNh6oCW/TJo4pWDz67DpVe3vgJ5w9rD2i19nB2PeschYlj/+TmXSpzvR6+CnaYPo2aaV2Y8hbF9SVgFRJ9J55ocDlJbDokl9GN7dAp8je7/QugtVmTZGq1WolrJE7wD/fB96tdyWrfqqKZ2EXqfnm5HfkFuSS2xWLHHZccRmxRKbFUtSXtIl93lFqyto69mWIPcggtyCCHQPJMgtiCD3IDwdPY3Bd1P7vpFAq4mw5UCrpluKm+svk2aptAhWToeDy7S/r34MbnhZmy8u4yR4t7dIS1ZhSRkj3tlCfHo+U65ux9ybu5v9GKJpee3XoyzcdII2Xi78NnOIZea2zEq8cF27esOPD0PMT9q6IbPhutkgYwTrpD6f/fkl+cRnxxOXHcfWxK2sPLGyXsdytXclyD0InU7H35l/A9q4sEndJnHrFbfi6+KLwclQp5bQxg7UJNBqIhoj0LrYxVdYWkh6YTrpBecfhemkFaRxKvsUK0+avmH0Oj1rb1vbJH5ptHj5GfDtRIj/U2u5GvUW9L23UQ79xtqjfLDxBAGezqyfeS0ezjKAtqXLLy5l2FubScoqZPqwjsy8sZPlD1peDhv+DVvf1v7ucafWumXvZPljNwOXk5y5ttawuQPnUlBaQFJeEom5iSTlJnEm7wwZhRl12q+D3gFfF1/jo7VLa+3/rr74OmvLdp/dzTt73qGcxmsYkECribB0oBX5dyRzo+aaJLJriNs63kYbjzbVlle+JCofq9blletTpWq1lavLMWrdTz3rUZdtzVqPixyjPs9BKaXlFDr+GxTlaF0n7YeAZ3Cj1COnoIQNx1JRStEvzJtAg3OdzrFZ69GQ1/oyXvf6XjfmqkddrrcG16Ou19wl6qGUIruwhNOZBeh00N7XDUd7fd331ZB6FGZBXqpW2t4FPPxBZ1fv6+1y6mHOa64+dbBmPUrLSylTZca/7XR26HX6GutQNbO9OTVGw0CzDbS2bNnCG2+8wZ49e0hKSuLHH3/klltuMa5XSvHSSy+xaNEiMjMz6d+/Px988AHdu1/ovigqKmLWrFksXbqUgoIChg0bxoIFC2jT5kIAkZmZyfTp01m5UmvVufnmm3nvvfdo1aqVsUxCQgLTpk1jw4YNuLi4MGHCBN58800cHR3r/Hws+ULVNnWDEEII0dx9NuIzIgIiLLb/ZptHKy8vj169enHvvfdy2223VVv/+uuv8/bbb7NkyRI6derEK6+8wo033sixY8fw8PAAYMaMGaxatYply5bh4+PDk08+yejRo9mzZw92dtrYgQkTJnD69Gl+/fVXAB588EEmTZrEqlWrACgrK2PUqFG0bt2arVu3kp6ezuTJk1FK8d577zXS2bi4hOyEGoMse509rV1b4+Psg6+LL+6O7ujQ+r8r94PHZ8dzMPUgCoUOHb1a96KtZ1uTMhXbVd228vLKat22lv2YbFuHY1Xd1lz1MNmPyX8tcC7qUDeTXRbmoCvIhNwUOPE7OlUGHkHQ8050Th6XPIZJnepZj8pl9iWc45foJBztdDxy3RV4nu8yrOtrVZcy9b5mGvI612U/Frjmq+2nAdebua5bc10vZ7IKjHciPnljJ3qHetVr+wbV41wCrPs/yE3W3hc3vATuAZCTCJ7B4Nba4p8LF92+DmUs/flUp/004HW6WLk1sWt4e8/bxnFhM/vMZFT7UXXaNjU/lTtX3Uk5pl2WIR6Wme2iMdhM16FOpzNp0VJKERQUxIwZM3jmmWcArfXK39+f1157jYceeoisrCxat27Nl19+ybhx4wA4c+YMISEhrF69mhEjRnDkyBG6devG9u3b6d+/PwDbt29n4MCBHD16lM6dO7NmzRpGjx7NqVOnCAoKAmDZsmVMmTKFlJSUOke3jd2iVd/m1GY3iXJzUzXTO0DXMXDrInC0QM6iWqTkFHLDW5vJLizlxdHduG9wWKMdWzQt81Yf4aMtJwn1dmXdE9daZmB8bXJTYOld2tQ9Orvz7xuFZJO3DQ35vrHGzVuW/P622WQ4sbGxJCcnM3z4cOMyJycnhgwZQlRUFAB79uyhpKTEpExQUBDh4eHGMtu2bcNgMBiDLIABAwZgMBhMyoSHhxuDLIARI0ZQVFTEnj17aq1jUVER2dnZJg9LqUhYp9dpL1nFxVefCzjALYCIgAgJsmxRTZne0cGIeY0aZAG88vMRsgtL6RFsYPLV7Rr12KJpeWxYR/w9nUjIyOejzScb9+DufjD5Z7jiRi0VRMVgoYps8lmJjVsfYaIh3zdjO45l7W1r+WzEZ6y9bW2Tv0PeZgOt5ORkAPz9/U2W+/v7G9clJyfj6OiIl5fXRcv4+flV27+fn59JmarH8fLywtHR0VimJvPmzcNgMBgfISGWbdpsbhefOE8p2LO4epJGlJY/qBFtOpbCygNn0Ovg1Vt7YKeXW+hF7dyd7Hl+VDcAFmw6zqkapuuxKEdXGPhY9eWqTEsPIZqs5tQwYLOBVoWq/cNKqUvm4Khapqbyl1OmqmeffZasrCzj49SpUxetlzk0p4tPoP3qXjoetrxRfZ3OTssj1EhiU/N46ocDAEy5OowebQyNdmzRdI3pGcjA9j4UlZbz8s8xjV8B3yu07sKq9n6uzf8phJXZbKAVEKAFElVblFJSUoytTwEBARQXF5OZmXnRMmfPnq22/9TUVJMyVY+TmZlJSUlJtZauypycnPD09DR5CFEn5eWw61P4oD/8tUZL3dDpH1pwBRgzvTfSRLrf7krg+rc2kZpTDEA7n8btrhRNl06n46V/dsder2N9zFk2Hktp3AoYgrUxWRXvnYpB1tHfa++vo6sbtz5CVGGzgVZYWBgBAQGsX39hFvLi4mI2b97M1VdfDUCfPn1wcHAwKZOUlMShQ4eMZQYOHEhWVhY7d+40ltmxYwdZWVkmZQ4dOkRS0oXpB9atW4eTkxN9+vSx6PMULVDacfh8NPwyE4pzILgvPPwHTFimzVU4+Wft30YazJuUVcDsyGiTDDsvrYohKaugUY4vmr5O/h7cO6gdAC+tPExRadnFNzC3q+658N554jBM+UVrDc5JgmXj4fsp2uB5IazAqukdcnNzOX78uPHv2NhY9u/fj7e3N6GhocyYMYNXX32Vjh070rFjR1599VVcXV2ZMGECAAaDgfvvv58nn3wSHx8fvL29mTVrFj169OCGG24AoGvXrtx000088MADfPTRR4CW3mH06NF07twZgOHDh9OtWzcmTZrEG2+8QUZGBrNmzeKBBx6QViphPmUlEPUubHoNyorAwRWGvQj9HgT9+V/jhuBGa8WqEJuaR9V7j8uUIi4tn0CDS6PWRTRdj9/QiRX7zxCXns/HW07y6PUdG7cCld87hmB4JAo2zYeo9+Dwj3Byk3ZzSa+7kOl7RGOyaqC1e/duhg4davx75syZAEyePJklS5bw9NNPU1BQwNSpU40JS9etW2fMoQXw3//+F3t7e+68805jwtIlS5YYc2gBfP3110yfPt14d+LNN9/M+++/b1xvZ2fHL7/8wtSpUxk0aJBJwlIhGiQrETJOQEmBNoVIcrS2vMP1MPod8Gpr1eoBxKbnVVtmp9PRzle6D0XdaQPju/L4sv28v/E4t/QOpo2XFa8hBxe48SXofiusfFR77/30MER/p7339Pbae9O7Q6P/uBEti83k0WoObHlSaWEFNeXFcvGCm+ZDz3E28au6sKSMYW9tJvFcATq0G+TtdDpeHRvOuIhQa1dPNDFKKe5atJ0dsRlc16k1Dw5pT5ivm/VbRstKtJatTfO11mQ7R22Z5N0S5zXbKXiaGwm0hFHKMVjQH9PJH3Xw0BYI7GmtWlXz7u9/8/b6vwg0OPPNvwaQnF1IO19X638xiibrWHION/1vi7E7Wq+DeWN72EbgnvY3RD4EZ6rkR9TZaWO8pGWrxWqRCUuFaJJykmH9HPj4OqrNsI2yqdvNz5wrYMEmbYzksyO7EtbajYEdfCTIEg3i6WJvcumXK3gu8pBt3Fzh2xGGzam+XJXB5tdt6v0pmg+rjtESotlI+xv+/B8c/BbKimsu08h5sS7ltV+PUlhSTkQ7L8b0DLR2dUQzEZuWV+0nhk3dXFGRd6tqguC9S+DQcuh7LwyYCp7ynhDmIS1aQjTEqZ2w7G54PwL2fakFWSEDYPwyGPOu1fJiXcruuAxW7D+DTgdzxnS/ZBJgIeoqzNeNqhMK6HXYzs0VVfNu6ezgyonQuouWbiXqXXinB6yYBqnHrFtX0SxIi5YQ9ZGVqOXByj4F+76ChG0X1nUeBYOmQ+iAC8uuuEGbCsS7vc0EWeXlipdWaRm8x/UNITxYMsAL8wk0uDBvbA+ei4ym7HzTVo9gg220ZlW46h7oMMz0vVleDn+vgz/f0d7X+77SHp1HwaDHwdBG7lIUl0UGw5uRDIZv5rYtgLXPYTIAxc5Ru4Pw6unQupPVqlYf3+06xdPLD+LhZM+GWdfR2sPJ2lUSzVBSVgEbj6bw3I+HAFj24AAGtPexcq3qKGGH1rJ19Ofq6+QuxWbJkt/f0qIlxMXkpsCRVXDwOzi1vcpKHdy3DoJ7W6VqlyOnsITX1x4FYPqwjhJkCYsJNLgwoX9bYpKy+Wp7AnNWHOaX6YOxt2sCI1ZC+0Po15D6F2yaB4cjL6xT5bByOhTnw5XjwVlahMXFNYErXohGln0GdnwEi0fCm520qXKqBVkACopzG716DfH+huOk5RYT5uvG5KvbWbs6ogWYNbwzXq4OHDubw5fb461dnfpp3UkbHF+Ngl+fgTeugK/vhH1fQ35Go1dPNA3SoiVaroqs7d4dtNu7Y1ZCzAo4vdO0XHAfCBuijd2ofKeSjd1FeCmxaXl89mcsAP83uiuO9vI7S1heK1dHZo3ozPM/HuLt9X8xplcQvu5NqCXVu0MNdynqwCsMMk/C32u1h95e+5zodjN0GQ2lRTKmSwASaImWqqas7ZWFDIBu/4SuY6BViLbMOwxWzdCCMhu7i7Au/vNLDCVliiGdWjO0s5+1qyNakLsiQlm6M4FDidm8/utRXr+9l7WrVHcVdylWfe9fdQ+kHIUj53+gnT0EJ37XHqtmYBzLqdNrU/70mWy1pyCsSwbDm5EMhm8ishLhnfDqQVabftDzTu3XaG05dLISbe4uwrrY/Fcqkz/bib1ex68zruUKP3drV0m0MHviM7ltYRQAP069mt6hXlauUT1d6r2fdhyOrIDo7yHlSPX17a+HDkMhdCAE9gJ7R8vXWdSZDIYXwpwyTtTckjXsRQi75uLbGoKbVIAFkJCRx7ORBwGYfHU7CbKEVfRp68Xtfdrww57TvLjiMD9NG4Rd1YRbtuxS733fK+CaJ6FNBHw+pvr6kxu0B4C9szYkIXSAFni1iQCXVtq6ykMamthnjaiZBFqi5alpzEUTG29VV9/uSmD28mhjQoo2XjaUy0i0OM/c1IW1h5KJTsziu92nGN/PBuY/NLcaP1/0MHim1tJ1ajvkp0P8n9pDKwB+3bRJ5+P/RCa7bl6k69CMpOuwCdn7Rc1jLpqRpKwCBs3fQHnltF86HVtnD7Wt5JGiRflsaywv/xyDl6sDG2ddRyvXZtiFdrHPF6Ug/biWFDVhu/Zvxsna99XpJgjuC35dtGDMqx3o7UzLSCtYg1ny+1sCLTOSQKuJaaLjrepq699pTPx0R7XlSx8YwMAOTSRxpGh2SsvKGfXuVo6dzWFs72Bu79uGMF+35hf81+fzJTcFdn+m5ey6FHtn8O2kBV1+XbSJ7Hcu0lrQpBXsskmg1URIoCVsyfw1R/hws+kvZWnRErZg24l0xn98ITedXgfzxvZgXEQz7Eqsq5pu0tHpYdAMyEmClBgtgWppwaX31fMu8O8GrULPP9qCqw9UndNUWsKMZDC8EKJeDp4+xyd/aDmzdDqtt8JOp+PVseESZAmrqzrBdLmC5yIPcW2n1i33+rxYGokK5WWQGQepR7XA6+QWiNtSfV8Hl1Vf5uBWKfAKhbxULS1FxXiwkW9BxH11r68EaXUmLVpmJC1awhbkF5cy+t2tnEzLY2SPAF4Y1ZX49ALa+bq23C8xYVOiTqQx4WPp1q5Rfboca2wF00HEg1CQCecS4Fy81iJWF86twDMI3P3APeD8v/7gEWC6LGYF/Dzj8rorbTRAkxYtIUSd/eeXI5xMy8Pf04lXb+1BK1dHglq5XnpDIRpJmK8beh0mN2rodNVbulqk+qSQqUsrGGhZ6rNOa0FXZjzEbYVDP1TfX+E57ZESU/f6Vsz9mLANDCHg4g2u3uf/9brwt5Mn7PvyQqLoyxlPZqNB2qVIoCVEM/L7kbN8vSMBgLfuuLJ53tElmrxAgwvzxvbguchDlJ3vVHHQ68grKrNyzZqgq+6BDsMu3gpm7wQ+HbQHQMfh2kTZVVPcTIwEyrXB+TnJ2r+5Z00fhVk1VELB/m8uUVG9tm/jJucDtNgt4BEIzp7gZAAnj/P/97zwr5OnloH/l5lNctC/dB2akXQdCmtKzSnipne2kJ5XzP2Dw/i/0d2sXSUhLiopq4ATKbm8vf4v9iaco5O/Oz9NG4Sro7QBWNzlprhJPwHv960+92O/B6G8RJtcuyATCjIg//y/Jfnmr7/ODmZEm61lS7oOhRAXpZTimeUHSc8rpkuAB0+N6GztKglxSYEGFwINLnQK8GD0u1v562yuNvn0nb3QVb1DTphXXVrCauLToW7dlZWVFGrdkZ8Mqx6gDZgGOrSWsqIcKMqGwmzt36Ic7f8ledX3qcq0ujeBLkQJtIRo4pKyCvh4y0k2HE3B0V7PO3ddibOD3aU3FMJG+Hk489743kz4ZAc/7kskop03E/q34FQPjeVypxSrb5Dm4AzBV9U/QKuQGQ/vXtlkZ/OQQEuIJuzbXQk8GxltHFQ8vJs/XQKk21o0Pf3b+/D0iM7MW3OUuSsP0yPYQI82BmtXS9TmcoK0y21F82pbc5DWBFqzQMZomZWM0RKNqeYpdmDr7OsljYNokpRSPPjlHtbHnCXE24VP74kgLa+oeWaOF/Vnwdk8ZIyWEKKag6ezTIIsgDIFcWn58qUkmiSdTsebd/Ri9Ht/cCqjgOHvaMk4JXO8AC6/q9PK9NaugBCi/lJyCpm/5ki15XY6neQiEk2awcWBf/8z3GRZReb4pKw6TD8jhI2RQEuIJubMuQLGfbSd2LR8PJ3t0Z+/OUum2BHNhaN99a+mMqWIS7NAmgAhLEy6DoVoQhLS85nwyXZOZxYQ3MqFbx7oj6O9nri0fJliRzQbNWWOB/DzcLJOhYRoAGnREqKJOJGay50fbeN0ZgHtfFz57uGBtPXRBgkP7OAjQZZoNioyx9tVyaU1fdk+zmYXWqlWQlweuevQjOSuQ2EJSVkFbPkrlflrjpKZX0JHP3e+/ld//DydrV01ISwqKauAuLR8copKeHZ5NOl5xQQanPl0cgTdguQzVpiP3HUoRAv17a4EZkdGU/FzKNDgzLIHB+DjLl0oovmryBwP0DXAkymLd3IiNY87Pozi/QlXMbSLn5VrKMSlSdehEDbqVEYes5dfCLIAzmYXUlxWXvtGQjRTId6uRD4yiKs7+JBXXMb9n+/i/Q1/E3UiTe5GFDZNAi0hbNDuuAzu/mQnVfv1y8/nyRKiJTK4OrDk3n7c0acN5QreXPcXEz7ewaD5G/h2V4K1qydEjSTQEsKGpOYU8eR3B7j9w20kZFQPqCRPlmjpHO31PHFjRyoPky9XMDsymoT0GiYfFsLKJNASwoqSsgqIOpHGqYw8lvwZy/VvbWL53tMAjOsbwv+N6mq880ryZAmhiUvPr9baqxTc/uE2VuxPpLxqXgghrEgGwwthJVUnhK7QI9jAy//sTu9QLwBG9gyUPFlCVFJbnq2UnCIeX7afDzef5OmbOnNdp9YkZxcSm5Yn8yUKq5H0DmYk6R1EXZ05l8+g1zZS9d339E2deejaDtjpdTVvKIQAtB8qz0Ueokwp7HQ65ozpRnZhCR9tPklOUSmgBWRx6XkoJfMliouT9A5CNBOnMvJZvvc03+xIqBZkAfQO8ZIgS4g6GBcRyrWdWldr7b27f1sWbDrOkqg4YtMujNkqV/BsZDTXdmotLVuiUUmgJYQFJGUVGLsr3J3sWROdzA97T7MzNqPWbWSguxD1UznPVgUvN0eeH9WNnm1a8djSfSbryhU8+d0B7hnYlus6++HsYNeY1RUtlARaQphZ1bFX9nY6Ssu0P3Q6GHyFL7dd1YacwlLmrjxs7PqQge5CmE/fdl41juOKOpFO1Il0XB3tuL6LH6N6BHJdZz/OFRTLWC5hETJGy4xkjFbLVVJWzuEz2Ww4cpZ3Nxyvtj7U24VxEaHc2juYoFYXPsQrphiRge5CmJ/pOC54aEgHikvLWXMomcRzF5KcOtjpKDn/Y0jGcrVMlvz+lkDLjCTQat4qdwd6ODuwLyGTXXGZ7IrNYN+pTApLas/YvvSB/gzs4NuItRVCQM0/ZpRSHDidxeroJFYeOENyVvWJqq/v4seA9t70atOK8GADbk721fYrLWDNhwyGF8IC6vpBmVtUyoebT/DBhuPG3D06qJbHp5WrA+FBBv48nmayTht75Wbm2gsh6qKmcVw6nY4rQ1pxZUgrruvUmgmf7Ki23YajKWw4mgJorVyd/D3o1aYVvUJacTa7kPc2/E253M0o6kACLdEiVR5HpdfBf27twTUdfTmRmsfJ1FxOpOZyMjWPE6m5nM0uqra9AgJbOTMwzIe+7byJaOdFh9bu6PW6aredy9grIWxXWOvqObn057sZY1PzOHD6HElZhRxNzuFocg7f7j5lsn25gtnLo8nMK6FniIF2Pm4EeDqjr3L3sLSAtVzSdWhG0nXYNCRlFTBo/oZqg2Tra+kDAxjYwafWY8jYKyGahpp+HFVuoTqbXcj+U+c4cOocm/9K5fCZ7Ivuz9FeT1tvV9r6uNHOx5W03CJWHDhjzOf16q09uKvfxVvAJDBrXDJGq4mQQKtpiDqRxoSPq3cV2Omhva87HVq70761Gx1au9PBzx1XRztuemeLSWBmp9OxdfZQ+QAUopmo64+jmn6o6YD+7b05m13EqYx8SuvwKy7A05kAgzOtPZzw83DCz+PC//efPseCjccvq2tSArTLI4FWEyGBVtNQ0welXgebn7qOEO+ax1Jd6hevEKLluNjnQWlZOWfOFRKXnkd8Rj7bjqex+lByg495hZ87rd2dMLg4YHBxoJWrA57n/zW4ONDKxZEdsel8cJkBWksngVYTIYFW03E5gZN0BwohKjSkBUyvg48m9aFcQWpOESk5RaTmFJKaU8TJ1DxOVspo3xDS8l53ctehEGZW2/QdF1PT3UtCiJaprp8HgQYX5o3tUe2H3Y3dAmosX1tg9tadV2Kn15GVX0xWQQnn8ku0fwu0f5POFXAqs8BkX2VKEZeWL59bViaBlmixJHASQjSG+vywqy0wu7V38EWPUVOAJtN62QYJtIQQQggLq88Pu8ttca8pQJMfk9YngZYQQghhYy6nxf1yAjRheRJoCSGEEM2EDImwPXprV0AIIYQQormSQEsIIYQQwkIk0BJCCCGEsBAJtIQQQgghLEQCLSGEEEIIC5FASwghhBDCQiTQEkIIIYSwEAm0hBBCCCEsRAItIYQQQggLkUBLCCGEEMJCJNASQgghhLAQmevQjJRSAGRnZ1u5JkIIIYSoq4rv7YrvcXOSQMuMcnJyAAgJCbFyTYQQQghRXzk5ORgMBrPuU6csEb61UOXl5Zw5cwYPDw90Ol29ts3OziYkJIRTp07h6elpoRqKyuScW4ec98Yn59w65Lw3vss950opcnJyCAoKQq8376gqadEyI71eT5s2bRq0D09PT3lDNjI559Yh573xyTm3Djnvje9yzrm5W7IqyGB4IYQQQggLkUBLCCGEEMJCJNCyEU5OTsyZMwcnJydrV6XFkHNuHXLeG5+cc+uQ8974bPGcy2B4IYQQQggLkRYtIYQQQggLkUBLCCGEEMJCJNASQgghhLAQCbSEEEIIISxEAi0bsGDBAsLCwnB2dqZPnz788ccf1q6STZo3bx4RERF4eHjg5+fHLbfcwrFjx0zKKKWYO3cuQUFBuLi4cN1113H48GGTMkVFRTz22GP4+vri5ubGzTffzOnTp03KZGZmMmnSJAwGAwaDgUmTJnHu3DmTMgkJCYwZMwY3Nzd8fX2ZPn06xcXFFnnutmLevHnodDpmzJhhXCbn3DISExOZOHEiPj4+uLq6cuWVV7Jnzx7jejnv5lVaWsoLL7xAWFgYLi4utG/fnpdffpny8nJjGTnnDbdlyxbGjBlDUFAQOp2On376yWS9rZ3j6OhohgwZgouLC8HBwbz88sv1nw9RCatatmyZcnBwUB9//LGKiYlRjz/+uHJzc1Px8fHWrprNGTFihFq8eLE6dOiQ2r9/vxo1apQKDQ1Vubm5xjLz589XHh4eavny5So6OlqNGzdOBQYGquzsbGOZhx9+WAUHB6v169ervXv3qqFDh6pevXqp0tJSY5mbbrpJhYeHq6ioKBUVFaXCw8PV6NGjjetLS0tVeHi4Gjp0qNq7d69av369CgoKUo8++mjjnAwr2Llzp2rXrp3q2bOnevzxx43L5ZybX0ZGhmrbtq2aMmWK2rFjh4qNjVW//fabOn78uLGMnHfzeuWVV5SPj4/6+eefVWxsrPr++++Vu7u7euedd4xl5Jw33OrVq9Xzzz+vli9frgD1448/mqy3pXOclZWl/P391V133aWio6PV8uXLlYeHh3rzzTfr9Zwl0LKyfv36qYcffthkWZcuXdTs2bOtVKOmIyUlRQFq8+bNSimlysvLVUBAgJo/f76xTGFhoTIYDOrDDz9USil17tw55eDgoJYtW2Ysk5iYqPR6vfr111+VUkrFxMQoQG3fvt1YZtu2bQpQR48eVUppHxZ6vV4lJiYayyxdulQ5OTmprKwsyz1pK8nJyVEdO3ZU69evV0OGDDEGWnLOLeOZZ55RgwcPrnW9nHfzGzVqlLrvvvtMlo0dO1ZNnDhRKSXn3BKqBlq2do4XLFigDAaDKiwsNJaZN2+eCgoKUuXl5XV+ntJ1aEXFxcXs2bOH4cOHmywfPnw4UVFRVqpV05GVlQWAt7c3ALGxsSQnJ5ucTycnJ4YMGWI8n3v27KGkpMSkTFBQEOHh4cYy27Ztw2Aw0L9/f2OZAQMGYDAYTMqEh4cTFBRkLDNixAiKiopMuneai2nTpjFq1ChuuOEGk+Vyzi1j5cqV9O3blzvuuAM/Pz969+7Nxx9/bFwv5938Bg8ezO+//85ff/0FwIEDB9i6dSsjR44E5Jw3Bls7x9u2bWPIkCEmyU9HjBjBmTNniIuLq/PzkkmlrSgtLY2ysjL8/f1Nlvv7+5OcnGylWjUNSilmzpzJ4MGDCQ8PBzCes5rOZ3x8vLGMo6MjXl5e1cpUbJ+cnIyfn1+1Y/r5+ZmUqXocLy8vHB0dm91rt2zZMvbu3cuuXbuqrZNzbhknT55k4cKFzJw5k+eee46dO3cyffp0nJycuOeee+S8W8AzzzxDVlYWXbp0wc7OjrKyMv7zn/8wfvx4QK71xmBr5zg5OZl27dpVO07FurCwsDo9Lwm0bIBOpzP5WylVbZkw9eijj3Lw4EG2bt1abd3lnM+qZWoqfzllmrpTp07x+OOPs27dOpydnWstJ+fcvMrLy+nbty+vvvoqAL179+bw4cMsXLiQe+65x1hOzrv5fPvtt3z11Vd88803dO/enf379zNjxgyCgoKYPHmysZycc8uzpXNcU11q27Y20nVoRb6+vtjZ2VX7hZKSklIt0hYXPPbYY6xcuZKNGzfSpk0b4/KAgACAi57PgIAAiouLyczMvGiZs2fPVjtuamqqSZmqx8nMzKSkpKRZvXZ79uwhJSWFPn36YG9vj729PZs3b+bdd9/F3t7e5NddZXLOGyYwMJBu3bqZLOvatSsJCQmAXOuW8NRTTzF79mzuuusuevTowaRJk3jiiSeYN28eIOe8MdjaOa6pTEpKClC91e1iJNCyIkdHR/r06cP69etNlq9fv56rr77aSrWyXUopHn30USIjI9mwYUO1ZtuwsDACAgJMzmdxcTGbN282ns8+ffrg4OBgUiYpKYlDhw4ZywwcOJCsrCx27txpLLNjxw6ysrJMyhw6dIikpCRjmXXr1uHk5ESfPn3M/+StZNiwYURHR7N//37jo2/fvtx9993s37+f9u3byzm3gEGDBlVLXfLXX3/Rtm1bQK51S8jPz0evN/1KtLOzM6Z3kHNuebZ2jgcOHMiWLVtMUj6sW7eOoKCgal2KF1XnYfPCIirSO3z66acqJiZGzZgxQ7m5uam4uDhrV83mPPLII8pgMKhNmzappKQk4yM/P99YZv78+cpgMKjIyEgVHR2txo8fX+OtwW3atFG//fab2rt3r7r++utrvDW4Z8+eatu2bWrbtm2qR48eNd4aPGzYMLV3717122+/qTZt2jSL268vpfJdh0rJObeEnTt3Knt7e/Wf//xH/f333+rrr79Wrq6u6quvvjKWkfNuXpMnT1bBwcHG9A6RkZHK19dXPf3008Yycs4bLicnR+3bt0/t27dPAertt99W+/btM6Y0sqVzfO7cOeXv76/Gjx+voqOjVWRkpPL09JT0Dk3RBx98oNq2bascHR3VVVddZUxXIEwBNT4WL15sLFNeXq7mzJmjAgIClJOTk7r22mtVdHS0yX4KCgrUo48+qry9vZWLi4saPXq0SkhIMCmTnp6u7r77buXh4aE8PDzU3XffrTIzM03KxMfHq1GjRikXFxfl7e2tHn30UZPbgJurqoGWnHPLWLVqlQoPD1dOTk6qS5cuatGiRSbr5bybV3Z2tnr88cdVaGiocnZ2Vu3bt1fPP/+8KioqMpaRc95wGzdurPFzfPLkyUop2zvHBw8eVNdcc41ycnJSAQEBau7cufVK7aCUUjql6pviVAghhBBC1IWM0RJCCCGEsBAJtIQQQgghLEQCLSGEEEIIC5FASwghhBDCQiTQEkIIIYSwEAm0hBBCCCEsRAItIYQQQggLkUBLCCGEEMJCJNASQjQ5S5YsoVWrVvXaJi4uDp1Ox/79+wHYtGkTOp2Oc+fOmb1+tuDYsWMEBASQk5PToP1EREQQGRlpploJ0fJIoCWEuKSKoKS2x9ChQ61dxXq7+uqrSUpKwmAwXLJsUwzKnn/+eaZNm4aHh0eD9vN///d/zJ492zi5shCifiTQEkJcUkVQUvXx0UcfodPpmDp1qrWrWG+Ojo4EBASg0+msXRWzO336NCtXruTee+9t8L5GjRpFVlYWa9euNUPNhGh5JNASQlxSRVBS+ZGZmclTTz3Fc889xx133GEsu3nzZvr164eTkxOBgYHMnj2b0tJS4/rrrruO6dOn8/TTT+Pt7U1AQABz5841Od7bb79Njx49cHNzIyQkhKlTp5Kbm1uvOu/cuZPevXvj7OxM37592bdvn8n6qq1U8fHxjBkzBi8vL9zc3OjevTurV68mLi7O2GLn5eWFTqdjypQpAPz6668MHjyYVq1a4ePjw+jRozlx4oTxGBXdlZGRkQwdOhRXV1d69erFtm3bTOry559/MmTIEFxdXfHy8mLEiBFkZmYCoJTi9ddfp3379ri4uNCrVy9++OGHiz737777jl69etGmTRvjsoru1p9//pnOnTvj6urK7bffTl5eHp9//jnt2rXDy8uLxx57jLKyMuN2dnZ2jBw5kqVLl9br/AshNBJoCSHq7dy5c9xyyy0MGTKEf//738bliYmJjBw5koiICA4cOMDChQv59NNPeeWVV0y2//zzz3Fzc2PHjh28/vrrvPzyy6xfv964Xq/X8+6773Lo0CE+//xzNmzYwNNPP13n+uXl5TF69Gg6d+7Mnj17mDt3LrNmzbroNtOmTaOoqIgtW7YQHR3Na6+9hru7OyEhISxfvhzQxj0lJSXxv//9z3icmTNnsmvXLn7//Xf0ej233nprtW62559/nlmzZrF//346derE+PHjjcHn/v37GTZsGN27d2fbtm1s3bqVMWPGGIOdF154gcWLF7Nw4UIOHz7ME088wcSJE9m8eXOtz2XLli307du32vL8/Hzeffddli1bxq+//sqmTZsYO3Ysq1evZvXq1Xz55ZcsWrSoWiDXr18//vjjj0ucdSFEjZQQQtRDWVmZ+sc//qG6du2qsrKyTNY999xzqnPnzqq8vNy47IMPPlDu7u6qrKxMKaXUkCFD1ODBg022i4iIUM8880ytx/zuu++Uj4+P8e/Fixcrg8FQa/mPPvpIeXt7q7y8POOyhQsXKkDt27dPKaXUxo0bFaAyMzOVUkr16NFDzZ07t8b9VS1bm5SUFAWo6OhopZRSsbGxClCffPKJsczhw4cVoI4cOaKUUmr8+PFq0KBBNe4vNzdXOTs7q6ioKJPl999/vxo/fnyt9ejVq5d6+eWXTZYtXrxYAer48ePGZQ899JBydXVVOTk5xmUjRoxQDz30kMm2K1asUHq93vgaCiHqTlq0hBD18txzz7Ft2zZWrFiBp6enybojR44wcOBAk3FPgwYNIjc3l9OnTxuX9ezZ02S7wMBAUlJSjH9v3LiRG2+8keDgYDw8PLjnnntIT08nLy+vTnU8cuQIvXr1wtXV1bhs4MCBF91m+vTpvPLKKwwaNIg5c+Zw8ODBSx7nxIkTTJgwgfbt2+Pp6UlYWBgACQkJJuUqP9/AwEAA4/OtaNGqSUxMDIWFhdx44424u7sbH1988YVJF2VVBQUFODs7V1vu6upKhw4djH/7+/vTrl073N3dTZZVfi0AXFxcKC8vp6ioqNZjCiFqJoGWEKLOvv32W958802WLVtGx44dq61XSlUbXK6UAjBZ7uDgYFJGp9MZu9vi4+MZOXIk4eHhLF++nD179vDBBx8AUFJSUqd6VhyzPv71r39x8uRJJk2aRHR0NH379uW999676DZjxowhPT2djz/+mB07drBjxw4AiouLTcpVfr4V56Hi+bq4uNS6/4oyv/zyC/v37zc+YmJiLjpOy9fX1zjGq7Z6VNTlYq9FhYyMDFxdXS9aVyFEzSTQEkLUyf79+7nvvvuYP38+I0aMqLFMt27diIqKMgl0oqKi8PDwIDg4uE7H2b17N6Wlpbz11lsMGDCATp06cebMmXrVtVu3bhw4cICCggLjsu3bt19yu5CQEB5++GEiIyN58skn+fjjjwHtZgDAZJB4eno6R44c4YUXXmDYsGF07dq1xuDmUnr27Mnvv/9e6/NwcnIiISGBK664wuQREhJS6z579+5NTExMvetSm0OHDnHVVVeZbX9CtCQSaAkhLiktLY1bbrmF6667jokTJ5KcnGzySE1NBWDq1KmcOnWKxx57jKNHj7JixQrmzJnDzJkz0evr9nHToUMHSktLee+99zh58iRffvklH374Yb3qO2HCBPR6Pffffz8xMTGsXr2aN99886LbzJgxg7Vr1xIbG8vevXvZsGEDXbt2BaBt27bodDp+/vlnUlNTyc3NxcvLCx8fHxYtWsTx48fZsGEDM2fOrFc9AZ599ll27drF1KlTOXjwIEePHmXhwoWkpaXh4eHBrFmzeOKJJ/j88885ceIE+/bt44MPPuDzzz+vdZ8jRoxg27ZtJoFhQ/zxxx8MHz7cLPsSoqWRQEsIcUm//PIL8fHxrF69msDAwGqPiIgIAIKDg1m9ejU7d+6kV69ePPzww9x///288MILdT7WlVdeydtvv81rr71GeHg4X3/9NfPmzatXfd3d3Vm1ahUxMTH07t2b559/ntdee+2i25SVlTFt2jS6du3KTTfdROfOnVmwYIHxeb300kvMnj0bf39/Hn30UfR6PcuWLWPPnj2Eh4fzxBNP8MYbb9SrngCdOnVi3bp1HDhwgH79+jFw4EBWrFiBvb09AP/+97958cUXmTdvHl27dmXEiBGsWrXKOB6sJiNHjsTBwYHffvut3vWpKjExkaioKLPk5BKiJdKpyxnMIIQQwqYtWLCAFStWNDjR6FNPPUVWVhaLFi0yU82EaFnsrV0BIYQQ5vfggw+SmZlJTk5Og6bh8fPzu2QOMiFE7aRFSwghhBDCQmSMlhBCCCGEhUigJYQQQghhIRJoCSGEEEJYiARaQgghhBAWIoGWEEIIIYSFSKAlhBBCCGEhEmgJIYQQQliIBFpCCCGEEBYigZYQQgghhIX8P6qoMmMNrlXLAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ds = xr.open_zarr(\"PeriodicParticle.zarr\")\n", "\n", "plt.plot(ds.lon.T, ds.lat.T, \".-\")\n", "plt.xlabel(\"Zonal distance (m)\")\n", "plt.ylabel(\"Meridional distance (m)\")\n", "plt.show()" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "We can see that the particles start at 20,000m, move eastward, and once they hit the boundary at 90,000m, they jump to the other side of the domain (the horizontal lines). So we have periodic boundary conditions!\n" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "As a note, one may ask \"why do we need the halo? Why can't we use simply the `PeriodicBC` kernel?\" This is because, if the particle is close to the edge of the fieldset (but still in it), `AdvectionRK4` will need to interpolate velocities that may lay outside the fieldset domain. With the halo, we make sure `AdvectionRK4` can access these values.\n" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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.11.6" } }, "nbformat": 4, "nbformat_minor": 1 }